Energy dependence on fractional charge for strongly interacting subsystems

The energies of a pair of strongly-interacting subsystems with arbitrary noninteger charges are examined from closed and open system perspectives. An ensemble representation of the charge dependence is derived, valid at all interaction strengths. Transforming from resonance-state ionicity to ensemble charge dependence imposes physical constraints on the occupation numbers in the strong-interaction limit. For open systems, the chemical potential is evaluated using microscopic and thermodynamic models, leading to a novel correlation between ground-state charge and an electronic temperature.Comment: 4 pages, 3 figs.; as accepted (Phys. Rev. Lett.

Improved tensor-product expansions for the two-particle density matrix

We present a new density-matrix functional within the recently introduced framework for tensor-product expansions of the two-particle density matrix. It performs well both for the homogeneous electron gas as well as atoms. For the homogeneous electron gas, it performs significantly better than all previous density-matrix functionals, becoming very accurate for high densities and outperforming Hartree-Fock at metallic valence electron densities. For isolated atoms and ions, it is on a par with previous density-matrix functionals and generalized gradient approximations to density-functional theory. We also present analytic results for the correlation energy in the low density limit of the free electron gas for a broad class of such functionals.Comment: 4 pages, 2 figure

A natural orbital functional for the many-electron problem

The exchange-correlation energy in Kohn-Sham density functional theory is expressed as a functional of the electronic density and the Kohn-Sham orbitals. An alternative to Kohn-Sham theory is to express the energy as a functional of the reduced first-order density matrix or equivalently the natural orbitals. In the former approach the unknown part of the functional contains both a kinetic and a potential contribution whereas in the latter approach it contains only a potential energy and consequently has simpler scaling properties. We present an approximate, simple and parameter-free functional of the natural orbitals, based solely on scaling arguments and the near satisfaction of a sum rule. Our tests on atoms show that it yields on average more accurate energies and charge densities than the Hartree Fock method, the local density approximation and the generalized gradient approximations

Long-term monitoring and experimentalmanipulation of a Chihuahuan Desert ecosystem near Portal, Arizona, USA

Desert ecosystems have long served as model systems in the study of ecological concepts (e.g., competition, resource pulses, top-down/bottom-up dynamics). However, the inherent variability of resource availability in deserts, and hence consumer dynamics, can also make them challenging ecosystems to understand. Study of a Chihuahuan desert ecosystem near Portal, Arizona, USA, began in 1977. At this site, 24 experimental plots were established in 1977 and divided among controls and experimental manipulations. Experimental manipulations over the years include removal of all or some rodent species, all or some ants, seed additions, and various alterations of the annual plant community. While some of these manipulations were discontinued early on, others (i.e., ant and rodent manipulations) have been maintained throughout the study. Monitoring of the composition and abundances of ants, plants, and rodents has occurred continuously on all 24 plots. From 1977 to 2002, individual-level data on rodents (i.e., species, sex, size, reproductive condition) were collected monthly for each plot. From 1983 to 2002, the species-level abundances of plants were sampled on permanent quadrats. From 1977 to 2002, the species-level abundance of ant colonies was recorded for each plot, and from 1988 to 2002 additional information on ant abundances were recorded. Finally, from 1980 to 2002 we recorded precipitation at the study site. These data have been used in a variety of publications documenting the effects of the experimental manipulations as well as the response of populations and communities to longterm changes in climate and habitat. Sampling is ongoing and this database will be periodically updated

Density-matrix functional theory of the Hubbard model: An exact numerical study

A density functional theory for many-body lattice models is considered in which the single-particle density matrix is the basic variable. Eigenvalue equations are derived for solving Levy's constrained search of the interaction energy functional W, which is expressed as the sum of Hartree-Fock energy and the correlation energy E_C. Exact results are obtained for E_C of the Hubbard model on various periodic lattices. The functional dependence of E_C is analyzed by varying the number of sites, band filling and lattice structure. The infinite one-dimensional chain and one-, two-, or three-dimensional finite clusters with periodic boundary conditions are considered. The properties of E_C are discussed in the limits of weak and strong electronic correlations, as well as in the crossover region. Using an appropriate scaling we observe a pseudo-universal behavior which suggests that the correlation energy of extended systems could be obtained quite accurately from finite cluster calculations. Finally, the behavior of E_C for repulsive (U>0) and attractive (U<0) interactions are contrasted.Comment: Phys. Rev. B (1999), in pres

Interaction energy functional for lattice density functional theory: Applications to one-, two- and three-dimensional Hubbard models

The Hubbard model is investigated in the framework of lattice density functional theory (LDFT). The single-particle density matrix $\gamma_{ij}$ with respect the lattice sites is considered as the basic variable of the many-body problem. A new approximation to the interaction-energy functional $W[\gamma]$ is proposed which is based on its scaling properties and which recovers exactly the limit of strong electron correlations at half-band filling. In this way, a more accurate description of $W$ is obtained throughout the domain of representability of $\gamma_{ij}$, including the crossover from weak to strong correlations. As examples of applications results are given for the ground-state energy, charge-excitation gap, and charge susceptibility of the Hubbard model in one-, two-, and three-dimensional lattices. The performance of the method is demonstrated by comparison with available exact solutions, with numerical calculations, and with LDFT using a simpler dimer ansatz for $W$. Goals and limitations of the different approximations are discussed.Comment: 25 pages and 8 figures, submitted to Phys. Rev.

An Empirical Charge Transfer Potential with Correct Dissociation Limits

The empirical valence bond (EVB) method [J. Chem. Phys. 52, 1262 (1970)] has always embodied charge transfer processes. The mechanism of that behavior is examined here and recast for use as a new empirical potential energy surface for large-scale simulations. A two-state model is explored. The main features of the model are: (1) Explicit decomposition of the total system electron density is invoked; (2) The charge is defined through the density decomposition into constituent contributions; (3) The charge transfer behavior is controlled through the resonance energy matrix elements which cannot be ignored; and (4) A reference-state approach, similar in spirit to the EVB method, is used to define the resonance state energy contributions in terms of "knowable" quantities. With equal validity, the new potential energy can be expressed as a nonthermal ensemble average with a nonlinear but analytical charge dependence in the occupation number. Dissociation to neutral species for a gas-phase process is preserved. A variant of constrained search density functional theory is advocated as the preferred way to define an energy for a given charge.Comment: Submitted to J. Chem. Phys. 11/12/03. 14 pages, 8 figure

Thermodynamics as an alternative foundation for zero-temperature density functional theory and spin density functional theory

Thermodynamics provides a transparent definition of the free energy of density functional theory (DFT), and of its derivatives - the potentials, at finite temperatures T. By taking the T to 0 limit, it is shown here that both DFT and spin-dependent DFT (for ground states) suffer from precisely the same benign ambiguities: (a) charge and spin quantization lead to "up to a constant" indeterminacies in the potential and the magnetic field respectively, and (b) the potential in empty subspaces is undetermined but irrelevant. Surprisingly, these simple facts were inaccessible within the standard formulation, leading to recent discussions of apparent difficulties within spin-DFT.Comment: RevTeX, to appear in Phys. Rev.