24小时热门版块排行榜    

查看: 2923  |  回复: 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文章,我:8O.5.5.1O.54,科目全,可十急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 05:45 by gy1nBQXYQJqL
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +4 gouxfjh 2026-08-28 4/200 2026-08-29 04:46 by alongwaytogo
[基金申请] 国自然评审意见 +11 wangmingqi 2026-08-28 15/750 2026-08-29 00:17 by fangyl2005
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
请填处理意见