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

Josiah Schwab

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

ResearcherUnclaimed source record
astro-ph.SRastro-ph.HEastro-ph.COastro-ph.GAastro-ph.IMgr-qcnucl-thphysics.flu-dyn

Catalog footprint

What is connected

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

14 published item(s)

preprint2021arXiv

Skye: A Differentiable Equation of State

Stellar evolution and numerical hydrodynamics simulations depend critically on access to fast, accurate, thermodynamically consistent equations of state. We present Skye, a new equation of state for fully-ionized matter. Skye includes the effects of positrons, relativity, electron degeneracy, Coulomb interactions, non-linear mixing effects, and quantum corrections. Skye determines the point of Coulomb crystallization in a self-consistent manner, accounting for mixing and composition effects automatically. A defining feature of this equation of state is that it uses analytic free energy terms and provides thermodynamic quantities using automatic differentiation machinery. Because of this, Skye is easily extended to include new effects by simply writing new terms in the free energy. We also introduce a novel thermodynamic extrapolation scheme for extending analytic fits to the free energy beyond the range of the fitting data while preserving desirable properties like positive entropy and sound speed. We demonstrate Skye in action in the MESA stellar evolution software instrument by computing white dwarf cooling curves.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Evolutionary Models for the Remnant of the Merger of Two Carbon-Oxygen Core White Dwarfs

We construct evolutionary models of the remnant of the merger of two carbon-oxygen (CO) core white dwarfs (WDs). With total masses in the range $1-2 {\rm M_\odot}$, these remnants may either leave behind a single massive WD or undergo a merger-induced collapse to a neutron star (NS). On the way to their final fate, these objects generally experience a $\sim 10$ kyr luminous giant phase, which may be extended if sufficient helium remains to set up a stable shell-burning configuration. The uncertain, but likely significant, mass loss rate during this phase influences the final remnant mass and fate (WD or NS). We find that the initial CO core composition of the WD is converted to oxygen-neon (ONe) in remnants with final masses $\gtrsim 1.05 {\rm M_\odot}$. This implies that the CO core / ONe core transition in single WDs formed via mergers occurs at a similar mass as in WDs descended from single stars, and thus that WD-WD mergers do not naturally provide a route to producing ultra-massive CO-core WDs. As the remnant contracts towards a compact configuration, it experiences a "bottleneck" that sets the characteristic total angular momentum that can be retained. This limit predicts single WDs formed from WD-WD mergers have rotational periods of $\approx 10-20$ min on the WD cooling track. Similarly, it predicts remnants that collapse can form NSs with rotational periods $\sim 10$ ms.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Laminar Flame Speeds in Degenerate Oxygen-Neon Mixtures

The collapse of degenerate oxygen-neon cores (i.e., electron-capture supernovae or accretion-induced collapse) proceeds through a phase in which a deflagration wave ("flame") forms at or near the center and propagates through the star. In models, the assumed speed of this flame influences whether this process leads to an explosion or to the formation of a neutron star. We calculate the laminar flame speeds in degenerate oxygen-neon mixtures with compositions motivated by detailed stellar evolution models. These mixtures include trace amounts of carbon and have a lower electron fraction than those considered in previous work. We find that trace carbon has little effect on the flame speeds, but that material with electron fraction $Y_e \approx 0.48-0.49$ has laminar flame speeds that are $\approx 2$ times faster than those at $Y_e = 0.5$. We provide tabulated flame speeds and a corresponding fitting function so that the impact of this difference can be assessed via full star hydrodynamical simulations of the collapse process.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Updated Parameter Estimates for GW190425 Using Astrophysical Arguments and Implications for the Electromagnetic Counterpart

The progenitor system of the compact binary merger GW190425 had a total mass of $3.4^{+0.3}_{-0.1}$ M$_\odot$ (90th-percentile confidence region) as measured from its gravitational wave signal. This mass is significantly different from the Milky Way (MW) population of binary neutron stars (BNSs) that are expected to merge in a Hubble time and from that of the first BNS merger, GW170817. Here we explore the expected electromagnetic signatures of such a system. We make several astrophysically motivated assumptions to further constrain the parameters of GW190425. By simply assuming that both components were NSs, we reduce the possible component masses significantly, finding $m_1 = 1.85^{+0.27}_{-0.19}$ M$_\odot$ and $m_2 = 1.47^{+0.16}_{-0.18}$ M$_\odot$. However if the GW190425 progenitor system was a NS-black hole merger, we find best-fitting parameters $m_1 = 2.19^{+0.21}_{-0.17}$ M$_\odot$ and $m_2 = 1.26^{+0.10}_{-0.08}$ M$_\odot$. For a well-motivated BNS system where the lighter NS has a mass similar to the mass of non-recycled NSs in MW BNS systems, we find $m_1 = 2.03^{+0.15}_{-0.14}$ M$_\odot$ and $m_2 = 1.35 \pm 0.09$ M$_\odot$, corresponding to only 7% mass uncertainties. For all scenarios, we expect a prompt collapse of the resulting remnant to a black hole. Examining detailed models with component masses similar to our best-fitting results, we find the electromagnetic counterpart to GW190425 is expected to be significantly redder and fainter than that of GW170817. We find that almost all reported observations used to search for an electromagnetic counterpart for GW190425 were too shallow to detect the expected counterpart. If the LIGO-Virgo Collaboration promptly provides the chirp mass, the astronomical community can adapt their observations to improve the likelihood of detecting a counterpart for similarly "high-mass" BNS systems. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Carbon Shell or Core Ignitions in White Dwarfs Accreting from Helium Stars

White dwarfs accreting from helium stars can stably burn at the accreted rate and avoid the challenge of mass loss associated with unstable Helium burning that is a concern for many Type Ia supernovae scenarios. We study binaries with helium stars of mass $1.25 M_\odot\le M_{\rm{He}} \le 1.8 M_\odot$, which have lost their hydrogen rich envelopes in an earlier common envelope event and now orbit with periods ($P_{\rm orb}$) of several hours with non-rotating $0.84$ and $1.0 M_\odot$ C/O WDs. The helium stars fill their Roche lobes (RLs) after exhaustion of central helium and donate helium on their thermal timescales (${\sim}10^5$yr). As shown by others, these mass transfer rates coincide with the steady helium burning range for WDs, and grow the WD core up to near the Chandrasekhar mass ($M_{\rm Ch}$) and a core carbon ignition. We show here, however, that many of these scenarios lead to an ignition of hot carbon ashes near the outer edge of the WD and an inward going carbon flame that does not cause an explosive outcome. For $P_{\rm orb} = 3$ hours, $1.0 M_\odot$ C/O WDs with donor masses $M_{\rm He}\gtrsim1.8 M_\odot$ experience a shell carbon ignition, while $M_{\rm He}\lesssim1.3 M_\odot$ will fall below the steady helium burning range and undergo helium flashes before reaching core C ignition. Those with $1.3 M_\odot \lesssim M_{\rm He} \lesssim 1.7 M_\odot$ will experience a core C ignition. We also calculate the retention fraction of accreted helium when the accretion rate leads to recurrent weak helium flashes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Neutronization During Carbon Simmering In Type Ia Supernova Progenitors

When a Type Ia supernova (SN Ia) progenitor first ignites carbon in its core, it undergoes ${\sim} \,10^{3}-10^{4} \,$yr of convective burning prior to the onset of thermonuclear runaway. This carbon simmering phase is important for setting the thermal profile and composition of the white dwarf. Using the \texttt{MESA} stellar evolution code, we follow this convective burning and examine the production of neutron-rich isotopes. The neutron content of the SN fuel has important consequences for the ensuing nucleosynthesis, and, in particular, for the production of secondary Fe-peak nuclei like Mn and stable Ni. These elements have been observed in the X-ray spectra of SN remnants like Tycho, Kepler, and 3C 397, and their yields can provide valuable insights into the physics of SNe Ia and the properties of their progenitors. We find that weak reactions during simmering can at most generate a neutron excess of ${\approx} \, 3 \times 10^{-4}$. This is ${\approx} \, 8 \times 10^{-4}$ lower than that found in previous studies that do not take the full density and temperature profile of the simmering region into account. Our results imply that the progenitor metallicity is the main contributor to the neutron excess in SN Ia fuel for $Z \gtrsim 1/3 \, Z_{\odot}$. Alternatively, at lower metallicities, this neutron excess provides a floor that should be present in any centrally-ignited SN~Ia scenario.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Evolution and Fate of Super-Chandrasekhar Mass White Dwarf Merger Remnants

We present stellar evolution calculations of the remnant of the merger of two carbon-oxygen white dwarfs (CO WDs). We focus on cases that have a total mass in excess of the Chandrasekhar mass. After the merger, the remnant manifests as an $L \sim 3 \times 10^4 L_\odot$ source for $\sim 10^4$ yr. A dusty wind may develop, leading these sources to be self-obscured and to appear similar to extreme AGB stars. Roughly $\sim 10$ such objects should exist in the Milky Way and M31 at any time. As found in previous work, off-center carbon fusion is ignited within the merger remnant and propagates inward via a carbon flame, converting the WD to an oxygen-neon (ONe) composition. By following the evolution for longer than previous calculations, we demonstrate that after carbon-burning reaches the center, neutrino-cooled Kelvin-Helmholtz contraction leads to off-center neon ignition in remnants with masses $\ge 1.35 M_\odot$. The resulting neon-oxygen flame converts the core to a silicon WD. Thus, super-Chandrasekhar WD merger remnants do not undergo electron-capture induced collapse as traditionally assumed. Instead, if the remnant mass remains above the Chandrasekhar mass, we expect that it will form a low-mass iron core and collapse to form a neutron star. Remnants that lose sufficient mass will end up as massive, isolated ONe or Si WDs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Thermal Runaway During the Evolution of ONeMg Cores towards Accretion-Induced Collapse

We study the evolution of degenerate electron cores primarily composed of the carbon burning products oxygen, neon, and magnesium (hereafter ONeMg cores) that are undergoing compression. Electron capture reactions on A=20 and A=24 isotopes reduce the electron fraction and heat the core. We develop and use a new capability of the Modules for Experiments in Stellar Astrophysics (MESA) stellar evolution code that provides a highly accurate implementation of these key reactions. These new accurate rates and the ability of MESA to perform extremely small spatial zoning demonstrates a thermal runaway in the core triggered by the temperature and density sensitivity of the Ne-20 electron capture reactions. Both analytics and numerics show that this thermal runaway does not trigger core convection, but rather leads to a centrally concentrated (r < km) thermal runaway that will subsequently launch an oxygen deflagration wave from the center of the star. We use MESA to perform a parameter study that quantifies the influence of the magnesium mass fraction, the central temperature, the compression rate, and uncertainties in the electron capture reaction rates on the ONeMg core evolution. This allows us to establish a lower limit on the central density at which the oxygen deflagration wave initiates of $ρ_c > 8.5 \times 10^9\, \textrm{g cm}^{-3}$. Based on previous work and order-of-magnitude calculations, we expect objects which ignite oxygen at or above these densities to collapse and form a neutron star. Calculations such as these are an important step in producing more realistic progenitor models for studies of the signature of accretion-induced collapse.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Turbulent Chemical Diffusion in Convectively Bounded Carbon Flames

It has been proposed that mixing induced by convective overshoot can disrupt the inward propagation of carbon deflagrations in super-asymptotic giant branch stars. To test this theory, we study an idealized model of convectively bounded carbon flames with 3D hydrodynamic simulations of the Boussinesq equations using the pseudospectral code Dedalus. Because the flame propagation timescale is much longer than the convection timescale, we approximate the flame as fixed in space, and only consider its effects on the buoyancy of the fluid. By evolving a passive scalar field, we derive a {\it turbulent} chemical diffusivity produced by the convection as a function of height, $D_{\rm t}(z)$. Convection can stall a flame if the chemical mixing timescale, set by the turbulent chemical diffusivity, $D_{\rm t}$, is shorter than the flame propagation timescale, set by the thermal diffusivity, $κ$, i.e., when $D_{\rm t}>κ$. However, we find $D_{\rm t}<κ$ for most of the flame because convective plumes are not dense enough to penetrate into the flame. Extrapolating to realistic stellar conditions, this implies that convective mixing cannot stall a carbon flame and that "hybrid carbon-oxygen-neon" white dwarfs are not a typical product of stellar evolution.

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
preprint2013arXiv

Type Ia Supernovae from Merging White Dwarfs II) Post-Merger Detonations

Merging carbon-oxygen (CO) white dwarfs are a promising progenitor system for Type Ia supernovae (SN Ia), but the underlying physics and timing of the detonation are still debated. If an explosion occurs after the secondary star is fully disrupted, the exploding primary will expand into a dense CO medium that may still have a disk-like structure. This interaction will decelerate and distort the ejecta. Here we carry out multi-dimensional simulations of ``tamped" SN Ia models, using both particle and grid-based codes to study the merger and explosion dynamics, and a radiative transfer code to calculate synthetic spectra and light curves. We find that post-merger explosions exhibit an hourglass-shaped asymmetry, leading to strong variations in the light curves with viewing angle. The two most important factors affecting the outcome are the scale-height of the disk, which depends sensitively on the binary mass ratio, and the total ${}^{56}$Ni yield, which is governed by the central density of the remnant core. The synthetic broadband light curves rise and decline very slowly, and the spectra generally look peculiar, with weak features from intermediate mass elements but relatively strong carbon absorption. We also consider the effects of the viscous evolution of the remnant, and show that a longer time delay between merger and explosion probably leads to larger ${}^{56}$Ni yields and more symmetrical remnants. We discuss the relevance of this class of aspherical ``tamped" SN Ia for explaining the class of ``super-Chandrasekhar'' SN Ia.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Viscous Evolution of White Dwarf Merger Remnants

The merger of two white dwarfs (WDs) creates a differentially rotating remnant which is unstable to magnetohydrodynamic instabilities. These instabilities can lead to viscous evolution on a time-scale short compared to the thermal evolution of the remnant. We present multi-dimensional hydrodynamic simulations of the evolution of WD merger remnants under the action of an $α$-viscosity. We initialize our calculations using the output of eight WD merger simulations from Dan et al. (2011), which span a range of mass ratios and total masses. We generically find that the merger remnants evolve towards spherical states on time-scales of hours, even though a significant fraction of the mass is initially rotationally supported. The viscous evolution unbinds only a very small amount of mass $(< 10^{-5} M_\odot)$. Viscous heating causes some of the systems we study with He WD secondaries to reach conditions of nearly dynamical burning. It is thus possible that the post-merger viscous phase triggers detonation of the He envelope in some WD mergers, potentially producing a Type Ia supernova via a double detonation scenario. Our calculations provide the proper initial conditions for studying the long-term thermal evolution of WD merger remnants. This is important for understanding WD mergers as progenitors of Type Ia supernovae, neutron stars, R Coronae Borealis stars and other phenomena.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Further Evidence for the Bimodal Distribution of Neutron Star Masses

We use a collection of 14 well-measured neutron star masses to strengthen the case that a substantial fraction of these neutron stars was formed via electron-capture supernovae (SNe) as opposed to Fe-core collapse SNe. The e-capture SNe are characterized by lower resultant gravitational masses and smaller natal kicks, leading to lower orbital eccentricities when the e-capture SN has led to the formation of the second neutron star in a binary system. Based on the measured masses and eccentricities, we identify four neutron stars, which have a mean post-collapse gravitational mass of ~1.25 solar masses, as the product of e-capture SNe. We associate the remaining ten neutron stars, which have a mean mass of 1.35 solar masses, with Fe-core collapse SNe. If the e-capture supernova occurs during the formation of the first neutron star, then this should substantially increase the formation probability for double neutron stars, given that more systems will remain bound with the smaller kicks. However, this does not appear to be the case for any of the observed systems, and we discuss possible reasons for this.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Galaxy-Scale Strong Lensing Tests of Gravity and Geometric Cosmology: Constraints and Systematic Limitations

Galaxy-scale strong gravitational lenses with measured stellar velocity dispersions allow a test of the weak-field metric on kiloparsec scales and a geometric measurement of the cosmological distance-redshift relation, provided that the mass-dynamical structure of the lensing galaxies can be independently constrained to a sufficient degree. We combine data on 53 galaxy-scale strong lenses from the Sloan Lens ACS Survey with a well-motivated fiducial set of lens-galaxy parameters to find (1) a constraint on the post-Newtonian parameter gamma = 1.01 +/- 0.05 and (2) a determination of Omega_Lambda = 0.75 +/- 0.17 under the assumption of a flat universe. These constraints assume that the underlying observations and priors are free of systematic error. We evaluate the sensitivity of these results to systematic uncertainties in (1) total mass-profile shape, (2) velocity anisotropy, (3) light-profile shape, and (4) stellar velocity dispersion. Based on these sensitivities, we conclude that while such strong-lens samples can in principle provide an important tool for testing general relativity and cosmology, they are unlikely to yield precision measurements of gamma and Omega_Lambda unless the properties of the lensing galaxies are independently constrained with substantially greater accuracy than at present.

0 citations0 reviewsNo code linkedOpen access