Structural anomaly and dynamic heterogeneity in cycloether/water binary mixtures: Signatures from composition dependent dynamic fluorescence measurements and computer simulations

Water Fraction Dependence of the Aggregation Behavior of Hydrophobic Fluorescent Solutes in Water-Tetrahydrofuran

This work investigates the water fraction dependence of the aggregation behavior of hydrophobic solutes in water-tetrahydrofuran (THF) and the elucidation of the role of THF using fluorescence microscopy, dynamic light scattering, neutron and X-ray scattering, and photoluminescence measurements. On the basis of the obtained results, the following model is proposed: hydrophobic molecules are molecularly dispersed in the low-water-content region (10-20 vol %), while they form mesoscopic particles upon increasing the water fraction to ∼30 vol %. This abrupt change is due to the composition fluctuation of the water-THF binary system to form hydrophobic areas in THF, followed by THF-rich droplets where hydrophobic solutes are incorporated and form loose aggregates. Further increasing the water content prompts the desolvation of THF, which decreases the particle size and generates tight aggregates of solute molecules. This model is consistent with the luminescence behavior of the solutes and will be helpful to control the aggregation state of hydrophobic solutes in various applications.

Does Viscosity Decoupling Guarantee Dynamic Heterogeneity? A Clue through an Excitation and Emission Wavelength-Dependent Time-Resolved Fluorescence Anisotropy Study

Traditionally, deviation from Stokes-Einstein-Debye (SED) relation in terms of viscosity dependence of medium dynamics, i.e., τ x ∝ ( η T ) p with p ≠ 1, is taken as a signature of dynamic heterogeneity. However, it does not guarantee medium heterogeneity, as the decoupling may also originate from the deviation of the basic assumption of SED. Here, we developed a method to find a stronger relation between viscosity decoupling (p ≠ 1) and dynamic heterogeneity in terms of rotational motion. Our approach exploited the fact that in heterogeneous media, a solvatochromic probe will be solvated to a different extent at different microdomains (subpopulations), and photoselection of these subpopulations can be achieved by excitation or emission wavelength-dependent measurements. We hypothesized that the dynamics of a homogeneous system might show viscosity decoupling, but the extent of decoupling at different excitations (or at different emissions) should not be different. On the other hand, in a heterogeneous medium, this extent of viscosity decoupling (p-value) should be different at different excitations (or at different emissions). As proof of concept, we investigated three versatile solvent media: squalane (viscous molecular liquid), 1-ethyle-3-methylimidazolium ethyl sulfate ionic liquid (IL), and [0.78 acetamide + 0.22 LiNO3] deep eutectic solvent (DES). We found that squalane is homogeneous, although it shows fractional viscosity dependence (p ≠ 1). Interestingly, mild heterogeneity in IL and significant heterogeneity in the DES were observed. Overall, we conclude that the difference in the p-value as a function of excitation (or emission) wavelength-dependent might be a superior way for the detection of dynamic heterogeneity.

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.

Dynamical Anomaly of Aqueous Amphiphilic Solutions: Connection to Solution H-Bond Fluctuation Dynamics?

We have investigated the possible connection between "dynamical anomaly" observed in time-resolved fluorescence measurements of reactive and nonreactive solute-centered relaxation dynamics in aqueous binary mixtures of different amphiphiles and the solution intra- and interspecies H-bond fluctuation dynamics. Earlier studies have connected the anomalous thermodynamic properties of binary mixtures at very low amphiphile concentrations to the structural distortion of water. This is termed as "structural anomaly." Interestingly, the abrupt changes in the composition-dependent average rates of solute relaxation dynamics occur at amphiphile mole fractions approximately twice as large as those where structural anomalies appear. We have investigated this anomalous solution dynamical aspect by considering (water + tertiary butanol) as a model system and performed molecular dynamics simulations at several tertiary butanol (TBA) concentrations covering the extremely dilute to the moderately concentrated regimes. The "dynamical anomaly" has been followed via monitoring the composition dependence of the intra- and interspecies H-bond fluctuations and reorientational relaxations of TBA and water molecules. Solution structural aspects have been followed via examining the tetrahedral order parameter, radial and spatial distribution functions, numbers of H bonds per water and TBA molecules, and the respective populations participating in H-bond formation. Our simulations reveal abrupt changes in the H-bond fluctuations and reorientational dynamics and tetrahedral order parameter at amphiphile concentrations differing approximately by a factor of 2 and corroborates well with the steady-state and the time-resolved spectroscopic measurements. This work therefore explains, following a uniform and cogent manner, both the experimentally observed structural and dynamical anomalies in microscopic terms.

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.

Effects of molecular shape on alcohol aggregation and water hydrogen bond network behavior in butanol isomer solutions

Despite butanol isomers such as n-butanol, sec-butanol, isobutanol and tert-butanol having the same chemical formula, their liquid-liquid phase diagrams are distinct. That is, tert-butanol is miscible in water at all concentrations, while the other three butanol isomers are partially miscible under ambient conditions. The molecular shape of tert-butanol is close to globular and differs from the other three butanol molecules with a relatively long carbon chain. By performing molecular dynamics simulations and graph theoretical analysis of the four water-butanol isomer mixtures at varying concentrations, we show how distinct butanol aggregates are formed which depend upon the molecular shape and affect the water H-bond network structure and phase diagram in the binary liquid. The three butanol isomers of n-butanol, sec-butanol and isobutanol at concentrated solutions form chain-like alcohol aggregates, but tert-butanol forms small aggregates due to the distinct packing behavior caused by its globular molecular shape. By employing the graph theoretical analysis such as the degree distribution and the eigenvalue spectrum from the adjacency matrix in the graphical representation of the alcohol H-bond network, we show that the tert-butanol aggregates have a different morphological structure from that of the other three butanol isomers in aqueous solution. The graph theoretically distinct butanol aggregates are categorized into two groups, water-compatible and water-incompatible, depending upon the interaction between the alcohol and water molecules. Based upon our observations, we propose that the water-incompatible networks of n-butanol, sec-butanol and isobutanol aggregates do not change the water structure significantly, forming two separate liquid phases that are alcohol-rich and water-rich. However, the water-compatible network of tert-butanol aggregates has a considerable interaction with the water molecules and causes significant disruption of the water H-bond network, forming a homogeneous solution. Understanding the alcohol aggregation behavior and water structure in butanol-water mixtures provides a critical clue in appreciating fundamental issues such as miscibility and phase separation in aqueous solution systems.

Molecular heterogeneity in aqueous cosolvent systems

Aqueous cosolvent systems (ACoSs) are mixtures of small polar molecules such as amides, alcohols, dimethyl sulfoxide, or ions in water. These liquids have been the focus of fundamental studies due to their complex intermolecular interactions as well as their broad applications in chemistry, medicine, and materials science. ACoSs are fully miscible at the macroscopic level but exhibit nanometer-scale spatial heterogeneity. ACoSs have recently received renewed attention within the chemical physics community as model systems to explore the relationship between intermolecular interactions and microscopic liquid-liquid phase separation. In this perspective, we provide an overview of ACoS spatial segregation, dynamic heterogeneity, and multiscale relaxation dynamics. We describe emerging approaches to characterize liquid microstructure, H-bond networks, and dynamics using modern experimental tools combined with molecular dynamics simulations and network-based analysis techniques.

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.

Microscopic origin of breakdown of Stokes-Einstein relation in binary mixtures: Inherent structure analysis

Aqueous binary mixtures often exhibit dramatic departure from the predicted hydrodynamic behavior when transport properties are plotted against composition. We show by inherent structure (IS) analysis that this sharp composition dependent breakdown of the Stokes-Einstein relation can be attributed to the non-monotonic variation in the average inherent structure energy of these mixtures. Further IS analysis reveals the existence of a unique ground state, stabilized by both the formation of an optimum number of H-bonds and a favorable hydrophobic interaction at this composition. The surprisingly sharp turnaround behavior observed in the effective hydrodynamic radius also owes its origin to the same combination of these two factors. Interestingly, the temperature dependence of isothermal compressibility shows a minimum at the particular composition. Extensive studies on water-dimethyl sulfoxide and water-ethanol mixtures using two different force-fields of water reveal many features that are nearly universal. A justification of this quasi-universal behavior is provided in terms of a mode-coupling theory (MCT) of viscosity, which can serve as the starting point of a remarkable correlation observed with the nearest neighbor structure, as captured by the first peaks of the radial distribution function, and the slowdown in the intermediate scattering function at intermediate wavenumbers. Therefore, the formation of the local structure captured through IS analysis can be correlated with the MCT.

Cloud Point Driven Dynamics in Aqueous Solutions of Thermoresponsive Copolymers: Are They Akin to Criticality Driven Solution Dynamics?

Cloud point driven interaction and relaxation dynamics of aqueous solutions of amphiphilic thermoresponsive copolymers were explored through picosecond resolved and steady state fluorescence measurements employing hydrophilic (coumarin 343, C343) and hydrophobic (coumarin 153, C153) solute probes of comparable sizes. These thermoresponsive random copolymers, with tunable cloud point temperatures (T_cp's) between 298 and 323 K, were rationally designed first and then synthesized via reversible addition-fragmentation chain transfer (RAFT) copolymerization of methyl methacrylate (MMA) and poly(ethylene glycol) monomethyl ether methacrylate (PEGMA). Subsequently, copolymers were characterized by NMR spectroscopy and size exclusion chromatography (SEC). A balance between the hydrophilic (PEGMA) and the hydrophobic (MMA) content dictates the critical aggregation concentration (CAC), with CAC ∼ 2-14 mg/L for these copolymers in aqueous media. No abrupt changes in the steady state spectral features of both C153 and C343 in the aqueous solutions of these polymers near but below the cloud point temperatures were observed. Interestingly, spectral properties of C153 in these solutions show the impact of hydrophobic/hydrophilic interaction balance but not by those of C343. More specifically, C153 reported a blue shift (relative to that in neat water) and heterogeneity in its local environment. This suggested different locations for the hydrophilic (C343) and the hydrophobic (C153) probes. In addition, the excited state fluorescence lifetime (⟨τ_life⟩) of C153 increased with the increase of hydrophobic (MMA) content in these copolymers. However, C343 reported no such variations, although fluorescence anisotropy decays for both solutes were significantly slowed down in these aqueous solutions compared to neat water. Anisotropy decays indicated bimodal time-dependent friction for these solutes in aqueous solutions of these copolymers but monomodal in neat water. A linear dependence of the average rotational relaxation rates (⟨k_rot⟩ = ⟨τ_rot⟩^-1) of the type ⟨k_rot⟩ ∝ (|T - T_cp|/T_cp)^γ with negative values for the exponent γ was observed for both solutes. No slowing down of the solute rotation with temperature approaching the T_cp was detected; rather, rotation became faster upon increasing the solution temperature, suggesting domination of the local friction.

Interaction and Dynamics in a Fully Biodegradable Glucose-Containing Naturally Abundant Deep Eutectic Solvent: Temperature-Dependent Time-Resolved Fluorescence Measurements

A new room-temperature deep eutectic solvent (DES) composed of glucose, urea, and water has been prepared and its relaxation dynamics explored via temperature-dependent time-resolved fluorescence measurements employing hydrophilic and hydrophobic solute probes. Differential scanning calorimetry measurements indicate a glass transition temperature (T_g) of ∼236 K. Measured viscosity coefficients (η) vary from ∼600 to ∼100 cP in the temperature range 318 ≤ T/K ≤ 343 and exhibit Arrhenius-type temperature dependence with an activation energy of ∼65 kJ mol^-1. Interestingly, this DES forms a stable liquid at ∼300 K but is too viscous to be accurately measured by us below 318 K. Temperature-dependent dynamic fluorescence anisotropy measurements using hydrophobic and hydrophilic solutes of similar sizes reveal bi-exponential kinetics and Arrhenius-type temperature dependence for solute rotation times (⟨τ_r⟩) but with significantly decreased activation energies, ∼31 kJ mol^-1 (hydrophobic) and ∼21 kJ mol^-1 (hydrophilic). Deviation from hydrodynamics is further reflected in the strong fractional viscosity dependence of ⟨τ_r⟩: ⟨τ_r⟩ ∝ (η/T)^p with p ≈ 0.3-0.5, indicating pronounced temporal heterogeneity in the relaxation dynamics. Dynamic fluorescence Stokes shift measurements (temporal resolution ∼85 ps) produce dynamic shifts of ∼500-700 cm^-1, bi-exponential solvation energy relaxation with time constants in the range ∼0.2 ns and ∼4 ns, and estimated missing amplitudes of ∼65-75%. Impact of the density difference between a nonpolar solvent and this DES on the estimated missing amplitudes is explored via measuring the temperature-dependent densities and refractive indices of this DES. Lifetime measurements suggest considerable temperature dependence for the hydrophobic solute but no such dependence for the hydrophilic one. Excitation energy dependence of fluorescence emission of various solutes with widely different lifetimes indicates mild spatial heterogeneity for this DES.

Formation and dissociation of synthetic hetero-double-helix complex in aqueous solutions: significant effect of water content on dynamics of structural change

A 1 : 1 mixture of the ethynylhelicene pseudoenantiomers (M)-tetramer and (P)-pentamer, which possess hydrophilic terminal tri(ethyleneglycol) (TEG) groups, changes their structures in the water-THF (10 μM) solvent system between dissociated random-coils and an associated hetero-double-helix upon heating and cooling. A small change in water content between 30 and 33% significantly affects the dynamics of structural changes. At 30% water content, heating to 60 °C causes rapid formation of random-coil and cooling to 10 °C causes the rapid formation of hetero-double-helix, accompanied by repeated changes in Δε at 369 nm between 0 and -2000 cm-1 M-1. Heating and cooling experiments at constant rates between 60 and 10 °C resulted in sigmoidal curves in Δε/temperature profiles, which indicate rapid structural changes. Different phenomena occurred at 33% water content. Heating to 60 °C and cooling to 0 °C initially induced changes in Δε between 0 and -2000 cm-1 M-1, and repeated cycles gradually reduced the range between 0 and -500 cm-1 M-1. Heating and cooling experiments at constant rates between 60 and 10 °C caused small changes in Δε, and repeated cycles at 10 °C gradually increased Δε to -500 cm-1 M-1. These phenomena involved rapid changes in molecular structure and slow structural changes in the water-THF solvent system. The sharp switching of the dynamics of structural changes at water content between 30 and 33% indicated discontinuous structural changes in the hydration of TEG and/or in water clusters in the vicinity of oligomer molecules.

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.