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.
Std. err. in the mean DMC energy (a.u.)<br />
0.0005<br />
0.0004<br />
0.0003<br />
0.0002<br />
0.0001<br />
0<br />
0 5 10 15 20<br />
Reblocking transformation number<br />
The file config.out contains the final positions <strong>of</strong> the electrons and the current state <strong>of</strong> the random<br />
number generator so that the VMC run can be continued if desired. To do this, set the input keyword<br />
newrun to F and rename the config.out file to config.in (in fact the runqmc script will normally<br />
ask if it can do this for you), then run the calculation for another vmc nstep steps. All the extra data<br />
will be put onto the end <strong>of</strong> vmc.hist and the error bar from the reblocking analysis will be smaller<br />
(in principle the error bar can be made as small as desired). As I already mentioned, you need to use<br />
the reblock utility for continued runs since the on-the-fly reblocking will no longer work.<br />
Finally, go and find a system with more than one electron and a Jastrow factor, such as the beryllium<br />
dimer in ~/<strong>CASINO</strong>/examples/molecule/be2/. Experiment with switching use jastrow between F<br />
(HFVMC) and T (VMC with a Jastrow factor), and see the energy- and error-bar-lowering effect <strong>of</strong><br />
the Jastrow factor (results shown below for 12 cores).<br />
Without Jastrow:<br />
E = -29.107(8) ; var = 6.6(6)<br />
Total <strong>CASINO</strong> CPU time ::: 4.0744 seconds<br />
With Jastrow:<br />
E = -29.215(2) ; var = 0.303(2)<br />
Total <strong>CASINO</strong> CPU time ::: 6.0221 seconds<br />
After any casino calculation, you can delete all the output that casino produces and restore the<br />
directory to the state it was in at the start <strong>of</strong> the calculation by typing clearup.<br />
6.3 Wave function optimization<br />
The values for the parameters in the Jastrow factor used in the last part <strong>of</strong> the previous section (and<br />
other adjustable parameters such as the coefficients in a multideterminant expansion, or backflow<br />
parameters) generally need to be optimized so that we have best possible trial wave function <strong>of</strong> the<br />
given functional form. There are various methods for doing this, including minimization <strong>of</strong> the variance(‘varmin’),<br />
minimization <strong>of</strong> the energy (‘emin’), and minimization <strong>of</strong> the mean absolute deviation<br />
<strong>of</strong> the local energies from the median (‘madmin’). Wave-function optimization is probably the most<br />
difficult part <strong>of</strong> QMC for beginners, and perseverance will help. The user might first like to take a<br />
look at Sec. 25.1, which contains a detailed summary <strong>of</strong> the theory and best practice. Here we simply<br />
summarize.<br />
17