RUMD is a general purpose, high-performance molecular dynamics (MD) simulation package running on graphical processing units (GPU's). RUMD addresses the challenge of utilizing the many-core nature of modern GPU hardware when simulating small to medium system sizes (roughly from a few thousand up to hundred thousand particles). It has a performance that is comparable to other GPU-MD codes at large system sizes and substantially better at smaller sizes. RUMD is open-source and consists of a library written in C++ and the CUDA extension to C, an easy-to-use Python interface, and a set of tools for set-up and post-simulation data analysis. The paper describes RUMD's main features, optimizations and performance benchmarks.

This review is focused on recent progress achieved by mesoscopic simulations of plastic deformation. The methods presently available for discretizing the dislocation lines are critically discussed with emphasis on a new lattice-based model. Progress in large-scale simulations is represented by a study on the influence of long range elastic stresses on the formation of dislocation patterns in fcc crystals. A hybrid discrete-continuum method that provides an exact treatment of the boundary conditions is described and illustrated by an investigation of the critical conditions for dislocation motion in the channels of γ /γ superalloys.

We propose a four-component united-atom molecular model of bitumen. The model includes realistic chemical constituents and introduces a coarse graining level that suppresses the highest frequency modes. Molecular dynamics simulations of the model are carried out using graphic-processor-units based software in time spans in order of microseconds, which enables the study of slow relaxation processes characterizing bitumen. This paper also presents results of the model dynamics as expressed through the mean-square displacement, the stress autocorrelation function, and rotational relaxation. The diffusivity of the individual molecules changes little as a function of temperature and reveals distinct dynamical time scales. Different time scales are also observed for the rotational relaxation. The stress autocorrelation function features a slow non-exponential decay for all temperatures studied. From the stress autocorrelation function, the shear viscosity and shear modulus are evaluated, showing a viscous response at frequencies below 100 MHz. The model predictions of viscosity and diffusivities are compared to experimental data, giving reasonable agreement. The model shows that the asphaltene, resin, and resinous oil tend to form nano-aggregates. The characteristic dynamical relaxation time of these aggregates is larger than that of the homogeneously distributed parts of the system, leading to strong dynamical heterogeneity.