[問題] 版上有人熟Quantum Espresso的嗎?

[peter308 (pete)]

 抱歉打擾各位~~

 我目前想重複一篇文章的數值

 可是我試了很久都重複不了, 我把文章所提共的設定用來跑QE

 但不論我怎麼試就試著不出來

 能否請教如果有熟QE或化學背景知識的, 幫我看一下

 說不定我是卡在很簡單的觀念上而不自知

 感謝!!!

 這是我的input script


 &control
    calculation = 'scf'
    restart_mode='from_scratch'
    title='Au3'
    pseudo_dir = "/opt/espresso-5.0.1/pseudo/"
    outdir = "/XXX/Au_Q=+1/"
    prefix='Au3'
    nstep = 1000
 /
&system
    ibrav = 1,
    celldm(1)= 40.0
    nat = 1
    ntyp = 1
    ecutwfc = 30.0
    ecutrho = 150.0
    occupations='smearing', smearing='gauss', degauss=0.0022,
    tot_charge = +1
    assume_isolated = 'makov-payne'
    nspin = 2
    starting_magnetization(1) = 0
    nbnd = 20
 /
&electrons
    conv_thr =  1.D-6
    mixing_beta = 0.7D0
 /
&ions
 /
ATOMIC_SPECIES
 Au  196.96 Au.pbe-van_ak.UPF
ATOMIC_POSITIONS (angstrom)
 Au                 0.00000000    0.00000000    1.73948706 0 0 0
K_POINTS Gamma



 我其實就只是想用QE 來算 Au (single atom) 的第一游離能和親合力

 可是我怎麼算就是算不出和paper一樣的值

 能否請各位幫忙看一下script並給一些建議?

 萬分感激!!!!!!


 能幫忙解決我的問題的 願給1萬p幣作為酬勞



http://tinyurl.com/jg4qvoj

文章pdf可在此下載



--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 36.231.170.153
※ 文章網址: https://www.ptt.cc/bbs/Chemistry/M.1472298575.A.E69.html
※ 編輯: peter308 (36.231.170.153), 08/27/2016 19:57:46
※ 編輯: peter308 (36.231.170.153), 08/27/2016 19:59:17
※ 編輯: peter308 (36.231.170.153), 08/27/2016 20:04:21
→ kkarthur: 考慮一下spin multiplicity的問題
→ kkarthur: starting_magnetization 建議查一下手冊相關資訊
→ peter308: 感謝1F, 我有考慮Multiplicity問題,但數值依舊沒改善
→ peter308: 我後來想到應該是舊版本的psedupotential basis set問題
→ peter308: 新的pseduopotential目前都可以從資料庫下載 但paper
→ peter308: 使用的非常可能是一個非相對論性的pseudopotential,我
→ peter308: 目前在想辦法用舊的QE的工具庫去製造無相對論性效應的
→ peter308: pseudopotential基組 看能否去解決此問題
→ peter308: 相對論性對於金這個材料是非常重要的,有或沒有差很多
→ fauna: 你算的和paper差多少kcal/mol?
主要是electronegativity. 我算出來和paper差0.35 eV
hardness 則是差不多一樣
→ touyalin: conv_thr =  1.D-6 這設定跟 paper 一樣嗎?
conv_thr 似乎影響不大, 我有調到1.D-8 但總能量算出來一樣
※ 編輯: peter308 (36.231.170.153), 08/28/2016 14:35:03
※ 編輯: peter308 (36.231.170.153), 08/28/2016 14:35:45
推 profya: 原文似乎沒有提到有考慮makov-payne修正
→ profya: 不知道是否是差在這邊
→ kancu: QE版本差異跟編譯用的數學函式庫都會有影響