We propose a novel strategy to study numerically the glass transition of molecular fluids. Our approach combines the construction of simple yet realistic models with the development of Monte Carlo algorithms to accelerate equilibration and sampling. Inspired by the well-studied ortho-terphenyl glass-former, we construct a molecular model with an analogous triangular geometry and construct a `flip' Monte Carlo algorithm. We demonstrate that the flip Monte Carlo algorithm achieves a sampling speedup of about 10⁹ at the experimental glass transition temperature Tg. This allows us to systematically analyze the equilibrium structure and molecular dynamics of the model over a temperature regime previously inaccessible. We carefully compare the observed physical behavior to earlier studies that used atomistic models. In particular, we find that the glass fragility and the departure from the Stokes-Einstein relation are much closer to experimental observations. We characterize the development and temperature evolution of spatial correlations in the relaxation dynamics, using both orientational and translational degrees of freedom. Excess wings emerge at intermediate frequencies in dynamic rotational spectra, and we directly visualize the corresponding molecular motion near Tg. Our approach can be generalized to a number of molecular geometries and paves the way to a complete understanding of how molecular details may affect more universal physical aspects characterizing molecular liquids approaching their glass transition.
Simon et al. (Thu,) studied this question.