The role of the photochemical fragmentation in laser ablation: a molecular dynamics study

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 × 10¹² 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.

Molecular dynamics investigation of desorption and ion separation following picosecond infrared laser (PIRL) ablation of an ionic aqueous protein solution

Molecular dynamics simulations were performed to characterize the ablation process induced by a picosecond infrared laser (PIRL) operating in the regime of desorption by impulsive vibrational excitation (DIVE) of a model peptide (lysozyme)/counter-ion system in aqueous solution. The simulations were performed for ablation under typical experimental conditions found within a time-of-flight mass spectrometer (TOF-MS), that is in vacuum with an applied electric field (E = ± 10⁷ V/m), for up to 2 ns post-ablation and compared to the standard PIRL-DIVE ablation condition (E = 0 V/m). Further, a simulation of ablation under an extreme field condition (E = 10¹⁰ V/m) was performed for comparison to extend the effective dynamic range of the effect of the field on charge separation. The results show that the plume dynamics were retained under a typical TOF-MS condition within the first 1 ns of ablation. Efficient desorption was observed with more than 90% of water molecules interacting with lysozyme stripped off within 1 ns post-ablation. The processes of ablation and desolvation of analytes were shown to be independent of the applied electric field and thus decoupled from the ion separation process. Unlike under the extreme field conditions, the electric field inside a typical TOF-MS was shown to modify the ions' motion over a longer time and in a soft manner with no enhancement to fragmentation observed as compared to the standard PIRL-DIVE. The study indicates that the PIRL-DIVE ablation mechanism could be used as a new, intrinsically versatile, and highly sensitive ion source for quantitative mass spectrometry.

A molecular dynamics study of the effects of the inclusion of dopants on ablation in polymethyl methacrylate

Molecular dynamics simulations are used to elucidate mechanisms of ablation in dopant-polymer systems. In one set of simulations, a uniform distribution of thermal absorbers are added to a polymethyl methacrylate substrate and are excited. Chemical decomposition occurs in the regions near the absorbers. Ejection of large pieces of substrate then follows the thermo-chemical breakdown of material. In another set of simulations, an absorbing cluster is embedded in the polymethyl methacrylate substrate at a depth of 50 or 250 A. Only the particles comprising the cluster are excited during the laser pulse. Ejection of material is initiated upon the fracture of the cluster and the cleavage of the surrounding polymer bonds with little chemical damage during the process. These two mechanisms of ejection suggest different pathways of ablation in doped polymer materials.

Coupled molecular dynamics-Monte Carlo model to study the role of chemical processes during laser ablation of polymeric materials

The coarse grained chemical reaction model is enhanced to build a molecular dynamics (MD) simulation framework with an embedded Monte Carlo (MC) based reaction scheme. The MC scheme utilizes predetermined reaction chemistry, energetics, and rate kinetics of materials to incorporate chemical reactions occurring in a substrate into the MD simulation. The kinetics information is utilized to set the probabilities for the types of reactions to perform based on radical survival times and reaction rates. Implementing a reaction involves changing the reactants species types which alters their interaction potentials and thus produces the required energy change. We discuss the application of this method to study the initiation of ultraviolet laser ablation in poly(methyl methacrylate). The use of this scheme enables the modeling of all possible photoexcitation pathways in the polymer. It also permits a direct study of the role of thermal, mechanical, and chemical processes that can set off ablation. We demonstrate that the role of laser induced heating, thermomechanical stresses, pressure wave formation and relaxation, and thermochemical decomposition of the polymer substrate can be investigated directly by suitably choosing the potential energy and chemical reaction energy landscape. The results highlight the usefulness of such a modeling approach by showing that various processes in polymer ablation are intricately linked leading to the transformation of the substrate and its ejection. The method, in principle, can be utilized to study systems where chemical reactions are expected to play a dominant role or interact strongly with other physical processes.

Prediction of solubility parameters using partial least square regression

The total solubility parameter (delta) values were effectively predicted by using computed molecular descriptors and multivariate partial least squares (PLS) statistics. The molecular descriptors in the derived models included heat of formation, dipole moment, molar refractivity, solvent-accessible surface area (SA), surface-bounded molecular volume (SV), unsaturated index (Ui), and hydrophilic index (Hy). The values of these descriptors were computed by the use of HyperChem 7.5, QSPR Properties module in HyperChem 7.5, and Dragon Web version. The other two descriptors, hydrogen bonding donor (HD), and hydrogen bond-forming ability (HB) were also included in the models. The final reduced model of the whole data set had R(2) of 0.853, Q(2) of 0.813, root mean squared error from the cross-validation of the training set (RMSEcv(tr)) of 2.096 and RMSE of calibration (RMSE(tr)) of 1.857. No outlier was observed from this data set of 51 diverse compounds. Additionally, the predictive power of the developed model was comparable to the well recognized systems of Hansen, van Krevelen and Hoftyzer, and Hoy.

Ion formation mechanisms in UV-MALDI

Matrix Assisted Laser Desorption/Ionization (MALDI) is a very widely used analytical method, but has been developed in a highly empirical manner. Deeper understanding of ionization mechanisms could help to design better methods and improve interpretation of mass spectra. This review summarizes current mechanistic thinking, with emphasis on the most common MALDI variant using ultraviolet laser excitation. A two-step framework is gaining acceptance as a useful model for many MALDI experiments. The steps are primary ionization during or shortly after the laser pulse, followed by secondary reactions in the expanding plume of desorbed material. Primary ionization in UV-MALDI remains somewhat controversial, the two main approaches are the cluster and pooling/photoionization models. Secondary events are less contentious, ion-molecule reaction thermodynamics and kinetics are often invoked, but details differ. To the extent that local thermal equilibrium is approached in the plume, the mass spectra may be straightforwardly interpreted in terms of charge transfer thermodynamics.

Molecular Dynamics Model of Ultraviolet Matrix-Assisted Laser Desorption/Ionization Including Ionization Processes

A molecular dynamics model of UV-MALDI including ionization processes is presented. In addition to the previously described breathing sphere approach developed for simulation of laser ablation/desorption of molecular systems, it includes radiative and nonradiative decay, exciton hopping, two pooling processes, and electron capture. The results confirm the main conclusions of the continuum model of Knochenmuss, Anal. Chem. 2003, 75, 2199, but provide a much more detailed description of the interaction between ablation/desorption and ionization processes in the critical early time regime. Both desorption and ablation regimes generate free ions, and yields are in accordance with experiment. The first molecular ions are emitted at high velocities shortly before neutral desorption begins, because of surface charging caused by electron escape from the top of the sample. Later ions are entrained and thermalized in the plume of neutral molecules and clusters. Clusters are found to be stable on a nanosecond time scale, so the ions in them will be released only slowly, if at all. Exciton hopping rate and the mean radius for ion recombination are shown to be key parameters that can have a significant effect on net ion yield.

Coarse-Grained Model of the Interaction of Light with Polymeric Material: Onset of Ablation

A coarse-grained model has been developed for molecular dynamics simulations of the interaction of light with polymeric materials. The photon energy can result in a vibrational excitation (photothermal process) or disruption of a chemical bond (photochemical process) in a polymer. In the latter case, the formation of active radial sites and the occurrence of chemical reactions have to be taken into consideration. The novel feature of this model is the incorporation of chemical reactions into the united atom approximate representation of the polymer structure, which permits the study of laser ablation, degradation, or the effect of various chemical reactions on large time and length scales. The chemical reactions are included in the model in a probabilistic manner as in the kinetic Monte Carlo method. This model adopts physically and experimentally known quantities such as enthalpies and probabilities of reactions. Properties such as laser irradiation time, laser fluence, and wavelength are explicitly included. Moreover, no chemically correct interaction potential is required to incorporate the effects of chemical reactions on the dynamics of the system after energy deposition. We find that the model provides a plausible description of the essential processes. The laser-induced pressure relaxation is the main mechanism responsible for the onset of polymer ablation. Since the pressure relaxation processes are slow, there is a delay in the onset of ablation after the end of the laser pulse as is observed experimentally. The vaporization processes are not efficient for material removal, and the effect is minimal for both photochemical and photothermal processes. A lower fluence is needed for the onset of ablation with photochemical processes than photothermal processes.