共轭梯度算法计算团簇的最稳定结构.docVIP

共轭梯度算法计算团簇的最稳定构型 摘要：采用共轭梯度算法（conjugate gradient method），从初始给定的位置出发，计算了具有Lennard-Jones势的团簇（一般为惰性元素，如Ne、Ar、Kr等）的基态最稳定构型。结果表明，团簇系统的能量存在多个极小点（鞍点），对应不同的能量，选取能量最小的，即给出了N=7，13，19时的最佳稳定构型。 关键词： 共轭梯度算法，零点算法，Lennard-Jones势，稳定构型 算法原理： (1) Lennard-Jones势的团簇能量 Lennard-Jones势一般适用于惰性气体原子间的Van der Waals键，即 1 其中，和是Weyl系数，对具体的元素是常数。 对势能求一阶导数，令其等于0得， 1式可简化为 2 此处，为计算方便，令得 3 求导得 4 团簇系统的Hamiltonian的形式为