Molecular dynamics simulation of the laser disintegration of aerosol particles

Physical Regimes and Mechanisms of Picosecond Laser Fragmentation of Gold Nanoparticles in Water from X-ray Probing and Atomistic Simulations

Laser fragmentation in liquids has emerged as a promising green chemistry technique for changing the size, shape, structure, and phase composition of colloidal nanoparticles, thus tuning their properties to the needs of practical applications. The advancement of this technique requires a solid understanding of the mechanisms of laser-nanoparticle interactions that lead to the fragmentation. While theoretical studies have made impressive practical and mechanistic predictions, their experimental validation is required. Hence, using the picosecond laser fragmentation of Au nanoparticles in water as a model system, the transient melting and fragmentation processes are investigated with a combination of time-resolved X-ray probing and atomistic simulations. The direct comparison of the diffraction profiles predicted in the simulations and measured in experiments has revealed a sequence of several nonequilibrium processes triggered by the laser irradiation. At low laser fluences, in the regime of nanoparticle melting and resolidification, the results provide evidence of a transient superheating of crystalline nanoparticles above the melting temperature. At fluences about three times the melting threshold, the fragmentation starts with evaporation of Au atoms and their condensation into small satellite nanoparticles. As fluence increases above five times the melting threshold, a transition to a rapid (explosive) phase decomposition of superheated nanoparticles into small liquid droplets and vapor phase atoms is observed. The transition to the phase explosion fragmentation regime is signified by prominent changes in the small-angle X-ray scattering profiles measured in experiments and calculated in simulations. The good match between the experimental and computational diffraction profiles gives credence to the physical picture of the cascade of thermal fragmentation regimes revealed in the simulations and demonstrates the high promise of the joint tightly integrated computational and experimental efforts.

Molecular Modeling in Anion Exchange Membrane Research: A Brief Review of Recent Applications

Anion Exchange Membrane (AEM) fuel cells have attracted growing interest, due to their encouraging advantages, including high power density and relatively low cost. AEM is a polymer matrix, which conducts hydroxide (OH-) ions, prevents physical contact of electrodes, and has positively charged head groups (mainly quaternary ammonium (QA) groups), covalently bound to the polymer backbone. The chemical instability of the quaternary ammonium (QA)-based head groups, at alkaline pH and elevated temperature, is a significant threshold in AEMFC technology. This review work aims to introduce recent studies on the chemical stability of various QA-based head groups and transportation of OH- ions in AEMFC, via modeling and simulation techniques, at different scales. It starts by introducing the fundamental theories behind AEM-based fuel-cell technology. In the main body of this review, we present selected computational studies that deal with the effects of various parameters on AEMs, via a variety of multi-length and multi-time-scale modeling and simulation methods. Such methods include electronic structure calculations via the quantum Density Functional Theory (DFT), ab initio, classical all-atom Molecular Dynamics (MD) simulations, and coarse-grained MD simulations. The explored processing and structural parameters include temperature, hydration levels, several QA-based head groups, various types of QA-based head groups and backbones, etc. Nowadays, many methods and software packages for molecular and materials modeling are available. Applications of such methods may help to understand the transportation mechanisms of OH- ions, the chemical stability of functional head groups, and many other relevant properties, leading to a performance-based molecular and structure design as well as, ultimately, improved AEM-based fuel cell performances. This contribution aims to introduce those molecular modeling methods and their recent applications to the AEM-based fuel cells research community.

A GPU-Accelerated Fast Multipole Method for GROMACS: Performance and Accuracy

An important and computationally demanding part of molecular dynamics simulations is the calculation of long-range electrostatic interactions. Today, the prevalent method to compute these interactions is particle mesh Ewald (PME). The PME implementation in the GROMACS molecular dynamics package is extremely fast on individual GPU nodes. However, for large scale multinode parallel simulations, PME becomes the main scaling bottleneck as it requires all-to-all communication between the nodes; as a consequence, the number of exchanged messages scales quadratically with the number of involved nodes in that communication step. To enable efficient and scalable biomolecular simulations on future exascale supercomputers, clearly a method with a better scaling property is required. The fast multipole method (FMM) is such a method. As a first step on the path to exascale, we have implemented a performance-optimized, highly efficient GPU FMM and integrated it into GROMACS as an alternative to PME. For a fair performance comparison between FMM and PME, we first assessed the accuracies of the methods for various sets of input parameters. With parameters yielding similar accuracies for both methods, we determined the performance of GROMACS with FMM and compared it to PME for exemplary benchmark systems. We found that FMM with a multipole order of 8 yields electrostatic forces that are as accurate as PME with standard parameters. Further, for typical mixed-precision simulation settings, FMM does not lead to an increased energy drift with multipole orders of 8 or larger. Whereas an ≈50 000 atom simulation system with our FMM reaches only about a third of the performance with PME, for systems with large dimensions and inhomogeneous particle distribution, e.g., aerosol systems with water droplets floating in a vacuum, FMM substantially outperforms PME already on a single node.

Finite element simulation of infrared laser ablation for mass spectrometry

RATIONALE

Laser ablation is widely used in conjunction with ambient ionization techniques, and a fundamental understanding of the mechanism of material removal is important to its optimal use in mass spectrometry. Finite element analysis simulates the laser material interaction on larger time and distance scales than atomistic approaches. Here, a two-dimensional finite element model was developed to simulate infrared laser irradiation of glycerol using a wavelength-tunable infrared (IR) laser.

METHODS

The laser fluence used for the simulations was varied from 1000 to 6000 J/m(2), the wavelength was varied from 2.7 to 3.7 µm, and both flat-top and Gaussian shape laser profiles were studied.

RESULTS

Phase explosion conditions were found for laser wavelengths near 3 µm (which corresponds to the OH stretch absorption of glycerol) and fluences above 2000 J/m(2). This suggests that laser ablation of glycerol is driven by phase explosion in the OH stretch region. The Gaussian profile generated regions of higher glycerol temperature, whereas the flat-top profile heated a larger volume of material above the phase explosion temperature.

CONCLUSIONS

These results suggest that the best performance for pulsed IR laser sample irradiation is in the wavelength range from 2.9 to 3.1 µm for materials with a strong OH stretch absorption.

Ion formation mechanism in laser desorption ionization of individual nanoparticles

Covariance mapping is used to study ion formation mechanisms in laser desorption ionization of individual 50 or 220 nm diameter particles having compositions similar to ambient aerosol. Single particle mass spectra are found to vary substantially from particle to particle. This variation is systematic--the energetically preferred ions (e.g., lowest ionization energy, highest electron affinity) are positively correlated with each other and negatively correlated with less preferred ions. For the compositions studied, the average positive ion yield is two to five times greater than the negative ion yield, indicating that free electrons are the main negatively charged species. For many particles, typically 20% to 40% of those analyzed, only positive ions are detected. Smaller particles give fewer negative ions, presumably because the plume is less dense and electron capture is less likely. The results suggest that ion formation occurs by a two stage process. In the first stage, photoionization of laser desorbed neutrals gives cations and free electrons. In the second stage, collisions in the plume cause electron capture and competitive charge transfer. When the particle ablates in a manner giving a dense plume with many collisions, the energetically preferred positive and negative ions are dominant. When the particle ablates in a manner giving a less dense plume with fewer collisions, the less preferred ions are able to survive and the energetically preferred ions constitute a lower fraction of the total ion signal. Systematic particle to particle variations of relative signal intensities can complicate ambient particle classification efforts by spreading a single particle composition over several classes.

The design of single particle laser mass spectrometers

This review explores some of the design choices made with single particle mass spectrometers. Different instruments have used various configurations of inlets, particle sizing techniques, ionization lasers, mass spectrometers, and other components. Systematic bias against non-spherical particles probably exceeds a factor of 2 for all instruments. An ionization laser tradeoff is the relatively poor beam quality and reliability of an excimer laser versus the longer wavelengths and slower response time of an Nd-YAG laser. Single particle instruments can make special demands on the speed and dynamic range of the mass spectrometers. This review explains some of the choices made for instruments that were developed for different types of measurements in the atmosphere. Some practical design notes are also given from the author's experience with each section of the instrument.

Time-Resolved Micro Liquid Desorption Mass Spectrometry: Mechanism, Features, and Kinetic Applications

Liquid water beam desorption mass spectrometry is an intriguing technique to isolate charged molecular aggregates directly from the liquid phase and to analyze them employing sensitive mass spectrometry. The liquid phase in this approach consists of a 10 µm diameter free liquid filament in vacuum which is irradiated by a focussed infrared laser pulse resonant with the OH-stretch vibration of bulk water. Depending upon the laser wavelength, charged (e.g. protonated) macromolecules are isolated from solution through a still poorly characterized mechanism. After the gentle liquid-to-vacuum transfer the low-charge-state aggregates are analyzed using time-of-flight mass spectrometry. A recent variant of the technique uses high performance liquid chromatography valves for local liquid injections of samples in the liquid carrier beam, which enables very low sample consumption and high speed sample analysis. In this review we summarize recent work to characterize the ‘desorption’ or ion isolation mechanism in this type of experiment. A decisive and interesting feature of micro liquid beam desorption mass spectrometry is that — under certain conditions — the gas-phase mass signal for a large number of small as well as supramolecular systems displays a surprisingly linear response on the solution concentration over many orders of magnitude, even for mixtures and complex body fluids. This feature and the all-liquid state nature of the technique makes this technique a solution-type spectroscopy that enables real kinetic studies involving (bio)polymers in solution without the need for internal standards. Two applications of the technique monitoring enzyme digestion of proteins and protein aggregation of an amyloid model system are highlighted, both displaying its potential for monitoring biokinetics in solution.

Photochemical Interaction of Polystyrene Nanospheres with 193 nm Pulsed Laser Light

The photochemical interaction of 193 nm light with polystyrene nanospheres is used to produce particles with a controlled size and morphology. Laser fluences from 0 to 0.14 J/cm2 at 10 and 50 Hz photofragment nearly monodisperse 110 nm spherical polystyrene particles. The size distributions before and after irradiation are measured with a scanning mobility particle sizer (SMPS), and the morphology of the irradiated particles is examined with a transmission electron microscope (TEM). The results show that the irradiated particles have a smaller mean diameter ( approximately 25 nm) and a number concentration more than an order of magnitude higher than nonirradiated particles. The particles are formed by nucleation of gas-phase species produced by photolytic decomposition of nanospheres. A nondimensional parameter, the photon-to-atom ratio (PAR), is used to interpret the laser-particle interaction energetics.

Improved sensitivity and mass range in time-of-flight bioaerosol mass spectrometry using an electrostatic ion guide

Bioearosol mass spectrometry (BAMS) analyzes single particles in real time from ambient air, placing strict demands on instrument sensitivity. Modeling of the BAMS reflectron time of flight (TOF) with SIMION revealed design limitations associated with ion transmission and instrument sensitivity at higher masses. Design and implementation of a BAMS linear TOF with electrostatic ion guide and delayed extraction capabilities has greatly increased the sensitivity and mass range relative to the reflectron design. Initial experimental assessment of the new instrument design revealed improved sensitivity at high masses as illustrated when using standard particles of cytochrome C (m/z approximately 12,000), from which the compound's monomer, dimer (m/z approximately 24,000) and trimer (m/z approximately 36,000) were readily detected.

Calibration Effects for Laser-Induced Breakdown Spectroscopy of Gaseous Sample Streams: Analyte Response of Gas-Phase Species versus Solid-Phase Species

The effects of analyte phase on the calibration response for laser-induced breakdown spectroscopy is investigated for a range of carbon species. Significant differences in the atomic emission signal from carbon were observed when comparing calibration streams of gas-phase and submicrometer-sized solid-phase carbon species. The resulting calibration curve slopes varied by a factor of 8 over a comparable range of atomic carbon concentrations for five different analyte sources, while the plasma electron density and temperature remained essentially constant. The current findings challenge a widely held assumption that complete dissociation of constituent species within a highly energetic laser-induced plasma results in independence of the analyte atomic emission signal on the analyte source. A physical model of the plasma-analyte interaction is proposed that provides a framework to account for the observed dependence on the physical state of the analyte.

Initial velocity distributions of ions generated by in-flight laser desorption/ionization of individual polystyrene latex microparticles as studied by the delayed ion extraction method

The delayed ion extraction method has been used to study characteristics of the initial velocity distributions of positive and negative ions produced simultaneously by laser desorption/ionization (LDI) from non-impacted single aerosol polymeric particles, using a bipolar time-of-flight (TOF) instrument (LAMPAS 2). Due to the geometry of the setup and the characteristics of the ablation process, only the projections of the velocities on the axis of the mass spectrometer can be directly studied. Additionally, since the mean initial velocity under these conditions should be close to zero, it was necessary to extend the method by taking into account higher order contributions of the velocity distribution. Theoretical expressions for these higher order terms are presented and discussed. The bipolar characteristics of the instrument permit evaluation and treatment of a possible instrumental artifact caused by small inclinations of the ionizing laser with respect to the ideal incidence direction. Results of a number of experiments are presented and discussed in relation to the theoretical expressions presented, and to possible ablation scenarios. Evidence pointing out that, under our experimental conditions, only partial ablation of the latex particles occurs was obtained. The variance of the distribution of the projection of the initial velocities can be directly estimated from these results. By assuming that the total initial velocities of the ions are developed completely according to a single-temperature adiabatic expansion mechanism, temperatures of approximately 50 K/Da can be assigned to the ion clouds from the variance estimations. If a two-temperature model is used, a radial temperature of about 100 K/Da results. These values are in reasonable agreement with results for polymer ablation from the literature.

Time-resolved studies of the interactions between pulsed lasers and aerosols

Studies of the interaction between a pulsed CO2 laser and micrometer-sized aqueous and organic particles by use of light-scattering methods and step-scan Fourier-transform infrared (FTIR) spectroscopy are reported. Visible two-color extinction experiments indicate primary particle shattering, accompanied by a high fraction of vaporization, followed by secondary particle evaporation. The extent of the latter depends on the pulse intensity and particle composition. Angle-resolved light-scattering investigations provide insight into the aerosol size distribution and temperature following the pulsed heating event. The time dependence of the vapor plume, monitored with step-scan FTIR spectroscopy, confirms that a large fraction of the initial particle is quickly evaporated during the shattering event, followed by secondary fragment evaporation and thermal expansion.

Depth profiling of heterogeneously mixed aerosol particles using single-particle mass spectrometry

Infrared laser evaporation of single aerosol particles in a vacuum followed by vacuum ultraviolet (VUV) laser ionization and time-of-flight mass spectroscopy of the resulting vapor provides a depth profile of the particle's composition. Analyzing glycerol particles coated with 60-150-nm coatings of oleic acid using either a CO2 laser or a tunable optical parametric oscillator as an evaporation laser results in mass spectra that depend on the IR laser power. Low infrared laser powers incompletely vaporize particles and preferentially probe the composition of the surface layers of the particle, but high laser powers evaporate the entire particle and produce spectra representative of the particle's total composition. In the limit of low laser power, the fraction of oleic acid in the mass spectra is as much as 50 times greater than the fraction of oleic acid in the particle, providing a surface-layer-specific characterization. The OPO laser provides even more surface specificity, producing an [oleic acid]/[glycerol] ratio as much as four times larger (for a 60-nm coating) than that obtained using the CO2 laser. The infrared laser power required to sample the core of the particle increases with the thickness of the coating and is sensitive to changes in the coating thickness on the order of 10 nm. In contrast to these intuitively appealing results, high CO2 laser powers (approximately 90 mJ/pulse) produce mass spectra that, at short delays between the CO2 and VUV lasers, show enrichment of the core material rather than the coating. Likewise, tuning the OPO to frequencies that are resonant with the core material but transparent to the coating also results in selective detection of the core. The results suggest that a shattering mechanism dominates the vaporization dynamics in these situations.

Enhancing the detection of sulfate particles for laser ablation aerosol mass spectrometry

A major limitation to the application of laser ablation aerosol mass spectrometry for the detection of particles less than 200 nm in diameter is a low ablation efficiency for sulfate particles. (Ablation efficiency is the probability that an ablated particle produces a detectable ion signal.) A method is described here to enhance the ablation efficiency of sulfate particles by coating them with a UV-absorbing compound. The method can be applied in-line with the aerosol mass spectrometer in a manner that does not significantly alter the aerosol size distribution. It is shown that a 12-nm coating of 1-naphthyl acetate increases the ablation efficiency of 136-nm ammonium sulfate particles by at least a factor of 20, while similar coatings on oleic acid and ammonium nitrate particles do not significantly alter the ablation efficiency. The results suggest that "undetected" particles, presumably sulfate, in ambient aerosol can be assessed.