关于Cu过渡金属的f,g函数的系数问题

TanJing

各位老师好，我的催化剂体系是含Cu金属的配合物，这个时候老师需要我们采用SDD+f,g函数基组进行优化，课题组往届只有关于Ru金属的函数式，如下；我查找了相关资料没有发现关于改内容的描述，请问这个地方的系数应该怎么修改呢
Ru 0 S   3 1.00       0.000000000000
      0.7936570000D+01 -0.1119665531D+01
      0.5984245000D+01  0.1445329211D+01
      0.4882220000D+01  0.6261652615D+00
S   1 1.00       0.000000000000
      0.1144624000D+01  0.1000000000D+01
S   1 1.00       0.000000000000
      0.5230170000D+00  0.1000000000D+01
S   1 1.00       0.000000000000
      0.1175730000D+00  0.1000000000D+01
S   1 1.00       0.000000000000
      0.4805000000D-01  0.1000000000D+01
S   1 1.00       0.000000000000
      0.1600000000D-01  0.1000000000D+01
P   2 1.00       0.000000000000
      0.3754609000D+01 -0.4722656422D+01
      0.2916571000D+01  0.4990908318D+01
P   2 1.00       0.000000000000
      0.1048675000D+01  0.7285466549D+00
      0.5073200000D+00  0.3039042812D+00
P   1 1.00       0.000000000000
      0.2673980000D+00  0.1000000000D+01
P   1 1.00       0.000000000000
      0.6974800000D-01  0.1000000000D+01
P   1 1.00       0.000000000000
      0.2292700000D-01  0.1000000000D+01
D   4 1.00       0.000000000000
      0.6009913000D+01 -0.3271599963D-01
      0.2104280000D+01  0.2657391970D+00
      0.9215000000D+00  0.4812397945D+00
      0.3885980000D+00  0.4099777953D+00
D   1 1.00       0.000000000000
      0.1528360000D+00  0.1000000000D+01
D   1 1.00       0.000000000000
      0.5100000000D-01  0.1000000000D+01
F   1 1.00       0.000000000000
      0.4780000000D+00  0.1000000000D+01
F   1 1.00       0.000000000000
      0.1666000000D+01  0.1000000000D+01
G   1 1.00       0.000000000000
      0.1057000000D+01  0.1000000000D+01
****

zjxitcc

完全不需要修改，直接用def2-TZVP基组即可，该基组对前四周期元素是全电子基组，从第五周期开始自带赝势的定义。如果课题组算力很强，可以考虑在高精度单点时考虑使用def2-TZVPP, def2-QZVP, def2-QZVPP等基组。gjf文件中的书写格式可以参考Sob老师这篇博文
详解Gaussian中混合基组、自定义基组和赝势基组的输入
http://sobereva.com/60

“Ru金属的函数式”不知所云。如果你自己在SDD基础上添加一些来路不明的高角动量基函数，这就像“五金店里有卖成型上好的电锯，而且是免费可获取的，你说你不要，你偏要自己去手工制作一把斧头，要体验古代人工砍树”。

sobereva

SDD额外加极化函数是很不推荐的做法，放着现成的高质量的def2系列不用，效仿某些文献手动加额外极化函数又麻烦、质量又没法保证。虽然def2系列对Cu是全电子基组，用SDD赝势基组原理上能节约时间，但只有区区一个Cu（假定是单核过渡金属配合物），耗时也差不了哪去

另外，如果你的目的是几何优化，加g那么高角动量基函数纯粹白浪费时间，更何况DFT本身对极化函数的角动量要求又不高。SDD标配的赝势基组用于几何优化就足够了，别被以讹传讹

TanJing

zjxitcc 发表于 2026-6-26 17:35
完全不需要修改，直接用def2-TZVP基组即可，该基组对前四周期元素是全电子基组，从第五周期开始自带赝势的 ...

好的 谢谢老师 请问这个对于H元素而言 没有加极化函数这个影响大吗

TanJing

sobereva 发表于 2026-6-26 22:09
SDD额外加极化函数是很不推荐的做法，放着现成的高质量的def2系列不用，效仿某些文献手动加额外极化函数又 ...

好的 谢谢老师

zjxitcc

TanJing 发表于 2026-6-29 12:22
好的 谢谢老师 请问这个对于H元素而言 没有加极化函数这个影响大吗

如果你想对H元素用SDD，并且还想再SDD基础上再添加极化函数，那不要这么做，会让人觉得你很奇怪或不专业。如果你用现成的6-31G(d,p)，括号里的p就是对H的极化函数。

如果你研究的体系涉及氢转移、氢键等因素，应当对H加极化函数。即使不涉及氢转移、氢键等因素，在计算量允许的情况下，一般也会对H使用极化函数（注意动词是使用，不是自己手动加，因为很多现成的基组都对H有极化函数，根据需求选用即可）

TanJing

zjxitcc 发表于 2026-6-29 13:17
如果你想对H元素用SDD，并且还想再SDD基础上再添加极化函数，那不要这么做，会让人觉得你很奇怪或不专业 ...

好的谢谢 我采用的是6-31g*  没有对H加极化函数 这里应该要重算了