Symplectic split-operator method for the time-dependent unitary Tavis-Cummings model

Andrii G. Sotnikov; Denys I. Bondar; Kurt Jacobs; Roman Ovsiannikov

arxiv: 2604.21778 · v2 · submitted 2026-04-23 · 🪐 quant-ph

Symplectic split-operator method for the time-dependent unitary Tavis-Cummings model

Roman Ovsiannikov , Kurt Jacobs , Andrii G. Sotnikov , Denys I. Bondar This is my paper

Pith reviewed 2026-05-09 22:42 UTC · model grok-4.3

classification 🪐 quant-ph

keywords methodmodelbasiscavitymodesystemtavis-cummingstri-diagonal

0 comments

The pith

A fast, memory-efficient symplectic split-operator method for the time-dependent Tavis-Cummings model that achieves linear scaling by re-indexing the Hamiltonian to tri-diagonal form.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The Tavis-Cummings model describes how multiple spins or atoms interact with a single mode of light in a cavity. When the spin frequencies or other parameters change over time, exact simulation becomes hard because the quantum state evolves according to a complex Hamiltonian. The new approach notices that a simple reordering of the mathematical basis states turns the interaction terms into a tri-diagonal structure. This allows a split-operator technique, which splits the evolution into easy-to-compute pieces, to run efficiently while keeping the quantum evolution unitary. Truncating the number of light quanta considered keeps memory and time costs linear in the total number of states. The same trick applies to any closed quantum system whose Hamiltonian can be made tri-diagonal by reordering basis elements.

Core claim

The method exploits the fact that, while the Tavis-Cummings model is not tri-diagonal, it can be brought into tri-diagonal form by a change of basis that can be implemented purely by re-indexing (permuting basis elements), which is a fast operation. By truncating the Fock basis of the cavity mode, the computational complexity of the method is linear in the total dimension of the coupled system, both in time and memory.

Load-bearing premise

That the Tavis-Cummings Hamiltonian (and similar systems) can always be made exactly tri-diagonal by a pure permutation of basis elements without additional transformations or approximations that would affect accuracy or unitarity.

Figures

Figures reproduced from arXiv: 2604.21778 by Andrii G. Sotnikov, Denys I. Bondar, Kurt Jacobs, Roman Ovsiannikov.

read the original abstract

We present a fast, memory-efficient, unitarity-preserving numerical method beyond the rotating-wave approximation for the closed Tavis-Cummings model in which a multilevel spin system interacts with a cavity mode. This model can describe the interaction of an ensemble of spins with a cavity mode in which the spin frequency and other parameters are time-dependent. The method exploits the fact that, while the Tavis-Cummings model is not tri-diagonal, it can be brought into tri-diagonal form by a change of basis that can be implemented purely by re-indexing (permuting basis elements), which is a fast operation. By truncating the Fock basis of the cavity mode, the computational complexity of the method is linear in the total dimension of the coupled system, both in time and memory. The method can be employed to simulate any closed quantum system whose Hamiltonian terms can be brought into tri-diagonal form.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The re-indexing trick for tridiagonalizing the Tavis-Cummings Hamiltonian doesn't work for N > 1, undermining the linear scaling claim.

read the letter

Your colleague should know right away that this paper's proposed symplectic split-operator method for the time-dependent Tavis-Cummings model hinges on re-indexing the basis to get a tridiagonal matrix, but that re-indexing is not feasible once you have more than one spin in the system. The new part is the combination of the split-operator approach with handling time-dependent parameters in the full unitary model, without relying on the rotating wave approximation. They emphasize that the re-indexing is just a permutation, which would be cheap, and then the integration becomes linear in the dimension after truncating the cavity Fock states. That would be a nice practical tool if it panned out, especially for larger ensembles in cavity QED. What works in the paper is the motivation: simulating these systems accurately with time-varying frequencies or couplings is relevant for experiments, and preserving unitarity in the numerical scheme is a solid goal. The abstract lays out the complexity claims clearly. The soft spot is the assumption that any Tavis-Cummings Hamiltonian can be made tri-diagonal by permutation alone. As the stress-test note explains, the interaction graph for N spins has maximum degree 2N, which exceeds 2 for N greater than 1. A tridiagonal matrix after permutation similarity requires the graph to be a collection of paths. This means the method as described can't deliver exact linear scaling without some other approximation that isn't mentioned. The paper doesn't provide the derivation or benchmarks that would let us check if they have a workaround or if the claim is limited to N=1. This kind of paper would interest researchers who need efficient solvers for quantum optical Hamiltonians in control or simulation contexts. A reader already familiar with split-operator methods might pick up the re-indexing idea, but they'd have to verify the tri-diagonal claim themselves. I wouldn't recommend sending this to peer review yet. The core technical step needs to be fixed or clarified first, otherwise referees will spend time on the same graph theory point.

Referee Report

2 major / 2 minor

Summary. The paper proposes a symplectic split-operator numerical method for the time-dependent closed Tavis-Cummings model (N spins coupled to a cavity mode, beyond the rotating-wave approximation). It claims that a change of basis implemented purely by re-indexing (permuting basis elements) renders the Hamiltonian tri-diagonal, enabling a linear-complexity algorithm in both time and memory (after Fock-space truncation) that exactly preserves unitarity. The method is asserted to apply to any closed quantum system whose Hamiltonian terms can be made tri-diagonal.

Significance. If the central algorithmic claim held, the approach would offer a scalable, structure-preserving integrator for large-N spin-cavity systems with time-dependent parameters, filling a gap between exact diagonalization (exponential cost) and approximate methods such as mean-field or truncated Wigner. The emphasis on symplectic/unitary preservation and the absence of free parameters in the core integrator would be genuine strengths for long-time coherent simulations.

major comments (2)

[Abstract, §2] Abstract and §2 (method description): the assertion that the Tavis-Cummings Hamiltonian 'can be brought into tri-diagonal form by a change of basis that can be implemented purely by re-indexing (permuting basis elements)' is incorrect for N ≥ 2. In the product basis |{m_i}, n⟩ each state couples to up to 2N distinct states via the rotating and counter-rotating terms; the underlying undirected graph therefore has maximum degree 2N > 2. A matrix is permutation-similar to a tridiagonal matrix only if its graph is a disjoint union of paths (maximum degree ≤ 2). No such permutation exists, so the claimed O(D) complexity (D = total dimension) cannot be achieved by re-indexing alone without additional approximations that would violate the stated exact unitarity preservation.
[§3, §4] §3 (split-operator construction) and §4 (numerical examples): the linear-complexity claim and the unitarity proof rest on the tridiagonal property obtained in the preceding step. Because that property does not hold, the subsequent symplectic splitting and error analysis do not apply to the Tavis-Cummings model as stated; any numerical demonstrations for N > 1 must therefore rely on an unstated approximation or truncation whose effect on unitarity and accuracy is not quantified.

minor comments (2)

The manuscript supplies no derivation of the claimed re-indexing permutation, no explicit matrix elements after re-ordering, and no operation-count table confirming O(D) scaling; these omissions make it impossible to verify the complexity statements even if the graph-theoretic objection were set aside.
Notation for the time-dependent parameters (ω(t), g(t), …) is introduced without a clear table or equation reference, complicating reproduction of the numerical tests.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the domain assumption that the Tavis-Cummings Hamiltonian can be exactly tri-diagonalized by basis permutation alone and that Fock-space truncation preserves the essential dynamics for the intended simulations.

axioms (1)

domain assumption The Tavis-Cummings model Hamiltonian can be brought into tri-diagonal form by a change of basis implemented purely by re-indexing (permuting basis elements).
This is the key exploitation stated in the abstract that enables the linear scaling and split-operator approach.

pith-pipeline@v0.9.0 · 5464 in / 1367 out tokens · 25422 ms · 2026-05-09T22:42:26.171538+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

52 extracted references · 52 canonical work pages · 1 internal anchor

[1]

J. R. Johansson, P. D. Nation, and F. Nori, QuTiP: An open- source Python framework for the dynamics of open quantum systems, Computer Physics Communications183, 1760 (2012)

work page 2012
[2]

J. R. Johansson, P. D. Nation, and F. Nori, QuTiP 2: A Python framework for the dynamics of open quantum systems, Com- puter Physics Communications184, 1234 (2013)

work page 2013
[3]

Blanes, F

S. Blanes, F. Casas, and A. Murua, Splitting methods for differ- ential equations, Acta Numerica33, 1 (2024)

work page 2024
[4]

Schneider and H

R. Schneider and H. Gharibnejad, Numerical methods for the time-dependent schrödinger equation: Beyond short-time prop- agators, Atoms13, 70 (2025)

work page 2025
[5]

H. F. Trotter, On the product of semi-groups of operators, Pro- ceedings of the American Mathematical Society10, 545 (1959)

work page 1959
[6]

M. Suzuki, Generalized Trotter’s formula and systematic ap- proximants of exponential operators and inner derivations with applications to many-body problems, Communications in Mathematical Physics51, 183 (1976)

work page 1976
[7]

Strang, On the construction and comparison of difference schemes, SIAM Journal on Numerical Analysis5, 506 (1968)

G. Strang, On the construction and comparison of difference schemes, SIAM Journal on Numerical Analysis5, 506 (1968)

work page 1968
[8]

M. D. Feit, J. A. Fleck, and A. D. Steiger, Solution of the Schrödinger equation by a spectral method, Journal of Com- putational Physics47, 412 (1982)

work page 1982
[9]

S. K. Gray and D. E. Manolopoulos, Symplectic integrators tai- lored to the time-dependent Schrödinger equation, The Journal of Chemical Physics104, 7099 (1996)

work page 1996
[10]

D. I. Bondar and O. Steuernagel, Symplectic split-operator propagators from tridiagonalized multi-mode bosonic hilbert spaces for bose-hubbard hamiltonians (2026), arXiv:2603.25639 [quant-ph]

work page arXiv 2026
[11]

E. T. Jaynes and F. W. Cummings, Comparison of quantum and semiclassical radiation theories with application to the beam maser, Proceedings of the IEEE51, 89 (1963)

work page 1963
[12]

Tavis and F

M. Tavis and F. W. Cummings, Exact solution for an N- molecule-radiation-field Hamiltonian, Physical Review170, 379 (1968)

work page 1968
[13]

Tavis and F

M. Tavis and F. W. Cummings, Approximate solutions for an N-molecule-radiation-field Hamiltonian, Physical Review188, 692 (1969)

work page 1969
[14]

B. W. Shore and P. L. Knight, The Jaynes–Cummings model, Journal of Modern Optics40, 1195 (1993)

work page 1993
[15]

M. W. Doherty, N. B. Manson, P. Delaney, F. Jelezko, J. Wrachtrup, and L. C. L. Hollenberg, The nitrogen-vacancy colour centre in diamond, Physics Reports528, 1 (2013)

work page 2013
[16]

Amsüss, C

R. Amsüss, C. Koller, T. Nöbauer, S. Putz, S. Rotter, K. Sand- ner, S. Schneider, M. Schramböck, G. Steinhauser, H. Ritsch, J. Schmiedmayer, and J. Majer, Cavity QED with magnetically coupled collective spin states, Physical Review Letters107, 060502 (2011)

work page 2011
[17]

Y . Kubo, F. R. Ong, P. Bertet, D. Vion, V . Jacques, D. Zheng, A. Dréau, J.-F. Roch, A. Auffeves, F. Jelezko, J. Wrachtrup, M. F. Barthe, P. Bergonzo, and D. Esteve, Strong coupling of a spin ensemble to a superconducting resonator, Physical Review Letters105, 140502 (2010)

work page 2010
[18]

J. D. Breeze, E. Salvadori, J. Sathian, N. M. Alford, and C. W. M. Kay, Room-temperature cavity quantum electrody- namics with strongly coupled Dicke states, npj Quantum Infor- mation3, 40 (2017)

work page 2017
[19]

Imamo ˘glu, Cavity QED based on collective magnetic dipole coupling: Spin ensembles as hybrid two-level systems, Phys

A. Imamo ˘glu, Cavity QED based on collective magnetic dipole coupling: Spin ensembles as hybrid two-level systems, Phys. Rev. Lett.102, 083602 (2009)

work page 2009
[20]

D. I. Schuster, A. P. Sears, E. Ginossar, L. DiCarlo, L. Frunzio, J. J. L. Morton, H. Wu, G. A. D. Briggs, B. B. Buckley, D. D. Awschalom, and R. J. Schoelkopf, High-cooperativity coupling of electron-spin ensembles to superconducting cavities, Phys. Rev. Lett.105, 140501 (2010)

work page 2010
[21]

E. R. Eisenach, J. F. Barry, M. F. O’Keeffe,et al., Cavity- enhanced microwave readout of a solid-state spin sensor, Na- ture Communications12, 1357 (2021)

work page 2021
[22]

Zhang, Q

Y . Zhang, Q. Wu, S.-L. Su, Q. Lou, C. Shan, and K. Mølmer, Cavity quantum electrodynamics effects with nitrogen vacancy center spins coupled to room temperature microwave res- onators, Physical Review Letters128, 253601 (2022)

work page 2022
[23]

W. Ng, H. Wu, and M. Oxborrow, Quasi-continuous cooling of a microwave mode on a benchtop using hyperpolarized NV- diamond, Applied Physics Letters119, 234001 (2021)

work page 2021
[24]

D. P. Fahey, K. Jacobs, M. J. Turner, H. Choi, J. E. Hoffman, D. Englund, and M. E. Trusheim, Steady-state microwave mode 6 cooling with a diamond n-vensemble, Phys. Rev. Appl.20, 014033 (2023)

work page 2023
[25]

Zhang, Q

Y . Zhang, Q. Wu, H. Wu, X. Yang, S.-L. Su, C. Shan, and K. Mølmer, Microwave mode cooling and cavity quantum elec- trodynamics effects at room temperature with optically cooled nitrogen-vacancy center spins, npj Quantum Information8, 125 (2022)

work page 2022
[26]

Q. Wu, Y . Zhang, X. Yang,et al., A superradiant maser with nitrogen-vacancy center spins, Science China Physics, Mechan- ics & Astronomy65, 217311 (2022)

work page 2022
[27]

Kersten, N

W. Kersten, N. de Zordo, O. Diekmann, E. S. Redchenko, et al., Self-induced superradiant masing, Nature Physics22, 158 (2026)

work page 2026
[28]

H. Wang, K. L. Tiwari, K. Jacobs, M. Judy, X. Zhang, D. R. Englund, and M. E. Trusheim, A spin-refrigerated cavity quan- tum electrodynamic sensor, Nature Communications15, 10320 (2024)

work page 2024
[29]

J. C. Hermann, R. Rizzato, F. Bruckmaier, R. D. Allert, A. Blank, and D. B. Bucher, Extending radiowave frequency detection range with dressed states of solid-state spin ensem- bles, npj Quantum Information10, 103 (2024)

work page 2024
[30]

Z. Yin, J. J. Welter, C. A. Hart, P. V . Petruzzi,et al., High- resolution and wide-frequency-range magnetic spectroscopy with solid-state spin ensembles, npj Quantum Information11, 190 (2025)

work page 2025
[31]

H. Wang, S. Wu, K. Jacobs, Y . S. Duan, D. R. Englund, and M. E. Trusheim, Cavity-enhanced solid-state nuclear spin gy- roscope, Physical Review Letters134, 183603 (2025)

work page 2025
[32]

H. Wang, K. Jacobs, D. R. Englund, and M. E. Trusheim, Ultralow-power microwave frequency comb at a bistable phase transition, Phys. Rev. X16, 021005 (2026)

work page 2026
[33]

H. Wang, K. Jacobs, D. Fahey, Y . Hu, D. R. Englund, and M. E. Trusheim, Exceptional sensitivity near the bistable tran- sition point of a hybrid quantum system, Nature Physics 10.1038/s41567-026-03217-3 (2026)

work page doi:10.1038/s41567-026-03217-3 2026
[34]

Ovsiannikov, K

R. Ovsiannikov, K. Jacobs, A. G. Sotnikov, M. E. Trusheim, and D. I. Bondar, Hybrid non-degenerate parametric amplifier for a microwave cavity mode and an NV ensemble (2026), arXiv:2601.03407 [quant-ph]

work page arXiv 2026
[35]

NV-ensemble enabled microwave/NV parametric amplifier with optimal driving

R. Ovsiannikov, K. Jacobs, A. G. Sotnikov, M. E. Trusheim, and D. I. Bondar, NV-ensemble enabled microwave/NV para- metric amplifier with optimal driving (2026), arXiv:2604.11643 [quant-ph]

work page internal anchor Pith review Pith/arXiv arXiv 2026
[36]

Holstein and H

T. Holstein and H. Primakoff, Field dependence of the intrinsic domain magnetization of a ferromagnet, Physical Review58, 1098 (1940)

work page 1940
[37]

L. G. Gunderman, T. W. Borneman, and D. G. Cory, Thermal state structure in the Tavis–Cummings model and rapid simu- lations in mesoscopic quantum ensembles, Phys. Rev. Research 7, 023155 (2025)

work page 2025
[38]

D. I. Bondar, G. McCaul, and A. Sotnikov, Observables in Mo- tion: A guide to simulating classical and quantum dynamics (2025)

work page 2025
[39]

Lambert, E

N. Lambert, E. Giguère, P. Menczel, B. Li, P. Hopf, G. Suarez, M. Gali Labarias, J. Lishman, R. Gadhvi, R. Agarwal, A. Gali- cia, N. Shammah, P. Nation, J. Johansson, S. Ahmed, S. Cross, A. Pitchford, and F. Nori, Qutip 5: The quantum toolbox in python, Physics Reports1153, 1 (2026)

work page 2026
[40]

R. H. Dicke, Coherence in spontaneous radiation processes, Physical Review93, 99 (1954)

work page 1954
[41]

For example, in Python,scipy.linalg.eigh_tridiagonal

work page
[42]

Crank and P

J. Crank and P. Nicolson, A practical method for numerical evaluation of solutions of partial differential equations of the heat-conduction type, Mathematical Proceedings of the Cam- bridge Philosophical Society43, 50 (1947)

work page 1947
[43]

Goldberg, H

A. Goldberg, H. M. Schey, and J. L. Schwartz, Computer- generated motion pictures of one-dimensional quantum- mechanical transmission and reflection phenomena, American Journal of Physics35, 177 (1967)

work page 1967
[44]

L. H. Thomas,Elliptic Problems in Linear Difference Equations over a Network, Tech. Rep. (Watson Scientific Computing Lab- oratory, Columbia University, New York, 1949)

work page 1949
[45]

Le Bris and P

C. Le Bris and P. Rouchon, Low-rank numerical approxima- tions for high-dimensional Lindblad equations, Physical Re- view A87, 022125 (2013)

work page 2013
[46]

Finazzi, A

S. Finazzi, A. Le Boité, F. Storme, A. Baksic, and C. Ciuti, Corner-Space Renormalization Method for Driven-Dissipative Two-Dimensional Correlated Systems, Physical Review Letters 115, 080604 (2015)

work page 2015
[47]

McCaul, K

G. McCaul, K. Jacobs, and D. I. Bondar, Fast computation of dissipative quantum systems with ensemble rank truncation, Physical Review Research3, 013017 (2021)

work page 2021
[48]

Y .-T. Chen, C. Farquhar, and R. M. Parrish, Low-rank density- matrix evolution for noisy quantum circuits, npj Quantum In- formation7, 61 (2021)

work page 2021
[49]

Donatella, Z

K. Donatella, Z. Denis, A. Le Boité, and C. Ciuti, Continuous- time dynamics and error scaling of noisy highly entangling quantum circuits, Physical Review A104, 062407 (2021)

work page 2021
[50]

Suzuki, Fractal decomposition of exponential operators with applications to many-body theories and Monte Carlo simula- tions, Phys

M. Suzuki, Fractal decomposition of exponential operators with applications to many-body theories and Monte Carlo simula- tions, Phys. Lett. A146, 319 (2002)

work page 2002
[51]

A. D. Bandrauk and H. Shen, Higher order exponential split operator method for solving time-dependent Schrödinger equa- tions, Canadian Journal of Chemistry70, 555 (1992)

work page 1992
[52]

Ovsiannikov,https : / / github

R. Ovsiannikov,https : / / github . com / dibondar / NewPropagators(2026), gitHub repository, accessed 2026-04-22

work page 2026

[1] [1]

J. R. Johansson, P. D. Nation, and F. Nori, QuTiP: An open- source Python framework for the dynamics of open quantum systems, Computer Physics Communications183, 1760 (2012)

work page 2012

[2] [2]

J. R. Johansson, P. D. Nation, and F. Nori, QuTiP 2: A Python framework for the dynamics of open quantum systems, Com- puter Physics Communications184, 1234 (2013)

work page 2013

[3] [3]

Blanes, F

S. Blanes, F. Casas, and A. Murua, Splitting methods for differ- ential equations, Acta Numerica33, 1 (2024)

work page 2024

[4] [4]

Schneider and H

R. Schneider and H. Gharibnejad, Numerical methods for the time-dependent schrödinger equation: Beyond short-time prop- agators, Atoms13, 70 (2025)

work page 2025

[5] [5]

H. F. Trotter, On the product of semi-groups of operators, Pro- ceedings of the American Mathematical Society10, 545 (1959)

work page 1959

[6] [6]

M. Suzuki, Generalized Trotter’s formula and systematic ap- proximants of exponential operators and inner derivations with applications to many-body problems, Communications in Mathematical Physics51, 183 (1976)

work page 1976

[7] [7]

Strang, On the construction and comparison of difference schemes, SIAM Journal on Numerical Analysis5, 506 (1968)

G. Strang, On the construction and comparison of difference schemes, SIAM Journal on Numerical Analysis5, 506 (1968)

work page 1968

[8] [8]

M. D. Feit, J. A. Fleck, and A. D. Steiger, Solution of the Schrödinger equation by a spectral method, Journal of Com- putational Physics47, 412 (1982)

work page 1982

[9] [9]

S. K. Gray and D. E. Manolopoulos, Symplectic integrators tai- lored to the time-dependent Schrödinger equation, The Journal of Chemical Physics104, 7099 (1996)

work page 1996

[10] [10]

D. I. Bondar and O. Steuernagel, Symplectic split-operator propagators from tridiagonalized multi-mode bosonic hilbert spaces for bose-hubbard hamiltonians (2026), arXiv:2603.25639 [quant-ph]

work page arXiv 2026

[11] [11]

E. T. Jaynes and F. W. Cummings, Comparison of quantum and semiclassical radiation theories with application to the beam maser, Proceedings of the IEEE51, 89 (1963)

work page 1963

[12] [12]

Tavis and F

M. Tavis and F. W. Cummings, Exact solution for an N- molecule-radiation-field Hamiltonian, Physical Review170, 379 (1968)

work page 1968

[13] [13]

Tavis and F

M. Tavis and F. W. Cummings, Approximate solutions for an N-molecule-radiation-field Hamiltonian, Physical Review188, 692 (1969)

work page 1969

[14] [14]

B. W. Shore and P. L. Knight, The Jaynes–Cummings model, Journal of Modern Optics40, 1195 (1993)

work page 1993

[15] [15]

M. W. Doherty, N. B. Manson, P. Delaney, F. Jelezko, J. Wrachtrup, and L. C. L. Hollenberg, The nitrogen-vacancy colour centre in diamond, Physics Reports528, 1 (2013)

work page 2013

[16] [16]

Amsüss, C

R. Amsüss, C. Koller, T. Nöbauer, S. Putz, S. Rotter, K. Sand- ner, S. Schneider, M. Schramböck, G. Steinhauser, H. Ritsch, J. Schmiedmayer, and J. Majer, Cavity QED with magnetically coupled collective spin states, Physical Review Letters107, 060502 (2011)

work page 2011

[17] [17]

Y . Kubo, F. R. Ong, P. Bertet, D. Vion, V . Jacques, D. Zheng, A. Dréau, J.-F. Roch, A. Auffeves, F. Jelezko, J. Wrachtrup, M. F. Barthe, P. Bergonzo, and D. Esteve, Strong coupling of a spin ensemble to a superconducting resonator, Physical Review Letters105, 140502 (2010)

work page 2010

[18] [18]

J. D. Breeze, E. Salvadori, J. Sathian, N. M. Alford, and C. W. M. Kay, Room-temperature cavity quantum electrody- namics with strongly coupled Dicke states, npj Quantum Infor- mation3, 40 (2017)

work page 2017

[19] [19]

Imamo ˘glu, Cavity QED based on collective magnetic dipole coupling: Spin ensembles as hybrid two-level systems, Phys

A. Imamo ˘glu, Cavity QED based on collective magnetic dipole coupling: Spin ensembles as hybrid two-level systems, Phys. Rev. Lett.102, 083602 (2009)

work page 2009

[20] [20]

D. I. Schuster, A. P. Sears, E. Ginossar, L. DiCarlo, L. Frunzio, J. J. L. Morton, H. Wu, G. A. D. Briggs, B. B. Buckley, D. D. Awschalom, and R. J. Schoelkopf, High-cooperativity coupling of electron-spin ensembles to superconducting cavities, Phys. Rev. Lett.105, 140501 (2010)

work page 2010

[21] [21]

E. R. Eisenach, J. F. Barry, M. F. O’Keeffe,et al., Cavity- enhanced microwave readout of a solid-state spin sensor, Na- ture Communications12, 1357 (2021)

work page 2021

[22] [22]

Zhang, Q

Y . Zhang, Q. Wu, S.-L. Su, Q. Lou, C. Shan, and K. Mølmer, Cavity quantum electrodynamics effects with nitrogen vacancy center spins coupled to room temperature microwave res- onators, Physical Review Letters128, 253601 (2022)

work page 2022

[23] [23]

W. Ng, H. Wu, and M. Oxborrow, Quasi-continuous cooling of a microwave mode on a benchtop using hyperpolarized NV- diamond, Applied Physics Letters119, 234001 (2021)

work page 2021

[24] [24]

D. P. Fahey, K. Jacobs, M. J. Turner, H. Choi, J. E. Hoffman, D. Englund, and M. E. Trusheim, Steady-state microwave mode 6 cooling with a diamond n-vensemble, Phys. Rev. Appl.20, 014033 (2023)

work page 2023

[25] [25]

Zhang, Q

Y . Zhang, Q. Wu, H. Wu, X. Yang, S.-L. Su, C. Shan, and K. Mølmer, Microwave mode cooling and cavity quantum elec- trodynamics effects at room temperature with optically cooled nitrogen-vacancy center spins, npj Quantum Information8, 125 (2022)

work page 2022

[26] [26]

Q. Wu, Y . Zhang, X. Yang,et al., A superradiant maser with nitrogen-vacancy center spins, Science China Physics, Mechan- ics & Astronomy65, 217311 (2022)

work page 2022

[27] [27]

Kersten, N

W. Kersten, N. de Zordo, O. Diekmann, E. S. Redchenko, et al., Self-induced superradiant masing, Nature Physics22, 158 (2026)

work page 2026

[28] [28]

H. Wang, K. L. Tiwari, K. Jacobs, M. Judy, X. Zhang, D. R. Englund, and M. E. Trusheim, A spin-refrigerated cavity quan- tum electrodynamic sensor, Nature Communications15, 10320 (2024)

work page 2024

[29] [29]

J. C. Hermann, R. Rizzato, F. Bruckmaier, R. D. Allert, A. Blank, and D. B. Bucher, Extending radiowave frequency detection range with dressed states of solid-state spin ensem- bles, npj Quantum Information10, 103 (2024)

work page 2024

[30] [30]

Z. Yin, J. J. Welter, C. A. Hart, P. V . Petruzzi,et al., High- resolution and wide-frequency-range magnetic spectroscopy with solid-state spin ensembles, npj Quantum Information11, 190 (2025)

work page 2025

[31] [31]

H. Wang, S. Wu, K. Jacobs, Y . S. Duan, D. R. Englund, and M. E. Trusheim, Cavity-enhanced solid-state nuclear spin gy- roscope, Physical Review Letters134, 183603 (2025)

work page 2025

[32] [32]

H. Wang, K. Jacobs, D. R. Englund, and M. E. Trusheim, Ultralow-power microwave frequency comb at a bistable phase transition, Phys. Rev. X16, 021005 (2026)

work page 2026

[33] [33]

H. Wang, K. Jacobs, D. Fahey, Y . Hu, D. R. Englund, and M. E. Trusheim, Exceptional sensitivity near the bistable tran- sition point of a hybrid quantum system, Nature Physics 10.1038/s41567-026-03217-3 (2026)

work page doi:10.1038/s41567-026-03217-3 2026

[34] [34]

Ovsiannikov, K

R. Ovsiannikov, K. Jacobs, A. G. Sotnikov, M. E. Trusheim, and D. I. Bondar, Hybrid non-degenerate parametric amplifier for a microwave cavity mode and an NV ensemble (2026), arXiv:2601.03407 [quant-ph]

work page arXiv 2026

[35] [35]

NV-ensemble enabled microwave/NV parametric amplifier with optimal driving

R. Ovsiannikov, K. Jacobs, A. G. Sotnikov, M. E. Trusheim, and D. I. Bondar, NV-ensemble enabled microwave/NV para- metric amplifier with optimal driving (2026), arXiv:2604.11643 [quant-ph]

work page internal anchor Pith review Pith/arXiv arXiv 2026

[36] [36]

Holstein and H

T. Holstein and H. Primakoff, Field dependence of the intrinsic domain magnetization of a ferromagnet, Physical Review58, 1098 (1940)

work page 1940

[37] [37]

L. G. Gunderman, T. W. Borneman, and D. G. Cory, Thermal state structure in the Tavis–Cummings model and rapid simu- lations in mesoscopic quantum ensembles, Phys. Rev. Research 7, 023155 (2025)

work page 2025

[38] [38]

D. I. Bondar, G. McCaul, and A. Sotnikov, Observables in Mo- tion: A guide to simulating classical and quantum dynamics (2025)

work page 2025

[39] [39]

Lambert, E

N. Lambert, E. Giguère, P. Menczel, B. Li, P. Hopf, G. Suarez, M. Gali Labarias, J. Lishman, R. Gadhvi, R. Agarwal, A. Gali- cia, N. Shammah, P. Nation, J. Johansson, S. Ahmed, S. Cross, A. Pitchford, and F. Nori, Qutip 5: The quantum toolbox in python, Physics Reports1153, 1 (2026)

work page 2026

[40] [40]

R. H. Dicke, Coherence in spontaneous radiation processes, Physical Review93, 99 (1954)

work page 1954

[41] [41]

For example, in Python,scipy.linalg.eigh_tridiagonal

work page

[42] [42]

Crank and P

J. Crank and P. Nicolson, A practical method for numerical evaluation of solutions of partial differential equations of the heat-conduction type, Mathematical Proceedings of the Cam- bridge Philosophical Society43, 50 (1947)

work page 1947

[43] [43]

Goldberg, H

A. Goldberg, H. M. Schey, and J. L. Schwartz, Computer- generated motion pictures of one-dimensional quantum- mechanical transmission and reflection phenomena, American Journal of Physics35, 177 (1967)

work page 1967

[44] [44]

L. H. Thomas,Elliptic Problems in Linear Difference Equations over a Network, Tech. Rep. (Watson Scientific Computing Lab- oratory, Columbia University, New York, 1949)

work page 1949

[45] [45]

Le Bris and P

C. Le Bris and P. Rouchon, Low-rank numerical approxima- tions for high-dimensional Lindblad equations, Physical Re- view A87, 022125 (2013)

work page 2013

[46] [46]

Finazzi, A

S. Finazzi, A. Le Boité, F. Storme, A. Baksic, and C. Ciuti, Corner-Space Renormalization Method for Driven-Dissipative Two-Dimensional Correlated Systems, Physical Review Letters 115, 080604 (2015)

work page 2015

[47] [47]

McCaul, K

G. McCaul, K. Jacobs, and D. I. Bondar, Fast computation of dissipative quantum systems with ensemble rank truncation, Physical Review Research3, 013017 (2021)

work page 2021

[48] [48]

Y .-T. Chen, C. Farquhar, and R. M. Parrish, Low-rank density- matrix evolution for noisy quantum circuits, npj Quantum In- formation7, 61 (2021)

work page 2021

[49] [49]

Donatella, Z

K. Donatella, Z. Denis, A. Le Boité, and C. Ciuti, Continuous- time dynamics and error scaling of noisy highly entangling quantum circuits, Physical Review A104, 062407 (2021)

work page 2021

[50] [50]

Suzuki, Fractal decomposition of exponential operators with applications to many-body theories and Monte Carlo simula- tions, Phys

M. Suzuki, Fractal decomposition of exponential operators with applications to many-body theories and Monte Carlo simula- tions, Phys. Lett. A146, 319 (2002)

work page 2002

[51] [51]

A. D. Bandrauk and H. Shen, Higher order exponential split operator method for solving time-dependent Schrödinger equa- tions, Canadian Journal of Chemistry70, 555 (1992)

work page 1992

[52] [52]

Ovsiannikov,https : / / github

R. Ovsiannikov,https : / / github . com / dibondar / NewPropagators(2026), gitHub repository, accessed 2026-04-22

work page 2026