雑多な記録

User Tools

Site Tools

Trace: • 擬ポテンシャルの作成
quantumespresso:gipaw:擬ポテンシャルの作成

This is an old revision of the document!

Table of Contents

概要

GIPAW計算用の擬ポテンシャルを作成します。

  • pslibrary.1.0.0 を使用するので、以下のサイトからダウンロードしてください。
  • 各原子の擬ポテンシャル生成用の入力ファイルは、次のファイルに含まれています
    • us_ps_high.job
    • us_ps_low.job
    • paw_ps_high.job
    • paw_ps_low.job
      • 上2つがウルトラソフト型で、下2つがPAW型です。
      • highというのはsemicore状態が多く含まれていて、計算の精度は上がりますが計算時間もかかります
      • lowはsemicore状態が少なく、計算時間は少なくて済みますが、必要な軌道が含まれていない場合もあります
        • 基本的にはlowに収録されている擬ポテンシャルで問題ありませんが、油断していると軌道が足りなくて実は変な結果が出ていたということもあるので、highを使った計算も早い段階で試しておくのが無難です。
  • このファイルの他に us_ps_collection.job, us_ps_collection_alt.job, us_lan_collection.job と、これらのusをpawに変えたものも同梱されていますが、これは古いバージョンのpslibraryのデータなので基本的に無視してください。
  • ノルム保存型も収録されていますが (nc_ps_collection.job)、テストが十分に行われていないそうです。
  • その他、PSLibraryを使う前に https://dalcorso.github.io/pslibrary/PP_list.html の Known PPs problems for pslibrary.1.0.0 を必ず一読してください。

擬ポテンシャルの作成

Ceを例に擬ポテンシャル（PBE, Scalar relativistic）を作成します。

ウルトラソフト型

pslibraryの中にある us_ps_high.job から、Ceに該当する部分をコピーしてきて、Ce_us.inというファイルを作成します。

Ce_us.in
 &input
   title='Ce',
   zed=58.0,
   rel=1,
   config='[Xe] 4f0.5 5d1.5 6s2.0 6p0.0',
   iswitch=3,
   dft='pbe'
 /
&inputp
   lpaw=.false.,
   pseudotype=3,
   file_pseudopw='Ce.pbe-spdfn-rrkjus_mypseudo.UPF',
   author='K.S.',
   lloc=-1,
   rcloc=2.0,
   which_augfun='PSQ',
   rmatch_augfun_nc=.true.,
   nlcc=.true.,
   new_core_ps=.true.,
   rcore=1.2,
   tm=.true.
   lgipaw_reconstruction = .true.
/
8
5S  1  0  2.00  0.00  0.90  1.60  0.0
6S  2  0  2.00  0.00  0.90  1.60  0.0
5P  2  1  6.00  0.00  1.00  1.90  0.0
6P  3  1  0.00  0.30  1.00  1.90  0.0
4D  3  2 10.00  0.00  0.90  1.00  0.0
5D  4  2  1.50  0.00  0.95  1.10  0.0
4F  4  3  0.50  0.00  1.00  1.30  0.0
4F  4  3  0.00  1.00  1.00  1.30  0.0
&test
/
6
5S  1  0  2.00  0.00  0.90  1.60  0.0
6S  2  0  2.00  0.00  0.90  1.60  0.0
5P  2  1  6.00  0.00  1.00  1.90  0.0
4D  3  2 10.00  0.00  0.90  1.00  0.0
5D  4  2  1.50  0.00  0.95  1.10  0.0
4F  4  3  0.50  0.00  1.00  1.30  0.0
  • &input
    • Scalar relativistic の場合 rel = 1 と指定
    • 交換相関汎関数としてPBE型を使うので dft='pbe' と変更
  • &inputp
    • file_pseudopw で擬ポテンシャルファイル名を変更
    • author で擬ポテンシャルの作者名を変更
    • GIPAW用に擬ポテンシャルを再構成するため lgipaw_reconstruction = .true. と use_paw_as_gipaw=.true. を追加
  • &test
    • ファイルの最後に &test 以下を追加
    • &testの軌道一覧は&inputpのものをコピーしてきますが、四列目がゼロあるいは負になっているものは削除してください。また軌道の数もそれに合わせて修正します（今の例では8→6）。

ファイルが用意できたら擬ポテンシャルを生成します 

$ ld1.x < Ce_us.in > Ce_us.out

PAW型

pslibraryの中にある paw_ps_high.job から、Ceに該当する部分をコピーしてきて、Ce_paw.inというファイルを作成します。

Ce_paw.in
 &input
   title='Ce',
   zed=58.0,
   rel=1,
   config='[Xe] 4f0.5 5d1.5 6s2.0 6p0.0',
   iswitch=3,
   dft='$gfun'
 /
 &inputp
   lpaw=.true.,
   pseudotype=3,
   file_pseudopw='Ce.pbe-spdfn-kjpaw_mypseudo.UPF',
   author='K.S.',
   lloc=-1,
   rcloc=2.0,
   which_augfun='PSQ',
   rmatch_augfun_nc=.true.,
   nlcc=.true.,
   new_core_ps=.true.,
   rcore=1.2,
   tm=.true.
   lgipaw_reconstruction=.true.
   use_paw_as_gipaw=.true.
 /
8
5S  1  0  2.00  0.00  0.90  1.60  0.0
6S  2  0  2.00  0.00  0.90  1.60  0.0
5P  2  1  6.00  0.00  1.00  1.90  0.0
6P  3  1  0.00  0.30  1.00  1.90  0.0
4D  3  2 10.00  0.00  0.90  1.00  0.0
5D  4  2  1.50  0.00  0.95  1.10  0.0
4F  4  3  0.50  0.00  1.00  1.30  0.0
4F  4  3  0.00  1.00  1.00  1.30  0.0
EOF
  • &input
    • Scalar relativistic の場合 rel = 1 と指定
    • 交換相関汎関数としてPBE型を使うので dft='pbe' と変更
  • &inputp
    • file_pseudopw で擬ポテンシャルファイル名を変更
    • author で擬ポテンシャルの作者名を変更
    • GIPAW用に擬ポテンシャルを再構成するため lgipaw_reconstruction = .true. と use_paw_as_gipaw=.true. を追加

ファイルが用意できたら擬ポテンシャルを生成します 

$ ld1.x < Ce_paw.in > Ce_paw.out
quantumespresso/gipaw/擬ポテンシャルの作成.1593162842.txt.gz · Last modified: 2021/06/27 21:59 (external edit)