作者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/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
※ 发信站: 批踢踢实业坊(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