24小时热门版块排行榜    

查看: 841  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 环境工程调剂 +9 大可digkids 2026-03-16 9/450 2026-03-20 17:38 by 醉在风里
[考研] 287求调剂 +4 晨昏线与星海 2026-03-19 5/250 2026-03-20 16:42 by Iveryant
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +21 rare12345 2026-03-18 21/1050 2026-03-20 14:31 by 无懈可击111
[考博] 招收博士1-2人 +3 QGZDSYS 2026-03-18 3/150 2026-03-20 11:58 by 呱呱呱呱叫
[考研] 085700资源与环境308求调剂 +10 墨墨漠 2026-03-18 11/550 2026-03-20 11:57 by 学员8dgXkO
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 14/700 2026-03-20 10:46 by EBSD
[考研] 材料学硕318求调剂 +5 February_Feb 2026-03-19 5/250 2026-03-19 23:51 by 23Postgrad
[考研] 0703化学调剂 +10 妮妮ninicgb 2026-03-15 14/700 2026-03-19 22:59 by 学员8dgXkO
[考研] 307求调剂 +9 冷笙123 2026-03-17 9/450 2026-03-19 22:44 by 学员8dgXkO
[考研] 一志愿天大材料与化工(085600)总分338 +5 蔡大美女 2026-03-13 5/250 2026-03-19 10:44 by 是小刘呀~
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +10 Liwangman 2026-03-15 10/500 2026-03-19 10:25 by 无际的草原
[考研] 330求调剂 +3 小材化本科 2026-03-18 3/150 2026-03-18 21:55 by 无懈可击111
[考研] 311求调剂 +4 冬十三 2026-03-18 4/200 2026-03-18 21:47 by 尽舜尧1
[考研] 一志愿武理材料305分求调剂 +5 想上岸的鲤鱼 2026-03-18 6/300 2026-03-18 17:53 by 无际的草原
[考研] 311求调剂 +11 冬十三 2026-03-15 12/600 2026-03-18 14:36 by 星空星月
[考研] 0854可跨调剂,一作一项核心论文五项专利,省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 297求调剂 +8 戏精丹丹丹 2026-03-17 8/400 2026-03-18 14:30 by laoshidan
[考研] 一志愿,福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[基金申请] 今年的国基金是打分制吗? 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 321求调剂 +5 大米饭! 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
信息提示
请填处理意见