The extent of molecular orientation at liquid/vapor interface of pyridine and its alkyl derivatives by molecular dynamics simulation

In this paper, molecular dynamics simulation was performed to investigate the liquid∕vapor interfacial structure of neat polar liquids. Large-scale ensembles of liquid pyridine and its alkyl derivatives, 4-methylpyridine and 4-ethylpyridine, were simulated by classical molecular dynamics at 298 K. For the liquid system of low polarity, the surface density profile of the atoms meet exactly at the middle of interfacial region, and atoms of hydrophobic nature can be hardly discriminated from hydrophilic ones in either vapor or liquid sides. For a liquid system of high polarity, the density profile of atoms with different nature is highly discriminated all over the interfacial region, and as the polarity increases, a dense region of atomic density is clearly developed in the subsurface region. The recognized bivariate method was also used to study the molecular orientational distribution quantitatively. Orientational analysis of the three liquid systems indicates that the pyridine ring plane in the outmost surface tends to be vertical. Its tendency in the innermost interfacial region is parallel. The orientational states available to 4-ethylpyridine and pyridine are discriminated by predicting the possibility of a bisector-wise tumbling for the ring plane in pyridine and a side-wise tumbling in 4-ethylpyridine. The orientational distribution maps explain the trend of experimental surface tension and surface entropy. As the dipole moment of these liquids increases with the alkyl chain length, the surface structural profile changes from a regular definite one to a surface of complex atomic structure involving a dense phase near the interface. The development of dense region in alkyl derivatives is the result of segregation of molecules due to the alkyl group, which is captured and discriminated by molecular dynamics simulation even when the length of a short alkyl chain is increased by one carbon atom.     

Spectroscopic studies of solvated hydrogen and hydroxide ions at aqueous surfaces

Measuring the molecular properties of the surface of acidic and basic aqueous solutions is essential to understanding a wide range of important biological, chemical, and environmental processes on our planet. In the present studies, vibrational sum-frequency spectroscopy (VSFS) is employed in combination with isotopic dilution experiments at the vapor/water interface to elucidate the interfacial water structure as the pH is varied with HCl and NaOH. In acidic solutions, solvated proton species are seen throughout the interfacial region, and they alter the hydrogen bonding between water molecules in ways that reflect their depth in the interfacial region. At the higher frequencies of the OH stretch region, there is spectral evidence for solvated proton species residing in the topmost layers of the interfacial region. As reported in previous VSF studies, more strongly bound solvated proton species are observed at lower OH stretching frequencies. The solvated proton species that have stronger hydrogen bonding are similar in structure to those found in bulk acid solutions and likely reside somewhat deeper in the interfacial region. There is also evidence of OH stretching from solvated protons and relatively strong hydrogen bonding in the solvation sphere that is similar to other solvated ions. In contrast, water molecules solvating OH(-) ions show relatively weak hydrogen bonding and significantly less interfacial order. VSF spectra are acquired under multiple polarizations to provide crucial information for the interpretation of the spectra and for the determination of interfacial structure.     

Interfacial properties of semifluorinated alkane diblock copolymers

The liquid-vapor interfacial properties of semifluorinated linear alkane diblock copolymers of the form F(3)C(CF(2))(n-1)(CH(2))(m-1)CH(3) are studied by fully atomistic molecular dynamics simulations. The chemical composition and the conformation of the molecules at the interface are identified and correlated with the interfacial energies. A modified form of the Optimized Parameter for Liquid Simulation All-Atom (OPLS-AA) force field of Jorgensen and co-workers [J. Am. Chem. Soc. 106, 6638 (1984); 118, 11225 (1996); J. Phys. Chem. A 105, 4118 (2001)], which includes specific dihedral terms for H-F blocks-and corrections to the H-F nonbonded interaction, is used together with a new version of the exp-6 force field developed in this work. Both force fields yield good agreement with the available experimental liquid density and surface tension data as well as each other over significant temperature ranges and for a variety of chain lengths and compositions. The interfacial regions of semifluorinated alkanes are found to be rich in fluorinated groups compared to hydrogenated groups, an effect that decreases with increasing temperature but is independent of the fractional length of the fluorinated segments. The proliferation of fluorine at the surface substantially lowers the surface tension of the diblock copolymers, yielding values near those of perfluorinated alkanes and distinct from those of protonated alkanes of the same chain length. With decreasing temperatures within the liquid state, chains are found to preferentially align perpendicular to the interface, as previously seen.     

Probing molecular structure at interfaces for comparison with bulk solution behavior: water/2-propanol mixtures monitored by vibrational sum frequency spectroscopy

The orientation of the isopropyl group at the liquid/vapor interface in 2-propanol/water binary mixtures was studied by vibrational sum frequency spectroscopy. The CH(3) stretch modes of the two methyl groups were used to determine the molecule's orientation by employing a novel united atom approach to model the (CH(3))(2)X moiety. For this purpose, the changes in the molecular susceptibility of the isopropyl group stretches were derived in the laboratory frame as a function of the tilt and twist angles. The results indicated that the methyl groups lay down on the surface at low alcohol mole fraction and gradually twisted with increasing mole fraction. At the azeotrope, x(iso) = 0.68, one of the methyl groups aligned approximately parallel to the surface normal, whereas the other was nearly parallel with the liquid/vapor interface. When the mole fraction of 2-propanol was higher than 0.68, the orientation of 2-propanol remained almost constant. The change in the alcohol's orientation with 2-propanol mole fraction closely tracked changes in its bulk activity coefficient. Such results lead to a picture in which the surface structure and bulk properties of the system are closely linked.     

Adsorption and mixing behavior of ethers and alkanes at the solid/liquid interface

The behavior of binary mixtures of linear symmetrical ethers and alkanes adsorbed to a graphite surface from the bulk liquid mixtures is described on the basis of differential scanning calorimetry (DSC) data. Both the ethers and the alkanes are found to form solid monolayers when adsorbed from the liquid. In addition, the monolayer mixing behavior is addressed. The results indicate that there is good, essentially ideal, mixing in the monolayers for ethers and alkanes of the same overall chain length, where the chain length is equal to the total number of carbon and oxygen atoms in the molecule. However, a difference in chain length of more than one atom results in a variation of mixing behavior from nonideal mixing (for long pairs) to phase separation (for short pairs) on the graphite surface. Hence, we conclude that it is the relative chain lengths that control mixing behavior. The results are quantified using a regular solution model with a correction for preferential adsorption. The phase behavior of the mixed monolayers is also compared to the behavior of the bulk. Interestingly, we observe mixtures where the bulk and monolayer behavior are quite different, for example, phase separation in the bulk but essentially ideal mixing in the monolayer for mixtures of ethers and alkanes with the same chain lengths. At present, we attribute this mixing in the monolayer to dilution of the unfavorable ether oxygen-ether oxygen lone pair interactions by the coadsorbed alkanes. In addition, we find evidence for the preferential adsorption of the alkane over the ether. For example, heptane is preferentially adsorbed over dibutyl ether even though it contains two fewer atoms in the molecular chain. This contrasts with the preferential adsorption of alcohols over alkanes reported previously (Messe, L.; Perdigon, A.; Clarke, S. M.; Inaba, A.; Arnold, T. Langmuir 2005, 21, 5085-5093).     

Adsorption and reaction of CO2 and SO2 at a water surface

The orientation and hydrogen bonding of water molecules in the vapor/water interfacial region in the presence of SO2 and CO2 gas are examined using vibrational sum-frequency spectroscopy (VSFS) to gain insight into the adsorption and reactions of these gases in atmospheric aerosols. The results show that an SO2 surface complex forms when the water surface is exposed to an atmosphere of SO2 gas. Reaction of SO2 with interfacial water leads to other spectral changes that are examined by studying the VSF spectra and surface tension isotherms of several salts added to the aqueous phase, specifically NaHSO3, NaHCO3, Na2SO3, Na2CO3, Na2SO4, and NaHSO4. The results are compared with similar studies of CO2 adsorption and reaction at the surface. A weakly bound surface complex is not observed with CO2.     

Is there any group additive rules in the calculation of electron correlation energies of long straight chain alkane molecules?

According to the definition in the text, the correlation energy of 1s2C of carbon atoms, the primary and secondary C-H bonding electron pairs in some CH3, CH2 fragments and CH3(CH2)mCH3 (m=1-5) linear alkane molecules are calculated and analyzed. The transferability of the correlation energies of these electron pairs in the linear alkanes is investigated. The results indicate that the correlation energy of 1s2C is perfectly transferable in the respective methyl and methylene groups, while the correlation energies of the primary and secondary C-H bonding electron pairs are approximately transferable in methyl and methylene groups. The analysis of the results of group correlation energy shows that both of the correlation energies of methyl and methylene groups are transferable in these linear alkanes. The correlation energies of methylene group in CH3(CH2)mCH3 (m=1-5) molecules are slightly decreasing showing a converging trend to a "standard" methylene group in linear alkanes. The excellent fitting relationship between the total correlation energy and the number of methylene groups of the linear alkanes shows that the total correlation energy is a linear function of the number of methylene groups, which means that the total correlation energies of large linear alkanes can be reproduced and predicted by counting the numbers of methylene groups. In this way, total correlation energy of large linear alkane molecule can be approximately calculated using this simple group additive scheme with substantial saving in computational time.     

Air-liquid interfaces of aqueous solutions containing ammonium and sulfate: spectroscopic and molecular dynamics studies

Investigations of the air-liquid interface of aqueous salt solutions containing ammonium (NH(4)(+)) and sulfate (SO(4)(2-)) ions were carried out using molecular dynamics simulations and vibrational sum frequency generation spectroscopy. The molecular dynamics simulations show that the predominant effect of SO(4)(2-) ions, which are strongly repelled from the surface, is to increase the thickness of the interfacial region. The vibrational spectra reported are in the O-H stretching region of liquid water. Isotropic Raman and ATR-FTIR (attenuated total reflection Fourier transform infrared) spectroscopies were used to study the effect of ammonium and sulfate ions on the bulk structure of water, whereas surface sum frequency generation spectroscopy was used to study the effect of these ions on the interfacial structure of water. Analysis of the interfacial and bulk vibrational spectra reveal that aqueous solutions containing SO(4)(2-) perturb the interfacial water structure differently than the bulk and, consistent with the molecular dynamics simulations, reveal an increase in the thickness of the interfacial region.     

Alkane/Alcohol mixed monolayers at the solid/liquid interface

In this work, we present the behavior of solid monolayers of binary mixtures of alkanes and alcohols adsorbed on the surface of graphite from their liquid mixtures. We demonstrate that solid monolayers form for all the combinations investigated here. Differential scanning calorimetry (DSC) is used to identify the surface phase behavior of these mixtures, and elastic neutron incoherent scattering has been used to determine the composition of the mixed monolayers inferred by the calorimetry. The mixing behavior of the alcohol/alkane monolayer mixtures is compared quantitatively with alkane/alkane and alcohol/alcohol mixtures using a regular solution approach to model the incomplete mixing in the solid monolayer with preferential adsorption determining the surface composition. This analysis indicates the preferential adsorption of alcohols over alkanes of comparable alkyl chain length and even preferential adsorption of shorter alcohols over longer alkanes, which contrasts strongly with mixtures of alkane/alkane and alcohol/alcohol of different alkyl chain lengths where the longer homologue is always found to preferentially adsorb over the shorter. The alcohol/alkane mixtures are all found to phase separate to a significant extent in the adsorbed layer mixtures even when molecules are of a similar size. Again, this contrasts strongly with alkane/alkane and alcohol/alcohol mixtures where, although phase separation is found for molecules of significantly different size, good mixing is found for similar size species.     

Probing the protic ionic liquid surface using X-ray reflectivity

The structure of the free liquid surface of three protic ionic liquids, ethylammonium nitrate (EAN), propylammonium nitrate (PAN), and ethylammonium formate (EAF), has been elucidated using X-ray reflectivity. The results show all three liquids have an extended interfacial region, spanning at least five ion pairs, which can be divided into two parts. Adjacent to the gas phase are aggregates consisting of multiple cations and anions. Below this are layers oriented parallel to the macroscopic surface that are alternately enriched and depleted in cation alkyl chains and polar domains of cation ammonium groups and their anions, gradually decaying to the isotropic sponge-like bulk structure. The most pronounced layering is observed for PAN, driven by strong solvophobic interactions, while reduced hydrogen bonding in EAF results in the least structured and least extensive interfacial region.     

Self-assembly of normal alkanes on the Au (111) surfaces

The self-assembled monolayers (SAMs) of normal alkanes (n-C(n)H(2n+2)) with different carbon chain lengths (n=14-38) in the interfaces between alkane solutions (or liquids), and the reconstructed Au (111) surfaces have been systematically studied by means of scanning tunneling microscopy (STM). In contrast to previous studies, which concluded that some n-alkanes (n=18-26) can not form well-ordered structures on Au (111) surfaces, we observed SAM formations for all these n-alkanes without any exceptions. We find that gold reconstruction plays a critical role in the SAM formation. The alkane monolayers adopt a lamellar structure in which the alkane molecules are packed side-by-side, to form commensurate structures with respect to the reconstructed Au (111) surfaces. The carbon skeletons are found to lie flat on the surfaces, which is consistent with the infrared spectroscopic studies. Interestingly, we find that two-dimensional chiral lamellar structures form for alkanes with an even carbon number due to the specific packing of alkane molecules in a tilted lamella. Furthermore, we find that the orientation of alkane molecules deviates from the exact [011] direction, because of the intermolecular interactions among the terminal methyl groups of neighboring lamellae; this results in differences of molecular orientation between mirror structures of adjacent zigzag alkane lamellae. Structural models have been proposed, that shed new light on monolayer formation.     

Interfacial properties of free-standing poly(3-hexylthiophene) films

Using full atomistic classical molecular dynamics simulations, the interfacial properties of free-standing poly(3-hexylthiophene) (P3HT) films have been investigated. The orientations of different parts of the P3HT chain and the surface tensions of the films were calculated in a temperature range of 540 K-600 K. At the liquid/vacuum interface, the P3HT chain shows ordering by exposing hexyl groups at the interface, while the chain backbone lays flat with the thiophene ring preferentially tilt toward the surface. At the interface, the terminal methyl groups of hexyl side chains are in excess compared to the methylene groups or thiophene rings. The surface tension of P3HT in its melt state shows similar temperature dependence to that of polymers that have long alkyl side chains. The surface tension values are comparable to those polymers that expose methyl or methylene groups on the surface. The surface tension values determined for the melt state are lower than the experimental reported values for crystalline P3HT films, as expected.     

Local structure at the air/liquid interface of room-temperature ionic liquids probed by infrared-visible sum frequency generation vibrational spectroscopy: 1-alkyl-3-methylimidazolium tetrafluoroborates

The air/liquid interface of 1-alkyl-3-methylimidazolium tetrafluoroborates with the general formula [C(n)mim]BF(4) (n = 4-11) was studied using infrared-visible sum frequency generation (SFG) vibrational spectroscopy. The probability of the gauche defect per CH2-CH2 bond in the alkyl chain decreases as the number of carbon atoms in the alkyl chain increases. This observation suggests that the interaction between the alkyl chains is enhanced as the alkyl chain length becomes longer. The frequencies of the C-H stretching vibrational modes observed in the SFG spectra are higher than those of the corresponding peak positions observed in the infrared spectra of the bulk liquids. This shift is consistent with a structure in which the alkyl chain protrudes from the bulk liquid into the air. A local structure, which originates from the intermolecular interaction between the ionic liquid molecules, is proposed to explain these observations.     

Adsorption of organosulfur species at aqueous surfaces: molecular bonding and orientation

Four sulfur-containing compounds important to tropospheric chemistry have been examined at the vapor/H2O and vapor/D2O interfaces. These adsorbates, DMS, DMSO, DMSO2, and DMSO3, were studied by surface tension and vibrational sum-frequency spectroscopy (VSFS). Each adsorbate is surface active and each orients with the hydrophobic methyl groups pointed out of the plane of the interface. Their influence on the interfacial water structure is adsorbate dependent. Strong and weak interactions with surface water are observed as well as reorientation of subsurface water molecules, resulting in an increase in interfacial thickness.     

Molecular shape and orientational order. Effects in the energy of cavity formation in liquids

A recently proposed method for calculating the energy of cavity formation in liquids is presented in which the cavity formation process is described as work against the surface forces of the solvent, at the microscopic scale. The energy involved in the cavity formation process is, on the other hand, viewed as a strictly interaction potential energy and the reference cavity, which has the size and the shape of the space occupied by each molecule in the liquid, is considered as having short-range orientational order characteristic of the pure liquid. The method is successfully applied to binary alkane mixtures at infinite dilution whose components have different chemical structure (linear, cyclic and branched alkanes). The importance of the changes in the molecular order of the solute and the solvent occuring in the mixing process is emphasized.     

Structural inhomogeneity of interfacial water at lipid monolayers revealed by surface-specific vibrational pump-probe spectroscopy

We report vibrational lifetime measurements of the OH stretch vibration of interfacial water in contact with lipid monolayers, using time-resolved vibrational sum frequency (VSF) spectroscopy. The dynamics of water in contact with four different lipids are reported and are characterized by vibrational relaxation rates measured at 3200, 3300, 3400, and 3500 cm(-1). We observe that the water molecules with an OH frequency ranging from 3300 to 3500 cm(-1) all show vibrational relaxation with a time constant of T(1) = 180 ± 35 fs, similar to what is found for bulk water. Water molecules with OH groups near 3200 cm(-1) show distinctly faster relaxation dynamics, with T(1) < 80 fs. We successfully model the data by describing the interfacial water containing two distinct subensembles in which spectral diffusion is, respectively, rapid (3300-3500 cm(-1)) and absent (3200 cm(-1)). We discuss the potential biological implications of the presence of the strongly hydrogen-bonded, rapidly relaxing water molecules at 3200 cm(-1) that are decoupled from the bulk water system.     

Fully atomistic molecular‐mechanical model of liquid alkane oils: Computational validation

Fully atomistic molecular dynamics simulations were performed on liquid n‐pentane, n‐hexane, and n‐heptane to derive an atomistic model for middle‐chain‐length alkanes. All simulations were based on existing molecular‐mechanical parameters for alkanes. The computational protocol was optimized, for example, in terms of thermo‐ and barostat, to reproduce properly the properties of the liquids. The model was validated by comparison of thermal, structural, and dynamic properties of the normal alkane liquids to experimental data. Two different combinations of temperature and pressure coupling algorithms were tested. A simple differential approach was applied to evaluate fluctuation‐related properties with sufficient accuracy. Analysis of the data reveals a satisfactory representation of the hydrophobic systems behavior. Thermodynamic parameters are close to the experimental values and exhibit correct temperature dependence. The observed intramolecular geometry corresponds to extended conformations domination, whereas the intermolecular structure demonstrates all characteristics of liquid systems. Cavity size distribution function was calculated from coordinates analysis and was applied to study the solubility of gases in hexane and heptane oils. This study provides a platform for further in‐depth research on hydrophobic solutions and multicomponent systems. © 2014 Wiley Periodicals, Inc.     

A systematic molecular simulation study of ionic liquid surfaces using intrinsic analysis methods

In this paper, we apply novel intrinsic analysis methods, coupled with bivariate orientation analysis, to obtain a detailed picture of the molecular-level structure of ionic liquid surfaces. We observe pronounced layering at the interface, alternating non-polar with ionic regions. The outermost regions of the surface are populated by alkyl chains, which are followed by a dense and tightly packed layer formed of oppositely charged ionic moieties. We then systematically change the cation chain length, the anion size, the temperature and the molecular model, to examine the effect of each of these parameters on the interfacial structure. Increasing the cation chain length promotes orientations in which the chain is pointing into the vapor, thus increasing the coverage of the surface with alkyl groups. Larger anions promote a disruption of the dense ionic layer, increasing the orientational freedom of cations and increasing the amount of free space. The temperature had a relatively small effect on the surface structure, while the effect of the choice of molecular model was clearly significant, particularly on the orientational preferences at the interface. Our study demonstrates the usefulness of molecular simulation methods in the design of ionic liquids to suit particular applications.     

Formation and structure of self-assembled monolayers of alkanethiolates on palladium

The adsorption of n-alkanethiols onto polycrystalline thin films of palladium containing a strong (111) texture produces well-organized, self-assembled monolayers. The organization of the alkane chains in the monolayer and the nature of the bonding between the palladium and the thiol were studied by contact angle measurements, optical ellipsometry, reflection absorption infrared spectroscopy (RAIRS), and X-ray photoelectron spectroscopy (XPS). The XPS data reveals that a compound palladium-sulfide interphase is present at the surface of the palladium film. The RAIR spectra, ellipsometry data, and wetting properties show that the palladium-sulfide phase is terminated with an organized, methyl-terminated monolayer of alkanethiolates. The local molecular environment of the alkane chains transitions from a conformationally disordered, liquidlike state to a mostly all-trans, crystalline-like structure with increasing chain length (n = 8-26). The intensities and dichroism of the methylene and methyl stretching modes support a model for the average orientation of an ensemble of all-trans-conformer chains with a tilt angle of approximately 14-18 degrees with respect to the surface normal and a twist angle of the CCC plane relative to the tilt plane of approximately 45 degrees. The SAMs are stable in air, although the sulfur present at the surface oxidizes in air over a period of 2-5 days at room temperature. The differences in chain organization between SAMs formed by microcontact printing and by solution deposition are also examined by RAIRS and XPS.     

Molecular level structures of poly(n-alkyl methacrylate)s with different side chain lengths at the polymer/air and polymer/water interfaces

Sum frequency generation (SFG) vibrational spectroscopy has been successfully applied to study molecular structures of several poly(n-alkyl methacrylate)s (PAMAs) with different side chain lengths at the PAMA/air and PAMA/water interfaces. We have observed that the ester side chains from all PAMAs always dominate the interface, but the orientation information of the methyl end group on the side chains varies, depending on the length of the side chain. The contributions from methylene groups on the side chains have been evaluated, and the surface structures have been related to the surface tension of these polymers. Different water restructuring behaviors have been observed for different PAMAs. This phenomenon and its reversibility are strongly dependent on the glass transition temperature of each polymer, which is influenced by the side chain length. Detailed data fitting and analysis has been discussed.