计算马德隆常数
- 格式:docx
- 大小:51.95 KB
- 文档页数:6
二维三角离子晶体马德隆常数的计算一、概述二维三角离子晶体的研究在固体物理学领域具有重要意义,其中马德隆常数是一个关键的物理量。
马德隆常数是用来描述晶体中离子间相互作用的强度和类型的参数,对于研究晶体的热力学性质、电子结构和声学性能等方面都具有重要意义。
本文将对二维三角离子晶体马德隆常数的计算进行深入探讨,以期为读者提供有价值的信息。
二、二维三角离子晶体的结构和性质在二维三角离子晶体中,正负电荷的离子按照一定的规则排列在平面上,形成一种特殊的晶体结构。
这种结构具有较强的几何规则性和周期性,因此在研究中具有一定的优势。
二维三角离子晶体的性质受到离子间相互作用的影响很大,而马德隆常数则是描述这种相互作用的重要参数。
三、马德隆常数的定义和物理意义马德隆常数是由英国物理学家马德隆在20世纪初提出的,它是描述晶体中离子间相互作用的参数。
马德隆常数的大小和符号决定了晶体的稳定性和各种物理性质,因此在材料科学研究中具有重要意义。
对于二维三角离子晶体来说,计算马德隆常数可以帮助我们更好地理解晶体结构和性质。
四、马德隆常数的计算方法在实际研究中,计算马德隆常数通常需要借助于复杂的计算方法和理论模型。
一般来说,我们可以通过经典力学或量子力学的理论来建立晶格模型,并在此基础上进行数值计算。
对于二维三角离子晶体来说,由于结构的特殊性,计算方法可能会更加复杂和繁琐。
五、个人观点和理解在我看来,二维三角离子晶体马德隆常数的计算是一个具有挑战性但又具有重要意义的课题。
通过对马德隆常数的准确计算,我们可以更好地理解晶体的结构和性质,为材料科学研究提供重要的参考。
未来,我希望能够进一步深入研究这一课题,为相关领域的发展贡献自己的力量。
六、总结和回顾通过本文的阐述,我们对二维三角离子晶体马德隆常数的计算有了更深入的了解。
马德隆常数作为描述晶体中离子间相互作用的重要参数,对于研究晶体的性质具有重要意义。
未来,我们需要进一步深入研究这一领域,为材料科学的发展做出更大的贡献。
计算马德隆常数《固体物理》马德隆常数的计算学院:物理学院学号:2011012643姓名:刘娴雅马德隆常数的计算摘要:通过分析马德隆常数的三种计算方法和其相应的使用范围,得出不同晶体结构下相应的计算方法和使用范围.关键字:马德隆常数离子晶体在固体物理学中,当计算离子晶体的结合能时,需知道马德隆常数的值, 因此,马德隆常数在离子晶体的理论研究和科学实验中占有十分重要的地位.该值一般由实验确定。
马德隆常数是描述离子晶体结构的常数,其定义公式为:n1、n2、n3为离子晶体中任一离子相对于中心离子的坐标,?为求和遍及晶体中所有离子。
由于离子晶体为数目巨大的多粒子系统,因此马德隆常数一般情况下由实验确定。
nnn123,,,,,1,,,, 222nnnn,n,n123123离子晶体结合的性质比较简单,在近代微观理论发展初期,计算离子晶体的结合能获得很好的结果,对于验证理论起到了重要作用,所用的方法和概念在处+-以NaCl为例,由于Na和Cl都是满壳层的结构,具有球对理许多问题中还常用到.称性,考虑库仑作用时,可以看做点电荷.先考虑一个正离子的平均库仑能.如果令r表示相邻离子的距离,该能量可表示为nnn1232,,1q,1,,, (1) 2222222nnn4,,(nr,nr,nr)123012312222222(nr,nr,nr)如果以所考虑的正离子为原点, 可以表示其他各123 离子所占格点的距离,并且对于所有负离子格点,n1+n2+n3=奇数,所有正离子格n1+n2+n3点,n1+n2+n3=偶数.考虑到正负离子电荷的差别,引入因子(-1),一个原胞的能量为,,nnn22123,,,q,1q,,, (2) 222222,,4r4,,rnnn00(nr,nr,nr)123123 nnn123,,,,,1,,, (3) ,222nnnn,n,n123123α为一无量纲的数,完全决定于晶体结构,称之为马德隆常数.在具体计算中发现,求和时既有正项,又有负项,如果逐项相加,并不能得到收敛的结果.对于一维情况,其级数求和很容易计算,如两种一价离子组成的一维晶格的马德隆常数,利用定义很容易计算出α=2ln2,但对于三维情况,其级数收敛很慢. 1918年Madelung首先计算这种级数和,他先将晶体中点阵视为一系列中性平面点阵组成即该平面内点阵由一系列中性直线点阵组成,其上正负电荷相等且按格点周期分布.由此将电势展开成傅里叶级数并用了享克尔函数(Hankel function),进而求出马德隆常数.这种方法对于计算像氯化钠那样简单的离子晶体取得了成功.但对大多数离子晶体而言并不适用。
K3C60的马德隆常数计算081120… riverman一、理论分析当选取晶胞较小时,非中性外表面的面电荷对中心离子的马德隆常数影响较大。
外表面类型相同的晶胞马德隆常数收敛,而外表面类型不同的晶胞马德龙常数相差较大。
应将带不同类型面电荷的晶胞的马德隆常数取平均,以快速计算较大晶胞的马德隆常数。
对于K3C60,当计算K(1)离子的马德隆常数时,应将外表面为C+K(1)、K(2)、K(1)+C、K(2)的四个连续扩大的晶胞的马德龙常数取平均;计算C离子的马德隆常数同K(1);计算K(2)的马德隆常数时,应将外表面为C+K(1)和K(2)的马德隆常数取平均,以计算无限大晶体的马德隆常数。
二、K(1)离子的马德隆常数:程序运行结果:程序代码:#include "stdio.h"#include "math.h"#define R 4int convertflag(long x) /*K(1)与C原离子变换*/ {int flag;if(x%2==1) flag=1;else flag=-3;return flag;}void main(){int n,m,l,i,j,k,flag,n1,n2,n3; /*m为胞长,l半胞长,flag为K(1)和C的电荷数*/long x;/*最外层边长m>=4,正半轴长l>=2*/double r,a=0,ratio=1,b[100]; /*a为马德隆常数*/printf("max length of cell n=");scanf("%d",&n);printf("note:m is the length of cell\tset crystal constant=4\n");for(l=2;l<=n/2;l++){m=2*l;if(l%2==0) /*最外层C+K(1)*/ {/*计算与C+K(1)相互作用*/if(l%4==0) /*最外层顶点同号*/{flag=1; /*顶点电荷*/x=0; /*顶点翻转*/}else /*最外层顶点异号*/{flag=-3; /*顶点电荷*/x=1; /*顶点翻转*/}for(i=-l;i<=l;){n1=i;i+=2;for(j=-l;j<=l;){n2=j;j+=2;for(k=-l;k<=l;){n3=k;k+=2;if(n1*n1+n2*n2+n3*n3==0){flag=convertflag(x++);continue;}r=sqrt(n1*n1+n2*n2+n3*n3);if(n1==l||n1==-l) ratio=ratio/2;if(n2==l||n2==-l) ratio=ratio/2;if(n3==l||n3==-l) ratio=ratio/2;a+=-flag/r*ratio*R;/*printf("(%d,%d,%d)\tflag=%d\tratio=%f\ta=%f\n",n1,n2,n3,flag,ratio,a);*/flag=convertflag(x++);ratio=1;}}}/*计算与K(2)相互作用*/for(i=-l+1;i<=l-1;){n1=i;i+=2;for(j=-l+1;j<=l-1;){n2=j;j+=2;for(k=-l+1;k<=l-1;){n3=k;k+=2;r=sqrt(n1*n1+n2*n2+n3*n3);a+=-1/r*R;/*printf("(%d,%d,%d)\ta=%f\n",n1,n2,n3,a);*/}}}}else /*最外层K(2)*/ {/*与计算K(2)相互作用*/for(i=-l;i<=l;){n1=i;i+=2;for(j=-l;j<=l;){n2=j;j+=2;for(k=-l;k<=l;){n3=k;k+=2;r=sqrt(n1*n1+n2*n2+n3*n3);if(n1==l||n1==-l) ratio=ratio/2;if(n2==l||n2==-l) ratio=ratio/2;if(n3==l||n3==-l) ratio=ratio/2;a+=-1/r*ratio*R;/*printf("(%d,%d,%d)\tratio=%f\ta=%f\n",n1,n2,n3,ratio,a);*/ratio=1;}}}/*计算与C+K(1)相互作用*/if((l-1)%4==0) /*次外层顶点同号*/ {flag=1;x=0;}else /*次外层顶点异号*/ {flag=-3;x=1;}for(i=-l+1;i<=l-1;){n1=i;i+=2;for(j=-l+1;j<=l-1;){n2=j;j+=2;for(k=-l+1;k<=l-1;){n3=k;k+=2;if(n1*n1+n2*n2+n3*n3==0){flag=convertflag(x++);continue;}r=sqrt(n1*n1+n2*n2+n3*n3);a+=-flag/r*R;/*printf("(%d,%d,%d)\tflag=%d\ta=%f\n",n1,n2,n3,flag,a);*/flag=convertflag(x++);ratio=1;}}}}…………若干}getch();}三、C离子的马德隆常数程序运行结果:程序代码:程序与计算K(1)类似,省略相同部分,划线部分为代码差别之处。
高等固体物理作业题 目: 马德隆常数的计算方法及实例计算 学生姓名: 学 院:理学院 专 业:物理电子学 指导教师:2013 年 12 月 7日学校代码:10128 学 号:摘要在固体物理学中,当计算离子晶体的结合能、晶格能、表面能等时,需知道马德隆常数的值,该值一般由实验确定。
马德隆常数是描述离子晶体结构的常数,是晶体结构的一个重要的特征参数,为一无量纲的数,只取决于晶体结构,在离子晶体的研究中占有重要的地位。
本文概述了晶体马德隆常数的几种计算方法及其使用范围,并举例简述了一维离子链,二维正方离子格子,以及三维Nacl离子晶体实例的马德隆常数的计算方法。
关键词:离子晶体;马德隆常数;计算方法;实例AbstractIn solid state physics, when calculate the combined energy, attice energy, surface energy, etc. of the ionic crystals, we need to know the Madelung constant value, which is generally determined by experiment. Madelung constant is used to describe structure of ionic crystal. Madelung constant is an important feature of the crystal structure parameters. Madelung constant is a dimensionless number that only depends on the crystal structure, and plays an important role in the study of ionic crystals. This article outlines several crystal Madelung constant calculation methods and its scope of application, and an example calculation methods outlined Madelung constant one-dimensional ion chains, two-dimensional square lattice ions, as well as three-dimensional Nacl ionic crystals instance.Keywords: ionic crystals; Madelung constant; calculation methods; examples目录引言 (1)1 晶体马德隆常数的几种计算方法 (2)1.1 定义法 (2)1.2 Evjen晶胞法 (2)1.3 计算晶格静电能法 (3)1.4 小结 (4)2 马德隆常数的实例计算 (5)2.1 一维离子链的马德隆常数计算 (5)2.2 二维正方离子格子的马德隆常数计算 (6)2.3 三维离子晶体(Nacl)的马德隆常数计算 (7)参考文献 (10)引 言马德隆(Madelung)常数α是晶体结构中的一个重要的特征参数,是描述粒子晶体结构的常数。
马德隆常数计算公式
马德隆常数是一种重要的数学常数,常用于计算圆周率π的近似值。
这个常数由英国数学家约翰·马德隆(John Madelaune)于1706年推导出来,因此得名。
马德隆常数的计算公式相对简单,可以通过以下式子得出:M = ∑(4/((2n-1)(2n+1))), 其中n从1开始递增,一直累加到无穷大为止。
通过这个公式,我们可以逼近得到π的值。
然而,马德隆常数不同于其他常数的地方在于它的收敛速度非常缓慢。
换句话说,要通过计算马德隆常数来得到π的精确值,需要进行无数次的累加运算。
这使得计算过程非常繁琐,甚至对于现代计算机而言也是一个巨大的挑战。
尽管马德隆常数在计算圆周率方面的应用存在一定的困难,但它在数值计算和近似计算领域仍然有着重要的意义。
它可以作为一种近似计算π的方法,尤其适用于没有高精度计算设备的时代。
在缺乏现代计算工具的情况下,数学家们通过手工计算马德隆常数,得到了种种与π相关的近似值,为后来的研究和计算提供了重要的基础。
此外,马德隆常数还引发了数学上的一些有趣问题和研究方向。
例如,我们可以探究该常数的收敛速度,研究如何加速计算马德隆常数以及优化近似π的方法等等。
这些问题激发了众多数学家的兴趣,并催生了一系列关于马德隆常数的研究成果。
总而言之,马德隆常数是一个重要而有趣的数学常数,虽然计算过程繁琐但在近似计算π和数值计算中仍然有一定的应用价值。
它不仅是对历史上数学家们智慧的体现,也激发了更多关于近似计算和数值分析的研究。
固体物理》
马德隆常数的计算
学院:物理学院
学号:
2011012643 姓
名:刘娴雅
马德隆常数的计算
摘要:通过分析马德隆常数的三种计算方法和其相应的使用范围,得出不同晶 体结构下相应的计算方法和使用范围•
关键字:马德隆常数离子晶体
在固体物理学中,当计算离子晶体的结合能时,需知道马德隆常数的值,因 此,马德隆常数在离子晶体的理论研究和科学实验中占有十分重要的地位 .该值 一般由实验确定。
马德隆常数是描述离子晶体结构的常数,其定义公式为:n1、n2、 n3为离子晶体中任一离子相对于中心离子的坐标 ,刀为求和遍及晶体中所有离 子。
由于离子晶体为数目巨大的多粒子系统,因此马德隆常数一般情况下由实验 确定。
1 n 1 n
2 n
3 门1门2门3 电 n 12 n ; n :
离子晶体结合的性质比较简单,在近代微观理论发展初期,计算离子晶体的 结合能获得很好的结果,对于验证理论起到了重要作用,所用的方法和概念在处 理许多问题中还常用到.以NaCI 为例,由于NS 和Cl -都是满壳层的结构,具有球对 称性,考虑库仑作用时,可以看做点电荷•先考虑一个正离子的平均库仑能•如果 令r 表示相邻离子的距离,该能量可表示为
2 n i n 2 n 3
q i ,2 2 口 r
离子所占格点的距离,并且对于所有负离子格点,n1+ n2+门3=奇数,所有正离子格 点, n1+ n2+n3=禺数.考虑到正负离子电荷的差别,引入因子(-1) n1+n2+n3, 一个原胞的 能量为
a 为一无量纲的数,完全决定于晶体结构,称之为马德隆常数.在具体计算中 发现,求和时既有正项,又有负项,如果逐项相加,并不能得到收敛的结果.对于一 维情况,其级数求和很容易计算,如两种一价离子组成的一维晶格的马德隆常数, 利用定义很容易计算出a =2ln2,但对于三维情况,其级数收敛很慢.1918年 Madelung 首先计算这种级数和,他先将晶体中点阵视为一系列中性平面点阵组 成即(1)
2 2 n 2 r 如果以所考虑的正离子为原点,(m 1
n 3r 2)2
可以表示其他各 2
q 1 n n 2 r b 2 q 4 0r ngrb
.. (n :r 2 n|r 2 nfr 2)
4 o r (2)
(3)
n 3 n 1 n 2
该平面内点阵由一系列中性直线点阵组成,其上正负电荷相等且按格点周期分布.由此将电势展开成傅里叶级数并用了享克尔函数(Hankel function), 进而求出马德隆常数.这种方法对于计算像氯化钠那样简单的离子晶体取得了成功. 但对大多数离子晶体而言并不适用。
因为不能满足上述中性平面的条件.因此该方法仅有历史价值.
1932年,Evjen认为把(1)式级数中的各项合理安排使其正项和负项的贡献几乎互相抵消,使级数迅速收敛,由此提出了计算马德隆常数的方法,其基本思想是:把晶体看成是由Evjen晶胞构成,Evjen晶胞内所有离子的代数和为零,把这些中性晶胞对参考离子的库仑能量的贡献加起来,若离子在这个中性立方体的面上、棱上或角上,其贡献取1/2、1/4或1/8,进而计算马德隆常数.以NaCI晶体为例,采用Evjen 方法,其收敛速度为1.456,1.752,1.747,计算到第90个Evjen 晶胞时,其马德隆常数为1.747564595,可见其计算是精确的.采用Evjen晶胞方法计算NaCl晶体马德隆常数,是一个很成功的例子,但对CsCl晶体结构,当Evjen晶胞最外层离子与参考离子同号时计算的马德隆常数,与当Evjen晶胞最外层离子与参考离子异号时计算的马德隆常数迥然不同.可见利用Evjen晶胞的方法计算马德隆常数,不便于推广使用,尤其对于复杂的离子晶体,以参考离子为中心构造一个比一个更大的Evjen晶胞并确定相应立方体边上、面上、棱上的正负离子数比较困难.因此,此法只适用于一些简单立方晶系的离子晶体马德隆数的计算,而不能计算复杂离子晶体的马德隆常数.
计算马德隆常数的目的是计算晶格静电能,因此不妨从晶格静电能出发计算马德隆常数.一摩尔离子晶体的晶格能UT是指晶体内各离子间静电相互吸引能
UC和玻尔排斥能UB,即卩L T"U C+=U C(1-1/n)=1/2N A U C分子(1-1/n)=1/2N A*U C晶胞/m 式中NA uC分子、uC晶胞分别为阿伏伽德罗常数和分子、晶胞的静电能,m、n分别为晶胞内分子数和玻恩指数;1/2是计算相互作用能时为避免重复计算而引入的设晶体一个晶胞涉及有k个正离子和w个负离子,则
u
c晶胞P j u Cj q j u Cj
1 j k k j (k w)
1
Uc B晶胞u c晶胞(1-—)
n
式中j=1,2,…,k对应于晶体内一个参考晶胞所涉及的k个正离子的编号;j=k+1,k+2,…,k+w对应于晶体内一个参考晶胞所涉及的w个负离子的编号. 对于立方晶胞,若离子处在界面上或棱边上或顶角上,则对应的p j和q分别取1/2 或1/4.其他晶胞与此类同;u+cj,u-cj, a+, a -j分别为参考晶胞中第j个离子静电能和相应的马德隆常数,其表达式为
Z-Z f +
叱=
4托切心)
jr
==—心p ( p fnjZ-^/Z-匸l( 1
J
Z+Z—/ _
wr.=
_4叭心
——一1/Z -/ Z +
竹柑孑 f jj辰” (盘f/j
式中i=1,2,3,…,k对应于计算时所及晶胞内正离子的编号;而i=k+1,k+2,…,k+w对应于计算时所及晶胞内负离子的编号m和l i的取值方法与P j和q相同;Z+和乙分别为正负离子价电子数;n i,门2小3分别为晶胞沿x,y,z方向堆积数.计算时须排除离子自身相互作用,r ij为计算时所及晶胞内第i离子与参考晶胞内第j离子的间距•若用马德隆常数a来表示晶体的结合能,则有
f J MZ + Z 一呵]1
+ B —
4応1 n j
由此可得晶体的马德隆常数a⑴为
注意a (1)只能用于二元化合物晶体,当晶体为二元以上化合物晶体时,须引入诸离子价电荷数Zj(j=1,2,…,(k+w))之间的最大公因子Z.由此,参考晶胞中第j 个离子静电能和相应的马德隆常数为U+Cj,U-Cj, a +j, a -j分别为
(Z"
4nE
n 7?
n
Z jr ( Z
-i~ )/( Z -I-X
乙(二nt,
"十it? iWfWJt
(Z E F-
4JT£O/?O J
2/i I
空口+3 Fy J
(/ =女+ 1 M + 2* …•佐+ R)
相应的晶格能和马德隆常数a (2分别为
% = /?()I S c s
其马德隆常数a⑵比a⑴具有普遍性.也可以不引入最大公因子Z来定义马德隆常数a(3),即
上述三种马德隆常数的关系为
利用这种方法可以计算出各种晶体结构的马德隆常数,如CsCI: a - 1.76266466,与文献值a =1.76268比较接近,可见这种方法的精确性.
综上所述,对于简单的离子晶体,可米用定义法直接对离子晶体计算马德隆常数;对于简单立方晶系离子晶体马德隆常数的计算,宜采用Evjen晶胞的方法. 对于复杂离子晶体,应使用计算晶格静电能法计算其马德隆常数.
参考文献:
【1】黄昆,韩汝琦.固体物理学[M].北京:高等教育出版社,2001.
【2】令狐荣锋,陈明能.马德隆常数的数值确定.贵州师范大学理学院, 1004—5570(2005)01 -0095 -02
【3】王矜奉,等.采用双Evjen晶胞计算离子晶体的马德隆常数[J].四川师范大学学报,2001,24(5).
【4】张维佳,等.复杂离子晶体马德隆常数研究[J].物理学报,2005,54(2).。