基于分子动力学模拟探究氢原子在碳纳米管及单壁碳纳米管的储氢行为.docxVIP

基于分子动力学模拟探究氢原子在碳纳米管及单壁碳纳米管的储氢行为

一、引言

二、理论基础与研究方法

2.1分子动力学模拟基本原理

2.1.1牛顿运动方程

分子动力学模拟的核心是基于牛顿运动方程来描述原子的运动。在一个由N个原子组成的体系中，对于第i个原子，其运动遵循牛顿第二定律：

F_{i}=m_{i}a_{i}=m_{i}\frac{d^{2}r_{i}}{dt^{2}}

其中，F_{i}是作用在第i个原子上的力，m_{i}是该原子的质量，a_{i}为加速度，r_{i}表示原子的位置矢量，t代表时间。体系中原子间的相互作用力F_{i}可通过势能函数U(r_{1},r_{2},\cdots,r_{N})对位置矢量求梯度得到：

F_{i}=-\nabla_{i}U(r_{1},r_{2},\cdots,r_{N})

通过对牛顿运动方程进行数值积分，就能得到每个原子在不同时刻的位置和速度，进而模拟出整个体系随时间的演化过程。在实际模拟中，原子间的相互作用是复杂的，势能函数的准确选择至关重要，它直接影响模拟结果的准确性和可靠性。例如，在研究氢原子与碳纳米管的相互作用时，需要选择能够准确描述碳-氢原子间相互作用的势能函数。

2.1.2时间积分算法

在分子动力学模拟中，由于牛顿运动方程是连续的微分方程，需要通过数值方法将其离散化进行求解，时间积分算法就起到了这个作用。Verlet