No Arabic abstract
We developed a new population synthesis code for groups of massive stars, where we model the emission of different forms of energy and matter from the stars of the association. In particular, the ejection of the two radioactive isotopes 26Al and 60Fe is followed, as well as the emission of hydrogen ionizing photons, and the kinetic energy of the stellar winds and supernova explosions. We investigate various alternative astrophysical inputs and the resulting output sensitivities, especially effects due to the inclusion of rotation in stellar models. As the aim of the code is the application to relatively small populations of massive stars, special care is taken to address their statistical properties. Our code incorporates both analytical statistical methods applicable to small populations, as well as extensive Monte Carlo simulations. We find that the inclusion of rotation in the stellar models has a large impact on the interactions between OB associations and their surrounding interstellar medium. The emission of 26Al in the stellar winds is strongly enhanced, compared to non-rotating models with the same mass-loss prescription. This compensates the recent reductions in the estimates of mass-loss rates of massive stars due to the effects of clumping. Despite the lower mass-loss rates, the power of the winds is actually enhanced for rotating stellar models. The supernova power (kinetic energy of their ejecta) is decreased due to longer lifetimes of rotating stars, and therefore the wind power dominates over supernova power for the first 6 Myr after a burst of star-formation. For populations typical of nearby star-forming regions, the statistical uncertainties are large and clearly non-Gaussian.
Employing the the stellar evolution code (Modules for Experiments in Stellar Astrophysics), we calculate yields of heavy elements from massive stars via stellar wind and core-collapse supernovae (CCSN) ejecta to interstellar medium (ISM). In our models, the initial masses ($M_{rm ini}$) of massive stars are taken from 13 to 80 $M_odot$, their initial rotational velocities (V) are 0, 300 and 500 km s$^{-1}$, and their metallicities are [Fe/H] = -3, -2, -1, and 0. The yields of heavy elements coming from stellar winds are mainly affected by the stellar rotation which changes the chemical abundances of stellar surfaces via chemically homogeneous evolution, and enhances mass-loss rate. We estimate that the stellar wind can produce heavy element yields of about $10^{-2}$ (for low metallicity models) to several $M_odot$ (for low metallicity and rapid rotation models) mass. The yields of heavy element produced by CCSN ejecta also depend on the remnant mass of massive mass which is mainly determined by the mass of CO-core. Our models calculate that the yields of heavy elements produced by CCSN ejecta can get up to several $M_odot$. Compared with stellar wind, CCSN ejecta has a greater contribution to the heavy elements in ISM. We also compare the $^{56}$Ni yields by calculated in this work with observational estimate. Our models only explain the $^{56}$Ni masses produced by faint SNe or normal SNe with progenitor mass lower than about 25 $M_odot$, and greatly underestimate the $^{56}$Ni masses produced by stars with masses higher than about 30 $M_odot$.
For understanding the process of star formation it is essential to know how many stars are formed as singles or in multiple systems, as a function of environment and binary parameters. This requires a characterization of the primordial binary population, which we define as the population of binaries that is present just after star formation has ceased, but before dynamical and stellar evolution have significantly altered its characteristics. In this article we present the first results of our adaptive optics survey of 200 (mainly) A-type stars in the nearby OB association Sco OB2. We report the discovery of 47 new candidate companions of Sco OB2 members. The next step will be to combine these observations with detailed simulations of young star clusters, in order to find the primordial binary population.
The Interstellar Medium (ISM) comprises gases at different temperatures and densities, including ionized, atomic, molecular species, and dust particles. The neutral ISM is dominated by neutral hydrogen and has ionization fractions up to 8%. The concentration of chemical elements heavier than helium (metallicity) spans orders of magnitudes in Galactic stars, because they formed at different times. Instead, the gas in the Solar vicinity is assumed to be well mixed and have Solar metallicity in traditional chemical evolution models. The ISM chemical abundances can be accurately measured with UV absorption-line spectroscopy. However, the effects of dust depletion, which removes part of the metals from the observable gaseous phase and incorporates it into solid grains, have prevented, until recently, a deeper investigation of the ISM metallicity. Here we report the dust-corrected metallicity of the neutral ISM measured towards 25 stars in our Galaxy. We find large variations in metallicity over a factor of 10 (with an average 55 +/- 7% Solar and standard deviation 0.28 dex) and including many regions of low metallicity, down to ~17% Solar and possibly below. Pristine gas falling onto the disk in the form of high-velocity clouds can cause the observed chemical inhomogeneities on scales of tens of pc. Our results suggest that this low-metallicity accreting gas does not efficiently mix into the ISM, which may help us understand metallicity deviations in nearby coeval stars.
It has been commonly conjectured that all massive >10 Msun stars are born in OB associations or clusters. Many O and B stars in the Galaxy or the Magellanic Clouds appear to exist in isolation, however. While some of these field OB stars have been ejected from their birthplaces, some are too far away from massive star forming regions to be runaways. Can massive stars form in isolation? The Spitzer survey of the Large Magellanic Cloud (aka SAGE) provides a unique opportunity for us to investigate and characterize the formation sites of massive stars for an entire galaxy. We have identified all massive young stellar objects (YSOs) in the Large Magellanic Cloud. We find that ~85% of the massive YSOs are in giant molecular clouds and ~65% are in OB associations. Only ~7% of the massive YSOs are neither in OB associations nor in giant molecular clouds. This fraction of isolated massive stars in the Large Magellanic Cloud is comparable to the 5-10% found in the Galaxy.
We assemble a census of the most massive stars in Orion, then use stellar isochrones to estimate their masses and ages, and use these results to establish the stellar content of Orions individual OB associations. From this, our new population synthesis code is utilized to derive the history of the emission of UV radiation and kinetic energy of the material ejected by the massive stars, and also follow the ejection of the long-lived radioactive isotopes 26Al and 60Fe. In order to estimate the precision of our method, we compare and contrast three distinct representations of the massive stars. We compare the expected outputs with observations of 26Al gamma-ray signal and the extent of the Eridanus cavity. We find an integrated kinetic energy emitted by the massive stars of 1.8(+1.5-0.4)times 10^52 erg. This number is consistent with the energy thought to be required to create the Eridanus superbubble. We also find good agreement between our model and the observed 26Al signal, estimating a mass of 5.8(+2.7-2.5) times 10^-4 Msol of 26Al in the Orion region. Our population synthesis approach is demonstrated for the Orion region to reproduce three different kinds of observable outputs from massive stars in a consistent manner: Kinetic energy as manifested in ISM excavation, ionization as manifested in free-free emission, and nucleosynthesis ejecta as manifested in radioactivity gamma-rays. The good match between our model and the observables does not argue for considerable modifications of mass loss. If clumping effects turn out to be strong, other processes would need to be identified to compensate for their impact on massive-star outputs. Our population synthesis analysis jointly treats kinematic output and the return of radioactive isotopes, which proves a powerful extension of the methodology that constrains feedback from massive stars.