第一章 晶体结构后期编辑：霍团长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）面心立方的正格子基矢（固体物理学原胞基矢）：123()2()2()2a a j k a a i k a a i j ⎧=+⎪⎪⎪=+⎨⎪⎪=+⎪⎩由倒格子基矢的定义：1232()b a a π=⨯Ω31230,,22(),0,224,,022a aaa a a a a a a Ω=⋅⨯==，223,,,0,()224,,022i j ka a a a a i j k aa ⨯==-++ 213422()()4ab i j k i j k a aππ∴=⨯⨯-++=-++同理可得：232()2()b i j k ab i j k aππ=-+=+-即面心立方的倒格子基矢与体心立方的正格基矢相同。

所以，面心立方的倒格子是体心立方。

（2）体心立方的正格子基矢（固体物理学原胞基矢）：123()2()2()2a a i j k a a i j k a a i j k ⎧=-++⎪⎪⎪=-+⎨⎪⎪=+-⎪⎩由倒格子基矢的定义：1232()b a a π=⨯Ω3123,,222(),,2222,,222a a a a a a a a a a a a a -Ω=⋅⨯=-=-，223,,,,()2222,,222i j k a a a a a a j k a a a ⨯=-=+-213222()()2a b j k j k a aππ∴=⨯⨯+=+同理可得：232()2()b i k a b i j aππ=+=+即体心立方的倒格子基矢与面心立方的正格基矢相同。

所以，体心立方的倒格子是面心立方。

1.5、证明倒格子矢量112233G h b h b h b =++垂直于密勒指数为123()h h h 的晶面系。

证明：因为33121323,a aa a CA CB h h h h =-=-，112233G h b h b h b =++ 利用2i j ij a b πδ⋅=，容易证明12312300h h h h h h G CA G CB ⋅=⋅=所以，倒格子矢量112233G h b h b h b =++垂直于密勒指数为123()h h h 的晶面系。

1.6、对于简单立方晶格，证明密勒指数为(,,)h k l 的晶面系，面间距d 满足：22222()d a h k l =++，其中a 为立方边长；并说明面指数简单的晶面，其面密度较大，容易解理。

解：简单立方晶格：123a a a ⊥⊥，123,,a ai a aj a ak === 由倒格子基矢的定义：2311232a a b a a a π⨯=⋅⨯，3121232a a b a a a π⨯=⋅⨯，1231232a a b a a a π⨯=⋅⨯倒格子基矢：123222,,b i b j b k a a aπππ=== 倒格子矢量：123G hb kb lb =++，222G h i k j l k a a aπππ=++晶面族()hkl 的面间距：2d Gπ=2221()()()h k l a a a=++22222()a d h k l =++ 面指数越简单的晶面，其晶面的间距越大，晶面上格点的密度越大，单位表面的能量越小，这样的晶面越容易解理。

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

解：(111) (111)1、(111)面与(100)面的交线的AB ，AB 平移，A 与O 点重合，B 点位矢：B R aj ak =-+， (111)面与(100)面的交线的晶向AB aj ak =-+，晶向指数[011]。

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

第二章 固体结合2.1、两种一价离子组成的一维晶格的马德隆常数（2ln 2=α）和库仑相互作用能，设离子的总数为2N 。

＜解＞ 设想一个由正负两种离子相间排列的无限长的离子键，取任一负离子作参考离子（这样马德隆常数中的正负号可以这样取，即遇正离子取正号，遇负离子取负号），用r 表示相邻离子间的距离，于是有(1)11112[ (234)ij rr r r r rα±'==-+-+∑ 前边的因子2是因为存在着两个相等距离i r 的离子，一个在参考离子左面，一个在其右面，故对一边求和后要乘2，马德隆常数为234(1) (34)n x x x x x x +=-+-+ 当X=1时，有1111 (2234)n-+-+= 2.3、若一晶体的相互作用能可以表示为 ()mnu r r r αβ=-+试求：（1）平衡间距0r ；（2）结合能W （单个原子的）； （3）体弹性模量；（4）若取02,10,3,4m n r A W eV ====，计算α及β的值。

解：（1）求平衡间距r 0由0)(0==r r drr du ，有：mn nm n m m n n m r r n r m --++⎪⎭⎫⎝⎛=⎪⎪⎭⎫ ⎝⎛=⇒=-1101.0100αββαβα结合能：设想把分散的原子（离子或分子）结合成为晶体，将有一定的能量释放出来，这个能量称为结合能（用w 表示）（2）求结合能w （单个原子的）题中标明单个原子是为了使问题简化，说明组成晶体的基本单元是单个原子，而非原子团、离子基团，或其它复杂的基元。

显然结合能就是平衡时，晶体的势能，即U min即：n m rrr U W 00)(βα-+=-= （可代入r 0值，也可不代入）（3）体弹性模量由体弹性模量公式：0220209r r U V r k ⎪⎪⎭⎫ ⎝⎛∂∂=1112[1...]234α=-+-+22n α∴=（4）m = 2，n = 10，A r 30=， w = 4eV ，求α、β818105210⎪⎭⎫ ⎝⎛=⎪⎭⎫ ⎝⎛=αβαβr ①)5(54)(802010.20代入αβαβα=-=+-=r r r rr UeV r r U W 454)(200==-=⇒α② 将A r 30=，J eV 1910602.11-⨯=代入①②211523810459.910209.7mN m N ⋅⨯=⋅⨯=⇒--βα 详解：（1）平衡间距r 0的计算 晶体内能()()2m n N U r r rαβ=-+ 平衡条件0r r dUdr==，11000m n m n r r αβ++-+=，10()n m n r m βα-= （2）单个原子的结合能01()2W u r =-，00()()m n r r u r r r αβ==-+，10()n m n r m βα-= 1(1)()2mn m m n W n m βαα--=-（3）体弹性模量0202()V UK V V ∂=⋅∂ 晶体的体积3V NAr =，A 为常数，N 为原胞数目 晶体内能()()2m n N U r r rαβ=-+ U U r V r V ∂∂∂=∂∂∂1121()23m n N m n r r NArαβ++=- 221121[()]23m n U N r m n V V r r r NArαβ++∂∂∂=-∂∂∂ 022222000001[]29m n m n V V U N m n m n V V r r r r αβαβ=∂=-+-+∂ 由平衡条件1120001()023m n V V U N m n Vr r NAr αβ++=∂=-=∂，得00m n m n r r αβ=222220001[]29m n V V U N m n V V r r αβ=∂=-+∂ 02220001[]29m nV V U N m n m n V V r r αβ=∂=-+∂2000[]29m n N nm V r r αβ=--+ 000()2m n N U r r αβ=-+ 020220()9V V Umn U V V =∂=-∂ 体弹性模量009mnK U V = （4）若取02,10,3,4m n r A W eV ====10()n mn r m βα-=，1(1)()2m n m m n W n m βαα--=- 1002W r β=，20100[2]r W r βα=+ -95101.210eV m β=⨯⋅，1929.010eV m α-=⨯⋅2.7、对于2H ，从气体的测量得到Lennard —Jones 参数为65010, 2.96.J A εσ-=⨯=计算fcc 结构的2H 的结合能[以KJ/mol 单位)，每个氢分子可当做球形来处理．结合能的实验值为0.751kJ ／mo1，试与计算值比较．＜解＞ 以2H 为基团，组成fcc 结构的晶体，如略去动能，分子间按Lennard —Jones 势相互作用，则晶体的总相互作用能为：1261262.ij ij i j U N P P R R σσε--⎡⎤⎛⎫⎛⎫''=-⎢⎥ ⎪ ⎪⎝⎭⎝⎭⎢⎥⎣⎦∑∑61214.45392;12.13188,ijij jiP P --''==∑∑16235010, 2.96, 6.02210/.erg A N mol εσ-=⨯==⨯()()12628162.96 2.962602210/501012.1314.45 2.55/.3.16 3.16U U mol erg KJ mol -⎡⎤⎛⎫⎛⎫=⨯⨯⨯⨯⨯-≈-⎢⎥ ⎪ ⎪⎝⎭⎝⎭⎢⎥⎣⎦0将R 代入得到平衡时的晶体总能量为。