• A study published on arXiv introduces tensor network methodologies to resolve the chronic memory expansion problem in radical-pair spin dynamics simulations. • Utilizing Matrix Product States (MPS) and Matrix Product Density Operators (MPDO), the researchers accurately calculated open quantum dynamics for systems containing up to 60 spins. • Applied to the flavin-tryptophan (FAD-Trp) system, which is central to avian magnetoreception, the framework precisely captured anisotropic spin dynamics and electron-transfer pathways. [Quantum Biology Society] Simulating the quantum dynamics of radical pairs (pairs of spin qubits) has long faced a significant technological barrier in the fields of quantum biology and spin chemistry. When attempting to include all surrounding nuclear spins in calculations, the Hilbert space expands exponentially, quickly exceeding the memory limits of current computing systems. However, a groundbreaking study has recently been published that fundamentally overcomes this "curse of dimensionality" by introducing tensor network techniques previously used for many-body quantum systems. A paper titled "Introduction to modelling radical pair quantum spin dynamics with tensor networks," released on the preprint repository arXiv, presents a revolutionary solution to these massive computational limits. According to the original publication, the development of this framework was led by a joint research team including Kentaro Hino and Yuki Kurashige from Kyoto University (Japan), Damyan S. Frantzov from the University of Oxford (UK), and Lewis M. Antill from Sungkyunkwan University (Republic of Korea). ■ Memory Reduction via Tensor Networks To efficiently compress the explosively growing Hilbert and Liouville spaces, the research team introduced representations based on Matrix Product States (MPS) and Matrix Product Density Operators (MPDO). This enabled the simulation of full open quantum dynamics while explicitly accounting for dozens of nuclear spins—a feat previously considered impossible with traditional computing methods. The team successfully benchmarked the methodology on systems with up to 60 spins, proving the overwhelming efficiency of the tensor network approach. ■ Demonstration in the Flavin-Tryptophan System To verify the algorithm's validity in a real biological context, the researchers applied it to the flavin-tryptophan (FAD-Trp) radical-pair model, a key biological component linked to the avian magnetoreception mechanism. The analysis successfully identified minute dependencies on the direction of an external magnetic field and captured anisotropic spin dynamics and electron-transfer pathways with high precision. This effectively reproduced the complex quantum dynamics of an actual biological environment, where numerous nuclear spin hyperfine interactions are intricately entangled, without distortion. ■ New Horizons for Quantum Biology Computing This study demonstrates that tensor networks are not just theoretical tools for physics but serve as a highly powerful and practical computing framework that enables ultra-large-scale quantum information processing in spin chemistry and quantum biology. As this technology extends into fields like drug discovery and quantum sensing, it is expected to drastically accelerate the analysis of precise spin dynamics within large protein complexes—tasks that were once impossible due to computational constraints.
• A study published on arXiv introduces tensor network methodologies to resolve the chronic memory expansion problem in radical-pair spin dynamics simulations. • Utilizing Matrix Product States (MPS) and Matrix Product Density Operators (MPDO), the researchers accurately calculated open quantum dynamics for systems containing up to 60 spins. • Applied to the flavin-tryptophan (FAD-Trp) system, which is central to avian magnetoreception, the framework precisely captured anisotropic spin dynamics and electron-transfer pathways. [Quantum Biology Society] Simulating the quantum dynamics of radical pairs (pairs of spin qubits) has long faced a significant technological barrier in the fields of quantum biology and spin chemistry. When attempting to include all surrounding nuclear spins in calculations, the Hilbert space expands exponentially, quickly exceeding the memory limits of current computing systems. However, a groundbreaking study has recently been published that fundamentally overcomes this "curse of dimensionality" by introducing tensor network techniques previously used for many-body quantum systems. A paper titled "Introduction to modelling radical pair quantum spin dynamics with tensor networks," released on the preprint repository arXiv, presents a revolutionary solution to these massive computational limits. According to the original publication, the development of this framework was led by a joint research team including Kentaro Hino and Yuki Kurashige from Kyoto University (Japan), Damyan S. Frantzov from the University of Oxford (UK), and Lewis M. Antill from Sungkyunkwan University (Republic of Korea). ■ Memory Reduction via Tensor Networks To efficiently compress the explosively growing Hilbert and Liouville spaces, the research team introduced representations based on Matrix Product States (MPS) and Matrix Product Density Operators (MPDO). This enabled the simulation of full open quantum dynamics while explicitly accounting for dozens of nuclear spins—a feat previously considered impossible with traditional computing methods. The team successfully benchmarked the methodology on systems with up to 60 spins, proving the overwhelming efficiency of the tensor network approach. ■ Demonstration in the Flavin-Tryptophan System To verify the algorithm's validity in a real biological context, the researchers applied it to the flavin-tryptophan (FAD-Trp) radical-pair model, a key biological component linked to the avian magnetoreception mechanism. The analysis successfully identified minute dependencies on the direction of an external magnetic field and captured anisotropic spin dynamics and electron-transfer pathways with high precision. This effectively reproduced the complex quantum dynamics of an actual biological environment, where numerous nuclear spin hyperfine interactions are intricately entangled, without distortion. ■ New Horizons for Quantum Biology Computing This study demonstrates that tensor networks are not just theoretical tools for physics but serve as a highly powerful and practical computing framework that enables ultra-large-scale quantum information processing in spin chemistry and quantum biology. As this technology extends into fields like drug discovery and quantum sensing, it is expected to drastically accelerate the analysis of precise spin dynamics within large protein complexes—tasks that were once impossible due to computational constraints.
Radical pairs are short-lived, spin-correlated intermediates that underpin processes in chemistry, biology, and emerging quantum technologies. Their behavior is governed by coupled electron-nuclear spin dynamics and is sensitive to weak magnetic fields, but full quantum treatments have been obstructed by the extreme computational cost of modeling many interacting spins. Here, this barrier is removed, demonstrating that open-system radical-pair dynamics can be resolved at nuclear-spin scales previously intractable, explicitly reaching regimes with tens of coupled nuclei and validated up to 60 spins. In biologically relevant flavin-tryptophan radical pairs, electron-transfer pathways and magnetic-field anisotropy are found to reshape spin evolution and, in turn, spin-selective reaction yields. The resulting directional responses expose a strong mechanistic link between the nuclear environment, magnetic geometry, and chemical outcome-a relationship central to hypotheses of avian magnetoreception and other magnetic-field effects in biology. This establishes a widely applicable simulation framework that removes a long-standing barrier in spin chemistry, quantum biology, and spin-based device science.
• Published in AIP Publishing's APL Quantum, fully applying digital quantum simulation (DQS) to spin biology. • Successfully implemented the quantum spin dynamics of the radical pair mechanism utilizing the 'Trotterization' technique. • Analyzed the noise limits of actual quantum hardware and demonstrated the potential of utilizing quantum computing for the future development of new drugs targeting complex biomolecules. [Quantum Biology Society] In the emerging field of quantum biology, simulating spin-chemistry systems with actual biological complexity has posed computational limitations that are difficult to overcome with conventional classical computers. This is because as the number of elements constituting the spin system increases, the dimension of the Hilbert space required for computation grows exponentially. Recently, a study simulating biological phenomena using an actual digital quantum computer was published to resolve this problem, drawing the attention of the academic community. A paper titled "Simulating spin biology using a digital quantum computer: Prospects on a near-term quantum hardware emulator" was published on August 19 in APL Quantum, a journal of AIP Publishing. This research was led by a multinational joint research team including Pedro H. Alvarez and Clarice D. Aiello. ■ Mapping the Radical Pair Mechanism via Trotterization The radical pair mechanism (RPM), represented by the magnetic field sensing (magnetoreception) ability of migratory birds, centers on the quantum entanglement and coherence phenomena of electron spins in vivo. The research team introduced the Trotterization technique to directly apply gate-based digital quantum simulation (DQS) to biological systems. Through this technique, they successfully mapped the coherent spin dynamics of radical pairs onto quantum circuits. In ideal, noise-free quantum circuit simulations, they obtained accurate spin dynamics results consistent with the conventional classical Haberkorn master equation approach using only about 15 Trotter steps. ■ Analysis of Quantum Hardware Noise Limits and Future Prospects Furthermore, going beyond ideal environments, the research team conducted an in-depth analysis of the impact of realistic error models on simulations using a near-term (NISQ) quantum hardware emulator that mimics the noise profile of an IBM quantum device (e.g., ibmq_jakarta). The analysis revealed that realistic limitations were observed due to the constrained performance and noise (such as phase damping and amplitude damping) of current quantum hardware, including a shift in singlet yield and a decreased ability to sense magnetic field anisotropy. However, the researchers showed that by fitting an analytical model to the measured noisy data, the overall dynamical shape depending on the magnetic field angle could be successfully derived. In addition, they presented a practical upper limit for the number of Trotter steps that can be directly run on quantum devices if error mitigation techniques are introduced in the future. Unlike conventional classical computing methods that face exponential limits, the DQS approach possesses the potential to efficiently model biological spin systems simply through linear resource scaling by adding qubits. This study is evaluated as a groundbreaking first step demonstrating that digital quantum simulation can be utilized as an innovative tool for new drug development targeting complex biomolecules and ultra-precise biosensing modeling by directly calculating the entanglement and coherence of biological phenomena.
This paper explores the application of tensor network methods for simulating molecular dynamics on quantum computers. Traditional molecular dynamics simulations are computationally expensive, particularly for complex systems exhibiting strong quantum correlations. Tensor networks provide a powerful framework for representing and manipulating high-dimensional quantum states, potentially offering a pathway to efficiently simulate molecular dynamics by directly leveraging quantum computational resources. This work details the theoretical framework, outlining how tensor network techniques, specifically Matrix Product States (MPS) and Tree Tensor Networks (TTN), can be adapted for this purpose. We present a simplified model of a diatomic molecule to illustrate the core concept: the construction of a tensor network representation of the molecular wavefunction, followed by the application of time evolution operators to simulate dynamics. The focus is on demonstrating the potential of this approach rather than a fully optimized simulation, highlighting key challenges and future research directions. The efficiency gains compared to classical methods are discussed, acknowledging the significant technological hurdles involved in realizing this approach on near-term quantum devices. We investigate the scalability of the method, considering the limitations imposed by the bond dimension of the tensor network. The ultimate goal is to demonstrate how tensor networks can provide a viable strategy for tackling computationally intractable molecular dynamics problems, paving the way for a new era of quantum-enhanced molecular simulations.
The study of interacting quantum many-body systems is one of the central challenges in modern physics. Interestingly, by understanding the physical properties of such systems from a microscopic description, one can explain many complex quantum phenomena such as the high-temperature superconductivity, fractional quantum Hall effect, etc. However, due to high complexity of these systems, theoretical investigations such as analytical and perturbation theories are limited to studying only a small subset of problems. In the past few decades, numerical simulation has emerged as an indispensable method for studying many phenomena, and has led to immense success and progress in recent years. The main obstacle in addressing quantum many-body systems is that the number of parameters required to represent a state of the system scales exponentially with the system size. However, one hardly needs the full Hilbert space but only a small subspace to describe physical properties of the system. From the point of view in quantum information theory, quantum entanglement is one of the main factors that is responsible for the complexity of studying quantum many-body systems. In particular, the amount of quantum entanglement in such system quantifies the correlation between constituents and therefore it determines how big the subspace is required to represent the system. In this thesis, we are interested in using tensor network states, a new emerging numerical technique to study dynamical and ground state properties of the system. The key ingredient is that a tensor network state which consists of a fairly small number of parameters can represent very well the complete wave function of a system. In one spatial dimension, Matrix Product State (MPS) has become a standard tool, combined with related algorithms such as Time Evolving Block Decimation algorithm (TEBD) to study dynamical properties of a locally-perturbed quantum lattice system. In this thesis, we introduce a technique called ``infinite boundary conditions'' that allows us to use a finite-size MPS to represent the state of an infinite-size system. As a result, one can simulate dynamics of a system for a longer time with a significant computational cost reduction compared to traditional methods. In two spatial dimensions, the tensor network known as Projected Entangled Pair States (PEPS) is emerging as a useful approach. However algorithms for optimizing the matrix elements of these tensors are complex and not very efficient. In particular, the efficiency of the optimization depends critically on the quality of the so-called ``environment'', which is expensive to compute. We have introduced a new method to accelerate the convergence of calculating the PEPS ground state by recycling the environment tensors for many optimization steps, vastly reducing the number of times that the full environment needs to be calculated.
Abstract We present a collection of methods to simulate entangled dynamics of open quantum systems governed by the Lindblad master equation with tensor network methods. Tensor network methods using matrix product states (MPSs) have been proven very useful to simulate many-body quantum systems and have driven many innovations in research. Since the MPS design is tailored for closed one-dimensional systems governed by the Schrödinger equation, the next step for many-body quantum dynamics is the simulation of one-dimensional open quantum systems. We review the three dominant approaches to the simulation of open quantum systems via the Lindblad master equation: quantum trajectories, matrix product density operators (MPDOs), and locally purified tensor networks. Selected examples guide possible applications of the methods and serve moreover as a benchmark between the techniques. These examples include the finite temperature states of the transverse quantum Ising model, the dynamics of an exciton traveling under the influence of spontaneous emission and dephasing, and a double-well potential simulated with the Bose–Hubbard model including dephasing. We analyze which approach is favorable leading to the conclusion that a complete set of all three methods is most beneficial, pushing the limits of different scenarios. The convergence studies using analytical results for macroscopic variables and exact diagonalization methods as comparison, show, for example, that MPDOs are favorable for the exciton problem in our study. All three methods access the same library, i.e., the software package Open Source Matrix Product States , allowing us to have a meaningful comparison between the different approaches based on the selected examples. For example, tensor operations are accessed from the same subroutines and with the same optimization eliminating one possible bias in a comparison of such numerical methods.
This thesis treats the classical simulation of strongly-interacting many-body quantum-mechanical systems in more than one dimension using matrix product states and the more general tensor product states. Contrary to classical systems, quantum many-body systems possess an exponentially larger number of degrees of freedom, thereby significantly complicating their numerical treatment on a classical computer. For this thesis two different representations of quantum many-body states were employed. The first, the so-called matrix product states (MPS) form the basis for the extremely successful density matrix renormalization group (DMRG) algorithm. While originally conceived for one-dimensional systems, MPS are in principle capable of describing arbitrary quantum many-body states. Using concepts from quantum information theory it is possible to show that MPS provide a representation of one-dimensional quantum systems that scales polynomially in the number of particles, therefore allowing an efficient simulation of one-dimensional systems on a classical computer. One of the key results of this thesis is that MPS representations are indeed efficient enough to describe even large systems in two dimensions, thereby enabling the simulation of such systems using DMRG. As a demonstration of the power of the DMRG algorithm, it is applied to the Heisenberg antiferromagnet with spin S=1/2 on the kagome lattice. This model's ground state has long been under debate, with proposals ranging from static spin configurations to so-called quantum spin liquids, states where quantum fluctuations destroy conventional order and give rise to exotic quantum orders. Using a fully SU(2)-symmetric implementation allowed us to handle the exponential growth of entanglement and to perform a large-scale study of this system, finding the ground state for cylinders of up to 700 sites. Despite employing a one-dimensional algorithm for a two-dimensional system, we were able to compute the spin gap (i.e. the energy gap to the first spinful excitation) and study the ground state properties, such as the decay of correlation functions, the static spin structure factors, and the structure and distribution of the nearest-neighbor spin-spin correlations. Additionally, by applying a new tool from quantum information theory, the topological entanglement entropy, we could also with high confidence demonstrate the ground state of this model to be the elusive gapped Z2 quantum spin liquid with topological order. To complement this study, we also considered the extension of MPS to higher dimensions, known as tensor product states (TPS). We implemented an optimization algorithm exploiting symmetries for this class of states and applied it to the bilinear-biquadratic-bicubic Heisenberg model with spin S=3/2 on the z=3 Bethe lattice. By carefully analyzing the simulation data we were able to determine the presence of both conventional and symmetry-protected topological order in this model, thereby demonstrating the analytically predicted existence of the Haldane phase in higher dimensions within an extended region of the phase diagram. Key properties of this symmetry-protected topological order include a doubling of the levels in the entanglement spectrum and the presence of edge spins, both of which were confirmed in our simulations. This finding simultaneously validated the applicability of the novel TPS algorithms to the search for exotic order.