卡开发发 发表于 2019-8-29 16:01 嗯 但用vaspkit出来的数据纵坐标会很大（数据都超万了），而实验文献给出的图纵坐标就很小（0~1.2左右），那您有没有听过optics.sh 这个脚本，我在网上也搜到了可以用optics.sh 脚本获得光学性质的一些数据 |
小螃蟹 发表于 2019-8-29 15:44 可能直接用vaspkit会简单一些吧。 |
小螃蟹 发表于 2019-8-27 21:37 前面提到过扩胞之后对称性是有变化的，高对称点和路径其实不该用同一个路径，所以看DOS是好的选择。 |
granvia 发表于 2019-8-27 21:54 好 非常谢谢两位的帮助  |
小螃蟹 发表于 2019-8-27 21:37 正如卡开发发建议的，你可以先比较下DOS，看看是否有很大差别。同时看看用DOS确定的带隙是否与从能带上看的一致 |
granvia 发表于 2019-8-27 12:23 嗯 “相同参数”是指进行结构优化及静态自洽计算时除了K点设置不同，就是K点会根据胞的比例进行改变外其它的参数一致，然后两者在计算能带时的路径是一样的。 |
小螃蟹 发表于 2019-8-27 10:54 扩胞前后的第一布里渊区发生了变化，而且1x3x2的不均匀比例扩胞，使得第一布里渊区的形状也发生了变化。那么，你确定前后两个不同布里渊区里的k点路径是一样的吗？ 你提到相同参数，也指k点采样数目吗？按常规方法，扩胞后的计算k点数目按相应缩放比例减少 |
granvia 发表于 2019-8-26 20:07 嗯 非常感谢你们的解答。但有些奇怪的是我用原胞计算带隙和扩成1X3X2的胞用相同的参数及路径算出来的能隙值差别还挺大的，原胞2.75 eV vs. 扩胞3.4 eV。另外，我是想计算有机分子吸附在TiO2表面后整体的能带，可以直观地看出分子吸附后TiO2的CB值与吸附前的变化，但原子个数一旦增多后，计算能带那一步就总是自己跳出，系统提示是内存太大，不知道你们有没有计算过大体系的能带，出现过这种情况吗？ |
granvia 发表于 2019-8-26 22:45 其实也没这个必要去折腾，直接单胞算就行了，又不是啥掺杂或者缺陷体系。 |
卡开发发 发表于 2019-8-26 22:29 对于他这个问题，完全可以将扩胞后的能带unfold成原胞的能带结构（算法应该不难，我记得有这方面的文章），然后比照一下直接计算出的原胞能带，就能找到hint了 |
granvia 发表于 2019-8-26 22:20 道理上确实是这样，对于gap这一个量来说其实路径对不对无所谓，只要k经过VBM和CBM就行。k网格的尺寸只要k做好收敛性测试就行。 |
卡开发发 发表于 2019-8-26 20:28 扩胞使得第一布里渊区变小，但不会改变k空间大小（其也是周期性重复的）。只要选取足够多有代表性的k点路径就比较安全了 另一可能原因是原胞计算时k点数目不够，因为理论上说k矢是连续的，实际计算都是根据玻恩冯卡曼近似用的有限数目k点 |
小螃蟹 发表于 2019-8-26 18:20 1、楼上的说法是没啥问题。扩胞的话其实本不会影响到带隙的结果。但如果随意挑选一个路径正好全都不经过价带顶或导带底的话就有可能给出来的结果有问题，这只是可能性一种。扩胞之后你的格子本身的对称性其实是有变化的，比如一个cubic扩胞1x2x1之后就不是cubic了。 2、和一些细节的处理还是有关系，比如优化晶格的方式等等，有些体系晶格参数对能带是有些影响的。因此参数做好测试是很重要的。 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2025-8-16 18:26 , Processed in 0.176336 second(s), 25 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites