---
title: Gibbs-Preserving Operations in Quantum Thermodynamics
url: https://www.emergentmind.com/topics/gibbs-preserving-operations
type: topic
---

# Gibbs-Preserving Operations in Quantum Thermodynamics

Searching arXiv for recent papers on Gibbs-preserving operations and related frameworks.
Gibbs-preserving operations, also called Gibbs-preserving maps, are completely positive trace-preserving maps that leave a thermal Gibbs state invariant. For a finite-dimensional system with Hamiltonian $H$ at inverse temperature $\beta$, the Gibbs state is
\[
\gamma_\beta=\frac{e^{-\beta H}}{Z},\qquad Z=\operatorname{Tr}[e^{-\beta H}],
\]
and a channel $\Phi$ is Gibbs-preserving when $\Phi(\gamma_\beta)=\gamma_\beta$. In quantum thermodynamics, this condition abstracts the equilibrium-preservation aspect of the second law and defines a broad class of free processes. The subject has developed into a resource-theoretic framework connecting thermal operations, thermo-majorization, $D$-majorization, covariant dynamics, and the operational role of coherence. A central structural fact is that Gibbs-preserving operations coincide with thermal operations only in the energy-diagonal regime; in the genuinely quantum regime they are strictly more powerful [1406.3618].

## 1. Formal setting and basic definitions

The standard setting fixes a finite-dimensional quantum system and a temperature. In the single-system formulation, a Gibbs-preserving operation is a CPTP map $\Phi$ satisfying $\Phi(D)=D$ for a full-rank Gibbs state
\[
D:=\frac{e^{-\beta H}}{\operatorname{Tr}[e^{-\beta H}]}.
\]
In the more general input-output formulation, a channel $\Lambda:S\to S'$ is Gibbs-preserving if $\Lambda(\tau_S)=\tau_{S'}$, where $\tau_X=e^{-\beta H_X}/\operatorname{Tr}(e^{-\beta H_X})$ denotes the Gibbs state of system $X$ [2003.04164][2404.03479].

At infinite temperature, $\beta\to 0$, the Gibbs state becomes maximally mixed, and Gibbs preservation reduces to unitality. In this sense, Gibbs-preserving operations generalize unital channels from ordinary majorization theory to a thermodynamic setting with a nontrivial fixed point $D$ [2003.04164].

In the quasi-classical or energy-diagonal regime, a Gibbs-preserving operation is represented by a Gibbs-stochastic matrix $M$ acting on probability vectors, with nonnegative entries, column sums equal to one, and
\[
Mg=g,
\]
where $g_i=e^{-\beta E_i}/Z$ is the Gibbs distribution. In the unnormalized convention $d_i:=e^{-\beta E_i}$, the condition becomes $Ad=d$ [2212.04305].

## 2. Relation to thermal operations and covariance

Thermal operations are defined microscopically by coupling the system to a thermal environment, applying an energy-conserving unitary, and tracing out an environment subsystem:
\[
\Lambda(\rho)=\mathrm{Tr}_{E'}\bigl[U(\rho\otimes \tau_E)U^\dagger\bigr],\qquad [U,H_{\rm tot}]=0.
\]
They are automatically Gibbs-preserving, so
\[
O_{\rm TO}\subseteq O_{\rm GP}.
\]
A more restrictive intermediate class is given by covariant Gibbs-preserving operations, also called enhanced thermal operations, which satisfy both Gibbs preservation and time-translation covariance [1406.3618][2406.06234].

| Class | Defining condition | Inclusion |
|---|---|---|
| GPO/GPM | $\Phi(\gamma_\beta)=\gamma_\beta$ | Broadest class here |
| CGPO / EnTO | Gibbs-preserving and covariant | TO $\subset$ CGPO $\subseteq$ GPO |
| TO | Thermal ancilla + energy-conserving unitary + discard | Physically motivated subclass |

The crucial dynamical distinction is covariance. Thermal operations satisfy
\[
\Phi\!\big(e^{-iHt}\rho\,e^{iHt}\big)=e^{-iHt}\Phi(\rho)e^{iHt}
\]
for all $t\in\mathbb{R}$. Consequently, for nondegenerate spectra, they cannot create coherence between distinct energy eigenspaces from an energy-diagonal input. An energy eigenstate is invariant under the system’s free evolution, and covariance forces the output to be time-invariant as well; for nondegenerate Hamiltonians, that means energy-diagonal. Gibbs-preserving operations impose no such covariance requirement: they need only fix the Gibbs state [1406.3618].

This separation is not merely formal. The 2014 note established that thermal operations are strictly contained in Gibbs-preserving maps in the quantum regime, while the 2024 coherence-cost analysis showed that some Gibbs-preserving operations remain unattainable even if thermal operations are assisted by any finite amount of quantum coherence [1406.3618][2404.03479].

## 3. Classical regime: $D$-majorization and thermo-majorization

For matrices or states $X,Y$, the preorder induced by Gibbs-preserving CPTP maps is
\[
Y \preceq_D X \quad \Longleftrightarrow \quad \exists\,\Phi\ \text{CPTP such that}\ \Phi(D)=D,\ \Phi(X)=Y.
\]
In the diagonal case $D=\operatorname{diag}(d)$, $X=\operatorname{diag}(x)$, $Y=\operatorname{diag}(y)$, this reduces to vector $d$-majorization, written $y\prec_d x$ [2003.04164].

For commuting states, this is the standard thermo-majorization criterion. One characterization uses $\beta$-ordered Lorenz curves:
\[
\sigma \preceq_D \rho \;\; \text{iff} \;\; L_\sigma^{(g)}(x) \le L_\rho^{(g)}(x) \;\; \forall x,
\]
where the Gibbs weights $g_i\propto e^{-\beta E_i}$ determine the ordering. In this abelian regime, Gibbs-preserving operations and thermal operations coincide at the level of state transitions, and thermo-majorization completely characterizes reachability [1406.3618][2003.04164].

The order-theoretic structure is richer than the diagonal criterion alone suggests. The relation $\preceq_D$ is a preorder—reflexive and transitive—but not a partial order in general. For trace-constrained Hermitian matrices, $D$ is the unique minimal element, while a rank-one projector supported on the minimal coordinate of $d$ is maximal in the positive cone under suitable trace normalization. The reachable-set operator
\[
M_D(S):=\bigcup_{\rho\in S}\{X:\;X\preceq_D \rho\}
\]
has strong regularity properties: $M_D(X)$ is convex; $M_D(P)$ is compact for compact $P$; on sets of states it is star-shaped with respect to the Gibbs state $D/\operatorname{Tr}(D)$; and on compact subsets it is non-expansive in the Hausdorff metric [2003.04164].

In the quasi-classical theory, these reachable sets become thermomajorization polytopes. For a Gibbs vector $d$ and initial vector $y$, the polytope
\[
M_d(y)=\{Ay: A\in s_d(n)\}
\]
admits both a halfspace description and an extreme-point description. The geometry depends sensitively on the Gibbs spectrum. “Stable” Gibbs states are precisely those for which global cyclic state transfers are impossible in the quasi-classical regime, while “well-structured” Gibbs states allow degeneracies of extreme points to witness equilibrium subspaces [2212.04305].

## 4. Quantum regime: strict separation and explicit coherence-generating maps

The decisive quantum feature is energy coherence, meaning off-diagonal matrix elements between distinct energy eigenspaces. In this regime,
\[
\mathrm{TO}\subsetneq \mathrm{GPM}.
\]
The strict inclusion is witnessed by explicit Gibbs-preserving channels that transform an energy eigenstate into an arbitrary target state, including coherent superpositions, while preserving the Gibbs state [1406.3618].

For a qubit with Gibbs state
\[
\gamma=p_0|0\rangle\langle 0|+p_1|1\rangle\langle 1|,
\]
fix any target state $\rho$, possibly coherent in the energy basis. Define
\[
\Phi(X)=\langle 0|X|0\rangle\,\sigma+\langle 1|X|1\rangle\,\rho,
\]
with
\[
\sigma=p_0^{-1}(\gamma-p_1\rho).
\]
This channel is a measurement in the energy basis followed by conditional state preparation. It is CPTP, satisfies $\Phi(\gamma)=\gamma$, and obeys $\Phi(|1\rangle\langle 1|)=\rho$. If $\rho$ contains coherence across distinct energies, the transition is impossible under thermal operations by covariance, but possible under Gibbs preservation alone [1406.3618].

The construction generalizes to $n$-level systems by taking a highest-energy level $|n\rangle$ and setting
\[
\Phi(X)=\operatorname{Tr}\!\big[(\mathbb{I}-|n\rangle\langle n|)X\big]\sigma+\operatorname{Tr}\!\big[|n\rangle\langle n|X\big]\rho,
\qquad
\sigma=\frac{\gamma-p_n\rho}{1-p_n}.
\]
Again, $\Phi$ is CPTP, maps $|n\rangle\langle n|$ to $\rho$, and preserves $\gamma$ [1406.3618].

A common misconception is that classical equivalence implies full quantum equivalence. The literature distinguishes these statements sharply. The classical result concerns state transitions between block-diagonal states; it does not imply that the channel sets themselves coincide. Even when diagonal-state convertibility agrees, the quantum channel sets differ because thermal operations retain symmetry restrictions that Gibbs-preserving maps do not [1406.3618].

## 5. Monotones and free-energy formulations

The basic monotone under any Gibbs-preserving operation is the relative entropy to the Gibbs state,
\[
D(\rho\Vert\tau)=\operatorname{Tr}[\rho(\log\rho-\log\tau)].
\]
If $\mathcal{E}(\tau)=\tau$, then the data-processing inequality gives
\[
D(\rho\Vert\tau)\ge D(\mathcal{E}(\rho)\Vert\tau).
\]
Equivalently, the nonequilibrium free energy
\[
F(\rho):=\operatorname{Tr}[\rho H]-kT\,S(\rho)
\]
satisfies
\[
F(\rho)-F(\tau)=kT\,D(\rho\Vert\tau),
\]
so $F(\rho)-F(\tau)$ cannot increase under Gibbs-preserving maps. The same logic applies to Rényi relative entropies and the associated generalized free energies [1406.3618].

These monotones constrain both Gibbs-preserving operations and thermal operations, but they do not capture the full distinction between them. The missing ingredient is symmetry: thermal operations satisfy time-translation covariance, whereas Gibbs preservation alone does not encode the asymmetry constraints associated with coherence [1406.3618].

Recent work on covariant Gibbs-preserving operations sharpened the free-energy picture in a different direction. For finite-dimensional systems with commensurate energy gaps, if the initial state is coherent, distillable, and has shortest period $2\pi/\Delta$, then correlated-catalytic state convertibility under covariant Gibbs-preserving operations is fully characterized by the single monotone
\[
F(\rho):=D(\rho\Vert\gamma_\beta).
\]
Under these assumptions, $\rho\to \sigma$ is possible arbitrarily well if and only if
\[
D(\rho\Vert\gamma_\beta)\ge D(\sigma\Vert\gamma_\beta),
\]
and if the inequality is strict and $\sigma$ is full-rank, there exists an exact correlated-catalytic implementation. In this correlated-catalytic regime, imposing covariance does not change convertibility for coherent initial states [2406.06234].

This result does not erase the earlier separation between thermal operations and Gibbs-preserving operations. Rather, it identifies a catalytic regime in which the additional covariance constraint ceases to affect state convertibility, provided the initial coherence already supplies the relevant timing resource [2406.06234].

## 6. Coherence cost and operational significance

The operational status of Gibbs-preserving operations is complicated by implementation cost. Coherence in this setting means time-translation asymmetry, and the 2024 analysis quantifies it by the quantum Fisher information
\[
I(\rho)=2\sum_{i,j}\frac{(\lambda_i-\lambda_j)^2}{\lambda_i+\lambda_j}\; \big|\langle e_i|H|e_j\rangle\big|^2
\]
for $\rho=\sum_i\lambda_i|e_i\rangle\langle e_i|$. The coherence cost of implementing a channel via thermal operations assisted by an ancillary state is defined as the minimum ancilla QFI required for exact or approximate realization [2404.03479].

For pairwise reversible Gibbs-preserving operations, the key quantity is
\[
C(\Lambda,P)\coloneqq\bigl\|\sqrt{\rho_1}\bigl(H_S-\Lambda^\dagger(H_{S'})\bigr)\sqrt{\rho_2}\bigr\|_2,
\]
where $P=\{\rho_1,\rho_2\}$ is a reversible pair. The main lower bound is
\[
\sqrt{C_c^\epsilon(\Lambda)}\ \ge\ \frac{C(\Lambda,P)}{\epsilon}\ -\ \Delta(H_S)\ -\ 3\,\Delta(H_{S'}),
\]
so if $C(\Lambda,P)>0$, the required coherence diverges as $\epsilon\to 0$. Exact implementation then requires an infinite amount of coherence [2404.03479].

A concrete qubit example is the coherence-detecting channel
\[
\Lambda(\rho)=\langle +|\rho|+\rangle\,|0\rangle\langle 0|+\langle -|\rho|-\rangle\,|1\rangle\langle 1|,
\]
with $|\pm\rangle=(|0\rangle\pm|1\rangle)/\sqrt{2}$. For a qubit input Hamiltonian $H_S=|1\rangle\langle 1|$ and trivial output Hamiltonian, this channel is Gibbs-preserving, the reversible pair is $\{|+\rangle,|-\rangle\}$, and
\[
C(\Lambda,P)=\frac{1}{2}.
\]
Accordingly, no finite amount of coherence suffices for exact implementation by thermal operations. The same paper shows that there are uncountably many Gibbs-preserving operations with unbounded coherence cost, obtained from continuous families of coherent-basis measure-and-prepare channels [2404.03479].

The contrast with earlier examples is important. Coherence-creating Gibbs-preserving measure-and-prepare maps, such as the qubit channel that sends an excited energy eigenstate to a coherent target state, can have finite coherence cost, whereas coherence-detecting Gibbs-preserving maps can require infinite coherence. This suggests that Gibbs preservation alone is too weak to guarantee operational affordability [1406.3618][2404.03479].

## 7. Broader hierarchies, limitations, and open problems

A broader 2025 hierarchy of thermodynamically consistent quantum operations places Gibbs-preservation-like fixed-point conditions in a larger framework. In that hierarchy, Class III channels—intended to encode consistency with both the second and third laws—must be rank non-decreasing and must not perturb a strictly positive state. More precisely, any Class III system channel has at least one faithful fixed point $\rho_*>0$. This generalizes the fixed-point aspect of Gibbs-preserving operations, but does not identify the fixed point with a preassigned Gibbs state unless additional structure, such as energy conservation and a thermal ancilla, is imposed. Thermal operations therefore lie inside Class III, but Class III is in general broader than Gibbs-preserving operations for a fixed $(H,\beta)$ [2505.23360].

Several technical limitations recur across the literature. Most results are formulated for finite-dimensional systems and finite inverse temperature. The sharp classical equivalence between thermal operations and Gibbs-preserving operations concerns only block-diagonal state transitions. The qubit case admits special simplifications: $D$-majorization is completely characterized by the family of trace-norm inequalities
\[
\|A-tD\|_1\le \|B-tD\|_1\qquad \forall t\in\mathbb{R},
\]
but this characterization fails in dimensions $d\ge 3$, where the transpose-based counterexample shows that simple trace-norm families are not complete [2003.04164].

Open problems identified in the cited works include a full characterization of $\preceq_D$ beyond qubits, efficient algorithmic tests for noncommuting states, a complete family of monotones for general Gibbs-preserving operations in higher dimensions, tighter relations between thermal operations, enhanced thermal operations, and Gibbs-preserving operations in catalytic regimes, and the extension of coherence-cost bounds to infinite-dimensional systems with unbounded Hamiltonians [2003.04164][2404.03479][2406.06234].

Taken together, these developments delineate a precise picture. Gibbs-preserving operations provide a mathematically natural and thermodynamically meaningful fixed-point condition. In the quasi-classical regime they reproduce thermo-majorization and the reachability theory of Gibbs-stochastic maps. In the quantum regime they strictly exceed thermal operations because Gibbs preservation does not enforce time-translation covariance. Recent work further shows that this broader power can conceal unbounded coherence requirements, so the choice between Gibbs-preserving operations, covariant Gibbs-preserving operations, and thermal operations is not merely formal; it determines which thermodynamic resources are treated as genuinely free [1406.3618][2404.03479][2406.06234].

Source: https://www.emergentmind.com/topics/gibbs-preserving-operations