Stellar Evolution Code Research Articles

Long-term double synchronization in close-in gas giant planets

ABSTRACT Hot Jupiters, orbiting their host stars at extremely close distances, undergo tidal evolution, with some being engulfed by their stars due to angular momentum exchanges induced by tidal forces. However, achieving double synchronization can prolong their survival. Using the mesa stellar evolution code, combined with the magnetic braking model of Matt et al. (2015), we calculate 25 000 models with different metallicity and study how to attain the conditions that trigger the long-term double synchronization. Our results indicate that massive planets orbiting stars with lower convective turnover time are easier to achieve long-term double synchronization. The rotation angular velocity at the equilibrium point (Ωsta) is almost equal to orbital angular velocity of planet (n) for the majority of the main sequence lifetime if a system has undergone a long-term double synchronization, regardless of their state at this moment. We further compared our results with known parameters of giant planetary systems and found that those systems with larger planetary masses and lower convective turnover time seem to be less sensitive to changes in the tidal quality factor $Q^{\prime }_{_*}$. We suggest that for systems that fall on the state of Ωsta ≈ n, regardless of their current state, the synchronization will persist for a long time if orbital synchronization occurs at any stage of their evolution. Our results can be applied to estimate whether a system has experienced long-term double synchronization in the past or may experience it in the future.

The Influence of Stellar Rotation in Binary Systems on Core-collapse Supernova Progenitors and Multimessenger Signals

The detailed structure of core-collapse supernova progenitors is crucial for studying supernova explosion engines and the corresponding multimessenger signals. In this paper, we investigate the influence of stellar rotation on binary systems consisting of a 30M ⊙ donor star and a 20M ⊙ accretor using the MESA stellar evolution code. We find that through mass transfer in binary systems, fast-rotating red- and blue-supergiant progenitors can be formed within a certain range of the initial orbital periods, although the correlation is not linear. We also find that even with the same initial mass ratio of the binary system, the resulting final masses of the collapsars, the iron core masses, the compactness parameters, and the final rotational rates can vary widely and are sensitive to the initial orbital periods. For instance, the progenitors with strong convection form a thinner Si shell and a wider O shell compared to those in single-star systems. In addition, we conduct 2D self-consistent core-collapse supernova simulations with neutrino transport for these rotating progenitors derived from binary stellar evolution. We find that the neutrino and gravitational-wave signatures of these binary progenitors could exhibit significant variations. Progenitors with larger compactness parameters produce more massive proto-neutron stars, have higher mass accretion rates, and emit brighter neutrino luminosity and louder gravitational emissions. Finally, we observe stellar-mass black hole formation in some of our failed exploding models.

Effect of nuclear reactions rates and core boundary mixing on the seismology of red clump stars

Context. Modeling of the stars in the red clump (RC), that is, core helium-burning stars that have gone through a He flash, is challenging because of the uncertainties associated with the physical processes in their core and during the helium flash. By probing the internal stellar structure, asteroseismology allows us to constrain the core properties of RC stars and eventually, to improve our understanding of this evolutionary phase. Aims. We aim to quantify the impact on the seismic properties of the RC stars of the two main core modeling uncertainties: core boundary mixing, and helium-burning nuclear reaction rates. Methods. Using the MESA stellar evolution code, we computed models with different core boundary mixing as well as different 3α and 12C(α, γ)16O nuclear reaction rates. We investigated the impact of these parameters on the period spacing ΔΠ, which is a probe of the region around the core. Results. We find that different core boundary mixing schemes yield significantly different period spacings, with differences of 30 s between the maximum ΔΠ value computed with semiconvection and maximal overshoot. We show that an increased rate of 12C(α, γ)16O lengthens the core helium-burning phase, which extends the range of ΔΠ covered by the models during their evolution. This results in a difference of 10 s between the models computed with a nominal rate and a rate multiplied by 2, which exceeds the observational uncertainties. The effect of changing the 3α reaction rate is comparatively small. Conclusions. The core boundary mixing is the main source of uncertainty in the seismic modeling of RC stars. Moreover, the effect of the 12C(α, γ)16O is non-negligible, even though it is difficult to distinguish from the effect of the mixing. This degeneracy could be seen more frequently in the future in the new seismic data from the PLATO mission and through theoretical constraints from numerical simulations.

The Measurement of Masses of OB-type Stars from LAMOST DR5

The measurements of masses and luminosities of massive stars play an important role in understanding the formation and evolution of their host galaxies. In this work, we present the measurement of masses and luminosities of 2946 OB-type stars, including 78 O-type stars and 2868 B-type stars, based on their stellar parameters (effective temperature, surface gravity, and metallicity) and the Padova and Trieste Stellar Evolution Code isochrones model. Our results show that the median mass and luminosity of the 2946 OB-type stars are 5.4 M ⊙ and log(L/L ⊙) = 3.2 with median relative errors of 21.4% and 71.1%, respectively. A good agreement between our results estimated by using our method and those derived by using the orbital motions of binary stars from the literature is found for some B-type stars. In addition, we also fit the mass–luminosity relation of B-type stars by using our derived mass and the luminosity from Gaia Data Release 3.

Convective-core overshooting and the final fate of massive stars

A massive star can explode in powerful supernova (SN) and form a neutron star, but it may also collapse directly into a black hole. Understanding and predicting the final fate of such stars is increasingly important, for instance, in the context of gravitational-wave astronomy. The interior mixing of stars (in general) and convective boundary mixing (in particular) remain some of the largest uncertainties in their evolution. Here, we investigate the influence of convective boundary mixing on the pre-SN structure and explosion properties of massive stars. Using the 1D stellar evolution code MESA, we modeled single, non-rotating stars of solar metallicity, with initial masses of 5 − 70 M⊙ and convective core step-overshooting of 0.05 − 0.50 pressure scale heights. Stars were evolved until the onset of iron core collapse and the pre-SN models were exploded using a parametric, semi-analytic SN code. We used the compactness parameter to describe the interior structure of stars at core collapse and we found a pronounced peak in compactness at carbon-oxygen core masses of MCO ≈ 7 M⊙, along with generally high compactness at MCO ≳ 14 M⊙. Larger convective core overshooting will shift the location of the compactness peak by 1 − 2 M⊙ to higher MCO. These core masses correspond to initial masses of 24 M⊙ (19 M⊙) and ≳40 M⊙ (≳30 M⊙), respectively, in models with the lowest (highest) convective core overshooting parameter. In both high-compactness regimes, stars are found to collapse into black holes. As the luminosity of the pre-supernova progenitor is determined by MCO, we predict black hole formation for progenitors with luminosities of 5.35 ≤ log(L/L⊙)≤5.50 and log(L/L⊙)≥5.80. The luminosity range of black hole formation from stars in the compactness peak is in good agreement with the observed luminosity of the red supergiant star N6946 BH1, which disappeared without a bright supernova, indicating that it had likely collapsed into a black hole. While some of our models in the luminosity range of log(L/L⊙) = 5.1 − 5.5 do indeed collapse to form black holes, this does not fully explain the lack of observed SN IIP progenitors at these luminosities. This case specifically refers to the “missing red supergiant” problem. The amount of convective boundary mixing also affects the wind mass loss of stars, such that the lowest black hole masses are 15 M⊙ and 10 M⊙ in our models, with the lowest and highest convective core overshooting parameter, respectively. The compactness parameter, central specific entropy, and iron core mass describe a qualitatively similar landscape as a function of MCO, and we find that entropy is a particularly good predictor of the neutron-star masses in our models. We find no correlation between the explosion energy, kick velocity, and nickel mass production with the convective core overshooting value, but we do see a tight relation with the compactness parameter. Furthermore, we show how convective core overshooting affects the pre-supernova locations of stars in the Hertzsprung–Russell diagram (HRD) and the plateau luminosity and duration of SN IIP light curves.

Betelgeuse as a Merger of a Massive Star with a Companion

We investigate the merger between a 16M ⊙ star, on its way to becoming a red supergiant (RSG), and a 4M ⊙ main-sequence companion. Our study employs three-dimensional hydrodynamic simulations using the state-of-the-art adaptive mesh refinement code Octo-Tiger. The initially corotating binary undergoes interaction and mass transfer, resulting in the accumulation of mass around the companion and its subsequent loss through the second Lagrangian point (L2). The companion eventually plunges into the envelope of the primary, leading to its spin-up and subsequent merger with the helium core. We examine the internal structural properties of the post-merger star, as well as the merger environment and the outflow driven by the merger. Our findings reveal the ejection of approximately ∼0.6 M ⊙ of material in an asymmetric and somewhat bipolar outflow. We import the post-merger stellar structure into the MESA stellar evolution code to model its long-term nuclear evolution. In certain cases, the post-merger star exhibits persistent rapid equatorial surface rotation as it evolves in the H–R diagram toward the observed location of Betelgeuse. These cases demonstrate surface rotation velocities of a similar magnitude to those observed in Betelgeuse, along with a chemical composition resembling that of Betelgeuse. In other cases, efficient rotationally induced mixing leads to slower surface rotation. This pioneering study aims to model stellar mergers across critical timescales, encompassing dynamical, thermal, and nuclear evolutionary stages.

On the origin of accretion bursts in FU Ori

ABSTRACT Accretion luminosity of young star FU Ori (FUOR) increased from undetectable levels to hundreds of L⊙ in 1937 and remains nearly as high at the present time. In a recent paper, we showed how Extreme Evaporation (EE) of a young gas giant planet that migrated to a ∼10 d orbit around the star may power FUOR. However, our model assumed a power-law mass–radius relation for the evaporating planet. Here, we employ a stellar evolution code to model mass losing planets. We find that adiabatic planets expand rapidly, which results in runaway FUOR bursts. Superadiabatic planets contract while losing mass; their outbursts are dimming with time. Long steadily declining bursts such as FUOR require relatively fine-tuned internal planetary structure, which may be rare. More commonly we find that superadiabatic planets contract too rapidly and their EE falters, leading to FUOR burst stutter. This stutter allows a single planet to produce many short repeating bursts, which may be relevant to bursts observed in V346 Nor, V899, and V1647. We compute broad-band spectra of our best-fitting scenario for FUOR. Since the outburst is triggered behind the planet location, the mid-infrared (mid-IR) emission rises many months before the optical, similar to bursts in Gaia-17bpi and Gaia-18dvy. We show that in outbursts powered by the classic thermal instability, mid-IR lags the optical, whereas the dead zone activation models predict mid-IR light precede the optical burst by many years to decades. We comment on the stellar flyby scenario for FUOR.

Tayler–Spruit dynamo simulations for the modeling of radiative stellar layers

Context. Maxwell stresses exerted by dynamo-generated magnetic fields have been proposed as an efficient mechanism to transport angular momentum in radiative stellar layers. Numerical simulations are still needed to understand its trigger conditions and the saturation mechanisms. Aims. The present study follows up on a recent paper where we reported on the first simulations of Tayler-Spruit dynamos. Here we extend the parameter space explored to assess in particular the influence of stratification on the dynamo solutions. We also present numerical verification of theoretical assumptions made previously that were instrumental in deriving the classical prescription for angular momentum transport implemented in stellar evolution codes. Methods. A simplified radiative layer is modeled numerically by considering the dynamics of a stably stratified, differentially rotating, magnetized fluid in a spherical shell. Results. Our simulations display a diversity of magnetic field topologies and amplitudes depending on the flow parameters, including hemispherical solutions. The Tayler-Spruit dynamos reported here are found to satisfy magnetostrophic equilibrium and achieve efficient turbulent transport of angular momentum, following Spruit’s heuristic prediction.

Scalable stellar evolution forecasting

Many astrophysical applications require efficient yet reliable forecasts of stellar evolution tracks. One example is population synthesis, which generates forward predictions of models for comparison with observations. The majority of state-of-the-art rapid population synthesis methods are based on analytic fitting formulae to stellar evolution tracks that are computationally cheap to sample statistically over a continuous parameter range. The computational costs of running detailed stellar evolution codes, such as MESA, over wide and densely sampled parameter grids are prohibitive, while stellar-age based interpolation in-between sparsely sampled grid points leads to intolerably large systematic prediction errors. In this work, we provide two solutions for automated interpolation methods that offer satisfactory trade-off points between cost-efficiency and accuracy. We construct a timescale-adapted evolutionary coordinate and use it in a two-step interpolation scheme that traces the evolution of stars from zero age main sequence all the way to the end of core helium burning while covering a mass range from 0.65 to 300 M⊙. The feedforward neural network regression model (first solution) that we train to predict stellar surface variables can make millions of predictions, sufficiently accurate over the entire parameter space, within tens of seconds on a 4-core CPU. The hierarchical nearest-neighbor interpolation algorithm (second solution) that we hard-code to the same end achieves even higher predictive accuracy, the same algorithm remains applicable to all stellar variables evolved over time, but it is two orders of magnitude slower. Our methodological framework is demonstrated to work on the MESA ISOCHRONES AND STELLAR TRACKS (Choi et al. 2016) data set, but is independent of the input stellar catalog. Finally, we discuss the prospective applications of these methods and provide guidelines for generalizing them to higher dimensional parameter spaces.

ESTIMATING GAMMA-RAY FLUX FROM MILLISECOND PULSARS ORIGINATING IN GLOBULAR CLUSTERS NEAR THE GALACTIC CENTER

In this study, we investigate the contribution of millisecond pulsars (MSPs) to the gamma-ray excess observed in the Galactic Center by analyzing data from high-resolution direct N-body simulations of six globular clusters (GCs) that experience close encounters with the nuclear star cluster. Using the φ-GPU code, we tracked the orbits of individual neutron stars (NSs) formed during the simulations, assuming a fraction of these NSs evolve into MSPs. Our model includes state-of-the-art single stellar evolution code including prescription for neutron star formation. We estimated the gamma-ray flux from these MSPs, considering known values for their gamma-ray emission. Our results show that MSPs originating from the six modeled GCs contribute a small but non-negligible fraction of the observed gamma-ray flux. This finding suggests that the actual gamma-ray flux from MSPs could be much higher when considering the entire population of GCs, potentially significantly contributing to the gamma-ray excess. This study highlights the importance of considering MSPs in the Galactic Center, originating from nearby globular clusters, as a potential source of the observed gamma-ray excess. Future work will involve more sophisticated simulations incorporating binary stellar evolution and comparing the fraction of MSPs in observed GCs to refine our models and improve the accuracy of our estimates.

Thermal pulses with mesa: resolving the third dredge-up

ABSTRACT The Thermally Pulsing- (Super) Asymptotic Giant Branch is a late stage in the evolution of low- and intermediate-mass stars. These stars undergo strong wind mass-loss and diverse nucleosynthesis. Third dredge-up events, that occur following thermal pulses, are responsible for enriching the surfaces of Asymptotic Giant Branch stars, hence an understanding of this process is crucial for constraining galactic chemical evolution. Using a custom numerical scheme, we investigate the temporal and spatial resolution required to resolve the third dredge-up in the 1D stellar evolution code mesa. With mesa’s default controls, the third dredge-up efficiency is underestimated by as much as $\approx 76~{{\ \rm per\ cent}}$. In stars that undergo hot third dredge-up (M ≳ 6 M⊙), the third dredge-up efficiency is overestimated by $\approx 55~{{\ \rm per\ cent}}$. The Thermally Pulsing- (Super) Asymptotic Giant Branch (TP-(S)AGB) evolution is computed for models with initial masses 1 ≤ Mi/M⊙ ≤ 8 at Solar metallicity (Z = 0.014). The minimum initial mass for carbon stars falls in the range 1.5–1.75 M⊙, compatible with observations. The use of mesa for TP-(S)AGB evolution is validated by comparison to the widely used monash models which show good agreement in the maximum third dredge-up efficiency at initial masses Mi > 2 M⊙. We also compare the third dredge-up efficiency in models produced using two independent stellar evolution codes, fruity and aton, which were computed with various differences in input physics including mass-loss, and which exhibit weaker third dredge-up episodes.

The mass range of hot subdwarf B stars from MESA simulations

ABSTRACT Hot subdwarf B (sdB) stars are helium core burning stars that have lost almost their entire hydrogen envelope due to binary interaction. Their assumed canonical mass of ${M_{\mathrm{sdB}}}\sim 0.47$ M⊙ has recently been debated given a broad range found both from observations as well as from the simulations. Here, we revise and refine the mass range for sdBs derived two decades ago with the Eggleton code, using the stellar evolution code mesa, and discuss the effects of metallicity and the inclusion of core overshooting during the main sequence. We find an excellent agreement for low-mass progenitors, up to ∼2.0 M⊙. For stars more massive than ∼2.5 M⊙, we obtain a wider range of sdB masses compared to the simulations from the literature. Our mesa models for the lower metallicity predict, on average, slightly more massive sdBs. Finally, we show the results for the sdB lifetime as a function of sdB mass and discuss the effect this might have in the comparison between simulations and observational samples. This study paves the way for reproducing the observed Galactic mass distribution of sdB binaries.

Neutrinos from pre-supernova in the framework of TQRPA method

Abstract We propose a new method for calculating spectra and luminosities for (anti)neutrinos produced in the pre-supernova environment by weak processes with hot nuclei. It is based on the thermal quasiparticle random phase approximation (TQRPA), that allows microscopic thermodynamically consistent calculations of the weak-interaction response of nuclei at finite temperatures. For realistic representative pre-supernova conditions from the stellar evolution code MESA, we compute (anti)neutrino luminosities and spectra arising from neutral- and charged-current weak reactions with hot 56Fe and compare them with the contribution of thermal processes. We find that the TQRPA approach produces not only a higher total luminosity of electron neutrinos (mainly born in the electron capture reaction), compared to the standard technique based on the large-scale shell model (LSSM) weak-interaction rates, but also a harder neutrino spectrum. Besides, applying the TQRPA and LSSM, we find that in the context of electron antineutrino generation, the neutral-current nuclear de-excitation (ND) process via neutrino-antineutrino pair emission is at least as important as the electron-positron pair annihilation process. We also show that flavor oscillations enhance the high-energy contribution of the ND process to the electron antineutrino flux. This could potentially be important for pre-supernova antineutrino registration by the Earth’s detectors.

A self-synthesized origin for heavy metals in hot subdwarf stars

Context. A number of He-rich hot subdwarf stars present high abundances for trans-iron elements, such as Sr, Y, Zr, and Pb. Diffusion processes are important in hot subdwarf stars and it is generally believed that the high abundances of heavy elements in these peculiar stars are due to the action of radiative levitation. However, during the formation of He-rich hot subdwarf stars, hydrogen can be ingested into the convective zone driven by the He-core flash. It is known that episodes of protons being ingested into He-burning convective zones can lead to neutron-capture processes and the formation of heavy elements. Aims. In this work, we aim to explore, for the first time, whether neutron-capture processes can occur in late He-core flashes taking place in the cores of the progenitors of He-rich hot subdwarfs. We aim to explore the possibility of a self-synthesized origin for the heavy elements observed in some He-rich hot subdwarf stars. Methods. We computed a detailed evolutionary model for a stripped red-giant star using a stellar evolution code with a nuclear network comprising 32 isotopes. Then we post-processed the stellar models in the phase of helium and hydrogen burning using a post-processing nucleosynthesis code with a nuclear network of 1190 species, which allowed us to follow the neutron-capture processes in detail. Results. We find the occurrence of neutron-capture processes in our model, with neutron densities reaching a value of ∼5 × 1012 cm−3. We determined that the trans-iron elements are enhanced in the surface by 1 to 2 dex, as compared to initial compositions. Moreover, the relative abundance pattern [Xi/Fe] produced by neutron-capture processes closely resembles those observed in some He-rich hot subdwarf stars, hinting at a possible self-synthesized origin for the heavy elements in these stars. Conclusions. We conclude that intermediate neutron-capture processes can occur during a proton ingestion event in the He-core flash of stripped red-giant stars. This mechanism offers a natural channel for the production of the heavy elements observed in certain He-rich hot subdwarf stars.

Dimming the Lights: 2D Simulations of Deflagrations of Hybrid C/O/Ne White Dwarfs Using FLASH

The dimmest and most numerous outlier of the Type Ia supernova population, the Type Iax event, is increasingly being found in the results of observational campaigns. There is currently no single accepted model to describe these events. This 2D study explores the viability of modeling Type Iax events as a hybrid C/O/Ne white dwarf progenitor undergoing a deflagration using the multiphysics software FLASH. This hybrid was created using the stellar evolution code MESA, and its C-depleted core and mixed structure have demonstrated lower yields than traditional C/O progenitors in previous deflagration-to-detonation studies. To generate a sample, 30 “realizations” of this simulation were performed, the only difference being the shape of the initial match head used to start the deflagration. Consistent with earlier work, these realizations produce the familiar hot dense bound remnant surrounded by sparse ejecta. Our results indicate that the majority of the star remains unburned (∼70%) and bound (>90%). Our realizations produce total ejecta yields on the order of 10−2–10−1 M ☉, ejected 56Ni yields on the order of 10−4–10−2 M ☉, and ejecta kinetic energies on the order of 1048–1049 erg. Compared to yields inferred from recent observations of the dimmest Type Iax events—SN 2007qd, SN 2008ha, SN 2010ae, SN 2019gsc, SN 2019muj, SN 2020kyg, and SN 2021fcg—our simulation produces comparable 56Ni yields but too-small total yields and kinetic energies. Reignition of the remnant is also seen in some realizations.

Very massive star models

Context. In addition to being spectacular objects, very massive stars (VMSs) are suspected to have a tremendous impact on their environment and on cosmic evolution in general. The nucleosynthesis both during their advanced stages and their final explosion may contribute greatly to the overall enrichment of the Universe. Their resulting supernovae are candidates for the most superluminous events possible and their extreme conditions also lead to very important radiative and mechanical feedback effects, from local to cosmic scale. Aims. We explore the impact of rotation and metallicity on the evolution of VMSs over cosmic time. Methods. With the recent implementation of an equation of state in the GENEC stellar evolution code, which is appropriate for describing the conditions in the central regions of very massive stars in their advanced phases, we present new results on VMS evolution from Population III to solar metallicity. Results. Low-metallicity VMS models are highly sensitive to rotation, while the evolution of higher-metallicity models is dominated by mass-loss effects. The mass loss strongly affects their surface velocity evolution, breaking quickly at high metallicity while reaching the critical velocity for low-metallicity models. Comparison to observed VMSs in the LMC shows that the mass-loss prescriptions used for these models are compatible with observed mass-loss rates. In our framework for modeling rotation, our models of VMS need a high initial velocity in order to reproduce the observed surface velocities. The surface enrichment of these VMSs is difficult to explain with only one initial composition, and could suggest multiple populations in the R136 cluster. At a metallicity typical of R136, only our non- or slowly rotating VMS models may produce pair-instability supernovae. The most massive black holes that can be formed are less massive than about 60 M⊙. Conclusions. Direct observational constraints on VMS are still scarce. Future observational campaigns will hopefully gather more pieces of information to guide the theoretical modeling of these objects, whose impacts can be very important. VMS tables are available at the CDS.

The Effect of Irradiation on the Spin of Millisecond Pulsars

A millisecond pulsar (MSP) is an old neutron star (NS) that has accreted material from its companion star, causing it to spin up, which is known as the recycling scenario. During the mass transfer phase, the system manifests itself as an X-ray binary. PSR J1402+13 is an MSP with a spin period of 5.89 ms and a spin period derivative of . These properties make it a notable object within the pulsar population, as MSPs typically exhibit low spin period derivatives. In this paper, we aim to explain how an MSP can possess a high spin period derivative by binary evolution. By utilizing the stellar evolution code MESA, we examine the effects of irradiation on the companion star and the propeller effect on the NS during binary evolution. We demonstrate that irradiation can modify the spin period and mass of an MSP, resulting in a higher spin period derivative. These results suggest that the irradiation effect may serve as a key factor in explaining MSPs with high spin period derivatives.

Exploring the Red Supergiant wind kink

Red supergiants (RSG) are key objects in studying the evolution of massive stars and their endpoints, but uncertainties related to their underlying mass-loss mechanism have stood in the way of an appropriate framework for massive star evolution thus far. In this work, we analyse a recently uncovered empirical mass-loss “kink” feature and we highlight its similarity to hot star radiation-driven wind models and observations at the optically thin-to-thick transition point. We motivate a new RSG mass-loss prescription that depends on the Eddington factor, Γ, (including both a steep luminosity, L, dependence and an inverse steep mass, Mcur, dependence). We subsequently implement this new RSG mass-loss prescription in the stellar evolution code MESA. We find that our physically motivated mass-loss behaviour naturally reproduces the Humphreys-Davidson limit without the need for any ad hoc tweaks. It also resolves the RSG supernova “problem”. We argue that a universal behaviour that is seen for radiation-driven winds across the HR diagram, independent of the exact source of opacity, is a key feature of the evolution of the most massive stars.

Evolutionary Origins of Binary Neutron Star Mergers: Effects of Common Envelope Efficiency and Metallicity

The formation histories of compact binary mergers, especially stellar-mass binary black hole mergers, have recently come under increased scrutiny and revision. We revisit the question of the dominant formation channel and efficiency of forming binary neutron star (BNS) mergers. We use the stellar and binary evolution code MESA and implement a detailed method for common envelope and mass transfer. We perform simulations for donor masses between 7 M ⊙ and 20 M ⊙ with a neutron star (NS) companion of 1.4 M ⊙ and 2.0 M ⊙ at two metallicities, using varying common envelope efficiencies and two different prescriptions to determine if the donor undergoes core collapse or electron capture, given their helium and carbon–oxygen cores. In contrast to the case of binary black hole mergers, for an NS companion of 1.4 M ⊙, all BNS mergers are formed following a common envelope phase. For an NS mass of 2.0 M ⊙, we identify a small subset of mergers following only stable mass transfer if the NS receives a natal kick sampled from a Maxwellian distribution with velocity dispersion σ = 265 km s−1. Regardless of the supernova prescription, we find more BNS mergers at subsolar metallicity compared to solar.

Comparative Study of Low-temperature Opacities with GARSTEC Models

We present a comparative study of the effect of low-temperature opacities on stellar models up to the Red Giant branch (RGB), computed with the GARching STellar Evolution Code. We have used two sets of low-temperature opacities; ÆSOPUS (Æ) from the University of Padova and those from the Wichita State University group (F05). In the relevant range of temperatures for this study, < . Therefore, to compare stellar evolutionary tracks, we performed a solar calibration of the α mlt, for each set of low-temperature opacities. After carrying such a calibration, we find that stellar evolutionary tracks are almost unaffected by the choice of low-temperature opacities, with largest variations of 25–30 K at the latest evolutionary stages of the RGB phase.