雑多な記録

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状態が少なく、計算時間は少なくて済みますが、必要な軌道が含まれていない場合もあります
        • 物質によって必要な軌道が変わるので、状況に応じて使い分けてください
  • このファイルの他に us_ps_collection.job, us_lan_collection.job, us_ps_collection_alt.job と、これらのusをpawに変えたものも同梱されていますが、これは古いバージョンのpslibraryのデータなので基本的に無視してください。
  • ノルム保存型も収録されていますが (nc_ps_collection.job)、テストが十分に行われていないそうです。

擬ポテンシャルの作成

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/擬ポテンシャルの作成.1593161151.txt.gz · Last modified: 2021/06/27 21:59 (external edit)