Stripped-envelope supernovae (Type IIb, Ib, Ic) showing little or no hydrogen are one of the main classes of explosions of massive stars. Their origin and the evolution of their progenitors are not fully understood as yet. Very massive single stars s
tripped by their own winds ($gtrsim 25-30 M_{odot}$ at solar metallicity) are considered viable progenitors of these events. However, recent 1D core-collapse simulations show that some massive stars may collapse directly onto black holes after a failed explosion, with weak or no visible transient. In this letter, we estimate the effect of direct collapse onto a black hole on the rates of stripped-envelope supernovae that arise from single stars. For this, we compute single star MESA models at solar metallicity and map their final state to their core-collapse outcome following prescriptions commonly used in population synthesis. According to our models, no single stars that have lost their entire hydrogen-rich envelope are able to explode, and only a fraction of progenitors with a thin hydrogen envelope left (IIb progenitor candidates) do, unless we invoke increased wind mass-loss rates. This result increases the existing tension between the single-star scenario for stripped-envelope supernovae and their observed rates and properties. At face value, our results point towards an even higher contribution of binary progenitors for stripped-envelope supernovae. Alternatively, they may suggest inconsistencies in the common practice of mapping different stellar models to core-collapse outcomes and/or higher overall mass loss in massive stars.
The origin of nebular HeII emission, which is frequently observed in low-metallicity (O/H) star-forming galaxies, remains largely an unsolved question. Using the observed anticorrelation of the integrated X-ray luminosity per unit of star formation r
ate ($L_X/{rm SFR}$) of an X-ray binary population with metallicity and other empirical data from the well-studied galaxy I Zw 18, we show that the observed HeII 4686 intensity and its trend with metallicity is naturally reproduced if the bulk of He$^+$ ionizing photons are emitted by the X-ray sources. We also show that a combination of X-ray binary population models with normal single and/or binary stellar models reproduces the observed $I(4686)/I(Hbeta)$ intensities and its dependency on metallicity and age. We conclude that both empirical data and theoretical models suggest that high-mass X-ray binaries are the main source of nebular HeII emission in low-metallicity star-forming galaxies.
We present the largest-scale comparison to date between observed extragalactic X-ray binary (XRB) populations and theoretical models of their production. We construct observational X-ray luminosity functions (oXLFs) using Chandra observations of 12 l
ate-type galaxies from the Spitzer Infrared Nearby Galaxy Survey (SINGS). For each galaxy, we obtain theoretical XLFs (tXLFs) by combining XRB synthetic models, constructed with the population synthesis code StarTrack, with observational star formation histories (SFHs). We identify highest-likelihood models both for individual galaxies and globally, averaged over the full galaxy sample. Individual tXLFs successfully reproduce about half of oXLFs, but for some galaxies we are unable to find underlying source populations, indicating that galaxy SFHs and metallicities are not well matched and/or XRB modeling requires calibration on larger observational samples. Given these limitations, we find that best models are consistent with a product of common envelope ejection efficiency and central donor concentration ~=0.1, and a 50% uniform -- 50% twins initial mass-ratio distribution. We present and discuss constituent subpopulations of tXLFs according to donor, accretor and stellar population characteristics. The galaxy-wide X-ray luminosity due to low-mass and high-mass XRBs, estimated via our best global model tXLF, follows the general trend expected from the L_x - star formation rate and L_x - stellar mass relations of Lehmer et al 2010. Our best models are also in agreement with modeling of the evolution both of XRBs over cosmic time and of the galaxy X-ray luminosity with redshift.
In black hole X-ray binaries, a misalignment between the spin axis of the black hole and the orbital angular momentum can occur during the supernova explosion that forms the compact object. In this letter we present population synthesis models of Gal
actic black hole X-ray binaries, and study the probability density function of the misalignment angle, and its dependence on our model parameters. In our modeling, we also take into account the evolution of misalignment angle due to accretion of material onto the black hole during the X-ray binary phase. The major factor that sets the misalignment angle for X-ray binaries is the natal kick that the black hole may receive at its formation. However, large kicks tend to disrupt binaries, while small kicks allow the formation of XRBs and naturally select systems with small misalignment angles. Our calculations predict that the majority (>67%) of Galactic field BH XRBs have rather small (>10 degrees) misalignment angles, while some systems may reach misalignment angles as high as ~90 degrees and even higher. This results is robust among all population synthesis models. The assumption of small small misalignment angles is extensively used to observationally estimate black hole spin magnitudes, and for the first time we are able to confirm this assumption using detailed population synthesis calculations.
We propose a physically motivated and self-consistent prescription for the modeling of transient neutron star (NS) low-mass X-ray binary (LMXB) properties, such as duty cycle (DC), outburst duration and recurrence time. We apply this prescription to
the population synthesis (PS) models of field LMXBs presented by Fragos et al. (2008), and compare the transient LMXB population to the Chandra X-ray survey of the two elliptical galaxies NGC 3379 and NGC 4278, which revealed several transient sources (Brassington et al., 2008, 2009). We are able to exclude models with a constant DC for all transient systems, while models with a variable DC based on the properties of each system are consistent with the observed transient populations. We predict that the majority of the observed transient sources in these two galaxies are LMXBs with red giant donors. Our comparison suggests that LMXBs formed through evolution of primordial field binaries are dominant in globular cluster (GC) poor elliptical galaxies, while they still have a significant contribution in GC rich ones.
In recent years, an increasing number of proper motions have been measured for Galactic X-ray binaries. When supplemented with accurate determinations of the component masses, orbital period, and donor effective temperature, these kinematical constra
ints harbor a wealth of information on the systems past evolution. Here, we consider all this available information to reconstruct the full evolutionary history of the black hole X-ray binary XTE J1118+480, assuming that the system originated in the Galactic disk and the donor has solar metallicity. This analysis accounts for four evolutionary phases: mass transfer through the ongoing X-ray phase, tidal evolution before the onset of Roche-lobe overflow, motion through the Galactic potential after the formation of the black hole, and binary orbital dynamics due to explosive mass loss and possibly a black hole natal kick at the time of core collapse. We find that right after black hole formation, the system consists of a ~6.0-10.0 solar masses black hole and a ~1.0-1.6 solar masses main-sequence star. We also find that that an asymmetric natal kick is not only plausible but required for the formation of this system, and derive a lower and upper limit on the black hole natal kick velocity magnitude of 80 km/s and 310 km/s, respectively.
We present theoretical models for the formation and evolution of populations of low-mass X-ray binaries (LMXB) in the two elliptical galaxies NGC 3379 and NGC 4278. The models are calculated with the recently updated StarTrack code (Belczynski et al.
, 2006), assuming only a primordial galactic field LMXB population. StarTrack is an advanced population synthesis code that has been tested and calibrated using detailed binary star calculations and incorporates all the important physical processes of binary evolution. The simulations are targeted to modeling and understanding the origin of the X-ray luminosity functions (XLF) of point sources in these galaxies. For the first time we explore the population XLF down to luminosities of 3X10^36 erg/s, as probed by the most recent observational results (Kim et al., 2006). We consider models for the formation and evolution of LMXBs in galactic fields with different CE efficiencies, stellar wind prescriptions, magnetic braking laws and initial mass functions. We identify models that produce an XLF in excellent agreement with the observations both in shape and absolute normalization. We also find that the treatment of the outburst luminosity of transient systems remains a crucial factor for the determination of the XLF since the modeled populations are dominated by transient X-ray systems.
In recent years, an increasing number of proper motions have been measured for Galactic X-ray binaries. When supplemented with accurate determinations of the component masses, orbital period, and donor luminosity and effective temperature, these kine
matical constraints harbor a wealth of information on the systems past evolution. The constraints on compact object progenitors and kicks derived from this are of immense value for understanding compact object formation and exposing common threads and fundamental differences between black hole and neutron star formation. Here, we present the results of such an analysis for the black hole X-ray binary XTE J1118+480. We present results from modeling the mass transfer phase, following the motion in the Galaxy back to the birth site of the black hole, and examining the dynamics of symmetric and asymmetric core-collapses of the black hole progenitor.
We revisit the formation of PSR J0737-3039, taking into account the most recent observational constraints. We show that the most likely kick velocity and progenitor parameters depend strongly on the consideration of the full five-dimensional PDF for
the magnitude and direction of the kick velocity imparted to pulsar B at birth, the mass of pulsar Bs pre-supernova helium star progenitor, and the pre-supernova orbital separation, and on the adopted prior assumptions. The priors consist of the transverse systemic velocity, the age of the system, and the treatment of the unknown radial velocity. Since the latter cannot be determined from observation, we adopt a statistical approach and use theoretical radial-velocity distributions obtained from population synthesis calculations for coalescing double neutron stars. We find that the prior assumptions about the pre-supernova helium star mass affect the derived most likely parameters significantly: when the minimum helium star mass required for neutron star formation is assumed to be 2.1Msun, the most likely kick velocity ranges from 70-180km/s; when masses lower than 2.1Msun are assumed to allow neutron star formation, the most likely kick velocity can be as low as a few km/s, although the majority of the considered models still yield most likely kick velocities of 50-170km/s. We also show that the proximity of the double pulsar to the Galactic plane and the small proper motion do not pose stringent constraints on the kick velocity and progenitor mass of pulsar B. Instead, the constraints imposed by the orbital dynamics of asymmetric supernova explosions turn out to be much more restrictive. We conclude that the currently available observational constraints cannot be used to favor a specific core-collapse and neutron star formation mechanism. (abridged)
