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

24小时热门版块排行榜    

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

virtualzx

木虫 (著名写手)
当你扫描势能曲线时，轨道特性会发生不连续改变。设计或许空间时不只需要知道轨道不可约表象，还需要知道每个轨道具体的电子结构。比如你的sigma轨道数超出价层轨道总数了，最高的几个轨道就非常容易改变。

不管你如何限制对称性，扫描全局势能面时不连续几乎是无法避免的。

例如你如果平均3个态，那么一旦第三和第四态相交时就会出现不连续。这只有通过动态权重调整才能实现。

轨道特征在某个点发生突变也很常见。例如一般键长较短时，低激发态会有rydberg state，一旦键长拉长它们就会被其他同对称性轨道取代。如果可能的话，活性空间应该只用来描述静态相关，动态相关应该由CI部分实现。也就是尽量不要在活性空间里包含不被任何一个被平均的态使用的轨道。如果有占据数接近零的轨道，它们会非常容易发生突变的

发自小木虫IOS客户端
赞一下(1人) 回复此楼
8楼2016-01-05 01:53:50
已阅   回复此楼   关注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 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭