In this work we introduce a phase-space description based on the positive P representation for bosonic fields interacting with a system of quantum emitters. The formalism is applicable to collective light-matter interactions and open quantum systems with decoherence. Conservation of particle numbers is considered, and a Jordan-Schwinger transformation enables the representation of multi-level quantum emitters. The evolution of the phase-space description of the combined system of emitters and field is formulated in terms of stochastic trajectories and we derive the rules of mapping from traditional quantum mechanics to this stochastic formalism. The resulting equations of motion encode deterministic, classical evolution with quantum effects incorporated by stochastic noise terms. The framework's equations and properties are provided without specifying the Hamiltonian, aiming for broad applicability in diverse research domains. A potential future application is the quantum mechanical description of collective spontaneous emission of an incoherently pumped ensemble of atoms.
In this work we introduce a phase-space description based on the positive P representation for bosonic fields interacting with a system of quantum emitters. The formalism is applicable to collective light-matter interactions and open quantum systems with decoherence. Conservation of particle numbers is considered, and a Jordan-Schwinger transformation enables the representation of multi-level quantum emitters. The evolution of the phase-space description of the combined system of emitters and field is formulated in terms of stochastic trajectories and we derive the rules of mapping from traditional quantum mechanics to this stochastic formalism. The resulting equations of motion encode deterministic, classical evolution with quantum effects incorporated by stochastic noise terms. The framework's equations and properties are provided without specifying the Hamiltonian, aiming for broad applicability in diverse research domains. A potential future application is the quantum mechanical description of collective spontaneous emission of an incoherently pumped ensemble of atoms.
We critically explore the applicability of a recently proposed framework to sample the quantum dynamics of a many-body quantum system interacting with light by stochastic trajectories, applying it to the closed and open Tavis-Cummings model (TCM). The stochastic differential equations (SDEs) sample the positive P phase-space representation by analog complex-valued dynamical variables that are linked to the quantum operators. Statistical average over the stochastic trajectories yields the evolution of the quantum mechanical expectation values. However, numerical implementation of these SDEs for the TCM indicates divergent solutions, also known from other phase-space methods. This limits the applicability of the framework to finite propagation times, that are strongly dependent on the physical parameters and initial conditions of the system. We outline the underlying mathematical reason for these divergences and show that their contribution to the averages are, however, essential. To attempt to regularize the divergences, we transform the SDEs to an equivalent set of SDEs with different noise realisations, thereby pushing the valid time boundary. Quantum collapse and revival of the TCM, however, cannot be recovered by the stochastic trajectory approach, pointing to the general difficulty of the applicability of stochastic phase-space sampling methods to systems with strong quantum features.
Classical Brownian motion has well been investigated since the pioneering work of Einstein, which inspired mathematicians to lay the theoretical foundation of stochastic processes. A stochastic formulation for quantum dynamics of dissipative systems described by the system-plus-bath model has been developed and found many applications in chemical dynamics, spectroscopy, quantum transport, and other fields. This article provides a tutorial review of the stochastic formulation for quantum dissipative dynamics. The key idea is to decouple the interaction between the system and the bath by virtue of the Hubbard-Stratonovich transformation or Itô calculus so that the system and the bath are not directly entangled during evolution, rather they are correlated due to the complex white noises introduced. The influence of the bath on the system is thereby defined by an induced stochastic field, which leads to the stochastic Liouville equation for the system. The exact reduced density matrix can be calculated as the stochastic average in the presence of bath-induced fields. In general, the plain implementation of the stochastic formulation is only useful for short-time dynamics, but not efficient for long-time dynamics as the statistical errors go very fast. For linear and other specific systems, the stochastic Liouville equation is a good starting point to derive the master equation. For general systems with decomposable bath-induced processes, the hierarchical approach in the form of a set of deterministic equations of motion is derived based on the stochastic formulation and provides an effective means for simulating the dissipative dynamics. A combination of the stochastic simulation and the hierarchical approach is suggested to solve the zero-temperature dynamics of the spin-boson model. This scheme correctly describes the coherent-incoherent transition (Toulouse limit) at moderate dissipation and predicts a rate dynamics in the overdamped regime. Challenging problems such as the dynamical description of quantum phase transition (local- ization) and the numerical stability of the trace-conserving, nonlinear stochastic Liouville equation are outlined.
We propose a stochastic Schrödinger equation in which the momentum is coupled to a white Gaussian noise. In the stochastic representation, the kinetic energy representing the self-interaction of momentum is reduced to a linear term of momentum. As such, the quantum evolution operator factorizes into two contributions due to the momentum and the potential, respectively. The exact quantum propagator thereby becomes an expectation of the stochastic one in which the amplitude results from the potential with a fluctuating position: The particle moves with a constant velocity, subjected to a complex Brownian motion. We demonstrate that the stochastic Schrödinger equation can be feasibly used to derive the quantum propagators for the linear potential and the harmonic oscillator system. Novel semiclassical and other approximations may be developed from the new representation of quantum mechanics.
Abstract This work introduces a theoretical framework to model the collective dynamics of quantum emitters in highly non-Markovian environments, interacting through the exchange of photons with significant retardations. The formalism consists on a set of coupled delay differential equations for the emitter’s raising/lowering operators σ i ± , supplemented by input–output relations that describe the field mediating the interactions. These equations capture the dynamics of both linear (bosonic) and nonlinear (two-level) emitter arrays. It is exact in some limits—e.g. bosonic emitters or generic systems with up to one collective excitation—and can be integrated to provide accurate results for larger numbers of photons. These equations support a study of collective spontaneous emission of emitter arrays in open waveguide-QED environments. This study uncovers an effect we term cascaded super- and sub-radiance, characterized by light-cone-limited propagation and increasingly correlated photon emission across distant emitters. The collective nature of this dynamics for two-level systems is evident both in the enhancement of collective emission rates, as well as in a superradiant burst with a faster than linear growth. While these effects should be observable in existing circuit QED devices or slight generalizations thereof, the formalism put forward in this work can be extended to model other systems such as network of quantum emitters or the generation of correlated photon states.
The stochastic Schrödinger equation (SSE) provides a trajectory-level route to simulate the dynamics of open quantum systems with applications ranging from molecular processes to quantum technologies. We study a colored-noise extension of the SSE based on an Ornstein-Uhlenbeck (OU) noise drive and benchmark its ensemble-averaged dynamics against the standard white-noise SSE and against a fluctuating OU random Hamiltonian. When the environment exhibits a finite correlation time, averaging over pure-state trajectories yields master equations that are generally open-form and not of Lindblad type, yet remain positive by construction. By considering the differential of the OU process, we define an effective correlated noise, whose properties we analyze and use to formulate an SSE unraveling of its associated open-form quantum master equation. We show that the averaged dissipator separates into a Lindblad contribution stemming from the white-noise component and additional correlation terms arising from the fluctuations of the OU Hamiltonian. To obtain a practical closed description and physical intuition, we introduce a Redfield-inspired perturbative closure for these correlation terms, providing an effective master equation for the colored SSE. For a two-level system, the resulting dynamics exhibit long-lived coherences, nontrivial stationary (including oscillatory) states, and multi-timescale relaxation, rationalized through the components of a time-dependent Redfield tensor.
We consider a class of open quantum many-body Lindblad dynamics characterized by an all-to-all coupling Hamiltonian and by dissipation featuring collective ``state-dependent" rates. The latter encodes local incoherent transitions that depend on average properties of the system. This type of open quantum dynamics can be seen as a generalization of classical (mean-field) stochastic Markov dynamics, in which transitions depend on the instantaneous configuration of the system, to the quantum domain. We study the time evolution in the limit of infinitely large systems, and we demonstrate the exactness of the mean-field equations for the dynamics of average operators. We further derive the effective dynamical generator governing the time evolution of (quasi-)local operators. Our results allow for a rigorous and systematic investigation of the impact of quantum effects on paradigmatic classical models, such as quantum generalized Hopfield associative memories or (mean-field) kinetically-constrained models.