Polar nanoregions in water: A study of the dielectric properties of TIP4P/2005, TIP4P/2005f and TTM3F

Role of Dielectric Drag in Circumventing the Solubility-Diffusivity Trade-off in Zwitterionic Copolymer Membranes

Recent experiments have revealed that random zwitterionic amphiphilic copolymer (r-ZAC) membranes exhibit excellent Cl-/F- permselectivity circumventing the solubility-diffusivity trade-off. We conducted molecular dynamics simulations to investigate the origin of the experimental results on the transport of sodium halides in r-ZAC membranes. Our results indicate that the enhancement of Cl-/F- diffusivity selectivity in r-ZAC membranes (relative to that in bulk water) stems from the increase in dielectric drag dominating over the increase in Stokes drag, zwitterionic group-induced steric hindrance, and ion-polymer interactions. The importance of dielectric drag is further demonstrated by showing that reduction in ionic charges leads to a complete reversal of the diffusivity selectivity trends. We conclude that leveraging the impact of hydrophilic nanoconfinement on the dynamics of water can be utilized as a strategy to simultaneously augment solubility selectivity and diffusivity selectivity for separations, wherein the flux of the larger ionic species is desired over that of the smaller.

Isochores and Heat Capacity of Liquid Water in Terms of the Ion–Molecular Model

Thermodynamics of liquid water in terms of a non-standard approach—the ion–molecular model—is considered. Water is represented as a dense gas of neutral H2O molecules and single charged H3O+ and OH− ions. The molecules and ions perform thermal collisional motion and interconvert due to ion exchange. The energy-rich process—vibrations of an ion in a hydration shell of molecular dipoles—well known to spectroscopists with its dielectric response at 180 cm−1 (5 THz), is suggested to be key for water dynamics. Taking into account this ion–molecular oscillator, we compose an equation of state of liquid water to obtain analytical expressions for the isochores and heat capacity.

Weakly Anisotropic Dielectric Properties of Water Droplets at the Nanoscale

The slab-confined water at the nanoscale exhibits anomalous dielectric properties compared to bulk water, for example, significantly low dielectric constant. In this work, we study the dielectric properties of nanoscale water droplets at room temperature using molecular dynamics simulations. We find that the nanoscale water droplets feature weakly anisotropic dielectric constant: the radial component of dielectric constants is distinctly smaller than the tangential component although they both decrease with reducing droplet size in a similar way. Such dielectric behavior is closely related to the orientational preference of water molecules near the convex surface. The molecular dipole prefers to slightly orientate toward the interior of droplets in contrast to the out-of-plane preference for free-standing water films and slab-confined water, which suppresses the fluctuation of dipole moments in the radial direction. Meanwhile, it facilitates the formation of the open hydrogen-bond network in the surface layer and ultimately leads to the relatively weak suppression of tangential fluctuations. The differential suppression is responsible for the anisotropic dielectric constant of water droplets. This anisotropic characteristic is also found in dielectric relaxation: both the radial and the tangential relaxation are consistently slowed down upon approaching surface but the latter is universally slower.

Dielectric Profile and Electromelting of a Monolayer of Water Confined in Graphene Slit Pore

A monolayer of water confined between two parallel graphene sheets exists in many different phases and exhibits fascinating dielectric properties that have been studied in experiments. In this work, we use molecular dynamics simulations to study how the dielectric properties of a confined monolayer of water is affected by its structure. We consider six of the popular nonpolarizable water models-SPC/E, SPC/Fw, TIP3P, TIP3P_M (modified), TIP4P-2005, and TIP4P-2005f-and find that the in-plane structure of the water molecules at ambient temperature and pressure is strongly dependent on the water model: all the 3-point water models considered here show square ice formation, whereas no such structural ordering is observed for the 4-point water models. This allows us to investigate the role of the in-plane structure of the water monolayer on its dielectric profile. Our simulations show an anomalous perpendicular dielectric constant compared to the bulk, and the models that do not exhibit ice formation show very different dielectric response along the channel width compared to models that exhibit square ice formation. We also demonstrate the occurrence of electromelting of the in-plane ordered water under the application of a perpendicular electric field and find that the critical field for electromelting strongly depends on the water model. Together, we have shown the dependence of confined water properties on the different water structures that it may take when sandwiched between bilayer graphene. These remarkable properties of confined water can be exploited in various nanofluidic devices, artificial ion channels, and molecular sieving.

Nonlocal Dielectric Response of Water in Nanoconfinement

Recent experiments reporting a very low dielectric permittivity for nanoconfined water have renewed the interest in the structure and dielectric properties of water in narrow gaps. Here, we describe such systems with a minimal Landau-Ginzburg field theory composed of a nonlocal bulk-determined term and a local water-surface interaction term. We show how the interplay between the boundary conditions and intrinsic bulk correlations encodes the dielectric properties of confined water. Our theoretical analysis is supported by molecular dynamics simulations and comparison with the experimental data.

Dielectric ordering of water molecules arranged in a dipolar lattice

Intermolecular hydrogen bonds impede long-range (anti-)ferroelectric order of water. We confine H₂O molecules in nanosized cages formed by ions of a dielectric crystal. Arranging them in channels at a distance of ~5 Å with an interchannel separation of ~10 Å prevents the formation of hydrogen networks while electric dipole-dipole interactions remain effective. Here, we present measurements of the temperature-dependent dielectric permittivity, pyrocurrent, electric polarization and specific heat that indicate an order-disorder ferroelectric phase transition at T₀ ≈ 3 K in the water dipolar lattice. Ab initio molecular dynamics and classical Monte Carlo simulations reveal that at low temperatures the water molecules form ferroelectric domains in the ab-plane that order antiferroelectrically along the channel direction. This way we achieve the long-standing goal of arranging water molecules in polar order. This is not only of high relevance in various natural systems but might open an avenue towards future applications in biocompatible nanoelectronics.

Despite the apparent simplicity of a H2O molecule, the mutual ferroelectric ordering of the molecules is unresolved. Here, the authors realize a macroscopic ferroelectric phase transition in a network of dipole-dipole coupled water molecules located in nanopores of gemstone.

Exclusion Zone Phenomena in Water-A Critical Review of Experimental Findings and Theories

The existence of the exclusion zone (EZ), a layer of water in which plastic microspheres are repelled from hydrophilic surfaces, has now been independently demonstrated by several groups. A better understanding of the mechanisms which generate EZs would help with understanding the possible importance of EZs in biology and in engineering applications such as filtration and microfluidics. Here we review the experimental evidence for EZ phenomena in water and the major theories that have been proposed. We review experimental results from birefringence, neutron radiography, nuclear magnetic resonance, and other studies. Pollack theorizes that water in the EZ exists has a different structure than bulk water, and that this accounts for the EZ. We present several alternative explanations for EZs and argue that Schurr's theory based on diffusiophoresis presents a compelling alternative explanation for the core EZ phenomenon. Among other things, Schurr's theory makes predictions about the growth of the EZ with time which have been confirmed by Florea et al. and others. We also touch on several possible confounding factors that make experimentation on EZs difficult, such as charged surface groups, dissolved solutes, and adsorbed nanobubbles.

Nuclear Quantum Effects from the Analysis of Smoothed Trajectories: Pilot Study for Water

Nuclear quantum effects have significant contributions to thermodynamic quantities and structural properties; furthermore, very expensive methods are necessary for their accurate computation. In most calculations, these effects, for instance, zero-point energies, are simply neglected or only taken into account within the quantum harmonic oscillator approximation. Herein, we present a new method, Generalized Smoothed Trajectory Analysis, to determine nuclear quantum effects from molecular dynamics simulations. The broad applicability is demonstrated with the examples of a harmonic oscillator and different states of water. Ab initio molecular dynamics simulations have been performed for ideal gas up to the temperature of 5000 K. Classical molecular dynamics have been carried out for hexagonal ice, liquid water, and vapor at atmospheric pressure. With respect to the experimental heat capacity, our method outperforms previous calculations in the literature in a wide temperature range at lower computational cost than other alternatives. Dynamic and structural nuclear quantum effects of water are also discussed.

Electrically induced liquid-liquid phase transition in water at room temperature

In this work we expand on findings previously reported [Wexler et al., Phys. Chem. Chem. Phys., 2016, 18, 16281] on the experimental observation of a phase transition in a hydrogen bonded liquid manifesting in long range dipole-dipole interactions. The studied system, liquid water stressed by an electric field, exhibits collective oscillations brought about through spontaneous breakdown of symmetry. Raman spectroscopy identifies the primary excitation of the emergent phase as transverse optically active phonon-like sidebands that appear on the hydrogen bonded asymmetric stretch mode. The phase transition is observed throughout the entire volume of liquid. The system also exhibits a self-similarity relation between the scattered Raman intensity and the electric field strength which further supports the conclusion that collective behavior persists against thermal disruption. The experimental findings are discussed in terms of a quantum field theory for macroscopic quantum systems.

Enhancement of the local asymmetry in the hydrogen bond network of liquid water by an ultrafast electric field pulse

Condensed phase electron decomposition analysis based on density functional theory has recently revealed an asymmetry in the hydrogen-bond network in liquid water, in the sense that a significant population of water molecules are simultaneously donating and accepting one strong hydrogen-bond and another substantially weaker one. Here we investigate this asymmetry, as well as broader structural and energetic features of water's hydrogen-bond network, following the application of an intense electric field square pulse that invokes the ultrafast reorientation of water molecules. We find that the necessary field-strength required to invoke an ultrafast alignment in a picosecond time window is on the order of 108 Vm-1. The resulting orientational anisotropy imposes an experimentally measurable signature on the structure and dynamics of the hydrogen-bond network, including its asymmetry, which is strongly enhanced. The dependence of the molecular reorientation dynamics on the field-strength can be understood by relating the magnitude of the water dipole-field interaction to the rotational kinetic energy, as well as the hydrogen-bond energy.

The effect of polymorphism on the structural, dynamic and dielectric properties of plastic crystal water: A molecular dynamics simulation perspective

We have employed molecular dynamics simulations based on the TIP4P/2005 water model to investigate the local structural, dynamical, and dielectric properties of the two recently reported body-centered-cubic and face-centered-cubic plastic crystal phases of water. Our results reveal significant differences in the local orientational structure and rotational dynamics of water molecules for the two polymorphs. The probability distributions of trigonal and tetrahedral order parameters exhibit a multi-modal structure, implying the existence of significant local orientational heterogeneities, particularly in the face-centered-cubic phase. The calculated hydrogen bond statistics and dynamics provide further indications of the existence of a strongly heterogeneous and rapidly interconverting local orientational structural network in both polymorphs. We have observed a hindered molecular rotation, much more pronounced in the body-centered-cubic phase, which is reflected by the decay of the fourth-order Legendre reorientational correlation functions and angular Van Hove functions. Molecular rotation, however, is additionally hindered in the high-pressure liquid compared to the plastic crystal phase. The results obtained also reveal significant differences in the dielectric properties of the polymorphs due to the different dipolar orientational correlation characterizing each phase.

Quantitative High-Dynamic-Range Electron Diffraction of Polar Nanodomains in Pb2 ScTaO6

Highly B-site ordered Pb₂ ScTaO₆ crystals are studied as a function of temperature via dielectric spectroscopy and in situ high-dynamic-range electron diffraction. The degree of ordering is examined on the local and macroscopic scale and is determined to be 76%. Novel analysis of the electron diffraction patterns provides structural information with two types of antiferroelectric displacements determined to be present in the polar structure. It is then found that a low-temperature transition occurs on cooling at ≈210 K that is not present on heating. This phenomenon is discussed in terms of the freezing of dynamic polar nanodomains where a high density of domain walls creates a metastable state.

Protein-protein interactions within photosystem II under photoprotection: the synergy between CP29 minor antenna, subunit S (PsbS) and zeaxanthin at all-atom resolution

The assembly and disassembly of protein complexes within cells are crucial life-sustaining processes. In photosystem II (PSII) of higher plants, there is a delicate yet obscure balance between light harvesting and photo-protection under fluctuating light conditions, that involves protein-protein complexes. Recent breakthroughs in molecular dynamics (MD) simulations are combined with new approaches herein to provide structural and energetic insight into such a complex between the CP29 minor antenna and the PSII subunit S (PsbS). The microscopic model involves extensive sampling of bound and dissociated states at atomic resolution in the presence of photo-protective zeaxanthin (Zea), and reveals well defined protein-protein cross-sections. The complex is placed within PSII, and macroscopic connections are emerging (PsbS-CP29-CP24-CP47) along the energy transfer pathways from the antenna to the PSII core. These connections explain macroscopic observations in the literature, while the previously obscured atomic scale details are now revealed. The implications of these findings are discussed in the context of the Non-Photochemical Quenching (NPQ) of chlorophyll fluorescence, the down-regulatory mechanism of photosynthesis, that enables the protection of PSII against excess excitation load. Zea is found at the PsbS-CP29 cross-section and a pH-dependent equilibrium between PsbS dimer/monomers and the PsbS-CP29 dissociation/association is identified as the target for engineering tolerant plants with increased crop and biomass yields. Finally, the new MD based approaches can be used to probe protein-protein interactions in general, and the PSII structure provided can initiate large scale molecular simulations of the photosynthetic apparatus, under NPQ conditions.

The Photosystem II Subunit S under Stress

Nonphotochemical quenching is the protective mechanism against overexcitation of photosystem II, triggered by excess ΔpH in photosynthetic membranes. The light-harvesting complexes (LHCs), the de-epoxidation of violaxanthin to zeaxanthin, and the photosystem II subunit S (PsbS) work in synergy for an optimized multilevel response. Understanding the fine details of this synergy has proven challenging to scientific research. Here, we employ large-scale, all-atom molecular simulations and beyond experimental insight, we proceed a step further in identifying the PsbS dynamics that could possibly be associated with this synergy. For the first time, to our knowledge, we probe the distinct behavior of PsbS under ΔpH that probes the details of the potential dimer-to-monomer transition, and in a violaxanthin/zeaxanthin-rich membrane, at an all-atom resolution. We propose that the lumen-exposed residues, threonine 162 and glutamic acid 173, form stabilizing hydrogen bonds between the PsbS monomers only at high lumen pH, whereas at low pH (excess ΔpH) this interaction is lost, and leads to higher flexibility of the protein and potentially to the dimer-to-monomer transition. Lastly, we discuss how conformational changes under the presence of ΔpH/zeaxanthin are related to the PsbS role in the current nonphotochemical quenching model in the literature. For the latter, we probe a PsbS-monomeric LHCII association. The association is proposed to potentially alter the monomeric LHCII sensitivity to ΔpH by changing the pKa values of interacting LHCII residues. This serves as an example where protonation-ligation events enhance protein-protein interactions fundamental to many life processes.

How Water's Properties Are Encoded in Its Molecular Structure and Energies

How are water's material properties encoded within the structure of the water molecule? This is pertinent to understanding Earth's living systems, its materials, its geochemistry and geophysics, and a broad spectrum of its industrial chemistry. Water has distinctive liquid and solid properties: It is highly cohesive. It has volumetric anomalies-water's solid (ice) floats on its liquid; pressure can melt the solid rather than freezing the liquid; heating can shrink the liquid. It has more solid phases than other materials. Its supercooled liquid has divergent thermodynamic response functions. Its glassy state is neither fragile nor strong. Its component ions-hydroxide and protons-diffuse much faster than other ions. Aqueous solvation of ions or oils entails large entropies and heat capacities. We review how these properties are encoded within water's molecular structure and energies, as understood from theories, simulations, and experiments. Like simpler liquids, water molecules are nearly spherical and interact with each other through van der Waals forces. Unlike simpler liquids, water's orientation-dependent hydrogen bonding leads to open tetrahedral cage-like structuring that contributes to its remarkable volumetric and thermal properties.

A pathway for protective quenching in antenna proteins of Photosystem II

Photosynthesis is common in nature, converting sunlight energy into proton motive force and reducing power. The increased spectral range absorption of light exerted by pigments (i.e. chlorophylls, Chls) within Light Harvesting Complexes (LHCs) proves an important advantage under low light conditions. However, in the exposure to excess light, oxidative damages and ultimately cell death can occur. A down-regulatory mechanism, thus, has been evolved (non-photochemical quenching, NPQ). The mechanistic details of its major component (qE) are missing at the atomic scale. The research herein, initiates on solid evidence from the current NPQ state of the art, and reveals a detailed atomistic view by large scale Molecular Dynamics, Metadynamics and ab initio Simulations. The results demonstrate a complete picture of an elaborate common molecular design. All probed antenna proteins (major LHCII from spinach-pea, CP29 from spinach) show striking plasticity in helix-D, under NPQ conditions. This induces changes in Qy bands in excitation and absorption spectra of the near-by pigment pair (Chl613-614) that could emerge as a new quenching site. Zeaxanthin enhances this plasticity (and possibly the quenching) even at milder NPQ conditions.

The origin of the Debye relaxation in liquid water and fitting the high frequency excess response

Debye relaxation in the microwave/THz region of the dielectric spectra of water is related to propagation of defects.

The hydrogen-bond network of water supports propagating optical phonon-like modes

The local structure of liquid water as a function of temperature is a source of intense research. This structure is intimately linked to the dynamics of water molecules, which can be measured using Raman and infrared spectroscopies. The assignment of spectral peaks depends on whether they are collective modes or single-molecule motions. Vibrational modes in liquids are usually considered to be associated to the motions of single molecules or small clusters. Using molecular dynamics simulations, here we find dispersive optical phonon-like modes in the librational and OH-stretching bands. We argue that on subpicosecond time scales these modes propagate through water's hydrogen-bond network over distances of up to 2 nm. In the long wavelength limit these optical modes exhibit longitudinal–transverse splitting, indicating the presence of coherent long-range dipole–dipole interactions, as in ice. Our results indicate the dynamics of liquid water have more similarities to ice than previously thought.

How the local structure of water varies as a function of temperature is a long-studied topic, which is still under debate. Here, the authors show that dielectric susceptibility measurements might be used to probe and identify propagating optical phonon-like modes in the hydrogen-bond network of water.

Dielectric constant and low-frequency infrared spectra for liquid water and ice Ih within the E3B model

Two intrinsic difficulties in modeling condensed-phase water with conventional rigid non-polarizable water models are: reproducing the static dielectric constants for liquid water and ice Ih, and generating the peak at about 200 cm(-1) in the low-frequency infrared spectrum for liquid water. The primary physical reason for these failures is believed to be the missing polarization effect in these models, and consequently various sophisticated polarizable water models have been developed. However, in this work we pursue a different strategy and propose a simple empirical scheme to include the polarization effect only on the dipole surface (without modifying a model's intermolecular interaction potential). We implement this strategy for our explicit three-body (E3B) model. Our calculated static dielectric constants and low-frequency infrared spectra are in good agreement with experiment for both liquid water and ice Ih over wide temperature ranges, albeit with one fitting parameter for each phase. The success of our modeling also suggests that thermal fluctuations about local minima and the energy differences between different proton-disordered configurations play minor roles in the static dielectric constant of ice Ih. Our analysis shows that the polarization effect is important in resolving the two difficulties mentioned above and sheds some light on the origin of several features in the low-frequency infrared spectra for liquid water and ice Ih.

Fluctuations of local electric field and dipole moments in water between metal walls

We examine the thermal fluctuations of the local electric field Ek (loc) and the dipole moment μk in liquid water at T = 298 K between metal walls in electric field applied in the perpendicular direction. We use analytic theory and molecular dynamics simulation. In this situation, there is a global electrostatic coupling between the surface charges on the walls and the polarization in the bulk. Then, the correlation function of the polarization density pz(r) along the applied field contains a homogeneous part inversely proportional to the cell volume V. Accounting for the long-range dipolar interaction, we derive the Kirkwood-Fröhlich formula for the polarization fluctuations when the specimen volume v is much smaller than V. However, for not small v/V, the homogeneous part comes into play in dielectric relations. We also calculate the distribution of Ek (loc) in applied field. As a unique feature of water, its magnitude |Ek (loc)| obeys a Gaussian distribution with a large mean value E0 ≅ 17 V/nm, which arises mainly from the surrounding hydrogen-bonded molecules. Since |μk|E0 ∼ 30kBT, μk becomes mostly parallel to Ek (loc). As a result, the orientation distributions of these two vectors nearly coincide, assuming the classical exponential form. In dynamics, the component of μk(t) parallel to Ek (loc)(t) changes on the time scale of the hydrogen bonds ∼5 ps, while its smaller perpendicular component undergoes librational motions on time scales of 0.01 ps.

Long-range orientation correlation in water

Strong short-range intermolecular interactions result in position and orientation correlations between nearest neighbour molecules in isotropic liquids, but it is generally assumed that such correlations extend at most a few molecular diameters. Results from second-harmonic light scattering experiments presented here reveal long-range molecular orientation correlations in liquid water, where the molecular dipole orientation distribution has the form of a nearly pure transverse vector field. Spatial scales in the range 200-2000 nm are probed by the angle-dependent measurements and the observed correlations are thought to result from rotation-translation coupling in acoustic phonons in the liquid.