ピーク関数のパラメーターなど

1. ：RIETAN-2000でHERMESデータを解析する場合、ピーク形状関数としては、NPRFN=2を選択と比較的一致がよいようですし、また扱いやすいようです。
   
   
   
2. ：RIETAN-97b、あるいはRIETAN-2000でプロファイル緩和を用いずに解析を行う場合、以下のパラメーターは装置に強く依存します。
   
   

   #Gaussian profile parameters, U, V, W, and P
   # Lorentzian profile parameters, X, Xe, Y, and Ye
   # Asymmetry parameter, As
   

   このうち、P，Xe，Yeはふつう0に固定しておきます。XeとYeはプロファイルが異方的に拡がっている場合に精密化しますが、ほとんどの場合、動かす必要はありません。
   
   
   
   
3. ：解析を行う時、U, V, W,X,YおよびAsなどの分解能パラメーターは、はじめはサンプルファイルの値で固定しておき、他のパラメーターでfittingする事をおすすめします。充分あってきたところで、これらのパラメーターをひとつづつ回し始めるのが良いと思います。
   なぜなら、これらのパラメーターは互いに強く相関があるので、いきなりいっぺんにまわし始めると、爆発したり、変な極小値に落ち込んでしまう可能性があるからです。
   
   
   
4. ：うごかす順番は、U, UV, UVW、ここで安定したらこの後はUVWをいっぺんに動かしても大丈夫です。UVWを動かしながら、X,YおよびAsはひとつつつ注意して動かします。つまり、Xだけ、Yだけ、Asだけを動かして安定して収束する事を確認してから、XYを同時にうごかし、さらにXYAsと動かします。
   
   
   
5. ：計算が正常にいっているならば、UとWは正、Vは負になります。XとYは負の値に収束する事もありますが、負の値を初期値とするとエラーを起こす場合がときどきあります。ちゃんと動いていたはずなのに、突然計算がスタートしなくなった場合は、X,Yの符号を確認し、負の場合は、とりあえず符号を正にしてまわします。ちゃんと動くはずです。
   
   
   
6. ：標準試料のデータを用いて得られたこれらのパラメーターはこちら-->HERMESパラメーター
7. ：通常は最小二乗法として修正Marquardt法を使用するので、NLESQ = 0とします。この場合、各Cycleの値をチェックしてください。正常に終了しているようにみえても、時として、
   

   Too large Marquardt parameter

   と表示されている場合があります。このときはパラメーターが変な値に落ちている可能性があります。動かすパラメーターを減らしてみましょう。どうしてもこのメッセージが消えないときは共役方向法(NLESQ = 2)を使うとうまくいくはずです。ただし、共役方向法は収束が遅いので何度もまわす必要があります。
   
   
   

プロファイル緩和

1. ：プロファイル緩和を使っても必ず一致が向上するわけではありません。プロファイル緩和を使わずによく一致するなら、無理に使う必要はありません。プロファイル緩和が有効なのは、プロファイル緩和なしで積分強度がほぼ説明できていても、形状があっていないピークがある場合です。このときは、プロファイル緩和をトライしてみます。
   
   
   
2. ：RIETAN-2000でプルファイル緩和を用いる場合、はじめはプロファイル緩和を用いずに全体の分解能パラメーターでfittingしておいて、おおよそあってきてから、単独ピークを１本つつえらんで、プロファイル緩和を行っていきます。いきなりたくさんのピークを加えると、爆発するかもしれません。プロファイル緩和を用いるとき、ひとつのピークには４，５個のパラメーターがありますが、これは一度に動かしても大丈夫なようです。
   
   
   
3. ：プロファイル緩和を行う場合、単独ピークのうちでプロファイルの不一致が一番目立つものからはじめます。もともと良くあっているピークにプロファイル緩和を用いても、S値は向上しないでしょう。
   
   
   
4. ：プロファイル緩和を用いる時は、他のピークが重なっていない単独のピークのみを選びます。複数のピークが重なって一本に見えているピークを使っては、プロファイル緩和の意味がないし、計算が爆発するかもしれません。
   単独ピークかどうかは、Igorグラフ中のブラッグ位置のマークの本数をチェックするか、ブラッグ指数とブラッグ角、強度の表を使うとよいでしょう。NPRINT = 1として 標準的なプリンター出力にすると、計算の最後にブラッグ指数とブラッグ角、強度の表が出力されます。（一度表を作ってしまえば、それ以降はNPRINT=0としたほうがよいでしょう。）
   
   
   
5. ：プロファイル緩和を用いるかどうかを決めるパラメーターはありません。プロファイル緩和を行うピークを指定するだけで、選ばれたピークのみでプロファイル緩和を行います。他の部分は通常のRIETANと同じ計算をします。
   
   
   
6. ：プロファイル緩和を行うピークの指定は、insファイル中の以下の部分で行います。
   

   # プロファイル緩和した反射において精密化するPPPは次の通り:
   # NPRFN = 1 (分割型pseudo-Voigt関数を適用): W, A, ηL, ηH.
   # NPRFN = 2 (拡張分割型pseudo-Voigt関数を適用): W1, W2, A, ηL, ηH.
   # NPRFN = 3 (分割型Pearson VII関数を適用): W, A, mL, mH.
   
   PPP1_2.1.0 0.450262 0.369632 1.13682 3.1259E-8 7.59959E-9 11111
   

   上の例では、2 1 0ブラッグピークに対し,ピーク形状関数の５つのパラメーターを全部動かしてプロファイル緩和を行う事になります。このPPP1_2.1.0 の赤い部分を書き換える事でピークの指定ができます。マイナスの指数も可能です。
   
   
   
7. ：プロファイル緩和を用いた単独ピークはIgorグラフ中ではブラッグ位置マークが紫に変わりますので、選択が正しいかグラフで確認しましょう。プロファイル緩和を行わなかったピークは通常どおり緑でマークされます。
   
   
   
8. ：さらに、１つピークを加えるたびにかならず計算結果のグラフをチェックして、そのピークでの一致が改善されてていくのを確認してください。なお、プロファイル緩和を行う際は全体の分解能パラメーター（U,V,Wなど）は固定しておいたほうがよいでしょう。プロファイル緩和によって安定して収束するのを確認してから、改めて分解能パラメーターも加えて計算するのが賢明です。
   
   
   
9. ：RIETAN-2000でも、計算が正常にいっているならばUとWは正、Vは負になります。（たぶん）
   標準試料のデータを用いて得られたこれらのパラメーターはこちら-->HERMESパラメーター
   
   
   
10. ：プルファイル緩和を用いれば多くのS値がかなり向上しますが、かならず向上するとは限りません。むしろSが悪くなる場合や、ピーク形状が変になる事すらあります。（例えばピークが分裂する場合もあります。）プルファイル緩和を用いたら、グラフ上の残差とピーク形状を必ず確認して、異常があればプロファイル緩和からははずします。

原子変位パラメーターBの扱い

# Label/species, g, x, y, z, B, and refinement identifiers (ID)
La/La 0.997093 0.0 0.0 0.0 0.344725 00001
B11/B11 1.0 0.361365 0.861365 0.5 0.425061 00001
C/C 0.989889 0.160048 0.660048 0.5 0.425061 00002
}
%
A(C,B)=A(B11,B)

一般的な注意点

# Peak-shift parameters Z, Ds, and Ts

# Label/species, g, x, y, z, B, and refinement identifiers (ID)
La/La 0.997093 0.0 0.0 0.0 0.344725 00001
B11/B11 1.0 0.361365 0.861365 0.5 0.425061 01201
C/C 0.989889 0.160048 0.660048 0.5 0.425061 00001
}
%
A(B11,y)=A(B11,x)+0.5