计算化学公社

标题: 用Gaussview计算opt+freq，选用哪个基组最为合适？ [打印本页]

作者
Author: lindaiyu 时间: 2025-3-19 17:40
标题: 用Gaussview计算opt+freq，选用哪个基组最为合适？
本帖最后由 lindaiyu 于 2025-3-19 23:05 编辑

各位老师好，本人最近在计算一种含C、H、O、N、Te元素的小分子，结构如下：
[attach]108704[/attach]
图片涉密已删除

利用def2-TZVP基组，B3LYP方法（高斯语言为：# opt freq B3LYP/def2TZVP pop=npa em=gd3bj）计算出的结构模型如下：
[attach]108707[/attach]
结果得出与吡啶上N相连的C原子不与吡啶共面，应该是基组选择的问题。
同结构类型的含硒小分子利用6-311g基组，B3LYP方法（高斯语言为：# opt freq B3LYP/6-311g(d,p) pop=npa em=gd3bj）计算得出的模型结构就很理想。结构模型如下：
[attach]108705[/attach]
图中可以看出与吡啶上N相连的C原子是与吡啶共面的。
后面我尝试了用6-311g*计算（#p b3lyp/6-311G* pop=npa em=gd3bj opt freq），同样地报错超范围了。
Standard basis: 6-311G(d) (5D, 7F)
Atomic number out of range for 6-311G basis.
Error termination via Lnk1e in /share/apps/g09/l301.exe at Wed Mar 19 17:10:16 2025.

所以由于硒和碲都是第IV主族元素，它们构成的这种小分子的结构也应该近似，想请教各位大佬对于这个含碲的小分子的结构计算是否合理？不合理应该选用哪种基组哪种方法计算最精确，效果最好？如果不是基组方法的问题，后续进行opt+freq计算应该加上何种限制？本人量子化学小白，对重原子的弥散作用还不是很懂，一时半会也了解不深刻。还请各位大佬费事费心指点一二。感谢！

作者
Author: sobereva 时间: 2025-3-19 18:14
那叫关键词（keyword），不叫语言（language）

基组名严谨写完整，不得把6-311g(d,p)说成“6-311g基组”，认真看
计算化学中的一些常见不良写法和用词
http://sobereva.com/298（http://bbs.keinsci.com/thread-1358-1-1.html）

对于这类分子结构优化，def2-TZVP不存在任何问题，精度绰绰有余；倘若结果有什么不合理之处（分清楚是你自己的直觉错误还是计算结果不合理），那也绝对不是基组的锅。既然def2-TZVP你算得动，就完全没理由考虑明显更差的6-311g(d,p)。若6-311g(d,p)结果和def2-TZVP差异不可忽略，正说明6-311g(d,p)此时表现烂。

把这些都仔细看了
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）
谈谈赝势基组的选用
http://sobereva.com/373（http://bbs.keinsci.com/thread-5625-1-1.html）
谈谈量子化学研究中什么时候用B3LYP泛函优化几何结构是适当的
http://sobereva.com/557（http://bbs.keinsci.com/thread-17899-1-1.html）

作者
Author: lindaiyu 时间: 2025-3-19 23:19

sobereva 发表于 2025-3-19 18:14
那叫关键词（keyword），不叫语言（language）

基组名严谨写完整，不得把6-311g(d,p)说成“6-311g基组” ...

好的，非常感谢！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)