计算化学公社
标题: 关于计算振动分辨的电子光谱 [打印本页]
作者Author: 小范范1989 时间: 2014-11-17 07:43
标题: 关于计算振动分辨的电子光谱
Sobereva你好，
我想咨询一下您，就是我看你写过的一个帖子，振动分辨的电子光谱，最后帖子的末尾说：
可以计算发射光谱，通过添加emission。我学习了这个帖子，我想在问一下几个问题，
1：帖子中计算吸收光谱，首先是优化了激发态，同时保存了他的震动信息（通过freq=saveNM），然后计算了基态的震动，同时和第一激发态的耦合（freq=FC，saveNM）。
我的理解是：基态是处于v=0，吸收跃迁到第一激发态的v=0，v=1，v=2…..，所以先计算激发态的，再计算基态的。
2：我想计算从s1到s0的振动分辨的发射光谱，
我首先计算了基态的振动信息，# hf/6-31G* opt freq=saveNM 这样写对不对？
同时，我又计算了s1的振动和优化，# cis(root=1)/6-31G* opt freq=(FC,saveNM,emission)，在最后加上了s0的chk的路径。
我这样改写不知道对不对？我看高斯手册，里面的emission写的是freq=(readFC,FC,Emission
3：这是我做出来的‘emission振动光谱’，还请sob查阅
4：还有一个小问题，就是计算的时候，这些初始的坐标是哪来的？就是我们高斯view画出来保存的吗？我的理解是，首先我们先优化s0，得到坐标，为后面的计算使用。
这样，我们在计算激发态的时候，里面的坐标是s0的坐标。我们再计算s0和s1振动耦的时候，用的也是我们早优化出来的s0的坐标？
5：溶剂中的怎样计算？就是在命令行加上scrf=（）？？这个吗？我看高斯守则中，还写得是freq=(readfc,fc,emission,readFCHT)?不理解，还请sob指点。
6：如果说想计算电离光谱，那么就可以把激发态改为阳离子。也就是改写它的电荷数是吧？电离一个电子，在输入文件中的电荷数就写1，电离两个就写2？
通过看您的帖子还有您的回复，学到了很多，谢谢您。Thanks。
作者Author: 小范范1989 时间: 2014-11-17 07:46
这事模拟出来的，也和您的不太一样
作者Author: sobereva 时间: 2014-11-17 08:39
1 和这没必然关系。这种计算顺序纯粹是因为高斯程序这么要求的。
2 吸收的时候怎么写发射的时候还怎么写，计算顺序也完全一样，先激发态后基态。唯一区别就是S0计算的时候加上emission关键词即可。
3 做出来应该是这样，HWHM设的是10。0-0是重合的，其余部分存在较明显的左右对称。
(, 下载次数 Times of downloads: 44)
4 初始坐标无所谓，反正还得做优化
5 多写上scrf=PCM即可
6 对
用帖子里的输入文件，严格按帖子里的步骤去做，肯定能做出帖子里的图，注意坐标轴范围。
作者Author: 小范范1989 时间: 2014-11-17 08:54
本帖最后由 小范范1989 于 2015-3-25 11:23 编辑
谢谢sobereva，我重新计算试试，有问题还需要麻烦您。
ps，后来我知道了。我计算文件写错了，所以，一定要认真。。。
作者Author: 小范范1989 时间: 2015-5-4 21:43
sob老师好，这个帖子这么久了，我现在计算这个发射，想做一个振动分辨的荧光光谱，所以看了看这个发射，我计算了那个帖子中的例子，得到的发射光谱和你的不太一样，你是调节了其他的参数吗？吸收光谱我做出来符合的很好。
附件是我计算的激发态的log和基态的log，请老师过目
作者Author: sobereva 时间: 2015-5-5 01:48
我还写了readFCHT，最后写上SPECHWHM=10
当时用的输入输出文件：
作者Author: 小范范1989 时间: 2015-5-5 09:13
谢谢sob老师，已经弄出来了，老师休息的有点晚，哈哈哈。
作者Author: janihyhubo 时间: 2017-2-14 09:17
您好，我现在做这个也出现您上面图片所示的错误了，但一直没找到错在哪，您方便说一下您之前是怎么改正的吗？谢谢
作者Author: 小范范1989 时间: 2017-2-14 10:06
我记得应该是我的输入文件弄错了，好久了，有点遗忘。但是看当时回复，应该是输入文件错误，
可能是不细心，完全按照sob老师的来，你看看哪个地方你应该是遗漏了。
作者Author: janihyhubo 时间: 2017-2-14 10:27
我就是直接复制的博文中的例子，但是做出来图还是您上边那个图，附件是输入输出文件，http://www.sobereva.com/223 链接是博文，当局者迷，可能旁观者更能清楚的看明白问题出在哪，麻烦您抽时间帮我看看哈，谢谢
作者Author: janihyhubo 时间: 2017-2-14 15:36
谢谢，我后来又测试了一下，问题出在作图，必须把Y的数据中的D改成E才行，origin不识别D
作者Author: liaorongbao 时间: 2019-7-27 17:22
本帖最后由 liaorongbao 于 2019-7-28 11:26 编辑
用您的anisole.log中出现的两列结果数据输入origin作图(并把第二列数据D改成E)。做出来的图和小范范1989给出来的第一个图一样的。发现似乎是绘图软件origin的纵坐标乱来，曲线与数据对不上。但不知为啥。调整了两三个小时，现在终于正常了。重新建立了一个origin的Project就好了。之前那个origin的Project是我复制后删除数据重新写上数据的，可能保留有记忆。
作者Author: wangxm 时间: 2021-6-22 19:22
你好，我现在也在做振动分辨的光谱计算。我看到很多人用FCclasses程序，初学有点吃力，看不太懂怎么操作。然后看了卢老师的帖子，比较简单。我想请教一下对于一些FC和HT等方面的考虑，高斯计算和FCclasses有什么区别吗？这些在高斯中如何体现
| 欢迎光临 计算化学公社 (http://bbs.keinsci.com/) |
Powered by Discuz! X3.3 |