24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 857  |  回复: 3

mawenfeng415

金虫 (小有名气)

[求助] molpro 运行出错

为什么求出的曲线是震荡的?
输入文件如下:
       ***,LiCl
       memory,60,m
       gprint,basis,orbitals=4,civector
       geometry={Li;Cl,Li,rLiCl}
       rvec=[1.0,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.8,1.9,2.0,2.5,2.6,3.0,3.5,4.0,4.5,5.0,5.5,5.6,5.8,5.9,6.0,6.1,6.2,6.3,6.4,6.5,6.6,6.7,6.8,6.9,7.0,8.0,9.0,10.0,11.0,12.0,13.0,14.0,15.0,16.0,18.0] ang
       do i=1,#rvec
       rLiCl=rvec(i)
       basis=avqz
       hf;
       ehf(i)=energy(1)
       {multi;
       occ,8,3,3,0;closed,5,1,1,0;
       wf,20,3,0;Lquant,1;}
       e1cas(i)=energy(1);
       e2cas(i)=energy(2);
       {mrci,check=0;orbital,ignore_error;
       occ,8,3,3,0;closed,5,1,1,0;core,5,1,1,0;
       wf,20,3,0;option,nstati=4;}
        e1mrci(i)=energy(1);
        e1mrda(i)=energd(1);
       table,rvec,ehf,e1mrda
       save,LiCl.tab
       title,Results for LiCl,basis $basis
       enddo
输出结果如下:
   1.0   -466.1294593   -465.9816344
     1.1   -466.4132956   -466.2702298
     1.2   -466.6244971   -466.1455674
     1.3   -466.7753238   -466.2313596
     1.4   -466.8797649   -466.7443837
     1.5   -466.9501947   -466.8617340
     1.6   -466.9963508   -466.6741125
     1.7   -467.0254537   -466.9447424
     1.8   -467.0426934   -466.4603611
     1.9   -467.0517418   -466.9876306
     2.0   -467.0551767   -464.6108939
     2.5   -467.0353506   -467.0251917
     2.6   -467.0287454   -467.0237818
     3.0   -467.0023228   -467.0254933
     3.5   -466.9739480   -466.8886113
     4.0   -466.9518099   -466.1772926
     4.5   -466.9347448   -466.8553181
     5.0   -466.9214207   -466.9896837
     5.5   -466.9108053   -467.0269363
     5.6   -466.9089381   -465.9699087
     5.8   -466.9054203   -466.9320624
     5.9   -466.9037614   -467.0223690
     6.0   -466.9021639   -466.6611873
     6.1   -466.9006246   -466.9791568
     6.2   -466.8991402   -467.0138725
     6.3   -466.8977076   -467.0328475
     6.4   -466.8963242   -466.5475732
     6.5   -466.8949874   -466.4958467
     6.6   -466.8936947   -467.0358892
     6.7   -466.8924438   -467.0273590
     6.8   -466.8912328   -466.5344676
     6.9   -466.8900595   -467.0287771
     7.0   -466.8889223   -466.2926032
     8.0   -466.8792002   -467.0283780
     9.0   -466.8717214   -467.0094948
    10.0   -466.8657739   -464.6418484
    11.0   -466.8609245   -467.0285566
    12.0   -466.8568915   -465.3473465
    13.0   -466.8534841   -467.0279868
    14.0   -466.8505658   -467.0299572
    15.0   -466.8480386   -467.0231719
    16.0   -466.8458286   -467.0283384
    18.0   -466.8421472   -466.1336679
做出的图没有任何的物理意义,问题出在哪儿?应该怎么去设置?
请高人详细指点一下,感激不尽!
回复此楼

» 猜你喜欢

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

在生活中亮化自己
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

先保证casscf能量曲线是正确的,否则ci,ignore_error得到的都是垃圾
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
2楼2012-04-18 00:16:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mawenfeng415

金虫 (小有名气)

能不能给我说的更加详细一点呢?谢谢!
在生活中亮化自己
3楼2012-04-18 13:54:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mawenfeng415

金虫 (小有名气)

引用回帖:
2楼: Originally posted by beefly at 2012-04-18 00:16:51:
先保证casscf能量曲线是正确的,否则ci,ignore_error得到的都是垃圾

但是如果不用ci,ignore_error的话总是出错算不下去啊。
在生活中亮化自己
4楼2012-04-18 14:31:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 mawenfeng415 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料工程281还有调剂机会吗 +42 xaw. 2026-04-11 43/2150 2026-04-15 09:57 by ssssjqrr
[考研] 291分调剂 +8 上岸小莹加油 2026-04-09 9/450 2026-04-15 00:41 by zzzggc
[考研] 22408 312求调剂 +16 门路摸摸 2026-04-14 17/850 2026-04-14 23:59 by Xurambo2014
[考研] 药学305求调剂 +10 玛卡巴卡boom 2026-04-10 10/500 2026-04-14 15:55 by zs92450
[考研] 279求调剂 +12 张番茄不炒蛋 2026-04-11 12/600 2026-04-14 15:38 by zs92450
[考研] 考研求调剂 +6 ban班小七 2026-04-11 6/300 2026-04-14 14:06 by 哆啦A梦只是个梦
[考研] 农学0904 312求调剂 +4 Say Never 2026-04-11 4/200 2026-04-14 09:10 by zs92450
[考研] 297工科,求调剂? +13 河南农业大学-能 2026-04-12 13/650 2026-04-13 14:12 by dingyanbo1
[考研] 求调剂,一志愿材料科学与工程985,365分, +8 材化李可 2026-04-11 10/500 2026-04-12 08:42 by 852137818
[考研] 求调剂 +11 翩翩一书生 2026-04-09 11/550 2026-04-11 19:57 by 逆水乘风
[考研] 求调剂,一志愿大连理工大学354分 +5 雨声余生 2026-04-11 6/300 2026-04-11 16:12 by 雨声余生
[考研] 求调剂 +3 胃痉挛累了 2026-04-11 5/250 2026-04-11 14:13 by luhong1990
[考研] 085410 273分调剂 +4 X1999 2026-04-09 4/200 2026-04-11 13:05 by pies112
[考研] 0859,337求调剂 +4 研s. 2026-04-10 4/200 2026-04-11 11:34 by caotw2020
[考研] 本科211 工科085400 280分求调剂 可跨专业 +11 LZH(等待调剂中 2026-04-10 11/550 2026-04-11 08:39 by zhq0425
[考研] 284求调剂 +9 让我上岸吧阿西 2026-04-09 11/550 2026-04-10 19:18 by 靖jing
[考研] 青岛科技大学材料学院,环境学院调剂补录4月10日以前都可以 +3 1青科大。 2026-04-09 5/250 2026-04-10 09:58 by 翩翩一书生
[考研] 0703化学求调剂 +21 不知名的小卅 2026-04-08 21/1050 2026-04-09 18:55 by l_paradox
[考研] 086004 求调剂 309 +7 Yin DY 2026-04-08 7/350 2026-04-09 13:59 by Delta2012
[考研] 296求调剂 +3 汪!?! 2026-04-08 3/150 2026-04-08 22:00 by zhouyuwinner
信息提示
请填处理意见