
Based on its condensedmatter properties, crystal structure, and metallurgy,
which includes a phase diagram with six allotropic phases, plutonium is one of
the most complicated pure elements in its solid state. Its anomalous
properties, which are indicative of a very strongly correlated state, are
related to its special position in the periodic table, which is at the boundary
between the light actinides that have itinerant 5$f$ electrons and the heavy
actinides that have localized 5$f$ electrons. As a foundational study to probe
the role of local electronic correlations in Pu, we use the localdensity
approximation together with a continuoustime quantum Monte Carlo simulation to
investigate the electronic structure of a single Pu atom that is either
substitutionally embedded in the bulk and or adsorbed on the surface of a Th
host. This is a simpler case than the solid phases of Pu metal, which must also
include the interactions between Pu 5$f$ electrons on different Pu atoms. For
the Pu impurity atom we have found a Kondo resonance peak, which is an
important signature of electronic correlations, in the local density of states
around the Fermi energy. Furthermore, we show that the peak width of this
resonance is narrower for Pu atoms at the surface of Th than for those in the
bulk due to a weakened Pu 5$f$ligand hybridization at the surface.

An understanding of the phase diagram of elemental plutonium (Pu) must
include both the effects of the strong directional bonding and the high density
of states of the Pu 5f electrons, as well as how that bonding weakens under the
influence of strong electronic correlations. We present for the first time
electronicstructure calculations of the full 16atom per unit cell alphaphase
structure within the framework of density functional theory (DFT) together with
dynamical meanfield theory (DMFT). Our calculations demonstrate that Pu atoms
sitting on different sites within the alphaPu crystal structure have a
strongly varying site dependence of the localizationdelocalization correlation
effects of their 5f electrons and a corresponding effect on the bonding and
electronic properties of this complicated metal. In short, alphaPu has the
capacity to simultaneously have multiple degrees of electron
localization/delocalization of Pu 5f electrons within a pure singleelement
material.

Full charge selfconsistence (CSC) over the electron density has been
implemented into the local density approximation plus dynamical meanfield
theory (LDA+DMFT) scheme based on a fullpotential linear muffintin orbital
method (FPLMTO). Computational details on the construction of the electron
density from the density matrix are provided. The method is tested on the
prototypical chargetransfer insulator NiO using a simple static HartreeFock
approximation as impurity solver. The spectral and ground state properties of
bcc Fe are then addressed, by means of the spinpolarized Tmatrix fluctuation
exchange solver (SPTF). Finally the permanent magnet SmCo$_5$ is studied using
multiple impurity solvers, SPTF and Hubbard I, as the strength of the local
Coulomb interaction on the Sm and Co sites are drastically different. The
developed CSCDMFT method is shown to in general improve on materials
properties like magnetic moments, electronic structure and the materials
density.

Since their discovery nearly a decade ago, plutoniumbased superconductors
have attracted considerable interest, which is now heightened by the latest
discovery of superconductivity in PuCoIn5. In the framework of density
functional theory (DFT) within the generalized gradient approximation (GGA)
together with dynamical meanfield theory (DMFT), we present a comparative
study of the electronic structure of PuCoIn5 with the related material,
PuCoGa5. Overall, a similar GGAbased electronic structure, including the
density of states, energy dispersion, and Fermi surface topology, was found for
both compounds. The GGA Pu 5f band was narrower in PuCoIn5 than in PuCoGa5,
resulting in an effective reduction of Kondo screening in the former system, as
also shown by DMFT calculations. This phenomenon is due to the expanded lattice
for PuCoIn5.

By combining the local density approximation (LDA) with dynamical mean field
theory (DMFT), we report a systematic analysis of the spectral properties of
$\delta$plutonium with varying $5f$ occupancy. The LDA Hamiltonian is
extracted from a tightbinding (TB) fit to fullpotential linearized augmented
planewave (FPLAPW) calculations. The DMFT equations are solved by the exact
quantum Monte Carlo (QMC) method and the HubbardI approximation. We have shown
for the first time the strong sensitivity of the spectral properties to the
$5f$ occupancy, which suggests using this occupancy as a fitting parameter in
addition to the Hubbard $U$. By comparing with PES data, we conclude that the
``open shell'' $5f^{5}$ configuration gives the best agreement, resolving the
controversy over $5f$ ``open shell'' versus ``close shell'' atomic
configurations in $\delta$Pu.

Heat capacity measurements of the deltaphase stabilized alloy PuAl suggest
that strong electronphonon coupling is required to explain the moderate
renormalization of the electronic density of states near the Fermi energy. We
calculate the heat capacity contributions from the lattice and electronic
degrees of freedom as well as from the electronlattice coupling term and find
good overall agreement between experiment and theory assuming a dimensionless
electronphonon coupling parameter of order unity, lambda ~ 0.8. This large
electronphonon coupling parameter is comparable to reported values in other
superconducting metals with facecentered cubic crystal structure, for example,
Pd (lambda ~ 0.7) and Pb (lambda ~ 1.5). Further, our analysis shows evidence
of a sizable residual lowtemperature entropy contribution, S_{res} ~ 0.4 k_B
(per atom). We can fit the residual specific heat to a twolevel system.
Therefore, we speculate that the observed residual entropy originates from
crystalelectric field effects of the Pu atoms or from selfirradiation induced
defects frozen in at low temperatures.

We propose using an equationofmotion approach as an impurity solver for
dynamical mean field theory. As an illustration of this technique, we consider
a finite$U$ Hubbard model defined on the Bethe lattice with infinite
connectivity at arbitrary filling. Depending on the filling, the spectra that
is obtained exhibits a quasiparticle peak, and lower and upper Hubbard bands as
typical features of strongly correlated materials. The results are also
compared and in good agreement with exact diagonalization. We also find a
different picture of the spectral weight transfer than the iterative
perturbation theory.

The most popular electronic structure method, the linear muffintin orbital
method (LMTO), in its fullpotential (FP) and relativistic forms has been
extended to calculate the spectroscopic properties of materials form first
principles, i.e, optical spectra, xray magnetic circular dichroism (XMCD) and
magnetooptical kerr effect (MOKE). The paper describes an overview of the
FPLMTO basis set and the calculation of the momentum matrix elements. Some
applications concerning the computation of optical properties of semiconductors
and XMCD spectra of transition metal alloys are reviewed.

In a recent paper by Jones et al., it is argued, based on FPLAPW
bandstructure calculations, that previous calculations of groundstate
properties for actinides, using the FPLMTO method implemented by J. M. Wills,
are in error. We demonstrate in this paper that the conclusions of Jones et.
al. are unfounded. Calculations using the FPLMTO method are compared with
calculations performed with the FPLAPW method and, in contradiction to
statements by Jones et. al., they agree very well. The limitation of the
present implementation of the spinorbit coupling is discussed, where
especially the 6p states pose a difficulty. The observed discrepancy (010%)
between the FPLMTO calculations and the recent FPLAPW calculation of Jones
et. al. is shown to be due to the choice of muffintin radius in the
calculations. We argue that the choice of a constant muffintin radius, common
for all calculated volumes, is less satisfactory compared to the choice of a
muffintin radius that scales with the volume, especially when a large volume
interval and open packed structures are considered. The conclusion of Jones et
al. that alphaPu has partially delocalized 5f electrons is argued to be
erroneous.

The macroscopic dielectric function in the randomphaseapproximation without
local field effect has been implemented using the local density approximation
with an all electron, fullpotential linear muffintin orbital basisset. This
method is used to investigate the optical properties of the semiconductors Si,
Ge, and GaAs under hydrostatic pressure. The pressure dependence of the
effective dielectric function is compared to the experimental data of Go\~ni
and coworkers, and an excellent agreement is found when the so called
``scissorsoperator'' shift (SOS) is used to account for the correct band gap
at $\Gamma$. The effect of the $3d$ semicore states in the interband
transitions hardly changes the static dielectric function, $\epsilon_\infty$;
however, their contribution to the intensity of absorption for higher photon
energies is substantial. The spinorbit coupling has a significant effect on
$\epsilon_\infty$ of Ge and GaAs, but not of Si. The $E_1$ peak in the
dynamical dielectric function is strongly underestimated for Si, but only
slightly for Ge and GaAs, suggesting that excitonic effects might be important
only for Si.

This is a reply to J. T. Gammel's Comment (9404047) on Tuning of the CDW in
MX Chains. We show that while the harmonic elastic energy can fit a small
region around zero dimerization, it is insufficient to treat the large
amplitude CDW of the MX systems which was the main physics of our Letter (Phys.
Rev. Lett. {\bf 71}, 1415 (1993)). We also stress that \underline{anharmonic}
elastic potentials are {\it essential} to explain the chargedensitywave
systematics of an entire class of MX material.

Localdensityapproximation calculations are used to show that the
metalmetal distance along the chains controls the chargedensity wave (CDW) in
halogenbridged transitionmetal linearchain (MX) compounds. The strength of
the CDW can be understood in terms of a two band SuSchriefferHeeger model if
a hardcore ionion repulsion potential is also added. We predict a
secondorder phase transition from an insulating to a semimetallic ground
state and explain trends in dimerization, bondlength ratios, band gaps, and
Raman breathing modes in terms of the metalmetal distance. Three "ps" figures
are included at the end of the file.