1
General Questions and Answers / Re: Formation energy and choice of reference states
« on: May 4, 2018, 15:43 »
Yes, thank you. I meant to say SG15 pseudo-potentials which are implemented in ATK. Thanks for the correction.
Show Posts
This section allows you to view all posts made by this member. Note that you can only see posts made in areas you currently have access to.
Pages: [1]
2
General Questions and Answers / Re: Formation energy and choice of reference states
« on: May 4, 2018, 15:30 »
Thank you for the quick reply! I am relatively new to DFT and therefore didn't know that the codes could be so different even with similar parameters.
To the point of question 2, I have only found cif files for example from the Crystallgraphic Open DB and other databases, no bulk configuration in a database in QW itself (leading to some ambiguity of which bulk structure to use). Otherwise, there is no bulk via File>Add>From Database... in the builder.
To the point of question 2, I have only found cif files for example from the Crystallgraphic Open DB and other databases, no bulk configuration in a database in QW itself (leading to some ambiguity of which bulk structure to use). Otherwise, there is no bulk via File>Add>From Database... in the builder.
3
General Questions and Answers / Formation energy and choice of reference states
« on: May 4, 2018, 14:16 »
Hello,
I am trying to calculate the energy of formation for a set of perovskites. I have calculated the total energy for each compound, and I understand that I need to now calculate the energies of the constituents in their standard phases, following dHformation = Etot - SUM(x*Etot,i), for each element i appearing x times in the composition.
These would typically be the chemical potentials of the elements in standard states, such as in the OQMD (or so I thought):
--> Elemental ground-state structures and chemical potentials predicted by DFT at 0 K
--> https://www.nature.com/articles/npjcompumats201510#t1
However the total energies output by ATK are very large (~20,000 eV for 20 atoms), and therefore do not give the correct results for dH. Following the tutorial on the subject (https://docs.quantumwise.com/tutorials/formation_energies/formation_energies.html), I see that for GaAs for example, one needs to calculate the bulk Ga and bulk As total energies. I therefore have 2 questions:
1. If the lowest energy phase should be chosen as the reference state, why can't the potentials in the above link be used? (i.e. why is the output so different than these values derived from VASP - I suspect this is about the PAW approach vs. GGA-PBE I am using).
2. How are the reference states to be chosen? I can download several different cif files of Ga and As phases from Crystallography Open DB. What reference states are actually used to get this number in the tutorial? (I have read the pointed to reference, and they also do not explicitly say how they calculate these numbers.)
It seems that the standard phases for all constituents must be calculated in the same ATK environment to keep the same DFT parameters (rather than compared with external values), and if this is true, I need to understand how to choose the correct phase. Thank you for the help!
I am trying to calculate the energy of formation for a set of perovskites. I have calculated the total energy for each compound, and I understand that I need to now calculate the energies of the constituents in their standard phases, following dHformation = Etot - SUM(x*Etot,i), for each element i appearing x times in the composition.
These would typically be the chemical potentials of the elements in standard states, such as in the OQMD (or so I thought):
--> Elemental ground-state structures and chemical potentials predicted by DFT at 0 K
--> https://www.nature.com/articles/npjcompumats201510#t1
However the total energies output by ATK are very large (~20,000 eV for 20 atoms), and therefore do not give the correct results for dH. Following the tutorial on the subject (https://docs.quantumwise.com/tutorials/formation_energies/formation_energies.html), I see that for GaAs for example, one needs to calculate the bulk Ga and bulk As total energies. I therefore have 2 questions:
1. If the lowest energy phase should be chosen as the reference state, why can't the potentials in the above link be used? (i.e. why is the output so different than these values derived from VASP - I suspect this is about the PAW approach vs. GGA-PBE I am using).
2. How are the reference states to be chosen? I can download several different cif files of Ga and As phases from Crystallography Open DB. What reference states are actually used to get this number in the tutorial? (I have read the pointed to reference, and they also do not explicitly say how they calculate these numbers.)
It seems that the standard phases for all constituents must be calculated in the same ATK environment to keep the same DFT parameters (rather than compared with external values), and if this is true, I need to understand how to choose the correct phase. Thank you for the help!
Pages: [1]