求助molpro的casscf计算C2分子的垂直激发能

原子分子物理 · 发表于 Post on 2020-12-4 10:37:38

本帖最后由原子分子物理于 2020-12-7 17:00 编辑

大家好，最近刚开始学习molpro的使用，想用molpro的casscf计算C2分子的垂直激发能，在论坛上也找了molpro manual和moloro的quickstart,刚刚熟悉了一些最简单的操作，quickstart上有一个甲醛的casscf 计算的例子，还有一个关于O2分子的态平均的计算例子（附图如下），我在关于C2分子的文献里看到C2分子的casscf计算是CASSCF(8,8)/aug-cc-pvtz,请问大家1.我可以模仿这两个例子计算C2分子的垂直激发能吗？2.怎么在输入卡中体现出活性空间是CAS(8,8)呢？如果这样不可以计算垂直激发能也请大家多多指教我应该怎么写出正确的输入卡，如果大家有学习molpro 更好的资料也可以分享给我，谢谢大家。

wzkchem5 · 发表于 Post on 2020-12-4 11:41:48

1.可以
2.需要仔细看文献，看文献选的8个轨道分别属于什么不可约表示，然后照你贴的图的示例写输入文件就行了

原子分子物理 · 发表于 Post on 2020-12-4 15:25:30

wzkchem5 发表于 2020-12-4 11:41
1.可以
2.需要仔细看文献，看文献选的8个轨道分别属于什么不可约表示，然后照你贴的图的示例写输入文件就 ...

好的，谢谢。我先尝试一下，有问题再来请教您

原子分子物理 · 发表于 Post on 2020-12-7 17:04:39

本帖最后由原子分子物理于 2020-12-7 17:09 编辑

wzkchem5 发表于 2020-12-4 11:41
1.可以
2.需要仔细看文献，看文献选的8个轨道分别属于什么不可约表示，然后照你贴的图的示例写输入文件就 ...

老师您好，我仿照文献中和quickstart中O2态平均的例子计算了C2的CASSCF，活性空间选取的是全部的价轨道，我想要的垂直激发能的数据在哪看呢（不好意思，刚刚入门），输入文件和输出文件已上传，麻烦老师指导一下，谢谢。

wzkchem5 · 发表于 Post on 2020-12-7 20:01:57

原子分子物理发表于 2020-12-7 17:04
老师您好，我仿照文献中和quickstart中O2态平均的例子计算了C2的CASSCF，活性空间选取的是全部的价轨道， ...

需要手动把激发态和基态的能量相减。比如可以从输出文件里看到
!MCSCF STATE 1.1 Energy -75.623333598242
!MCSCF STATE 1.2 Energy -75.555939556315
那么第一激发能就是(-75.555939556315) - (-75.623333598242) au

原子分子物理 · 发表于 Post on 2020-12-7 20:06:39

wzkchem5 发表于 2020-12-7 20:01
需要手动把激发态和基态的能量相减。比如可以从输出文件里看到
!MCSCF STATE 1.1 Energy ...

那激发态的电子态怎么判断啊，老师不好意思，初次上手，好多都不懂，我也没找到更多的资料

wzkchem5 · 发表于 Post on 2020-12-7 20:18:46

原子分子物理发表于 2020-12-7 20:06
那激发态的电子态怎么判断啊，老师不好意思，初次上手，好多都不懂，我也没找到更多的资料

可以从态的名字上看出来，比如STATE 1.2就是第2个不可约表示的第一个态。不同点群下各个不可约表示的编号在molpro手册里有。
另外，你现在每个不可约表示只算一个态的话，理论上只能知道基态的不可约表示，万一第一激发态和基态的不可约表示相同，你就把第一激发态漏掉了。除非文献已经告诉你，你需要的激发态每个不可约表示至多有一个，否则每个不可约表示应该多算几个态。

原子分子物理 · 发表于 Post on 2020-12-7 20:26:38

本帖最后由原子分子物理于 2020-12-7 20:28 编辑

wzkchem5 发表于 2020-12-7 20:18
可以从态的名字上看出来，比如STATE 1.2就是第2个不可约表示的第一个态。不同点群下各个不可约表示的编号 ...

好的，谢谢老师，因为还不是很熟练，所以就每个只算一个态试试，我再好好研究研究，尽快熟练起来，非常感谢老师的指导，另外，老师您有其他关于molpro学习的资料吗？我在网上找了好久，只找到了quickstart和manual.如果老师有其他与之相关的资料可以分享给我吗？

wzkchem5 · 发表于 Post on 2020-12-7 21:25:24

原子分子物理发表于 2020-12-7 20:26
好的，谢谢老师，因为还不是很熟练，所以就每个只算一个态试试，我再好好研究研究，尽快熟练起来，非常感 ...

我基本上也就是看manual，有的地方写得不是特别明晰，但是自己试着跑一两个输入文件就懂了

原子分子物理 · 发表于 Post on 2020-12-8 08:38:34

wzkchem5 发表于 2020-12-7 21:25
我基本上也就是看manual，有的地方写得不是特别明晰，但是自己试着跑一两个输入文件就懂了

好的，谢谢老师

，那我再研究研究手册吧

原子分子物理 · 发表于 Post on 2020-12-10 15:32:02

本帖最后由原子分子物理于 2020-12-10 16:46 编辑

原子分子物理发表于 2020-12-8 08:38
好的，谢谢老师，那我再研究研究手册吧

原子分子物理 · 发表于 Post on 2020-12-10 15:42:04

每个不可约表示下的两个态的计算的输出文件

原子分子物理 · 发表于 Post on 2020-12-10 16:46:07

wzkchem5 发表于 2020-12-7 21:25
我基本上也就是看manual，有的地方写得不是特别明晰，但是自己试着跑一两个输入文件就懂了

老师，不好意思又来打扰您，如您之前所说：!MCSCF STATE 1.1 Energy             -75.623333598242
!MCSCF STATE 1.2 Energy             -75.555939556315
那么第一激发能就是(-75.555939556315) - (-75.623333598242) au，那么第二激发能是不是就是!MCSCF STATE 1.3 Energy             -75.535139690849与!MCSCF STATE 1.1 Energy             -75.601109764769的差？以此类推，第三激发能是!MCSCF STATE 1.4 Energy             -75.495997872160与!MCSCF STATE 1.1 Energy             -75.601109764769的差？如果我每个不可约表示下了算了两个态，又该如何相减呢，请老师多多指教，谢谢老师

wzkchem5 · 发表于 Post on 2020-12-10 16:59:42

原子分子物理发表于 2020-12-10 16:46
老师，不好意思又来打扰您，如您之前所说：!MCSCF STATE 1.1 Energy -75.623333598242
!MC ...

“第一激发能就是(-75.555939556315) - (-75.623333598242) au”有个前提条件，就是你的STATE 1.2确实是第一激发态，而这不一定成立。第二到第四激发态同理。所以必须每个不可约表示多算几个态，然后把所有的能量从低到高排序，最低的是基态，第二低的是第一激发态，以此类推，然后把第n激发态的能量和基态能量相减，就是第n激发能。

原子分子物理 · 发表于 Post on 2020-12-10 17:34:06

wzkchem5 发表于 2020-12-10 16:59
“第一激发能就是(-75.555939556315) - (-75.623333598242) au”有个前提条件，就是你的STATE 1.2确实是 ...

老师您这里面说的基态能量是指state1.1吗？因为我每个不可约表示算了两个态，所以有state1.1和state1.2，那这两个哪个是基态能量呢?主要是这一点不是很懂

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

[Molpro] 求助molpro的casscf计算C2分子的垂直激发能

浏览过的版块