Researcher profile

Alexander Heger

Alexander Heger contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

14 published item(s)

preprint2024arXiv

Spectacular nucleosynthesis from early massive stars

Stars formed with initial mass over 50 Msun are very rare today, but they are thought to be more common in the early universe. The fates of those early, metal-poor, massive stars are highly uncertain. Most are expected to directly collapse to black holes, while some may explode as a result of rotationally powered engines or the pair-creation instability. We present the chemical abundances of J0931+0038, a nearby low-mass star identified in early followup of SDSS-V Milky Way Mapper, which preserves the signature of unusual nucleosynthesis from a massive star in the early universe. J0931+0038 has relatively high metallicity ([Fe/H] = -1.76 +/- 0.13) but an extreme odd-even abundance pattern, with some of the lowest known abundance ratios of [N/Fe], [Na/Fe], [K/Fe], [Sc/Fe], and [Ba/Fe]. The implication is that a majority of its metals originated in a single extremely metal-poor nucleosynthetic source. An extensive search through nucleosynthesis predictions finds a clear preference for progenitors with initial mass > 50 Msun, making J0931+0038 one of the first observational constraints on nucleosynthesis in this mass range. However the full abundance pattern is not matched by any models in the literature. J0931+0038 thus presents a challenge for the next generation of nucleosynthesis models and motivates study of high-mass progenitor stars impacted by convection, rotation, jets, and/or binary companions. Though rare, more examples of unusual early nucleosynthesis in metal-poor stars should be found in upcoming large spectroscopic surveys.

preprint2022arXiv

Long-term evolution of post-explosion Helium-star Companions of Type Iax Supernovae

Supernovae of Type Iax (SNe Iax) are an accepted faint subclass of hydrogen-free supernovae. Their origin, the nature of the progenitor systems, however, is an open question. Recent studies suggest that the weak deflagration explosion of a near-Chandrasekhar-mass white dwarf in a binary system with a helium star donor could be the origin of SNe Iax. In this scenario, the helium star donor is expected to survive the explosion. We use the one-dimensional stellar evolution codes \textsc{MESA} and \textsc{Kepler} to follow the post-impact evolution of the surviving helium companion stars. The stellar models are based on our previous hydrodynamical simulations of ejecta-donor interaction, and we explore the observational characteristics of these surviving helium companions. We find that the luminosities of the surviving helium companions increase significantly after the impact: They could vary from $2\mathord,500\,\mathrm{L_{\odot}}$ to $16\mathord,000\,\mathrm{L_{\odot}}$ for a Kelvin-Helmholtz timescale of about $10^{4}\,\mathrm{yr}$. After the star reaches thermal equilibrium, it evolves as an O-type hot subdwarf (sdO) star and continues its evolution along the evolutionary track of a normal sdO star with the same mass. Our results will help to identify the surviving helium companions of SNe Iax in future observations and to place new constraints on their progenitor models.

preprint2022arXiv

The Impact of the New $^{65\!}$As(p,$γ$)$^{66\!}$Se Reaction Rate on the Two-Proton Sequential Capture of $^{64}\!$Ge, Weak GeAs Cycles, and Type-I X-Ray Bursts such as the Clocked Burster GS 1826$-$24

We re-assess $^{65}$As(p,$γ$)$^{66}$Se reaction rates based on a set of proton thresholds of $^{66}$Se, $S_\mathrm{p}$($^{66}$Se), estimated from the experimental mirror nuclear masses, theoretical mirror displacement energies, and full $pf$-model space shell-model calculation. The self-consistent relativistic Hartree-Bogoliubov theory is employed to obtain the mirror displacement energies with much reduced uncertainty, and thus reducing the proton-threshold uncertainty up to 161 keV compared to the AME2020 evaluation. Using the simulation instantiated by the one-dimensional multi-zone hydrodynamic code, KEPLER, that closely reproduces the observed GS 1826$-$24 clocked bursts, the present forward and reverse $^{65}$As(p,$γ$)$^{66}$Se reaction rates based on a selected $S_\mathrm{p}$($^{66}$Se) = 2.469$\pm$0.054 MeV, and the latest $^{22}$Mg($α$,p)$^{25}$Al, $^{56}$Ni(p,$γ$)$^{57}$Cu(p,$γ$)$^{58}$Zn, $^{55}$Ni(p,$γ$)$^{56}$Cu, and $^{64}$Ge(p,$γ$)$^{65}$As reaction rates, we find that though the GeAs cycles is weakly established in the rapid-proton capture process path, the $^{65}$As(p,$γ$)$^{66}$Se reaction still strongly characterizes the burst tail end due to the two-proton sequential capture on $^{64}$Ge, not found by Cyburt et al. (2016) sensitivity study. The $^{65}$As(p,$γ$)$^{66}$Se reaction influences the abundances of nuclei $A$ = 64, 68, 72, 76, and 80 up to a factor of 1.4. The new $S_\mathrm{p}$($^{66}$Se) and the inclusion of the updated $^{22}$Mg($α$,p)$^{25}$Al reaction rate increases the production of $^{12}$C up to a factor of $4.5$ that is not observable and could be the main fuel for superburst. The waiting point status of and two-proton sequential capture on $^{64}$Ge, weak-cycle feature of GeAs at region heavier than $^{64}$Ge, and impact of other possible $S_\mathrm{p}$($^{66}$Se) are also discussed.

preprint2022arXiv

The Regulated NiCu Cycles with the new $^{57}$Cu(p,$γ$)$^{58}$Zn reaction rate and the Influence on Type-I X-Ray Bursts: GS 1826$-$24 Clocked Burster

During the X-ray bursts of GS 1826$-$24, "clocked burster", the nuclear reaction flow that surges through the rapid-proton capture process path has to pass through the NiCu cycles before reaching the ZnGa cycles that moderate the further extent of hydrogen burning in the region above germanium and selenium isotopes. The $^{57}$Cu(p,$γ$)$^{58}$Zn reaction located in the NiCu cycles plays an important role in influencing the burst light curves as found by Cyburt et al. (2016). We deduce the $^{57}$Cu(p,$γ$)$^{58}$Zn reaction rate based on the experimentally determined important nuclear structure information, isobaric-multiplet-mass equation, and large-scale shell model calculations. Based on the isobaric-multiplet-mass equation, we propose a possible order of $1^+_1$ and $2^+_3$ dominant resonance states and constrain the resonance energy of the $1^+_2$ state. The latter reduces the contribution of the $1^+_2$ dominant resonance state. The new reaction rate is up to a factor of four lower than the Forstner et al. (2001) rate recommended by JINA REACLIB v2.2 at the temperature regime sensitive to clocked bursts of GS 1826$-$24. Using the simulation from the one-dimensional implicit hydrodynamic code, KEPLER, to model the thermonuclear X-ray bursts of GS 1826$-$24 clocked burster, we find that the new $^{57}$Cu(p,$γ$)$^{58}$Zn coupled with the latest $^{56}$Ni(p,$γ$)$^{57}$Cu and $^{55}$Ni(p,$γ$)$^{56}$Cu reaction rates redistributes the reaction flow in the NiCu cycles and strongly influences the burst ash composition, whereas the $^{59}$Cu(p,$α$)$^{56}$Ni and $^{59}$Cu(p,$γ$)$^{60}$Zn reactions suppress the influence of the $^{57}$Cu(p,$γ$)$^{58}$Zn reaction and diminish the impact of nuclear reaction flow that by-passes the important $^{56}$Ni waiting point induced by the $^{55}$Ni(p,$γ$)$^{56}$Cu reaction on burst light curve.

preprint2022arXiv

The Role of the Hadron-Quark Phase Transition in Core-Collapse Supernovae

The hadron-quark phase transition in quantum chromodyanmics has been suggested as an alternative explosion mechanism for core-collapse supernovae. We study the impact of three different hadron-quark equations of state (EoS) with first-order (DD2F\_SF, STOS-B145) and second-order (CMF) phase transitions on supernova dynamics by performing 97 simulations for solar- and zero-metallicity progenitors in the range of $14\texttt{-}100\,\text{M}_\odot$. We find explosions only for two low-compactness models ($14 \text{M}_\odot$ and $16\,\text{M}_\odot$) with the DD2F\_SF EoS, both with low explosion energies of $\mathord{\sim}10^{50}\,\mathrm{erg}$. These weak explosions are characterised by a neutrino signal with several mini-bursts in the explosion phase due to complex reverse shock dynamics, in addition to the typical second neutrino burst for phase-transition driven explosions. The nucleosynthesis shows significant overproduction of nuclei such as $^{90}\mathrm{Zr}$ for the $14\,\text{M}_\odot$ zero-metallicity model and $^{94}\mathrm{Zr}$ for the $16\,\text{M}_\odot$ solar-metallicity model, but the overproduction factors are not large enough to place constraints on the occurrence of such explosions. Several other low-compactness models using the DD2F\_SF EoS and two high-compactness models using the STOS EoS end up as failed explosions and emit a second neutrino burst. For the CMF EoS, the phase transition never leads to a second bounce and explosion. For all three EoS, inverted convection occurs deep in the core of the proto-compact star due to anomalous behaviour of thermodynamic derivatives in the mixed phase, which heats the core to entropies up to $4k_\text{B}/\text{baryon}$ and may have a distinctive gravitational wave signature, also for a second-order phase transition.

preprint2021arXiv

On the origin of nitrogen at low metallicity

Understanding the evolution of the N/O ratio in the interstellar medium (ISM) of galaxies is essential if we are to complete our picture of the chemical evolution of galaxies at high redshift, since most observational calibrations of O/H implicitly depend upon the intrinsic N/O ratio. The observed N/O ratio, however, shows large scatter at low O/H, and is strongly dependent on galactic environment. We show that several heretofore unexplained features of the N/O distribution at low O/H can be explained by the N seen in metal-poor galaxies being mostly primary nitrogen that is returned to the ISM via pre-supernova winds from rapidly rotating massive stars ($M \gtrsim 10$ M$_\odot$, $v/v_{\rm crit} \gtrsim 0.4$). This mechanism naturally produces the observed N/O plateau at low O/H. We show that the large scatter in N/O at low O/H also arises naturally from variations in star-formation efficiency. By contrast, models in which the N and O come primarily from supernovae provide a very poor fit to the observed abundance distribution. We propose that the peculiar abundance patterns we observe at low O/H are a signature that dwarf galaxies retain little of their SN ejecta, leaving them with abundance patterns typical of winds.

preprint2020arXiv

A Minimum Dilution Scenario for Supernovae and Consequences for Extremely Metal-Poor Stars

To date no metal-free stars have been identified by direct observations. The most common method of constraining their properties is searching the spectra of the most metal-poor stars for the chemical elements created in the first stars and their supernova. In this approach, modelled supernova yields are compared to the observed abundance patterns in extremely metal-poor stars. The method typically only uses the abundance ratios, i.e., the yields are diluted to the observed level. Following the usual assumption of spherical symmetry we compute a simple lower limit of the mass a supernova can mix with and find that it is consistent with all published simulations of early chemical enrichment in the interstellar medium. For three different cases, we demonstrate that this dilution limit can change the conclusions from the abundance fitting. There is a large discrepancy between the dilution found in simulations of SN explosions in minihaloes and the dilution assumed in many abundance fits. Limiting the dilution can significantly alter the likelihood of which supernovae are possible progenitors of observed CEMP-no stars. In particular, some of the faint, very low-yield SNe, which have been suggested as models for the abundance pattern of SMSS0313-6708, cannot explain the measured metal abundances, as their predicted metal yields are too small by two orders of magnitude. Altogether, the new dilution model presented here emphasizes the need to better understand the mixing and dilution behaviour of aspherical SNe.

preprint2020arXiv

Helium and Nitrogen Enrichment in Massive Main Sequence Stars: Mechanisms and Implications for the Origin of WNL Stars

The evolutionary paths taken by massive stars with $M \gtrsim 60 \, \mathrm{M}_\odot$ remain substantially uncertain. They begin their lives as main sequence (MS) O-stars. Depending on their masses, rotation rates, and metallicities, they can then encounter a wide range of evolutionary states with an equally broad set of possible surface compositions and spectral classifications. We present a new grid of calculations for the evolution of such stars that covers a broad range in mass, M/M$_\odot = 60$ to $150$, rotation rate, $v \, / \, v_{\rm crit} = 0$ to $0.6$, metallicity, $[\mathrm{Fe}/\mathrm{H}] = -4$ to $0$, and $α$-element enhancement, $[α/\mathrm{Fe}] = 0$ to $0.4$. We show that rotating stars undergo rotationally-induced dredge-up of nucleosynthetic products, mostly He and N, to their surfaces while still on the MS. Non-rotating metal-rich stars also reveal the products of nucleosynthesis on their surfaces because even modest amounts of mass loss expose their "fossil" convective cores: regions that are no longer convective, but which were part of the convective core at an early stage in the star's evolution. Thus surface enhancement of He and N is expected for rotating stars at all metallicities, and for non-rotating stars if they are relatively metal-rich. We calculate a stellar atmosphere for a representative model from our grid, properly accounting for He- and N-enhancement, and show that the resulting spectrum provides a good match to observed WNL stars, strongly suggesting that the physical mechanisms we have identified are the ultimate cause of the WNL phase.

preprint2020arXiv

Multi-epoch X-ray burst modelling: MCMC with large grids of 1D simulations

Type-I X-ray bursts are recurring thermonuclear explosions on the surface of accreting neutron stars. Matching observed bursts to computational models can help to constrain system properties, such as the neutron star mass and radius, crustal heating rates, and the accreted fuel composition, but systematic parameter studies to date have been limited. We apply Markov chain Monte Carlo methods to 1D burst models for the first time, and obtain system parameter estimations for the `Clocked Burster', GS 1826$-$238, by fitting multiple observed epochs simultaneously. We explore multiple parameters which are often held constant, including the neutron star mass, crustal heating rate, and hydrogen composition. To improve the computational efficiency, we precompute a grid of 3840 KEPLER models - the largest set of 1D burst simulations to date - and by interpolating over the model grid, we can rapidly sample burst predictions. We obtain estimates for a CNO metallicity of $Z_\mathrm{CNO} = 0.010^{+0.005}_{-0.004}$, a hydrogen fraction of $X_0 = 0.74^{+0.02}_{-0.03}$, a distance of $d \sqrt{ξ_\mathrm{b}} = 6.5^{+0.4}_{-0.6}\, \mathrm{kpc}$, and a system inclination of $i = {69^{+2}_{-3}}^{\circ}$.

preprint2020arXiv

On Monolithic Supermassive Stars

Supermassive stars have been proposed as the progenitors of the massive ($\sim 10^{9}\,\rm{M}_{\odot}$) quasars observed at $z\sim7$. Prospects for directly detecting supermassive stars with next-generation facilities depend critically on their intrinsic lifetimes, as well as their formation rates. We use the 1D stellar evolution code Kepler to explore the theoretical limiting case of zero-metallicity, non-rotating stars, formed monolithically with initial masses between $10\,\rm{kM}_{\odot}$ and $190\,\rm{kM}_{\odot}$. We find that stars born with masses between $\sim60\,\rm{kM}_{\odot}$ and $\sim150\,\rm{kM}_{\odot}$ collapse at the end of the main sequence, burning stably for $\sim1.5\,\rm{Myr}$. More massive stars collapse directly through the general relativistic instability after only a thermal timescale of $\sim3\,\rm{kyr}$--$4\,\rm{kyr}$. The expected difficulty in producing such massive, thermally-relaxed objects, together with recent results for currently preferred rapidly-accreting formation models, suggests that such ``truly direct'' or ``dark'' collapses may not be typical for supermassive objects in the early Universe. We close by discussing the evolution of supermassive stars in the broader context of massive primordial stellar evolution and the possibility of supermassive stellar explosions.

preprint2020arXiv

Production of Lithium in Primordial Supernovae

The first generation of stars is quite unique. The absence of metals likely affects their formation, with current models suggesting a much more top-heavy initial mass fraction than what we observe today, and some of their other properties, such as rotation rates and binarity, are largely unknown or constrained by direct observations. But even non-rotation single stars of a given mass will evolve quite differently due to the absence of the metals: the stars will mostly remain much more compact until their death, with the hydrogen-rich later reaching down ten teems deeper in radius then in modern stars. When they explode as supernovae, the exposure to the supernova neutrino flux is much enhanced, allowing for copious production of lithium. This production will not be constant for all stars but largely vary across the mass range. Such production even more challenges the presence of the Spite Plateau.

preprint2020arXiv

The Chemical Evolution of Iron-Peak Elements with Hypernovae

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

preprint2020arXiv

The impact of fallback on the compact remnants and chemical yields of core-collapse supernovae

Fallback in core-collapse supernovae plays a crucial role in determining the properties of the compact remnants and of the ejecta composition. We perform three-dimensional simulations of mixing and fallback for selected non-rotating supernova models to study how explosion energy and asymmetries correlate with the remnant mass, remnant kick, and remnant spin. We find that the strongest kick and spin is imparted by partial fallback in an asymmetric explosion. Black hole (BH) kicks of several hundred $\mathrm{km}\,\mathrm{s}^{-1}$ and spin parameters of $\mathord{\sim}0.25$ can be obtained in this scenario. If the initial explosion energy barely exceeds the envelope binding energy, stronger fallback results, and the remnant kick and spin remain small. If the explosion energy is high with respect to the envelope binding energy, there is little fallback with a small effect on the remnant kick, but the spin-up by fallback can be substantial. For a non-rotating $12\,\mathrm{M}_\odot$ progenitor, we find that the neutron star (NS) is spun up to millisecond periods. The high specific angular momentum of the fallback material can also lead to disk formation around black holes. Fallback may thus be a pathway towards millisecond-magnetar or collapsar-type engines for hypernovae and gamma-ray bursts that does not require rapid progenitor rotation. Within our small set of simulations, none reproduced the peculiar layered fallback necessary to explain the metal-rich iron-poor composition of many carbon-enhanced metal-poor (CEMP) stars. Models with different explosion energy and different realisations of asymmetries may, however, be compatible with CEMP abundance patterns.

preprint2019arXiv

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

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