  在分子动力学中，分子被描述为一系列带电的点（原子），通过化学键相互连接。 为了描述键长、键角、扭转角的时间演化，以及原子之间的非键范德华力和静电相互作用, 人们使用力场来描述这些参数，力场通过参数化这些能量项的公式和常数，将原子之间的相互作用转化为力和能量的数学表示。这些参数可以通过实验数据、量子化学计算或经验拟合等方法获得。 通过力场，可以计算分子系统中原子之间的相互作用力，从而模拟分子的结构、动力学和热力学行为。力场在分子模拟和计算化学研究中起着重要的作用，能够预测分子的构象、热力学性质、反应机理等重要信息。不同的力场可能适用于不同类型的分子和研究问题，选择合适的力场对于获得准确和可靠的模拟结果至关重要。 分子动力学力场通常由以下文件组成： 1. 参数文件（Parameter files）：参数文件包含了力场所需的参数信息，如原子质量、电荷、键长、键角、二面角等。这些参数是通过实验数据、量子力学计算或经验拟合获得的。 2. 拓扑文件（Topology files）：拓扑文件定义了分子系统的拓扑结构，包括原子之间的连接方式、键类型、键长、键角等信息。它通常与参数文件一起使用。 3. 势能函数文件（Potential function files）：势能函数文件包含了力场所使用的势能函数的具体形式和参数。它定义了键能、角能、二面角能和非键相互作用能的计算方式。 4. 输入文件（Input files）：输入文件是包含了系统结构、参数文件和模拟参数的文件。它指定了模拟的起始结构、模拟步长、温度等参数，并指定使用的力场和相应的文件。 势能函数在分子动力学力场中由几个主要的组成部分构成： 1. 键能（Bond Energy）：描述共价键的能量。它基于键的长度和键的弹性常数来计算，用于模拟原子之间的共价键振动。 2. 角能（Angle Energy）：描述角度的能量。它基于角度的大小和角度的弹性常数来计算，用于模拟原子之间的键角变化。 3. 二面角能（Dihedral Energy）：描述二面角（扭转角）的能量。它基于二面角的大小和周期性函数来计算，用于模拟原子之间的扭转和旋转。 4. 非键相互作用能（Nonbonded Interaction Energy）：描述非共价相互作用的能量，包括静电相互作用和范德华力。它通常由库仑势和范德华势来表示，用于模拟原子之间的静电吸引和排斥以及范德华力的作用。 - 库仑势能（Coulombic Potential）：用于描述带电粒子之间的静电相互作用。它与电荷之间的距离成反比，通过库仑常数和电荷值来计算。 - 范德华势能（Van der Waals Potential）：用于描述非极性相互作用，如范德华力。它基于原子之间的距离和范德华力常数来计算，通常采用Lennard-Jones势函数或其修正形式。  Ubond=共价震动能量 Uangle=三个原子组成的键角震动能量 Udihedrals=四个原子组成的二面角旋转势能 Unonbond=非键相互作用能量，包括静电相互作用和范德华势能； 原子之间的键结相互作用：  键长与键能之间的相互关系如下图所示:
键角和不当角（improper terms）具有类似的二次形式，但弹性系数较小。这些力常数可以通过分子的振动分析获得。 二面角势能与二面角大小的相互关系如下图所示：  其中n表示是周期性项的整数值，表示旋转的周期性； 静电势能和原子之间距离的关系下图所示：
范德华势能和距离的相互关系如下图所示：  下面是一些常见的分子动力学力场文件： 1. AMBER力场：AMBER（Assisted Model Building with Energy Refinement）力场是广泛应用的力场之一。它提供了多个版本，如AMBER94、AMBER99和AMBER ff03等。AMBER力场主要用于蛋白质和核酸的模拟研究。 2. CHARMM力场：CHARMM（Chemistry at Harvard Molecular Mechanics）力场是另一个常用的力场。它包括多个版本，如CHARMM22和CHARMM36等。CHARMM力场广泛应用于蛋白质、核酸、脂质和碳水化合物等生物分子的模拟。 3. GROMOS力场：GROMOS（Groningen Molecular Simulation）力场是由Groningen大学开发的力场。它适用于溶液中的生物分子模拟，并提供了多个版本，如GROMOS53A6和GROMOS96等。 4. OPLS力场：OPLS（Optimized Potentials for Liquid Simulations）力场是用于液体模拟的力场。它提供了多个版本，包括OPLS-AA（全原子力场）和OPLS-UA（United Atom力场）。OPLS力场广泛应用于脂质、多肽和蛋白质等系统的模拟。 5. AMBER-ff系列力场：AMBER-ff（AMBER Force Field）系列力场是AMBER力场的衍生版本，针对特定类别的分子进行了参数化。例如，AMBER-ff99SB力场专门用于蛋白质模拟，而AMBER-ff99SB-ILDN力场用于蛋白质折叠和去折叠的模拟。 这些力场文件通常以特定的格式（如AMBER格式、CHARMM格式或GROMOS格式）存储，并包含描述分子的拓扑信息和相互作用参数。此外，还有其他一些力场，如MMFF（Merck Molecular Force Field）、UFF（Universal Force Field）等，它们适用于小分子和药物化学等领域的模拟研究。 参考资料： https://www.ks./Training/Workshop/Boston/emad_lecture 关注我获得 |
|
|
