5.5.3 Standard Quadrature Grids

Four different “standard grids" are available in Q-Chem, designated SG-$n$, for $n=0,1,2$, or 3; both quality and the computational cost of these grids increases with $n$. These grids are constructed starting from a “parent” grid ($N_r,N_{\mathrm{\Omega }}$) consisting of $N_r$ radial spheres with $N_{\mathrm{\Omega }}$ angular (Lebedev) grid points on each, then systematically pruning the number of angular points in regions where sophisticated angular quadrature is not necessary, such as near the nuclei where the charge density is nearly spherically symmetric and at long distance from the nuclei where it varies slowly. A large number of points are retained in the valence region where angular accuracy is critical. The SG-$n$ grids are summarized in Table 5.3. While many electronic structure programs use some kind of procedure to delete unnecessary grid points in the interest of computational efficiency, Q-Chem’s SG-$n$ grids are notable in that the complete grid specifications are available in the peer-reviewed literature, ³⁸⁶ Gill P. M. W., Johnson B. G., Pople J. A.
Chem. Phys. Lett.
(1993), 209, pp. 506. Link ^{, 207} Chien S.-H., Gill P. M. W.
J. Comput. Chem.
(2006), 27, pp. 730. Link ^{, 260} Dasgupta S., Herbert J. M.
J. Comput. Chem.
(2017), 38, pp. 869. Link to facilitate reproduction of Q-Chem DFT calculations using other electronic structure programs. Just as computed energies may vary quite strongly with the choice of basis set, so too in DFT they may vary strongly with the choice of quadrature grid. In publications, users should always specify the grid that is used, and it is suggested to cite the appropriate literature reference from Table 5.3.

The SG-0 and SG-1 grids are designed for calculations on large molecules using GGA functionals. SG-1 affords integration errors on the order of $\sim$0.2 kcal/mol for medium-sized molecules and GGA functionals, including for demanding test cases such as reaction enthalpies for isomerizations. (Integration errors in total energies are no more than a few $\mu$hartree, or $\sim$0.01 kcal/mol.) The SG-0 grid was derived in similar fashion, and affords a root-mean-square error in atomization energies of 72 $\mu$hartree with respect to SG-1, while relative energies are reproduced well. ²⁰⁷ Chien S.-H., Gill P. M. W.
J. Comput. Chem.
(2006), 27, pp. 730. Link In either case, errors of this magnitude are typically considerably smaller than the intrinsic errors in GGA energies, and hence acceptable. As seen in Table 5.3, SG-1 retains % of the grid points of its parent grid, which translates directly into cost savings.

Both SG-0 and SG-1 were optimized so that the integration error in the energy falls below a target threshold, but derivatives of the energy (including such properties as (hyper)polarizabilities ¹⁸³ Castet F., Champagne B.
J. Chem. Theory Comput.
(2012), 8, pp. 2044. Link ) are often more sensitive to the quality of the integration grid. Special care is required, for example, when imaginary vibrational frequencies are encountered, as low-frequency (but real) vibrational frequencies can manifest as imaginary if the grid is sparse. If imaginary frequencies are found, or if there is some doubt about the frequencies reported by Q-Chem, the recommended procedure is to perform the geometry optimization and vibrational frequency calculations again using a higher-quality grid. (The optimization should converge quite quickly if the previously-optimized geometry is used as an initial guess.)

SG-1 was the default DFT integration grid for all density functionals for Q-Chem versions 3.2–4.4. Beginning with Q-Chem v. 4.4.2, however, the default grid is functional-dependent, as summarized in Table 5.3. This is a reflection of the fact that although SG-1 is adequate for energy calculations using most GGA and hybrid functionals (although care must be taken for some other properties, as discussed below), it is not adequate to integrate many functionals developed since $\sim$2005. These include meta-GGAs, which are more complicated due to their dependence on the kinetic energy density (${\tau }_{\sigma }$ in Eq. (5.10)) and/or the Laplacian of the density (${\nabla }^2{\rho }_{\sigma }$). Functionals based on B97, along with the Minnesota suite of functionals, ¹³⁵⁶ Zhao Y., Truhlar D. G.
Theor. Chem. Acc.
(2008), 120, pp. 215. Link ^{, 1357} Zhao Y., Truhlar D. G.
Chem. Phys. Lett.
(2011), 502, pp. 1. Link contain relatively complicated expressions for the exchange inhomogeneity factor, and are therefore also more sensitive to the quality of the integration grid. ¹²⁷⁰ Wheeler S. E., Houk K. N.
J. Chem. Theory Comput.
(2010), 6, pp. 395. Link ^{, 787} Mardirossian N., Head-Gordon M.
Phys. Chem. Chem. Phys.
(2014), 16, pp. 9904. Link ^{, 260} Dasgupta S., Herbert J. M.
J. Comput. Chem.
(2017), 38, pp. 869. Link To integrate these modern density functionals, the SG-2 and SG-3 grids were developed, ²⁶⁰ Dasgupta S., Herbert J. M.
J. Comput. Chem.
(2017), 38, pp. 869. Link which are pruned versions of the medium-quality (75, 302) and high-quality (99, 590) integration grids, respectively. Tests of properties known to be highly sensitive to the quality of the integration grid, such as vibrational frequencies, hyper-polarizabilities, and potential energy curves for non-bonded interactions, demonstrate that SG-2 is usually adequate for meta-GGAs and B97-based functionals, and in many cases is essentially converged with respect to an unpruned (250, 974) grid. ²⁶⁰ Dasgupta S., Herbert J. M.
J. Comput. Chem.
(2017), 38, pp. 869. Link The Minnesota functionals are more sensitive to the grid, and while SG-3 is often adequate, it is not completely converged in the case of non-bonded interactions. ²⁶⁰ Dasgupta S., Herbert J. M.
J. Comput. Chem.
(2017), 38, pp. 869. Link

Note:  1. SG-0 was re-optimized for Q-Chem v. 3.0, so results may differ slightly as compared to older versions of the program. 2. The SG-2 and SG-3 grids use a double-exponential radial quadrature, ²⁶⁰ Dasgupta S., Herbert J. M.
J. Comput. Chem.
(2017), 38, pp. 869. Link whereas a general grid (selected by setting XC_GRID = $X\mathbf{}Y$, as described in Section 5.4) uses an Euler-MacLaurin radial quadrature. As such, absolute energies cannot be compared between, e.g., SG-2 and XC_GRID = 000075000302, even though SG-2 uses a pruned (75, 302) grid. However, energy differences should be quite similar between the two. 3. As noted in Ref.  260 Dasgupta S., Herbert J. M.
J. Comput. Chem.
(2017), 38, pp. 869. Link , for Minnesota functionals some wiggles in potential energy surfaces may persist with the SG-3 grid, especially for longer-range non-bonded interactions. Although these are rarely problematic for energy differences, if the user wants to eliminate these oscillations then we recommend an unpruned (99, 590) grid, i.e., XC_GRID = 000099000590.