Some formal properties of ensemble density functionals

Formal properties of ensemble density functionals are examined. Expressions for the difference between energy functionals where the particle number differs by one are constructed in terms of their first functional derivatives for the universal energy functional, the electron–electron repulsion energy functional, and the interacting kinetic energy functional. Equations that must be satisfied by second and higher order functional derivatives are derived. It is also shown that the shape of ${\delta V_{ee}[\rho]\over\delta\rho({\bf r})}$ and ${\delta K[\rho]\over\delta\rho({\bf r})}$ , the functional derivatives of the mutual electron–electron repulsion, and kinetic energy, respectively, are separately particle number independent for particle numbers between successive integers. © 2013 Wiley Periodicals, Inc.     

Exact on-event expressions for discrete potential systems

The properties of systems composed of atoms interacting though discrete potentials are dictated by a series of events which occur between pairs of atoms. There are only four basic event types for pairwise discrete potentials and the square-well/shoulder systems studied here exhibit them all. Closed analytical expressions are derived for the on-event kinetic energy distribution functions for an atom, which are distinct from the Maxwell-Boltzmann distribution function. Exact expressions are derived that directly relate the pressure and temperature of equilibrium discrete potential systems to the rates of each type of event. The pressure can be determined from knowledge of only the rate of core and bounce events. The temperature is given by the ratio of the number of bounce events to the number of disassociation/association events. All these expressions are validated with event-driven molecular dynamics simulations and agree with the data within the statistical precision of the simulations.     

Exact density functionals for two-electron systems in an external magnetic field

In principle, the extension of density functional theory (DFT) to Coulombic systems in a nonvanishing magnetic field is via current DFT (CDFT). Though CDFT is long established formally, relatively little is known with respect to any generally applicable, reliable approximate E(XC) and A(XC) functionals analogous with the workhorse approximate functionals (local density approximation and generalized gradient approximation) of ordinary DFT. Progress can be aided by having benchmark studies on a solvable correlated system. At zero field, the best-known finite system for such purposes is Hooke's atom. Recently we extended the exact ground state solutions for this two-electron system to certain combinations of nonzero external magnetic fields and confinement strengths. From those exact solutions, as well as high-accuracy numerical results for other field and confinement combinations, we construct the correlated electron density and paramagnetic current density, the exact Kohn-Sham orbitals, and the exact DFT and CDFT exchange-correlation energies and potentials. We compare with results from several widely used approximate functionals, all of which exhibit major qualitative failures, whether in CDFT or in naive application of ordinary DFT. We also illustrate how the CDFT vorticity variable nu is a computationally difficult quantity which may not be appropriate in practice to describe the external B field effects on E(XC) and A(XC).     

Density functionals for exchange and correlation energies: Exact conditions and comparison of approximations

A set of exact conditions is compiled for the purpose of developing and testing approximations for the exchange-correlation energy as a functional of the electron density. Special emphasis is placed upon recently developed density-scaling relationships. Commonly used generalized gradient approximations are compared against several of these conditions. A direct tabular comparison of these functionals (not of calculated properties) with one another is also made. © 1994 John Wiley & Sons, Inc.     

Exact analytical expressions for the potential of electrical double layer interactions for a sphere-plate system

Exact, closed-form analytical expressions are presented for evaluating the potential energy of electrical double layer (EDL) interactions between a sphere and an infinite flat plate for three different types of interactions: constant potential, constant charge, and an intermediate case as given by the linear superposition approximation (LSA). By taking advantage of the simpler sphere-plate geometry, simplifying assumptions used in the original Derjaguin approximation (DA) for sphere-sphere interaction are avoided, yielding expressions that are more accurate and applicable over the full range of κa. These analytical expressions are significant improvements over the existing equations in the literature that are valid only for large κa because the new equations facilitate the modeling of EDL interactions between nanoscale particles and surfaces over a wide range of ionic strength.     

Single particle versus ensemble average: from power-law intermittency of a single quantum dot to quasistretched exponential fluorescence decay of an ensemble

Light-induced diffusion-controlled electron transfer is proposed as an underlying mechanism for the intermittency (power law and breakdown) of a single quantum dot and ensemble-averaged fluorescence decay. The intensity decay can be approximated to a stretched exponential expression. The physical links to the free energy gap, reorganization energy, electronic coupling, and diffusion correlation times are discussed. A procedure is described for extracting these molecular-based parameters from experiments and is demonstrated with examples using existing data.     

Comparison of the accuracy of aerosol refractive index measurements from single particle and ensemble techniques

The ability of two techniques, aerosol cavity ring down spectroscopy (A-CRDS) and optical tweezers, to retrieve the refractive index of atmospherically relevant aerosol was compared through analysis of supersaturated sodium nitrate at a range of relative humidities. Accumulation mode particles in the diameter range 300-600 nm were probed using A-CRDS, with optical tweezer measurements performed on coarse mode particles several micrometers in diameter. A correction for doubly charged particles was applied in the A-CRDS measurements. Both techniques were found to retrieve refractive indices in good agreement with previously published results from Tang and Munkelwitz, with a precision of ±0.0012 for the optical tweezers and ±0.02 for the A-CRDS technique. The coarse mode optical tweezer measurements agreed most closely with refractive index predictions made using a mass-weighted linear mixing rule. The uncertainty in the refractive index retrieved by the A-CRDS technique prevented discrimination between predictions using both mass-weighted and volume-weighted linear mixing rules. No efflorescence or kinetic limitations on water transport between the particle and the gas phase were observed at relative humidities down to 14%. The magnitude of the uncertainty in refractive index retrieved using the A-CRDS technique reflects the challenges in determining particle optical properties in the accumulation mode, where the extinction efficiency varies steeply with particle size.     

Approximate analytic expressions for the diffusiophoretic velocity of a spherical colloidal particle

The general expression is derived for the diffusiophoretic velocity of a spherical colloidal particle of radius a in a concentration gradient of symmetrical electrolyte. On the basis of this expression, simple approximate analytic expressions for the diffusiophoretic velocity correct up to the order of 1/κa is derived, where κ is the Debye-Hückel parameter. It is found that the approximate expression correct to order unity can be applied for κa ≥ 50 with negligible errors, while the approximate expression correct to order 1/κa can be applied for κa ≥ 20 with negligible errors.     

Extension of the cumulant definition for low particle number systems

Cumulants are usually interpreted as the connected components of density matrices, but this interpretation fails and practical problems arise when the rank n of cumulants is larger than the number of particles (N) in the system. In that case, cumulants defined in the traditional way become disconnected. To solve this problem, the definition of cumulants is extended by introducing a simple modulation factor. The modified cumulants reduce to the conventional definition, but they vanish when N < n. Using the modified definition also eliminates the error in the approximation of density matrices by low‐rank cumulants, when N < n. The problem assumes a slightly different form when we work with active space–based theories, and it can be solved by a similar approach. Another problem with cumulants, due to spin coupling [Herbert, Int J Quantum Chem 2007, 107, 703], can be solved via the introduction of a similar modulation factor. A related yet more serious issue, termed as the local particle number constraint problem, is also discussed. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Diminished gradient dependence of density functionals: constraint satisfaction and self-interaction correction

The Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation for the exchange-correlation energy functional has two nonempirical constructions, based on satisfaction of universal exact constraints on the hole density or on the energy. We show here that, by identifying one possible free parameter in exchange and a second in correlation, we can continue to satisfy these constraints while diminishing the gradient dependence almost to zero (i.e., almost recovering the local spin density approximation or LSDA). This points out the important role played by the Perdew-Wang 1991 nonempirical hole construction in shaping PBE and later constructions. Only the undiminished PBE is good for atoms and molecules, for reasons we present, but a somewhat diminished PBE could be useful for solids; in particular, the surface energies of solids could be improved. Even for atoms and molecules, a strongly diminished PBE works well when combined with a scaled-down self-interaction correction (although perhaps not significantly better than LSDA). This shows that the undiminished gradient dependence of PBE and related functionals works somewhat like a scaled-down self-interaction correction to LSDA.     

Exact expressions of mean first-passage times and splitting probabilities for random walks in bounded rectangular domains

We give exact and explicit expressions of mean first-passage times for random walks in a rectangular domain in both cases of reflecting boundary conditions and periodic boundary conditions. The situations with one or two absorbing targets are considered.     

Proof of the generalized expressions for the number of perfect matchings of polycube graphs

The number of perfect rnatchings for the linear 2 × 2 ×n cubic lattice was analytically derived by diagonalizing the skew—symmetric 4n × 4n determinant, whose non—zero off—diagonal elements are either ±1 or ±i (pure imaginary number). The basic formulation invoking the matrix manipulation follows that of Kasteleyn, but the result obtained in this paper is the first example of the analytical solution for a special case of the three-dimensional Ising model.received by the Publisher 20 September 1989     

Simulating prescribed particle densities in the grand canonical ensemble using iterative algorithms

We present two efficient iterative Monte Carlo algorithms in the grand canonical ensemble with which the chemical potentials corresponding to prescribed (targeted) partial densities can be determined. The first algorithm works by always using the targeted densities in the kT log(rho(i)) (ideal gas) terms and updating the excess chemical potentials from the previous iteration. The second algorithm extrapolates the chemical potentials in the next iteration from the results of the previous iteration using a first order series expansion of the densities. The coefficients of the series, the derivatives of the densities with respect to the chemical potentials, are obtained from the simulations by fluctuation formulas. The convergence of this procedure is shown for the examples of a homogeneous Lennard-Jones mixture and a NaCl-CaCl(2) electrolyte mixture in the primitive model. The methods are quite robust under the conditions investigated. The first algorithm is less sensitive to initial conditions.     

Legendre-transform functionals for spin-density-functional theory

We provide a rigorous proof that the Hohenberg-Kohn theorem holds for spin densities by extending Lieb's Legendre-transform formulation to spin densities. The resulting spin-density-functional theory resolves several troublesome issues. Most importantly, the present paper provides an explicit construction for the spin potentials at any point along the adiabatic connection curve, thus providing a formal basis for the use of exchange-correlation functionals of the spin density in the Kohn-Sham density-functional theory (DFT). The practical implications of this result for unrestricted Kohn-Sham DFT calculations is considered, and the existence of holes below the Fermi level is discussed. We argue that an orbital's energy tends to increase as its occupation number increases, which provides the basis for a computational algorithm for determining the occupation numbers in Kohn-Sham DFT and helps explain the origin of Hund's rules and holes below the Fermi level.     

Improved density functionals for water

The accuracy of existing density functional methods for describing the noncovalent interaction energies in small water clusters is investigated by testing 25 density functionals against a data set of 28 water dimers and 8 water trimers whose structures are taken from the literature and from simulations. The most accurate functionals are found to be PW6B95 with a mean unsigned error of 0.13 kcal/mol and MPWB1K and B98 with mean unsigned errors of 0.15 kcal/mol; the best functional with no Hartree-Fock exchange is mPWLYP, which is a GGA with a mean unsigned error of 0.28 kcal/mol. In comparison, the most popular GGA functionals, PBE and BLYP, have mean unsigned errors of 0.52 and 1.03 kcal/mol, respectively. Since GGAs are very cost efficient for both condensed-phase simulations and electronic structure calculations on large systems, we optimized four new GGAs for water. The best of these, PBE1W and MPWLYP1W, have mean unsigned errors of 0.12 and 0.17 kcal/mol, respectively. These new functionals are well suited for use in condensed-phase simulations of water and ice.