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

Evan O'Connor

Evan O'Connor appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEgr-qcastro-ph.SRnucl-thhep-exastro-ph.IMhep-phnucl-ex

Catalog footprint

What is connected

23works
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

23 published item(s)

preprint2023arXiv

Impact of Rotation on the Multimessenger Signatures of a Hadron-quark Phase Transition in Core-collapse Supernovae

We study the impact of rotation on the multimessenger signals of core-collapse supernovae (CCSNe) with the occurrence of a first-order hadron-quark phase transition (HQPT). We simulate CCSNe with the \texttt{FLASH} code starting from a 20~$M_\odot$ progenitor with different rotation rates, and using the RDF equation of state from \textit{Bastian} 2021 that prescribes the HQPT. Rotation is found to delay the onset of the HQPT and the resulting dynamical collapse of the protocompact star (PCS) due to the centrifugal support. All models with the HQPT experience a second bounce shock which leads to a successful explosion. The oblate PCS as deformed by rotation gives rise to strong gravitational-wave (GW) emission around the second bounce with a peak amplitude larger by a factor of $\sim10$ than that around the first bounce. The breakout of the second bounce shock at the neutrinosphere produces a $\barν_e$-rich neutrino burst with a luminosity of serveral 10$^{53}$~erg~s$^{-1}$. In rapidly rotating models the PCS pulsation following the second bounce generates oscillations in the neutrino signal after the burst. In the fastest rotating model with the HQPT, a clear correlation is found between the oscillations in the GW and neutrino signals immediately after the second bounce. In addition, the HQPT-induced collapse leads to a jump in the ratio of rotational kinetic energy to gravitational energy ($β$) of the PCS, for which persistent GW emission may arise due to secular nonaxisymmetric instabilities.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Measuring the distance and mass of galactic core-collapse supernovae using neutrinos

Neutrinos from a Galactic core-collapse supernova will be measured by neutrino detectors minutes to days before an optical signal reaches Earth. We present a novel calculation showing the ability of current and near-future neutrino detectors to make fast predictions of the progenitor distance and place constraints on the zero-age main sequence mass in order to inform the observing strategy for electromagnetic follow-up. We show that for typical Galactic supernovae, the distance can be constrained with an uncertainty of $\sim$5\% using IceCube or Hyper-K and, furthermore, the zero-age main sequence mass can be constrained for extremal values of compactness.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Neutrino Echos following Black Hole Formation in Core-Collapse Supernovae

During a failed core-collapse supernova, the protoneutron star eventually collapses under its own gravitational field and forms a black hole. This collapse happens quickly, on the dynamical time of the protoneutron star, $\lesssim$0.5 ms. During this collapse, barring any excessive rotation, the entire protoneutron star is accreted into the newly formed black hole. The main source of neutrinos is now removed and the signal abruptly shuts off over this formation timescale. However, while the source of neutrinos is turned off, the arrival times at an Earth-based detector will depend on the neutrino path. We show here that a modest amount of neutrinos, emitted just prior to the black hole forming, scatter on the infalling material into our line of sight and arrive after the formation of the black hole, up to 15 ms in our model. This neutrino echo, which we characterize with Monte Carlo simulations and analytic models, has a significantly higher average energy (upwards of $\sim$ 50 MeV) compared to the main neutrino signal, and for the canonical failed supernova explored here, is likely detectable in $\mathcal{O}$(10 kT) supernova neutrino detectors for Galactic failed supernovae. The presence of this signal is important to consider if using black hole formation as a time post for triangulation or the post black hole timing profile for neutrino mass measurements. On its own, it can also be used to characterize or constrain the structure and nature of the accretion flow.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Impact of neutrino pair-production rates in Core-Collapse Supernovae

In this paper, we present a careful study on the impact of neutrino pair-production on core-collapse supernovae via spherically-symmetric, general-relativistic simulations of two different massive star progenitors with energy-dependent neutrino transport. We explore the impact and consequences of both the underlying microphysics and the implementation in the radiation transport algorithms on the supernova evolution, neutrino signal properties, and the explosion dynamics. We consider the two dominant neutrino pair-production processes found in supernovae, electron-positron annihilation as well as nucleon-nucleon bremsstrahlung in combination with both a simplified and a complete treatment of the processes in the radiation transport algorithms. We find that the use of the simplified prescription quantitatively impacts the neutrino signal at the 10\% level and potentially the supernova dynamics, as we show for the case of a 9.6M$_\odot$ progenitor. We also show that the choice of nucleon-nucleon bremsstrahlung interaction can also have a quantitative impact on the neutrino signal. A self-consistent treatment with state-of-the-art microphysics is suggested for precision simulations of core collapse, however the simplified treatment explored here is both computationally less demanding and results in a qualitatively similar evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Black Hole-Neutron Star Mergers with a Hot Nuclear Equation of State: Outflow and Neutrino-Cooled Disk for a Low-Mass, High-Spin Case

Neutrino emission significantly affects the evolution of the accretion tori formed in black hole-neutron star mergers. It removes energy from the disk, alters its composition, and provides a potential power source for a gamma-ray burst. To study these effects, simulations in general relativity with a hot microphysical equation of state and neutrino feedback are needed. We present the first such simulation, using a neutrino leakage scheme for cooling to capture the most essential effects and considering a moderate mass (1.4 M_{\odot} neutron star, 5.6 M_{\odot} black hole), high spin (black hole J/M^2=0.9) system with the K_0=220 MeV Lattimer-Swesty equation of state. We find that about 0.08 M_{\odot} of nuclear matter is ejected from the system, while another 0.3 M_{\odot} forms a hot, compact accretion disk. The primary effects of the escaping neutrinos are (i) to make the disk much denser and more compact, (ii) to cause the average electron fraction Y_e of the disk to rise to about 0.2 and then gradually decrease again, and (iii) to gradually cool the disk. The disk is initially hot (T~6 MeV) and luminous in neutrinos (L_ν~10^{54} erg s^{-1}), but the neutrino luminosity decreases by an order of magnitude over 50 ms of post-merger evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Impact of an improved neutrino energy estimate on outflows in neutron star merger simulations

Binary neutron star mergers are promising sources of gravitational waves for ground-based detectors such as Advanced LIGO. Neutron-rich material ejected by these mergers may also be the main source of r-process elements in the Universe, while radioactive decays in the ejecta can power bright electromagnetic post-merger signals. Neutrino-matter interactions play a critical role in the evolution of the composition of the ejected material, which significantly impacts the outcome of nucleosynthesis and the properties of the associated electromagnetic signal. In this work, we present a simulation of a binary neutron star merger using an improved method for estimating the average neutrino energies in our energy-integrated neutrino transport scheme. These energy estimates are obtained by evolving the neutrino number density in addition to the neutrino energy and flux densities. We show that significant changes are observed in the composition of the polar ejecta when comparing our new results with earlier simulations in which the neutrino spectrum was assumed to be the same everywhere in optically thin regions. In particular, we find that material ejected in the polar regions is less neutron rich than previously estimated. Our new estimates of the composition of the polar ejecta make it more likely that the color and timescale of the electromagnetic signal depend on the orientation of the binary with respect to an observer's line-of-sight. These results also indicate that important observable properties of neutron star mergers are sensitive to the neutrino energy spectrum, and may need to be studied through simulations including a more accurate, energy-dependent neutrino transport scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Low mass binary neutron star mergers : gravitational waves and neutrino emission

Neutron star mergers are among the most promising sources of gravitational waves for advanced ground-based detectors. These mergers are also expected to power bright electromagnetic signals, in the form of short gamma-ray bursts, infrared/optical transients, and radio emission. Simulations of these mergers with fully general relativistic codes are critical to understand the merger and post-merger gravitational wave signals and their neutrinos and electromagnetic counterparts. In this paper, we employ the SpEC code to simulate the merger of low-mass neutron star binaries (two $1.2M_\odot$ neutron stars) for a set of three nuclear-theory based, finite temperature equations of state. We show that the frequency peaks of the post-merger gravitational wave signal are in good agreement with predictions obtained from simulations using a simpler treatment of gravity. We find, however, that only the fundamental mode of the remnant is excited for long periods of time: emission at the secondary peaks is damped on a millisecond timescale in the simulated binaries. For such low-mass systems, the remnant is a massive neutron star which, depending on the equation of state, is either permanently stable or long-lived. We observe strong excitations of l=2, m=2 modes, both in the massive neutron star and in the form of hot, shocked tidal arms in the surrounding accretion torus. We estimate the neutrino emission of the remnant using a neutrino leakage scheme and, in one case, compare these results with a gray two-moment neutrino transport scheme. We confirm the complex geometry of the neutrino emission, also observed in previous simulations with neutrino leakage, and show explicitly the presence of important differences in the neutrino luminosity, disk composition, and outflow properties between the neutrino leakage and transport schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Supernova Physics at DUNE

The DUNE/LBNF program aims to address key questions in neutrino physics and astroparticle physics. Realizing DUNE's potential to reconstruct low-energy particles in the 10-100 MeV energy range will bring significant benefits for all DUNE's science goals. In neutrino physics, low-energy sensitivity will improve neutrino energy reconstruction in the GeV range relevant for the kinematics of DUNE's long-baseline oscillation program. In astroparticle physics, low-energy capabilities will make DUNE's far detectors the world's best apparatus for studying the electron-neutrino flux from a supernova. This will open a new window to unrivaled studies of the dynamics and neutronization of a star's central core in real time, the potential discovery of the neutrino mass hierarchy, provide new sensitivity to physics beyond the Standard Model, and evidence of neutrino quantum-coherence effects. The same capabilities will also provide new sensitivity to `boosted dark matter' models that are not observable in traditional direct dark matter detectors.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Unequal mass binary neutron star mergers and multimessenger signals

We study the merger of binary neutron stars with different mass ratios adopting three different realistic, microphysical nuclear equations of state, as well as incorporating neutrino cooling effects. In particular, we concentrate on the influence of the equation of state on the gravitational wave signature and also on its role, in combination with neutrino cooling, in determining the properties of the resulting hypermassive neutron star, of the neutrinos produced, and of the ejected material. The ejecta we find are consistent with other recent studies that find that small mass ratios produce more ejecta than equal mass cases (up to some limit) and this ejecta is more neutron rich. This trend indicates the importance with future kilonovae observations of measuring the individual masses of an associated binary neutron star system, presumably from concurrent gravitational wave observations, in order to be able to extract information about the nuclear equation of state

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An Open-Source Neutrino Radiation Hydrodynamics Code for Core-Collapse Supernovae

We present an open-source update to the spherically-symmetric, general-relativistic hydrodynamics, core-collapse supernova (CCSN) code GR1D (O'Connor & Ott, 2010, CQG, 27, 114103). The source code is available at http://www.GR1Dcode.org. We extend its capabilities to include a general relativistic treatment of neutrino transport based on the moment formalisms of Shibata et al., 2011, PTP, 125, 1255 and Cardall et al., 2013, PRD, 87 103004. We pay special attention to implementing and testing numerical methods and approximations that lessen the computational demand of the transport scheme by removing the need to invert large matrices. This is especially important for the implementation and development of moment-like transport methods in two and three dimensions. A critical component of neutrino transport calculations are the neutrino-matter interaction coefficients that describe the production, absorption, scattering, and annihilation of neutrinos. In this article we also describe our open-source, neutrino interaction library NuLib (available at http://www.nulib.org). We believe that an open-source approach to describing these interactions is one of the major steps needed to progress towards robust models of CCSNe and robust predictions of the neutrino signal. We show, via comparisons to full Boltzmann neutrino transport simulations of CCSNe, that our neutrino transport code performs remarkably well. Furthermore, we show that the methods and approximations we employ to increase efficiency do not decrease the fidelity of our results. We also test the ability of our general relativistic transport code to model failed CCSN by evolving a 40 solar-mass progenitor to the onset of collapse to a black hole.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Monte Carlo Neutrino Transport Through Remnant Disks from Neutron Star Mergers

We present Sedonu, a new open source, steady-state, special relativistic Monte Carlo (MC) neutrino transport code, available at bitbucket.org/srichers/sedonu. The code calculates the energy- and angle-dependent neutrino distribution function on fluid backgrounds of any number of spatial dimensions, calculates the rates of change of fluid internal energy and electron fraction, and solves for the equilibrium fluid temperature and electron fraction. We apply this method to snapshots from two-dimensional simulations of accretion disks left behind by binary neutron star mergers, varying the input physics and comparing to the results obtained with a leakage scheme for the case of a central black hole and a central hypermassive neutron star. Neutrinos are guided away from the densest regions of the disk and escape preferentially around 45 degrees from the equatorial plane. Neutrino heating is strengthened by MC transport a few scale heights above the disk midplane near the innermost stable circular orbit, potentially leading to a stronger neutrino-driven wind. Neutrino cooling in the dense midplane of the disk is stronger when using MC transport, leading to a globally higher cooling rate by a factor of a few and a larger leptonization rate by an order of magnitude. We calculate neutrino pair annihilation rates and estimate that an energy of 2.8e46 erg is deposited within 45 degrees of the symmetry axis over 300 ms when a central BH is present. Similarly, 1.9e48 erg is deposited over 3 s when an HMNS sits at the center, but neither estimate is likely to be sufficient to drive a GRB jet.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Post-merger evolution of a neutron star-black hole binary with neutrino transport

We present a first simulation of the post-merger evolution of a black hole-neutron star binary in full general relativity using an energy-integrated general relativistic truncated moment formalism for neutrino transport. We describe our implementation of the moment formalism and important tests of our code, before studying the formation phase of a disk after a black hole-neutron star merger. We use as initial data an existing general relativistic simulation of the merger of a neutron star of 1.4 solar mass with a black hole of 7 solar mass and dimensionless spin a/M=0.8. Comparing with a simpler leakage scheme for the treatment of the neutrinos, we find noticeable differences in the neutron to proton ratio in and around the disk, and in the neutrino luminosity. We find that the electron neutrino luminosity is much lower in the transport simulations, and that the remnant is less neutron-rich. The spatial distribution of the neutrinos is significantly affected by relativistic effects. Over the short timescale evolved, we do not observe purely neutrino-driven outflows. However, a small amount of material (3e-4Msun) is ejected in the polar region during the circularization of the disk. Most of that material is ejected early in the formation of the disk, and is fairly neutron rich. Through r-process nucleosynthesis, that material should produce high-opacity lanthanides in the polar region, and could thus affect the lightcurve of radioactively powered electromagnetic transients. We also show that by the end of the simulation, while the bulk of the disk is neutron-rich, its outer layers have a higher electron fraction. As that material would be the first to be unbound by disk outflows on longer timescales, the changes in Ye experienced during the formation of the disk could have an impact on the nucleosynthesis outputs from neutrino-driven and viscously-driven outflows. [Abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Sensitivity of Core-Collapse Supernovae to Nuclear Electron Capture

A weak-rate library aimed at investigating the sensitivity of astrophysical environments to variations of electron-capture rates on medium-heavy nuclei has been developed. With this library, the sensitivity of the core-collapse and early post-bounce phases of core-collapse supernovae to nuclear electron-capture is examined by systematically and statistically varying electron-capture rates of individual nuclei. The rates are adjusted by factors consistent with uncertainties indicated by comparing theoretical rates to those deduced from charge-exchange and $β$-decay measurements. To ensure a model independent assessment, sensitivity studies across a comprehensive set of progenitors and equations of state are performed. In our systematic study, we find a +16/-4 % range in the mass of the inner-core at the time of shock formation and a $\pm$20% range of peak νe-luminosity during the deleptonization burst. These ranges are each five times as large as those seen from a separate progenitor study in which we evaluated the sensitivity of these parameters to 32 presupernova stellar models. It is also found that the simulations are more sensitive to a reduction in the electron-capture rates than an enhancement, and in particular to the reduction in the rates for neutron-rich nuclei near the N = 50 closed neutron-shell. As measurements for medium-heavy (A > 65) and neutron-rich nuclei are sparse, and because accurate theoretical models which account for nuclear structure considerations on the individual nucleus level are not readily available, rates for these nuclei may be overestimated. If more accurate estimates confirm this, results from this study indicate that significant changes to the core-collapse trajectory can be expected. For this reason, experimental and theoretical efforts should focus in this region of the nuclear chart.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Magnetized Neutron Stars With Realistic Equations of State and Neutrino Cooling

We incorporate realistic, tabulated equations of state into fully relativistic simulations of magnetized neutron stars along with a neutrino leakage scheme which accounts for cooling via neutrino emission. Both these improvements utilize open-source code (GR1D) and tables from http://stellarcollapse.org. Our implementation makes use of a novel method for the calculation of the optical depth which simplifies its use with distributed adaptive mesh refinement. We present various tests with and without magnetization and preliminary results both from single stars and from the merger of a binary system.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Neutron star-black hole mergers with a nuclear equation of state and neutrino cooling: Dependence in the binary parameters

We present a first exploration of the results of neutron star-black hole mergers using black hole masses in the most likely range of $7M_\odot-10M_\odot$, a neutrino leakage scheme, and a modeling of the neutron star material through a finite-temperature nuclear-theory based equation of state. In the range of black hole spins in which the neutron star is tidally disrupted ($χ_{\rm BH}\gtrsim 0.7$), we show that the merger consistently produces large amounts of cool ($T\lesssim 1\,{\rm MeV}$), unbound, neutron-rich material ($M_{\rm ej}\sim 0.05M_\odot-0.20M_\odot$). A comparable amount of bound matter is initially divided between a hot disk ($T_{\rm max}\sim 15\,{\rm MeV}$) with typical neutrino luminosity $L_ν\sim 10^{53}\,{\rm erg/s}$, and a cooler tidal tail. After a short period of rapid protonization of the disk lasting $\sim 10\,{\rm ms}$, the accretion disk cools down under the combined effects of the fall-back of cool material from the tail, continued accretion of the hottest material onto the black hole, and neutrino emission. As the temperature decreases, the disk progressively becomes more neutron-rich, with dimmer neutrino emission. This cooling process should stop once the viscous heating in the disk (not included in our simulations) balances the cooling. These mergers of neutron star-black hole binaries with black hole masses $M_{\rm BH}\sim 7M_\odot-10M_\odot$ and black hole spins high enough for the neutron star to disrupt provide promising candidates for the production of short gamma-ray bursts, of bright infrared post-merger signals due to the radioactive decay of unbound material, and of large amounts of r-process nuclei.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Arduous Journey to Black-Hole Formation in Potential Gamma-Ray Burst Progenitors

We present a quantitative study on the properties at death of fast-rotating massive stars evolved at low-metallicity, objects that are proposed as likely progenitors of long-duration gamma-ray bursts (LGRBs). We perform 1D+rotation stellar-collapse simulations on the progenitor models of Woosley & Heger (2006) and critically assess their potential for the formation of a black hole and a Keplerian disk (namely a collapsar) or a proto-magnetar. We note that theoretical uncertainties in the treatment of magnetic fields and the approximate handling of rotation compromises the accuracy of stellar-evolution models. We find that only the fastest rotating progenitors achieve sufficient compactness for black-hole formation while the bulk of models possess a core density structure typical of garden-variety core-collapse supernova (SN) progenitors evolved without rotation and at solar metallicity. Of the models that do have sufficient compactness for black-hole formation, most of them also retain a large amount of angular momentum in the core, making them prone to a magneto-rotational explosion, therefore preferentially leaving behind a proto-magnetar. A large progenitor angular-momentum budget is often the sole criterion invoked in the community today to assess the suitability for producing a collapsar. This simplification ignores equally important considerations such as the core compactness, which conditions black-hole formation, the core angular momentum, which may foster a magneto-rotational explosion preventing black-hole formation, or the metallicity and the residual envelope mass which must be compatible with inferences from observed LGRB/SNe. Our study suggests that black-hole formation is non trivial, that there is room for accommodating both collapsars and proto-magnetars as LGRB progenitors, although proto-magnetars seem much more easily produced by current stellar-evolutionary models.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Progenitor Dependence of the Preexplosion Neutrino Emission in Core-Collapse Supernovae

We perform spherically-symmetric general-relativistic simulations of core collapse and the postbounce preexplosion phase in 32 presupernova stellar models of solar metallicity with zero-age-main-sequence masses of 12 M_{sun} to 120 M_{sun}. Using energy-dependent three-species neutrino transport in the two-moment approximation with an analytic closure, we show that the emitted neutrino luminosities and spectra follow very systematic trends that are correlated with the compactness (~M/R) of the progenitor star's inner regions via the accretion rate in the preexplosion phase. We find that these qualitative trends depend only weakly on the nuclear equation of state, but quantitative observational statements will require independent constraints on the equation of state and the rotation rate of the core as well as a more complete understanding of neutrino oscillations. We investigate the simulated response of water Cherenkov detectors to the electron antineutrino fluxes from our models and find that the large statistics of a galactic core collapse event may allow robust conclusions on the inner structure of the progenitor star.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A New Monte Carlo Method for Time-Dependent Neutrino Radiation Transport

Monte Carlo approaches to radiation transport have several attractive properties such as simplicity of implementation, high accuracy, and good parallel scaling. Moreover, Monte Carlo methods can handle complicated geometries and are relatively easy to extend to multiple spatial dimensions, which makes them potentially interesting in modeling complex multi-dimensional astrophysical phenomena such as core-collapse supernovae. The aim of this paper is to explore Monte Carlo methods for modeling neutrino transport in core-collapse supernovae. We generalize the Implicit Monte Carlo photon transport scheme of Fleck & Cummings and gray discrete-diffusion scheme of Densmore et al. to energy-, time-, and velocity-dependent neutrino transport. Using our 1D spherically-symmetric implementation, we show that, similar to the photon transport case, the implicit scheme enables significantly larger timesteps compared with explicit time discretization, without sacrificing accuracy, while the discrete-diffusion method leads to significant speed-ups at high optical depth. Our results suggest that a combination of spectral, velocity-dependent, Implicit Monte Carlo and discrete-diffusion Monte Carlo methods represents a robust approach for use in neutrino transport calculations in core-collapse supernovae. Our velocity-dependent scheme can easily be adapted to photon transport.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Charged current neutrino interactions in core-collapse supernovae in a virial expansion

Core-collapse supernovae may depend sensitively on charged current neutrino interactions in warm, low density neutron rich matter. A proton in neutron rich matter is more tightly bound than is a neutron. This energy shift ΔU increases the electron energy in ν_e + n --> p + e, increasing the available phase space and absorption cross section. Likewise ΔU decreases the positron energy in \bar ν_e + p --> n + e^+, decreasing the phase space and cross section. We have calculated ΔU using a model independent virial expansion and we find ΔU is much larger, at low densities, than the predictions of many mean field models. Therefore ΔU could have a significant impact on charged current neutrino interactions in supernovae. Preliminary simulations of the accretion phase of core-collapse supernovae find that ΔU increases \bar ν_e energies and decreases the ν_e luminosity.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Black Hole Formation in Failing Core-Collapse Supernovae

We present results of a systematic study of failing core-collapse supernovae and the formation of stellar-mass black holes (BHs). Using our open-source general-relativistic 1.5D code GR1D equipped with a three-species neutrino leakage/heating scheme and over 100 presupernova models, we study the effects of the choice of nuclear equation of state (EOS), zero-age main sequence (ZAMS) mass and metallicity, rotation, and mass-loss prescription on BH formation. We find that the outcome, for a given EOS, can be estimated, to first order, by a single parameter, the compactness of the stellar core at bounce. By comparing protoneutron star (PNS) structure at the onset of gravitational instability with solutions of the Tolman-Oppenheimer-Volkof equations, we find that thermal pressure support in the outer PNS core is responsible for raising the maximum PNS mass by up to 25% above the cold NS value. By artificially increasing neutrino heating, we find the critical neutrino heating efficiency required for exploding a given progenitor structure and connect these findings with ZAMS conditions, establishing, albeit approximately, for the first time based on actual collapse simulations, the mapping between ZAMS parameters and the outcome of core collapse. We also study the effect of progenitor rotation and find that the dimensionless spin of nascent BHs may be robustly limited below a^* = Jc/GM^2 = 1 by the appearance of nonaxisymmetric rotational instabilities.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Thermal Effects on Black Hole Formation in Failed Core-Collapse Supernovae

We investigate several aspects of black hole formation in failing core-collapse supernovae using 1D general-relativistic hydrodynamic simulations. We use the open-source code GR1D and incorporate into it nucleon-nucleon Bremsstrahlung, a crucial neutrino pair-production channel. We focus on how various thermal effects can influence the postbounce supernova evolution towards black hole formation. By performing simulations with and without nucleon-nucleon Bremsstrahlung, we investigate the sensitivity of black hole formation to thermal support in the protoneutron star. We also investigate delayed black hole formation by artificially driving explosions in an extreme model where the protoneutron star is initially thermally supported above the maximum baryonic cold neutron star mass but then collapses to a black hole after neutrino cooling removes sufficient thermal support.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A New Open-Source Code for Spherically-Symmetric Stellar Collapse to Neutron Stars and Black Holes

We present the new open-source spherically-symmetric general-relativistic (GR) hydrodynamics code GR1D. It is based on the Eulerian formulation of GR hydrodynamics (GRHD) put forth by Romero-Ibanez-Gourgoulhon and employs radial-gauge, polar-slicing coordinates in which the 3+1 equations simplify substantially. We discretize the GRHD equations with a finite-volume scheme, employing piecewise-parabolic reconstruction and an approximate Riemann solver. GR1D is intended for the simulation of stellar collapse to neutron stars and black holes and will also serve as a testbed for modeling technology to be incorporated in multi-D GR codes. Its GRHD part is coupled to various finite-temperature microphysical equations of state in tabulated form that we make available with GR1D. An approximate deleptonization scheme for the collapse phase and a neutrino-leakage/heating scheme for the postbounce epoch are included and described. We also derive the equations for effective rotation in 1D and implement them in GR1D. We present an array of standard test calculations and also show how simple analytic equations of state in combination with presupernova models from stellar evolutionary calculations can be used to study qualitative aspects of black hole formation in failing rotating core-collapse supernovae. In addition, we present a simulation with microphysical EOS and neutrino leakage/heating of a failing core-collapse supernova and black hole formation in a presupernova model of a 40 solar mass zero-age main-sequence star. We find good agreement on the time of black hole formation (within 20%) and last stable protoneutron star mass (within 10%) with predictions from simulations with full Boltzmann neutrino radiation hydrodynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Studies of Stellar Collapse and Black Hole Formation with the Open-Source Code GR1D

We discuss results from simulations of black hole formation in failing core-collapse supernovae performed with the code GR1D, a new open-source Eulerian spherically-symmetric general-relativistic hydrodynamics code. GR1D includes rotation in an approximate way (1.5D), comes with multiple finite-temperature nuclear equations of state (EOS), and treats neutrinos in the post-core-bounce phase via a 3-flavor leakage scheme and a heating prescription. We chose the favored K_0=220 MeV-variant of the Lattimer & Swesty (1990) EOS and present collapse calculations using the progenitor models of Limongi & Chieffi (2006). We show that there is no direct (or ``prompt'') black hole formation in the collapse of ordinary massive stars (8 M_Sun ~< M_ZAMS ~< 100 M_Sun) and present first results from black hole formation simulations that include rotation.

0 citations0 reviewsNo code linkedOpen access