Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
\WarningFilter

revtex4-2Repair the float

Enhancing Gaussian quantum metrology with position-momentum correlations

João C. P. Porto [Uncaptioned image] Departamento de Física, Universidade Federal do Piauí, Campus Ministro Petrônio Portela, CEP 64049-550, Teresina, PI, Brazil    Lucas S. Marinho [Uncaptioned image] lucas.marinho@ufpi.edu.br Departamento de Física, Universidade Federal do Piauí, Campus Ministro Petrônio Portela, CEP 64049-550, Teresina, PI, Brazil    Pedro R. Dieguez[Uncaptioned image] pedro.dieguez@ug.edu.pl International Centre for Theory of Quantum Technologies, University of Gdańsk, Jana Bażyńskiego 1A, 80-309 Gdańsk, Poland    Irismar G. da Paz [Uncaptioned image] Departamento de Física, Universidade Federal do Piauí, Campus Ministro Petrônio Portela, CEP 64049-550, Teresina, PI, Brazil    Carlos H. S. Vieira [Uncaptioned image] carloshsv09@gmail.com Centro de Ciências Naturais e Humanas, Universidade Federal do ABC, Avenida dos Estados 5001, 09210-580 Santo André, São Paulo, Brazil.
(August 23, 2024)
Abstract

Quantum metrology offers significant improvements in several quantum technologies. In this work, we propose a Gaussian quantum metrology protocol assisted by initial position-momentum correlations (PM). We employ a correlated Gaussian wave packet as a probe to examine the dynamics of Quantum Fisher Information (QFI) and purity based on PM correlations to demonstrate how to estimate the PM correlations and, more importantly, to unlock its potential applications such as a resource to enhance quantum thermometry. In the low-temperature regime, we find an improvement in the thermometry of the surrounding environment when the original system exhibits a non-null initial correlation (correlated Gaussian state). In addition, we explore the connection between the loss of purity and the gain in QFI during the process of estimating the effective environment coupling and its effective temperature.

preprint: APS/123-QED

I Introduction

The search for optimal measurement of properties during classical or quantum noisy processes i.e. determining the ultimate precision in which a parameter can be estimated, after eliminating all technical noise, using ideally accurate instruments, and repeatedly preparing the system in the same state, is crucial for the development of quantum technologies [1]. Quantum metrology offers significant applications such as superresolution imaging [2], high-precision clocks [3], estimation of proper times and accelerations in quantum field theory [4], navigation devices [5], magnetometry [6], optical and gravitational-wave interferometry [7, 8], and thermometry [9].

Currently, the capacity to measure low temperatures accurately in quantum systems is important in a wide range of proof of principles experiments and has attracted significant interest recently due to its critical role in optimizing the performance of quantum technologies [10, 11]. For instance, effective two-level atoms can be employed to minimize the undesired disturbance on the sample, and an optimum quantum probe (a small controllable quantum system) can act as a thermometer with maximal thermal sensitivity [12]. The scaling of the temperature estimation precision with the number of quantum probes was investigated [13], and the impact of initial system-environment correlations was examined to enhance the estimation of environment parameters within spin-spin model at low temperatures [14]. Still in the low-temperature limit, as a fundamental aspect of the third law of thermodynamics, the unattainability principle of reaching absolute zero also implies the impossibility of precisely measuring temperatures near absolute zero [15, 16, 17]. A general approach to low temperature quantum thermometry considers restrictions from both the sample and the measurement process [15]. In continuous-variable systems, a non-Markovian quantum thermometer was proposed to measure the temperature of a quantum reservoir, effectively avoiding error divergence in the low-temperature [18]. Additionally, optimizing the interaction time in a bipartite Gaussian state allows for precise estimation of the local temperature of trapped ions [19].

Employing quantum resources such as coherence and entanglement, it is possible to improve, until the Heisenberg limit, the measurement precision beyond the classical shot-noise limit [20, 21, 22]. Recently, it was presented a method for low-temperature measurement that improves thermal range and sensitivity by generating quantum coherence in a thermometer probe [23]. Optimal quantum metrology strategies usually employ entanglement in initialization and readout stages to significantly enhance measurement precision [24, 25]. However, in practice, producing entangled states or measurements with high precision, in high-dimensional or continuous-variable systems is not a trivial task. Relatively small errors are known to ruin the success of such tasks when compared to the optimal classical one, and the utilization of a hybrid quantum-classical approach to automatically optimize the controls has been taken into consideration [26].

In this work, we employ a position-momentum (PM) correlated Gaussian state in quantum metrology. In order to work with it, a real parameter can be controlled such that, when it is null, the state recovers the standard uncorrelated Gaussian wave packet form. PM correlations were originally investigated by considering the quantized operators x^^𝑥\hat{x}over^ start_ARG italic_x end_ARG and p^^𝑝\hat{p}over^ start_ARG italic_p end_ARG, with [x^,p^]=i^𝑥^𝑝𝑖Planck-constant-over-2-pi[\hat{x},\hat{p}]=i\hbar[ over^ start_ARG italic_x end_ARG , over^ start_ARG italic_p end_ARG ] = italic_i roman_ℏ [27]. These correlations generalize Glauber’s coherent states [28, 29] and minimize the Robertson-Schrödinger uncertainty relation. Practically, this parameter arises from atomic beam propagation along a transverse harmonic potential, acting as a thin lens and causing a quadratic phase shift in the initial state [30]. Gaussian correlated packets have applications in various fields, including quantum optics and double-slit matter-wave interferometers [31, 32, 33, 34, 35]. However, practical applications face challenges due to errors in tuning the parameter, caused by non-ideal incoherent sources of matter waves and the focalization method used to produce it [36, 37, 38].

Here, we propose a scheme to estimate PM correlations. Our protocol consists of three stages, initialization, interaction, and estimation. Moreover, we further apply our protocol to estimate the effective environment coupling and its temperature (thermometry), and we show that the correlated Gaussian system can outperform the standard (uncorrelated) one. More specifically, we estimated the effective temperature of a Markovian bath via the environmental coupling with our quantum thermometer based on a single-mode correlated Gaussian system. Such correlations impact the Quantum Fisher Information (QFI) after the evolution of the initial state through a Markovian bath and establish the conditions to improve our quantum thermometer model assisted by PM correlations.

The dynamics of the QFI and purity will be analyzed to better estimate this initial position-momentum correlation using the Classical Fisher Information (CFI). Purity, like other measures such as entanglement [39], quantum discord [40, 41, 42], and quantum coherence [43], is also a resource, where it is quantified by deviations from the maximally mixed state [44, 39, 45, 46, 43, 22, 47, 48]. Purity can be operationally interpreted as the maximum coherence achievable by unitary operations [48]. It bounds the maximum amount of entanglement and quantum discord that unitary operations can generate, acting then as a fundamental resource for quantum information processing [48]. Unlike other measures, purity is easily accessible in experiments, providing experimental bounds to other quantum quantifiers [49, 50, 51]. In the following, we present our main results. First, we describe our general estimation protocol in Sec. II.1. Then, we investigate the noise environmental effect on PM correlations estimation in Sec. II.2. Finally, in Sec. II.3 we analyze the role of PM correlation in improving the thermal sensitivity in the low-temperature regime, and in Sec. III we discuss our results.

II Results

II.1 Estimation protocol

Refer to caption
Fig. 1: Theoretical protocol to describe the estimation of the parameters γ𝛾\gammaitalic_γ (PM correlations) and ΛΛ\Lambdaroman_Λ (effective environmental coupling). The initialization step is the source production of fullerene wave packets with momentum k𝑘\vec{k}over→ start_ARG italic_k end_ARG and initial coherent length 0subscript0\ell_{0}roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. Then, the focalization step takes place, generating the initial position-momentum correlations γ𝛾\gammaitalic_γ (see appendix B). After, the system interacts with the Markovian thermal bath, characterized by the effective environmental coupling ΛΛ\Lambdaroman_Λ. In the readout stage, parameter estimations are done.

Our quantum estimation protocol schematized in Fig. 1 is represented in three parts: initialization, interaction, and estimation. The initialization step is characterized by the generation of fullerene Gaussian wave packets with momentum k𝑘\vec{k}over→ start_ARG italic_k end_ARG and an initial coherent length 0subscript0\ell_{0}roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. In the sequence, the parameterization process of unknown parameters is separated into two parts. The first one, here denoted as focalization, describes the unitary parameterization associated with the encoding of the PM correlation parameter γ𝛾\gammaitalic_γ employing a real-noisy source (details about this process are discussed in Appendix B). The second one, denoted as interaction, represents a non-unitary parameterization for which the probe interacts with a Markovian bath producing an encoding of the effective environment coupling parameter ΛΛ\Lambdaroman_Λ and in turn its temperature T𝑇Titalic_T in the correlated Gaussian state. Finally, the estimation procedure of the unknown parameters is carried out during the estimation stage. This stage includes (I) PM correlations and (II) thermometry estimations.

We consider as the initial state the following correlated Gaussian state of transverse width σ0subscript𝜎0\sigma_{0}italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT

ψ0(xi)=1σπexp[xi22σ02+iγxi22σ02],subscript𝜓0subscript𝑥𝑖1𝜎𝜋subscriptsuperscript𝑥2𝑖2subscriptsuperscript𝜎20𝑖𝛾subscriptsuperscript𝑥2𝑖2subscriptsuperscript𝜎20\psi_{0}(x_{i})=\frac{1}{\sqrt{\sigma\sqrt{\pi}}}\exp\left[-\frac{{x^{2}_{i}}}% {{2\sigma^{2}_{0}}}+\frac{i\gamma x^{2}_{i}}{2\sigma^{2}_{0}}\right],italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) = divide start_ARG 1 end_ARG start_ARG square-root start_ARG italic_σ square-root start_ARG italic_π end_ARG end_ARG end_ARG roman_exp [ - divide start_ARG italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG start_ARG 2 italic_σ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG + divide start_ARG italic_i italic_γ italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG start_ARG 2 italic_σ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ] , (1)

that represents a PM-correlated Gaussian state [52]. The real parameter γ𝛾\gammaitalic_γ ensures that the initial state is correlated  [52, 53]. Considering the the initial state, Eq. (1), the variance in position becomes σxx=σ02/2subscript𝜎𝑥𝑥superscriptsubscript𝜎022\sigma_{xx}=\sigma_{0}^{2}/2italic_σ start_POSTSUBSCRIPT italic_x italic_x end_POSTSUBSCRIPT = italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / 2, whereas the variance in momentum is σpp=(1+γ2)2/2σ02subscript𝜎𝑝𝑝1superscript𝛾2superscriptPlanck-constant-over-2-pi22superscriptsubscript𝜎02\sigma_{pp}=(1+\gamma^{2})\hbar^{2}/2\sigma_{0}^{2}italic_σ start_POSTSUBSCRIPT italic_p italic_p end_POSTSUBSCRIPT = ( 1 + italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) roman_ℏ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / 2 italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT, and the initial correlations PM is σxp=(x^p^+p^x^)/2x^p^=γ/2subscript𝜎𝑥𝑝delimited-⟨⟩^𝑥^𝑝^𝑝^𝑥2delimited-⟨⟩^𝑥delimited-⟨⟩^𝑝Planck-constant-over-2-pi𝛾2\sigma_{xp}=\langle(\hat{x}\hat{p}+\hat{p}\hat{x})/2\rangle-\langle\hat{x}% \rangle\langle\hat{p}\rangle=\hbar\gamma/2italic_σ start_POSTSUBSCRIPT italic_x italic_p end_POSTSUBSCRIPT = ⟨ ( over^ start_ARG italic_x end_ARG over^ start_ARG italic_p end_ARG + over^ start_ARG italic_p end_ARG over^ start_ARG italic_x end_ARG ) / 2 ⟩ - ⟨ over^ start_ARG italic_x end_ARG ⟩ ⟨ over^ start_ARG italic_p end_ARG ⟩ = roman_ℏ italic_γ / 2. Notice that, for γ=0𝛾0\gamma=0italic_γ = 0 we have a simple uncorrelated Gaussian wave packet, i.e., σxp=0subscript𝜎𝑥𝑝0\sigma_{xp}=0italic_σ start_POSTSUBSCRIPT italic_x italic_p end_POSTSUBSCRIPT = 0. Exploring the Pearson correlation coefficient between x^^𝑥\hat{x}over^ start_ARG italic_x end_ARG and p^^𝑝\hat{p}over^ start_ARG italic_p end_ARG, i.e., r=σxp/σxxσpp(1r1)𝑟subscript𝜎𝑥𝑝subscript𝜎𝑥𝑥subscript𝜎𝑝𝑝1𝑟1r=\sigma_{xp}/\sqrt{\sigma_{xx}\sigma_{pp}}\,(-1\leq r\leq 1)italic_r = italic_σ start_POSTSUBSCRIPT italic_x italic_p end_POSTSUBSCRIPT / square-root start_ARG italic_σ start_POSTSUBSCRIPT italic_x italic_x end_POSTSUBSCRIPT italic_σ start_POSTSUBSCRIPT italic_p italic_p end_POSTSUBSCRIPT end_ARG ( - 1 ≤ italic_r ≤ 1 ), the γ𝛾\gammaitalic_γ parameter turns out to be γ=r/1r2(γ)𝛾𝑟1superscript𝑟2𝛾\gamma=r/\sqrt{1-r^{2}}\,(-\infty\leq\gamma\leq\infty)italic_γ = italic_r / square-root start_ARG 1 - italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ( - ∞ ≤ italic_γ ≤ ∞ ), illustrating the physical meaning of the γ𝛾\gammaitalic_γ as a parameter that encoded the initial correlations between position and momentum for the initial state.

In the standard quantum phase estimation protocols, it is usual to estimate the physical parameter that implements a relative phase shift [54, 55]. Here, our main goal will be to explore whether initial PM correlations can work as a resource to improve quantum metrology under noise. To estimate the initial PM correlations, we consider the production of beam particles from a non-ideal or partially incoherent source (see Fig. 1), where each particle has momentum px=kxsubscript𝑝𝑥Planck-constant-over-2-pisubscript𝑘𝑥p_{x}=\hbar k_{x}italic_p start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT = roman_ℏ italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT in the x𝑥xitalic_x-direction within a certain range σkxsubscript𝜎subscript𝑘𝑥\sigma_{k_{x}}italic_σ start_POSTSUBSCRIPT italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT end_POSTSUBSCRIPT of wavenumber kxsubscript𝑘𝑥k_{x}italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT and can be modeled by the following initial density matrix [56, 37, 57]

ργ(x,x0)=𝑑kxp(kx)ψ0kx(x0)ψ0kx(x0),subscript𝜌𝛾𝑥superscriptsubscript𝑥0differential-dsubscript𝑘𝑥𝑝subscript𝑘𝑥subscript𝜓0subscript𝑘𝑥subscript𝑥0superscriptsubscript𝜓0subscript𝑘𝑥superscriptsubscript𝑥0\rho_{\gamma}(x,x_{0}^{\prime})=\int dk_{x}p(k_{x})\psi_{0k_{x}}(x_{0})\psi_{0% k_{x}}^{*}(x_{0}^{\prime}),italic_ρ start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT ( italic_x , italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) = ∫ italic_d italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT italic_p ( italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT ) italic_ψ start_POSTSUBSCRIPT 0 italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) italic_ψ start_POSTSUBSCRIPT 0 italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) , (2)

with ψ0kx(x0)=ψ0(x0)eikxx0subscript𝜓0subscript𝑘𝑥subscript𝑥0subscript𝜓0subscript𝑥0superscript𝑒𝑖subscript𝑘𝑥subscript𝑥0\psi_{0k_{x}}(x_{0})=\psi_{0}(x_{0})e^{ik_{x}x_{0}}italic_ψ start_POSTSUBSCRIPT 0 italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) = italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) italic_e start_POSTSUPERSCRIPT italic_i italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT. In this model, we consider wave effects only in x𝑥xitalic_x-direction since the energy associated with the momentum of the particles in the z𝑧zitalic_z-direction is high enough such that the momentum component pzsubscript𝑝𝑧p_{z}italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT is sharply defined, i.e., Δpzpzmuch-less-thanΔsubscript𝑝𝑧subscript𝑝𝑧\Delta p_{z}\ll p_{z}roman_Δ italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ≪ italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT. Then, the behavior in the z𝑧zitalic_z-direction can be considered classical [37]. The geometry of the collimation device is such that it selects only particles with a transverse wave number kxsubscript𝑘𝑥k_{x}italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT inside a specific range defined by the width σkxsubscript𝜎𝑘𝑥\sigma_{kx}italic_σ start_POSTSUBSCRIPT italic_k italic_x end_POSTSUBSCRIPT, being described by a classical Maxwell–Boltzmann distribution, p(kx)=e(kx2/2σkx2)/2πσkx𝑝subscript𝑘𝑥superscript𝑒superscriptsubscript𝑘𝑥22superscriptsubscript𝜎𝑘𝑥22𝜋subscript𝜎𝑘𝑥p(k_{x})=e^{-(k_{x}^{2}/2\sigma_{kx}^{2})}/\sqrt{2\pi}\sigma_{kx}italic_p ( italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT ) = italic_e start_POSTSUPERSCRIPT - ( italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / 2 italic_σ start_POSTSUBSCRIPT italic_k italic_x end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT / square-root start_ARG 2 italic_π end_ARG italic_σ start_POSTSUBSCRIPT italic_k italic_x end_POSTSUBSCRIPT. After some algebraic modifications, the initial density matrix of the coherent correlated Gaussian wave packet becomes

ργ(x0,x0)=𝒩0exp{𝒜0x020x02+𝒞0x0x0},subscript𝜌𝛾subscript𝑥0superscriptsubscript𝑥0subscript𝒩0subscript𝒜0superscriptsubscript𝑥02subscript0superscriptsubscript𝑥02subscript𝒞0subscript𝑥0superscriptsubscript𝑥0\rho_{\gamma}(x_{0},x_{0}^{\prime})=\mathcal{N}_{0}\exp\left\{-\mathcal{A}_{0}% x_{0}^{2}-\mathcal{B}_{0}x_{0}^{\prime 2}+\mathcal{C}_{0}x_{0}x_{0}^{\prime}% \right\},italic_ρ start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) = caligraphic_N start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT roman_exp { - caligraphic_A start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - caligraphic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ 2 end_POSTSUPERSCRIPT + caligraphic_C start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT } , (3)

where 𝒩0=1/σ0πsubscript𝒩01subscript𝜎0𝜋\mathcal{N}_{0}=1/\sigma_{0}\sqrt{\pi}caligraphic_N start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1 / italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT square-root start_ARG italic_π end_ARG, 𝒜0=1/202+(1iγ)/2σ02subscript𝒜012superscriptsubscript021𝑖𝛾2superscriptsubscript𝜎02\mathcal{A}_{0}=1/2\ell_{0}^{2}+(1-i\gamma)/2\sigma_{0}^{2}caligraphic_A start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1 / 2 roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( 1 - italic_i italic_γ ) / 2 italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT, 0=1/202+(1+iγ)/2σ02subscript012superscriptsubscript021𝑖𝛾2superscriptsubscript𝜎02\mathcal{B}_{0}=1/2\ell_{0}^{2}+(1+i\gamma)/2\sigma_{0}^{2}caligraphic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1 / 2 roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( 1 + italic_i italic_γ ) / 2 italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT, 𝒞0=1/02subscript𝒞01superscriptsubscript02\mathcal{C}_{0}=1/\ell_{0}^{2}caligraphic_C start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1 / roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT, and 0=σkx1subscript0superscriptsubscript𝜎𝑘𝑥1\ell_{0}=\sigma_{kx}^{-1}roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = italic_σ start_POSTSUBSCRIPT italic_k italic_x end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT is related with the coherence level of the initial state. It is important to note that, in the limit σkx0subscript𝜎𝑘𝑥0\sigma_{kx}\rightarrow 0italic_σ start_POSTSUBSCRIPT italic_k italic_x end_POSTSUBSCRIPT → 0 (ideal collimation), one has 0subscript0\ell_{0}\rightarrow\inftyroman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT → ∞ (𝒞00subscript𝒞00\mathcal{C}_{0}\rightarrow 0caligraphic_C start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT → 0), meaning that the initial state (3) is completely coherent or pure. On the other hand, if σkxsubscript𝜎𝑘𝑥\sigma_{kx}\rightarrow\inftyitalic_σ start_POSTSUBSCRIPT italic_k italic_x end_POSTSUBSCRIPT → ∞, then one has 00subscript00\ell_{0}\rightarrow 0roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT → 0 (𝒞0subscript𝒞0\mathcal{C}_{0}\rightarrow\inftycaligraphic_C start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT → ∞), indicating an completely incoherent or mixed state.

According to the protocol depicted in Fig. 1, after the focalization process, the system interacts with an environment that is considered to be a Markovian bath [58] and its evolution is given by

ρΛ(x,x,t)=𝑑x0𝑑x0KΛ(x,x,t;x0,x0,0)ργ(x0,x0),subscript𝜌Λ𝑥superscript𝑥𝑡differential-dsubscript𝑥0differential-dsubscriptsuperscript𝑥0subscript𝐾Λ𝑥superscript𝑥𝑡subscript𝑥0subscriptsuperscript𝑥00subscript𝜌𝛾subscript𝑥0subscriptsuperscript𝑥0\rho_{\Lambda}(x,x^{\prime},t)=\int\int dx_{0}dx^{\prime}_{0}K_{\Lambda}(x,x^{% \prime},t;x_{0},x^{\prime}_{0},0)\rho_{\gamma}(x_{0},x^{\prime}_{0}),italic_ρ start_POSTSUBSCRIPT roman_Λ end_POSTSUBSCRIPT ( italic_x , italic_x start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_t ) = ∫ ∫ italic_d italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_d italic_x start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_K start_POSTSUBSCRIPT roman_Λ end_POSTSUBSCRIPT ( italic_x , italic_x start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_t ; italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_x start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , 0 ) italic_ρ start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_x start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) , (4)

where

KΛ(x,x,t;x0,x0,0)=m2πtexp{im2t[(xx0)2+(xx0)2]Λt3[(xx)2+(x0x0)2+(xx)(x0x0)]},subscript𝐾Λ𝑥superscript𝑥𝑡subscript𝑥0subscriptsuperscript𝑥00𝑚2𝜋Planck-constant-over-2-pi𝑡𝑖𝑚2Planck-constant-over-2-pi𝑡delimited-[]superscript𝑥subscript𝑥02superscriptsuperscript𝑥subscript𝑥02Λ𝑡3delimited-[]superscript𝑥superscript𝑥2superscriptsubscript𝑥0subscriptsuperscript𝑥02𝑥superscript𝑥subscript𝑥0subscriptsuperscript𝑥0\displaystyle K_{\Lambda}(x,x^{\prime},t;x_{0},x^{\prime}_{0},0)=\frac{m}{2\pi% \hbar t}\;\exp{\Big{\{}\frac{im}{2\hbar t}\left[(x-x_{0})^{2}+(x^{\prime}-x_{0% })^{2}\right]-\frac{\Lambda t}{3}\left[(x-x^{\prime})^{2}+(x_{0}-x^{\prime}_{0% })^{2}+(x-x^{\prime})(x_{0}-x^{\prime}_{0})\right]\Big{\}}},italic_K start_POSTSUBSCRIPT roman_Λ end_POSTSUBSCRIPT ( italic_x , italic_x start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_t ; italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_x start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , 0 ) = divide start_ARG italic_m end_ARG start_ARG 2 italic_π roman_ℏ italic_t end_ARG roman_exp { divide start_ARG italic_i italic_m end_ARG start_ARG 2 roman_ℏ italic_t end_ARG [ ( italic_x - italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( italic_x start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT - italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] - divide start_ARG roman_Λ italic_t end_ARG start_ARG 3 end_ARG [ ( italic_x - italic_x start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_x start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( italic_x - italic_x start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) ( italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_x start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ] } , (5)

is the propagator which includes the environmental effect [37, 59]. We consider the propagation of fullerene molecules and the decoherence effect produced by air molecules scattering [37]. In this case, the effective scattering constant is given by Λ=(8/32)2πmair(kBT)3/2Nw2Λ83superscriptPlanck-constant-over-2-pi22𝜋subscript𝑚airsuperscriptsubscript𝑘𝐵𝑇32𝑁superscript𝑤2\Lambda=(8/3\hbar^{2})\sqrt{2\pi m_{\text{air}}}(k_{B}T)^{3/2}Nw^{2}roman_Λ = ( 8 / 3 roman_ℏ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) square-root start_ARG 2 italic_π italic_m start_POSTSUBSCRIPT air end_POSTSUBSCRIPT end_ARG ( italic_k start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT italic_T ) start_POSTSUPERSCRIPT 3 / 2 end_POSTSUPERSCRIPT italic_N italic_w start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT, where N𝑁Nitalic_N is the total number density of the air, mairsubscript𝑚airm_{\text{air}}italic_m start_POSTSUBSCRIPT air end_POSTSUBSCRIPT is the mass of the air molecule, w𝑤witalic_w is the size of the molecule of the quantum system, kBsubscript𝑘𝐵k_{B}italic_k start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT is the Boltzmann constant, and T𝑇Titalic_T is the bath temperature [59]. In the propagator described by Eq. (5), it is assumed that the center-of-mass state of the system remains undisturbed by the scattering events, a valid condition for fullerene molecules (mmairmuch-greater-than𝑚subscript𝑚airm\gg m_{\text{air}}italic_m ≫ italic_m start_POSTSUBSCRIPT air end_POSTSUBSCRIPT[37]. The environmental scattering decoherence is a ubiquitous effect constantly monitoring the position of quantum systems in the cosmos, that can be caused by air molecules, light (optical photons), solar neutrinos, cosmic muons, background radioactivity, and even the universe’s 3 K cosmic background radiation [59]. It is also essential to emphasize that in the derivation of this model, it was also assumed the following conditions [59, 60]: the system and the surroundings are not initially correlated; when the composite object–environment system is translated, the scattering interaction remains unchanged; the rate of scattering is much faster than the characteristic rate of change in the state of the system; and the distribution of momenta of the particles in the environment is isotropic and obeys a Maxwell–Boltzmann distribution. In physical terms, the ΛΛ\Lambdaroman_Λ parameter quantifies the rate at which spatial coherence over a given distance ΔxΔ𝑥\Delta xroman_Δ italic_x is suppressed, motivating the introduction of a decoherence timescale given by [59]

τdec=1Λ(Δx)2.subscript𝜏dec1ΛsuperscriptΔ𝑥2\tau_{\text{dec}}=\frac{1}{\Lambda(\Delta x)^{2}}.italic_τ start_POSTSUBSCRIPT dec end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG roman_Λ ( roman_Δ italic_x ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG . (6)

This equation will be useful for interpreting how PM correlation is related to loss of purity and gain of Fisher information. After integration and manipulating Eq. (4), we get

ρΛ(x,x,t)=𝒩texp{𝒜tx2tx2+𝒞txx},subscript𝜌Λ𝑥superscript𝑥𝑡subscript𝒩𝑡subscript𝒜𝑡superscript𝑥2subscript𝑡superscript𝑥2subscript𝒞𝑡𝑥superscript𝑥\rho_{\Lambda}(x,x^{\prime},t)=\mathcal{N}_{t}\exp\left\{-\mathcal{A}_{t}x^{2}% -\mathcal{B}_{t}x^{\prime 2}+\mathcal{C}_{t}xx^{\prime}\right\},italic_ρ start_POSTSUBSCRIPT roman_Λ end_POSTSUBSCRIPT ( italic_x , italic_x start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_t ) = caligraphic_N start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT roman_exp { - caligraphic_A start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - caligraphic_B start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT ′ 2 end_POSTSUPERSCRIPT + caligraphic_C start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_x italic_x start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT } , (7)

where 𝒜tsubscript𝒜𝑡\mathcal{A}_{t}caligraphic_A start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT, tsubscript𝑡\mathcal{B}_{t}caligraphic_B start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT, 𝒞tsubscript𝒞𝑡\mathcal{C}_{t}caligraphic_C start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT and 𝒩tsubscript𝒩𝑡\mathcal{N}_{t}caligraphic_N start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT are parameters that include the interaction with the Markovian bath (see Appendix C).

Employing Eq. (7), we can cast the purity of the state μ(γ,Λ,t)=tr[ρΛ2(x,x,t)]𝜇𝛾Λ𝑡trdelimited-[]subscriptsuperscript𝜌2Λ𝑥superscript𝑥𝑡\mu(\gamma,\Lambda,t)=\text{tr}\left[\rho^{2}_{\Lambda}(x,x^{\prime},t)\right]italic_μ ( italic_γ , roman_Λ , italic_t ) = tr [ italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_Λ end_POSTSUBSCRIPT ( italic_x , italic_x start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_t ) ] in the concise form

μ(γ,Λ,t)=𝒜t+t𝒞t4𝒜RR+𝒞R2,𝜇𝛾Λ𝑡subscript𝒜𝑡subscript𝑡subscript𝒞𝑡4subscript𝒜𝑅subscript𝑅superscriptsubscript𝒞𝑅2\mu(\gamma,\Lambda,t)=\sqrt{\frac{\mathcal{A}_{t}+\mathcal{B}_{t}-\mathcal{C}_% {t}}{4\mathcal{A}_{R}\mathcal{B}_{R}+\mathcal{C}_{R}^{2}}},italic_μ ( italic_γ , roman_Λ , italic_t ) = square-root start_ARG divide start_ARG caligraphic_A start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT + caligraphic_B start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - caligraphic_C start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_ARG start_ARG 4 caligraphic_A start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT caligraphic_B start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT + caligraphic_C start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG , (8)

where 𝒜R=Re[𝒜t]subscript𝒜𝑅Redelimited-[]subscript𝒜𝑡\mathcal{A}_{R}=\text{Re}[\mathcal{A}_{t}]caligraphic_A start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT = Re [ caligraphic_A start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ], R=Re[𝒜t]subscript𝑅Redelimited-[]subscript𝒜𝑡\mathcal{B}_{R}=\text{Re}[\mathcal{A}_{t}]caligraphic_B start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT = Re [ caligraphic_A start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ] and 𝒞R=Re[𝒜t]subscript𝒞𝑅Redelimited-[]subscript𝒜𝑡\mathcal{C}_{R}=\text{Re}[\mathcal{A}_{t}]caligraphic_C start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT = Re [ caligraphic_A start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ]. This result produces μ=1𝜇1\mu=1italic_μ = 1 when 0subscript0\ell_{0}\rightarrow\inftyroman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT → ∞ (𝒞t0subscript𝒞𝑡0\mathcal{C}_{t}\rightarrow 0caligraphic_C start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT → 0), i.e., for a completely coherent source and no environment effect Λ=0Λ0\Lambda=0roman_Λ = 0. Outside this regime, we always have μ(γ,Λ,t)<1𝜇𝛾Λ𝑡1\mu(\gamma,\Lambda,t)<1italic_μ ( italic_γ , roman_Λ , italic_t ) < 1, typical of a noisy scenario represented by mixed states (details in Appendix C).

In the following subsections, we present the dynamics of the QFI and the purity in each part of the parameterization process illustrated in our estimating protocol in Fig. 1. We begin estimating PM correlations in a scenery where the quantum system unavoidably interacts with its surrounding environments and then we explore how these correlations can enhance noisy quantum metrology.

II.2 PM correlations estimation

From the mixed Gaussian state described in Eq. (7), it is easy to check that the first moments x^delimited-⟨⟩^𝑥\langle\hat{x}\rangle⟨ over^ start_ARG italic_x end_ARG ⟩ and p^delimited-⟨⟩^𝑝\langle\hat{p}\rangle⟨ over^ start_ARG italic_p end_ARG ⟩ are null and the dimensionless second moments are given by

σxx=t2τ02[12+2(τ02t+γ2)2+σ0202+2Λtσ023];σpp=(1+γ2)2+σ0202+2tΛσ02;formulae-sequencesubscript𝜎𝑥𝑥superscript𝑡2superscriptsubscript𝜏02delimited-[]122superscriptsubscript𝜏02𝑡𝛾22superscriptsubscript𝜎02superscriptsubscript022Λ𝑡superscriptsubscript𝜎023subscript𝜎𝑝𝑝1superscript𝛾22superscriptsubscript𝜎02superscriptsubscript022𝑡Λsuperscriptsubscript𝜎02\sigma_{xx}=\frac{t^{2}}{\tau_{0}^{2}}\left[\frac{1}{2}+2\left(\frac{\tau_{0}}% {2t}+\frac{\gamma}{2}\right)^{2}+\frac{\sigma_{0}^{2}}{\ell_{0}^{2}}+\frac{2% \Lambda t\sigma_{0}^{2}}{3}\right];\;\;\;\;\;\;\sigma_{pp}=\frac{(1+\gamma^{2}% )}{2}+\frac{\sigma_{0}^{2}}{\ell_{0}^{2}}+2t\Lambda\sigma_{0}^{2};italic_σ start_POSTSUBSCRIPT italic_x italic_x end_POSTSUBSCRIPT = divide start_ARG italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_τ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG [ divide start_ARG 1 end_ARG start_ARG 2 end_ARG + 2 ( divide start_ARG italic_τ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 2 italic_t end_ARG + divide start_ARG italic_γ end_ARG start_ARG 2 end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + divide start_ARG italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + divide start_ARG 2 roman_Λ italic_t italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 3 end_ARG ] ; italic_σ start_POSTSUBSCRIPT italic_p italic_p end_POSTSUBSCRIPT = divide start_ARG ( 1 + italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG 2 end_ARG + divide start_ARG italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + 2 italic_t roman_Λ italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ; (9)
σxp=γ2+t2τ0(1+γ2)+tσ02τ002+t2Λσ02τ0.subscript𝜎𝑥𝑝𝛾2𝑡2subscript𝜏01superscript𝛾2𝑡superscriptsubscript𝜎02subscript𝜏0superscriptsubscript02superscript𝑡2Λsuperscriptsubscript𝜎02subscript𝜏0\displaystyle\sigma_{xp}=\frac{\gamma}{2}+\frac{t}{2\tau_{0}}(1+\gamma^{2})+% \frac{t\sigma_{0}^{2}}{\tau_{0}\ell_{0}^{2}}+\frac{t^{2}\Lambda\sigma_{0}^{2}}% {\tau_{0}}.italic_σ start_POSTSUBSCRIPT italic_x italic_p end_POSTSUBSCRIPT = divide start_ARG italic_γ end_ARG start_ARG 2 end_ARG + divide start_ARG italic_t end_ARG start_ARG 2 italic_τ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ( 1 + italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) + divide start_ARG italic_t italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_τ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + divide start_ARG italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_Λ italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_τ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG . (10)

Here, τ0=(mσ02)/subscript𝜏0𝑚superscriptsubscript𝜎02Planck-constant-over-2-pi\tau_{0}=(m\sigma_{0}^{2})/\hbaritalic_τ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = ( italic_m italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) / roman_ℏ is the time at which the distance of the order of the wave packet extension is traversed with a speed corresponding to the dispersion in velocity [61]. Then, the covariance matrix 𝝈𝝈\boldsymbol{\sigma}bold_italic_σ and its inverse 𝝈1superscript𝝈1\boldsymbol{\sigma}^{-1}bold_italic_σ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT can be written as

𝝈(γ,0,Λ,t)=(σxxσxpσpxσpp);𝝈1(γ,0,Λ,t)=μ2adj(𝝈),formulae-sequence𝝈𝛾subscript0Λ𝑡subscript𝜎𝑥𝑥subscript𝜎𝑥𝑝subscript𝜎𝑝𝑥subscript𝜎𝑝𝑝superscript𝝈1𝛾subscript0Λ𝑡superscript𝜇2adj𝝈\boldsymbol{\sigma}(\gamma,\ell_{0},\Lambda,t)=\left(\begin{array}[]{cc}\sigma% _{xx}&\sigma_{xp}\\ \sigma_{px}&\sigma_{pp}\end{array}\right);\;\;\;\;\;\;\;\boldsymbol{\sigma}^{-% 1}(\gamma,\ell_{0},\Lambda,t)=\mu^{2}\text{adj}(\boldsymbol{\sigma}),bold_italic_σ ( italic_γ , roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , roman_Λ , italic_t ) = ( start_ARRAY start_ROW start_CELL italic_σ start_POSTSUBSCRIPT italic_x italic_x end_POSTSUBSCRIPT end_CELL start_CELL italic_σ start_POSTSUBSCRIPT italic_x italic_p end_POSTSUBSCRIPT end_CELL end_ROW start_ROW start_CELL italic_σ start_POSTSUBSCRIPT italic_p italic_x end_POSTSUBSCRIPT end_CELL start_CELL italic_σ start_POSTSUBSCRIPT italic_p italic_p end_POSTSUBSCRIPT end_CELL end_ROW end_ARRAY ) ; bold_italic_σ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ( italic_γ , roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , roman_Λ , italic_t ) = italic_μ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT adj ( bold_italic_σ ) , (11)

where σpx=σxpsubscript𝜎𝑝𝑥subscript𝜎𝑥𝑝\sigma_{px}=\sigma_{xp}italic_σ start_POSTSUBSCRIPT italic_p italic_x end_POSTSUBSCRIPT = italic_σ start_POSTSUBSCRIPT italic_x italic_p end_POSTSUBSCRIPT, adj(𝝈)adj𝝈\text{adj}(\boldsymbol{\sigma})adj ( bold_italic_σ ) is the adjugate matrix of 𝝈𝝈\boldsymbol{\sigma}bold_italic_σ and we have used that, for Gaussian states, det(𝝈)=μ2det𝝈superscript𝜇2\text{det}(\boldsymbol{\sigma})=\mu^{-2}det ( bold_italic_σ ) = italic_μ start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT [62]. Therefore, we can write the Quantum Fisher Information (QFI) for single-mode Gaussian states (see Eq.(24) in Appendix A) as following

Θ Q(Θ)=μ42(1+μ2)Tr{[adj(𝝈)Θ𝝈]2}+2(Θμ)21μ4.superscriptsubscriptΘ QΘsuperscript𝜇421superscript𝜇2Trsuperscriptdelimited-[]adj𝝈subscriptΘ𝝈22superscriptsubscriptΘ𝜇21superscript𝜇4\mathcal{F}_{\Theta}^{\text{\;Q}}(\Theta)=\frac{\mu^{4}}{2(1+\mu^{2})}\text{Tr% }\{[\text{adj}(\boldsymbol{\sigma})\partial_{\Theta}\boldsymbol{\sigma}]^{2}\}% +2\frac{(\partial_{\Theta}\mu)^{2}}{1-\mu^{4}}.caligraphic_F start_POSTSUBSCRIPT roman_Θ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT Q end_POSTSUPERSCRIPT ( roman_Θ ) = divide start_ARG italic_μ start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG start_ARG 2 ( 1 + italic_μ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG Tr { [ adj ( bold_italic_σ ) ∂ start_POSTSUBSCRIPT roman_Θ end_POSTSUBSCRIPT bold_italic_σ ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT } + 2 divide start_ARG ( ∂ start_POSTSUBSCRIPT roman_Θ end_POSTSUBSCRIPT italic_μ ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 1 - italic_μ start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG . (12)

The above expression will be useful for interpreting our results, where we can explicitly understand the role of purity and its derivative as a resource for estimating unknown parameters in the presence of noise. Furthermore, we also investigate the dynamics of Classical Fisher Information (CFI) obtained from the density distribution in position space as

Θ C(Θ)=1P(xi|Θ)(P(xi|Θ)Θ)2𝑑xi,superscriptsubscriptΘ CΘ1𝑃conditionalsubscript𝑥𝑖Θsuperscript𝑃conditionalsubscript𝑥𝑖ΘΘ2differential-dsubscript𝑥𝑖\displaystyle\mathcal{F}_{\Theta}^{\text{\;C}}(\Theta)=\int\frac{1}{P(x_{i}|% \Theta)}\Bigg{(}\frac{\partial P(x_{i}|\Theta)}{\partial\Theta}\Bigg{)}^{2}dx_% {i},caligraphic_F start_POSTSUBSCRIPT roman_Θ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT C end_POSTSUPERSCRIPT ( roman_Θ ) = ∫ divide start_ARG 1 end_ARG start_ARG italic_P ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT | roman_Θ ) end_ARG ( divide start_ARG ∂ italic_P ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT | roman_Θ ) end_ARG start_ARG ∂ roman_Θ end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , (13)

where P(xi|Θ)=ρΛ(x,x=x,t)P(x_{i}|\Theta)=\rho_{\Lambda}(x,x^{\prime}=x,t)italic_P ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT | roman_Θ ) = italic_ρ start_POSTSUBSCRIPT roman_Λ end_POSTSUBSCRIPT ( italic_x , italic_x start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = italic_x , italic_t ). In what follows, the estimated parameters will be the initial position-momentum correlation Θ=γΘ𝛾\Theta=\gammaroman_Θ = italic_γ and the effective environmental coupling constant Θ=ΛΘΛ\Theta=\Lambdaroman_Θ = roman_Λ.

We employ the QFI and CFI for the estimation of the initial correlations. They were calculated from Eqs. (12) and (13) by making Θ=γΘ𝛾\Theta=\gammaroman_Θ = italic_γ. After some simplifications, we obtain

γQ(γ,t,Λ)=μ42(1+μ2)Φγ(γ,t,Λ)+2(γμ)21μ4,superscriptsubscript𝛾Q𝛾𝑡Λsuperscript𝜇421superscript𝜇2subscriptΦ𝛾𝛾𝑡Λ2superscriptsubscript𝛾𝜇21superscript𝜇4\mathcal{F}_{\gamma}^{\;\text{Q}}(\gamma,t,\Lambda)=\frac{\mu^{4}}{2(1+\mu^{2}% )}\;\Phi_{\gamma}(\gamma,t,\Lambda)+2\frac{(\partial_{\gamma}\mu)^{2}}{1-\mu^{% 4}},caligraphic_F start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT Q end_POSTSUPERSCRIPT ( italic_γ , italic_t , roman_Λ ) = divide start_ARG italic_μ start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG start_ARG 2 ( 1 + italic_μ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG roman_Φ start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT ( italic_γ , italic_t , roman_Λ ) + 2 divide start_ARG ( ∂ start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT italic_μ ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 1 - italic_μ start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG , (14)

and

γC(γ,t,Λ)=18σ04B4(mt+γσ02)2,superscriptsubscript𝛾C𝛾𝑡Λ18superscriptsubscript𝜎04superscript𝐵4superscript𝑚Planck-constant-over-2-pi𝑡𝛾superscriptsubscript𝜎022\mathcal{F}_{\gamma}^{\;\text{C}}(\gamma,t,\Lambda)=\frac{1}{8\sigma_{0}^{4}B^% {4}}\left(\frac{m}{\hbar t}+\frac{\gamma}{\sigma_{0}^{2}}\right)^{2},caligraphic_F start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT C end_POSTSUPERSCRIPT ( italic_γ , italic_t , roman_Λ ) = divide start_ARG 1 end_ARG start_ARG 8 italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_B start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG ( divide start_ARG italic_m end_ARG start_ARG roman_ℏ italic_t end_ARG + divide start_ARG italic_γ end_ARG start_ARG italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , (15)
Refer to caption
Fig. 2: Curves for QFI (γ Qsuperscriptsubscript𝛾 Q\mathcal{F}_{\gamma}^{\text{\;Q}}caligraphic_F start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT Q end_POSTSUPERSCRIPT), CFI (γ Csuperscriptsubscript𝛾 C\mathcal{F}_{\gamma}^{\text{\;C}}caligraphic_F start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT C end_POSTSUPERSCRIPT), purity μ𝜇\muitalic_μ, and its absolute relative derivative 1μ|γμ|1𝜇subscript𝛾𝜇\frac{1}{\mu}|\partial_{\gamma}\mu|divide start_ARG 1 end_ARG start_ARG italic_μ end_ARG | ∂ start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT italic_μ | as a function of the initial correlation γ𝛾\gammaitalic_γ for t=1.0𝑡1.0t=1.0italic_t = 1.0 μ𝜇\muitalic_μs, and crescents environmental effect coupling constants: (a) Λ=0Λ0\Lambda=0roman_Λ = 0 (T=0𝑇0T=0italic_T = 0 K), (b) Λ=1020Λsuperscript1020\Lambda=10^{20}roman_Λ = 10 start_POSTSUPERSCRIPT 20 end_POSTSUPERSCRIPT m-2s-1 (T=952𝑇952T=952italic_T = 952 K), (c) Λ=1022Λsuperscript1022\Lambda=10^{22}roman_Λ = 10 start_POSTSUPERSCRIPT 22 end_POSTSUPERSCRIPT m-2s-1 (T=20.5×103𝑇20.5superscript103T=20.5\times 10^{3}italic_T = 20.5 × 10 start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT K), and (d) Λ=1023Λsuperscript1023\Lambda=10^{23}roman_Λ = 10 start_POSTSUPERSCRIPT 23 end_POSTSUPERSCRIPT m-2s-1 (T=95.2×103𝑇95.2superscript103T=95.2\times 10^{3}italic_T = 95.2 × 10 start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT K). In panels (a) and (b), we observe that when the environment has a small effect (low-temperature regime), the QFI exhibits a behavior similar to the purity. In contrast, panels (c) and (d) show that the QFI exhibits a behavior related to the absolute purity variation when the environment has a strong effect (high-temperature regime). The advantage of employing PM correlations can be realized for small values of purity.

where the parameter B𝐵Bitalic_B is explicitly defined in Appendix C and Φγ(γ,t,Λ)subscriptΦ𝛾𝛾𝑡Λ\Phi_{\gamma}(\gamma,t,\Lambda)roman_Φ start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT ( italic_γ , italic_t , roman_Λ ) is a monotonic increasing function of the parameters γ𝛾\gammaitalic_γ, t𝑡titalic_t and ΛΛ\Lambdaroman_Λ (see Appendix D). Since this function is multiplied by the fourth power of μ𝜇\muitalic_μ, the purity behavior determines the first term in QFI. This can be illustrated using the following parameters: fullerene mass m=1.2×1024𝑚1.2superscript1024m=1.2\times 10^{-24}italic_m = 1.2 × 10 start_POSTSUPERSCRIPT - 24 end_POSTSUPERSCRIPT kg, molecular size w=7𝑤7w=7italic_w = 7 Å̊𝐴\mathring{A}over̊ start_ARG italic_A end_ARG, width of the initial wave packet σ0=7.8subscript𝜎07.8\sigma_{0}=7.8italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 7.8 nm, mass of the air molecule mair=5.0×1026subscript𝑚air5.0superscript1026m_{\text{air}}=5.0\times 10^{-26}italic_m start_POSTSUBSCRIPT air end_POSTSUBSCRIPT = 5.0 × 10 start_POSTSUPERSCRIPT - 26 end_POSTSUPERSCRIPT kg and density of the air molecules N=1012𝑁superscript1012N=10^{12}italic_N = 10 start_POSTSUPERSCRIPT 12 end_POSTSUPERSCRIPT molecules/m3 [38, 56, 37]. It is assumed that the collimator apparatus selects wave numbers in the Ox𝑂𝑥Oxitalic_O italic_x direction, with a transverse wavenumber dispersion of σkx107m1subscript𝜎subscript𝑘𝑥superscript107superscriptm1\sigma_{k_{x}}\approx 10^{7}\;\mathrm{m^{-1}}italic_σ start_POSTSUBSCRIPT italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT end_POSTSUBSCRIPT ≈ 10 start_POSTSUPERSCRIPT 7 end_POSTSUPERSCRIPT roman_m start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT. This corresponds to an initial coherence length of 0=σkx150nmsubscript0superscriptsubscript𝜎subscript𝑘𝑥150nm\ell_{0}=\sigma_{k_{x}}^{-1}\approx 50\;\mathrm{nm}roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = italic_σ start_POSTSUBSCRIPT italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ≈ 50 roman_nm. Parameters of this order of magnitude were previously used in experiments with fullerene molecules by Zeilinger [63, 64]. With the other parameters fixed, the variation in ΛΛ\Lambdaroman_Λ corresponds to the variation in the environmental temperature such that, ΛT3/2proportional-toΛsuperscript𝑇32\Lambda\propto T^{3/2}roman_Λ ∝ italic_T start_POSTSUPERSCRIPT 3 / 2 end_POSTSUPERSCRIPT. For instance, ΛΛ\Lambdaroman_Λ of the order of 1019superscript101910^{19}10 start_POSTSUPERSCRIPT 19 end_POSTSUPERSCRIPT m-2s-1, corresponds to an effective temperature around 205 K. The value 3.2×10153.2superscript10153.2\times 10^{15}3.2 × 10 start_POSTSUPERSCRIPT 15 end_POSTSUPERSCRIPT m-2s-1 for the scattering by air molecules estimated in [37] for the experiment with fullerene molecules corresponds to a temperature of 300 K and a density of the air molecules of N=1.8×108𝑁1.8superscript108N=1.8\times 10^{8}italic_N = 1.8 × 10 start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT molecules/m3, smaller than the densities we are considering. In other words, the air molecule density and temperature range we are exploring here are experimentally feasible with the current technology.

To analyze the role of each term in the QFI (12) and the influence of the environmental effect, we present in Fig. 2 the behavior of the QFI (γ Qsuperscriptsubscript𝛾 Q\mathcal{F}_{\gamma}^{\text{\;Q}}caligraphic_F start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT Q end_POSTSUPERSCRIPT), CFI (γ Csuperscriptsubscript𝛾 C\mathcal{F}_{\gamma}^{\text{\;C}}caligraphic_F start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT C end_POSTSUPERSCRIPT), purity μ𝜇\muitalic_μ, and its absolute relative derivative 1μ|γμ|1𝜇subscript𝛾𝜇\frac{1}{\mu}|\partial_{\gamma}\mu|divide start_ARG 1 end_ARG start_ARG italic_μ end_ARG | ∂ start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT italic_μ | as a function of the initial correlation γ𝛾\gammaitalic_γ for an characteristic interaction time t=1.0𝑡1.0t=1.0italic_t = 1.0 μ𝜇\muitalic_μs, and different environmental effect coupling constants. Panels (a) and (b) demonstrate that the QFI behaves similarly to the purity when the environment has a relatively small effect (low-temperature regime). In contrast, panels (c) and (d) demonstrate that the QFI behaves according to the absolute purity variation when the environment effect is stronger. Also, we see that only in the regime of strong environmental effect (quantum to classical transition in the high-temperature regime) does the CFI become comparable with the QFI.

One of the most intriguing features here is that besides the suppression of the QFI when the noise increases, we can find PM correlations that enhance the corresponding QFI compared to the usual Gaussian state. It is also important to note that Fig. 2(d) shows that even a maximal purity value (null PM correlations) does not provide the maximum value for the QFI. This happens because, in the regime of strong environmental effect, the change in purity (second term in Eq. (14)) is dominant over the purity value itself that is close to zero and will be elevated to the fourth power (first term in Eq. (14)). To comprehend in more detail how the QFI behavior is governed by either the purity or by the purity variation, and how this transition occurs, we present in Appendix D these quantities where we fixed ΛΛ\Lambdaroman_Λ and plotted them as a function of the PM correlations and the propagation time.

II.3 Quantum thermometry assisted by PM correlations

Refer to caption
Fig. 3: Curves for QFI (Λ2Λ QsuperscriptΛ2superscriptsubscriptΛ Q\Lambda^{2}\mathcal{F}_{\Lambda}^{\text{\;Q}}roman_Λ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT caligraphic_F start_POSTSUBSCRIPT roman_Λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT Q end_POSTSUPERSCRIPT), CFI (Λ2Λ CsuperscriptΛ2superscriptsubscriptΛ C\Lambda^{2}\mathcal{F}_{\Lambda}^{\text{\;C}}roman_Λ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT caligraphic_F start_POSTSUBSCRIPT roman_Λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT C end_POSTSUPERSCRIPT), purity μ𝜇\muitalic_μ, and its absolute relative variation 1μ|γμ|1𝜇subscript𝛾𝜇\frac{1}{\mu}|\partial_{\gamma}\mu|divide start_ARG 1 end_ARG start_ARG italic_μ end_ARG | ∂ start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT italic_μ | as a function of initial correlation γ𝛾\gammaitalic_γ for t=50.0𝑡50.0t=50.0italic_t = 50.0 μ𝜇\muitalic_μs and considering (a) weak environmental effect (Λ=1.0×1015Λ1.0superscript1015\Lambda=1.0\times 10^{15}roman_Λ = 1.0 × 10 start_POSTSUPERSCRIPT 15 end_POSTSUPERSCRIPT m-2s-1 or T=442𝑇442T=442italic_T = 442 mK) and (b) strong environmental effect (Λ=1.0×1021Λ1.0superscript1021\Lambda=1.0\times 10^{21}roman_Λ = 1.0 × 10 start_POSTSUPERSCRIPT 21 end_POSTSUPERSCRIPT m-2s-1 or T=4.4×103𝑇4.4superscript103T=4.4\times 10^{3}italic_T = 4.4 × 10 start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT K). We observe in panel (b) [(a)] that the gain in Fisher information is related to purity (relative variation in purity) when the environmental effect is strong (weak).

Here, we analyze how the PM correlations affect the estimation of the effective environmental coupling constant, showing a specific regime in which the correlated Gaussian state provides a quantum advantage over the standard one in estimating the effective parameter related to an environmental interaction. Using Eqs. (12) and (13) and setting Θ=ΛΘΛ\Theta=\Lambdaroman_Θ = roman_Λ, the QFI and CFI reads

ΛQ(γ,t,Λ)=μ42(1+μ2)ΦΛ(γ,t,Λ)+2(Λμ)21μ4,superscriptsubscriptΛQ𝛾𝑡Λsuperscript𝜇421superscript𝜇2subscriptΦΛ𝛾𝑡Λ2superscriptsubscriptΛ𝜇21superscript𝜇4\mathcal{F}_{\Lambda}^{\;\text{Q}}(\gamma,t,\Lambda)=\frac{\mu^{4}}{2(1+\mu^{2% })}\;\Phi_{\Lambda}(\gamma,t,\Lambda)+2\frac{(\partial_{\Lambda}\mu)^{2}}{1-% \mu^{4}},caligraphic_F start_POSTSUBSCRIPT roman_Λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT Q end_POSTSUPERSCRIPT ( italic_γ , italic_t , roman_Λ ) = divide start_ARG italic_μ start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG start_ARG 2 ( 1 + italic_μ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG roman_Φ start_POSTSUBSCRIPT roman_Λ end_POSTSUBSCRIPT ( italic_γ , italic_t , roman_Λ ) + 2 divide start_ARG ( ∂ start_POSTSUBSCRIPT roman_Λ end_POSTSUBSCRIPT italic_μ ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 1 - italic_μ start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG , (16)

and

ΛC(γ,Λ,t)=t212σ04B4,superscriptsubscriptΛC𝛾Λ𝑡superscript𝑡212superscriptsubscript𝜎04superscript𝐵4\mathcal{F}_{\Lambda}^{\;\text{C}}(\gamma,\Lambda,t)=\frac{t^{2}}{12\sigma_{0}% ^{4}B^{4}},caligraphic_F start_POSTSUBSCRIPT roman_Λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT C end_POSTSUPERSCRIPT ( italic_γ , roman_Λ , italic_t ) = divide start_ARG italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 12 italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_B start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG , (17)

where ΦΛ(γ,t,Λ)subscriptΦΛ𝛾𝑡Λ\Phi_{\Lambda}(\gamma,t,\Lambda)roman_Φ start_POSTSUBSCRIPT roman_Λ end_POSTSUBSCRIPT ( italic_γ , italic_t , roman_Λ ) is another monotonic increasing function of the parameters γ𝛾\gammaitalic_γ, t𝑡titalic_t and ΛΛ\Lambdaroman_Λ (see Appendix E). Similarly, this function does not determine the profile of QFI in the first term.

In Fig. 3, we exhibit the curves for QFI (Λ2Λ QsuperscriptΛ2superscriptsubscriptΛ Q\Lambda^{2}\mathcal{F}_{\Lambda}^{\text{\;Q}}roman_Λ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT caligraphic_F start_POSTSUBSCRIPT roman_Λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT Q end_POSTSUPERSCRIPT), CFI (Λ2Λ CsuperscriptΛ2superscriptsubscriptΛ C\Lambda^{2}\mathcal{F}_{\Lambda}^{\text{\;C}}roman_Λ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT caligraphic_F start_POSTSUBSCRIPT roman_Λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT C end_POSTSUPERSCRIPT), purity μ𝜇\muitalic_μ, and its absolute relative variation 1μ|γμ|1𝜇subscript𝛾𝜇\frac{1}{\mu}|\partial_{\gamma}\mu|divide start_ARG 1 end_ARG start_ARG italic_μ end_ARG | ∂ start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT italic_μ | as a function of initial correlation γ𝛾\gammaitalic_γ for t=50.0𝑡50.0t=50.0italic_t = 50.0 μ𝜇\muitalic_μs and considering (a) weak environmental effect (Λ=1.0×1015Λ1.0superscript1015\Lambda=1.0\times 10^{15}roman_Λ = 1.0 × 10 start_POSTSUPERSCRIPT 15 end_POSTSUPERSCRIPT m-2s-1 or T=442𝑇442T=442italic_T = 442 mK) and (b) strong environmental effect (Λ=1.0×1021Λ1.0superscript1021\Lambda=1.0\times 10^{21}roman_Λ = 1.0 × 10 start_POSTSUPERSCRIPT 21 end_POSTSUPERSCRIPT m-2s-1 or T=4.4×103𝑇4.4superscript103T=4.4\times 10^{3}italic_T = 4.4 × 10 start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT K). Again, we observe two characteristic regimes, the gain in QFI is related either to purity or relative variation in purity when the environmental effect is strong (high-temperature regime) or weak (low-temperature regime), respectively. Such regimes are associated with the first and second term in Eq. (16), respectively. One of the more interesting results of this work, rarely explored in the literature, is investigating the relation between loss of purity and gain in Fisher information. Also, we can observe that in the low-temperature regime, the PM can improve the QFI, and this behavior paves the way for exploring the enhancement of noisy quantum metrology by position-momentum correlations that will be discussed in the following. On the other hand, these PM correlations quantum advantage ceases to be valid in the limit of strong environmental effect (Λ=1.0×1021Λ1.0superscript1021\Lambda=1.0\times 10^{21}roman_Λ = 1.0 × 10 start_POSTSUPERSCRIPT 21 end_POSTSUPERSCRIPT m-2s-1 or T=4.4×103𝑇4.4superscript103T=4.4\times 10^{3}italic_T = 4.4 × 10 start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT K), with the standard uncorrelated Gaussian state γ=0𝛾0\gamma=0italic_γ = 0 providing the best sensitivity for thermometry in this regime [see Fig. 3 (b)].

Also, it is important to mention that the effect of increasing the propagation time t𝑡titalic_t or changing the environmental effect ΛΛ\Lambdaroman_Λ is almost equivalent from a theoretical point of view, they both decrease the purity and coherence of the quantum state. However, from an experimental point of view, the variation of these two parameters is very different. The first is varied just by moving the position of the detector to regions further away from the source. In contrast, the second can be varied by changing, for example, the environmental temperature or the air pressure.

Refer to caption
Fig. 4: (a) QFI (Λ2ΛQsuperscriptΛ2superscriptsubscriptΛQ\Lambda^{2}\mathcal{F}_{\Lambda}^{\text{Q}}roman_Λ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT caligraphic_F start_POSTSUBSCRIPT roman_Λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT Q end_POSTSUPERSCRIPT) and (b) purity μ𝜇\muitalic_μ as a function of time t𝑡titalic_t and for different values of the initial correlation γ𝛾\gammaitalic_γ and considering Λ=1.0×1015Λ1.0superscript1015\Lambda=1.0\times 10^{15}roman_Λ = 1.0 × 10 start_POSTSUPERSCRIPT 15 end_POSTSUPERSCRIPT m-2s-1 (T=442𝑇442T=442italic_T = 442 mK). We can observe a quantum enhancement due to γ𝛾\gammaitalic_γ. Also, note the relation exists between the time of information saturation and the correspondent maximum temporal rate of purity variation [inset of panel(b)].

In Fig. 4 we investigate the time dependence of (a) QFI (Λ2ΛQsuperscriptΛ2superscriptsubscriptΛQ\Lambda^{2}\mathcal{F}_{\Lambda}^{\text{Q}}roman_Λ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT caligraphic_F start_POSTSUBSCRIPT roman_Λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT Q end_POSTSUPERSCRIPT) and (b) purity μ𝜇\muitalic_μ for different values of the initial correlation γ𝛾\gammaitalic_γ and considering Λ=1.0×1015Λ1.0superscript1015\Lambda=1.0\times 10^{15}roman_Λ = 1.0 × 10 start_POSTSUPERSCRIPT 15 end_POSTSUPERSCRIPT m-2s-1 (T=442𝑇442T=442italic_T = 442 mK). We observe a temporal quantum enhancement due to PM correlation γ𝛾\gammaitalic_γ. Therefore, employing an initial correlated (position-momentum) Gaussian state can reduce the time required to obtain the maximum amount of information until it is saturated. Moreover, note the relationship [inset of panel (b)] between the maximum temporal rate of purity variations and the time of information saturation.

Refer to caption
Fig. 5: Temporal Gain of Information (TGI) as a function of initial PM correlation γ𝛾\gammaitalic_γ. The solid orange line corresponds to the approximate gain calculated from Eq. (42). The full circles correspond to some values of TGI obtained from the exact value that maximizes the function 1μ|tμ|1𝜇subscript𝑡𝜇\frac{1}{\mu}|\partial_{t}\mu|divide start_ARG 1 end_ARG start_ARG italic_μ end_ARG | ∂ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_μ | without approximation in the purity expression. These points are summarized in Tab. 1.

Besides this, to quantify such an enhancement in the time required to saturate the QFI due to PM correlation γ𝛾\gammaitalic_γ, we introduce the Temporal Gain of Information (TGI), defining it on the decibel scale as follows

TGI(γ)=(10dB)log[τmax(Λ,γ)τmax(Λ,γ=0)],TGI𝛾10dBsubscript𝜏maxΛ𝛾subscript𝜏maxΛ𝛾0\text{TGI}(\gamma)=-(10\;\text{dB})\log\Bigg{[}\frac{\tau_{\text{max}}(\Lambda% ,\gamma)}{\tau_{\text{max}}(\Lambda,\gamma=0)}\Bigg{]},TGI ( italic_γ ) = - ( 10 dB ) roman_log [ divide start_ARG italic_τ start_POSTSUBSCRIPT max end_POSTSUBSCRIPT ( roman_Λ , italic_γ ) end_ARG start_ARG italic_τ start_POSTSUBSCRIPT max end_POSTSUBSCRIPT ( roman_Λ , italic_γ = 0 ) end_ARG ] , (18)

where the time τmaxsubscript𝜏max\tau_{\text{max}}italic_τ start_POSTSUBSCRIPT max end_POSTSUBSCRIPT corresponds to the maximization point of the relative purity variation, i.e., the maximum of 1μ|tμ|1𝜇subscript𝑡𝜇\frac{1}{\mu}|\partial_{t}\mu|divide start_ARG 1 end_ARG start_ARG italic_μ end_ARG | ∂ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_μ | [see inset of Fig. 4(b)]. In this definition, the reduction in the time to saturate the QFI is compared with the correspondent time needed for the standard uncorrelated Gaussian state (γ=0𝛾0\gamma=0italic_γ = 0), such that the uncorrelated state is always associated with a null temporal gain. After some approximations (see more details in Appendix F) the TGI can be explicitly given by

TGIapprox(γ)=(10dB)log(1+γ2)1/3.\text{TGI}_{\text{approx}}(\gamma)=(10\;\text{dB})\log(1+\gamma^{2})^{1/3}.TGI start_POSTSUBSCRIPT approx end_POSTSUBSCRIPT ( italic_γ ) = ( 10 dB ) roman_log ( 1 + italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 1 / 3 end_POSTSUPERSCRIPT . (19)

In Fig. 5, we examine the enhancement in TGI as a function of PM correlation γ𝛾\gammaitalic_γ. The solid orange line corresponds to the approximate gain calculated from Eq. (19). The full circles correspond to particular TGI values determined by a precise number that optimizes the function 1μ|tμ|1𝜇subscript𝑡𝜇\frac{1}{\mu}|\partial_{t}\mu|divide start_ARG 1 end_ARG start_ARG italic_μ end_ARG | ∂ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_μ | without any approximation in the purity expression. These points are obtained directly by probe inspection in the inset of Fig 4. (b) and are summarized in Tab. 1. Moreover, this procedure was applied in order to obtain an approximate and simple expression for the TGI and its dependence on the PM correlations γ𝛾\gammaitalic_γ. Note that for PM correlation on the order of 150, the gain achievable is almost 15 dB over the standard uncorrelated Gaussian state. In experimental terms, this parameter can be controlled, for example, by changing the laser frequency or intensity in the initialization step (see Appendix B).

III Discussion

Concepts such as purity distillation and purity cost, analogous to entanglement distillation and entanglement cost, were introduced to quantify this resource [65, 66]. Transformations from pure to mixed states are usually associated with information loss and irreversibility [65]. However, this work explores the relationship between Quantum Fisher Information (QFI) and purity rate, focusing on a regime where purity loss primarily drives QFI behavior. It examines how position-momentum (PM)-correlated probes can enhance noisy quantum metrology of the environmental thermal bath and demonstrates that PM-correlated Gaussian states can outperform standard uncorrelated Gaussian states in estimating the effective environmental coupling constant.

We started presenting the probe initialization procedure and how the PM correlation can be controlled and produced at this stage. We then calculated the dynamics of QFI and purity for the probe as it evolves through a decoherence channel, representing the scattering model, one of the most common sources of decoherence in physical systems [59, 60]. This decoherence model describes how the entanglement of the system with the environmental particles (fermionic or bosonic bath), caused by countless scattering events, delocalizes local phase relations between spatially separated wave-function components, leading to decoherence in position space (i.e., to localization) [59]. Our findings demonstrate the role of purity and purity “velocity” (Θμ)subscriptΘ𝜇(\partial_{\Theta}\mu)( ∂ start_POSTSUBSCRIPT roman_Θ end_POSTSUBSCRIPT italic_μ ) as a resource for estimating unknown parameters in the presence of noise [see Eq. (12)].

We observed that, in estimating the PM correlation, when the decoherence factor is relatively small (i.e., low-temperature regime), the QFI and purity curves have similar behaviors. However, for high temperatures, the QFI is dependent on the rate of change of purity. In this configuration, maximum purity does not necessarily imply maximum values for the QFI due to the dynamics of the purity rate of change term for regimes with an effective strong coupling with the Markovian bath. As expected, QFI tends to match CFI under conditions of quantum-classical transition (high-temperature regime). Furthermore, we found scenarios where the QFI can be asymptotically estimated by the value of CFI by a suitable preparation of the correlated γ𝛾\gammaitalic_γ initial state [see, for example, Fig. 2(a)]. From a practical standpoint, this result is both surprising and useful, since obtaining the QFI is challenging because it requires maximizing all possible POVMs to ensure that estimations are independent of the chosen POVM and to achieve optimal measurement precision.

Regarding the estimation of the effective environmental coupling constant, we observed an improvement in the thermometry of the surrounding environment when the probe was prepared in a non-null initial correlated Gaussian state. Specifically, in this part, we explored the acceleration in the loss of purity, due to correlations, as a resource for increasing the QFI [second term in Eq. (16)]. Once the parameter ΛΛ\Lambdaroman_Λ quantifies the decrease rate in spatial coherence, and since PM correlation values imply that high momentum values are correlated with high spatial values, as a result, we have an extended wavepacket for the probe state. This consequently causes an increase in the decoherence rate τdec1superscriptsubscript𝜏dec1\tau_{\text{dec}}^{-1}italic_τ start_POSTSUBSCRIPT dec end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT [see Eq. (6)], reducing the purity of the state. In our model, we considered the long-wavelength limit [59], where the wavelengths of scattering air molecules are much longer than the wave packet extension of the probe. Consequently, a large number of scattering events are needed to resolve the fullerene position by encoding substantial which-path information in the environment, ensuring weak coupling and resulting in an effective Markovian bath. Once we consider that the momentum of environment particles obeys a Maxwell–Boltzmann distribution, the temperature regime investigated here is still more than one million times greater than the ultra-low temperatures obtained, for example, in a Bose-Einstein condensate on the order of 170 nK [67]. Similarly, for example, the amplitude-damping channel [68], which models physical processes such as spontaneous emission, occurs at T=0𝑇0T=0italic_T = 0 K, and even so, Markovianity is maintained since the coupling of the system with the infinite vacuum modes is weak.

The findings discussed here can be extended to more general scenarios. For instance, recently it was shown that incorporating additional degrees of freedom in molecules, such as rotations and vibrations, can significantly enhance sensitivity when using molecular sensors as probes [69, 70], creating new opportunities for their application in quantum protocols. Building on this, one could integrate these degrees of freedom into our estimation protocol to fully explore the advantages of molecular probes. Additionally, our findings could also be expanded to include more complex dynamics through a generalized channel that accounts for non-Markovian memory processes. This would offer deeper insights into non-Markovian decoherence features in structured reservoirs, such as decoherence-free states [71]. By pursuing this approach, one could also advance in the exploration at even lower temperatures by relaxing the Born-Markov approximation [18, 72].

Acknowledgements.
J. C. P. Porto acknowledges Fundação de Amparo à Pesquisa do Estado do Piauí (FAPEPI) for financial support. L.S.M. acknowledges the Federal University of Piauí for providing the workspace. P.R.D acknowledges support from the NCN Poland, ChistEra-2023/05/Y/ST2/00005 under the project Modern Device Independent Cryptography (MoDIC). I.G.P. acknowledges Grant No. 306528/2023-1 from CNPq. C.H.S.V acknowledges the São Paulo Research Foundation (FAPESP), Grant. No. 2023/13362-0, for financial support and the Federal University of ABC (UFABC) to provide the workspace.

Appendix A Fisher information

Here, we review the definitions and differences between classical and quantum Fisher information, as well as the definition of purity and its relation to resource theory. Our focus is on describing valid relations specifically for Gaussian states.

A.1 Classical Fisher Information

The Fisher information is a measure of information that an observable random variable carries about an unknown parameter [73]. It expresses the level of uncertainty of a measured physical quantity. Let P(xi|Θ)𝑃conditionalsubscript𝑥𝑖ΘP(x_{i}|\Theta)italic_P ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT | roman_Θ ) be the conditioned probability of measuring data xisubscript𝑥𝑖x_{i}italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT given a specific value of ΘΘ\Thetaroman_Θ. The Classical Fisher Information (CFI) for estimation of the parameter ΘΘ\Thetaroman_Θ is defined by

Θ C(Θ)=iP(xi|Θ)(ln[P(xi|Θ)]Θ)2=i1P(xi|Θ)(P(xi|Θ)Θ)2.superscriptsubscriptΘ CΘsubscript𝑖𝑃conditionalsubscript𝑥𝑖Θsuperscript𝑃conditionalsubscript𝑥𝑖ΘΘ2subscript𝑖1𝑃conditionalsubscript𝑥𝑖Θsuperscript𝑃conditionalsubscript𝑥𝑖ΘΘ2\displaystyle\mathcal{F}_{\Theta}^{\text{\;C}}(\Theta)=\sum_{i}P(x_{i}|\Theta)% \Bigg{(}\frac{\partial\ln[P(x_{i}|\Theta)]}{\partial\Theta}\Bigg{)}^{2}=\sum_{% i}\frac{1}{P(x_{i}|\Theta)}\Bigg{(}\frac{\partial P(x_{i}|\Theta)}{\partial% \Theta}\Bigg{)}^{2}.caligraphic_F start_POSTSUBSCRIPT roman_Θ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT C end_POSTSUPERSCRIPT ( roman_Θ ) = ∑ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_P ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT | roman_Θ ) ( divide start_ARG ∂ roman_ln [ italic_P ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT | roman_Θ ) ] end_ARG start_ARG ∂ roman_Θ end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = ∑ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT divide start_ARG 1 end_ARG start_ARG italic_P ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT | roman_Θ ) end_ARG ( divide start_ARG ∂ italic_P ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT | roman_Θ ) end_ARG start_ARG ∂ roman_Θ end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT . (20)

If P(xi|Θ)𝑃conditionalsubscript𝑥𝑖ΘP(x_{i}|\Theta)italic_P ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT | roman_Θ ) exhibits a peak as a response to variations in ΘΘ\Thetaroman_Θ, then, the data xisubscript𝑥𝑖x_{i}italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT provides the information needed to estimate the parameter ΘΘ\Thetaroman_Θ. On the other hand, when P(xi|Θ)𝑃conditionalsubscript𝑥𝑖ΘP(x_{i}|\Theta)italic_P ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT | roman_Θ ) is flat, numerous samples of xisubscript𝑥𝑖x_{i}italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT would be required to estimate the value of ΘΘ\Thetaroman_Θ, which would be determined only by using the complete sample population. In this context, the Cramér-Rao inequality [74] establishes the lower bound on standard deviation ΔΘ=Θ2Θ2ΔΘdelimited-⟨⟩superscriptΘ2superscriptdelimited-⟨⟩Θ2\Delta\Theta=\sqrt{\langle\Theta^{2}\rangle-\langle\Theta\rangle^{2}}roman_Δ roman_Θ = square-root start_ARG ⟨ roman_Θ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ⟩ - ⟨ roman_Θ ⟩ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG of the estimated parameter

ΔΘ1NΘ C(Θ),ΔΘ1𝑁superscriptsubscriptΘ CΘ\Delta\Theta\geq\frac{1}{\sqrt{N\mathcal{F}_{\Theta}^{\text{\;C}}(\Theta)}},roman_Δ roman_Θ ≥ divide start_ARG 1 end_ARG start_ARG square-root start_ARG italic_N caligraphic_F start_POSTSUBSCRIPT roman_Θ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT C end_POSTSUPERSCRIPT ( roman_Θ ) end_ARG end_ARG , (21)

where N𝑁Nitalic_N is the number of times that the experiment is repeated. Therefore, the Fisher information determines the reachable accuracy of the estimated quantity and represents the figure of merit in parameter estimation problems. Import to mention, that this inequality is valid only for unbiased estimators, i.e., estimators for which Θest=Θrealdelimited-⟨⟩subscriptΘ𝑒𝑠𝑡subscriptΘ𝑟𝑒𝑎𝑙\langle\Theta_{est}\rangle=\Theta_{real}⟨ roman_Θ start_POSTSUBSCRIPT italic_e italic_s italic_t end_POSTSUBSCRIPT ⟩ = roman_Θ start_POSTSUBSCRIPT italic_r italic_e italic_a italic_l end_POSTSUBSCRIPT.

In quantum theory, we use a set of positive operator-values measure (POVM) {A^(κ)}^𝐴𝜅\{\hat{A}(\kappa)\}{ over^ start_ARG italic_A end_ARG ( italic_κ ) }, parameterized by κ𝜅\kappaitalic_κ, to describe the measurement procedure. These operators are positive and satisfy κA^(κ)=I^subscript𝜅^𝐴𝜅^𝐼\sum_{\kappa}\hat{A}(\kappa)=\hat{I}∑ start_POSTSUBSCRIPT italic_κ end_POSTSUBSCRIPT over^ start_ARG italic_A end_ARG ( italic_κ ) = over^ start_ARG italic_I end_ARG to ensure normalization. The probability can be expressed as P(κ|Θ)=Tr[ρ^{A^(κ)}]𝑃conditional𝜅ΘTrdelimited-[]^𝜌^𝐴𝜅P(\kappa|\Theta)=\text{Tr}[\hat{\rho}\{\hat{A}(\kappa)\}]italic_P ( italic_κ | roman_Θ ) = Tr [ over^ start_ARG italic_ρ end_ARG { over^ start_ARG italic_A end_ARG ( italic_κ ) } ]. Using P(κ|Θ)𝑃conditional𝜅ΘP(\kappa|\Theta)italic_P ( italic_κ | roman_Θ ) as the probability distribution, the CFI is defined as

ΘC(Θ;A^)=κ1P(κ|Θ)(P(κ|Θ)Θ)2,superscriptsubscriptΘCΘ^𝐴subscript𝜅1𝑃conditional𝜅Θsuperscript𝑃conditional𝜅ΘΘ2\displaystyle\mathcal{F}_{\Theta}^{\;\text{C}}(\Theta;\hat{A})=\sum_{\kappa}% \frac{1}{P(\kappa|\Theta)}\Bigg{(}\frac{\partial P(\kappa|\Theta)}{\partial% \Theta}\Bigg{)}^{2},caligraphic_F start_POSTSUBSCRIPT roman_Θ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT C end_POSTSUPERSCRIPT ( roman_Θ ; over^ start_ARG italic_A end_ARG ) = ∑ start_POSTSUBSCRIPT italic_κ end_POSTSUBSCRIPT divide start_ARG 1 end_ARG start_ARG italic_P ( italic_κ | roman_Θ ) end_ARG ( divide start_ARG ∂ italic_P ( italic_κ | roman_Θ ) end_ARG start_ARG ∂ roman_Θ end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , (22)

where ρ^^𝜌\hat{\rho}over^ start_ARG italic_ρ end_ARG is the density operator describing the system. Naturally, optimal POVMs are those characterized by a statistical distribution of measurement results that is maximally sensitive to changes in the parameter ΘΘ\Thetaroman_Θ [75].

A.2 Quantum Fisher Information

The Quantum Fisher information (QFI), here denoted by Θ Qsubscriptsuperscript QΘ\mathcal{F}^{\text{\;Q}}_{\Theta}caligraphic_F start_POSTSUPERSCRIPT Q end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_Θ end_POSTSUBSCRIPT, is defined by maximizing Eq. (22) over all possible POVMs {A^(κ)}^𝐴𝜅\{\hat{A}(\kappa)\}{ over^ start_ARG italic_A end_ARG ( italic_κ ) } as follows [76, 77, 78, 79]

Θ Q(Θ)=maxA^Θ C(Θ;A^).superscriptsubscriptΘ QΘsubscript^𝐴superscriptsubscriptΘ CΘ^𝐴\mathcal{F}_{\Theta}^{\text{\;Q}}(\Theta)=\max_{\hat{A}}\mathcal{F}_{\Theta}^{% \text{\;C}}(\Theta;\hat{A}).caligraphic_F start_POSTSUBSCRIPT roman_Θ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT Q end_POSTSUPERSCRIPT ( roman_Θ ) = roman_max start_POSTSUBSCRIPT over^ start_ARG italic_A end_ARG end_POSTSUBSCRIPT caligraphic_F start_POSTSUBSCRIPT roman_Θ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT C end_POSTSUPERSCRIPT ( roman_Θ ; over^ start_ARG italic_A end_ARG ) . (23)

The quantum estimation is then related to the optimal possible measurement precision. Note that this maximization procedure turns the QFI an upper bound for the classical one, i.e., ΘQ(Θ)ΘC(Θ)superscriptsubscriptΘQΘsuperscriptsubscriptΘCΘ\mathcal{F}_{\Theta}^{\;\text{Q}}(\Theta)\geq\mathcal{F}_{\Theta}^{\;\text{C}}% (\Theta)caligraphic_F start_POSTSUBSCRIPT roman_Θ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT Q end_POSTSUPERSCRIPT ( roman_Θ ) ≥ caligraphic_F start_POSTSUBSCRIPT roman_Θ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT C end_POSTSUPERSCRIPT ( roman_Θ ) [78]. Quantum metrology is not the only application of quantum Fisher information, alternatives include quantum cloning [80], entanglement detection  [81, 82], and quantum phase transition [83]. Feng and Wei  [84] have highlighted the relationship between quantum coherence and Quantum Fisher Information (QFI), demonstrating that QFI is useful for quantifying quantum coherence. This relationship is established because QFI satisfies monotonicity under typical incoherent operations and convexity when quantum states are mixed. In Ref. [85], it was demonstrated how QFI can be used to identify quantum correlations such as steering, showing that such correlations can be useful for quantum-enhanced measurement protocols.

For a single-mode Gaussian state, with covariance matrice 𝝈𝝈\boldsymbol{\sigma}bold_italic_σ and first moments 𝒅𝒅\boldsymbol{d}bold_italic_d, the QFI is given by [62, 86]

ΘQ(Θ)=Tr[(𝝈𝟏Θ𝝈)2]2(1+μ2)+2(Θμ)21μ4+2(Θ𝒅)T(𝝈𝟏)(Θ𝒅),superscriptsubscriptΘQΘTrdelimited-[]superscriptsuperscript𝝈1subscriptΘ𝝈221superscript𝜇22superscriptsubscriptΘ𝜇21superscript𝜇42superscriptsubscriptΘ𝒅Tsuperscript𝝈1subscriptΘ𝒅\mathcal{F}_{\Theta}^{\;\text{Q}}(\Theta)=\frac{\text{Tr}[(\boldsymbol{\sigma^% {-1}}\partial_{\Theta}\boldsymbol{\sigma})^{2}]}{2(1+\mu^{2})}+2\frac{(% \partial_{\Theta}\mu)^{2}}{1-\mu^{4}}+2(\partial_{\Theta}\boldsymbol{d})^{% \text{T}}(\boldsymbol{\sigma^{-1}})(\partial_{\Theta}\boldsymbol{d}),caligraphic_F start_POSTSUBSCRIPT roman_Θ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT Q end_POSTSUPERSCRIPT ( roman_Θ ) = divide start_ARG Tr [ ( bold_italic_σ start_POSTSUPERSCRIPT bold_- bold_1 end_POSTSUPERSCRIPT ∂ start_POSTSUBSCRIPT roman_Θ end_POSTSUBSCRIPT bold_italic_σ ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] end_ARG start_ARG 2 ( 1 + italic_μ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG + 2 divide start_ARG ( ∂ start_POSTSUBSCRIPT roman_Θ end_POSTSUBSCRIPT italic_μ ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 1 - italic_μ start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG + 2 ( ∂ start_POSTSUBSCRIPT roman_Θ end_POSTSUBSCRIPT bold_italic_d ) start_POSTSUPERSCRIPT T end_POSTSUPERSCRIPT ( bold_italic_σ start_POSTSUPERSCRIPT bold_- bold_1 end_POSTSUPERSCRIPT ) ( ∂ start_POSTSUBSCRIPT roman_Θ end_POSTSUBSCRIPT bold_italic_d ) , (24)

where μ=1/det(𝝈)𝜇1det𝝈\mu=1/\sqrt{\text{det}(\boldsymbol{\sigma})}italic_μ = 1 / square-root start_ARG det ( bold_italic_σ ) end_ARG is the purity of the quantum state and ΘsubscriptΘ\partial_{\Theta}∂ start_POSTSUBSCRIPT roman_Θ end_POSTSUBSCRIPT denotes differentiation with respect to the parameter ΘΘ\Thetaroman_Θ. The first term is associated with the dynamical dependence of the covariance matrix with the ΘΘ\Thetaroman_Θ parameter. The second one is the dynamic of purity under ΘΘ\Thetaroman_Θ variation, and the third is the contribution of the moments dynamics of the Gaussian state for the estimated parameter. For example, this result has been utilized in metrological protocols that take advantage of the superradiant phase transition in the Rabi model [87]. For our system under investigation, we will apply the result (24) throughout the work to quantify the QFI.

Appendix B How generate PM-correlated Gaussian states?

We review the procedure that can generate correlated Gaussian states for pedagogical reasons and to clarify the work from an experimental point of view. We will explicitly show the parameters determining the initial position-momentum correlation γ𝛾\gammaitalic_γ in Eq. (1). The traditional method known as atom lens [88, 89] can be applied to produce these initial position-momentum correlations. This method uses the atom-light interactions to create spatially dependent AC-Stark shift of the electronic ground state of the atoms which is induced in the vicinity of an intensity maximum of a sub-resonant standing wave laser field [30]. The standing wave is generated using a cavity perpendicular to the atom direction (see Fig. 1). If the atoms are modeled as two-level systems, the atom-laser interaction’s Rabi frequency is given by

Ω(x,z)=Ω0cos(2πx/λ)eπz2/(VCMtint)2Ω𝑥𝑧subscriptΩ02𝜋𝑥𝜆superscript𝑒𝜋superscript𝑧2superscriptsubscript𝑉CMsubscript𝑡int2\Omega(x,z)=\Omega_{0}\cos(2\pi x/\lambda)e^{-\pi z^{2}/(V_{\text{CM}}t_{\text% {int}})^{2}}roman_Ω ( italic_x , italic_z ) = roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT roman_cos ( 2 italic_π italic_x / italic_λ ) italic_e start_POSTSUPERSCRIPT - italic_π italic_z start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / ( italic_V start_POSTSUBSCRIPT CM end_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT int end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT (25)

to describe a standing wave with period λ𝜆\lambdaitalic_λ (light wavelength) in the x𝑥xitalic_x-direction and a Gaussian beam profile in the z𝑧zitalic_z-direction. Here, VCMsubscript𝑉CMV_{\text{CM}}italic_V start_POSTSUBSCRIPT CM end_POSTSUBSCRIPT is the atomic center-of-mass velocity, and tintsubscript𝑡intt_{\text{int}}italic_t start_POSTSUBSCRIPT int end_POSTSUBSCRIPT is the effective interaction time, assumed smaller than the radiative lifetime of the excited state of the atoms so that spontaneous emission can be neglected. This approach also applies to the focusing of molecules, since we can consider their vibrational states practically degenerate as long as the energy difference between them is much smaller than the excitation energy of the light beam inside the optical cavity. Due to the Stark shift, the optical potential for incoming atoms in the ground state is given by

U(x,z)=12Ω(x,z)2+δ2,𝑈𝑥𝑧12Ωsuperscript𝑥𝑧2superscript𝛿2U(x,z)=-\frac{1}{2}\sqrt{\Omega(x,z)^{2}+\delta^{2}},italic_U ( italic_x , italic_z ) = - divide start_ARG 1 end_ARG start_ARG 2 end_ARG square-root start_ARG roman_Ω ( italic_x , italic_z ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_δ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG , (26)

where δ=ωω0𝛿𝜔subscript𝜔0\delta=\omega-\omega_{0}italic_δ = italic_ω - italic_ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, with ω=2πc/λ𝜔2𝜋𝑐𝜆\omega=2\pi c/\lambdaitalic_ω = 2 italic_π italic_c / italic_λ is the laser frequency and ω0subscript𝜔0\omega_{0}italic_ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is the resonance transition [30]. In the vicinity of a maximum intensity region, the atoms roughly feel the harmonic potential U(x)14Ω0k2x2(1+δ2/Ω02)𝑈𝑥14subscriptΩ0superscript𝑘2superscript𝑥21superscript𝛿2superscriptsubscriptΩ02U(x)\approx\frac{1}{4}\Omega_{0}k^{2}x^{2}(1+\delta^{2}/\Omega_{0}^{2})italic_U ( italic_x ) ≈ divide start_ARG 1 end_ARG start_ARG 4 end_ARG roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( 1 + italic_δ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ), with photon momentum k=2π/λ𝑘2𝜋𝜆k=2\pi/\lambdaitalic_k = 2 italic_π / italic_λ. For a short interaction time, the harmonic potential changes the initial state to

ψ0(x)eimVCMx22R(z)ψ0(x),subscript𝜓0𝑥superscript𝑒𝑖𝑚subscript𝑉CMsuperscript𝑥22Planck-constant-over-2-pi𝑅𝑧subscript𝜓0𝑥\psi_{0}(x)\rightarrow e^{\frac{imV_{\text{CM}}x^{2}}{2\hbar R(z)}}\psi_{0}(x),italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) → italic_e start_POSTSUPERSCRIPT divide start_ARG italic_i italic_m italic_V start_POSTSUBSCRIPT CM end_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 roman_ℏ italic_R ( italic_z ) end_ARG end_POSTSUPERSCRIPT italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) , (27)

where R(z)𝑅𝑧R(z)italic_R ( italic_z ) is the curvature radius of the wavefronts associated with the beam propagation. Then, it acts as a thin lens with focal length [30]

f=λ2πΩ0tintλdB(1+δ2Ω02)1/2,𝑓superscript𝜆2𝜋subscriptΩ0subscript𝑡𝑖𝑛𝑡subscript𝜆dBsuperscript1superscript𝛿2superscriptsubscriptΩ0212f=\frac{\lambda^{2}}{\pi\Omega_{0}t_{int}\lambda_{\text{dB}}}\Bigg{(}1+\frac{% \delta^{2}}{\Omega_{0}^{2}}\Bigg{)}^{1/2},italic_f = divide start_ARG italic_λ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_π roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT italic_i italic_n italic_t end_POSTSUBSCRIPT italic_λ start_POSTSUBSCRIPT dB end_POSTSUBSCRIPT end_ARG ( 1 + divide start_ARG italic_δ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT , (28)

where λdB=h/(mVCM)subscript𝜆dB𝑚subscript𝑉CM\lambda_{\text{dB}}=h/(mV_{\text{CM}})italic_λ start_POSTSUBSCRIPT dB end_POSTSUBSCRIPT = italic_h / ( italic_m italic_V start_POSTSUBSCRIPT CM end_POSTSUBSCRIPT ) is the de Broglie wavelength of the atoms. We see therefore from this analysis which experimental parameters determine the initial position correlation characterized by the γ𝛾\gammaitalic_γ parameter encoded in the quadratic phase of the initial state (1). From an experimental point of view, the focal length f𝑓fitalic_f can be varied, for example, by modulating the laser power. Also, note that positive (negative) values of γ𝛾\gammaitalic_γ are associated with a diverging (converging) beam with curvature radius R>0𝑅0R>0italic_R > 0 (R<0𝑅0R<0italic_R < 0).

Appendix C Parameters of density matrix for a correlated Gaussian state and purity

Here, we describe the parameters obtained through the action of the time evolution propagator over the density matrix shown in Eq. (7)

𝒜t=A1+A2iA3,t=A1A2iA3,𝒞t=2iA3,formulae-sequencesubscript𝒜𝑡subscript𝐴1subscript𝐴2𝑖subscript𝐴3formulae-sequencesubscript𝑡subscript𝐴1subscript𝐴2𝑖subscript𝐴3subscript𝒞𝑡2𝑖subscript𝐴3\displaystyle\mathcal{A}_{t}=A_{1}+A_{2}-iA_{3},\quad\mathcal{B}_{t}=A_{1}-A_{% 2}-iA_{3},\quad\mathcal{C}_{t}=2iA_{3},caligraphic_A start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - italic_i italic_A start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT , caligraphic_B start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - italic_i italic_A start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT , caligraphic_C start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = 2 italic_i italic_A start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT , (29)
A1=m282t2σ02B2,B2=14σ04+1202σ02+(m2t+γ2σ02)2+Λt3σ02,formulae-sequencesubscript𝐴1superscript𝑚28superscriptPlanck-constant-over-2-pi2superscript𝑡2superscriptsubscript𝜎02superscript𝐵2superscript𝐵214superscriptsubscript𝜎0412superscriptsubscript02superscriptsubscript𝜎02superscript𝑚2Planck-constant-over-2-pi𝑡𝛾2superscriptsubscript𝜎022Λ𝑡3superscriptsubscript𝜎02\displaystyle A_{1}=\frac{m^{2}}{8\hbar^{2}t^{2}\sigma_{0}^{2}B^{2}},\quad B^{% 2}=\frac{1}{4\sigma_{0}^{4}}+\frac{1}{2\ell_{0}^{2}\sigma_{0}^{2}}+\left(\frac% {m}{2\hbar t}+\frac{\gamma}{2\sigma_{0}^{2}}\right)^{2}+\frac{\Lambda t}{3% \sigma_{0}^{2}},italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = divide start_ARG italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 8 roman_ℏ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_B start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG , italic_B start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = divide start_ARG 1 end_ARG start_ARG 4 italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG + divide start_ARG 1 end_ARG start_ARG 2 roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + ( divide start_ARG italic_m end_ARG start_ARG 2 roman_ℏ italic_t end_ARG + divide start_ARG italic_γ end_ARG start_ARG 2 italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + divide start_ARG roman_Λ italic_t end_ARG start_ARG 3 italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG , (30)
A2=m242t2B2(1202+Λt)+Λt12σ02B2(Λt+12σ02+202)+mΛγ4σ02B2+Λtγ212σ04B2,subscript𝐴2superscript𝑚24superscriptPlanck-constant-over-2-pi2superscript𝑡2superscript𝐵212superscriptsubscript02Λ𝑡Λ𝑡12superscriptsubscript𝜎02superscript𝐵2Λ𝑡12superscriptsubscript𝜎022superscriptsubscript02𝑚Λ𝛾4Planck-constant-over-2-pisuperscriptsubscript𝜎02superscript𝐵2Λ𝑡superscript𝛾212superscriptsubscript𝜎04superscript𝐵2\displaystyle A_{2}=\frac{m^{2}}{4\hbar^{2}t^{2}B^{2}}\left(\frac{1}{2\ell_{0}% ^{2}}+\Lambda t\right)+\frac{\Lambda t}{12\sigma_{0}^{2}B^{2}}\left(\Lambda t+% \frac{1}{2\sigma_{0}^{2}}+\frac{2}{\ell_{0}^{2}}\right)+\frac{m\Lambda\gamma}{% 4\hbar\sigma_{0}^{2}B^{2}}+\frac{\Lambda t\gamma^{2}}{12\sigma_{0}^{4}B^{2}},italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = divide start_ARG italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 4 roman_ℏ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_B start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ( divide start_ARG 1 end_ARG start_ARG 2 roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + roman_Λ italic_t ) + divide start_ARG roman_Λ italic_t end_ARG start_ARG 12 italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_B start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ( roman_Λ italic_t + divide start_ARG 1 end_ARG start_ARG 2 italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + divide start_ARG 2 end_ARG start_ARG roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) + divide start_ARG italic_m roman_Λ italic_γ end_ARG start_ARG 4 roman_ℏ italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_B start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + divide start_ARG roman_Λ italic_t italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 12 italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_B start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG , (31)

and

A3=m4tσ02B2(Λt+12σ02+102)+mγ8tσ02B2(mt+γσ02),𝒩t=2A1π.formulae-sequencesubscript𝐴3𝑚4Planck-constant-over-2-pi𝑡superscriptsubscript𝜎02superscript𝐵2Λ𝑡12superscriptsubscript𝜎021superscriptsubscript02𝑚𝛾8Planck-constant-over-2-pi𝑡superscriptsubscript𝜎02superscript𝐵2𝑚Planck-constant-over-2-pi𝑡𝛾superscriptsubscript𝜎02subscript𝒩𝑡2subscript𝐴1𝜋\displaystyle A_{3}=\frac{m}{4\hbar t\sigma_{0}^{2}B^{2}}\left(\Lambda t+\frac% {1}{2\sigma_{0}^{2}}+\frac{1}{\ell_{0}^{2}}\right)+\frac{m\gamma}{8\hbar t% \sigma_{0}^{2}B^{2}}\left(\frac{m}{\hbar t}+\frac{\gamma}{\sigma_{0}^{2}}% \right),\;\;\;\;\;\;\;\mathcal{N}_{t}=\sqrt{\frac{2A_{1}}{\pi}}.italic_A start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT = divide start_ARG italic_m end_ARG start_ARG 4 roman_ℏ italic_t italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_B start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ( roman_Λ italic_t + divide start_ARG 1 end_ARG start_ARG 2 italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + divide start_ARG 1 end_ARG start_ARG roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) + divide start_ARG italic_m italic_γ end_ARG start_ARG 8 roman_ℏ italic_t italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_B start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ( divide start_ARG italic_m end_ARG start_ARG roman_ℏ italic_t end_ARG + divide start_ARG italic_γ end_ARG start_ARG italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) , caligraphic_N start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = square-root start_ARG divide start_ARG 2 italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_π end_ARG end_ARG . (32)

Besides this, employing the density matrix, Eq. (7), and the parameters (C1-C4), the purity Eq. (8), after some algebraic simplifications, can be explicitly written as follows

μ(γ,Λ,t)=[1+2σ0202+4σ02Λt+4γΛmt2+4Λ(γ2+1+2σ0202)3τ0mt3+4Λ223m2t4]1/2.𝜇𝛾Λ𝑡superscriptdelimited-[]12superscriptsubscript𝜎02superscriptsubscript024superscriptsubscript𝜎02Λ𝑡4𝛾ΛPlanck-constant-over-2-pi𝑚superscript𝑡24Planck-constant-over-2-piΛsuperscript𝛾212superscriptsubscript𝜎02superscriptsubscript023subscript𝜏0𝑚superscript𝑡34superscriptΛ2superscriptPlanck-constant-over-2-pi23superscript𝑚2superscript𝑡412\displaystyle\mu(\gamma,\Lambda,t)=\Bigg{[}1+\frac{2\sigma_{0}^{2}}{\ell_{0}^{% 2}}+4\sigma_{0}^{2}\Lambda t+\frac{4\gamma\Lambda\hbar}{m}t^{2}+\frac{4\hbar% \Lambda(\gamma^{2}+1+\frac{2\sigma_{0}^{2}}{\ell_{0}^{2}})}{3\tau_{0}m}t^{3}+% \frac{4\Lambda^{2}\hbar^{2}}{3m^{2}}t^{4}\Bigg{]}^{-1/2}.italic_μ ( italic_γ , roman_Λ , italic_t ) = [ 1 + divide start_ARG 2 italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + 4 italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_Λ italic_t + divide start_ARG 4 italic_γ roman_Λ roman_ℏ end_ARG start_ARG italic_m end_ARG italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + divide start_ARG 4 roman_ℏ roman_Λ ( italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 1 + divide start_ARG 2 italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) end_ARG start_ARG 3 italic_τ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_m end_ARG italic_t start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT + divide start_ARG 4 roman_Λ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_ℏ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 3 italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG italic_t start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT ] start_POSTSUPERSCRIPT - 1 / 2 end_POSTSUPERSCRIPT . (33)

Appendix D QFI parameters for correlation estimation

In this appendix, we treat the parameters associated with the function on Eq.(14), representing the QFI for correlation estimation.

Φγ(γ,t,Λ)=17202τ04k=02𝒞kΛk,subscriptΦ𝛾𝛾𝑡Λ172superscriptsubscript02superscriptsubscript𝜏04superscriptsubscript𝑘02subscript𝒞𝑘superscriptΛ𝑘\displaystyle\Phi_{\gamma}(\gamma,t,\Lambda)=\frac{1}{72\ell_{0}^{2}\tau_{0}^{% 4}}\sum_{k=0}^{2}\mathcal{C}_{k}\Lambda^{k},roman_Φ start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT ( italic_γ , italic_t , roman_Λ ) = divide start_ARG 1 end_ARG start_ARG 72 roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_τ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG ∑ start_POSTSUBSCRIPT italic_k = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT caligraphic_C start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT roman_Λ start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT , (34)

where 𝒞k=𝒞k(γ,t,Λ)subscript𝒞𝑘subscript𝒞𝑘𝛾𝑡Λ\mathcal{C}_{k}=\mathcal{C}_{k}(\gamma,t,\Lambda)caligraphic_C start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT = caligraphic_C start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ( italic_γ , italic_t , roman_Λ ), are expressed as follows

𝒞0=902τ04[1+2(σ00)2],𝒞1=12σ0202τ02t3[(2(σ00)2+(γ2+1))+3γτ0t+3(τ0t)2],formulae-sequencesubscript𝒞09superscriptsubscript02superscriptsubscript𝜏04delimited-[]12superscriptsubscript𝜎0subscript02subscript𝒞112superscriptsubscript𝜎02superscriptsubscript02superscriptsubscript𝜏02superscript𝑡3delimited-[]2superscriptsubscript𝜎0subscript02superscript𝛾213𝛾subscript𝜏0𝑡3superscriptsubscript𝜏0𝑡2\displaystyle\mathcal{C}_{0}=9\ell_{0}^{2}\tau_{0}^{4}\left[1+2\left(\frac{% \sigma_{0}}{\ell_{0}}\right)^{2}\right],\;\;\;\;\mathcal{C}_{1}=12\sigma_{0}^{% 2}\ell_{0}^{2}\tau_{0}^{2}t^{3}\left[\left(2\left(\frac{\sigma_{0}}{\ell_{0}}% \right)^{2}+(\gamma^{2}+1)\right)+3\gamma\frac{\tau_{0}}{t}+3\left(\frac{\tau_% {0}}{t}\right)^{2}\right],caligraphic_C start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 9 roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_τ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT [ 1 + 2 ( divide start_ARG italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] , caligraphic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = 12 italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_τ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_t start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT [ ( 2 ( divide start_ARG italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 1 ) ) + 3 italic_γ divide start_ARG italic_τ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_t end_ARG + 3 ( divide start_ARG italic_τ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_t end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] , (35)

and

𝒞2=32σ0402t6[γ2+3γτ0t+218(τ0t)2].subscript𝒞232superscriptsubscript𝜎04superscriptsubscript02superscript𝑡6delimited-[]superscript𝛾23𝛾subscript𝜏0𝑡218superscriptsubscript𝜏0𝑡2\displaystyle\mathcal{C}_{2}=32\sigma_{0}^{4}\ell_{0}^{2}t^{6}\left[\gamma^{2}% +3\gamma\frac{\tau_{0}}{t}+\frac{21}{8}\left(\frac{\tau_{0}}{t}\right)^{2}% \right].caligraphic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = 32 italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_t start_POSTSUPERSCRIPT 6 end_POSTSUPERSCRIPT [ italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 3 italic_γ divide start_ARG italic_τ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_t end_ARG + divide start_ARG 21 end_ARG start_ARG 8 end_ARG ( divide start_ARG italic_τ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_t end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] . (36)

We show in Fig. 6 (a) QFI (γ Qsuperscriptsubscript𝛾 Q\mathcal{F}_{\gamma}^{\text{\;Q}}caligraphic_F start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT Q end_POSTSUPERSCRIPT), (b) purity μ𝜇\muitalic_μ and (c) its absolute relative derivative 1μ|γμ|1𝜇subscript𝛾𝜇\frac{1}{\mu}|\partial_{\gamma}\mu|divide start_ARG 1 end_ARG start_ARG italic_μ end_ARG | ∂ start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT italic_μ | as a function of initial correlation γ𝛾\gammaitalic_γ and the propagation time t𝑡titalic_t, considering Λ=1022Λsuperscript1022\Lambda=10^{22}roman_Λ = 10 start_POSTSUPERSCRIPT 22 end_POSTSUPERSCRIPT m-2s-1 (T=20.5×103𝑇20.5superscript103T=20.5\times 10^{3}italic_T = 20.5 × 10 start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT K). Observe that the QFI behavior is directly influenced by the purity (absolute variation in purity) in the regime of small (larger) propagation times t𝑡titalic_t, allowing us to describe how the transition occurs between a regime determined by purity or its rate of variation in the behavior of the QFI.

Refer to caption
Refer to caption
Refer to caption
Fig. 6: (a) QFI (γ Qsuperscriptsubscript𝛾 Q\mathcal{F}_{\gamma}^{\text{\;Q}}caligraphic_F start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT Q end_POSTSUPERSCRIPT), (b) purity μ𝜇\muitalic_μ and (c) its absolute relative derivative 1μ|γμ|1𝜇subscript𝛾𝜇\frac{1}{\mu}|\partial_{\gamma}\mu|divide start_ARG 1 end_ARG start_ARG italic_μ end_ARG | ∂ start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT italic_μ | as a function of initial correlation γ𝛾\gammaitalic_γ and the propagation time t𝑡titalic_t, considering Λ=1022Λsuperscript1022\Lambda=10^{22}roman_Λ = 10 start_POSTSUPERSCRIPT 22 end_POSTSUPERSCRIPT m-2s-1 (T=20.5×103𝑇20.5superscript103T=20.5\times 10^{3}italic_T = 20.5 × 10 start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT K). We can observe that, for small (larger) propagation times t𝑡titalic_t, the QFI behavior is governed by the purity (absolute variation in purity).

Appendix E QFI parameters for the effective environmental coupling estimation

In this appendix, let us show the constants associated with the estimation of the environment coupling through the QFI [ΛC(γ,Λ,t)superscriptsubscriptΛC𝛾Λ𝑡\mathcal{F}_{\Lambda}^{\text{C}}(\gamma,\Lambda,t)caligraphic_F start_POSTSUBSCRIPT roman_Λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT C end_POSTSUPERSCRIPT ( italic_γ , roman_Λ , italic_t )].

ΦΛ(γ,t,Λ)=11804τ04k=02𝒵¯kΛk,subscriptΦΛ𝛾𝑡Λ118superscriptsubscript04superscriptsubscript𝜏04superscriptsubscript𝑘02subscript¯𝒵𝑘superscriptΛ𝑘\displaystyle\Phi_{\Lambda}(\gamma,t,\Lambda)=\frac{1}{18\ell_{0}^{4}\tau_{0}^% {4}}\sum_{k=0}^{2}\bar{\mathcal{Z}}_{k}\Lambda^{k},roman_Φ start_POSTSUBSCRIPT roman_Λ end_POSTSUBSCRIPT ( italic_γ , italic_t , roman_Λ ) = divide start_ARG 1 end_ARG start_ARG 18 roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_τ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG ∑ start_POSTSUBSCRIPT italic_k = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over¯ start_ARG caligraphic_Z end_ARG start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT roman_Λ start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT , (37)

where 𝒵¯k=𝒵¯k(γ,t,Λ)subscript¯𝒵𝑘subscript¯𝒵𝑘𝛾𝑡Λ\bar{\mathcal{Z}}_{k}=\bar{\mathcal{Z}}_{k}(\gamma,t,\Lambda)over¯ start_ARG caligraphic_Z end_ARG start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT = over¯ start_ARG caligraphic_Z end_ARG start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ( italic_γ , italic_t , roman_Λ ), expressed as follows

𝒵¯0=204σ04t6[Γ2+6γτ0tΓ+15(τ0t)2(35(σ00)2+(γ2+310))+18γ(τ0t)3+9(τ0t)4],subscript¯𝒵02superscriptsubscript04superscriptsubscript𝜎04superscript𝑡6delimited-[]superscriptΓ26𝛾subscript𝜏0𝑡Γ15superscriptsubscript𝜏0𝑡235superscriptsubscript𝜎0subscript02superscript𝛾231018𝛾superscriptsubscript𝜏0𝑡39superscriptsubscript𝜏0𝑡4\displaystyle\bar{\mathcal{Z}}_{0}=2\ell_{0}^{4}\sigma_{0}^{4}t^{6}\left[% \Gamma^{2}+6\gamma\frac{\tau_{0}}{t}\Gamma+15\left(\frac{\tau_{0}}{t}\right)^{% 2}\left(\frac{3}{5}\left(\frac{\sigma_{0}}{\ell_{0}}\right)^{2}+\left(\gamma^{% 2}+\frac{3}{10}\right)\right)+18\gamma\left(\frac{\tau_{0}}{t}\right)^{3}+9% \left(\frac{\tau_{0}}{t}\right)^{4}\right],over¯ start_ARG caligraphic_Z end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 2 roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_t start_POSTSUPERSCRIPT 6 end_POSTSUPERSCRIPT [ roman_Γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 6 italic_γ divide start_ARG italic_τ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_t end_ARG roman_Γ + 15 ( divide start_ARG italic_τ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_t end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG 3 end_ARG start_ARG 5 end_ARG ( divide start_ARG italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + divide start_ARG 3 end_ARG start_ARG 10 end_ARG ) ) + 18 italic_γ ( divide start_ARG italic_τ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_t end_ARG ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT + 9 ( divide start_ARG italic_τ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_t end_ARG ) start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT ] , (38)
𝒵¯1=4σ0604t7[2(σ00)2+(γ2+1)+3γτ0t+3(τ0t)2],𝒵¯2=4σ0804t8,andΓ=(2(σ00)2+(γ2+1)).formulae-sequencesubscript¯𝒵14superscriptsubscript𝜎06superscriptsubscript04superscript𝑡7delimited-[]2superscriptsubscript𝜎0subscript02superscript𝛾213𝛾subscript𝜏0𝑡3superscriptsubscript𝜏0𝑡2formulae-sequencesubscript¯𝒵24superscriptsubscript𝜎08superscriptsubscript04superscript𝑡8andΓ2superscriptsubscript𝜎0subscript02superscript𝛾21\displaystyle\bar{\mathcal{Z}}_{1}=4\sigma_{0}^{6}\ell_{0}^{4}t^{7}\left[2% \left(\frac{\sigma_{0}}{\ell_{0}}\right)^{2}+\left(\gamma^{2}+1\right)+3\gamma% \frac{\tau_{0}}{t}+3\left(\frac{\tau_{0}}{t}\right)^{2}\right],\;\;\;\bar{% \mathcal{Z}}_{2}=4\sigma_{0}^{8}\ell_{0}^{4}t^{8},\;\;\text{and}\;\;\;\;\Gamma% =\left(2\left(\frac{\sigma_{0}}{\ell_{0}}\right)^{2}+(\gamma^{2}+1)\right).over¯ start_ARG caligraphic_Z end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = 4 italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 6 end_POSTSUPERSCRIPT roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_t start_POSTSUPERSCRIPT 7 end_POSTSUPERSCRIPT [ 2 ( divide start_ARG italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 1 ) + 3 italic_γ divide start_ARG italic_τ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_t end_ARG + 3 ( divide start_ARG italic_τ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_t end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] , over¯ start_ARG caligraphic_Z end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = 4 italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_t start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT , and roman_Γ = ( 2 ( divide start_ARG italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 1 ) ) . (39)

In the Fig. 7, we analyze the dependence of the QFI (Λ2ΛQsuperscriptΛ2superscriptsubscriptΛQ\Lambda^{2}\mathcal{F}_{\Lambda}^{\;\text{Q}}roman_Λ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT caligraphic_F start_POSTSUBSCRIPT roman_Λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT Q end_POSTSUPERSCRIPT) and purity μ𝜇\muitalic_μ as a function of environmental effect ΛΛ\Lambdaroman_Λ. The inset in Fig. 7(b) shows how the absolute relative variation 1μ|Λμ|1𝜇subscriptΛ𝜇\frac{1}{\mu}|\partial_{\Lambda}\mu|divide start_ARG 1 end_ARG start_ARG italic_μ end_ARG | ∂ start_POSTSUBSCRIPT roman_Λ end_POSTSUBSCRIPT italic_μ | directly affects the way the QFI behaves. Furthermore, note that initial correlated states with γ0𝛾0\gamma\neq 0italic_γ ≠ 0 have more QFI when compared with the standard uncorrelated Gaussian state. Also, although for the standard Gaussian state (γ=0𝛾0\gamma=0italic_γ = 0), the purity remains almost unchanged, preserving the value of this quantity does not translate into any gain in QFI. Therefore, the PM-correlations can be employed to enhance the thermal sensitivity of the probe in the regime of small environmental effects (low-temperature regime). Moreover, the range of thermal precision improvement is broad, albeit decreases with temperature growth.

Refer to caption
Fig. 7: (a) QFI (Λ2Λ QsuperscriptΛ2superscriptsubscriptΛ Q\Lambda^{2}\mathcal{F}_{\Lambda}^{\text{\;Q}}roman_Λ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT caligraphic_F start_POSTSUBSCRIPT roman_Λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT Q end_POSTSUPERSCRIPT) and (b) purity μ𝜇\muitalic_μ as a function of environmental effect ΛΛ\Lambdaroman_Λ for t=50.0𝑡50.0t=50.0italic_t = 50.0 μ𝜇\muitalic_μs, and considering the following values for the initial correlation: γ=10𝛾10\gamma=-10italic_γ = - 10, γ=0𝛾0\gamma=0italic_γ = 0 and γ=5𝛾5\gamma=5italic_γ = 5. The inset in (b) shows the absolute relative variation 1μ|Λμ|1𝜇subscriptΛ𝜇\frac{1}{\mu}|\partial_{\Lambda}\mu|divide start_ARG 1 end_ARG start_ARG italic_μ end_ARG | ∂ start_POSTSUBSCRIPT roman_Λ end_POSTSUBSCRIPT italic_μ |, where we can see how it directly affects the way the QFI behaves. Also, note that initial correlated states with γ0𝛾0\gamma\neq 0italic_γ ≠ 0 show a more resilient QFI when compared with the standard uncorrelated Gaussian state, however, this does not necessarily result into a gain in the QFI.

Appendix F Temporal Gain of Information Gain (TGI)

Here we introduce some approximations to obtain a simplified and closed expression for the TGI as a function of PM-correlations γ𝛾\gammaitalic_γ. Considering the propagation times in the order of some microsecond and since (σ0/0)20.024superscriptsubscript𝜎0subscript020.024(\sigma_{0}/\ell_{0})^{2}\approx 0.024( italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT / roman_ℓ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ≈ 0.024 the cubic term in Eq.(33) becomes dominant, and therefore the purity can be approximately expressed as

μapprox(γ,Λ,t)[1+4Λ(γ2+1)3τ0mt3]12.subscript𝜇approx𝛾Λ𝑡superscriptdelimited-[]14Planck-constant-over-2-piΛsuperscript𝛾213subscript𝜏0𝑚superscript𝑡312\displaystyle\mu_{\text{approx}}(\gamma,\Lambda,t)\approx\left[1+\frac{4\hbar% \Lambda(\gamma^{2}+1)}{3\tau_{0}m}t^{3}\right]^{-\frac{1}{2}}.italic_μ start_POSTSUBSCRIPT approx end_POSTSUBSCRIPT ( italic_γ , roman_Λ , italic_t ) ≈ [ 1 + divide start_ARG 4 roman_ℏ roman_Λ ( italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 1 ) end_ARG start_ARG 3 italic_τ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_m end_ARG italic_t start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT ] start_POSTSUPERSCRIPT - divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT . (40)
Table 1: Values with the time improvement in the estimation environmental effect for Λ=1015Λsuperscript1015\Lambda=10^{15}roman_Λ = 10 start_POSTSUPERSCRIPT 15 end_POSTSUPERSCRIPT m-2s-1.
γ𝛾\gammaitalic_γ τmaxsubscript𝜏max\tau_{\text{max}}italic_τ start_POSTSUBSCRIPT max end_POSTSUBSCRIPT (μ𝜇\muitalic_μs) μ(τmax)𝜇subscript𝜏max\mu(\tau_{\text{max}})italic_μ ( italic_τ start_POSTSUBSCRIPT max end_POSTSUBSCRIPT ) 1μ|dμdt|(τmax)1𝜇𝑑𝜇𝑑𝑡subscript𝜏max\frac{1}{\mu}|\frac{d\mu}{dt}|(\tau_{\text{max}})divide start_ARG 1 end_ARG start_ARG italic_μ end_ARG | divide start_ARG italic_d italic_μ end_ARG start_ARG italic_d italic_t end_ARG | ( italic_τ start_POSTSUBSCRIPT max end_POSTSUBSCRIPT ) (s1)superscript𝑠1(s^{-1})( italic_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ) Λ2ΛQ(τmax)superscriptΛ2superscriptsubscriptΛQsubscript𝜏max\Lambda^{2}\mathcal{F}_{\Lambda}^{\text{Q}}(\tau_{\text{max}})roman_Λ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT caligraphic_F start_POSTSUBSCRIPT roman_Λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT Q end_POSTSUPERSCRIPT ( italic_τ start_POSTSUBSCRIPT max end_POSTSUBSCRIPT ) TGI (dB)
-50.0 17.1 0.563 58 488 0.246 11.28
-25.0 27.2 0.563 36 861 0.247 9.24
-1.0 183.2 0.563 5 472 0.247 0.97
0.0 228.4 0.563 4 377 0.247 0
35.0 21.7 0.563 46 117 0.247 10.22
70.0 13.7 0.563 73 191 0.248 12.23
150.0 8.2 0.563 121 646 0.247 14.45

Then, the time τmaxsubscript𝜏max\tau_{\text{max}}italic_τ start_POSTSUBSCRIPT max end_POSTSUBSCRIPT that maximizes the relative purity variation, i.e., the maximization point of 1μapproxtμapprox1subscript𝜇approxsubscript𝑡subscript𝜇approx\frac{1}{\mu_{\text{approx}}}\partial_{t}\mu_{\text{approx}}divide start_ARG 1 end_ARG start_ARG italic_μ start_POSTSUBSCRIPT approx end_POSTSUBSCRIPT end_ARG ∂ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_μ start_POSTSUBSCRIPT approx end_POSTSUBSCRIPT is given by

τmax(Λ,γ)=[3τ022(1+γ2)Λσ02]1/3.subscript𝜏maxΛ𝛾superscriptdelimited-[]3superscriptsubscript𝜏0221superscript𝛾2Λsuperscriptsubscript𝜎0213\tau_{\text{max}}(\Lambda,\gamma)=\Bigg{[}\frac{3\tau_{0}^{2}}{2(1+\gamma^{2})% \Lambda\sigma_{0}^{2}}\Bigg{]}^{1/3}.italic_τ start_POSTSUBSCRIPT max end_POSTSUBSCRIPT ( roman_Λ , italic_γ ) = [ divide start_ARG 3 italic_τ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 ( 1 + italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) roman_Λ italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ] start_POSTSUPERSCRIPT 1 / 3 end_POSTSUPERSCRIPT . (41)

Note that this time can be reduced by increasing the initial correlation γ𝛾\gammaitalic_γ. In this way, using τmax(Λ,γ)subscript𝜏maxΛ𝛾\tau_{\text{max}}(\Lambda,\gamma)italic_τ start_POSTSUBSCRIPT max end_POSTSUBSCRIPT ( roman_Λ , italic_γ ) from Eq. (41), the approximate TGI(γ𝛾\gammaitalic_γ) can be written explicitly as

TGIapprox(γ)=(10dB)log[τmax(Λ,γ)τmax(Λ,γ=0)]=(10dB)log(1+γ2)1/3.\text{TGI}_{\text{approx}}(\gamma)=-(10\;\text{dB})\log\Bigg{[}\frac{\tau_{% \text{max}}(\Lambda,\gamma)}{\tau_{\text{max}}(\Lambda,\gamma=0)}\Bigg{]}=(10% \;\text{dB})\log(1+\gamma^{2})^{1/3}.TGI start_POSTSUBSCRIPT approx end_POSTSUBSCRIPT ( italic_γ ) = - ( 10 dB ) roman_log [ divide start_ARG italic_τ start_POSTSUBSCRIPT max end_POSTSUBSCRIPT ( roman_Λ , italic_γ ) end_ARG start_ARG italic_τ start_POSTSUBSCRIPT max end_POSTSUBSCRIPT ( roman_Λ , italic_γ = 0 ) end_ARG ] = ( 10 dB ) roman_log ( 1 + italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 1 / 3 end_POSTSUPERSCRIPT . (42)

In Tab. 1, we illustrate the achievable values of the TGI that can be obtained by employing the PM correlations. These values are plotted in Fig. 5.

References