An approach for carrying out vibrationally adiabatic mixed quantum-classical molecular dynamics simulations is presented. An appropriate integration scheme is described for the vibrationally adiabatic equations of motion of a diatomic solute in a monatomic solvent and an approach for calculating the adiabatic energy levels is presented. Specifically, an iterative Lanczos algorithm with full reorthogonalization is used to solve for the lowest few vibrational eigenvalues and eigenfunctions. The eigenfunctions at one time step in a mixed quantum-classical trajectory are used to initiate the Lanczos calculation at the next time step. The basis set size is reduced by using a potential-optimized discrete variable representation. As a demonstration the problem of a homonuclear diatomic molecule in a rare gas fluid (N-2 in Ar) has been treated. The approach is shown to be efficient and accurate. An important advantage of this approach is that it can be straightforwardly applied to polyatomic solutes that have multiple vibrational degrees-of-freedom that must be quantized. (C) 2003 American Institute of Physics.