高锰酸钾的UV计算

发条鸟

想请教一下，我想用gauss计算一下高锰酸钾的紫外吸收图谱，用来验证和标定之前相关体系的计算结果，但计算出来的结果总是和实验结果不符，期望是在525nm处有吸收峰，但是算出来的结果虽然在525nm处有激发波长，但是并没有出现吸收峰，也考虑了溶剂效应。希望各位老师抽空指教，谢谢
%mem=2000MB
%rwf=BecpKMnO4-spin1-td-dft.rwf
%nosave
%nprocshared=4
%chk=BecpKMnO4-spin1-td-dft.chk
#p b3lyp/genecp nosymm iop(5/13=1,2/16=3) td(nstate=10,direct) scrf(pcm,solvent=water)

Title Card Required

0 1
Mn                -0.83266700   -0.31001000    0.68450100
O                 -0.69412400    0.87191900    1.76036000
O                 -0.22562600    0.24253700   -0.69378200
O                 -0.11543700   -1.62265600    1.15771400
O                 -2.40231700   -0.53655300    0.44297100
K                 -2.19722400    2.18473100   -0.21569700

O K 0
6-31+g**
*********
Mn 0
lanl2dz
*********

Mn 0
lanl2dz

sobereva

2/16毫无意义
万恶的5/13=1更是切记不能加。这个最糟的关键词这两年简直就像瘟疫一样传播。

看看525处的振子强度是否精确为0，是的话可能是对称禁阻，这在考虑振动耦合时实际还是可以观测到的。

发条鸟

sobereva 发表于 2015-9-23 17:33
2/16毫无意义
万恶的5/13=1更是切记不能加。这个最糟的关键词这两年简直就像瘟疫一样传播。

不是精确为0
Excited State   1:      Singlet-?Sym    2.3606 eV  525.23 nm  f=0.0003  <S**2>=0.000
      31 -> 34         0.30435
      32 -> 35         0.30448
      33 -> 34         0.51152
      33 -> 35         0.21731
把这两个关键词去掉了又算一遍，还是不对。。

sobereva

用FCHT方式考虑振动耦合看看结果
振动分辨的电子光谱的计算
http://sobereva.com/223

发条鸟

sobereva 发表于 2015-9-23 20:42
用FCHT方式考虑振动耦合看看结果
振动分辨的电子光谱的计算
http://sobereva.com/223

真是麻烦sob老师了，昨天晚上按照帖子里面的方法进行了计算，激发态的振动分析优化好了，但是基态的报错了，麻烦sob老师再指导下怎么解决这个
这是两次的输入文件：
激发态的计算
%mem=2000MB
%rwf=BecpKMnO4-spin1-exc.rwf
%nosave
%nprocshared=4
%chk=BecpKMnO4-spin1-exc.chk
#p opt cis(root=1)/genecp freq=saveNM

Title Card Required

0 1
Mn                -0.63978200   -0.29525900    0.72417600
O                 -0.29015000    0.70076000    1.93294300
O                  0.19840300    0.11084700   -0.56501700
O                 -0.27314200   -1.79420300    1.15277700
O                 -2.42884686   -0.19152292    0.35352015
K                 -3.03387680    2.29934527   -0.46233183

O K 0
6-31+g**
*********
Mn 0
lanl2dz
*********

Mn 0
lanl2dz

基态的：

%mem=2000MB
%rwf=BecpKMnO4-spin1-base.rwf
%nosave
%nprocshared=4
%chk=BecpKMnO4-spin1-base.chk
#p opt mp2/genecp freq=FCHT nosymm

Title Card Required

0 1
Mn                -0.63978200   -0.29525900    0.72417600
O                 -0.29015000    0.70076000    1.93294300
O                  0.19840300    0.11084700   -0.56501700
O                 -0.27314200   -1.79420300    1.15277700
O                 -2.42884686   -0.19152292    0.35352015
K                 -3.03387680    2.29934527   -0.46233183

O K 0
6-31+g**
*********
Mn 0
lanl2dz
*********

Mn 0
lanl2dz

BecpKMnO4-spin1-exc.chk

这是报错输出中的最后几段：

Approx. of the electronic transition dipole moment: FCHT
Type of transition requested: ONE-PHOTON ABSORPTION

Data for initial state taken from current calculation.
Passed-in normal modes used.
Data for   final state taken from checkpoint file "BecpKMnO4-spin1-exc.chk"
Normal modes recovered from file.
Transition dipole derivatives not found.
Error termination via Lnk1e in /home/msc/g09/l716.exe at Wed Sep 23 21:22:14 2015.
Job cpu time:  0 days  0 hours 17 minutes  8.9 seconds.
File lengths (MBytes):  RWF=   2097 Int=      0 D2E=      0 Chk=      2 Scr=      1
                                                                         21063,2       Bot

sobereva

注意第四节第一段的末尾“G09里CIS和TDDFT都能用于FC，但是HT或FCHT只能用CIS，而且做CIS时必须用freq=numer，因为只有这样才会产生跃迁偶极矩导数信息。”

另外当前情况基态最好用HF来算

发条鸟

sobereva 发表于 2015-9-24 21:41
注意第四节第一段的末尾“G09里CIS和TDDFT都能用于FC，但是HT或FCHT只能用CIS，而且做CIS时必须用freq=nume ...

sober老师，
基态的输入文件
%mem=2000MB
%rwf=BecpKMnO4-spin1-base.rwf
%nosave
%nprocshared=4
%chk=BecpKMnO4-spin1-base.chk
#p opt HF/genecp freq=FCHT nosymm

Title Card Required

0 1
Mn                -0.63978200   -0.29525900    0.72417600
O                 -0.29015000    0.70076000    1.93294300
O                  0.19840300    0.11084700   -0.56501700
O                 -0.27314200   -1.79420300    1.15277700
O                 -2.42884686   -0.19152292    0.35352015
K                 -3.03387680    2.29934527   -0.46233183

O K 0
6-31+g**
*********
Mn 0
lanl2dz
*********

Mn 0
lanl2dz

BecpKMnO4-spin1-exc.chk


激发态的输入文件：
%mem=2000MB
%rwf=BecpKMnO4-spin1-exc.rwf
%nosave
%nprocshared=4
%chk=BecpKMnO4-spin1-exc.chk
#p opt cis(root=1)/genecp freq=numer

Title Card Required

0 1
Mn                -0.63978200   -0.29525900    0.72417600
O                 -0.29015000    0.70076000    1.93294300
O                  0.19840300    0.11084700   -0.56501700
O                 -0.27314200   -1.79420300    1.15277700
O                 -2.42884686   -0.19152292    0.35352015
K                 -3.03387680    2.29934527   -0.46233183

O K 0
6-31+g**
*********
Mn 0
lanl2dz
*********

Mn 0
lanl2dz

报错文件中最后几行是：
ERROR: The Franck-Condon factor corresponding to the overlap integral
  between both vibrational ground states is too small:
  |<0'|0">|^2 <  0.13152E-56
那sober老师，这个是不是说明这个体系不适合用FC方法来计算？

sobereva

激发态优化完了你得到的结构是什么样的？如果和基态优化的结构偏离很大就没法算了。