User Tools

Site Tools


wien2k:電場勾配

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
Next revision Both sides next revision
wien2k:電場勾配 [2020/03/31 20:42]
koudai [注意]
wien2k:電場勾配 [2020/04/01 09:17]
koudai [SCF計算後]
Line 43: Line 43:
  
  
-電荷の収束を厳し目にして、通常のSCF計算を実行します+収束を厳し目にして、通常のSCF計算を実行します
  
 <code> <code>
-$ init_lapw -b -numk 1000 -rkmax 7.5 +$ init_lapw -b -numk 10000 -rkmax 7.5 
-$ run_lapw -ec 0.00001 -cc 0.00001 -i 100 -p+$ run_lapw -ec 0.000001 -cc 0.000001 -i 100 -p
 </code> </code>
  
-====== 電場勾配 ======+===== 計算結果 =====
  
 SCF計算が終了したら case.scf0 を開いてください。 SCF計算が終了したら case.scf0 を開いてください。
Line 76: Line 76:
 が1番目(今の例ではCa)の原子のEFG計算の結果です。 が1番目(今の例ではCa)の原子のEFG計算の結果です。
  
-===== 電場勾配 =====+==== 電場勾配 ====
  
   :EFG001:                        EFG             0.65471   *10**21  V / m**2   :EFG001:                        EFG             0.65471   *10**21  V / m**2
Line 97: Line 97:
  
  
-===== 電気四重極テンソル =====+==== 電気四重極テンソル ====
  
                 0.01496    0.44715    0.00000       -0.29757    0.00000    0.00000                 0.01496    0.44715    0.00000       -0.29757    0.00000    0.00000
Line 117: Line 117:
 右側が主軸変換(対角化)後の結果で、固有値の絶対値が大きい順に$|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   1.0000  0.6989  0.0000            MAIN DIRECTIONS OF THE EFG   1.0000  0.6989  0.0000
Line 127: Line 127:
  
  
-===== 非対称パラメータ =====+==== 非対称パラメータ ====
  
   :ETA001:                         ASYMM. ETA =     0.09099   :ETA001:                         ASYMM. ETA =     0.09099
Line 136: Line 136:
 ====== 注意 ====== ====== 注意 ======
  
-  * 収束は非常に遅いので<code> +  * SCF計算各ステップの電場勾配の値は<code> 
-$ grep ":EFG001:" case.scf</code>必ず収束を確認してください +$ grep ":EFG001:" case.scf</code>で確認できます。収束が十分でないようでしたら、収束条件をより厳しくして再計算してください。 
-  * 収束してないようでしたら、run_lapwでオプション-ccを使って、電荷の収束をより厳しい条件で再計算してください。 +  * k点数に対しても収束が遅いので、必ずk点数に関する収束確認してください 
-  * また、k点数に対しても収束が遅いので、必ずk点数を変化させて収束確認してください +  * 周辺の原子位置に敏感な量ですので、実験と合わないときは構造最適化を行ってから計算してみるのも一つの手です。
-  * 周辺の原子位置に敏感な量ですので、実験と合わないときは構造最適化を行って計算してみるのもありです。+
  
wien2k/電場勾配.txt · Last modified: 2021/06/27 22:04 (external edit)