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

Rodrigo Fernández

Rodrigo Fernández appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEastro-ph.SRnucl-thgr-qcastro-phastro-ph.GAastro-ph.IMhep-ph

Catalog footprint

What is connected

24works
8topics
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

24 published item(s)

preprint2025arXiv

Collapsar Disk Outflows III: Detectable Neutrino and Gravitational Wave Signatures

We investigate the neutrino and gravitational wave (GW) signals from accretion disks formed during the failed collapse of a rotating massive star (a collapsar). Following black hole formation, a neutrino-cooled, shocked accretion disk forms, which displays non-spherical oscillations for a period of seconds before becoming advective and exploding the star. We compute the neutrino and GW signals (matter quadrupole, frequencies $\lesssim 100$ Hz) from collapsar disks using global axisymmetric, viscous hydrodynamic simulations. The neutrino signal with typical energies of O$(10)$ MeV is maximal during the neutrino-cooled (NDAF) phase that follows shock formation. This phase lasts for a few seconds and is easily detectable within O$(10-100)$ kpc by the IceCube Neutrino Telescope. Additional neutrino signatures from a precursor equatorial shock and by stochastic accretion plumes during the advective phase are detectable within the galaxy. The GW signal during the NDAF phase is detectable in the galaxy by current and next-generation ground-based observatories. The explosion (memory) GW signal is similar to that of standard core-collapse supernovae and can be probed with a deci-Hertz space-based detector. Shock oscillations during the NDAF phase impart time variations with frequency O$(10-100)$ Hz to the neutrino and GW signals, encoding information about the shock dynamics and inner disk. These time variations can be detectable in neutrinos by IceCube within O$(1-10)$ kpc depending on progenitor model, flavor transformation scenario, and detailed properties of the angular momentum transport mechanism.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Long-term 3D-MHD Simulations of Black Hole Accretion Disks formed in Neutron Star Mergers

We examine the long-term evolution of accretion tori around black hole (BH) remnants of compact object mergers involving at least one neutron star, to better understand their contribution to kilonovae and the synthesis of r-process elements. To this end, we modify the unsplit magnetohydrodynamic (MHD) solver in FLASH4.5 to work in non-uniform three-dimensional spherical coordinates, enabling more efficient coverage of a large dynamic range in length scales while exploiting symmetries in the system. This modified code is used to perform BH accretion disk simulations that vary the initial magnetic field geometry and disk compactness, utilizing a physical equation of state, a neutrino leakage scheme for emission and absorption, and modeling the BH's gravity with a pseudo-Newtonian potential. Simulations run for long enough to achieve a radiatively-inefficient state in the disk. We find robust mass ejection with both poloidal and toroidal initial field geometries, and suppressed outflow at high disk compactness. With the included physics, we obtain bimodal velocity distributions that trace back to mass ejection by magnetic stresses at early times, and to thermal processes in the radiatively-inefficient state at late times. The electron fraction distribution of the disk outflow is broad in all models, and the ejecta geometry follows a characteristic hourglass shape. We test the effect of removing neutrino absorption or nuclear recombination with axisymmetric models, finding $\sim 50\%$ less mass ejection and more neutron-rich composition without neutrino absorption, and a subdominant contribution from nuclear recombination. Tests of the MHD and neutrino leakage implementations are included.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Probing magnetar emission mechanisms with spectropolarimetry

Over the next year, a new era of observations of compact objects in X-ray polarization will commence. Among the key targets for the upcoming Imaging X-ray Polarimetry Explorer mission, will be the magnetars 4U 0142+61 and 1RXS J170849.0-400910. Here we present the first detailed predictions of the expected polarization from these sources that incorporate realistic models of emission physics at the surface (gaseous or condensed), the temperature distribution on the surface, general relativity, quantum electrodynamics and scattering in the magnetosphere, and also account for the broadband spectral energy distribution of these sources from below 1 keV to nearly 100 keV. We find that either atmospheres or condensed surfaces can account for the emission at a few keV; in both cases either a small hot polar cap or scattering is required to account for the emission at 5-10 keV, and above 10 keV scattering by a hard population of electrons can account for the rising power in the hard X-rays observed in many magnetars in quiescence. Although these different scenarios result in very similar spectral energy distributions, they generate dramatically different polarization signatures from 2-10 keV, which is the range of sensitivity of the Imaging X-ray Polarimetry Explorer. Observations of these sources in X-ray polarization will therefore probe the emission from magnetars in an essentially new way.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Deep CFHT Optical Search for a Counterpart to the Possible Neutron Star -- Black Hole Merger GW190814

We present a wide-field optical imaging search for electromagnetic counterparts to the likely neutron star - black hole (NS-BH) merger GW190814/S190814bv. This compact binary merger was detected through gravitational waves by the LIGO/Virgo interferometers, with masses suggestive of a NS-BH merger. We imaged the LIGO/Virgo localization region using the MegaCam instrument on the Canada-France-Hawaii Telescope. We describe our hybrid observing strategy of both tiling and galaxy-targeted observations, as well as our image differencing and transient detection pipeline. Our observing campaign produced some of the deepest multi-band images of the region between 1.7 and 8.7 days post-merger, reaching a 5sigma depth of g > 22.8 (AB mag) at 1.7 days and i > 23.1 and i > 23.9 at 3.7 and 8.7 days, respectively. These observations cover a mean total integrated probability of 67.0% of the localization region. We find no compelling candidate transient counterparts to this merger in our images, which suggests that either the lighter object was tidally disrupted inside of the BH's innermost stable circular orbit, the transient lies outside of the observed sky footprint, or the lighter object is a low-mass BH. We use 5sigma source detection upper limits from our images in the NS-BH interpretation of this merger to constrain the mass of the kilonova ejecta to be Mej < 0.015Msun for a 'blue' (kappa = 0.5 cm^2 g^-1) kilonova, and Mej < 0.04Msun for a 'red' (kappa = 5-10 cm^2 g^-1) kilonova. Our observations emphasize the key role of large-aperture telescopes and wide-field imagers such as CFHT MegaCam in enabling deep searches for electromagnetic counterparts to gravitational wave events.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The landscape of disk outflows from black hole - neutron star mergers

We investigate mass ejection from accretion disks formed in mergers of black holes (BHs) and neutron stars (NSs). The third observing run of the LIGO/Virgo interferometers provided BH-NS candidate events that yielded no electromagnetic (EM) counterparts. The broad range of disk configurations expected from BH-NS mergers motivates a thorough exploration of parameter space to improve EM signal predictions. Here we conduct 27 high-resolution, axisymmetric, long-term hydrodynamic simulations of the viscous evolution of BH accretion disks that include neutrino emission/absorption effects and post-processing with a nuclear reaction network. In the absence of magnetic fields, these simulations provide a lower-limit to the fraction of the initial disk mass ejected. We find a nearly linear inverse dependence of this fraction on disk compactness (BH mass over initial disk radius). The dependence is related to the fraction of the disk mass accreted before the outflow is launched, which depends on the disk position relative to the innermost stable circular orbit. We also characterize a trend of decreasing ejected fraction and decreasing lanthanide/actinide content with increasing disk mass at fixed BH mass. This trend results from a longer time to reach weak freezout and an increasingly dominant role of neutrino absorption at higher disk masses. We estimate the radioactive luminosity from the disk outflow alone available to power kilonovae over the range of configurations studied, finding a spread of two orders of magnitude. For most of the BH-NS parameter space, the disk outflow contribution is well below the kilonova mass upper limits for GW190814.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

LRP2020: The cosmic origin and evolution of the elements

The origin of many elements of the periodic table remains an unsolved problem. While many nucleosynthetic channels are broadly understood, significant uncertainties remain regarding certain groups of elements such as the intermediate and rapid neutron-capture processes, the p-process, or the origin of odd-Z elements in the most metal-poor stars. Canada has a long tradition of leadership in nuclear astrophysics, dating back to the work of Alastair Cameron in the 1950s. Recent faculty hires have further boosted activity in the field, including transient observation and theory, survey science on galactic nucleosynthesis, and nuclear experiments. This white paper contains a brief overview of recent activity in the community, highlighting strengths in each sub-field, and provides recommendations to improve interdisciplinary collaboration. Sustaining Canadian leadership in the next decade will require, on the observational side, access to transient and non-transient surveys like LSST, SKA, or MSE, support for target-of-opportunity observing in current and future Canadian telescopes, and participation in next-generation X-ray telescopes such as ATHENA. State-of-the-art theoretical predictions will require an ambitious succession plan for the Niagara supercomputer to support large parallel jobs. We propose a funding instrument for postdoctoral training that reflects the interdisciplinary nature of nuclear astrophysics research, and the creation of a national collaborative funding program that allows for joint projects and workshop organization.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Production of the entire range of r-process nuclides by black hole accretion disc outflows from neutron star mergers

We consider $r$-process nucleosynthesis in outflows from black hole accretion discs formed in double neutron star and neutron star -- black hole mergers. These outflows, powered by angular momentum transport processes and nuclear recombination, represent an important -- and in some cases dominant -- contribution to the total mass ejected by the merger. Here we calculate the nucleosynthesis yields from disc outflows using thermodynamic trajectories from hydrodynamic simulations, coupled to a nuclear reaction network. We find that outflows produce a robust abundance pattern around the second $r$-process peak (mass number $A \sim 130$), independent of model parameters, with significant production of $A < 130$ nuclei. This implies that dynamical ejecta with high electron fraction may not be required to explain the observed abundances of $r$-process elements in metal poor stars. Disc outflows reach the third peak ($ A \sim 195$) in most of our simulations, although the amounts produced depend sensitively on the disc viscosity, initial mass or entropy of the torus, and nuclear physics inputs. Some of our models produce an abundance spike at $A = 132$ that is absent in the Solar system $r$-process distribution. The spike arises from convection in the disc and depends on the treatment of nuclear heating in the simulations. We conclude that disc outflows provide an important -- and perhaps dominant -- contribution to the $r$-process yields of compact binary mergers, and hence must be included when assessing the contribution of these systems to the inventory of $r$-process elements in the Galaxy.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Electromagnetic Signatures of Neutron Star Mergers in the Advanced LIGO Era

The mergers of binaries containing neutron stars and stellar-mass black holes are the most promising sources for direct detection in gravitational waves by the interferometers Advanced LIGO and Virgo over the next few years. The concurrent detection of electromagnetic emission from these events would greatly enhance the scientific return of these discoveries. Here we review the state of the art in modeling the electromagnetic signal of neutron star binary mergers across different phases of the merger and multiple wavelengths. We focus on those observables which provide the most sensitive diagnostics of the merger physics and the contribution to the synthesis of rapid neutron capture ($r$-process) elements in the Galaxy. We also outline expected future developments on the observational and theoretical sides of this rapidly evolving field.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The interplay of disk wind and dynamical ejecta in the aftermath of neutron star - black hole mergers

We explore the evolution of the different ejecta components generated during the merger of a neutron star (NS) and a black hole (BH). Our focus is the interplay between material ejected dynamically during the merger, and the wind launched on a viscous timescale by the remnant accretion disk. These components are expected to contribute to an electromagnetic transient and to produce r-process elements, each with a different signature when considered separately. Here we introduce a two-step approach to investigate their combined evolution, using two- and three-dimensional hydrodynamic simulations. Starting from the output of a merger simulation, we identify each component in the initial condition based on its phase space distribution, and evolve the accretion disk in axisymmetry. The wind blown from this disk is injected into a three-dimensional computational domain where the dynamical ejecta is evolved. We find that the wind can suppress fallback accretion on timescales longer than ~100 ms. Due to self-similar viscous evolution, the disk accretion at late times nevertheless approaches a power-law time dependence $\propto t^{-2.2}$. This can power some late-time GRB engine activity, although the available energy is significantly less than in traditional fallback models. Inclusion of radioactive heating due to the r-process does not significantly affect the fallback accretion rate or the disk wind. We do not find any significant modification to the wind properties at large radius due to interaction with the dynamical ejecta. This is a consequence of the different expansion velocities of the two components.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Three-Dimensional Simulations of SASI- and Convection-Dominated Core-Collapse Supernovae

We investigate the effect of dimensionality on the transition to explosion in neutrino-driven core-collapse supernovae. Using parameterized hydrodynamic simulations of the stalled supernova shock in one-, two- (2D), and three spatial dimensions (3D), we systematically probe the extent to which hydrodynamic instabilities alone can tip the balance in favor of explosion. In particular, we focus on systems that are well into the regimes where the Standing Accretion Shock Instability (SASI) or neutrino-driven convection dominate the dynamics, and characterize the difference between them. We find that SASI-dominated models can explode with up to ~20% lower neutrino luminosity in 3D than in 2D, with the magnitude of this difference decreasing with increasing resolution. This improvement in explosion conditions is related to the ability of spiral modes to generate more non-radial kinetic energy than a single sloshing mode, increasing the size of the average shock radius, and hence generating better conditions for the formation of large-scale, high-entropy bubbles. In contrast, convection-dominated explosions show a smaller difference in their critical heating rate between 2D and 3D (<8%), in agreement with previous studies. The ability of our numerical implementation to maintain arbitrary symmetries is quantified with a set of SASI-based tests. We discuss implications for the diversity of explosion paths in a realistic supernova environment.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Angular momentum redistribution by SASI spiral modes and consequences for neutron star spins

In the collapsing core of massive stars, the standing accretion shock instability (SASI) can drive spiral modes that efficiently redistribute angular momentum. This process can impart a spin to the forming neutron star even when the progenitor star is non-rotating. Here we develop the first analytical description of the angular momentum redistribution driven by a spiral mode of the SASI. Our analysis, valid in the limit of small mode amplitude, shows that the angular momentum separation is driven by the Reynolds stress generated by the spiral mode. The resulting solutions compare favorably with previous three-dimensional hydrodynamic simulations of the SASI in the linear and weakly non-linear phases. Reasonable agreement is also found when extrapolating the solutions into the fully non-linear phase. A Reynolds-decomposition of the flow is performed in the saturated state of these simulations, showing that outward angular momentum transport by the Reynolds stress and the fluctuating component of the mass flux balance inward transport by advection. We derive an approximate analytic expression for the maximum angular momentum deposited in the neutron star as a function of the mass accretion rate, shock radius, shock compression ratio, and amplitude of the spiral mode at the time of explosion. Implications for the birth spin periods of neutron stars are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Characterizing SASI- and Convection-Dominated Core-Collapse Supernova Explosions in Two Dimensions

The success of the neutrino mechanism of core-collapse supernovae relies on the supporting action of two hydrodynamic instabilities: neutrino-driven convection and the Standing Accretion Shock Instability (SASI). Depending on the structure of the stellar progenitor, each of these instabilities can dominate the evolution of the gain region prior to the onset of explosion, with implications for the ensuing asymmetries. Here we examine the flow dynamics in the neighborhood of explosion by means of parametric two-dimensional, time-dependent hydrodynamic simulations for which the linear stability properties are well understood. We find that systems for which the convection parameter is sub-critical (SASI-dominated) develop explosions once large-scale, high-entropy bubbles are able to survive for several SASI oscillation cycles. These long-lived structures are seeded by the SASI during shock expansions. Finite-amplitude initial perturbations do not alter this outcome qualitatively, though they can lead to significant differences in explosion times. Supercritical systems (convection-dominated) also explode by developing large-scale bubbles, though the formation of these structures is due to buoyant activity. Non-exploding systems achieve a quasi-steady state in which the time-averaged flow adjusts itself to be convectively sub-critical. We characterize the turbulent flow using a spherical Fourier-Bessel decomposition, identifying the relevant scalings and connecting temporal and spatial components. Finally, we verify the applicability of these principles on the general relativistic, radiation-hydrodynamic simulations of Mueller, Janka, & Heger (2012), and discuss implications for the three-dimensional case.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Outflows from accretion disks formed in neutron star mergers: effect of black hole spin

The accretion disk that forms after a neutron star merger is a source of neutron-rich ejecta. The ejected material contributes to a radioactively-powered electromagnetic transient, with properties that depend sensitively on the composition of the outflow. Here we investigate how the spin of the black hole remnant influences mass ejection on the thermal and viscous timescales. We carry out two-dimensional, time-dependent hydrodynamic simulations of merger remnant accretion disks including viscous angular momentum transport and approximate neutrino self-irradiation. The gravity of the spinning black hole is included via a pseudo-Newtonian potential. We find that a disk around a spinning black hole ejects more mass, up to a factor of several, relative to the non-spinning case. The enhanced mass loss is due to energy release by accretion occurring deeper in the gravitational potential, raising the disk temperature and hence the rate of viscous heating in regions where neutrino cooling is ineffective. The mean electron fraction of the outflow increases moderately with BH spin due to a highly-irradiated (though not neutrino-driven) wind component. While the bulk of the ejecta is still very neutron-rich, thus generating heavy r-process elements, the leading edge of the wind contains a small amount of Lanthanide-free material. This component can give rise to a ~1 day blue optical `bump' in a kilonova light curve, even in the case of prompt BH formation, which may facilitate its detection.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Red or blue? A potential kilonova imprint of the delay until black hole formation following a neutron star merger

Mergers of binary neutron stars (NSs) usually result in the formation of a hypermassive neutron star (HMNS). Whether- and when this remnant collapses to a black hole (BH) depends primarily on the equation of state and on angular momentum transport processes, both of which are uncertain. Here we show that the lifetime of the merger remnant may be directly imprinted in the radioactively powered kilonova emission following the merger. We employ axisymmetric, time-dependent hydrodynamic simulations of remnant accretion disks orbiting a HMNS of variable lifetime, and characterize the effect of this delay to BH formation on the disk wind ejecta. When BH formation is relatively prompt (~ 100 ms), outflows from the disk are sufficiently neutron rich to form heavy r-process elements, resulting in ~ week-long emission with a spectral peak in the near-infrared (NIR), similar to that produced by the dynamical ejecta. In contrast, delayed BH formation allows neutrinos from the HMNS to raise the electron fraction in the polar direction to values such that potentially Lanthanide-free outflows are generated. The lower opacity would produce a brighter, bluer, and shorter-lived ~ day-long emission (a `blue bump') prior to the late NIR peak from the dynamical ejecta and equatorial wind. This new diagnostic of BH formation should be useful for events with a signal to noise lower than that required for direct detection of gravitational waveform signatures.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Delayed outflows from black hole accretion tori following neutron star binary coalescence

Expulsion of neutron-rich matter following the merger of neutron star (NS) binaries is crucial to the radioactively-powered electromagnetic counterparts of these events and to their relevance as sources of r-process nucleosynthesis. Here we explore the long-term (viscous) evolution of remnant black hole accretion disks formed in such mergers by means of two-dimensional, time-dependent hydrodynamical simulations. The evolution of the electron fraction due to charged-current weak interactions is included, and neutrino self-irradiation is modeled as a lightbulb that accounts for the disk geometry and moderate optical depth effects. Over several viscous times (~1s), a fraction ~10% of the initial disk mass is ejected as a moderately neutron-rich wind (Y_e ~ 0.2) powered by viscous heating and nuclear recombination, with neutrino self-irradiation playing a sub-dominant role. Although the properties of the outflow vary in time and direction, their mean values in the heavy-element production region are relatively robust to variations in the initial conditions of the disk and the magnitude of its viscosity. The outflow is sufficiently neutron-rich that most of the ejecta forms heavy r-process elements with mass number A >130, thus representing a new astrophysical source of r-process nucleosynthesis, distinct from that produced in the dynamical ejecta. Due to its moderately high entropy, disk outflows contain a small residual fraction ~1% of helium, which could produce a unique spectroscopic signature.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Head-on collisions of white dwarfs in triple systems could explain type Ia supernova

Type Ia supernovae (SNe Ia), thermonuclear explosions of carbon-oxygen white dwarfs (CO-WDs), are currently the best cosmological "standard candles", but the triggering mechanism of the explosion is unknown. It was recently shown that the rate of head-on collisions of typical field CO-WDs in triple systems may be comparable to the SNe Ia rate. Here we provide evidence supporting a scenario in which the majority of SNe Ia are the result of such head-on collisions of CO-WDs. In this case, the nuclear detonation is due to a well understood shock ignition, devoid of commonly introduced free parameters such as the deflagration velocity or transition to detonation criteria. By using two-dimensional hydrodynamical simulations with a fully resolved ignition process, we show that zero-impact-parameter collisions of typical CO-WDs with masses $0.5-1\,M_{\odot}$ result in explosions that synthesize $^{56}$Ni masses in the range of $\sim0.1-1\,M_{\odot}$, spanning the wide distribution of yields observed for the majority of SNe Ia. All collision models yield the same late-time ($>60$ days since explosion) bolometric light curve when normalized by $^{56}$Ni masses (to better than $30\%$), in agreement with observations. The calculated widths of the $^{56}$Ni-mass-weighted-line-of-sight velocity distributions are correlated with the calculated $^{56}$Ni yield, agreeing with the observed correlation. The strong correlation, shown here for the first time, between $^{56}$Ni yield and total mass of the colliding CO-WDs (insensitive to their mass ratio), is suggestive as the source for the continuous distribution of observed SN Ia features, possibly including the Philips relation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Nonlinear Evolution of the Radiation-Driven Magneto-Acoustic Instability (RMI)

We examine the nonlinear development of unstable magnetosonic waves driven by a background radiative flux -- the Radiation-Driven Magneto-Acoustic Instability (RMI, a.k.a. the "photon bubble" instability). The RMI may serve as a persistent source of density, radiative flux, and magnetic field fluctuations in stably-stratified, optically-thick media. The conditions for instability are present in a variety of astrophysical environments, and do not require the radiation pressure to dominate or the magnetic field to be strong. Here we numerically study the saturation properties of the RMI, covering three orders of magnitude in the relative strength of radiation, magnetic field, and gas energies. Two-dimensional, time-dependent radiation-MHD simulations of local, stably-stratified domains are conducted with Zeus-MP in the optically-thick, highly-conducting limit. Our results confirm the theoretical expectations of Blaes and Socrates (2003) in that the RMI operates even in gas pressure-dominated environments that are weakly magnetized. The saturation amplitude is a monotonically increasing function of the ratio of radiation to gas pressure. Keeping this ratio constant, we find that the saturation amplitude peaks when the magnetic pressure is comparable to the radiation pressure. We discuss the implications of our results for the dynamics of magnetized stellar envelopes, where the RMI should act as a source of sub-photospheric perturbations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Hydrodynamics of Core-Collapse Supernovae at the Transition to Explosion. I. Spherical Symmetry

We study the transition to runaway expansion of an initially stalled core-collapse supernova shock. The neutrino luminosity, mass accretion rate, and neutrinospheric radius are all treated as free parameters. In spherical symmetry, this transition is mediated by a global non-adiabatic instability that develops on the advection time and reaches non-linear amplitude. Here we perform high-resolution, time-dependent hydrodynamic simulations of stalled supernova shocks with realistic microphysics to analyze this transition. We find that radial instability is a sufficient condition for runaway expansion if the neutrinospheric parameters do not vary with time and if heating by the accretion luminosity is neglected. For a given unstable mode, transition to runaway occurs when fluid in the gain region reaches positive specific energy. We find approximate instability criteria that accurately describe the behavior of the system over a wide region of parameter space. The threshold neutrino luminosities are in general different than the limiting value for a steady-state solution. We hypothesize that multidimensional explosions arise from the excitation of unstable large-scale modes of the turbulent background flow, at threshold luminosities that are lower than in the laminar case.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Nuclear Dominated Accretion Flows in Two Dimensions. I. Torus Evolution with Parametric Microphysics

We explore the evolution of radiatively inefficient accretion disks in which nuclear reactions are dynamically important (`Nuclear Dominated Accretion Flows', or NuDAFs). Examples of such disks are those generated by the merger of a white dwarf with a neutron star or black hole, or by the collapse of a rotating star. Here we present two-dimensional hydrodynamic simulations that systematically explore the effect of adding a single nuclear reaction to a viscous torus. The equation of state, anomalous shear stress, and nuclear reactions are given parametric forms. Our results point to the existence of two qualitatively different regimes of NuDAF evolution: (1) steady accretion with quiescent burning; or (2) detonation of the disk. These outcomes are controlled primarily by the ratio of the nuclear energy released to the enthalpy at the burning radius. Disks detonate if this ratio exceeds a critical value (~1), and if burning occurs in regions where neutrino cooling is unimportant. Thermonuclear runaways are seeded by the turbulent mixing of hot ash with cold fuel at the burning front. Subcritical disks do not explode, but instead power a persistent collimated outflow of unbound material composed primarily of ash, with a mass-loss rate that increases with the nuclear binding energy released. We discuss the implications of our results for supernova-like counterparts from astrophysical events in the NuDAF regime. In particular, detonations following a white dwarf - neutron star merger could account for some subluminous Type Ia supernovae, such as the class defined by SN 2002cx.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The X-Ray Polarization Signature of Quiescent Magnetars: Effect of Magnetospheric Scattering and Vacuum Polarization

In the magnetar model, the quiescent non-thermal soft X-ray emission from Anomalous X-ray Pulsars and Soft-Gamma Repeaters is thought to arise from resonant comptonization of thermal photons by charges moving in a twisted magnetosphere. Robust inference of physical quantities from observations is difficult, because the process depends strongly on geometry and current understanding of the magnetosphere is not very deep. The polarization of soft X-ray photons is an independent source of information, and its magnetospheric imprint remains only partially explored. In this paper we calculate how resonant cyclotron scattering would modify the observed polarization signal relative to the surface emission, using a multidimensional Monte Carlo radiative transfer code that accounts for the gradual coupling of polarization eigenmodes as photons leave the magnetosphere. We employ a globally-twisted, self-similar, force-free magnetosphere with a power-law momentum distribution, assume a blackbody spectrum for the seed photons, account for general relativistic light deflection close to the star, and assume that vacuum polarization dominates the dielectric properties of the magnetosphere. The latter is a good approximation if the pair multiplicity is not much larger than unity. Phase-averaged polarimetry is able to provide a clear signature of the magnetospheric reprocessing of thermal photons and to constrain mechanisms generating the thermal emission. Phase-resolved polarimetry, in addition, can characterize the spatial extent and magnitude of the magnetospheric twist angle at ~100 stellar radii, and discern between uni- or bidirectional particle energy distributions, almost independently of every other parameter in the system. We discuss prospects for detectability with GEMS.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Spiral Modes of the Standing Accretion Shock Instability

A stalled spherical accretion shock, such as that arising in core-collapse supernovae, is unstable to non-spherical perturbations. In three dimensions, this Standing Accretion Shock Instability (SASI) can develop spiral modes that spin-up the protoneutron star. Here we study these non-axisymmetric modes by combining linear stability analysis and three-dimensional, time-dependent hydrodynamic simulations with Zeus-MP, focusing on characterizing their spatial structure and angular momentum content. We do not impose any rotation on the background accretion flow, and use simplified microphysics with no neutrino heating or nuclear dissociation. Spiral modes are examined in isolation by choosing flow parameters such that only the fundamental mode is unstable for a given polar index, leading to good agreement with linear theory. We find that any superposition of sloshing modes with non-zero relative phases survives in the nonlinear regime and leads to angular momentum redistribution. It follows that the range of perturbations required to obtain spin-up is broader than that needed to obtain the limiting case of a phase shift of pi/2. The bulk of the angular momentum redistribution occurs during a phase of exponential growth, and arises from internal torques that are second order in the perturbation amplitude. This redistribution gives rise to at least two counter rotating regions, with the maximum angular momentum of a given sign approaching a significant fraction of the mass accretion rate times the shock radius squared (Mdot * r_shock^2 ~ 1E+47 g/cm^2/s, spin period ~60 ms). Nonlinear mode coupling at saturation causes the angular momentum to fluctuate in all directions with much smaller amplitudes.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Dynamics of a Spherical Accretion Shock with Neutrino Heating and Alpha-Particle Recombination

We investigate the effects of neutrino heating and alpha-particle recombination on the hydrodynamics of core-collapse supernovae. Our focus is on the non-linear dynamics of the shock wave that forms in the collapse, and the assembly of positive energy material below it. To this end, we perform time-dependent hydrodynamic simulations with FLASH2.5 in spherical and axial symmetry. These generalize our previous calculations by allowing for bulk neutrino heating and for nuclear statistical equilibrium between n, p and alpha. The heating rate is freely tunable, as is the starting radius of the shock relative to the recombination radius of alpha-particles. An explosion in spherical symmetry involves the excitation of an overstable mode, which may be viewed as the L=0 version of the `Standing Accretion Shock Instability'. In 2D simulations, non-spherical deformations of the shock are driven by plumes of material with positive Bernoulli parameter, which are concentrated well outside the zone of strong neutrino heating. The non-spherical modes of the shock reach a large amplitude only when the heating rate is also high enough to excite convection below the shock. The critical heating rate that causes an explosion depends sensitively on the initial position of the shock relative to the recombination radius. Weaker heating is required to drive an explosion in 2D than in 1D, but the difference also depends on the size of the shock. Forcing the infalling heavy nuclei to break up into n and p below the shock only causes a slight increase in the critical heating rate, except when the shock starts out at a large radius. This shows that heating by neutrinos (or some other mechanism) must play a significant role in pushing the shock far enough out that recombination heating takes over.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Stability of a Spherical Accretion Shock with Nuclear Dissociation

We examine the stability of a standing shock wave within a spherical accretion flow onto a gravitating star, in the context of core-collapse supernova explosions. Our focus is on the effect of nuclear dissociation below the shock on the linear growth, and non-linear saturation, of non-radial oscillations of the shocked fluid. We combine two-dimensional, time-dependent hydrodynamic simulations using FLASH2.5 with a solution to the linear eigenvalue problem, and demonstrate the consistency of the two approaches. Previous studies of this `Standing Accretion Shock Instability' (SASI) have focused either on zero-energy accretion flows without nuclear dissociation, or made use of a detailed finite-temperature nuclear equation of state and included strong neutrino heating. Our main goal in this and subsequent papers is to introduce equations of state of increasing complexity, in order to isolate the various competing effects. In this work we employ an ideal gas equation of state with a constant rate of nuclear dissociation below the shock, and do not include neutrino heating. We find that a negative Bernoulli parameter below the shock significantly lowers the real frequency, growth rate, and saturation amplitude of the SASI. A decrease in the adiabatic index has similar effects. The non-linear development of the instability is characterized by an expansion of the shock driven by turbulent kinetic energy at nearly constant internal energy. Our results also provide further insight into the instability mechanism: the rate of growth of a particular mode is fastest when the radial advection time from the shock to the accretor overlaps with the period of a standing lateral sound wave. The fastest-growing mode can therefore be modified by nuclear dissociation.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Rotochemical Heating in Millisecond Pulsars. Formalism and Non-superfluid case

Rotochemical heating originates in a departure from beta equilibrium due to spin-down compression in a rotating neutron star. The main consequence is that the star eventually arrives at a quasi-equilibrium state, in which the thermal photon luminosity depends only on the current value of the spin-down power, which is directly measurable. Only in millisecond pulsars the spin-down power remains high long enough for this state to be reached with a substantial luminosity. We report an extensive study of the effect of this heating mechanism on the thermal evolution of millisecond pulsars, developing a general formalism in the slow-rotation approximation of general relativity that takes the spatial structure of the star fully into account, and using a sample of realistic equations of state to solve the non-superfluid case numerically. We show that nearly all observed millisecond pulsars are very likely to be in the quasi-equilibrium state. Our predicted quasi-equilibrium temperatures for PSR J0437-4715 are only 20% lower than inferred from observations. Accounting for superfluidity should increase the predicted value.

0 citations0 reviewsNo code linkedOpen access