关于计算振动分辨的电子光谱

小范范1989 · 发表于 Post on 2014-11-17 07:43:22

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。

小范范1989 · 发表于 Post on 2014-11-17 07:46:39

这事模拟出来的，也和您的不太一样

sobereva · 发表于 Post on 2014-11-17 08:39:33

1 和这没必然关系。这种计算顺序纯粹是因为高斯程序这么要求的。
2 吸收的时候怎么写发射的时候还怎么写，计算顺序也完全一样，先激发态后基态。唯一区别就是S0计算的时候加上emission关键词即可。
3 做出来应该是这样，HWHM设的是10。0-0是重合的，其余部分存在较明显的左右对称。

4 初始坐标无所谓，反正还得做优化
5 多写上scrf=PCM即可
6 对

用帖子里的输入文件，严格按帖子里的步骤去做，肯定能做出帖子里的图，注意坐标轴范围。

小范范1989 · 发表于 Post on 2014-11-17 08:54:34

本帖最后由小范范1989 于 2015-3-25 11:23 编辑

sobereva 发表于 2014-11-17 08:39
1 和这没必然关系。这种计算顺序纯粹是因为高斯程序这么要求的。
2 吸收的时候怎么写发射的时候还怎么写， ...

谢谢sobereva，我重新计算试试，有问题还需要麻烦您。
ps，后来我知道了。我计算文件写错了，所以，一定要认真。。。

小范范1989 · 发表于 Post on 2015-5-4 21:43:40

sobereva 发表于 2014-11-17 08:39
1 和这没必然关系。这种计算顺序纯粹是因为高斯程序这么要求的。
2 吸收的时候怎么写发射的时候还怎么写， ...

sob老师好，这个帖子这么久了，我现在计算这个发射，想做一个振动分辨的荧光光谱，所以看了看这个发射，我计算了那个帖子中的例子，得到的发射光谱和你的不太一样，你是调节了其他的参数吗？吸收光谱我做出来符合的很好。
附件是我计算的激发态的log和基态的log，请老师过目

sobereva · 发表于 Post on 2015-5-5 01:48:02

我还写了readFCHT，最后写上SPECHWHM=10
当时用的输入输出文件：

小范范1989 · 发表于 Post on 2015-5-5 09:13:30

sobereva 发表于 2015-5-5 01:48
我还写了readFCHT，最后写上SPECHWHM=10
当时用的输入输出文件：

谢谢sob老师，已经弄出来了，老师休息的有点晚，哈哈哈。

janihyhubo · 发表于 Post on 2017-2-14 09:17:59

小范范1989 发表于 2014-11-17 08:54
谢谢sobereva，我重新计算试试，有问题还需要麻烦您。
ps，后来我知道了。我计算文件写错了，所以，一 ...

您好，我现在做这个也出现您上面图片所示的错误了，但一直没找到错在哪，您方便说一下您之前是怎么改正的吗？谢谢

小范范1989 · 发表于 Post on 2017-2-14 10:06:32

janihyhubo 发表于 2017-2-14 09:17
您好，我现在做这个也出现您上面图片所示的错误了，但一直没找到错在哪，您方便说一下您之前是怎么改正的 ...

我记得应该是我的输入文件弄错了，好久了，有点遗忘。但是看当时回复，应该是输入文件错误，
可能是不细心，完全按照sob老师的来，你看看哪个地方你应该是遗漏了。

janihyhubo · 发表于 Post on 2017-2-14 10:27:44

我就是直接复制的博文中的例子，但是做出来图还是您上边那个图，附件是输入输出文件，http://www.sobereva.com/223 链接是博文，当局者迷，可能旁观者更能清楚的看明白问题出在哪，麻烦您抽时间帮我看看哈，谢谢

janihyhubo · 发表于 Post on 2017-2-14 15:36:26

谢谢，我后来又测试了一下，问题出在作图，必须把Y的数据中的D改成E才行，origin不识别D

liaorongbao · 发表于 Post on 2019-7-27 17:22:35

本帖最后由 liaorongbao 于 2019-7-28 11:26 编辑

janihyhubo 发表于 2017-2-14 15:36
谢谢，我后来又测试了一下，问题出在作图，必须把Y的数据中的D改成E才行，origin不识别D

用您的anisole.log中出现的两列结果数据输入origin作图(并把第二列数据D改成E)。做出来的图和小范范1989给出来的第一个图一样的。发现似乎是绘图软件origin的纵坐标乱来，曲线与数据对不上。但不知为啥。调整了两三个小时，现在终于正常了。重新建立了一个origin的Project就好了。之前那个origin的Project是我复制后删除数据重新写上数据的，可能保留有记忆。

wangxm · 发表于 Post on 2021-6-22 19:22:53

你好，我现在也在做振动分辨的光谱计算。我看到很多人用FCclasses程序，初学有点吃力，看不太懂怎么操作。然后看了卢老师的帖子，比较简单。我想请教一下对于一些FC和HT等方面的考虑，高斯计算和FCclasses有什么区别吗？这些在高斯中如何体现

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] 关于计算振动分辨的电子光谱

浏览过的版块