quantumespresso:gipaw:磁気遮蔽テンソル
Differences
This shows you the differences between two versions of the page.
| Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
| quantumespresso:gipaw:磁気遮蔽テンソル [2020/03/15 20:12] – [出力ファイルの見方] koudai | quantumespresso:gipaw:磁気遮蔽テンソル [2020/03/15 21:54] (current) – 削除 koudai | ||
|---|---|---|---|
| Line 1: | Line 1: | ||
| - | ====== 概要 ====== | ||
| - | 石英 (quartz, SiO2) の磁気感受率と化学シフトの計算を行います。 | ||
| - | |||
| - | |||
| - | ====== 準備 ====== | ||
| - | |||
| - | ===== 計算 ===== | ||
| - | |||
| - | 擬ポテンシャルはGIPAW計算に対応したものを使用します。 | ||
| - | 擬ポテンシャルファイルを開いて、GIPAWの記載があれば確実でしょう。 | ||
| - | |||
| - | 今回、擬ポテンシャルは次のサイトからダウンロードしました | ||
| - | * https:// | ||
| - | |||
| - | <file - quartz.scf.in> | ||
| - | & | ||
| - | calculation = ' | ||
| - | prefix | ||
| - | pseudo_dir | ||
| - | outdir | ||
| - | / | ||
| - | &system | ||
| - | ibrav = 0, | ||
| - | celldm(1) = 4.6415377, | ||
| - | nat = 9, | ||
| - | ntyp = 2, | ||
| - | ecutwfc | ||
| - | nosym = .true. | ||
| - | / | ||
| - | & | ||
| - | conv_thr = 1.0d-10 | ||
| - | / | ||
| - | |||
| - | ATOMIC_SPECIES | ||
| - | Si 28.086 | ||
| - | | ||
| - | |||
| - | ATOMIC_POSITIONS crystal | ||
| - | Si 0.4701 | ||
| - | Si 0.0000 | ||
| - | Si | ||
| - | O | ||
| - | O | ||
| - | O | ||
| - | O 0.2674 | ||
| - | O 0.1465 | ||
| - | O | ||
| - | |||
| - | CELL_PARAMETERS alat | ||
| - | 1.0000000 | ||
| - | 1.0000000 | ||
| - | 0.0000000 | ||
| - | |||
| - | K_POINTS automatic | ||
| - | 4 4 4 1 1 1 | ||
| - | </ | ||
| - | |||
| - | * パラメタの意味は次のとおりです | ||
| - | |||
| - | ^変数^初期値^説明^ | ||
| - | |nosym|.false.|falseにすると、結晶の対称性に応じて等価なk点での計算を省略して計算量を減らす。現在のところGIPAWは一部の結晶構造でこのような計算に対応していない| | ||
| - | |||
| - | $ pw.x < quartz.scf.in > quartz.scf.out | ||
| - | |||
| - | |||
| - | ===== NMR計算 ===== | ||
| - | |||
| - | <file - quartz.nmr.in> | ||
| - | & | ||
| - | job = ' | ||
| - | prefix = ' | ||
| - | tmp_dir = ' | ||
| - | q_gipaw = 0.01 | ||
| - | use_nmr_macroscopic_shape = .true. | ||
| - | / | ||
| - | </ | ||
| - | |||
| - | ^変数^初期値^説明^ | ||
| - | |job|nmr|計算の種類。nmrとした場合は磁気感受率と化学シフトの計算を行う。| | ||
| - | |q_gipaw|0.01|NMRのパラメタは線形応答理論から計算するが、これは有限の波数qで定式化される。実際の測定は一様磁場 (q=0) で行うので、数値計算が破綻しない程度の小さいqで計算してq-> | ||
| - | |use_nmr_macroscopic_shape|.false.|サンプルのマクロな形状を考慮して化学シフトの計算をするかどうか。NMR計算ではtrueにすることを推奨。なお、形状はnmr_macroscopic_shapeの値により指定する（デフォルトではサンプルが球体であると仮定） http:// | ||
| - | |||
| - | |||
| - | 実行は次のようにする。 | ||
| - | |||
| - | $ gipaw.x < quartz.nmr.in > quartz.nmr.out | ||
| - | |||
| - | NMRの計算には時間がかかるので注意。 | ||
| - | |||
| - | ===== 出力ファイルの見方 ===== | ||
| - | |||
| - | quartz.nmr.outに結果が出力されます。 | ||
| - | |||
| - | ==== 磁気感受率 ==== | ||
| - | |||
| - | <file - quartz.nmr.out> | ||
| - | End of magnetic susceptibility calculation | ||
| - | |||
| - | f-sum rule (1st term): | ||
| - | | ||
| - | -0.0000 | ||
| - | | ||
| - | |||
| - | f-sum rule (2nd term): | ||
| - | -1.0457 | ||
| - | -0.0000 | ||
| - | | ||
| - | |||
| - | f-sum rule (should be | ||
| - | | ||
| - | -0.0000 | ||
| - | | ||
| - | </ | ||
| - | * f-sum rule (should be | ||
| - | |||
| - | <file - quartz.nmr.out> | ||
| - | | ||
| - | | ||
| - | -0.002133 | ||
| - | | ||
| - | |||
| - | | ||
| - | -0.002133 | ||
| - | | ||
| - | |||
| - | | ||
| - | | ||
| - | -0.001493 | ||
| - | | ||
| - | |||
| - | | ||
| - | -0.001493 | ||
| - | | ||
| - | |||
| - | | ||
| - | | ||
| - | -0.0101 | ||
| - | | ||
| - | |||
| - | | ||
| - | | ||
| - | -0.0071 | ||
| - | | ||
| - | </ | ||
| - | * （電子間の相互作用を考慮しない）磁気感受率 $\chi^{0}_{ij}$ が chi_bare pGv と chi_bare vGv に出力されます。実際の値は、これらの値の範囲内に入ります。 | ||
| - | |||
| - | ==== 磁気遮蔽テンソル ==== | ||
| - | |||
| - | <file - quartz.nmr.out> | ||
| - | End of magnetic susceptibility calculation | ||
| - | |||
| - | f-sum rule (1st term): | ||
| - | | ||
| - | -0.0000 | ||
| - | | ||
| - | |||
| - | f-sum rule (2nd term): | ||
| - | -1.0457 | ||
| - | -0.0000 | ||
| - | | ||
| - | |||
| - | f-sum rule (should be | ||
| - | | ||
| - | -0.0000 | ||
| - | | ||
| - | |||
| - | | ||
| - | | ||
| - | -0.002133 | ||
| - | | ||
| - | |||
| - | | ||
| - | -0.002133 | ||
| - | | ||
| - | |||
| - | | ||
| - | | ||
| - | -0.001493 | ||
| - | | ||
| - | |||
| - | | ||
| - | -0.001493 | ||
| - | | ||
| - | |||
| - | | ||
| - | | ||
| - | -0.0101 | ||
| - | | ||
| - | |||
| - | | ||
| - | | ||
| - | -0.0071 | ||
| - | | ||
| - | </ | ||
| - | |||
| - | 主軸変換した結果も出力されます。 | ||
| - | 座標軸は $|V_{xx}| \leq |V_{yy}| \leq |V_{zz}|$ となるように選ばれます。 | ||
| - | <file - quartz.efg.out> | ||
| - | | ||
| - | | ||
| - | | ||
| - | | ||
| - | | ||
| - | |||
| - | | ||
| - | | ||
| - | | ||
| - | | ||
| - | |||
| - | | ||
| - | | ||
| - | | ||
| - | | ||
| - | |||
| - | | ||
| - | | ||
| - | | ||
| - | | ||
| - | |||
| - | | ||
| - | | ||
| - | | ||
| - | | ||
| - | |||
| - | | ||
| - | | ||
| - | | ||
| - | | ||
| - | |||
| - | | ||
| - | | ||
| - | | ||
| - | | ||
| - | |||
| - | | ||
| - | | ||
| - | | ||
| - | | ||
| - | |||
| - | | ||
| - | | ||
| - | | ||
| - | | ||
| - | </ | ||
| - | * axis ... 主軸変換後の座標軸の向き。主軸がz軸になるように選ばれます。 | ||
| - | * $C_q = \frac{e^2 Q q}{h}$ ... 四重極結合定数。$V_{zz} = eq$ である | ||
| - | * $\eta = \frac{V_{xx} - V_{yy}}{V_{zz}}$ ... 異方性パラメタ | ||
| - | Siは核スピン1/ | ||
| - | |||
| - | なお、酸素は $I=5/2$ の核スピンを持つので、エネルギー固有値は3次方程式 $x^3 -7(3+\eta^2) x -20(1-\eta^2) = 0$ の解$x$を使って $E = \frac{e^2 Q q}{20} x$ で表されます（導出はNQRの教科書を見てください）。 | ||
| - | NQR周波数は、この3つのエネルギー固有値の差をプランク定数$h$で割ったもので表されるので、3つ出てきます。 | ||
quantumespresso/gipaw/磁気遮蔽テンソル.1584270736.txt.gz · Last modified: 2021/06/27 21:59 (external edit)