雑多な記録

User Tools

Site Tools

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

Differences

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

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision		Previous revision
quantumespresso:gipaw:擬ポテンシャルの作成 [2020/06/26 17:46] – [概要] koudaiquantumespresso:gipaw:擬ポテンシャルの作成 [2021/06/27 22:04] (current) – external edit 127.0.0.1
Line 13: Line 13:
         * highというのはsemicore状態が多く含まれていて、計算の精度は上がりますが計算時間もかかります         * highというのはsemicore状態が多く含まれていて、計算の精度は上がりますが計算時間もかかります
         * lowはsemicore状態が少なく、計算時間は少なくて済みますが、必要な軌道が含まれていない場合もあります         * lowはsemicore状態が少なく、計算時間は少なくて済みますが、必要な軌道が含まれていない場合もあります
-          * 物質によって必要な軌道が変わるので、状況に応じて使い分けてください+          * 基本的にはlowに収録されている擬ポテンシャルで問題ありませんが、油断していると軌道が足りなくて実は変な結果が出ていたということもあるので、highを使った計算も早い段階で試しておくのが無難です。
   * このファイルの他に us_ps_collection.job, us_ps_collection_alt.job, us_lan_collection.job と、これらのusをpawに変えたものも同梱されていますが、これは古いバージョンのpslibraryのデータなので基本的に無視してください。     * このファイルの他に us_ps_collection.job, us_ps_collection_alt.job, us_lan_collection.job と、これらのusをpawに変えたものも同梱されていますが、これは古いバージョンのpslibraryのデータなので基本的に無視してください。  
   * ノルム保存型も収録されていますが (nc_ps_collection.job)、テストが十分に行われていないそうです。    * ノルム保存型も収録されていますが (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）を作成します。+酸素Oを例に擬ポテンシャル（PBE, Scalar relativistic）を作成します。
  
 ===== ウルトラソフト型 ===== ===== ウルトラソフト型 =====
  
-pslibraryの中にある us_ps_high.job から、Ceに該当する部分をコピーしてきて、Ce_us.inというファイルを作成します。+pslibraryの中にある us_ps_high.job から、Oに該当する部分をコピーしてきて、Ce_us.inというファイルを作成します。
  
-<file - Ce_us.in>+<file - O_us.in>
  &input  &input
-   title='Ce', +   title='O', 
-   zed=58.0,+   zed=8.,
    rel=1,    rel=1,
-   config='[Xe] 4f0.5 5d1.5 6s2.0 6p0.0',+   config='[He] 2s2 2p4',
    iswitch=3,    iswitch=3,
    dft='pbe'    dft='pbe'
  /  /
-&inputp+ &inputp
    lpaw=.false.,    lpaw=.false.,
    pseudotype=3,    pseudotype=3,
-   file_pseudopw='Ce.pbe-spdfn-rrkjus_mypseudo.UPF',+   file_pseudopw='O.pbe-nl-rrkjus_psl.1.0.0_gipaw.UPF',
    author='K.S.',    author='K.S.',
    lloc=-1,    lloc=-1,
-   rcloc=2.0,+   rcloc=1.1
    which_augfun='PSQ',    which_augfun='PSQ',
    rmatch_augfun_nc=.true.,    rmatch_augfun_nc=.true.,
    nlcc=.true.,    nlcc=.true.,
    new_core_ps=.true.,    new_core_ps=.true.,
-   rcore=1.2,+   rcore=0.7,
    tm=.true.    tm=.true.
-   lgipaw_reconstruction = .true. +   lgipaw_reconstruction=.true. 
-/ + / 
-8 +4 
-5S  1  0  2.00  0.00  0.90  1.60  0.0 +2S  1  0  2.00  0.00  1.00  1.30  0.0 
-6S  2  0  2.00  0.00  0.90  1.60  0.0 +2S  1  0  0.00  1.00  1.00  1.30  0.0 
-5P  2  1  6.00  0.00  1.00  1.90  0.0 +2P  2  1  4.00  0.00  0.90  1.45  0.0 
-6P  3  1  0.00  0.30  1.00  1.90  0.0 +2P  2  1  0.00  0.05  0.90  1.45  0.0 
-4D  3  2 10.00  0.00  0.90  1.00  0.0 + &test 
-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 +2 
-4F  4  3  0.00  1.00  1.00  1.30  0.0 +2S  1  0  2.00  0.00  1.00  1.30  0.0 
-&test +2P  2  1  4.00  0.00  0.90  1.45  0.0
-/ +
-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+
 </file> </file>
  
Line 76: Line 69:
     * file_pseudopw で擬ポテンシャルファイル名を変更     * file_pseudopw で擬ポテンシャルファイル名を変更
     * author で擬ポテンシャルの作者名を変更     * author で擬ポテンシャルの作者名を変更
-    * GIPAW用に擬ポテンシャルを再構成するため lgipaw_reconstruction = .true. と use_paw_as_gipaw=.true. を追加 +    * GIPAW用に擬ポテンシャルを再構成するため lgipaw_reconstruction=.true. を追加 
-  * &test  +  * &test の項目を追加（中身はなし） 
-    * ファイルの最後に &test 以下を追加 +  * 軌道一覧を&testの下に追加 
-    * &testの軌道一覧は&inputpのものをコピーしてきますが、四列目がゼロあるいは負になっているものは削除してください。また軌道の数もそれに合わせて修正します（今の例では8->6）。+    * &inputpのものをコピーしてきますが、四列目がゼロあるいは負になっているものは削除してください。また軌道の数もそれに合わせて修正します（今の例では4->2）。 
 +    * いろいろ試した結果、こうしないといけないことがわかったのですが、どうしてかは謎です。 
  
 ファイルが用意できたら擬ポテンシャルを生成します ファイルが用意できたら擬ポテンシャルを生成します
  
 <code> <code>
-$ ld1.x < Ce_us.in > Ce_us.out+$ ld1.x < O_us.in > O_us.out
 </code> </code>
 +
 +注意：酸素は今の所 QE6.4 以降でエラーが出るので、6.3で作成してください。
  
  
 ===== PAW型 ===== ===== PAW型 =====
  
-pslibraryの中にある paw_ps_high.job から、Ceに該当する部分をコピーしてきて、Ce_paw.inというファイルを作成します。+pslibraryの中にある paw_ps_high.job から、Oに該当する部分をコピーしてきて、O_paw.inというファイルを作成します。
  
-<file - Ce_paw.in>+<file - O_paw.in>
  &input  &input
-   title='Ce', +   title='O', 
-   zed=58.0,+   zed=8.,
    rel=1,    rel=1,
-   config='[Xe] 4f0.5 5d1.5 6s2.0 6p0.0',+   config='[He] 2s2 2p4',
    iswitch=3,    iswitch=3,
-   dft='$gfun'+   dft='pbe'
  /  /
  &inputp  &inputp
    lpaw=.true.,    lpaw=.true.,
    pseudotype=3,    pseudotype=3,
-   file_pseudopw='Ce.pbe-spdfn-kjpaw_mypseudo.UPF',+   file_pseudopw='O.pbe-nl-kjpaw_psl.1.0.0_gipaw.UPF',
    author='K.S.',    author='K.S.',
    lloc=-1,    lloc=-1,
-   rcloc=2.0,+   rcloc=1.1
    which_augfun='PSQ',    which_augfun='PSQ',
    rmatch_augfun_nc=.true.,    rmatch_augfun_nc=.true.,
    nlcc=.true.,    nlcc=.true.,
    new_core_ps=.true.,    new_core_ps=.true.,
-   rcore=1.2,+   rcore=0.7,
    tm=.true.    tm=.true.
    lgipaw_reconstruction=.true.    lgipaw_reconstruction=.true.
    use_paw_as_gipaw=.true.    use_paw_as_gipaw=.true.
  /  /
-8 +4 
-5S  1  0  2.00  0.00  0.90  1.60  0.0 +2S  1  0  2.00  0.00  1.00  1.30  0.0 
-6S  2  0  2.00  0.00  0.90  1.60  0.0 +2S  1  0  0.00  1.00  1.00  1.30  0.0 
-5P  2  1  6.00  0.00  1.00  1.90  0.0 +2P  2  1  4.00  0.00  0.90  1.45  0.0 
-6P  3  1  0.00  0.30  1.00  1.90  0.0 +2P  2  1  0.00  0.05  0.90  1.45  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+
 </file> </file>
  
Line 140: Line 131:
  
 <code> <code>
-$ ld1.x < Ce_paw.in > Ce_paw.out+$ ld1.x < O_paw.in > O_paw.out
 </code> </code>
quantumespresso/gipaw/擬ポテンシャルの作成.1593161207.txt.gz · Last modified: 2021/06/27 21:59 (external edit)