*＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊ * * FEFF 8.20 * feff.inp 見本ファイル * * 中山将伸 * * 書き間違い、解釈間違い等についてはなんら中山は責任を負えません。 * あしからず・・・・・（間違いやアドバイスなどは大歓迎です。） * masanobu@o.cc.titech.ac.jp====del-this-part====まで *＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊ *=========================================================================== * MAIN CONTROL (除: ATOMS) *=========================================================================== *タイトル TITLE Fd3m LiMn2O4 *モジュールコントロール************************************************* * pot xsph fms paths genfmt ff2chi CONTROL 1 1 1 1 1 1 PRINT 1 0 0 0 0 0 * pot 散乱ポテンシャルの計算(SCF or atomic overlap (Mattheiss) prescrpt) * xsph 位相シフト, 双極子行列, LDOS(副量子数別のDOS）散乱断面積、 * fms 多重散乱計算 * paths 散乱パスの計算 * genfmt 有効散乱振幅とその他のXAFSパラメータの計算 * ff2chi 得られた各パスの散乱信号のコンバイン * Print は p.9参照 * ff2chi=3でfeffNNNN.datファイル、 * xsph=1でldosNNNN.datファイル *格子の膨張収縮（１％アップなら、1.01)*********************************** * RMULTIPLIER 1.00 *ADVANCED OPTION***************** *CFAVERAGE(p.8) *OVERLAP(p.10) *=========================================================================== * POT: Scattering potentials *=========================================================================== *最初に孤立原子のポテンシャルを計算（Dirac−Fock atom code.) *求めた電子密度を重ねてマフィン-ティン-ポテンシャル近似(Mattheiss prescrip- *tion)またはSCF法で散乱ポテンシャル算出。平均自由行程はinterstial potential *の虚数部から求まる。 * AFOLP 1.0で、オーバーラップさせない ************ AFOLP * エッジ指定, (K, L1, L2, L3)またはNO (no hole) **************** EDGE K * S02が0.1以下で自動計算 ************************** S02 0.0 *POTENTIALS****************************************** POTENTIALS * ipot, Z, element: ipot=0 は吸収元素, Zは参照 elementは計算時参照しない * scmt セルフコンシステント fms 完全多重散乱（副量子数,副量子数 0:S,1:L,2:D) * stoichio: Fermiレベルの決定に使う 吸収元素は結晶性物質の場合0.01が推奨) * ipot Z element l_scmt l_fms stoichiometry 0 25 Mn 2 2 0.01 1 3 Li 1 1 8 2 25 Mn 2 2 16 3 8 O 1 1 32 *EXCHANGE******************************************** * ixc [ Vr Vi ix0] EXCHANGE 0 4 0.5 * Vr はFermiレベルの補正(約1eV) * Vi はpure imaginary "optical" potential * ixc0は、バックグラウンド関数 *NOHOLE********************************************** *NOOHOLE *コアホールがないと考えて計算する。バンド計算との対比に便利 *dDOSやL2,L3吸収ではむしろこちらのほうがお奨め。 *RGRID*********************************************** *RGRID *計算時のグリッド（動径部） デフォルト0.05 *場合によっては0.03,0.01がお奨め。特に高エネルギー側の場合 *SCF************************************************** * r_scf [ l_scf n_scf ca nmix] SCF 4.0 0 *** Radius of small cluster for self-consistency calculation *** r_scf クラスター半径（最低第２近接まで含める） *** l_scf 0:固体, 1:分子(def:0) *** n_scf scfでまわす数、10あれば十分(def:30) *** ca 収束加速因子(def:0.2), ρnext=ca*ρnew + (1-ca)ρold (f原子の場合0.05) *** nmix n_scfのサブでまわす数、(f原子の場合caと同時に10で指定） *UNFREEZEF******************************************** *UNFREEZEF *FEFF8.2では、f電子は最初のfree原子のそれをデフォルトで代用する *SCF計算したい場合はこのオプションが必須 *ADVANCED OPTIONS****** *INTERSTITIAL *ION *SPIN *=========================================================================== * XSPH: Cross-section and phase shifts *=========================================================================== *ELLIPTICITY (POLARIZATIONとともに用いる) *POLARIZATION(単結晶＆偏向ビーム） *MULTIPOLE le2 [l2lp]******************************************** *MULTIPOLE *双極子のときはle2=0, 双極子＋磁極子le2=1, 双極子＋四重極子のときはle2=2 *l2lpはL->L+1, L->L-1などの寄与が分かる *LDOS emin emax eimag******************************************** * emin emax eimag LDOS -30 20 0.1 * eimag ローレンツ関数でどの程度ブロードにさせるか（単位eV) *** Energy grid over which to calculate DOS functions (fix =84) *EXAFS xkmax ************************************************ *EXAFS 20 *デフォルト20A-1 *XANES [xkmax xkstep estep]************************************ XANES 4.0 *デフォルトでは計算されない！ *絶対エネルギーは、小さいZで数eV, 大きなZでは数百eVに及ぶ *** Upper limit of XANES calculation. *** This *must* be uncommented to *** make Feff calculate full multiple *** scattering rather than a path expansion *ADVANCED OPTIONS****** *DANES *FPRIME *XES *XMCD, XNCD *RPHASE (虚数フェーズシフト使わない） *RSIGMA (虚数エネルギーを使わない)＝＞XANES計算では過剰に見積もる傾向ある *=========================================================================== * FMS: Full multiple scattering *=========================================================================== *FMS: fMS計算をする******************************** * r_fms l_fms [minv toler1 toler2 rdirec] FMS 4.30925 0 * r_fmsはfMS計算をする球範囲（SCF計算ではない！） * rfmsはPOT以外のモジュール（LDOS, FMSなど）で用いるクラスター半径 * クラスター中の原子数は87以下がデフォルト(変えるにはコンパイル前のファイル * を書き換える必要がある)87以上の原子がある場合には、適切なクラスター半径 * になるようにr_fmsの値を自動的に下げる * EXAFSの場合は、FMS精度が高エネルギー領域で落ちるので、FMSとXANESオプ * ションは使わないほうがいい。ただし、SCFは用いたほうがいい。 * l_fms 0は固体、1は分子、（デフォルト0) * minv FMS計算のアルゴリズム選択。minv=2は試してみる価値あり * のこりはアルゴリズムを制御するパラメーター（よく分からない） * RPATHについては、rfms値より小さくすることを推奨する。MSの多重度増加は * XANES領域の計算発散を招きやすい。この場合はNLEGで2回散乱だけを使用した * ほうがいい *DEBYE temp thetad [idwopt] : 熱振動を考慮する ************************ * DEBYE * temp デバイ温度 *=========================================================================== * PATHS: Path enumeration *=========================================================================== *散乱パスで重要な寄与のみをフィルターリングし、パス長を長くする *ただし、FMS領域にあるクラスター以遠の領域を対象としている。 *RPATH rpath:上述の計算をする最大クラスター半径********************* RPATH 0.1 *XANESの場合はFMS計算のみで、high-order-MS計算を加えないほうがいい *EXAFSの場合は6が目安 *デフォルトは、第一近接×2.2となる *ADVANCED OPTION************* *PCRITERIA *SS (OVERLAPと使う, 結晶構造が分からないとき) *=========================================================================== * GENFMT: XAFS parameters *=========================================================================== * 散乱強度と位相フェーズなどの実際の計算ルーチン * FMS指定している範囲はCRITERIA, NLEGともに意味はない * CRITERIA critcw critpw ***************************** CRITERIA 0 0 * 重要でない（振幅の小さい）散乱の寄与を無視する * critcw は球面波の計算のcutoff(%) デフォルト4% cwが4%いかで計算する * critpw は平面波の場合　デフォルトは2.5% *NLEG nleg ******************************************* * NLEG 8 * 多重散乱の散乱回数のリミット（デフォルト８） *ADVANCED OPTION************* *NSTAR *=========================================================================== * FF2CHI: XAFS spectum *=========================================================================== *得られた散乱スペクトルの結合を行う *DEBYE temperature Debye-temperature [idwopt]********************** *DWの計算[単位K] *idwopt は計算法の選択 *CORRECTIONS real-energy-shift imaginary-energy-shift************** * CORRECTIONS 3.0 1.0 *real partはE0補正のために、imaginary partはピークのブロード化のため *SIG2 sig2 ******************************************************** *グローバルDWパラメーター *SIG2 0.001 *=========================================================================== * 座標 *=========================================================================== ATOMS * this list contains 83 atoms * x y z ipot tag distance 0.00000 0.00000 0.00000 0 Mn1 0.00000 0 1.95288 0.10712 0.10712 3 O1 1.95875 1 0.10712 1.95288 0.10712 3 O1 1.95875 2 0.10712 0.10712 1.95288 3 O1 1.95875 3 -1.95288 -0.10712 -0.10712 3 O1 1.95875 4 -0.10712 -1.95288 -0.10712 3 O1 1.95875 5 -0.10712 -0.10712 -1.95288 3 O1 1.95875 6 2.06000 2.06000 0.00000 2 Mn1 2.91328 7 -2.06000 -2.06000 0.00000 2 Mn1 2.91328 8 2.06000 0.00000 2.06000 2 Mn1 2.91328 9 0.00000 2.06000 2.06000 2 Mn1 2.91328 10 -2.06000 0.00000 -2.06000 2 Mn1 2.91328 11 0.00000 -2.06000 -2.06000 2 Mn1 2.91328 12 1.95288 1.95288 1.95288 3 O1 3.38249 13 -1.95288 -1.95288 -1.95288 3 O1 3.38249 14 -3.09000 1.03000 1.03000 1 Li1 3.41612 15 1.03000 -3.09000 1.03000 1 Li1 3.41612 16 -1.03000 -1.03000 3.09000 1 Li1 3.41612 17 -1.03000 3.09000 -1.03000 1 Li1 3.41612 18 3.09000 -1.03000 -1.03000 1 Li1 3.41612 19 1.03000 1.03000 -3.09000 1 Li1 3.41612 20 -2.16712 -2.16712 1.95288 3 O1 3.63408 21 -1.95288 2.16712 2.16712 3 O1 3.63408 22 2.16712 -1.95288 2.16712 3 O1 3.63408 23 2.16712 2.16712 -1.95288 3 O1 3.63408 24 -2.16712 1.95288 -2.16712 3 O1 3.63408 25 1.95288 -2.16712 -2.16712 3 O1 3.63408 26 -4.01288 -2.16712 0.10712 3 O1 4.56192 27 -2.16712 -4.01288 0.10712 3 O1 4.56192 28 -0.10712 4.01288 2.16712 3 O1 4.56192 29 4.01288 -0.10712 2.16712 3 O1 4.56192 30 -0.10712 2.16712 4.01288 3 O1 4.56192 31 2.16712 -0.10712 4.01288 3 O1 4.56192 32 4.01288 2.16712 -0.10712 3 O1 4.56192 33 2.16712 4.01288 -0.10712 3 O1 4.56192 34 -4.01288 0.10712 -2.16712 3 O1 4.56192 35 0.10712 -4.01288 -2.16712 3 O1 4.56192 36 -2.16712 0.10712 -4.01288 3 O1 4.56192 37 0.10712 -2.16712 -4.01288 3 O1 4.56192 38 -2.16712 4.22712 0.10712 3 O1 4.75147 39 4.22712 -2.16712 0.10712 3 O1 4.75147 40 -4.22712 -0.10712 2.16712 3 O1 4.75147 41 -0.10712 -4.22712 2.16712 3 O1 4.75147 42 -2.16712 0.10712 4.22712 3 O1 4.75147 43 0.10712 -2.16712 4.22712 3 O1 4.75147 44 -4.22712 2.16712 -0.10712 3 O1 4.75147 45 2.16712 -4.22712 -0.10712 3 O1 4.75147 46 4.22712 0.10712 -2.16712 3 O1 4.75147 47 0.10712 4.22712 -2.16712 3 O1 4.75147 48 -0.10712 2.16712 -4.22712 3 O1 4.75147 49 2.16712 -0.10712 -4.22712 3 O1 4.75147 50 -2.06000 4.12000 2.06000 2 Mn1 5.04595 51 4.12000 -2.06000 2.06000 2 Mn1 5.04595 52 -4.12000 -2.06000 2.06000 2 Mn1 5.04595 53 -2.06000 -4.12000 2.06000 2 Mn1 5.04595 54 -2.06000 2.06000 4.12000 2 Mn1 5.04595 55 2.06000 -2.06000 4.12000 2 Mn1 5.04595 56 4.12000 2.06000 -2.06000 2 Mn1 5.04595 57 -4.12000 2.06000 -2.06000 2 Mn1 5.04595 58 2.06000 4.12000 -2.06000 2 Mn1 5.04595 59 2.06000 -4.12000 -2.06000 2 Mn1 5.04595 60 -2.06000 2.06000 -4.12000 2 Mn1 5.04595 61 2.06000 -2.06000 -4.12000 2 Mn1 5.04595 62 5.15000 1.03000 1.03000 1 Li1 5.35204 63 1.03000 5.15000 1.03000 1 Li1 5.35204 64 3.09000 3.09000 3.09000 1 Li1 5.35204 65 1.03000 1.03000 5.15000 1 Li1 5.35204 66 -5.15000 -1.03000 -1.03000 1 Li1 5.35204 67 -1.03000 -5.15000 -1.03000 1 Li1 5.35204 68 -3.09000 -3.09000 -3.09000 1 Li1 5.35204 69 -1.03000 -1.03000 -5.15000 1 Li1 5.35204 70 4.12000 4.12000 0.00000 2 Mn1 5.82656 71 -4.12000 4.12000 0.00000 2 Mn1 5.82656 72 4.12000 -4.12000 0.00000 2 Mn1 5.82656 73 -4.12000 -4.12000 0.00000 2 Mn1 5.82656 74 4.12000 0.00000 4.12000 2 Mn1 5.82656 75 -4.12000 0.00000 4.12000 2 Mn1 5.82656 76 0.00000 4.12000 4.12000 2 Mn1 5.82656 77 0.00000 -4.12000 4.12000 2 Mn1 5.82656 78 4.12000 0.00000 -4.12000 2 Mn1 5.82656 79 -4.12000 0.00000 -4.12000 2 Mn1 5.82656 80 0.00000 4.12000 -4.12000 2 Mn1 5.82656 81 0.00000 -4.12000 -4.12000 2 Mn1 5.82656 82 END