Firefly and PC GAMESS-related discussion club


 
Learn how to ask questions correctly  
 
 
We are NATO-free zone
 



XMCQDPT calculation, excited states in ISTSYM=2, CSFS-Problem

Johannes
johannes.heuser@pci.uni-heidelberg.de


Hello,
I have the following problem:
I try to calculate excited state energies for a Polyene (C2h) in the Bu symmetry.
My input works fine for the Ag states, but if I switch ISTSYM=1 to ISTSYM=2, I get the following error message:

********** ROUTINE MQCONF **********************************************
###   RESULTS   ###
NUMBER OF ORBITAL CONFIGURATIONS =      1107
MAXIMUM NUMBER OF SPIN FUNCTIONS =        14
NUMBER OF CSFS                   =      1764

CPU        TIME:   STEP =      0.01 ,  TOTAL =       49.4 SECONDS (    0.8 MIN)
WALL CLOCK TIME:   STEP =      0.01 ,  TOTAL =       51.2 SECONDS (    0.9 MIN)
CPU UTILIZATION:   STEP =     77.97%,  TOTAL =      96.63%
***** CPU TIME FOR MQCONF =          0.010 SEC.

********** ROUTINE MQCIG1 **********************************************

CPU        TIME:   STEP =      0.11 ,  TOTAL =       49.5 SECONDS (    0.8 MIN)
WALL CLOCK TIME:   STEP =      0.11 ,  TOTAL =       51.3 SECONDS (    0.9 MIN)
CPU UTILIZATION:   STEP =    101.32%,  TOTAL =      96.64%
***** CPU TIME FOR MQCIG1 =          0.110 SEC.

********** ROUTINE MQSYM  **********************************************
###   RESULTS   ###
# OF SYMMETRY ADAPTED CSFS =         0
*** ERROR STOP IN SUB.MQSYM ***
# OF SYMMETRY ADAPTED CSFS MUST BE GREATER THAN 0.

ADDRESS 0x0921F945 HAS INITIATED PROGRAM ABORT BECAUSE OF FATAL ERROR(S)

My Input is as follows:

$CONTRL SCFTYP=MCSCF MPLEVL=2 RUNTYP=ENERGY
UNITS=ANGS COORD=UNIQUE ICHARG=0 MULT=1 MAXIT=50 $END
$SYSTEM TIMLIM=30000 MEMORY=8000000 $END
$BASIS extfil=.t. GBASIS=cc-pVDZ $END
$DATA
ccd XMCQDPT Ag Au Bu Bg cc-pVDZ
Cnh 2
...
$END
$GUESS  GUESS=MOREAD  NORB=170 NORDER=1
IORDER(31)=31,37,39,32,33,34,35,36,38 $END
$MCSCF CISTEP=ALDET ntrack=10 $END
$DET NCORE=25 NACT=8 NELS=8 NSTATE=10 WSTATE(1)=1,1,1,-0 GROUP=C2h ISTSYM=2 $END  
$XMCQDPT KSTATE(1)=1,1,1,-0 ISTSYM=2 INORB=0 WSTATE(1)=1,1,1,-0 $END
$VEC
...
$END

Please advice on how to solve this problem. If there is anything wrong with my input that would also apply to the Ag symmetry, please let me know. The basis set was taken from the EMSL library

Also I'd like to report that ISTSYM=2 corresponds to the Bu symmetry unlike stated in the manual (according to which it should be ISTSYM=4 or ISTSYM=3 (in the example section))

Any help would be appreciated.


[ Previous ] [ Next ] [ Index ]           Wed Aug 27 '14 2:13pm
[ Reply ] [ Edit ] [ Delete ]           This message read 442 times