Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2512.23678

Symbolic recursion method for strongly correlated fermions in two and three dimensions

Igor Ermakov, Oleg Valentinovich Lychkovskiy

Published Dec 29, 2025• 0 citations
arXivHTMLPDF

Abstract

We present a symbolic implementation of the recursion method for dynamical correlations and transport in fermionic systems on one-, two-, and three-dimensional lattices. The implementation is applicable in the strongly correlated regime, yields results directly in the thermodynamic limit, and covers all time scales, from short times through the thermalization time to the late-time asymptotics. Focusing on two paradigmatic models -- interacting spinless fermions and the Hubbard model -- we confirm the universal operator growth hypothesis in the fermionic case, compute infinite-temperature current-current autocorrelation functions, and determine the charge diffusion constant and the high-temperature conductivity. The diffusion constant is obtained in both the perturbative and the nonperturbative regime of the interaction strength, and we locate the boundary between them. We compare our approach with the Majorana propagation method, which converges up to intermediate times and agrees with ours where converged. Our results highlight a symbolic computational paradigm in which the most expensive step is performed once, producing reusable symbolic output that yields physical observables for arbitrary model parameters.

DOI
10.48550/arxiv.2512.23678

Related papers

Quantum dynamics in one and two dimensions via the recursion method

Filipp Gennadevich Uskov, Oleg Valentinovich Lychkovskiy

We report an implementation of the recursion method that addresses quantum many-body dynamics in the nonperturbative regime. The method essentially amounts to constructing a Lanczos basis in the space of operators and solving coupled Heisenberg equations in this basis. The reported implementation has two key ingredients: a computer-algebraic routine for the symbolic calculation of nested commutators and a procedure to extrapolate the sequence of Lanczos coefficients according to the universal operator growth hypothesis. We apply the method to calculate infinite-temperature correlation functions for spin-$1/2$ systems on one- and two-dimensional lattices. In two dimensions the accessible timescale is large enough to essentially embrace the relaxation to equilibrium. The method allows one to accurately calculate transport coefficients. As an illustration, we compute the diffusion constant for the transverse-field Ising model on a square lattice.

arXiv:2401.17211•Physical review. B./Physical review. B•Apr 1, 2024•19 citations

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

Oleg Lychkovskiy, Viacheslav Khrushchev, Ilya Shirokov

We present QCommute, a software tool implemented in C++ for symbolic computation of nested commutators between a Hamiltonian and local observables in quantum many-body spin-1/2 systems on one-, two-, and three-dimensional hypercubic lattices. The computation is performed algebraically directly in the thermodynamic limit, and the Hamiltonian parameters are kept symbolic. Importantly, this way the entire parameter space is covered in a single run. The implementation supports extensive parallelization to achieve high computational performance. QCommute can serve as a computational backend for Heisenberg-picture approaches to quantum dynamics in strongly correlated regimes, ranging from direct Taylor expansion in time to advanced techniques such as the recursion method.

cond-mat.str-elCross-listedarXiv:2604.04778•SciPost Physics Codebases•Apr 7, 2026•1 citation

Nonequilibrium Green Functions Approach to Strongly Correlated Fermions in Lattice Systems

Niclas Schlünzen, M. Bönitz

Abstract Quantum dynamics in strongly correlated systems are of high current interest in many fields including dense plasmas, nuclear matter and condensed matter and ultracold atoms. An important model case are fermions in lattice systems that is well suited to analyze, in detail, a variety of electronic and magnetic properties of strongly correlated solids. Such systems have recently been reproduced with fermionic atoms in optical lattices which allow for a very accurate experimental analysis of the dynamics and of transport processes such as diffusion. The theoretical analysis of such systems far from equilibrium is very challenging since quantum and spin effects as well as correlations have to be treated non‐perturbatively. The only accurate method that has been successful so far are density matrix renormalization group (DMRG) simulations. However, these simulations are presently limited to one‐dimensional (1D) systems and short times. Extension of quantum dynamics simulations to two and three dimensions is commonly viewed as one of the major challenges in this field. Recently we have reported a breakthrough in this area [N. Schlünzen et al., Phys. Rev. B (2016)] where we were able to simulate the expansion dynamics of strongly correlated fermions in a Hubbard lattice following a quench of the confinement potential in 1D, 2D and 3D. The results not only exhibited excellent agreement with the experimental data but, in addition, revealed new features of the short‐time dynamics where correlations and entanglement are being build up. The method used in this work are nonequilibrium Green functions (NEGF) which are found to be very powerful in the treatment of fermionic lattice systems filling the gap presently left open by DMRG in 2D and 3D. In this paper we present a detailed introduction in the NEGF approach and its application to inhomogeneous Hubbard clusters. In detail we discuss the proper strong coupling approximation which is given by T ‐matrix selfenergies that sum up two‐particle scattering processes to infinite order. The efficient numerical implemen‐tation of the method is discussed in detail as it has allowed us to achieve dramatic performance gains. This has been the basis for the treatment of more than 100 particles over large time intervals. The numerical results presented in this paper concentrate on the diffusion in 1D to 3D lattices. We find that the expansion dynamics consist of three different phases that are linked with the build‐up of correlations. In the long time limit, a universal scaling with the particle number is revealed. By extrapolating the expansion velocities to the macroscopic limit, the obtained results show excellent agreement with recent experiments on ultracold fermions in optical lattices. Moreover we present results for the site‐resolved behavior of correlations and entanglement that can be directly compared with experiments using the recently developed atomic microscope technique. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)

Contributions to Plasma Physics•Jan 1, 2016•55 citations

Recursion method for quench dynamics: Strengths and limitations

Ilya Shirokov, Viacheslav Khrushchev, Filipp Gennadevich Uskov +3 more

The recursion method, which solves coupled Heisenberg equations in a Lanczos operator basis, has recently emerged as a powerful nonperturbative tool for computing dynamical correlation functions in strongly correlated two- and three-dimensional quantum many-body systems. Motivated by this success, we investigate whether the method can be extended to expectation values of observables following a quantum quench. We find that such an extension encounters an obstacle absent in the computation of dynamical correlation functions. The latter are fully determined by the Lanczos coefficients b_n b n , which in generic systems exhibit universal behavior, enabling reliable extrapolation from the first few dozens of explicitly computed coefficients. In contrast, quench dynamics additionally requires “quench coefficients” c_n c n , defined as overlaps of Lanczos basis operators with the initial state. We show that, unlike the Lanczos coefficients, the quench coefficients exhibit no universal structure and cannot be reliably extrapolated, thereby limiting the time up to which the method yields accurate results. The behavior of quench coefficients is highly state-dependent, ranging from decaying to irregular or even growing sequences; typically, the less regular the sequence c_n c n , the shorter the accessible timescale. Nevertheless, for favorable initial states, the method remains competitive with state-of-the-art approaches. Moreover, its symbolic implementation allows a single computation to be reused across different Hamiltonian parameters and initial states, making it particularly advantageous in studies requiring extensive scans over Hamiltonian parameters or initial states.

arXiv:2503.24362•SciPost Physics Core•Aug 19, 2026•0 citations

Codebase release 1.0 for QCommute

Oleg Valentinovich Lychkovskiy, Viacheslav Khrushchev, Ilya Shirokov

We present QCommute, a software tool implemented in C++ for symbolic computation of nested commutators between a Hamiltonian and local observables in quantum many-body spin- 1/2 1 / 2 systems on one-, two-, and three-dimensional hypercubic lattices. The computation is performed algebraically directly in the thermodynamic limit, and the Hamiltonian parameters are kept symbolic. Importantly, this way the entire parameter space is covered in a single run. The implementation supports extensive parallelization to achieve high computational performance. QCommute can serve as a computational backend for Heisenberg-picture approaches to quantum dynamics in strongly correlated regimes, ranging from direct Taylor expansion in time to advanced techniques such as the recursion method.

SciPost Physics Codebases•Oct 20, 2026•0 citations

Correlation functions of an interacting spinless fermion model at finite temperature

Kohei Motegi, Kazumitsu Sakai

We formulate correlation functions for a one-dimensional interacting spinless fermion model at finite temperature. By combination of a lattice path integral formulation for thermodynamics with the algebraic Bethe ansatz for fermion systems, the equal-time one-particle Green's function at arbitrary particle density is expressed as a multiple-integral form. Our formula reproduces previously known results in the following three limits: the zero-temperature, the infinite-temperature and the free fermion limits.

arXiv:0712.1399•Journal of Statistical Mechanics Theory and Experiment•Feb 11, 2008•0 citations
We present a symbolic implementation of the recursion method for dynamical correlations and transport in fermionic systems on one-, two-, and three-dimensional lattices. The implementation is applicable in the strongly correlated regime, yields results directly in the thermodynamic limit, and covers all time scales, from short times through the thermalization time to the late-time asymptotics. Focusing on two paradigmatic models -- interacting spinless fermions and the Hubbard model -- we confirm the universal operator growth hypothesis in the fermionic case, compute infinite-temperature current-current autocorrelation functions, and determine the charge diffusion constant and the high-temperature conductivity. The diffusion constant is obtained in both the perturbative and the nonperturbative regime of the interaction strength, and we locate the boundary between them. We compare our approach with the Majorana propagation method, which converges up to intermediate times and agrees with ours where converged. Our results highlight a symbolic computational paradigm in which the most expensive step is performed once, producing reusable symbolic output that yields physical observables for arbitrary model parameters.
We report an implementation of the recursion method that addresses quantum many-body dynamics in the nonperturbative regime. The method essentially amounts to constructing a Lanczos basis in the space of operators and solving coupled Heisenberg equations in this basis. The reported implementation has two key ingredients: a computer-algebraic routine for the symbolic calculation of nested commutators and a procedure to extrapolate the sequence of Lanczos coefficients according to the universal operator growth hypothesis. We apply the method to calculate infinite-temperature correlation functions for spin-1/21/21/2 systems on one- and two-dimensional lattices. In two dimensions the accessible timescale is large enough to essentially embrace the relaxation to equilibrium. The method allows one to accurately calculate transport coefficients. As an illustration, we compute the diffusion constant for the transverse-field Ising model on a square lattice.
We present QCommute, a software tool implemented in C++ for symbolic computation of nested commutators between a Hamiltonian and local observables in quantum many-body spin-1/2 systems on one-, two-, and three-dimensional hypercubic lattices. The computation is performed algebraically directly in the thermodynamic limit, and the Hamiltonian parameters are kept symbolic. Importantly, this way the entire parameter space is covered in a single run. The implementation supports extensive parallelization to achieve high computational performance. QCommute can serve as a computational backend for Heisenberg-picture approaches to quantum dynamics in strongly correlated regimes, ranging from direct Taylor expansion in time to advanced techniques such as the recursion method.
Abstract Quantum dynamics in strongly correlated systems are of high current interest in many fields including dense plasmas, nuclear matter and condensed matter and ultracold atoms. An important model case are fermions in lattice systems that is well suited to analyze, in detail, a variety of electronic and magnetic properties of strongly correlated solids. Such systems have recently been reproduced with fermionic atoms in optical lattices which allow for a very accurate experimental analysis of the dynamics and of transport processes such as diffusion. The theoretical analysis of such systems far from equilibrium is very challenging since quantum and spin effects as well as correlations have to be treated non‐perturbatively. The only accurate method that has been successful so far are density matrix renormalization group (DMRG) simulations. However, these simulations are presently limited to one‐dimensional (1D) systems and short times. Extension of quantum dynamics simulations to two and three dimensions is commonly viewed as one of the major challenges in this field. Recently we have reported a breakthrough in this area [N. Schlünzen et al., Phys. Rev. B (2016)] where we were able to simulate the expansion dynamics of strongly correlated fermions in a Hubbard lattice following a quench of the confinement potential in 1D, 2D and 3D. The results not only exhibited excellent agreement with the experimental data but, in addition, revealed new features of the short‐time dynamics where correlations and entanglement are being build up. The method used in this work are nonequilibrium Green functions (NEGF) which are found to be very powerful in the treatment of fermionic lattice systems filling the gap presently left open by DMRG in 2D and 3D. In this paper we present a detailed introduction in the NEGF approach and its application to inhomogeneous Hubbard clusters. In detail we discuss the proper strong coupling approximation which is given by T ‐matrix selfenergies that sum up two‐particle scattering processes to infinite order. The efficient numerical implemen‐tation of the method is discussed in detail as it has allowed us to achieve dramatic performance gains. This has been the basis for the treatment of more than 100 particles over large time intervals. The numerical results presented in this paper concentrate on the diffusion in 1D to 3D lattices. We find that the expansion dynamics consist of three different phases that are linked with the build‐up of correlations. In the long time limit, a universal scaling with the particle number is revealed. By extrapolating the expansion velocities to the macroscopic limit, the obtained results show excellent agreement with recent experiments on ultracold fermions in optical lattices. Moreover we present results for the site‐resolved behavior of correlations and entanglement that can be directly compared with experiments using the recently developed atomic microscope technique. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
The recursion method, which solves coupled Heisenberg equations in a Lanczos operator basis, has recently emerged as a powerful nonperturbative tool for computing dynamical correlation functions in strongly correlated two- and three-dimensional quantum many-body systems. Motivated by this success, we investigate whether the method can be extended to expectation values of observables following a quantum quench. We find that such an extension encounters an obstacle absent in the computation of dynamical correlation functions. The latter are fully determined by the Lanczos coefficients b_n b n , which in generic systems exhibit universal behavior, enabling reliable extrapolation from the first few dozens of explicitly computed coefficients. In contrast, quench dynamics additionally requires “quench coefficients” c_n c n , defined as overlaps of Lanczos basis operators with the initial state. We show that, unlike the Lanczos coefficients, the quench coefficients exhibit no universal structure and cannot be reliably extrapolated, thereby limiting the time up to which the method yields accurate results. The behavior of quench coefficients is highly state-dependent, ranging from decaying to irregular or even growing sequences; typically, the less regular the sequence c_n c n , the shorter the accessible timescale. Nevertheless, for favorable initial states, the method remains competitive with state-of-the-art approaches. Moreover, its symbolic implementation allows a single computation to be reused across different Hamiltonian parameters and initial states, making it particularly advantageous in studies requiring extensive scans over Hamiltonian parameters or initial states.
We present QCommute, a software tool implemented in C++ for symbolic computation of nested commutators between a Hamiltonian and local observables in quantum many-body spin- 1/2 1 / 2 systems on one-, two-, and three-dimensional hypercubic lattices. The computation is performed algebraically directly in the thermodynamic limit, and the Hamiltonian parameters are kept symbolic. Importantly, this way the entire parameter space is covered in a single run. The implementation supports extensive parallelization to achieve high computational performance. QCommute can serve as a computational backend for Heisenberg-picture approaches to quantum dynamics in strongly correlated regimes, ranging from direct Taylor expansion in time to advanced techniques such as the recursion method.
We formulate correlation functions for a one-dimensional interacting spinless fermion model at finite temperature. By combination of a lattice path integral formulation for thermodynamics with the algebraic Bethe ansatz for fermion systems, the equal-time one-particle Green's function at arbitrary particle density is expressed as a multiple-integral form. Our formula reproduces previously known results in the following three limits: the zero-temperature, the infinite-temperature and the free fermion limits.