24小时热门版块排行榜    

查看: 862  |  回复: 7

1989shaonan

金虫 (小有名气)


[交流] PCM溶剂化中优化构型出错,L301错误

在PCM模型中优化12个水分子组成的环状团簇,溶剂半径radii = UAKS,结果出现L301错误,输入和出错界面如下,请问这是什么原因导致的错误?该如何解决?谢谢各位大侠!

输入文件
# opt freq rb3lyp/6-311+g(d,p) scrf=(solvent=water,read,pcm) nosymm
pop=npa

12w-pcm-20150423

0 1
H                  0.68424400    3.67313200    1.39072600
O                  1.60129400    3.33449100    1.34403400
H                  2.09676300    3.76387300    2.06550200
H                  4.33478000   -0.57290500    1.97189600
O                  3.66287500   -0.72317000    1.28082900
H                  4.06265700   -0.44720100    0.43300200
H                  2.41882800    3.30312900   -0.43894600
O                  2.45200900    2.81506900   -1.28457500
H                  2.65511800    3.46878900   -1.97954900
H                  3.51968600    1.24952500   -1.38173100
O                  3.68964900    0.28652400   -1.33973400
H                  4.31222700    0.07910300   -2.06059100
H                 -0.67607600   -3.66330900   -1.36918500
O                 -1.59751500   -3.33476900   -1.32826600
H                 -2.08416300   -3.77183900   -2.05103200
O                 -3.66745900    0.71679400   -1.29231800
H                 -4.33182100    0.57026300   -1.99142500
H                 -4.07684300    0.43899100   -0.45005000
H                 -2.66118600   -3.46655600    1.98757900
O                 -2.45310900   -2.81364400    1.29338800
H                 -2.41530900   -3.30173100    0.44779300
H                 -4.30816400   -0.07988900    2.05451100
O                 -3.68576700   -0.28647400    1.33328100
H                 -3.50947600   -1.24831700    1.38026300
H                 -2.22344100    3.69506600   -2.05405300
O                 -2.09399300    3.05110700   -1.33370300
H                 -2.84350300    2.42336600   -1.37879200
H                  1.66695300   -4.04087600   -1.97987200
O                  1.20684400   -3.53194800   -1.28650600
H                  1.64995300   -3.74083600   -0.44111200
H                 -1.65057000    3.75172800    0.45122400
O                 -1.20753900    3.53259500    1.29361100
H                 -1.66914300    4.03212100    1.99286200
H                  2.22055800   -3.69557800    2.05648500
O                  2.09212800   -3.05358200    1.33421900
H                  2.83957400   -2.42296900    1.37981100

radii=uaks

输出文件出错界面:

Polarizable Continuum Model (PCM)
=================================
Model                : PCM.
Atomic radii         : UAKS (United Atom Topological Model).
Polarization charges : Total charges.
Charge compensation  : None.
Solution method      : Matrix inversion.
Cavity               : GePol (RMin=0.200 OFac=0.890).
                        Default sphere list used, NSphG=   12.
                        Tesserae with average area of 0.200 Ang**2.
1st derivatives      : Analytical V*U(x)*V algorithm (CHGder, D1EAlg=0).
                        Cavity 1st derivative terms included.
Solvent              : Water, Eps     =  78.390000
                               Eps(inf)=   1.776000
                               RSolv   =   1.385000 Ang.
------------------------------------------------------------------------------
AdVTs1: ISph=  887 is engulfed by JSph=  888 but Ae(  887) is not yet zero!
Error termination via Lnk1e in /home/user/opt/g03/l301.exe at Sat Feb  9 01:31:50 2013.
Job cpu time:  1 days  0 hours 10 minutes 20.2 seconds.
File lengths (MBytes):  RWF=   1490 Int=      0 D2E=      0 Chk=     16 Scr=      1
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kong890109

金虫 (小有名气)



小木虫: 金币+0.5, 给个红包,谢谢回帖
Problem is related to building of the cavity in solvent calculations. One possible solution is to change the cavity model (default in g03 is UAO, can be changed by adding RADII keyword in section below coordinates in the .inp file, e.g. RADII=UFF, see SCRF keyword).
链接地址:http://blog.sina.com.cn/s/blog_8be57be70100xabz.html
8楼2015-04-28 09:10:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
syhorchid2楼
2015-04-28 00:17   回复  
1989shaonan(金币+1): 谢谢参与
dmbb3楼
2015-04-28 08:28   回复  
1989shaonan(金币+1): 谢谢参与
2015-04-28 08:31   回复  
1989shaonan(金币+1): 谢谢参与
mtw9225楼
2015-04-28 08:34   回复  
1989shaonan(金币+1): 谢谢参与
2015-04-28 08:38   回复  
1989shaonan(金币+1): 谢谢参与
2015-04-28 08:53   回复  
相关版块跳转 我要订阅楼主 shaonan 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见