Computer simulations in the study of gold nanowires: the effect of impurities

Suspended gold nanowires have recently been made in an ultra-high vacuum ambient and were imaged by electron microscopy. Two puzzles were presented: one atom thick wires are produced and some of the atomic distances between these atoms before their breaking were too large. Simulations using realistic molecular dynamics method were able to unveil some processes to explain the mechanisms of formation, evolution, and breaking of these atomically thin Au nanowires under stress. The calculations showed how defects induce the formation of constrictions that eventually will form the one-atom chains. Atomically thin chains, five atoms long were obtained, before breaking. The results were in excellent agreement with experimental results except for the large Au-Au distances. In fact no theoretical calculation of pure gold nanowires have been able to produce such large distances. Light impurities that cannot be imaged in these experiments may be responsible for these large Au-Au distances. Using ab initio total energy calculations based on the density functional theory, we have studied the effect of H, C, O, N, B, S, CH, CH2, and H-2 impurities on the nanowire's electronic and structural properties, in particular how they affect the rupture of the nanowire. We find that the impurities tend to locally increase the nanowire's strength, in such a way that its rupture always occurs at an Au-Au bond and never at an Au-X bond (X being an impurity). In particular, oxygen seems to form very stable bonds that may be used to pull longer Au chains. Regarding the observed large Au-Au bond lengths, it was found, based on quasi-static calculations, that the best candidate to explain the large distances is H. However, some particular experimental conditions may lead to different results.