作者peter308 (pete)
看板PhD
標題Fw: [問題] 版上有人熟Quantum Espresso的嗎?
時間Sat Aug 27 20:09:04 2016
※ [本文轉錄自 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://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
※ 發信站: 批踢踢實業坊(ptt.cc)
※ 轉錄者: peter308 (36.231.170.153), 08/27/2016 20:09:04
1F：推 pttnowash: 我熟Espresso 每天來一杯 08/28 12:11
2F：推 concil: 可惡，想喝。 08/28 15:56
3F：推 JustOneJFS: 沒算過游離能 不過這裡你假設你算的是一整個crystal 08/30 06:57
4F：→ JustOneJFS: 你確定會得到跟single atom一樣的特性嗎？ 08/30 06:57
5F：→ JustOneJFS: 又 你設定這個晶體是正方晶格 unit cell只有一顆原子 08/30 06:58
6F：→ JustOneJFS: 那atomic positions理應在0 0 0 怎麼會在0 0 z? 08/30 06:59
7F：→ JustOneJFS: 最後就算金屬而言這個degauss似乎設得有點低 不過如果 08/30 06:59
8F：→ JustOneJFS: 計算能夠收斂那也許沒什麼關係 08/30 07:00
9F：→ Casimir: 你是只做這個帶+1 價的計算，還是計算+1，0，-1的間的能 08/31 22:03
10F：→ Casimir: 差？ 08/31 22:03
11F：推 seraphwind: What else 09/03 20:48