nextuppreviouscontents
MOLCAS manual:

Next: 6.19 genano Up: 6. Programs Previous: 6.17 ffpt

Subsections



6.18 gateway

The Gateway module collects information about molecular system (geometry, basis sets, symmetry) to be used for future calculations.

Gateway module is a subset of seward. All keywords for this module can also appear as an input for SEWARD, however, for clearity the information about molecular system can be placed as an input for this module. Note, that gateway module does not compute any integral, and so must be followed by run of SEWARD module.

GATEWAY destroys the communication file RUNFILE, if it is used in a combination with geometry optimization it should run outside the optimization loop.

6.18.1 Input

This sections will describe the various possible input blocks in Gateway. These control
  • the molecular structure (coordinates, symmetry and basis sets),
  • explicit auxiliary basis sets in terms of CD basis sets (aCD and acCD) or external auxiliary basis sets,
  • parameters for reaction field calculations, i.e. parameters for the Kirkwood model or the PCM model and options for Pauli repulsion integral and external field integrals,
  • options for finite nuclear charge distribution models in association with relativistic calculations, and
  • the option to use the Saddle method to locate transitions state geometries.

The Gateway input section always starts with the program reference:

  &GATEWAY

6.18.1.1 General keywords

KeywordMeaning
TITLEThe keyword followed by a title.

BASDIRThe keyword allows to set up an extra location for basis set files. The value can be either an absolute path (started from /) or relative to submit directory, e.g. BASDIR=. In order to use a local copy of a basis set file with name FOO - place this file into directory specified in BASDIR

BASLIBThe keyword followed by the absolute path to the basis set library directory. The default is the $MOLCAS/basis_library directory. Note that this directory must also be host to local copies of the .tbl files.

RTRNMax number of atoms for which bond lengths, angles and dihedral angles are listed, and the radius defining the maximum length of a bond follows on the next line. The latter is used as a threshold when printing out angles and dihedral angles. The length can be followed by Bohr or Angstrom which indicates the unit in which the length was specified, the default is Bohr. The default values are 15 and 3.0 au.

6.18.1.2 Molecular structure: coordinates, symmetry and basis sets

There are three different ways to specify the molecular structure, symmetry and the basis sets in Gateway:
  • the so-called native input (old MOLCAS standard),
  • XYZ input and
Note that only XYZ input for Gateway is supported by Graphical User interface. Gateway makes a decision about the type of the input based on keywords. If Coord is used, it assumes that the input is in XYZ format.

The three different modes will be described below.

6.18.1.2.1 Native input

If the geometry is specified in a native MOLCAS format, only symmetry inequivalent atoms should be specified. The default units are atomic units. By default, symmetry is not used in the calculation.

KeywordMeaning
SYMMetrySymmetry specification follows on next line. There may be up to three different point group generators specified on that line. The generators of a point group is the minimal set of symmetry operators which is needed to generate all symmetry operators of a specific point group. A generator is in the input represented as a sequence of up to three of the characters x, y, and z. The order within a given sequence is arbitrary and the generators can be given in any sequence. Observe that the order of the irreps is defined by the order of the generators as ( E, g1, g2, g1g2, g3, g1g3, g2g3, g1g2g3)! Note that E is always assumed and should never be specified.

Below is listed the possible generators.

  • x -- Reflection in the yz-plane.
  • y -- Reflection in the xz-plane.
  • z -- Reflection in the xy-plane.
  • xy -- Twofold rotation around the z-axis.
  • xz -- Twofold rotation around the y-axis.
  • yz -- Twofold rotation around the x-axis.
  • xyz -- Inversion through the origin.
The default is no symmetry.
BASIs SetThis notes the start of a basis set definition. The next line always contains a basis set label. The basis set definition is alway terminated with the "End of Basis" keyword. For the definitions of basis set labels see the subsequent sections. Below follows a description of the options associated with the basis set definition.

  • Label [/ option] - The label is a specification of a specific basis set, e.g. C.ANO$\ldots$4s3p2d., which is an ANO basis set. If no option is specified GATEWAY will look for the basis set in the default basis directory. If an option is specified it could either be the name of an alternative basis directory or the wording ``Inline'' which defines that the basis set will follow in the current input file. For the format of the Inline option see the section `Basis set format'. Observe that the label is arbitrary for this option and will not be decoded. The Label card is mandatory.
  • Name x, y, z (Angstrom or Bohr) - This card specifies an arbitrary (see next sentence!) name for a symmetry distinct center and its Cartesian coordinates. Observe, that the name "DBAS" is restricted to assign the center of the diffuse basis functions required to model the continuum orbitals in R-matrix calculations. The label is truncated to four characters. Observe that this label must be unique to each center. The coordinate unit can be specified as an option. The default unit is Bohr. There should at least be one card of this type in a basis set definition.
  • Charge - The real entry on the subsequent line defines the charge associated with this basis set. This will override the default which is defined in the basis set library. The option can be used to put in ghost orbitals as well as to augment the basis sets of the library. The Charge card is optional.
  • Spherical [option] - Specifying which shells will be in real spherical Gaussians. Valid options are "all" or a list of the shell characters separated by a blank. The shell characters are s, p, d, f, etc. All shells after p are by default in real spherical Gaussians, except for the d-functions in the 6-31G family of basis sets which are in Cartesian. The Spherical card is optional. The s and p shells and the d-functions of the 6-31G family of basis sets are by default in Cartesian Gaussians.
  • Cartesian [option] - Specifying which shells will be in a Cartesian Gaussian representation. For syntax consult the corresponding Spherical keyword.
  • Contaminant [option] - Specifying for which shells the contaminant will be kept. The contaminants are functions of lower rank which are generated when a Cartesian shell is transformed to a spherical representation (e.g. r2=x2+y2+z2 for d-shells, p contaminants for f-shells, s and d contaminants for g-shells, etc). Valid options are the same as for the Spherical keyword. The default is no contaminant in any shell. The Contaminant card is optional.
  • End of Basis set - Marks the end of the basis set specification. This card is mandatory.

Example of an input in native MOLCAS format:


  &GATEWAY
Title
formaldehyde
SYMMETRY
X  Y
Basis  set
H.STO-3G....
H1  0.000000  0.924258  -1.100293  /Angstrom
End  of  basis

Basis  set
C.STO-3G....
C3  0.000000  0.000000  -0.519589  /Angstrom
End  of  basis

Basis  set
O.STO-3G....
O  0.000000  0.000000  0.664765  /Angstrom
End  of  basis

End  of  input

6.18.1.2.2 XYZ input

If the geometry is specified in XYZ format, all atoms should be specified. The default units are Ångstroms. By default, maximum possible symmetry is used.

'Molcas XYZ' file format is an extension of plain XYZ format.

First line of this file contains the number of atoms.
Second line (a comment line) can contain 'a.u.' or 'bohr' to use atomic units, instead of default Ångstroms. Also this line can contain keyword TRANS, followed by 3 numbers, and/or ROTATE, followed by 9 numbers (in this case coordinates will be Translated by specified vector, and/or Rotated).
Remaining lines are used to specify Element and cartesian coordinates.

Element name might be optionally followed by a Number (e.g. H7), a Label (separated by _ sign: e.g. $H\_INNER$), or Basis Set (separated by . , e.g. H.STO-3G)

KeywordMeaning
COORD

The keyword followed on the next line by the name of XYZ file, or inline coordinates in XYZ format. If the file is located in the same directory, where MOLCAS job was submitted there is no need to specify the PATH to this file. The keyword may appear several times. In this case all coordinate files will be concatenated, and considered as individual fragments.

BASISThe keyword can be used to specify global basis set for all atoms, or for a group of atoms. The keyword followed by a label of basis set, or by coma separated list of basis sets for individual atoms.

Note! The basis set definition in XYZ mode does not allow to use inline basis set.

Example:

COORD
4

C 0.00000 0.00000 0.00000
H 1.00000 0.00000 0.00000
H 0.00000 1.00000 0.00000
H 0.00000 0.00000 1.00000
BASIS
STO-3G, H.6-31G*

In this example, the C atom (in the origin) will have the basis set STO-3G and the H atoms 6-31G*.

If keyword BASIS never appears in the input, the default basis, ANO-S-MB, will be used.

GROUP

The keyword can be used to specify the symmetry of the molecule.

The keyword must be followed by one of:

  • FULL (default) - use maximum possible subgroup of D2h
  • NOSYM (same as E, or C1)
  • space separated list of generators: e.g. X XY (for more details see SYMMETRY keyword)

Limitations: in the current implementation atom labels, and basis sets are ignored during symmetry recognition.

If XYZ input has been used in gateway, a file with native MOLCAS input will be produced and stored in working directory under the name findsym.std.

Advanced keywords:
KeywordMeaning
SYMThresholdfollowed by a real number - threshold for symmetry recognition (default is 0.1)

MOVEallow to translate and rotate molecule in order to find highest possible symmetry. (this is a default for all groups, except of C1)
NOMOVEdo not allow to transform coordinates while searching for highest group (default for C1 group)

BSSEfollowed by an integer. Indicates which XYZ-file that should be treated like ghost atoms.

VARTSpecifies that the energy should not be considered invariant to translations. Translational variance is detected automatically, but sometimes it may be useful to enforce it.

VARRSpecifies that the energy should not be considered invariant to rotations. Rotational variance is detected automatically, but sometimes it may be useful to enforce it.



Example:
&GATEWAY
COORD
water.xyz
BASIS
STO-3G

or, in short EMIL notation:

&GATEWAY
COORD=water.xyz; BASIS=STO-3G

6.18.1.3 Constraints

In case of optimizations with constraints these are defined in the GATEWAY input. For a complete description of this keyword see the section [*].

KeywordMeaning
CONStraintsThis marks the start of the definition of the constraints which the optimization is subject to. This section is always ended by the keyword End of Constraints. This option can be used in conjunction with any definition of the internal coordinates.

6.18.1.4 Explicit auxiliary basis sets

The so-called Resolution of Identity (RI) technique (also called Density Fitting, DF) is implemented in the MOLCAS package. This option involves the use of an auxiliary basis set in the effective computation of the 2-electron integrals. MOLCAS incorporates both the use of conventionally computed, externally provided, auxiliary basis sets (RIJ, RIJK, and RIC types), and on-the-fly generated auxiliary basis sets. The latter are atomic CD (aCD) or the atomic compact CD (acCD) basis sets, based on the Cholesky decomposition method. The externally provided auxiliary basis sets are very compact, since they are tailored for special wave function methods. However, they are not provided for all available valence basis sets. The aCD or acCD RI auxiliary basis sets are a more general option and provides auxiliary basis sets for any wave function model and valence basis set.

KeywordMeaning
RIJUse the RI-J basis in the density fitting (DF) approach to treat the two-electron integrals. Note that the valence basis set must have a supporting auxiliary basis set for this to work.
RIJKUse the RI-JK auxiliary basis in the density fitting (DF) approach to treat the two-electron integrals. Note that the valence basis set must have a supporting auxiliary basis set for this to work.
RICUse the RI-C auxiliary basis in the density fitting (DF) approach to treat the two-electron integrals. Note that the valence basis set must have a supporting auxiliary basis set for this to work.
RICDUse the aCD or acCD approach [7] to treat the two-electron integrals. This procedure will use an on-the-fly generated auxiliary basis set.
CDTHresholdThreshold for on-the-fly generation of aCD or acCD auxiliary basis sets for RI calculations (default value 1.0d-4).
SHACSkip high angular combinations à la Turbomole when creating on-the-fly basis sets (default of).
KHACKeep high angular combinations when creating on-the-fly basis sets (default on).
aCD basisGenerate an atomic CD (aCD) auxiliary basis sets (default off).
acCD basisGenerate an atomic compact CD (acCD) auxiliary basis sets (default on).


6.18.1.5 Reaction field calculations

The effect of the solvent on the quantum chemical calculations has been introduced in MOLCAS through the reaction field created by the surrounding environment, represented by a polarizable dielectric continuum outside the boundaries of a cavity containing the solute molecule. MOLCAS-4 supports Self Consistent Reaction Field (SCRF) and Multi Configurational Self Consistent Reaction Field (MCSCRF) calculations within the framework of the SCF and the RASSCF programs. The reaction field, computed in a self-consistent fashion, can be later added as a constant perturbation for the remaining programs, as for example CASPT2.

The purpose of this facility is to incorporate the effect of the environment (a solvent or a solid matrix) on the studied molecule. The utility itself it is not a program, but requires an additional input which has to be provided to the GATEWAY program. Two methods are available for SCRF calculations: one is based on the Kirkwood model, the other is the so called Polarizable Continuum Model (PCM). The reaction field is computed as the response of a dielectric medium polarized by the solute molecule: the solute is placed in a ``cavity'' surrounded by the dielectric. In Kirkwood model the cavity is always spherical, whereas in PCM the cavity is modeled on the actual solute shape.

The possible set of parameters controlled by input are:

  • the Kirkwood model,
  • the PCM model, and
  • one-electron integrals representing Pauli repulsion and external fields.
First a brief presentation of the Kirkwood and the PCM models.


6.18.1.5.1 The Kirkwood Model

The Kirkwood model is an expansion of the so-called Onsager model where the surrounding will be characterized by its dielectric permitivity and a radius describing a spherical cavity, indicating where the dielectric medium starts. (Note that all atoms in the studied molecule must be inside the spherical cavity.) The Pauli repulsion due to the medium can be introduced by use of the spherical well integrals which are generated by SEWARD. The charge distribution of the molecule will introduce an electric field acting on the dielectric medium. This reaction field will interact with the charge distribution of the molecule. This interaction will manifest itself as a perturbation to the one-electron Hamiltonian. The perturbation will be automatically computed in a direct fashion (no multipole integrals are stored on disk) and added to the one-electron Hamiltonian. Due to the direct way in which this contribution is computed rather high terms in the multipole expansion of the charge can be afforded.


6.18.1.5.2 The Polarizable Continuum Model, PCM

The PCM has been developed in order to describe the solvent reaction field in a more realistic way, basically through the use of cavities of general shape, modeled on the solute. The cavity is built as the envelope of spheres centered on solute atoms or atomic groups (usually, hydrogen atoms are included in the same sphere as the heavy atoms they are bonded to). The reaction field is described by means of apparent charges (solvation charges) spread on the cavity surface, designed to reproduced the electrostatic potential due to the polarized dielectric inside the cavity. Such charges are used both to compute solute-solvent interactions (modifying the total energy of the solute), and to perturb the molecular Hamiltonian through a suitable operator (thus distorcing the solute wave-function, and affecting all the electronic properties). The PCM operator contains both one- and two-electron terms: it is computed using atomic integrals already present in the program, through a ``geometry matrix'' connecting different points lying on the cavity surface. It can be shown that with this approach the SCF and RASSCF variational procedures lead to the free energy of the given molecule in solution: this is the thermodynamic meaning of the SCF or CI energy provided by the program. More precisely, this is the solute-solvent electrostatic contribution to the free energy (of course, other terms depending on solute atomic motions, like vibrational and rotational free energies, should be included separately); it can be used to get a good approximation of the solvation free energy, by subtracting the SCF or CI energy computed in vacuo, and also to compute directly energy surfaces and reaction paths in solution. On the other hand, the solute wave-function perturbed by the reaction field can be used to compute any electronic property in solution.

Also other quantities can be computed, namely the cavitation free energy (due the the work spent to create the cavity in the dielectric) and the dispersion-repulsion free energy: these terms affect only the total free energy of the molecule, and not its electronic distribution. They are collectively referred to as non-electrostatic contributions.

Note that two other keywords are defined for the RASSCF program: they refer to the CI root selected for the calculation of the reaction field (RFROOT), and to the possibility to perform a non-equilibrium calculation (NONEQ) when vertical electronic transitions are studied in solution. These keywords are referenced in the RASSCF section. To include the reaction field perturbation in a SCF, RASSCF, CASPT2 or RASSI calculation, another keyword must be specified (RFPERT), as explained in the respective program sections.

Complete and detailed examples of how to add a reaction field, through the Kirkwood or the PCM model, into quantum chemical calculations in MOLCAS is presented in section [*] of the examples manual. The user is encouraged to read that section for further details.

6.18.1.5.3 Input for the Kirkwood and PCM models


6.18.1.5.3.1 Files

The reaction field calculations will store the information in the following files, which will be used by the following programs

FileContents
ONEINTOne-electron integral file used to store the Pauli repulsion integrals
RUNFILECommunications file. The last computed self-consistent reaction field (SCF or RASSCF) will be stored here to be used by following programs
GV.offInput file for the external program ``geomview'' (see Tutorial section ``Solvent models''), for the visualization of PCM cavities

6.18.1.5.3.2 Input
Below follows a description of the input to the reaction field utility in the GATEWAY program. The RASSCF program has its own keywords to compute reaction fields for excited states.

Compulsory keywords
KeywordMeaning
RF-InputActivate reaction field options.

END Of RF-InputThis marks the end of the input to the reaction field utility.

Optional keywords for the Kirkwood Model
KeywordMeaning
REACtion FieldThis command is exclusive to the Kirkwood model. It indicates the beginning of the specification of the reaction field parameters. The subsequent line will contain the dielectric constant of the medium, the radius of the cavity in Bohrs (the cavity is always centered around the origin), and the angular quantum number of the highest multipole moment used in the expansion of the change distribution of the molecule (only charge is specified as 0, charge and dipole moments as 1, etc.). The input specified below specifies that a dielectric permitivity of 80.0 is used, that the cavity radius is 14.00 a.u., and that the expansion of the charge distribution is truncated after l=4, i.e hexadecapole moments are the last moments included in the expansion. Optionally a fourth argument can be added giving the value of the dielectric constant of the fast component of the solvent (default value 1.0).
Sample input for the reaction field part (Kirkwood model)



RF-Input
Reaction  field
80.0  14.00  4
End  Of  RF-Input

Sample input for a complete reaction field calculation using the Kirkwood model. The SCF computes the reaction field in a self consistent manner while the MRCI program adds the effect as a constant perturbation.



  &GATEWAY
Title  =  HF  molecule
Symmetry
X  Y
Basis  set
F.ANO-S...3S2P.
F  0.00000  0.00000  1.73300
End  of  basis
Basis  set
H.ANO-S...2S.
H  0.00000  0.00000  0.00000
End  of  basis
Well  integrals
  4
  1.0  5.0  6.75
  1.0  3.5  7.75
  1.0  2.0  9.75
  1.0  1.4  11.75

RF-Input
Reaction  field
  80.0  4.75  4
End  of  RF-Input
  &SEWARD

  &SCF
Occupied  =  3  1  1  0

  &MOTRA
LumOrb
Frozen  =  1  0  0  0
RFPert

  &GUGA
Electrons  =  8
Spin  =  1
Inactive  =  2  1  1  0
Active  =  0  0  0  0
CiAll  =  1

  &MRCI
SDCI

Optional keywords for the PCM Model
KeywordMeaning
PCM-modelIf no other keywords are specified, the program will execute a standard PCM calculation with water as solvent. The solvent reaction field will be included in all the programs (SCF, RASSCF, CASPT2, etc) invoked after SEWARD: note that in some cases additional keywords are required in the corresponding program sections. Some PCM parameters can be changed through the following keywords.
SOLVentUsed to indicate which solvent is to be simulated. The name of the requested solvent must be written in the line below this keyword. Find implemented solvents in the PCM model below this section.
DIELectric constantDefines a different dielectric constant for the selected solvent; useful to describe the system at temperatures other that 298 K, or to mimic solvent mixtures. The value is read in the line below the keyword. An optional second value might be added on the same line which defines a different value for the infinite frequency dielectric constant for the selected solvent (this is used in non-equilibrium calculations; by default it is defined for each solvent at 298 K).
CONDuctor versionIt requires a PCM calculation where the solvent is represented as a polarized conductor: this is an approximation to the dielectric model which works very well for polar solvents (i. e. dielectric constant greater than about 5), and it has some computational advantages being based on simpler equations. It can be useful in cases when the dielectric model shows some convergence problems.
AAREaIt is used to define the average area (in Å2) of the small elements on the cavity surface where solvation charges are placed; when larger elements are chosen, less charges are defined, what speeds up the calculation but risks to worsen the results. The default value is 0.4 Å2 (i. e. 60 charges on a sphere of radius 2 Å). The value is read in the line below the keyword.
R-MInIt sets the minimum radius (in Å) of the spheres that the program adds to the atomic spheres in order to smooth the cavity surface (default 0.2 Å). For large solute, if the programs complains that too many sphere are being created, or if computational times become too high, it can be useful to enlarge this value (for example to 1 or 1.5 Å), thus reducing the number of added spheres. The value is read in the line below the keyword.
PAULingIt invokes the use of Pauling's radii to build the solute cavity: in this case, hydrogens get their own sphere (radius 1.2 Å).
SPHEre radiusIt is used to provide sphere radii from input: for each sphere given explicitly by the user, the keyword ``Sphere radius'' is required, followed by a line containing two numbers: an integer indicating the atom where the sphere has to be centered, and a real indicating its radius (in Å). For example, ``Sphere radius'' followed by ``3 1.5'' indicates that a sphere of radius 1.5 Å is placed around atom #3; ``Sphere radius'' followed by ``4 2.0'' indicates that another sphere of radius 2 Å is placed around atom #4 and so on.

Solvents implemented in the PCM model are

Name  Dielectric Name  Dielectric Name  Dielectric
   constant    constant    constant
water  78.39 dichloroethane  10.36 toluene  2.38
dimethylsulfoxide  46.70 quinoline  9.03 benzene  2.25
nitromethane  38.20 methylenchloride  8.93 carbontetrachloride  2.23
acetonitrile  36.64 tetrahydrofuran  7.58 cyclohexane  2.02
methanol  32.63 aniline  6.89 heptane  1.92
ethanol  24.55 chlorobenzene  5.62 xenon  1.71
acetone  20.70 chloroform  4.90 krypton  1.52
isoquinoline  10.43 ethylether  4.34 argon  1.43

Sample input for the reaction field part (PCM model): the solvent is water, a surface element average area of 0.2 Å2 is requested.



RF-input
PCM-model
Solvent
water
AAre
0.2
End  of  RF-input

Sample input for a standard PCM calculation in water. The SCF and RASSCF programs compute the reaction field self consistently and add its contribution to the Hamiltonian. The RASSCF is repeated twice: first the ground state is determined, then a non-equilibrium calculation on the first excited state is performed.



  &GATEWAY
Coord
4
formaldehyde
O  0.000000  0.000000  -1.241209
C  0.000000  0.000000  0.000000
H  0.000000  0.949585  0.584974
H  0.000000  -0.949585  0.584974

Basis  =  STO-3G
Group  =  C1
RF-input
PCM-model
solvent  =  water
End  of  RF-input

  &SEWARD  ;  &SCF

  &RASSCF
nActEl  =  4  0  0
Symmetry  =  1
Inactive  =  6
Ras2  =  3
CiRoot
1  1
1
LumOrb

  &RASSCF
nActEl  =  4  0  0
Symmetry  =  1
Inactive  =  6
Ras2  =  3
CiRoot
2  2
1  2
0  1
JOBIPH
NonEq
RFRoot  =  2
Again the user is recommended to read section [*] of the examples manual for further details.

6.18.1.5.4 Keywords associated to one-electron integrals

KeywordMeaning
WELL integralsRequest computation of Pauli repulsion integrals for dielectric cavity reaction field calculations. The first line specifies the total number of primitive well integrals in the repulsion integral. Then follows a number of lines, one for each well integral, specifying the coefficient of the well integral in the linear combination of the well integrals which defines the repulsion integral, the exponent of the well integral, and the distance of the center of the Gaussian from the origin. In total three entries on each line. All entries in atomic units. If zero or a negative number is specified for the number of well integrals a standard set of 3 integrals with their position adjusted for the radius of the cavity will be used. If the distance of the center of the Gaussian from the origin is negative displacements relative to the cavity radius is assumed.
XFIEld integralsRequest the presence of an external electric field represented by a number of partial charges and dipoles. Optionally, polarisabilities may be specified whose induced dipoles are determined self-consistently during the SCF iteration. The first line may contain, apart from the first integer [nXF] (number of centers), up to four additional integers. The second integer [nOrd] specifies the maximum multipole order, or -1 signifying no permanent multipoles. Default is 1 (charges and dipoles). The third integer [p] specifies the type of external polarisabilities: 0 (default) no polarisabilities, 1 (isotropic), or 2 (anisotropic). The fourth integer [nFrag] specifies the number of fragments one multipole may contribute to (relevant only if polarisabilities are present). The default is 0, meaning that each permanent multipole is only excluded in the calculation of the field at its own polarisability, 1 means that one gives a fragment number to each multipole and that the static multipoles do not contribute to the polarising field within the same fragment, whereas 2 can be used in more complex situations, e.g. polymers, allowing you to specify a second fragment number so that junction atoms does not contribute to either of the neighbouring fragments. Finally, the fifth and last integer [nRead] (relevant only if Langevin dipoles are used) may be 0 or 1 (where 0 is default), specifying wheather an element number (e.g. 8 for oxygen) should be read for each multipole. In that case the default radius for that element is used to determine which Langevin grid points should be annihilated. A negative element number signifies that a particular radius should be used for that multipole, in thousands of a Bohr (-1400 meaning 1.4 Bohr). Then follows nXF lines, one for each center. On each line is first nFrag+nRead (which may equal 0) integers, specifying the fragments that the multipole should not contribute to (the first fragment is taken as the fragment that the polarisability belongs to) and the element number. Then follows the three coordinates of the center, followed by the multipoles and polarisabilities. The number of multipole entries is 0 for nOrd=-1, 1 for nOrd=0, 4 for nOrd=1, and 10 for nOrd=2. The number of polarisability entries are 0 for p=0, 1 for p=1, and 6 for p=2. The order of quadrupole moment and anisotropic polarisability entries is xx, xy, xz, yy, yz, zz. If default is used, i.e. only specifying the number of centers on the first line, each of these lines will contain 7 entries (coordinates, charge, and dipole vector). All entries are in atomic units, if not otherwise requested by the Angstrom keyword that must be placed between nXF and nOrd. All these data can be stored in a separate file whose name must be passed as an argument of the XField keyword.
ANGMSupplement ONEINT for transition angular momentum calculations. Entry which specifies the angular momentum origin (in au).
AMPRRequest the computation of angular momentum product integrals. The keyword is followed by values which specifies the angular momentum origin (in au).
DSHDRequests the computation of diamagnetic shielding integrals. The first entry specifies the gauge origin. Then follows an integer specifying the number of points at which the diamagnetic shielding will be computed. If this entry is zero, the diamagnetic shielding will be computed at each nucleus. If nonzero, then the coordinates (in au) for each origin has to be supplied, one entry for each origin.

EPOTAn integer follows which represents the number of points for which the electric potential will be computed. If this number is zero, the electric field acting on each nucleus will be computed. If nonzero, then the coordinates (in au) for each point have to be supplied, one entry for each point. This keyword is mutually exclusive with EFLD and FLDG.
EFLDAn integer follows which represents the number of points for which the electric potential and electric field will be computed. If this number is zero, the electric field acting on each nucleus will be computed. If nonzero, then the coordinates (in au) for each point have to be supplied, one entry for each point. This keyword is mutually exclusive with EPOT and FLDG.
FLDGAn integer required which represents the number of points for which the electric potential, electric field and electric field gradient will be computed. If this number is zero, the electric field gradient acting on each nucleus will be computed. If nonzero, then either the coordinates (in au) for each point or labels for each atom center have to be supplied, one entry for each point. In case a label is supplied it must match one of those given previous in the input during specification of the coordinates of the atom centers. Using a label instead of a coordinate can e.g. be useful in something like a geometry optimization where the coordinate isn't known when the input is written. This keyword is mutually exclusive with EPOT and EFLD.
EMPCUse point charges specified by the keyword XField when calculating the Orbital-Free Embedding potential.
RF-InputSpecification of reaction field parameters, consult the reaction field section of this manual.

6.18.1.5.5 Keywords associated with nuclear charge distribution models

Input parameters associated with different models of the nuclear charge distribution. The default is to use a point charge representation.
KeywordMeaning
FINIteRequest a finite center representation of the nuclei by a single exponent s-type Gaussian.
MGAUSsianRequest a finite center representation of the nuclei by a modified Gaussian.

6.18.1.6 The Saddle method for transition state optimization

The Saddle method [45] is a method to locate transition states (TS). The method, in practice, can be viewed as a series of constrained optimization along the reaction path, which connects two starting structure (could be the reactants and products of a reaction), to locate the region of the TS and a subsequent unconstrained optimization to locate the TS. The only data needed for the procedure are the energies and coordinates of the two structures. Note that this option will overwrite the coordinates which have already been specified with the normal input of the molecular geometry. However, this does not make that input section redundant and should always be included.

KeywordMeaning
RP-CoordinatesThis activates the Saddle method for TS geometry optimization. The line is followed by an integer specifying the number of symmetry unique coordinates to be specified. This is followed by two sets of input - one line with the energy and then the Cartesian coordinates in bohr - for each of the two starting structures of the Saddle method. Note that the order of the coordinates must always match the order specified with the conventional input of the coordinates of the molecular system. Alternatively, two lines with the filenames containing the coordinates of reactants and products, respectively, (in XYZ format) can be given.

NOALignBy default, the two starting structures are aligned to minimize the root mean square distance (RMSD) between them, in particular, the first structure is moved and the second structure remains fixed. If this keyword is given, the starting structures are used as given.

ALIGn onlyThe two starting structures are aligned, but nothing more is done. An input block for seward is still needed, but no integrals are computed.

WEIGhtsRelative weights of each atom to use for the alignment and for the calculations of the ``distance'' between structures. The possibilities are:

MASS. This is the default. Each atom is given a weight proportional to its mass. Equivalent to mass-weighted coordinates.

EQUAL. All atoms have an equal weight.

HEAVY. Only heavy atoms are considered, with equal weights. Hydrogens are given zero weight.

A list of N numbers can also be provided, and they will be used as weights for the N symmetry-unique atoms. For example:



WEIGht
0 0 0 0 0 0 1 1 1 1 1 1 0 0 0 0

will align only atoms 7-12 out of 16.

Note that, in any case, weights of 0 are likely to cause problems with constraints, and they will be increased automatically.

SADDleStep size reduction for each macro iteration of the saddle method. The value is given in weighted coordinates, divided by the square root of the total weight (see the WEIGHTS keyword). Default value is 0.1 au.

6.18.1.7 Geometry Optimization using constrained internal coordinates

These keyword are used together with the geo to optimize the relative position of two or more rigid fragments. The starting geometry can either be defined by supplying an xyz-file for each fragment using the keyword coord or by placing a file named $project.zmt in a directory named $project.GEO . The z-matrix should be in the following format:



H
O 0.982011 0 1
H 0.982013 0 104.959565 0 2 1
H 1.933697 1 107.655494 1 114.496053 1 2 3 1
O 0.988177 0 173.057942 1 -56.200750 1 4 2 3
H 0.979890 0 104.714572 0 179.879745 1 5 4 2

where the three columns of real numbers are internal coordinates, and the last three columns of integers indicate which other atoms that are used to define the coordinate. The type of coordinates from left to right are bond distances, bond angles and dihedral angels, both for the coordinates and the link. The column of integers just to the right of each coordinate indicate if this coordinate should be optimized or not (1 = optimize, 0 = do not optimize).

There are also two utility-keywords used to create a z-matrix or to write out a constraint-definition for slapaf and keywords to rotate and translate fragments. (See documentation for GEO for more details)

KeywordMeaning
HYPERThis keyword is used to specify that a geometry optimization with constrained internal coordinates shall be performed later, a z-matrix and a set of displaced geometries are therefore constructed. The keyword should be followed by three real numbers defining the maximum displacement for each coordinate type. The order from left to right is bond distances, bond angles, and dihedral angles. Default values of 0.15 Å, 2.5 degrees, and 2.5 degrees are used for the maximum displacement of bond distances, bond angles and dihedral angles respectively.

OLDZThis keyword is used both to start a new calculation from a user-defined z-matrix and to restart calculations. When using the keyword for a new calculation a directory $project.GEO must exist and contain a file called $project.zmt with a z-matrix in the format defined above. The directory must not contain any files with the suffix .info when performing a fresh calculation since these files contain restart information.

ZONLYThis keyword is used to construct a z-matrix from a set of xyz-files (fragments) and store it in the directory $project.GEO. The optimization parameters of the resulting z-matrix are set so that only coordinates linking fragments are set to 1 (= optimize coordinate).

ZCONSThis keyword is used to define constraints from a set of xyz-files (fragments) on a form that could be supplied to the slapaf in order to keep the fragments rigid. The resulting constraints-file is named $project.cns and stored in the directory $project.GEO. The atom-numbers in this constraint-file will not match those of your original xyz-file and should not be used together with this. Instead a new xyz-file named cons.xyz is created and placed into the directory $project.GEO, this has the proper numbering to use together with the constraints.

ORIGINThis keyword is used to translate and rotate a set of xyz-files. The keyword must be entered before the xyz-files is entered with coord. The keyword should be followed by one line with the number of fragments and then one line for each fragment that should be translated. This row should contain 13 numbers. One integer defining which fragment that should be moved, (the fragments are numbered based on order of appearance in the input-file from top to bottom), 3 real numbers defining a translation (x, y, z) and 9 numbers defining a rotation (row1, row2, row3 of rotation matrix). The keyword origin is mutually exclusive with the keyword frgm which is an alternative way to express the same rotations and translations.

FRGMThis keyword is used together with the keywords rot and trans to define rotation and translation of a specific fragment. Frgm defines an active fragment (each xyz-file is considered a fragment, the files are numbered based on order of appearance in the input from top to bottom). The keyword must be entered before the xyz-file it is supposed to modify is entered with coord. Each occurence of frgm should be followed by either one of or both of the keywords rot and trans to define rotation and translation. This keyword is mutually exclusive with the keyword orgin

ROTThis keyword should be followed by nine real numbers defining the rotation for the fragment defined by the preceeding frgm. The numbers should be the nine elements of a rotation matrix listed with one full row at the time.

TRANSThis keyword should be followed by three real numbers defining the translation for the fragment defined by the preceeding frgm. The numbers should be the x, y and z coordinates of the translation in that order.

Example of an input:

  &GATEWAY
Title
Water  Dimer
frgm=2
trans=3.0  0.0  0.0
Coord=water_monomer.xyz
Coord=water_monomer.xyz
Group=c1
basis=cc-pVTZ
hyper
0.2  3.0  3.0

In this example a water dimer is constructed from a single monomer by translating it with the keyword trans. For more details see the manual entry for the module geo.


next up previous contents
Next: 6.19 genano Up: 6. Programs Previous: 6.17 ffpt