Chemical Reaction Dynamics of the Liquid/Vapor Interface Studied by Mass Spectrometry

This thesis presents investigations of chemical reactions occurring at the liquid/vapor interface studied using novel sampling methodologies coupled with detection by mass spectrometry. Chapters 2 and 3 utilize the recently developed technique of field-induced droplet ionization mass spectrometry (FIDI-MS), in which the application of a strong electric field to a pendant microliter droplet results in the ejection of highly charged progeny droplets from the liquid surface. In Chapter 2, this method is employed to study the base-catalyzed dissociation of a surfactant molecule at the liquid/vapor interface upon uptake of ammonia from the gas phase. This process is observed to occur without significant modulation of the bulk solution pH, suggesting a transient increase in surface pH following the uptake of gaseous ammonia. Chapter 3 presents real-time studies of the oxidation of the model tropospheric organic compound glycolaldehyde by photodissociation of iron (III) oxalate complexes. The oxidation products of glycolaldehyde formed in this process are identified, and experiments in a deoxygenated environment identify the role of oxygen in the oxidation pathway and in the regeneration of iron (III) following photo-initiated reduction. Chapter 4 explores alternative methods for the study of heterogeneous reaction processes by mass spectrometric sampling from liquid surfaces. Bursting bubble ionization (BBI) and interfacial sampling with an acoustic transducer (ISAT) generate nanoliter droplets from a liquid surface that can be sampled via the atmospheric pressure interface of a mass spectrometer. Experiments on the oxidation of oleic acid by ozone using ISAT are also presented. Chapters 5 and 6 detail mechanistic studies and applications of free-radical-initiated peptide sequencing (FRIPS), a technique employing gas-phase free radical chemistry to the sequencing of peptides and proteins by mass spectrometry. Chapter 5 presents experimental and theoretical studies on the anomalous mechanism of dissociation observed in the presence of serine and threonine residues in peptides. Chapter 6 demonstrates the combination of FRIPS with ion mobility-mass spectrometry (IM-MS) for the separation of isomeric peptides.

The Convective Instabilities in a Liquid-Vapor System with a Non-equilibrium Evaporation Interface

Rayleigh-Marangoni-B,nard instability in a system consisting of a horizontal liquid layer and its own vapor has been investigated. The two layers are separated by a deformable evaporation interface. A linear stability analysis is carried out to study the convective instability during evaporation. In previous works, the interface is assumed to be under equilibrium state. In contrast with previous works, we give up the equilibrium assumption and use Hertz-Knudsen's relation to describe the phase change under non-equilibrium state. The influence of Marangoni effect, gravitational effect, degree of non-equilibrium and the dynamics of the vapor on the instability are discussed.

Premelting-Induced Smoothening of the Ice-Vapor Interface

We perform computer simulations of the quasiliquid layer of ice formed at the ice-vapor interface close to the ice Ih-liquid-vapor triple point of water. Our study shows that the two distinct surfaces bounding the film behave at small wavelengths as atomically rough and independent ice-water and water-vapor interfaces. For long wavelengths, however, the two surfaces couple, large scale parallel fluctuations are inhibited, and the ice-vapor interface becomes smooth. Our results could help explain the complex morphology of ice crystallites.

Premelting-Induced Smoothening of the Ice-Vapor Interface

We perform computer simulations of the quasi-liquid layer of ice formed at the ice/vapor interface close to the ice Ih/liquid/vapor triple point of water. Our study shows that the two distinct surfaces bounding the film behave at small wave-lengths as atomically rough and independent ice/water and water/vapor interfaces. For long wave-lengths, however, the two surfaces couple, large scale parallel fluctuations are inhibited and the ice/vapor interface becomes smooth. Our results could help explaining the complex morphology of ice crystallites.

Manipulating double-decker molecules at the liquid-solid interface

We have used a scanning tunneling microscope to manipulate heteroleptic phthalocyaninato, naphthalocyaninato, porphyrinato double-decker molecules at the liquid/solid interface between 1-phenyloctane solvent and graphite. We employed nano-grafting of phthalocyanines with eight octyl chains to place these molecules into a matrix of heteroleptic double-decker molecules; the overlayer structure is epitaxial on graphite. We have also used nano-grafting to place double-decker molecules in matrices of single-layer phthalocyanines with octyl chains. Rectangular scans with a scanning tunneling microscope at low bias voltage resulted in the removal of the adsorbed doubledecker molecular layer and substituted the double-decker molecules with bilayer-stacked phthalocyanines from phenyloctane solution. Single heteroleptic double-decker molecules with lutetium sandwiched between naphthalocyanine and octaethylporphyrin were decomposed with voltage pulses from the probe tip; the top octaethylporphyrin ligand was removed and the bottom naphthalocyanine ligand remained on the surface. A domain of decomposed molecules was formed within the double-decker molecular domain, and the boundary of the decomposed molecular domain self-cured to become rectangular. We demonstrated a molecular “sliding block puzzle” with cascades of double-decker molecules on the graphite surface.

Dynamic structure factor across the liquid-solid interface: appearance of a delta-function elastic peak

We present a theoretical calculation of the dynamic structure factor, S(k, ω), at the liquid-solid interface for large values of the wavevector k. An analytic expression is derived which shows the evolution of the elastic peak as the solid surface is approached from the liquid side.

Waring and Riedel Functions for the Liquid-Vapor Coexistence Curve

There exists a minimum in the Waring function, psi(T) = -d(ln p)/d(1/T), and in the Riedel function, alpha(T) = d(ln p)/d(In T), in the liquid-vapor coexistence curve for most fluids. By analyzing National Institute of Standards and Technology data for the molar enthalpy of vaporization and the compressibility variation at the liquid-vapor phase change of 105 fluids, we find that the temperatures of these minima are linearly correlated with the critical temperature, T-c. Using reduced coordinates, we also demonstrate that the minima are well-correlated with the acentric factor. These correlations are used for testing four well-known vapor pressure equations in the Pitzer corresponding states scheme.

The curvature of the liquid-vapor reduced pressure curve and its relation with the critical region

For most fluids, there exist a maximum and a minimum in the curvature of the reduced vapor pressure curve, p(r) = p(r)(T-r) (with p(r) = p/p(c) and T-r = T/T-c, p(c) and T-c being the pressure and temperature at the critical point). By analyzing National Institute of Standards and Technology (NIST) data on the liquid-vapor coexistence curve for 105 fluids, we find that the maximum occurs in the reduced temperature range 0.5 <= T-r <= 0.8 while the minimum occurs in the reduced temperature range 0.980 <= T-r <= 0.995. Vapor pressure equations for which d(2)p(r)/dT(r)(2) diverges at the critical point present a minimum in their curvature. Therefore, the point of minimum curvature can be used as a marker for the critical region. By using the well-known Ambrose-Walton (AW) vapor pressure equation we obtain the reduced temperatures of the maximum and minimum curvature in terms of the Pitzer acentric factor. The AW predictions are checked against those obtained from NIST data. (C) 2013 Elsevier Ltd. All rights reserved.

Interface capturing simulations of acoustically driven bubble dynamics in and near tissue

Tissue injury during therapeutic ultrasound or lithotripsy is thought, in cases, to be due to the action of cavitation bubbles. Assessing this and mitigating it is challenging since bubble dynamics in the complex confinement of tissues or in small blood vessels are challenging to predict. Simulations tools require specialized algorithms to simultaneously represent strong acoustic waves and shocks, topologically complex liquid‐vapor phase boundaries, and the complex viscoelastic material dynamics of tissue. We discuss advances in a simulation tool for such situations. A single‐mesh Eulerian solver is used to solve the governing equations. Special sharpening terms maintain the liquid‐vapor interface in face of the finite numerical dissipation included in the scheme to accurately capture shocks. A recent enhancement to this formulation has significantly improved this interface capturing procedure, which is demonstrated for simulation of the Rayleigh collapse of a bubble. The solver also transports elastic stresses and can thus be used to assess the effects of elastic properties on bubble dynamics. A shock‐induced bubble collapse adjacent to a model elastic tissue is used to demonstrate this and draw some conclusions regarding the injury suppressing role that tissue elasticity might play.

Spatial modulation of the composition of a binary liquid near a repulsive wall

When a binary liquid is confined by a strongly repulsive wall, the local density is depleted near the wall and an interface similar to that between the liquid and its vapor is formed. This analogy suggests that the composition of the binary liquid near this interface should exhibit spatial modulation similar to that near a liquid-vapor interface even if the interactions of the wall with the two components of the liquid are the same. The Guggenheim adsorption relation quantifies the concentrations of two components of a binary mixture near a liquid-vapor interface and qualitatively states that the majority (minority) component enriches the interface for negative (positive) mixing energy if the surface tensions of the two components are not very different. From molecular dynamics simulations of binary mixtures with different compositions and interactions we find that the Guggenheim relation is qualitatively satisfied at wall-induced interfaces for systems with negative mixing energy at all state points considered. For systems with positive mixing energy, this relation is found to be qualitatively valid at low densities, while it is violated at state points with high density where correlations in the liquid are strong. This observation is validated by a calculation of the density profiles of the two components of the mixture using density functional theory with the Ramakrishnan-Yussouff free-energy functional. Possible reasons for the violation of the Guggenheim relation are discussed.

Maximum Spreading of Liquid Drops Impacting on Groove-Textured Surfaces: Effect of Surface Texture

Maximum, spreading of liquid drops impacting on solid surfaces textured with unidirectional parallel grooves is studied for drop Weber number in the range 1-100 focusing on the role of texture geometry and wettability. The maximum spread factor of impacting drops measured perpendicular to grooves; beta(m,perpendicular to) is seen to be less than, that:measured parallel to grooves, beta(m,perpendicular to).The difference between beta(m,perpendicular to), and beta(m,parallel to) increases with drop impact velocity. This deviation of beta(m,perpendicular to) from beta(m,parallel to) is analyzed by considering the possible mechanisms, correspond, ing to experimental observations (1) impregnation of drop into the grooves, (2) convex shape of liquid vapor interface near contact line at maximum spreading, and (3) contact line pinning of spreading drop at the pillar edges by incorporating them into an energy conservation-based model. The analysis reveals that contact line pinning offers a physically meaningful justification of the observed: deviation of beta(m,perpendicular to) from beta(m,parallel to) compared to other possible candidates. A unified model, incorporating all the above-mentioned mechanisms, is formulated, which predicts beta(m,perpendicular to) on several groove-textured surfaces made of intrinsically hydrophilic and hydrophobic materials with an average error of 8.3%. The effect of groove-texture geometrical parameters,on maximum drop spreading is explained using this unified model. A special case of the unified model, with contact line pinning, absent, predicts beta(m,parallel to) with an average error of 6.3%.

Marangoni-Benard instability with the exchange of evaporation at liquid-vapour interface

A new two-sided model rather than the one-sided model in previous works is put forward. The linear instability analysis is performed on the Marangoni-Benard convection in the two-layer system with an evaporation interface. We define a new evaporation Biot number which is different from that in the one-sided model, and obtain the curves of critical Marangoni number versus wavenumber. The influence of evaporation velocity and Biot number on the system is discussed and a new phenomenon uninterpreted before is now explained from our numerical results.

Electrochemical detector at the liquid/liquid and liquid/membrane interface

Electrochemistry at the liquid/liquid and liquid/membrane interface has been developed recently. It can be used for a new kind of electrochemical detector and sensor. Some new aspects of the field in this laboratory are presented.

Study of adsorption of methylene blue and new methylene blue in liquid-solid interface by slab optical waveguide spectroscopy

A slab optical waveguide (SOWG) has been used for study of adsorption of both methylene blue (MB) and new methylene blue (NMB) in liquid-solid interface. Adsorption characteristics of MB and NMB on both bare SOWG and silanized SOWG by octadecyltrichlorosilane (ODS) were compared. Effect of pH on adsorption on MB and NMB was investigated. Binding rate constant analysis showed that both MB and NMB on bare SOWG demonstrates larger association constants than those on ODS-SOWG. Interactions of NIB and NMB on bare SOWG and ODS-SOWG were analyzed by molecular mechanics calculation method. The binding energy change was in the following order: ENMB-bare > EMB-bare > ENMB-ODS > EMB-ODS. (c) 2004 Elsevier B.V. All rights reserved.

Freezing of 4He and its liquid-solid interface from density functional theory

We show that, at high densities, fully variational solutions of solidlike types can be obtained from a density functional formalism originally designed for liquid 4He . Motivated by this finding, we propose an extension of the method that accurately describes the solid phase and the freezing transition of liquid 4He at zero temperature. The density profile of the interface between liquid and the (0001) surface of the 4He crystal is also investigated, and its surface energy evaluated. The interfacial tension is found to be in semiquantitative agreement with experiments and with other microscopic calculations. This opens the possibility to use unbiased density functional (DF) methods to study highly nonhomogeneous systems, like 4He interacting with strongly attractive impurities and/or substrates, or the nucleation of the solid phase in the metastable liquid.