Prediction of Contact Angles and Density Profiles of Sessile Droplets Using Classical Density Functional Theory Based on the PCP-SAFT Equation of State

Predicting Solvation Free Energies from the Minnesota Solvation Database Using Classical Density Functional Theory Based on the PC-SAFT Equation of State

We critically assess the capabilities of classical density functional theory (DFT) based on the perturbed-chain statistical associating fluid theory (PC-SAFT) equation of state to predict the solvation free energies of small molecules in various hydrocarbon solvents. We compare DFT results with experimental data from the Minnesota solvation database and utilize statistical methods to analyze the accuracy of our approach, as well as its weaknesses. The mean absolute error of the solvation free energies is 3.7 kJ mol-1 for n-alkane solvents, ranging from pentane to hexadecane, with 473 solute-solvent systems. For solvents consisting of cyclic hydrocarbons (cyclohexane, benzene, toluene, and ethylbenzene) with 245 solute-solvent systems, we report a slightly larger mean absolute error of 4.2 kJ mol-1. We identify three possible sources of errors: (i) the neglect of solute-solvent and solvent-solvent Coulomb interactions, which limits the applicability of PC-SAFT DFT to nonpolar and weakly polar molecules; (ii) the solute's Lennard-Jones parameters supplied by the general AMBER force field, which are not parametrized toward solvation free energies; and (iii) the application of the Lorentz-Berthelot combining rules to the dispersive interactions between a segment of the PC-SAFT solvent and a Lennard-Jones interaction site of the solute. The approach is more accurate than standard implementations of phenomenological models in common chemistry software packages, which exhibit mean absolute errors larger than 9.12 kJ mol-1, even though newer phenomenological models achieve a mean absolute error of about 2 kJ mol-1. PC-SAFT DFT is more computationally efficient than state of the art explicit molecular simulations in combination with free energy perturbation methods. It is predictive with respect to solvation free energies, i.e., the input for the model is the (element-specific) molecular force field, the solute configuration from molecular dynamics simulations, and the (substance-specific) PC-SAFT parameters. The PC-SAFT parametrization uses pure-component data and does not require experimental solvation free energies. The PC-SAFT equation of state, without applying a DFT formalism, can also be used to calculate solvation free energies, provided that the PC-SAFT parameters for the solute are available. A large number of substances was recently parametrized by members of our group (Esper, T.; Bauer, G.; Rehner, P.; Gross, J. Ind. Eng. Chem. Res. 2023, 62), which enables a comparison to the DFT approach for 103 substances. Accurate results are obtained from the PC-SAFT equation of state with an MAE below 2.51 kJ mol-1. The DFT approach does not require PC-SAFT parameters for the solute and can be applied to all solutes that can be represented by the molecular force field.

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

We present a classical density functional theory (DFT) for fluid mixtures that is based on a third-order thermodynamic perturbation theory of Feynman-Hibbs-corrected Mie potentials. The DFT is developed to study the interfacial properties of hydrogen, helium, neon, deuterium, and their mixtures, i.e., fluids that are strongly influenced by quantum effects at low temperatures. White Bear fundamental measure theory is used for the hard-sphere contribution of the Helmholtz energy functional, and a weighted density approximation is used for the dispersion contribution. For mixtures, a contribution is included to account for non-additivity in the Lorentz-Berthelot combination rule. Predictions of the radial distribution function from DFT are in excellent agreement with results from molecular simulations, both for pure components and mixtures. Above the normal boiling point and 5% below the critical temperature, the DFT yields surface tensions of neon, hydrogen, and deuterium with average deviations from experiments of 7.5%, 4.4%, and 1.8%, respectively. The surface tensions of hydrogen/deuterium, para-hydrogen/helium, deuterium/helium, and hydrogen/neon mixtures are reproduced with a mean absolute error of 5.4%, 8.1%, 1.3%, and 7.5%, respectively. The surface tensions are predicted with an excellent accuracy at temperatures above 20 K. The poor accuracy below 20 K is due to the inability of Feynman-Hibbs-corrected Mie potentials to represent the real fluid behavior at these conditions, motivating the development of new intermolecular potentials. This DFT can be leveraged in the future to study confined fluids and assess the performance of porous materials for hydrogen storage and transport.

Wettability of net C, net W and net Y: a molecular dynamics simulation study

The experimental synthesis of biphenylene, a two-dimensional carbon allotrope, theoretically predicted in 1997, took place in 2021. Biphenylene is also called net C. Two close relatives of this structure, known as net W and net Y, have not yet been experimentally synthesized. In this article, the wettability properties of these three carbon allotropes are investigated, using molecular dynamics simulation. The electronic and mechanical properties of these allotropes have been extensively studied, but their wettability properties are unknown. The chemical structure of the three allotropes is similar and contain four, six, and eight carbon membered rings. The results of molecular dynamics calculations with reactive potential show that net C, net W and net Y are hydrophobic substrates with contact angles of 122.3° ± 1.3°, 126.2° ± 1.3° and 127.8° ± 1.2°, respectively. The droplets on the above-mentioned substrates have a completely layered structure. That is, the water molecules inside the droplet are completely placed in certain layers. Calculating the order parameter for water molecules shows that the degree of water molecules' tetrahedrality on all three substrates is exactly the same. In terms of hydrogen bonding at the interface, the three substrates act identically and show almost the same effect. The droplet displacement is the highest on net W and the lowest on net Y. Furthermore, the van der Waals potential on all three substrates has been scanned. It is demonstrated that the amount of droplet displacement on the surface is inversely related to the surface density of the potential peaks.

Hydrodynamic density functional theory for mixtures from a variational principle and its application to droplet coalescence

Dynamic density functional theory (DDFT) allows the description of microscopic dynamical processes on the molecular scale extending classical DFT to non-equilibrium situations. Since DDFT and DFT use the same Helmholtz energy functionals, both predict the same density profiles in thermodynamic equilibrium. We propose a molecular DDFT model, in this work also referred to as hydrodynamic DFT, for mixtures based on a variational principle that accounts for viscous forces as well as diffusive molecular transport via the generalized Maxwell-Stefan diffusion. Our work identifies a suitable expression for driving forces for molecular diffusion of inhomogeneous systems. These driving forces contain a contribution due to the interfacial tension. The hydrodynamic DFT model simplifies to the isothermal multicomponent Navier-Stokes equation in continuum situations when Helmholtz energies can be used instead of Helmholtz energy functionals, closing the gap between micro- and macroscopic scales. We show that the hydrodynamic DFT model, although not formulated in conservative form, globally satisfies the first and second law of thermodynamics. Shear viscosities and Maxwell-Stefan diffusion coefficients are predicted using an entropy scaling approach. As an example, we apply the hydrodynamic DFT model with a Helmholtz energy density functional based on the perturbed-chain statistical associating fluid theory equation of state to droplet and bubble coalescence in one dimension and analyze the influence of additional components on coalescence phenomena.

Predicting solvation free energies in non-polar solvents using classical density functional theory based on the PC-SAFT equation of state

We propose a predictive Density Functional Theory (DFT) for the calculation of solvation free energies. Our approach is based on a Helmholtz free-energy functional that is consistent with the Perturbed-Chain Statistical Associating Fluid Theory (PC-SAFT) equation of state. This allows for a coarse-grained description of the solvent based on an inhomogeneous density of PC-SAFT segments. The solute, on the other hand, is described in full detail by atomistic Lennard-Jones interaction sites. The approach is entirely predictive as it only takes the PC-SAFT parameters of the solvent and the force-field parameters of the solute as input. No adjustable parameters or empirical corrections are involved. The framework is applied to study self-solvation of n-alkanes and to the calculation of residual chemical potentials in binary solvent mixtures. Our DFT approach accurately predicts solvation free energies of small molecular solutes in three different non-polar solvents, namely n-hexane, cyclohexane, and benzene. Additionally, we show that the calculated solvation free energies agree well with those obtained by molecular dynamics simulations and with the residual chemical potential calculated by the bulk PC-SAFT equation of state. We observe higher deviations for the solvation free energy of systems with significant solute-solvent Coulomb interactions.

Free-Energy-Averaged Potentials for Adsorption in Heterogeneous Slit Pores Using PC-SAFT Classical Density Functional Theory

This study analyzes the adsorption behavior in two-dimensional heterogeneous slit pores using nonlocal density functional theory based on the perturbed-chain statistical associating fluid theory (PC-SAFT) equation of state. Both chemical heterogeneity and surface roughness on small atomistic scales are investigated. The solid structure is considered as individual solid interaction sites whereby chemical heterogeneity is introduced through the presence of different solid-fluid sites and molecular roughness by varying the position of the interaction sites in the first solid layers. The effect of both forms of heterogeneity on the adsorption behavior is assessed individually. Effective one-dimensional solid-fluid potentials provide a way to reduce the dimensionality and computational demand of the density functional theory (DFT) calculations. We determine one-dimensional free-energy-averaged (FEA) solid-fluid potentials of methane and n-butane in the low-density limit for solid systems with molecular roughness and chemical heterogeneity. Using this effective one-dimensional solid-fluid potential at any density, we find excellent agreement of adsorption isotherms for both solid descriptions in systems with homogeneous slit pores. Subcritical adsorption isotherms of n-butane in slit pores with surface roughness show deviations at higher pressures due to the formation of fluid layers in the one-dimensional FEA potential. Chemical heterogeneity introduces a shift of the capillary condensation pressure below the saturation pressure of the bulk liquid, which is well described by the free-energy-averaged system.

Prediction of Adsorption Isotherms and Selectivities: Comparison between Classical Density Functional Theory Based on the Perturbed-Chain Statistical Associating Fluid Theory Equation of State and Ideal Adsorbed Solution Theory

This study gives an assessment of the predictive capability of classical density functional theory (DFT) for adsorption processes of pure substances and mixtures of spherical and nonspherical molecular species. A Helmholtz energy functional based on the perturbed-chain statistical associating fluid theory (PC-SAFT) is applied to calculate isotherms and selectivities of multicomponent adsorption. In order to unambiguously assess the accuracy of the DFT model, we conduct molecular simulations. Monte Carlo (MC) simulations are performed in the grand canonical ensemble using the transition matrix. Two types of systems are studied: a model system, where fluid-fluid and solid-fluid interactions are defined as (single-site) Lennard-Jones interactions, and a more realistic methane-n-butane mixture in a graphite-like pore. Differences between a slit-shaped and a cylindrical pore geometry are examined for the model system. Adsorption isotherms and selectivities obtained from DFT calculations and MC simulations are found in very good agreement, particularly at high pressures. Capillary condensation observed along adsorption isotherms containing n-butane was accurately predicted, both, in equilibrium pressure and in density-increase. Comparisons with results from the ideal adsorbed solution theory are presented, confirming powerful predictions of the DFT approach.

A heuristic approach for nanodrops on a smooth solid surface

A heuristic approach is developed to obtain a simple equation for the contact angle of a nanodrop on a smooth planar solid surface. First, nanodrops of various fluids in contact with various solid surfaces are considered on the basis of nonlocal density functional theory (DFT). Along with the traditional (apparent) contact angle, θ_a, which the drop profile makes with the solid surface, another one, θ_d, formed by the smooth part of the drop profile and the horizontal plane separating that part from the oscillatory part of the profile was examined. For each of the contact angles, a separate simple equation resembling the Young equation for the macroscopic drops but containing, instead of surface tensions, the microscopic parameters of intermolecular interactions, temperature, and average density of the fluid was hypothesized and the parameters of this equation were determined using the results of DFT calculations. It was shown that predictions of these equations coincide with the results provided by DFT.

An analog to Bond number for pendant nanodrops

A new dimensionless number is introduced which characterizes the shape and stability of a pendant nanodrop.