Source author record

A. Heger

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

31works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

31 published item(s)

preprint2020arXiv

Properties of gamma-ray decay lines in 3D core-collapse supernova models, with application to SN 1987A and Cas A

Comparison of theoretical line profiles to observations provides important tests for supernova explosion models. We study the shapes of radioactive decay lines predicted by current 3D core-collapse explosion simulations, and compare these to observations of SN 1987A and Cas A. Both the widths and shifts of decay lines vary by several thousand kilometers per second depending on viewing angle. The line profiles can be complex with multiple peaks. By combining observational constraints from 56Co decay lines, 44Ti decay lines, and Fe IR lines, we delineate a picture of the morphology of the explosive burning ashes in SN 1987A. For M_ZAMS=15-20 Msun progenitors exploding with ~1.5 *10^51 erg, ejecta structures suitable to reproduce the observations involve a bulk asymmetry of the 56Ni of at least ~400 km/s and a bulk velocity of at least ~1500 km/s. By adding constraints to reproduce the UVOIR bolometric light curve of SN 1987A up to 600d, an ejecta mass around 14 Msun is favoured. We also investigate whether observed decay lines can constrain the neutron star (NS) kick velocity. The model grid provides a constraint V_NS > V_redshift, and applying this to SN 1987A gives a NS kick of at least 500 km/s. For Cas A, our single model provides a satisfactory fit to the NuSTAR observations and reinforces the result that current neutrino-driven core-collapse SN models can achieve enough bulk asymmetry in the explosive burning material. Finally, we investigate the internal gamma-ray field and energy deposition, and compare the 3D models to 1D approximations.

preprint2020arXiv

The Formation of a 70 Msun Black Hole at High Metallicity

A 70Msun BH was discovered in Milky Way disk in a long period and almost circular detached binary system (LB-1) with a high metallicity 8Msun B star companion. Current consensus on the formation of BHs from high metallicity stars limits the black hole mass to be below 20Msun. Using simple evolutionary model, we show that the formation of a 70Msun BH in high metallicity environment is possible if stellar wind mass loss rates are reduced by factor of 5. As observations indicate, a fraction of massive stars have surface magnetic fields which may quench the wind mass-loss, independently of stellar mass and metallicity. We also computed detailed stellar evolution models and we confirm such a scenario. A non-rotating 85Msun model at Z=0.014 with decreased winds ends up as a 71Msun star prior core-collapse with a 32Msun helium core and a 28Msun CO core. Such star avoids pair-instability pulsation supernova mass loss and may form a 70Msun BH in the direct collapse. Stars that can form such BHs expand to significant size with radius of R>600Rsun, exceeding the size of LB-1 orbit. Therefore, we can explain the formation of BHs upto 70Msun at high metallicity and this result is independent from LB-1. However, if LB-1 hosts a massive BH we are unable to explain how such a binary star system could have formed without invoking some exotic scenarios.

preprint2020arXiv

Three-dimensional Models of Core-collapse Supernovae From Low-mass Progenitors With Implications for Crab

We present 3D full-sphere supernova simulations of non-rotating low-mass (~9 Msun) progenitors, covering the entire evolution from core collapse through bounce and shock revival, through shock breakout from the stellar surface, until fallback is completed several days later. We obtain low-energy explosions [~(0.5-1.0)x 10^{50} erg] of iron-core progenitors at the low-mass end of the core-collapse supernova (LMCCSN) domain and compare to a super-AGB (sAGB) progenitor with an oxygen-neon-magnesium core that collapses and explodes as electron-capture supernova (ECSN). The onset of the explosion in the LMCCSN models is modelled self-consistently using the Vertex-Prometheus code, whereas the ECSN explosion is modelled using parametric neutrino transport in the Prometheus-HOTB code, choosing different explosion energies in the range of previous self-consistent models. The sAGB and LMCCSN progenitors that share structural similarities have almost spherical explosions with little metal mixing into the hydrogen envelope. A LMCCSN with less 2nd dredge-up results in a highly asymmetric explosion. It shows efficient mixing and dramatic shock deceleration in the extended hydrogen envelope. Both properties allow fast nickel plumes to catch up with the shock, leading to extreme shock deformation and aspherical shock breakout. Fallback masses of <~5x10^{-3} Msun have no significant effects on the neutron star (NS) masses and kicks. The anisotropic fallback carries considerable angular momentum, however, and determines the spin of the newly-born NS. The LMCCSNe model with less 2nd dredge-up results in a hydrodynamic and neutrino-induced NS kick of >40 km/s and a NS spin period of ~30 ms, both not largely different from those of the Crab pulsar at birth.

preprint2019arXiv

Gamma-ray Emission of 60Fe and 26Al Radioactivities in our Galaxy

The isotopes $^{60}$Fe and $^{26}$Al originate from massive stars and their supernovae, reflecting ongoing nucleosynthesis in the Galaxy. We studied the gamma-ray emission from these isotopes at characteristic energies 1173, 1332, and 1809 keV with over 15 years of SPI data, finding a line flux in $^{60}$Fe combined lines of $(0.31\pm 0.06) \times 10^{-3}$ ph cm$^{-2}$ s$^{-1}$ and the $^{26}$Al line flux of $(16.8\pm 0.7) \times 10^{-4}$ ph cm$^{-2}$ s$^{-1}$ above the background and continuum emission for the whole sky. Based on the exponential-disk grid maps, we characterise the emission extent of $^{26}$Al to find scale parameters $R_0 =7.0^{+1.5}_{-1.0}$ kpc and $z_0=0.8^{+0.3}_{-0.2}$ kpc, however the $^{60}$Fe lines are too weak to spatially constrain the emission. Based on a point source model test across the Galactic plane, the $^{60}$Fe emission would not be consistent with a single strong point source in the Galactic center or somewhere else, providing a hint for a diffuse nature. We carried out comparisons of emission morphology maps using different candidate-source tracers for both $^{26}$Al and $^{60}$Fe emissions, and suggests that the $^{60}$Fe emission is more likely to be concentrated towards the Galactic plane. We determine the $^{60}$Fe/$^{26}$Al $γ$-ray flux ratio at $(18.4\pm4.2)\,\%$ , when using a parameterized spatial morphology model. Across the range of plausible morphologies, it appears possible that $^{26}$Al and $^{60}$Fe are distributed differently in the Galaxy. Using the best fitting maps for each of the elements, we constrain flux ratios in the range 0.2--0.4. We discuss its implications for massive star models and their nucleosynthesis.

preprint2016arXiv

A Simple Approach to the Supernova Progenitor-Explosion Connection

We present a new approach to understand the landscape of supernova explosion energies, ejected nickel masses, and neutron star birth masses. In contrast to other recent parametric approaches, our model predicts the properties of neutrino-driven explosions based on the pre-collapse stellar structure without the need for hydrodynamic simulations. The model is based on physically motivated scaling laws and simple differential equations describing the shock propagation, the contraction of the neutron star, the neutrino emission, the heating conditions, and the explosion energetics. Using model parameters compatible with multi-D simulations and a fine grid of thousands of supernova progenitors, we obtain a variegated landscape of neutron star and black hole formation similar to other parameterised approaches and find good agreement with semi-empirical measures for the "explodability" of massive stars. Our predicted explosion properties largely conform to observed correlations between the nickel mass and explosion energy. Accounting for the coexistence of outflows and downflows during the explosion phase, we naturally obtain a positive correlation between explosion energy and ejecta mass. These correlations are relatively robust against parameter variations, but our results suggest that there is considerable leeway in parametric models to widen or narrow the mass ranges for black hole and neutron star formation and to scale explosion energies up or down. Our model is currently limited to an all-or-nothing treatment of fallback and there remain some minor discrepancies between model predictions and observational constraints.

preprint2016arXiv

Dependence of X-Ray Burst Models on Nuclear Reaction Rates

X-ray bursts are thermonuclear flashes on the surface of accreting neutron stars and reliable burst models are needed to interpret observations in terms of properties of the neutron star and the binary system. We investigate the dependence of X-ray burst models on uncertainties in (p,$γ$), ($α$,$γ$), and ($α$,p) nuclear reaction rates using fully self-consistent burst models that account for the feedbacks between changes in nuclear energy generation and changes in astrophysical conditions. A two-step approach first identified sensitive nuclear reaction rates in a single-zone model with ignition conditions chosen to match calculations with a state-of-the-art 1D multi-zone model based on the {\Kepler} stellar evolution code. All relevant reaction rates on neutron deficient isotopes up to mass 106 were individually varied by a factor of 100 up and down. Calculations of the 84 highest impact reaction rate changes were then repeated in the 1D multi-zone model. We find a number of uncertain reaction rates that affect predictions of light curves and burst ashes significantly. The results provide insights into the nuclear processes that shape X-ray burst observables and guidance for future nuclear physics work to reduce nuclear uncertainties in X-ray burst models.

preprint2016arXiv

NuGrid stellar data set I. Stellar yields from H to Bi for stars with metallicities Z = 0.02 and Z = 0.01

We provide a set of stellar evolution and nucleosynthesis calculations that applies established physics assumptions simultaneously to low- and intermediate-mass and massive star models. Our goal is to provide an internally consistent and comprehensive nuclear production and yield data base for applications in areas such as pre-solar grain studies. Our non-rotating models assume convective boundary mixing where it has been adopted before. We include 8 (12) initial masses for $Z = 0.01$ ($0.02$). Models are followed either until the end of the asymptotic giant branch phase or the end of Si burning, complemented by a simple analytic core-collapse supernova models with two options for fallback and shock velocities. The explosions show which pre-supernova yields will most strongly be effected by the explosive nucleosynthesis. We discuss how these two explosion parameters impacts the light elements and the $s$ and $p$ process. For low- and intermediate-mass models our stellar yields from H to Bi include the effect of convective boundary mixing at the He-intershell boundaries and the stellar evolution feedback of the mixing process that produces the $^{13}$C pocket. All post-processing nucleosynthesis calculations use the same nuclear reaction rate network and nuclear physics input. We provide a discussion of the nuclear production across the entire mass range organized by element group. All our stellar nucleosynthesis profile and time evolution output is available electronically, and tools to explore the data on the NuGrid VOspace hosted by the Canadian Astronomical Data Centre are introduced.

preprint2016arXiv

The Effect of Pair-Instability Mass Loss on Black Hole Mergers

Mergers of two stellar origin black holes are a prime source of gravitational waves and are under intensive investigations. One crucial ingredient in their modeling has so far been neglected. Pair-instability pulsation supernovae with associated severe mass loss may suppress formation of massive black holes, decreasing black hole merger rates for the highest black hole masses. The mass loss associated with pair-instability pulsation supernovae limits the Population I/II stellar-origin black hole mass to 50 Msun, in tension with earlier predictions that the maximum black hole mass could be as high as 100 Msun. Suppression of double black hole merger rates by pair-instability pulsation supernovae is negligible for our evolutionary channel. Our standard evolutionary model with inclusion of pair-instability pulsation supernovae and with pair-instability supernovae is fully consistent with the LIGO detections of black hole mergers: GW150914, GW151226 and LVT151012. The LIGO observations seem to exclude high (>400 km/s) BH natal kicks. We predict the detection of several up to 60 BH-BH mergers with total mass 10--150 Msun (most likely range: 20--80 Msun) in the forthcoming 60 effective days of LIGO O2 observations.

preprint2016arXiv

Uncertainties in the production of $p$ nuclei in massive stars obtained from Monte Carlo variations

Nuclear uncertainties in the production of $p$ nuclei in massive stars have been quantified in a Monte Carlo procedure. Bespoke temperature-dependent uncertainties were assigned to different types of reactions involving nuclei from Fe to Bi. Their simultaneous impact was studied in postprocessing explosive trajectories for three different stellar models. It was found that the grid of mass zones in the model of a 25 $M_\odot$ star, which is widely used for investigations of $p$ nucleosynthesis, is too crude to properly resolve the detailed temperature changes required for describing the production of $p$ nuclei. Using models with finer grids for 15 $M_\odot$ and 25 $M_\odot$ stars with initial solar metallicity, it was found that most of the production uncertainties introduced by nuclear reaction uncertainties are smaller than a factor of two. Since a large number of rates were varied at the same time in the Monte Carlo procedure, possible cancellation effects of several uncertainties could be taken into account. Key rates were identified for each $p$ nucleus, which provide the dominant contribution to the production uncertainty. These key rates were found by examining correlations between rate variations and resulting abundance changes. This method is superior to studying flow patterns, especially when the flows are complex, and to individual, sequential variation of a few rates.

preprint2015arXiv

Carbon production on accreting neutron stars in a new regime of stable nuclear burning

Accreting neutron stars exhibit Type I X-ray bursts from both frequent hydrogen/helium flashes as well as rare carbon flashes. The latter (superbursts) ignite in the ashes of the former. Hydrogen/helium bursts, however, are thought to produce insufficient carbon to power superbursts. Stable burning could create the required carbon, but this was predicted to only occur at much larger accretion rates than where superbursts are observed. We present models of a new steady-state regime of stable hydrogen and helium burning that produces pure carbon ashes. Hot CNO burning of hydrogen heats the neutron star envelope and causes helium to burn before the conditions of a helium flash are reached. This takes place when the mass accretion rate is around 10% of the Eddington limit: close to the rate where most superbursts occur. We find that increased heating at the base of the envelope sustains steady-state burning by steepening the temperature profile, which increases the amount of helium that burns before a runaway can ensue.

preprint2015arXiv

Constraints on explosive silicon burning in core-collapse supernovae from measured Ni/Fe ratios

Measurements of explosive nucleosynthesis yields in core-collapse supernovae provide tests for explosion models. We investigate constraints on explosive conditions derivable from measured amounts of nickel and iron after radioactive decays using nucleosynthesis networks with parameterized thermodynamic trajectories. The Ni/Fe ratio is for most regimes dominated by the production ratio of 58Ni/(54Fe + 56Ni), which tends to grow with higher neutron excess and with higher entropy. For SN 2012ec, a supernova that produced a Ni/Fe ratio of $3.4\pm1.2$ times solar, we find that burning of a fuel with neutron excess $η\approx 6\times 10^{-3}$ is required. Unless the progenitor metallicity is over 5 times solar, the only layer in the progenitor with such a neutron excess is the silicon shell. Supernovae producing large amounts of stable nickel thus suggest that this deep-lying layer can be, at least partially, ejected in the explosion. We find that common spherically symmetric models of $M_{\rm ZAMS} \lesssim 13$ Msun stars exploding with a delay time of less than one second ($M_{\rm cut} < 1.5$ Msun) are able to achieve such silicon-shell ejection. Supernovae that produce solar or sub-solar Ni/Fe ratios, such as SN 1987A, must instead have burnt and ejected only oxygen-shell material, which allows a lower limit to the mass cut to be set. Finally, we find that the extreme Ni/Fe value of 60-75 times solar derived for the Crab cannot be reproduced by any realistic-entropy burning outside the iron core, and neutrino-neutronization obtained in electron-capture models remains the only viable explanation.

preprint2015arXiv

Nebular spectra of pair-instability supernovae

If very massive stars (M >~ 100 Msun) can form and avoid too strong mass loss during their evolution, they are predicted to explode as pair-instability supernovae (PISNe). One critical test for candidate events is whether their nucleosynthesis yields and internal ejecta structure, being revealed through nebular-phase spectra at t >~ 1 yr, match those of model predictions. Here we compute theoretical spectra based on model PISN ejecta at 1-3 years post-explosion to allow quantitative comparison with observations. The high column densities of PISNe lead to complete gamma-ray trapping for t >~ 2 years which, combined with fulfilled conditions of steady state, leads to bolometric supernova luminosities matching the 56Co decay. Most of the gamma-rays are absorbed by the deep-lying iron and silicon/sulphur layers. The ionization balance shows a predominantly neutral gas state, which leads to emission lines of Fe I, Si I, and S I. For low-mass PISNe the metal core expands slowly enough to produce a forest of distinct lines, whereas high-mass PISNe expand faster and produce more featureless spectra. Line blocking is complete below ~5000 A for several years, and the model spectra are red. The strongest line is typically [Ca II] 7291,7323, one of few lines from ionized species. We compare our models with proposed PISN candidates SN 2007bi and PTF12dam, finding discrepancies for several key observables and thus no support for a PISN interpretation. We discuss distinct spectral features predicted by the models, and the possibility of detecting pair-instability explosions among non-superluminous supernovae.

preprint2015arXiv

The first direct measurement of 12C(12C,n)23Mg at stellar energies

Neutrons produced by the carbon fusion reaction 12C(12C,n)23Mg play an important role in stellar nucleosynthesis. However, past studies have shown large discrepancies between experimental data and theory, leading to an uncertain cross section extrapolation at astrophysical energies. We present the first direct measurement that extends deep into the astrophysical energy range along with a new and improved extrapolation technique based on experimental data from the mirror reaction 12C(12C,p)23Na. The new reaction rate has been determined with a well-defined uncertainty that exceeds the precision required by astrophysics models. Using our constrained rate, we find that 12C(12C,n)23Mg is crucial to the production of Na and Al in Pop-III Pair Instability Supernovae. It also plays a non-negligible role in the production of weak s-process elements as well as in the production of the important galactic gamma-ray emitter 60Fe.

preprint2015arXiv

The GALAH Survey: Scientific Motivation

The GALAH survey is a large high-resolution spectroscopic survey using the newly commissioned HERMES spectrograph on the Anglo-Australian Telescope. The HERMES spectrograph provides high-resolution (R ~28,000) spectra in four passbands for 392 stars simultaneously over a 2 degree field of view. The goal of the survey is to unravel the formation and evolutionary history of the Milky Way, using fossil remnants of ancient star formation events which have been disrupted and are now dispersed throughout the Galaxy. Chemical tagging seeks to identify such dispersed remnants solely from their common and unique chemical signatures; these groups are unidentifiable from their spatial, photometric or kinematic properties. To carry out chemical tagging, the GALAH survey will acquire spectra for a million stars down to V~14. The HERMES spectra of FGK stars contain absorption lines from 29 elements including light proton-capture elements, alpha-elements, odd-Z elements, iron-peak elements and n-capture elements from the light and heavy s-process and the r-process. This paper describes the motivation and planned execution of the GALAH survey, and presents some results on the first-light performance of HERMES.

preprint2014arXiv

A single low-energy, iron-poor supernova as the source of metals in the star SMSS J 031300.36-670839.3

The element abundance ratios of four low-mass stars with extremely low metallicities indicate that the gas out of which the stars formed was enriched in each case by at most a few, and potentially only one low-energy, supernova. Such supernovae yield large quantities of light elements such as carbon but very little iron. The dominance of low-energy supernovae is surprising, because it has been expected that the first stars were extremely massive, and that they disintegrated in pair-instability explosions that would rapidly enrich galaxies in iron. What has remained unclear is the yield of iron from the first supernovae, because hitherto no star is unambiguously interpreted as encapsulating the yield of a single supernova. Here we report the optical spectrum of SMSS J031300.36- 670839.3, which shows no evidence of iron (with an upper limit of 10^-7.1 times solar abundance). Based on a comparison of its abundance pattern with those of models, we conclude that the star was seeded with material from a single supernova with an original mass of ~60 Mo (and that the supernova left behind a black hole). Taken together with the previously mentioned low-metallicity stars, we conclude that low-energy supernovae were common in the early Universe, and that such supernovae yield light element enrichment with insignificant iron. Reduced stellar feedback both chemically and mechanically from low-energy supernovae would have enabled first-generation stars to form over an extended period. We speculate that such stars may perhaps have had an important role in the epoch of cosmic reionization and the chemical evolution of early galaxies.

preprint2014arXiv

Reaction Rate and Composition Dependence of the Stability of Thermonuclear Burning on Accreting Neutron Stars

The stability of thermonuclear burning of hydrogen and helium accreted onto neutron stars is strongly dependent on the mass accretion rate. The burning behavior is observed to change from Type I X-ray bursts to stable burning, with oscillatory burning occurring at the transition. Simulations predict the transition at a ten times higher mass accretion rate than observed. Using numerical models we investigate how the transition depends on the hydrogen, helium, and CNO mass fractions of the accreted material, as well as on the nuclear reaction rates of triple alpha and the hot-CNO breakout reactions 15O(a,g)19Ne and 18Ne(a,p)21Na. For a lower hydrogen content the transition is at higher accretion rates. Furthermore, most experimentally allowed reaction rate variations change the transition accretion rate by at most 10%. A factor ten decrease of the 15O(a,g)19Ne rate, however, produces an increase of the transition accretion rate of 35%. None of our models reproduce the transition at the observed rate, and depending on the true 15O(a,g)19Ne reaction rate, the actual discrepancy may be substantially larger. We find that the width of the interval of accretion rates with marginally stable burning depends strongly on both composition and reaction rates. Furthermore, close to the stability transition, our models predict that X-ray bursts have extended tails where freshly accreted fuel prolongs nuclear burning.

preprint2014arXiv

The Large Observatory For x-ray Timing

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

preprint2013arXiv

New Primary Mechanisms for the Synthesis of Rare 9Be in Early Supernovae

We present two new primary mechanisms for the synthesis of the rare nucleus 9Be, both triggered by nu-induced production of 3H followed by 4He(3H,gamma)7Li in the He shells of core-collapse supernovae. For progenitors of ~8M_sun, 7Li(3H,n_0)9Be occurs during the rapid expansion of the shocked He shell. Alternatively, for ultra-metal-poor progenitors of ~11-15M_sun, 7Li(n,gamma)8Li(n,gamma)9Li(e^-anti-nu_e)9Be occurs with neutrons produced by 4He(anti-nu_e,e^+n)3H, assuming a hard effective anti-nu_e spectrum from oscillations (which also leads to heavy element production through rapid neutron capture) and a weak explosion (so the 9Be survives shock passage). We discuss the associated production of 7Li and 11B, noting patterns in LiBeB production that might distinguish the new mechanisms from others.

preprint2013arXiv

Production of carbon-rich presolar grains from massive stars

About a year after core collapse supernova, dust starts to condense in the ejecta. In meteorites, a fraction of C-rich presolar grains (e.g., silicon carbide (SiC) grains of Type-X and low density graphites) are identified as relics of these events, according to the anomalous isotopic abundances. Several features of these abundances remain unexplained and challenge the understanding of core-collapse supernovae explosions and nucleosynthesis. We show, for the first time, that most of the measured C-rich grain abundances can be accounted for in the C-rich material from explosive He burning in core-collapse supernovae with high shock velocities and consequent high temperatures. The inefficiency of the $^{12}$C($α$,$γ$)$^{16}$O reaction relative to the rest of the $α$-capture chain at $T > 3.5\times10^8 \mathrm{K}$ causes the deepest He-shell material to be carbon rich and silicon rich, and depleted in oxygen. The isotopic ratio predictions in part of this material, defined here as the C/Si zone, are in agreement with the grain data. The high-temperature explosive conditions that our models reach at the bottom of the He shell, can also be representative of the nucleosynthesis in hypernovae or in the high-temperature tail of a distribution of conditions in asymmetric supernovae. Finally, our predictions are consistent with the observation of large $^{44}$Ca/$^{40}$Ca observed in the grains. This is due to the production of $^{44}$Ti together with $^{40}$Ca in the C/Si zone, and/or to the strong depletion of $^{40}$Ca by neutron captures.

preprint2012arXiv

LOFT: the Large Observatory For X-ray Timing

The LOFT mission concept is one of four candidates selected by ESA for the M3 launch opportunity as Medium Size missions of the Cosmic Vision programme. The launch window is currently planned for between 2022 and 2024. LOFT is designed to exploit the diagnostics of rapid X-ray flux and spectral variability that directly probe the motion of matter down to distances very close to black holes and neutron stars, as well as the physical state of ultra-dense matter. These primary science goals will be addressed by a payload composed of a Large Area Detector (LAD) and a Wide Field Monitor (WFM). The LAD is a collimated (<1 degree field of view) experiment operating in the energy range 2-50 keV, with a 10 m^2 peak effective area and an energy resolution of 260 eV at 6 keV. The WFM will operate in the same energy range as the LAD, enabling simultaneous monitoring of a few-steradian wide field of view, with an angular resolution of <5 arcmin. The LAD and WFM experiments will allow us to investigate variability from submillisecond QPO's to year-long transient outbursts. In this paper we report the current status of the project.

preprint2012arXiv

New Two-Dimensional Models of Supernova Explosions by the Neutrino-Heating Mechanism: Evidence for Different Instability Regimes in Collapsing Stellar Cores

The neutrino-driven explosion mechanism for core-collapse supernovae in its modern flavor relies on the additional support of hydrodynamical instabilities in achieving shock revival. Two possible candidates, convection and the so-called standing accretion shock instability (SASI), have been proposed for this role. In this paper, we discuss new successful simulations of supernova explosions that shed light on the relative importance of these two instabilities. While convection has so far been observed to grow first in self-consistent hydrodynamical models with multi-group neutrino transport, we here present the first such simulation in which the SASI grows faster while the development of convection is initially inhibited. We illustrate the features of this SASI-dominated regime using an explosion model of a 27 solar mass progenitor, which is contrasted with a convectively-dominated model of an 8.1 solar mass progenitor with subsolar metallicity, whose early post-bounce behavior is more in line with previous 11.2 and 15 solar mass explosion models. We analyze the conditions discriminating between the two different regimes, showing that a high mass-accretion rate and a short advection time-scale are conducive for strong SASI activity. We also briefly discuss some important factors for capturing the SASI-driven regime, such as general relativity, the progenitor structure, a nuclear equation of state leading to a compact proto-neutron star, and the neutrino treatment. Finally, we evaluate possible implications of our findings for 2D and 3D supernova simulations. Our results show that a better understanding of the SASI and convection in the non-linear regime is required.

preprint2012arXiv

Superburst Models for Neutron Stars with Hydrogen and Helium-Rich Atmospheres

Superbursts are rare day-long Type I X-ray bursts due to carbon flashes on accreting neutron stars in low-mass X-ray binaries. They heat the neutron star envelope such that the burning of accreted hydrogen and helium becomes stable, and the common shorter X-ray bursts are quenched. Short bursts reappear only after the envelope cools down. We study multi-zone one-dimensional models of the neutron star envelope, in which we follow carbon burning during the superburst, and we include hydrogen and helium burning in the atmosphere above. We investigate both the case of a solar composition and a helium-rich atmosphere. This allows us to study for the first time a wide variety of thermonuclear burning behavior as well as the transitions between the different regimes in a self-consistent manner. For solar composition, burst quenching ends much sooner than previously expected. This is because of the complex interplay between the 3-alpha, hot CNO, and CNO breakout reactions. Stable burning of hydrogen and helium transitions via marginally stable burning (mHz QPOs) to less energetic bursts with short recurrence times. We find a short-lived bursting mode where weaker and stronger bursts alternate. Eventually the bursting behavior changes back to that of the pre-superburst bursts. Because of the scarcity of observations, this transition has not been directly detected after a superburst. Using the MINBAR burst catalog we identify the shortest upper limit on the quenching time for 4U 1636-536, and derive further constraints on the time scale on which bursts return.

preprint2011arXiv

Multi-Zone Models of Superbursts from Accreting Neutron Stars

Superbursts are rare and energetic thermonuclear carbon flashes observed to occur on accreting neutron stars. We create the first multi-zone models of series of superbursts using a stellar evolution code. We self-consistently build up the fuel layer at different rates, spanning the entire range of observed mass accretion rates for superbursters. For all models light curves are presented. They generally exhibit a shock breakout, a precursor burst due to shock heating, and a two-component power-law decay. Shock heating alone is sufficient for a bright precursor, that follows the shock breakout on a short dynamical time scale due to the fall-back of expanded layers. Models at the highest accretion rates, however, lack a shock breakout, precursor, and the first power law decay component. The ashes of the superburst that form the outer crust are predominantly composed of iron, but a superburst leaves a silicon-rich layer behind in which the next one ignites. Comparing the model light curves to an observed superburst from 4U 1636-53, we find for our accretion composition the best agreement with a model at three times the observed accretion rate. We study the dependence on crustal heating of observables such as the recurrence time and the decay time scale. It remains difficult, however, to constrain crustal heating, if there is no good match with the observed accretion rate, as we see for 4U 1636-53.

preprint2010arXiv

Multi-Instrument X-ray Observations of Thermonuclear Bursts with Short Recurrence Times

Type I X-ray bursts from low-mass X-ray binaries result from a thermonuclear runaway in the material accreted onto the neutron star. Although typical recurrence times are a few hours, consistent with theoretical ignition model predictions, there are also observations of bursts occurring as promptly as ten minutes or less after the previous event. We present a comprehensive assessment of this phenomenon using a catalog of 3387 bursts observed with the BeppoSAX/WFCs and RXTE/PCA X-ray instruments. This catalog contains 136 bursts with recurrence times of less than one hour, that come in multiples of up to four events, from 15 sources. Short recurrence times are not observed from so-called ultra-compact binaries, indicating that hydrogen burning processes play a crucial role. As far as the neutron star spin frequency is known, these sources all spin fast at over 500 Hz; the rotationally induced mixing may explain burst recurrence times of the order of 10 min. Short recurrence time bursts generally occur at all mass accretion rates where normal bursts are observed, but for individual sources the short recurrence times may be restricted to a smaller interval of accretion rate. The fraction of such bursts is roughly 30%. We also report the shortest known recurrence time of 3.8 minutes.

preprint2003arXiv

Neutrino Nucleosynthesis

We study neutrino process nucleosynthesis in massive stars using newly calculated cross sections, an expanded reaction network, and complete and self-consistent models of the progenitor star. We reevaluate the production of light isotopes from abundant progenitors as well as that of rare, heavy, proton-rich isotopes. In particular, new results are given for B11, F19, La138, and Ta180. The production of these isotopes places limits on neutrino spectrum and oscialltions.

preprint2003arXiv

Sensitivity of the C and O production on the 3-alpha rate

We investigate the dependence of the carbon and oxygen production in stars on the 3-alpha rate by varying the energy of the 02+-state of 12C and determine the resulting yields for a selection of low-mass, intermediate-mass, and massive stars. The yields are obtained using modern stellar evolution codes that follow the entire evolution of massive stars, including the supernova explosion, and consider in detail the 3rd dredge-up process during the thermally pulsating asymptotic giant branch of low-mass and intermediate-mass stars. Our results show that the C and O production in massive stars depends strongly on the initial mass, and that it is crucial to follow the entire evolution. A rather strong C production during the He-shell flashes compared to quiescent He burning leads to a lower sensitivity of the C and O production in low-mass and intermediate-mass stars on the 3-alpha-rate than predicted in our previous work. In particular, the C production of intermediate-mass stars seems to have a maximum close to the actual value of the 02+ energy level of 12C.

preprint2001arXiv

Formation of High Mass X-ray Black Hole Binaries

The discrepancy in the past years of many more black-hole soft X-ray transients (SXTs), of which a dozen have now been identified, had challenged accepted wisdom in black hole evolution. The sheer number of inferred SXTs requires that many, if not most, stars of ZAMS masses $20-35\msun$ end up as black holes (Portegies Zwart et al. 1997; Ergma & van den Heuvel 1998). In this paper we show that this can be understood by challenging the accepted wisdom that the result of helium core burning in a massive star is independent of whether the core is covered by a hydrogen envelope, or "naked" while it burns. The latter case occurs in binaries when the envelope of the more massive star is transferred to the companion by Roche Lobe overflow while in either main sequence or red giant stage. For solar metallicity, whereas the helium cores which burn while naked essentially never go into high-mass black holes, those that burn while clothed do so, beginning at ZAMS mass $\sim 20\msun$, the precise mass depending on the $^{12}C(α,γ)^{16}$O rate as we outline. In this way the SXTs can be evolved, provided that the H envelope of the massive star is removed only following the He core burning. We also show that in order to evolve a black hole of mass $\gsim 10\msun$ such as observed in Cyg X-1, even employing extremely massive progenitors of ZAMS mass $\gsim 60\msun$ for the black hole, the core must be covered by hydrogen during a substantial fraction of the core burning. In other words, the progenitor must be a WNL star. We evolve Cyg X-1 in an analogous way to which the SXTs are evolved, the difference being that the companion in Cyg X-1 is more massive than those in the SXTs, so that Cyg X-1 shines continuously.

preprint2001arXiv

Massive Star Evolution: Nucleosynthesis and Nuclear Reaction Rate Uncertainties

We present a nucleosynthesis calculation of a 25 solar mass star of solar composition that includes all relevant isotopes up to polonium. In particular, all stable isotopes and necessary nuclear reaction rates are covered. We follow the stellar evolution from hydrogen burning till iron core collapse and simulate the explosion using a ``piston'' approach. We discuss the influence of two key nuclear reaction rates, C12(a,g) and Ne22(a,n), on stellar evolution and nucleosynthesis. The former significantly influences the resulting core sizes (iron, silicon, oxygen) and the overall presupernova structure of the star. It thus has significant consequences for the supernova explosion itself and the compact remnant formed. The later rate considerably affects the s-process in massive stars and we demonstrate the changes that different currently suggested values for this rate cause.

preprint2000arXiv

Pair-Instability Supernovae, Gravity Waves, and Gamma-Ray Transients

Growing theoretical evidence suggests that the first generation of stars may have been quite massive (~100-300 solar masses). If they retain their high mass until death, such stars will, after about 3Myr, make pair-instability supernovae. We consider the complete evolution of two zero-metallicity stars of 250 and 300 solar masses. Explosive oxygen and silicon burning cause the 130 solar mass helium core to explode, but explosive burning is unable to drive an explosion in the 300 solar mass star and it collapses to a black hole. For this star, the calculated angular momentum in the presupernova model is sufficient to delay black hole formation and the star initially forms a 50 solar mass, 1000km core within which neutrinos are trapped. Although the star does not become dynamically unstable, the calculated growth time of secular rotational instabilities is shorter than the black hole formation time, and such instabilities may develop. We estimate the energy and amplitude of the gravitational waves emitted during this collapse. After the black hole forms, accretion continues through a disk. Although the disk is far too large and cool to transport energy efficiently to the rotational axis by neutrino annihilation, it has ample potential energy to produce a 1e54erg jet driven by magnetic fields. The interaction of this jet with surrounding circumstellar gas may produce an energetic gamma-ray transient, but given the redshift and time scale, this is probably not a model for typical gamma-ray bursts.

preprint1997arXiv

B[e] supergiants: What is their evolutionary status?

In this paper, we investigate the evolutionary status of B[e]~stars from the point of view of stellar evolution theory. We try to answer to the question of how massive hot supergiants --- i.e. evolved stars --- can be capable of producing a circumstellar disk. We find and discuss three possibilities: very massive evolved main sequence stars close to critical rotation due to their proximity to their Eddington-limit, blue supergiants which have just left the red supergiant branch, and single star merger remnants of a close binary system. While the latter process seems to be required to understand the properties of the spectroscopic binary R4 in the LMC, the other two scenarios may be capable of explaining the distribution of the B[e] stars in the HR~diagram. The three scenarios make different predictions about the duration of the B[e]~phase, the time integrated disk mass and the stellar properties during the B[e]~phase, which may ultimately allow to distinguish them observationally.