Parametric Feedback Cooling

• Parametric feedback cooling is an active method that modulates trap stiffness at twice the resonance frequency to add damping and lower effective temperature.
• It has been applied in levitated optomechanics, cavity QED, and pendulum systems, achieving significant temperature reductions and high quality factors.
• Advanced implementations combine digital feedback, adaptive phase control, and optimal estimation to enhance cooling efficiency and approach quantum limits.

Searching arXiv for the provided ids to ground the article in the cited literature. Parametric feedback cooling is an active cooling method for mechanical motion in which the stiffness of a trap or resonator is modulated at approximately twice the mechanical resonance frequency so that the modulation extracts energy rather than injects it. In its standard form, the method acts on the center-of-mass or librational motion of a trapped object by engineering an additional damping term through phase-controlled modulation of the restoring potential, with the effective temperature reduced according to the ratio of intrinsic to total damping (Gieseler et al., 2012). The technique has been developed in levitated optomechanics, cavity QED, rotational nanomechanics, and low-frequency pendulum systems, with implementations ranging from analog sinusoidal modulation to square-wave digital control, phase-adaptive protocols, and measurement-based optimal quantum feedback (Vovrosh et al., 2016).

1. Fundamental mechanism

For a harmonically trapped degree of freedom $x(t)$, the uncooled dynamics are commonly written as

$$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$

with $F_{\rm th}(t)$ a Langevin force satisfying $$\langle F_{\rm th}(t)F_{\rm th}(t')\rangle = 2m\Gamma_0 k_B T_0 \delta(t-t')$$ (Gieseler et al., 2012). Parametric feedback introduces a time-dependent spring constant $k(t)=k_0[1+\epsilon(t)]$, so that the oscillator equation becomes

$$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2[1+\epsilon(t)]x = F_{\rm th}(t).$$

A convenient choice is $\epsilon(t)=\epsilon_0\cos(2\Omega_0 t+\phi)$, namely stiffness modulation at $2\Omega_0$ with a controllable phase $\phi$ (Gieseler et al., 2012).

Under a rotating-wave or slow-amplitude analysis, the $2\Omega_0$ modulation yields an effective friction term $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$0, with

$$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$1

If $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$2, the modulation adds damping; if $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$3, it amplifies the motion (Gieseler et al., 2012). Closely related formulations appear in cavityless levitated optomechanics, where the control signal $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$4 modulates the trap frequency in the Hamiltonian $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$5, and in librational systems, where the restoring potential is written as $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$6 (Ferialdi et al., 2019); (Gao et al., 2024).

The physical picture is consistent across platforms. When the particle is near a turning point, the trap is stiffened; when it passes through the origin, the trap is softened. If the modulation is kept in the appropriate phase relation with the motion, the net effect is energy extraction over each cycle rather than parametric excitation (Penny et al., 2021). In experimental practice, this phase relation can be implemented by using $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$7, by PLL-based phase tracking, or by continuous quadrature estimation (Vovrosh et al., 2016); (Penny et al., 2021).

2. Dynamical descriptions and temperature reduction

The basic thermal consequence of the added damping is expressed by

$$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$8

so that the center-of-mass temperature is reduced by increasing the parametric damping while the thermal force remains tied to the intrinsic damping $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$9 (Gieseler et al., 2012). The same structure appears in several analyses. In the direct comparison study of levitated oscillators, the effective parametric damping rate is identified as $F_{\rm th}(t)$0, giving

$F_{\rm th}(t)$1

when detection noise and back-action are neglected in the analytic model (Penny et al., 2021). In mechanical pendulum cooling, the corresponding extra damping is

$F_{\rm th}(t)$2

derived from a suspension-point modulation at $F_{\rm th}(t)$3 (Hartwig et al., 2020).

The frequency-domain description is equally central. For translational motion of a nanoparticle, the displacement power spectral density can be written as

$F_{\rm th}(t)$4

so that the spectral area is proportional to $F_{\rm th}(t)$5 (Gieseler et al., 2012). In the parabolic-mirror-trap implementation, the spectral density is written as a shifted, broadened Lorentzian,

$F_{\rm th}(t)$6

with $F_{\rm th}(t)$7 and $F_{\rm th}(t)$8, thereby making explicit that parametric actuation changes both damping and resonance frequency (Vovrosh et al., 2016).

Several variants modify the feedback waveform while retaining the same principle. In digital parametric feedback, the trap stiffness is switched between two levels through

$F_{\rm th}(t)$9

leading to an extra damping

$$\langle F_{\rm th}(t)F_{\rm th}(t')\rangle = 2m\Gamma_0 k_B T_0 \delta(t-t')$$0

and a temperature

$$\langle F_{\rm th}(t)F_{\rm th}(t')\rangle = 2m\Gamma_0 k_B T_0 \delta(t-t')$$1

For the square-wave case, optimal cooling corresponds to $$\langle F_{\rm th}(t)F_{\rm th}(t')\rangle = 2m\Gamma_0 k_B T_0 \delta(t-t')$$2, unlike the sinusoidal analog, where the optimum phase is tied to the sign of $$\langle F_{\rm th}(t)F_{\rm th}(t')\rangle = 2m\Gamma_0 k_B T_0 \delta(t-t')$$3 (Zheng et al., 2019). This difference is a frequent source of confusion when comparing implementations.

3. Experimental realizations across platforms

The earliest levitated-nanoparticle demonstrations emphasized that a laser-trapped nanoparticle is entirely isolated from the thermal bath and lacks a clamping mechanism, enabling robust decoupling from internal vibrations and cooling in all degrees of freedom by means of a single laser beam (Gieseler et al., 2012). In one implementation, a single 1064 nm laser of approximately 100 mW was focused by an NA $$\langle F_{\rm th}(t)F_{\rm th}(t')\rangle = 2m\Gamma_0 k_B T_0 \delta(t-t')$$4 objective into a vacuum chamber to trap a fused-silica nanosphere of radius $$\langle F_{\rm th}(t)F_{\rm th}(t')\rangle = 2m\Gamma_0 k_B T_0 \delta(t-t')$$5 nm and mass $$\langle F_{\rm th}(t)F_{\rm th}(t')\rangle = 2m\Gamma_0 k_B T_0 \delta(t-t')$$6 kg. Three balanced photodetectors monitored forward-scattered light in $$\langle F_{\rm th}(t)F_{\rm th}(t')\rangle = 2m\Gamma_0 k_B T_0 \delta(t-t')$$7, $$\langle F_{\rm th}(t)F_{\rm th}(t')\rangle = 2m\Gamma_0 k_B T_0 \delta(t-t')$$8, and $$\langle F_{\rm th}(t)F_{\rm th}(t')\rangle = 2m\Gamma_0 k_B T_0 \delta(t-t')$$9, with a phase-sensitive interferometric noise floor of approximately $k(t)=k_0[1+\epsilon(t)]$0 pm/$k(t)=k_0[1+\epsilon(t)]$1 (Gieseler et al., 2012). A related parabolic-mirror geometry used a single-mode 1550 nm fibre laser, a single photodiode, and a single beam for trapping, position detection, and cooling of all three dimensions, with backscattered self-homodyne detection reaching a position sensitivity of approximately $k(t)=k_0[1+\epsilon(t)]$2 fm/$k(t)=k_0[1+\epsilon(t)]$3 (Vovrosh et al., 2016).

Feedback-loop architectures differ substantially. The standard analog loop in nanoparticle trapping differentiates each detector signal, multiplies $k(t)=k_0[1+\epsilon(t)]$4 to generate a component at $k(t)=k_0[1+\epsilon(t)]$5, applies a phase shift, sums the three axes’ $k(t)=k_0[1+\epsilon(t)]$6 signals, and drives a Pockels cell or AOM so that the trapping-laser power becomes $k(t)=k_0[1+\epsilon(t)]$7 (Gieseler et al., 2012). Digital parametric feedback replaces continuous phase-shifters by FPGA processing: the position signals are digitized, Kalman-filtered, delay-compensated, converted into synchronous square waves, and combined through a “majority voting” logic for tri-axial cooling (Zheng et al., 2019). In cavityless levitated optomechanics, the control signal is derived from an LQG formulation and applied by an AOM through $k(t)=k_0[1+\epsilon(t)]$8, while the state estimator explicitly incorporates the time-dependent stiffness $k(t)=k_0[1+\epsilon(t)]$9 to avoid lock loss at large modulation depth (Ferialdi et al., 2019).

Outside center-of-mass levitation, the same principle has been adapted to rotational and atomic degrees of freedom. In optically levitated libration, a backward-scattering heterodyne scheme with a 1550 nm trapping beam and an LO offset by $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2[1+\epsilon(t)]x = F_{\rm th}(t).$$0 Hz provided linear signals for all three libration modes, and the PLL phase was shifted by $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2[1+\epsilon(t)]x = F_{\rm th}(t).$$1 to approximate $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2[1+\epsilon(t)]x = F_{\rm th}(t).$$2 before driving an EOM at $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2[1+\epsilon(t)]x = F_{\rm th}(t).$$3 (Gao et al., 2024). For a single atom in an optical cavity, parametric feedback combined trap-depth modulation at $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2[1+\epsilon(t)]x = F_{\rm th}(t).$$4 with fast repetitive measurements of the atomic position derived from cavity transmission and processed on an FPGA through an IQ demodulator (Sames et al., 2018). In pendulum-based gravity experiments, the restoring stiffness was modulated mechanically via vertical suspension-point motion $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2[1+\epsilon(t)]x = F_{\rm th}(t).$$5, with a digital phase tracker maintaining the actuator at $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2[1+\epsilon(t)]x = F_{\rm th}(t).$$6 and near-optimal phase (Hartwig et al., 2020).

4. Demonstrated performance and limiting mechanisms

The 2012 nanoparticle experiment reported bare trap frequencies $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2[1+\epsilon(t)]x = F_{\rm th}(t).$$7 kHz, $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2[1+\epsilon(t)]x = F_{\rm th}(t).$$8 kHz, and $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2[1+\epsilon(t)]x = F_{\rm th}(t).$$9 kHz, with $\epsilon(t)=\epsilon_0\cos(2\Omega_0 t+\phi)$0 in the free-molecular regime and $\epsilon(t)=\epsilon_0\cos(2\Omega_0 t+\phi)$1 mHz at $\epsilon(t)=\epsilon_0\cos(2\Omega_0 t+\phi)$2 mBar, corresponding to $\epsilon(t)=\epsilon_0\cos(2\Omega_0 t+\phi)$3 (Gieseler et al., 2012). With feedback engaged, all three axes cooled; in the lowest-noise run, the reported values were $\epsilon(t)=\epsilon_0\cos(2\Omega_0 t+\phi)$4 mK at $\epsilon(t)=\epsilon_0\cos(2\Omega_0 t+\phi)$5 mBar, $\epsilon(t)=\epsilon_0\cos(2\Omega_0 t+\phi)$6 K at $\epsilon(t)=\epsilon_0\cos(2\Omega_0 t+\phi)$7 mBar, and $\epsilon(t)=\epsilon_0\cos(2\Omega_0 t+\phi)$8 K at the same pressure, with overall compression factors $\epsilon(t)=\epsilon_0\cos(2\Omega_0 t+\phi)$9 up to $2\Omega_0$0 (Gieseler et al., 2012). In the parabolic-mirror trap, room-temperature motion was cooled to as little as $2\Omega_0$1 mK in all three axes, with $2\Omega_0$2 mHz at $2\Omega_0$3 mbar and estimated $2\Omega_0$4 (Vovrosh et al., 2016). Digital square-wave cooling achieved center-of-mass temperatures of approximately $2\Omega_0$5 mK at $2\Omega_0$6 mbar with modulation depth up to $2\Omega_0$7 on the $2\Omega_0$8-axis (Zheng et al., 2019).

Librational cooling has reached lower occupancies in selected modes. Using backward-scattering detection, all three libration degrees of freedom were cooled below $2\Omega_0$9 mK, with one mode at $\phi$0 mK corresponding to $\phi$1, and the measurement efficiency for the $\phi$2 mode was reported as approximately $\phi$3 (Gao et al., 2024). In single-atom cavity QED, parametric cooling of a radial mode at $\phi$4 kHz increased the average storage time by a factor of $\phi$5 to more than $\phi$6 s, while a $\phi$7 kHz axial mode exhibited a $\phi$8 extension of trap lifetime within microseconds of cooling (Sames et al., 2018). For pendulum-based gravity experiments, a proof-of-principle demonstration achieved a damping factor of $\phi$9, reducing the effective mode temperature from $2\Omega_0$0 K to approximately $2\Omega_0$1 K in the seismic-noise-dominated regime (Hartwig et al., 2020).

The limiting mechanisms are platform-dependent but structurally similar. Reported limits include reheating by gas collisions, measurement noise floor, recoil heating, finite PLL bandwidth, phase jitter, slow trap-frequency drift, cross-coupling to other modes, and actuator or electronics noise (Gieseler et al., 2012); (Penny et al., 2021); (Gao et al., 2024). In the original nanoparticle experiment, the measurement noise estimate $2\Omega_0$2 pm in $2\Omega_0$3 Hz implied a theoretical $2\Omega_0$4 $2\Omega_0$5K in the absence of back-action, while photon-recoil heating was estimated to give only one event per $2\Omega_0$6 oscillations for $2\Omega_0$7 nm and $2\Omega_0$8 nm (Gieseler et al., 2012). In the direct comparison study, two bounds were emphasized: a stability-bandwidth limit $2\Omega_0$9 and a loop-SNR limit $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$00, with the practical minimum approximated by $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$01 (Penny et al., 2021).

5. Optimal, adaptive, and quantum formulations

A major development beyond fixed-phase feedback is the explicit use of optimal control and adaptive phase updates. In cavityless levitated optomechanics, parametric cooling has been cast as an LQG-type problem with performance index

$$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$02

leading through Pontryagin’s Minimum Principle and the Riccati equation to the feedback law

$$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$03

or $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$04 in the infinite-horizon limit (Ferialdi et al., 2019). The essential implementation point is that the Kalman filter or PLL must incorporate the known modulation $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$05 through the time-dependent matrix $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$06; this allows stable tracking at modulation depths up to $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$07, whereas the conventional double-phase scheme cannot be pushed much beyond $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$08 without losing lock (Ferialdi et al., 2019). The reported consequence is up to $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$09 faster cooling and $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$10 lower $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$11 than standard double-phase schemes (Ferialdi et al., 2019).

A distinct adaptive strategy updates the modulation phase directly from measured quadratures. In phase-adaptive parametric cooling, the phase at the $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$12-th update is set to

$$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$13

with repeated stroboscopic updates yielding purely exponential energy decay,

$$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$14

The classical and quantum steady-state occupancies are then

$$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$15

for $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$16 (Ghosh et al., 2022). This result explicitly distinguishes adaptive phase control from plain fixed-phase parametric cooling, which is described there as producing a non-exponential, slow energy decay (Ghosh et al., 2022).

Quantum analyses impose stronger constraints. A quantum calculation in the shot-noise-dominant regime derived a stochastic master equation with measurement efficiency $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$17, shot-noise heating per oscillation $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$18, and parametric-feedback Hamiltonian $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$19 (Zhong et al., 2017). In that treatment, the minimum occupation under parametric feedback scales approximately as

$$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$20

and ground-state cooling requires substantially higher $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$21 than force feedback (Zhong et al., 2017). By contrast, an optimal quantum protocol based on heterodyne measurement, resonant parametric modulation, and conditional choice of both phase and duration derived the phase relation

$$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$22

the optimal squeezing duration

$$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$23

and a steady-state fixed point $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$24, even for an isolated oscillator $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$25 (Manikandan et al., 2022). This suggests that measurement-conditioned parametric modulation can, in principle, cool below one quantum without sideband or linear-feedback cooling (Manikandan et al., 2022).

6. Comparisons, misconceptions, and scope

Parametric feedback cooling is often discussed alongside velocity damping, cold damping, and cavity sideband cooling, but the methods are not interchangeable. In a direct comparison performed on the same levitated particle and with the same detection system, velocity damping cooled the oscillator to $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$26 mK, almost an order of magnitude below the best parametric-feedback result of approximately $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$27 mK, and was reported to be more resilient to imperfect experimental conditions (Penny et al., 2021). The reasons given include lower effective back-action, automatic tracking of frequency drift, and reduced sensitivity to cross-couplings and phase-estimation errors (Penny et al., 2021). Likewise, in the shot-noise-dominant quantum calculation, force feedback reached lower occupation than parametric feedback at fixed measurement efficiency (Zhong et al., 2017). A common misconception is therefore that parametric feedback is universally the lowest-temperature route; the comparative literature does not support that claim.

Another misconception is that parametric feedback acts identically on all mechanical coordinates. In rigid-body nanodumbbells trapped in a linearly polarized laser beam, standard parametric feedback can extract energy from two of the five rotational degrees of freedom, but the dynamics after feedback are characterized by a normal mode describing precession about the laser polarization axis together with spin about the nanoparticle’s symmetry axis (Seberson et al., 2018). Full cooling of the librational coordinates requires an asymmetry in the librational frequencies and feedback modulation containing both rotational frequencies (Seberson et al., 2018). This is not a minor technicality; it shows that the modal structure of the underlying Hamiltonian can obstruct naïve extensions of center-of-mass protocols to rotational motion.

At the same time, parametric feedback remains attractive in geometries with limited optical access or restricted actuator choices. In the single-atom cavity experiment, only one optical mode was used for both measurement and actuation, and the method remained effective for a $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$28 kHz oscillation mode within microseconds (Sames et al., 2018). In levitated optomechanics, the single-beam implementations and digital FPGA realizations show that all-optical or minimally invasive control of several motional degrees of freedom is feasible (Vovrosh et al., 2016); (Zheng et al., 2019). Applications identified across the literature include ultrasensitive force sensing with $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$29 N/$$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$30, tests of quantum mechanics with mesoscopic objects, searches for nonstandard forces, quantum-state engineering, damping of pendulum modes in gravitational-wave detectors, and preparation of rotational states relevant to matter-wave interferometry (Gieseler et al., 2012); (Hartwig et al., 2020).

The long-term outlook is framed by the competition among damping, measurement imprecision, recoil or technical heating, and phase-tracking fidelity. For laser-trapped nanoparticles with $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$31 kHz, the condition $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$32 implies $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$33 pK, and extrapolation of $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$34 to $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$35 mBar together with higher feedback gain was identified as a route toward $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$36 pK (Gieseler et al., 2012). More recent analyses instead emphasize that reaching the near-ground-state regime requires not only high $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$37 and ultralow gas damping, but also measurement and estimator designs that remain accurate under strong modulation (Ferialdi et al., 2019); (Ghosh et al., 2022); (Manikandan et al., 2022). A plausible implication is that the future of parametric feedback cooling lies less in the elementary $$m\ddot x + m\Gamma_0 \dot x + m\Omega_0^2 x = F_{\rm th}(t),$$38 loop by itself than in hybrid architectures that combine parametric modulation with state estimation, adaptive phase control, or complementary cooling channels.