Wolf-Rayet Galaxies in SDSS-IV MaNGA. II. Metallicity Dependence of the High-mass Slope of the Stellar Initial Mass Function

As hosts of living high-mass stars, Wolf-Rayet (WR) regions or WR galaxies are ideal objects for constraining the high-mass end of the stellar initial mass function (IMF). We construct a large sample of 910 WR galaxies/regions that cover a wide range of stellar metallicity (from Z ∼ 0.001 to 0.03) by combining three catalogs of WR galaxies/regions previously selected from the SDSS and SDSS-IV/MaNGA surveys. We measure the equivalent widths of the WR blue bump at ∼4650 Å for each spectrum. They are compared with predictions from stellar evolutionary models Starburst99 and BPASS, with different IMF assumptions (high-mass slope α of the IMF ranging from 1.0 to 3.3). Both singular evolution and binary evolution are considered. We also use a Bayesian inference code to perform full spectral fitting to WR spectra with stellar population spectra from BPASS as fitting templates. We then make a model selection among different α assumptions based on Bayesian evidence. These analyses have consistently led to a positive correlation of the IMF high-mass slope α with stellar metallicity Z, i.e., with a steeper IMF (more bottom-heavy) at higher metallicities. Specifically, an IMF with α = 1.00 is preferred at the lowest metallicity (Z ∼ 0.001), and an Salpeter or even steeper IMF is preferred at the highest metallicity (Z ∼ 0.03). These conclusions hold even when binary population models are adopted.

Pre-explosion Properties of Helium Star Donors to Thermonuclear Supernovae

Helium star–carbon-oxygen white dwarf (CO WD) binaries are potential single-degenerate progenitor systems of thermonuclear supernovae. Revisiting a set of binary evolution calculations using the stellar evolution code MESA, we refine our previous predictions about which systems can lead to a thermonuclear supernova and then characterize the properties of the helium star donor at the time of explosion. We convert these model properties to near-UV/optical magnitudes assuming a blackbody spectrum and support this approach using a matched stellar atmosphere model. These models will be valuable to compare with pre-explosion imaging for future supernovae, though we emphasize the observational difficulty of detecting extremely blue companions. The pre-explosion source detected in association with SN 2012Z has been interpreted as a helium star binary containing an initially ultra-massive WD in a multiday orbit. However, extending our binary models to initial CO WD masses of up to 1.2 M ⊙, we find that these systems undergo off-center carbon ignitions and thus are not expected to produce thermonuclear supernovae. This tension suggests that, if SN 2012Z is associated with a helium star–WD binary, then the pre-explosion optical light from the system must be significantly modified by the binary environment and/or the WD does not have a carbon-rich interior composition.

The Maximum Accreted Mass of Recycled Pulsars

The maximum mass of neutron stars (NSs) is of great importance for constraining equations of state of NSs and understanding the mass gap between NSs and stellar-mass black holes. NSs in X-ray binaries increase in mass by accreting material from their companions (known as the recycling process), and the uncertainties in the accretion process make studying the NS mass at birth a challenge. In this work, we investigate the NS accreted mass while considering the effect of NS spin evolution and provide the maximum accreted mass for NSs in the recycling process. By exploring a series of binary evolution calculations, we obtain the final NS mass and the maximum accreted mass for a given birth mass of an NS and a mass transfer efficiency. Our results show that NSs can accrete relatively more material for binary systems with donor masses in the range of 1.8 ∼ 2.4 M ⊙, NSs accrete relatively more mass when the remnant WD mass is in the range of ∼ 0.25–0.30 M ⊙, and the maximum accreted mass is positively correlated with the initial NS mass. For a 1.4 M ⊙ NS at birth with a moderate mass transfer efficiency of 0.3, the maximum accreted mass could be 0.27 M ⊙. The results can be used to estimate the minimum birth mass for systems with massive NSs in observations.

Most “Young” α-rich Stars Have High Masses but are Actually Old

Recent observations have revealed a population of α-element abundances, enhanced giant stars with unexpected high masses (≳1 M ⊙) from asteroseismic analysis and spectroscopy. Assuming single-star evolution, their masses imply young ages (τ < 6 Gyr) incompatible with the canonical Galactic chemical evolution scenario. Here we study the chemistry and kinematics of a large sample of such α-rich, high-mass red giant branch (RGB) stars drawn from the LAMOST spectroscopic surveys. Using LAMOST and Gaia, we found these stars share the same kinematics as the canonical high-α old stellar population in the Galactic thick disk. The stellar abundances show that these high-α massive stars have α- and iron-peak element abundances similar to those of the high-α old thick-disk stars. However, a portion of them exhibit higher [(N+C)/Fe] and [Ba/Fe] ratios, which implies they have gained C- and Ba-rich materials from extra sources, presumably asymptotic giant branch (AGB) companions. The results support the previous suggestion that these RGB stars are products of binary evolution. Their high masses thus mimic “young” single stars, yet in fact they belong to an intrinsic old stellar population. To fully explain the stellar abundance patterns of our sample stars, a variety of binary evolution channels, such as main-sequence (MS) + RGB, MS + AGB, RGB + RGB, and RGB + AGB, are required, pointing to diverse formation mechanisms of these seemly rejuvenated cannibals. With this larger sample, our results confirm earlier findings that most, if not all, α-rich stars in the Galactic disk seem to be old.

Simulation of Formation of Binary Compact Objects in Globular Cluster

Binary compact objects are the special binary systems, which were composed of compact objects (i.e., white dwarf, neutron star and black hole). They contribute a lot to the sources of gravitational waves. The study of binary compact objects in star clusters and galaxies can provide a theoretical guidance for gravitational wave detection and improve the probability of detection effectively. We simulate the formation and fraction of binary compact objects in a small globular cluster, via NBODY6++GPU, an efficient N-body simulation code. We obtain the fractions of white dwarf binaries and neutron star binaries at different ages. The results show that the dynamic interactions among stars can also result in some binary compact objects besides binary evolution.

Population synthesis of accreting white dwarfs: Rates and evolutionary pathways of H and He novae

Novae are some of the most commonly detected optical transients and have the potential to provide valuable information about binary evolution. Binary population synthesis codes have emerged as the most effective tool for modelling populations of binary systems, but such codes have traditionally employed greatly simplified nova physics, precluding detailed study. In this work, we implement a model treating H and He novae as individual events into the binary population synthesis code binary_c. This treatment of novae represents a significant improvement on the ‘averaging’ treatment currently employed in modern population synthesis codes. We discuss the evolutionary pathways leading to these phenomena and present nova event rates and distributions of several important physical parameters. Most novae are produced on massive white dwarfs, with approximately 70 and 55 per cent of nova events occurring on O/Ne white dwarfs for H and He novae respectively. Only 15 per cent of H-nova systems undergo a common-envelope phase, but these systems are responsible for the majority of H nova events. All He-accreting He-nova systems are considered post-common-envelope systems, and almost all will merge with their donor star in a gravitational-wave driven inspiral. We estimate the current annual rate of novae in M31 (Andromeda) to be approximately 41 ± 4 for H novae, underpredicting the current observational estimate of $65^{+15}_{-16}$, and 0.14 ± 0.015 for He novae. When varying common-envelope parameters, the H nova rate varies between 20 and 80 events per year.