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.     

H NMR studies of poly(N-isopropylacrylamide) gels near the phase transition

The phase transition and critical phenomenon of equilibrium swollen poly(N-isopropylacrylamide) (NIPA) hydrogels were studied by ¹H NMR spectroscopy in liquid solution mode. The quantitative NMR observation shows that the peak height and line width of polymer proton and of the HOD proton, and relaxation times of HOD proton all transitionally change as the temperature approaches the transition temperature. The relaxation times of water protons are also measured quantitatively, which shows that the temperature dependence of relaxation times of HOD on temperature before the transition is not consistent with relaxation theory based on the assumption of dominated dipolar interaction between like-spin nuclei and isotropic rotational motion. To explain the surprising relaxation behavior of HOD, we suggest that the amount of bound water in gels increases gradually with temperature at the approach of the phase transition. The pulsed-gradient spin-echo NMR experiments of NIPA gel confirm this suggestion. We believe that these results have important implications concerning the mechanism of the phase transition of NIPA hydrogels.     

Phase diagram and structural evolution of Ag-Au bimetallic nanoparticles: molecular dynamics simulations

We studied the structural evolution of a 270-atom Ag-Au bimetallic nanoparticle (2 nm in size) with varying composition and temperature. The liquid to solid transition region and the solid-state structure were investigated using molecular dynamics simulations. To determine the exact transition temperature region, we applied the mean square displacement and structure deviation methods, as well as the generally used caloric curve of potential energy versus temperature. The results showed that a complete solid-solution phase diagram of the binary Ag-Au system was obtained. Irrespective of the composition, the freezing temperature of a Ag-Au bimetallic nanoparticle was lower than that of the bulk state by a margin of several hundred degrees, and three different solid-state structures are proposed in relation to the Au composition. Our phase diagram offers guidance for the application of Ag-Au nanoparticles.     

An off-lattice Wang-Landau study of the coil-globule and melting transitions of a flexible homopolymer

The Wang-Landau Monte Carlo approach is applied to the coil-globule and melting transitions of off-lattice flexible homopolymers. The solid-liquid melting point and coil-globule transition temperatures are identified by their respective peaks in the heat capacity as a function of temperature. The melting and theta points are well separated, indicating that the coil-globule transition occurs separately from melting even in the thermodynamic limit. We also observe a feature in the heat capacity between the coil-globule and melting transitions which we attribute to a transformation from a low-density liquid globule to a high-density liquid globule.     

Waterlike glass polyamorphism in a monoatomic isotropic Jagla model

We perform discrete-event molecular dynamics simulations of a system of particles interacting with a spherically-symmetric (isotropic) two-scale Jagla pair potential characterized by a hard inner core, a linear repulsion at intermediate separations, and a weak attractive interaction at larger separations. This model system has been extensively studied due to its ability to reproduce many thermodynamic, dynamic, and structural anomalies of liquid water. The model is also interesting because: (i) it is very simple, being composed of isotropically interacting particles, (ii) it exhibits polyamorphism in the liquid phase, and (iii) its slow crystallization kinetics facilitate the study of glassy states. There is interest in the degree to which the known polyamorphism in glassy water may have parallels in liquid water. Motivated by parallels between the properties of the Jagla potential and those of water in the liquid state, we study the metastable phase diagram in the glass state. Specifically, we perform the computational analog of the protocols followed in the experimental studies of glassy water. We find that the Jagla potential calculations reproduce three key experimental features of glassy water: (i) the crystal-to-high-density amorphous solid (HDA) transformation upon isothermal compression, (ii) the low-density amorphous solid (LDA)-to-HDA transformation upon isothermal compression, and (iii) the HDA-to-very-high-density amorphous solid (VHDA) transformation upon isobaric annealing at high pressure. In addition, the HDA-to-LDA transformation upon isobaric heating, observed in water experiments, can only be reproduced in the Jagla model if a free surface is introduced in the simulation box. The HDA configurations obtained in cases (i) and (ii) are structurally indistinguishable, suggesting that both processes result in the same glass. With the present parametrization, the evolution of density with pressure or temperature is remarkably similar to the corresponding experimental measurements on water. Our simulations also suggest that the Jagla potential may reproduce features of the HDA-VHDA transformations observed in glassy water upon compression and decompression. Snapshots of the system during the HDA-VHDA and HDA-LDA transformations reveal a clear segregation between LDA and HDA but not between HDA and VHDA, consistent with the possibility that LDA and HDA are separated by a first order transformation as found experimentally, whereas HDA and VHDA are not. Our results demonstrate that a system of particles with simple isotropic pair interactions, a Jagla potential with two characteristic length scales, can present polyamorphism in the glass state as well as reproducing many of the distinguishing properties of liquid water. While most isotropic pair potential models crystallize readily on simulation time scales at the low temperatures investigated here, the Jagla potential is an exception, and is therefore a promising model system for the study of glass phenomenology.     

Liquid-liquid phase transitions in supercooled water studied by computer simulations of various water models

Liquid-liquid and liquid-vapor coexistence regions of various water models were determined by Monte Carlo (MC) simulations of isotherms of density fluctuation-restricted systems and by Gibbs ensemble MC simulations. All studied water models show multiple liquid-liquid phase transitions in the supercooled region: we observe two transitions of the TIP4P, TIP5P, and SPCE models and three transitions of the ST2 model. The location of these phase transitions with respect to the liquid-vapor coexistence curve and the glass temperature is highly sensitive to the water model and its implementation. We suggest that the apparent thermodynamic singularity of real liquid water in the supercooled region at about 228 K is caused by an approach to the spinodal of the first (lowest density) liquid-liquid phase transition. The well-known density maximum of liquid water at 277 K is related to the second liquid-liquid phase transition, which is located at positive pressures with a critical point close to the maximum. A possible order parameter and the universality class of liquid-liquid phase transitions in one-component fluids are discussed.     

Dynamics near a liquid-liquid phase transition in a non-tetrahedral liquid: the case of gallium

We use molecular simulation to analyze liquid dynamics in the vicinity of the liquid-liquid phase transition (LLPT) recently discovered in the modified embedded-atom model for elemental gallium. For this purpose we analyze the diffusive behavior in terms of the mean-squared displacement and self-intermediate scattering functions for two systems obtained by cooling the stable liquid through the LLPT at different cooling rates. The results show a pronounced heterogeneity of the dynamics upon the onset of the LLPT. Furthermore, it is found that this heterogeneity is closely correlated to the structural properties of the 9-fold coordinated high-density and 8-fold coordinated low-density liquid forms involved in the transition, showing a mixture of domains with very different diffusion time scales. The dynamics of the low-density liquid is found to be much more sluggish than that of the high-density form. Analysis of the energetics suggests that the origin of this difference is rooted in the fact that the cohesion in the former is significantly stronger than that in the latter.     

Equilibrium statistics of irregularly branched and dendrimeric polymeric liquid crystals

In this paper we develop a self-consistent model for the equilibrium statistics of nematic branched polymeric liquid crystals in the mean-field approximation. We have solved the resulting system of equations numerically and find a nematic-isotropic phase transition. We find that the order-disorder transition temperature scales as a function of the bond continuation probability, or equivalently the molecular weight, with an exponent that depends on the interaction potential. These results are compared with the experimentally observed behaviour.     

Equilibrium statistics of irregularly branched and dendrimeric polymeric liquid crystals

In this paper we develop a self-consistent model for the equilibrium statistics of nematic branched polymeric liquid crystals in the mean-field approximation. We have solved the resulting system of equations numerically and find a nematic-isotropic phase transition. We find that the order-disorder transition temperature scales as a function of the bond continuation probability, or equivalently the molecular weight, with an exponent that depends on the interaction potential. These results are compared with the experimentally observed behaviour.     

Phase behavior and local structure of a binary mixture in pores: mean-field lattice model calculations for analyzing neutron scattering data

We investigate the phase behavior of an asymmetric binary liquid A-W mixture confined between two planar homogenous substrates (slit pore). Molecules of species W interact preferentially with the solid walls via a long-range potential. Assuming nearest-neighbor attractions between the liquid molecules, we employ a lattice-gas model and a mean-field approximation for the grand potential. Minimization of this potential yields the density profiles of thermodynamically stable phases for fixed temperature, chemical potentials of both species, pore width and strengths of attraction. This model is used to analyze experimental small-angle neutron-scattering (SANS) data on the microscopic structure of the binary system isobutyric acid (iBA)+heavy water (D2O) inside a mesoscopic porous matrix (controlled-pore glass of about 10 nm mean pore width). Confinement-independent model parameters are adjusted so that the theoretical liquid-liquid coexistence curve in the bulk matches its experimental counterpart. By choosing appropriate values of the pore width and the attraction strength between substrates and water we analyze the effect of confinement on the phase diagram. In addition to a depression of the liquid-liquid critical point we observe surface induced phase transitions as well as water-film adsorption near the walls. The temperature dependence of the structure of water-rich and iBA-rich phases of constant composition are discussed in detail. The theoretical predictions are consistent with results of the SANS study and assist their interpretation.     

Does transparent nematic phase exist in 5CB∕DDAB∕water microemulsions? From the viewpoint of temperature dependent dielectric spectroscopy

Liquid crystal colloids have received tremendous attention because of its great potential both in the understanding of the liquid crystalline phase and in searching for new application of liquid crystals. Inverse microemulsion composed of 4-cyano-4-n-pentylbiphenyl (5CB), didodecyl dimethyl ammonium bromide, and water was investigated by means of broadband dielectric spectroscopy in this study. Based on the understanding of previous investigations on the same system, the isotropic phase was taken into account to quantitatively characterize the bulklike relaxations after the isotropic-to-nematic phase transition. Analogous results concerning the phase transition and phase composition to other investigations were obtained. In addition to bulklike relaxations, a new relaxation was observed at the frequency range about two orders lower than bulklike relaxations. This new relaxation shows abnormal temperature dependence, suggesting that superstructures composed of water droplets and confined 5CB molecules exist. This superstructure possibly possesses a confined nanoscaled liquid crystal ordering and may correspond to the notion of the transparent nematic phase.     

Effect of temperature on the interfacial behavior of a polystyrene-b-poly(methyl methacrylate) diblock copolymer at the air/water interface

Monolayers of a polystyrene-poly(methyl methacrylate) (PS-PMMA) diblock copolymer at the air-water interface were studied by measuring the surface pressure-area isotherms at several temperatures. Langmuir film balance experiments and atomic force microscopy showed that the diblock copolymer molecules formed surface micelles. In the plot of the surface pressure versus surface area per repeating unit, the monolayer changed from the gas phase to the liquid expanded phase at lower surface pressure for systems at low temperature compared to those at high temperature. In addition, a plateau, corresponding to the transition from the liquid expanded to liquid condensed phase, appeared in that plot at lower surface pressure for systems with a higher subphase (water) temperature. Hysteresis was observed in the compression-expansion cycle process. Increasing the subphase temperature alleviated this hyteresis gap, especially at low surface pressures. The minimum in the plot of the surface pressure versus surface area per repeating unit in the expansion process (which arises from the transition) and the transition plateau appeared more vividly at higher water temperature. These dynamic experimental results show that PS-PMMA diblock copolymers, in which both blocks are insoluble in water, do not form complicated entanglements in two-dimensional space. Although higher water temperature provided more entropy to the chains, and thus more conformational freedom, it did not change the surface morphology of the condensed film because both blocks of PS-PMMA are insoluble in water.     

The phase behavior, structure, and dynamics of rodlike mesogens with various flexibility using dissipative particle dynamics simulation

We present a systematic dissipative particle dynamics (DPD) study on the phase behavior, structure, and dynamics of rodlike mesogens. In addition to a rigid fused-bead-chain model with RATTLE constraint method, we also construct a semirigid model in which the flexibility is controlled by the bending constant of k(φ). Using this notation, the rigid model has an infinite bending constant of k(φ)=∞. Within the parameter space studied, both two kinds of models exhibit the nematic and smectic-A phases in addition to the isotropic and solid phases. All of the phase transitions are accompanied by the discontinuities in the thermodynamical, structural, and dynamical quantities and the hysteresis around the transition points, and are therefore first order. Note that the obtained solid state exhibits an in-layer tetragonal packing due to the high density. For the rigid model, the simulations show that the liquid crystal phases can be observed for mesogens with at least five beads and the nematic phase is the first one to appear. More importantly, the phase diagram of seven-bead-chain models is obtained as a function of k(φ) and temperature. It is found that decreasing the value of k(φ) reduces the anisotropy of molecular shape and the orientational ordering, and thereby shifts the liquid crystal phases to the lower temperature end of the phase diagram. Due to the different k(φ) dependence of phase transition temperatures, the nematic phase range exhibits a more marked narrowing than the smectic-A phase as k(φ) is reduced, implying that the flexibility has a destabilizing effect on the nematic and smectic-A phases. We also have investigated the anisotropic translational diffusion in liquid crystal phases and its temperature and flexibility dependence. In our study, we find that the phases formed, their statical and dynamic properties, as well as the transition properties are in close accord with those observations in real thermotropic liquid crystals. It is clear that both the rigid and semirigid models we used are valuable models with which to study the behavior of thermotropic liquid crystals using DPD algorithm.     

Dynamics of water confined in single- and double-wall carbon nanotubes

Using high-resolution quasielastic neutron scattering, we investigated the temperature dependence of single-particle dynamics of water confined in single- and double-wall carbon nanotubes with the inner diameters of 14+/-1 and 16+/-3 A, respectively. The temperature dependence of the alpha relaxation time for water in the 14 A nanotubes measured on cooling down from 260 to 190 K exhibits a crossover at 218 K from a Vogel-Fulcher-Tammann law behavior to an Arrhenius law behavior, indicating a fragile-to-strong dynamic transition in the confined water. This transition may be associated with a structural transition from a high-temperature, low-density (<1.02 gcm(3)) liquid to a low-temperature, high-density (>1.14 gcm(3)) liquid found in molecular dynamics simulation at about 200 K. However, no such dynamic transition in the investigated temperature range of 240-195 K was detected for water in the 16 A nanotubes. In the latter case, the dynamics of water simply follows a Vogel-Fulcher-Tammann law. This suggests that the fragile-to-strong crossover for water in the 16 A nanotubes may be shifted to a lower temperature.     

Glass transition of low-density amorphous water and related structures

In this work, the glass transition of water was studied with density functional theory. The transition temperature was determined by measuring the heat capacity Cp of low-density amorphous water during rapid heating. This technique ensures that all measurements were implemented without crystallization occurring, which is difficult to be achieved experimentally. The results showed that the glass transition occurs at 171 K, which is much higher than the reported value of 136 K. In addition, the triply hydrogen-bonded water molecules were found when T > 180 K, demonstrating the existence of the liquid structure at the higher temperature range.     

Binary mass diffusion in smectic C and CA phases as observed by forced Rayleigh scattering

We have observed the diffusion constants of a dye in several liquid crystals by forced Rayleigh scattering. In a liquid crystal which has a standard phase sequence of N-S_A-S_C, the diffusion anisotropy changes at the N-S_A phase transition and increases with decreasing temperature in S_A and S_C phases. The diffusion constants exhibit a rather smooth decrease with decreasing temperature except an anomaly at the S_A-S_c phase transition. In a liquid crystal which has the antiferroelectric S_CA phase, however, the diffusion constants show discontinuous increase and decrease at the S_A-S_c and the S_c-S_cA phase transition temperatures, respectively: the diffusion constant in S_C is larger than that in the higher temperature S_A phase. Anomalous signal increase and profile were observed at the phase transition temperatures, and were ascribed to the complementary transient grating due to the coexistence of two phases.     

Orientational dynamics for an amphiphilic-solvent solution

In this work, we performed Monte Carlo simulations on a lattice model for spontaneous amphiphilic aggregation, in order to study the orientational and hydrogen-bonding dynamics of water on different regions inside the micellar solution. We employed an associating lattice gas model that mimics the aqueous solvent, which presents a rich phase diagram with first- and second-order transition lines. Even though this is a simplified model, it makes possible to investigate the orientational dynamics of water in an equilibrium solution of amphiphiles, as well as the influence of the different phases of the solvent in the interfacial and bulk water dynamics. By means of extensive simulations, we showed that, at high temperatures, the behavior of the orientational relaxation and hydrogen bonding of water molecules in the bulk, first, and second hydration shells are considerable different. We observe the appearance of a very slow component for water molecules in the first hydration shell of micelles when the system reaches a high-density phase, consistent with previous theoretical and experimental studies concerning biological water. Also, at high temperatures, we find that water molecules in the second hydration shell of micelles have an orientational decay similar to that of bulk water, but with a generally slower dynamics. Otherwise, at low temperatures, we have two components for the orientational relaxation of bulk water in the low density liquid phase, and only a single component in the high density liquid (HDL) phase, which reflect the symmetry properties of the different phases of the solvent model. In the very dense region of water molecules in the first hydration shell of micelles at low temperatures, we find two components for the orientational relaxation on both liquid phases, one of them much slower than that in the single component of bulk water in the HDL phase. This happens even though our model does not present any hindrance to the water rotational freedom caused by the presence of the amphiphiles.     

Influence of non-mesogenic impurities on a nematic to isotropic phase transition

Non-mesomorphic solutes depress the normal nematic-isotropic transition temperature in liquid crystals. When non-mesomorphic solutes are added to a nematic liquid crystal, the nematic-isotropic transition temperature is depressed and a two phase region is formed due to the presence of impurities of the solutes. The present paper explains the formation of this two phase region by the Landau-de Gennes phenomenological theory, which agrees fairly well with the experimental observations. We also note that this two phase region indicates the tricritical behaviour of the nematic-isotropic phase transition and the phase diagram near the tricritical point is also obtained.     

Effect of Nanoscale Confinement on Freezing of Modified Water at Room Temperature and Ambient Pressure

Understanding the phase behavior of confined water is central to fields as diverse as heterogeneous catalysis, corrosion, nanofluidics, and to emerging energy technologies. Altering the state points (temperature, pressure, etc.) or introduction of a foreign surface can result in the phase transformation of water. At room temperature, ice nucleation is a very rare event and extremely high pressures in the GPa–TPa range are required to freeze water. Here, we perform computer experiments to artificially alter the balance between electrostatic and dispersion interactions between water molecules, and demonstrate nucleation and growth of ice at room temperature in a nanoconfined environment. Local perturbations in dispersive and electrostatic interactions near the surface are shown to provide the seed for nucleation (nucleation sites), which lead to room temperature liquid–solid phase transition of confined water. Crystallization of water occurs over several tens of nanometers and is shown to be independent of the nature of the substrate (hydrophilic oxide vs. hydrophobic graphene and crystalline oxide vs. amorphous diamond‐like carbon). Our results lead us to hypothesize that the freezing transition of confined water can be controlled by tuning the relative dispersive and electrostatic interaction.     

Direct Raman evidence for a weak continuous phase transition in liquid water

This paper presents the Raman depolarization ratio of degassed ultrapure water as a function of temperature, in the range 303.4-314.4 K (30.2-41.2 degrees C). The pressure of the sample was the vapor pressure of water at the measurement temperature. The data provide a direct indication of the existence of a phase transition in the liquid at 307.7 K, 5.8 kPa (34.6 degrees C, 0.057 atm). The minimum in the heat capacity, C(p)(), of water occurs at 34.5 degrees C, 1.0 atm (J. Res. Natl. Bur. Stand. 1939, 23, 197(1)). The minimum in C(p)() is shallow, and the transition is a weak-continuous phase transition. The pressure coefficient of the viscosity of water changes sign as pressure increases for temperatures below 35 degrees C (Nature 1965, 207, 620(2)). The viscosity minimum tracks the liquid phase transition in the P, T plane where it connects with the minimum in the freezing point of pure water in the same plane (Proc. Am. Acad. Arts Sci. 1911-12, 47, 441(3)). Previously we argued (J. Chem. Phys. 1998, 109, 7379(4)) that the minimum in the pressure coefficient of viscosity signaled the elimination of three-dimensional connectivity in liquid water. These observations coupled with recent measurements of the coordination shell of water near 300 K (Science 2004, 304, 995(5)) suggest that the structural component that changes during the phase transition is tetrahedrally coordinated water. At temperatures above the transition, there is no tetrahedrally coordinated water in the liquid and locally planar water structures dominate the liquid structure. Water is a structured liquid with distinct local structures that vary with temperature. Furthermore, liquid water has a liquid-liquid phase transition near the middle of the normal liquid range.