Influence of thermostats and carrier gas on simulations of nucleation.

We investigate the influence of carrier gas and thermostat on molecular dynamics (MD) simulations of nucleation. The task of keeping the temperature constant in MD simulations is not trivial and an inefficient thermalization may have a strong influence on the results. Different thermostating mechanisms have been proposed and used in the past. In particular, we analyze the efficiency of velocity rescaling, Nose-Hoover, and a carrier gas (mimicking the experimental situation) by extensive MD simulations. Since nucleation is highly sensitive to temperature, one would expect that small variations in temperature might lead to differences in nucleation rates of up to several orders of magnitude. Surprisingly, the results indicate that the choice of the thermostating method in a simulation does not have--at least in the case of Lennard-Jones argon--a very significant influence on the nucleation rate. These findings are interpreted in the context of the classical theory of Feder et al. [Adv. Phys. 15, 111 (1966)] by analyzing the temperature distribution of the nucleating clusters. We find that the distribution of cluster temperatures is non-Gaussian and that subcritically sized clusters are colder while postcritically sized clusters are warmer than the bath temperature. However, the average temperature of all clusters is found to be always higher than the bath temperature.