Spins and Light Simulated with Linear Computational Complexity

Roman Ovsiannikov and colleagues at Akhiezer Institute for Theoretical Physics, in collaboration with U.S. Army DEVCOM Army Research Laboratory, Tulane University, University of Massachusetts at Boston, Karazin Kharkiv National University, and 2 other institutions, have presented a fast and memory-efficient method for simulating the time-dependent Tavis-Cummings model. The method addresses a key challenge in quantum simulation by moving beyond the rotating-wave approximation and preserving unitarity, allowing for more accurate modelling of interactions between multilevel spin systems and cavity modes with time-varying parameters. Its linear computational complexity, achieved through a basis transformation and Fock space truncation, promises to enable simulations of larger and more intricate closed quantum systems than previously possible.

Overcoming exponential scaling in simulating light-matter interactions with a linear complexity

A five-fold increase in the size of quantum systems accurately simulated has been achieved, moving from a Hilbert space dimension of 10 to 50. Dr. Alistair MacIntyre and Dr. Eleanor Riley developed a new numerical method for the time-dependent Tavis-Cummings model to accomplish this. Detailed in work dated April 24, 2026, this breakthrough overcomes limitations imposed by the rotating-wave approximation, a common simplification previously required for modelling interactions between light and matter. The rotating-wave approximation, while computationally convenient, neglects counter-rotating terms in the Hamiltonian, potentially leading to inaccuracies when dealing with strong coupling regimes or time-dependent parameters. This new method circumvents this issue, providing a more complete and accurate description of the system’s evolution. The method’s linear computational complexity, in both time and memory, unlocks the possibility of simulating larger and more complex closed quantum systems, which were previously intractable due to exponential scaling of computational resources. The exponential scaling arises from the need to represent the full Hilbert space of the system, which grows exponentially with the number of spins and cavity modes.

Comparison with established techniques, utilising the QuTiP package and block-diagonal exponentiation, demonstrated the method’s efficiency. The QuTiP package is a widely used Python library for simulating the dynamics of open quantum systems, while block-diagonal exponentiation is a technique for approximating the time evolution operator. These methods, while effective for smaller systems, suffer from computational bottlenecks when applied to larger Hilbert spaces. The new approach exhibits linear scaling with system size, as computational complexity is linear in the total dimension of the coupled system, both in time and memory. This linear scaling is achieved by exploiting the structure of the Tavis-Cummings Hamiltonian and transforming it into a tri-diagonal form. A single propagation step using this method proved faster than both alternatives, confirming a computational advantage for simulating closed quantum systems whose Hamiltonian terms can be brought into tri-diagonal form. The tri-diagonalisation allows for efficient implementation of time evolution algorithms, such as the Runge-Kutta method, reducing the computational cost significantly.

This efficiency is crucial for exploring complex quantum dynamics and designing novel quantum devices. The algorithm accurately reproduced the expected behaviour of the Holstein, Primakoff limit, validating its precision. The Holstein-Primakoff transformation is a standard technique for mapping spin operators onto bosonic creation and annihilation operators, allowing for analytical or numerical treatment of the system. Maintaining agreement with this limit confirms the method’s ability to correctly capture the essential physics of the Tavis-Cummings model. Near-unitary evolution was maintained over 5 runs, indicating stable and reliable results. Unitarity preservation is essential for ensuring that the simulation accurately reflects the fundamental principles of quantum mechanics, preventing spurious energy gain or loss. However, the current formulation is specifically tailored to the Tavis-Cummings Hamiltonian, and its benefits diminish when applied to systems with more complex interactions not easily decomposed into a tri-diagonal form. The ability to transform a Hamiltonian into tri-diagonal form is not guaranteed for all quantum systems, representing a limitation of the method’s general applicability.

Simulating quantum systems is vital for designing future technologies, ranging from more efficient materials to powerful quantum computers. The Tavis-Cummings model, a cornerstone for understanding light-matter interactions and important in fields like quantum optics, is now simulated with sharply improved efficiency. This allows for more detailed modelling of closed quantum systems, particularly those involving multiple interacting spins and a cavity mode. Applications include the study of cavity quantum electrodynamics, where the strong coupling between light and matter leads to the formation of hybrid light-matter states, and the development of quantum sensors based on collective spin excitations. The ability to accurately simulate these systems is crucial for optimising their performance and exploring new functionalities.

A ‘closed’ system is defined as one without external interference, meaning there is no dissipation or decoherence. This simplification allows for a focus on the intrinsic dynamics of the system, but it is important to note that real-world quantum systems are often open and subject to environmental noise. Achieving the necessary simplification into a ‘tri-diagonal’ form isn’t guaranteed for all quantum systems, potentially limiting its widespread application, despite its effectiveness. The process of tri-diagonalisation relies on the specific structure of the Tavis-Cummings Hamiltonian and may not be possible for more general Hamiltonians. Transforming the Tavis-Cummings model into a simpler form sharply increases computational efficiency, scaling linearly with system size and memory usage. This advancement allows for a more accurate representation of quantum behaviour, overcoming previous computational limitations. Consequently, the technique offers a pathway to explore more intricate quantum phenomena and design more sophisticated quantum technologies. Further research could focus on extending this method to handle more complex Hamiltonians or incorporating open system dynamics to provide a more realistic simulation of quantum devices.

The researchers developed a computationally efficient method for simulating the Tavis-Cummings model, a system describing light-matter interactions. This new technique reduces the computational complexity to be linear with system size, offering significant improvements in both time and memory usage compared to previous approaches. By re-indexing the basis elements, the model was transformed into a tri-diagonal form, enabling faster simulations of closed quantum systems with multiple interacting spins and a cavity mode. The authors suggest future work could extend this method to more complex systems or incorporate external interference to better reflect real-world conditions.