24小时热门版块排行榜    

查看: 829  |  回复: 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的回帖

mawenfeng415

金虫 (小有名气)

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

beefly

专家顾问 (职业作家)

地沟油冶炼专家

先保证casscf能量曲线是正确的,否则ci,ignore_error得到的都是垃圾
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
2楼2012-04-18 00:16:51
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 311求调剂 +6 南迦720 2026-02-28 6/300 2026-02-28 17:17 by 刘茉莉
[考研] 高分子化学与物理调剂 +3 好好好1233 2026-02-28 4/200 2026-02-28 17:13 by sandychj
[考研] 化工专硕348,一志愿985求调剂 +3 弗格个 2026-02-28 5/250 2026-02-28 17:04 by sandychj
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
[考研] 0856材料求调剂 +8 hyf hyf hyf 2026-02-28 9/450 2026-02-28 16:27 by etapple
[考研] 285求调剂 +3 满头大汗的学生 2026-02-28 3/150 2026-02-28 16:22 by 无际的草原
[考研] 312求调剂 +4 吃宵夜1 2026-02-28 5/250 2026-02-28 16:11 by gjm133
[考博] 26申博 +3 想申博! 2026-02-26 3/150 2026-02-28 16:07 by nxgogo
[考研] 材料类求调剂 +3 wana_kiko 2026-02-28 3/150 2026-02-28 15:03 by lature00
[考研] 295求调剂 +4 19171856320 2026-02-28 4/200 2026-02-28 13:39 by ms629
[考研] 290求调剂 +4 材料专硕调剂; 2026-02-28 5/250 2026-02-28 13:32 by houyaoxu
[考研] 寻找调剂 +3 LYidhsjabdj 2026-02-28 3/150 2026-02-28 12:59 by miniwendy
[考研] 304求调剂 +5 曼殊2266 2026-02-28 6/300 2026-02-28 12:44 by 迷糊CCPs
[硕博家园] 博士自荐 +6 科研狗111 2026-02-26 9/450 2026-02-28 12:32 by seaskyy
[考研] 272求调剂 +3 田智友 2026-02-28 3/150 2026-02-28 12:31 by 王加浩to
[考研] 298求调剂 +4 axyz3 2026-02-28 4/200 2026-02-28 11:21 by wang_dand
[基金申请] 什么是人一生最重要的? +10 瞬息宇宙 2026-02-21 10/500 2026-02-27 08:46 by tfang
[基金申请] 面上可以超过30页吧? +12 阿拉贡aragon 2026-02-22 13/650 2026-02-26 22:09 by Hahaxia
[教师之家] 版面费该交吗 +14 苹果在哪里 2026-02-22 17/850 2026-02-26 11:55 by AGI智创机器人
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
请填处理意见