Intrinsic-Density Functionals

The Hohenberg-Kohn theorem and Kohn-Sham procedure are extended to functionals of the localized intrinsic density of a self-bound system such as a nucleus. After defining the intrinsic-density functional, we modify the usual Kohn-Sham procedure slightly to evaluate the mean-field approximation to the functional, and carefully describe the construction of the leading corrections for a system of fermions in one dimension with a spin-degeneracy equal to the number of particles N. Despite the fact that the corrections are complicated and nonlocal, we are able to construct a local Skyrme-like intrinsic-density functional that, while different from the exact functional, shares with it a minimum value equal to the exact ground-state energy at the exact ground-state intrinsic density, to next-to-leading order in 1/N. We briefly discuss implications for real Skyrme functionals.Comment: 15 page

Schiff Screening of Relativistic Nucleon Electric-Dipole Moments by Electrons

We show, at leading-order in the multipole expansion of the electron-nucleus interaction, that nucleon electric-dipole moments are completely shielded by electrons so that they contribute nothing to atomic electric-dipole moments, even when relativity in the nucleus is taken into account. It is well known that relativistic electron motion, by contrast, leads to dipole moments that are not screened; we discuss the reasons for the difference.Comment: 4 pages, typeset by REVTeX, submitted to PR

Self-consistent description of multipole strength: systematic calculations

We use the quasiparticle random phase approximation with a few Skyrme density functionals to calculate strength functions in the Jpi = 0+, 1-, and 2+ channels for even Ca, Ni, and Sn isotopes, from the proton drip line to the neutron drip line. We show where and how low-lying strength begins to appear as N increases. We also exhibit partial energy-weighted sums of the transition strength as functions of N for all nuclei calculated, and transition densities for many of the interesting peaks. We find that low-energy strength increases with N in all multipoles, but with distinctive features in each. The low-lying 0+ strength near the neutron at large N barely involves protons at all, with the strength coming primarily from a single two-quasineutron configuration with very large spatial extent. The low-lying 1- strength is different, with protons contributing to the transition density in the nuclear interior together with neutrons at large radii. The low-lying 2+ transition strength goes largely to more localized states. The three Skyrme interactions we test produce similar results, differing most significantly in their predictions for the location of the neutron drip line, the boundaries of deformed regions, energies of and transition strengths to the lowest 2+ states between closed shells, and isovector energy-weighted sum rules.Comment: 43 pages, 48 figures, 1 tabl

Limit on T-violating P-conserving rhoNN interaction from the gamma decay of Fe-57

We use the experimental limit on the interference of M1 and E2 multipoles in the γ decay of 57Fe to bound the time-reversal-violating parity-conserving ρNN vertex. Our approach is a large-basis shell-model calculation of the interference. We find an upper limit on the parameter g¯ρ, the relative strength of the T-violating ρNN vertex, of close to 10^(-2), a value similar to the best limits from other experiments

Quasielastic neutrino scattering from oxygen and the atmospheric neutrino problem

We examine several phenomena beyond the scope of Fermi-gas models that affect the quasielastic scattering (from oxygen) of neutrinos in the 0.1 -- 3.0 GeV range. These include Coulomb interactions of outgoing protons and leptons, a realistic finite-volume mean field, and the residual nucleon-nucleon interaction. None of these effects are accurately represented in the Monte Carlo simulations used to predict event rates due to $\mu$ and $e$ neutrinos from cosmic-ray collisions in the atmosphere. We nevertheless conclude that the neglected physics cannot account for the anomalous $\mu$ to $e$ ratio observed at Kamiokande and IMB, and is unlikely to change absolute event rates by more than 10--15\%. We briefly mention other phenomena, still to be investigated in detail, that may produce larger changes.Comment: In Revtex version 2. 14 pages, 3 figures (available on request from J. Engel, tel. 302-831-4354, [email protected]

Nuclear time-reversal violation and the Schiff moment of 225Ra

We present a comprehensive mean-field calculation of the Schiff moment of the nucleus 225Ra, the quantity which determines the static electric dipole moment of the corresponding atom if time-reversal (T) invariance is violated in the nucleus. The calculation breaks all possible intrinsic symmetries of the nuclear mean field and includes, in particular, both exchange and direct terms from the full finite-range T-violating nucleon-nucleon interaction, and the effects of short-range correlations. The resulting Schiff moment, which depends on three unknown T-violating pion-nucleon coupling constants, is much larger than in 199Hg, the isotope with the best current experimental limit on its atomic electric-dipole moment.Comment: 4 pages, 2 figures; this version (references added) to be published in PR

Effective Operators for Double-Beta Decay

We use a solvable model to examine double-beta decay, focusing on the neutrinoless mode. After examining the ways in which the neutrino propagator affects the corresponding matrix element, we address the problem of finite model-space size in shell-model calculations by projecting our exact wave functions onto a smaller subspace. We then test both traditional and more recent prescriptions for constructing effective operators in small model spaces, concluding that the usual treatment of double-beta-decay operators in realistic calculations is unable to fully account for the neglected parts of the model space. We also test the quality of the Quasiparticle Random Phase Approximation and examine a recent proposal within that framework to use two-neutrino decay to fix parameters in the Hamiltonian. The procedure eliminates the dependence of neutrinoless decay on some unfixed parameters and reduces the dependence on model-space size, though it doesn't eliminate the latter completely.Comment: 10 pages, 8 figure

Time-Reversal Violating Schiff Moment of 225Ra

We use the Skyrme-Hartree-Fock method, allowing all symmetries to be broken, to calculate the time-reversal-violating nuclear Schiff moment (which induces atomic electric dipole moments) in the octupole-deformed nucleus 225Ra. Our calculation includes several effects neglected in earlier work, including self consistency and polarization of the core by the last nucleon. We confirm that the Schiff moment is large compared to those of reflection-symmetric nuclei, though ours is generally a few times smaller than recent estimates.Comment: Typos corrected, references added, minor changesin text. Version to appear in PRC. 10 pages, 4 figure

Self-consistent Skyrme QRPA for use in axially-symmetric nuclei of arbitrary mass

We describe a new implementation of the quasiparticle random phase approximation (QRPA) in axially-symmetric deformed nuclei with Skyrme and volume-pairing energy-density functionals. After using a variety of tests to demonstrate the accuracy of the code in ^{24,26}Mg and ^{16}O, we report the first fully self-consistent application of the Skyrme QRPA to a heavy deformed nucleus, calculating strength distributions for several K^pi in ^{172}Yb. We present energy-weighted sums, properties of gamma-vibrational and low-energy K^pi=0^+ states, and the complete isovector E1 strength function. The QRPA calculation reproduces the properties of the low-lying 2^+ states as well or better than it typically does in spherical nuclei.Comment: 5 pages, 6 figure