用6-31G(d,p)和6-311G(d,p)做TDDFT计算，耗时差异巨大

zhangzb

老师好，我在用Gaussian做TDDFT计算的时候用的是B3LYP-D3/6-311G(d,p)算了20个states, 算了24个小时都没有结束。然后我试着把基组换成B3LYP-D3/6-31G(d,p)其他参数不变，1h20min就结束了。我想请问：
1. 在用Gaussian做TDDFT计算的时候，为何6-31G(d,p)和6-311G(d,p)耗时的差异如此之大。因为我之前对比过他们在做opt freq的时候，6-311G(d,p)是比6-31G(d,p)用时要多，但没有多出这么多；
2. 我看到Sob老师有讲过，计算TDDFT的时候，6-311G*对比6-31G*提升其实比较有限。因为我在做opt freq的时候是用的B3LYP-D3/6-311G(d,p)，如果这里用B3LYP-D3/6-31G(d,p)做TDDFT是否合理，到时候审稿人会不会有意见。
谢谢！

sobereva

监控输出文件，弄清楚davidson迭代都做了多少轮，不会判断就把输出文件压缩后上传
没有体系特征没法告诉你是否有必要用6-311G**，体系大小都不明确。

zhangzb

sobereva 发表于 2022-10-15 17:15
监控输出文件，弄清楚davidson迭代都做了多少轮，不会判断就把输出文件压缩后上传
没有体系特征没法告诉你 ...

谢谢老师！中性有机分子，总共100个原子，含有C,H,S元素。使用6-311g(d,p)计算的已经迭代了39轮，还没有一个state converge。而使用6-31g(d,p)计算的迭代了9轮，所有的states就converge 了。

wzkchem5

B3LYP-D3/6-311G(d,p)的计算是卡在了SCF步骤还是TDDFT步骤？如果是卡在了TDDFT，是因为每步TDDFT迭代时间变长了还是总迭代次数变多了？

zhangzb

wzkchem5 发表于 2022-10-15 17:40
B3LYP-D3/6-311G(d,p)的计算是卡在了SCF步骤还是TDDFT步骤？如果是卡在了TDDFT，是因为每步TDDFT迭代时间变 ...

SCF很快就结束了，卡在了TDDFT, 总的迭代次数变多了。

wzkchem5

zhangzb 发表于 2022-10-15 10:32
谢谢老师！中性有机分子，总共100个原子，含有C,H,S元素。使用6-311g(d,p)计算的已经迭代了39轮，还没有 ...

是不是给的内存太少了？如果内存刚好够6-31G(d,p)的计算但是不够6-311G(d,p)的计算，那么程序就会限制Davidson方法的展开空间的维度来节省内存，导致收敛变慢

zhangzb

wzkchem5 发表于 2022-10-15 18:51
是不是给的内存太少了？如果内存刚好够6-31G(d,p)的计算但是不够6-311G(d,p)的计算，那么程序就会限制Dav ...

6-31G(d,p)和6-311G(d,p)用的是同一台服务器，设置相同，都是%nprocshared=16，%mem=128GB，因为这台服务器只有16个核心。现在使用6-311G(d,p)计算的还在迭代，但依然没有看到有state converge。
我用另外一台服务器计算相同的分子，也用6-311G(d,p)，设置成%nprocshared=32，%mem=128GB，现在经过3h20min迭代到29轮的时候，可以看到有2个states converge.
PS:这两台服务器的单核性能差不多。我下次可以把内存加大些。谢谢您！

zhangzb

wzkchem5 发表于 2022-10-15 18:51
是不是给的内存太少了？如果内存刚好够6-31G(d,p)的计算但是不够6-311G(d,p)的计算，那么程序就会限制Dav ...

老师您好！我按照您的建议把内存增加到512GB, %nprocshared=16，%mem=512GB,结果还是一样，16h之后还在迭代，没有一个state收敛。
然后我按照Sob老师的方法（http://sobereva.com/348）用6-311G(d,p), %nprocshared=16，%mem=128GB进行分批计算看看他们耗时是怎么样的。先计算states(1~5), 用时1h23min结束。接着算states(6~10), 1h28min结束。后面的states(11~15)正在进行中。