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/04/04 16:59]
koudai [例:α-CaAlF5]
wien2k:電場勾配 [2020/07/05 21:34]
koudai [注意]
Line 66: Line 66:
 ===== SCF計算 ===== ===== SCF計算 =====
  
-  - 収束を厳し目にして、SCF計算を実行します<code>+  - SCF計算を実行します<code>
 $ init_lapw -b -numk 100 -rkmax 7.0 $ init_lapw -b -numk 100 -rkmax 7.0
-$ run_lapw -ec 0.000001 -cc 0.000001 -i 100 -p+$ run_lapw -ec 0.00001 -cc 0.0001 -i 100 -p
 </code> </code>
 ===== 計算結果 ===== ===== 計算結果 =====
Line 201: Line 201:
   * 電気四重極モーメント一覧 https://doi.org/10.1016/j.adt.2005.04.001   * 電気四重極モーメント一覧 https://doi.org/10.1016/j.adt.2005.04.001
   * 例えば27Alですと $Q=0.15~\mathrm{[10^{-28}m^2]}$ 程度です   * 例えば27Alですと $Q=0.15~\mathrm{[10^{-28}m^2]}$ 程度です
-  * bはバーンという単位で、$1~\mathrm{[b]} = 10^{-28}~\mathrm{[m^2]}$です+  * bはバーンという単位で、$1~\mathrm{[b]} = 10^{-28}~\mathrm{[m^2]}$です。散乱断面積に使われる単位ですが、NQRでも同じくらいのオーダーの量を扱うために借用されることが多いです。 
 +  * 四重極結合定数からNQR周波数を求める方法は核スピンの大きさで変わってくるので、教科書で確認してください。
  
  
Line 248: Line 249:
 $ grep ":EFG001:" case.scf</code>で確認できます。収束が十分でないようでしたら、収束条件をより厳しくして再計算してください。 $ grep ":EFG001:" case.scf</code>で確認できます。収束が十分でないようでしたら、収束条件をより厳しくして再計算してください。
   * k点数に対しても収束が遅いので、必ずk点数に関する収束も確認してください。特に金属の場合はかなり多くのk点をとらないと収束しません。   * k点数に対しても収束が遅いので、必ずk点数に関する収束も確認してください。特に金属の場合はかなり多くのk点をとらないと収束しません。
-  * 周辺の原子位置に敏感な量ですので、実験と合わないときは構造最適化を行っから計算してみるとよいでしょう+  * 周辺の原子位置に敏感な量ですので、基本的に内部座標の構造最適化と組み合わせて計算します。run_lapw -p -fc 0.1 -min ...
  
wien2k/電場勾配.txt · Last modified: 2021/06/27 22:04 (external edit)