The technique of H-2 MAS NMR spectroscopy is presented for the investigation of multiaxial molecular dynamics. To evaluate the effects of discrete random reorientation a Lie algebraic formalism based on the stochastic Liouville-von Neumann equation is developed. The solution to the stochastic Liouville-von Neumann equation is obtained both in the presence and absence of rf irradiation. This allows effects of molecular dynamics to be evaluated during rf pulses and extends the applicability of the formalism to arbitrary multiple pulse experiments. Theoretical methods are presented for the description of multiaxial dynamics with particular emphasis ion the application of vector parameters to represent molecular rotations. Numerical time and powder integration algorithms are presented that are both efficient and easy to implement computationally. The applicability of H-2 MAS NMR spectroscopy for investigating molecular dynamics is evaluated from theoretical spectra. To demonstrate the potential of the technique the dynamics of thiourea-H-2(4) is investigated experimentally. From a series of variable temperature MAS and quadrupole echo spectra it has been found that the dynamics can be described by composite rotation about the CS and CN bonds. Both experiments are sensitive to the fast CS rotation which is shown to be described by the Arrhenius parameters E-CS = 46.4 +/- 2.3 kJ mol(-1) and ln(A(CS)) = 32.6 +/- 0.9. The MAS experiment represents a significant improvement by simultaneously allowing the dynamics of the slow CN rotation to be fully characterized in terms of E-CN = 56.3 +/- 3.4 kJ mol(-1) and ln(A(CN))= 25.3 +/- 1.1. (C) 1998 Elsevier Science B.V. All rights reserved.