雑多な記録

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/01 21:08]
koudai 		wien2k:電場勾配 [2020/04/02 13:32]
koudai [例：CaGa4]
Line 4: Line 4:
  
  
-====== 例：CaGa4 ======+====== 例：α-CaAlF5 ======
  
 以下の論文の結晶構造から作成しました。 以下の論文の結晶構造から作成しました。
  
-  * https://www.sciencedirect.com/science/article/pii/0022508889901719+  * https://scripts.iucr.org/cgi-bin/paper?du0231
  
-<file - CaGa4.struct>+<file - alpa-CaAlF5.struct>
 blebleble                                                                        blebleble                                                                       
-CXZ LATTICE,NONEQUIV.ATOMS   3   12 C2/m    +CXZ LATTICE,NONEQUIV.ATOMS   5   15 C2/c    
 MODE OF CALC=RELA unit=bohr MODE OF CALC=RELA unit=bohr
- 11.680586 11.559833 11.584210 90.000000 90.000000118.942000 + 16.463294 13.887597 11.937400 90.000000 90.000000115.040000 
-ATOM  -1: X=0.00000000 Y=0.00000000 Z=0.00000000 +ATOM  -1: X=0.22110000 Y=0.61300000 Z=0.01780000 
-          MULT= 1          ISPLIT=15 +          MULT= 4          ISPLIT=15 
-Ca         NPT=  781  R0=.000050000 RMT= 2.50000     Z:  20.00000+      -1: X=0.77890000 Y=0.38700000 Z=0.98220000 
 +      -1: X=0.77890000 Y=0.88700000 Z=0.01780000 
 +      -1: X=0.22110000 Y=0.11300000 Z=0.98220000 
 +F          NPT=  781  R0=.000100000 RMT=   2.00000   Z:   9.00000
 LOCAL ROT MATRIX:    1.0000000 0.0000000 0.0000000 LOCAL ROT MATRIX:    1.0000000 0.0000000 0.0000000
                      0.0000000 1.0000000 0.0000000                      0.0000000 1.0000000 0.0000000
                      0.0000000 0.0000000 1.0000000                      0.0000000 0.0000000 1.0000000
-ATOM  -2: X=0.40503000 Y=0.76923000 Z=0.00000000+ATOM  -2: X=0.48860000 Y=0.02630000 Z=0.21690000 
 +          MULT= 4          ISPLIT=15 
 +      -2: X=0.51140000 Y=0.97370000 Z=0.78310000 
 +      -2: X=0.51140000 Y=0.47370000 Z=0.21690000 
 +      -2: X=0.48860000 Y=0.52630000 Z=0.78310000 
 +F          NPT=  781  R0=.000100000 RMT=   2.00000   Z:   9.00000 
 +LOCAL ROT MATRIX:    1.0000000 0.0000000 0.0000000 
 +                     0.0000000 1.0000000 0.0000000 
 +                     0.0000000 0.0000000 1.0000000 
 +ATOM  -3: X=0.00000000 Y=0.25000000 Z=0.05780000
           MULT= 2          ISPLIT=15           MULT= 2          ISPLIT=15
-      -2: X=0.59497000 Y=0.23077000 Z=0.00000000 +      -3: X=0.00000000 Y=0.75000000 Z=0.94220000 
-Ga         NPT=  781  R0=.000050000 RMT= 2.33        Z:  31.00000+F          NPT=  781  R0=.000100000 RMT=   2.00000   Z:   9.00000
 LOCAL ROT MATRIX:    1.0000000 0.0000000 0.0000000 LOCAL ROT MATRIX:    1.0000000 0.0000000 0.0000000
                      0.0000000 1.0000000 0.0000000                      0.0000000 1.0000000 0.0000000
                      0.0000000 0.0000000 1.0000000                      0.0000000 0.0000000 1.0000000
-ATOM  -3: X=0.00000000 Y=0.50000000 Z=0.23283000+ATOM  -4: X=0.00000000 Y=0.25000000 Z=0.45760000
           MULT= 2          ISPLIT=15           MULT= 2          ISPLIT=15
-      -3: X=0.00000000 Y=0.50000000 Z=0.76717000 +      -4: X=0.00000000 Y=0.75000000 Z=0.54240000 
-Ga         NPT=  781  R0=.000050000 RMT= 2.33        Z:  31.00000+Ca         NPT=  781  R0=.000050000 RMT=   2.00000   Z:  20.00000 
 +LOCAL ROT MATRIX:    1.0000000 0.0000000 0.0000000 
 +                     0.0000000 1.0000000 0.0000000 
 +                     0.0000000 0.0000000 1.0000000 
 +ATOM  -5: X=0.00000000 Y=0.00000000 Z=0.00000000 
 +          MULT= 2          ISPLIT=15 
 +      -5: X=0.00000000 Y=0.50000000 Z=0.00000000 
 +Al         NPT=  781  R0=.000100000 RMT=   2.00000   Z:  13.00000
 LOCAL ROT MATRIX:    1.0000000 0.0000000 0.0000000 LOCAL ROT MATRIX:    1.0000000 0.0000000 0.0000000
                      0.0000000 1.0000000 0.0000000                      0.0000000 1.0000000 0.0000000
Line 38: Line 57:
 </file> </file>
  
-この物質のWIEN2kによるNQRの研究は、以下の論文の5.1節が参考になります+この物質のWIEN2kによるNQRの研究は、以下の論文が参考になります
  
-  * https://d-nb.info/1007751142/34+  * https://doi.org/10.1021/jp0740696
  
  
-収束を厳し目にして、通常のSCF計算を実行します +  - 内部構造最適化を行います<code> 
- +$ init_lapw -b -numk 1000 -rkmax 7.0 
-<code> +$ run_lapw -cc 0.001 -ec 0.0001 -fc 0.5 -i 100 -p -min 
-$ init_lapw -b -numk 20000 -rkmax 7.5+</code> 
 +  - 収束を厳し目にして、SCF計算を再実行します<code> 
 +$ init_lapw -b -numk 20000 -rkmax 7.5 -lvns 5
 $ run_lapw -ec 0.000001 -cc 0.000001 -i 100 -p $ run_lapw -ec 0.000001 -cc 0.000001 -i 100 -p
 </code> </code>
Line 57: Line 78:
 ... ...
  
 +:VKCOUL :  VK-COUL convergence: 0.127E-07
        ATOMNUMBER=  1 Ca         VCOUL-ZERO =  0.22527E+00        ATOMNUMBER=  1 Ca         VCOUL-ZERO =  0.22527E+00
 :EFG001:                        EFG         =     0.57135   *10**21  V / m**2 :EFG001:                        EFG         =     0.57135   *10**21  V / m**2
Line 120: Line 142:
 がEFG計算の結果です。 がEFG計算の結果です。
  
-以下では1番目（今の例ではCa）の原子を使って見方を解説します。+以下では1番目の原子（今の例ではCa）を使って見方を解説します。
  
 ==== 電場勾配 ==== ==== 電場勾配 ====
Line 162: Line 184:
 の順に並んでいます。 の順に並んでいます。
  
-右側が主軸変換（対角化）後の結果で、固有値の絶対値が大きい順に$|V_{zz}| \geq |V_{yy}| \geq |V_{xx}|$と定義されますので、$V_{zz} = 0.57135, $V_{yy} = -0.35308$, , $V_{xx} = -0.21827$です。+右側が主軸変換（対角化）後の結果で、固有値の絶対値が大きい順に$|V_{zz}| \geq |V_{yy}| \geq |V_{xx}|$と定義されますので、$V_{zz} = 0.57135$, $V_{yy} = -0.35308$, , $V_{xx} = -0.21827$です。
  
 ==== 主軸回転 ==== ==== 主軸回転 ====
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