Ultrafast Physical Ablation

• Ultrafast physical ablation is a process using sub-picosecond laser pulses to non-thermally remove material with minimal heat diffusion.
• It enables precision micro- and nanomachining across metals, semiconductors, insulators, and 2D materials through nonlinear absorption and energy localization.
• Quantitative metrics like ablation thresholds and incubation effects inform process optimization and defect engineering for deterministic material removal.

Ultrafast physical ablation is a regime of material removal driven by intense, sub-picosecond laser pulses that deposit energy on timescales shorter than thermal diffusion or lattice relaxation, leading to highly localized, non-thermal, and often deterministic elimination of material. This process underpins a diverse set of applications in micro- and nanomachining, device prototyping, direct-write lithography, and precision surface engineering across metals, semiconductors, insulators, and atomically thin materials. The following sections survey the principal mechanisms, quantitative metrics, material-specific behaviors, spatial characteristics, and engineering implications of ultrafast physical ablation, synthesizing numerical and experimental results from recent literature.

1. Physical Mechanisms and Timescales

Ultrafast ablation is initiated by the nonlinear absorption of laser photons (typically in the femtosecond–picosecond regime), rapidly driving the electronic subsystem of the target far from equilibrium. Key stages include:

• Nonlinear Ionization and Carrier Excitation: At intensities exceeding 10¹¹–10¹² W/cm², multiphoton and tunneling ionization dominate, generating dense electron–hole plasmas on 10–100 fs timescales in materials like graphene, semiconductors, and metals (Vasquez et al., 2019, McDonnell et al., 2020, Werner et al., 2019, Klein et al., 5 Feb 2026).
• Carrier–Lattice Energy Transfer: Electron–phonon coupling transmits energy to the lattice on 1–10 ps timescales. In metals and degenerate conductors, this process controls the subsequent lattice heating; in covalent solids and 2D materials, strong excitation-induced bond weakening alters or bypasses thermal pathways (McDonnell et al., 2020, Klein et al., 5 Feb 2026).
• Non-Thermal Bond Disruption and Material Removal: At high electronic temperatures or carrier densities, energy localization leads to non-thermal melting, electron-blast-driven lattice deformation, or even antibonding-states-driven bond scission. These regimes enable evaporation, spallation, or void formation without significant heating of the bulk lattice (McDonnell et al., 2020, Klein et al., 5 Feb 2026, Werner et al., 2019).
• Plasma Expansion, Shock, and Cavitation: In liquids and transparent materials, the rapid conversion of deposited energy produces explosive vaporization, pressure waves, and shock fronts (Sedov–Taylor expansion), observable on 10 ps–10 ns timescales (Hernandez-Rueda et al., 2018, Zhang et al., 17 Apr 2025).

The fundamental advantage of ultrafast ablation lies in its ability to spatially confine material modification and restrict the heat-affected zone to sub-micron scales, leveraging timescales too short for significant thermal diffusion.

2. Ablation Thresholds: Quantitative Metrics and Material Dependence

The ablation threshold is the minimum incident fluence required to cause irreversible material removal under a specified set of conditions (pulse duration, wavelength, sample geometry, and, for 2D materials, substrate structure).

Material/System	Pulse Duration	Wavelength	Ablation Threshold F_th	Notes
Graphene (monolayer)	100 fs	819 nm	~9.2 mJ/cm² (multi-shot, 15s)	F_th(1 pulse) ≈16.7 mJ/cm², reduces to ≈9.2 mJ/cm² at N≫10⁶ (Vasquez et al., 2019)
ITO on glass (20 nm)	500 fs	1030 nm	7.50 mJ/cm² (absorbed)	Dominated by non-thermal electron-blast ejection (McDonnell et al., 2020)
MoS₂ (monolayer)	160 fs	800 nm	F₁=111 mJ/cm² (single pulse)	Intrinsic substrate-independent F_nt = 66 mJ/cm² (Solomon et al., 2021)
Silicon (single crystal)	200 fs	2.75–4.15 μm	0.52–0.87 J/cm² (ablation)	Ponderomotive-driven, non-thermal ablation in MIR (Werner et al., 2019)
Water (liquid interface)	150 fs	800 nm	8.1 J/cm²	Threshold for explosive vaporization (Hernandez-Rueda et al., 2018)

Thresholds are established via area–fluence fitting (Liu method) or via intrinsic-field normalization in atomically thin systems. Multi-shot incubation typically lowers F_th with increasing pulse count N, but robust materials like monolayer MoS₂ show remarkably small and rapidly saturating reductions (Solomon et al., 2021, Vasquez et al., 2019).

3. Distinct Regimes: Non-Thermal, Thermal, and Hybrid Ablation

Ablation manifests as a function of fluence relative to F_th, pulse duration, and material electron–phonon coupling:

• Non-Thermal/Electron-Blast Regimes: Below thresholds associated with lattice melting, grain ejection, surface evaporation, and defect formation are dominated by ultrafast electronic excitation. Characteristic outcomes include minimal melt, preservation of substrate or contiguous layers, and the absence of resolidification features (e.g., in ITO and Ti/Au thin films) (McDonnell et al., 2020, Kim et al., 2020).
• Thermal/Vaporization Regimes: At fluences exceeding ~1 J/cm² in most materials (excluding 2D systems), lattice temperatures reach melt or vaporization points, and classical phase-explosion or melt ejection is observed (Werner et al., 2019, Hernandez-Rueda et al., 2018).
• Hybrid Regimes and Phase Diagram Shifts: In covalent semiconductors (Si, 2D TMDs), excitation-dependent bond weakening shifts melting and ablation boundaries to lower lattice temperatures as electronic temperature increases, enabling non-thermal processes even when T_l is far below T_m at equilibrium (Klein et al., 5 Feb 2026, Solomon et al., 2021).

Distinct sub-ablation damage thresholds—where irreversible point or line defects, lattice disorder, or nano-voids are induced without full removal—are quantitatively observed in monolayer MoS₂ (78% of F₁) and graphene (~75% F_th), providing windows for defect engineering (Solomon et al., 2021, Vasquez et al., 2019).

4. Incubation Effects and Cumulative Damage

Multi-pulse irradiation can lower the effective ablation threshold (incubation), driven by both defect-seeded absorption enhancement and bond softening. Empirical and microscopic models reveal:

• Power-law or Exponential Saturation: F_th(N) ≈ F_∞ + (F_1 – F_∞)*N⁻ᵏ or F_th(N)=F_th(∞)+[F_th(1)–F_th(∞)] exp(–γN).
• Physical Origins: Absorption rise (ΔA) and critical surface energy reduction (ΔG′) due to accumulated vacancies or disorder (Solomon et al., 2021).
• Material-Specific Saturation: MoS₂ displays ΔG′/G₀′ ≈ –0.013, ΔA/A₀ ≈ 0.31, leading to R = F_∞/F₁ ≈ 0.75 (only ~25% reduction and N_sat ≈ 10 shots). By contrast, graphene and WS₂ have larger reductions and slower saturation, with N_sat ≫ 10³ (Solomon et al., 2021, Vasquez et al., 2019).

This behavior underlies the robustness of certain 2D materials in high-repetition or high-throughput laser patterning contexts.

5. Patterning Resolution and Spatial Characteristics

Advances in beam engineering, material–optic coupling, and control of energy localization govern the achievable spatial resolution and aspect ratios:

• Diffraction-Limited and Beyond: In monolayer crystals, ablation feature size is set primarily by the focused spot, with MoS₂ trenches as narrow as 0.52 µm at NA = 0.55 (Solomon et al., 2021). In transparent solids, use of Bessel–Gauss beams with back-surface processing and self-seeded near-field enhancement allows nanochannel structuring down to 7 nm—exceeding the λ/100 limit and delivering aspect ratios >1,000 (Zhang et al., 17 Apr 2025).
• Radial Defect Zones: For graphene and other single-layer targets, a defect “halo” (e.g., width ~2 µm, L_D ≈ 58 nm) surrounds cleanly ablated holes, with gradations in disorder controlled by local fluence distributions (Vasquez et al., 2019).
• Near-Field/Plasmonic Enhancement: Substrate choice and optical stack tuning (e.g., DBRs, plasmonic films) can dramatically lower the required incident fluence via increased local field intensity (η), enabling low-power ablation at high speeds (Solomon et al., 2021).

6. Process Engineering: Pulse Structuring, Efficiency, and Application Domains

Pulse temporal structuring (single-pulse, burst, biburst) and beam-size optimization substantially affect ablation efficiency, removal selectivity, and collateral impact:

• Efficiency Maximization: For metals (copper, steel), ablation efficiency peaks in MHz-burst (N ≈ 3) regimes (η ≃ 8.8 μm³/μJ for copper drilling). Pure GHz bursts reduce efficiency due to plasma shielding, while “biburst” (MHz-in-GHz) can provide modest gains and mitigate odd-even oscillations (Žemaitis et al., 2020).
• Selective Material Removal: In metal multilayers, incorporation of a high electron–phonon coupling transition metal interlayer (e.g., Ti with G ≈ 2.6×10¹⁸ W m⁻³ K⁻¹ vs. Au with G ≈ 2.6×10¹⁶ W m⁻³ K⁻¹) enables selective ultrafast vaporization and ultra-clean lift-off of noble metal films at fluences ~3.2–3.9 J/cm², without substrate damage (Kim et al., 2020).
• Soft Tissue and Fiber Delivery: In biomedical ablation, burst-mode femtosecond pulses (e.g., N = 4, E_p = 3–8 μJ, intra-burst rate 22.3 MHz) reveal ablation rates of ~1 mm³/min with negligible collateral heating, compatible with flexible fiber delivery (Kerse et al., 2014).
• Adiabatic 2D Material Structuring: Due to negligible phonon coupling to substrates on the few-ps timescale, the ablation of 2D van der Waals systems is fundamentally substrate-adiabatic, governed by the internal field (optical etalon effect), allowing for predictable, uniform removal independent of substrate thermal properties (Solomon et al., 2021, Solomon et al., 2021).

7. Modeling, Predictive Frameworks, and Future Directions

Modern treatment of ultrafast ablation couples continuum two-temperature models (TTM) or thermal-spike models (TSM) with MD using excitation-dependent, parameterized bond-order potentials (Klein et al., 5 Feb 2026, McDonnell et al., 2020, Kim et al., 2020). For light–matter interactions in 2D materials, analytic closed-form treatments employing the zero-thickness approximation and etalon field enhancement allow explicit prediction of thresholds and field profiles (Solomon et al., 2021).

Ongoing research elucidates the universality of non-thermal mechanisms—such as excitation-mediated bond weakening and void-driven spallation—in crystals beyond silicon (e.g., diamond, chalcogenides), and pushes the ultimate limits of resolution, pattern fidelity, and process flexibility (Klein et al., 5 Feb 2026, Zhang et al., 17 Apr 2025). A plausible implication is that future application realms, including atomically precise device engineering, metamaterial fabrication, and in-vivo bioprocessing, will leverage physical ablation regimes precisely tailored by temporal pulse structuring, optical environment design, and atomistic modeling.