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:32]
koudai [例:CaGa4]
wien2k:電場勾配 [2020/03/31 20:42]
koudai [注意]
Line 3: Line 3:
   * WIEN2kではSCF計算の際に、NQR周波数を求めるのに使われる電場勾配 (EFG) の計算も同時に行ってくれます。   * WIEN2kではSCF計算の際に、NQR周波数を求めるのに使われる電場勾配 (EFG) の計算も同時に行ってくれます。
  
-====== 注意 ====== 
- 
-  * 収束は非常に遅いので<code> 
-$ grep ":EFG001:" case.scf</code>で必ず収束を確認してください。 
-  * 収束してないようでしたら、run_lapwでオプション-ccを使って、電荷の収束をより厳しい条件で再計算してください。 
-  * また、k点数に対しても収束が遅いので、必ずk点数を変化させて収束を確認してください 
-  * 周辺の原子位置に敏感な量ですので、構造最適化を行った後に計算することをおすすめします。 
  
 ====== 例:CaGa4 ====== ====== 例:CaGa4 ======
Line 54: Line 47:
 <code> <code>
 $ init_lapw -b -numk 1000 -rkmax 7.5 $ init_lapw -b -numk 1000 -rkmax 7.5
-$ run_lapw -cc 0.0001 -ec 0.00001 -i 100 -p+$ run_lapw -ec 0.00001 -cc 0.00001 -i 100 -p
 </code> </code>
  
Line 81: Line 74:
 :ETA001:                         ASYMM. ETA =     0.09099 :ETA001:                         ASYMM. ETA =     0.09099
 </file> </file>
-が1番目の原子のEFG計算の結果です。+が1番目(今の例ではCa)の原子のEFG計算の結果です。
  
 ===== 電場勾配 ===== ===== 電場勾配 =====
Line 139: Line 132:
  
 非対称パラメータ $\eta = \frac{V_{xx} - V_{yy}}{V_{zz}}$ の値です。 非対称パラメータ $\eta = \frac{V_{xx} - V_{yy}}{V_{zz}}$ の値です。
 +
 +
 +====== 注意 ======
 +
 +  * 収束は非常に遅いので<code>
 +$ grep ":EFG001:" case.scf</code>で必ず収束を確認してください。
 +  * 収束してないようでしたら、run_lapwでオプション-ccを使って、電荷の収束をより厳しい条件で再計算してください。
 +  * また、k点数に対しても収束が遅いので、必ずk点数を変化させて収束を確認してください
 +  * 周辺の原子位置に敏感な量ですので、実験と合わないときは構造最適化を行って計算してみるのもありです。
 +
wien2k/電場勾配.txt · Last modified: 2021/06/27 22:04 (external edit)