Skip to content

Quantum Heat Engines and Refrigerators

A quantum heat engine is a thermodynamic machine whose working medium is described quantum mechanically. The working medium may be a two-level system, harmonic oscillator, spin ensemble, atom-cavity system, quantum dot, superconducting circuit, or another small system coupled to reservoirs and controls.

The word “quantum” does not suspend the ordinary thermodynamic questions. One still asks:

  • which reservoirs provide heat;
  • which degree of freedom stores or receives work;
  • whether the process is cyclic or steady;
  • whether the heat currents are defined by a thermodynamically consistent model;
  • whether coherence, measurement, or finite-time driving changes power or irreversibility.

Quantum heat engines are most useful as precise models of microscopic energy conversion. They are not loopholes around Carnot’s bound. When the reservoirs are genuine thermal reservoirs and the accounting is complete, the usual second-law bounds reappear.

This chapter uses QαQ_\alpha for heat transferred into the working medium from reservoir α\alpha. For a cyclic engine between a hot bath hh and cold bath cc,

ΔEcycle=0.\Delta E_{\mathrm{cycle}} = 0 .

Let

Qh>0,Qc<0,Q_h\gt0, \qquad Q_c\lt0,

where QhQ_h is heat absorbed from the hot bath and QcQ_c is heat dumped into the cold bath. Positive extracted work is

Wext=Qh+Qc.W_{\mathrm{ext}} = Q_h+Q_c .

The efficiency is

η=WextQh=1+QcQh.\eta = \frac{W_{\mathrm{ext}}}{Q_h} = 1+\frac{Q_c}{Q_h}.

The same process described using work supplied to the system would use Wsup=−WextW_{\mathrm{sup}}=-W_{\mathrm{ext}}. Many sign mistakes in heat-engine papers are just convention changes that were not stated clearly.

For an ideal cycle coupled only to two thermal reservoirs at temperatures ThT_h and TcT_c, with Th>TcT_h\gt T_c, the working medium returns to its initial state. The total entropy production is therefore

Σ=−QhTh−QcTc≥0.\Sigma = - \frac{Q_h}{T_h} - \frac{Q_c}{T_c} \ge0 .

Since Qh>0Q_h\gt0, this implies

Qc≤−TcThQh.Q_c \le - \frac{T_c}{T_h}Q_h .

Substituting into the efficiency gives

η≤1−TcTh≡ηC.\eta \le 1-\frac{T_c}{T_h} \equiv \eta_{\mathrm C}.

This derivation is deliberately compact. It shows what must be true for the bound to apply: the reservoirs must be thermal, all entropy changes outside the working medium must be included, and the cycle must really close. If a squeezed reservoir, measurement record, chemical potential, or battery is present, it must be counted as an additional resource rather than hidden inside the words “hot bath.”

A refrigerator uses supplied work to move heat from a cold reservoir to a hot reservoir. It is often clearer to use positive quantities

Qcin>0,Win>0,Q_c^{\mathrm{in}}\gt0, \qquad W_{\mathrm{in}}\gt0,

where QcinQ_c^{\mathrm{in}} is heat extracted from the cold bath and WinW_{\mathrm{in}} is work supplied by the external source. The heat delivered to the hot bath is

Qhout=Qcin+Win.Q_h^{\mathrm{out}} = Q_c^{\mathrm{in}}+W_{\mathrm{in}} .

The coefficient of performance for refrigeration is

COPR=QcinWin.\mathrm{COP}_{\mathrm R} = \frac{Q_c^{\mathrm{in}}} {W_{\mathrm{in}}}.

For two ideal thermal reservoirs,

COPR≤TcTh−Tc.\mathrm{COP}_{\mathrm R} \le \frac{T_c} {T_h-T_c}.

The same device can be called a heat pump when the useful output is heat delivered to the hot reservoir. Then

COPHP=QhoutWin≤ThTh−Tc.\mathrm{COP}_{\mathrm{HP}} = \frac{Q_h^{\mathrm{out}}} {W_{\mathrm{in}}} \le \frac{T_h} {T_h-T_c}.

Again, the bound assumes ideal thermal reservoirs and complete resource accounting.

Quantum heat engines are commonly grouped into two architectural classes.

Stroke engines separate the tasks into time intervals:

  • thermalization with a hot reservoir;
  • isolated unitary driving;
  • thermalization with a cold reservoir;
  • another isolated unitary stroke.

The quantum Otto engine is the standard example. Its clean separation makes it useful pedagogically and experimentally, especially in trapped ions, spins, superconducting circuits, and oscillator platforms.

Autonomous engines operate continuously. Different transitions or subsystems are coupled to different reservoirs, and work is stored in a mode, load, field, or battery. The three-level maser is the classic autonomous model.

The two architectures can share the same thermodynamic efficiency in ideal limits, but they distribute irreversibility differently. Stroke engines emphasize finite-time driving and quantum friction. Autonomous engines emphasize steady-state currents, spectral filtering, and load modeling.

In weak-coupling Markovian models, the working medium may obey

dρdt=−iℏ[H(t),ρ]+∑αLαt(ρ),\frac{d\rho}{dt} = - \frac{i}{\hbar}[H(t),\rho] + \sum_\alpha \mathcal L_\alpha^t(\rho),

where each dissipator Lαt\mathcal L_\alpha^t represents a reservoir. A heat current into the system is often defined as

Q˙α=Tr⁡[H(t)Lαt(ρ)].\dot Q_\alpha = \operatorname{Tr} \left[ H(t)\mathcal L_\alpha^t(\rho) \right].

The power supplied by explicit driving is

W˙sup=Tr⁡[ρdHdt].\dot W_{\mathrm{sup}} = \operatorname{Tr} \left[ \rho\frac{dH}{dt} \right].

The instantaneous first-law balance is then

ddtTr⁡[ρH]=W˙sup+∑αQ˙α.\frac{d}{dt} \operatorname{Tr}[\rho H] = \dot W_{\mathrm{sup}} + \sum_\alpha \dot Q_\alpha .

For a periodically driven steady cycle, the average internal-energy change over one period vanishes, so

W˙‾ext=∑αQ˙‾α,W˙‾ext=−W˙‾sup.\overline{\dot W}_{\mathrm{ext}} = \sum_\alpha \overline{\dot Q}_\alpha, \qquad \overline{\dot W}_{\mathrm{ext}} = - \overline{\dot W}_{\mathrm{sup}} .

This bookkeeping becomes thermodynamics only if the dissipators are physically tied to reservoirs. For thermal reservoirs, the rates should satisfy detailed balance with respect to the Hamiltonian used in the heat current. See Thermal Master Equations and Detailed Balance.

Let the working medium be a two-level system with Hamiltonian

H(ω)=ℏω∣e⟩⟨e∣.H(\omega) = \hbar\omega |e\rangle\langle e|.

The Otto cycle uses two gaps:

ωh>ωc.\omega_h\gt\omega_c .

The four ideal strokes are:

  • hot isochore at gap ωh\omega_h;
  • isolated expansion from ωh\omega_h to ωc\omega_c;
  • cold isochore at gap ωc\omega_c;
  • isolated compression from ωc\omega_c to ωh\omega_h.

Let the excited-state populations after ideal hot and cold thermalization be

ph=11+eβhℏωh,pc=11+eβcℏωc.p_h = \frac{1} {1+e^{\beta_h\hbar\omega_h}}, \qquad p_c = \frac{1} {1+e^{\beta_c\hbar\omega_c}}.

The heat absorbed from the hot bath is

Qh=ℏωh(ph−pc).Q_h = \hbar\omega_h(p_h-p_c).

The heat dumped into the cold bath is

Qc=ℏωc(pc−ph).Q_c = \hbar\omega_c(p_c-p_h).

The net extracted work is

Wext=Qh+Qc=ℏ(ωh−ωc)(ph−pc).W_{\mathrm{ext}} = Q_h+Q_c = \hbar(\omega_h-\omega_c)(p_h-p_c).

The device acts as an engine when

ph>pc,p_h\gt p_c,

equivalently

βhωh<βcωc.\beta_h\omega_h \lt \beta_c\omega_c .

The ideal efficiency is

ηOtto=WextQh=1−ωcωh.\eta_{\mathrm{Otto}} = \frac{W_{\mathrm{ext}}}{Q_h} = 1-\frac{\omega_c}{\omega_h}.

The engine condition ensures that this does not exceed the Carnot efficiency:

1−ωcωh≤1−TcTh.1-\frac{\omega_c}{\omega_h} \le 1-\frac{T_c}{T_h}.

Finite-time driving can create coherence in the instantaneous energy basis. That coherence is often called quantum friction because it raises the energetic cost of closing the cycle unless it is controlled, dephased, or used coherently by another stroke.

For a harmonic oscillator,

H(ω)=ℏω(a†a+12).H(\omega) = \hbar\omega \left( a^\dagger a+\frac{1}{2} \right).

The thermal mean occupation is

nˉ(β,ω)=1eβℏω−1.\bar n(\beta,\omega) = \frac{1} {e^{\beta\hbar\omega}-1}.

In an ideal adiabatic Otto cycle, the occupation number is unchanged during frequency ramps. If

nˉh=nˉ(βh,ωh),nˉc=nˉ(βc,ωc),\bar n_h = \bar n(\beta_h,\omega_h), \qquad \bar n_c = \bar n(\beta_c,\omega_c),

then

Qh=ℏωh(nˉh−nˉc),Qc=ℏωc(nˉc−nˉh).Q_h = \hbar\omega_h(\bar n_h-\bar n_c), \qquad Q_c = \hbar\omega_c(\bar n_c-\bar n_h).

The zero-point contribution cancels over the complete cycle. The extracted work is

Wext=ℏ(ωh−ωc)(nˉh−nˉc),W_{\mathrm{ext}} = \hbar(\omega_h-\omega_c) (\bar n_h-\bar n_c),

and the ideal efficiency is again

ηOtto=1−ωcωh.\eta_{\mathrm{Otto}} = 1-\frac{\omega_c}{\omega_h}.

The oscillator version is common in ion traps, optomechanics, nanomechanical systems, and circuit resonators. It also makes clear why “higher temperature” is not enough by itself; the relevant comparison is the occupation at the two different gaps.

The three-level maser is the classic autonomous quantum heat engine. Let the levels be

E1=0,E2=ℏωc,E3=ℏωh,E_1=0, \qquad E_2=\hbar\omega_c, \qquad E_3=\hbar\omega_h,

with

ωh>ωc,ωw=ωh−ωc.\omega_h\gt\omega_c, \qquad \omega_w=\omega_h-\omega_c .

The hot bath couples the 1↔31\leftrightarrow3 transition, the cold bath couples the 1↔21\leftrightarrow2 transition, and the work output is associated with the 3↔23\leftrightarrow2 transition. Each hot quantum of energy ℏωh\hbar\omega_h can be converted into one work quantum ℏωw\hbar\omega_w plus one cold-bath quantum ℏωc\hbar\omega_c.

The ideal efficiency is therefore

η=ℏωwℏωh=1−ωcωh.\eta = \frac{\hbar\omega_w} {\hbar\omega_h} = 1-\frac{\omega_c}{\omega_h}.

The engine regime requires population inversion on the work transition:

p3>p2.p_3\gt p_2 .

If the hot and cold transitions are locally thermalized, this corresponds to

e−βhℏωh>e−βcℏωc,e^{-\beta_h\hbar\omega_h} \gt e^{-\beta_c\hbar\omega_c},

or

βhωh<βcωc.\beta_h\omega_h \lt \beta_c\omega_c .

The same inequality guarantees that the maser efficiency stays below the Carnot value. In modern language, the output field or load must be modeled explicitly if one wants power, fluctuations, or backaction rather than only an efficiency ratio.

Efficiency is not the only figure of merit. A reversible engine can approach Carnot efficiency only with vanishing entropy production, which usually means vanishing power in ideal quasistatic limits. Practical microscopic engines are judged by several quantities:

  • output power;
  • efficiency or coefficient of performance;
  • stability of the output load;
  • work fluctuations;
  • entropy production;
  • control cost;
  • robustness to parameter uncertainty;
  • heat leaks and parasitic transitions.

For periodically driven engines, average power is

Pout=Wextτcyc,P_{\mathrm{out}} = \frac{W_{\mathrm{ext}}}{\tau_{\mathrm{cyc}}},

where τcyc\tau_{\mathrm{cyc}} is the cycle time. For steady autonomous engines, power is an energy current into the work storage degree of freedom.

Quantum coherence can increase or decrease power depending on the model. It can enable coherent energy exchange with a load, but it can also create nonadiabatic excitations that reduce useful work. A trustworthy claim must specify the Hamiltonian, reservoirs, control protocol, and load.

Quantum heat-engine models differ from classical textbook engines in several concrete ways:

  • the working medium can have discrete spectra and few levels;
  • heat currents may be transition-resolved;
  • coherence can appear during finite-time strokes;
  • measurement backaction can affect feedback engines;
  • work is a process variable, not a one-time observable;
  • strong coupling can make the system-bath boundary thermodynamically relevant;
  • noise and fluctuations are comparable to mean values in small systems.

These features change modeling practice. They do not by themselves violate the second law. When apparent violations appear, the usual diagnosis is an omitted resource: a nonthermal reservoir, measurement record, coherent drive, hidden battery, changing coupling energy, or inconsistent local dissipator.

The same quantum machine can often operate in several modes depending on parameters. For the two-level Otto model, changing frequency ratios and bath temperatures can switch among:

  • engine operation, with Qh>0Q_h\gt0 and Wext>0W_{\mathrm{ext}}\gt0;
  • refrigerator operation, with heat removed from the cold bath while work is supplied;
  • heat pump operation, where useful output is heat delivered to the hot bath;
  • dissipator operation, where supplied work is simply degraded into heat.

For absorption refrigerators, no externally driven work stroke is required. Instead, a third reservoir or transition plays the role of a work resource. Such machines are autonomous, and their performance must be bounded using entropy production for all reservoirs, not only a hot and cold pair.

  • Quoting an efficiency without specifying heat-current sign conventions.
  • Calling a squeezed, inverted, or chemically biased reservoir a thermal bath without accounting for the extra resource.
  • Using local dissipators in an interacting working medium without checking thermodynamic consistency.
  • Treating coherence as a universal advantage rather than a protocol-dependent resource.
  • Ignoring the work storage device or load in an autonomous engine.
  • Comparing finite-time power to quasistatic Carnot efficiency as if they were optimized by the same protocol.
  • Forgetting that refrigeration has a coefficient of performance, not an efficiency.
  • Counting heat leaks as useful heat input.

Before trusting a quantum heat-engine model, identify:

  • the working medium and its Hamiltonian;
  • all reservoirs and their temperatures or nonequilibrium resources;
  • the work repository, load, drive, or battery;
  • the sign convention for heat and work;
  • whether the model is cyclic, periodically driven, or autonomous steady state;
  • whether each dissipator satisfies detailed balance for the Hamiltonian being used;
  • whether strong coupling or finite reservoirs alter the energy boundary;
  • whether measurements, feedback, or records enter the entropy balance;
  • the operating regime: engine, refrigerator, heat pump, or dissipator.

This checklist is mundane in the best way: it prevents impressive formulas from being attached to incomplete machines.

Assume a cyclic engine absorbs Qh>0Q_h\gt0 from a hot bath and has Qc<0Q_c\lt0 for the cold bath, with no other reservoirs. Using

Σ=−QhTh−QcTc≥0,\Sigma = - \frac{Q_h}{T_h} - \frac{Q_c}{T_c} \ge0,

derive the Carnot efficiency bound.

Solution

The entropy-production inequality gives

−QcTc≥QhTh.- \frac{Q_c}{T_c} \ge \frac{Q_h}{T_h}.

Multiplying by TcT_c,

−Qc≥TcThQh.-Q_c \ge \frac{T_c}{T_h}Q_h .

Thus

Qc≤−TcThQh.Q_c \le - \frac{T_c}{T_h}Q_h .

The extracted work is Wext=Qh+QcW_{\mathrm{ext}}=Q_h+Q_c, so

η=WextQh=1+QcQh≤1−TcTh.\eta = \frac{W_{\mathrm{ext}}}{Q_h} = 1+\frac{Q_c}{Q_h} \le 1-\frac{T_c}{T_h}.

For the two-level Otto engine, show that the engine condition ph>pcp_h\gt p_c is equivalent to

βhωh<βcωc.\beta_h\omega_h \lt \beta_c\omega_c .
Solution

The excited-state population of a thermal two-level system is

p(β,ω)=11+eβℏω.p(\beta,\omega) = \frac{1} {1+e^{\beta\hbar\omega}}.

This function decreases monotonically with βω\beta\omega. Therefore

ph>pcp_h\gt p_c

is equivalent to

βhℏωh<βcℏωc.\beta_h\hbar\omega_h \lt \beta_c\hbar\omega_c .

Since ℏ>0\hbar\gt0, this is the same as

βhωh<βcωc.\beta_h\omega_h \lt \beta_c\omega_c .

Using

Qh=ℏωh(ph−pc),Qc=ℏωc(pc−ph),Q_h=\hbar\omega_h(p_h-p_c), \qquad Q_c=\hbar\omega_c(p_c-p_h),

derive the ideal Otto efficiency.

Solution

The extracted work is

Wext=Qh+Qc=ℏ(ωh−ωc)(ph−pc).W_{\mathrm{ext}} = Q_h+Q_c = \hbar(\omega_h-\omega_c)(p_h-p_c).

The efficiency is

η=WextQh=ℏ(ωh−ωc)(ph−pc)ℏωh(ph−pc).\eta = \frac{W_{\mathrm{ext}}}{Q_h} = \frac{ \hbar(\omega_h-\omega_c)(p_h-p_c) } { \hbar\omega_h(p_h-p_c) }.

In the engine regime ph≠pcp_h\ne p_c, so

η=1−ωcωh.\eta = 1-\frac{\omega_c}{\omega_h}.

Assume the hot transition gives p3/p1=e−βhℏωhp_3/p_1=e^{-\beta_h\hbar\omega_h} and the cold transition gives p2/p1=e−βcℏωcp_2/p_1=e^{-\beta_c\hbar\omega_c}. Derive the condition for population inversion on the work transition.

Solution

Population inversion on the work transition means

p3>p2.p_3\gt p_2 .

Using the two ratios,

p3p1>p2p1\frac{p_3}{p_1} \gt \frac{p_2}{p_1}

becomes

e−βhℏωh>e−βcℏωc.e^{-\beta_h\hbar\omega_h} \gt e^{-\beta_c\hbar\omega_c}.

Taking logarithms and reversing the sign after multiplying by −1-1 gives

βhωh<βcωc.\beta_h\omega_h \lt \beta_c\omega_c .

This is the same condition that appears in the two-level Otto engine.

  • H. E. D. Scovil and E. O. Schulz-DuBois, “Three-level masers as heat engines,” Physical Review Letters 2, 262-263 (1959).
  • J. E. Geusic, E. O. Schulz-DuBois, R. W. De Grasse, and H. E. D. Scovil, “Quantum equivalent of the Carnot cycle,” Physical Review 156, 343-351 (1967).
  • R. Alicki, “The quantum open system as a model of the heat engine,” Journal of Physics A: Mathematical and General 12, L103-L107 (1979).
  • R. Kosloff, “Quantum thermodynamics: A dynamical viewpoint,” Entropy 15, 2100-2128 (2013).
  • A. Levy and R. Kosloff, “Quantum absorption refrigerator,” Physical Review Letters 108, 070604 (2012).
  • R. Kosloff and A. Levy, “Quantum heat engines and refrigerators: Continuous devices,” Annual Review of Physical Chemistry 65, 365-393 (2014).
  • D. Gelbwaser-Klimovsky, W. Niedenzu, and G. Kurizki, “Thermodynamics of quantum systems under dynamical control,” Advances in Atomic, Molecular, and Optical Physics 64, 329-407 (2015).
  • J. Goold, M. Huber, A. Riera, L. del Rio, and P. Skrzypczyk, “The role of quantum information in thermodynamics: A topical review,” Journal of Physics A: Mathematical and Theoretical 49, 143001 (2016).
  • M. T. Mitchison, “Quantum thermal absorption machines: Refrigerators, engines and clocks,” Contemporary Physics 60, 164-187 (2019).