用SHELXTL程序进行晶体结构分析的方法 单晶结构分析电子教案.ppt
《用SHELXTL程序进行晶体结构分析的方法 单晶结构分析电子教案.ppt》由会员分享,可在线阅读,更多相关《用SHELXTL程序进行晶体结构分析的方法 单晶结构分析电子教案.ppt(98页珍藏版)》请在三一办公上搜索。
1、第五章,用SHELXTL程序进行结构分析的方法,new open copy archive deleteexit,edit edit.res edit.ins edit.cif edit.prp edit.lst copy.res.ins,XlXH,一、SHELXTL文件,1.文件名,一般,同一结构,所有文件都用相同的名(不能超过个字符),只是扩展名不同,2.两个必要文件(由XPREP程序产生),*.hkl文件:所有的衍射点,每一点一行。格式为:h k l F2(F2),-4 2-2 21.189 3.050 4-2 2 19.539 2.120-4-2-2 17.129 2.710-4-2-
2、2 20.459 3.120.,*.ins 文件:组成上可分成五部分,TITL 041203e in P2(1)/c CELL 0.71073 11.5870 19.3080 17.2144 90.000 108.471 90.000 ZERR 4.00 0.0069 0.0112 0.0102 0.000 0.009 0.000 LATT 1 SYMM-X,0.5+Y,0.5-Z SFAC C H N O Co S UNIT 108 104 40 20 8 16,信息区,SIZE 0.47 0.43 0.35HTAB 2 L.S.8ACTABOND FMAP 2PLAN 20,晶胞参数,X光
3、波长,标题,空间群,单胞中分子数目,晶胞参数的标准偏差,晶格类型,对称操作码,原子类型,单胞中原子数目,命令区,可输入各种让XS和XL执行的命令,DFIX 1.43 0.02 o5 c27 WGHT 0.074900 1.631200 FVAR 0.169240TEMP 25,Co1 5 0.639436 0.182296 0.778299 11.000000 0.037950=0.043120 0.052330-0.002120 0.007180-0.002370Co2 5 0.132299 0.088237 0.782543 11.000000 0.061270=0.080610 0.08
4、2660-0.004140 0.011300 0.005080N1 3 0.816545 0.181404 0.857527 11.000000 0.048210=0.044720 0.042580-0.001460 0.010450-0.001390N2 3 0.641818 0.098803 0.864065 11.000000 0.055810=0.050870 0.070080 0.000700 0.013560-0.009990.,命令区,HKLF 4,END,原子表坐标、占有率、温度因子,衍射点文件类型,文件结束命令,3.其它文件,reslstpltciffcfpcftex,xs、
5、xl、refine产生的文件,记录xs、xl、refine过程和结果的文件,XP中做的图形文件,晶体学信息文件,结构因子文件,晶体结构报表文件,记录仪器型号、晶体外观等的文件,4.INS文件的建立和更新,结构解析和精修的过程,是ins文件建立和不断更新的过程,这主要是下列过程实现的:xprep、xshellrefine、xl、xp、edit、copy,SHELXTL的主要子程序和文件,二、常用的XS和XL指令,三、SHELXTL结构分析的步骤,1.项目的设立,打开SHELXTL程序:点project new,输入文件名,查找到hkl 文件并打开,2.结构的解析,1)结构解析的基本原理,XS用直
6、接法或Patterson法解决相角问题,试验性找出部分原子或重原子的位置(坐标),*所谓直接法(direct methods),就是运用数学的方法,利用不同衍射点的关系,从大量强度数据中,直接找出各个衍射点的相角,从而达到解析晶体结构的目的。其过程概括如下:,a 将|Fo|转化为归一化结构因子|Eo|,b 建立可以利用正切公式的三相角及四相角关系,c 赋于起始相角,d 利用正切公式精修相角,e 计算诊断指标,判断各套相角的质量,f 采用诊断指标最佳的相角数据计算解析电子密度图,即E图,*Patterson法是其本人1934年提出,通常只用来解析含有重原子的结构,用这种方法时,首先利用重原子的特
7、征峰,即Harker峰,求出重原子坐标,再通过Fourier合成获得其它原子的坐标,Harker峰,或Harker截面,就是同一套等效的原子组成的Patterson峰,由于平方效应,重原子的Harker峰会显得十分突出,寻找起来一般比较容易。因此,可以轻松地从分析重原子的Harker峰,得到重原子的坐标,2)用XS和XSHELL程序定出结构雏形(初始套),一般先试直接法,再试Pattson法,用什么方法是通过改变Edit.ins文件中的指令来实现的,直接法:,TITL 020908b in C2/c CELL 0.710730 30.1927 8.5175 13.9108 90.0000 95
8、.1300 90.0000 ZERR 8.00 0.0146 0.0042 0.0071 0.0000 0.0100 0.0000LATT 7 SYMM-X,Y,0.5-Z SFAC C H N O Cr UNIT 144 112 24 56 8TEMP 25TREFHKLF 4END,Patterson法:,TITL 020908b in C2/c CELL 0.710730 30.1927 8.5175 13.9108 90.0000 95.1300 90.0000 ZERR 8.00 0.0146 0.0042 0.0071 0.0000 0.0100 0.0000LATT 7 SYMM
9、-X,Y,0.5-Z SFAC C H N O Cr UNIT 144 112 24 56 8TEMP 25PATTHKLF 4END,继续尝试直接法:,TITL 020908b in C2/c CELL 0.710730 30.1927 8.5175 13.9108 90.0000 95.1300 90.0000 ZERR 8.00 0.0146 0.0042 0.0071 0.0000 0.0100 0.0000LATT 7 SYMM-X,Y,0.5-Z SFAC C H N O Cr UNIT 144 112 24 56 8TEMP 25ESEL 1.5 3.5TREFHKLF 4END
10、,归一化结构因子E的 下、上限值,可试探性设置,设置好ins文件,点击XS,就可自动进行计算,XS计算结果的评估,#直接法,RE越小越好,一般大于0.3,就预示着不成功,#打开Xshell,检查计算出的原子或(Q)峰位置是否构成合理的化学结构,#精修(refint)一次后,R1值是否大小于0.5;或继续精修R1值是否很快降低,*如果得到的原子或峰的位置能够构成合理的化学结构,就要通过删除不合理的原子或(Q)峰,并把Q 峰选择后定为相应的原子。操作方法是:,选择原子或(Q)峰的方法:,a 光标指在欲选择位置,点S 键。可依次选择所有的要选择的原子或Q 峰,c 用Select菜单(或背景右手键点“
11、Select”)可选择部分和所有的原子或Q 峰,顺序与原顺序相同,b 光标指在欲选择位置,右手键点“Select”,可依次选择原子,d 按住“Shift”键,左手键拖出一个选择框并同时轻开,可选择一个或多个原子,删除原子或(Q)峰的方法:,a 光标指在欲删除位置(原子或键),点K 键,b 光标指在欲删除位置,右手键点Bonds and Angles,再在对话框中点“Delete”,e 点refine菜单,系统会提示有Q 峰没命名,并选择所有Q 峰,顺序与原来顺序相同,e 在ins文件删除,c 光标指在欲删除位置,右手键点“Delete”,d 先选择,再点Select菜单中的Kill selec
12、ted,恢复删错的原子或(Q)峰的方法:,a 先选择,再用Lablels菜单中的Group label,命名(标记)原子或(Q)峰的方法:,点“EditRestore Killed Atom”,先选择要恢复的原子(或Q峰),再点“Restore”,b 光标指在要命名的位置,用右手键的Edit,c 光标指在欲命名位置,右手键点Bonds and Angles,再在对话框中点“Rename”,*如果仍得不到较为合理的初始结构,则应试验用三种方法中的另一法,或试验用不同的ESEL值,*如果试验用各种方法,仍不行,那就要考虑是不是出现了如下的问题:,常用 的小技巧:,a 解析和精修结构 的初期,可把C
13、、O、N都定为C,b 如果独立单元中分子较多,可光标指在 一个分子中的某一原子上,再用右手键菜单中的UNIQ指令把其分离出来处理,c 如果独立单元只是分子的一部分,原子不好取舍,可用右手键菜单中的GROW指令把分子长全处理,结构的解析中的若干问题,A 晶胞的错误确定,用CCD探测器,仪器可以自动确定晶胞(分别在matrix、sadabs和xprep中),但有时会不正确,因而晶系和空间群会因之错误,解决方法:在xprep中重新确定晶胞;实在不行,需用收集数据重新做晶胞,B 空间群的错误指认,主要的三种情况:,a 晶胞的错误确定引起的,b 仪器只根据Rint和CFOM(诊断因子)值 来确定,有时会
14、降低晶体的对称性,c 由系统消光规律的错误判断引起,解决方法:在xprep中重新确定空间群;重新做晶胞,*XPREP的主要功能和应用,单击进入XPREP程序,根据程序的提示输入晶胞参数,选择可能的晶格,程序则显示以下菜单:,D Read,Modify or Merge DATDSETSP Contour PATTERSON SecionsH Search for HIGHER mertric symmetryS Determine or input SPACE GROUP,读入、更改、合并衍射数据,计算显示Patterson截面,寻找更高的对称性,确定或输入已知的空间群,A Apply ABS
15、ORPTION correctionsM Test for MEROHEDRAL TWINNINGL Reset LATTICE type of Original CellC Define unit-cell CONTENTSF Setup SHELXTL FILESR RECIPROCAL Space DisplaysU UNIT-CELL transformationsT Change TOLERANCESO Self-rotaion functionQ Quit Program,吸收校正,重设原始晶胞的晶格类型,孪晶缺面试验,定义单胞的化学组成,建立计算指令文件,显示倒易空间,转换晶胞,
16、改变一些变量的容忍值,自旋函数,退出程序,Select option H:,程序显示目前的晶胞参数和其它的可能选择,可认同程序的选择,但有时需自行选择,Select option A:,认同程序的选择后,程序会提示下一选项S,Select option S:,可直接输入已知空间群,也可测定,Select option S:,则还要进一步选择晶系,Select option M:,Select option P:,进一步选择晶格类型,程序显示目前的空间群和其它的可能选择,可认同程序的选择,但有时需自行选择,认同程序的选择后,程序会提示下一选项D,Select option A:,Select o
17、ption D:,程序让显示统计强度的选项S,程序让显示合并等效点的选项A,Select option S:,Select option A:,程序显示出退出主菜单的选项E,Select option E:,认同程序的选择后,程序会提示下一个定义单胞的化学组成选项C,Select option C:,认同原有的或更改或重新输入后,程序又显示出退出主菜单的选项E,Select option E:,认同程序的选择后,程序会提示下一个建立SHELXTL计算指令文件选项F,Select option F:,认同程序的选择后,程序会提示输入文件名,输入并回車后,程序会问是否写hkl文件,同意后,程序又提
18、示下一个选项Q,文件名,可任意取,Select option Q:D,程序又显示出退出主菜单的选项E,Select option E:H,输入0.84以使为25,程序又显示出退出主菜单的选项E,这里也可不限制,解或精修结构时用OMIT指令限制,Select option E:W,再输入文件名后,两个回车后又显示如下选项:,应与前面的文件名相同!,Select option E:M,便 会显示如下选项,Select option A:I,就可一直回车选项退出XPREP程序,C 严重给错分子式,不同的原子对X-射线散射能力不同,如果给错了分子式,尤其是重原子错,一般都 会导致解析不了结构(需仔细核
19、实分子式),D 赝对称结构,就是晶体结构中含有某些较重的原子,它们处于特殊位置时,使得自身的对称性高于轻原子及整体的对称性(如具有额外的平移和中心对称性),这时就会造成在电子密度图上轻原子的多解,即轻原子除了在正确位置出现个,还在一个或多个错误位置出现(要仔细的重解,仔细地确定原子),E 衍射数据质量差,在没有以上错误,仍解不出结构时,可能是数据质量太差,是由以下原因造成的:,其一、衍射强度数据太弱;数据还原程序统计得到的平均信噪比I/(I)太小,而解结构程序统计的平均信噪比的参数Rsigma太大,其二、吸收校正没做好,所有等效点的平均差别Rint较大,解决方法:重新培养较大的单晶;重收数据,
20、解决方法:重新用其它方法做吸收校正,如果没有以上问题,还是解不出结构,还可试用先将有心改为无心空间群的方法,解出初始结构,再转换成有心空间群,3.结构精修的基本过程,在得到初始结构后,所有的Refine(或XL)计算都 是由计算机完成的,对于非晶体学专家,该过程有如:“暗箱操作”,但差值Fourier电子云密度图的分析和结构模型的逐步完善必须人工参与完成,A 读入晶胞参数、空间群、分子式、衍射数据等,B 程序利用已知的部分原子的坐标和种类,计算 出结构的部分正确的相角,C 程序结合计算出结构的的部分正确的相角和测 的衍射强度,模拟出“观察”到的电子云密度图,D 程序把已经确认的原子从“观察”到
21、的电子云密 度中扣除,计算出差值Fourier电子云密度图,E 在差值Fourier图中,找出分子中其它原子的位置,F 程序根据衍射强度,对原子坐标和位移参数等结 构参数进行精修,G 重复CF步骤至分子中所有原子都被指认、精修 收敛,*结构精修的典型步骤,A 在XSHELL中用Refine精修各向同性的非氢原子,定出基本正确的结构模型后,在前几轮精修(各向同性)中,可以通过检查R1和U值,来检查原子确定的是否正确,并逐步指认出非氢原子,R1应小于0.4(wR2为0.50.7),越小越好,U应在0.1以下,如果在精修过程中,U值急剧增大,说明该原子并非存在,如果U偏大,该位置可能是一个更轻的原子
22、,如果U偏小,则可能是一个更重的原子,还要注意残余峰的大小,如果其大于1就应考虑是不是还有未发现的非氢原子,在Refine中选择Aniso则可进行各向异性精修(也可在Edit.ins文件中加ANIS指令,用XL进行),这样可更确切地表达电子云的真实分布情况,所以R1值会明显下降,B 对所有的非氢原子进行各向异性精修,这时如果是非中心对称的空间群,应检查Flack参数是否接近于0,否则要进行绝对构型转换,指令:,MOVE 1 1 1-1,对于具有确定立体结构的有机基团,可采用理论加氢;对于无法理论加氢的,如水分子,可从差值Fourier图中找出氢原子,参加精修的情况视 数据的质量而定,C 计算/
23、指认并精修氢原子,理论加氢的方法,退出XSHELL,在Edit-.ins文件的命令区,添加指令:,HFIX mn 需加氢原子名(产生AFIX固定),m 是一或两位数,指定氢的类型:=1 叔-H,=2 仲-H,=3(或13)伯-H,=4 芳-H,=8(或14)X-O-H,=9 X=CH2或X-NH2,=15 笼状B-H,n 是一位数,指定固定的类型:=1 坐标、占有率、位移因子固定,=2 占有率、位移因子固定,=3(或7)坐标固定,=4 同3,但允许修正X-H的键长(方向固定),然后进行XL精修,完成后进入XSHELL继续精修,也可在XSHELL中先点Atoms-Hybridize All 再点
24、Atoms-Calculate Hydrogens,然后检查H加的是否合理,如不合理,可打开Edit-Edit Current File修改不合理的部分,再用Refine继续精修,还可在XSHELL中先选择同类原子,再用Select Atoms Set HFIX分别进行加氢,D 进一步精修更多的晶体学参数,先在ins文件中写入要精修的参数指令,再进行精修,例如:,SIZE 0.38 0.32 0.24,HTAB 2,CONF,CONF X1 X2 X3 X4,BOND$H,BIND X1 X2 X3 X4_$1,MPLA n X1 X8,EQIV$1 X-1,Y+1/2,-Z,OMIT 1-2
25、 4,输入晶体的尺寸,计算氢键并写入lst文件,计算所有的扭转角,计算这四个原子间的扭转角,加上$H,可计算包括H在内的各种健长和键角,计算输入的原子对的健长和键角,计算输入的前n个原子的平面,并给出所有输入原子偏离该平面的距离,还能给出与前一个平面间的夹角,提供所键合原子的对称操作码,MOVE xi yi zi,删除h、k、l 分别为1、-2、4的坏点,移动坐标,EXTI x,精修消光效应,x 是程序自动产生的,这些计算一般在作数据表时进行,找出对称操作码的方法:,而用ENVI命令会得到有关的用四位数字表示的对称操作码(如:2746、3654、4545等),可表示为:SXYZ,其中,S为对称
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 用SHELXTL程序进行晶体结构分析的方法 单晶结构分析电子教案 SHELXTL 程序 进行 晶体结构 分析 方法 结构 电子 教案
链接地址:https://www.31ppt.com/p-2999286.html