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.

Any fracture process ultimately involves the rupture of atomic bonds. Processes at the atomic scale therefore critically influence the toughness and overall fracture behavior of materials. Atomistic simulation methods including large-scale molecular dynamics simulations with classical potentials, density functional theory calculations and advanced concurrent multiscale methods have led to new insights e.g. on the role of bond trapping, dynamic effects, crack-microstructure interactions and chemical aspects on the fracture toughness and crack propagation patterns in metals and ceramics. This review focuses on atomistic aspects of fracture in crystalline materials where significant advances have been achieved over the last ten years and provides an outlook on future perspectives for atomistic modelling of fracture.

Mechanical response of metal nanowires has recently attracted a lot of interest due to their ultra-high strengths and unique deformation behaviours. Atomistic simulations have predicted that face-centered cubic metal nanowires deform in different modes depending on the orientation between wire axis and loading direction. Here we report, by combination of in situ transmission electron microscopy and molecular dynamic simulation, the conditions under which particular deformation mechanisms take place during the uniaxial loading of [110]-oriented Au nanowires. Furthermore, by performing cyclic uniaxial loading, we show reversible plastic deformation by twinning and consecutive detwinning in tension and compression, respectively. Molecular dynamics simulations rationalize the observed behaviours in terms of the orientation-dependent resolved shear stress on the leading and trailing partial dislocations, their potential nucleation sites and energy barriers. This reversible twinning–detwinning process accommodates large strains that can be beneficially utilized in applications requiring high ductility in addition to ultra-high strength.