在Gaussian16中计算水团簇这种氢键体系时泛函和基组的选择问题

hzz_xizhe

各位老师好，我的组研究的体系是HX(X=F，Cl，Br，I)分子和水分子组成的HX(H2O)n的氢键团簇体系，其中HX会发生解离。对这个体系进行结构优化，计算单点能和加电场优化，泛函和基组的选择一直比较迷茫（计算软件为Molclus+Gaussian16）。师兄师姐们研究的是HX(X=F，Cl，Br)分子和水分子组成的HX(H2O)n团簇，我研究的是HI(H2O)n团簇，一共有下面三个问题：
1. 结构优化：师兄师姐在研究HX(X=F，Cl，Br)分子和水分子组成的HX(H2O)n团簇时，采用B3LYP-D3(BJ)/6-311+G**进行结构优化同时进行振动频率分析，我不理解师兄师姐这边为什么要在基组上加弥散？（我想可能是因为体系是弱相互作用体系）在阅读sob老师的谈谈弥散函数和“月份”基组 - 思想家公社的门口：量子化学·分子模拟·二次元后，我觉得结构优化的时候并不需要加弥散，所以我在结构优化和振动频率分析HI-(H2O)n时，对H和O的基组采用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碳环的超强调控作用 - 思想家公社的门口：量子化学·分子模拟·二次元中提到电子会因为外电场的作用而被严重极化，因此需要加弥散函数。决定基组对H和O采用6-311+G**,对I采用ma-TZVP，方法依旧是B3LYP-D3(BJ)，这样的选择是否合理？

请老师们不吝赐教。

wzkchem5

你需要多高精度？2kcal/mol的精度和0.5kcal/mol的精度，答案是截然不同的

hzz_xizhe

wzkchem5 发表于 2025-10-2 13:04
你需要多高精度？2kcal/mol的精度和0.5kcal/mol的精度，答案是截然不同的

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