Cohesive energy of Albite from Materials Project for LAMMPS comparison


I am doing some work using LAMMPS to calculate the cohesive energy of various silicates. I notice when I do albite I am getting a much lower cohesive energy than expected so I am looking for some kind of a comparison.

Does materials project have the cohesive energy of albite listed? I see final energy/atom, but am not sure if this is the cohesive energy? Further, should cohesive energy be for the mineral (all atoms in the system), or for each individual atom type in the system. For example, the cohesive energy of Na in albite is …

If anyone has DFT data or materials project data for the cohesive energy of albite it would be much appreciated.



Mat Proj has pure albite if that’s what you’re looking for: mp-721988: NaAlSi3O8 (triclinic, P-1, 2)

1 Like

Thanks, but where do I actually find the cohesive energy?

It is not directly listed.

E_{cohesive} = E_{bulk} - \sum_{i \in bulk} n_i E_{i}

Where n is the number of atoms of type i in the bulk structure, and E_i is the energy of isolated atom i. Usually if you want to analyze phase stability, you can use the formation energy, which defined with the pure elemental crystals as the reference states in MP. If you want to calculate the cohesive energy itself, you should use VASP to calculate the isolated atom energies. You can find these values listed in the pseudopotential file, which we cannot give out publically, you need to have a vasp license.

You might want to check the MPDS: MPDS materials platform however we don’t have that property from ab initio just yet (only ML). ETA is August 2021.

MP hosts a collection of calculated isolated-atom energies and facilitates calculation of cohesive energy given a material ID via the pymatgen.ext.matproj.MPRester.get_cohesive_energy method. You can see from the linked method code how to directly fetch the isolated-atom calculation results if you wish to compute cohesive energy in an alternative way. We never quite got these up on the website. Example:

from pymatgen.ext.matproj import MPRester
mpr = MPRester()
mpr.get_cohesive_energy('mp-721988', per_atom=True) # 6.73 eV/atom

Hi all,

Thanks very much for the help. I will go through these and message back if it isn’t working.