Hydrogen bonds in aqueous electrolyte solutions: statistics and dynamics based on both geometric and energetic criteria

Enhanced fluidity of water in superhydrophobic nanotubes: estimating viscosity using jump-corrected confined Stokes-Einstein approach

Accurately predicting the viscosity of water confined within nanotubes is vital for various technological applications. Traditional methods have failed in this regard, necessitating a novel approach. We introduced the jump-corrected confined Stokes-Einstein (JCSE) method and now employ the same to estimate the viscosity and diffusion in superhydrophobic nanotubes. Our study covers a temperature range of 230-300 K and considers three nanotube diameters. Results show that water inside superhydrophobic nanotubes exhibits a significantly lower viscosity and higher diffusion than those inside hydrophobic nanotubes. Narrower nanotubes and lower temperatures accentuate these effects. Furthermore, water inside superhydrophobic nanotubes display a lower viscosity than bulk water, with the difference increasing at lower temperatures. This reduction is attributed to weaker water-water interactions caused by a lower water density in the interfacial region. These findings highlight the importance of interfacial water density and its influence on nanotube viscosity, shedding light on nanoscale fluid dynamics and opening avenues for diverse applications.

Interaction of Nitrite Ions with Hydrated Portlandite Surfaces: Atomistic Computer Simulation Study

The nitrite admixtures in cement and concrete are used as corrosion inhibitors for steel reinforcement and also as anti-freezing agents. The characterization of the protective properties should account for the decrease in the concentration of free NO₂^- ions in the pores of cement concretes due to their adsorption. Here we applied the classical molecular dynamics computer simulation approach to quantitatively study the molecular scale mechanisms of nitrite adsorption from NaNO₂ aqueous solution on a portlandite surface. We used a new parameterization to model the hydrated NO₂^- ions in combination with the recently upgraded ClayFF force field (ClayFF-MOH) for the structure of portlandite. The new NO₂^- parameterization makes it possible to reproduce the properties of hydrated NO₂^- ions in good agreement with experimental data. In addition, the ClayFF-MOH model improves the description of the portlandite structure by explicitly taking into account the bending of Ca-O-H angles in the crystal and on its surface. The simulations showed that despite the formation of a well-structured water layer on the portlandite (001) crystal surface, NO₂^- ions can be strongly adsorbed. The nitrite adsorption is primarily due to the formation of hydrogen bonds between the structural hydroxyls on the portlandite surface and both the nitrogen and oxygen atoms of the NO₂^- ions. Due to that, the ions do not form surface adsorption complexes with a single well-defined structure but can assume various local coordinations. However, in all cases, the adsorbed ions did not show significant surface diffusional mobility. Moreover, we demonstrated that the nitrite ions can be adsorbed both near the previously-adsorbed hydrated Na⁺ ions as surface ion pairs, but also separately from the cations.

Vibrational signature of hydrated protons confined in MXene interlayers

The hydration structure of protons has been studied for decades in bulk water and protonated clusters due to its importance but has remained elusive in planar confined environments. Two-dimensional (2D) transition metal carbides known as MXenes show extreme capacitance in protic electrolytes, which has attracted attention in the energy storage field. We report here that discrete vibrational modes related to protons intercalated in the 2D slits between Ti₃C₂T_x MXene layers can be detected using operando infrared spectroscopy. The origin of these modes, not observed for protons in bulk water, is attributed to protons with reduced coordination number in confinement based on Density Functional Theory calculations. This study therefore demonstrates a useful tool for the characterization of chemical species under 2D confinement.

Ionic Dynamics and Vibrational Spectral Diffusion of a Protic Alkylammonium Ionic Salt through Intrinsic Cationic N-H Vibrational Probe from FPMD Simulations

We employed density functional theory (DFT)-based molecular dynamics simulations to explore the structure, dynamics, and spectral properties of the protic ionic entity trimethylammonium chloride (TMACl). Structural investigations include calculating the site-site radial distribution functions (RDFs), the distribution of constituent cations and anions in three-dimensional space, and combined distribution functions of the hydrogen-bonded pair RDF versus angle, revealing the structural characteristics of the ionic solvation and the intermolecular interactions within ions. Further, we determined the instantaneous vibrational stretching frequencies of the intrinsic N-H stretch probe modes by applying the time-series wavelet method. The associated ionic dynamics within the protic ionic compound were investigated by calculating the time-evolution of the fluctuating frequencies and the frequency-time correlation functions (FFCFs). The time scale related to the local structural relaxation process and the average hydrogen bond lifetime, ion cage dynamics, and mean squared displacement were investigated. The faster decay component of the FFCFs, depicting the intermolecular motion of intact hydrogen bonds in TMACl, is 0.07 ps for the Perdew-Burke-Ernzerhof (PBE)-based simulation and 0.06 ps for the PBE-D2 representation. The slower time scale of the longer picosecond decay time component of PBE and PBE-D2 representations are 3.13 and 2.87 ps, respectively. These picosecond time scales represent more significant fluctuations of the hydrogen-bonding partners in the ionic entity and hydrogen-bond jump events accompanied by large angular jumps. The longest picosecond time scales represent structural relaxation, including large angular jumps and ion-pair dynamics. Also, ion cage lifetimes correlate with the slowest time scale of the associated dynamics of vibrational spectral diffusion despite the type of DFT functional. This study benchmarks DFT treatments of the exchange-correlation functional with and without the van der Waals (vdW) dispersion correction scheme. The inclusion of vdW interactions to the PBE functional represents a less structured state of the ionic entity and faster dynamics of the molecular motions relative to the one predicted by the PBE system. All the results illustrate the necessity of accurately describing the Coulomb interactions, vdW dispersive interactive forces, and localized hydrogen bonds required to sustain the energetic balance in this ionic salt.

Hydration in aqueous NaCl

Atomistic details about the hydration of ions in aqueous solutions are still debated due to the disordered and statistical nature of the hydration process. However, many processes from biology, physical chemistry to materials sciences rely on the complex interplay between solute and solvent. Oxygen K-edge X-ray excitation spectra provide a sensitive probe of the local atomic and electronic surrounding of the excited sites. We used ab initio molecular dynamics simulations together with extensive spectrum calculations to relate the features found in experimental oxygen K-edge spectra of a concentration series of aqueous NaCl with the induced structural changes upon solvation of the salt and distill the spectral fingerprints of the first hydration shells around the Na⁺- and Cl^--ions. By this combined experimental and theoretical approach, we find the strongest spectral changes to indeed result from the first hydration shells of both ions and relate the observed shift of spectral weight from the post- to the main-edge to the origin of the post-edge as a shape resonance.

A computational approach on the stereoselective binding of peptides from aqueous medium with endo-functionalized molecular tubes

The need to obtain enantiomerically pure isomers of amino acids and peptides is often realized in the field of biology and in the pharmaceutical industry. Research is underway to devise simple methods for the chiral resolution of amino acids from their racemic mixtures. Inspired by this objective, in our present work, we have computationally shown the possibility of chiral separation of the enantiomeric pairs of two model peptides, namely, (D,L)-aspargine and (D,L)-phenylalanine, in the presence of water. For this purpose, we have used two synthetic supramolecular receptors named host-1a and host-1b, respectively. Molecular dynamics simulations and quantum chemical methods are employed to analyze the structural features and the energy aspects involved in the separation process. The information obtained at the molecular level helps us gain better insights into the key interactions that operate to produce such enantioselectivity. We have also investigated the dynamics and changes in the water structure in the vicinity of the host molecules, both in the presence and absence of the model peptides. The D- and L-isomers of the same peptide undergo complexation with a particular host molecule registering a difference of more than 1.5 kcal mol^-1 (obtained from PMF and MM-PBSA analyses) in their respective energies. This indicates that the chiral separation of the peptides with the help of these endo-fuctionalized molecular tube receptors may be energetically feasible. The connection between the peptide stereochemistry and its interaction with the endo-functionalized hosts would be instrumental in designing novel segregation techniques that can be further extended to separate larger model peptides or proteins.

Dynamics of Ionic Liquid through Intrinsic Vibrational Probes Using the Dispersion-Corrected DFT Functionals

First principles molecular dynamics simulations have been utilized to study the spectral properties of the protic ionic liquid, methylammonium formate (MAF). All simulations were performed using density functional theory (DFT) and various van der Waals-corrected exchange-correlation functionals. We calculated the vibrational stretch frequency distributions, determined the time-frequency correlations of the intrinsic vibrational probes, the N-H and C-O modes in MAF, and the frequency-structure correlations. We also estimated the average hydrogen-bond lifetimes and orientation dynamics to capture the ultrafast spectral response. The spectroscopic signature of the N-H stretching vibrations using the Becke-Lee-Yang-Parr (BLYP) and Perdew-Burke-Ernzerhof (PBE) functionals displays a spectral shift in the lower frequency side, suggesting stronger hydrogen-bonding interactions represented by the gradient approximation functionals than the van der Waals (vdW)-corrected simulations. The carboxylate frequency profiles with the dispersion-corrected representations are almost similar without a significant difference in the normalized distributions. Besides, the COO stretching frequencies at the peak maxima positions of the PBE functionals exhibit a lesser deviation from the experimental data. Spectral diffusion dynamics of the intrinsic vibrational probes on the cationic and anionic sites of the ionic liquid proceed through a short time relaxation of the intact hydrogen bonds followed by an intermediate time constant and a longer time decay indicating the switchover of hydrogen bonds. Dispersion-corrected atom-centered one-electron potential (DCACP) correction added to the BLYP system slows down the picosecond time scales of frequency correlation and the time constants of rotational motion, lengthening the overall system dynamics. The observed trends in the time-dependent decays of frequency fluctuations and the orientation autocorrelation functions correlate with the structural interactions in liquid MAF and hydrogen-bond dynamics. In this study, we examine the predictions made by different density functional treatments comparing the results of the uncorrected BLYP and PBE representations with the semiempirical vdW methods of Grimme and matching our calculated data with the experimental observations.

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.

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.

Role of hydrogen bond capacity of solvents in reactions of amines with CO2: A computational study

Various computational methods were employed to investigate the zwitterion formation, a critical step for the reaction of monoethanolamine with CO₂, in five solvents (water, monoethanolamine, propylamine, methanol and chloroform) to probe the effect of hydrogen bond capacity of solvents on the reaction of amine with CO₂ occurring in the amine-based CO₂ capture process. The results indicate that the zwitterion can be formed in all considered solvents except chloroform. For two pairs of solvents (methanol and monoethanolamine, propylamine and chloroform) with similar dielectric constant but different hydrogen bond capacity, the solvents with higher hydrogen bond capacity (monoethanolamine and propylamine) facilitate the zwitterion formation. More importantly, kinetics parameters such as activation free energy for the zwitterion formation are more relevant to the hydrogen bond capacity than to dielectric constant of the considered solvents, clarifying the hydrogen bond capacity could be more important than dielectric constant in determining the kinetics of monoethanolamine with CO₂.

Insight into Cellulose Dissolution with the Tetrabutylphosphonium Chloride-Water Mixture using Molecular Dynamics Simulations

All-atom molecular dynamics simulations are utilized to determine the properties and mechanisms of cellulose dissolution using the ionic liquid tetrabutylphosphonium chloride (TBPCl)–water mixture, from 63.1 to 100 mol % water. The hydrogen bonding between small and large cellulose bundles with 18 and 88 strands, respectively, is compared for all concentrations. The Cl, TBP, and water enable cellulose dissolution by working together to form a cooperative mechanism capable of separating the cellulose strands from the bundle. The chloride anions initiate the cellulose breakup, and water assists in delaying the cellulose strand reformation; the TBP cation then more permanently separates the cellulose strands from the bundle. The chloride anion provides a net negative pairwise energy, offsetting the net positive pairwise energy of the peeling cellulose strand. The TBP–peeling cellulose strand has a uniquely favorable and potentially net negative pairwise energy contribution in the TBPCl–water solution, which may partially explain why it is capable of dissolving cellulose at moderate temperatures and high water concentrations. The cellulose dissolution declines rapidly with increasing water concentration as hydrogen bond lifetimes of the chloride–cellulose hydroxyl hydrogens fall below the cellulose’s largest intra-strand hydrogen bonding lifetime.

Molecular dynamics simulations of plutonium binding and its decorporation from the binding-cleft of human serum transferrin

The possibility of plutonium (Pu) intake by radiation workers can not be ruled out. Transportation of Pu(IV) to various organs/cells is mainly carried through iron-carrying protein, serum transferrin (sTf), by receptor-mediated endocytosis. Understanding the Pu-sTf interaction is a primary step toward future design of its decorporating agents. We report MD simulations of Pu(IV) binding with sTf and look out for its decorporation at extracellular pH using suitable ligands. MD simulations were carried out in polarizable water environment at different protonation states of the protein. Results unravel the binding motif of Pu(IV): (1) sTf binds the ion in closed conformation at extracellular serum pH with carbonate as synergistic anions, (2) change in protonation state of dilysine (K206 and K296)-trigger and that of the carbonate ion at acidic endosomal pH is found to cause conformational changes of protein, conducive for the heavy ion to be released, although; (3) strong electrostatic interaction between D63 in the binding-cleft and Pu(IV) is found not to ever set free the ion. In an endeavour to decorporate Pu(IV), fragmented molecular form of hydroxypyridinone (HOPO) and catechol (CAM)-based ligands are docked at the binding site (BS) of the protein and metadynamics simulations are conducted. Pu(IV) binding at BS is found to be so strong that it was not detached from BS with the docked HOPO. However, for the identical set of simulation parameters, CAM is found to facilitate dislodging the heavy ion from the protein's binding influence. Differential behaviour of the two chelators is further explored. Fragmented molecular form of hydroxy-pyridinone (HOPO) and catecholamide (CAM) ligands were docked at the binding-site (BS) of human serum transferrin (sTf) to explore their feasibility as plausible Pu(IV) decorporating agents by employing metadynamics method. CAM was found to dislodge Pu from the sTf BS, while HOPO could not.

Thermodynamic and Transport Properties of Tetrabutylphosphonium Hydroxide and Tetrabutylphosphonium Chloride-Water Mixtures via Molecular Dynamics Simulation

Thermodynamic, structural, and transport properties of tetrabutylphosphonium hydroxide (TBPH) and tetrabutylphosphonium chloride (TBPCl)–water mixtures have been investigated using all-atom molecular dynamics simulations in response to recent experimental work showing the TBPH–water mixtures capability as a cellulose solvent. Multiple transitional states exist for the water—ionic liquid (IL) mixture between 70 and 100 mol% water, which corresponds to a significant increase in water hydrogen bonds. The key transitional region, from 85 to 92.5 mol% water, which coincides with the mixture’s maximum cellulose solubility, reveals small and distinct water veins with cage structures formed by the TBP⁺ ions, while the hydroxide and chloride ions have moved away from the P atom of TBP⁺ and are strongly hydrogen bonded to the water. The maximum cellulose solubility of the TBPH–water solution at approximately 91.1 mol% water, appears correlated with the destruction of the TBP’s interlocking structure in the simulations, allowing the formation of water veins and channeling structures throughout the system, as well as changing from a subdiffusive to a near-normal diffusive regime, increasing the probability of the IL’s interaction with the cellulose polymer. A comparison is made between the solution properties of TBPH and TBPCl with those of alkylimidazolium-based ILs, for which water appears to act as anti-solvent rather than a co-solvent.

Tuning Water Networks via Ionic Liquid/Water Mixtures

Water in nanoconfinement is ubiquitous in biological systems and membrane materials, with altered properties that significantly influence the surrounding system. In this work, we show how ionic liquid (IL)/water mixtures can be tuned to create water environments that resemble nanoconfined systems. We utilize molecular dynamics simulations employing ab initio force fields to extensively characterize the water structure within five different IL/water mixtures: [BMIM+][BF4−], [BMIM+][PF6−], [BMIM+][OTf−], [BMIM+][NO3−] and [BMIM+][TFSI−] ILs at varying water fraction. We characterize water clustering, hydrogen bonding, water orientation, pairwise correlation functions and percolation networks as a function of water content and IL type. The nature of the water nanostructure is significantly tuned by changing the hydrophobicity of the IL and sensitively depends on water content. In hydrophobic ILs such as [BMIM+][PF6−], significant water clustering leads to dynamic formation of water pockets that can appear similar to those formed within reverse micelles. Furthermore, rotational relaxation times of water molecules in supersaturated hydrophobic IL/water mixtures indicate the close-connection with nanoconfined systems, as they are quantitatively similar to water relaxation in previously characterized lyotropic liquid crystals. We expect that this physical insight will lead to better design principles for incorporation of ILs into membrane materials to tune water nanostructure.

The role of urea in the solubility of cellulose in aqueous quaternary ammonium hydroxide

We examine the role of water and urea in cellulose solubility in tetrabutylammonium hydroxide (TBAH). Molecular dynamics simulations were performed for several different solvent compositions with a fixed cellulose fraction. For each composition, two simulations were carried out with cellulose fixed in each of the crystalline and the dissolved states. From the enthalpy and the entropy of the two states, the difference in Gibbs free energy (ΔG) and hence the spontaneity is determined. A comparison with solubility experiments showed a strong correlation between the calculated ΔG and the experimental measurements. A breakdown of the enthalpic and entropic contributions reveals the roles of water and urea in solubility. At high water concentration, a drop in solubility is attributed to both increased enthalpy and decreased entropy of dissolution. Water displaces strong IL–cellulose interactions for weaker water–cellulose interactions, resulting in an overall enthalpy increase. This is accompanied by a strong decrease in entropy, which is primarily attributed to both water and the entropy of mixing. Adding urea to TBAH(aq) increases solubility by an addition to the mixing term and by reducing losses in solvent entropy upon dissolution. In the absence of urea, the flexible [TBA]⁺ ions lose substantial degrees of freedom when they interact with cellulose. When urea is present, it partially replaces [TBA]⁺ and to a lesser extent OH⁻ near cellulose, losing less entropy because of its rigid structure. This suggests that one way to boost the dissolving power of an ionic liquid is to limit the number of degrees of freedom from the outset.

We examine the role of water and urea in cellulose solubility in tetrabutylammonium hydroxide (TBAH).

An ab initio molecular dynamics study of benzene in water at supercritical conditions: Structure, dynamics, and polarity of hydration shell water and the solute

Anisotropic structure and dynamics of the hydration shell of a benzene solute in supercritical water are investigated by means of ab initio molecular dynamics simulations. The polarity and structural distortion of the benzene solute in supercritical water are also investigated in this study. Calculations are done at 673 K for three different densities of the solvent. The simulations are carried out using the Becke-Lee-Yang-Parr (BLYP) and also the Becke-Lee-Yang-Parr functional including dispersion corrections of Grimme (BYLP-D). The structural anisotropy is found to exist even at supercritical conditions as elucidated by the radial distribution functions of different conical regions and also by angular and spatial distribution functions. The benzene-water πH-bond and also the water-water hydrogen bonds are found to exist even at the supercritical temperature of 673 K. However, the numbers of these hydrogen bonds are reduced substantially with a decrease in water density. The water molecules in the axial region of benzene are found to be preferably oriented with one OH vector pointing toward the benzene ring, whereas the water molecules located in the equatorial region are found to orient their dipoles mostly parallel to the ring plane. The orientational distributions, however, are found to be rather broad at the supercritical temperature due to thermal fluctuations. Although the water molecules have faster dynamics at these supercritical conditions, a slight difference is observed in the dynamics of the solvation shell and bulk molecules. The conformational flexibility of the ring is found to be enhanced which causes an increase in polarity of the benzene solute in water under supercritical conditions.

Dynamical crossover of water confined within the amphiphilic nanocores of aggregated amyloid β peptides

It is believed that the self-assembly of amyloid beta (Aβ) peptides in the brain is the cause of Alzheimer's disease. Atomistic molecular dynamics simulations of aqueous solutions of Aβ protofilaments of different sizes at room temperature have been carried out to explore the dynamic properties of water confined within the core and at the exterior surface of the protofilaments. Attempts have been made to understand how the non-uniform distortion of the protofilaments associated with their structural crossover influences the diffusivity and the hydrogen bonding environment of the confined water molecules. In contrast to the homogeneous solvent dynamical environment at the exterior surface, the calculations revealed heterogeneously restricted motions of water confined within the distorted cores of the protofilaments. Importantly, it is demonstrated that the structural crossover of the aggregates observed for the decamer is associated with a dynamical transition of water confined within its core. A direct one-to-one correlation between the heterogeneously restricted core water motions and the kinetics of the breaking and formation of hydrogen bonds quantitatively demonstrated that a modified hydrogen bond arrangement within the cores of higher order Aβ protofilaments is the origin behind the crossover in core water mobility. A fraction of the water molecules forming short-lived water-water hydrogen bonds within the core of the crossover protofilament decamer are believed to diffuse away easily from the core and thus play a crucial role in further growth of the protofilament by facilitating the binding of new peptide monomers.

Nuclear quantum effects on the vibrational dynamics of liquid water

Based on quantum-mechanical path-integral molecular dynamics simulations, the impact of nuclear quantum effects on the vibrational and hydrogen bond dynamics in liquid water is investigated. The instantaneous fluctuations in the frequencies of the O-H stretch modes are calculated using the wavelet method of time-series analysis, while the time scales of the vibrational spectral diffusion are determined from frequency-time correlation functions, joint probability distributions, and the slope of three-pulse photon echo. We find that the inclusion of nuclear quantum effects leads not only to a redshift of the vibrational frequency distribution by around 130 cm^-1 but also to an acceleration of the vibrational dynamics by as much as 30%. In addition, quantum fluctuations also entail a significantly faster decay of correlation in the initial diffusive regime, which is in agreement with recent vibrational echo experiments.

On the hydrogen-bond network and the non-Arrhenius transport properties of water

We study the structural and dynamic transformations of SPC/E water with temperature, through molecular dynamics (MD), and discuss the non-Arrhenius behavior of the transport properties and orientational dynamics, and the magnitude of the breakdown of the Stokes-Einstein (SE) and the Stokes-Einstein-Debye (SED) relations, in the light of these transformations. Our results show that deviations from Arrhenius behavior of the self-diffusion at low temperatures cannot be exclusively explained by the reduction of water defects (interstitial waters) and the increase of the local tetrahedrality, thus, suggesting the importance of the slowdown of collective rearrangements. Interestingly we find that at high temperatures (T  ⩾  340 K) water defects lead to a slight increase of the tetrahedrality and a decrease of the self-diffusion, opposite to water at low temperatures. The relative magnitude of the breakdown of the SE and the SED relations is found to be in accord with recent experiments (Dehaoui et al 2015 Proc. Natl Acad. Sci. USA 112 12020) resolving the discrepancy with previous MD results. Further, we show that SPC/E hydrogen-bond (HB) lifetimes deviate from Arrhenious behaviour at low temperatures in contrast with some previous MD studies. This deviation is nevertheless much smaller than that observed for the orientational dynamics and the transport properties of water, consistent with the relaxation times measured by several experimental methods. The HB acceptor exchange dynamics defined here by the acceptor switch and reform (librational dynamics) frequencies exhibit similar Arrhenius deviations, thus explaining to some extent the non-Arrhenius behavior of the transport properties and of the orientational dynamics of water. Our results also show that the fraction of HB switches through a bifurcated pathway follow a power law with the temperature decrease. Thus, at low temperatures HB acceptor switches are less frequent but occur on a faster time scale consistent with the temperature dependence of the ratio of the rotational relaxation times for the different Legendre polynomial ranks.

Probing the dynamics of N-methylacetamide in methanol via ab initio molecular dynamics

Two-dimensional infrared (2D IR) spectroscopy of amide 1 vibrational bands provides a valuable probe of proteins as well as molecules such as N-methylacetamide (NMA), which present peptide-like H-bonding possibilities to a solvent.

Anion Dependent Dynamics and Water Solubility Explained by Hydrogen Bonding Interactions in Mixtures of Water and Aprotic Heterocyclic Anion Ionic Liquids

Molecular dynamics simulations were used to compare water solubilities and the effects of water on the structure and dynamics of ionic liquids (ILs) composed of phosphonium cations paired with azolide and phenolate anions. The addition of water decreases ordering of the ions compared to the dry ILs with the exception of anion-anion ordering in the phenolate IL. The result is that the dynamics of the azolide ionic liquids increase significantly upon addition of water, whereas the phenolate IL dynamics show little change. The relative water solubilities were compared through calculation of Henry's law constants. Water is much more soluble in the phenolate IL due to strong hydrogen bonding interactions between water and the phenolate oxygen atom. Anions can therefore be selected to control IL-water hydrogen bonding for optimal performance in applications such as CO₂ separation.

Anisotropic structure and dynamics of the solvation shell of a benzene solute in liquid water from ab initio molecular dynamics simulations

The anisotropic structure and dynamics of the hydration shell of a benzene solute in liquid water have been investigated by means of ab initio molecular dynamics simulations using the BLYP (Becke–Lee–Yang–Parr) and dispersion corrected BLYP-D functionals.

Mechanisms of hydrogen bond formation between ionic liquids and cellulose and the influence of water content

We explore the complex network of transitions occurring between different hydrogen bonding states within ionic liquids and cellulose.

SINGLE-PARTICLE AND PAIR DYNAMICAL PROPERTIES OF ACETONE–METHANOL MIXTURES CONTAINING CHARGED AND NEUTRAL SOLUTES: A MOLECULAR DYNAMICS STUDY

The dynamical properties of acetone–methanol mixtures containing either an ionic or a neutral hydrophobic solute are investigated by means of a series of molecular dynamics simulations. The primary goal has been to study how the solute and solvent dynamical properties change with variation of composition of the mixture ranging from pure acetone to pure methanol. The variations of structure and energetics of the mixture with composition are also calculated. The diffusion coefficients of both ionic and neutral solutes are found to show nonlinear variation with composition of the mixture, although the extent of nonlinearity in the diffusion of the neutral solute is much weaker. Calculations of appropriate solute-solvent distribution functions reveal the extent and nature of selective solvation of these solute species which play a role in determining the nonideal dynamical characteristics of these solutes. The free energies of solvation of the ionic solutes are also calculated and the results are discussed in the context of their dynamical behavior. The hydrogen bond statistics and dynamics of these mixtures are also calculated over their entire composition range. The energies and lifetimes of hydrogen bonds between an acetone and a methanol molecule or between two methanol molecules are found to increase with increase of acetone mole fraction of the mixture. Residence times of methanol molecules in solvation shells of acetone and methanol are also found to follow the same trend as relaxation times. However, these pair dynamical properties show essentially linear dependence on composition, thus behave almost ideally with respect to changes in composition of the mixture.