里特沃尔德精修(Rietveld refinement),简称里特沃尔德法(Rietveld method)[1],是由荷兰晶体学家雨果·里特沃尔德提出的一种晶体材料表征技术,是全谱拟合法的一种。根据实验测得的多晶衍射数据全谱信息,利用物质的晶体结构参数(晶胞参数、原子坐标等)和非结构参数(峰宽、择优取向因子等)信息模拟计算并利用非线性最小二乘法不断拟合接近实验谱图,从而得到晶体结构信息的方法[2][3]。
里特沃尔德精修还有里特沃尔德全谱拟合(Rietveld whole pattern fitting)、里特沃尔德图谱拟合(Rietveld profile fitting)、里特沃尔德结构精修(Rietveld structure refinement)等其他名称[4][5]。
历史
1969年,荷兰晶体学家里特沃尔德(Hugo M. Rietveld)首次提出“全谱拟合”的概念,并应用到中子粉末衍射技术中[6]。里特沃尔德先后用ALGOL语言和Fortran语言编写出相关程序并将其自由分享,得到当时晶体结构研究者内使用[7]。
1977年,奇达姆(A.K. Cheetham)[8]对中子衍射图谱的拟合法进行了归纳总结。在此之前全谱拟合技术仅局限于中子衍射技术中[9]。
1979年,库珀(M. J. Cooper)将结构参数和峰型参数同时拟合改进为先拟合峰型参数,后拟合结构参数,使得拟合结果更合理[10]。 杨格(Robert Alan Young)等人将里特沃尔德法应用到了晶体X射线衍射中,包括同步辐射X射线衍射技术[11][12]。
随后里特沃尔德精修范围从简单到化学体系不断扩大到复杂体系的分析,如硅酸盐复合水泥的物相分析,蛋白质结构确定[7]等。经过不断完善与发展,里特沃尔德法已经发展成为获得晶体结构信息的重要方法[2][9],并诞生一系列相关软件进行里特沃尔德精修[13][14][15][16]。
原理
实验测得的粉末衍射图案常常取决于一系列因素,包括:
因此粉末衍射图谱可以由晶体结构参数和仪器因素(峰型参数)确定,并通过计算叠加得到完整图谱。
里特沃尔德精修是基于散射能量守恒的原理,在假设晶体结构模型和结构参数基础上,结合峰形函数来计算多晶衍射谱,并用最小二乘法不断调整结构参数与峰形参数,使计算得到衍射谱不断逼近实验测得的谱图[2],使拟合偏差最小,从而得出相关晶体结构信息。
精修过程
里特沃尔德精修主要经过全谱计算、拟合逼近、拟合评价三个步骤进行。
在给定假设晶体模型下,衍射谱可以由布拉格方程、衍射强度公式、本底函数计算叠加得出。
在给定晶体结构和入射束波长λ下,根据布拉格方程可以确定各衍射晶面对应的衍射角。
根据衍射积分强度公式[18]:
其中:
- : 尺寸因子
- : 多重因子
- : 洛伦兹因子
- : 偏振因子
- :吸收因子
- : 择优取向因子
- : 消光因子
- :晶胞结构因子
可以得到晶面衍射峰对应积分强度。随后加上归一化峰型函数,就得到了第个衍射峰处的实测强度
设定背景强度函数,将各衍射峰强度叠加即可得到全谱各点处的计算实测强度分布
- 其中,代表第个衍射峰。
根据计算得到的衍射谱强度分布,不断调整相关参数,利用非线性最小二乘法使得与实验实际测得的衍射谱强度分布的方差最小:
其中为权重因子。
为了确定精修拟合参数的优劣性,常用可信度因子来判定参数的正确性:
- 谱图残差,又称可靠因子(Profile residual,Reliability factor)
- 加权谱图残差(Weighted profile residual)
- 布拉格残差或强度残差(Bragg residual,Intensive residual)
- 期望谱图残差(Expected profile residual)。其为期望值
- 拟合优度,或 (Goodness of fitting)
和是根据峰实测强度得来,其中最能反应拟合的优劣。
是根据积分强度得到,与结构模型高度相关,用于判断结构模型和合理性。
拟合函数
完美理想衍射条件下(晶粒无穷小,完全随机取向,仪器几何系统完全准直,无发散纯单色光),得到的衍射图谱符合理想衍射强度公式,且衍射峰为线型[18]。然而衍射测定中不可能达到理想条件,故需要引入相关拟合函数进行修正[2][9]。
峰型函数用于模拟实际测定中衍射峰偏移理想衍射条件引起的峰变形。
里特沃尔德最早采用高斯函数对中子多晶衍射进行精修[6]。洛伦兹函数也是早期常用的峰型函数[11]。由于两种函数均不能很好描述峰型,后续在两种函数基础上进行卷积、线性组合或改进得到一些列新的峰型函数[9]。
函数名称 | 简写 | 备注 |
---|---|---|
高斯函数(Gaussian) | 对于峰型,峰顶太宽,峰尾太窄 | |
洛伦兹函数(Lorentzian) | 对于峰型,峰顶太窄,峰尾太宽 | |
居间洛伦兹函数(Intermediate Lorentzian) | ||
变形洛伦兹函数(Modified Lorentzian) | ||
皮尔逊VII函数(Pearson VII) | 由和的线性组合得到 | |
福格特函数(Voigt) | 由个和个卷积得到 | |
赝福格特函数(Pseudo-Voigt) | , 其中 | |
余弦洛伦兹函数(Cosine-Lorentzian) | 自变量为衍射角的余弦值的洛伦兹函数形式 | |
学生t-分布(Student's t-distribution) | 高角度和低角度下所用函数系数不同 |
峰宽函数用于模拟实际测定中衍射峰偏移理想衍射条件引起的峰宽化,常用半峰宽表示。针对不同的衍射条件提出了不同的峰宽函数形式。
本底函数用于模拟衍射谱中的背景。包括X射线荧光、非相干散射以及仪器噪声等机械误差;分子热振动引起的漫散射;以及非晶体成分造成的背景。
择优取向函数用于校正非球形(特别是棒状和片状)晶粒在制样过程中朝向分布不均造成的衍射强度分布不均。
最早采用指数分布分布形式模拟晶体的择优取向分布[6],后面又发展出更多择优取向函数形式来适应复杂的情况[2]。
已经提出的择优取向函数:
其中为拟合参数,为择优取向与衍射晶面之间夹角。
应用与局限
里特沃尔德精修常用于用于晶体结构分析,包括薄膜材料,纳米材料等低维材料[2]以及生物材料[7]的分析:
里特沃尔德精修是全谱拟合中属于需要相关已知晶体结构的数据的一类,另一类则不需要相关结构数据,需要相关纯物质标准谱。因此里特沃尔德法不能测定未知结构的晶体粉末样品,而且需要一个较为准确的初始晶体模型,计算过程相对复杂[2][9]。
在X射线衍射谱图精修中,还有Pawley法,Le Bail法等变种。 Pawley法不需要结构模型,但参数多,耗时,误差较大,难以获得晶体结构内层原子信息。Le Bail法参数少,收敛快,准确度较高。但不能区分位置较近的衍射峰[23][24][15]。
相关条目
参考资料
参考书籍
外部链接
Wikiwand in your browser!
Seamless Wikipedia browsing. On steroids.
Every time you click a link to Wikipedia, Wiktionary or Wikiquote in your browser's search results, it will show the modern Wikiwand interface.
Wikiwand extension is a five stars, simple, with minimum permission required to keep your browsing private, safe and transparent.