Hydrogen-bond dynamics of water in presence of an amphiphile, tetramethylurea: signature of confinement-induced effects

Temperature-dependent dielectric relaxation measurements of (acetamide + K/Na SCN) deep eutectic solvents: Decoding the impact of cation identity via computer simulations

The impact of successive replacement of K+ by Na+ on the megahertz-gigahertz polarization response of 0.25[fKSCN + (1 - f)NaSCN] + 0.75CH3CONH2 deep eutectic solvents (DESs) was explored via temperature-dependent (303 ≤ T/K ≤ 343) dielectric relaxation (DR) measurements and computer simulations. Both the DR measurements (0.2 ≤ ν/GHz ≤ 50) and the simulations revealed multi-Debye relaxations accompanied by a decrease in the solution static dielectric constant (ɛs) upon the replacement of K+ by Na+. Accurate measurements of the DR response of DESs below 100 MHz were limited by the well-known one-over-frequency divergence for conducting solutions. This problem was tackled in simulations by removing the zero frequency contributions arising from the ion current to the total simulated DR response. The temperature-dependent measurements revealed a much stronger viscosity decoupling of DR times for Na+-containing DES than for the corresponding K+ system. The differential scanning calorimetry measurements indicated a higher glass transition temperature for Na+-DES (∼220 K) than K+-DES (∼200 K), implying more fragility and cooperativity for the former (Na+-DES) than the latter. The computer simulations revealed a gradual decrease in the average number of H bonds (⟨nHB⟩) per acetamide molecule and increased frustrations in the average orientational order upon the replacement of K+ by Na+. Both the measured and simulated ɛs values were found to decrease linearly with ⟨nHB⟩. Decompositions of the simulated DR spectra revealed that the cation-dependent cross interaction (dipole-ion) term contributes negligibly to ɛs and appears in the terahertz regime. Finally, the simulated collective single-particle reorientational relaxations and the structural H-bond fluctuation dynamics revealed the microscopic origin of the cation identity dependence shown by the measured DR relaxation times.

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.

Investigating the Vapor-Phase Adsorption of Aroma Molecules on the Water-Vapor Interface using Molecular Dynamics Simulations

Surfactants are amphiphilic additives primarily used to reduce the surface tension of water and manipulate its wettability on various surfaces. Recent reports suggest that volatile surfactants, such as aroma molecules, diffuse more quickly to the interface from the vapor-phase than conventional surfactants typically used in the aqueous phase. The ability to adsorb from the vapor phase, in addition to their use as cosurfactants, expands the potential applications of volatile surfactants, particularly in situations where adding surfactants from the liquid phase is difficult. Here, we present a molecular level understanding of the adsorption kinetics of linalool, a common aroma molecule, on the water interface using molecular dynamics simulations. We note that the value of surface tension while adsorption from vapor and liquid phases is dependent only on the surface coverage. A minimum surface tension of 32 ± 1.8 mN/m is obtained in both cases at a maximum surface coverage of 4.88 μmol/m2 at 300 K. We observe the extent of decrease of the H-bonds between linalool-water and linalool-linalool molecules at various surface coverages to explain the mechanism of surface tension reduction. We solve Gibb's adsorption equation to establish a correlation between the surface coverage of linalool and the corresponding bulk concentration in experiments. We investigate the free energy profile of linalool's adsorption behavior at different surface coverages and temperatures. Our report suggests that linalool adsorption onto the water interface is an enthalpy-driven process primarily dependent on the strength of the interaction between the hydroxyl group of linalool and water molecules. These insights are crucial for selecting a suitable aroma molecule for various applications that target the vapor-phase adsorption mechanism.

Non-monotonic composition dependence of the breakdown of Stokes-Einstein relation for water in aqueous solutions of ethanol and 1-propanol: explanation using translational jump-diffusion approach

A series of experimental and simulation studies examined the validity of the Stokes-Einstein relationship (SER) of water in binary water/alcohol mixtures of different mixture compositions. These studies revealed a strong non-monotonic composition dependence of the SER with maxima at the specific alcohol mole fraction where the non-idealities of the thermodynamic and transport properties are observed. The translational jump-diffusion (TJD) approach elucidated the breakdown of the SER in pure supercooled water as caused by the jump translation of molecules. The breakdown of SER in the supercooled water/methanol binary mixture was successfully explained using the same TJD approach. To further generalize the picture, here we focus on the non-monotonic composition dependence of SER breakdown of water in two water/alcohol mixtures (water/ethanol and water/propanol) for a broad temperature range. In agreement with previous studies, maximum breakdown of SER is observed for the mixture with alcohol mole fraction x = 0.2. Diffusion of the water molecules at the maximum SER breakdown point is largely contributed by jump-diffusion. The residual-diffusion, obtained by subtracting the jump-diffusion from the total diffusion, approximately follows the SER for different compositions and temperatures. We also performed hydrogen (H-)bond dynamics and observed that the contribution of jump-diffusion is proportional to the total free energy of activation of breaking all H-bonds that exist around a molecule. This study, therefore, suggests that the more a molecule is trapped by H-bonding, the more likely it is to diffuse through the jump-diffusion mechanism, eventually leading to an increasing degree of SER breakdown.

Translational Jump-Diffusion of Hydroxide Ion in Anion Exchange Membrane: Deciphering the Nature of Vehicular Diffusion

Earlier, ab initio and reactive force-field-based molecular dynamics (MD) simulation studies suggested an overwhelming contribution of the vehicular diffusion in the total diffusion of hydroxide ions rather than structural diffusion. But does the vehicular diffusion occur via small-step displacement? This question is important to have an understanding of the real characteristics of vehicular diffusion. To answer this question, we perform a classical molecular dynamics simulation of a system containing a hydroxide ion exchange membrane polymer and hydroxide ion at different hydration levels and temperatures using the same molecular force field (Dubey, V. Chem. Phys. Lett. 2020, 755, 137802), which successfully captured the microscopic structure and dynamics of the system. We use the translational jump-diffusion approach, used previously in supercooled water for understanding the origin of breakdown of the Stokes-Einstein relation, to calculate the jump-diffusion coefficient of hydroxide ion and water in the anion exchange membrane. We have seen a significant role of hydration level and temperature in the mechanism of vehicular diffusion. In overhydrated membrane, both hydroxide ions and water molecules diffuse via both small- and large-step displacement. With decreasing hydration level and temperature, the diffusion is increasingly governed by the jump-diffusion mechanism. The larger contribution of jump-diffusion comes from the stronger caging of the diffusing species by the solvent at lower hydration levels and temperature. These results, therefore, suggest that the hydration level and temperature of the hydroxide ion exchange membrane determine the detailed mechanism of the vehicular diffusion of hydroxide ion, especially whether the diffusion follows hydrodynamics or not.

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.

Breakdown of the Stokes-Einstein relation in supercooled water: the jump-diffusion perspective

Although water is the most ubiquitous liquid it shows many thermodynamic and dynamic anomalies. Some of the anomalies further intensify in the supercooled regime. While many experimental and theoretical studies have focused on the thermodynamic anomalies of supercooled water, fewer studies explored the dynamical anomalies very extensively. This is due to the intricacy of the experimental measurement of the dynamical properties of supercooled water. Violation of the Stokes-Einstein relation (SER), an important relation connecting the diffusion of particles with the viscosity of the medium, is one of the major dynamical anomalies. In absence of experimentally measured viscosity, researchers used to check the validity of SER indirectly using average translational relaxation time or α-relaxation time. Very recently, the viscosity of supercooled water was accurately measured at a wide range of temperatures and pressures. This allowed direct verification of the SER at different temperature-pressure thermodynamic state points. An increasing breakdown of the SER was observed with decreasing temperature. Increasing pressure reduces the extent of breakdown. Although some well-known theories explained the above breakdown, a detailed molecular mechanism was still elusive. Recently, a translational jump-diffusion (TJD) approach has been able to quantitatively explain the breakdown of the SER in pure supercooled water and an aqueous solution of methanol. The objective of this article is to present a detailed and state-of-the-art analysis of the past and present works on the breakdown of SER in supercooled water with a specific focus on the new TJD approach for explaining the breakdown of the SER.

Breakdown of the Stokes-Einstein Relation in Supercooled Water/Methanol Binary Mixtures: Explanation Using the Translational Jump-Diffusion Approach

A recent experiment has directly checked the validity of the Stokes-Einstein (SE) relation for pure water, pure methanol, and their binary mixtures of three different compositions at different temperatures. The effect of composition on the nature of breakdown of the SE relation is interesting. While in the majority of the systems, an increasing SE breakdown is observed with decreasing temperature, the breakdown is already significant at higher temperatures for the equimolar mixture. Violations of the SE relation in pure supercooled water at different temperatures and pressures have been previously explained using the translational jump-diffusion (TJD) approach, which provides a fundamental molecular basis, by directly connecting the SE breakdown with jump-diffusion of the molecules. We have used the same TJD approach for explaining the SE breakdown for the methanol/water binary mixtures of compositions studied in the experiment over a wide range of temperatures between 220 K and 300 K. We have understood that the jump-diffusion is the key responsible factor for the SE breakdown. The maximum jump-diffusion contribution gives rise to the early SE breakdown observed for the equimolar mixture observed in the experiment. This study, therefore, provides molecular insight into the SE breakdown for the supercooled water/methanol binary mixture, as found in the experiment.

Collective dynamic dipole moment and orientation fluctuations, cooperative hydrogen bond relaxations, and their connections to dielectric relaxation in ionic acetamide deep eutectics: Microscopic insight from simulations

The paper reports a detailed simulation study on collective reorientational relaxation, cooperative hydrogen bond (H-bond) fluctuations, and their connections to dielectric relaxation (DR) in deep eutectic solvents made of acetamide and three uni-univalent electrolytes, lithium nitrate (LiNO3), lithium bromide (LiBr), and lithium perchlorate (LiClO4). Because cooperative H-bond fluctuations and ion migration complicate the straightforward interpretation of measured DR timescales in terms of molecular dipolar rotations for these conducting media which support extensive intra- and inter-species H-bonding, one needs to separate out the individual components from the overall relaxation for examining the microscopic origin of various timescales. The present study does so and finds that reorientation of ion-complexed acetamide molecules generates relaxation timescales that are in sub-nanosecond to nanosecond range. This explains in molecular terms the nanosecond timescales reported by recent giga-Hertz DR measurements. Interestingly, the simulated survival timescale for the acetamide-Li(+) complex has been found to be a few tens of nanosecond, suggesting such a cation-complexed species may be responsible for a similar timescale reported by mega-Hertz DR measurements of acetamide/potassium thiocyanate deep eutectics near room temperature. The issue of collective versus single particle relaxation is discussed, and jump waiting time distributions are determined. Dependence on anion-identity in each of the cases has been examined. In short, the present study demonstrates that assumption of nano-sized domain formation is not required for explaining the DR detected nanosecond and longer timescales in these media.

How Heterogeneous Are Trehalose/Glycerol Cryoprotectant Mixtures? A Combined Time-Resolved Fluorescence and Computer Simulation Investigation

Heterogeneity and molecular motions in representative cryoprotectant mixtures made of trehalose and glycerol are investigated in the temperature range 298 ≤ T (K) ≤ 353, via time-resolved fluorescence Stokes shift and anisotropy measurements, and molecular dynamics simulations of four-point density-time correlations and H-bond relaxations. Mixtures containing 5 and 20 wt % of trehalose along with neat glycerol are studied. Viscosity coefficients for these systems lie in the range 0.30 < η (P) < 23. Measured solute (Coumarin 153) rotation and solvation times reveal a substantial departure from the hydrodynamic viscosity dependence, suggesting the strong microheterogeneous nature of these systems. Fluorescence anisotropy decays are highly nonexponential, reflecting a non-Markovian character of the medium friction. A complete missing of the Stokes shift dynamics in these systems at 298 K but partial detection of it at other higher temperatures (shift magnitude being ∼400-600 cm^-1) indicates rigid solute environments. An amorphous solid-like feature emerges in the simulated radial distribution functions at these temperatures. Analyses of mean squared displacements reveal rattling-in-a-cage motion, non-Gaussian displacement distributions, and strong dynamic heterogeneity features. Simulated dynamic structure factors and four-point correlations hint, respectively, at very long α-relaxation and correlated time scales at 298 K. This explains the long solute rotation times (∼80-200 ns) measured at 298 K. Stretched exponential decay of the simulated H-bond relaxations with long time scales further highlights the strong temporal heterogeneity and slow dynamics inherent to these systems. In summary, this work provides the first insight into the molecular motions and interspecies interaction in a representative cryoprotectant mixture, and stimulates further study to investigate the interconnection between cryoprotection and dynamic heterogeneity.

Microscopic origin of temporal heterogeneities in translational dynamics of liquid water

Liquid water is known to reorient via a combination of large angular jumps (due to exchange of hydrogen bonding (H-bond) partners) and diffusive orientations. Translation of the molecule undergoing the orientational jump and its initial and final H-bond acceptors plays a key role in the microscopic reorientation process. Here, we partition the translational dynamics into those occurring during intervals when rotating water molecules (and their initial and final H-bonding partners) undergo orientational jump and those arising when molecules wait between consecutive orientational jumps. These intervals are chosen in such a way that none of the four possible H-bonds involving the chosen water molecule undergo an exchange process within its duration. Translational dynamics is analysed in terms of the distribution of particle displacements, van Hove functions, and its moments. We observe that the translational dynamics, calculated from molecular dynamics simulations of liquid water, is fastest during the orientational jumps and slowest during periods of waiting. The translational dynamics during all temporal intervals shows an intermediate behaviour. This is the microscopic origin of temporal dynamic heterogeneity in liquid water, which is mild at 300 K and systematically increases with supercooling. Study of such partitioned dynamics in supercooled water shows increased disparity in dynamics occurring in the two different types of intervals. Nature of the distribution of particle displacements in supercooled water is investigated and it reveals signatures non-Gaussian behaviour.