1 2####################################################################### 3# TEST DESCRIPTION 4####################################################################### 5cat > cc3_noddy_hyp.info <<%EOF% 6 cc3_noddy_hyp 7 ------------- 8 Molecule: LIH 9 Theory Level: CC3 / 6-31G 10 Test Purpose: static and freq.-dep. first hyperpolarizability 11 as computed with the noddy test code 12%EOF% 13 14####################################################################### 15# INTEGRAL INPUT 16####################################################################### 17cat > cc3_noddy_hyp.mol <<%EOF% 18BASIS 196-31G 20blubber 21blubber 22 2 0 X Y 1.00D-18 23 3.0 1 1 1 24Li 0.0000000000000000 0.000000000 1.30000000 25 1.0 1 1 1 26H 0.0000000000000000 0.000000000 0.00000000 27%EOF% 28 29####################################################################### 30# DALTON INPUT 31####################################################################### 32cat > cc3_noddy_hyp.dal <<%EOF% 33**DALTON 34.RUN WAVE FUNCTIONS 35**INTEGRAL 36.DIPLEN 37**WAVE FUNCTION 38.CC 39*SCF INPUT 40.THRESH 41 1.0D-12 42.NODIIS 43*ORBITAL INPUT 44*CC INP 45.CC3 46.PRINT 47 1 48.THRENR 491.0D-14 50.THRLEQ 511.0D-14 52.MXLRV 53 1000 54.MAX IT 55 1000 56*CCLR 57.OPERATOR 58ZDIPLEN ZDIPLEN 59.FREQUE 60 3 610.00 0.10 0.20 62*CCQR 63.NOBMAT 64.OPERATOR 65ZDIPLEN ZDIPLEN ZDIPLEN 66.ORFREQ 67 2 68 0.00 0.10 69.SHGFRE 70 1 71 0.10 72*NODDY 73.XI 74.ETA 75.FMAT 76**END OF DALTON 77%EOF% 78####################################################################### 79 80 81####################################################################### 82# CHECK SCRIPT 83####################################################################### 84echo '#!/bin/ksh 85log=$1 86 87if [ `uname` = Linux ]; then 88 GREP="egrep -a" 89else 90 GREP="egrep" 91fi 92 93# SCF, MP2 and CC3 energies: 94CRIT1=`$GREP "Total SCF energy: * -7\.5744842494" $log | wc -l` 95CRIT2=`$GREP "Total MP2 energy: * -7\.5922238176" $log | wc -l` 96CRIT3=`$GREP "Total CC3 energy: * -7\.6000272811" $log | wc -l` 97TEST[1]=`expr $CRIT1 \+ $CRIT2 \+ $CRIT3` 98CTRL[1]=6 99ERROR[1]="ENERGIES NOT CORRECT" 100 101 102# CC3 polarizability : 103CRIT1=`$GREP "ZDIPLEN *.unrel.. *..0000 *ZDIPLEN *.unrel.. *..0000 * 12\.684665" $log | wc -l` 104CRIT2=`$GREP "...1000 * ..1000 * 230.56897" $log | wc -l` 105CRIT3=`$GREP "...2000 * ..2000 * 8.8105639" $log | wc -l` 106TEST[2]=`expr $CRIT1 \+ $CRIT2 \+ $CRIT3` 107CTRL[2]=3 108ERROR[2]="CC3 POLARIZABILITY NOT CORRECT" 109 110# CC3 hyperpolarizability : 111CRIT1=`$GREP "ZDIPLEN .unrel.. *..0000 ZDIPLEN .unrel.. *..0000 ZDIPLEN .unrel.. *..0000 * -1050.3968" $log | wc -l` 112CRIT2=`$GREP "..0000 * ..1000 * ...1000 * -401811.17" $log | wc -l` 113CRIT3=`$GREP "..2000 * ..1000 * ...1000 * -115.90926" $log | wc -l` 114TEST[3]=`expr $CRIT1 \+ $CRIT2 \+ $CRIT3` 115CTRL[3]=3 116ERROR[3]="CC3 HYPERPOLARIZABILITY NOT CORRECT" 117 118 119PASSED=1 120for i in 1 2 3 121do 122 if [ ${TEST[i]} -ne ${CTRL[i]} ]; then 123 echo "${ERROR[i]} ( test = ${TEST[i]}; control = ${CTRL[i]} ); " 124 PASSED=0 125 fi 126done 127 128if [ $PASSED -eq 1 ] 129then 130 echo TEST ENDED PROPERLY 131 exit 0 132else 133 echo THERE IS A PROBLEM 134 exit 1 135fi 136 137' > cc3_noddy_hyp.check 138####################################################################### 139