分子动力学方法MOLECULAR DYNAMICS SIMULATION

分子动力学方法MOLECULARDYNAMICSSIMULATION—原理、算法与应用Contents目录分子动力学模拟方法的理论基础、核心算法与前沿应用全景导览。01基本概念与发展历史02理论基础与势能函数03核心算法与系综方法04模拟流程与软件工具05应用领域与前沿进展CHAPTER01基本概念与发展历史从硬球模型到生物大分子：分子动力学的演进之路MolecularDynamics什么是分子动力学模拟分子动力学（MD）模拟是一类基于经典力学原理的计算机模拟技术，通过数值求解牛顿运动方程来追踪体系中每个原子随时间的运动轨迹，从而在微观层面揭示物质的结构、动力学行为与宏观性质之间的内在联系。分子动力学模拟原子运动轨迹可视化01经典力学框架MD模拟基于经典力学框架，将原子视为遵循牛顿第二定律的粒子，通过力场描述原子间相互作用势能02时间步长推进模拟过程在给定初始条件后逐步推进时间步长，记录每个时刻所有原子的坐标与速度，生成完整的运动轨迹03宏观性质推导最终目标是从微观原子运动出发，通过统计力学方法推导出体系的宏观热力学性质与动力学参数HISTORY分子动力学发展简史分子动力学自1950年代末诞生以来，经历了从简单硬球模型到复杂生物大分子体系的跨越式发展。每一次关键突破都与势能函数的改进、计算硬件的升级以及算法创新密切相关，如今已成为物理、化学、生物学和材料科学不可或缺的研究工具。1950s—1970s奠基时期硬球模型开创—1957年Alder与Wainwright首次引入硬球模型进行MD模拟，研究简单液体的相变行为真实势能引入—1964年Rahman使用Lennard-Jones势模拟液态氩，首次将真实势能函数引入MD溶剂化模拟奠基—1974年Rahman与Stillinger完成液态水首次MD模拟，奠定方法论基础1977—2000s生物大分子时期蛋白质模拟先河—1977年McCammon等完成BPTI首次蛋白质MD模拟，开创生物分子模拟经典力场建立—AMBER、CHARMM、GROMOS等力场相继建立，大幅提升模拟精度与可靠性长程静电突破—粒子网格Ewald（PME）方法使长程静电相互作用计算效率获根本性改善2010s—至今现代高性能时期GPU与专用硬件—GPU加速使性能提升数十倍，Anton超算实现毫秒级蛋白质折叠模拟增强采样方法—元动力学、副本交换等方法突破时间尺度瓶颈，拓展可研究过程范围机器学习力场—ANI、DeePMD等将量子力学精度与经典MD效率结合，开启模拟新范式METHODOLOGYCOMPARISONMD与其他计算模拟方法的比较分子动力学在计算模拟方法谱系中占据独特的'介观'位置：相比蒙特卡洛方法它能提供真实的动力学演化信息，相比量子力学方法它能处理更大体系（百万原子级别），相比连续介质方法它保留了原子级分辨率，是连接微观机理与宏观性质的关键桥梁。COMPARISON蒙特卡洛（MC）MC通过随机采样探索构型空间，擅长计算热力学平衡态性质，但无法提供体系随时间演化的真实动力学轨迹MD基于确定性牛顿力学方程积分，能同时获取热力学性质与输运系数、时间相关函数等动力学信息Advantage动力学轨迹COMPARISON量子力学（QM）QM从电子层面求解薛定谔方程，精度高但计算代价为O(N³–N⁴)，通常限于数百原子以内的小体系MD使用经验力场近似原子间作用力，计算代价接近O(N)，可处理数十万至百万原子的大规模体系Advantage百万原子COMPARISON连续介质方法有限元等方法将材料视为连续体，适合宏观工程问题，但无法揭示原子尺度的微观机理MD保留每个原子的自由度，能直接观察位错运动、相变成核、蛋白质折叠等原子级动态过程Advantage原子级分辨率CHAPTER02理论基础与势能函数牛顿运动方程、力场参数化与势能面的数学描述CoreEquations牛顿运动方程与力学基础分子动力学模拟的数学核心是牛顿第二定律：体系势能函数对原子坐标的负梯度给出每个原子所受的力，由此计算加速度并通过数值积分获得原子运动轨迹。这一简洁的物理框架使得从微观相互作用推导宏观性质成为可能。01运动方程对体系中第i个原子，运动方程为mᵢaᵢ=Fᵢ=-∇ᵢV(r₁,r₂,...,rN)，其中V为体系总势能函数F=-∇V02势能函数势能函数V包含键合相互作用（键长、键角、二面角）和非键相互作用（范德华力、静电作用）两大部分键合·非键03数值积分在已知初始位置和速度的条件下，通过数值积分方法（如Verlet算法）逐步更新原子坐标与速度，生成运动轨迹Verlet04宏观桥梁轨迹数据结合统计力学公式可计算温度、压力、扩散系数、热容等宏观可观测量，实现微观到宏观的桥梁T·P·CvMolecularDynamics·ForceField力场的基本概念与组成力场是描述分子体系势能面的参数化数学模型，由键合项和非键项两部分构成，其选择直接决定模拟精度。键合相互作用项01键伸缩势：谐振子模型E=½kb(r−r₀)²描述化学键偏离平衡键长时的能量惩罚02键角弯曲势：谐振模型E=½kθ(θ−θ₀)²描述三个相邻原子键角偏离平衡值的能量变化03二面角扭转势：余弦级数E=ΣVn[1+cos(nφ−γ)]描述绕化学键旋转时的能量周期性变化非键相互作用项01范德华作用：Lennard-Jones12-6势E=4ε[(σ/r)¹²−(σ/r)⁶]描述原子间短程排斥与长程色散吸引包含排斥项（r⁻¹²，近距离主导）和吸引项（r⁻⁶，色散作用），参数ε为势阱深度，σ为原子直径02静电作用：库仑定律E=q₁q₂/(4πε₀r)描述带电或极性原子之间的远程静电相互作用随距离衰减较慢（r⁻¹），长程效应显著，通常采用Ewald求和或粒子网格方法处理周期性边界条件MOLECULARDYNAMICS·FORCEFIELDS主流分子力场体系概览主流分子力场经过数十年的发展已形成多个成熟体系，各自针对特定类型的分子体系进行了系统性的参数化与验证。选择力场时需考虑研究体系的化学类型、关注的物理性质以及力场在该领域的验证充分度，不恰当的力场选择是导致模拟结果不可靠的最常见原因之一。主流分子力场对比力场名称开发团队主要适用体系核心特点AMBERKollman(UCSF)蛋白质、核酸生物大分子模拟标杆，参数体系完善CHARMMKarplus(Harvard)蛋白质、脂质、核酸覆盖体系最广，支持QM/MM混合模拟GROMOSvanGunsteren(ETH)蛋白质、溶液体系联合原子模型，计算效率高OPLS-AAJorgensen(Yale)有机分子、液态体系液态热力学性质重现精度突出GAFF/GAFF2Kollman/Wang通用有机小分子AMBER兼容，可自动化参数生成GLYCAMWoods(UGA)糖类、糖蛋白专为碳水化合物设计，糖苷键参数精确六大主流力场各有侧重，AMBER和CHARMM在生物大分子领域应用最广，OPLS-AA在液态性质上表现突出POTENTIALFUNCTIONLennard-Jones势与范德华相互作用Lennard-Jones12-6势函数是描述范德华相互作用最广泛使用的经验模型，通过r⁻¹²排斥项和r⁻⁶吸引项的组合精确捕捉原子间短程排斥与长程色散吸引的物理本质。参数σ和ε通过混合规则（如Lorentz-Berthelot）可推广到异种原子对的相互作用计算。Lennard-Jones12-6势能函数曲线示意01势能公式：E(r)=4ε[(σ/r)¹²−(σ/r)⁶]，其中σ为零势能距离、ε为势阱深度，分别控制排斥范围和吸引强度02物理本质：r⁻¹²项近似描述电子云重叠时的泡利排斥效应，r⁻⁶项精确描述瞬时偶极-诱导偶极的伦敦色散吸引力03混合规则：Lorentz-Berthelot规则σᵢⱼ=(σᵢ+σⱼ)/2，εᵢⱼ=√(εᵢεⱼ)，将同种原子参数推广到异种原子对04截断处理：采用截断半径（通常1.0–1.4nm）限制计算范围，配合切换函数或长程色散修正减少截断误差LONG-RANGEELECTROSTATICS静电相互作用与长程处理方法库仑静电作用随1/r缓慢衰减，属于典型的长程相互作用，简单截断会引入严重的能量不连续和人为伪影。Ewald求和及PME方法通过实空间与倒空间分解配合FFT实现O(NlogN)高效精确计算，已成为现代MD模拟的标准做法。库仑势的长程性E=q₁q₂/(4πε₀r)随距离1/r缓慢衰减，在周期性边界条件下直接求和条件收敛，简单截断会产生严重伪影。1/rEwald求和分解将静电势分为实空间（短程、快速收敛）和倒空间（长程、傅里叶级数求和）两部分，确保总能量精确收敛。实空间+倒空间PME方法粒子网格Ewald将倒空间电荷映射到网格，利用FFT加速，复杂度从O(N²)降至O(NlogN)。O(NlogN)简化替代方案反应场方法和截断加移位方法是PME的简化替代，计算更快但精度和物理一致性不如PME严格。精度权衡Chapter03核心算法与系综方法数值积分、统计系综与温度/压力控制算法NumericalIntegrationVerlet系列数值积分算法Verlet系列算法是分子动力学模拟中最核心的数值积分方法，通过对牛顿运动方程进行泰勒展开近似实现位置与速度的时间步进更新。VelocityVerlet因其同时给出位置和速度、良好的能量守恒性和辛结构保持特性，已成为现代MD程序中最广泛采用的标准积分器。基本Verlet算法利用泰勒展开：r(t+Δt)=2r(t)−r(t−Δt)+a(t)Δt²，仅需存储当前位置和前一步位置，内存占用少优点为时间可逆且辛结构保持，长时间模拟能量漂移极小；缺点是不直接给出速度，需额外计算辛结构VelocityVerlet算法分两步更新：先用位置公式更新位置，再用v(t+Δt)=v(t)+½[a(t)+a(t+Δt)]Δt更新速度同时提供精确的位置和速度信息，便于计算动能和温度，是GROMACS、LAMMPS等主流软件的首选积分器首选积分器Leap-frog蛙跳算法速度在半点t+½Δt处更新，位置在整点更新：r(t+Δt)=r(t)+v(t+½Δt)Δt，两者交替计算计算量略低于VelocityVerlet，但速度不与位置同步，需额外插值获取同时刻速度半点步进MolecularDynamics周期性边界条件与最小镜像约定周期性边界条件（PBC）通过将模拟盒子在三维空间无限复制来消除表面效应，使有限的N粒子体系能近似模拟宏观块体材料的行为。配合最小镜像约定，每对粒子的相互作用仅取最近镜像距离计算，既保证了物理合理性又控制了计算复杂度，是现代MD模拟的基础设施级技术。01PBC将中心模拟盒子在三个维度无限平移复制，粒子离开一个面时从对面重新进入，保持粒子数N恒定并消除自由表面N=const02最小镜像约定：任意两粒子间的距离取其在所有镜像中最近的那个，确保每对粒子只计算一次相互作用，截断半径需≤盒长的一半rc≤L/203PBC隐含体系具有周期性结构的假设，对均匀液体和晶体是良好近似，但对界面、缺陷等非均匀体系可能引入人为周期性伪影ARTIFACTRISK04模拟盒子形状可选立方体、长方体、菱形十二面体或截断八面体，后两者在球形溶质体系中能减少约30%的溶剂分子数≈30%↓周期性边界条件中模拟盒子与三维镜像复制示意图OPTIMIZATIONSTRATEGIES邻居列表与计算优化策略Verlet与Cell列表将力计算复杂度从O(N²)降至O(N)，使百万原子级长时间模拟成为可能01Verlet邻居列表为每个原子构建邻居清单（截断半径rc+缓冲层rskin），缓冲层允许列表在10–20步内复用而无需重建02Cell空间列表将模拟盒子划分为边长≥rc的小格子，仅需搜索所在格子及26个相邻格子，复杂度降为O(N)03双重优化架构现代MD程序（GROMACS、LAMMPS）将Cell列表用于快速构建Verlet列表，再用Verlet列表加速每步力计算04GPU与并行化GPU加速与域分解并行化进一步优化大规模体系计算效率，使百万原子体系纳秒级日产出成为常规条件STATISTICALENSEMBLE统计系综的概念与分类统计系综定义了分子动力学模拟中体系的热力学约束条件，直接决定了模拟所采样的概率分布和可获得的热力学量。NVE、NVT、NPT三种系综覆盖了绝大多数模拟需求，正确选择系综是确保模拟结果具有明确物理意义的前提条件。MD模拟常用统计系综对比系综恒定物理量物理对应实现方法NVE（微正则）粒子数N、体积V、能量E完全孤立体系，能量守恒基本Verlet积分，无额外控制NVT（正则）粒子数N、体积V、温度T与恒温热浴平衡的体系thermostat（热浴耦合）NPT（等温等压）粒子数N、压力P、温度T最接近常规实验条件thermostat+barostat双重耦合μVT（巨正则）化学势μ、体积V、温度T与粒子库交换粒子的开放体系蒙特卡洛粒子插入/删除NVE适合能量守恒验证，NPT最接近实验条件，NVT适用于固定密度研究，μVT用于相平衡与吸附研究ThermostatMethods温度控制算法（Thermostat）温度控制算法通过调节原子速度使体系温度维持在目标值附近，不同thermostat在收敛速度、系综采样的严格性和对动力学信息的干扰程度上各有优劣。Berendsen热浴通过指数弛豫机制按比例缩放速度：vnew=λ·vold，使温度以特征时间常数τT向目标温度T₀指数趋近收敛速度快且稳定，适合初始平衡阶段；不能严格产生正则系综概率分布，可能导致"飞冰棒"效应λ·vNosé-Hoover热浴引入虚拟热浴变量s作为额外自由度，通过扩展拉格朗日量严格实现正则系综采样，理论基础最为扎实小体系或刚性体系可能出现非遍历性问题，可采用Nosé-Hoover链（NHC）引入多个串联热浴变量解决NHC速度重采样（V-rescale）Bussi等人提出的改良方法，在Berendsen缩放基础上加入随机项，既保持快速收敛又能产生严格正则系综分布实现简单、数值稳定、采样正确，已成为GROMACS默认thermostat，适用于绝大多数模拟场景GROMACSBarostat·NPTEnsemble压力控制算法（Barostat）压力控制算法通过动态调节模拟盒子的体积（和形状）使体系压力维持在目标值，是NPT系综模拟的核心组件。不同barostat在各向同性/各向异性压力耦合、系综严格性和数值稳定性方面存在差异，正确的使用策略是先用简单方法预平衡再切换到严格方法正式采样。Berendsen压控按比例缩放盒子尺寸使压力以特征时间τ_P向目标值指数趋近，收敛快速但不产生严格NPT系综分布。τ_P指数趋近Parrinello-Rahman方法将盒子向量作为动力学变量引入扩展拉格朗日量，允许盒子形状变化，适合各向异性体系和相变研究。各向异性耦合Martyna-Tobias-Klein（MTK）在Nosé-Hoover框架下实现严格等温等压采样，理论基础最完善，是NAMD等软件的首选barostat。严格NPT系综推荐工作流程平衡阶段使用Berendsen快速收敛体积，采样阶段切换为Parrinello-Rahman或MTK以保证系综采样的正确性。预平衡→正式采样Chapter04模拟流程与软件工具从体系构建到轨迹分析的完整工作流程WORKFLOWMD模拟完整工作流程分子动力学模拟是一个多步骤、有严格先后依赖关系的系统化流程：从初始结构获取到力场参数化，从溶剂化到能量最小化，从NVT/NPT平衡到正式采样，每一步都必须达到特定的收敛判据后才能进入下一步。规范的操作流程是获得可靠模拟结果的基本保障。01体系准备从PDB数据库获取实验结构或通过同源建模、从头预测构建初始原子坐标，检查缺失残基和原子，确保结构完整性02力场分配为体系中每种分子选择力场并生成拓扑文件，包含原子类型、电荷、键合参数和非键参数等完整信息03溶剂化与中和将溶质置于模拟盒子中，填充水分子（如TIP3P、SPC/E模型），加入Na⁺/Cl⁻离子使体系电中性并设定离子浓度04能量最小化使用最速下降法或共轭梯度法消除初始结构中的原子重叠和高能构型，使体系势能降至合理范围05NVT平衡在位置限制下加热体系至目标温度（如300K），验证温度稳定波动后解除限制进入下一步06NPT平衡与生产模拟在恒温恒压下平衡体系密度，确认能量、温度、压力均稳定后开始正式轨迹采集OPTIMIZATIONMETHODS能量最小化方法能量最小化通过优化算法将初始结构调整到最近的势能面局部极小点，消除原子重叠、不合理键长等高能冲突，为后续动力学模拟提供物理上合理的起始构型。最速下降法沿势能面负梯度方向移动原子，步长由线搜索确定，每步独立计算方向，无需存储历史信息，内存占用极低鲁棒性强、能处理严重原子重叠；接近极小点时收敛极慢（锯齿形振荡），适合前期粗优化负梯度搜索共轭梯度法利用当前与历史梯度构造共轭搜索方向，避免最速下降法的锯齿问题，收敛速度显著更快需存储上一步梯度信息，适合初步消除冲突后的精细化优化，对二次势能面可达超线性收敛共轭方向L-BFGS拟牛顿法用有限步历史梯度近似构造Hessian矩阵的逆，无需显式存储完整Hessian即获得接近二阶方法的收敛速度对大分子体系效率优异，是GROMACS等现代MD软件中大体系能量最小化的推荐方法GROMACS推荐SolventModels&Solvation溶剂模型与溶剂化策略溶剂环境的处理方式对MD模拟的精度和效率有决定性影响，选择时需权衡溶剂细节依赖与计算资源。显式水模型01TIP3P：三站点刚性模型，计算效率高，AMBER和CHARMM力场默认搭配，水扩散系数偏高02SPC/E：改进偶极矩和液态密度，适用于GROMOS力场，自洽极化修正提升热力学精度03TIP4P/2005：四站点模型，虚拟电荷精确重现水的密度异常、相图和热容等物理性质隐式溶剂模型01广义Born（GB）模型：解析公式近似静电溶剂化能，速度快10–100倍，适合大规模构象采样和自由能估算。通过有效玻恩半径简化溶剂-溶质边界处理，在蛋白质折叠和分子对接中广泛应用02Poisson-Boltzmann（PB）方程：数值求解连续介质静电学方程，精度高于GB，常用于结合自由能后处理。通过有限差分或有限元方法求解非线性方程，可精确描述复杂分子表面的静电势分布精度-效率权衡：显式模型适合研究溶剂化层结构和氢键网络，隐式模型适合探索长时尺度动力学和构象空间TOOLS&PLATFORMS主流MD模拟软件工具分子动力学模拟软件经过数十年发展已形成多个成熟的开源和商业平台，各自在计算效率、力场支持、并行扩展性和易用性等方面具有不同优势。选择合适的软件需综合考虑研究体系类型、所需功能、可用计算资源和团队的熟悉程度，实际研究中多种软件配合使用也十分常见。主流MD软件对比一览6款平台·4维度软件开发机构核心优势许可GROMACSKTH/Uppsala等计算效率极高，GPU优化出色，社区活跃开源免费AMBERStonyBrook/UCSF等生物分子力场标杆，分析工具丰富商业授权NAMDUIUC(Schulten)并行扩展性优异，与VMD深度集成学术免费LAMMPSSandia国家实验室材料科学全能，支持多种势函数开源免费OpenMMStanford(Pande)Python接口灵活，适合方法开发开源免费CHARMMHarvard(Karplus)QM/MM混合模拟支持完善商业授权GROMACS和LAMMPS是开源领域的主力，AMBER和CHARMM在生物分子领域深耕，NAMD和OpenMM各有特色AlgorithmOptimization键长约束算法与时间步长优化键长约束算法通过数学方法将含氢键等高频振动自由度固定在其平衡值，从而消除限制时间步长的最快运动模式。这一技术可将时间步长从0.5-1fs扩大到2fs甚至更大，使模拟效率提升2-4倍，是现代MD模拟中不可或缺的效率优化手段。SHAKE与RATTLE算法SHAKE在Verlet积分的位置更新步骤后，通过迭代拉格朗日乘子法将被约束键长调整回目标值，精度由收敛阈值控制RATTLE是SHAKE在VelocityVerlet框架下的对应版本，同时约束位置和速度，保证被约束键上无速度分量迭代求解LINCS算法线性约束求解器通过矩阵展开近似直接求解约束方程，避免了SHAKE的迭代过程，计算效率更高特别适合大规模并行计算和大体系模拟，是GROMACS的默认约束方法，通常允许2fs时间步长且能量守恒良好GROMACSSETTLE算法专为水分子设计的解析约束方法，一步直接求解三个约束方程（两个O-H键和H-H距离），无需迭代由于水分子数量在体系中占比极大，SETTLE的解析求解显著降低了溶剂约束的计算开销解析求解CHAPTER05应用领域与前沿进展从药物设计到材料科学的多维应用场景MOLECULARDYNAMICS·CHAPTER27蛋白质折叠与构象动力学分子动力学模拟为蛋白质折叠研究提供了实验手段难以获得的原子级时间分辨信息。D.E.Shaw团队使用Anton超级计算机完成的毫秒级全原子模拟首次直接观察了蛋白质自发折叠的完整动力学过程，揭示了折叠漏斗中多路径并行的复杂图景，深化了对蛋白质能量景观理论的理解。01首次观测：Anton超级计算机完成了多个小型蛋白质（如FiP35WW结构域）的毫秒级全原子模拟，首次直接观测到自发折叠事件02多路径假说：MD模拟揭示了折叠路径的多样性，同一蛋白质可通过多条不同路径到达天然态，支持能量漏斗理论03增强采样：副本交换MD、元动力学等方法通过克服自由能壁垒加速构象空间探索，弥补常规MD时间尺度的不足04AI互补：AlphaFold2给出静态终点结构，MD揭示到达终点的动态路径和中间态信息，两者形成研究互补蛋白质三维分子结构的科学可视化渲染MolecularDynamics计算机辅助药物设计中的应用分子动力学模拟在药物设计中弥补了分子对接方法忽略蛋白质柔性和溶剂效应的不足，通过模拟药物-靶蛋白复合物在真实溶液环境中的动态行为，提供更可靠的结合模式评估和更精确的结合自由能预测，已成为现代药物研发流程中从先导化合物优化到候选药物评价的关键计算工具。结合模式评估与优化对接结果经MD模拟松弛后，蛋白质侧链和骨架的柔性调整使结合模式更接近真实构象，减少假阳性预测。水分子介导的氢键网络和疏水口袋的动态变化可通过MD直接观察，指导药物化学家设计更好的取代基团。构象优化结合自由能计算MM-PBSA/GBSA方法从MD轨迹中提取能量项进行结合自由能估算，计算成本低，适合大批量候选物的初步筛选排序。自由能微扰（FEP）和热力学积分（TI）方法提供更高精度的相对结合自由能，能预测0.5kcal/mol以内的活性变化。FEP/TI变构效应与选择性分析MD模拟可揭示药物结合引起的远端变构效应，帮助理解变构调控机制并设计变构药物。通过比较药物与同家族不同亚型蛋白的MD模拟，分析选择性来源，指导高选择性药物设计。变构调控MolecularDynamics膜蛋白与脂质双分子层模拟膜蛋白模拟是MD领域最具技术挑战性的方向之一，需要在磷脂双分子层、水溶液和离子的复杂多相环境中模拟蛋白质行为。全原子模拟能揭示离子通道的选择性机制和受体激活的构象变化，粗粒化模型则将可模拟的时间和空间尺度拓展数倍，两者结合构成了研究膜蛋白功能的完整计算框架。膜蛋白嵌入磷脂双分子层的分子动力学模拟可视化01膜蛋白模拟需构建磷脂双分子层环境（如POPC、DPPC），体系通常包含10万-100万原子，对计算资源要求远高于溶液蛋白模拟02离子通道选择性机制：MD模拟揭示了钾离子通道KcsA中羰基氧排列形成的选择性滤器如何精确区分K⁺和Na⁺03GPCR激活机制：长时间MD模拟观察到受体从非活性态到活性态的构象转变路径，揭示了胞内信号转导的分子开关04MARTINI粗粒化力场将约4个重原子映射为1个珠子，可模拟微秒到毫秒尺度的膜重组、蛋白质聚集等大尺度膜过程APPLICATIONS材料科学中的MD模拟应用分子动力学在材料科学中的应用涵盖了从金属合金到聚合物、从纳米材料到陶瓷的广泛领域。通过选择适当的势函数，MD能揭示材料的微观变形机理与相变动力学。金属与合金EAM嵌入原子势精确描述金属键合特性，可模拟位错成核与运动、晶界迁移、裂纹扩展等塑性变形微观机制高温蠕变和辐照损伤模拟揭示了缺陷的产生、演化和湮灭动力学，为核材料设计提供关键数据EAM·位错·辐照聚合物与软物质全原子和粗粒化MD模拟可研究聚合物链段运动、玻璃化转变温度预测和力学性能评估高分子共混物的相分离动力学、嵌段共聚物的自组装过程可通过大规模MD模拟直接观察链段·相分离·自组装纳米材料与界面碳纳米管和石墨烯的拉伸力学、热导率等性质可通过MD精确预测，与实验结果高度吻合纳米颗粒-基体界面、液-固界面的分子级结构和输运性质可通过MD揭示，指导纳米复合材料设计CNT·石墨烯·界面FREEENERGYMETHODS自由能计算方法自由能计算将分子动力学模拟从定性的结构动态描述提升为定量的热力学预测工具。PATH-BASED热力学积分（TI）定义耦合参数λ连接初末态，计算每个λ窗口的⟨∂H/∂λ⟩并数值积分得到ΔG，精度高但需要多个窗口。⟨∂H/∂λ⟩PATH-BASED自由能微扰（FEP）基于Zwanzig公式，适合相邻状态间的小扰动计算，是药物优化的核心工具。ΔG=−kT·ln⟨e−ΔU/kT⟩ENHANCEDSAMPLING伞形采样沿反应坐标分窗口施加偏置势增强采样，用WHAM或MBAR重构无偏自由能面（PMF），适合一维反应坐标。PMFENHANCEDSAMPLING元动力学在CV空间逐步沉积高斯势填平自由能低谷，迫使体系逃离局部极小并探索新区域，可在线估计自由能面。CVComputationalMethodQM/MM混合量子-经典方法QM/MM方法将体系划分为量子力学区域和分子力学区域，使含数万原子的复杂体系反应研究成为可能，Warshel与L