---
title: Classical Ergotropy Extraction Driving
url: https://www.emergentmind.com/topics/classical-ergotropy-extraction-driving
type: topic
---

# Classical Ergotropy Extraction Driving

Searching arXiv for recent and foundational papers on classical ergotropy extraction and related formulations.
Classical ergotropy extraction driving is the problem of identifying a time-dependent perturbation that extracts the maximal amount of work from a thermally isolated classical Hamiltonian system prepared in a nonequilibrium distribution. In the formulation developed for classical systems, the maximal extractable work is the classical ergotropy, defined as the gap between the initial mean energy and the minimal final energy reachable by volume-preserving dynamics. Under an ergodic assumption, the optimal driving is given by a two-step quench–adiabat protocol: an instantaneous quench to a Hamiltonian for which the initial distribution is already passive, followed by an adiabatic return to the original Hamiltonian. This construction provides a classical analogue of the quantum passive-state rearrangement problem and supports a coherent–incoherent decomposition of ergotropy closely paralleling the quantum case [2508.12797].

## 1. Problem statement and definition of classical ergotropy

A classical system is specified by canonical coordinates $z=(q,p)$ and an unperturbed Hamiltonian $H_0(z)$. The initial condition is a nonequilibrium phase-space distribution $\rho_0(z)\ge 0$ with $\int dz\,\rho_0(z)=1$, and the system is then isolated from heat baths. One applies a time-dependent perturbation $V(z,t)$ over a finite interval $t\in[0,\tau]$, so that the full Hamiltonian is
$$
H(z,t)=H_0(z)+V(z,t).
$$
Because the dynamics is Hamiltonian, Liouville’s theorem enforces phase-space volume preservation. The extracted work is the decrease in the average energy with respect to the original Hamiltonian:
$$
W = E_i - E_f = \int dz\,\rho_0(z)\,H_0(z)-\int dz\,\rho_\tau(z)\,H_0(z).
$$
The central optimization problem is therefore to determine which driving $V(z,t)$ maximizes $W$ [2508.12797].

Under very general arguments identified with Gardner’s theorem, the minimal reachable final energy $\breve E_c$ over all volume-preserving maps coincides with the Gardner free energy. The classical ergotropy is then
$$
\mathcal E_c \equiv \max_{\rm drivings} W
= \int dz\,\rho_0(z)\,H_0(z)-\breve E_c.
$$
This definition places classical ergotropy on the same conceptual footing as quantum ergotropy, where the maximum cyclic-unitary work is determined by rearranging a fixed spectrum into a passive ordering [2508.12797].

A complementary line of work formulates classical ergotropy in terms of classical relative entropies. In that approach, for an initial distribution on an energy shell $A$, a final distribution on shell $B$, and a Gibbs reference distribution $p_B^{\rm eq}$, one defines
$$
\Omega_{\rm cl}
=\frac{1}{\beta}\Bigl[
D\!\bigl(p_{B|A}\|p_B^\mathrm{eq}\bigr)
-
D\!\bigl(p_A\|p_B^\mathrm{eq}\bigr)
\Bigr].
$$
This makes explicit that classical ergotropy quantifies the additional nonequilibrium structure produced by Liouville evolution relative to thermal equilibrium on the target shell [2103.10850]. A plausible implication is that the Gardner-style rearrangement picture and the relative-entropy picture are complementary descriptions of the same optimization problem.

## 2. Passive rearrangement and the Gardner construction

The optimal classical rearrangement is expressed through two cumulative phase-space functions. The first is
$$
R(\sigma)=\int dz\,\Theta\bigl(\rho_0(z)-\sigma\bigr),
$$
which is strictly decreasing in $\sigma$. The second is
$$
\Omega_0(E)=\int dz\,\Theta\bigl(E-H_0(z)\bigr),
$$
which is strictly increasing in $E$. Gardner’s construction then identifies the minimal-energy reachable distribution as a function of the unperturbed energy alone,
$$
\rho_1(z)=g\!\bigl(H_0(z)\bigr), \qquad
g^{-1}(\sigma)=E_0\bigl(R(\sigma)\bigr),
$$
with energy
$$
\breve E_c
= \int dz\,\rho_1(z)\,H_0(z)
= \int_0^\infty d\Phi\,E_0(\Phi)\,R^{-1}(\Phi).
$$
Accordingly,
$$
\mathcal E_c
=
\int dz\,\rho_0(z)\,H_0(z)
-\int_0^\infty d\Phi\,E_0(\Phi)\,R^{-1}(\Phi).
$$
Equivalently,
$$
\mathcal E_c
=
\int dz\,\rho_0(z)\,\bigl[H_0(z)-H_0^{\rm pass}(z)\bigr],
$$
where $H_0^{\rm pass}(z)$ is the one-dimensional rearrangement of $H_0$ that makes the final distribution a decreasing function of energy [2508.12797].

This is the classical passive-state construction. Its significance is structural: the minimizing state is not obtained by thermalization, but by a measure-preserving sorting of phase-space weight so that larger values of $\rho_0$ are assigned to lower-energy regions of $H_0$. That is the direct classical counterpart of the quantum passive state, where the largest eigenvalues are aligned with the lowest eigenenergies [2508.12797].

The relative-entropy formulation sharpens the same point in a different language. In that framework, any deviation from uniformity on a constant-energy surface stores active energy, and one can write
$$
\Omega_{\rm cl}
= \int d\Gamma'\,\varphi_B(\Gamma')\,E_B(\Gamma'),
$$
where $\varphi_B(\Gamma')$ is the inhomogeneity that develops on the final shell. This suggests that passive rearrangement and energy-shell homogenization are dual descriptions of the same resource structure [2103.10850].

## 3. The quench–adiabat driving protocol

The central result for classical ergotropy extraction is the quench–adiabat, or QA, protocol. It consists of two steps [2508.12797].

First, one performs an instantaneous quench at $t=0$,
$$
H_0(z)\longrightarrow H_1(z)=f\!\bigl(\rho_0(z)\bigr),
\qquad f\ \text{strictly decreasing}.
$$
Since the quench is instantaneous, the phase-space distribution remains $\rho_0$. However, by construction, $\rho_0$ becomes a decreasing function of the new Hamiltonian $H_1$, so it is passive relative to $H_1$.

Second, one performs an adiabatic return from $H_1$ to $H_0$ over $t\in[0,\tau]$. Under the ergodic hypothesis for each frozen Hamiltonian $H(z,t)$, the enclosed phase-space volume
$$
\Omega_t(E)=\int dz\,\Theta\bigl(E-H(z,t)\bigr)
$$
is an adiabatic invariant. Each phase-space point is therefore transported from a level set of $H_1$ to a level set of $H_0$ with the same volume label. The result is exactly the Gardner passive rearrangement, so the final energy is $\breve E_c$ and the extracted work is maximal [2508.12797].

The protocol has a direct operational interpretation. The quench “locks in” the initial distribution as a passive state of an auxiliary Hamiltonian, and the adiabatic return uses Liouville volume preservation together with adiabatic invariance to map that passive ordering back onto the original Hamiltonian. The paper explicitly states that this is the classical analogue of the quantum “reorder-and-adiabatic-ramp” proof [2508.12797].

This protocol should be distinguished from work on classically driven quantum systems, where a classical field controls the ergotropy of a quantum emitter or battery. In those settings, “classical drive” denotes an externally applied field acting on a quantum system, whereas in classical ergotropy extraction driving the system itself is classical and the objective is the optimal Hamiltonian deformation implementing the passive rearrangement [2410.23589].

## 4. Why the QA protocol is optimal

The optimality proof proceeds by introducing the phase-space volume function of the quenched Hamiltonian,
$$
\Omega_1(E)=\int dz\,\Theta\bigl(E-H_1(z)\bigr),
$$
and the shell probability density
$$
P_1(\Phi)
=
\int dz\,\rho_0(z)\,\delta\!\bigl[\Phi-\Omega_1\bigl(H_1(z)\bigr)\bigr].
$$
Because $\rho_0=f^{-1}(H_1)$, one obtains
$$
P_1(\Phi)=f^{-1}\!\bigl(E_1(\Phi)\bigr),
\qquad E_1=\Omega_1^{-1}.
$$
Adiabatic invariance of $\Omega$ implies that this same $P_1(\Phi)$ becomes the probability density over shells of the original Hamiltonian after the return. The final mean energy is therefore
$$
E_f=\int_0^\infty d\Phi\,P_1(\Phi)\,E_0(\Phi).
$$
Choosing the special monotone $f=R$ makes $\Omega_1(x)=x$ and $P_1(\Phi)=R^{-1}(\Phi)$, which yields
$$
E_f
=
\int_0^\infty d\Phi\,R^{-1}(\Phi)\,E_0(\Phi)
=
\breve E_c.
$$
The QA protocol thus exactly attains the Gardner minimum [2508.12797].

The same result can be re-expressed in action–angle variables when the unperturbed dynamics admits coordinates $(I,\theta)$ with $H_0=H_0(I)$. Adiabatic invariance then states that the action is approximately conserved during the slow stage, so the final map permutes action “rings” in the way required to minimize the average energy. This is again described as the classical counterpart of the quantum reorder-and-ramp argument [2508.12797].

A plausible implication is that the proof isolates the genuinely dynamical content of classical ergotropy extraction: the problem is not merely combinatorial sorting of probabilities, but the physical implementation of that sorting by Hamiltonian flows constrained by Liouville’s theorem and adiabatic invariance.

## 5. Coherent and incoherent parts of classical ergotropy

The 2025 treatment shows that classical ergotropy splits into coherent and incoherent contributions, just as in the quantum case [2508.12797]. The relevant classical dephasing, or homogenization, operator acts on the energy shells of $H_0$:
$$
\mathcal D[\rho](z)
=
\int dz'\,\rho(z')\,
\delta\!\bigl[
\Omega_0\!\bigl(H_0(z')\bigr)-\Omega_0\!\bigl(H_0(z)\bigr)
\bigr].
$$
This map preserves the mean energy while erasing the fine structure of $\rho$ along each iso-$H_0$ shell. One then defines
$$
\mathcal E_c^{\rm inc}[\rho_0]
=
\mathcal E_c\!\bigl[\mathcal D[\rho_0]\bigr],
$$
and
$$
\mathcal E_c^{\rm coh}[\rho_0]
=
\mathcal E_c[\rho_0]-\mathcal E_c^{\rm inc}[\rho_0].
$$
The incoherent part is therefore the ergotropy stored in radial shell populations, while the coherent part is the additional ergotropy contained in angular inhomogeneities within each shell [2508.12797].

The paper also gives an information-theoretic identity:
$$
\beta\,\mathcal E_c^{\rm coh}[\rho_0]
=
C[\rho_0]
+
D\!\bigl[\mathcal P(\mathcal D[\rho_0])\|\rho_\beta\bigr]
-
D\!\bigl[\mathcal P(\rho_0)\|\rho_\beta\bigr],
$$
where $\mathcal P(\rho)$ is the passive rearrangement of $\rho$, $D[p\|q]$ is the Kullback–Leibler divergence, and
$$
C[\rho_0]=D[\rho_0\|\mathcal D[\rho_0]]
$$
measures deviation from the dephased distribution, identified as “classical coherence” [2508.12797].

This decomposition is closely related to the earlier relative-entropy formulation of classical ergotropy, in which energy-shell inhomogeneity is the active resource and the classical limit emerges from the geometric formulation of quantum mechanics when probability densities are supported on energy-eigenstate poles [2103.10850]. It is also consistent with the distinction between incoherent and coherent ergotropy extraction in isolated quantum settings, where incoherent extraction corresponds to level permutations and coherent extraction addresses off-diagonal structure in the energy basis [2111.03116]. The classical case replaces Hilbert-space coherence by phase-space inhomogeneity on constant-energy surfaces.

## 6. Examples, topology, and the role of ergodicity

A canonical example is the harmonic oscillator with displaced Gaussian initial distribution
$$
H_0(q,p)=\tfrac12(p^2+q^2),\qquad
\rho_0(q,p)=\frac{1}{2\pi\sigma}
\exp\!\Bigl[-\frac{(q-q_0)^2+p^2}{2\sigma}\Bigr].
$$
Under the unperturbed flow, this Gaussian rotates and does not become passive. The QA protocol first quenches to
$$
H_1(q,p)=\tfrac12\bigl[p^2+(q-q_0)^2\bigr]
=f\bigl(\rho_0(q,p)\bigr),
$$
with
$$
f(x)=-\sigma\ln(2\pi\sigma\,x),
$$
so that the state is stationary and passive relative to $H_1$. One then adiabatically moves the potential minimum from $q=q_0$ back to $q=0$, producing the final passive Gibbs-form distribution
$$
\rho_1(q,p)=\frac{1}{2\pi\sigma}e^{-(p^2+q^2)/(2\sigma)},
$$
and the extracted work is
$$
\mathcal E_c
=
\langle H_0\rangle_{\rho_0}-\langle H_0\rangle_{\rho_1}
=
\frac{q_0^2}{2}.
$$
Thus the ergotropy equals the potential-energy drop of the displaced center [2508.12797].

The analysis becomes subtler for topologically nontrivial initial distributions, such as a ring-shaped or “Mexican-hat” distribution with a central valley. In that case, the monotone $f(\rho_0)$ generates an auxiliary Hamiltonian $H_1$ with level sets that are not simply connected. A purely adiabatic global return cannot permute all such structures in a single step. The proposed resolution is to choose an implementable $H_1$ sharing the same level-set topology, perform the adiabatic return, and then supplement it with a finite sequence of small non-adiabatic “ring-permutation” drivings that permute subregions and approach the true passive ordering arbitrarily closely. The paper states that this combination of QA plus controlled mixing fully recovers $\mathcal E_c$ [2508.12797].

The role of ergodicity is central. The adiabatic invariance of enclosed phase-space volume $\Omega_t$ relies on ergodic motion on each energy shell, in the sense discussed through Khinchin’s ergodic-theorem-based proofs and extensions by Hertz, Jarzynski, and Ott [2508.12797]. If the instantaneous Hamiltonians are not ergodic, volume invariance may fail for a global adiabatic deformation. The same work nonetheless emphasizes two extensions: one may often select an $H_1$ with the same shell topology as $H_0$, and when shells are multiply connected one can replace the purely adiabatic step by topologically targeted field pulses that permute small subregions. This suggests that ergodicity is sufficient for the clean QA construction, but not always necessary for near-optimal or fully optimal extraction in more structured phase spaces.

## 7. Relation to quantum ergotropy and broader significance

The classical theory is explicitly framed as an analogue of the solved quantum ergotropy problem. In finite-dimensional quantum systems, ergotropy is the maximal cyclic work
$$
\mathcal E(\rho,H)=\operatorname{Tr}[\rho H]-\operatorname{Tr}[P_\rho H],
$$
where the passive state $P_\rho$ is obtained by assigning the largest eigenvalues of $\rho$ to the lowest energy levels of $H$ [2202.05050]. The classical QA protocol implements the corresponding rearrangement through Hamiltonian deformation rather than unitary diagonalization and eigenvalue sorting [2508.12797].

This correspondence has several consequences. First, it strengthens the claim that there is a unified passive-state principle underlying both classical and quantum ergotropy. Second, it clarifies that coherent versus incoherent contributions are not uniquely quantum; the classical setting admits an analogous decomposition once “coherence” is interpreted as inhomogeneity on energy shells [2508.12797; 2103.10850]. Third, it places classical ergotropy extraction alongside work on practical ergotropy estimation and control in quantum platforms, where passive states are reached through feedback or externally driven protocols, such as FQErgo for unknown quantum states [2409.04087]. The conceptual parallel is that both cases seek constructive procedures that attain the passive rearrangement rather than merely proving its existence.

A common misconception is to identify “classical ergotropy extraction driving” with the use of a classical field to control a quantum battery. The literature distinguishes these notions sharply. In one case, the controlled object is a classical Hamiltonian system and the optimal driving is the QA protocol [2508.12797]. In the other, a classical field acts on a quantum system, such as a two-level emitter or battery, and the analysis concerns the evolution of quantum ergotropy under continuous or pulsed driving [2410.23589; 2302.12279]. The shared term “classical” refers to the controller in the latter case, but to the dynamical regime of the working medium in the former.

Within this broader landscape, classical ergotropy extraction driving establishes that the maximal-work problem for isolated classical systems can be solved constructively under the ergodic hypothesis, while also providing a framework for more general non-ergodic and topologically constrained cases. The result is a precise classical counterpart to passive-state engineering in quantum thermodynamics and a basis for practical energy-recovery protocols in classical Hamiltonian systems [2508.12797].

Source: https://www.emergentmind.com/topics/classical-ergotropy-extraction-driving