Gravitational-wave (GW) detections of binary black hole (BH) mergers have begun to sample the cosmic BH mass distribution. The evolution of single stellar cores predicts a gap in the BH mass distribution due to pair-instability supernovae (PISNe). Determining the upper and lower edges of the BH mass gap can be useful for interpreting GW detections of merging BHs. We use
We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (
- NSF-PAR ID:
- 10398798
- Author(s) / Creator(s):
- ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more »
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astrophysical Journal Supplement Series
- Volume:
- 265
- Issue:
- 1
- ISSN:
- 0067-0049
- Format(s):
- Medium: X Size: Article No. 15
- Size(s):
- ["Article No. 15"]
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract MESA to evolve single, nonrotating, massive helium cores with a metallicity ofZ = 10−5, until they either collapse to form a BH or explode as a PISN, without leaving a compact remnant. We calculate the boundaries of the lower BH mass gap for S-factors in the range S(300 keV) = (77,203) keV b, corresponding to the ±3σ uncertainty in our high-resolution tabulated12C(α ,γ )16O reaction rate probability distribution function. We extensively test temporal and spatial resolutions for resolving the theoretical peak of the BH mass spectrum across the BH mass gap. We explore the convergence with respect to convective mixing and nuclear burning, finding that significant time resolution is needed to achieve convergence. We also test adopting a minimum diffusion coefficient to help lower-resolution models reach convergence. We establish a new lower edge of the upper mass gap asM lower≃M ⊙from the ±3σ uncertainty in the12C(α ,γ )16O rate. We explore the effect of a larger 3α rate on the lower edge of the upper mass gap, findingM lower≃M ⊙. We compare our results with BHs reported in the Gravitational-Wave Transient Catalog. -
more » « less
Abstract Feedback from massive stars plays an important role in the formation of star clusters. Whether a very massive star is born early or late in the cluster formation timeline has profound implications for the star cluster formation and assembly processes. We carry out a controlled experiment to characterize the effects of early-forming massive stars on star cluster formation. We use the star formation software suite
Torch , combining self-gravitating magnetohydrodynamics, ray-tracing radiative transfer,N -body dynamics, and stellar feedback, to model four initially identical 104M ⊙giant molecular clouds with a Gaussian density profile peaking at 521.5 cm−3. Using theTorch software suite through theAMUSE framework, we modify three of the models, to ensure that the first star that forms is very massive (50, 70, and 100M ⊙). Early-forming massive stars disrupt the natal gas structure, resulting in fast evacuation of the gas from the star-forming region. The star formation rate is suppressed, reducing the total mass of the stars formed. Our fiducial control model, without an early massive star, has a larger star formation rate and total efficiency by up to a factor of 3, and a higher average star formation efficiency per freefall time by up to a factor of 7. Early-forming massive stars promote the buildup of spatially separate and gravitationally unbound subclusters, while the control model forms a single massive cluster. -
more » « less
ABSTRACT Planetary engulfment events can occur while host stars are on the main sequence. The addition of rocky planetary material during engulfment will lead to refractory abundance enhancements in the host star photosphere, but the level of enrichment and its duration will depend on mixing processes that occur within the stellar interior, such as convection, diffusion, and thermohaline mixing. We examine engulfment signatures by modelling the evolution of photospheric lithium abundances. Because lithium can be burned before or after the engulfment event, it produces unique signatures that vary with time and host star type. Using mesa stellar models, we quantify the strength and duration of these signatures following the engulfment of a 1, 10, or 100 M⊕ planetary companion with bulk Earth composition, for solar-metallicity host stars with masses ranging from 0.5 to 1.4 M⊙. We find that lithium is quickly depleted via burning in low-mass host stars ($\lesssim 0.7 \, {\rm M}_\odot$) on a time-scale of a few hundred Myrs, but significant lithium enrichment signatures can last for Gyrs in G-type stars ($\sim \! 0.9 \, {\rm M}_{\odot }$). For more massive stars (1.3−1.4 M⊙), engulfment can enhance internal mixing and diffusion processes, potentially decreasing the surface lithium abundance. Our predicted signatures from exoplanet engulfment are consistent with observed lithium-rich solar-type stars and abundance enhancements in chemically inhomogeneous binary stars.
-
more » « less
Abstract Stellar evolution and numerical hydrodynamics simulations depend critically on access to fast, accurate, thermodynamically consistent equations of state. We present Skye, a new equation of state for fully ionized matter. Skye includes the effects of positrons, relativity, electron degeneracy, Coulomb interactions, nonlinear mixing effects, and quantum corrections. Skye determines the point of Coulomb crystallization in a self-consistent manner, accounting for mixing and composition effects automatically. A defining feature of this equation of state is that it uses analytic free energy terms and provides thermodynamic quantities using automatic differentiation machinery. Because of this, Skye is easily extended to include new effects by simply writing new terms in the free energy. We also introduce a novel
thermodynamic extrapolation scheme for extending analytic fits to the free energy beyond the range of the fitting data while preserving desirable properties like positive entropy and sound speed. We demonstrate Skye in action in theMESA stellar evolution software instrument by computing white dwarf cooling curves. -
more » « less
Abstract We present the discovery of a new double-detonation progenitor system consisting of a hot subdwarf B (sdB) binary with a white dwarf companion with a
P orb= 76.34179(2) minutes orbital period. Spectroscopic observations are consistent with an sdB star during helium core burning residing on the extreme horizontal branch. Chimera light curves are dominated by ellipsoidal deformation of the sdB star and a weak eclipse of the companion white dwarf. Combining spectroscopic and light curve fits, we find a low-mass sdB star,M sdB= 0.383 ± 0.028M ⊙with a massive white dwarf companion,M WD= 0.725 ± 0.026M ⊙. From the eclipses we find a blackbody temperature for the white dwarf of 26,800 K resulting in a cooling age of ≈25 Myr whereas ourMESA model predicts an sdB age of ≈170 Myr. We conclude that the sdB formed first through stable mass transfer followed by a common envelope which led to the formation of the white dwarf companion ≈25 Myr ago. Using theMESA stellar evolutionary code we find that the sdB star will start mass transfer in ≈6 Myr and in ≈60 Myr the white dwarf will reach a total mass of 0.92M ⊙with a thick helium layer of 0.17M ⊙. This will lead to a detonation that will likely destroy the white dwarf in a peculiar thermonuclear supernova. PTF1 J2238+7430 is only the second confirmed candidate for a double-detonation thermonuclear supernova. Using both systems we estimate that at least ≈1% of white dwarf thermonuclear supernovae originate from sdB+WD binaries with thick helium layers, consistent with the small number of observed peculiar thermonuclear explosions.
