计算化学公社

标题: 高斯09D01版本进行配合物激发态频率计算，为什么一直循环计算，无法产生频率数据？ [打印本页]

作者
Author: zhficcas 时间: 2020-8-12 10:56
标题: 高斯09D01版本进行配合物激发态频率计算，为什么一直循环计算，无法产生频率数据？
本帖最后由 zhficcas 于 2020-8-12 10:57 编辑

各位好！我用高斯09.D.01版对一个配合物激发态进行频率计算（已优化好结构），输入文件如下。为什么计算一直在L502—L914—L1002几个模块间循环运算，无法产生频率数据？谢谢！09版的无法进行TDDFT 的激发态频率计算吗？谢谢！

输入文件如下：
%nproc=50
%mem=60gb
%chk=s1-Fix-freq.chk
#P freq td pbepbe/def2svp iop(3/76=1000003333,3/77=0666706667)

11

0 1
Cu             0 -0.01514900 0.90378607 -0.04731300
P             -1 1.92337114 -0.48072903 0.03724100
P             -1 -1.98407914 -0.45097403 -0.10357301
O             0 -0.23875902 -1.88311513 1.54247711
N             -1 0.27922102 2.35255417 -1.51248211
N             -1 -0.20594401 3.58701826 -1.25215009
N             -1 -0.08778301 2.22980716 1.53355411
N             -1 -0.53242904 3.48021925 1.27821509
C             0 1.01323107 2.44671617 -2.61900119
C             0 1.00665207 3.75612227 -3.09933622
C             0 0.21569202 4.44164432 -2.18690316
C             0 0.27522702 2.19187216 2.81169320
C             0 0.07845601 3.43892225 3.40809025
C             0 -0.43883903 4.22245330 2.38543817
C             0 3.20988023 0.80807706 0.02547700
C             0 3.43720324 1.56885111 1.17933309
C             0 4.37022931 2.59653619 1.16973909
C             0 5.08542737 2.88135021 0.00943900
C             0 4.86535235 2.13094416 -1.14186508
C             0 3.93099828 1.10266008 -1.13706308
C             0 2.26824216 -1.49037010 -1.43715410
C             0 3.41362825 -2.30885016 -1.48792711
C             0 2.63415519 -3.30895524 -3.54551525
......................................................................................
......................................................................................
......................................................................................
......................................................................................

作者
Author: pwzhou 时间: 2020-8-12 11:22
可以计算，因为是数值的，所以需要计算3N-6个TDDFT，N等于你的原子数，你可以统计一下算一个需要多长时间，你现在算了多少个了，剩下的还需要多长时间。

如果有16就用16吧，解析梯度要快很多很多。

作者
Author: zhficcas 时间: 2020-8-12 11:30

pwzhou 发表于 2020-8-12 11:22
可以计算，因为是数值的，所以需要计算3N-6个TDDFT，N等于你的原子数，你可以统计一下算一个需要多长时间， ...

哦，原来如此，谢谢。我87个原子喔，要算243个循环，现在才算到第五个! 16我也算过了，但我需要用DUSHIN软件来分析频率数据，只能用09版的，DUSHIN软件不支持16版的，谢谢！

作者
Author: pwzhou 时间: 2020-8-12 11:47

zhficcas 发表于 2020-8-12 11:30
哦，原来如此，谢谢。我87个原子喔，要算243个循环，现在才算到第五个! 16我也算过了，但我需要用DUSHIN ...

需要的循环次数我说错了，应该是6N次，你重新估算一下吧

作者
Author: ABetaCarw 时间: 2020-8-12 12:13

pwzhou 发表于 2020-8-12 11:47
需要的循环次数我说错了，应该是6N次，你重新估算一下吧

更绝望了

作者
Author: zhficcas 时间: 2020-8-12 13:47

pwzhou 发表于 2020-8-12 11:47
需要的循环次数我说错了，应该是6N次，你重新估算一下吧

天哪，没法算了！

作者
Author: hebrewsnabla 时间: 2020-8-12 14:21

zhficcas 发表于 2020-8-12 11:30
哦，原来如此，谢谢。我87个原子喔，要算243个循环，现在才算到第五个! 16我也算过了，但我需要用DUSHIN ...

dushin可以分析g16的结果，参见http://bbs.keinsci.com/thread-7574-1-1.html

作者
Author: sobereva 时间: 2020-8-13 08:14
用g09的半数值TDDFT做振动分析，计算一个中小分子通常就得一、两个礼拜。

作者
Author: zhficcas 时间: 2020-8-13 17:01

sobereva 发表于 2020-8-13 08:14
用g09的半数值TDDFT做振动分析，计算一个中小分子通常就得一、两个礼拜。

谢谢Sob老师关注留言，但如果要通过Dushin来分析S1态的频率，只能用09高斯，用16的有时候dushin识别不了。今天我尝试用16计算的频率，把关键数据拷贝到09文件里，到是可以解决用16计算频率，通过Dushin软件来分析。不过当用在限制性优化和相应的频率分析这种类型时，貌似dushin软件无法识别，正在尝试。

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