Saturday, December 28, 2013

wavefunction01 原子核对电子作用的库仑势,当rR 时,库仑势可表示-7/r

2010 年11 月四川大学学报自然科学版)

第47 卷第6 期Journal of Sichuan University (Natural Science Edition)



doi: 103969/j. issn. 0490-6756. 2010. 06. 032

用类氟原子波函数对~ 团簇的

结构与能量的研究

赵茂娟1 , 2 朱博严永刚朱俊1

(1.四川大学物理科学与技术学院,成都610064;

2. 成都大学电子信息工程学院,成都610106)

Nov. 2010

Vol. 47 No.6



摘要:在单中心球模型近似下,选用类氟原子解析波函数,用变分法计算了m 团簇正八

面中心结构与能量.结果表明,当中心氢原子核到顶角氢原子核之间的距离 = 3.16α。时,

体系能量有一极小值=-4.153ho C α。=0. 529177X10-10 m , ho =27.2 eV) ,这表明H~ 团簇正八面中心体结构是稳定的结构, H~ 团菜是可能存在的.



关键词~ 团簇;羊中心球模型近似;变分法

中图分类号: 0469 文献标识码A 文章编号: 0490-6756(2010)06-1378-03



Calculation of the energy curve of body-centered regular octahedral structure

。f H~- cluster by one-center spherical symmetry wave f unction


ZHAOMao-]uan1

•2 , ZHU B01 , YAN Yong-Gang 1 , ZHU ]un 1
(1. College of Physics Science and Technology , Sichuan University, Chengdu 610065 , China;

2. College of Electronic and Information Engineering , Chengdu University , Chengdu 610106 , China)

Absìracì: The simple fluorine wave functions are chosen as the electronic wave functions of the H~- cluster



with one center spherical model approximation. The energy curve of the body-centered regular octa

hedral structure of the H~- cluster is calculated by means of variational method. The result of the calcu

lation shows that the cure has a minimal energy -4. 153ho at R=3. 16 αo. This means that the body-cen

tered regular octahedral structure of the H~- cluster may exist stably.

Key words: H~- cluster , one center spherical model approximation , variational method



1 引

对原子团簇的结构与能量的研究,多年来,我

们探索和总结出了计算量小、物理思想明确且易于

推广的两种不同的方法即改进了的排列通道量子

力学方法CMACQM) 和球对称模型下的变分法.我

们先后对原子团簇Li3' He;j、Het 、Hei 、Het 、

Het 、H';- 、Ht 、H g , H7 , H4 , Ht 、Hf 、H g H7 、

收稿日期: 2010-01-12



Hi+ 、H 俨、Hi+ 、的结构与能量进行了系统的理论研究[1 叫,理论计算结果与实验结果符合得较好,

表明我们采用的上述两方法是有效可靠的.本文用

球对称模型,选用类氟原子解析波函数,用变分法

对町的正八面体中心结构与能量进行了变分计算,得到了H? 团簇可能以以键长=3.16冉的



正八面体中心结构稳定存在,其能量值为=-4.

153 hoCao=O. 529177X10-10 m , ho=27. 2 eV).



基金项目:国家科技支撑计划重点项目(2007BAE13BOO) ;国家自然科学基金项目(1 0576020)

作者简介:赵茂娟(1 963 一) ,女,副教授,四川大学访问学者,主要从事原子与分子团簇研究

通讯作者·朱博. E-mail:zhubosw@gmail. com

第6 期赵茂娟等:用类氟原子泼函数对l 团簇的结构与能量的研究1379

Hl 正八面体中心结构Fig. 1 The body-centered regular octahedral structure

of the Hï- Cluster

2 理论与方法



H~ 团簇正八面体中心结构如图 所示,可

视中心的氢核在球心,距该中心为 且位于正八

面体顶角的其余6 个氢核均匀分布在以R 为半径

的球面上, 9 个电子在这种结构中运动在球内,电子受到7 个氢核的作用可视为 个在球心处的氢核电荷加上球面上的6 个氢核电荷的作用.在球



外,可视 个氢核都集中在球心处的作用.因此各

原子核对电子作用的库仑势,当r<R 时,库仑势可表示为一l/r - 6/R , 当r二?::R 时,库仑势可表示



为/r. 在此模型近似下,町团簇的电子运动可

用类氟原子波函数来近似描述.

根据以上分析,在单中心球模型近似下,可把

H? 团簇体系的哈密顿算符H 写为



H=一斗三:听- ~(于十~6 )十

'"飞 i f> I

三;于十二字且(1)

i<j f ij 卢 、4



上式中的第二项表示各电子与核之间的库仑势之

和;当r>R 时, r>=ri' 当r~;;R 时, r>=R. 上式的

第三项表示各电子之间的排斥势之和,最后一项表

示各核之间的排斥势之和.

H 2 的电子波函数可近似地用类氟原子波函



数:

φ = (SSS 头 兴7 PxPxPyPyPz) (2)



取S , S 兴,且,尺, Pz 为文献[3J 中的解析波函数

S. 伊1 (r) = N1 e 阳S 头伊2 (r) = N汇μr)e 卢e 阳

Px: 伊3(r) = N3 (卢cos8e-pr ~ (3)

Py :CP4 (r) = N4 (p )s in8e 甲卢Pz :cp5(r) = N5 (μr )sin8e 甲卢其中μ 、"为变分参数, N 为正交化系数, N1 、风、



风、陀、凡为归一化常数.

上述单电子波函数乘以自旋函数,并考虑到多

电子体系总波函数是反对称的,应取为later 行列

式形式:

φDet 1 伊1 (1 )ψ1 (2)cp2 (3) ρ( 的机(5)

件(6) 伊4 (7)伊4(8) 伊5(9) 1 (4)



其中伊z 与伊z 表示自旋相反波函数.

根据变分原理,团体系的能量可由下式确定:

(5)

若引人符号:

T i =- ~ f cpt V2 cpi dτ ,

V i = f伊t (二十去) CPi

V ij = f Jψ(r1)77(rJ|2dhdτ2

_ r t (r1 )cp; (r2) ψ;C r2 'j (r1) 't'l ,. 1 ./ 't'J ,.;:;./ "j/l




"


;:: ./ 't'J " 1./ d Ll dτ2 '1 J r12

(5) 式的能量E 可以展开写成:

E= 红卫十T2 ) 十5T3 (2V1 2V2

5V3) + V ll + V22 V33 V44 2(2V12

5V13 5V23 3V34 ) - 2X12 - 5 (X13

X 23 + X 34 ) + V NN (6)

其中VNN 为各核之间的排斥势之和.

3 计算结果与讨论(6) 式可知町的能量E 是参变量α.μ



键长 的函数即:

E = E(α μ R) (7)

应用变分原理于(7) 式,可求出平衡结构的键长R

与能量E 的数值用数值法先确定α μ 值,然后

计算出E R 变化的关系.当选取α= 4.98 μ=

3. 95 时,计算出E 随R 变化关系的数值列于附表

中.由此可得到当R=3.16矶时,能量有一极小值E=-4.153ho . 这一结果表明H~ 团簇可能以键R=3.16 α口的正八面体中心结构稳定存在.在

1380 四川大学学报自然科学版) 第47 卷

表1 酣一八四面体中心结构与能量计算值Tab. 1 The energy calculation value of the body- centered regular octahedral structure of the Hj- cluster




R/ao E/ho R/ao E/ho


2.30 0.229 3.00 4.109

2.40 l. 158 3.10 4.129

2.50 l. 715 3.11 4.135



2.60 2.433 3.12 4.141

2. 70 2.946 3.13 4.145

2.80 3.325 3.14 4.148

2.90 3.728 3.15 4. 151

其它文献中,我们没有查阅到相关或类似的理论和

实验结果.

参考文献:

[1 J Yang J , Li p , Gou Q Q. Formation mechanism and



the binding energy of the body-centred regular tetra

hedral structure of Li3 [J]. Commun Theor Phys ,

2005 , 44 (3): 525.

[2J Zhang J p , Gou Q Q , Li P. Formation mechanism



and binding energy for icosahedral central structure of

He;t; cluster [j]. Commun Theor Phys , 2004 , 42



(5): 763.

[3J Zhang J p , Gou Q Q , Li P. Formation mechanism



and binding energy for body-centered cubic structure

of Het Ciuster[J]. Commun Theor Phys , 2004 , 41

(1): 10 1.

[4J Zhang J p , Gou Q Q , Li P. The formation mechanism



and binding energy for the octahedral central

structure of the Hei cluster [J]. Commun Theor

Phys , 2003 , 40 (2): 203.

[5J Li p , Xiong y , Gou Q Q , et al. The formation



mechanism and the binding energy of the body-cen

tred regular tetrahedral structure of Het [J]. Chinese

Physics , 2002 , 11 (10): 1018.

[6J Gou Q Q , Zhang J p , Li P. Formation mechanism



and binding energy for equilateral triangle structure of

Het cluster [j]. Commun Theor Phys , 2002 , 38



(5): 597.

[7J Li p , Gou Q Q , Zhang J p , etal. The MACQM calculation



for the binding energy of the equilateral tri

angle structure of H3 cluster [J J. Commun Theor

Phys , 2001 , 35 (3): 327.

[8J Li p , Yin H , Gou Q Q , etal. The MACQM calcula



tion of the total energy curve for the equilateral trian

gle structure of Ht cluster [j]. Commun Theor




R/ao E/ho R/ao E/ho


3.16 4. 153 3.50 2.901

3.17 4. 150 3.60 2. 793

3.18 4. 146 3. 70 2.602

3.19 4. 141 3.80 2.523

3.20 3. 137 3.90 2.420

3.30 3. 121 4.00 2.318

3.40 3. 102 4.10 2.256

Phys , 1999 , 32 (3): 467.

[9J 曹乐千,于开宣,朱俊.用类 原子波函数对比团簇的体心立方结构与能量的研究[J].四川大学学

报自然科学版, 2003 , 40(3) :485.

[10J 朱俊,卿勇.用单中心球对称波函数对,团簇正八



面体中心结构与能量的计算[J].四川大学学报自

然科学版, 1998 , 35( 的: 854.

[11J 蒲利春,朱俊,张雪峰.比氢原子团簇结构与能量



的计算[JJ.四川大学学报自然科学版, 2005 , 42

(5) :993.

[12J 朱俊,曹乐千.用类。原子波函数对t 团簇的体



心立方结构与能量的研究[J].原子与分子物理学

报, 200 1, 18(2) :238.

[13J 陈纲,朱俊.用类碳原子波函数对i 团簇的正八



面体中心结构与能量的研究[J].四川大学学报自

然科学版, 2005 , 42(2) :340.

[14J 苟清泉,朱俊.用单中心球对称波函数计算H g 团簇体心立方结构的能量曲线JJ.原子与分子物理

学报, 1998 , 15(3) :288.

[15J 朱俊,卿勇,苟清泉.用类氧原子波函数计算H,团簇正八面体中心结构的能量曲线JJ.原子与分子

物理学报, 1999 , 16( 1) :69.

[16J 汪志刚,朱俊,陈艳用类i 原子波函数对l+团簇的正四面体中心结构与能量的研究[J].原子与



分子物理学报, 2007 , 24(2):303.

[17J 汪志刚,朱俊,程艳用类 原子波函数对四十团簇的正八面体中心结构与能量的研究[JJ.原子与



分子物理学报, 2006 , 23(4):734.

[18J 汪志刚,朱俊,陈刚用类 原子波函数对;+ 团簇的体心立方结构与能量的研究[J].四川大学学

报自然科学版, 2006 , 43 (4) : 852. .

[19J 苟清泉,黄树勋原子结构与变分计算[M]. 成都:



成都科技大学出版社, 1989.

[责任编辑陈向荣

No comments:

Post a Comment