| 查看: 3086 | 回复: 7 | |||
didi5158金虫 (小有名气)
|
[交流]
【求助】优化的结果，scf计算该结构，force为什么仍然很大？ 已有7人参与
|
|
vc-relax计算晶体的晶格常数，用计算所得到的结构再进行自洽计算，而所得结果中force，stress，total energy都差别较大，这个是为什么？请大家指教。 vc-relax计算结果： --------------------------------------------------------------------------------------------------------- total energy = -588.70418379 Ry Harris-Foulkes estimate = -588.70418379 Ry estimated scf accuracy < 3.9E-10 Ry Total force = 0.000201 Total SCF correction = 0.000059 SCF correction compared to forces is too large, reduce conv_thr entering subroutine stress ... total stress (Ry/bohr**3) (kbar) P= 0.06 0.00000076 0.00000000 0.00000000 0.11 0.00 0.00 0.00000000 0.00000076 0.00000000 0.00 0.11 0.00 0.00000000 0.00000000 -0.00000024 0.00 0.00 -0.03 bfgs converged in 14 scf cycles and 12 bfgs steps End of BFGS Geometry Optimization Final enthalpy = -588.7041837859 Ry new unit-cell volume = 657.99457 a.u.^3 ( 97.50473 Ang^3 ) CELL_PARAMETERS (alat) 0.781116722 0.000000000 0.000000000 -0.390558361 0.676466924 0.000000000 0.000000000 0.000000000 1.486696131 ATOMIC_POSITIONS (crystal) Li 0.333333333 0.666666667 0.250000000 Li -0.333333333 -0.666666667 -0.250000000 F 0.420499391 -0.010919477 0.122948838 F 0.010919477 0.431418866 0.122948838 F -0.431418866 -0.420499390 0.122948838 F 0.010919477 -0.420499390 0.377051161 F -0.431418866 -0.010919477 0.377051161 F 0.420499390 0.431418866 0.377051161 F -0.420499390 0.010919478 -0.122948838 F -0.010919477 -0.431418866 -0.122948838 F 0.431418866 0.420499390 -0.122948838 F -0.010919477 0.420499390 0.622948838 F 0.431418866 0.010919477 0.622948838 F -0.420499390 -0.431418866 0.622948838 Al 0.666666667 0.333333333 0.250000000 Al -0.666666667 -0.333333333 -0.250000000 Ca 0.000000000 0.000000000 0.000000000 Ca 0.000000000 0.000000000 0.500000000 ---------------------------------------------------------------------------------------------------------- 用上面所得到结构，scf计算结果： ---------------------------------------------------------------------------------------------------------- total energy = -571.49242692 Ry Harris-Foulkes estimate = -571.49242692 Ry estimated scf accuracy < 2.4E-11 Ry Total force = 0.262955 Total SCF correction = 0.000025 entering subroutine stress ... total stress (Ry/bohr**3) (kbar) P=-4722.80 -0.03210576 0.00000000 0.00000000 -4722.92 0.00 0.00 0.00000000 -0.03210576 0.00000000 0.00 -4722.92 0.00 0.00000000 0.00000000 -0.03210340 0.00 0.00 -4722.57 ---------------------------------------------------------------------------------------------------------- |
回复此楼» 猜你喜欢
估计是周四
已经有3人回复
-
基金啊基金
已经有3人回复
-
科研孤儿太难了
已经有15人回复
-
让我中一个面上吧！
已经有10人回复
-
言之凿凿，明天官网会通知放榜？
已经有17人回复
-
时间戳又变了
已经有17人回复
-
今天系统多次维护，明天很可能放榜！
已经有11人回复
-
时间戳变了，能看出什么问题？
已经有23人回复
-
只有每年这种时候来逛逛小木虫
已经有22人回复
-
放榜前的不淡定
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 为什么采用CASTEP计算CaTiO3的光学性质其中的静态介电常数与实验值相差很大
已经有14人回复
- 求助结构解析。。。。
已经有16人回复
- 关于碳材料的内部结构
已经有5人回复
- vasp计算过渡态问题
已经有13人回复
- 优化和性质计算成功，但是分析结果时键结构图看不了，怎么办？
已经有3人回复
- 几何优化
已经有7人回复
- 结构收敛与scf收敛的关系
已经有7人回复
- 【求助】VASP中，“结构优化”和“电子自洽计算”的关系？
已经有8人回复
- 【求助】用gaussian view打开优化好的结构后,将断裂键连上,是否影响后面计算
已经有9人回复
- 【求助】castep结构优化与性质计算
已经有5人回复
- 【求助】为什么优化结果总是有虚频
已经有5人回复
goldenfisher
金虫 (著名写手)
- 1ST强帖: 3
- 应助: 55 (初中生)
- 贵宾: 0.153
- 金币: 5779.1
- 散金: 2908
- 红花: 21
- 帖子: 1552
- 在线: 351.5小时
- 虫号: 117764
- 注册: 2005-11-26
- 性别: GG
- 专业: 原子和分子物理
2楼2010-09-21 16:35:58
dxcharlary
专家顾问 (职业作家)
-
专家经验: +83 - 1ST强帖: 2
- 应助: 251 (大学生)
- 金币: 10507.9
- 散金: 155
- 红花: 65
- 帖子: 3141
- 在线: 758.7小时
- 虫号: 936718
- 注册: 2010-01-04
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 第一性原理
3楼2010-09-22 10:03:40
qphll
金虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 2708.7
- 散金: 2294
- 红花: 8
- 帖子: 554
- 在线: 227.1小时
- 虫号: 89654
- 注册: 2005-08-29
- 性别: GG
- 专业: 碳素材料与超硬材料
4楼2010-09-24 09:09:20
gemucai
木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 1774.6
- 散金: 157
- 红花: 4
- 帖子: 693
- 在线: 201.7小时
- 虫号: 816478
- 注册: 2009-07-27
- 专业: 凝聚态物性 II ：电子结构
5楼2011-11-17 10:53:48
hanyanli0475
铜虫 (小有名气)
- 应助: 11 (小学生)
- 金币: 173.1
- 散金: 32
- 红花: 3
- 帖子: 296
- 在线: 130.5小时
- 虫号: 1160645
- 注册: 2010-12-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
6楼2011-11-17 17:03:55
7楼2015-12-22 11:10:15
siestafonda
至尊木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 20250.2
- 散金: 5
- 红花: 6
- 帖子: 1486
- 在线: 131.5小时
- 虫号: 958984
- 注册: 2010-03-01
- 专业: 凝聚态物性 II ：电子结构
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
| One possibility is that a fixed set of G-vector is selected using the initial cell based on the specified `ecutwfc`. This G-vector set is kept fixed during a vc-relax run (which can modify effectively the stored G-vectors through changes in the unit cell). As such, at the end of vc-relax, the planewave cutoff can differ from direct scf calculation if the unit cell undergoes sizable strain. One solution to this problem is to use the approach by M. Bernasconi et al, [J. Phys. Chem. Solids 56, 501 (1995), doi:10.1016/0022-3697(94)00228-2]. Relevant QE input parameters include: ecfixed, qcutz, and q2sigma see: https://www.quantum-espresso.org/Doc/INPUT_PW.html#idm406 |
8楼2022-06-19 23:18:42
40