计算化学公社

标题: 高锰酸钾的UV计算 [打印本页]

作者
Author: 发条鸟 时间: 2015-9-23 17:01
标题: 高锰酸钾的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

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

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

作者
Author: 发条鸟 时间: 2015-9-23 20:08

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
把这两个关键词去掉了又算一遍，还是不对。。

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

作者
Author: 发条鸟 时间: 2015-9-24 09:24

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

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

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

作者
Author: 发条鸟 时间: 2015-9-25 12:54

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方法来计算？

作者
Author: sobereva 时间: 2015-9-25 18:11
激发态优化完了你得到的结构是什么样的？如果和基态优化的结构偏离很大就没法算了。

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