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.
that the number <strong>of</strong> configuration moves required to achieve a given error bar on the total energy scales<br />
as O(N); hence, in practice, the CPU time for ‘linear-scaling’ QMC calculations scales as O(N 2 ).<br />
27.2 Using <strong>CASINO</strong> to carry out ‘linear-scaling’ QMC calculations<br />
27.2.1 Generation <strong>of</strong> Bloch orbitals<br />
At present the Bloch orbitals must be represented in a plane-wave basis. A plane-wave DFT code<br />
should be used to generate a pwfn.data file as though an ordinary QMC calculation with a plane-wave<br />
basis were to be carried out. Note that only one k point may be used: the Γ point. This is not a<br />
severe restriction since, for large systems, one would usually only carry out calculations at Γ anyway.<br />
27.2.2 Generation <strong>of</strong> localized orbitals<br />
The localizer code should be used to carry out the linear transformation to a localized set <strong>of</strong> orbitals.<br />
The code has the following features:<br />
• localizer implements the method described in Ref. [70].<br />
• localizer requires a pwfn.data file holding the Bloch orbitals represented in a plane-wave basis<br />
and a centres.dat file <strong>of</strong> format:<br />
Number <strong>of</strong> centres<br />
<br />
Display coefficients <strong>of</strong> linear transformation (0=NO; 1=YES)<br />
<br />
Use spherical (1) or parallelepiped (2) localization regions<br />
<br />
x,y & z coords <strong>of</strong> centres ; radius ; no. orbs on centre (up & dn)<br />
<br />
...<br />
<br />
where ‘N’ is the total number <strong>of</strong> localization centres. If ‘iprint’ is 1 then the coefficients <strong>of</strong> the<br />
linear combination are written to stdout; otherwise they are not. ‘icut’ can take values 1 or<br />
2, specifying that the localization regions are spherical or parallelepiped-shaped, respectively.<br />
‘pos(i,j)’ is the ith Cartesian component <strong>of</strong> the position vector <strong>of</strong> the jth centre. ‘radius(j)’ is<br />
the cut<strong>of</strong>f radius for the jth localization centre. ‘norbs up(j)’ and ‘norbs dn(j)’ are the number<br />
<strong>of</strong> spin-up and spin-down orbitals to be localized on the jth centre, respectively. Note that if a<br />
parallelepiped-shaped localization region is used then the shape <strong>of</strong> the parallelepiped is defined<br />
by the lattice vectors, but the distance to each face is given by the cut<strong>of</strong>f radius.<br />
• The choice <strong>of</strong> localization centres requires some chemical intuition. See the discussion in Ref.<br />
[70]. At some point in the future, the optimization <strong>of</strong> the localization centres will be enabled.<br />
Note that in some highly symmetric molecules, symmetric choices <strong>of</strong> localization centres can<br />
lead to two localized orbitals being identical. This problem can be avoided by breaking the<br />
symmetry <strong>of</strong> the localization centres.<br />
• Note that the orbitals will become linearly dependent as two centres approach one another. In<br />
the limit that two centres are located in the same place, the orbitals localized on those centres<br />
will be identical. (One should instead define a single centre and increase the number <strong>of</strong> orbitals<br />
localized on that centre.)<br />
• localizer produces a pwfn.data.localized file that contains the localized orbitals represented<br />
in the same plane-wave basis as that specified in pwfn.data.<br />
• localizer can only work with Bloch orbitals at Γ. Since it is intended for use in large systems,<br />
this should not be a serious restriction.<br />
• One can only choose different numbers <strong>of</strong> localized orbitals for spin-up and spin-down electrons<br />
if pwfn.data contains spin-polarized data.<br />
167