24小时热门版块排行榜    

查看: 2228  |  回复: 13
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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
已阅   回复此楼   关注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这一步完成
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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 4/200 2026-09-20 07:24 by 9i7PCSZyMzot
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 4/200 2026-09-20 07:24 by 9i7PCSZyMzot
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 E4V9bALS54sj 2026-09-19 4/200 2026-09-20 07:12 by 9i7PCSZyMzot
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 E4V9bALS54sj 2026-09-19 4/200 2026-09-20 07:12 by 9i7PCSZyMzot
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 E4V9bALS54sj 2026-09-19 4/200 2026-09-20 07:01 by 9i7PCSZyMzot
[基金申请] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 E4V9bALS54sj 2026-09-19 4/200 2026-09-20 07:01 by 9i7PCSZyMzot
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 4/200 2026-09-20 06:55 by 9i7PCSZyMzot
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 3/150 2026-09-20 00:15 by 9i7PCSZyMzot
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-19 3/150 2026-09-20 00:03 by 9i7PCSZyMzot
[基金申请] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:39 by 9i7PCSZyMzot
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:27 by 9i7PCSZyMzot
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:27 by 9i7PCSZyMzot
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:03 by 9i7PCSZyMzot
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:03 by 9i7PCSZyMzot
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-18 7/350 2026-09-19 21:03 by 9i7PCSZyMzot
[教师之家] 现代”学阀”该如何界定 +3 iaeyuan 2026-09-19 3/150 2026-09-19 20:52 by 没昵称呀
[考博] 申博发邮件 +8 Lmengk 2026-09-15 9/450 2026-09-19 09:27 by 中间体超市
[有机交流] 各位大神,目前国内有哪些比较好用的逆合成软件? 20+3 zapen 2026-09-15 9/450 2026-09-18 18:17 by zapen
[考博] 上海工程技术大学激光智能制造课题组|2027级博士研究生招生公告 +7 水士口 2026-09-16 8/400 2026-09-16 13:47 by kkklens
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +8 两三岁ss 2026-09-14 8/400 2026-09-16 09:46 by 水士口
信息提示
请填处理意见