Chengping Wu

Atomistic modeling of nanoparticle generation in short pulse laser ablation of thin metal films in water.

Laser ablation in liquids is actively used for generation of clean colloidal nanoparticles with unique shapes and functionalities. The fundamental mechanisms of the laser ablation in liquids and the key processes that control the nanoparticle structure, composition, and size distribution, however, are not yet fully understood. In this paper, we report the results of first atomistic simulations of laser ablation of metal targets in liquid environment. A model combining a coarse-grained representation of the liquid environment (parameterized for water), a fully atomistic description of laser interactions with metal targets, and acoustic impedance matching boundary conditions is developed and applied for simulation of laser ablation of a thin silver film deposited on a silica substrate. The simulations, performed at two laser fluences in the regime of phase explosion, predict a rapid deceleration of the ejected ablation plume and the formation of a dense superheated molten layer at the water-plume interface. The water in contact with the hot metal layer is brought to the supercritical state and transforms into an expanding low density metal-water mixing region that serves as a precursor for the formation of a cavitation bubble. Two distinct mechanisms of the nanoparticle formation are predicted in the simulations: (1) the nucleation and growth of small (mostly ⩽10nm) nanoparticles in the metal-water mixing region and (2) the formation of larger (tens of nm) nanoparticles through the breakup of the superheated molten metal layer triggered by the emergence of complex morphological features attributed to the Rayleigh-Taylor instability of the interface between at the superheated metal layer and the supercritical water. The first mechanism is facilitated by the rapid cooling of the growing nanoparticles in the supercritical water environment, resulting in solidification of the nanoparticles located in the upper part of the mixing region on the timescale of nanoseconds. The computational prediction of the two mechanisms of nanoparticle formation yielding nanoparticles with different characteristic sizes is consistent with experimental observations of two distinct nanoparticle populations appearing at different stages of the ablation process.

Nanoparticle generation and transport resulting from femtosecond laser ablation of ultrathin metal films: Time-resolved measurements and molecular dynamics simulations

The synthesis of metal nanoparticles by ultrafast laser ablation of nanometers-thick metal films has been studied experimentally and computationally. Near-threshold backside laser ablation of 2–20 nm-thick Pt films deposited on fused silica substrates was found to produce nanoparticles with size distributions that were bimodal for the thicker films, but collapsed into a single mode distribution for the thinnest film. Time-resolved imaging of blackbody emission from the Pt nanoparticles was used to reveal the nanoparticle propagation dynamics and estimate their temperatures. The observed nanoparticle plume was compact and highly forward-directed with a well-defined collective velocity that permitted multiple rebounds with substrates to be revealed. Large-scale molecular dynamics simulations were used to understand the evolution of compressive and tensile stresses in the thicker melted liquid films that lead to their breakup and ejection of two groups of nanoparticles with different velocity and size distributions. Ultrafast laser irradiation of ultrathin (few nm) metal films avoids the splitting of the film and appears to be a method well-suited to cleanly synthesize and deposit nanoparticles from semitransparent thin film targets in highly directed beams.

Growth Twinning and Generation of High-Frequency Surface Nanostructures in Ultrafast Laser-Induced Transient Melting and Resolidification

The structural changes generated in surface regions of single crystal Ni targets by femtosecond laser irradiation are investigated experimentally and computationally for laser fluences that, in the multipulse irradiation regime, produce sub-100 nm high spatial frequency surface structures. Detailed experimental characterization of the irradiated targets combining electron back scattered diffraction analysis with high-resolution transmission electron microscopy reveals the presence of multiple nanoscale twinned domains in the irradiated surface regions of single crystal targets with (111) surface orientation. Atomistic- and continuum-level simulations performed for experimental irradiation conditions reproduce the generation of twinned domains and establish the conditions leading to the formation of growth twin boundaries in the course of the fast transient melting and epitaxial regrowth of the surface regions of the irradiated targets. The observation of growth twins in the irradiated Ni(111) targets provides strong evidence of the role of surface melting and resolidification in the formation of high spatial frequency surface structures. This also suggests that the formation of twinned domains can be used as a sensitive measure of the levels of liquid undercooling achieved in short pulse laser processing of metals.

Atomistic simulation study of short pulse laser interactions with a metal target under conditions of spatial confinement by a transparent overlayer

The distinct characteristics of short pulse laser interactions with a metal target under conditions of spatial confinement by a solid transparent overlayer are investigated in a series of atomistic simulations. The simulations are performed with a computational model combining classical molecular dynamics (MD) technique with a continuum description of the laser excitation, electron-phonon equilibration, and electronic heat transfer based on two-temperature model (TTM). Two methods for incorporation of the description of a transparent overlayer into the TTM-MD model are designed and parameterized for Ag-silica system. The material response to the laser energy deposition is studied for a range of laser fluences that, in the absence of the transparent overlayer, covers the regimes of melting and resolidification, photomechanical spallation, and phase explosion of the overheated surface region. In contrast to the irradiation in vacuum, the spatial confinement by the overlayer facilitates generation of sustain...

Generation of Subsurface Voids, Incubation Effect, and Formation of Nanoparticles in Short Pulse Laser Interactions with Bulk Metal Targets in Liquid: Molecular Dynamics Study

The ability of short pulse laser ablation in liquids to produce clean colloidal nanoparticles and unusual surface morphology has been employed in a broad range of practical applications. In this paper, we report the results of large-scale molecular dynamics simulations aimed at revealing the key processes that control the surface morphology and nanoparticle size distributions by pulsed laser ablation in liquids. The simulations of bulk Ag targets irradiated in water are performed with an advanced computational model combining a coarse-grained representation of liquid environment and an atomistic description of laser interaction with metal targets. For the irradiation conditions that correspond to the spallation regime in vacuum, the simulations predict that the water environment can prevent the complete separation of the spalled layer from the target, leading to the formation of large subsurface voids stabilized by rapid cooling and solidification. The subsequent irradiation of the laser-modified surface is found to result in a more efficient ablation and nanoparticle generation, thus suggesting the possibility of the incubation effect in multipulse laser ablation in liquids. The simulations performed at higher laser fluences that correspond to the phase explosion regime in vacuum reveal the accumulation of the ablation plume at the interface with the water environment and the formation of a hot metal layer. The water in contact with the metal layer is brought to the supercritical state and provides an environment suitable for nucleation and growth of small metal nanoparticles from metal atoms emitted from the hot metal layer. The metal layer itself has limited stability and can readily disintegrate into large (tens of nanometers) nanoparticles. The layer disintegration is facilitated by the Rayleigh–Taylor instability of the interface between the higher density metal layer decelerated by the pressure from the lighter supercritical water. The nanoparticles emerging from the layer disintegration are rapidly cooled and solidified due to the interaction with water environment, with a cooling rate of ∼2 × 1012 K/s observed in the simulations. The computational prediction of two distinct mechanisms of nanoparticle formation yielding nanoparticles with different characteristic sizes provides a plausible explanation for the experimental observations of bimodal nanoparticle size distributions in laser ablation in liquids. The ultrahigh cooling and solidification rates suggest the possibility for generation of nanoparticles featuring metastable phases and highly nonequilibrium structures.

Atomic Movies of Laser-Induced Structural and Phase Transformations from Molecular Dynamics Simulations

Molecular dynamics (MD) simulations of laser-materials interactions are playing an important role in investigation of complex and highly non-equilibrium processes involved in short pulse laser processing and surface modification. This role is defined by the ability of MD simulations to reveal in-depth information on the structural and phase transformations induced by the laser excitation and, at the same time, to provide clear visual representations, or “atomic movies,” of laser-induced dynamic processes. This chapter provides a brief overview of recent progress in the description of laser coupling and relaxation of photo-excited states in metals, semiconductors, insulators and molecular systems within the general framework of the classical MD technique and presents several examples of MD simulations of laser melting, generation of crystal defects, photomechanical spallation, explosive boiling and molecular entrainment in laser ablation. Possible directions of future progress in atomistic modeling of laser-materials interactions and the potential role of MD simulations in the design of an integrated multiscale computational model capable of accounting for interrelations between processes occurring at different time- and length-scales are discussed.

Molecular Dynamics Simulations of Laser-Materials Interactions: General and Material-Specific Mechanisms of Material Removal and Generation of Crystal Defects

Molecular dynamics simulations of laser-materials interactions are capable of providing detailed information on the complex processes induced by the fast laser energy deposition and can help in the advancement of laser-driven applications. This chapter provides a brief overview of recent progress in the atomic- and molecular-level modeling of laser-materials interactions and presents several examples of the application of atomistic simulations for investigation of laser melting and resolidification, generation of crystal defects, photomechanical spallation, and ablation of metals and molecular targets. A particular focus of the analysis of the computational results is on revealing the general and material-specific phenomena in laser-materials interactions and on making connections to experimental observations.

Mechanism of acoustically induced diffusional structuring of surface adatoms

Physical mechanisms of time-averaged structuring of adatoms induced by a standing surface acoustic wave (SAW) on a solid substrate are studied. Despite some similarity with conventional mechanisms based on averaging of fast oscillation-type motion or radiation-pressure effects, we demonstrate that, for diffusional (i.e., strongly damped) adatom motion, the origin of time-averaged structuring is essentially different. The proposed analytical model and kinetic Monte–Carlo (kMC) simulations reveal several distinct structuring regimes and directly relate them to the transient modification of diffusion barriers and adiabatic temperature variations induced by SAW strains.

VIBRATIONAL CONTRIBUTION TO THERMAL CONDUCTIVITY OF SILICON NEAR SOLID-LIQUID TRANSITION

Although thermal transport in silicon is dominated by phonons in the solid state, electrons also participate as the system approaches, and exceeds, its melting point. Thus, the contribution from both phonons and electrons must be considered in any model for the thermal conductivity, k, of silicon near the melting point. In this paper, equilibrium molecular dynamics simulations measure the vibration mediated thermal conductivity in Stillinger-Weber silicon at temperatures ranging from 1400 to 2000 K — encompassing the solid-liquid phase transition. Non-equilibrium molecular dynamics is also employed as a confirmatory study. The electron contribution may then be estimated by comparing these results to experimental measurements of k. The resulting relationship may provide a guide for the modeling of heat transport under conditions realized in high temperature applications, such as laser irradiation or rapid thermal processing of silicon substrates.Copyright