Specifying Occupancies For Parallel Finite-difference Frequency Calculations
One always needs to be aware that use of the OCCUPATION keyword implies some symmetry (or, most precisely, the number of symmetry species in the point group used in the calculation). More often than not, the SCF procedure will converge to the lowest state, but there are sometimes problems; there are sometimes cases where one wants to calculate some state other than the ground state, etc. Now, when one does harmonic frequency calculations with VIB=FINDIF, this generally leads to the situation in which the displaced geometries (used in the energy or gradient calculations) sometimes have lower symmetry than the reference structure. That is, whatever OCCUPATION specification works for the equilibrium structure will not work for all of the displaced points, except in the special case that the molecule has no symmetry at all. Have a means to run such a calculation is useful when 1) the default SCF doesn't hit the correct state (at whatever geometry); and 2) special cases like Renner-Teller or Jahn-Teller molecules.
It is possible to specify different occupations for different points when one uses VIB=FINDIF in conjunction with FREQ_ALGORITHM=PARALLEL, this sort of calculation now being amply documented on this website. To do this, one must use a ZMAT file with an additional %fdocc record, similar to that below:
you make a mistake, you can simply edit your zmatxxx files and insert the correct OCCUPATION information the old fashioned way...