wien2k:電場勾配
Differences
This shows you the differences between two versions of the page.
| Both sides previous revision Previous revision Next revision | Previous revision | ||
|
wien2k:電場勾配 [2020/04/02 19:36] koudai [計算結果] |
wien2k:電場勾配 [2020/07/10 22:01] koudai [主軸回転] |
||
|---|---|---|---|
| Line 7: | Line 7: | ||
| プロジェクト名を alpha-CaAlF5とします。 | プロジェクト名を alpha-CaAlF5とします。 | ||
| - | 以下の論文の結晶構造から作成しました。 | + | 以下の論文の結晶構造からstructファイルを作成しました。 |
| * https:// | * https:// | ||
| Line 61: | Line 61: | ||
| * https:// | * https:// | ||
| + | * スライド：http:// | ||
| * この論文によると、構造最適化して計算したほうが良いということですが、今回はチュートリアルですので構造最適化前の結果を示します。 | * この論文によると、構造最適化して計算したほうが良いということですが、今回はチュートリアルですので構造最適化前の結果を示します。 | ||
| ===== SCF計算 ===== | ===== SCF計算 ===== | ||
| - | - 収束を厳し目にして、SCF計算を再実行します< | + | - SCF計算を実行します< |
| - | $ init_lapw -b -numk 2000 -rkmax 7.0 | + | $ init_lapw -b -numk 100 -rkmax 7.0 |
| - | $ run_lapw -ec 0.000001 -cc 0.000001 -i 100 -p | + | $ run_lapw -ec 0.00001 -cc 0.0001 -i 100 -p |
| </ | </ | ||
| ===== 計算結果 ===== | ===== 計算結果 ===== | ||
| Line 77: | Line 78: | ||
| :VKCOUL : VK-COUL convergence: | :VKCOUL : VK-COUL convergence: | ||
| - | | + | |
| - | : | + | : |
| - | | + | |
| - | | + | |
| - | V22M TOT/ | + | V22M TOT/ |
| - | | + | |
| - | V21M TOT/ | + | V21M TOT/ |
| - | 6.90223 0.06182 0.92145 6.98229 | + | 6.90150 0.06185 0.92134 6.98155 |
| - | 0.06182 -3.25365 -0.20537 | + | 0.06185 -3.25329 -0.20536 |
| - | 0.92145 -0.20537 -3.64858 | + | 0.92134 -0.20536 -3.64821 |
| MAIN DIRECTIONS OF THE EFG | MAIN DIRECTIONS OF THE EFG | ||
| Line 97: | Line 98: | ||
| - | | + | |
| - | : | + | : |
| - | | + | |
| - | | + | |
| - | V22M TOT/ | + | V22M TOT/ |
| - | | + | |
| - | V21M TOT/ | + | V21M TOT/ |
| - | -5.20181 -0.04398 0.15546 -5.20427 | + | -5.20136 -0.04393 0.15543 -5.20382 |
| - | -0.04398 -0.49895 -0.72140 | + | -0.04393 -0.49886 -0.72130 |
| - | 0.15546 -0.72140 5.70076 | + | 0.15543 -0.72130 5.70022 |
| - | MAIN DIRECTIONS OF THE EFG | + | MAIN DIRECTIONS OF THE EFG |
| 0.0072 -0.1149 | 0.0072 -0.1149 | ||
| | | ||
| : | : | ||
| - | : | + | : |
| - | | + | |
| - | : | + | : |
| - | | + | |
| - | | + | |
| V22M TOT/ | V22M TOT/ | ||
| | | ||
| V21M TOT/ | V21M TOT/ | ||
| - | -7.05536 -0.01953 | + | -7.05479 -0.01953 |
| - | | + | |
| - | 0.00000 | + | 0.00000 |
| MAIN DIRECTIONS OF THE EFG | MAIN DIRECTIONS OF THE EFG | ||
| Line 137: | Line 138: | ||
| - | | + | |
| - | : | + | : |
| | | ||
| - | | + | |
| - | V22M TOT/ | + | V22M TOT/ |
| | | ||
| V21M TOT/ | V21M TOT/ | ||
| - | -2.25240 0.15886 | + | -2.25237 0.15884 |
| - | 0.15886 0.24454 | + | 0.15884 0.24450 |
| - | 0.00000 | + | 0.00000 |
| MAIN DIRECTIONS OF THE EFG | MAIN DIRECTIONS OF THE EFG | ||
| Line 154: | Line 155: | ||
| : | : | ||
| - | : | + | : |
| - | | + | |
| - | : | + | : |
| - | | + | |
| - | | + | |
| - | V22M TOT/ | + | V22M TOT/ |
| - | | + | |
| - | V21M TOT/ | + | V21M TOT/ |
| - | -2.85833 -0.03752 0.27084 -2.88523 | + | -2.85854 -0.03740 0.27082 -2.88543 |
| - | -0.03752 2.95110 0.23493 | + | -0.03740 2.95115 0.23490 |
| - | 0.27084 0.23493 -0.09276 | + | 0.27082 0.23490 -0.09260 |
| - | MAIN DIRECTIONS OF THE EFG | + | MAIN DIRECTIONS OF THE EFG |
| - | 0.0104 -0.0762 | + | 0.0103 -0.0762 |
| - | -0.0979 | + | -0.0978 |
| : | : | ||
| - | : | + | : |
| ... | ... | ||
| Line 185: | Line 185: | ||
| ==== 電場勾配 ==== | ==== 電場勾配 ==== | ||
| - | | + | |
| は電場勾配 $eq = V_{zz}~\mathrm{[10^{21} V/m^2]}$ の値です。 | は電場勾配 $eq = V_{zz}~\mathrm{[10^{21} V/m^2]}$ の値です。 | ||
| Line 200: | Line 200: | ||
| * 電気四重極モーメント一覧 https:// | * 電気四重極モーメント一覧 https:// | ||
| - | * 例えば27Alですと $Q=0.15~\mathrm{[\times | + | * 例えば27Alですと $Q=0.15~\mathrm{[10^{-28}m^2]}$ 程度です |
| - | * bはバーンという単位で、$1~\mathrm{[b]} = 10^{-28}~\mathrm{[m^2]}$です | + | * bはバーンという単位で、$1~\mathrm{[b]} = 10^{-28}~\mathrm{[m^2]}$です。散乱断面積に使われる単位ですが、NQRでも同じくらいのオーダーの量を扱うために借用されることが多いです。 |
| + | * 四重極結合定数からNQR周波数を求める方法は核スピンの大きさで変わってくるので、教科書で確認してください。 | ||
| ==== 電気四重極テンソル ==== | ==== 電気四重極テンソル ==== | ||
| - | 0.13321 | + | < |
| - | 0.39243 0.21987 0.00000 | + | -2.85854 -0.03740 0.27082 -2.88543 |
| - | 0.00000 0.00000 -0.35308 | + | -0.03740 2.95115 0.23490 |
| + | 0.27082 0.23490 -0.09260 | ||
| + | </ | ||
| - | の左側の3×3行列は電気四重極テンソルで、 | + | の左側の3×3行列は、元の座標系における電気四重極テンソルで、 |
| \begin{equation} | \begin{equation} | ||
| Line 223: | 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} = 2.96921$, $V_{yy} = -2.88543$, , $V_{xx} = -0.08379$です。 |
| ==== 主軸回転 ==== | ==== 主軸回転 ==== | ||
| - | MAIN DIRECTIONS OF THE EFG | + | MAIN DIRECTIONS OF THE EFG |
| - | | + | |
| - | 0.0000 0.0000 | + | -0.0978 1.0000 |
| 主軸回転後の軸の向きです。 | 主軸回転後の軸の向きです。 | ||
| - | 主軸回転後のz軸の方向（主軸方向）は (0.8957, 1, 0) 、y軸の方向は(0,0,1)、x軸の方向は (1, -0.8957, 0) です。 | + | 列ごとにみてください。 |
| + | 主軸回転後のz軸の方向（主軸方向）は (-0.0029, 1, 0.0765) 、y軸の方向は(1, 0.0103, -0.0978)、x軸の方向は (0.0986, -0.0762, 1) です。 | ||
| ==== 非対称パラメータ ==== | ==== 非対称パラメータ ==== | ||
| - | :ETA001: | + | :ETA005: |
| 非対称パラメータ $\eta = \frac{V_{xx} - V_{yy}}{V_{zz}}$ の値です。 | 非対称パラメータ $\eta = \frac{V_{xx} - V_{yy}}{V_{zz}}$ の値です。 | ||
| Line 248: | Line 250: | ||
| $ grep ": | $ grep ": | ||
| * k点数に対しても収束が遅いので、必ずk点数に関する収束も確認してください。特に金属の場合はかなり多くのk点をとらないと収束しません。 | * k点数に対しても収束が遅いので、必ずk点数に関する収束も確認してください。特に金属の場合はかなり多くのk点をとらないと収束しません。 | ||
| - | * 周辺の原子位置に敏感な量ですので、実験と合わないときは構造最適化を行ってから計算してみるとよいでしょう | + | * 周辺の原子位置に敏感な量ですので、基本的には内部座標の構造最適化と組み合わせて計算します。run_lapw -p -fc 0.1 -min ... |
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
