| 查看: 2436 | 回复: 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点设置的时候如何确定要不要平移啊？ |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有18人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
中青基了要发朋友圈吗？
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
国社科申报系统有变化
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 急求解决-计算Fe的cohesive energy
已经有3人回复
- 【求助】cohesive energy如何用CASTEP计算
已经有17人回复
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
zhaoyue_ucla
铁虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 104.9
- 帖子: 14
- 在线: 18.7小时
- 虫号: 3465223
- 注册: 2014-10-10
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
5楼2014-11-18 15:06:30
liqizuiyang
木虫 (著名写手)
- 应助: 59 (初中生)
- 金币: 6859.5
- 散金: 261
- 红花: 94
- 帖子: 2093
- 在线: 708.9小时
- 虫号: 2366686
- 注册: 2013-03-21
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
★ ★
ljw4010: 金币+2, 谢谢指导解答！ 2014-11-20 20:38:43
ljw4010: 金币+2, 谢谢指导解答！ 2014-11-20 20:38:43
|
1. 这个可以在站内找下优化晶格常数的教程. 2.QE有两种指定晶体结构的方法，第一种基于晶胞，用内置的十四种布拉伐格子，W是体心立方，所以是2。另一种是用原胞，自己手输基矢坐标。体心立方的基矢不是(a,0,0)(0,a,0)(0,0,a)，你原先是按照(a,0,0)(0,a,0)(0,0,a)来理解的。 3.http://www.quantum-espresso.org/ ... 0&origin_id去这里找标着pz/pw和scalar relativistic的赝势。 如果计算资源允许，推荐用ELK或exciting。这两个程序不使用赝势，并且使用缀加平面波，算W这种重元素可能更准确。 |
6楼2014-11-18 15:15:44
huai
木虫 (职业作家)
- 应助: 4 (幼儿园)
- 金币: 1130.9
- 散金: 118
- 红花: 6
- 帖子: 4537
- 在线: 554小时
- 虫号: 420388
- 注册: 2007-07-09
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
7楼2017-03-08 11:46:55
投票:
30