雑多な記録

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 18:03] – [概要] koudaiquantumespresso:gipaw:擬ポテンシャルの作成 [2021/06/27 22:04] (current) – external edit 127.0.0.1
Line 16: Line 16:
   * このファイルの他に 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/擬ポテンシャルの作成.1593162192.txt.gz · Last modified: 2021/06/27 21:59 (external edit)