Hydration structure of water confined between mica surfaces

We report further molecular dynamics simulations on the structure of bound hydration layers under extreme confinement between mica surfaces. We find that the liquid phase of water is maintained down to 2 monolayer (ML) thick, whereas the structure of the K(+) ion hydration shell is close to the bulk structure even under D = 0.92 nm confinement. Unexpectedly, the density of confined water remains approximately the bulk value or less, whereas the diffusion of water molecules decreases dramatically. Further increase in confinement leads to a transition to a bilayer ice, whose density is much less than that of ice Ih due to the formation of a specific hydrogen-bonding network.     

Electrofreezing of confined water

We report results from molecular dynamics simulations of the freezing transition of TIP5P water molecules confined between two parallel plates under the influence of a homogeneous external electric field, with magnitude of 5 V/nm, along the lateral direction. For water confined to a thickness of a trilayer we find two different phases of ice at a temperature of T=280 K. The transformation between the two, proton-ordered, ice phases is found to be a strong first-order transition. The low-density ice phase is built from hexagonal rings parallel to the confining walls and corresponds to the structure of cubic ice. The high-density ice phase has an in-plane rhombic symmetry of the oxygen atoms and larger distortion of hydrogen bond angles. The short-range order of the two ice phases is the same as the local structure of the two bilayer phases of liquid water found recently in the absence of an electric field [J. Chem. Phys. 119, 1694 (2003)]. These high- and low-density phases of water differ in local ordering at the level of the second shell of nearest neighbors. The results reported in this paper, show a close similarity between the local structure of the liquid phase and the short-range order of the corresponding solid phase. This similarity might be enhanced in water due to the deep attractive well characterizing hydrogen bond interactions. We also investigate the low-density ice phase confined to a thickness of 4, 5, and 8 molecular layers under the influence of an electric field at T=300 K. In general, we find that the degree of ordering decreases as the distance between the two confining walls increases.     

Breakdown of hydration repulsion between charged surfaces in aqueous Cs+ solutions

Using a surface force balance, we have measured the normal and shear forces between mica surfaces across aqueous caesium salt solutions (CsNO(3) and CsCl) up to 100 mM concentrations. In contrast to all other alkali metal ions at these concentrations, we find no evidence of hydration repulsion between the mica surfaces on close approach: the surfaces appear to be largely neutralized by condensation of the Cs ions onto the charged lattice sites, and are attracted on approach into adhesive contact. The contact separation at adhesion indicates that the condensed Cs ions protrude by 0.3 +/- 0.2 nm from each surface, an observation supported both by the relatively weak adhesion energies between the surfaces, and the relatively weak frictional yield stress when they are made to slide past each other. These observations show directly that the hydration shells about the Cs(+) ions are removed as the ions condense into the charged surface lattice. This effect is attributed to the low energies-resulting from their large ionic radius-required for dehydration of these ions.     

Adding salt to an aqueous solution of t-butanol: is hydrophobic association enhanced or reduced?

Recent neutron scattering experiments on aqueous salt solutions of amphiphilic t-butanol by Bowron and Finney [Phys. Rev. Lett. 89, 215508 (2002); J. Chem. Phys. 118, 8357 (2003)] suggest the formation of t-butanol pairs, bridged by a chloride ion via O-H...Cl- hydrogen bonds, leading to a reduced number of intermolecular hydrophobic butanol-butanol contacts. Here we present a joint experimental/theoretical study on the same system, using a combination of molecular dynamics (MD) simulations and nuclear magnetic relaxation measurements. Both MD simulation and experiment clearly support the more classical scenario of an enhanced number of hydrophobic contacts in the presence of salt, as it would be expected for purely hydrophobic solutes. [T. Ghosh et al., J. Phys. Chem. B 107, 612 (2003)]. Although our conclusions arrive at a structurally completely distinct scenario, the molecular dynamics simulation results are within the experimental error bars of the Bowron and Finney data.     

Normal and frictional forces between surfaces bearing polyelectrolyte brushes

Normal and shear forces were measured as a function of surface separation, D, between hydrophobized mica surfaces bearing layers of a hydrophobic-polyelectrolytic diblock copolymer, poly(methyl methacrylate)- block-poly(sodium sulfonated glycidyl methacrylate) copolymer (PMMA- b-PSGMA). The copolymers were attached to each hydrophobized surface by their hydrophobic PMMA moieties with the nonadsorbing polyelectrolytic PSGMA tails extending into the aqueous medium to form a polyelectrolyte brush. Following overnight incubation in 10 (-4) w/v aqueous solution of the copolymer, the strong hydrophobic attraction between the hydrophobized mica surfaces across water was replaced by strongly repulsive normal forces between them. These were attributed to the osmotic repulsion arising from the confined counterions at long-range, together with steric repulsion between the compressed brush layers at shorter range. The corresponding shear forces on sliding the surfaces were extremely low and below our detection limit (+/-20-30 nN), even when compressed down to a volume fraction close to unity. On further compression, very weak shear forces (130 +/- 30 nN) were measured due to the increase in the effective viscous drag experienced by the compressed, sliding layers. At separations corresponding to pressures of a few atmospheres, the shearing motion led to abrupt removal of most of the chains out of the gap, and the surfaces jumped into adhesive contact. The extremely low frictional forces between the charged brushes (prior to their removal) is attributed to the exceptional resistance to mutual interpenetration displayed by the compressed, counterion-swollen brushes, together with the fluidity of the hydration layers surrounding the charged, rubbing polymer segments.     

Forces between mica surfaces, prepared in different ways, across aqueous and nonaqueous liquids confined to molecularly thin films

We have measured normal and lateral interactions across a range of different liquids between mica surfaces using a surface force balance (SFB). The mica surfaces were prepared either by melt cutting using Pt wire and standard procedures in our laboratories or by tearing sheets (that had not been exposed to Pt) off from a freshly cleaved sheet of mica. AFM micrographs revealed the substantial absence of Pt nanoparticles on the melt cut and torn-off mica surfaces. Normal-force versus surface-separation (D) profiles and shear force versus D measurements for purified water (no added salt), for concentrated aqueous NaCl solutions, and for cyclohexane revealed that in all cases the behavior of the highly confined liquids between melt-cut and between torn-off mica sheets was identical within experimental scatter. These results demonstrate directly that interactions measured between melt-cut mica surfaces as routinely prepared using established procedures in our laboratories and in other laboratories are free of the effect of any Pt contamination.     

A facile approach for assembling lipid bilayer membranes on template-stripped gold

Lipid vesicles are designed with functional chemical groups to promote vesicle fusion on template-stripped gold (TS Au) surfaces that does not spontaneously occur on unfunctionalized Au surfaces. Three types of vesicles were exposed to TS Au surfaces: (1) vesicles composed of only 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipids; (2) vesicles composed of lipid mixtures of 2.5 mol % of 1,2-distearoyl-sn-glycero-3-phosphoethanolamine-N-poly(ethylene glycol)-2000-N-[3-(2-pyridyldithio)propionate] (DSPE-PEG-PDP) and 97.5 mol % of POPC; and (3) vesicles composed of 2.5 mol % of 1,2-distearoyl-sn-glycero-3-phosphoethanolamine-N-[methoxy(poly(ethylene glycol))-2000] (DSPE-PEG) and 97.5 mol % POPC. Atomic force microscopy (AFM) topography and force spectroscopy measurements acquired in a fluid environment confirmed tethered lipid bilayer membrane (tLBM) formation only for vesicles composed of 2.5 mol % DSPE-PEG-PDP/97.5 mol % POPC, thus indicating that the sulfur-containing PDP group is necessary to achieve tLBM formation on TS Au via Au-thiolate bonds. Analysis of force-distance curves for 2.5 mol % DSPE-PEG-PDP/97.5 mol % POPC tLBMs on TS Au yielded a breakthrough distance of 4.8 ± 0.4 nm, which is about 1.7 nm thicker than that of POPC lipid bilayer membrane formed on mica. Thus, the PEG group serves as a spacer layer between the tLBM and the TS Au surface. Fluorescence microscopy results indicate that these tLBMs also have greater mechanical stability than solid-supported lipid bilayer membranes made from the same vesicles on mica. The described process for assembling stable tLBMs on Au surfaces is compatible with microdispensing used in array fabrication.     

Nanorheology of dioctyl phthalate confined between surfaces coated with long alkyl chains

To shed light on the mechanism related to the high viscosity of a precipitated calcium carbonate (PCC) dispersion in dioctyl phthalate (DOP), the viscosity of DOP in a nanometer space was investigated using the shear resonance measurement. We used mica surfaces modified with dioctadecyldimethylammonium bromide (DODA) as a model of a PCC surface which bears long alkyl chains and bare mica surfaces as a reference. We found that a resonance peak for DOP of high water content (1164 ppm) confined between the DODA-modified surfaces gradually decreased upon approach from a distance of 320.6 nm and disappeared at ca. 57 nm. This indicated that the viscosity of DOP increased with decreasing distance, and a highly ordered state appeared at a large separation of 57 nm. This highly ordered state, however, cannot be observed for DOP of low water content (469 ppm) confined between the DODA-modified surfaces. The resonance peak for DOP of high water content (1164 ppm) between bare mica surfaces gradually decreased upon approach from a distance of 358.4 nm and disappeared at 10.7 nm, which is much shorter than the value for the DODA-modified surfaces. These results could be correlated with the high macroscopic viscosity of a dispersion of PCC coated with long alkyl chains in DOP. We propose the interpenetration of the long alkyl chains and the hydrogen-bonding network of DOP molecules mediated by water molecules as a plausible model for interpreting the high viscosity of coated PCC dispersions.     

Selective adsorption of poly(ethylene oxide) onto a charged surface mediated by alkali metal ions

Using a surface force balance, we have measured normal and shear interactions between mica surfaces across pure water and across 0.1 M aqueous solutions of LiNO3, NaNO3, KNO3, and CsNO3, both prior to adding polymer and following addition of 1.5 x 10(-4) w/w poly(ethylene oxide) (PEO, Mw = 170 kD) and overnight incubation. Our results reveal that while the PEO adsorbs strongly from the KNO3 and CsNO3 solutions, unexpectedly it does not adsorb at all from the LiNO3 and NaNO3 salt solutions. We attribute this to the different nature of the hydration layers about the alkali metal ions: these favor liganding to the negatively charged mica surface of the etheric -O- group on the ethylene oxide monomer for the case of the more weakly hydrated K+ and Cs+, but not for the case of Na+ or Li+ with their more strongly bound water. A simple model relating the electrostatic energy changes occurring upon such liganding to the experimentally measured hydration energies of the different alkali metal ions supports this attribution.     

Microhydration of cytosine and its radical anion: cytosine.(H2O)n (n=1-5)

Microhydration effects on cytosine and its radical anion have been investigated theoretically, by explicitly considering various structures of cytosine complexes with up to five water molecules. Each successive water molecule (through n=5) is bound by 7-10 kcal mol(-1) to the relevant cytosine complex. The hydration energies are uniformly higher for the analogous anion systems. While the predicted vertical detachment energy (VDE) of the isolated cytosine is only 0.48 eV, it is predicted to increase to 1.27 eV for the lowest-lying pentahydrate of cytosine. The adiabatic electron affinity (AEA) of cytosine was also found to increase from 0.03 to 0.61 eV for the pentahydrate, implying that the cytosine anion, while questionable in the gas phase, is bound in aqueous solution. Both the VDE and AEA values for cytosine are smaller than those of uracil and thymine for a given hydration number. These results are in qualitative agreement with available experimental results from photodetachment-photoelectron spectroscopy studies of Schiedt et al. [Chem. Phys. 239, 511 (1998)].     

Surface relaxation in water clusters: evidence from theoretical analysis of the oxygen 1s photoelectron spectrum

We present a theoretical interpretation of the oxygen 1s photoelectron spectrum published by Ohrwall et al. [J. Chem. Phys. 123, 054310 (2005)]. A water cluster that contains 200 molecules was simulated at 215 K using the polarizable AMOEBA force field. The force field predicts longer O...O distances at the cluster surface than in the bulk. Comparisons to ab initio molecular dynamics (MD) simulations indicate that the force field underestimates the degree of surface relaxation. By comparing cluster lineshape models, computed from MD simulations, to the experimental spectrum we find further evidence of surface relaxation.     

Relaxation and jump dynamics of water at the mica interface

The orientational relaxation dynamics of water confined between mica surfaces is investigated using molecular dynamics simulations. The study illustrates the wide heterogeneity that exists in the dynamics of water adjacent to a strongly hydrophilic surface such as mica. Analysis of the survival probabilities in different layers is carried out by normalizing the corresponding relaxation times with bulk water layers of similar thickness. A 10-fold increase in the survival times is observed for water directly in contact with the mica surface and a non-monotonic variation in the survival times is observed moving away from the mica surface to the bulk-like interior. The orientational relaxation time is highest for water in the contact layer, decreasing monotonically away from the surface. In all cases the ratio of the relaxation times of the 1st and 2nd rank Legendre polynomials of the HH bond vector is found to lie between 1.5 and 1.9 indicating that the reorientational relaxation in the different water layers is governed by jump dynamics. The orientational dynamics of water in the contact layer is particularly novel and is found to undergo distinct two-dimensional hydrogen bond jump reorientational dynamics with an average waiting time of 4.97 ps. The waiting time distribution is found to possess a long tail extending beyond 15 ps. Unlike previously observed jump dynamics in bulk water and other surfaces, jump events in the mica contact layer occur between hydrogen bonds formed by the water molecule and acceptor oxygens on the mica surface. Despite slowing down of the water orientational relaxation near the surface, life-times of water in the hydration shell of the K(+) ion are comparable to that observed in bulk salt solutions.     

Time-dependent conformational changes in fibrinogen measured by atomic force microscopy

Tapping-mode atomic force microscopy was used to study the time-dependent changes in the structure of fibrinogen under aqueous conditions following adsorption on two model surfaces: hydrophobic graphite and hydrophilic mica. Fibrinogen was observed in the characteristic trinodular form, and the dimensions of the adsorbed molecules were consistent with previously reported values for these surfaces. On the basis of the differences in the relative heights of the D and the E domains, four orientation states were observed for fibrinogen adsorbed on both the surfaces. On graphite, the initial asymmetric orientation states disappeared with spreading over time. Some small lateral movements of the adsorbed proteins were observed on mica during repeated scanning, whereas no such movement was observed on graphite, indicating strong adhesion of fibrinogen to a hydrophobic surface. Spreading kinetics of fibrinogen on the two surfaces was determined by measuring the heights of the D and E domains over a time period of approximately 2 h. On graphite, the heights of both the D and E domains decreased with time to a lower plateau value of 1.0 nm. On mica, the heights of both the D and E domains showed an increase, rising to an upper plateau value of approximately 2.1 nm. The spreading of the D and E domains on graphite was analyzed using an 'exponential-decay-of-height' model. A spreading rate constant of approximately 4.7 x 10(-4) s(-1) was observed for the whole fibrinogen molecule adsorbed on graphite, corresponding to a free energy of unfolding of approximately 37 kT. Extrapolation of the exponential curve in the model to t = 0 yielded values of 2.3 and 2.2 nm for the heights of the D and the E domains at the time of contact with the hydrophobic graphite substrate, significantly less than their free solution diameters. A two-step spreading model is proposed to explain this observation.     

Molecular nitrogen-N2 properties: the intermolecular potential and the equation of state

Quantum mechanical (QM) high precision calculations were used to determine N(2)-N(2) intermolecular interaction potential. Using QM numerical data the anisotropic potential energy surface was obtained for all orientations of the pair of the nitrogen molecules in the rotation invariant form. The new N(2)-N(2) potential is in reasonably good agreement with the scaled potential obtained by van der Avoird et al. using the results of Hartree-Fock calculations [J. Chem. Phys. 84, 1629 (1986)]. The molecular dynamics (MD) of the N(2) molecules has been used to determine nitrogen equation of state. The classical motion of N(2) molecules was integrated in rigid rotor approximation, i.e., it accounted only translational and rotational degrees of freedom. Fincham [Mol. Simul. 11, 79 (1993)] algorithm was shown to be superior in terms of precision and energy stability to other algorithms, including Singer [Mol. Phys. 33, 1757 (1977)], fifth order predictor-corrector, or Runge-Kutta, and was therefore used in the MD modeling of the nitrogen pressure [S. Krukowski and P. Strak, J. Chem. Phys. 124, 134501 (2006)]. Nitrogen equation of state at pressures up to 30 GPa (300 kbars) and temperatures from the room temperature to 2000 K was obtained using MD simulation results. Results of MD simulations are in very good agreement (the error below 1%) with the experimental data on nitrogen equation of state at pressures below 1 GPa (10 kbars) for temperatures below 1800 K [R. T. Jacobsen et al., J. Phys. Chem. Ref. Data 15, 735 (1986)]. For higher temperatures, the deviation is slightly larger, about 2.5% which still is a very good agreement. The slightly larger difference may be attributed to the vibrational motion not accounted explicitly by rigid rotor approximation, which may be especially important at high temperatures. These results allow to obtain reliable equation of state of nitrogen for pressures up to 30 GPa (300 kbars), i.e., close to molecular nitrogen stability limit, determined by Nellis et al. [Phys. Rev. Lett. 53, 1661 (1984)].     

Hydration force between mica surfaces in aqueous KCl electrolyte solution

Liquid-vapor molecular dynamics simulations are performed to study the interaction forces between two mica surfaces in an aqueous KCl electrolyte solution. Strong repulsive hydration force is obtained within a distance of ~2 nm between the two mica surfaces, which cannot be explained by the continuum theory of double-layer repulsion. We find that this short-range repulsive hydration force is much stronger than the double-layer force between mica surfaces. Whereas the simulation system is much smaller than the surface force measurement system, fundamental mechanisms of repulsive hydration force are revealed. In particular, important features of the step-like force oscillatory behavior during normal compression and force hysteresis during retraction are observed. Detailed analysis of the ionic density distributions shows that the "forced adsorption" of diffusive K(+) ions onto mica surfaces during compression and the subsequent "slow desorption" of the absorbed K(+) ions from mica surfaces upon retraction are responsible for the hysteresis phenomenon. From a mechanics point of view, we attribute the load bearing capacity of the dense electrolyte to the very hard hydration shells of K(+) metal ions under confinement. We find that the hydrated K(+) ions and Cl(-) co-ions remain very diffusive in the aqueous film. Water molecules in the hydration layer are also very fluidic, in the sense that the diffusion constant of water molecules is less than its bulk value by at most 3 orders of magnitude under the extreme confinement.     

The influence of ion binding and ion specific potentials on the double layer pressure between charged bilayers at low salt concentrations

Measurements of surface forces between double-chained cationic bilayers adsorbed onto molecularly smooth mica surfaces across different millimolar salt solutions have revealed a large degree of ion specificity [Pashley et al., J. Phys. Chem. 90, 1637 (1986)]. This has been interpreted in terms of highly specific anion binding to the adsorbed bilayers. We show here that inclusion in the double layer theory of nonspecific ion binding and ion specific nonelectrostatic potentials acting between ions and the two surfaces can account for the phenomenon. It also gives the right Hofmeister series for the double layer pressure.     

Na(110)和Na2O(111)表面弛豫的研究

We have studied the relaxation of Na(110) and Na_2O(111) surfaces with LEED, and found that the best agreement between theory and experiment in such a structure, the surface layer distance d=0.291±0.01 nm (contracted 0.34%±0.01 nm) [for Na(110)]; and the three atomic surface layer distances d_1=0.117±0.01 nm (expanded 4.63%±0.01 nm), d_2=0.161±0.01 nm(expanded 0.063%±0.01 nm) and d_3=0.086±0.01 nm (expanded 0.75%±0.01 nm) [for Na_2O(111)].     

Molecular dynamics in the isothermal-isobaric ensemble: the requirement of a "shell" molecule. I. Theory and phase-space analysis

Current constant pressure molecular-dynamics (MD) algorithms are not consistent with the recent reformulation of the isothermal-isobaric (NpT) ensemble. The NpT ensemble partition function requires the use of a "shell" molecule to identify uniquely the volume of the system, thereby avoiding the redundant counting of configurations [e.g., G. J. M. Koper and H. Reiss, J. Phys. Chem. 100, 422 (1996); D. S. Corti, Phys. Rev. E, 64, 016128 (2001)]. So far, only the NpT Monte Carlo method has been updated to allow the system volume to be defined by a shell particle [D. S. Corti, Mol. Phys. 100, 1887 (2002)]. A shell particle has yet to be incorporated into MD simulations. The proper modification of the NpT MD algorithm is therefore the subject of this paper. Unlike Andersen's method [H. C. Andersen, J. Chem. Phys. 72, 2384 (1980)] where a piston of unknown mass serves to control the response time of volume fluctuations, the newly proposed equations of motion impose a constant external pressure via the introduction of a shell particle of known mass. Hence, the system itself sets the time scales for pressure and volume fluctuations. The new algorithm is subject to a number of fundamentally rigorous tests to ensure that the equations of motion sample phase space correctly. We also show that the Hoover NpT algorithm [W. G. Hoover, Phys. Rev. A. 31, 1695 (1985); 34, 2499 (1986)] does sample phase correctly, but only when periodic boundary conditions are employed.     

Critical fluctuations of the micellar triethylene glycol monoheptyl ether-water system

Using the equal volume criterion and also the pseudospinodal conception the critical demixing point of the triethylene glycol monoheptyl ether/water system (C7E3H2O) has been determined as Ycrit=0.1 and Tcrit=296.46 K (Y, mass fraction of surfactant). From density measurements the critical micelle concentration (cmc) followed as Ycmc=0.007 at 288.15 K and Ycmc=0.0066 at 298.15 K. The (static) shear viscosity etas and the mutual diffusion coefficient D of the C7E3H2O mixture of critical composition have been evaluated to yield their singular and background parts. From a combined treatment of both quantities the relaxation rate Gamma of order parameter fluctuations has been derived. Gamma follows power law with universal critical exponent and amplitude Gamma0=3.1 x 10(9) s(-1). Broadband ultrasonic spectra of C7E3H2O mixtures exhibit a noncritical relaxation, reflecting the monomer exchange between micelles and the suspending phase, and a critical term due to concentration fluctuations. The former is subject to a relaxation time distribution that broadens when approaching the critical temperature. The latter can be well represented with the aid of the dynamic scaling model by Bhattacharjee and Ferrell (BF) [Phys. Rev. A. 31, 1788 (1985)]. The half-attenuation frequency in the scaling function of the latter model is noticeably smaller (Omega12 (BF) approximately 1) than the theoretically predicted value Omega12 (BF)=2.1. This result has been taken as an indication of a coupling between the fluctuations in the local concentration and the kinetics of micelle formation, in correspondence with the idea of a fluctuation controlled monomer exchange [T. Telgmann and U. Kaatze, Langmuir 18, 3068 (2002)].     

Sound attenuation, shear viscosity, and mutual diffusivity behavior in the nitroethane-cyclohexane critical mixture

The shear viscosity eta(s), mutual diffusion coefficient D, and ultrasonic attenuation spectra of the nitroethane-cyclohexane mixture of critical composition have been measured at various temperatures near the critical temperature T(c). The relaxation rate of order parameter fluctuations resulting from a combined evaluation of the eta(s) and D data follows power law behavior with the theoretical exponent and with the large amplitude Gamma(o)=(156+/-2)x10(9) s(-1). The ultrasonic spectra have been evaluated in terms of a critical contribution and a noncritical background contribution. The amplitude of the former exhibits a temperature dependence, in conformity with a temperature dependence in the adiabatic coupling constant (|g| = 0.064 near T(c) and 0.1 at T-T(c)=3 K). If the variation of the critical amplitude with T is taken into account the experimental attenuation coefficient data display a scaling function which nicely fits to the theoretical prediction from the Bhattacharjee-Ferrell dynamic scaling model [R. A. Ferrell and J. K. Bhattacharjee, Phys. Rev. A 31, 1788 (1985)].