作者peter308 (pete)
看板Chemistry
標題[問題] 版上有人熟Quantum Espresso的嗎?
時間Sat Aug 27 19:49:32 2016
抱歉打擾各位~~
我目前想重複一篇文章的數值
可是我試了很久都重複不了, 我把文章所提共的設定用來跑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://webptt.com/m.aspx?n=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
1F：→ kkarthur: 考慮一下spin multiplicity的問題 08/28 02:35
2F：→ kkarthur: starting_magnetization 建議查一下手冊相關資訊 08/28 02:37
3F：→ peter308: 感謝1F, 我有考慮Multiplicity問題,但數值依舊沒改善 08/28 09:18
4F：→ peter308: 我後來想到應該是舊版本的psedupotential basis set問題 08/28 09:18
5F：→ peter308: 新的pseduopotential目前都可以從資料庫下載 但paper 08/28 09:19
6F：→ peter308: 使用的非常可能是一個非相對論性的pseudopotential,我 08/28 09:20
7F：→ peter308: 目前在想辦法用舊的QE的工具庫去製造無相對論性效應的 08/28 09:21
8F：→ peter308: pseudopotential基組 看能否去解決此問題 08/28 09:22
9F：→ peter308: 相對論性對於金這個材料是非常重要的,有或沒有差很多 08/28 09:23
10F：→ fauna: 你算的和paper差多少kcal/mol? 08/28 12:44
主要是electronegativity. 我算出來和paper差0.35 eV
hardness 則是差不多一樣
11F：→ touyalin: conv_thr = 1.D-6 這設定跟 paper 一樣嗎? 08/28 13:54
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
12F：推 profya: 原文似乎沒有提到有考慮makov-payne修正 08/29 01:02
13F：→ profya: 不知道是否是差在這邊 08/29 01:03
14F：→ kancu: QE版本差異跟編譯用的數學函式庫都會有影響 09/07 13:26