AOIntegrals - Maple Help
For the best experience, we recommend viewing online help using Google Chrome or Microsoft Edge.

Online Help

All Products    Maple    MapleSim


QuantumChemistry

  

AOIntegrals

  

compute electronic integrals in an atomic-orbital basis set

 

Calling Sequence

Description

Options

Examples

Calling Sequence

AOIntegrals(molecule, options)

Parameters

molecule

-

list of lists; each list has 4 elements, the string of an atom's symbol and atom's x, y, and z coordinates

options

-

(optional) equation(s) of the form option = value where option is one of integral_type, basis, spin, charge, symmetry, unit, max_memory

Description

• 

The AOIntegrals command generates electronic integrals in an atomic-orbital basis set.

• 

The following types of integrals, specified by the integral_type keyword, are available:
   
        - "overlap" (default): overlap matrix (S)
        - "kinetic": kinetic energy matrix (T)
        - "nuclear": nuclear attraction matrix (V)
        - "dipole": dipole moment integrals (x, y, z)
        - "eri": Electron repulsion integrals (U)
        - "quadrupole": Nuclear quadrupole moment integrals (xx, yy, zz, xy, xz, yz)
        - "giao": gauge-including atomic orbital (GIAO) integrals for NMR shielding (3 components)
        - "efg": electric-field gradient integrals (xx, yy, zz, xy, xz, yz)
.

• 

The nuclear attraction integrals, kinetic energy integrals, and overlap integrals are one-electron integrals, and hence, each is returned as an r×r Matrix where r is the number of orbitals.

• 

The electron repulsion integrals are returned as an r×r×r×r Array where r is the number of orbitals.

• 

Integrals with multiple components (e.g., dipole, quadrupole, giao, efg) are returned are 3-index Arrays with the first index being the component.

• 

On the Windows operating system the AOIntegrals command requires the installation of Microsoft's Windows Subsystem for Linux (WSL).  For Windows 10 (version 2004 and higher) and Windows 11 you can install the WSL by opening the Command Prompt in administrator mode and entering the command: wsl --install -d Ubuntu  For additional details, please refer to: https://learn.microsoft.com/en-us/windows/wsl/install

Options

• 

basis = string -- name of the basis set.  See Basis for a list of available basis sets.  Default is "sto-3g".

• 

spin = nonnegint -- twice the total spin S (= 2S). Default is 0.

• 

charge = nonnegint -- net charge of the molecule. Default is 0.

• 

symmetry = string/boolean -- is the Schoenflies symbol of the abelian point-group symmetry which can be one of the following:  D2h, C2h, C2v, D2, Cs, Ci, C2, C1. true finds the appropriate symmetry while false (default) does not use symmetry.

• 

ghost = list of lists -- each list has the string of an atom's symbol and the atom's x, y, and z coordinates.  See Ghost Atoms.

• 

unit = string -- "Angstrom" or "Bohr". Default is "Angstrom".

• 

max_memory = posint -- allowed memory in MB. Default is 4000.

See Also

QuantumChemistry

Examples

withQuantumChemistry:

Consider the hydrogen fluoride HF molecule

molecule  H,0,0,0,F,0,0,0.95;

moleculeH,0,0,0,F,0,0,0.95000000

(1)

Compute the AO overlap integrals

overlap_hf  AOIntegralsmolecule,integral_type=overlap;

overlap_hf1.000000000.045738070.427030840.0.−0.337699110.045738071.000000000.237989880.0.0.0.427030840.237989881.000000000.0.0.0.0.0.1.000000000.0.0.0.0.0.1.000000000.−0.337699110.0.0.0.1.00000000

(2)

Compute the AO kinetic integrals

kinetic_hf  AOIntegralsmolecule,integral_type=kinetic;

kinetic_hf0.76003188−0.002051800.111616900.0.−0.24905820−0.0020518036.98457674−0.212857500.0.0.0.11161690−0.212857501.037995410.0.0.0.0.0.3.248014690.0.0.0.0.0.3.248014690.−0.249058200.0.0.0.3.24801469

(3)

Compute the AO nuclear attraction integrals

nuclear_hf  AOIntegralsmolecule,integral_type=nuclear;

nuclear_hf−6.05238151−1.84242279−3.971805670.0.2.78232116−1.84242279−77.56660830−9.309375880.0.0.01402913−3.97180567−9.30937588−12.094839930.0.0.169172160.0.0.−11.966162810.0.0.0.0.0.−11.966162810.2.782321160.014029130.169172160.0.−12.07363412

(4)

Compute the AO electron repulsion integrals (ERI)

eri_hf  AOIntegralsmolecule,integral_type=eri;

Compute the AO electric dipole integrals

dipole_hf  AOIntegralsmolecule,integral_type=dipole;

 

 


Download Help Document