The effect of confinement on water structure

Structure-Thermodynamic Relationship of a Polysaccharide Gel (Alginate) as a Function of Water Content and Counterion Type (Na vs Ca)

Biofilms are the predominant mode of microbial life on Earth, and so a deep understanding of microbial communities─and their impacts on environmental processes─requires a firm understanding of biofilm properties. Because of the importance of biofilms to their microbial inhabitants, microbes have evolved different ways of engineering and reconfiguring the matrix of extracellular polymeric substances (EPS) that constitute the main non-living component of biofilms. This ability makes it difficult to distinguish between the biotic and abiotic origins of biofilm properties. An important route toward establishing this distinction has been the study of simplified models of the EPS matrix. This study builds on such efforts by using atomistic simulations to predict the nanoscale (≤10 nm scale) structure of a model EPS matrix and the sensitivity of this structure to interpolymer interactions and water content. To accomplish this, we use replica exchange molecular dynamics (REMD) simulations to generate all-atom configurations of ten 3.4 kDa alginate polymers at a range of water contents and Ca-Na ratios. Simulated systems are solvated with explicitly modeled water molecules, which allows us to capture the discrete structure of the hydrating water and to examine the thermodynamic stability of water in the gels as they are progressively dehydrated. Our primary findings are that (i) the structure of the hydrogels is highly sensitive to the identity of the charge-compensating cations, (ii) the thermodynamics of water within the gels (specific enthalpy and free energy) are, surprisingly, only weakly sensitive to cation identity, and (iii) predictions of the differential enthalpy and free energy of hydration include a short-ranged enthalpic term that promotes hydration and a longer-ranged (presumably entropic) term that promotes dehydration, where short and long ranges refer to distances shorter or longer than ∼0.6 nm between alginate strands.

Rotational dynamics of proteins in nanochannels: role of solvent's local viscosity

Viscosity variation of solvent in local regions near a solid surface, be it a biological surface of a protein or an engineered surface of a nanoconfinement, is a direct consequence of intermolecular interactions between the solid body and the solvent. The current coarse-grained molecular dynamics study takes advantage of this phenomenon to investigate the anomaly in a solvated protein's rotational dynamics confined using a representative solid matrix. The concept of persistence time, the characteristic time of structural reordering in liquids, is used to compute the solvent's local viscosity. With an increase in the degree of confinement, the confining matrix significantly influences the solvent molecule's local viscosity present in the protein hydration layer through intermolecular interactions. This effect contributes to the enhanced drag force on protein motion, causing a reduction in the rotational diffusion coefficient. Simulation results suggest that the direct matrix-protein non-bonded interaction is responsible for the occasional jump and discontinuity in orientational motion when the protein is in very tight confinement.

Imaging of nuclear magnetic resonance spin-lattice relaxation activation energy in cartilage

Samples of human and bovine cartilage have been examined using magnetic resonance imaging to determine the proton nuclear magnetic resonance spin-lattice relaxation time, T₁, as a function of depth within through the cartilage tissue. T₁ was measured at five to seven temperatures between 8 and 38°C. From this, it is shown that the T₁ relaxation time is well described by Arrhenius-type behaviour and the activation energy of the relaxation process is quantified. The activation energy within the cartilage is approximately 11 ± 2 kJ mol^-1 with this notably being less than that for both pure water (16.6 ± 0.4 kJ mol^-1) and the phosphate-buffered solution in which the cartilage was immersed (14.7 ± 1.0 kJ mol^-1). It is shown that this activation energy increases as a function of depth in the cartilage. It is known that cartilage composition varies with depth, and hence, these results have been interpreted in terms of the structure within the cartilage tissue and the association of the water with the macromolecular constituents of the cartilage.

Confinement Effects on the Benzene Orientational Structure

Liquids under confinement exhibit different properties compared with their corresponding bulk phases, for example, miscibility, phase transitions, and diffusion. The underlying cause is the local ordering of molecules, which is usually only studied using pure simulation methods. Herein, we derive experimentally the structure of benzene confined in MCM-41 using total neutron scattering measurements. The study reveals a layering of molecules across a pore, and four concentric cylindrical shells can be distinguished for a pore with the radius of 18 Å. The nanoscale confinement of the liquid has a major effect on the spatial and orientational correlations observed between the molecules, when compared with the structure of the bulk liquid. These differences are most marked for molecules in parallel configurations, and this suggests differences in chemical reactivity between the confined and bulk liquids.

Counterion-Regulated Dynamics of Water Confined in Lyotropic Liquid Crystalline Morphologies

The dynamics of confined water is of fundamental and long-standing interest. In technologically important forms of confinement, such as proton-exchange membranes, electrostatic interactions with the confining matrix and counterions play significant roles on the properties of water. There has been recent interest on the dynamics of water confined to the lyotropic liquid crystalline (LLC) morphologies of Gemini dicarboxylate surfactants. These systems are exciting because the nature of confinement, for example, size and curvature of channels and surface functionality is dictated by the chemistry of the self-assembling surfactant molecules. Quasielastic neutron scattering experiments have shown an interesting dependence of the water self-diffusion constant, D_α, on the identity (denoted α) of the counterion: at high hydration, the magnitude of the water self-diffusion constant is in the order D_TMA < D_Na < D_K, where TMA, Na, and K refer to tetramethyl ammonium, sodium, and potassium counterions, respectively. This sequence is similar to what is seen in bulk electrolyte solutions. At low hydrations, however, the order of water self-diffusion is different, that is, D_Na < D_TMA < D_K. In this work, we present molecular dynamics simulations for the dynamics of water in the LLC phases of dicarboxylate Gemini surfactants. The simulations reproduce the trends seen in experiments. From an analysis of the trajectories, we hypothesize that two competing factors play a role: the volume accessible to the water molecules and the correlations between the water and the counterion. The excluded volume effect is the largest with TMA⁺, and the electrostatic correlation is the strongest with Na⁺. The observed trend is a result of which of these two effects is dominant at a given water to surfactant ratio.

Time-Resolved Light Scattering Study on the Kinetics of the Liquid-Liquid Transition in Triphenyl Phosphite

There is experimental evidence suggesting the existence of a liquid-liquid transition (LLT) in a single-component liquid. However, none of this evidence is free from controversy, including the case of a molecular liquid, triphenyl phosphite, which we study here. Furthermore, the kinetics of LLT has been largely unexplored. Here we study the phase-transition dynamics of triphenyl phosphite in a supercooled liquid state by means of time-resolved polarized and depolarized small-angle light scattering to clarify whether the transition is a liquid-liquid transition (LLT) or merely nanocrystal formation. A part of this study was recently reported in another of our papers [Shimizu, R.; Kobayashi, M.; Tanaka, H. Phys. Rev. Lett. 2014, 112, 125702]. A detailed analysis of our experimental results of light scattering and the comparison with heat evolution during LLT have revealed the following facts. The polarized scattering from domains has a finite (nonzero) intensity in the low-wavenumber limit, and the time evolution of its average intensity is almost proportional to the square of the heat-releasing rate. The depolarized scattering intensity monotonically increases in the process of LLT during isothermal annealing above the spinodal temperature TSD but exhibits a peak below TSD. On the basis of these results, we suggest that the primary process is LLT, whose order parameter is of a nonconserved nature, but accompanies nanocrystal formation. In the NG-type LLT, the sharp interface between liquid II droplets and the liquid I matrix promotes nanocrystal formation there, whereas much less nanocrystal formation is induced in the SD-type LLT due to the lack of such sharp interfaces.

Microscopic identification of the order parameter governing liquid-liquid transition in a molecular liquid

A liquid-liquid transition (LLT) in a single-component substance is an unconventional phase transition from one liquid to another. LLT has recently attracted considerable attention because of its fundamental importance in our understanding of the liquid state. To access the order parameter governing LLT from a microscopic viewpoint, here we follow the structural evolution during the LLT of an organic molecular liquid, triphenyl phosphite (TPP), by time-resolved small- and wide-angle X-ray scattering measurements. We find that locally favored clusters, whose characteristic size is a few nanometers, are spontaneously formed and their number density monotonically increases during LLT. This strongly suggests that the order parameter of LLT is the number density of locally favored structures and of nonconserved nature. We also show that the locally favored structures are distinct from the crystal structure and these two types of orderings compete with each other. Thus, our study not only experimentally identifies the structural order parameter governing LLT, but also may settle a long-standing debate on the nature of the transition in TPP, i.e., whether the transition is LLT or merely microcrystal formation.

Radical re-appraisal of water structure in hydrophilic confinement

The structure of water confined in MCM41 silica cylindrical pores is studied to determine whether confined water is simply a version of the bulk liquid which can be substantially supercooled without crystallisation. A combination of total neutron scattering from the porous silica, both wet and dry, and computer simulation using a realistic model of the scattering substrate is used. The water in the pore is divided into three regions: core, interfacial and overlap. The average local densities of water in these simulations are found to be about 20% lower than bulk water density, while the density in the core region is below, but closer to, the bulk density. There is a decrease in both local and core densities when the temperature is lowered from 298 K to 210 K. The radical proposal is made here that water in hydrophilic confinement is under significant tension, around -100 MPa, inside the pore.

Bond orientational order in liquids: Towards a unified description of water-like anomalies, liquid-liquid transition, glass transition, and crystallization: Bond orientational order in liquids

There are at least three fundamental states of matter, depending upon temperature and pressure: gas, liquid, and solid (crystal). These states are separated by first-order phase transitions between them. In both gas and liquid phases a complete translational and rotational symmetry exist, whereas in a solid phase both symmetries are broken. In intermediate phases between liquid and solid, which include liquid crystal and plastic crystal phases, only one of the two symmetries is preserved. Among the fundamental states of matter, the liquid state is the most poorly understood. We argue that it is crucial for a better understanding of liquids to recognize that a liquid generally has the tendency to have a local structural order and its presence is intrinsic and universal to any liquid. Such structural ordering is a consequence of many-body correlations, more specifically, bond angle correlations, which we believe are crucial for the description of the liquid state. We show that this physical picture may naturally explain difficult unsolved problems associated with the liquid state, such as anomalies of water-type liquids (water, Si, Ge, ...), liquid-liquid transition, liquid-glass transition, crystallization and quasicrystal formation, in a unified manner. In other words, we need a new order parameter representing a low local free-energy configuration, which is a bond orientational order parameter in many cases, in addition to a density order parameter for the physical description of these phenomena. Here we review our two-order-parameter model of liquid and consider how transient local structural ordering is linked to all of the above-mentioned phenomena. The relationship between these phenomena is also discussed.

Liquid-liquid transition without macroscopic phase separation in a water-glycerol mixture

The existence of more than two liquid states in a single-component substance and the ensuing liquid-liquid transitions (LLTs) has attracted considerable attention because of its counterintuitive nature and its importance in the fundamental understanding of the liquid state. Here we report direct experimental evidence for a genuine (isocompositional) LLT without macroscopic phase separation in an aqueous solution of glycerol. We show that liquid I transforms into liquid II by way of two types of kinetics: nucleation and growth, and spinodal decomposition. Although liquid II is metastable against crystallization, we could access both its static and dynamical properties experimentally. We find that liquids I and II differ in density, refractive index, structure, hydrogen bonding state, glass transition temperature and fragility, and that the transition between the two liquids is mainly driven by the local structuring of water rather than of glycerol, suggesting a link to a plausible LLT in pure water.

Density profile of water confined in cylindrical pores in MCM-41 silica

Recently, water absorbed in the porous silica material MCM-41-S15 has been used to demonstrate an apparent fragile to strong dynamical crossover on cooling below ∼220 K, and also to claim that the density of confined water reaches a minimum at a temperature around 200 K. Both of these behaviours are purported to arise from the crossing of a Widom line above a conjectured liquid-liquid critical point in bulk water. Here it is shown that traditional estimates of the pore diameter in this porous silica material (of order 15 Å) are too small to allow the amount of water that is observed to be absorbed by these materials (around 0.5 g H(2)O/g substrate) to be absorbed only inside the pore. Either the additional water is absorbed on the surface of the silica particles and outside the pores, or else the pores are larger than the traditional estimates. In addition the low Q Bragg intensities from a sample of MCM-41-S15 porous silica under different dry and wet conditions and with different hydrogen isotopes are simulated using a simple model of the water and silica density profile across the pore. It is found the best agreement of these intensities with experimental data is shown by assuming the much larger pore diameter of 25 Å (radius 12.5 Å). Qualitative agreement is found between these simulated density profiles and those found in recent empirical potential structure refinement simulations of the same data, even though the latter data did not specifically include the Bragg peaks in the structure refinement. It is shown that the change in the (100) peak intensity on cooling from 300 to 210 K, which previously has been ascribed to a change in density of the confined water on cooling, can equally be ascribed to a change in density profile at constant average density. It is further pointed out that, independent of whether the pore diameter really is as large as 25 Å or whether a significant amount of water is absorbed outside the pore, the earlier reports of a dynamic crossover in supercooled confined water could in fact be a crystallization transition in the larger pore or surface water.

Density minimum of confined water at low temperatures: a combined study by small-angle scattering of X-rays and neutrons

A simple explanation is given for the low-temperature density minimum of water confined within cylindrical pores of ordered nanoporous materials of different pore size. The experimental evidence is based on combined data from in-situ small-angle scattering of X-rays (SAXS) and neutrons (SANS), corroborated by additional wide-angle X-ray scattering (WAXS). The combined scattering data cannot be described by a homogeneous density distribution of water within the pores, as was originally suggested from SANS data alone. A two-step density model reveals a wall layer covering approximately two layers of water molecules with higher density than the residual core water in the central part of the pores. The temperature-induced changes of the scattering signal from both X-rays and neutrons are consistent with a minimum of the average water density. We show that the temperature at which this minimum occurs depends monotonically on the pore size. Therefore we attribute this minimum to a liquid-solid transition of water influenced by confinement. For water confined in the smallest pores of only 2 nm in diameter, the density minimum is explained in terms of a structural transition of the surface water layer closest to the hydrophilic pore walls.