#include <gaussbas.h>
Inheritance diagram for sc::GaussianBasisSet:
Public Types | |
enum | UnitType { Unit } |
This can be given to a CTOR to construct a unit basis function. | |
Public Member Functions | |
GaussianBasisSet (const Ref< KeyVal > &) | |
The KeyVal constructor. | |
GaussianBasisSet (UnitType) | |
This can be given GaussianBasisSet::Unit to construct a basis set with a single basis function that is one everywhere. | |
GaussianBasisSet (StateIn &) | |
Ref< GaussianBasisSet > | operator+ (const Ref< GaussianBasisSet > &B) |
Returns a GaussianBasisSet object that consists of the basis functions for each atom in this followed by the basis functions in B for the corresponding atom. | |
void | save_data_state (StateOut &) |
Save the base classes (with save_data_state) and the members in the same order that the StateIn CTOR initializes them. | |
const char * | name () const |
Return the name of the basis set (is nonnull only if keyword "name" was provided). | |
const char * | label () const |
Return the label of the basis set. | |
Ref< Molecule > | molecule () const |
Return the Molecule object. | |
Ref< SCMatrixKit > | matrixkit () |
Returns the SCMatrixKit that is to be used for AO bases. | |
Ref< SCMatrixKit > | so_matrixkit () |
Returns the SCMatrixKit that is to be used for SO bases. | |
RefSCDimension | basisdim () |
Returns the SCDimension object for the dimension. | |
int | ncenter () const |
Return the number of centers. | |
int | nshell () const |
Return the number of shells. | |
int | nshell_on_center (int icenter) const |
Return the number of shells on the given center. | |
int | shell_on_center (int icenter, int shell) const |
Return an overall shell number, given a center and the shell number on that center. | |
int | shell_to_center (int ishell) const |
Return the center on which the given shell is located. | |
int | shell_to_primitive (int ishell) const |
Return the overall index of the first primitive from the given shell. | |
int | nbasis () const |
Return the number of basis functions. | |
int | nbasis_on_center (int icenter) const |
Return the number of basis functions on the given center. | |
int | nprimitive () const |
Return the number of primitive Gaussians. | |
int | has_pure () const |
Return true if basis contains solid harmonics Gaussians. | |
int | max_nfunction_in_shell () const |
Return the maximum number of functions that any shell has. | |
int | max_ncartesian_in_shell (int aminc=0) const |
Return the maximum number of Cartesian functions that any shell has. | |
int | max_nprimitive_in_shell () const |
Return the maximum number of primitive Gaussian that any shell has. | |
int | max_angular_momentum () const |
Return the highest angular momentum in any shell. | |
int | max_ncontraction () const |
Return the maximum number of Gaussians in a contraction in any shell. | |
int | max_am_for_contraction (int con) const |
Return the maximum angular momentum found in the given contraction number for any shell. | |
int | max_cartesian () const |
Return the maximum number of Cartesian functions in any shell. | |
int | shell_to_function (int i) const |
Return the number of the first function in the given shell. | |
int | function_to_shell (int i) const |
Return the shell to which the given function belongs. | |
const GaussianShell & | operator() (int i) const |
Return a reference to GaussianShell number i. | |
GaussianShell & | operator() (int i) |
Return a reference to GaussianShell number i. | |
const GaussianShell & | operator[] (int i) const |
Return a reference to GaussianShell number i. | |
GaussianShell & | operator[] (int i) |
Return a reference to GaussianShell number i. | |
const GaussianShell & | shell (int i) const |
Return a reference to GaussianShell number i. | |
GaussianShell & | shell (int i) |
Return a reference to GaussianShell number i. | |
const GaussianShell & | operator() (int icenter, int ishell) const |
Return a reference to GaussianShell number ishell on center icenter. | |
GaussianShell & | operator() (int icenter, int ishell) |
Return a reference to GaussianShell number ishell on center icenter. | |
const GaussianShell & | shell (int i, int j) const |
Return a reference to GaussianShell number j on center i. | |
GaussianShell & | shell (int i, int j) |
Return a reference to GaussianShell number j on center i. | |
double | r (int icenter, int xyz) const |
The location of center icenter. | |
int | values (const SCVector3 &r, ValueData *, double *basis_values) const |
Compute the values for this basis set at position r. | |
int | grad_values (const SCVector3 &r, ValueData *, double *g_values, double *basis_values=0) const |
Like values(. | |
int | hessian_values (const SCVector3 &r, ValueData *, double *h_values, double *g_values=0, double *basis_values=0) const |
Like values(. | |
int | shell_values (const SCVector3 &r, int sh, ValueData *, double *basis_values) const |
Compute the values for the given shell functions at position r. | |
int | grad_shell_values (const SCVector3 &r, int sh, ValueData *, double *g_values, double *basis_values=0) const |
Like values(. | |
int | hessian_shell_values (const SCVector3 &r, int sh, ValueData *, double *h_values, double *g_values=0, double *basis_values=0) const |
Like values(. | |
int | equiv (const Ref< GaussianBasisSet > &b) |
Returns true if this and the argument are equivalent. | |
void | print_brief (std::ostream &=ExEnv::out0()) const |
Print a brief description of the basis set. | |
void | print (std::ostream &=ExEnv::out0()) const |
Print a detailed description of the basis set. | |
Protected Member Functions | |
GaussianBasisSet (const GaussianBasisSet &) | |
virtual void | set_matrixkit (const Ref< SCMatrixKit > &) |
Classes | |
class | ValueData |
This holds scratch data needed to compute basis function values. More... |
Inputs for common basis sets are included in the MPQC distribution. They have been obtained from the EMSL Basis Set Database and translated into the MPQC format. The citation for this database is below. The technical citation for each basis set is listed in the individual basis set data files, in MPQC's lib/basis
directory.
Following is a table with available basis sets listing the supported elements for each basis and the number of basis functions for H, , first row,
, and second row,
, atoms. Basis sets with non-alpha-numerical characters in their name must be given in quotes.
Basis Set | Elements | ![]() | ![]() | ![]() |
STO-2G | H-Ca | 1 | 5 | 9 |
STO-3G | H-Kr | 1 | 5 | 9 |
STO-3G* | H-Ar | 1 | 5 | 14 |
STO-6G | H-Kr | 1 | 5 | 9 |
MINI (Huzinaga) | H-Ca | 1 | 5 | 9 |
MINI (Scaled) | H-Ca | 1 | 5 | 9 |
MIDI (Huzinaga) | H-Na, Al-K | 2 | 9 | 13 |
DZ (Dunning) | H, Li, B-Ne, Al-Cl | 2 | 10 | 18 |
DZP (Dunning) | H, Li, B-Ne, Al-Cl | 5 | 16 | 24 |
DZP + Diffuse (Dunning) | H, B-Ne | 6 | 19 | |
3-21G | H-Kr | 2 | 9 | 13 |
3-21G* | H-Ar | 2 | 9 | 19 |
3-21++G | H-Ar | 3 | 13 | 17 |
3-21++G* | H-Ar | 3 | 13 | 23 |
4-31G | H-Ne, P-Cl | 2 | 9 | 13 |
6-31G | H-Zn | 2 | 9 | 13 |
6-31G* | H-Zn | 2 | 15 | 19 |
6-31G** | H-Zn | 5 | 15 | 19 |
6-31+G* | H-Ar | 2 | 19 | 23 |
6-31++G | H-Ca | 3 | 13 | 17 |
6-31++G* | H-Ar | 3 | 19 | 23 |
6-31++G** | H-Ar | 6 | 19 | 23 |
6-311G | H-Ca, Ga-Kr | 3 | 13 | 21 |
6-311G* | H-Ca, Ga-Kr | 3 | 18 | 26 |
6-311G** | H-Ca, Ga-Kr | 6 | 18 | 26 |
6-311G(2df,2pd) | H-Ne, K, Ca | 14 | 30 | |
6-311++G** | H-Ne | 7 | 22 | |
6-311++G(2d,2p) | H-Ca | 10 | 27 | 35 |
6-311++G(3df,3pd) | H-Ar | 18 | 39 | 47 |
cc-pVDZ | H-Ar, Ca, Ga-Kr | 5 | 14 | 18 |
cc-pVTZ | H-Ar, Ca, Ga-Kr | 14 | 30 | 34 |
cc-pVQZ | H-Ar, Ca, Ga-Kr | 30 | 55 | 59 |
cc-pV5Z | H-Ar, Ca, Ga-Kr | 55 | 91 | 95 |
cc-pV6Z | H, He, B-Ne, Al-Ar | 91 | 140 | 144 |
aug-cc-pVDZ | H, He, B-Ne, Al-Ar, Ga-Kr | 9 | 23 | 27 |
aug-cc-pVTZ | H, He, B-Ne, Al-Ar, Ga-Kr | 23 | 46 | 50 |
aug-cc-pVQZ | H, He, B-Ne, Al-Ar, Ga-Kr | 46 | 80 | 84 |
aug-cc-pV5Z | H, He, B-Ne, Al-Ar, Ga-Kr | 80 | 127 | 131 |
aug-cc-pV6Z | H, He, B-Ne, Al-Ar | 127 | 189 | 193 |
cc-pCVDZ | Li, B-Ar | 18 | 27 | |
cc-pCVTZ | Li, B-Ar | 43 | 59 | |
cc-pCVQZ | Li, B-Ar | 84 | 109 | |
cc-pCV5Z | B-Ne | 145 | ||
aug-cc-pCVDZ | B-F, Al-Ar | 27 | 36 | |
aug-cc-pCVTZ | B-Ne, Al-Ar | 59 | 75 | |
aug-cc-pCVQZ | B-Ne, Al-Ar | 109 | 134 | |
aug-cc-pCV5Z | B-F | 181 | ||
NASA Ames ANO | H, B-Ne, Al, P, Ti, Fe, Ni | 30 | 55 | 59 |
pc-0 | H, C-F, Si-Cl | 2 | 9 | 13 |
pc-1 | H, C-F, Si-Cl | 5 | 14 | 18 |
pc-2 | H, C-F, Si-Cl | 14 | 30 | 34 |
pc-3 | H, C-F, Si-Cl | 34 | 64 | 64 |
pc-4 | H, C-F, Si-Cl | 63 | 109 | 105 |
pc-0-aug | H, C-F, Si-Cl | 3 | 13 | 17 |
pc-1-aug | H, C-F, Si-Cl | 9 | 23 | 27 |
pc-2-aug | H, C-F, Si-Cl | 23 | 46 | 50 |
pc-3-aug | H, C-F, Si-Cl | 50 | 89 | 89 |
pc-4-aug | H, C-F, Si-Cl | 88 | 145 | 141 |
All basis sets except for the pc-n and pc-n-aug basis sets were obtained from the Extensible Computational Chemistry Environment Basis Set Database, Version 12/03/03, as developed and distributed by the Molecular Science Computing Facility, Environmental and Molecular Sciences Laboratory which is part of the Pacific Northwest Laboratory, P.O. Box 999, Richland, Washington 99352, USA, and funded by the U.S. Department of Energy. The Pacific Northwest Laboratory is a multi-program laboratory operated by Battelle Memorial Institute for the U.S. Department of Energy under contract DE-AC06-76RLO 1830. Contact David Feller or Karen Schuchardt for further information.
The pc-n and pc-n-aug basis sets are the polarization consistent basis sets of Frank Jensen. See J. Chem. Phys. 115 (2001) 9113; J. Chem. Phys. 116 (2002) 7372; J. Chem. Phys. 117 (2002) 9234; and J. Chem. Phys. 121 (2004) 3463.
The KeyVal constructor.
molecule
puream
name
basis
element
basisdir
basisfiles
matrixkit
Several files in various directories are checked for basis set data. First, basis sets can be given by the user in the basis section at the top level of the main input file. Next, if a path is given with the basisdir keyword, then all of the files given with the basisfiles keyword are read in after appending their names to the value of basisdir. Basis sets can be given in these files in the basis section at the top level as well. If the named basis set still cannot be found, then GaussianBasisSet will try convert the basis set name to a file name and check first in the directory given by basisdir. Next it checks for the environment variable SCLIBDIR. If it is set it will look for the basis file in $SCLIBDIR/basis. Otherwise it will look in the source code distribution in the directory SC/lib/basis. If the executable has changed machines or the source code has be moved, then it may be necessary to copy the library files to your machine and set the SCLIBDIR environmental variable.
The basis set itself is also given in the ParsedKeyVal format. There are two recognized formats for basis sets:
type
am
puream
exp
coef
basis: ( carbon: "STO-2G": [ (type: [(am = s)] { exp coef:0 } = { 27.38503303 0.43012850 4.87452205 0.67891353 }) (type: [(am = p) (am = s)] { exp coef:1 coef:0 } = { 1.13674819 0.04947177 0.51154071 0.28830936 0.96378241 0.61281990 }) ] )
am
and nprim
, which specify the angular momentum and the number of shells in each block of even-tempered primitives. In addition, one must provide any two of the following vectors:
first_exp
last_exp
exp_ratio
Here's an example of a basis set composed of 2 blocks of even-tempered s-functions and 1 block of even-tempered p-functions.
basis: ( neon: "20s5s13p":(
am = [ 0 0 1 ] nprim = [ 20 5 13 ] first_exp = [ 1000.0 0.1 70.0 ] last_exp = [ 1.0 0.01 0.1 ]
) )
sc::GaussianBasisSet::GaussianBasisSet | ( | UnitType | ) |
This can be given GaussianBasisSet::Unit to construct a basis set with a single basis function that is one everywhere.
This can be used with integral evaluators to compute certain classes of integrals, with limitations.
Ref<GaussianBasisSet> sc::GaussianBasisSet::operator+ | ( | const Ref< GaussianBasisSet > & | B | ) |
Returns a GaussianBasisSet object that consists of the basis functions for each atom in this followed by the basis functions in B for the corresponding atom.
The Molecule object for the two basis sets must be identical.
void sc::GaussianBasisSet::save_data_state | ( | StateOut & | ) | [virtual] |
Save the base classes (with save_data_state) and the members in the same order that the StateIn CTOR initializes them.
This must be implemented by the derived class if the class has data.
Reimplemented from sc::SavableState.
const char* sc::GaussianBasisSet::label | ( | ) | const [inline] |
int sc::GaussianBasisSet::shell_on_center | ( | int | icenter, | |
int | shell | |||
) | const |
Return an overall shell number, given a center and the shell number on that center.
int sc::GaussianBasisSet::max_ncartesian_in_shell | ( | int | aminc = 0 |
) | const |
Return the maximum number of Cartesian functions that any shell has.
The optional argument is an angular momentum increment.
int sc::GaussianBasisSet::max_am_for_contraction | ( | int | con | ) | const |
Return the maximum angular momentum found in the given contraction number for any shell.
double sc::GaussianBasisSet::r | ( | int | icenter, | |
int | xyz | |||
) | const |
The location of center icenter.
The xyz argument is 0 for x, 1 for y, and 2 for z.
int sc::GaussianBasisSet::values | ( | const SCVector3 & | r, | |
ValueData * | , | |||
double * | basis_values | |||
) | const |
Compute the values for this basis set at position r.
The basis_values argument must be vector of length nbasis.
int sc::GaussianBasisSet::grad_values | ( | const SCVector3 & | r, | |
ValueData * | , | |||
double * | g_values, | |||
double * | basis_values = 0 | |||
) | const |
Like values(.
..), but computes gradients of the basis function values, too. The g_values argument must be vector of length 3*nbasis. The data will be written in the order bf1_x, bf1_y, bf1_z, ...
int sc::GaussianBasisSet::hessian_values | ( | const SCVector3 & | r, | |
ValueData * | , | |||
double * | h_values, | |||
double * | g_values = 0 , |
|||
double * | basis_values = 0 | |||
) | const |
Like values(.
..), but computes first and second derivatives of the basis function values, too. h_values must be vector of length 6*nbasis. The data will be written in the order bf1_xx, bf1_yx, bf1_yy, bf1_zx, bf1_zy, bf1_zz, ...
int sc::GaussianBasisSet::shell_values | ( | const SCVector3 & | r, | |
int | sh, | |||
ValueData * | , | |||
double * | basis_values | |||
) | const |
Compute the values for the given shell functions at position r.
See the other values(...) members for more information.
int sc::GaussianBasisSet::grad_shell_values | ( | const SCVector3 & | r, | |
int | sh, | |||
ValueData * | , | |||
double * | g_values, | |||
double * | basis_values = 0 | |||
) | const |
Like values(.
..), but computes gradients of the shell function values, too. See the other grad_values(...) members for more information.
int sc::GaussianBasisSet::hessian_shell_values | ( | const SCVector3 & | r, | |
int | sh, | |||
ValueData * | , | |||
double * | h_values, | |||
double * | g_values = 0 , |
|||
double * | basis_values = 0 | |||
) | const |
Like values(.
..), but computes first and second derivatives of the shell function values, too. See the other hessian_values(...) members for more information.