Hydrogen Bonding and Related Properties in Liquid Water: A Car–Parrinello Molecular Dynamics Simulation Study

Terahertz spectroscopy as a method for investigation of hydration shells of biomolecules

The hydration of biomolecules is one of the fundamental processes underlying the construction of living matter. The formation of the native conformation of most biomolecules is possible only in an aqueous environment. At the same time, not only water affects the structure of biomolecules, but also biomolecules affect the structure of water, forming hydration shells. However, the study of the structure of biomolecules is given much more attention than their hydration shells. A real breakthrough in the study of hydration occurred with the development of the THz spectroscopy method, which showed that the hydration shell of biomolecules is not limited to 1-2 layers of strongly bound water, but also includes more distant areas of hydration with altered molecular dynamics. This review examines the fundamental features of the THz frequency range as a source of information about the structural and dynamic characteristics of water that change during hydration. The applied approaches to the study of hydration shells of biomolecules based on THz spectroscopy are described. The data on the hydration of biomolecules of all main types obtained from the beginning of the application of THz spectroscopy to the present are summarized. The emphasis is placed on the possible participation of extended hydration shells in the realization of the biological functions of biomolecules and at the same time on the insufficient knowledge of their structural and dynamic characteristics.

Hydrogen-Bond Configurations of Hydration Water around Glycerol Investigated by HOH Bending and OH Stretching Analysis

Toward a comprehensive understanding of the mechanism of glycerol as a moisturizer, studies on the hydrogen-bond (HB) structure of hydration water, which is known to be disordered by glycerol, are insufficient. To this aim, we evaluated the HB configurations based on the HOH bending and OH stretching spectra of the hydration water from those of glycerol/water mixtures by subtracting the contributions of bulk water and glycerol using dielectric relaxation spectroscopy. Analysis of the HOH bending band showed that hydration water-donating HBs lose the intermolecular bending coupling with increasing glycerol by replacing the water-water HBs with water-glycerol HBs. The OH stretching band provided more detailed insight into the HB configuration, indicating that the double-donor double-acceptor and double-donor single-acceptor configurations in bulk water change to a predominantly double-donor single-acceptor configuration in hydration water around glycerol. The formation of more donor HBs than acceptor HBs may be due to the steric constrains by glycerol and/or differences in the partial charge on the oxygen atom between water and glycerol.

Broadband Dielectric Spectroscopic Analysis toward Characterization of the Hydration State and Bioprotective Superiority of Trehalose

Alteration of the hydrogen-bond (H-bond) network by trehalose is acknowledged as a bioprotective agent. However, most studies exploring the hydration superiority of the trehalose structure are limited structure are limited by the computational cost or a narrow-range spectrum. In the present study, the structural and dynamical behaviors of the H-bond network of trehalose and maltose solutions were observed and compared with a broadband dielectric spectrum (100 MHz-18 THz) to investigate the influence of the trehalose structure on the bioprotective function. From the relaxation time, the reorientation cooperativity, resonant frequency, and damping constant of water-water vibration, the symmetric structure of trehalose allowed a more significant H-bond strengthening effect and homogeneous aqueous environment. In contrast, the difference in the hydration number between trehalose and maltose was negligible. Thus, the enhanced H-bond strengthening effect and homogeneous aqueous environment owing to the symmetric structure are the essential factors that contribute to the remarkable bioprotective effect of trehalose.

The Proton in Biochemistry: Impacts on Bioenergetics, Biophysical Chemistry, and Bioorganic Chemistry

The proton is the smallest atomic particle, and in aqueous solution it is the smallest hydrated ion, having only two waters in its first hydration shell. In this article we survey key aspects of the proton in chemistry and biochemistry, starting with the definitions of pH and pK _a and their application inside biological cells. This includes an exploration of pH in nanoscale spaces, distinguishing between bulk and interfacial phases. We survey the Eigen and Zundel models of the structure of the hydrated proton, and how these can be used to explain: a) the behavior of protons at the water-hydrophobic interface, and b) the extraordinarily high mobility of protons in bulk water via Grotthuss hopping, and inside proteins via proton wires. Lastly, we survey key aspects of the effect of proton concentration and proton transfer on biochemical reactions including ligand binding and enzyme catalysis, as well as pH effects on biochemical thermodynamics, including the Chemiosmotic Theory. We find, for example, that the spontaneity of ATP hydrolysis at pH ≥ 7 is not due to any inherent property of ATP (or ADP or phosphate), but rather to the low concentration of H⁺. Additionally, we show that acidification due to fermentation does not derive from the organic acid waste products, but rather from the proton produced by ATP hydrolysis.

Highly Reversible Aqueous Zinc Batteries enabled by Zincophilic-Zincophobic Interfacial Layers and Interrupted Hydrogen-Bond Electrolytes

Aqueous Zn batteries promise high energy density but suffer from Zn dendritic growth and poor low-temperature performance. Here, we overcome both challenges by using an eutectic 7.6 m ZnCl₂ aqueous electrolyte with 0.05 m SnCl₂ additive, which in situ forms a zincophilic/zincophobic Sn/Zn₅ (OH)₈ Cl₂ ⋅H₂ O bilayer interphase and enables low temperature operation. Zincophilic Sn decreases Zn plating/stripping overpotential and promotes uniform Zn plating, while zincophobic Zn₅ (OH)₈ Cl₂ ⋅H₂ O top-layer suppresses Zn dendrite growth. The eutectic electrolyte has a high ionic conductivity of ≈0.8 mS cm^-1 even at -70 °C due to the distortion of hydrogen bond network by solvated Zn²⁺ and Cl^- . The eutectic electrolyte enables Zn∥Ti half-cell a high Coulombic efficiency (CE) of >99.7 % for 200 cycles and Zn∥Zn cell steady charge/discharge for 500 h with a low overpotential of 8 mV at 3 mA cm^-2 . Practically, Zn∥VOPO₄ batteries maintain >95 % capacity with a CE of >99.9 % for 200 cycles at -50 °C, and retain ≈30 % capacity at -70 °C of that at 20 °C.

Aqueous Affinity and Interfacial Dynamics of Anisotropic Buckled Black Phosphorous

The structure of black phosphorous (BP) is similar to the honeycomb arrangement of graphene, but the layered BP is found to be buckled and highly anisotropic. The buckled surface structure affects interfacial molecule mobility and plays a vital role in various nanomaterial applications. The BP is also known for wettability, droplet formation, stability, and hydrophobicity in the aqueous environment. However, there is a gap concerning the structural and dynamical behavior of water molecules, which is available in abundance for other monoatomic and polyatomic two-dimensional (2D) materials. Motivated by the technological importance, we try to bridge the gap by explaining the surface anisotropy-facilitated behavior of water molecules on bilayer BP using classical and first principles molecular dynamics (MD) simulations. From our classical MD study, we find three distinct layers of water molecules. The water layer closest to the interface is L1, followed by L2 and L3/bulk perpendicular to the BP surface. Water molecules in the L1 layer experience some structural disintegration in hydrogen bond (HB) phenomena compared to the bulk. There is a loss of HB donor-acceptor count per water molecule. The average HB count decreases because of an elevated rate of HB formation and deformation; this would affect the dynamic properties in terms of HB lifetime. Therefore, we observe the reduced lifetime of HB in the layer in close contact with BP, which again complements our finding on the diffusion coefficient of water molecules in distinct layers. Water diffuses relatively faster with diffusion coefficient 3.25 × 10^-9 m² s^-1 in L1, followed by L2 and L3. The BP layer shows moderate hydrophobic nature. Our results also indicate the anisotropic behavior as the diffusion along the x-direction is faster than that along the y-direction. The gap in the slope of the x and y components of mean-squared displacement (MSD) complements the pinning effect in an aqueous environment. We observe blue-shifted and red-shifted libration and O-H stretching modes from the calculated power spectra for the L1 water molecules compared to the L2 and L3 molecules from first principles MD simulations. Our analysis may help understand the physical phenomena that occur during the surface wetting of the predroplet formation process observed experimentally.

Graphene Oxide as a Pb(II) Separation Medium: Has Part of the Story Been Overlooked?

A key problem associated with the design of graphene oxide (GO) materials and their tuning for nanoscale separations is how specific functional groups influence the competitive adsorption of solvated ions and water at liquid/graphene interfaces. Computation accompanied by experiment shows that OH and COOH exert an influence on water adsorption properties stronger than that of O and H functional groups. The COO- anions, following COOH deprotonation, stabilize Pb(II) through strong electrostatic interactions. This suggests that, among the functional groups under study, COOH offers the best Pb(II) adsorption capacity and the ability to regenerate the sorbent through a pH swing. In line with computation, striking experimental observations revealed that a substantial increase in Pb(II) adsorption occurs with increasing pH. Our findings provide a systematic framework for controlled design and implementation of regenerable C-based sorbents used in separations and desalination.

Conformation-induced vibrational spectral dynamics of hydrogen peroxide and vicinal water molecules

We studied the conformation-induced spectral response of water molecules due to site-specific structural alterations of solvated hydrogen peroxide (H₂O₂) employing DFT-based first principles molecular dynamics (FPMD) simulations. Wavelet transform was used to determine the time-dependent frequencies of the hydroxyls of water molecules and the O-H stretch modes of H₂O₂. Shifts in the vibrational frequency of the hydrogen-bonded hydroxyls inside the solvation shell of H₂O₂ support multiple distinctive hydrogen bonding environments. This paper classifies two distinct hydrogen bond types inside the O-O_W solvation shell of H₂O₂, and the dynamical calculations provide a quantitative estimation of the relative hydrogen bond strength. We ascertain the reason for not observing the escape of water molecules from the hydrogen peroxide hydration shell, unlike the solvation shell of ionic solutions and neutral solutes. Besides, we provide a comprehensive analysis of the spectral shifts in the normalized frequency distribution, the time-dependent decay of frequency-frequency correlation functions, and the hydrogen bond length scale fluctuations. We also quantify the relative contribution of the cisoid and transoid conformers affecting the vibrational spectral signature of the vicinal water molecules. While the transoid conformers promote the hydrogen bonding interactions through the oxygen site (OH_W), the cisoid conformers facilitate hydrogen peroxide-water hydrogen bond formation through the hydrogen site (HO_W). These non-identical hydrogen bond associations stabilize hydrogen peroxide in water.

Dynamical Collectivity and Nuclear Quantum Effects on the Intermolecular Stretching Mode of Liquid Water

This study investigated the broadband terahertz and low-frequency Raman spectroscopy of liquid water (H₂O, D₂O, and H₂¹⁸O) over 2 decades of frequency to address long-standing challenges regarding the interpretation of the intermolecular stretching mode at around 5 THz. We experimentally demonstrated that the intermolecular stretching mode of liquid water obtained via terahertz spectroscopy is significantly redshifted and broadened compared with that via Raman. This result was rationalized by the enhanced dynamical collectivity probed by terahertz spectroscopy, although both have a common origin in the kinetic motion. Their temperature and isotope dependences emphasize the significance of oscillation mass in determining the intermolecular stretching lineshape, while quantum effects cannot be overlooked in both terahertz and low-frequency Raman spectra.

Hydrogen-Bond Structure and Low-Frequency Dynamics of Electrolyte Solutions: Hydration Numbers from ab Initio Water Reorientation Dynamics and Dielectric Relaxation Spectroscopy

We present an atomistic simulation scheme for the determination of the hydration number (h) of aqueous electrolyte solutions based on the calculation of the water dipole reorientation dynamics. In this methodology, the time evolution of an aqueous electrolyte solution generated from ab initio molecular dynamics simulations is used to compute the reorientation time of different water subpopulations. The value of h is determined by considering whether the reorientation time of the water subpopulations is retarded with respect to bulk-like behavior. The application of this computational protocol to magnesium chloride (MgCl2 ) solutions at different concentrations (0.6-2.8 mol kg-1 ) gives h values in excellent agreement with experimental hydration numbers obtained using GHz-to-THz dielectric relaxation spectroscopy. This methodology is attractive because it is based on a well-defined criterion for the definition of hydration number and provides a link with the molecular-level processes responsible for affecting bulk solution behavior. Analysis of the ab initio molecular dynamics trajectories using radial distribution functions, hydrogen bonding statistics, vibrational density of states, water-water hydrogen bonding lifetimes, and water dipole reorientation reveals that MgCl2 has a considerable influence on the hydrogen bond network compared with bulk water. These effects have been assigned to the specific strong Mg-water interaction rather than the Cl-water interaction.

Solvation structure and dynamics of the dimethylammonium cation diluted in liquid water: A molecular dynamics approach

Classical molecular dynamics simulation techniques were employed to investigate the local solvation structure and related dynamics of the dimethylammonium cation diluted in liquid water at ambient conditions. The translational and orientational order around the dimethylammonium cation was investigated in terms of the corresponding radial and angular distribution functions. The results obtained revealed that the first solvation shell of the dimethylammonium consists mainly of two and, less frequently, three water molecules. The two nearest water neighbors form hydrogen bonds with the ammonium hydrogen atoms of the cation, whereas the third neighbor interacts with the methyl hydrogen atoms as well. The distribution of the trigonal order parameter exhibits a bimodal behavior, signifying the existence of local orientational heterogeneities in the solvation shell of the dimethylammonium cation. The calculated continuous and intermittent residence and hydrogen bond lifetimes for the cation-water pairs have also been found to be longer in comparison with the water-water ones. The very similar self-diffusion coefficients of the dimethylammonium cation and the water molecules in the bulk dilute solution indicate that the translational motions of the cation are mainly controlled by the translational mobility of the surrounding water molecules.

Water Distribution and Clustering on the Lyophilized IgG1 Surface: Insight from Molecular Dynamics Simulations

Water has a critical role in the stability of the higher-order structure of proteins. In addition, it is considered to be a major destabilization factor for the physical and chemical stability of freeze-dried proteins and peptides. Physical and chemical aspects of protein/water relationships are commonly studied with the use of water vapor sorption isotherms for amorphous lyophilized proteins, which, in turn, are commonly analyzed using the Brunauer-Emmett-Teller (BET) equation to obtain the parameters, W_m and C_B. The parameter W_m is generally referred to as the "monolayer limit of adsorption" and has a narrow range of 6-8% for most proteins. In this study, the water distribution on an IgG1 surface is investigated by molecular dynamics (MD) simulations at different water contents. The monolayer of water molecules was found to have limited coverage of the protein surface, and the true monolayer coverage of the protein globule actually occurs at a hydration level above 30%. The distribution of water molecules on the IgG1 surface is also highly heterogeneous, and the heterogeneity is not considered in the BET theory. In this study, a mechanistic model has been developed to describe the water vapor sorption isotherm. This model is based on the analysis of the hydrogen bonding network extracted from the MD simulations. The model is consistent with the experimental Type-II isotherm, which is usually observed for proteins. The physical meaning of the BET monolayer was redefined as the onset of water cluster formation. A simple model to calculate the onset water level, W_m, is proposed based on the hydration of different amino acids, as determined from the MD simulations.

Computational Analysis of Energy Landscapes Reveals Dynamic Features That Contribute to Binding of Inhibitors to CFTR-Associated Ligand

The CFTR-associated ligand PDZ domain (CALP) binds to the cystic fibrosis transmembrane conductance regulator (CFTR) and mediates lysosomal degradation of mature CFTR. Inhibition of this interaction has been explored as a therapeutic avenue for cystic fibrosis. Previously, we reported the ensemble-based computational design of a novel peptide inhibitor of CALP, which resulted in the most binding-efficient inhibitor to date. This inhibitor, kCAL01, was designed using osprey and evinced significant biological activity in in vitro cell-based assays. Here, we report a crystal structure of kCAL01 bound to CALP and compare structural features against iCAL36, a previously developed inhibitor of CALP. We compute side-chain energy landscapes for each structure to not only enable approximation of binding thermodynamics but also reveal ensemble features that contribute to the comparatively efficient binding of kCAL01. Finally, we compare the previously reported design ensemble for kCAL01 vs the new crystal structure and show that, despite small differences between the design model and crystal structure, significant biophysical features that enhance inhibitor binding are captured in the design ensemble. This suggests not only that ensemble-based design captured thermodynamically significant features observed in vitro, but also that a design eschewing ensembles would miss the kCAL01 sequence entirely.

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.

A Water Solvation Shell Can Transform Gold Metastable Nanoparticles in the Fluxional Regime

Solvated gold nanoparticles have been modeled in the fluxional regime by density functional theory including dispersion forces for an extensive set of conventional morphologies. The study of isolated adsorption of one water molecule shows that the most stable adsorption forms are similar (corners and edges) regardless of the nanoparticle shape and size, although the adsorption strength differs significantly (0.15 eV). When a complete and explicit water solvation shell interacts with gold nanoclusters, metastable in vacuum and presenting a predominance of (100) square facets (ino-decahedra Au₅₅ and Au₁₄₇), these nanoparticles are found unstable and transform into the closest morphologies exhibiting mainly (111) triangular facets and symmetries. The corresponding adsorption strength per water molecule becomes independent of shape and size and is enhanced by the formation of two hydrogen bonds on average. For applications in radiotherapy, this study suggests that the shapes of small gold nanoparticles should be homogenized by interacting with the biological environment.

Understanding the effect of nanoconfinement on the structure of water hydrogen bond networks

Dynamic hydrogen bond trails in water confined between two phospholipid membranes traced by the information flow model.

Using a monomer potential energy surface to perform approximate path integral molecular dynamics simulation of ab initio water at near-zero added cost

It is now established that nuclear quantum motion plays an important role in determining water's hydrogen bonding, structure, and dynamics. Such effects are important to include in density functional theory (DFT) based molecular dynamics simulation of water. The standard way of treating nuclear quantum effects, path integral molecular dynamics (PIMD), multiplies the number of energy/force calculations by the number of beads required. In this work we introduce a method whereby PIMD can be incorporated into a DFT simulation with little extra cost and little loss in accuracy. The method is based on the many body expansion of the energy and has the benefit of including a monomer level correction to the DFT energy. Our method calculates intramolecular forces using the highly accurate monomer potential energy surface developed by Partridge-Schwenke, which is cheap to evaluate. Intermolecular forces and energies are calculated with DFT only once per timestep using the centroid positions. We show how our method may be used in conjunction with a multiple time step algorithm for an additional speedup and how it relates to ring polymer contraction and other schemes that have been introduced recently to speed up PIMD simulations. We show that our method, which we call "monomer PIMD", correctly captures changes in the structure of water found in a full PIMD simulation but at much lower computational cost.

Reconsideration of the relaxational and vibrational line shapes of liquid water based on ultrabroadband dielectric spectroscopy

The Debye relaxation function is widely used to describe the large dielectric dispersion of ambient water around 20 GHz. However, from a theoretical point of view, this function is supposed to give incorrect predictions at high frequencies owing to the inappropriate assumption that inertial effects and intermolecular interactions do not affect the relaxation dynamics. Our ultrabroadband spectroscopy investigation of liquid water ranging from 500 MHz to 400 THz did demonstrate that the Debye function is inaccurate far above the microwave region. As an alternative, we tried a stochastic frequency modulation (SFM) model assuming instantaneous modification of the line shapes by the correlation with the surrounding system. The SFM relaxation model reproduced the experimental dielectric spectra up to 400 THz, showing that the hydrogen-bond dynamics are associated with the inertial effect that causes the non-exponential relaxation behaviour in a very short time (typically 25 fs). Within the framework of this relaxation model, the hindered translation modes are able to be approximated as fast modulation (homogeneous) line shapes because the interaction time with frequency modulation is too short. Compared with them, the libation mode is found to have a relatively slow modulation (inhomogeneous) origin, where disturbance of water hydrogen bonds induced by the hindered translations leads to fluctuations in the libration frequency.

Characterization of the hydrogen-bond network of water around sucrose and trehalose: Microwave and terahertz spectroscopic study

Modification of the water hydrogen bond network imposed by disaccharides is known to serve as a bioprotective agent in living organisms, though its comprehensive understanding is still yet to be reached. In this study, aiming to characterize the dynamical slowing down and destructuring effect of disaccharides, we performed broadband dielectric spectroscopy, ranging from 0.5 GHz to 12 THz, of sucrose and trehalose aqueous solutions. The destructuring effect was examined in two ways (the hydrogen bond fragmentation and disordering) and our result showed that both sucrose and trehalose exhibit an obvious destructuring effect with a similar strength, by fragmenting hydrogen bonds and distorting the tetrahedral-like structure of water. This observation strongly supports a chaotropic (structure-breaking) aspect of disaccharides on the water structure. At the same time, hydration water was found to exhibit slower dynamics and a greater reorientational cooperativity than bulk water because of the strengthened hydrogen bonds. These results lead to the conclusion that strong disaccharide-water hydrogen bonds structurally incompatible with native water-water bonds lead to the rigid but destructured hydrogen bond network around disaccharides. Another important finding in this study is that the greater dynamical slowing down of trehalose was found compared with that of sucrose, at variance with the destructuring effect where no solute dependent difference was observed. This discovery suggests that the exceptionally greater bioprotective impact especially of trehalose among disaccharides is mainly associated with the dynamical slowing down (rather than the destructuring effect).

Peculiar Molecular Shape and Size Dependence of the Dynamics of Fluids Confined in a Small-Pore Metal-Organic Framework

Force-field-based molecular dynamics simulations were deployed to systematically explore the dynamics of confined molecules of different shapes and sizes, that is, linear (CO₂ and N₂) and spherical (CH₄) fluids, in a model small pore system, that is, the metal-organic framework SIFSIX-2-Cu-i. These computations unveil an unprecedented molecular symmetry dependence of the translational and rotational dynamics of fluids confined in channel-like nanoporous materials. In particular, this peculiar behavior is reflected by the extremely slow decay of the Legendre reorientational correlation functions of even-parity order for the linear fluids, which is associated with jump-like orientation flips, while the spherical fluid shows a very fast decay taking place on a subpicosecond time scale. Such a fundamental understanding is relevant to diverse disciplines such as in chemistry, physics, biology, and materials science, where diatomic or polyatomic molecules of different shapes/sizes diffuse through nanopores.

First-Principles Simulations of Liquid Water Using a Dielectric-Dependent Hybrid Functional

We carried out first-principles simulations of liquid water under ambient conditions using a dielectric-dependent hybrid functional, where the fraction of exact exchange is set equal to the inverse of the high-frequency dielectric constant of the liquid. We found excellent agreement with experiment for the oxygen-oxygen partial correlation function at the experimental equilibrium density and 311 ± 3 K. Other structural and dynamical properties, such as the diffusion coefficient, molecular dipole moments, and vibrational spectra, are also in good agreement with experiment. Our results, together with previous findings on electronic properties of the liquid with the same functional, show that the dielectric-dependent hybrid functional accurately describes both the structural and electronic properties of liquid water.

Molecular investigation of the wettability of rough surfaces using molecular dynamics simulation

In the present study, a computational investigation on the effect of surface roughness on the wettability behavior of water nanodroplets has been performed via molecular dynamics simulation.

Quantum mechanical force fields for condensed phase molecular simulations

Molecular simulations are powerful tools for providing atomic-level details into complex chemical and physical processes that occur in the condensed phase. For strongly interacting systems where quantum many-body effects are known to play an important role, density-functional methods are often used to provide the model with the potential energy used to drive dynamics. These methods, however, suffer from two major drawbacks. First, they are often too computationally intensive to practically apply to large systems over long time scales, limiting their scope of application. Second, there remain challenges for these models to obtain the necessary level of accuracy for weak non-bonded interactions to obtain quantitative accuracy for a wide range of condensed phase properties. Quantum mechanical force fields (QMFFs) provide a potential solution to both of these limitations. In this review, we address recent advances in the development of QMFFs for condensed phase simulations. In particular, we examine the development of QMFF models using both approximate and ab initio density-functional models, the treatment of short-ranged non-bonded and long-ranged electrostatic interactions, and stability issues in molecular dynamics calculations. Example calculations are provided for crystalline systems, liquid water, and ionic liquids. We conclude with a perspective for emerging challenges and future research directions.

Structure and dynamics of liquid CS2: Going from ambient to elevated pressure conditions

Molecular dynamics simulation studies were performed to investigate the structural and dynamic properties of liquid carbon disulfide (CS₂) from ambient to elevated pressure conditions. The results obtained have revealed structural changes at high pressures, which are related to the more dense packing of the molecules inside the first solvation shell. The calculated neutron and X-ray structure factors have been compared with available experimental diffraction data, also revealing the pressure effects on the short-range structure of the liquid. The pressure effects on the translational, reorientational, and residence dynamics are very strong, revealing a significant slowing down when going from ambient pressure to 1.2 GPa. The translational dynamics of the linear CS₂ molecules have been found to be more anisotropic at elevated pressures, where cage effects and librational motions are reflected on the shape of the calculated time correlation functions and their corresponding spectral densities.

A quantum-mechanics molecular-mechanics scheme for extended systems

We introduce and discuss a hybrid quantum-mechanics molecular-mechanics (QM-MM) approach for Car-Parrinello DFT simulations with pseudopotentials and planewaves basis, designed for the treatment of periodic systems. In this implementation the MM atoms are considered as additional QM ions having fractional charges of either sign, which provides conceptual and computational simplicity by exploiting the machinery already existing in planewave codes to deal with electrostatics in periodic boundary conditions. With this strategy, both the QM and MM regions are contained in the same supercell, which determines the periodicity for the whole system. Thus, while this method is not meant to compete with non-periodic QM-MM schemes able to handle extremely large but finite MM regions, it is shown that for periodic systems of a few hundred atoms, our approach provides substantial savings in computational times by treating classically a fraction of the particles. The performance and accuracy of the method is assessed through the study of energetic, structural, and dynamical aspects of the water dimer and of the aqueous bulk phase. Finally, the QM-MM scheme is applied to the computation of the vibrational spectra of water layers adsorbed at the TiO2 anatase (1 0 1) solid-liquid interface. This investigation suggests that the inclusion of a second monolayer of H2O molecules is sufficient to induce on the first adsorbed layer, a vibrational dynamics similar to that taking place in the presence of an aqueous environment. The present QM-MM scheme appears as a very interesting tool to efficiently perform molecular dynamics simulations of complex condensed matter systems, from solutions to nanoconfined fluids to different kind of interfaces.

A flow cell for measuring X-ray Compton scattering of liquid at temperatures up to 623 K and pressures up to 20 MPa

A flow-type cell was developed for measuring Compton scattering spectra of heat-sensitive aqueous solution. Compton scattering spectra of water and ethanol were measured in the region from ambient conditions to 623 K and 20 MPa. Compton profiles derived from measurement with the flow-type cell were comparable with those in the literature. Results obtained from the flow-type cell showed that delocalization of electronic charge density of water and ethanol at high temperatures occurred. Delocalization of the electronic charge density of ethanol was greater than that of water at high temperature, which is consistent with the prior works that use proton NMR chemical shifts to describe hydrogen bonding.