求助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 更好的资料也可以分享给我，谢谢大家。

zjxitcc · 发表于 Post on 2020-12-17 20:57:39

陈小北发表于 2020-12-17 10:19
啊是这样

但是老师的链接跳转到的帖子需要权限40才能看，也就是level 4, 你是很久之前下的吗？那时候 ...

quickstart.zip (263.97 KB, 下载次数 Times of downloads: 7)

陈小北 · 发表于 Post on 2020-12-17 20:42:17

原子分子物理发表于 2020-12-17 10:38
不好意思，我记错了，我是在小*虫下的，你可以去那里下，我现在没有那个链接了，你可以去搜一下

好的，我去找一下

谢谢！

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

陈小北发表于 2020-12-17 10:19
啊是这样

但是老师的链接跳转到的帖子需要权限40才能看，也就是level 4, 你是很久之前下的吗？那时候 ...

不好意思，我记错了，我是在小*虫下的，你可以去那里下，我现在没有那个链接了，你可以去搜一下

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

陈小北发表于 2020-12-17 10:19
啊是这样

但是老师的链接跳转到的帖子需要权限40才能看，也就是level 4, 你是很久之前下的吗？那时候 ...

这样啊，那我给你发一份吧，我也是刚学MOLPRO，半个月前刚下的，可能不是从老师这个链接下的

陈小北 · 发表于 Post on 2020-12-17 10:19:31

原子分子物理发表于 2020-12-17 09:11
我也是在这个论坛上找的，下面有老师发的链接，你可以下载一下。

啊是这样

但是老师的链接跳转到的帖子需要权限40才能看，也就是level 4, 你是很久之前下的吗？那时候这个帖子还没有被锁吗？

原子分子物理 · 发表于 Post on 2020-12-17 09:11:50

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

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

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

zjxitcc · 发表于 Post on 2020-12-16 21:21:03

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

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

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

陈小北 · 发表于 Post on 2020-12-16 20:57:29

楼主好！

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

非常感谢！！

zjxitcc · 发表于 Post on 2020-12-16 13:04:45

本帖最后由 zjxitcc 于 2020-12-16 13:30 编辑

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

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

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

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

原子分子物理 · 发表于 Post on 2020-12-16 11:36:18

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

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

zjxitcc · 发表于 Post on 2020-12-15 19:23:07

本帖最后由 zjxitcc 于 2020-12-16 13:31 编辑

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

%mem=8GB
%nprocshared=8
#p CASSCF(8,8)/aug-cc-pVTZ
mokit{casscf_prog=molpro}
0 3
C 0.0 0.0 0.0
C 0.0 0.0 2.0

复制代码

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

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

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

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

原子分子物理 · 发表于 Post on 2020-12-11 15:03:27

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

感谢，感谢。

pwzhou · 发表于 Post on 2020-12-11 12:58:44

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

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

jiangning198511 · 发表于 Post on 2020-12-11 11:36:48

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

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

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

浏览过的版块