| 查看: 2445 | 回复: 6 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
zhaoyue_ucla铁虫 (初入文坛)
|
[求助]
钨（tungsten）cohesive energy计算
|
|
|
本人QE刚入门小白，W的scf计算input 如下，请各位大侠帮忙看看问题（我算的结果5.8ev跟实际结果8.5ev差距较大。如果有可能的话，有QE高手愿意语音或者视频帮助解决问题，可以给30~50RMB/per hour的酬劳，不甚感激。 &CONTROL calculation = 'scf' , pseudo_dir = './' , prefix = 'W_exc' , lkpoint_dir = .false. , / &SYSTEM ibrav = 1, celldm(1) = 5.964, nat = 2, ntyp = 1, ecutwfc = 300 , occupations = 'smearing' , degauss = 0.06 , smearing = 'marzari-vanderbilt' , / &ELECTRONS mixing_beta = 0.7 , / ATOMIC_SPECIES W 183.84000 W.pw-mt_fhi.UPF ATOMIC_POSITIONS alat W 0.000000000 0.000000000 0.000000000 W 0.500000000 0.500000000 0.500000000 K_POINTS automatic 6 6 6 1 1 1 对ecutwfc的测试： cut(Ry) 30 50 150 300 E（Ry） -129.05979680 -148.60944687 -209.35569942 -210.14896228 350 370 -210.15691865 -210.16924485 为啥我都测到370Ry, 总能还不能完全收敛啊？（变化>0.01Ry) 以及K点设置的时候如何确定要不要平移啊？ |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有19人回复
-
国社科又开始会评了，不知道这次命运如何
已经有9人回复
-
能否申诉？
已经有5人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 急求解决-计算Fe的cohesive energy
已经有3人回复
- 【求助】cohesive energy如何用CASTEP计算
已经有17人回复
huai
木虫 (职业作家)
- 应助: 4 (幼儿园)
- 金币: 1137.9
- 散金: 118
- 红花: 6
- 帖子: 4538
- 在线: 554小时
- 虫号: 420388
- 注册: 2007-07-09
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
7楼2017-03-08 11:46:55
liqizuiyang
木虫 (著名写手)
- 应助: 59 (初中生)
- 金币: 6859.5
- 散金: 261
- 红花: 94
- 帖子: 2093
- 在线: 708.9小时
- 虫号: 2366686
- 注册: 2013-03-21
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2014-11-18 14:24:15
liqizuiyang
木虫 (著名写手)
- 应助: 59 (初中生)
- 金币: 6859.5
- 散金: 261
- 红花: 94
- 帖子: 2093
- 在线: 708.9小时
- 虫号: 2366686
- 注册: 2013-03-21
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2014-11-18 14:25:36
zhaoyue_ucla
铁虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 104.9
- 帖子: 14
- 在线: 18.7小时
- 虫号: 3465223
- 注册: 2014-10-10
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
|
非常感谢你的回复！ 还有问题想请教你一下 1.celldm(1)我是去wikipedia上找的晶格常数5.964,请问你说的优化怎么进行呀，我刚接触这个方向不是很懂。 2.问一个小白问题，ibrav=2,nat= 1,ntyp=1,W 0 0 0和 ibrav = 1,nat = 2,ntyp=1，W 0 0 0 , W 0.5 0.5 0.5等价吗？我用xcryden画出来两个结果是一样的，但是总能量E却相差一倍，有点困惑。 3.关于赝势的选择。我导师的文章上用的是siesta算的，他让我用QE来算，导师文章上是这么说的：The ab initio calculations are performed using the SIESTA [22] code, which carries out self-consistent density functional theory calculations with norm-conserving pseudopotentials and a flexible linear combination of numerical atomic orbital basis. The core electrons are replaced by Troullier–Martins pseudopotentials [24] and the valence electrons are described by a double-zeta polar- ized orbital basis set [22]. For the exchange–correlation part of the energy functional, we used the local density approximation (LDA) parameterized by [4]. The conver- gence of the energy is achieved for a plane-wave cut-off of 300 Ry. Throughout the calculations, both local and volumetric relaxations are carried out using the Broyden conjugate gradient algorithm and a force tolerance of 0.01 eV A ̊ 􏰛。其实我完全不知道怎么去选这个赝势，我就是在QE官网上找了个跟文章上最接近的赝势:W.pw-mt_fhi.UPF Pseudopotential type: NORMCONS Method: Martins-Troullier Functional type: Perdew-Wang LDA non relativistic 请问我应该选择官网上给的哪一个赝势才能算的比较准呀？或者我是不是得自己做个赝势出来？不甚感激~ |
4楼2014-11-18 15:01:43
投票:
30