| 查看: 1387 | 回复: 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
» 猜你喜欢
招贤纳博（已结束）
已经有28人回复
-
表征测试机构后台私聊我
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有194人回复
-
需要一个圆偏振光诱导化学反应，请问使用圆偏振光作为光源，应该需要哪些设备搭建
已经有0人回复
-
Coordination Chemistry Reviews
已经有1人回复
-
请问现在还有电池材料方向的博导招人吗
已经有2人回复
» 本主题相关价值贴推荐，对您同样有帮助:
请教用molpro打印氢分子离子的重叠积矩阵、单电子积分矩阵、偶极矩阵
已经有3人回复- 【求助】Molpro如何执行旋轨耦合计算 如何执行态平均计算?
已经有10人回复
3楼2013-11-02 08:52:44
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +457 - QC强帖: 12
- 应助: 86 (初中生)
- 金币: 12717.9
- 散金: 27307
- 红花: 300
- 帖子: 4663
- 在线: 925.3小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
2楼2013-11-01 23:47:42
4楼2013-11-02 10:01:50
5楼2013-11-02 15:31:58