The calculation of Energy vs. lattice constant is supported by the following keywords:
MD.Type EvsLC # MD.EvsLC.Step 0.4 # default=0.4% MD.maxIter 32 # default=1When 'MD.Type' is set to 'EvsLC', the total energy is calculated step by step by changing unit cell vectors, a, b, and c. The change of unit cell vectors is done uniformly by expanding them by a percentage, where the reference is the initial vectors, specified with 'MD.EvsLC.Step'. The number of steps is specified by the keyword 'MD.maxIter'. After the calculation, you will obtain a file '*.EvsLC', where '*' is 'System.Name'. The columns in the file '*.EvsLC' are arranged in order of , , , , , , , , in Å, and the total energy in Hartree, where , , and are x-, y-, and z-coordinates of the a(b,c) vector, respectively. As an example, calculation of Energy vs. lattice for the fcc Mn bulk is shown in Fig. 42, where the equilibrium lattice constant and bulk modulus were evaluated by fitting the data to the Murnaghan equation of state with a code 'murn.f' provided on the web site [88].