14.09.2014 Views

CASINO manual - Theory of Condensed Matter

CASINO manual - Theory of Condensed Matter

CASINO manual - Theory of Condensed Matter

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!