Kirkwood-Buff integrals of aqueous alcohol binary mixtures

Kirkwood-Buff Analysis of Binary and Ternary Systems Consisting of Alcohols (Methanol, Ethanol, 1-Propanol, and 2-Propanol), Water, and n-Hexane to Understand the Formation of Surfactant-Free Microemulsions

Surfactant-free microemulsion (SFME) represents a class of fluid mixtures that can form microheterogeneous structures without detergents, offering an environmentally benign alternative to traditional microemulsions. However, the formation mechanism is still elusive. This work applies the Kirkwood-Buff theory to mixtures of alcohols, water, and n-hexane to elucidate the SFME formation mechanism. To ensure robust calculation of the Kirkwood-Buff integrals (KBIs), we construct a data set of densities and excess free energies of binary and ternary systems. Multiple excess Gibbs free energy models are assessed against this data set to select the most suitable model reproducing the experimental results. In addition, we introduce statistical methods to determine the optimal polynomial order of the Redlich-Kister correlation for the excess volume data. We first validate our methodology in binary systems. Then, we extend the calculation method to ternary mixtures. The KBI calculation results reveal that the alcohol-hexane and water-hexane interactions do not significantly affect SFME formation. In contrast, the interplay among water-water, water-alcohol, and alcohol-alcohol interactions critically influences the ability of a liquid mixture to form SFME structures. SFME systems exhibit the facile formation of water aggregates enveloped by alcohols, whereas non-SFME systems demonstrate homogeneous alcohol/water droplets dispersed in an oil continuous medium.

Concentration Fluctuation/Microheterogeneity Duality Illustrated with Aqueous 1,4-Dioxane Mixtures

The structural properties of aqueous 1-4 dioxane mixtures are studied by computer simulations of different water and dioxane force field models, from the perspective of illustrating the link between structural properties at the molecular level and measurable properties such as radiation scattering intensities and Kirkwood-Buff integrals (KBIs). A strategy to consistently correct the KBI obtained from simulations is proposed, which allows us to obtain the genuine KBI corresponding to a given pair of molecular species, in the entire concentration range, and without necessitating excessively large system sizes. The application of this method to the aqueous dioxane mixtures, with an all-atom CHARMM dioxane model and 2 water models, namely, SPC/E and TIP3P, allows one to understand the differences in the structure of the corresponding mixtures at the molecular level, particularly concerning the role of the water aggregates and its model dependence. This study allows us to characterize the dual role played by the concentration fluctuations and the domain segregation, particularly in what concerns the calculated X-ray spectra.

Smart Polymers for Soft Materials: From Solution Processing to Organic Solids

Polymeric materials are ubiquitous in our everyday life, where they find a broad range of uses-spanning across common household items to advanced materials for modern technologies. In the context of the latter, so called "smart polymers" have received a lot of attention. These systems are soluble in water below their lower critical solution temperature Tℓ and often exhibit counterintuitive solvation behavior in mixed solvents. A polymer is known as smart-responsive when a slight change in external stimuli can significantly change its structure, functionm and stability. The interplay of different interactions, especially hydrogen bonds, can also be used for the design of lightweight high-performance organic solids with tunable properties. Here, a general scheme for establishing a structure-property relationship is a challenge using the conventional simulation techniques and also in standard experiments. From the theoretical side, a broad range of all-atom, multiscale, generic, and analytical techniques have been developed linking monomer level interaction details with macroscopic material properties. In this review, we briefly summarize the recent developments in the field of smart polymers, together with complementary experiments. For this purpose, we will specifically discuss the following: (1) the solution processing of responsive polymers and (2) their use in organic solids, with a goal to provide a microscopic understanding that may be used as a guiding tool for future experiments and/or simulations regarding designing advanced functional materials.

Camel back shaped Kirkwood-Buff Integrals

Some binary mixtures, such as specific alcohol-alkane mixtures, or even water-tbutanol, exhibit two humps "camel back" shaped KBI. This is in sharp contrast with usual KBI of binary mixtures having a single extremum. This extremum is interpreted as the region of maximum concentration fluctuations, and usually occurs in binary mixtures presenting appreciable micro-segregation, and corresponds to where the mixture exhibit a percolation of the two species domains. In this paper, it is shown that two extrema occur in binary mixtures when one species forms "meta-particle" aggregates, the latter which act as a meta-species, and have their own concentration fluctuations, hence their own KBI extremum. This "meta-extremum" occurs at low concentration of the aggregate-forming species (such as alcohol in alkane), and is independant of the other usual extremum observed at mid volume fraction occupancy. These systems are a good illustration of the concept of the duality between concentration fluctuations and micro-segregation.

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.

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.

On the Solute-Induced Structure-Making/Breaking Effect: Rigorous Links among Microscopic Behavior, Solvation Properties, and Solution Non-Ideality

We studied the solute-induced perturbation of the solvent environment around a solute species from a microscopic viewpoint and propose a novel approach to the understanding of the structure-making/breaking process, regardless of the type and nature of the solute-solvent interactions. Based on the Kirkwood-Buff fluctuation formalism, we present a rigorous statistical mechanics description of the evolution of the solvent structure around the solute, analyze its response to small perturbations of the ( TP) state conditions and composition of the system, and make direct connections between a few equivalent micro- and macroscopic manifestations as probes for, and targets of, experimental measurements. We illustrate the analysis with theoretical results from integral equation calculations of model fluids and experimental evidence from available data for a variety of aqueous electrolyte and nonelectrolyte real fluid solutions. Finally, we provide a critical discussion about the inadequacy underlying a widely used de facto criterion for the classification of structure-making/breaking solutes.

Thermodynamic Analysis of n-Hexane-Ethanol Binary Mixtures Using the Kirkwood-Buff Theory

A complete thermodynamic analysis of mixtures consisting of molecules with complex chemical constitution can be rather demanding. The Kirkwood-Buff theory of solutions allows the estimation of thermodynamic properties, which cannot be directly extracted from atomistic simulations, such as the Gibbs energy of mixing (Δ_mix G). In this work, we perform molecular dynamics simulations of n-hexane-ethanol binary mixtures in the liquid state under two temperature-pressure conditions and at various mole fractions. On the basis of the recently published methodology of Galata [ Fluid Phase Equilib. 2018 , 470 , 25 - 37 ] , we first calculate the Kirkwood-Buff integrals in the isothermal-isobaric ( NpT) ensemble, identifying how system size affects their estimation. We then extract the activity coefficients, excess Gibbs energy, excess enthalpy, and excess entropy for the n-hexane-ethanol binary mixtures we simulate. We employ two approaches for quantifying composition fluctuations: one based on counting molecular centers of mass and a second one based on counting molecular segments. Results from the two approaches are practically indistinguishable. We compare our results against predictions of vapor-liquid equilibria obtained in a previous simulation work using the same force field, as well as with experimental data, and find very good agreement. In addition, we develop a simple methodology to identify the hydrogen bonds between ethanol molecules and analyze their effects on mixing properties.

Time-Resolved Structured Illumination Microscopy for Phase Separation Dynamics of Water and 2-Butoxyethanol Mixtures: Interpretation of "Early Stage" Involving Micelle-Like Structures

Phase separation dynamics of a water/2-butoxyethanol (2BE) mixture was studied with newly developed time-resolved structured illumination microscopy (SIM). Interestingly, an employed hydrophobic fluorescent probe for SIM showed spectral shifts up to 500 ns after a laser-induced temperature jump, which suggests 2BE micellar-like aggregates become more hydrophobic at the initial stage of phase separation. This hydrophobic environment in 2BE aggregates, probably due to the ejection of water molecules, continued up to at least 10 μs. Time-resolved SIM and previously reported light scattering data clearly showed that the size of a periodic structure remained constant (ca. 300 nm) from 3 to 10 μs, and then the growth of periodic structures having the self-similarity started. We think that the former and the latter processes correspond to "early stage" (concentration growth) and "late stage" (size growth), respectively, in phase separation dynamics. Here we suggest that, in the early stage, the entity to bear 2BE phase be water-poor 2BE aggregates, and the number density of these aggregates would simply increase in time.

Depleted depletion drives polymer swelling in poor solvent mixtures

Establishing a link between macromolecular conformation and microscopic interaction is a key to understand properties of polymer solutions and for designing technologically relevant "smart" polymers. Here, polymer solvation in solvent mixtures strike as paradoxical phenomena. For example, when adding polymers to a solvent, such that all particle interactions are repulsive, polymer chains can collapse due to increased monomer-solvent repulsion. This depletion induced monomer-monomer attraction is well known from colloidal stability. A typical example is poly(methyl methacrylate) (PMMA) in water or small alcohols. While polymer collapse in a single poor solvent is well understood, the observed polymer swelling in mixtures of two repulsive solvents is surprising. By combining simulations and theoretical concepts known from polymer physics and colloidal science, we unveil the microscopic, generic origin of this collapse-swelling-collapse behavior. We show that this phenomenon naturally emerges at constant pressure when an appropriate balance of entropically driven depletion interactions is achieved.

Kirkwood-Buff integrals for hard-core Yukawa fluids

The Kirkwood-Buff (KB) theory of solution is employed to investigate several macroscopic properties of the one-component hard-core Yukawa (HCY) fluid, where the key physical quantities are the KB integrals (KBIs). For both repulsive and attractive HCY fluids, the radial distribution functions are calculated by using the classical density functional theory, and then the corresponding KBIs are carried out. Since the local structure and global properties of a fluid can be related by KBI, we presented the isothermal compressibility and the derivative of the chemical potential with respect to bulk density for both repulsive and attractive HCY fluids. It is found that a transition of the affinity of particles in an attractive HCY fluid exists. The corresponding phase diagrams on the affinity are illustrated, which consist of repulsive and attractive regions with the boundary line of KBIs being zero. These results show that the aggregated structure of a HCY fluid can be effectively regulated by the screening parameter, bulk density and interaction energy, while KBIs can provide a quantitative reliable description on the properties of HCY fluids.

Molecular emulsions: from charge order to domain order

Aqueous mixtures of small molecules, such as lower n-alkanols for example, are known to be micro-segregated, with domains in the nano-meter range.

On the nature of the molecular ordering of water in aqueous DMSO mixtures

Computer simulation studies of aqueous dimethyl sulfoxyde (DMSO) mixtures show micro-heterogeneous structures, just like aqueous alcohol mixtures. However, there is a marked difference in the aggregate structure of water between the two types of systems. While water molecules form multiconnected globular clusters in alcohols, we report herein that the typical water aggregates in aqueous DMSO mixtures are linear, favouring a 2 hydrogen bond structure per water molecule, and for all DMSO mole fractions ranging from 0.1 to 0.9. This linear-aggregate structure produces a particular signature in the water site-site structure factors, in the form of a pre-peak at k ≈ 0.2-0.8 Å(-1), depending on DMSO concentration. This pre-peak is either absent in other aqueous mixtures, such as aqueous methanol mixtures, or very difficult to see through computer simulations, such as in aqueous-t-butanol mixtures. This difference in the topology of the aggregates explains why the Kirkwood-Buff integrals of aqueous-DMSO mixture look nearly ideal, in contrast with those of aqueous alcohol mixtures, suggesting a connection between the shape of the water aggregates, its fluctuations, and the concentration fluctuations. In order to further study this discrepancy between aqueous DMSO and aqueous alcohol mixture, two models of pseudo-DMSO are introduced, where the size of the sulfur atom is increased by a factor 1.6 and 1.7, respectively, hence increasing the hydrophobicity of the molecule. The study shows that these mixtures become closer to the emulsion type seen in aqueous alcohol mixtures, with more globular clustering of the water molecules, long range domain oscillations in the water-water correlations and increased water-water Kirkwood-Buff integrals. It demonstrates that the local ordering of the water molecules is influenced by the nature of the solute molecules, with very different consequences for structural properties and related thermodynamic quantities. This study illustrates the unique plasticity of water in presence of different types of solutes.

Microstructures of negative and positive azeotropes

Azeotropes famously impose fundamental restrictions on distillation processes, yet their special thermodynamic properties make them highly desirable for a diverse range of industrial and technological applications. Using neutron diffraction, this study provides first insights into the microstructures of azeotropes.

Simple and complex disorder in binary mixtures with benzene as a common solvent

Substituting benzene for water in computer simulations of binary mixtures, allows one to study the various forms of disorder, without the complications often encountered in aqueous mixtures.

Microstructure of neat alcohols

Formation of microstructure in homogeneous associated liquids is analyzed through the density-density pair correlation functions, both in direct and reciprocal space, as well as an effective local one-body density function. This is illustrated through a molecular dynamics study of two neat alcohols, namely, methanol and tert-butanol, which have a rich microstructure: chainlike molecular association for the former and micellelike for the latter. The relation to hydrogen bonding interaction is demonstrated. The apparent failure to find microstructure in water--a stronger hydrogen bonding liquid--with the same tools is discussed.

Ion-induced multiply reentrant liquid-liquid transitions and the nature of criticality in ethanol-water mixture

We report a quite unusual feature of four liquid-liquid reentrant transitions in ethanol (E)+water (W)+ammonium sulfate mixture by meticulous tuning of the ammonium sulfate concentration in a narrow range, as a function of temperature, at atmospheric pressure. Detailed exploration of the intricate phase behavior in terms of E/W sections shows that the range of triple reentrance shrinks with increasing E/W. The behavior of osmotic susceptibility is investigated by light scattering, near the critical point, in the one-phase region by varying the temperature at fixed concentration of the components, in a particular E/W section. The critical exponent of susceptibility (gamma) and correlation length (nu) are observed to have Fisher renormalized Ising values [Phys. Rev. 176, 237 (1968)], with gamma(r)=1.41 and nu(r)=0.718. The effective susceptibility exponent, gamma(eff), exhibits a sharp, nonmonotonic crossover from Ising to mean-field critical behavior, which is completed outside the critical regime. The amplitude of the correlation length, xi(o)(=21.2+/-0.4 A), deduced from light scattering experiment, is an order of magnitude larger than the typical values in usual aqueous electrolyte systems. This value of xi(o) is further verified from small-angle x-ray scattering (SAXS) experiments and found to be consistent. SAXS experiments on the critical sample reveal the presence of long-ranged intermolecular correlations, leading to supramolecular structuring, at a temperature far away from the critical point. These results convincingly demonstrate that the finite length scale arising due to the structuring competes with the diverging correlation length of critical concentration fluctuations, which influences the nonasymptotic critical behavior in this aqueous electrolyte system. The sulphate ions play a dominant role in both structuring and the complex phase behavior.