Source author record

Bernhard Müller

Bernhard Müller 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

13works
9topics
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

13 published item(s)

preprint2022arXiv

Inferring Astrophysical Parameters of Core-Collapse Supernovae from their Gravitational-Wave Emission

Nearby core-collapse supernovae (CCSNe) are powerful multi-messenger sources for gravitational-wave, neutrino and electromagnetic telescopes as they emit gravitational waves in the ideal frequency band for ground based detectors. Once a CCSN gravitational-wave signal is detected, we will need to determine the parameters of the signal, and understand how those parameters relate to the source's explosion, progenitor and remnant properties. This is a challenge due to the stochastic nature of CCSN explosions, which is imprinted on their time series gravitational waveforms. In this paper, we perform Bayesian parameter estimation of CCSN signals using an asymmetric chirplet signal model to represent the dominant high-frequency mode observed in spectrograms of CCSN gravitational-wave signals. We use design sensitivity Advanced LIGO noise and CCSN waveforms from four different hydrodynamical supernova simulations with a range of different progenitor stars. We determine how well our model can reconstruct time-frequency images of the emission modes, and show how well we can determine parameters of the signal such as the frequency, amplitude, and duration. We show how the parameters of our signal model may allow us to place constraints on the proto-neutron star mass and radius, the turbulent kinetic energy onto the proto-neutron star, and the time of shock revival.

preprint2021arXiv

Explodability fluctuations of massive stellar cores enable asymmetric compact object mergers such as GW190814

The first three observing runs with Advanced LIGO and Virgo have resulted in the detection of binary black hole mergers (BBH) with highly unequal mass components, which are difficult to reconcile with standard formation paradigms. The most representative of these is GW190814, a highly asymmetric merger between a 23 M$_{\odot}$ black hole and a 2.6 M$_{\odot}$ compact object. Here, we explore recent results suggesting that a sizeable fraction of stars with pre-collapse carbon-oxygen core masses above 10 M$_{\odot}$, and extending up to at least 30 M$_{\odot}$, may produce objects inside the so-called lower mass gap that bridges the division between massive pulsars and BHs in Galactic X-ray binaries. We demonstrate that such an explosion landscape would naturally cause a fraction of massive binaries to produce GW190814-like systems instead of symmetric-mass BBHs. We present examples of specific evolutionary channels leading to the formation of GW190814 and GW200210, a 24+2.8 M$_{\odot}$ merger discovered during the O3b observing run. We estimate the merger-rate density of these events in our scenario to be $\mathcal{O}$(5%) of the total BBH merger rate. Finally, we discuss the broader implications of this formation channel for compact object populations, and its possible relevance to less asymmetric merger events such as GW200105 and GW200115

preprint2021arXiv

Gravitational Wave Signals from Two-Dimensional Core-Collapse Supernova Models with Rotation and Magnetic Fields

We investigate the impact of rotation and magnetic fields on the dynamics and gravitational wave emission in 2D core-collapse supernova simulations with neutrino transport. We simulate 17 different models of $15\,M_\odot$ and $39\,M_\odot$ progenitor stars with various initial rotation profiles and initial magnetic fields strengths up to $10^{12}\, \mathrm{G}$, assuming a dipolar field geometry in the progenitor. Strong magnetic fields generally prove conducive to shock revival, though this trend is not without exceptions. The impact of rotation on the post-bounce dynamics is more variegated, in line with previous studies. A significant impact on the time-frequency structure of the gravitational wave signal is found only for rapid rotation or strong initial fields. For rapid rotation, the angular momentum gradient at the proto-neutron star surface can appreciably affect the frequency of the dominant mode, so that known analytic relations for the high-frequency emission band no longer hold. In case of two magnetorotational explosion models, the deviation from these analytic relations is even more pronounced. One of the magnetorotational explosions has been evolved to more than half a second after the onset of the explosion and shows a subsidence of high-frequency emission at late times. Its most conspicuous gravitational wave signature is a high-amplitude tail signal. We also estimate the maximum detection distances for our waveforms. The magnetorotational models do not stick out for higher detectability during the post-bounce and explosion phase.

preprint2020arXiv

A 3D Simulation of a Neutrino-Driven Supernova Explosion Aided By Convection and Magnetic Fields

We study the impact of a small-scale dynamo in core-collapse supernovae using a 3D neutrino magnetohydrodynamics simulation of a $15 M_\odot$ progenitor. The weak seed field is amplified exponentially in the gain region once neutrino-driven convection develops, and remains dominated by small-scale structures. About $250\, \mathrm{ms}$ after bounce, the field energy in the gain region reaches $\mathord{\sim} 50\%$ of kinetic equipartition. This supports the development of a neutrino-driven explosion with modest global anisotropy, which does not occur in a corresponding model without magnetic fields. Our results suggest that magnetic fields may play a beneficial subsidiary role in neutrino-driven supernovae even without rapid progenitor rotation. Further investigation into the nature of magnetohydrodynamic turbulence in the supernova core is required.

preprint2020arXiv

Defining mass transfer in a capillary wave micro-bioreactor

For high-throughput cell culture and associated analytics, droplet-based cultivation systems open up the opportunities for parallelization and rapid data generation. In contrast to microfluidics with continuous flow, sessile droplet approaches enhance the flexibility for fluid manipulation with usually less operational effort. Generating biologically favorable conditions and promoting cell growth in a droplet, however, is particularly challenging due to mass transfer limitations, which has to be solved by implementing an effective mixing technique. Here, capillary waves induced by vertical oscillation are used to mix inside a sessile droplet micro-bioreactor (MBR) system avoiding additional moving parts inside the fluid. Depending on the excitation frequency, different patterns are formed on the oscillating liquid surface, which are described by a model of a vibrated sessile droplet. Analyzing mixing times and oxygen transport into the liquid, a strong dependency of mass transfer on the oscillation parameters, especially the excitation frequency, is demonstrated. Oscillations at distinct capillary wave resonant frequencies lead to rapid homogenization with mixing times of 2 s and volumetric liquid-phase mass transfer coefficients of more than 340 h-1. This shows that the mass transfer in a droplet MBR can be specifically controlled via capillary waves, what is subsequently demonstrated for cultivations of Escherichia coli BL21 cells. Therefore, the presented MBR in combination with vertical oscillation mixing for intensified mass transfer is a promising tool for highly parallel cultivation and data generation.

preprint2020arXiv

The Chemical Evolution of Iron-Peak Elements with Hypernovae

We calculate the mean evolution of the iron-peak abundance ratios [(Cr,Mn,Co,Zn)/Fe] in the Galaxy, using modern supernova and hypernova chemical yields and a Galactic Chemical Evolution code that assumes homogeneous chemical evolution. We investigate a range of hypernova occurrence rates and are able to produce a chemical composition that is a reasonable fit to the observed values in metal-poor stars. This requires a hypernova occurence rate that is large (50%) in the early Universe, decreasing throughout evolution to a value that is within present day observational constraints (>~ 1%). A large hypernova occurence rate is beneficial to matching the high [Zn/Fe] observed in the most metal-poor stars, although including hypernovae with progenitor mass >= 60 solar masses is detrimental to matching the observed [(Mn,Co)/Fe] evolution at low [Fe/H]. A significant contribution from HNe seems to be critical for producing supersolar [(Co,Zn)/Fe] at low metallicity, though more work will need to be done in order to match the most extreme values. We also emphasise the need to update models for the enrichment sources at higher metallicity, as the satisfactory recovery of the solar values of [(Cr,Mn,Co,Zn)/Fe] still presents a challenge.

preprint2020arXiv

Three-dimensional core-collapse supernova simulations of massive and rotating progenitors

We present three-dimensional simulations of the core-collapse of massive rotating and non-rotating progenitors performed with the general relativistic neutrino hydrodynamics code CoCoNuT-FMT and analyse their explosion properties and gravitational-wave signals. The progenitor models include Wolf-Rayet stars with initial helium star masses of $39\,M_{\odot}$ and $20\,M_{\odot}$, and an $18\,M_{\odot}$ red supergiant. The $39\,M_{\odot}$ model is a rapid rotator, whereas the two other progenitors are non-rotating. Both Wolf-Rayet models produce healthy neutrino-driven explosions, whereas the red supergiant model fails to explode. By the end of the simulations, the explosion energies have already reached $1.1\times 10^{51}\,\mathrm{erg}$ and $0.6\times 10^{51}\,\mathrm{erg}$ for the $39\,M_{\odot}$ and $20\,M_{\odot}$ model, respectively. The explosions produce neutron stars of relatively high mass, but with modest kicks. Due to the alignment of the bipolar explosion geometry with the rotation axis, there is a relatively small misalignment of $30^\circ$ between the spin and the kick in the $39\,M_{\odot}$ model. In terms of gravitational-wave signals, the massive and rapidly rotating $39\,M_{\odot}$ progenitor stands out by large gravitational-wave amplitudes that would make it detectable out to almost 2 Mpc by the Einstein Telescope. For this model, we find that rotation significantly changes the dependence of the characteristic gravitational-wave frequency of the f-mode on the proto-neutron star parameters compared to the non-rotating case. The other two progenitors have considerably smaller detection distances, despite significant low-frequency emission in the most sensitive frequency band of current gravitational-wave detectors due to the standing accretion shock instability in the $18\,M_{\odot}$ model.

preprint2019arXiv

Large-Scale Mixing in a Violent Oxygen-Neon Shell Merger Prior to a Core-Collapse Supernova

We present a seven-minute long $4π$-3D simulation of a shell merger event in a non-rotating $18.88\, M_\odot$ supernova progenitor before the onset of gravitational collapse. The key motivation is to capture the large-scale mixing and asymmetries in the wake of the shell merger before collapse using a self-consistent approach. The $4π$ geometry is crucial as it allows us to follow the growth and evolution of convective modes on the largest possible scales. We find significant differences between the kinematic, thermodynamic and chemical evolution of the 3D and the 1D model. The 3D model shows vigorous convection leading to more efficient mixing of nuclear species. In the 3D case the entire oxygen shell attains convective Mach numbers of $\mathord{\approx}\, 0.1$, whereas in the 1D model, the convective velocities are much lower and there is negligible overshooting across convective boundaries. In the 3D case, the convective eddies entrain nuclear species from the neon (and carbon) layers into the deeper part of the oxygen burning shell, where they burn and power a violent convection phase with outflows. This is a prototypical model of a convective-reactive system. Due to the strong convection and the resulting efficient mixing, the interface between the neon layer and the silicon-enriched oxygen layer disappears during the evolution, and silicon is mixed far out into merged oxygen/neon shell. Neon entrained inwards by convective downdrafts burns, resulting in lower neon mass in the 3D model compared to the 1D model at time of collapse. In addition, the 3D model develops remarkable large-scale, large-amplitude asymmetries, which may have important implications for the impending gravitational collapse and the subsequent explosion.

preprint2016arXiv

Progenitor-dependent Explosion Dynamics in Self-consistent, Axisymmetric Simulations of Neutrino-driven Core-collapse Supernovae

We present self-consistent, axisymmetric core-collapse supernova simulations performed with the Prometheus-Vertex code for 18 pre-supernova models in the range of 11-28 solar masses, including progenitors recently investigated by other groups. All models develop explosions, but depending on the progenitor structure, they can be divided into two classes. With a steep density decline at the Si/Si-O interface, the arrival of this interface at the shock front leads to a sudden drop of the mass-accretion rate, triggering a rapid approach to explosion. With a more gradually decreasing accretion rate, it takes longer for the neutrino heating to overcome the accretion ram pressure and explosions set in later. Early explosions are facilitated by high mass-accretion rates after bounce and correspondingly high neutrino luminosities combined with a pronounced drop of the accretion rate and ram pressure at the Si/Si-O interface. Because of rapidly shrinking neutron star radii and receding shock fronts after the passage through their maxima, our models exhibit short advection time scales, which favor the efficient growth of the standing accretion-shock instability. The latter plays a supportive role at least for the initiation of the re-expansion of the stalled shock before runaway. Taking into account the effects of turbulent pressure in the gain layer, we derive a generalized condition for the critical neutrino luminosity that captures the explosion behavior of all models very well. We validate the robustness of our findings by testing the influence of stochasticity, numerical resolution, and approximations in some aspects of the microphysics.

preprint2016arXiv

The Last Minutes of Oxygen Shell Burning in a Massive Star

We present the first 3D simulation of the last minutes of oxygen shell burning in an 18 solar mass supernova progenitor up to the onset of core collapse. A moving inner boundary is used to accurately model the contraction of the silicon and iron core according to a 1D stellar evolution model with a self-consistent treatment of core deleptonization and nuclear quasi-equilibrium. The simulation covers the full solid angle to allow the emergence of large-scale convective modes. Due to core contraction and the concomitant acceleration of nuclear burning, the convective Mach number increases to ~0.1 at collapse, and an l=2 mode emerges shortly before the end of the simulation. Aside from a growth of the oxygen shell from 0.51 to 0.56 solar masses due to entrainment from the carbon shell, the convective flow is reasonably well described by mixing length theory, and the dominant scales are compatible with estimates from linear stability analysis. We deduce that artificial changes in the physics, such as accelerated core contraction, can have precarious consequences for the state of convection at collapse. We argue that scaling laws for the convective velocities and eddy sizes furnish good estimates for the state of shell convection at collapse and develop a simple analytic theory for the impact of convective seed perturbations on shock revival in the ensuing supernova. We predict a reduction of the critical luminosity for explosion by 12--24% due to seed asphericities for our 3D progenitor model relative to the case without large seed perturbations.

preprint2015arXiv

The Dynamics of Neutrino-Driven Supernova Explosions after Shock Revival in 2D and 3D

We study the growth of the explosion energy after shock revival in neutrino-driven explosions in two and three dimensions (2D/3D) using multi-group neutrino hydrodynamics simulations of an $11.2 M_\odot$ star. The 3D model shows a faster and steadier growth of the explosion energy and already shows signs of subsiding accretion after one second. By contrast, the growth of the explosion energy in 2D is unsteady, and accretion lasts for several seconds as confirmed by additional long-time simulations of stars of similar masses. Appreciable explosion energies can still be reached, albeit at the expense of rather high neutron star masses. In 2D, the binding energy at the gain radius is larger because the strong excitation of downward-propagating g-modes removes energy from the freshly accreted material in the downflows. Consequently, the mass outflow rate is considerably lower in 2D than in 3D. This is only partially compensated by additional heating by outward-propagating acoustic waves in 2D. Moreover, the mass outflow rate in 2D is reduced because much of the neutrino energy deposition occurs in downflows or bubbles confined by secondary shocks without driving outflows. Episodic constriction of outflows and vertical mixing of colder shocked material and hot, neutrino-heated ejecta due to Rayleigh-Taylor instability further hamper the growth of the explosion energy in 2D. Further simulations will be necessary to determine whether these effects are generic over a wider range of supernova progenitors.

preprint2014arXiv

Experimental and computational studies of water drops falling through model oil with surfactant and subjected to an electric field

The behaviour of a single sub-millimetre-size water drop falling through a viscous oil while subjected to an electric field is of fundamental importance to industrial applications such as crude oil electrocoalescers. Detailed studies, both experimental and computational, have been performed previously, but an often challenging issue has been the characterization of the fluids. As numerous authors have noted, it is very difficult to have a perfectly clean water-oil system even for very pure model oils, and the presence of trace chemicals may significantly alter the interface behaviour. In this work, we consider a well- characterized water-oil system where controlled amounts of a surface active agent (Span 80) have been added to the oil. This addition dominates any trace contaminants in the oil, such that the interface behaviour can also be well-characterized. We present the results of experiments and corresponding two-phase- flow simulations of a falling water drop covered in surfactant and subjected to a monopolar square voltage pulse. The results are compared and good agreement is found for surfactant concentrations below the critical micelle concentration.

preprint2014arXiv

Extending a serial 3D two-phase CFD code to parallel execution over MPI by using the PETSc library for domain decomposition

To leverage the last two decades' transition in High-Performance Computing (HPC) towards clusters of compute nodes bound together with fast interconnects, a modern scalable CFD code must be able to efficiently distribute work amongst several nodes using the Message Passing Interface (MPI). MPI can enable very large simulations running on very large clusters, but it is necessary that the bulk of the CFD code be written with MPI in mind, an obstacle to parallelizing an existing serial code. In this work we present the results of extending an existing two-phase 3D Navier-Stokes solver, which was completely serial, to a parallel execution model using MPI. The 3D Navier-Stokes equations for two immiscible incompressible fluids are solved by the continuum surface force method, while the location of the interface is determined by the level-set method. We employ the Portable Extensible Toolkit for Scientific Computing (PETSc) for domain decomposition (DD) in a framework where only a fraction of the code needs to be altered. We study the strong and weak scaling of the resulting code. Cases are studied that are relevant to the fundamental understanding of oil/water separation in electrocoalescers.