字数 2266，阅读大约需 12 分钟

《Molecular Dynamics Simulation of Nanocomposites》

本书聚焦纳米复合材料分子动力学（MD）模拟，系统讲解Materials Studio、LAMMPS、GROMACS三大主流软件的使用方法、模拟原理与实战案例，覆盖金属基、聚合物基、陶瓷基纳米复合材料的性能预测与机理分析，是纳米材料模拟领域的实操指南。

第1章 分子动力学导论

1.1 核心模拟方法

• • 分子动力学（MD）：基于牛顿运动方程，模拟粒子随时间的运动，可研究平衡/非平衡物理现象，是本书核心方法。
• • 蒙特卡洛（MC）模拟：通过概率法则生成微观状态，仅适用于热力学平衡体系，无法体现时间维度的动态过程。
• • 布朗动力学（BD）：将溶剂视为连续介质，模拟分散粒子的布朗运动，忽略溶剂分子细节。
• • 耗散粒子动力学（DPD）：将溶剂分子簇视为虚拟流体粒子，适配大尺度、长时程模拟。
• • 格子玻尔兹曼法（LBM）：基于虚拟流体粒子的格子运动，侧重流体体系模拟。

1.2 基础核心概念

1. 1. 力场：描述粒子势能的数学函数，分键合（键伸缩、键角、二面角）与非键合（范德华力、静电作用）项，主流力场包括AMBER、CHARMM、COMPASS、CVFF。
2. 2. 势能函数：
• • Tersoff/Brenner势：适配碳基材料（碳纳米管、石墨烯）的多体键序势。
• • 莫尔斯势：描述化学键的非简谐振动，可模拟键断裂。
• • 伦纳德-琼斯势：经典范德华作用势，广泛用于非键合相互作用。
3. 3. 系综：
• • 微正则系综（NVE）：粒子数、能量、体积恒定；
• • 正则系综（NVT）：粒子数、体积、温度恒定；
• • 等温等压系综（NPT）：粒子数、压力、温度恒定；
• • 巨正则系综（μVT）：化学势、体积、温度恒定。
4. 4. 控温/控压方法：
• • 控温：Andersen、Berendsen、Nose-Hoover温控器；
• • 控压：Berendsen、Parrinello-Rahman控压器。
5. 5. 边界条件：
• • 周期性边界条件（PBC）：模拟无限大体系，消除边界效应；
• • Lees-Edwards边界：适配简单剪切流的非平衡模拟。

1.3 MD模拟流程

1. 1. 设定粒子初始位置与速度；
2. 2. 计算粒子间作用力；
3. 3. 采用Velocity Verlet/Leapfrog算法更新粒子位置与速度；
4. 4. 循环迭代，实现体系弛豫与动态模拟。

第2章 Materials Studio、LAMMPS、GROMACS概述

2.1 BIOVIA Materials Studio（商业可视化软件）

1. 1. 核心模块：
• • Materials Visualizer：建模、可视化、结构编辑；
• • Amorphous Cell：构建无定形聚合物、复合材料模型；
• • Forcite：经典力学计算、几何优化、MD模拟、力学性能预测；
• • COMPASS：凝聚相材料专用力场，适配聚合物、金属、陶瓷。
2. 2. 模拟策略：构建碳纳米管（CNT）→ 聚合单体建模 → 无定形晶胞包裹聚合物与增强相 → 几何优化 → 动力学弛豫 → 力学/热学性能计算。
3. 3. 案例：CNT/聚合物、CNT/金属基复合材料建模与性能模拟。

2.2 LAMMPS（开源并行MD软件）

1. 1. 核心特性：开源、并行计算，支持百万级粒子模拟，适配金属、聚合物、纳米复合材料。
2. 2. 模拟流程：
• • 定义体系构型与边界条件；
• • 编写输入脚本（初始化、原子定义、参数设置、运行指令）；
• • 体系弛豫（高温弛豫→降温至目标温度）；
• • 力学/结构性能计算（应力应变、玻璃化转变温度、应力松弛）。
3. 3. 输入/输出：输入脚本控制模拟参数，输出轨迹文件、日志文件，支持后处理分析。
4. 4. 案例：聚丙烯/CNT纳米复合材料的拉伸、应力松弛模拟。

2.3 GROMACS（开源生物/材料MD软件）

1. 1. 核心特性：高效计算、开源，适配生物分子、聚合物、界面体系模拟。
2. 2. 工作原理：构建拓扑结构→设定初始条件→计算作用力→积分运动方程→输出轨迹；支持周期性边界、控温/控压、表面张力耦合。
3. 3. 算法：Leapfrog、Velocity Verlet积分器，Berendsen/Nose-Hoover控温，支持非平衡模拟（热传导、剪切流）。
4. 4. 应用：生物芯片、生物分子建模、聚合物界面、离子传输模拟。
5. 5. 案例：磷脂膜、酶-配体结合、聚合物/电解质界面模拟。

第3章 金属基纳米复合材料的MD模拟

3.1 Materials Studio模拟石墨烯/SiC增强铝基复合材料

1. 1. 研究对象：SiC纳米颗粒、石墨烯增强铝（Al）基复合材料。
2. 2. 模拟方法：
• • 构建纯Al、SiC/Al、石墨烯/SiC/Al模型，采用EAM势（Al）、Tersoff势（SiC/石墨烯）、莫尔斯势（界面）；
• • NPT弛豫（300K）→ 单轴拉伸模拟，预测弹性模量、拉伸强度。
3. 3. 核心结论：
• • SiC体积分数提升，复合材料弹性模量、强度升高，塑性下降；
• • 小尺寸SiC颗粒增强效果更优（界面更多，载荷传递效率高）；
• • 石墨烯引入改善界面结合，进一步提升强度，降低塑性。

3.2 LAMMPS模拟石墨烯/铜纳米层状复合材料

1. 1. 研究对象：铜（Cu）-石墨烯交替层状复合材料。
2. 2. 模拟方法：
• • 构建不同厚度的Cu/石墨烯层状模型，采用EAM势（Cu）、Tersoff势（石墨烯）、LJ势（界面）；
• • NPT弛豫→ 单轴拉伸，分析应力应变、弹性模量、变形机理。
3. 3. 核心结论：
• • 石墨烯层显著提升Cu基复合材料强度、刚度，Cu层越薄增强效果越明显；
• • 石墨烯阻碍位错运动，抑制滑移，实现强化；
• • 层状结构提升材料韧性与失效应变。

3.3 GROMACS模拟锂金属/聚合物电解质界面

1. 1. 研究对象：聚氧化乙烯（PEO）-LiTFSI聚合物电解质与锂（Li）金属界面。
2. 2. 模拟方法：
• • 构建Li壁/Li slab与电解质界面模型，采用OPLS-AA力场、LJ势（Li）；
• • NPT弛豫，分析径向分布函数（RDF）、配位数、扩散系数。
3. 3. 核心结论：
• • 界面区域离子、聚合物链有序度更高，形成双电层；
• • 界面处离子扩散速率显著低于本体，界面阻力影响离子传输；
• • Li⁺与PEO氧原子强配位，界面处离子配对更显著。

第4章 聚合物基纳米复合材料的MD模拟

1. 1. 研究对象：碳纳米管（CNT）、功能化CNT增强聚合物（PMMA、PmPV、天然橡胶）复合材料。
2. 2. 模拟方法：
• • 构建SWCNT、聚合物单体→ 无定形晶胞包裹→ 几何优化+动力学弛豫；
• • 预测力学性能（弹性模量、阻尼）、热导率、摩擦学性能。
3. 3. 核心结论：
• • CNT增强可显著提升聚合物刚度、强度、热导率；
• • CNT功能化改善界面结合，提升载荷传递效率；
• • 复合材料阻尼、摩擦性能与CNT-聚合物界面作用强相关。

第5章 陶瓷基纳米复合材料的MD模拟

1. 1. 研究对象：CNT增强SiC陶瓷、Al/Al₂O₃金属陶瓷、BN/CNT同轴复合材料。
2. 2. 模拟方法：
• • 构建陶瓷基+纳米增强相模型，适配陶瓷体系的势函数；
• • 几何优化、动力学模拟，预测力学性能、界面作用。
3. 3. 核心结论：
• • CNT、BNNT可有效提升陶瓷韧性、强度，抑制脆性断裂；
• • 陶瓷-纳米相界面结合决定复合材料整体性能；
• • 同轴纳米管（BN/CNT）展现优异的结构稳定性与力学性能。

第6章 分子动力学脚本编写

1. 1. Materials Studio脚本：建模、模拟、后处理自动化，支持应力应变、热导率、玻璃化转变温度计算脚本。
2. 2. LAMMPS脚本：输入脚本编写（体系定义、势函数、弛豫、拉伸指令），包含空位形成能、纳米线变形、热导率计算脚本。
3. 3. GROMACS脚本：拓扑生成、控温/控压、轨迹处理、性能分析脚本。

第7章 三大软件的跨领域应用

1. 1. Materials Studio：量子力学、介观模拟、材料结晶、高分子设计、催化模拟；
2. 2. LAMMPS：金属塑性变形、纳米摩擦、软物质流变、纳米压痕；
3. 3. GROMACS：生物分子、药物输运、生物芯片、纳米流体、聚合物界面。

全书核心总结

1. 1. 三大软件定位：
• • Materials Studio：可视化、易上手，适合材料建模与快速模拟；
• • LAMMPS：开源、并行强，适合大尺度金属/聚合物体系；
• • GROMACS：高效、界面优，适合生物/聚合物/电解质界面。
2. 2. 模拟核心逻辑：力场选择→模型构建→弛豫优化→动态模拟→性能预测→机理分析。
3. 3. 纳米复合材料强化机理：增强相（CNT、石墨烯、陶瓷颗粒）通过载荷传递、位错阻碍、界面结合提升基体力学/热学性能。