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

看板PhD作者 (pete)時間9年前 (2016/08/27 20:09), 9年前編輯推噓4(407)
留言11則, 5人參與, 最新討論串1/1
※ [本文轉錄自 Chemistry 看板 #1NmNvFvf ] 作者: 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://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 ※ 發信站: 批踢踢實業坊(ptt.cc) ※ 轉錄者: peter308 (36.231.170.153), 08/27/2016 20:09:04

08/28 12:11, , 1F
我熟Espresso 每天來一杯
08/28 12:11, 1F

08/28 15:56, , 2F
可惡,想喝。
08/28 15:56, 2F

08/30 06:57, , 3F
沒算過游離能 不過這裡你假設你算的是一整個crystal
08/30 06:57, 3F

08/30 06:57, , 4F
你確定會得到跟single atom一樣的特性嗎?
08/30 06:57, 4F

08/30 06:58, , 5F
又 你設定這個晶體是正方晶格 unit cell只有一顆原子
08/30 06:58, 5F

08/30 06:59, , 6F
那atomic positions理應在0 0 0 怎麼會在0 0 z?
08/30 06:59, 6F

08/30 06:59, , 7F
最後就算金屬而言這個degauss似乎設得有點低 不過如果
08/30 06:59, 7F

08/30 07:00, , 8F
計算能夠收斂那也許沒什麼關係
08/30 07:00, 8F

08/31 22:03, , 9F
你是只做這個帶+1 價的計算,還是計算+1,0,-1的間的能
08/31 22:03, 9F

08/31 22:03, , 10F
差?
08/31 22:03, 10F

09/03 20:48, , 11F
What else
09/03 20:48, 11F
文章代碼(AID): #1NmOBY9s (PhD)