Density functional theory of inhomogeneous liquids. II. A fundamental measure approach

Previously, it has been shown that the direct correlation function for a Lennard-Jones fluid could be modeled by a sum of that for hard-spheres, a mean-field tail, and a simple linear correction in the core region constructed so as to reproduce the (known) bulk equation of state of the fluid [Lutsko, J. Chem. Phys. 127, 054701 (2007)]. Here, this model is combined with ideas from the fundamental measure theory to construct a density functional theory for the free energy. The theory is shown to accurately describe a range of inhomogeneous conditions including the liquid vapor interface, the fluid in contact with a hard wall, and a fluid confined in a slit pore. The theory gives quantitatively accurate predictions for the surface tension, including its dependence on the potential cutoff. It also obeys two important exact conditions: That relating the direct correlation function to the functional derivative of the free energy with respect to density and the wall theorem.     

Density functional theory of inhomogeneous liquids. I. The liquid-vapor interface in Lennard-Jones fluids

A simple model is proposed for the direct correlation function (DCF) for simple fluids consisting of a hard-core contribution, a simple parametrized core correction, and a mean-field tail. The model requires as input only the free energy of the homogeneous fluid, obtained, e.g., from thermodynamic perturbation theory. Comparison to the DCF obtained from simulation of a Lennard-Jones fluid shows this to be a surprisingly good approximation for a wide range of densities. The model is used to construct a density functional theory for inhomogeneous fluids which is applied to the problem of calculating the surface tension of the liquid-vapor interface. The numerical values found are in good agreement with simulation.     

Density functional theory for inhomogeneous associating chain fluids

We propose a nonlocal density functional theory for associating chain molecules. The chains are modeled as tangent spheres, which interact via Lennard-Jones (12,6) attractive interactions. A selected segment contains additional, short-ranged, highly directional interaction sites. The theory incorporates an accurate treatment of the chain molecules via the intramolecular potential formalism and should accurately describe systems with strongly varying external fields, e.g., attractive walls. Within our approach we investigate the structure of the liquid-vapor interface and capillary condensation of a simple model of associating chains with only one associating site placed on the first segment. In general, the properties of inhomogeneous associating chains depend on the association energy. Similar to the bulk systems we find the behavior of associating chains of a given length to be in between that for the nonassociating chains of the same length and that for the nonassociating chains twice as large.     

Crystal growth and classical nucleation theory

Calculations were performed of the crystal growth rates in lithium disilicate glass in the low-temperature regime where homogeneous nucleation is observed. The computations were executed using the gain-loss (Becker–Doring) equations that form the framework of Classical Nucleation Theory (CNT). The growth rates were obtained in several different ways, using various choices for the kinetic model, the generalized diffusion coefficient, and the physical input data. The results of these calculations are compared with recently obtained experimental values of the growth rates.     

Exploiting classical nucleation theory for reverse self-assembly

In this paper we introduce a new method to design interparticle interactions to target arbitrary crystal structures via the process of self-assembly. We show that it is possible to exploit the slope of the crystal nucleation free-energy barrier to sample and select optimal interparticle interactions for self-assembly into a desired structure. We apply this method to find interactions to target two simple crystal structures: a crystal with simple cubic symmetry and a two-dimensional plane with square symmetry embedded in a three-dimensional space. Finally, we discuss the potential and limits of our method and propose a general model by which a functionally infinite number of different interaction geometries may be constructed and to which our reverse self-assembly method could in principle be applied.     

Test of classical nucleation theory via molecular-dynamics simulation

A direct test of classical nucleation theory (CNT) is made using molecular-dynamics simulations. The relation between critical nucleus size and undercooling temperature is extracted and the result yields the solid-liquid interfacial energy. It is shown that the CNT, within the assumptions made for spherical nucleus in supercooled liquid, is valid in the critical regime of nucleation for a large range of undercooling and nucleus size.     

Electrochemical nucleation and the classical theory: Overpotential and temperature dependence of the nucleation rate

A comparison between the activation energies for critical nuclei formation obtained from both the classical nucleation theory and Arrhenius kinetic analysis has been proposed. Results obtained indicate that such comparison is hampered by temperature-induced changes of the surface energies, modifying the contact angles and correspondingly the activation energies for nucleation. Published in Russian in Elektrokhimiya, 2008, Vol. 44, No. 6, pp. 704–711. The text was submitted by the authors in English.     

Simple correction to the classical theory of homogeneous nucleation

Formation of the new disperse phase via homogeneous nucleation plays a fundamental role wherever the first-order phase transitions occur. Inconsistent temperature dependence of the nucleation rates and poor agreement of theoretical critical supersaturations with experimental data for a number of substances are fundamental problems of the classical nucleation theory (CNT). Here we show that these problems can be solved with a simple empirical correction to CNT. Despite its simplicity, the corrected CNT (CCNT) accurately predicts temperature dependences and absolute values of the critical supersaturations for both organic and inorganic substances with widely varying properties at different ambient conditions and it works surprisingly well in a wide size range down to few molecules. The difference in predictions of CCNT and other versions of the classical nucleation theory commonly used in analyzing experimental data is discussed. It has been found that CCNT consistently gives better agreement with experimental data than other versions of classical nucleation theory.     

Homogeneous nucleation of n-nonane and n-propanol mixtures: a comparison of classical nucleation theory and experiments

The homogeneous nucleation rates for n-nonane-n-propanol vapor mixtures have been calculated as a function of vapor-phase activities at 230 K using the classical nucleation theory (CNT) with both rigorous and approximate kinetic prefactors and compared to previously reported experimental data. The predicted nucleation rates resemble qualitatively the experimental results for low n-nonane gas phase activity. On the high nonane activity side the theoretical nucleation rates are about three orders of magnitude lower than the experimental data when using the CNT with the approximate kinetics. The accurate kinetics improves the situation by reducing the difference between theory and experiments to two orders of magnitude. Besides the nucleation rate comparison and the experimental and predicted onset activities, the critical cluster composition is presented. The total number of molecules is approximated by CNT with reasonable accuracy. Overall, the classical nucleation theory with rigorous kinetic prefactor seems to perform better. The thermodynamic parameters needed to calculate the nucleation rates are revised extensively. Up-to-date estimates of liquid phase activities using universal functional activity coefficient Dortmund method are presented together with the experimental values of surface tensions obtained in the present study.     

Density functional theory and molecular clusters

Density functional theory (DFT) methods, including nonlocal density gradient terms in the exchange and correlation energy functionals, were applied to various types of molecular clusters: H-bonded, ionic, electrostatic, and London. Reliable results on the structure and stabilization energy were obtained for the first two types of cluster as long as Becke3LYP and Becke3P86 functionals and basis sets of at least DZ + P quality were used. DFT methods with currently available functionals failed completely, however, for London-type clusters, for which no minimum was found on the potential energy surfaces. DFT interaction energy exhibits the same basis set extension dependence as the Hartree-Fock (HF) interaction energy. Therefore, the Boys-Bernardi function counterpoise procedure should be employed for elimination of the DFT basis set superposition error. © 1995 John Wiley & Sons, Inc.     

Comprehensive theory for star-like polymer micelles; combining classical nucleation and polymer brush theory

A comprehensive theory is proposed that combines classical nucleation and polymer brush theory to describe star-like polymer micelles. With a minimum of adjustable parameters, the model predicts properties such as critical micelle concentrations and micellar size distributions. The validity of the present theory is evidenced in direct comparison to experiments; this revealed that the proportionality constant in the Daoud-Cotton model is of the order of unity and that the star-limit is valid down to relatively short corona chains. Furthermore, we show that the predicted saddle points in the free energy correspond to those solutions that are accessible with self-consistent field methods for self-assembly.     

Test of classical nucleation theory on deeply supercooled high-pressure simulated silica

We test classical nucleation theory (CNT) in the case of simulations of deeply supercooled, high density liquid silica, as modeled by the van Beest-Kramer-van Santen potential. We find that at density rho=4.38 gcm(3), spontaneous nucleation of crystalline stishovite occurs in conventional molecular dynamics simulations at temperature T=3000 K, and we evaluate the nucleation rate J directly at this T via "brute force" sampling of nucleation events in numerous independent runs. We then use parallel, constrained Monte Carlo simulations to evaluate DeltaG(n), the free energy to form a crystalline embryo containing n silicon atoms, at T=3000, 3100, 3200, and 3300 K. By comparing the form of DeltaG(n) to CNT, we test the ability of CNT to reproduce the observed behavior as we approach the regime where spontaneous nucleation occurs on simulation time scales. We find that the prediction of CNT for the n dependence of DeltaG(n) fits reasonably well to the data at all T studied. Deltamu, the chemical potential difference between bulk liquid and stishovite, is evaluated as a fit parameter in our analysis of the form of DeltaG(n). Compared to directly determined values of Deltamu extracted from previous work, the fitted values agree only at T=3300 K; at lower T the fitted values increasingly overestimate Deltamu as T decreases. We find that n(*), the size of the critical nucleus, is approximately ten silicon atoms at T=3300 K. At 3000 K, n(*) decreases to approximately 3, and at such small sizes methodological challenges arise in the evaluation of DeltaG(n) when using standard techniques; indeed even the thermodynamic stability of the supercooled liquid comes into question under these conditions. We therefore present a modified approach that permits an estimation of DeltaG(n) at 3000 K. Finally, we directly evaluate at T=3000 K the kinetic prefactors in the CNT expression for J, and find physically reasonable values; e.g., the diffusion length that Si atoms must travel in order to move from the liquid to the crystal embryo is approximately 0.2 nm. We are thereby able to compare the results for J at 3000 K obtained both directly and based on CNT, and find that they agree within an order of magnitude. In sum, our work quantifies how certain predictions of CNT (e.g., for Deltamu) break down in this deeply supercooled limit, while others [the n dependence of DeltaG(n)] are not as adversely affected.     

Density functional theory investigation of decamethyldizincocene

A density functional investigation into the structure and vibrational properties of the recently synthesized, novel, Zn(I)-containing species decamethyldizincocene has been performed. Our analysis is in agreement with the general structural properties of the experimental results. We have corroborated the experimental geometry as a true minimum on the global molecular energy surface, confirmed the experimental hypothesis that the Zn atoms are in a Zn(I) state, and provided a detailed analysis of the experimentally undefined Zn-dominant IR and Raman spectral bands of this unusual Zn(I) species.     

Density functional theory: Approximating quasiparticle density functional

A constructive approach for deriving the approximating quasiparticle energy density functional is proposed. As a matter of fact, the proposed approach is the direct development of the Kohn–Sham quasiparticle concept and the Levy–Valone approach. The approach presented takes into account a pseudopotential character of the exchange-correlation part of the density functional and results in a system of functional equations to obtain ground-state energies of many-electron systems.     

Density functional theory for Yukawa fluids

We develop an approximate field theory for particles interacting with a generalized Yukawa potential. This theory improves and extends a previous splitting field theory, originally developed for counterions around a fixed charge distribution. The resulting theory bridges between the second virial approximation, which is accurate at low particle densities, and the mean-field approximation, accurate at high densities. We apply this theory to charged, screened ions in bulk solution, modeled to interact with a Yukawa potential; the theory is able to accurately reproduce the thermodynamic properties of the system over a broad range of conditions. The theory is also applied to "dressed counterions," interacting with a screened electrostatic potential, contained between charged plates. It is found to work well from the weak coupling to the strong coupling limits. The theory is able to reproduce the counterion profiles and force curves for closed and open systems obtained from Monte Carlo simulations.     

Density functional theory for triplet superconductors

The density functional theory of superconductivity is extended to triplet superconductors and superfluid helium 3. We prove a Hohenberg-Kohn-type theorem for these systems and derive effective single-particle equations. The latter include exchange and correlations in a formally exact way and allow the treatment of both electronic and phonon-induced superconductivity. The relation of this approach to the Bogolubov-de Gennes mean-field theory and to phenomenological theories based on Ginzburg-Landau functionals is discussed. © 1997 John Wiley & Sons, Inc.     

Density functional theory on phase space

Forty‐five years after the point de départ [Hohenberg and Kohn, Phys Rev, 1964, 136, B864] of density functional theory, its applications in chemistry and the study of electronic structures keep steadily growing. However, the precise form of the energy functional in terms of the electron density still eludes us—and possibly will do so forever [Schuch and Verstraete, Nat Phys, 2009, 5, 732]. In what follows we examine a formulation in the same spirit with phase space variables. The validity of Hohenberg–Kohn–Levy‐type theorems on phase space is recalled. We study the representability problem for reduced Wigner functions, and proceed to analyze properties of the new functional. Along the way, new results on states in the phase space formalism of quantum mechanics are established. Natural Wigner orbital theory is developed in depth, with the final aim of constructing accurate correlation‐exchange functionals on phase space. A new proof of the overbinding property of the Müller functional is given. This exact theory supplies its home at long last to that illustrious ancestor, the Thomas–Fermi model. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012