Monte Carlo simulation study of the effects of nonequilibrium chemical reactions during pulsed laser desorption

Planar Gas Expansion under Intensive Nanosecond Laser Evaporation into Vacuum as Applied to Time-of-Flight Analysis

A computational investigation of the dynamics of gas expansion due to intense nanosecond laser evaporation into vacuum has been carried out. The problem is solved in a one-dimensional approximation, which simplifies calculations and at the same time allows one to analyze the main features of the expansion dynamics. For analysis we use three different approaches. Two of them are based on kinetic analysis via the direct simulation Monte Carlo (DSMC) method and numerical solution of the model Bhatnagar-Gross-Krook (BGK) equation. The third one focuses on derivation of an analytical continuum solution. Emphasis is placed on the analysis of the velocity distribution function and the average energy of particles passing through the time-of-flight detector on the normal to the evaporation surface, which is important for interpreting experimental measurements. The formulated problem is quite difficult as the considered flow is time-dependent, contains discontinuities in boundary conditions and involves large variations of local Knudsen numbers as well as steep gradients of the velocity distribution function. Data were obtained on the particle energy in the time-of-flight distribution for the range of regimes from the free molecular flow to continuum one. The maximum attainable average energy of particles in the time-of-flight distribution is determined. The non-monotonicity of the energy increase was found, which is explained based on analysis of the velocity distribution of particles.

Interpretation of time-of-flight distributions for neutral particles under pulsed laser evaporation using direct Monte Carlo simulation

A theoretical study of the time-of-flight (TOF) distributions under pulsed laser evaporation in vacuum has been performed. A database of TOF distributions has been calculated by the direct simulation Monte Carlo (DSMC) method. It is shown that describing experimental TOF signals through the use of the calculated TOF database combined with a simple analysis of evaporation allows determining the irradiated surface temperature and the rate of evaporation. Analysis of experimental TOF distributions under laser ablation of niobium, copper, and graphite has been performed, with the evaluated surface temperature being well agreed with results of the thermal model calculations. General empirical dependences are proposed, which allow indentifying the regime of the laser induced thermal ablation from the TOF distributions for neutral particles without invoking the DSMC-calculated database.

Computer Simulations of Laser Ablation, Plume Expansion and Plasma Formation

An overview is given of different modeling work that has been carried out, and is currently going on in our research group, in the field of modeling for laser ablation (LA). Most emphasis will be put on nanosecond (ns) LA, more specifically describing the laser-solid interaction, leading to heating, melting and vaporization of the target, by a heat conduction model, the expansion of the evaporated plume in vacuum or in a background gas by a set of conservation equations, and the plasma formation in the plume, assuming local thermal equilibrium. Some first results for nanoparticle formation in the expanding plume will be presented as well. Also, the process of target heating in the case of femtosecond (fs) LA will be described by means of a two-temperature model, and phase transitions, more specifically evaporation, will be illustrated by means of molecular dynamics simulations.

Direction-selective free expansion of laser-produced plasmas from planar targets

Direction-selective expansion of laser-produced plasmas from planar slab targets of Al, Ni, Mo, and Ta are reported. Angular distributions of the particles emitted from the targets, produced by a 130 mJ, 5 nsec, Nd:YAG laser, were obtained by means of a retarding potential analyzer and a quartz crystal. It was observed that the angular distributions of the particles had mainly three characteristics. For a given laser energy and a given target element, the angular distribution showed more preferential focusing toward the target normal as the value of the focal spot size B increased. Second, for a given laser energy and a given focal spot size, the focusing was more pronounced toward the target normal as the atomic mass number of the target materials increased. Third, for a given energy, a given focal spot size and a given element, the particles with higher ionization states were much more focused toward the target normal. Our experimental results confirm the Monte Carlo simulation results of the earlier works taking into account collision and recombination processes.

Laser-generated plasma plume expansion: combined continuous-microscopic modeling

The physical phenomena involved in the interaction of a laser-generated plasma plume with a background gas are studied numerically. A three-dimensional combined model is developed to describe the plasma plume formation and its expansion in vacuum or into a background gas. The proposed approach takes advantages of both continuous and microscopic descriptions. The simulation technique is suitable for the simulation of high-rate laser ablation for a wide range of background pressure. The model takes into account the mass diffusion and the energy exchange between the ablated and background species, as well as the collective motion of the ablated species and the background-gas particles. The developed approach is used to investigate the influence of the background gas on the expansion dynamics of the plume obtained during the laser ablation of aluminum. At moderate pressures, both plume and gas compressions are weak and the process is mainly governed by the diffusive mixing. At higher pressures, the interaction is determined by the plume-gas pressure interplay, the plume front is strongly compressed, and its center exhibits oscillations. In this case, the snowplough effect takes place, leading to the formation of a compressed gas layer in front of the plume. The background pressure needed for the beginning of the snowplough effect is determined from the plume and gas density profiles obtained at various pressures. Simulation results are compared with experimentally measured density distributions. It is shown that the calculations suggest localized formation of molecules during reactive laser ablation.