An algorithm for enhanced configurational sampling in molecular dynamics simulations based on stochastic scaling of the intermolecular potential is described. It is implemented in simulations of homogeneous nucleation of a supercooled molecular liquid, carbon dioxide, resulting in significant savings of computer time. This method preserves a canonical form for the potential energy distribution and does not alter average dynamical properties for equilibrium systems. [S1063-651X(99)06603-9].