We use a combination of tandem ion mobility spectrometry (IMS–IMS, with differential mobility analyzers), molecular dynamics (MD) simulations, and analytical models to examine both neutral solvent (H2O) and ion (solvated Na+) evaporation from aqueous sodium chloride nanodrops. For experiments, nanodrops were produced via electrospray ionization (ESI) of an aqueous sodium chloride solution. Two nanodrops were examined in MD simulations: a 2500 water molecule nanodrop with 68 Na+ and 60 Cl− ions (an initial net charge of z = +8), and (2) a 1000 water molecule nanodrop with 65 Na+ and 60 Cl− ions (an initial net charge of z = +5). Specifically, we used MD simulations to examine the validity of a model for the neutral evaporation rate incorporating both the Kelvin (surface curvature) and Thomson (electrostatic) influences, while both MD simulations and experimental measurements were compared to predictions of the ion evaporation rate equation of Labowsky et al. [Anal. Chim. Acta, 2000, 406, 105–118]. Within a single fit parameter, we find excellent agreement between simulated and modeled neutral evaporation rates for nanodrops with solute volume fractions below 0.30. Similarly, MD simulation inferred ion evaporation rates are in excellent agreement with predictions based on the Labowsky et al. equation. Measurements of the sizes and charge states of ESI generated NaCl clusters suggest that the charge states of these clusters are governed by ion evaporation, however, ion evaporation appears to have occurred with lower activation energies in experiments than was anticipated based on analytical calculations as well as MD simulations. Several possible reasons for this discrepancy are discussed.
Accurate calculation of the ion–ion recombination rate coefficient has been of long-standing interest as it controls the ion concentration in gas phase systems and in aerosols. We describe the development of a hybrid continuum-molecular dynamics (MD) approach to determine the ion–ion recombination rate coefficient. This approach is based on the limiting sphere method classically used for transition regime collision phenomena in aerosols. When ions are sufficiently far from one another, the ion–ion relative motion is described by diffusion equations, while within a critical distance, MD simulations are used to model ion–ion motion. MD simulations are parameterized using the Assisted Model Building with Energy Refinement force-field as well as by considering partial charges on atoms. Ion–neutral gas collisions are modeled in two mutually exclusive cubic domains composed of 103 gas atoms each, which remain centered on the recombining ions throughout calculations. Example calculations are reported for NH4+ recombination with NO2− in He, across a pressure range from 10 kPa to 10 000 kPa. Excellent agreement is found in comparison with calculations to literature values for the 100 kPa recombination rate coefficient (1.0 × 10−12 m3 s−1) in He. We also recover the experimentally observed increase in the recombination rate coefficient with pressure at sub-atmospheric pressures, and the observed decrease in the recombination rate coefficient in the high pressure continuum limit. We additionally find that non-dimensionalized forms of rate coefficients are consistent with recently developed equations for the dimensionless charged particle–ion collision rate coefficient based on Langevin dynamics simulations.
We conducted intratracheal instillations of polyacrylic acid (PAA) with crosslinking and non-crosslinking into rats in order to examine what kinds of physicochemical characteristics of acrylic-acid-based polymers affect responses in the lung. F344 rats were intratracheally exposed to similar molecular weights of crosslinked PAA (CL-PAA) (degree of crosslinking: ~0.1%) and non-crosslinked PAA (Non-CL-PAA) at low and high doses. Rats were sacrificed at 3 days, 1 week, 1 month, 3 months, and 6 months post-exposure. Both PAAs caused increases in neutrophil influx, cytokine-induced neutrophil chemoattractants (CINC) in the bronchoalveolar lavage fluid (BALF), and heme oxygenase-1 (HO-1) in the lung tissue from 3 days to 6 months following instillation. The release of lactate dehydrogenase (LDH) activity in the BALF was higher in the CL-PAA-exposed groups. Histopathological findings of the lungs demonstrated that the extensive fibrotic changes caused by CL-PAA were also greater than those in exposure to the Non-CL- PAA during the observation period. CL-PAA has more fibrogenicity of the lung, suggesting that crosslinking may be one of the physicochemical characteristic factors of PAA-induced lung disorder.
Background: We conducted intratracheal instillations of different molecular weights of polyacrylic acid (PAA) into rats in order to examine what kinds of physicochemical characteristics of acrylic acid-based polymer affect responses in the lung. Methods: F344 rats were intratracheally exposed to a high molecular weight (HMW) of 598 thousand g/mol or a low molecular weight (LMW) of 30.9 thousand g/mol PAA at low and high doses. Rats were sacrificed at 3 days, 1 week, 1 month, 3 months and 6 months post exposure. Results: HMW PAA caused persistent increases in neutrophil influx, cytokine-induced neutrophil chemoattractants (CINC) in the bronchoalveolar lavage fluid (BALF), and heme oxygenase-1 (HO-1) in the lung tissue from 3 days to 3 months and 6 months following instillation. On the other hand, LMW PAA caused only transient increases in neutrophil influx, CINC in BALF, and HO-1 in the lung tissue from 3 days to up to 1 week or 1 month following instillation. Histopathological findings of the lungs demonstrated that the extensive inflammation and fibrotic changes caused by the HMW PAA was greater than that in exposure to the LMW PAA during the observation period. Conclusion: HMW PAA induced persistence of lung disorder, suggesting that molecular weight is a physicochemical characteristic of PAA-induced lung disorder.
We performed a molecular dynamics simulation to calculate the self-diffusion coefficients of 1-Butyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide and water in a water–ionic liquid mixture. We then compared the simulated self-diffusion coefficients of cation, anion and water molecules with experimental data and with simulated data from the literature. Although the simulation overestimated the self-diffusion coefficients of ions, the simulated results qualitatively reproduced the enhancement of the self-diffusion coefficients of water as the water molar fraction increased. We also calculated the radial distribution functions to investigate the solution structure, i.e. the clustering of water molecules. The clustering of water in ionic liquid was found to play an important role in the enhancement of the diffusion of water molecules in the ionic liquid.