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

wzkchem5

原子分子物理 发表于 2020-12-10 17:34
老师您这里面说的基态能量是指state1.1吗？因为我每个不可约表示算了两个态，所以有state1.1和state1.2， ...

基态永远是能量最低的那个态，你把所有态的能量比一下就知道了。一般而言，基态可能但不一定是state 1.1，应该完全按能量的高低来判断

原子分子物理

wzkchem5 发表于 2020-12-10 18:23
基态永远是能量最低的那个态，你把所有态的能量比一下就知道了。一般而言，基态可能但不一定是state 1.1 ...

好的老师，我想我这次是真的理解了，真的谢谢老师，特别感谢您。

原子分子物理

wzkchem5 发表于 2020-12-10 18:23
基态永远是能量最低的那个态，你把所有态的能量比一下就知道了。一般而言，基态可能但不一定是state 1.1 ...

老师，不好意思，还有一个问题不是很清楚，文献中对C2分子的活化空间有的是CAS(8,8),有的是CAS(4,4),我之前算的是选的全部的价轨道，完整的输入卡是这样的：*** , C2
print ,basis , orbitals
geometry ={
c1
c2 ,c1 ,r
}

r=2.3 bohr
basis =avdz

{hf
wf ,12 ,1 ,0  
occ ,2,1,1,0,2,0,0,0
}

{ casscf
wf ,12 ,1 ,0 ;state ,1
wf ,12 ,2 ,0 ;state ,2
wf ,12 ,3 ,0 ;state ,2
wf ,12 ,4 ,0 ;state ,2
wf ,12 ,5 ,0 ;state ,2
wf ,12 ,6 ,0 ;state ,2
wf ,12 ,7 ,0 ;state ,2}
如果我想用CAS(4,4)这个活化空间，occ和closed应该怎么写啊，怎么在输入卡上体现出我用的活化空间是CAS(4,4)呢？MOLPRO手册上的态平均的例子都是选的全部价轨道，所以这个地方就有点疑问，请老师多多指教，感谢。

wzkchem5

原子分子物理 发表于 2020-12-10 21:39
老师，不好意思，还有一个问题不是很清楚，文献中对C2分子的活化空间有的是CAS(8,8),有的是CAS(4,4),我之 ...

非态平均的例子应该有不选全部价轨道的例子吧？你参考非态平均的例子，再加上state, xxx指定态平均就行了

原子分子物理

wzkchem5 发表于 2020-12-11 11:03
非态平均的例子应该有不选全部价轨道的例子吧？你参考非态平均的例子，再加上state, xxx指定态平均就行了

首先感谢老师一直很耐心的解答我的疑问。有这样的例子，我已经尝试过CAS(6,6) ，但是算出来的的数值比文献中的数值偏高，老师我还想问一下，我想用CAS(4,4),但是需要交换轨道比如说：Occupied  (SGG) (SGU) (SGG) (SGU) (PIU) (PIU)
                                                      Virtual   (SGG) (PIG) (PIG) (SGU) (SGG) (PIU) (PIU) (SGG)
在这里面我想选用四个PI轨道，即PIU PIU PIG PIG 所以需要把空轨道中的SGG和PIG交换，我查到的资料是，比如：rotate,1,4,5.4,90，后面这个为什么是90，我在其他地方看到这个数值有写0的,这是有什么规则吗？还有就是这个命令应该是写在哪里啊，是写在{hf}这一部分还是{casscf}?

jiangning198511

MOLCAS手册里面有C2的计算例子，虽然两个软件有所不同，但你仔细分析一下相关内容估计会对你理解计算有所帮助
http://www.molcas.org/documentat ... 6112000000000000000

pwzhou

原子分子物理 发表于 2020-12-11 11:29
首先感谢老师一直很耐心的解答我的疑问。有这样的例子，我已经尝试过CAS(6,6) ，但是算出来的的数值比文 ...

交换轨道是90
交换需要写在casscf部分，因为是通过读取hf轨道并对其进行交换

原子分子物理

jiangning198511 发表于 2020-12-11 11:36
MOLCAS手册里面有C2的计算例子，虽然两个软件有所不同，但你仔细分析一下相关内容估计会对你理解计算有所帮 ...

感谢，感谢。

原子分子物理

pwzhou 发表于 2020-12-11 12:58
交换轨道是90
交换需要写在casscf部分，因为是通过读取hf轨道并对其进行交换

感谢老师，我也猜测是写在CASSCF部分，虽然提交程序没有报错，但是心里不确定到底对不对，我的课题跟师兄师姐们的不一样，他们也不用molpro,所以都是自己摸索猜测，特别感谢老师这几天一直耐心的指导，谢谢老师。

zjxitcc

担心自己结果不对、不会交换轨道、不知道如何选活性轨道，可以用 自动做多参考计算的程序MOKIT (https://gitlab.com/jxzou/mokit) 试算一下，以三重态C2为例，输入文件如下

1. %mem=8GB
   
2. %nprocshared=8
   
3. #p CASSCF(8,8)/aug-cc-pVTZ
   
4. 
   
5. mokit{casscf_prog=molpro}
   
6. 
   
7. 0 3
   
8. C   0.0   0.0   0.0
   
9. C   0.0   0.0   2.0
   
10. 
    

复制代码

整个过程是自动的，不用你看、选轨道，只要你机子上装有Gaussian, PySCF（免费开源）, GAMESS（免费开源）, Molpro等软件，（安装MOKIT后）会自动调用相应模块，8核8GB下算这个例子仅需23s时间。而人工手动一步步来耗时长多了。结果是
E(CASSCF) =       -75.48488843 a.u.
仅供参考。当然了，这只是基态的单点，但你可以用算完的molpro文件读取之、做后续计算，就不用再看、调轨道了。

原子分子物理

zjxitcc 发表于 2020-12-15 19:23
担心自己结果不对、不会交换轨道、不知道如何选活性轨道，可以用 自动做多参考计算的程序MOKIT (https://gi ...

MOLPRO在服务器里也可以调用吗？

zjxitcc

原子分子物理 发表于 2020-12-16 11:36
MOLPRO在服务器里也可以调用吗？

如果你指的是 在主节点提交，作业调度系统（如bsub, qsub等）会将其分配到子节点计算 这种形式的话，你可以在主节点提交MOKIT作业，MOKIT任务会被分配到子节点，它会在子节点里调用molpro，不用你做啥。

如果你指的服务器没有作业调度系统的话、当前节点直接可以计算，那更简单了，你只要能执行which molpro查看到molpro路径，那MOKIT自然也能调用。

另外，MOKIT没有用对称性，如果你要用对称性，可能要考虑读取算完的molpro输出文件、波函数文件等开了对称性再续算，或者自己仍手动算、不用MOKIT。

陈小北

楼主好！

想问一下您的moloro的quickstart可以分享一下吗？最近也在学molpro,想要做一些简单的练手

非常感谢！！

zjxitcc

陈小北 发表于 2020-12-16 20:57
楼主好！

想问一下您的moloro的quickstart可以分享一下吗？最近也在学molpro,想要做一些简单的练手

quickstart.pdf在解压后的doc文件夹中，可自行下载
http://bbs.keinsci.com/thread-15060-1-1.html

原子分子物理

陈小北 发表于 2020-12-16 20:57
楼主好！

想问一下您的moloro的quickstart可以分享一下吗？最近也在学molpro,想要做一些简单的练手

我也是在这个论坛上找的，下面有老师发的链接，你可以下载一下。