| 查看: 1484 | 回复: 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
» 猜你喜欢
求催化作用基础第三版 （甄开吉) 课件
已经有5人回复
-
26年材料/化学相关专业博士研究生招生 温州大学化学与材料工程学院(26年秋季入学)
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有273人回复
-
求助电催化剂制备，怎么在碳布表面生长PdCo合金
已经有0人回复
-
求助电催化剂制备，怎么在碳布表面生长PdCo合金
已经有1人回复
-
304求调剂
已经有3人回复
-
上海华东理工大学机械与动力工程学院先进材料与制造研究组招聘公告
已经有0人回复
-
山西大学资源与环境工程研究所化工分离/环境功能材料小组招收2026年博士生
已经有2人回复
-
304求调剂
已经有6人回复
-
304求调剂
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教用molpro打印氢分子离子的重叠积矩阵、单电子积分矩阵、偶极矩阵
已经有3人回复
- 【求助】Molpro如何执行旋轨耦合计算 如何执行态平均计算?
已经有10人回复
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +458 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13001.7
- 散金: 27307
- 红花: 300
- 帖子: 4672
- 在线: 926.9小时
- 虫号: 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楼2013-11-02 15:31:58
