Multi-Frequency Cooling Algorithms

• Multi-frequency cooling algorithms are techniques that leverage multiple spectral components to drive efficient energy extraction in quantum systems.
• They encompass diverse methods such as heat-bath algorithmic cooling, measurement-based cooling, engineered dissipation, and optical frequency comb strategies.
• These protocols rely on spectral matching and optimization of control parameters to balance speed, fidelity, and asymptotic cooling performance.

Multi-frequency cooling algorithm denotes a class of cooling procedures in which entropy extraction or energy relaxation is organized by multiple spectral components, multiple resonant channels, or a nontrivial spectrum of decay modes. In the literature surveyed here, the term spans several distinct constructions: fixed-map heat-bath algorithmic cooling whose relaxation is resolved into discrete eigenmodes, measurement-based cooling timed by dominant or collective Rabi frequencies, engineered dissipation using ancilla splittings sampled from a frequency window, quantum digital cooling via sweeps over ancilla energies, and optical schemes based on frequency combs or multiple red-detuned MOT components (Raeisi et al., 2019, Yan et al., 2022, Schlömer et al., 1 May 2026, Polla et al., 2019, Ahmad et al., 2012, Hopton et al., 25 Apr 2026). This suggests that the expression functions less as the name of a single canonical protocol than as a spectral-design paradigm for cooling dynamics.

1. Scope of the concept

The surveyed works use “multi-frequency” in several technically distinct senses. In heat-bath algorithmic cooling, the relevant “frequencies” are the eigenmodes of a time-homogeneous Markov chain on populations; in measurement-based resonator cooling they are dominant JC Rabi frequencies such as $\Omega_d$, $\Omega_{d+1}$, and the collective thermal Rabi frequency $\Omega_{\rm th}$; in engineered dissipation for interacting fermions they are the ancilla splittings $\omega_\ell$ or bath frequencies $\Delta_r$; and in laser-cooling realizations they are literal optical frequencies, comb teeth, or multiple red-detuned spectral components (Raeisi et al., 2019, Yan et al., 2022, Schlömer et al., 1 May 2026, Jayich et al., 2016, Hopton et al., 25 Apr 2026).

Domain	Multi-frequency object	Representative mechanism
HBAC / fixed-map AC	Eigenvalues $\lambda_k$ of a population-transfer matrix	Repeated fixed compression map damps all $|\lambda_k|<1$ and leaves the OAS (Raeisi et al., 2019)
Measurement-based cooling	$\Omega_d,\Omega_{d+1},\Omega_{\rm th}$	Measurement interval chosen from dominant or collective Rabi scales (Yan et al., 2022)
Multimode ancilla cooling	$\Omega_{\rm th}=\sqrt{\sum_k g_k^2\bar n_k}$	One ancilla cools several resonators simultaneously (Yan et al., 2022)
Dissipative many-body cooling	Random $\omega_\ell$ or averaged $\Omega_{d+1}$0	Broadband or multi-frequency ancilla-assisted energy extraction (Schlömer et al., 1 May 2026, Molpeceres et al., 31 Mar 2025)
Frequency-comb and MOT cooling	Comb teeth or multiple red-detuned optical components	Different transitions or velocity classes are addressed in parallel (Ahmad et al., 2012, Buhin et al., 2020, Hopton et al., 25 Apr 2026)

A recurrent structural theme is spectral matching. Cooling is efficient when the control spectrum overlaps energy-lowering transitions while suppressing reverse processes. The main design variables are therefore frequency placement, linewidth or effective broadening, cycle time, reset strategy, and the locality or symmetry properties of the coupling operators.

2. Algorithmic and digital cooling on spin and qubit registers

In HBAC, the system is partitioned into computation qubits and reset qubits. Standard protocols iterate an entropy-compression unitary followed by rethermalization of the reset qubit to a state with polarization $\Omega_{d+1}$1. The Partner Pairing Algorithm sorts diagonal populations at each iteration, so the required compression depends on the instantaneous state. The two-sort algorithmic cooling protocol (TSAC) replaces this state-dependent compression with a single fixed unitary $\Omega_{d+1}$2 that swaps every pair of neighboring populations except the first and last, followed by the reset map

$\Omega_{d+1}$3

The induced population dynamics is a stochastic map $\Omega_{d+1}$4 with a unique fixed point equal to the optimal asymptotic state

$\Omega_{d+1}$5

and eigenvalues

$\Omega_{d+1}$6

plus $\Omega_{d+1}$7. The protocol therefore reaches the HBAC limit asymptotically, but the mixing time scales as $\Omega_{d+1}$8, while the per-iteration gate complexity is $\Omega_{d+1}$9 and the total gate complexity near the OAS is $\Omega_{\rm th}$0. The paper also proves a monotonicity statement: if the first computation qubit is below the HBAC limit, each TSAC iteration increases its polarization (Raeisi et al., 2019). In this setting, “multi-frequency” refers to the spectrum of the Markov chain rather than to explicit physical drive tones.

A separate line of work formulates cooling digitally with a resettable ancilla qubit of splitting $\Omega_{\rm th}$1, coupled through

$\Omega_{\rm th}$2

For a one-qubit toy model, resonant cooling gives $\Omega_{\rm th}$3, while reheating depends on detuning and coupling strength. Two scalable protocols follow from this analysis. LogSweep uses weak coupling and logarithmically spaced ancilla energies $\Omega_{\rm th}$4 to sweep an interval of transition energies, thereby implementing an explicit multi-frequency energy scan. BangBang uses strong coupling and locally chosen frequencies $\Omega_{\rm th}$5 to maximize rapid de-excitation with shallow depth. LogSweep was tested on the 1D transverse-field Ising model and yields approximate ground-state preparation with an error that can be made polynomially small in the computation time; BangBang is faster but does not promise long-time convergence (Polla et al., 2019).

Earlier algorithmic-cooling work on spins provides a multiscale antecedent. Semi-optimal practicable AC introduces a fixed number $\Omega_{\rm th}$6 of compression cycles per recursive level, with mPAC obeying

$\Omega_{\rm th}$7

in the small-bias regime, and the level-dependent $\Omega_{\rm th}$8PAC generalization giving

$\Omega_{\rm th}$9

This suggests a multi-frequency interpretation in the timescale hierarchy: different recursion levels are “driven” with different cycle counts $\omega_\ell$0, interpolating between PAC and exhaustive all-bonacci-type cooling (Elias et al., 2011).

3. Measurement-based and collective-mode formulations

Measurement-based cooling of a resonator coupled to a qubit via the JC Hamiltonian

$\omega_\ell$1

makes the multi-frequency structure explicit. The relevant mode-resolved Rabi frequencies are

$\omega_\ell$2

For conditional measurement (CM), the Fock population $\omega_\ell$3 is multiplied by

$\omega_\ell$4

For unconditional measurement (UM), the update rule is

$\omega_\ell$5

The central analytical result is the frequency-dependent optimal UM interval

$\omega_\ell$6

where $\omega_\ell$7 and $\omega_\ell$8, together with the CM rule

$\omega_\ell$9

A PPO actor-critic agent then chooses between UM and CM, optimizing the cooperative cooling performance

$\Delta_r$0

For the representative parameter set $\Delta_r$1, $\Delta_r$2, $\Delta_r$3, and $\Delta_r$4, the RL-optimized sequence reduces the average population by four orders of magnitude in 16 rounds, achieves $\Delta_r$5, and keeps $\Delta_r$6; pure CM gives deeper cooling but $\Delta_r$7, whereas pure UM preserves $\Delta_r$8 but cools weakly (Yan et al., 2022).

The same logic extends to simultaneous cooling of several modes by repeated projection of a single ancilla onto its ground state. For two resonators coupled to a $\Delta_r$9-type qutrit, the nonunitary propagator is diagonal in the joint Fock basis, and under the resonant condition the relevant double-mode Rabi frequency is

$\lambda_k$0

For $\lambda_k$1 modes, the collective thermal Rabi frequency becomes

$\lambda_k$2

The protocol scales to more resonators and is robust to moderate frequency mismatch. In the reported simulations, about 20 optimized measurements reduce the average populations by six orders of magnitude, and a 5-mode example reaches $\lambda_k$3 for all modes after 20 measurements (Yan et al., 2022). Here the “multi-frequency” label refers to simultaneous cooling of non-degenerate modes through a collective timing rule rather than to multiple external drive frequencies.

4. Engineered dissipation for many-body fermions

For strongly interacting fermions, the most literal multi-frequency implementation is randomized dissipative cooling with ancilla qubits whose splittings are sampled from a window $\lambda_k$4 at every cycle. The joint Hamiltonian is

$\lambda_k$5

with $\lambda_k$6 and local, symmetry-preserving operators $\lambda_k$7. A cycle consists of ancilla initialization, random frequency sampling, unitary evolution, and ancilla reset. In a Fermi’s-golden-rule picture, a transition $\lambda_k$8 is resonantly enhanced when $\lambda_k$9, with probability

$|\lambda_k|<1$0

Because $|\lambda_k|<1$1 is redrawn from a continuous distribution, the effective bath is broadband and model-agnostic. Benchmarks on the spinless $|\lambda_k|<1$2–$|\lambda_k|<1$3–$|\lambda_k|<1$4 chain, the 1D $|\lambda_k|<1$5–$|\lambda_k|<1$6 chain, a mixed-dimensional ladder, and 2D spinless fermions show monotonic energy relaxation, concentration of spectral weight at low energies, and stabilization of correlated order, including charge order, pairing, and superconducting correlations (Schlömer et al., 1 May 2026).

A complementary solvable case study analyzes a quadratic fermion chain cooled by repeated coupling to a reset bath with one or several frequencies. In the weak-coupling limit, the single-frequency protocol yields mode-resolved jump operators with coefficients

$|\lambda_k|<1$7

so heating and cooling compete. Randomizing the cycle time turns the discrete interference structure into Lorentzian rates, and using $|\lambda_k|<1$8 bath frequencies gives averaged rates

$|\lambda_k|<1$9

Multi-frequency and randomized cycles significantly enhance cooling and circumvent the accidental heating resonances of the fixed-time, single-frequency protocol. The paper further shows that with $\Omega_d,\Omega_{d+1},\Omega_{\rm th}$0, $\Omega_d,\Omega_{d+1},\Omega_{\rm th}$1, and $\Omega_d,\Omega_{d+1},\Omega_{\rm th}$2, constant-fidelity near-ground-state cooling can be achieved in total time $\Omega_d,\Omega_{d+1},\Omega_{\rm th}$3. In the depolarizing-noise model, the rates acquire an additive term,

$\Omega_d,\Omega_{d+1},\Omega_{\rm th}$4

and resonant cooling remains feasible when $\Omega_d,\Omega_{d+1},\Omega_{\rm th}$5 (Molpeceres et al., 31 Mar 2025). Relative to dissipative state preparation in the same model, cooling generally reaches lower energies and is more resilient to noise (Molpeceres et al., 31 Mar 2025).

5. Frequency-comb, dual-species, and MOT realizations

In laser cooling, the term “multi-frequency” is literal. A coherent pulse train produces a frequency comb with teeth

$\Omega_d,\Omega_{d+1},\Omega_{\rm th}$6

or, in standard notation,

$\Omega_d,\Omega_{d+1},\Omega_{\rm th}$7

The “see-saw” protocol for a three-level $\Omega_d,\Omega_{d+1},\Omega_{\rm th}$8 system switches the carrier-envelope phase offset between values $\Omega_d,\Omega_{d+1},\Omega_{\rm th}$9 and $\Omega_{\rm th}=\sqrt{\sum_k g_k^2\bar n_k}$0, thereby moving comb teeth in and out of resonance with $\Omega_{\rm th}=\sqrt{\sum_k g_k^2\bar n_k}$1 and $\Omega_{\rm th}=\sqrt{\sum_k g_k^2\bar n_k}$2. The switching period $\Omega_{\rm th}=\sqrt{\sum_k g_k^2\bar n_k}$3 is optimized numerically; for the reported examples, $\Omega_{\rm th}=\sqrt{\sum_k g_k^2\bar n_k}$4, meaning a switch every pulse. Time evolution of the simulated velocity distribution shows both deceleration and narrowing after 25,000 and 100,000 pulses, illustrating true cooling rather than mere slowing (Ahmad et al., 2012).

Direct comb cooling of atoms on a two-photon transition uses a different frequency-space object: the effective two-photon comb has teeth

$\Omega_{\rm th}=\sqrt{\sum_k g_k^2\bar n_k}$5

with a two-photon Rabi frequency

$\Omega_{\rm th}=\sqrt{\sum_k g_k^2\bar n_k}$6

For the Rb $\Omega_{\rm th}=\sqrt{\sum_k g_k^2\bar n_k}$7 demonstration, the comb-based two-photon Doppler limit is

$\Omega_{\rm th}=\sqrt{\sum_k g_k^2\bar n_k}$8

with an ideal prediction of $\Omega_{\rm th}=\sqrt{\sum_k g_k^2\bar n_k}$9, a broadened-linewidth estimate of $\omega_\ell$0, and an experimental value of $\omega_\ell$1 (Jayich et al., 2016). The same comb logic scales naturally to simultaneous multi-transition or multi-species addressing. In the dual-species $\omega_\ell$2 experiment, two distinct comb modes are simultaneously red detuned from the two cycling transitions; with $\omega_\ell$3, the relevant modes differ by 14 teeth, reflecting the 1126 MHz isotopic splitting, and the two isotopes cool simultaneously to the Doppler temperature without degrading each other’s cooling characteristics (Buhin et al., 2020).

A vapor-cell MOT implements yet another multi-frequency construction by adding several closely spaced red-detuned components to the cooling light of a standard $\omega_\ell$4Rb MOT. The additional components are typically placed around $\omega_\ell$5 to $\omega_\ell$6 MHz, i.e. roughly one linewidth apart, and formed into an annular beam that overlaps the outer trapping region but avoids the dense core. The mechanism is straightforward: different detunings decelerate different velocity classes, extending the force–velocity curve and increasing the capture velocity $\omega_\ell$7, while the loading rate from a thermal vapor scales as $\omega_\ell$8. Experimentally, this yields up to a factor-of-4 increase in loading rate and a doubling of steady-state atom number; the reported maxima are $\omega_\ell$9 and $\Omega_{d+1}$00 trapped atoms, obtained without additional slowing stages (Hopton et al., 25 Apr 2026).

6. Common design principles, trade-offs, and limitations

Across these implementations, efficient cooling requires a unique low-energy fixed point together with suppression of reverse processes. In TSAC this condition is encoded in a stochastic matrix with one eigenvalue at 1 and all others in $\Omega_{d+1}$01 (Raeisi et al., 2019). In measurement-based cooling it appears as interval choices such as $\Omega_{d+1}$02 or $\Omega_{d+1}$03, which maximize de-excitation while protecting the vacuum (Yan et al., 2022, Yan et al., 2022). In digital and dissipative fermionic cooling it is achieved by sweeping or randomizing ancilla energies so that many energy-lowering transitions are resonantly addressed while reheating remains off-resonant (Polla et al., 2019, Schlömer et al., 1 May 2026, Molpeceres et al., 31 Mar 2025). In optical implementations it is realized by distributing spectral weight over comb teeth or red-detuned MOT components so that different internal transitions or velocity classes are cooled in parallel (Ahmad et al., 2012, Hopton et al., 25 Apr 2026).

The main trade-offs are equally consistent. State-independent maps simplify control but may converge slowly: TSAC has a fixed circuit and asymptotic optimality, yet its saturation can require exponentially many iterations (Raeisi et al., 2019). Narrow-band, weak-coupling strategies suppress reheating but demand many cycles or deep circuits, as in LogSweep and in the weak-coupling fermion analyses (Polla et al., 2019, Molpeceres et al., 31 Mar 2025). Strong-coupling or hybrid protocols cool faster but sacrifice asymptotic guarantees or success probability, as seen in BangBang and in CM-dominated measurement sequences (Polla et al., 2019, Yan et al., 2022). In optical systems, adding frequencies is beneficial only if blue-detuned components, deleterious interference, or collision-enhanced losses are controlled; the MOT study attributes earlier limitations largely to symmetric broadband spectra and to the overlap of highly red-detuned light with the dense trapped core (Hopton et al., 25 Apr 2026).

Several common misconceptions are corrected by the literature. One is that “multi-frequency” must refer to multiple externally applied tones. In fact, TSAC uses the phrase only indirectly through a spectral decomposition of the cooling channel, and simultaneous ancilla-based resonator cooling encodes the relevant frequencies in collective Rabi scales rather than in distinct drives (Raeisi et al., 2019, Yan et al., 2022). Another is that more frequencies always improve performance. The solvable fermion case shows that gains from increasing $\Omega_{d+1}$04 eventually saturate or reverse when the weak-coupling condition breaks down, and the MOT analysis shows that spectral asymmetry and spatial separation are as important as the raw number of components (Molpeceres et al., 31 Mar 2025, Hopton et al., 25 Apr 2026).

A plausible synthesis is that multi-frequency cooling algorithms are best understood as spectral-engineering schemes. Their distinguishing feature is not a particular hardware platform or update rule, but the deliberate shaping of a control spectrum—whether an eigenvalue spectrum, a set of ancilla splittings, a ladder of measurement intervals, or an optical comb—to make all relevant non-equilibrium modes relax while preserving the desired cold manifold.