24小时热门版块排行榜

>论坛更新日志 (2862)
>虫友互识 (636)
>文献求助 (144)
>导师招生 (129)
>休闲灌水 (126)
>考博 (121)
>教师之家 (77)
>找工作 (70)
>博后之家 (68)
>硕博家园 (68)
>论文投稿 (66)
>论文道贺祈福 (63)
>基金申请 (58)
>考研 (52)
>公派出国 (48)
>招聘信息布告栏 (31)

返回列表

屠户

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 88.8
帖子: 28
在线: 72小时
虫号: 221604
注册: 2006-03-20
专业: 无机材料化学

[求助] 【turbomole】练习mpgrad所遇的问题，请求帮助，非常感谢

初学tm，谢谢关注，不管怎样，非常感谢！

借助define，使用mpgrad计算benzene的能量（并行在一个节点12个procs）：

完成dscf后（convergence criteria satisfied），在命令行中输入：

mp2prep -e -m memory -p discspace

出现如下错误：

ERROR! grep of disk storage requirement in /tmp/mp2prep.15822.mpout failed

之后运行mpgrad，end abnormally。以下是master文件最后几行的出错信息：
MODTRACE: no modules on stack
fatal error in trfdrv.
mpgrad ended abnormally

--------------下面是control文件，有点长------------------------------
$title
$operating system unix
$symmetry d6h
$redundant file=coord
$coord file=coord
$intdef file=coord
$user-defined bonds file=coord
$atoms
c  1-6                                                                      \
basis =c TZV                                                             \
cbas  =c cc-pVTZ
h  7-12                                                                      \
basis =h TZV                                                             \
cbas  =h cc-pVTZ
$basis file=basis
$rundimensions
dim(fock,dens)=5307
natoms=12
nshell=66
nbf(CAO)=102
nbf(AO)=102
dim(trafo[SAO<-->AO/CAO])=676
rhfshells=1
$scfmo file=mos
$closed shells
a1g    1-3                                  ( 2 )
e1g    1                                     ( 2 )
e2g    1-3                                  ( 2 )
a2u    1                                     ( 2 )
b1u    1                                     ( 2 )
b2u    1-2                                  ( 2 )
e1u    1-3                                  ( 2 )
$scfiterlimit    30
$scfconv       7
$thize    0.10000000E-04
$thime       5
$scfdamp start=0.300  step=0.050  min=0.100
$scfdump
$scfintunit
unit=30    size=0       file=twoint
$scfdiis
$scforbitalshift  automatic=.1
$drvopt
cartesian  on
basis    off
global    off
hessian on
dipole    on
nuclear polarizability
$interconversion  off
qconv=1.d-7
maxiter=25
$optimize
internal on
redundant  off
cartesian  off
global    off
basis    off logarithm
$coordinateupdate
dqmax=0.3
interpolate  on
statistics 5
$forceupdate
ahlrichs numgeo=0  mingeo=3 maxgeo=4 modus= dynamic fail=0.3
threig=0.005  reseig=0.005  thrbig=3.0  scale=1.00  damping=0.0
$forceinit on
diag=default
$energy file=energy
$grad file=gradient
$forceapprox file=forceapprox
$lock off
$denconv    0.10000000E-06
$cbas file=auxbasis
$actual step    mpgrad
$2e-ints_shell_statistics file=metastase
$orbital_max_rnorm 0.23869415437551E-05
$last SCF energy change = -230.67052
$dipole from dscf
  x    0.00000000000002 y    0.00000000000001 z    0.00000000000000 a.u.
| dipole | = 0.0000000000  debye
$scratch files
mpgrad  ddens       ./ddens_14697
mpgrad  mdens       ./mdens_14697
mpgrad  quasy       ./quasy_14697
mpgrad  GRD          ./GRD_14697
mpgrad  DIPOL       ./DIPOL_14697
mpgrad  restartg    ./restartg_14697
mpgrad  restarte    ./restarte_14697
mpgrad  restart.mpgrad ./restart.mpgrad_14697
mpgrad  recstr       ./recstr_14697
mpgrad  recstb       ./recstb_14697
mpgrad  indices       ./indices_14697
$mp2energy
$mointunit
type=intermed unit=61    size=2       file=halfint
type=1112    unit=63    size=0       file=moint#1
type=1122    unit=64    size=0       file=moint#j
type=1212    unit=65    size=1       file=moint#k
type=1212a unit=70    size=0       file=moint#a
type=gamma#1  unit=71    size=0       file=gamma#1
type=gamma#2  unit=72    size=0       file=gamma#2
type=dtdb#1 unit=76    size=0       file=dtdb#1
type=dtdb#2 unit=77    size=0       file=dtdb#2
$traloop 12
$statistics  off
$end
------------------------------以上是control文件--------------------

在tm-forum上查了下，有人说是HF计算没收敛，但是dscf显示的是收敛的。
不知道是不是自己在环境变量的设置上出问题（按照pwzhou兄的建议进行的），该贴链接：http://muchong.com/bbs/viewthread.php?tid=2058578&fpage=1

另外，还想问一下：data group $actual step is not empty
                           due to the abend of mpgrad
这个是什么意思？

对您所给予的关注和解答，不胜感激。
o了