[molpro-user] frequency calc crashes at start of hessian
Benj FitzPatrick
benjfitz at uchicago.edu
Thu Jul 16 06:20:56 BST 2009
Hello all,
I'm trying to calculate the ub3lyp geometry and frequencies of the first
electronically excited state of CH3SO2 (the molecule has Cs symmetry, X = A', A
= A"). The frequency calculation (input, output, and molpro.log are below)
crashes during the first calculation of the hessian (either when I try to
calculation the frequency by itself or at the end of the geometry optimization).
I've tried 2006.1 (patch level 82) and 2008.1 (patch level 31), and both
exhibit the same problem. Any suggestions would be greatly appreciated.
Thanks,
Benj FitzPatrick
University of Chicago
---.inp file---
***,trying to look @ 1st excited state, A"***
memory,200,M
gthresh,oneint=1.d-12,twoint=1.d-12,zero=1.d-12
CH2 = 1.09646054
SC3 = 1.79992022
SCH3 = 107.87658656
OS4 = 1.54328295
OSC4 = 103.43622747
DIH4 = -47.86474777
HC6 = 1.09889138
HCS6 = 107.89395182
DIH6 = -72.85245290
geometry={
ang;
h;
c, 1, ch2;
s, 2, sc3, 1, sch3;
o, 3, os4, 2, osc4, 1, dih4;
o, 3, os4, 2, osc4, 1, -dih4;
h, 2, hc6, 3, hcs6, 5, dih6;
h, 2, hc6, 3, hcs6, 4, -dih6;
}
basis=6-31G*;
{uhf;
wf,41,2,1}
basis=6-311G**;
{uhf;
wf,41,2,1}
basis,default=aug-cc-pvdz,s=avtz+d;
{uhf;
wf,41,2,1}
{uks,b3lyp;
wf,41,2,1}
{freq;
save,5300.2}
put,molden,ch3so2-A-opt-ub3lyp-avtdz-d2.molden
---.out file---
PROGRAM * HESSIAN
*** Long output written to logfile /comp/ab/benjfitz/molpro.log ***
Computing numerical hessian using default procedure for command UKS-SCF
Symmetry elements for unique displacements: Y
Using no symmetry in wavefunction calculations
Running default procedure: UKS-SCF003
Recomputing wavefunction at reference geometry
?ERROR IN SCFOCC: NO ORBITALS IN SYMMETRY 2
GLOBAL ERROR fehler on processor 0
0:0:fehler on processor 0:: 0
(rank:0 hostname:benjfitz pid:460):ARMCI DASSERT fail. armci.c:ARMCI_Error():260
cond:0
0: ARMCI aborting 0 (0).
tmp =
/home/benjfitz/pdir//cluster/molprop-2008.1_10_4/molprop_2008_1_Linux_x86_64_i8.exe.p
Creating: host=benjfitz, user=benjfitz,
file=/cluster/molprop-2008.1_10_4/molprop_2008_1_Linux_x86_64_i8.exe,
port=57515
---molpro.log---
PROGRAM * UKS-SCF (Unrestricted Kohn-Sham) Authors: W. Meyer, H.-J. Werner,
P.J. Knowles, 1993
NUMBER OF ELECTRONS: 21+ 20- SPACE SYMMETRY=2 SPIN SYMMETRY=Doublet
CONVERGENCE THRESHOLDS: 1.00E-05 (Density) 1.00E-06 (Energy)
MAX. NUMBER OF ITERATIONS: 60
INTERPOLATION TYPE: DIIS
INTERPOLATION STEPS: 2 (START) 1 (STEP)
LEVEL SHIFTS: 0.30 (ALPHA) 0.30 (BETA)
Use fine grid record 1800.2 for energy threshold 1.0E-06
Molecular orbitals read from record 2203.2 Type=UKS/ALPHA (state 1.2)
Present basis size: 151
Previous basis size: 91 60
Previous occupation: 14 7
Previous closed-shells: 14 6
Previous core-orbitals: 6 2
Present occupation: 21
Present closed-shells: 20
Present core-orbitals: 8
Wavefunction symmetry: 1 Doublet
Orbital overlap <old|new>: < 1.1| 1.1> = 1.000000 < 2.1| 2.1> = 1.000000 <
3.1| 3.1> = 1.000000 < 4.1| 4.1> = 1.000000
< 5.1| 5.1> = 1.000000 < 6.1| 6.1> = 1.000000 <
1.2| 7.1> = 1.000000 < 2.2| 8.1> = 1.000000
< 7.1| 9.1> = 1.000000 < 8.1|10.1> = 1.000000 <
9.1|11.1> = 1.000000 <10.1|12.1> = 1.000000
<11.1|13.1> = 1.000000 <12.1|14.1> = 1.000000
<13.1|15.1> = 1.000000 <14.1|16.1> = 1.000000
< 3.2|17.1> = 1.000000 < 4.2|18.1> = 1.000000 <
5.2|19.1> = 1.000000 < 6.2|20.1> = 1.000000
< 7.2|21.1> = 1.000000
More information about the Molpro-user
mailing list