Panwang Zhou
pwzhou@gmail.com
Recently, I performed a MCSCF(16,14) calculation then following a XMC-QDPT2 calculation on the constructed LIIC curve among the optimized cis-HBDI, CI and trans-HBDI structures, the obtained MCSCF curve is smooth and continuous whereas the obtained XMC-QDPT2 curve is discontinuous at the CI point (see attachment). I'm not clear why the XMC-QDPT2 curve is discontinuous and how to resolve this problem?
The following is the part of input file:
$CONTRL SCFTYP=MCSCF RUNTYP=ENERGY EXETYP=RUN MAXIT=50 ICHARG=-1
MULT=1 FSTINT=.T. GENCON=.T. INTTYP=HONDO NOSYM=1
ICUT=11 ITOL=30 WIDE=.T. MPLEVL=2 $END
$SYSTEM MWORDS=480 TIMLIM=60000.0 KDIAG=0 NOJAC=100 $END
$SYSTEM MKLNP=1 NP=12 $END
$SMP SMPPAR=.T. HTTNP=1 $END
$SCF DIRSCF=.T. FDIFF=.F. NCONV=8 $END
$P2P P2P=.T. DLB=.T. $END
$TRANS MPTRAN=2 DIRTRF=.T. AOINTS=DIST ALTPAR=.T. MODE=112 $END
$MCSCF CISTEP=ALDET FULLNR=.F. SOSCF=.T. MAXIT=100 $END
$MCSCF IFORB=.T. $END
$DET NCORE=49 NACT=14 NELS=16 NSTATE=6 WSTATE(1)=1,1 DISTCI=12 $END
$XMCQDPT NSTATE=2 EDSHFT=0.02 THRGEN=1D-12 MXBASE=90 $END
$XMCQDPT HALLOC=.T. $END
$XMCQDPT IFORB(1)=-1,1,1 WSTATE(1)=1,1,-0 AVECOE(1)=1,1,-0 $END
$BASIS GBASIS=N31 NGAUSS=6 NDFUNC=1 NPFUNC=1 DIFFSP=.TRUE. $END
$GUESS GUESS=MOREAD NORB=359 $END
$MCQFIT $END
I have tried to increas the NSTATE in XMCQDPT group, however, the required memory is larger than 480 MWORDS, thus the calculation failed.
Best Regards!
Panwang Zhou
This message contains the 68 kb attachment [ LIIC_curves.zip ] MCSCF and XMC-QDPT2 Curves |