里特沃爾德精修(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.