Scattering of Ne from the liquid–vapor interface of glycerol: A molecular dynamics study

Energy- and Angle-Resolved Scattering of Ne from Dodecane Liquid Surfaces: Theory Corroborating Experiment

Motivated by recent experimental work by the Neumark group, we present here an all-atom molecular dynamics study of Ne scattering from a dodecane liquid surface with the objective of elucidating the fundamental aspects of gas-liquid dynamics. Using a fine-tuned force field, the GPU-accelerated simulations reproduced semiquantitatively the energy- and angle-resolved experimental results. The branching ratio between the impulsive scattering (IS) and thermal desorption (TD) channels exhibits a clear correlation with the incidence energy (Ei) and angle. Ne atoms with lower Ei values are more likely to be trapped, yielding an increased TD ratio. For a given Ei, a large incidence angle led to a higher IS ratio. The energy transfer between Ne atoms and liquid dodecane was found to be more sensitive to the deflection angle than to the incidence or reflection angle. With an increasing deflection angle, the fractional energy loss increases, suggesting that more kinetic energy is transferred to the liquid.

Competing Segregation of Br- and Cl- to a Surface Coated with a Cationic Surfactant: Direct Measurements of Ion and Solvent Depth Profiles

Ion-surface scattering experiments can be used to measure elemental depth profiles on the angstrom scale in complex liquid mixtures. We employ NICISS (neutral impact collision ion scattering spectroscopy) to measure depth profiles of dissolved ions and solvent in liquid glycerol containing the cationic surfactant tetrahexylammonium bromide (THA⁺/Br^-) at 0.013 M and mixtures of NaBr + NaCl at 0.4 M total concentration. The experiments reveal that Br^- outcompetes Cl^- in its attraction to surface THA⁺, and that THA⁺ segregates more extensively when more Br^- ions are present. Intriguingly, the depths spanned by THA⁺, Br^-, and Cl^- ions generally increase with Br^- bulk concentration, expanding from ∼10 to ∼25 Å for both Br^- and Cl^- depth profiles. This broadening likely occurs because of an increasing pileup of THA⁺ ions in a multilayer region that spreads the halide ions over a wider depth. The experiments indicate that cationic surfactants enhance Br^- and Cl^- concentrations in the surface region far beyond their bulk-phase values, making solutions coated with these surfactants potentially more reactive toward gases that can oxidize the halide ions.

Experimental Depth Profiles of Surfactants, Ions, and Solvent at the Angstrom Scale: Studies of Cationic and Anionic Surfactants and Their Salting Out

Neutral impact ion scattering spectroscopy (NICISS) is used to measure the depth profiles of ionic surfactants, counterions, and solvent molecules on the angstrom scale. The chosen surfactants are 0.010 m tetrahexylammonium bromide (THA⁺/Br^-) and 0.0050 m sodium dodecyl sulfate (Na⁺/DS^-) in the absence and presence of 0.30 m NaBr in liquid glycerol. NICISS determines the depth profiles of the elements C, O, Na, S, and Br through the loss in energy of 5 keV He atoms that travel into and out of the liquid, which is then converted into depth. In the absence of NaBr, we find that THA⁺ and its Br^- counterion segregate together because of charge attraction, forming a narrow double layer that is 10 Å wide and 150 times more concentrated than in the bulk. With the addition of NaBr, THA⁺ is "salted out" to the surface, increasing the interfacial Br^- concentration by 3-fold and spreading the anions over a ∼30 Å depth. Added NaBr similarly increases the interfacial concentration of DS^- ions and broadens their positions. Conversely, the dissolved Br^- ions are significantly depleted over a depth of 0-40 Å from the surface because of charge repulsion from DS^- ions within the interfacial region. These different interfacial Br^- propensities correlate with previously measured gas-liquid reactivities: gaseous Cl₂ readily reacts with Br^- ions in the presence of THA⁺ but drops 70-fold in the presence of DS^-, demonstrating that surfactant headgroup charge controls the reactivity of Br^- through changes in its depth profile.

Salty glycerol versus salty water surface organization: bromide and iodide surface propensities

Salty NaBr and NaI glycerol solution interfaces are examined in the OH stretching region using broadband vibrational sum frequency generation (VSFG) spectroscopy. Raman and infrared (IR) spectroscopy are used to further understand the VSFG spectroscopic signature. The VSFG spectra of salty glycerol solutions reveal that bromide and iodide anions perturb the interfacial glycerol organization in a manner similar as that found in aqueous halide salt solutions, thus confirming the presence of bromide and iodide anions at the glycerol surface. Surface tension measurements are consistent with the surface propensity suggested by the VSFG data and also show that the surface excess increases with increasing salt concentration, similar to that of water. In addition, iodide is shown to have more surface prevalence than bromide, as has also been determined from aqueous solutions. These results suggest that glycerol behaves similarly to water with respect to surface activity and solvation of halide anions at its air/liquid interface.

Molecular dynamics simulation study of glycerol-water liquid mixtures

To study the effects of water on conformational dynamics of polyalcohols, Molecular Dynamics simulations of glycerol-water liquid mixtures have been carried out at different concentrations: 42.9 and 60.0 wt % of glycerol, respectively. On the basis of the analysis of backbone conformer distributions, it is found that the surrounding water molecules have a large impact on the populations of the glycerol conformers. While the local structure of water in the liquid mixture is surprisingly close to that in pure liquid water, the behavior of glycerols can be divided into three different categories where roughly 25% of them occur in a structure similar to that in pure liquid of glycerol, ca. 25% of them exist as monomers, solvated by water, and the remaining 50% of glycerols in the mixture form H-bonded strings as remains of the glycerol H-bond network. The typical glycerol H-bond network still exists even at the lower concentration of 40 wt % of glycerol. The microheterogeneity of water-glycerol mixtures is analyzed using time-averaged distributions of the sizes of the water aggregates. At 40 wt % of glycerol, the cluster sizes from 3 to 10 water molecules are observed. The increase of glycerol content causes a depletion of clusters leading to smaller 3-5 molecule clusters domination. Translational diffusion coefficients have been calculated to study the dynamical behavior of both glycerol and water molecules. Rotational-reorientational motion is studied both in overall and in selected substructures on the basis of time correlation functions. Characteristic time scales for different motional modes are deduced on the basis of the calculated correlation times. The general conclusion is that the presence of water increases the overall mobility of glycerol, while glycerol slows the mobility of water.

A stochastic, local mode study of neon-liquid surface collision dynamics

Equations of motion for a fast, light rare gas atom passing over a liquid surface are derived and used to infer the dynamics of neon collisions with squalane and perfluorinated polyether surfaces from experimental data. The equations incorporate the local mode model of a liquid surface via a stochastic process and explicitly account for impulsive collisional energy loss to the surface. The equations predict angular distributions for scattering of neon that are in good quantitative agreement with experimental data. Our key dynamical conclusions are that experimental angular distributions derive mainly from local mode surface topography rather than from structural features of individual surface molecules, and that the available data for these systems can be accounted for almost exclusively by single collisions between neon atoms and the liquid surface.

A study on orientation and absorption spectrum of interfacial molecules by using continuum model

In this work, a numerical procedure based on the continuum model is developed and applied to the solvation energy for ground state and the spectral shift against the position and the orientation of the interfacial molecule. The interface is described as a sharp boundary separating two bulk media. The polarizable continuum model (PCM) allows us to account for both electrostatic and nonelectrostatic solute-solvent interactions when we calculate the solvation energy. In this work we extend PCM to the interfacial system and the information about the position and orientation of the interfacial molecule can be obtained. Based on the developed expression of the electrostatic free energy of a nonequilibrium state, the numerical procedure has been implemented and used to deal with a series of test molecules. The time-dependent density functional theory (TDDFT) associated with PCM is used for the electron structure and the spectroscopy calculations of the test molecules in homogeneous solvents. With the charge distribution of the ground and excited states, the position- and orientation-dependencies of the solvation energy and the spectrum have been investigated for the interfacial systems, taking the electrostatic interaction, the cavitation energy, and the dispersion-repulsion interaction into account. The cavitation energy is paid particular attention, since the interface portion cut off by the occupation of the interfacial molecule contributes an extra part to the stabilization for the interfacial system. The embedding depth, the favorable orientational angle, and the spectral shift for the interfacial molecule have been investigated in detail. From the solvation energy calculations, an explanation has been given on why the interfacial molecule, even if symmetrical in structure, tends to take a tilting manner, rather than perpendicular to the interface.

Collisions of DCl with Pure and Salty Glycerol: Enhancement of Interfacial D → H Exchange by Dissolved NaI

The fate of DCl molecules striking pure glycerol and a 2.6 M NaI-glycerol solution is investigated using scattering, uptake, and residence time measurements. We find that dissolved Na+ and I- ions alter every gas-liquid pathway from the moment of contact of DCl with the surface to its eventual emergence as HCl. In particular, the salt enhances both trapping-desorption of DCl and interfacial DCl --> HCl exchange at the expense of DCl entry into the bulk solution. The reduced entry and enhanced desorption of thermalized DCl molecules are interpreted by assuming that Na+ and I- ions bind to interfacial OH groups and tie up surface sites that would otherwise capture incoming DCl molecules. These ion-glycerol interactions may also be responsible for enhancing interfacial D --> H exchange by disrupting the interfacial hydrogen bond network that carries the newly formed H+ ion away from its Cl- pair. This disruption may increase the fraction of interfacial Cl- and H+ that recombine and desorb immediately as HCl before the ions separate and diffuse deeply into the bulk.

Molecular Dynamics Simulations of Atmospheric Oxidants at the Air−Water Interface: Solvation and Accommodation of OH and O3

A comparative study of OH, O3, and H2O equilibrium aqueous solvation and gas-phase accommodation on liquid water at 300 K is performed using a combination of ab initio calculations and molecular dynamics simulations. Polarizable force fields are developed for the interaction potential of OH and O3 with water. The free energy profiles for transfer of OH and O3 from the gas phase to the bulk liquid exhibit a pronounced minimum at the surface, but no barrier to solvation in the bulk liquid. The calculated surface excess of each oxidant is comparable to calculated and experimental values for short chain, aliphatic alcohols. Driving forces for the surface activity are discussed in terms of the radial distribution functions and dipole orientation distributions for each molecule in the bulk liquid and at the surface. Simulations of OH, O3, and H2O impinging on liquid water with a thermal impact velocity are used to calculate thermal accommodation (S) and mass accommodation (alpha) coefficients. The values of S for OH, O3, and H2O are 0.95, 0.90, and 0.99, respectively. The approaching molecules are accelerated toward the liquid surface when they are approximately 5 angstroms above it. The molecules that reach thermal equilibrium with the surface do so within 2 ps of striking the surface, while those that do not scatter into the gas phase with excess translational kinetic energy in the direction perpendicular to the surface. The time constants for absorption and desorption range from approximately 35 to 140 ps, and the values of alpha for OH, O3, and H2O are 0.83, 0.047, and 0.99, respectively. The results are consistent with previous formulations of gas-phase accommodation from simulations, in which the process occurs by rapid thermal and structural equilibration followed by diffusion on the free energy profile. The implications of these results with respect to atmospheric chemistry are discussed.

Classical trajectory study of collisions of Ar with alkanethiolate self-assembled monolayers: Potential-energy surface effects on dynamics

We have investigated collisions between Ar and alkanethiolate self-assembled monolayers (SAMs) using classical trajectory calculations with several potential-energy surfaces. The legitimacy of the potential-energy surfaces is established through comparison with molecular-beam data and ab initio calculations. Potential-energy surfaces used in previous work overestimate the binding of Ar to the SAM, leading to larger energy transfer than found in the experiments. New calculations, based on empirical force fields that better reproduce ab initio calculations, exhibit improved agreement with the experiments. In particular, polar-angle-dependent average energies calculated with explicit-atom potential-energy surfaces are in excellent agreement with the experiments. Polar- and azimuthal-angle-dependent product translational energies are examined to gain deeper insight into the dynamics of Ar+SAM collisions.

MOLECULAR BEAM STUDIES OF GAS-LIQUID INTERFACES

Molecular beam scattering experiments provide a way to disentangle the elementary steps involved in energy transfer and chemical reactions between gases and liquids. After surveying the history and recent progress in this field, we review studies of the kinematics of gas-liquid collisions and proton exchange of HCl, DCl, and HBr with supercooled sulfuric acid and liquid glycerol. These experiments help to clarify the role of the surface region in controlling trapping and interfacial- and bulk-phase reactions.

Adsorption and solvation of ethanol at the water liquid-vapor interface: a molecular dynamics study

The free energy profiles of methanol and ethanol at the water liquid-vapor interface at 310K were calculated using molecular dynamics computer simulations. Both alcohols exhibit a pronounced free energy minimum at the interface and, therefore, have positive adsorption at this interface. The surface excess was computed from the Gibbs adsorption isotherm and was found to be in good agreement with experimental results. Neither compound exhibits a free energy barrier between the bulk and the surface adsorbed state. Scattering calculations of ethanol molecules from a gas phase thermal distribution indicate that the mass accommodation coefficient is 0.98, and the molecules become thermalized within 10 ps of striking the interface. It was determined that the formation of the solvation structure around the ethanol molecule at the interface is not the rate-determining step in its uptake into water droplets. The motion of an ethanol molecule in a water lamella was followed for 30 ns. The time evolution of the probability distribution of finding an ethanol molecule that was initially located at the interface is very well described by the diffusion equation on the free energy surface.