Temperature-dependent structure of methanol-water mixtures on cooling: X-ray and neutron diffraction and molecular dynamics simulations

Second inflection point of supercooled water surface tension induced by hydrogen bonds: A molecular-dynamics study

Surface tension of supercooled water is a fundamental property in various scientific processes. In this study, we perform molecular dynamics simulations with the TIP4P-2005 model to investigate the surface tension of supercooled water down to 220 K. Our results show a second inflection point (SIP) in the surface tension at temperature TSIP ≈ 267.5 ± 2.3 K. Using an extended IAPWS-E functional fit for the water surface tension, we calculate the surface excess internal-energy and entropy terms of the excess Helmholtz free energy. Similar to prior studies [Wang et al., Phys. Chem. Chem. Phys. 21, 3360 (2019); Gorfer et al., J. Chem. Phys. 158, 054503 (2023)], our results show that the surface tension is governed by two driving forces: a surface excess entropy change above the SIP and a surface excess internal-energy change below it. We study hydrogen-bonding near the SIP because it is the main cause of water's anomalous properties. With decreasing temperature, our results show that the entropy contribution to the surface tension reaches a maximum slightly below the SIP and then decreases. This is because the number of hydrogen bonds increases more slowly below the SIP. Moreover, the strengths and lifetimes of the hydrogen bonds also rise dramatically below the SIP, causing the internal-energy term to dominate the excess surface free energy. Thus, the SIP in the surface tension of supercooled TIP4P-2005 water is associated with an increase in the strengths and lifetimes of hydrogen bonds, along with a decrease in the formation rate (#/K) of new hydrogen bonds.

The liquidus temperature curve of aqueous methanol mixtures: a numerical simulation study

The liquidus temperature curve that characterizes the boundary between the liquid methanol/water mixture and its coexistence with ice Ih is determined using the direct-coexistence method. Several methanol concentrations and pressures of 0.1 MPa, 50 MPa, and 100 MPa are considered. In this study, we used the TIP4P/Ice model for water and two different models for methanol: OPLS and OPLS/2016, using the geometric rule for the Lennard-Jones cross interactions. We compared our simulation results with available experimental data and found that this combination of models reproduces reasonably well the liquidus curve for methanol mole fractions up to xm=0.3 at p=0.1 MPa. The freezing point depression of these mixtures is calculated and compared to experimental results. We also analyzed the effect of pressure on the liquidus curve, and we found that both models also reproduce qualitatively well the experimental decreasing of the liquidus temperatures as the pressure increases.

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.

Properties of Hydrogen-Bonded Networks in Ethanol-Water Liquid Mixtures as a Function of Temperature: Diffraction Experiments and Computer Simulations

New X-ray and neutron diffraction experiments have been performed on ethanol–water mixtures as a function of decreasing temperature, so that such diffraction data are now available over the entire composition range. Extensive molecular dynamics simulations show that the all-atom interatomic potentials applied are adequate for gaining insight into the hydrogen-bonded network structure, as well as into its changes on cooling. Various tools have been exploited for revealing details concerning hydrogen bonding, as a function of decreasing temperature and ethanol concentration, like determining the H-bond acceptor and donor sites, calculating the cluster-size distributions and cluster topologies, and computing the Laplace spectra and fractal dimensions of the networks. It is found that 5-membered hydrogen-bonded cycles are dominant up to an ethanol mole fraction x_eth = 0.7 at room temperature, above which the concentrated ring structures nearly disappear. Percolation has been given special attention, so that it could be shown that at low temperatures, close to the freezing point, even the mixture with 90% ethanol (x_eth = 0.9) possesses a three-dimensional (3D) percolating network. Moreover, the water subnetwork also percolates even at room temperature, with a percolation transition occurring around x_eth = 0.5.

Pressure-Dependent Structure of Methanol-Water Mixtures up to 1.2 GPa: Neutron Diffraction Experiments and Molecular Dynamics Simulations

Total scattering structure factors of per-deuterated methanol and heavy water, CD₃OD and D₂O, have been determined across the entire composition range as a function of pressure up to 1.2 GPa, by neutron diffraction. The largest variations due to increasing pressure were observed below a scattering variable value of 5 Å⁻¹, mostly as shifts in terms of the positions of the first and second maxima. Molecular dynamics computer simulations, using combinations of all-atom potentials for methanol and various water force fields, were conducted at the experimental pressures with the aim of interpreting neutron diffraction results. The peak-position shifts mentioned above could be qualitatively reproduced by simulations, although in terms of peak intensities, the accord between neutron diffraction and molecular dynamics was much less satisfactory. However, bearing in mind that increasing pressure must have a profound effect on repulsive forces between neighboring molecules, the agreement between experiment and computer simulation can certainly be termed as satisfactory. In order to reveal the influence of changing pressure on local intermolecular structure in these “simplest of complex” hydrogen-bonded liquid mixtures, simulated structures were analyzed in terms of hydrogen bond-related partial radial distribution functions and size distributions of hydrogen-bonded cyclic entities. Distinct differences between pressure-dependent structures of water-rich and methanol-rich composition regions were revealed.

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.