课程导入：学习目标、课程定位与分子动力学模拟的价值

分子动力学（Molecular Dynamics, MD）模拟是连接微观原子运动与宏观材料性能的关键桥梁。在材料科学中，实验手段往往受限于时间尺度（难以捕捉皮秒级过程）和空间尺度（难以直接观察原子迁移路径），而MD模拟通过数值求解牛顿运动方程，能够以原子级分辨率重现材料在特定条件下的演化过程。本课程面向材料科学与工程、凝聚态物理、计算化学等方向的研究生及科研工程师，旨在帮助学习者建立从"力场选择—模型构建—系综设置—轨迹分析"的完整知识体系，并能够独立开展针对金属、陶瓷、聚合物、纳米材料等体系的模拟研究。

1.1 学习目标：从"会跑代码"到"懂物理、会分析、能设计"

本课程的学习目标分为四个递进层次，每个层次对应不同的能力要求。第一层次是概念理解：掌握分子动力学的基本假设（经典近似、Born-Oppenheimer近似、势函数可迁移性），理解系综（NVE、NVT、NPT、NPH）的物理意义与适用场景。第二层次是操作技能：能够独立完成LAMMPS、GROMACS或OpenMD等主流软件的输入文件编写、力场参数选择、模拟流程搭建，并正确处理单位制转换（如LJ约化单位与真实单位的换算）。第三层次是分析能力：熟练使用VMD、OVITO、MDAnalysis等工具进行轨迹可视化与定量分析，包括计算径向分布函数（RDF）、均方位移（MSD）、速度自相关函数（VACF）、应力张量等。第四层次是研究思维：能够针对具体材料问题（如位错滑移、晶界扩散、聚合物玻璃化转变）设计合理的模拟方案，识别模拟结果的可靠性边界，并将模拟结论与实验现象进行对照。

1.2 课程定位：在计算材料学中的坐标

分子动力学模拟处于计算材料学"多尺度方法"的中间层。在空间尺度上，它介于第一性原理计算（DFT，处理数百个原子）与有限元模拟（FEM，处理宏观构件）之间，典型模拟体系为10³～10⁸个原子，对应尺寸从纳米到亚微米。在时间尺度上，MD可覆盖飞秒至微秒级过程，恰好填补了DFT（皮秒级）与相场/动力学蒙特卡洛（毫秒级以上）之间的空白。这种"承上启下"的定位决定了MD的核心价值：它能够以原子分辨率为宏观本构模型提供参数输入，同时又能验证和补充DFT在有限温度下的预测。

课程内容设计遵循"理论—方法—应用"的螺旋上升结构。前五章夯实物理基础（统计力学要点、系综理论、势函数分类），中间八章聚焦核心技术（积分算法、长程作用力处理、控温控压技术、约束动力学），后十七章面向材料专题（晶体缺陷、非晶结构、界面行为、力学响应、热输运、相变动力学等）。每章均配备实际案例，例如在"位错运动"章节中，将对比分析使用EAM势函数与MEAM势函数模拟FCC铝中刃型位错滑移的差异。

1.3 分子动力学模拟在材料科学中的核心价值

MD模拟在材料科学中的价值可以归纳为三个维度：机制揭示、参数预测、实验指导。在机制揭示方面，MD能够追踪单个原子在纳秒时间尺度内的运动轨迹，直接观察到位错形核、裂纹扩展、晶界迁移等微观过程的原子路径。例如，在纳米压痕模拟中，可以清晰看到压头下方材料发生相变（如硅的体心立方→β锡相变）的原子重排过程，这是实验难以直接观测的。在参数预测方面，MD可以计算材料的弹性常数、热导率、扩散系数、界面能等关键参数，这些参数往往难以通过实验直接测量（如极端条件下的材料性能）。在实验指导方面，MD模拟结果可以帮助实验人员解释异常现象，例如在合金设计中，模拟预测的短程有序结构可以解释实验X射线衍射中出现的额外峰。

以锂离子电池硅负极材料为例：实验发现硅在嵌锂过程中体积膨胀超过300%，导致容量快速衰减。通过MD模拟，研究者发现非晶硅在锂化过程中形成Li-Si合金的原子机制——锂原子优先占据硅网络中的间隙位，随后引发Si-Si键断裂和网络重构。这一模拟结果不仅解释了实验观察到的应力分布特征，还指导了"预锂化"工艺的优化方向。该案例充分体现了MD模拟在"实验现象解释—机理假设验证—工艺参数优化"中的闭环价值。

1.4 整体学习路径与课程结构

本课程共30章，分为四大模块。模块一（第1-5章）为"基础与准备"，涵盖统计力学基础、力场分类与参数化、积分算法原理、周期性边界条件与邻居列表技术。模块二（第6-10章）为"模拟技术进阶"，包括系综控制（Nose-Hoover恒温器、Parrinello-Rahman恒压器）、长程静电力处理（Ewald求和、PPPM方法）、约束动力学（SHAKE算法）、增强采样方法（伞形采样、元动力学）。模块三（第11-18章）为"材料行为分析"，聚焦力学性能（弹性常数、屈服强度、断裂韧性）、热学性能（热导率、比热容）、扩散与输运、相变与成核。模块四（第19-30章）为"专题应用与前沿"，涵盖纳米材料（纳米线、纳米晶）、非晶合金、聚合物与生物材料、界面与摩擦、辐照损伤、电化学界面等方向。

学习路径建议采用"纵向深入+横向拓展"策略。纵向方面，建议学习者从第2章开始，每章完成后完成对应的上机练习，确保每个技术点都能独立操作。横向方面，鼓励学习者选择一种自己感兴趣的材料体系（如高熵合金、钙钛矿太阳能电池材料），在课程进行到第15章后，尝试将所学技术综合应用于该体系的研究。课程最终项目要求提交一份完整的模拟研究报告，包括研究问题定义、模型构建依据、模拟参数选择、结果分析与实验对照。

1.5 常见误区与学习风险提示

初学者在进入MD模拟领域时，常存在以下误区。误区一："力场是通用的"。实际上，每个力场都有其适用边界，例如CHARMM力场针对生物分子优化，EAM势函数适用于金属体系，而Tersoff势函数适用于共价晶体。将金属力场用于聚合物模拟，会得到完全失真的结果。误区二："模拟时间越长越好"。虽然延长模拟时间可以提高统计精度，但受限于计算资源，需要在"模拟时长"与"体系尺寸"之间做权衡。对于扩散系数计算，通常需要模拟时间超过相关弛豫时间两个数量级；而对于玻璃化转变温度测定，则需要足够长的降温速率（通常要求低于1K/ns）。误区三："忽略有限尺寸效应"。当模拟体系尺寸小于相关长度（如位错核心宽度、晶界厚度）时，结果会显著偏离真实材料行为。建议在正式模拟前，进行体系尺寸收敛性测试。

另一个重要风险是对模拟结果的过度解读。MD模拟基于经典力学，无法描述化学键断裂与形成（除非使用反应力场），也无法处理电子极化效应（除非使用极化力场）。因此，当研究涉及电荷转移、化学吸附或光激发过程时，需要谨慎评估MD的适用性。课程将在第5章和第20章分别讨论力场选择策略和MD与DFT的耦合方法，帮助学习者建立"何时该用MD、何时不该用MD"的判断力。

最后，需要强调的是，MD模拟不是"一键出结果"的工具。一个典型的研究级模拟项目，从模型构建到结果分析，通常需要数周时间。其中，模型构建（包括缺陷引入、界面搭建、随机共混）往往占据40%以上的工作量。课程将在第7章专门讲解复杂模型构建技巧，包括使用Atomsk、Moltemplate等工具生成初始构型，以及如何通过"退火-淬火"循环获得合理的非晶结构。掌握这些技能，将显著提升模拟研究的效率与可靠性。