分子动力学简介
分子动力学（Molecular Dynamics, MD）是一种数值模拟方法，用于研究分子系统如何随时间演化。这些模拟为我们提供了关于原子和分子如何互相作用、移动和分布的深入洞察。
基础知识
1. 原理 分子动力学是基于经典物理学，尤其是牛顿的第二定律来模拟原子和分子的动态行为的方法。在这种模拟中，原子或分子被视为质点，仅考虑其质量和位置，忽略其内部结构。它们的运动由键合和非键合相互作用的力决定，这些力来自于它们与其他原子或分子的相互作用，从而通过这种方法可以预测和理解复杂系统中的微观行为。
2. 力场 在分子动力学模拟中，原子和分子之间的相互作用决定了其运动。这些相互作用由“力场”来描述。 力场分为 键合 和 非键合相互作用 。键合相互作用包括 键长 （两键合原子间的理想距离）和 键角 （三键合原子间的夹角）。 非键合相互作用涉及 范德华力 （原子间的吸引或排斥）和 静电相互作用 （带电的原子或分子之间的吸引或排斥）。常见的力场如 AMBER 主要用于蛋白质和核酸， CHARMM 适用于各种生物分子， OPLS 主要针对有机分子，而 GROMOS 常用于蛋白质模拟。 C-H-O体系，在有机化学、生物化学和材料科学中均有广泛应用，因此有多个力场适用于其模拟。常用的有AMBER的GAFF版本、 CHARMM、OPLS-AA、GROMOS、MM3/MM4、UFF、DREIDING和COMPASS。其中，AMBER和CHARMM主要用于生物大分子，而OPLS-AA、MM系列和DREIDING则广泛应用于有机分子。 UFF是一个全元素通用力场，而COMPASS主要针对材料科学。 选择合适的力场需考虑研究问题和系统特性，确保模拟的准确性。
3. 时间步 分子动力学模拟是一种计算方法，旨在模拟和预测分子系统随时间的演化。为了达到这个目的，我们必须对系统进行“离散化”，也就是将连续的时间分解成一系列小的时间间隔或“时间步”。每个时间步中，我们都会根据物理定律计算原子和分子的位置和速度的变化。 选择合适的时间步非常关键。由于原子和分子在高频率下发生快速振动，特别是在高温和高能状态下，我们需要足够小的时间步来确保捕捉到这些运动。通常，这些时间步的大小是飞秒级别的，即10^-15秒。例如，使用1飞秒或10飞秒作为时间步可以确保在模拟中捕获到化学键的振动，如碳氢键。 为什么要这么小的时间步呢？考虑一个简单的弹簧振子模型，如果我们采用的时间步太大，可能会错过振子的某些振动周期，导致能量不守恒，进而影响模拟的准确性。这种现象在真实的分子系统中也同样存在。原子的快速振动，特别是在固体和液体中，需要足够细致的时间分辨率来观察。 此外，太大的时间步可能会导致数值不稳定，使模拟系统在计算过程中“爆炸”，这意味着原子可能会被突然推到非常远的地方，导致模拟失效。
4. 初始条件 在分子动力学模拟启动之初，为系统内的每个原子分配初始位置和速度是至关重要的步骤。原子的初始位置通常是基于实验结构数据，如X射线晶体学或核磁共振结构。如果这些数据不可用，可能还会参考先前模拟的结构或通过量子力学计算得出的预测结构。而关于速度，它们通常基于所考虑温度的玻尔兹曼分布随机初始化。确保原子有适当的初始速度非常重要，因为它决定了模拟的动态行为和系统的初始热化。如果速度选择不恰当，系统可能不会正确地模拟，或需要额外的时间来达到平衡态。 因此，从合适的来源为原子选择初步位置和速度是确保模拟质量的关键。在为分子动力学模拟设定原子的初始位置和速度时，通常可以采用从实验如X射线、NMR获取的结构，或基于已知结构的模型构建，或使用计算方法进行预测。速度则可从玻尔兹曼分布中根据指定温度随机选取，或继承自先前的模拟。此外，分子动力学软件包如GROMACS、AMBER和NAMD提供了设定这些初始条件的工具和程序。选择具体方法取决于所研究的系统和可用的数据。
5. 边界条件 在分子动力学模拟中，为了避免由于有限系统大小而产生的边缘效应，并且更好地模拟真实的宏观条件，常用的边界条件是周期性边界条件(PBC)。PBC的核心思想是将一个小的、有限的模拟盒子在各个方向上无限次复制，从而创造了一个似乎是无限扩展的宏观系统。当模拟中的一个原子移动并穿过模拟盒子的一个面时，另一个与其相同的原子会从相对的另一面进入，确保系统内的原子数保持不变。这种方法有效地模拟了一个连续、无限的材料，同时在计算上仍然是可行的，因为仅计算原始模拟盒中的相互作用。
6. 热力学系综 分子动力学模拟通常在特定的热力学条件下进行，例如恒定的温度和压力（NPT系综）或恒定的粒子数、体积和温度（NVT系综）。此外，还有NVE系综，也称为微正则系综。在NVE系综中，系统的粒子数（N）、体积（V）和能量（E）都是恒定的。这意味着在整个模拟过程中，粒子的总数是不变的，没有粒子进入或离开系统。体积也是固定的，而总能量（包括势能和动能）是恒定的，这意味着没有能量从系统中流入或流出。NVE系综是最简单的模拟系综之一，因为它不需要任何形式的热浴或压力耦合。但是，由于系统能量是恒定的，因此系统的温度可能会因为动能的变化而波动。在现实的实验室条件下，恒定能量的情况是不太常见的，但NVE系综在研究一些特定问题时是非常有用的，例如研究系统的基本动力学或在没有外部干扰的情况下观察系统的自然演化。
7. 应用 分子动力学广泛应用于化学、物理、生物和材料科学中，用于研究液体的性质、生物大分子的构象变化、材料的机械性质等。
8. 局限性 由于计算能力的限制，模拟的时间尺度通常限制在纳秒到微秒的范围内。而且，模拟的准确性高度依赖于所选的力场。
9. 工具和软件 有许多软件可用于进行分子动力学模拟，包括GROMACS、LAMMPS、NAMD和AMBER。
10. 分析和可视化 模拟结束后，可以对轨迹数据进行各种分析，例如原子的平均位置、速度分布、能量分布等。常见的可视化工具包括VMD、PyMOL和Chimera。
进阶主题
1. 增强采样技术 增强采样技术致力于更高效地探索系统的配置空间，特别是那些难以直接访问的高能区域。例如，Metadynamics方法通过在势能表面上不断添加Gaussian势垒，使得系统能够越过能垒并探索新的相空间区域。Umbrella Sampling则是在反应坐标上分段采样，并使用偏置势来确保均匀地探索整个反应路径。最后，加权直方图分析方法 (WHAM) 是一种重新加权多个模拟中的数据以获得统计上准确的自由能曲线的方法。
2. 非等温-非等压系综 为了模拟在特定温度和压力下的系统，采用了不同的恒温和恒压算法。Berendsen热统和压力统为系统提供了快速的弛豫，是为了快速地达到所需的温度或压力；Nose-Hoover热统是基于Hamiltonian的方法，为系统提供一个更严格的温度控制；Parrinello-Rahman方法则调整模拟盒的形状和体积以维持所需的压力。
3. 多尺度模拟 多尺度模拟结合了不同描述水平的方法，从量子到经典，从原子到介观。例如，某些关键的电子结构变化可以通过量子化学方法来描述，而模拟的其余部分则使用原子级的经典描述，同时还可以利用介观级的描述来模拟更大的系统尺度，如微米尺度的流体流动。
4. 粗粒化模型 粗粒化方法是一种降维策略，其中多个原子或分子被表示为一个效果粒子或块。这种简化使得模拟可以扩展到更大的时间和空间尺度，但牺牲了一些详细信息。这对于研究如膜、聚合物和生物大分子的组装和动态尤为有用。
5. 反应动力学 反应动力学模拟主要关注原子尺度上化学反应的速率和机制。转换态理论基于能量表面上的鞍点来估计反应速率，而动态隧道效应考虑了量子效应，尤其是在低温下的反应速率增强。
6. 非平衡分子动力学 (NEMD) NEMD是一种模拟非热平衡状态下的系统的方法。通过施加温度梯度、剪切率或其他非平衡驱动力，可以研究如热传导、粘度和其他输运性质的系统响应。
7. 量子分子动力学 量子分子动力学结合了量子力学描述电子的波函数和经典描述原子核的运动方程。这种方法特别适用于电子结构和原子运动之间有强烈耦合的系统，例如在化学反应中的键断裂和形成。
8. 自由能计算 自由能与系统在特定条件下的稳定性直接相关。自由能摄动 (FEP) 和热力学积分 (TI) 是常用的计算相对自由能变化的方法，例如在分子改造或配体结合中的自由能差异。
9. 拉伸模拟 在拉伸模拟中，通过施加外部力来逐渐改变分子的形态，这有助于了解分子在外部应力下的稳定性、机械性质和可能的变形路径。
10. 蛋白质-配体结合模拟 这种模拟专门研究药物小分子与生物大分子（如蛋白质或核酸）之间的相互作用。通过观察结合过程和计算结合自由能，为药物设计和优化提供了有价值的洞察。