Hydrophobic Segregation, Phase Transitions and the Anomalous Thermodynamics of Water/Methanol Mixtures

Solvation shell thermodynamics of extended hydrophobic solutes in mixed solvents

The ability of various cosolutes and cosolvents to enhance or quench solvent density fluctuations at solute–water interfaces has crucial implications on the conformational equilibrium of macromolecules such as polymers and proteins. Herein, we use an extended hydrophobic solute as a model system to study the effect of urea and methanol on the density fluctuations in the solute’s solvation shell and the resulting thermodynamics. On strengthening the solute–water/cosolute repulsive interaction, we observe distinct trends in the mutual affinities between various species in, and the thermodynamic properties of, the solvation shell. These trends strongly follow the respective trends in the preferential adsorption of urea and methanol: solute–water/cosolute repulsion strengthens, urea accumulation decreases, and methanol accumulation increases. Preferential accumulation of urea is found to quench the density fluctuations around the extended solute, leading to a decrease in the compressibility of the solvation shell. In contrast, methanol accumulation enhances the density fluctuations, leading to an increase in the compressibility. The mode of action of urea and methanol seems to be strongly coupled to their hydration behavior. The observations from this simple model is discussed in relation to urea driven swelling and methanol induced collapse of some well-known thermo-responsive polymers.

Dynamics of molecular associates in methanol/water mixtures

The dynamics of molecular associates in a methanol/water mixture was investigated using quasielastic neutron scattering. By measuring the signal from four methanol/water samples differing only by their isotopic composition, the relative motion of the water to methanol molecules, i.e. their mutual dynamics, was determined at the nanoscale. The thus obtained nanoscopic mutual diffusion coefficient signals a significantly slower process than the single particle diffusion of either methanol or water in the system as well as their macroscopic mutual diffusion. The data do not provide any indication of microsegregation in this preeminent alcohol/water mixture; however, they do indicate the existence of long lived but dynamic molecular associates of water and methanol molecules. Analysis of the structural relaxation shows that the lifetime of molecular association through hydrogen bonding determines the fact that viscosity of the mixtures at intermediate concentrations is higher than that of both pure components.

Thermodynamic properties of a molecular dipolar liquid using the two-phase thermodynamic approach

A modified version of the original two-phase thermodynamic approach has been extended to evaluate the thermodynamic properties of molecular systems. Its basic assumption states that the density of states can be decomposed into solid-like and gas-like components. The solid part has been approximated by that of a set of harmonic oscillators, whereas a subset composed of rough hard spheres has been considered for the gas part. In this new approach, molecules have been modelled as rotating hard spheres that experience elastic collisions. The technique has been tested on a system made up of dipolar diatomic molecules, and it leads to very good results for total entropy, potential energy reference and heat capacity. Translation and rotation solid components of the overall spectra have been compared to the real part of the instantaneous normal mode vibrational densities of states. Similarities between them reinforce the validity of the two-phase thermodynamic approach.

Methane Clathrate Formation is Catalyzed and Kinetically Inhibited by the Same Molecule: Two Facets of Methanol

Here, we perform molecular dynamics simulations to provide atomic-level insights into the dual roles of methanol in enhancing and delaying the rate of methane clathrate hydrate nucleation. Consistent with experiments, we find that methanol slows clathrate hydrate nucleation above 250 K but promotes clathrate formation at temperatures below 250 K. We show that this behavior can be rationalized by the unusual temperature dependence of the methane-methanol interaction in an aqueous solution, which emerges due to the hydrophobic effect. In addition to its antifreeze properties at temperatures above 250 K, methanol competes with water to interact with methane prior to the formation of clathrate nuclei. Below 250 K, methanol encourages water to occupy the space between methane molecules favoring clathrate formation and it may additionally promote water mobility.

Molecular-scale origins of solution nanostructure and excess thermodynamic properties in a water/amphiphile mixture

The molecular and nanoscale origins of nonideality in excess thermodynamic properties are essential to understanding cosolvent mixtures, yet they remain challenging to determine. Here, we consider a binary mixture of water and an amphiphile, N,N,N',N'-tetramethylmalonamide (TMMA), which is characterized by strong hydrogen bonding between the two components and no hydrogen bonding between amphiphiles. Using molecular dynamics simulation, validated with excess volume measurements and X-ray scattering, we identify three distinct solution regimes across the composition range of the binary mixture and find that the transition between two of these regimes, marked by the water percolation threshold, is closely correlated with minima in the excess volume and excess enthalpy. Structural analysis of the simulations reveals an interplay between local interactions and solution nanostructure, determined by the relative strength of the water-water and water-amphiphile hydrogen bonding interactions. By comparison with other amphiphiles, such as linear alcohols, the relative strength of like and unlike interactions between water and amphiphile affects the relationship between thermodynamics and structural regimes. This provides insight into how molecular forces of mutual solvation interact across length scales and how they manifest in excess thermodynamic properties.

Determination of the Free Energies of Mixing of Organic Solutions through a Combined Molecular Dynamics and Bayesian Statistics Approach

As new generations of thin-film semiconductors are moving toward solution-based processing, the development of printing formulations will require information pertaining to the free energies of mixing of complex mixtures. From the standpoint of in silico material design, this move necessitates the development of methods that can accurately and quickly evaluate these formulations in order to maximize processing speed and reproducibility. Here, we make use of molecular dynamics (MD) simulations, in combination with the two-phase thermodynamic (2PT) model, to explore the free energy of mixing surfaces for a series of halogenated solvents and high-boiling point solvent additives used in the development of thin-film organic semiconductors. Although the combined methods generally show good agreement with available experimental data, the computational cost to traverse the free-energy landscape is considerable. Hence, we demonstrate how a Bayesian optimization scheme, coupled with the MD and 2PT approaches, can drastically reduce the number of simulations required, in turn shrinking both the computational cost and time.

Molecular Mechanism of Water Reorientation Dynamics in Dimethyl Sulfoxide Aqueous Mixtures

Nonmonotonic composition dependence is often observed for numerous properties in the aqueous mixtures of small amphiphilic molecules. The molecular picture underlying this structure-activity relationship, however, remains largely elusive. We herein studied water reorientation dynamics in the aqueous mixture of dimethyl sulfoxide (DMSO), which has a significant nonmonotonic composition dependence, using molecular dynamic simulation and an extended molecular jump model. The analysis indicates that this nonideal behavior is driven by the collective frame diffusion component of water reorientation, which decelerates in the water-rich regime because of the strengthened hydrogen bonds and accelerates in the water-poor regime as the hydrogen bonding network is broken into smaller aggregates. The current work therefore connects the microheterogeneity in the solvation structure of DMSO-water with its nonmonotonic hydration dynamics and sheds new light on how microsegregation leads to the multiscale hydration nonideality in general.

Exploring the hydrogen bond kinetics of methanol–water solutions using Raman scattering

Stimulated Raman scattering was used to clearly show the hydrogen bond kinetics of water–methanol mixed solutions.

Accurate calculation of zero point energy from molecular dynamics simulations of liquids and their mixtures

The two-phase thermodynamic (2PT) method is used to compute the zero point energy (ZPE) of several liquids and their mixtures. The 2PT method uses the density of states (DoS), which is computed from the velocity autocorrelation (VAC) function obtained from a short classical molecular dynamics trajectory. By partitioning the VAC and the DoS of a fluid into solid and gaslike components, quantum mechanical corrections to thermodynamical properties can be computed. The ZPE is obtained by combining the partition function of the quantum harmonic oscillator with the vibrational part of the solidlike DoS. The resulting ZPE is found to be in excellent agreement with both experimental and ab initio results. Solvent effects such as hydrogen bonding and polarization can be included by the utilization of ab initio density functional theory based molecular dynamics simulations. It is found that these effects significantly influence the DoS of water molecules. The obtained results demonstrate that the 2PT model is a powerful method for efficient ZPE calculations, in particular, to account for solvent effects and polarization.

Transferability of Local Density-Assisted Implicit Solvation Models for Homogeneous Fluid Mixtures

The application of bottom-up coarse grained (CG) models to study the equilibrium mixing behavior of liquids is rather challenging, since these models can be significantly influenced by the density or the concentration of the state chosen during parametrization. This dependency leads to low transferability in density/concentration space and has been one of the major limitations in bottom-up coarse graining. Recent approaches proposed to tackle this shortcoming range from the addition of thermodynamic constraints, to an extended ensemble parametrization, to the addition of supplementary terms to the system's Hamiltonian. To study fluid phase equilibria with bottom-up CG models, the application of local density (LD) potentials appears to be a promising approach, as shown in previous work by Sanyal and Shell [T. Sanyal, M. S. Shell, J. Phys. Chem. B, 2018, 122, 5678]. Here, we want to further explore this method and test its ability to model a system which contains structural inhomogeneities only on the molecular scale, namely, solutions of methanol and water. We find that a water-water LD potential improves the transferability of an implicit-methanol CG model toward high water concentration. Conversely, a methanol-methanol LD potential does not significantly improve the transferability of an implicit-water CG model toward high methanol concentration. These differences appear due to the presence of cooperative interactions in water at high concentrations that the LD potentials can capture. In addition, we compare two different approaches to derive our CG models, namely, relative entropy optimization and the Inverse Monte Carlo method, and formally demonstrate under which analytical and numerical assumptions these two methods yield equivalent results.

Molecular dynamics simulations of nano-confined methanol and methanol-water mixtures between infinite graphite plates: Structure and dynamics

Molecular dynamics simulations are used to investigate microscopic structures and dynamics of methanol and methanol-water binary mixture films confined between hydrophobic infinite parallel graphite plate slits with widths, H, in the range of 7-20 Å at 300 K. The initial geometric densities of the liquids were chosen to be the same as bulk methanol at the same temperature. For the two narrowest slit widths, two smaller initial densities were also considered. For the nano-confined system with H = 7 Å and high pressure, a solid-like hexagonal arrangement of methanol molecules arranged perpendicular to the plates is observed which reflects the closest packing of the molecules and partially mirrors the structure of the underlying graphite structure. At lower pressures and for larger slit widths, in the contact layer, the methanol molecules prefer having the C-O bond oriented parallel to the walls. Layered structures of methanol parallel to the wall were observed, with contact layers and additional numbers of central layers depending on the particular slit width. For methanol-water mixtures, simulations of solutions with different composition were performed between infinite graphite slits with H = 10 and 20 Å at 300 K. For the nanoslit with H = 10 Å, in the solution mixtures, three layers of molecules form, but for all mole fractions of methanol, methanol molecules are excluded from the central fluid layer. In the nanopore with H = 20 Å, more than three fluid layers are formed and methanol concentrations are enhanced near the confining plates walls compared to the average solution stoichiometry. The self-diffusion coefficients of methanol and water molecules in the solution show strong dependence on the solution concentration. The solution mole fractions with minimal diffusivity are the same in confined and non-confined bulk methanol-water mixtures.

Synthesis and EPR-spectroscopic characterization of the perchlorotriarylmethyl tricarboxylic acid radical (PTMTC) and its 13C labelled analogue (13C-PTMTC)

A hydrophilic tris(tetrachlorotriaryl)methyl (tetrachloro-TAM) radical labelled 50% with ¹³C at the central carbon atom was prepared. The mixture of isotopologue radicals was characterised by continuous wave and pulsed X-band electron paramagnetic spectroscopy (EPS). For the pharmaceutical and medical applications planned, the quantitative influence of oxygen, viscosity, temperature and pH on EPR line widths was studied in aqueous buffer, DMSO, water-methanol and water-glycerol mixtures. Under in vivo conditions, pH can be disregarded. There is a clear oxygen dependence of the width of the ¹²C isotopologue single EPR line in aqueous solutions while changes in rotational motion (viscosity) are observable only in the doublet lines of the central carbon of the ¹³C isotopologue. The tetrachloro-TAM proved to be very stable as a solid. Its thermal decay was determined quantitatively by thermal annealing. Towards ascorbic acid as a reducing agent and towards an oocyte cell extract it had a half-life of approx. 60 and 10 min. Thus for in vivo applications, 50% ¹³C tetrachloro-TAMs are suitable for selective and simultaneous oxygen and macroviscosity measurements in a formulation, e.g. nanocapsules.

Accurate schemes for calculation of thermodynamic properties of liquid mixtures from molecular dynamics simulations

We explore different schemes for improved accuracy of entropy calculations in aqueous liquid mixtures from molecular dynamics (MD) simulations. We build upon the two-phase thermodynamic (2PT) model of Lin et al. [J. Chem. Phys. 119, 11792 (2003)] and explore new ways to obtain the partition between the gas-like and solid-like parts of the density of states, as well as the effect of the chosen ideal "combinatorial" entropy of mixing, both of which have a large impact on the results. We also propose a first-order correction to the issue of kinetic energy transfer between degrees of freedom (DoF). This problem arises when the effective temperatures of translational, rotational, and vibrational DoF are not equal, either due to poor equilibration or reduced system size/time sampling, which are typical problems for ab initio MD. The new scheme enables improved convergence of the results with respect to configurational sampling, by up to one order of magnitude, for short MD runs. To ensure a meaningful assessment, we perform MD simulations of liquid mixtures of water with several other molecules of varying sizes: methanol, acetonitrile, N, N-dimethylformamide, and n-butanol. Our analysis shows that results in excellent agreement with experiment can be obtained with little computational effort for some systems. However, the ability of the 2PT method to succeed in these calculations is strongly influenced by the choice of force field, the fluidicity (hard-sphere) formalism employed to obtain the solid/gas partition, and the assumed combinatorial entropy of mixing. We tested two popular force fields, GAFF and OPLS with SPC/E water. For the mixtures studied, the GAFF force field seems to perform as a slightly better "all-around" force field when compared to OPLS+SPC/E.

Water-Methanol Mixtures: Simulations of Mixing Properties over the Entire Range of Mole Fractions

Numerous experimental and theoretical investigations have been devoted to the hydrogen bond in pure liquids and mixtures. Among the different theoretical approaches, molecular dynamics (MD) simulations are predominant in obtaining detailed information, on the molecular level, simultaneously on the structure and the dynamics. Water and methanol are the two most prominent hydrogen-bonded liquids, and they and their mixtures have consequently been the subject of many studies; we revisit here the problem of the mixtures. An important first step is to check whether a classical potential model, the components of which are deemed to be satisfactory for the pure liquids, is able to reproduce the known thermodynamic excess properties of the mixtures sufficiently well. We have used the available BJH (water) and PHH (methanol) flexible models because they are by construction mutually compatible and also well suited to study, in a second step, some dynamic property characteristic of hydrogen-bonded liquids. In this article we show that these models, after a slight reparametrization for use in NpT simulations, reproduce the essential features of the excess mixing and molar properties of water-methanol mixtures. Furthermore, in the pure liquids, the agreement of the radial distribution functions with experiment remains as satisfactory as before. Similarly, the translation self-diffusion coefficients D are modified by less than 10%. In the mixtures, they evolve nonmonotonously as a function of mole fraction.

Optimizing the oxygen evolution reaction for electrochemical water oxidation by tuning solvent properties

Electrochemical water-based energy cycles provide a most promising alternative to fossil-fuel sources of energy. However, current electrocatalysts are not adequate (high overpotential, lack of selectivity toward O2 production, catalyst degradation). We propose here mechanistic guidelines for experimental examination of modified catalysts based on the dependence of kinetic rates on the solvent dielectric constant. To illustrate the procedure we consider the fcc(111) platinum surface and show that the individual steps for the oxygen evolution reaction (OER) and oxygen reduction reaction (ORR) change systematically with the polarizability of the medium. Thus changing this environmental variable can be used to tune the rate determining steps and the barriers, providing a means for screening and validating new systems to optimize the rate determining steps for the ORR and OER reaction pathways.

High pressure Raman spectroscopy of H2O-CH3OH mixtures

Complex intra-molecular interactions and the hydrogen-bonding network in H2O-volatile mixtures play critical roles in many dynamics processes in physical chemistry, biology, and Earth and planetary sciences. We used high pressure Raman spectroscopy to study the pressure evolution of vibrational frequencies and bonding behavior in H2O-CH3OH mixtures. We found that the presence of low CH3OH content in H2O increases the transition pressure where water crystallizes to ice VI, but does not significantly change the pressure where ice VI transforms to ice VII. Furthermore, the stiffening rates of C-H stretching frequencies dω/dP in CH3OH significantly decrease upon the crystallization of water, and the softening rates of the O-H stretching frequencies of ice VII are suppressed over a narrow pressure range, after which the frequencies of these modes shift with pressure in ways similar to pure CH3OH and ice VII, respectively. Such complex pressure evolution of Raman frequencies along with pronounced variations in Raman intensities of CH3OH within the sample, and the hysteresis of the water-ice VI phase transition suggest pressure-induced segregation of low content CH3OH from ice VII. These findings indicate the significant influence of volatiles on the crystallization of sub-surface ocean and thermal evolution within large icy planets and satellites.

Structural properties of methanol–water binary mixtures within the quantum cluster equilibrium model

The Quantum Cluster Equilibrium (QCE) method computes cluster distributions and thermodynamic properties of binary methanol–water mixtures in agreement with experiments.

Experiments on tracer diffusion in aqueous and non-aqueous solvent combinations

Forced Rayleigh scattering is used to study the tracer diffusion of an azobenzene in binary combinations of polar solvents, including water. In the absence of water, the tracer diffusion coefficient D in the mixture lies between the diffusion coefficients within the pure solvents, on a curve that is reasonably close to the prediction of free-volume theory. If water is present, on the other hand, the diffusion coefficient displays a minimum that is less than the smaller of the two pure-solvent values. We attempt to understand the different behavior in water by concentrating on the fairly hydrophobic nature of the solute, leading to a first solvent shell that is hydrophobic on the inside and hydrophilic on the outside. We also believe that clusters of amphiphiles explain the observation that, in aqueous combinations, D is nearly constant above a certain amphiphile mole fraction.

Dramatic increase in the oxygen reduction reaction for platinum cathodes from tuning the solvent dielectric constant

Hydrogen fuel cells (FC) are considered essential for a sustainable economy based on carbon-free energy sources, but a major impediment are the costs. First-principles quantum mechanics (density functional theory including solvation) is used to predict how the energies and barriers for the mechanistic steps of the oxygen reduction reaction (ORR) over the fcc(111) platinum surface depend on the dielectric constant of the solvent. The ORR kinetics can be strongly accelerated by decreasing the effective medium polarizability from the high value it has in water. Possible ways to realize this experimentally are suggested. The calculated volcano structure for the dependence of rate on solvent polarization is considered to be general, and should be observed in other electrochemical systems.

Low-temperature water dynamics in an aqueous methanol solution

An aqueous methanol solution (x(MeOH) = 0.30) has been studied by quasielastic neutron scattering. The single-particle water dynamics were effectively isolated by employing deuterated methanol. A smooth dynamic transition to a sub-Arrhenius temperature dependence has been observed in the relaxation times. We associate this behavior with the formation of small crystallites in the system. These findings are compared with molecular dynamics simulations and previous nuclear magnetic resonance measurements. We discuss possible dynamic signatures of structuring in the mixture.

Methanol incorporation in clathrate hydrates and the implications for oil and gas pipeline flow assurance and icy planetary bodies

One of the best-known uses of methanol is as antifreeze. Methanol is used in large quantities in industrial applications to prevent methane clathrate hydrate blockages from forming in oil and gas pipelines. Methanol is also assigned a major role as antifreeze in giving icy planetary bodies (e.g., Titan) a liquid subsurface ocean and/or an atmosphere containing significant quantities of methane. In this work, we reveal a previously unverified role for methanol as a guest in clathrate hydrate cages. X-ray diffraction (XRD) and NMR experiments showed that at temperatures near 273 K, methanol is incorporated in the hydrate lattice along with other guest molecules. The amount of included methanol depends on the preparative method used. For instance, single-crystal XRD shows that at low temperatures, the methanol molecules are hydrogen-bonded in 4.4% of the small cages of tetrahydrofuran cubic structure II hydrate. At higher temperatures, NMR spectroscopy reveals a number of methanol species incorporated in hydrocarbon hydrate lattices. At temperatures characteristic of icy planetary bodies, vapor deposits of methanol, water, and methane or xenon show that the presence of methanol accelerates hydrate formation on annealing and that there is unusually complex phase behavior as revealed by powder XRD and NMR spectroscopy. The presence of cubic structure I hydrate was confirmed and a unique hydrate phase was postulated to account for the data. Molecular dynamics calculations confirmed the possibility of methanol incorporation into the hydrate lattice and show that methanol can favorably replace a number of methane guests.