Dalton

From CUC3
Jump to navigation Jump to search

Dalton is an ab initio package installed on destiny. Comprehensive documentation can be found at [1].

AJWT has written a number of scripts to simplify running Dalton, and in particular extracting the four-index integrals required as input to NECI. These are available at destiny:~ajwt3/dalton

run.pl

run.pl [DALTON]
Parse an input file substituting a placeholder for a series of values, and then call ./submit.sh.2 for each value with the appropriate specification for dirroot (placeholders changed etc.) and queue. tag is a tag to append to the output filename. submit.sh.2 runs Dalton, and then the postprocessing scripts to produce an FCIDUMP file. As well as the file.mol containing the placeholder, the filename can also contain it.

Usage:
run.pl file.mol placeholder start end step dirroot tag frozen queue
or
run.pl file.mol placeholder listofvalues 0 0 dirroot tag frozen queue


molrun.pl

molrun.pl [DALTON]
Perform a series of runs for a set of molecules specified in a MOLECULES file.

Usage:
./molrun.pl BASIS MOLLIST FROZEN QUEUE QSIZE

DALTEMPLATE
specifies a template .dal file
BASIS
specification of basis (e.g. cc-pVDZ). Must be as specified for MOLPRO
TAG
tag for root of directory and FCIDUMP FILE
MOLLIST
optional list of specific molecules rather than all. e.g. "HOH NN"
FROZEN
0 if not freezing cores, or 1 if freezing cores
QUEUE
queue. If queue has a q prepended, jobs are batched.

MOLECULES file

A file containing a list of molecules and geometries to run. Routines to read and interpret it are in mols.pm. The format is as follows

HEADER

MOLECULE1

MOLECULE2

...

HEADER

The header contains a list of input parameters to the script, which may be the following:

SCALE=s1,s2
This gives two possible values for a scale parameter, s. All distance measurements are multiplied by s. If s1=s2 then only the single value is used.
DISTVAR=dpc,npts
This allows a scan over distance variables. For each distance variable, d, choose npts points ranging from <math>d*(1-dpc)</math> to <math>d*(1+dpc)</math>.
ANGVAR=apc,npts
This allows a scan over angle variables. For each angle variable, a, choose npts points ranging from <math>a*(1-apc)</math> to <math>a*(1+apc)</math>.
LSHS
If present, use low-spin values at scale s1, and high-spin values at scale s2.

MOLECULEn

Each molecule specification consists of the following lines:

NAME NEL NCORE HSSPIN HSSYM
ZMAT1;ZMAT2;ZMAT3;...;ZMATn
var1=val1 [;CALC]
...
varn=valn [;CALC]

These have the following meanings

NAME
A moniker given to the molecule. This is used to tag all FCIDUMP files and run directories.
NEL
The total number of electrons in the molecule
NCORE
The number of core electrons in the molecule (used when frozen-core calculations are performed)
HSSPIN
The value of 2*Ms for a high-spin state (used with LSHS) (untested with Dalton])
HSSYM
The MOLPRO symmetry specification for a high-spin state (used with LSHS) (untested with Dalton)
ZMATn
The nth line of the Z-matrix for the molecule. This is comma-separated, and can contain variables and references to atoms on previous lines.
varn=valn [;CALC]
Variable definitions. Distance variables begin with r, and angular begin with a. valn is an evaluable expression comprising of previously defined variables. All distance and angle variable are liable to a scan if DISTVAR or ANGVAR are specified in the header except if ;CALC is appended to the definition. Distances must be specified in Angstroms.