| 查看: 1484 | 回复: 3 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
pgi1706木虫 (小有名气)
|
[交流]
晶格常数和能量是这种形式？ 已有2人参与
|
||
|
晕了，算出来的晶格常数和能量是这种形式？ a= 10.891 1 F= -.22064566E+03 E0= -.22064566E+03 d E =-.220646E+03 a= 10.892 1 F= -.22064519E+03 E0= -.22064519E+03 d E =-.220645E+03 a= 10.893 1 F= -.22064626E+03 E0= -.22064626E+03 d E =-.220646E+03 a= 10.894 1 F= -.22064570E+03 E0= -.22064570E+03 d E =-.220646E+03 a= 10.895 1 F= -.22064509E+03 E0= -.22064509E+03 d E =-.220645E+03 a= 10.896 1 F= -.22064517E+03 E0= -.22064517E+03 d E =-.220645E+03 a= 10.897 1 F= -.22064448E+03 E0= -.22064448E+03 d E =-.220644E+03 a= 10.898 1 F= -.22064530E+03 E0= -.22064530E+03 d E =-.220645E+03 a= 10.899 1 F= -.22064452E+03 E0= -.22064452E+03 d E =-.220645E+03 a= 10.900 1 F= -.22064591E+03 E0= -.22064591E+03 d E =-.220646E+03 所谓能量此起彼伏，何为最优结构？ 这是采用4＊4＊4K点优化的结果，我想得到比较精确的能量。 原来采用单G点优化时，a＝10.894时是最低能量，但是采用4＊4＊4后，能量却是现在这个模样？ 可有人遇到这种情况，能否交流一下  ？谢谢 |
回复此楼» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 计算Fe的超晶胞晶格常数与能量
已经有4人回复
- 关于势函数的求助
已经有23人回复
- 晶格常数优化
已经有7人回复
- 【500金币本版活动】“经验力场分子模拟+第一性原理计算问题讨论互助及资源共享”
已经有82人回复
- 【求助】拟合晶格常数
已经有14人回复
- 【求助】晶格常数
已经有6人回复
acridine
木虫 (正式写手)
- 1ST强帖: 11
- 应助: 13 (小学生)
- 金币: 4901.6
- 红花: 13
- 帖子: 706
- 在线: 389.5小时
- 虫号: 305478
- 注册: 2006-12-10
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
贺仪(金币+5): 多谢指教！ 2011-09-09 17:03:41
小木虫(金币+0.5):给个红包，谢谢回帖
贺仪(金币+5): 多谢指教！ 2011-09-09 17:03:41
|
说两点： 1，k点可能不够，要做收敛测试，也就是不断的增加k点数目进行晶格常数的计算，得到一系列的晶格常数值，当k点足够多时，直到得到的晶格常数趋近稳定，这个时候的k数目才是足够的。同样，也应该对Ecut做类似的测试。 2，也是更重要的，要对自己计算精度有一定的认识，这也是和上面说的收敛测试密切相关的。一般来说，能量收敛到meV的量级，晶格常数收敛到0.01埃以下，就已经是比较不错的精度了。从k点数目来看，你使用的参数并非精度极高的计算吧，如果a是原子单位的话，你所追究的这些0.001的差别可能已经超出了你的计算参数所能达到的精度极限了，能量也是同理。因此你所说的这些”能量此起彼伏“可能都是在误差范围之内的。 |
3楼2011-09-09 16:19:18
laplacetalor
银虫 (小有名气)
- 应助: 6 (幼儿园)
- 金币: 690
- 散金: 542
- 帖子: 248
- 在线: 41.1小时
- 虫号: 593064
- 注册: 2008-09-05
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2011-09-09 08:30:15
pgi1706
木虫 (小有名气)
- 1ST强帖: 1
- 应助: 0 (幼儿园)
- 金币: 2895.2
- 散金: 500
- 红花: 1
- 帖子: 116
- 在线: 140.4小时
- 虫号: 604222
- 注册: 2008-09-16
- 专业: 凝聚态物性 II ：电子结构
|
K点我已经优化了，如下 1 Found 1 irreducible k-points: free energy TOTEN = -219.066683 eV 2 Found 1 irreducible k-points: free energy TOTEN = -220.632965 eV 3 Found 4 irreducible k-points: free energy TOTEN = -220.644897 eV 4 Found 4 irreducible k-points: free energy TOTEN = -220.645697 eV 5 Found 10 irreducible k-points: free energy TOTEN = -220.645609 eV 所以我采用4＊4＊4的MP－Kpoints 来计算的，难道是Ecut，之前1＊1＊1的晶胞优化后的Ecut＝400eV，2＊2＊2的超晶胞，还要优化Ecut？ |
4楼2011-09-10 21:41:01
