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

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

PERMIT DEN SYMM (Logical) If this flag is set to T then he symmetry <strong>of</strong> the SCF charge<br />

density (in mpc.data) will be imposed on the QMC charge-density data used in the MPC<br />

interaction (with the justification that imposing an exact condition on the charge density can’t<br />

hurt) and also when writing the QMC density to expval.data. It is possible however that<br />

DMC will break the symmetry <strong>of</strong> the SCF calculation; in this case the user should turn <strong>of</strong>f<br />

permit den symm.<br />

PLOT BACKFLOW (Block) This block allows a plot <strong>of</strong> the backflow transformation to be made<br />

(see Sec. 23) just after VMC equilibration. The block should contain 2 lines, plus an optional<br />

line for plotting the Φ term: (1) ‘0’ or ‘1’ to (de-)activate this facility; (2) ‘kspin’, ‘knumber’<br />

and ‘zposition’; (3) value <strong>of</strong> fixed electron–nucleus distance r iI . This will produce various files:<br />

bfconfig.dat (reference config), bfconfigx.dat (associated quasi-particle config), bfions.dat<br />

(coordinates <strong>of</strong> nuclei for which backflow terms exist), bfeta 〈s〉.dat (η vs. r ij for each spinpair<br />

type 〈s〉), bfmu 〈s〉 〈set〉.dat (µ vs r i for each spin type 〈s〉 in each set 〈set〉), bfphi.dat<br />

[contribution <strong>of</strong> Φ to the 3D backflow displacement on electron j vs. 2D projection <strong>of</strong> r jI on the<br />

plane defined by electron j, ion I in set 〈set〉, and electron i at distance r iI from the nucleus with<br />

spin such that 〈s〉 is the spin-pair type <strong>of</strong> (i, j)], and bffield.dat (3D backflow displacement<br />

on electron (kspin,knumber) vs. its 2D position on the plane z = zposition).<br />

PLOT EXPVAL (Block) This block allows the plot expval utility to be told about the geometrical<br />

region over which a particular expectation value is to be plotted (which can be a line AB / plane<br />

AB-AC / volume AB-AC-AD). The structure <strong>of</strong> the block is as follows:<br />

Line 1: dimensionality d (1, 2 or 3);<br />

Line 2: No. <strong>of</strong> points along the d directions;<br />

Line 3: x, y and z coordinates <strong>of</strong> point A;<br />

Line 4: x, y and z coordinates <strong>of</strong> point B;<br />

Line 5: x, y and z coordinates <strong>of</strong> point C (if required);<br />

Line 6: x, y and z coordinates <strong>of</strong> point D (if required).<br />

The data will be plotted in a format suitable for xmgrace in the file lineplot.data or in<br />

a format suitable for gnuplot in 2Dplot.dat or 3Dplot.dat which can be quickly visualized<br />

using the plot 2D utility. This block is ignored by casino itself.<br />

POPSTATS (Logical) If popstats is T then the variance <strong>of</strong> the local energies sampled in DMC will<br />

be evaluated.<br />

POSTFIT VMC (Logical) If postfit vmc is set to T then a VMC calculation will be performed<br />

after the final optimization phase when runtype=‘vmc opt’ or ‘opt vmc’. This should normally<br />

be done.<br />

POSTFIT KEEP CFG (Logical) If postfit keep cfg is set to T then the configurations generated<br />

in the post-fit VMC calculation will be written to config.out.<br />

PRIMITIVE CELL (Block) Sometimes the ‘primitive lattice vectors’ in the xwfn.data file do not<br />

correspond to the true primitive cell. If the actual primitive lattice vectors are needed (e.g., for<br />

accumulation <strong>of</strong> the charge density), then override values can be supplied in the primitive cell<br />

block.<br />

PRINTGSCREENING (Logical) Before doing a periodic Gaussian calculation, casino prepares<br />

lists <strong>of</strong> potentially significant (primitive) cells and sites in each such cell which could contain<br />

Gaussians having a nonzero value in a reference primitive cell centred on the origin. Zero is<br />

defined as 10 −gautol . Turning on the printgscreening flag prints out the important information<br />

about this screening.<br />

PSI S (Text) If psi s is ‘none’ then the Slater wave function is set to one; if psi s is ‘slater’ (default)<br />

then an expansion in one or more Slater determinants is used; if psi s is ‘exmol’ then a specially<br />

crafted wave function for excitonic and positronic molecules is used.<br />

QMC DENSITY MPC (Logical) If this flag is set to T then the QMC charge data at the end <strong>of</strong><br />

the expval.data file will be used to compute the MPC interaction (see Sec. 19.4.4), rather than<br />

the default SCF density in the mpc.data file. This is likely to be useful in cases such as the<br />

Wigner crystal where the Hartree-Fock charge density is very different to the true charge density<br />

(it is too localized) as opposed to, say, the Fermi fluid where the Hartree-Fock charge density<br />

is exact. Note that using this option is likely to increase the time taken to evaluate the MPC<br />

interaction; in both DFT and QMC cases, the code counts backwards from the end <strong>of</strong> the list<br />

50

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

Saved successfully!

Ooh no, something went wrong!