Martin
Schnedlitz
a,
Maximilian
Lasserus
a,
Daniel
Knez
b,
Andreas W.
Hauser
*a,
Ferdinand
Hofer
b and
Wolfgang E.
Ernst
*a
aInstitute of Experimental Physics, Graz University of Technology, Petersgasse 16, A-8010 Graz, Austria. E-mail: andreas.w.hauser@gmail.com; wolfgang.ernst@tugraz.at; Fax: +43 (316) 873-108140; Tel: +43 (316) 873-8157 Tel: +43 (316) 873-8140
bInstitute for Electron Microscopy and Nanoanalysis & Graz Centre for Electron Microscopy, Graz University of Technology, Steyrergasse 17, A-8010 Graz, Austria
First published on 10th March 2017
We present time-resolved transmission electron microscopy studies of the degradation of Au, Ag, Cu and Ni nanowires deposited on a heated support. The wires are grown under fully inert conditions in superfluid helium droplets and deposited onto amorphous carbon. The inherent stability of these pristine metal nanowires with diameters below 10 nm is investigated in the absence of any stabilizers, templates or solvents. The phenomenon of Rayleigh-breakup, a consequence of diffusion processes along the wire surfaces, is analysed in situ via scans over time and support temperature. Our experimental efforts are combined with simulations based on a novel model featuring a cellular automaton to emulate surface diffusion. Based on this model, correlations between the material parameters and actual breakup behaviour are studied.
In a recent publication we studied the Rayleigh breakup of Ag nanowires,14 exploiting superfluid He nanodroplets (HeN) as nanolabs for the production of one-dimensional metallic structures.15,16 The doping of HeN with particles is well established in spectroscopy,17–19 and has recently been adapted for the controlled production and structure preserving soft deposition of metal clusters.20–24 For the creation of metallic nanowires we take advantage of a vortex formation, a quantum effect25,26 occurring for larger droplets with diameters of a few μm. As was originally observed in bulk superfluid helium (He II),27–29 these vortices attract dopants immersed in the helium due to a pressure gradient around the vortex core,30,31 causing a preference for one-dimensional growth which leads to the formation of nanowires.32–35
The experimentally challenging technique offers the benefit of complete inertness during wire production, deposition and measurement. Typically, nanowires in the sub 10 nm diameter range are stabilized by templates,36–38 which introduce a large bias in stability measurements.
When comparing the structures obtained for different metals after deposition on amorphous carbon via transmission electron microscopy (TEM), it became obvious that silver behaved significantly differently from the other elements. Whereas continuous wires could be observed for Au, Cu and Ni under standard conditions, the typical outcome for silver deposition would be an accumulation of segmented chains of silver clusters. This phenomenon was originally thought to be a consequence of the formation within the He droplets before deposition. However, a detailed study of the wire structures as a function of the substrate temperature revealed that this segmentation was caused by Rayleigh breakup, a temperature-dependent outcome of surface diffusion processes along the wires.14 This behaviour is well captured by an adaptation of Rayleigh's original model for the breakup of a fluid jet39 performed by Nichols and Mullins.40–42 While this has been postulated earlier by Gordon et al.,33,43 it could only be proven by a series of TEM images taken at increasing substrate temperatures. The study revealed that silver nanowires are particularly sensitive to heating with breakup temperatures as low as 263 K.
The current article aims to deepen our understanding of this diffusion process with regard to material-specific constants. In order to do this, we perform a combined approach involving theory and experiment.
On the experimental side, we broaden our spectrum of material investigation by detailed, temperature-programmed diffusion studies, originally only performed for Ag, now also for the elements Au, Ni and Cu. These metals have been chosen in order to cover a wide stability range. In addition to previous studies, the current experimental setup combines a very efficient and reproducible method for the production of metallic nanowires with a novel method for highly resolved in situ measurements of surface diffusion processes in ultrathin nanowires. This unique combination of apparatus and analysis tools offers the first in situ observations of the very nature of the formation and surface diffusion of metallic clusters.
On the theoretical side, we present a novel approach of describing the experimentally observed Rayleigh breakup phenomenon. Following ref. 44, we previously applied a two-dimensional computational method based on the Mullins equation,45 which states a relation between the flux of surface atoms and the gradient of the local surface curvature. With this method we could predict the location and the order (with respect to time evolution) of future breakups for initial structures of the arbitrary shape.14 Despite its success, this approach suffers from severe limitations, both in terms of computational efficiency as well as in its ability to grasp the real-world process of material redistribution: the actual motion of atoms in the diffusion process is represented by an artificial force acting on the contour line of the metallic structure. Therefore, we replaced the contour-line method by a cellular automaton (CA) ansatz on a three-dimensional grid. This new method, which is discussed in greater detail in the theory section of this article, together with additional data from the experiment, allows us to address several new questions regarding surface diffusion on metallic nanostructures:
(1) The higher computational efficiency allows the simulation of larger nanostructures and takes their actual height on the substrate into consideration.
(2) Material properties such as bond lengths and the binding energy enter the simulations. This way, metal-specific behaviour can be simulated and concrete changes of the structure can be reproduced in the time- and temperature domain.
(3) The broader choice of metals allows for a first analysis of trends in the periodic system and of tentative correlations between the breakup behaviour and the related material properties.
Our article is structured as follows: first, we give a brief overview of our experimental setup in Section 2. We discuss the He-droplet-mediated creation of metallic nanowires and their follow-up measurement. In Section 3 we provide details of the computational approach. A comparison of the simulated results and experimental nature is given in Section 4. Correlations between the actual behaviour of various metals and their corresponding atomic properties or material constants are analysed.
To investigate the temperature dependent diffusion processes for room-temperature-stable metal nanowires, scanning transmission electron microscope (STEM) high-angle annular dark-field (HAADF) images are recorded using a probe-corrected FEI Titan3 G2 60-300. The microscope is equipped with a Gatan Quantum energy filter for electron energy loss spectroscopy (EELS) and a four-quadrant energy-dispersive X-ray spectroscopy (EDX) detector (FEI Super-X), for the determination of material compositions. The temperature of the heatable carbon grid can be controlled from room temperature up to ≈1300 °C with heating rates of 200 °C ms−1. Since some nanowires show breakups below ambient temperature conditions14 the second substrate holder is cooled to liquid nitrogen (N2,liq) temperature. After the deposition of nanowires the cooled substrates are detached from the sample holder and immersed in N2,liq within ≈5 s. During this time the sample holder is flooded with dry nitrogen to prevent contact with ambient air. Afterwards the sample is analysed in a N2,liq environment via a FEI Tecnai 12 TEM, equipped with a Bio Scan Camera (Gatan, Model 792), where it is possible to manipulate the substrate temperature from 77 K to 363 K via the cryo-TEM holder (Gatan model No. 613).
The advantage over previous measurements is the possibility of in situ recording of breakup events as a function of the substrate temperature. To observe temperature dependent breakup events, recordings are taken during the heat-up process using a small temperature ramp of 2 K min−1 and at constant substrate temperatures. For each of the Au, Ag, Ni and Cu wires, one particular segment of the observable area was kept constant for a continuous observation throughout the whole heating time.
In the literature, experiments were described on nanowires made from e.g. Ag,47 Cu,48,49 Au,50,51 Sn52 and Pt,53 all with diameters above 20 nm. Breakup temperatures well below the bulk melting temperature were reported in all cases, but still far above room temperature. Only for Ag nanofractals, with initial branch diameters of ≈16 nm, structural changes near room temperature were observed.54,55 We note that a breakup of Ag nanowires, due to surface effects, has also been observed by Gordon et al. in bulk He II.32,43
Our novel ansatz is fundamentally different from this previous approach and much closer to the actual physical process of particle migration. We employ a cellular automaton (CA), which uses the information obtained from the TEM data to create a three dimensional reconstruction of the nanowires in the form of a three-dimensional logical grid.
Given a change in the detected intensity with the material thickness, and assuming a circular cross section for each nanostructure, we can translate image contrast data into a three-dimensional object. As TEM intensity information is provided in the form of a two-dimensional matrix, the wires are created layer by layer, starting with a lower threshold I0 for the central layer, and counting only those pixels with intensities above I0. For the next layers this threshold is continuously increased, which adds additional voxels above the starting layer and results in a three-dimensional reconstruction of the total structure.
Each element of the grid represents a ‘voxel’, whose size is given by the pixel size of the TEM image. Following the basic idea of a CA, such a voxel is either filled with material or remains empty. For our data set, one voxel measures (0.95 nm)3, which corresponds to approximately 40 atoms per filled voxel.
The dynamics of the CA is defined by the choice of update rules which apply to all cells or voxels of the system. These updates are executed simultaneously for all cells in each time step.56 A simulation of the thermally induced surface diffusion process can now be obtained by choosing update rules which describe a particle hopping on the grid. This process depends on the local environment of a voxel, but must also account for the stochastic character of a temperature-driven particle migration. This can be achieved by employing an ansatz from reaction kinematics, the Eyring equation,57 defining jumping rates via
(1) |
As k is strongly dependent on the temperature, the reaction rate will drastically increase as soon as the thermal energy is high enough to overcome the energy barrier ΔE and detach one voxel from its neighbour cells. After detaching, the voxel is free to change its position within its nearest neighbourhood. The number of nearest and next nearest neighbours determines, together with the probabilities derived for possible migrations, to which position the voxel will most likely be transferred. With the ongoing simulation time, this algorithm leads to preferably spherical structures, where the number of nearest neighbours in total reaches a maximum. Note that this model does not take any substrate interaction into account due to the extremely low wetting behaviour on the chosen substrate. This is particularly obvious for the interaction of gold and carbon in comparison to the intermetallic bond energy (EAu–C = 0.03 eV per atom,59EAu–Au = 3.81 eV per atom58). For a two-dimensional study of surface diffusion including the interaction with various substrates we refer to the work of Solov'yov et al. on Ag nanoparticles.60 Diffusion time is not directly linked to the in situ observed diffusion process. Therefore, a fit between the observed time dependent TEM surface diffusion process at constant temperature and the stepwise iteration process is necessary for realistic predictions of the behaviour over time. The simple linear dependence of the pre-exponential factor in eqn (1) is derived from statistical mechanics under the assumption of a single degree of freedom being responsible for reaction barrier transmission. The fact that not every attempt in this degree of motion leads to a successful hopping is usually taken into consideration by a proportionality constant κ referred to as the transmission coefficient. Although not in the atomic picture, we can follow this nomenclature for our discussion of voxel migration and include all unknown scaling factors which occur in the coarse-grained model into this constant κ. With F = E − TS we can rewrite eqn (1) as
(2) |
(3) |
Fig. 1 displays graphically the situation of a selected surface voxel and its initial neighbourhood. In this setup, the energy can be minimized by moving the green voxel in a way which maximizes its number of neighbour voxels. The stochastic character of the process of breaking and forming bonds enters viaeqn (1), which yields the probabilities for hopping events.
With our technique we were able to perform in situ observations of the very nature of this thermally induced breakup down to wire diameters of 1 nm covering a total observation area of (3 mm)2.
Particularly interesting is the finding that the breakup temperatures are far below the corresponding bulk melting temperatures. Therefore, a special cryo-TEM experiment for silver nanowires had to be designed, which showed that breakup starts at temperatures as low as ≈258 K (−15 °C).14 In the current extended study, all metals are analysed with respect to different diameters of the initially stable nanowires. The measurement confirmed the known relation between the breakup temperature of nanowires and their corresponding bulk melting temperature given by58
(4) |
When heating above the breakup-temperature, a fragmentation takes place as it is predicted by the CA simulations. The resulting nanostructures show elliptical or spherical shapes. Details of this process and a comparison with our simulations are given below.
The following features and their temperature dependence are well predicted by our simulations for all materials tested:
˙ Smoothing of the surface of the wires as a result of the induced surface diffusion.
˙ Rayleigh breakup and fragmentation of the nanowire.
˙ Contraction of nanostructures by minimization of the number of surface atoms
˙ The exact breakup points and the final shape of the wires
Our CA ansatz further allows the prediction of the breakup-sequence of constriction points. We note that the diffusion process is simulated in three dimensions, although measurement is limited to a “projected”, two-dimensional view. The dynamics during the simulations matches the structural changes observed in the TEM nicely. The agreement between simulations and measured structures after heating can be seen in Fig. 4. Even for very complex forms we are able to accurately predict the exact surface diffusion process. However, a correction factor has to be introduced to the Eyring equation in the case of the more reactive metals Cu and Ni, which accounts for the formation of an oxide layer on the surface of those clusters. This oxide layer is assumed to quench particle migration, leading to an inhibited or delayed breakup.
For the case of gold nanowires we also present a more detailed study of time dependence which includes a fit of the transmission coefficient κ(T). Its deviation is based on a series of TEM images taken at a constant substrate temperature of 150 °C. The surface diffusion of gold nanowires in real time is studied and compared to the structure evolution simulated with the CA. A selection of nanoparticles and their corresponding constriction points in time are used to perform a linear fit of the real time evolution as a function of the number of update cycles. For this particular combination of materials (Au nanowires on a non-wetting support) and the given temperature we find that 1500 real time seconds of TEM observation correspond to 100 iteration steps. One example of a Rayleigh breakup as a function of time and iteration step is presented in Fig. 5.
For the non-noble metals like Ni and Cu an oxide layer covering the metal cluster is observed (see Fig. 3) which results in a higher breakup temperature as the oxide has to evaporate before the diffusion process can continue. As a result, larger deviations from the linear behaviour are observed.
Cu shows a particularly interesting behaviour, as the surface diffusion for pure Cu takes place at room temperature. However, the exposure to oxygen at least during the transfer process of the specimen stops the diffusion process, which is continued at much higher temperatures after the oxygen layer is evaporated inside the TEM.
This journal is © the Owner Societies 2017 |