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.

 

lattice constraints
 

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.