What ab initio methods are better suited to modelling disorder in materials?

Matter Modeling Asked on November 7, 2021

In contrast to the perfect periodic bulk materials, now computational approaches are moving towards modeling ‘real’ materials or solid solutions with dopants, dislocations, grain boundaries, and interfaces [1]. For example, the special quasirandom structures (SQS) approach is one way of modeling random alloys with statistical site occupations [2].

So my question is, what are other ab initio/theoretical approaches of modeling disorder in materials and what disorder are they good at modeling?

References

[1] Yang, Y., Chen, C., Scott, M. et al. Deciphering chemical order/disorder and material properties at the single-atom level. Nature 542, 75–79 (2017). https://doi.org/10.1038/nature21042

[2] Alex Zunger, S.-H. Wei, L. G. Ferreira, and James E. Bernard Special quasirandom structures Phys. Rev. Lett. 65, 353 – Published 16 July 1990

Coherent Potential Approximation (CPA) Method.

This method is implemented successfully in the first-principles matter modeling packages Questaal and EMTO.

You can find a general introduction about CPA on Wikipedia.

Answered by Jack on November 7, 2021

A very neat recent development for the modelling of disorder is the so-called localization landscape theory, first introduced here.

Consider a general Hamiltonian $$hat{H}=-frac{hbar^2}{2m}nabla^2+V$$. The landscape function $$u$$ is defined as the solution to this equation:

$$hat{H}u=1.$$

This innocent-looking landscape function has been shown to encode a great deal of information for studying the complex energy landscapes associated with disorder. For example, it can be shown that (i) the eigenstates $$psi(mathbf{r})$$ and eigenenergies $$E$$ of the Hamiltonian obey $$|psi(mathbf{r})|leq Eu(mathbf{r})$$, so that the localization landscape function $$u$$ delimits localization regions in space, or that (ii) the function $$W(mathbf{r})=1/u(mathbf{r})$$ acts as a confining potential that controls quantities like the exponential decay of Anderson localized states.

The localization landscape theory has recently been adopted in materials simulations of disorder. This series of three papers provides a very comprehensive introduction to the theory, and a demonstration of what is possible:

1. Basic theory and modelling: paper.
2. Urbach tails of disordered InGaN alloy quantum well layers: paper.
3. Carrier transport and recombination in light emitting diodes: paper.

Answered by ProfM on November 7, 2021

Better late than never for a partial answer.

I have worked with modeling the disorder of coverage on surfaces of materials. Recently we published a paper where we give a very general approach to modeling adsorbate-adsorbate interactions which we have seen is a problem in current literature. By using non-ideal buildups of coverage (where each adsorbate is not equally spaced), we have been able to identify better models of NO coverage on Pt3Sn. You can see below that at accessible high coverages, the structure is heavily disordered and would be incredibly difficult to find.

A large problem with modeling disorder is getting any useful information out of it. This problem needs to be solved before very disordered models will give better results than ideal models. In contrast to the comment by Camps, I don't feel making a good disordered model is difficult, the hard part is understanding what part of the disorder is actually relevant to the result.

Answered by Tristan Maxson on November 7, 2021

Related Questions

What are typical RPA capabilities of plane-wave codes?

1  Asked on August 19, 2021 by michael-f-herbst

Quadrature over three Euler Angles for orientation averaging

1  Asked on August 19, 2021 by emil-zak

What are some open-source all-electron DFT alternatives to Wien2K?

4  Asked on August 19, 2021

Real space projection vs reciprocal space projection in DFT calculations

1  Asked on August 19, 2021

How to calculate the Fock matrix in the molecular orbital basis PySCF?

1  Asked on August 19, 2021 by wychh

Bond Order: When and how is it used today?

2  Asked on August 19, 2021

Is it possible to calculate/estimate the value of the J parameter to be used in the Heisenberg/Ising Hamiltonians?

1  Asked on August 19, 2021

What are some methods for modeling bulk phase infrared spectra?

2  Asked on August 19, 2021

What is the closest thing we have to “the” universal density functional?

4  Asked on August 19, 2021

Molden installation on macOS?

1  Asked on August 19, 2021 by dmitry-eremin

What are the types of bond orders?

4  Asked on August 19, 2021

What is band inversion and how to recognize it in band structure?

1  Asked on August 19, 2021

What are the types of charge analysis?

10  Asked on August 19, 2021 by dgkang

What are the types of DFT?

5  Asked on August 19, 2021

Supercomputers around the world

14  Asked on August 19, 2021

What are the types of pseudopotentials?

2  Asked on August 19, 2021

What are the types of ab initio Molecular Dynamics?

3  Asked on August 19, 2021 by etienne-palos

Are there any MD packages that do proper free energy sampling in the NPT ensemble?

1  Asked on August 19, 2021 by b-kelly

VASP vs MedeA-VASP

1  Asked on August 19, 2021

Derivation of Slater-Koster equations

2  Asked on August 19, 2021