Source author record

Martin Obergaulinger

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

10works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

10 published item(s)

preprint2022arXiv

R-process viable outflows are suppressed in global alpha-viscosity models of collapsar disks

Collapsar disks have been proposed to be rich factories of heavy elements, but the major question of whether their outflows are neutron-rich, and could therefore represent significant sites of the rapid neutron-capture (r-) process, or dominated by iron-group elements remains unresolved. We present the first global models of collapsars that start from a stellar progenitor and self-consistently describe the evolution of the disk, its composition, and its outflows in response to the imploding stellar mantle, using energy-dependent M1 neutrino transport and an alpha-viscosity to approximate turbulent angular-momentum transport. We find that a neutron-rich, neutrino-dominated accretion flow (NDAF) is established only marginally--either for short times or relatively low viscosities--because the disk tends to disintegrate into an advective disk (ADAF) already at relatively high mass-accretion rates, launching powerful outflows but preventing it from developing a hot, dense, and therefore neutron-rich core. Viscous outflows disrupt the star within ~100s with explosion energies close to that of hypernovae. If viscosity is neglected, a stable NDAF with disk mass of about 1Msun is formed but is unable to release neutron-rich ejecta, while it produces a relatively mild explosion powered by a neutrino-driven wind blown off its surface. With ejecta electron fractions close to 0.5, all models presumably produce large amounts of Ni56. Our results suggest that collapsar models based on the alpha-viscosity are inefficient r-process sites and that genuinely magnetohydrodynamic effects may be required to generate neutron-rich outflows. A relatively weak effective viscosity generated by magnetohydrodynamic turbulence would improve the prospects for obtaining neutron-rich ejecta.

preprint2021arXiv

Nucleosynthesis in magneto-rotational supernovae

We present the nucleosynthesis of magneto-rotational supernovae (MR-SNe) including neutrino-driven and magneto-rotational-driven ejecta based, for the first time, on 2D simulations with accurate neutrino transport. The models analysed here have different rotation and magnetic fields, allowing us to explore the impact of these two key ingredients. The accurate neutrino transport of the simulations is critical to analyse the slightly neutron-rich and proton-rich ejecta that are similar to the, also neutrino-driven, ejecta in standard supernovae. In the model with strong magnetic field, the r-process produces heavy elements up to the third r-process peak ($A\sim 195$), in agreement with previous works. This model presents a jet-like explosion with proton-rich jets surrounded by neutron-rich material where the r-process occurs. We have estimated a lower limit for $^{56}$Ni of $2.5\times10^{-2} M_\odot$, which is still well below the expected hypernova value. Longer simulations including the accretion disc evolution are required to get a final prediction. In addition, we have found that the late evolution is critical in a model with weak magnetic field in which late-ejected neutron-rich matter produces elements up to the second r-process peak. Even if we cannot yet provide conclusions for hypernova nucleosynthesis, our results agree with observations of old stars and radioactive isotopes in supernova remnants. This makes MR-SNe a good additional scenario to neutron star mergers for the synthesis of heavy elements and brings us closer to understand their origin and the role of MR-SNe in the early Galaxy nucleosynthesis.

preprint2016arXiv

Neutron-star merger ejecta as obstacles to neutrino-powered jets of gamma-ray bursts

We present the first special relativistic, axisymmetric hydrodynamic simulations of black hole-torus systems (approximating general relativistic gravity) as remnants of binary-neutron star (NS-NS) and neutron star-black hole (NS-BH) mergers, in which the viscously driven evolution of the accretion torus is followed with self-consistent energy-dependent neutrino transport and the interaction with the cloud of dynamical ejecta expelled during the NS-NS merging is taken into account. The modeled torus masses, BH masses and spins, and the ejecta masses, velocities, and spatial distributions are adopted from relativistic merger simulations. We find that energy deposition by neutrino annihilation can accelerate outflows with initially high Lorentz factors along polar low-density funnels, but only in mergers with extremely low baryon pollution in the polar regions. NS-BH mergers, where polar mass ejection during the merging phase is absent, provide sufficiently baryon-poor environments to enable neutrino-powered, ultrarelativistic jets with terminal Lorentz factors above 100 and considerable dynamical collimation, favoring short gamma-ray bursts (sGRBs), although their typical energies and durations might be too small to explain the majority of events. In the case of NS-NS mergers, however, neutrino emission of the accreting and viscously spreading torus is too short and too weak to yield enough energy for the outflows to break out from the surrounding ejecta shell as highly relativistic jets. We conclude that neutrino annihilation alone cannot power sGRBs from NS-NS mergers.

preprint2016arXiv

On the maximum magnetic field amplification by the magnetorotational instability in core-collapse supernovae

Whether the magnetorotational instability (MRI) can amplify initially weak magnetic fields to dynamically relevant strengths in core collapse supernovae is still a matter of active scientific debate. Recent numerical studies have shown that the first phase of MRI growth dominated by channel flows is terminated by parasitic instabilities of the Kelvin-Helmholtz type that disrupt MRI channel flows and quench further magnetic field growth. However, it remains to be prop- erly assessed by what factor the initial magnetic field can be amplified and how it depends on the initial field strength and the amplitude of the perturbations. Different termination criteria leading to different estimates of the amplification factor were proposed within the parasitic model. To determine the amplification factor and test which criterion is a better predictor of the MRI termination, we perform three-dimensional shearing-disc and shearing-box simula- tions of a region close to the surface of a differentially rotating proto-neutron star in non-ideal MHD with two different numerical codes. We find that independently of the initial magnetic field strength, the MRI channel modes can amplify the magnetic field by, at most, a factor of 100. Under the conditions found in proto-neutron stars a more realistic value for the mag- netic field amplification is of the order of 10. This severely limits the role of the MRI channel modes as an agent amplifying the magnetic field in proto-neutron stars starting from small seed fields. A further amplification should therefore rely on other physical processes, such as for example an MRI-driven turbulent dynamo.

preprint2016arXiv

Termination of the magnetorotational instability via parasitic instabilities in core-collapse supernovae

The magnetorotational instability (MRI) can be a powerful mechanism amplifying the magnetic field in core-collapse supernovae. Whether initially weak magnetic fields can be amplified by this instability to dynamically relevant strengths is still a matter of debate. One of the main uncertainties concerns the process that terminates the growth of the instability. Parasitic instabilities of both Kelvin-Helmholtz and tearing-mode type have been suggested to play a crucial role in this process, disrupting MRI channel flows and quenching magnetic field amplification. We perform two-dimensional and three-dimensional sheering-disc simulations of a differentially rotating protoneutron star layer in non-ideal magnetohydrodynamics with unprecedented high numerical accuracy, finding that Kelvin-Helmholtz parasitic modes dominate tearing modes in the regime of large hydrodynamic and magnetic Reynolds numbers, as encountered close to the surface of protoneutron stars. They also determine the maximum magnetic field stress achievable during the exponential growth of the MRI. Our results are consistent with the theory of parasitic instabilities based on a local stability analysis. To simulate the Kelvin-Helmholtz instabilities properly, a very high numerical resolution is necessary. Using ninth-order spatial reconstruction schemes, we find that at least eight grid zones per MRI channel are necessary to simulate the growth phase of the MRI and reach an accuracy of ~10 per cent in the growth rate, while more than ~60 zones per channel are required to achieve convergent results for the value of the magnetic stress at MRI termination.

preprint2015arXiv

A new multidimensional, energy-dependent two-moment transport code for neutrino-hydrodynamics

We present the new code ALCAR developed to model multidimensional, multi energy-group neutrino transport in the context of supernovae and neutron-star mergers. The algorithm solves the evolution equations of the 0th- and 1st-order angular moments of the specific intensity, supplemented by an algebraic relation for the 2nd-moment tensor to close the system. The scheme takes into account frame-dependent effects of order O(v/c) as well as the most important types of neutrino interactions. The transport scheme is significantly more efficient than a multidimensional solver of the Boltzmann equation, while it is more accurate and consistent than the flux-limited diffusion method. The finite-volume discretization of the essentially hyperbolic system of moment equations employs methods well-known from hydrodynamics. For the time integration of the potentially stiff moment equations we employ a scheme in which only the local source terms are treated implicitly, while the advection terms are kept explicit, thereby allowing for an efficient computational parallelization of the algorithm. We investigate various problem setups in one and two dimensions to verify the implementation and to test the quality of the algebraic closure scheme. In our most detailed test, we compare a fully dynamic, one-dimensional core-collapse simulation with two published calculations performed with well-known Boltzmann-type neutrino-hydrodynamics codes and we find very satisfactory agreement.

preprint2014arXiv

Magnetic field amplification and magnetically supported explosions of collapsing, non-rotating stellar cores

We study the amplification of magnetic fields in the collapse and the post-bounce evolution of the core of a non-rotating star of 15 solar masses in axisymmetry. To this end, we solve the coupled equations of magnetohydrodynamics and neutrino transport in the two-moment approximation. The pre-collapse magnetic field is strongly amplified by compression in the infall. Initial fields of the order of 1010 G translate into proto-neutron star fields similar to the ones observed in pulsars, while stronger initial fields yield magnetar-like final field strengths. After core bounce, the field is advected through the hydrodynamically unstable neutrino-heating layer, where non-radial flows due to convection and the standing accretion shock instability amplify the field further. Consequently, the resulting amplification factor of order five is the result of the number of small-eddy turnovers taking place within the time scale of advection through the post-shock layer. Due to this limit, most of our models do not reach equipartition between kinetic and magnetic energy and, consequently, evolve similarly to the non-magnetic case, exploding after about 800 ms when a single or few high-entropy bubbles persist over several dynamical time scales. In the model with the strongest initial field we studied, \$10^{12}\$ G, for which equipartition between flow and field is achieved, the magnetic tension favours a much earlier development of such long-lived high-entropy bubbles and enforces a fairly ordered large-scale flow pattern. Consequently, this model, after exhibiting very regular shock oscillations, explodes much earlier than non-magnetic ones.

preprint2013arXiv

Gravitational wave signatures in black-hole-forming core collapse

We present numerical simulations in general relativity of collapsing stellar cores. Our initial model consists of a low metallicity rapidly-rotating progenitor which is evolved in axisymmetry with the latest version of our general relativistic code CoCoNuT, which allows for black hole formation and includes the effects of a microphysical equation of state (LS220) and a neutrino leakage scheme to account for radiative losses. The motivation of our study is to analyze in detail the emission of gravitational waves in the collapsar scenario of long gamma-ray bursts. Our simulations show that the phase during which the proto-neutron star (PNS) survives before ultimately collapsing to a black hole is particularly optimal for gravitational wave emission. The high-amplitude waves last for several seconds and show a remarkable quasi-periodicity associated with the violent PNS dynamics, namely during the episodes of convection and the subsequent nonlinear development of the standing-accretion shock instability (SASI). By analyzing the spectrogram of our simulations we are able to identify the frequencies associated with the presence of g-modes and with the SASI motions at the PNS surface. We note that the gravitational waves emitted reach large enough amplitudes to be detected with third-generation detectors as the Einstein Telescope within a Virgo cluster volume at rates ~< 0.1 /y.

preprint2012arXiv

Core-Collapse Supernovae: Reflections and Directions

Core-collapse supernovae are among the most fascinating phenomena in astrophysics and provide a formidable challenge for theoretical investigation. They mark the spectacular end of the lives of massive stars and, in an explosive eruption, release as much energy as the sun produces during its whole life. A better understanding of the astrophysical role of supernovae as birth sites of neutron stars, black holes, and heavy chemical elements, and more reliable predictions of the observable signals from stellar death events are tightly linked to the solution of the long-standing puzzle how collapsing stars achieve to explode. In this article our current knowledge of the processes that contribute to the success of the explosion mechanism are concisely reviewed. After a short overview of the sequence of stages of stellar core-collapse events, the general properties of the progenitor-dependent neutrino emission will be briefly described. Applying sophisticated neutrino transport in axisymmetric (2D) simulations with general relativity as well as in simulations with an approximate treatment of relativistic effects, we could find successful neutrino-driven explosions for a growing set of progenitor stars. First results of three-dimensional (3D) models have been obtained, and magnetohydrodynamic simulations demonstrate that strong initial magnetic fields in the pre-collapse core can foster the onset of neutrino-powered supernova explosions even in nonrotating stars. These results are discussed in the context of the present controversy about the value of 2D simulations for exploring the supernova mechanism in realistic 3D environments, and they are interpreted against the background of the current disagreement on the question whether the standing accretion shock instability (SASI) or neutrino-driven convection is the crucial agency that supports the onset of the explosion.

preprint2011arXiv

Magnetic field amplification in collapsing, non-rotating stellar cores

Context. The influence of magnetic fields on stellar core collapse and explosion is not well explored. It depends on the possibility to amplify the pre-collapse fields. Without rotation this can happen by compression, convection, the standing accretion shock instability, and the accumulation and growth of Alfven waves in the accretion flows. Aims. We investigate such amplification mechanisms of the magnetic field during the collapse and post-bounce evolution of the core of a non-rotating 15 solar mass star with varied initial field strengths, taking into account the microphysical equation of state and neutrino physics crucial for supernova cores. Methods. We perform simulations of ideal MHD with neutrino transport in axisymmetry. The transport of electron neutrinos and antineutrinos is treated with a new scheme that solves the energy-dependent set of radiation energy and momentum equations in 2d by using an analytic closure relation. Results. The magnetic field undergoes amplification by turbulent flows. We find indications for amplification by interacting waves in accretion streams. The fields can reach up to equipartition with the velocity field. Very high magnetic field strengths require very strong pre-collapse fields and are able to shape the post-bounce flow, leading to a pattern dominated by low-order multipoles. Such models are closest to a successful explosion. Conclusions. Magnetic fields can build up to interesting strengths even in non-rotating collapsing stellar cores. Starting with fields in the pre-collapse core as predicted by present stellar evolution models, typical neutron star fields emerge, whereas stronger progenitor fields lead to fields of magnetar strength. Only in the latter case the fields have dynamical effects on the flows in the supernova core. However, in none of our 2D simulations, we find an explosion until 500 ms post-bounce.