CASINO manual - Theory of Condensed Matter
CASINO manual - Theory of Condensed Matter
CASINO manual - Theory of Condensed Matter
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
for all-electron calculations scales rather badly with atomic number Z; this scaling is improved by<br />
using pseudopotentials.<br />
casino reads pseudopotential data from a file called xx pp.data, where xx is the chemical symbol<br />
<strong>of</strong> the element in lower-case letters. The file contains a logarithmic radial grid and the values <strong>of</strong> the<br />
different angular momentum components <strong>of</strong> the pseudopotential at each grid point.<br />
On the rare occasions when you might want to use two or more different pseudopotentials for atoms<br />
with the same atomic number (say in a surface, and in an atom or molecule absorbed on that surface),<br />
then you may use additional pseudopotentials renamed as, e.g., xx2 pp.data. Different types <strong>of</strong><br />
pseudoatom are flagged in xwfn.data by adding multiples <strong>of</strong> 1000 to the original atomic number,<br />
e.g., atomic numbers 12, 1012 and 2012 correspond to atoms using pseudopotentials mg pp.data,<br />
mg2 pp.data and mg3 pp.data, etc.<br />
A library <strong>of</strong> pseudopotentials suitable for use with casino is available at:<br />
http://vallico.net/casinoqmc/pplib/<br />
An important point is that exactly the same pseudopotential should be used in the DFT/HF calculation<br />
that generates the trial wave function and in casino. Other programs do not in general understand<br />
the casino pseudopotential format, and so the information must somehow be transformed so<br />
that they do.<br />
For programs using Gaussian basis sets such as gaussian9X/0X and crystal, the pseudopotential<br />
must be re-expanded in Gaussian functions multiplied by powers <strong>of</strong> r. If you are using the Cambridge<br />
pseudopotentials, such expansions (in formats suitable for these two programs) are included in the<br />
on-line library. There is a casino utility—ptm—which can manipulate pseudopotential files on grids<br />
and their Gaussian expansions in various useful ways.<br />
As for plane-wave programs, castep understands the casino grid format and can read such<br />
files directly. Other plane-wave programs require conversion utilities, which are included in the<br />
~/<strong>CASINO</strong>/utils/pseudo converters/ directory.<br />
For atomic calculations on radial grids, awfn.data files generated with the Cambridge pseudopotentials<br />
are available in the on-line library.<br />
6.1.3 The input file<br />
Having prepared a trial wave-function file and (perhaps) a pseudopotential file, you need to tell casino<br />
exactly what to do with them. casino takes its instructions from a file called input, which contains<br />
a flexible list <strong>of</strong> keywords. These control the behaviour <strong>of</strong> the calculation, switch on and <strong>of</strong>f various<br />
options and so on. Take a moment to examine the various input files in ~/<strong>CASINO</strong>/examples/ to get<br />
a feel for what they look like.<br />
A complete list <strong>of</strong> input keywords, together with their definitions, is given in Sec. 7.3. Further details,<br />
including default values, may be found by using the casinohelp utility. This tends to be more up to<br />
date than the <strong>manual</strong>, since it interrogates casino directly. Type casinohelp all to get a list <strong>of</strong> all<br />
keywords that casino knows about, or casinohelp keyword for detailed help on a particular keyword.<br />
Type casinohelp search text to search for the string ‘text’ in all the keyword descriptions.<br />
Although there are many keywords, the beginner can play around by changing only a few <strong>of</strong> them.<br />
Here’s a (very) rough guide. Advice on good values to use will be given in the subsequent sections<br />
explaining how to do VMC, DMC and optimization calculations.<br />
• General (system-dependent) keywords:<br />
NEU, NED Number <strong>of</strong> electrons <strong>of</strong> up and down spin;<br />
PERIODIC Whether the system is periodic or not;<br />
NPCELL Array <strong>of</strong> primitive cells making up the simulation cell (not required for finite systems).<br />
• Other vital keywords:<br />
RUNTYPE Type <strong>of</strong> QMC calculation: ‘vmc’, ‘vmc dmc’, ‘vmc opt’, etc.;<br />
12