Showing 1–50 of 77 results for author: Schwenk, A

Neutron matter from local chiral EFT interactions at large cutoffs

Authors: I. Tews, R. Somasundaram, D. Lonardoni, H. Göttling, R. Seutin, J. Carlson, S. Gandolfi, K. Hebeler, A. Schwenk

Abstract: Neutron matter is an important many-body system that provides valuable constraints for the equation of state (EOS) of neutron stars. Neutron-matter calculations employing chiral effective field theory (EFT) interactions have been extensively used for this purpose. Among the various many-body methods, quantum Monte Carlo (QMC) methods stand out due to their nonperturbative nature and the achievable… ▽ More Neutron matter is an important many-body system that provides valuable constraints for the equation of state (EOS) of neutron stars. Neutron-matter calculations employing chiral effective field theory (EFT) interactions have been extensively used for this purpose. Among the various many-body methods, quantum Monte Carlo (QMC) methods stand out due to their nonperturbative nature and the achievable precision. However, QMC methods require local interactions as input, which leads to the appearance of stronger regulator artifacts as compared to non-local interactions. To circumvent this, we employ large-cutoff interactions derived within chiral EFT (400 MeV $\leq Λ_c \leq$ 700 MeV) for studies of pure neutron matter. These interactions have been adjusted to nucleon-nucleon scattering phase shifts, the triton binding energy, as well as the triton beta-decay half life. We find that regulator artifacts significantly decrease with increasing cutoff, leading to a significant reduction of uncertainties in the neutron-matter EOS. We discuss implications for the symmetry energy and demonstrate how our new calculations lead to a reduction in the theoretical uncertainty of predicted neutron-star radii by up to 30% for low-mass stars. △ Less

Submitted 12 July, 2024; originally announced July 2024.

Comments: 7 pages, 3 figures

Report number: LA-UR-23-32163

Constraining the dense matter equation of state with new NICER mass-radius measurements and new chiral effective field theory inputs

Authors: Nathan Rutherford, Melissa Mendes, Isak Svensson, Achim Schwenk, Anna L. Watts, Kai Hebeler, Jonas Keller, Chanda Prescod-Weinstein, Devarshi Choudhury, Geert Raaijmakers, Tuomo Salmi, Patrick Timmerman, Serena Vinciguerra, Sebastien Guillot, James M. Lattimer

Abstract: Pulse profile modeling of X-ray data from NICER is now enabling precision inference of neutron star mass and radius. Combined with nuclear physics constraints from chiral effective field theory ($χ$EFT), and masses and tidal deformabilities inferred from gravitational wave detections of binary neutron star mergers, this has lead to a steady improvement in our understanding of the dense matter equa… ▽ More Pulse profile modeling of X-ray data from NICER is now enabling precision inference of neutron star mass and radius. Combined with nuclear physics constraints from chiral effective field theory ($χ$EFT), and masses and tidal deformabilities inferred from gravitational wave detections of binary neutron star mergers, this has lead to a steady improvement in our understanding of the dense matter equation of state (EOS). Here we consider the impact of several new results: the radius measurement for the 1.42$\,M_\odot$ pulsar PSR J0437$-$4715 presented by Choudhury et al. (2024), updates to the masses and radii of PSR J0740$+$6620 and PSR J0030$+$0451, and new $χ$EFT results for neutron star matter up to 1.5 times nuclear saturation density. Using two different high-density EOS extensions -- a piecewise-polytropic (PP) model and a model based on the speed of sound in a neutron star (CS) -- we find the radius of a 1.4$\,M_\odot$ (2.0$\,M_\odot$) neutron star to be constrained to the 95\% credible ranges $12.28^{+0.50}_{-0.76}\,$km ($12.33^{+0.70}_{-1.34}\,$km) for the PP model and $12.01^{+0.56}_{-0.75}\,$km ($11.55^{+0.94}_{-1.09}\,$km) for the CS model. The maximum neutron star mass is predicted to be $2.15^{+0.14}_{-0.16}\,$$M_\odot$ and $2.08^{+0.28}_{-0.16}\,$$M_\odot$ for the PP and CS model, respectively. We explore the sensitivity of our results to different orders and different densities up to which $χ$EFT is used, and show how the astrophysical observations provide constraints for the pressure at intermediate densities. Moreover, we investigate the difference $R_{2.0} - R_{1.4}$ of the radius of 2$\,M_\odot$ and 1.4$\,M_\odot$ neutron stars within our EOS inference. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: 23 pages, 15 figures; accepted for publication in ApJ Letters. This paper was made to be as similar to the ApJ Letters version as possible

Neutron star matter as a dilute solution of protons in neutrons

Authors: J. Keller, K. Hebeler, C. J. Pethick, A. Schwenk

Abstract: Neutron stars contain neutron-rich matter with around 5% protons at nuclear saturation density. In this Letter, we consider equilibrium between bulk phases of matter based on asymmetric nuclear matter calculations using chiral effective field theory interactions rather than, as has been done in the past, by interpolation between the properties of symmetric nuclear matter and pure neutron matter. N… ▽ More Neutron stars contain neutron-rich matter with around 5% protons at nuclear saturation density. In this Letter, we consider equilibrium between bulk phases of matter based on asymmetric nuclear matter calculations using chiral effective field theory interactions rather than, as has been done in the past, by interpolation between the properties of symmetric nuclear matter and pure neutron matter. Neutron drip (coexistence of nuclear matter with pure neutrons) is well established, but from earlier work it is unclear whether proton drip (equilibrium between two phases, both of which contain protons and neutrons) is possible. We find that proton drip is a robust prediction of any physically reasonable equation of state, but that it occurs over a limited region of densities and proton fractions. An analytical model based on expanding the energy in powers of the proton density, rather than the neutron excess, is able to account for these features of the phase diagram. △ Less

Submitted 10 June, 2024; v1 submitted 24 January, 2024; originally announced January 2024.

Comments: 6 pages, 4 figures, published version

Report number: NORDITA 2024-006

Journal ref: Phys. Rev. Lett. 132, 232701 (2024)

Effects of nuclear matter properties in neutron star mergers

Authors: Maximilian Jacobi, Federico Maria Guercilena, Sabrina Huth, Giacomo Ricigliano, Almudena Arcones, Achim Schwenk

Abstract: The dynamics in mergers of binary neutron star (BNS) systems depend sensitively on the equation of state (EOS) of dense matter. This has profound implications on the emission of gravitational waves (GWs) and the ejection of matter in the merger and post-merger phases and is thus of high interest for multi-messenger astronomy. Today, a variety of nuclear EOSs are available with various underlying m… ▽ More The dynamics in mergers of binary neutron star (BNS) systems depend sensitively on the equation of state (EOS) of dense matter. This has profound implications on the emission of gravitational waves (GWs) and the ejection of matter in the merger and post-merger phases and is thus of high interest for multi-messenger astronomy. Today, a variety of nuclear EOSs are available with various underlying microphysical models. This calls for a study to focus on EOS effects from different physical nuclear matter properties and their influence on BNS mergers. We perform simulations of equal-mass BNS mergers with a set of 9 different EOSs based on Skyrme density functionals. In the models, we systematically vary the effective nucleon mass, incompressibility, and symmetry energy at saturation density. This allows us to investigate the influence of specific nuclear matter properties on the dynamics of BNS mergers. We analyze the impact of these properties on the merger dynamics, the fate of the remnant, disk formation, ejection of matter, and gravitational wave emission. Our results indicate that some aspects of the merger are sensitive to the EOS around saturation density while others are sensitive to the behavior towards higher densities, e.g., characterized by the slope of the pressure as a function of density. The detailed density dependence of the EOS thus needs to be taken into account to describe its influence on BNS mergers. △ Less

Submitted 21 December, 2023; v1 submitted 17 September, 2023; originally announced September 2023.

Comments: 17 pages, 19 figures, published version

Journal ref: MNRAS, 527, 8812 (2023)

Symmetry energy and neutron star properties constrained by chiral effective field theory calculations

Authors: Yeunhwan Lim, Achim Schwenk

Abstract: We investigate the nuclear symmetry energy and neutron star properties using a Bayesian analysis based on constraints from different chiral effective field theory calculations using new energy density functionals that allow for large variations at high densities. Constraints at high densities are included from observations of GW170817 and from NICER. In particular, we show that both NICER analyses… ▽ More We investigate the nuclear symmetry energy and neutron star properties using a Bayesian analysis based on constraints from different chiral effective field theory calculations using new energy density functionals that allow for large variations at high densities. Constraints at high densities are included from observations of GW170817 and from NICER. In particular, we show that both NICER analyses lead to very similar posterior results for the symmetry energy and neutron star properties when folded into our equation-of-state framework. Using the posteriors, we provide results for the symmetry energy and the slope parameter, as well as for the proton fraction, the speed of sound, and the central density in neutron stars. Moreover, we explore correlations of neutron star radii with the pressure and the speed of sound in neutron stars. Our 95\% credibility ranges for the symmetry energy $S_v$, the slope parameter $L$, and the radius of a 1.4$\msun$ neutron star, $R_{1.4}$, are $S_v=(30.6\text{--}33.9)$\,MeV, $L=(43.7\text{--}70.0)$\,MeV, and $R_{1.4}=(11.6\text{--}13.2)$\,km. Our analysis for the proton fraction shows that larger and/or heavier neutron stars are more likely to cool rapidly via the direct Urca process. Within our equation-of-state framework a maximum mass of neutron stars $M_{\rm max}>2.1\msun$ indicates that the speed of sound needs to exceed the conformal limit. △ Less

Submitted 10 March, 2024; v1 submitted 8 July, 2023; originally announced July 2023.

Comments: Published version, 12 pages, 13 figures

Journal ref: Phys. Rev. C 109, 035801 (2024)

Constraints on strong phase transitions in neutron stars

Authors: Tyler Gorda, Kai Hebeler, Aleksi Kurkela, Achim Schwenk, Aleksi Vuorinen

Abstract: We study current bounds on strong first-order phase transitions (PTs) along the equation of state (EOS) of dense strongly interacting matter in neutron stars, under the simplifying assumption that on either side of the PT the EOS can be approximated by a simple polytropic form. We construct a large ensemble of possible EOSs of this form, anchor them to chiral effective field theory calculations at… ▽ More We study current bounds on strong first-order phase transitions (PTs) along the equation of state (EOS) of dense strongly interacting matter in neutron stars, under the simplifying assumption that on either side of the PT the EOS can be approximated by a simple polytropic form. We construct a large ensemble of possible EOSs of this form, anchor them to chiral effective field theory calculations at nuclear density and perturbative QCD at high densities, and subject them to astrophysical constraints from high-mass pulsars and gravitational-wave observations. Within this setup, we find that a PT permits neutron-star solutions with larger radii, but only if the transition begins below twice nuclear saturation density. We also identify a large parameter space of allowed PTs currently unexplored by numerical-relativity studies. Additionally, we locate a small region of parameter space allowing twin-star solutions, though we find them to only marginally pass the current astrophysical constraints. Finally, we find that sizeable cores of high-density matter beyond the PT may be located in the centers of some stable neutron stars, primarily those with larger masses. △ Less

Submitted 8 August, 2023; v1 submitted 20 December, 2022; originally announced December 2022.

Comments: 10 pages, 11 figures; version accepted for publication in ApJ. Substantial change from v1 in treatment of matching to chiral EFT in a continuous manner. Figures and discussions updated. Conclusions qualitatively unchanged

Report number: HIP-2022-34/TH

Horizons: Nuclear Astrophysics in the 2020s and Beyond

Authors: H. Schatz, A. D. Becerril Reyes, A. Best, E. F. Brown, K. Chatziioannou, K. A. Chipps, C. M. Deibel, R. Ezzeddine, D. K. Galloway, C. J. Hansen, F. Herwig, A. P. Ji, M. Lugaro, Z. Meisel, D. Norman, J. S. Read, L. F. Roberts, A. Spyrou, I. Tews, F. X. Timmes, C. Travaglio, N. Vassh, C. Abia, P. Adsley, S. Agarwal , et al. (140 additional authors not shown)

Abstract: Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilit… ▽ More Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilities across an ever growing number of disciplines and subfields that need to be integrated. We take a holistic view of the field discussing the unique challenges and opportunities in nuclear astrophysics in regards to science, diversity, education, and the interdisciplinarity and breadth of the field. Clearly nuclear astrophysics is a dynamic field with a bright future that is entering a new era of discovery opportunities. △ Less

Submitted 16 May, 2022; originally announced May 2022.

Comments: 96 pages. Submitted to Journal of Physics G

Report number: LA-UR-22-23997

Nuclear equation of state for arbitrary proton fraction and temperature based on chiral effective field theory and a Gaussian process emulator

Authors: J. Keller, K. Hebeler, A. Schwenk

Abstract: We calculate the equation of state of asymmetric nuclear matter at finite temperature based on chiral effective field theory interactions to next-to-next-to-next-to-leading order. Our results assess the theoretical uncertainties from the many-body calculation and the chiral expansion. Using a Gaussian process emulator for the free energy, we derive the thermodynamic properties of matter through co… ▽ More We calculate the equation of state of asymmetric nuclear matter at finite temperature based on chiral effective field theory interactions to next-to-next-to-next-to-leading order. Our results assess the theoretical uncertainties from the many-body calculation and the chiral expansion. Using a Gaussian process emulator for the free energy, we derive the thermodynamic properties of matter through consistent derivatives and use the Gaussian process to access arbitrary proton fraction and temperature. This enables a first nonparametric calculation of the equation of state in beta equilibrium, and of the speed of sound and the symmetry energy at finite temperature. Moreover, our results show that the thermal part of the pressure decreases with increasing densities. △ Less

Submitted 20 February, 2023; v1 submitted 29 April, 2022; originally announced April 2022.

Comments: 7 pages, 5 figures, minor updates and improvements, published version

Journal ref: Phys. Rev. Lett. 130, 072701 (2023)

Coherent elastic neutrino-nucleus scattering: Terrestrial and astrophysical applications

Authors: M. Abdullah, H. Abele, D. Akimov, G. Angloher, D. Aristizabal-Sierra, C. Augier, A. B. Balantekin, L. Balogh, P. S. Barbeau, L. Baudis, A. L. Baxter, C. Beaufort, G. Beaulieu, V. Belov, A. Bento, L. Berge, I. A. Bernardi, J. Billard, A. Bolozdynya, A. Bonhomme, G. Bres, J-. L. Bret, A. Broniatowski, A. Brossard, C. Buck , et al. (250 additional authors not shown)

Abstract: Coherent elastic neutrino-nucleus scattering (CE$ν$NS) is a process in which neutrinos scatter on a nucleus which acts as a single particle. Though the total cross section is large by neutrino standards, CE$ν$NS has long proven difficult to detect, since the deposited energy into the nucleus is $\sim$ keV. In 2017, the COHERENT collaboration announced the detection of CE$ν$NS using a stopped-pion… ▽ More Coherent elastic neutrino-nucleus scattering (CE$ν$NS) is a process in which neutrinos scatter on a nucleus which acts as a single particle. Though the total cross section is large by neutrino standards, CE$ν$NS has long proven difficult to detect, since the deposited energy into the nucleus is $\sim$ keV. In 2017, the COHERENT collaboration announced the detection of CE$ν$NS using a stopped-pion source with CsI detectors, followed up the detection of CE$ν$NS using an Ar target. The detection of CE$ν$NS has spawned a flurry of activities in high-energy physics, inspiring new constraints on beyond the Standard Model (BSM) physics, and new experimental methods. The CE$ν$NS process has important implications for not only high-energy physics, but also astrophysics, nuclear physics, and beyond. This whitepaper discusses the scientific importance of CE$ν$NS, highlighting how present experiments such as COHERENT are informing theory, and also how future experiments will provide a wealth of information across the aforementioned fields of physics. △ Less

Submitted 14 March, 2022; originally announced March 2022.

Comments: contribution to Snowmasss 2021. Contact authors: P. S. Barbeau, R. Strauss, L. E. Strigari

A Next-Generation Liquid Xenon Observatory for Dark Matter and Neutrino Physics

Authors: J. Aalbers, K. Abe, V. Aerne, F. Agostini, S. Ahmed Maouloud, D. S. Akerib, D. Yu. Akimov, J. Akshat, A. K. Al Musalhi, F. Alder, S. K. Alsum, L. Althueser, C. S. Amarasinghe, F. D. Amaro, A. Ames, T. J. Anderson, B. Andrieu, N. Angelides, E. Angelino, J. Angevaare, V. C. Antochi, D. Antón Martin, B. Antunovic, E. Aprile, H. M. Araújo , et al. (572 additional authors not shown)

Abstract: The nature of dark matter and properties of neutrinos are among the most pressing issues in contemporary particle physics. The dual-phase xenon time-projection chamber is the leading technology to cover the available parameter space for Weakly Interacting Massive Particles (WIMPs), while featuring extensive sensitivity to many alternative dark matter candidates. These detectors can also study neut… ▽ More The nature of dark matter and properties of neutrinos are among the most pressing issues in contemporary particle physics. The dual-phase xenon time-projection chamber is the leading technology to cover the available parameter space for Weakly Interacting Massive Particles (WIMPs), while featuring extensive sensitivity to many alternative dark matter candidates. These detectors can also study neutrinos through neutrinoless double-beta decay and through a variety of astrophysical sources. A next-generation xenon-based detector will therefore be a true multi-purpose observatory to significantly advance particle physics, nuclear physics, astrophysics, solar physics, and cosmology. This review article presents the science cases for such a detector. △ Less

Submitted 4 March, 2022; originally announced March 2022.

Comments: 77 pages, 40 figures, 1262 references

Report number: INT-PUB-22-003

Journal ref: J. Phys. G: Nucl. Part. Phys. 50 (2023) 013001

Constraining Neutron-Star Matter with Microscopic and Macroscopic Collisions

Authors: S. Huth, P. T. H. Pang, I. Tews, T. Dietrich, A. Le Fèvre, A. Schwenk, W. Trautmann, K. Agarwal, M. Bulla, M. W. Coughlin, C. Van Den Broeck

Abstract: 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 exper… ▽ More 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. △ Less

Submitted 24 June, 2022; v1 submitted 13 July, 2021; originally announced July 2021.

Comments: 7 pages, 2 figures, Supplemental Material, minor changes

Report number: LA-UR-21-22072

Journal ref: Nature 606, 276-280 (2022)

A Detailed Examination of Astrophysical Constraints on the Symmetry Energy and the Neutron Skin of $^{208}$Pb with Minimal Modeling Assumptions

Authors: Reed Essick, Philippe Landry, Achim Schwenk, Ingo Tews

Abstract: The symmetry energy and its density dependence are pivotal 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\pm0.071$ fm for the neutron-skin thickness of $^{208}$Pb, $R_{\rm skin}^{^{208}\text{Pb}}$, implying a sy… ▽ More The symmetry energy and its density dependence are pivotal 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\pm0.071$ fm for the neutron-skin thickness of $^{208}$Pb, $R_{\rm skin}^{^{208}\text{Pb}}$, implying a symmetry-energy slope parameter $L$ of $106\pm37$ 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_{\rm skin}^{^{208}\text{Pb}}$ directly from observations of neutron stars with minimal modeling assumptions. The resulting astrophysical constraints from heavy pulsar masses, LIGO/Virgo, and NICER favor smaller values of the neutron skin and $L$, as well as negative symmetry incompressibilities. Combining astrophysical data with chiral effective field theory ($χ$EFT) and PREX-II constraints yields $S_0 = 33.0^{+2.0}_{-1.8}$ MeV, $L=53^{+14}_{-15}$ MeV, and $R_{\rm skin}^{^{208}\text{Pb}} = 0.17^{+0.04}_{-0.04}$ fm. We also examine the consistency of several individual $χ$EFT calculations with astrophysical observations and terrestrial experiments. We find that there is only mild tension between $χ$EFT, astrophysical data, and PREX-II's $R_\mathrm{skin}^{^{208}\mathrm{Pb}}$ measurement ($p$-value $= 12.3\%$) and that there is excellent agreement between $χ$EFT, astrophysical data, and other nuclear experiments. △ Less

Submitted 15 December, 2021; v1 submitted 12 July, 2021; originally announced July 2021.

Comments: 18 pages, 12 figures, 1 table

Report number: LA-UR-21-21363

Journal ref: Phys. Rev. C 104, 065804 (2021)

Constraints on the dense matter equation of state and neutron star properties from NICER's mass-radius estimate of PSR J0740+6620 and multimessenger observations

Authors: G. Raaijmakers, S. K. Greif, K. Hebeler, T. Hinderer, S. Nissanke, A. Schwenk, T. E. Riley, A. L. Watts, J. M. Lattimer, W. C. G. Ho

Abstract: 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… ▽ More 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. △ Less

Submitted 13 June, 2021; v1 submitted 14 May, 2021; originally announced May 2021.

Comments: 17 pages, 8 figures; accepted for publication in the Astrophysical Journal Letters

Journal ref: The Astrophysical Journal Letters, Volume 918 (2021), L29

Astrophysical Constraints on the Symmetry Energy and the Neutron Skin of $^{208}$Pb with Minimal Modeling Assumptions

Authors: Reed Essick, Ingo Tews, Philippe Landry, Achim Schwenk

Abstract: 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 \pm 0.071$ fm for the neutron-skin thickness of $^{208}$Pb, implying a slope parameter… ▽ More 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 \pm 0.071$ fm for the neutron-skin thickness of $^{208}$Pb, implying a slope parameter $L = 106 \pm 37$ 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^{+2.0}_{-1.8}$ MeV, $L=53^{+14}_{-15}$ MeV, and $R_\mathrm{skin}^{^{208}\mathrm{Pb}}=0.17^{+0.04}_{-0.04}$ fm. △ Less

Submitted 3 November, 2021; v1 submitted 19 February, 2021; originally announced February 2021.

Comments: 11 pages, 6 figures, 2 tables

Report number: LA-UR-21-20527

Neutron matter at finite temperature based on chiral effective field theory interactions

Authors: J. Keller, C. Wellenhofer, K. Hebeler, A. Schwenk

Abstract: 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 inc… ▽ More 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. △ Less

Submitted 17 May, 2021; v1 submitted 11 November, 2020; originally announced November 2020.

Comments: 14 pages, 14 figures, minor changes, published version

Journal ref: Phys. Rev. C 103, 055806 (2021)

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

Authors: S. Huth, C. Wellenhofer, A. Schwenk

Abstract: 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 simul… ▽ More 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. △ Less

Submitted 12 February, 2021; v1 submitted 18 September, 2020; originally announced September 2020.

Comments: 21 pages, 19 figures, minor changes, published version

Journal ref: Phys. Rev. C 103, 025803 (2021)

Equation of state constraints from nuclear physics, neutron star masses, and future moment of inertia measurements

Authors: S. K. Greif, K. Hebeler, J. M. Lattimer, C. J. Pethick, A. Schwenk

Abstract: 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 s… ▽ More 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. △ Less

Submitted 7 October, 2020; v1 submitted 28 May, 2020; originally announced May 2020.

Comments: 11 pages, 10 figures, published version

Report number: NORDITA 2020-046

Journal ref: Astrophys. J. 901, 155 (2020)

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

Authors: G. Raaijmakers, S. K. Greif, T. E. Riley, T. Hinderer, K. Hebeler, A. Schwenk, A. L. Watts, S. Nissanke, S. Guillot, J. M. Lattimer, R. M. Ludlam

Abstract: 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 s… ▽ More 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. △ Less

Submitted 17 April, 2020; v1 submitted 23 December, 2019; originally announced December 2019.

Comments: 18 pages, 8 figures. Accepted for publication in ApJ Letters

Journal ref: The Astrophysical Journal Letters, 893, L21 (2020)

A NICER view of PSR J0030+0451: Implications for the dense matter equation of state

Authors: G. Raaijmakers, T. E. Riley, A. L. Watts, S. K. Greif, S. M. Morsink, K. Hebeler, A. Schwenk, T. Hinderer, S. Nissanke, S. Guillot, Z. Arzoumanian, S. Bogdanov, D. Chakrabarty K. C. Gendreau, W. C. G. Ho, J. M. Lattimer, R. M. Ludlam, M. T. Wolff

Abstract: Both the mass and radius of the millisecond pulsar PSR J0030+0451 have been inferred via pulse-profile modeling of X-ray data obtained by NASA's NICER mission. In this Letter we study the implications of the mass-radius inference reported for this source by Riley et al. (2019) for the dense matter equation of state (EOS), in the context of prior information from nuclear physics at low densities. U… ▽ More Both the mass and radius of the millisecond pulsar PSR J0030+0451 have been inferred via pulse-profile modeling of X-ray data obtained by NASA's NICER mission. In this Letter we study the implications of the mass-radius inference reported for this source by Riley et al. (2019) for the dense matter equation of state (EOS), in the context of prior information from nuclear physics at low densities. Using a Bayesian framework we infer central densities and EOS properties for two choices of high-density extensions: a piecewise-polytropic model and a model based on assumptions of the speed of sound in dense matter. Around nuclear saturation density these extensions are matched to an EOS uncertainty band obtained from calculations based on chiral effective field theory interactions, which provide a realistic description of atomic nuclei as well as empirical nuclear matter properties within uncertainties. We further constrain EOS expectations with input from the current highest measured pulsar mass; together, these constraints offer a narrow Bayesian prior informed by theory as well as laboratory and astrophysical measurements. The NICER mass-radius likelihood function derived by Riley et al. (2019) using pulse-profile modeling is consistent with the highest-density region of this prior. The present relatively large uncertainties on mass and radius for PSR J0030+0451 offer, however, only a weak posterior information gain over the prior. We explore the sensitivity to the inferred geometry of the heated regions that give rise to the pulsed emission, and find a small increase in posterior gain for an alternative (but less preferred) model. Lastly, we investigate the hypothetical scenario of increasing the NICER exposure time for PSR J0030+0451. △ Less

Submitted 11 December, 2019; originally announced December 2019.

Comments: Appears in ApJ Letters Focus Issue on NICER Constraints on the Dense Matter Equation of State, 17 pages, 5 figures

Journal ref: ApJ Letters, 887, L22 (2019)

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

Authors: I. Tews, A. Schwenk

Abstract: 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… ▽ More 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. △ Less

Submitted 24 March, 2020; v1 submitted 7 August, 2019; originally announced August 2019.

Comments: 10 pages, 3 figures, published version

Report number: LA-UR-19-27912

Journal ref: Astrophys. J. 892, 14 (2020)

${\it Ab~initio}$ constraints on thermal effects of the nuclear equation of state

Authors: Arianna Carbone, Achim Schwenk

Abstract: 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… ▽ More 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. △ Less

Submitted 22 August, 2019; v1 submitted 1 April, 2019; originally announced April 2019.

Comments: 11 pages, 7 figures, added references

Journal ref: Phys. Rev. C 100, 025805 (2019)

Equation of state sensitivities when inferring neutron star and dense matter properties

Authors: S. K. Greif, G. Raaijmakers, K. Hebeler, A. Schwenk, A. L. Watts

Abstract: 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… ▽ More 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. △ Less

Submitted 6 March, 2019; v1 submitted 19 December, 2018; originally announced December 2018.

Comments: to appear in MNRAS

Dense matter with eXTP

Authors: Anna L. Watts, Wenfei Yu, Juri Poutanen, Shu Zhang, Sudip Bhattacharyya, Slavko Bogdanov, Long Ji, Alessandro Patruno, Thomas E. Riley, Pavel Bakala, Altan Baykal, Federico Bernardini, Ignazio Bombaci, Edward Brown, Yuri Cavecchi, Deepto Chakrabarty, Jérôme Chenevez, Nathalie Degenaar, Melania Del Santo, Tiziana Di Salvo, Victor Doroshenko, Maurizio Falanga, Robert D. Ferdman, Marco Feroci, Angelo F. Gambino , et al. (51 additional authors not shown)

Abstract: 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, b… ▽ More 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. △ Less

Submitted 10 December, 2018; originally announced December 2018.

Comments: Accepted for publication on Sci. China Phys. Mech. Astron. (2019)

Equation of state effects in core-collapse supernovae

Authors: Hannah Yasin, Sabrina Schäfer, Almudena Arcones, Achim Schwenk

Abstract: We investigate the impact of different properties of the nuclear equation of state in core-collapse supernovae, with a focus on the proto-neutron-star contraction and its impact on the shock evolution. To this end, we introduce a range of equations of state that vary the nucleon effective mass, incompressibility, symmetry energy, and nuclear saturation point. This allows us to point to the differe… ▽ More We investigate the impact of different properties of the nuclear equation of state in core-collapse supernovae, with a focus on the proto-neutron-star contraction and its impact on the shock evolution. To this end, we introduce a range of equations of state that vary the nucleon effective mass, incompressibility, symmetry energy, and nuclear saturation point. This allows us to point to the different effects in changing these properties from the Lattimer and Swesty to the Shen et al. equations of state, the two most commonly used equations of state in simulations. In particular, we trace the contraction behavior to the effective mass, which determines the thermal nucleonic contributions to the equation of state. Larger effective masses lead to lower pressures at nuclear densities and a lower thermal index. This results in a more rapid contraction of the proto-neutron star and consequently higher neutrino energies, which aids the shock evolution to a faster explosion. △ Less

Submitted 12 March, 2020; v1 submitted 5 December, 2018; originally announced December 2018.

Comments: 6 pages, 4 figures, minor changes, published version

Journal ref: Phys. Rev. Lett. 124, 092701 (2020)

First results on the scalar WIMP-pion coupling, using the XENON1T experiment

Authors: E. Aprile, J. Aalbers, F. Agostini, M. Alfonsi, L. Althueser, F. D. Amaro, M. Anthony, V. C. Antochi, F. Arneodo, L. Baudis, B. Bauermeister, M. L. Benabderrahmane, T. Berger, P. A. Breur, A. Brown, A. Brown, E. Brown, S. Bruenner, G. Bruno, R. Budnik, C. Capelli, J. M. R. Cardoso, D. Cichon, D. Coderre, A. P. Colijn , et al. (107 additional authors not shown)

Abstract: We present first results on the scalar WIMP-pion coupling from 1 t$\times$yr of exposure with the XENON1T experiment. This interaction is generated when the WIMP couples to a virtual pion exchanged between the nucleons in a nucleus. In contrast to most non-relativistic operators, these pion-exchange currents can be coherently enhanced by the total number of nucleons, and therefore may dominate in… ▽ More We present first results on the scalar WIMP-pion coupling from 1 t$\times$yr of exposure with the XENON1T experiment. This interaction is generated when the WIMP couples to a virtual pion exchanged between the nucleons in a nucleus. In contrast to most non-relativistic operators, these pion-exchange currents can be coherently enhanced by the total number of nucleons, and therefore may dominate in scenarios where spin-independent WIMP-nucleon interactions are suppressed. Moreover, for natural values of the couplings, they dominate over the spin-dependent channel due to their coherence in the nucleus. Using the signal model of this new WIMP-pion channel, no significant excess is found, leading to an upper limit cross section of $6.4\times10^{-46}$ cm$^2$ (90 % confidence level) at 30 GeV/c$^2$ WIMP mass. △ Less

Submitted 22 February, 2019; v1 submitted 29 November, 2018; originally announced November 2018.

Report number: INT-PUB-18-058

Journal ref: Phys. Rev. Lett. 122, 071301 (2019)

Discriminating WIMP-nucleus response functions in present and future XENON-like direct detection experiments

Authors: A. Fieguth, M. Hoferichter, P. Klos, J. Menéndez, A. Schwenk, C. Weinheimer

Abstract: 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, whic… ▽ More 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. △ Less

Submitted 31 May, 2018; v1 submitted 12 February, 2018; originally announced February 2018.

Comments: 10 pages, 7 figures

Report number: INT-PUB-18-006

Journal ref: Phys. Rev. D 97, 103532 (2018)

Chiral interactions up to next-to-next-to-next-to-leading order and nuclear saturation

Authors: C. Drischler, K. Hebeler, A. Schwenk

Abstract: We present an efficient Monte Carlo framework for perturbative calculations of infinite nuclear matter based on chiral two-, three-, and four-nucleon interactions. The method enables the incorporation of all many-body contributions in a straightforward and transparent way, and makes it possible to extract systematic uncertainty estimates by performing order-by-order calculations in the chiral expa… ▽ More We present an efficient Monte Carlo framework for perturbative calculations of infinite nuclear matter based on chiral two-, three-, and four-nucleon interactions. The method enables the incorporation of all many-body contributions in a straightforward and transparent way, and makes it possible to extract systematic uncertainty estimates by performing order-by-order calculations in the chiral expansion as well as the many-body expansion. The versatility of this new framework is demonstrated by applying it to chiral low-momentum interactions, exhibiting a very good many-body convergence up to fourth order. Following these benchmarks, we explore new chiral interactions up to next-to-next-to-next-to-leading order (N$^3$LO). Remarkably, simultaneous fits to the triton and to saturation properties can be achieved, while all three-nucleon low-energy couplings remain natural. The theoretical uncertainties of nuclear matter are significantly reduced when going from next-to-next-to-leading order to N$^3$LO. △ Less

Submitted 1 February, 2019; v1 submitted 23 October, 2017; originally announced October 2017.

Comments: published version, incl. supplemental material

Journal ref: Phys. Rev. Lett. 122, 042501 (2019)

Improved limits for Higgs-portal dark matter from LHC searches

Authors: Martin Hoferichter, Philipp Klos, Javier Menéndez, Achim Schwenk

Abstract: Searches for invisible Higgs decays at the Large Hadron Collider constrain dark matter Higgs-portal models, where dark matter interacts with the Standard Model fields via the Higgs boson. While these searches complement dark matter direct-detection experiments, a comparison of the two limits depends on the coupling of the Higgs boson to the nucleons forming the direct-detection nuclear target, typ… ▽ More Searches for invisible Higgs decays at the Large Hadron Collider constrain dark matter Higgs-portal models, where dark matter interacts with the Standard Model fields via the Higgs boson. While these searches complement dark matter direct-detection experiments, a comparison of the two limits depends on the coupling of the Higgs boson to the nucleons forming the direct-detection nuclear target, typically parameterized in a single quantity $f_N$. We evaluate $f_N$ using recent phenomenological and lattice-QCD calculations, and include for the first time the coupling of the Higgs boson to two nucleons via pion-exchange currents. We observe a partial cancellation for Higgs-portal models that makes the two-nucleon contribution anomalously small. Our results, summarized as $f_N=0.308(18)$, show that the uncertainty of the Higgs-nucleon coupling has been vastly overestimated in the past. The improved limits highlight that state-of-the-art nuclear physics input is key to fully exploiting experimental searches. △ Less

Submitted 20 October, 2017; v1 submitted 7 August, 2017; originally announced August 2017.

Comments: 6 pages, 3 figures; version to appear in PRL

Report number: INT-PUB-17-031

Journal ref: Phys. Rev. Lett. 119, 181803 (2017)

Neutrino-nucleon scattering in supernova matter from the virial expansion

Authors: C. J. Horowitz, O. L. Caballero, Zidu Lin, Evan O'Connor, A. Schwenk

Abstract: We generalize our virial approach to study the neutral current neutrino response of nuclear matter at low densities. In the long-wavelength limit, the virial expansion makes model-independent predictions for neutrino-nucleon scattering rates and the density S_V and spin S_A responses. We find S_A is significantly reduced from one even at low densities. We provide a simple fit S_A^f(n,T,Y_p) of the… ▽ More We generalize our virial approach to study the neutral current neutrino response of nuclear matter at low densities. In the long-wavelength limit, the virial expansion makes model-independent predictions for neutrino-nucleon scattering rates and the density S_V and spin S_A responses. We find S_A is significantly reduced from one even at low densities. We provide a simple fit S_A^f(n,T,Y_p) of the axial response as a function of density n, temperature T and proton fraction Y_p. This fit reproduces our model independent virial results at low densities and reproduces the Burrows and Sawyer random phase approximation (RPA) results at high densities. Our fit can be incorporated into supernova simulations in a straight forward manner. Preliminary one dimensional supernova simulations suggest that the reduction in the axial response may enhance neutrino heating rates in the gain region during the accretion phase of a core-collapse supernovae. △ Less

Submitted 15 November, 2016; originally announced November 2016.

Comments: 7 pages, 3 figures

Journal ref: Phys. Rev. C 95, 025801 (2017)

Pairing in neutron matter: New uncertainty estimates and three-body forces

Authors: C. Drischler, T. Krüger, K. Hebeler, A. Schwenk

Abstract: We present solutions of the BCS gap equation in the channels ${}^1S_0$ and ${}^3P_2-{}^3F_2$ in neutron matter based on nuclear interactions derived within chiral effective field theory (EFT). Our studies are based on a representative set of nonlocal nucleon-nucleon (NN) plus three-nucleon (3N) interactions up to next-to-next-to-next-to-leading order (N$^3$LO) as well as local and semilocal chiral… ▽ More We present solutions of the BCS gap equation in the channels ${}^1S_0$ and ${}^3P_2-{}^3F_2$ in neutron matter based on nuclear interactions derived within chiral effective field theory (EFT). Our studies are based on a representative set of nonlocal nucleon-nucleon (NN) plus three-nucleon (3N) interactions up to next-to-next-to-next-to-leading order (N$^3$LO) as well as local and semilocal chiral NN interactions up to N$^2$LO and N$^4$LO, respectively. In particular, we investigate for the first time the impact of subleading 3N forces at N$^3$LO on pairing gaps and also derive uncertainty estimates by taking into account results for pairing gaps at different orders in the chiral expansion. Finally, we discuss different methods for obtaining self-consistent solutions of the gap equation. Besides the widely-used quasi-linear method by Khodel et al. we demonstrate that the modified Broyden method is well applicable and exhibits a robust convergence behavior. In contrast to Khodel's method it is based on a direct iteration of the gap equation without imposing an auxiliary potential and is straightforward to implement. △ Less

Submitted 17 October, 2016; originally announced October 2016.

Journal ref: Phys. Rev. C 95, 024302 (2017)

Neutron matter from chiral two- and three-nucleon calculations up to N$^3$LO

Authors: C. Drischler, A. Carbone, K. Hebeler, A. Schwenk

Abstract: Neutron matter is an ideal laboratory for nuclear interactions derived from chiral effective field theory since all contributions are predicted up to next-to-next-to-next-to-leading order (N$^3$LO) in the chiral expansion. By making use of recent advances in the partial-wave decomposition of three- nucleon (3N) forces, we include for the first time N$^3$LO 3N interactions in many-body perturbation… ▽ More Neutron matter is an ideal laboratory for nuclear interactions derived from chiral effective field theory since all contributions are predicted up to next-to-next-to-next-to-leading order (N$^3$LO) in the chiral expansion. By making use of recent advances in the partial-wave decomposition of three- nucleon (3N) forces, we include for the first time N$^3$LO 3N interactions in many-body perturbation theory (MBPT) up to third order and in self-consistent Green's function theory (SCGF). Using these two complementary many-body frameworks we provide improved predictions for the equation of state of neutron matter at zero temperature and also analyze systematically the many-body convergence for different chiral EFT interactions. Furthermore, we present an extension of the normal-ordering framework to finite temperatures. These developments open the way to improved calculations of neutron-rich matter including estimates of theoretical uncertainties for astrophysical applications. △ Less

Submitted 9 December, 2016; v1 submitted 19 August, 2016; originally announced August 2016.

Comments: minor changes, published version

Journal ref: Phys. Rev. C 94, 054307 (2016)

Impact of Nucleon-Nucleon Bremsstrahlung Rates Beyond One-Pion Exchange

Authors: Alexander Bartl, Robert Bollig, Hans-Thomas Janka, Achim Schwenk

Abstract: Neutrino-pair production and annihilation through nucleon-nucleon bremsstrahlung is included in current supernova simulations by rates that are based on the one-pion-exchange approximation. Here we explore the consequences of bremsstrahlung rates based on a modern nuclear interactions for proto-neutron star cooling and the corresponding neutrino emission. We find that despite a reduction of the br… ▽ More Neutrino-pair production and annihilation through nucleon-nucleon bremsstrahlung is included in current supernova simulations by rates that are based on the one-pion-exchange approximation. Here we explore the consequences of bremsstrahlung rates based on a modern nuclear interactions for proto-neutron star cooling and the corresponding neutrino emission. We find that despite a reduction of the bremsstrahlung emission by a factor of 2-5 in the neutrinospheric region, models with the improved treatment exhibit only $\lesssim$5% changes of the neutrino luminosities and an increase of $\lesssim$0.7 MeV of the average energies of the radiated neutrino spectra, with the largest effects for the antineutrinos of all flavors and at late times. Overall, the proto-neutron star cooling evolution is slowed down modestly by $\lesssim$0.5-1 s. △ Less

Submitted 14 October, 2016; v1 submitted 17 August, 2016; originally announced August 2016.

Comments: 11 pages, 7 figures, minor changes and additions, to appear in Phys. Rev. D

Report number: INT-PUB-16-022

Journal ref: Phys. Rev. D 94, 083009 (2016)

Analysis strategies for general spin-independent WIMP-nucleus scattering

Authors: Martin Hoferichter, Philipp Klos, Javier Menéndez, Achim Schwenk

Abstract: We propose a formalism for the analysis of direct-detection dark-matter searches that covers all coherent responses for scalar and vector interactions and incorporates QCD constraints imposed by chiral symmetry, including all one- and two-body WIMP-nucleon interactions up to third order in chiral effective field theory. One of the free parameters in the WIMP-nucleus cross section corresponds to st… ▽ More We propose a formalism for the analysis of direct-detection dark-matter searches that covers all coherent responses for scalar and vector interactions and incorporates QCD constraints imposed by chiral symmetry, including all one- and two-body WIMP-nucleon interactions up to third order in chiral effective field theory. One of the free parameters in the WIMP-nucleus cross section corresponds to standard spin-independent searches, but in general different combinations of new-physics couplings are probed. We identify the interference with the isovector counterpart of the standard spin-independent response and two-body currents as the dominant corrections to the leading spin-independent structure factor, and discuss the general consequences for the interpretation of direct-detection experiments, including minimal extensions of the standard spin-independent analysis. Fits for all structure factors required for the scattering off xenon targets are provided based on state-of-the-art nuclear shell-model calculations. △ Less

Submitted 8 August, 2016; v1 submitted 25 May, 2016; originally announced May 2016.

Comments: 20 pages, 7 figures; v2 includes the two-body current related to the trace anomaly of the QCD energy-momentum tensor

Report number: INT-PUB-16-012

Journal ref: Phys. Rev. D 94, 063505 (2016)

Measuring the neutron star equation of state using X-ray timing

Authors: Anna L. Watts, Nils Andersson, Deepto Chakrabarty, Marco Feroci, Kai Hebeler, Gianluca Israel, Frederick K. Lamb, M. Coleman Miller, Sharon Morsink, Feryal Özel, Alessandro Patruno, Juri Poutanen, Dimitrios Psaltis, Achim Schwenk, Andrew W. Steiner, Luigi Stella, Laura Tolos, Michiel van der Klis

Abstract: One of the primary science goals of the next generation of hard X-ray timing instruments is to determine the equation of state of the matter at supranuclear densities inside neutron stars, by measuring the radius of neutron stars with different masses to accuracies of a few percent. Three main techniques can be used to achieve this goal. The first involves waveform modelling. The flux we observe f… ▽ More One of the primary science goals of the next generation of hard X-ray timing instruments is to determine the equation of state of the matter at supranuclear densities inside neutron stars, by measuring the radius of neutron stars with different masses to accuracies of a few percent. Three main techniques can be used to achieve this goal. The first involves waveform modelling. The flux we observe from a hotspot on the neutron star surface offset from the rotational pole will be modulated by the star's rotation, giving rise to a pulsation. Information about mass and radius is encoded into the pulse profile via relativistic effects, and tight constraints on mass and radius can be obtained. The second technique involves characterising the spin distribution of accreting neutron stars. The most rapidly rotating stars provide a very clean constraint, since the mass-shedding limit is a function of mass and radius. However the overall spin distribution also provides a guide to the torque mechanisms in operation and the moment of inertia, both of which can depend sensitively on dense matter physics. The third technique is to search for quasi-periodic oscillations in X-ray flux associated with global seismic vibrations of magnetars (the most highly magnetized neutron stars), triggered by magnetic explosions. The vibrational frequencies depend on stellar parameters including the dense matter equation of state. We illustrate how these complementary X-ray timing techniques can be used to constrain the dense matter equation of state, and discuss the results that might be expected from a 10m$^2$ instrument. We also discuss how the results from such a facility would compare to other astronomical investigations of neutron star properties. [Modified for arXiv] △ Less

Submitted 2 February, 2016; originally announced February 2016.

Comments: To appear in Reviews of Modern Physics as a Colloquium, 23 pages, 9 figures

Asymmetric nuclear matter based on chiral two- and three-nucleon interactions

Authors: Christian Drischler, Kai Hebeler, Achim Schwenk

Abstract: We calculate the properties of isospin-asymmetric nuclear matter based on chiral nucleon-nucleon (NN) and three-nucleon (3N) interactions. To this end, we develop an improved normal-ordering framework that allows to include general 3N interactions starting from a plane-wave partial-wave-decomposed form. We present results for the energy per particle for general isospin asymmetries based on a set o… ▽ More We calculate the properties of isospin-asymmetric nuclear matter based on chiral nucleon-nucleon (NN) and three-nucleon (3N) interactions. To this end, we develop an improved normal-ordering framework that allows to include general 3N interactions starting from a plane-wave partial-wave-decomposed form. We present results for the energy per particle for general isospin asymmetries based on a set of different Hamiltonians, study their saturation properties, the incompressibility, symmetry energy, and also provide an analytic parametrization for the energy per particle as a function of density and isospin asymmetry. △ Less

Submitted 9 December, 2016; v1 submitted 22 October, 2015; originally announced October 2015.

Comments: minor changes, published version

Journal ref: Phys. Rev. C 93, 054314 (2016)

Bose-Einstein condensates in neutron stars

Authors: C. J. Pethick, Thomas Schaefer, A. Schwenk

Abstract: 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 star… ▽ More 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. △ Less

Submitted 21 July, 2015; originally announced July 2015.

Comments: 22 pages, 2 figures, to appear in "Universal Themes of Bose-Einstein Condensation", edited by D.W. Snoke, N.P. Proukakis and P.B. Littlewood (Cambridge University Press)

Chiral power counting of one- and two-body currents in direct detection of dark matter

Authors: Martin Hoferichter, Philipp Klos, Achim Schwenk

Abstract: 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 interact… ▽ More 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. △ Less

Submitted 20 May, 2015; v1 submitted 16 March, 2015; originally announced March 2015.

Comments: 8 pages, 1 table; journal version

Nuclear structure aspects of spin-independent WIMP scattering off xenon

Authors: L. Vietze, P. Klos, J. Menéndez, W. C. Haxton, A. Schwenk

Abstract: We study the structure factors for spin-independent WIMP scattering off xenon based on state-of-the-art large-scale shell-model calculations, which are shown to yield a good spectroscopic description of all experimentally relevant isotopes. Our results are based on the leading scalar one-body currents only. At this level and for the momentum transfers relevant to direct dark matter detection, the… ▽ More We study the structure factors for spin-independent WIMP scattering off xenon based on state-of-the-art large-scale shell-model calculations, which are shown to yield a good spectroscopic description of all experimentally relevant isotopes. Our results are based on the leading scalar one-body currents only. At this level and for the momentum transfers relevant to direct dark matter detection, the structure factors are in very good agreement with the phenomenological Helm form factors used to give experimental limits for WIMP-nucleon cross sections. In contrast to spin-dependent WIMP scattering, the spin-independent channel, at the one-body level, is less sensitive to nuclear structure details. In addition, we explicitly show that the structure factors for inelastic scattering are suppressed by ~ 10^{-4} compared to the coherent elastic scattering response. This implies that the detection of inelastic scattering will be able to discriminate clearly between spin-independent and spin-dependent scattering. Finally, we provide fits for all calculated structure factors. △ Less

Submitted 19 February, 2015; v1 submitted 5 December, 2014; originally announced December 2014.

Comments: 10 pages, 24 figures, minor changes, published version

Journal ref: Physical Review D 91, 043520 (2015)

Neutrino-pair bremsstrahlung from nucleon-$α$ versus nucleon-nucleon scattering

Authors: Rishi Sharma, Sonia Bacca, Achim Schwenk

Abstract: We study the impact of the nucleon-$α$ P-wave resonances on neutrino-pair bremsstrahlung. Because of the non-central spin-orbit interaction, these resonances lead to an enhanced contribution to the nucleon spin structure factor for temperatures $T \lesssim 4$ MeV. If the $α$-particle fraction is significant and the temperature is in this range, this contribution is competitive with neutron-neutron… ▽ More We study the impact of the nucleon-$α$ P-wave resonances on neutrino-pair bremsstrahlung. Because of the non-central spin-orbit interaction, these resonances lead to an enhanced contribution to the nucleon spin structure factor for temperatures $T \lesssim 4$ MeV. If the $α$-particle fraction is significant and the temperature is in this range, this contribution is competitive with neutron-neutron bremsstrahlung. This may be relevant for neutrino production in core-collapse supernovae or other dense astrophysical environments. Similar enhancements are expected for resonant non-central nucleon-nucleus interactions. △ Less

Submitted 12 November, 2014; originally announced November 2014.

Comments: 6 pages, 4 figures

Report number: TIFR/TH/14-22

The Large Observatory For x-ray Timing

Authors: M. Feroci, J. W. den Herder, E. Bozzo, D. Barret, S. Brandt, M. Hernanz, M. van der Klis, M. Pohl, A. Santangelo, L. Stella, A. Watts, J. Wilms, S. Zane, M. Ahangarianabhari, C. Albertus, M. Alford, A. Alpar, D. Altamirano, L. Alvarez, L. Amati, C. Amoros, N. Andersson, A. Antonelli, A. Argan, R. Artigue , et al. (320 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 29 August, 2014; v1 submitted 27 August, 2014; originally announced August 2014.

Comments: Proc. SPIE 9144, Space Telescopes and Instrumentation 2014: Ultraviolet to Gamma Ray, 91442T

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

Authors: T. Krüger, K. Hebeler, A. Schwenk

Abstract: 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 o… ▽ More 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. △ Less

Submitted 25 March, 2015; v1 submitted 18 August, 2014; originally announced August 2014.

Comments: 5 pages, 4 figures, minor changes, published version

Report number: INT-PUB-14-029

Journal ref: Phys. Lett. B 744, 18 (2015)

Charged-current reactions in the supernova neutrino-sphere

Authors: Ermal Rrapaj, J. W. Holt, Alexander Bartl, Sanjay Reddy, A. Schwenk

Abstract: 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 interac… ▽ More 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. △ Less

Submitted 1 April, 2015; v1 submitted 14 August, 2014; originally announced August 2014.

Report number: INT-PUB-14-028

Journal ref: Phys. Rev. C 91, 035806 (2015)

Pairing and superfluidity of nucleons in neutron stars

Authors: A. Gezerlis, C. J. Pethick, A. Schwenk

Abstract: 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 inspi… ▽ More 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. △ Less

Submitted 13 April, 2015; v1 submitted 23 June, 2014; originally announced June 2014.

Comments: 48 pages, 13 figures; as published in "Novel Superfluids, Volume 2", edited by K. H. Bennemann and J. B. Ketterson (Oxford University Press)

Report number: NORDITA-2014-77

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

Authors: A. Bartl, C. J. Pethick, A. Schwenk

Abstract: 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 scat… ▽ More 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. △ Less

Submitted 3 September, 2014; v1 submitted 17 March, 2014; originally announced March 2014.

Comments: 5 pages, 4 figures, published version, NORDITA-2014-28

Journal ref: Phys. Rev. Lett. 113, 081101 (2014)

Microscopic calculations and energy expansions for neutron-rich matter

Authors: C. Drischler, V. Soma, A. Schwenk

Abstract: 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 matt… ▽ More 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. △ Less

Submitted 5 March, 2014; v1 submitted 21 October, 2013; originally announced October 2013.

Comments: 15 pages, 9 figures, published version

Journal ref: Phys. Rev. C 89, 025806 (2014)

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

Authors: I. Tews, T. Krüger, A. Gezerlis, K. Hebeler, A. Schwenk

Abstract: 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. 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. △ Less

Submitted 14 October, 2013; originally announced October 2013.

Comments: 10 pages, 7 figures, proceedings of the International Conference on Nuclear Theory in the Supercomputing Era 2013, Iowa State University, May 13-17, 2013

Signatures of Dark Matter Scattering Inelastically Off Nuclei

Authors: L. Baudis, G. Kessler, P. Klos, R. F. Lang, J. Menendez, S. Reichard, A. Schwenk

Abstract: 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 ef… ▽ More 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. △ Less

Submitted 8 November, 2013; v1 submitted 3 September, 2013; originally announced September 2013.

Comments: 9 pages, 7 figures. Matches resubmitted version to Phys. Rev. D. One figure added; supplemental material (fits to the structure functions) added as an Appendix

Journal ref: Phys. Rev. D 88, 115014 (2013)

The chiral condensate in neutron matter

Authors: T. Krüger, I. Tews, B. Friman, K. Hebeler, A. Schwenk

Abstract: 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 restorati… ▽ More 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. △ Less

Submitted 30 September, 2013; v1 submitted 8 July, 2013; originally announced July 2013.

Comments: published version, 6 pages, 4 figures

Journal ref: Phys. Lett. B 726, 412 (2013)

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

Authors: P. Klos, J. Menéndez, D. Gazit, A. Schwenk

Abstract: 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 s… ▽ More 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. △ Less

Submitted 23 January, 2015; v1 submitted 29 April, 2013; originally announced April 2013.

Comments: 25 pages, 13 figures. Version v4 with small typos corrected in Eqs.(8),(10),(20) and in the Appendix

Journal ref: Physical Review D 88, 083516 (2013)

Neutron matter from chiral effective field theory interactions

Authors: T. Krüger, I. Tews, K. Hebeler, A. Schwenk

Abstract: 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… ▽ More 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. △ Less

Submitted 12 August, 2013; v1 submitted 8 April, 2013; originally announced April 2013.

Comments: published version; 21 pages, 11 figures, 5 tables

Journal ref: Phys. Rev. C 88, 025802 (2013)