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/02 13:42]
koudai
wien2k:電場勾配 [2020/04/02 15:03]
koudai [注意]
Line 64: Line 64:
 ===== SCF計算 ===== ===== SCF計算 =====
  
-  - 内部構造最適化を行います<code>+  - (任意)内部構造最適化を行います<code>
 $ init_lapw -b -numk 1000 -rkmax 7.0 $ init_lapw -b -numk 1000 -rkmax 7.0
 $ run_lapw -ec 0.0001 -cc 0.001 -fc 0.5 -i 100 -p -min $ run_lapw -ec 0.0001 -cc 0.001 -fc 0.5 -i 100 -p -min
Line 71: Line 71:
 $ run_lapw -ec 0.000001 -cc 0.000001 -i 100 -p -NI $ run_lapw -ec 0.000001 -cc 0.000001 -i 100 -p -NI
 </code> </code>
 +    * 構造最適化のフラグ-minをつけて厳しい収束条件にすると、いつまでたっても計算が終わらなくなります。
 ===== 計算結果 ===== ===== 計算結果 =====
  
Line 201: Line 201:
 電気四重極モーメントは原子核の種類によって決まります。 電気四重極モーメントは原子核の種類によって決まります。
  
-  * 電気四重極モーメント一覧 https://www.sciencedirect.com/science/article/pii/S0092640X05000239 +  * 電気四重極モーメント一覧 https://doi.org/10.1016/j.adt.2005.04.001 
-  * 例えば43Caですと $Q=4.3~\mathrm{[\times 10^{-28}m^2]}$ です+  * 例えば27Alですと $Q=0.15~\mathrm{[\times 10^{-28}m^2]}$ 程度です
   * bはバーンという単位で、$1~\mathrm{[b]} = 10^{-28}~\mathrm{[m^2]}$です   * bはバーンという単位で、$1~\mathrm{[b]} = 10^{-28}~\mathrm{[m^2]}$です
  
Line 250: Line 250:
 $ grep ":EFG001:" case.scf</code>で確認できます。収束が十分でないようでしたら、収束条件をより厳しくして再計算してください。 $ grep ":EFG001:" case.scf</code>で確認できます。収束が十分でないようでしたら、収束条件をより厳しくして再計算してください。
   * k点数に対しても収束が遅いので、必ずk点数に関する収束も確認してください。特に金属の場合はかなり多くのk点をとらないと収束しません。   * k点数に対しても収束が遅いので、必ずk点数に関する収束も確認してください。特に金属の場合はかなり多くのk点をとらないと収束しません。
-  * 周辺の原子位置に敏感な量ですので、実験と合わないときは構造最適化を行ってから計算してみるのも一つの手す。+  * 周辺の原子位置に敏感な量ですので、実験と合わないときは構造最適化を行ってから計算してみるとよいしょう
  
wien2k/電場勾配.txt · Last modified: 2021/06/27 22:04 (external edit)