加电场VASP和Dmol3计算的结果不一致

菩城

图片1是看到的一篇文献，体系是MoS2上CO2的吸附，首先是C原子靠近S原子物理吸附，作者用Dmol3把电场加到0.0032 a.u.时，发现C原子和S原子相互靠近，变成化学吸附，我打算用VASP重复他的结果。根据sob老师的帖子http://bbs.keinsci.com/thread-189-1-1.html
1 a.u. =51.423 V/Angstrom
那么我在VASP中需要把电场加到0.0032*51.432=0.16458 V 可以得到跟他相似的结果，但结果并不是，我从0.1 V开始加，都加到了0.8 V还是没有发现C原子和S原子成键，很奇怪，附上我0.8 V的结构（图片2）以及INCAR关于电场的设置，麻烦前辈们帮忙分析一下，万分感谢！


#EEF
   EFIELD = 0.8
   LDIPOL=.TRUE.
   IDIPOL= 3

卡开发发

不好说，做几点猜测看看是否有可能解决：
1、检查下每步SCF是否都收敛了，如果没有那有可能是这个原因导致力和能量不一定正确，那么跑出来的结构可能不正常。加电场和偶极修正一般都可能让SCF收敛变困难，所以确实有这种可能性。
2、DMol3外加电场的时候我印象中是没办法进行偶极修正。偶极修正之后就会把加上去的外场抵消掉（我自己私下测试过），也许你可以试试看相同真空层、同样不做偶极修正做一下看看。
3、也有可能是DMol3的赝势和VASP的PAW对于这种情况有差异，但这个就难说了。

菩城

卡开发发 发表于 2019-3-22 13:32
不好说，做几点猜测看看是否有可能解决：
1、检查下每步SCF是否都收敛了，如果没有那有可能是这个原因导致 ...

感谢老师回复！您给出的建议我已参考，具体结果如下：1.计算过程中每步的SCF都收敛。不过我有个疑问，对于复杂体系，有时候前几个离子步里面的电子步并没有完全收敛，达到60步就接着下一个离子了，随着离子步进行最后电子步也会收敛很好，这样的结果有问题吗？
2.已测试在VASP中不加偶极修正，结果仍然跟文献不一致。
3.不清楚Dmol3与VASP计算电场的区别，不过如果两个软件在大概相同的电场下，结果却是不一样的，这感觉有些问题呀。
期待您的回复，再次感谢！

卡开发发

菩城 发表于 2019-3-22 17:04
感谢老师回复！您给出的建议我已参考，具体结果如下：1.计算过程中每步的SCF都收敛。不过我有个疑问，对 ...

1、这种情况最好做个单点看看和优化的时候最后一步离子步看看能量差别多大，如果差的很大的话那就是有问题了，这种情况VASP确实有。
2、3、这就很难说了，涉及到很细节的东西。

菩城

卡开发发 发表于 2019-3-22 22:01
1、这种情况最好做个单点看看和优化的时候最后一步离子步看看能量差别多大，如果差的很大的话那就是有问 ...

好的，感谢老师！

吕小强强

楼主，你好，你们计算吸附能的时候，吸附小分子是否也需要单独加上电场优化一次得到能量？跟不加电场区别大吗？

卡开发发

吕小强强 发表于 2019-5-23 15:35
楼主，你好，你们计算吸附能的时候，吸附小分子是否也需要单独加上电场优化一次得到能量？跟不加电场区别大 ...

应当认为表面和吸附后的结构是在电场作用下的，而分子从无穷远处不受电场作用的气相中过来发生吸附。

一无是处lzc

楼主，你好！请问这篇文献方便给DOI之类的信息吗？我最近也想做类似的计算，想参考一下文章的做法。。。