| 查看: 837 | 回复: 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的点，如果取这个点的话事实上是不正确的，因为后面的变化又变大了。 |
回复此楼» 猜你喜欢
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
假如你的研究生提出不合理要求
已经有6人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
不自信的我
已经有12人回复
-
所感
已经有4人回复
-
论文终于录用啦！满足毕业条件了
已经有28人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
找到一些相关的精华帖子，希望有用哦~
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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人