计算化学公社

标题: 在Gaussian16中计算水团簇这种氢键体系时泛函和基组的选择问题 [打印本页]

作者
Author:
hzz_xizhe    时间: 2025-10-2 11:22
标题: 在Gaussian16中计算水团簇这种氢键体系时泛函和基组的选择问题
本帖最后由 hzz_xizhe 于 2025-10-2 11:22 编辑

各位老师好,我的组研究的体系是HX(X=FClBrI)分子和水分子组成的HX(H2O)n的氢键团簇体系,其中HX会发生解离。对这个体系进行结构优化,计算单点能和加电场优化,泛函和基组的选择一直比较迷茫(计算软件为Molclus+Gaussian16)。师兄师姐们研究的是HX(X=FClBr)分子和水分子组成的HX(H2O)n团簇,我研究的是HI(H2O)n团簇,一共有下面三个问题:
1.
结构优化:师兄师姐在研究HX(X=FClBr)分子和水分子组成的HX(H2O)n团簇时,采用B3LYP-D3(BJ)/6-311+G**进行结构优化同时进行振动频率分析,我不理解师兄师姐这边为什么要在基组上加弥散?(我想可能是因为体系是弱相互作用体系)在阅读sob老师的谈谈弥散函数和“月份”基组 - 思想家公社的门口:量子化学·分子模拟·二次元后,我觉得结构优化的时候并不需要加弥散,所以我在结构优化和振动频率分析HI-(H2O)n时,对HO的基组采用6-311G**,对I采用Def2-TZVP的赝势基组。想问一下各位老师我的理解是否有误以及我这样的选择是否合理?
2.单点能计算:师兄师姐们在进行单点能计算时采用的是MP2/jul-cc-pVTZ,但是该方法算的很慢还很占内存。在阅读简谈量子化学计算中DFT泛函的选择 - 思想家公社的门口:量子化学·分子模拟·二次元后,我选用了revDSD-PBEP86-D3(BJ)双杂化泛函。因为体系涉及到弱相互作用,因此选取了带了弥散函数的基组。决定用revDSD-PBEP86-D3(BJ)/ jul-cc-pVTZ,我这样的选择是否合理?此外,考虑到后面可能会算偶极矩,用revDSD-PBEP86-D3(BJ)/Def2-TZVPD可以吗?还是用更好的Def2-TZVPPD
3.
加电场结构优化: sob老师的一篇文章深入揭示外电场对18碳环的超强调控作用 - 思想家公社的门口:量子化学·分子模拟·二次元中提到电子会因为外电场的作用而被严重极化,因此需要加弥散函数。决定基组对HO采用6-311+G**,I采用ma-TZVP,方法依旧是B3LYP-D3(BJ),这样的选择是否合理?
请老师们不吝赐教。


作者
Author:
wzkchem5    时间: 2025-10-2 13:04
你需要多高精度?2kcal/mol的精度和0.5kcal/mol的精度,答案是截然不同的
作者
Author:
hzz_xizhe    时间: 2025-10-2 15:55
wzkchem5 发表于 2025-10-2 13:04
你需要多高精度?2kcal/mol的精度和0.5kcal/mol的精度,答案是截然不同的

结构优化的话,结果能用就行,不用太高精度
单点能希望精度高一点




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3