CASINO manual - Theory of Condensed Matter
CASINO manual - Theory of Condensed Matter
CASINO manual - Theory of Condensed Matter
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
DTDMC (Real) Time step for DMC run (atomic units). The DMC time step must be small, as<br />
the DMC Green’s function is only exact in the limit <strong>of</strong> zero time step: see Sec. 13. Typically<br />
the DMC time step is about two orders <strong>of</strong> magnitude smaller than the VMC time step, and<br />
the DMC move-acceptance ratio should be about 99.9%. For accurate work, one must always<br />
investigate time-step bias by plotting the DMC energy against the value <strong>of</strong> dtvmc. Provided<br />
the time step is sufficiently small that the root-mean-square distance diffused by each particle at<br />
each step is much less than the shortest physically relevant length scale, one can expect to find<br />
the time-step bias in the DMC energy to be linear; hence it is straightforward to extrapolate<br />
to zero time step. The extrapolate tau utility exists to help with the extrapolation to zero<br />
time step.<br />
DTVMC (Real) Time step for VMC run (atomic units). The form <strong>of</strong> the VMC transition-probability<br />
density used by casino is discussed in Sec. 12. As described in Sec. 12.4, dtvmc should be<br />
chosen so that the overall move acceptance ratio is close to 50%. In most normal systems, the<br />
appropriate value is between 0.1 a.u. and 0.6 a.u. If dtvmc is given a sensible starting value<br />
(and, for normal systems, anything in this range is sensible), setting opt dtvmc to 1 will cause<br />
the time step to be optimized automatically. For very low-density systems, a much larger value<br />
<strong>of</strong> dtvmc is appropriate.<br />
DTVMCS (Block) Use this keyword to specify a VMC time step for each particle family explicitly,<br />
as well as to determine whether to optimize each <strong>of</strong> them individually. The contents <strong>of</strong> this<br />
block override the values <strong>of</strong> dtvmc and opt dtvmc. One line is to be written for each ‘family’<br />
<strong>of</strong> particles, the format <strong>of</strong> each line being:<br />
DTVMC OPT DTVMC<br />
where ‘DTVMC’ is the value <strong>of</strong> the time step, and ‘OPT DTVMC’ can be 0 or 1, indicating<br />
whether to optimize the corresponding time step or not.<br />
DTVMC SHIFT (Real) dtvmc shift is an optional shift in the VMC transition probability, which<br />
can be used to ‘encourage’ electrons to be more mobile. dtvmc shift is expressed in units <strong>of</strong><br />
the square root <strong>of</strong> dtvmc.<br />
E OFFSET (Physical) This keyword gives a constant shift E <strong>of</strong>fset in the total energy per electron<br />
such that the final result will be E = E calc − E <strong>of</strong>fset . The default is zero. This allows the user<br />
to add any constant contributions to the total energy that are not calculated within casino.<br />
EBEST AV WINDOW (Integer) Averaging window for calculating the ground-state energy during<br />
equilibration. During DMC equilibration the best estimate <strong>of</strong> the ground-state energy is taken<br />
to be the average local energy over the last ebest av window moves. The default <strong>of</strong> 25 is<br />
usually sufficient.<br />
EDIST BY ION, EDIST BY IONTYPE (Block) The edist by ion block allows fine control<br />
<strong>of</strong> the initial distribution <strong>of</strong> the electrons before equilibration starts. The standard algorithm<br />
shares out the electrons amongst the various ions weighted by the pseudo-charge/atomic number<br />
<strong>of</strong> the ion. Each electron is placed randomly on the surface <strong>of</strong> a sphere surrounding its parent<br />
ion. There are certain situations, for example a simple crystal with a very large lattice constant,<br />
where the standard algorithm in the points routine may give a bad initial distribution, which<br />
cannot be undone by equilibrating for a reasonable amount <strong>of</strong> time. This keyword allows a<br />
user-defined set <strong>of</strong> electron/ion associations to be supplied. The syntax is to supply N ion lines<br />
within the block which look like, e.g., 1 4 4, where the three numbers are: the ion sequence<br />
number; the number <strong>of</strong> up-spin electrons associated with this ion; the number <strong>of</strong> down-spin<br />
electrons associated with this ion. Alternatively one may use the edist by iontype keyword<br />
block, where you replace the ion sequence number with the ion type sequence number and the<br />
information is supplied only for each particular type <strong>of</strong> ion.<br />
EMIN MIN ENERGY (Physical) This keyword sets a minimum energy threshold for energy minimization,<br />
used to reject low-quality wave functions which produce spurious low VMC energy<br />
estimates. The value <strong>of</strong> emin min energy should (ideally) be set slightly below the groundstate<br />
energy. The ground-state energy is <strong>of</strong>ten not known, in which case a good estimate can<br />
sometimes be supplied. If this keyword is not set <strong>manual</strong>ly, casino will supply an automatic<br />
guess derived from the preceding VMC run. See Sec. 25.3.6 for more details.<br />
EMIN XI VALUE (Real) This keyword sets the value <strong>of</strong> the ξ parameter used to control semiorthogonalization<br />
in energy minimization. It should rarely be changed by users. See Sec. 25.3<br />
for details.<br />
41