Structure and dynamics of water confined in silica nanopores

We report the results of molecular simulation of water in silica nanopores at full hydration and room temperature. The model systems are approximately cylindrical pores in amorphous silica, with diameters ranging from 20 to 40 ?. The filled pores are prepared using grand canonical Monte Carlo simulation and molecular dynamics simulation is used to calculate the water structure and dynamics. We found that water forms two distinct molecular layers at the interface and exhibits uniform, but somewhat lower than bulk liquid, density in the core region. The hydrogen bond density profile follows similar trends, with lower than bulk density in the core and enhancements at the interface, due to hydrogen bonds between water and surface non-bridging oxygens and OH groups. Our studies of water dynamics included translational mean squared displacements, orientational time correlations, survival probabilities in interfacial shells, and hydrogen bond population relaxation. We found that the radial-axial anisotropy in translational motion largely follows the predictions of a model of free diffusion in a cylinder. However, both translational and rotational water mobilities are strongly dependent on the proximity to the interface, with pronounced slowdown in layers near the interface. Within these layers, the effects of interface curvature are relatively modest, with only a small increase in mobility in going from the 20 to 40 ? diameter pore. Hydrogen bond population relaxation is nearly bulk-like in the core, but considerably slower in the interfacial region.     

Structure and dynamics of liquid methanol confined within functionalized silica nanopores

Molecular dynamics simulations have been carried out to investigate the structure and dynamics of liquid methanol confined in 3.3 nm diameter cylindrical silica pores. Three cavities differing in the characteristics of the functional groups at their walls have been examined: (i) smooth hydrophobic pores in which dispersive forces prevail, (ii) hydrophilic cavities with surfaces covered by polar silanol groups, and (iii) a much more rugged pore in which 60% of the previous interfacial hydroxyl groups were replaced by the bulkier trimethylsilyl ones. Confinement promotes a considerable structure at the vicinity of the pore walls which is enhanced in the case of hydroxylated surfaces. Moreover, in the presence of the trimethylsilyl groups, the propagation of this interface-induced spatial ordering extends down to the central region of the pore. Concerning the dynamical modes, we observed an overall slowdown in both the translational and rotational motions. An analysis of these mobilities from a local perspective shows that the largest retardations operate at the vicinity of the interfaces. The gross features of the rotational dynamics were analyzed in terms of contributions arising from bulk and surface states. Compared to the bulk dynamical behavior, the characteristic timescales associated with the rotational motions show the most dramatic increments. A dynamical analysis of hydrogen bond formation and breaking processes is also included.     

Solvated calcium ions in charged silica nanopores

Hydroxyl surface density in porous silica drops down to nearly zero when the pH of the confined aqueous solution is greater than 10.5. To study such extreme conditions, we developed a model of slit silica nanopores where all the hydrogen atoms of the hydroxylated surface are removed and the negative charge of the resulting oxygen dangling bonds is compensated by Ca(2+) counterions. We employed grand canonical Monte Carlo and molecular dynamics simulations to address how the Ca(2+) counterions affect the thermodynamics, structure, and dynamics of confined water. While most of the Ca(2+) counterions arrange themselves according to the so-called "Stern layer," no diffuse layer is observed. The presence of Ca(2+) counterions affects the pore filling for strong confinement where the surface effects are large. At full loading, no significant changes are observed in the layering of the first two adsorbed water layers compared to nanopores with fully hydroxylated surfaces. However, the water structure and water orientational ordering with respect to the surface is much more disturbed. Due to the super hydrophilicity of the Ca(2+)-silica nanopores, water dynamics is slowed down and vicinal water molecules stick to the pore surface over longer times than in the case of hydroxylated silica surfaces. These findings, which suggest the breakdown of the linear Poisson-Boltzmann theory, provide important information about the properties of nanoconfined electrolytes upon extreme conditions where the surface charge and ion concentration are large.     

Water adsorption in ion-bearing nanopores

Grand canonical Monte Carlo simulations are used to examine the adsorption of water into cylindrical nanopores containing single ions. The isotherms for water adsorbing into nanopores with radii of 0.44, 0.54, 0.64, and 0.74 nm and containing Na+, K+, Ca2+, Cl-, or F- at 298 K are computed. In all cases the nanopores are found to fill at reservoir chemical potentials below the chemical potential of saturated water vapor at 298 K. The threshold chemical potential is found to be sensitive to both the size of the channel and the ion species, with the anion-bearing pores filling at lower chemical potentials. Additionally, the filling threshold chemical potential is found to decrease as the radius of the pores is decreased. Pores with K+ and Cl- are compared, and the Cl- pores are found to exhibit higher water densities in the filled states and a more energetically favorable water structure while yielding lower per particle entropies. Sample simulation configurations are also examined and indicate that at low chemical potentials, the adsorbed water forms a cluster around the ion. Finally, the influence of the choice of water model on the adsorption isotherms is examined.     

Size-selective molecular transport through silica colloidal nanopores

Diffusion rate of dye-labelled PAMAM dendrimers through free-standing silica colloidal crystals was studied as a function of the dendrimer generation and nanopore size to determine the transport selectivity.     

Aqueous electrolytes confined within functionalized silica nanopores

Molecular dynamics simulations have been carried out to investigate structural and dynamical characteristics of NaCl aqueous solutions confined within silica nanopores in contact with a "bulk-like" reservoir. Two types of pores, with diameters intermediate between 20 A? and 37.5 A?, were investigated: The first one corresponded to hydrophobic cavities, in which the prevailing wall-solution interactions were of the Lennard-Jones type. In addition, we also examined the behavior of solutions trapped within hydrophilic cavities, in which a set of unsaturated O-sites at the wall were transformed in polar silanol Si-OH groups. In all cases, the overall concentrations of the trapped electrolytes exhibited important reductions that, in the case of the narrowest pores, attained 50% of the bulk value. Local concentrations within the pores also showed important fluctuations. In hydrophobic cavities, the close vicinity of the pore wall was coated exclusively by the solvent, whereas in hydrophilic pores, selective adsorption of Na(+) ions was also observed. Mass and charge transport were also investigated. Individual diffusion coefficients did not present large modifications from what is perceived in the bulk; contrasting, the electrical conductivity exhibited important reductions. The qualitative differences are rationalized in terms of simple geometrical considerations.     

Melting and freezing of water in cylindrical silica nanopores

Freezing and melting of H(2)O and D(2)O in the cylindrical pores of well-characterized MCM-41 silica materials (pore diameters from 2.5 to 4.4 nm) was studied by differential scanning calorimetry (DSC) and (1)H NMR cryoporometry. Well-resolved DSC melting and freezing peaks were obtained for pore diameters down to 3.0 nm, but not in 2.5 nm pores. The pore size dependence of the melting point depression DeltaT(m) can be represented by the Gibbs-Thomson equation when the existence of a layer of nonfreezing water at the pore walls is taken into account. The DSC measurements also show that the hysteresis connected with the phase transition, and the melting enthalpy of water in the pores, both vanish near a pore diameter D* approximately equal to 2.8 nm. It is concluded that D* represents a lower limit for first-order melting/freezing in the pores. The NMR spin echo measurements show that a transition from low to high mobility of water molecules takes place in all MCM-41 materials, including the one with 2.5 nm pores, but the transition revealed by NMR occurs at a higher temperature than indicated by the DSC melting peaks. The disagreement between the NMR and DSC transition temperatures becomes more pronounced as the pore size decreases. This is attributed to the fact that with decreasing pore size an increasing fraction of the water molecules is situated in the first and second molecular layers next to the pore wall, and these molecules have slower dynamics than the molecules in the core of the pore.     

Proton quantum coherence observed in water confined in silica nanopores

Deep inelastic neutron scattering measurements of water confined in nanoporous xerogel powders, with average pore diameters of 24 and 82 A, have been carried out for pore fillings ranging from 76% to nearly full coverage. DINS measurements provide direct information on the momentum distribution n(p) of protons, probing the local structure of the molecular system. The observed scattering is interpreted within the framework of the impulse approximation and the longitudinal momentum distribution determined using a model independent approach. The results show that the proton momentum distribution is highly non-Gaussian. A bimodal distribution appears in the 24 A pore, indicating coherent motion of the proton over distances d of approximately 0.3 A. The proton mean kinetic energy W of the confined water molecule is determined from the second moment of n(p). The W values, higher than in bulk water, are ascribed to changes of the proton dynamics induced by the interaction between interfacial water and the confining surface.     

Water dynamics in graphite oxide investigated with neutron scattering

Graphite oxide is an inorganic multilayer system that preserves the layered structure of graphite but not the conjugated bond structure. In the past few years, detailed studies of the static structure of graphite oxide were carried out. This was mainly done by NMR investigations and led to a new structural model of graphite oxide. The layer distance of graphite oxide increases with increasing humidity level, giving rise to different spacings of the carbon layers in the range from 6 to 12 A. As a consequence, different types of motions of water and functional groups appear. Information about the mobility of the water molecules is not yet complete but is crucial for the understanding of the structure of the carbon layers as well as the intercalation process. In this paper, the hydration- and temperature-dependent dynamic behavior of graphite oxide will be investigated by quasielastic neutron scattering using the time-of-flight spectrometer NEAT at the Hahn-Meitner-Institut Berlin. The character of the embedded water does not change over a wide range of hydration levels. Especially the interlayer water remains tightly bound and does not show any translational motion. In samples with excess water, however, the water is also distributed in noninterlayer voids, leading to the observation of additional motions of bulklike or confined water. The dynamic behavior of hydrated graphite oxide can be described by a consistent model that combines two two-site jump motions for the motions of the water molecules and the motions of OH groups.     

Structural characterization of surfactant aggregates adsorbed in cylindrical silica nanopores

The self-assembly of nonionic surfactants in the cylindrical pores of SBA-15 silica with a pore diameter of 8 nm was studied by small-angle neutron scattering (SANS) at different solvent contrasts. The alkyl ethoxylate surfactants C(10)E(5) and C(12)E(5) exhibit strong aggregative adsorption in the pores as indicated by the sigmoidal shape of the adsorption isotherms. The SANS intensity profiles can be represented by a sum of two terms, one accounting for diffuse scattering from surfactant aggregates in the pores and the other for Bragg scattering from the pore lattice of the silica matrix. The Bragg reflections are analyzed with a form factor model in which the radial density profile of the surfactant in the pore is approximated by a two-step function. Diffuse scattering is represented by a Teubner-Strey-type scattering function which indicates a preferred distance between adsorbed surface aggregates in the pores. Our results suggest that adsorption starts with formation of discrete surface aggregates which increase in number and eventually merge to interconnected patches as the plateau value of the adsorption isotherm is approached. A grossly different behavior, viz. formation of micelles as in solution, is found for the maltoside surfactant C(10)G(2), in agreement with the observed weak adsorption of this surfactant in SBA-15.     

Solidification and melting of cetane confined in the nanopores of silica gel

Synthetic silica gels with six different effective diameters of nano-pores (3–30 nm) were loaded with n-hexadecane (cetane) after the elimination of adsorbed water. Kinetics of the solidification and melting of cetane was studied by differential scanning calorimetry (DSC) above the room temperature. Two thermodynamically different states of cetane were found in the samples: the free (bulk)-cetane state and the confined-cetane state. As suspected, the third state of cetane can be amorphous. This has been indicated by the small total transformation heat. The complex crystallization effect of cetane has been found to obey the nucleation-and-growth kinetics and also to depend on the dimensions of confining pores of silica gel. The melting of cetane seems to vary only with the average diameter of silica gel pores, which satisfies the Gibbs–Thompson relation. The presented results validate the applicability of the DSC technique for the porometry. The cetane-medium calibration curve for the silica gel nano-thermoporometry has been determined.     

Langevin dynamics simulations of polymer translocation through nanopores

We investigate the dynamics of polymer translocation through a nanopore using two-dimensional Langevin dynamics simulations. In the absence of an external driving force, we consider a polymer which is initially placed in the middle of the pore and study the escape time tau(e) required for the polymer to completely exit the pore on either side. The distribution of the escape times is wide and has a long tail. We find that tau(e) scales with the chain length N as tau(e) approximately N(1+2nu), where nu is the Flory exponent. For driven translocation, we concentrate on the influence of the friction coefficient xi, the driving force E, and the length of the chain N on the translocation time tau, which is defined as the time duration between the first monomer entering the pore and the last monomer leaving the pore. For strong driving forces, the distribution of translocation times is symmetric and narrow without a long tail and tau approximately E(-1). The influence of xi depends on the ratio between the driving and frictional forces. For intermediate xi, we find a crossover scaling for tau with N from tau approximately N(2nu) for relatively short chains to tau approximately N(1+nu) for longer chains. However, for higher xi, only tau approximately N(1+nu) is observed even for short chains, and there is no crossover behavior. This result can be explained by the fact that increasing xi increases the Rouse relaxation time of the chain, in which case even relatively short chains have no time to relax during translocation. Our results are in good agreement with previous simulations based on the fluctuating bond lattice model of polymers at intermediate friction values, but reveal additional features of dependency on friction.     

Water dynamics in hardened ordinary Portland cement paste or concrete: from quasielastic neutron scattering

Portland cement reacts with water to form an amorphous paste through a chemical reaction called hydration. In concrete the formation of pastes causes the mix to harden and gain strength to form a rock-like mass. Within this process lies the key to a remarkable peculiarity of concrete: it is plastic and soft when newly mixed, strong and durable when hardened. These qualities explain why one material, concrete, can build skyscrapers, bridges, sidewalks and superhighways, houses, and dams. The character of the concrete is determined by the quality of the paste. Creep and shrinkage of concrete specimens occur during the loss and gain of water from cement paste. To better understand the role of water in mature concrete, a series of quasielastic neutron scattering (QENS) experiments were carried out on cement pastes with water/cement ratio varying between 0.32 and 0.6. The samples were cured for about 28 days in sealed containers so that the initial water content would not change. These experiments were carried out with an actual sample of Portland cement rather than with the components of cement studied by other workers. The QENS spectra differentiated between three different water interactions: water that was chemically bound into the cement paste, the physically bound or "glassy water" that interacted with the surface of the gel pores in the paste, and unbound water molecules that are confined within the larger capillary pores of cement paste. The dynamics of the "glassy" and "unboud" water in an extended time scale, from a hundred picoseconds to a few nanoseconds, could be clearly differentiated from the data. While the observed motions on the picosecond time scale are mainly stochastic reorientations of the water molecules, the dynamics observed on the nanosecond range can be attributed to long-range diffusion. Diffusive motion was characterized by diffusion constants in the range of (0.6-2) 10(-9) m(2)/s, with significant reduction compared to the rate of diffusion for bulk water. This reduction of the water diffusion is discussed in terms of the interaction of the water with the calcium silicate gel and the ions present in the pore water.     

Water structure and dynamics at a silica surface: Pake doublets in 1H NMR spectra

Detailed knowledge about the dynamics and structure of liquids in the vicinity of a solid surface is important in several fields of research. In this study a homogeneous model system of colloidal and nonporous silica particles with a narrow particle size distribution was used to examine such properties of adsorbed water and 1-heptanol. Doublet (1)H water resonances ("Pake doublets") indicate a preferred spatial orientation for the water molecules, as well as a lower molecular density in the surface-induced water structures compared to bulk water. These surface-induced structures are found to extend at least 8 nm from the silica surface. T(1) relaxation measurements at several temperatures indicate weaker H-bonding in the adsorbed water compared to bulk water. T(2) relaxation measurements at several temperatures reveal the presence of two water phases and give quantitative information on the mobility of water molecules and proton exchange processes. The presence of 1-heptanol changes the water characteristics, primarily in the water phase closer to the surface, where water molecules experience decreased translational and increased rotational freedom. In the absence of water, adsorbed 1-heptanol forms surface aggregates encompassing several molecular layers, where the first adsorbed layer shows severe restrictions in mobility and subsequent layers are more mobile.     

Langevin dynamics simulations of ds-DNA translocation through synthetic nanopores

We have implemented a coarse-grained model to study voltage-driven as-DNA translocation through nanopores located in synthetic membranes. The simulated trajectory of the DNA through the nanopores was calculated using Langevin dynamics. We present the results based on more than 120,000 individual translocations. We are particularly interested in this work in probing the physical basis of various experimentally observed--yet poorly understood--phenomena. Notably, we observe in our simulations the formation of ds-DNA hairpins, widely suspected to be the basis for quantized blockage. We study the translocation time, a measurable quantity crucially important in polyelectrolyte characterization, as a function of hairpin vertex location along the polymer backbone, finding that this behavior can be tuned to some degree by simulation parameters. We also study the voltage dependence of the tendency of hairpins to serve as the initiators of translocation events. Surprisingly, we find that the resulting probability depends vitally upon whether the events counted are ultimately successful or not. Further details lead us to propose that failed attempts in experimental translocation studies may be more common--and deceptive--than is generally recognized. We find the time taken by successful single file translocations to be directly proportional to the ratio of chain length to the applied voltage. Finally, we address a common yet puzzling phenomenon in translocation experiments: translocation events in which the current through the pore is highly, yet incompletely, blocked. We present the findings that offer a new explanation for such events.     

Water structure in nanopores of agarose gel by Raman spectroscopy

The evolution of water structure during the gelation process is examined in aqueous solution of agarose using Raman spectroscopy of the O-H stretching band. The measurements have been performed at room temperature for different concentrations of agarose, which yields different dimensions of nanopores in the network of the created gel. Our results show that water confined in the gel pores exhibits evident changes in the local order of molecules in comparison with bulk water and water in the sol state. During the sol-gel transition the number of molecules that participate in the regular tetrahedral H-bond structure increases, and the effect is stronger for higher concentration of the biopolymer.     

Structures of alkali metals in silica gel nanopores: new materials for chemical reductions and hydrogen production

Alkali metals and their alloys can be protected from spontaneous reaction with dry air by intercalation (with subsequent heating) into the pores of silica gel (SG) at loadings up to 40 wt %. The resulting loose, black powders are convenient materials for chemical reduction of organic compounds and the production of clean hydrogen. The problem addressed in this paper is the nature of the reducing species present in these amorphous materials. The atomic pair distribution function (PDF), which considers both Bragg and diffuse scattering components, was used to examine their structures. Liquid Na-K alloys added to silica gel at room temperature (stage 0) or heated to 150 degrees C (stage I) as well as stage I Na-SG, retain the overall pattern of pure silica gel. Broad oscillations in the PDF show that added alkali metals remain in the pores as nanoscale metal clusters. 23Na MAS NMR studies confirm the presence of Na(0) and demonstrate that Na+ ions are formed as well. The relative amounts of Na(0) and Na(+) depend on both the overall metal loading and the average pore size. The results suggest that ionization occurs near or in the SiO2 walls, with neutral metal present in the larger cavities. The fate of the electrons released by ionization is uncertain, but they may add to the silica gel lattice, or form an "electride-like plasma" near the silica gel walls. A remaining mystery is why the stage I material does not show a melting endotherm of the encapsulated metal and does not react with dry oxygen. Na-SG when heated to 400 degrees C (stage II) yields a dual-phase reaction product that consists of Na(4)Si(4) and Na(2)SiO(3).     

Function of membrane protein in silica nanopores: incorporation of photosynthetic light-harvesting protein LH2 into FSM

A high amount of functional membrane protein complex was introduced into a folded-sheet silica mesoporous material (FSM) that has nanometer-size pores of honeycomb-like hexagonal cylindrical structure inside. The photosynthetic light-harvesting complex LH2, which is a typical membrane protein, has a cylindrical structure of 7.3 nm diameter and contains 27 bacteriochlorophyll a and nine carotenoid molecules. The complex captures light energy in the anoxygenic thermophilic purple photosynthetic bacterium Thermochromatium tepidum. The amount of LH2 adsorbed to FSM was determined optically and by the adsorption isotherms of N2. The FSM compounds with internal pore diameters of 7.9 and 2.7 nm adsorbed LH2 at 1.11 and 0.24 mg/mg FSM, respectively, suggesting the high specific affinity of LH2 to the interior of the hydrophobic nanopores with a diameter of 7.9 nm. The LH2 adsorbed to FSM showed almost intact absorption bands of bacteriochlorophylls, and was fully active in the capture and transfer of excitation energy. The LH2 complex inside the FSM showed increased heat stability of the exciton-type absorption band of bacteriochlorophylls (B850), suggesting higher circular symmetry. The environment inside the hydrophobic silica nanopores can be a new matrix for the membrane proteins to reveal their functions. The silica-membrane protein adduct will be useful for the construction of new probes and reaction systems.