Classical density functional theory for interfacial properties of hydrogen, helium, deuterium, neon, and their mixtures

Integrative studies of ionic liquid interface layers: bridging experiments, theoretical models and simulations

Ionic liquids (ILs) are a class of salts existing in the liquid state below 100 °C, possessing low volatility, high thermal stability as well as many highly attractive solvent and electrochemical capabilities, etc., making them highly tunable for a great variety of applications, such as lubricants, electrolytes, and soft functional materials. In many applications, ILs are first either physi- or chemisorbed on a solid surface to successively create more functional materials. The functions of ILs at solid surfaces can differ considerably from those of bulk ILs, mainly due to distinct interfacial layers with tunable structures resulting in new ionic liquid interface layer properties and enhanced performance. Due to an almost infinite number of possible combinations among the cations and anions to form ILs, the diversity of various solid surfaces, as well as different external conditions and stimuli, a detailed molecular-level understanding of their structure-property relationship is of utmost significance for a judicious design of IL-solid interfaces with appropriate properties for task-specific applications. Many experimental techniques, such as atomic force microscopy, surface force apparatus, and so on, have been used for studying the ion structuring of the IL interface layer. Molecular Dynamics simulations have been widely used to investigate the microscopic behavior of the IL interface layer. To interpret and clarify the IL structure and dynamics as well as to predict their properties, it is always beneficial to combine both experiments and simulations as close as possible. In another theoretical model development to bridge the structure and properties of the IL interface layer with performance, thermodynamic prediction & property modeling has been demonstrated as an effective tool to add the properties and function of the studied nanomaterials. Herein, we present recent findings from applying the multiscale triangle "experiment-simulation-thermodynamic modeling" in the studies of ion structuring of ILs in the vicinity of solid surfaces, as well as how it qualitatively and quantitatively correlates to the overall ILs properties, performance, and function. We introduce the most common techniques behind "experiment-simulation-thermodynamic modeling" and how they are applied for studying the IL interface layer structuring, and we highlight the possibilities of the IL interface layer structuring in applications such as lubrication and energy storage.

Revised Enskog theory for Mie fluids: Prediction of diffusion coefficients, thermal diffusion coefficients, viscosities, and thermal conductivities

Since the 1920s, the Enskog solutions to the Boltzmann equation have provided a route to predicting the transport properties of dilute gas mixtures. At higher densities, predictions have been limited to gases of hard spheres. In this work, we present a revised Enskog theory for multicomponent mixtures of Mie fluids, where the Barker-Henderson perturbation theory is used to calculate the radial distribution function at contact. With parameters of the Mie-potentials regressed to equilibrium properties, the theory is fully predictive for transport properties. The presented framework offers a link between the Mie potential and transport properties at elevated densities, giving accurate predictions for real fluids. For mixtures of noble gases, diffusion coefficients from experiments are reproduced within ±4%. For hydrogen, the predicted self-diffusion coefficient is within 10% of experimental data up to 200 MPa and at temperatures above 171 K. Binary diffusion coefficients of the CO2/CH4 mixture from simulations are reproduced within 20% at pressures up to 14.7 MPa. Except for xenon in the vicinity of the critical point, the thermal conductivity of noble gases and their mixtures is reproduced within 10% of the experimental data. For other molecules than noble gases, the temperature dependence of the thermal conductivity is under-predicted, while the density dependence appears to be correctly predicted. Predictions of the viscosity are within ±10% of the experimental data for methane, nitrogen, and argon up to 300 bar, for temperatures ranging from 233 to 523 K. At pressures up to 500 bar and temperatures from 200 to 800 K, the predictions are within ±15% of the most accurate correlation for the viscosity of air. Comparing the theory to an extensive set of measurements of thermal diffusion ratios, we find that 49% of the model predictions are within ±20% of the reported measurements. The predicted thermal diffusion factor differs by less than 15% from the simulation results of Lennard-Jones mixtures, even at densities well exceeding the critical density.

Free energy of critical droplets-from the binodal to the spinodal

Arguably, the main challenge of nucleation theory is to accurately evaluate the work of formation of a critical embryo in the new phase, which governs the nucleation rate. In Classical Nucleation Theory (CNT), this work of formation is estimated using the capillarity approximation, which relies on the value of the planar surface tension. This approximation has been blamed for the large discrepancies between predictions from CNT and experiments. In this work, we present a study of the free energy of formation of critical clusters of the Lennard-Jones fluid truncated and shifted at 2.5σ using Monte Carlo simulations, density gradient theory, and density functional theory. We find that density gradient theory and density functional theory accurately reproduce molecular simulation results for critical droplet sizes and their free energies. The capillarity approximation grossly overestimates the free energy of small droplets. The incorporation of curvature corrections up to the second order with the Helfrich expansion greatly remedies this and performs very well for most of the experimentally accessible regions. However, it is imprecise for the smallest droplets and largest metastabilities since it does not account for a vanishing nucleation barrier at the spinodal. To remedy this, we propose a scaling function that uses all relevant ingredients without adding fitting parameters. The scaling function reproduces accurately the free energy of the formation of critical droplets for the entire metastability range and all temperatures examined and deviates from density gradient theory by less than one k_BT.