24小时热门版块排行榜    

查看: 2891  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 不应该看fileCode +6 且听虎啸 2026-08-12 8/400 2026-08-13 12:07 by 小木虫6752397
[基金申请] FileCode能看出啥? +8 要乐观耀哥 2026-08-10 26/1300 2026-08-13 12:01 by Tide man
[基金申请] 重要来源:本周末出结果 +9 瞬息宇宙 2026-08-12 9/450 2026-08-13 10:12 by Vivilian
[基金申请] 好奇怪的filecode +5 布布和一二 2026-08-08 6/300 2026-08-12 16:11 by 云上清扬
[基金申请] filecode +13 documentary 2026-08-10 15/750 2026-08-12 15:56 by 云上清扬
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +8 Tide man 2026-08-10 9/450 2026-08-11 20:39 by beefly
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +3 majunge000 2026-08-11 4/200 2026-08-11 20:13 by lch2012
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
信息提示
请填处理意见