English

LREI: A fast numerical solver for quantum Landau-Lifshitz equations

Quantum Physics 2025-09-01 v1 Mesoscale and Nanoscale Physics Materials Science Numerical Analysis Numerical Analysis Computational Physics

Abstract

We develop LREI (Low-Rank Eigenmode Integration), a memory- and time-efficient scheme for solving quantum Landau-Lifshitz (q-LL) and quantum Landau-Lifshitz-Gilbert (q-LLG) equations, which govern spin dynamics in open quantum systems. Although system size grows exponentially with the number of spins, our approach exploits the low-rank structure of the density matrix and the sparsity of Hamiltonians to avoid full matrix computations. By representing density matrices via low-rank factors and applying Krylov subspace methods for partial eigendecompositions, we reduce the per-step complexity of Runge-Kutta and Adams-Bashforth schemes from O(N3)\mathcal{O}(N^3) to O(r2N)\mathcal{O}(r^2N), where N=2nN = 2^n is the Hilbert space dimension for nn spins and rNr \ll N the effective rank. Similarly, memory costs shrink from O(N2)\mathcal{O}(N^2) to O(rN)\mathcal{O}(rN), since no full N×NN\times N matrices are formed. A key advance is handling the invariant subspace of zero eigenvalues. By using Householder reflectors built for the dominant eigenspace, we perform the solution entirely without large matrices. For example, a time step of a twenty-spin system, with density matrix size over one million, now takes only seconds on a standard laptop. Both Runge-Kutta and Adams-Bashforth methods are reformulated to preserve physical properties of the density matrix throughout evolution. This low-rank algorithm enables simulations of much larger spin systems, which were previously infeasible, providing a powerful tool for comparing q-LL and q-LLG dynamics, testing each model validity, and probing how quantum features such as correlations and entanglement evolve across different regimes of system size and damping.

Keywords

Cite

@article{arxiv.2508.21200,
  title  = {LREI: A fast numerical solver for quantum Landau-Lifshitz equations},
  author = {Davoud Mirzaei and Behnam Hashemi and Vahid Azimi-Mousolou},
  journal= {arXiv preprint arXiv:2508.21200},
  year   = {2025}
}

Comments

23 pages and 12 figures