The coupled Schrodinger equations that describe nonadiabatic dynamics are recast using the Bohm formulation of quantum mechanics. The resulting coupled Bohm equations are solved numerically for two scattering models, giving results that are essentially identical to wave-packet solution of the original coupled Schrodinger equations. The classical limit of the set of coupled Bohm equations is then described, producing a mixed quantum-classical theory incorporating classical-like motion on each potential energy surface accompanied by quantum transitions between the quantum states. Numerical tests of the mixed quantum-classical method are in excellent agreement with the accurate full-quantum results for the model problems. The method is contrasted with the related surface-hopping approach. It is shown that computing the dynamics of a distribution of classical particles is more consistent and more accurate than computing the motion of independent point particles as with surface hopping. (C) 2000 American Institute of Physics. [S0021-9606(00)00814-X].