多尺度分子动力学模拟

很多材料问题横跨好几个数量级的时间和空间：电子层面决定键的属性，原子层面决定局部运动，介观层面决定相行为，宏观层面决定构件响应。单一方法根本覆盖不了。我这个项目客户的高分子复合材料要算力学松弛，从单体到宏观构件全尺度都要管，于是我搭了一套 多尺度分子动力学模拟 流程，把 QM → 全原子 MD → 粗粒化 MD → 连续介质 FEM 串起来。

尺度衔接从最底层开始。第一层用 DFT 算单体（重复单元）的电荷分布、二面角旋转势垒，把得到的量子精度信息转成全原子力场的二面角参数。这一步保证了原子层力场在”电子层面”是自洽的，而不是随便套个通用力场。第二层全原子 MD 跑 100 ns，输出径向分布函数（RDF）、密度、扩散系数等结构量，作为下一层的参考。第三层用映射法把若干个原子粗粒化成（bead）一个 CG 珠子，用 IBI 或力匹配拟合 CG 势，跑 10 μs 介观模拟，覆盖全原子够不到的时间尺度。第四层把 CG 得到的应力-应变关系作为本构模型，喂给宏观 FEM 算构件变形。

这个流程里最关键、也最容易被做成”概念游戏”的，是层与层之间的传递量。我给每一层都定义了明确的输出和输入：DFT → 力场参数；全原子 → RDF/密度/扩散系数；CG → 黏弹本构；FEM → 宏观响应。没有这些硬传递量，多尺度就是四张互不相关的图贴在一起。

验证环节用 RDF 一致性。右图对比了全原子 MD 和粗粒化 MD 的径向分布函数，第一峰位置和高度吻合得很好，说明 CG 力场的结构描述是可靠的。如果 RDF 对不上，说明映射方案或 CG 势有问题，必须回退到全原子数据重新拟合——我这个项目里第一版 CG 势的 RDF 第二峰偏高，排查发现是 bead 内的键长约束没设对，修正后吻合。

时间尺度是另一个重点。全原子 MD 能稳稳覆盖 ns–μs，但很多高分子松弛要 ms 以上，必须靠 CG 才能摸到。我这个项目 CG 模拟跑了 10 μs，才看到应力松弛的完整衰减曲线；如果只跑全原子，连第一个松弛平台都够不到。

踩坑点：第一，力场传递的保真度。DFT 算的二面角参数如果不做温度/相态适用性检查，直接用于液态 MD 会出问题，我做了液相构象抽样验证。第二，CG 映射非唯一。同一分子可以有多种 bead 方案，映射比太大丢化学细节（氢键、π-π 堆积），太小又失去加速意义，我这个项目用 5:1 映射兼顾。第三，跨尺度边界。FEM 用的本构若来自单一 CG 构型，可能漏掉多轴应力状态，我补做了多应变率 CG 扫描生成完整本构曲面。第四，验证不全。只验 RDF 不够，还需验扩散系数、玻璃化转变温度等，我列了 4 个一致性指标。

从我的工程经验看，多尺度分子动力学模拟 是”用对的方法算对的尺度”，而不是”算得越细越好”。我给客户的建议是：先明确要回答的工程问题在哪个尺度，再决定流程下探到哪一层，别为炫技做全 QM。更多跨尺度案例，见 [多尺度分子动力学模拟](https://www.keyanxueshu.com/category/md/)；需要搭建流程，[多尺度分子动力学模拟](https://www.keyanxueshu.com/) 上有服务入口。