**INTEGRALS

General directives

.NUCMOD

Specify nuclear model.

Point nucleus:

.NUCMOD
 1

Gaussian charge distribution (default):

.NUCMOD
 2

The point nucleus model is useful to compare Lévy-Leblond type calculations with regular nonrelativistic calculations done with another code, e.g. Dalton. The two methods should give precisely the same energies.

For the Gaussian charge distribution the default exponents are in accordance with values proposed by Visscher and Dyall [Visscher1997b].

.SELECT

Restrict range of nuclei in one-electron integrals involving single atomic centers, for example electric field gradients.

Example: Restrict range to four nuclei (first number), the nuclei 1, 3, 7, and 8:

.SELECT
 4
 1
 3
 7
 8

.MAGCOR

Print the symmetrized nuclear magnetic moments. This corresponds to taking symmetry combinations of rotations, not coordinates, at each nuclear center. The numbering is used in labels of various magnetic integrals.

.PRINT

General print level in the integral module. Default:

.PRINT
 1

*ONEINT

One-electron integrals

This subsection gives directives for the generation of one-electron integrals. Based on input in the other section, the program will determine what integrals to calculate.

.PRINT

Print level in one-electron integral routines. By default print level is taken from **INTEGRALS.

*READIN

The mol file

This subsection allows changes of defaults in the reading of the mol file.

.UNCONTRACT

Decontract basis sets specified as contracted (when being read from the library, for example). In the case of using non-relativistically contracted basis set the decontraction is necessary for heavy elements. Two-component quasirelativistic Hamiltonians (like X2C) work only with decontracted basis set. By default only the small component is decontracted.

.PRINT

Print level in the reading of the mol file. By default print level is taken from **INTEGRALS.

.MAXPRI

Maximum number of primitive functions in a given block in basis file. Default:

.MAXPRI
 22

*TWOINT

Two-electron integrals

This subsection gives directives for the generation of two-electron integrals. It also gives directives for the construction of Fock matrices, such as screening.

.SCREEN

Screening threshold for integral direct calculations of Fock matrices [Saue1997]. Default:

.SCREEN
 1.0D-12

Note that the screening threshold may influence the convergence. In general, the screening threshold must be about three orders of magnitude smaller than the desired norm of the electronic gradient at convergence.

Choosing a negative value for the screening turns the integral screening completely off.

.ICEDIF

Separate screening of Coulomb and exchange contributions [Saue1997]. Useful for fine regulation of the convergence process. Default: Coulomb and exchange on (1 = on; off = 0):

.ICEDIF
 1 1

.THRFAC

Adjust the integral thresholds for SL and SS integrals. For conventional integral calculations only integrals above the threshold given in the mol file are written to disk. The thresholds for the SL and SS integrals are divided by the factors given here. Default:

.THRFAC
 1.0 1.0

.AOFOCK

Do direct Fock matrix construction in non-symmetry-adapted basis (AO basis).

The direct Fock matrix construction is performed in AO basis using the skeleton matrix approach. This may give better screening, and does give more tasks for better parallelization, but AOFOCK is more memory intensive than SOFOCK. Default is AOFOCK if 25 or more MPI nodes.

.SOFOCK

Do direct Fock matrix construction in symmetry-adapted basis (SO basis). This is the default setting.

AOFOCK may give better screening, and does give more tasks for better parallelization, but AOFOCK is more memory intensive than SOFOCK. Default is SOFOCK if at most 24 MPI nodes.

.PRINT

Set the print level in two-electron integral routines for the calculation of a particular shell quadruplet. The print level is changed only for the given shell quadruplet. A zero matches all shells, thus:

.PRINT
 4 0 0 0 0

or just:

.PRINT
 4

sets the print level to 4 for all shell quadruplets.

Use with care to avoid massive output! At print level 15 the individual integrals are printed.

.TIME

Give detailed timing for integral calculation.