-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
Three-Nucleon Forces: Implementation and Applications to Atomic Nuclei and Dense Matter
Authors:
Kai Hebeler
Abstract:
Recent advances in nuclear structure theory have significantly enlarged the accessible part of the nuclear landscape via ab initio many-body calculations. These developments open new ways for microscopic studies of light, medium-mass and heavy nuclei as well as nuclear matter and represent an important step toward a systematic and comprehensive understanding of atomic nuclei across the nuclear cha…
▽ More
Recent advances in nuclear structure theory have significantly enlarged the accessible part of the nuclear landscape via ab initio many-body calculations. These developments open new ways for microscopic studies of light, medium-mass and heavy nuclei as well as nuclear matter and represent an important step toward a systematic and comprehensive understanding of atomic nuclei across the nuclear chart. While remarkable agreement has been found between different many-body methods for a given nuclear Hamiltonian, the comparison with experiment and the understanding of theoretical uncertainties are still important open questions. The observed discrepancies to experiment indicate deficiencies in presently used nuclear interactions and operators. Chiral effective field theory (EFT) allows to systematically derive contributions to nucleon-nucleon (NN), three-nucleon (3N) and higher-body interactions including estimates of theoretical uncertainties. While the treatment of NN interactions in many-body calculations is well established, the calculation of 3N interactions and their incorporation in ab initio frameworks is still a frontier. This work reviews in detail recent and current developments on the derivation and implementation of improved 3N interactions and provides a comprehensive introduction to fundamental methods for their practical calculation and representation. We further give an overview of novel and established methods that facilitate the inclusion and treatment of 3N interactions in ab initio nuclear structure frameworks and present a selection of the latest calculations of atomic nuclei as well as nuclear matter based on state-of-the-art nuclear NN and 3N interactions derived within chiral EFT. Finally, we discuss ongoing efforts, open questions and future directions.
△ Less
Submitted 6 January, 2021; v1 submitted 21 February, 2020;
originally announced February 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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
Equation of state and neutron star properties constrained by nuclear physics and observation
Authors:
K. Hebeler,
J. M. Lattimer,
C. J. Pethick,
A. Schwenk
Abstract:
Microscopic calculations of neutron matter based on nuclear interactions derived from chiral effective field theory, combined with the recent observation of a 1.97 +- 0.04 M_sun neutron star, constrain the equation of state of neutron-rich matter at sub- and supranuclear densities. We discuss in detail the allowed equations of state and the impact of our results on the structure of neutron stars,…
▽ More
Microscopic calculations of neutron matter based on nuclear interactions derived from chiral effective field theory, combined with the recent observation of a 1.97 +- 0.04 M_sun neutron star, constrain the equation of state of neutron-rich matter at sub- and supranuclear densities. We discuss in detail the allowed equations of state and the impact of our results on the structure of neutron stars, the crust-core transition density, and the nuclear symmetry energy. In particular, we show that the predicted range for neutron star radii is robust. For use in astrophysical simulations, we provide detailed numerical tables for a representative set of equations of state consistent with these constraints.
△ Less
Submitted 22 July, 2013; v1 submitted 19 March, 2013;
originally announced March 2013.
-
New applications of Renormalization Group methods to nuclear matter
Authors:
Kai Hebeler
Abstract:
We give an overview of recent results for the nuclear equation of state and properties of neutron stars based on microscopic two- and three-nucleon interactions derived within chiral effective field theory (EFT). It is demonstrated that the application of Renormalization Group (RG) transformations allows efficient, controlled and simplified calculations of nuclear matter.
We give an overview of recent results for the nuclear equation of state and properties of neutron stars based on microscopic two- and three-nucleon interactions derived within chiral effective field theory (EFT). It is demonstrated that the application of Renormalization Group (RG) transformations allows efficient, controlled and simplified calculations of nuclear matter.
△ Less
Submitted 21 August, 2012;
originally announced August 2012.
-
Neutron matter at next-to-next-to-next-to-leading order in chiral effective field theory
Authors:
I. Tews,
T. Krueger,
K. Hebeler,
A. Schwenk
Abstract:
Neutron matter presents a unique system for chiral effective field theory (EFT), because all many-body forces among neutrons are predicted to next-to-next-to-next-to-leading order (N3LO). We present the first complete N3LO calculation of the neutron matter energy. This includes the subleading three-nucleon (3N) forces for the first time and all leading four-nucleon (4N) forces. We find relatively…
▽ More
Neutron matter presents a unique system for chiral effective field theory (EFT), because all many-body forces among neutrons are predicted to next-to-next-to-next-to-leading order (N3LO). We present the first complete N3LO calculation of the neutron matter energy. This includes the subleading three-nucleon (3N) forces for the first time and all leading four-nucleon (4N) forces. We find relatively large contributions from N3LO 3N forces. Our results provide constraints for neutron-rich matter in astrophysics with controlled theoretical uncertainties.
△ Less
Submitted 4 February, 2013; v1 submitted 31 May, 2012;
originally announced June 2012.
-
Equation-of-state dependence of the gravitational-wave signal from the ring-down phase of neutron-star mergers
Authors:
A. Bauswein,
H. -T. Janka,
K. Hebeler,
A. Schwenk
Abstract:
Neutron-star (NS) merger simulations are conducted for 38 representative microphysical descriptions of high-density matter in order to explore the equation-of-state dependence of the postmerger ring-down phase. The formation of a deformed, oscillating, differentially rotating very massive NS is the typical outcome of the coalescence of two stars with 1.35 $M_{\odot}$ for most candidate EoSs. The o…
▽ More
Neutron-star (NS) merger simulations are conducted for 38 representative microphysical descriptions of high-density matter in order to explore the equation-of-state dependence of the postmerger ring-down phase. The formation of a deformed, oscillating, differentially rotating very massive NS is the typical outcome of the coalescence of two stars with 1.35 $M_{\odot}$ for most candidate EoSs. The oscillations of this object imprint a pronounced peak in the gravitational-wave (GW) spectra, which is used to characterize the emission for a given model. The peak frequency of this postmerger GW signal correlates very well with the radii of nonrotating NSs, and thus allows to constrain the high-density EoS by a GW detection. In the case of 1.35-1.35 $M_{\odot}$ mergers the peak frequency scales particularly well with the radius of a NS with 1.6 $M_{\odot}$, where the maximum deviation from this correlation is only 60 meters for fully microphysical EoSs which are compatible with NS observations. Combined with the uncertainty in the determination of the peak frequency it appears likely that a GW detection can measure the radius of a 1.6 $M_{\odot}$ NS with an accuracy of about 100 to 200 meters. We also uncover relations of the peak frequency with the radii of nonrotating NSs with 1.35 $M_{\odot}$ or 1.8 $M_{\odot}$, with the radius or the central energy density of the maximum-mass Tolman-Oppenheimer-Volkoff configuration, and with the pressure or sound speed at a fiducial rest-mass density of about twice nuclear saturation density. Furthermore, it is found that a determination of the dominant postmerger GW frequency can provide an upper limit for the maximum mass of nonrotating NSs. The prospects for a detection of the postmerger GW signal and a determination of the dominant GW frequency are estimated to be in the range of 0.015 to 1.2 events per year with the upcoming Advanced LIGO detector.
△ Less
Submitted 18 July, 2012; v1 submitted 9 April, 2012;
originally announced April 2012.
-
Constraints on the symmetry energy and neutron skins from experiments and theory
Authors:
M. B. Tsang,
J. R. Stone,
F. Camera,
P. Danielewicz,
S. Gandolfi,
K. Hebeler,
C. J. Horowitz,
Jenny Lee,
W. G. Lynch,
Z. Kohley,
R. Lemmon,
P. Moller,
T. Murakami,
S. Riordan,
X. Roca-Maza,
F. Sammarruca,
A. W. Steiner,
I. Vidaña,
S. J. Yennello
Abstract:
The symmetry energy contribution to the nuclear Equation of State (EoS) impacts various phenomena in nuclear astrophysics, nuclear structure, and nuclear reactions. Its determination is a key objective of contemporary nuclear physics with consequences for the understanding of dense matter within neutron stars. We examine the results of laboratory experiments that have provided initial constraints…
▽ More
The symmetry energy contribution to the nuclear Equation of State (EoS) impacts various phenomena in nuclear astrophysics, nuclear structure, and nuclear reactions. Its determination is a key objective of contemporary nuclear physics with consequences for the understanding of dense matter within neutron stars. We examine the results of laboratory experiments that have provided initial constraints on the nuclear symmetry energy and its density dependence at and somewhat below normal nuclear matter density. Some of these constraints have been derived from properties of nuclei. Others have been derived from the nuclear response to electroweak and hadronic probes. We also examine the most frequently used theoretical models that predict the symmetry energy and its slope. By comparing existing constraints on the symmetry pressure to theories, we demonstrate how the contribution of the three-body force, an essential ingredient in neutron matter models, can be determined.
△ Less
Submitted 2 April, 2012;
originally announced April 2012.
-
Constraints on neutron star radii based on chiral effective field theory interactions
Authors:
K. Hebeler,
J. M. Lattimer,
C. J. Pethick,
A. Schwenk
Abstract:
We show that microscopic calculations based on chiral effective field theory interactions constrain the properties of neutron-rich matter below nuclear densities to a much higher degree than is reflected in commonly used equations of state. Combined with observed neutron star masses, our results lead to a radius R = 9.7 - 13.9 km for a 1.4 M_{solar} star, where the theoretical range is due, in abo…
▽ More
We show that microscopic calculations based on chiral effective field theory interactions constrain the properties of neutron-rich matter below nuclear densities to a much higher degree than is reflected in commonly used equations of state. Combined with observed neutron star masses, our results lead to a radius R = 9.7 - 13.9 km for a 1.4 M_{solar} star, where the theoretical range is due, in about equal amounts, to uncertainties in many-body forces and to the extrapolation to high densities.
△ Less
Submitted 26 October, 2010; v1 submitted 10 July, 2010;
originally announced July 2010.
-
Chiral three-nucleon forces and neutron matter
Authors:
K. Hebeler,
A. Schwenk
Abstract:
We calculate the properties of neutron matter and highlight the physics of chiral three-nucleon forces. For neutrons, only the long-range 2 pi-exchange interactions of the leading chiral three-nucleon forces contribute, and we derive density-dependent two-body interactions by summing the third particle over occupied states in the Fermi sea. Our results for the energy suggest that neutron matter is…
▽ More
We calculate the properties of neutron matter and highlight the physics of chiral three-nucleon forces. For neutrons, only the long-range 2 pi-exchange interactions of the leading chiral three-nucleon forces contribute, and we derive density-dependent two-body interactions by summing the third particle over occupied states in the Fermi sea. Our results for the energy suggest that neutron matter is perturbative at nuclear densities. We study in detail the theoretical uncertainties of the neutron matter energy, provide constraints for the symmetry energy and its density dependence, and explore the impact of chiral three-nucleon forces on the S-wave superfluid pairing gap.
△ Less
Submitted 23 June, 2010; v1 submitted 2 November, 2009;
originally announced November 2009.
-
Dependence of the BCS 1S0 superfluid pairing gap on nuclear interactions
Authors:
K. Hebeler,
A. Schwenk,
B. Friman
Abstract:
We study in detail the dependence of the 1S0 superfluid pairing gap on nuclear interactions and on charge-independence breaking at the BCS level. Starting from chiral effective-field theory and conventional nucleon-nucleon (NN) interactions, we use the renormalization group to generate low-momentum interactions V_{low k} with sharp and smooth regulators. The resulting BCS gaps are well constrain…
▽ More
We study in detail the dependence of the 1S0 superfluid pairing gap on nuclear interactions and on charge-independence breaking at the BCS level. Starting from chiral effective-field theory and conventional nucleon-nucleon (NN) interactions, we use the renormalization group to generate low-momentum interactions V_{low k} with sharp and smooth regulators. The resulting BCS gaps are well constrained by the NN scattering phase shifts, and the cutoff dependence is very weak for sharp or sufficiently narrow smooth regulators with cutoffs Lambda > 1.6 fm^{-1}. It is therefore likely that the effect of three-nucleon interactions on 1S0 superfluidity is small at the BCS level. The charge dependence of nuclear interactions has a 10% effect on the pairing gap.
△ Less
Submitted 5 March, 2007; v1 submitted 8 November, 2006;
originally announced November 2006.