Abstract: Uniform-acceptance force-bias Monte Carlo (fbMC) methods have been shown to be a powerful technique to access longer timescales in atomistic simulations allowing, for example, phase transitions and growth. Recently, a new fbMC method, the time-stamped force-bias Monte Carlo (tfMC) method, was derived with inclusion of an estimated effective timescale; this timescale, however, does not seem able to explain some of the successes the method. In this contribution, we therefore explicitly quantify the effective timescale tfMC is able to access for a variety of systems, namely a simple single-particle, one-dimensional model system, the Lennard-Jones liquid, an adatom on the Cu(100) surface, a silicon crystal with point defects and a highly defected graphene sheet, in order to gain new insights into the mechanisms by which tfMC operates. It is found that considerable boosts, up to three orders of magnitude compared to molecular dynamics, can be achieved for solid state systems by lowering of the apparent activation barrier of occurring processes, while not requiring any system-specific input or modifications of the method. We furthermore address the pitfalls of using the method as a replacement or complement of molecular dynamics simulations, its ability to explicitly describe correct dynamics and reaction mechanisms, and the association of timescales to MC simulations in general.

Abstract: We calculate bubble nucleation rates in a Lennard-Jones fluid through explicit molecular dynamics simulations. Our approach-based on a recent free energy method (dubbed reweighted Jarzynski sampling), transition state theory, and a simple recrossing correction-allows us to probe a fairly wide range of rates in several superheated and cavitation regimes in a consistent manner. Rate predictions from this approach bridge disparate independent literature studies on the same model system. As such, we find that rate predictions based on classical nucleation theory, direct brute force molecular dynamics simulations, and seeding are consistent with our approach and one another. Published rates derived from forward flux sampling simulations are, however, found to be outliers. This study serves two purposes: First, we validate the reliability of common modeling techniques and extrapolation approaches on a paradigmatic problem in materials science and chemical physics. Second, we further test our highly generic recipe for rate calculations, and establish its applicability to nucleation processes.

Abstract: We report on the ordering and dynamics of interacting colloidal particles confined by a parabolic potential. By means of Brownian dynamics simulations, we find that by varying the magnitude of the trap stiffness, it is possible to control the dimension of the system and, thus, explore both the structural transitions and the long-time self-diffusion coefficient as a function of the degree of confinement. We particularly study the structural ordering in the directions perpendicular and parallel to the confinement. Further analysis of the local distribution of the first-neighbors layer allows us to identify the different structural phases induced by the parabolic potential. These results are summarized in a structural state diagram that describes the way in which the colloidal suspension undergoes a structural re-ordering while increasing the confinement. To fully understand the particle dynamics, we take into account hydrodynamic interactions between colloids; the parabolic potential constricts the available space for the colloids, but it does not act on the solvent. Our findings show a non-linear behavior of the long-time self-diffusion coefficient that is associated to the structural transitions induced by the external field. (C) 2015 AIP Publishing LLC.

Abstract: We study the structure and phase behavior of fullerene-cubane C60·C8H8 by Monte Carlo simulation. Using a simple potential model capturing the icosahedral and cubic symmetries of its molecular constituents, we reproduce the experimentally observed phase transition from a cubic to an orthorhombic crystal lattice and the accompanying rotational freezing of the C60 molecules. We elaborate a scheme to identify the low-temperature orientations of individual molecules and to detect a pattern of orientational ordering similar to the arrangement of C60 molecules in solid C60. Our configuration of orientations supports a doubled periodicity along one of the crystal axes.