What's in the database?
|Full Heusler||Half Heusler||Inverse Heusler|
|A||Cr, Mn, Fe, and Co||Cr, Mn, Fe, Ru, Co, Rh, and Ni||Cu, Ni, Co, Fe, Mn, Cr, V, Ti, and Sc|
|B||Ti, V, Cr, Mn, and Fe||Ti, V, Cr, Mn, Fe, and Ni||Zn, Cu, Ni, Co, Fe, Mn, Cr, V, and Ti|
|C||Al, Ga, In, Si, Ge, Sn, P, As, and Sb||Al, Ga, In, Si, Ge, Sn, P, As, and Sb||Al, Ga, In, Si, Ge, Sn, P, As and Sb|
Details About DFT Calculations
The calculations were performed using the Vienna Ab-Initio Simulation Package (VASP) which uses a plane wave basis set and projector-augmented wave (PAW) based pseudo-potentials. A uniform cut-off energy of 520eV was chosen for all calculations. The Perdew-Burke-Ernzerhof (PBE) version of the generalized gradient approximation (GGA) to the exchange correlation functional of DFT was adopted. The pseudo-potential version for each element is listed below. These pseudo-potentials were chosen to be consistent with those used to construct the OQMD database of calculated energies of alloy phases.
|Zn: PAW_PBE Zn 06Sep2000 1|
|Cu: PAW_PBE Cu 05Jan200|
|Ni: PAW_PBE Ni 06Sep2000||Al: PAW_PBE Al 04Jan2001|
|Co: PAW_PBE Co 06Sep2000||Ga: PAW_PBE Ga_d 06Sep2000|
|Fe: PAW_PBE Fe 06Sep2000||In: PAW_PBE In_d 06Sep2000|
|Mn: PAW_PBE Mn 06Sep2000||Si: PAW_PBE Si 05Jan2001|
|Cr: PAW_PBE Cr 06Sep2000||Ge: PAW_PBE Ge_d 06Sep2000|
|V: PAW_PBE V 08Apr2002||Sn: PAW_PBE Sn_d 06Sep2000|
|Ti: PAW_PBE Ti 08Apr2002||P: PAW_PBE P 17Jan2003|
|Sc: PAW_PBE Sc_sv 07Sep2000||As: PAW_PBE As 06Sep2000|
|Ru: PAW_PBE Ru 06Sep2000||Sb: PAW_PBE Sb 06Sep2000|
|Rh: PAW_PBE Rh 06Sep2000|
The integration over the irreducible part of the Brillouin zone was done using the “automatic” 50 Monkhorst-Pack k-point grid and linear tetrahedron method with Blöchl corrections . The superlattices were relaxed without any constraint on shape or volume, using the conjugate-gradient algorithm. When ionic relaxations were performed, the stopping-criterion for structural optimization was an energy change less than 10-5 eV.
Determination of Lattice Constants
For a given A, B and C combination, we performed full ionic relaxations within a tetragonal cell for all of the Heusler alloys. Some of these alloys were found to remain in the cubic L21, C1b and XA structures, while some relaxed to a tetragonal structure; i.e. c/a is not equal to 1. Some of the half Heusler alloys were neither cubic nor tetragonal. All relaxations started from the cubic structure with small displacements to avoid vanishing of the net force on each atom due to symmetry.
Heusler alloys are prone to having multiple DFT solutions with different magnetic configurations. Therefore all energy versus lattice constant calculations were performed using multiple initial magnetic moment configurations, no moment, parallel moments and anti-parallel moments of the transition metal atoms.
The ground state of the Heusler alloy was determined by analyzing the results of the relaxation and cubic studies. For CoMnAl (Fig. 1), the Heusler alloy is stable when its lattice constants is 5.46 Å and the magnetic configuration of the A and B atoms is ferrimagnetic.
Formation Energy Per Atom
The formation energy per atom for each of the Heusler alloys was calculated by subtracting the energies of the individual elements in pure form from the alloy's total energy. For full Heusler A2BC, the formation energy per atoms is:
Eform = (EA2BC - 2EA - EB - EC)/4
where EA2BC is the total energy of the full Heusler unit cell containing 4 atoms (two A, one B, one C), EA, EB and EC stand for the energies when A, B, C, X, Y and Z elements crystallize in pure metals, respectively. If the formation energy per atoms is positive, the alloy is unstable and cannot be formed in equilibrium conditions. Negative formation energy per atom suggests that such an alloy may be grown. However, negative formation energy does not, for example, guarantee against formation of binary or ternary phases comprising any two or three of the available three elements. The reference states for the elements were chosen to be consistent with the OQMD database so that formation energies in this database can be compared to those of competing phases in the OQMD database.