雑多な記録

User Tools

Site Tools

Trace:
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:44]
koudai [電場勾配] 		wien2k:電場勾配 [2020/04/02 19:46]
koudai [SCF計算]
Line 61: Line 61:
  
   * https://doi.org/10.1021/jp0740696   * https://doi.org/10.1021/jp0740696
 +  * この論文によると、構造最適化して計算したほうが良いということですが、今回はチュートリアルですので構造最適化前の結果を示します。
  
 ===== SCF計算 ===== ===== SCF計算 =====
  
-  - 内部構造最適化を行います<code> 
-$ init_lapw -b -numk 1000 -rkmax 7.0 
-$ run_lapw -ec 0.0001 -cc 0.001 -fc 0.5 -i 100 -p -min 
-</code> 
   - 収束を厳し目にして、SCF計算を再実行します<code>   - 収束を厳し目にして、SCF計算を再実行します<code>
-$ run_lapw -ec 0.000001 -cc 0.000001 -i 100 -p -NI+$ init_lapw -b -numk 100 -rkmax 7.0 
 +$ run_lapw -ec 0.000001 -cc 0.000001 -i 100 -p
 </code> </code>
- 
 ===== 計算結果 ===== ===== 計算結果 =====
  
Line 79: Line 76:
 ... ...
  
-:VKCOUL :  VK-COUL convergence: 0.306E-06 +:VKCOUL :  VK-COUL convergence: 0.307E-06 
-       ATOMNUMBER=  1 F          VCOUL-ZERO =  0.28719E+00 +       ATOMNUMBER=  1 F          VCOUL-ZERO =  0.28637E+00 
-:EFG001:                        EFG         =     6.97758   *10**21  V / m**2 +:EFG001:                        EFG         =     6.98229   *10**21  V / m**2 
-                               V20  TOT/SRF=    -3.15744    -0.61722 +                               V20  TOT/SRF=    -3.15977    -0.61904 
-                               V22  TOT/SRF=     5.07464     1.07082 +                               V22  TOT/SRF=     5.07794     1.07422 
-                               V22M TOT/SRF=     0.06173     0.00930 +                               V22M TOT/SRF=     0.06182     0.00939 
-                               V21  TOT/SRF=     0.92073     0.14769 +                               V21  TOT/SRF=     0.92145     0.14812 
-                               V21M TOT/SRF=    -0.20544    -0.01225+                               V21M TOT/SRF=    -0.20537    -0.01217
  
-              6.89759    0.06173    0.92073        6.97758    0.00000    0.00000 +              6.90223    0.06182    0.92145        6.98229    0.00000    0.00000 
-              0.06173   -3.25169   -0.20544        0.00000   -3.80538    0.00000 +              0.06182   -3.25365   -0.20537        0.00000   -3.80799    0.00000 
-              0.92073   -0.20544   -3.64590        0.00000    0.00000   -3.17220+              0.92145   -0.20537   -3.64858        0.00000    0.00000   -3.17431
  
-         MAIN DIRECTIONS OF THE EFG   1.0000 -0.0882  0.0285 +         MAIN DIRECTIONS OF THE EFG   1.0000 -0.0882  0.0284 
-                                      0.0043  0.3809  1.0000 +                                      0.0043  0.3803  1.0000 
-                                      0.0866  1.0000 -0.3784+                                      0.0866  1.0000 -0.3778
 :ANG001:  ANGLE WITH OLD X-AXIS =         5.3 :ANG001:  ANGLE WITH OLD X-AXIS =         5.3
  
-:ETA001:                         ASYMM. ETA =     0.09074+:ETA001:                         ASYMM. ETA =     0.09076
  
  
-       ATOMNUMBER=  2 F          VCOUL-ZERO =  0.24371E+00 +       ATOMNUMBER=  2 F          VCOUL-ZERO =  0.24339E+00 
-:EFG002:                        EFG         =     5.78076   *10**21  V / m**2 +:EFG002:                        EFG         =     5.78591   *10**21  V / m**2 
-                               V20  TOT/SRF=     4.93250     0.80157 +                               V20  TOT/SRF=     4.93700     0.80506 
-                               V22  TOT/SRF=    -2.35224    -0.30915 +                               V22  TOT/SRF=    -2.35143    -0.30831 
-                               V22M TOT/SRF=    -0.04393    -0.00921 +                               V22M TOT/SRF=    -0.04398    -0.00926 
-                               V21  TOT/SRF=     0.15501     0.04408 +                               V21  TOT/SRF=     0.15546     0.04446 
-                               V21M TOT/SRF=    -0.72116    -0.22754+                               V21M TOT/SRF=    -0.72140    -0.22776
  
-             -5.20002   -0.04393    0.15501       -5.20247    0.00000    0.00000 +             -5.20181   -0.04398    0.15546       -5.20427    0.00000    0.00000 
-             -0.04393   -0.49554   -0.72116        0.00000    5.78076    0.00000 +             -0.04398   -0.49895   -0.72140        0.00000    5.78591    0.00000 
-              0.15501   -0.72116    5.69556        0.00000    0.00000   -0.57829+              0.15546   -0.72140    5.70076        0.00000    0.00000   -0.58164
  
-         MAIN DIRECTIONS OF THE EFG   1.0000  0.0146 -0.0056 +         MAIN DIRECTIONS OF THE EFG   1.0000  0.0146 -0.0057 
-                                      0.0072 -0.1150  1.0000 +                                      0.0072 -0.1149  1.0000 
-                                     -0.0137  1.0000  0.1151+                                     -0.0138  1.0000  0.1150
 :ANG002:  ANGLE WITH OLD X-AXIS =         0.9 :ANG002:  ANGLE WITH OLD X-AXIS =         0.9
  
-:ETA002:                         ASYMM. ETA =     0.79993+:ETA002:                         ASYMM. ETA =     0.79894
  
  
-       ATOMNUMBER=  3 F          VCOUL-ZERO =  0.19796E+00 +       ATOMNUMBER=  3 F          VCOUL-ZERO =  0.19692E+00 
-:EFG003:                        EFG         =     8.98492   *10**21  V / m**2 +:EFG003:                        EFG         =     8.98942   *10**21  V / m**2 
-                               V20  TOT/SRF=    -1.67223    -0.39784 +                               V20  TOT/SRF=    -1.67493    -0.39934 
-                               V22  TOT/SRF=    -8.01943    -1.13060 +                               V22  TOT/SRF=    -8.02238    -1.13225 
-                               V22M TOT/SRF=    -0.01957     0.00012+                               V22M TOT/SRF=    -0.01953     0.00017
                                V21  TOT/SRF=     0.00000     0.00000                                V21  TOT/SRF=     0.00000     0.00000
                                V21M TOT/SRF=     0.00000     0.00000                                V21M TOT/SRF=     0.00000     0.00000
  
-             -7.05397   -0.01957    0.00000       -7.05400    0.00000    0.00000 +             -7.05536   -0.01953    0.00000       -7.05538    0.00000    0.00000 
-             -0.01957    8.98490    0.00000        0.00000    8.98492    0.00000 +             -0.01953    8.98940    0.00000        0.00000    8.98942    0.00000 
-              0.00000    0.00000   -1.93093        0.00000    0.00000   -1.93093+              0.00000    0.00000   -1.93404        0.00000    0.00000   -1.93404
  
          MAIN DIRECTIONS OF THE EFG   1.0000 -0.0012  0.0000          MAIN DIRECTIONS OF THE EFG   1.0000 -0.0012  0.0000
Line 137: Line 134:
 :ANG003:  ANGLE WITH OLD X-AXIS =         0.1 :ANG003:  ANGLE WITH OLD X-AXIS =         0.1
  
-:ETA003:                         ASYMM. ETA =     0.57019+:ETA003:                         ASYMM. ETA =     0.56971
  
  
-       ATOMNUMBER=  4 Ca         VCOUL-ZERO = -0.24919E+00 +       ATOMNUMBER=  4 Ca         VCOUL-ZERO = -0.24761E+00 
-:EFG004:                        EFG         =    -2.26265   *10**21  V / m**2 +:EFG004:                        EFG         =    -2.26246   *10**21  V / m**2 
-                               V20  TOT/SRF=     1.73744     0.03173 +                               V20  TOT/SRF=     1.73886     0.03329 
-                               V22  TOT/SRF=    -1.24948    -0.03004 +                               V22  TOT/SRF=    -1.24847    -0.03071 
-                               V22M TOT/SRF=     0.15884     0.00427+                               V22M TOT/SRF=     0.15886     0.00430
                                V21  TOT/SRF=     0.00000     0.00000                                V21  TOT/SRF=     0.00000     0.00000
                                V21M TOT/SRF=     0.00000     0.00000                                V21M TOT/SRF=     0.00000     0.00000
  
-             -2.25260    0.15884    0.00000       -2.26265    0.00000    0.00000 +             -2.25240    0.15886    0.00000       -2.26246    0.00000    0.00000 
-              0.15884    0.24637    0.00000        0.00000    0.25643    0.00000 +              0.15886    0.24454    0.00000        0.00000    0.25460    0.00000 
-              0.00000    0.00000    2.00623        0.00000    0.00000    2.00623+              0.00000    0.00000    2.00786        0.00000    0.00000    2.00786
  
-         MAIN DIRECTIONS OF THE EFG   1.0000  0.0633  0.0000 +         MAIN DIRECTIONS OF THE EFG   1.0000  0.0634  0.0000 
-                                     -0.0633  1.0000  0.0000+                                     -0.0634  1.0000  0.0000
                                       0.0000  0.0000  1.0000                                       0.0000  0.0000  1.0000
 :ANG004:  ANGLE WITH OLD X-AXIS =         3.6 :ANG004:  ANGLE WITH OLD X-AXIS =         3.6
  
-:ETA004:                         ASYMM. ETA =     0.77334+:ETA004:                         ASYMM. ETA =     0.77493
  
  
-       ATOMNUMBER=  5 Al         VCOUL-ZERO = -0.77946E+00 +       ATOMNUMBER=  5 Al         VCOUL-ZERO = -0.78227E+00 
-:EFG005:                        EFG         =     2.97056   *10**21  V / m**2 +:EFG005:                        EFG         =     2.96917   *10**21  V / m**2 
-                               V20  TOT/SRF=    -0.07974     0.00494 +                               V20  TOT/SRF=    -0.08033     0.00488 
-                               V22  TOT/SRF=    -2.90640    -0.12276 +                               V22  TOT/SRF=    -2.90471    -0.11829 
-                               V22M TOT/SRF=    -0.03759    -0.00159 +                               V22M TOT/SRF=    -0.03752    -0.00157 
-                               V21  TOT/SRF=     0.27118     0.01021 +                               V21  TOT/SRF=     0.27084     0.00970 
-                               V21M TOT/SRF=     0.23527     0.02798+                               V21M TOT/SRF=     0.23493     0.02708
  
-             -2.86036   -0.03759    0.27118       -2.88730    0.00000    0.00000 +             -2.85833   -0.03752    0.27084       -2.88523    0.00000    0.00000 
-             -0.03759    2.95244    0.23527        0.00000   -0.08326    0.00000 +             -0.03752    2.95110    0.23493        0.00000   -0.08394    0.00000 
-              0.27118    0.23527   -0.09208        0.00000    0.00000    2.97056+              0.27084    0.23493   -0.09276        0.00000    0.00000    2.96917
  
          MAIN DIRECTIONS OF THE EFG   1.0000  0.0987 -0.0029          MAIN DIRECTIONS OF THE EFG   1.0000  0.0987 -0.0029
-                                      0.0104 -0.0763  1.0000 +                                      0.0104 -0.0762  1.0000 
-                                     -0.0979  1.0000  0.0766+                                     -0.0979  1.0000  0.0765
 :ANG005:  ANGLE WITH OLD X-AXIS =         5.6 :ANG005:  ANGLE WITH OLD X-AXIS =         5.6
  
-:ETA005:                         ASYMM. ETA =     0.94394+:ETA005:                         ASYMM. ETA =     0.94346 
  
 ... ...
Line 201: Line 199:
 電気四重極モーメントは原子核の種類によって決まります。 電気四重極モーメントは原子核の種類によって決まります。
  
-  * 電気四重極モーメント一覧 https://www.sciencedirect.com/science/article/pii/S0092640X05000239 +  * 電気四重極モーメント一覧 https://doi.org/10.1016/j.adt.2005.04.001 
-  * 例えば27Alですと $Q=14.9~\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 248:
 $ 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)

Page Tools

Except where otherwise noted, content on this wiki is licensed under the following license: CC Attribution-Share Alike 4.0 International
CC Attribution-Share Alike 4.0 International Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki