Surface vibrational spectroscopic study of surface melting of ice

Sublimation: A Two-Step Phase Transition

It is shown that the existing notions of the triple point phase equilibrium and the sublimation process need revision. To understand these phenomena a different approach is proposed, substantiated with the principle of least time for phase transitions. It is deduced that sublimation is a two-step phase transition through an intermediate (liquid) state.The new approach is consistent with the modern viewpoint of the first-order phase transitions. Some experimental data support the proposed approach.

Radiation-induced premelting of ice at silica interfaces

The existence of surface and interfacial melting of ice below 0 degrees C has been confirmed by many different experimental techniques. Here we present a high-energy x-ray reflectivity study of the interfacial melting of ice as a function of both temperature and x-ray irradiation dose. We found a clear increase of the thickness of the quasiliquid layer with the irradiation dose. By a systematic x-ray study, we have been able to unambiguously disentangle thermal and radiation-induced premelting phenomena. We also confirm the previously announced very high water density (1.25 g/cm(3)) within the emerging quasiliquid layer.

Interfacial structures of acidic and basic aqueous solutions

Phase-sensitive sum-frequency vibrational spectroscopy was used to study water/vapor interfaces of HCl, HI, and NaOH solutions. The measured imaginary part of the surface spectral responses provided direct characterization of OH stretch vibrations and information about net polar orientations of water species contributing to different regions of the spectrum. We found clear evidence that hydronium ions prefer to emerge at interfaces. Their OH stretches contribute to the "ice-like" band in the spectrum. Their charges create a positive surface field that tends to reorient water molecules more loosely bonded to the topmost water layer with oxygen toward the interface, and thus enhances significantly the "liquid-like" band in the spectrum. Iodine ions in solution also like to appear at the interface and alter the positive surface field by forming a narrow double-charge layer with hydronium ions. In NaOH solution, the observed weak change of the "liquid-like" band and disappearance of the "ice-like" band in the spectrum indicates that OH(-) ions must also have excess at the interface. How they are incorporated in the interfacial water structure is, however, not clear.

Organic contaminant amplification during snowmelt

The release of organic contaminants from melting snow poses risks to aquatic and terrestrial organisms and to humans who rely on drinking water and food production from regions that are seasonally snow-covered. Measured and model-predicted spring peak concentrations in waters receiving snowmelt motivate a thorough investigation of organic contaminant behaviour during melting. On the basis of the current understanding of snow metamorphosis, snowmelt hydrology and chemical partitioning in snow, this critical review aims to provide a qualitative picture of the processes involved in the release of organic contaminants from a melting snowpack. The elution sequence of organic substances during snowmelt is strongly dependent on their environmental partitioning properties and the physical properties of the snowpack. Water-soluble organic contaminants can be discharged in greatly elevated concentrations at an early stage of melting, while the bulk of the hydrophobic chemicals attached to particles is often released at the end of the melt period. Melting of a highly metamorphosed and deep snowpack promotes such shock load releases, whereas a shallow snow cover over a relatively warm ground experiencing irregular melting over the winter season is unlikely to generate notable peak releases of organic substances. Meltwater runoff over frozen ground directly transfers contaminant shock loads into receiving water bodies, while permeable soils buffer and dilute the contaminants. A more quantitative understanding of the behaviour of organic contaminants in varying snowmelt scenarios will depend on controlled laboratory studies combined with field investigations. Reliable numerical process descriptions will need to be developed to integrate water quality and contaminant fate models.

Characterization of vibrational resonances of water-vapor interfaces by phase-sensitive sum-frequency spectroscopy

Phase-sensitive sum-frequency spectroscopy provides correct characterization of vibrational resonances of water-vapor interfaces and allows better identification of interfacial water species contributing to different parts of the spectra. Iodine ions emerging at an interface create a surface field that tends to reorient the more loosely bonded water molecules below the topmost layer.

The single-crystal, basal face of ice I(h) investigated with sum frequency generation

Sum frequency generation spectroscopy has been used to investigate the hydrogen-bonded region of single-crystal, hexagonal ice in the temperature range of 113-178 K. The temperature and polarization dependences of the signal are used in conjunction with a recent theoretical model to suggest an interpretation of the bluest and reddest of the hydrogen-bonded peaks. The reddest feature is associated with strong hydrogen bonding; the dynamic polarizability of this feature is primarily parallel to the surface. It is assigned to a cooperative motion among the companion to the free-OH and four-coordinate oscillators hydrogen bonded to dangling lone-pair molecules on the surface. The bluest hydrogen-bonded feature is similarly assigned to a cooperative motion of the OH stretch of dangling lone-pair molecules and of four-coordinate molecules in the lower half bilayer that are hydrogen bonded to free-OH molecules. Reconstruction induced strain is present at as low as 113 K. These results provide a richer picture of the ice surface than has heretofore been possible.

Excess electron relaxation dynamics at water/air interfaces

We have performed mixed quantum-classical molecular dynamics simulations of the relaxation of a ground state excess electron at interfaces of different phases of water with air. The investigated systems included ambient water/air, supercooled water/air, Ih ice/air, and amorphous solid water/air interfaces. The present work explores the possible connections of the examined interfacial systems to finite size cluster anions and the three-dimensional infinite, fully hydrated electron. Localization site analyses indicate that in the absence of nuclear relaxation the electron localizes in a shallow potential trap on the interface in all examined systems in a diffuse, surface-bound (SB) state. With relaxation, the weakly bound electron undergoes an ultrafast localization and stabilization on the surface with the concomitant collapse of its radius. In the case of the ambient liquid interface the electron slowly (on the 10 ps time scale) diffuses into the bulk to form an interior-bound state. In each other case, the excess electron persists on the interface in SB states. The relaxation dynamics occur through distinct SB structures which are easily distinguishable by their energetics, geometries, and interactions with the surrounding water bath. The systems exhibiting the most stable SB excess electron states (supercooled water/air and Ih ice/air interfaces) are identified by their characteristic hydrogen-bonding motifs which are found to contain double acceptor-type water molecules in the close vicinity of the electron. These surface states correlate reasonably with those extrapolated to infinite size from simulated water cluster anions.

On the nature of ions at the liquid water surface

A qualitatively new understanding of the nature of ions at the liquid water surface is emerging. Traditionally, the characterization of liquid surfaces has been limited to macroscopic experimental techniques such as surface tension and electrostatic potential measurements, wherein the microscopic picture then has to be inferred by applying theoretical models. Because the surface tension of electrolyte solutions generally increases with ion concentration, all inorganic ions have been thought to be repelled from the air-water interface, leaving the outermost surface layer essentially devoid of ions. This oversimplified picture has recently been challenged: first by chemical kinetics measurements, then by theoretical molecular dynamics simulations using polarizable models, and most recently by new surface sensitive experimental observations. Here we present an overview of the nature of the interfacial structure of electrolyte solutions and give a detailed description of the new picture that is emerging.

Photolysis of Polycyclic Aromatic Hydrocarbons on Water and Ice Surfaces

Laser-induced fluorescence detection was used to measure photolysis rates of anthracene and naphthalene at the air-ice interface, and the kinetics were compared to those observed in water solution and at the air-water interface. Direct photolysis proceeds much more quickly at the air-ice interface than at the air-water interface, whereas indirect photolysis due to the presence of nitrate or hydrogen peroxide appears to be suppressed at the ice surface with respect to the liquid water surface. Both naphthalene and anthracene self-associate readily on the ice surface, but not on the water surface. The increase in photolysis rates observed on ice surfaces is not due to this self-association, however. The wavelength dependence of the photolysis indicates that it is due to absorption by the PAH. No dependence of the rate on temperature is seen, either at the liquid water surface or at the ice surface. Molecular oxygen appears to play a complex role in the photolytic loss mechanism, increasing or decreasing the photolysis rate depending on its concentration.

Dynamics of TIP5P and TIP4P/ice potentials

The dynamics of a thin film of ice Ih deposited on MgO (001) is studied through molecular dynamics simulations performed with two new potential models of ice. This system is chosen because it is possible to compare the results of the simulations to incoherent neutron quasielastic scattering experiments performed few years ago and to previous molecular dynamics simulations using the TIP4P potential model. The present simulations are performed to determine the evolution of the translational and orientational order parameters of the ice film upon temperature increase in the 250-280 K range. They are also used to calculate the translational and orientational diffusion coefficients of the water molecules in the supported film as a function of the temperature. When using the TIP5P potential, the present results show a better agreement with experimental data than those calculated with the TIP4P potential, especially regarding the temperature above which significant changes are obtained in the dynamics of the water film. Similar conclusions are obtained when using the TIP4P/ice potential, although this latter potential clearly underestimates the translational diffusion coefficients.

Surface properties of water ice at 150–191K studied by elastic helium scattering

A highly surface sensitive technique based on elastic scattering of low-energy helium atoms has been used to probe the conditions in the topmost molecular layer on ice in the temperature range of 150-191 K. The elastically scattered intensity decreased slowly as the temperature was increased to about 180 K, followed by a rapid decrease at higher temperatures. An effective surface Debye temperature of 185+/-10 K was calculated from the data below 180 K. The changes in the ice surface above 180 K are interpreted as the onset of an anomalous enhancement of the mean square vibrational amplitude for the surface molecules and/or the onset of a limited amount of disorder in the ice surface. The interpretation is consistent with earlier experimental studies and molecular dynamics simulations. The observed changes above 180 K can be considered as the first sign of increased mobility of water molecules in the ice surface, which ultimately leads to the formation of a quasiliquid layer at higher temperatures. A small shift and broadening of the specular peak was also observed in the range of 150-180 K and the effect is explained by the inherent corrugation of the crystalline ice surface. The peak shift became more pronounced with increasing temperature, which indicates that surface corrugation increases as the temperature approaches 180 K. The results have implications for the properties and surface chemistry of atmospheric ice particles, and may contribute to the understanding of solvent effects on the internal molecular motion of hydrated proteins and other organic structures such as DNA.

Surface freezing in binary alkane-alcohol mixtures

Surface freezing was detected and studied in mixtures of alcohol and alkane molecules, using surface tensiometry and surface-specific x-ray scattering methods. Considering that surface freezing in pure alkanes forms an ordered monolayer and in alcohols it forms an ordered bilayer, the length mismatch repulsion was minimized by varying the carbon number of the alkane component around 2n, where n is the carbon number of the alcohol molecule. A solutionlike behavior was found for all mixtures, where the ideal liquid mixture phase-separates upon freezing both in the bulk and the surface. The solid exhibits a herringbone crystalline phase below an alkane mole fraction phi(t) approximately 0.8 and a rotator phase above it. The surface frozen film below phi(t) is an alkane monolayer exhibiting a next-nearest neighbor molecular tilt of a composition-dependent magnitude. Above phi(t), no diffraction peaks were observed. This could be explained by the intrinsically shorter-range order of the rotator phase and a possible proliferation of defects.

CORRELATED ELECTRONIC STRUCTURE NONLINEAR RESPONSE METHODS FOR STRUCTURED ENVIRONMENTS

▪ Abstract  This contribution concerns a brief outline of structural environment models where correlated electronic structure response methods are utilized for the determination of nonlinear optical properties of molecules. The presentation provides theory and applications of a heterogeneous dielectric media model and a quantum mechanical–classical mechanical model at the level of correlated electronic structure response methods. The correlated electronic structure response methods include the multiconfigurational self-consistent field method.

Polarization and experimental configuration analyses of sum frequency generation vibrational spectra, structure, and orientational motion of the air/water interface

Here we report a detailed study on spectroscopy, structure, and orientational distribution, as well as orientational motion, of water molecules at the air/water interface, investigated with sum frequency generation vibrational spectroscopy (SFG-VS). Quantitative polarization and experimental configuration analyses of the SFG data in different polarizations with four sets of experimental configurations can shed new light on our present understanding of the air/water interface. Firstly, we concluded that the orientational motion of the interfacial water molecules can only be in a limited angular range, instead of rapidly varying over a broad angular range in the vibrational relaxation time as suggested previously. Secondly, because different vibrational modes of different molecular species at the interface has different symmetry properties, polarization and symmetry analyses of the SFG-VS spectral features can help the assignment of the SFG-VS spectra peaks to different interfacial species. These analyses concluded that the narrow 3693 cm(-1) and broad 3550 cm(-1) peaks belong to C(infinityv) symmetry, while the broad 3250 and 3450 cm(-1) peaks belong to the symmetric stretching modes with C2v symmetry. Thus, the 3693 cm(-1) peak is assigned to the free OH, the 3550 cm(-1) peak is assigned to the singly hydrogen-bonded OH stretching mode, and the 3250 and 3450 cm(-1) peaks are assigned to interfacial water molecules as two hydrogen donors for hydrogen bonding (with C2v symmetry), respectively. Thirdly, analysis of the SFG-VS spectra concluded that the singly hydrogen-bonded water molecules at the air/water interface have their dipole vector directed almost parallel to the interface and is with a very narrow orientational distribution. The doubly hydrogen-bonded donor water molecules have their dipole vector pointing away from the liquid phase.

Nonlinear Optical Studies of the Agricultural Antibiotic Morantel Interacting with Silica/Water Interfaces

It is now known that the untreated discharge of pharmaceuticals into the environment can impact human health and development and lead to increased drug resistance in biota. Here, we present the first direct interface-specific studies that address the mobility of the widely used agricultural antibiotic morantel, which is commonly present in farm runoff. Surface-bound morantel was spectroscopically identified using second harmonic generation (SHG) via a two-photon resonance of its n-pi* transition and in the C-H stretching region by vibrational sum frequency generation (VSFG). Resonantly enhanced SHG adsorption isotherm measurements carried out at the silica/water interface between 6 x 10(-7) and 5 x 10(-5) M morantel concentration result in a free energy of adsorption of 42(2) kJ/mol at pH 7. Finally, real-time tracking of morantel interaction with the silica/water interface shows that the binding events are fully reversible, consistent with its high mobility in silica-rich soil environments. This work thus indicates that pharmaceuticals discharged into the environment can enter the groundwater supply of municipal water systems, at which point their removal is challenging. In addition, the high mobility of morantel in silica-rich soil environments could lead to developing increased interaction of this antibiotic with target organisms, which could respond by increased drug resistance.

Thermodynamic model of quasiliquid formation on H2O ice: Comparison with experiment

We have developed a new thermodynamic theory of the quasiliquid layer, which has been shown to be effective in modeling the phenomenon in a number of molecular systems. Here we extend our analysis to H(2)O ice, which has obvious implications for environmental and atmospheric chemistry. In the model, the liquid layer exists in contact with an ice defined as a two-dimensional lattice of sites. The system free energy is defined by the bulk free energies of ice I(h) and liquid water and is minimized in the grand canonical ensemble. An additional configurational entropy term arises from the occupation of the lattice sites. Furthermore, the theory predicts that the layer thickness as a function of temperature depends only on the liquid activity. Two additional models are derived, where slightly different approximations are used to define the free energy. With these two models, we illustrate the connection between the quasiliquid phenomenon and multilayer adsorption and the possibility of a two-dimensional phase transition connecting a dilute low coverage phase of adsorbed H(2)O and the quasiliquid phase. The model predictions are in agreement with a subset of the total suite of experimental measurements of the liquid thickness on H(2)O ice as a function of temperature. The theory indicates that the quasiliquid layer is actually equivalent to normal liquid water, and we discuss the impact of such an identification. In particular, observations of the liquid layer to temperatures as low as 200 K indicate the possibility that the quasiliquid is, in fact, an example of deeply supercooled normal water. Finally, we briefly discuss the obvious extension of the pure liquid theory to a thermodynamic theory of interfacial solutions on ice in the environment.

Kinetic Studies of Chromium (VI) Binding to Carboxylic Acid- and Methyl Ester-Functionalized Silica/Water Interfaces Important in Geochemistry

Real-time kinetic measurements of hexavalent chromium binding to fused silica surfaces functionalized with carboxylic acid and methyl ester terminal groups are performed in situ using resonantly enhanced surface second harmonic generation (SHG) at pH 7 and 300 K. These functional groups were chosen because of their high abundance in humic acids and related biopolymers. Kinetic measurements are conducted in the submonolayer regime using chromate solution concentrations ranging from 1 x 10(-6) to 2 x 10(-5) M. The adsorption rates were analyzed using the standard Langmuir model and the Frumkin-Fowler-Guggenheim model. The desorption kinetics are consistent with a first-order process. These results indicate that hexavalent chromium mobility in carboxylic acid- and ester-rich soil environments increases with decreasing chromate concentrations. Based on the measured half-lives of the adsorbed Cr(VI) species, remobilization of bound hexavalent chromium due to natural or anthropogenic events that lower the chromate concentration in the aqueous phase can occur within minutes.

Interfacial melting of ice in contact with SiO(2)

The physical behavior of condensed matter can be drastically altered in the presence of interfaces. Using a high-energy x-ray transmission-reflection scheme, we have studied ice-SiO2 model interfaces. We observed the formation of a quasiliquid layer below the bulk melting temperature and determined its thickness and density as a function of temperature. The quasiliquid layer has stronger correlations than water and a large density close to rho(HDA)=1.17 g/cm(3) of high-density amorphous ice suggesting a structural relationship with the postulated high-density liquid phase of water.

Sum frequency vibrational spectroscopy of leucine molecules adsorbed at air–water interface

Sum frequency vibrational spectroscopy was used to study adsorption of leucine molecules at air-water interface from solutions with different concentrations and pH values. The surface density and the orientation of the isopropyl head group of the adsorbed leucine molecules could be deduced from the measurements. It was found that the orientation depends on the surface density, but only weakly on bulk pH value at the saturated surface density. The vibrational spectra of the interfacial water molecules appeared to be strongly affected by the charge state of the adsorbed leucine molecules. Enhancement and inversion of polar orientation of interfacial water molecules by surface charges or field controllable by the bulk pH value were observed.

Surface freezing in binary mixtures of chain molecules. I. Alkane mixtures

X-ray surface scattering and surface tension measurements are used to study surface freezing in molten mixtures of alkanes. These binary mixtures consist of protonated and deuterated alkanes, as well as of alkanes of different lengths. As for pure alkanes, a crystalline monolayer is formed at the surface a few degrees above the bulk freezing temperature. The structure of the monolayer has been determined on an angstrom scale. A simple theoretical approach is used to account for the thermodynamical observations at the surface and in the bulk. The model is based on a competition between entropic mixing and a repulsive interaction due to chain-length mismatch. The surface and bulk liquid phases are treated as ideal mixtures, while the solid phases are treated as regular mixtures. The theory is found to account well for all the mixtures studied, both hydrogenated-hydrogenated and hydrogenated-deuterated. The repulsive interaction and its dependence on the chain lengths of the components are determined from fits to the measured data.

Air-snow interactions and atmospheric chemistry

The presence of snow greatly perturbs the composition of near-surface polar air, and the higher concentrations of hydroxyl radicals (OH) observed result in a greater oxidative capacity of the lower atmosphere. Emissions of nitrogen oxides, nitrous acid, light aldehydes, acetone, and molecular halogens have also been detected. Photolysis of nitrate ions contained in the snow appears to play an important role in creating these perturbations. OH formed in the snowpack can oxidize organic matter and halide ions in the snow, producing carbonyl compounds and halogens that are released to the atmosphere or incorporated into snow crystals. These reactions modify the composition of the snow, of the interstitial air, and of the overlying atmosphere. Reconstructing the composition of past atmospheres from ice-core analyses may therefore require complex corrections and modeling for reactive species.