什么是均方位移MSD？从定义到误差来源全解析

说明：本文华算科技主要介绍均方位移（MSD）的统计定义、时间标度、方向与集体分解，以及有限轨迹中的坐标和误差来源。

MSD 是什么，怎样由轨迹构造？

分子动力学轨迹记录原子或分子在离散时刻的位置。某个粒子 i 在起始时刻 t₀ 与稍后时刻 t₀+τ 之间的位移为 Δr_i=r_i(t₀+τ)−r_i(t₀)。MSD 统计的是位移矢量模长的平方，不是轨迹长度，也不是相邻帧步长的简单累加。

对含 N 个等价粒子的体系，常用定义写作 MSD(τ)=(1/N)Σ_i⟨|r_i(t₀+τ)−r_i(t₀)|²⟩_t0。式中的 滞后时间 τ 是两帧之间的时间间隔，横坐标并非整条模拟已经运行的时长。每一个 τ 都会对应一组起止帧配对。

一条长度为 T 的轨迹可为较短的 τ 提供许多起始时刻。τ 接近 T 时，可用配对快速减少。多个时间起点与多个等价粒子的平均降低了单次随机位移造成的波动，但重叠时间片共享大量坐标，所得样本仍含相关。

平方运算消除了方向符号：向 +x 与 −x 移动相同距离，会给出相同的平方位移。它还会提高长位移事件的权重。MSD 是位移分布的二阶矩；两种轨迹即使拥有相同 MSD，也可能分别由连续小步或少量长跳组成，运动机制仍要由位移分布、van Hove 函数和原子轨迹补充。

时间标度怎样区分振动/困笼/扩散？

许多运动区间可近似写成 MSD(τ)∝τ^α。双对数坐标中的局部标度指数 α(τ)=d ln MSD/d ln τ记录曲线在不同时间尺度上的形状。它与线性坐标中用于求扩散系数的斜率承担不同任务：前者识别幂律状态，后者在正常扩散区间给出 D。

惯性尚未耗散的极短时间内，粒子近似保持速度，α≈2；局部热振动或邻近粒子形成的困笼会使曲线弯折，平台附近 α 接近 0；无偏、长程随机游走进入 Fick 扩散后，α≈1。0αα≤2 对应超扩散。定向漂移在长时间处也可产生接近 2 的指数；它来自持续平均速度，短时弹道区则来自惯性记忆，两者对应的速度相关函数不同。

平台有多种物理来源。晶体中的原子可在势阱附近振动；玻璃中的粒子可在邻居形成的困笼内往复运动；孔隙或界面法向位移受几何范围限制；轨迹过短时，偶发跳跃尚未进入统计窗口。平台高度给出受限运动的长度尺度，而平台后的再次上升才对应逃离局部环境后的长程迁移。

图2中的胶体由稀疏表面向致密内部过渡。靠近表面的曲线接近线性，内部曲线逐渐出现亚扩散与平台；自中间散射函数 F_s(q,t) 的衰减也随深度变慢。MSD 与相关函数共同限定困笼持续时间，避免把低斜率笼统写成扩散系数减小。

同一深度仍可同时存在高位移与低位移粒子。图3用颜色编码单粒子位移，并标出相邻高位移粒子形成的团簇。总体 MSD 把它们压缩成一个平均数，粒子分辨的平方位移分布则保留迁移非均一和协同重排的空间信息。

粒子/方向/集体位移分别对应什么？

按元素求平均时，Li⁺ 的 MSD 回答“一个 Li⁺ 平均离开自身初始位置多远”；逐离子曲线回答哪些 Li⁺ 参与了长程迁移；按晶格位点或界面区域分组，则比较不同局部环境的停留与跳跃。粒子集合的定义直接改变平均值，跨组交换的原子还会使固定空间区域的统计出现成员变化。

三轴分量满足 MSD_3D=MSD_x+MSD_y+MSD_z。层状导体、晶界和纳米孔中的运动常有明显方向差异，三轴总和会隐藏某一方向的受限状态。分量 MSD 的斜率构成扩散张量，坐标轴应与晶向、界面法向或孔轴对应。

α-Li₃N 晶界中的 Li⁺ 位移主要沿 xy 平面增长，z 分量接近平台。若按三维 Einstein 关系直接除以 6τ，会把面内迁移稀释到三个方向。正常扩散区间的表达式为 MSD=2dDτ，其中自由运动方向数 d应与实际统计分量一致：一轴取 1，平面取 2，三轴取 3。

单粒子 MSD 给出示踪扩散 D^*。电流却由全部带电粒子的净位移产生，对同号离子可写成 ⟨|Σ_iΔr_i|²⟩；展开平方后，除各粒子的自项外，还出现 Δr_i·Δr_j。交叉项记录离子运动相关，同向协同与反向补偿会改变电荷扩散 D_σ。

图5中 Na⁺ 自项较接近线性，Na⁺–Cl⁻ 交叉项在不同独立体系间分散很宽，平均后呈负值。由 D^* 换算电导率时忽略该项，相当于采用无相关近似。Haven 比常写作 H_R=D^*/D_σ；部分文献使用它的倒数，数值比较时应连同定义式一并记录。

有限轨迹会把哪些误差写进 MSD？

周期性模拟盒把越过盒面的粒子映射回主盒后，坐标序列会在盒面处形成接近盒长的突变，连续迁移也会被折回有限范围。MSD 采用跨盒展开坐标，相邻帧间还应满足位移小于半个盒长，避免高速粒子在输出间隔内跨越多个周期像。

NPT 轨迹还含盒长涨落。图6的 TOR 方案累积主盒中的最小位移，保持普通扩散轨迹的统计特征；LAT 方案追踪涨落晶格中的周期像，远离原点后会放大盒缩放噪声。坐标展开算法属于 MSD 定义的一部分。长时间 NPT 轨迹含盒缩放项，沿用 NVT 坐标处理会遗漏该项。

体系整体漂移会在 MSD 中产生近似线性或二次增长，去除质心运动时，参照原子集合应固定为不参与目标扩散的结构骨架。研究固态电解质中的 Li⁺ 迁移时，可用承载骨架定义整体平移；若把可迁移 Li⁺ 自身纳入漂移参照，它们的净位移会被部分扣除。漂移修正对象与扩散对象使用不同原子集合，两者重合会改变物理量。

长 τ 处的起始时刻较少，方差随之增大；相邻 MSD 点又共享轨迹片段，误差彼此相关且大小不同。普通最小二乘把各点视作同方差独立数据，回归公式给出的斜率误差常偏小。分块轨迹、独立速度种子或带协方差的回归可把扩散系数写成分布或置信区间，而非没有误差条的单值。

图7a汇总了 192 个 LLZO 有效模拟得到的 D^* 点估计，红色直方图给出运行间分布，灰线来自完整数据的贝叶斯回归，粉色横线标出标准差。图7b中的蓝色直方图记录单次有效模拟给出的方差估计，粉色竖线是 192 次结果的样本方差。单次轨迹的误差估计本身也有分布，随意延长到噪声占主导的曲线末端，并不会自动增加扩散信息。

高温 AIMD 常用来增加跳跃数，但温度升高可能改变晶相、缺陷浓度或局部配位。若承载骨架已经熔化，所得线性 MSD 描述液态扩散；若 20 ps 内仅出现一次 Li⁺ 跳跃，斜率主要受单个事件控制。相结构、有效跳跃数与分块斜率分布限定该段 MSD 可对应的材料状态。