BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Luke F. Roberts

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

ResearcherUnclaimed source record
astro-ph.HEgr-qcnucl-thastro-ph.COastro-ph.GAhep-phMachine Learningmath.NANumerical Analysisphysics.comp-ph

Catalog footprint

What is connected

11works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

11 published item(s)

preprint2022arXiv

HARM3D+NUC: A new method for simulating the post-merger phase of binary neutron star mergers with GRMHD, tabulated EOS and neutrino leakage

The first binary neutron star merger has already been detected in gravitational waves. The signal was accompanied by an electromagnetic counterpart including a kilonova component powered by the decay of radioactive nuclei, as well as a short $γ$-ray burst. In order to understand the radioactively-powered signal, it is necessary to simulate the outflows and their nucleosynthesis from the post-merger disk. Simulating the disk and predicting the composition of the outflows requires general relativistic magnetohydrodynamical (GRMHD) simulations that include a realistic, finite-temperature equation of state (EOS) and self-consistently calculating the impact of neutrinos. In this work, we detail the implementation of a finite-temperature EOS and the treatment of neutrinos in the GRMHD code HARM3D+NUC, based on HARM3D. We include formal tests of both the finite-temperature EOS and the neutrino leakage scheme. We further test the code by showing that, given conditions similar to those of published remnant disks following neutron star mergers, it reproduces both recombination of free nucleons to a neutron-rich composition and excitation of a thermal wind.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Machine learning moment closure models for the radiative transfer equation I: directly learning a gradient based closure

In this paper, we take a data-driven approach and apply machine learning to the moment closure problem for radiative transfer equation in slab geometry. Instead of learning the unclosed high order moment, we propose to directly learn the gradient of the high order moment using neural networks. This new approach is consistent with the exact closure we derive for the free streaming limit and also provides a natural output normalization. A variety of benchmark tests, including the variable scattering problem, the Gaussian source problem with both periodic and reflecting boundaries, and the two-material problem, show both good accuracy and generalizability of our machine learning closure model.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Wave heating from proto-neutron star convection and the core-collapse supernova explosion mechanism

Our understanding of the core-collapse supernova explosion mechanism is incomplete. While the favoured scenario is delayed revival of the stalled shock by neutrino heating, it is difficult to reliably compute explosion outcomes and energies, which depend sensitively on the complex radiation hydrodynamics of the post-shock region. The dynamics of the (non-)explosion depend sensitively on how energy is transported from inside and near the proto-neutron star (PNS) to material just behind the supernova shock. Although most of the PNS energy is lost in the form of neutrinos, hydrodynamic and hydromagnetic waves can also carry energy from the PNS to the shock. We show that gravity waves excited by core PNS convection can couple with outgoing acoustic waves that present an appreciable source of energy and pressure in the post-shock region. Using one-dimensional simulations, we estimate the gravity wave energy flux excited by PNS convection and the fraction of this energy transmitted upward to the post-shock region as acoustic waves. We find wave energy fluxes near $10^{51}\,\mathrm{erg}\,\mathrm{s}^{-1}$ are likely to persist for $\sim1\,\mathrm{s}$ post-bounce. The wave pressure on the shock may exceed $10\%$ of the thermal pressure, potentially contributing to shock revival and, subsequently, a successful and energetic explosion. We also discuss how future simulations can better capture the effects of waves, and more accurately quantify wave heating rates.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Exciting Prospects for Detecting Late-Time Neutrinos from Core-Collapse Supernovae

The importance of detecting neutrinos from a Milky Way core-collapse supernova is well known. An under-studied phase is proto-neutron star cooling. For SN 1987A, this seemingly began at about 2 s, and is thus probed by only 6 of the 19 events (and only the $\barν_e$ flavor) in the Kamiokande-II and IMB detectors. With the higher statistics expected for present and near-future detectors, it should be possible to measure detailed neutrino signals out to very late times. We present the first comprehensive study of neutrino detection during the proto-neutron star cooling phase, considering a variety of outcomes, using all flavors, and employing detailed detector physics. For our nominal model, the event yields (at 10 kpc) after 10 s -- the approximate duration of the SN 1987A signal -- far exceed the entire SN 1987A yield, with $\simeq$250 $\barν_e$ events (to 50 s) in Super-Kamiokande, $\simeq$110 $ν_e$ events (to 40 s) in DUNE, and $\simeq$10 $ν_μ, ν_τ, \barν_μ, \barν_τ$ events (to 20 s) in JUNO. These data would allow unprecedented probes of the proto-neutron star, including the onset of neutrino transparency and hence its transition to a neutron star. If a black hole forms, even at very late times, this can be clearly identified. But will the detectors fulfill their potential for this perhaps once-ever opportunity for an all-flavor, high-statistics detection of a core collapse? Maybe. Further work is urgently needed, especially for DUNE to thoroughly investigate and improve its MeV capabilities.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dynamical Mass Ejection from Binary Neutron Star Mergers

We present fully general-relativistic simulations of binary neutron star mergers with a temperature and composition dependent nuclear equation of state. We study the dynamical mass ejection from both quasi-circular and dynamical-capture eccentric mergers. We systematically vary the level of our treatment of the microphysics to isolate the effects of neutrino cooling and heating and we compute the nucleosynthetic yields of the ejecta. We find that eccentric binaries can eject significantly more material than quasi-circular binaries and generate bright infrared and radio emission. In all our simulations the outflow is composed of a combination of tidally- and shock-driven ejecta, mostly distributed over a broad $\sim 60^\circ$ angle from the orbital plane, and, to a lesser extent, by thermally driven winds at high latitudes. Ejecta from eccentric mergers are typically more neutron rich than those of quasi-circular mergers. We find neutrino cooling and heating to affect, quantitatively and qualitatively, composition, morphology, and total mass of the outflows. This is also reflected in the infrared and radio signatures of the binary. The final nucleosynthetic yields of the ejecta are robust and insensitive to input physics or merger type in the regions of the second and third r-process peaks. The yields for elements on the first peak vary between our simulations, but none of our models is able to explain the Solar abundances of first-peak elements without invoking additional first-peak contributions from either neutrino and viscously-driven winds operating on longer timescales after the mergers, or from core-collapse supernovae.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

How loud are neutron star mergers?

We present results from the first large parameter study of neutron star mergers using fully general relativistic simulations with finite-temperature microphysical equations of state and neutrino cooling. We consider equal and unequal-mass binaries drawn from the galactic population and simulate each binary with three different equations of state. Our focus is on the emission of energy and angular momentum in gravitational waves in the postmerger phase. We find that the emitted gravitational-wave energy in the first $\sim$$10\,\mathrm{ms}$ of the life of the resulting hypermassive neutron star (HMNS) is about twice the energy emitted over the entire inspiral history of the binary. The total radiated energy per binary mass is comparable to or larger than that of nonspinning black hole inspiral-mergers. About $0.8-2.5\%$ of the binary mass-energy is emitted at kHz frequencies in the early HMNS evolution. We find a clear dependence of the postmerger GW emission on binary configuration and equation of state and show that it can be encoded as a broad function of the binary tidal coupling constant $κ^T_2$. Our results also demonstrate that the dimensionless spin of black holes resulting from subsequent HMNS collapse are limited to $\lesssim0.7-0.8$. This may significantly impact the neutrino pair annihilation mechanism for powering short gamma-ray bursts (sGRB).

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Compact Stellar Binary Assembly in the First Nuclear Star Clusters and r-Process Synthesis in the Early Universe

Investigations of elemental abundances in the ancient and most metal deficient stars are extremely important because they serve as tests of variable nucleosynthesis pathways and can provide critical inferences of the type of stars that lived and died before them. The presence of r-process elements in a handful of carbon-enhanced metal-poor (CEMP-r) stars, which are assumed to be closely connected to the chemical yield from the first stars, is hard to reconcile with standard neutron star mergers. Here we show that the production rate of dynamically assembled compact binaries in high-z nuclear star clusters can attain a sufficient high value to be a potential viable source of heavy r-process material in CEMP-r stars. The predicted frequency of such events in the early Galaxy, much lower than the frequency of Type II supernovae but with significantly higher mass ejected per event, can naturally lead to a high level of scatter of Eu as observed in CEMP-r stars.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A New Code for Proto-Neutron Star Evolution

A new code for following the evolution and emissions of proto-neutron stars during the first minute of their lives is developed and tested. The code is one dimensional, fully implicit, and general relativistic. Multi-group, multi-flavor neutrino transport is incorporated that makes use of variable Eddington factors obtained from a formal solution of the static general relativistic Boltzmann equation with linearized scattering terms. The timescales of neutrino emission and spectral evolution obtained using the new code are broadly consistent with previous results. Unlike other recent calculations, however, the new code predicts that the neutrino-driven wind will be characterized, at least for part of its existence, by a neutron excess. This change, potentially consequential for nucleosynthesis in the wind, is due to an improved treatment of the charged-current interactions of electron flavored neutrinos and anti-neutrinos with nucleons. A comparison is also made between the results obtained using either variable Eddington factors or simple equilibrium flux-limited diffusion. The latter approximation, which has been frequently used in previous studies of proto-neutron star cooling, accurately describes the total neutrino luminosities (to within 10%) for most of the evolution, until the proto-neutron star becomes optically thin.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Medium modification of the charged current neutrino opacity and its implications

Previous work on neutrino emission from proto-neutron stars which employed full solutions of the Boltzmann equation showed that the average energies of emitted electron neutrinos and antineutrinos are closer to one another than predicted by older, more approximate work. This in turn implied that the neutrino driven wind is proton rich during its entire life, precluding $r$-process nucleosynthesis and the synthesis of Sr, Y, and Zr. This work relied on charged current neutrino interaction rates that are appropriate for a free nucleon gas. Here, it is shown in detail that the inclusion of the nucleon potential energies and collisional broadening of the response significantly alters this conclusion. Iso-vector interactions, which give rise to the nuclear symmetry energy, produce a difference between the neutron and proton single-particle energies $ΔU=U_n-U_p$ and alter the kinematics of the charged current reactions. In neutron-rich matter, and for a given neutrino/antineutrino energy, the rate for $ν_e+n\rightarrow e^-+p$ is enhanced while $ \barν_e+p\rightarrow n+e^+$ is suppressed because the $Q$ value for these reactions is altered by $\pmΔU$, respectively. In the neutrino decoupling region, collisional broadening acts to enhance both $ν_e$ and $\barν_e$ cross-sections and RPA corrections decrease the $ν_e$ cross-section and increase the $\bar ν_e$ cross-section, but mean field shifts have a larger effect. Therefore, electron neutrinos decouple at lower temperature than when the nucleons are assumed to be free and have lower average energies. The change is large enough to allow for a reasonable period of time when the neutrino driven wind is predicted to be neutron rich. It is also shown that the electron fraction in the wind is influenced by the nuclear symmetry energy.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Electromagnetic Transients Powered by Nuclear Decay in the Tidal Tails of Coalescing Compact Binaries

The possibility that long tidal tails formed during compact object mergers may power optical transients through the decay of freshly synthesized r-process material is investigated. Precise modeling of the merger dynamics allows for a realistic determination of the thermodynamic conditions in the ejected debris. The results of hydrodynamic and full nuclear network calculations are combined to calculate the resultant r-process abundances and the heating of the material by their decays. The subsequent homologous structure is mapped into a radiative transfer code to synthesize emergent model light curves and determine how their properties (variability and color evolution) depend on the mass ratio and orientation of the merging binary. The radiation emanating from the ejected debris, though less spectacular than a typical supernova, should be observable in transient surveys and we estimate the associated detection rates. The case for (or against) compact object mergers as the progenitors of short gamma-ray bursts can be tested if such electromagnetic transients are detected (or not) in coincidence with some bursts, although they may be obscured by on-axis afterglows.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Proto-Neutron Star Cooling with Convection: The Effect of the Symmetry Energy

We model neutrino emission from a newly born neutron star subsequent to a supernova explosion to study its sensitivity to the equation of state, neutrino opacities, and convective instabilities at high baryon density. We find the time period and spatial extent over which convection operates is sensitive to the behavior of the nuclear symmetry energy at and above nuclear density. When convection ends within the proto-neutron star, there is a break in the predicted neutrino emission that may be clearly observable.

0 citations0 reviewsNo code linkedOpen access