求助计算对二苯硫醇分子拉曼光谱出错

CohenTang

我想用gaussian计算BDT分子的拉曼光谱，但是计算总是出错：No CIS Freq=Raman. Error termination via Lnk1e
输入文件为：
%chk=BDTO.chk
# freq=raman td cam-b3lyp/6-311g(d,p) scrf=(solvent=water)
geom=connectivity

BDT-R

0 1
C                  1.43190800    0.04006100   -0.00004800
C                  0.66291000    1.24734300   -0.00028600
C                 -0.70786500    1.23122100    0.00001000
C                 -1.40043500   -0.00198300    0.00006200
C                 -0.66772200   -1.21063000   -0.00018400
C                  0.70457000   -1.19251400   -0.00023400
H                  1.19487800    2.19030900   -0.00063800
H                 -1.26068100    2.16357800    0.00034100
H                 -1.19542800   -2.15772700   -0.00049800
H                  1.26598000   -2.11563400   -0.00027700
S                 -3.14775900   -0.11140500    0.00012000
S                  3.13691000    0.11119100    0.00012400
H                 -3.40019500    1.21314900    0.00016600
H                  3.42883500   -1.97123900    0.00108700

1 2 1.5 6 1.5 12 1.0
2 3 2.0 7 1.0
3 4 1.5 8 1.0
4 5 1.5 11 1.0
5 6 2.0 9 1.0
6 10 1.0
7
8
9
10
11 13 1.0
12 14 1.0
13


输出文件：
Gaussian 16:  ES64L-G16RevA.03 25-Dec-2016
                18-Aug-2021
******************************************
%chk=BDTO.chk
----------------------------------------------------------------------
# freq=raman td cam-b3lyp/6-311g(d,p) scrf=(solvent=water) geom=connec
tivity
----------------------------------------------------------------------
No CIS Freq=Raman.
Error termination via Lnk1e in /txk/g16/l1.exe at Wed Aug 18 11:30:53 2021.
Job cpu time:       0 days  0 hours  0 minutes  0.4 seconds.
Elapsed time:       0 days  0 hours  0 minutes  0.1 seconds.
File lengths (MBytes):  RWF=      6 Int=      0 D2E=      0 Chk=      8 Scr=      1


不知道怎么回事，求大神指导！！

CohenTang

snljty 发表于 2021-8-18 12:05
G16 B.01才开始支持CIS/TD-DFT下的激发态Raman。从你的描述我也看不出你想做激发态Raman的想法。
这个体系 ...

谢谢，明白了

snljty

G16 B.01才开始支持CIS/TD-DFT下的激发态Raman。从你的描述我也看不出你想做激发态Raman的想法。
这个体系没什么理由用CAM-B3LYP。用B3LYP即可。（考虑校正因子）
你的初始结构有问题，有个氢位置明显不对，至少先优化。
Raman光谱计算极化率导数原理上应该加弥散函数，不然结果误差很大。

CohenTang

已解决，method里面误选TD-scf，是用来算紫外的，去掉TD就好了