24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

didi5158

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 482.4
帖子: 250
在线: 89.9小时
虫号: 273449
注册: 2006-08-20
性别: GG
专业: 无机材料化学

[交流] 【求助】优化的结果，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
----------------------------------------------------------------------------------------------------------

回复此楼

» 猜你喜欢

职称评审没过，求安慰已经有41人回复
回收溶剂求助已经有7人回复
硝基苯如何除去已经有3人回复
A期刊撤稿已经有4人回复
垃圾破二本职称评审标准已经有17人回复
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有22人回复
EST投稿状态问题已经有7人回复
毕业后当辅导员了，天天各种学生超烦已经有4人回复
求助文献已经有3人回复
三无产品还有机会吗已经有6人回复

» 本主题相关价值贴推荐，对您同样有帮助:

为什么采用CASTEP计算CaTiO3的光学性质其中的静态介电常数与实验值相差很大已经有14人回复
求助结构解析。。。。已经有16人回复
关于碳材料的内部结构已经有5人回复
vasp计算过渡态问题已经有13人回复
优化和性质计算成功，但是分析结果时键结构图看不了，怎么办？已经有3人回复
几何优化已经有7人回复
结构收敛与scf收敛的关系已经有7人回复
【求助】VASP中，“结构优化”和“电子自洽计算”的关系？已经有8人回复
【求助】用gaussian view打开优化好的结构后,将断裂键连上,是否影响后面计算已经有9人回复
【求助】castep结构优化与性质计算已经有5人回复
【求助】为什么优化结果总是有虚频已经有5人回复

1楼 2010-09-21 10:50:11

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

hanyanli0475

铜虫 (小有名气)

应助: 11 (小学生)
金币: 173.1
散金: 32
红花: 3
帖子: 296
在线: 130.5小时
虫号: 1160645
注册: 2010-12-02
性别: GG
专业: 凝聚态物性 II ：电子结构

★
小木虫(金币+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，并且保证两个输入文件是一样的，这种情况应该就不会发生。