Efficient $d$-dimensional molecular dynamics simulations for studies of the glass-jamming transition

Abstract

We develop an algorithm suitable for parallel molecular dynamics simulations in $d$ spatial dimensions and describe its implementation in C++. All routines work in arbitrary $d$; the maximum simulated $d$ is limited only by available computing resources. These routines include several that are particularly useful for studies of the glass/jamming transition, such as SWAP Monte Carlo and FIRE energy minimization. Scaling of simulation runtimes with the number of particles $N$ and number of simulation threads $n_{\mathrm{threads}}$ is comparable to popular MD codes such as LAMMPS. The efficient parallel implementation allows simulation of systems that are much larger than those employed in previous high-dimensional glass-transition studies. As a demonstration of the code's capabilities, we show that supercooled $d=6$ liquids can possess dynamics that are substantially more heterogeneous and experience a breakdown of the Stokes-Einstein relation that is substantially stronger than previously reported, owing at least in part to the much smaller system sizes employed in earlier simulations.

References in corpus (11)

• Theoretical perspective on the glass transition and amorphous materials
• Strong scaling of general-purpose molecular dynamics simulations on GPUs
• Fractal free energy landscapes in structural glasses
• Universal nature of particle displacements close to glass and jamming transitions
• Packing Hyperspheres in High-Dimensional Euclidean Spaces
• Hopping and the Stokes-Einstein relation breakdown in simple glass formers
• Finite size effects in the dynamics of glass-forming liquids
• Dimensional Study of the Caging Order Parameter at the Glass Transition
• Glass transition of hard spheres in high dimensions
• Geometrical Frustration: A Study of 4d Hard Spheres
• The role of the attractive forces in a supercooled liquid