Structure and dynamics in aqueous mixtures of glycerol: insights from molecular dynamics simulations

Aqueous glycerol mixtures are investigated over the whole concentration range of glycerol xGLY = 0.1-0.9 via molecular dynamics (MD) simulations at ambient pressure and temperature. Two glycerol force fields are used: an all-atom (AA) and a united-atom (UA) model. Structural changes upon different mixing ratios are discussed through the site-site radial distribution functions (RDFs), coordination numbers and cluster analysis. As both species are hydrogen bonded, they form an almost perfect H-bonded network, with no observed clusters. There are, however, noticeable changes in the RDFs. Glycerol correlations grow stronger with increasing glycerol content, as do water correlations. There is significant transformation in dynamics as well, as evidenced by the self-diffusion coefficients, the velocity autocorrelation functions and the rotational autocorrelation functions. Diffusion of both species slows down with increasing glycerol content. Rotational relaxation is also altered depending on the mixture composition and there is a slow-down at the lower end of glycerol content.

Nanoclustering in non-ideal ethanol/heptane solutions alters solvation dynamics

Alcohol/alkane solutions widely used in chemical synthesis and as transportation fuels are highly non-ideal due to the nanoscale clustering of the amphiphilic alcohol molecules within the nonpolar alkanes. Besides impacting reactivity, such as combustion, non-ideal solutions are likely to exhibit unusual solvation dynamics on ultrafast time scales arising from the structurally heterogeneous nature of molecular-scale association. Using a convenient transition metal carbonyl vibrational probe [(C5H5)Mn(CO)3, CMT], linear absorption and nonlinear two-dimensional infrared (2D-IR) spectroscopy reveal composition-dependent solvation dynamics as reported by the frequency fluctuation correlation function in a series of ethanol/heptane solutions. Slow spectral diffusion with dilute ethanol indicates preferential solvation of the polar solute by the alcohol with a mechanism largely dominated by solvent exchange. Comparison with an ethanol/acetonitrile solution series yields no substantial preferential solvation or solvent exchange signatures in the linear or 2D-IR spectra. In ethanol/heptane solutions, increasing the ethanol concentration speeds up the solvation dynamics, which is largely consistent with a model that includes solvent exchange and single-solvent spectral diffusion. Detailed analysis of the deviation from the experimental time constants from the model's optimal parameters yields a remarkable resemblance of the concentration-weighted Kirkwood-Buff integrals for ethanol/heptane solutions. This trend indicates that solution non-ideality alters the spectral diffusion dynamics of the probe solute. Given that nanoscale clustering drives the non-ideality, these experiments reveal a dynamical consequence of nanoscale heterogeneity on the ultrafast dynamics of the solution. Refined understanding of the structural and dynamical aspects of mixed solvents will be necessary for predictive solution strategies in chemistry.

Unraveling the interplay of temperature with molecular aggregation and miscibility in TEA-water mixtures

In the phase diagram of binary liquid mixtures, a miscibility gap is found with the concomitant liquid-liquid phase separation, wherein temperature is a key parameter in modulating the phase behavior. This includes critical temperatures such as the lower critical solution temperature (LCST) and upper critical solution temperature (UCST). Using a comprehensive approach including molecular dynamics (MD) simulation, graph theoretical analysis and spatial inhomogeneity measurement in an LCST-type mixture, we attempt to establish the relationship between the molecular aggregation pattern and phase behavior in TEA-water mixtures. At lower temperatures of binary liquid mixtures, TEA molecules tend to aggregate while simultaneously interacting with water forming a homogeneous solution. As the temperature increases, these TEA aggregates tend to self-associate by minimizing the interaction with water, which facilitates formation of two distinct liquid phases in the binary liquid. The spatial distribution analysis also reveals that the TEA aggregates compatible with water promote uniform distribution of water molecules, maintaining a homogeneous solution, while the water-incompatible ones generate isolation of water H-bond aggregates, leading to liquid-liquid phase separation in the binary system. This current study on temperature-induced molecular aggregation behavior is anticipated to contribute to a critical understanding of the phase behavior in binary liquid mixtures, including UCST, LCST, and reentrant phase behavior.

Elucidating the Quenching Mechanism of Tris(2,2'-bipyridyl)ruthenium(II) Complex in the Water-Glycerol Binary System Based on the Microscopic Structure of the Media

Kinetic studies on the photochemical quenching reaction of the tris(2,2'-bipyridyl) ruthenium(II) complex ([Ru(bpy)3]2+) in water-glycerol binary media were conducted based on the Einstein-Smoluchowski (E-S) theory. Dynamic and static quenching behaviors were analyzed by comparing results from time-resolved spectroscopy and emission spectroscopy. While the dynamic quenching reaction aligns well with the E-S theory, static quenching was observed, leading to a notable increase in the overall photoquenching reaction rate constant. Employing chromatography and infrared spectroscopy, we correlated the microscopic molecular structure of the binary solvent system and the solvation environment around the emitters with the reaction mechanism. This correlation was found to correspond to ion pair formation and the confinement effect of the emitter, respectively.

On the dual behaviour of water in octanol-rich aqueous n-octanol mixtures: an X-ray scattering and computer simulation study

Aqueous n-octanol (n = 1, 2, 3, and 4) mixtures from the octanol rich side are studied by X-ray scattering and computer simulation, with a focus on structural changes, particularly in what concerns the hydration of the hydroxyl-group aggregated chain-like structures, under the influence of various branching of the alkyl tails. Previous studies have indicated that hydroxyl-group chain-cluster formation is hindered in proportion to the branching number. Here, water mole fractions up to x = 0.2 are examined, i.e. up to the miscibility limit. It is found that water molecules within the hydroxyl-chain domains participate in the chain formations in a different manner for 1-octanol and the branched octanols. The hydration of the octanol hydroxyl chains is confirmed by the shifting of the scattering pre-peak position kPP to smaller values, both from measured and simulated X-ray scattering intensities, which corresponds to an increased size of the clusters. Experimental pre-peak amplitudes are seen to increase with increasing water content for 1-octanol, while this trend is reversed in all branched octanols, with the amplitudes decreasing with the increase of the branching number. Conjecturing that the amplitudes of pre-peaks are related to the density of the corresponding aggregates, these results are interpreted as water breaking large OH hydroxyl chains in 1-octanol, hence increasing the density of aggregates, while enhancing hydroxyl aggregates in branched alcohols by inserting itself into the OH chains. The analysis of the cluster distributions from computer simulations provide more details on the role of water. For cluster sizes smaller than dc = 2π/kPP, water is found to always play the role of a structure enforcer for all n-octanols, while for clusters of size dc water is always a destructor. For cluster sizes larger than dc, the role of water differs from 1-octanol and the branched ones: it acts as a structure maker or breaker in inverse proportion to the hindering of OH hydroxyl chain structures arising from the topology of the alkyl tails (branched or not).

Molecular Aggregation Behavior and Microscopic Heterogeneity in Binary Osmolyte-Water Solutions

Osmolytes, small organic compounds, play a key role in modulating the protein stability in aqueous solutions, but the operating mechanism of the osmolyte remains inconclusive. Here, we attempt to clarify the mode of osmolyte action by quantitatively estimating the microheterogeneity of osmolyte-water mixtures with the aid of molecular dynamics simulation, graph theoretical analysis, and spatial distribution measurement in the four osmolyte solutions of trimethylamine-N-oxide (TMAO), tetramethylurea (TMU), dimethyl sulfoxide, and urea. TMAO, acting as a protecting osmolyte, tends to remain isolated with no formation of osmolyte aggregates while preferentially interacting with water, but there is a strong aggregation propensity in the denaturant TMU solution, characterized by favored hydrophobic interactions between TMU molecules. Taken together, the mechanism of osmolyte action on protein stability is proposed as a comprehensive one that encompasses the direct interactions between osmolytes and proteins and indirect interactions through the regulation of water properties in the osmolyte-water mixtures.

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.

Polymer Dynamics in Glycerol-Water Mixtures

Velocity correlation spectra (VAS) in binary mixtures of water and glycerol (G/W), obtained by measurements using the modulated gradient spin echo (MGSE) NMR method, were explained by the interactions of water molecules with clusters formed around the hydrophilic glycerol molecule, which drastically change the molecular dynamics and rheology of the mixture. It indicates a thickening of the shear viscosity, which could affect the dynamics of submerged macromolecules. The calculation of the polymer dynamics with the Langevin equations according to the Rouse model, where the friction was replaced by the memory function of the retarded friction, gave the dependence of the dynamics of the polymer on the rate of shear viscous properties of the solvent. The obtained formula was used to calculate the segmental VAS of the polymer when immersed in pure water and in a G/W mixture with 33 vol% glycerol content, taking into account the inverse proportionality between the solvent VAS and friction. The spectrum shows that in the G/W mixture, the fast movements of the polymer segments are strongly inhibited, which creates the conditions for slow processes caused by the internal interaction between the polymer segments, such as interactions that cause disordered polypeptides to spontaneously fold into biologically active protein molecules when immersed in such a solvent.

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.

Solvent Effects in Ligand Stripping Behavior of Colloidal Nanoparticles

Inorganic colloidal nanoparticle (NP) properties can be tuned by stripping stabilizing ligands using a poor solvent. However, the mechanism behind ligand stripping is poorly understood, in part because in situ measurements of ligand stripping are challenging at the nanoscale. Here, we investigate ethanol solvent-mediated oleylamine ligand stripping from magnetite (Fe₃O₄) NPs in different compositions of ethanol/hexane mixtures using atomistic molecular dynamics (MD) simulations and thermogravimetric analysis (TGA). Our study elucidates a complex interplay of ethanol interactions with system components and indicates the existence of a threshold concentration of ∼34 vol % ethanol, above which ligand stripping saturates. Moreover, hydrogen bonding between ethanol and stripped ligands inhibits subsequent readsorption of the ligands on the NP surface. A proposed modification of the Langmuir isotherm explains the role of the enthalpy of mixing of the ligands and solvents on the ligand stripping mechanism. A good agreement between the MD predictions and TGA measurements of ligand stripping from Fe₃O₄ NPs validates the simulation observations. Our findings demonstrate that the ligand coverage of NPs can be controlled by using a poor solvent below the threshold concentration and highlight the importance of ligand-solvent interactions that modulate the properties of colloidal NPs. The study also provides an approach for a detailed in silico study of ligand stripping and exchange from colloidal NPs that are crucial for applications of NPs spanning self-assembly, optoelectronics, nanomedicine, and catalysis.

Detecting Liquid-Liquid Phase Separations Using Molecular Dynamics Simulations and Spectral Clustering

A stringent test of the accuracy of empirical force fields is reproducing the phase diagram of bulk phases and mixtures. Exploring the phase diagram of mixtures requires the detection of phase boundaries and critical points. In contrast to most solid-liquid transitions, in which a global order parameter (average density) can be used to discriminate between two phases, some demixing transitions entail relatively subtle changes in the local environment of each molecule. In such cases, finite sampling errors and finite-size effects make the identification of trends in local order parameters extremely challenging. Here we analyze one such example, namely a methanol/hexane mixture, and compute several local and global structural properties. We simulate the system at various temperatures and study the structural changes associated with demixing. We show that despite a seemingly continuous transformation between mixed and demixed states, the topological properties of the H-bond network change abruptly as the system crosses the demixing line. In particular, by using spectral clustering, we show that the distribution of cluster sizes develops a fat tail (as expected from percolation theory) in the vicinity of the critical point. We illustrate a simple criterion to identify this behavior, which results from the emergence of large system-spanning clusters from a collection of aggregates. We further tested the spectral clustering analysis on a Lennard-Jones system as a standard example of a system with no H-bonds, and also, in this case, we were able to detect the demixing transition.

The impact of H/D exchange on the thermal and structural properties as well as high-pressure relaxation dynamics of melatonin

In this paper, thermal properties, atomic-scale structure, and molecular dynamics (at ambient and high pressure) of native melatonin (MLT) and its partially-deuterated derivative (MLT-d₂) have been investigated. Based on infrared spectroscopy, it was shown that treating MLT with D₂O causes the replacement of hydrogen atoms attached to the nitrogen by deuterium. The degree of such substitution was very high (> 99%) and the deuterated sample remained stable after exposure to the air as well as during the melting and vitrification processes. Further calorimetric studies revealed the appearance of a peculiar thermal event before the melting of crystalline MLT-d₂, which was assigned by the X-ray diffraction to a local negative thermal expansion of the unit cell. Finally, the high-pressure dielectric experiments indicated a few interesting findings, including the variation in the shape of the structural relaxation peak during compression, the difference in the pressure evolution of the glass transition temperature, and the temperature dependence of activation volume for both MLT species. The variations in these parameters manifest a different impact of the compression/densification on the dynamics of hydrogen and deuterium bonds in the native and partially-deuterated MLT, respectively.

Supramolecular structures of self-assembled oligomers under confinement

We study the molecular origin of a prepeak (PP) observed at low q values in the structure factors of three oligomers in a bulk (poly(mercaptopropyl)methylsiloxane, PMMS, poly(methylmercaptopropyl)-grafted-hexylmethacrylate, PMMS-g-HMA, and poly(methylphenyl)siloxane, PMPS) in order to understand the lowering of the PP intensity detected for oligomers confined in cylindrical pores with low diameter. For this purpose, we use a combination of X-ray diffraction measurements and coarse-grained bead-spring molecular dynamics simulations. Our molecular modelling demonstrated that the planarity of the pendant groups triggers the self-association of oligomers into nanoaggregates. However, the formation of oligomeric nanodomains is not sufficient for building-up the PP. The latter requires spatial disturbance in the arrangement of the side groups of oligomers within clusters. Importantly, our numerical analysis revealed that the increasing degree of the confinement of oligomers limits their aggregation and consequently lowers the amplitude of the PP observed in the experimental data.

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.

Supramolecular Structure of Phenyl Derivatives of Butanol Isomers

Wide-angle X-ray scattering patterns were recorded for a series of aliphatic butanol isomers (n-, iso-, sec-, tert-butanol) and their phenyl derivatives (4-phenyl-1-butanol, 2-methyl-3-phenyl-1-propanol, 4-phenyl-2-butanol, and 2-methyl-1-phenyl-2-propanol, respectively) to determine their atomic-scale structure with particular emphasis on the formation of supramolecular clusters. In addition, molecular dynamics simulations were carried out and yielded good agreement with experimental data. The combination of experimental and theoretical results allowed clarification of the origin of the pre-peak appearing at low scattering angles for the aliphatic butanols and its absence for their phenyl counterparts. It was demonstrated that the location of the hydroxyl group in the molecule of alkyl butanol, its geometry, and rigidity determine the morphology of the supramolecular clusters, while the addition of the aromatic moiety causes more disordered organization of molecules. The phenyl group significantly decreases the number of hydrogen bonds and size of the supramolecular clusters formed via the O–H···O scheme. The lower association ability of phenyl alcohols via H-bonds is additionally attenuated by the appearance of competing π–π configurations evidenced by the structural models.

Dynamical Anomaly of Aqueous Amphiphilic Solutions: Connection to Solution H-Bond Fluctuation Dynamics?

We have investigated the possible connection between "dynamical anomaly" observed in time-resolved fluorescence measurements of reactive and nonreactive solute-centered relaxation dynamics in aqueous binary mixtures of different amphiphiles and the solution intra- and interspecies H-bond fluctuation dynamics. Earlier studies have connected the anomalous thermodynamic properties of binary mixtures at very low amphiphile concentrations to the structural distortion of water. This is termed as "structural anomaly." Interestingly, the abrupt changes in the composition-dependent average rates of solute relaxation dynamics occur at amphiphile mole fractions approximately twice as large as those where structural anomalies appear. We have investigated this anomalous solution dynamical aspect by considering (water + tertiary butanol) as a model system and performed molecular dynamics simulations at several tertiary butanol (TBA) concentrations covering the extremely dilute to the moderately concentrated regimes. The "dynamical anomaly" has been followed via monitoring the composition dependence of the intra- and interspecies H-bond fluctuations and reorientational relaxations of TBA and water molecules. Solution structural aspects have been followed via examining the tetrahedral order parameter, radial and spatial distribution functions, numbers of H bonds per water and TBA molecules, and the respective populations participating in H-bond formation. Our simulations reveal abrupt changes in the H-bond fluctuations and reorientational dynamics and tetrahedral order parameter at amphiphile concentrations differing approximately by a factor of 2 and corroborates well with the steady-state and the time-resolved spectroscopic measurements. This work therefore explains, following a uniform and cogent manner, both the experimentally observed structural and dynamical anomalies in microscopic terms.

All-Atom Simulation Method for Zeeman Alignment and Dipolar Assembly of Magnetic Nanoparticles

Magnetic nanoparticles (MNPs) can organize into novel structures in solutions with excellent order and unique geometries. However, studies of the self-assembly of smaller MNPs are challenging due to a complicated interplay between external magnetic fields and van der Waals, electrostatic, dipolar, steric, and hydrodynamic interactions. Here, we present a novel all-atom molecular dynamics simulation method to enable detailed studies of the dynamics, self-assembly, structure, and properties of MNPs as a function of core sizes and shapes, ligand chemistry, solvent properties, and external field. We demonstrate the use and effectiveness of the model by simulating the self-assembly of oleic acid ligand-functionalized magnetite (Fe₃O₄) nanoparticles, with spherical and cubic shapes, into rings, lines, chains, and clusters under a uniform external magnetic field. We found that the long-range electrostatic interactions can favor the formation of a chain over a ring, the ligands promote MNP cluster growth, and the solvent can reduce the rotational diffusion of the MNPs. The algorithm has been parallelized to take advantage of multiple processors of a modern computer and can be used as a plugin for the popular simulation software LAMMPS to study the behavior of small MNPs and gain insights into the physics and chemistry of different magnetic assembly processes with atomistic details.

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.

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.

Universal features in the lifetime distribution of clusters in hydrogen-bonding liquids

Hydrogen-bonding liquids, typically water and alcohols, are known to form labile structures (network, chains, etc.); hence, the lifetime of these structures is an important microscopic parameter, which can be calculated via computer simulations. Since these cluster entities are mostly statistical in nature, one would expect that, in the short-timescale regime, their lifetime distribution would be a broad Gaussian-like function of time, with a single maximum representing their mean lifetime, and be weakly dependent on criteria such as the bonding distance and angle, much similar to non-hydrogen-bonding simple liquids, while the long-timescale regime is known to have some power law dependence. Unexpectedly, all the hydrogen-bonding liquids studied herein, namely water and alcohols, display three highly hierarchical specific lifetimes, in the sub-picosecond range 0-0.5 ps. The dominant lifetime depends very strongly on the bonding-distance criterion and is related to hydrogen-bonded pairs. This mode is absent in non-H-bonding simple liquids. The secondary and tertiary mean lifetimes are related to clusters and are nearly independent of the bonding criterion. Of these two lifetimes, only the first one can be related to that of simple liquids, which poses the question of the nature of the third lifetime. The study of alcohols reveals that this third lifetime is related to the topology of the H-bonded clusters and that its distribution may also be affected by the alkyl tail surrounding the "bath". This study shows that hydrogen-bonding liquids have a universal hierarchy of hydrogen-bonding lifetimes with a timescale regularity across very different types, and which depend on the topology of the cluster structures.

Mesostructuring in Liquid-Liquid Extraction Organic Phases Originating from Critical Points

Organic phase structure plays an important role in solute extraction energetics and phase behavior of liquid-liquid extraction (LLE) systems. For a binary extractant (amphiphile)/solvent mixture of relevance to LLE, we find that the organic phase mesostructuring is consistent with extractant concentration fluctuations as the compositional isotherm traverses the Widom line above its liquid-liquid critical point. This reveals a different mechanism for the well-documented heterogeneities in LLE organic phases that are typically attributed to micellization.

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.

Length-Scale Effects in Hydrophobic Polymer Collapse Transitions

The study of molecular mechanisms for cosolvent-driven hydrophobic polymer collapse transitions in water is of pivotal importance in the field of smart responsive materials. Computational studies together with complementary experimental data have led to the discovery and understanding of new phenomena in recent years. However, elementary mechanisms, generally contributing to polymer coil-globule transitions in different classes of cosolvent-water systems, remain elusive due to compensating energy-entropy effects. Herein, I discuss the role of length scales in polymer solubility problems. New ideas on surfactant mechanisms are discussed based on examples in which these mechanisms drive polymer swelling or collapse.

Peptoids as a Chiral Stationary Phase for Liquid Chromatography: Insights from Molecular Dynamics Simulations

Peptoids are peptide regioisomers with attractive structural tunability in terms of sequence and three-dimensional arrangement. Peptoids are foreseen to have a great potential for many diverse applications, including their utilization as a chiral stationary phase in chromatography. To achieve chiral recognition, a chiral side chain is required to allow specific interactions with a given enantiomer from a racemic mixture. One of the most studied chiral stationary phases, built with (S)-N-1-phenylethyl (Nspe) units, was shown to be successful in resolving racemic mixtures of binaphthyl derivatives. However, there is currently no description at the atomic scale of the factors favoring its enantioselectivity. Here, we take advantage of steered molecular dynamics simulations to mimic the elution process at the atomic scale and present evidence that the predominantly right-handed helical conformation of Nspe peptoids and their ability to form stronger hydrogen bonds with the (S) enantiomer are responsible for the chiral recognition of the popular chiral probe 2,2'-bihydroxy-1,1'-binaphthyl.

Cluster Identification Using Modularity Optimization to Uncover Chemical Heterogeneity in Complex Solutions

Structural heterogeneity is commonly manifested in solutions and liquids that feature competition of different interparticle forces. Identifying and characterizing heterogeneity across different length scales requires multimodal experimental measurement and/or the application of new techniques for the interrogation of atomistic simulation data. Within the latter, the parsing of networks of interparticle interactions (chemical networks) has been demonstrated to be a valuable tool for identifying subensembles of chemical environments. However, chemical networks can adopt a wide variety of topologies that challenge generalizable methods for identifying heterogeneous behavior, and few network analysis algorithms have been proposed for multiscale resolution. In this study, we apply a method of partitioning using the graph theoretic concept of clusters and communities. Using a modularity optimization algorithm, the cluster partition creates subgraphs based on their relative internal and external connectivities. The methodology is tested on two soft matter systems that have significantly different network topologies so as to probe its ability to identify multiple scale features and its generalizability. A binary Lennard-Jones fluid is first examined, where one component causes subgraphs that have high internal network connectivity yet are still connected to the rest of the interparticle network of interactions. The impact of connectivity and edge weighting on the cluster partition is investigated. In the second system, hierarchically organized molecular structures comprised of hydrogen bonded water molecules are identified at a liquid/liquid interface. These structures have a much more sparse network with significantly varied internal connectivity that is a challenge to differentiate from the background hydrogen bonding network of water molecules at the instantaneous interface. The organized macrostructures are effectively isolated from the background network using the cluster partition, and a time-dependent implementation allows us to reveal their reactivity. These studies indicate that cluster partitioning based upon intermolecular network connectivity patterns is broadly generalizable, depending only on user-defined intermolecular connectivity, is operable across different length scales, and is extensible to the study of dynamic phenomena.

Diffusion of the carbon dioxide-ethanol mixture in the extended critical region

The effect of traces of ethanol in supercritical carbon dioxide on the mixture's thermodynamic properties is studied by molecular simulations and Taylor dispersion measurements. This mixture is investigated along the isobar p = 10 MPa in the temperature range between T = 304 and 343 K. Along this path, the mixture undergoes two transitions: First, the Widom line is crossed, marking the transition from liquid-like to gas-like conditions. A second transition occurs from the supercritical gas-like domain to a subcritical gas. The Widom line crossover entails inflection points for most of the studied properties, i.e. density, enthalpy, shear viscosity, Maxwell-Stefan and intradiffusion coefficients. On the other hand, the transition between the super- and subcritical regions is found to be generally smooth, an observation that is qualitatively confirmed by experimental Taylor dispersion measurements. Dedicated atomistic simulations show the presence of microheterogeneities due to ethanol self-association along the investigated path, which lead to the mixture's anomalous behavior in its extended critical region.

Theoretical and Experimental Study of the Excess Thermodynamic Properties of Highly Nonideal Liquid Mixtures of Butanol Isomers + DBE

Binary alcohol + ether liquid mixtures are of significant importance as potential biofuels or additives for internal combustion engines and attract considerable fundamental interest as model systems containing one strongly H-bonded self-associating component (alcohol) and one that is unable to do so (ether), but that can interact strongly as a H-bond acceptor. In this context, the excess thermodynamic properties of these mixtures, specifically the excess molar enthalpies and volumes (H^E and V^E), have been extensively measured. Butanol isomer + di-n-butyl ether (DBE) mixtures received significant attention because of interesting differences in their V^E, changing from negative (1- and isobutanol) to positive (2- and tert-butanol) with increasing alkyl group branching. With the aim of shedding light on the differences in alcohol self-association and cross-species H-bonding, considered responsible for the observed differences, we studied representative 1- and 2-butanol + DBE mixtures by molecular dynamics simulations and experimental excess property measurements. The simulations reveal marked differences in the self-association of the two isomers and, while supporting the existing interpretations of the H^E and V^E in a general sense, our results suggest, for the first time, that subtle changes in H-bonded topologies may contribute significantly to the anomalous volumetric properties of these mixtures.

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.

Self-Assembled Solute Networks in Crowded Electrolyte Solutions and Nanoconfinement of Charged Redoxmer Molecules

Redoxmers are electrochemically active organic molecules storing charge and energy in electrolyte fluids circulating through redox flow batteries (RFBs). Such molecules typically have solvent-repelling cores and solvent-attracting pendant groups introduced to increase solubility in liquid electrolytes. These two features can facilitate nanoscale aggregation of the redoxmer molecules in crowded solutions. In some cases, this aggregation leads to the emergence of continuous networks of solute molecules in contact, and the solution becomes microscopically heterogeneous. Here, we use small-angle X-ray scattering (SAXS) and molecular dynamics modeling to demonstrate formation of such networks and examine structural factors controlling this self-assembly. We also show that salt ions become excluded from these solute aggregates into small pockets of electrolytes, where these ions strongly associate. This confinement by exclusion is also likely to occur to charged redoxmer molecules in a "sea" of neutral precursors coexisting in the same solution. Here, we demonstrate that the decay lifetime of the confined charged molecules in such solutions can increase several fold compared to dilute solutions. We attribute this behavior to a "microreactor effect" on reverse reactions of the confined species during their decomposition.

Comparative analysis of ethanol dynamics in aqueous and non-aqueous solutions

In this study, we compare the results for vibrational, reorientational and hydrogen bond dynamics of ethanol in water and in hexane across the whole concentration range. Water and hexane are both commonly used as solvents, but so far, it has been unclear to what extent they modify the solute dynamics. Ethanol is chosen as the solute because it is an aliphatic molecule that is miscibile with both solvents. It is known that ethanol forms micelle-like domains in water and cyclic clusters resembling loops in hexane. This structural micro-heterogeneity is well known both in experiments and in simulations. The main question we raise here is: is there a signature of micro-heterogeneity in the dynamical quantities of ethanol? We focus on quantities such as the vibrational spectra, the reorientational correlation functions, the self-diffusion coefficients, the ethanol-ethanol hydrogen bond correlation functions and the corresponding hydrogen bond histograms. For the first time ever, we compute the van Hove functions to reveal the dynamical variations of spatial correlations in these systems. All these results complement each other and provide a unifying dynamical description of ethanol in binary mixtures.

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).

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.