We present a matrix-product-state-based numerical approach for simulating systems composed of several qubits and a common one-dimensional waveguide. In the presented approach, the one-dimensional waveguide is modeled in real space. Thus, one can use the advantage of matrix-product states that are suited for simulating low-entangled one-dimensional systems. The price to pay is that the vacuum of the waveguide in this modeling becomes the Bogoliubov vacuum, and one has to consider a not-so-small local Hilbert space for bosonic degrees of freedom. To manage the large local Hilbert space, we adopt the recently proposed single-site schemes. We demonstrate the potential of the presented approach by simulating superradiant phenomena within the Hamiltonian dynamics.
We present a matrix-product-state-based numerical approach for simulating systems composed of several qubits and a common one-dimensional waveguide. In the presented approach, the one-dimensional waveguide is modeled in real space. Thus, one can use the advantage of matrix-product states that are suited for simulating low-entangled one-dimensional systems. The price to pay is that the vacuum of the waveguide in this modeling becomes the Bogoliubov vacuum, and one has to consider a not-so-small local Hilbert space for bosonic degrees of freedom. To manage the large local Hilbert space, we adopt the recently proposed single-site schemes. We demonstrate the potential of the presented approach by simulating superradiant phenomena within the Hamiltonian dynamics.
A large, or even infinite, local Hilbert space dimension poses a significant computational challenge for simulating quantum systems. In this work, we present a matrix product state (MPS) based method for simulating one-dimensional quantum systems with a large local Hilbert space dimension, an example being bosonic systems with a large on-site population. To this end, we split the local Hilbert space corresponding to one site into two sites, each with a smaller Hilbert space dimension. An advantage of this method is that it can be easily integrated into MPS based techniques such as time-dependent variational principle (TDVP) without changing their standard algorithmic structure. Here, we implement our method using the TDVP to simulate the dynamics of the spin-boson model, a prototypical model of a spin interacting with a large bath of bosonic modes. We benchmark our method against and find excellent agreement with previous studies.
Simulating quantum many-body systems (QMBS) is one of the long-standing, highly non-trivial challenges in condensed matter physics and quantum information due to the exponentially growing size of the system's Hilbert space. To date, tensor networks have been an essential tool for studying such quantum systems, owing to their ability to efficiently capture the entanglement properties of the systems they represent. One of the well-known tensor network architectures, namely matrix product states (MPS), is the standard method for simulating one-dimensional QMBS. Here, we propose a simple, yet efficient, method to augment the already available MPS algorithms to simulate the dynamics of time-dependent Hamiltonians with better accuracy and a faster convergence rate, giving a second-order convergence compared to the first-order convergence of the standard method. We apply our proposed method to simulate the dynamics of a chain of single spins associated with nitrogen-vacancy color centers in diamonds, which has potential applications for practical and scalable quantum technologies, and find that our method improves the average error for a system of few NV centers by a factor of about 1000 for moderate step sizes. Our work paves the way for efficient simulation of QMBS under the influence of time-dependent Hamiltonians.
QwaveMPS is an open-source Python library for simulating one-dimensional quantum many-body waveguide systems using matrix product states (MPS). It provides a user-friendly interface for constructing, evolving, and analyzing quantum states and operators, facilitating studies in quantum physics and quantum information with waveguide-QED systems. This approach enables efficient, scalable simulations by focusing computational resources on the most relevant parts of the quantum system. Thus, one can study a wide range of complex dynamical interactions, including time-delayed feedback effects in the non-Markovian regime and deeply nonlinear systems, at a highly reduced computational cost compared to full Hilbert space approaches, making it both practical and convenient to model a variety of open waveguide-QED systems (in Markovian and non-Markovian regimes), treating quantized atoms and quantized photons on an equal footing. The QwaveMPS package automates the waveguide-QED-specific tensor-network operations, including construction of the time-bin evolution operator from physical parameters, the swap network required for delayed interactions, orthogonality-center and truncation management, the matrix-product factorization of N-photon Fock-state inputs, and the evaluation of single-time and two-time field observables (such as the second-order quantum correlation function) at chosen points along the waveguide.
In this thesis, the numerical simulation of strongly-interacting many-body quantum-mechanical systems using matrix product states (MPS) is considered. Compared to classical systems, quantum many-body systems possess an exponentially enlarged number of degrees of freedom, significantly complicating a simulation on a classical computer. Matrix-Product-States are a novel representation of arbitrary quantum many-body states. Using quantum information theory, it is possible to show that Matrix-Product-States provide a polynomial-sized representation of one-dimensional quantum systems, thus allowing an efficient simulation of one-dimensional quantum system on classical computers. Matrix-Product-States form the conceptual framework of the density-matrix renormalization group (DMRG). Based upon this connection, deeper understanding of the density matrix renormalization group can be obtained. After a general introduction in the first chapter of this thesis, the second chapter deals with Matrix-Product-States, focusing on the development of fast and stable algorithms. It is possible to extend the Matrix-Product-States approach to be able to represent arbitrary operators, the so-called Matrix-Product-Operators, which allows a fast and flexible calculation of arbitrary expectation values. To obtain algorithms to efficiently calculate groundstates, the density-matrix renormalization group is reformulated using the Matrix-Product-States framework. Further, time-dependent problems are considered. Two different algorithms are presented, one based on a Trotter decomposition of the time-evolution operator, the other one on Krylov subspaces. Finally, the evaluation of dynamical spectral functions is discussed, and a correction vector-based method is presented. In the following chapters, the methods presented in the second chapter, are applied to a number of different physical problems. The third chapter deals with the existence of chiral phases in isotropic one-dimensional quantum spin systems. A preceding analytical study based on a mean-field approach indicated the possible existence of those phases in an isotopic Heisenberg model with a frustrating zig-zag interaction and a magnetic field. In this thesis, the existence of the chiral phases will be shown numerically by using Matrix-Product-States-based algorithms. A key effect of interacting one-dimensional quantum-mechanical many-body systems is the spin-charge separation. However, up to now only signs of the spin-charge separation have been observed in experiments. In the fourth chapter, we propose an experiment using ultracold atomic gases in optical lattices, which allows a well controlled observation of the spin-charge separation (of different hyperfine states of the ultracold atoms) with current state of the art experimental techniques. Ultracold atoms in optical lattices are well described by (Bose)-Hubbard models. In order to support this proposal, we present numerical results for realistic system parameters. Matrix-Product-States are an excellent tool for the simulation of one-dimensional quantum systems, however, they are not well suited for the simulation of higher dimensional systems. For strongly-correlated systems, for instance cuprates-based high-temperature superconductors, quantum fluctuations play an essential role. Classical mean-field theories neglect any kind of fluctuations, thus they are not suitable to describe strongly-correlated systems. The dynamical mean-field theory (DMFT) fully takes local quantum fluctuations into account but neglects any kind of spatial fluctuations. The many-body problem on the lattice is mapped onto an impurity problem, which needs to be solved self-consistently. In the last chapter of this thesis, Matrix-Product-States-based algorithms are used to solve the impurity problem of the dynamical mean-field theory. We present results for a Hubbard model on a one-dimensional lattice and on a Bethe lattice obtained by the dynamical mean-field and compare them with exact results.
Ultrashort pulses propagating in nonlinear nanophotonic waveguides can simultaneously leverage both temporal and spatial field confinement, promising a route towards single-photon nonlinearities in an all-photonic platform. In this multimode quantum regime, however, faithful numerical simulations of pulse dynamics naïvely require a representation of the state in an exponentially large Hilbert space. Here, we employ a time-domain, matrix product state (MPS) representation to enable efficient simulations by exploiting the entanglement structure of the system. To extract physical insight from these simulations, we develop an algorithm to unravel the MPS quantum state into constituent temporal supermodes, enabling, e.g., access to the phase-space portraits of arbitrary pulse waveforms. As a demonstration, we perform exact numerical simulations of a Kerr soliton in the quantum regime. We observe the development of non-classical Wigner-function negativity in the solitonic mode as well as quantum corrections to the semiclassical dynamics of the pulse. A similar analysis of χ ( 2 ) simultons reveals a unique entanglement structure between the fundamental and second harmonics. Our approach is also readily compatible with quantum trajectory theory, allowing full quantum treatment of propagation loss and decoherence. We expect this work to establish the MPS technique as part of a unified engineering framework for the emerging field of broadband quantum photonics.
Modeling open quantum systems---quantum systems coupled to a bath---is of value in condensed-matter theory, cavity quantum electrodynamics, nanosciences, and biophysics. The real-time simulation of open quantum systems was advanced significantly by the recent development of chain mapping techniques and the use of matrix product states that exploit the intrinsic entanglement structure in open quantum systems. The computational cost of simulating open quantum systems, however, remains high when the bath is excited to high-lying quantum states. We develop an approach to reduce the computational costs in such cases. The interaction representation for the open quantum system is used to distribute excitations among the bath degrees of freedom so that the occupation of each bath oscillator is ensured to be low. The interaction picture also causes the matrix dimensions to be much smaller in a matrix product state of a chain-mapped open quantum system than in the Schr\"odinger picture. Using the interaction representation accelerates the calculations by one to two orders of magnitude over the existing matrix-product-state method. In the regime of strong system-bath coupling and high temperatures, the speedup can be as large as three orders of magnitude. The approach developed here is especially promising to simulate the dynamics of open quantum systems in high-temperature and strong-coupling regimes.