作者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/cn.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