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.

First-Principles Simulations of CuCl in High-Temperature Water Vapor.

Experimental data suggest that the solubility of copper in high-temperature water vapor is controlled by the formation of hydrated clusters of the form CuCl(H2O)n, where the average number of water molecules in the cluster generally increases with increasing density [Migdisov, A. A.; et al. Geochim. Cosmochim. Acta 2014, 129, 33-53]. However, the precise nature of these clusters is difficult to probe experimentally. Moreover, there are some discrepancies between experimental estimates of average cluster size and prior simulation work [Mei, Y. Geofluids 2018, 2018, 4279124]. We have performed first-principles Monte Carlo (MC) and molecular dynamics (MD) simulations to explore these clusters in finer detail. We find that molecular dynamics is not the most appropriate technique for studying aggregation in vapor phases, even at relatively high temperatures. Specifically, our MD simulations exhibit substantial problems in adequately sampling the equilibrium cluster size distribution. In contrast, MC simulations with specialized cluster moves are able to accurately sample the phase space of hydrogen-bonding vapors. At all densities, we find a stable, slightly distorted linear H2O-Cu-Cl structure, which is in agreement with the earlier simulations, surrounded by a variable number of water molecules. The surrounding water molecules do not form a well-defined second solvation shell but rather a loose network of hydrogen-bonded water with molecular CuCl on the outside edge of the water cluster. We also find a broad distribution of hydration numbers, especially at higher densities. In contrast to previous simulation work but in agreement with experimental data, we find that the average hydration number substantially increases with increasing density. Moreover, the value of the hydration number depends on the choice of cluster definition.

Calculation of self-diffusion coefficients in supercritical carbon dioxide using mean force kinetic theory.

This paper presents an application of mean force kinetic theory (MFT) to the calculation of the self-diffusivity of CO2 in the supercritical fluid regime. Two modifications to the typical application of MFT are employed to allow its application to a system of molecular species. The first is the assumption that the inter-particle potential of mean force can be obtained from the molecule center-of-mass pair correlation function, which in the case of CO2 is the C-C pair correlation function. The second is a new definition of the Enskog factor that describes the effect of correlations at the surface of the collision volume. The new definition retains the physical picture that this quantity represents a local density increase, resulting from particle correlations, relative to that in the zero density homogeneous fluid limit. These calculations are facilitated by the calculation of pair correlation functions from molecular dynamics (MD) simulations using the FEPM2 molecular CO2 model. The self-diffusivity calculated from theory is in good agreement with that from MD simulations up to and slightly beyond the density at the location of the Frenkel line. The calculation is compared with and is found to perform similarly well to other commonly used models but has a greater potential for application to systems of mixed species and to systems of particles with long range interatomic potentials due to electrostatic interactions.

NaCl aggregation in water at elevated temperatures and pressures: Comparison of classical force fields.

The properties of water vary dramatically with temperature and density. This can be exploited to control its effectiveness as a solvent. Thus, supercritical water is of keen interest as solvent in many extraction processes. The low solubility of salts in lower density supercritical water has even been suggested as a means of desalination. The high temperatures and pressures required to reach supercritical conditions can present experimental challenges during collection of required physical property and phase equilibria data, especially in salt-containing systems. Molecular simulations have the potential to be a valuable tool for examining the behavior of solvated ions at these high temperatures and pressures. However, the accuracy of classical force fields under these conditions is unclear. We have, therefore, undertaken a parametric study of NaCl in water, comparing several salt and water models at 200 bar-600 bar and 450 K-750 K for a range of salt concentrations. We report a comparison of structural properties including ion aggregation, hydrogen bonding, density, and static dielectric constants. All of the force fields qualitatively reproduce the trends in the liquid phase density. An increase in ion aggregation with decreasing density holds true for all of the force fields. The propensity to aggregate is primarily determined by the salt force field rather than the water force field. This coincides with a decrease in the water static dielectric constant and reduced charge screening. While a decrease in the static dielectric constant with increasing NaCl concentration is consistent across all model combinations, the salt force fields that exhibit more ionic aggregation yield a slightly smaller dielectric decrement.

Thermodynamics, dynamics, and structure of supercritical water at extreme conditions.

Molecular dynamics (MD) simulations to understand the thermodynamic, dynamic, and structural changes in supercritical water across the Frenkel line and the melting line have been performed. The two-phase thermodynamic model [J. Phys. Chem. B, 2010, 114(24), 8191-8198] and the velocity autocorrelation functions are used to locate the Frenkel line and to calculate the thermodynamic and dynamic properties. The Frenkel lines obtained from the two-phase thermodynamic model and the velocity autocorrelation criterion do not agree with each other. Structural characteristics and the translational diffusion dynamics of water suggest that this inconsistency could arise from the two oscillatory modes in water, which are associated with the bending of hydrogen bonds and intermolecular collisions inside the first coordination shell. The overall results lead us to conclude that the universality of the Frenkel line as a dynamic crossover line from rigid to nonrigid fluids is preserved in water.

Universality, Scaling, and Collapse in Supercritical Fluids.

Supercritical fluid (SCF) is known to exhibit salient dynamic and thermodynamic crossovers and an inhomogeneous molecular distribution. However, the question as to what basic physics underlies these microscopic and macroscopic anomalies remains open. Here, using an order parameter extracted by machine learning, the fraction of gas-like (or liquid-like) molecules, we find simplicity and universality in SCF: First, all isotherms of a given fluid collapse onto a single master curve described by a scaling relation. The observed power law holds from the high-temperature and -pressure regime down to the critical point where it diverges. Second, phase diagrams of different compounds collapse onto their master curves by the same scaling exponent, thereby demonstrating a putative law of corresponding supercritical states in simple fluids. The reported results support a model of the SCF as a mixture of two interchangeable microstates, whose spatiotemporal dynamics gives rise to unique macroscopic properties.

Electrical conductivity, ion pairing, and ion self-diffusion in aqueous NaCl solutions at elevated temperatures and pressures.

We have performed classical molecular dynamics (MD) simulations of aqueous sodium chloride (NaCl) solutions from 298 to 674 K at 200 bars to understand the influence of ion pairing and ion self-diffusion on electrical conductivity in high-temperature/high-pressure salt solutions. Conductivity data obtained from the MD simulation highlight an apparent anomaly, namely, a conductivity maximum as temperature increases along an isobar, which has been also observed in experimental studies. By examining both velocity autocorrelation and cross-correlation terms of the Green-Kubo integral, we quantitatively demonstrate that the conductivity anomaly arises mainly from a competition between the single-ion self-diffusion and the contact ion pair formation. The velocity autocorrelation function in conjunction with structural analysis suggests that diffusive motion of ions is suppressed at high temperatures due to the persistence of an inner hydration shell. The contribution of velocity cross-correlation functions between oppositely charged ions becomes significant at the onset of the conductivity decrease. Structural analysis based on Voronoi tessellation and pair correlation functions indicates that the fraction of contact ion pairs increases as temperature increases. Spatial decomposition of the electrical conductivity also indicates that the formation of contact ion pairs significantly decreases the electrical conductivity compared to Nernst-Einstein conductivity, but the contribution of distant opposite charges cannot be ignored except at the highest temperature due to unscreened long-range interactions.

Topological extension of the isomorph theory based on the Shannon entropy.

Isomorph theory is one of the promising theories for understanding the quasiuniversal relationship between thermodynamic, dynamic, and structural characteristics. Based on the hidden scale invariance of the inverse power law potentials, it rationalizes the excess entropy scaling law of dynamic properties. This work aims to show that this basic idea of isomorph theory can be extended by examining the microstructural features of the system. Using the topological framework in conjunction with the entropy calculation algorithm, we demonstrate that Voronoi entropy, a measure of the topological diversity of single atoms, provides a scaling law for the transport properties of soft-sphere fluids, which is comparable to the frequently used excess entropy scaling. By examining the relationship between the Voronoi entropy and the solidlike fraction of simple fluids, we suggest that the Frenkel line, a rigid-nonrigid crossover line, be a topological isomorphic line at which the scaling relation qualitatively changes.

Topological generalization of the rigid-nonrigid transition in soft-sphere and hard-sphere fluids.

A fluid particle changes its dynamics from diffusive to oscillatory as the system density increases up to the melting density. Hence the notion of the Frenkel line was introduced to demarcate the fluid region into rigid and nonrigid liquid subregions based on the collective particle dynamics. In this work, we apply a topological framework to locate the Frenkel lines of the soft-sphere and the hard-sphere models relying on the system configurations. The topological characteristics of the ideal gas and the maximally random jammed state are first analyzed, then the classification scheme designed in our earlier work is applied to soft-sphere and hard-sphere fluids. The dependence of the classification result on the bulk density is understood based on the theory of fluid polyamorphism. The percolation behavior of solid-like clusters is described based on the fraction of solid-like molecules in an integrated manner. The crossover densities are obtained by examining the percolation of solid-like clusters. The resultant crossover densities of soft-sphere fluids converge to that of hard-sphere fluid. Hence the topological method successfully highlights the generality of the Frenkel line.

A corresponding-state framework for the structural transition of supercritical fluids across the Widom delta.

This work proposes a classification algorithm based on the radical Voronoi tessellation to define the Widom delta, supercritical gas-liquid coexistence region, of polyatomic molecules. Specifically, we use a weighted mean-field classification method to classify a molecule into either gaslike or liquidlike. Classical percolation theory methods are adopted to understand the generality of the structural transition and to locate the Widom delta. A structural analysis on various supercritical fluids shows that the proposed method detects the influence of the attractive interaction on the structural transition of supercritical fluids. Moreover, we demonstrate that the supercritical gas-liquid coexistence region of water overlaps with the ridges of the response function maxima. From the pressure-temperature relation, a three-parameter corresponding state theorem is derived, which states that the fraction of gaslike molecules of a substance is equal to that of another if their reduced pressure, reduced temperature, and the critical compressibility factor are the same.

Topological Characterization of Rigid-Nonrigid Transition across the Frenkel Line.

The dynamics of supercritical fluids, a state of matter beyond the gas-liquid critical point, changes from diffusive to oscillatory motions at high pressure. This transition is believed to occur across a locus of thermodynamic states called the Frenkel line. The Frenkel line has been extensively investigated from the viewpoint of the dynamics, but its structural meaning is still not well-understood. This Letter interprets the mesoscopic picture of the Frenkel line entirely based on a topological and geometrical framework. This discovery makes it possible to understand the mechanism of rigid-nonrigid transition based not on the dynamics of individual atoms but on their instantaneous configurations. The topological classification method reveals that the percolation of solid-like structures occurs above the rigid-nonrigid crossover densities.

Fabrication and Characterization of Multiscale PLA Structures Using Integrated Rapid Prototyping and Gas Foaming Technologies.

Multiscale structured polymers have been considered as a promising category of functional materials with unique properties. We combined rapid prototyping and gas foaming technologies to fabricate multiscale functional materials of superior mechanical and thermal insulation properties. Through scanning electron microscope based morphological characterization, formation of multiscale porous structure with nanoscale cellular pores was confirmed. Improvement in mechanical strength is attributed to rearrangement of crystals within CO2 saturated grid sample. It is also shown that a post-foaming temperature higher than the glass transition temperature deteriorates mechanical strength, providing process guidelines. Thermal decomposition of filament material sets the upper limit of temperature for 3D printed features, characterized by simultaneous differential scanning calorimetry and thermogravimetric analysis. Porosity of the fabricated 3D structured polylactic acid (PLA) foam is controllable by suitable tuning of foaming conditions. The fabricated multiscale 3D structures have potential for thermal insulation applications with lightweight and reasonable mechanical strength.

"Two-Phase" Thermodynamics of the Frenkel Line.

The Frenkel line, a crossover line between rigid and nonrigid dynamics of fluid particles, has recently been the subject of intense debate regarding its relevance as a partitioning line of the supercritical phase, where the main criticism comes from the theoretical treatment of collective particle dynamics. From an independent point of view, this Letter suggests that the two-phase thermodynamics model may alleviate this contentious situation. The model offers new criteria for defining the Frenkel line in the supercritical region and builds a robust connection among the preexisting, seemingly inconsistent definitions. In addition, one of the dynamic criteria locates the rigid-nonrigid transition of the soft-sphere and the hard-sphere models. Hence, we suggest the Frenkel line be considered as a dynamic rigid-nonrigid fluid boundary, without any relation to gas-liquid transition. These findings provide an integrative viewpoint combining fragmentized definitions of the Frenkel line, allowing future studies to be carried out in a more reliable manner.

Probabilistic characterization of the Widom delta in supercritical region.

We present a probabilistic classification algorithm to understand the structural transition of supercritical Lennard-Jones (LJ) fluid. The classification algorithm is designed based on the exploratory data analysis on the nearest Voronoi neighbors of subcritical vapor and liquid. The algorithm is tested and applied to LJ type fluids modeled with the truncated and shifted potential and the Weeks-Chandler-Andersen potential. The algorithm makes it available to locate the Widom delta, which encloses the supercritical gas-liquid boundary and the percolation transition loci in a geometrical manner, and to conjecture the role of attractive interactions on the structural transition of supercritical fluids. Thus, the designed algorithm offers an efficient and comprehensible method to understand the phase behavior of a supercritical mesophase.

Widom Delta of Supercritical Gas-Liquid Coexistence.

Density fluctuations and the Widom line are of great importance in understanding the critical phenomena and the behaviors of supercritical fluids (SCFs). We report on the direct classification of liquid-like and gas-like molecules coexisting in the SCF, identified by machine learning analysis on simulation data. The deltoid coexistence region encloses the Widom line and may therefore be termed the Widom delta. Number fractions of gas-like and liquid-like particles are found to undergo continuous transition across the delta, following a simplified two-state model. These fractions are closely related to the magnitude of supercritical anomaly, which originates from the fluctuation between the two types. This suggests a microscopic view of the SCF as a mixture of liquid-like and gas-like structures, providing an integrative explanation to the anomalous behaviors near the critical point and the Widom line.