黄昆版固体物理学课后答案解析答案

黄昆 原着 韩汝琦改编 （陈志远解答，仅供参考）

第一章 晶体结构

1.1、

解：实验表明，很多元素的原子或离子都具有或接近于球形对称结构。因此，可以把这些原子或离子构成的晶体看作是很多刚性球紧密堆积而成。这样，一个单原子的晶体原胞就可以看作是相同的小球按点阵排列堆积起来的。它的空间利用率就是这个晶体原胞所包含的点的数目n和小球体积V所得到的小球总体积nV与晶体原胞体积Vc之比，即：晶体原胞的空间利用率， x43nV Vc（1）对于简立方结构：（见教材P2图1-1） a=2r， V=r3，Vc=a3，n=1

4343rr33∴x330.52 6a8r（2）对于体心立方：晶胞的体对角线BG=3a4ran=2, Vc=a3

442r32r3333∴x0.68 8a3433(r)343x 3（3）对于面心立方：晶胞面对角线BC=2a4r,a22r n=4，Vc=a3

（4）对于六角密排：a=2r晶胞面积：S=6SABO6晶胞的体积：V=SC16123328aa32a3242r3 23aasin60332=a 22n=121223=6个

（5）对于金刚石结构，晶胞的体对角线BG=3a42ra1.2、试证：六方密排堆积结构中()1/21.633

证明：在六角密堆积结构中，第一层硬球A、B、O的中心联线形成一个边长a=2r的正三角

形，第二层硬球N位于球ABO所围间隙的正上方并与这三个球相切，于是： NA=NB=NO=a=2R.

即图中NABO构成一个正四面体。…

ca838r3 n=8, Vc=a3

1.3、证明：面心立方的倒格子是体心立方；体心立方的倒格子是面心立方。

aa12(jk)a证明：（1）面心立方的正格子基矢（固体物理学原胞基矢）：a2(ik)

2aa32(ij)由倒格子基矢的定义：b12(a2a3) 0,a1(a2a3)a,2a,2a,20,a,2ai,2aaa3，a2a3,224a,02j,0,a,2kaa2(ijk) 2402(ijk)a同理可得：即面心立方的倒格子基矢与体心立方的正格基矢相同。

2b3(ijk)ab2所以，面心立方的倒格子是体心立方。

aa12(ijk)a（2）体心立方的正格子基矢（固体物理学原胞基矢）：a2(ijk)

2aa32(ijk)由倒格子基矢的定义：b12(a2a3) aaa,,i,j,k222aaaa2aaaa3(jk) ，a2a3,,a1(a2a3),,22222222aaaaaa,,,,2222222(ik)a同理可得：即体心立方的倒格子基矢与面心立方的正格基矢相同。

2b3(ij)ab2所以，体心立方的倒格子是面心立方。

1.5、证明倒格子矢量Gh1b1h2b2h3b3垂直于密勒指数为(h1h2h3)的晶面系。

证明：因为CA

a1a3aa,CB23，Gh1b1h2b2h3b3 h1h3h2h3利用aibj2ij，容易证明

Gh1h2h3CA0Gh1h2h3CB0

所以，倒格子矢量Gh1b1h2b2h3b3垂直于密勒指数为(h1h2h3)的晶面系。

1.6、对于简单立方晶格，证明密勒指数为(h,k,l)的晶面系，面间距d满足：

d2a2(h2k2l2)，其中a为立方边长；并说明面指数简单的晶面，其面密度较大，容

易解理。 解：简单立方晶格：a1a2a3，a1ai,a2aj,a3ak 由倒格子基矢的定义：b12倒格子基矢：b1a2a3a3a1aa，b22，b3212

a1a2a3a1a2a3a1a2a3222i,b2j,b3k aaa222倒格子矢量：Ghb1kb2lb3，Ghikjlk

aaa晶面族(hkl)的面间距：d2G1

hkl()2()2()2aaa面指数越简单的晶面，其晶面的间距越大，晶面上格点的密度越大，单位表面的能量越小，

这样的晶面越容易解理。

1.9、画出立方晶格（111）面、（100）面、（110）面，并指出（111）面与（100）面、（111）面与（110）面的交线的晶向。

解：(111)

1、(111)面与(100)面的交线的AB，AB平移，A与O点重合，B点位矢：RBajak， (111)面与(100)面的交线的晶向ABajak，晶向指数[011]。

2、(111)面与(110)面的交线的AB，将AB平移，A与原点O重合，B点位矢：RBaiaj，(111)面与(110)面的交线的晶向ABaiaj，晶向指数[110]。

第二章 固体结合

2.1、两种一价离子组成的一维晶格的马德隆常数（2ln2）和库仑相互作用能，设离子的总数为2N。

＜解＞ 设想一个由正负两种离子相间排列的无限长的离子键，取任一负离子作参考离子（这样马德隆常数中的正负号可以这样取，即遇正离子取正号，遇负离子取负号），用r表示相邻离子间的距离，于是有

rj(1)11112[...] rijr2r3r4r前边的因子2是因为存在着两个相等距离ri的离子，一个在参考离子左面，一个在其右面，故对一边求和后要乘2，马德隆常数为 111

2[1...]2342xx3x4... n(1x)xx34111当X=1时，有1...n2 2n22342.3、若一晶体的相互作用能可以表示为 试求：（1）平衡间距r0；

（2）结合能W（单个原子的）；

（3）体弹性模量；

（4）若取m2,n10,r03A,W4eV，计算及的值。 解：（1）求平衡间距r0

由

du(r)dr0，有：

rr0结合能：设想把分散的原子（离子或分子）结合成为晶体，将有一定的能量释放出

来，这个能量称为结合能（用w表示） （2）求结合能w（单个原子的）

题中标明单个原子是为了使问题简化，说明组成晶体的基本单元是单个原子，而非

原子团、离子基团，或其它复杂的基元。

显然结合能就是平衡时，晶体的势能，即Umin

即：WU(r0)（3）体弹性模量

r0mr0n （可代入r0值，也可不代入）

r02由体弹性模量公式：k9V02Ur2 r0（4）m = 2，n = 10，r03A， w = 4eV，求α、β

10 r021858 ① 44eV ② 5r021WU(r0)将r03A，1eV1.6021019J代入①②

（1）平衡间距r0的计算

N晶体内能U(r)(mn)

2rrdU平衡条件

drrr01nnmn)m 0，m1n10，r0(r0r0m（2）单个原子的结合能

1nn1)m Wu(r0)，u(r0)(mn)，r0(m2rrrr02U（3）体弹性模量K(2)V0V0

V晶体的体积VNAr3，A为常数，N为原胞数目 晶体内能U(r)由平衡条件

N(mn) 2rrUVVV0mnNmn1n ，得(m1n1)0m2r0r02r0r03NAr0体弹性模量KU0mn 9V0（4）若取m2,n10,r03A,W4eV

1mnn1mnr0()m，W(1)()nm

m2nmW10r0，r02[102W]

r021.210-95eVm10，9.01019eVm2

2.6、bcc和fcc Ne的结合能，用林纳德—琼斯(Lennard—Jones)势计算Ne在bcc和fcc

结构中的结合能之比值．

1＜解＞u(r)4()12()6,u(r)N(4)An()12Al()6

r2rrr2.7、对于H2，从气体的测量得到Lennard—Jones参数为50106J,2.96A.计算fcc结构的H2的结合能[以KJ/mol单位)，每个氢分子可当做球形来处理．结合能的实验值为0.751kJ／mo1，试与计算值比较．

＜解＞ 以H2为基团，组成fcc结构的晶体，如略去动能，分子间按Lennard—Jones势相互作用，则晶体的总相互作用能为：

jP614.45392;Pij1212.13188,iji501016erg,2.96A,N6.0221023/mol.将R0代入U得到平衡时的晶体总能量为U26。02210/mol50102816126因此，2.962.96erg12.1314.452.55KJ/mol.3.163.16计算得到的H2晶体的结合能为2．55KJ／mol，远大于实验观察值0.75lKJ／mo1．对于H2的晶体，量子修正是很重要的，我们计算中没有考虑零点能的量子修正，这正是造成理论和实验值之间巨大差别的原因．

第三章 固格振动与晶体的热学性质

3.1、已知一维单原子链，其中第j个格波，在第n个格点引起的位移为，

njajsin(jt_naqjj)，j为任意个相位因子，并已知在较高温度下每个格波的平均能量为，具体计算每个原子的平方平均位移。

＜解＞任意一个原子的位移是所有格波引起的位移的叠加，即 nnjajsin(jtnaqjj) （1）

jj由于njnj数目非常大为数量级，而且取正或取负几率相等，因此上式得第2项与第一项

22相比是一小量，可以忽略不计。所以n njj由于nj是时间t的周期性函数，其长时间平均等于一个周期内的时间平均值为

2j1T0T00a2jsin(jtnaqjj)dt12aj （2） 2已知较高温度下的每个格波的能量为KT，nj的动能时间平均值为

其中L是原子链的长度，使质量密度，T0为周期。 所以Tnj112w2LaKT （3） jj422因此将此式代入（2）式有njKT PL2j22所以每个原子的平均位移为 nnjjjKTKT1 22PLjPLjj3.2、讨论N个原胞的一维双原子链（相邻原子间距为a），其2N个格波解，当M= m时与一维单原子链的结果一一对应。

解：质量为M的原子位于2n-1， 2n+1， 2n+3 ……；质量为m的原子位于2n， 2n+2， 2n+4 ……。 牛顿运动方程

m2n(22n2n12n1)M2n1(22n12n22n)

N个原胞，有2N个独立的方程 设方程的解

2nAei[t(2na)q]2n1Bei[t(2n1)aq]，代回方程中得到

2m22cosaqA、B有非零解，0，则

2cosaq2M21(mM)4mM2{1[1sinaq]2}2mM(mM)22两种不同的格波的色散关系

(mM)4mM{1[1sin2aq]}2mM(mM)12

一个q对应有两支格波：一支声学波和一支光学波.总的格波数目为2N.

当Mm时

4aqcosm24aqsinm2，

两种色散关系如图所示： 长波极限情况下q0，sin(qaqa)， 22(2m)q与一维单原子晶格格波的色散关系一致.

3.3、考虑一双子链的晶格振动，链上最近邻原子间的力常数交错地为和10，两种原子质量相等，且最近邻原子间距为a2。试求在q0,qa处的(q)，并粗略画出色散关

系曲线。此问题模拟如H2这样的双原子分子晶体。

答：（1）

浅色标记的原子位于2n-1， 2n+1， 2n+3 ……；深色标记原子位于2n， 2n+2， 2n+4 ……。 第2n个原子和第2n＋1个原子的运动方程： 体系N个原胞，有2N个独立的方程 方程的解：

2nAe1i[t(2n)aq]21i[t(2n1)aq]222/m，将解代入上述方程得： ，令121/m,22n1BeA、B有非零的解，系数行列式满足：

c10c222因为1、210，令0得到 12,2100mm2两种色散关系：20(1120cosqa101)

当q0时，(11121)，

2202200

当qa时，(1181)，

22020020

（2）色散关系图：

3.7、设三维晶格的光学振动在q=0附近的长波极限有(q)0Aq2 求证：f()V11/2,0;f()0,0. 042A3/21212＜解＞0时，0Aq20f()0,00Aq2qA0 依据q(q)2Aq,f()V23ds，并带入上边结果有

q(q)3.8、有N个相同原子组成的面积为S的二维晶格，在德拜近似下计算比热，并论述在低温极限比热正比与T2。

证明：在k到kdk间的独立振动模式对应于平面中半径n到ndn间圆环的面积2ndn，

L253s且2ndnkdkkdk即d则 2222v3skBTE012v2233sE2v2m2de/kBT0DD3d3skTkTkTBBB/kBTe12v222xDDx2dx， ex13.9、写出量子谐振子系统的自由能，证明在经典极限下，自由能为FU0kBTqq nkTB1q证明:量子谐振子的自由能为FUkBTq2kBTqkBTn1e 经典极限意味着（温度较高）kBT应用ex1xx2... 所以eg

qkBT1q...

kBTkBTq21因此FUqkBTq2q其中U0UqqU0kBTn11kTBq nkTB1q 23.10、设晶体中每个振子的零点振动能为

1，使用德拜模型求晶体的零点振动能。 2证明:根据量子力学零点能是谐振子所固有的，与温度无关，故T=0K时振动能E0就是各振动模零点能之和。E0E0m0E0gd将E03V1和g232代入积分有

2vs23V994N，由于k得ENkBD mmmBD02316vs88一股晶体德拜温度为~102K，可见零点振动能是相当大的，其量值可与温升数百度所需热能相比拟．

3.11、一维复式格子m51.671024g,00A,min（1），光学波max，声学波max。

M4,1.5101N/m(即1.51104dyn/cm),求m（2）相应声子能量是多少电子伏。

（3）在300k时的平均声子数。

0（4）与max相对应的电磁波波长在什么波段。

＜解＞（1），Axam221.5104dyn/cm3113.0010s, 24M451.6710Amax6.5810165.991013s11.97102eVo6.5810166.701013s14.41102eV （2）maxomin6.5810163.001013s13.95102eVAn（3）max1eAmax/kBT1O0.873,nmax1eOmax/kBT10.221

（4）2c28.1m

第四章 能带理论

4.1、根据ka状态简并微扰结果，求出与E及E相应的波函数及?，并说明它

2们的特性．说明它们都代表驻波，并比较两个电子云分布说明能隙的来源(假设Vn=Vn*)。

＜解＞令ka，ka，简并微扰波函数为Ak0(x)Bk0(x)

0 VnAEkEB0 取EE

带入上式，其中EE0(k)Vn V(x)<0,Vn0,从上式得到B= -A,于是

AA(x)(x)L0k0knxixina2Anasinx ee=aL取EE，EE0(k)Vn VnAVnB,得到AB

AA(x)(x)L0k0knxixina2Anacosx ee=aL 由教材可知，及均为驻波． 在驻波状态下，电子的平均速度(k)为零．产生

n22a时，电子波的波长，恰好满足布拉格发射条件，这时

akn电子波发生全反射，并与反射波形成驻波由于两驻波的电子分布不同，所以对应不同代入能量。

4.2、写出一维近自由电子近似，第n个能带(n=1，2，3)中，简约波数k的0级波函

2a数。

驻波因为电子波矢k221mxximxi(m)x1ikx1ikxi2a1i21＜解＞(x)eeeeaeaea4

LLLL*kx1i2a第一能带：m0,m0,(x)e

2aL*k23ixixi221*第二能带：bb则bb,m,即m1,（ea=e2a）k(x)e2a

aaL2xixx221i21i5*aa2a第三能带：cc,m,即m1,k(x)eee

aaLL4.3、电子在周期场中的势能．

V(x) 0 ， 当(n-1)a+bxnab

其中d＝4b，是常数．试画出此势能曲线，求其平均值及此晶体的第一个和第二个禁带

度．

＜解＞(I)题设势能曲线如下图所示． (2)势能的平均值：由图可见，V(x)是个以a为周期的周期函数，所以

题设a4b，故积分上限应为ab3b，但由于在b,3b区间内V(x)0，故只需在

b,b区间内积分．这时，n0，于是

1bm2b2m22VV(x)dx(bx)dxbba2a2a2bxbb1x33bb12。 mb6（3），势能在[-2b,2b]区间是个偶函数，可以展开成傅立叶级数

u2利用积分公式u2cosmudu2musinmu2cosmum3sinmu得 mEg116m23b2第二个禁带宽度Eg22V2,以m2代入上式,代入上式

b22m2Eg2b0(bx)cosxbdx再次利用积分公式有Eg22m22b2

4.4、

解：我们求解面心立方，同学们做体心立方。

（1）如只计及最近邻的相互作用，按照紧束缚近似的结果，晶体中S态电子的能量可表示成：

在面心立方中，有12个最近邻，若取Rm0，则这12个最近邻的坐标是：

aaaa①(1,1,0),(1,1,0),(1,1,0),(1,1,0) 2222aaaa②(0,1,1),(0,1,1),(0,1,1),(0,1,1) 2222aaaa③(1,0,1)(1,0,1),(1,0,1),(1,0,1) 2222由于S态波函数是球对称的，在各个方向重叠积分相同，因此J(RS)有相同的值，简单表示为J1=J(RS)。又由于s态波函数为偶宇称，即s(r)s(r)

∴在近邻重叠积分J(Rs)i*(Rs)U()V(Rs)i()d中，波函数的贡献为正 ∴J1＞0。

于是，把近邻格矢RS代入Es(RS)表达式得到：

aaa(kxky)i(kxky)i(kxky)i(kxky)ia222=SJ0J1e eee2ai(kykz)2ai(kykz)2ai(kykz)2ai(kykz)2eeee+

eai(kxkz)2eai(kxkz)2eai(kxkz)2eai(kxkz)2 aaaa=SJ02J1cos(kxky)cos(kxky)cos(kykz)cos(kykz)

2222aaaaaa=sJ04J1coskxcoskycoskycoskzcoskzcoskx

222222（2）对于体心立方：有8个最近邻，这8个最近邻的坐标是：

4.7、有一一维单原子链，间距为a，总长度为Na。求（1）用紧束缚近似求出原子s态能级对应的能带E(k)函数。（2）求出其能态密度函数的表达式。（3）如果每个原子s态只有

00一个电子，求等于T=0K的费米能级EF及EF处的能态密度。

＜解＞(1),E(k)sJ0J1(eikaeika)sJ02J1coskaE02J1coska (2) ,N(E)20kFLdk2Na1N 22dE2J1asinkaJ1sinka(3), N002NakFNa002(k)2dk22kFkF

22a4.8、证明一个自由简单晶格在第一布里渊区顶角上的一个自由电子动能比该区一边中点大

2倍．(b)对于一个简单立力晶格在第一布里渊区顶角上的一个自由电子动能比该区面心上大多少？(c)(b)的结果对于二价金属的电导率可能会产生什么影响7

＜解＞（a）二维简单正方晶格的晶格常数为a，倒格子晶格基矢A2ˆ2ˆi,Bj aa第一布里渊区如图所示

a 0

区边中点的波矢为KA2ˆˆˆi,角顶B点的波矢为KBij.aaa2x2KyKz2,自由电子能量2K2mA点能量AKx2,

2m2ma2ma222222B点能量BKxKy2m2,所以B/A2 2maa2ma2222222b)简单立方晶格的晶格常数为a，倒格子基矢为Aaˆ2ˆ2i,Bj,Caaˆk, 第一布里渊区如图7—2所示．

所以B/A3

(c)如果二价金属具有简单立方品格结构，布里渊区如图7—2所示．根据自由电子理论，自由电子的能量为4切于4点的内切球的体积

3K2m322x2KyKz2，FerM面应为球面．由(b)可知，内

，于是在K空间中，内切球内能容纳的电子数为a43V2N1.047N 其中VNa3 3a233二价金属每个原子可以提供2个自由电子，内切球内只能装下每原子1.047个电子，余下的0.953个电子可填入其它状态中．如果布里渊区边界上存在大的能量间隙，则余下的电子只能填满第一区内余下的所有状态(包括B点)．这样，晶体将只有绝缘体性质．然而由(b)可知，B点的能员比A点高很多，从能量上看，这种电子排列是不利的．事实上，对于二价金属，布里渊区边界上的能隙很小，对于三维晶体，可出现

一区、二区能带重迭．这样，处于第一区角顶附近的高能态的电子可以“流向”第二区中的能量较低的状态，并形成横跨一、二区的球形Ferm面．因此，一区中有空态存在，而二区中有电子存在，从而具有导电功能．实际上，多数的二价金届具有六角密堆和面心立方结构，能带出现重达，所以可以导电． 4.10、

解：设晶体中有N个Cu原子，向其中掺入x个锌原子。则晶体中电子的总数为： (N-x)+2x=N+x

N由于Cu是面心立方，每一个原胞中含4个电子。因此：晶体中包含的原胞数为：

4其倒格子为体心立方，倒格子的边长为：

434，对角线的长度为：

aa1433于是：布里渊区边界到原点的距离为：

4aa即：当Fermi球与第一布里渊区边界相切时，kF又由：23 a43kFNx 323V于是有：

Nx3Nx33 N3aN4a4即：当锌原子与铜原子之比为0.56时，Fermi球与第一布里渊区边界相接触。

2x2ycos4.12、正方晶格．设有二维正方晶格，晶体势为Ux,y4Ucos. aa用基本方程，近似求出布里渊区角,处的能隙．

aaˆ,bˆˆ,ˆj表示位置矢量的单位矢量，以b ＜解＞以i12表示倒易矢量的单位矢量，则有，

ˆGbˆ2gˆˆ,g,g为整数。ˆyiˆ,GG rxib11221b1g2b212a2x2ycos晶体势能Ux,y4Ucos. aai2xi2xi2yi2yiG11UrUeeeeUG11eG11其中UG11U，而其他势能傅氏系数UG10UG20...0。这样基本方程

kCKUGG(KG)0变为

GKCKUG11CKG11UG11CKG11UG11CKG11UG11CKG110求布

111里渊区角顶,，即kG(,)G11处的能隙，可利用双项平面波近似

222aaC(K)eiKrC(KG)ei(KG)r来处理。

11当KG11,KG11时依次有

2211KG11G11,KG11G11而其他的KG11，

22KG11G11，所以在双项平面波近似下上式中只有

11CG11,CKG11CG11;22 11CG11,CKG11CG11;22 u 1G112 =0，因为

第五章 晶体中电子在电场和磁场中的运动

5.1、设有一维晶体的电子能带可写成 E(k)常数，m是电子的质量。 试求（1）能带宽度；

（2）电子在波矢k状态的速度； （3）带顶和带底的电子有效质量。 解：（1） E(k)71(coskacos2ka) 2ma88271(coskacos2ka)， 其中a为晶格2ma882127 =－coska+(2cos2ka－1)] 288ma2 ＝

4ma2(coska－2)2－1

当ka＝(2n+1)时,n=0,1,2… 当ka=2n时, Emin(k)0 能带宽度＝EmaxEmin （2）22 ma2

1dE(k)1(sinkasin2ka) dkma421 (3) m2Em(coskacos2ka)1

2k2* 当k0时,带底，m*2m 当k5.5、

2时,带顶，m*m a31解：（1）电子的运动速度：vkE(k)

dvd11dE (kE)kdtdtdt由于单位时间内能量的增加=力在单位时间内作的功

dEds1即：FFvkEF

dtdt∴ 加速度：∴

dv11EEE2k[kEF]2k[F1F2F3] dtk1k2k3写成分量的形式：

112E其中： （i,j=1,2,3） 2mijkikj22k32k12k2由题知：E 2m12m22m321111112E1,容易得出： 22 同理：mmmmm11k1m122233311111112E0 20 同理：

m13m21m23m31m31m12k1k2dv1m1dtF1dv故运动方程为：m22F2

dtdv3m3dtF3（2）当存在磁场B作用时，电子将受到洛仑兹力作用

当B相对于椭球主轴的方向余弦为,,时，电子的运动方程可写成： ∴ 电子的运动方程可写成：

dv12vB3)v23v3m1dtF1e(vBdv∴ m22F2e(vB3VB1)v31v 1dtdv3v21VB1)v12m3dtF3e(vB其中：1eB,2eB,3eB

231由于电子在磁场B作用下作周期性运动，故可设试探解：

v1v10eitv2v20eit 代入上述方程组可得： v3v30eitim1v12v33v2im1v13v22v30im2v23v11v3 即3v1im2v21v30 imvvvvvimv033122123321v1,v2,v3有非零解的条件是

222222m11m22m3322m1m2m3即：=eB

m2m2m3m1m2m32mm2m3e2B2∴=*2 m*22212 222m1m2m3m2

mmmeB即： 其中：m*221222322 证毕

m*m1m2m31/2第六章 金属电子论

第七章 半导体电子论

7.1、InSb电子有效质量me0.015m，介电常数18，晶格常数a6.49A。试计算；(1)施主的电离能；(2)基态轨道的半径；(3)施主均匀分布，相邻杂质原于的轨道之间将产生交叠时掺有的施主浓度应该高于多少？

＜解＞（1）由于施主电离能ED是氢原子电离能Ei的

m*倍， m02420m0170.5228a(A)6.3110(A)6.3110(m) （2），a002m*em*0.014（3），如果施主的电子与类氢基态轨道发生重叠，则均匀分布于InSb中施主杂质浓度ND就

一定满足(2a)3ND1,ND(131)4.981020(m2) 832a(26.3110)第十二章 晶体中的缺陷和扩散

例1．假设把一个钠原子从钠晶体内部移到边界上所需要的能量为1ev(1)，试计算室温（300k）时，sckottky空位的浓度？（已知：Na=0.97克/原米3，原子量为23）

解：（1）设N为单位体积内的Na原子数，则在温度T时，schottky定位的浓度n可写成：

nNeuKBT

由题知：u=1ev=1.602 ×10-19J

Na=0.97克/厘米3每cm3含Na的mol数为：

0.97 23每cmNa中所含的原子数为： N=

3

0.976.021023 23u0.9723236.0210e1.3810300 231.61019于是：nNeKBT例2. 如果u代表形成一个Frenkel缺陷所需的能量，证明在温度T时，达到热平衡的晶

u体中，Frenkel缺陷的数目为：nNN'e2KBT.

解：达到热平衡时，在N个原子的晶体中形成n个空位的可能方式数为： 这n个原子排列在N＇个间隙位置上的可能方式数为：

这样，从N个原子中取出n个原子并把它们排n＇个间隙位置上的总方式数为：

=12=

N!N＇!

(Nn)!n!(N＇n)!n!由此引起的熵的增量为：

利用斯特令公式:luN!=NlnN-N

得skBNlnN(Nn)ln(Nn)nlnnkBN'lnN'(N'n)ln(N'n)nlnn 系统的自

由能改变：FnuTs

U为形成一个Frenkel缺陷所需的能量。

F由热平衡条件：0有：

nTun2eKBT 即：

(N'n)(Nn)证毕

例3．铁一碳合金是面心立方结构，晶体常数为3.61A，设碳占的比重为1.7%试计算当碳是填隙式或替代式渗入时，合金的密度各等于多少？

解：C的mol质量：12g，C：含量为：1.7%

Fe的mol质量：55.85g,Fe：应为：98.3%

因此：在100克合余中，如果C是替代式溶入，则

100g合金的总mol数=0.1415+1.76=1.9015(mol)

在1mol的合金中：C的mol数=Fe的 mol数=

0.14150.0745mol

1.90151.760.9255mol.

1.9015于是替代式的合金密度：

=7.43克/厘米3。

如果C以填隙方式溶入铁的晶格中

7.4555.851292.554236.0210=7.88克/厘米3 则833.6110显然，C以不同的方式存在于Fe中，其密度的差别较大，因此，可以用实验的方法确定碳原子在合余中的存在方式。

例4．设t=0时，在晶体表面每单位面积上集中着扩散粒子N，试求径过t时间后，扩散粒子的均方位移？x22Dt

解：本题属于有源扩散问题。t时刻，x处扩散粒子的浓度是： 扩散粒子的均方位移可按下式求出：

例5．已知c在Fe中的扩散系数D与温度关系的实验数据为：当温度为200℃时，扩散

系数D200℃=10-11cm2/秒。温度为760℃时，D760℃=10-6cm2/秒，试求扩散过程的激活能Q（千焦耳/摩尔）。 气体常数R= 8.31J/mol·K.[华工99年研究生入学考试三（15分）] 解：由D=Doe-Q/RT有：D0DeQ/RT =83.38kJ/mol. 作业：

（1）对于铜，形成一个肖脱基空位的能量为1.2ev，在接近熔点时（1300K），估算空位的浓度。形成一个填隙原子所需的能量约为4ev，估算接近熔点时填隙原子的浓度。试比较这两种浓度的数量级差多少？

（2）在离子晶体中，由于电中性的要求，Schottky缺陷都是成对地产生的。令n代表正、负离子空位的对数，E是形成一对缺陷所需要的能量，N为整个离子晶体中正、负离子对的数目；证明：

n≈Ne

E2KBT.

(3)已知C在铁中的扩散系数D和温度T的关系一组实验数据如下 T/℃ -10 27 200 760 D/（cm2/s） 5×10-19 10-16 10-11 10-6 求1mol原子的激活能Q=N0和频率因子D0.