[Paper Review] Nonlinear Dynamics of Metallic Nanofabrication

J. Bürki|arXiv (Cornell University)|Jun 16, 2004

Nanofabrication and Lithography Techniques3 citations

TL;DR

This paper proposes a nonlinear dynamical model for metallic nanowire evolution based on a semiclassical energy functional incorporating electron-shell effects and surface self-diffusion. The model predicts universal equilibrium shapes consisting of cylindrical segments connected to unduloid-like leads, explaining the experimentally observed formation of long, stable, nearly perfect cylindrical nanowires via kink-mediated thinning and shape relaxation.

ABSTRACT

Using concepts from fluid dynamics, a partial differential equation for the shape evolution of a metallic nanowire is derived from a semiclassical energy functional that includes electron-shell effects. A rich dynamics, involving movement and interaction of kinks connecting locally stable radii, leads to the formation of a wire whose equilibrium shape is universal and consists of a cylindrical part connected to unduloid-like leads. The universality of the equilibrium shape may provide an explanation for the formation of cylindrical nanowires observed in recent experiments.

Motivation & Objective

• To explain the experimental observation of long, stable, nearly cylindrical Au and Ag nanowires formed in TEM after electron-beam-induced thinning.
• To resolve the apparent paradox of Rayleigh instability suppression in long nanowires with high surface-to-volume ratios.
• To develop a dynamical model that captures the evolution of initially random nanowire shapes toward universal equilibrium configurations.
• To demonstrate the role of electron-shell effects and surface self-diffusion in stabilizing cylindrical nanowires and enabling kink-mediated thinning.

Proposed method

• Derives a partial differential equation for nanowire shape evolution from a semiclassical grand-canonical energy functional including surface tension and mesoscopic electron-shell potential.
• Models surface self-diffusion using Fick’s law and ionic mass conservation, with surface current driven by chemical potential gradients.
• Incorporates a Gutzwiller-type trace formula to compute the electron-shell potential as a function of wire radius and temperature, capturing quantum shell effects.
• Uses periodic boundary conditions and axisymmetric geometry to simplify the dynamics while preserving key features like kink formation and propagation.
• Performs numerical simulations starting from random initial wire profiles to mimic experimental conditions, tracking evolution over time.
• Analyzes equilibrium shapes by rescaling and comparing to Delaunay unduloids to reveal universal geometric features.

Experimental results

Research questions

• RQ1 Why do long metallic nanowires formed via electron-beam thinning in TEM exhibit near-perfect cylindrical shapes despite Rayleigh instability?
• RQ2 How do quantum shell effects stabilize specific wire radii and influence the dynamics of shape evolution?
• RQ3 What is the role of kink formation and propagation in the thinning process of nanowires?
• RQ4 Is the equilibrium shape of a nanowire universal across different initial conditions and system sizes?
• RQ5 How do electron-shell effects and surface diffusion interact to produce stable, cylindrical nanowires with unduloid-like leads?

Key findings

• The equilibrium shape of a nanowire is universal and consists of a cylindrical segment connected to unduloid-like leads, independent of initial conditions.
• The shape of the leads is well approximated by Delaunay unduloids, whose curvature is determined by the ratio of the cylindrical radius to the maximum radius.
• Electron-shell effects stabilize the unduloid shape by pinning it at the junction with the cylindrical part; without these effects, wires break apart due to Rayleigh instability.
• Kink-mediated thinning is observed in simulations, with kinks nucleating near one end and moving at constant speed toward the other end, matching experimental observations.
• The conductance of the final wire decreases by one unit of conductance quantum (G₀) during kink motion, consistent with experimental data (e.g., from G=8G₀ to G=6G₀).
• Simulations of 14 different random initial wires with conductances from 1 to 200G₀ and lengths between 200kFL and 600kFL all converge to the same universal shape upon relaxation.