24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1502  |  回复: 1

henanlin

新虫 (初入文坛)

[求助] ?CORE ORBITAL RECORD 2140.2 CORRESPONDS TO DIFFERENT GEOMETRY急求遇到过相同问题

恳请各位大神支招
  ***,TiC molecule
  ***,TiC molecule
***,TiC molecule
       memory,5000,m
  
  
     rTiC=[1.3,1.31,1.33,1.35,1.402,1.443,1.475,1.503,1.542,1.583,1.615,1.641,1.673,1.702,1.712,1.723,1.732,1.738,1.743,1.748,1.752,1.756,1.760,1.763,1.766,1.770,1.775,1.780,1.785,1.790,1.795,1.801,1.806,1.812,1.818,1.826,1.842,1.863,1.885,1.902,2.003,2.201,2.405,2.906,3.432,3.823,4.512,5.123,5.745,6.324,7.0]
ang
  
   geometry={
       Ti
       C, Ti, r
       }

  do i=1,#rTiC
       r=rTiC(i)
  
      if(i.eq.1) then
         hf
  
        end if
  
  
      {casscf;frozen,6,2,2,0;occ,12,5,5,1;
  
         wf,28,1,2;
           wf,28,2,2;wf,28,3,2;
           maxit,500,500;
              }
       e1(i)=energy(1)
       dip1(i)=dmz(1)
        e2(i)=energy(2)
       dip2(i)=dmz(2)
        e3(i)=energy(3)
       dip3(i)=dmz(3)
  
  
      {mrci;core,6,2,2,0;occ,12,5,5,1;wf,28,3,2;save,2351.2;maxiter,500,500;orbital,ignore_error;}
        e1(i)=energy(1)
         ed1(i)=energd(1)
         dip1(i)=dmz(1)
  
       {mrci;core,6,2,2,0;occ,12,5,5,1;wf,28,1,2;save,2352.2;maxiter,500,500;orbital,ignore_error;}
        e2(i)=energy(1)
          ed2(i)=energd(1)
          dip2(i)=dmz(1)
  
     {ci;trans,2351.2,2352.2,dm;}
  
     t1dipx(i)=trdmx;
     t1dipy(i)=trdmy;
     t1dipz(i)=trdmz;
  
         table,rTiC,e1,dip1,e2,dip2,e3,dip3
        table,rTiC,e1,ed1,dip1,e2,ed2,dip2,t1dipx,t1dipy,t1dipz
       end do
基组跟着文献走的:ANO(Ti)/cc-pVQZ(C)  活性空间是3d4s4p(Ti)/2s2p(C)    第一个点能正常计算,到第二个点就报错:?CORE ORBITAL RECORD 2140.2 CORRESPONDS TO DIFFERENT GEOMETRY

后面参考MOLPRO使用手册:如果指定了NOCHECK,就会跳过一些寻找正确轨道的一致性检查,并且忽略“ORBITALS CORRESPOND TO DIFFERENT GEOMETRY”一类的错误信息。在CASSCF、MRCI中添加了:orbital,nocheck;        任然是一样的错误

再次修改了轨道:{casscf;frozen,4,1,1,0;closed,6,2,2,0;occ,12,4,4,1;          设置了闭壳层,还是报同样的错误

跪求各位大神指点啊!!!
回复此楼

» 猜你喜欢

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

850487336

新虫 (初入文坛)

geometry裡面r改成l把if(i.eq.1)刪去

发自小木虫Android客户端
2楼2017-05-25 20:05:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 henanlin 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然评分 +3 无名者登山 2026-05-20 4/200 2026-05-20 13:29 by vito刘
[基金申请] 提交了我也来说说感想 +4 fummck 2026-05-20 4/200 2026-05-20 12:56 by Equinoxhua
[基金申请] 面上本子正文33页,违规吗?会被低分嘛? +9 1234567wang 2026-05-17 11/550 2026-05-20 12:37 by zhuzg0628
[基金申请] 评审有感 +13 popular289 2026-05-18 22/1100 2026-05-20 11:28 by 水和泥不是水泥
[有机交流] 反应很差,大量原料没有反应 5+3 Mr.Zot 2026-05-19 6/300 2026-05-20 10:36 by xtlilibin
[有机交流] 如何实现卤原子转化 +3 BT20230424 2026-05-15 8/400 2026-05-20 09:32 by xtlilibin
[教师之家] 上海大学实验技术岗位非升即走 +12 嘻嘻哈哈乐呵呵 2026-05-15 13/650 2026-05-20 08:34 by xli1984
[文学芳草园] 献血感触 +6 呀呀好傻 2026-05-19 11/550 2026-05-19 22:26 by 而立得乐
[基金申请] 2026山东省优青 +3 luckinging 2026-05-14 4/200 2026-05-19 16:15 by 52wls
[考博] 26/27申博自荐-锂/钠电池方向 5+3 狗头军师. 2026-05-15 4/200 2026-05-19 09:10 by moonboat
[基金申请] 别被青基扩招骗了!26年科研内卷才刚刚开始 +3 国自然国社科中 2026-05-14 4/200 2026-05-19 08:48 by archvillain
[考博] 2026博士还有哪些学校有名额 +7 小王求读研 2026-05-15 8/400 2026-05-19 08:27 by zhyzzh
[考博] 博士申请 +5 星…… 2026-05-18 6/300 2026-05-18 23:49 by 糊糊涂涂好
[基金申请] 国自然上会要求 +5 无名者登山 2026-05-18 9/450 2026-05-18 17:50 by BlakeReary
[文学芳草园] 半夜喝咖啡 +3 myrtle 2026-05-15 5/250 2026-05-18 01:03 by 小沈2018
[考博] 光量子物理方向 博士招生 1人(2026.09) +3 sandyworld 2026-05-15 4/200 2026-05-17 14:38 by sandyworld
[有机交流] 求助2,4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 6/300 2026-05-16 19:46 by Equinoxhua
[硕博家园] 申请博士 +3 呃?呃 2026-05-15 3/150 2026-05-16 11:01 by a4742549
[考博] 26应届毕业生考博求助 +3 wo一定上岸 2026-05-13 3/150 2026-05-14 21:47 by 明海天涯
[考博] 材料类只有一篇综述能申博么 +4 乐逍遥谷 2026-05-13 4/200 2026-05-14 12:05 by zhyzzh
信息提示
请填处理意见