December 22, 2025
Accurate temperature estimation in the quantum and cryogenic regimes remains a fundamental challenge. Here, we investigate nonequilibrium quantum thermometry using a single-qubit probe coupled to a bosonic bath through noncommuting interaction operators, which unify pure dephasing and dissipative dynamics within a spin–boson model. We show that the interference between these two coupling channels induces strong non-Markovian feedback between populations and coherences, leading to coherence trapping and enhanced thermal sensitivity. Remarkably, by tuning the coupling structure, the probe’s temperature sensitivity exhibits a quadratic low-temperature scaling, even under weak coupling. Moreover, while coherence-based measurements are formally suboptimal, they become the most informative in the early nonequilibrium regime, where memory effects dominate. Our findings identify noncommutative system–bath couplings as a practical and tunable resource for achieving high-precision quantum thermometry in realistic open-system architectures.
High-precision temperature estimation at the nanoscale holds major importance, ranging applications from fundamental natural sciences to the rapidly developing field of quantum technologies [1]–[3]. Indeed, recent advances in microscopic physics shed light on quantum thermodynamics in which the precise control and measurement of thermodynamic quantities at the quantum level necessitates quantum-based methods, since classical approaches frequently prove ineffective [4]–[11]. In this inspiration, the so-called quantum thermometry aims to accurately determine the temperature of quantum systems. Indeed, it represents a general strategy for probing the temperature of a quantum reservoir via placing a quantum thermometer in thermal contact with it. Hence, according to their interaction, one can encode thermal information into the state of the thermometer. Then, the temperature can be inferred either from the equilibrium thermal state [12]–[14] or from the nonequilibrium dynamical state [15]–[19] of the thermometer by measuring suitable observables.
Quantum thermometry can surpass classical approaches by exploiting quantum coherence [20]–[22], strong coupling [23]–[25], quantum correlations [26], [27], and periodic driving [28]–[30]. However, a major challenge in most schemes reflect the sensing error typically diverges as the temperature decreases. Indeed, strategies such as strong coupling, periodic driving, and finite measurement resolution can decrease this divergence [31], where reliable performance in the ultra-low temperature regime remains elusive. In the thermalized case, universal results link thermometric precision to the heat capacity via the Gibbs ensemble. But, for very low temperatures, quantum probes often fail to thermalize and instead reach model-dependent nonthermal steady states. For equilibrium probes, the estimation error diverges exponentially as \(T \to 0\), whereas nonequilibrium probes can reduce this to a polynomial divergence, though the error can still be substantial.
In this work, we investigate the process of temperature estimation in a bosonic thermal bath using a single-qubit thermometer evolving under non-Markovian dynamics. Our main goal is to elucidate how the structure of the system–bath interaction, implemented through noncommuting coupling operators, affects thermometric sensitivity. This framework unifies the limiting cases of pure dephasing and purely dissipative interactions and allows one to explore the full range of intermediate couplings, where interference between these channels becomes operational. We show that such noncommutative couplings generate a dynamical feedback between populations and coherences of the probe, leading to long-lived coherence trapping and an enhancement of temperature sensitivity. In particular, we find that in the low-temperature limit, \(T \to 0\), the sensitivity exhibits a quadratic scaling, even under weak coupling, whereas conventional Born–Markov treatments predict an exponential suppression. Finally, we demonstrate that coherence-based measurements, while formally suboptimal, become the most informative at early nonequilibrium regime, regimes that cannot be captured by Markovian master equations [32], where non-Markovian memory effects are most pronounced. This highlights noncommutative couplings as a resource for enhancing the precision and robustness of quantum thermometry in experimentally realistic open-system settings.
The paper is organized as follows. Section 2 describes the model, its dynamical properties, and provides a brief overview of the temperature estimation formalism. Section 3 presents and discusses the results. Section 4 outlines the experimental outlook of our proposed model. Finally, Section 5 summarizes the main conclusions of the paper.
The model— We consider a single qubit, i.e., a thermometer (we will use the terms interchangeably) of the following Hamiltonian: \[\hat{H}_S = \dfrac{\hbar~\varepsilon}{2}\sigma_z,\] where \(\varepsilon\) is the spin splitting. Whereas, the sample to be probed is a bath of a non-interacting bosonic system that rests at thermal equilibrium at a temperature \(T\). Then, the Hamiltonian of the bath reads as: \[\hat{H}_{B}=\sum_{k}\hbar\omega_k \hat{b}_k^{\dagger} \hat{b}_k,\] where \(\omega_k\) is the frequency of the reservoir modes, \(\hat{b}_k^{\dagger}(\hat{b}_k)\) defines the bosonic creation (annihilation) operator for mode \(k\). However, the thermometer is coupled to the sample using the following linear interaction: \[\hat{H}_{int}=\hat{\sigma}_{\alpha}\otimes\sum_{k}\hbar~g_k( \hat{b}_k^{\dagger}+ \hat{b}_k)=\hbar\hat{\sigma}_{\alpha}\otimes B,\] where \(g_k\) denotes the strength coupling of each mode with the thermometer. Moreover, it is common practice to incorporate the frequency dependence of the interaction strengths into the spectral density, which is defined as \(J_{\omega} \equiv \sum_k g_k^2\, \delta(\omega - \omega_k)\) [33], [34]. However, in this work, we adopt an Ohmic spectral density of the following form: \[J_{\omega} = \eta\, \omega\, e^{\omega/\omega_c},\] where \(\eta\) defines a dimensionless coupling constant that carries the order of magnitude of the couplings \(g_k\). While, \(\omega_c\) is the cutoff frequency.
We consider the role of noncommutative coupling [35]–[37] in the interaction between the thermometer and the sample by introducing a parameter \(\alpha\) that defines the structure of the coupling operator. Specifically, we define the coupling operator as: \[\label{int-oper} \hat{\sigma}_{\alpha} = (1 - \alpha)\, \hat{\sigma}_\mathrm{z}+ \alpha\, \hat{\sigma}_\mathrm{x},\tag{1}\] which allows interpolation between purely diagonal and purely off-diagonal interactions in the energy basis of the thermometer. In fact, for \(\alpha = 0\), the coupling reduces to a purely dephasing interaction mediated by \(\hat{\sigma}_z\), which preserves the populations of the thermometer’s energy levels while affecting the coherences. Conversely, for \(\alpha = 1\), the interaction is governed by \(\hat{\sigma}_x\), corresponding to a purely dissipative or energy-exchange coupling. This parametrization enables the exploration of intermediate regimes, where both diagonal and off-diagonal coupling effects coexist. This offers a flexible framework to investigate how the nature of the system–probe interaction influences the performance of quantum thermometry.
Since we are interested only in the dynamics of the thermometer, we employ an exact non-Markovian master equation derived using Zwanzig’s projection-operator method [33], [34], [38] to trace out the sample degrees of freedom from the Liouville equation of motion, retaining terms up to second order in the thermometer–sample coupling. \[\begin{align} \label{NME} \dot{\rho}_{S}&=&i/\hbar~[ \rho_{S}, H_S]\\ &+&\intop_{0}^{t}d\tau\{\Phi_T(t-\tau)[\sigma_{\alpha}(t,\tau)~\rho_{S},\sigma_{\alpha}]+h.c\}\nonumber, \end{align}\tag{2}\] where \(\Phi_T(t)=\langle e^{i H_B t} B e^{-i H_B t} B \rangle_{\rho_{B}}\) is the bath correlation functions, with \(\langle \dots \rangle_{\rho_{B}}=\mathrm{Tr}\{\dots~\rho_{B}\}\). \(\rho_B\) is the thermal equilibrium state of the reservoir at temperature \(T\). Indeed, it takes the following form: \[\rho_B = \prod_k (1 - e^{\beta \omega_k}) e^{-\beta \omega_k b_k^\dagger b_k},\] where \(\beta = \hbar / k_B T\) is the inverse temperature (\(\hbar=k_B=1\)). The time-dependent interaction operator is expressed as \(\hat{\sigma}_{\alpha}(t,\tau) = \alpha ~ \hat{\sigma} _-~ e^{i \varepsilon (t-\tau )}+\alpha ~ \hat{\sigma} _+~ e^{-i \varepsilon (t-\tau )}+(1-\alpha ) ~\hat{\sigma} _z\). Additionally, we assume that the initial state of the thermometer-sample system as \(\rho_{tot}(0) = \rho_{S}(0) \otimes \rho_{B}\), where the initial state of the thermometer is chosen as \(\rho_{S}(0)=\ket{+}\bra{+}\), such that \(\ket{+}=\dfrac{1}{\sqrt{2}}(\ket{0}+\ket{1})\).
It is worth emphasizing that the master equation (2 ) is exact, where no Markov and even no secular approximations have been applied. Hence, the density matrix of the thermometer by means of a single qubit can be conveniently rewritten in the Bloch representation [39] as: \(\rho_S = \frac{1}{2} \left( \mathbb{I}_2 + \mathbf{\Delta} \cdot \hat{\boldsymbol{\sigma}} \right)\), where \(\mathbb{I}_2\) is the \(2 \times 2\) identity matrix. While, \(\mathbf{\Delta} = (\Delta_x, \Delta_y, \Delta_z)\) denotes the real Bloch vector characterizing the state of the qubit and \(\hat{\boldsymbol{\sigma}} = (\hat{\sigma}_x, \hat{\sigma}_y, \hat{\sigma}_z)\) denotes the vector of Pauli matrices. Therefore, one can obtain the generalized Bloch equations for the components of the thermometer’s state as:
\[\begin{align} \tag{3} \dot{\Delta}_x &=& -\varepsilon\, \Delta_y - 4 \alpha (\alpha - 1) G_t - 4 \alpha (\alpha - 1)\, \Delta_z\, \mathcal{K}_t - 4 (\alpha - 1)^2\, \Delta_x\, R_t, \tag{4} \\ \dot{\Delta}_y &=& \Delta_x \left(\varepsilon + 4 \alpha^2 X_t \right) + 4 (\alpha - 1) \alpha (F_t - L_t) - 4 \Delta_y \left(\alpha^2 \mathcal{K}_t + (\alpha - 1)^2 R_t \right) - 4 (\alpha - 1) \alpha\, \Delta_z\, X_t, \tag{5}\\ \dot{\Delta}_z &=& -4 \alpha^2 G_t - 4 \alpha^2 \Delta_z\, \mathcal{K}_t - 4 (\alpha - 1) \alpha\, \Delta_x\, R_t, \tag{6} \end{align}\]
where the time-dependent coefficients in the above equation are expressed as:
\[\begin{align} \label{Coeff} R_t&=&\int_0^{\infty } \frac{J_{\omega} \sin (t \omega ) \coth \left(\frac{\omega }{2 T}\right)}{\omega } \, d\omega,~~~~\mathcal{K}_t=\int_0^{\infty } \frac{J_{\omega} \coth \left(\frac{\omega }{2 T}\right) (\varepsilon \sin (\varepsilon t) \cos (t \omega )-\omega \cos (\varepsilon t) \sin (t \omega ))}{\varepsilon^2-\omega ^2} \, d\omega,\nonumber\\ L_t&=&\int_0^{\infty } \frac{J_{\omega} (1-\cos (t \omega ))}{\omega } \, d\omega,~~~X_t=\int_0^{\infty } \frac{J_{\omega} \coth \left(\frac{\omega }{2 T}\right) (-\omega \sin (\varepsilon t) \sin (t \omega )+\varepsilon (-\cos (\varepsilon t)) \cos (t \omega )+\varepsilon)}{\varepsilon^2-\omega ^2} \, d\omega,\nonumber\\ F_t&=&\int_0^{\infty } \frac{J_{\omega} (\varepsilon \sin (\varepsilon t) \sin (t \omega )+\omega \cos (\varepsilon t) \cos (t \omega )-\omega )}{\varepsilon^2-\omega ^2} \, d\omega,~~~G_t=\int_0^{\infty } \frac{J_{\omega} (\omega \sin (\varepsilon t) \cos (t \omega )-\varepsilon \cos (\varepsilon t) \sin (t \omega ))}{\varepsilon^2-\omega ^2} \, d\omega.\nonumber \end{align}\tag{7}\]
The structure of the generalized Bloch equations Eqs. (4 –6 ) provides a transparent physical picture of the thermometer’s dynamics under noncommutative couplings. The time-dependent coefficients \(R_t, \mathcal{K}_t, X_t, F_t, L_t\), and \(G_t\) represent dynamical memory kernels induced by the structured bath, capturing how past system–bath correlations influence the instantaneous evolution of the probe. Crucially, the terms proportional to \(\alpha(1-\alpha)\) in Eqs. (4 –6 ) are direct manifestations of quantum interference between the dephasing \((\hat{\sigma}_\mathrm{z})\) and dissipative \((\hat{\sigma}_\mathrm{x})\) interaction pathways. For intermediate values of \(\alpha\), these noncommutative contributions encode a cross-coupling between the population \((\Delta_z)\) and coherence \((\Delta_x)\) dynamics. Meaning the bath simultaneously monitors the qubit’s energy (a quantum non-demolition process linked to dephasing) while inducing transitions between its levels (a dissipative process). The concurrent action of these two channels establishes a dynamical feedback, where coherences affect population flow and vice versa.
Information backflow— As it has already been mentioned, we aim to employ a non-Markovian master equation to investigate the role of the coupling mixing parameter \(\alpha\) on the thermometric performance. Indeed, we shall study the influence of \(\alpha\) on the emergence of non-Markovian dynamics in this model. In this regard, note that various definitions of non-Markovianity have been proposed in the literature [40]–[44]. In this work, we adopt a coherence-based approach, where non-Markovianity is witnessed through re-coherence [45]–[47], that is, the temporary revival of quantum coherence during the system’s evolution. This revival indicates a reversal in the decoherence process, which gives rise to back-flow of information from the environment to the open system. Accordingly, the presence of re-coherence serves as a qualitative signature of non-Markovian behavior in the dynamics [41]. \[\begin{align} \label{NM} \mathcal{N}_C = \int_{\dot{C}(t) > 0} \dot{C}(t)\,dt, \end{align}\tag{8}\] where \(C(t)\) reflects the quantum coherence [48], [49] of the system, quantified by means of the so-called \(l_1\)-norm of coherence. Note that in Bloch representation, the \(l_1\)-norm of coherence can be expressed as \(C(t) = \sqrt{\Delta_x^2 + \Delta_y^2}\), where \(\Delta_x\) and \(\Delta_y\) are the transverse components of the Bloch vector obtained from Eq. (3 ). In general, if the coherence at \(t=0\) satisfies \(C(0) = 1\), then the state corresponds to a pure state on the equator of the Bloch sphere. However, the re-coherence is identified for distinct intervals when \(\dot{C}(t) > 0\) [46], where one can integrate successfully these revivals over time to obtain the total non-Markovianity witness, namely \(\mathcal{N}\).
Figure (1) illustrates the delicate interplay between the qubit’s coherence dynamics and the degree of non-Markovianity. The top panel displays how the coupling structure, controlled by the parameter \(\alpha\), determines the trajectory of the probe’s coherence on the Bloch sphere. As anticipated, when \(\alpha = 1/2\), the probe maintains quantum coherence for a longer duration and eventually stabilizes at a nonzero steady-state value. In contrast, for the case of \(\alpha = 0\) (dephasing) and \(\alpha = 1\) (dissipative), coherence decays faster to zero. This behavior stems from the interference between diagonal and off-diagonal couplings, encoded in the cross terms proportional to \(\alpha(1-\alpha)\) in Eqs. (4 –6 ). These terms mediate a continuous exchange of information between the population and coherence sectors, generating a feedback mechanism in which the environment both induces transitions and partially restores lost coherence. Consequently, the coexistence of the two coupling channels suppresses decoherence and gives rise to steady-state or trapped coherence, a signature of non-Markovian memory effects [50]–[52]. The bottom panel further supports this interpretation: the non-Markovianity measure \(\mathcal{N}\) exhibits pronounced values that coincide with the re-coherence regions in the dynamics. The residual long-time coherence \(|\Delta_x(\infty)|\) thus quantifies the persistent memory retained by the probe, indicating that noncommutative coupling can effectively stabilize coherence against environmental noise by exploiting bath-induced backflow and interference effects.
Thermometric performance— It has been shown that quantum sensing of the sample’s temperature is fundamentally bounded by means of the quantum Cramér–Rao bound (QCRB), which sets the minimum achievable uncertainty in the estimation of any unknown parameter, temperature in our case. The QCRB is given by [53], [54] \[\begin{align} \label{QCRB} \delta T \geq \frac{1}{\sqrt{M \mathcal{F}_Q(T)}}, \end{align}\tag{9}\] where \(\delta T\) denotes the standard deviation in the temperature estimate, \(M\) is the number of independent measurements and \(\mathcal{F}_Q(T)\) defines the quantum Fisher information (QFI) associated with the parameter \(T\). Therefore, maximizing the QFI plays a central role in optimizing the sensitivity of quantum thermometry. The QFI is defined as \(\mathcal{F}_Q(T) = \mathrm{Tr} \{ \hat{L}_T^2 \, \rho_T \}\) [54], where the symmetric logarithmic derivative (SLD) operator \(\hat{L}_T\) satisfies \(\partial_T \rho_T = \frac{1}{2} ( \hat{L}_T \rho_T + \rho_T \hat{L}_T ).\) This quantity quantifies the maximum information that can be extracted about the temperature \(T\) from the quantum state \(\rho_T\).
For the adopted model and by using the master equation given in Eq. (2 ), one can compute the QFI with respect to the temperature \(T\) as \(\mathcal{F}_Q = \left( \frac{\partial \vec{\mathbf{\Delta}}}{\partial T} \right)^{\!T} \cdot \mathbf{M}^{-1} \cdot \left( \frac{\partial \vec{\mathbf{\Delta}}}{\partial T} \right)\), where \(\mathbf{M}^{-1}\) defines the inverse of the QFI metric tensor for a qubit, that is: \(\mathbf{M}^{-1} = \mathbb{I}_3 + \frac{\vec{\mathbf{\Delta}} \vec{\mathbf{\Delta}}^T}{1 - \lvert \vec{\mathbf{\Delta}} \rvert^2}\), such that \(\mathbb{I}_3\) is the \(3 \times 3\) identity matrix.
Note that when performing the numerical calculation of the QFI, one needs to evaluate the first-order derivatives of \(\partial_T \boldsymbol{\Delta}\). In this paper, the derivative of an arbitrary temperature-dependent function \(\boldsymbol{f}(T)\) is computed numerically using the following finite-difference method [55]: \[\begin{align} \frac{d\boldsymbol{f(T)}}{dT} \approx \frac{-\boldsymbol{f}(T + 2\delta) + 8\boldsymbol{f}(T + \delta) - 8\boldsymbol{f}(T - \delta) + \boldsymbol{f}(T - 2\delta)}{12\delta},\nonumber\\ \end{align}\] where \(\delta = 10^{-7} \times T\), which provides a good accuracy.
The central aim of this work is to investigate how the structure of the probe–sample interaction, characterized by the mixing parameter \(\alpha\) (Eq. 1 ), influences the performance of quantum thermometry in a nonequilibrium setting. Indeed, by varying \(\alpha\), we interpolate between purely dephasing thermometer [56]–[58] and purely dissipative thermometer [30], [59]–[61], thereby exploring a broad class of hybrid interactions that combine coherence-preserving and energy-exchange effects. Our main goal is to evaluate the QFI associated with the probe’s state. In what follows, we analyze the behavior of QFI as a function of \(\alpha\) at various fixed evolution times. This allows us to identify how nonequilibrium effects and the nature of the coupling jointly affect the sensitivity of the probe.
Figure (2) exhibits the QFI as a function of the mixing parameter \(\alpha\) for various evolution times \(\omega_c t\). The top panel highlights the behavior at short to intermediate times. At a short interval of time, the temperature information is predominantly encoded in the nonequilibrium state of the probe via pure dephasing interactions (\(\alpha \to 0\)). However, as \(\alpha\) increases toward unity, the sensitivity of the probe to temperature variations decreases, indicating that the off-diagonal interactions are less effective in this regime. However, at intermediate times, increasing \(\alpha\) leads to a noticeable improvement in sensitivity, up to a saturation point where the probe reaches thermal equilibrium with the environment, and no further information can be extracted.
In contrast, the lower panel of Fig. 2 reveals that the QFI exhibits a clear enhancement at intermediate values of \(\alpha\), corresponding to noncommutative coupling configurations. This behavior indicates that the coexistence of dephasing and dissipative channels facilitates a more efficient extraction of temperature information from the bath. This enhancement originates from the interference between population and coherence dynamics, which allows the probe to simultaneously access both the energy exchange and phase fluctuation sectors of the bath. In this regime, the system operates far from thermal equilibrium, and the backflow of information characteristic of non-Markovian dynamics enables partial recovery of lost sensitivity during evolution. Such nonequilibrium feedback amplifies the temperature dependence of the populations and coherences, resulting in improved thermometric precision beyond the limits attainable in purely dephasing or purely dissipative interactions. This mechanism is reminiscent of recent findings in open quantum transport [36], where noncommutative couplings and structured environments in a generalized spin-boson model are shown to significantly boost energy transport and thermal rectificationvia purely quantum effects. Therefore, the observed enhancement in QFI at intermediate \(\alpha\) reflects an interplay between coherence preservation and energy exchange, establishing noncommutative coupling as a resource for nonequilibrium quantum thermometry.
In what follows, we focus on the early-time evolution of the probe and the scaling behavior at low temperatures. Measuring the temperature of a sample in the low-temperature regime is particularly challenging due to the vanishing of the QFI, which leads to a diverging quantum Cramér-Rao bound, i.e., \(\lim_{T \to 0} \mathcal{F}_Q(T) = 0\). However, under the assumptions of weak coupling and the Born-Markov approximation, it has been shown that the QCRB diverges exponentially at low temperatures as [31], [62] \[\delta T^2 \geq \frac{2 T^4 e^{\varepsilon /T}}{M \varepsilon ^2}.\] The above expression clearly indicates that for \(T \to 0\), the uncertainty in estimating the temperature increases rapidly. These limitations motivate the exploration of beyond-Markovian dynamics, where memory effects can alter the scaling behavior [63]. In particular, we will examine how our non-Markovian master equation can enhance the QFI in the low-temperature regime.
The top panel of Fig. (3) shows the non-equilibrium QFI versus the sample temperature for different probing times, in the presence of both diagonal and off-diagonal couplings (\(\alpha=1/2\)). We observe that non-Markovian dynamics can enhance the low-temperature scaling, where a non-vanishing QFI persists as \(T \to 0\) and scales as \(\mathcal{F}_Q \propto T^2\), in contrast to the Markovian case, where the QFI vanishes at low temperatures (\(T \to 0\), black-dotted line). Interestingly enough, this favorable scaling behavior is achieved at early times [63] and under weak coupling, i.e., through nondestructive probing. The \(\mathcal{F}_Q\propto T^2\) low-temperature scaling has already been proven in a system of a fermionic tight-binding chain with access to only two sites [31], and for a quantum Brownian probe strongly coupled to an infinite bath of harmonic oscillators [64]. It is worth mentioning that in this setup, the low-temperature scaling of the QFI also holds for a gapless thermometer. In Ref. [65], it has been shown that for \(\varepsilon = 0\) and \(T/\omega_c \ll 1\), the QFI diverges as \(\mathcal{F}_Q \propto 1/T^2\).
Now, let us examine the Fisher information to identify the (sub)optimal thermometer observable \(\hat{\mathcal{O}}\) with the associated temperature uncertainty given in [54], [66] \[\delta T(\hat{\mathcal{O}}) \geq \frac{1}{\sqrt{M \mathcal{F}_C(\hat{\mathcal{O}})}},\] where \(\mathcal{F}_C(\hat{\mathcal{O}})\) reflects the Fisher information associated with the performed measurement on the observable \(\hat{\mathcal{O}}\). The Fisher information serves as a precision quantifier for a specific measurement. Whereas the QFI determines the fundamental lower bound on the estimation error, regardless of the measurement performed. For a two-level system, the Fisher information associated with a measurement reads [67] \[\mathcal{F}_C(\hat{\mathcal{O}})=\dfrac{1}{\braket{\Delta \hat{\mathcal{O}}^2}}\bigg(\dfrac{\partial\braket{\hat{\mathcal{O}}}}{\partial T}\bigg)^2,\] where \(\braket{\hat{\mathcal{O}}}\) and \(\braket{\Delta \hat{\mathcal{O}}^2}\) represent the mean and variance of the measured observable \(\hat{\mathcal{O}}\). In this work, we will focus on comparing measurements of coherences and populations-difference, i.e., measuring \(\hat{\sigma}_\mathrm{x}\) and \(\hat{\sigma}_\mathrm{z}\), receptively. Thus, Fisher information reads \(\mathcal{F}_C(\hat{\sigma}_\mathrm{x})=(\partial_T \Delta_x)^2/(1 - \Delta_x^2)\) and for populations-difference (\(\hat{\sigma}_\mathrm{z}\)), we simple replace \(\Delta_x\) with \(\Delta_z\).
In Fig. 3 (bottom), we compare the Fisher information obtained from coherence-based (solid) and population-based (dashed) measurements in the low-temperature and short-time regimes. Both observables exhibit a quadratic scaling
as \(T \to 0\), consistent with the behavior of the QFI. The superiority of coherence-based measurements at early times arises because quantum coherences encode the immediate response of the probe to temperature-induced
phase fluctuations, before dissipative energy exchange with the bath becomes dominant. In this transient regime, off-diagonal elements of the density matrix carry direct information about the noise spectrum and its temperature dependence, allowing a fast
and sensitive estimation process. In contrast, population-based observables respond on longer relaxation timescales, reflecting slower thermalization dynamics.
Moreover, short-time non-equilibrium thermometry can be both highly informative and minimally invasive [57], since
measurements performed before full equilibration disturb the system less while still capturing the essential thermal signatures. In this sense, the enhanced sensitivity of coherence-based measurements observed here not only reflects their access to fast
non-Markovian correlations but also aligns with the principle of optimal, non-invasive quantum thermometry, where precision and minimal disturbance coexist in the early-time regime. This highlights the operational relevance of coherence measurements and
provides a practical route to experimentally feasible, nonequilibrium quantum thermometers.
We note that at sufficiently low temperatures, the (quantum) Fisher information deviates from the Markovian result, even for weak coupling. This highlights that faithfully extracting temperature in this regime requires accounting for system–environment correlations, which are absent in simple dissipation models such as those described by the GKSL master equation [68], [69].
From an experimental perspective, the coupling structure embodied in \(\hat{\sigma}_{\alpha}\) is not purely theoretical. For example, in superconducting-circuit architectures an inductively shunted transmon qubit has been demonstrated with flux-tunable switching between transverse (\(\hat{\sigma}_\mathrm{x}\)) and longitudinal (\(\hat{\sigma}_\mathrm{z}\)) coupling to a harmonic mode [70]. Moreover, hybrid spin-photon systems have achieved parametric longitudinal coupling (\(\hat{\sigma}_\mathrm{z}\)) between a semiconductor spin qubit and a high-impedance resonator, underscoring the feasibility of engineered qubit–bath interactions beyond standard transverse (\(\hat{\sigma}_\mathrm{x}\)) coupling [71]. In such platforms, the coupling strength and orientation can be tuned through drive amplitude, flux bias, or coupling-axis control, and qubit state tomography can be performed [72] to extract both populations and coherences. By sweeping the sample temperature in the cryogenic regime [73] and comparing the achieved temperature estimation error for different coupling configurations, one could directly test the Fisher-information enhancement and the influence of non-Markovian dynamics.
The present work introduces a tunable system–bath coupling operator \(\hat{\sigma}_{\alpha} = (1 - \alpha) \hat{\sigma}_z + \alpha \hat{\sigma}_x\) that continuously interpolates between pure dephasing and pure dissipation, thereby unifying the two major coupling paradigms within quantum thermometry. This unified treatment reveals qualitatively new physical insights: in particular, we demonstrate that intermediate values of \(\alpha\) enable enhanced temperature sensitivity, manifested via non-Markovian memory effects and a low-temperature scaling of the QFI that sharply contrasts with the exponential suppression predicted under standard Born–Markov approaches. By revealing optimal coupling structures for thermometric precision, our work goes beyond formal extension of prior models and uncovers controllable resources (non-commutativity of couplings + non-Markovianity) for quantum sensing.
Moreover, we analyzed the scaling behavior of the QFI at low temperatures and short times for the co-existing diagonal and off-diagonal couplings. In this regime, we concluded that \(\mathcal{F}_Q\propto T^2\), which is consistent with fundamental thermodynamic limits [31], [63], [64]. At sufficiently low temperatures, deviations from the Markovian prediction appeared even for the weak coupling regime. Further, we showed that in this regime a measurement based on the coherence [74], [75] of the thermometer is suboptimal to extract thermal information from the sample.
Beyond the theoretical implications of our model, the present framework provides a physically grounded route toward experimentally relevant quantum thermometry. The generalized noncommutative coupling operator can be realized in controllable settings where both diagonal and off-diagonal noise channels are routinely engineered. In such environments, the interference between dephasing and dissipative processes identified in this work could be tuned through external controls [70] enabling the observation of coherence trapping and enhanced thermometric sensitivity in nonequilibrium regimes. These findings indicate that noncommutative system–bath couplings represent a versatile resource for improving temperature estimation in realistic quantum devices operating at cryogenic or mesoscopic scales [73].