Skip to content

Ergotropy and Passive States

Ergotropy is the maximum work that can be extracted from a quantum state by a cyclic unitary process. The Hamiltonian at the beginning and end is the same. The controller may act coherently on the system, but there is no heat bath, no measurement feedback, and no change in the system Hamiltonian after the cycle closes.

The concept is useful because it separates a sharp finite-system question from broader thermodynamic work statements. Given only a pair (ρ,H)(\rho,H) and full coherent control, which part of the state’s energy is ordered enough to be converted into work?

It is also easy to overinterpret. Ergotropy is not the same as nonequilibrium free energy, not the same as the random work distribution from a two-point measurement protocol, and not automatically available under thermal operations that lack an external phase reference.

Let a system have time-independent Hamiltonian

H=∑nEn∣En⟩⟨En∣,H = \sum_n E_n |E_n\rangle\langle E_n|,

with energy levels ordered as

E1≤E2≤⋯ .E_1\le E_2\le\cdots .

A cyclic unitary protocol begins with state ρ\rho, applies a unitary UU, and ends with the same Hamiltonian HH. The final state is

ρ′=UρU†.\rho' = U\rho U^\dagger .

Using the convention that positive work is extracted from the system, the extracted work is

Wext(U)=Tr⁡(ρH)−Tr⁡(UρU†H).W_{\mathrm{ext}}(U) = \operatorname{Tr}(\rho H) - \operatorname{Tr}(U\rho U^\dagger H).

The unitary does not change the eigenvalues of ρ\rho. It can only rearrange the state’s eigenvectors relative to the energy eigenbasis. Therefore the minimization problem is:

minimize Tr⁡(UρU†H)over all unitaries U.\text{minimize } \operatorname{Tr}(U\rho U^\dagger H) \quad \text{over all unitaries } U .

The ergotropy of (ρ,H)(\rho,H) is the resulting maximum:

W(ρ,H)=Tr⁡(ρH)−min⁡UTr⁡(UρU†H).\mathcal W(\rho,H) = \operatorname{Tr}(\rho H) - \min_U \operatorname{Tr}(U\rho U^\dagger H).

Here W≥0\mathcal W\ge0. It vanishes exactly when ρ\rho is already passive with respect to HH.

Write the spectral decomposition of the state as

ρ=∑krk∣rk⟩⟨rk∣,r1≥r2≥⋯ .\rho = \sum_k r_k |r_k\rangle\langle r_k|, \qquad r_1\ge r_2\ge\cdots .

For a nondegenerate Hamiltonian with ordered energies, the unitary that minimizes the final energy maps the largest eigenvalue of ρ\rho to the lowest energy eigenstate, the next largest eigenvalue to the next lowest energy eigenstate, and so on. The corresponding passive state is

πρ=∑krk∣Ek⟩⟨Ek∣.\pi_\rho = \sum_k r_k |E_k\rangle\langle E_k|.

Thus

W(ρ,H)=Tr⁡(ρH)−Tr⁡(πρH).\mathcal W(\rho,H) = \operatorname{Tr}(\rho H) - \operatorname{Tr}(\pi_\rho H).

This is the quantum version of sorting. The spectrum of ρ\rho is fixed by unitary evolution, and the final energy is smallest when probability weight is placed as low in the energy ladder as the spectrum permits.

Degeneracies require a small convention. Inside a degenerate energy subspace, unitary rotations do not change the energy. Inside a degenerate eigenspace of ρ\rho, the eigenvectors are not unique. The passive rearrangement is therefore unique only up to such rotations, but its energy and ergotropy are well defined.

A state σ\sigma is passive with respect to HH if no cyclic unitary can extract work:

Tr⁡(σH)≤Tr⁡(UσU†H)for every unitary U.\operatorname{Tr}(\sigma H) \le \operatorname{Tr}(U\sigma U^\dagger H) \qquad \text{for every unitary } U .

For a nondegenerate Hamiltonian, this is equivalent to two conditions:

  • σ\sigma is diagonal in the energy eigenbasis;
  • the populations are nonincreasing with energy.

Explicitly, if

σ=∑npn∣En⟩⟨En∣,\sigma = \sum_n p_n |E_n\rangle\langle E_n|,

then passivity requires

Ei≤Ej⟹pi≥pj.E_i\le E_j \quad\Longrightarrow\quad p_i\ge p_j .

If a higher energy level is more populated than a lower one, a unitary swap lowers the mean energy and extracts work. Such a population inversion is active. If the state has energy-basis coherence, a suitable unitary can usually lower the energy even when the diagonal populations alone look passive.

Passivity is a statement about a specified Hamiltonian. The same density operator can be passive for one Hamiltonian and active for another.

Passivity of one copy is weaker than thermodynamic equilibrium. A state σ\sigma is completely passive if every tensor power

σ⊗N\sigma^{\otimes N}

is passive with respect to the noninteracting Hamiltonian

H(N)=∑j=1NI⊗⋯⊗H⊗⋯⊗I.H^{(N)} = \sum_{j=1}^N I\otimes\cdots\otimes H\otimes\cdots\otimes I .

The distinction matters because correlations among many copies can unlock work even when each copy is individually passive. In finite-dimensional systems under standard assumptions, Gibbs states

τβ=e−βHZ,β≥0,\tau_\beta = \frac{e^{-\beta H}} {Z}, \qquad \beta\ge0,

are completely passive. Conversely, complete passivity singles out thermal equilibrium states, with appropriate qualifications for infinite systems and zero-temperature sectors.

This is one reason ergotropy should not be read as the entire second law. A single passive nonthermal state may contain no one-copy cyclic-unitary work, while many copies can still contain an extractable resource.

At fixed temperature TT, the nonequilibrium free energy is

FT(ρ)=Tr⁡(ρH)−TS(ρ),F_T(\rho) = \operatorname{Tr}(\rho H) - T S(\rho),

where

S(ρ)=−kBTr⁡(ρln⁡ρ).S(\rho) = -k_B\operatorname{Tr}(\rho\ln\rho).

The difference

FT(ρ)−FT(τβ)=kBT D(ρ∥τβ)F_T(\rho)-F_T(\tau_\beta) = k_B T\,D(\rho\Vert\tau_\beta)

measures the maximum reversible work available in ideal bath-assisted transformations, under the assumptions behind equilibrium thermodynamics.

Ergotropy asks a different question. It allows a cyclic unitary on the system but no heat bath and no entropy disposal. Because unitary evolution preserves S(ρ)S(\rho) and the eigenvalues of ρ\rho, the only available work comes from arranging those eigenvalues more favorably relative to HH.

A thermal state has zero ergotropy:

W(τβ,H)=0(β≥0).\mathcal W(\tau_\beta,H) = 0 \qquad (\beta\ge0).

This does not mean the state has zero thermodynamic free energy. It means that no work can be extracted from that state alone by a cyclic unitary.

Energy-basis coherence can contribute to ergotropy if the controller is allowed to implement arbitrary unitaries. For example, a pure state

∣ψ⟩=1−p ∣g⟩+eiϕp ∣e⟩|\psi\rangle = \sqrt{1-p}\,|g\rangle + e^{i\phi}\sqrt p\,|e\rangle

of a qubit with Hamiltonian H=ϵ∣e⟩⟨e∣H=\epsilon |e\rangle\langle e| has mean energy ϵp\epsilon p. Since any pure state can be unitarily mapped to the ground state, its ergotropy is

W(∣ψ⟩⟨ψ∣,H)=ϵp.\mathcal W(|\psi\rangle\langle\psi|,H) = \epsilon p .

The phase ϕ\phi affects which unitary performs the extraction, even though the maximum extracted work depends only on the mean energy for this pure qubit example.

Operational restrictions change the statement. If allowed operations must commute with total energy, or if no phase reference is available, coherence between distinct energy eigenspaces is not freely convertible into mechanical work. In resource-theoretic thermodynamics, population imbalance and coherence are often separate resources. See Thermal Operations Preview for that viewpoint.

The useful rule is:

  • ergotropy with arbitrary cyclic unitaries counts coherent control as part of the work-extraction apparatus;
  • thermal operations without an external coherence reference do not treat energy-basis coherence as automatically extractable work.

Let

H=ϵ∣e⟩⟨e∣,ϵ>0,H = \epsilon |e\rangle\langle e|, \qquad \epsilon>0,

with state

ρ=pg∣g⟩⟨g∣+pe∣e⟩⟨e∣,pg+pe=1.\rho = p_g |g\rangle\langle g| + p_e |e\rangle\langle e|, \qquad p_g+p_e=1 .

The mean energy is

E(ρ)=ϵpe.E(\rho) = \epsilon p_e .

If pe≤pgp_e\le p_g, the state is passive: the larger population is already on the lower energy level. Its ergotropy is zero.

If pe>pgp_e>p_g, the state is inverted. A unitary swap places population pep_e in ∣g⟩|g\rangle and pgp_g in ∣e⟩|e\rangle. The final energy is ϵpg\epsilon p_g, so

W(ρ,H)=ϵ(pe−pg)=ϵ(2pe−1).\mathcal W(\rho,H) = \epsilon(p_e-p_g) = \epsilon(2p_e-1).

Combining the two cases,

W(ρ,H)=ϵmax⁡(0,2pe−1).\mathcal W(\rho,H) = \epsilon\max(0,2p_e-1).

This simple formula is a good check on sign conventions: population inversion is active, a positive-temperature thermal qubit is passive.

Consider a three-level system with

E1≤E2≤E3E_1\le E_2\le E_3

and a diagonal state whose populations are

(p1,p2,p3)=(0.2,0.5,0.3).(p_1,p_2,p_3) = (0.2,0.5,0.3).

The state is not passive because p2>p1p_2>p_1 while E2≥E1E_2\ge E_1. The passive rearrangement sorts the same probabilities in decreasing order over increasing energy:

(0.5,0.3,0.2).(0.5,0.3,0.2).

The ergotropy is

W=∑n=13pnEn−(0.5E1+0.3E2+0.2E3).\mathcal W = \sum_{n=1}^3 p_n E_n - \left( 0.5E_1+0.3E_2+0.2E_3 \right).

No entropy has been removed. The work comes only from undoing the population disorder relative to the energy ordering.

A harmonic oscillator thermal state is passive:

τβ=e−βℏωa†aTr⁡(e−βℏωa†a).\tau_\beta = \frac{e^{-\beta\hbar\omega a^\dagger a}} {\operatorname{Tr}(e^{-\beta\hbar\omega a^\dagger a})}.

If the oscillator is squeezed, the resulting state can have the same entropy as a thermal state but a larger mean energy. Under suitable coherent control, the squeezing can be undone and part of that excess energy can be extracted as work. This is why squeezed reservoirs and squeezed states are often described as carrying an ordered, work-like resource in quantum thermodynamic models.

The exact ergotropy depends on which Hamiltonian defines the energy and which unitaries the controller can implement. In open-system settings, be careful to distinguish a squeezed bath as an engineered nonequilibrium reservoir from a thermal bath at a higher temperature.

Ergotropy is a state function of the pair (ρ,H)(\rho,H) under an allowed control class. It is a maximum extractable average work.

A work distribution is different. In a two-point measurement protocol, one samples initial and final energies and assigns

W=Emτ−En0.W = E_m^\tau-E_n^0 .

That random variable depends on projective measurements, driving, and a sign convention. It can obey fluctuation relations such as Jarzynski and Crooks identities when their assumptions hold. See Two-Point Measurement Scheme, Work Distributions, and Jarzynski and Crooks Relations.

The two languages meet in carefully designed protocols, but they should not be collapsed into one slogan. Ergotropy optimizes over unitaries for a given state; fluctuation relations compare probability distributions generated by specified protocols.

  • Calling every nonthermal state active. Some nonthermal states are passive for one copy.
  • Forgetting that passivity is defined relative to a Hamiltonian.
  • Confusing ergotropy with nonequilibrium free energy.
  • Treating energy-basis coherence as freely extractable when the allowed operations lack a phase reference.
  • Ignoring degeneracies when identifying the passive rearrangement.
  • Assuming one-copy passivity implies complete passivity.
  • Using ergotropy formulas for open driven systems without specifying which degrees of freedom are included in the cyclic unitary.

When using ergotropy in a calculation, specify:

  • the Hamiltonian HH defining energy;
  • whether the Hilbert space is finite-dimensional or requires domain care;
  • the allowed control unitaries;
  • whether the process is cyclic in the Hamiltonian;
  • whether positive work means work extracted or work supplied;
  • whether the state has energy-basis coherence and whether a phase reference is available;
  • whether one copy or many copies are being considered.

These choices determine whether a number called “extractable work” has a well-defined operational meaning.

For the diagonal qubit state

ρ=(1−p)∣g⟩⟨g∣+p∣e⟩⟨e∣,H=ϵ∣e⟩⟨e∣,\rho = (1-p)|g\rangle\langle g| + p|e\rangle\langle e|, \qquad H=\epsilon |e\rangle\langle e|,

derive

W(ρ,H)=ϵmax⁡(0,2p−1).\mathcal W(\rho,H) = \epsilon\max(0,2p-1).
Solution

The eigenvalues of ρ\rho are pp and 1−p1-p. The passive state places the larger eigenvalue on the ground state and the smaller eigenvalue on the excited state.

If p≤1/2p\le1/2, then 1−p≥p1-p\ge p and the state is already passive. Thus W=0\mathcal W=0.

If p>1/2p>1/2, the passive state has excited-state population 1−p1-p. The initial and passive energies are

Ein=ϵp,Epass=ϵ(1−p).E_{\mathrm{in}} = \epsilon p, \qquad E_{\mathrm{pass}} = \epsilon(1-p).

Therefore

W=ϵp−ϵ(1−p)=ϵ(2p−1).\mathcal W = \epsilon p-\epsilon(1-p) = \epsilon(2p-1).

Combining both regimes gives the stated formula.

Let

∣ψ⟩=1−p ∣g⟩+eiϕp ∣e⟩,H=ϵ∣e⟩⟨e∣.|\psi\rangle = \sqrt{1-p}\,|g\rangle + e^{i\phi}\sqrt p\,|e\rangle, \qquad H=\epsilon |e\rangle\langle e|.

Assuming arbitrary cyclic unitaries are available, find the ergotropy.

Solution

The state is pure, so its eigenvalues are 11 and 00. A unitary can map ∣ψ⟩|\psi\rangle to the ground state ∣g⟩|g\rangle. The passive state is therefore ∣g⟩⟨g∣|g\rangle\langle g|, whose energy is zero.

The initial energy is

⟨ψ∣H∣ψ⟩=ϵp.\langle\psi|H|\psi\rangle = \epsilon p.

Thus

W=ϵp.\mathcal W = \epsilon p.

The phase ϕ\phi affects the extracting unitary, but not the maximum value in this two-level pure-state example.

Show that a finite-dimensional Gibbs state

τβ=e−βHZ\tau_\beta = \frac{e^{-\beta H}}{Z}

with β≥0\beta\ge0 is passive.

Solution

Let H∣En⟩=En∣En⟩H|E_n\rangle=E_n|E_n\rangle, with energies ordered increasingly. The Gibbs populations are

pn=e−βEnZ.p_n = \frac{e^{-\beta E_n}}{Z}.

If Ei≤EjE_i\le E_j and β≥0\beta\ge0, then

e−βEi≥e−βEj.e^{-\beta E_i} \ge e^{-\beta E_j}.

Therefore pi≥pjp_i\ge p_j. The Gibbs state is diagonal in the energy basis and its populations decrease with increasing energy, so it is passive.

Consider a three-level system with equally spaced energies 0,ϵ,2ϵ0,\epsilon,2\epsilon and diagonal populations

(p0,p1,p2)=(0.5,0.3,0.2).(p_0,p_1,p_2) = (0.5,0.3,0.2).

Is the state passive? Is it necessarily a Gibbs state?

Solution

The populations decrease with energy, so the state is passive for one copy.

It is not necessarily Gibbs. For a Gibbs state with equally spaced levels, the ratios would obey

p1p0=p2p1=e−βϵ.\frac{p_1}{p_0} = \frac{p_2}{p_1} = e^{-\beta\epsilon}.

Here

p1p0=0.6,p2p1=23.\frac{p_1}{p_0} = 0.6, \qquad \frac{p_2}{p_1} = \frac{2}{3}.

The ratios are unequal, so this specific state is passive but not Gibbs. This illustrates why passivity and thermal equilibrium are different statements.

  • W. Pusz and S. L. Woronowicz, “Passive states and KMS states for general quantum systems,” Communications in Mathematical Physics 58, 273-290 (1978).
  • A. Lenard, “Thermodynamical proof of the Gibbs formula for elementary quantum systems,” Journal of Statistical Physics 19, 575-586 (1978).
  • A. E. Allahverdyan, R. Balian, and Th. M. Nieuwenhuizen, “Maximal work extraction from finite quantum systems,” Europhysics Letters 67, 565-571 (2004).
  • R. Alicki and M. Fannes, “Entanglement boost for extractable work from ensembles of quantum batteries,” Physical Review E 87, 042123 (2013).
  • K. V. Hovhannisyan, M. Perarnau-Llobet, M. Huber, and A. Acin, “Entanglement generation is not necessary for optimal work extraction,” Physical Review Letters 111, 240401 (2013).
  • 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).
  • F. Binder, L. A. Correa, C. Gogolin, J. Anders, and G. Adesso, eds., Thermodynamics in the Quantum Regime: Fundamental Aspects and New Directions, Springer (2018).