78Ni revealed as a doubly magic stronghold against nuclear deformation

High-Precision Mass Measurements of Neutron Deficient Silver Isotopes Probe the Robustness of the N=50 Shell Closure

High-precision mass measurements of exotic ^{95-97}Ag isotopes close to the N=Z line have been conducted with the JYFLTRAP double Penning trap mass spectrometer, with the silver ions produced using the recently commissioned inductively heated hot cavity catcher laser ion source at the Ion Guide Isotope Separator On-Line facility. The atomic mass of ^{95}Ag was directly determined for the first time. In addition, the atomic masses of β-decaying 2^{+} and 8^{+} states in ^{96}Ag have been identified and measured for the first time, and the precision of the ^{97}Ag mass has been improved. The newly measured masses, with a precision of ≈1  keV/c^{2}, have been used to investigate the N=50 neutron shell closure, confirming it to be robust. Empirical shell-gap and pairing energies determined with the new ground-state mass data are compared with the state-of-the-art ab initio calculations with various chiral effective field theory Hamiltonians. The precise determination of the excitation energy of the ^{96m}Ag isomer in particular serves as a benchmark for ab initio predictions of nuclear properties beyond the ground state, specifically for odd-odd nuclei situated in proximity to the proton dripline below ^{100}Sn. In addition, density functional theory calculations and configuration-interaction shell-model calculations are compared with the experimental results. All theoretical approaches face challenges to reproduce the trend of nuclear ground-state properties in the silver isotopic chain across the N=50 neutron shell and toward the proton dripline.

First Exploration of Monopole-Driven Shell Evolution above the N=126 Shell Closure: New Millisecond Isomers in ^{213}Tl and ^{215}Tl

Isomer spectroscopy of heavy neutron-rich nuclei beyond the N=126 closed shell has been performed for the first time at the Radioactive Isotope Beam Factory of the RIKEN Nishina Center. New millisecond isomers have been identified at low excitation energies, 985.3(19) keV in ^{213}Tl and 874(5) keV in ^{215}Tl. The measured half-lives of 1.34(5) ms in ^{213}Tl and 3.0(3) ms in ^{215}Tl suggest spins and parities 11/2^{-} with the single proton-hole configuration πh_{11/2} as leading component. They are populated via E1 transitions by the decay of higher-lying isomeric states with proposed spin and parity 17/2^{+}, interpreted as arising from a single πs_{1/2} proton hole coupled to the 8^{+} seniority isomer in the ^{A+1}Pb cores. The lowering of the 11/2^{-} states is ascribed to an increase of the πh_{11/2} proton effective single-particle energy as the second νg_{9/2} orbital is filled by neutrons, owing to a significant reduction of the proton-neutron monopole interaction between the πh_{11/2} and νg_{9/2} orbitals. The new ms isomers provide the first experimental observation of shell evolution in the almost unexplored N>126 nuclear region below doubly magic ^{208}Pb.

Nuclear structure opportunities with GeV radioactive beams at FAIR

The Facility for Antiproton and Ion Research (FAIR) is in its final construction stage next to the campus of the Gesellschaft für Schwerionenforschung Helmholtzzentrum for heavy-ion research in Darmstadt, Germany. Once it starts its operation, it will be the main nuclear physics research facility in many basic sciences and their applications in Europe for the coming decades. Owing to the ability of the new fragment separator, Super-FRagment Separator, to produce high-intensity radioactive ion beams in the energy range up to about 2 GeV/nucleon, these can be used in various nuclear reactions. This opens a unique opportunity for various nuclear structure studies across a range of fields and scales: from low-energy physics via the investigation of multi-neutron systems and halos to high-density nuclear matter and the equation of state, following heavy-ion collisions, fission and study of short-range correlations in nuclei and hypernuclei. The newly developed reactions with relativistic radioactive beams (R3B) set up at FAIR would be the most suitable and versatile for such studies. An overview of highlighted physics cases foreseen at R3B is given, along with possible future opportunities, at FAIR. This article is part of the theme issue 'The liminal position of Nuclear Physics: from hadrons to neutron stars'.

Demonstration of nuclear gamma-ray polarimetry based on a multi-layer CdTe Compton camera

To detect and track structural changes in atomic nuclei, the systematic study of nuclear levels with firm spin-parity assignments is important. While linear polarization measurements have been applied to determine the electromagnetic character of gamma-ray transitions, the applicable range is strongly limited due to the low efficiency of the detection system. The multi-layer Cadmium-Telluride (CdTe) Compton camera can be a state-of-the-art gamma-ray polarimeter for nuclear spectroscopy with the high position sensitivity and the detection efficiency. We demonstrated the capability to operate this detector as a reliable gamma-ray polarimeter by using polarized 847-keV gamma rays produced by the [Formula: see text]([Formula: see text]) reaction. By combining the experimental data and simulated calculations, the modulation curve for the gamma ray was successfully obtained. A remarkably high polarization sensitivity was achieved, compatible with a reasonable detection efficiency. Based on the obtained results, a possible future gamma-ray polarimetery is discussed.

Further Evidence for Shape Coexistence in ^{79}Zn^{m} near Doubly Magic ^{78}Ni

Isomers close to doubly magic _{28}^{78}Ni_{50} provide essential information on the shell evolution and shape coexistence near the Z=28 and N=50 double shell closure. We report the excitation energy measurement of the 1/2^{+} isomer in _{30}^{79}Zn_{49} through independent high-precision mass measurements with the JYFLTRAP double Penning trap and with the ISOLTRAP multi-reflection time-of-flight mass spectrometer. We unambiguously place the 1/2^{+} isomer at 942(10) keV, slightly below the 5/2^{+} state at 983(3) keV. With the use of state-of-the-art shell-model diagonalizations, complemented with discrete nonorthogonal shell-model calculations which are used here for the first time to interpret shape coexistence, we find low-lying deformed intruder states, similar to other N=49 isotones. The 1/2^{+} isomer is interpreted as the bandhead of a low-lying deformed structure akin to a predicted low-lying deformed band in ^{80}Zn, and points to shape coexistence in ^{79,80}Zn similar to the one observed in ^{78}Ni. The results make a strong case for confirming the claim of shape coexistence in this key region of the nuclear chart.

First observation of 28O

Subjecting a physical system to extreme conditions is one of the means often used to obtain a better understanding and deeper insight into its organization and structure. In the case of the atomic nucleus, one such approach is to investigate isotopes that have very different neutron-to-proton (N/Z) ratios than in stable nuclei. Light, neutron-rich isotopes exhibit the most asymmetric N/Z ratios and those lying beyond the limits of binding, which undergo spontaneous neutron emission and exist only as very short-lived resonances (about 10-21 s), provide the most stringent tests of modern nuclear-structure theories. Here we report on the first observation of 28O and 27O through their decay into 24O and four and three neutrons, respectively. The 28O nucleus is of particular interest as, with the Z = 8 and N = 20 magic numbers1,2, it is expected in the standard shell-model picture of nuclear structure to be one of a relatively small number of so-called 'doubly magic' nuclei. Both 27O and 28O were found to exist as narrow, low-lying resonances and their decay energies are compared here to the results of sophisticated theoretical modelling, including a large-scale shell-model calculation and a newly developed statistical approach. In both cases, the underlying nuclear interactions were derived from effective field theories of quantum chromodynamics. Finally, it is shown that the cross-section for the production of 28O from a 29F beam is consistent with it not exhibiting a closed N = 20 shell structure.

Nuclear shell-model simulation in digital quantum computers

The nuclear shell model is one of the prime many-body methods to study the structure of atomic nuclei, but it is hampered by an exponential scaling on the basis size as the number of particles increases. We present a shell-model quantum circuit design strategy to find nuclear ground states by exploiting an adaptive variational quantum eigensolver algorithm. Our circuit implementation is in excellent agreement with classical shell-model simulations for a dozen of light and medium-mass nuclei, including neon and calcium isotopes. We quantify the circuit depth, width and number of gates to encode realistic shell-model wavefunctions. Our strategy also addresses explicitly energy measurements and the required number of circuits to perform them. Our simulated circuits approach the benchmark results exponentially with a polynomial scaling in quantum resources for each nucleus. This work paves the way for quantum computing shell-model studies across the nuclear chart and our quantum resource quantification may be used in configuration-interaction calculations of other fermionic systems.

First Evidence of Axial Shape Asymmetry and Configuration Coexistence in ^{74}Zn: Suggestion for a Northern Extension of the N=40 Island of Inversion

The excited states of N=44 ^{74}Zn were investigated via γ-ray spectroscopy following ^{74}Cu β decay. By exploiting γ-γ angular correlation analysis, the 2_{2}^{+}, 3_{1}^{+}, 0_{2}^{+}, and 2_{3}^{+} states in ^{74}Zn were firmly established. The γ-ray branching and E2/M1 mixing ratios for transitions deexciting the 2_{2}^{+}, 3_{1}^{+}, and 2_{3}^{+} states were measured, allowing for the extraction of relative B(E2) values. In particular, the 2_{3}^{+}→0_{2}^{+} and 2_{3}^{+}→4_{1}^{+} transitions were observed for the first time. The results show excellent agreement with new microscopic large-scale shell-model calculations, and are discussed in terms of underlying shapes, as well as the role of neutron excitations across the N=40 gap. Enhanced axial shape asymmetry (triaxiality) is suggested to characterize ^{74}Zn in its ground state. Furthermore, an excited K=0 band with a significantly larger softness in its shape is identified. A shore of the N=40 "island of inversion" appears to manifest above Z=26, previously thought as its northern limit in the chart of the nuclides.

Nuclear moments of indium isotopes reveal abrupt change at magic number 82

In spite of the high-density and strongly correlated nature of the atomic nucleus, experimental and theoretical evidence suggests that around particular 'magic' numbers of nucleons, nuclear properties are governed by a single unpaired nucleon^1,2. A microscopic understanding of the extent of this behaviour and its evolution in neutron-rich nuclei remains an open question in nuclear physics^3-5. The indium isotopes are considered a textbook example of this phenomenon⁶, in which the constancy of their electromagnetic properties indicated that a single unpaired proton hole can provide the identity of a complex many-nucleon system^6,7. Here we present precision laser spectroscopy measurements performed to investigate the validity of this simple single-particle picture. Observation of an abrupt change in the dipole moment at N = 82 indicates that, whereas the single-particle picture indeed dominates at neutron magic number N = 82 (refs. ^2,8), it does not for previously studied isotopes. To investigate the microscopic origin of these observations, our work provides a combined effort with developments in two complementary nuclear many-body methods: ab initio valence-space in-medium similarity renormalization group and density functional theory (DFT). We find that the inclusion of time-symmetry-breaking mean fields is essential for a correct description of nuclear magnetic properties, which were previously poorly constrained. These experimental and theoretical findings are key to understanding how seemingly simple single-particle phenomena naturally emerge from complex interactions among protons and neutrons.

Electric Monopole Transition from the Superdeformed Band in ^{40}Ca

The electric monopole (E0) transition strength ρ^{2} for the transition connecting the third 0^{+} level, a "superdeformed" band head, to the "spherical" 0^{+} ground state in doubly magic ^{40}Ca is determined via e^{+}e^{-} pair-conversion spectroscopy. The measured value ρ^{2}(E0;0_{3}^{+}→0_{1}^{+})=2.3(5)×10^{-3} is the smallest ρ^{2}(E0;0^{+}→0^{+}) found in A<50 nuclei. In contrast, the E0 transition strength to the ground state observed from the second 0^{+} state, a band head of "normal" deformation, is an order of magnitude larger ρ^{2}(E0;0_{2}^{+}→0_{1}^{+})=25.9(16)×10^{-3}, which shows significant mixing between these two states. Large-scale shell-model (LSSM) calculations are performed to understand the microscopic structure of the excited states and the configuration mixing between them; experimental ρ^{2} values in ^{40}Ca and neighboring isotopes are well reproduced by the LSSM calculations. The unusually small ρ^{2}(E0;0_{3}^{+}→0_{1}^{+}) value is due to destructive interference in the mixing of shape-coexisting structures, which are based on several different multiparticle-multihole excitations. This observation goes beyond the usual treatment of E0 strengths, where two-state shape mixing cannot result in destructive interference.

Evidence of a sudden increase in the nuclear size of proton-rich silver-96

Understanding the evolution of the nuclear charge radius is one of the long-standing challenges for nuclear theory. Recently, density functional theory calculations utilizing Fayans functionals have successfully reproduced the charge radii of a variety of exotic isotopes. However, difficulties in the isotope production have hindered testing these models in the immediate region of the nuclear chart below the heaviest self-conjugate doubly-magic nucleus ¹⁰⁰Sn, where the near-equal number of protons (Z) and neutrons (N) lead to enhanced neutron-proton pairing. Here, we present an optical excursion into this region by crossing the N = 50 magic neutron number in the silver isotopic chain with the measurement of the charge radius of ⁹⁶Ag (N = 49). The results provide a challenge for nuclear theory: calculations are unable to reproduce the pronounced discontinuity in the charge radii as one moves below N = 50. The technical advancements in this work open the N = Z region below ¹⁰⁰Sn for further optical studies, which will lead to more comprehensive input for nuclear theory development.

Quasifree Neutron Knockout Reaction Reveals a Small s-Orbital Component in the Borromean Nucleus ^{17}B

A kinematically complete quasifree (p,pn) experiment in inverse kinematics was performed to study the structure of the Borromean nucleus ^{17}B, which had long been considered to have a neutron halo. By analyzing the momentum distributions and exclusive cross sections, we obtained the spectroscopic factors for 1s_{1/2} and 0d_{5/2} orbitals, and a surprisingly small percentage of 9(2)% was determined for 1s_{1/2}. Our finding of such a small 1s_{1/2} component and the halo features reported in prior experiments can be explained by the deformed relativistic Hartree-Bogoliubov theory in continuum, revealing a definite but not dominant neutron halo in ^{17}B. The present work gives the smallest s- or p-orbital component among known nuclei exhibiting halo features and implies that the dominant occupation of s or p orbitals is not a prerequisite for the occurrence of a neutron halo.

Ab Initio Neutrinoless Double-Beta Decay Matrix Elements for ^{48}Ca, ^{76}Ge, and ^{82}Se

We calculate basis-space converged neutrinoless ββ-decay nuclear matrix elements for the lightest candidates: ^{48}Ca, ^{76}Ge, and ^{82}Se. Starting from initial two- and three-nucleon forces, we apply the ab initio in-medium similarity renormalization group to construct valence-space Hamiltonians and consistently transformed ββ-decay operators. We find that the tensor component is non-negligible in ^{76}Ge and ^{82}Se, and the resulting nuclear matrix elements are overall 25%-45% smaller than those obtained from the phenomenological shell model. While a final matrix element with uncertainties still requires substantial developments, this work nevertheless opens a path toward a true first-principles calculation of neutrinoless ββ decay in all nuclei relevant for ongoing large-scale searches.

The impact of nuclear shape on the emergence of the neutron dripline

Atomic nuclei are composed of a certain number of protons Z and neutrons N. A natural question is how large Z and N can be. The study of superheavy elements explores the large Z limit^1,2, and we are still looking for a comprehensive theoretical explanation of the largest possible N for a given Z-the existence limit for the neutron-rich isotopes of a given atomic species, known as the neutron dripline³. The neutron dripline of oxygen (Z = 8) can be understood theoretically as the result of single nucleons filling single-particle orbits confined by a mean potential, and experiments confirm this interpretation. However, recent experiments on heavier elements are at odds with this description. Here we show that the neutron dripline from fluorine (Z = 9) to magnesium (Z = 12) can be predicted using a mechanism that goes beyond the single-particle picture: as the number of neutrons increases, the nuclear shape assumes an increasingly ellipsoidal deformation, leading to a higher binding energy. The saturation of this effect (when the nucleus cannot be further deformed) yields the neutron dripline: beyond this maximum N, the isotope is unbound and further neutrons 'drip' out when added. Our calculations are based on a recently developed effective nucleon-nucleon interaction⁴, for which large-scale eigenvalue problems are solved using configuration-interaction simulations. The results obtained show good agreement with experiments, even for excitation energies of low-lying states, up to the nucleus of magnesium-40 (which has 28 neutrons). The proposed mechanism for the formation of the neutron dripline has the potential to stimulate further thinking in the field towards explaining nucleosynthesis with neutron-rich nuclei.

Absence of Low-Energy Shape Coexistence in ^{80}Ge: The Nonobservation of a Proposed Excited 0_{2}^{+} Level at 639 keV

The ^{80}Ge structure was investigated in a high-statistics β-decay experiment of ^{80}Ga using the GRIFFIN spectrometer at TRIUMF-ISAC through γ, β-e, e-γ, and γ-γ spectroscopy. No evidence was found for the recently reported 0_{2}^{+} 639-keV level suggested as evidence for low-energy shape coexistence in ^{80}Ge. Large-scale shell model calculations performed in ^{78,80,82}Ge place the 0_{2}^{+} level in ^{80}Ge at 2 MeV. The new experimental evidence combined with shell model predictions indicate that low-energy shape coexistence is not present in ^{80}Ge.

Sequential Nature of (p,3p) Two-Proton Knockout from Neutron-Rich Nuclei

Twenty-one two-proton knockout (p,3p) cross sections were measured from neutron-rich nuclei at ∼250  MeV/nucleon in inverse kinematics. The angular distribution of the three emitted protons was determined for the first time, demonstrating that the (p,3p) kinematics are consistent with two sequential proton-proton collisions within the projectile nucleus. Ratios of (p,3p) over (p,2p) inclusive cross sections follow the trend of other many-nucleon removal reactions, further reinforcing the sequential nature of (p,3p) in neutron-rich nuclei.

Metastable States of ^{92,94}Se: Identification of an Oblate K Isomer of ^{94}Se and the Ground-State Shape Transition between N=58 and 60

Here we present new information on the shape evolution of the very neutron-rich ^{92,94}Se nuclei from an isomer-decay spectroscopy experiment at the Radioactive Isotope Beam Factory at RIKEN. High-resolution germanium detectors were used to identify delayed γ rays emitted following the decay of their isomers. New transitions are reported extending the previously known level schemes. The isomeric levels are interpreted as originating from high-K quasineutron states with an oblate deformation of β∼0.25, with the high-K state in ^{94}Se being metastable and K hindered. Following this, ^{94}Se is the lowest-mass neutron-rich nucleus known to date with such a substantial K hindrance. Furthermore, it is the first observation of an oblate K isomer in a deformed nucleus. This opens up the possibility for a new region of K isomers at low Z and at oblate deformation, involving the same neutron orbitals as the prolate orbitals within the classic Z∼72 deformed hafnium region. From an interpretation of the level scheme guided by theoretical calculations, an oblate deformation is also suggested for the ^{94}Se_{60} ground-state band.

Global Sensitivity Analysis of Bulk Properties of an Atomic Nucleus

We perform a global sensitivity analysis of the binding energy and the charge radius of the nucleus ^{16}O to identify the most influential low-energy constants in the next-to-next-to-leading order chiral Hamiltonian with two- and three-nucleon forces. For this purpose, we develop a subspace-projected coupled-cluster method using eigenvector continuation [Frame D. et al., Phys. Rev. Lett. 121, 032501 (2018)PRLTAO0031-900710.1103/PhysRevLett.121.032501]. With this method, we compute the binding energy and charge radius of ^{16}O at more than 10^{6} different values of the 16 low-energy constants in one hour on a standard laptop computer. For relatively small subspace projections, the root-mean-square error is about 1% compared to full-space coupled-cluster results. We find that 58(1)% of the variance in energy can be apportioned to a single contact term in the ^{3}S_{1} wave, whereas the radius depends sensitively on several low-energy constants and their higher-order correlations. The results identify the most important parameters for describing nuclear saturation and help prioritize efforts for uncertainty reduction of theoretical predictions. The achieved acceleration opens up an array of computational statistics analyses of the underlying description of the strong nuclear interaction in nuclei across the Segrè chart.

Quasifree Neutron Knockout from ^{54}Ca Corroborates Arising N=34 Neutron Magic Number

Exclusive cross sections and momentum distributions have been measured for quasifree one-neutron knockout reactions from a ^{54}Ca beam striking on a liquid hydrogen target at ∼200  MeV/u. A significantly larger cross section to the p_{3/2} state compared to the f_{5/2} state observed in the excitation of ^{53}Ca provides direct evidence for the nature of the N=34 shell closure. This finding corroborates the arising of a new shell closure in neutron-rich calcium isotopes. The distorted-wave impulse approximation reaction formalism with shell model calculations using the effective GXPF1Bs interaction and ab initio calculations concur our experimental findings. Obtained transverse and parallel momentum distributions demonstrate the sensitivity of quasifree one-neutron knockout in inverse kinematics on a thick liquid hydrogen target with the reaction vertex reconstructed to final state spin-parity assignments.

Multiple Shape Coexistence in ^{110,112}Cd

From detailed spectroscopy of ^{110}Cd and ^{112}Cd following the β^{+}/electron-capture decay of ^{110,112}In and the β^{-} decay of ^{112}Ag, very weak decay branches from nonyrast states are observed. The transition rates determined from the measured branching ratios and level lifetimes obtained with the Doppler-shift attenuation method following inelastic neutron scattering reveal collective enhancements that are suggestive of a series of rotational bands. In ^{110}Cd, a γ band built on the shape-coexisting intruder configuration is suggested. For ^{112}Cd, the 2^{+} and 3^{+} intruder γ-band members are suggested, the 0_{3}^{+} band is extended to spin 4^{+}, and the 0_{4}^{+} band is identified. The results are interpreted using beyond-mean-field calculations employing the symmetry conserving configuration mixing method with the Gogny D1S energy density functional and with the suggestion that the Cd isotopes exhibit multiple shape coexistence.

Probing the nuclear structure in the vicinity of 78Ni

Theoretical and experimental studies of neutron-rich nuclei have shown that the general concept of shell structure is not as robust and universal as earlier thought, but can exhibit significant changes as a function of neutron excess. New magic numbers appear and some other conventional ones disappear mainly because of a different ordering of the single-particle orbitals. In the present contribution, recent experimental studies of neutron-rich Cu isotopes, performed at RIKEN using β decay and one-proton knockout reactions, will be discussed. Neutron-rich nuclei near 78Ni were populated through in-flight fission of 238U on thick 9Be targets in both experiments. In the β-decay study, 75,77Ni nuclei were implanted into the WAS3ABi silicon array, while γ rays from excited states in 75,77Cu emitted after β decay of the implanted ions were detected with the EURICA Ge detector array that was surrounding the active stopper. In a second experiment within the SEASTAR campaign at RIKEN, the same 75,77Cu nuclei were produced in (p,2p) knockout reactions from 76,78Zn beam particles at around 250 MeV/nucleon impinging onto the MINOS liquid hydrogen target. In the latter experiment the DALI2 NaI array was used to detect de-excitation γ rays measured in coincidence with Cu nuclei identified in the Zero Degree Spectrometer. Both studies are complimentary and greatly contribute to our understanding on the nuclear structure in the 78Ni region.