应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Molpro/Molcas » 用MRCI计算CO+ 势能曲线以及跃迁振幅 的速度问题
51/1返回列表
查看: 1909  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

ideallzc

铜虫 (小有名气)

[求助] 用MRCI计算CO+ 势能曲线以及跃迁振幅 的速度问题

输入文件如下：

***COplus_potential
memory,100,m;
symmetry,x,y
geometry={c;o,c,r(i)}

!basis=aug-cc-pV5Z
basis=AVQZ
!basis=cc-pVQZ

distances=[1.5,1.6,1.7,1.8,1.9,2.0,2.1,2.2,2.3,2.4,2.5,2.6,2.65,2.7,2.75,2.8,2.9,3.0,3.2,3.4,3.6,3.8,3.9,4.0,4.2,4.6,5.0,5.5,6.0]
do i=1,#distances
r(i)=distances(i)

{hf;wf,13,1,1;}

{casscf; ! do casscf calculaion
orbital,ignore_error;
OCC,7,2,2,0; ! specify occupied space
!CLOSED,0,0,0,0; ! specify closed-shell (inactive) orbitals
FROZEN,2,0,0,0; ! specify Frozen core orbitals
WF,13,1,1;state,3;lquant,0,0,0; ! specifies number of electrons, point group symmetry and spin
WF,13,2,1;state,2;lquant,1,1;
} ! define wavefunction symmetry


{mrci;
orbital,ignore_error;
maxiter,100,100;
occ,7,2,2,0;core,2,0,0,0;
wf,13,1,1;state,6;}
sE1(i)=energy(1);
sE2(i)=energy(2);
sE3(i)=energy(3);
sE4(i)=energy(4);
sE5(i)=energy(5);
sE6(i)=energy(6);

enddo

{table,r,sE1,sE2,sE3,sE4,sE5,sE6
head,r,sE1,sE2,sE3,sE4,sE5,sE6
save,co+_pes.tab
title,MRCI PES Results for CO+, basis $basis
}

想问下有没有办法可以加快这个程序的计算速度，就是有没有什么优化的地方。 其次，我用MRCI计算A1对应的态的时候发现没法区分开 Σ和△这两个态，导致算的PES 从一个态跑到另一个态上。不知道有没有类似于家如lquant（CASSCF）的方法进行解决。  
还想问下有没有介绍CO+不同Σ跃迁矩阵的文章？

谢谢，回复有重谢~~
回复此楼

» 猜你喜欢
ideal-neil
1楼 2015-12-30 20:29:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luqing6879

木虫 (著名写手)
1) lquant most time does not work. Just give it up.
2) The states with same symmetry should not cross. Get confused of your above descriptions.
3) You may try to calculate the energies of the Delta state in symmetry 4, rather than symmetry 1, since they are degenerate.
赞一下 回复此楼
cishan.goodweb.cn/慈善点击，点滴帮助。
9楼2016-01-05 06:07:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答

beefly

专家顾问 (职业作家)

地沟油冶炼专家
加快速度的办法就是每个点或者n个点放入一个文件，分成几个任务算

区分Sigma、Delta应该在casscf这一步完成
赞一下(1人) 回复此楼
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
2楼2015-12-30 21:58:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ideallzc

铜虫 (小有名气)
引用回帖:
2楼: Originally posted by beefly at 2015-12-30 21:58:08
加快速度的办法就是每个点或者n个点放入一个文件，分成几个任务算

区分Sigma、Delta应该在casscf这一步完成

加快速度的方法好像可行，我可以试下。 区分Sigma和Delta，我是在casscf用lquant进行区分，但是在MRCI还是会出现 当sigma和delta两条势能曲线相交的时候，算出来的结果（因为按照能量排序），势能曲线并不相交，可能前半段是sigma轨道，交点之后是delta轨道，我不知道这种情况怎么处理。
赞一下 回复此楼
ideal-neil
3楼2015-12-31 09:20:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家
引用回帖:
3楼: Originally posted by ideallzc at 2015-12-31 09:20:29
加快速度的方法好像可行，我可以试下。 区分Sigma和Delta，我是在casscf用lquant进行区分，但是在MRCI还是会出现 当sigma和delta两条势能曲线相交的时候，算出来的结果（因为按照能量排序），势能曲线并不相交，可 ...

加上swap以后，ci结果按照casscf参考态的顺序排序
ci,swap;
赞一下 回复此楼
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
4楼2015-12-31 10:08:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭