We describe an efficient numerical method for simulating the dynamics of interacting spin ensembles in the presence of dephasing and decay. The method builds on the discrete truncated Wigner approximation for isolated systems, which combines the mean-field dynamics of a spin ensemble with a Monte Carlo sampling of discrete initial spin values to account for quantum correlations. Here we show how this approach can be generalized for dissipative spin systems by replacing the deterministic mean-field evolution by a stochastic process, which describes the decay of coherences and populations while preserving the length of each spin. We demonstrate the application of this technique for simulating nonclassical spin-squeezing effects or the dynamics and steady states of cavity QED models with ${10}^{5}$ interacting two-level systems. This opens up the possibility to perform accurate real-scale simulations of a diverse range of experiments in quantum optics or with solid-state spin ensembles under realistic laboratory conditions.We describe an efficient numerical method for simulating the dynamics of interacting spin ensembles in the presence of dephasing and decay. The method builds on the discrete truncated Wigner approximation for isolated systems, which combines the mean-field dynamics of a spin ensemble with a Monte Carlo sampling of discrete initial spin values to account for quantum correlations. Here we show how this approach can be generalized for dissipative spin systems by replacing the deterministic mean-field evolution by a stochastic process, which describes the decay of coherences and populations while preserving the length of each spin. We demonstrate the application of this technique for simulating nonclassical spin-squeezing effects or the dynamics and steady states of cavity QED models with ${10}^{5}$ interacting two-level systems. This opens up the possibility to perform accurate real-scale simulations of a diverse range of experiments in quantum optics or with solid-state spin ensembles under realistic laboratory conditions.

arXiv:2105.00004Physical Review AJan 24, 202234 citations

An efficient method of spin dynamics for long-range interacting systems

Taichi Hinokihara, M. Sasaki, Seiji Miyashita

Bulletin of the American Physical SocietyMar 2, 20200 citations

First-principles simulation of spin diffusion in static solids using dynamic mean-field theory

Timo Gräßer, Götz S. Uhrig, Matthias Ernst

The dynamics of disordered nuclear spin ensembles are the subject of nuclear magnetic resonance studies. Due to the through-space long-range dipolar interaction generically many spins are involved in the time evolution, so that exact brute force calculations are impossible. The recently established spin dynamic mean-field theory (spinDMFT) represents an efficient and unbiased alternative to overcome this challenge. The approach only requires the dipolar couplings as input and the only prerequisite for its applicability is that each spin interacts with a large number of other spins. In this article, we show that spinDMFT can be used to describe spectral spin diffusion in static samples and to simulate zero-quantum line shapes which eluded an efficient quantitative simulation so far to the best of our knowledge. We perform benchmarks for two test substances that establish an excellent match with published experimental data. As spinDMFT combines low computational effort with high accuracy, we strongly suggest to use it for large-scale simulations of spin diffusion, which are important in various areas of magnetic resonance.The dynamics of disordered nuclear spin ensembles are the subject of nuclear magnetic resonance studies. Due to the through-space long-range dipolar interaction generically many spins are involved in the time evolution, so that exact brute force calculations are impossible. The recently established spin dynamic mean-field theory (spinDMFT) represents an efficient and unbiased alternative to overcome this challenge. The approach only requires the dipolar couplings as input and the only prerequisite for its applicability is that each spin interacts with a large number of other spins. In this article, we show that spinDMFT can be used to describe spectral spin diffusion in static samples and to simulate zero-quantum line shapes which eluded an efficient quantitative simulation so far to the best of our knowledge. We perform benchmarks for two test substances that establish an excellent match with published experimental data. As spinDMFT combines low computational effort with high accuracy, we strongly suggest to use it for large-scale simulations of spin diffusion, which are important in various areas of magnetic resonance.

cond-mat.mtrl-sciCross-listedarXiv:2512.15572Dec 18, 2025

Monte Carlo Simulation of Spin Models with Long-Range Interactions

Erik Luijten

arXiv:cond-mat/9906219Springer proceedings in physicsJan 1, 20001 citation

A generalized phase space approach for solving quantum spin dynamics

Bihui Zhu, Ana Maria Rey, Johannes Schachenmayer

Abstract Numerical techniques to efficiently model out-of-equilibrium dynamics in interacting quantum many-body systems are key for advancing our capability to harness and understand complex quantum matter. Here we propose a new numerical approach which we refer to as generalized discrete truncated Wigner approximation (GDTWA). It is based on a discrete semi-classical phase space sampling and allows to investigate quantum dynamics in lattice spin systems with arbitrary S ≥ 1/2. We show that the GDTWA can accurately simulate dynamics of large ensembles in arbitrary dimensions. We apply it for S > 1/2 spin-models with dipolar long-range interactions, a scenario arising in recent experiments with magnetic atoms. We show that the method can capture beyond mean-field effects, not only at short times, but it also can correctly reproduce long time quantum-thermalization dynamics. We benchmark the method with exact diagonalization in small systems, with perturbation theory for short times, and with analytical predictions made for models which feature quantum-thermalization at long times. We apply our method to study dynamics in large S > 1/2 spin-models and compute experimentally accessible observables such as Zeeman level populations, contrast of spin coherence, spin squeezing, and entanglement quantified by single-spin Renyi entropies. We reveal that large S systems can feature larger entanglement than corresponding S = 1/2 systems. Our analyses demonstrate that the GDTWA can be a powerful tool for modeling complex spin dynamics in regimes where other state-of-the art numerical methods fail.Abstract Numerical techniques to efficiently model out-of-equilibrium dynamics in interacting quantum many-body systems are key for advancing our capability to harness and understand complex quantum matter. Here we propose a new numerical approach which we refer to as generalized discrete truncated Wigner approximation (GDTWA). It is based on a discrete semi-classical phase space sampling and allows to investigate quantum dynamics in lattice spin systems with arbitrary S ≥ 1/2. We show that the GDTWA can accurately simulate dynamics of large ensembles in arbitrary dimensions. We apply it for S > 1/2 spin-models with dipolar long-range interactions, a scenario arising in recent experiments with magnetic atoms. We show that the method can capture beyond mean-field effects, not only at short times, but it also can correctly reproduce long time quantum-thermalization dynamics. We benchmark the method with exact diagonalization in small systems, with perturbation theory for short times, and with analytical predictions made for models which feature quantum-thermalization at long times. We apply our method to study dynamics in large S > 1/2 spin-models and compute experimentally accessible observables such as Zeeman level populations, contrast of spin coherence, spin squeezing, and entanglement quantified by single-spin Renyi entropies. We reveal that large S systems can feature larger entanglement than corresponding S = 1/2 systems. Our analyses demonstrate that the GDTWA can be a powerful tool for modeling complex spin dynamics in regimes where other state-of-the art numerical methods fail.

arXiv:1905.08782New Journal of PhysicsJul 25, 201959 citations

{10}^{5}

interacting two-level systems. This opens up the possibility to perform accurate real-scale simulations of a diverse range of experiments in quantum optics or with solid-state spin ensembles under realistic laboratory conditions.