Differences

This shows you the differences between two versions of the page.

--- wien2k:電場勾配 [2020/04/02 10:53]
koudai [例：CaGa4]
+++ wien2k:電場勾配 [2020/04/02 13:33]
koudai [例：α-CaAlF5]
@@ Line 4: / Line 4: @@
-====== 例：CaGa4 ======
+====== 例：α-CaAlF5 ======
+プロジェクト名を alpha-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
-CXZ LATTICE,NONEQUIV.ATOMS   3   12 C2/m
+CXZ LATTICE,NONEQUIV.ATOMS   5   15 C2/c
 MODE OF CALC=RELA unit=bohr
-.680586 11.559833 11.584210 90.000000 90.000000118.942000
+.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
 .0000000 1.0000000 0.0000000
 .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
+.0000000 1.0000000 0.0000000
+.0000000 0.0000000 1.0000000
+ATOM  -3: X=0.00000000 Y=0.25000000 Z=0.05780000
           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
 .0000000 1.0000000 0.0000000
 .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
-      -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
+.0000000 1.0000000 0.0000000
+.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
 .0000000 1.0000000 0.0000000
@@ Line 38: / Line 58: @@
 </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
+</code>
+  - 収束を厳し目にして、SCF計算を再実行します<code>
 $ init_lapw -b -numk 20000 -rkmax 7.5 -lvns 5
 $ run_lapw -ec 0.000001 -cc 0.000001 -i 100 -p
 </code>
-  * オプション -lvns 5 は、これがないと init_lapw の lstart 実行時に文句を言われるのでつけました。
-  * 金属ですので、大量のk点数が必要です。
-  * 実際に計算するときは、少ないk点数から初めて徐々に増やしていきましょう。
 ===== 計算結果 =====

雑多な記録

User Tools

Site Tools

Differences

Page Tools