Calculation of 'Energy vs. lattice constant'

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=1

When 'MD.Type' is set to 'EvsLC', the total enery 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 $a_x$, $a_y$, $a_z$, $b_x$, $b_y$, $b_z$, $c_x$, $c_y$, $c_z$ in Å, and the total energy in Hartree, where $a(b,c)_x$, $a(b,c)_y$, and $a(b,c)_z$ 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. 41, 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 [83].

Figure 41: Total energy vs. lattice constant for the fcc Mn bulk calculated by the keyword 'EvsLC'. The input file used for the calculation is 'Mnfcc-EvsLC.dat' in the direcoty 'work'.