| 查看: 883 | 回复: 0 | |||
[求助]
新手计算TiO2求助
|
|
最近在用VASP对体相金红石结构的TiO2做动能截断测试的时候发现：文献里面大多采用400eV,并且晶胞参数和实验值符合的非常好（差别为0.01埃以下）。但是我测试的时候发现400eV并不能保证能量收敛，真正的收敛在800-850eV左右。 下面是一组测试结果。体系为2*2*2大小的超胞，k点密度4*4*6，paw_pbe赝势，ISMEAR=0。其他赝势和k点密度，ISMEAR都测试过，差不多相似的结果。 300 -433.557008 350 -430.159627 400 -429.264379 450 -428.753864 500 -428.606147 550 -428.616493 600 -428.679304 650 -428.739563 700 -428.7912 750 -428.814651 800 -428.829391 850 -428.831697  对此，我有几个疑问: 1 能量趋势表现为先上升后下降趋向于收敛，这是什么原因？截断能比较小的时候，为什么能量很低？ 2 paw_PBE这种赝势是不是k点尽量取偶数个？我看到很多资料都讲k点尽量取以gamma点为中心对称分布，这样各方向一般为奇数。但是文献中使用这种赝势的大多是偶数值。 3 我检查过400eV的输出，晶胞参数确实计算的非常好，但是看我的结果，能量并不收敛啊。而若取800eV截断，晶胞参数和能量都会收敛，但是晶胞参数比实验值变化了0.06埃，虽然不大，但是没有文献报道的那么好，而且计算时间增加太多。 4 即便是将上述体系总能量除以原子数48，能量变化的趋势仍然没有变。这样在能量从上升到下降转变的地方肯定能找到一个“能量/原子”变化小于0.001eV的点，如果取这个点的话事实上是不正确的，因为后面的变化又变大了。 |
回复此楼» 猜你喜欢
2027广东省杰青
已经有5人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有72人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
-
招聘 金属平磨液，抛光液研发工程师
已经有3人回复
-
时间戳又变了8-15
已经有23人回复
-
filecode=后面第一个是大写字母
已经有7人回复
找到一些相关的精华帖子，希望有用哦~
VASP计算TiO2表面能问题 求助
已经有8人回复
- XPS和EDS表征TiO2
已经有20人回复
- 用castep计算TiO2表面能和吸附能
已经有4人回复
- 作为一个初学者，lz表示完全hold不住cp2k。。。求助方法、资料~~
已经有18人回复
- XRD计算的问题
已经有10人回复
- 金红石TiO2(110)面密度是多少？
已经有4人回复
- 【求助】金红石结构TiO2（110）表面体系能带结构计算得到的结果只有一半
已经有5人回复
- 【求助】怎样通过TiO2的吸收光谱 确定能带宽度
已经有17人回复
- 【求助】ＸＲＤ表征数据 计算TiO2锐钛矿相与红晶石相比例
已经有16人回复
- 【求助】如何计算禁带宽度
已经有3人回复
- 【求助】请问从这张紫外图中怎么算 band gap ？
已经有17人回复
- 【请教】Cu2O 和TiO2 那些晶面能够晶格匹配
已经有7人回复
- 【求助】怎样用紫外吸收光谱图来求TiO2的带隙能，急！
已经有13人回复
- 【求助】怎么根据紫外漫反射光谱图来计算TiO2的的带隙能
已经有9人回复
- 【求助】关于 TiO2的禁带宽度的计算
已经有9人回复
- 【求助】TiO2是离子晶体为什么不溶于水？
已经有4人回复
- 【求助】计算出的费米能级为什么和价带顶重合
已经有16人回复
- 【求助成功】请教：castep计算结果中2*Integrated Spin Density的意义？
已经有10人回复
- 【求助】关于钛酸四丁酯制备TiO2
已经有37人回复
- 【求助】用vasp4.6计算GaAs的能带结构
已经有18人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
