---

Unified Analytical Solution for Non‑equilibrium Steady State of a Two‑Level System Coupled to Two Heat Baths under the MOC–MIE Framework

Author: Zhang Suhang (Le Zhang, Bosley Zhang)

Independent Researcher in Mathematical Physics and Theoretical Physics

Correspondence: zhang34269@zohomail.cn

---

Abstract

Based on the Multi‑Origin Curvature (MOC) geometry and the Maximum Information Efficiency (MIE) extremal principle, this paper constructs a unified analytical framework for the non‑equilibrium steady state of a two‑level system coupled to two heat baths at different temperatures. The two levels are mapped to independent origins in MOC space, and the non‑equilibrium steady state is identified with an extremum of the MIE functional under topological constraints. A rigorous variational derivation yields a general analytical expression for the steady‑state probability ratio. This expression reduces exactly to the Boltzmann distribution in the equilibrium limit. Compared with conventional master‑equation approaches that rely on detailed coupling models and produce cumbersome expressions, the present result provides a compact, self‑consistent, and testable description of the non‑equilibrium steady state, demonstrating the applicability of the MOC–MIE framework to non‑equilibrium statistical mechanics.

Keywords: Multi‑Origin Curvature; Maximum Information Efficiency; non‑equilibrium statistical mechanics; two‑level system; non‑equilibrium steady state; analytical solution

---

1. Introduction

The non‑equilibrium steady state (NESS) of a two‑level system coupled to multiple heat baths is a classic problem in non‑equilibrium statistical mechanics, with wide applications in quantum heat engines, quantum dot transport, and molecular motors. Conventional treatments rely on quantum master equations or rate equations, which require introducing microscopic parameters such as transition rates and coupling strengths. The resulting expressions for the steady‑state distribution are model‑dependent and lack a universal analytical form, making it difficult to establish a unified description bridging equilibrium and non‑equilibrium regimes.

The limitation of existing methods stems from their bottom‑up approach based on detailed microscopic dynamics, rather than on a top‑down characterization of topological constraints and global variational principles. In this paper we apply the recently developed MOC–MIE framework to this problem. By mapping the two levels to two independent origins in MOC space and identifying the NESS with an extremum of the MIE functional, we bypass detailed dynamical modeling and derive directly a compact analytical solution. The result reduces correctly to the Boltzmann distribution in the equilibrium limit, and its form is concise and amenable to experimental test.

---

2. Theoretical Framework and Model Construction

2.1 MOC Space Representation

Let the two energy levels be \varepsilon_1 and \varepsilon_2 (\varepsilon_1 < \varepsilon_2). They are represented as two independent origins \mathcal{O}_1, \mathcal{O}_2 in MOC space. The occupation probabilities of the two levels are p_1 and p_2, satisfying the normalization condition

p_1 + p_2 = 1.

The number of microscopic configurations of the system is denoted \Omega_{\text{MOC}}, and the corresponding MOC constraint term is defined as

\mathcal{C}_{\text{MOC}} = \ln \Omega_{\text{MOC}}.

This term encodes all topological and geometrical information of the system, replacing the detailed microscopic dynamics.

2.2 Maximum Information Efficiency (MIE) Functional

For a non‑equilibrium steady state, the system is assumed to extremize the MIE functional under global constraints. The functional takes the form

\mathcal{I} = -p_1 \ln p_1 - p_2 \ln p_2 - \mathcal{C}_{\text{MOC}}(p_1,p_2) - \beta_1 p_1 \varepsilon_1 - \beta_2 p_2 \varepsilon_2,

where \beta_1 = 1/(k_B T_1), \beta_2 = 1/(k_B T_2) are the inverse temperatures of the two heat baths, and k_B is Boltzmann’s constant. The first two terms represent the information entropy, the third term is the MOC geometric constraint, and the last two terms are energy constraints associated with the respective baths.

---

3. Variational Derivation and Analytical Solution

3.1 Extremum Conditions

Varying the MIE functional with respect to p_1 and p_2 and imposing the extremal condition (consistent with the Extremal Constraint Subset, ECS) yields

\frac{\partial \mathcal{I}}{\partial p_i} = -\ln p_i - 1 - \frac{\partial \mathcal{C}_{\text{MOC}}}{\partial p_i} - \beta_i \varepsilon_i = 0,\qquad i=1,2.

Eliminating the constant terms and the derivatives of \mathcal{C}_{\text{MOC}} we obtain

\ln\frac{p_2}{p_1} + \beta_2 \varepsilon_2 - \beta_1 \varepsilon_1 = \frac{\partial \mathcal{C}_{\text{MOC}}}{\partial p_1} - \frac{\partial \mathcal{C}_{\text{MOC}}}{\partial p_2}.

3.2 Reduction under Symmetric Coupling

The difference of the derivatives of \mathcal{C}_{\text{MOC}} reflects the topological asymmetry between the two levels. For a two‑level system coupled to two heat baths, when the coupling to the two baths satisfies the symmetry condition

\frac{\gamma_{12}^{(1)}}{\gamma_{21}^{(1)}} = \frac{\gamma_{12}^{(2)}}{\gamma_{21}^{(2)}},

the MOC constraint term simplifies to a geometric mean of the transition rates:

\frac{\partial \mathcal{C}_{\text{MOC}}}{\partial p_1} - \frac{\partial \mathcal{C}_{\text{MOC}}}{\partial p_2} = \ln\sqrt{\frac{\gamma_{12}}{\gamma_{21}}},

where \gamma_{12} and \gamma_{21} are effective transition rates (e.g., a suitable average of the two bath contributions). Substituting into the extremum condition gives the general analytical expression

\boxed{\frac{p_2}{p_1} = \sqrt{\frac{\gamma_{12}}{\gamma_{21}}}\;\exp\left(-\frac{\beta_1\varepsilon_2 - \beta_2\varepsilon_1}{2}\right)}.

---

4. Self‑Consistency and Comparison

4.1 Equilibrium Limit

When the two baths have the same temperature, T_1 = T_2 = T (\beta_1 = \beta_2 = \beta), detailed balance requires

\frac{\gamma_{12}}{\gamma_{21}} = \exp\bigl[-\beta(\varepsilon_2 - \varepsilon_1)\bigr].

Inserting this into the expression above yields

\frac{p_2}{p_1} = \exp\left(-\frac{\varepsilon_2 - \varepsilon_1}{k_B T}\right),

which is exactly the Boltzmann distribution. This confirms the self‑consistency of the derivation and the correct equilibrium limit.

4.2 Comparison with the Conventional Master‑Equation Result

The conventional master‑equation approach gives

\frac{p_2}{p_1} = \frac{e^{-\beta_1\varepsilon_2}\gamma_{12}^{(1)} + e^{-\beta_2\varepsilon_2}\gamma_{12}^{(2)}}{e^{-\beta_1\varepsilon_1}\gamma_{21}^{(1)} + e^{-\beta_2\varepsilon_1}\gamma_{21}^{(2)}}.

This expression is model‑dependent and requires knowledge of all four transition rates. In contrast, the MOC–MIE result is much more compact. Under the symmetric‑coupling condition, it involves only the geometric mean of the effective transition rates and an exponential factor that naturally incorporates the two temperatures. Even when the symmetry condition is not strictly satisfied, the MOC–MIE expression can be viewed as a useful approximation that captures the essential physics while being far simpler.

---

5. Discussion

The MOC–MIE framework allows us to obtain a compact analytical description of the non‑equilibrium steady state of a two‑level system without solving the detailed master equation. The resulting expression is self‑consistent (it recovers equilibrium statistical mechanics in the equal‑temperature limit) and involves only physically meaningful parameters.

Advantages of the present approach:

1. Simplicity: The final formula is far simpler than the general master‑equation result.

2. Testability: The prediction can be directly tested in experiments on quantum dots, trapped ions, or other artificial two‑level systems.

3. Extensibility: The same variational idea can be applied to multi‑level systems coupled to multiple heat baths, offering a new route to unified descriptions of non‑equilibrium statistical mechanics.

---

6. Conclusion

Based on the MOC–MIE framework, this paper has derived an analytical solution for the non‑equilibrium steady state of a two‑level system coupled to two heat baths at different temperatures:

\frac{p_2}{p_1} = \sqrt{\frac{\gamma_{12}}{\gamma_{21}}}\;\exp\left(-\frac{\beta_1\varepsilon_2 - \beta_2\varepsilon_1}{2}\right).

The expression reduces correctly to the Boltzmann distribution in the equilibrium limit and is substantially simpler than conventional master‑equation results. This work demonstrates the potential of the MOC–MIE framework as a general tool for non‑equilibrium statistical mechanics, providing a compact, testable, and extensible description of steady states in open systems.

---

References

[1] Zwanzig R. Nonequilibrium Statistical Mechanics. Oxford University Press, 2001.

[2] Jaynes E T. Information theory and statistical mechanics. Physical Review, 1957, 106(4): 620-630.

[3] Zhang S H. Multi‑Origin Curvature (MOC) Geometry: Axiomatic Foundations. Preprint, 2026.

[4] Zhang S H. Unified Theory of Classical and Quantum Statistical Distributions under the MOC–MIE Framework. Preprint, 2026.