k-points分别尝试了1*2*1,2*3*1,2*4*2,都会出现报错。经过多次尝试,发现当把k-points设置为2*2*2时不再报错,这是为何?这样的结果可靠么?(模型是非正交体系,其中a的长度几乎是b的两倍,怎么会k值取一样呢) |
卡开发发 发表于 2023-3-26 22:32 好的,谢谢老师指点 |
卡卡西 发表于 2023-3-24 20:31 如果只是表面薄层的话,可以看成个面,不过未必能和单质Co的有所对应,有可能会重构。如果是团簇就只好直接贴。实际上得看你实验测下来的情况,这个没办法臆断。 |
Deepast 发表于 2023-3-24 23:47 嗯嗯好的,谢谢老师 |
如果是固体表面一个个Co原子分开吸附的话,常见做法是研究不同吸附位点和覆盖度 如果是金属簇的话,通常是先找一个簇的构象吸附在表面(可以是文献中稳定的构型,可以是自己构象搜索的,实在不行就从Co晶体截一部分下来)。放上去后做一个AIMD,从结果中找比较稳定的构象进行后续研究。 |
卡开发发 发表于 2023-3-24 11:47 学到了,谢谢老师。还想请教两个问题,我先对CoO切面扩胞,超胞优化之后想在其表面添加少量Co(<10原子,不会完全覆盖CoO表面),以此构建一个金属氧化物表面负载单质的模型,请问在这个模型中Co可以设定为特定晶面么?还是说只能手动添加这些Co原子,没办法构建晶面。 如果上述的特定晶面无法构建或者无需构建,那么我在原有的CoO表面添加几个Co原子形成的团簇,作为Co/CoO催化剂模型,来计算其催化过程,请问这样合理么。 刚开始接触VASP,问题有点多,感谢老师的解答 |
jiangning198511 发表于 2023-3-24 11:03 先说k-spacing的问题,k-spacing选取和这个没有直接关系,就是保证总能和力可靠即可,吸附能比这个收敛实际会更快。这样的情况下,覆盖度θ→0的情况下吸附能肯定是可以比较的。 如果要比较不同覆盖度的情况,那得看要讨论的问题是什么,如果关注单位面积上的吸附能,那就折合到单位面积,如果关注的是1mol分子的情况则折合到1mol分子,得看具体讨论的问题。 |
k点建议取样密度相同(当然,如果绝对能量收敛了差别也不大),直接用vaspkit生成kpoints就行 吸附能的话,比较相同覆盖度就可以了 |
卡开发发 发表于 2023-3-24 10:29 不同晶面的面积大小一般不同,可能导致吸附一个分子后导致其覆盖度不同,这样对比吸附能可行吗? |
卡卡西 发表于 2023-3-23 16:21 不需要,k值都测到收敛就行了,不同体系k收敛快慢不一定一样。 上面的错误来源是三个轴向的k-spacing差异过大,程序才出这个警告。 |
jiangning198511 发表于 2023-3-24 08:45 好的,谢谢老师 |
卡卡西 发表于 2023-3-23 16:21 这个确实不好直接对比,比如面积不同可能覆盖度都有差异,你可以找找金属的不同表面(如111,100等)吸附分子的计算文章看看他们是如何处理这个问题的 |
jiangning198511 发表于 2023-3-23 09:06 老师,再问一个问题,我想对比Co不同晶面对小分子的吸附能力,切出来的两个晶面的晶胞大小不完全一样,比如002晶面对应abc分别为9.8、8.4、22.9,100晶面分别为12.6、8.5、21.9,K-POINTS全取的2*2*1,这样对比两个晶胞的吸附能合适么?我好像之前在一个教程里看说需要晶胞*K值相等才可以比较,但这个很难实现。麻烦老师给解答一下吧 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2024-11-26 22:17 , Processed in 0.177666 second(s), 26 queries , Gzip On.