sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Farmer, R.

Normalized to: Farmer, R.

22 article(s) in total. 215 co-authors, from 1 to 10 common article(s). Median position in authors list is 3,0.

[1]  [pdf] - 2077415
Polluting the pair-instability mass gap for binary black holes through super-Eddington accretion in isolated binaries
Comments: submitted
Submitted: 2020-04-10
The theory for single stellar evolution predicts a gap in the mass distribution of black holes (BHs) between approximately 45--130\Msun, the so-called ``the 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 90\Msun. We find no merging binary BH systems with a total mass exceeding 100\Msun. 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.
[2]  [pdf] - 2089142
On Stellar Evolution In A Neutrino Hertzsprung-Russell Diagram
Comments: 11 pages, 5 figures, accepted for publication in the astrophysical journal
Submitted: 2020-03-12
We explore the evolution of a select grid of solar metallicity stellar models from their pre-main sequence phase to near their final fates in a neutrino Hertzsprung-Russell diagram, where the neutrino luminosity replaces the traditional photon luminosity. Using a calibrated \MESA\ solar model for the solar neutrino luminosity ($L_{\nu,\odot}$ = 0.02398 $\cdot$ $L_{\gamma,\odot}$ = 9.1795 $\times$ 10$^{31}$ erg s$^{-1}$) as a normalization, we identify $\simeq$ 0.3 MeV electron neutrino emission from helium burning during the helium flash (peak $L_{\nu} / L_{\nu,\odot} \simeq$ 10$^4$, flux $\Phi_{\nu, {\rm He \ flash}} \simeq$ 170 (10 pc/$d$)$^{2}$ cm$^{-2}$ s$^{-1}$ for a star located at a distance of $d$ parsec, timescale $\simeq$ 3 days) and the thermal pulse (peak $L_{\nu} / L_{\nu,\odot} \simeq$ 10$^9$, flux $\Phi_{\nu, {\rm TP}} \simeq$ 1.7$\times$10$^7$ (10 pc/$d$)$^{2}$ cm$^{-2}$ s$^{-1}$, timescale $\simeq$ 0.1 yr) phases of evolution in low mass stars as potential probes for stellar neutrino astronomy. We also delineate the contribution of neutrinos from nuclear reactions and thermal processes to the total neutrino loss along the stellar tracks in a neutrino Hertzsprung-Russell diagram. We find, broadly but with exceptions, that neutrinos from nuclear reactions dominate whenever hydrogen and helium burn, and that neutrinos from thermal processes dominate otherwise.
[3]  [pdf] - 2089128
The expansion of stripped-envelope stars: consequences for supernovae and gravitational-wave progenitors
Comments: Main text 17 pages, 7 figures, accepted for publication in Astronomy & Astrophysics
Submitted: 2020-03-02
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.
[4]  [pdf] - 2061748
Sensitivity of the lower-edge of the pair instability black hole mass gap to the treatment of time dependent convection
Comments: 7 pages + 1 appendix, accepted in MNRAS, online results at
Submitted: 2020-02-19
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.
[5]  [pdf] - 2047677
Predictions for the hydrogen-free ejecta of pulsational pair-instability supernovae
Comments: 17 pages + 3 appendixes, submitted to A&A, simulation input and output will be available at acceptance
Submitted: 2020-02-12
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.
[6]  [pdf] - 2057890
Laminar Flame Speeds in Degenerate Oxygen-Neon Mixtures
Comments: 12 pages, 12 figures; Accepted to ApJ
Submitted: 2020-01-21
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.
[7]  [pdf] - 2026162
Mind the gap: The location of the lower edge of the pair instability supernovae black hole mass gap
Comments: 17 pages, 5 figures, 1 Table, Accepted ApJ
Submitted: 2019-10-28
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(\alpha,\gamma\right)^{16}O$ reaction rate within its $1\sigma$ 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.
[8]  [pdf] - 1979640
The diverse lives of progenitors of hydrogen-rich core-collapse supernovae: the role of binary interaction
Comments: The paper has been accepted for publication in Astronomy and Astrophysics
Submitted: 2019-07-15
Hydrogen-rich supernovae, known as Type II (SNe II), are the most common class of explosions observed following the collapse of the core of massive stars. We use analytical estimates and population synthesis simulations to assess the fraction of SNe II progenitors that are expected to have exchanged mass with a companion prior to explosion. We estimate that 1/3 to 1/2 of SN II progenitors have a history of mass exchange with a binary companion before exploding. The dominant binary channels leading to SN II progenitors involve the merger of binary stars. Mergers are expected to produce a diversity of SN II progenitor characteristics, depending on the evolutionary timing and properties of the merger. Alternatively, SN II progenitors from interacting binaries may have accreted mass from their companion, and subsequently been ejected from the binary system after their companion exploded. We show that the overall fraction of SN II progenitors that are predicted to have experienced binary interaction is robust against the main physical uncertainties in our models. However, the relative importance of different binary evolutionary channels is affected by changing physical assumptions. We further discuss ways in which binarity might contribute to the observed diversity of SNe II by considering potential observational signatures arising from each binary channel. For supernovae which have a substantial H-rich envelope at explosion (i.e., excluding Type IIb SNe), a surviving non-compact companion would typically indicate that the supernova progenitor star was in a wide, non-interacting binary. We argue that a significant fraction of even Type II-P SNe are expected to have gained mass from a companion prior to explosion.
[9]  [pdf] - 1916820
Modules for Experiments in Stellar Astrophysics (MESA): Pulsating Variable Stars, Rotation, Convective Boundaries, and Energy Conservation
Comments: 57 pages, 57 figures; Accepted to ApJS
Submitted: 2019-03-04, last modified: 2019-05-16
We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). RSP is a new functionality in MESAstar that models the non-linear radial stellar pulsations that characterize RR Lyrae, Cepheids, and other classes of variable stars. We significantly enhance numerical energy conservation capabilities, including during mass changes. For example, this enables calculations through the He flash that conserve energy to better than 0.001 %. To improve the modeling of rotating stars in MESA, we introduce a new approach to modifying the pressure and temperature equations of stellar structure, and a formulation of the projection effects of gravity darkening. A new scheme for tracking convective boundaries yields reliable values of the convective-core mass, and allows the natural emergence of adiabatic semiconvection regions during both core hydrogen- and helium-burning phases. We quantify the parallel performance of MESA on current generation multicore architectures and demonstrate improvements in the computational efficiency of radiative levitation. We report updates to the equation of state and nuclear reaction physics modules. We briefly discuss the current treatment of fallback in core-collapse supernova models and the thermodynamic evolution of supernova explosions. We close by discussing the new MESA Testhub software infrastructure to enhance source-code development.
[10]  [pdf] - 1867947
Massive runaways and walkaway stars
Comments: accepted version: added sec. 8.1 on other possible contributions to the observed velocities of runaway stars
Submitted: 2018-04-24, last modified: 2019-03-04
Anticipating the kinematic constraints from the Gaia mission, we perform an extensive numerical study of the evolution of massive binary systems to predict the peculiar velocities that stars obtain when their companion collapses and disrupts the system. Our aim is to (1) identify which predictions are robust against model uncertainties and assess their implications, (2) investigate which physical processes leave a clear imprint and may therefore be constrained observationally and (3) provide a suite of publicly available model predictions. We find that $22_{-8}^{+26}$% of all massive binary systems merge prior to the first core collapse in the system. Of the remainder, $86_{-9}^{+11}$% become unbound because of the core-collapse. Remarkably, this rarely produce runaway stars (i.e., stars with velocities above 30 km/s). These are outnumbered by more than an order of magnitude by slower unbound companions, or "walkaway stars". This is a robust outcome of our simulations and is due to the reversal of the mass ratio prior to the explosion and widening of the orbit, as we show analytically and numerically. We estimate a $10^{+5}_{-8}$% of massive stars to be walkaways and only $0.5^{+1.0}_{-0.4}$% to be runaways, nearly all of which have accreted mass from their companion. Our findings are consistent with earlier studies, however the low runaway fraction we find is in tension with observed fractions 10%. If Gaia confirms these high fractions of massive runaway stars resulting from binaries, it would imply that we are currently missing physics in the binary models. Finally, we show that high end of the mass distributions of runaway stars is very sensitive to the assumed black hole natal kicks and propose this as a potentially stringent test for the explosion mechanism. We discuss companions remaining bound which can evolve into X-ray and gravitational wave sources.
[11]  [pdf] - 1958055
Pulsational pair-instability supernovae in very close binaries
Comments: 20 pages, 15 figures, including 4 pages of appendices
Submitted: 2018-10-31
Pair-instability and pulsational pair-instability supernovae (PPISN) have not been unambiguously observed so far. They are, however, promising candidates for the progenitors of the heaviest binary black hole (BBH) mergers detected. Motivated by this, we discuss the implications of a PPISN happening with a close binary companion, and what impact these explosions have on the formation of merging BBHs through field binary evolution. For this, we have computed a set of models of metal-poor $(Z_\odot/10)$ helium stars using the MESA software instrument. For PPISN progenitors with pre-explosion masses $> 50M_\odot$ we find that, after a pulse, heat deposited throughout the layers of the star that remain bound cause it to expand to more than $100R_\odot$ for periods of $10^2-10^4$ yrs depending on the mass of the progenitor. This results in long-lived phases of Roche-lobe overflow or even common-envelope events if there is a close binary companion, leading to additional electromagnetic transients associated to PPISN eruptions. If we ignore the effect of these interactions, we find that mass loss from PPISNe reduces the final black hole spin by $\sim 30\%$, induces eccentricities that can be detected by the LISA observatory, and can produce a double-peaked distribution of measured chirp masses in BBH mergers observed by ground-based detectors.
[12]  [pdf] - 1637552
Modules for Experiments in Stellar Astrophysics (MESA): Convective Boundaries, Element Diffusion, and Massive Star Explosions
Comments: 64 pages, 61 figures; Accepted to AAS Journals
Submitted: 2017-10-23, last modified: 2018-01-03
We update the capabilities of the software instrument Modules for Experiments in Stellar Astrophysics (MESA) and enhance its ease of use and availability. Our new approach to locating convective boundaries is consistent with the physics of convection, and yields reliable values of the convective core mass during both hydrogen and helium burning phases. Stars with $M<8\,{\rm M_\odot}$ become white dwarfs and cool to the point where the electrons are degenerate and the ions are strongly coupled, a realm now available to study with MESA due to improved treatments of element diffusion, latent heat release, and blending of equations of state. Studies of the final fates of massive stars are extended in MESA by our addition of an approximate Riemann solver that captures shocks and conserves energy to high accuracy during dynamic epochs. We also introduce a 1D capability for modeling the effects of Rayleigh-Taylor instabilities that, in combination with the coupling to a public version of the STELLA radiation transfer instrument, creates new avenues for exploring Type II supernovae properties. These capabilities are exhibited with exploratory models of pair-instability supernova, pulsational pair-instability supernova, and the formation of stellar mass black holes. The applicability of MESA is now widened by the capability of importing multi-dimensional hydrodynamic models into MESA. We close by introducing software modules for handling floating point exceptions and stellar model optimization, and four new software tools -- MESAWeb, MESA-Docker, pyMESA, and -- to enhance MESA's education and research impact.
[13]  [pdf] - 1630310
The Impact of Nuclear Reaction Rate Uncertainties On The Evolution of Core-Collapse Supernova Progenitors
Comments: Accepted for publication in ApJS, 33 pages, 21 figures
Submitted: 2017-12-17
We explore properties of core-collapse supernova progenitors with respect to the composite uncertainties in the thermonuclear reaction rates by coupling the reaction rate probability density functions provided by the STARLIB reaction rate library with $\texttt{MESA}$ stellar models. We evolve 1000 15 $M_{\odot}$ models from the pre main-sequence to core O-depletion at solar and subsolar metallicities for a total of 2000 Monte Carlo stellar models. For each stellar model, we independently and simultaneously sample 665 thermonuclear reaction rates and use them in a $\texttt{MESA}$ in situ reaction network that follows 127 isotopes from $^{1}$H to $^{64}$Zn. With this framework we survey the core mass, burning lifetime, composition, and structural properties at five different evolutionary epochs. At each epoch we measure the probability distribution function of the variations of each property and calculate Spearman Rank-Order Correlation coefficients for each sampled reaction rate to identify which reaction rate has the largest impact on the variations on each property. We find that uncertainties in $^{14}$N$(p,\gamma)^{15}$O, triple-$\alpha$, $^{12}$C$(\alpha,\gamma)^{16}$O, $^{12}$C($^{12}$C,$p$)$^{23}$Na, $^{12}$C($^{16}$O,$p$)$^{27}$Al, $^{16}$O($^{16}$O,$n$)$^{31}$S, $^{16}$O($^{16}$O,$p$)$^{31}$P, and $^{16}$O($^{16}$O,$\alpha$)$^{28}$Si reaction rates dominate the variations of the properties surveyed. We find that variations induced by uncertainties in nuclear reaction rates grow with each passing phase of evolution, and at core H-, He-depletion are of comparable magnitude to the variations induced by choices of mass resolution and network resolution. However, at core C-, Ne-, and O-depletion, the reaction rate uncertainties can dominate the variation causing uncertainty in various properties of the stellar model in the evolution towards iron core-collapse.
[14]  [pdf] - 1604897
Neutrinos from beta processes in a presupernova: probing the isotopic evolution of a massive star
Comments: 14 pages, 5 figures
Submitted: 2017-09-06
We present a new calculation of the neutrino flux received at Earth from a massive star in the $\sim 24$ hours of evolution prior to its explosion as a supernova (presupernova). Using the stellar evolution code MESA, the neutrino emissivity in each flavor is calculated at many radial zones and time steps. In addition to thermal processes, neutrino production via beta processes is modeled in detail, using a network of 204 isotopes. We find that the total produced $\nu_{e}$ flux has a high energy spectrum tail, at $E \gtrsim 3 - 4$ MeV, which is mostly due to decay and electron capture on isotopes with $A = 50 - 60$. In a tentative window of observability of $E \gtrsim 0.5$ MeV and $t < 2$ hours pre-collapse, the contribution of beta processes to the $\nu_{e}$ flux is at the level of $\sim90\%$ . For a star at $D=1$ kpc distance, a 17 kt liquid scintillator detector would typically observe several tens of events from a presupernova, of which up to $\sim 30\%$ due to beta processes. These processes dominate the signal at a liquid argon detector, thus greatly enhancing its sensitivity to a presupernova.
[15]  [pdf] - 1560973
Presupernova neutrinos: realistic emissivities from stellar evolution
Comments: 17 pages, 6 figures
Submitted: 2015-11-09, last modified: 2017-04-10
We present a new calculation of neutrino emissivities and energy spectra from a massive star going through the advanced stages of nuclear burning (presupernova) in the months before becoming a supernova. The contributions from beta decay and electron capture, pair annihilation, plasmon decay, and the photoneutrino process are modeled in detail, using updated tabulated nuclear rates. We also use realistic conditions of temperature, density, electron fraction and nuclear isotopic composition of the star from the state of the art stellar evolution code MESA. Results are presented for a set of progenitor stars with mass between 15 $M_\odot$ and 30 $M_\odot$. It is found that beta processes contribute substantially to the neutrino emissivity above realistic detection thresholds of few MeV, at selected positions and times in the evolution of the star.
[16]  [pdf] - 1530264
Modules for Experiments in Stellar Astrophysics (MESA): Binaries, Pulsations, and Explosions
Comments: Published in the ApJ Supplement Series; updated to account for published erratum
Submitted: 2015-06-09, last modified: 2017-01-08
We substantially update the capabilities of the open-source software instrument Modules for Experiments in Stellar Astrophysics (MESA). MESA can now simultaneously evolve an interacting pair of differentially rotating stars undergoing transfer and loss of mass and angular momentum, greatly enhancing the prior ability to model binary evolution. New MESA capabilities in fully coupled calculation of nuclear networks with hundreds of isotopes now allow MESA to accurately simulate advanced burning stages needed to construct supernova progenitor models. Implicit hydrodynamics with shocks can now be treated with MESA, enabling modeling of the entire massive star lifecycle, from pre-main sequence evolution to the onset of core collapse and nucleosynthesis from the resulting explosion. Coupling of the GYRE non-adiabatic pulsation instrument with MESA allows for new explorations of the instability strips for massive stars while also accelerating the astrophysical use of asteroseismology data. We improve treatment of mass accretion, giving more accurate and robust near-surface profiles. A new MESA capability to calculate weak reaction rates "on-the-fly" from input nuclear data allows better simulation of accretion induced collapse of massive white dwarfs and the fate of some massive stars. We discuss the ongoing challenge of chemical diffusion in the strongly coupled plasma regime, and exhibit improvements in MESA that now allow for the simulation of radiative levitation of heavy elements in hot stars. We close by noting that the MESA software infrastructure provides bit-for-bit consistency for all results across all the supported platforms, a profound enabling capability for accelerating MESA's development.
[17]  [pdf] - 1532436
On Variations Of Pre-Supernova Model Properties
Comments: 32 pages, 1 table, 25 figures, Accepted ApJS
Submitted: 2016-11-03
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.
[18]  [pdf] - 1410146
Properties of Carbon-Oxygen White Dwarfs From Monte Carlo Stellar Models
Comments: Accepted for publication in The Astrophysical Journal; 19 Pages, 23 Figures, 5 Tables
Submitted: 2016-03-21
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 $\Delta M_{{\rm 1TP}}$ $\approx$ 0.019 M$_{\odot}$ for the core mass at the first thermal pulse, $\Delta$$t_{\rm{1TP}}$ $\approx$ 12.50 Myr for the age, $\Delta \log(T_{{\rm c}}/{\rm K}) \approx$ 0.013 for the central temperature, $\Delta \log(\rho_{{\rm c}}/{\rm g \ cm}^{-3}) \approx$ 0.060 for the central density, $\Delta Y_{\rm{e,c}} \approx$ 2.6$\times$10$^{-5}$ for the central electron fraction, $\Delta X_{\rm c}(^{22}\rm{Ne}) \approx$ 5.8$\times$10$^{-4}$, $\Delta X_{\rm c}(^{12}\rm{C}) \approx$ 0.392, and $\Delta X_{\rm c}(^{16}\rm{O}) \approx$ 0.392. Uncertainties in the experimental $^{12}$C($\alpha,\gamma)^{16}\rm{O}$, triple-$\alpha$, and $^{14}$N($p,\gamma)^{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.
[19]  [pdf] - 1245835
On Carbon Burning in Super Asymptotic Giant Branch Stars
Comments: Accepted ApJ; 23 pages, 21 figures, 5 tables
Submitted: 2015-06-02
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, $\rho_{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 $\Delta M_{\rm ZAMS}$/$\Delta 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.
[20]  [pdf] - 1179670
The PLATO 2.0 Mission
Rauer, H.; Catala, C.; Aerts, C.; Appourchaux, T.; Benz, W.; Brandeker, A.; Christensen-Dalsgaard, J.; Deleuil, M.; Gizon, L.; Goupil, M. -J.; Güdel, M.; Janot-Pacheco, E.; Mas-Hesse, M.; Pagano, I.; Piotto, G.; Pollacco, D.; Santos, N. C.; Smith, A.; -C., J.; Suárez; Szabó, R.; Udry, S.; Adibekyan, V.; Alibert, Y.; Almenara, J. -M.; Amaro-Seoane, P.; Eiff, M. Ammler-von; Asplund, M.; Antonello, E.; Ball, W.; Barnes, S.; Baudin, F.; Belkacem, K.; Bergemann, M.; Bihain, G.; Birch, A. C.; Bonfils, X.; Boisse, I.; Bonomo, A. S.; Borsa, F.; Brandão, I. M.; Brocato, E.; Brun, S.; Burleigh, M.; Burston, R.; Cabrera, J.; Cassisi, S.; Chaplin, W.; Charpinet, S.; Chiappini, C.; Church, R. P.; Csizmadia, Sz.; Cunha, M.; Damasso, M.; Davies, M. B.; Deeg, H. J.; DÍaz, R. F.; Dreizler, S.; Dreyer, C.; Eggenberger, P.; Ehrenreich, D.; Eigmüller, P.; Erikson, A.; Farmer, R.; Feltzing, S.; Fialho, F. de Oliveira; Figueira, P.; Forveille, T.; Fridlund, M.; García, R. A.; Giommi, P.; Giuffrida, G.; Godolt, M.; da Silva, J. Gomes; Granzer, T.; Grenfell, J. L.; Grotsch-Noels, A.; Günther, E.; Haswell, C. A.; Hatzes, A. P.; Hébrard, G.; Hekker, S.; Helled, R.; Heng, K.; Jenkins, J. M.; Johansen, A.; Khodachenko, M. L.; Kislyakova, K. G.; Kley, W.; Kolb, U.; Krivova, N.; Kupka, F.; Lammer, H.; Lanza, A. F.; Lebreton, Y.; Magrin, D.; Marcos-Arenal, P.; Marrese, P. M.; Marques, J. P.; Martins, J.; Mathis, S.; Mathur, S.; Messina, S.; Miglio, A.; Montalban, J.; Montalto, M.; Monteiro, M. J. P. F. G.; Moradi, H.; Moravveji, E.; Mordasini, C.; Morel, T.; Mortier, A.; Nascimbeni, V.; Nelson, R. P.; Nielsen, M. B.; Noack, L.; Norton, A. J.; Ofir, A.; Oshagh, M.; Ouazzani, R. -M.; Pápics, P.; Parro, V. C.; Petit, P.; Plez, B.; Poretti, E.; Quirrenbach, A.; Ragazzoni, R.; Raimondo, G.; Rainer, M.; Reese, D. R.; Redmer, R.; Reffert, S.; Rojas-Ayala, B.; Roxburgh, I. W.; Salmon, S.; Santerne, A.; Schneider, J.; Schou, J.; Schuh, S.; Schunker, H.; Silva-Valio, A.; Silvotti, R.; Skillen, I.; Snellen, I.; Sohl, F.; Sousa, S. G.; Sozzetti, A.; Stello, D.; Strassmeier, K. G.; Švanda, M.; Szabó, Gy. M.; Tkachenko, A.; Valencia, D.; van Grootel, V.; Vauclair, S. D.; Ventura, P.; Wagner, F. W.; Walton, N. A.; Weingrill, J.; Werner, S. C.; Wheatley, P. J.; Zwintz, K.
Comments: 63 pages, 17 figures, submitted to Experimental Astronomy (ExA)
Submitted: 2013-10-02, last modified: 2014-03-04
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.
[21]  [pdf] - 1203290
Prospects for detecting asteroseismic binaries in Kepler data
Comments: 21 pages, 5 figures, ApJ Letters, accepted
Submitted: 2014-02-11
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.
[22]  [pdf] - 676242
The true stellar parameters of the Kepler target list
Comments: Submitted to MNRAS, 14 pages with 9 figures; v2 Corrected typo; v3 Replaced with version accepted by MNRAS 15 pages with 9 figures
Submitted: 2013-02-12, last modified: 2013-05-07
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 \Delta Teff=+500K and \Delta log g =-0.2dex for main-sequence stars, and \Delta Teff=+50K and \Delta 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 \Delta Teff and \Delta 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%.