Water-hydrocarbon interfaces: effect of hydrocarbon branching on interfacial structure

Molecular dynamics simulation are performed for the water/hydrocarbon system to study the effect of hydrocarbon branching on interfacial properties. The following two series of hydrocarbons are considered: (1) n-pentane, 2-methyl pentane, and 2,2,4-trimethyl pentane (constant chain length) and (2) n-octane, 2-methyl heptane, and 2,2,4-trimethyl pentane (constant molecular mass). With a simple algorithm for identification of surface sites and mapping nonsurface sites to these surface sites, intrinsic profiles were constructed with respect to the surface layer. Intrinsic density profiles for water and hydrocarbons with respect to the hydrocarbon and water surface, respectively, resemble density profiles of liquids in the presence of a wall. Order parameters were used to study orientation of molecules with respect to the surface normal and the hydrogen bond network was characterized in terms of the number of hydrogen bonds per water molecule and percentage of hydrogen bonded molecules in the first coordination shell. The corresponding intrinsic profiles were obtained. The O-H bond for surface water was found to have two preferential orientations, pointing toward the hydrocarbon phase and parallel to the interface. Hydrocarbon molecules in series 1 orient along the interface with the more branched molecule better aligned. For molecules in series 2, the larger molecular length reduces the alignment of molecules along the interface.     

Structure of the nonionic surfactant triethoxy monooctylether C8E3 adsorbed at the free water surface, as seen from surface tension measurements and Monte Carlo simulations

The adsorption layer of the nonionic surfactant triethoxy monooctylether C8E3 has been investigated at the free water surface by means of both experimental and computer simulation methods. The surface tension of the aqueous solution of C8E3 has been measured by pendant drop shape analysis in the entire concentration range in which C8E3 is soluble in water. The data obtained from these measurements are used to derive the adsorption isotherm. The critical micellar concentration and the surface excess concentration of the saturated adsorption layer are found to be 7.48 mM and 4.03 micromol/m2, respectively, the latter value corresponding to the average area per molecule of 41 A2. In order to analyze the molecular level structure of the unsaturated adsorption layer, Monte Carlo simulations have been performed at four different surface concentration values, i.e., 0.68, 1.36, 2.04, and 2.72 micromol/m2, respectively. It has been found that the water surface is already almost fully covered at the lowest surface density value investigated, and the adsorbed molecules show a strong preference for lying parallel with the interface in elongated conformations. No sign of the penetration of the hydrophilic triethoxy headgroups into the aqueous phase to any extent has been observed. With increasing surface densities the preferential orientation of the apolar octyl tails gradually turns from lying parallel with the interface to pointing toward the vapor phase by their CH3 end, whereas the conformation of the adsorbed molecules becomes gradually less elongated. Both of these changes lead to the increase of the number of C8E3 molecules being in a direct contact (i.e., forming hydrogen bonds) with water. However, the increasing number of the C8E3 molecules hydrogen bonded to water is found to be accompanied by the weakening of this binding, i.e., the decrease of both the number of hydrogen bonds a bound C8E3 molecule forms with water and the magnitude of the average binding energy of the adsorbed C8E3 molecules.     

Mechanism of OH radical hydration: a comparative computational study of liquid and supercritical solvent

Flexible models of the radical and water molecules including short-range interaction of hydrogen atoms have been employed in molecular dynamic simulation to understand mechanism of (●)OH hydration in aqueous systems of technological importance. A key role of H-bond connectivity patterns of water molecules has been identified. The behavior of (●)OH(aq) strongly depends on water density and correlates with topological changes in the hydrogen-bonded structure of water driven by thermodynamic conditions. Liquid and supercritical water above the critical density exhibit the radical localization in cavities existing in the solvent structure. A change of mechanism has been found at supercritical conditions below the critical density. Instead of cavity localization, we have identified accumulation of water molecules around (●)OH associated with the formation of a strong H-donor bond and diminution of non-homogeneity in the solvent structure. For all the systems investigated, the computed hydration number and the internal energy of hydration Δ(h)U showed approximately linear decrease with decreasing density of the solvent but a degree of radical-water hydrogen bonding exhibited non-monotonic dependence on density. The increase in the number of radical-water H-acceptor bonds is associated with diminution of extended nets of four-bonded water molecules in compressed solution at ~473 K. Up to 473 K, the isobaric heat of hydration in compressed liquid water remains constant and equal to -40 ± 1 kJ mol(-1).     

Molecular dynamics of phenol at the liquid-vapor interface of water

Molecular dynamics results are presented for phenol at the water liquid-vapor interface at 300 K. The calculated excess free energy of phenol at the interface is -2.8 +/- 0.4 kcal/mol, in good agreement with the recent experimental results of Eisenthal and co-workers. The most probable orientation of the phenol molecule at the surface is such that the aromatic ring is perpendicular to the interface and the OH group is fully immersed in water. The hydroxyl substituent has a preferred orientation which is similar to the orientation of OH bonds of water at the pure water liquid-vapor interface. The transition between interfacial and bulk-like behavior of phenol is abrupt and occurs when the center of mass of the solute is located about 6 angstroms from the Gibbs surface of water. In this region the para carbon atom of the hydrophobic benzene ring can reach the interface and become partially dehydrated. This result suggests that the width of the interfacial region in which the behavior of a simple amphiphilic solute in water is influenced by the presence of the surface depends primarily on the size of its hydrophobic part. The role of the OH substituent was investigated by comparing phenol at the interface with two model systems: benzene with and without partial charges on carbon and hydrogen atoms. It is shown that in the absence of the hydrophilic substituent the solute is located further away from the liquid phase and is more likely to be oriented parallel to the interface. However, when the center of mass of the solute is moved into the interfacial region where the density of water approaches that of the bulk solvent, all three molecules become oriented perpendicularly to the surface. In this orientation the work of cavity formation needed to accommodate the hydrophobic ring in aqueous solvent is minimized.     

金纳米颗粒周围水的结构和动力学性质的分子动力学模拟

采用经典的分子动力学模拟方法系统地研究了在常温条件下金纳米颗粒周围水的结构与动力学性质. 结果表明, 水分子在纳米颗粒附近形成了明显的多层结构. 同时随着径向距离的减小, 水分子的空间取向也从无序排列趋向于有序排列. 通过分析界面处不同水层中的均方位移及停留时间分布, 发现紧贴颗粒表面的第一和第二水层中的水分子表现出很低的扩散系数, 而第三和第四水层中的水分子则能够轻易地离开界面区域而进入主体相区域. 此外, 在界面处的每个水分子的氢键平均数要高于在主体相的平均值.     

Aqueous Heterogeneity at the Air/Water Interface Revealed by 2D‐HD‐SFG Spectroscopy

Water molecules interact strongly with each other through hydrogen bonds. This efficient intermolecular coupling causes strong delocalization of molecular vibrations in bulk water. We study intermolecular coupling at the air/water interface and find intermolecular coupling 1) to be significantly reduced and 2) to vary strongly for different water molecules at the interface—whereas in bulk water the coupling is homogeneous. For strongly hydrogen‐bonded OH groups, coupling is roughly half of that of bulk water, due to the lower density in the near‐surface region. For weakly hydrogen‐bonded OH groups that absorb around 3500 cm^?1, which are assigned to the outermost, yet hydrogen‐bonded OH groups pointing towards the liquid, coupling is further reduced by an additional factor of 2. Remarkably, despite the reduced structural constraints imposed by the interfacial hydrogen‐bond environment, the structural relaxation is slow and the intermolecular coupling of these water molecules is weak.     

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.     

Quantum Chemical Study of Water Adsorption on the Surfaces of SrTiO3 Nanotubes

We have studied the adsorption of water molecules on the inner and outer surfaces of nanotubes generated by rolling (001) layers of SrTiO₃ cubic crystals. The stability and the atomic and electronic structures of the adsorbed layers are determined by using hybrid density functional theory. The absorption energy and the preferred adsorbate structure are essentially governed by the nature of the surface of the nanotube. Dissociative adsorption prevails on the outer nanotube surfaces. The stability of the adsorbed layers on the inner surfaces is related to the possibility of the formation of hydrogen bonds between water molecules and surface oxygen atoms, and depends on the surface curvature. The presence of water molecules on the inner surface of the nanotubes leads to an increase of the electronic band gap. Externally TiO₂‐terminated nanotubes could be used for the photocatalytic decomposition of water by ultraviolet radiation.     

Structural characterization of interfacial n-octanol and 3-octanol using molecular dynamic simulations

Structurally isomeric octanol interfacial systems, water/vapor, 3-octanol/vapor, n-octanol/vapor, 3-octanol/water, and n-octanol/water are investigated at 298 K using molecular dynamics simulation techniques. The present study is intended to investigate strongly associated liquid/liquid interfaces and probe the atomistic structure of these interfaces. The octanol and water molecules were initially placed randomly into a box and were equilibrated using constant pressure techniques to minimize bias within the initial conditions as well as to fully sample the structural conformations of the interface. An interface formed via phase separation during equilibration and resulted in a slab geometry with a molecularly sharp interface. However, some water molecules remained within the octanol phase with a mole fraction of 0.12 after equilibration. The resulting "wet" octanol interfaces were analyzed using density profiles and orientational order parameters. Our results support the hypothesis of an ordered interface only 1 or 2 molecular layers deep before bulk properties are reached for both the 3-octanol and water systems. However, in contrast to most other interfacial systems studied by molecular dynamics simulations, the n-octanol interface extends for several molecular layers. The octanol hydroxyl groups form a hydrogen-bonding network with water which orders the surface molecules toward a preferred direction and produces a hydrophilic/hydrophobic layering. The ordered n-octanol produces an oscillating low-high density of oxygen atoms out of phase with a high-low density of carbon atoms, consistent with an oscillating dielectric. In contrast, the isomeric 3-octanol has only a single carbon-rich layer directly proximal to the interface, which is a result of the different molecular topology. Both 3-octanol and n-octanol roughen the water interface with respect to the water/vapor interface. The "wet" octanol phases, in the octanol/water systems reach bulk properties in a shorter distance than the "dry" octanol/vapor interfaces.     

Adsorption of octyl cyanide at the free water surface as studied by Monte Carlo simulation

Monte Carlo simulations of the adsorption layer of octyl cyanide have been performed on the canonical (N, V, T) ensemble at 300 K. The systems simulated cover the range of octyl cyanide surface densities from 0.27 to 7.83 mumol/m2. The surface density value at which the saturation of the adsorption layer occurs is estimated to be 1.7 mumol/m2. At low surface densities, the main driving force of the adsorption is found to be the formation of hydrogen bonds between the water and octyl cyanide molecules, whereas at higher surface concentrations, the dipole-dipole attraction between the neighboring adsorbed octyl cyanide molecules becomes more important. At low surface concentrations, the water-octyl cyanide hydrogen bonds prefer tilted alignments relative to the interface; however, in the case of the saturated adsorption layer, the number of such hydrogen bonds is maximized, leading to the preference of these bonds for the orientation perpendicular to the interface. Contrary to nonionic surfactants of multiple hydrogen bonding abilities (e.g., 1-octanol, C8E3), the increasing surface concentration of octyl cyanide was not found to lead to considerable competition of the molecules for positions of optimal arrangement. As a consequence, the energy and geometry of the water-octyl cyanide hydrogen bonds are found to be insensitive to the octyl cyanide surface concentration.     

Molecular dynamics study of the hydration of the hydroxyl radical at body temperature

Classical molecular dynamics (MD) simulation of ˙OH in liquid water at 37 °C has been performed using flexible models of the solute and solvent molecules. We derived the Morse function describing the bond stretching of the radical and the potential for ˙OH-H(2)O interactions, including short-range interactions of hydrogen atoms. Scans of the potential energy surface of the ˙OH-H(2)O complex have been performed using the DFT method with the B3LYP functional and the 6-311G(d,p) basis set. The DFT-derived partial charges, ±0.375e, and the equilibrium bond-length, 0.975 ?, of ˙OH resulted in the dipole moment of 1.76 D. The radical-water radial distribution functions revealed that ˙OH is not built into the solvent structure but it rather occupies distortions or cavities in the hydrogen-bonded network. The solvent structure at 37 °C has been found to be the same as that of pure water. The hydration cage of the radical comprises 13-14 water molecules. The estimated hydration enthalpy -42 ± 5 kJ mol(-1) is comparable with the experimental value -39 ± 6 kJ mol(-1) for 25 °C. Inspection of hydrogen bonds showed the importance of short-range interaction of hydrogen atoms and indicated that neglect of the angular condition greatly overestimates the number of the H-acceptor radical-water bonds. The mean number ?n = 0.85 of radical-water H-bonds has been calculated using geometric definition of H-bond and ?n = 0.62 has been obtained when the energetic condition, E(da)≤-8 kJ mol(-1), was additionally considered. The continuous lifetimes of 0.033 ps and 0.024 ps have been estimated for the radical H-donor and the H-acceptor bonds, respectively. Within statistical uncertainty the radical self-diffusion coefficient, (2.9 ± 0.6) × 10(-9) m(2) s(-1), is the same as (3.1 ± 0.5) × 10(-9) m(2) s(-1) calculated for water in solution and in pure solvent. To the best of our knowledge, this is the first study of the ˙OH(aq) properties at a biologically relevant body temperature.     

Molecular Structure and Dynamics of Water at the Water–Air Interface Studied with Surface‐Specific Vibrational Spectroscopy

Water interfaces provide the platform for many important biological, chemical, and physical processes. The water–air interface is the most common and simple aqueous interface and serves as a model system for water at a hydrophobic surface. Unveiling the microscopic (<1 nm) structure and dynamics of interfacial water at the water–vapor interface is essential for understanding the processes occurring on the water surface. At the water interface the network of very strong intermolecular interactions, hydrogen‐bonds, is interrupted and the density of water is reduced. A central question regarding water at interfaces is the extent to which the structure and dynamics of water molecules are influenced by the interruption of the hydrogen‐bonded network and thus differ from those of bulk water. Herein, we discuss recent advances in the study of interfacial water at the water–air interface using laser‐based surface‐specific vibrational spectroscopy.     

Molecular dynamics study of free energy of transfer of alcohol and amine from water phase to the micelle by thermodynamic integration method

Free energy of transfer of methylamine, octylamine, methanol, and octanol from water phase to sodium dodecyl sulfate (SDS) micelle has been calculated using thermodynamic integration method combined with molecular dynamics calculations. Together with the results for alkanes obtained in our previous study [K. Fujimoto, N. Yoshii, and S. Okazaki, J. Chem. Phys. 133, 074511 (2010)], the effect of polar group on the partition of hydrophilic solutes between water phase and the micelle has been investigated in detail at a molecular level. The calculations showed that the molecules with octyl group are more stable in the SDS micelle than in the water phase due to their hydrophobicity of long alkyl chain. In contrast, methanol and methylamine are stable in the water phase as well as in the micelle because of their high hydrophilicity. The spatial distribution of methylamine, octylamine, methanol, and octanol has also been evaluated as a function of the distance, R, from the center of mass of SDS micelle to the solutes. The distribution shows that the methylamine molecule is adsorbed on the SDS micelle surface, while the methanol molecule is delocalized among the whole system, i.e., in the water phase, on the surface of the micelle, and in the hydrophobic core of the micelle. The octylamine and octanol molecules are solubilized in the SDS micelle with palisade layer structure and are not found in the water phase.     

Spectroscopic characterization of microscopic hydrogen-bonding disparities in supercritical water

The local hydrogen-bonding environment in supercritical water (380 degrees C, 300 bars, density 0.54 gcm3) was studied by x-ray Raman scattering at the oxygen K edge. The spectra are compared to those of the gas phase, liquid surface, bulk liquid, and bulk ice, as well as to calculated spectra. The experimental model systems are used to assign spectral features and to quantify specific local hydrogen-bonding situations in supercritical water. The first coordination shell of the molecules is characterized in more detail with the aid of the calculations. Our analysis suggests that approximately 65% of the molecules in supercritical water are hydrogen bonded in configurations that are distinctly different from those in liquid water and ice. In contrast to liquid water the bonded molecules in supercritical water have four intact hydrogen bonds and in contrast to ice large variations of bond angles and distances are observed. The remaining approximately 35% of the molecules exhibit two free O-H bonds and are thus either not involved in hydrogen bonding at all or have one or two hydrogen bonds on the oxygen side. We determine an average O-O distance of 3.1+/-0.1 A in supercritical water for the H bonded molecules at the conditions studied here. This and the corresponding hydrogen bond lengths are shown to agree with neutron- and x-ray-diffraction data at similar conditions. Our results on the local hydrogen-bonding environment with mainly two disparate hydrogen-bonding configurations are consistent with an extended structural model of supercritical water as a heterogeneous system with small patches of bonded molecules in various tetrahedral configurations and surrounding nonbonded gas-phase-like molecules.     

Water adsorption on ZnO(1010): from single molecules to partially dissociated monolayers

Static and dynamic density functional calculations have been used to study the structure and energetics of water adsorbed on the main cleavage plane of ZnO. In the single molecule limit we find that molecular adsorption is strongly preferred. The water binding energy increases for higher coverages due to an almost isotropic attractive water-water interaction which leads to clustering and formation of monolayer islands in the low water coverage regime. A thermodynamic analysis further shows that the full water monolayer is clearly the most stable phase until water starts to desorb. The water monolayer is even more stabilized by a partial dissociation of the water molecules, yielding as most stable configuration a (2x1) superstructure where every second water molecule is cleaved. The dissociation barrier for this process is very small which allows for an auto-dissociation of the water molecules even at low temperatures as observed experimentally. Finally we find that the energy cost involved to form [1210]-oriented domain boundaries between (2x1) patches with different orientation is almost negligible which explains the abundance of such domain boundaries in STM images.     

A Theoretical Study on Hydrogen‐Bonded Complex of Proflavine Cation and Water: The Site‐dependent Feature of Hydrogen Bond Strengthening and Weakening

The intermolecular hydrogen‐bonds between proflavine cation (PC) and water molecules are investigated by density functional theory (DFT) and time‐dependent density functional theory (TDDFT) methods. The ground‐state geometry optimizations, electronic excitation energies and corresponding oscillation strengths of the low‐lying electronically excited states for the isolated proflavine cation, the hydrogen‐bonded PC–H2O dimer and PC–(H2O)2 trimer are calculated. Intermolecular hydrogen bonds at the central site of proflavine molecule are found to be stronger than the peripheral site. The hydrogen bond N–H???O for the hydrogen‐bonded dimer are indicated to be weakened in the excited states, since the excitation energy is increased slightly comparing to the monomer. Hydrogen bonds of PC–(H2O)2 trimer with the same type as the dimer are strengthened in the excited state, which is demonstrated by the decrease of the excited energies. Thus, hydrogen bond strengthening and weakening are observed to reveal site dependent feature in proflavine molecule. Furthermore, the hydrogen bond at central site induces the blue‐shift of the absorption spectrum, while the ones at peripheral site induce red‐shift. Hydrogen bonds with the same type at peripheral and central sites of proflavine molecule provide different effects on the photochemical and photophysical properties of proflavine.     

Simulation of water cluster assembly on a graphite surface

The assembly of small water clusters (H2O)n, n = 1-6, on a graphite surface is studied using a density functional tight-binding method complemented with an empirical van der Waals force correction, with confirmation using second-order M?ller-Plesset perturbation theory. It is shown that the optimized geometry of the water hexamer may change its original structure to an isoenergy one when interacting with a graphite surface in some specific orientation, while the smaller water cluster will maintain its cyclic or linear configurations (for the water dimer). The binding energy of water clusters interacting with graphite is dependent on the number of water molecules that form hydrogen bonds, but is independent of the water cluster size. These physically adsorbed water clusters show little change in their IR peak position and leave an almost perfect graphite surface.     

Order-Disorder Transition of Carboxyl Terminated Chains in Polydiacetylenes Vesicles Probed by Second Harmonic Generation and Two-Photon Fluorescence

To understand and control the interfacial properties of polydiacetylenes (PDAs) vesicles with π-conjugated backbone is very important for their colorimetric sensing of chemical and biological targets. In this work, we adopted 10, 12-pentacosadiynoic acid (PCDA) as the model molecule to prepare PDAs vesicles in aqueous solution with different forms (from monomer to blue-to-purple-to-red phase) by controlling the UV irradiation dose. The variations of the interfacial conformation of PDAs vesicles during chromatic transitions were inspected by the adsorption behaviors of probe molecules (4-(4-diethylaminostyry)-1-methylpyridinium iodide, D289) on vesicle surface with surface-specific second harmonic generation (SHG) and zeta potential measurements. Resonant SHG signal from D289 adsorbed on vesicle surface attenuated sharply, and the adsorption free energy as well as the corresponding two-photon fluorescence signal decreased slightly in chromatic transitions. While, the change in the surface density of the adsorbed D289 molecules for PDAs vesicles with different forms was relatively small as estimated from zeta potential measurements. The attenuation of the SHG intensity was thus attributed to the overall order-disorder transition and the changed orientation of D289 molecules caused by the gradual distortion of carboxyl head group driven by backbone perturbation.     

Surface Charges at the CaF2/Water Interface Allow Very Fast Intermolecular Vibrational‐Energy Transfer

We investigate the dynamics of water in contact with solid calcium fluoride, where at low pH, localized charges can develop upon fluorite dissolution. We use 2D surface‐specific vibrational spectroscopy to quantify the heterogeneity of the interfacial water (D₂O) molecules and provide information about the sub‐picosecond vibrational‐energy‐relaxation dynamics at the buried solid/liquid interface. We find that strongly H‐bonded OD groups, with a vibrational frequency below 2500 cm^?1, display very rapid spectral diffusion and vibrational relaxation; for weakly H‐bonded OD groups, above 2500 cm^?1, the dynamics slows down substantially. Atomistic simulations based on electronic‐structure theory reveal the molecular origin of energy transport through the local H‐bond network. We conclude that strongly oriented H‐bonded water molecules in the adsorbed layer, whose orientation is pinned by the localized charge defects, can exchange vibrational energy very rapidly due to the strong collective dipole, compensating for a partially missing solvation shell.     

Quasi chiral phase separation in a two-dimensional orientationally disordered system: 6-nitrospiropyran on Au(111)

The adsorption and chiral expression of 6-nitrospiropyran (SP6) molecules on a Au(111) surface are studied by scanning tunneling microscopy (STM) in combination with density functional theory (DFT) calculations. Both the chirality and the adsorption orientation of each adsorbed SP6 molecule are determined. The racemic mixture of SP6 enantiomers forms two-dimensional (2D) domains with same close packed positional orders but different internal orientational structures due to the random distribution of two adsorption orientations in each domain. However, all these orientationally disordered 2D domains undergo spontaneous quasi chiral phase separation; the 2D SP6 domains separate into 1D homochiral chains in which the SP6 molecules adopt two orientations randomly. This novel phenomenon is attributed to the preferential formation of the energetic favorable configurations with both the C-H...O weak hydrogen bonds and the pi-stacking of the two moieties of each SP6 molecule.