固体物理第二章答案
- 格式:doc
- 大小:1.35 MB
- 文档页数:16
《固体物理学》习题解答黄昆 原着 韩汝琦改编 (陈志远解答,仅供参考)第一章 晶体结构1.1、解:实验表明,很多元素的原子或离子都具有或接近于球形对称结构。
因此,可以把这些原子或离子构成的晶体看作是很多刚性球紧密堆积而成。
这样,一个单原子的晶体原胞就可以看作是相同的小球按点阵排列堆积起来的。
它的空间利用率就是这个晶体原胞所包含的点的数目n 和小球体积V 所得到的小球总体积nV 与晶体原胞体积Vc 之比,即:晶体原胞的空间利用率, VcnVx = (1)对于简立方结构:(见教材P2图1-1)a=2r , V=3r 34π,Vc=a 3,n=1∴52.06r 8r34a r 34x 3333=π=π=π= (2)对于体心立方:晶胞的体对角线BG=x 334a r 4a 3=⇒= n=2, Vc=a 3∴68.083)r 334(r 342a r 342x 3333≈π=π⨯=π⨯= (3)对于面心立方:晶胞面对角线BC=r 22a ,r 4a 2=⇒= n=4,Vc=a 3(4)对于六角密排:a=2r 晶胞面积:S=6260sin a a 6S ABO ⨯⨯=⨯∆=2a 233 晶胞的体积:V=332r 224a 23a 38a 233C S ==⨯=⨯ n=1232126112+⨯+⨯=6个 (5)对于金刚石结构,晶胞的体对角线BG=3r 8a r 24a 3=⇒⨯= n=8, Vc=a 31.2、试证:六方密排堆积结构中633.1)38(a c 2/1≈= 证明:在六角密堆积结构中,第一层硬球A 、B 、O 的中心联线形成一个边长a=2r 的正三角形,第二层硬球N 位于球ABO 所围间隙的正上方并与这三个球相切,于是: NA=NB=NO=a=2R.即图中NABO 构成一个正四面体。
…1.3、证明:面心立方的倒格子是体心立方;体心立方的倒格子是面心立方。
证明:(1)面心立方的正格子基矢(固体物理学原胞基矢):123()2()2()2a a j k a a i k a a i j ⎧=+⎪⎪⎪=+⎨⎪⎪=+⎪⎩由倒格子基矢的定义:1232()b a a π=⨯Ω31230,,22(),0,224,,022a aa a a a a a a a Ω=⋅⨯==,223,,,0,()224,,022i j ka a a a a i j k a a ⨯==-++ 同理可得:232()2()b i j k ab i j k aππ=-+=+-即面心立方的倒格子基矢与体心立方的正格基矢相同。
第一章 金属自由电子气体模型习题及答案1. 你是如何理解绝对零度时和常温下电子的平均动能十分相近这一点的?[解答] 自由电子论只考虑电子的动能。
在绝对零度时,金属中的自由(价)电子,分布在费米能级及其以下的能级上,即分布在一个费米球内。
在常温下,费米球内部离费米面远的状态全被电子占据,这些电子从格波获取的能量不足以使其跃迁到费米面附近或以外的空状态上,能够发生能态跃迁的仅是费米面附近的少数电子,而绝大多数电子的能态不会改变。
也就是说,常温下电子的平均动能与绝对零度时的平均动能十分相近。
2. 晶体膨胀时,费米能级如何变化?[解答] 费米能级3/222)3(2πn mE o F= , 其中n 单位体积内的价电子数目。
晶体膨胀时,体积变大,电子数目不变,n 变小,费密能级降低。
3. 为什么温度升高,费米能反而降低?[解答] 当K T 0≠时,有一半量子态被电子所占据的能级即是费米能级。
除了晶体膨胀引起费米能级降低外,温度升高,费米面附近的电子从格波获取的能量就越大,跃迁到费米面以外的电子就越多,原来有一半量子态被电子所占据的能级上的电子就少于一半,有一半量子态被电子所占据的能级必定降低,也就是说,温度生高,费米能反而降低。
4. 为什么价电子的浓度越大,价电子的平均动能就越大?[解答] 由于绝对零度时和常温下电子的平均动能十分相近,我们讨论绝对零度时电子的平均动能与电子的浓度的关系。
价电子的浓度越大,价电子的平均动能就越大,这是金属中的价电子遵从费米—狄拉克统计分布的必然结果。
在绝对零度时,电子不可能都处于最低能级上,而是在费米球中均匀分布。
由式3/120)3(πn k F =可知,价电子的浓度越大费米球的半径就越大,高能量的电子就越多,价电子的平均动能就越大。
这一点从3/2220)3(2πn m E F=和3/222)3(10353πn mE E oF ==式看得更清楚。
电子的平均动能E 正比于费米能o F E ,而费米能又正比于电子浓度32l n。
固体物理习题参考答案(部分)第一章 晶体结构1.氯化钠:复式格子,基元为Na +,Cl -金刚石:复式格子,基元为两个不等价的碳原子 氯化钠与金刚石的原胞基矢与晶胞基矢如下:原胞基矢)ˆˆ()ˆˆ()ˆˆ(213212211j i a a i k a a k j a a +=+=+= , 晶胞基矢 ka a j a a ia a ˆˆˆ321===2. 解:31A A O ':h:k;l;m==-11:211:11:111:1:-2:1 所以(1 1 2 1) 同样可得1331B B A A :(1 1 2 0); 5522A B B A :(1 1 0 0);654321A A A A A A :(0 0 0 1)3.简立方: 2r=a ,Z=1,()63434r 2r a r 3333πππ===F体心立方:()πππ833r4r 342a r 3422a 3r 4a r 4a 33333=⨯=⨯=∴===F Z ,,则面心立方:()πππ622r 4r 34434442r 4a r 4a 233ar 33=⨯=⨯=∴===F Z ,,则 六角密集:2r=a, 60sin 2c a V C = a c 362=,πππ622336234260sin 34223232=⨯⨯⨯=⨯=⎪⎭⎫ ⎝⎛a a c a r F a金刚石:()πππ163r 38r 348a r 3488Z r 8a 33333=⨯=⨯===F ,, 4. 解:'28109)31arccos(312323)ˆˆˆ()ˆˆˆ(cos )ˆˆˆ()ˆˆˆ(021*******12211=-=-=++-⋅+-=⋅=++-=+-=θθa a k j i a k j i a a a a a kj i a a kj i a a 5.解:对于(110)面:2a 2a a 2S =⋅=所包含的原子个数为2,所以面密度为22a2a22=对于(111)面:2a 2323a 22a 2S =⨯⨯= 所包含的原子个数为2,所以面密度为223a34a 232=8.证明:ABCD 是六角密堆积结构初基晶胞的菱形底面,AD=AB=a 。
黄昆版固体物理学课后答案《固体物理学》习题解答黄坤原著韩汝琦改编(陈志远答案,仅供参考)第一章晶体结构1.1、解:实验表明,很多元素的原子或离子都具有或接近于球形对称结构。
因此,可以把这些原子或离子构成的晶体看作是很多刚性球紧密堆积而成。
这样,一个单原子的晶体原胞就可以看作是相同的小球按点阵排列堆积起来的。
它的空间利用率就是这个晶体原胞所包含的点的数目n和小球体积v所得到的小球总体积nv与晶体原胞体积vc之比,即:晶体原胞的空间利用率,x?(1)对于简立方结构:(见教材p2图1-1)nvvc43?r、 vc=a3,n=134343?RR33∴十、0.526a38r3a=2r,v=(2)对于体心立方:晶胞的体对角线bg=3a?4r?a?n=2,vc=a343x32?∴十、434? r2??r33330.688a3433(R)3(3)对于面心立方:单元面对角线BC=2A?4r,?A.22rn=4,vc=a3444??r34??r3233x0.74336a(22r)(4)对于六角密排:a=2r晶胞面积:s=6?s?abo?6?晶胞的体积:v=s?c?a?asin60332a=223328a?a?32a3?242r323n=1212?11?2??3=6个6246??r323x0.7436242r(5)对于金刚石结构,晶胞的体对角线bg=3a?4?2r?a?8r3n=8,vc=a3一448??r38??r33?33x0.346a3833r33c81。
2.测试:在六角形紧密堆积结构中?()1/2? 一点六三三a3证明:在六角密堆积结构中,第一层硬球a、b、o的中心联线形成一个边长a=2r 的正三角形,第二层硬球n位于球abo所围间隙的正上方并与这三个球相切,于是:na=nb=no=a=2r.也就是说,图中的Nabo形成一个正四面体1.3、证明:面心立方的倒格子是体心立方;体心立方的倒格子是面心立方。
A.a1?2(j?k)A.证明了:(1)面心立方的法向晶格基向量(固体物理的原胞基向量):?a2?(i?k)2aa3?2(i?j)??2(a2?a3)由倒格子基矢的定义:b1??0,aa1?(a2?a3)?,2a,2a,20,a,2a?i,2aa3a???,a2?a3?,242a0,2?j,0,a,2?kaa2(?i?j?k)240?4a2???2?b1?2??3?(?i?j?k)?(?i?j? k)a4a?2.b2?(I?J?K)a类似地:2即面心立方的倒格子基矢与体心立方的正格基矢相同。
《固体物理学》习题解答黄昆 原着 韩汝琦改编 (陈志远解答,仅供参考)第一章 晶体结构、解:实验表明,很多元素的原子或离子都具有或接近于球形对称结构。
因此,可以把这些原子或离子构成的晶体看作是很多刚性球紧密堆积而成。
这样,一个单原子的晶体原胞就可以看作是相同的小球按点阵排列堆积起来的。
它的空间利用率就是这个晶体原胞所包含的点的数目n 和小球体积V 所得到的小球总体积nV 与晶体原胞体积Vc 之比,即:晶体原胞的空间利用率, VcnVx = (1)对于简立方结构:(见教材P2图1-1) a=2r , V=3r 34π,Vc=a 3,n=1∴52.06r8r34a r 34x 3333=π=π=π=(2)对于体心立方:晶胞的体对角线BG=x 334a r 4a 3=⇒= n=2, Vc=a 3∴68.083)r 334(r 342a r 342x 3333≈π=π⨯=π⨯= (3)对于面心立方:晶胞面对角线BC=r 22a ,r 4a 2=⇒= n=4,Vc=a 3(4)对于六角密排:a=2r 晶胞面积:S=6260sin a a 6S ABO ⨯⨯=⨯∆=2a 233 晶胞的体积:V=332r 224a 23a 38a 233C S ==⨯=⨯ n=1232126112+⨯+⨯=6个(5)对于金刚石结构,晶胞的体对角线BG=3r 8a r 24a 3=⇒⨯= n=8, Vc=a 3、试证:六方密排堆积结构中633.1)38(ac 2/1≈=证明:在六角密堆积结构中,第一层硬球A 、B 、O 的中心联线形成一个边长a=2r 的正三角形,第二层硬球N 位于球ABO 所围间隙的正上方并与这三个球相切,于是:NA=NB=NO=a=2R.即图中NABO 构成一个正四面体。
…、证明:面心立方的倒格子是体心立方;体心立方的倒格子是面心立方。
证明:(1)面心立方的正格子基矢(固体物理学原胞基矢):123()2()2()2a a j k a a i k a a i j ⎧=+⎪⎪⎪=+⎨⎪⎪=+⎪⎩r r r r r rr r r由倒格子基矢的定义:1232()b a a π=⨯Ωr r r31230,,22(),0,224,,022a a a a a a a a a a Ω=⋅⨯==r r rQ ,223,,,0,()224,,022i j ka a a a a i j k a a ⨯==-++r rr r r r r r同理可得:232()2()b i j k ab i j k aππ=-+=+-r rr r r r r r 即面心立方的倒格子基矢与体心立方的正格基矢相同。
第一章 晶体结构1.1、 如果将等体积球分别排成下列结构,设x 表示钢球所占体积与总体积之比,证明:结构 X简单立方52.06=π体心立方68.083≈π 面心立方74.062≈π 六角密排74.062≈π 金刚石34.063≈π解:实验表明,很多元素的原子或离子都具有或接近于球形对称结构。
因此,可以把这些原子或离子构成的晶体看作是很多刚性球紧密堆积而成。
这样,一个单原子的晶体原胞就可以看作是相同的小球按点阵排列堆积起来的。
它的空间利用率就是这个晶体原胞所包含的点的数目n 和小球体积V 所得到的小球总体积nV 与晶体原胞体积Vc 之比,即:晶体原胞的空间利用率, VcnVx = (1)对于简立方结构:(见教材P2图1-1) a=2r , V=3r 34π,Vc=a 3,n=1 ∴52.06834343333====πππrra r x (2)对于体心立方:晶胞的体对角线BG=x 334a r 4a 3=⇒= n=2, Vc=a 3∴68.083)334(3423423333≈=⨯=⨯=πππr r a r x (3)对于面心立方:晶胞面对角线BC=r 22a ,r 4a 2=⇒= n=4,Vc=a 374.062)22(3443443333≈=⨯=⨯=πππr r a r x (4)对于六角密排:a=2r 晶胞面积:S=6260sin a a 6S ABO ⨯⨯=⨯∆=2a 233晶胞的体积:V=332r 224a 23a 38a 233C S ==⨯=⨯ n=1232126112+⨯+⨯=6个 74.062)22(3443443333≈=⨯=⨯=πππr r a r x (5)对于金刚石结构,晶胞的体对角线BG=3r 8a r 24a 3=⇒⨯= n=8, Vc=a 334.06333834834833333≈=⨯=⨯=πππr r a r x 1.3、证明:面心立方的倒格子是体心立方;体心立方的倒格子是面心立方。
《固体物理学》习题解答黄昆 原著 韩汝琦改编 (陈志远解答,仅供参考)第一章 晶体结构1.1、解:实验表明,很多元素的原子或离子都具有或接近于球形对称结构。
因此,可以把这些原子或离子构成的晶体看作是很多刚性球紧密堆积而成。
这样,一个单原子的晶体原胞就可以看作是相同的小球按点阵排列堆积起来的。
它的空间利用率就是这个晶体原胞所包含的点的数目n 和小球体积V 所得到的小球总体积nV 与晶体原胞体积Vc 之比,即:晶体原胞的空间利用率, VcnVx = (1)对于简立方结构:(见教材P2图1-1)a=2r , V=3r 34π,Vc=a 3,n=1 ∴52.06r8r34a r 34x 3333=π=π=π= (2)对于体心立方:晶胞的体对角线BG=x 334a r 4a 3=⇒= n=2, Vc=a 3∴68.083)r 334(r 342a r 342x 3333≈π=π⨯=π⨯= (3)对于面心立方:晶胞面对角线BC=r 22a ,r 4a 2=⇒= n=4,Vc=a 374.062)r 22(r 344a r 344x 3333≈π=π⨯=π⨯= (4)对于六角密排:a=2r 晶胞面积:S=6260sin a a 6S ABO ⨯⨯=⨯∆=2a 233 晶胞的体积:V=332r 224a 23a 38a 233C S ==⨯=⨯ n=1232126112+⨯+⨯=6个 74.062r224r 346x 33≈π=π⨯= (5)对于金刚石结构,晶胞的体对角线BG=3r 8a r 24a 3=⇒⨯= n=8, Vc=a334.063r 338r 348a r 348x 33333≈π=π⨯=π⨯=1.2、试证:六方密排堆积结构中633.1)38(a c 2/1≈= 证明:在六角密堆积结构中,第一层硬球A 、B 、O 的中心联线形成一个边长a=2r 的正三角形,第二层硬球N 位于球ABO 所围间隙的正上方并与这三个球相切,于是: NA=NB=NO=a=2R.即图中NABO 构成一个正四面体。
第2章晶体的结合思考题1.是否有与库仑力无关的晶体结合类型?[解答]共价结合中, 电子虽然不能脱离电负性大的原子, 但靠近的两个电负性大的原子可以各出一个电子, 形成电子共享的形式, 即这一对电子的主要活动范围处于两个原子之间, 通过库仑力, 把两个原子连接起来. 离子晶体中, 正离子与负离子的吸引力就是库仑力. 金属结合中, 原子实依靠原子实与电子云间的库仑力紧紧地吸引着. 分子结合中, 是电偶极矩把原本分离的原子结合成了晶体. 电偶极矩的作用力实际就是库仑力. 氢键结合中, 氢先与电负性大的原子形成共价结合后, 氢核与负电中心不在重合, 迫使它通过库仑力再与另一个电负性大的原子结合. 可见, 所有晶体结合类型都与库仑力有关.2.如何理解库仑力是原子结合的动力?[解答]晶体结合中, 原子间的排斥力是短程力, 在原子吸引靠近的过程中, 把原本分离的原子拉近的动力只能是长程力, 这个长程吸引力就是库仑力. 所以, 库仑力是原子结合的动力.3.晶体的结合能, 晶体的内能, 原子间的相互作用势能有何区别?[解答]自由粒子结合成晶体过程中释放出的能量, 或者把晶体拆散成一个个自由粒子所需要的能量, 称为晶体的结合能.原子的动能与原子间的相互作用势能之和为晶体的内能.在0K时, 原子还存在零点振动能. 但零点振动能与原子间的相互作用势能的绝对值相比小得多. 所以, 在0K时原子间的相互作用势能的绝对值近似等于晶体的结合能.4.原子间的排斥作用取决于什么原因?[解答]相邻的原子靠得很近, 以至于它们内层闭合壳层的电子云发生重叠时, 相邻的原子间便产生巨大排斥力. 也就是说, 原子间的排斥作用来自相邻原子内层闭合壳层电子云的重叠.5.原子间的排斥作用和吸引作用有何关系? 起主导的范围是什么?[解答]在原子由分散无规的中性原子结合成规则排列的晶体过程中, 吸引力起到了主要作用. 在吸引力的作用下, 原子间的距离缩小到一定程度, 原子间才出现排斥力. 当排斥力与吸引力相等时, 晶体达到稳定结合状态. 可见, 晶体要达到稳定结合状态, 吸引力与排斥力缺一不可. 设此时相邻原子间的距离为0r, 当相邻原子间的距离r>0r时, 吸引力起主导作用; 当相邻原子间的距离r<0r时, 排斥力起主导作用.6.共价结合为什么有“饱和性”和“方向性”?[解答]设N为一个原子的价电子数目, 对于IV A、V A、VI A、VII A族元素,价电子壳层一共有8个量子态, 最多能接纳(8- N)个电子, 形成(8- N)个共价键. 这就是共价结合的“饱和性”.共价键的形成只在特定的方向上, 这些方向是配对电子波函数的对称轴方向, 在这个方向上交迭的电子云密度最大. 这就是共价结合的 “方向性”.7. 共价结合, 两原子电子云交迭产生吸引, 而原子靠近时, 电子云交迭会产生巨大的排斥力, 如何解释?[解答]共价结合, 形成共价键的配对电子, 它们的自旋方向相反, 这两个电子的电子云交迭使得体系的能量降低, 结构稳定. 但当原子靠得很近时, 原子内部满壳层电子的电子云交迭, 量子态相同的电子产生巨大的排斥力, 使得系统的能量急剧增大. 8.试解释一个中性原子吸收一个电子一定要放出能量的现象.[解答]当一个中性原子吸收一个电子变成负离子, 这个电子能稳定的进入原子的壳层中, 这个电子与原子核的库仑吸引能的绝对值一定大于它与其它电子的排斥能. 但这个电子与原子核的库仑吸引能是一负值. 也就是说, 当中性原子吸收一个电子变成负离子后, 这个离子的能量要低于中性原子原子的能量. 因此, 一个中性原子吸收一个电子一定要放出能量. 9.如何理解电负性可用电离能加亲和能来表征?[解答]使原子失去一个电子所需要的能量称为原子的电离能, 电离能的大小可用来度量原子对价电子的束缚强弱. 一个中性原子获得一个电子成为负离子所释放出来的能量称为电子亲和能. 放出来的能量越多, 这个负离子的能量越低, 说明中性原子与这个电子的结合越稳定. 也就是说, 亲和能的大小也可用来度量原子对电子的束缚强弱. 原子的电负性大小是原子吸引电子的能力大小的度量. 用电离能加亲和能来表征原子的电负性是符合电负性的定义的.10.为什么许多金属为密积结构?[解答]金属结合中, 受到最小能量原理的约束, 要求原子实与共有电子电子云间的库仑能要尽可能的低(绝对值尽可能的大). 原子实越紧凑, 原子实与共有电子电子云靠得就越紧密, 库仑能就越低. 所以, 许多金属的结构为密积结构. 11.何为杂化轨道?[解答]为了解释金刚石中碳原子具有4个等同的共价键, 1931年泡林(Pauling)和斯莱特(Slater)提出了杂化轨道理论. 碳原子有4个价电子, 它们分别对应s 2ϕ、xp 2ϕ、yp 2ϕ、zp 2ϕ量子态, 在构成共价键时, 它们组成了4个新的量子态).(21),(21),(21),(2122221222212222122221z y x z y x z y x z y x p p p s p p p s p p p s p p p s ϕϕϕϕψϕϕϕϕψϕϕϕϕψϕϕϕϕψ+−−=−+−=−−+=+++=,4个电子分别占据1ψ、2ψ、3ψ、4ψ新轨道, 在四面体顶角方向(参见图1.18)形成4个共价键.12.你认为固体的弹性强弱主要由排斥作用决定呢, 还是吸引作用决定?[解答]如上图所示, 0r 附近的力曲线越陡, 当施加一定外力, 固体的形变就越小. 0r 附近力曲线的斜率决定了固体的弹性性质. 而0r 附近力曲线的斜率主要取决于排斥力. 因此, 固体的弹性强弱主要由排斥作用决定. 13.固体呈现宏观弹性的微观本质是什么?[解答]固体受到外力作用时发生形变, 外力撤消后形变消失的性质称为固体的弹性. 设无外力时相邻原子间的距离为0r , 当相邻原子间的距离r >0r 时, 吸引力起主导作用; 当相邻原子间的距离r <0r 时, 排斥力起主导作用. 当固体受挤压时, r <0r , 原子间的排斥力抗击着这一形变. 当固体受拉伸时, r >0r , 原子间的吸引力抗击着这一形变. 因此, 固体呈现宏观弹性的微观本质是原子间存在着相互作用力, 这种作用力既包含着吸引力, 又包含着排斥力.14.你是如何理解弹性的, 当施加一定力, 形变大的弹性强呢, 还是形变小的强?[解答]对于弹性形变, 相邻原子间的距离在0r 附近变化. 令r r r ∆+=0, 则有).1(),1()1()(0000000r rnr r r rmr r rr r r r n n m m m m m ∆∆∆∆−≈−≈+=+=−−−−−−−因为0/r r ∆是相对形变, 弹性力学称为应变, 并计作S , 所以原子间的作用力.)(000000S r Bn r Am r BnS r AmS r B r A r B r A f n m n m n m n m −=−++−=+−=再令c r Bnr Am nm =−00, cS f =.可见, 当施加一定力, 形变S 大的固体c 小, 形变S 小的固体c 大. 固体的弹性是固体的属性, 它与外力和形变无关. 弹性常数c 是固体的属性, 它的大小可作为固体弹性强弱的度量. 因此, 当施加一定力, 形变大的弹性弱, 形变小的强. 从这种意义上说, 金刚石的弹性最强.15.拉伸一长棒, 任一横截面上的应力是什么方向? 压缩时, 又是什么方向?[解答]如上图所示, 在长棒中取一横截面, 长棒被拉伸时, 从截面的右边看, 应力向右, 但从截面的左边看, 应力向左. 压缩时, 如下图所示, 应力方向与拉伸时正相反. 可见, 应16.固体中某一面积元两边的应力有何关系?[解答以上题为例, 在长棒中平行于横截面取一很薄的体积元, 拉伸时体积元两边受的应力如图所示.压缩时体积元两边受的应力如下图所示.当体积元无限薄, 体积元将变成面积元. 从以上两图可以看出, 面积元两边的应力大小相等方向相反.17.沿某立方晶体一晶轴取一细长棒做拉伸实验, 忽略宽度和厚度的形变, 由此能否测出弹性劲度常数11c ?[解答]立方晶体c b a , ,轴是等价的, 设长棒方向为x (a , 或b , 或c )轴方向, 做拉伸实验时若忽略宽度和厚度的形变, 则只有应力1T 应变1S 不为0, 其它应力应变分量都为0. 由(2.55)可得 1111S c T =. 设长棒的横截面积为A , 长度为L , 拉伸力为F , 伸长量为L ∆, 则有: L L S A F T / ,/11∆==. 于是, L A FL c ∆/11=.18.若把上题等价成弹簧的形变, 弹簧受的力kx F −=, k 与11c 有何关系?[解答]上题中长棒受的力L c L AF ∆11=,长棒的伸长量L ∆即是弹簧的伸长量x . 因此,.11c L A k =可见, 弹簧的弹性系数k 与弹性劲度常数的量纲是不同的.19.固体中的应力与理想流体中的压强有何关系?[解答]固体受挤压时, 固体中的正应力321 , ,T T T 与理想流体中的压强是等价的, 但654 , ,T T T 不同于理想流体中的压强概念. 因为压强的作用力与所考虑截面垂直, 而654 , ,T T T 与所考虑截面平行. 也就是说, 理想流体中不存在与所考虑截面平行的作用力.这是因为理想流体分子间的距离比固体原子间距大得多, 流层与流层分子间不存在切向作用力.20.固体中的弹性波与理想流体中的传播的波有何差异? 为什么?[解答]理想流体中只能传播纵波. 固体中不仅能传播纵波, 还能传播切变波. 这是因为理想流体分子间距离大, 分子间不存在切向作用力, 只存在纵向作用力;而固体原子间距离小, 原子间不仅存在纵向作用力, 还存在切向作用力.。
第21. 有一晶体,平衡时体积为 0V , 原子间相互作用势为0.如果相距为 r 的两原子互作用势为 ()n m r r a r u β+-= 证明(1) 体积弹性模量为 K=.90V mnU (2) 求出体心立方结构惰性分子的体积弹性模量.[解答]设晶体共含有 N 个原子,则总能量为U(r)=()∑∑i jij r u '21. 由于晶体表面层的原子数目与晶体内原子数目相比小得多,因此可忽略它们之间的基异,于是上式简化为 U=().2'∑jijr u N设最近邻原子间的距离为R 则有j ij a r =R再令 A ,1'∑=j m j m a A ,1'∑=jn j n a 得到 U=.200⎪⎪⎭⎫ ⎝⎛+-n n m m R A R A N βα 平衡时R=R 0,则由已知条件U(R 0) = 0U 得0002U R A R A N n n m m =⎪⎪⎭⎫⎝⎛+-βα 由平衡条件 0)(0=R dRR dU得021010=⎪⎪⎭⎫⎝⎛-++n nm m R A n R A m N βα. 由(1),(2)两式可解得.)(2,)(20000n n m m nR n m N U A nR n m N U A -=-=βα利用体积弹性模量公式[参见《固体物理教程》(2.14)式]K=0220209R R U V R ⎪⎪⎭⎫ ⎝⎛∂∂得K= ⎥⎦⎤⎢⎣⎡+++-n n m m R A n n R A m m N V 000)1()1(291βα = ⎥⎦⎤⎢⎣⎡-++-+-)(2)1()(2)1(2910000000n m N mR U R n n n m N nR U R m m N V nnm m = .900V mn U - 由于,00<U 因此,00U U -= 于是 K= .90V mnU (1) 由《固体物理教程》(2.18)式可知,一对惰性气体分子的互作用能为.)(126r B r A r u +-=若令 61,42⎪⎭⎫⎝⎛==A B B A σε,则N 个惰性气体分子的互作用势能可表示为⎥⎥⎦⎤⎢⎢⎣⎡⎪⎭⎫ ⎝⎛-⎪⎭⎫ ⎝⎛=6612122)(R A R A N r U σσε.由平衡条件0)(0=R dRR dU 可得 R .2616120⎪⎪⎭⎫ ⎝⎛=A A σ进一步得 .2)(122600A A N R U U ε-==代入K=.900V mn U 并取 m =6,n =12,V 300334R N =得 K=5126123233⎪⎪⎭⎫⎝⎛A A A σε.对体心立方晶体有 A .11.9,25.12126==A 于是.1.703σε=K 2. 一维原子链,正负离子间距为a ,试证:马德隆常数为2=μ1n2. [解答] 相距ij r 的两个离子间的互作用势能可表示成.4)(2n ijij ij r br q r u +=πμ设最近邻原子间的距离为R 则有 R a r j ij =, 则总的离子间的互作用势能 U=()∑∑∑-⎪⎪⎭⎫ ⎝⎛±-=jn jn j j j ij a bRa R q N r u N ''0'114[22πε. 基中 jja 1'±=∑μ 为离子晶格的马德隆常数,式中+;- 号分别对应于与参考离子相异和相同的离子.任选一正离子作为参考离子,在求和中对负离子到正号,对正离子取负号,考虑到对一维离子两边的离子是正负对称分布的,则有.413121112)1('⎥⎦⎤⎢⎣⎡+-+-=±=∑Λj ja μ利用正面的展开式 1n(1+x ),432432Λ+-+-x x x x 并令 1=x 得Λ+-+-41312111=1n(1+1)=1n2.于是,一维离子链的马德常数为2=μ1n23. 计算面心立方面简单格子的6A 和12A(1) 只计最近邻; (2) 计算到次近邻; (3) 计算到次近邻.[解答]图2.26示出了面心立方简单格子的一个晶胞.角顶O 原子周围有8个这样的晶胞,标号为1的原子是原子O 的最近邻标号为2的原子是O 原子的最近邻,标号为3的原子是O 原子的次次近邻.由此得到,面心立方简单格子任一原子有12个最近邻,6个次近邻及24个次次近邻.以最近邻距离度量,其距离分别为:.3,2,1===j j j a a a 由 .1,112'126'6⎪⎪⎭⎫⎝⎛=⎪⎪⎭⎫ ⎝⎛=∑∑jj j j a A a A图2.6 面心立方晶胞得(1) 只计最近邻时1211*12)1(66=⎪⎭⎫⎝⎛=A , 1211*12)1(1212\=⎪⎭⎫⎝⎛=A .(2) 计算到次近邻时.094.1221*611*12)2(,750.1221*611*12)2(121212666=⎪⎭⎫ ⎝⎛+⎪⎭⎫⎝⎛==⎪⎭⎫ ⎝⎛+⎪⎭⎫⎝⎛=A A(3) 计算到次次近邻时.127.12033.0094.1231*2421*611*12)3(,639.13899.0750.1231*2421*611*12)3(121212126666=+=⎪⎪⎭⎫ ⎝⎛+⎪⎭⎫ ⎝⎛+⎪⎭⎫⎝⎛==+=⎪⎪⎭⎫ ⎝⎛+⎪⎭⎫ ⎝⎛+⎪⎭⎫⎝⎛=A A 由以上可以看出,由于12A中的幂指数较大,12A 收敛得很快,而6A 中的幂指数较小,因此 6A 收敛得较慢,通常所采用的面心立方简单格子的 6A 和 12A 的数值分别是14.45与12.13.4. 用埃夫琴方法计算二维正方离子(正负两种)格子的马德隆常数. [解答]马德隆常数的定义式为 jja 1'±=∑μ,式中+、-号分别对应于与参考离子相异和相同的离子,二维正方离子(正负两种)格子,实际是一个面心正方格子,图 2.7示出了一个埃夫琴晶胞.设参考离子O 为正离子,位于边棱中点的离子为负离子,它们对晶胞的贡献为4*(1/2).对参考离子库仑能的贡献为图2.7二维正方离子晶格.121*4顶角上的离子为正离子,它们对晶胞的贡献为4*(1/4), 对参考离子库仑能的贡献为 .241*4-因此通过一个埃夫琴晶胞算出的马德隆常数为 .293.1241*4121*4=-=ν再选取422=个埃夫琴晶胞作为考虑对象,这时离子O 的最的邻,次近邻均在所考虑的范围内,它们对库仑能的贡献为,2414-而边棱上的离子对库仑能的贡献为 ,521*8221*4+- 顶角上的离子对为库仑能的贡献为 ,841*4-这时算出的马德隆常数为图 2.8 4个埃夫琴晶胞同理对932=个埃夫琴晶胞进行计算,所得结果为611.11841*41321*81021*8321*48458242414=⎪⎪⎪⎪⎭⎫ ⎝⎛-+-+⎪⎪⎭⎫ ⎝⎛-+-+⎪⎭⎫ ⎝⎛-=μ 对 1642=个埃夫琴晶胞进行计算,所得结果为614.13241*42521*81721*81021*8421*4184138108348458242414=⎪⎪⎪⎪⎭⎫ ⎝⎛-+-+-+⎪⎪⎭⎫ ⎝⎛-+-+⎪⎪⎭⎫ ⎝⎛-+-+⎪⎭⎫ ⎝⎛-=μ当选取 n 2个埃夫琴晶胞来计算二维正方离子(正负两种)格子的马德隆常数,其计算公式(参见刘策军,二维NaC1 晶体马德隆常数计算,《大学物理》,Vo1.14,No.12,1995.)为 [][].1,8411>+++=--n D C B A n n n n μ其中 ,21)1(,1)1(11111nB t A n n n t t n +-=+--=-=∑,1)1(1)1()2()1(1)1()1(2112212221112122122222222221⎪⎪⎪⎪⎪⎭⎫⎝⎛+--+--+-+-+--++⎪⎪⎭⎫⎝⎛+++-+⎪⎪⎭⎫ ⎝⎛+---n n n n n C n n ΛΛ.121)1()1(2181222222+-+-++++-=n n n n n D n n Λ5. 用埃夫琴方法计算CsCl 型离子晶体的马德隆常数(1) 只计最近邻 (2) 取八个晶胞 [解答](1) 图2.29是CsCl 晶胸结构,即只计及最近邻的最小埃夫琴晶胞,图2.29()a 是将Cs +双在体心位置的结构,图2.9(a)是将 Cl -取在体心位置的结构,容易求得在只计及最近邻情况下,马德隆常数为1.图2.29 (a )Cs 取为体心的CsC1晶胞,(b) C1取为体心的CsC1晶胞(2)图2.10是由8个CsCl 晶胞构成的埃夫琴晶胞,8个最近邻在埃夫琴晶胞内,每个离子对晶胞的贡献为1,它们与参考离子异号,所以这8个离子对马德隆常数的贡献为8埃夫琴晶胞6个面上的离子与参考离子同号,它们对埃夫琴晶胞的贡献是21,它们与参考离子的距离为32R 它们对马德隆常数的贡献为-()3/2*621图 2.10 8个CsCl 晶胞构成的一个埃夫琴晶胞埃夫琴晶胞楞上的12个离子,与参考离子同号,它们对埃夫琴晶胞的贡献是41它们与参考离子的距离为322R 它们对马德隆常数的贡献为-()3224/1*12埃夫琴晶胞角顶上的 8个离子,与参考离子同号,它们对埃夫琴晶胞的贡献是81它们与参考离子的距离为2R 它们对马德隆常数的贡献为 -()281*8,由8个CsCl 晶胞构成的埃夫琴晶胞计算的马德隆常数.064806.32)8/1(*8322)4/1(*123/2)2/1(*68=---=μ 为了进一步找到马德常数的规律,我们以计算了由27个CsCl 晶胞构成的埃夫琴晶胞的马德隆常数,结果发现,由27个CsCl 晶胞构成的埃夫琴晶胞的马德隆常数是0.439665.马德隆常数的不收敛,说明CsCl 晶胞的结构的马德隆常数不能用传统的埃夫琴方法计算.为了找出合理的计算方法,必须首先找出采用单个埃夫琴晶胞时马德隆常数不收敛的原因.为了便于计算,通常取参考离子处于埃夫琴晶胞的中心.如果以Cs +作参考离子,由于埃夫琴晶胞是电中性的要求,则边长为pa 2(p 是大于或等于1的整数)的埃夫琴晶胞是由(2p )3个CsCl 晶胞所构成,埃夫琴晶胞最外层的离子与参考离子同号,而边长为(2p +1)的埃夫琴晶胞是由(2p +1)3 个 CsCl 晶胞所构成,但埃夫琴晶胞的最外层离子与参考离子异号,如果以C1-作参考离子也有同样的规律,设参考离子处于坐标原点O ,沿与晶胞垂直的方向(分别取为x,y,z 图2.11示出了z 轴)看去,与参考郭同号的离子都分布在距O 点ia 的层面上,其中i 是大于等于 1的整数,与 O 点离子异号的离子都分布在距O 点(i -0.5)a 的层面上,图 2.11(a) 示出了同号离子层,图2.11(b)示出了异号离子层.图2.11 离子层示意图(a)表示同号离子层, O 离子所在层与 O '离子所在层相距ia(b)表示异号离子层, O 离子所在层和O ' 离子所在层相距(i -0.5)a当 CsCl 埃夫琴晶胞边长很大时,晶胞最外层的任一个离子对参考离子的库仑能都变得很小,但它们对参考离子总的库仑能不能忽略.对于由(2p )3个CsCl 晶胞所构成的埃夫琴晶胞来说,最外层有6*(2p )2个与参考离子同号的离子,它们与参考离子的距离为(1/2)pa ~(23)pa ,它们与参考离子的库仑能为a pe 024πε量级,这是一个相对大的正值.对于由(2p +1)3个CsCl 晶胞所构成的埃夫琴晶胞来说,离外层有6*(2p +1)2个与参考离子异号的离子,它们与参考离子的库仑能为a pe 024πε-量级,这是一个绝对值相对大的负值,因此,由(2p )3个CsCl 晶胞构成的埃夫琴晶胞所计算的库仑能,与由(2p +1)3个CsCl 晶胞构成的埃夫琴晶胞所计算的库仑能会有较大的差异.即每一情况计算的库仑能都不能代表CsCl 晶体离子间相互作用的库仑能.因此这两种情况所计算的马德隆常数也必定有较大的差异,由1个CsCl 晶胞、8个CsCl 晶胞和27个CsCl 晶胞构成的埃夫琴晶胞的计算可知, CsCl 埃夫琴晶胞体积不大时,这种现象已经存在.为了克服埃夫琴方法在计算马德隆常数时的局限性,可采取以下方法,令由 (2p )3个CsCl 晶胞构成的埃夫琴晶胞计算的库仑能为1U ,由(2p +1)3个CsCl 晶胞构成的埃夫琴晶胞所计算的库仑能为1U ,则CsCl 晶体离子间相互作用的库仑能可近似取作 )(2121U U U +=(1) 因子1/2 的引入是考虑除了(2p +1)3个CsCl 晶胞构成的埃夫琴晶胞最外层离子外,其他离子间的库仑能都累计了两偏,计算1U 和2U 时要选取体积足够大的埃夫琴晶胞,此时埃夫琴晶胞最外层离子数与晶胞内的离子数相比是个很小的数,相应的马德隆常数应为 )(2121μμμ+=(2) 其中:=1μ⎪⎪⎭⎫ ⎝⎛±∑i ja 1'是由(2p )3个CsC1晶胞构成的埃夫琴晶胞计算的值; =1μ⎪⎪⎭⎫ ⎝⎛±∑i ja 1'由 (2p +1)3 个CsC1晶胞构成的埃夫琴晶胞所计算成本的值.为简化计算,特选取晶胞边长a 为计算单位,由于,32a R =所以,23'μμ= ⎪⎪⎭⎫ ⎝⎛±=∑'''1i i a μ (3) 其中'i a 是某一离子到参点的距离与a 的比值.考虑到对称性,对选定的埃夫琴晶胞,把晶胞的离子看成分布在一个个以参考离子为对称心的正六面体的六个面上,体积不同的正六面六个面上的离子分别计算.由(2p )3个CsC1晶胞构成埃夫琴晶胞时,由分析整理可得,231111⎪⎪⎭⎫ ⎝⎛++=∑∑=-=p pi i p i i C B A μ (4) 由(2p +1)3个 CsC1 晶胸构成埃夫琴晶胞时,,231112⎪⎪⎭⎫ ⎝⎛++=∑∑=-=p pi i p i i D B A μ (5)其中:),1(''''22'2'p i i y x k A i x iy y x i <≤++-=∑∑(6)i A 表示与 O 点距离为ia 的6个面上所有的离子对马德隆常数的面贡献,因为这些离子与参考离子同号,故到负号.'x 、'y 是离子在平面 '''y x o 上的坐标, ''y x k 代表 6个面上等价离子的个数,其取值规则为:(1) 在角上(如E 点),即'x =i 且 'y = i. 时, ''y x k =8;(2) 在棱与坐标轴的交点(如 F 点),'x =i 且'y = 0或 'x =0且'y = 0时, ''y x k =6 (3) 在棱上的其他点(如H 、I 点)即不满足上述条件,且'x =i 或'y = i.时, ''y x k =12 (4) 在'O 点,即'x =0且'y = 0时, ''y x k =6(5) 在除'O 点外的面上的点(如J 点),即不满足上述条件时,''y x k =24.),1()5.0(5.05.05.05.022'2''''''p i i y x k B i x i y yx i ≤≤-++=∑∑-=-=(7)i B 代表距O 点距离为(i -0.5)a 的6个面上的离子对马德隆常数的贡献,因为这种些离子与参考离子异号,故取正号. 'x ,'y 是离子在平面'''y x o 上的坐标, '''y x k 代表这6个面上等价离子的个数,其取值规则为:(1) 在角上(如K 点),即'x =i 且 'y = i.时, '''y x k =8;(2) 在棱下(如L 、M 点),即不满足不述条件,且'x =i 或'y = i 时,'''y x k =12; (3) 在面上(如N 点)好不满足上述条件时, '''y x k =24.),(0022'2'"''''p i i y x k C i x iy i yx =++-=∑∑==i C 表示在边长为2pa 的晶胞最外层,即与参考离子相距pa 的6个面上的离子对马德隆常数的贡献,应取负号,与iA 的不同在于"''y x k的取值: (1) 在角上, "''y x k =''y x k /8; (2) 在棱上, "''y x k =''y x k /4; (3) 在面上, "''y x k=''y x k /2.),()5.0(5.05.05.05.022'2''''''''p i i y x k D i x i y yx i =-++=∑∑-=-=i D 表示在边长为2a p )1(+的晶胞最外层,即与参考离子相距(p +0.5)a 的离子层对马德隆常数的贡献,应取正号,与i B 的不同在于'''''yx k 的取值: (1) 在角上, '''''y x k ='''y x k /8; (2) 在棱上, '''''y x k ='''y x k /4; (3) 在面上, '''''y x k ='''y x k /2.表2.1给出了计算结果,给出的μ是由分别对应2p 和2p+1的1μ和2μ求得的,实际上, 1μ和2μ只需对应边长相近的埃夫琴晶胞即可,如取对应2p 和2p-1的埃夫琴晶胞也可得到一样的收敛结果,由以上数据可见,马德隆常数μ随晶胞边长的增大而迅速收敛.该方法适用于NaC1结构以外离子晶体马德隆常数的计算.6.只计及最近邻间的排斥作用时,一离子晶体离子间的互作用势为⎪⎪⎩⎪⎪⎨⎧±-=-)2(,)1(,)(22r e R e e r u R ρλ(1)最近邻(2)最近邻以外 式中ρλ,是常数,R 是最近邻距离,求晶体平衡时,原子间总的互作用势.[解 答]设离子数目为2N,以j ij a r =R 表示第j 个离子到参考离子i 的距离,忽略表面效应,则总的相互作用能可表示为U =N ⎥⎥⎦⎤⎢⎢⎣⎡+⎪⎪⎭⎫ ⎝⎛±-∑∑-ρλR j j e R a e 2' (∑表示最近邻)=N ,2⎥⎦⎤⎢⎣⎡+--ρλμR e Z R e其中⎪⎪⎭⎫⎝⎛±=∑j ia 1'μ 为马德隆常数,+号对应于异号离子,-号对应于同号离子;Z 为任一离子的最近邻数目,设平衡时R=R 0 ,由平衡条件,02020=⎥⎦⎤⎢⎣⎡+=-ρρλμR R e Z R e N dRdU 得.0202ρλμρR e Z R e -=平衡时的总相互作用为.1)(0020200⎪⎪⎭⎫⎝⎛-=⎥⎦⎤⎢⎣⎡+-=-R R e N e Z R e N R U R ρμλμρ 7. 设离子晶体中,离子间的互作用势为⎪⎪⎩⎪⎪⎨⎧±+-=最近邻以外最近邻,,)(22re R b R e r u m(1) 求晶体平衡时,离子间总的相互作用势能)(0R U (2) 证明: )(0R U 11-⎪⎪⎭⎫⎝⎛∝m mZ μ其中μ是马德隆常数,Z 是晶体配位数 [解答](1)设离子数目为2N , 以j ij a r =R 表示第j 个离子到参考离子i 的距离,忽略表面效应,则总的相互作用能可表示U =N ⎥⎥⎦⎤⎢⎢⎣⎡+⎪⎪⎭⎫ ⎝⎛±-∑∑m j j R b R a e 2'(∑表示最近邻) =N ,2⎥⎦⎤⎢⎣⎡+-m R b Z Re μ其中⎪⎪⎭⎫ ⎝⎛±=∑j i a 1'μ,为马德隆常数,+号对应于异号离子,-号对应于同号离子.Z 为任一离子的最近邻数目,设平衡时R=R 0由平衡条件,0102020=⎥⎦⎤⎢⎣⎡-=+m R R Zmb R e N drdUμ得10-m R Zmb=2e μ即1120-⎪⎪⎭⎫ ⎝⎛=m e Zmb R μ.于是,晶体平衡时离子间总的相互作用势能0U =).1(000--=⎥⎦⎤⎢⎣⎡+-m R NZbR b Z R Zmb N m m m(2)晶体平衡时离子间的相互作用势能可进一步化为0U =.)()()1()1(1111121211--------=⎪⎪⎭⎫ ⎝⎛--m m m m mm m m m m mb Ze Nbm e Zmb ZNbm μμ由上式可知 .110-⎪⎪⎭⎫⎝⎛∝m mZ U μ8.一维离子链,其上等间距载有正负2N 个离子,设离子间的泡利排斥只出现在最近邻离子之间,且为b/R n,b,n 是常R 是两最近邻离子的间距,设离子电荷为q ,(1) 试证明平衡间距下 )(0R U =;114212002⎪⎭⎫⎝⎛--n R n Nq πε(2) 令晶体被压缩,使0R )1(0δ-→R , 试证明在晶体被压缩单位长度的过程中外力作功的主项为c 2δ其中c=;21)1(02R n q n -(3) 求原子链被压缩了2)1(0<<e e NR δδ时的外力[解答](1) 因为离子间是等间距的,且都等于R ,所以认定离子与第j 个离子的距离j r 总可表示成为R a r j j =ja 是一整数,于是离子间总的互作用势能⎥⎥⎦⎤⎢⎢⎣⎡-⎪⎪⎭⎫ ⎝⎛±-=⎥⎥⎦⎤⎢⎢⎣⎡+=∑∑n i in j j j R b a R q N r b r q N R U 214242)('202'0πεπεμ,其中+、-分别对应相异离子和相同离子的相互作用.一维离子晶格的马德隆常数(参见本章习题2)为=⎪⎪⎭⎫ ⎝⎛±∑i ia 1'21n2. 利用平衡条件0)(0=R dRR dU得到b=nq 01-n 0241n2R πε,)(R U =⎪⎪⎭⎫ ⎝⎛---n n nR R R Nq 102141n22πε. 在平衡间距下⎪⎭⎫⎝⎛--n R Nq R U 1141n22)(0020πε.(2) 将互作用势能在平衡间距附近展成级数Λ+-⎪⎪⎭⎫ ⎝⎛+-⎪⎭⎫⎝⎛+=202200)(21)()()(0R R dR U d R R dR dU R U R U R R 由外力作的功等于晶体内能的增量,可得外力作功的主项为W=20220)(21)()(0R R dR U d R U R U R-⎪⎪⎭⎫ ⎝⎛=-, 其中利用平衡条件,将R=R )1(0δ- ,代入上式,得到W=δδπε)2(421)1(2102002NR R n q n ⎥⎥⎦⎤⎢⎢⎣⎡-. 晶体被压缩单位长度的过程中,外力作的功的主项δ02W NR =δπε⎥⎥⎦⎤⎢⎢⎣⎡-2002421)1(21R n q n 令c=202421)1(R n q n πε-(CGS)得到在晶体被压缩单位长度的过程中,外力作的功的主项为2δc . (3)设e δδ=时外力为F e ,由于在弹性范围内,外力与晶格的形变成正比,所以 F= )2(0δαNR , F e = )2(0e NR δα,其中α为比例系数离子链被压缩e NR δ02过程中外力作的功W e =δδαδδd NR NR Fdx e eNR e 020002)]2([0⎰⎰== e e e F NR NR δδα022022121)2(=.由于 W e =)2(20e eNR c δδ,所以离子链被压缩了e NR δ02时的外力为F e =202)1(21R n n q c ee δδ-=.9.设泡利排斥项的形式不变,讨论电荷加倍对NaC1晶格常数,体积弹性模量以及结合能的影响。
固体物理学答案朱建国版3HUA system office room 【HUA16H-TTMS2A-HUAS8Q8-HUAH1688】固体物理学·习题指导配合《固体物理学(朱建国等编着)》使用2022年4月28日第1章晶体结构 0第2章晶体的结合 (11)第3章晶格振动和晶体的热学性质 (17)第4章晶体缺陷 (26)第5章金属电子论 (30)第1章 晶体结构1.1 有许多金属即可形成体心立方结构,也可以形成面心立方结构。
从一种结构转变为另一种结构时体积变化很小.设体积的变化可以忽略,并以R f 和R b 代表面心立方和体心立方结构中最近邻原子间的距离,试问R f /R b 等于 多少?答:由题意已知,面心、体心立方结构同一棱边相邻原子的距离相等,都设为a :对于面心立方,处于 面心的原子与顶角原子的距离为:R f =2a对于体心立方,处于体心的原子与顶角原子的距离为:R b那么,RfRb =31.2 晶面指数为(123)的晶面ABC 是离原点O 最近的晶面,OA 、OB 和OC 分别与基失a 1,a 2和a 3重合,除O 点外,OA ,OB 和OC 上是否有格点若ABC 面的指数为(234),情况又如何答:晶面族(123)截a 1,a 2,a 3分别为1,2,3等份,ABC 面是离原点O 最近的晶面,OA 的长度等于a 1的长度,OB 的长度等于a 2长度的1/2,OC 的长度等于a 3长度的1/3,所以只有A 点是格点。
若ABC 面的指数为(234)的晶面族,则A 、B 和C 都不是格点。
1.3 二维布拉维点阵只有5种,试列举并画图表示之。
答:二维布拉维点阵只有五种类型,两晶轴b a 、,夹角 ,如下表所示。
1 简单斜方2 简单正方3 简单六角4 简单长方5 有心长方二维布拉维点阵1.4 在六方晶系中,晶面常用4个指数(hkil )来表示,如图所示,前3个指数表示晶面族中最靠近原点的晶面在互成120°的共平面轴a 1,a 2,a 3上的截距a 1/h ,a 2/k ,a 3/i ,第四个指数表示该晶面的六重轴c 上的截距c/l.证明:i=-(h+k ) 并将下列用(hkl )表示的晶面改用(hkil )表示:(001)(133)(110)(323)(100)(010)(213) 答:证明设晶面族(hkil )的晶面间距为d ,晶面法线方向的单位矢量为n °。
黄昆版固体物理学课后答案解析答案Prepared on 24 November 2020《固体物理学》习题解答黄昆 原着 韩汝琦改编 (陈志远解答,仅供参考)第一章 晶体结构、解:实验表明,很多元素的原子或离子都具有或接近于球形对称结构。
因此,可以把这些原子或离子构成的晶体看作是很多刚性球紧密堆积而成。
这样,一个单原子的晶体原胞就可以看作是相同的小球按点阵排列堆积起来的。
它的空间利用率就是这个晶体原胞所包含的点的数目n 和小球体积V 所得到的小球总体积nV 与晶体原胞体积Vc 之比,即:晶体原胞的空间利用率, VcnV x =(1)对于简立方结构:(见教材P2图1-1) a=2r , V=3r 34π,Vc=a 3,n=1∴52.06r8r 34a r 34x 3333=π=π=π=(2)对于体心立方:晶胞的体对角线BG=x 334a r 4a 3=⇒= n=2, Vc=a 3∴68.083)r 334(r 342a r 342x 3333≈π=π⨯=π⨯=(3)对于面心立方:晶胞面对角线BC=r 22a ,r 4a 2=⇒=n=4,Vc=a 3(4)对于六角密排:a=2r 晶胞面积:S=6260sin a a 6S ABO ⨯⨯=⨯∆=2a 233晶胞的体积:V=332r 224a 23a 38a 233C S ==⨯=⨯ n=1232126112+⨯+⨯=6个(5)对于金刚石结构,晶胞的体对角线BG=3r 8a r 24a 3=⇒⨯= n=8, Vc=a 3、试证:六方密排堆积结构中633.1)38(ac 2/1≈=证明:在六角密堆积结构中,第一层硬球A 、B 、O 的中心联线形成一个边长a=2r 的正三角形,第二层硬球N 位于球ABO 所围间隙的正上方并与这三个球相切,于是: NA=NB=NO=a=2R.即图中NABO 构成一个正四面体。
…、证明:面心立方的倒格子是体心立方;体心立方的倒格子是面心立方。
固体物理(胡安)第⼆版课后习题答案__完整版_校核版Word 版完整版校核版第⼀章晶体的结构及其对称性1.1⽯墨层中的碳原⼦排列成如图所⽰的六⾓⽹状结构,试问它是简单还是复式格⼦。
为什么?作出这⼀结构所对应的两维点阵和初基元胞。
解:⽯墨层中原⼦排成的六⾓⽹状结构是复式格⼦。
因为如图点A 和点B 的格点在晶格结构中所处的地位不同,并不完全等价,平移A →B,平移后晶格结构不能完全复原所以是复式格⼦。
1.2在正交直⾓坐标系中,若⽮量k l j l i l R l321 ,i ,j ,k 为单位向量。
3,2,1 i l i 为整数。
问下列情况属于什么点阵?(a )当i l为全奇或全偶时;(b )当i l之和为偶数时。
解: 112233123l R l a l a l a l i l j l kr r r r r r r...2,1,0,,321 l l l 当l 为全奇或全偶时为⾯⼼⽴⽅结构点阵,当321l l l 之和为偶数时是⾯⼼⽴⽅结构 1.3 在上题中若321l l l 奇数位上有负离⼦,321l l l 偶数位上有正离⼦,问这⼀离⼦晶体属于什么结构?解:是离⼦晶体,属于氯化钠结构。
1.4 (a )分别证明,⾯⼼⽴⽅(fcc )和体⼼⽴⽅(bcc )点阵的惯⽤初基元胞三基⽮间夹⾓相等,对fcc 为,对bcc 为(b )在⾦刚⽯结构中,作任意原⼦与其四个最近邻原⼦的连线。
证明任意两条线之间夹⾓θ均为'1cos 109273arc o '1cos 109273arco解:(1)对于⾯⼼⽴⽅ 12a a j k r r r 22a a i k r r r32a a i j r r r13222a a a a r r r1212121602a a COS a a a a o r rr r2323231602a a COS a a a a o r rr r1360COS a a o r r(2)对于体⼼⽴⽅ 12a a i j k r r r r 22a a i j k r r r r32a a i j k r r r r12332a a a a r r r12'12121129273a a COS a a a a o r rr r'1313131129273a a COS a a a a o r rr r r r'2312927COS a a o r r(3)对于⾦刚⽯晶胞134a i j k rr r r234a i j k r r r r2212122122314934a COS a r rr r1.5 证明:在六⾓晶系中密勒指数为(h,k,l )的晶⾯族间距为212222234c l a k hk h d证明: a b a r r元胞基⽮的体积a ai r rcos60cos301322b a i j ai ajo o r r r r rc ck r r20033022200a a a a c c倒格⼦基⽮ )33(2][2j i a c b ajaa c b334][2k c b a c2][2倒格⽮:***hkl G ha kb lc r r r r晶⾯间距***222cl b k a h Gd hklhkl2222222222ha kb lch a k b l c hk a b kl b c hl a cr r r r r r r r r 22423a a r 22423b a r 222c cr 2223a b ar r 0b c r r 0a cr r 122222222122222242424242333343hkld h k l hk a a a a h k kl l a c1.6 证明:底⼼正交的倒点阵仍为底⼼正交的。
第一章 晶体结构1.1、 如果将等体积球分别排成下列结构,设x 表示钢球所占体积与总体积之比,证明:结构 X简单立方52.06=π体心立方68.083≈π 面心立方74.062≈π 六角密排74.062≈π 金刚石34.063≈π解:实验表明,很多元素的原子或离子都具有或接近于球形对称结构。
因此,可以把这些原子或离子构成的晶体看作是很多刚性球紧密堆积而成。
这样,一个单原子的晶体原胞就可以看作是相同的小球按点阵排列堆积起来的。
它的空间利用率就是这个晶体原胞所包含的点的数目n 和小球体积V 所得到的小球总体积nV 与晶体原胞体积Vc 之比,即:晶体原胞的空间利用率, VcnVx = (1)对于简立方结构:(见教材P2图1-1) a=2r , V=3r 34π,Vc=a 3,n=1 ∴52.06834343333====πππrra r x (2)对于体心立方:晶胞的体对角线BG=x 334a r 4a 3=⇒= n=2, Vc=a 3∴68.083)334(3423423333≈=⨯=⨯=πππr r a r x (3)对于面心立方:晶胞面对角线BC=r 22a ,r 4a 2=⇒= n=4,Vc=a 374.062)22(3443443333≈=⨯=⨯=πππr r a r x (4)对于六角密排:a=2r 晶胞面积:S=6260sin a a 6S ABO ⨯⨯=⨯∆=2a 233晶胞的体积:V=332r 224a 23a 38a 233C S ==⨯=⨯ n=1232126112+⨯+⨯=6个 74.062)22(3443443333≈=⨯=⨯=πππr r a r x (5)对于金刚石结构,晶胞的体对角线BG=3r 8a r 24a 3=⇒⨯= n=8, Vc=a 334.06333834834833333≈=⨯=⨯=πππr r a r x 1.3、证明:面心立方的倒格子是体心立方;体心立方的倒格子是面心立方。
固体物理第二章习题参考答案1.已知某晶体两相邻原子间的互作用能可表示成 n m rb r a r U +-=)( (1) 求出晶体平衡时两原子间的距离;(2) 平衡时的二原子间的互作用能;(3) 若取m=2,n=10,两原子间的平衡距离为3Å,仅考虑二原子间互作用则离解能为4ev ,计算a 及b 的值; (4) 若把互作用势中排斥项b/r n 改用玻恩-梅叶表达式λexp(-r/p),并认为在平衡时对互作用势能具有相同的贡献,求n 和p 间的关系。
解:(1)平衡时 01100=-=∂∂----n m r bnr amr r u 得 am bn r mn =-0m n ambn r -=1)(0 (2)平衡时 把r 0表示式代入u(r) u(r 0)=-m n n m n m ambn b am bn a --+)()(=-m n n m b amn a bn m -+)()( (3)由r 0表示式得: 81)5(10310ab =⨯-若理解为互作用势能为二原子平衡时系统所具有的能量,由能量最小原理,平衡时系统能量具有极小值,且为负值;离解能和结合能为要把二原子拉开,外力所作的功,为正值,所以,离解能=结合能=-互作用势能,由U(r)式的负值,得101021019)103()103(106.14---⨯-⨯+=⨯⨯ba 化简为 80101039104.6-⨯-+=⨯ba 略去第二项 a=5.76⨯102 上式代入a 值得 b=7.55⨯10-75 (4)由题意得 λex ρ(-r 0/ρ)=br -n *ln λ-r 0/ρ=lnb -nlnr 0 nlnr o =r 0/ρ+lnb/λln ln 0r n b r λρ+=又解:*式两边对r 0求导,得:λ/ρ×ex ρ(-r 0/ρ)=bnr -n+1, 与*式比较得: n/r 0 =1/ρ 得:r 0 = n ρ2.N 对离子组成的Nacl 晶体相互作用势能为 ⎥⎦⎤⎢⎣⎡-=R e RB N R U n 024)(πεα(1) 证明平衡原子间距为 n e BR n 2014απε=- (2) 证明平衡时的互作用势能为 )11(4)(0020nR Ne R U --=πεα(3) 若试验试验测得Nacl 晶体的结合能为765kj/mol,晶格常数为5.63⨯10-10m ,计算Nacl晶体的排斥能的幂指数n ,已知Nacl 晶体的马德隆常数是α=1.75证: (1)⎥⎦⎤⎢⎣⎡---=---2021)1(4)(R e e n B N dr du n πεα)4(1202+--=n RBnR e N πεα令00==R R dRdu得20104e BnR n απε=- 证毕 (2)把以上结果代入U(R)式,并把R 取为R 0⎥⎥⎦⎤⎢⎢⎣⎡-=12012020)(4)()(402440e Bn e Bn e Bn e B N R U απεαπεαπεπεα =-N )11(4002nR e -πεα若认为结合能与互作用能符号相反,则上式乘“-” 证毕(3)由(2)之结论 整理可得)(400022R U R e N e N n πεαα+= 式中 阿氏常数N =6.0⨯1023 电子电量 e=1.6⨯10-19库仑 真空介电常数 ε0=8.85⨯10-12法/米若题中R 0为异种原子的间矩,R 0=0.5×5.63⨯10-10m U(R 0)=-765000j/mol (平衡时互作用势能取极小值,且为负,而结合能为正值) 马德隆常数α=1.7520)(411e N R U R n απε-=8.811510121056.275.1100.61065,71082.21085.814.34≈-=--⨯⨯⨯⨯⨯⨯⨯⨯⨯⨯⨯。
固体物理学·习题指导配合《固体物理学(朱建国等编著)》使用2020年6月21日第1章晶体结构 0第2章晶体的结合 (13)第3章晶格振动和晶体的热学性质 (22)第4章晶体缺陷 (35)第5章金属电子论 (39)第1章晶体结构有许多金属即可形成体心立方结构,也可以形成面心立方结构。
从一种结构转变为另一种结构时体积变化很小.设体积的变化可以忽略,并以R f和R b代表面心立方和体心立方结构中最近邻原子间的距离,试问R f/R b等于多少答:由题意已知,面心、体心立方结构同一棱边相邻原子的距离相等,都设为a:对于面心立方,处于面心的原子与顶角原子的距离为:R f=2a对于体心立方,处于体心的原子与顶角原子的距离为:R b=2a那么,Rf Rb晶面指数为(123)的晶面ABC是离原点O最近的晶面,OA、OB和OC分别与基失a1,a2和a3重合,除O点外,OA,OB和OC上是否有格点若ABC面的指数为(234),情况又如何答:晶面族(123)截a1,a2,a3分别为1,2,3等份,ABC面是离原点O最近的晶面,OA的长度等于a1的长度,OB的长度等于a2长度的1/2,OC的长度等于a3长度的1/3,所以只有A 点是格点。
若ABC面的指数为(234)的晶面族,则A、B和C都不是格点。
二维布拉维点阵只有5种,试列举并画图表示之。
答:二维布拉维点阵只有五种类型,两晶轴ba、,夹角ϕ,如下表所示。
4长方2,πϕ=≠ba简单长方(图中4所示)有心长方(图中5所示)1mm,2mm1 简单斜方2 简单正方3 简单六角4 简单长方5 有心长方二维布拉维点阵在六方晶系中,晶面常用4个指数(hkil)来表示,如图所示,前3个指数表示晶面族中最靠近原点的晶面在互成120°的共平面轴a1,a2,a3上的截距a1/h,a2/k,a3/i,第四个指数表示该晶面的六重轴c上的截距c/l.证明:i=-(h+k)并将下列用(hkl)表示的晶面改用(hkil)表示:(001)(133)(110)(323)(100)(010)(213)答:证明设晶面族(hkil)的晶面间距为d,晶面法线方向的单位矢量为n°。
第一章 晶体结构1.1、 如果将等体积球分别排成下列结构,设x 表示钢球所占体积与总体积之比,证明:结构 X简单立方52.06=π体心立方68.083≈π 面心立方74.062≈π 六角密排74.062≈π 金刚石34.063≈π解:实验表明,很多元素的原子或离子都具有或接近于球形对称结构。
因此,可以把这些原子或离子构成的晶体看作是很多刚性球紧密堆积而成。
这样,一个单原子的晶体原胞就可以看作是相同的小球按点阵排列堆积起来的。
它的空间利用率就是这个晶体原胞所包含的点的数目n 和小球体积V 所得到的小球总体积nV 与晶体原胞体积Vc 之比,即:晶体原胞的空间利用率, VcnVx = (1)对于简立方结构:(见教材P2图1-1) a=2r , V=3r 34π,Vc=a 3,n=1 ∴52.06834343333====πππrra r x (2)对于体心立方:晶胞的体对角线BG=x 334a r 4a 3=⇒= n=2, Vc=a 3∴68.083)334(3423423333≈=⨯=⨯=πππr r a r x (3)对于面心立方:晶胞面对角线BC=r 22a ,r 4a 2=⇒= n=4,Vc=a 374.062)22(3443443333≈=⨯=⨯=πππr r a r x (4)对于六角密排:a=2r 晶胞面积:S=6260sin a a 6S ABO ⨯⨯=⨯∆=2a 233晶胞的体积:V=332r 224a 23a 38a 233C S ==⨯=⨯ n=1232126112+⨯+⨯=6个 74.062)22(3443443333≈=⨯=⨯=πππr r a r x (5)对于金刚石结构,晶胞的体对角线BG=3r 8a r 24a 3=⇒⨯= n=8, Vc=a 334.06333834834833333≈=⨯=⨯=πππr r a r x 1.3、证明:面心立方的倒格子是体心立方;体心立方的倒格子是面心立方。
第一章 晶体结构1.1、解:实验表明,很多元素的原子或离子都具有或接近于球形对称结构。
因此,可以把这些原子或离子构成的晶体看作是很多刚性球紧密堆积而成。
这样,一个单原子的晶体原胞就可以看作是相同的小球按点阵排列堆积起来的。
它的空间利用率就是这个晶体原胞所包含的点的数目n 和小球体积V 所得到的小球总体积nV 与晶体原胞体积Vc 之比,即:晶体原胞的空间利用率, VcnVx =(1)对于简立方结构:(见教材P2图1-1)a=2r , V=3r 34π,Vc=a 3,n=1∴52.06r8r34a r 34x 3333=π=π=π= (2)对于体心立方:晶胞的体对角线BG=x 334a r 4a 3=⇒= n=2, Vc=a 3∴68.083)r 334(r 342a r 342x 3333≈π=π⨯=π⨯= (3)对于面心立方:晶胞面对角线BC=r 22a ,r 4a 2=⇒= n=4,Vc=a 374.062)r 22(r 344a r 344x 3333≈π=π⨯=π⨯= (4)对于六角密排:a=2r 晶胞面积:S=6260sin a a 6S ABO ⨯⨯=⨯∆=2a 233 晶胞的体积:V=332r 224a 23a 38a 233C S ==⨯=⨯ n=1232126112+⨯+⨯=6个 74.062r224r 346x 33≈π=π⨯= (5)对于金刚石结构,晶胞的体对角线BG=3r 8a r 24a 3=⇒⨯= n=8, Vc=a 334.063r 338r 348a r 348x 33333≈π=π⨯=π⨯=1.2、试证:六方密排堆积结构中633.1)38(a c 2/1≈= 证明:在六角密堆积结构中,第一层硬球A 、B 、O 的中心联线形成一个边长a=2r 的正三角形,第二层硬球N 位于球ABO 所围间隙的正上方并与这三个球相切,于是: NA=NB=NO=a=2R.即图中NABO 构成一个正四面体。
第21. 有一晶体,平衡时体积为 0V , 原子间相互作用势为0.如果相距为 r 的两原子互作用势为 ()n m r r a r u β+-= 证明(1) 体积弹性模量为 K=.90V mnU (2) 求出体心立方结构惰性分子的体积弹性模量.[解答]设晶体共含有 N 个原子,则总能量为U(r)=()∑∑i jij r u '21. 由于晶体表面层的原子数目与晶体内原子数目相比小得多,因此可忽略它们之间的基异,于是上式简化为 U=().2'∑jijr u N设最近邻原子间的距离为R 则有j ij a r =R再令 A ,1'∑=j m j m a A ,1'∑=jn j n a 得到 U=.200⎪⎪⎭⎫ ⎝⎛+-n n m m R A R A N βα 平衡时R=R 0,则由已知条件U(R 0) = 0U 得0002U R A R A N n n m m =⎪⎪⎭⎫⎝⎛+-βα 由平衡条件 0)(0=R dRR dU得021010=⎪⎪⎭⎫⎝⎛-++n nm m R A n R A m N βα. 由(1),(2)两式可解得.)(2,)(20000n n m m nR n m N U A nR n m N U A -=-=βα利用体积弹性模量公式[参见《固体物理教程》(2.14)式]K=0220209R R U V R ⎪⎪⎭⎫ ⎝⎛∂∂得K= ⎥⎦⎤⎢⎣⎡+++-n n m m R A n n R A m m N V 000)1()1(291βα = ⎥⎦⎤⎢⎣⎡-++-+-)(2)1()(2)1(2910000000n m N mR U R n n n m N nR U R m m N V nnm m = .900V mn U - 由于,00<U 因此,00U U -= 于是 K= .90V mnU (1) 由《固体物理教程》(2.18)式可知,一对惰性气体分子的互作用能为.)(126r B r A r u +-=若令 61,42⎪⎭⎫⎝⎛==A B B A σε,则N 个惰性气体分子的互作用势能可表示为⎥⎥⎦⎤⎢⎢⎣⎡⎪⎭⎫ ⎝⎛-⎪⎭⎫ ⎝⎛=6612122)(R A R A N r U σσε.由平衡条件0)(0=R dRR dU 可得 R .2616120⎪⎪⎭⎫ ⎝⎛=A A σ进一步得 .2)(122600A A N R U U ε-==代入K=.900V mn U 并取 m =6,n =12,V 300334R N =得 K=5126123233⎪⎪⎭⎫⎝⎛A A A σε.对体心立方晶体有 A .11.9,25.12126==A 于是.1.703σε=K 2. 一维原子链,正负离子间距为a ,试证:马德隆常数为2=μ1n2. [解答] 相距ij r 的两个离子间的互作用势能可表示成.4)(2n ijij ij r br q r u +=π设最近邻原子间的距离为R 则有 R a r j ij =, 则总的离子间的互作用势能 U=()∑∑∑-⎪⎪⎭⎫ ⎝⎛±-=jn jn j j j ij a bRa R q N r u N ''0'114[22πε. 基中 jja 1'±=∑μ 为离子晶格的马德隆常数,式中+;- 号分别对应于与参考离子相异和相同的离子.任选一正离子作为参考离子,在求和中对负离子到正号,对正离子取负号,考虑到对一维离子两边的离子是正负对称分布的,则有.413121112)1('⎥⎦⎤⎢⎣⎡+-+-=±=∑ j ja μ利用正面的展开式 1n(1+x ),432432 +-+-x x x x 并令 1=x 得 +-+-41312111=1n(1+1)=1n2.于是,一维离子链的马德常数为2=μ1n23. 计算面心立方面简单格子的6A 和12A(1) 只计最近邻; (2) 计算到次近邻; (3) 计算到次近邻.[解答]图2.26示出了面心立方简单格子的一个晶胞.角顶O 原子周围有8个这样的晶胞,标号为1的原子是原子O 的最近邻标号为2的原子是O 原子的最近邻,标号为3的原子是O 原子的次次近邻.由此得到,面心立方简单格子任一原子有12个最近邻,6个次近邻及24个次次近邻.以最近邻距离度量,其距离分别为:.3,2,1===j j j a a a 由 .1,112'126'6⎪⎪⎭⎫⎝⎛=⎪⎪⎭⎫ ⎝⎛=∑∑jj j j a A a A图2.6 面心立方晶胞得(1) 只计最近邻时1211*12)1(66=⎪⎭⎫⎝⎛=A , 1211*12)1(1212\=⎪⎭⎫⎝⎛=A .(2) 计算到次近邻时.094.1221*611*12)2(,750.1221*611*12)2(121212666=⎪⎭⎫ ⎝⎛+⎪⎭⎫⎝⎛==⎪⎭⎫ ⎝⎛+⎪⎭⎫⎝⎛=A A(3) 计算到次次近邻时.127.12033.0094.1231*2421*611*12)3(,639.13899.0750.1231*2421*611*12)3(121212126666=+=⎪⎪⎭⎫ ⎝⎛+⎪⎭⎫ ⎝⎛+⎪⎭⎫⎝⎛==+=⎪⎪⎭⎫ ⎝⎛+⎪⎭⎫ ⎝⎛+⎪⎭⎫⎝⎛=A A 由以上可以看出,由于12A中的幂指数较大,12A 收敛得很快,而6A 中的幂指数较小,因此 6A 收敛得较慢,通常所采用的面心立方简单格子的 6A 和 12A 的数值分别是14.45与12.13.4. 用埃夫琴方法计算二维正方离子(正负两种)格子的马德隆常数. [解答]马德隆常数的定义式为 jja 1'±=∑μ,式中+、-号分别对应于与参考离子相异和相同的离子,二维正方离子(正负两种)格子,实际是一个面心正方格子,图 2.7示出了一个埃夫琴晶胞.设参考离子O 为正离子,位于边棱中点的离子为负离子,它们对晶胞的贡献为4*(1/2).对参考离子库仑能的贡献为图2.7二维正方离子晶格.121*4顶角上的离子为正离子,它们对晶胞的贡献为4*(1/4), 对参考离子库仑能的贡献为 .241*4-因此通过一个埃夫琴晶胞算出的马德隆常数为 .293.1241*4121*4=-=ν再选取422=个埃夫琴晶胞作为考虑对象,这时离子O 的最的邻,次近邻均在所考虑的范围内,它们对库仑能的贡献为,2414-而边棱上的离子对库仑能的贡献为 ,521*8221*4+- 顶角上的离子对为库仑能的贡献为 ,841*4-这时算出的马德隆常数为图 2.8 4个埃夫琴晶胞同理对932=个埃夫琴晶胞进行计算,所得结果为611.11841*41321*81021*8321*48458242414=⎪⎪⎪⎪⎭⎫ ⎝⎛-+-+⎪⎪⎭⎫ ⎝⎛-+-+⎪⎭⎫ ⎝⎛-=μ 对 1642=个埃夫琴晶胞进行计算,所得结果为614.13241*42521*81721*81021*8421*4184138108348458242414=⎪⎪⎪⎪⎭⎫ ⎝⎛-+-+-+⎪⎪⎭⎫ ⎝⎛-+-+⎪⎪⎭⎫ ⎝⎛-+-+⎪⎭⎫ ⎝⎛-=μ当选取 n 2个埃夫琴晶胞来计算二维正方离子(正负两种)格子的马德隆常数,其计算公式(参见刘策军,二维NaC1 晶体马德隆常数计算,《大学物理》,Vo1.14,No.12,1995.)为 [][].1,8411>+++=--n D C B A n n n n μ其中 ,21)1(,1)1(11111nB t A n n n t t n +-=+--=-=∑,1)1(1)1()2()1(1)1()1(2112212221112122122222222221⎪⎪⎪⎪⎪⎭⎫⎝⎛+--+--+-+-+--++⎪⎪⎭⎫⎝⎛+++-+⎪⎪⎭⎫ ⎝⎛+---n n n n n C n n.121)1()1(2181222222+-+-++++-=n n n n n D n n5. 用埃夫琴方法计算CsCl 型离子晶体的马德隆常数(1) 只计最近邻 (2) 取八个晶胞 [解答](1) 图2.29是CsCl 晶胸结构,即只计及最近邻的最小埃夫琴晶胞,图2.29()a 是将Cs +双在体心位置的结构,图2.9(a)是将 Cl -取在体心位置的结构,容易求得在只计及最近邻情况下,马德隆常数为1.图2.29 (a )Cs 取为体心的CsC1晶胞,(b) C1取为体心的CsC1晶胞(2)图2.10是由8个CsCl 晶胞构成的埃夫琴晶胞,8个最近邻在埃夫琴晶胞内,每个离子对晶胞的贡献为1,它们与参考离子异号,所以这8个离子对马德隆常数的贡献为8埃夫琴晶胞6个面上的离子与参考离子同号,它们对埃夫琴晶胞的贡献是21,它们与参考离子的距离为32R 它们对马德隆常数的贡献为-()3/2*621图 2.10 8个CsCl 晶胞构成的一个埃夫琴晶胞埃夫琴晶胞楞上的12个离子,与参考离子同号,它们对埃夫琴晶胞的贡献是41它们与参考离子的距离为322R 它们对马德隆常数的贡献为-()3224/1*12埃夫琴晶胞角顶上的 8个离子,与参考离子同号,它们对埃夫琴晶胞的贡献是81它们与参考离子的距离为2R 它们对马德隆常数的贡献为 -()281*8,由8个CsCl 晶胞构成的埃夫琴晶胞计算的马德隆常数.064806.32)8/1(*8322)4/1(*123/2)2/1(*68=---=μ 为了进一步找到马德常数的规律,我们以计算了由27个CsCl 晶胞构成的埃夫琴晶胞的马德隆常数,结果发现,由27个CsCl 晶胞构成的埃夫琴晶胞的马德隆常数是0.439665.马德隆常数的不收敛,说明CsCl 晶胞的结构的马德隆常数不能用传统的埃夫琴方法计算.为了找出合理的计算方法,必须首先找出采用单个埃夫琴晶胞时马德隆常数不收敛的原因.为了便于计算,通常取参考离子处于埃夫琴晶胞的中心.如果以Cs +作参考离子,由于埃夫琴晶胞是电中性的要求,则边长为pa 2(p 是大于或等于1的整数)的埃夫琴晶胞是由(2p )3个CsCl 晶胞所构成,埃夫琴晶胞最外层的离子与参考离子同号,而边长为(2p +1)的埃夫琴晶胞是由(2p +1)3 个 CsCl 晶胞所构成,但埃夫琴晶胞的最外层离子与参考离子异号,如果以C1-作参考离子也有同样的规律,设参考离子处于坐标原点O ,沿与晶胞垂直的方向(分别取为x,y,z 图2.11示出了z 轴)看去,与参考郭同号的离子都分布在距O 点ia 的层面上,其中i 是大于等于 1的整数,与 O 点离子异号的离子都分布在距O 点(i -0.5)a 的层面上,图 2.11(a) 示出了同号离子层,图2.11(b)示出了异号离子层.图2.11 离子层示意图(a)表示同号离子层, O 离子所在层与 O '离子所在层相距ia(b)表示异号离子层, O 离子所在层和O ' 离子所在层相距(i -0.5)a当 CsCl 埃夫琴晶胞边长很大时,晶胞最外层的任一个离子对参考离子的库仑能都变得很小,但它们对参考离子总的库仑能不能忽略.对于由(2p )3个CsCl 晶胞所构成的埃夫琴晶胞来说,最外层有6*(2p )2个与参考离子同号的离子,它们与参考离子的距离为(1/2)pa ~(23)pa ,它们与参考离子的库仑能为a pe 024πε量级,这是一个相对大的正值.对于由(2p +1)3个CsCl 晶胞所构成的埃夫琴晶胞来说,离外层有6*(2p +1)2个与参考离子异号的离子,它们与参考离子的库仑能为a pe 024πε-量级,这是一个绝对值相对大的负值,因此,由(2p )3个CsCl 晶胞构成的埃夫琴晶胞所计算的库仑能,与由(2p +1)3个CsCl 晶胞构成的埃夫琴晶胞所计算的库仑能会有较大的差异.即每一情况计算的库仑能都不能代表CsCl 晶体离子间相互作用的库仑能.因此这两种情况所计算的马德隆常数也必定有较大的差异,由1个CsCl 晶胞、8个CsCl 晶胞和27个CsCl 晶胞构成的埃夫琴晶胞的计算可知, CsCl 埃夫琴晶胞体积不大时,这种现象已经存在.为了克服埃夫琴方法在计算马德隆常数时的局限性,可采取以下方法,令由 (2p )3个CsCl 晶胞构成的埃夫琴晶胞计算的库仑能为1U ,由(2p +1)3个CsCl 晶胞构成的埃夫琴晶胞所计算的库仑能为1U ,则CsCl 晶体离子间相互作用的库仑能可近似取作 )(2121U U U +=(1) 因子1/2 的引入是考虑除了(2p +1)3个CsCl 晶胞构成的埃夫琴晶胞最外层离子外,其他离子间的库仑能都累计了两偏,计算1U 和2U 时要选取体积足够大的埃夫琴晶胞,此时埃夫琴晶胞最外层离子数与晶胞内的离子数相比是个很小的数,相应的马德隆常数应为 )(2121μμμ+=(2) 其中:=1μ⎪⎪⎭⎫ ⎝⎛±∑i ja 1'是由(2p )3个CsC1晶胞构成的埃夫琴晶胞计算的值; =1μ⎪⎪⎭⎫ ⎝⎛±∑i ja 1'由 (2p +1)3 个CsC1晶胞构成的埃夫琴晶胞所计算成本的值.为简化计算,特选取晶胞边长a 为计算单位,由于,32a R =所以,23'μμ= ⎪⎪⎭⎫ ⎝⎛±=∑'''1i i a μ (3) 其中'i a 是某一离子到参点的距离与a 的比值.考虑到对称性,对选定的埃夫琴晶胞,把晶胞的离子看成分布在一个个以参考离子为对称心的正六面体的六个面上,体积不同的正六面六个面上的离子分别计算.由(2p )3个CsC1晶胞构成埃夫琴晶胞时,由分析整理可得,231111⎪⎪⎭⎫ ⎝⎛++=∑∑=-=p pi i p i i C B A μ (4) 由(2p +1)3个 CsC1 晶胸构成埃夫琴晶胞时,,231112⎪⎪⎭⎫ ⎝⎛++=∑∑=-=p pi i p i i D B A μ (5)其中:),1(''''22'2'p i i y x k A i x iy y x i <≤++-=∑∑(6)i A 表示与 O 点距离为ia 的6个面上所有的离子对马德隆常数的面贡献,因为这些离子与参考离子同号,故到负号.'x 、'y 是离子在平面 '''y x o 上的坐标, ''y x k 代表 6个面上等价离子的个数,其取值规则为:(1) 在角上(如E 点),即'x =i 且 'y = i. 时, ''y x k =8;(2) 在棱与坐标轴的交点(如 F 点),'x =i 且'y = 0或 'x =0且'y = 0时, ''y x k =6 (3) 在棱上的其他点(如H 、I 点)即不满足上述条件,且'x =i 或'y = i.时, ''y x k =12 (4) 在'O 点,即'x =0且'y = 0时, ''y x k =6(5) 在除'O 点外的面上的点(如J 点),即不满足上述条件时,''y x k =24.),1()5.0(5.05.05.05.022'2''''''p i i y x k B i x i y yx i ≤≤-++=∑∑-=-=(7)i B 代表距O 点距离为(i -0.5)a 的6个面上的离子对马德隆常数的贡献,因为这种些离子与参考离子异号,故取正号. 'x ,'y 是离子在平面'''y x o 上的坐标, '''y x k 代表这6个面上等价离子的个数,其取值规则为:(1) 在角上(如K 点),即'x =i 且 'y = i.时, '''y x k =8;(2) 在棱下(如L 、M 点),即不满足不述条件,且'x =i 或'y = i 时,'''y x k =12; (3) 在面上(如N 点)好不满足上述条件时, '''y x k =24.),(0022'2'"''''p i i y x k C i x iy i yx =++-=∑∑==i C 表示在边长为2pa 的晶胞最外层,即与参考离子相距pa 的6个面上的离子对马德隆常数的贡献,应取负号,与iA 的不同在于"''y x k的取值: (1) 在角上, "''y x k =''y x k /8; (2) 在棱上, "''y x k =''y x k /4; (3) 在面上, "''y x k=''y x k /2.),()5.0(5.05.05.05.022'2''''''''p i i y x k D i x i y yx i =-++=∑∑-=-=i D 表示在边长为2a p )1(+的晶胞最外层,即与参考离子相距(p +0.5)a 的离子层对马德隆常数的贡献,应取正号,与i B 的不同在于'''''yx k 的取值: (1) 在角上, '''''y x k ='''y x k /8; (2) 在棱上, '''''y x k ='''y x k /4; (3) 在面上, '''''y x k ='''y x k /2.表2.1给出了计算结果,给出的μ是由分别对应2p 和2p+1的1μ和2μ求得的,实际上, 1μ和2μ只需对应边长相近的埃夫琴晶胞即可,如取对应2p 和2p-1的埃夫琴晶胞也可得到一样的收敛结果,由以上数据可见,马德隆常数μ随晶胞边长的增大而迅速收敛.该方法适用于NaC1结构以外离子晶体马德隆常数的计算.6.只计及最近邻间的排斥作用时,一离子晶体离子间的互作用势为⎪⎪⎩⎪⎪⎨⎧±-=-)2(,)1(,)(22r e R e e r u R ρλ(1)最近邻(2)最近邻以外 式中ρλ,是常数,R 是最近邻距离,求晶体平衡时,原子间总的互作用势.[解 答]设离子数目为2N,以j ij a r =R 表示第j 个离子到参考离子i 的距离,忽略表面效应,则总的相互作用能可表示为U =N ⎥⎥⎦⎤⎢⎢⎣⎡+⎪⎪⎭⎫ ⎝⎛±-∑∑-ρλR j j e R a e 2' (∑表示最近邻)=N ,2⎥⎦⎤⎢⎣⎡+--ρλμR e Z R e其中⎪⎪⎭⎫⎝⎛±=∑j ia 1'μ 为马德隆常数,+号对应于异号离子,-号对应于同号离子;Z 为任一离子的最近邻数目,设平衡时R=R 0 ,由平衡条件,02020=⎥⎦⎤⎢⎣⎡+=-ρρλμR R e Z R e N dRdU 得.0202ρλμρR e Z R e -=平衡时的总相互作用为.1)(0020200⎪⎪⎭⎫⎝⎛-=⎥⎦⎤⎢⎣⎡+-=-R R e N e Z R e N R U R ρμλμρ 7. 设离子晶体中,离子间的互作用势为⎪⎪⎩⎪⎪⎨⎧±+-=最近邻以外最近邻,,)(22re R b R e r u m(1) 求晶体平衡时,离子间总的相互作用势能)(0R U (2) 证明: )(0R U 11-⎪⎪⎭⎫⎝⎛∝m mZ μ其中μ是马德隆常数,Z 是晶体配位数 [解答](1)设离子数目为2N , 以j ij a r =R 表示第j 个离子到参考离子i 的距离,忽略表面效应,则总的相互作用能可表示U =N ⎥⎥⎦⎤⎢⎢⎣⎡+⎪⎪⎭⎫ ⎝⎛±-∑∑m j j R b R a e 2'(∑表示最近邻) =N ,2⎥⎦⎤⎢⎣⎡+-m R b Z Re μ其中⎪⎪⎭⎫ ⎝⎛±=∑j i a 1'μ,为马德隆常数,+号对应于异号离子,-号对应于同号离子.Z 为任一离子的最近邻数目,设平衡时R=R 0由平衡条件,0102020=⎥⎦⎤⎢⎣⎡-=+m R R Zmb R e N drdUμ得10-m R Zmb=2e μ即1120-⎪⎪⎭⎫ ⎝⎛=m e Zmb R μ.于是,晶体平衡时离子间总的相互作用势能0U =).1(000--=⎥⎦⎤⎢⎣⎡+-m R NZbR b Z R Zmb N m m m(2)晶体平衡时离子间的相互作用势能可进一步化为0U =.)()()1()1(1111121211--------=⎪⎪⎭⎫ ⎝⎛--m m m m mm m m m m mb Ze Nbm e Zmb ZNbm μμ由上式可知 .110-⎪⎪⎭⎫⎝⎛∝m mZ U μ8.一维离子链,其上等间距载有正负2N 个离子,设离子间的泡利排斥只出现在最近邻离子之间,且为b/R n,b,n 是常R 是两最近邻离子的间距,设离子电荷为q ,(1) 试证明平衡间距下 )(0R U =;114212002⎪⎭⎫⎝⎛--n R n Nq πε(2) 令晶体被压缩,使0R )1(0δ-→R , 试证明在晶体被压缩单位长度的过程中外力作功的主项为c 2δ其中c=;21)1(02R n q n -(3) 求原子链被压缩了2)1(0<<e e NR δδ时的外力[解答](1) 因为离子间是等间距的,且都等于R ,所以认定离子与第j 个离子的距离j r 总可表示成为R a r j j =ja 是一整数,于是离子间总的互作用势能⎥⎥⎦⎤⎢⎢⎣⎡-⎪⎪⎭⎫ ⎝⎛±-=⎥⎥⎦⎤⎢⎢⎣⎡+=∑∑n i in j j j R b a R q N r b r q N R U 214242)('202'0πεπε,其中+、-分别对应相异离子和相同离子的相互作用.一维离子晶格的马德隆常数(参见本章习题2)为=⎪⎪⎭⎫ ⎝⎛±∑i ia 1'21n2. 利用平衡条件0)(0=R dRR dU得到b=nq 01-n 0241n2R πε,)(R U =⎪⎪⎭⎫ ⎝⎛---n n nR R R Nq 102141n22πε. 在平衡间距下⎪⎭⎫⎝⎛--n R Nq R U 1141n22)(0020πε.(2) 将互作用势能在平衡间距附近展成级数+-⎪⎪⎭⎫ ⎝⎛+-⎪⎭⎫⎝⎛+=202200)(21)()()(0R R dR U d R R dR dU R U R U R R 由外力作的功等于晶体内能的增量,可得外力作功的主项为W=20220)(21)()(0R R dR U d R U R U R-⎪⎪⎭⎫ ⎝⎛=-, 其中利用平衡条件,将R=R )1(0δ- ,代入上式,得到W=δδπε)2(421)1(2102002NR R n q n ⎥⎥⎦⎤⎢⎢⎣⎡-. 晶体被压缩单位长度的过程中,外力作的功的主项δ02W NR =δπε⎥⎥⎦⎤⎢⎢⎣⎡-2002421)1(21R n q n 令c=202421)1(R n q n πε-(CGS)得到在晶体被压缩单位长度的过程中,外力作的功的主项为2δc . (3)设e δδ=时外力为F e ,由于在弹性范围内,外力与晶格的形变成正比,所以 F= )2(0δαNR , F e = )2(0e NR δα,其中α为比例系数离子链被压缩e NR δ02过程中外力作的功W e =δδαδδd NR NR Fdx e eNR e 020002)]2([0⎰⎰== e e e F NR NR δδα022022121)2(=.由于 W e =)2(20e eNR c δδ,所以离子链被压缩了e NR δ02时的外力为F e =202)1(21R n n q c ee δδ-=.9.设泡利排斥项的形式不变,讨论电荷加倍对NaC1晶格常数,体积弹性模量以及结合能的影响。