next up previous
Next: Examples of input command Up: users_guide Previous: users_guide

Description of the command structure

The input command file consists of separate lines each containing

-
a label
-
a label followed by a character, integer and/or fix-point or
-
a character, integer and/or fix-point constant.

Note that

-
labels can be in upper or lower case,

-
compulsory labels, i.e. these that must be included in the input command file, are marked ``$ \bullet$'' and the optional ones ``$ \circ$'',

-
the compulsory labels must follow the order given below; the optional ones can be inserted anywhere between title and stop labels,

-
optional parameters are enclosed in square brackets,

-
$ r$ denotes a fix-point constant, $ i$ - an integer, $ c$ - a string of characters,

-
exclamation mark placed anywhere in an input line starts a comment; what follows "!" is ignored.

The following labels can be used:

 $&bull#bullet;$
TITLE
Format:
title
$ c$
$ c$ is any string of up to 80 characters describing the current case. This string is used as a header of disk output files containing orbitals and potentials.

 $&bull#bullet;$
METHOD
Format:
method $ c$
Select the type of calculation.
$ c$:
HF - Hartree-Fock method
$ c$:
HFS - Hartree-Fock-Slater method
$ c$:
OED - One Electron Diatomic states

 $&bull#bullet;$
NUCLEI
Format:
nuclei $ Z_A \;\;Z_B\;\; R \;\;[\; c\;]$
Set the nuclei charges and the bond length.
$ Z_A$:
nuclear charge of centre A (integer)
$ Z_B$:
nuclear charge of centre B (integer)
$ R$:
bond length (real)
$ c$:
$ angstrom$ - the internuclear separation can be given in angstrom units if this string is included (the conversion factor 0.529177249 is used)
If $ \vert Z_A-Z_B\vert<10^{-6}$ then the molecule is considered to be a homonuclear one (this threshold can be changed by redefining HOMOLEVL variable in blk_data.inc.raw).

$ \circ$
HOMO
Format:
homo
This label is used to impose explicitly $ D_{\infty h}$ symmetry upon orbitals of homonuclear molecules in order to improve scf/sor convergence.

$ \circ$
BREAK
Format:
break
When this label is present homonuclear molecules are calculated in $ C_{\infty v}$ symmetry and the $ D_{\infty h}$ symmetry labels (u or g) are superfluous (see below).

 $&bull#bullet;$
CONFIG
Format:
config $ i$
$ i$:
the total charge of a system

The following cards define molecular orbitals and their occupation. Note that the last orbital description card must contain an extra label $ end$.

The possible formats are:

Format:
$ i\;\;\;c$
$ i$:
number of fully occupied orbitals of a given irreducible representation (irrep) of the $ C_{\infty v}$ group; 2 electrons make $ \sigma$ orbitals fully occupied and 4 electrons are needed for filling the orbitals belonging to the other irreps
$ c$:
symbol of the $ C_{\infty v}$ irrep to which the orbitals belong (sigma, pi, delta or phi)

Format:
$ i\;\;\;c_1 \;\;c_2$
$ i$:
number of fully occupied orbitals of a given irrep of the $ D_{\infty h}$ group
$ c_1$:
symbol of the $ C_{\infty v}$ irrep to which the orbitals belong (sigma, pi, delta or phi)
$ c_2$:
symbol for the inversion symmetry of the $ D_{\infty h}$ irreps (u or g)
Use this format for a homonuclear molecule unless break command is included.

Format:
$ i\;\;\;c_1\;\;c_2\;\;[c_3\;\;[\;c_4\;\;[\;c_5\;]\;]\;]$
$ i$:
number of orbitals of a given irrep of the $ C_{\infty v}$ group
$ c_1$:
symbol for the $ C_{\infty v}$ irreps to which the orbitals belong (sigma, pi, delta, phi)
$ c_2$-$ c_5$:
$ +,\;-$ or . (a dot); $ +/-$ denotes spin up/down electron and . denotes an unoccupied spin orbital

Format:
$ i\;\;\;c_1\;\;c_2\;\;c_3\;\;[\;c_4\;\;[\;c_5\;\;[\;c_6\;]\;]\;]$
$ i$:
number of orbitals of a given irrep of the $ D_{\infty h}$ group
$ c_1$:
symbol for the $ C_{\infty v}$ irrep to which the orbitals belong (sigma, pi, delta, phi)
$ c_2$:
symbol for the inversion symmetry of the $ D_{\infty h}$ irrep (u or g)
$ c_3$-$ c_6$:
$ +,\;-$ or . (a dot); $ +/-$ denotes spin up/down electron and . denotes an unoccupied spin orbital

$ \circ$
GRID

Two possible formats are (the second one is retained for the backward compatibility):

Format:
grid $ n_{\nu}$ $ R_{\infty}$
An integer and a real define a single 2d grid.
$ n_{\nu}$:
the number of grid points in $ \nu$ variable
$ R_{\infty}$:
the practical infinity
$ n_{\mu}$ is calculated so as to make the step size in $ \mu$ variable equal to the stepsize in $ \nu$ variable. $ n_{\nu}$ and $ n_{\mu}$ have to meet special conditions. If the conditions are not fulfilled the nearest (but smaller) appropriate values are used.

Format:
grid $ n_{\nu}$ $ n_{\mu}$ $ R_{\infty}$
Two integers and one real define a single 2d grid.
$ n_{\nu}$:
the number of grid points in $ \nu$ variable
$ n_{\mu}$:
the number of grid points in $ \mu$ variable
$ R_{\infty}$:
the practical infinity
$ n_{\nu}$ and $ n_{\mu}$ have to meet special conditions. If the conditions are not fulfilled the nearest (but smaller) appropriate values are used.

$ \circ$
SUBGRID
Format:
subgrid $ n_{grids}$
$ n_{\nu}$ $ n_{\mu}^{(1)} \;[\; n_{\mu}^{(2)} \;
[\;n_{\mu}^{(3)}\;]\;]$
$ h_{\mu}^{(1)} \;[\; s^{(2)}[\; s^{(3)} ]\;]$
Define up to 3 subgrids.
$ n_{grids}$:
number of subgrids (integer)
$ n_{\nu}$:
the number of grid points in $ \nu$ variable (integer)
$ n_{\mu}^{(1)}\ldots n_{\mu}^{(3)}$:
number of points in the $ \mu$ direction for each of the subgrids (integers)
$ h_{\mu}^{(1)}$:
the step size in $ \mu$ variable for the first subgrid (real)
$ s^{(2)}$ $ s^{(3)}$:
step sizes (reals) for other subgrids are defined as $ h_{\mu}^{(2)}=s^{(2)} h_{\mu}^{(1)}$, $ h_{\mu}^{(3)}=s^{(3)} h_{\mu}^{(1)}$

Note that either GRID of SUBGRID label must be included in the input command file.

$ \circ$
INTERP
Format:
interp
Use this label to change the grid between separate runs of the program. The restriction is that only the number of grid points in one of the variables or $ R_{\infty}$ can be changed at a time.

 $&bull#bullet;$
INITIAL
Format:
initial $ i_1\;\;i_2\;[\;i_3\;]$
$ i_1$:
determine the initial source of orbitals and potentials:
  • $ i_1=1$ - molecular orbitals are formed as a linear combination of hydrogenic functions on centres $ A$ and $ B$; in the case of HF or HFS calculations Coulomb (exchange) potentials are approximated as a linear combination of Thomas-Fermi ($ 1/r$) potentials at the two centres; if method OED is chosen the potential function is approximated as a linear combination of $ Z_A/r_1$ and $ Z_B/r_2$ terms and the exchange potentials are set to zero

  • $ i_1=2$ - GAUSSIAN94 output is used to retrieve exponents and expansion coefficients of (uncontracted) molecular orbitals (it is assumed that the output is contained in gauss94.out and gauss94.pun files) and Coulomb and exchange potentials are initialized as in $ i_1=1$ case; see routine PREPG94 for more details

  • $ i_1=3$ - GAUSSIAN94 modified output is used to retrieve exponents and expansion coefficients of molecular orbitals (it is assumed that the output is contained in gauss94l.out file) and Coulomb and exchange potentials are initialized as in $ i_1=1$ case; see routine PREPG94L for details

  • $ i_1=5$ - initial orbitals and potentials are retrieved from disk files created in a previous run

  • $ i_1=6$ - orbitals and Coulomb potentials are retrieved from disk files and exchange potentials are initialized as in $ i_1=1$ case (convenient when going from HFS to HF calculations)

$ i_2$:
specifies how exchange potentials will be read/written and manipulated (stored in memory). The program always keeps all orbitals and Coulomb potentials in core. The exchange potentials can also be all kept in core (if there is enough memory). However, during a relaxation of a particular orbital only a fraction of them is in fact needed. Thus all exchange potentials can be kept on disk as separate files (named fort.31, fort.32, ... during a run) and only relevant ones are being retrieved when necessary.1
  • $ i_2=0$ - read exchange potentials as separate files and write them back as separate files
  • $ i_2=1$ - read all exchange potentials in a file but write them out as separate files
  • $ i_2=2$ - read all exchange potentials separately but write them out as a single file
  • $ i_2=3$ - read and write exchange potentials in the form of a single file

$ i_3$:
if $ i_1=1$ then this parameter must be set to 1 or 2. In such a case the initialization of each of the orbitals has to be defined in terms of the linear combination of atom centered hydrogen-like functions. For each orbital include a card of the following format (the order of orbitals should match the order specified under the config label):
Format:
$ c_A\;\;n_A \;\;l_A \;\;\zeta_A \;\;\;\;c_B\;\;n_B\;\;l_B\;\;
\zeta_B\;\;\;\;i_1\;\;\;\;[\;i_2\;]$
where
$ c_A$ - mixing coefficient for a hydrogenic orbital on the $ Z_A$ centre (real),
$ n_A$ - its principle quantum number (integer)
$ l_A$ - its orbital quantum number (integer)
$ \zeta_A$ - the effective nuclear charge if $ i_3=1$ or a screening parameter if $ i_3=2$ (real)
$ c_B$ - mixing coefficient for a hydrogenic orbital on the $ Z_B$ centre (real),
$ n_B$ - its principle quantum number (integer)
$ l_B$ - its orbital quantum number (integer)
$ \zeta_B$ - the effective nuclear charge if $ i_3=1$ or a screening parameter if $ i_3=2$ (real)
$ i_1$ - set to 1 to freeze the orbital during scf; otherwise set to 0 (integer)
$ i_2$ - a number of successive over-relaxations for a given orbital (integer); if omitted is set to 10

$ \circ$
FEFIELD
Format:
fefield $ r$
$ r$:
a strength of an external static electric field directed along the internuclear axis (in atomic units)

$ \circ$
MULTIPOL
Format:
multipol $ r$
$ r$:
if $ r>0$ multipole moment expansion coefficients are recalculated when the maximum error in orbital energy is reduced by $ r$ (the default value is 1.15). To suppress recalculation of the coefficients set $ r$ to a negative real number. This is useful when generating potentials from a set of fixed orbitals, e.g from GAUSSIAN94 orbitals.

$ \circ$
SCF
Format:
scf $ [\;i_1\;[\;i_2\;[\;i_3\;[\;i_4\;[i_5]\;]\;]\;]\;]$
$ i_1$:
maximum number of scf iterations (default 1000); to skip the scf step set $ i_1$ to a negative integer,
$ i_2$:
every $ i_2$ scf iterations orbitals and potentials are saved on disk (default 100). If $ i_2=0$ functions are saved on disk upon completion of the scf process. If $ i_2<0$ functions are never written to disk,
$ i_3$:
if the maximum error in orbital energy is less than $ 10^{-i_3}$ than the scf process is terminated (the default value is 10),
$ i_4$:
if the maximum error in orbital norm is less than $ 10^{-i_4}$ than scf process is terminated (the default is 10),
$ i_5$:
the level of output during scf process
  • $ i_5=1$ - the orbital convergence rate, orbital energy and normalization of every orbital is printed in every scf iteration
  • $ i_5=2$ - the orbital convergence rate, orbital energy and normalization of the worst converged orbital is printed in every scf iteration (default)
  • $ i_5=3$ - the orbital convergence rate, orbital energy and normalization of the worst converged orbital is printed every $ i_2$ iterations. Printing of ``... multipole moment expansion coefficients (re)calculated ...'' communique is suppressed
Total energy is printed every $ i_2$ iterations.

$ \circ$
FIX
Format:
fix $ [\;i_1\;[\;i_2\;[\;i_3\;]\;]\;]$
If $ i_1$, $ i_2$ or $ i_3$ are set to 1 then orbitals, Coulomb potentials or exchange potentials, respectively, are kept frozen during the scf/sor process (the respective default values are 0, 0 and 2). If $ i_3=2$ then exchange potentials are relaxed only once during an scf cycle. $ i_2$ and $ i_3$ cannot be set to 1 if hydrogenic orbitals are used to initiate the orbitals.

$ \circ$
XALPHA
Format:
xalpha $ alpha$
If this label is present the Slater exchange approximation, i.e. $ V_x=-{3\over 2}\alpha({3 \over \pi} \rho)^{1/3}$, is used and the parameter $ alpha$ can be modified. This approximation is useful to quickly improve the initial HF orbitals. $ alpha$ is a real number and 0.7 is its default value.

$ \circ$
SOR
Format:
sor $ [\;i_1\;[\;i_2\;[\;i_3\;]\;]\;]$
Change default values of sor parameters.
$ i_1$:
the number of (MC)SOR iterations (over each subgrid) for a given function being relaxed in a single SCF cycle (the default value is 10)

$ i_2$:
a scaling factor used to change $ i_1$ when relaxing potentials (see RELCOUL1/2 and RELEXCH1/2). The new value of (MC)SOR iterations is calculalated as $ i_1$/$ i_2$; the default value of $ i_2$ is 1.

$ i_3$:
if $ i_1=1 (2)$ SOR (MCSOR) method is used to solve Poisson equations for orbitals and potentials (the default value is 1); if $ i_1=3$ SOR method is used to solve Poisson equations for orbitals and MCSOR - for potentials

 $&bull#bullet;$
OMEGA
Format:
omega
$ \omega_{orb}^{(1)}\;[\;\omega_{orb}^{(2)}\;
[\;\omega_{orb}^{(3)}\;]\;]$
$ \omega_{pot}^{(1)}\;[\;\omega_{pot}^{(2)}\;
[\;\omega_{pot}^{(3)}\;]\;]$
Up to three real numbers in each line setting over-relaxation parameter $ \omega$ for relaxation of orbitals and potentials for each subgrid. If the second or the third parameter is omitted, the value for the first subgrid is used.

Note that a semiempirical formula can be used to calculate a near-optimal value of $ \omega_{pot}$ by specifying $ \omega_{pot}^{(1)}$ as a negative real number.

$ \circ$
ORDER
Format:
order $ [\;i_1\;[\;i_2\;[\;i_3\;]\;]\;]$
Up to three integers defining the ordering of mesh points on subgrids.
$ i_k=1$:
natural ordering
$ i_k=2$:
'middle' type of sweep (the default)

$ \circ$
FERMI
Format:
fermi $ r_A$ $ r_B$
When this label is present, the Fermi nuclear charge distribution is used. $ r_A$ and $ r_B$ define the atomic masses (in amu) of nuclei A and B.

$ \circ$
GAUSS
Format:
gauss $ r_A$ $ r_B$
When this label is present, the Gauss nuclear charge distribution is used. $ r_A$ and $ r_B$ define the atomic masses (in amu) of nuclei A and B.

$ \circ$
DEBUG
Format:
debug $ i_1\ldots i_{20}$
Up to 20 different debug flags can be set at a time. If the integer $ i_k$ is encountered the debug flag $ i_k$ is set, i.e. idbg $ (i_k)=1\;\; (1 \leq i_k<600)$.
 $&bull#bullet;$
STOP
Format:
stop
This label indicates the end of input data.


next up previous
Next: Examples of input command Up: users_guide Previous: users_guide
Jacek Kobus 2003-01-17