国科大XRD作业参考答案Word文档下载推荐.docx
《国科大XRD作业参考答案Word文档下载推荐.docx》由会员分享，可在线阅读，更多相关《国科大XRD作业参考答案Word文档下载推荐.docx（17页珍藏版）》请在冰豆网上搜索。
SrVO3:
cubic,3.841Å
注2：
各原子或离子的散射因子拟合参数参见网络课堂上的上传国际表格，也可通过网络搜索得到；
目前O2－的参数只能用O－的参数代替（x）。
解一：
取用Ca1-xSrxTiO3固溶体
钙钛矿具有ABO3结构，Sr和Ca占据（0,0,0）；
O2-占据（0,1/2,1/2）,（1/2,0,1/2）,（1/2,1/2,0）；
Ti4+占据（1/2,1/2,1/2）。
要让（hkl）=（001）峰消光，即让I∝|F001|2=0
首先计算其结构散射因子为：
（1）
其中:
=2
（2）
（001）方向即h=k=0，l=1；
将
（2）式具体展开可得:
F001=
=
（3）
要使（001）峰消失只要令（3）式为0即可，为求解该方程，可先利用直通光情况（Q=0）的散射因子（等于核外电子数）作近似计算：
36x+（1-x）18-18-10=0
x=5/9=0.5556
CaTiO3、SrTiO3都是立方钙钛矿结构，晶格常数相近，且容忍因子也都接近于1，所以能够互相溶解且保持立方相结构，根据Vegard定律可得到其近似晶格常数:
c=（1-x）×
3.827+x×
3.905（4）
代入x=5/9,可得c=3.8703Å
=0.0166898（5）
原子结构因子的计算公式为:
（6）
通过查阅InternationalTableofCrystallography可获得
利用式（6）和式（5）的结果，编程求解各元素的原子结构因子可得：
=33.13058,
=16.28962,
=7.23750,
=16.78898
再将该结果回代到公式（3）中可解得x=0.4594。
利用该x值重新计算固熔体的晶格常数c和
C=3.8628，
=0.01669264
重复上述步骤，得到x=0.4593,基本收敛。
所以材料配比为
解二：
取用Sr1-xCaxVO3固溶体
求解步骤与上面类似：
晶胞中各粒子的坐标：
Sr2+（Ca2+）:
（0,0,0）
V5+:
（½
½
）
O2-:
0）,（½
0,½
）,（0,½
a=3.8135Å
（Average）（or3.806Å
accordingtoorthorhombicCaVO3）
q=2/a=1.6476Å
-1
（q/4）2=0.0171906
A（q）=
[（1-x）fSr+xfCa+fVexp（iq（1/2,1/2,1/2））+fO[exp（iq（0,1/2,½
））+exp（iq（1/2,0,½
））+exp（iq（1/2,1/2,0））]
[（1-x）fSr+xfCa-fV-fO]
（1-x）fSr+xfCa-fV-fO=0
（1-x）*{18.0874*exp（-1.4907*0.0171906）+8.1373*exp（-12.6963*0.0171906）+2.5654*exp（-24.5651*0.0171906）-34.193*exp（0.0138*0.0171906）+41.4025}+x*{15.6348*exp（0.0074*0.0171906）+7.9518*exp（-0.6089*0.0171906）+8.4372*exp（-10.3116*0.0171906）+0.8537*exp（-25.9905*0.0171906）-14.875}-{15.6887*exp（-0.679003*0.0171906）+8.14208*exp（-5.40135*0.0171906）+2.03081*exp（-9.97278*0.0171906）-9.576*exp（0.940464*0.0171906）+1.7143}-{4.1916*exp（-12.8573*0.0171906）+1.63969*exp（-4.17236*0.0171906）+1.52673*exp（-47.0179*0.0171906）-20.307*exp（0.01404*0.0171906）+21.9412}=0
解方程得x=0.55
用此x值将晶格常数修正为
0.45*3.841+0.55*3.786=3.8108
所以
（q/4）2=0.0172
重复上述步骤得到新的x值
x=0.55
该值已经收敛，所以最后的化合物为Sr0.45Ca0.55VO3（如果不考虑V3+成分）。
2.推导出双分子气体Br2气的散射强度I与散射矢量Q的关系式，写出详细的论证与解决步骤，并使用任何你喜欢的计算机语言编写程序，画出Br2气的I－Q曲线。
改变Br2分子的结构（此即化学键长），I－Q曲线有何变化？
每个Br2分子内两个原子之间的散射相干，要振幅叠加；
但不同Br2分子之间无固定相位差，非相干，要强度叠加。
结构因子：
其中
为原子散射因子。
散射强度为：
O2中氧原子间距|r|固定，但不同氧分子的r的取向不同，并且随时间变化，因此应对r的方向进行平均才得到测量强度
r的取向呈球分布
所以，测量强度与Q的关系为:
图2-1I-Q曲线
若再考虑极化因子则：
，单位为
，所以
，设λ=1.54Å
（CuKα辐射）
程序为：
r=1.52;
Q=0.01:
0.04:
6;
%unit:
2*pi/r
c=（Q/（2*r））.^2;
f=3.0485*exp（-13.2771*c）+2.2868*exp（-5.7011*c）+1.5463*exp（-0.3239*c）+0.867*exp（-32.9089*c）+0.2508;
Intensity=2*f.^2.*（1+sin（Q*r）./（Q*r））;
%Unit:
r02
ssita=Q./（2*r）*1.54;
%sin（sita）
ssita_2=sqrt（1-ssita.^2）.*2.*ssita;
%sin（2sita）
fa=（1+（1-2.*ssita.^2）.^2）/2./ssita_2;
%计算洛仑兹因子和极化因子
Intensity=Intensity.*fa;
%semilogy（Q,Intensity）
plot（Q,Intensity）
xlabel（'
Q（2\pi/r）'
'
FontWeight'
'
Bold'
FontSize'
15）
ylabel（'
Intensity'
改变Br2分子中的原子间距，可以看出散射峰之间的距离随原子间距的增大而减小，这符合正空间的材料结构与其散射花样之间是倒易关系（傅立叶变换）。
3.EuTiO3具有理想的立方钙钛矿结构，其晶格常数与SrTiO3的相等，均为3.905Å
，这二者可以组成一个比较理想无应力体系，用于研究复杂结构的氧化物薄膜无应力生长的机理。
a）从课程网站上下载实验数据，使用你喜欢的任何软件画出EuTiO3粉末的XRD谱，然后计算并标出各个衍射峰的指数（写出计算过程和结果）。
b）*（选做）写一个电脑程序计算EuTiO3的粉末衍射谱，5<
2<
80,=1.54Å
（CuK辐射）。
可以用任何你喜欢的语言。
解决方案中应包括程序代码，两列数据（角度和强度）和相关图线，其中一张图谱是实验数据和计算数据的比较。
计算程序需要考虑的因素包括:
1）原子形状因子
2）结构因子
3）Debye-Waller因子
4）Compton散射
5）热漫散射
6）有限尺寸因子
7）多重性因子
8）仪器分辨率
关于仪器分辨率，X主光束的发散角为0.12°
，探测器的前狭缝处在距样品17cm远处，缝宽为0.1mm,因此张角为0.03°
根据实验测量结果，二者的卷积大约为0.13°
。
（1）衍射峰的位置由布拉格方程决定，其衍射峰的强度由以上因素决定：
立方结构的面间距：
c为晶格常数。
（100）方向，即h=1,k=l=0。
=11.372deg。
其它衍射峰的位置同理可得，结果列表如下：
入射角
（deg）
11.372
16.192
19.970
23.226
26.162
28.881
33.898
36.267
指标化（HKL）
100
110
111
200
210
211
220
300
%布拉格方程计算衍射峰的位置,其中第7项不存在,因为h^2+k^+l^!
=7
lamda=1.54;
%入射x光波长
c=3.905;
%晶格常数
fori=1:
1:
10
d=2*c/sqrt（i）;
sita（i）=asin（lamda/d）;
end
sita/pi*180%转化为角度
接下来以（100）方向为例计算衍射强度：
粉末衍射积累强度的表达式为：
各参数的意义为：
I：
衍射的强度；
I0：
入射X射线的强度；
R0：
样品到衍射环的距离；
re=e2/4pe0mc2，电子的经典半径；
N0：
单位体积样品中晶胞的数目；
：
X射线的波长；
Vcr：
样品被X射线照射到的面积；
F：
为结构因子；
Mf：
多重性因子；
为温度因子即Debye-Waller因子。
（1+cos22qB）/sin2qBcosqB：
角因子，也称为洛仑兹-偏振因子。
A为吸收因子。
现已知晶格常数，和入射X射线波长，再假定样品的体积一定，则可合并一些常数则累计强度为：
A
又因为A和
都受
变化的影响，且两者变化方向相反，虽然数值不完全相等，但大致可抵消[1]，故用相对强度表示为：
1.考虑原子形状因子：
（6）
=11.372，
=1.54
代入上式可得：
f_Eu=56.6490f_O=7.2621f_Ti=16.8091
2.结构因子：
其中:
（100）方向即h=k=0，l=1；
|F100|2=1061.3。
3.多重性因子
（100）方向的多重性因子是6。
这里给出（100），（110），（111），（200），（210）方向的衍射峰。
图3-1衍射峰的衍射谱，从左到右分别为（100），（110），（111），（200），（210）
源程序如下：
clear;
si=[11.37216.19219.97023.22626.162];
%各衍射峰的出现位置,这里计算前5个
%入射X线波长
sita=si./180.*pi;
%角度转化为弧度
Q=sin（sita）./lamda;
%计算国际晶体学表中的
%计算原子结构因子
O_1=[4.1916012.85731.639694.172361.5267347.0179-20.307-0.0140421.9412];
%氧离子-1
Ti_4=[19.51140.178848.234736.670182.01341-0.292631.5208012.9464-13.280];
%正四价钛离子
Eu_2=[24.00632.2778319.95040.17353011.803411.60963.8724326.51561.36389];
a=zeros（1,4）;
b=zeros（1,4）;
c=0;
4
a（i）=Eu_2（2*i-1）;
b（i）=Eu_2（2*i）;
c=Eu_2（2*i+1）;
f_Eu=form_factor（a,b,c,Q）
a（i）=Ti_4（2*i-1）;
b（i）=Ti_4（2*i）;
c=Ti_4（2*i+1）;
f_Ti=form_factor（a,b,c,Q）
a（i）=O_1（2*i-1）;
b（i）=O_1（2*i）;
c=O_1（2*i+1）;
f_O=form_factor（a,b,c,Q）
%原子结构因子
fstr=zeros（1,5）;
jiaodu=（1+cos（2.*sita）.^2）./sin（sita）./sin（sita）./cos（sita）;
%角因子
fstr
（1）=（f_Eu
（1）-f_Ti
（1）-f_O
（1））.^2;
%100面的晶胞结构因子
fstr
（2）=（f_Eu
（2）+f_Ti
（2）-f_O
（2））.^2;
%110
fstr（3）=（f_Eu（3）-f_Ti（3）+3*f_O（3））.^2;
%111
fstr（4）=（f_Eu（4）+f_Ti（4）+3*f_O（4））.^2;
%200
fstr（5）=（f_Eu（5）-f_Ti（5）-f_O（5））.^2;
%210
Mf=[6,12,8,6,24];
%多重性因子
In=fstr.*Mf.*jiaodu;
%最后的衍射峰的相对强度
stem（2*si,In）;
2\theta/（deg）'
FontName'
Times'
15,'
）;
Intensify'
axis（[10,55,0,12e5]）;
4.判断CsCl结构的X射线衍射谱中，（100）和（200）衍射谱线的强度那个大，为什么？
【参考答案】
如果是单晶体，则可根据单晶体反射强度公式计算（100）和（200）衍射峰的强度比值
是大于1还是小于1来判断这两个衍射峰强度的大小。
其中（100）和（200）的多重性因子相等。
查阅CsCl的晶格常数，假设所用X射线是CuK辐射，则可计算上面单晶和多晶的I100和I200的比值。
经过计算可知，CsCl的（200）峰强度要高于（100）峰的强度。
注意：
有的固体物理学教材上的解释－高指数晶面的原子面密度小（如方俊鑫、陆栋的1980版《固体物理学》上册P26：
“又由于面上的原子密度大，对射线的散射强，因而密勒指数简单的晶面族，在X射线衍射中，往往为照片中的浓黑斑点所对应”）是错误的。
5.为什么实际测到的衍射峰都来自低指数晶面的衍射？
首先，无论是衍射仪还是面探测器，其探测的最大衍射角2是180，根据布拉格公式
可知，
，即实际探测到的衍射面的晶面间距不能小于波长的一半。
因此，高指数hkl的晶面不能满足布拉格方程，也就是说，实际测量到的衍射峰都来自低指数晶面。
第二个方面的原因是，高指数晶面的衍射峰出现在高衍射角位置，但原子的散射因子随角度的增加而快速衰减，所以受此影响高角处的衍射峰强度很小，大于一定的角度衍射峰强度就会衰减到难以测量的程度，所以高指数晶面的衍射不出现。
其他方面的原因还有Lorentz因子，德拜温度因子也对高角处的衍射峰有限制。
6.（选作）阅读文章[QunShen,StefanKycia,Determinationofinterfacialstraindistributioninquantum-wirestructuresbysynchrotronx-rayscattering,Phys.Rev.B55,15791（1997）]，
从中体会对特殊的问题如何应用X射线衍射的运动学理论，并根据你所学运动学理论方法，写一篇对该论文的评论（Review）文章，要求讲明：
1）所要解决的问题；
2）所采用的方法的原理，必要时将论文中省去的强度表达式的推导写出来；
3）所获得的结果及其意义；
4）对解决方法的评论。
第六章衍射实验方法
问题1：
我现在在挪威奥斯陆大学交流学习，7月初就回去了。
他们新买了一台衍射仪，但是还没有调试好。
用这台仪器测试ZnO/c-sapphire样品，分别在31.14°
34.52°
37.56°
41.73°
存在很明显的峰位，其中第二个和第四个峰位应该是ZnO和蓝宝石，但是第一个和第三个不知道来自哪儿，会不会是和光源的单色性有关？
可以看出，所用的靶是铜靶，可以推算出，第一个和第三个峰是来自于Cu的Kβ辐射，衍射峰仍然是ZnO（002）和Al2O3（006）。
这意味着狭缝卡光源的位置不准，没有靶Kβ卡去。
下面是我计算的各个峰位和波长的对应关系：
实验：
31.14°
34.52°
37.56°
41.73°
CuKα1：
Al2O3（005）34.4926
Al2O3（006）41.682
ZnO（002）34.4775
CuKα2：
Al2O3（005）41.5808
Al2O3（006）41.7902
ZnO（002）34.5055
CuKα平均
Al2O3（005）34.5224
Al2O3（006）41.7186
CuKβ
Al2O3（005）31.1
Al2O3（006）37.52
ZnO（002）31.02
有几个问题，还想向王老师请教一下：
（1）如果探测器探测到的直通光不是高斯分布，这对测试结果影响有多大呢？
直通光不是高斯分布的原因，我想到的有：
a，X射线源光斑本身存在强度分布，由于狭缝原因，卡到的光斑不是高斯对称分布。
这种情况，只要光源单色性好，应该是不影响测试结果的。
b，由于探测器之前的狭缝影响，这种原因肯定会影响探测结果。
只要薄膜和衬底的峰不重合，则不影响分峰，但影响峰位的精细拟合；
如果要根据峰形变化分析缺陷的类型，则有致命的影响，不能使用。
但是固定探测器做rockingcurve，峰形比较对称，这是不是就可以说明探测器的狭缝没有问题？
不能说明。
（2）我们在测试theta-2thate时，总是要先通过固定探测器进行omega扫描和固定theta进行探测器扫描，寻找到最优化的theta及2theta值，这样做的目的是为了消除晶面斜切的影响。
对于同步辐射站上，我们可以根据优化后的结果重新定义theta及2theta值；
但是对于衍射仪，这个该怎么定义？
2theta都不能重新定义，theta定义方式相同，也可以不定义。
（3）所有角度的定义是以什么为参照呢？
如果这个参照不对，是不是会影响衍射峰的位置（也就是测试结果与实际结果存在偏差）？
2theta以直通光为参照，theta开始时以样品表面为参考，得到衍射角2theta之后，可以以Bragg角为参考，也可以仍然以表面为参考。
问题2：
、在劳厄法中，如果采用X射线管或X光机作光源，为什么降低光源的加速电压会导致透射式劳厄衍射图样中靠近中心的衍射点消失，但不会导致反射式劳厄衍射图样中的衍射点消失？
参考：
由E=h=hc/min=eV，得出min=hc/（eV），所以，如果加速电压V降低，则X光的最短波长min就会增大，从而满足劳厄条件的衍射点Q＝Ghkl＝4sin>
4minsin，因而低角度处的－即靠近中心的－衍射点会消失。
问题3、在粉末的X射线衍射中，具有同样面间距的衍射晶面是否具有同样的系统消光条件？
为什么？
参考答案：
是。
因为（写出晶胞结构因子表达式，指出，晶面间距相等，则Q＝Ghkl相等，各原子的位置rj不变，所以晶胞结构因子相等，因此计算结果相等，所以如果其一系统消逛，则另外一个也系统消光，即系统消光条件相同。
可以面心立方晶系的（300
升级会员 