| 查看: 2594 | 回复: 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 ---------------------------------------------------------------------------------------------------------- |
» 猜你喜欢
不自信的我
已经有11人回复
北核录用
已经有3人回复
要不要辞职读博?
已经有6人回复
实验室接单子
已经有3人回复
磺酰氟产物,毕不了业了!
已经有8人回复
求助:我三月中下旬出站,青基依托单位怎么办?
已经有10人回复
26申博(荧光探针方向,有机合成)
已经有4人回复
论文终于录用啦!满足毕业条件了
已经有26人回复
2026年机械制造与材料应用国际会议 (ICMMMA 2026)
已经有4人回复
Cas 72-43-5需要30g,定制合成,能接单的留言
已经有8人回复
» 本主题相关价值贴推荐,对您同样有帮助:
为什么采用CASTEP计算CaTiO3的光学性质其中的静态介电常数与实验值相差很大
已经有14人回复
求助结构解析。。。。
已经有16人回复
关于碳材料的内部结构
已经有5人回复
vasp计算过渡态问题
已经有13人回复
优化和性质计算成功,但是分析结果时键结构图看不了,怎么办?
已经有3人回复
几何优化
已经有7人回复
结构收敛与scf收敛的关系
已经有7人回复
【求助】VASP中,“结构优化”和“电子自洽计算”的关系?
已经有8人回复
【求助】用gaussian view打开优化好的结构后,将断裂键连上,是否影响后面计算
已经有9人回复
【求助】castep结构优化与性质计算
已经有5人回复
【求助】为什么优化结果总是有虚频
已经有5人回复
siestafonda
至尊木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 20139.2
- 散金: 5
- 红花: 6
- 帖子: 1465
- 在线: 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












回复此楼