We propose two methods for computing the large deviations of the first-passage-time statistics in general open quantum systems. The first method determines the region of convergence of the joint Laplace transform and the z transform of the first-passage time distribution by solving an equation of poles with respect to the z-transform parameter. The scaled cumulant generating function, which is the logarithm of the boundary values within this region, is subsequently obtained. The theoretical basis is that the dynamics of open quantum systems can be unraveled into a piecewise deterministic process and that a tilted Liouville master equation exists in Hilbert space. The second method uses a simulation-based approach that is built on the wave-function cloning algorithm. To validate both methods, we derive analytical expressions for the scaled cumulant generating functions in field-driven two-level and three-level systems. In addition, we present numerical results and cloning simulations for a field-driven system comprising two interacting two-level atoms.
We propose two methods for computing the large deviations of the first-passage-time statistics in general open quantum systems. The first method determines the region of convergence of the joint Laplace transform and the z transform of the first-passage time distribution by solving an equation of poles with respect to the z-transform parameter. The scaled cumulant generating function, which is the logarithm of the boundary values within this region, is subsequently obtained. The theoretical basis is that the dynamics of open quantum systems can be unraveled into a piecewise deterministic process and that a tilted Liouville master equation exists in Hilbert space. The second method uses a simulation-based approach that is built on the wave-function cloning algorithm. To validate both methods, we derive analytical expressions for the scaled cumulant generating functions in field-driven two-level and three-level systems. In addition, we present numerical results and cloning simulations for a field-driven system comprising two interacting two-level atoms.
In a specific class of open quantum systems with finite and fixed numbers of collapsed quantum states, the semi-Markov process method is used to calculate the large deviations of the first passage time statistics. The core formula is an equation of poles, which is also applied in determining the scaled generating functions (SCGFs) of the counting statistics. For simple counting variables, the SCGFs of the first passage time statistics are derived by finding the largest modulus of the roots of this equation with respect to the z-transform parameter and then calculating its logarithm. The procedure is analogous to that of solving for the SCGFs of the counting statistics. However, for current-like variables, the method generally fails unless the equation of pole is simplified to a quadratic form. The fundamental reason for this lies in the nonuniqueness between the roots and the region of convergence for the joint transform. We illustrate these results via a resonantly driven two-level quantum system, where for several counting variables the solutions to the SCGFs of the first passage time are analytically obtained. Furthermore, we apply these functions to investigate quantum violations of the classical kinetic and thermodynamic uncertainty relations.
We use a semi-Markov process method to calculate large deviations of counting statistics for three open quantum systems, including a resonant two-level system and resonant three-level systems in the Λ and V configurations. In the first two systems, radical solutions to the scaled cumulant generating functions are obtained. Although this is impossible in the third system, since a general sixth-degree polynomial equation is present, we still obtain asymptotically large deviations of the complex system. Our results show that, in these open quantum systems, the large deviation rate functions at zero current are equal to two times the largest nonzero real parts of the eigenvalues of operator -iH[over ̂], where H[over ̂] is a non-Hermitian Hamiltonian, while at a large current these functions possess a unified formula.
A semi-Markov process method for obtaining general counting statistics for open quantum systems is extended to the scenario of resetting. The simultaneous presence of random resets and wave function collapses means that the quantum jump trajectories are no longer semi-Markov. However, focusing on trajectories and using simple probability formulas, general counting statistics can still be constructed from reset-free statistics. An exact tilted matrix equation is also obtained. The inputs of these methods are the survival distributions and waiting-time density distributions instead of quantum operators. In addition, a continuous-time cloning algorithm is introduced to simulate the large-deviation properties of open quantum systems. Several quantum optics systems are used to demonstrate these results.
It is well known that the state operator of an open quantum system can be generically represented as the solution of a time-local equation — a quantum master equation. Unraveling in quantum trajectories offers a picture of open system dynamics dual to solving master equations. In the unraveling picture, physical indicators are computed as Monte Carlo averages over a stochastic process valued in the Hilbert space of the system. This approach is particularly adapted to simulate systems in large Hilbert spaces. We show that the dynamics of an open quantum system generically admits an unraveling in the Hilbert space of the system described by a Markov process generated by ordinary stochastic differential equations for which rigorous concentration estimates are available. The unraveling can be equivalently formulated in terms of norm-preserving state vectors or in terms of linear “ostensible” processes trace preserving only on average. We illustrate the results in the case of a two level system in a simple boson environment. Next, we derive the state-of-the-art form of the Diósi-Gisin-Strunz Gaussian random ostensible state equation in the context of a model problem. This equation provides an exact unraveling of open systems in Gaussian environments. We compare and contrast the two unravelings and their potential for applications to quantum error mitigation.
We analyze dynamical large deviations of quantum trajectories in Markovian open quantum systems in their full generality. We derive a quantum level-2.5 large deviation principle for these systems, which describes the joint fluctuations of time-averaged quantum jump rates and of the time-averaged quantum state for long times. Like its level-2.5 counterpart for classical continuous-time Markov chains (which it contains as a special case), this description is both explicit and complete, as the statistics of arbitrary time-extensive dynamical observables can be obtained by contraction from the explicit level-2.5 rate functional we derive. Our approach uses an unraveled representation of the quantum dynamics which allows these statistics to be obtained by analyzing a classical stochastic process in the space of pure states. For quantum reset processes we show that the unraveled dynamics is semi-Markovian and derive bounds on the asymptotic variance of the number of quantum jumps which generalize classical thermodynamic uncertainty relations. We finish by discussing how our level-2.5 approach can be used to study large deviations of nonlinear functions of the state, such as measures of entanglement.
Open quantum systems have drawn attention over decades due to its applicability in the foundation of theoretical physics, e.g. statistical mechanics, quantum mechanics and condensed matter physics. The dynamics of open quantum systems has been described as separate entities from their surrounding environment that consists of a very large number of modes, somehow coupled to the mode of the system. Even though the exact solution of the dynamical behavior of the system is impossible to calculate, we obtain a tentative solution using the crucial Markov approximation. The input-output formalism of the quantum Langevin equation (QLE) has been considered as a useful tool which provides a semi-classical description of the dynamics of the system, whereas the master equation provides a complete picture of the dynamics of the system expressed in terms of density matrix. While studying the dynamics of nonlinear system/environment coupling using QLE, we see, for a small value of external field, that the steady state system field does not change much from the steady state field obtained in the absence of nonlinear dissipation. However, in a case where a stronger external field is applied, we see that the deviation becomes substantial from the solution of linear system. We also see that the nonlinear coupling introduces significant difference in the cavity fluctuation spectrum. The description, therefore, provides a potential explanation of parametric effects in terms of nonlinear dissipation phenomena associated with the nonlinear coupling. Even though the theories developed in the context of open quantum systems have proven to be powerful tools, they do not provide a satisfactory platform to be implemented on non-linear Hamiltonians. We often approximate it by linearizing over nonlinear steady state field amplitude, and therefore, the interesting effects are often overlooked. The limitation of the analytics provokes us to simulate open quantum dynamics numerically. The numerical method consists of transformation of the environmental degrees of freedom to a one-dimensional many-body chain, and the computational technique includes numerical diagonalization and renormalization process. The time-adaptive density matrix renormalisation group (t-DMRG) is known as one of the most powerful techniques for the simulation of strongly-correlated many-body quantum systems. In this thesis, along with the theoretical modeling, we implement DMRG numerical scheme for the simulation of canonical S/B model by mapping it to one-dimensional harmonic chain with nearest neighbor interactions, and use the method to investigate the dynamics of the free dissipative system. The thermalization of open quantum systems is also studied by generating minimally entangled typical thermal states (METTS) through imaginary time evolution, and real-time evolving an empty system in the presence of the thermal bath. Further, we simulate coherently driven free dissipative Kerr nonlinear system numerically using Euler's method by solving Heisenberg equation of motion and t-DMRG algorithm, and demonstrate how the numerical results are analogous to classical bistability. By comparing with analytics, we see that the DMRG numerics is analogous to the quantum-mechanical exact solution obtained by mapping the equation of motion of the density matrix of the system to a Fokker-Plank equation. The comparison between two different numerical techniques shows that the semi-classical Euler's method determines the dynamics of the system field of one among two coherent branches, whereas DMRG numerics gives the superposition of both of them. Hence, DMRG-determined time dynamics undergoes generating non-classical states. Our approach of dealing with nonlinearity represents an important contribution in the developments of technique to study the dynamical and steady-state behavior of open quantum systems, which is a fundamental aspect of quantum physics.