Firstly, I calculated bulk modulus of diamond by cp2k with QS METHOD, the lattice constant varied from 3.537 to 3.607, then the bulk modulus was calculated to 434.85 GPa by fitting the correlation between volume and energy (Birch–Murnaghan equation of state). And the result matches the real experiment.
Question1 : Why do their results have a big difference?(CP2K and LAMMPS) And is my input file wrong?
Question2 : I want to find a classical force field of diamond for QMMM, do you have a better advice?
The input files and results are attached.
Read more here: Source link