Generalized Cluster-Correlation Expansion (gCCE)

• gCCE is an extension of the standard CCE that includes the central spin in every cluster Hamiltonian for refined modeling of quantum spin bath dynamics.
• It has been successfully applied to pure dephasing, multi-level coherence, and clock transition regimes, delivering quantitatively accurate predictions in controlled settings.
• However, the multiplicative ansatz of gCCE can fail for dipolar relaxation due to overlapping energy-transfer pathways, highlighting the need for additive reorganization in these cases.

Searching arXiv for the cited papers to ground the article in the current literature. Generalized Cluster-Correlation Expansion (gCCE) is an extension of the Cluster-Correlation Expansion (CCE) for central-spin dynamics in a quantum spin bath that retains the central spin inside each cluster Hamiltonian, so that the reduced density matrix elements of the central system are assembled from irreducible cluster contributions computed on the joint central-spin-plus-cluster Hilbert space. In the literature, gCCE has been used in distinct senses: as a population-based extension of CCE for longitudinal relaxation and other driven population-transfer problems (Yang et al., 2018, Fazio et al., 2024); as a multi-level coherence formalism for coupled-spin qubits in the pure-dephasing regime (Chen et al., 2024); and as a modified construction for clock transitions that absorbs quadratic Overhauser effects into configuration-dependent central-spin energies to restore cluster locality (Zhang et al., 2020). A recent mathematical analysis of dipolar relaxation shows that, in its standard product form, gCCE is capable of simulating the transfer of energy from the central spin into the bath in principle, but is insufficient for providing even a qualitatively accurate description of spin-spin relaxation in long-range dipolar baths (Ryan et al., 16 Feb 2026).

1. Definition and scope

The original CCE reorganizes the central-spin coherence $L(t)$ into a product over irreducible bath-cluster contributions, and for pure dephasing it can project the central spin out and treat bath dynamics conditional on the central-spin state (Ryan et al., 16 Feb 2026). In that setting, the factorization over bath clusters captures the additivity of phase fluctuations in the exponent and delivers excellent quantitative predictions for $T_2$ in many platforms (Ryan et al., 16 Feb 2026).

The generalized construction keeps the central spin explicitly in each cluster’s dynamics. In the formulation emphasized for dipolar central-spin problems, the cluster-reduced matrix element is

$$\rho_{ij,C}(t) = \langle i| e^{-i\hat H_C t} \rho_{S+C}(0) e^{+i\hat H_C t} |j\rangle,$$

and the irreducible contribution is defined by inclusion–exclusion,

$$\tilde\rho_{ij,C}(t) = \rho_{ij,C}(t) \Big/ \prod_{C' \subset C} \tilde\rho_{ij,C'}(t),$$

so that the truncated approximation is

$$\rho_{ij}^{(M)}(t) = \prod_{|C|\le M} \tilde\rho_{ij,C}(t).$$

This construction is stated explicitly in the dipolar-relaxation analysis (Ryan et al., 16 Feb 2026), in the two-qubit pure-dephasing formalism (Chen et al., 2024), and in PyCCE’s implementation-oriented overview (Onizhuk et al., 2021).

In population-based applications, the same product logic is applied to a survival probability or target-state population rather than to an off-diagonal coherence. For the NV-center longitudinal-relaxation problem, the observable is the survival probability $P(t)$ of $|0\rangle$, written as

$$P^{(M)}(t)=\prod_{|\mathfrak c|\le M}\tilde P_{\mathfrak c}(t),$$

with the NV spin included in every cluster (Yang et al., 2018). For STIRAP in a five-state system coupled to a nuclear bath, the generalized product expansion is written as

$$P(t)=\tilde P_{\{0\}}\prod_i \tilde P_{\{i\}}\prod_{i,j}\tilde P_{\{i,j\}}\cdots,$$

with time-dependent cluster evolution under a restricted driven Hamiltonian (Fazio et al., 2024).

A separate, modified/generalized CCE was introduced for clock transitions, where the difficulty is not energy exchange but the failure of standard CCE near quadratic-noise points. There the central-spin Hamiltonian is diagonalized for each bath eigenstate of the longitudinal hyperfine operator, and the long-range hyperfine-mediated terms are absorbed into configuration-dependent central energies $E_M(\delta)$ rather than treated as explicit bath-bath couplings (Zhang et al., 2020). This suggests that “gCCE” in the literature is not a single unique recipe, but a family of central-spin-inclusive or central-spin-conditioned cluster expansions adapted to regimes where standard bath-only CCE is inadequate.

2. Hamiltonian structure and observables

A general central-spin Hamiltonian used in the dipolar-relaxation analysis is decomposed as

$T_2$0

with tensor form

$T_2$1

where $T_2$2 and $T_2$3 are zero-field splitting and quadrupole tensors, $T_2$4 is the central-bath coupling tensor, $T_2$5 is the bath-bath dipolar coupling tensor, and $T_2$6 are gyromagnetic tensors (Ryan et al., 16 Feb 2026). A representative magnetic dipolar interaction between spins at separation $T_2$7 is

$T_2$8

and in a strong static field one ხშირად uses the secular form

$T_2$9

(Ryan et al., 16 Feb 2026).

For a bath cluster $$\rho_{ij,C}(t) = \langle i| e^{-i\hat H_C t} \rho_{S+C}(0) e^{+i\hat H_C t} |j\rangle,$$0, the generalized cluster Hamiltonian is

$$\rho_{ij,C}(t) = \langle i| e^{-i\hat H_C t} \rho_{S+C}(0) e^{+i\hat H_C t} |j\rangle,$$1

so that non-cluster spins are included at mean-field level via static averages $$\rho_{ij,C}(t) = \langle i| e^{-i\hat H_C t} \rho_{S+C}(0) e^{+i\hat H_C t} |j\rangle,$$2 (Ryan et al., 16 Feb 2026). PyCCE presents the same structural form for a general spin Hamiltonian and notes that for a fully mixed bath the mean-field terms vanish because $$\rho_{ij,C}(t) = \langle i| e^{-i\hat H_C t} \rho_{S+C}(0) e^{+i\hat H_C t} |j\rangle,$$3 (Onizhuk et al., 2021).

The relevant observable depends on regime. For relaxation, the central spin is initially polarized along $$\rho_{ij,C}(t) = \langle i| e^{-i\hat H_C t} \rho_{S+C}(0) e^{+i\hat H_C t} |j\rangle,$$4, for example $$\rho_{ij,C}(t) = \langle i| e^{-i\hat H_C t} \rho_{S+C}(0) e^{+i\hat H_C t} |j\rangle,$$5, and the population $$\rho_{ij,C}(t) = \langle i| e^{-i\hat H_C t} \rho_{S+C}(0) e^{+i\hat H_C t} |j\rangle,$$6 is followed (Ryan et al., 16 Feb 2026). For dephasing, the central spin is initialized in a superposition and the observable is the coherence $$\rho_{ij,C}(t) = \langle i| e^{-i\hat H_C t} \rho_{S+C}(0) e^{+i\hat H_C t} |j\rangle,$$7, $$\rho_{ij,C}(t) = \langle i| e^{-i\hat H_C t} \rho_{S+C}(0) e^{+i\hat H_C t} |j\rangle,$$8 (Ryan et al., 16 Feb 2026). In the pure-dephasing projection, when $$\rho_{ij,C}(t) = \langle i| e^{-i\hat H_C t} \rho_{S+C}(0) e^{+i\hat H_C t} |j\rangle,$$9, the total Hamiltonian reduces to

$$\tilde\rho_{ij,C}(t) = \rho_{ij,C}(t) \Big/ \prod_{C' \subset C} \tilde\rho_{ij,C'}(t),$$0

and

$$\tilde\rho_{ij,C}(t) = \rho_{ij,C}(t) \Big/ \prod_{C' \subset C} \tilde\rho_{ij,C'}(t),$$1

for a bath basis $$\tilde\rho_{ij,C}(t) = \rho_{ij,C}(t) \Big/ \prod_{C' \subset C} \tilde\rho_{ij,C'}(t),$$2 (Ryan et al., 16 Feb 2026).

Two-qubit gCCE adopts the same pure-dephasing logic in a multilevel central system. For an off-diagonal reduced-density-matrix element $$\tilde\rho_{ij,C}(t) = \rho_{ij,C}(t) \Big/ \prod_{C' \subset C} \tilde\rho_{ij,C'}(t),$$3, the normalized coherence is

$$\tilde\rho_{ij,C}(t) = \rho_{ij,C}(t) \Big/ \prod_{C' \subset C} \tilde\rho_{ij,C'}(t),$$4

which is independent of the initial central state provided the initial coherence is nonzero (Chen et al., 2024).

3. Product expansion, inclusion–exclusion, and cluster truncation

The algebraic core of gCCE is an inclusion–exclusion or Möbius-type construction over bath-spin subsets. For standard coherence CCE, one writes

$$\tilde\rho_{ij,C}(t) = \rho_{ij,C}(t) \Big/ \prod_{C' \subset C} \tilde\rho_{ij,C'}(t),$$5

and truncation by cluster size $$\tilde\rho_{ij,C}(t) = \rho_{ij,C}(t) \Big/ \prod_{C' \subset C} \tilde\rho_{ij,C'}(t),$$6 gives CCE-$$\tilde\rho_{ij,C}(t) = \rho_{ij,C}(t) \Big/ \prod_{C' \subset C} \tilde\rho_{ij,C'}(t),$$7 (Ryan et al., 16 Feb 2026). The generalized version replaces the bath-only cluster quantity by a reduced central-spin matrix element computed with the central spin retained in the cluster Hamiltonian (Ryan et al., 16 Feb 2026).

This same recursive structure appears in several applications. In two-qubit gCCE,

$$\tilde\rho_{ij,C}(t) = \rho_{ij,C}(t) \Big/ \prod_{C' \subset C} \tilde\rho_{ij,C'}(t),$$8

with practical truncation

$$\tilde\rho_{ij,C}(t) = \rho_{ij,C}(t) \Big/ \prod_{C' \subset C} \tilde\rho_{ij,C'}(t),$$9

(Chen et al., 2024). In the STIRAP application, each cluster’s irreducible contribution is

$$\rho_{ij}^{(M)}(t) = \prod_{|C|\le M} \tilde\rho_{ij,C}(t).$$0

where $$\rho_{ij}^{(M)}(t) = \prod_{|C|\le M} \tilde\rho_{ij,C}(t).$$1 is the target-state population computed with only bath spins in $$\rho_{ij}^{(M)}(t) = \prod_{|C|\le M} \tilde\rho_{ij,C}(t).$$2 retained (Fazio et al., 2024). For NV longitudinal relaxation, one likewise defines

$$\rho_{ij}^{(M)}(t) = \prod_{|C|\le M} \tilde\rho_{ij,C}(t).$$3

and forms $$\rho_{ij}^{(M)}(t) = \prod_{|C|\le M} \tilde\rho_{ij,C}(t).$$4 from all clusters up to size $$\rho_{ij}^{(M)}(t) = \prod_{|C|\le M} \tilde\rho_{ij,C}(t).$$5 (Yang et al., 2018).

In practice, exact dynamics is recovered only when the bath is taken as one cluster, $$\rho_{ij}^{(M)}(t) = \prod_{|C|\le M} \tilde\rho_{ij,C}(t).$$6, but one usually truncates at $$\rho_{ij}^{(M)}(t) = \prod_{|C|\le M} \tilde\rho_{ij,C}(t).$$7 (Ryan et al., 16 Feb 2026). Computational cost then grows combinatorially with cluster size, so applications rely on convergence with relatively small clusters. In the two-qubit pure-dephasing study, convergence for off-diagonal elements is typically achieved at $$\rho_{ij}^{(M)}(t) = \prod_{|C|\le M} \tilde\rho_{ij,C}(t).$$8–$$\rho_{ij}^{(M)}(t) = \prod_{|C|\le M} \tilde\rho_{ij,C}(t).$$9 over microsecond-and-shorter windows in free induction, and may require $P(t)$0–$P(t)$1 under dynamical decoupling where bath-bath flip-flops dominate (Chen et al., 2024). In the NV longitudinal-relaxation study, 4-CCE and 5-CCE are indistinguishable for the tested $P(t)$2 bath, indicating that 4-CCE suffices there (Yang et al., 2018). In the SiV$P(t)$3 STIRAP study, convergence already occurs at gCCE1 because the laser amplitudes dominate the hyperfine and nuclear-dipolar couplings (Fazio et al., 2024).

A practical implementation perspective is given by PyCCE, which exposes both conventional CCE and gCCE, uses connected clusters under a dipolar-connectivity cutoff, and evaluates small-cluster evolution exactly with matrix exponentials while applying control pulses as explicit unitary rotations (Onizhuk et al., 2021). This indicates that gCCE is best viewed not merely as a formal series, but as a computational framework whose validity depends on both the physical regime and the truncation strategy.

4. Regimes of success: pure dephasing, multi-level coherence, and clock transitions

The most robust domain of gCCE is pure dephasing. In the dipolar-relaxation analysis, the same 9-spin model that shows pathological population dynamics under truncated gCCE exhibits clean, monotonic convergence of $P(t)$4 with increasing cluster order toward the exact result (Ryan et al., 16 Feb 2026). The short-time expansion within a cluster is

$P(t)$5

with $P(t)$6, and the irreducible coefficient obeys

$P(t)$7

(Ryan et al., 16 Feb 2026). Because phase fluctuations from independent clusters add in the exponent, the CCE product is physically correct for dephasing (Ryan et al., 16 Feb 2026).

The two-qubit study extends this logic to a four-level central system. Two electron-spin qubits coupled by isotropic exchange or magnetic dipolar interaction interact with a disordered nuclear bath, and gCCE is used to calculate the decay of off-diagonal elements of the two-qubit reduced density matrix under free evolution and dynamical decoupling (Chen et al., 2024). The formalism captures clusters that couple simultaneously to both electrons, thereby incorporating cross-qubit bath correlations that are absent in single-qubit formulations (Chen et al., 2024). In the exchange-coupled singlet–triplet manifold, the off-diagonal elements separate into three behaviors in the pure-dephasing window: four elements with identical decay envelopes, one element decaying fastest with $P(t)$8 exactly half that of the first group, and one element not decaying to leading order (Chen et al., 2024). These behaviors are traced to effective two-level Hamiltonians derived in the relevant subspaces (Chen et al., 2024).

Near clock transitions, standard CCE may fail because the quadratic dependence of the central splitting on the Overhauser field generates effective long-range interactions $P(t)$9 between all bath-spin pairs (Zhang et al., 2020). The modified/generalized CCE introduced for that regime avoids a perturbative expansion of

$|0\rangle$0

by diagonalizing the central-spin Hamiltonian for each bath product eigenstate $|0\rangle$1, and then defining the conditioned bath Hamiltonians

$|0\rangle$2

(Zhang et al., 2020). The resulting product expansion for fixed $|0\rangle$3,

$|0\rangle$4

absorbs the hyperfine-mediated long-range effects into diagonal, configuration-dependent central energies, leaving the residual interactions as the physical short-range dipolar couplings in $|0\rangle$5 (Zhang et al., 2020). Modified CCE-2 then matches exact dynamics in the examples studied, while the original CCE does not converge near the clock transition (Zhang et al., 2020).

These successes define the part of gCCE that is least controversial: when decoherence is dominated by conditional phase accumulation rather than by real energy transfer, and when cluster contributions behave like additive phase fluctuations in an exponent, the product structure is physically well matched to the problem.

5. Population dynamics, relaxation, and the breakdown of the standard product ansatz

The central controversy surrounding gCCE concerns relaxation driven by spin-spin dipolar interactions. A detailed analysis of a free electron coupled to a bath of eight free electron spins, all dipolarly interacting and initialized with the bath maximally mixed and the central spin in $|0\rangle$6, shows that exact dynamics yields an approximately thermal equilibrium with $|0\rangle$7, whereas truncated gCCE shows no systematic convergence (Ryan et al., 16 Feb 2026). Even orders, CCE-2, 4, and 6, exhibit unphysical dynamics where $|0\rangle$8 exceeds $|0\rangle$9; odd orders, CCE-1, 3, 5, and 7, yield overdamping to $$P^{(M)}(t)=\prod_{|\mathfrak c|\le M}\tilde P_{\mathfrak c}(t),$$0, which does not approach a sensible steady state (Ryan et al., 16 Feb 2026).

The mathematical origin of this failure appears already in the short-time Dyson expansion. For a cluster $$P^{(M)}(t)=\prod_{|\mathfrak c|\le M}\tilde P_{\mathfrak c}(t),$$1,

$$P^{(M)}(t)=\prod_{|\mathfrak c|\le M}\tilde P_{\mathfrak c}(t),$$2

where $$P^{(M)}(t)=\prod_{|\mathfrak c|\le M}\tilde P_{\mathfrak c}(t),$$3 projects onto central-spin states orthogonal to $$P^{(M)}(t)=\prod_{|\mathfrak c|\le M}\tilde P_{\mathfrak c}(t),$$4 (Ryan et al., 16 Feb 2026). Writing

$$P^{(M)}(t)=\prod_{|\mathfrak c|\le M}\tilde P_{\mathfrak c}(t),$$5

the irreducible coefficient satisfies

$$P^{(M)}(t)=\prod_{|\mathfrak c|\le M}\tilde P_{\mathfrak c}(t),$$6

by Möbius inversion on the subset poset (Ryan et al., 16 Feb 2026).

Two generic outcomes then follow. If overlapping pathways make $$P^{(M)}(t)=\prod_{|\mathfrak c|\le M}\tilde P_{\mathfrak c}(t),$$7, one gets $$P^{(M)}(t)=\prod_{|\mathfrak c|\le M}\tilde P_{\mathfrak c}(t),$$8, and multiplying many such factors drives the population outside the physical interval $$P^{(M)}(t)=\prod_{|\mathfrak c|\le M}\tilde P_{\mathfrak c}(t),$$9 (Ryan et al., 16 Feb 2026). If $$P(t)=\tilde P_{\{0\}}\prod_i \tilde P_{\{i\}}\prod_{i,j}\tilde P_{\{i,j\}}\cdots,$$0 for most clusters, then each factor is $$P(t)=\tilde P_{\{0\}}\prod_i \tilde P_{\{i\}}\prod_{i,j}\tilde P_{\{i,j\}}\cdots,$$1, and their product drives the population toward zero, producing overdamped relaxation rather than equilibration with the bath (Ryan et al., 16 Feb 2026). In either case, the problem is structural rather than merely numerical.

The analysis identifies the failure as a mismatch between additive relaxation pathways and multiplicative cluster probabilities. Population transfer pathways are not statistically independent: resonant central-bath flip-flops generate nonlocal, overlapping processes across clusters, particularly in long-range dipolar baths (Ryan et al., 16 Feb 2026). For relaxation, distinct pathways should combine additively at the level of rates or amplitudes, not multiplicatively as independent survival probabilities (Ryan et al., 16 Feb 2026). This is reinforced by a direct test at CCE-2: if one chooses disjoint cluster sets, avoiding overlaps between order-1 and order-2 clusters, the product remains bounded and appears convergent; reintroducing overlapping order-2 clusters with an order-1 cluster immediately reintroduces divergence (Ryan et al., 16 Feb 2026).

Earlier work on NV longitudinal relaxation had already shown that gCCE can numerically describe nearly resonant cross relaxation in a $$P(t)=\tilde P_{\{0\}}\prod_i \tilde P_{\{i\}}\prod_{i,j}\tilde P_{\{i,j\}}\cdots,$$2C bath, with population of $$P(t)=\tilde P_{\{0\}}\prod_i \tilde P_{\{i\}}\prod_{i,j}\tilde P_{\{i,j\}}\cdots,$$3 decaying toward zero on $$P(t)=\tilde P_{\{0\}}\prod_i \tilde P_{\{i\}}\prod_{i,j}\tilde P_{\{i,j\}}\cdots,$$4 time scales near resonance and only small damped oscillations at large detuning (Yang et al., 2018). The later deconstruction does not invalidate those calculations as numerical observations, but it shows that for long-range dipolar relaxation the standard multiplicative population ansatz has no guarantee of qualitative correctness at finite cluster order (Ryan et al., 16 Feb 2026). A plausible implication is that apparent success in specific small or weakly overlapping regimes should be interpreted cautiously unless convergence to the full bath is demonstrated.

The same caution is stated in the two-qubit pure-dephasing study, which remarks on “the application and limitations of gCCE in simulating nuclear-spin induced two-qubit relaxation processes” while emphasizing that its actual simulations remain in the pure-dephasing regime where populations are approximately constant (Chen et al., 2024). Likewise, the STIRAP study applies a population-based gCCE to a regime in which bath couplings are negligible compared with the drive, so that gCCE0 and gCCE1 are nearly indistinguishable (Fazio et al., 2024). This suggests that population-based gCCE can be benign when bath-induced relaxation is perturbatively negligible, but not that the standard product construction is generally reliable for dipolar $$P(t)=\tilde P_{\{0\}}\prod_i \tilde P_{\{i\}}\prod_{i,j}\tilde P_{\{i,j\}}\cdots,$$5 physics.

6. Applications, implementations, and proposed directions for resolution

Published applications of gCCE span several distinct physical settings. For NV centers in diamond, gCCE was introduced to simulate longitudinal relaxation induced by cross relaxation with a $$P(t)=\tilde P_{\{0\}}\prod_i \tilde P_{\{i\}}\prod_{i,j}\tilde P_{\{i,j\}}\cdots,$$6C nuclear bath, thereby extending CCE beyond pure dephasing and into nearly resonant population dynamics (Yang et al., 2018). For two coupled spin qubits, gCCE has been used to compute off-diagonal reduced-density-matrix elements and gate fidelities under free evolution and dynamical decoupling in a disordered nuclear environment, with parameter trends governed by nuclear-free radius, bath density, electron separation, exchange coupling, and magnetic field orientation (Chen et al., 2024). For clock transitions, a modified generalized CCE has been used for nitrogen-vacancy centers near zero magnetic field and for singlet–triplet transitions in double quantum dots, where CCE-2 already converges because the problematic long-range mediated terms are absorbed into bath-configuration-dependent central energies (Zhang et al., 2020). For driven population transfer, gCCE has been applied to STIRAP in the SiV$$P(t)=\tilde P_{\{0\}}\prod_i \tilde P_{\{i\}}\prod_{i,j}\tilde P_{\{i,j\}}\cdots,$$7 defect in 4H-SiC, with a five-state driven system coupled to a spherical nuclear bath of $$P(t)=\tilde P_{\{0\}}\prod_i \tilde P_{\{i\}}\prod_{i,j}\tilde P_{\{i,j\}}\cdots,$$8Si and $$P(t)=\tilde P_{\{0\}}\prod_i \tilde P_{\{i\}}\prod_{i,j}\tilde P_{\{i,j\}}\cdots,$$9C (Fazio et al., 2024).

The SiV$E_M(\delta)$0 STIRAP study is a particularly clear illustration of a regime where generalized population expansion is technically applicable but physically unchallenging. The protocol uses Gaussian pump and Stokes pulses with $E_M(\delta)$1, $E_M(\delta)$2, $E_M(\delta)$3, $E_M(\delta)$4 MHz, and $E_M(\delta)$5 MHz, over a total duration of $E_M(\delta)$6, and finds that for bath radii $E_M(\delta)$7 nm or $E_M(\delta)$8 nm the bath has no appreciable effect on STIRAP, with gCCE converging already at gCCE1 (Fazio et al., 2024). The stated reason is that

$E_M(\delta)$9

so that the coherent optical drive dominates the spin-bath couplings (Fazio et al., 2024).

PyCCE provides a software implementation of both conventional and generalized CCE for realistic spin baths, including bath generation, hyperfine couplings from first principles or point-dipole approximations, cluster generation via connectivity rules, and explicit pulse-controlled cluster evolution (Onizhuk et al., 2021). In PyCCE’s formulation, gCCE is particularly useful near avoided crossings, for non-secular couplings, and for strong driving and time-dependent control sequences, because it evolves the full central-spin-plus-cluster Hamiltonian and applies pulses as explicit rotation operators (Onizhuk et al., 2021). This implementation perspective is consistent with the broader literature: retaining the full central system in the cluster is valuable whenever perturbative projection onto a fixed qubit subspace becomes inaccurate.

The most explicit proposals for fixing the relaxation problem are set out in the dipolar-relaxation analysis. The required change is a reorganization of the expansion so that irreducible contributions to relaxation add rather than multiply (Ryan et al., 16 Feb 2026). Suggested directions include: an additive expansion for a rate kernel $T_2$00 or for transition amplitudes at the level of a $T_2$01-matrix, where cluster-resolved amplitudes are summed before squaring; diagrammatic or linked-cluster resummations that subtract overlaps at the level appropriate for rates rather than survival probabilities; time-dependent projection-operator approaches such as Nakajima–Zwanzig or time-convolutionless master equations with cluster-resolved memory kernels; and partial resummations or renormalizations across strongly overlapping clusters (Ryan et al., 16 Feb 2026). Alternative controlled approaches named in that analysis are tensor-network real-time simulations of spin-bath dynamics and HEOM-like non-Markovian frameworks, which address the additive-rate and memory-kernel structure absent from the standard gCCE product (Ryan et al., 16 Feb 2026).

The present consensus emerging from these works is therefore differentiated rather than uniform. gCCE remains a controlled and convergent tool for central-spin coherences in pure-dephasing regimes, including multi-level and driven settings where central-spin mixing or explicit pulse action matter (Chen et al., 2024, Onizhuk et al., 2021). Modified gCCE constructions can also restore convergence at clock transitions by redefining the unperturbed problem (Zhang et al., 2020). By contrast, for dipolar relaxation driven by overlapping energy-exchange pathways, the standard multiplicative population ansatz fails at the level of principle, and any future relaxation-capable “gCCE” will need a relaxation-specific reorganization of cluster contributions (Ryan et al., 16 Feb 2026).