- ベストアンサー
X線解析で
試料をX線回折で解析したとき、Kα1とKα2二つの波の重ねあわせでピークが出ますが、この時のKα2を除去する方法(計算法)を教えてください。 新しいX線回折装置を購入したのですが、その計算ソフトの信頼性を検討するため手計算することを考えています。よろしくお願いします。
- titanuser
- お礼率40% (18/44)
- その他(学問・教育)
- 回答数7
- ありがとう数1
- みんなの回答 (7)
- 専門家の回答
質問者が選んだベストアンサー
> このやりかただとKα1とKα2のピークが同じ位置に出ないんですか? 質問の意味が不明です。ブラッグの式はご存知ないのでしょうか? 横に引き伸ばす式は以下のとおりです。 θ2 = Arcsin [ (λ2 / λ1) ( sinθ1 ) ] 縦軸に関してはそのまま強度比を使えばいいはずです。少しは手を動かしながら考えましょう。
その他の回答 (6)
>分布関数が何になるのか、正規分布でいいんですかね? 覚えていません。 >とありますが、中心は単純にピークトップでいいのですか? 「ピークトップ」の意味を知らないので回答不能。 >計算方法はできればエクセルでできたら助かります。 当時.エクセルは販売されておらず.マルチプラン互換ソフトを使用していました。このソフトの製造・販売、対応機器の製造・販売が停止され10年以上が経ちました。よって使い物になりません。
- 38endoh
- ベストアンサー率53% (264/494)
先ほど帰りの電車の中で,もっとエレガントな算出方法がひらめきました。早速答えを書くと次のようになります。 1.Kα1 と Kα2 が混ざった実測の XRD を A とする。 2.これから求めたい,Kα1 のみからなる B とする。 3.XRD を,Kα1 と Kα2 の波長の差を元に横に引き伸ばし,強度比を元に縦に潰す操作を [X} とする。 このとき, B = A - [X] A' + [X]^2 A' - [X]^3 A' + … が答えです。つまり,実測の XRD から一回変換したものを差し引き,次に二回変換したものを足し,と繰り返せばよいのです。実際に図示すれば,式の意味が直感的にも分かるのではないでしょうか? なお,A' は A からベースラインを差し引き,平滑化したものとします。理想的な XRD では A = A' ですが,実際はそうはいきません。 この級数をどこまで展開するかは,実測の XRD の P/N 比(ピーク・ノイズ比,XRD ではこういいますよね)から判断すると良いと思います。恐らく,ノイズが多ければ第二項までで十分,P/N 比が良くても第三項までとれば十分でしょう。このような数値計算では多く展開すれば,それだけ数値計算に伴う丸め誤差などが生じ,かえって精度が落ちてしまいます。 って,意味分かります?
補足
またまた回答ありがとうございます。 >>XRD を,Kα1 と Kα2 の波長の差を元に横に引き伸ばし,強度比を元に縦に潰す操作を [X} とする。 ですが、このやりかただとKα1とKα2のピークが同じ位置に出ないんですか? 横に引き伸ばすのは移動する要素が入ってくるんですかね? 数学弱くてすみません。
- 38endoh
- ベストアンサー率53% (264/494)
1.Kα1 と Kα2 が混ざった実測の XRD を A とする。 2.これから求めたい,Kα1 のみの寄与を B とする。 3.何もないところから B を作るのは大変なので,初期値の B として A を用いる。この B は Kα1 のみからできていると思い込む。 4.B を横に引き伸ばし,縦に潰し,Kα2 の場合をシミュレートする。これを C とする。 5.B と C を足す。ただし,横に引き伸ばされた C は角度の間隔が広いため,C を平滑化したりラグランジュ補間したりして B と合う角度の部分を内挿によって求める。そして B と C が足されたものを D とする。この D がもし A と完全に一致していれば,B は Kα1 のみからできていると考えられる。 6.各々の角度について,A と D の差を求め,これを二乗し,そして足す(二乗平方和)。 7.この二乗平方和が D と A との乖離を示すので,これが最小値になるよう,B を微調整する。 8.新しい B を用いて,4.からもう一度やりなおす。二乗平方和が最小値になったら終わり。 初期値として Kα1 と Kα2 が混ざったものを使うのは,解がそれから大きく外れないからです。なお,7.における B の微調整を効率よく行うのはとても難しいです。…考えてください。
- 38endoh
- ベストアンサー率53% (264/494)
edogawaranpo さんの方法は,本質的には No.1 における私の方法と同等だと思います。 Kα2 成分の計算は,あまり深く考えなくても,Kα1 と Kα2 の波長の差を考えて横方向に引き伸ばし,強度比を考えて縦方向に圧縮するだけでよいと思います。ただ測定結果は離散的なデータですので,2 つの XRD の和を取るには,平滑化して平均してみるなどの何らかの補完処理は必要でしょう。 > 計算方法はできればエクセルでできたら助かります。 手入力で trial and error でもいいのですが,Excel には「ソルバーアドイン」なるものがありまして,これを用いれば最適値を簡単に計算することができます。残差二乗和のセルを用意すれば非線形最小二乗法を行うこともできます(参考 URL を見てください)。 もし余力がありましたら,Excel には VBA というプログラミング言語が内蔵されておりますので,この VBA からソルバーアドインを呼び出せば,所望の処理であれば十数行程度のプログラミングできますよ。
補足
2度も回答ありがとうございます。 >>その(実測の)XRD が Kα1 のみからなると仮定し,Kα2 のみの場合の XRD をシミュレートする… ということですが、その場合Kα2の分布はどうなるのでしょう?実測は両方の重ね合わせで成り立っていてKα1のピークは実測の低角度側に、Kα2は高角度側にくると思うのですが。 ちょっと、いやかなり混乱してきました-。-)
ASTMカード(今の名称はどなたかが回答しているので見つけてください)を見て.2つの波長と強度を入手します。 (名称忘却)分布関数にしたがって.中心から左右に別れたピークが観察されますから.強度にあわせた分布関数を計算します(岩波数学辞典.後ろの付録に計算式が載っています)。 2つのピークの和が測定値と重なるように.強度(と横方向への広がり方)を適当に変化させて計算します。 これは.試しに1%程度ずらして一致したら....と何回か計算してください(手早く計算するのであれば.非線型最小二乗法のヤコビ法なりなんなりを適当に選択して.もっとも.ソフトを造る手間を考えたらば.適当な数値を入れて「再計算」が楽です)。エクセル等表計算ソフトをつかえば簡単だと思います。 そのうち.ある程度のところで発散するなり振動しますから.ここで打ち切ります。 「ある程度のところ」は.計数したX線の数から計算できますね。 以上20年前にやってはいるのですが.すべて.忘却。いまさら.できません。
補足
回答ありがとうございます。 2つの波長と強度比は分かってます。 分布関数が何になるのか、正規分布でいいんですかね? あと >>中心から左右に別れたピーク... とありますが、中心は単純にピークトップでいいのですか? 計算方法はできればエクセルでできたら助かります。 申し訳ありませんがどなたか補足お願いします。
- 38endoh
- ベストアンサー率53% (264/494)
最も単純な方法として,以下のアルゴリズムでは如何でしょうか? 計算量は多くなりますが,確実に解は求まります。 1.初期値として Kα2 を含んだ実測の XRD を用いる。 2.その XRD が Kα1 のみからなると仮定し,Kα2 のみの場合の XRD をシミュレートする。 3.初期値に,2.でシミュレートした結果を足す。 4.実測値と3.との残差平方和を計算する。 5.残差平方和が最小となるよう,各角度における強度を求める。 6.残差平方和が一定値に収束するまで,2~5を繰り返す。 Kα1 と Kα2 の強度比,両者の波長については文献等で調べてください。 ただ,もっとちゃんと数学的な方法があると思いますので,他の方の回答を待ったほうが良いと思います(私も興味があります)。
関連するQ&A
- 流動解析ソフトの状況
何年か前にプラ型の流動解析ソフトを導入しようと考え調査しましたが、信頼性がまだまだ不十分と言うことと、高価であることから断念しました。 最近変形やガス・ウエルドといった問題が多く、再度検討したいと考えています。 最近の開発信頼状況は、どのようなレベルなのでしょうか? また、解析ソフトはどんなものがあるのでしょうか? お教えいただければ幸いです。
- 締切済み
- CAE
- X線を利用した回折実験についての質問です。
X線を利用した回折実験についての質問です。 以下のような回折角を得られた場合、外挿法を用いて格子定数を求める方法を教えてほしいです。 また、実験に用いた試料の結晶系は岩塩構造です。 お願いします。 hkl 特性X線 λ[nm] (400) (420) (422) (511) Cu-Kα 0.15418 94.17 109.90 127.44 144.04 Cu-Kβ1 0.13922 82.85 95.37 108.12 118.44
- 締切済み
- 物理学
- 解析用PCのスペックについて
FEM解析(有限要素法解析)用のPCの購入を検討しています。 解析用のため、スペックの高いものが必要だと思うのですが、下記項目を満たしたい場合、具体的にどのスペックが重要となるのでしょうか?(CPU?、メモリ?、他) ・解析時間をできるだけ短くしたい。 ・解析ソフトをサクサク操作したい。 (解析専用のため、他のソフトウェアは立ち上げません。) PCについてはあまり詳しくないので、初歩的な質問ですがよろしくお願いします。
- ベストアンサー
- 科学
- 応答曲面解析について
応答曲面解析は統計的手法の実験計画で直交表を用いない少ない実験回数でも交互作用を含めた2次曲面モデルが推定可能な解析法という事らしいのですが、業務に適用できるか検討したく、具体的な解析ソフトや専門書がありましたら、教えていただけませんか?
- 締切済み
- 数学・算数
- X線回折の管球の選択について
X線回折装置について質問いたします。 X線にCu-Kα線を使用している場合、Fe化合物の回折ができないと聞いたことがあるのですが、この理由をご存知なら教えてください。 使用管球の選択法もできれば教えてください。
- ベストアンサー
- 化学
- 台形公式を使って解析解を、。
2 ∫ X^5dx の定積分の解析解を求めるのに台形公式法を使って 0 求めよ。 なんですが、 全然回りの答えと合いません・・ 10.~ に出来るだけ近い数値が解です。 どの様にして求めればいいのでしょうか。 プログラムではなく手計算での求め方お願いします。
- 締切済み
- 数学・算数
- ある濃度範囲における検量線の解析について
例えば,0.2~50ng/mLの濃度範囲において,内部標準法にて8ポイント程度の 濃度を決めて(0.2,0.4,2,4,10,20,40,50ng/mL)検量線を作成したとします。 検量線を解析する際に, 原点を通る,通さない,含んで計算の三通りが解析ソフトで可能なのですが, それら三つの解析方法のメリットとデメリットは何があるのでしょうか? 理由として, (原点を通る式を使用することは,ほぼないのですが,) 通さないと含んでの計算の2パターンでサンプルの濃度にズレが出てきます。 どちらかというと,含んで計算が理論値に近いです。
- 締切済み
- 化学
- 真空雰囲気での熱流体解析について
半導体製造装置の設計しているものです。 現在、熱流体解析ソフトの導入し、チャンバ形状やヒータの熱解析等を検討しようと考えています。 ただ一点気になっているのは、半導体製造にはチャンバ内を真空状態で使用することが多くあるのですが、市販の熱流体解析ソフトは、雰囲気が粘性流域⇔中間流域⇔分子流域と変化しても対応しているものでしょうか? 熱流体解析ソフトは使ったことが無いので、何方か御教授下さい。 よろしくお願いします。
- 締切済み
- CAE
- 吸収スペクトルの波形分離解析ソフトについて
FT-IRの測定で得られたスペクトルに対してガンマ関数やローレンツ関数を用いての波形分離解析を考えています。 導入された測定装置がこれらの処理には未対応でしたが、ASCIIファイルでのデーター出力はできるようなので、外部の解析ソフトでの処理を検討しています。 ただ、予算的に高額な市販ソフトの購入が厳しそうなので、適当なフリーウェアでもないか探しているところです。ちゃんと山の分離ができれば、とりあえずの方向性検討には充分だと思っているのですが、どなたか適当なソフト、ご存知ないでしょうか。
- ベストアンサー
- 化学
お礼
うそです。 θ2ですね。
補足
>>少しは手を動かしながら考えましょう。 すみません、頑張ります。 >>θ2 = Arcsin [ (λ2 / λ1) ( sinθ1 ) ] 左辺はsinθ2ですか?