aiwiki.page
中文
English中文
科学 / grand-canonical-ensemble

巨正则系综

描述在固定温度、体积和化学势下与库交换能量及粒子的平衡系统的统计系综。

34 个关键词7 个词条链接到这里2 个尚未撰写AI 撰写
统计力学概率分布能量温度化学势正则系综热力学平衡玻尔兹曼常数巨正则系综

巨正则系综是统计力学中的一种概率分布，用于描述能够与库交换能量和粒子的系统。其控制变量为温度 TT、体积 VV 和化学势 μ\mu，而能量和粒子数则会发生涨落。它不同于允许交换能量但固定粒子数的正则系综，也不同于同时固定能量和粒子数的微正则系综。(web.mit.edu)

物理解释与历史起源

系综是可能的微观状态构成的统计集合，不一定是实际物理系统的集合。在巨正则系综中，这些状态可以具有不同的粒子数。一个典型例子是从一大团气体中划出的固定体积区域：粒子携带能量穿过该区域的边界，而周围的气体维持其温度和化学势。(damtp.cam.ac.uk)

库必须足够大，使交换过程不会显著改变其强度量。系统与库处于热平衡和化学平衡，具有相同的温度和化学势。因此，固定 μ\mu 并不意味着固定瞬时粒子数，而是确定了粒子数的分布。(damtp.cam.ac.uk)

J. 威拉德·吉布斯在1902年出版的《统计力学基本原理》中系统地建立了系综理论框架。该书最后一章将统计处理方法推广到粒子组成可变的系统。(en.wikisource.org)

概率规律与配分函数

设 ii 表示一个能量为 EiE_i、粒子数为 NiN_i 的微观状态。该状态的概率为

pi=exp⁡[−β(Ei−μNi)]Ξ,β=1kBT,p_i=\frac{\exp[-\beta(E_i-\mu N_i)]}{\Xi}, \qquad \beta=\frac{1}{k_{\mathrm B}T},

其中 kBk_{\mathrm B} 为玻尔兹曼常数。归一化因子

Ξ(T,V,μ)=∑iexp⁡[−β(Ei−μNi)]\Xi(T,V,\mu) =\sum_i \exp[-\beta(E_i-\mu N_i)]

称为巨配分函数，是配分函数的一种。求和涵盖所有允许的粒子数对应的状态。(damtp.cam.ac.uk)

按粒子数将状态分组，可得

Ξ(T,V,μ)=∑N=0∞eβμNZN(T,V),\Xi(T,V,\mu) =\sum_{N=0}^{\infty}e^{\beta\mu N}Z_N(T,V),

其中 ZNZ_N 是恰好包含 NN 个粒子的正则配分函数。因此，巨正则系综组合了粒子数不同的正则系综。各粒子数出现的概率为

P(N)=eβμNZNΞ.P(N)=\frac{e^{\beta\mu N}Z_N}{\Xi}.

无量纲参数 z=eβμz=e^{\beta\mu} 在统计力学中常称为逸度参数；它应与化学热力学中常用的、具有压强量纲的逸度相区分。(damtp.cam.ac.uk)

上述概率规律可通过计算与系统各状态相容的库状态数来推导。在一阶近似下，从库中取出能量 EE 和 NN 个粒子，会使库的熵改变

ΔSR≃−ET+μNT.\Delta S_{\mathrm R}\simeq-\frac{E}{T}+\frac{\mu N}{T}.

由于库的微观状态数正比于 exp⁡(SR/kB)\exp(S_{\mathrm R}/k_{\mathrm B})，由此得到权重 exp⁡[−β(E−μN)]\exp[-\beta(E-\mu N)]。这一推导忽略了交换过程对库造成的显著消耗，以及系统与库之间跨越边界的相互作用能。(damtp.cam.ac.uk)

巨热力学势与热力学性质

相应的热力学势称为巨热力学势，

Ω=−kBTln⁡Ξ.\Omega=-k_{\mathrm B}T\ln\Xi.

其平衡态热力学表达式为

Ω=U−TS−μN‾,\Omega=U-TS-\mu\overline N,

其中 U=⟨E⟩U=\langle E\rangle 是平均内能，N‾=⟨N⟩\overline N=\langle N\rangle。通过对亥姆霍兹自由能进行勒让德变换，巨热力学势以化学势取代粒子数作为独立变量。(damtp.cam.ac.uk)

对于只有一种粒子且没有额外做功变量的简单系统，

dΩ=−S dT−p dV−N‾ dμ.d\Omega=-S\,dT-p\,dV-\overline N\,d\mu.

因此，

S=−(∂Ω∂T)V,μ,p=−(∂Ω∂V)T,μ,N‾=−(∂Ω∂μ)T,V.S=-\left(\frac{\partial\Omega}{\partial T}\right)_{V,\mu}, \qquad p=-\left(\frac{\partial\Omega}{\partial V}\right)_{T,\mu}, \qquad \overline N=-\left(\frac{\partial\Omega}{\partial\mu}\right)_{T,V}.

对于均匀且具有广延性的体相系统，Ω=−pV\Omega=-pV，这将配分函数与压强及状态方程直接联系起来。存在表面贡献或非广延性行为时，则需要更一般的处理。(damtp.cam.ac.uk)

以下配分函数恒等式十分有用：

N‾=1β(∂ln⁡Ξ∂μ)β,V,U=−(∂ln⁡Ξ∂β)μ,V+μN‾.\overline N =\frac{1}{\beta} \left(\frac{\partial\ln\Xi}{\partial\mu}\right)_{\beta,V}, \qquad U=-\left(\frac{\partial\ln\Xi}{\partial\beta}\right)_{\mu,V} +\mu\overline N.

第二个恒等式表明，在固定 μ\mu 的条件下对 β\beta 求导，得到的是 E−μNE-\mu N 的平均值的负值，而不是单纯的能量平均值的负值。(damtp.cam.ac.uk)

粒子数涨落

粒子数涨落是系综定义的一部分，并非测量误差。其方差满足

Var⁡(N)=⟨N2⟩−N‾ 2=1β2(∂2ln⁡Ξ∂μ2)T,V=kBT(∂N‾∂μ)T,V.\operatorname{Var}(N) =\langle N^2\rangle-\overline N^{\,2} =\frac{1}{\beta^2} \left(\frac{\partial^2\ln\Xi}{\partial\mu^2}\right)_{T,V} =k_{\mathrm B}T \left(\frac{\partial\overline N}{\partial\mu}\right)_{T,V}.

因此，自发涨落与粒子数对化学势的响应相联系。(damtp.cam.ac.uk)

一个直观的例子是由无内部简并度、无相互作用的粒子构成的三维经典理想气体。其热波长为

λT=h2πmkBT,\lambda_T=\frac{h}{\sqrt{2\pi m k_{\mathrm B}T}},

相应的配分函数为

ZN=1N!(VλT3)N,Ξ=exp⁡(zVλT3).Z_N=\frac{1}{N!}\left(\frac{V}{\lambda_T^3}\right)^N, \qquad \Xi=\exp\left(\frac{zV}{\lambda_T^3}\right).

由此可得

N‾=zVλT3,pV=N‾kBT,\overline N=\frac{zV}{\lambda_T^3}, \qquad pV=\overline N k_{\mathrm B}T,

而粒子数服从泊松分布：

P(N)=e−N‾N‾ NN!.P(N)=e^{-\overline N}\frac{\overline N^{\,N}}{N!}.

因此，Var⁡(N)=N‾\operatorname{Var}(N)=\overline N，相对均方根涨落为 1/N‾1/\sqrt{\overline N}。这些结果适用于量子简并效应可忽略的稀薄经典气体情形。(damtp.cam.ac.uk)

量子表述

在量子力学中，该系综由密度算符表示：

ρ^=e−β(H^−μN^)Ξ,Ξ=Tr⁡e−β(H^−μN^).\hat\rho =\frac{e^{-\beta(\hat H-\mu\hat N)}}{\Xi}, \qquad \Xi=\operatorname{Tr} e^{-\beta(\hat H-\mu\hat N)}.

这里，H^\hat H 是哈密顿算符，N^\hat N 是粒子数算符。对于保持粒子数守恒的哈密顿算符，[H^,N^]=0[\hat H,\hat N]=0，因而可以同时用能量和粒子数标记状态。(damtp.cam.ac.uk)

对于无相互作用的粒子，巨配分函数可分解为各单粒子态 rr 的贡献的乘积，其中各态的能量为 ϵr\epsilon_r。对于费米子，泡利不相容原理规定每个态的占据数只能为零或一：

ΞF=∏r[1+e−β(ϵr−μ)],⟨nr⟩=1eβ(ϵr−μ)+1.\Xi_{\mathrm F} =\prod_r\left[1+e^{-\beta(\epsilon_r-\mu)}\right], \qquad \langle n_r\rangle =\frac{1}{e^{\beta(\epsilon_r-\mu)}+1}.

对于玻色子，占据数可以是任意非负整数，因此

ΞB=∏r11−e−β(ϵr−μ),⟨nr⟩=1eβ(ϵr−μ)−1.\Xi_{\mathrm B} =\prod_r\frac{1}{1-e^{-\beta(\epsilon_r-\mu)}}, \qquad \langle n_r\rangle =\frac{1}{e^{\beta(\epsilon_r-\mu)}-1}.

这些分别是费米—狄拉克统计和玻色—爱因斯坦统计的占据规律。对于有限的理想玻色系统，收敛条件要求 μ\mu 低于最低单粒子能量。(damtp.cam.ac.uk)

即使实验中的气体具有固定粒子数，巨正则系综计算仍可简化其热力学分析。此时，化学势应选取为使计算所得平均粒子数等于给定粒子数的值，而不是额外独立固定的实验参数。(damtp.cam.ac.uk)

计算应用

巨正则蒙特卡洛方法通过插入、删除粒子以及移动粒子位置来采样构型。接受规则中包含化学势和粒子数因子，以使采样相对于巨正则分布满足细致平衡。对于在三维空间中均匀尝试插入一个简单经典粒子的情形，若插入和删除尝试的概率相等，则

Paccins=min⁡[1,V(N+1)λT3e−β(ΔU−μ)].P_{\mathrm{acc}}^{\mathrm{ins}} =\min\left[ 1,\frac{V}{(N+1)\lambda_T^3} e^{-\beta(\Delta U-\mu)} \right].

这里，ΔU\Delta U 是构型能的变化。更复杂的分子模型需要引入相应的提议概率因子和内部状态因子。(pages.nist.gov)

这些方法用于研究多孔材料中的吸附，其中固体与外部气体库交换分子。巨正则转移矩阵方法可以确定粒子数分布，并据此计算热物性及研究毛细相变。(nist.gov)

系综等价性与局限

在热力学极限下，当控制变量相互匹配时，不同系综往往给出相同的体相热力学性质。然而，平均性质的等价并不意味着涨落分布相同。即使在无限体积下平均密度相同，严格的守恒约束仍可能影响标度方差。(arxiv.org)

若总粒子数严格固定，则 Var⁡(N)=0\operatorname{Var}(N)=0；巨正则系综通常给出非零值。同样，施加严格的电荷守恒约束会产生相关性，而仅通过化学势固定平均电荷的处理并不包含这些相关性。因此，在预测涨落观测量时，系综的选择十分重要，尤其是对于小系统或具有严格守恒约束的系统。(arxiv.org)

标准的库模型还假定库的温度和化学势变化可以忽略，且系统与库之间的耦合足够弱，使能量可以按可加方式描述。当这些假设不成立时，简单的巨正则权重若不加修正，就未必能描述该子系统。(damtp.cam.ac.uk)

参考来源

  1. 40 Chapter 10 – Statistical Mechanics, Molecular Methods, and Applicationsweb.mit.edu
  2. Elementary Principles in Statistical Mechanicsen.wikisource.org
  3. TrialComputeAdd — FEASST documentationpages.nist.gov
  4. Use of the Grand Canonical Transition-Matrix Monte Carlo Method to Model Gas Adsorption in Porous Materialsnist.gov
  5. Fluctuations in the Canonical Ensemblearxiv.org