Open Access Article

This Open Access Article is licensed under a

Creative Commons Attribution 3.0 Unported Licence

DOI: 10.1039/C9CE01893J
(Paper)
CrystEngComm, 2020, Advance Article

Vyacheslav N. Gorshkov^{a},
Vladimir V. Tereshchuk^{a} and
Pooya Sareh*^{b}
^{a}National Technical University of Ukraine, Igor Sikorsky Kyiv Polytechnic Institute, 37 Prospect Ceremony, Kiev 03056, Ukraine
^{b}School of Engineering, University of Liverpool, London Campus, 33 Finsbury Square, London EC2A 1AG, UK. E-mail: pooya.sareh@liverpool.ac.uk

Received
29th November 2019
, Accepted 12th February 2020

First published on 25th March 2020

We have analysed the expressed manifestation of the anisotropy of surface energy density in the dynamics of ultrathin nanowires, which break up into disjointed clusters when annealed below their melting temperature. The breakup process is studied for different temperatures and orientations of the nanowire axis relative to its internal crystal structure using the Monte Carlo kinetic method. We have also presented an approximate analytical model of the instability of nanowires. Generally, the interpretation of experimental results refers to the theoretical model developed by Nichols and Mullins, which is based on conceptions about the Rayleigh instability of liquid jets. In both cases, the theories – which do not take into account the anisotropy of surface energy density – predict the breakup of a nanowire/liquid jet with radius r into fragments with an average length Λ = 9r. However, the observed value, Λ/r, often deviates from 9 either to lower values or to substantially greater ones (up to 24–30). Our results explain various observed features of the breakup and the significant variations in the values of its parameter Λ/r depending on experimental conditions. In particular, the ambiguous role of exchange by atoms of the surface of a nanowire with a surrounding layer of free atoms formed as a result of their rather intensive sublimation, which occurs in a number of cases, has been investigated. We have shown that this exchange can lead both to a decrease and to a significant increase in the parameter Λ/r. The obtained results could be potentially useful in applications such as the development of optical waveguides based on ordered nanoparticle chains.

The existing analytical model of the nanowire breakup by Nichols and Mullins^{24} is very close to the model of Plateau–Rayleigh instability in liquid inviscid threads^{25} because it is constructed under the assumption that the surface energy density, σ, of the nanowire is isotropic. Despite the fact that internal flows are absent in the case of nanowire breakup and according to the statement in ref. 24, only the surface diffusion of atoms is responsible for nanowire dynamics, the wavelength of periodic perturbation of its radius with maximal growth increment, λ_{max}, is the same as that in cylindrical liquid streams, i.e. λ_{max} ≈ 9r_{0}, where r_{0} is the initial radius of the nanowire. The origin of this coincidence is that, in both cases, changes in total surface energy and the displacement of center of mass contained in a half-wavelength segment (between necking and widening) for the supposed sinusoidal modulations of radius are similar.

The correspondence of the predictions of the theoretical model^{24} with the experimentally measured average sizes of nanodroplets and the distances between them, Λ, was noted in a number of previous studies performed for Cu, Ag, Au, and Pt-nanowires.^{26–30} The discrepancies in the experimental and theoretical values, which are approximately 10–16% as noted by the authors, are quite expected, taking into account the assumptions of the model itself and the measurement errors in the experiments as well as the role of uncontrolled factors. However, there are also known cases of extremely large values of the breakup parameter, Λ/r_{0} ∼ 25–30.^{31} Moreover, significant periodic modulations of the nanowire radius (Δr_{max} ∼ 0.5r_{0}) can be observed in a high-resolution TEM photo^{32} with a wavelength of λ_{max} ≈ 5.7r_{0}, which is below the classical threshold for the development of instability,^{24} λ_{max} > 2πr_{0}. Therefore, the physical mechanisms responsible for such dispersion in the parameters of this instability require careful future research.

The observed deviations in the average distance between the centers of neighboring nanodroplets, Λ, which was associated with the wavelength λ_{max}, must be interpreted as a result of the anisotropy of surface energy density. The first discussion of this factor^{33,34} resulted in a clear qualitative estimation: if σ decreases/increases at the necking regions relative to the widening regions when the surface perturbations arise, then this decrease/increase leads to acceleration/inhibition in the development of the instability. Less obvious and more complex is the answer to the questions how σ inhomogeneity affects the value Λ and why in some experiments a strong dependency of Λ on temperature, T, is observed.^{26,28} This is the fact that the well-known model^{24} cannot explain, as well as the processes that are responsible for giant wavelength perturbations, Λ ∼ 30r_{0}, as reported in the mentioned studies^{28,31} (see Fig. 1(a)). Finally, we would like to note a series of publications^{35–38} in which the physical processes concerning the nanowire dynamics were studied in the frameworks of stochastic atomistic modeling and molecular dynamics. However, the physical mechanisms determining the experimentally observed diversity of the breakup parameters, Λ and λ_{max}, (in general, the origin of these values is different^{39,40}) are not sufficiently clear in detail.

Fig. 1 (a) Experimental configuration obtained for an Au-nanowire at temperature T = 400 °C adapted from the study of Vigonski et al.^{31} (b) Result of our numerical experiments with a nanowire (L/r = 124 and Λ/r ≈ 25) oriented along the [110]-direction; the initial nanowire shape is also outlined. |

The effects associated with σ(r) anisotropy are expected to manifest when the axis of the nanowire does not coincide with the axis of symmetry of the crystal lattice. This is why these ‘asymmetrical’ orientations were studied in our previous work.^{39} But it turned out, as shown in our present work, that the most striking deviations in the dynamics of the nanowire, in comparison with the model,^{24} are manifested precisely in the axisymmetric orientations [100], [111], and [110] (the nanowire axis is the axis of symmetry of the fourth-, third-, and second-order, respectively). A variety of scenarios in terms of the physical factors determining the breakup parameter, Λ, is based on that in these orientations the nanowire lateral surface is mainly composed of faces of (100), (111), and (110) types with the lowest surface energy densities, σ(r). Therefore, the formation of the nanowire constrictions is accompanied by the appearance of a noticeable imbalance between the values, σ(r), at the narrowing and bulging areas of the nanowire in one direction or another, with the total free-energy decreasing in any case.

It is to be noted that, in addition to the surface diffusion of atoms, detachment–reattachment processes can play an important role in the development of nanowire instabilities. We have established that the role of the exchange by atoms strongly depends on the orientation of a nanowire and can lead either to a significant reduction in the wavelengths of excited surface perturbations or to their considerable increase (up to λ_{max} ≈ (25–30) × r_{0}; see Fig. 1(b) and section 4), as compared with the results obtained in cases in which the detachment processes are blocked. Also, the approximate analytical model of the disintegration process is given taking into account the anisotropy of the surface energy density (see section 3), the results of which predict the effects observed in the experiments.

All the results of our numerical experiments were interpreted on the basis of clear qualitative physical concepts that represent in detail the features of breakup scenarios, which can be realized in the FCC-nanowire dynamics depending on experimental conditions. The exploited Monte Carlo method, operating with two dimensionless parameters, has been successfully used earlier in a series of articles dedicated to modeling the surface- and nanotube-templated growth^{41,42} of nanoparticles and other nanostructures, as well as the sintering of the former.^{43,44} The detailed description of this method can be found in our previous study.^{39} However, its main concepts will be briefly outlined below (see section 2).

The free atoms diffuse by hopping with random angles of scattering and fixed length, , of the hops (, the lattice spacing, a, is taken as the unit of length in sequence). They can be reattached to the nanowire surface (at unoccupied lattice sites that are the nearest-neighbors to the nanostructure atoms) if they hop into a cube with edge length a/2, centered at such a site. The nanowire axis is also the axis of a cylindrical container with walls that reflect the free atoms. The radius of this container is much greater than the nanowire radius.

In order to evaluate the role of detachment–reattachment processes in the nanowire break-up, we also conducted numerical experiments in which the evaporation from the nanowire surface was reduced. For this, an additional coefficient P_{filter} < 1 sometimes was used in the model to decrease the detachment probability of bound atoms by 1/P_{filter} times.

To avoid the shrinkage of the finitely-long nanowire and early generation of nanodroplets at its ends because of the so-called “end-effect”,^{39,40} we consider long nanowires with ‘frozen’ ends. In this case, 5–7 atomic planes of moveless/frozen atoms, which are oriented perpendicular to the nanowire axis, connect the nanowire butt-ends with the end-walls of the container. This approach is an analogue of the periodic boundary condition.

In accordance with prior studies of the morphology dynamics of FCC nanoclusters,^{39–47} the present mesoscale model operates by two dimensionless thermodynamic parameters,

α = |ε|/kT and p = exp(−Δ/kT), | (1) |

(2) |

(3) |

The calculation of ΔE_{surf} is performed for a cylinder with a variable radius along its axis, X, as follows

r(x, t = 0) = r_{0} − (ε^{2}/4r_{0}) − εcos(kx), k = 2π/λ, ε ≪ r_{0},
| (4) |

(5) |

(6) |

(7) |

(8a) |

(8b) |

(9) |

If we take the surface energy density on a facet of (111) type equal to one, i.e. σ_{(111)} = 1, then according to the broken-bond rule for the surface energies of noble metals,^{45} we will have σ_{(100)} = 1.155 and σ_{(110)} = 1.225. The results presented in ref. 31 for Au give slightly different values for the surface energy density, but in all cases the relation σ_{(111)} < σ_{(100)} < σ_{(110)} holds, which allows us to estimate the value of λ_{max}/r_{0} for nanowires with different orientations on the basis of eqn (9). Such an estimation enables one to predict only the deviations of the parameter λ_{max}/r_{0} from the classical value either in a positive or negative direction. As will be shown in the next section, in some cases these deviations can significantly exceed 9 (in agreement with experiments^{28,30,31}), because the presented model does not take into account some important physical factors.

(10) |

We begin our analysis of the numerical results with the most interesting case in which the wire is oriented along the [110] axis. As noted in many previous studies (see, e.g.,^{28,30,31,48}), in this condition the nanowire is the most resistant to breakup and is characterized by an anomalously large value of the parameter Λ/r_{0} ∼ 20–30. The fact that the wavelength of perturbations, λ_{max}, of the nanowire radius should exceed 9r_{0} immediately before the start of breakup follows from eqn (9) and is based on the considerations explained below.

It is known that the surface of a cylindrical nanowire is rearranged by minimizing free energy at the initial stages of evolution. Fig. 2 shows the structure of the lateral surface of some nanowires as a result of such transformation. In the case of the [110] orientation (Fig. 2(a)), the lateral surface of the nanowire is represented mainly by the (111) faces with the minimum surface energy density, σ. On the (100) faces, σ is also relatively low. The arising instability is followed by the formation of necking regions with slopes that are formed by facets with greater values of σ. Thus, the average value of the parameter β in eqn (9) is positive, i.e. β > 0, that corresponds to the excitation of long-wave modulations of the nanowire radius: λ_{max}/r_{0} > 9. The results of numerical experiments (Fig. 3) confirm this preliminary estimate. The data presented in Fig. 3 were obtained by blocking the evaporation of atoms from the surface of a nanowire (P_{filter} = 0). That is, the development of instability is associated only with surface diffusion, which was assumed in the derivation of eqn (9).

Fig. 2 Configurations of the lateral surface of a nanowire when its axis is oriented in the (a) [110]-, (b) [111]-, and (c) [100]-directions (images show fragments of long nanowires). |

One can see that the breakup parameter, Λ/r_{0}, monotonically decreases with increasing temperature. According to the results presented in Fig. 3, the ratio, Λ/r_{0}, is close to 9 at high temperatures (Fig. 3C; Λ/r_{0} ≈ 10) and significantly exceeds this value in the cold regime (Fig. 3A; Λ/r_{0} ≈ 14.5). The physical interpretation of this experimentally observed fact^{26} is straightforward in our numerical model. The jump frequency of bound atoms – locally determined by the value of p_{jump} = p^{m0} = exp(−m_{0} × Δ/kt) (see section 2) – is inhomogeneous along a nanowire surface (p_{jump} is minimal at widening regions). At low temperatures, this inhomogeneity becomes more pronounced. Naturally, the higher surface activity of atoms in necking regions leads to the accumulation of atoms in front of the widening zones and results in increasing the wavelength of surface perturbations.

The results presented in Fig. 3 show that under long-wave perturbations, the number of nanodroplets after the breakup of the nanowire coincides with the number of modulation periods of its cross section, i.e. Λ = λ_{max}. Such a coincidence is violated in the case of the [100]-orientation (see below), when, at the stage of the breakup of the nanowire, the developed shorter wavelength perturbations merge into separate nanoclusters, which results in inequality Λ > λ_{max}.

Note that significant variations in the breakup parameter, Λ/r_{0}, as a function of temperature are observed in previously reported experiments with Cu (ref. 26) when the temperature increases from 400 °C (Λ/r_{0} ≈ 13–15) to 500 °C. In the second case, the diameter of the formed nanodroplets is such that the gap between them is approximately equal to their diameter. For this geometric ratio the breakup parameter is Λ/r_{0} ≈ 9–10. However, in experiments with Au reported in ref. 31 and 28, the observed ratio Λ/r_{0} is up to 20–24 (at T ≈ 500 °C) and demonstrates a record high value Λ/r_{0} ≈ 30–35 at T ≈ 400 °C. We obtain similar results if the evaporation of bound atoms from the nanowire surface is taken into account (P_{filter} = 1) (see Fig. 4). The effect of atom exchange between the nanowire surface and the surface layer of free atoms increases the parameter Λ/r_{eff} from 10 to 16 in the hot regime and from 14.5 to 24–30 in the cold regime. The effective nanowire radius is

(11) |

A significant increase in the Λ/r_{eff} parameter (compare the results of Fig. 3 and 4) and the sensitivity of this parameter to temperature changes are inherent to only the [110]-orientation. The breakup characteristics of a nanowire with a [111]-orientation are completely opposite. If the evaporation of atoms from the surface of the nanowire is blocked (P_{filter} = 0), then in the warm mode Λ/r_{0} ≈ 12.5 (see Fig. 5D). Accounting for detachment–reattachment processes (P_{filter} = 1) does not increase (as it occurs in the [110]-orientation), but reduces the breakup parameter to the classical Λ/r_{eff} ≈ 9 regardless of temperature (see Fig. 5A–C).

The reasons for this striking difference can be explained on the basis of the structural features of the side surface of the nanowire in the cases under consideration (in the [110] and [111]-orientations), as discussed below.

The flux density of atomic evaporation from the nanowire surface, Φ^{(−)}, is a local characteristic that, in a rough approximation, depends on the average surface curvature. The deposition of free atoms (the density of their reverse flow, Φ^{(+)}) is associated with the distribution of the surface curvature along the entire nanowire and is its integral result. The balance of these flows (Φ^{(+)} and Φ^{(−)}) is inhomogeneous along the nanowire surface (a detailed analysis of this balance was carried out in our work^{47} for Si; see the ESI of ref. 47). For the sinusoidal perturbations of the wire radius, ε ≪ r_{0}, in eqn (4), the value Φ^{(+)} in the narrowing regions is less than that in the widening regions if the wavelength is λ ≲ 10r_{0}. In this case, conditions arise for the development of instability, which is similar to the unstable growth of the shell during the diffusion deposition of the material on the nanowire^{23} (this type of instability may be qualitatively represented by the model developed by Mullins and Sekerka^{49}). It is this mechanism that manifests itself in the dynamics of the nanowire in the [111] orientation. A very important factor is that the surface of the nanowire before the breakup has a scaly structure in the widening zones (see Fig. 2b) composed of fragments of [111] faces (a result of thermal roughening^{50}). Free atoms deposited on such a surface are captured by numerous corner vacancies, which contribute to the further redistribution of atoms to widening regions and the development of instability (see Fig. 5A–C). Turning off this mechanism of atomic transfer in the regime of blocked evaporation (P_{filter} = 0) leads to an expected result that may be obtained on the basis of eqn (9).

Note that the [111] scaly lateral surface (Fig. 2b) corresponds to a rather low surface energy density. At the same time, the slopes of constrictions arising at the developed stage of instability are formed by the fragments of the [111] and [100] faces^{39} with a larger value of σ, which corresponds to β > 0 in eqn (9) and, therefore, λ_{max} > 9. Such a qualitative prediction is in agreement with the data of Fig. 5D: Λ/r_{0} ≈ 12.5.

Returning to the analysis of [110]-nanowires, note that, in contrast to [111]-nanowires, the effect of exchange by atoms between the surface and the near-surface layer leads to a significant increase in the parameter Λ/r_{eff} (see Fig. 4). In this case, the side surface of the nanowire is mainly composed of [111] facets (see Fig. 2(a) and the lower inset in Fig. 4C). Atoms deposited on the mainly flat facets [111] have only three bonds and a high coefficient of surface diffusion. If the length between neighboring necking regions is not sufficient, then these deposited atoms will not have enough time to form new clusters on the side facets. As a result, the surface diffusion flows of bonded atoms directed to the necking regions from left to right (see the lower inset in Fig. 4C, red square) can lead to its shortening, which diminishes the intensity of the detachment process at this region. Thus, the instability may arise only for more elongated necking and widening regions, than in the case of blocked evaporation.

Finally, we analyze the breakup of a [100]-nanowire. Results presented in Fig. 6 show that there is a weak response of the breakup parameter, Λ/r_{eff}, to temperature variations and the effects of surface-vapor exchange in this case. The structure of the lateral nanowire surface before breakup is shown in Fig. 2(c). It is composed of four [100]- and four roughened [110]-facets. Taking into account that the slopes of the necking regions are mainly formed by flat [111] facets^{39} with the lowest surface energy density, we may claim that the value β < 0 in eqn (9), and the length of the excited perturbations must be less than 9r_{0}.

A numerical estimate of λ_{max} in this case can be obtained on the basis of the following assumptions. Since the lateral surface of the nanowire is bounded by (100)- and scaly (110)-facets, the value of σ_{0} in eqn (5) can be set equal to σ_{(100)}. The angle, γ, between the normal to the axis of the nanowire and the normal to the cutting planes (111), which form the slopes of the necking regions, is equal to γ = 0.62 rad. Then

(12) |

This prediction is in agreement with results shown in Fig. 6B at some regions of the nanowire where the deviations Δ = − 1 may be associated with the linear stage of the surface dynamics. Moreover, the similar short-wave, λ_{max}/r_{0} ≈ 5.7, modulations, ε, of the nanowire radius (ε ≈ 0.5r_{0}) were observed in ref. 32. In the final stage of instability (see Fig. 6B), when the nanowire radius is rather visible, λ_{max} ≈ 6.3 according to (x)-Fourier transform. However, these short-wave radius modulations cannot result in single droplets of the corresponding size at the nonlinear stage of breakup, so they either merge into larger-scale clusters (in all the cases under consideration: Λ/r_{eff} ≈ 8.3–9.5) or remain in the metastable state for a while.^{32}

Now, we briefly discuss the reasons for the insensitivity of the [100]-nanowire breakup parameter, Λ/r_{eff}, to temperature changes. For sinusoidal-like short-wave perturbations (λ/r_{eff} ≲ 8 and ε ≲ 0.4r_{0}), the flow density Φ^{(+)} is the largest in the widening regions,^{47} although the surface curvature in these zones, κ_{w}, may be lower than that in the narrowing zones, κ_{n}. If the sublimation flux density from the nanowire surface, Φ^{(−)}, is associated with the local curvature, κ(x), of the surface, then in the narrowing zones Φ^{(−)} > Φ^{(+)}. Additional atomic transport in the near-surface layer accelerates the breakup processes. This effect is realized in the case of the [100]-orientation with the maximal order of rotational symmetry and the sinusoidal-like surface perturbations (see the data presented in Fig. 6A and B).

The slight increase in the wavelength, λ_{max}, in Fig. 6A, compared to Fig. 6B, is a result of the thermal roughening of the nanowire surface caused by its exchange with the surrounding vapor of free atoms, which decreases the level of the anisotropy of σ. Upon decreasing the temperature, this exchange is suppressed, and λ_{max} takes the previous value (see Fig. 6B and C).

As stated above, the exchange by free atoms between the nanowire surface and the near-surface layer may play an ambiguous role. Qualitative analysis of the relevant results presented in Fig. 3, 5D and 6B (the sublimation of atoms from the nanowire surface is absent) is the following.

The development of instability should be accompanied by a decrease in the free energy of the nanowire. In the case of an isotropic density of surface energy, β = 0, the area of the lateral surface, S, should decrease, i.e. ΔS < 0. With the sinusoidal perturbations of the nanowire radius, this happens only for λ > 2πr. The value |ΔS| increases with increasing wavelength and this is reflected in eqn (8a) and (8b) by the factor [1 − ^{2}] on their right-hand side. However, it should be noted that with increasing wavelength, the mass of the transported substance also increases, which is reflected by the factor 1/^{2} on the left-hand side of eqn (8a) and (8b). The optimal ratio between the mentioned values, corresponding to the maximum instability increment, γ, is achieved for λ ≈ 9r (the result is well known in theory^{24,25}). For β ≠ 0, two options are realized.

If the surface energy density, σ, decreases on the slopes of necking regions with an increase in the angle of inclination (β < 0 in eqn (5)), then for the same amplitude ε(t), the effect of an additional decrease in surface energy is more pronounced for shorter wavelength disturbances, which is a result of the σ-anisotropy. The maximum of the increment, γ, shifts toward shorter waves, λ < 9r, and the value of γ itself increases due to the corresponding decrease in the mass of matter transported. Note that for a large value |β|, the optimal ratio between the change in free energy and the mass transferred (see eqn (3)) can be realized for very short wavelengths, when changes in the area of the lateral surface of the nanowire are minimal. However, the excitation of such short waves at the linear stage of instability can lead to the phenomenon that at the nonlinear stage, neighboring proto-droplets will merge into larger clusters. The described effects are demonstrated by the results presented in Fig. 6. Noticeable modulations of the number of atoms in atomic layers in Fig. 6B are developed with a wavelength of λ ≈ 6.3r, which is substantially less than 9r. At the nonlinear stage of breakup, the number of drops in all variants in Fig. 6 is such that Λ ≈ 9.3r. Thus, the value of the breakup parameter, Λ/r, determined in experiments, hides the real dynamics of the process associated with the pronounced anisotropy (β < 0) of σ.

Completely different decay scenarios are realized in cases when β > 0 (in the [111]- and [110]-orientations). The most evident are the features of the breakup of a nanowire with the [110]-orientation. As noted above, its lateral surface is characterized by a minimum average surface energy density (see Fig. 2a). On the surface of the arising necks, the value of σ is higher and in all the presented temperature regimes (see Fig. 3) 9r < λ_{max} ≈ Λ (merging of neighboring drops practically does not occur). An increase in these decay parameters with decreasing temperature (up to Λ/r ≈ 14.5) is associated with the inhomogeneous surface mobility of atoms, which leads to a noticeable deviation of the developing perturbations from the sinusoidal shape and pronounced localization of the breakup regions. This feature of the [110]-orientation was observed in the results of experiments with copper.^{26}

Depending on the material, diameter, and temperature, the formation of a shell of free atoms around a nanowire can result in significant modifications to the dynamics of its decay. The nature of these modifications depends on the intensity of surface diffusion in the narrowing and broadening regions, the sublimation intensity in these regions, and the distribution of diffusion fluxes of free atoms from the surface layer onto the nanowire surface.

Obviously, the dominance of sublimation fluxes, Φ^{(−)}, in the narrowing regions over the reverse flow of free atoms, Φ^{(+)}, contributes to the development of instability that is caused by surface diffusion. However, the realization of such a ratio with the corresponding self-tuning of the surface can lead to a substantial rearrangement of the spectrum of excited perturbations, as the results of effects that are analogous to the roughening transition.^{50} For sinusoidal radius modulations, a dominance of the back-flows, Φ^{(+)}, at the widening regions (as compared to that at the necking regions) can be realized either for short-wave perturbations (λ_{max} < 8r and the shadow effect in diffusion becomes apparent) or for longer perturbations (λ_{max} ≳ 11r, when Φ^{(+)} ∼ 1/r(x)).

The first option is realized in the case of nanowires with the [111] orientation, when the slow surface diffusion along the scaly-surface structure does not suppress the growth of short-scale bulges, and the breakup parameter Λ/r_{eff} is reduced from 12.5 (see Fig. 5B) to the classical value of 9, independent of the temperature regime (see Fig. 5A and B). Thus, the effects of the exchange by atoms between the surface of a nanowire and the enveloping near-surface layer hide the manifestation of anisotropy σ in real experiments.

The second option is clearly manifested in the nanowires with the [110] orientation, the lateral surface of which is mainly formed by (111) faces. The relatively high diffusion rate on such planes inhibits the formation of nanoclusters (new layers in the widening regions) on them. Free atoms settling on these faces drift to the narrowing regions and inhibit their extension along the nanowire. Instabilities arise only with increased distances between necks, which is demonstrated by the results presented in Fig. 4. Note that the effects of sublimation noticeably increase the breakup parameters of nanowires even at a relatively low intensity (see Fig. 4B, configurations (e) and (e′); P_{filter} = 0.3).

Hence, the manifestation of the anisotropy effects of the physical properties of the surface can be clearly demonstrated by a significant deviation of the breakup parameter from its classical value, 9. However, the absence of the noticeable deviation from this value, 9, can also be a result of the complex interaction of various anisotropy factors at different stages of breakup.

- K. Dawson and A. O'Riordan, Towards Nanowire (Bio) Sensors, J. Phys.: Conf. Ser., 2011, 307, 012004 CrossRef.
- J. Cao, T. Sun and K. Grattan, Gold nanorod-based localized surface plasmon resonance biosensors: A review, Sens. Actuators, B, 2014, 195, 332–351 CrossRef CAS.
- K. Dawson, J. Strutwolf, K. Rodgers, G. Herzog, D. Arrigan, A. Quinn and A. O'Riordan, Single Nanoskived Nanowires for Electrochemical Applications, Anal. Chem., 2011, 83(14), 5535–5540 CrossRef CAS PubMed.
- W. Liu, Tunable light absorption of Si nanowires by coating non-absorbing dielectric shells for photovoltaic applications, Optoelectron. Adv. Mater., Rapid Commun., 2011, 5, 916–919 CAS.
- X. Hua, Y. Zeng and W. Shen, Semiconductor/dielectric half-coaxial nanowire arrays for large-area nanostructured photovoltaics, J. Appl. Phys., 2014, 115(12), 124309 CrossRef.
- R. Day, M. Mankin, R. Gao, Y. No, S. Kim, D. Bell, H. Park and C. Lieber, Plateau–Rayleigh crystal growth of periodic shells on one-dimensional substrates, Nat. Nanotechnol., 2015, 10(4), 345–352 CrossRef CAS PubMed.
- Z. Xue, M. Xu, Y. Zhao, J. Wang, X. Jiang, L. Yu, J. Wang, J. Xu, Y. Shi, K. Chen and P. Roca i Cabarrocas, Engineering island-chain silicon nanowires via a droplet mediated Plateau-Rayleigh transformation, Nat. Commun., 2016, 7(1), 12836 CrossRef CAS PubMed.
- I. Goldthorpe, A. Marshall and P. McIntyre, Synthesis and Strain Relaxation of Ge-Core/Si-Shell Nanowire Arrays, Nano Lett., 2008, 8(11), 4081–4086 CrossRef CAS PubMed.
- P. Namdari, H. Daraee and A. Eatemadi, Recent Advances in Silicon Nanowire Biosensors: Synthesis Methods, Properties, and Applications, Nanoscale Res. Lett., 2016, 11(1), 406 CrossRef PubMed.
- S. Dong, P. Hu, X. Zhang, Y. Cheng, C. Fang, J. Xu and G. Chen, Facile synthesis of silicon nitride nanowires with flexible mechanical properties and with diameters controlled by flow rate, Sci. Rep., 2017, 7(1), 45538 CrossRef CAS PubMed.
- Y. Chu, S. Jing, X. Yu and Y. Zhao, High-Temperature Plateau–Rayleigh Growth of Beaded SiC/SiO2 Nanochain Heterostructures, Cryst. Growth Des., 2018, 18(5), 2941–2947 CrossRef CAS.
- Z. Xue, M. Xu, X. Li, J. Wang, X. Jiang, X. Wei, L. Yu, Q. Chen, J. Wang, J. Xu, Y. Shi, K. Chen and P. Roca i Cabarrocas, In-Plane Self-Turning and Twin Dynamics Renders Large Stretchability to Mono-Like Zigzag Silicon Nanowire Springs, Adv. Funct. Mater., 2016, 26(29), 5352–5359 CrossRef CAS.
- S. Hwang, C. Lee, H. Cheng, J. Jeong, S. Kang, J. Kim, J. Shin, J. Yang, Z. Liu, G. Ameer, Y. Huang and J. Rogers, Biodegradable Elastomers and Silicon Nanomembranes/Nanoribbons for Stretchable, Transient Electronics, and Biosensors, Nano Lett., 2015, 15(5), 2801–2808 CrossRef CAS PubMed.
- J. Durham and Y. Zhu, Fabrication of Functional Nanowire Devices on Unconventional Substrates Using Strain-Release Assembly, ACS Appl. Mater. Interfaces, 2012, 5(2), 256–261 CrossRef PubMed.
- W. Wang, C. Chen, K. Lin, Y. Fang and C. Lieber, Label-free detection of small-molecule-protein interactions by using nanowire nanosensors, Proc. Natl. Acad. Sci. U. S. A., 2005, 102(9), 3208–3212 CrossRef CAS PubMed.
- E. Stern, A. Vacic, N. Rajan, J. Criscione, J. Park, B. Ilic, D. Mooney, M. Reed and T. Fahmy, Label-free biomarker detection from whole blood, Nat. Nanotechnol., 2009, 5(2), 138–142 CrossRef PubMed.
- J. Chua, R. Chee, A. Agarwal, S. Wong and G. Zhang, Label-Free Electrical Detection of Cardiac Biomarker with Complementary Metal-Oxide Semiconductor-Compatible Silicon Nanowire Sensor Arrays, Anal. Chem., 2009, 81(15), 6266–6271 CrossRef CAS PubMed.
- H. Li, J. Biser, J. Perkins, S. Dutta, R. Vinci and H. Chan, Thermal stability of Cu nanowires on a sapphire substrate, J. Appl. Phys., 2008, 103(2), 024315 CrossRef.
- M. Rauber, F. Muench, M. Toimil-Molares and W. Ensinger, Thermal stability of electrodeposited platinum nanowires and morphological transformations at elevated temperatures, Nanotechnology, 2012, 23(47), 475710 CrossRef CAS PubMed.
- A. Volk, D. Knez, P. Thaler, A. Hauser, W. Grogger, F. Hofer and W. Ernst, Thermal instabilities and Rayleigh breakup of ultrathin silver nanowires grown in helium nanodroplets, Phys. Chem. Chem. Phys., 2015, 17(38), 24570–24575 RSC.
- S. Gill, Controlling the Rayleigh instability of nanowires, Appl. Phys. Lett., 2013, 102(14), 143108 CrossRef.
- S. Maier, P. Kik, H. Atwater, S. Meltzer, E. Harel, B. Koel and A. Requicha, Local detection of electromagnetic energy transport below the diffraction limit in metal nanoparticle plasmon waveguides, Nat. Mater., 2003, 2(4), 229–232 CrossRef CAS PubMed.
- Y. Yin, Y. Sun, M. Yu, X. Liu, T. Jiang, B. Yang, D. Liu, S. Liu and W. Cao, ZnO Nanorod Array Grown on Ag Layer: A Highly Efficient Fluorescence Enhancement Platform, Sci. Rep., 2015, 5(1), 8152 CrossRef CAS PubMed.
- F. A. Nichols and W. W. Mullins, Surface- (interface-) and volume-diffusion contributions to morphological changes driven by capillarity, Trans. Metall. Soc. AIME, 1965, 233, 1840–1848 CAS.
- J. Plateau, Experimental and Theoretical Statics of Liquids Subject to Molecular Forces Only, Gauthier-Villars, Paris, 1873, vol. 1 Search PubMed.
- M. Toimil Molares, A. Balogh, T. Cornelius, R. Neumann and C. Trautmann, Fragmentation of nanowires driven by Rayleigh instability, Appl. Phys. Lett., 2004, 85(22), 5337–5339 CrossRef CAS.
- C. Bréchignac, P. Cahuzac, F. Carlier, C. Colliex, J. Leroux, A. Masson, B. Yoon and U. Landman, Instability Driven Fragmentation of Nanoscale Fractal Islands, Phys. Rev. Lett., 2002, 88(19), 196103 CrossRef PubMed.
- S. Karim, M. Toimil-Molares, A. Balogh, W. Ensinger, T. Cornelius, E. Khan and R. Neumann, Morphological evolution of Au nanowires controlled by Rayleigh instability, Nanotechnology, 2006, 17(24), 5954–5959 CrossRef CAS.
- S. Karim, M. Toimil-Molares, W. Ensinger, A. Balogh, T. Cornelius, E. Khan and R. Neumann, Influence of crystallinity on the Rayleigh instability of gold nanowires, J. Phys. D: Appl. Phys., 2007, 40(12), 3767–3770 CrossRef CAS.
- M. Rauber, F. Muench, M. Toimil-Molares and W. Ensinger, Thermal stability of electrodeposited platinum nanowires and morphological transformations at elevated temperatures, Nanotechnology, 2012, 23(47), 475710 CrossRef CAS PubMed.
- S. Vigonski, V. Jansson, S. Vlassov, B. Polyakov, E. Baibuz, S. Oras, A. Aabloo, F. Djurabekova and V. Zadin, Au nanowire junction breakup through surface atom diffusion, Nanotechnology, 2017, 29(1), 015704 CrossRef PubMed.
- S. Xu, P. Li and Y. Lu, In situ atomic-scale analysis of Rayleigh instability in ultrathin gold nanowires, Nano Res., 2017, 11(2), 625–632 CrossRef.
- K. Gurski and G. McFadden, The effect of anisotropic surface energy on the Rayleigh instability, Proc. R. Soc. London, Ser. A, 2003, 459(2038), 2575–2598 CrossRef.
- J. Cahn, Stability of rods with anisotropic surface free energy, Scr. Metall., 1979, 13(11), 1069–1071 CrossRef.
- D. Min and H. Wong, Rayleigh's instability of Lennard-Jones liquid nanothreads simulated by molecular dynamics, Phys. Fluids, 2006, 18(2), 024103 CrossRef.
- C. Zhao, J. Sprittles and D. Lockerby, Revisiting the Rayleigh–Plateau instability for the nanoscale, J. Fluid Mech., 2019, 861, R3 CrossRef CAS.
- K. Lai and J. Evans, Reshaping and sintering of 3D fcc metal nanoclusters: Stochastic atomistic modeling with realistic surface diffusion kinetics, Phys. Rev. Mater., 2019, 3(2), 026001 CrossRef CAS.
- P. García-Mochales, R. Paredes, S. Peláez and P. Serena, Statistical Molecular Dynamics Study of (111) and (100) Ni Nanocontacts: Evidences of pentagonal Nanowires, J. Nanomater., 2008, 1–9 Search PubMed.
- V. Gorshkov, P. Sareh, V. Tereshchuk and A. Soleiman-Fallah, Dynamics of Anisotropic Break-Up in Nanowires of FCC Lattice Structure, Adv. Theory Simul., 2019, 1900118 CrossRef.
- V. Gorshkov and V. Privman, Kinetic Monte Carlo model of breakup of nanowires into chains of nanoparticles, J. Appl. Phys., 2017, 122(20), 204301 CrossRef.
- V. Privman, V. Gorshkov and O. Zavalov, Formation of nanoclusters and nanopillars in nonequilibrium surface growth for catalysis applications: growth by diffusional transport of matter in solution synthesis, Heat Mass Transfer, 2013, 50(3), 383–392 CrossRef.
- V. Privman, V. Gorshkov and Y. Yaish, Kinetics modeling of nanoparticle growth on and evaporation off nanotubes, J. Appl. Phys., 2017, 121(1), 014301 CrossRef.
- V. Gorshkov, V. Kuzmenko and V. Privman, Nonequilibrium kinetic modeling of sintering of a layer of dispersed nanocrystals, CrystEngComm, 2014, 16(45), 10395–10409 RSC.
- V. Gorshkov, V. Kuzmenko and V. Privman, Mechanisms of interparticle bridging in sintering of dispersed nanoparticles, J. Coupled Syst. Multiscale Dyn., 2014, 2(2), 91–99 CrossRef.
- I. Galanakis, G. Bihlmayer, V. Bellini, N. Papanikolaou, R. Zeller, S. Blügel and P. Dederichs, Broken-bond rule for the surface energies of noble metals, Europhys. Lett., 2002, 58(5), 751–757 CrossRef CAS.
- G. Grochola, S. Russo and I. Snook, On fitting a gold embedded atom method potential using the force matching method, J. Chem. Phys., 2005, 123(20), 204719 CrossRef PubMed.
- V. Gorshkov, V. Tereshchuk and P. Sareh, Restructuring and breakup of nanowires with the diamond cubic crystal structure into nanoparticles, Mater. Today Commun., 2019, 100727 Search PubMed.
- H. Wang, M. Upmanyu and C. Ciobanu, Morphology of Epitaxial Core−Shell Nanowires, Nano Lett., 2008, 8(12), 4305–4311 CrossRef CAS PubMed.
- W. Mullins and R. Sekerka, Stability of a Planar Interface During Solidification of a Dilute Binary Alloy, J. Appl. Phys., 1964, 35(2), 444–451 CrossRef.
- J. Heyraud, J. Métois and J. Bermond, The roughening transition of the Si{113} and Si{110} surfaces – an in situ, real time observation, Surf. Sci., 1999, 425(1), 48–56 CrossRef CAS.
- D. Chatterjee, S. Shetty, K. Müller-Caspary, T. Grieb, F. Krause, M. Schowalter, A. Rosenauer and N. Ravishankar, Ultrathin Au-Alloy Nanowires at the Liquid–Liquid Interface, Nano Lett., 2018, 18(3), 1903–1907 CrossRef CAS PubMed.

This journal is © The Royal Society of Chemistry 2020 |