> ./fp2sflip.py BaSO4.ffoこれでできるBaSO4.inflipファイルはまだ完全ではないが、とりあえずはそのままでSuperflipを実行できるようになっている。直す部分としては組成とZが分かっているなら、キーワードcompositionを正しく与えてやる必要がある。BaSO4の場合は、Z=4なので、"composition Ba4 S4 O16"となる。また先頭の#を削除する。さらにhistogram compositionの先頭の#を削除する。もし組成自体が不確定な場合はcompositionとhistogram行は#を先頭につけたままにする。このテンプレートは標準的なセッティングなので、必要に応じてパラメーターなどを変える必要があるが、自分の好みでコード自体を直せばよい。
> ./superflip BaSO4.inflipで計算を実行させることができる。計算が終了したらSuccessがあるかチェック。Successがゼロでも実際には正解の電子密度が得られていることもある。また、何度か繰り返して計算させることも必要。計算で出来たxplorファイルをダブルクリックして、Vestaで電子密度を表示させる。RIETAN-FP付属のBaSO4データを解析した場合、Ba, Sは明瞭に見えるが、酸素は最初はっきりしないが、VestaのObjects->Properties->Isosurfaces...の設定からIsosurface levelを下げてやるとSの周りに4個電子分布が見えるので、それらが酸素であろう。
> EDMA BaSO4.edmaで実行する。cifファイルが作られるので、それをVestaで読む。BaSO4の結果では酸素を含めて全ての元素が検出された。原子種の間違いなどはVesta上で編集すればよい。RIETAN-FP付属のfluoro-apatiteの例ではCa, P, Fは明瞭に見えるが、酸素位置は求められなかった。求まったCa, P, Fの位置を使って、RIETAN-FPで部分構造を入れたLe Bail法を行えばいいはずだが、今のところうまくいっていない。