Source author record

A. Schwenk

A. Schwenk appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

83works
13topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

83 published item(s)

preprint2022arXiv

Constraining Neutron-Star Matter with Microscopic and Macroscopic Collisions

Interpreting high-energy, astrophysical phenomena, such as supernova explosions or neutron-star collisions, requires a robust understanding of matter at supranuclear densities. However, our knowledge about dense matter explored in the cores of neutron stars remains limited. Fortunately, dense matter is not only probed in astrophysical observations, but also in terrestrial heavy-ion collision experiments. In this work, we use Bayesian inference to combine data from astrophysical multi-messenger observations of neutron stars and from heavy-ion collisions of gold nuclei at relativistic energies with microscopic nuclear theory calculations to improve our understanding of dense matter. We find that the inclusion of heavy-ion collision data indicates an increase in the pressure in dense matter relative to previous analyses, shifting neutron-star radii towards larger values, consistent with recent NICER observations. Our findings show that constraints from heavy-ion collision experiments show a remarkable consistency with multi-messenger observations and provide complementary information on nuclear matter at intermediate densities. This work combines nuclear theory, nuclear experiment, and astrophysical observations, and shows how joint analyses can shed light on the properties of neutron-rich supranuclear matter over the density range probed in neutron stars.

preprint2022arXiv

Importance truncation for the in-medium similarity renormalization group

Ab initio nuclear many-body frameworks require extensive computational resources, especially when targeting heavier nuclei. Importance-truncation (IT) techniques allow to significantly reduce the dimensionality of the problem by neglecting unimportant contributions to the solution of the many-body problem. In this work, we apply IT methods to the nonperturbative in-medium similarity renormalization group (IMSRG) approach and investigate the induced errors for ground-state energies in different mass regimes based on different nuclear Hamiltonians. We study various importance measures, which define the IT selection, and identify two measures that perform best, resulting in only small errors to the full IMSRG(2) calculations even for sizable compression ratios. The neglected contributions are accounted for in a perturbative way and serve as an estimate of the IT-induced error. Overall we find that the IT-IMSRG(2) performs well across all systems considered, while the largest compression ratios for a given error can be achieved when using soft Hamiltonians and for large single-particle bases.

preprint2022arXiv

Least-square approach for singular value decompositions of scattering problems

It was recently observed that chiral two-body interactions can be efficiently represented using matrix factorization techniques such as the singular value decomposition. However, the exploitation of these low-rank structures in a few- or many-body framework is nontrivial and requires reformulations that explicitly utilize the decomposition format. In this work, we present a general least-square approach that is applicable to different few- and many-body frameworks and allows for an efficient reduction to a low number of singular values in the least-square iteration. We verify the feasibility of the least-square approach by solving the Lippmann-Schwinger equation in factorized form. The resulting low-rank approximations of the $T$ matrix are found to fully capture scattering observables. Potential applications of the least-square approach to other frameworks with the goal of employing tensor factorization techniques are discussed.

preprint2021arXiv

Comparing different density-matrix expansions for long-range pion exchange

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.

preprint2021arXiv

Natural orbitals for many-body expansion methods

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.

preprint2021arXiv

New equations of state constrained by nuclear physics, observations, and QCD calculations of high-density nuclear matter

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.

preprint2021arXiv

Nuclear Charge Radii of the Nickel Isotopes $^{58-68,70}$Ni

Collinear laser spectroscopy is performed on the nickel isotopes $^{58-68,70}$Ni, using a time-resolved photon counting system. From the measured isotope shifts, nuclear charge radii $R_c$ are extracted and compared to theoretical results. Three ab initio approaches all employ, among others, the chiral interaction NNLO$_{\rm sat}$, which allows an assessment of their accuracy. We find agreement with experiment in differential radii $δ\left\langle r_\mathrm{c}^2 \right\rangle$ for all employed ab initio methods and interactions, while the absolute radii are consistent with data only for NNLO$_{\rm sat}$. Within nuclear density functional theory, the Skyrme functional SV-min matches experiment more closely than the Fayans functional Fy($Δr$,HFB).

preprint2020arXiv

$\boldsymbol{N=32}$ shell closure below calcium: Low-lying structure of $^{50}$Ar

Low-lying excited states in the $N=32$ isotope $^{50}$Ar were investigated by in-beam $γ$-ray spectroscopy following proton- and neutron-knockout, multi-nucleon removal, and proton inelastic scattering at the RIKEN Radioactive Isotope Beam Factory. The energies of the two previously reported transitions have been confirmed, and five additional states are presented for the first time, including a candidate for a 3$^-$ state. The level scheme built using $γγ$ coincidences was compared to shell-model calculations in the $sd-pf$ model space, and to ab initio predictions based on chiral two- and three-nucleon interactions. Theoretical proton- and neutron-knockout cross sections suggest that two of the new transitions correspond to $2^+$ states, while the previously proposed $4^+$ state could also correspond to a $2^+$ state.

preprint2020arXiv

Ab initio limits of atomic nuclei

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.

preprint2020arXiv

Ab initio short-range-correlation scaling factors from light to medium-mass nuclei

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.

preprint2020arXiv

Charge radius of the short-lived $^{68}$Ni and correlation with the dipole polarizability

We present the first laser spectroscopic measurement of the neutron-rich nucleus $^{68}$Ni at the \mbox{$N=40$} subshell closure and extract its nuclear charge radius. Since this is the only short-lived isotope for which the dipole polarizability $α_{\rm 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.

preprint2020arXiv

Constraining the dense matter equation of state with joint analysis of NICER and LIGO/Virgo measurements

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.

preprint2020arXiv

Dilute Fermi gas at fourth order in effective field theory

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

preprint2020arXiv

First glimpse of the $N=82$ shell closure below $Z=50$ from masses of neutron-rich cadmium isotopes and isomers

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}$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}$Sn. Using the recently implemented phase-imaging ion-cyclotron-resonance method, the ordering of the low-lying isomers in $^{129}$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.

preprint2020arXiv

Spin-polarized neutron matter, the maximum mass of neutron stars, and GW170817

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.

preprint2020arXiv

Testing ab initio nuclear structure in neutron-rich nuclei: lifetime measurements of second 2+ states in 16C and 20O

To test the predictive power of ab initio nuclear structure theory, the lifetime of the second 2+ state in neutron-rich 20O, tau(2+_2 ) = 150(+80-30) fs, and an estimate for the lifetime of the second 2+ state in 16C have been obtained, for the first time. The results were achieved via a novel Monte Carlo technique that allowed us to measure nuclear state lifetimes in the tens-to-hundreds femtoseconds range, by analyzing the Doppler-shifted gamma-transition line shapes of products of low-energy transfer and deep-inelastic processes in the reaction 18O (7.0 MeV/u) + 181Ta. The requested sensitivity could only be reached owing to the excellent performances of the AGATA gamma-tracking array, coupled to the PARIS scintillator array and to the VAMOS++ magnetic spectrometer. The experimental lifetimes agree with predictions of ab initio calculations using two- and three-nucleon interactions, obtained with the valence-space in-medium similarity renormalization group for 20O, and with the no-core shell model for 16C. The present measurement shows the power of electromagnetic observables, determined with high-precision gamma spectroscopy, to assess the quality of first-principles nuclear structure calculations, complementing common benchmarks based on nuclear energies. The proposed experimental approach will be essential for short lifetimes measurements in unexplored regions of the nuclear chart, including r-process nuclei, when intense ISOL-type beams become available.

preprint2019arXiv

Electromagnetic properties of $^{21}$O for benchmarking nuclear Hamiltonians

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 $τ_{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 $γ$-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.

preprint2016arXiv

Chiral Three-Nucleon Interactions in Light Nuclei, Neutron-$α$ Scattering, and Neutron Matter

We present quantum Monte Carlo calculations of light nuclei, neutron-$α$ scattering, and neutron matter using local two- and three-nucleon ($3N$) interactions derived from chiral effective field theory up to next-to-next-to-leading order (N$^2$LO). The two undetermined $3N$ low-energy couplings are fit to the $^4$He binding energy and, for the first time, to the spin-orbit splitting in the neutron-$α$ $P$-wave phase shifts. Furthermore, we investigate different choices of local $3N$-operator structures and find that chiral interactions at N$^2$LO are able to simultaneously reproduce the properties of $A=3,4,5$ systems and of neutron matter, in contrast to commonly used phenomenological $3N$ interactions.

preprint2016arXiv

Ground and excited states of doubly open-shell nuclei from ab initio valence-space Hamiltonians

We present ab initio predictions for ground and excited states of doubly open-shell fluorine and neon isotopes based on chiral two- and three-nucleon interactions. We use the in-medium similarity renormalization group, in both flow-equation and Magnus formulations, to derive mass-dependent sd valence-space Hamiltonians. The experimental ground-state energies are reproduced through neutron number N=14, beyond which a new targeted normal-ordering procedure improves agreement with data and large-scale multi-reference calculations. For spectroscopy, we focus on neutron-rich 23-26F and 24-26Ne isotopes near N=14,16 magic numbers. In all cases we find an agreement with experiment competitive with established phenomenology. Moreover, yrast states are well described in 20Ne and 24Mg, providing an ab initio description of deformation in the medium-mass region.

preprint2016arXiv

Quantum Monte Carlo calculations of neutron matter with chiral three-body forces

Chiral effective field theory (EFT) enables a systematic description of low-energy hadronic interactions with controlled theoretical uncertainties. For strongly interacting systems, quantum Monte Carlo (QMC) methods provide some of the most accurate solutions, but they require as input local potentials. We have recently constructed local chiral nucleon-nucleon (NN) interactions up to next-to-next-to-leading order (N$^2$LO). Chiral EFT naturally predicts consistent many-body forces. In this paper, we consider the leading chiral three-nucleon (3N) interactions in local form. These are included in auxiliary field diffusion Monte Carlo (AFDMC) simulations. We present results for the equation of state of neutron matter and for the energies and radii of neutron drops. In particular, we study the regulator dependence at the Hartree-Fock level and in AFDMC and find that present local regulators lead to less repulsion from 3N forces compared to the usual nonlocal regulators.

preprint2016arXiv

Quantum Monte Carlo calculations of two neutrons in finite volume

Ab initio calculations provide direct access to the properties of pure neutron systems that are challenging to study experimentally. In addition to their importance for fundamental physics, their properties are required as input for effective field theories of the strong interaction. In this work, we perform auxiliary-field diffusion Monte Carlo calculations of the ground and first excited state of two neutrons in a finite box, considering a simple contact potential as well as chiral effective field theory interactions. We compare the results against exact diagonalizations and present a detailed analysis of the finite-volume effects, whose understanding is crucial for determining observables from the calculated energies. Using the Lüscher formula, we extract the low-energy S-wave scattering parameters from ground- and excited-state energies for different box sizes.

preprint2016arXiv

Unexpected distribution of $\nu1f_{7/2}$ strength in the calcium isotopes at $N$=30

The calcium isotopes have emerged as an important testing ground for new microscopically derived shell-model interactions, and a great deal of focus has been directed toward this region. We investigate the relative spectroscopic strengths associated with $1f_{7/2}$ neutron hole states in $^{47, 49}$Ca following one-neutron knockout reactions from $^{48,50}$Ca. The observed reduction of strength populating the lowest 7/2$^{-}_{1}$ state in $^{49}$Ca, as compared to $^{47}$Ca, is consistent with the description given by shell-model calculations based on two- and three-nucleon forces in the neutron $pf$ model space, implying a fragmentation of the $l$=3 strength to higher-lying states. The experimental result is inconsistent with both the GXPF1 interaction routinely used in this region of the nuclear chart and with microscopic calculations in an extended model space including the $\nu1g_{9/2}$ orbital.

preprint2016arXiv

Unexpectedly large charge radii of neutron-rich calcium isotopes

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.

preprint2015arXiv

Bose-Einstein condensates in neutron stars

In the two decades since the appearance of the book "Bose-Einstein Condensation" in 1995, there have been a number of developments in our understanding of dense matter. After a brief overview of neutron star structure and the Bose-Einstein condensed phases that have been proposed, we describe selected topics, including neutron and proton pairing gaps, the physics of the inner crust of neutron stars, where a neutron fluid penetrates a lattice of nuclei, meson condensates, and pairing in dense quark matter. Especial emphasis is placed on basic physical effects and on connections to the physics of cold atomic gases.

preprint2015arXiv

Charge, neutron, and weak size of the atomic nucleus

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.

preprint2015arXiv

Charged-current reactions in the supernova neutrino-sphere

We calculate neutrino absorption rates due to charged-current reactions $ν_e+n \rightarrow e^- + p $ and $\barν_e+p \rightarrow e^+ + n $ in the outer regions of a newly born neutron star called the neutrino-sphere. To improve on recent work which has shown that nuclear mean fields enhance the $ν_e$ cross-section and suppress the $\barν_e$ cross-section, we employ realistic nucleon-nucleon interactions that fit measured scattering phase shifts. Using these interactions we calculate the momentum-, density-, and temperature-dependent nucleon self-energies in the Hartree-Fock approximation. A potential derived from chiral effective field theory and a pseudo-potential constructed to reproduce nucleon-nucleon phase shifts at the mean-field level are used to study the equilibrium proton fraction and the charged-current rates are studied in detail. We compare our results to earlier calculations obtained using phenomenological mean-field models and to those obtained in the virial expansion valid at low density. We find that for typical ambient conditions in the neutrino-sphere, $T=5-10$ MeV and $ρ=10^{11}-10^{13}$ g/cm$^3$, the difference between the $ν_{e}$ and $\barν_e$ absorption rates is much larger than predicted earlier. Our results have important implications for heavy-element nucleosynthesis in supernovae and for supernova neutrino detection.

preprint2015arXiv

Exploring sd-shell nuclei from two- and three-nucleon interactions with realistic saturation properties

We study ground- and excited-state properties of all sd-shell nuclei with neutron and proton numbers 8 <= N,Z <= 20, based on a set of low-resolution two- and three-nucleon interactions that predict realistic saturation properties of nuclear matter. We focus on estimating the theoretical uncertainties due to variation of the resolution scale, the low-energy couplings, as well as from the many-body method. The experimental two-neutron and two-proton separation energies are reasonably well reproduced, with an uncertainty range of about 5 MeV. The first excited 2+ energies also show overall agreement, with a more narrow uncertainty range of about 500 keV. In most cases, this range is dominated by the uncertainties in the Hamiltonian.

preprint2015arXiv

Ground-State Electromagnetic Moments of Calcium Isotopes

High-resolution bunched-beam collinear laser spectroscopy was used to measure the optical hyperfine spectra of the $^{43-51}$Ca isotopes. The ground state magnetic moments of $^{49,51}$Ca and quadrupole moments of $^{47,49,51}$Ca were measured for the first time, and the $^{51}$Ca ground state spin $I=3/2$ was determined in a model-independent way. Our results provide a critical test of modern nuclear theories based on shell-model calculations using phenomenological as well as microscopic interactions. The results for the neutron-rich isotopes are in excellent agreement with predictions using interactions derived from chiral effective field theory including three-nucleon forces, while lighter isotopes illustrate the presence of particle-hole excitations of the $^{40}$Ca core in their ground state.

preprint2015arXiv

Large-scale nuclear structure calculations for spin-dependent WIMP scattering with chiral effective field theory currents

We perform state-of-the-art large-scale shell-model calculations of the structure factors for elastic spin-dependent WIMP scattering off 129,131Xe, 127I, 73Ge, 19F, 23Na, 27Al, and 29Si. This comprehensive survey covers the non-zero-spin nuclei relevant to direct dark matter detection. We include a pedagogical presentation of the formalism necessary to describe elastic and inelastic WIMP-nucleus scattering. The valence spaces and nuclear interactions employed have been previously used in nuclear structure calculations for these mass regions and yield a good spectroscopic description of these isotopes. We use spin-dependent WIMP-nucleus currents based on chiral effective field theory (EFT) at the one-body level and including the leading long-range two-body currents due to pion exchange, which are predicted in chiral EFT. Results for all structure factors are provided with theoretical error bands due to the nuclear uncertainties of WIMP currents in nuclei.

preprint2015arXiv

Nuclear forces and their impact on neutron-rich nuclei and neutron-rich matter

We review the impact of nuclear forces on matter at neutron-rich extremes. Recent results have shown that neutron-rich nuclei become increasingly sensitive to three-nucleon forces, which are at the forefront of theoretical developments based on effective field theories of quantum chromodynamics. This includes the formation of shell structure, the spectroscopy of exotic nuclei, and the location of the neutron dripline. Nuclear forces also constrain the properties of neutron-rich matter, including the neutron skin, the symmetry energy, and the structure of neutron stars. We first review our understanding of three-nucleon forces and show how chiral effective field theory makes unique predictions for many-body forces. Then, we survey results with three-nucleon forces in neutron-rich oxygen and calcium isotopes and neutron-rich matter, which have been explored with a range of many-body methods. Three-nucleon forces therefore provide an exciting link between theoretical, experimental and observational nuclear physics frontiers.

preprint2015arXiv

Nuclear structure studies of $^{24}$F

The structure of the $^{24}$F nucleus has been studied at GANIL using the $β$ decay of $^{24}$O and the in-beam $γ$-ray spectroscopy from the fragmentation of projectile nuclei. Combining these complementary experimental techniques, the level scheme of $^{24}$F has been constructed up to 3.6 Mev by means of particle-$γ$ and particle-$γγ$ coincidence relations. Experimental results are compared to shell-model calculations using the standard USDA and USDB interactions as well as ab-initio valence-space Hamiltonians calculated from the in-medium similarity renormalization group based on chiral two- and three-nucleon forces. Both methods reproduce the measured level spacings well, and this close agreement allows unidentified spins and parities to be consistently assigned.

preprint2015arXiv

Pairing and superfluidity of nucleons in neutron stars

We survey the current status of understanding of pairing and superfluidity of neutrons and protons in neutron stars from a theoretical perspective, with emphasis on basic physical properties. During the past two decades, the blossoming of the field of ultracold atomic gases and the development of quantum Monte Carlo methods for solving the many-body problem have been two important sources of inspiration, and we shall describe how these have given insight into neutron pairing gaps. The equilibrium properties and collective oscillations of the inner crust of neutron stars, where neutrons paired in a $^1$S$_0$ state coexist with a lattice of neutron-rich nuclei, are also described. While pairing gaps are well understood at densities less than one tenth of the nuclear saturation density, significant uncertainties exist at higher densities due to the complicated nature of nucleon-nucleon interactions, the difficulty of solving the many-body problem under these conditions, and the increasing importance of many-nucleon interactions. We also touch more briefly on the subject of pairing of neutrons in other angular momentum states, specifically the $^3$P$_2$ state, as well as pairing of protons.

preprint2015arXiv

The In-Medium Similarity Renormalization Group: A Novel Ab Initio Method for Nuclei

We present a comprehensive review of the In-Medium Similarity Renormalization Group (IM-SRG), a novel ab inito method for nuclei. The IM-SRG employs a continuous unitary transformation of the many-body Hamiltonian to decouple the ground state from all excitations, thereby solving the many-body problem. Starting from a pedagogical introduction of the underlying concepts, the IM-SRG flow equations are developed for systems with and without explicit spherical symmetry. We study different IM-SRG generators that achieve the desired decoupling, and how they affect the details of the IM-SRG flow. Based on calculations of closed-shell nuclei, we assess possible truncations for closing the system of flow equations in practical applications, as well as choices of the reference state. We discuss the issue of center-of-mass factorization and demonstrate that the IM-SRG ground-state wave function exhibits an approximate decoupling of intrinsic and center-of-mass degrees of freedom, similar to Coupled Cluster (CC) wave functions. To put the IM-SRG in context with other many-body methods, in particular many-body perturbation theory and non-perturbative approaches like CC, a detailed perturbative analysis of the IM-SRG flow equations is carried out. We conclude with a discussion of ongoing developments, including IM-SRG calculations with three-nucleon forces, the multi-reference IM-SRG for open-shell nuclei, first non-perturbative derivations of shell- model interactions, and the consistent evolution of operators in the IM-SRG. We dedicate this review to the memory of Gerry Brown, one of the pioneers of many-body calculations of nuclei.

preprint2015arXiv

To which densities is spin-polarized neutron matter a weakly interacting Fermi gas?

We study the properties of spin-polarized neutron matter at next-to-next-to-next-to-leading order in chiral effective field theory, including two-, three-, and four-neutron interactions. The energy of spin-polarized neutrons is remarkably close to a non-interacting system at least up to saturation density, where interaction effects provide less than 10% corrections. This shows that the physics of neutron matter is similar to a unitary gas well beyond the scattering-length regime. Implications for energy-density functionals and for a possible ferromagnetic transition in neutron stars are discussed. Our predictions can be tested with lattice QCD, and we present results for varying pion mass.

preprint2014arXiv

Breakdown of the Isobaric Multiplet Mass Equation for the A = 20 and 21 Multiplets

Using the Penning trap mass spectrometer TITAN, we performed the first direct mass measurements of 20,21Mg, isotopes that are the most proton-rich members of the A = 20 and A = 21 isospin multiplets. These measurements were possible through the use of a unique ion-guide laser ion source, a development that suppressed isobaric contamination by six orders of magnitude. Compared to the latest atomic mass evaluation, we find that the mass of 21Mg is in good agreement but that the mass of 20Mg deviates by 3σ. These measurements reduce the uncertainties in the masses of 20,21Mg by 15 and 22 times, respectively, resulting in a significant departure from the expected behavior of the isobaric multiplet mass equation in both the A = 20 and A = 21 multiplets. This presents a challenge to shell model calculations using either the isospin non-conserving USDA/B Hamiltonians or isospin non-conserving interactions based on chiral two- and three-nucleon forces.

preprint2014arXiv

Local chiral effective field theory interactions and quantum Monte Carlo applications

We present details of the derivation of local chiral effective field theory interactions to next-to-next-to-leading order, and show results for nucleon-nucleon phase shifts and deuteron properties for these potentials. We then perform systematic auxiliary-field diffusion Monte Carlo calculations for neutron matter based on the developed local chiral potentials at different orders. This includes studies of the effects of the spectral-function regularization and of the local regulators. For all orders, we compare the quantum Monte Carlo results with perturbative many-body calculations and find excellent agreement for low cutoffs.

preprint2014arXiv

Microscopic calculations and energy expansions for neutron-rich matter

We investigate asymmetric nuclear matter with two- and three-nucleon interactions based on chiral effective field theory, where three-body forces are fit only to light nuclei. Focusing on neutron-rich matter, we calculate the energy for different proton fractions and include estimates of the theoretical uncertainty. We use our ab-initio results to test the quadratic expansion around symmetric matter with the symmetry energy term, and confirm its validity for highly asymmetric systems. Our calculations are in remarkable agreement with an empirical parametrization for the energy density. These findings are very useful for astrophysical applications and for developing new equations of state.

preprint2014arXiv

Nonperturbative shell-model interactions from the in-medium similarity renormalization group

We present the first ab initio construction of valence-space Hamiltonians for medium-mass nuclei based on chiral two- and three-nucleon interactions using the in-medium similarity renormalization group. When applied to the oxygen isotopes, we find experimental ground-state energies are well reproduced, including the flat trend beyond the drip line at 24O. Similarly, natural-parity spectra in 21,22,23,24O are in agreement with experiment, and we present predictions for excited states in 25,26O. The results exhibit a weak dependence on the harmonic-oscillator basis parameter and reproduce spectroscopy within the standard sd valence space.

preprint2014arXiv

Quantum Monte Carlo Calculations of Light Nuclei Using Chiral Potentials

We present the first Green's function Monte Carlo calculations of light nuclei with nuclear interactions derived from chiral effective field theory up to next-to-next-to-leading order. Up to this order, the interactions can be constructed in a local form and are therefore amenable to quantum Monte Carlo calculations. We demonstrate a systematic improvement with each order for the binding energies of $A=3$ and $A=4$ systems. We also carry out the first few-body tests to study perturbative expansions of chiral potentials at different orders, finding that higher-order corrections are more perturbative for softer interactions. Our results confirm the necessity of a three-body force for correct reproduction of experimental binding energies and radii, and pave the way for studying few- and many-nucleon systems using quantum Monte Carlo methods with chiral interactions.

preprint2014arXiv

Supernova matter at subnuclear densities as a resonant Fermi gas: Enhancement of neutrino rates

At low energies nucleon-nucleon interactions are resonant and therefore supernova matter at subnuclear densities has many similarities to atomic gases with interactions dominated by a Feshbach resonance. We calculate the rates of neutrino processes involving nucleon-nucleon collisions and show that these are enhanced in mixtures of neutrons and protons at subnuclear densities due to the large scattering lengths. As a result, the rate for neutrino pair bremsstrahlung and absorption is significantly larger below 10^{13} g cm^{-3} compared to rates used in supernova simulations.

preprint2014arXiv

The Large Observatory For x-ray Timing

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

preprint2014arXiv

The neutron polaron as a constraint on nuclear density functionals

We study the energy of an impurity (polaron) that interacts strongly in a sea of fermions when the effective range of the impurity-fermion interaction becomes important, thereby mapping the Fermi polaron of condensed matter physics and ultracold atoms to strongly interacting neutrons. We present Quantum Monte Carlo results for this neutron polaron, and compare these with effective field theory calculations that also include contributions beyond the effective range. We show that state-of-the-art nuclear density functionals vary substantially and generally underestimate the neutron polaron energy. Our results thus provide constraints for adjusting the time-odd components of nuclear density functionals to better characterize polarized systems.

preprint2014arXiv

Three-nucleon forces and spectroscopy of neutron-rich calcium isotopes

We study excited-state properties of neutron-rich calcium isotopes based on chiral two- and three-nucleon interactions. We first discuss the details of our many-body framework, investigate convergence properties, and for two-nucleon interactions benchmark against coupled-cluster calculations. We then focus on the spectroscopy of 47-56Ca, finding that with both 3N forces and an extended pfg9/2 valence space, we obtain a good level of agreement with experiment. We also study electromagnetic transitions and find that experimental data are well described by our calculations. In addition, we provide predictions for unexplored properties of neutron-rich calcium isotopes.

preprint2013arXiv

Chiral three-nucleon forces and bound excited states in neutron-rich oxygen isotopes

We study the spectra of neutron-rich oxygen isotopes based on chiral two- and three-nucleon interactions. First, we benchmark our many-body approach by comparing ground-state energies to coupled-cluster results for the same two-nucleon interaction, with overall good agreement. We then calculate bound excited states in 21,22,23O, focusing on the role of three-nucleon forces, in the standard sd shell and an extended sdf7/2p3/2 valence space. Chiral three-nucleon forces provide important one- and two-body contributions between valence neutrons. We find that both these contributions and an extended valence space are necessary to reproduce key signatures of novel shell evolution, such as the N = 14 magic number and the low-lying states in 21O and 23O, which are too compressed with two-nucleon interactions only. For the extended space calculations, this presents first work based on nuclear forces without adjustments. Future work is needed and open questions are discussed.

preprint2013arXiv

Constraints on Skyrme Equations of State from Properties of Doubly Magic Nuclei and Ab-Initio Calculations of Low-Density Neutron Matter

We use properties of doubly-magic nuclei and ab-initio calculations of low-density neutron matter to constrain Skyrme equations of state for neutron-rich conditions. All of these properties are consistent with a Skyrme functional form and a neutron-matter equation of state that depends on three parameters. With a reasonable range for the neutron-matter effective mass, the values of the two other Skyrme parameters are well constrained. This leads to predictions for other quantities. The neutron skins for $^{208}$Pb and $^{48}$Ca are predicted to be 0.182(10) fm and 0.173(5) fm, respectively. Other results including the dipole polarizability are discussed.

preprint2013arXiv

Equation of state and neutron star properties constrained by nuclear physics and observation

Microscopic calculations of neutron matter based on nuclear interactions derived from chiral effective field theory, combined with the recent observation of a 1.97 +- 0.04 M_sun neutron star, constrain the equation of state of neutron-rich matter at sub- and supranuclear densities. We discuss in detail the allowed equations of state and the impact of our results on the structure of neutron stars, the crust-core transition density, and the nuclear symmetry energy. In particular, we show that the predicted range for neutron star radii is robust. For use in astrophysical simulations, we provide detailed numerical tables for a representative set of equations of state consistent with these constraints.

preprint2013arXiv

Neutron matter at next-to-next-to-next-to-leading order in chiral effective field theory

Neutron matter presents a unique system for chiral effective field theory (EFT), because all many-body forces among neutrons are predicted to next-to-next-to-next-to-leading order (N3LO). We present the first complete N3LO calculation of the neutron matter energy. This includes the subleading three-nucleon (3N) forces for the first time and all leading four-nucleon (4N) forces. We find relatively large contributions from N3LO 3N forces. Our results provide constraints for neutron-rich matter in astrophysics with controlled theoretical uncertainties.

preprint2013arXiv

Neutron matter from chiral effective field theory interactions

The neutron-matter equation of state constrains the properties of many physical systems over a wide density range and can be studied systematically using chiral effective field theory (EFT). In chiral EFT, all many-body forces among neutrons are predicted to next-to-next-to-next-to-leading order (N3LO). We present details and additional results of the first complete N3LO calculation of the neutron-matter energy, which includes the subleading three-nucleon as well as the leading four-nucleon forces, and provides theoretical uncertainties. In addition, we discuss the impact of our results for astrophysics: for the supernova equation of state, the symmetry energy and its density derivative, and for the structure of neutron stars. Finally, we give a first estimate for the size of the N3LO many-body contributions to the energy of symmetric nuclear matter, which shows that their inclusion will be important in nuclear structure calculations.

preprint2013arXiv

Neutron matter with chiral EFT interactions: Perturbative and first QMC calculations

Neutron matter presents a unique system in chiral effective field theory (EFT), because all many-body forces among neutrons are predicted to next-to-next-to-next-to-leading order (N3LO). We discuss perturbative and first Quantum Monte Carlo (QMC) calculations of neutron matter with chiral EFT interactions and their astrophysical impact for the equation of state and neutron stars.

preprint2013arXiv

Nuclear astrophysics and electron beams

Electron beams provide important probes and constraints for nuclear astrophysics. This is especially exciting at energies within the regime of chiral effective field theory (EFT), which provides a systematic expansion for nuclear forces and electroweak operators based on quantum chromodynamics. This talk discusses some recent highlights and future directions based on chiral EFT, including nuclear structure and reactions for astrophysics, the neutron skin and constraints for the properties of neutron-rich matter in neutron stars and core-collapse supernovae, and the dark matter response of nuclei.

preprint2013arXiv

Quantum Monte Carlo Calculations with Chiral Effective Field Theory Interactions

We present the first quantum Monte Carlo (QMC) calculations with chiral effective field theory (EFT) interactions. To achieve this, we remove all sources of nonlocality, which hamper the inclusion in QMC calculations, in nuclear forces to next-to-next-to-leading order. We perform auxiliary-field diffusion Monte Carlo (AFDMC) calculations for the neutron matter energy up to saturation density based on local leading-order, next-to-leading order, and next-to-next-to-leading order nucleon-nucleon interactions. Our results exhibit a systematic order-by-order convergence in chiral EFT and provide nonperturbative benchmarks with theoretical uncertainties. For the softer interactions, perturbative calculations are in excellent agreement with the AFDMC results. This work paves the way for QMC calculations with systematic chiral EFT interactions for nuclei and nuclear matter, for testing the perturbativeness of different orders, and allows for matching to lattice QCD results by varying the pion mass.

preprint2013arXiv

Role of three-nucleon forces and many-body processes in nuclear pairing

We present microscopic valence-shell calculations of pairing gaps in the calcium isotopes, focusing on the role of three-nucleon (3N) forces and many-body processes. In most cases, we find a reduction in pairing strength when the leading chiral 3N forces are included, compared to results with low-momentum two-nucleon (NN) interactions only. This is in agreement with a recent energy density functional study. At the NN level, calculations that include particle-particle and hole-hole ladder contributions lead to smaller pairing gaps compared with experiment. When particle-hole contributions as well as the normal-ordered one- and two-body parts of 3N forces are consistently included to third order, we find reasonable agreement with experimental three-point mass differences. This highlights the important role of 3N forces and many-body processes for pairing in nuclei. Finally, we relate pairing gaps to the evolution of nuclear structure in neutron-rich calcium isotopes and study the predictions for the 2+ excitation energies, in particular for 54Ca.

preprint2013arXiv

Signatures of Dark Matter Scattering Inelastically Off Nuclei

Direct dark matter detection focuses on elastic scattering of dark matter particles off nuclei. In this study, we explore inelastic scattering where the nucleus is excited to a low-lying state of 10-100 keV, with subsequent prompt de-excitation. We calculate the inelastic structure factors for the odd-mass xenon isotopes based on state-of-the-art large-scale shell-model calculations with chiral effective field theory WIMP-nucleon currents. For these cases, we find that the inelastic channel is comparable to or can dominate the elastic channel for momentum transfers around 150 MeV. We calculate the inelastic recoil spectra in the standard halo model, compare these to the elastic case, and discuss the expected signatures in a xenon detector, along with implications for existing and future experiments. The combined information from elastic and inelastic scattering will allow to determine the dominant interaction channel within one experiment. In addition, the two channels probe different regions of the dark matter velocity distribution and can provide insight into the dark halo structure. The allowed recoil energy domain and the recoil energy at which the integrated inelastic rates start to dominate the elastic channel depend on the mass of the dark matter particle, thus providing a potential handle to constrain its mass.

preprint2013arXiv

The chiral condensate in neutron matter

We calculate the chiral condensate in neutron matter at zero temperature based on nuclear forces derived within chiral effective field theory. Two-, three- and four-nucleon interactions are included consistently to next-to-next-to-next-to-leading order (N3LO) of the chiral expansion. We find that the interaction contributions lead to a modest increase of the condensate, thus impeding the restoration of chiral symmetry in dense matter and making a chiral phase transition in neutron-rich matter unlikely for densities that are not significantly higher than nuclear saturation density.

preprint2013arXiv

Three-body forces and proton-rich nuclei

We present the first study of three-nucleon (3N) forces for proton-rich nuclei along the N=8 and N=20 isotones. Our results for the ground-state energies and proton separation energies are in very good agreement with experiment where available, and with the empirical isobaric multiplet mass equation. We predict the spectra for all N=8 and N=20 isotones to the proton dripline, which agree well with experiment for 18Ne, 19Na, 20Mg and 42Ti. In all other cases, we provide first predictions based on nuclear forces. Our results are also very promising for studying isospin symmetry breaking in medium-mass nuclei based on chiral effective field theory.

preprint2013arXiv

Three-body forces: From cold atoms to nuclei

It is often assumed that few- and many-body systems can be accurately described by considering only pairwise two-body interactions of the constituents. We illustrate that three- and higher-body forces enter naturally in effective field theories and are especially prominent in strongly interacting quantum systems. We focus on three-body forces and discuss examples from atomic and nuclear physics. In particular, we highlight the importance and the challenges of three-nucleon forces for nuclear structure and reactions, including applications to astrophysics and fundamental symmetries.

preprint2013arXiv

Three-nucleon forces and nuclei at the extremes

Neutron-rich nuclei become increasingly sensitive to three-nucleon forces. These components of nuclear forces are at the forefront of theoretical developments based on effective field theories of quantum chromodynamics. We discuss our understanding of three-nucleon forces and their impact on exotic nuclei, and show how new measurements test and constrain them. Three-nucleon forces therefore provide an exciting link between theoretical and experimental nuclear physics frontiers.

preprint2012arXiv

Equation-of-state dependence of the gravitational-wave signal from the ring-down phase of neutron-star mergers

Neutron-star (NS) merger simulations are conducted for 38 representative microphysical descriptions of high-density matter in order to explore the equation-of-state dependence of the postmerger ring-down phase. The formation of a deformed, oscillating, differentially rotating very massive NS is the typical outcome of the coalescence of two stars with 1.35 $M_{\odot}$ for most candidate EoSs. The oscillations of this object imprint a pronounced peak in the gravitational-wave (GW) spectra, which is used to characterize the emission for a given model. The peak frequency of this postmerger GW signal correlates very well with the radii of nonrotating NSs, and thus allows to constrain the high-density EoS by a GW detection. In the case of 1.35-1.35 $M_{\odot}$ mergers the peak frequency scales particularly well with the radius of a NS with 1.6 $M_{\odot}$, where the maximum deviation from this correlation is only 60 meters for fully microphysical EoSs which are compatible with NS observations. Combined with the uncertainty in the determination of the peak frequency it appears likely that a GW detection can measure the radius of a 1.6 $M_{\odot}$ NS with an accuracy of about 100 to 200 meters. We also uncover relations of the peak frequency with the radii of nonrotating NSs with 1.35 $M_{\odot}$ or 1.8 $M_{\odot}$, with the radius or the central energy density of the maximum-mass Tolman-Oppenheimer-Volkoff configuration, and with the pressure or sound speed at a fiducial rest-mass density of about twice nuclear saturation density. Furthermore, it is found that a determination of the dominant postmerger GW frequency can provide an upper limit for the maximum mass of nonrotating NSs. The prospects for a detection of the postmerger GW signal and a determination of the dominant GW frequency are estimated to be in the range of 0.015 to 1.2 events per year with the upcoming Advanced LIGO detector.

preprint2012arXiv

In-Medium Similarity Renormalization Group for Open-Shell Nuclei

We present a nonperturbative derivation of effective valence-shell Hamiltonians in the framework of the recently developed in-medium similarity renormalization group (IM-SRG). As a first application, we calculate the spectra of $p$- and $sd$-shell nuclei, $^6$Li and $^{18}$O, based on evolved chiral nucleon-nucleon interactions. For $^6$Li, the spectrum is in very good agreement with ab-initio results. For $^{18}$O, the IM-SRG provides a new method for the shell model to systematically go beyond effective interaction techniques based on diagrammatic expansions.

preprint2012arXiv

In-Medium Similarity Renormalization Group with Chiral Two- Plus Three-Nucleon Interactions

We use the recently proposed In-Medium Similarity Renormalization Group (IM-SRG) to carry out a systematic study of closed-shell nuclei up to $\nuc{Ni}{56}$, based on chiral two- plus three-nucleon interactions. We analyze the capabilities of the IM-SRG by comparing our results for the ground-state energy to Coupled Cluster calculations, as well as to quasi-exact results from the Importance-Truncated No-Core Shell Model. Using chiral two- plus three-nucleon Hamiltonians whose resolution scales are lowered by free-space SRG evolution, we obtain good agreement with experimental binding energies in $\nuc{He}{4}$ and the closed-shell oxygen isotopes, while the calcium and nickel isotopes are somewhat overbound.

preprint2012arXiv

LOFT: the Large Observatory For X-ray Timing

The LOFT mission concept is one of four candidates selected by ESA for the M3 launch opportunity as Medium Size missions of the Cosmic Vision programme. The launch window is currently planned for between 2022 and 2024. LOFT is designed to exploit the diagnostics of rapid X-ray flux and spectral variability that directly probe the motion of matter down to distances very close to black holes and neutron stars, as well as the physical state of ultra-dense matter. These primary science goals will be addressed by a payload composed of a Large Area Detector (LAD) and a Wide Field Monitor (WFM). The LAD is a collimated (<1 degree field of view) experiment operating in the energy range 2-50 keV, with a 10 m^2 peak effective area and an energy resolution of 260 eV at 6 keV. The WFM will operate in the same energy range as the LAD, enabling simultaneous monitoring of a few-steradian wide field of view, with an angular resolution of <5 arcmin. The LAD and WFM experiments will allow us to investigate variability from submillisecond QPO's to year-long transient outbursts. In this paper we report the current status of the project.

preprint2012arXiv

Matter and charge radius of 6He in the hyperspherical-harmonics approach

We present ab-initio calculations of the binding energy and radii of the two-neutron halo nucleus 6He using two-body low-momentum interactions based on chiral effective field theory potentials. Calculations are performed via a hyperspherical harmonics expansion where the convergence is sped up introducing an effective interaction for non-local potentials. The latter is essential to reach a satisfactory convergence of the extended matter radius and of the point-proton radius. The dependence of the results on the resolution scale is studied. A correlation is found between the radii and the two-neutron separation energy. The importance of three-nucleon forces is pointed out comparing our results and previous calculations to experiment.

preprint2012arXiv

Neutrino processes in partially degenerate neutron matter

We investigate neutrino processes for conditions reached in simulations of core-collapse supernovae. Where neutrino-matter interactions play an important role, matter is partially degenerate, and we extend earlier work that addressed the degenerate regime. We derive expressions for the spin structure factor in neutron matter, which is a key quantity required for evaluating rates of neutrino processes. We show that, for essentially all conditions encountered in the post-bounce phase of core-collapse supernovae, it is a very good approximation to calculate the spin relaxation rates in the nondegenerate limit. We calculate spin relaxation rates based on chiral effective field theory interactions and find that they are typically a factor of two smaller than those obtained using the standard one-pion-exchange interaction alone.

preprint2012arXiv

New precision mass measurements of neutron-rich calcium and potassium isotopes and three-nucleon forces

We present precision Penning-trap mass measurements of neutron-rich calcium and potassium isotopes in the vicinity of neutron number N=32. Using the TITAN system the mass of $^{51}$K was measured for the first time, and the precision of the $^{51,52}$Ca mass values were improved significantly. The new mass values show a dramatic increase of the binding energy compared to those reported in the atomic mass evaluation. In particular, $^{52}$Ca is more bound by 1.74 MeV, and the behavior with neutron number deviates substantially from the tabulated values. An increased binding was predicted recently based on calculations that include three-nucleon (3N) forces. We present a comparison to improved calculations, which agree remarkably with the evolution of masses with neutron number, making neutron-rich calcium isotopes an exciting region to probe 3N forces at neutron-rich extremes.

preprint2012arXiv

Spin-dependent WIMP scattering off nuclei

Chiral effective field theory (EFT) provides a systematic expansion for the coupling of WIMPs to nucleons at the momentum transfers relevant to direct cold dark matter detection. We derive the currents for spin-dependent WIMP scattering off nuclei at the one-body level and include the leading long-range two-body currents, which are predicted in chiral EFT. As an application, we calculate the structure factor for spin-dependent WIMP scattering off 129,131Xe nuclei, using nuclear interactions that have been developed to study nuclear structure and double-beta decays in this region. We provide theoretical error bands due to the nuclear uncertainties of WIMP currents in nuclei.

preprint2011arXiv

Chiral two-body currents in nuclei: Gamow-Teller transitions and neutrinoless double-beta decay

We show that chiral effective field theory (EFT) two-body currents provide important contributions to the quenching of low-momentum-transfer Gamow-Teller transitions, and use chiral EFT to predict the momentum-transfer dependence that is probed in neutrinoless double-beta decay. We then calculate for the first time the neutrinoless double-beta decay operator based on chiral EFT currents and study the nuclear matrix elements at successive orders. The contributions from chiral two-body currents are significant and should be included in all calculations.

preprint2011arXiv

First direct mass-measurement of the two-neutron halo nucleus 6He and improved mass for the four-neutron halo 8He

The first direct mass-measurement of $^{6}$He has been performed with the TITAN Penning trap mass spectrometer at the ISAC facility. In addition, the mass of $^{8}$He was determined with improved precision over our previous measurement. The obtained masses are $m$($^{6}$He) = 6.018 885 883(57) u and $m$($^{8}$He) = 8.033 934 44(11) u. The $^{6}$He value shows a deviation from the literature of 4$σ$. With these new mass values and the previously measured atomic isotope shifts we obtain charge radii of 2.060(8) fm and 1.959(16) fm for $^{6}$He and $^{8}$He respectively. We present a detailed comparison to nuclear theory for $^6$He, including new hyperspherical harmonics results. A correlation plot of the point-proton radius with the two-neutron separation energy demonstrates clearly the importance of three-nucleon forces.

preprint2011arXiv

Improved nuclear matter calculations from chiral low-momentum interactions

We present new nuclear matter calculations based on low-momentum interactions derived from chiral effective field theory potentials. The current calculations use an improved treatment of the three-nucleon force contribution that includes a corrected combinatorial factor beyond Hartree-Fock that was omitted in previous nuclear matter calculations. We find realistic saturation properties using parameters fit only to few-body data, but with larger uncertainty estimates from cutoff dependence and the 3NF parametrization than in previous calculations.

preprint2011arXiv

In-Medium Similarity Renormalization Group for Nuclei

We present a new ab-initio method that uses similarity renormalization group (SRG) techniques to continuously diagonalize nuclear many-body Hamiltonians. In contrast with applications of the SRG to two- and three-nucleon interactions in free space, we perform the SRG evolution "in medium" directly in the $A$-body system of interest. The in-medium approach has the advantage that one can approximately evolve $3,...,A$-body operators using only two-body machinery based on normal-ordering techniques. The method is nonperturbative and can be tailored to problems ranging from the diagonalization of closed-shell nuclei to the construction of effective valence shell-model Hamiltonians and operators. We present first results for the ground-state energies of $^4$He, $^{16}$O and $^{40}$Ca, which have accuracies comparable to coupled-cluster calculations.

preprint2011arXiv

Three-body interactions in Fermi systems

We show that the contributions of three-quasiparticle interactions to normal Fermi systems at low energies and temperatures are suppressed by n_q/n compared to two-body interactions, where n_q is the density of excited or added quasiparticles and n is the ground-state density. For finite Fermi systems, three-quasiparticle contributions are suppressed by the corresponding ratio of particle numbers N_q/N. This is illustrated for polarons in strongly interacting spin-polarized Fermi gases and for valence neutrons in neutron-rich calcium isotopes.

preprint2010arXiv

Chiral three-nucleon forces and neutron matter

We calculate the properties of neutron matter and highlight the physics of chiral three-nucleon forces. For neutrons, only the long-range 2 pi-exchange interactions of the leading chiral three-nucleon forces contribute, and we derive density-dependent two-body interactions by summing the third particle over occupied states in the Fermi sea. Our results for the energy suggest that neutron matter is perturbative at nuclear densities. We study in detail the theoretical uncertainties of the neutron matter energy, provide constraints for the symmetry energy and its density dependence, and explore the impact of chiral three-nucleon forces on the S-wave superfluid pairing gap.

preprint2010arXiv

Constraints on neutron star radii based on chiral effective field theory interactions

We show that microscopic calculations based on chiral effective field theory interactions constrain the properties of neutron-rich matter below nuclear densities to a much higher degree than is reflected in commonly used equations of state. Combined with observed neutron star masses, our results lead to a radius R = 9.7 - 13.9 km for a 1.4 M_{solar} star, where the theoretical range is due, in about equal amounts, to uncertainties in many-body forces and to the extrapolation to high densities.

preprint2010arXiv

From low-momentum interactions to nuclear structure

We present an overview of low-momentum two-nucleon and many-body interactions and their use in calculations of nuclei and infinite matter. The softening of phenomenological and effective field theory (EFT) potentials by renormalization group (RG) transformations that decouple low and high momenta leads to greatly enhanced convergence in few- and many-body systems while maintaining a decreasing hierarchy of many-body forces. This review surveys the RG-based technology and results, discusses the connections to chiral EFT, and clarifies various misconceptions.

preprint2010arXiv

How should one formulate, extract, and interpret `non-observables' for nuclei?

Nuclear observables such as binding energies and cross sections can be directly measured. Other physically useful quantities, such as spectroscopic factors, are related to measured quantities by a convolution whose decomposition is not unique. Can a framework for these nuclear structure `non-observables' be formulated systematically so that they can be extracted from experiment with known uncertainties and calculated with consistent theory? Parton distribution functions in hadrons serve as an illustrative example of how this can be done. A systematic framework is also needed to address questions of interpretation, such as whether short-range correlations are important for nuclear structure.

preprint2010arXiv

Lowest-order contributions of chiral three-nucleon interactions to pairing properties of nuclear ground states

We perform a systematic study of the odd-even mass staggering generated using a pairing interaction computed at first order in low-momentum interactions. Building on previous work including the (nuclear plus Coulomb) two-nucleon interaction only, we focus here on the first-order contribution from chiral three-nucleon forces. We observe a significant repulsive effect from the three-nucleon interaction. The combined contribution from two- and three-nucleon interactions accounts for approximately 70% of the experimental pairing gaps. This leaves room for higher-order contributions to the pairing kernel and the normal self-energy that need to be computed consistently.

preprint2009arXiv

Isospin-symmetry-breaking corrections to superallowed Fermi beta decay: Radial excitations

Based on an exact formalism, we study the effects of isospin-symmetry-breaking interactions on superallowed 0+ -> 0+ transitions. We calculate the second-order renormalization of the Fermi matrix element due to radial contributions and show that radial excitations neglected in the treatment of Towner and Hardy are significant. These are estimated to decrease the isospin-symmetry-breaking corrections. Our results provide a correction term that can be included in existing approaches.

preprint2009arXiv

Nuclei and chiral dynamics

Nuclear theory has entered an exciting era. This is due to advances on many fronts, including the development of effective field theory and the renormalization group for nuclear forces, advances in ab-initio methods for nuclear structure, an effort to develop a universal density functional based on microscopic interactions, and the application of large-scale computing resources. I discuss their impact, recent highlights and the frontiers in understanding and predicting nuclei and the structure of strongly-interacting matter based on chiral interactions.

preprint2009arXiv

Spin response of a normal Fermi liquid with noncentral interactions

We consider the spin response of a normal Fermi liquid with noncentral interactions under conditions intermediate between the collisionless and hydrodynamic regimes. This problem is of importance for calculations of neutrino properties in dense matter. By expressing the deviation of the quasiparticle distribution function from equilibrium in terms of eigenfunctions of the transport equation under the combined influence of collisions and an external field, we derive a closed expression for the spin-density--spin-density response function and compare its predictions with that of a relaxation time approximation. Our results indicate that the relaxation time approximation is reliable for neutrino properties under astrophysically relevant conditions.

preprint2008arXiv

Three-nucleon interactions: A frontier in nuclear structure

Three-nucleon interactions are a frontier in understanding and predicting the structure of strongly-interacting matter in laboratory nuclei and in the cosmos. We present results and discuss the status of first calculations with microscopic three-nucleon interactions beyond light nuclei. This coherent effort is possible due to advances based on effective field theory and renormalization group methods in nuclear physics.