Computing Large Deviations of First-Passage-Time Statistics in Open Quantum Systems: Two Methods

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 is then obtained as the logarithm of the boundary values within this region. The theoretical basis lies in the facts that the dynamics of the open quantum systems can be unraveled into a piecewise deterministic process and there exists a tilted Liouville master equation in Hilbert space. The second method uses a simulation-based approach 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 alongside cloning simulations for a field-driven system comprising two interacting two-level atoms.