Search

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

Quantum state preparation by adiabatic evolution is currently rendered ineffective by the long implementation times of the underlying quantum circuits, comparable to the decoherence time of present and near-term quantum devices. These implementation times can be significantly reduced by realizing the evolution with a minimal number of customized gates. Employing a realistic model of a two-qubit processor, we carried out classical device-level simulations of the adiabatic evolution of a two-spin system implemented with customized two-qubit gates. These device-level simulations were compared with (experimental) ones solving the same problem on IBMQ systems. When used to emulate the IBMQ quantum circuit, our device-level simulations reached state fidelities ranging from 65% to 85%, similar to the actual performance of a diverse set of IBMQ devices.When we reduced the implementation time by using a minimal number of customized gates, however, the loss of fidelity was reduced by at least a factor of four, allowing us to accurately extract the energy of the target state. This improvement is enough to render adiabatic evolution useful for quantum state preparation for small systems or as a preconditioner for other state preparation methods., Comment: 12 pages, 10 figures

Combining laser spectroscopy in a Versatile Arc Discharge and Laser Ion Source, with Penning-trap mass spectrometry at the CERN-ISOLDE facility, this work reports on mean-square charge radii of neutron-rich mercury isotopes across the $N = 126$ shell closure, the electromagnetic moments of $^{207}$Hg and more precise mass values of $^{206-208}$Hg. The odd-even staggering (OES) of the mean square charge radii and the kink at $N = 126$ are analyzed within the framework of covariant density functional theory (CDFT), with comparisons between different functionals to investigate the dependence of the results on the underlying single-particle structure. The observed features are defined predominantly in the particle-hole channel in CDFT, since both are present in the calculations without pairing. However, the magnitude of the kink is still affected by the occupation of the $1i_{11/2}$ and $2g_{9/2}$ orbitals with a dependence on the relative energies as well as pairing., Comment: 20 pages, 12 figures, Phys. Rev. C in press

The ECHo experiment has been designed for the determination of the effective electron neutrino mass by means of the analysis of the end-point region of the Ho-163 electron capture spectrum. Metallic magnetic calorimeters enclosing Ho-163 are used for the high energy resolution calorimetric measurement of the Ho-163 spectrum. For the first phase of the experiment, ECHo-1k, a 72-pixel MMC array has been developed. The single-pixel design has been optimised to reach 100% stopping power for the radiation emitted in the Ho-163 electron capture process (besides the electron neutrino) and an energy resolution < 10 eV FWHM. We describe the design of the ECHo-1k detector chip, the fabrication steps and the characterisation at room temperature, at 4 K and at the final operation temperatures. In particular, a detailed analysis of the results from these tests allowed to define a quality check protocol based on parameters measurable at room temperature. We discuss the performance achieved with the two ECHo-1k detector chips - the first one with Ho-163 implanted in gold and the second one with Ho-163 implanted in silver - which have been used for the high statistics measurement of the ECHo-1k experiment. An average activity per pixel of 0.81 Bq and 0.71 Bq and an average energy resolution of 6.07 eV FWHM and 5.55 eV FWHM have been achieved with these two detectors, fulfilling the requirements for the first phase of the ECHo experiment.

Isotope shifts of $^{223-226,228}$Ra$^{19}$F were measured for different vibrational levels in the electronic transition $A^{2}{}{\Pi}_{1/2}\leftarrow X^{2}{}{\Sigma}^{+}$. The observed isotope shifts demonstrate the particularly high sensitivity of radium monofluoride to nuclear size effects, offering a stringent test of models describing the electronic density within the radium nucleus. Ab initio quantum chemical calculations are in excellent agreement with experimental observations. These results highlight some of the unique opportunities that short-lived molecules could offer in nuclear structure and in fundamental symmetry studies., Comment: The paper was accepted for publication as a Letter in Physical Review Letters

The mean-square charge radii of $^{207,208}$Hg ($Z=80, N=127,128$) have been studied for the first time and those of $^{202,203,206}$Hg ($N=122,123,126$) remeasured by the application of in-source resonance-ionization laser spectroscopy at ISOLDE (CERN). The characteristic \textit{kink} in the charge radii at the $N=126$ neutron shell closure has been revealed, providing the first information on its behavior below the $Z=82$ proton shell closure. A theoretical analysis has been performed within relativistic Hartree-Bogoliubov and non-relativistic Hartree-Fock-Bogoliubov approaches, considering both the new mercury results and existing lead data. Contrary to previous interpretations, it is demonstrated that both the kink at $N=126$ and the odd-even staggering (OES) in its vicinity can be described predominately at the mean-field level, and that pairing does not need to play a crucial role in their origin. A new OES mechanism is suggested, related to the staggering in the occupation of the different neutron orbitals in odd- and even-$A$ nuclei, facilitated by particle-vibration coupling for odd-$A$ nuclei., Comment: 7 pages, 5 figures, Physical Review C, in press

The mesoscopic nature of the atomic nucleus gives rise to a wide array of macroscopic and microscopic phenomena. The size of the nucleus is a window into this duality: while the charge radii globally scale as $A^{1/3}$, their evolution across isotopic chains reveals unanticipated structural phenomena [1-3]. The most ubiquitous of these is perhaps the Odd-Even Staggering (OES) [4]: isotopes with an odd number of neutrons are usually smaller in size than the trend of their even-neutron neighbours suggests. This OES effect varies with the number of protons and neutrons and poses a significant challenge for nuclear theory [5-7]. Here, we examine this problem with new measurements of the charge radii of short-lived copper isotopes up to the very exotic $^{78}$Cu $(Z=29, N=49)$, produced at only 20 ions/s, using the highly-sensitive Collinear Resonance Ionisation Spectroscopy (CRIS) method at ISOLDE-CERN. Due to the presence of a single proton outside of the closed Z=28 shell, these measurements provide crucial insights into the single-particle proton structure and how this affects the charge radii. We observe an unexpected reduction in the OES for isotopes approaching the $N=50$ shell gap. To describe the data, we applied models based on nuclear Density Functional Theory [2,8] (DFT) and ab-initio Valence-Space In-Medium Similarity Renormalization Group (VS-IMSRG) theory [9,10]. Through these comparisons, we demonstrate a relation between the global behavior of charge radii and the saturation density of nuclear matter, and show that the local charge radii variations, which reflect the many-body polarization effects due to the odd neutron, naturally emerge from the VS-IMSRG calculations.

With increasing demand for accurate calculation of isotope shifts of atomic systems for fundamental and nuclear structure research, an analytic energy derivative approach is presented in the relativistic coupled-cluster theory framework to determine the atomic field shift and mass shift factors. This approach allows the determination of expectation values of atomic operators, overcoming fundamental problems that are present in existing atomic physics methods, i.e. it satisfies the Hellmann-Feynman theorem, does not involve any non-terminating series, and is free from choice of any perturbative parameter. As a proof of concept, the developed analytic response relativistic coupled-cluster theory has been applied to determine mass shift and field shift factors for different atomic states of indium. High-precision isotope-shift measurements of $^{104-127}$In were performed in the 246.8-nm (5p $^2$P$_{3/2}$ $\rightarrow$ 9s $^2$S$_{1/2}$) and 246.0-nm (5p $^2$P$_{1/2}$ $\rightarrow$ 8s $^2$S$_{1/2}$) transitions to test our theoretical results. An excellent agreement between the theoretical and measured values is found, which is known to be challenging in multi-electron atoms. The calculated atomic factors allowed an accurate determination of the nuclear charge radii of the ground and isomeric states of the $^{104-127}$In isotopes, providing an isotone-independent comparison of the absolute charge radii., Comment: Submitted

The study of molecular systems provides exceptional opportunities for the exploration of the fundamental laws of nature and for the search for physics beyond the Standard Model of particle physics. Measurements of molecules composed of naturally occurring nuclei have provided the most stringent upper bounds to the electron electric dipole moment to date, and offer a route to investigate the violation of fundamental symmetries with unprecedented sensitivity. Radioactive molecules - where one or more of their atoms possesses a radioactive nucleus - can contain heavy and deformed nuclei, offering superior sensitivity for EDM measurements as well as for other symmetry-violating effects. Radium monofluoride, RaF, is of particular interest as it is predicted to have an appropriate electronic structure for direct laser cooling. Furthermore, some Ra isotopes are known to be octupole deformed, thereby resulting in a large enhancement of their symmetry-violating nuclear moments. Until now,however, no experimental measurements of RaF have been performed, and their study is impeded by major experimental challenges, as no stable isotopes of radium exist. Here, we present a novel experimental approach to study short-lived radioactive molecules using the highly sensitive collinear resonance ionisation method. With this technique we have measured, for the first time, the energetically low-lying electronic states for each of the isotopically pure RaF molecules at the ISOLDE-CERN. Our results provide strong evidence of the existence of a suitable laser-cooling scheme for these molecules and constitute a pivotal step towards high-precision studies in these systems. Our findings open up new opportunities in the synthesis, manipulation and study of short-lived radioactive molecules, which will have a direct impact in many-body physics, astrophysics, nuclear structure, and fundamental physics research., Comment: Submitted

Neutron-deficient $^{177-185}$Hg isotopes were studied using in-source laser resonance-ionization spectroscopy at the CERN-ISOLDE radioactive ion-beam facility, in an experiment combining different detection methods tailored to the studied isotopes. These include either alpha-decay tagging or Multi-reflection Time-of-Flight gating to identify the isotopes of interest. The endpoint of the odd-even nuclear shape staggering in mercury was observed directly by measuring for the first time the isotope shifts and hyperfine structures of $^{177-180}$Hg. Changes in the mean-square charge radii for all mentioned isotopes, magnetic dipole and electric quadrupole moments of the odd-A isotopes and arguments in favor of $I = 7/2$ spin assignment for $^{177,179}$Hg were deduced. Experimental results are compared with Density Functional Theory (DFT) and Monte-Carlo Shell Model (MCSM) calculations. DFT calculations with several Skyrme parameterizations predict a large jump in the charge radius around the neutron $N = 104$ mid shell, with an odd-even staggering pattern related to the coexistence of nearly-degenerate oblate and prolate minima. This near-degeneracy is highly sensitive to many aspects of the effective interaction, a fact that renders perfect agreement with experiment out of reach for current functionals. Despite this inherent diffculty, the SLy5s1 and a modified UNEDF1^{SO} parameterization predict a qualitatively correct staggering that is off by two neutron numbers. MCSM calculations of states with the experimental spins and parities show good agreement for both electromagnetic moments and the observed charge radii. A clear mechanism for the origin of shape staggering within this context is identified: a substantial change in occupancy of the proton $\pi h_{9/2}$ and neutron $\nu i_{13/2}$ orbitals., Comment: 19 pages, 18 figures

$\beta$-decay, a process that changes a neutron into a proton (and vice versa), is the dominant decay mode of atomic nuclei. This decay offers a unique window to physics beyond the standard model, and is at the heart of microphysical processes in stellar explosions and the synthesis of the elements in the Universe. For 50 years, a central puzzle has been that observed $\beta$-decay rates are systematically smaller than theoretical predictions. This was attributed to an apparent quenching of the fundamental coupling constant $g_A \simeq $ 1.27 in the nucleus by a factor of about 0.75 compared to the $\beta$-decay of a free neutron. The origin of this quenching is controversial and has so far eluded a first-principles theoretical understanding. Here we address this puzzle and show that this quenching arises to a large extent from the coupling of the weak force to two nucleons as well as from strong correlations in the nucleus. We present state-of-the-art computations of $\beta$-decays from light to heavy nuclei. Our results are consistent with experimental data, including the pioneering measurement for $^{100}$Sn. These theoretical advances are enabled by systematic effective field theories of the strong and weak interactions combined with powerful quantum many-body techniques. This work paves the way for systematic theoretical predictions for fundamental physics problems. These include the synthesis of heavy elements in neutron star mergers and the search for neutrino-less double-$\beta$-decay, where an analogous quenching puzzle is a major source of uncertainty in extracting the neutrino mass scale., Comment: 20 pages, 18 figures

Nuclear spins and precise values of the magnetic dipole and electric quadrupole moments of the ground-states of neutron-rich $^{76-78}$Cu isotopes were measured using the Collinear Resonance Ionization Spectroscopy (CRIS) experiment at ISOLDE, CERN. The nuclear moments of the less exotic $^{73,75}$Cu isotopes were re-measured with similar precision, yielding values that are consistent with earlier measurements. The moments of the odd-odd isotopes, and $^{78}_{29}$Cu ($N=49$) in particular, are used to investigate excitations of the assumed doubly-magic $^{78}$Ni core through comparisons with large-scale shell-model calculations. Despite the narrowing of the $Z=28$ shell gap between $N\sim45$ and $N=50$, the magicity of $Z=28$ and $N=50$ is restored towards $^{78}$Ni. This is due to weakened dynamical correlations, as clearly probed by the present moment measurements.

Long-lived isotopes of plutonium were studied using two complementary techniques, high-resolution resonance ionisation spectroscopy (HR-RIS) and collinear laser spectroscopy (CLS). Isotope shifts have been measured on the $5f^67s^2\ ^7F_0 \rightarrow 5f^56d^27s\ (J=1)$ and $5f^67s^2\ ^7F_1 \rightarrow 5f^67s7p\ (J=2)$ atomic transitions using the HR-RIS method and the hyperfine factors have been extracted for the odd mass nuclei $^{239,241}$Pu. Collinear laser spectroscopy was performed on the $5f^67s\ ^8F_{1/2} \rightarrow J=1/2\; (27523.61\text{cm}^{-1})$ ionic transition with the hyperfine $A$ factors measured for $^{239}$Pu. Changes in mean-squared charge radii have been extracted and show a good agreement with previous non-optical methods, with an uncertainty improvement by approximately one order of magnitude. Plutonium represents the heaviest element studied to date using collinear laser spectroscopy., Comment: accepted by Phys. Rev. A (2017-03-02)

Despite being a complex many-body system, the atomic nucleus exhibits simple structures for certain "magic" numbers of protons and neutrons. The calcium chain in particular is both unique and puzzling: evidence of doubly-magic features are known in 40,48Ca, and recently suggested in two radioactive isotopes, 52,54Ca. Although many properties of experimentally known Ca isotopes have been successfully described by nuclear theory, it is still a challenge to predict their charge radii evolution. Here we present the first measurements of the charge radii of 49,51,52Ca, obtained from laser spectroscopy experiments at ISOLDE, CERN. The experimental results are complemented by state-of-the-art theoretical calculations. The large and unexpected increase of the size of the neutron-rich calcium isotopes beyond N = 28 challenges the doubly-magic nature of 52Ca and opens new intriguing questions on the evolution of nuclear sizes away from stability, which are of importance for our understanding of neutron-rich atomic nuclei., Comment: Submitted version. See original publication (doi:10.1038/nphys3645) for final version

We present a comprehensive review of keV-scale sterile neutrino Dark Matter, collecting views and insights from all disciplines involved - cosmology, astrophysics, nuclear, and particle physics - in each case viewed from both theoretical and experimental/observational perspectives. After reviewing the role of active neutrinos in particle physics, astrophysics, and cosmology, we focus on sterile neutrinos in the context of the Dark Matter puzzle. Here, we first review the physics motivation for sterile neutrino Dark Matter, based on challenges and tensions in purely cold Dark Matter scenarios. We then round out the discussion by critically summarizing all known constraints on sterile neutrino Dark Matter arising from astrophysical observations, laboratory experiments, and theoretical considerations. In this context, we provide a balanced discourse on the possibly positive signal from X-ray observations. Another focus of the paper concerns the construction of particle physics models, aiming to explain how sterile neutrinos of keV-scale masses could arise in concrete settings beyond the Standard Model of elementary particle physics. The paper ends with an extensive review of current and future astrophysical and laboratory searches, highlighting new ideas and their experimental challenges, as well as future perspectives for the discovery of sterile neutrinos., Comment: v2: 257 pages, 57 figures, content matches published version [JCAP01(2017)025]; over 100 authors from several different communities

We develop interactions from chiral effective field theory (EFT) that are tailored to the harmonic oscillator basis. As a consequence, ultraviolet convergence with respect to the model space is implemented by construction and infrared convergence can be achieved by enlarging the model space for the kinetic energy. In oscillator EFT, matrix elements of EFTs formulated for continuous momenta are evaluated at the discrete momenta that stem from the diagonalization of the kinetic energy in the finite oscillator space. By fitting to realistic phase shifts and deuteron data we construct an effective interaction from chiral EFT at next-to-leading order. Many-body coupled-cluster calculations of nuclei up to 132Sn exhibit a fast convergence of ground-state energies and radii in feasible model spaces., Comment: 16 pages, 6 figures; published version

What is the size of the atomic nucleus? This deceivably simple question is difficult to answer. While the electric charge distributions in atomic nuclei were measured accurately already half a century ago, our knowledge of the distribution of neutrons is still deficient. In addition to constraining the size of atomic nuclei, the neutron distribution also impacts the number of nuclei that can exist and the size of neutron stars. We present an ab initio calculation of the neutron distribution of the neutron-rich nucleus $^{48}$Ca. We show that the neutron skin (difference between radii of neutron and proton distributions) is significantly smaller than previously thought. We also make predictions for the electric dipole polarizability and the weak form factor; both quantities are currently targeted by precision measurements. Based on ab initio results for $^{48}$Ca, we provide a constraint on the size of a neutron star.

Chiral effective field theory (chiEFT) provides a systematic approach to describe low-energy nuclear forces. Moreover, chiEFT is able to provide well-founded estimates of statistical and systematic uncertainties -- although this unique advantage has not yet been fully exploited. We fill this gap by performing an optimization and statistical analysis of all the low-energy constants (LECs) up to next-to-next-to-leading order. Our optimization protocol corresponds to a simultaneous fit to scattering and bound-state observables in the pion-nucleon, nucleon-nucleon, and few-nucleon sectors, thereby utilizing the full model capabilities of chiEFT. We study the effect on other observables by demonstrating error-propagation methods that can easily be adopted by future works. We employ mathematical optimization and implement automatic differentiation to attain efficient and machine-precise first- and second-order derivatives of the objective function with respect to the LECs. We use power-counting arguments to estimate the systematic uncertainty that is inherent to chiEFT and we construct chiral interactions at different orders with quantified uncertainties. Statistical error propagation is compared with Monte Carlo sampling showing that statistical errors are in general small compared to systematic ones. In conclusion, we find that a simultaneous fit to different sets of data is critical to (i) identify the optimal set of LECs, (ii) capture all relevant correlations, (iii) reduce the statistical uncertainty, and (iv) attain order-by-order convergence in chiEFT. Furthermore, certain systematic uncertainties in the few-nucleon sector are shown to get substantially magnified in the many-body sector; in particlar when varying the cutoff in the chiral potentials. The methodology and results presented in this Paper open a new frontier for uncertainty quantification in ab initio nuclear theory., Comment: Accepted for publication in Physical Review X. This version includes Supplemental Material

We precisely determine the infrared (IR) length scale of the no-core shell model (NCSM). In the NCSM, the $A$-body Hilbert space is truncated by the total energy, and the IR length can be determined by equating the intrinsic kinetic energy of $A$ nucleons in the NCSM space to that of $A$ nucleons in a $3(A-1)$-dimensional hyper-radial well with a Dirichlet boundary condition for the hyper radius. We demonstrate that this procedure indeed yields a very precise IR length by performing large-scale NCSM calculations for $^{6}$Li. We apply our result and perform accurate IR extrapolations for bound states of $^{4}$He, $^{6}$He, $^{6}$Li, $^{7}$Li. We also attempt to extrapolate NCSM results for $^{10}$B and $^{16}$O with bare interactions from chiral effective field theory over tens of MeV., Comment: 8 pages, 4 figures. Updated to address issues fixed during publication. We have prepared a python (v2.7) script that will compute the IR scales (both full scale and intrinsic scale). This script depends only on numpy and scipy and can be found at http://github.com/KyleWendt/ncsm_ho_ir_scale

With the goal of developing predictive ab-initio capability for light and medium-mass nuclei, two-nucleon and three-nucleon forces from chiral effective field theory are optimized simultaneously to low-energy nucleon-nucleon scattering data, as well as binding energies and radii of few-nucleon systems and selected isotopes of carbon and oxygen. Coupled-cluster calculations based on this interaction, named NNLO$_{\rm sat}$, yield accurate binding energies and radii of nuclei up to $^{40}$Ca, and are consistent with the empirical saturation point of symmetric nuclear matter. In addition, the low-lying collective $J^\pi=3^-$ states in $^{16}$O and $^{40}$Ca are described accurately, while spectra for selected $p$- and $sd$-shell nuclei are in reasonable agreement with experiment.