| 查看: 1538 | 回复: 4 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
molpro,casscf作多个态的平均，然后mrci要计算五个态的势能与偶极矩？
|
|||
|
大家好，我正学molpro，我看文献来计算CaO基态的势能曲线和偶极矩，但是最低的五个态势能曲线靠的很近，所以要用casscf来做态平均，对不？然后用mrci计算势能和偶极矩，对吗？但是有两个Σ态，都是单重的怎么区分? 下面是输入，哪有问题？ ***,CaO molecule memory,800,m rcao=[2.8,3.0,3.2,3.3,3.41,3.42,3.43,3.432,3.434,3.436,3.438,3.44,3.442,3.443,3.444,3.446,3.45,3.46,3.48,3.5,3.55,3.6,3.66,3.7,3.76,3.80,3.84,3.88,3.9,3.92,3.94,3.96,3.98,4.0,4.02,4.05,4.1,4.15,4.2,4.25,4.3,4.35,4.4,4.45,4.5,4.55,4.6,4.65,4.7,4.75,4.8,4.85,4.9,4.92,4.94,4.96,4.98,5.0,5.02,5.04,5.08,5.1,5.15,5.2,5.25,5.3,5.35,5.4,5.45,5.5,5.55,5.6,5.65,5.7,5.75,5.8,5.85,5.9,5.92,5.94,5.96,5.98,6.0,6.02,6.04,6.08,6.1,6.15,6.2,6.5,7.0,7.5,8.0,8.5,9.0,10] geometry={ O Ca, O, r } basis={ spdfgh,O,aug-cc-pv5z spdfgh,Ca,cc-pcv5z; } do i=1,#rcao r=rcao(i) if(i.eq.1) then hf end if {casscf;closed,7,2,2,0;occ,12,5,5,1; wf,28,1,0;state,2; wf,28,1,2; wf,28,2,0; wf,28,2,2; wf,28,3,0; wf,28,3,2; } ecasscf(i)=energy dipcasscf(i)=dmz {mrci;closed,6,2,2,0;occ,10,4,4,0; wf,28,1,0;state,2;weight,1,0;orbital,ignore_error; } e1(i)=energy dip1(i)=dmz {mrci;closed,6,2,2,0;occ,10,4,4,0; wf,28,1,0;state,2;weight,0,1;orbital,ignore_error; } e2(i)=energy dip2(i)=dmz table,rcao,ecasscf,dipcasscf table,rcao,e1,dip1, end do 附文献 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : CaO.pdf
2013-11-01 21:51:55, 565.37 K
» 猜你喜欢
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有132人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有1人回复
-
实际气体内压强的微观机制或将改写--------从一个基本的力学佯谬说起
已经有1人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有5人回复
-
招收0703开头的学硕调剂生 线上复试
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教用molpro打印氢分子离子的重叠积矩阵、单电子积分矩阵、偶极矩阵
已经有3人回复
- 【求助】Molpro如何执行旋轨耦合计算 如何执行态平均计算?
已经有10人回复
5楼2013-11-02 15:31:58
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +458 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13087.2
- 散金: 27307
- 红花: 300
- 帖子: 4673
- 在线: 927.3小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
2楼2013-11-01 23:47:42
3楼2013-11-02 08:52:44
4楼2013-11-02 10:01:50
5