A normal mode coordinate-based molecular optimization algorithm was implemented and its performance tested against other optimization techniques. In certain cases the method was found to be computationally simpler and numerically more stable than the optimization algorithms based on Cartesian or internal valence coordinates. The usual redundant/internal coordinate scheme provided fastest convergence for compact covalently bonded molecules, while the normal mode method was found to be more suitable for more weakly bonded molecular complexes. For constrained optimizations use of the normal coordinates allows one to naturally separate the lower-energy modes from those more typically studied with vibrational spectroscopy. Thus, it provides an appropriate tool for simulations of IR and Raman spectra of larger molecules and complex systems when specific conformations are desired.