Gamesscomsol求解不收敛敛,求帮助

小木虫 --- 600万学术达人喜爱的学术科研平台
热门搜索:
&&Gamess不收敛,求帮助
Gamess不收敛,求帮助
输出文件:
...... END OF RHF CALCULATION ......
CPU& &&&0: STEP CPU TIME=& &&&1.13 TOTAL CPU TIME=& && &&&1.7 (& & 0.0 MIN)
TOTAL WALL CLOCK TIME=& && &&&3.2 SECONDS, CPU UTILIZATION IS&&52.78%
SCF DID NOT CONVERGE...NO MPLEVL=2 CALCULATION
EXECUTION OF GAMESS TERMINATED -ABNORMALLY- AT Tue Jul&&7 16:47:56 2015
& && && && && &6784206&&WORDS OF DYNAMIC MEMORY USED
CPU& &&&0: STEP CPU TIME=& &&&0.00 TOTAL CPU TIME=& && &&&1.7 (& & 0.0 MIN)
TOTAL WALL CLOCK TIME=& && &&&3.2 SECONDS, CPU UTILIZATION IS&&52.62%
DDI Process 0: error code 911
ddikick.x: application process 0 quit unexpectedly.
ddikick.x: Fatal error detected.
The error is most likely to be in the application, so check for
input errors, disk space, memory needs, application bugs, etc.
ddikick.x will now clean up all processes, and exit...
DDI Process 10: terminated upon request. DDI Process 23: terminated upon request.
ddikick.x: Sending kill signal to DDI processes.
ddikick.x: Execution terminated due to error(s).
unset echo
----- accounting info -----
Files used on the master node node15 were:
-rw-rw-r-- 1 menglp menglp& &&&1533 Jul&&7 16:47 /tmp/gamesstmp_11889/1-b-1.F05
-rw-rw-r-- 1 menglp menglp
Jul&&7 16:47 /tmp/gamesstmp_11889/1-b-1.F08
-rw-rw-r-- 1 menglp menglp
Jul&&7 16:47 /tmp/gamesstmp_11889/1-b-1.F08.001
-rw-rw-r-- 1 menglp menglp
Jul&&7 16:47 /tmp/gamesstmp_11889/1-b-1.F08.002
-rw-rw-r-- 1 menglp menglp
Jul&&7 16:47 /tmp/gamesstmp_11889/1-b-1.F08.003
-rw-rw-r-- 1 menglp menglp
Jul&&7 16:47 /tmp/gamesstmp_11889/1-b-1.F08.004
-rw-rw-r-- 1 menglp menglp
Jul&&7 16:47 /tmp/gamesstmp_11889/1-b-1.F08.005
-rw-rw-r-- 1 menglp menglp
Jul&&7 16:47 /tmp/gamesstmp_11889/1-b-1.F08.006
-rw-rw-r-- 1 menglp menglp
Jul&&7 16:47 /tmp/gamesstmp_11889/1-b-1.F08.007
-rw-rw-r-- 1 menglp menglp
Jul&&7 16:47 /tmp/gamesstmp_11889/1-b-1.F08.008
-rw-rw-r-- 1 menglp menglp
Jul&&7 16:47 /tmp/gamesstmp_11889/1-b-1.F08.009
-rw-rw-r-- 1 menglp menglp
Jul&&7 16:47 /tmp/gamesstmp_11889/1-b-1.F08.010
-rw-rw-r-- 1 menglp menglp
Jul&&7 16:47 /tmp/gamesstmp_11889/1-b-1.F08.011
-rw-rw-r-- 1 menglp menglp
Jul&&7 16:47 /tmp/gamesstmp_11889/1-b-1.F08.012
-rw-rw-r-- 1 menglp menglp
Jul&&7 16:47 /tmp/gamesstmp_11889/1-b-1.F08.013
-rw-rw-r-- 1 menglp menglp
Jul&&7 16:47 /tmp/gamesstmp_11889/1-b-1.F08.014
-rw-rw-r-- 1 menglp menglp
Jul&&7 16:47 /tmp/gamesstmp_11889/1-b-1.F08.015
-rw-rw-r-- 1 menglp menglp
Jul&&7 16:47 /tmp/gamesstmp_11889/1-b-1.F08.016
-rw-rw-r-- 1 menglp menglp
Jul&&7 16:47 /tmp/gamesstmp_11889/1-b-1.F08.017
-rw-rw-r-- 1 menglp menglp
Jul&&7 16:47 /tmp/gamesstmp_11889/1-b-1.F08.018
-rw-rw-r-- 1 menglp menglp
Jul&&7 16:47 /tmp/gamesstmp_11889/1-b-1.F08.019
-rw-rw-r-- 1 menglp menglp
Jul&&7 16:47 /tmp/gamesstmp_11889/1-b-1.F08.020
-rw-rw-r-- 1 menglp menglp
Jul&&7 16:47 /tmp/gamesstmp_11889/1-b-1.F08.021
-rw-rw-r-- 1 menglp menglp
Jul&&7 16:47 /tmp/gamesstmp_11889/1-b-1.F08.022
-rw-rw-r-- 1 menglp menglp
Jul&&7 16:47 /tmp/gamesstmp_11889/1-b-1.F08.023
-rw-rw-r-- 1 menglp menglp
Jul&&7 16:47 /tmp/gamesstmp_11889/1-b-1.F10
ls: No match.
ls: No match.
ls: No match.
Tue Jul&&7 16:47:58 CST 2015
0.137u 0.296s 0:06.74 6.2%& & & & 0+0k 8+240io 0pf+0w
输入文件:
$CONTRL SCFTYP=RHF RUNTYP=eda ICHARG=0 ispher=1 mplevl=2
&&MULT=1 $end
!$guess guess=hcore $end
$BASIS GBASIS=accd $END
$LMOEDA MATOM(1)=7 12 MCHARG(1)=0 0 MMULT(1)=1 1 $END
!$SCF DIIS=.t. SOSCF=.f. DIRSCF=.T. FDIFF=.F. NPUNCH=1 $END
$scf maxdii=100 $end
$SYSTEM TIMLIM= MWORDS=2000 memddi=2000 $END
optimized by b3lyp/ccd
P& & 15.0& && && && & 0.& & 0.& & 1.
O& & 8.0& && && && & -1.& &-0.& & 2.
O& & 8.0& && && && &&&0.& & 1.& & 2.
O& & 8.0& && && && &&&1.& &-1.& & 2.
H& & 1.0& && && && &&&1.& & 0.& & 2.
H& & 1.0& && && && &&&1.& & 1.& & 2.
H& & 1.0& && && && &&&0.& &-1.& & 2.
C& & 6.0& && && && &-0.& & 1.& &-2.&&
C& & 6.0& && && && & -1.& & 0.& &-2.
C& & 6.0& && && && & -1.& &-0.& &-2.
C& & 6.0& && && && &&&0.& &-1.& &-2.
C& & 6.0& && && && &&&1.& &-0.& &-2.
C& & 6.0& && && && &&&1.& & 0.& &-2.
H& & 1.0& && && && &&&0.& & 2.& &-2.
H& & 1.0& && && && & -2.& & 1.& &-2.
H& & 1.0& && && && & -2.& &-1.& &-2.
H& & 1.0& && && && &&&0.& &-2.& &-2.
H& & 1.0& && && && &&&2.& &-1.& &-2.
H& & 1.0& && && && &&&2.& & 1.& &-2.& &
初学GAMESS,遇到不收敛不知道什么原因……求各位大神帮助!
学术必备与600万学术达人在线互动!
扫描下载送金币GAMESS自洽场(SCF)出问题或者不收敛时的解决办法
最近做课题的时候需要需要用GAMESS取积分,但是在指定的构型下GAMESS的HF-SCF收敛的能量和GAUSSIAN或者molpro收敛的能量
不一致。记得以前听GY说过GAMESS的SCF经常出现这种问题,此时可以以GAUSSIAN的初始轨道作为初值猜测,使GAMESS的SCF也可以收
敛到更加合适的位置。
比如,以我使用的是2为氢分子体系,如下图所示(晕 图片传不上)
&H ... H ... H ... H ... H ... H
&H ... H ... H ... H ... H ... H
&H ... H ... H ... H ... H ... H
&H ... H ... H ... H ... H ... H
对应的GAMESS输入文件如下& ( HF-SCF能量为 -10.)
&$CONTRL SCFTYP=RHF RUNTYP=ENERGY ISPHER=1
MAXIT=100 $END
&$SYSTEM MWORDS=1480 $END
&$BASIS& GBASIS=n311 NGAUSS=6
而此体系GAUSSIAN或者molpro能量为& -10.
为了使GAMESS的结果与GAUSSIAN或者molpro一致,我们要做的是用GAMESS来读入GAUSSIAN的分子轨道(MOs)作为初值。GAUSSIAN的MOs可以由以下输入文件得到:
(注意这里为了与GAMESS构型保持一致,采用了原子单位以及NOSYM来固定坐标)
%chk=H4X6.chk
#HF/6-311G UNIT=AU pop=full nosym
-2.&&&&&&&
-6.&&&&&&&
-10.&&&&&&&
-10.&&&&&&
10.&&&&&&&
-10.&&&&&&&
-10.&&&&&&
10.&&&&&&&
-6.&&&&&&&
-2.&&&&&&&
计算完后,将*.chk转化为*.fchk,同时读出分子轨道并写成GAMESS的输入形式(见附),并在GAMESS输入里面添加MOREAD的部分,如下
&$CONTRL SCFTYP=RHF RUNTYP=ENERGY ISPHER=1
MAXIT=100 $END
&$SYSTEM MWORDS=1480 $END
&$BASIS& GBASIS=n311 NGAUSS=6
&$GUESS& GUESS=MOREAD norb=72
0. 0. 0. 0.
0. 0. 0. 0.
0. 0. 0. 0.
0. 0. 0. 0.
之后,GAMESS也可以得出跟GAUSSIAN或者Molpro一样的SCF能量,以及和GAUSSIAN一样的MOs。可以放心的提取分子轨道下的积分。
acknowledge:
GY&& LY(subroutine)
=========================================================================
附,转换fchk到vec(GAMESS)的程序(subroutine)
Program Read_CS
& implicit none
!parameter
integer :: natom=72
integer :: file_id
character*50 :: file_in='H4X6.fchk'
character*50 :: file_out='H4X6.vec'
character*8 :: flag='Alpha MO'
integer :: error
character*8 :: buffer
double precision,allocatable :: NBO(:)
character*100 :: string
integer :: i,j,k
allocate (NBO(natom*natom),stat=error)
if(error /= 0) then
open(unit=101,file=file_in,status='OLD',iostat=error)
if(error/=0) then
do while(.TRUE.)
read(101,"(A8)") buffer
if(buffer==flag) then
&&&&&&&&&&&
do i=1,natom**2
&&&&&&&&&&&&&
read(101,"(E16.8)",advance='NO') NBO(i)
&&&&&&&&&&&&&
if(mod(i,5)==0) then
&&&&&&&&&&&&&&&
read(101,*)
&&&&&&&&&&&&&
&&&&&&&&&&&
&&&&&&&&&&&
close(101)
open(unit=102,file=file_out,status="REPLACE",iostat=error)
if(error/=0) then
write(102,"(A5)",advance='NO') ' $VEC'
do i=1,natom
do j=1,natom
&&&&&&&&&&&
if(mod(j,5)==1) then
&&&&&&&&&&&&&
write(102,*)
&&&&&&&&&&&&&
&&&&&&&&&&&&&
do while (k & 99)
&&&&&&&&&&&&&&&
&&&&&&&&&&&&&
&&&&&&&&&&&&&
write(102,"(I2,I3)",advance='NO'),k,ceiling(j/.5d1)
&&&&&&&&&&&
&&&&&&&&&&&
write(102,"(E15.8)",advance='NO') NBO((i-1)*natom+j)
write(102,*)
write(102,"(A5)") ' $END'
close(102)
open(unit=105,file="VEC.txt")
do i=1,natom
&&&&&&&&&&&
do j=1,natom
&&&&&&&&&&&&&
write(105,*)i,j,NBO((i-1)*natom+j)
&&&&&&&&&&&
close(105)
deallocate (NBO)
已投稿到:
以上网友发言只代表其个人观点,不代表新浪网的观点或立场。}

我要回帖

更多关于 comsol求解不收敛 的文章

更多推荐

版权声明:文章内容来源于网络,版权归原作者所有,如有侵权请点击这里与我们联系,我们将及时删除。

点击添加站长微信