Are N-methyl groups of Tetramethylurea (TMU) Hydrophobic? A composition and temperature-dependent fluorescence spectroscopic investigation of TMU/water binary mixtures

Identical Diffusion Distributions and Co-Cluster Formation Dictate Azeotrope Formation: Microscopic Evidences and Experimental Signatures

What selects azeotropic pairs and governs the azeotropic conditions (composition and temperature) is an open and intriguing question. A combined simulation and experimental work presented here investigates this by considering ethanol-water mixtures. We find identical distributions of center-of-mass diffusion coefficients for ethanol and water molecules under the azeotropic condition (95.5 wt % ethanol +4.5 wt % water, Tazeo = 351.1K). Moreover, the particle displacements show strong interspecies correlations at Tazeo. Interestingly, simulated reorientation time distributions become identical at Tazeo but at a composition different from that at which the translational diffusion distributions overlapped. Cluster analyses indicate that solutions at Tazeo with xwater ≤ 15 wt % are more microheterogeneous than those with higher water content, although no anomaly in the composition-dependent solution structural properties was detected. Ethanol-water and ethanol-ethanol interaction energies show pronounced nonideal composition dependence, but the size of the relative fluctuations in them remained small (∼0.5kBT). Rare water-water H-bonding, predominant water-ethanol H-bonding, and a sizable population of "free" water molecules characterize the azeotropic solutions. The red edge excitation spectroscopic (REES) measurements with a dissolved anionic fluorescent dye, coumarin343 (C343), support the predicted solution microheterogeneity by showing a nonmonotonic composition dependence of the excitation energy-induced changes in the fluorescence emission spectral frequencies and bandwidths, the largest changes being under the azeotropic condition. Subsequent dynamic anisotropy measurements reveal a nonmonotonic composition dependence of C343 rotation times with a peak under the azeotropic condition. In summary, equalization of the component translational diffusion coefficients and solution microheterogeneity with regular composition dependence of the solution structure appear to characterize the ethanol-water azeotrope.

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.

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.

Non-ideal behavior of binary aqueous mixtures of some urea derivatives and their capacity to induce lysozyme gelation

The urea derivatives, namely, ethylurea (EU), 1,3 dimethylurea (1,3-DMU) and 1,1 diethylurea (1,1-DEU), in the limiting regions of their solubilities in water, and tetramethylurea (TMU) at w≥0.65 were investigated in relation to their capacity of inducing hen egg white lysozyme (HEWL) physical (non-covalent) gelation. Protein transparent gels were generated out of TMU/H₂O and 1,1-DEU/H₂O, respectively, whereas an intensively turbid gel resulted from sol-gel transition taking place in EU/H₂O. Oscillatory rheology revealed distinctions in the gels' structural and dynamic characteristics. Hydration patterns of the derivatives in solution, sizes of their non-polar domains and supramolecular symmetry features played a central role in their capacity of gel formation and in the gels' rheological behavior and morphology. Effects on gel characteristics of distinctively positioned ions in the Hofmeister series showed that SCN^- disrupted water H-bonding interconnectivity in TMU lysozyme gel, strengthening gel structure, yet maintaining gel transparency. Citrate enhanced system elasticity albeit causing intense turbidity and leading to phase separation. Larger values of the storage modulus, G', were verified for gels generated from binary mixtures containing urea derivatives with higher dipole moments.

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.