Physical Properties and Hydrogen-Bonding Network of Water-Ethanol Mixtures from Molecular Dynamics Simulations

Spidroins under the Influence of Alcohol: Effect of Ethanol on Secondary Structure and Molecular Level Solvation of Silk-Like Proteins

Future sustainable materials based on designer biomolecules require control of the solution assembly, but also interfacial interactions. Alcohol treatments of protein materials are an accessible means to this, making understanding of the process at the molecular level of seminal importance. We focus here on the influence of ethanol on spidroins, the main proteins of silk. By large-scale atomistically detailed molecular dynamics (MD) simulations and interconnected experiments, we characterize the protein aggregation, secondary structure changes, molecular level origins of them, and solvation environment changes for the proteins, as induced by ethanol as a solvation additive. The MD and circular dichoroism (CD) findings jointly show that ethanol promotes ordered structure in the protein molecules, leading to an increase of helix content and turns but also increased aggregation, as revealed by dynamic light scattering (DLS) and light microscopy. The structural changes correlate at the molecular level with increased intramolecular hydrogen bonding. The simulations reveal that polar amino acids, such as glutamine and serine, are most influenced by ethanol, whereas glycine residues are most prone to be involved in the ethanol-induced secondary structure changes. Furthermore, ethanol engages in interactions with the hydrophobic alanine-rich regions of the spidroin, significantly decreasing the hydrophobic interactions of the protein with itself and its surroundings. The protein solutes also change the microstructure of water/ethanol mixtures, essentially decreasing the level of larger local clustering. Overall, the work presents a systematic characterization of ethanol effects on a widely used, common protein type, spidroins, and generalizes the findings to other intrinsically disordered proteins by pinpointing the general features of the response. The results can aid in designing effective alcohol treatments for proteins, but also enable design and tuning of protein material properties by a relatively controllable solvation handle, the addition of ethanol.

Effects of solvent phase recycling on microalgae liquefaction in ethanol: Bio-oil production and nitrogen transformation

Liquefaction of microalgae in ethanol offers an eco-friendly bio-oil alternative, but solvent recycling is crucial for sustainability due to extra costs. In this work, Chlorella vulgaris was liquefied in supercritical ethanol at 260 °C, and the solvent phase (SP) separated from bio-oil was recovered and reused. Five liquefaction cycles were performed at identical temperature and pressure conditions to investigate the effects on oil production and nitrogen transformation. The findings demonstrated a gradual increase in water content in recycled SP. Ethanol-water co-solvent as the reaction medium promoted the decomposition and re-polymerization of protein in raw material, thus increasing the bio-oil yield (76.84 %) and higher heating value (33.53 MJ/kg) to some extent. Simultaneously, the relative nitrogen content of bio-oil rose from 8.03 % to 8.52 %, predominantly in the form of nitrogen heterocycles. The potential pathway for nitrogen conversion was revealed, which establishes a theoretical basis for the subsequent denitrification of bio-oil.

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.

Radial distribution and hydrogen bonded network graphs of alcohol-aniline binary mixture

CONTEXT

Hydrogen bonds play a vital role in the stability and functioning of biomolecules. Suitable binary liquids are often used as prototypes for the study of biologically significant hydrogen bond studies and their intricate networks. Often, such systems show deviations in their physico-chemical properties from ideal conditions. As a continuation of our research interest in biologically important hydrogen-bonded systems, this paper reports the classical molecular dynamic studies on mixtures of aniline with 8 primary alcohols (C_RH_2R+1-OH, R = 1 to 8) for the complete concentration range. The energetics results indicate the predominance of OH--O interactions over other hydrogen bonds. Structures in the network are analyzed using radial distribution function (RDF), hydrogen bond statistics, and graph theoretical analysis (GTA). Coordination numbers, hydrogen bond statistics, and GTA show a bunching of alcohol-alcohol hydrogen bonds for lower aniline concentrations, while the aniline-aniline interactions are not affected by changes in the concentration.

METHODS

Interaction energies are calculated using B3LYP/6-311G++(d, p) density functional theory using Gaussian-09. The molecular dynamics simulations are carried out using GROMACS (V 2020.6) with the OPLS/AA force field and the simulation box is visualized using VMD. The NetworkX Python package is used for GTA calculation.

An atomistic explanation of the ethanol-water azeotrope

Ethanol and water form an azeotropic mixture at an ethanol molecular percentage of ∼91% (∼96% by volume), which prohibits ethanol from being further purified via distillation. Aqueous solutions at different concentrations in ethanol have been studied both experimentally and theoretically. We performed cylindrical micro-jet photoelectron spectroscopy, excited by synchrotron radiation, 70 eV above C1s ionization threshold, providing optimal atomic-scale surface-probing. Large model systems have been employed to simulate, by molecular dynamics, slabs of the aqueous solutions and obtain an atomistic description of both bulk and surface regions. We show how the azeotropic behaviour results from an unexpected concentration-dependence of the surface composition. While ethanol strongly dominates the surface and water is almost completely depleted from the surface for most mixing ratios, the different intermolecular bonding patterns of the two components cause water to penetrate to the surface region at high ethanol concentrations. The addition of surface water increases its relative vapour pressure, giving rise to the azeotropic behaviour.

Geometries of hydrogen bonds in water-ethanol mixtures from ab initio molecular dynamics simulations

We outline a simple procedure to determine the geometry of hydrogen bonds between different molecular species in binary mixtures from ab initio molecular dynamics (AIMD) trajectories. Here we determine the geometry of the hydrogen bonds arising from intermolecular OH⋯O interactions between different H-bonded pairs, water-water, ethanol-ethanol and water-ethanol in water-alcohol mixtures at different compositions by plotting the intermolecular non-bonded OH⋯O and O⋯O distances, and the ∠HO⋯O (θ) angles for each of the possible pairs in the ensemble. Two regions separate out in each of the scatter-plots; the one with short OH⋯O and O⋯O intermolecular distances and almost linear ∠HO⋯O angles may be identified as the region where the intermolecular OH⋯O geometry would be favorable for hydrogen bonding. Using the different geometric criteria for each of the three possible H-bonded pairs we estimate the average number of water and ethanol molecules that are hydrogen bonded to a water molecule, and to an ethanol molecule, respectively, at different mole fractions of the mixture. We validate the results from values of the chemical shift of the two OH resonances (water and ethanol) in the proton NMR spectra of the mixtures at different concentrations as these values are known to be sensitive to the local chemical environment of the resonating nuclei.

Optical Fiber Sensor for Monitoring the Evaporation of Ethanol–Water Mixtures

An inline optical fiber sensor is proposed to monitor in real time the evaporation process of ethanol–water binary mixtures. The sensor presents two interferometers, a cladding modal interferometer (CMI) and a Mach–Zehnder interferometer (MZI). The CMI is used to acquire the variations in the external medium refractive index, presenting a maximum sensitivity of 387 nm/RIU, and to attain the variation in the sample concentration profile, while the MZI monitors temperature fluctuations. For comparison purposes, an image analysis is also conducted to obtain the droplet profile. The sensor proposed in this work is a promising alternative in applications where a rigorous measurement of volatile organic compound concentrations is required, and in the study of chemical and physical properties related to the evaporation process.

Hydrophobic Cluster Formation in Aqueous Ethanol Solutions Probed by Soft X-ray Absorption Spectroscopy

Hydrophobic cluster structures in aqueous ethanol solutions at different concentrations have been investigated by soft X-ray absorption spectroscopy (XAS). In the O K-edge XAS, we have found that hydrogen bond structures among water molecules are enhanced in the middle-concentration region by the hydrophobic interaction of the ethyl groups in ethanol. In the C K-edge XAS, the lower energy features arise from a transition from the terminal methyl C 1s electron to an unoccupied orbital of 3s Rydberg character, which is sensitive to the nearest-neighbor intermolecular interactions. From the comparison of C K-edge XAS with the inner-shell calculations, we have found that ethanol clusters are easily formed in the middle-concentration region due to the hydrophobic interaction of the ethyl group in ethanol, resulting in the enhancement of the hydrogen bond structures among water molecules. This behavior is different from aqueous methanol solutions, where the methanol-water mixed clusters are more predominant in the middle-concentration region due to the relatively weak hydrophobic interactions of the methyl group in methanol.

Insights on choline chloride-based deep eutectic solvent (reline) + primary alcohol mixtures: a molecular dynamics simulation study

Deep eutectic solvents (DESs) emerged as green solvents for new generation technologies owing to their high chemical and thermal stability. Addition of restricted amount of organic solvents into the DESs plays a significant role in the improvement of thermodynamic and the transport properties to work as a potential solvent in process industries. In this paper, molecular dynamics (MD) simulations were performed to understand the thermophysical and transport properties of choline chloride-based DES (reline) and primary alcohol (methanol and ethanol) mixture in relation to microscopic structure. Density, radial distribution function, coordination number, average number of H-bond, diffusion coefficient and spatial distribution function was calculated in order to understand the structure and involvement of H-bond network at an atomic level. H-bond and spatial distribution function analyses revealed that the chloride ion prefers to be spatially distributed around hydroxyl group of alcohol and found to be more pronounced upon increase in alcohol concentration. As a consequence, it was observed that the H-bonds between Cl^- and urea decreases overall with the loading of alcohol and effect is more pronounced beyond a concentration of 0.4. Self-diffusion values for choline, Cl^- and urea were found to be increased significantly upon increase in concentration of alcohol from 0.6 to 0.8. Overall, our simulation points to the interplay and interactions between the chloride ions and the solvents in determining the structural and transport properties of choline chloride-based DES.

Rigid Cations Induce Enhancement of Microheterogeneity and Exhibit Anomalous Ion Diffusion in Water-Ethanol Mixtures

Because of the amphiphilic nature of ethanol in the aqueous solution, ions cause an interesting microheterogeneity where the water molecules and the hydroxy groups of ethanol preferentially solvate the ions, while the ethyl groups tend to occupy the intervening space. Using computer simulations, we study the dynamics of rigid monovalent cations (Li⁺, Na⁺, K⁺, and Cs⁺) in aqueous ethanol solutions with chloride as the counterion. We vary both the size of the ions and the composition of the mixture to explore size- and composition-dependent ion diffusion. The relative stability of enhanced microheterogeneous configurations makes ion diffusion slower than what would be surmised by using the bulk properties of the mixture, using the Stokes-Einstein relation. We study the structure through partial radial distribution functions and the stability through coordination number fluctuations. The ion diffusion coefficient exhibits sharp re-entrant behavior when plotted against viscosity varied by composition. Our studies reveal multiple anomalous features of ion motion in this mixture. We formulate a mode-coupling theory (MCT) that takes into account the interaction between different dynamical components; MCT can incorporate the effects of heterogeneous dynamics and nonlinearity in composition dependence that arise from the feedback between mutually dependent ion-solvent dynamics.

Correlating solvation with conformational pathways of proteins in alcohol-water mixtures: a THz spectroscopic insight

Water, being an active participant in most of the biophysical processes, is important to trace how protein solvation changes as its conformation evolves in the presence of solutes or co-solvents. In this study, we investigate how the secondary structures of two diverse proteins - lysozyme and β-lactoglobulin - change in the aqueous mixtures of two alcohols - ethanol and 2,2,2-trifluoroethanol (TFE) using circular dichroism measurements. We observe that these alcohols change the secondary structures of these proteins and the changes are protein-specific. Subsequently, we measure the collective solvation dynamics of these two proteins both in the absence and in the presence of alcohols by measuring the frequency-dependent absorption coefficient (α(ν)) in the THz (0.1-1.2 THz) frequency domain. The alcohol-water mixtures exhibit a non-ideal behaviour with the highest absorption difference (Δα) obtained at X_alcohol = 0.2. The protein solvation in the presence of the alcohols shows an oscillating behaviour in which Δα_protein changes with X_alcohol. Such an oscillatory behaviour of protein solvation results from a delicate interplay between the protein-water, protein-alcohol and water-alcohol associations. We attempt to correlate the various structural conformations of the proteins with the associated solvation.

H2O- and ethanol concentration-responsive polymer/gel inverse opal photonic crystal

Responsive photonic crystals have attracted much attention due to their strong capability to manipulate the propagation of light in the visible region, but it is still a big challenge to invisibility and mechanical stability. Here, the novel Poly(ether sulfone)/Poly(acrylic acid) inverse opal photonic crystals, which have high mechanical stability and can release visible patterns after wetting with water, are discussed. The Poly(ether sulfone)/Poly(acrylic acid) inverse opal photonic crystals are also responsive to the concentration of ethanol, and the structural color response times increase with increasing ethanol concentration. This design uses the selective infiltration, hydrogen bonding and capillary action of solvent to realize the spectral diversity of reflectance. Owing to the high polarity and hydrogen bonding ability of carboxyl groups, water molecules are adsorbed easily by the poly(acrylic acid) gel. Subsequently, the encrypted information is decrypted due to the redshift of the structural color. Because of its lower polarity and hydrogen bonding ability relative to water, ethanol can impede the absorption of solvent by gel. Therefore, the ethanol concentration can be identified based on the structural color response time. Furthermore, reliable information decryption methods make Poly(ether sulfone)/Poly(acrylic acid) inverse opal photonic crystals potentially uesful as trusted encryption devices.

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.

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.

The Structures of ZnCl2-Ethanol Mixtures, a Spectroscopic and Quantum Chemical Calculation Study

We report in this article the structural properties, spectral behavior and heterogeneity of ZnCl₂-ethanol (EtOH) mixtures in a wide-composition range (1:3 to 1:14 in molar ratios), using ATR-FTIR spectroscopy and quantum chemical calculations. To improve the resolution of the initial IR spectra, excess spectroscopy and two-dimensional correlation spectroscopy were employed. The transformation process was suggested to be from EtOH trimer and EtOH tetramer to EtOH monomer, EtOH dimer and ZnCl₂-3EtOH complex upon mixing. The theoretical findings showed that increasing the content of EtOH was accompanied with the flow of negative charge to ZnCl₂. This led to reinforcement of the Zn←O coordination bonds, increase of the ionic character of Zn‒Cl bond and weakening and even dissociation of the Zn‒Cl bond. It was found that in some of the ZnCl₂-EtOH complexes optimized at the gas phase or under the solvent effect, there existed hydroxyls with a very special interactive array in the form of Cl‒Zn⁺←O‒H^…Cl^-, which incredibly red-shifted to wavenumbers <3000 cm^-1. This in-depth study shows the physical insights of the respective electrolyte alcoholic solutions, particularly the solvation process of the salt, help to rationalize the reported experimental results, and may shed light on understanding the properties of the deep eutectic solvents formed from ZnCl₂ and an alcohol.

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.

Molecular Origin of the Prepeak in the Structure Factor of Alcohols

Prepeak in the structure factor of alcohols is known for a half century and was attributed to one of two mechanisms (i) self-assembly in aggregates and (ii) existence of spatial heterogeneity. Although both explnations are often argued the molecular origin is yet unclear. In this work, molecular dynamics simulation of neat alcohols and their mixtures in the presence of an apolar liquid in bulk and in confined phases is performed to unveil and to clarify the origin of the prepeak at the molecular scale. Unambiguously, we show that the existence of the prepeak is the result of the self-assembly in clusters leading to long-range correlations rather than the spatial heterogeneity. We also establish that the confinement of neat liquids at the nanoscale does not erase the clustering and the prepeak but strongly reduce the spatial heterogeneity. Regarding the binary alcohol/toluene mixtures, we highlight the possibility to erase the clustering and the spatial heterogeneity from nanoconfinement inducing the formation of a core-shell structure. By tuning the interfacial chemistry and pore size, we shed light on the possibility to control the spatial heterogeneity, the self-assembly, and the microphase separation.

Molecular modelling techniques for predicting liquid-liquid interfacial properties of methanol plus alkane (n-hexane, n-heptane, n-octane) mixtures

In this work, the liquid-liquid interfacial properties of methanol plus n-alkane (n-hexane, n-heptane, n-octane) mixtures are investigated at atmospheric pressure by two complementary molecular modelling techniques; namely, molecular dynamic simulations (MD) and density gradient theory (DGT) coupled with the PC-SAFT (perturbed-chain statistical associating fluid theory) equation of state. Furthermore, two molecular models of methanol are used, which are based on a non-polarisable three site approach. On the one hand, is the original (flexible) TraPPE-UA model force field. On the other hand, is the rigid approximation denoted as OPLS/2016. In both cases, n-alkanes are modelled using the TraPPE-UA model. Simulations are performed using the direct coexistence technique in the ensemble. Special attention is paid to the comparison between the estimations obtained from different methanol models, the available experimental data and theoretical calculations. In all cases, the rigid model is capable of predicting the experimental phase equilibrium and interfacial properties accurately. Unsurprisingly, the methanol-rich density and interfacial tension are overestimated using the TraPPE model combined with Lorentz-Berthelot mixing rules for predicting the mixture behaviour. Accurate comparison between MD and DGT plus PC-SAFT requires consideration of the cross-interactions between individual influence parameters and fitting the βij values. This latter aspect is particularly important because it allows the exploitation of the link between the EOS model and the direct molecular simulation of the corresponding fluid. At the same time, it was demonstrated that the key property defining the interfacial tension value is the absolute concentration of methanol in the methanol-rich phase. This behaviour indicates that there are more hydrogens bonded with each other, and they interact favourably with an increasing number of carbon atoms in the alkane.

Force-Field Simulations of a Hydrated Lanreotide-Based Derivative: Hydration, Dynamics, and Numerical Evidence of Self-Assembly in Dimers

Recently, self-organization of the cyclic octapeptide lanreotide and lanreotide-based derivatives in a nanotube to from a dimer structure has been experimentally evidenced. While the nature of the interactions between both monomers has been strongly investigated no molecular details of the hydration of the monomer and the formation of the dimer have been provided. Using molecular dynamics simulations, this work focuses on the structure, hydration, and dynamics of water and an analog of lanreotide. To do so, several models of monomers based on different schemes of partial charges and electrostatic interaction calculations are considered. By comparison with the experiments, we show that the model based on the combination of the AMBER force-field, CHELPG charge calculation, Ewald sum is the most relevant. Additionally, by mapping the interfacial hydration of the lanreotide monomer we evidence a heterogeneous surface in terms of hydrophilicity involving heterogeneous hydration. Furthermore, we show a slowdown in the translational dynamics of water molecules located close to the lanreotide surface. We also provide the molecular details of the self-assembly in the dimer in terms of structure, hydration, and energy.

Looking for the Azeotrope: A Computational Study of (Ethanol)6-Water, (Methanol)6-Water, (Ethanol)7, and (Methanol)7 Heptamers

Considering that a molecular-level understanding of the azeotropic ethanol-water system can contribute to the search of new methodologies and/or modifications of industrial separation methods, this study tries to provide some clues to understand why azeotropes should be expected for ethanol, but not for methanol. Our exploration of the potential energy surface of (ethanol)₆-water heteroheptamers, carried out at the B3LYP-D3/6-311++G(d,p) level, shows these heteroclusters to exhibit a cyclic structure where the cooperativity effects between the OH···O HBs is a fundamental ingredient. An analysis of this cooperativity clearly indicates that ethanol-water systems will exhibit a similarly high stability as the heterocluster size approaches the azeotrope. However, a similar behavior should not be expected for the methanol-containing analogues. A comparison between (ethanol)₇, (ethanol)₆-water, (methanol)₇, and (methanol)₆-water shows the ethanol-containing systems to be significantly more stable than the methanol-containing analogues. This result is probably due to the fact that the OH···O HBs are weaker than those found between ethanol molecules. However, our atoms in molecule (AIM) and noncovalent interaction (NCI) analyses unambiguously show that important contributors to the enhanced stability of the ethanol-containing clusters are the secondary van der Waals interactions between ethyl groups, which are not observed between methyl groups. Hence, while the formation of stable azeotropes is expected for the case of ethanol, for the methanol-containing analogues, the relative stability of the clusters is significantly smaller, and its formation is accompanied by an increase of the free energy.

Understanding alcohol aggregates and the water hydrogen bond network towards miscibility in alcohol solutions: graph theoretical analysis

Under ambient conditions, methanol and ethanol are miscible in water at all concentrations, while n-butanol is partially miscible. This is the first study to quantitatively examine the miscibility of butanol and compare with miscible alcohols by employing molecular dynamics simulations and graph theoretical analysis of three water-alcohol mixtures at various concentrations. We show how distinct alcohol aggregates are formed, thereby affecting the water structure, which established the relationship between the morphological structure of the aggregates and the miscibility of the alcohol in aqueous solution. The aggregates of methanol and ethanol in highly concentrated solutions form an extended H-bond network that intertwines well with the H-bond network of water. n-Butanol tends to self-associate and form large aggregates, while such aggregates are segregated from water. Graph theoretical analysis revealed that the alcohol aggregates of methanol and ethanol solutions have a morphological structure different from that of n-butanol, although there is no significant difference in morphology between the three pure alcohols. These two distinct alcohol aggregates are classified as water-compatible and water-incompatible depending upon their interaction with the water H-bond network, and their effect on the water structure was investigated. Our study reveals that the water-compatible network of alcohol aggregates in methanol and ethanol solutions disrupts the water H-bond networks, while the water-incompatible network of n-butanol aggregates does not considerably alter the water structure, which is consistent with the experimental results. Furthermore, we propose that miscible alcohols form water-compatible networks in binary aqueous systems while partially miscible alcohols form water-incompatible networks. The bifurcating hypothesis on the alcohol aggregation behavior in liquid water is of critical use to understand the fundamental issues such as solubility and phase separation in solution systems.

Decelerated Liquid Dynamics Induced by Component-Dependent Supercooling in Hydrogen and Deuterium Quantum Mixtures

Isotopic mixtures of p-H₂ and o-D₂ molecules have been an attractive binary system because they include two kinds of purely isotopic molecules which possess the same electronic potential but the twice different mass inducing differently pronounced nuclear quantum effects (NQEs). Accessing details of structures and dynamics in such quantum mixtures combining complex molecular dynamics with NQEs of different strengths remains a challenging problem. Taking advantage of the nonempirical molecular dynamics method which describes p-H₂ and o-D₂ molecules, we found that the liquid dynamics slows down at a specific mixing ratio, which can be connected to the observed anomalous slowdown of crystallization in the quantum mixtures. We attributed the decelerated dynamics to the component-dependent supercooling of p-H₂ taking place in the mixtures, demonstrating that there is an optimal mixing ratio to hinder crystallization. The obtained physical insights will help in experimentally controlling and achieving unknown quantum mixtures including superfluid.

Hydrogen-Bonding Interactions in Hybrid Aqueous/Nonaqueous Electrolytes Enable Low-Cost and Long-Lifespan Sodium-Ion Storage

Although "water-in-salt" electrolytes have opened a new pathway to expand the electrochemical stability window of aqueous electrolytes, the electrode instability and irreversible proton co-insertion caused by aqueous media still hinder the practical application, even when using exotic fluorinated salts. In this study, an accessible hybrid electrolyte class based on common sodium salts is proposed, and crucially an ethanol-rich media is introduced to achieve highly stable Na-ion electrochemistry. Here, ethanol exerts a strong hydrogen-bonding effect on water, simultaneously expanding the electrochemical stability window of the hybridized electrolyte to 2.5 V, restricting degradation activities, reducing transition metal dissolution from the cathode material, and improving electrolyte-electrode wettability. The binary ethanol-water solvent enables the impressive cycling of sodium-ion batteries based on perchlorate, chloride, and acetate electrolyte salts. Notably, a Na_0.44MnO₂ electrode exhibits both high capacity (81 mAh g^-1) and a remarkably long cycle life >1000 cycles at 100 mA g^-1 (a capacity decay rate per cycle of 0.024%) in a 1 M sodium acetate system. The Na_0.44MnO₂/Zn full cells also show excellent cycling stability and rate capability in a wide temperature range. The gained understanding of the hydrogen-bonding interactions in the hybridized electrolyte can provide new battery chemistry guidelines in designing promising candidates for developing low-cost and long-lifespan batteries based on other (Li⁺, K⁺, Zn^2+, Mg²⁺, and Al³⁺) systems.

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.

Molecular Mechanism for Azeotrope Formation in Ethanol/Benzene Binary Mixtures through Gibbs Ensemble Monte Carlo Simulation

Azeotropes have been studied for decades due to the challenges they impose on separation processes but fundamental understanding at the molecular level remains limited. Although molecular simulation has demonstrated its capability of predicting mixture vapor-liquid equilibrium (VLE) behaviors, including azeotropes, its potential for mechanistic investigation has not been fully exploited. In this study, we use the united atom transferable potentials for phase equilibria (TraPPE-UA) force field to model the ethanol/benzene mixture, which displays a positive azeotrope. Gibbs ensemble Monte Carlo (GEMC) simulation is performed to predict the VLE phase diagram, including an azeotrope point. The results accurately agree with experimental measurements. We argue that the molecular mechanism of azeotrope formation cannot be fully understood by studying the mixture liquid-state stability at the azeotrope point alone. Rather, azeotrope occurrence is only a reflection of the changing relative volatility between the two components over a much wider composition range. A thermodynamic criterion is thus proposed on the basis of the comparison of partial excess Gibbs energy between the components. In the ethanol/benzene system, molecular energetics shows that with increasing ethanol mole fraction, its volatility initially decreases but later plateaus, while benzene volatility is initially nearly constant and only starts to decrease when its mole fraction is low. Analysis of the mixture liquid structure, including a detailed investigation of ethanol hydrogen-bonding configurations at different composition levels, reveals the underlying molecular mechanism for the changing volatilities responsible for the azeotrope.

Self-assembly of small molecules at hydrophobic interfaces using group effect

Although common in nature, the self-assembly of small molecules at sold-liquid interfaces is difficult to control in artificial systems. The high mobility of dissolved small molecules limits their residence at the interface, typically restricting the self-assembly to systems under confinement or with mobile tethers between the molecules and the surface. Small hydrogen-bonding molecules can overcome these issues by exploiting group-effect stabilization to achieve non-tethered self-assembly at hydrophobic interfaces. Significantly, the weak molecular interactions with the solid makes it possible to influence the interfacial hydrogen bond network, potentially creating a wide variety of supramolecular structures. Here we investigate the nanoscale details of water and alcohols mixtures self-assembling at the interface with graphite through group-effect. We explore the interplay between inter-molecular and surface interactions by adding small amounts of foreign molecules able to interfere with the hydrogen bond network and systematically varying the length of the alcohol hydrocarbon chain. The resulting supramolecular structures forming at room temperature are then examined using atomic force microscopy with insights from computer simulations. We show that the group-based self-assembly approach investigated here is general and can be reproduced on other substrates such as molybdenum disulphide and graphene oxide, potentially making it relevant for a wide variety of systems.

Anomalous dynamics of water at the octopeptide lanreotide surface

This work reports the study of water dynamics close to the cyclic octapeptide lanreotide from atomistic simulations of hydrated lanreotide, a cyclic octapeptide. Calculation of the hydrogen bonds between water molecules allows mapping of the hydrophilic regions of lanreotide. Whereas a super-diffusivity of the interfacial water molecules is established, a slowdown in rotational dynamics is observed, involving a decoupling between both processes. Acceleration in translation dynamics is connected to the hopping process between hydrophilic zones. Microscopically, this is correlated with the weakness of the interfacial hydrogen bonding network due to a hydrophobic interface at the origin of the interfacial sliding of water molecules. Heterogeneous rotational dynamics of water molecules close the lanreotide surface is evidenced and connected to heterogeneous hydration.

Molecular dynamics simulations of a hydrated mutated lanreotide, a cyclic octapeptide, were carried out to characterize its hydration state. We studied the water dynamics close to the peptide using atomistic simulations.

Diffusion in Binary Aqueous Solutions of Alcohols by Molecular Simulation

Based on the molecular dynamics method, the calculations for diffusion coefficients were carried out in binary aqueous solutions of three alcohols: ethanol, isopropanol, and tert-butanol. The intermolecular potential TIP4P/2005 was used for water; and five force fields were analyzed for the alcohols. The force fields providing the best accuracy of calculation were identified based on a comparison of the calculated self-diffusion coefficients of pure alcohols with the experimental data for internal (Einstein) diffusion coefficients of alcohols in solutions. The temperature and concentration dependences of the interdiffusion coefficients were determined using Darken’s Equation. Transport (Fickian) diffusion coefficients were calculated using a thermodynamic factor determined by the non-random two-liquid (NRTL) and Willson models. It was demonstrated that for adequate reproduction of the experimental data when calculating the transport diffusion coefficients, the thermodynamic factor has to be 0.64. Simple approximations were obtained, providing satisfactory accuracy in calculating the concentration and temperature dependences of the transport diffusion coefficients in the studied mixtures.

Microstructure of nonideal methanol binary liquid mixtures

The nonideality of binary mixtures is often related to the nature of the interactions between both liquids and of the heterogeneity at the nanoscale-named microstructure. When one of the liquids is a hydrogen bonds former and the second is aprotic, the progressive diluting of the hydrogen-bonding network leads to a clustering and nanophases. By considering two mixtures, toluene-methanol and cyclohexane-methanol, the nonideality and its connection with the structure at the nanoscale and the intermolecular interactions are numerically investigated. Contrary to the toluene that is fully miscible in methanol, cyclohexane presents a high range of immiscibility which makes it a relevant system to study the nucleation (local segregation) and its propagation. In both mixtures, the deviation from the ideal behavior is observed. In the case of the toluene-methanol mixture, the initial hydrogen-bonding network corresponding to a homogenous structure is locally broken due to the favorable toluene-methanol interactions leading to the spatial heterogeneity at the origin of the nonideality. In the range of miscibility of the cyclohexane-methanol mixtures, the formation of hydrophobic nanophases of larger size is observed due to the unfavorable interactions between both components leading to a self-organizing of cyclohexane molecules. The immiscibility of cyclohexane and methanol are then correlated to the formation of nanophases and their propagation, which are also at the origin of the spatial heterogeneity. In the pure methanol, we highlight the disconnection between the clustering and the heterogeneity. We shed light on the fact that the prepeak observed in the structure factor is independent of the degree of heterogeneity, but is connected to the presence of cyclic clusters.

Highly Selective Supported Graphene Oxide Membranes for Water-Ethanol Separation

A polyethersulfone (PES)-supported graphene oxide (GO) membrane has been developed by a simple casting approach. This stable membrane is applied for ethanol/water separation at different temperatures. The 5.0 µm thick GO film coated on PES support membrane showed a long-term stability over a testing period of one month and excellent water/ethanol selectivity at elevated temperatures. The water/ethanol selectivity is dependent on ethanol weight percentage in water/ethanol feed mixtures and on operating temperature. The water/ethanol selectivity was enhanced with an increase of ethanol weight percentage in water/ethanol mixtures, from below 100 at RT to close to 874 at a 90 °C for 90% ethanol/10% water mixture. Molecular dynamics simulation of water-ethanol mixtures in graphene bilayers, that are considered to play a key role in transport, revealed that molecular transport is negligible for layer spacing below 1 nm. The differences in the diffusion of ethanol and water in the bilayer are not consistent with the large selectivity value experimentally observed. The entry of water and ethanol into the interlayer space may be the crucial step controlling the selectivity.

Water fluctuation in methanol, ethanol, and 1-propanol aqueous-mixture probed by Brownian motion

The origin of the driving force in Brownian motion is the collision between the colloidal particle and the molecules of the surrounding fluid. Therefore, Brownian motion contains information on the local solvent structures of the surrounding colloid. The mean square displacement in a water-ethanol mixture is greater than that anticipated from the macroscopic shear viscosity, indicating that the microscopic movement of Brownian motion involves the local information on the water-ethanol mixture on a molecular level, i.e., an inhomogeneity in the Brownian particle size (∼1 μm). Here, the Brownian motion of mixtures of water and methanol, ethanol, and 1-propanol are systematically investigated. Similar discrepancies between the microscopic and macroscopic viscosities are observed at low alcohol molar concentrations, for all the alcohol mixtures. This means that inhomogeneity with water fluctuation is important in explanation of the unusual Brownian diffusions of alcohol aqueous solutions. The Brownian motion also reveals a thermal energy conversion mechanism between translation and rotation.