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.

The influence of charge ordering in the microscopic structure of monohydroxy alcohols

While radiation scattering data provides insight inside the microstructure of liquids, the Debye relation relating the scattering intensityI(k) to the atom-atom structure factorsSab(k)shows that, ultimately, it is these individual structure correlation functions which contain the relevant information about the micro-structure. However, these quantities are not observables, except in few cases where one can invert the Debye relation to obtain the structure functions. In the majority of other cases, the need for model dependent computer simulations is unavoidable. The resulting calculations reveal that the scattering pre-peak is the result of cancellations between positive pre-peaks and negative anti-peaks contributions from the atom-atom structure factors. What of systems where this cancellation is such that it entirely suppresses the scattering pre-peak? One would be tempted to falsely conclude that there is no underlying micro-heterogeneity. Hence, the structure functions appear as hidden variables, and it is important to understand the relation between their features and the micro-structure of the system. Through the computer simulation study of various mono-ols, ranging from methanol to 1-nonanol, as well as the branched octanols, we show how the features of the atom-atom pair correlation functiongab(r)affect that of the structure factorsSab(k), and reveal that the micro-structure is ultimately the result of the charge ordering between different atoms in the system.

High-Pressure and Temperature Effects on the Clustering Ability of Monohydroxy Alcohols

This study examined the clustering behavior of monohydroxy alcohols, where hydrogen-bonded clusters of up to a hundred molecules on the nanoscale can form. By performing X-ray diffraction experiments at different temperatures and under high pressure, we investigated how these conditions affect the ability of alcohols to form clusters. The pioneering high-pressure experiment performed on liquid alcohols contributes to the emerging knowledge in this field. Implementation of molecular dynamics simulations yielded excellent agreement with the experimental results, enabling the analysis of theoretical models. Here we show that at the same global density achieved either by alteration of pressure or temperature, the local aggregation of molecules at the nanoscale may significantly differ. Surprisingly, high pressure not only promotes the formation of hydrogen-bonded clusters but also induces the serious reorganization of molecules. This research represents a milestone in understanding association under extreme thermodynamic conditions in other hydrogen bonding systems such as water.

Experimental and Computational Approach to Studying Supramolecular Structures in Propanol and Its Halogen Derivatives

A series of four alcohols, n-propanol and its halogen (Cl, Br, and I) derivatives, were selected to study the effects of variation in polarity and halogen-driven interactions on the hydrogen bonding pattern and supramolecular structure by means of experimental and theoretical methods. It was demonstrated on both grounds that the average strength of H-bonds remains the same but dissociation enthalpy, the size of molecular nanoassemblies, as well as long-range correlations between dipoles vary with the molecular weight of halogen atom. Further molecular dynamics simulations indicated that it is connected to the variation in the molecular order introduced by specific halogen-based hydrogen bonds and halogen-halogen interactions. Our results also provided important experimental evidence supporting the assumption of the transient chain model on the molecular origin of the structural process in self-assembling alcohols.

Role of Charge Ordering in the Dynamics of Cluster Formation in Associated Liquids

Liquids are archetypes of disordered systems, yet liquids of polar molecules are locally more ordered than nonpolar molecules, due to the Coulomb interaction based charge ordering phenomenon. Hydrogen bonded liquids, such as water or alcohols, for example, represent a special type of polar liquids, in that they form labile clustered local structures. For water, in particular, hydrogen bonding and the related local tetrahedrality, play an important role in the various attempts to understand this liquid. However, labile structures imply dynamics, and it is not clear how it affects the understanding of this type of liquids from purely static point of view. Herein, we propose to reconsider hydrogen bonding as a charge ordering process. This concept allows us to demonstrate the insufficiency of the analysis of the microscopic structure based solely on static pair correlation functions, and the need for dynamical correlation functions, both in real and reciprocal space. The subsequent analysis allows to recover several aspects of our understanding of hydrogen bonded liquids, but from the charge order perspective. For water, it confirms the jump rotation picture found recently, and it allows to rationalize the contradicting pictures that arise when following the interpretations based on hydrogen bonding. For alcohols, it allows to understand the dynamical origin of the scattering prepeak, which does not exist for water, despite the fact that both these liquids have very similar hydroxyl group chain clusters. The concept of charge ordering complemented by the analysis of dynamical correlation functions appear as a promising way to understand microheterogeneity in complex liquids and mixtures from kinetics point of view.

Non-simple flow behavior in a polar van der Waals liquid: Structural relaxation under scrutiny

The non-exponential character of the structural relaxation is considered one of the hallmarks of the glassy dynamics, and in this context, the relatively narrow shape observed by dielectric techniques for polar glass formers has attracted the attention of the community for long time. This work addresses the phenomenology and role of specific non-covalent interactions in the structural relaxation of glass-forming liquids by the study of polar tributyl phosphate. We show that dipole interactions can couple to shear stress and modify the flow behavior, preventing the occurrence of the simple liquid behavior. We discuss our findings in the general framework of glassy dynamics and the role of intermolecular interactions.

Structure of water-in-salt and water-in-bisalt electrolytes

We report a systematic diffraction study of two "water-in-salt" electrolytes and a "water-in-bisalt" electrolyte combining high-energy X-ray diffraction (HEXRD) with polarized and unpolarized neutron diffraction (ND) on both H₂O and D₂O solutions. The measurements provide three independent combinations of correlations between the different pairs of atom types that reveal the short- and intermediate-range order in considerable detail. The ND interference functions show pronounced peaks around a scattering vector Q ∼ 0.5 Å^-1 that change dramatically with composition, indicating significant rearrangements of the water network on a length scale around 12 Å. The experimental results are compared with two sets of Molecular Dynamics (MD) simulations, one including polarization effects and the other based on a non-polarizable force field. The two simulations reproduce the general shapes of the experimental structure factors and their changes with concentration, but differ in many detailed respects, suggesting ways in which their force fields might be modified to better represent the actual systems.

Spontaneous Ouzo Emulsions Coexist with Pre-Ouzo Ultraflexible Microemulsions

Even in the absence of surfactants, polymers, or particles, spontaneous emulsions produced by dilution with water can be stable over days. This "Ouzo effect" used by the industry is obtained by rapid dilution from an identified "pre-Ouzo" domain of composition where weak aggregates are present: nanometer-sized clusters covered by a surface layer enriched in a hydrotrope such as ethanol. In these systems, Ostwald ripening is not an effective destabilizing mechanism. Using in situ autodilution small-angle X-ray scattering (SAXS), we follow the morphological transitions occurring in a ternary mixture of water/n-octanol/ethanol throughout the monophasic and biphasic regions. This allows for the first time an online characterization of the multiscale coexisting microstructures. Small-angle neutron scattering (SANS) profiles on metastable emulsions as well as phase-separated samples complete the SAXS data, taking advantage of contrast variation via isotopic substitution. After crossing the phase boundary into the two-phase region, coexisting phases are both ternary solutions structured at the nanometer scale when the emulsion is stable. The transition from single phase to two phases is asymmetric around the plait point. When the initial concentration of the hydrotrope is below the minimum hydrotrope concentration (MHC), emulsification failure occurs, i.e., emulsions cream within seconds. Beyond MHC, the low interfacial tension between coexisting ternary fluids results in a Laplace pressure below 100 Pa, explaining the puzzling resilience of spontaneous emulsion to the universal mechanism of Ostwald ripening.

On the X-ray Scattering Pre-peak of Linear Mono-ols and the Related Microstructure from Computer Simulations

The X-ray scattering intensities (I(k)) of linear alkanols OH(CH₂)_n-1CH₃ obtained from experiments (methanol to 1-undecanol) and computer simulations (methanol to 1-nonanol) of different force field models are comparatively studied particularly in order to explain the origin and the properties of the scattering pre-peak in the k-vector range 0.3-1 Å^-1. The experimental I(k) values show two apparent features: the pre-peak position k_P decreases with increasing n, and more intriguingly, the amplitude A_P goes through a maximum at 1-butanol (n = 4). The first feature is well reproduced by all force-field models, while the second shows strong model dependence. The simulations reveal various shapes of clusters of the hydroxyl head-group from n>2. k_P is directly related to the size of the meta-objects corresponding to such clusters surrounded by their alkyl tails. The explanation of the A_P turnover at n = 4 is more involved in terms of cancellations of atom-atom structure factor S(k) contributions related to domain ordering. The flexibility of the alkyl tails tends to reduce the cross contributions, thus revealing the crucial importance of this parameter in the models. Force fields with all-atom representation are less successful in reproducing the pre-peak features for smaller alkanols, n<6, possibly because they blur the charge ordering process since all atoms bear partial charges. The analysis clearly shows that it is not possible to obtain a model-free explanation of the features of I(k).

Coupling between Structural and Dielectric Relaxations of Methanol and Ethanol Studied by Molecular Dynamics Simulation

The microscopic origin of the fast dielectric relaxation modes and the integrated dielectric relaxation times of methanol and ethanol was investigated by means of cross-correlation analysis of molecular dynamics simulation. Random force on the fluctuation of collective dipole moment was correlated with the two-body density mode in both real and reciprocal spaces. A strong coupling was observed with the OH alternation mode at 30 nm^-1, suggesting that alternating switching of the hydrogen bond within a hydrogen-bonding chain is the principal origin of the retarded friction on the collective dipole moment. The relaxation of the coupling was much slower than that of the partial intermediate scattering functions at the corresponding wavenumber, which suggests the breakdown of the factorization approximation employed in the mode-coupling theory. Although the prepeak structure is strongly coupled to the viscoelastic relaxation, its coupling with the dielectric relaxation is relatively weak. The difference between the viscoelastic and the dielectric relaxations was discussed in terms of the different symmetries of the shear stress tensor and the collective dipole moment.

Modeling micro-heterogeneity in mixtures: The role of many body correlations

A two-component interaction model is introduced herein, which allows us to describe macroscopic miscibility with various modes of tunable micro-segregation, ranging from phase separation to micro-segregation, and is in excellent agreement with structural quantities obtained from simulations and the liquid state hypernetted-chain like integral equation theory. The model is based on the conjecture that the many-body correlation bridge function term in the closure relation can be divided into one part representing the segregation effects, which are modeled herein, and the usual part representing random many body fluctuations. Furthermore, the model allows us to fully neglect these second contributions, thus increasing the agreement between the simulations and the theory. The analysis of the retained part of the many body correlations gives important clues about how to model the many body bridge functions for more realistic systems exhibiting micro-segregation, such as aqueous mixtures.

Microscopic origin of the scattering pre-peak in aqueous propylamine mixtures: X-ray and neutron experiments versus simulations

The structure of aqueous propylamine mixtures is investigated through X-ray and neutron scattering experiments, and the scattered intensities compared with computer simulation data.

Analysis of shear viscosity and viscoelastic relaxation of liquid methanol based on molecular dynamics simulation and mode-coupling theory

The role of the prepeak structure of liquid methanol in determining its shear viscosity was studied by means of molecular dynamics (MD) simulation and mode-coupling theory (MCT). The autocorrelation function of the shear stress and the intermediate scattering functions at both the prepeak and the main peak were calculated from the MD trajectories. Their comparison based on MCT suggests that the viscoelastic relaxation in the ps regime is affected by the slow structural dynamics at the prepeak. On the other hand, the MCT for molecular liquids based on the interaction-site model (site-site MCT) fails to describe the coupling between the prepeak dynamics and shear stress. The direct evaluation of the coupling between the two-body density and the shear stress reveals that the viscoelastic relaxation is actually affected by the prepeak dynamics, although the coupling is not captured by the site-site MCT. The site-site MCT works well for a model methanol without partial charges, suggesting that the failure of the site-site MCT originates from the existence of a hydrogen-bonding network structure.

Multimodal character of shear viscosity response in hydrogen bonded liquids

Non-simple viscosity response of 2E1H alcohol forming supramolecular aggregates.

Dynamic heterogeneity in aqueous ionic solutions

It is well known that supercooled liquids have heterogeneous dynamics, but it is still unclear whether dynamic heterogeneity also exists in aqueous ionic solutions at room or even higher temperatures.

Analysis of Prepeak Structure of Concentrated Organic Lithium Electrolyte by Means of Neutron Diffraction with Isotopic Substitution and Molecular Dynamics Simulation

The prepeak structure of a 3 mol/kg solution of LiClO₄ in propylene carbonate (PC) was studied by both neutron diffraction with isotopic substitution (NDIS) and molecular dynamics (MD) simulation. The NDIS data showed that the intensity of the prepeak decreases experimentally with an increase in the scattering length of the lithium atom from ⁷Li to ⁶Li in PC-d₆. On the other hand, although the prepeak was observed in solutions of both PC-d₆ and PC-h₆, it disappears when the 1:1 mixture of PC-d₆ and PC-h₆ was used as the solvent. The prepeak structure and its variation with the isotope substitution were reproduced well by MD simulation, and they were explained in terms of the contrast of the scattering length densities of the ionic and nonpolar domains.

Lifshitz phase: the microscopic structure of aqueous and ethanol mixtures of 1,n-diols

We study binary mixtures of ethylene glycol and 1,3-propandiol with water or ethanol using computer simulations.

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.