Solid-state combustion synthesis of ceramics and alloys in reduced gravity

Possible microgravity effects are explored in the combustion synthesis of ceramics and alloys from their constituent elements. Molten intermediates are typically present during the combustion process, thereby offering the chance for natural convection to take place. Numerical simulations suggest that the combustion front in concert with gravity may act as a partial zone-refinement mechanism which is attempting to sweep out porosity in the sample. Contrary to suggestions by dimensional analysis, no effects on the combustion rate are seen. An analytical model of the combustion velocity as a function of the gravitational field and the spreading rate of molten material gives the correct order of magnitude of the gravity effect as measured by centrifuge experiments

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

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

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

Molecular dynamics simulations and thermochemistry of reactive ion etching of silicon by chlorine, chlorine dimer, bromine, and bromine dimer cations

Simulations of Cl plasma etch of Si surfaces with MD techniques agree reasonably well with the available experimental information on yields and surface morphologies. This information has been supplied to a Monte Carlo etch profile resulting in substantial agreement with comparable inputs provided through controlled experiments. To the extent that more recent measurements of etch rates are more reliable than older ones, preliminary MD simulations using bond-order corrections to the atomic interactions between neighboring Si atoms on the surface improves agreement with experiment through an increase in etch rate and improved agreement with XPS measurements of surface stoichiometry. Thermochemical and geometric analysis of small Si-Br molecules is consistent with the current notions of the effects of including brominated species in etchant gases

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.

Atomic interactions between plutonium and helium.

An essential issue in gallium (Ga)-stabilized fcc-phase plutonium ({delta}-Pu) is the formation of helium (He) voids and bubbles emanating from the radiolytic decay of the Pu. The rate of formation of He voids and bubbles is related to the He-defect formation energies and their associated migration barriers. The size and shape distributions of the bubbles are coupled to these critical migration processes. The values of the defect formation energies, internal pressure, and migration barriers can be estimated from atomistic calculations. Complicating this picture is the destruction of He-filled voids and bubbles by subsequent radiolytic decay events. The present study concerns the construction of the necessary potential energy surfaces for the Pu-He and He-He interactions within the modified embedded atom method (MEAM). Once fully tested, the potentials will be used to estimate the He-defect formation energies and barriers to the migration of these defects for both interstitial and substitutional He on an fcc Pu lattice. The He-He interactions are modeled from ab initio electronic structure calculations for the He{sub 2} dimer and the equilateral He, trimer. The experimental data and the electronic structure calculations on He{sub 2} agree very well. These data were fit to a Rose function fn{sub R}(x) = A P({alpha}x) exp(-{alpha}x), where P is a polynomial, x = R/R{sub 0}-1, R is the bond length, and R{sub 0} is its equilibrium value. The fits are very satisfactory. Both linear (P = 1+{alpha}x, zeroth-order Rose) and rational (P = 1+{alpha}x+a{sub 3} ({alpha}x){sup 3}/(1+x) first-order Rose) polynomials in the Rose function were tried. The more flexible rational form does improve the fit, but only marginally. Only the linear form was used thereafter. The resulting MEAM potential was used to predict the behavior of the linear trimer and the fcc cold compression curve. The results are shown in Fig. 2 and appear to be satisfactory. The compression regions of the curves are of particular interest for several reasons. First, an octahedral interstitial He atom in anfcc Pu lattice with a lattice constant of 4.64 {angstrom} has a nearest Pu neighbor distance of 2.32 {angstrom}. This distance is in the compressive region of the potential energy curve. Second, the compressive region will partially determine the internal pressure of He-filled voids and bubbles. Third, the shape of the He-filled voids will be influenced by the compression region of the potential. The Pu-He interactions are also modeled from ab initio electronic structure calculations, this time only for the PuHe dimer. The lowest-energy spin state of the dimer appears to be the S=7/2 state with a 'Stuttgart small-core RECP/6-31 g' basis. Two electronic structure methods were tried which would bound the extremes of the Pu-He interaction. One was the local density approximation (LDA), which tends to overestimate binding strength. It gives a well depth of 0.08 eV and a bond length of 3.6 {angstrom}. The other used the Becke-3-Lee-Yang-Parr (B3LYP) exchange-correlation energy functional, which tends to underestimate binding strength. It predicts no binding at any separation. For purposes of fitting to the Rose functional form, a well depth of 0.03 eV and bond length of 4.8 {angstrom} was used. The bond length exceeds the cutoff distance that will be used in future simulations to limit the maximum range of the atomic interactions and is effectively purely repulsive. Furthermore, the dimer information is sufficient to determine the painvise part of the Pu-He MEAM potential, but not the effective electron density that determines the many-body part of the potential, the embedding functions F{sub He} and F{sub Pu} in Fig. 3. The effective electron density, as well as determining which of the two dimer curves (LDA or B3LYP) is preferable, will be decided by comparing simulation results to known information about He bubble formation rates at elevated temperatures and estimates of Me bubble sizes. Initial simulations suggest that an interstitial He defect, based on either the LDA or the B3LYP dimer curve with a He:Pu density ratio of 0.04, will not remain at an octahedral site as in other fcc metals such as nickel. The He defect may also form a split interstitial with a Pu atom. The details remain to be determined

The influence of defects on magnetic properties of fcc-Pu

The influence of vacancies and interstitial atoms on magnetism in Pu has been considered in frames of the Density Functional Theory (DFT). The relaxation of crystal structure arising due to different types of defects was calculated using the molecular dynamic method with modified embedded atom model (MEAM). The LDA+U+SO (Local Density Approximation with explicit inclusion of Coulomb and spin-orbital interactions) method in matrix invariant form was applied to describe correlation effects in Pu with these types of defects. The calculations show that both vacancies and interstitials give rise to local moments in $f$-shell of Pu in good agreement with experimental data for annealed Pu. Magnetism appears due to destroying of delicate balance between spin-orbital and exchange interactions.Comment: 13 pages, 4 figure

Settling Decisions and Heterospecific Social Information Use in Shrikes

Animals often settle near competitors, a behavior known as social attraction, which belies standard habitat selection theory. Two hypotheses account for these observations: individuals obtain Allee benefits mediated by the physical presence of a competitor, or they use successfully settled individual as a source of information indicating the location of high quality habitat. We evaluated these hypotheses experimentally in two species of shrikes. These passerine birds with a raptor-like mode of life impale prey to create larders that serve as an indicator of male/habitat quality. Thus, two forms of indirect information are available in our system: a successfully settled shrike and its larder. Typically these two cues are associated with each other, however, our experimental treatment created an unnatural situation by disassociating them. We manipulated the presence of larders of great grey shrikes and examined the settling decisions of red-backed shrikes within and outside the great grey shrike territories. Male red-backed shrikes did not settle sooner on plots with great grey shrikes compared to plots that only contained artificial larders indicating that red-backed shrikes do not use the physical presence of a great grey shrike when making settling decisions which is inconsistent with the Allee effect hypothesis. In contrast, for all plots without great grey shrikes, red-backed shrikes settled, paired and laid clutches sooner on plots with larders compared to plots without larders. We conclude that red-backed shrikes use larders of great grey shrikes as a cue to rapidly assess habitat quality