STO双中心势能积分

势能积分的计算公式

由公式得到两个java计算方法

VabA( double R ,int zn, double za,double zb, int n1,int L1 ,int m1 ,int n2 ,int L2 ,int m2 )

VabB( double R ,int zn, double za,double zb, int n1,int L1 ,int m1 ,int n2 ,int L2 ,int m2 )

zn是原子核的核电荷数，za，zb分别是ζ，ζ‘。公式里是假设两个原子的核电荷数都是1。

验算这两个方法，参考了两篇文章

1.Calculation of the One-Electron Two-Center Integrals Over Slater-Type Orbitals by Means of the Ellipsoidal Coordinates Method by SIDI MOHAMED MEKELLECHE, ABDELLATIF BABA-AHMED（1996）

2..Evaluation of Two-center One- and Two-electron Integrals over Slater Type Orbitals by YAKAR, Yusuf , ÖZMEN, Ayhan，ATAV, Ülfet （2006）

得到表格

	n	l	n'	l'	λ	ζ	ζ'	R	YAKAR, Yusuf	计算值	计算值/参照值
zna=1	2	1	3	2	1	4.2	4.4	2.5	0.011498097921007800	0.011498097921203000	1.000000000016980000
	3	2	2	1	0	2.3	2.3	2.5	0.154172007331042000	0.154172007331022000	0.999999999999870000
	3	2	3	2	2	5.8	4.5	0.2	-1.471238017764440000	-1.471238018706040000	1.000000000640010000
	3	2	3	2	1	2	4.4	8.7	0.000000857400767270	0.000000857431091410	1.000035367522470000
	4	3	5	2	1	1.5	0.5	15	0.002572464174355450	0.002572464177822200	1.000000001347640000
	5	3	6	4	3	1.5	1	0.1	0.006949396708185720	-270.000000000000000000
	4	3	8	5	2	3.6	3.7	0.1	-0.007197120074513900		0.000000000000000000
	6	4	12	4	4	3.5	0.5	20	-0.000003811906504442		0.000000000000000000
	5	3	8	5	3	2.5	2	2.2	-0.182136029883179000		0.000000000000000000
	10	5	5	4	2	0.3	0.1	10	-0.006308043501843390		0.000000000000000000
	15	12	10	8	3	4	2.5	1.5	-0.034376370777494700		0.000000000000000000
	22	20	24	18	15	1.5	3.5	0.8	-0.001059509498481520		0.000000000000000000

									SIDI MOHAMED MEKELLECHE	计算值
zna=35	3	2	3	2	0	6.5197	2.0387	3.75803	-0.486277949158820000	-0.486277948407194000	0.999999998454328000
	3	2	3	2	1	6.5197	2.0387	3.75803	0.265448562634640000	0.265448562222838000	0.999999998448656000
	3	2	3	2	2	6.5197	2.0387	3.75803	-0.041387530907970000	-0.041387530988898000	1.000000001955370000

znb=17	3	2	3	2	0	6.5197	2.0387	3.75803	-0.086569711360165000	-0.086569711326935500	0.999999999616153000
	3	2	3	2	1	6.5197	2.0387	3.75803	0.039704646900045000	0.039704646872093500	0.999999999296014000
	3	2	3	2	2	6.5197	2.0387	3.75803	-0.005394666015387200	-0.005394666000166480	0.999999997178561000

结果表明用这个程序计算势能积分，主量子数不能大于5.

代入比如第一组

double d2= VabA( 2.5 , 1 , 4.2, 4.4 , 2 , 1 , 1 , 3 , 2 , 1 ); // 0.011498097921203

最后一组

double d2= VabB( 3.75803 , 17 , 6.5197 , 2.0387 , 3 , 2 , 2 , 3 , 2 , 2 );// -0.00539466600016648

对应的java程序为

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;}

STO双中心势能积分相关推荐

双中心积分计算氟化氢HF的动能和势能
<F 1s | H 1s>,<F 2s | H 1s>,<F 2p0 | H 1s>,氟化氢中有这三种情况是双中心积分, 列成表格 n l n' l' λ ζ ζ' ...
双中心单电子积分计算氢分子离子H2+的轨道能量
双中心用椭圆坐标用线性组合原子轨道LCAO方法把电子波函数设为假设这个轨道是两个轨道的线性组合,φA,φB用氢原子的1s轨道代入轨道能量E 将E的表达式展开因为cA,cB是实系数,因此可以进一 ...
浅谈同城双中心的网络部署模型
企业建设数据中心时,出于灾备的考虑,会建设两个甚至多个数据中心.例如我们经常提到的"两地三中心",即同城双中心+异地中心. 同城双中心是指在同城或邻近城市建立两个可独立承担业务的数 ...
ES双中心数据稽核（同步）
数据稽核场景基于ES双中心使用场景,需要近实时周期性验证2个ES索引双写的数据差异,并采取措施保证数据的一致性. 数据稽核方案由于ES一般承载的数据量较大,无法直接查询到内存中进行明细比对,采用的 ...
整流十—双二阶广义积分锁相
前言前面<整流专题>中描述了电网平衡下的 1.锁相环-单同步锁相环(获得电网基波旋转角度) 以及电网不平衡下情况下的锁相环. 2.锁相环-双同步坐标系锁相(获得电网基本正序分量旋转角度) ...
智慧党建/党员双报到/党建积分系统/智慧党建小程序/党建云平台
功能亮点 1. 基层问题得到解决 2. 党员.单位组织建立双报到机制 3. 基础群众反馈有渠道 4. 党建资源的开放实现效果 1. 提升基层治理水平,提高党员干部能力素质,满足人民群众的需求 2. ...
验算双中心重叠积分程序
用前面得到的计算勒让德多项式系数的新方法改写了计算重叠积分的程序,这次验证这个程序的性能,主要参考了 1.Approximate Molecular Orbital Theory by Pople J ...
[会员积分运营了解]各大主流电商平台会员及积分体系概况集合！
一.淘宝·天猫 (一)会员名称: 普通会员:淘气值>=300且<1000分: 超级会员:淘气值>=1000分: APASS会员:是属于主动邀请制,不支持任何形式的会员主动申请,邀请朋 ...
前阿里巴巴产品专家明卿：当下做通用积分还是个好生意吗？
1. 准备建立一个垂直社群,致力于专家与群友的主题分享与线下沙龙.请问在初期需要建立严格的会员系统吗? 问题描述:垂直社群,以学习和分享为初期发展目标.在社群建立初期是否应该建立严格的会员系统,还是只 ...

STO双中心势能积分

STO双中心势能积分相关推荐

最新文章

热门文章