Differences

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

--- quantumespresso:スピン軌道相互作用 [2018/02/25 18:18] – 作成 koudai
+++ quantumespresso:スピン軌道相互作用 [2024/03/31 00:01] (current) – koudai
@@ Line 8: / Line 8: @@
 スピン軌道相互作用を含めた電子状態を計算する場合、相対論的効果を含んだ擬ポテンシャルをダウンロードします。
 擬ポテンシャルに相対論的効果が含まれて、スピン軌道相互作用の計算に使えるかどうかは、擬ポテンシャルの説明に「full relativistic」という説明があるかどうかで判断できます。
-scalar relativisticですと、擬ポテンシャルの作成にスピン軌道相互作用は入っていますが、電子状態計算にスピン軌道相互作用を入れることができません。
+scalar relativisticですと、擬ポテンシャルの作成にスピン軌道相互作用は入っていますが、SCF計算にスピン軌道相互作用を入れることができません。
 ここではPt.rel-pz-n-rrkjus.UPFという擬ポテンシャルを使うことにします。
 単体のプラチナの結晶構造は面心立方格子です。 次の入力ファイルを作成します。
+systemのところにlspinorb, noncolin, starting_magnetizationが追加されます。
-<file - pt.scf.in>
+<file - Pt.scf.in>
 &control
   calculation = 'scf',
-  prefix='pt',
+  prefix='Pt',
   pseudo_dir = './',
   outdir='./',
@@ Line 41: / Line 42: @@
 </file>
-磁化も指定しないといけませんが、プラチナは磁性体ではないのでstarting_magnetization=0.0としておきます。
+スピン軌道結合の計算では磁化を必ず指定しないといけませんが、プラチナは磁性体ではないので starting_magnetization(1) = 0.0 としておきます。
-見慣れないのは&systemにあるlspinorbとnoncolinの2つだと思います。 それぞれ次を意味します。
+また lspinorb と noncolin を .true. にしておく必要があります。
+それぞれ次を意味します。
-|lspinorb|.true.にするとスピン軌道相互作用を入れた計算を実行する。full relativisticの擬ポテンシャルを使う場合、これがないとエラーになる|
+^変数^説明^
-|noncolin|デフォルトでは、スピンはz軸方向に向くもののみしか考慮されない。これを.true.にすることで任意の方向にスピンを向けることができるようになる。ただし計算時間が増大する。|
+|noncolin|デフォルトでは、スピンはz軸方向に向くもののみしか考慮されない。これを.true.にすることで様々な方向にスピンを向けることができるようになる（ノンコリニア）。ただし計算時間が増大する。|
+|lspinorb|.true.にするとスピン軌道相互作用を入れた計算を実行する。noncolin=.true. のときのみ指定可能。full relativisticの擬ポテンシャルを使う場合は.true.にしておかないと正しい結果が得られない|
 実行は次のようにします
-  $ pw.x < pt.scf.in > pt.scf.out
+  $ pw.x < Pt.scf.in > Pt.scf.out
 バンドを書いたりするのは今までと同様に計算できます。
@@ Line 58: / Line 61: @@
 鉄Feは強磁性体ですのでstarting_magnetizationを有限の値にする必要があります。 擬ポテンシャルにはFe.rel-pbe-kjpaw.UPFを使います。
-<file - fe.scf.in>
+<file - Fe.scf.in>
 &control
   calculation = 'scf'
   pseudo_dir = './',
-  outdir='./',
+  outdir='./tmp/',
 /
 &system
@@ Line 88: / Line 91: @@
 </file>
-  $ pw.x < fe.scf.in > fe.scf.out
+  $ pw.x < Fe.scf.in > Fe.scf.out