next up previous contents index
Next: One-electron integrals: *ONEINT Up: **INTEGRALS directives Previous: End of input: *END   Contents   Index


General: **INTEGRALS

General-purpose directives are given in the **INTEGRALS section. This mainly includes requests for different atomic integrals, as well as some directives affecting the outcome of such an integral evaluation. Note that although not explicitly stated, none of the test options work with symmetry.

For all atomic integrals, the proper expression for the integral is given, together with the labels written on the file AOPROPER, for reference in later stages of a DALTON calculation (like for instance in during the evaluation of dynamic response properties, or for non-DALTON programs).

We also note that as long as any single atomic property integral is requested in this module, the overlap integrals will also be calculated. Note also, that unless the Hückel starting guess is turned off, this overlap matrix will not only be calculated for the requested basis set, but also for a ``ghost'' ano-4 basis set appended to the original set in order to do the Hückel starting guess.

.1ELPOT
One-electron potential energy integrals.

.AD2DAR

READ (LUCMD,*) DARFAC

Add two-electron Darwin integrals to the standard electron-repulsion integrals with a perturbation factor DARFAC.

.ANGLON
Contribution to the one-electron contribution of the magnetic moment using London orbitals arising from the differentiation of London-orbital transformed Hamiltonian, see Ref. [115].

.ANGMOM
Angular momentum around the molecular origin. This can be adjusted by changing the gauge origin through the use of the .GAUGEO keyword.

.AUXBAS
An auxiliary basis is used. Basis sets must be identified as orbital basis and auxiliary basis in the MOLECULE.INP file (line 5).

.CARMOM

READ (LUCMD,*) IORCAR

Cartesian multipole integrals to order IORCAR. Read one more line specifying order. See also the keyword .SPHMOM.

where $ii+jj+kk =$IORDER, and where ii = (i/10)*10+mod(i,10).

.CM-1

READ (LUCMD, '(A7)') FIELD1

First derivative of the electric dipole operator with respect to an external magnetic field due to differentiation of the London phase factors, see Ref. [92]. Read one more line giving the direction of the electric field (A7). These include X-FIELD, Y-FIELD, and Z-FIELD.

where $D$ is the direction of the applied electric field as specified in the input.

.CM-2

READ (LUCMD, '(A7)') FIELD2

Second derivative electric dipole operator with respect to an external magnetic field due to differentiation of the London phase factors, see Ref. [92]. Read one more line giving the direction of the electric field (A7). These include X-FIELD, Y-FIELD, and Z-FIELD.

where $D$ is the direction of the applied electric field as specified in the input.

.DARWIN
One-electron Darwin integrals [116].

.DCCR12
Only required for interfaces to other implementations of the R12 approach. Obsolete, do not use.

.DEROVL
Geometrical first derivatives of overlap integrals.

where $xyz$ is the symmetry adapted nuclear coordinate.

.DERHAM
Geometrical first derivatives of the one-electron Hamiltonian matrix.

where $xyz$ is the symmetry adapted nuclear coordinate.

.DIASUS
Diamagnetic magnetizability integrals, as calculated with London atomic orbitals, see Ref. [115]. It is calculated as the sum of the three contributions DSUSLH, DSUSLL, and DSUSNL.

.DIPGRA

Calculate dipole gradient integrals, that is, the geometrical first derivatives of the dipole length integrals.

where $abc$ is the symmetry adapted nuclear coordinate, and $d$ the direction (x/y/z) of the dipole moment.

.DIPLEN
Dipole length integrals.

.DIPORG

READ (LUCMD, *) (DIPORG(I), I = 1, 3)

Specify the dipole origin to be used in the calculation. Read one more line containing the three Cartesian components (*). Default is (0,0,0).

.DIPVEL
Dipole velocity integrals.

.DNS-KE
Kinetic-energy correction to the diamagnetic contribution to nuclear shielding constants with a common gauge origin, see Ref. [].

.DPTOVL
DPT (Direct Perturbation Theory) integrals: Small-component one-electron overlap integrals.

.DPTPOT
DPT (Direct Perturbation Theory) integrals: Small-component one-electron potential energy integrals.

.DPTPXP
DPT (Direct Perturbation Theory) integrals: Small-component dipole length integrals for Direct Perturbation Theory.

.DSO
Diamagnetic spin-orbit integrals. These are calculated using Gaussian quadrature as described in Ref. [117]. The number of quadrature point is controlled by the keyword .POINTS.

.DSO-KE
Kinetic energy correction to the diamagnetic spin-orbit integrals. These are calculated using Gaussian quadrature as described in Ref. [117]. The number of quadrature point is controlled by the keyword .POINTS. Please note that this integral has not been extensively tested, and the use of this integral is at the risk of the user.

.DSUSLH
The contribution to diamagnetic magnetizability integrals from the differentiation of the London orbital phase-factors, see Ref. [115].

.DSUSLL
The contribution to the diamagnetic magnetizability integrals from mixed differentiation on the Hamiltonian and the London orbital phase factors, see Ref. [115].

.DSUSNL
The contribution to the diamagnetic magnetizability integrals using London orbitals but with contributions from the differentiation of the Hamiltonian only, see Ref. [115].

.DSUTST
Test of the diamagnetic magnetizability integrals with London atomic orbitals. Mainly for debugging purposes.

.EFGCAR
Cartesian electric field gradient integrals.

.EFGSPH
Spherical electric field gradient integrals. Obtained by transforming the Cartesian electric-field gradient integrals (see .EFGCAR) to spherical basis.

.ELGDIA
Diamagnetic one-electron spin-orbit integrals without London orbitals.

.ELGDIL
Diamagnetic one-electron spin-orbit integrals with London orbitals.

.EXPIKR

READ (LUCMD, *) (EXPKR(I), I = 1, 3)

Cosine and sine integrals. Read one more line containing the wave numbers in the three Cartesian directions. The center of expansion is always (0,0,0).

.FC
Fermi-contact integrals, see Ref. [61].

.FC-KE
Kinetic energy correction to Fermi-contact integrals, see Ref. [].

.FINDPT

READ (LUCMD, *) DPTFAC

A direct relativistic perturbation is added to the Hamiltonian and metric with the perturbation parameter DPTFAC, where the actually applied perturbation is DPTFAC*$\alpha_{fs}^2$.

.GAUGEO

READ (LUCMD, *) (GAGORG(I), I = 1, 3)

Specify the gauge origin to be used in the calculation. Read one more line containing the three Cartesian components (*). Default is (0,0,0).

.HBDO
Symmetric combination of half-differentiated overap matrix with respect to an external magnetic field perturbation when London orbitals are used.

.HDO
Symmetrized, half-differentiated overlap integrals with respect to geometric distortions, see Ref. [118]. Differentiation on the ket-vector.

.HDOBR
Geometric half-differentiated overlap matrix differentiated once more on the ket-vector with respect to an external magnetic field, see Ref. [78].

.HDOBRT
Test the calculation of the .HDOBR integral. Mainly for debugging purposes.

.INPTES
Test the correctness of the **INTEGRALS-input. Mainly for debugging purposes, but also a good option to check if the MOLECULE input has been typed in correctly.

.KINENE
Kinetic energy integrals. Note however, that the kinetic energy integrals used in the wave function optimization is generated in the *ONEINT section.

.LONMOM
Contribution to the London magnetic moment from the differentiation with respect to magnetic field on the London orbital phase factors, see Ref. [115].

.MAGMOM
One-electron contribution to the magnetic moment around the nuclei to which the atomic orbitals are attached. This is the London atomic orbital magnetic moment as defined in Eq. (35) of Ref. [56]. The integral is calculated as the sum of .LONMOM and .ANGLON.

.MASSVE
Mass-velocity integrals.

.MGMO2T
Test of two-electron integral contribution to magnetic moment.

.MGMOMT
Test the calculation of the .MAGMOM integrals.

.MGMTHR

READ (LUCMD, *) PRTHRS

Set the threshold for which two-electron integrals should be tested with the keyword .MGMO2T. Default is 10$^{-10}$.

.MNF-SO
Calculates the atomic mean-field spin-orbit integrals as described in Ref. [103]. As the calculation of these integrals require a proper description of the atomic states, reliable results can only be expected for generally contracted basis sets such as the ANO sets, and in some cases also the correlation-consistent basis sets ((aug-)cc-p(C)VXZ)).

.NELFLD
Nuclear electric field integrals.

.NO HAM
Do not calculate ordinary one- and two-electron Hamiltonian integrals.

.NO2SO
Do not calculate two-electron contribution to spin-orbit integrals.

.NOPICH
Do not add direct perturbation theory correction to Hamiltonian integral, see keyword  .FINDPT.

.NOSUP
Do not calculate the supermatrix integral file. This may be required in order to reduce the amount of disc space used in the calculation (to approximately one-third before entering the evaluation of molecular properties). Note however, that this will increase the time used for the evaluation of the wave function significantly in ordinary Hartree-Fock runs. It is default for direct and parallel calculations.

.NOTV12
Obsolete keyword, do not use.

.NOTWO
Only calculate the one-electron part of the Hamiltonian integrals. It is default for direct and parallel calculations.

.NPOTST
Test of the nuclear potential integrals calculated with the keyword .NUCPOT. Mainly for debugging purposes.

.NSLTST
Test of the integrals calculated with the keyword .NSTLON. Mainly for debugging purposes.

.NSNLTS
Test of the integrals calculated with the keyword .NSTNOL. Mainly for debugging purposes.

.NST
Calculate the one-electron contribution to the diamagnetic nuclear shielding tensor integrals using London atomic orbitals, see Ref. [115]. It is calculated as the sum of NSTLON and NSTNOL.

.NSTCGO
Calculate the diamagnetic nuclear shielding tensor integrals without using London atomic orbitals. Note that the gauge origin is controlled by the keyword .GAUGEO.

.NSTLON
Calculate the contribution to the London orbital nuclear shielding tensor from the differentiation of the London orbital phase-factors, see Ref. [115].

.NSTNOL
Calculate the contribution to the nuclear shielding tensor when using London atomic orbitals from the differentiation of the Hamiltonian alone, see Ref. [115].

.NSTTST
Test the calculation of the one-electron diamagnetic nuclear shielding tensor using London atomic orbitals.

.NUCMOD

READ (LUCMD, *) INUC

Choose nuclear model. A 1 corresponds to a point nucleus (which is the default), and 2 corresponds to a Gaussian distribution model.

.NUCPOT
Calculate the nuclear potential energy. Currently this keyword can only be used in calculations not employing symmetry.

.OCTGRA
Calculate octupole gradient integrals, that is, the geometrical first derivatives of the third moment integrals ( .THIRDM).

where ab is the symmetry adapted nuclear coordinate, and cde the component (x/y/z) of the third moment tensor. Currently, this integral does not work with symmetry.

.OZ-KE
Calculates the kinetic energy correction to the orbital Zeeman operator, see Ref. [].

.PHASEO

READ (LUCMD, *) (ORIGIN(I), I = 1, 3)

Set the origin appearing in the London atomic orbital phase-factors. Read one more line containing the Cartesian components of this origin (*). Default is (0,0,0).

.POINTS

READ (LUCMD,*) NPQUAD

Read the number of quadrature points to be used in the evaluation of the diamagnetic spin-orbit integrals, as requested by the keyword .DSO. Read one more line containing the number of quadrature points. Default is 40.

.PRINT

READ (LUCMD,*) IPRDEF

Set default print level during the integral evaluation. Read one more line containing print level. Default is the value of IPRDEF from the general input module for DALTON.

.PROPRI
Print all one-electron property integrals requested.

.PSO
Paramagnetic spin-orbit integrals, see Ref. [61].

.PSO-KE
Kinetic energy correction to the paramagnetic spin-orbit integrals, see Ref. [].

.PSO-OZ
Orbital-Zeeman correction to the paramagnetic spin-orbit integrals, see Ref. [].

.PVP
Calculate the $pVp$ integrals that appear in the Douglas-Kroll-Heß transformation [119].

.PVIOLA
Parity-violating electroweak interaction.

.QDBINT

READ (LUCMD,'(A7)') FIELD3

London orbital corrections arising from the second-moment of charge operator in finite-perturbation calculations involving an external electric field gradient. Possible values for the perturbation (FIELD3) may be XX/XY/XZ/YY/YZ/ZZ-FGRD.

.QDBTST
Test of the .QDBINT integrals, mainly for debugging purposes.

.QUADRU
Quadrupole moment integrals. For traceless quadrupole moment integrals as defined by Buckingham [43], see the keyword .THETA.

.QUAGRA

Calculate quadrupole gradient integrals, that is, the geometrical first derivatives of the second moment integrals (i.e. .SECMOM, note: it is NOT the gradient of the .QUADRU integrals).

where abc is the symmetry adapted nuclear coordinate, and de the component (xx/xy/xz/yy/yz/zz) of the second moment tensor. Currently symmetry can not be used with these integrals.

.QUASUM
Calculate all atomic integrals as square matrices, irrespective of their inherent Hermiticity or anti-Hermiticity.

.RANGMO
Calculate the diamagnetic magnetizability integrals using the CTOCD-DZ method, see Ref. [58,59]. The gauge origin is, as default, in the center of mass.

.R12
Perform integral evaluation as required by the R12 method. One-electron integrals for the R12 method (Cartesian multipole integrals up to order 2) are precomputed and stored on the file AOPROPER. Two-electron integrals are computed in direct mode.

.R12EXP

READ (LUCMD,*) GAMMAC

Same as .R12 but with Gaussian-damped linear $r_{12}$ terms of the form $r_{12}\exp{-\gamma r_{12}^2}$. The value of $\gamma$ is read from the input line.

.R12INT
Calculation of two-electron integrals over r12.

.ROTSTR
Rotational strength integrals in the mixed representation [120].

.RPSO
Calculate the diamagnetic nuclear shielding tensor integrals using the CTOCD-DZ method, see Ref. [58,59,60]. The gauge origin is, as default, at the center of mass. Setting the gauge origin somewhere else will give wrong results in calculations using symmetry.

.S1MAG
Calculate the first derivative overlap matrix with respect to an external magnetic field by differentiation of the London phase factors, see Ref. [115].

.S1MAGL
Calculate the first magnetic half-differentiated overlap matrix with respect to an external magnetic field as needed with the natural connection, see Ref. [57]. Differentiated on the bra-vector.

.S1MAGR
Calculate the first magnetic half-differentiated overlap matrix with respect to an external magnetic field as needed with the natural connection, see Ref. [57]. Differentiated on the ket-vector.

.S1MAGT
Test the integrals calculated with the keyword .S1MAG. Mainly for debugging purposes.

.S1MLT
Test the integrals calculated with the keyword .S1MAGL. Mainly for debugging purposes.

.S1MRT
Test the integrals calculated with the keyword .S1MAGR. Mainly for debugging purposes.

.S2MAG
Calculate the second derivative of the overlap matrix with respect to an external magnetic field by differentiation of the London phase factors, see Ref. [115].

.S2MAGT
Test the integrals calculated with the keyword .S2MAG. Mainly for debugging purposes.

.SD
Spin-dipole integrals, see Ref. [61].

.SD+FC
Calculate the sum of the spin-dipole and Fermi-contact integrals.

.SD-KE
Kinetic energy correction to spin-dipole integrals, see Ref. [].

.SECMOM
Second-moment integrals.

.SELECT

READ (LUCMD, *) NPATOM
READ (LUCMD, *) (IPATOM(I), I = 1, NPATOM

Select which atoms for which a given atomic integral is to be calculated. This applies mainly to property integrals for which there exist a set of integrals for each nucleus. Read one more line containing the number of atoms selected, and then another line containing the numbers of the atoms selected. Most useful when calculating diamagnetic spin-orbit integrals, as this is a rather time-consuming calculation. The numbering is of symmetry-independent nuclei.

.SOFIEL
External magnetic-field dependence of the spin-orbit operator integrals [121].

.SOMAGM
Nuclear magnetic moment dependence of the spin-orbit operator integrals [122].

.SORT I
Requests that the two-electron integrals should be sorted for later use in SIRIUS. See also keywords .PRESORT in the **DALTON and *TRANSFORMATION input sections.

.SOTEST
Test the calculation of spin-orbit integrals as requested by the keyword .SPIN-O.

.SPHMOM

READ (LUCMD,*) IORSPH

Spherical multipole integrals to order IORSPH. Read one more line specifying order. See also the keyword .CARMOM.

where $i+j+k =$IORDER, and where ii = (i/10)*10+mod(i,10).

.SPIN-O
Spatial spin-orbit integrals, see Ref. [123]. Both the one- and the two-electron integrals are calculated, the latter is stored on the file AO2SOINT.

.SQHDOR
Square, non-symmetrized half-differentiated overlap integrals with respect to geometric distortions, see Ref. [118]. Differentiation on the ket-vector.

.SUPONL
Only calculate the supermatrix. Requires the presence of the two-electron integral file.

.SUSCGO
Diamagnetic magnetizability integrals calculated without the use of London atomic orbitals. The choice of gauge origin can be controlled by the keyword .GAUGEO.

.THETA
Traceless quadrupole moment integrals as defined by Buckingham [43].

.THIRDM
Third-moment integrals.

.U12INT
Calculation of two-electron integrals over $\left[T1,r_{12}\right]$.

.U21INT
Calculation of two-electron integrals over $\left[T2,r_{12}\right]$.

.WEINBG

READ (LUCMD,*) BGWEIN
Read in the square of the sin of the Weinberg angle appearing in the definition of parity-violating integrals, see .PVIOLA. The Weinberg angle factor will if this keyword is used be set to [1-4*BGWEIN].

.XDDXR3
Direct perturbation theory paramagnetic spin-orbit like integrals.


next up previous contents index
Next: One-electron integrals: *ONEINT Up: **INTEGRALS directives Previous: End of input: *END   Contents   Index
Dalton Manual - Release 1.2.1