24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1109  |  回复: 13
当前主题已经存档。

liuxy861120

新虫 (初入文坛)

[交流] 【求助】请问高斯计算中出现这样的错误怎么处理

我在计算的过程中出现了这样的问题。虽然最大力和最大位移等都收敛了,但是在计算频率的时候出现了这样的东西,不知道是什么意思,Defaulting to unpruned grid for atomic number  47
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

easygeton

铜虫 (小有名气)


liuxy861120(金币+5): 2010-03-08 15:58
yjcmwgk(金币+1): 2010-03-08 16:00
这不是错误。意思是对于47号原子使用“未修剪”的积分密度,不用管它。
2楼2010-03-08 14:37:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuxy861120

新虫 (初入文坛)

引用回帖:
Originally posted by easygeton at 2010-03-08 14:37:19:
这不是错误。意思是对于47号原子使用“未修剪”的积分密度,不用管它。

那为什么后面的频率计算停止了呢?
3楼2010-03-08 16:00:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiewei

专家顾问 (著名写手)

小才

引用回帖:
Originally posted by liuxy861120 at 2010-03-08 16:00:54:

那为什么后面的频率计算停止了呢?

你把最后那段文字贴出来 我帮你看看
我想做个科学家,可我又喜欢写诗;我想做一个诗人,可我却读的理工科..
4楼2010-03-08 18:38:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuxy861120

新虫 (初入文坛)

**********************************************************************

            Population analysis using the SCF density.

**********************************************************************

Orbital symmetries:
Alpha Orbitals:
       Occupied  (A1) (A1) (B2) (B2) (A1) (B2) (A1) (B1) (A1) (A1)
                 (B2) (B2) (B1) (A2) (A1) (B2) (A1) (B2) (A2) (B1)
                 (A1) (B2) (A1) (A1) (A2) (A1) (B1) (B2) (B1) (A2)
                 (B2) (B1) (A1) (B2) (A2) (B1) (A1) (B1) (B2) (A1)
                 (A2) (A2) (B2) (B2) (A1) (A1) (B2) (A1)
       Virtual   (A1) (B2) (B1) (A1) (B2) (A2) (A1) (B1) (B2) (A1)
                 (A1) (B1) (B2) (A2) (B2) (A1) (B2) (B1) (A1) (A1)
                 (B2) (A2) (B1) (B2) (B1) (A2) (A1) (A1) (A1) (B2)
                 (B2) (B2) (A1) (B1) (A2) (B1) (A2) (B2) (A1) (A1)
                 (B1) (A1) (B2) (B1) (A1) (A2) (B2) (B1) (A2) (A1)
                 (B2) (A1) (B2) (A2) (A1) (B2) (A1) (B2) (B2) (A1)
                 (A1) (B2)
Beta  Orbitals:
       Occupied  (A1) (A1) (B2) (B2) (A1) (B2) (A1) (B1) (A1) (B2)
                 (A1) (B2) (B1) (A2) (A1) (B2) (A1) (B2) (A2) (B1)
                 (A1) (B2) (A1) (A1) (A2) (A1) (B1) (B2) (B1) (A2)
                 (B2) (B1) (A1) (B2) (A2) (B1) (A1) (B1) (B2) (A1)
                 (A2) (A2) (B2) (B2) (A1) (A1) (B2)
       Virtual   (A1) (A1) (B2) (B1) (A1) (B2) (A2) (A1) (B1) (B2)
                 (A1) (A1) (B1) (B2) (A2) (B2) (A1) (B2) (B1) (A1)
                 (A1) (B2) (A2) (B1) (B2) (B1) (A2) (A1) (A1) (A1)
                 (B2) (B2) (B2) (A1) (A2) (B1) (B1) (A2) (B2) (A1)
                 (A1) (B1) (A1) (B2) (B1) (A1) (A2) (B2) (B1) (A2)
                 (A1) (B2) (A1) (B2) (A2) (A1) (B2) (A1) (B2) (B2)
                 (A1) (A1) (B2)
The electronic state is 2-A1.
Alpha  occ. eigenvalues --   -3.61404  -3.60889  -3.60877  -3.60350  -3.60350
Alpha  occ. eigenvalues --   -2.25475  -2.25360  -2.25012  -2.24966  -2.24826
Alpha  occ. eigenvalues --   -2.24821  -2.24743  -2.24549  -2.24530  -2.24383
Alpha  occ. eigenvalues --   -2.24342  -2.24178  -2.24171  -2.23979  -2.23978
Alpha  occ. eigenvalues --   -0.37989  -0.37468  -0.37000  -0.36079  -0.35750
Alpha  occ. eigenvalues --   -0.35606  -0.35439  -0.35033  -0.34889  -0.33816
Alpha  occ. eigenvalues --   -0.33488  -0.33196  -0.33125  -0.32883  -0.32675
Alpha  occ. eigenvalues --   -0.32577  -0.32348  -0.32129  -0.31977  -0.31954
Alpha  occ. eigenvalues --   -0.31744  -0.31573  -0.31398  -0.30664  -0.30523
Alpha  occ. eigenvalues --   -0.26626  -0.20329  -0.16792
Alpha virt. eigenvalues --   -0.08837  -0.07999  -0.05235  -0.03497  -0.01833
Alpha virt. eigenvalues --   -0.01044  -0.00326   0.00809   0.01682   0.03311
Alpha virt. eigenvalues --    0.03814   0.04146   0.04472   0.04554   0.05191
Alpha virt. eigenvalues --    0.06469   0.08161   0.08195   0.08507   0.09387
Alpha virt. eigenvalues --    0.09803   0.10142   0.11352   0.12809   0.14483
Alpha virt. eigenvalues --    0.14816   0.16325   0.19275   0.21415   0.21910
Alpha virt. eigenvalues --    0.25187   0.31936   0.44766   0.48530   0.48660
Alpha virt. eigenvalues --    0.49386   0.54855   0.55376   0.55700   0.58182
Alpha virt. eigenvalues --    0.62655   0.62677   0.62909   0.69712   0.71748
Alpha virt. eigenvalues --    0.72545   0.77657   0.77948   0.78115   0.78954
Alpha virt. eigenvalues --    0.80728   0.80998   0.85626   0.86348   0.87366
Alpha virt. eigenvalues --    0.98587   1.00234   1.01277   1.35123   1.44084
Alpha virt. eigenvalues --    1.70543   1.84619
  Beta  occ. eigenvalues --   -3.61330  -3.60831  -3.60819  -3.60313  -3.60313
  Beta  occ. eigenvalues --   -2.25370  -2.25293  -2.24987  -2.24883  -2.24793
  Beta  occ. eigenvalues --   -2.24793  -2.24665  -2.24527  -2.24511  -2.24348
  Beta  occ. eigenvalues --   -2.24306  -2.24140  -2.24130  -2.23963  -2.23962
  Beta  occ. eigenvalues --   -0.37833  -0.37342  -0.36833  -0.35842  -0.35659
  Beta  occ. eigenvalues --   -0.35395  -0.35321  -0.34935  -0.34737  -0.33723
  Beta  occ. eigenvalues --   -0.33389  -0.33095  -0.33004  -0.32793  -0.32574
  Beta  occ. eigenvalues --   -0.32482  -0.32220  -0.31993  -0.31881  -0.31826
  Beta  occ. eigenvalues --   -0.31659  -0.31505  -0.31195  -0.30564  -0.30303
  Beta  occ. eigenvalues --   -0.25814  -0.19905
  Beta virt. eigenvalues --   -0.12459  -0.08487  -0.07119  -0.04947  -0.02493
  Beta virt. eigenvalues --   -0.01622  -0.00935   0.00179   0.01247   0.01982
  Beta virt. eigenvalues --    0.03505   0.03853   0.04205   0.04578   0.04648
  Beta virt. eigenvalues --    0.05350   0.06702   0.08218   0.08258   0.08752
  Beta virt. eigenvalues --    0.09627   0.09909   0.10224   0.11787   0.12997
  Beta virt. eigenvalues --    0.14668   0.15140   0.16617   0.19830   0.21951
  Beta virt. eigenvalues --    0.22138   0.25684   0.32178   0.45022   0.48928
  Beta virt. eigenvalues --    0.48961   0.49827   0.55178   0.55698   0.56042
  Beta virt. eigenvalues --    0.58607   0.62934   0.63054   0.63225   0.70117
  Beta virt. eigenvalues --    0.72247   0.72840   0.77966   0.78395   0.78530
  Beta virt. eigenvalues --    0.79310   0.80940   0.81177   0.85852   0.86634
  Beta virt. eigenvalues --    0.87572   0.99029   1.00642   1.01638   1.35417
  Beta virt. eigenvalues --    1.44423   1.70736   1.84763
          Condensed to atoms (all electrons):
              1          2          3          4          5
     1  Ag  18.665319   0.155592   0.030642   0.006616   0.142897
     2  Ag   0.155592  18.665319   0.030642   0.142897   0.006616
     3  Ag   0.030642   0.030642  18.493252   0.150805   0.150805
     4  Ag   0.006616   0.142897   0.150805  18.770045   0.000497
     5  Ag   0.142897   0.006616   0.150805   0.000497  18.770045
Mulliken atomic charges:
              1
     1  Ag  -0.001067
     2  Ag  -0.001067
     3  Ag   0.143854
     4  Ag  -0.070860
     5  Ag  -0.070860
Sum of Mulliken charges=   0.00000
Atomic charges with hydrogens summed into heavy atoms:
              1
     1  Ag  -0.001067
     2  Ag  -0.001067
     3  Ag   0.143854
     4  Ag  -0.070860
     5  Ag  -0.070860
Sum of Mulliken charges=   0.00000
          Atomic-Atomic Spin Densities.
              1          2          3          4          5
     1  Ag   0.229017   0.114717  -0.037058  -0.009063  -0.032681
     2  Ag   0.114717   0.229017  -0.037058  -0.032681  -0.009063
     3  Ag  -0.037058  -0.037058   0.147191   0.041857   0.041857
     4  Ag  -0.009063  -0.032681   0.041857   0.153722   0.002841
     5  Ag  -0.032681  -0.009063   0.041857   0.002841   0.153722
Mulliken atomic spin densities:
              1
     1  Ag   0.264931
     2  Ag   0.264931
     3  Ag   0.156789
     4  Ag   0.156674
     5  Ag   0.156674
Sum of Mulliken spin densities=   1.00000
APT atomic charges:
              1
     1  Ag   0.015980
     2  Ag   0.015980
     3  Ag   0.124204
     4  Ag  -0.078082
     5  Ag  -0.078082
Sum of APT charges=   0.00000
APT Atomic charges with hydrogens summed into heavy atoms:
              1
     1  Ag   0.015980
     2  Ag   0.015980
     3  Ag   0.124204
     4  Ag  -0.078082
     5  Ag  -0.078082
Sum of APT charges=   0.00000
Electronic spatial extent (au):  =  2009.5392
Charge=     0.0000 electrons
Dipole moment (field-independent basis, Debye):
    X=     0.0000    Y=     0.0000    Z=    -0.0440  Tot=     0.0440
Quadrupole moment (field-independent basis, Debye-Ang):
   XX=  -101.0660   YY=   -92.0080   ZZ=   -97.8312
   XY=     0.0000   XZ=     0.0000   YZ=     0.0000
Traceless Quadrupole moment (field-independent basis, Debye-Ang):
   XX=    -4.0976   YY=     4.9604   ZZ=    -0.8628
   XY=     0.0000   XZ=     0.0000   YZ=     0.0000
Octapole moment (field-independent basis, Debye-Ang**2):
  XXX=     0.0000  YYY=     0.0000  ZZZ=    -2.2932  XYY=     0.0000
  XXY=     0.0000  XXZ=    -0.2561  XZZ=     0.0000  YZZ=     0.0000
  YYZ=     8.3734  XYZ=     0.0000
Hexadecapole moment (field-independent basis, Debye-Ang**3):
XXXX=  -131.6555 YYYY= -2289.7987 ZZZZ= -1001.9358 XXXY=     0.0000
XXXZ=     0.0000 YYYX=     0.0000 YYYZ=     0.0000 ZZZX=     0.0000
ZZZY=     0.0000 XXYY=  -428.8712 XXZZ=  -188.2528 YYZZ=  -539.8184
XXYZ=     0.0000 YYXZ=     0.0000 ZZXY=     0.0000
N-N= 5.911020269201D+02 E-N=-2.748726949070D+03  KE= 2.700353844003D+02
Symmetry A1   KE= 9.218909490852D+01
Symmetry A2   KE= 4.499323091149D+01
Symmetry B1   KE= 4.991836504031D+01
Symmetry B2   KE= 8.293469353998D+01
  Exact polarizability: 126.914   0.000 372.507   0.000   0.000 243.499
Approx polarizability: 213.994   0.000 751.059   0.000   0.000 469.779
                          Isotropic Fermi Contact Couplings
        Atom                 a.u.       MegaHertz       Gauss      10(-4) cm-1
     1  Ag(107             0.00002      -0.00411      -0.00147      -0.00137
     2  Ag(107             0.00002      -0.00411      -0.00147      -0.00137
     3  Ag(107             0.00002      -0.00279      -0.00100      -0.00093
     4  Ag(107             0.00002      -0.00278      -0.00099      -0.00093
     5  Ag(107             0.00002      -0.00278      -0.00099      -0.00093
--------------------------------------------------------
       Center         ----  Spin Dipole Couplings  ----
                      3XX-RR        3YY-RR        3ZZ-RR
--------------------------------------------------------
     1   Atom       -0.008925      0.029080     -0.020155
     2   Atom       -0.008925      0.029080     -0.020155
     3   Atom       -0.013666      0.027090     -0.013424
     4   Atom       -0.004099      0.010909     -0.006810
     5   Atom       -0.004099      0.010909     -0.006810
--------------------------------------------------------
                        XY            XZ            YZ
--------------------------------------------------------
     1   Atom        0.000000      0.000000      0.007726
     2   Atom        0.000000      0.000000     -0.007726
     3   Atom        0.000000      0.000000      0.000000
     4   Atom        0.000000      0.000000      0.013438
     5   Atom        0.000000      0.000000     -0.013438
--------------------------------------------------------


---------------------------------------------------------------------------------
              Anisotropic Spin Dipole Couplings in Principal Axis System
---------------------------------------------------------------------------------

       Atom             a.u.   MegaHertz   Gauss  10(-4) cm-1        Axes

              Baa    -0.0213     0.463     0.165     0.154  0.0000 -0.1515  0.9885
     1 Ag(107 Bbb    -0.0089     0.194     0.069     0.065  1.0000  0.0000  0.0000
              Bcc     0.0303    -0.657    -0.234    -0.219  0.0000  0.9885  0.1515

              Baa    -0.0213     0.463     0.165     0.154  0.0000  0.1515  0.9885
     2 Ag(107 Bbb    -0.0089     0.194     0.069     0.065  1.0000  0.0000  0.0000
              Bcc     0.0303    -0.657    -0.234    -0.219  0.0000  0.9885 -0.1515

              Baa    -0.0137     0.297     0.106     0.099  1.0000  0.0000  0.0000
     3 Ag(107 Bbb    -0.0134     0.291     0.104     0.097  0.0000  0.0000  1.0000
              Bcc     0.0271    -0.588    -0.210    -0.196  0.0000  1.0000  0.0000

              Baa    -0.0140     0.305     0.109     0.102  0.0000 -0.4741  0.8805
     4 Ag(107 Bbb    -0.0041     0.089     0.032     0.030  1.0000  0.0000  0.0000
              Bcc     0.0181    -0.394    -0.140    -0.131  0.0000  0.8805  0.4741

              Baa    -0.0140     0.305     0.109     0.102  0.0000  0.4741  0.8805
     5 Ag(107 Bbb    -0.0041     0.089     0.032     0.030  1.0000  0.0000  0.0000
              Bcc     0.0181    -0.394    -0.140    -0.131  0.0000  0.8805 -0.4741


---------------------------------------------------------------------------------

   4 Symmetry operations used in ECPInt.
ECPInt:  NShTT=     820 NPrTT=    4175 LenC2=     821 LenP2D=    4175.
LDataN:  DoStor=F MaxTD1= 7 Len=  274
LDataN:  DoStor=T MaxTD1= 7 Len=  274
Defaulting to unpruned grid for atomic number  47.

这就是后面的输出内容,你帮我看看啊,谢谢

你把最后那段文字贴出来 我帮你看看 [/quote]
5楼2010-03-09 14:30:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiewei

专家顾问 (著名写手)

小才

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
yjcmwgk(金币+1): 2010-03-12 13:21
引用回帖:
Originally posted by liuxy861120 at 2010-03-09 14:30:17:
**********************************************************************

            Population analysis using the SCF density.

**************************************************************** ...

你建立输入文件的时候,是不是增加对称性了? 你对称性改为 loose 或者 nosymm试试,  你的回复我这里咋收不到消息呢,我还是今天无事的时候点击我的回复话题才开见的。。。
我想做个科学家,可我又喜欢写诗;我想做一个诗人,可我却读的理工科..
6楼2010-03-09 14:41:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuxy861120

新虫 (初入文坛)

引用回帖:
Originally posted by jiewei at 2010-03-09 14:41:47:
哦 谢谢你啊 我试试看 我也没有增加对称性啊 这是别人算的东西,我试着用别人的集合构型算的 算了之后就出现了问题 我也从没遇见过这种情况。至于你说的没有消息的事情,我也不是很清楚,估计我回复你的过程中出现了,问题吧
你建立输入文件的时候,是不是增加对称性了? 你对称性改为 loose 或者 nosymm试试,  你的回复我这里咋收不到消息呢,我还是今天无事的时候点击我的回复话题才开见的。。。

Sample Text
7楼2010-03-09 16:53:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangmt

至尊木虫 (著名写手)

我叫MT

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
yjcmwgk(金币+3): 2010-03-09 18:09
你没加对称性,但是gaussian计算之前会先判断对称性,虽然你没有给出坐标信息,但从数据处理上看,显然这是一个高度对称的体系。
首先拿收敛的结构出来看看,是不是明显的不合理,
如果自认为合理,那么加大基组大概就可以解决问题,
如果还不行,尝试加nosymm或者grid=ultrafine
说白了,很大可能就是体系跑散了,或者距离太近了,总之结构不合理造成
没法进行轨道重叠的猜测了——如果还就是合理的,那就是你要另外加极化或弥散轨道了。

[ Last edited by zhangmt on 2010-3-9 at 18:12 ]
一群自以为正义凛然的年轻人将一切不能以科学解释的事情定性为封建迷信并大刀阔斧地进行消灭,其实这是修养不足学识浅薄的一种体现,也是可恶的偏执和愚蠢的自以
8楼2010-03-09 18:07:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiewei

专家顾问 (著名写手)

小才

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
yjcmwgk(金币+3): 2010-03-10 10:42
引用回帖:
Originally posted by zhangmt at 2010-03-09 18:07:11:
你没加对称性,但是gaussian计算之前会先判断对称性,虽然你没有给出坐标信息,但从数据处理上看,显然这是一个高度对称的体系。
首先拿收敛的结构出来看看,是不是明显的不合理,
如果自认为合理,那么加大基组 ...

其实我个人经验来看:可以用CHEM画好分子结构后,放入CHEM3D初步优化,这个时候可以根据经验来判断分子大致在什么结构比较合理,然后放入高斯 不加任何对称性情况下优化,我一般是这么做的,结构一般也不赖,所以到现在为止基本没有出现过楼主这样的情况。
经验之谈,呵呵 仅供参考,当然这个需要你比较多实际经验,比如看到一个结构,你大体就能知道它的优化结构在怎么样的
我想做个科学家,可我又喜欢写诗;我想做一个诗人,可我却读的理工科..
9楼2010-03-09 19:09:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuxy861120

新虫 (初入文坛)

引用回帖:
Originally posted by jiewei at 2010-03-09 19:09:43:
谢谢你的建议 我试试看 呵呵
其实我个人经验来看:可以用CHEM画好分子结构后,放入CHEM3D初步优化,这个时候可以根据经验来判断分子大致在什么结构比较合理,然后放入高斯 不加任何对称性情况下优化,我一般是这么做的,结构一般也不赖,所 ...

Sample Text
10楼2010-03-09 21:57:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 liuxy861120 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见