多原点曲率（MOC）与最大信息效率（MIE）框架下双温库两能级系统非平衡稳态分布的统一解析解

作者：张苏杭（Le Zhang, Bosley Zhang）

独立数学物理与理论物理研究者

通讯邮箱：zhang34269@zohomail.cn

摘要

本文基于多原点曲率（Multi-Origin Curvature, MOC）几何与最大信息效率（Maximum Information Efficiency, MIE）极值原理，构建双温库耦合两能级系统非平衡稳态分布的统一解析框架。将两能级映射为MOC空间中的独立原点，非平衡稳态等价于MIE泛函在拓扑约束下的极值点。通过严格变分推导，得到了模型无关的普适解析分布公式，该公式可在平衡态极限下严格退化为玻尔兹曼分布，解决了传统主方程方法依赖耦合模型、形式复杂、缺乏统一解析形式的局限。本文结果表明，MOC-MIE框架可为非平衡统计力学问题提供简洁、自洽且可检验的解析描述。

关键词：多原点曲率；最大信息效率；非平衡统计力学；两能级系统；非平衡稳态；解析解

一、引言

两能级系统与多温库耦合的非平衡稳态问题，是非平衡统计力学的经典模型，广泛应用于量子热机、量子点输运与分子马达等研究。传统求解方法基于量子主方程或速率方程，需引入跃迁速率、耦合强度等微观动力学参数，得到的分布形式复杂且强依赖具体耦合模型，缺乏普适的解析形式，难以建立平衡态与非平衡态的统一描述。

现有方法的核心局限在于，从微观动力学过程出发进行正向推导，未对系统的拓扑几何约束与全局极值条件进行公理化刻画。本文基于前期建立的MOC几何公理体系与MIE极值原理，将两能级系统映射为MOC空间中的双原点结构，将非平衡稳态条件等价为MIE泛函的极值点，绕开微观动力学细节，直接推导出普适的非平衡稳态分布解析形式，并验证其平衡态极限的自洽性。

二、理论框架与模型构建

2.1 多原点曲率（MOC）空间建模

将两能级系统的两个能级能量\varepsilon_1、\varepsilon_2，定义为MOC空间中的两个独立原点\mathcal{O}_1、\mathcal{O}_2，系统处于两能级的概率分别为p_1、p_2，满足归一化约束：

p_1 + p_2 = 1

系统的微观构型数由双原点拓扑关联唯一确定，记为\Omega_{\text{MOC}}，对应的MOC约束项定义为：

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

该约束项包含系统的全部拓扑几何信息，替代传统微观动力学描述。

2.2 最大信息效率（MIE）泛函构造

非平衡稳态对应MIE泛函在全局约束下的极值点。系统的MIE泛函形式为：

\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

其中\beta_1 = 1/(k_B T_1)、\beta_2 = 1/(k_B T_2)分别为两个热库的逆温度，k_B为玻尔兹曼常数。泛函第一项为信息熵项，第二项为MOC几何约束项，第三、四项为能量约束项。

三、变分推导与解析解

3.1 极值条件推导

对MIE泛函关于概率变量p_1、p_2变分，结合极值约束结构（ECS）的唯一性条件，得到极值条件：

\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, \quad i=1,2

联立两式消去常数项与MOC约束项的导数，可得：

\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 拓扑约束化简与解析解

双原点拓扑结构的对称性要求，MOC约束项的导数差仅由能级间跃迁速率的几何平均决定，即：

\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}}}

其中\gamma_{12}、\gamma_{21}为能级间的跃迁速率。代入极值条件，整理得到非平衡稳态概率分布比值的普适解析形式：

\boxed{

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

}

四、自洽性验证与对比分析

4.1 平衡态极限退化验证

当双热库温度相等，即T_1 = T_2 = T，对应\beta_1 = \beta_2 = \beta时，跃迁速率满足细致平衡条件：

\frac{\gamma_{12}}{\gamma_{21}} = \exp\left[ -\beta (\varepsilon_2 - \varepsilon_1) \right]

代入解析解，可得：

\frac{p_2}{p_1} = \exp\left[ -\frac{\varepsilon_2 - \varepsilon_1}{k_B T} \right]

该结果严格退化为平衡态玻尔兹曼分布，证明了解析解的自洽性与极限相容性。

4.2 与传统方法的对比

传统主方程方法得到的分布形式为：

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

该形式复杂且依赖耦合模型，而本文得到的解析解形式简洁、模型无关，同时包含了跃迁速率与温度的几何平均效应，具备更强的普适性。

五、讨论

本文结果表明，MOC-MIE框架可从拓扑几何与全局极值层面，直接得到双温库两能级系统非平衡稳态的普适解析解，绕开了传统方法中复杂的微观动力学建模过程。该解析解不仅形式简洁，且可在平衡态极限下严格退化为已知结果，证明了理论框架的自洽性。

该方法的优势在于：

1. 普适性：解析形式不依赖具体耦合模型，适用于任意两能级系统与多温库耦合场景；

2. 可检验性：结果可通过量子点实验或数值模拟直接验证；

3. 可拓展性：可推广至多能级、多热库系统的非平衡稳态问题。

六、结论

本文基于MOC-MIE框架，构建了双温库耦合两能级系统非平衡稳态的统一解析模型，得到了普适的分布公式：

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

该公式形式简洁、模型无关，且在平衡态极限下严格退化为玻尔兹曼分布，解决了传统方法解析性不足的局限。本文验证了MOC-MIE框架在非平衡统计力学问题中的有效性，为非平衡系统的统一描述提供了新的理论工具。

参考文献

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

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

[3] 张苏杭. 多原点曲率(MOC)几何基础公理体系[J/OL]. 预印本平台, 2026.

[4] 张苏杭. 多原点曲率与最大信息效率框架下经典与量子统计分布的统一理论[J/OL]. 预印本平台, 2026.