This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision Next revision Both sides next revision | ||
wien2k:電場勾配 [2020/03/28 15:41] koudai [電場勾配] |
wien2k:電場勾配 [2020/04/01 14:03] koudai [注意] |
||
---|---|---|---|
Line 3: | Line 3: | ||
* WIEN2kではSCF計算の際に、NQR周波数を求めるのに使われる電場勾配 (EFG) の計算も同時に行ってくれます。 | * WIEN2kではSCF計算の際に、NQR周波数を求めるのに使われる電場勾配 (EFG) の計算も同時に行ってくれます。 | ||
- | ====== 注意 ====== | ||
- | * 収束は非常に遅いので<code> | + | ====== 例:CaGa4 ====== |
- | $ grep ":EFG001:" case.scf</code>で必ず収束を確認してください。 | + | |
- | * 収束してないようでしたら、run_lapwでオプション-ccを使って、電荷の収束をより厳しい条件で再計算してください。 | + | 以下の論文の結晶構造から作成しました。 |
- | * また、k点数に対しても収束が遅いので、必ずk点数を変化させて収束を確認してください | + | |
- | * 周辺の原子位置に敏感な量ですので、構造最適化を行った後に計算することをおすすめします。 | + | * https:// |
+ | |||
+ | <file - CaGa4.struct> | ||
+ | blebleble | ||
+ | CXZ LATTICE, | ||
+ | MODE OF CALC=RELA unit=bohr | ||
+ | | ||
+ | ATOM -1: X=0.00000000 Y=0.00000000 Z=0.00000000 | ||
+ | MULT= 1 ISPLIT=15 | ||
+ | Ca | ||
+ | LOCAL ROT MATRIX: | ||
+ | | ||
+ | | ||
+ | ATOM -2: X=0.40503000 Y=0.76923000 Z=0.00000000 | ||
+ | MULT= 2 ISPLIT=15 | ||
+ | -2: X=0.59497000 Y=0.23077000 Z=0.00000000 | ||
+ | Ga | ||
+ | LOCAL ROT MATRIX: | ||
+ | | ||
+ | | ||
+ | ATOM -3: X=0.00000000 Y=0.50000000 Z=0.23283000 | ||
+ | MULT= 2 ISPLIT=15 | ||
+ | -3: X=0.00000000 Y=0.50000000 Z=0.76717000 | ||
+ | Ga | ||
+ | LOCAL ROT MATRIX: | ||
+ | | ||
+ | | ||
+ | | ||
+ | </file> | ||
+ | |||
+ | この物質のWIEN2kによるNQRの研究は、以下の論文の5.1節が参考になります | ||
+ | |||
+ | * https:// | ||
+ | |||
+ | |||
+ | 収束を厳し目にして、通常のSCF計算を実行します | ||
+ | |||
+ | < | ||
+ | $ init_lapw -b -numk 10000 -rkmax 7.5 | ||
+ | $ run_lapw -ec 0.000001 -cc 0.000001 -i 100 -p | ||
+ | </ | ||
- | ====== 電場勾配 ====== | + | ===== 計算結果 |
- | case.scf0 を開いてください。 | + | SCF計算が終了したら CaGe4.scf0 を開いてください |
- | 以下はCaGa4の例です。 | + | |
<file - CaGe4.scf0> | <file - CaGe4.scf0> | ||
Line 35: | Line 73: | ||
: | : | ||
</ | </ | ||
- | が1番目の原子のEFG計算の結果です。 | + | が1番目(今の例ではCa)の原子のEFG計算の結果です。 |
- | ===== 電場勾配 | + | ==== 電場勾配 ==== |
: | : | ||
Line 58: | Line 96: | ||
- | ===== 電気四重極テンソル | + | ==== 電気四重極テンソル ==== |
0.01496 | 0.01496 | ||
Line 78: | Line 116: | ||
右側が主軸変換(対角化)後の結果で、固有値の絶対値が大きい順に$|V_{zz}| \geq |V_{yy}| \geq |V_{xx}|$と定義されますので、$V_{zz} = 0.65471$, $V_{yy} = -0.35714$, , $V_{xx} = -0.29757$です。 | 右側が主軸変換(対角化)後の結果で、固有値の絶対値が大きい順に$|V_{zz}| \geq |V_{yy}| \geq |V_{xx}|$と定義されますので、$V_{zz} = 0.65471$, $V_{yy} = -0.35714$, , $V_{xx} = -0.29757$です。 | ||
- | ===== 主軸回転 | + | ==== 主軸回転 ==== |
MAIN DIRECTIONS OF THE EFG | MAIN DIRECTIONS OF THE EFG | ||
Line 88: | Line 126: | ||
- | ===== 非対称パラメータ | + | ==== 非対称パラメータ ==== |
: | : | ||
非対称パラメータ $\eta = \frac{V_{xx} - V_{yy}}{V_{zz}}$ の値です。 | 非対称パラメータ $\eta = \frac{V_{xx} - V_{yy}}{V_{zz}}$ の値です。 | ||
+ | |||
+ | |||
+ | ====== 注意 ====== | ||
+ | |||
+ | * SCF計算の各ステップでの電場勾配の値は< | ||
+ | $ grep ": | ||
+ | * k点数に対しても収束が遅いので、必ずk点数に関する収束も確認してください。特に金属の場合はかなり多くのk点をとらないと収束しません。 | ||
+ | * 周辺の原子位置に敏感な量ですので、実験と合わないときは構造最適化を行ってから計算してみるのも一つの手です。 | ||
+ |