应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Molpro/Molcas » 【求助】molpro 计算问题，请高手指教！！！！
131/2返回列表12下一页
查看: 2371  |  回复: 12
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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这个关键词怎么用？请大家给些宝贵意见，谢谢，很急切！！！！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2010-07-05 18:31:52
已阅   回复此楼   关注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]

另外，默认的活性空间很大，你可能需要自己定义活性空间
赞一下(2人) 回复此楼
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可以试试态平均。
赞一下(1人) 回复此楼
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
赞一下(1人) 回复此楼
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挑出来。
赞一下(1人) 回复此楼
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 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭