***, hydroperoxyl HO2 splitting instanton
if(NPROC_MPP.gt.1) then
skipped
endif
FILE,2,ho2inst.wfu,new

proc beadpot={
{uhf;wf,sym=1,spin=1}
}

proc beadgrad={
beadpot
force
}

proc beadfreq1={
frequencies
}
proc beadfreq3={
beadpot
frequencies
}

geometry={ ! optimized geometry at well minimum
    3
 UHF-SCF002/VDZ  ENERGY=-150.16631226 ! comment line
H 0.904294157899999895 0.919502897200000113 0
O 0.681281911699999965 -0.00702956160000003938 0
O -0.683576069499999939 0.00706386439999998061 0
}

mass,isotope,print
basis=vdz

instanton, splitting, input='splitting_initial.xyz', \\
   output='splitting_final.xyz', maxit=200
table,action,fluctratio,tunsplit*1e9