  We study the 38-atom Lennard-Jones cluster with parallel tempering Monte
Carlo methods in the microcanonical and molecular dynamics ensembles. A new
Monte Carlo algorithm is presented that samples rigorously the molecular
dynamics ensemble for a system at constant total energy, linear and angular
momenta. By combining the parallel tempering technique with molecular dynamics
methods, we develop a hybrid method to overcome quasi-ergodicity and to extract
both equilibrium and dynamical properties from Monte Carlo and molecular
dynamics simulations. Several thermodynamic, structural and dynamical
properties are investigated for LJ$_{38}$, including the caloric curve, the
diffusion constant and the largest Lyapunov exponent. The importance of
insuring ergodicity in molecular dynamics simulations is illustrated by
comparing the results of ergodic simulations with earlier molecular dynamics
simulations.
