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.

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

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

Saved successfully!

Ooh no, something went wrong!