We introduce a simple local atomic structure optimization algorithm which is significantly faster than standard implementations of the conjugate gradient method and often competitive with more sophisticated quasi-Newton schemes typically used in ab initio calculations. It is based on conventional molecular dynamics with additional velocity modifications and adaptive time steps. The surprising efficiency and especially the robustness and versatility of the method is illustrated using a variety of test cases from nanoscience, solid state physics, materials research, and biochemistry.

Graphite and other lamellar materials are used as dry lubricants for macroscale metallic sliding components and high-pressure contacts. It has been shown experimentally that monolayer graphene exhibits higher friction than multilayer graphene and graphite, and that this friction increases with continued sliding, but the mechanism behind this remains subject to debate. It has long been conjectured that the true contact area between two rough bodies controls interfacial friction. The true contact area, defined for example by the number of atoms within the range of interatomic forces, is difficult to visualize directly but characterizes the quantity of contact. However, there is emerging evidence that, for a given pair of materials, the quality of the contact can change, and that this can also strongly affect interfacial friction. Recently, it has been found that the frictional behaviour of two-dimensional materials exhibits traits unlike those of conventional bulk materials. This includes the abovementioned finding that for few-layer two-dimensional materials the static friction force gradually strengthens for a few initial atomic periods before reaching a constant value. Such transient behaviour, and the associated enhancement of steady-state friction, diminishes as the number of two-dimensional layers increases, and was observed only when the two-dimensional material was loosely adhering to a substrate. This layer-dependent transient phenomenon has not been captured by any simulations. Here, using atomistic simulations, we reproduce the experimental observations of layer-dependent friction and transient frictional strengthening on graphene. Atomic force analysis reveals that the evolution of static friction is a manifestation of the natural tendency for thinner and less-constrained graphene to re-adjust its configuration as a direct consequence of its greater flexibility. That is, the tip atoms become more strongly pinned, and show greater synchrony in their stick-slip behaviour. While the quantity of atomic-scale contacts (true contact area) evolves, the quality (in this case, the local pinning state of individual atoms and the overall commensurability) also evolves in frictional sliding on graphene. Moreover, the effects can be tuned by pre-wrinkling. The evolving contact quality is critical for explaining the time-dependent friction of configurationally flexible interfaces.

The atomic-scale details during melting of a surface-free Lennard-Jones crystal were monitored using molecular dynamics simulations. Melting occurs when the superheated crystal spontaneously generates a sufficiently large number of spatially correlated destabilized particles that simultaneously satisfy the Lindemann and Born instability criteria. The accumulation and coalescence of these internal local lattice instabilities constitute the primary mechanism for homogeneous melt nucleation inside the crystal, in lieu of surface nucleation for equilibrium melting. The vibrational and elastic lattice instability criteria as well as the homogeneous nucleation theory all coincide in determining the superheating limit.