Gibbs-State Limit Cycle in Quantum Otto Machines

• GSLC is defined as the ideal limit cycle where each isochoric stroke instantaneously equilibrates the system to the Gibbs state of the interacting Hamiltonian.
• Internal coupling alters effective energy gaps and induces coherence, broadening the operational regime beyond traditional Otto cycle bounds.
• The global GKSL master equation validates GSLC behavior in the infinite interaction-time limit, providing a benchmark against finite-time NELC and ELC dynamics.

The Gibbs-state limit cycle (GSLC) is the idealized limit cycle of a quantum Otto cycle in which the working medium is assumed to thermalize essentially instantly during each isochoric stroke, so that after each hot or cold contact the state is the corresponding Gibbs state of the interacting system Hamiltonian. In the treatment of the internally coupled Otto cycle studied in "Quantum Thermal Machines Improved by Internal Coupling: From Equilibrium to Non-equilibrium Limit Cycles" (Gao et al., 2 Mar 2026), the GSLC serves as the equilibrium benchmark against which the equilibrating limit cycle (ELC) and non-equilibrating limit cycle (NELC) are compared. Within that framework, internal coupling broadens the operational regime, can turn parameter regions where the uncoupled cycle does not operate as any thermal machine into engine or refrigerator regimes, and can enhance efficiency or coefficient of performance (COP) beyond the standard Otto bounds while remaining below the Carnot limit (Gao et al., 2 Mar 2026).

1. Conceptual definition and placement among limit cycles

In the paper, the GSLC is the limit cycle in which one ignores the influence of the interaction time and assumes that during each isochoric process the system rapidly reaches equilibrium. The end-of-stroke states are therefore Gibbs states of the interacting Hamiltonian. The paper’s main conceptual point is that GSLC is the equilibrium benchmark, ELC is the long-time dynamical realization, and NELC is the finite-time nonequilibrium realization.

The three limit cycles are distinguished entirely by the bath-contact time and by whether equilibration is complete.

Limit cycle	Bath-contact regime	End-of-stroke state
GSLC	Equilibrium assumed immediately	Gibbs state of the interacting Hamiltonian
ELC	Very long time	Steady state of the open dynamics; Gibbs-like equilibrium state of the global master equation
NELC	Finite/short time	Periodic non-equilibrium steady cycle with distinct hot- and cold-stroke states

A central result is that ELC approaches GSLC when the global Gorini--Kossakowski--Sudarshan--Lindblad (GKSL) description is used. This establishes GSLC as the reference equilibrium construction and ELC as its open-system long-time realization, while NELC captures the finite-time regime in which the system does not fully equilibrate.

2. Interacting Hamiltonian and effective level structure

For the uncoupled cycle, the system Hamiltonian during the hot and cold strokes is

$$H_S^\alpha = \omega_\alpha \sigma_S^+ \sigma_S^- = \begin{pmatrix} 0 & 0 \ 0 & \omega_\alpha \end{pmatrix}, \qquad \alpha=h,c.$$

The central modification is the inclusion of an internal coupling $g_\alpha$ between the ground and excited levels,

$$H_S^\alpha = \begin{pmatrix} 0 & g_\alpha \ g_\alpha & \omega_\alpha \end{pmatrix}.$$

This coupling is present during both isochoric strokes, with generally different values $g_h$ and $g_c$. It is not a system-bath coupling; rather, it is an internal coherence-generating coupling in the working medium.

The baths are bosonic,

$$H_B^\alpha = \sum_k \hat{b}_\alpha^\dagger \hat{b}_\alpha,$$

and the system-bath interaction is taken in the standard form

$$H_I^\alpha = \sum_{k,\alpha} V_{k,\alpha}\,\sigma_S^x(\hat{b}_\alpha+\hat{b}_\alpha^\dagger).$$

Diagonalizing $H_S^\alpha$ yields the eigenvalues

$$\epsilon_\alpha^\pm = \frac{1}{2}\left(\omega_\alpha \pm \sqrt{4g_\alpha^2+\omega_\alpha^2}\right),$$

with splitting

$$\tilde{\omega}_\alpha = \epsilon_\alpha^+ - \epsilon_\alpha^- = \sqrt{4g_\alpha^2+\omega_\alpha^2}.$$

The diagonalizing angle is

$g_\alpha$0

and the diagonal basis states $g_\alpha$1 define the Gibbs state. The physical role of $g_\alpha$2 is twofold: it changes the effective energy gaps and creates off-diagonal coherence in the original basis. According to the paper, this is the mechanism behind broader operation regimes, enhanced efficiency and COP, and the possibility of engine or refrigerator operation in regions where the uncoupled Otto cycle would fail.

3. Otto-cycle strokes and Gibbs-state construction

The cycle has four strokes:

1. Hot isochore $g_\alpha$3: system contacts hot bath at $g_\alpha$4, Hamiltonian fixed at $g_\alpha$5.
2. Adiabatic stroke $g_\alpha$6: Hamiltonian changes from $g_\alpha$7 to $g_\alpha$8, state frozen.
3. Cold isochore $g_\alpha$9: system contacts cold bath at $$H_S^\alpha = \begin{pmatrix} 0 & g_\alpha \ g_\alpha & \omega_\alpha \end{pmatrix}.$$0, Hamiltonian fixed at $$H_S^\alpha = \begin{pmatrix} 0 & g_\alpha \ g_\alpha & \omega_\alpha \end{pmatrix}.$$1.
4. Adiabatic stroke $$H_S^\alpha = \begin{pmatrix} 0 & g_\alpha \ g_\alpha & \omega_\alpha \end{pmatrix}.$$2: Hamiltonian changes back from $$H_S^\alpha = \begin{pmatrix} 0 & g_\alpha \ g_\alpha & \omega_\alpha \end{pmatrix}.$$3 to $$H_S^\alpha = \begin{pmatrix} 0 & g_\alpha \ g_\alpha & \omega_\alpha \end{pmatrix}.$$4, state frozen.

For GSLC, the defining assumption is that during each isochoric stroke the system quickly equilibrates to the Gibbs state of the interacting Hamiltonian. In the eigenbasis, the state is

$$H_S^\alpha = \begin{pmatrix} 0 & g_\alpha \ g_\alpha & \omega_\alpha \end{pmatrix}.$$5

In the original basis, the Gibbs state contains coherence:

$$H_S^\alpha = \begin{pmatrix} 0 & g_\alpha \ g_\alpha & \omega_\alpha \end{pmatrix}.$$6

with

$$H_S^\alpha = \begin{pmatrix} 0 & g_\alpha \ g_\alpha & \omega_\alpha \end{pmatrix}.$$7

The off-diagonal term is the explicit manifestation of interaction-induced coherence in the original basis. This is why the coupled cycle is not governed solely by the bare frequency ratio $$H_S^\alpha = \begin{pmatrix} 0 & g_\alpha \ g_\alpha & \omega_\alpha \end{pmatrix}.$$8.

4. Operational regime, work and heat, and performance criteria

Without coupling, the standard Otto expressions are

$$H_S^\alpha = \begin{pmatrix} 0 & g_\alpha \ g_\alpha & \omega_\alpha \end{pmatrix}.$$9

with engine operation when

$g_h$0

and refrigerator operation when

$g_h$1

A key special case emphasized in the paper is

$g_h$2

For the uncoupled cycle, this gives zero work and no machine operation. With internal coupling, however, the cycle can still function (Gao et al., 2 Mar 2026).

In the GSLC, the operational regime depends on the ratios $g_h$3 and $g_h$4, and on

$g_h$5

For the engine, a higher efficiency than Otto occurs when

$g_h$6

together with $g_h$7, $g_h$8, $g_h$9.

For the refrigerator, an improved COP occurs when

$g_c$0

together with $g_c$1, $g_c$2, $g_c$3.

The paper identifies two special cases of particular importance. If $g_c$4, the uncoupled cycle cannot work, but the coupled one can, controlled entirely by $g_c$5 and $g_c$6. If $g_c$7, the uncoupled cycle has $g_c$8, but coupling can produce nonzero work. In that sense, the internal coupling acts as a thermodynamic resource.

With coupling, the heats are written as

$g_c$9

$$H_B^\alpha = \sum_k \hat{b}_\alpha^\dagger \hat{b}_\alpha,$$0

and the adiabatic-stroke work contributions are

$$H_B^\alpha = \sum_k \hat{b}_\alpha^\dagger \hat{b}_\alpha,$$1

$$H_B^\alpha = \sum_k \hat{b}_\alpha^\dagger \hat{b}_\alpha,$$2

The total work is

$$H_B^\alpha = \sum_k \hat{b}_\alpha^\dagger \hat{b}_\alpha,$$3

For engine operation,

$$H_B^\alpha = \sum_k \hat{b}_\alpha^\dagger \hat{b}_\alpha,$$4

where $$H_B^\alpha = \sum_k \hat{b}_\alpha^\dagger \hat{b}_\alpha,$$5 is a coupling-dependent factor built from $$H_B^\alpha = \sum_k \hat{b}_\alpha^\dagger \hat{b}_\alpha,$$6 and the coherence function $$H_B^\alpha = \sum_k \hat{b}_\alpha^\dagger \hat{b}_\alpha,$$7.

For refrigerator operation,

$$H_B^\alpha = \sum_k \hat{b}_\alpha^\dagger \hat{b}_\alpha,$$8

where $$H_B^\alpha = \sum_k \hat{b}_\alpha^\dagger \hat{b}_\alpha,$$9 is another coupling-dependent factor.

The standard Otto bounds are recovered when $$H_I^\alpha = \sum_{k,\alpha} V_{k,\alpha}\,\sigma_S^x(\hat{b}_\alpha+\hat{b}_\alpha^\dagger).$$0, because then $$H_I^\alpha = \sum_{k,\alpha} V_{k,\alpha}\,\sigma_S^x(\hat{b}_\alpha+\hat{b}_\alpha^\dagger).$$1. With coupling, $$H_I^\alpha = \sum_{k,\alpha} V_{k,\alpha}\,\sigma_S^x(\hat{b}_\alpha+\hat{b}_\alpha^\dagger).$$2 or $$H_I^\alpha = \sum_{k,\alpha} V_{k,\alpha}\,\sigma_S^x(\hat{b}_\alpha+\hat{b}_\alpha^\dagger).$$3 may become larger or smaller than their Otto counterparts. The paper’s explanation is that internal coupling changes the effective energy gaps, introduces coherence in the Gibbs state, and modifies how much heat is drawn from each bath relative to work output. Despite this enhancement, the performance remains below the Carnot limit, because the cycle remains a thermal machine obeying the first and second laws.

5. GSLC, ELC, NELC, and the global GKSL description

A major result is the comparison between the ideal GSLC and the open-system dynamical limits (Gao et al., 2 Mar 2026). The paper stresses that the local master equation neglects internal coupling in the dissipator and gives a steady state that is not the correct Gibbs state for the interacting Hamiltonian. In this setting, it is therefore thermodynamically inconsistent.

By contrast, the global master equation diagonalizes the interacting system first and yields the equilibrium state

$$H_I^\alpha = \sum_{k,\alpha} V_{k,\alpha}\,\sigma_S^x(\hat{b}_\alpha+\hat{b}_\alpha^\dagger).$$4

which is Gibbs-like in the eigenbasis.

In the infinite interaction-time limit, the system equilibrates during each stroke and reaches the ELC. The paper shows that the ELC has the same properties as the GSLC, which is taken as validation of the global GKSL approach.

For finite short interaction times, the system does not fully equilibrate and instead converges to a periodic non-equilibrium cycle. The authors prove the existence of this NELC using a primitive CPTP map and the quantum Perron--Frobenius theorem. Its hot- and cold-stroke states satisfy

$$H_I^\alpha = \sum_{k,\alpha} V_{k,\alpha}\,\sigma_S^x(\hat{b}_\alpha+\hat{b}_\alpha^\dagger).$$5

and the cycle is obtained stroboscopically from repeated application of the one-cycle map.

The comparison is therefore structurally sharp: GSLC and ELC are equilibrium-like and exhibit higher efficiency or COP, whereas NELC is nonequilibrium and has reduced efficiency or COP. The global GKSL master equation reproduces the ELC and agrees with GSLC in the long-time limit.

6. Finite-time operation and the power--efficiency trade-off

The finite-time analysis isolates the role of interaction time in the transition from NELC to ELC. For short interaction time, energy exchange is limited, so efficiency and COP are low, but power can be high because the cycle time is short. For long interaction time, the system thermalizes more completely, so efficiency and COP increase, but power decreases because the cycle duration becomes long. In the infinite interaction-time limit, the NELC converges to ELC; power goes to zero, while efficiency and COP approach their equilibrium values.

For engine operation, the power is defined as

$$H_I^\alpha = \sum_{k,\alpha} V_{k,\alpha}\,\sigma_S^x(\hat{b}_\alpha+\hat{b}_\alpha^\dagger).$$6

The paper explicitly observes that increasing $$H_I^\alpha = \sum_{k,\alpha} V_{k,\alpha}\,\sigma_S^x(\hat{b}_\alpha+\hat{b}_\alpha^\dagger).$$7 makes NELC approach ELC, that efficiency rises toward its equilibrium maximum, and that power falls toward zero. This is the standard power--efficiency trade-off: approaching reversible equilibrium improves performance metrics while reducing output rate.

Taken together, these results position the GSLC as the equilibrium reference model for the internally coupled Otto machine. In that role, it clarifies how internal coupling creates coherence, modifies effective energy gaps, enlarges the operational regime, and enables performance beyond the standard Otto bounds without surpassing Carnot, while the ELC and NELC supply the corresponding long-time and finite-time dynamical realizations.