Supported configuration files and their contents

Pycasino can read input files in well-known Casino format.

Input files are supported

  • input is the main input parameter file.

  • correlation.data file contains all optimizable parameters together with accompanying data (for example, the parameters used to define a Jastrow factor or backflow function).

  • gwfn.data file contains the data that defines the geometry, the gaussian-type orbitals and, if appropriate, the determinant expansion coefficients.

  • stowfn.data file contains the data that defines the geometry and the slater-type orbitals.

  • x_pp.data (where x is the chemical symbol of an element in lower-case letters.) This file contains the pseudopotential data for the corresponding element.

Input parameter are supported

General keywords

  • NEU, NED - number of electrons of up and down spin

  • RUNTYPE - type of QMC calculation: vmc, vmc_dmc, vmc_opt

  • TESTRUN - if this flag is T then read input files, print information and stop

  • ATOM_BASIS_TYPE - the type of orbitals to be used: gaussian, slater-type

VMC keywords

  • VMC_EQUIL_NSTEP - number of equilibration steps

  • VMC_NSTEP - number of VMC energy-evaluation steps

  • VMC_DECORR_PERIOD - number of steps between VMC energy-evaluation moves

  • VMC_NCONFIG_WRITE - number of VMC configurations stored for later use in DMC or optimization

  • VMC_NBLOCK - number of blocks into which the total VMC run is divided post-equilibration

  • DTVMC - VMC time step (size of trial steps in random walk)

  • VMC_METHOD - (1) - EBES (work in progress), (3) - CBCS.

Optimization keywords

  • OPT_CYCLES - number of optimization VMC cycles to perform

  • OPT_METHOD - optimization method to use: varmin, emin

  • OPT_JASTROW - optimize the Jastrow factor in wave-function optimization

  • OPT_BACKFLOW - optimize backflow parameters in wave-function optimization

  • OPT_DET_COEFF - optimize the coefficients of the determinants in wave-function optimization

  • OPT_MAXEVAL - maximum number of evaluations of the variance during variance minimization (default 50)

  • OPT_PLAN - allows specifying different parameters for each optimization cycle

  • VM_REWEIGHT - if set then the reweighted variance-minimization algorithm will be used, else the unreweighted algorithm will be used Unreweighted variance minimization is recommended

  • EMIN_METHOD - energy minimization method to use: newton, linear (default), reconf

DMC keywords

  • DMC_TARGET_WEIGHT - target number of configurations in DMC

  • DMC_EQUIL_NSTEP - number of DMC steps in equilibration

  • DMC_STATS_NSTEP - number of DMC steps in statistics accumulation

  • DMC EQUIL NBLOCK - number of blocks into which the DMC equilibration phase is divided

  • DMC STATS NBLOCK - number of blocks into which the DMC statistics accumulation phase is divided

  • DTDMC - DMC time step

  • DMC_METHOD - (1) - EBES, (2) - CBCS

  • LIMDMC - set modifications to Green’s function in DMC. Only (4) Umrigar mods to drift velocity, Zen–Sorella–Alfè mods to energy

  • ALIMIT - parameter required by DMC drift-velocity- and energy-limiting schemes

  • NUCLEUS_GF_MODS - this keyword is the switch for enabling the use of the modifications to the DMC Green’s function for the presence of bare nuclei

  • EBEST_AV_WINDOW - averaging window for calculating the ground-state energy during equilibration (work in progress)

WFN definition keywords

  • BACKFLOW - turns on backflow corrections. Backflow parameters are read from correlation.data

  • USE_JASTROW - use a wave function of the Slater-Jastrow form. The Jastrow factor is read from correlation.data

  • USE_GJASTROW - use gjastrow Jastrow factor. This Jastrow factor is defined in a parameters.casl file (work in progress)

Cusp correction keywords

  • CUSP_CORRECTION - when the cusp correction flag is activated, the s-type Gaussian basis functions centred on each atom are replaced within a small sphere by a function which ensures that the electron–nucleus cusp condition is obeyed

  • CUSP_INFO - if set then information about how cusp correction is done will be printed to the log-file

  • CUSP_THRESHOLD - if the magnitude of the s-component of a Gaussian orbital is less than this threshold then it will not be cusp corrected

Pseudopotential keywords

  • NON_LOCAL_GRID - selects the grid for nonlocal integration, can take values between 1 and 7, the default being 4

  • LCUTOFFTOL - this is used to define the cutoff radius for the local part of the pseudopotential, the default being 10-5

  • NLCUTOFFTOL - this is used to define the cutoff radius for the nonlocal parts of the pseudopotential, the default being 10-5

Configs read/write

Config loader is represented by the casino.readers.CasinoConfig class. The CasinoConfig instance represents the relevant information in various attributes:

  • config.input - data from input file

  • config.wfn - data from gwfn.data or stowfn.data files

  • config.mdet - multideterminant data from correlation.data files

  • config.jastrow - jastrow related data from correlation.data files

  • config.backflow - backflow related data from correlation.data files

Config files can be read:

from casino.readers import CasinoConfig

config_path = <path to a directory containing input file>
config = CasinoConfig(config_path)
config.read()

modified:

config.title = 'New config title'
config.mdet.title = 'New mdet section title'
config.jastrow.title = 'New jastrow section title'
config.backflow.title = 'New backflow section title'

and written to a new destination:

config.write('.', version=0)