Martin Obergaulinger

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, \

NEUTRON-STAR MERGER EJECTA AS OBSTACLES TO NEUTRINO-POWERED JETS OF GAMMA-RAY BURSTS

10^{12}\

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

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.

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

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.

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

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.

Hydrodynamic simulations of the interaction of supernova shock waves with a clumpy environment: the case of the RX J0852.0−4622 (Vela Jr) supernova remnant

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 di erentially rotating proto-neutron 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 proto-neutron 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 9th order spatial reconstruction schemes, we find that at least 8 grid zones per MRI channel are necessary to simulate the growth phase of the MRI and reach an accuracy of 10% 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.

Protomagnetar and black hole formation in high-mass stars

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.

Estimation of the mechanical properties of the eye through the study of its vibrational modes

Observations in all electromagnetic bands show that many supernova remnants (SNRs) have a very aspherical shape. This can be the result of asymmetries in the supernova explosion or a clumpy circumstellar medium. We study the generation of inhomogeneities and the mixing of elements arising from these two sources in multidimensional hydrodynamic simulations of the propagation of a supernova blast wave into a cloudy environment. We model a specific SNR, Vela Jr (RX J0852.0-4622). By comparing our results with recent observations, we can constrain the properties of the explosion. We find that a very energetic explosion of several 10^{51} erg occurring roughly about 800 years ago is consistent with the shape and emission of the SNR, as well as a supernova with an energy closer to the canonical value of 10^{51} erg a few thousand years ago.

On the measurements of numerical viscosity and resistivity in Eulerian MHD codes

Using axisymmetric simulations coupling special relativistic MHD, an approximate post-Newtonian gravitational potential and two-moment neutrino transport, we show different paths for the formation of either protomagnetars or stellar mass black holes. The fraction of prototypical stellar cores which should result in collapsars depends on a combination of several factors, among which the structure of the progenitor star and the profile of specific angular momentum are probably the foremost. Along with the implosion of the stellar core, we also obtain supernova-like explosions driven by neutrino heating and hydrodynamic instabilities or by magneto-rotational effects in cores of high-mass stars. In the latter case, highly collimated, mildly relativistic outflows are generated. We find that after a rather long post-collapse phase (lasting >~ 1 sec) black holes may form in cases both of successful and failed supernovalike explosions. A basic trend is that cores with a specific angular momentum smaller than that obtained by standard, one-dimensional stellar evolution calculations form black holes (and eventually collapsars). Complementary, protomagnetars result from stellar cores with the standard distribution of specific angular momentum obtained from prototypical stellar evolution calculations including magnetic torques and moderate to large mass loss rates.

Core-collapse supernova simulations in one and two dimensions: comparison of codes and approximations

Measuring the eye’s mechanical properties in vivo and with minimally invasive techniques can be the key for individualized solutions to a number of eye pathologies. The development of such techniques largely relies on a computational modelling of the eyeball and, it optimally requires the synergic interplay between experimentation and numerical simulation. In Astrophysics and Geophysics the remote measurement of structural properties of the systems of their realm is performed on the basis of (helio-)seismic techniques. As a biomechanical system, the eyeball possesses normal vibrational modes encompassing rich information about its structure and mechanical properties. However, the integral analysis of the eyeball vibrational modes has not been performed yet. Here we develop a new finite difference method to compute both the spheroidal and, specially, the toroidal eigenfrequencies of the human eye. Using this numerical model, we show that the vibrational eigenfrequencies of the human eye fall in the interval 100 Hz–10 MHz. We find that compressible vibrational modes may release a trace on high frequency changes of the intraocular pressure, while incompressible normal modes could be registered analyzing the scattering pattern that the motions of the vitreous humour leave on the retina. Existing contact lenses with embebed devices operating at high sampling frequency could be used to register the microfluctuations of the eyeball shape we obtain. We advance that an inverse problem to obtain the mechanical properties of a given eye (e.g., Young’s modulus, Poisson ratio) measuring its normal frequencies is doable. These measurements can be done using non-invasive techniques, opening very interesting perspectives to estimate the mechanical properties of eyes in vivo. Future research might relate various ocular pathologies with anomalies in measured vibrational frequencies of the eye.