Researcher profile

S. Justham

S. Justham contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
4topics
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

11 published item(s)

preprint2022arXiv

Discovery of a highly magnetic He-sdO star from a double-degenerate binary merger

Helium-rich hot subdwarf stars of spectral type O (He-sdO) are considered prime candidates for stellar merger remnants. Such events should lead to the generation of strong magnetic fields. However, no magnetic He-sdO has yet been unambiguously discovered despite the high magnetic rate (20%) among white dwarf stars, the progeny of hot subdwarfs. Here we present the discovery of a strong magnetic field (B = 353 $\pm$ 10 kG) from Zeeman-split hydrogen, helium, and metal lines in the optical X-SHOOTER spectrum of an He-sdO and present the first spectroscopic analysis of any magnetic hot subdwarf. For this we used line-blanketed Tlusty non-local thermodynamic equilibrium models and assumed a simple homogeneous magnetic field. The derived atmospheric parameters $T_\mathrm{eff}$ = 44900 $\pm$ 1000 K and log g = 5.93 $\pm$ 0.15 are typical for He-sdO stars, while the star is less hydrogen-poor than most He-sdOs at log n(He)/n(H) = +0.28 $\pm$ 0.10. The star is a slow rotator ($v_\mathrm{rot}\sin i$ < 40 km s$^{-1}$). Its chemical composition is N-rich and C- and O-poor, and the Si and S abundances are close to solar. Combining the atmospheric parameters with Gaia parallax and photometry, the stellar radius and luminosity are found to be typical for He-sdOs and place the star on the helium main sequence in the Hertzsprung-Russell diagram. Its mass of $0.93^{+0.44}_{-0.30}$ $M_\odot$, although uncertain, appears to be remarkably high. The strong magnetic field along with the atmospheric parameters and metal abundances provide overwhelming evidence for the double-degenerate merger scenario.

preprint2022arXiv

Properties and applications of a predicted population of runaway He-sdO/B stars ejected from single degenerate He-donor SNe

This study builds on previous works, producing the most extensive prediction of the properties of such a hypothetical population to date, taking into account both Chandrasekhar and non-Chandrasekhar mass events. These results are then used to define criteria for membership of this population and characterise putative subpopulations This study contains 6x10^6 individual ejection trajectories out of the Galactic plane calculated with the stellar kinematics framework SHyRT, which are analysed with regard to their bulk observational properties. These are then put into context with the only previously identified population member US\,708 and applied to a number of other possible candidate objects. We find that two additional previously observed objects possess properties to warrant a designation as candidate objects. Characterisation of these object with respect to the predicted population finds all of them to be extreme in at least one astrometric observable. We find that current observations support a Galactic SN rate on the order of ~3x10^-7/yr to ~2x10^-6/yr, three orders of magnitude below the inferred Galactic SN Ia rate and two orders of magnitude below the formation rate of predicted He-donor progenitors. The number of currently observed population members suggests that the He-donor scenariois not a dominant contributor to the number of observed SNe Ia. However, we find that, even at the low event rate suggested, the majority of possibly detectable population members is still undetected. The extreme nature of current population members suggests that a still larger number of objects has simply evaded detection up to this point, hinting at a higher contribution than is currently supported by observation. - abridged -

preprint2022arXiv

Rejuvenated accretors have less bound envelopes: Impact of Roche lobe overflow on subsequent common envelope events

Common-envelope (CE) evolution is an outstanding open problem in stellar evolution, critical to the formation of compact binaries including gravitational-wave sources. In the ``classical&#39;&#39; isolated binary evolution scenario for double compact objects, the CE is usually the second mass transfer phase. Thus, the donor star of the CE is the product of a previous binary interaction, often stable Roche-lobe overflow (RLOF). Because of the accretion of mass during the first RLOF, the main-sequence core of the accretor star grows and is ``rejuvenated&#39;&#39;. This modifies the core-envelope boundary region and decreases significantly the envelope binding energy for the remaining evolution. Comparing accretor stars from self-consistent binary models to stars evolved as single, we demonstrate that the rejuvenation can lower the energy required to eject a CE by $\sim 42-96\%$ for both black hole and neutron star progenitors, depending on the evolutionary stage and final orbital separation. Therefore, binaries experiencing first stable mass transfer may more easily survive subsequent CE events and result in possibly wider final separations compared to current predictions. Despite their high mass, our accretors also experience extended ``blue loops&#39;&#39;, which may have observational consequences for low-metallicity stellar populations and asteroseismology.

preprint2022arXiv

The redshift evolution of the binary black hole merger rate: a weighty matter

Gravitational wave detectors are starting to reveal the redshift evolution of the binary black hole (BBH) merger rate, $R_{\mathrm{BBH}}(z)$. We make predictions for $R_{\mathrm{BBH}}(z)$ as a function of black hole mass for systems originating from isolated binaries. To this end, we investigate correlations between the delay time and black hole mass by means of the suite of binary population synthesis simulations, COMPAS. We distinguish two channels: the common envelope (CE), and the stable Roche-lobe overflow (RLOF) channel, characterised by whether the system has experienced a common envelope or not. We find that the CE channel preferentially produces BHs with masses below about $30\rm{M}_{\odot}$ and short delay times ($t_{\rm delay} \lesssim 1$Gyr), while the stable RLOF channel primarily forms systems with BH masses above $30\rm{M}_{\odot}$ and long delay times ($t_{\rm delay} \gtrsim 1$Gyr). We provide a new fit for the metallicity specific star-formation rate density based on the Illustris TNG simulations, and use this to convert the delay time distributions into a prediction of $R_{\mathrm{BBH}}(z)$. This leads to a distinct redshift evolution of $R_{\mathrm{BBH}}(z)$ for high and low primary BH masses. We furthermore find that, at high redshift, $R_{\mathrm{BBH}}(z)$ is dominated by the CE channel, while at low redshift it contains a large contribution ($\sim 40\%$) from the stable RLOF channel. Our results predict that, for increasing redshifts, BBHs with component masses above $30\rm{M}_{\odot}$ will become increasingly scarce relative to less massive BBH systems. Evidence of this distinct evolution of $R_{\mathrm{BBH}}(z)$ for different BH masses can be tested with future detectors.

preprint2021arXiv

The cosmic carbon footprint of massive stars stripped in binary systems

The cosmic origin of carbon, a fundamental building block of life, is still uncertain. Yield predictions for massive stars are almost exclusively based on single star models, even though a large fraction interact with a binary companion. Using the MESA stellar evolution code, we predict the carbon ejected in the winds and supernovae of single and binary-stripped stars at solar metallicity. We find that binary-stripped stars are twice as efficient at producing carbon (1.5-2.6 times, depending on choices on the slope of the initial mass function and black hole formation). We confirm that this is because the convective helium core recedes in stars that have lost their hydrogen envelope, as noted previously. The shrinking of the core disconnects the outermost carbon-rich layers created during the early phase of helium burning from the more central burning regions. The same effect prevents carbon destruction, even when the supernova shock wave passes. The yields are sensitive to the treatment of mixing at convective boundaries, specifically during carbon-shell burning (variations up to 40%) and improving upon this should be a central priority for more reliable yield predictions. The yields are robust (variations less than 0.5%) across our range of explosion assumptions. Black hole formation assumptions are also important, implying that the stellar graveyard now explored by gravitational-wave detections may yield clues to better understand the cosmic carbon production. Our findings also highlight the importance of accounting for binary-stripped stars in chemical yield predictions and motivates further studies of other products of binary interactions.

preprint2020arXiv

Polluting the pair-instability mass gap for binary black holes through super-Eddington accretion in isolated binaries

The theory for single stellar evolution predicts a gap in the mass distribution of black holes (BHs) between approximately 45-130M$_{\odot}$, the so-called &#34;pair-instability mass gap&#34;. We examine whether BHs can pollute the gap after accreting from a stellar companion. To this end, we simulate the evolution of isolated binaries using a population synthesis code, where we allow for super-Eddington accretion. Under our most extreme assumptions, we find that at most about 2% of all merging binary BH systems contains a BH with a mass in the pair-instability mass gap, and we find that less than 0.5% of the merging systems has a total mass larger than 90M$_{\odot}$. We find no merging binary BH systems with a total mass exceeding 100M$_{\odot}$. We compare our results to predictions from several dynamical pathways to pair-instability mass gap events and discuss the distinguishable features. We conclude that the classical isolated binary formation scenario will not significantly contribute to the pollution of the pair-instability mass gap. The robustness of the predicted mass gap for the isolated binary channel is promising for the prospective of placing constraints on (i) the relative contribution of different formation channels, (ii) the physics of the progenitors including nuclear reaction rates, and (iii), tentatively, the Hubble parameter.

preprint2020arXiv

Predictions for the hydrogen-free ejecta of pulsational pair-instability supernovae

Present time-domain astronomy efforts will unveil a variety of rare transients. We focus here on pulsational pair-instability evolution, which can result in signatures observable with electromagnetic and gravitational waves. We simulate grids of bare helium stars to characterize the resulting black hole (BH) masses and ejecta composition, velocity, and thermal state. The stars do not react &#34;elastically&#34; to the thermonuclear explosion: there is not a one-to-one correspondence between pair-instability driven ignition and mass ejections, causing ambiguity in what is an observable pulse. In agreement with previous studies, we find that for carbon-oxygen core masses 28Msun< M_CO<30.5Msun the explosions are not strong enough to affect the surface. With increasing mass, they first cause large radial expansion (30.5Msun<M_CO<31.4Msun), and finally, also mass ejection episodes (M_CO>31.4Msun). The lowest mass to be fully disrupted in a pair-instability supernova is M_CO=57Msun. Models with M_CO>121Msun reach the photodisintegration regime, resulting in BHs with M_BH>125Msun. If the pulsating models produce BHs via (weak) explosions, the previously-ejected material might be hit by the blast wave. We characterize the H-free circumstellar material from the pulsational pair-instability of helium cores assuming simply that the ejecta maintain a constant velocity after ejection. Our models produce He-rich ejecta with mass 10^{-3}Msun<M_CSM<40Msun. These ejecta are typically launched at a few thousand \kms and reach distances of ~10^{12}-10^{15} cm before core-collapse. The delays between mass ejection events and the final collapse span a wide and mass-dependent range (from sub-hour to 10^4 years), and the shells ejected can also collide with each other. The range of properties we find suggests a possible connection with (some) type Ibn supernovae.

preprint2020arXiv

Sensitivity of the lower-edge of the pair instability black hole mass gap to the treatment of time dependent convection

Gravitational-wave detections are now probing the black hole (BH) mass distribution, including the predicted pair-instability mass gap. These data require robust quantitative predictions, which are challenging to obtain. The most massive BH progenitors experience episodic mass ejections on timescales shorter than the convective turn-over timescale. This invalidates the steady-state assumption on which the classic mixing-length theory relies. We compare the final BH masses computed with two different versions of the stellar evolutionary code \texttt{MESA}: (i) using the default implementation of \cite{paxton:18} and (ii) solving an additional equation accounting for the timescale for convective deceleration. In the second grid, where stronger convection develops during the pulses and carries part of the energy, we find weaker pulses. This leads to lower amounts of mass being ejected and thus higher final BH masses of up to $\sim$\,$5\,M_\odot$. The differences are much smaller for the progenitors which determine the maximum mass of BHs below the gap. This prediction is robust at $M_{\rm BH, max}\simeq 48\,M_\odot$, at least within the idealized context of this study. This is an encouraging indication that current models are robust enough for comparison with the present-day gravitational-wave detections. However, the large differences between individual models emphasize the importance of improving the treatment of convection in stellar models, especially in the light of the data anticipated from the third generation of gravitational wave detectors.

preprint2020arXiv

The expansion of stripped-envelope stars: consequences for supernovae and gravitational-wave progenitors

Massive binaries that merge as compact objects are the progenitors of gravitational-wave sources. Most of these binaries experience one or more phases of mass transfer, during which one of the stars loses part or all of its outer envelope and becomes a stripped-envelope star. The evolution of the size of these stripped stars is crucial in determining whether they experience further interactions and their final fate. We present new calculations of stripped-envelope stars based on binary evolution models computed with MESA. We use these to investigate their radius evolution as a function of mass and metallicity. We further discuss their pre-supernova observable characteristics and potential consequences of their evolution on the properties of supernovae from stripped stars. At high metallicity we find that practically all of the hydrogen-rich envelope is removed, in agreement with earlier findings. Only progenitors with initial masses below 10\Msun expand to large radii (up to 100\Rsun), while more massive progenitors stay compact. At low metallicity, a substantial amount of hydrogen remains and the progenitors can, in principle, expand to giant sizes (> 400\Rsun), for all masses we consider. This implies that they can fill their Roche lobe anew. We show that the prescriptions commonly used in population synthesis models underestimate the stellar radii by up to two orders of magnitude. We expect that this has consequences for the predictions for gravitational-wave sources from double neutron star mergers, in particular for their metallicity dependence.

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.

preprint2019arXiv

Mind the gap: The location of the lower edge of the pair instability supernovae black hole mass gap

Gravitational-wave detections are now starting to probe the mass distribution of stellar-mass black holes (BHs). Robust predictions from stellar models are needed to interpret these. Theory predicts the existence of a gap in the BH mass distribution because of pair-instability supernova. The maximum BH mass below the gap is the result of pulsational mass loss. We evolve massive helium stars through their late hydrodynamical phases of evolution using the open-source MESA stellar evolution code. We find that the location of the lower edge of the mass gap at 45$M_\odot$ is remarkably robust against variations in the metallicity ($\approx 3M_\odot$), the treatment of internal mixing ($\approx 1M_\odot$), stellar wind mass loss ($\approx 4M_\odot$), making it the most robust predictions for the final stages of massive star evolution. The reason is that the onset of the instability is dictated by the near-final core mass, which in turn sets the resulting BH mass. However, varying $^{12}C\left(α,γ\right)^{16}O$ reaction rate within its $1σ$ uncertainties shifts the location of the gap between $40M_\odot$ and $56M_\odot$. We provide updated analytic fits for population synthesis simulations. Our results imply that the detection of merging BHs can provide constraints on nuclear astrophysics. Furthermore, the robustness against metallicity suggests that there is a universal maximum for the location of the lower edge of the gap, which is insensitive to the formation environment and redshift for first-generation BHs. This is promising for the possibility to use the location of the gap as a &#34;standard siren&#34; across the Universe.