双中心积分计算氟化氢HF的动能和势能

<F 1s | H 1s>,<F 2s | H 1s>,<F 2p0 | H 1s>,氟化氢中有这三种情况是双中心积分，

列成表格

	n	l	n'	l'	λ	ζ	ζ'	R	Pople John A	计算值	参考值/计算值
<F 1s \| H 1s>	1	0	1	0	0	8.7	1	1.733	-2.1798	-2.179823778	0.999989092
<F 2s \| H 1s>	2	0	1	0	0	2.6	1	1.733	-4.3051	-4.305097499	1.000000581
<F 2p0 \| H 1s>	2	1	1	0	0	2.6	1	1.733	-2.2822	-2.282193291	1.00000294

利用前面得到Tab，VabA和VabB方法得到一个新的方法


public static  double cou( double R ,int zn1,  int zn2,double za,double zb, int n1,int L1 ,int m1 ,int n2 ,int L2 ,int m2  ) throws IOException, ParseException, InterruptedException {double d= Tab (R , za , zb ,   n1,  L1 , m1 , n2 , L2 , m2  );//0.00523319906316793double d1= VabA( R  , zn1 , za , zb ,   n1,  L1 , m1 , n2 , L2 , m2 );//-2.152403725820715 double d2= VabB( R , zn2 , za , zb ,   n1,  L1 , m1 , n2 , L2 , m2 );//-0.03265325106338999double f=d+d1+d2;return f;   }

直接代入这个新方法就可以得到这三个积分，

double d=cou( 1.733 ,9,1, 8.7 , 1 , 1, 0 , 0 , 1 , 0 , 0 ); //-2.1798237778209373

double d=cou( 1.733 ,9,1, 2.6 , 1 , 2, 0 , 0 , 1 , 0 , 0 ); //-4.305097498533459

double d=cou( 1.733 ,9,1, 2.6 , 1 , 2, 1 , 0 , 1 , 0 , 0 );//-2.282193290607004

完整的java程序

package greenwood;import java.io.IOException;
import java.text.ParseException;
import udp.HFoverlapA;public class kin {public static  double FACT( double n ) throws IOException, ParseException {double prodt=1.0;for(int a=1 ;a<n+1 ;a++){prodt=prodt*a;}if(n<0){System.out.println( prodt+   "  负数阶乘 "  );}return prodt;}public static  double Tab( double R,double za,double zb, int n1,int L1 ,int m1 ,int n2 ,int L2 ,int m2 ) throws IOException, ParseException, InterruptedException {double f=0;double f1=0;double f2=0;double f3=0;f1= ( -zb*zb/2)*HFoverlapA.Sab ( R,za, zb,  n1,  L1 ,  m1 ,  n2 ,  L2  );System.out.println(  f1 +" f1 Tab "  );double d1=( ( 2*n2*zb*zb) /Math.pow ( 2*n2*(2*n2-1) ,0.5 ) );double d2=HFoverlapA.Sab (R, za, zb,  n1,  L1 ,  m1 ,  (n2-1) ,  L2  );f2=d1*d2;System.out.println(  f2 +" f2 Tab   " +d1+"  "+d2   );     double d5=0;double d6=0;if(n2>1){double d4=2*zb*zb*(L2*(L2+1)-n2*(n2-1) );if(Math.abs(d4)>1e-6){d5=Math.pow ( 2*n2*(2*n2-1)*(2*n2-2)*(2*n2-3) ,0.5);d6= HFoverlapA.Sab (R, za, zb,  n1,  L1 ,  m1 ,  (n2-2) ,  L2  );f3=(d4/d5)*d6;}System.out.println(  f3 +" f3 Tab " +d4+"  "+d5+"  "+d6 );}f=f1+f2+f3;System.out.println( f+"  "+ f1+"  "+f2+"  "+f3+" f Tab "  );return f;}public static  double VabA( double R ,int zn, double za,double zb, int n1,int L1 ,int m1 ,int n2 ,int L2 ,int m2 ) throws IOException, ParseException, InterruptedException {double f1=-zn* (2*za/ Math.pow(2*n1*(2*n1-1),0.5 ) )*HFoverlapA .Sab ( R, za, zb ,    n1-1, L1 , m1 ,n2 ,L2 ); return f1;}public static  double VabB( double R ,int zn, double za,double zb, int n1,int L1 ,int m1 ,int n2 ,int L2 ,int m2 ) throws IOException, ParseException, InterruptedException {double f1=-zn* (2*zb/ Math.pow(2*n2*(2*n2-1),0.5 ) )*HFoverlapA .Sab ( R, za, zb ,    n1, L1 , m1 ,n2-1 ,L2 ); return f1;}public static  double cou( double R ,int zn1,  int zn2,double za,double zb, int n1,int L1 ,int m1 ,int n2 ,int L2 ,int m2  ) throws IOException, ParseException, InterruptedException {double d= Tab (R , za , zb ,   n1,  L1 , m1 , n2 , L2 , m2  );//0.00523319906316793double d1= VabA( R  , zn1 , za , zb ,   n1,  L1 , m1 , n2 , L2 , m2 );//-2.152403725820715 double d2= VabB( R , zn2 , za , zb ,   n1,  L1 , m1 , n2 , L2 , m2 );//-0.03265325106338999double f=d+d1+d2;return f;   }public static  void res(   ) throws IOException, ParseException, InterruptedException {// double d=cou( 1.733 ,9,1, 8.7 , 1 ,   1,  0 , 0 , 1 , 0 , 0 ); //-2.1798237778209373// double d=cou( 1.733 ,9,1, 2.6 , 1 ,   2,  0 , 0 , 1 , 0 , 0 ); //-4.305097498533459 double d=cou( 1.733 ,9,1, 2.6 , 1 ,   2,  1 , 0 , 1 , 0 , 0 );//-2.282193290607004System.out.println(   d+"     res "   );}public static void main(String[] args) throws IOException, ParseException, InterruptedException {res(   );}}

双中心积分计算氟化氢HF的动能和势能相关推荐

双中心单电子积分计算氢分子离子H2+的轨道能量
双中心用椭圆坐标用线性组合原子轨道LCAO方法把电子波函数设为假设这个轨道是两个轨道的线性组合,φA,φB用氢原子的1s轨道代入轨道能量E 将E的表达式展开因为cA,cB是实系数,因此可以进一 ...
浅谈同城双中心的网络部署模型
企业建设数据中心时,出于灾备的考虑,会建设两个甚至多个数据中心.例如我们经常提到的"两地三中心",即同城双中心+异地中心. 同城双中心是指在同城或邻近城市建立两个可独立承担业务的数 ...
ES双中心数据稽核（同步）
数据稽核场景基于ES双中心使用场景,需要近实时周期性验证2个ES索引双写的数据差异,并采取措施保证数据的一致性. 数据稽核方案由于ES一般承载的数据量较大,无法直接查询到内存中进行明细比对,采用的 ...
STO双中心势能积分
势能积分的计算公式 *2 由公式得到两个java计算方法 VabA( double R ,int zn, double za,double zb, int n1,int L1 ,int m1 ,int ...
验算双中心重叠积分程序
用前面得到的计算勒让德多项式系数的新方法改写了计算重叠积分的程序,这次验证这个程序的性能,主要参考了 1.Approximate Molecular Orbital Theory by Pople J ...
计算约化重叠积分氟化氢HF
计算F原子的(2,1,0)电子和H原子(1,0,0)电子的重叠积分 F原子核外电子的角量子数不为0,这个重叠积分的计算方法为其中ζa和ζb是斯莱特规则下的核电荷数,R是两个原子核距离为1.733,α ...
LCAOSCF自洽场氟化氢HF斯莱特函数
设椭圆坐标得到 F原子的轨道是 H原子的轨道是因为和斯莱特规则则 z1=8.7 z2=2.6 z3=1 R=1.733 F原子的斯莱特轨道是 fx1=(z1)**(1.5)*sympy.exp ...
“双碳政策”下关于数据中心网络可视化管理的一些思考
站在新一年的起点,回顾2021全年,疫情持续下全球经济动荡,国内新基建大潮方兴未艾,双碳与能耗双控牵动千行百业,数据中心产业迎来转折,由快速发展向规范发展过渡.回顾数据中心行业近几年的发展方向,我们可 ...
盘点2021：数据中心行业迎来转折点
站在新一年的起点,回顾2021全年,疫情持续下全球经济动荡,国内新基建大潮方兴未艾,双碳与"能耗双控"牵动千行百业,数据中心产业迎来转折,由快速发展向规范发展过渡. 回望数据中心行 ...

双中心积分计算氟化氢HF的动能和势能

双中心积分计算氟化氢HF的动能和势能相关推荐

最新文章

热门文章