Source author record

R. Farmer

R. Farmer 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

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

14 published item(s)

preprint2022arXiv

Pair-instability mass loss for top-down compact object mass calculations

Population synthesis relies on semi-analytic formulae to determine masses of compact objects from the (helium or carbon-oxygen) cores of collapsing stars. Such formulae are combined across mass ranges that span different explosion mechanisms, potentially introducing artificial features in the compact object mass distribution. Such artifacts impair the interpretation of gravitational-wave observations. We propose a "top-down" remnant mass prescription where we remove mass from the star for each possible mass-loss mechanism, instead of relying on the fallback onto a "proto-compact-object" to get the final mass. For one of these mass-loss mechanisms, we fit the metallicity-dependent mass lost to pulsational-pair instability supernovae from numerical simulations. By imposing no mass loss in the absence of pulses, our approach recovers the existing compact object masses prescription at the low mass end and ensures continuity across the core-collapse/pulsational-pair-instability regime. Our remnant mass prescription can be extended to include other mass-loss mechanisms at the final collapse.

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.

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

Laminar Flame Speeds in Degenerate Oxygen-Neon Mixtures

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

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

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.

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.

preprint2016arXiv

On Variations Of Pre-Supernova Model Properties

We explore the variation in single star 15-30 $\rm{M}_{\odot}$, non-rotating, solar metallicity, pre-supernova MESA models due to changes in the number of isotopes in a fully-coupled nuclear reaction network and adjustments in the mass resolution. Within this two-dimensional plane we quantitatively detail the range of core masses at various stages of evolution, mass locations of the main nuclear burning shells, electron fraction profiles, mass fraction profiles, burning lifetimes, stellar lifetimes, and compactness parameter at core-collapse for models with and without mass loss. Up to carbon burning we generally find mass resolution has a larger impact on the variations than the number of isotopes, while the number of isotopes plays a more significant role in determining the span of the variations for neon, oxygen and silicon burning. Choice of mass resolution dominates the variations in the structure of the intermediate convection zone and secondary convection zone during core and shell hydrogen burning respectively, where we find a minimum mass resolution of $\approx$ 0.01 $\rm{M}_{\odot}$ is necessary to achieve convergence in the helium core mass at the $\approx$5% level. On the other hand, at the onset of core-collapse we find $\approx$30% variations in the central electron fraction and mass locations of the main nuclear burning shells, a minimum of $\approx$127 isotopes is needed to attain convergence of these values at the $\approx$10% level.

preprint2016arXiv

Properties of Carbon-Oxygen White Dwarfs From Monte Carlo Stellar Models

We investigate properties of carbon-oxygen white dwarfs with respect to the composite uncertainties in the reaction rates using the stellar evolution toolkit, Modules for Experiments in Stellar Astrophysics (MESA) and the probability density functions in the reaction rate library STARLIB. These are the first Monte Carlo stellar evolution studies that use complete stellar models. Focusing on 3 M$_{\odot}$ models evolved from the pre main-sequence to the first thermal pulse, we survey the remnant core mass, composition, and structure properties as a function of 26 STARLIB reaction rates covering hydrogen and helium burning using a Principal Component Analysis and Spearman Rank-Order Correlation. Relative to the arithmetic mean value, we find the width of the 95\% confidence interval to be $ΔM_{\rm 1TP}$ $\approx$ 0.019 M$_{\odot}$ for the core mass at the first thermal pulse, $Δ$$t_{\rm{1TP}}$ $\approx$ 12.50 Myr for the age, $Δ\log(T_{\rm c}/{\rm K}) \approx$ 0.013 for the central temperature, $Δ\log(ρ_{\rm c}/{\rm g \ cm}^{-3}) \approx$ 0.060 for the central density, $ΔY_{\rm{e,c}} \approx$ 2.6$\times$10$^{-5}$ for the central electron fraction, $ΔX_{\rm c}(^{22}\rm{Ne}) \approx$ 5.8$\times$10$^{-4}$, $ΔX_{\rm c}(^{12}\rm{C}) \approx$ 0.392, and $ΔX_{\rm c}(^{16}\rm{O}) \approx$ 0.392. Uncertainties in the experimental $^{12}$C($α,γ)^{16}\rm{O}$, triple-$α$, and $^{14}$N($p,γ)^{15}\rm{O}$ reaction rates dominate these variations. We also consider a grid of 1 to 6 M$_{\odot}$ models evolved from the pre main-sequence to the final white dwarf to probe the sensitivity of the initial-final mass relation to experimental uncertainties in the hydrogen and helium reaction rates.

preprint2015arXiv

On Carbon Burning in Super Asymptotic Giant Branch Stars

We explore the detailed and broad properties of carbon burning in Super Asymptotic Giant Branch (SAGB) stars with 2755 MESA stellar evolution models. The location of first carbon ignition, quenching location of the carbon burning flames and flashes, angular frequency of the carbon core, and carbon core mass are studied as a function of the ZAMS mass, initial rotation rate, and mixing parameters such as convective overshoot, semiconvection, thermohaline and angular momentum transport. In general terms, we find these properties of carbon burning in SAGB models are not a strong function of the initial rotation profile, but are a sensitive function of the overshoot parameter. We quasi-analytically derive an approximate ignition density, $ρ_{ign} \approx 2.1 \times 10^6$ g cm$^{-3}$, to predict the location of first carbon ignition in models that ignite carbon off-center. We also find that overshoot moves the ZAMS mass boundaries where off-center carbon ignition occurs at a nearly uniform rate of $ΔM_{\rm ZAMS}$/$Δf_{\rm{ov}}\approx$ 1.6 $M_{\odot}$. For zero overshoot, $f_{\rm{ov}}$=0.0, our models in the ZAMS mass range $\approx$ 8.9 to 11 $M_{\odot}$ show off-center carbon ignition. For canonical amounts of overshooting, $f_{\rm{ov}}$=0.016, the off-center carbon ignition range shifts to $\approx$ 7.2 to 8.8 $M_{\odot}$. Only systems with $f_{\rm{ov}}$ $\geq 0.01$ and ZAMS mass $\approx$ 7.2-8.0 $M_{\odot}$ show carbon burning is quenched a significant distance from the center. These results suggest a careful assessment of overshoot modeling approximations on claims that carbon burning quenches an appreciable distance from the center of the carbon core.

preprint2014arXiv

Prospects for detecting asteroseismic binaries in Kepler data

Asteroseismology may in principle be used to detect unresolved stellar binary systems comprised of solar-type stars and/or red giants. This novel method relies on the detection of the presence of two solar-like oscillation spectra in the frequency spectrum of a single lightcurve. Here, we make predictions of the numbers of systems that may be detectable in data already collected by the NASA Kepler Mission. Our predictions, which are based upon TRILEGAL and BiSEPS simulations of the Kepler field of view, indicate that as many as 200 or more ``asteroseismic binaries'' may be detectable in this manner. Most of these binaries should be comprised of two He-core-burning red giants. Owing largely to the limited numbers of targets with the requisite short-cadence Kepler data, we expect only a small number of detected binaries containing solar-type stars. The predicted yield of detections is sensitive to the assumed initial mass ratio distribution of the binary components and therefore represents a sensitive calibration of the much debated initial mass ratio distribution near mass ratio unity.

preprint2014arXiv

The PLATO 2.0 Mission

PLATO 2.0 has recently been selected for ESA's M3 launch opportunity (2022/24). Providing accurate key planet parameters (radius, mass, density and age) in statistical numbers, it addresses fundamental questions such as: How do planetary systems form and evolve? Are there other systems with planets like ours, including potentially habitable planets? The PLATO 2.0 instrument consists of 34 small aperture telescopes (32 with 25 sec readout cadence and 2 with 2.5 sec candence) providing a wide field-of-view (2232 deg2) and a large photometric magnitude range (4-16 mag). It focusses on bright (4-11 mag) stars in wide fields to detect and characterize planets down to Earth-size by photometric transits, whose masses can then be determined by ground-based radial-velocity follow-up measurements. Asteroseismology will be performed for these bright stars to obtain highly accurate stellar parameters, including masses and ages. The combination of bright targets and asteroseismology results in high accuracy for the bulk planet parameters: 2%, 4-10% and 10% for planet radii, masses and ages, respectively. The planned baseline observing strategy includes two long pointings (2-3 years) to detect and bulk characterize planets reaching into the habitable zone (HZ) of solar-like stars and an additional step-and-stare phase to cover in total about 50% of the sky. PLATO 2.0 will observe up to 1,000,000 stars and detect and characterize hundreds of small planets, and thousands of planets in the Neptune to gas giant regime out to the HZ. It will therefore provide the first large-scale catalogue of bulk characterized planets with accurate radii, masses, mean densities and ages. This catalogue will include terrestrial planets at intermediate orbital distances, where surface temperatures are moderate. Coverage of this parameter range with statistical numbers of bulk characterized planets is unique to PLATO 2.0.

preprint2013arXiv

The true stellar parameters of the Kepler target list

Using population synthesis tools we create a synthetic Kepler Input Catalogue (KIC) and subject it to the Kepler Stellar Classification Program (SCP) method for determining stellar parameters such as the effective temperature Teff and surface gravity g. We achieve a satisfactory match between the synthetic KIC and the real KIC in the log g vs log Teff diagram, while there is a significant difference between the actual physical stellar parameters and those derived by the SCP of the stars in the synthetic sample. We find a median difference ΔTeff=+500K and Δlog g =-0.2dex for main-sequence stars, and ΔTeff=+50K and Δlog g =-0.5dex for giants, although there is a large variation across parameter space. For a MS star the median difference in g would equate to a ~3% increase in stellar radius and a consequent ~3% overestimate of the radius for any transiting exoplanet. We find no significant difference between ΔTeff and Δlog g for single stars and the primary star in a binary system. We also re-created the Kepler target selection method and found that the binary fraction is unchanged by the target selection. Binaries are selected in similar proportions to single star systems; the fraction of MS dwarfs in the sample increases from about 75% to 80%, and the giant star fraction decreases from 25% to 20%.