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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Molpro/Molcas » 用MRCI计算CO+ 势能曲线以及跃迁振幅 的速度问题
51/1返回列表
查看: 1926  |  回复: 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的回帖

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的回帖
查看全部 14 个回答

beefly

专家顾问 (职业作家)

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

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

virtualzx

木虫 (著名写手)
你的活性空间奇怪。哪些轨道在里面？
mcscf和ci不连续未必是对称性的问题，也可能轨道或被平均的态发生改变了，这很常见

Molpro会通过块对角化重叠积分发现更高对称性的，但是mcscf时具体什么命令指定我就记不得了，应该是sym吧

发自小木虫IOS客户端
赞一下 回复此楼
5楼2015-12-31 22:45:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭