14.5 Primitive set definition

A group of basis functions is defined by a data card specifying a set of primitive gaussians, optionally followed by one or more cards specifying particular contractions of primitives to be included in the final basis. When all contraction definitions have been read (delimited by the next data card other than a contraction definition), the remaining primitives in the set which have not been included in any contraction set are added uncontracted to the basis set.

There are four different input forms, as explained below under a) to d):

a) type,atom,key,scale,nprim;

Load basis named key from the library with angular symmetry type (S, P, D, F, G, H, or I). This basis is added from all atoms with number i=atom on the A cards (see above). If scale is present, all exponents are scaled by scale**2. If nprim is specified, the first nprim exponents only are taken from the library. If nprim is negative, the last $\vert nprim\vert$ basis functions from the library set are deleted. Associated with the library basis may be a set of default contraction coefficients which may be accessed in subsequent contraction cards. type can include several types, e.g., SPD or DF. This usually makes sense only with default contractions, i.e., such cards should be followed only by ``C'' without any other specifications for contractions.

b) type,atom,exp1,exp2,$\ldots$expn;expn+1,$\ldots$;

General specification of exponents; continuation onto subsequent cards (separated by semicolon) is permitted as shown (the first card can hold up to 19 exponents (cray 13), each following card 20 exponents (cray 15).

The exponents (and other numerical parameters described below such as numbers of functions, and contraction coefficients) can be given as general input expressions, possibly involving variables. It is important to note, however, that these expressions are evaluated typically just once, at the same time as the complete basis set is parsed. This generally happens the first time that the basis set is required, perhaps before the first SCF calculation can be done. If the variables on which the basis depends are altered, this will not be noticed by the program, and the new basis set will not be used for subsequent stages of the computation. If, however, a new basis block is presented in the input, then the program marks as outdated any quantities such as integrals that have been calculated with the old basis set; subsequent job steps will then use the new basis. Such behaviour can be forced at any stage by issuing the command

INT

and this is the recommended way of ensuring that a new basis set is adopted; note, however, that INT will cause atomic orbital integrals to be evaluated, and this overhead should of course not be incurred unnecessarily in non-direct calculations.

c) type,atom,EVEN,nprim,ratio,centre,dratio

Generates a generalized even tempered set of functions. The number of functions $n$ is specified by nprim, their geometric mean $c$ by centre, the mean ratio of successive exponents $r$ by ratio, and the variation of this ratio, $d$, by dratio. If centre is not given, the previous basis of the same type is extended by diffuse functions. If in this case ratio is not given, $r$ is determined from the exponents of the last two function of the previous basis. If this is not possible, the default $r=2.5$ is adopted. $d=1$ (the default) specifies a true even-tempered set, but otherwise the ratio between successive exponents changes linearly; the exponents are given explicitly by

$$\log e_i = \log c + ((n+1)/2-i) \,\log r + \frac12 ((n+1)/2-i)^2 \,\log d \quad i=1,2,\dots,n$$

Example 1

SP,1,VTZ;C;SP,1,EVEN,1;
generates the generally contracted $s$ and $p$ triple-zeta basis sets for atom 1 and extends these by one diffuse function.

Example 2

SPD,1,VTZ,,$-1$;C;SP,1,EVEN,2,2.5;
generates the generally contracted $s$, $p$ triple-zeta basis sets for atom 1. Two energy optimized $d$-functions of Dunning are included. The last $s$ and $p$ functions are deleted and replaced by two even tempered functions with ratio 2.5.

d) type,atom,EVENR,nprim,aa,ap,bb,bp

Generates an even tempered set of nprim functions according to the ``regular'' prescription described in M W Schmidt and K Ruedenberg, J. Chem. Phys. 71 (1970) 3951. If any of the parameters aa, ap, bb, bp is zero or omitted, the values are taken from table III of the above.