---
title: 'THC-sRI-CC2: Hybrid Reduced-Scaling CC2 Method'
url: https://www.emergentmind.com/topics/thc-sri-cc2
type: topic
---

# THC-sRI-CC2: Hybrid Reduced-Scaling CC2 Method

THC-sRI-CC2 is a noise-reduced formulation of stochastic resolution of identity applied to CC2 in which tensor hypercontraction (THC) is used to decouple the expensive Coulomb term and stochastic resolution of the identity (sRI) is used to treat the time-determining exchange term, while preserving an overall $O(N^3)$ scaling with system size $N$ [2509.21885]. The method is designed to mitigate the substantial stochastic errors that conventional sRI-CC2 can introduce in energy derivatives, particularly in contexts such as molecular dynamics simulations, and benchmarks report greater accuracy and markedly reduced stochastic noise than conventional sRI-CC2 at identical computational samplings [2509.21885].

## 1. Definition and computational rationale

The sRI approximation reduces the computational scaling of CC2 from $O(N^5)$ to $O(N^3)$, where $N$ is a measure of system size, but its inherent stochastic noise can introduce substantial errors in energy derivatives [2509.21885]. THC-sRI-CC2 addresses this limitation by combining two distinct low-scaling decompositions within the same CC2 workflow.

In this formulation, the expensive Coulomb term, which scales as $O(N^4)$, is decoupled via THC, while the time-determining exchange term with an $O(N^5)$ cost is addressed through the sRI scheme, collectively yielding an overall $O(N^3)$ scaling [2509.21885]. The central design choice is therefore not a uniform stochastic treatment of all contractions, but a hybrid partitioning in which the deterministic THC factorization removes the dominant source of noise from the Coulomb contribution and leaves stochastic sampling only in the exchange contribution.

This suggests that the defining feature of THC-sRI-CC2 is not merely reduced asymptotic cost, but a redistribution of approximation error: deterministic low-rank compression for the Coulomb structure and tunable stochastic estimation for the exchange structure.

## 2. CC2 equations and the key contraction

The method is built from the CC2 working equations. With $|HF\rangle$ denoting the Hartree–Fock determinant, the CC cluster operator is defined as
$$
T = T_1 + T_2,
$$
with
$$
T_1 = \sum_{ia} t_i^a\,E_{ai},
$$
and
$$
T_2 = \frac{1}{4}\sum_{ijab} t_{ij}^{ab}\,E_{ai}\,E_{bj}.
$$

The CC2 ground-state energy and amplitude equations follow from the similarity-transformed Hamiltonian
$$
\bar H = e^{-T_1} H e^{T_1}.
$$
Projecting onto singles and doubles yields
$$
E_{CC2} = \langle HF|\,\bar H + [\bar H,T_2]\,|HF\rangle,
$$
$$
0 = \Omega_{i}^{a} \equiv \langle a_i|\,\bar H + [\bar H,T_2]\,|HF\rangle,
$$
and
$$
0 = \Omega_{ij}^{ab} \equiv \langle a_i b_j|\,\bar H + [F,T_2]\,|HF\rangle.
$$

Once the doubles amplitudes $t_{ij}^{ab}$ are found by solving the doubles equation, one back-substitutes to obtain the singles amplitudes and then the CC2 energy [2509.21885]. In practical form, the key tensor contraction appearing in $\Omega_{ij}^{ab}$ or in excited-state Jacobian contractions is
$$
A_{ik} = \sum_{abj} t_{ij}^{ab}\,(2\,(ai|bj) - (bi|aj))\,(bj|ak).
$$

This contraction separates naturally into a Coulomb-like part and an exchange-like part. The hybrid structure of THC-sRI-CC2 follows directly from this decomposition: the Coulomb contribution is factorized deterministically, while the exchange contribution is estimated stochastically.

## 3. Tensor hypercontraction for the Coulomb contribution

Tensor HyperContraction factorizes a four-index electron-repulsion quantity $(\alpha\beta|\gamma\delta)$ by introducing four low-rank interpolation matrices and a small core matrix:
$$
(\alpha\beta|\gamma\delta) \approx \sum_{K,L} A_{\alpha K}\,B_{\beta K}\,M_{KL}\,C_{\gamma L}\,D_{\delta L}.
$$
A common choice, identified as Interpolative Separable Density Fitting, is
$$
A_{\alpha K}=\phi^K_\alpha,\quad B_{\beta K}=\phi^K_\beta,\quad C_{\gamma L}=\phi^L_\gamma,\quad D_{\delta L}=\phi^L_\delta,
$$
so that
$$
(\alpha\beta|\gamma\delta) \approx \sum_{K,L} \phi^K_\alpha\,\phi^K_\beta\,M_{KL}\,\phi^L_\gamma\,\phi^L_\delta.
$$

All five tensor factors $\{\phi^K_\mu,M_{KL}\}$ scale as $O(N^2)$ elements, and $N_{IP}\sim c\cdot N$ [2509.21885]. Inserting this factorization into the Coulomb part of the key CC2 contraction and regrouping yields an $O(N^3)$ form:
$$
C_{ik}^C = 2\sum_{K,N} X_{KN}\left(\sum_a \phi^K_a\,\phi^N_a\right)\phi^K_i\,\phi^N_k,
$$
where
$$
X_{KN} = \sum_{L,M} \left(\sum_b \phi^L_b\,\phi^M_b\right)\left(\sum_j \phi^L_j\,\phi^M_j\right)M_{KL}\,M_{MN},
$$
and $i,k$ run over occupied orbitals while $a,b$ and $j$ are active indices.

The stated implementation precomputes the $N_{IP}\times N_{IP}$ intermediates
$$
B_{LM}=\sum_b \phi^L_b\,\phi^M_b,\qquad J_{LM}=\sum_j \phi^L_j\,\phi^M_j,
$$
then
$$
Y_{KN}=\sum_L M_{KL}\,J_{LM},\qquad X_{KN}=\sum_M Y_{KM}\,B_{MN},
$$
so that all sums reduce to $O(N^3)$ steps [2509.21885]. The significance of this deterministic treatment is that the Coulomb part no longer contributes stochastic variance.

## 4. Stochastic resolution of the identity for the exchange contribution

The exchange contraction in the key CC2 tensor expression involves two distinct four-index integrals and is formally $O(N^5)$ [2509.21885]. THC-sRI-CC2 treats this part with stochastic resolution of the identity.

The sRI procedure draws $N_s$ random auxiliary-space vectors $\theta^\xi\in\{\pm 1\}^{N_{\mathrm aux}}$, for $\xi=1\ldots N_s$, satisfying
$$
E_\xi[\theta^\xi(\theta^\xi)^T] = I_{\mathrm aux}.
$$
For each $\xi$, one builds the stochastic RI tensor
$$
R^{\xi}_{\alpha\beta} = \sum_{P,Q} (\alpha\beta|P)\,V^{-1/2}_{PQ}\,\theta^\xi_Q,
$$
so that
$$
(\alpha\beta|\gamma\delta) \approx E_\xi[\,R^{\xi}_{\alpha\beta}\,R^{\xi}_{\gamma\delta}\,]
= \frac{1}{N_s}\sum_\xi R^{\xi}_{\alpha\beta}\,R^{\xi}_{\gamma\delta}.
$$

The exchange contribution is then written as
$$
C^X_{ik} = - E_{\xi,\xi'} \left[ \sum_{bj} (R^{\xi}_{bi}\,R^{\xi}_{aj}) (R^{\xi'}_{bj}\,R^{\xi'}_{ak}) \right],
$$
which, after rearranging inner sums over $b,j,a,i,k$, is realized in $O(N_s\,N^3)$ [2509.21885].

Variance analysis is given explicitly. For each fixed index tuple $(b,i,a,j)$,
$$
\mathrm{Var}[\,R^{\xi}_{bi}\,R^{\xi}_{aj}\,] = E[(R_{bi}R_{aj})^2] - (bi|aj)^2,
$$
and by standard sampling theory
$$
\mathrm{Var}[\,C^X_{ik}\,] = O(1/N_s).
$$
Thus the per-element stochastic error is $\sim O(N_s^{-1/2})$ and is tunable by $N_s$ [2509.21885]. In this respect, THC-sRI-CC2 retains the controllable-error character of sRI while confining it to the exchange sector.

## 5. Scaling and noise characteristics

The asymptotic cost analysis is stated in terms of a system-size measure $N$, such as the number of basis functions or electrons [2509.21885]. For the THC steps associated with the Coulomb term, building $B_{LM}$ and $J_{LM}$ scales as
$$
O(N_{\mathrm occ}\,N_{IP}^2 + N_{\mathrm vir}\,N_{IP}^2)=O(N^3),
$$
building $Y_{KN}$ and $X_{KN}$ scales as
$$
O(N_{IP}^3)=O(N^3),
$$
and the final contraction scales as
$$
O(N^3).
$$

For the sRI steps associated with the exchange term, forming $R^\xi_{\alpha\beta}$ costs
$$
O(N_s\,N_{\mathrm aux}\,N_{\alpha\beta})=O(N^3)
$$
since $N_s\sim\mathrm{const}$, and the exchange contraction scales as
$$
O(N_s\,N_{\mathrm occ}^2\,N_{\mathrm vir})=O(N^3).
$$
Because $N_s$, $N_{\mathrm aux}$, and $N_{IP}$ are all proportional to $N$ and treated as linear in $N$, every stage grows as $N^3$, giving overall $O(N^3)$ scaling [2509.21885].

The noise-reduction argument is correspondingly explicit. In pure sRI-CC2, both Coulomb and exchange carry stochastic error $\sim O(N_s^{-1/2})$. In THC-sRI-CC2, the Coulomb part is entirely deterministic, so the only fluctuation stems from the much smaller exchange contribution. Denoting
$$
\sigma_C^2 = \mathrm{Var}_{\mathrm Coulomb},\qquad \sigma_X^2 = \mathrm{Var}_{\mathrm Exchange},
$$
one has in sRI-CC2
$$
\sigma_{\mathrm total}^2\approx \sigma_C^2+\sigma_X^2,
$$
whereas in THC-sRI-CC2
$$
\sigma_{\mathrm total}^2\approx \sigma_X^2.
$$
Numerically $\sigma_C\gg\sigma_X$, so eliminating the Coulomb noise often reduces the total $\sigma$ by an order of magnitude. Quantitatively,
$$
\epsilon_{\mathrm total}\equiv \mathrm{std.dev.}(\mathrm{energy}) \simeq \frac{\sigma_X}{\sqrt{N_s}} \ll \frac{\sqrt{\sigma_C^2+\sigma_X^2}}{\sqrt{N_s}}.
$$

A plausible implication is that the principal practical gain of the method is not only asymptotic efficiency, but also improved stability of derivative-level quantities at fixed sampling.

## 6. Benchmarks and stated scope

The reported benchmarks compare conventional sRI-CC2, THC-sRI-CC2, and RI-CC2 [2509.21885].

| Quantity | sRI-CC2 | THC-sRI-CC2 |
|---|---:|---:|
| H$_2$O excitation energy, $N_s=5\,000$ | $8.0500\pm0.1986$ eV | $8.0965\pm0.0330$ eV |
| Benzene excitation energy, $N_s=5\,000$ | $5.3048\pm0.0951$ eV | $5.3768\pm0.0288$ eV |
| H$_2$O $\Delta_{\max}$ ground-state gradient, $N_s=50\,000$ | $0.0028$ | $0.0017$ |
| Benzene $\Delta_{\max}$ ground-state gradient, $N_s=50\,000$ | $0.1200$ | $0.0147$ |

For the same excitation-energy benchmark, the RI-CC2 reference values are $8.0924$ eV for H$_2$O and $5.3668$ for benzene [2509.21885]. The benchmark summary further states “$\sigma$ reduced by $\sim 5\times$” and “abs. error $<0.01$ eV” for excitation energies. For ground-state gradients, the summary states “$5$–$10\times$ noise reduction,” and excited-state gradients and derivative couplings are reported to exhibit similar $5$–$10\times$ drops in stochastic error at identical $N_s$ [2509.21885].

Scaling timings on all-E olefins are reported as
$$
\text{RI-CC2: } t\propto N^{4.3},\qquad
\text{sRI-CC2: } t\propto N^{2.8},\qquad
\text{THC-sRI-CC2: } t\propto N^{2.8},
$$
with crossover around $N_e\approx 70$–$120$ electrons, beyond which the cubic methods become faster by factors of $2$–$5$ at $N_e\gg 100$ [2509.21885].

The stated implications for application are direct. By combining the mathematical rigor of THC for Coulomb decoupling with the flexibility of sRI for exchange, THC-sRI-CC2 delivers full $O(N^3)$ scaling for energies, gradients, and derivative couplings; noise levels comparable to pure RI-CC2 with $5$–$10\times$ fewer stochastic samples; and wall-clock times $2$–$5\times$ faster than RI-CC2 for $N_e\gtrsim 100$ [2509.21885]. The same source states that this opens the door to routine CC2-quality excited-state and nonadiabatic dynamics in systems of several hundred atoms, with systematically controllable stochastic error and no loss in formal accuracy, and further establishes a general THC-sRI hybrid strategy for the development of reduced-scaling electronic structure methods.

Source: https://www.emergentmind.com/topics/thc-sri-cc2