统计力学是物理学的一个分支，将粒子的微观行为与物质的宏观性质联系起来。它为热力学提供微观基础，通过力学定律和概率解释压强、热容和熵等物理量。统计力学并不逐一追踪每个组分，而是研究各种可能微观状态的分布。其方法可用于描述气体、液体、固体、辐射，以及许多包含大量相互作用自由度的其他系统。(damtp.cam.ac.uk)
历史发展
统计力学起源于19世纪通过分子运动解释热和气体行为的尝试。詹姆斯·克拉克·麦克斯韦建立了分子速度的统计描述，路德维希·玻尔兹曼则发展了气体动理论，并将熵与微观构型的数目联系起来。1877年，玻尔兹曼提出了熵的这一统计解释。随后，约西亚·威拉德·吉布斯以系综为核心构建了统计力学体系；他于1902年出版的《统计力学基本原理》为平衡理论提供了系统框架。此后，量子统计将这些方法推广到微观状态遵循量子规律的系统。(damtp.cam.ac.uk)
微观状态与宏观描述
微观状态规定了系统的微观状况。在经典力学中，这通常指系统中各粒子的位置和动量，可用相空间中的一个点表示。在量子力学中，状态和可观测量需要用量子理论描述，统计预测则由概率或密度算符得出。宏观状态由少数可测量的物理量确定，例如体积、温度和总能量。许多不同的微观状态可以对应同一个宏观状态。(damtp.cam.ac.uk)
当各个可及微观状态等概率出现时，熵可由玻尔兹曼关系式表示：
[ S=k_{\mathrm B}\ln\Omega, ]
其中，(\Omega) 为这些状态的数目，(k_{\mathrm B}) 为玻尔兹曼常数。对于各状态概率为 (p_i) 的离散分布，吉布斯表达式为
[ S=-k_{\mathrm B}\sum_i p_i\ln p_i. ]
当每个可及状态的概率均为 (1/\Omega) 时，这一表达式就化为玻尔兹曼表达式。因此，熵具有精确的统计意义，而不只是对“无序程度”的定性度量。(ocw.mit.edu)
平衡系综
统计系综是满足指定宏观约束的系统的假想副本所构成的集合，并为这些副本赋予概率分布。其中有三种系综尤为重要：
- **微正则系综：**描述能量、体积和粒子数均固定的孤立系统。在其标准平衡表述中，指定能量区间内的各个可及状态被赋予相同的概率。
- **正则系综：**描述与热库交换能量、但体积和粒子数保持不变的系统。
- **巨正则系综：**描述在指定温度、体积和化学势下，既可与库交换能量，也可交换粒子的系统。(damtp.cam.ac.uk)
这些系综采用不同的概率分布，但在热力学极限下，往往会给出相同的宏观热力学预测；所谓热力学极限，是指在密度固定的条件下，粒子数和体积趋于无穷大。这种等价性需要满足适当条件，并非普遍成立，对于具有长程相互作用的系统尤其如此。(damtp.cam.ac.uk)
配分函数与可观测量
在正则系综中，能量为 (E_i) 的微观状态具有概率
[ p_i=\frac{e^{-\beta E_i}}{Z}, \qquad \beta=\frac{1}{k_{\mathrm B}T}, \qquad Z=\sum_i e^{-\beta E_i}. ]
配分函数 (Z) 使概率归一化，并包含系统的平衡热力学信息。若可观测量 (A) 在各微观状态中的取值为 (A_i)，其期望值为 (\langle A\rangle=\sum_i p_iA_i)。在经典表述中，对状态的求和由经过适当归一化的相空间积分替代。(damtp.cam.ac.uk)
对于不随温度变化的能谱，常用关系式包括
[ U=-\frac{\partial\ln Z}{\partial\beta}, \qquad F=-k_{\mathrm B}T\ln Z, ]
其中，(U) 为平均内能，(F) 为亥姆霍兹自由能。能量涨落将微观统计与可测量的响应联系起来：
[ \operatorname{Var}(E)=k_{\mathrm B}T^2C_V. ]
因此，能量的方差决定了定容热容 (C_V)。(damtp.cam.ac.uk)
经典统计与量子统计
对于相互作用可忽略的稀薄经典气体，统计力学给出状态方程 (PV=Nk_{\mathrm B}T)，也能推导出麦克斯韦—玻尔兹曼速度分布。粒子间的作用力可以解释气体行为对理想气体的偏离，而经典能量均分定理指出，能量中的每个独立二次项所对应的平均能量为 (k_{\mathrm B}T/2)。量子效应限制了这一结论的适用范围。(damtp.cam.ac.uk)
全同量子粒子需要采用不同的计数规则。玻色子遵循玻色—爱因斯坦统计，允许多个粒子占据同一个单粒子态。费米子遵循费米—狄拉克统计和泡利不相容原理。这些区别解释了电子简并、黑体辐射以及玻色—爱因斯坦凝聚态的形成等现象。在稀薄、非简并的极限下，这两种量子分布都趋近于经典的麦克斯韦—玻尔兹曼分布。(damtp.cam.ac.uk)
集体行为与相变
相互作用可以产生单个粒子所不具备的集体性质。统计力学通过考察平衡状态及其自由能，描述凝结和磁有序化等相变。伊辛模型由相互作用的离散自旋组成，为研究磁有序和临界行为提供了一个便于分析的模型。序参量用于区分不同的相；例如，磁化强度可用于表征铁磁有序。(damtp.cam.ac.uk)
在连续相变附近，涨落所涉及的空间尺度越来越大。与重整化密切相关的重整化群方法，解释了微观细节不同的系统为何能够具有相同的临界指数和标度行为。这种普适性使相变可以划分为不同类别，各类别由维数、对称性和相互作用范围等特征决定。(ocw.mit.edu)
非平衡理论与计算方法
非平衡统计力学研究弛豫、输运以及受外界驱动的系统。布朗运动和扩散将微观涨落与可观测的运动和扩散过程联系起来。动力学方程、随机过程和响应理论描述了随时间变化的行为中的一些方面，而这些并不能仅由平衡系综确定。(ocw.mit.edu)
精确计算通常只适用于相对简单的模型。分子动力学通过数值计算追踪粒子运动，而马尔可夫链蒙特卡洛则无需枚举所有状态，就能对统计分布进行采样。当对状态的求和或积分难以解析求解时，这些相互补充的方法可用于研究相互作用粒子、材料性质和相行为。(arxiv.org)