24小时热门版块排行榜

南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知（长期有效）

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

长风_破浪

捐助贵宾 (小有名气)

应助: 1 (幼儿园)
金币: 720.1
散金: 30
红花: 2
帖子: 139
在线: 134.2小时
虫号: 2499503
注册: 2013-06-07
性别: MM
专业: 光学

[求助] 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

» 猜你喜欢

华南师范大学先进光电子研究院电子纸团队诚聘2026级博士研究生已经有0人回复
QE利用声子计算的一系列lambda和Tc如何选已经有8人回复
物理化学论文润色/翻译怎么收费? 已经有97人回复
PbS量子点紫外吸收已经有0人回复
理论计算合作已经有2人回复
小木虫的论文辅导靠谱吗？有没有用过的同学？已经有1人回复
QE计算声子ph.out的Γ点出现虚频（-30cm-1）vasp计算没有已经有0人回复
40-50万/年，中国散裂中子源诚聘计算模拟方向博士后已经有81人回复
2026博士申请—化学能源方向，本人熟练分子动力学模拟的使用已经有4人回复
比利时鲁汶大学与国家留学基金委共同资助博士研究生CSC-KU Leuven PhD Scholarship 已经有0人回复
书籍求助：汽车市场营销理论与实务（电子版）——章小平已经有0人回复

» 本主题相关价值贴推荐，对您同样有帮助:

请教用molpro打印氢分子离子的重叠积矩阵、单电子积分矩阵、偶极矩阵已经有3人回复
【求助】Molpro如何执行旋轨耦合计算如何执行态平均计算? 已经有10人回复

1楼2013-11-01 21:53:10

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

长风_破浪

捐助贵宾 (小有名气)

应助: 1 (幼儿园)
金币: 720.1
散金: 30
红花: 2
帖子: 139
在线: 134.2小时
虫号: 2499503
注册: 2013-06-07
性别: MM
专业: 光学

引用回帖:

2楼: Originally posted by beefly at 2013-11-01 23:47:42
两个单重Σ+态，放在同一个mrci中计算就可以了。如果分开算虽然会更快，但是容易出其它问题

请问是这样？
  {mrci;closed,6,2,2,0;occ,10,4,4,0;
   wf,28,1,0;state,2;orbital,ignore_error;
      }
  e1(i)=energy
dip1(i)=dmz
那怎么分别输出它们两个的势能和偶极矩呢？
哦，前面的态平均对吗？是不是还要多加几个态来平均效果会好？活性空间对吗？
还有用casscf优化时怎么能输出每个态的势能呢？very thank you !
求帮忙