This module controls the main features of the property calculation, that is,
which properties is to be calculated. In addition it includes
directives affecting the performance of several of the program
sections. It should be noted, however, that the specification of what
kind of walk (minimization , location of
transition states , dynamical
walks ) is given in the *WALK
or *MINIMI
submodules in the
general input module. See also Chapter .
READ (LUCMD,*) (CAVORG(ICOOR), ICOOR = 1, 3)
Reads the origin to be used for the cavity during a solvent calculation. By this default this is chosen to be the center of mass . Should by used with care, as it has to correspond to the center used in the evaluation of the undifferentiated solvent integrals in the HERMIT section.
READ (LUCMD, *) (DIPORG(ICOOR), ICOOR = 1, 3)
Reads in a user defined dipole origin . This may
affect properties in
which changes in dipole origin is canceled by
similar changes in the
nuclear part. It should also be used with care, as the same dipole
origin must be used during the integral evaluation sections , in
particular if one is doing numerical
differentiation with respect to electric field perturbations. For such
finite-field calculations we refer to Chapter
, which deals with finite field calculations. It's
use is mainly for debugging.
READ (LUCMD, *) (GAGORG(ICOOR), ICOOR = 1, 3)
Reads in a user defined gauge origin and overwrite both the .NOCMC option, as well as the default value of center-of-mass coordinates. Not that an unsymmetric position of the gauge origin will lead to wrong results in calculations employing symmetry.
READ (LUCMD,*) NIS
READ (LUCMD, *) (ISOTOP(IS), IS = 1, NIS)
Read one more line indicating the number of nuclei for which the isotope will be explicitly given. On the next line read the isotope chosen for each nuclei .
By default the isotope chosen for each nucleus is the most abundant
one. Although one may choose the number of nuclei for which an
isotopic substitution is to be given explicitly, one should notice
that this ordering follows the ordering of the atoms in the input.
Thus, if only the last atom is to be an isotope different from the
most common, all atoms has to be specified (and this applies to all
nuclei, not only the symmetry independent one, se also Sec. ).
The ordering of the isotopes for each nuclei is in the order of natural abundance. Thus deuterium will be isotope 2 of hydrogen, while tritium will be isotope 3.
The choice of isotopic substituted molecules will affect the gauge origin , as well as dynamical walks and walks using mass-weighted normal modes, which depend on the choice of isotopic substitution. We notice that for the vibrational analysis , the isotopic substitution may be introduced rather late, there is a similar .ISOTOP | keyword in the *VIBANA input module. This can also be used to study the vibrational and rotational structure of several isotopic substituted species.
Furthermore, the natural connection (Ref. [34, 50]) is default in order to ensure numerically stable results. The natural connection can be turned off and instead use the symmetric connection by the keyword .NODIFC .
The gauge origin are chosen to be the center of mass of the molecule. This origin can be changed by the two keywords .GAUGEO and .NOCMC . This will of course not affect the total magnetizability, only the size of the dia- and paramagnetic terms.
By definition the gauge origin of the molecular
g-factor is to be the
center of mass of the molecule, and although the
gauge origin can be
changed through the keywords .NOCMC
and .GAUGEO
, this
is not recommended, and may give erroneous results.
As the symmetric connection may give numerical inaccurate results, it's use is not recommended for other than comparisons with other programs.
READ (LUCMD, *) (ORIGIN(ICOOR), ICOOR = 1, 3)
Changes the origin of the phase-factors entering the London atomic orbitals. This will change the value of all of the contributions to the different magnetic field dependent properties when using London atomic orbitals. To be used for debugging purposes only.
READ (LUCMD, *) IPRDEF
Set default print level for the calculation. Read one
more line containing print level. Default print level is the
value of IPRDEF
from the general input module.
READ (LUMCD, *) NREPS
READ (LUMCD, *) (IDOSYM(I),I = 1, NREPS)
Consider perturbations of selected symmetries only. Read one more line specifying how many symmetries, then one line listing the desired symmetries. This option is currently only implemented for geometric perturbations.
READ (LUCMD,*) NPERT
READ (LUCMD, *) (IPOINT(I),I=1,NPERT)
Select which nuclear geometric perturbations are to be considered. Read one more line specifying how many perturbations, then on a new line the sequence of perturbations to be considered. By default, all perturbations are to be considered, but by invoking this keyword, only those perturbations specified in the sequence will be considered.
The perturbation ordering follows the ordering of the symmetrized
nuclear coordinates. This ordering can be obtained by setting the
print level in the *READIN
module to 20 or higher.
.NOLOND
.
Furthermore, the natural connection
(Ref. [34, 50]) is default in order to ensure
numerically stable results as well as physically interpretable
results for the paramagnetic and diamagnetic terms. The natural
connection can be turned off and instead use the symmetric
connection
by the keyword .NODIFC
.
The gauge origin are chosen to be the center of
mass of the molecule.
This origin can be changed by the two keywords .GAUGEO
and
.NOCMC
. This choice of gauge origin will of course not affect
the final shieldings, only the size of the dia- and paramagnetic
contributions.
.NOLOND
.
By definition the gauge origin of the
spin-rotation constant is to be the
center of mass of the molecule, and although the
gauge origin can be
changed through the keywords .NOCMC
and .GAUGEO
, this
is not recommended, and may give erroneous results.
In the current implementation, symmetry dependent nuclei cannot be used during the calculation of spin-rotation constants.
As this is a very time consuming property, as well as requiring MCSCF
wave functions in order to get reliable results, it is recommended to
consult the chapter describing the calculation of NMR-parameters
(Ch. ). The main control of which
contributions and which nuclei to calculate spin-spin couplings
between is done in the
*SPIN-S
module.
.NOLOND
.
Furthermore, the natural connection
(Ref. [34, 50]) is default in order to ensure
numerically stable results. The natural
connection can be turned off and instead use the symmetric
connection
by the keyword .NODIFC
.
The gauge origin are chosen to be the center of
mass of the molecule.
This origin can be changed by the two keywords .GAUGEO
and
.NOCMC
. This will of course not affect the final VCD results,
only the size of the contributing mechanisms.
In the current implementation, the keyword .NOCMC will be set true in calculations of Vibrational Circular Dichroism, that is, the coordinate system origin will be used as gauge origin.
.DIPGRA
the IR intensities
will be calculated as well .