QCommute: a tool for symbolic computation of nested commutators in quantum many-body spin-1/2 systems

1. Problem Statement

Quantitatively describing quantum many-body dynamics in non-perturbative regimes is a significant challenge, particularly in two and three dimensions where Matrix Product State (MPS) techniques are less effective. While methods operating in the Heisenberg picture (such as the recursion method and sparse Pauli dynamics) have shown promise, their efficiency relies heavily on the core computational primitive: the evaluation of nested commutators $[H, [H, \dots, [H, A]\dots]]$ between a many-body Hamiltonian ($H$) and an observable ($A$).

Existing tools face several limitations:

• Finite Lattices: Most software operates on finite lattices, introducing boundary effects that require large system sizes to approximate the thermodynamic limit, leading to high memory overhead.
• Numerical Specificity: Many tools compute results for specific numerical values of Hamiltonian parameters, requiring separate runs for different parameter sets.
• Memory Bottlenecks: The number of terms in nested commutators grows exponentially with the nesting order ($n$), often causing out-of-memory (OOM) crashes before reaching high orders necessary for accurate dynamics.

2. Methodology

The authors present QCOMMUTE, a C++17 software tool designed to overcome these limitations through three core methodological pillars:

A. Symbolic Computation in the Thermodynamic Limit

• Translation Invariance: The tool assumes translation-invariant Hamiltonians and observables on infinite hypercubic lattices ($D \in \{1, 2, 3\}$). It utilizes a linear superoperator $T$ to distribute local "densities" (anchored Pauli strings) over the infinite lattice.
• Canonical Form: Operators are represented by their "densities" (finite sums of anchored Pauli strings). The algorithm computes commutators on these densities, collapsing translated terms into single coefficients. This allows the computation to be performed directly in the thermodynamic limit without boundary effects.
• Symbolic Parameters: Hamiltonian parameters (e.g., magnetic fields, coupling constants) are treated as symbolic variables. The coefficients of Pauli strings are computed as polynomials in these parameters with integer coefficients, rather than floating-point numbers. This allows a single run to cover the entire parameter space.

B. Algorithmic Framework

The algorithm iteratively computes the sequence of nested commutators $A^{(n)} = [H, A^{(n-1)}]$:

1. Input: Densities of the Hamiltonian ($h$) and observable ($a$).
2. Commutation: The Commute routine calculates $[H, A^{(n-1)}]$ by iterating over overlapping Pauli strings.
3. Canonicalization & Simplification: Non-anchored strings are translated to the anchor site, and like terms are collected.
4. Arbitrary-Precision Arithmetic: To handle the rapid growth of integer coefficients, the code uses the GMP (GNU Multiple Precision) library, ensuring mathematically exact results regardless of nesting order.

C. Memory Management via Disk Swapping

To address the exponential growth of the number of Pauli strings (which can exceed RAM capacity):

• Chunking: The program does not keep the entire expression in memory.
• Disk Caching: Intermediate results are dynamically split into chunks and written to binary .comm files on the disk (NVMe storage).
• Optimization: This allows the computation of deeply nested commutators (high $n$) on standard hardware by offloading data to disk, avoiding OOM errors.

3. Key Contributions

• First C++ Tool for Infinite Lattices: Unlike Julia or Python packages (e.g., PAULISTRINGS.JL, SPD) that typically operate on finite lattices, QCOMMUTE operates natively in the thermodynamic limit, eliminating boundary effects and reducing memory overhead associated with finite-size scaling.
• Symbolic Parameter Handling: It is one of the few tools (alongside PAULIPROPAGATION.JL) that retains Hamiltonian parameters symbolically. This enables the generation of analytical expressions for dynamical quantities across the entire parameter space in a single execution.
• High-Performance Architecture: The combination of C++17, GMP for exact arithmetic, and optimized disk-swapping strategies allows for the computation of significantly higher nesting orders than previous methods.
• Rigorous Bounds: By providing exact moments ($\mu_{2n}$) and quench coefficients, the tool enables the construction of rigorous upper and lower bounds for dynamical quantities (like autocorrelation functions) using truncated Taylor expansions.

4. Results and Performance

The authors benchmarked QCOMMUTE on the Ising model in 1D, 2D, and 3D:

• Scalability:
  • 1D: Reached nesting order $n_{max} = 48$ on a 1 TB RAM server without disk swapping, generating $\sim 2.4 \times 10^9$ unique Pauli strings.
  • 2D/3D: Successfully handled $n_{max} = 23$ (2D) and $n_{max} = 17$ (3D) by utilizing disk caching. Standard RAM-only approaches would crash around $n \approx 20$ in 2D.
• Comparison: Table 1 in the paper shows QCOMMUTE achieving higher maximal nesting orders than previous state-of-the-art studies (e.g., Parker et al., Noh, De et al.) for the same models.
• Dynamics Application:
  • Quench Dynamics: The tool was used to compute post-quench expectation values via Taylor expansion. Results for the 1D mixed-field Ising model matched Exact Diagonalization (ED) perfectly. In 2D and 3D, results agreed with high-precision Sparse Pauli Dynamics (SPD) and iPEPS data.
  • Autocorrelation: The tool generated tight two-sided bounds for autocorrelation functions. In 1D, these bounds converged to the exact result; in higher dimensions, they provided rigorous constraints where other methods showed discrepancies.

5. Significance

QCOMMUTE represents a significant advancement in the study of quantum many-body dynamics:

1. Benchmarks for New Methods: It provides a high-precision, boundary-free benchmark for testing other approximate methods (like neural networks or tensor networks) in the short-time regime.
2. Non-Perturbative Insights: By enabling high-order expansions in the thermodynamic limit, it allows researchers to probe strongly correlated regimes inaccessible to perturbative approaches.
3. Analytical Power: The symbolic output allows for the derivation of universal properties and the identification of critical points without re-running simulations for every parameter set.
4. Future Directions: The authors note that while currently limited to hypercubic lattices and product states, the framework is a strong foundation for extending to more complex lattices and arbitrary initial states, as well as implementing full Pauli propagation algorithms directly within the package.

In summary, QCOMMUTE bridges the gap between exact analytical methods and numerical simulations, offering a unique capability to compute exact, parameter-independent dynamical quantities for infinite quantum spin systems.