本征微观态

理论背景

系综（统计系综）

如果我们用概率来描述热系统，采用的方法一般是，设想一次又一次地重复某一个实验，来测量某一个系统的某一种性质，因为我们无法控制它们的微观性质（由系统的微观态所描述），在尝试表述这个方法时，乔赛亚·威拉德·吉布斯（Josiah Willard Gibbs）在1902年提出了系综的概念 。这是一种理想化的方法，在该方法中他考虑对系统进行大量想象的“影印”，其中每一个都代表了该系统所处的一个可能状态。

在物理学，特别是在统计物理学中，在统计物理中，系综代表一定条件下，一个体系的大量可能状态的集合。换句话说，系综是系统状态的一个概率分布。对一相同性质的体系，其微观状态（比如每个粒子的位置和速度）仍然可以大不相同。更进一步地说，统计系综是统计力学中用来描述单一系统的一组粒子系统 。

下面举一个例子来说明这样的表述。考虑抛一枚硬币的实验，这样一个简单的实验只有两种可能的结果，“正面”或“反面”。原则上，如果我们能够确切地知道硬币是如何被抛出的，以及与硬币和桌子相互作用力等等信息，那么只要根据经典力学的理论进行一定的计算，实验的结果应该是完全可以预测的。实际上，关于这个实验详细的、精确的信息是无法获取的。所以对于某一次实验结果，我们不可能作出唯一的预测，可是实验的统计表述却是比较简单的。

为了展示抛掷一枚硬币的概率，我们考虑由N枚相同的硬币组成的系综（N是一个非常大的的值），当每一枚硬币都被抛掷后，系综的“面貌”在就在这张图中被展示出来了。

我们只要考虑由很大数目，[math]\displaystyle{ N }[/math]枚相似的硬币组成的一个系综，当这些硬币以同样的方式抛出，我们可以数出结果中硬币正反面的个数，进而得到正面的概率[math]\displaystyle{ p }[/math]和反面的概率[math]\displaystyle{ q }[/math]。统计理论希望能够预测这些概率。

现在考虑稍微复杂一点的掷N枚硬币的实验，由于抛掷任何一枚硬币都有两个可能的结果，那么掷N枚硬币就可以出现[math]\displaystyle{ 2×2×2×…×2=2^N }[/math]个可能结果中的任何一个。如果不是只讨论一组[math]\displaystyle{ N }[/math]枚硬币，而是考虑[math]\displaystyle{ N }[/math]个这样的组（每组有[math]\displaystyle{ N }[/math]枚硬币）所组成的系综，每组都以相似的方式抛掷硬币，那么值得我们探究的问题便是，[math]\displaystyle{ 2^N }[/math]个可能结果中，任何一个特殊的结果在系综中出现的概率为多大。

“平衡系综”的概念对统计系综的许多应用至关重要。如果每一时刻体系的统计系综中，呈现任一特殊事件的体系数目是一样的（或等价地表示为：如果这个系综中任一特殊事件出现的概率与时间无关），那么就说这个系综是与时间无关的。换句话说，尽管一个机械系统肯定会随着时间的推移而演变，但系综并不一定会发生演变。事实上，如果系综包含了系统的所有过去和未来阶段，那么它就不会演变。这样的统计描述就为平衡提供一个非常清楚的定义：如果孤立宏观体系的一个统计系综是与时间无关的，那么这样一个体系就称为处于平衡。

吉布斯定义了三种主要的系综：

（1）微正则系综（microcanonical ensemble）：系综里的每个体系具有相同的能量。

（2）正则系综（canonical ensemble）：系综里的各体系可以和外界环境交换能量，这种能量交换将确定（并且定义）了系统的温度。

（3）巨正则系综（grand canonical ensemble）：是正则系综的推广，各体系可以和外界环境交换能量和粒子，但系综内各个体系有相同的温度和化学势。

在系综中，物理量的变化范围与其本身大小的比值会随着体系的增大而减小。那么对于一个宏观体系来说，从各种系综计算出的物理量的差异将趋于零。

相变

在化学、热力学和其他许多相关领域，相变（或相位变化）是指介质的一种状态（由某些参数确定）和另一种状态（参数值不同）之间的过渡物理过程。这个术语通常用于指固体、液体和气体，等离子体（少数情况下）等物质的基本状态之间的变化。

此图显示了不同相变的命名法

举例来说，热力学系统的一个阶段和物质状态具有统一的物理特性。在某一介质的相变过程中，由于外部条件的变化，如温度、压强或其他条件的变化，介质的某些属性发生的变化，往往是不连续的。比如液体在加热到沸点时，可能变成气体，导致体积的突然变化。这种对发生转变的外部条件的测度，称为相变。相变通常发生在自然界中，并且今天在许多技术中也常常被使用。

保罗·埃伦费斯特（Paul Ehrenfest）提出了如下关于相变的一种分类 ：相变的级（order）是在[math]\displaystyle{ T_C }[/math]处吉布斯自由能[math]\displaystyle{ G }[/math]（或者[math]\displaystyle{ \mu }[/math]）的微分显示不连续性的最低的阶数。那么一级相变（rstorder phase transition）包含潜热，因为熵[math]\displaystyle{ S }[/math]（[math]\displaystyle{ G }[/math]的一阶微分）此时显示出不连续性，体积[math]\displaystyle{ V }[/math]（也是[math]\displaystyle{ G }[/math]的一阶微分）也显示出一个并不连续的跳跃 。热容[math]\displaystyle{ C }[/math]是[math]\displaystyle{ G }[/math]的二阶微分，并且也显示出一个尖锐的跳跃，压缩系数[math]\displaystyle{ \beta }[/math]也是如此。一级相变的例子包括固液相变、固气相变以及气液相变。

根据埃伦费斯特的分类，二级相变（second-order phase transition）没有潜热，因为熵[math]\displaystyle{ S }[/math]没有显示出不连续性，并且体积[math]\displaystyle{ V }[/math]也没有不连续性，但是如热容[math]\displaystyle{ C }[/math]和压缩系数[math]\displaystyle{ \beta }[/math]等量却依然不连续 。二级相变的例子包括超导相变，或者[math]\displaystyle{ \beta }[/math]黄铜中的有序无序转变。

然而，在研究相变时我们所使用的这种方法，存在一个很大的问题。我们在热力学中常所作这样一个近似，粒子的数目很大，以至于平均性质（如压强和密度）能够被很好地定义，但这种近似在相变时将会失效。在接近相变时，涨落增加，所以在非常接近相变温度时，系统的行为并不遵循我们所预期的分析。在所有的长度标度下，这个临界区域会由涨落所表征。

不符合埃伦费斯特分类的相变的第一个例子是1944年由拉尔斯·昂萨格（Lars Onsager）发现的Ising模型的精确解。在接下来的几十年里，埃伦费斯特分类法被一个简化的分类方案所取代，这个方案不在此过多赘述，感兴趣的读者可以参见相变。

临界点与临界现象

在热力学中，一个临界点（或临界状态）就是相平衡曲线的终点。最突出的例子是液-汽临界点，它是压力-温度曲线的终点，它指明了液体和其蒸汽可以共存的条件。在较高的温度下，气体不能仅靠压力来液化。在由临界温度[math]\displaystyle{ T_C }[/math]和临界压力[math]\displaystyle{ p_C }[/math]定义的临界点，相边界消失。其他例子包括混合物中的液-液临界点，以及在没有外部磁场的情况下的铁磁体-准磁体转变 。

该图为典型的相图。其中虚线部分表示了水的反常行为

为使表述简单明晰，临界点的一般概念我们通过液体-蒸汽临界点来介绍。

右图显示了纯物质的P-T图。众所周知，固相、液相和汽相通过相边界分离，即两相可以共存的压力-温度组合。在三相点，所有三个相可以共存。然而，在临界温度[math]\displaystyle{ T_C }[/math]和临界压力[math]\displaystyle{ p_C }[/math]时，液-汽边界终止于一个点，也就是临界点。

在临界点附近，液体和蒸汽的物理性质发生了巨大变化，两相变得越来越相似。在临界点，只有一个相存在。这意味着在临界点：

[math]\displaystyle{ \left(\frac{\partial p}{\partial V}\right)_T = 0, }[/math]

[math]\displaystyle{ \left(\frac{\partial^2p}{\partial V^2}\right)_T = 0. }[/math]

由右图可见，液相和气相区域之间的边界线在临界点终结，因此有可能通过选取一条穿过相图、避免不连续变化的路径“逃脱”清晰的相变。对于高于临界温度（如水为647K）的温度，气态和液态仅可以通过它们的密度来区分。气态和液态之间的转变不涉及对称性的变化，所以可通过围绕临界端点改变状态，从而避免明显的相变。与此相比较，固液相变涉及对称性的变化，因此熔化曲线没有临界点。关于更多涉及对称性的变化的定义与示例，参见对称性破缺。

平衡态与非平衡态

“平衡”这个词意味着一种平衡的状态。平衡热力学，通常被认为起源于对卡诺循环的分析。在这个典型的系统里，比如一缸气体，最初处于其自身的内部热力学平衡状态，通过燃烧反应的热输入而失去平衡。然后通过一系列的步骤，随着系统进入最终的平衡状态，功被释放出来。

在平衡态下，系统内的热力学势或驱动力处于完全平衡状态。平衡热力学的一个核心目标是：给定一个处于明确的热力学平衡初始状态的系统，受到被精确规定的约束，计算当约束被外部施加的干预改变时，一旦系统达到新的平衡，系统的状态将会是什么。平衡态在数学上是通过寻求热力学势函数的极值来确定的，其性质取决于施加在系统上的约束。例如，在恒定温度和压力下的化学反应，将在吉布斯自由能[math]\displaystyle{ G }[/math]最小和熵[math]\displaystyle{ S }[/math]最大时达到平衡 。

平衡热力学与非平衡热力学不同，对于后者而言，被分析系统的状态通常不是均匀的，而是在能量[math]\displaystyle{ E }[/math]、熵[math]\displaystyle{ S }[/math]和温度[math]\displaystyle{ T }[/math]分布上会发生局部变化。相比之下，在平衡热力学中，系统的状态将被认为是均匀的，由温度、压力或体积等量宏观定义的。系统是根据从一个平衡状态到另一个平衡状态的变化来研究的，这样的变化被称为热力学过程。

几乎所有在自然界中发现的系统都不是在热力学平衡中，因为它们正在随着时间变化或者可以被触发而发生变化，并且不断地和其他系统交换物质和能量以及参与化学反应。然而，某些系统和过程在某些程度上足够接近于热力学平衡态，允许目前已有的非平衡态热力学对其进行可行的精确描述。然而，在许多自然系统和过程中，由于非变分动力学的存在，使得自由能的概念并不存在，因此总是远远超出非平衡热力学方法的范围 。

对于平衡态的研究，使用了“准静态过程”的理想化概念。准静态过程是一种概念上沿着热力学平衡状态连续路径的平滑数学过程。它只可能是微分几何的练习，而不是现实中可能发生的过程 。此外，非平衡态热力学试图描述连续的时间进程，这需要它的状态变量与平衡态热力学的状态变量之间有非常密切的联系。这极大地限制了非平衡态热力学的范围，并对其概念框架提出了严格的要求 。

不同于前文所述的平衡态，非平衡态要复杂得多，它可能存在更多广延量的涨落。边界条件施加给它们某些强度量，如温度梯度或形变集体运动（剪切运动、涡旋等），通常称为热力学力。在这里，我们必须强调，这里并没有像平衡态热力学中熵的热力学第二定律那样定义能量的静态非平衡性质的一般定律。这就是为什么在这种情况下，应该考虑一个更一般的Legendre变换（也就是拓展的马休势，extended Massieu function）。

微观态与宏观态

我们不妨用一个简单的例子来介绍微观态与宏观态。

设想一个大盒子里，放有100个完全相同的硬币。盖上盒盖后，用力并持续足够长的时间摇晃盒子，随后打开盒盖朝里看，有些硬币正面朝上，有些硬币反面朝上 ，有大量可以获得的可能组态（准确地说，应该是[math]\displaystyle{ 2^{100} }[/math]种，大约为[math]\displaystyle{ 10^{30} }[/math])。这里我们假定，这些不同组态中的每一种，均是等可能出现的。因此每种可能组态出现的概率约为[math]\displaystyle{ 10^{-30} }[/math]。我们称上述每一种特定的组态（configuration），为该系统的一个微观态（microstate）。这些微观态的某一个例子是：“一号硬币正面朝上，二号硬币正面朝上，三号硬币反面朝上，……，一百号硬币反面朝上”。为了辨别一个微观态，我们可能需要单独地辨别每一个硬币，这的确令人烦躁，但是这毕竟只是在简单地数数：有多少硬币正面朝上，有多少硬币反面朝上（例如，有53枚正面朝上，47枚反面朝上)。这样的分类称为该系统的一个宏观态（macrostate）。但值得注意的是，每个宏观态并不是等可能出现的。例如，在约为[math]\displaystyle{ 10^{30} }[/math]个可能的组态（微观态）中，50枚硬币正面朝上，50枚硬币反面朝上的组态数为[math]\displaystyle{ \frac{100!}{50!×50!}≈4×10^{27} }[/math]；53枚硬币正面朝上，47枚硬币反面朝上的组态数为[math]\displaystyle{ \frac{100!}{53!×47!}≈3×10^{27} }[/math]；100枚硬币正面朝上，0枚硬币反面朝上的组态数为1。

这样看来，100枚硬币正面全朝上的结果是不太可能发生的，因为这个宏观态只含有一个微观态。当然，有53枚正面和47枚反面的一个特定微观态也同样也不太可能发生，这是因为还有将近[math]\displaystyle{ 3×10^{27} }[/math]个有53枚正面和47枚反面的、看上去与之极端相似的其他微观态存在。

这个简单的例子说明了两个关键点：

（1）可以用数量巨大的、同等可能的微观态描述一个系统；

（2）实际测量的是系统宏观态的一个性质，但各个宏观态并不是同等可能出现的，因为不同宏观态对应不同数量的微观态。

系统最可能所处的宏观态，是对应于最多微观态数的宏观态。

热力学系统的行为，与上面考虑的例子中的情况非常类似。要指定一个热力学系统的一个微观态，就需要给出系统中每一个原子的微观位形 （configuration，可能是位置和速度或者能量）。但值得注意的是，我们通常不可能测量出系统处在哪一个微观态。另一方面，仅仅给出一个热力学系统的宏观性质（如压强，总能量或者体积），我们就能指定该系统的个宏观态。举例来说，在体积[math]\displaystyle{ 1m^3 }[/math]中压强为105Pa的气体的一个宏观位形，通常会与大量的微观态相联系。

本征微观态及其演化

对于一个由[math]\displaystyle{ N }[/math]个主体组成的复杂系统，我们可以从实验测量或计算机模拟中获得主体的状态。依次使用时间[math]\displaystyle{ t=1, 2, … ,M }[/math]的状态，我们可以得到主体[math]\displaystyle{ i=1,2, … , N }[/math]的状态序列[math]\displaystyle{ S_i(t) }[/math]。

一个主体[math]\displaystyle{ i }[/math]的平均状态是 [math]\displaystyle{ \begin{eqnarray}\langle {S}_{i}\rangle =\displaystyle \frac{1}{M}\displaystyle \sum _{t=1}^{M}{S}_{i}(t)\end{eqnarray} }[/math]。在某个时间[math]\displaystyle{ t }[/math]，主体[math]\displaystyle{ i }[/math]有一个波动 [math]\displaystyle{ \begin{eqnarray}\delta {S}_{i}(t)={S}_{i}(t)-\langle {S}_{i}\rangle\end{eqnarray} }[/math]。这里定义一个具有所有主体的波动的微观态，它由一个[math]\displaystyle{ N }[/math]维矢量表示： [math]\displaystyle{ \begin{eqnarray}\delta {\boldsymbol{S}}(t)=\left[\begin{array}{c}\delta {S}_{1}(t)\\ \delta {S}_{2}(t)\\ \vdots \\ \delta {S}_{N}(t)\end{array}\right]\end{eqnarray} }[/math]。

有了[math]\displaystyle{ M }[/math]个微观态，我们可以组成一个复杂系统的统计系综。这个系综由一个[math]\displaystyle{ N×M }[/math]的矩阵[math]\displaystyle{ \boldsymbol{A} }[/math]来描述，其元素为[math]\displaystyle{ \begin{eqnarray}{A}_{{it}}=\displaystyle \frac{\delta {S}_{i}(t)}{\sqrt{{C}_{0}}}\end{eqnarray} }[/math]，其中 [math]\displaystyle{ {C}_{0}={\sum }_{t=1}^{M}{\sum }_{i=1}^{N}\delta {S}_{i}^{2}(t) }[/math]，[math]\displaystyle{ \boldsymbol{A} }[/math]的列序与微观态的演变相一致。

[math]\displaystyle{ t }[/math]和[math]\displaystyle{ t^{\prime} }[/math]的微观态之间的相关性由它们的矢量乘积定义：[math]\displaystyle{ \begin{eqnarray}{C}_{{tt}^{\prime} }=\delta {\boldsymbol{S}}{\left(t\right)}^{{\rm{T}}}\cdot \delta {\boldsymbol{S}}(t^{\prime} )=\displaystyle \sum _{i=1}^{N}\delta {S}_{i}(t)\delta {S}_{i}(t^{\prime} )\end{eqnarray} }[/math] 。 以[math]\displaystyle{ {C}_{tt}^{\prime} }[/math]作为其元素，我们可以得到一个[math]\displaystyle{ M×M }[/math]的微观态相关矩阵：[math]\displaystyle{ \begin{eqnarray}{\boldsymbol{C}}={C}_{0}{{\boldsymbol{A}}}^{{\rm{T}}}\cdot {\boldsymbol{A}}\end{eqnarray} }[/math]，其轨迹[math]\displaystyle{ Tr\boldsymbol{C}={\sum }_{t=1}^{M}{C}_{tt}={C}_{0} }[/math] 。相关矩阵[math]\displaystyle{ \boldsymbol{C} }[/math]有[math]\displaystyle{ M }[/math]个特征向量[math]\displaystyle{ V_J }[/math]，其中[math]\displaystyle{ J = 1, 2, ⋯ , M }[/math]，我们可以用它们组成一个[math]\displaystyle{ M×M }[/math]的单元矩阵 [math]\displaystyle{ \boldsymbol{V}=\left[\boldsymbol{V}_1 \boldsymbol{V}_2 \ldots \boldsymbol{V}_M\right] }[/math]

此外，我们在这里考虑动态微观态[math]\displaystyle{ δS_i }[/math]和[math]\displaystyle{ δS_j }[/math]之间的相关性：[math]\displaystyle{ \begin{eqnarray}{{\boldsymbol{S}}}_{j}^{{{T}}}=\displaystyle \sum _{t=1}^{M}\delta {S}_{i}(t)\delta {S}_{j}(t)\end{eqnarray} }[/math]。以[math]\displaystyle{ K_{ij} }[/math]为元素，我们可以得到一个[math]\displaystyle{ N×N }[/math]的动态微观态的相关矩阵： [math]\displaystyle{ \begin{eqnarray}{\boldsymbol{K}}={C}_{0}{\boldsymbol{A}}\cdot {\boldsymbol{A}}^{T}\end{eqnarray} }[/math]，其轨迹[math]\displaystyle{ Tr\boldsymbol{K}={\sum }_{t=1}^{N}{K}_{ii}={C}_{0} }[/math]。相关矩阵[math]\displaystyle{ \boldsymbol{K} }[/math]有[math]\displaystyle{ N }[/math]个特征向量[math]\displaystyle{ U_I }[/math]，其中[math]\displaystyle{ I = 1, 2, ⋯ , N }[/math]，我们可以用它们组成一个[math]\displaystyle{ N×N }[/math]的单元矩阵 [math]\displaystyle{ \boldsymbol{U}=\left[\boldsymbol{U}_1 \boldsymbol{U}_2 \ldots \boldsymbol{U}_N\right] }[/math]。

根据奇异值分解，系综矩阵[math]\displaystyle{ \boldsymbol{A} }[/math]可以被分解为[math]\displaystyle{ \begin{eqnarray}{\boldsymbol{A}}={\boldsymbol{U}}\cdot {\boldsymbol{\Sigma }}\cdot {{\boldsymbol{V}}}^{{\rm{T}}}，\end{eqnarray} }[/math]

对[math]\displaystyle{ \boldsymbol{A} }[/math]进行奇异值分解，统计系综可被分解为：

[math]\displaystyle{ A=\sum_{i=1}^r \sigma_i U_i \otimes V_i }[/math]

这里的[math]\displaystyle{ r=min (N, M) }[/math]，其中，[math]\displaystyle{ U_i }[/math]为本征微观态，[math]\displaystyle{ V_i }[/math]为该本征微观态遵循的时间演化，⊗表示克罗内克积，其中[math]\displaystyle{ \sigma_i }[/math]表示[math]\displaystyle{ U_i }[/math]在系综中的概率，且满足归一化条件

[math]\displaystyle{ \sum_{i=1}^r \sigma_i^2=1 }[/math]

用这样的方法，可以将原来相互关联的微观态转变为相互独立的本征微观态，就可以将最初的微观态用本征微观态进行线性组合，这个线性组合的大小与本征值[math]\displaystyle{ \sigma_i }[/math]（即权重因子）成正比，即[math]\displaystyle{ \sigma_i }[/math]越大，占比越多。

演化的傅里叶谱分析

第[math]\displaystyle{ I }[/math]个本征微观态的时间演化由[math]\displaystyle{ {{\boldsymbol{V}}}_{I}^{{\rm{T}}}=\left[{V}_{1I},{V}_{2I},\cdots ,{V}_{{MI}}\right] }[/math]来进行描述。为了获得更多关于本征微观态演变的信息，这里对[math]\displaystyle{ V }[/math]进行傅里叶级数展开，有[math]\displaystyle{ V_{t I}=\sum_{n=0}^{M-1} b_{n I} \mathrm{e}^{\mathrm{i} \frac{2 \pi}{M} n t} }[/math]，其中傅里叶系数为[math]\displaystyle{ b_{n I}=\frac{1}{M} \sum_{t=0}^{M-1} V_{t I} \mathrm{e}^{-\mathrm{i} \frac{2 \pi}{M} n t} }[/math]。

为了表征不同频率的比例，这里引入了傅里叶功率谱密度（Fourier power spectrum density）[math]\displaystyle{ I_n=\frac{\left|b_n^I\right|^2}{\sum_{n=0}^{M-1}\left|b_n^I\right|^2} }[/math]，这将可以展示其周期[math]\displaystyle{ T_n=M/n }[/math]。

本征微观态和相变的凝聚

一般来说，在一个无序的复杂系统中并不存在一个主要的本征值和微观态，但当我们将微观态转变为本征微观态时，原来的每个微观态都可以看作一系列的本征微观态的线性求和。

如果一个概率振幅[math]\displaystyle{ σ }[/math]成为主导，以至于它在[math]\displaystyle{ M → ∞ }[/math]和[math]\displaystyle{ N → ∞ }[/math]时能够取得一个有限的极限，那么在统计集合中就会出现本征微观态[math]\displaystyle{ {\boldsymbol{A}}_{I}^{e} }[/math]的凝聚，这类似于玻色气体里的玻色-爱因斯坦凝聚。一个本征微观态的凝聚，表现为有限比例的无限本征微观态共享这个本征微观态。一个新的相将出现在系统中，而新相的特征是本征微观态[math]\displaystyle{ U_I }[/math]，新相的演变，可以通过[math]\displaystyle{ V_I }[/math]来描述。

应用

在平衡系统中

伊辛模型（Ising Model），是用来解释铁磁系统相变的一个简单模型。在相当长时间的探索中，人们逐渐认识到了它作为相变模型的普适性，并且发现伊辛模型可以用来对几乎所有有趣的热力学现象进行建模，甚至包括在物理之外的其他学科中。毫不夸张地说，伊辛模型是统计物理中迄今为止唯一的一个同时具备表述简单、内涵丰富、应用广泛这三种优点的模型。

我们考虑有这么一个模型：具有方形几何结构的、具有周期性边界条件的三维伊辛模型（Ising model）。[math]\displaystyle{ N = L × L × L }[/math]的哈密顿方程是 [math]\displaystyle{ \begin{eqnarray}H=-J\displaystyle \sum _{\langle i，j\rangle }^{N}{S}_{i}{S}_{j}\end{eqnarray} }[/math]， 其中[math]\displaystyle{ S_i=±1 }[/math]，[math]\displaystyle{ ⟨i, j⟩ }[/math]表示对最近邻域的所有自旋进行求和。这个模型的系统大小为[math]\displaystyle{ L }[/math]，还原温度[math]\displaystyle{ T = k*BT/J }[/math]。我们可以使用Wolff算法来获得系统的微观态，这个算法可以翻转一簇自旋而不是单一自旋。

在有限伊辛系统中，不存在对称性破缺。因此[math]\displaystyle{ ⟨Si⟩ = 0 }[/math]，[math]\displaystyle{ δSi(t) = Si(t) }[/math]。在时间[math]\displaystyle{ t处 }[/math]的微观态是 [math]\displaystyle{ \delta \boldsymbol{S}(t)=\left[\begin{array}{c} \delta S_1(t) \\ \delta S_2(t) \\ \vdots \\ \delta S_N(t) \end{array}\right] }[/math]。利用这些微观态，我们得到系综矩阵[math]\displaystyle{ \boldsymbol{A} }[/math]，其元素为[math]\displaystyle{ A_{i t}=\frac{\delta S_i(t)}{\sqrt{C_0}} }[/math]，其中[math]\displaystyle{ C_0 = M·N }[/math]。

在（A）[math]\displaystyle{ T^{*}=5.5116 }[/math]、（B）[math]\displaystyle{ T^{*}=4.5116 }[/math]、（C）[math]\displaystyle{ T^{*}=3.5116 }[/math]时，与[math]\displaystyle{ U1,U2,U3 }[/math]等距的四张截面图

下面我们继续通过研究本征微观态[math]\displaystyle{ U1 }[/math]的空间分布，来观察相变的特性。对于一个在三维空间的本征微观态，我们给出了四个等距的截面图来展示其空间分布。其中，在

[math]\displaystyle{ T^{*}=5.5116 }[/math]时（高于临界温度），从图中可以看到，这些本征微观态中的自旋簇都具有微小尺寸，且在空间中随机分布。在[math]\displaystyle{ T^{*}=4.5116 }[/math]时（接近临界温度）以及在[math]\displaystyle{ T^{*}=3.5116 }[/math]时（低于临界温度），最大的本征微观态EM1将产生一个凝聚。这表明，当本征微观态的概率变得有限时，就会出现一个铁磁相变。

在（A）[math]\displaystyle{ T^{*}=5.5116 }[/math]、（B）[math]\displaystyle{ T^{*}=4.5116 }[/math]、（C）[math]\displaystyle{ T^{*}=3.5116 }[/math]时，三个最大本征微观态随时间的演化图

右图显示了临界点周围的本征微观态的演变，这与Wolff算法 中阐明的动力学相关。可以发现，[math]\displaystyle{ V_{tI} }[/math]随着时间[math]\displaystyle{ t }[/math]在[math]\displaystyle{ 0 }[/math]附近波动。

在地球系统中

气候问题是人类社会得以拥有稳定发展的最重要因素之一，显著的气候变迁可能直接或间接地对物种进化、种族迁徙、社会变革等产生影响。能否良好地分析与解释地球气候，确乎成为关乎人类福祉的重大命题。所有复杂系统都由无数不同的、相互作用的个体组成，而这些个体又可能受大量不同、随机的因素支配。这些复杂的因素甚至可能是混沌的，如地球气候系统，初值的微小差异可能带来后续的完全不同的演化进程。但是，复杂系统往往具有很好的普适性（Universal-Examples of Behavior-sity），这使得刻画系统中所有可能状态的集合成为可能。

地球系统科学起源于20世纪前期，经过近百年的发展，行星的观测运动、地球系统模型的健全，为综合性科学的发展提供了动力。对于地球系统的稳定性问题，尤其是气候问题，是涉及整个地球系统的非线性、突变与临界点的相互作用的一大挑战。2021年诺贝尔物理学奖被颁发给真锅淑郎（Syukuro Manabe）和克劳斯·哈塞尔曼（Klaus Hasselmann），他们对复杂地球物理系统进行量化气候变化的物理建模，并且给出了对地球系统整体的可靠的趋势预测。

在某个时间点[math]\displaystyle{ t }[/math]上，第[math]\displaystyle{ i }[/math]个格点的地球表面气温（SAT）是[math]\displaystyle{ T_i（t） }[/math]。在一个时间段[math]\displaystyle{ M }[/math]中，网格[math]\displaystyle{ i }[/math]的平均SAT可以被计算为[math]\displaystyle{ \left\langle T_i\right\rangle=\frac{1}{M} \sum_{t=1}^M\gt T_i(t) }[/math]，在第[math]\displaystyle{ i }[/math]个格点处的SAT的波动是[math]\displaystyle{ \delta T_i(t)=T_i(t)-\left\langle T_i\right\rangle, t=1,2, \cdots, M }[/math]，均方根偏差为[math]\displaystyle{ \Delta_i=\sqrt{\frac{1}{M} \sum_{t=1}^M \delta T_i(t)^2} }[/math]。那么，用下式可以很好地描述SAT的波动特征： [math]\displaystyle{ \delta S_i(t)=\delta T_i(t) / \Delta_i }[/math]。

我们现在引入地球的微观态为 [math]\displaystyle{ \delta \boldsymbol{S}(t)=\left[\begin{array}{c} \delta S_1(t) \\ \delta S_2(t) \\ \vdots \\ \delta S_N(t) \end{array}\right] }[/math]，使用从1950年1月1日到2018年12月31日的数据集，我们可以得到一个[math]\displaystyle{ N×M }[/math]的集合矩阵[math]\displaystyle{ A }[/math]，元素为[math]\displaystyle{ A_{i t}=\frac{\delta S_i(t)}{\sqrt{C_0}} }[/math]，其中[math]\displaystyle{ C_0 = M · N }[/math]。

参考资料

https://mp.weixin.qq.com/s/ehGrm1FPlkq7b6KATmIhTw