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.

No abstract

Diamond is the hardest material on Earth. Nevertheless, polishing diamond is possible with a process that has remained unaltered for centuries and is still used for jewellery and coatings: the diamond is pressed against a rotating disc with embedded diamond grit. When polishing polycrystalline diamond, surface topographies become non-uniform because wear rates depend on crystal orientations. This anisotropy is not fully understood and impedes diamond's widespread use in applications that require planar polycrystalline films, ranging from cutting tools to confinement fusion. Here, we use molecular dynamics to show that polished diamond undergoes an sp(3)-sp(2) order-disorder transition resulting in an amorphous adlayer with a growth rate that strongly depends on surface orientation and sliding direction, in excellent correlation with experimental wear rates. This anisotropy originates in mechanically steered dissociation of individual crystal bonds. Similarly to other planarization processes, the diamond surface is chemically activated by mechanical means. Final removal of the amorphous interlayer proceeds either mechanically or through etching by ambient oxygen.

The ultrasmoothness of diamond-like carbon coatings is explained by an atomistic/continuum multiscale model. At the atomic scale, carbon ion impacts induce downhill currents in the top layer of a growing film. At the continuum scale, these currents cause a rapid smoothing of initially rough substrates by erosion of hills into neighboring hollows. The predicted surface evolution is in excellent agreement with atomic force microscopy measurements. This mechanism is general, as shown by similar simulations for amorphous silicon. It explains the recently reported smoothing of multilayers and amorphous transition metal oxide films and underlines the general importance of impact-induced downhill currents for ion deposition, polishing, and nanopattering.

Materials performance in structural applications is often restricted by a transition from ductile response to brittle fracture with decreasing temperature. This transition is currently viewed as being controlled either by dislocation mobility or by the nucleation of dislocations. Fracture experiments on tungsten single crystals reported here provide evidence for the importance of dislocation nucleation for the fracture toughness in the semibrittle regime. However, it is shown that the transition itself, in general, is controlled by dislocation mobility rather than by nucleation.

No abstract