6.12.11 Cholesky Decomposition with CC (CD-CC)

Two-electron integrals can be decomposed using Cholesky decomposition ³⁶⁶ Epifanovsky E. et al.
J. Chem. Phys.
(2013), 139, pp. 134105. Link ^{, 382} Feng X. et al.
J. Chem. Phys.
(2019), 151, pp. 014110. Link giving rise to the same representation as in RI and substantially reducing the cost of integral transformation, disk storage requirements, and improving parallel performance:

The rank of Cholesky decomposition, $M$, is typically 3-10 times larger than the number of basis functions $N$ (Ref.  49 Aquilante F., Pedersen T. B., Lindh R.
Theor. Chem. Acc.
(2009), 124, pp. 1. Link ); it depends on the decomposition threshold $\delta$ and is considerably smaller than the full rank of the matrix, $N(N+1)/2$ (Refs.  49 Aquilante F., Pedersen T. B., Lindh R.
Theor. Chem. Acc.
(2009), 124, pp. 1. Link , 103 Beebe N. H. F., Linderberg J.
Int. J. Quantum Chem.
(1977), 12, pp. 683. Link , 1478 Wilson S.
Comput. Phys. Commun.
(1990), 58, pp. 71–81. Link ). Cholesky decomposition removes linear dependencies in product densities $(\mu \nu |$, ⁴⁹ Aquilante F., Pedersen T. B., Lindh R.
Theor. Chem. Acc.
(2009), 124, pp. 1. Link allowing one to obtain compact approximation to the original matrix with accuracy, in principle, up to machine precision.

Decomposition threshold $\delta$ is the only parameter that controls accuracy and the rank of the decomposition. Cholesky decomposition is invoked by specifying CHOLESKY_TOL that defines the accuracy with which decomposition should be performed. For most calculations tolerance of $\delta ={10}^{-3}$ gives a good balance between accuracy and compactness of the rank. Tolerance of $\delta ={10}^{-2}$ can be used for exploratory calculations and $\delta ={10}^{-4}$ for high-accuracy calculations. Similar to RI, Cholesky-decomposed integrals can be transformed back, into the canonical MO form, using CC_DIRECT_RI keyword.

The default execution uses the new two-step algorithm for the decomposition. ¹²¹⁵ Schnack-Petersen A. K. et al.
J. Chem. Phys.
(2022), 156, pp. 244111. Link ^{, 937} Melega L. et al.
J. Chem. Phys.
(2026), 130, pp. 2097. Link The old algorithm ³⁶⁶ Epifanovsky E. et al.
J. Chem. Phys.
(2013), 139, pp. 134105. Link ^{, 382} Feng X. et al.
J. Chem. Phys.
(2019), 151, pp. 014110. Link can be invoked by specifying CHOLESKY_TWOSTEP = FALSE. The two-step algorithm uses additional parameter, CD_TS_SPAN_FACTOR, which is defaulted to 2. A larger value may be used to increase accuracy, as illustrated by Example 6.12.13.

Note:  Cholesky decomposition is available for all CCMAN2 methods, including energy, analytic gradients, and properties calculations. For maximum computational efficiency, combine with FNO (see Sections 6.15.1 and 7.9.14) when appropriate.

Examples 6.12.13, 6.12.13, and 6.12.13 illustrate the usage of this feature.