计算化学公社

标题: 如何利用Gaussian16计算分子共振拉曼光谱（Resonance Raman Spectra）？ [打印本页]
作者
Author: 我爱问路    时间: 2023-10-18 19:05
标题: 如何利用Gaussian16计算分子共振拉曼光谱（Resonance Raman Spectra）？
本帖最后由 我爱问路 于 2023-10-18 19:07 编辑

如题，本人欲采用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

  6. gs freq for C1Hdd

  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

  6. s1 opt and freq for C1Hdd

  8. 0 1
复制代码
分子第一激发态结构收敛，且无虚频。

第三步，计算共振拉曼光谱。问题就出在此步。
参考Gaussian16官网的范例，我尝试了两个不同的输入文件，但计算均未成功。
1）第一种输入文件如下
  1. %Chk=gs_freq
  2. # Freq=(FC,ReadFC,ReadFCHT) Geom=AllCheck

  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.)

  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

  4. RR spectrum

  6. 0 1

  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)

  14. s1_freq.chk

  16. 1.000D0 1.000D0 0.00D0

  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.
复制代码


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


作者
Author: liurui    时间: 2023-12-28 16:27
亲亲，我也遇到了这个问题，想问解决了吗
作者
Author: liurui    时间: 2023-12-28 17:20
我刚刚仔细研究了说明书，对于第二种输入，貌似是你的倒数两行有问题哇，跃迁电偶极矩，倒数第二行是三个数，倒数第一行的入射波波数两个，这两个应该是相互对应的，然后geom=allcheck 试试呢？



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3