应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【求助】优化的结果，scf计算该结构，force为什么仍然很大？
81/1返回列表
查看: 2552  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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
----------------------------------------------------------------------------------------------------------
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-09-21 10:50:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zxzj05(金币+1):谢谢回帖交流 2010-09-21 16:52:39
是否scf和relax用的某些参数不一样，比如k点，比如截断等。
这就证明你的relax并没有收敛，还应该按照更精确的参数算一算
赞一下(3人) 回复此楼
2楼2010-09-21 16:35:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zxzj05(金币+1):谢谢回帖交流 2010-09-22 10:47:42
一般只有在K点和，ecut设置不一样的时候·会出现这样的情况。仔细的对比一下，两个输入文件。
赞一下(2人) 回复此楼
3楼2010-09-22 10:03:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qphll

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
学习了.
回复此楼
Life, Love, Laugh.
4楼2010-09-24 09:09:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
uuv2010(金币+1): 多谢提示 2011-11-17 11:39:38
一般就算K点和截断能完全一致也会出现很大差别。所以每次我做完vc-relax都要重新算一遍自洽计算。
赞一下(2人) 回复此楼
5楼2011-11-17 10:53:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hanyanli0475

铜虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
1楼: Originally posted by didi5158 at 2010-09-21 10:50:11:
vc-relax计算晶体的晶格常数，用计算所得到的结构再进行自洽计算，而所得结果中force，stress，total energy都差别较大，这个是为什么？请大家指教。

vc-relax计算结果：
---------------------------------- ...

个人建议应该进行收敛测试后（ecut k都要测试），在一定收敛精度范围内，再去做relax和scf，并且保证两个输入文件是一样的，这种情况应该就不会发生。
赞一下(1人) 回复此楼
什么时候能算完！！！！
6楼2011-11-17 17:03:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zqclyyq

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
5楼: Originally posted by gemucai at 2011-11-17 10:53:48
一般就算K点和截断能完全一致也会出现很大差别。所以每次我做完vc-relax都要重新算一遍自洽计算。

vc-relax最后它自己不是会进行一次scf吗，如果我把vc最后得到的原子位置拿出来进行一次scf，这两次scf的区别在哪里，老师你知道吗?
赞一下 回复此楼
7楼2015-12-22 11:10:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

siestafonda

至尊木虫 (著名写手)
★
小木虫: 金币+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
赞一下(1人) 回复此楼
8楼2022-06-19 23:18:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 didi5158 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭