Search

The symmetry energy and its density dependence are crucial inputs for many nuclear physics and astrophysics applications, as they determine properties ranging from the neutron-skin thickness of nuclei to the crust thickness and the radius of neutron stars. Recently, PREX-II reported a value of $0.283\ifmmode\pm\else\textpm\fi{}0.071\text{ }\text{ }\mathrm{fm}$ for the neutron-skin thickness of $^{208}\mathrm{Pb}$, implying a slope parameter $L=106\ifmmode\pm\else\textpm\fi{}37\text{ }\text{ }\mathrm{MeV}$, larger than most ranges obtained from microscopic calculations and other nuclear experiments. We use a nonparametric equation of state representation based on Gaussian processes to constrain the symmetry energy ${S}_{0}$, $L$, and ${R}_{\mathrm{skin}}^{^{208}\mathrm{Pb}}$ directly from observations of neutron stars with minimal modeling assumptions. The resulting astrophysical constraints from heavy pulsar masses, LIGO/Virgo, and NICER clearly favor smaller values of the neutron skin and $L$, as well as negative symmetry incompressibilities. Combining astrophysical data with PREX-II and chiral effective field theory constraints yields ${S}_{0}=33.{0}_{\ensuremath{-}1.8}^{+2.0}\text{ }\text{ }\mathrm{MeV}$, $L=5{3}_{\ensuremath{-}15}^{+14}\text{ }\text{ }\mathrm{MeV}$, and ${R}_{\mathrm{skin}}^{^{208}\mathrm{Pb}}=0.1{7}_{\ensuremath{-}0.04}^{+0.04}\text{ }\text{ }\mathrm{fm}$.

Understanding the dense matter equation of state at extreme conditions is an important open problem. Astrophysical observations of neutron stars promise to solve this, with NICER poised to make precision measurements of mass and radius for several stars using the waveform modelling technique. What has been less clear, however, is how these mass-radius measurements might translate into equation of state constraints and what are the associated equation of state sensitivities. We use Bayesian inference to explore and contrast the constraints that would result from different choices for the equation of state parametrization; comparing the well-established piecewise polytropic parametrization to one based on physically motivated assumptions for the speed of sound in dense matter. We also compare the constraints resulting from Bayesian inference to those from simple compatibility cuts. We find that the choice of equation of state parametrization and particularly its prior assumptions can have a significant effect on the inferred global mass-radius relation and the equation of state constraints. Our results point to important sensitivities when inferring neutron star and dense matter properties. This applies also to inferences from gravitational wave observations., to appear in MNRAS

Steven Weinberg's seminal papers from 1990-92 initiated the use of effective field theories (EFTs) for nuclei. We summarize progress, priorities, and open questions for nuclear EFT developments based on the 2019 INT program "Nuclear Structure at the Crossroads.", Comment: 15 pages, Contribution to Special Issue: Celebrating 30 years of Steven Weinberg's papers on Nuclear Forces from Chiral Lagrangians; published version

We discuss high-order calculations in perturbative effective field theory for fermions at low energy scales. The Fermi-momentum or $k_{\rm F} a_s$ expansion for the ground-state energy of the dilute Fermi gas is calculated to fourth order, both in cutoff regularization and in dimensional regularization. For the case of spin one-half fermions we find from a Bayesian analysis that the expansion is well-converged at this order for ${| k_{\rm F} a_s | \lesssim 0.5}$. Further, we show that Pad{\'e}-Borel resummations can improve the convergence for ${| k_{\rm F} a_s | \lesssim 1}$. Our results provide important constraints for nonperturbative calculations of ultracold atoms and dilute neutron matter., Comment: 18 pages, 5 figures, minor changes, published version

A direct measurement of the decay width of the excited ${0}_{1}^{+}$ state of $^{6}\mathrm{Li}$ using the relative self-absorption technique is reported. Our value of ${\mathrm{\ensuremath{\Gamma}}}_{\ensuremath{\gamma},{0}_{1}^{+}\ensuremath{\rightarrow}{1}_{1}^{+}}=8.17(14{)}_{\mathrm{stat}.}(11{)}_{\mathrm{syst}.}\text{ }\text{ }\mathrm{eV}$ provides sufficiently low experimental uncertainties to test modern theories of nuclear forces. The corresponding transition rate is compared to the results of ab initio calculations based on chiral effective field theory that take into account contributions to the magnetic dipole operator beyond leading order. This enables a precision test of the impact of two-body currents that enter at next-to-leading order.

Over the past decade the in-medium similarity renormalization group (IMSRG) approach has proven to be a powerful and versatile ab initio many-body method for studying medium-mass nuclei. So far, the IMSRG was limited to the approximation in which only up to two-body operators are incorporated in the renormalization group flow, referred to as the IMSRG(2). In this work, we extend the IMSRG(2) approach to fully include three-body operators yielding the IMSRG(3) approximation. We use a perturbative scaling analysis to estimate the importance of individual terms in this approximation and introduce truncations that aim to approximate the IMSRG(3) at a lower computational cost. The IMSRG(3) is systematically benchmarked for different nuclear Hamiltonians for ${}^{4}\text{He}$ and ${}^{16}\text{O}$ in small model spaces. The IMSRG(3) systematically improves over the IMSRG(2) relative to exact results. Approximate IMSRG(3) truncations constructed based on computational cost are able to reproduce much of the systematic improvement offered by the full IMSRG(3). We also find that the approximate IMSRG(3) truncations behave consistently with expectations from our perturbative analysis, indicating that this strategy may also be used to systematically approximate the IMSRG(3)., 23 pages, 10 figures, published version

We present new equations of state for applications in core-collapse supernova and neutron star merger simulations. We start by introducing an effective mass parametrization that is fit to recent microscopic calculations up to twice saturation density. This is important to capture the predicted thermal effects, which have been shown to determine the proto-neutron star contraction in supernova simulations. The parameter range of the energy-density functional underlying the equation of state is constrained by chiral effective field theory results at nuclear densities as well as by functional renormalization group computations at high densities based on QCD. We further implement observational constraints from measurements of heavy neutron stars, the gravitational wave signal of GW170817, and from the recent NICER results. Finally, we study the resulting allowed ranges for the equation of state and for properties of neutron stars, including the predicted ranges for the neutron star radius and maximum mass., Comment: 21 pages, 19 figures, minor changes, published version

The nuclear many-body problem for medium-mass systems is commonly addressed using wave-function expansion methods that build upon a second-quantized representation of many-body operators with respect to a chosen computational basis. While various options for the computational basis are available, perturbatively constructed natural orbitals recently have been shown to lead to significant improvement in many-body applications yielding faster model-space convergence and lower sensitivity to basis set parameters in large-scale no-core shell model diagonalizations. This work provides a detailed comparison of single-particle basis sets and a systematic benchmark of natural orbitals in nonperturbative many-body calculations using the in-medium similarity renormalization group approach. As a key outcome we find that the construction of natural orbitals in a large single-particle basis enables for performing the many-body calculation in a reduced space of much lower dimension, thus offering significant computational savings in practice that help extend the reach of ab initio methods towards heavier masses and higher accuracy., Comment: 15 pages, 9 figures, published version

In recent years our understanding of the dense matter equation of state (EOS) of neutron stars has significantly improved by analyzing multimessenger data from radio/X-ray pulsars, gravitational wave events, and from nuclear physics constraints. Here we study the additional impact on the EOS from the jointly estimated mass and radius of PSR J0740+6620, presented in Riley et al. (2021) by analyzing a combined dataset from X-ray telescopes NICER and XMM-Newton. We employ two different high-density EOS parameterizations: a piecewise-polytropic (PP) model and a model based on the speed of sound in a neutron star (CS). At nuclear densities these are connected to microscopic calculations of neutron matter based on chiral effective field theory interactions. In addition to the new NICER data for this heavy neutron star, we separately study constraints from the radio timing mass measurement of PSR J0740+6620, the gravitational wave events of binary neutron stars GW190425 and GW170817, and for the latter the associated kilonova AT2017gfo. By combining all these, and the NICER mass-radius estimate of PSR J0030+0451 we find the radius of a 1.4 solar mass neutron star to be constrained to the 95% credible ranges 12.33^{+0.76}_{-0.81} km (PP model) and 12.18^{+0.56}_{-0.79} km (CS model). In addition, we explore different chiral effective field theory calculations and show that the new NICER results provide tight constraints for the pressure of neutron star matter at around twice saturation density, which shows the power of these observations to constrain dense matter interactions at intermediate densities., 17 pages, 8 figures; accepted for publication in the Astrophysical Journal Letters

100Sn is of singular interest for nuclear structure. Its closed-shell proton and neutron configuration exhibit exceptional binding and 100Sn is the heaviest nucleus comprising protons and neutrons in equal number, a feature that enhances the contribution of the short-range, proton-neutron pairing interaction and strongly influences its decay via the weak interaction. Decays studies in the region of 100Sn have attempted to prove its doubly magic character but few have studied it from the ab initio theoretical perspective and none have addressed the odd-proton nuclear forces. Here we present, the first direct measurement of the exotic odd-proton nuclide 100In - the beta-decay daughter of 100Sn - and 99In, only one proton below 100Sn. The most advanced mass spectrometry techniques were used to measure 99In, produced at a rate of only a few ions per second, and to resolve the ground and isomeric states in 101In. The experimental results are confronted with new ab initio many-body approaches. The 100-fold improvement in precision of the 100In mass value exarcebates a striking discrepancy in the atomic mass values of 100Sn deduced from recent beta-decay results., 12 pages, 4 figures, 1 Table

We develop a method that uses truncation-order-dependent re-expansions constrained by generic strong-coupling information to extrapolate perturbation series to the nonperturbative regime. The method is first benchmarked against a zero-dimensional model field theory and then applied to the dilute Fermi gas in one and three dimensions. Overall, our method significantly outperforms Pad\'e and Borel extrapolations in these examples. The results for the ground-state energy of the three-dimensional Fermi gas are robust with respect to changes of the form of the re-expansion and compare well with quantum Monte Carlo simulations throughout the BCS regime and beyond., Comment: minor improvements, matches published version; 11 pages, 8 figures

We study the equation of state of neutron matter at finite temperature based on two- and three-nucleon interactions derived within chiral effective field theory to next-to-next-to-next-to-leading order. The free energy, pressure, entropy, and internal energy are calculated using many-body perturbation theory including terms up to third order around the self-consistent Hartree-Fock solution. We include contributions from three-nucleon interactions without employing the normal-ordering approximation and provide theoretical uncertainty estimates based on an order-by-order analysis in the chiral expansion. Our results demonstrate that thermal effects can be captured remarkably well via a thermal index and a density-dependent effective mass. The presented framework provides the basis for studying the dense matter equation of state at general temperatures and proton fractions relevant for core-collapse supernovae and neutron star mergers., 14 pages, 14 figures, minor changes, published version

The structure of exotic nuclei provides valuable tests for state-of-the-art nuclear theory. In particular electromagnetic transition rates are more sensitive to aspects of nuclear forces and many-body physics than excitation energies alone. We report the first lifetime measurement of excited states in $^{21}$O, finding $\tau_{1/2^+}=420^{+35}_{-32}\text{(stat)}^{+34}_{-12}\text{(sys)}$\,ps. This result together with the deduced level scheme and branching ratio of several $\gamma$-ray decays are compared to both phenomenological shell-model and ab initio calculations based on two- and three-nucleon forces derived from chiral effective field theory. We find that the electric quadrupole reduced transition probability of $\rm B(E2;1/2^+ \rightarrow 5/2^+_{g.s.}) = 0.71^{+0.07\ +0.02}_{-0.06\ -0.06}$~e$^2$fm$^4$, derived from the lifetime of the $1/2^+$ state, is smaller than the phenomenological result where standard effective charges are employed, suggesting the need for modifications of the latter in neutron-rich oxygen isotopes. We compare this result to both large-space and valence-space ab initio calculations, and by using multiple input interactions we explore the sensitivity of this observable to underlying details of nuclear forces., Comment: 23 pages, 3 figures

We explore constraints on the equation of state (EOS) of neutron-rich matter based on microscopic calculations up to nuclear densities and observations of neutron stars. In a previous work we showed that predictions based on modern nuclear interactions derived within chiral effective field theory and the observation of two-solar-mass neutron stars result in a robust uncertainty range for neutron star radii and the EOS over a wide range of densities. In this work we extend this study, employing both the piecewise polytrope extension from Hebeler et al. as well as the speed of sound model of Greif et al., and show that moment of inertia measurements of neutron stars can significantly improve the constraints on the EOS and neutron star radii., 11 pages, 10 figures, published version

We present a determination of optical potentials using the double-folding method based on chiral effective field theory nucleon-nucleon interactions at next-to-next-to-leading order combined with dispersion relations to constrain the imaginary part. This approach is benchmarked on O16-O16 collisions, and extended to the C12-C12 and C12-O16 cases. Predictions derived from these potentials are compared to data for elastic scattering at energies up to 1000 MeV, as well as for fusion at low energy. Without adjusting parameters, excellent agreement with experiment is found. In addition, we study the sensitivity of the corresponding cross sections to the nucleon-nucleon interactions and nuclear densities used., SCOPUS: ar.j, info:eu-repo/semantics/published

The strength of the $N=28$ magic number in neutron-rich argon isotopes is examined through high-precision mass measurements of $^{46\text{--}48}\mathrm{Ar}$, performed with the ISOLTRAP mass spectrometer at ISOLDE/CERN. The new mass values are up to 90 times more precise than previous measurements. While they suggest the persistence of the $N=28$ shell closure for argon, we show that this conclusion has to be nuanced in light of the wealth of spectroscopic data and theoretical investigations performed with the SDPF-U phenomenological shell model interaction. Our results are also compared with ab initio calculations using the valence space in-medium similarity renormalization group and the self-consistent Green's function approaches. Both calculations provide a very good account of mass systematics at and around $Z=18$ and, generally, a consistent description of the physics in this region. This combined analysis indicates that $^{46}\mathrm{Ar}$ is the transition between the closed-shell $^{48}\mathrm{Ca}$ and collective $^{44}\mathrm{S}$.

The NICER collaboration recently published a joint estimate of the mass and the radius of PSR J0030+0451, derived via X-ray pulse-profile modeling. Raaijmakers et al. (2019) explored the implications of this measurement for the dense matter equation of state (EOS) using two parameterizations of the high-density EOS: a piecewise-polytropic model, and a model based on the speed of sound in neutron stars. In this work we obtain further constraints on the EOS following this approach, but we also include information about the tidal deformability of neutron stars from the gravitational wave signal of the compact binary merger GW170817. We compare the constraints on the EOS to those set by the recent measurement of a 2.14 solar mass pulsar, included as a likelihood function approximated by a Gaussian, and find a small increase in information gain. To show the flexibility of our method, we also explore the possibility that GW170817 was a neutron star-black hole merger, which yields weaker constraints on the EOS., 18 pages, 8 figures. Accepted for publication in ApJ Letters

We present the first laser spectroscopic measurement of the neutron-rich nucleus ^{68}Ni at the N=40 subshell closure and extract its nuclear charge radius. Since this is the only short-lived isotope for which the dipole polarizability α_{D} has been measured, the combination of these observables provides a benchmark for nuclear structure theory. We compare them to novel coupled-cluster calculations based on different chiral two- and three-nucleon interactions, for which a strong correlation between the charge radius and dipole polarizability is observed, similar to the stable nucleus ^{48}Ca. Three-particle-three-hole correlations in coupled-cluster theory substantially improve the description of the experimental data, which allows to constrain the neutron radius and neutron skin of ^{68}Ni.

We probe the $N=82$ nuclear shell closure by mass measurements of neutron-rich cadmium isotopes with the ISOLTRAP spectrometer at ISOLDE-CERN. The new mass of $^{132}\mathrm{Cd}$ offers the first value of the $N=82$, two-neutron shell gap below $Z=50$ and confirms the phenomenon of mutually enhanced magicity at $^{132}\mathrm{Sn}$. Using the recently implemented phase-imaging ion-cyclotron-resonance method, the ordering of the low-lying isomers in $^{129}\mathrm{Cd}$ and their energies are determined. The new experimental findings are used to test large-scale shell-model, mean-field, and beyond-mean-field calculations, as well as the ab initio valence-space in-medium similarity renormalization group.

The cross section for coherent elastic neutrino-nucleus scattering (CE$\nu$NS) depends on the response of the target nucleus to the external current, in the Standard Model (SM) mediated by the exchange of a $Z$ boson. This is typically subsumed into an object called the weak form factor of the nucleus. Here, we provide results for this form factor calculated using the large-scale nuclear shell model for a wide range of nuclei of relevance for current CE$\nu$NS experiments, including cesium, iodine, argon, fluorine, sodium, germanium, and xenon. In addition, we provide the responses needed to capture the axial-vector part of the cross section, which does not scale coherently with the number of neutrons, but may become relevant for the SM prediction of CE$\nu$NS on target nuclei with nonzero spin. We then generalize the formalism allowing for contributions beyond the SM. In particular, we stress that in this case, even for vector and axial-vector operators, the standard weak form factor does not apply anymore, but needs to be replaced by the appropriate combination of the underlying nuclear structure factors. We provide the corresponding expressions for vector, axial-vector, but also (pseudo-)scalar, tensor, and dipole effective operators, including two-body-current effects as predicted from chiral effective field theory. Finally, we update the spin-dependent structure factors for dark matter scattering off nuclei according to our improved treatment of the axial-vector responses., Comment: 28 pages, 11 figures; journal version

7 pags., 4 figs., 1 tab., Excited states in the N=40 isotone Ti were populated via the V(p,2p)Ti reaction at ∼200 MeV/nucleon at the Radioactive Isotope Beam Factory and studied using γ-ray spectroscopy. The energies of the 2 →0 and 4 →2 transitions, observed here for the first time, indicate a deformed Ti ground state. These energies are increased compared to the neighboring Cr and Fe isotones, suggesting a small decrease of quadrupole collectivity. The present measurement is well reproduced by large-scale shell-model calculations based on effective interactions, while ab initio and beyond mean-field calculations do not yet reproduce our findings. The shell-model calculations for Ti show a dominant configuration with four neutrons excited across the N=40 gap. Likewise, they indicate that the N=40 island of inversion extends down to Z=20, disfavoring a possible doubly magic character of the elusive Ca., We thank the RIKEN Nishina Center accelerator staff and the Bi-gRIPS team for the stable operation of the high-intensity Zn beam and for the preparation of the secondary beam setting. K.O. ac-knowledges the support by Grant-in-Aid for Scientific Research of the Japan Society for the Promotion of Science (JSPS) JP16K05352. A.P. is supported in part by the Ministerio de Ciencia, Innovación y Universidades (Spain), Severo Ochoa Programme SEV-2016-0597 and grant PGC-2018-94583. F.B. is supported by the RIKEN Spe-cial Postdoctoral Researcher Program. L.X.C. and B.D.L. would like to thank the Vietnam Ministry of Science and Technology (MOST) for its support through the Physics Development Program Grant No. ÐTÐLCN.25/18. I.G. has been supported by HIC for FAIR and Croatian Science Foundation under projects no. 1257 and 7194. D. So. was supported by the the European Regional Develop-ment Fund contract No. GINOP-2.3.3-15-2016-00034 and the National Research, Development and Innovation Fund of Hungary via Project No. K128947. V.V. acknowledges support from the Span-ish Ministerio de Economía y Competitividad under Contract No. FPA2017-84756-C4-2-P. K.I.H., D.K. and S.Y.P. acknowledge the sup-port from the National Research Foundation of Korea grant No. 2018R1A5A1025563 and 2019M7A1A1033186. The development of MINOS was supported by the European Research Council through the ERC Grant No. MINOS-258567. This work was also supported by the JSPS KAKENHI Grant No. 18K03639, MEXT as “Priority is-sue on post-K computer” (Elucidation of the fundamental laws and evolution of the universe), the Joint Institute for Computational Fundamental Science (JICFuS), the CNS-RIKEN joint project for large-scale nuclear structure calculations, Natural Sciences and Engineering Research Council (NSERC) of Canada, the Deutsche Forschungsgemeinschaft – Projektnummer 279384907 – SFB 1245, the PRISMA Cluster of Excellence, and the BMBF under Contracts No. 05P18RDFN1 and 05P19RDFN1. TRIUMF receives funding via a contribution through the National Research Council Canada. Com-putations were performed at the Jülich Supercomputing Center (JURECA)

Empirical energy density functionals (EDFs) are generally successful in describing nuclear properties across the table of nuclides. But their limitations motivate using the density-matrix expansion (DME) to embed long-range pion interactions into a Skyrme functional. Recent results on the impact of the pion were both encouraging and puzzling, necessitating a careful re-examination of the DME implementation. Here we take the first steps, focusing on two-body scalar terms in the DME. Exchange energies with long-range one-pion contributions are well approximated by all DME implementations considered, with preference for variants that do not truncate at two derivatives in every EDF term. The use of the DME for chiral pion contributions is therefore supported by this investigation. For scalar-isovector energies it is important to treat neutrons and protons separately. The results are found to apply under broad conditions, although self-consistency is not yet tested., Comment: 21 pages, 15 figures, including Supplemental Material, published version

First principles calculations of atomic nuclei based on microscopic nuclear forces derived from chiral effective field theory (EFT) have blossomed in the past years. A key element of such ab initio studies is the understanding and quantification of systematic and statistical errors arising from the omission of higher-order terms in the chiral expansion as well as the model calibration. While there has been significant progress in analyzing theoretical uncertainties for nucleon-nucleon scattering observables, the generalization to multi-nucleon systems has not been feasible yet due to the high computational cost of evaluating observables for a large set of low-energy couplings. In this Letter we show that a new method called eigenvector continuation (EC) can be used for constructing an efficient and accurate emulator for nuclear many-body observables, thereby enabling uncertainty quantification in multi-nucleon systems. We demonstrate the power of EC emulation with a proof-of-principle calculation that lays out all correlations between bulk ground-state observables in the few-nucleon sector. On the basis of ab initio calculations for the ground-state energy and radius in 4He, we demonstrate that EC is more accurate and efficient compared to established methods like Gaussian processes., Comment: 8 pages, 6 figures, Python code and input files provided as ancillary material, published version

Quantum many-body theory has witnessed tremendous progress in various fields, ranging from atomic and solid-state physics to quantum chemistry and nuclear structure. Due to the inherent computational burden linked to the ab initio treatment of microscopic fermionic systems, it is desirable to obtain accurate results through low-order perturbation theory. In atomic nuclei however, effects such as strong short-range repulsion between nucleons can spoil the convergence of the expansion and make the reliability of perturbation theory unclear. Mathematicians have devised an extensive machinery to overcome the problem of divergent expansions by making use of so-called resummation methods. In large-scale many-body applications such schemes are often of limited use since no a priori analytical knowledge of the expansion is available. We present here eigenvector continuation as an alternative resummation tool that is both efficient and reliable because it is based on robust and simple mathematical principles., 6 pages, 2 figures

We exploit the many-body self-consistent Green's function method to analyze finite-temperature properties of infinite nuclear matter and to explore the behavior of the thermal index used to simulate thermal effects in equations of state for astrophysical applications. We show how the thermal index is both density and temperature dependent, unlike often considered, and we provide an error estimate based on our ${\it ab~initio}$ calculations. The inclusion of many-body forces is found to be critical for the density dependence of the thermal index. We also compare our results to a parametrization in terms of the density dependence of the nucleon effective mass. Our study questions the validity of predictions made for the gravitational-wave signal from neutron-star merger simulations with a constant thermal index., 11 pages, 7 figures, added references

We study ground-state energies and charge radii of closed-shell medium-mass nuclei based on novel chiral nucleon-nucleon (NN) and three-nucleon (3N) interactions, with a focus on exploring the connections between finite nuclei and nuclear matter. To this end, we perform in-medium similarity renormalization group (IM-SRG) calculations based on chiral interactions at next-to-leading order (NLO), N$^2$LO, and N$^3$LO, where the 3N interactions at N$^2$LO and N$^3$LO are fit to the empirical saturation point of nuclear matter and to the triton binding energy. Our results for energies and radii at N$^2$LO and N$^3$LO overlap within uncertainties, and the cutoff variation of the interactions is within the EFT uncertainty band. We find underbound ground-state energies, as expected from the comparison to the empirical saturation point. The radii are systematically too large, but the agreement with experiment is better. We further explore variations of the 3N couplings to test their sensitivity in nuclei. While nuclear matter at saturation density is quite sensitive to the 3N couplings, we find a considerably weaker dependence in medium-mass nuclei. In addition, we explore a consistent momentum-space SRG evolution of these NN and 3N interactions, exhibiting improved many-body convergence. For the SRG-evolved interactions, the sensitivity to the 3N couplings is found to be stronger in medium-mass nuclei., 10 pages, 11 figures, published version

We investigate how a phase transition from neutron-star matter to spin-polarized neutron matter affects the equation of state and mass-radius relation of neutron stars. While general extension schemes for the equation of state allow for high pressures inside neutron stars, we find that a phase transition to spin-polarized neutron matter excludes extreme regimes. Hence, such a transition limits the maximum mass of neutron stars to lie below 2.6$-$$2.9 \, M_{\odot}$, depending on the microscopic nuclear forces used, while significantly larger masses could be reached without these constraints. These limits are in good agreement with recent constraints extracted from the neutron-star merger GW170817 and its electromagnetic counterpart. Assuming the description in terms of spin-polarized neutron matter to be valid in the center of neutron stars, we find that stars with a large spin-polarized domain in their core are ruled out by GW170817., 10 pages, 3 figures, published version

We reply to a Comment on our Letter [Phys. Rev. Lett. 118, 232501 (2017), arXiv:1612.01502] by A. Deltuva and R. Lazauskas [Phys. Rev. Lett 123, 069201 (2019), arXiv:1904.00925].

Chiral effective field theory allows one to calculate the response of few-nucleon systems to external currents, both for currents that can be probed in the Standard Model and ones that only exist in Standard-Model extensions. In combination with state-of-the-art many-body methods, the constraints from chiral symmetry can then be implemented in nuclear structure factors that describe the response of atomic nuclei in direct-detection searches for dark matter. We review the present status of this approach, including the role of coherently enhanced two-body currents, the discrimination of dark matter candidates based on the nuclear response functions, and limits on Higgs-portal dark matter., 10 pages, 6 figures, proceedings for The 9th International Workshop on Chiral Dynamics, Durham, North Carolina, USA, September 17-21, 2018

$\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

The first γ-ray spectroscopy of ^{52}Ar, with the neutron number N=34, was measured using the ^{53}K(p,2p) one-proton removal reaction at ∼210 MeV/u at the RIBF facility. The 2_{1}^{+} excitation energy is found at 1656(18) keV, the highest among the Ar isotopes with N>20. This result is the first experimental signature of the persistence of the N=34 subshell closure beyond ^{54}Ca, i.e., below the magic proton number Z=20. Shell-model calculations with phenomenological and chiral-effective-field-theory interactions both reproduce the measured 2_{1}^{+} systematics of neutron-rich Ar isotopes, and support a N=34 subshell closure in ^{52}Ar.

High-energy scattering processes, such as deep inelastic scattering (DIS) and quasielastic (QE) scattering provide a wealth of information about the structure of atomic nuclei. The remarkable discovery of the empirical linear relationship between the slope of the European Muon Collaboration (EMC) effect in DIS and the short-range-correlation (SRC) scaling factors $a_2$ in QE kinematics is naturally explained in terms of scale separation in effective field theory. This explanation has powerful consequences, allowing us to calculate and predict SRC scaling factors from ab initio low-energy nuclear theory. We present ab initio calculations of SRC scaling factors for a nucleus $A$ relative to the deuteron $a_2(A/d)$ and relative to $^3\rm He$ $a_2(A/^3\rm He)$ in light and medium-mass nuclei. Our framework further predicts that the EMC effect and SRC scaling factors have minimal or negligible isovector corrections., Comment: 25 pages, 6 figures, 2 tables

In this White Paper we present the potential of the Enhanced X-ray Timing and Polarimetry (eXTP) mission for determining the nature of dense matter; neutron star cores host an extreme density regime which cannot be replicated in a terrestrial laboratory. The tightest statistical constraints on the dense matter equation of state will come from pulse profile modelling of accretion-powered pulsars, burst oscillation sources, and rotation-powered pulsars. Additional constraints will derive from spin measurements, burst spectra, and properties of the accretion flows in the vicinity of the neutron star. Under development by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s., Comment: Accepted for publication on Sci. China Phys. Mech. Astron. (2019)

Nuclear magic numbers, which emerge from the strong nuclear force based on quantum chromodynamics, correspond to fully occupied energy shells of protons, or neutrons inside atomic nuclei. Doubly magic nuclei, with magic numbers for both protons and neutrons, are spherical and extremely rare across the nuclear landscape. While the sequence of magic numbers is well established for stable nuclei, evidence reveals modifications for nuclei with a large proton-to-neutron asymmetry. Here, we provide the first spectroscopic study of the doubly magic nucleus $^{78}$Ni, fourteen neutrons beyond the last stable nickel isotope. We provide direct evidence for its doubly magic nature, which is also predicted by ab initio calculations based on chiral effective field theory interactions and the quasi-particle random-phase approximation. However, our results also provide the first indication of the breakdown of the neutron magic number 50 and proton magic number 28 beyond this stronghold, caused by a competing deformed structure. State-of-the-art phenomenological shell-model calculations reproduce this shape coexistence, predicting further a rapid transition from spherical to deformed ground states with $^{78}$Ni as turning point.

The hyperfine coupling constants of neutron deficient $^{37}\mathrm{Ca}$ were deduced from the atomic hyperfine spectrum of the $4s\phantom{\rule{0.16em}{0ex}}^{2}S_{1/2}\phantom{\rule{0.16em}{0ex}}\ensuremath{\leftrightarrow}\phantom{\rule{0.16em}{0ex}}4p\phantom{\rule{0.16em}{0ex}}^{2}P_{3/2}$ transition in Ca ii, measured using the collinear laser spectroscopy technique. The ground-state magnetic-dipole and spectroscopic electric-quadrupole moments were determined for the first time as $\ensuremath{\mu}=+0.7453(72){\ensuremath{\mu}}_{N}$ and $Q=\ensuremath{-}15(11)\phantom{\rule{0.16em}{0ex}}{e}^{2}\phantom{\rule{0.16em}{0ex}}{\mathrm{fm}}^{2}$, respectively. The experimental values agree well with nuclear shell-model calculations using the universal sd model-space Hamiltonians versions A and B (USDA/B) in the $sd$-model space with a 95% probability of the canonical nucleon configuration. It is shown that the magnetic moment of $^{39}\mathrm{Ca}$ requires a larger non-$sd$-shell component than that of $^{37}\mathrm{Ca}$ for good agreement with the shell-model calculation, indicating a more robust closed subshell structure of $^{36}\mathrm{Ca}$ at the neutron number $N=16$ than $^{40}\mathrm{Ca}$. The results are also compared to valence-space in-medium similarity renormalization group calculations based on chiral two- and three-nucleon interactions.

We predict the limits of existence of atomic nuclei, the proton and neutron drip lines, from the light through medium-mass regions. Starting from a chiral two- and three-nucleon interaction with good saturation properties, we use the valence-space in-medium similarity renormalization group to calculate ground-state and separation energies from helium to iron, nearly 700 isotopes in total. We use the available experimental data to quantify the theoretical uncertainties for our ab initio calculations towards the drip lines. Where the drip lines are known experimentally, our predictions are consistent within the estimated uncertainty. For the neutron-rich sodium to chromium isotopes, we provide predictions to be tested at rare-isotope beam facilities., Comment: 6 pages, 3 figures, supplementary material included; accepted for publication in Physical Review Letters

Using effective field theory methods, we calculate for the first time the complete fourth-order term in the Fermi-momentum or $k_{\rm F} a_s$ expansion for the ground-state energy of a dilute Fermi gas. The convergence behavior of the expansion is examined for the case of spin one-half fermions and compared against quantum Monte-Carlo results, showing that the Fermi-momentum expansion is well-converged at this order for $| k_{\rm F} a_s | \lesssim 0.5$., 6 pages, 2 figures; v2: minor improvements in text, references updated, matches published version

We present nuclear structure factors that describe the generalized spin-independent coupling of weakly interacting massive particles (WIMPs) to nuclei. Our results are based on state-of-the-art nuclear structure calculations using the large-scale nuclear shell model. Starting from quark- and gluon-level operators, we consider all possible coherently enhanced couplings of spin-1/2 and spin-0 WIMPs to one and two nucleons up to third order in chiral effective field theory. This includes a comprehensive discussion of the structure factors corresponding to the leading two-nucleon currents covering, for the first time, the contribution of spin-2 operators. We provide results for the most relevant nuclear targets considered in present and planned dark matter direct detection experiments: fluorine, silicon, argon, and germanium, complementing our previous work on xenon. All results are also publicly available in a Python notebook., 24 pages, 19 figures, Python notebook available at https://theorie.ikp.physik.tu-darmstadt.de/strongint/ChiralEFT4DM.html; further details on nuclear structure added, journal version

In the two decades since the appearance of the book \Bose{Einstein Con- densation" in 1995, there have been a number of developments in our under- standing of dense matter. After a brief overview of neutron star structure and the Bose{Einstein condensed phases that have been proposed, we de- scribe selected topics, including neutron and proton pairing gaps, the physics of the inner crust of neutron stars, where a neutron uid penetrates a lattice

We construct valence-space Hamiltonians for use in shell-model calculations, where the residual two-body interaction is based on symmetry principles and the low-momentum expansion from chiral effective field theory. In addition to the usual free-space contact interactions, we also include novel center-of-mass--dependent operators that arise due to the Galilean invariance breaking by in-medium effects. We fitted the low-energy constants to 441 ground- and excited-state energies in the sd shell and obtained a root-mean-square derivation of 1.8 MeV at leading order and of 0.5 MeV at next-to-leading order, with natural low-energy constants in all cases. The developed chiral shell-model interactions enable order-by-order uncertainty estimates and show promising predictions for neutron-rich isotopes beyond the fitted data set., 19 pages, 12 figures, published version

We study systems of bosons and fermions in finite periodic boxes and show how the existence and properties of few-body resonances can be extracted from studying the volume dependence of the calculated energy spectra. Using a plane-wave-based discrete variable representation to conveniently implement periodic boundary conditions, we establish that avoided level crossings occur in the spectra of up to four particles and can be linked to the existence of multi-body resonances. To benchmark our method we use two-body calculations, where resonance properties can be determined with other methods, as well as a three-boson model interaction known to generate a three-boson resonance state. Finding good agreement for these cases, we then predict three-body and four-body resonances for models using a shifted Gaussian potential. Our results establish few-body finite-volume calculations as a new tool to study few-body resonances. In particular, the approach can be used to study few-neutron systems, where such states have been conjectured to exist., 13 pages, 10 figures, 2 tables, published version

The determination of nucleus–nucleus potentials is important not only to describe the properties of the colliding system, but also to extract nuclear-structure information and for modelling nuclear reactions for astrophysics. We present the first determination of double-folding potentials based on chiral effective field theory at leading, next-to-leading, and next-to-next-to-leading order. To this end, we construct new soft local chiral effective field theory interactions. We benchmark this approach in the 16O–16O system, and present results for cross sections computed for elastic scattering up to 700 MeV in energy, as well as for the astrophysical S-factor of the fusion reaction., SCOPUS: ar.j, info:eu-repo/semantics/published

The neutron-rich isotopes $^{58-63}$Cr were produced for the first time at the ISOLDE facility and their masses were measured with the ISOLTRAP spectrometer. The new values are up to 300 times more precise than those in the literature and indicate significantly different nuclear structure from the new mass-surface trend. A gradual onset of deformation is found in this proton and neutron mid-shell region, which is a gateway to the second island of inversion around \emph{N}=40. In addition to comparisons with density-functional theory and large-scale shell-model calculations, we present predictions from the valence-space formulation of the \emph{ab initio} in-medium similarity renormalization group, the first such results for open-shell chromium isotopes.

Interactions from chiral effective field theory have been successfully employed in a broad range of \textit{ab initio} calculations of nuclei and nuclear matter, but it has been observed that most results of few- and many-body calculations experience a substantial residual regulator and cutoff dependence. In this work, we investigate the behavior of local chiral potentials at leading order under variation of the cutoff scale for different local regulators. When varying the cutoff, we require that the resulting interaction produces no spurious bound states in the deuteron channel. We find that, for a particular choice of leading-order operators, nucleon-nucleon phase shifts and the deuteron ground-state energy converge to cutoff-independent plateaus, for all regulator functions we investigate. This observation may enable improved calculations with chiral Hamiltonians that also include three-nucleon interactions., 10 pages, 6 figures, published version

The standard interpretation of direct-detection limits on dark matter involves particular assumptions of the underlying WIMP-nucleus interaction, such as, in the simplest case, the choice of a Helm form factor that phenomenologically describes an isoscalar spin-independent interaction. In general, the interaction of dark matter with the target nuclei may well proceed via different mechanisms, which would lead to a different shape of the corresponding nuclear structure factors as a function of the momentum transfer $q$. We study to what extent different WIMP-nucleus responses can be differentiated based on the $q$-dependence of their structure factors (or "form factors"). We assume an overall strength of the interaction consistent with present spin-independent limits and consider an exposure corresponding to XENON1T-like, XENONnT-like, and DARWIN-like direct detection experiments. We find that, as long as the interaction strength does not lie too much below current limits, the DARWIN settings allow a conclusive discrimination of many different response functions based on their $q$-dependence, with immediate consequences for elucidating the nature of dark matter., Comment: 10 pages, 7 figures

This erratum fixes an error in the use of the low-energy constant in the two-body currents. The corrected version of one of the equations and of some of the figures is reported.

We link the structure of nuclei around $^{100}$Sn, the heaviest doubly magic nucleus with equal neutron and proton numbers ($N=Z=50$), to nucleon-nucleon ($NN$) and three-nucleon ($NNN$) forces constrained by data of few-nucleon systems. Our results indicate that $^{100}$Sn is doubly magic, and we predict its quadrupole collectivity. We present precise computations of $^{101}$Sn based on three-particle--two-hole excitations of $^{100}$Sn, and reproduce the small splitting between the lowest $J^\pi=7/2^+$ and $5/2^+$ states. Our results are consistent with the sparse available data., Comment: 8 pages, 4 figures

We present a common chiral power-counting scheme for vector, axial-vector, scalar, and pseudoscalar WIMP-nucleon interactions, and derive all one- and two-body currents up to third order in the chiral expansion. Matching our amplitudes to non-relativistic effective field theory, we find that chiral symmetry predicts a hierarchy amongst the non-relativistic operators. Moreover, we identify interaction channels where two-body currents that so far have not been accounted for become relevant., Comment: 8 pages, 1 table; journal version