Equilibrating Limit Cycle in Quantum Systems

• Equilibrating Limit Cycle (ELC) is a construction where open or switched dynamics fully relax into stationary periodic cycles, providing a precise equilibrium benchmark in quantum thermodynamics.
• It employs the global GKSL master equation to derive dressed thermal fixed points, capturing the effects of internal coupling and dissipative processes.
• ELC serves as a reference to compare finite-time nonequilibrium cycles (NELC) against ideal Gibbs-state limit cycles, thus informing performance trade-offs in quantum Otto engines.

Searching arXiv for the cited primary paper and closely related limit-cycle literature to ground the article. An equilibrating limit cycle (ELC) is a limit-cycle construction in which the relevant open or switched dynamics are allowed to relax fully on the dissipative or constrained portions of the evolution, so that each cycle is built from the exact stationary state or stationary motion selected by the underlying dynamics rather than from an imposed finite-time trajectory. In the quantum thermal-machine setting, the term has a precise microscopic meaning: for a quantum Otto cycle with internal coupling, the ELC is obtained by taking infinitely long hot and cold isochores under the global Gorini–Kossakowski–Sudarshan–Lindblad (GKSL) master equation, so that the system reaches the true GKSL fixed point on each isochore (Gao et al., 2 Mar 2026). In other arXiv literature, the same phrase is either used directly or is a natural interpretive label for long-time periodic attractors generated by drive–dissipation balance, hysteretic switching, or noisy relaxation toward a stable cycle (Chan et al., 2015, Makarenkov, 2017, Louca, 2015, Ríos-Monje et al., 2024, Makarenkov et al., 4 Nov 2025).

1. Quantum-thermodynamic definition

In the quantum Otto model of a coupled two-level working medium, the system Hamiltonian on stroke $\alpha \in \{h,c\}$ is

$$H_S^{\alpha}(\lambda_\alpha) = \begin{pmatrix} 0 & g_\alpha \ g_\alpha & \omega_\alpha \end{pmatrix},$$

where $\omega_\alpha$ is the bare energy gap and $g_\alpha$ is the internal coupling strength operative on that stroke. The Otto cycle consists of a hot isochore, a compression stroke, a cold isochore, and an expansion stroke. During the isochores the Hamiltonian is fixed and the system exchanges heat with a hot bath at $T_h$ or a cold bath at $T_c$; during the adiabatic strokes the parameters switch between $(\omega_h,g_h)$ and $(\omega_c,g_c)$ while the state is kept unchanged (Gao et al., 2 Mar 2026).

Because $g_\alpha \neq 0$, the internal coupling mixes the bare ground and excited states. Diagonalization gives dressed eigenvalues

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

and dressed transition gap

$$H_S^{\alpha}(\lambda_\alpha) = \begin{pmatrix} 0 & g_\alpha \ g_\alpha & \omega_\alpha \end{pmatrix},$$0

The associated mixing angle is

$$H_S^{\alpha}(\lambda_\alpha) = \begin{pmatrix} 0 & g_\alpha \ g_\alpha & \omega_\alpha \end{pmatrix},$$1

These dressed quantities are the basic objects entering the open-system dynamics and the thermodynamic performance formulas (Gao et al., 2 Mar 2026).

Within this framework, the ELC is defined by modeling each isochore with the global GKSL master equation and then taking the infinite interaction-time limit. The state on each isochore relaxes to the true GKSL fixed point, so the cycle becomes a periodic steady operation built from those stationary states. Under the assumptions used in the model—Ohmic bath, secular approximation, and weak coupling—the fixed point coincides with the thermal state in the diagonalized basis, and the ELC reproduces the Gibbs-state limit cycle (GSLC) after transforming back to the original basis (Gao et al., 2 Mar 2026).

2. Microscopic dynamics and the stationary state

The open dynamics on each isochore are governed in the diagonalized basis by the global GKSL equation

$$H_S^{\alpha}(\lambda_\alpha) = \begin{pmatrix} 0 & g_\alpha \ g_\alpha & \omega_\alpha \end{pmatrix},$$2

For the interaction operator obtained from $$H_S^{\alpha}(\lambda_\alpha) = \begin{pmatrix} 0 & g_\alpha \ g_\alpha & \omega_\alpha \end{pmatrix},$$3, the only dissipators that survive under the Ohmic assumption $$H_S^{\alpha}(\lambda_\alpha) = \begin{pmatrix} 0 & g_\alpha \ g_\alpha & \omega_\alpha \end{pmatrix},$$4 are those associated with transitions between the dressed eigenstates. The jump operators are

$$H_S^{\alpha}(\lambda_\alpha) = \begin{pmatrix} 0 & g_\alpha \ g_\alpha & \omega_\alpha \end{pmatrix},$$5

with effective dressed rates

$$H_S^{\alpha}(\lambda_\alpha) = \begin{pmatrix} 0 & g_\alpha \ g_\alpha & \omega_\alpha \end{pmatrix},$$6

where the rates satisfy detailed balance at temperature $$H_S^{\alpha}(\lambda_\alpha) = \begin{pmatrix} 0 & g_\alpha \ g_\alpha & \omega_\alpha \end{pmatrix},$$7 (Gao et al., 2 Mar 2026).

In the ELC limit, the state is the zero-eigenvalue eigenmatrix of the global Liouvillian on each isochore:

$$H_S^{\alpha}(\lambda_\alpha) = \begin{pmatrix} 0 & g_\alpha \ g_\alpha & \omega_\alpha \end{pmatrix},$$8

Transforming back gives

$$H_S^{\alpha}(\lambda_\alpha) = \begin{pmatrix} 0 & g_\alpha \ g_\alpha & \omega_\alpha \end{pmatrix},$$9

In the present model this is exactly the dressed thermal fixed point, but the same paper notes that with structured or frequency-filtered baths the global GKSL steady state may instead be a generalized thermal state rather than exactly $\omega_\alpha$0 (Gao et al., 2 Mar 2026).

A central point is that the ELC is not defined by assuming a Gibbs state a priori. It is derived as the equilibrium limit of the microscopic global GKSL dynamics. This distinction matters because the paper contrasts the global construction with a local GKSL description built from bare $\omega_\alpha$1 dissipators; the local approach neglects internal coupling in the dissipative part and yields a stationary state that is not the Gibbs state of $\omega_\alpha$2, whereas the global approach correctly incorporates the coupling and recovers the dressed thermal fixed point (Gao et al., 2 Mar 2026).

3. Thermodynamic observables and the role of internal coupling

With the convention that heat absorbed by the system is positive and work is positive when transferred into the system during the unitary strokes, the heat and work over the Otto cycle are

$\omega_\alpha$3

$\omega_\alpha$4

with the first law

$\omega_\alpha$5

Engine efficiency and refrigerator coefficient of performance are

$\omega_\alpha$6

with $\omega_\alpha$7 for a refrigerator (Gao et al., 2 Mar 2026).

For GSLC and ELC, the heat and work admit closed forms in terms of

$\omega_\alpha$8

The uncoupled Otto bounds are

$\omega_\alpha$9

while the Carnot bounds are

$g_\alpha$0

The paper states that internal coupling shifts and enlarges the operational regime and enhances the efficiency and the coefficient of performance, allowing performance to exceed the standard Otto bounds while remaining below the Carnot limit (Gao et al., 2 Mar 2026).

The physical mechanism is described there as spectral engineering. Internal coupling modifies the dressed transition gaps $g_\alpha$1, mixes the eigenstates through $g_\alpha$2, and rescales the dissipative rates by $g_\alpha$3. In addition, equilibrium in the original basis carries off-diagonal coherence. For GSLC,

$g_\alpha$4

Accordingly, when heat and work are expressed in the original basis, coherence contributes explicitly to the thermodynamic balance (Gao et al., 2 Mar 2026).

Operationally, engine action requires $g_\alpha$5, $g_\alpha$6, $g_\alpha$7, while refrigeration requires $g_\alpha$8, $g_\alpha$9, $T_h$0. The paper shows that these conditions can be satisfied in parameter regimes where the uncoupled Otto cycle fails. In the special case $T_h$1 with $T_h$2, there is nonzero work and the performance depends purely on coupling:

$T_h$3

The engine operates for $T_h$4, and the refrigerator for $T_h$5 (Gao et al., 2 Mar 2026).

4. ELC, GSLC, and NELC

The relevant paper distinguishes three limit-cycle notions according to how the isochores are treated (Gao et al., 2 Mar 2026).

Concept	Isochoric treatment	Steady object
GSLC	Rapid or ideal equilibration to Gibbs state	Gibbs state of $T_h$6
ELC	Global GKSL dynamics with infinite interaction time	Fixed point of each isochore’s GKSL semigroup
NELC	Finite interaction times on isochores	Fixed point of the one-cycle CPTP map

In this classification, GSLC is a thermodynamic idealization, ELC is the exact equilibrium limit of the microscopic open-system dynamics, and NELC is the finite-time nonequilibrium periodic state. The paper validates the global GKSL approach by showing that ELC reproduces GSLC heat flows, work, and performance in the infinite-time limit of the isochores (Gao et al., 2 Mar 2026).

For NELC, the one-cycle stroboscopic map is primitive CPTP and admits a unique periodic limit cycle by the quantum Perron–Frobenius theorem. The power is

$T_h$7

Numerically, as $T_h$8, heat flows and work vanish, so $T_h$9 and $T_c$0; as $T_c$1 increases, efficiency and coefficient of performance increase toward the ELC values, while power eventually decreases and vanishes as $T_c$2. The paper identifies this behavior as a power–efficiency trade-off: ELC gives higher efficiency and COP, whereas NELC gives finite power at finite interaction time (Gao et al., 2 Mar 2026).

The representative parameter set $T_c$3, $T_c$4, $T_c$5, $T_c$6, $T_c$7, $T_c$8 is used to show convergence of $T_c$9, $(\omega_h,g_h)$0, and $(\omega_h,g_h)$1 from NELC to ELC/GSLC. Phase diagrams further show that NELC operational regions are narrower than ELC regions and expand toward them as $(\omega_h,g_h)$2 increases (Gao et al., 2 Mar 2026).

5. Assumptions, caveats, and validity domain

The ELC construction in the coupled quantum Otto problem is conditioned on weak system–bath coupling, Markovian dynamics, a secular approximation in the global eigenbasis of $(\omega_h,g_h)$3, and an Ohmic spectral density $(\omega_h,g_h)$4 that removes the zero-frequency dephasing channel. Within those assumptions, the global master equation is thermodynamically consistent in the presence of internal coupling and its stationary state coincides with the dressed thermal state (Gao et al., 2 Mar 2026).

The same source also states several limitations. Strong internal coupling or strong system–bath coupling can invalidate weak-coupling Markovian GKSL and the secular approximation. Structured baths, non-Ohmic spectra, or frequency filtering can produce steady states that differ from a simple Gibbs form, so ELC then reflects the true GKSL fixed point rather than $(\omega_h,g_h)$5. Finite-time non-secular corrections can modify NELC behavior and quantitatively affect the power–efficiency trade-off (Gao et al., 2 Mar 2026).

These caveats are important for the status of the term itself. In the narrow quantum-thermodynamic sense, ELC is a sharply defined infinite-time construction. This suggests that the phrase should not be treated as a generic synonym for any periodic attractor; rather, its exact meaning depends on which dynamics are being equilibrated and what object is taken as the long-time limit.

6. Broader uses of the ELC idea

Beyond the quantum Otto setting, several arXiv papers support a broader but non-uniform use of the idea. In driven-dissipative spin-$(\omega_h,g_h)$6 systems, the long-time attractor can be a time-periodic nonequilibrium steady state of a Lindblad master equation in which collective magnetizations oscillate periodically in time. That paper does not use the phrase “Equilibrating Limit Cycle,” but it explicitly describes the interpretation as natural: the system relaxes toward a stable periodic attractor under drive and dissipation, with a spatio-temporal Goldstone mode and quasi-long-range limit-cycle ordering in $(\omega_h,g_h)$7 (Chan et al., 2015).

In switched planar systems, a regular asymptotically stable sliding equilibrium of the Filippov equation can yield an orbitally stable limit cycle when a small hysteresis band is introduced. The paper gives the first-order period

$(\omega_h,g_h)$8

and interprets the resulting orbit as an ELC that replaces the regulating sliding equilibrium when sliding is eliminated by hysteresis (Makarenkov, 2017). A related 2025 glacial-cycles analysis studies an attracting limit cycle born from a degenerate fold-fold singularity at a switched equilibrium; that paper also says it does not use the term “Equilibrating Limit Cycle,” but that the construction matches the ELC concept (Makarenkov et al., 4 Nov 2025).

In stochastic dynamical systems with a stable deterministic cycle perturbed by additive Gaussian white noise, one paper describes an ELC as a noisy limit cycle whose transverse deviations equilibrate to a stationary Gaussian law while the phase diffuses in a Brownian-like fashion. In the Hopf normal form, the transverse coordinate obeys an Ornstein–Uhlenbeck process and the phase diffusion coefficient is $(\omega_h,g_h)$9, yielding stationary autocovariances and broadened spectral peaks (Louca, 2015). In optimal-control work on the van der Pol and more general Liénard oscillators, ELC denotes finite-time synchronization to a stable limit cycle by a time-dependent force, with an exact speed-limit inequality relating connection time and non-conservative work (Ríos-Monje et al., 2024).

Taken together, these papers indicate that ELC is best understood as a family resemblance rather than a single universal formalism. Across quantum thermodynamics, driven-dissipative many-body physics, Filippov systems, noisy oscillators, and controlled nonlinear oscillators, the common element is relaxation toward a stable periodic attractor or its exact long-time surrogate; the precise mathematical object can be a GKSL fixed point on each stroke, a periodic Lindblad attractor, an orbitally stable hysteresis-induced cycle, a noisy cycle with stationary transverse statistics, or a finite-time controlled arrival on a stable limit cycle (Gao et al., 2 Mar 2026, Chan et al., 2015, Makarenkov, 2017, Louca, 2015, Ríos-Monje et al., 2024).

7. Significance

Within quantum thermal machines, the ELC provides an exact equilibrium benchmark for internally coupled working media. It isolates the infinite-time behavior of the global GKSL dynamics, establishes when GSLC idealizations are microscopically justified, and supplies a reference point against which finite-time NELC power and performance can be evaluated. In the specific model studied, internal coupling broadens operational domains, enables engine or refrigerator action where the uncoupled Otto cycle is non-operational, and enhances efficiency and COP beyond standard Otto bounds while staying below Carnot limits (Gao et al., 2 Mar 2026).

More broadly, the recurring appearance of ELC-like constructions across disparate systems suggests a common research role: they formalize how a system “equilibrates” not necessarily to a static fixed point, but to a stable periodic object determined by the governing open, stochastic, switched, or controlled dynamics. That implication is explicit in some papers and interpretive in others, but it captures why the concept is useful: it extends the language of equilibration from time-independent steady states to periodic steady operation (Gao et al., 2 Mar 2026, Chan et al., 2015, Makarenkov, 2017, Louca, 2015, Makarenkov et al., 4 Nov 2025).