When modeling and interpreting the spectral energy distributions (SEDs) of galaxies, the simple stellar population (SSP) model, star formation history (SFH), and dust attenuation law (DAL) are three ...of the most important components. However, each of them carries significant uncertainties that have seriously limited our ability to reliably recover the physical properties of galaxies from the analysis of their SEDs. In this paper, we present a Bayesian framework to deal with these uncertain components simultaneously. Based on the Bayesian evidence, a quantitative implement of the principle of Occam's razor, the method allows a more objective and quantitative discrimination among the different assumptions about these uncertain components. With a Ks-selected sample of 5467 low-redshift (mostly with z 1) galaxies in the COSMOS/UltraVISTA field and classified into passively evolving galaxies (PEGs) and star-forming galaxies (SFGs) with the UVJ diagram, we present a Bayesian discrimination of a set of 16 SSP models from five research groups (BC03 and CB07, M05, GALEV, Yunnan-II, BPASS V2.0), five forms of SFH (Burst, Constant, Exp-dec, Exp-inc, Delayed-τ), and four kinds of DAL (Calzetti law, MW, LMC, SMC). We show that the results obtained with the method are either obvious or understandable in the context of stellar/galaxy physics. We conclude that the Bayesian model comparison method, especially that for a sample of galaxies, is very useful for discriminating the different assumptions in the SED modeling of galaxies. The new version of the BayeSED code, which is used in this work, is publicly available at https://bitbucket.org/hanyk/bayesed/.
Abstract
He accretion on to carbon–oxygen white dwarfs (CO WDs) plays a fundamental role when studying the formation of Type Ia supernovae (SNe Ia). Employing the mesa stellar evolution code, we ...calculated the long-term evolution of He-accreting CO WDs. Previous studies usually supposed that a WD can grow in mass to the Chandrasekhar limit in the stable He burning region and finally produce an SN Ia. However, in this study, we find that off-centre carbon ignition occurs in the stable He burning region if the accretion rate is above a critical value (∼2.05 × 10−6 M⊙ yr−1), resulting in accretion-induced collapse rather than an SN Ia. If the accretion rate is below the critical value, explosive carbon ignition will eventually happen in the centre producing an SN Ia. Taking into account the possibility of off-centre carbon ignition, we have re-determined the initial parameter space that produces SNe Ia in the He star donor channel, one of the promising channels to produce SNe Ia in young populations. Since this parameter space is smaller than was found in the previous study of Wang et al. (2009), the SN Ia rates are also correspondingly smaller. We also determined the chemical abundance profile of the He-accreting WDs at the moment of explosive carbon ignition, which can be used as initial input for SN Ia explosion models.
Binary evolution is indispensable in stellar evolution to understand the formation and evolution of most peculiar and energetic objects, such as binary compact objects, Type Ia supernovae, X-ray ...binaries, cataclysmic variables, blue stragglers, hot subdwarfs, and central binaries in planetary nebulae. Mass transfer in binary stars can change the evolutionary path and fate of the components compared to what is expected from single stellar evolution. The critical mass ratio at which unstable mass transfer occurs is an unsolved fundamental problem in binary evolution. To resolve this issue, we construct the thermal equilibrium mass-loss model and derive critical mass ratios for both thermal-timescale mass transfer and unstable mass transfer, the latter of which occurs when the outer Lagrangian point, L2, is overfilled. Using several 3.2 M stellar models as examples, we study the stellar response to thermal equilibrium mass loss and present the thresholds for thermal-timescale mass transfer. We study the possible mass-transfer channels of binary systems containing a 3.2 M donor star, taking into account thermal-timescale mass transfer, unstable mass transfer through L2, and dynamical-timescale mass transfer. We repeat this simulation for a grid of donor stars with different masses (from 0.1 to 100 M with Z = 0.02) and at different evolutionary stages, and present our results. The results show that unstable mass transfer due to the overfilling of the outer Lagrangian point may also play an essential role in the formation of common envelopes for late red giant branch and asymptotic giant branch donors.
ABSTRACT In the limit of extremely rapid mass transfer, the response of a donor star in an interacting binary becomes asymptotically one of adiabatic expansion. We survey here adiabatic mass loss ...from Population I stars (Z = 0.02) of mass 0.10 M -100 M from the zero-age main sequence to the base of the giant branch, or to central hydrogen exhaustion for lower main sequence stars. The logarithmic derivatives of radius with respect to mass along adiabatic mass-loss sequences translate into critical mass ratios for runaway (dynamical timescale) mass transfer, evaluated here under the assumption of conservative mass transfer. For intermediate- and high-mass stars, dynamical mass transfer is preceded by an extended phase of thermal timescale mass transfer as the star is stripped of most of its envelope mass. The critical mass ratio qad (throughout this paper, we follow the convention of defining the binary mass ratio as q Mdonor/Maccretor) above which this delayed dynamical instability occurs increases with advancing evolutionary age of the donor star, by ever-increasing factors for more massive donors. Most intermediate- or high-mass binaries with nondegenerate accretors probably evolve into contact before manifesting this instability. As they approach the base of the giant branch, however, and begin developing a convective envelope, qad plummets dramatically among intermediate-mass stars, to values of order unity, and a prompt dynamical instability occurs. Among low-mass stars, the prompt instability prevails throughout main sequence evolution, with qad declining with decreasing mass, and asymptotically approaching qad = 2/3, appropriate to a classical isentropic n = 3/2 polytrope. Our calculated qad values agree well with the behavior of time-dependent models by Chen & Han of intermediate-mass stars initiating mass transfer in the Hertzsprung gap. Application of our results to cataclysmic variables, as systems that must be stable against rapid mass transfer, nicely circumscribes the range in qad as a function of the orbital period in which they are found. These results are intended to advance the verisimilitude of population synthesis models of close binary evolution.
Eclipsing binary millisecond pulsars (MSPs; the so-called black widows and redbacks) can provide important information about accretion history, pulsar irradiation of their companion stars, and the ...evolutionary link between accreting X-ray pulsars and isolated MSPs. However, the formation of such systems is not well understood, nor the difference in progenitor evolution between the two populations of black widows and redbacks. Whereas both populations have orbital periods between 0.1 and 1.0 days, their companion masses differ by an order of magnitude. In this paper, we investigate the formation of these systems via the evolution of converging low-mass X-ray binaries by employing the MESA stellar evolution code. Our results confirm that one can explain the formation of most of these eclipsing binary MSPs using this scenario. More notably, we find that the determining factor for producing either black widows or redbacks is the efficiency of the irradiation process, such that the redbacks absorb a larger fraction of the emitted spin-down energy of the radio pulsar (resulting in more efficient mass loss via evaporation) compared to that of the black widow systems. We argue that geometric effects (beaming) are responsible for the strong bimodality of these two populations. Finally, we conclude that redback systems do not evolve into black widow systems with time.
We present BayeSED, a general purpose tool for Bayesian analysis of spectral energy distributions (SEDs) using pre-existing model SED libraries or their linear combinations. The artificial neural ...networks, principal component analysis, and multimodal-nested sampling (MultiNest) techniques are employed to allow the highly efficient sampling of posterior distribution and the calculation of Bayesian evidence. As a demonstration, we apply this tool to a sample of hyperluminous infrared galaxies (HLIRGs). The Bayesian evidence obtained for a pure starburst, a pure active galactic nucleus (AGN), and a linear combination of starburst+AGN models show that the starburst+AGN model has the highest evidence for all galaxies in this sample. The Bayesian evidence for the three models and the estimated contributions of starbursts and AGNs to infrared luminosity show that HLIRGs can be classified into two groups: one dominated by starbursts and the other dominated by AGNs. Other parameters and corresponding uncertainties about starbursts and AGNs are also estimated using the model with the highest Bayesian evidence. We find that the starburst region of the HLIRGs dominated by starbursts tends to be more compact and has a higher fraction of OB stars than that of HLIRGs dominated by AGNs. Meanwhile, the AGN torus of the HLIRGs dominated by AGNs tends to be more dusty than that of HLIRGs dominated by starbursts. These results are consistent with previous researches, but need to be tested further with larger samples. Overall, we believe that BayeSED could be a reliable and efficient tool for exploring the nature of complex systems such as dust-obscured starburst-AGN composite systems by decoding their SEDs.
Binary population synthesis shows that mass transfer from a giant star to a main-sequence (MS) companion may account for some observed long-orbital-period blue stragglers. However, little attention ...is paid to this blue straggler formation scenario as dynamical instability often happens when the mass donor is a giant star. In this paper, we have studied the critical mass ratio, qc, for dynamically stable mass transfer from a giant star to a MS companion using detailed evolution calculations. The results show that a more evolved star is generally less stable for Roche lobe overflow. Meanwhile, qc almost linearly increases with the amount of the mass and angular momentum lost during mass transfer, but has little dependance on stellar wind. To conveniently use the result, we give a fit of qc as a function of the stellar radius at the onset of Roche lobe overflow and of the mass-transfer efficiency during the Roche lobe overflow. To examine the formation of blue stragglers from mass transfer between giants and MS stars, we have performed Monte Carlo simulations with various qc. The simulations show that some binaries with the mass donor on the first giant branch may contribute to blue stragglers with qc obtained in this paper but will not from previous qc. Meanwhile, from our qc, blue stragglers from the mass transfer between an asymptotic giant branch star and a MS companion may be more numerous and have a wider range of orbital periods than those from the other qc.
The asymptotic response of donor stars in interacting binary systems to very rapid mass loss is characterized by adiabatic expansion throughout their interiors. In this limit, energy generation and ...heat flow through the stellar interior can be neglected. We model this response by constructing model sequences, beginning with a donor star filling its Roche lobe at an arbitrary point in its evolution, holding its specific entropy and composition profiles fixed as mass is removed from the surface. The stellar interior remains in hydrostatic equilibrium. Luminosity profiles in these adiabatic models of mass-losing stars can be reconstructed from the specific entropy profiles and their gradients. These approximations are validated by comparison with time-dependent binary mass transfer calculations. We describe how adiabatic mass-loss sequences can be used to quantify threshold conditions for dynamical timescale mass transfer, and to establish the range of post-common envelope binaries that are allowed energetically. In dynamical timescale mass transfer, the adiabatic response of the donor star drives it to expand beyond its Roche lobe, leading to runaway mass transfer and the formation of a common envelope with its companion star. For donor stars with surface convection zones of any significant depth, this runaway condition is encountered early in mass transfer, if at all; but for main-sequence stars with radiative envelopes, it may be encountered after a prolonged phase of thermal timescale mass transfer, a so-called delayed dynamical instability. We identify the critical binary mass ratio for the onset of dynamical timescale mass transfer as that ratio for which the adiabatic response of the donor star radius to mass loss matches that of its Roche lobe at some point during mass transfer; if the ratio of donor to accretor masses exceeds this critical value, dynamical timescale mass transfer ensues. In common envelope evolution, the dissipation of orbital energy of the binary provides the energy to eject the common envelope; the energy budget for this process consists essentially of the initial orbital energy of the binary and the initial self-energies of the binary components. We emphasize that, because the stellar core and envelope contribute mutually to each other's gravitational potential energy, proper evaluation of the total energy of a star requires integration over the entire stellar interior, and not just over the ejected envelope alone as commonly assumed. We show that the change in total energy of the donor star, as a function of its remaining mass along an adiabatic mass-loss sequence, can be calculated either by integration over initial and final models, or by a path integral along the mass-loss sequence. That change in total energy of the donor star, combined with the requirement that both remnant donor and its companion star fit within their respective Roche lobes, then circumscribes energetically possible survivors of common envelope evolution.