Next: 18.4.4 Atom partitioning of Up: 18.4 Numerical integration grid Previous: 18.4.2 Radial integration grid

18.4.3 Angular integration grid (ANGULAR)

ANGULAR,method,acca,crowd
LMIN, $l^{\text{min}}_0,l^{\text{min}}_1,l^{\text{min}}_2,l^{\text{min}}_3$
LMAX, $l^{\text{max}}_0,l^{\text{max}}_1,l^{\text{max}}_2,l^{\text{max}}_3$

Specify the details of the angular quadrature scheme. The default choice for method is LEBEDEV (ie. as in A. D. Becke, J. Chem. Phys. 88 (1988) 2547) which provides angular grids of octahedral symmetry. The alternative choice for method is LEGENDRE which gives Gauss-Legendre quadrature in $\theta$ and simple quadrature in $\phi$, as defined by C. W. Murray, N. C. Handy and G. J. Laming, Mol. Phys. 78 (1993) 997.

Each type of grid specifies a family of which the various members are characterized by a single quantum number $l$; spherical harmonics up to degree $l$ are integrated exactly. $l{\text{min}}_i$ and $l{\text{max}}_i, i=0,1,2,3$ specify allowed ranges of $l$ for hydrogen/helium, first row, second row, and other elements respectively. For the Lebedev grids, if the value of $l$ is not one of the set implemented in MOLPRO (3, 5, 7, 9, 11, 13, 15, 17, 19, 23, 29, 41, 47, 53), then $l$ is increased to give the next largest angular grid available. In general, different radial points will have different $l$, and in the absence of any moderation described below, will be taken from $l^{\text{max}}_i$.

crowd is a parameter to control the reduction of the degree of quadrature close to the nucleus, where points would otherwise be unnecessarily close together; larger values of crowd mean less reduction thus larger grids. A very large value of this parameter, or, conventionally, setting it c;to zero, will switch off this feature.

acca is a target energy accuracy. It is used to reduce $l$ for a given radial point as far as possible below $l^{\text{max}}_i$ but not lower than $l^{\text{max}}_i$. The implementation uses the error in the angular integral of the kernel of the Slater-Dirac exchange functional using a sum of approximate atomic densities. If acca is zero, the global threshold is used instead, or else it is ignored.



Next: 18.4.4 Atom partitioning of Up: 18.4 Numerical integration grid Previous: 18.4.2 Radial integration grid

P.J. Knowles and H.-J. Werner
molpro@tc.bham.ac.uk
Jan 15, 2002