晶格畸变与晶格应变：参考结构、定量表征及第一性原理计算

说明：本文华算科技主要介绍晶格畸变与晶格应变的参考结构、数学变量、定量指标及第一性原理计算信号。

晶格应变与晶格畸变的基础概念怎样分开

一份结构文件同时包含晶格矢量和原子分数坐标。设参考晶胞矩阵为 H₀，目标晶胞矩阵为 H，两者满足 H=FH₀。变形梯度 F 记录晶胞整体的拉伸、压缩、剪切和转动，剔除刚体转动后得到应变张量。晶格应变由晶胞相对参考态的几何变化定义，单轴 2% 拉伸、双轴压缩和剪切加载都可写成 F 或 ε_ij 的分量。

原子实空间坐标还可写成 r_i=Fr_i⁰+u_i。第一项是参考原子随晶胞发生的仿射位移，第二项 u_i 是原子偏离仿射位置的内部位移。键长分裂、键角偏转、配位多面体拉长、八面体转动、阳离子偏心以及缺陷近邻重构都写在 u_i 中。晶格畸变描述的是原子内部自由度偏离选定高对称构型的幅度与形式，它没有与应变张量完全相同的统一标量。

缺陷、固溶体和相变模型必须为目标结构指定可追踪的参考态。缺陷超胞可采用同尺寸、同电荷设置的完整晶体作为 H₀ 和 r_i⁰；固溶体可把未弛豫的理想占位模型作为母相；相变研究则选高对称相并完成原子映射。不同化学组成、不同原子排序或不同超胞取向若缺少一一映射，u_i 便失去统一参照。周期镜像、整体平移和刚体转动应在比较前消除。

两类量的差别不宜简化成“应变属于整体、畸变属于局部”。GPA 可以给出纳米尺度的局域应变场；一个反铁畸变转动模也能在完整晶体中周期重复。区分依据是参考结构与结构变量：晶胞矩阵变化归入应变，仿射变换之后仍保留的内部坐标偏移归入畸变。材料文献有时用“lattice distortion”泛指所有偏离理想晶格的变化，计算报告需额外列出 ε_ij、u_i 或具体多面体指标。

应变张量与畸变指标分别怎样计算

应变由晶胞或局域基矢的相对变化给出

小形变下常用 ε=(F+F^T)/2−I；形变较大时可改用 Green–Lagrange 应变 E=(F^TF−I)/2。ε_xx、ε_yy、ε_zz 给出正应变，非对角项给出剪切分量。二维 slab 的真空厚度不属于材料厚度，面外应变应从真实原子层厚度、键长或有效厚度定义，避免把真空层变化计入晶体形变。

轴向百分比 ε=(a−a₀)/a₀ 适用于晶轴平行且剪切分量为零的比较。单斜、三斜晶胞或发生晶格转动的结构若逐项比较 a、b、c，基矢旋转和轴置换会混入标量。极分解 F=RU 可把刚体转动 R 与纯伸缩 U 分离，再由 U 构造有限应变。

原子级局域应变可在每个中心原子周围拟合 F_i，令参考近邻矢量经过 F_i 后尽量接近弛豫近邻矢量。最优仿射拟合的对称部分给出局域 ε_ij，拟合残差 D²_min 则记录无法由同一个线性变换吸收的非仿射运动。位错芯、剪切带、缺陷近邻和化学无序区域常出现较大的 D²_min，此数值更接近局域畸变或重排幅度。

周期 DFT 施加均匀应变时，参考态通常取零外压下完全弛豫的同一相。外延模型则以衬底约束后的面内晶格为控制量，面外晶格和内部坐标继续弛豫。实验局域应变图会选取一块晶格规整区域作为零点。参考态一旦改变，应变的符号和数值也随之改变，文章中应同步给出参考晶格、晶向、加载方式和是否允许晶胞弛豫。

畸变由内部坐标、键长、键角和对称模量化

配位多面体可用键长畸变指数 Δ=(1/n)Σ[(d_i−d̄)/d̄]²、二次伸长率 λ 和键角方差 σ²_θ 描述。Δ 关注同一多面体内的键长离散，σ²_θ 关注键角偏离理想值的程度。平均键长相同的两个八面体，仍可能具有完全不同的键长分布和转动模式；单独报告平均值会丢失局部几何。

复杂相变还可把弛豫结构投影到高对称母相的对称适配模，得到每个模的振幅 Q_α。八面体同相转动、反相转动、呼吸模、Jahn–Teller 模和极化位移由不同不可约表示标记。模振幅保留位移方向、相位关系和对称性降低方式，比“结构发生畸变”提供更多可复算数据。

应变为何会诱发畸变，畸变又何时产生自发应变

外加 ε 会改变键长、键角和配位空间，原子在固定晶胞中继续寻找受力为零的位置。“固定晶胞 + 弛豫内部坐标”把应变控制量和畸变响应放进同一次计算。钙钛矿薄膜受双轴应变后，面外晶格按泊松响应改变，氧八面体还可能通过转动、倾斜或极化偏心降低总能。应变是计算输入，弛豫得到的 Q_α、Δ 和 σ²_θ 是结构响应。

固定同一个 ε 时，高对称结构的力常数矩阵可能出现负本征值，对应声子谱中的虚频。沿该本征矢施加微小位移并弛豫，能量随模振幅形成双势阱，结构转入低对称极小值。虚频本征矢给出畸变方向，E(Q) 曲线给出高对称点到极小值的能量变化；应变只负责规定每条曲线所处的晶胞条件。

畸变也可能在零外加应变下自发出现。轨道简并解除、软声子凝聚、局域电荷排布或离子尺寸失配会推动内部位移；当位移模与 ε_ij 存在弹性耦合项，例如 gQ²ε，晶胞形状随 Q 的出现而改变，形成自发应变。因果顺序要由计算约束确定：固定应变扫描 Q 可研究应变诱导畸变，固定高对称晶胞释放软模可研究畸变驱动的晶胞响应。

两类量还可各自独立存在。单原子 Bravais 晶格的均匀仿射拉伸具有非零 ε，原胞内没有可弛豫的内部自由度；反铁畸变八面体转动可在平均晶格常数近乎不变时保持较大的 Q_R。高熵或缺陷超胞常保留接近立方的平均晶胞，位点周围的键长分布却明显展宽。宏观应变为零与局域畸变为零是两条独立陈述。

DFT 结果里怎样把两类结构变化分开报告

采用两组弛豫提取内部坐标贡献

同一 ε 下可保留两份结构：第一份把分数坐标固定在仿射位置，第二份固定晶胞并弛豫所有内部坐标。两者的能量差 ΔE_int、应力差 Δσ_ij 和原子位移 u_i 直接刻画内部弛豫贡献。clamped-ion 结果对应纯仿射响应，relaxed-ion 结果包含应变诱发的畸变。晶胞与原子一起全弛豫时，输出结构给出两类变量共同作用后的终态；因果分解依赖额外受限计算。

结构量与电子结构量分栏报告

结构部分应给出 H₀、H、ε_ij、残余应力、关键 u_i、键长分布、键角方差和主导 Q_α。电子部分再列带隙、PDOS、晶场分裂、磁矩或 COHP 的变化。能带移动和 DOS 峰位属于结构变化后的电子响应，它们不承担应变或畸变的几何定义。结构量写全后，电子信号才能对应到特定键、特定位移模和特定加载方向。

应变扫描宜逐点输出 E(ε)、σ_ij(ε)、Q_α(ε)、Δ(ε) 和关键键长。弹性区内 E(ε) 平滑且可作低阶拟合；Q_α 突跳、空间群改变或新虚频出现时，当前结构已跨入另一相，跨相数据不再用于同一组弹性常数拟合。有限温度 AIMD 还要分开时间平均晶胞与瞬时内部位移，热平均恢复高对称性时，单帧结构仍可保留较大的动态畸变。

高熵、缺陷、界面和非晶化前驱结构还应报告分布，而非把所有位点压缩成一个平均数。均值描述晶胞尺度，分布宽度描述位点离散，模振幅描述协同位移。对应计算文件至少保留参考结构、原子映射、周期镜像选择和刚体转动消除方式；否则 u_i 会混入平移、转动或错误的原子配对。一个可复算的结构结论应分别写出 ε_ij 与主导 Q_α，并列出承担畸变的键长、键角和原子编号。