Finite-frequency fluctuation-response bounds for open quantum systems

Abstract: We derive a finite-frequency fluctuation-response inequality for Markovian open quantum systems in an input-output setting. For any downstream measurement of the emitted field, the measured lock-in response-to-noise matrix is bounded by the output-field quantum Fisher information rate. For dissipative amplitude modulation with vacuum inputs, this information rate is further bounded by a frequency-independent signal-channel activity, which reduces for kinetic modulation to the stationary channel fluxes. The result is detector-facing but unraveling-independent: it applies after choosing a measurement record, while the information ceiling is set by the quantum field before any detection scheme or trajectory representation is selected. We formulate the bound for multiple signal channels and real finite-frequency quadratures, and illustrate it with a single-sided cavity, resonance fluorescence, and a truncated Kerr-parametric cat resonator.

• The paper establishes a finite-frequency inequality showing that detector response-to-noise precision is bounded by the emitted field’s quantum Fisher information and, for dissipative amplitude signals, by calibrated channel activity.
• The methodology combines lock-in response statistics, output-field data processing, Liouvillian resolvents, and QFI monotonicity to produce detector-facing yet unraveling-independent bounds across homodyne, heterodyne, and counting measurements.
• The results recover classical finite-frequency fluctuation-response inequalities and are validated for cavities, resonance fluorescence, and a Kerr-parametric cat resonator, while leaving Hamiltonian signals, non-Markovian baths, and incompatible multiparameter estimation open.

Overview

This paper establishes a finite-frequency fluctuation-response inequality (FRI) for Markovian open quantum systems formulated at the quantum input-output level. The central result is a chain of inequalities,

$$\mathsf R^{T}(\omega)\,[\mathsf S^{out}(\omega)]^{+}\,\mathsf R(\omega)\ \preceq\ F^{Q}_{out}(\omega)\ \preceq\ A_{\rm sig}\otimes I_2,$$

where $\mathsf R(\omega)$ is the lock-in response matrix of measured output currents, $\mathsf S^{out}(\omega)$ is the corresponding noise covariance matrix, $F^{Q}_{out}(\omega)$ is the quantum Fisher information (QFI) rate carried by the emitted output field at frequency $\omega$, and $A_{\rm sig}$ is a signal-channel activity built from calibrated dissipative coupling tangents and stationary channel fluxes. The structural novelty is the placement of the bound: the left-hand side involves only detector-facing classical spectral data, while the ceiling is set by the output field before any detection scheme or trajectory unraveling is chosen. The paper describes this as "detector-facing but unraveling-independent."

Setup and main theorems

The system is finite-dimensional with GKSL dynamics $\dot\rho = -i[H,\rho] + \sum_\mu D[L_\mu]\rho$, coupled to Markovian bosonic input channels via the Gardiner–Collett relation $b_{out} = b_{in} + L$. A key standing assumption is exponential mixing: the Liouvillian spectrum restricted to the traceless subspace lies in $\{z: \operatorname{Re} z \leq -\gamma_{\rm mix}\}$. This guarantees invertibility of $-i\omega - L$ on that subspace for all real nonzero $\mathsf R(\omega)$0, so that stationary spectra and Liouvillian resolvents are unambiguous.

The signal is modeled as weak sinusoidal modulation of coupling amplitudes, $\mathsf R(\omega)$1, expanded in unit-RMS cosine/sine modes at fixed frequency. Two theorems compose into the main chain:

Data-processing bound: for any downstream measurement producing currents, the measured response-to-noise matrix satisfies $\mathsf R(\omega)$2. The proof combines a Schur-complement/score-identity argument on the lock-in statistics with monotonicity of QFI under the POVM implementing the detector.

Activity bound: for purely dissipative amplitude tangents satisfying $\mathsf R(\omega)$3 with vacuum inputs, the directional QFI rate obeys $\mathsf R(\omega)$4, where $\mathsf R(\omega)$5. For kinetic modulation $\mathsf R(\omega)$6, this reduces to weighted stationary channel fluxes $\mathsf R(\omega)$7 — photon fluxes, jump rates, or tunneling rates depending on platform — making the right-hand side experimentally calibratable and frequency-independent.

A companion coherent-input version bounds displacement sensing by the input QFI rate itself: $\mathsf R(\omega)$8, since a displaced vacuum mode carries QFI rate exactly $\mathsf R(\omega)$9 per real quadrature in this normalization.

Operational interpretation

Three points deserve emphasis. First, the bound is not an operator fluctuation-dissipation theorem; it relates measured output-current spectra to output-field information rather than internal commutators to symmetrized correlators. Second, it applies uniformly across homodyne, heterodyne, photon-counting, inefficient, and adaptive detection, because data processing removes any dependence on the chosen POVM. Third, noncommutativity enters through the Lindblad evolution and the input-output relation even though the final record is classical — homodyne spectra of a driven qubit contain phase-sensitive structure absent from any classical jump process with identical mean rate. The authors note that the derivation requires care with normalization conventions: omitting the $\mathsf S^{out}(\omega)$0 factors in the real lock-in vector produces a spurious factor-of-two mismatch between complex Fourier and real-mode Fisher information conventions.

Examples

The single-sided cavity provides an analytically solvable Gaussian benchmark saturating the coherent-input bound. Because lossless single-port scattering preserves quadrature vacuum noise ($\mathsf S^{out}(\omega)$1, $\mathsf S^{out}(\omega)$2), phase-matched homodyne detection achieves $\mathsf S^{out}(\omega)$3 at every frequency. The implication is sharp: passive linear scattering can reshape and delay signal information but cannot amplify it beyond what the input tone carries.

For resonance fluorescence, kinetic modulation of the radiative coupling yields activity $\mathsf S^{out}(\omega)$4 equal to the steady fluorescence flux. On resonance ($\mathsf S^{out}(\omega)$5), the paper verifies the scalar bound by explicit substitution: $\mathsf S^{out}(\omega)$6 evaluates to a manifestly nonnegative rational function of frequency, confirming $\mathsf S^{out}(\omega)$7 for every homodyne phase, with equality approached only in the trivial undriven limit.

The truncated Kerr-parametric cat resonator validates the full multiparameter matrix inequality numerically. With two dissipative signals (external and internal loss channels), a $\mathsf S^{out}(\omega)$8 response matrix, and cutoff $\mathsf S^{out}(\omega)$9, the largest eigenvalue of the normalized matrix $F^{Q}_{out}(\omega)$0 remains below unity across the sampled frequency window. Notably, the theorem is invoked only at each finite cutoff; no convergence analysis of the infinite-dimensional limit is provided.

Relation to existing bounds

In the classical counting limit — diagonal Lindblad jump operators $F^{Q}_{out}(\omega)$1 with ideally monitored channels — the main inequality reproduces Dechant's finite-frequency FRI for Markov jump processes exactly, not merely analogously: the signal activity becomes the standard activity matrix $F^{Q}_{out}(\omega)$2. Relative to response kinetic uncertainty relations [Liu–Gu], the present work implements the same hierarchy (response precision ≤ Fisher information ≤ activity) at the level of the emitted field, extending coverage to finite frequencies and phase-sensitive measurements. Relative to quantum-trajectory fluctuation-response bounds [Van Vu], the key distinction is that trajectory-level inequalities presuppose a chosen unraveling, whereas here different detectors are different POVMs on one field and the ceiling precedes that choice.

Limitations and open questions

Several restrictions are stated explicitly. The activity bound is special to purely dissipative amplitude tangents; without the condition $F^{Q}_{out}(\omega)$3, a coherent channel tangent appears in sequential channel QFI and is not bounded by jump activity alone. Hamiltonian perturbations $F^{Q}_{out}(\omega)$4 lie outside the activity formalism; the data-processing half should still apply, but the appropriate information cost would be a Hamiltonian-tangent "quantum Fisher strength" not reducing to channel fluxes. The frequency independence of $F^{Q}_{out}(\omega)$5 relies on Markovianity — non-Markovian baths or pre-filtered signals generally yield frequency-dependent ceilings. The result is first-order in the perturbation; nonlinear response requires higher-order information inequalities. In multiparameter settings, the SLD-QFI matrix need not be jointly attainable when parameters are incompatible, so the matrix inequality must be read as an upper bound rather than a jointly saturable estimation bound; a finite-frequency Holevo-type formulation is identified as the needed refinement. Finally, tightness is architecture-dependent: unobserved losses, inefficiency, thermal noise, and internal nonlinearity generically make the inequality strict, and exact saturation should not be expected outside linear passive systems.

Conclusion

The paper formulates a finite-frequency fluctuation-response inequality whose left-hand side consists entirely of measurable output-current spectra and lock-in responses, bounded above by the output-field QFI rate and, for calibrated dissipative amplitude modulation, by a frequency-independent signal-channel activity. Its distinctive contribution is enforcing this constraint before unraveling selection while retaining direct experimental interpretability through calibratable fluxes. Whether the structure extends to Hamiltonian signals, non-Markovian environments, and incompatible multiparameter estimation at finite frequency remains open.