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.
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