Multi-Step Sputtering Mechanism

• Multi-step sputtering is a process where sequential sputtering steps—including controlled erosion, surface diffusion, and redeposition—produce hierarchical nanostructures for advanced nanofabrication and analytical applications.
• It utilizes distinct ion-beam orientations, such as oblique and normal incidences, to template and nucleate structures like ripple formations and Au nanobead arrays, leading to long-range order.
• The approach integrates theoretical models (BH, KS, and eKS) with numerical simulations to predict surface evolution, optimize pattern uniformity, and manage contamination in analytical instrumentation.

The multi-step sputtering mechanism encompasses a class of physical and technological processes in which material is removed, transported, and redeposited through a controlled sequence of sputtering steps. Each step—characterized by distinct incident particle energies, geometries, or temporal orderings—enables functions ranging from nanostructuring at surfaces to the formation of ultrathin films or control of contamination in analytic instrumentation. The underlying processes feature a coupling between curvature-dependent erosion, surface diffusion, nonlinear mass redistribution, and local redeposition, and they may be exploited in both single-sample fabrication and multi-sample analytical contexts.

1. Fundamental Principles of Multi-Step Sputtering

In multi-step sputtering, material removal and structuring involve not only the primary ion-induced ejection of atoms (sputter erosion), but also subsequent redeposition events or erosion under new boundary conditions. These steps can be realized as sequential exposure to differently oriented ion beams, alternations between oblique and normal incidence, or the use of sample rotation, among other protocols.

The physical driving forces can be mapped to linear instability (curvature-dependent roughening) as described by the Bradley–Harper (BH) model, as well as nonlinear stabilization and coarsening effects. When appropriately sequenced or superposed, these steps generate hierarchical or composite nanostructures unattainable in one-step processes. In analytical instrumentation, multi-step pathways may entail both the initial deposition on inert surfaces in the environment and subsequent reflective sputtering, resulting in unintended cross-contaminations or, if harnessed, engineered coatings (Kim et al., 2011, Scheithauer, 2013, Yasseri et al., 2011).

2. Sequential Surface Patterning and Hierarchical Self-Assembly

Hierarchical bottom-up self-assembly by multi-step sputtering has been functionally validated in the formation of one-dimensional arrays of Au nanobeads on Au(001) surfaces (Kim et al., 2011). The canonical protocol consists of:

1. Step 1: Ripple Template Formation
   • Oblique-incidence Ar⁺ sputtering (E = 2 keV, θ = 72° from surface normal, P_Ar = 1.2×10⁻⁴ Torr, T ≈ 300 K, ψ = 4,500 ions nm⁻²), producing ripples of wavelength λ_R ≃ 47 nm.
2. Step 2: Nanobead Nucleation and Growth
   • Normal-incidence sputtering on the pre-patterned rippled surface (same ion species/energy; ψ₂ = 356–1,781 ions nm⁻²) induces bead formation selectively at the ripple crests. As fluence increases, bead diameter/spacing λ_B evolves from ≃35 nm to ≃50 nm, amplitude A_B from ≃0.5 nm to 1.5 nm, and coherence length ℓ_B surpasses 1 µm.

The two-step mechanism leverages the anisotropy introduced in Step 1 as a spatial template dictating mass transport and bead nucleation in Step 2, achieving long-range order and uniformity unattainable by normal-incidence sputtering alone.

3. Theoretical and Numerical Frameworks

The evolution of surface topography under multi-step sputtering is classically described by continuum equations. The pivotal models are:

• Linear Bradley–Harper (BH) Model:

$\partial_t h = -\nu \nabla^2 h - K \nabla^4 h$

Predicts ripple formation with fastest-growing wavelength $\lambda^* = 2\pi\sqrt{2K/\nu}$.

• Kuramoto–Sivashinsky (KS) Model:

Adds nonlinear, non-conserved term $(\nabla h)^2$.

• Extended KS (eKS) Model with Conserved KPZ (cKPZ) Term:

$$\partial_t h = - \nu \nabla^2 h - K \nabla^4 h + \lambda^{(1)} (\nabla h)^2 + \lambda^{(2)} \nabla^2 (\nabla h)^2 + \eta(\mathbf{r}, t)$$

Here, $\lambda^{(2)} \nabla^2 (\nabla h)^2$ models local redeposition and surface-confined transport. The dimensionless parameter $r = (\nu\lambda^{(2)})/(K\lambda^{(1)})$ quantifies the relative strength of mass-conserving nonlinearity, crucial for order preservation (Kim et al., 2011).

Numerical integration confirms that only with a nonzero cKPZ term ($r \lesssim -0.5$) does the system reproduce long-lived, highly ordered nanobead arrays. Without this term, order is rapidly lost.

4. Multi-Step Sputtering in Analytical Instrumentation

In AES and XPS, multi-step sputtering explains cross-contamination between samples through a tripartite mechanism (Scheithauer, 2013):

1. Primary Sputtering: Energetic Ar⁺ ions eject atoms from a primary target, with yields $$Y_\mathrm{Au}(2\,\mathrm{keV},45^\circ) \simeq 2\text{–}3$$ atoms/ion.
2. Deposition on Chamber Surfaces: Sputtered atoms deposit on nearby walls and fixtures, forming a contaminant layer.
3. Secondary (Reflective) Sputtering: Subsequent high-energy Ar reflections or energetic sputtered atoms re-sputter this deposit, leading to secondary redeposition on neighboring samples. The coverage and film formation follow a rate equation driven by the secondary flux $J_\mathrm{sec} = R Y_\mathrm{prim} j_\mathrm{Ar}$. The process enables in-situ ultrathin coatings (e.g., 1–2 ML) for improved electrical conductivity and energy referencing.

Table: Stages of Reflective Sputtering Cross-Contamination

Step	Description	Key Quantities
Primary	Ion-beam sputtering of target	$Y_\mathrm{prim}$ (yield), $j_\mathrm{Ar}$ (flux)
Deposition	Build-up on chamber surfaces	Layer $\sim$35 nm (Al-foil window observed)
Secondary	Re-sputtering to new samples	$J_\mathrm{sec}$, $R$ (reflectivity)

5. Kinetic Modeling and Pattern Evolution

Kinetic Monte Carlo simulations provide mechanistic insight into pattern evolution under dual, sequential, and rotated beam geometries (Yasseri et al., 2011):

• Sequential Beam Sputtering (SIBS):

Upon rotation by Δφ, pre-existing ripples collapse quickly, with structural relaxation time $T_0 \approx 2.5$ ML. This transient leads to a roughness dip and ephemeral superposition of structures from both beam directions.

• Rotated Ion-Beam Sputtering (RIBS):

The roughness exhibits a pronounced minimum at rotation frequency $\omega \approx 1/T_0$. At low ω, ripples adiabatically follow, while at high ω, dot-like patterns emerge.

• Dual Beam Sputtering (DIBS):

Exactly balanced, orthogonal beams yield square patterns. Any imbalance or non-orthogonality leads to ripples oriented along predictable principal axes.

The continuum BH framework captures the orientation dynamics through an effective erosion-rate tensor, whose eigenvalues and vectors determine dominant morphologies.

6. Integrated Simulation Strategies

A multi-step sputtering process in device-scale environments requires chained multi-physics modeling (Gentile, 2012):

1. Plasma Initiation: Self-consistent plasma potentials, electron and ion densities established via PIC or fluid-MC solvers.
2. Ion Acceleration: Sheath transport and collision modeling produce an energy-angle distribution onto the target.
3. Sputter-Ejection: Binary collision MC and molecular dynamics yield energy–angle distributions of ejected atoms.
4. Gas-Phase Transport: DSMC modeling of neutral atom propagation, thermalization, and arrival fluxes.
5. Film Growth: Continuum/KMC models predict adsorption, diffusion, and mechanical/thermal properties of the film.

These modules exchange statistical distributions at well-defined interfaces. The decoupling of fast plasma and slower film evolution enables practical simulation of complex multi-step deposition sequences.

7. Applications, Limitations, and Future Directions

Multi-step sputtering enables the fabrication of hierarchically ordered nanostructures, in-situ calibration layers in spectroscopy, and tailored ultrathin coatings for electrical and analytical functions. The preservation of order and sharpness in these assemblies critically depends on local redeposition and surface-confined transport—the nonlinear conservation dynamics quantified in the eKS formalism (Kim et al., 2011). Reflective sputtering in AES/XPS, previously a source of contamination, now serves as a controlled route to monolayer-scale coverage (Scheithauer, 2013).

A plausible implication is that extending these approaches to additional materials, complex chamber geometries, or dynamic substrate topographies will require further advancements in coupled simulation methods and in-situ diagnostics for validation. Accurate integration of erosion, transport, and co-sputtering effects remains a central challenge for comprehensive predictive modeling (Gentile, 2012, Yasseri et al., 2011).