1
|
Xie JZ, Zhou XY, Jin B, Jiang H. Machine Learning Force Field-Aided Cluster Expansion Approach to Phase Diagram of Alloyed Materials. J Chem Theory Comput 2024; 20:6207-6217. [PMID: 38940547 DOI: 10.1021/acs.jctc.4c00463] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/29/2024]
Abstract
First-principles approaches based on density functional theory (DFT) have played important roles in the theoretical study of multicomponent alloyed materials. Considering the highly demanding computational cost of direct DFT-based sampling of the configurational space, it is crucial to build efficient and low-cost surrogate Hamiltonian models with DFT accuracy for efficient simulation of alloyed systems with configurational disorder. Recently, the machine learning force field (MLFF) method has been proposed to tackle complicated multicomponent disordered systems. However, the importance of integrating significant physical considerations, including, in particular, convex hull preservation, which is the prerequisite for the accurate prediction of phase diagrams, into the training process of the MLFF remains rarely addressed. In this work, a workflow is proposed to train a convex-hull-preserved (CHP) MLFF for binary alloy systems, based on which the order-disorder phase boundary is predicted by using the Wang-Landau Monte Carlo (WLMC) technique. The predicted values for order-disorder phase transition temperatures agree well with the experiment. The CHP-MLFF is further used to build CE models with the same accuracy as the MLFF and higher efficiency in sampling configurational space. Using the results obtained from the MLFF-based WLMC simulation as a reference, the performances of different schemes for constructing CE models were evaluated in a transparent manner, which revealed the close correlation between the prediction accuracy of ground-state configurations and that of the order-disorder phase transition temperature. This work clearly indicates the great importance of reproducing the convex hull and energetics of ground-state configurations when constructing surrogate Hamiltonians for the statistical modeling of alloyed systems.
Collapse
Affiliation(s)
- Jun-Zhong Xie
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory of Rare Earth Material Chemistry and Application, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, 100871 Beijing, China
| | - Xu-Yuan Zhou
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory of Rare Earth Material Chemistry and Application, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, 100871 Beijing, China
| | - Bin Jin
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory of Rare Earth Material Chemistry and Application, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, 100871 Beijing, China
| | - Hong Jiang
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory of Rare Earth Material Chemistry and Application, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, 100871 Beijing, China
| |
Collapse
|
2
|
Harada JK, Charles N, Poeppelmeier KR, Rondinelli JM. Heteroanionic Materials by Design: Progress Toward Targeted Properties. ADVANCED MATERIALS (DEERFIELD BEACH, FLA.) 2019; 31:e1805295. [PMID: 30861235 DOI: 10.1002/adma.201805295] [Citation(s) in RCA: 90] [Impact Index Per Article: 18.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/14/2018] [Revised: 01/16/2019] [Indexed: 05/16/2023]
Abstract
The burgeoning field of anion engineering in oxide-based compounds aims to tune physical properties by incorporating additional anions of different size, electronegativity, and charge. For example, oxychalcogenides, oxynitrides, oxypnictides, and oxyhalides may display new or enhanced responses not readily predicted from or even absent in the simpler homoanionic (oxide) compounds because of their proximity to the ionocovalent-bonding boundary provided by contrasting polarizabilities of the anions. In addition, multiple anions allow heteroanionic materials to span a more complex atomic structure design palette and interaction space than the homoanionic oxide-only analogs. Here, established atomic and electronic principles for the rational design of properties in heteroanionic materials are contextualized. Also described are synergistic quantum mechanical methods and laboratory experiments guided by these principles to achieve superior properties. Lastly, open challenges in both the synthesis and the understanding and prediction of the electronic, optical, and magnetic properties afforded by anion-engineering principles in heteroanionic materials are reviewed.
Collapse
Affiliation(s)
- Jaye K Harada
- Department of Materials Science and Engineering, Northwestern University, Evanston, IL, 60208, USA
| | - Nenian Charles
- Department of Materials Science and Engineering, Northwestern University, Evanston, IL, 60208, USA
| | | | - James M Rondinelli
- Department of Materials Science and Engineering, Northwestern University, Evanston, IL, 60208, USA
| |
Collapse
|
3
|
Wong ZM, Tan TL, Yang SW, Xu GQ. Optimizing special quasirandom structure (SQS) models for accurate functional property prediction in disordered 2D alloys. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2018; 30:485402. [PMID: 30406769 DOI: 10.1088/1361-648x/aae764] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
2D materials such as MXenes have garnered attention in a wide field of applications ranging from energy to environment to medical. Properties of 2D materials can be tailored via alloying and in some cases, solid-solutions (disordered alloys) are formed. To predict the disordered alloy properties via first-principles, the model structure needs to imitate the random arrangements of alloyants and yet remains computationally tractable. Using density functional theory and the cluster expansion method, we investigate the accuracy of using of special quasirandom structures (SQSs) for predicting disordered 2D alloy properties, evaluating the effect of SQS supercell size on the prediction quality of formation energies, elastic properties, and structural parameters. We illustrate the findings with 5 different disordered binary [Formula: see text] MXene alloy systems (where M = Ti and M' = Zr, Hf, V, Nb, or Ta), demonstrating that SQSs around 6-8 times the primitive cell (N = 6-8) are sufficient to attain convergence in the property predictions versus supercell size. For formation energies, SQSs with N > 4 are found to reproduce the formation energies of the fully disordered phase within ~2.5 meV. For the simulation of the experimentally-synthesized TiNbCO2, we find convergence in structural parameters and elastic tensors at N ~ 6. We traced the convergence of the predictions to the convergence in the band structure-related properties via analysis of the electronic densities-of-states and the projected crystal overlap Hamilton population. Our findings suggest that modest sized SQSs would reproduce the properties of disordered MXene alloys. The results should help guide the investigations of structure-property relationships in other disordered 2D materials as well.
Collapse
Affiliation(s)
- Zicong Marvin Wong
- Department of Chemistry, National University of Singapore, 3 Science Drive 3, Singapore 117543, Singapore. Institute of High Performance Computing, Agency for Science, Technology and Research, 1 Fusionopolis Way, #16-16 Connexis, Singapore 138632, Singapore
| | | | | | | |
Collapse
|
4
|
D'Arco P, Mustapha S, Ferrabone M, Noël Y, De La Pierre M, Dovesi R. Symmetry and random sampling of symmetry independent configurations for the simulation of disordered solids. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2013; 25:355401. [PMID: 23912029 DOI: 10.1088/0953-8984/25/35/355401] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/02/2023]
Abstract
A symmetry-adapted algorithm producing uniformly at random the set of symmetry independent configurations (SICs) in disordered crystalline systems or solid solutions is presented here. Starting from Pólya's formula, the role of the conjugacy classes of the symmetry group in uniform random sampling is shown. SICs can be obtained for all the possible compositions or for a chosen one, and symmetry constraints can be applied. The approach yields the multiplicity of the SICs and allows us to operate configurational statistics in the reduced space of the SICs. The present low-memory demanding implementation is briefly sketched. The probability of finding a given SIC or a subset of SICs is discussed as a function of the number of draws and their precise estimate is given. The method is illustrated by application to a binary series of carbonates and to the binary spinel solid solution Mg(Al,Fe)2O4.
Collapse
Affiliation(s)
- Philippe D'Arco
- UPMC (Université Pierre et Marie Curie) Université Paris 6, ISTEP UMR 7193, F-75005, Paris, France. philippe.d
| | | | | | | | | | | |
Collapse
|
5
|
Mustapha S, D'Arco P, De La Pierre M, Noël Y, Ferrabone M, Dovesi R. On the use of symmetry in configurational analysis for the simulation of disordered solids. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2013; 25:105401. [PMID: 23388579 DOI: 10.1088/0953-8984/25/10/105401] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/01/2023]
Abstract
The starting point for a quantum mechanical investigation of disordered systems usually implies calculations on a limited subset of configurations, generated by defining either the composition of interest or a set of compositions ranging from one end member to another, within an appropriate supercell of the primitive cell of the pure compound. The way in which symmetry can be used in the identification of symmetry independent configurations (SICs) is discussed here. First, Pólya's enumeration theory is adopted to determine the number of SICs, in the case of both varying and fixed composition, for colors numbering two or higher. Then, De Bruijn's generalization is presented, which allows analysis of the case where the colors are symmetry related, e.g. spin up and down in magnetic systems. In spite of their efficiency in counting SICs, neither Pólya's nor De Bruijn's theory helps in solving the difficult problem of identifying the complete list of SICs. Representative SICs are obtained by adopting an orderly generation approach, based on lexicographic ordering, which offers the advantage of avoiding the (computationally expensive) analysis and storage of all the possible configurations. When the number of colors increases, this strategy can be combined with the surjective resolution principle, which permits the efficient generation of SICs of a problem in |R| colors starting from the ones obtained for the (|R| - 1)-colors case. The whole scheme is documented by means of three examples: the abstract case of a square with C(4v) symmetry and the real cases of the garnet and olivine mineral families.
Collapse
Affiliation(s)
- Sami Mustapha
- Institut de Mathématiques de Jussieu (UMR 7586 UPMC-CNRS), UPMC, Sorbonne Universités, Paris, France
| | | | | | | | | | | |
Collapse
|
6
|
Dzyabchenko AV, Reformatskaya II. Method of enumeration of crystal superstructures based on fcc, bcc, and hexagonal close packings. RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY A 2011. [DOI: 10.1134/s0036024412010098] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
7
|
Liu JZ, Zunger A. Thermodynamic theory of epitaxial alloys: first-principles mixed-basis cluster expansion of (In, Ga)N alloy film. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2009; 21:295402. [PMID: 21828531 DOI: 10.1088/0953-8984/21/29/295402] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
Epitaxial growth of semiconductor alloys onto a fixed substrate has become the method of choice to make high quality crystals. In the coherent epitaxial growth, the lattice mismatch between the alloy film and the substrate induces a particular form of strain, adding a strain energy term into the free energy of the alloy system. Such epitaxial strain energy can alter the thermodynamics of the alloy, leading to a different phase diagram and different atomic microstructures. In this paper, we present a general-purpose mixed-basis cluster expansion method to describe the thermodynamics of an epitaxial alloy, where the formation energy of a structure is expressed in terms of pair and many-body interactions. With a finite number of first-principles calculation inputs, our method can predict the energies of various atomic structures with an accuracy comparable to that of first-principles calculations themselves. Epitaxial (In, Ga)N zinc-blende alloy grown on GaN(001) substrate is taken as an example to demonstrate the details of the method. Two (210) superlattice structures, (InN)(2)/(GaN)(2) (at x = 0.50) and (InN)(4)/(GaN)(1) (at x = 0.80), are identified as the ground state structures, in contrast to the phase-separation behavior of the bulk alloy.
Collapse
|
8
|
d'Avezac M, Zunger A. Finding the atomic configuration with a required physical property in multi-atom structures. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2007; 19:402201. [PMID: 22049098 DOI: 10.1088/0953-8984/19/40/402201] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
In many problems in molecular and solid state structures one seeks to determine the energy-minimizing decoration of sites with different atom types. In other problems, one is interested in finding a decoration with a target physical property (e.g. alloy band gap) within a certain range. In both cases, the sheer size of the configurational space can be horrendous. We present two approaches which identify either the minimum-energy configuration or configurations with a target property for a fixed underlying Bravais lattice. We compare their efficiency at locating the deepest minimum energy configuration of face centered cubic Au-Pd alloy. We show that a global-search genetic-algorithm approach with diversity-enhancing constraints and reciprocal-space mating can efficiently find the global optimum, whereas the local-search virtual-atom approach presented here is more efficient at finding structures with a target property.
Collapse
|
9
|
Liu JZ, Trimarchi G, Zunger A. Strain-minimizing tetrahedral networks of semiconductor alloys. PHYSICAL REVIEW LETTERS 2007; 99:145501. [PMID: 17930682 DOI: 10.1103/physrevlett.99.145501] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/13/2007] [Indexed: 05/25/2023]
Abstract
The atomic size mismatch between different binary semiconductors has been long known to limit their mutual solubility, leading instead to phase separation into incoherent phases, forming inhomogeneous mixtures that severely limit technological applications that rely on carrier transport. We show here that this atomic size mismatch can lead, under coherent conditions, to the formation of a homogeneous alloy with characteristic (201) two-monolayer ordering. This occurs because such specific layer arrangement corresponds to a unique strain-minimizing network in tetrahedral systems.
Collapse
|
10
|
Chen S, Gong XG, Wei SH. Superhard Pseudocubic BC2N superlattices. PHYSICAL REVIEW LETTERS 2007; 98:015502. [PMID: 17358485 DOI: 10.1103/physrevlett.98.015502] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/28/2006] [Indexed: 05/14/2023]
Abstract
It is currently under debate whether diamondlike BC2N may be harder than cubic BN (c-BN). Using the bond counting rule, we have performed an unconstrained search and identified a series of short period (111) superlattices that have much lower total energy than previously proposed structures. By examining the ideal strength of these pseudocubic boron-carbonitrides, we show that they are harder than c-BN. Our results are consistent with experimental findings, but in contrast with a recent theoretical study, which claimed that the BC2N is less hard than c-BN.
Collapse
Affiliation(s)
- Shiyou Chen
- Surface Science Laboratory (National Key) and Physics Department, Fudan University, Shanghai, China
| | | | | |
Collapse
|
11
|
Sanati M, Wang LG, Zunger A. Adaptive crystal structures: CuAu and NiPt. PHYSICAL REVIEW LETTERS 2003; 90:045502. [PMID: 12570430 DOI: 10.1103/physrevlett.90.045502] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/12/2002] [Indexed: 05/24/2023]
Abstract
We discover that Au-rich Cu1-xAux and Pt-rich Ni1-xPtx contain a composition range in which there is a quasicontinuum of stable, ordered "adaptive structures" made of (001) repeat units of simple structural motifs. This is found by searching approximately 3x10(6) different fcc configurations whose energies are parametrized via a "cluster expansion" of first-principles-calculated total energies of just a few structures. This structural adaptivity is explained in terms of an anisotropic, long-range strain energy.
Collapse
Affiliation(s)
- M Sanati
- National Renewable Energy Laboratory, Golden, Colorado 80401, USA
| | | | | |
Collapse
|
12
|
Hart GL, Zunger A. Origins of nonstoichiometry and vacancy ordering in Sc1-x squarexS. PHYSICAL REVIEW LETTERS 2001; 87:275508. [PMID: 11800897 DOI: 10.1103/physrevlett.87.275508] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/13/2001] [Indexed: 05/23/2023]
Abstract
Whereas nearly all compounds A(n)B(m) obey Dalton's rule of integer stoichiometry (n: m, both integer), there is a class of systems, exemplified by the rocksalt structure Sc1-x squarexS, that exhibits large deviations from stoichiometry via vacancies, even at low temperatures. By combining first-principles total energy calculations with lattice statistical mechanics, we scan an astronomical number of possible structures, identifying the stable ground states. Surprisingly, all have the same motifs: (111) planes with (112) vacancy rows arranged in (110) columns. Electronic structure calculations of the ground states (identified out of approximately 3x10(6) structures) reveal the remarkable origins of nonstoichiometry.
Collapse
Affiliation(s)
- G L Hart
- National Renewable Energy Laboratory, Golden, Colorado 80401, USA
| | | |
Collapse
|
13
|
Mäder KA, Zunger A. Short- and long-range-order effects on the electronic properties of III-V semiconductor alloys. PHYSICAL REVIEW. B, CONDENSED MATTER 1995; 51:10462-10476. [PMID: 9977741 DOI: 10.1103/physrevb.51.10462] [Citation(s) in RCA: 24] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/12/2023]
|
14
|
Ceder G, Garbulsky GD, Avis D, Fukuda K. Ground states of a ternary fcc lattice model with nearest- and next-nearest-neighbor interactions. PHYSICAL REVIEW. B, CONDENSED MATTER 1994; 49:1-7. [PMID: 10009249 DOI: 10.1103/physrevb.49.1] [Citation(s) in RCA: 34] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/12/2023]
|