分子动力学模拟器 返回
分子模拟

分子动力学模拟器

使用Lennard-Jones势和Verlet积分计算粒子间的引力与斥力。拖动温度滑块,实时观察气相、液相、固相的转移。

参数

操作

气相
计算结果
运动温度 T*
运动能
势能
总能量
粒子模拟(蓝=低温 → 红=高温)
能量历史
理论·主要公式
$$V(r)=4\varepsilon\left[\left(\frac{\sigma}{r}\right)^{12}-\left(\frac{\sigma}{r}\right)^{6}\right]$$

引力最小点 rmin = 21/6σ ≈ 1.122σ

什么是分子动力学模拟器

🙋
分子动力学模拟可以在计算机上再现教科书中的气体、液体、固体状态变化,是这样吗?
🎓
完全正确!这个模拟器计算了大量粒子(分子)的运动,并以实时方式展示温度或密度变化时会发生什么。例如,尝试将上面的"温度"滑块向下拖动。原本散乱运动的粒子会逐渐聚集并规则排列,即"气体→固体"转移会在眼前发生。
🙋
是吗!但粒子间的相互吸引和排斥力是怎样确定的?
🎓
这正是模拟的关键。粒子间的力用"Lennard-Jones势"公式来确定。参数"ε(epsilon)"表示引力强度,"σ(sigma)"表示粒子大小。增大ε会使粒子更容易粘连,增大σ能模拟更大的粒子。试着改变这些参数,找出液体形成的条件。
🙋
粒子数"粒子数"也可以改,增加粒子数会改变什么?只是计算变得繁重?
🎓
提出好问题!增加粒子数时,少数粒子看不到的"宏观"现象会更容易显现。例如,粒子少时只看起来像一个团块,但当粒子数增加到数百时,降低温度会自然形成整齐的结晶结构。相转移的过程会更清晰。

常见问题

温度滑块调整整个系统的运动能(粒子速度)。在低温下粒子规则排列形成固相(晶体),中温时形成流动的液相,高温时粒子自由飞动形成气相。您可以实时观察这些变化。
当前版本中粒子数是固定的。但通过改变初始配置和温度,可以用相同粒子数再现气体、液体、固体的不同状态。增加粒子数会增加计算负荷,因此出于性能考虑有所限制。
这是表示随距离变化的粒子间力的模型。近距离有强斥力(防碰撞),中距离有弱引力(聚集),远距离近似为零。这个简单公式可以模拟分子间的范德华力。
温度过高时粒子速度很大,相对于计算时间步长,移动量会过大。此时请降低温度滑块,或刷新浏览器返回初始状态。

现实应用

纳米材料设计:使用原子级模拟预测碳纳米管和石墨烯片的机械和热性能。LJ势用于模拟碳原子间的范德华相互作用。

药物设计与蛋白质折叠:通过包含水分子的大规模分子动力学模拟,分析药物分子如何与靶蛋白结合,以及蛋白质如何形成三维结构。LJ势是生物大分子相互作用的基础模型。

润滑剂与界面科学:揭示金属表面与润滑油分子间的摩擦和润滑机制。通过模拟分子如何排列以及如何抵抗剪切力,为开发低摩擦耐久润滑剂做出贡献。

气体储存与分离材料:模拟多孔材料(如MOF)内气体分子(氢气、二氧化碳等)的吸附与扩散行为。通过改变材料结构,可提前评估储气效率。

常见误解和注意点

首先,重要的是不要将"模拟结果视为直接实验数据"。本模拟器只是再现了"Lennard-Jones(LJ)势这个简化模型"的世界。例如,要严格再现水的相转移,需要考虑水分子的极性(正负偏差),需要更复杂的力模型。将LJ模型视为理解稀有气体原子行为的"入门"。

其次,在参数设置中容易陷入的陷阱是设置"极端值"。例如,一下子将温度降至接近0K(绝对零度),粒子几乎不动,模拟看起来"冻结"了。反之,将ε(引力强度)极端增大,粒子会紧密固着不动,反而看不清相转移的"过程"。诀窍是从默认值开始,逐步改变(例如温度每次改变10%,ε每次增加1.2倍),观察其影响。

最后,认为"粒子数越多越好"是误区。确实,数百至数千粒子会明确展现液体或固体的"团块"行为,但计算负荷相应增大。实际应用中,精度与计算资源的平衡很关键。例如,研究纳米粒子凝聚时,详细跟踪数十粒子的行为可能比数千粒子更能抓住本质。用这个工具改变粒子数,掌握捕捉现象本质需要的最少粒子数感觉。

使用指南

  1. 设置粒子数滑块(nMDNum)在2〜1000范围内。氩原子模拟通常从256个开始
  2. 调整温度滑块(tempMDNum)在50〜5000K范围内。观察低温固体状态、中温液体、高温气体相的转移
  3. 设置Lennard-Jones参数ε滑块(epsMDNum),单位为kJ/mol,定义分子间相互作用能
  4. 设置σ滑块(sigmaMDNum)为0.3〜0.5Å单位,确定粒子直径。使用Verlet积分以0.005ps时间步长执行计算

具体计算示例

重现氩(Ar)气液共存状态时,设置粒子数256个、温度94.4K、ε=119.8kJ/mol、σ=3.405Å。初始密度0.8g/cm³,运行10000步模拟,势能收敛至约−600kJ/mol。104K以下呈液相(配位数10〜11),120K以上呈气相(配位数2〜3),可定量确认相图三相点附近的物理化学行为

实际应用注意事项