Vapor+liquid equilibrium of water, carbon dioxide, and the binary system, water+carbon dioxide, from molecular simulation

NVT- and NpT-Gibbs ensemble Monte Carlo (GEMC) simulations were applied to describe the vapor–liquid equilibrium of water (between 323 and 573 K), carbon dioxide (between 230 and 290 K) and their binary mixtures (between 348 and 393 K). The properties of supercritical carbon dioxide were determined between 310 and 520 K by NpT-MC simulations. Literature data for the effective pair potentials (for water: the SPC-, SPC/E-, and TIP4P potential models; for carbon dioxide: the EPM2 potential model) were used to describe the properties of the pure substances. The vapor pressures of water and carbon dioxide are calculated. For water, the SPC- and TIP4P models give superior results for the vapor pressure when compared to the SPC/E model. The vapor–liquid equilibrium of the binary mixture, carbon dioxide–water, was predicted using the SPC- as well as the TIP4P model for water and the EPM2 model for carbon dioxide. The interactions between carbon dioxide and water were estimated from the pair potentials of the pure components using common mixing rules without any adjustable binary parameter. Agreement of the predicted data for the compositions of the coexisting phases in vapor–liquid equilibrium and experimental results is observed within the statistical uncertainties of the simulation results in the investigated range of state, i.e. at pressures up to about 20 MPa.     

Thermodynamic properties of methane/water interface predicted by molecular dynamics simulations

Molecular dynamics simulations have been performed to examine the thermodynamic properties of methane/water interface using two different water models, the TIP4P/2005 and SPC/E, and two sets of combining rules. The density profiles, interfacial tensions, surface excesses, surface pressures, and coexisting densities are calculated over a wide range of pressure conditions. The TIP4P/2005 water model was used, with an optimized combining rule between water and methane fit to the solubility, to provide good predictions of interfacial properties. The use of the infinite dilution approximation to calculate the surface excesses from the interfacial tensions is examined comparing the surface pressures obtained by different approaches. It is shown that both the change of methane solubilities in pressure and position of maximum methane density profile at the interface are independent of pressure up to about 2 MPa. We have also calculated the adsorption enthalpies and entropies to describe the temperature dependency of the adsorption.     

Tailoring porous silica films through supercritical carbon dioxide processing of fluorinated surfactant templates

The tailoring of porous silica thin films synthesized using perfluoroalkylpyridinium chloride surfactants as templating agents is achieved as a function of carbon dioxide processing conditions and surfactant tail length and branching. Well-ordered films with 2D hexagonal close-packed pore structure are obtained from sol-gel synthesis using the following cationic fluorinated surfactants as templates: 1-(3,3,4,4,5,5,6,6,7,7,8,8,8-tridecafluoro-octyl)pyridinium chloride (HFOPC), 1-(3,3,4,4,5,5,6,6,7,8,8,8-dodecafluoro-7-trifluoromethyl -octyl)pyridinium chloride (HFDoMePC), and 1-(3,3,4,4,5,5,6,6,7,7,8,8,9,9,10,10,10-heptadecafluoro-decyl)pyridinium chloride (HFDePC). Processing the sol-gel film with CO2 (69-172 bar, 25 and 45 degrees C) immediately after coating results in significant increases in pore diameter relative to the unprocessed thin films (increasing from 20% to 80% depending on surfactant template and processing conditions). Pore expansion increases with CO2 processing pressure, surfactant tail length, and surfactant branching. The varying degree of CO2 induced expansion is attributed to the solvation of the "CO2-philic" fluorinated tail and is interpreted from interfacial behavior of HFOPC, HFDoMePC, and HFDePC at the CO2-water interface.     

Surface reactions of carbon dioxide at the adsorbed water-iron oxide interface

Despite the fact that carbon dioxide is an abundant atmospheric gas with profound environmental implications, there is little information on the reaction of carbon dioxide at the adsorbed water-oxide interface. In this study, the chemistry of carbon dioxide at the adsorbed water-iron oxide interface is investigated with FTIR spectroscopy. As shown here, the thin water layer on the iron oxide surface plays an important role in the surface chemistry of carbon dioxide. In particular, adsorbed water enhances CO(2) uptake, undergoes isotope exchange with CO(2) in O(18)-labeled experiments, and influences the chemical nature of the predominant adsorbed product on the surface from bicarbonate to carbonate. The resultant thin water film is acidic in nature from the reaction of CO(2). The IR spectrum recorded of adsorbed carbonate at the adsorbed water-iron oxide interface is remarkably similar to that at the bulk liquid water-iron oxide interface. Since reactions in thin water films estimated to be approximately 2 layers will play a role in a number of environmental processes, it is essential to understand the chemistry of these "wet" interfaces with atmospheric gases.     

DLVO-based estimates of adsorbed water film thicknesses in geologic CO2 reservoirs

When supercritical carbon dioxide (scCO(2)) is injected into deep subsurface reservoirs, much of the affected volume consists of pores containing both water and scCO(2), with water films remaining as the mineral-wetting phase. Although water films can affect multiphase flow and mediate reactions at mineral surfaces, little is known about how film thicknesses depend on system properties. Here, the thicknesses of water films were estimated on the basis of considerations of capillary pressure needed for the entry of CO(2) and disjoining pressures in films resulting from van der Waals and electric double-layer interactions. Depth-dependent CO(2) and water properties were used to estimate Hamaker constants for water films on silica and smectite surfaces under CO(2) confinement. Dispersion interactions were combined with approximate solutions to the electric double layer film thickness-pressure relationship in a Derjaguin-Landau-Verwey-Overbeek (DLVO) analysis, with CO(2) as the confining fluid. Under conditions of elevated pressure, temperature, and salinity commonly associated with CO(2) sequestration, adsorbed water films in reservoir rock surfaces are typically predicted to be less than 10 nm in thickness. Decreased surface charge of silica under the acidic pH of CO(2)-equilibrated water and elevated salinity is predicted to compress the electric double layer substantially, such that the dispersion contribution to the film thickness is dominant. Relative to silica, smectite surfaces are predicted to support thicker water films under CO(2) confinement because of greater electrostatic and dispersion stabilization.     

Solvation forces between silica bodies in supercritical carbon dioxide

We report Monte Carlo simulations of the solvation pressure between two planar surfaces, which represent the interface of spherical silica nanoparticles in supercritical carbon dioxide. Carbon dioxide (CO2) was modeled as an atomistic dumbbell or a spherical Lennard-Jones particle. The interaction between CO2 molecules and silica surfaces was characterized by the standard Steele potential with energetic heterogeneities representing the hydrogen bonds. The parameters for the solid-fluid interaction potentials were obtained by fitting our simulations to the experimental isotherms of CO2 sorption on mesoporous siliceous materials. We studied the dependence of the solvation force on the distance between planar silica surfaces at T = 318 K, at equilibrium bulk pressures p(bulk) ranging from 69 to 200 atm. At 69 atm, we observed a long-range attraction between the two surfaces, and it vanished when the pressure was increased to 102 and then 200 atm. The results obtained with different fluid models were consistent with each other. According to our observations, energetic heterogeneities of the surface have negligible influence on the solvation pressure. Using the Derjaguin approximation, we calculated the solvation forces between spherical silica nanoparticles in supercritical CO2 from the solvation pressures between the planar surfaces.     

Enhanced stabilization of reverse micelles by compressed CO2

The effect of compressed CO2 on the solubilization capacity of water in reverse micelles of sodium bis(2-ethylhexyl) sulfosuccinate (AOT) in longer chain n-alkanes was studied at different temperatures and pressures. It was found that the amount of solubilized water is increased considerably by CO2 in a suitable pressure range. The suitable CO2 pressure range in which the solubilization capacity of water could be enhanced decreased with increasing W0 (water-to-AOT molar ratio). The microenvironments in the CO2-stabilized reverse micelles were investigated by UV/Vis adsorption spectroscopy with methyl orange (MO) as probe. The mechanism by which the reverse micelles are stabilized by CO2 is discussed in detail. The main reason is likely to be that CO2 has a much smaller molecular volume than the n-alkane solvents studied in this work. Therefore, it can penetrate the interfacial film of the reverse micelles and stabilize them by increasing the rigidity of the micellar interface and thus reducing the attractive interaction between the droplets. However, if the CO2 pressure is too high, the solvent strength of the solvents is reduced markedly, and this induces phase separation in the micellar solution.     

Capillary waves at the liquid-vapor interface and the surface tension of water

Capillary waves occurring at the liquid-vapor interface of water are studied using molecular dynamics simulations. In addition, the surface tension, determined thermodynamically from the difference in the normal and tangential pressure at the liquid-vapor interface, is compared for a number of standard three- and four-point water models. We study four three-point models (SPC/E, TIP3P, TIP3P-CHARMM, and TIP3P-Ew) and two four-point models (TIP4P and TIP4P-Ew). All of the models examined underestimate the surface tension; the TIP4P-Ew model comes closest to reproducing the experimental data. The surface tension can also be determined from the amplitude of capillary waves at the liquid-vapor interface by varying the surface area of the interface. The surface tensions determined from the amplitude of the logarithmic divergence of the capillary interfacial width and from the traditional thermodynamic method agree only if the density profile is fitted to an error function instead of a hyperbolic tangent function.     

聚苯乙烯-聚(4-乙烯基吡啶)嵌段共聚物LB膜结构形态的研究

采用氯仿作为铺展溶剂,将嵌段共聚物聚苯乙烯-聚(4-乙烯基吡啶)(PS-b-P4VP)稀溶液铺展于空气与水界面上,利用Langmuir-Blodgett(LB)膜技术转移至固体基底.研究了不同的嵌段比、表面压和小分子1-芘丁酸(PBA)的加入对嵌段共聚物气液界面聚集组装的影响.研究发现随着亲水段(P4VP)的增加,聚集组装结构由纳米片状、带状转变成纳米条状、纳米点状结构.表面压对纯PS-b-P4VP聚集组装产生影响,表面压增大,组装体排列紧密;随着表面压的继续增大,单层聚集结构遭到破坏,发生堆叠.加入PBA小分子后,PBA与PS-b-P4VP形成氢键,形态发生明显变化,原来的片状结构转变为条状或点状结构.     

Hydrolysis reaction analysis of L-alpha-distearoylphosphatidylcholine monolayer catalyzed by phospholipase A2 with polarization-modulated infrared reflection absorption spectroscopy

The hydrolysis reaction of L-alpha-distearoylphosphatidylcholine (DSPC) monolayers catalyzed by phospholipase A2 (PLA(2)) has been studied using polarization-modulated infrared reflection absorption spectroscopy (PM-IRRAS) with film balance measurements. The PM-IRRAS analysis provides quantitative information about the reaction efficiency at different surface pressures. It was found that the reaction efficiency of L-DSPC monolayer hydrolysis catalyzed by PLA(2) decreased with increasing surface pressure. At zero pressure (lift-off point), the hydrolysis reaction efficiency has the highest value of 45%. Increasing surface pressure leads to the decrease of the hydrolysis efficiency. Since the surface pressure is above 20 mN/m, the hydrolysis reaction nearly stopped. PM-IRRAS technique provides a powerful means to study the hydrolysis process catalyzed by phospholipase A2 at the air/water interface.     

Contact angle of water on polystyrene thin films: effects of CO(2) environment and film thickness

We examine the contact angle of water droplets on polystyrene (PS) thin films of varying thicknesses supported by silicon wafers under both air and pressurized carbon dioxide (CO2) environments. At 23 degrees C, the contact angle is found to increase upon increasing CO2 pressure in the vapor regime and then levels off in the liquid CO2 regime. A macroscopic model based on Young's equation and the geometric-mean method for interfacial tensions, and long-range van der Waals interactions, correctly predicts the trends and the magnitude of the contact angle dependence on pressure, although deviations occur at high CO2 activities. The contact angle was also found to depend on film thickness, h, when h was comparable to or smaller than 50 nm. Specifically, the contact angle decreases with decreasing PS film thickness. This behavior could be accounted for with the use of a model that incorporates the effects of film thickness, CO2 pressure, and the long-range van der Waals potential.     

Structural and shear characteristics of adsorbed sodium caseinate and monoglyceride mixed monolayers at the air-water interface

The structural and shear characteristics of mixed monolayers formed by an adsorbed Na-caseinate film and a spread monoglyceride (monopalmitin or monoolein) on the previously adsorbed protein film have been analyzed. Measurements of the surface pressure (pi)-area (A) isotherm and surface shear viscosity (eta(s)) were obtained at 20 degrees C and at pH 7 in a modified Wilhelmy-type film balance. The structural and shear characteristics of the mixed films depend on the surface pressure and on the composition of the mixed film. At surface pressures lower than the equilibrium surface pressure of Na-caseinate (at pipi(e)(CS) have important repercussions on the shear characteristics of the mixed films.     

水合物生成条件下甲烷在水中的吸附

采用悬滴法系统地测定了温度274.2 ~ 282.2 K、压力0.1 ~ 10.1 MPa下甲烷/纯水间界面张力。实验结果表明在恒定温度下界面张力随压力的增加而增大。在高压条件下,压力对界面张力有很大的影响。不同温度和压力下计算出的甲烷在水中的表面过剩浓度结果表明,压力越高,温度越低,甲烷在水溶液中的吸附浓度越高。同时,计算出的甲烷在水溶液中的表面吸附自由能结果表明,在水合物生成条件下,甲烷在水中的吸附比298.2 K更容易。     

Evolution of intermolecular structure and dynamics in supercritical carbon dioxide with pressure: an ab initio molecular dynamics study

The effect of pressure on supercritical carbon dioxide (scCO2) has been characterized by using Car-Parrinello molecular dynamics simulations. Structural and dynamical properties along an isotherm of 318.15 K and at pressures ranging from 190 to 5000 bar have been obtained. Intermolecular pair correlation functions and three-dimensional atomic probability density map calculations indicate that the local environment of a central CO2 molecule becomes more structured with increasing pressure. The closest neighbors are predominantly oriented in a distorted T-shaped geometry while neighbors separated by larger distances are likely oriented in a slipped parallel arrangement. The structure of scCO2 at high densities has been compared with that of crystalline CO2. The probability distributions of intramolecular distances narrow down with increasing pressure. A marginal but non-negligible effect of pressure on the instantaneous intramolecular OCO angle is observed, lending credence to the idea that intermolecular interactions between CO2 molecules in an inhomogeneous near neighbor environment could contribute to the observed instantaneous molecular dipole moment. The extent of deviation from a perfect linear geometry of the carbon dioxide molecule decreases with increasing pressure. Time constants derived from reorientational time correlation functions of the molecular backbone compare well with experimental data. Within the range of thermodynamic conditions explored here, no significant changes are observed in the frequencies of intramolecular vibrational modes. However, a blue shift is observed in the low-frequency cage rattling mode with increasing pressure.     

CO-induced embedding of Pt adatoms in a partially reduced FeO(x) film on Pt(111)

The reduction of a single-layer FeO film grown on Pt(111) by CO at elevated pressures and temperatures has been studied through an interplay of scanning tunneling microscopy, ambient-pressure X-ray photoelectron spectroscopy, and density functional theory calculations. Exposure of the FeO thin film to CO at pressures between 1 and 30 Torr and temperatures between 500 and 530 K leads to formation of a honeycomb-structured Fe(3)O(2) film with hollow sites occupied by single Pt atoms extracted from the substrate surface. The formation of these adatoms is driven by an increase in CO adsorption energy. In addition, the structure incorporates undercoordinated Fe centers, which are proposed to have substantial effects on the catalytic properties of the surface.     

Nanobubbles at the interface between water and a hydrophobic solid

A very thin layer (5-80 nm) of gas phase, consisting of discrete bubbles with only about 40 000 molecules, is quite stable at the interface between a hydrophobic solid and water. We prepare this gas phase from either ambient air or from CO(2)(g) through a solvent exchange method reported previously. In this work, we examine the interface using attenuated total internal reflection infrared spectroscopy. The presence of rotational fine structure in the spectrum of CO(2) and D(2)O proves that molecules are present in the gas phase at the interface. The air bubbles are stable for more than 4 days, whereas the CO(2) bubbles are only stable for 1-2 h. We determine the average gas pressure inside the CO(2) bubbles from the IR spectrum in two ways: from the width of the rotational fine structure (P(gas) < 2 atm) and from the intensity in the IR spectrum (P(gas) = 1.1 +/- 0.4 atm). The small difference in gas pressure between the bubbles and the ambient (1 atm) is consistent with the long lifetime. The dimensions and curvature of a set of individual bubbles was determined by atomic force microscopy. The pressures of individual bubbles calculated from the measured curvature using the Laplace equation fall into the range P(gas) = 1.0-1.7 atm, which is concordant with the average pressure measured from the IR spectrum. We believe that the difference in stability of the CO(2) bubbles and the air bubbles is due to a combination of the much lower pressure of CO(2) in the atmosphere and the greater solubility of CO(2) in water, compared to N(2) and O(2). As expected, smaller bubbles have a shorter average lifetime than larger bubbles, and the average pressure and the curvature of individual bubbles decreases with time. Surface plasmon resonance measurements provide supporting evidence that the film is in the gas state: the thin film has a lower refractive index than water, and there are few common contaminants that satisfy this condition. Interfacial gas bubbles are not ubiquitous on hydrophobic solids: bubble-free and bubble-decorated hydrophobic interfaces can be routinely prepared.     

Study of the diclofenac/phospholipid interactions with liposomes and monolayers

The interaction of diclofenac sodium (SD) with soya phosphatidylcholine (SPC) has been studied with floating Langmuir monolayers and liposomes. SD was either introduced into the subphase of SPC monolayers or co-spread with SPC on an aqueous subphase. In both cases, SD caused the surface pressure isotherm to become more expanded, thus demonstrating the affinity between SD and SPC. The incorporation of SD caused SPC liposomes to have a decreased diameter according to light scattering experiments. When SPC liposomes were injected into an aqueous subphase, their destruction yielding surface-active monomers could be monitored by changes in surface pressure. SD-loaded liposomes displayed a much faster kinetics when the surface density of surface-active monomers was plotted against time, with rate constants increasing significantly with the SD concentration. The kinetic profile can be quantitatively analyzed by plotting ln[1 - (gamma/gamma infinity)] versus t1/2.     

Characterisation of spread monolayers of functionalised polyisobutylenes with various chain architectures

Nine polyisobutylene (PIB) samples, consisting of three chain architectures (linear and branched), each of three molecular weights, were each derivatised at one, two or three chain ends with a hydrophilic head group based on diethyl azodicarboxylate. The surface pressure–area isotherms for spread monolayers of these polymers at the air/water interface were determined using a Langmuir trough, and the thicknesses of the films determined by ellipsometry. The isotherms showed that the area per molecule increased, whereas the collapse pressure decreased, with increasing molecular weight. The area per molecule and collapse pressures both increased with functionality, suggesting that the PIBs adsorb with all of the head groups at the interface. This was confirmed by the ellipsometry data which showed that the film thickness decreased with increasing functionality for similar molecular weight PIBs. Compared to other head groups which have been used to functionalise PIB, the azo head group has a much weaker affinity for water, leading to higher collapse areas and lower collapse pressures.     

Langmuir monolayer and Langmuir–Blodgett films of amphiphilic triblock copolymers with water-soluble middle block

Langmuir monolayers and Langmuir–Blodgett (LB) film morphology of amphiphilic triblock copolymers are studied using surface pressure-area measurements and atomic force microscopy (AFM), respectively. The triblock copolymers are composed of long water-soluble poly(ethylene oxide) (PEO) chains as middle block with very short poly(perfluorohexylethyl methacrylate) (PFMA) end blocks. The surface pressure-area isotherms show phase transitions in the brush regime. This phase transition is due to a rearrangement of PFMA block at the air–water interface. It becomes more significant with increasing PFMA content in the copolymer. LB films transferred at low surface pressures from the air–water interface to hydrophilic silicon substrates show surface micelles in the size range of 50–100 nm. A typical crystalline morphology of the corresponding PEO homopolymer is observed in LB films of copolymers with very short PFMA blocks, transferred in the brush region at high surface pressure. This crystallization is hindered with increasing PFMA content in the copolymer.     

Adsorption of carbon dioxide of 1-site and 3-site models in pillared clays: a Gibbs ensemble Monte Carlo simulation

Adsorption behavior of carbon dioxide confined in pillared clays is analyzed by using constant pressure Gibbs ensemble Monte Carlo (GEMC) method. In our simulation, 1-site and 3-site models are used to represent carbon dioxide. At the 1-site model, carbon dioxide is described as a Lennard-Jones (LJ) sphere, while at the 3-site model, carbon dioxide is modeled as a three-sites linear chain represented by EPM2 potential considering the quadrapolar effect. The potential model from Yi et al. for pillared clays is used to emphasize its quasi two-dimensional structure. Comparing the calculated results from the 1-site and the 3-site models at T=228.15 and 258.15 K, we observe that the adsorption amount from the two models is the same basically. However, the local density presents a significant difference, because the shoulder in the main peak near the wall from 3-site model can reflect the orientation of carbon dioxide. Accordingly, in the systematical investigation to explore the effect of porosity and pore width on the adsorption of carbon dioxide in pillared clays, the 3-site model was only used. We observe that for a narrow pore of H=1.02 nm, each isotherm shape displays type I curve, suggesting that it is not inflected by the porosity. However, for the larger pores of H=1.70 and H=2.38 nm, the increase of the porosity alters the shape of adsorption isotherms from a simple linear relation to the first order jump, indicating that the porosity is of very important factor to affect adsorption and phase behavior of fluids confined in pillared clays. The excess adsorptions of carbon dioxide at supercritical temperatures of T=323.15 and 348.15 K are also investigated. We find that the maximum exists for each excess isotherm, and the optimal pressure corresponding to the maximum increases with the pore width. However, the porosity has no significant effect on the optimal pressure.