1. Introduction
The phenomenology of the quantum-to-classical transition, which is the process that drives an otherwise quantum system toward a fully classical description of its physical configuration, is the object of an extensive body of research. Indeed, whether such a transition is due to new fundamental physics or not is a controversial matter [
1]. In particular, it is still under debate if the decoherence of a quantum system that grows in complexity and size can be ascribed to an intrinsic mechanism or only to the unavoidable presence of the surrounding environment [
2,
3].
Motivated by the fact that environmental decoherence does not provide a satisfactory solution to the measurement problem, and thus to the quantum-to-classical transition issue, collapse models embody an alternative theoretical framework [
4,
5]. By elevating the collapse of the wavefunction to a universal physical mechanism embedded in stochastic dynamics, collapse models explain the quantum-to-classical transition in a phenomenological fashion, thus embodying an instance of macrorealistic modifications of quantum mechanics.
Such modification is achieved through a stochastic Schödinger equation and the introduction of new fundamental parameters. When used to assess the dynamics of microscopic systems, the framework of collapse models recovers standard quantum mechanics. Moving toward larger systems, coherence is rapidly suppressed to prevent large spatial superpositions of macroscopically distinguishable states.
The continuous spontaneous localization (CSL) is one of the most well-studied collapse models [
6,
7]. It describes the loss of coherence in the position basis by way of an extra dissipative term entering the master equation of a quantum system. This means that an open quantum system subjected to the collapse mechanism should experience additional dissipation not ascribable to any of the other environmental noise sources. Testing this model is of current interest for the exploration of the limits of validity quantum mechanics [
8]. However, the predicted collapse effect for most of the systems currently used in quantum labs is very weak and thus challenging to detect and distinguish from other environmental noise effects.
Entering the regime where the collapse mechanism is dominant is a tall order. It requires extremely good isolation from environmental noises and ultra-sensitive devices. On the other hand, studies of statistical inference techniques, such as hypothesis testing and parameter estimation, can be employed to ease the requirements and smooth the path toward experimental tests [
9,
10,
11,
12]. Indeed, in Ref. [
10], we showed that a quantum hypothesis testing protocol applied to a mesoscopic optomechanical system—whose massive mechanical mode would be subjected to the CSL mechanism, if any—provides advantages with respect to comparable classical strategies and during the transient dynamics before the onset of a stationary state. In this work, we employ the same optomechanical set-up and look at the problem of parameter estimation rather than noise-source discrimination. The parameter of interest will be the CSL diffusion rate
encoding the two free phenomenological parameters of the collapse model (cf.
Section 3).
The remainder of this work is organized as follows. In
Section 2, we briefly recall elements of quantum parameter estimation theory relevant to our investigation. In
Section 3, we describe shortly the optomechanical set-up of interest. This is a two-cavity system for which quantum advantages stemming from the application of hypothesis testing techniques to collapse model dynamics have been previously shown [
10]. In
Section 4, we perform a dynamical analysis showing that an advantage, analogous to the one for the hypothesis testing, emerges during the transient also for quantum parameter estimation. In
Section 5, we then analyze the steady state of our optomechanical set-up and show that classical measurement strategies and input noise outperform the quantum strategies considered. We conclude in
Section 6 with a discussion of our findings.
2. Parameter Estimation Theory
In Ref. [
10], we considered the advantage that arises in a quantum hypothesis testing scenario aiming to test the presence of a collapse mechanism. While a hypothesis testing protocol allows us to determine, up to a certain confidence level, whether something akin to a collapse mechanism is acting upon a system [
13], one can further wonder how well the collapse parameters can be estimated in principle and which measurement strategies offer the best chances. Addressing these issues involves the use of quantum estimation theory tools.
In (local) quantum estimation theory [
14,
15], the quantum Cramer–Ráo bound [
16] defines the ultimate precision limit for the estimation of a parameter encoded in the state of the system. Indeed, in general, the parameter of interest (
) does not correspond to a directly measurable observable of the system, and its estimation has to be made indirectly, via the measurement of an observable of the parameter-dependent state
of the system. In classical estimation theory, the Fisher information
provides a lower bound to the mean square error of any estimator of the parameter
—in the following, we will consider unbiased estimators for which the mean square error coincides with the variance [
17]. This is known as the (classical) Cramer–Ráo bound and reads
where
n is the number of measurements. The Fisher information is defined as
where
is the conditional probability of obtaining the outcome
x when the parameter has value
. It is important to note that this quantity, and thus the classical Cramer–Ráo bound, depends on the measurement strategy that is adopted to extract information from the state of the system. This is encoded in the conditional probabilities that can be recast in the form of the Born rule
with
defining the POVM corresponding to the measurement strategy.
The ultimate bound to the precision for the estimation of a parameter can then be achieved by optimizing over all possible generalized measurement schemes. This optimization defines the quantum Fisher information [
18,
19]
which, in turn, gives us the aforementioned quantum Cramer–Ráo bound
. Here, we introduced the symmetric logarithmic derivative
, defined by
.
In the following, we will focus on Gaussian quantum systems [
20,
21]. For this particular class of systems, it is convenient to use a phase-space formalism that focuses solely on the first and second moments of the quadratures of the system [
9]. The latter is compactly represented by the covariance matrix
, whose elements are given by
in terms of the components of the quadrature operators vector
. Moreover, we will also restrict our considerations to the local Gaussian measurement on a single Gaussian mode. These are represented by a spectralization of the identity in terms of single-mode Gaussian states characterized by their covariance matrix
. Following [
9], the classical Fisher information can be written in this case as
where
is the sum of the covariance matrix of the state of the system
and the one that characterizes the measurement POVM. According to our previous discussion, through the Cramer–Ráo bound, this quantity bounds the precision on the estimation of the parameter
by any single mode Gaussian measurement characterized by the covariance matrix
. Analogously, the quantum Fisher information for single-mode Gaussian states can be written as [
9]
3. System and Collapse Mechanism
In the rest of this work, we will investigate the precision limit for parameter estimation in a specific optomechanical set-up affected by a collapse mechanism. We consider the two-cavity system shown in
Figure 1.
This same set-up was recently analyzed for the purpose to show transient quantum advantages in quantum hypothesis testing for collapse models [
10], and it is inspired by the quantum reading scheme in [
22]. The system consists of an optomechanical cavity and a second, normal cavity used as an auxiliary system. The optomechanical cavity is initially pumped with coherent light until it reaches its steady state. Afterward, an extra laser beam is injected in both cavities, and the output modes are measured either directly or after recombination through a beam-splitter. We study the driven dynamics by comparing classical and quantum sources of input light and local or EPR-like measurements of the output modes. We assume that the driving field is strong enough to allow for the linearization of the dynamics. This means that any operators can be split into two parts
, where
is a mean-field part that behaves classically, and
is the quantum fluctuation part. We are going to study the dynamics of the quantum fluctuations and, for simplicity, we will suppress the
when indicating the quantum fluctuations’ quadrature from here on.
In the linear approximation, the Hamiltonian of the system is at most quadratic in the quadratures of the system, and the noise we consider has at most linear jump operators. Thus, the open system dynamic is Gaussian [
23]. In order to use the Gaussian formalism for the parameter estimation framework, we also restrict our analysis to Gaussian measurements. This simply requires addressing the covariance matrix of the CSL-affected optomechanical system
, whose elements
are obtained from the zero-mean quantum fluctuations vector
. Here, the first two components
are the dimensionless quadratures for the mechanical mode of the optomechanical cavity, which is modeled as a harmonic oscillator with frequency
and damping rate
. The remaining quadratures,
, are the optical modes for both cavities
. The time evolution then is given by the Lyapunov-like equation
where
A is the drift matrix depending on the physical parameters of the system, and
D is the diffusion matrix. The latter accounts for the noises as
, where the quantum noise operators vector is
, where
is the cavity decay rate, assumed to be the same for both cavity for simplicity [
24].
Here, we consider the following noise sources:
The extra input light fields, given by the operators for each cavity .
The Brownian noise, described by the noise operator , characterized by the Markovian correlation functions . Here, is the Boltzmann constant, while T is the temperature of the surrounding thermal environment.
The CSL collapse model described by , acting as an extra source of decoherence.
The decoherence due to the collapse mechanism can be effectively ascribed to a stochastic force [
25], characterized by the two-point correlation function
. The associated diffusion rate
depends on the two fundamental CSL parameters, the rate of collapse
, and the decoherence length
; it involves also a mass-scaling factor
which can be written as
where
amu (atomic mass unit) and
is the Fourier transform of the mass density of the system subject to the CSL. In the following, we will concentrate on the quantum estimation of the diffusion parameter
affecting the dynamics of the mechanical mode in our set-up. It should be noticed that the two phenomenological parameters entering in the definition of
, i.e.,
, constitute the parameter space of the CSL collapse model. Considerable effort has been spent in the past few years to experimentally rule out portions of this parameter space; for a review of these efforts, we refer the reader to [
5,
8,
26,
27] and the references therein. Since the aim of this work is to study the impact of quantum resources and measurement schemes in the quantum estimation of the CSL parameter
, for our numerical results, we will set the value of
. This is consistent with the values used in Ref. [
10] and results from considering realistic values of
, for a spherical micrometer oscillator, and Adler’s collapse rate [
28]
.
The observable consequence of the collapse mechanism is an overheating of the system, mathematically represented by the additional contribution
to the stochastic Brownian force
. In the diffusion matrix, the collapse diffusion rate
enters the mechanical mode as an extra thermal constant, added to the Brownian contribution
where
is the 4 × 4 covariance matrix associated to the driving light input modes
, and
is a
matrix of zeroes. Here,
is the unknown parameter at the center of our parameter estimation effort.
In our set-up, the dynamics of the mechanical system, affected by the collapse mechanism, is indirectly monitored by measuring the cavities’ output modes. As already discussed, we restrict the detection of the optical output modes to local Gaussian POVM measurements characterized by the single-mode Gaussian states covariance matrix
. Here,
parametrizes the degree of squeezing of the POVM, i.e.,
corresponds to homodyne detection and
heterodyne detection, and the matrix
describes a rotation in phase-space in terms of the Pauli matrix
with
determining the direction along which the measurement is performed. Thus, the total covariance matrix entering the definition of the classical Fisher information (
4) is given by
Here,
is a
diagonal block of the evolved
covariance matrix obtained as a solution of Equation (
6) and it pertains to a single optical cavity mode—either the one of the first cavity or a linear combination of the two optical modes via a beam splitter as explicated in the following. This will be the only quantity needed to calculate both the classical and the quantum Fisher information. In the latter case, we also just need
.
4. Dynamical Analysis
Let us consider the dynamic evolution of the system before it reaches its steady state. The initial state of the system is chosen to be the product of the steady states obtained when only coherent light is pumped into the cavities. Thus, the optomechanical cavity will be in a steady state of the light field and the mechanical element, while the second cavity will simply be in its ground state. Once this initial state is reached, it is possible to drive the system by using additional laser light and the output modes of both cavities can be measured.
We compared two strategies, that we call
classical and
quantum according to the choice of input resources and type of measurement performed [
10]. The
classical strategy involves two independent thermal input noises as classical sources driving the dynamics. This is combined with a local measurement of the optical field of the first cavity
. Note that here, we refrain from explicitly accounting for the input–output relations needed when one considers the measurement of the output cavity field. This is a reasonable assumption, given the linearity of the input–output relations, that allows performing measurements on the internal cavity modes without interfering with the output modes [
29].
The initial covariance matrix for thermal states with a mean number of photons
and
, respectively, reads
The local measurement is performed on the first cavity optical mode and thus it concerns the
central diagonal block of the full
covariance matrix solution of Equation (
6).
The
quantum strategy, on the other hand, makes use of a two-mode squeezed (TMS) light field as correlated input noise and a quantum measurement of EPR-type quadratures obtained by combining the optical fields of the two cavities with a beam-splitter. TMS states are Gaussian states whose covariance matrix, entering Equation (
9), depends only on the squeezing amplitude
and the squeezing angle
and can be written as
where
The EPR-like measurements correspond to measuring a linear combination of the optical modes of the cavities obtained via a 50:50 beam-splitter, giving rise to modes with quadratures
In terms of covariance matrix elements, it means that the
submatrix
of the solution to Equation (
6), representing the covariance matrix of the two optical cavity modes at a generic time
t, has to go through a simplectic transformation describing the modes recombination via the beam-splitter [
23]
The operator
satisfies the equation
, where
is the symplectic matrix and
the beam-splitter Hamiltonian with
as the beam-splitter angle. Here,
(
are the creation and annihilation operators for the two optical cavity modes, respectively [
30].
These two measurement strategies are analogous to the ones used in the quantum hypothesis testing employing the same optomechanical set-up in [
10]. In this case, however, since we are interested in the precision limit to the estimation of the CLS parameter, we look at the (classical) Fisher information for the two strategies that we discussed. The classical Fisher information is obtained from Equation (
4). While
characterizes the measurement on the single optical mode,
is the CSL-affected covariance matrix obtained by (i) solving Equation (
6) with either the input noise from Equation (
11), for the classical strategy, or Equation (
12), for the quantum strategy, and (ii) either focusing on the
submatrix corresponding to the optical mode of the first cavity, for the classical strategy, or one the
covariance matrix of one of the optical modes emerging from the beam-splitter mixing the optical modes of the two cavities, for the quantum strategy.
In
Figure 2, we show the classical Fisher information for the two strategies in function of time. We observe that for early times, the
quantum scheme gives a higher value of the Fisher information than the
classical one. This translates in a lower bound, with respect to the classical scheme, on the precision of the estimation of the CSL parameter
. The precision to which we can estimate the parameter
using classical input states and measurements can be overcome at short times by using non-classical resources, namely TMS states and EPR measurements. We also observe that the quantum advantage is lost at later times. This is expected for systems subjected to decoherence arising from thermal noises [
31]. These results are in agreement with those obtained in [
10], where, considering the same set-up with the same choice of parameters, a
quantum advantage at short times was proven for (quantum) hypothesis testing aimed at certifying the presence of the CLS collapse mechanism.
5. Steady-State Analysis
Having considered the transient dynamics, in this section, we perform a steady-state analysis. In line with the previous discussion, we observe that the best performance is always obtained with a classical scheme.
Once the full system reaches a steady state, all memory about the initial state is lost. However, according to the kind of input noises we subjected the system to—either thermal or TMS light—the dynamics will drive the systems to different steady states. We compute both the classical and the quantum Fisher information at the steady state using the covariance matrix
obtained as the solution of Equation (
6) when setting the right-hand side to zero, i.e.,
In
Figure 3, we show the classical Fisher information at the steady-state reached with TMS-input-noise driven dynamics. The plots are in function of the squeezing parameter
r of the input light field. We compare the two measurement schemes, the
quantum one, i.e., EPR measurement, and the
classical one, that employs local measurements. Higher values of the Fisher information are obtained for lower values of the squeezing parameter. In particular, the maximum is obtained when
and for local measurements. In other words, when there is no 2-mode squeezing in the input noise and we only focus on the first cavity, we obtain the minimum error in the estimation of the CSL parameter
. This corresponds to the case in which we drive the cavities with just coherent light, which can be considered a classical input light field and, as a matter of fact, we completely neglect the second optical cavity. In particular, we see that neither input-noise 2-mode squeezing nor EPR-like measurements of the optical cavity modes can lead to an advantage in the estimation of the CSL parameter. Therefore, we conclude that the use of quantum measurements and input noise is not helpful in the estimation of the CSL parameter
at the steady state, where instead, local measurements and vacuum input noise lead to the best estimate.