如何利用Gaussian16计算分子共振拉曼光谱（Resonance Raman Spectra）？

我爱问路

如题，本人欲采用Gaussian16计算一个分子的共振拉曼光谱（Resonance Raman Spectra），但计算一直无法成功，具体计算步骤、输入文件及错误如下。

第一步，我优化了基态结构，并做了Freq计算，输入文件信息如下：

1. %chk=gs_freq.chk
   
2. %Mem=32gb
   
3. %nproc=64
   
4. # CAM-B3LYP/def2SVP Freq=(raman,savenormalmodes) EmpiricalDispersion=GD3 NOSYMM
   
5. 
   
6. gs freq for C1Hdd
   
7. 
   
8. 0 1

复制代码

分子结构收敛，且无虚频。

第二步，我优化了第一激发态结构，并坐做Freq计算，输入文件信息如下：

1. %chk=s1_freq.chk
   
2. %Mem=32gb
   
3. %nproc=64
   
4. # CAM-B3LYP/def2SVP TD=(NStates=5,Root=1) Opt=CalcFC Freq=(raman,savenormalmodes) NoSymm
   
5. 
   
6. s1 opt and freq for C1Hdd
   
7. 
   
8. 0 1

复制代码

分子第一激发态结构收敛，且无虚频。

第三步，计算共振拉曼光谱。问题就出在此步。
参考Gaussian16官网的范例，我尝试了两个不同的输入文件，但计算均未成功。
1）第一种输入文件如下

1. %Chk=gs_freq
   
2. # Freq=(FC,ReadFC,ReadFCHT) Geom=AllCheck
   
3. 
   
4. TimeIndependent
   
5. Spectroscopy=ResonanceRaman
   
6. Spectrum=(Lower=800.,Upper=1800.,Broadening=Stick)
   
7. Intermediate=Source=Chk
   
8. RR=(OmegaMin=17400.,OmegaMax=17500.,OmegaStep=100.)
   
9. 
   
10. s1_freq.chk
    

复制代码

出现如下错误

1. ERROR: The Franck-Condon factor corresponding to the overlap integral
   
2.         between both vibrational ground states is too small:
   
3.          |<0'|0">|^2 =  0.28972E-08
   
4. Error termination via Lnk1e in /cm/shared/uniol/sw/SYSTEM/gaussian/g16.c01/g16/l718.exe at Wed Oct 18 12:02:24 2023.

复制代码

2）第二种输入文件如下

1. %Chk=gs_freq
   
2. # Freq=(FC,ReadFC,ReadFCHT) Geom=Check
   
3. 
   
4. RR spectrum
   
5. 
   
6. 0 1
   
7. 
   
8. Spectroscopy=RR
   
9. TransProp=EDip=Input
   
10. RR=OmegaList
    
11. TD=(2NSTEP=12,2NSTEPWIN=12,Time=1.0d-12,GauHWHM=100)
    
12. Print=(Tensors,Matrix=JK)
    
13. 
    
14. s1_freq.chk
    
15. 
    
16. 1.000D0 1.000D0 0.00D0
    
17. 
    
18. 17400. 17500.
    

复制代码

出现的错误如下

1. ERROR: End-of-file reached while reading list of incident frequencies
   
2. Error termination via Lnk1e in /cm/shared/uniol/sw/SYSTEM/gaussian/g16.c01/g16/l718.exe at Wed Oct 18 12:09:37 2023.

复制代码

我知道应该是我的输入文件有问题，但搜遍全网都没找到修正的方法或者相关的经验，故在此发帖，向各位请教，非常期望能得到大佬的指导和帮助！
不胜感激

liurui

我刚刚仔细研究了说明书，对于第二种输入，貌似是你的倒数两行有问题哇，跃迁电偶极矩，倒数第二行是三个数，倒数第一行的入射波波数两个，这两个应该是相互对应的，然后geom=allcheck 试试呢？

liurui

亲亲，我也遇到了这个问题，想问解决了吗