搜索
X X 射线粉末衍射射线粉末衍射测定晶体结构方法和测定晶体结构方法和 PowderX PowderX 软件介绍软件介绍董成董成超导国家重点实验室超导国家重点实验室中国科学院物理研究所中国科学院物理研究所20042004年年22月月1212日凝聚态物理前沿讲座日凝聚态物理前沿讲座粉末衍射结构分析的科学意义粉末衍射结构分析的科学意义（（11））自然界存在的和人工合成的绝大多数固体材料是多自然界存在的和人工合成的绝大多数固体材料是多晶体晶体((如南开大学每年合成新络合物如南开大学每年合成新络合物100100－－200200种，仅有种，仅有几种能长出单晶几种能长出单晶))。。（（22））由于易生成孪晶、包晶、生长条件苛刻等原因由于易生成孪晶、包晶、生长条件苛刻等原因,,很多材料的单晶生长都不容易甚至不可能很多材料的单晶生长都不容易甚至不可能。。（（33））新材料研究中最先得到的大多是多晶。用多晶粉末新材料研究中最先得到的大多是多晶。用多晶粉末衍射数据测定出晶体结构衍射数据测定出晶体结构,,有助于研究结构和性能的有助于研究结构和性能的关系。关系。（（44）粉末衍射能）粉末衍射能方便进行高温、低温、强电磁场、高压方便进行高温、低温、强电磁场、高压下的实验下的实验,,研究物质的相变。研究物质的相变。国内外研究现状用粉末衍射数据测定晶体结构近年来是国际晶体学界研究的热点，因为：（1）衍射实验条件的改善，如同步辐射光源的发展；（2）新的结构分析算法的发展；（3）计算机计算能力的发展。目前国际上已经可能用粉末衍射数据和从头计算方法测定单胞体积为2500Å 3、有200个原子参数的化合物晶体结构。但现在的粉末衍射结构分析面临很多挑战性的难题，远不是常规工作。现国内外还很少用超级计算平台解决粉末法测定晶体结构问题。国内外研究现状国内外研究现状到到20012001年的粉末法测定的晶体结构统计数字年的粉末法测定的晶体结构统计数字((累计累计592)592)国内外研究现状国内外研究现状近年来代表性的粉末衍射结构分析软件程序名称作者算法EXPO2000:意大利A.Altomare等,直接法＋MontCarloEndeavour：德国H.Putz等,综合优化Rwp和系统势能DASH：英国WIFDavid，K.Shankland，SAESPOIR：法国ArmelLeBail，MontCarloSAEAGER：英国KennethD.M.Harris，GAPSSP：美国P.W.Stephens，SAPowderSolve：英国G.E.Engel等，SA[药]氯丙嗪, 冬眠硫磷 C 10 H 13 ClN 2 O 3 a=26.645, b=9.072, c=5.217Å,a=b=g=90, Space group=P2 1 2 1 2 1 样品制备粉末衍射数据处理指标化指标化结构模型结构模型结构精修研究目标和内容研究目标和内容研究目标和内容研究目标和内容（（1 1）研究国际上现有的结构测定算法和软件的优缺）研究国际上现有的结构测定算法和软件的优缺点，点，建立起有自主知识产权的粉末衍射法晶体结构建立起有自主知识产权的粉末衍射法晶体结构测定软件系统和发展新算法。测定软件系统和发展新算法。争取经过几年时间的争取经过几年时间的努力，赶上国际同行的水平并争取在某些方面有所努力，赶上国际同行的水平并争取在某些方面有所创新。创新。（（2 2）面向粉末衍射结构分析的难点和瓶颈问题上，）面向粉末衍射结构分析的难点和瓶颈问题上，解决一些现有程序和方法解决不好或解决不了的问解决一些现有程序和方法解决不好或解决不了的问题：题：指标化、重叠峰强度分离与求解初始结指标化、重叠峰强度分离与求解初始结构模型。构模型。研究目标和内容研究目标和内容（（3 3））粉末法测定晶体结构的超级计算平台软件系统要粉末法测定晶体结构的超级计算平台软件系统要力求算法先进、结果可靠、使用方便和界面友好。充力求算法先进、结果可靠、使用方便和界面友好。充分发挥超级计算平台的作用，争取能够使解决问题的分发挥超级计算平台的作用，争取能够使解决问题的规模和成功率比用规模和成功率比用PC PC机和工作站运算明显提高。机和工作站运算明显提高。（（4 4）考虑到工作的难度和探索性）考虑到工作的难度和探索性, ,准备分两个阶段实行。准备分两个阶段实行。在第一阶段（在第一阶段（2003 2003－－2005 2005年），首先解决从无到有的年），首先解决从无到有的问题，建立起在超级计算机上运行的粉末结构分析的问题，建立起在超级计算机上运行的粉末结构分析的基本程序和计算平台；基本程序和计算平台；2006 2006年以后几年继续发展完善年以后几年继续发展完善本软件系统，使其达到更高水平。本软件系统，使其达到更高水平。研究方法和技术路线－研究方法和技术路线－有所为，有所不为有所为，有所不为 (1) (1) 有成熟算法和免费共享源代码的，作适当改进后移有成熟算法和免费共享源代码的，作适当改进后移植到超级计算机上。植到超级计算机上。 (2) (2) 现有算法可行，但没有源代码时，就独立开发相应现有算法可行，但没有源代码时，就独立开发相应的源代码。的源代码。 (3) (3) 发展新创的结构解析算法及其源代码程序，是这项发展新创的结构解析算法及其源代码程序，是这项研究的关键和难点研究的关键和难点，主要的时间和精力要花在这一方，主要的时间和精力要花在这一方面。面。 (4) (4) 综合坐标空间和衍射空间循环迭代法结构求解和结综合坐标空间和衍射空间循环迭代法结构求解和结构精修方法构精修方法。。要注意充分发挥超级计算机的并行功要注意充分发挥超级计算机的并行功能，并确保各个模块之间的数据传输方便可靠。能，并确保各个模块之间的数据传输方便可靠。粉末衍射图研究方法和技术路线－研究方法和技术路线－模块化模块化衍射数据处理由晶体结构计算粉末衍射图衍射图指标化空间群分析峰形拟合提取衍射强度 Patterson 函数法直接法结构解析电子密度计算和寻峰坐标空间结构约束表示和变换坐标空间搜索模拟退火遗传算法化学结构信息结构模型 Rietveld精修粉末衍射实验方法粉末衍射实验方法研究方法和技术路线－研究方法和技术路线－数据处理数据处理（1）扣除Ka2（2）扣除背景（3）数据平滑（4）寻峰（5）校正系统误差（6）指标化（7）峰形拟合（8）数据格式转换研究方法和技术路线－研究方法和技术路线－指标化指标化多晶粉末衍射图的指标化方法是寻找在实验误差范围内满足以下方程的解，从衍射峰对应的晶面间距求解晶格参数 (a,b,c和a,b,g)，同时确定衍射峰的晶面指标。 (d hkl ) 2 =h 2 a* 2 +k 2 b *2 +l 2 c* 2 +2hka*··b*+2lhc*··a*+2klb*··c* 令Q i = (d hkl ) 2 , A=a* 2 , B=b *2 , C=c* 2 , D=2a*··b*, E=2b*··c* F=2c*··a*,有： n i h Fl l Ek k Dh Cl Bk Ah Q i i i i i i i i i i ,..., 1 2 2 2=+++++=常用的方法和计算程序是面指数尝试法(TREOR90)，晶带法(POWDER)和对分法(DICVOL91)。指标化方程数学上是多解的, 所以人们常用品质因数 (Figure of Merit) 来表征指标化结果的可靠性。研究方法和技术路线－研究方法和技术路线－指标化指标化尽管已有多种指标化计算程序，但指标化还是粉末衍射结构分析的瓶颈。近来又有一些新发展的指标化算法：全谱MonteCarlo和格点搜索方法，程序名MCMaille，作者LeBail。每秒试验20000次，仍然需要很多小时计算。考虑系统消光影响的X-Cell，作者MarcusA.Neumann和利用SVD方法的指标化方法(Bruker商业软件)等。我们准备在并行计算机上开发能够容忍少量杂相衍射线条的指标化程序；并改进现有的指标化算法并发展新的算法。近来我们提出用欧几里德算法进行高对称性晶系衍射图指标化，并成功地指标化了LeBail等发起的国际循环赛中立方晶系晶胞参数为18.28Å的化合物衍射图。方法的优点是在有杂相存在时也能得出正确结果，受杂线影响很小。研究方法和技术路线－研究方法和技术路线－空间群分析空间群分析根据晶体中存在带心点阵和螺旋轴及滑移面导致的系统消光，可以在指标化后判定晶体的点阵类型和所属的X射线衍射群。结合某些物理性质的测量和衍射强度统计规律可以帮助确定空间群。具有旋光性的晶体不具有对称中心；具有压电性的晶体不具有对称中心（O-432点群除外）；具有热电性的晶体含有一个极性轴；具有倍频效应的晶体不具有对称中心（O-432，D6-622和D4-422点群除外）在很多情况下，低对称性空间群可以包括高对称性的空间群。空间群的确定一般不构成粉末结构分析的主要障碍。空间群的确定一般不构成粉末结构分析的主要障碍。研究方法和技术路线－研究方法和技术路线－基本参数法峰形拟合基本参数法峰形拟合提取衍射强度提取衍射强度基本参数法是根据衍射仪的几何参数和光源特征参数通过卷积来精确计算衍射峰形的方法，可以得到如同在“理想”衍射仪上的精确峰位，对准确提取衍射强度也很有好处。但现在使用该方法的大多是商业软件（BrukerTopasTopas和BGMNBGMN），），没有源代码。我们准我们准备自己开发出基本参数法峰形备自己开发出基本参数法峰形拟合的源代码。拟合的源代码。研究方法和技术路线－研究方法和技术路线－由晶体结构由晶体结构计算粉末衍射图计算粉末衍射图其中：Y ik －反射k对净强度的贡献； s －标度因子，同一实验条件下为常数； m k －多重性因子，由晶体对称性决定的整数； P k －择优取向因子； L k － Lorentz 偏极化因子，是衍射角的函数； F ck －结构因子； f(x)－峰形函数。 F(hkl) = Σ g i f i exp(2πi[hx i +ky i +lz i )T研究方法和技术路线－研究方法和技术路线－直接法直接法直接法的基本原理根据实验衍射强度只能得到结构因子的模 结构振幅，而不能得到结构因子的相角。 F I 2µ但实验中可得到比未知结构参数多的衍射强度数据, 从而结构振幅和相角之间存在相互关联。 ) . 2 cos( sin ) ( ) . 2 cos( cos ) ( r r j j j H j j j H H f H F H f H Fppffåå==直接法是根据结构振幅与相角的相互关联, 推演出实验中无法测定的相角信息, 进而测定晶体结构的方法。直接法是公认最成功, 应用最广的晶体结构测定方法。由于直接法对人类认识原子水平上的微观世界的巨大贡献, 两位早期从事直接法研究的科学家H. Hauptman和 J. Karle 获得诺贝尔化学奖。研究方法和技术路线－研究方法和技术路线－直接法直接法古典的直接法古典的直接法基于基于 3 3 个基本的个基本的假设：假设： (1) (1) 电子密度在单胞总是正的。电子密度在单胞总是正的。 (2) (2)单胞由分离的原子组成。单胞由分离的原子组成。 (3) (3)原子的位置是在整个单胞内原子的位置是在整个单胞内均匀分布的随机变量。衍射均匀分布的随机变量。衍射强度数据个数大大超过未知强度数据个数大大超过未知原子坐标参数，可以用概率原子坐标参数，可以用概率论推导衍射的相角。论推导衍射的相角。归一化结构因子的模 The tangent formula The tangent formula S 1 关系：U U 2H 2H ³³22U U HH22--11三重相角关系：ff H H + + ff K K + + ff H H K K »»00研究方法和技术路线－研究方法和技术路线－直接法直接法1.数据还原,计算归一化结构因子|EH|。选出|EH|较大的衍射。2.规定三个确定原点的衍射的相角,并根据不等式关系，符号关系等得出(或根据随机