APP下载

一种新的双电子积分的计算方法

2019-04-03辽宁省锦州中学赵宇阔

数学大世界 2019年4期
关键词:计算精度变分计算方法

辽宁省锦州中学 赵宇阔

一般来说,在单中心的双电子积分的计算中,采用如下的近似方式是一种比较常见的解决方案,即:

但是,在这种展开方式中存在着一种截断误差,特别是低轨道的时候,它的展开精度相对比较低。例如,在1s 轨道的计算中,它的计算精度为1 阶,即:。因此,除了Hylleraas 变分方法之外,其他的计算方法若不经过试验数据的拟合,其总能量的相对误差都将大于1.8%。

因此,为了解决这个问题,本文提出了一种新的计算方法,并且可以向任意阶展开。所以,本文的计算精度将优于上面的展开方式。

一、基本原理

(1)。

那么根据式(1),可得:

那么,根据式(2)和麦克劳林级数,可得:

那么,根据式(3),可得:

其中:

在定积分的计算中,由于高斯-勒让德积分可以向任意阶展开,所以,本文先通过积分变换,把一个六维的双电子积分简化为一个五维的定积分问题,然后再根据泰勒展开和高斯-勒让德积分的特点,从而获得式(4)的积分形式。并且,相对于文献1 来说,本文的计算精度要更高一些,但是时间复杂度却高于文献1。因此,在双电子积分的计算中,若是低轨道的积分计算,则建议采用本文的方法。否则,建议采用文献1 的方法。同时,有关这两种方法的结合使用,将是本文后续关注的一个焦点问题。

猜你喜欢

计算精度变分计算方法
求解变分不等式和不动点问题的公共元的修正次梯度外梯度算法
槽道侧推水动力计算方法研究
浮力计算方法汇集
极限的计算方法研究
基于SHIPFLOW软件的某集装箱船的阻力计算分析
自反巴拿赫空间中方向扰动的广义混合变分不等式的可解性
基于变分水平集方法的数字图像分割研究
钢箱计算失效应变的冲击试验
一种伺服机构刚度计算方法
具有(H,η)单调算子的广义非线性变分包含