高斯计算碳材料吸附问题

刀刀

各位老师们好，新手想求助一下。最近在做掺氮碳材料的吸附计算，用的Gaussian，优化方法是B3LYP-D3/6-311g(d,p) 。在吸附H2O的结构优化计算时发现和文献的结构差异很大，文献用的Vasp计算。图片左边是我优化的结果，右边是文献的；附件是我的.gif文件。想问问可能是什么原因导致的，因为我换了很多位置吸附，结果基本是图片上的那样，想问问有什么办法调整吗？谢谢老师们。

Uus/pMeC6H4-/キ

.gjf文件里面明明有em=gd3的色散校正设置，为什么计算方法说成B3LYP而不是B3LYP-D3呢？而且.gjf里面的体系有氮原子，和截图不是一个体系，太容易造成困惑了。

当前计算和文献的计算在软件、计算级别上差别明显，很可能不具备可比性，参考下述社长博文。
谈谈重复不出来计算化学文献里的数据的可能原因
http://sobereva.com/678（http://bbs.keinsci.com/thread-38911-1-1.html）

而且基底只用蒄这么小的分子且不冻结边缘，按簇模型计算并不能真实表现二维材料的表面。本来表面吸附就合适按周期性体系用CP2K等计算，如果真要截成簇模型按孤立体系用Gaussian计算，也得像这两篇博文一样小心地建模。
18碳环（cyclo[18]carbon）与石墨烯的相互作用：基于簇模型的研究一例
http://sobereva.com/615（http://bbs.keinsci.com/thread-25306-1-1.html）
使用量子化学程序基于簇模型计算金属表面吸附问题
http://sobereva.com/540（http://bbs.keinsci.com/thread-16532-1-1.html）

刀刀

Uus/pMeC6H4-/キ 发表于 2025-1-6 14:28
.gjf文件里面明明有em=gd3的色散校正设置，为什么计算方法说成B3LYP而不是B3LYP-D3呢？而且.gjf里面的体系 ...

好的谢谢老师，我再仔细看看，.gif文件放错了，已更正

sobereva

又不分gif和gjf