リートベルト解析のやり方と手順・結晶構造精密化の基本

リートベルト解析のやり方を、XRD測定から構造モデル作成・RIETAN-FP操作・R因子確認まで徹底解説。金属加工現場での定量分析や品質管理への応用もわかりやすく紹介。あなたの職場でもすぐ使えるノウハウを詳しく知りたくありませんか?

リートベルト解析のやり方・手順を基礎から実践まで解説

パターンを眺めるだけでは、結晶相の「割合」はわかりません。


この記事でわかること
🔬
リートベルト解析とは何か

XRD測定で得られた回折パターンを使い、結晶構造パラメータを最小二乗法でフィッティングする解析手法です。格子定数・原子位置・定量比まで一度に得られます。

⚙️
RIETAN-FPを使ったやり方の手順

試料調製→XRD測定→CIFファイル取得→INSファイル編集→BAT実行→R因子確認という6ステップの具体的な進め方を解説します。

🏭
金属加工現場での活用ポイント

製造プロセス管理・多相試料の定量・劣化原因特定など、金属加工業界でリートベルト解析が実際に役立つ場面と注意点を紹介します。


リートベルト解析のやり方の前に知っておきたい基本原理



リートベルト解析(Rietveld analysis)は、オランダの結晶学者ヒューゴ・リートベルト(Hugo Rietveld)が1969年に考案した粉末X線回折(XRD)データの解析手法です。測定で得られた回折プロファイル全体を、結晶構造モデルから理論計算したパターンと最小二乗法でフィッティングすることで、結晶構造パラメータを精密化します。


従来の「ピーク位置だけを読み取る手法」と大きく異なるのは、ピーク強度・ピーク幅・ピーク形状・バックグラウンドを含むパターン全体を解析に使うという点です。このおかげで、ピークが重なり合う多相試料でも定量分析が可能になります。これが原則です。


具体的に得られる情報は次の通りです。格子定数・原子位置(xyz座標)・席占有率・温度因子・結晶相ごとの存在割合(質量比)・ゼロ点シフト・結晶子サイズ・格子歪みなど、非常に多くのパラメータを一度に得られます。たとえばAl₂O₃とZnOを混ぜた粉末試料を解析すると、Al₂O₃が50.56%・ZnOが49.44%といった具体的な質量比が±1%以内の精度で算出できます。


金属加工の現場では、熱処理や加工の前後で結晶構造がどう変化したかを数値化できるため、製造プロセスの最適化や材料劣化の原因特定に直結します。「パターンを眺めて似ている・似ていないで判断する」段階から、「数値として変化量を比較する」段階に移行できるのが最大のメリットです。


なお、解析には「非晶質(アモルファス)ではなく結晶性試料であること」「粒径が1〜5µm以下の粉末試料であること」「試料の結晶構造モデル(候補)があること」の3条件が必要です。この3点が条件です。


リートベルト解析(Rietveld analysis)とは | MST(材料技術研究所)- 概要・解析事例・得られる情報の詳細


リートベルト解析のやり方・試料調製とXRD測定の手順

リートベルト解析の精度は、解析ソフトの使い方よりも前段階の「試料調製」と「XRD測定条件」に左右されます。ここを疎かにすると、どれだけ丁寧にパラメータを調整しても正確な結果は得られません。試料調製が先決です。


まず粉砕・粒径調整から始めます。粉末法では試料中の結晶の向きがランダムに分布していることが前提で、直径50µm程度の球を直径10mm・長さ0.1mmの円筒に1万個以上入れるイメージの充填密度が必要です。実際には試料を乳鉢で細かく粉砕し、粒径を1〜5µm程度(髪の毛の直径は約70µmなので、それより一桁以上細かいレベル)に揃えます。


重要な注意点が「選択配向(テクスチャー)」の抑制です。結晶が特定の方向に揃ってしまうと、ある回折ピークだけが異常に強くなり、定量結果が大きく狂います。金属粉末の場合、乾式粉砕よりもエタノールを用いた湿式粉砕の方がメカノケミカル現象(粉砕中に化学反応が起きる現象)を抑えやすく、選択配向も出にくいとされています。


試料板への充填は、表面が平坦になるよう摺り切りまたはプレスで行います。表面が凸凹していると入射X線がはみ出し、特に低角側の強度データが狂います。厚みが薄い試料の場合は、X線が基板まで突き抜けてしまうこともあるので注意が必要です。


XRD測定条件では、回折ピーク強度の平均が1,000カウント以上になるよう計数時間を設定することがJIS R1640(窒化けい素の相組成分析方法)でも要求されています。カウント数が少ないままだと統計誤差が大きくなり、解析精度が下がります。測定後は、回折パターンをRIETAN-FP対応の`.int`ファイル形式に変換して保存します。


リートベルト法による回折データの精密解析(JFEテクノリサーチ)- 複相試料の定量分析事例・得られる情報の一覧


リートベルト解析のやり方・RIETAN-FPを使ったINSファイル設定の手順

日本で最も広く使われているリートベルト解析ソフトがRIETAN-FPです。物質・材料研究機構の泉富士夫博士が開発し、指定ルールに従えばインターネットから無料でダウンロード・利用できます。ここでは、実際の操作手順を順番に説明します。


ステップ1:CIFファイルの取得と構造モデルの準備


CIFファイル(Crystallographic Information File)は、結晶構造の格子定数・空間群・原子座標などを記述したテキストファイルです。物質・材料研究機構が提供する無機材料データベースAtomWork(MatNavi)は無料登録で利用でき、多くの金属・セラミックス材料のCIFファイルを入手できます。有料ですが信頼性の高いデータベースとしてICSD(無機結晶構造データベース)もあり、専門家の間では標準ツールとして使われています。


ステップ2:VESTAでINSファイルを自動生成する


VESTAは三次元可視化ソフトウェアで、RIETAN-FPと連携した入力ファイル作成に使えます。VESTAにCIFファイルを読み込み、メニューの「Utilities → Powder Diffraction Pattern」を実行すると、`VESTA\temp`フォルダ内にINSファイルが自動生成されます。これが初期入力ファイルになります。


ステップ3:INSファイルの主要パラメータを編集する


生成されたINSファイルをテキストエディタで開き、以下の箇所を修正します。


パラメータ 変更前(シミュレーション) 変更後(リートベルト解析) 役割
NMODE NMODE = 1(シミュレーション) NMODE = 0(リートベルト解析) 解析モードの切り替え
NUPDT NUPDT = 0 NUPDT = 1 INSファイルへのパラメータ更新
SCALE 初期値 0.0001 1(IDを動かすパラメータに設定) 計算値と実測値の全体比率合わせ
BKGD 全て 0.0 低角側のバックグラウンド実測値 バックグラウンド関数の初期値
FWHM3 仮の半値幅パラメータ 実測パターンから読み取ったW値 ピーク幅の角度依存性モデル
R12 0.5(Kα2あり) 0.0(Kα2なし・放射光の場合) Kα1/Kα2の強度比


INSファイルはtabキーを使ってはいけません。半角スペースで記述するのがRIETAN-FPのルールです。これは必須です。また、空間群の指定(HKLM1)を誤ると全く収束しないため、VESTAで確認した値をそのままコピー&ペーストすることが推奨されています。


ステップ4:BATファイルを実行して解析を開始する


`.bat`ファイルの`SET SAMPLE=`の部分をINSファイルのファイル名(拡張子なし)に書き換え、ダブルクリックして実行します。計算が完了すると、グラフビューア(RIETVIEW)に赤点(実測値)・青線(計算値)・緑線(差分)が表示されます。差分の緑線がほぼゼロになれば収束成功です。


多相試料の解析では、狭い散乱角範囲から合わせ込みを始め、少しずつ角度範囲を広げていくと収束しやすくなります。


KEK放射光利用技術入門コース「リートベルト法による粉末構造解析入門」- RIETAN-FPによるCeO₂・TiO₂の実習手順


リートベルト解析のやり方・R因子による収束判定と結果の読み方

パラメータを動かして解析を実行したら、次に確認するのがR因子(信頼度因子)です。R因子は「計算パターンと実測パターンがどれだけ一致しているか」を数値化したもので、解析の品質を判断する最も重要な指標です。R因子の確認が鍵です。


主要なR因子の種類と目安を以下にまとめます。


R因子 名称 目安 特徴
Rwp 重み付きプロファイルR因子 15%未満、10%以下が望ましい 最も一般的な収束指標。バックグラウンドを含む全データで計算
Rp プロファイルR因子 Rwpより低い値になる バックグラウンドの影響を受けやすい
RB ブラッグR因子 5%以下が理想 積分強度の一致度を評価。単結晶解析のR因子に相当
S(GooF) 適合度 1.0〜1.3が理想 Rwpを期待プロファイルR因子Reで割った値。1に近いほど良好


Rwpが10%以下になると、結晶相の存在割合(定量分析)が安定した結果として得られる傾向があります。三重県工業研究所の研究(酸化チタンを使ったリートベルト法精度向上に関する報告)でも、Rwpが10%を超えると結晶相割合の誤差が1割程度生じることが確認されています。


収束しない場合の原因として最も多いのは、INSファイルの入力ミスです。特に空間群の指定(HKLM1)の誤りと、格子定数の初期値のずれが主な原因です。次に多いのが「プロファイル関数の半値幅パラメータ(U, V, W)の設定」の問題です。収束しないときはFWHMのW値を増やすことが有効なことがあります。


R因子が改善されても、出力ファイル(.lstファイル)に記載された格子定数・席占有率・原子変位パラメータ(温度因子)の値が化学的に妥当かどうかを必ず確認してください。たとえば原子変位パラメータBが負の値になっている場合は、モデルに問題がある可能性があります。また、ゼロ点シフトパラメータと格子定数を同時に精密化すると、両者の相関で誤差が生じやすいため、段階的に固定・解放を行うのが現場での鉄則です。これは注意が必要です。


リートベルト解析のやり方・金属加工現場での実践的な活用と注意点

リートベルト解析は研究室だけのツールではありません。金属加工の製造現場でも、材料の品質管理・プロセス最適化・不良原因の特定に活用される場面が増えています。


製造プロセス管理への応用として代表的なのが、熱処理後の相変態の定量化です。たとえばステンレス系の材料を加熱処理した際に、オーステナイト相・フェライト相・マルテンサイト相がそれぞれ何%存在するかをリートベルト解析で定量できます。単純にXRDパターンを見るだけでは判断が難しい多相の比率を、±1%以内の精度で数値化できることが最大の強みです。


原料調達・不純物検査にも有効です。仕入れた金属粉末に意図しない結晶相が混入していないか、あるいは想定外の酸化物が生成していないかを、リートベルト解析によって定量的に確認できます。「似たような回折ピークがある」だけでは判断できない微量成分(たとえば質量比で5%程度の副相)の検出と定量が可能です。


金属加工現場でのリートベルト解析の限界と注意点も理解しておく必要があります。以下の4条件に当てはまるデータはリートベルト解析が使えません。


- 選択配向の影響が強く出ているデータ(プレス成形品の表面など)
- 粗大粒子の影響が強く出ているデータ(粉砕不足の試料)
- 入射X線が試料からはみ出しているデータ(通常の試料ホルダー充填なら問題なし)
- 試料が薄く、X線が試料ホルダーや基板まで突き抜けているデータ


つまり解析できないケースもあります。


また、リートベルト解析は「既存の結晶構造モデルがある物質」にしか適用できない点も重要です。まったく未知の結晶構造の決定には、別途「未知構造解析」が必要になります。疑問に思ったら専門の分析機関への相談も選択肢に入れると時間の節約になります。JFEテクノリサーチや材料技術研究所(MST)などは、リートベルト解析の受託サービスも提供しており、解析に慣れていない段階では外部委託を活用しながらノウハウを蓄積するやり方も現実的です。


解析スキルを体系的に学ぶには、書籍『RIETAN-FPで学ぶリートベルト解析』(坪田雅己・伊藤孝憲 著、情報機構)が初学者向けにRIETAN-FPのインストールからエラー対処法まで丁寧に解説しており、現場で使える一冊として高く評価されています。


リートベルト解析マスターが勧める 新規機能材料の開発に必須のICSD活用法(JAICI)- 実務でのCIFファイル活用・解析スキル習得の実体験


リートベルト解析・パターンを「眺めるだけ」でなく「数値化」する独自視点の解析戦略

多くの金属加工現場では、XRDパターンを取得しても「ピーク位置で相を同定して終わり」という運用が主流です。しかし、この「眺めるだけ」では回折パターンから引き出せる情報の1割も活用できていません。これが実態です。


リートベルト解析が強力なのは、格子定数の精密化により「結晶格子の微小なひずみ」まで検出できるからです。たとえば金属の加工硬化によって結晶格子が変形すると、Åの単位(1Å=0.0000001mm)で格子定数がわずかに変化します。この変化をリートベルト解析で数値化することで、加工前・加工後の材料変化を定量的に比較検討できます。


さらに一歩進んだ活用として、Williamson-Hallプロットがあります。これはピーク幅の回折角依存性から「結晶子サイズ」と「格子歪み」を分離して評価する手法で、RIETAN-FPのv.2.6以降に実装されています。結晶子サイズとは、X線が「単一の結晶」とみなせる領域の大きさで、研磨・プレス・熱処理によってどのように変化するかを評価できます。感覚で「硬くなった」「軟らかくなった」と判断していたものを、数nm〜数百nmスケールの数値で表現できるわけです。これは使えそうです。


さらに高度な活用例として、標準試料を少量(既知量)混合することで、非晶質成分の割合も定量できる点があります。金属製品の表面酸化膜や熱処理後の部分アモルファス化を定量評価したい場面では、この手法が力を発揮します。


解析結果の信頼性を上げるためのチェックポイントとして、最終的に結合長や結合角がVESTAで確認できる化学的に妥当な範囲に収まっているか、原子変位パラメータBが正の値をとっているかを確認する習慣をつけると良いです。結果の見直しが品質の分岐点です。


結晶子サイズとリートベルト解析の関係(ナノ分析テクノロジー)- ピーク重畳試料での定量と結晶子評価の解説






自転車 空気入れ 電動 電動空気入れ バイク 車 用 コードレス 6000mAh 大容量バッテリー 仏/英/米式バルブ対応 自動停止 低騒音 エアポンプ