Thermal screening and critical scaling of quantum energy teleportation in a harmonic chain

AI summary

Overview

Research area: Quantum information and quantum thermodynamics, specifically quantum energy teleportation (QET) in many-body bosonic systems; intersecting with quantum critical phenomena and finite-temperature Gaussian-state methods.

Technical level: Advanced. The paper assumes familiarity with Gaussian states, covariance matrices, coherent-state POVMs, Matsubara (finite-temperature) expansions, and the thermodynamic limit in lattice models.

Scope: A single-author theoretical study that derives exact asymptotic results for the energy that can be extracted by a measurement-and-feedback QET protocol in a one-dimensional harmonic chain, both at finite temperature and at the zero-temperature critical point.

What This Paper Is About

Quantum energy teleportation is a protocol in which a local measurement on one site injects energy into a correlated many-body system, and a conditional local operation on a distant site then extracts usable energy — using only the measurement outcome plus pre-existing quantum correlations. Prior work on this effect in the harmonic chain was mostly numerical and zero-temperature, and reported an apparent critical decay exponent without deriving it. This paper builds a finite-temperature Gaussian formulation of the protocol, reduces the entire single-site protocol to one thermal correlation function, and derives exact analytic laws for how the extractable energy decays with distance in the gapped, thermally screened, and critical regimes.

Key Contributions

1. A finite-temperature Gaussian-state formulation of QET. For translationally invariant Gibbs states of the harmonic chain, the optimized Gaussian measurement-feedback protocol is shown to reduce entirely to thermal two-point correlators. For an arbitrary Gaussian block measurement of m = 2ℓ + 1 sites, the optimized extracted energy takes the closed form in Eq. (21), combining momentum and position sectors.

2. Reduction of the single-site protocol to a single correlator. The paper proves the identity J_q(β, α; d) = h_d(β, α), valid for α < 1 and 0 < β ≤ ∞, which collapses the full distance dependence of the single-site protocol onto one thermal momentum correlator h_d. This makes the thermodynamic and near-critical limits analytically tractable.

3. Exact asymptotic laws in both limiting regimes. At fixed finite temperature the extracted energy is exponentially screened, controlled by an explicit thermal correlation length. In the zero-temperature critical limit taken after the thermodynamic limit, the paper derives the exact closed-form coupling and the asymptotic law E_ext ∼ d⁻⁴.

4. Resolution of the crossover and of squeezed-measurement effects. Numerical diagnostics explain the apparent drift of the effective decay exponent at intermediate distances, and an analysis of squeezed Gaussian POVMs shows that squeezing changes only the extraction prefactor, not the large-distance scaling class.

Main Findings

• Protocol collapses to thermal correlators: In the single-site setting, E_ext(β, α; d) is given by two terms built from h_d, J_q, h_0, and g_0, together with the measurement noise Ω/2 and 1/(2Ω). The usable Alice–Bob correlations are weighted by the inverse covariance of Alice's noisy measurement record.

• Key identity: J_q(β, α; d) = h_d(β, α) holds for all α < 1 and 0 < β ≤ ∞, before any critical limit, because inserting the integral form of g_n into the definition of J_q reproduces the integral for h_d exactly.

• Thermal screening at finite temperature: In the ordered critical limit α ↑ 1 at fixed β, g_0 diverges as 1/(β√(1−α²)) so the position sector of E_ext vanishes, leaving h_d^crit(β) = −(2/β) Σ_m tanh(κ_m) exp(−2d κ_m) with κ_m = arsinh(√2 π m / β). The m = 1 term dominates at large d.

• Thermal correlation length: The screening length for h_d is ξ_β = [2 arsinh(√2 π / β)]⁻¹. Because E_ext is quadratic in this correlator, the extracted energy decays on the shorter length ξ_E = ξ_β/2 = (4κ_1)⁻¹, with E_ext^crit(β; d) ∼ [2 tanh²κ_1 / (β²[h_0^crit(β) + Ω/2])] exp(−4d κ_1).

• Exact zero-temperature critical result: With β → ∞ first and then α ↑ 1, the coupling is h_d^crit(∞) = −√2 / (π(4d² − 1)) for d ≥ 1, with h_0^crit(∞) = √2/π. This gives the closed-form E_ext^crit(∞; d) = 1 / [π²(4d² − 1)²(√2/π + Ω/2)].

• The d⁻⁴ law: Expanding the closed form at large separation yields E_ext^crit(∞; d) ∼ [1 / (16π²(√2/π + Ω/2))] d⁻⁴, providing the analytic derivation that earlier numerical work lacked.

• Benchmarks against numerics: Figure 3 compares the analytic α = 1 curves with a finite chain of N = 2¹⁵ sites. Near-critical data at α = 1 − 10⁻⁷ track the critical prediction closely over the full distance range, while gapped data at α = 0.98 decay far more steeply. Finite-temperature QET at criticality can exceed gapped ground-state QET over a broad window of separations, showing that proximity to criticality can outweigh the thermal penalty.

• Crossover diagnostic: Three regimes are compared with N = 2000: near-critical low temperature (α = 1 − 10⁻⁹, β = 10⁵), near-critical finite temperature (α = 1 − 10⁻⁹, β = 70), and gapped ground state (α = 0.99, β = 10⁵). The algebraic case gives an effective logarithmic slope η_eff ≈ 4 and local attenuation rate κ_eff ≈ 4/n_B → 0; an exponential envelope instead gives η_eff ≈ n_B/ξ_E and κ_eff ≈ 1/ξ_E. The growing η_eff and finite κ_eff in the screened and gapped cases directly measure the finite length that cuts off critical correlations.

• Two competing length scales: Finite temperature introduces the Matsubara cutoff ξ_β = ξ_h in the critical limit, while a nonzero gap gives ξ_Δ = [arcosh(1/α)]⁻¹ ≈ v/Δ_gap with Δ_gap = √(1 − α) and v = √(α/2). The gapped ground-state behavior is qualitatively consistent with the general exponential decaying bound on QET extraction in finite-range gapped lattice systems.

• Squeezing changes only the prefactor: For q-squeezed measurements the added noises are ν_X = e^(−2r)/(2Ω) and ν_P = Ω e^(2r)/2; for p-squeezed measurements they are ν_X = e^(2r)/(2Ω) and ν_P = Ω e^(−2r)/2. In the critical limit the q-sector drops out, so squeezing acts only by replacing Ω/2 → ν_P. Hence p-squeezing enhances the extracted energy and q-squeezing suppresses it, while Figs. 5 and 6 show both approach the same asymptotic exponent close to −4.

Methodology in Plain English

The author starts from a chain of masses on springs with an onsite pinning potential, written in dimensionless position and momentum variables. A single parameter α controls the ratio of intersite coupling to pinning: 0 < α < 1 is a gapped chain, and α = 1 is the point where the gap closes. The same quadratic Hamiltonian, after a canonical rescaling, describes a chain of coupled linear microwave resonators, so the model doubles as a circuit-QED effective theory.

The initial state is a thermal Gibbs state. Because the Hamiltonian is quadratic, this state is Gaussian, meaning it is fully described by its two-point position and momentum correlations, which can be written as integrals over the normal-mode dispersion ω_α(θ) = √(1 − α cos θ).

Alice performs a Gaussian (coherent-state or squeezed coherent-state) measurement on a block of sites and sends the classical outcomes to Bob, who applies a conditional displacement to his site. Averaging over outcomes gives the post-measurement energy of Bob's local Hamiltonian as a quadratic function of Bob's feedback parameters. Minimizing this quadratic form yields the optimized extracted energy: a ratio of squared correlations to measurement-noise-dressed variances.

The central technical move is to specialize to a single measured site, where these matrices collapse to scalars and the whole protocol depends on one function h_d. That function is then evaluated using a Mittag-Leffler expansion of the thermal occupation factor into bosonic Matsubara frequencies, turning the cosine integrals into rational integrals that can be done by contour methods. This produces a series of decaying terms whose longest-range piece defines the thermal correlation length. Care is taken with the order of limits: α < 1 is held first, the thermodynamic limit is taken, the temperature regime is chosen, and only then is the gap closed at α ↑ 1 — since setting α = 1 at finite N would make the zero mode singular.

Why This Matters

Impact on research: The work moves QET in the harmonic chain from numerical observation to analytic control. It supplies the derivation of the critical decay exponent that earlier ground-state numerics could only estimate, and it separates two physically distinct screening mechanisms — thermal (Matsubara) screening and gap-induced screening — within one framework. It also shows that measurement design (squeezing) is an efficiency knob, not a scaling knob, which sharpens what one should expect from optimizing protocols.

Potential applications (as motivated by the paper's framing, not demonstrated experimentally here):

• Benchmarking bosonic QET implementations in circuit-QED architectures, since the same quadratic Hamiltonian describes coupled linear microwave resonators.
• Characterizing thermal and critical correlations in arrays of coupled resonators or phononic lattices.
• Quantum thermodynamics of many-body devices, where passive mixed states can still hold extractable resources via local measurements.
• Using QET-type protocols as operational probes of proximity to a critical point in Gaussian many-body systems.

Industry relevance: Superconducting-circuit and microwave-resonator platforms are leading candidates for quantum hardware, and the paper's identification of the chain with a circuit-QED bosonic lattice gives device engineers a concrete analytic benchmark. The result that finite-temperature critical QET can outperform gapped ground-state QET over a broad range of separations is directly relevant to operating hardware at non-zero temperature.

Future Directions

• Finite-temperature block protocols. The general block formula, Eq. (21), is exact but the analytic treatment is carried out mainly for a single measured site; extending the asymptotics to multi-site blocks at finite temperature is a natural next step.
• Multisite feedback on Bob's side. The paper suggests this may improve the extraction prefactor while preserving the same asymptotic scaling.
• Reconciling with earlier numerics. The crossover diagnostics explain how an effective decay exponent can drift at intermediate distances; a systematic comparison with the previously reported apparent d^(−3.6) behavior of ground-state numerics would clarify what finite-size data actually measures.
• Experimental or architectural translation. Because the model maps onto coupled microwave resonators, the derived laws could be tested against circuit-QED implementations, including the role of squeezing in the measurement stage.

Target Audience

Advanced graduate students and researchers in quantum information, quantum thermodynamics, and condensed-matter theory — particularly those working on Gaussian states, many-body QET protocols, or quantum criticality in bosonic lattices. It is also relevant to circuit-QED researchers seeking analytic benchmarks for bosonic lattice dynamics, though readers without a background in covariance-matrix methods and finite-temperature field theory will find the derivations demanding.

Authors’ abstract

We develop a finite-temperature Gaussian-state formulation of quantum energy teleportation in the one-dimensional harmonic chain. For Gibbs states, the optimized measurement-feedback protocol reduces to thermal two-point functions. In the single-site protocol, the extracted energy is governed by a single correlator, which makes the thermodynamic and near-critical limits analytically tractable. At fixed finite temperature, the extracted energy is exponentially screened with distance and is controlled by a thermal correlation length for both critical and non-critical $α$. In the zero-temperature critical limit taken after the thermodynamic limit, we derive the exact asymptotic law $E_{\mathrm{ext}}\sim d^{-4}$. Numerical results confirm both regimes and resolve the crossover responsible for the apparent drift of the effective decay exponent at intermediate distances. We also analyze squeezed Gaussian measurements and show that they modify the extraction prefactor: $p-$squeezing enhances the extracted energy, whereas $q-$squeezing suppresses it, without changing the large-distance scaling.