Massively parallel GPU-accelerated minimization of classical density functional theory

A Jacobian-free pseudo-arclength continuation method for phase transitions in inhomogeneous thermodynamic systems

Developing phase diagrams for inhomogeneous systems in thermodynamics is difficult, in part, due to the large phase space and the possibility of unstable and metastable solutions arising from first-order phase transitions. Pseudo-arclength continuation (PAC) is a method that allows one to trace out stable and unstable solutions of nonlinear systems. Typically, PAC utilizes the Jacobian in order to implement Newton (or quasi-Newton) steps. In this work, we present a Jacobian-free PAC method that is amenable to the usual workflows in inhomogeneous thermodynamics. We demonstrate our method in systems that have first-order phase transitions, including a novel example of polyelectrolyte complex coacervation in confinement, where multiple surface phase transitions occur and can overlap with one another.

Comparative study of force-based classical density functional theory

We reexamine results obtained with the recently proposed density functional theory framework based on forces (force-DFT) [S. M. Tschopp et al., Phys. Rev. E 106, 014115 (2022)2470-004510.1103/PhysRevE.106.014115]. We compare inhomogeneous density profiles for hard sphere fluids to results from both standard density functional theory and from computer simulations. Test situations include the equilibrium hard sphere fluid adsorbed against a planar hard wall and the dynamical relaxation of hard spheres in a switched harmonic potential. The comparison to grand canonical Monte Carlo simulation profiles shows that equilibrium force-DFT alone does not improve upon results obtained with the standard Rosenfeld functional. Similar behavior holds for the relaxation dynamics, where we use our event-driven Brownian dynamics data as benchmark. Based on an appropriate linear combination of standard and force-DFT results, we investigate a simple hybrid scheme which rectifies these deficiencies in both the equilibrium and the dynamical case. We explicitly demonstrate that although the hybrid method is based on the original Rosenfeld fundamental measure functional, its performance is comparable to that of the more advanced White Bear theory.

Inside and out: surface thermodynamics from positive to negative curvature

To explore the curvature dependence of solid-fluid interfacial thermodynamics, we calculate, using Grand Canonical Monte Carlo simulation, the surface free energy γ for a 2-d hard-disk fluid confined in a circular hard container of radius R as a function of the bulk packing fraction η and wall curvature C = −1/ R. (The curvature is negative because the surface is concave) Combining this with our previous data [J. Phys. Chem. B, 124, 7938 (2020)] for the positive curvature case (a hard-disk fluid at a circular wall, C = +1/ R), we obtain a full picture of surface thermodynamics in this system over the full range of positive and negative wall curvatures. Our results show that γ is linear in C with a slope that is the same for both positive and negative wall curvatures, with deviations seen only at high negative curvatures (strong confinement) and high density. This observation indicates that the surface thermodynamics of this system is consistent with the predictions of so-called Morphometric Thermodynamics at both positive and negative curvatures. In addition, we show that classical Density Functional Theory and a generalized scaled particle theory can be constructed that give excellent agreement with the simulation data over most of the range of curvatures and densities. For extremely high curvatures, where only one or two disks can occupy the container at maximum packing, it is possible to calculate γ exactly. In this limit the simulations and DFT calculations are in remarkable agreement with the exact results.

Force density functional theory in- and out-of-equilibrium

When a fluid is subject to an external field, as is the case near an interface or under spatial confinement, then the density becomes spatially inhomogeneous. Although the one-body density provides much useful information, a higher level of resolution is provided by the two-body correlations. These give a statistical description of the internal microstructure of the fluid and enable calculation of the average interparticle force, which plays an essential role in determining both the equilibrium and dynamic properties of interacting fluids. We present a theoretical framework for the description of inhomogeneous (classical) many-body systems, based explicitly on the two-body correlation functions. By consideration of local Noether-invariance against spatial distortion of the system we demonstrate the fundamental status of the Yvon-Born-Green (YBG) equation as a local force-balance within the fluid. Using the inhomogeneous Ornstein-Zernike equation we show that the two-body correlations are density functionals and, thus, that the average interparticle force entering the YBG equation is also a functional of the one-body density. The force-based theory we develop provides an alternative to standard density functional theory for the study of inhomogeneous systems both in- and out-of-equilibrium. We compare force-based density profiles to the results of the standard potential-based (dynamical) density functional theory. In-equilibrium, we confirm both analytically and numerically that the standard approach yields profiles that are consistent with the compressibility pressure, whereas the force-density functional gives profiles consistent with the virial pressure. For both approaches we explicitly prove the hard-wall contact theorem that connects the value of the density profile at the hard-wall with the bulk pressure. The structure of the theory offers deep insights into the nature of correlation in dense and inhomogeneous systems.

Gas adsorption and separation in metal-organic frameworks by PC-SAFT based density functional theory

In this work, we examine the theoretical performance of perturbed-chain statistical associating fluid theory based density functional theory (DFT) in predicting gas adsorption and separation in metal-organic frameworks by using simulation and experimental data as the benchmark. Adsorption isotherms of methane and ethane in pure gas and mixtures and selectivities for ethane/methane mixtures are calculated. The predicted isotherms by DFT are in excellent agreement with simulation and experimental data for pure methane and ethane, whereas for the mixture, DFT is in semi-quantitative accordance with simulation results. For fast and high-throughput screening of material purpose, three algorithms including Picard iteration with line search, Anderson mixing, and Picard-Anderson-hybrid algorithm are proposed to calculate the three dimensional density distribution of confined gases. The advantages and limitations of the three algorithms under various conditions are discussed.

Direct Correlation Function of a Crystalline Solid

Direct correlation functions (DCFs), linked to the second functional derivative of the free energy with respect to the one-particle density, play a fundamental role in a statistical mechanics description of matter. This holds, in particular, for the ordered phases: DCFs contain information about the local structure including defects and encode the thermodynamic properties of crystalline solids; they open a route to the elastic constants beyond low temperature expansions. Via a demanding numerical approach, we have explicitly calculated for the first time the DCF of a solid: based on the fundamental measure concept, we provide results for the DCF of a hard sphere crystal. We demonstrate that this function differs at coexistence significantly from its liquid counterpart-both in shape as well as in its order of magnitude-because it is dominated by vacancies. We provide evidence that the traditional use of liquid DCFs in functional Taylor expansions of the free energy is conceptually wrong and show that the emergent elastic constants are in good agreement with simulation-based results.

Statics and dynamics of a finite two-dimensional colloidal system with competing attractive critical Casimir and repulsive magnetic dipole interactions

We continue our theoretical study of a recently proposed two-dimensional colloidal system with attractive critical Casimir and repulsive magnetic dipole forces that can be tuned easily and independently from each other via the temperature and the strength of an external magnetic field, respectively [K. Marolt, M. Zimmermann, and R. Roth, Phys. Rev. E 100, 052602 (2019)2470-004510.1103/PhysRevE.100.052602]. Using this freedom, it is possible to construct a competing interaction potential that causes microphase separation featuring spatially inhomogeneous cluster, stripe, and bubble phases in the bulk, i.e., in an infinite system without an external potential. In the present work, we demonstrate by means of density functional theory that microphase separation can also occur in finite geometries. In a square cell with a side length of 20 or 30 colloid diameters, we observe the emergence of highly structured cluster and ring phases at intermediate bulk densities in addition to almost uniform fluid phases for lower and higher bulk densities. We then employ dynamic density functional theory to determine how the system reacts when the temperature and the magnetic field are altered over time, and we show how to induce a transition from the liquid to the cluster/ring phase and also from the cluster directly to the ring phase. We find that often a slowly varying and nontrivial path in parameter space is required to reach a stable state, whereas abrupt changes are prone to lead to metastable configurations.

A GPU implementation of classical density functional theory for rapid prediction of gas adsorption in nanoporous materials

Nanoporous materials are promising as the next generation of absorbents for gas storage and separation with ultrahigh capacity and selectivity. The recent advent of data-driven approaches in materials modeling provides alternative routes to tailor nanoporous materials for customized applications. Typically, a data-driven model requires a large amount of training data that cannot be generated solely by experimental methods or molecular simulations. In this work, we propose an efficient implementation of classical density functional theory with a graphic processing unit (GPU) for the fast yet accurate prediction of gas adsorption isotherms in nanoporous materials. In comparison to serial computing with the central processing unit, the massively parallelized GPU implementation reduces the computational cost by more than two orders of magnitude. The proposed algorithm renders new opportunities not only for the efficient screening of a large materials database for gas adsorption but it may also serve as an important stepping stone toward the inverse design of nanoporous materials tailored to desired applications.

Single-precision open-shell CCSD and CCSD(T) calculations on graphics processing units

It has been shown that coupled-cluster calculations with single-precision data are able to provide correlation energy with insignificant loss of accuracy.

Microphase separation in a two-dimensional colloidal system with competing attractive critical Casimir and repulsive magnetic dipole interactions

We propose and study theoretically a colloidal system in two dimensions with attractive critical Casimir and repulsive magnetic dipole forces, wherein the strength of attraction and repulsion can be easily and independently tuned by adjusting the temperature and an external magnetic field, respectively. We expect this setup to be experimentally accessible and are confident that it can serve to deepen our understanding of the mechanisms behind microphase separation due to competing interactions. We develop a density functional theory for our model and present first results of our calculations in the form of a phase diagram for fixed temperature, but varying magnetic fields and bulk densities. For certain values of these parameters, we are able to confirm the existence of thermodynamically stable inhomogeneous density profiles in the bulk, such as parallel lamellar stripes, as well as clusters and voids on a hexagonal lattice.

Density functional theory for the crystallization of two-dimensional dipolar colloidal alloys

Two-dimensional mixtures of dipolar colloidal particles with different dipole moments exhibit extremely rich self-assembly behaviour and are relevant to a wide range of experimental systems, including charged and super-paramagnetic colloids at liquid interfaces. However, there is a gap in our understanding of the crystallization of these systems because existing theories such as integral equation theory and lattice sum methods can only be used to study the high temperature fluid phase and the zero-temperature crystal phase, respectively. In this paper we bridge this gap by developing a density functional theory (DFT), valid at intermediate temperatures, in order to study the crystallization of one and two-component dipolar colloidal monolayers. The theory employs a series expansion of the excess Helmholtz free energy functional, truncated at second order in the density, and taking as input highly accurate bulk fluid direct correlation functions from simulation. Although truncating the free energy at second order means that we cannot determine the freezing point accurately, our approach allows us to calculate ab initio both the density profiles of the different species and the symmetry of the final crystal structures. Our DFT predicts hexagonal crystal structures for one-component systems, and a variety of superlattice structures for two-component systems, including those with hexagonal and square symmetry, in excellent agreement with known results for these systems. The theory also provides new insights into the structure of two-component systems in the intermediate temperature regime where the small particles remain molten but the large particles are frozen on a regular lattice.

Phase diagrams and crystal-fluid surface tensions in additive and nonadditive two-dimensional binary hard-disk mixtures

Using density functionals from fundamental measure theory, phase diagrams and crystal-fluid surface tensions in additive and nonadditive (Asakura-Oosawa model) two-dimensional binary hard-disk mixtures are determined for the whole range of size ratios q=smalldiameter/largediameter, assuming random disorder (lattice points or interstitial occupied by large or small disks at random) in the crystal phase. The fluid-crystal transitions are first order due to the assumption of a periodic unit cell in the density-functional calculations. Qualitatively, the shape of the phase diagrams is similar to the case of three-dimensional hard-sphere mixtures. For the nonadditive case, a broadening of the fluid-crystal coexistence region is found for small q, whereas for large q a vapor-fluid transition intervenes. In the additive case, we find a sequence of spindle-type, azeotropic, and eutectic phase diagrams upon lowering q from 1 to 0.6. The transition from azeotropic to eutectic is different from the three-dimensional case. Surface tensions in general become smaller (up to a factor 2) upon addition of a second species and they are rather small. The minimization of the functionals proceeds without restrictions and optimized graphics card routines are used.

Nonequilibrium phase transitions of sheared colloidal microphases: Results from dynamical density functional theory

By means of classical density functional theory and its dynamical extension, we consider a colloidal fluid with spherically symmetric competing interactions, which are well known to exhibit a rich bulk phase behavior. This includes complex three-dimensional periodically ordered cluster phases such as lamellae, two-dimensional hexagonally packed cylinders, gyroid structures, or spherical micelles. While the bulk phase behavior has been studied extensively in earlier work, in this paper we focus on such structures confined between planar repulsive walls under shear flow. For sufficiently high shear rates, we observe that microphase separation can become fully suppressed. For lower shear rates, however, we find that, e.g., the gyroid structure undergoes a kinetic phase transition to a hexagonally packed cylindrical phase, which is found experimentally and theoretically in amphiphilic block copolymer systems. As such, besides the known similarities between the latter and colloidal systems regarding the equilibrium phase behavior, our work reveals further intriguing nonequilibrium relations between copolymer melts and colloidal fluids with competing interactions.

Phase behavior and bulk structural properties of a microphase former with anisotropic competing interactions: A density functional theory study

Using classical density functional theory, we investigate systems exhibiting interactions where a short-range anisotropic attractive force competes with a long-range spherically symmetric repulsive force. The former is modelled within Wertheim's first-order perturbation theory for patchy particles, and the repulsive part is assumed to be a Yukawa potential which is taken into account via a mean-field approximation. From previous studies of systems with spherically symmetric competing interactions, it is well known that such systems can exhibit stable bulk cluster phases (microphase separation) provided that the repulsion is sufficiently weak compared to the attraction. For the present model system, we find rich phase diagrams including both reentrant clustering and liquid-gas binodals. In particular, the model predicts inhomogeneous bulk phases at extremely low packing fractions, which cannot be observed in systems with isotropic competing interactions.