24小时热门版块排行榜    

查看: 2955  |  回复: 12

tiantian1358

新虫 (初入文坛)

[交流] 【求助】molpro 计算问题,请高手指教!!!! 已有2人参与

ERROR DETECTED
                                                  **************

     NON-ZERO P-SPACE GRADIENT


NON-ZERO P-SPACE GRADIENT. LARGEST VALUE:    0.93D+02
Probably the p-space is not closed under symmetry operations. This can happen if
the initial orbitals are not symmetry adapted or if there is root-flipping between states of
different symmetry. Try different p-space threshold, or select p-space configurations manually.
The check can be bypassed by setting thresh,pspace=thr to a larger value.
Non-zero p-space gradient elements smaller than thr are ignored.
molpro 计算时出现错误,请大家帮忙!!!!怎样才能解决这个问题,p-space这个关键词怎么用?请大家给些宝贵意见,谢谢,很急切!!!!
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

beefly

专家顾问 (职业作家)

地沟油冶炼专家


heyo_123(金币+1):鼓励交流!! 2010-07-10 18:36:53
第一次casscf计算的初始轨道不好,因为键长从HF的1.06跳到9.0。
你可以在HF用9.0,或者rvec=[1.06,1.5,2.0,3.0,...,9.0]

另外,默认的活性空间很大,你可能需要自己定义活性空间
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
7楼2010-07-08 23:36:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家


tiantian1358(金币+2):金币不多,请原谅,希望得到你的帮助 2010-07-07 22:38:41
heyo_123(金币+1):鼓励交流!! 2010-07-10 18:36:31
通常有两个原因:

活性空间设置不合理。

来自HF或CAS的初始轨道有问题。可能是没收敛。HF还可以试试加减电荷变成闭壳层,CAS可以试试态平均。
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
2楼2010-07-05 22:59:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tiantian1358

新虫 (初入文坛)

能不能具体说一下,加什么关键词什么的来修正呢,活性空间怎么设置?初始轨道有问题怎么解决呢,现在问题还是没有解决,一头雾水,希望能详细指点,非常感谢!!!!
3楼2010-07-07 22:37:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

先描述一下你的分子吧
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
4楼2010-07-08 01:53:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tiantian1358

新虫 (初入文坛)


heyo_123(金币+1):鼓励交流!! 2010-07-10 18:36:42
我的分子是N2O,分子基态为3重态对称性取的A“ 我对该分子基态对称性也有点不确定,看了些文献说是A“, 也有的是A'我的输入文件是这样的
***NN-rNN
    memory,100,m;
    symmetry,x;
    geometry={O;
              q,O,R1
              N1,q,r/2,O,zeta
              N2,q,r/2,N1,180.0,O,0}


    r=1.064 bohr
    R1=20.0 bohr
    zeta=90.0 deg
    l=0

    basis=VDZ
    {hf; wf,22,2,2;}


    rvec=[9.0,8.0,7.0,6.0] bohr

   do i=1,#rvec

      r=rvec(i)


      {casscf; wf,22,2,2;maxiter,200,200;orbital,ignore_error;}
      ecas(i)=energy


      rhh(i)=r
   table,rhh,ecas

   enddo

   table,rhh,ecas
5楼2010-07-08 21:40:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tiantian1358

新虫 (初入文坛)

真的非常感谢您的帮助,我有没有其他方式可以联系你,QQ或这邮箱?非常感谢
6楼2010-07-08 21:42:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tiantian1358

新虫 (初入文坛)

请问活性空间怎么设定呢?我的初始键长是1.06,我从1.0一直递增到10.0,一开始还给算,算到3.0左右的时候就出现上面的错误了,怎么办呢?后面的那些点也需要算啊,真的是太麻烦你了,十分感谢!!!
8楼2010-07-10 01:38:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家


heyo_123(金币+1):鼓励交流!! 2010-07-10 18:37:02
说明有相同对称性的电子态交叉,需要做态平均计算。

此外,molpro默认用12个活性轨道(3个2s+3个2p),可以试试把2s放到closed空间,计算会快得多。可以查看一下hf轨道能量,把3个2s挑出来。
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
9楼2010-07-10 11:57:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tiantian1358

新虫 (初入文坛)

太感谢你了,我是初学者,你的话好像有点听不懂,我很想和你交朋友,以后也许有很多问题我们可以直接探讨,不知道有没有这个荣幸,我的邮箱lijing52316@163.com  如果看到的话希望你给我发邮件,
10楼2010-07-10 17:35:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tiantian1358 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ero8OE6tv9cu 2026-09-02 4/200 2026-09-03 03:14 by rM1TE0WVDIIY
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 ero8OE6tv9cu 2026-09-02 6/300 2026-09-02 23:59 by rM1TE0WVDIIY
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +14 Kittylucky 2026-08-27 15/750 2026-09-02 19:35 by liumei0601
[基金申请] 科研人应该花精力去思考如何解决问题,而不是去凝练问题 +7 瞬息宇宙 2026-09-01 14/700 2026-09-02 17:33 by ma0526
[基金申请] 要骂人了,新模版改版就是要淡化问题凝练这种虚的东西,结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 10/500 2026-09-02 11:35 by 大不刘6
[基金申请] 面上合作单位盖章 +6 ssyjh 2026-08-27 9/450 2026-09-01 20:01 by huagongfeihu
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 17/850 2026-09-01 19:51 by 超级无敌华子
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
信息提示
请填处理意见