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.

OPT MAXITER (Integer) Largest permitted number <strong>of</strong> nl2sol or global EMIN iterations (default:<br />

10).<br />

OPT METHOD (Text) There are currently four optimization methods implemented in casino:<br />

(1) variance minimization (‘varmin’); (2) minimization <strong>of</strong> the mean absolute deviation <strong>of</strong> the<br />

set <strong>of</strong> local energies from the median (‘madmin’); (3) energy minimization (‘emin’); (4) an<br />

accelerated variance minimization technique for parameters that appear linearly in the Jastrow<br />

factor (‘varmin linjas’). The first three methods are capable <strong>of</strong> optimizing the Jastrow factor,<br />

orbitals, backflow functions, and determinant coefficients. The fourth method can only be used<br />

to optimize linear parameters in the Jastrow factor. There are other keywords that affect the<br />

behaviour <strong>of</strong> each <strong>of</strong> these methods. The default value <strong>of</strong> opt method is ‘varmin’.<br />

OPT NOCTF CYCLES (Integer) Supplying a positive integer X for this keyword will cause all<br />

cut<strong>of</strong>f parameters (in the Jastrow factor, backflow function and orbital functions) to be fixed for<br />

the final X cycles <strong>of</strong> an optimization. Defaults to 0 for variance minimization, and to 50% <strong>of</strong><br />

total cycles for energy minimization.<br />

OPT ORBITALS (Logical) Optimize parameters in the orbitals in wave-function optimization (relevant,<br />

e.g., if use orbmods=T).<br />

OPT STRICT (Logical) Setting opt strict=T will cause casino to stop a vmc opt or opt vmc run<br />

if the VMC energies are incremented within a 99.7% confidence interval during two consecutive<br />

cycles. Prevents wastage <strong>of</strong> CPU time in times <strong>of</strong> scarcity. Default value is F.<br />

ORB NORM (Real) Allows user to change normalization <strong>of</strong> orbitals by multiplying all <strong>of</strong> them by<br />

this constant. Of course this should have no effect on the energy, but it can be useful if the<br />

Slater determinant starts going singular, as it might for some very low-density systems.<br />

ORBBUF (Logical) Setting orbbuf=T turns on orbital buffering in DMC. This is an efficiency<br />

device in which buffered copies <strong>of</strong> orbitals/gradients/Laplacians are kept for later reuse. This<br />

has a significant memory cost. Orbital buffering should always be used unless you start running<br />

out <strong>of</strong> memory; hence the ability to turn it <strong>of</strong>f.<br />

PAIR CORR (Logical) Set pair corr to T to accumulate the reciprocal-space pair-correlation function<br />

in the expval.data file. Currently restricted to periodic systems. Note that you also need<br />

to give the position and type <strong>of</strong> a fixed particle using the pcf rfix block (unless the density is<br />

homogeneous). See Sec. 33.<br />

PAIR CORR SPH (Logical) If pair corr sph is set to T then the spherically averaged real-space<br />

pair-correlation function will be accumulated in the expval.data file (via a process <strong>of</strong> ‘binning’<br />

the electron-electron separations). This currently works for periodic homogeneous systems and<br />

finite isotropic systems such as electron-hole biexcitons. For periodic systems with atoms you<br />

can use the pair corr keyword instead which gives you the full (non-spherically averaged) paircorrelation<br />

function accumulated in reciprocal space. See Sec. 33.<br />

PARTICLES (Block) Using the particles block the user can define quantum particles (other than<br />

electrons, which can be introduced using neu and ned) to be used in the QMC calculation.<br />

The format <strong>of</strong> each line is ‘〈i〉 〈charge/|e|〉 〈mass/m e 〉 〈spin/¯h〉 〈name〉’. A negative value <strong>of</strong><br />

the mass indicates that the following three lines give an anisotropic 3 × 3 mass tensor [currently<br />

unused]. casino decides whether each particle type is a fermion or a boson (based on the spin),<br />

and selects the appropriate way to combine the one-particle orbitals (symmetric combination<br />

[not currently implemented] or antisymmetric Slater determinants). The particles defined here<br />

can be assigned orbitals using the free particles block.<br />

PCF RFIX (Block) This block contains two lines. The first line gives the type <strong>of</strong> particle to be<br />

fixed during accumulation <strong>of</strong> the pair correlation function g(r, r ′ ); the second line gives the<br />

coordinates <strong>of</strong> the position r at which to fix it (in a.u.). This applies to the reciprocal-space<br />

PCF activated with the pair corr input keyword. It also applies in principle to the spherical<br />

real space PCF activated with the pair corr sph input keyword, in the sense that the format<br />

<strong>of</strong> expval.data allows it, but the accumulation <strong>of</strong> the spherical PCF with fixed particles has<br />

not yet been implemented. See Sec. 33.<br />

PERIODIC (Logical) T if and only if the system is periodic in either 1, 2 or 3 dimensions.<br />

49

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

Saved successfully!

Ooh no, something went wrong!