A theoretical framework to study linear and nonlinear Richtmyer-Meshkov instability (RMI) is presented. This instability typically develops when an incident shock crosses a corrugated material interface separating two fluids with different thermodynamic properties. Because the contact surface is rippled, the transmitted and reflected wavefronts are also corrugated, and some circulation is generated at the material boundary. The velocity circulation is progressively modified by the sound wave field radiated by the wavefronts, and ripple growth at the contact surface reaches a constant asymptotic normal velocity when the shocks/rarefactions are distant enough. The instability growth is driven by two effects: an initial deposition of velocity circulation at the material interface by the corrugated shock fronts and its subsequent variation in time due to the sonic field of pressure perturbations radiated by the deformed shocks. First, an exact analytical model to determine the asymptotic linear growth rate is presented and its dependence on the governing parameters is briefly discussed. Instabilities referred to as RM-like, driven by localized non-uniform vorticity, also exist; they are either initially deposited or supplied by external sources. Ablative RMI and its stabilization mechanisms are discussed as an example. When the ripple amplitude increases and becomes comparable to the perturbation wavelength, the instability enters the nonlinear phase and the perturbation velocity starts to decrease. An analytical model to describe this second stage of instability evolution is presented within the limit of incompressible and irrotational fluids, based on the dynamics of the contact surface circulation. RMI in solids and liquids is also presented via molecular dynamics simulations for planar and cylindrical geometries, where we show the generation of vorticity even in viscid materials.

By extending molecular dynamics simulation to the experimental micrometer length scale, we observed, for the rst time, the complete dynamics of gold lms subjected to ultrashort (1 ps) laser irradiation, culminating in cavitation then ablation of the melt at the front, and crack nucleation then spallation at the rear side of the sample. For thick lms (> 0.5 µm) ablation and spallation are spatially separated; whereas for thin lms, they merge into a unique damage process, which leads to a dependence of ablation threshold on lm thickness. This work provides new insights into thermomechanical dynamics of irradiated metal lms, including coupling-decoupling of ablation and spallation processes, and opens up a eld of direct comparison between experiment and atomistic simulations.

By decoupling time and length scales in moving window molecular dynamics shock-wave simulations, a new regime of shock-wave propagation is uncovered characterized by a two-zone elastic-plastic shock-wave structure consisting of a leading elastic front followed by a plastic front, both moving with the same average speed and having a fixed net thickness that can extend to microns. The material in the elastic zone is in a metastable state that supports a pressure that can substantially exceed the critical pressure characteristic of the onset of the well-known split-elastic-plastic, two-wave propagation. The two-zone elastic-plastic wave is a general phenomenon observed in simulations of a broad class of crystalline materials and is within the reach of current experimental techniques.

Nanoscale thermally assisted hydrodynamic melt perturbations induced by ultrafast laser energy deposition in noble-metal films produce irreversible nanoscale translative mass redistributions and results in formation of radially-symmetric frozen surface structures. We demonstrate that the final three-dimensional (3D) shape of the surface structures formed after re-solidification of the molten part of the film is governed by incident laser fluence and, more importantly, predicted theoretically via molecular dynamics modeling. Considering the underlying physical processes associated with laser-induced energy absorption, electron-ion energy exchange, acoustic relaxation and hydrodynamic flows, the theoretical approach separating "slow" and "fast" physical processes and combining hybrid analytical two-temperature calculations, scalable molecular-dynamics simulations, and a semi-analytical thin-shell model was shown to provide accurate prediction of the final nanoscale solidified morphologies, fully consistent with direct electron-microscopy visualization of nanoscale focused ion-beam cuts of the surface structures produced at different incident laser fluences. Finally, these results are in reasonable quantitative agreement with mass distribution profiles across the surface nanostructures, provided by their noninvasive and non-destructive nanoscale characterization based on energy-dispersive x-ray fluorescence spectroscopy, operating at variable electron-beam energies.

No abstract

Laser-produced surface nanostructures show considerable promise for many applications while fundamental questions concerning the corresponding mechanisms of structuring are still debated. Here, we present a simple physical model describing those mechanisms happened in a thin metal film on dielectric substrate irradiated by a tightly focused ultrashort laser pulse. The main ingredients included into the model are (i) the film–substrate hydrodynamic interaction, melting and separation of the film from substrate with velocity increasing with increase of absorbed fluence; (ii) the capillary forces decelerating expansion of the expanding flying film; and (iii) rapid freezing into a solid state if the rate of solidification is comparable or larger than hydrodynamic velocities. The developed model and performed simulations explain appearance of microbump inside the focal spot on the film surface. The model follows experimental findings about gradual transformation of the bump from small parabolic to a conical shape and to the bump with a jet on its tip with increasing fluence. Disruption of the bump as a result of thinning down the liquid film to a few interatomic distances or due to mechanical break-off of solid film is described together with the jetting and formation of one or many droplets. Developed theory opens door for optimizing laser parameters for intended nanostructuring in applications.

We report an experimental and numerical investigation of the fragmentation mechanisms of micrometer-sized metal droplet irradiated by ultrashort laser pulses. The results of the experiment show that the fast one-side heating of such a droplet may lead to either symmetric or asymmetric expansion followed by different fragmentation scenarios. To unveil the underlying processes leading to fragmentation we perform simulation of liquid-tin droplet expansion produced by the initial conditions similar to those in experiment using the smoothed particle hydrodynamics (SPH) method. Simulation demonstrates that a thin heated surface layer generates a ultrashort shock wave propagating from the frontal side to rear side of the droplet. Convergence of such shock wave followed by a rarefaction tale to the droplet center results in the cavitation of material inside the central region by the strong tensile stress. Reflection of the shock wave from the rear side of droplet produces another region of highly stretched material where the spallation may occur producing a thin spallation layer moving with a velocity higher than expansion of the central shell after cavitation. It is shown both experimentally and numerically that the threshold laser intensity necessary for the spallation is higher than the threshold required to induce cavitation in the central region of droplet. Thus, the regime of asymmetrical expansion is realized if the laser intensity exceeds the spallation threshold. The transverse and longitudinal expansion velocities obtained in SPH simulations of different regimes of expansion are agreed well with our experimental data. * grigorev@phystech.edu †

The evolution of orientation-dependent metastable states during shock-induced solid-liquid phase transitions in crystalline Al is followed using moving window molecular dynamics simulations. The orientation-dependent transition pathways towards an orientation-independent final state Hugoniot include both "cold melting" followed by recrystallization in [110]- and [111]-oriented shock waves and crystal overheating followed by melting in [100] shock waves. The orientation-dependent dynamics take place within a zone that can extend up to hundreds of nanometers behind the shock front.