说明：本文华算科技主要介绍玻璃化转变温度 Tg 在材料计算中的含义、常见取值方式、冷却轨迹中的图谱信号，以及用 Tg 比较非晶聚合物时应报告的模型条件与适用范围。
基本概念：Tg描述的是非晶材料的链段冻结
在非晶聚合物、无定形氧化物和金属玻璃中，Tg 对应结构松弛时间接近观测时间的温度区间。高于该区间，链段构象、自由体积和局域配位仍能在实验或模拟窗口内调整；低于该区间，原子或链段位移受到局域笼限制，密度、焓、模量和热膨胀系数的斜率会改变。
晶体熔化常伴随潜热和长程有序消失，玻璃化转变则发生在缺少长程周期结构的非平衡材料中。Tg 常写成一个温度值，实际更接近一段窄温区；冷却速率升高时，非晶链段在较高温度已来不及完成构象采样，计算得到的 Tg 常偏高。慢速 DSC、介电谱、黏度测试和 MD 冷却轨迹给出的数值可能相差几十 K。
图1. ULTEM 聚合物的密度-温度曲线和 525 K 交点。DOI：10.48550/arXiv.2001.08889。
ULTEM 的 ρ(T) 随 T 近似线性下降，低温玻璃态和高温橡胶态的斜率不同。两段拟合交点给出 525 K；该数值来自 1 atm NPT 冷却后的密度数据，随力场、冷却程序、平衡时长和拟合温区变化。若同一组模型采用相同程序，交点温度可用于比较聚合物结构改变后的相对耐热差别。
酰亚胺环、醚键、砜基、异构连接方式和端基会改变链段转动、堆积密度与自由体积分布。同一类聚醚酰亚胺也可能因为重复单元刚性不同而给出不同 Tg；计算文章若只给温度，读者无法判断数值来自链段刚化、密度变化还是拟合窗口选择。
图2. 五种聚醚酰亚胺的重复单元结构。DOI：10.48550/arXiv.2001.08889。
聚酰亚胺数据集中，Tg 主要分布在 500 到 580 K 附近，低温端和高温端仍有样本。材料家族内部的分散提醒我们，Tg 不宜写成某类聚合物的固定常数；重复单元、分子量、端基、交联、含水量和测量速率都会改变数值。
图3. 聚酰亚胺 Tg 数据的温度分布。DOI：10.48550/arXiv.2001.08889。
静态 DFT 给局域势能面、键合强弱、电子态和构型能量，宏观 Tg 还依赖长时间构象采样。AIMD 的 ps 级轨迹适合小模型短时热响应和局域结构演化；聚合物 Tg 的主力数据仍来自经典 MD、粗粒化 MD、DSC、DMA 或介电松弛谱。
Tg在冷却轨迹中怎样取值
ρ(T)、比容 v(T)、焓 H(T) 和势能 U(T) 都可记录冷却过程中的热响应。常用处理是把低温玻璃态数据和高温熔体数据分别线性拟合，交点对应模拟观测时间下的 Tg。拟合点数量、温度间隔、是否剔除未平衡数据，会直接改变交点位置。
温度点太少时，两段直线可能只是拟合算法的产物；温度点过密但每个点平衡不足，曲线会带有冷却历史。用于材料比较的数据表应列出每个温度的平均密度、误差条和采样时长，并给出低温段与高温段的拟合范围。
势能曲线给出的是热运动和构象采样的能量响应。受限聚合物薄膜的 Utot(T) 在低温和高温区各有近似线性段，箭头标出的交点位于约 0.64 到 0.67 ε/kB。同一薄膜模型内，密度交点和势能交点可互相校验；若两个交点相距较远，平衡时长、样品厚度和拟合区间应重新核对。
图4. 受限聚合物薄膜的 Utot(T) 两段拟合交点。DOI：10.48550/arXiv.2306.01560。
动力学信号给的是另一类观测时间。小分子在聚合物自由体积中的扩散系数 D 可由 MSD 长时间斜率得到，D(T) 在 Arrhenius 坐标中若出现斜率改变，说明扩散活化能在玻璃态和橡胶态之间发生变化。ULTEM 中氩扩散的 D(T) 交点约为 504 K，低于密度法得到的 525 K，两个数值反映的采样对象不同。
图5. 氩在 ULTEM 中的扩散系数随 1/T 的变化。DOI：10.48550/arXiv.2001.08889。
密度法关心体积热膨胀，扩散法关心自由体积连通与小分子迁移，势能法关心构象采样后的平均能量。三类图谱给出的 Tg 可相近，也可有差别。科研写作中宜写明物理量名称、统计窗口和拟合温区，否则读者无法区分“材料真的更耐热”和“取值方式变了”。
冷却速率、链长和薄膜厚度为什么会改变Tg
MD 的时间步通常是 fs，轨迹长度常在 ns 到 µs 范围，冷却速率远高于常规 DSC。快速冷却会把高温构象冻结到较高温区，Tg 交点随速率升高而上移。若要比较两个聚合物结构，冷却速率、退火时间、温度步长和压力控制方式应保持一致。聚合物 Tg 对非键相互作用很敏感，范德华参数会改变链间堆积，部分电荷会改变酰亚胺、砜基和醚氧附近的偶极相互作用，二面角势控制链段旋转能垒。
短链和长链在端基占比、缠结程度和构象约束上不同，Tg(0)(N) 随 N 增大而趋近 Tg∞。Fox-Flory 形式 Tg(N)=Tg∞-K/N 常用于描述体相聚合物的链长依赖；静态 DFT 或半经验量子化学可作为单体电荷和构型能量的参照，长链冷却曲线仍依赖可采样的大模型。
图6. 体相聚合物 Tg(0)(N) 随链长 N 的变化。DOI：10.48550/arXiv.2306.01560。
对于 h 接近或小于 Rg(0) 的自由薄膜和受限薄膜，表面或墙面附近链段的垂直构象受到约束，密度曲线的斜率和交点都会变化。Tg(N,h)/Tg(0)(N) 随 h 降低而下降，薄膜厚度改变的是链段可采样空间，薄膜数值不宜直接代替体相材料数值。
图7. 受限与自由薄膜中的密度曲线和相对 Tg。DOI：10.48550/arXiv.2306.01560。
模型尺度带来的偏差并非只来自原子数。端基处理、周期条件、墙面相互作用、孔隙、残余应力、含水量和交联程度都会改变非晶结构的松弛时间。同一材料比较 Tg 时，N、h、密度构建方式和退火程序应随数值一起给出。
Tg判断应核对哪些结构和动力学信号
密度曲线给热膨胀斜率，势能曲线给平均构象能量，D(T) 给扩散活化行为，RDF 给局域配位，二面角分布给链段旋转状态，MSD 或中间散射函数给松弛时间。Tg 的可信程度来自这些信号是否处在相近温区，而不是一个交点温度本身。
链构象数据把薄膜约束写成可量化变量：Re,|| 与 Rg,|| 在较厚薄膜中接近体相，Re,⊥ 与 Rg,⊥ 随 h/Re(0) 降低而明显下降；取向参数 S 和 kink 数目随厚度改变。这些结构量能解释薄膜 Tg 下降的来源，也能区分表面自由度变化与链长效应。
图8. 聚合物薄膜中链构象和取向参数随厚度的变化。DOI：10.48550/arXiv.2306.01560。
对非晶氧化物和金属玻璃，Tg 不由链段旋转控制，局域配位数、键角分布、结构因子 S(q)、MSD 平台和黏度外推常作为取值依据。不同材料类型应选择不同结构变量：聚合物侧重二面角和自由体积，氧化物侧重配位网络重排，金属玻璃侧重短程序和剪切模量变化。
实验 Tg 与模拟 Tg 对比时，报告项宜包含冷却速率、温度步长、平衡时长、统计系综、压力、力场或 DFT 泛函、模型尺寸、密度构建方式、拟合温区和误差条。DSC 的 Cp 拐点、DMA 的 tan δ 峰、介电 αrelax 峰对应不同观测频率，和 MD 密度交点并列时应标注测量频率或扫描速率。
用于筛选耐热聚合物时，密度交点、Utot(T) 交点、D(T) 斜率变化和链构象约束宜处在相近温区。若数值分散，优先检查冷却速率、平衡时长、样品尺寸、端基、孔隙和薄膜界面设定；材料比较中优先采用同一套冷却程序下的相对排序，绝对 Tg 数值只随模型条件一起使用。