Source author record

S. Justham

S. Justham 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

16works
5topics
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

16 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'' 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''. 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'', 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 "pair-instability mass gap". 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 "elastically" 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 "standard siren" across the Universe.

preprint2014arXiv

On the role of recombination in common-envelope ejections

The energy budget in common-envelope events (CEEs) is not well understood, with substantial uncertainty even over to what extent the recombination energy stored in ionised hydrogen and helium might be used to help envelope ejection. We investigate the reaction of a red-giant envelope to heating which mimics limiting cases of energy input provided by the orbital decay of a binary during a CEE, specifically during the post-plunge-in phase during which the spiral-in has been argued to occur on a time-scale longer than dynamical. We show that the outcome of such a CEE depends less on the total amount of energy by which the envelope is heated than on how rapidly the energy was transferred to the envelope and on where the envelope was heated. The envelope always becomes dynamically unstable before receiving net heat energy equal to the envelope's initial binding energy. We find two types of outcome, both of which likely lead to at least partial envelope ejection: "runaway" solutions in which the expansion of the radius becomes undeniably dynamical, and superficially "self-regulated" solutions, in which the expansion of the stellar radius stops but a significant fraction of the envelope becomes formally dynamically unstable. Almost the entire reservoir of initial helium recombination energy is used for envelope expansion. Hydrogen recombination is less energetically useful, but is nonetheless important for the development of the dynamical instabilities. However, this result requires the companion to have already plunged deep into the envelope; therefore this release of recombination energy does not help to explain wide post-common-envelope orbits.

preprint2014arXiv

SDSSJ001153.08-064739.2, a cataclysmic variable with an evolved donor in the period gap

Secondary stars in cataclysmic variables (CVs) follow a well defined period-density relation. Thus, canonical donor stars in CVs are generally low-mass stars of spectral type M. However, several CVs have been observed containing secondary stars which are too hot for their inferred masses. This particular configuration can be explained if the donor stars in these systems underwent significant nuclear evolution before they reached contact. In this paper we present SDSSJ001153.08-064739.2 as an additional example belonging to this peculiar type of CV and discuss in detail its evolutionary history. We perform spectroscopic and photometric observations and make use of available Catalina Real-Time Transient Survey photometry to measure the orbital period of SDSSJ001153.08-064739.2 as 2.4 hours and estimate the white dwarf (Mwd>0.65Msun) and donor star (0.21Msun<Mdon<0.45Msun) masses, the mass ratio (q = 0.32 +- 0.08), the orbital inclination (47 degrees < i < 70 degrees), derive an accurate orbital ephemeris (T0 = 2453383.578 + E x 0.10028081), and report the detection of an outburst. We show that SDSSJ001153.08-064739.2 is one of the most extreme cases in which the donor star is clearly too hot for its mass. SDSSJ001153.08-064739.2 is therefore not only a peculiar CV containing an evolved donor star but also an accreting CV within the period gap. Intriguingly, approximately half of the total currently-observed sample of these peculiar CVs are located in the period gap with nearly the same orbital period.

preprint2013arXiv

Identification of the Long-Sought Common-Envelope Events

Common-envelope events (CEEs), during which two stars temporarily orbit within a shared envelope, are believed to be vital for the formation of a wide range of close binaries. For decades, the only evidence that CEEs actually occur has been indirect, based on the existence of systems that could not be otherwise explained. Here we propose a direct observational signature of CEE arising from a physical model where emission from matter ejected in a CEE is controlled by a recombination front as the matter cools. The natural range of timescales and energies from this model, as well the expected colors, light-curve shapes, ejection velocities and event rate, match those of a recently-recognized class of red transient outbursts.

preprint2012arXiv

Common Envelope Evolution: Where we stand and how we can move forward

This work aims to present our current best physical understanding of common-envelope evolution (CEE). We highlight areas of consensus and disagreement, and stress ideas which should point the way forward for progress in this important but long-standing and largely unconquered problem. Unusually for CEE-related work, we mostly try to avoid relying on results from population synthesis or observations, in order to avoid potentially being misled by previous misunderstandings. As far as possible we debate all the relevant issues starting from physics alone, all the way from the evolution of the binary system immediately before CEE begins to the processes which might occur just after the ejection of the envelope. In particular, we include extensive discussion about the energy sources and sinks operating in CEE, and hence examine the foundations of the standard energy formalism. Special attention is also given to comparing the results of hydrodynamic simulations from different groups and to discussing the potential effect of initial conditions on the differences in the outcomes. We compare current numerical techniques for the problem of CEE and also whether more appropriate tools could and should be produced (including new formulations of computational hydrodynamics, and attempts to include 3D processes within 1D codes). Finally we explore new ways to link CEE with observations. We compare previous simulations of CEE to the recent outburst from V1309 Sco, and discuss to what extent post-common-envelope binaries and nebulae can provide information, e.g. from binary eccentricities, which is not currently being fully exploited.

preprint2007arXiv

Detection of circumstellar material in a normal Type Ia Supernova

Type Ia supernovae are thought to be thermonuclear explosions of accreting white dwarfs that reach a critical mass limit. Despite their importance as cosmological distance indicators, the nature of their progenitors has remained controversial. Here we report the detection of circumstellar material in a normal Type Ia supernova. The expansion velocities, densities and dimensions of the circumstellar envelope indicate that this material was ejected from the progenitor system. The relatively low expansion velocities appear to favor a progenitor system where a white dwarf accretes material from a companion star which is in the red-giant phase at the time of explosion.