[Paper Review] On the challenge of simulating dipolar contributions to spin relaxation with generalized cluster correlation expansion methods

C. Ryan, Alessandro Lunghi|arXiv (Cornell University)|Feb 16, 2026

Advanced NMR Techniques and Applications0 citations

TL;DR

The paper analyzes the generalized Cluster-Correlation Expansion (gCCE) for spin-relaxation due to dipolar bath interactions and shows it yields unphysical or overdamped results, unlike its performance for dephasing.

ABSTRACT

The study of spin decoherence is often performed by assuming that spin-phonon interactions lead to relaxation at high temperatures, and spin-spin dipolar interactions instead contribute to pure dephasing at low temperatures. This has resulted in the neglect of spin relaxation due to spin-spin dipolar interactions and its influence on decoherence at low temperatures. For a complete understanding of low temperature spin dynamics, it is then imperative to focus also on the latter mechanism. One such method which has shown great promise in the efficient calculation of central spin dynamics due to spin-spin dipolar interactions with a surrounding spin bath is the Cluster-Correlation Expansion (CCE). An extension of this method through the explicit inclusion of the central spin degrees of freedom, known as the generalized Cluster-Correlation Expansion (gCCE) is capable of simulating the transfer of energy from the central spin into the bath, and thus could have the potential to investigate spin relaxation in this setting. In this work, we show that gCCE, in its standard form, is insufficient for providing even a qualitatively accurate description of spin-spin relaxation. A full mathematical deconstruction of the underlying theory of gCCE clearly points to the origin of such a breakdown and provides a starting point for its potential future resolution.

Motivation & Objective

• Motivate a quantitative understanding of spin-spin dipolar relaxation at low temperatures.
• Assess the ability of gCCE to predict central-spin relaxation times (T1) from bath spin dynamics.
• Identify mathematical origins of failures in gCCE when modeling relaxation.
• Compare gCCE performance in relaxation versus pure dephasing regimes.

Proposed method

• Present the central-spin Hamiltonian including dipolar bath interactions (D, A_i, J_ij tensors).
• Formulate the central-spin density matrix evolution and its decomposition into bath-cluster contributions (CCE).
• Derive the irreducible cluster contributions via Möbius inversion to isolate genuine many-body relaxation pathways.
• Expand short-time dynamics to second order to define relaxation rates per cluster (alpha_C).
• Analyze convergence and overlap issues by examining product structure in Eq. (13) and the sign structure of alpha_C^irr.
• Contrast relaxation with dephasing by showing convergent behavior for coherence but not for population relaxation.

Figure 4: CCE Relaxation Convergence. Dynamical evolution of the spin-up state of a free electron interacting with a maximally mixed spin bath of 8 free electrons. The steady state of the dynamics either becomes unphysical or reaches an overdamped value where there is no longer any spin-up populatio

Experimental results

Research questions

• RQ1 Can gCCE reliably predict T1-like relaxation times from dipolar spin baths?
• RQ2 What mathematical factors cause unphysical or overdamped relaxation dynamics in gCCE?
• RQ3 Why does gCCE appear well-behaved for dephasing but fail for relaxation?

Key findings

• gCCE yields either unphysical populations (outside [0,1]) or overdamped relaxation in central-spin T1 calculations.
• The root cause is overlapping relaxation pathways among bath-spin clusters, leading to breakdown of the simple product form.
• Möbius inversion reveals irreducible cluster contributions can be negative, triggering nonphysical behavior in the central-spin populations.
• When clusters have strong overlaps, the product of cluster contributions can diverge or converge to an unrealistic steady state with no initial population.
• In the pure dephasing regime, gCCE does converge toward the correct result, indicating regime-dependent validity.

Figure 5: CCE Convergence/Divergence. Simulated dynamics for a 9 electron spin system at CCE order 2. Convergent behavior is calculated when there is no overlap between the order 2 clusters and the order 1 clusters (blue line), while the dynamics diverges once the order 2 clusters begin to overlap w