2021
DOI: 10.48550/arxiv.2111.14858
|View full text |Cite
Preprint
|
Sign up to set email alerts
|

A New Moment-Based General-Relativistic Neutrino-Radiation Transport Code: Methods and First Applications to Neutron Star Mergers

David Radice,
Sebastiano Bernuzzi,
Albino Perego
et al.

Abstract: We present a new moment-based neutrino transport code for neutron star merger simulations in general relativity. In the merger context, ours is the first code to include Doppler effects at all orders in 𝜐/𝑐, retaining all nonlinear neutrino-matter coupling terms. The code is validated with a stringent series of tests. We show that the inclusion of full neutrino-matter coupling terms is necessary to correctly capture the trapping of neutrinos in relativistically moving media, such as in differentially rotatin… Show more

Help me understand this report
View published versions

Search citation statements

Order By: Relevance

Paper Sections

Select...
3
1
1

Citation Types

2
5
0

Year Published

2022
2022
2023
2023

Publication Types

Select...
4
2

Relationship

0
6

Authors

Journals

citations
Cited by 6 publications
(7 citation statements)
references
References 123 publications
(209 reference statements)
2
5
0
Order By: Relevance
“…Prior research into the effectiveness of neutrino schemes have shown that the differences between twomoment (M1) and Monte Carlo (MC) schemes can result in a ∼ 20% uncertainty in neutrino luminosity, translating to a difference of 10% in outflow electron fraction (Foucart et al 2020). Comparison between M1 and the leakage scheme of Ardevol-Pulpillo et al (2019) (which our scheme is based on, see Appendix C) shows a further 10% uncertainty in neutrino luminosities, and a comparison of leakage+M0 to M1 schemes shows that leakage schemes tend to decrease the average Ye, but with a minimal effect on nucleosynthetic yields (Radice et al 2021). The exclusion of relativistic effects implies that we cannot accurately model jet formation, but the effects of these approximations on mass ejection and composition are likely minimal, since the relevant processes operate far from the BH.…”
Section: Summary and Discussionmentioning
confidence: 71%
See 1 more Smart Citation
“…Prior research into the effectiveness of neutrino schemes have shown that the differences between twomoment (M1) and Monte Carlo (MC) schemes can result in a ∼ 20% uncertainty in neutrino luminosity, translating to a difference of 10% in outflow electron fraction (Foucart et al 2020). Comparison between M1 and the leakage scheme of Ardevol-Pulpillo et al (2019) (which our scheme is based on, see Appendix C) shows a further 10% uncertainty in neutrino luminosities, and a comparison of leakage+M0 to M1 schemes shows that leakage schemes tend to decrease the average Ye, but with a minimal effect on nucleosynthetic yields (Radice et al 2021). The exclusion of relativistic effects implies that we cannot accurately model jet formation, but the effects of these approximations on mass ejection and composition are likely minimal, since the relevant processes operate far from the BH.…”
Section: Summary and Discussionmentioning
confidence: 71%
“…Leakage schemes have been shown to capture the dominant effects of neutrinos in post-merger tori around compact objects, especially for BH disks, for which they are subdominant energy sources (Foucart et al 2019;Fernández & Metzger 2013;Siegel & Metzger 2018;Fernández et al 2019a). However, significant differences appear when compared quantitatively to more advanced Monte-Carlo or two-moment (M1) schemes (Richers et al 2015;Foucart et al 2015;Perego et al 2016;Ardevol-Pulpillo et al 2019;Radice et al 2021). For this reason it is necessary to make impovements to the previous leakage-scheme implemented in FLASH (Fernández & Metzger 2013;Metzger & Fernández 2014), while retaining computational efficiency.…”
Section: Appendix C: Neutrino Leakage Schemementioning
confidence: 99%
“…multi-scale and multi-physics systems that combine disparate time and spatial scales, and which demand the use of subgrid-scale precision to accurately resolve the evolution of physical fields. This is a well known problem in multiple disciplines, including general relativistic simulations [3][4][5], weather forecasting [6], ab initio density functional theory simulations [7], among many other computational grand challenges.…”
Section: Introductionmentioning
confidence: 99%
“…At present, the Riemann solver and constrained transport scheme implemented in existing numerical relativity magnetohydrodynamics codes, e.g. [55,[60][61][62][63][64][65][66][67], are based on the HLLE solver [4,59,68]. (An exception is the SpECTRE [69] code, which is based on the discontinuous Galerkin method.)…”
Section: Introductionmentioning
confidence: 99%