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.

ATOM BASIS TYPE The basis used to represent the orbitals: ‘plane-wave’, ‘gaussian’,<br />

‘slater-type‘, ‘blip’, ‘numerical’, ‘dimer’ or ‘none’ (the last option is for HEGs, etc.);<br />

• Important VMC keywords:<br />

VMC EQUIL NSTEP Number <strong>of</strong> equilibration steps;<br />

VMC NSTEP Number <strong>of</strong> VMC energy-evaluation steps;<br />

VMC DECORR PERIOD Number <strong>of</strong> steps between VMC energy-evaluation moves.<br />

VMC NCONFIG WRITE Number <strong>of</strong> VMC configurations stored for later use in DMC or<br />

optimization;<br />

DTVMC VMC time step (size <strong>of</strong> trial steps in random walk);<br />

• Important optimization keywords:<br />

OPT CYCLES Number <strong>of</strong> optimization+VMC cycles to perform.<br />

OPT METHOD Optimization method to use: variance minimization (‘varmin’), energy minimization<br />

(‘emin’), etc.<br />

• Important DMC keywords:<br />

DMC TARGET WEIGHT Target number <strong>of</strong> configurations in DMC;<br />

DMC EQUIL NSTEP Number <strong>of</strong> DMC steps in equilibration;<br />

DMC STATS NSTEP Number <strong>of</strong> DMC steps in statistics accumulation;<br />

DTDMC DMC time step.<br />

6.1.4 Correlation parameter file<br />

If you run a casino VMC calculation using a trial wave function consisting <strong>of</strong> only a single determinant<br />

<strong>of</strong> orbitals (referred to as an ‘HFVMC’ calculation), then the result will be the HF energy. If the<br />

orbitals were generated using a HF calculation, then the HFVMC energy should agree with the<br />

HF energy from the generating code. This is a good check that everything is being done correctly.<br />

Obviously, if the determinant is made up <strong>of</strong> Kohn–Sham orbitals from a DFT calculation then the total<br />

energies will not agree, because the DFT program adds an XC energy deduced from the self-consistent<br />

charge density; however, the kinetic energies should still be in agreement.<br />

The full Slater-Jastrow trial function normally used in QMC requires the determinantal part <strong>of</strong> the<br />

wave function stored in xwfn.data to be multiplied by a separate ‘Jastrow factor’, which defines the<br />

functional form <strong>of</strong> explicit interparticle correlations. In a typical VMC calculation one might recover<br />

60–80% <strong>of</strong> the correlation energy using such a wave function. This is not really enough to be generally<br />

useful, and in practice the main use <strong>of</strong> VMC is to prepare an accurate trial wave function to be given<br />

as input to a DMC calculation. The DMC energy does not in principle depend on the Jastrow factor,<br />

since the Jastrow factor is positive definite and the DMC energy depends only on the nodal surface<br />

(the set <strong>of</strong> points in configuration space where the many-electron wave function is zero). However, it<br />

makes the calculation vastly more efficient, and in general Jastrow factors should always be used.<br />

The Jastrow factor is stored in a file called correlation.data. Again, you should look at the<br />

examples to see what these look like. The various parameters in the files are defined in Sec. 7.4.2.<br />

The adjustable parameters in the file must be optimized for a specific system, and this is the purpose<br />

<strong>of</strong> the variance-minimization procedure.<br />

The correlation.data files may also contain other optimizable parameters not contained in the<br />

Jastrow factor (for example, ‘backflow’ parameters, or the coefficients in a determinant expansion)<br />

but you don’t need to know about these yet.<br />

6.1.5 The MPC data file<br />

If, for a two- or three-dimensionally periodic system, you want to use the model periodic Coulomb<br />

(MPC) interaction to calculate the electron–electron energies instead <strong>of</strong> (or as well as) the standard<br />

Ewald interaction, then you need to generate an extra file for the given geometry before you start<br />

doing VMC/DMC calculations. This is called mpc.data and contains the Fourier transform <strong>of</strong> the 1/r<br />

interaction and the Fourier transform <strong>of</strong> the charge density. The mpc.data file should be prepared<br />

13

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

Saved successfully!

Ooh no, something went wrong!