We consider the estimation of parameters encoded in the measurement record of a continuously monitored quantum system in the jump unraveling, corresponding to a single-shot scenario, where information is continuously gathered. Here, it is generally difficult to assess the precision of the estimation procedure via the Fisher Information due to intricate temporal correlations and memory effects. In this paper we provide a full set of solutions to this problem. First, for multi-channel renewal processes we relate the Fisher Information to an underlying Markov chain and derive a easily computable expression for it. For non-renewal processes, we introduce a new algorithm that combines two methods: the monitoring operator method for metrology and the Gillespie algorithm which allows for efficient sampling of a stochastic form of the Fisher Information along individual quantum trajectories. We show that this stochastic Fisher Information satisfies useful properties related to estimation on a single run. Finally, we consider the case where some information is lost in data compression/post-selection and provide tools for computing the Fisher Information in this case. All scenarios are illustrated with instructive examples from quantum optics and condensed matter.
We consider the estimation of parameters encoded in the measurement record of a continuously monitored quantum system in the jump unraveling, corresponding to a single-shot scenario, where information is continuously gathered. Here, it is generally difficult to assess the precision of the estimation procedure via the Fisher Information due to intricate temporal correlations and memory effects. In this paper we provide a full set of solutions to this problem. First, for multi-channel renewal processes we relate the Fisher Information to an underlying Markov chain and derive a easily computable expression for it. For non-renewal processes, we introduce a new algorithm that combines two methods: the monitoring operator method for metrology and the Gillespie algorithm which allows for efficient sampling of a stochastic form of the Fisher Information along individual quantum trajectories. We show that this stochastic Fisher Information satisfies useful properties related to estimation on a single run. Finally, we consider the case where some information is lost in data compression/post-selection and provide tools for computing the Fisher Information in this case. All scenarios are illustrated with instructive examples from quantum optics and condensed matter.
Continuous monitoring is a powerful tool for analyzing quantum systems and is increasingly used as a non-demolition technique in quantum metrology. In this context, noisy data acquired through continuous monitoring can be used to precisely pinpoint unknown parameters of the system. However, extending the theoretical framework that connects these data to the underlying parameters beyond Markovian dynamics is notoriously difficult, primarily because such dynamics cannot be expressed as completely positive divisible maps. To overcome this, we propose a method based on the reaction coordinate mapping to extend these parameter-estimation techniques beyond the Markovian regime. Our approach specifically targets linear systems with non-Markovian dynamics undergoing Gaussian continuous measurements, such as homodyne detection. Within this framework, we analyze Bayesian estimation and provide an analytical expression for the Fisher information, alongside the asymptotic scaling of the estimation precision for arbitrary parameters. Finally, we demonstrate the efficacy of our method through the example of thermometry of a bosonic bath.
The jump unraveling of a quantum master equation decomposes the dynamics of an open quantum system into abrupt jumps, interspersed by periods of coherent dynamics when no jumps occur. Such open quantum systems are ubiquitous in quantum optics and mesoscopic physics, hence the need for efficient techniques for their stochastic simulation. Numerical simulation techniques fall into two main categories. The first splits the evolution into small time steps and determines stochastically for each step if a jump occurs or not. The second, known as Monte Carlo wave-function simulation, is based on the reduction of the norm of an initially pure state in the conditional no-jump evolution. It exploits the fact that the purity of the state is preserved by the finest unraveling of the master equation. In this paper, we present an alternative method for the simulation of the quantum jump unraveling, inspired by the classical Gillespie algorithm. The method is particularly well suited for situations in which a large number of trajectories is required for relatively small systems. It allows for non-purity-preserving dynamics, such as the ones generated by partial monitoring and channel merging. We describe the algorithm in detail and discuss relevant limiting cases. To illustrate it, we include four example applications of increasing physical complexity and discuss the performance of the algorithm across regimes of interest for open quantum systems simulation. Publicly available implementations of our code are provided in julia and mathematica.
Stochastic unravelings provide a useful way to represent open quantum system dynamics in terms of pure state realizations, and have been widely studied both from a fundamental and from a computational point of view. They were initially formulated for Markovian dynamics described by the Gorini-Kossakowski-Sudarshan-Lindblad master equation. However, due to recent technological and experimental development, most physical relevant dynamics present temporal correlations beyond the Markov approximation. Such correlations cause decay rates to turn temporarily negative, thus requiring the generalization of stochastic unravelings from Markovian to non-Markovian scenarios. Indeed, many unraveling techniques have been introduced in this regime, and a comprehensive review of the different jump methods is currently missing. In this work, we provide an overview of widely used quantum jump unraveling techniques for non-Markovian systems and also discuss them in terms of their numerical efficiency, divisibility requirements, Hilbert space extension, and measurement interpretation.
Abstract Many real-world tasks include some kind of parameter estimation, i.e. the determination of a parameter encoded in a probability distribution. Often, such probability distributions arise from stochastic processes. For a stationary stochastic process with temporal correlations, the random variables that constitute it are identically distributed but not independent. This is the case, for instance, for quantum continuous measurements. In this article, we derive the asymptotic Fisher information rate for a stationary process with finite Markov order. We give a precise expression for this rate which is determined by the process’ conditional distribution up to its Markov order. Second, we demonstrate with suitable examples that correlations may both enhance or hamper the metrological precision. Indeed, unlike for entropic information quantities, in general nothing can be said about the sub- or super-additivity of the joint Fisher information in the presence of correlations. To illustrate our results, we apply them to thermometry on an Ising spin chain, considering nearest-neighbour and next-to-nearest neighbour coupling. In this case, the asymptotic Fisher information rate is directly connected to the specific heat capacity of the spin chain. We observe that the presence of correlations strongly enhances the estimation precision in an anti-ferromagnetic chain, while in a ferromagnetic chain this is not the case.
As we enter the era of quantum technologies, quantum estimation theory provides an operationally motivating framework for determining high precision devices in modern technological applications. The aim of any estimation process is to extract information from an unknown parameter embedded in a physical system such as the estimation converges to the true value of the parameter. According to the Cramér-Rao inequality in mathematical statistics, the Fisher information in the case of single-parameter estimation procedures, and the Fisher information matrix in the case of multi-parameter estimation, are the key quantities representing the ultimate precision of the parameters specifying a given statistical model. In quantum estimation strategies, it is usually difficult to derive the analytical expressions of such quantities in a given quantum state. This review provides comprehensive techniques on the analytical calculation of the quantum Fisher information as well as the quantum Fisher information matrix in various scenarios and via several methods. Furthermore, it provides a mathematical transition from classical to quantum estimation theory applied to many freedom quantum systems. To clarify these results, we examine these developments using some examples. Other challenges, including their links to quantum correlations and saturating the quantum Cramér-Rao bound, are also addressed.
We study the problem of parameter estimation in time series stemming from general stochastic processes, where the outcomes may exhibit arbitrary temporal correlations. In particular, we address the question of how much Fisher information is lost if the stochastic process is compressed into a single histogram, known as the empirical distribution. As we show, the answer is nontrivial due to the correlations between outcomes. We derive practical formulas for the resulting Fisher information for various scenarios, from generic stationary processes to discrete-time Markov chains to continuous-time classical master equations. The results are illustrated with several examples.