This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision Next revision Both sides next revision | ||
wien2k:電場勾配 [2020/04/01 21:09] koudai [計算結果] |
wien2k:電場勾配 [2020/04/02 13:44] koudai [電場勾配] |
||
---|---|---|---|
Line 4: | Line 4: | ||
- | ====== 例:CaGa4 ====== | + | ====== 例:α-CaAlF5 |
+ | プロジェクト名を alpha-CaAlF5とします。 | ||
以下の論文の結晶構造から作成しました。 | 以下の論文の結晶構造から作成しました。 | ||
- | * https://www.sciencedirect.com/science/article/ | + | * https://scripts.iucr.org/cgi-bin/paper? |
- | <file - CaGa4.struct> | + | <file - alpa-CaAlF5.struct> |
blebleble | blebleble | ||
- | CXZ LATTICE, | + | CXZ LATTICE, |
MODE OF CALC=RELA unit=bohr | MODE OF CALC=RELA unit=bohr | ||
- | 11.680586 11.559833 | + | 16.463294 13.887597 |
- | ATOM -1: X=0.00000000 | + | ATOM -1: X=0.22110000 |
- | MULT= 1 | + | MULT= 4 |
- | Ca NPT= 781 R0=.000050000 | + | -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 | ||
LOCAL ROT MATRIX: | LOCAL ROT MATRIX: | ||
| | ||
| | ||
- | ATOM -2: X=0.40503000 | + | ATOM -2: X=0.48860000 |
+ | 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= 1.6800 | ||
+ | LOCAL ROT MATRIX: | ||
+ | | ||
+ | | ||
+ | ATOM -3: X=0.00000000 | ||
MULT= 2 ISPLIT=15 | MULT= 2 ISPLIT=15 | ||
- | -2: X=0.59497000 | + | -3: X=0.00000000 |
- | Ga NPT= 781 R0=.000050000 | + | F |
LOCAL ROT MATRIX: | LOCAL ROT MATRIX: | ||
| | ||
| | ||
- | ATOM -3: X=0.00000000 Y=0.50000000 | + | ATOM -4: X=0.00000000 Y=0.25000000 |
MULT= 2 ISPLIT=15 | MULT= 2 ISPLIT=15 | ||
- | -3: X=0.00000000 Y=0.50000000 | + | -4: X=0.00000000 Y=0.75000000 |
- | Ga NPT= 781 R0=.000050000 RMT= 2.33 | + | Ca NPT= 781 R0=.000050000 RMT= 2.0200 |
+ | LOCAL ROT MATRIX: | ||
+ | | ||
+ | | ||
+ | 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 | ||
LOCAL ROT MATRIX: | LOCAL ROT MATRIX: | ||
| | ||
Line 38: | Line 58: | ||
</ | </ | ||
- | この物質のWIEN2kによるNQRの研究は、以下の論文の5.1節が参考になります | + | この物質のWIEN2kによるNQRの研究は、以下の論文が参考になります |
- | * https://d-nb.info/1007751142/34 | + | * https://doi.org/10.1021/jp0740696 |
+ | ===== SCF計算 ===== | ||
- | 収束を厳し目にして、通常のSCF計算を実行します | + | - 内部構造最適化を行います< |
- | + | $ init_lapw -b -numk 1000 -rkmax 7.0 | |
- | < | + | $ run_lapw -ec 0.0001 -cc 0.001 -fc 0.5 -i 100 -p -min |
- | $ init_lapw -b -numk 20000 -rkmax 7.5 | + | </ |
- | $ run_lapw -ec 0.000001 -cc 0.000001 -i 100 -p | + | - 収束を厳し目にして、SCF計算を再実行します< |
+ | $ run_lapw -ec 0.000001 -cc 0.000001 -i 100 -p -NI | ||
</ | </ | ||
===== 計算結果 ===== | ===== 計算結果 ===== | ||
- | SCF計算が終了したら | + | SCF計算が終了したら |
- | <file - CaGe4.scf0> | + | <file - alpha-CaAlF5.scf0> |
... | ... | ||
- | ATOMNUMBER= | + | :VKCOUL : VK-COUL convergence: |
- | : | + | ATOMNUMBER= |
- | | + | : |
- | | + | |
- | V22M TOT/ | + | |
- | | + | V22M TOT/ |
- | V21M TOT/ | + | |
+ | V21M TOT/ | ||
- | | + | |
- | 0.39243 0.21987 | + | 0.06173 -3.25169 -0.20544 |
- | 0.00000 | + | 0.92073 -0.20544 -3.64590 |
- | MAIN DIRECTIONS OF THE EFG | + | MAIN DIRECTIONS OF THE EFG |
- | -0.8957 1.0000 0.0000 | + | 0.0043 0.3809 1.0000 |
- | 0.0000 0.0000 1.0000 | + | 0.0866 1.0000 -0.3784 |
- | : | + | : |
- | : | + | : |
- | | + | |
- | : | + | : |
- | | + | |
- | V22 | + | |
- | V22M TOT/ | + | V22M TOT/ |
+ | | ||
+ | V21M TOT/ | ||
+ | |||
+ | | ||
+ | | ||
+ | 0.15501 | ||
+ | |||
+ | MAIN DIRECTIONS OF THE EFG | ||
+ | 0.0072 -0.1150 | ||
+ | | ||
+ | : | ||
+ | |||
+ | : | ||
+ | |||
+ | |||
+ | | ||
+ | : | ||
+ | | ||
+ | V22 | ||
+ | V22M TOT/ | ||
| | ||
V21M TOT/ | V21M TOT/ | ||
- | 1.21144 -1.40063 | + | -7.05397 -0.01957 |
- | -1.40063 -2.09428 | + | -0.01957 8.98490 |
- | 0.00000 | + | 0.00000 |
- | MAIN DIRECTIONS OF THE EFG | + | MAIN DIRECTIONS OF THE EFG |
- | -0.3667 | + | 0.0012 |
0.0000 | 0.0000 | ||
- | :ANG002: ANGLE WITH OLD X-AXIS = 20.1 | + | :ANG003: ANGLE WITH OLD X-AXIS = 0.1 |
- | :ETA002: | + | :ETA003: |
- | | + | |
- | :EFG003: EFG | + | :EFG004: EFG |
- | | + | |
- | | + | |
- | V22M TOT/ | + | V22M TOT/ |
| | ||
V21M TOT/ | V21M TOT/ | ||
- | 2.15594 | + | -2.25260 0.15884 |
- | -2.22503 -3.63402 | + | 0.15884 0.24637 |
- | 0.00000 | + | 0.00000 |
- | MAIN DIRECTIONS OF THE EFG | + | MAIN DIRECTIONS OF THE EFG |
- | -0.3399 | + | -0.0633 |
0.0000 | 0.0000 | ||
- | :ANG003: ANGLE WITH OLD X-AXIS = | + | :ANG004: ANGLE WITH OLD X-AXIS = 3.6 |
+ | |||
+ | : | ||
+ | |||
+ | |||
+ | | ||
+ | : | ||
+ | | ||
+ | | ||
+ | V22M TOT/ | ||
+ | | ||
+ | V21M TOT/ | ||
+ | |||
+ | | ||
+ | | ||
+ | 0.27118 | ||
+ | |||
+ | MAIN DIRECTIONS OF THE EFG | ||
+ | 0.0104 -0.0763 | ||
+ | | ||
+ | : | ||
- | :ETA003: | + | :ETA005: |
... | ... | ||
Line 120: | Line 183: | ||
がEFG計算の結果です。 | がEFG計算の結果です。 | ||
- | 以下では1番目の原子(今の例ではCa)を使って見方を解説します。 | + | 以下では5番目の原子(今の例ではAl)を使って見方を解説します。 |
==== 電場勾配 ==== | ==== 電場勾配 ==== | ||
Line 139: | Line 202: | ||
* 電気四重極モーメント一覧 https:// | * 電気四重極モーメント一覧 https:// | ||
- | * 例えば43Caですと $Q=4.3~\mathrm{[\times 10^{-28}m^2]}$ です | + | * 例えば43Alですと $Q=14.9~\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 162: | Line 225: | ||
の順に並んでいます。 | の順に並んでいます。 | ||
- | 右側が主軸変換(対角化)後の結果で、固有値の絶対値が大きい順に$|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$です。 |
==== 主軸回転 ==== | ==== 主軸回転 ==== |