Open Access is an initiative that aims to make scientific research freely available to all. To date our community has made over 100 million downloads. It’s based on principles of collaboration, unobstructed discovery, and, most importantly, scientific progression. As PhD students, we found it difficult to access the research we needed, so we decided to create a new Open Access publisher that levels the playing field for scientists across the world. How? By making research easy to access, and puts the academic needs of the researchers before the business interests of publishers.
We are a community of more than 103,000 authors and editors from 3,291 institutions spanning 160 countries, including Nobel Prize winners and some of the world’s most-cited researchers. Publishing on IntechOpen allows authors to earn citations and find new collaborators, meaning more people see your work not only from your own field of study, but from other related fields too.
To purchase hard copies of this book, please contact the representative in India:
CBS Publishers & Distributors Pvt. Ltd.
www.cbspd.com
|
customercare@cbspd.com
INFIQC, Departamento de Matemática y Física, Argentina
Patricio Vélez
INFIQC, Departamento de Matemática y Física, Argentina
Sergio A. Dassie
Departamento de Fisicoquímica, Facultad de Ciencias Químicas, Universidad Nacional de Córdoba, Argentina
Ezequiel P. M. Leiva
INFIQC, Departamento de Matemática y Física, Argentina
*Address all correspondence to:
1. Introduction
“The era in which the number of transistors on a computer chip doubles at a constant rate is drawing to a close”. This is not the prophecy of an obscure mind, but is more or less the conclusion drawn by none other than the man who coined Moores´ law
Times Online September 19, 2007, http://technology.timesonline.co.uk/tol/news/tech_and_web/article2489053.ece
. In an interview held in 2007, Gordon Moore recognized that by about 2020, his law would come up against the laws of physics. Furthermore, he recognized a change in a paradigm: the replacement of the top-down approach currently used for building circuits by a bottom-up procedure, where chips would be assembled using individual atoms or molecules. This is nothing but the realm of nanotechnology, while there is some consensus that the elementary switches of these circuits should be molecules with some feature allowing for the on/off status required for the components of logical devices, many questions remain concerning their stability. In the case of micrometric circuit components temperature may be an issue, but in the case of single molecules thermal effects may be overwhelming, since current flow occurs across a single bond. The lifetime of this bond, will determine the lifetime of the circuit component. Under these conditions, circuit engineering will be coming unexpectedly close to chemical kinetics.
It still is far from clear which will be the technological procedure for the massive production of these molecular circuits. However, there are a number of experimental techniques for the study of their properties that are well established. These are shown schematically in Fig. 1. Fig. 1d shows a method devised to study the structure of monatomic nanowires (NWs). It has been developed by Kondo and Takayanagi (Kondo & Takayanagi, 1997) using High Resolution Transmission Electronic Microscopy (HRTEM) and allows the generation of suspended NWs. In this approach nanowires are generated in situ by focusing an electron beam on adjacent sites of a self-supported metal thin film (ca. 3 nm), making holes and allowing them to grow until a nanometric bridge is formed inside or between grains. The relatively long lifetimes of these metallic nanowires (of the order of seconds and even minutes) allows for a detailed characterization of their geometries. This method has even been used to study the structure of alloyed metal monatomic nanowires (Bettini et al., 2006). However, because of the preparation procedure, these nanowires are susceptible to contamination by light impurity atoms (Galvão et al., 2004, Legoas et al., 2002, Legoas et al., 2004), not observed directly in the HRTEM, that affect their structure upon insertion in the monatomic chain.
Figure 1.
Different methods employed to generate suspended nanowires: a) metallic jump-to-contact, b) mechanically controllable break-junction, c) electrochemical growth, d) electron beam punching, and e) molecular jump-to-contact.
Another procedure that has been used to study the properties of monatomic metal contacts is the so-called Mechanically Controllable Break-Junctions method (MCBJ, Fig. 1 b). In this technique, a metal wire is mechanically broken by mounting its two ends on a bendable substrate. In this way, the length of the contact can be adjusted by controlling the bending of the substrate. Lithographic designed MCBJ provides highly stable electrodes, with standard deviations of less than 1 pm over 24 h at low temperatures. A detailed discussion on this technique has been given by van Ruitenbeek et al. (van Ruitenbeek et al., 2005). The metal neck at the MCBJ can be also tuned by electrochemical methods as long as the metal involved can be deposited electrochemically. The main drawbacks of this technique are the lack of information on the atomic configuration, and the fact that both electrodes must be made of the same metal. The former situation was challenged recently by Yoshida et al. (Yoshida et al., 2007), who developed MCBJ inside a HRTEM chamber that allows for simultaneous performance of conductance measurements and electronic imaging of the atomic chain.
In the case of metallic Scanning Tunnelling Microscopy (STM) and conducting Atomic Force Microscopy (AFM) break junctions (Ohnishi et al., 1998, Rubio-Bollinger et al., 2001, Xu et al., 2003a, Xu & Tao, 2003), the sharp tip of the scanning probe microscope is approached to the surface and brought into contact with it (Fig. 1a). The surface may be made either of the same material as the tip or of a different one containing islands of the tip material. Upon retraction, an atomic chain of atoms occurs, whose properties are analyzed. When employed in an Ultra High Vacuum (UHV). The main drawbacks of this technique are the need for an intensive cleaning procedure and the implementation of thermal and mechanical stability conditions, as it is the case of STM. On the other hand, this methodology can be used in situ for electrochemical experiments. In fact, it has been employed by Tao and co-workers to perform conductance measurements with Au nanowires under a variety of experimental electrochemical conditions (He et al., 2002, Xu et al., 2003b), including the properties of the nanojunction when molecules are inserted into it (Xu et al., 2003a; Xu & Tao, 2003). Efficient electrochemical methods for the generation of nanowires have also been developed. Li and Tao(Li & Tao, 1998) managed to bridge the gap between an STM tip and a substrate by a suitable potential control of both the tip and the surface(Fig. 1c). In a similar procedure, the gap between two supported metal pieces was filled by a nanometric wire using an adequate feedback electronic setup (Li et al, 1999). The latter procedure has the advantage of removing thermal drift problems.
Finally, the method developed Haiss et al. moves more into the spirit of bottom up nanostructuring (Fig. 1e). In this procedure, a molecule bridges spontaneously the gap between an STM tip and a surface (Haiss et al., 2004, Haiss et al., 2006). The substrate-tip connection is verified by jumps in the tunnelling current measured.
2. Structure and stability of pure and contaminated metallic monatomic nanowires
2.1. Experimental measurements
The conductance measured through metallic monatomic nanowires is quantized in units of G0 = 2e2/h (where e is the charge of an electron and h stands for Planck’s constant). The force during the fabrication and breaking of a gold nanowire was measured using a specific STM supplemented with a force sensor at room temperature (Rubio-Bollinger et al., 2001). Force and conductace curves were thus obtained simultaneously. The later displayed a steplike behavior down to a value close to one conductance quantum (G0), which corresponds to a one-atom contact, while the force curve showed a sawtooth like signal decreasing in amplitude in a sequence of elastic stages separated by sudden force relaxations. In this experiments the one-atom contact of gold was further stretched a distance of about 1 nm while the conductance remained close to G0, which signals the formation of a chain of about four atoms long that finally breaks. This corresponds to a monatomic neck of ~4 gold atoms. Similar observations were made in STM experiments at 4 K (Yanson et al., 1998), and from direct observations by means of HRTEM measurements (Ohnishi et al., 1998, Rodrigues & Ugarte, 2001a, Rodrigues & Ugarte, 2002; Rodrigues et al., 2000).
Due to the inherent irreproducibility of the contacts formed, and therefore of the measured conductance curves, it is useful to construct histograms with a few hundred measured curves. A force histogram for the Au nanowire showed a narrow distribution centered at 1.5 ± 0.3 nN for the force needed to break one single bond in the chain (Rubio-Bollinger et al., 2001).
Two striking experimental features have drawn the attention of the researches in the last few years: first, the existence of unusually large Au – Au distances in the order of 3.6 – 4.0 Å (Kizuka, 2008, Koizumi et al., 2001, Ohnishi et al., 1998, Sato et al., 2006b), longer than the bulk distance of 2.88 Å. These experimental observations along with the simulation of HRTEM measurements suggest that the atomic-sized wires are complexed with light elements (such as H, C, S for example) (Kizuka, 2008, Koizumi et al., 2001).
Time sequence of atomic resolution images of the formation, elongation and fracture of a suspended chain of gold atoms: (a) 0 s; (b) 0.64 s; (c) 1.12 s; and (d) 3.72 s. Atomic positions appear as dark lines or dots. A schematic representation of the chain structure is shown in (e) (distances are marked in Å and the error bar is 0.1 Å); the letters A and B indicate the apex position in (c). Note that the chain is attached to the tips through a two atom structure. The double arrows in (a) and (b) have been drawn to indicate the movement (rotation) of the lower apex. Reprinted figure with permission from Rodrigues, V.; Ugarte, D. Phys. Rev. B, Vol. 63, No.7, 073405 (pp 1-4), 2001. Copyright 2001 by the American Physical Society.
The stability and breakdown of gold nanojunctions at different stretching rates has recently been measured by Huang et al. (Huang et al., 2007a). Information about the lifetime of the Au – Au bond is extracted from the length of the last plateau in the conductance traces of Au – Au point contacts. The most probable stretching distance for this last plateau (corresponding to an atom-sized contact between the electrodes as verified by its conductance at a G0 value) maintains a constant value of ~0.1 nm at low stretching rates (0.8 – 8.3 nm/s). At high stretching rates (45.9 – 344 nm/s) a maximum plateau of ~0.17 nm is reached. Between these two plateaus, a linearly proportional regime is observed between the most probable distance and the logarithm of the stretching rate. This three-phase regime (illustrated in Fig. 3) has also been observed in the breakdown of biological molecules measured by AFM (Auletta et al., 2004; Merkel et al., 1999, Schönherr et al., 2000, Zou et al., 2005). The linear increase in stretching distance is related to the stretching rate according to Eq. 1 (Evans & Ritchie, 1997, , Evans 1999, Evans, 2001).
L*=kBTksxβln(toffksxβkBT)+kBTksxβlnυE1
Figure 3.
Most probable stretching distance, L*, vs logarithm of stretching rate for a Au – Au point contact. The black line is the linear fit of L* with the logarithm of the stretching rate, according to Eq 1. Reprinted with permission from Huang, Z. F.; Chen, F.; Bennett, P. A.; Tao, N. J. J. Am. Chem. Soc., Vol. 129, No.43, 13225-13231, 2007. Copyright 2007 American Chemical Society.
In Eq. 1L* is the most probable stretching distance, ks is the effective spring constant of the bond, xβ is the average thermal bond length along the pulling direction until dissociation, toff is the natural lifetime, and ν is the stretching rate. By fitting the linear regime observed for the Au – Au point contacts to the above equation, a value of toff = 81 ms was obtained.
2.2. Computational simulations
Experimental observations such as those described above, have triggered numerous theoretical approaches to gain further insight into the structure and transmission properties of these metallic nanowires. The mechanical structure and evolution of a tip-surface contact has been modelled by means of different computational techniques. In early works, Molecular Dynamics (MD) and effective medium theory potentials have been used to simulate the mechanical deformation of atomic-sized metallic contacts under tensile strain for Au, Ag, Pt, Ni (Bahn & Jacobsen, 2001; Dreher et al., 2005, Rubio-Bollinger et al., 2001, Sørensen et al., 1998). The relative utility of different semiempirical potentials for MD simulations of stretched gold nanowires has been recently reviewed (Pu et al., 2007a). In this work the authors find that the second-moment approximation of the tight-binding potential reproduces well the energetics of finite Au clusters. Also, the calculated tensile force just before the nanowire breaks is around 1.5 nN, consistent with the experimental result The stretching of gold nanowires has also been simulated in the presence of solvent (Pu et al., 2007b).
Tight binding MD simulations using the Naval Research Laboratory potentials, along with ab-initio quantum calculations within the Density Functional Theory (DFT) framework, provided evidence that a one-atom thick 5-atom long necklace is formed for gold nanowires under a stretching force (da Silva et al., 2001; da Silva et al., 2004). Before breaking, relatively long Au – Au distances, of the order of 3.0 – 3.1 Å are obtained, in agreement with experiments.
The main drawback of the MD simulations mentioned above is that the stretching rate used (typically in the order of 1 -2 m/s) is around 9 - 10 orders of magnitude larger than the experimental values (in the order of a tenths to a few nm/s).
The use of ab-initio quantum methods to study the stability of monatomic metallic nanowires is usually performed using two different methodologies. One of these procedures involves stretching the nanowire in steps of a certain elongation, minimizing the energy of the system at each step, until the breakage of the nanowire is achieved (Bahn et al., 2002, da Silva et al., 2004; De Maria & Springborg, 2000; Häkkinen et al., 2000, Nakamura et al., 1999, Nakamura et al., 2001; Novaes et al., 2003, Okamoto & Takayanagi, 1999, Rubio-Bollinger et al., 2001, Sánchez-Portal et al., 1999, Skorodumova & Simak, 2003). This method has been successfully used to, for example, predict the pulling force necessary for breaking the nanowire (da Silva et al., 2004; Nakamura et al., 1999, Novaes et al., 2003, Rubio-Bollinger et al., 2001). An alternative procedure, and more costly from a computational point of view, consists in performing Ab-initio Molecular Dynamics (AIMD), which has been employed to obtain a detailed description of the elongation process of pure gold nanowires (Hobi et al., 2008, Torres et al., 1999), and gold NWs in the presence of organic molecules (Krüger et al., 2002), or light weight elements as contaminants (Hobi et al., 2008, Legoas et al., 2004). AIMD takes into account the thermal motion of the system, for which in a first approach seems a more appropriate approximation to obtain reliable information about the mechanical proprieties and the stability of metallic atom-sized wires. Nonetheless, the experimental elongation times are far larger than those accessible by an AIMD simulation, typically in the order of nanoseconds.
As it was mentioned above, one of the possible explanations for the rather long Au – Au distances observed in a monatomic neck is the presence of light weight elements, such as C, H, O or S, intercalated between gold atoms. This possibility has only been investigated so far by the means of computational simulations, since these light weight elements would have a low contrast against the much heavier Au atoms, and thus rendering their direct visualization by today’s electronic microscopes very difficult. Carbon is a frequent contaminant in bulk gold (Legoas et al., 2002), while hydrogen (O'Hanlon, 2001) and oxygen (Bahn et al., 2002) are impurities very difficult to extract even in the best UHV conditions. Although the possibility of hydrogen acting as a contaminant has been ruled out by an AIMD study (Legoas et al., 2004), this result has been subsequently challenged (Hobi et al., 2005).
Legoas et al. (Legoas et al., 2002) modelled monatomic gold chains contaminated with carbon by means of geometry optimization at DFT-LDA level. The authors used an isolated linear chain and no tension was exerted on the system. Their results showed that long Au – Au distances of around 4 – 4.5 Å could be explained by the presence of two consecutive carbon atoms (C2) inserted into the gold chain. Whereas, another set of anomalously long bond in the order of 3 – 3.7 Å could be a consequence of a mixture of pure Au – Au bonds with contamination of such a bond by a single carbon atom.
Skorodumova and Simak (Skorodumova & Simak, 2003) showed, using DFT-GGA calculations, that the unusual structural stability of monatomic gold wires could be explained in terms of hydrogen contamination. Stretching the nanowire, the authors observed that the chain takes a linear structure with hydrogen atoms intercalated and a Au – Au distance of 8.8 Å before the nanowire breaks. The cohesive energy of the contaminated gold wire was found to be 2-fold higher than a pure chain of gold atoms. This last result was attributed to a partial charge transfer from gold to hydrogen. Subsequently, the influence of carbon was explored using the same computational methodology (Skorodumova et al., 2007, Skorodumova & Simak, 2004), also finding that carbon can enhance the stability of linear gold chains yielding large interatomic distances.
Novaes et al. (Novaes et al., 2003) studied through ab- initio calculations the effect of H, B, C, N, O and S impurities on a gold nanowire electronic and structural properties. The authors find that the most likely candidates to explain the distances in the range of 3.6 Å and 4.8 Å are H and S impurity atoms, respectively.
The main drawback of the procedures presented so far is that the presence of the impurity is simply assumed, and no description is obtained of how, when, or with what probability it migrates to the position it was assumed to have. To overcome this limitation, an AIMD study of the formation and growth of gold chains with a variety of impurities (H, C, O, S), without any assumption of their initial positions was performed (Anglada et al., 2007). One or two impurity atoms were introduced randomly in an amorphous column of 50 – 150 gold atoms. These amorphous solid columns were stretched during 4 – 18 ns until they broke. Hydrogen was always found to evaporate before formation of the monatomic chain took place. Carbon and oxygen were found in the final chains with low probability (~ 10 %), while sulphur was found participating in it with a high probability (~ 90 %). The mean distances between gold atoms bridged by C, O and S were 3.3, 3.4 and 5.0 Å, respectively, in good agreement with experiments.
Inasmuch as this last study provides a level of accuracy and reliability superior to those mentioned before, the stretching rate is still much higher than those typically used in experiments, and the simulations last only a few nanoseconds, while experiments take place in the order of 0.001 – 1 second.
The study of the effect of impurities on the structure and stability of gold nanowires is an ongoing investigation topic. Some of the most recent theoretical work can be found in (Jelínek et al., 2008;, Novaes et al., 2006, Zhang et al., 2008). Although gold is by far the most prominent element of interest in the formation of nanowires, some recent computational studies have also involved copper nanowires (Amorim et al., 2007, Amorim et al., 2008; Sato et al., 2006a).
Building a device in which a single molecule bridges two metallic electrodes is of major interest, since one could easily tailor the nanojunction electronic properties by changing the molecule or even only a substituent in the molecule. This opens an enormous range of possibilities in the field of molecular electronics.
Xu et. al. reported the first electromechanical measurement of a molecular junction (Xu et al., 2003a). The authors determined simultaneously the conductance and the force under mechanical stretching for the octanodithiol (ODT) and 4,4’-bipyridine (BYP) nanojunctions. The quantum conductance for BYP resulted 40 times larger than that of ODT, while the force quantum was 0.8 ± 0.2 nN, considerably smaller than the 1.5 ± 0.2 nN determined for ODT. This last value is the same as that required to break a Au – Au bond (Rubio-Bollinger et al., 2001). Thus, the authors concluded that the breakdown of the ODT nanojunction involves a Au – Au bond rupture, whereas in the case of BYP, the lower breaking force would indicate that a Au – N bond is breaking. These results are in agreement with the notion of the Au- S bond being stronger than the Au-N bond (Stolberg et al., 1990). In a subsequent study, Huang et. al. found that the behavior of the ODT nanojunction as a function of the stretching rate is essentially identical to that of a pure gold point contacts (Huang et al., 2007a).
One useful experimental parameter to determine how strong is the molecule bonded to the metallic electrodes is the length that the junction formed by a single molecule can be stretched before it breaks. This allowed Kiguchi et. al. to establish the following order in binding energies for 1,4-disubstituted benzenes with Au and Pt electrodes: Au-NH2< Pt-NH2 ~ Au-S < Au-isoCN < Pt-isoCN ~ Pt-S (Kiguchi et al., 2006; Kiguchi et al., 2007; Kiguchi et al., 2008).
The statistical analysis of the stretching length was also used to establish the following order in the sense of increasing binding strength: Au-COOH < Au-NH2< Au-S (Chen et al., 2006).
3.2. Computational simulations
As mentioned above, a large proportion of the theoretical work on molecular nanowires involves the calculation of the conductance. Only a few of these have addressed some aspects of the thermodynamic stability of such nanocontacts.
As it respects to the mechanical properties of a molecular nanojunction, one or the first ab-initio studies was performed by Krüger et. al. (Krüger et al., 2003). A thiomethyl radical bonded to a 5-atom planar cluster was used as a model for a typical Au – S contact. The junction was elongated, with geometry optimization at each step, obtaining an isomerization of the cluster into a linear chain which finally breaks at a Au – Au bond with a force of ~ 1.5 nN. In a Car-Parrinello molecular dynamics the same authors showed that when ethylthiol attached to a gold surface is pulled, this leads to the formation of a monoatomic gold nanowire, followed by breaking a Au – Au bond with a rupture force of about 1.2 nN (Krüger et al., 2002).
In a different study, the stretching and breaking behavior of a benzene dithiol molecule sandwiched between two Au(111) slabs was studied using DFT calculations (Lorenz et al., 2006). It was found that breakage occurs through a dissociation of one of the Au – S bonds with a maximum force of 1.25 nN in the case when the molecule is directly attached to the surface, and of 1.9 nN when an adatom is placed between the sulphur and the gold slab.
4. Long term stability of metallic monatomic nanowires
As can be gathered from the summary of the computational work performed so far with the aim of obtaining valuable information about the stability of NWs, the main challenge remains to develop models that would allow to use this valuable computational information to extrapolate to the experimental time scale, and taking thermal motion into account. In the following sections a simple kinetic model based on the Transition State Theory (TST) is presented. The utility of this model along with DFT calculations and an exploration of the energy landscape with a suitable algorithm will be illustrated for pure and contaminated gold monatomic NWs, as well as for single molecule contacts.
4.1. The minimum energy path and the transition state theory
A common and important problem in theoretical chemistry and solid state physics is to identify the path with the lowest energy for the reorganization of a group of atoms from a stable configuration to another. This path is referred to as the Minimum Energy Path (MEP) and is often used to define a Reaction Coordinate (RC) for transitions of the type of chemical reactions, conformational changes in molecules or diffusion processes in solids. The maximum potential energy along the MEP is referred to as the saddle point, and provides the activation energy for the occurrence of the process. To calculate the transition rate constant is of central importance in the TST (Eyring, 1935; Vineyard, 1957), as will be discussed in section 4.2.
Different methods have been developed to find the reaction path and saddle points (Michael & Michael, 2007). We focus our attention on methods that make use of two boundary conditions: the initial and final configurations for the transition. These settings should normally correspond to two local minima in the multidimensional potential energy surface. These minima may be obtained from different energy minimization techniques such as the simulated annealing, conjugate gradient, etc.
These methods require only the calculation of first derivatives of the potential energy. They generate a chain of images or replicas of the system between the initial and final configurations. All the intermediate images are simultaneously optimized in some concerted way of the potential energy surface that should be as close as possible to guarantee the convergence to the MEP. The method called Nudged Elastic Band (NEB) (Henkelman et al., 2000; Henkelman & Jonsson, 2000; Mills & Jónsson, 1994) works in the scheme of these methods and its implementation is particularly simple. The NEB method has been successfully applied to a variety of problems, such as studies of diffusion on metal surfaces (Villarba & Jónsson, 1994), the dissociative adsorption of molecules on a surface (Mills & Jónsson, 1994), and the formation of a contact between a STM tip and a surface (Sørensen et al., 1996).
4.2. Kinetic model
As it was shown by Krüger et. al. in the computational simulation studies described above (Krüger et al., 2002; Krüger et al., 2003), the creation of a Au nanowire takes place in a number of steps, involving an isomerization process. Our approach will only deal with the final stage, where the NW breaks but the model could be extended to a multi-step process. We will assume that, for a given elongation of the wire, it may exist either in a broken (b) or unbroken (u) state (see Fig. 4). The related system energies are denoted with Eb and Eu, respectively. When shifting from state u to b, the system will find an energy barrier (activation energy) that we will denote with ΔE‡. We will neglect the reverse process in all the treatment we give below. In principle, wire reformation could be easily introduced in the model. However, we must take into account that in order to surmount the activation energies involved, the system must gain a considerable amount of energy, which will be released in the downhill stage after crossing through the maximum. This excess energy will rapidly take the systems to other more compact configurations of the final state. To be more illustrative, immediately after the rupture the wire (state b) will find itself in a situation where the Au atoms are in a very low coordination, so that they will stabilize by merging to
Figure 4.
Scheme of the different stages of the rupture of a Au NW. The states labelled as u, b and # correspond to the unbroken, broken and activated states, respectively. The state labelled as c indicates a situation where the broken wire has evolved towards a more compact state, in which the pieces of the broken NW are attached to the bulky pieces of the junction. Although the present figure is meant with illustrative purposes, the configurations and energies were obtained from calculations using the modified embedded atom method potentials (Baskes et al., 1989) and the nudged elastic method (Henkelman et al., 2000; Henkelman & Jonsson, 2000). Reprinted from Vélez, P.; Dassie, S. A.; Leiva, E. P. M. Chem. Phys. Lett., Vol. 460, No.1-3, 261-265, 2008 with permission from Elsevier.
a bulky piece of metal of one of the tips making the junction (state c). These qualitative assessments are based on semiempirical calculations
P. Vélez, S. A. Dassie, E. P. M. Leiva unpublished results.
.
Following the TST, we have for the forward process (wire breaking) the frequency νf given by Eq. 2.
νf=νuexp(−ΔE‡kBT)E2
where kB is Boltzmann’s constant, T the absolute temperature, νu is a factor showing a weak dependence on temperature that represents an average oscillation frequency and ΔE‡ is calculated as E# - Eu. Since ΔE‡ and νu are in principle functions of the wire elongation, this will also be the case of νf. Typical values of νu are 3.5 – 7.0 x 1012 Hz (Bürki et al., 2005; Todorov et al., 2001). We have found that the assumption of constant (elongation independent) values in this range leads essentially to the same qualitative and quantitative predictions we show below. However, in order to get a parameter free description of the problem, it would be desirable to get not only ΔE‡ but also νu from first-principles calculations. Sánchez-Portal et al. (Sánchez-Portal et al., 1999) have calculated the transversal and longitudinal phonon frequencies of a Au NW by means of the frozen phonon method at different wire lengths. In order to get an estimation of νu for our problem, we have parameterized the results of these authors for the transversal mode as a function of the wire elongation and introduced it in our equations for νu (Δz).
Let us now consider a differential elongation of the wire dΔz performed in a period of time dt. The number of crossings from state u to state b over the barrier in dt will then be given by Eq. 3.
dn=[νuexp(−ΔE‡kBT)]dtE3
dn being the average number of times that the system moves from the unbroken to the broken state in the time dt. If we want to calculate the number of possible crossings in a finite period of time we integrate between 0 and t, which is given Eq. 4.
Δn=∫t=0t=τ*[νuexp(−ΔE‡kBT)]dtE4
We can estimate the lifetime of the wire τ* by setting in Δn = 1 (Eq. 5).
∫t=0t=τ*[νuexp(−ΔE‡kBT)]dt=1E5
Note that in the present formulation no assumption has been made on how the wire is elongated. We turn now to consider two different possibilities:
Static rupture of the nanowire: In this case, rupture of the NW is studied at a constant elongation Δz. Under these conditions, the argument of the integral in Eq. 5 becomes independent of time and the breaking time can be straightforwardly obtained by solving Eq. 6, given below.
[νuexp(−ΔE‡kBT)]τ*=1E6
Rupture of the nanowire at a constant elongation rate: If we assume that the NW is stretched at a constant rate, dΔz/dt=νe, then we have Eq. 7, where Δz0 is the elongation at t=0, Δz* is the breaking elongation for the NW and we have written the activation energy as ΔE‡(Δz) to emphasize the dependence of the quantity on the elongation. νu was also considered to be elongation dependent and calculated as pointed out above. However, the dependence of the predictions on this parameter is rather weak. Once the dependence of ΔE‡(Δz) on Δz is given, Eq. 7 can be solved numerically to get the dependence of Δz* on the elongation rate νe.
In order to illustrate the method, we consider a system consisting of a Au nanowire made of a supercell containing 4 atoms which are periodically repeated in space (Au4 NW). Fig. 5 shows a scheme of the unit cell employed to simulate the Au4 NW considered here. We used 4 atoms because the elongation distances at which a NW of this size breaks are in the range between 0.11 to 0.14 nm, which is the value that Huang et al. (Huang et al., 2007a) have found experimentally (see Fig. 9a below). Although this is a rather small system, the rupture of the wire has been found to be the displacement of one of atoms perpendicular to the wire axis (Ke et al., 2007), so that the interaction of the breaking bond with the rest of the system should be minimal. Atom 1 is fixed and the length of the supercell is stretched. For a given stretching of the NW, the energy of the system can be minimized with respect to all the atomic coordinates. Let us denote with Emin (Δz) the minimum energy for a given Δz = L – L0, where L0 is the chain length at equilibrium for each system and L is the total chain length; at some Δz, say Δzmin, Emin will present a minimum, that we denote with Emin (Δz). In the following discussion, we refer all the stretching lengths to Δzmin and all the energies to Emin (Δz).
While the unbroken state is clearly defined, some uncertainty remains concerning the broken one. With this purpose, ab initio molecular dynamics simulations were performed near the point where the NW breaks, finding structures which are very similar to that proposed below in the broken state of Fig. 7, which was obtained by minimizing the energy of the systems, as were the remaining configurations. It is interesting to point out that the broken configuration is the same as one of the most stable Au4 clusters reported in (Bonacic-Koutecky et al., 2002).
Figure 5.
Schematic representation of the unit cell employed in the present calculations. L is the total chain length and 4d¯Au-Au is the average Au-Au separation distance. Reprinted from Vélez, P.; Dassie, S. A.; Leiva, E. P. M. Chem. Phys. Lett., Vol. 460, No.1-3, 261-265, 2008 with permission from Elsevier.
This configuration was then adopted to obtain the broken state for different elongations. This was achieved by compressing the system to a cell size corresponding to the desired length L, and then performing a conjugate gradient minimization to obtain the broken state for each system. The minimum energy path between the unbroken and broken states at each Δz, was then obtained by means of the NEB algorithm.
Curves for the energy of the unbroken state Eu, the broken state Eb, and the activated state E# as a function of the wire elongation are shown in Fig. 6a. The energy curves of the broken and unbroken state meet for an elongation of Δzu-b = 0.1565 nm. At this point the activation energies for the backward and forward reactions are the same and they are equal to 0.33 eV.
Figure 6.
System energy and force curves upon system elongation for the Au4 nanowire. a) Energy of the unbroken state, Eu, energy of the broken state, Eb, and energy of the activated state, E#, as a function of the stretching of the NW, Δz. In all cases the energy of the system is referred to the minimum value that the system takes for the set of Δz considered. The activation energy for NW rupture is ΔE‡ = E# - Eu. b) Force calculated according to Eq. 8 as a function of Δz. Reprinted from Vélez, P.; Dassie, S. A.; Leiva, E. P. M. Chem. Phys. Lett., Vol. 460, No.1-3, 261-265, 2008 with permission from Elsevier.
This would yield a switching frequency between the broken and unbroken state of the order of 106 Hz at room temperature, going down to the order of 1 Hz at 135 K. However, as we stated previously, after the bond is broken the system will immediately rearrange to other more stable configuration so that the forward reaction can be considered as irreversible under the usual stretching conditions. On the other hand, the energies of the unbroken and activated state meet at an elongation Δzu-# = 0.1995 nm. This represents a non-activated rupture process for the NW, and thus the Au – Au average distance at this point (ca. 2.88 Å) represents an absolute limit for wire stability.
Fig. 6b shows the force of the system, calculated as the component of the force applied on the NW at each step from Eq. 8.
FZ=−∂∂ΔzEu(Δz)E8
In Fig. 7 we can observe curves for the energies of the system between broken and unbroken states at some sample elongations. From these curves it can be appreciated that the transition state is shifted towards the broken state. The configurations of the system at the broken, unbroken, and activated state are also shown. It is remarkable that the structure of the broken state is similar to the Au fragment attached to an ethylthiolate molecule, as it was found in an AIMD simulation (Krüger et al., 2002).
Figure 7.
Energy of the system at the reaction path between broken and unbroken states at some sample elongations. The energy is in all cases referred to the unbroken state. The configurations of the system at the broken, unbroken, and activated states are depicted within each frame. Reprinted from Vélez, P.; Dassie, S. A.; Leiva, E. P. M. Chem. Phys. Lett., Vol. 460, No.1-3, 261-265, 2008 with permission from Elsevier.
Static rupture of the nanowire: A logarithmic plot for the static rupture of the NW calculated according to Eq. 6 is shown in Fig. 8. The lowest stretching rates employed in the experiments allow the rupture of NWs in times of the order of 0.1 s. Fig. 8 shows that in this order of times the wires should break at distances close to 0.105 nm. This is very close to the experimental value of 0.1 nm (Huang et al., 2007a) obtained at the lowest stretching rates. Thus, the present results also support the general idea that long Au – Au distances such as those found in the experiments of (Legoas et al., 2002; Ohnishi et al., 1998; Rodrigues & Ugarte, 2001a; Rodrigues & Ugarte, 2001b) at room temperature cannot occur for pure Au NWs. On the other hand, at 150 K and below, pure Au NWs can be considerably stretched beyond that point.
Figure 8.
Decimal logarithm of the lifetime τ* of a Au4 nanowire at different elongations calculated according to Eq. 6 for different temperatures. The dotted line denotes a constant lifetime of 0.1 s. Reprinted from Vélez, P.; Dassie, S. A.; Leiva, E. P. M. Chem. Phys. Lett., Vol. 460, No.1-3, 261-265, 2008 with permission from Elsevier.
Figure 9.
Dynamic rupture of a Au4 nanowire: a) Breaking distance as a function of the logarithm of the elongation rate νe. The full line shows the calculations according to the present model and the squares are data for gold point contacts taken from (Huang et al., 2007a). The segments show the width of the distribution reported there. No fitting parameters were considered. b) Logarithm of the lifetime of the wire, τ*, as a function of the decimal logarithm of the elongation rate, νe. The full line shows the calculations according to the present model and the squares are experimental data for the rupture of the Au – ODT taken from (Huang et al., 2007a). The segments show the width of the distribution reported there. No fitting parameters were considered. Reprinted from Vélez, P.; Dassie, S. A.; Leiva, E. P. M. Chem. Phys. Lett., Vol. 460, No.1-3, 261-265, 2008 with permission from Elsevier.
Rupture of the nanowire at a constant elongation rate: The dynamic rupture of a Au4 nanowire was studied according to Eq. 7. Fig. 9a shows the breaking distance as a function of ln (νe) for elongation rates between e-1 and e6 in comparison with the experimental date taken from Huang et al. (Huang et al., 2007a). It can be observed that the calculated results resemble the experimental trend in the general features.
A further point that can be analyzed through the present calculations is the lifetime of the wires as a function of the elongation rates. These results are given in Fig. 9b. There is no experimental data available for a straightforward comparison with our calculations. However, AIMD simulations (Krüger et al., 2002), and considerations based on experiments (Huang et al., 2006; Huang et al., 2007a; Li et al., 2006a) indicate that the rupture of a nanocontact made of an alkanethiol and Au contacts should break at a Au – Au bond, so that comparison between the present results may be made with the experiment of Huang et al. (Huang et al., 2007a), who have studied the rupture of single molecule junctions involving Au contacts and ODT. The experimental data of Huang et al. is included in Fig. 9b, where it is found that the calculated lifetimes closer resemble those from the experiment; specially taking into account that no fitting attempt was made seeking for agreement.
4.2. Results for contaminated metallic nanowires
Fig. 10 shows a scheme of the unit cell employed to simulate the pure (Au-p) and contaminated Au NWs considered here. The light weight elements taken into account as contaminants are a H atom (Au-H) and a C atom (Au-C). The grey circles represent Au atoms, which remain fixed at their positions during the relaxation processes. The latter consists in a local energy minimization procedure by means of the conjugate gradient method or a search of a MEP by means of the NEB method.
In the case of contaminated NWs the circle marked with an X represents the location the contaminant atom. This figure also shows the definition of the α1, α2 and α3 bond angles, determined by the atoms relevant for the analysis of the rupture of the NW.
Figure 10.
Schematic representation of the unit cell employed in the present calculations: a) the rectangle indicates the extension of the unit cell. The grey circles represent the Au atoms fixed during the simulation. α1, α2 and α3 are the bond angles between the Au atoms, and αX is the angle defined by the atoms Au1 – X – Au2. L is the total chain length and d¯Au-Au is the average Au – Au separation distance. b) Front view of the pyramid of Au atoms shown in a).
The length of the NW, L, is defined here as the distance between the atoms Au1 and Au5. We also define an average Au – Au separation, asd¯Au-Au=L/4. Fig. 10b shows the front view of one of the pyramid of Au atoms, represented as grey circles in Fig. 10a. The present systems contains as a whole 18 atoms. It must be emphasized that in the case of calculations for pure Au NWs, the present systems deliver results that are close to those presented in the previous section for only four atoms in the unit cell. This indicates that the mechanical properties of the NW appear to be quite local, with a rather slight dependence from the bulky atomic environment.
The atomic impurity was located between the atoms Au1 and Au2. This choice was made because in the literature we found first principles calculations, similar to those performed here, where the H and C were positioned at a similar place, as well as between the Au2 and Au3 atoms, with similar results (Novaes et al., 2003; Skorodumova et al., 2007). As it will be found below, the present results agree with those where the impurity was located at another sites of the chain.
Structure and energetics of pure and contaminated Au NWs. For all the systems we shall refer to the equilibrium state as that where the derivative of the energy with respect to the elongation is equal to zero (ie. the external force, or stress, acting on the system is null). The rupture force will be considered to be the value of the force Fz where it presents a maximum at long elongations, being the force defined in Eq. 8 (da Silva et al., 2004; Jelínek et al., 2008; Novaes et al., 2003; Novaes et al., 2006; Rubio-Bollinger et al., 2001; Vélez et al., 2005; Vélez et al., 2008). Accordingly, we refer in the following to the “at rupture” state as that where the coordinates of the atoms are such that ∂Fz/∂Δz=0 (maximum force). With this definition, we are trying to address the status of the system just at the point where the NW is breaking by further force application. The structural information at the equilibrium and “at rupture” states of the NWs is summarized in Table 1 for the three systems types considered.
Considering the length difference between the “at rupture” and equilibrium states of the different systems, it is found that the stretching lengths of the Au-p and Au-C systems are 1.9 Å and 1.6 Å respectively, while the elongation of the Au-H system is considerably larger (2.6 Å). This fact bears direct consequences for the force constant kz, as will be seen later on. Inspecting Table 1 it becomes clear that the equilibrium geometries of the contaminated
Equilibrium geometries
„At rupture“ geometries
Au-p
Au-H
Au-C
Au-p
Au-H
Au-C
Le or L* / Å
9.80
9.60
10.8
11.7
12.2
12.4
d (1,2) / Å
2.63
2.76
3.78
2.85
3.58
3.87
d (2,3) / Å
2.64
2.62
2.59
2.95
2.75
2.67
d (3,4) / Å
2.63
2.67
2.66
3.05
3.16
3.15
d (4,5) / Å
2.63
2.61
2.62
2.85
2.71
2.71
α1 / °
139.5
149.6
161.8
179.3
178.6
178.5
α2 / °
125.2
108.6
116.5
179.7
177.0
176.5
α3 / °
139.7
157.4
136.0
179.9
178.5
178.2
αX / °
---
102.5
168.5
---
179.2
179.5
Table 1.
Structural information of the pure Au (Au-p) and hydrogen (Au-H) or carbon contaminated (Au-C) NWs at the equilibrium and close to rupture situations. In the first row, Le is the equilibrium length of the NW, whereas L* corresponds to length just before breakage occurs. Refer to Fig. 11 for an illustration of each of the geometries analyzed. The d (1,2) bond distance and the αX bond angle of the impurity are in boldface.
Figure 11.
a) Equilibrium geometries for the pure gold NW (Au-p) and contaminated with H (Au-H) and C (Au-C). b) Geometries just before rupture for stretched NWs. The arrows indicate the bond at which the NW breaks.
systems are very different from each other. Close to equilibrium, the presence of the H atom changes only slightly the d (1,2) bond distance, while the carbon atom inserts itself almost completely in the middle of this bond. This fact can be appreciated in the characteristics of the equilibrium configuration (see αX bond angle and d (1,2) bond distance).
At the situation of the rupture, all three systems present a linear configuration and the distance of the Au3 – Au4 bond where it occurs is very similar (~3.1 Å). The contaminated systems arrive to the point of rupture with large Au – Au distances at the Au1 – Au2 bond; 3.58 Å for the Au-H system and 3.87 Å for the Au-C system. These figures indicate that our results are in a very good agreement with the Au – Au separation distances found experimentally (Kizuka, 2008; Kondo & Takayanagi, 2000; Legoas et al., 2002; Ohnishi et al., 1998; Rodrigues et al., 2000; Yanson et al., 1998) and with the first-principles calculations from other research groups (Novaes et al., 2003; Skorodumova et al., 2007).
Up to now, we have only explored the possibility of the incorporation of a single atomic impurity. In relation to this, it must be recognized that the extremely large Au – Au distances of 4 – 5 Å (Legoas et al., 2002), as well as consecutive distances of 3.5 – 4 Å (Kondo & Takayanagi, 2000; Ohnishi et al., 1998; Rodrigues & Ugarte, 2001a; Rodrigues et al., 2000) cannot be reproduced by our calculations. To tackle this point, more impurities should be considered, and probably other molecular species as proposed by other groups that performed first-principles calculations (Bahn et al., 2002; Galvão et al., 2004; Legoas et al., 2002; Novaes et al., 2003; Skorodumova & Simak, 2003; Skorodumova & Simak, 2004). However, the present approach is not devoted to predict all the Au – Au distances observed, but to understand the effect that an atomic impurity produces on a monatomic Au NW from energetic, geometrical, electronic and kinetic viewpoints.
Analysis of the time-stability of pure and contaminated NWs. In the following the kinetic aspects of the rupture process for the Au-p, Au-H and Au-C systems are considered. Fig. 12a shows the energy curves for the unbroken (u), broken (b) and activated (#) states as a function of the elongation (Δz) for Au-H. For Au-p the respective curves are similar to those shown in Fig. 6 for the 4-atom gold nanowire. The results obtained for Au-C have qualitatively the same features as those obtained for Au-p.
Figure 12.
a) Energy values corresponding to the unbroken state Eu, broken state Eb and activated state E# of the system Au-H as a function of the elongation. b) Continuous curve: z component of the force acting on the system, Fz, calculated according to Eq. 8, as a function of Δz for Au-H. Broken line: Linear fit to Fz between zero and the maximum force. The maximum force Fzmax is indicated in normal typeset, while the force constant, kz (slope of the broken lines), is denoted in italics.
In order to find candidate configurations for the broken (final in TST jargon) state at each elongation, molecular dynamics runs were performed at an elongation corresponding to the “at-rupture” configuration. This led to structures that looked like those illustrated for Au-H for the broken state configurations in Fig. 14 below.
We will perform a more detailed discussion of the reaction path for the rupture of the NWs below. We first consider the behaviour of the forces along the stretching procedure, since they are closely related to the energy curves of the unbroken state. In fact, the energy curves for the unbroken states in Fig. 12a may be used to calculate the longitudinal force Fz acting on this system according to Eq. 8. In Fig. 6, it was found that the stability limit of a monoatomic Au NW is reached at a force of 1.45 nN, which is close to the experimental value of 1.5 ± 0.3 nN (Kizuka, 2008; Rubio et al., 1996; Rubio-Bollinger et al., 2001). Fig. 12b shows Fz for the Au-H as a function of the elongation Δz. The maximum force Fzmax is also reported there. Table 2 compiles relevant information for this and the other systems under consideration in this study (Au-p, Au-C). For the Au-p system, the maximum force is 1.57 nN, also in agreement with the experiments and other theoretical values (da Silva et al., 2004; Rubio-Bollinger et al., 2001; Vélez et al., 2008). In the case of the contaminated systems, this value is somewhat lower (1.19 nN for Au-H and 1.14 nN for Au-C). First principles results from literature show the same trend (Novaes et al., 2003; Skorodumova et al., 2007).
Fzmax, nN
kz, N/m
ΔE‡(Δz), eV(Å)
Au-p
1.57
9.1
0.99 (0.0)
0.22 (1.8)
Au-H
1.19
5.0
1.33 (0.0)
1.02 (2.7)
Au-C
1.14
8.2
1.12 (0.0)
1.10 (1.8)
Table 2.
Maximum tensile force sustained before breakage (Fzmax, nN), force constant (kz, N/m) and activation energy (ΔE‡, eV) for the minimum energy path at the equilibrium and “at-rupture” geometries, for the systems Au-p, Au-H and Au-C. The activation energy for the rupture of the NW at each Δz is calculated as ΔE‡ = E# - Eu.
In their studies on the mechanic properties of monatomic Au NWs, (Rubio-Bollinger et al., 2001) found that these chains are five times harder than the massive electrodes. They evaluated from the experimental results the slopes of the force curves in the last stage of the elastic deformation, before the rupture of the NW, from a set of 200 experiments. The average value of the force constant for an average chain length was 8 N/m. The present calculations show that the behaviour of the force curves of the Au-p, Au-H and Au-C systems is quite elastic, as can be inferred from the linear fit of Fz between zero and the maximum force, shown in Fig. 12b for Au-H as a broken line. In the elastic deformation region, the force Fz may be written as:
Fz=kzΔzE9
where kz is the force constant of the system. The value of kz found from the present calculations for Au-p, 9.1 N/m, is in perfect agreement with the experimental value reported in references (Rubio-Bollinger et al., 2001; Xu et al., 2003a). The contaminated systems present a remarkably different behaviour when compared between each other. The Au-C system has kz =8.1 N/m, a very similar value to that of Au-p, but 1.6 larger than the value of this property for Au-H (kz=5.0 N/m). These are interesting predictions, since up to date no measurements have been performed comparing pure and contaminated systems. A suitable fitting of the experimental data of this property could help to shed light on the type of impurity present when a NW exhibits large Au – Au separation distances.
Fig. 13 shows the energy of the system, E, as a function of the normalized reaction coordinate (NRC) for the system Au-H for some sample elongations. The corresponding configurations are presented in Fig. 14.
Fig. 13 shows that the Au-H system presents an energy curve along the minimum energy path that is more complex than those for Au-p and Au-C (not shown here). The distinctive feature for of Au-H is that at short elongations the energy curves present a minimum. This behavior is related to the incorporation of the hydrogen atom into the wire that takes place as the elongation proceeds. However, for the Au-C system, the C atom is incorporated into the wire since the beginning of the elongation. The minimum in the minimum energy path profile of the Au-H system disappears gradually with increasing elongation of the NW, remaining a shoulder at long elongations (see Fig. 13b).
The images of monatomic Au chains showing long separation distances (3.6 – 4 Å) were obtained by TEM and HRTEM using electron beam lithography for the fabrication of the NWs (Kizuka, 2008; Kondo & Takayanagi, 2000; Legoas et al., 2002; Ohnishi et al., 1998; Rodrigues & Ugarte, 2001a; Rodrigues et al., 2000). Under these experimental conditions, the elongation rate is not controlled and in principle not known. However, an estimation can be made looking at the pictures provided in some of these publications, as the one illustrated in
Figure 13.
System energy, E, as a function of the normalized reaction coordinate (NRC) along the reaction path between unbroken (u) and broken (b) states for different elongations Δz of the NW. a), b) and c) correspond to the Au-p system, d), e) and f) to the Au-H system and g), h) and i) to the Au-C system. In all cases the energy of the system is referred to the energy of the unbroken state (u). The activation energies ΔE‡ and the corresponding elongations are given for each case.
Figure 14.
Configurations of the unbroken (u), activated (#) and broken (b) states of the NW for each of the elongations Δz whose energy curves along the reaction path are shown in Fig. 13 for Au-H.
Fig. 2 (Rodrigues & Ugarte, 2001a), and is found to be very slow, of the order of 0.1 nm/s in the final stage, allowing for a complete equilibration of the system at the atomic scale all along the elongation process. Thus, the rupture of the NWs in these experiments takes place in the timescale of a second, so that it can be inferred that the activation barriers must be quite high. In the static limit, the lifetime τ* of a NW can be estimated from Eq 6. In this equation, we have used νu(Δz) = νu = 3 x 1012 Hz and T = 300 K. Similar values were used in references (Bürki et al., 2005; Rubio et al., 1996; Todorov et al., 2001; Vélez et al., 2008). Fig. 15 shows a plot of the decimal logarithm of the lifetime τ* of the NW, as a function of the elongation force Fz for the Au-p, Au-H and Au-C systems.
It becomes evident from this figure that the contaminated NWs live considerably longer than Au-p NWs for all elongations. It is found that pure Au NW become unstable in the experimental time scale from an elongation which corresponds to dAu-Au = 2.74 Å / atom on. In fact, for larger Fz (dAu-Au), pure NWs should live less than 0.1 s. The remarkable behaviour of the Au-H system is due to the shape of the activation energy, discussed above in Fig. 13. Fig. 15 also supports the idea that impurities are responsible for the large Au – Au distances observed experimentally. These impurities, as stated above, modify geometrically the structure of the NWs, giving place to the occurrence of the anomalous large Au – Au distances. Furthermore, they modify the NW chemically by charge transfer and rearrangement, which in turn changes the potential energy surfaces so that for each elongation the contaminated NWs present higher activation barriers for the rupture than the pure NWs. Fig. 15 makes also plausible that stretched (contaminated) NWs may have lifetimes of the order of the second, sometimes even of the order of minutes (Ohnishi et al., 1998).
Figure 15.
Decimal logarithm of the lifetimes of the NWs corresponding to the systems Au-p (), Au-H (---) and Au-C (••••) as a function of the elongation force Fz, calculated according to Eq. 6.
5. Long term stability of molecular nanowires. Effect of substituents on the Au – N nanojunction.
The formalism proposed in the previous sections could, in principle, be applied to study the mechanical and kinetic stability of single molecule nanowires. To the best of our knowledge, no such study has been yet attempted. One possible approximation to the long term stability problem is to estimate the kinetic barrier for the rupture of the nanowire by means of the binding energy of the system (Eb).
The mechanical properties of ortho substituted pyrazines bonded to planar gold clusters of 6 and 7 atoms have been studied using the same methodology described for the monatomic metallic nanowires. In Table 3, we summarize the values obtained for the binding energies (Eb) and rupture forces (Frup) for the substituted pyrazines bonded to the 6-atom (Au6) and the 7-atom (Au6+1) gold clusters.
The systems studied are illustrated in Fig. 16. These systems are taken as a model for the last stage in the stretching of a Au – bipyridine nanojunction, and allow an assessment of the effect of a substituent in the molecule on the stability of the nanojunction.
The Eb values summarized in Table 3 were obtained as the energy difference between the minimum energy structure for each system and the final broken structure after stretching the nanojunction. In all cases the rupture occurred at the Au – N bond. Frup values were obtained using Eq. 8.
Figure 16.
Illustration of the models adopted to study the effect of ortho substituents on the mechanical stability of a Au – N nanojunction. a) Au6 + X-pyrazine. b) Au6+1 + X-pyrizane.
Au6 + X-pyrazine
Au6+1 + X-pyrazine
X
σ
Eb, eV
Frup, nN
Eb, eV
Frup, nN
NH2
-0.66
0.795
1.16
0.913
1.26
CH3
-0.17
0.731
1.07
0.828
1.13
H
0.00
0.671
1.00
0.778
1.09
Cl
0.23
0.596
0.82
0.689
0.98
CONH2
0.36
0.549
0.77
---
---
CN
0.66
0.497
0.77
0.605
0.88
NO2
0.78
0.409
0.48
0.560
0.73
Table 3.
Hammett’s σ (values taken from (Hansch et al., 1991)), binding energies (Eb) and rupture forces (Frup) obtained from stretching experiments for Au6 + X-pyrazines and Au6+1 X-pyrazines.
Using the Eb values as a lower limit for the activation energy to break the nanojunction, and using Eq. 6 with νu = 1012 Hz (a typical value for these systems as shown above), time constant values, τ*, can be calculated for the different nanojunctions. The values for τ so obtained represent a lower limit for the lifetime of the nanojunctions. Interestingly, there is a an excellent correlation between the logarithmic lifetime and a pure empirical parameter describing the electronic nature of the substituent such as Hammett’s σ (Hansch et al., 1991), as illustrated in Fig. 17a.
Fig. 17b shows that a good correlation is also observed with respect to the calculated rupture force for the nanojunction. This is interesting since Frup values can be obtained experimentally.
Note that the electronic properties of the substituent has a marked effect on the lifetime of the nanojunction, which ranges from 10-6 seconds for X = NO2 to about 103 seconds for NH2. On the other hand, only a slight effect of the substituent was found on the experimentally measured conductance of a structurally similar system (Venkataraman et al., 2007).
When the molecule binds to a less coordinated gold atom, such as in Au6+1 – X-pyrazine systems, the strength of the nanojunction increases, as reflected by the higher values obtained for both Eb and Frup. This also causes an increase in the lifetime of the molecular junction, as illustrated in Fig. 17.
Figure 17.
Correlation between log τ, calculated as the lower limit for the nanojunction’s lifetime and a) Hammett’s σ and b) calculated rupture force (Frup). In both plots, open circles represent the Au6 + X-pyrazine systems, whereas filled triangles are for the Au6+1 – X-pyrazine systems.
A comprehensive revision of the most recent advances with respect to the experimental techniques and computational simulations focused on the stability and mechanical properties of monatomic metallic and single molecule nanowires is presented.
We have established a simple model based on the Transition State Theory and using the Minimum Energy Path to study the long-term stability of nanowires in a time scale corresponding to that of the experimental observations. The utility of this model has been demonstrated for the quantitative evaluation of the effect of impurities in gold nanowires, as well as for a qualitative assessment of the effect exerted by a change in the electronic properties of a molecule in the temporal stability of a metal – molecule junction.
Some of the future work within this topic includes the analysis of other systems, including different metals and / or other chemical linkers for molecular nanowires. It would also be of great interest to evaluate, if there is any, the correlation between the stability of the nanowire and its electronic properties, such as conductance.
References
1.AmorimE. P. M.daSilva. A. J. R.daSilva. E. Z.2008 Computer Simulations of Copper and Gold Nanowires and Single-wall Nanowires, J. Phys. Chem. C, 11239 15241-15246, 1932-7447
2.AmorimE. P. M.daSilva. A. J. R.FazzioA.daSilva. E. Z.2007 Short linear atomic chains in copper nanowires, Nanotechnology, 1814 145701, 0957-4484
3.AngladaE.TorresJ. A.YndurainF.SolerJ. M.2007 Formation of Gold Nanowires with Impurities: A First-Principles Molecular Dynamics Simulation, Phys. Rev. Lett., 989 096102-096104, 0031-9007
4.AulettaT.de JongM. R.MulderA.van VeggelF. C. J. M.HuskensJ.ReinhoudtD. N.ZouS.ZapotocznyS.SchönherrH.VancsoG. J.KuipersL.2004 β-Cyclodextrin Host−Guest Complexes Probed under Thermodynamic Equilibrium: Thermodynamics and AFM Force Spectroscopy, J. Am. Chem. Soc., 1265 1577-1584, 0002-7863
6.BahnS. R.JacobsenK. W.2001 Chain Formation of Metal Atoms, Phys. Rev. Lett., 8726 266101-1-266101-4, 0031-9007
7.BaskesM. I.NelsonJ. S.WrightA. F.1989 Semiempirical modified embedded-atom potentials for silicon and germanium, Phys. Rev. B, 409 6085-6100, 1098-0121
8.BatistaR. J. C.OrdejonP.ChachamH.ArtachoE.2007 Resistive and rectifying effects of pulling gold atoms at thiol-gold nanocontacts, Phys. Rev. B, 754 041402-041404, 1098-0121
9.BeebeJ. M.EngelkesV. B.MillerL. L.FrisbieC. D.2002 Contact Resistance in Metal-Molecule-Metal Junctions Based on Aliphatic SAMs: Effects of Surface Linker and Metal Work Function, J. Am. Chem. Soc., 12438 11268-11269, 0002-7863
10.BettiniJ.SatoF.CouraP. Z.DantasS. O.GalvãoD. S.UgarteD.2006 Experimental realization of suspended atomic chains composed of different atomic species, Nat Nano, 13 182-185, 1748-3387
11.Bonacic-KouteckyV.BurdaJ.MitricR.GeM.ZampellaG.FantucciP.2002 Density functional study of structural and electronic properties of bimetallic silver--gold clusters: Comparison with pure gold and silver clusters, J. Chem. Phys., 1177 3120-3131, 0021-9606
12.BürkiJ.StaffordC. A.SteinD. L.2005 Theory of Metastability in Simple Metal Nanowires, Phys. Rev. Lett., 959 090601, 0031-9007
13.ChenF.LiX.HihathJ.HuangZ.TaoN.2006 Effect of Anchoring Groups on Single-Molecule Conductance: Comparative Study of Thiol-, Amine-, and Carboxylic-Acid-Terminated Molecules, J. Am. Chem. Soc., 12849 15874-15881, 0002-7863
15.daSilva. E. Z.daSilva. A. J. R.FazzioA.2001 How Do Gold Nanowires Break?, Phys. Rev. Lett., 8725 256102, 0031-9007
16.daSilva. E. Z.NovaesF. D.daSilva. A. J. R.FazzioA.2004 Theoretical study of the formation, evolution, and breaking of gold nanowires, Phys. Rev. B, 6911 115411 (111 ), 1098-0121
17.De MariaL.SpringborgM.2000 Electronic structure and dimerization of a single monatomic gold wire, Chemical Physics Letters, 3233-4 , 293-299, 0009-2614
18.DreherM.PaulyF.HeurichJ.CuevasJ. C.ScheerE.NielabaP.2005 Structure and conductance histogram of atomic-sized Au contacts, Phys. Rev. B, 727 075435-11, 1098-0121
20.EvansE.1999 Introductory Lecture Energy landscapes of biomolecular adhesion and receptor anchoring at interfaces explored with dynamic force spectroscopy, Faraday Discussions, 111 1-16, 1364-5498
21.EvansE.2001 Probing the Relation Between Force-Lifetime and Chemistry in Single Molecular Bonds, Annual Review of Biophysics and Biomolecular Structure, 301 105-128, 1056-8700
22.EyringH.1935 The Activated Complex in Chemical Reactions, J. Chem. Phys., 32 107-115, 0021-9606
23.GalvãoD. S.RodriguesV.UgarteD.LegoasS. B.2004 The role of carbon contamination in metallic nanowires, Materials Research, 7 339-342, 1516-1439
24.HaissW.MartínS.LearyE.ZalingeH. v.HigginsS. J.BouffierL.NicholsR. J.2009 Impact of Junction Formation Method and Surface Roughness on Single Molecule Conductance, The Journal of Physical Chemistry C, 11314 5823-5833, 1932-7447
25.HaissW.NicholsR. J.van ZalingeH.HigginsS. J.BethelD.SchiffrinD. J.2004 Measurement of single molecule conductivity using the spontaneous formation of molecular wires, Phys. Chem. Chem. Phys., 6 4330-4337, 1463-9076
26.HaissW.WangC.GraceI.BatsanovA. S.SchiffrinD. J.HigginsS. J.BryceM. R.LambertC. J.NicholsR. J.2006 Precision control of single-molecule electrical junctions, Nature Materials, 512 995-1002, 1476-1122
27.HaissW.WangC.JitchatiR.GraceI.MartínS.BatsanovA. S.HigginsS.BryceM. R.LambertC. J.JensenP. S.NicholsR. J.2008 Variable contact gap single-molecule conductance determination for a series of conjugated molecular bridges, Journal of Physics: Condensed Matter, 2037 374119 (9pp), 0953-8984
28.HäkkinenH.BarnettR. N.ScherbakovA. G.LandmanU.2000 Nanowire Gold Chains: Formation Mechanisms and Conductance, J. Phys. Chem. B, 10439 9063-9066, 1520-6106
29.HanschC.LeoA.TaftR. W.1991 A Survey of Hammett Substituent Constants and Resonance and Field Parameters, Chem. Rev., 912 165-195, 0009-2665
30.HeH. X.ShuC.LiC. Z.TaoN. J.2002 Adsorbate effect on the mechanical stability of atomically thin metallic wires, Journal of Electroanalytical Chemistry, 5221 26-32, 0022-0728
31.HenkelmanG.JonssonH.2000 Improved tangent estimate in the nudged elastic band method for finding minimum energy paths and saddle points, J. Chem. Phys., 11322 9978-9985, 0021-9606
32.HenkelmanG.UberuagaB. P.JonssonH.2000 A climbing image nudged elastic band method for finding saddle points and minimum energy paths, J. Chem. Phys., 11322 9901-9904, 0021-9606
33.HobiE.daSilva. A. J. R.NovaesF. D.daSilva. E. Z.FazzioA.2005 Comment on "Contaminants in Suspended Gold Chains: An Ab Initio Molecular Dynamics Study", Phys. Rev. Lett., 9516 169601 (1 , 0031-9007
34.HobiE.Jr FazzioA.daSilva. A. J. R.2008 Temperature and Quantum Effects in the Stability of Pure and Doped Gold Nanowires, Phys. Rev. Lett., 1005 056104 (14 ), 0031-9007
35.HouS.NingJ.ShenZ.ZhaoX.XueZ.2006 Influences of the molecule-electrode interface structure on the conducting characteristics of the gold-4,4 bipyridine-gold molecular junction, Chem. Phys., 3271 1-9, 0301-0104
36.HouS.ZhangJ.LiR.NingJ.HanR.ShenZ.ZhaoX.XueZ.WuQ.2005 First-principles calculation of the conductance of a single 4,4 bipyridine molecule, Nanotechnology, 162 239-244, 0957-4484
37.HuangZ. F.ChenF.BennettP. A.TaoN. J.2007a Single Molecule Junctions Formed via Au-Thiol Contact: Stability and Breakdown Mechanism, J. Am. Chem. Soc., 12943 13225-13231, 0002-7863
38.HuangZ. F.XuB. Q.ChenY. C.VentraM. D.TaoN. J.2006 Measurement of Current-Induced Local Heating in a Single Molecule Junction, Nano Lett., 66 1240-1244, 1530-6984
39.HuangZ.ChenF.D’agostaR.BennettP. A.Di VentraM.TaoN.2007b Local ionic and electron heating in single-molecule junctions, Nature Nanotechnology, 211 698-703, 1748-3387
41.JelínekP.PérezR.OrtegaJ.FloresF.2008 Ab initio study of evolution of mechanical and transport properties of clean and contaminated Au nanowires along the deformation path, Phys. Rev. B, 7711 115447 (112 ), 1098-0121
42.KamenetskaM.KoentoppM.WhalleyA. C.ParkY. S.SteigerwaldM. L.NuckollsC.HybertsenM. S.VenkataramanL.2009 Formation and Evolution of Single-Molecule Junctions, Phys. Rev. Lett., 10212 126803-126804, 0031-9007
43.KeL.KotaniT.van SchilfgaardeM.BennettP. A.2007 Breakdown of a gold nanowire between electrodes, Nanotechnology, 1842 424002, 0957-4484
44.KiguchiM.MiuraS.HaraK.SawamuraM.MurakoshiK.2006 Conductance of a single molecule anchored by an isocyanide substituent to gold electrodes, Appl. Phys. Lett., 8921 213104-3, 0003-6951
45.KiguchiM.MiuraS.HaraK.SawamuraM.MurakoshiK.2007 Conductance of single 1,4-disubstituted benzene molecules anchored to Pt electrodes, Appl. Phys. Lett., 915 053110-053113, 0003-6951
46.KiguchiM.MiuraS.TakahashiT.HaraK.SawamuraM.MurakoshiK.2008 Conductance of Single 1,4-Benzenediamine Molecule Bridging between Au and Pt Electrodes, J. Phys. Chem. C, 11235 13349-13352, 1932-7447
47.KiguchiM.MurakoshiK.2006 Fabrication of stable Pd nanowire assisted by hydrogen in solution, Appl. Phys. Lett., 8825 253112-253113, 0003-6951
48.KimB. S.BeebeJ. M.JunY.ZhuX. Y.FrisbieC. D.2006a Correlation between HOMO Alignment and Contact Resistance in Molecular Junctions: Aromatic Thiols versus Aromatic Isocyanides, J. Am. Chem. Soc., 12815 4970-4971, 0002-7863
49.KimY. H.Tahir-KheliJ.SchultzP. A.GoddardW. A.I. I. I.2006b First-principles approach to the charge-transport characteristics of monolayer molecular-electronics devices: Application to hexanedithiolate devices, Phys. Rev. B, 7323 235419-10, 1098-0121
50.KizukaT.2008 Atomic configuration and mechanical and electrical properties of stable gold wires of single-atom width, Phys. Rev. B, 7715 155401 (111 ), 1098-0121
51.KoizumiH.OshimaY.KondoY.TakayanagiK.2001 Quantitative high-resolution microscopy on a suspended chain of gold atoms, Ultramicroscopy, 881 17-24, 0304-3991
53.KondoY.TakayanagiK.2000 Synthesis and Characterization of Helical Multi-Shell Gold Nanowires, Science, 2895479 606-608, 0036-8075
54.KrügerD.FuchsH.RousseauR.MarxD.ParrinelloM.2002 Pulling Monatomic Gold Wires with Single Molecules: An Ab Initio Simulation, Phys. Rev. Lett., 8918 186402, 0031-9007
55.KrügerD.RousseauR.FuchsH.MarxD.2003 Towards "Mechanochemistry": Mechanically Induced Isomerizations of Thiolate-Gold Clusters, Angew. Chem. Int. Ed. Eng., 4220 2251-2253, 1433-7851
56.LegoasS. B.GalvaoD. S.RodriguesV.UgarteD.2002 Origin of Anomalously Long Interatomic Distances in Suspended Gold Chains, Phys. Rev. Lett., 887 076105 (14 ), 0031-9007
57.LegoasS. B.RodriguesV.UgarteD.GalvãoD. S.2004 Contaminants in Suspended Gold Chains: An Ab Initio Molecular Dynamics Study, Phys. Rev. Lett., 9321 216103 (14 ), 0031-9007
58.Liet.alC. Z.1999 Fabrication of stable metallic nanowires with quantized conductance, Nanotechnology, 102 221, 0957-4484
59.LiC. Z.TaoN. J.1998 Quantum transport in metallic nanowires fabricated by electrochemical deposition/dissolution, Appl. Phys. Lett., 728 894-896, 0003-6951
60.LiC.PobelovI.WandlowskiT.BagretsA.ArnoldA.EversF.2007b Charge Transport in Single Au | Alkanedithiol | Au Junctions: Coordination Geometries and Conformational Degrees of Freedom, J. Am. Chem. Soc., 1301 318-326, 0002-7863
61.LiC.PobelovI.WandlowskiT.BagretsA.ArnoldA.EversF.2007a Charge Transport in Single Au | Alkanedithiol | Au Junctions: Coordination Geometries and Conformational Degrees of Freedom, J. Am. Chem. Soc., 1301 318-326, 0002-7863
62.LiQ.WuX.HuangJ.YangJ.2005 Electronic transport property of 4,4’-bipyridine molecular junction, Ultramicroscopy, 1051-4 , 293-298, 0304-3991
63.LiX.HeJ.HihathJ.XuB.LindsayS. M.TaoN.2006a Conductance of Single Alkanedithiols: Conduction Mechanism and Effect of Molecule-Electrode Contacts, J. Am. Chem. Soc., 1286 2135-2141, 0002-7863
64.LiZ. L.ZouB.WangC. K.LuoY.2006b Effects of Electrode Distances on Geometric Structure and Electronic Transport Properties of Molecular 4,4’-Bipyridine Junction, Journal of Physics: Conference Series, 29 110-114, 1742-6596
65.LiZ. L.ZouB.WangC. K.LuoY.2006c Electronic transport properties of molecular bipyridine junctions: Effects of isomer and contact structures, Phys. Rev. B, 737 075326-075327, 1098-0121
66.LorenzR.GeorgH.MathisG.Jean-LucB.EgbertZ.2006 Stretching and Breaking of a Molecular Junction, Small, 212 1468-1475, 1613-6829
67.MartínS.HaissW.HigginsS.CeaP.LopezM. C.NicholsR. J.2008 A Comprehensive Study of the Single Molecule Conductance of α,ω-Dicarboxylic Acid-Terminated Alkanes, J. Phys. Chem. C, 11210 3941-3948, 1932-7447
68.MerkelR.NassoyP.LeungA.RitchieK.EvansE.1999 Energy landscapes of receptor-ligand bonds explored with dynamic force spectroscopy, Nature, 3976714 50-53, 0028-0836
69.MichaelL. M.MichaelP.2007 Computing Reaction Pathways on Molecular Potential Energy Surfaces, In: Reviews in Computational Chemistry 4 Kenny, B. L. (Ed.), 35-65,
70.MillsG.JónssonH.1994 Quantum and thermal effects in H2 dissociative adsorption: Evaluation of free energy barriers in multidimensional quantum systems, Phys. Rev. Lett., 727 1124, 0031-9007
71.NakamuraA.BrandbygeM.HansenL. B.JacobsenK. W.1999 Density Functional Simulation of a Breaking Nanowire, Phys. Rev. Lett., 827 1538, 0031-9007
72.NakamuraJ.KobayashiN.AonoM.2001 Electronic states and structural stability of gold nanowires, RIKEN Review, 200137 17-20, 0919-3405
73.NovaesF. D.daSilva. A. J. R.daSilva. E. Z.FazzioA.2003 Effect of Impurities in the Large Au-Au Distances in Gold Nanowires, Phys. Rev. Lett., 903 036101 (14 ), 0031-9007
74.NovaesF. D.daSilva. A. J. R.daSilva. E. Z.FazzioA.2006 Oxygen Clamps in Gold Nanowires, Phys. Rev. Lett., 961 016104 (14 ), 0031-9007
75.O’HanlonJ. F.2001A User’s Guide to Vacuum Technology, John Wiley & Sons, Inc, 978-0-47127-052-2 Hoboken, NJ.
76.OhnishiH.KondoY.TakayanagiK.1998 Quantized conductance through individual rows of suspended gold atoms, Nature, 3956704 780-783, 0028-0836
77.OkamotoM.TakayanagiK.1999 Structure and conductance of a gold atomic chain, Phys. Rev. B, 6011 7808-7811, 1098-0121
78.ParkY. S.WhalleyA. C.KamenetskaM.SteigerwaldM. L.HybertsenM. S.NuckollsC.VenkataramanL.2007 Contact Chemistry and Single-Molecule Conductance: A Comparison of Phosphines, Methyl Sulfides, and Amines, J. Am. Chem. Soc., 12951 15768-15769, 0002-7863
79.ParkY. S.WidawskyJ. R.KamenetskaM.SteigerwaldM. L.HybertsenM. S.NuckollsC.VenkataramanL.2009 Frustrated Rotations in Single-Molecule Junctions, J. Am. Chem. Soc., 13131 10820-10821, 0002-7863
80.PatroneL.PalacinS.BourgoinJ. P.2003a Direct comparison of the electronic coupling efficiency of sulfur and selenium alligator clips for molecules adsorbed onto gold electrodes, Applied Surface Science, 212-213 , 446-451, 0169-4332
81.PatroneL.PalacinS.CharlierJ.ArmandF.BourgoinJ. P.TangH.GauthierS.2003b Evidence of the Key Role of Metal-Molecule Bonding in Metal-Molecule-Metal Transport Experiments, Phys. Rev. Lett., 919 096802, 0031-9007
82.PaulssonM.KragC.FrederiksenT.BrandbygeM.2008 Conductance of Alkanedithiol Single-Molecule Junctions: A Molecular Dynamics Study, Nano Lett., 91 117-121, 1530-6984
83.Perez-JimenezA.2005 Uncovering Transport Properties of 4,4’-Bipyridine/Gold Molecular Nanobridges, J. Phys. Chem. B, 10920 10052-10060, 1520-6106
84.PuQ.LengY.TsetserisL.ParkH. S.PantelidesS. T.CummingsP. T.2007a Molecular dynamics simulations of stretched gold nanowires: The relative utility of different semiempirical potentials, J. Chem. Phys., 12614 144707-6, 0021-9606
85.PuQ.LengY.ZhaoX.CummingsP. T.2007b Molecular simulations of stretching gold nanowires in solvents, Nanotechnology, 1842 424007, 0957-4484
88.QuinnJ. R.FossF. W.VenkataramanL.HybertsenM. S.BreslowR.2007 Single-Molecule Junction Conductance through Diaminoacenes, J. Am. Chem. Soc., 12921 6714-6715, 0002-7863
89.RodriguesV.BettiniJ.RochaA. R.RegoL. G. C.UgarteD.2002 Quantum conductance in silver nanowires: Correlation between atomic structure and transport properties, Phys. Rev. B, 6515 153402, 1098-0121
90.RodriguesV.UgarteD.2001a Real-time imaging of atomistic process in one-atom-thick metal junctions, Phys. Rev. B, 637 073405 (14 ), 1098-0121
91.RodriguesV.UgarteD.2001b Structural and electronic properties of gold nanowires, Eur. Phys. J. D, 161-3 , 395-398, 1434-6060
92.RodriguesV.UgarteD.2002 Quantum conductance properties of metal nanowires, Materials Science and Engineering B, 962 188-192, 0921-5107
93.RodriguesV.BettiniJ.SilvaP. C.UgarteD.2003 Evidence for Spontaneous Spin-Polarized Transport in Magnetic Nanowires, Phys. Rev. Lett., 919 096801, 0031-9007
94.RodriguesV.FuhrerT.UgarteD.2000 Signature of Atomic Structure in the Quantum Conductance of Gold Nanowires, Phys. Rev. Lett., 8519 4124-4127, 0031-9007
96.Rubio-BollingerG.BahnS. R.AgraïtN.JacobsenK. W.VieiraS.2001 Mechanical Properties and Formation Mechanisms of a Wire of Single Gold Atoms, Phys. Rev. Lett., 872 026101 (14 ), 0031-9007
97.Sánchez-PortalD.ArtachoE.JunqueraJ.OrdejónP.GarcíaA.SolerJ. M.1999 Stiff Monatomic Gold Wires with a Spinning Zigzag Geometry, Phys. Rev. Lett., 8319 3884-3887, 0031-9007
98.SatoF.MoreiraA. S.BettiniJ.CouraP. Z.DantasS. O.UgarteD.GalvaoD. S.2006a On the Formation of Copper Linear Atomic Suspended Chains, Condensed Matter, arXiv:cond-mat/06020921
99.SatoF.MoreiraA. S.BettiniJ.CouraP. Z.DantasS. O.UgarteD.GalvaoD. S.2006b Transmission electron microscopy and molecular dynamics study of the formation of suspended copper linear atomic chains, Phys. Rev. B, 7419 193401-193404, 1098-0121
100.SchönherrH.BeulenM. W. J.BüglerJ.HuskensJ.van VeggelF. C. J. M.ReinhoudtD. N.VancsoG. J.2000 Individual Supramolecular Host−Guest Interactions Studied by Dynamic Single Molecule Force Spectroscopy, J. Am. Chem. Soc., 12220 4963-4967, 0002-7863
101.SkorodumovaN. V.SimakS. I.2003 Stability of gold nanowires at large Au-Au separations, Phys. Rev. B, 6712 121404 (14 ), 1098-0121
102.SkorodumovaN. V.SimakS. I.2004 Stabilization of monoatomic gold wires by carbon impurities, Solid State Communications, 13011 755-757, 0038-1098
103.SkorodumovaN. V.SimakS. I.KochetovA. E.JohanssonB.2007 Ab initio study of electronic and structural properties of gold nanowires with light-element impurities, Phys. Rev. B, 7523 235440 (14 ), 1098-0121
104.SørensenM. R.BrandbygeM.JacobsenK. W.1998 Mechanical deformation of atomic-scale metallic contacts: Structure and mechanisms, Phys. Rev. B, 576 3283, 1098-0121
105.SørensenM. R.JacobsenK. W.JónssonH.1996 Thermal Diffusion Processes in Metal-Tip-Surface Interactions: Contact Formation and Adatom Mobility, Phys. Rev. Lett., 7725 5067, 0031-9007
106.StadlerR.ThygesenK. S.JacobsenK. W.2005 Forces and conductances in a single-molecule bipyridine junction, Phys. Rev. B, 7224 241401-241404, 1098-0121
107.StolbergL.LipkowskiJ.IrishD. E.1990 Adsorption of pyridine at the Au(110)- solution interface, Journal of Electroanalytical Chemistry, 2961 171-189, 0022-0728
108.TakaiY.KawasakiT.KimuraY.IkutaT.ShimizuR.2001 Dynamic Observation of an Atom-Sized Gold Wire by Phase Electron Microscopy, Phys. Rev. Lett., 8710 106105, 0031-9007
110.TorresJ. A.TosattiE.DalCorso. A.ErcolessiF.KohanoffJ. J.Di TollaF. D.SolerJ. M.1999 The puzzling stability of monatomic gold wires, Surface Science, 4263 L441-L446, 0039-6028
111.UlrichJ.EsrailD.PontiusW.VenkataramanL.MillarD.DoerrerL. H.2006 Variability of Conductance in Molecular Junctions, J. Phys. Chem. B, 1106 2462-2466, 1520-6106
112.van RuitenbeekJ.ScheerE.WeberH.2005 Contacting Individual Molecules Using Mechanically Controllable Break Junctions, In: Introducing Molecular Electronics, Cuniberti, G., Fagas, G., and Richter, K. (Ed.), 253274 , Springer Berlin Heidelberg, 16166361
113.VélezP.DassieS. A.LeivaE. P. M.2005 First Principles Calculations of Mechanical Properties of 4,4’-Bipyridine Attached to Au Nanowires, Phys. Rev. Lett., 954 045503-045504, 0031-9007
114.VélezP.DassieS. A.LeivaE. P. M.2008 When do nanowires break? A model for the theoretical study of the long-term stability of monoatomic nanowires, Chem. Phys. Lett., 4601-3 , 261-265, 0009-2614
115.VélezP.DassieS. A.LeivaE. P. M.2007 Understanding the properties of nanocontacts. Recent advances in experiment and theory, In: Recent Advances in Nanoscience, Mariscal, M. M. and Dassie, S. A. (Ed.), 138 , Research Signpost, 978-8-13080-207-7 Kerala, India.
117.VenkataramanL.ParkY. S.WhalleyA. C.NuckollsC.HybertsenM. S.SteigerwaldM. L.2007 Electronics and Chemistry: Varying Single-Molecule Junction Conductance Using Chemical Substituents, Nano Lett., 72 502-506, 1530-6984
118.VenkataramanL.KlareJ. E.NuckollsC.HybertsenM. S.SteigerwaldM. L.2006b Dependence of single-molecule junction conductance on molecular conformation, Nature, 4427105 904-907, 0028-0836
119.VillarbaM.JónssonH.1994 Diffusion mechanisms relevant to metal crystal growth: Pt/Pt(111), Surface Science, 3171-2 , 15-36, 0039-6028
120.VineyardG. H.1957 Frequency factors and isotope effects in solid state rate processes, Journal of Physics and Chemistry of Solids, 31-2 , 121-127, 0022-3697
121.WuX.LiQ.HuangJ.YangJ.2005 Nonequilibrium electronic transport of 4,4’-bipyridine molecular junction, J. Chem. Phys., 12318 184712-184716, ISSN:
122.XuB.XiaoX.TaoN. J.2003a Measurements of Single-Molecule Electromechanical Properties, J. Am. Chem. Soc., 12552 16164-16165, 0002-7863
123.XuB. Q.LiX. L.XiaoX. Y.SakaguchiH.TaoN. J.2005 Electromechanical and Conductance Switching Properties of Single Oligothiophene Molecules, Nano Lett., 57 1491-1495, 1530-6984
124.XuB.HeH.BoussaadS.TaoN. J.2003b Electrochemical properties of atomic-scale metal wires, Electrochimica Acta, 4820-22 , 3085-3091, 0013-4686
125.XuB.TaoN. J.2003 Measurement of Single-Molecule Resistance by Repeated Formation of Molecular Junctions, Science, 3015637 1221-1223, 0021-9606
126.YansonA. I.BollingerG. R.van denBrom. H. E.AgraïtN.van RuitenbeekJ. M.1998 Formation and manipulation of a metallic wire of single gold atoms, Nature, 3956704 783-785, 0028-0836
127.YasudaS.YoshidaS.SasakiJ.OkutsuY.NakamuraT.TaninakaA.TakeuchiO.ShigekawaH.2006 Bond Fluctuation of S/Se Anchoring Observed in Single-Molecule Conductance Measurements using the Point Contact Method with Scanning Tunneling Microscopy, J. Am. Chem. Soc., 12824 7746-7747, 0002-7863
128.YoshidaM.KuruiY.OshimaY.TakayanagiK.2007 In situ Observation of the Electrical Conductance of the Bridged Single Fullerene using TEM-STM, Microscopy and Microanalysis, 13 No.Supplement S02, 748-749, 1431-9276
129.ZhangC.BarnettR. N.LandmanU.2008 Bonding, Conductance, and Magnetization of Oxygenated Au Nanowires, Phys. Rev. Lett., 1004 046801-046804, 0031-9007
130.ZoloffMichoff. M. E.VélezP.LeivaE. P. M.2009 Substituent Effect on the Mechanical Properties of Au−N Nanojunctions, J. Phys. Chem. C, 1139 3850-3854, 1932-7447
131.ZouS.SchonherrH.VancsoG. J.2005 Force Spectroscopy of Quadruple H-Bonded Dimers by AFM: Dynamic Bond Rupture and Molecular Time−Temperature Superposition, J. Am. Chem. Soc., 12732 11230-11231, 0002-7863
Notes
Times Online September 19, 2007, http://technology.timesonline.co.uk/tol/news/tech_and_web/article2489053.ece
P. Vélez, S. A. Dassie, E. P. M. Leiva unpublished results.
Written By
Martin E. Zoloff Michoff, Patricio Vélez, Sergio A. Dassie and Ezequiel P. M. Leiva