Two semiclassical, initial value representation (IVR) treatments are presented for the correlation function , where psi((i)) and psi((f )), are energy eigenfunctions of a "zero-order" Hamiltonian describing an arbitrary, integrable, vibrational system. These wave functions are treated semiclassically so that quantum calculations and numerical integrations over these states are unnecessary. While one of the new approximations describes the correlation function as an integral over all phase space variables of the system, in a manner similar to most existing IVR treatments, the second approximation describes the correlation function as an integral over only half of the phase space variables (i.e., the angle variables for the initial system). The relationship of these treatments to the conventional Herman-Kluk approximation for correlation functions is discussed. The accuracy and convergence of these treatments are tested by calculations of absorption spectra for model systems having up to 18 degrees of freedom, using Monte Carlo techniques to perform the multidimensional phase space integrations. Both treatments are found to be capable of producing spectra of excited, anharmonic states that agree well with quantum results. Although generally less accurate than full phase space or Herman-Kluk treatments, the half phase space method is found to require far fewer trajectories to achieve convergence. In addition, this number is observed to increase much more slowly with the system size than it does for the former methods, making the half-phase space technique a very promising method for the treatment of large systems. (C) 2004 American Institute of Physics.