Search

The thermal plasma filling the early universe generated a stochastic gravitational wave background that peaks in the microwave frequency range today. If the graviton production rate is expressed as a series in a fine-structure constant, $\alpha$, and the temperature over the Planck mass, $T^2_{ } / m_{\rm pl}^2$, then the lowest-order contributions come from single ($\sim \alpha T^2_{ }/m_{\rm pl}^2$) and double ($\sim T^4_{ }/m_{\rm pl}^4$) graviton production via $2\to 2$ scatterings. We show that in the Standard Model, single-graviton production dominates if the maximal temperature is smaller than $4\times 10^{18}_{ }$ GeV. This justifies previous calculations which relied solely on single-graviton production. We mention Beyond the Standard Model scenarios in which the single and double-graviton contributions could be of comparable magnitudes. Finally, we elaborate on what these results imply for the range of applicability of General Relativity as an effective theory., Comment: 16 pages. v2: clarifications added

The observation of neutrinoless double-beta ($0\nu\beta\beta$) decay would offer proof of lepton number violation, demonstrating that neutrinos are Majorana particles, while also helping us understand why there is more matter than antimatter in the Universe. If the decay is driven by the exchange of the three known light neutrinos, a discovery would, in addition, link the observed decay rate to the neutrino mass scale through a theoretical quantity known as the nuclear matrix element (NME). Accurate values of the NMEs for all nuclei considered for use in $0\nu\beta\beta$ experiments are therefore crucial for designing and interpreting those experiments. Here, we report the first comprehensive ab initio uncertainty quantification of the $0\nu\beta\beta$-decay NME, in the key nucleus $^{76}$Ge. Our method employs nuclear strong and weak interactions derived within chiral effective field theory and recently developed many-body emulators. Our result, with a conservative treatment of uncertainty, is an NME of $2.60^{+1.28}_{-1.36}$, which, together with the best-existing half-life sensitivity and phase-space factor, sets an upper limit for effective neutrino mass of $187^{+205}_{-62}$ meV. The result is important for designing next-generation germanium detectors aiming to cover the entire inverted hierarchy region of neutrino masses., Comment: 7 pages, 1 table, and 2 figures

This whitepaper presents the research priorities decided on by attendees of the 2022 Town Meeting for Fundamental Symmetries, Neutrons and Neutrinos, which took place December 13-15, 2022 in Chapel Hill, NC, as part of the Nuclear Science Advisory Committee (NSAC) 2023 Long Range Planning process. A total of 275 scientists registered for the meeting. The whitepaper makes a number of explicit recommendations and justifies them in detail.

This White Paper, prepared for the Fundamental Symmetries, Neutrons, and Neutrinos Town Meeting related to the 2023 Nuclear Physics Long Range Plan, makes the case for double beta decay as a critical component of the future nuclear physics program. The major experimental collaborations and many theorists have endorsed this white paper., Comment: white paper submitted for the Fundamental Symmetries, Neutrons, and Neutrinos Town Meeting in support of the US Nuclear Physics Long Range Planning Process

In numerous astrophysical scenarios, such as core-collapse supernovae and neutron star mergers, as in well as heavy-ion collision experiments, transitions between thermally populated nuclear excited states have been shown to play an important role. Due to its simplicity and excellent extrapolation ability, the finite-temperature quasiparticle random phase approximation (FT-QRPA) presents itself as an efficient method to study the properties of hot nuclei. The statistical ensembles in the FT-QRPA make the theory much richer than its zero-temperature counterpart, but also obscure the meaning of various physical quantities. In this work, we clarify several aspects of the FT-QRPA, including notations seen in the literature, and demonstrate how to extract physical quantities from the theory. To exemplify the correct treatment of finite-temperature transitions, we place special emphasis on the charge-exchange transitions described within the proton-neutron FT-QRPA (FT-PNQRPA). With the FT-PNQRPA built on the nuclear energy-density functional theory, we obtain solutions using a relativistic matrix approach and also the non-relativistic finite amplitude method. We show that the Ikeda sum rule is fulfilled with the proper treatment of de-excitations from thermally populated excited states. Additionally, we demonstrate the impact of these transitions on stellar electron capture (EC) rates in ${}^{58,78}$Ni. While their inclusion does not influence the EC rates in ${}^{58}$Ni, the rates in ${}^{78}$Ni are dominated by de-excitations for temperatures $T > 0.5$ MeV. In systems with a large negative $Q$-value, the inclusion of de-excitations within the FT-QRPA is necessary for a complete description of reaction rates at finite temperature., Comment: 19 pages, 4 figures, submitted for publication

The rapid neutron capture process (r-process) is one of the main mechanisms whereby elements heavier than iron are synthesized, and is entirely responsible for the natural production of the actinides. Kilonova emissions are modeled as being largely powered by the radioactive decay of species synthesized via the r -process. Given that the r -process occurs far from nuclear stability, unmeasured beta decay rates play an essential role in setting the time scale for the r -process. In an effort to better understand the sensitivity of kilonova modeling to different theoretical global beta-decay descriptions, we incorporate these into nucleosynthesis calculations. We compare the results of these calculations and highlight differences in kilonova nuclear energy generation and light curve predictions, as well as final abundances and their implications for nuclear cosmochronometry. We investigate scenarios where differences in beta decay rates are responsible for increased nuclear heating on time scales of days that propagates into a significantly increased average bolometric luminosity between 1-10 days post-merger. We identify key nuclei, both measured and unmeasured, whose decay rates are directly impact nuclear heating generation on timescales responsible for light curve evolution. We also find that uncertainties in beta decay rates significantly impact ages estimates from cosmochronometry.

We present the results of a National Science Foundation (NSF) Project Scoping Workshop, the purpose of which was to assess the current status of calculations for the nuclear matrix elements governing neutrinoless double-beta decay and determine if more work on them is required. After reviewing important recent progress in the application of effective field theory, lattice quantum chromodynamics, and ab initio nuclear-structure theory to double-beta decay, we discuss the state of the art in nuclear-physics uncertainty quantification and then construct a road map for work in all these areas to fully complement the increasingly sensitive experiments in operation and under development. The road map contains specific projects in theoretical and computational physics as well as an uncertainty-quantification plan that employs Bayesian Model Mixing and an analysis of correlations between double-beta-decay rates and other observables. The goal of this program is a set of accurate and precise matrix elements, in all nuclei of interest to experimentalists, delivered together with carefully assessed uncertainties. Such calculations will allow crisp conclusions from the observation or non-observation of neutrinoless double-beta decay, no matter what new physics is at play., Comment: This Project Scoping Workshop report is focused on the US context for the theory of neutrinloess double beta decay. Its authors plan to produce a journal article that addresses similar issues, but is more inclusive as regards non-US efforts on this problem. We would be happy to receive further input that will help us refine our text before it is submitted to the journal

We use chiral interactions and several {\em ab initio} methods to compute the nuclear matrix elements (NMEs) for ground-state to ground-state double Gamow-Teller transitions in a range of isotopes, and explore the correlation of these NMEs with those for neutrinoless double beta decay produced by the exchange of a light Majorana neutrino. When all the NMEs of both isospin-conserving and isospin-changing transitions from the {\em ab initio} calculations are considered, the correlation is strong. For the experimentally relevant isospin-changing transitions by themselves, however, the correlation is weaker and may not be helpful for reducing the uncertainty in the NMEs for neutrinoless double-beta decay., Comment: 16 pages with 19 figures, Phys. Rev. C (in press)

We use the Lipkin-Meshkov-Glick (LMG) model and the valence-space nuclear shell model to examine the likely performance of variational quantum eigensolvers in nuclear-structure theory. The LMG model exhibits both a phase transition and spontaneous symmetry breaking at the mean-field level in one of the phases, features that characterize collective dynamics in medium-mass and heavy nuclei. We show that with appropriate modifications, the ADAPT-VQE algorithm, a particularly flexible and accurate variational approach, is not troubled by these complications. We treat up to 12 particles and show that the number of quantum operations needed to approach the ground-state energy scales linearly with the number of qubits. We find similar scaling when the algorithm is applied to the nuclear shell model with realistic interactions in the $sd$ and $pf$ shells. Although most of these simulations contain no noise, we use a noise model from real IBM hardware to show that for the LMG model with four particles, weak noise has no effect on the efficiency of the algorithm., Comment: 10 pages, 7 figures. Identical in content to published version. Now incudes analysis of noise

In light and medium-mass nuclei, two-body weak currents from chiral effective field theory account for a significant portion of the phenomenological quenching of Gamow-Teller transition matrix elements. Here we examine the systematic effects of two-body axial currents on Gamow-Teller strength and $\beta$-decay rates in heavy nuclei within energy-density functional theory. Using a Skyrme functional and the charge-changing finite amplitude method, we add the contributions of two-body currents to the usual one-body linear response in the Gamow-Teller channel, both exactly and though a density-matrix expansion. The two-body currents, as expected, usually quench both summed Gamow-Teller strength and decay rates, but by an amount that decreases as the neutron excess grows. In addition, they can enhance individual low-lying transitions, leading to decay rates that are quite different from those that an energy-independent quenching would produce, particularly in neutron-rich nuclei. We show that both these unexpected effects are related to changes in the total nucleon density as the number of neutrons increases., Comment: 13 pages, 8 figures

The temperature dependence of stellar electron-capture (EC) rates is investigated, with a focus on nuclei around $N=50$, just above $Z=28$, which play an important role during the collapse phase of core-collapse supernovae (CCSN). Two new microscopic calculations of stellar EC rates are obtained from a relativistic and a non-relativistic finite-temperature quasiparticle random-phase approximation approaches, for a conventional grid of temperatures and densities. In both approaches, EC rates due to Gamow-Teller transitions are included. In the relativistic calculation contributions from first-forbidden transitions are also included, and add strongly to the EC rates. The new EC rates are compared with large-scale shell model calculations for the specific case of $^{86}$Kr, providing insight into the finite-temperature effects on the EC rates. At relevant thermodynamic conditions for core-collapse, the discrepancies between the different calculations of this work are within about one order of magnitude. Numerical simulations of CCSN are performed with the spherically-symmetric GR1D simulation code to quantify the impact of such differences on the dynamics of the collapse. These simulations also include EC rates based on two parametrized approximations. A comparison of the neutrino luminosities and enclosed mass at core bounce shows that differences between simulations with different sets of EC rates are relatively small ($\approx 5\%$), suggesting that the EC rates used as inputs for these simulations have become well constrained.

The generator coordinate method begins with the variational construction of a set of non-orthogonal mean-field states that span a subspace of the full many-body Hilbert space. These states are then often projected onto states with good quantum numbers to restore symmetries, leading to a set with members that can be similar to one another, and it is sometimes possible to reduce this set without greatly affecting results. Here we propose a greedy algorithm that we call the energy-transition-orthogonality procedure (ENTROP) to select subsets of important states. As applied here, the approach selects on the basis of diagonal energy, orthogonality, and contribution to the matrix element that governs neutrinoless double-$\beta$ decay. We present both shell-model and preliminary ab initio calculations of this matrix element for the decay of $^{76}$Ge, with quadrupole deformation parameters and the isoscalar pairing strength as generator coordinates. ENTROP converges quickly, reducing significantly the number of basis states needed for an accurate calculation., Comment: 7 pages, 9 figures, authors added, some details clarified

We present 3D calculations for dielectric haloscopes such as the currently envisioned MADMAX experiment. For ideal systems with perfectly flat, parallel and isotropic dielectric disks of finite diameter, we find that a geometrical form factor reduces the emitted power by up to $30\,\%$ compared to earlier 1D calculations. We derive the emitted beam shape, which is important for antenna design. We show that realistic dark matter axion velocities of $10^{-3} c$ and inhomogeneities of the external magnetic field at the scale of $10\,\%$ have negligible impact on the sensitivity of MADMAX. We investigate design requirements for which the emitted power changes by less than $20\,\%$ for a benchmark boost factor with a bandwidth of $50\,{\rm MHz}$ at $22\,{\rm GHz}$, corresponding to an axion mass of $90\,\mu{\rm eV}$. We find that the maximum allowed disk tilt is $100\,\mu{\rm m}$ divided by the disk diameter, the required disk planarity is $20\,\mu{\rm m}$ (min-to-max) or better, and the maximum allowed surface roughness is $100\,\mu{\rm m}$ (min-to-max). We show how using tiled dielectric disks glued together from multiple smaller patches can affect the beam shape and antenna coupling., Comment: 25 pages, 12 figures; added detail, updated references, matches published version

We use coupled-cluster theory and nuclear interactions from chiral effective field theory to compute the nuclear matrix element for the neutrinoless double-beta decay of $^{48}$Ca. Benchmarks with the no-core shell model in several light nuclei inform us about the accuracy of our approach. For $^{48}$Ca we find a relatively small matrix element. We also compute the nuclear matrix element for the two-neutrino double-beta decay of $^{48}$Ca with a quenching factor deduced from two-body currents in recent ab-initio calculation of the Ikeda sum-rule in $^{48}$Ca [Gysbers et al., Nature Physics 15, 428-431 (2019)]., Comment: 14 pages, 13 figures; Version accepted for publication, Supplemental material also updated

We develop a fully self-consistent subtracted second random-phase approximation for charge-exchange processes with Skyrme energy-density functionals. As a first application, we study Gamow-Teller excitations in the doubly-magic nucleus $^{48}$Ca, the lightest double-$\beta$ emitter that could be used in an experiment, and in $^{78}$Ni, the single-beta-decay rate of which is known. The amount of Gamow-Teller strength below 20 or 30 MeV is considerably smaller than in other energy-density-functional calculations and agrees better with experiment in $^{48}$Ca, as does the beta-decay rate in $^{78}$Ni. These important results, obtained without \textit{ad hoc} quenching factors, are due to the presence of two-particle -- two-hole configurations. Their density progressively increases with excitation energy, leading to a long high-energy tail in the spectrum, a fact that may have implications for the computation of nuclear matrix elements for neutrinoless double-$\beta$ decay in the same framework.

We use the finite amplitude method (FAM), an efficient implementation of the quasiparticle random phase approximation, to compute beta-decay rates with Skyrme energy-density functionals for 3983 nuclei, essentially all the medium-mass and heavy isotopes on the neutron rich side of stability. We employ an extension of the FAM that treats odd-mass and odd-odd nuclear ground states in the equal filling approximation. Our rates are in reasonable agreement both with experimental data where available and with rates from other global calculations., Comment: 11 pages, 7 figures, supplemental material. Submitted to Phys. Rev. C

In this report we present the status of the MAgnetized Disk and Mirror Axion eXperiment (MADMAX), the first dielectric haloscope for the direct search of dark matter axions in the mass range of 40 to 400 $\mu$eV. MADMAX will consist of several parallel dielectric disks, which are placed in a strong magnetic field and with adjustable separations. This setting is expected to allow for an observable emission of axion induced electromagnetic waves at a frequency between 10 and 100 GHz corresponding to the axion mass. The present document orignated from a status report to the DESY PRC in 2019.

The Magnificent CE$\nu$NS Workshop (2018) was held November 2 & 3 of 2018 on the University of Chicago campus and brought together theorists, phenomenologists, and experimentalists working in numerous areas but sharing a common interest in the process of coherent elastic neutrino-nucleus scattering (CE$\nu$NS). This is a collection of abstract-like summaries of the talks given at the meeting, including links to the slides presented. This document and the slides from the meeting provide an overview of the field and a snapshot of the robust CE$\nu$NS-related efforts both planned and underway., Comment: The Magnificent CEvNS Workshop (2018), Nov 2-3, 2018; Chicago, IL, USA; 44 contributions

We compute nuclear matrix elements of neutrinoless double-beta decay mediated by light Majorana-neutrino exchange in the A = 6 system. The goal is to benchmark two many-body approaches, the No-Core Shell Model and the Multi-Reference In-Medium Similarity Renormalization Group. We use the SRG-evolved chiral N3LO-EM500 potential for the nuclear interaction, and make the approximation that isospin is conserved. We compare the results of the two approaches as a function of the cutoff on the many-body basis space. Although differences are seen in the predicted nuclear radii, the ground-state energies and neutrinoless double-beta decay matrix elements produced by the two approaches show significant agreement. We discuss the implications for calculations in heavier nuclei., Comment: 12 pages, 3 figures, 2 tables

Working with Hamiltonians from chiral effective field theory, we develop a novel framework for describing arbitrary deformed medium-mass nuclei by combining the in-medium similarity renormalization group with the generator coordinate method. The approach leverages the ability of the first method to capture dynamic correlations and the second to include collective correlations without violating symmetries. We use our scheme to compute the matrix element that governs the neutrinoless double beta decay of $^{48}$Ca to $^{48}$Ti, and find it to have the value $0.61$, near or below the predictions of most phenomenological methods. The result opens the door to ab initio calculations of the matrix elements for the decay of heavier nuclei such as $^{76}$Ge, $^{130}$Te, and $^{136}$Xe., Comment: 6 pages, 4 figures and 1 table. supplementary material included. version to be published

In the late stages of stellar core-collapse, prior to core bounce, electron captures on medium-heavy nuclei drive deleptonization and simulations require the use of accurate reaction rates. Nuclei with neutron number near $N=50$, just above atomic number $Z=28$, play an important role, but rates used in astrophysical simulations rely primarily on a relatively simple single-state approximation. In order to improve the accuracy of astrophysical simulations, experimental data are needed to test the electron-capture rates and to guide the development of better theoretical models. This work presents the results of the $^{86}$Kr($t$,$^{3}$He+$\gamma$) experiment at the NSCL, from which an upper limit for the Gamow-Teller strength up to an excitation energy in $^{86}$Br of 5 MeV is extracted. The derived upper limit for the electron-capture rate on $^{86}$Kr indicates that the rate estimated through the single-state approximation is too high and that rates based on Gamow-Teller strengths estimated in shell-model and QRPA calculations are more accurate. The QRPA calculations tested in this manner were used for estimating the electron capture rates for 78 isotopes near $N=50$ and above $Z=28$. The impact of using these new electron-capture rates in simulations of supernovae instead of the rates based on the single-state approximation is investigated, indicating a significant reduction in the deleptonization that affects multi-messenger signals, such as the emission of neutrinos and gravitational waves., Comment: 15 pages, 10 figures

The Gamow-Teller strength distribution from ${}^{88}$Sr was extracted from a $(t,{}^{3}\text{He}+\gamma)$ experiment at 115 MeV/$u$ to constrain estimates for the electron-capture rates on nuclei around $N=50$, between and including $^{78}$Ni and $^{88}$Sr, which are important for the late evolution of core-collapse supernovae. The observed strength below an excitation energy of 8 MeV was consistent with zero and below 10 MeV amounted to $0.1\pm0.05$. Except for a very-weak transition that could come from the 2.231-MeV $1^{+}$ state, no $\gamma$ lines that could be associated with the decay of known $1^{+}$ states were identified. The derived electron-capture rate from the measured strength distribution is more than an order of magnitude smaller than rates based on the single-state approximation presently used in astrophysical simulations for most nuclei near $N=50$. Rates based on shell-model and quasiparticle random-phase approximation calculations that account for Pauli blocking and core-polarization effects provide better estimates than the single-state approximation, although a relatively strong transition to the first $1^{+}$ state in $^{88}$Rb is not observed in the data. Pauli unblocking effects due to high stellar temperatures could partially counter the low electron-capture rates. The new data serves as a zero-temperature benchmark for constraining models used to estimate such effects.

The axion emerges in extensions of the Standard Model that explain the absence of CP violation in the strong interactions. Simultaneously, it can provide naturally the cold dark matter in our universe. Several searches for axions and axion-like particles (ALPs) have constrained the corresponding parameter space over the last decades but no unambiguous hints of their existence have been found. The axion mass range below 1 meV remains highly attractive and a well motivated region for dark matter axions. In this White Paper we present a description of a new experiment based on the concept of a dielectric haloscope for the direct search of dark matter axions in the mass range of 40 to 400 $\mu$eV. This MAgnetized Disk and Mirror Axion eXperiment (MADMAX) will consist of several parallel dielectric disks, which are placed in a strong magnetic field and with adjustable separations. This setting is expected to allow for an observable emission of axion induced electromagnetic waves at a frequency between 10 to 100 GHz corresponding to the axion mass.

We use a reference state based on symmetry-restored states from deformed mean-field or generator-coordinate-method (GCM) calculations in conjunction with the in-medium similarity-renormalization group (IMSRG) to compute spectra and matrix elements for neutrinoless double-beta ($0\nu\beta\beta$) decay. Because the decay involves ground states from two nuclei, we use evolved operators from the IMSRG in one nucleus in a subsequent GCM calculation in the other. We benchmark the resulting IMSRG+GCM method against complete shell-model diagonalization for both the energies of low-lying states in $^{48}$Ca and $^{48}$Ti and the $0\nu\beta\beta$ matrix element for the decay of $^{48}$Ca, all in a single valence shell. Our approach produces better spectra than either the IMSRG with a spherical-mean-field reference or GCM calculations with unevolved operators. For the $0\nu\beta\beta$ matrix element the improvement is slight, but we expect more significant effects in full ab-initio calculations., Comment: 12 pages, 5 figures, 1 tables

Background: The interacting boson model (IBM) has been used extensively to calculate the matrix elements governing neutrinoless double-beta decay. Studies within other models indicate that a good description of neutron-proton pairing is essential for accurate calculations of those matrix elements. The usual interacting boson model is based only on like-particle pairs, however, and the extent to which it captures neutron-proton pairing is not clear. Purpose: To determine whether neutron-proton pairing should be explicitly included as neutron-proton bosons in IBM calculations of neutrinoless double-beta decay matrix elements. Method: An isospin-invariant version of the nucleon-pair shell model is applied to carry out shell-model calculations in a large space and in a collective subspace, and to define effective operators in the latter. A democratic mapping is then used to define corresponding boson operators for the IBM, with and without an isoscalar neutron-proton pair boson. Results: IBM calculations with and without the isoscalar boson are carried out for nuclei near the beginning of the $pf$ shell, with a realistic shell-model Hamiltonian and neutrinoless double-beta-decay operator as the starting point. Energy spectra and double-beta matrix elements are compared to those obtained in the underlying shell model. Conclusions: The isoscalar boson does not improve energy spectra but does improve double-beta matrix elements. To be useful at the level of precision we need, the mapping procedure must be further developed to better determine the dependence of the boson Hamiltonian and decay operator on particle number and isospin. But the benefits provided by the isoscalar boson suggest that through an appropriate combination of mappings and fitting, it would make IBM matrix elements more accurate for the heavier nuclei used in experiments., Comment: 17 pages, 8 figures

We use the generator-coordinate method with realistic shell-model interactions to closely approximate full shell-model calculations of the matrix elements for the neutrinoless double-beta decay of $^{48}$Ca, $^{76}$Ge, and $^{82}$Se. We work in one major shell for the first isotope, in the $f_{5/2}pg_{9/2}$ space for the second and third, and finally in two major shells for all three. Our coordinates include not only the usual axial deformation parameter $\beta$, but also the triaxiality angle $\gamma$ and neutron-proton pairing amplitudes. In the smaller model spaces our matrix elements agree well with those of full shell-model diagonalization, suggesting that our Hamiltonian-based GCM captures most of the important valence-space correlations. In two major shells, where exact diagonalization is not currently possible, our matrix elements are only slightly different from those in a single shell., Comment: 8 pages, 7 figures

After identifying the nuclei in the regions near A =80 and A = 160 for which beta-decay rates have the greatest effect on weak and main r-process abundance patterns, we apply the finite-amplitude method (FAM) with Skyrme energy-density functionals (EDFs) to calculate beta-decay half-lives of those nuclei in the quasiparticle random-phase approximation (QRPA). We use the equal filling approximation to extend our implementation of the charge-changing FAM, which incorporates pairing correlations and allows axially symmetric deformation, to odd-A and odd-odd nuclei. Within this framework we find differences of up to a factor of seven between our calculated beta-decay half-lives and those of previous efforts. Repeated calculations with nuclei near A = 160 and multiple EDFs show a spread of two to four in beta-decay half-lives, with differences in calculated Q values playing an important role. We investigate the implications of these results for r-process simulations., Comment: 18 pages, 16 figures

We present a generator-coordinate calculation, based on a relativistic energy-density functional, of the low-lying spectra in the isotopes $^{150}$Nd and $^{150}$Sm and of the nuclear matrix element that governs the neutrinoless double-beta decay of the first isotope to the second. We carefully examine the impact of octupole correlations on both nuclear structure and the double-beta decay matrix element. Octupole correlations turn out to reduce quadrupole collectivity in both nuclei. Shape fluctuations, however, dilute the effects of octupole deformation on the double-beta decay matrix element, so that the overall octupole-induced quenching is only about 7\%., Comment: 7 pages, 5 figures, to be published in Phys. Rev. C as a regular article

We investigate the extent to which theories of collective motion can capture the physics that determines the nuclear matrix elements governing neutrinoless double-beta decay. To that end we calculate the matrix elements for a series of isotopes in the full $pf$ shell, omitting no spin-orbit partners. With the inclusion of isoscalar pairing, a separable collective Hamiltonian that is derived from the shell model effective interaction reproduces the full shell-model matrix elements with good accuracy. A version of the generator coordinate method that includes the isoscalar pairing amplitude as a coordinate also reproduces the shell model results well, an encouraging result for theories of collective motion, which can include more single-particle orbitals than the shell model. We briefly examine heavier nuclei relevant for experimental double-beta decay searches, in which shell-model calculations with all spin-orbit partners are not feasible; our estimates suggest that isoscalar pairing also plays a significant role in these nuclei, though one we are less able to quantify precisely., Comment: 9 pages, 7 figures, published version

We use the finite amplitude method for computing charge-changing Skyrme-QRPA transition strengths in axially-deformed nuclei together with a modern Skyrme energy-density functional to fit several previously unconstrained parameters in the charge-changing time-odd part of the functional. With the modified functional we then calculate rates of beta-minus decay for all medium-mass and heavy even-even nuclei between the valley of stability and the neutron drip line. We fit the Skyrme parameters to a limited set of beta-decay rates, a set of Gamow-Teller resonance energies, and a set of spin-dipole resonance energies, in both spherical and deformed nuclei. Comparison to available experimental beta-decay rates shows agreement at roughly the same level as in other global QRPA calculations. We estimate the uncertainty in our rates all the way to the neutron drip line through a construction that extrapolates the errors of known beta-decay rates in nuclei with intermediate Q values to less stable isotopes with higher Q values., Comment: 11 pages, 7 figures

We make use of a subtraction procedure, introduced to overcome double--counting problems in beyond--mean--field theories, in the second random--phase--approximation (SRPA) for the first time. This procedure guarantees the stability of SRPA (so that all excitation energies are real). We show that the method fits perfectly into nuclear density--functional theory. We illustrate applications to the monopole and quadrupole response and to low--lying $0^+$ and $2^+$ states in the nucleus $^{16}$O. We show that the subtraction procedure leads to: (i) results that are weakly cutoff dependent; (ii) a considerable reduction of the SRPA downwards shift with respect to the random--phase approximation (RPA) spectra (systematically found in all previous applications). This implementation of the SRPA model will allow a reliable analysis of the effects of 2 particle--2 hole configurations ($2p2h$) on the excitation spectra of medium--mass and heavy nuclei., Comment: 1 tex, 16 figures

We describe and apply a version of the finite amplitude method for obtaining the charge-changing nuclear response in the quasiparticle random phase approximation. The method is suitable for calculating strength functions and beta-decay rates, both allowed and forbidden, in axially-deformed open-shell nuclei. We demonstrate the speed and versatility of the code through a preliminary examination of the effects of tensor terms in Skyrme functionals on beta decay in a set of spherical and deformed open-shell nuclei. Like the isoscalar pairing interaction, the tensor terms systematically increase allowed beta-decay rates. This finding generalizes previous work in semimagic nuclei and points to the need for a comprehensive study of time-odd terms in nuclear density functionals., Comment: 11 pages, 8 figures, submitted to Physical Review C

We test the effects of an approximate treatment of two-body contributions to the axial-vector current on the QRPA matrix elements for neutrinoless double-beta decay in a range of isotopes. The form and strength of the two-body terms come from chiral effective-field theory. The two-body currents typically reduce the matrix elements by about 20%, not as much as in shell-model calculations. One reason for the difference is that standard practice in the QRPA is to adjust the strength of the isoscalar pairing interaction to reproduce two-neutrino double-beta decay lifetimes. Another may be the larger QRPA single-particle space. Whatever the reasons, the effects on neutrinoless decay are significantly less than those on two-neutrino decay, both in the shell model and the QRPA., Comment: 6 pages, 2 figures