Source author record

Zhanwen Han

Zhanwen Han appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

89works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

89 published item(s)

preprint2026arXiv

Pre-Supernova Eruptions Triggered by Sudden Energy Deposition in Low-Mass Core-Collapse Supernova Progenitors

In low-mass core-collapse supernova (CCSN) progenitors, nuclear burning beyond oxygen can become explosive under degenerate conditions, triggering eruptive mass loss before the final explosion. We investigate such pre-SN eruptions using \texttt{SNEC} hydrodynamic simulations and realistic stellar models, parameterizing the nuclear energy deposition as a fraction of the binding energy of the combined He layer and H-rich envelope. For the lowest-mass model (9 $M_\odot$), the ejecta mass ($M_{\rm ej}$) scales with the energy gained by the H-rich envelope via a power law (index$\sim$3.5). Across 9-10 $M_\odot$, this relation shows limited scatter within a factor of $\sim$2.6, enabling an estimation of the gained energy from $M_{\rm ej}$. The shock passage also flattens the bound envelope, which can affect the SN light curve morphology and provide another diagnostic for the eruption. Then, we compute the associated precursor light curves for the 9 $M_\odot$ model with the multi-group radiative-transfer code \texttt{STELLA}. These signals are typically faint, with bolometric luminosities of $\sim10^{39}$ erg s$^{-1}$ lasting hundreds of days. Their cool black-body spectra make them brighter in the infrared, yet several magnitudes fainter than observed pre-SN precursors at the threshold for full envelope ejection. To aid future studies, we make our post-eruption stellar profiles and precursor light curves publicly available.

preprint2025arXiv

Introduction to the Chinese Space Station Survey Telescope (CSST)

The Chinese Space Station Survey Telescope (CSST) is an upcoming Stage-IV sky survey telescope, distinguished by its large field of view (FoV), high image quality, and multi-band observation capabilities. It can simultaneously conduct precise measurements of the Universe by performing multi-color photometric imaging and slitless spectroscopic surveys. The CSST is equipped with five scientific instruments, i.e. Multi-band Imaging and Slitless Spectroscopy Survey Camera (SC), Multi-Channel Imager (MCI), Integral Field Spectrograph (IFS), Cool Planet Imaging Coronagraph (CPI-C), and THz Spectrometer (TS). Using these instruments, CSST is expected to make significant contributions and discoveries across various astronomical fields, including cosmology, galaxies and active galactic nuclei (AGN), the Milky Way and nearby galaxies, stars, exoplanets, Solar System objects, astrometry, and transients and variable sources. This review aims to provide a comprehensive overview of the CSST instruments, observational capabilities, data products, and scientific potential.

preprint2025arXiv

The Formation of Electron-capture Supernovae: A Review

It is generally believed that the electron-capture reactions happen when the oxygen-neon (ONe) cores grow in masses close to the Chandrasekhar limit, leading to the formation of neutron stars (NSs) via electron-capture supernovae (EC-SNe). EC-SNe are predicted to be the most likely short-lived and faint optical transients, and a small ejecta mass is expected during the collapse. This kind of SNe provide an alternative channel for producing isolated NSs and NS systems, especially for the formation of X-ray binaries and double NSs. Although EC-SNe were proposed ~45 yr ago, there are still some uncertainties for the origin of EC-SNe. In this article, we review recent studies on the two classic progenitor channels of EC-SNe, i.e., the single star channel and the binary star channel. In the single star channel, EC-SNe can happen in super asymptotic giant branch stars or He stars, whereas in the binary star channel EC-SNe can occur in He stars in binaries (involving He star+MS systems and NS+He star systems) or accretion-induced collapse in white dwarf binaries (involving the single-degenerate scenario and the double-degenerate scenario). Recent progress on the two progenitor channels is discussed, including the initial parameter range for EC-SNe, the evolutionary paths to EC-SNe, related objects and some observational constraints, etc. We also make some discussions on the possible candidates for EC-SNe in this article, and the impacts of EC-SNe on some research fields, e.g., the properties of NSs, double NS population and chemical products, etc. We also discuss the differences between EC-SNe and ultra-stripped SNe in this article. Research on EC-SNe is at a pivotal stage, with key theoretical uncertainties and observational challenges requiring integrated modeling and multi-wavelength observations for robust identification.

preprint2023arXiv

Influence of a mass transfer stability criterion on double white dwarf populations

Mass transfer stability is an essential issue in binary evolution. Ge et al. studied critical mass ratios for dynamically stable mass transfer by establishing adiabatic mass loss model and found that the donor stars on the giant branches tend to be more stable than that based on the composite polytropic stellar model. We would investigate the influence of mass transfer stability on the formation and properties of DWD populations. We performed a series of binary population synthesis, where the critical mass ratios from adiabatic mass loss model (Ge's model) and that from the composite polytropic model are adopted, respectively. For Ge's model, most of the DWDs are produced from the stable non-conservative Roche lobe overflow plus common envelope (CE) ejection channel (RL+CE channel) regardless of the CE ejection efficiency $α_{CE}$. While the results of the polytropic model strongly depend on the adopted value of $α_{ CE}$. We find DWDs produced from the RL+CE channel have comparable WD masses and the mass ratio distribution peaks at around 1. Based on the magnitude-limited sample of DWDs, the space densities for the detectable DWDs and those with extremely low-mass WD (ELM WD) companions in Ge's model is $1347$ and $473 kpc^{-3}$, respectively, close to observations. While the polytropic model overpredicts space density of DWDs by a factor of about $2-3$. We also find that the results of DWD merger rate distribution in Ge's model reproduce the observations better than that of the polytropic model, and the merger rate of DWDs with ELM WD companions in the Galaxy is about $1.8\times 10^{-3} yr^{-1}$ in Ge's model, which is comparable to the observation estimation. We confirm that the mass transfer stability plays important roles in the formation and properties of DWD populations, and then in the progenitors of SNe Ia and detectable GW sources.

preprint2022arXiv

A Roche Lobe-filling hot Subdwarf and White Dwarf Binary: possible detection of an ejected common envelope

Binaries consisting of a hot subdwarf star and an accreting white dwarf (WD) are sources of gravitational wave radiation at low frequencies and possible progenitors of type Ia supernovae if the WD mass is large enough. Here, we report the discovery of the third binary known of this kind: it consists of a hot subdwarf O (sdO) star and a WD with an orbital period of 3.495 hours and an orbital shrinkage of 0.1 s in 6 yr. The sdO star overfills its Roche lobe and likely transfers mass to the WD via an accretion disk. From spectroscopy, we obtain an effective temperature of $T_{\mathrm{eff}}=54\,240\pm1\,840$ K and a surface gravity of $\log{g}=4.841\pm0.108$ for the sdO star. From the light curve analysis, we obtain a sdO mass of $M_{\mathrm{sdO}}=0.55$ ${\mathrm{M_{\odot}}}$ and a mass ratio of $q=M_{\mathrm{WD}}/M_{\mathrm{sdO}}=0.738\pm0.001$. Also, we estimate that the disk has a radius of $\sim 0.41R_\odot$ and a thickness of $\sim 0.18R_\odot$. The origin of this binary is probably a common envelope ejection channel, where the progenitor of the sdO star is either an RGB star or, more likely, an early AGB star; the sdO star will subsequently evolve into a WD and merge with its WD companion, likely resulting in an R CrB star. The outstanding feature in the spectrum of this object is strong Ca H&K lines, which are blueshifted by $\sim$200 km/s and likely originate from the recently ejected common envelope, and we estimated that the remnant CE material in the binary system has a density $\sim 6\times 10^{-10} {\rm g/cm^3}$.

preprint2022arXiv

Evolution of AM CVn binaries with WD donors

The evolution and stability of mass transfer of CO+He WD binaries are not well understood. Observationally they may emerge as AM CVn binaries and are important gravitational wave (GW) emitters. In this work, we have modeled the evolution of double WD binaries with accretor masses of $0.50 - 1.30\;M_{\odot}$ and donor masses of $0.17\; - 0.45\;M_{\odot}$ using the detailed stellar evolution code MESA. We find that the evolution of binaries with same donor masses but different accretor masses is very similar and binaries with same accretor masses but larger He donor masses have larger maximum mass transfer rates and smaller minimum orbital periods. We also demonstrate that the GW signal from AM CVn binaries can be detected by space-borne GW observatories, such as LISA, TianQin. And there is a linear relation between the donor mass and gravitational wave frequency during mass transfer phase. In our calculation, all binaries can have dynamically stable mass transfer, which is very different from previous studies. The threshold donor mass of Eddington-limited mass transfer for a given accretor WD mass is lower than previous studies. Assuming that a binary may enter common envelope if the mass transfer rate exceeds the maximum stable burning rate of He, we provide a new criterion for double WDs surviving mass transfer, which is below the threshold of Eddington-limit. Finally, we find that some systems with ONe WDs in our calculation may evolve into detached binaries consisting of neutron stars (NSs) and extremely low mass He WDs and further ultra-compact X-ray binaries.

preprint2022arXiv

Formation of black widows through ultra-compact X-ray binaries with He star companions

Black widows (BWs) are a type of eclipsing millisecond pulsars (MSPs) with companion masses $M_2$ $\lesssim 0.05\,\rm M_\odot$, which can be used to study the accretion history and the radiation of pulsars, as well as the origin of isolated MSPs. Recent observations indicate that there are two sub-types of BWs. One is the BWs with $M_2$ $\sim$ $0.01-0.05\,\rm M_\odot$, whereas another with $M_2$ $\lesssim 0.01\,\rm M_\odot$. However, the origin of the latter is still highly uncertain. In this paper, we investigated the formation of BWs with $M_2$ $\lesssim 0.01\,\rm M_\odot$ through ultra-compact X-ray binaries (UCXBs) with He star companions, in which a neutron star (NS) accretes material from a He star through Roche-lobe overflow. By considering different He star masses and evaporation efficiencies with the stellar evolution code Modules for Experiments in Stellar Astrophysics (MESA), we evolved a series of NS+He star systems that can undergo UCXB stage. We found that this channel can explain the formation of BWs with $M_2$ $\lesssim 0.01\,\rm M_\odot$ within the Hubble time, especially three widely studied BWs, i.e. PSRs J1719-1438, J2322-2650 and J1311-3430. We also found that X-ray irradiation feedback does not affect the evolutionary tracks of evaporation process. Our simulations indicate that the UCXB channel with He star companions are the potential progenitors of isolated MSPs, and that the origin of BWs with $M_2$ $\lesssim 0.01\,\rm M_\odot$ is different with another sub-type of BWs. In addition, the present work suggests that the BWs with $M_2$ $\lesssim 0.01\,\rm M_\odot$ may not be produced by redback systems.

preprint2022arXiv

Identification of new classical Be stars from the LAMOST MRS survey

Be stars are B-type main-sequence stars that display broad Balmer emission lines in their spectra. Identification of Be population is essential to further examine the formation and evolutionary models. We report the detection of classical Be (CBe) stars from observations with the Large sky Area Multi-Object fiber Spectroscopic Telescope Medium Resolution Survey of Date Release 7 (LAMOST MRS DR7). We used a deep convolutional neural network, the ResNet, with an 18-layer module to examine the morphology of the H alpha profile. We identified 1,162 candidate Be stars from the collection of 2,260,387 spectra for 789,918 stars in the database. The ResNet network achieves a Be star classification accuracy of 99.5%. Among the detections, 151 of these are prior known Be stars cross-matched from the literature. By applying a three-step test, we identified 183 new CBe stars. We find that 41 CBe stars are members of known open clusters. Based upon an investigation of the kinematics of the identified CBe stars from the Gaia EDR3 astrometric solutions, we identified 16 new runaways. These new identifications will provide a reference for future follow-ups to further investigate their physical properties.

preprint2022arXiv

Investigating the stability of mass transfer in neutron star-helium white dwarf binaries

Neutron star-helium white dwarf (NS+He WD) binaries are important evolutionary products of close-orbit binary star systems. They are often observed as millisecond pulsars and may continue evolving into ultra-compact X-ray binaries (UCXBs) and continuous gravitational wave (GW) sources that will be detected by space-borne GW observatories, such as LISA, TianQin and Taiji. Nevertheless, the stability of NS+He WD binaries undergoing mass transfer is not well studied and still under debate. In this paper, we model the evolution of NS+He WD binaries with WD masses ranging from 0.17-0.45 $M_{\odot}$, applying the detailed stellar evolution code mesa. Contrary to previous studies based on hydrodynamics, we find that apparently all NS+He WD binaries undergo stable mass transfer. We find for such UCXBs that the larger the WD mass, the larger the maximum mass-transfer rate and the smaller the minimum orbital period during their evolution. Finally, we demonstrate numerically and analytically that there is a tight correlation between WD mass and GW frequency for UCXBs, independent of NS mass.

preprint2022arXiv

Long-term evolution of post-explosion Helium-star Companions of Type Iax Supernovae

Supernovae of Type Iax (SNe Iax) are an accepted faint subclass of hydrogen-free supernovae. Their origin, the nature of the progenitor systems, however, is an open question. Recent studies suggest that the weak deflagration explosion of a near-Chandrasekhar-mass white dwarf in a binary system with a helium star donor could be the origin of SNe Iax. In this scenario, the helium star donor is expected to survive the explosion. We use the one-dimensional stellar evolution codes \textsc{MESA} and \textsc{Kepler} to follow the post-impact evolution of the surviving helium companion stars. The stellar models are based on our previous hydrodynamical simulations of ejecta-donor interaction, and we explore the observational characteristics of these surviving helium companions. We find that the luminosities of the surviving helium companions increase significantly after the impact: They could vary from $2\mathord,500\,\mathrm{L_{\odot}}$ to $16\mathord,000\,\mathrm{L_{\odot}}$ for a Kelvin-Helmholtz timescale of about $10^{4}\,\mathrm{yr}$. After the star reaches thermal equilibrium, it evolves as an O-type hot subdwarf (sdO) star and continues its evolution along the evolutionary track of a normal sdO star with the same mass. Our results will help to identify the surviving helium companions of SNe Iax in future observations and to place new constraints on their progenitor models.

preprint2022arXiv

Mass-Ratio Distribution of Binaries From the LAMOST-MRS Survey

Binary evolution leads to the formation of important objects crucial to the development of astrophysics, but the statistical properties of binary populations are still poorly understood. The LAMOST-MRS has provided a large sample of stars to study the properties of binary populations, especially for the mass ratio distributions and the binary fractions. We have devised a Peak Amplitude Ratio (PAR) approach to derive the mass ratio of a binary system based on results obtained from its spectrum. By computing a cross-correlation function (CCF), we established a relationship between the derived mass ratio and the PARs of the binary systems. By utilizing spectral observations obtained from LAMSOT DR6 & DR7, we applied the PAR approach to form distributions of the derived mass ratio of the binary systems to the spectral types. We selected the mass ratio within the range of $0.6-1.0$ for investigating the mass-ratio distribution. Through a power-law fitting, we obtained the power index $γ$ values of $-0.42\pm0.27$, $0.03\pm0.12$, and $2.12\pm0.19$ for A-, F-, and G-type stars identified in the sample, respectively. The derived $γ$-values display an increasing trend toward lower primary star masses, and G-type binaries tend to be more in twins. The close binary fractions (for $P\lesssim 150\,{\rm d}$ and $q\gtrsim 0.6$) in our sample for A, F and G binaries are $7.6\pm 0.5 \%$, $4.9\pm 0.2 \%$ and $3.7 \pm 0.1 \%$, respectively. Note that the PAR approach can be applied to large spectroscopic surveys of stars.

preprint2022arXiv

On the rotation properties of a post-explosion helium-star companion in Type Iax supernovae

Recent studies have suggested that type Iax supernovae (SNe Iax) are likely to result from a weak deflagration explosion of a Chandrasekhar-mass white dwarf in a binary system with a helium (He)-star companion. Assuming that most SNe Iax are produced from this scenario, in this work we extend our previous work on the three-dimensional hydrodynamical simulation of ejecta-companion interaction by taking the orbital and spin velocities of the progenitor system into account. We then follow the post-impact evolution of a surviving He-star companion by using the one-dimensional stellar evolution code \textsc{MESA}. We aim to investigate the post-explosion rotation properties of a He-star companion in SNe Iax. We find that the He-star companion spins down after the impact due to the angular-momentum loss and expansion caused by the mass-stripping and shock heating during the interaction. This leads to the situation where the surface rotational speed of the surviving companion can drop to one-third of its pre-explosion value when it expands to a maximum radius a few years after the impact. Subsequently, the star shrinks and spins up again once the deposited energy is released. This spin-switching feature of the surviving He-star companions of SNe Iax may be useful for the identification of such objects in future observations.

preprint2022arXiv

Overview of the LAMOST survey in the first decade

The Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST), also known as the Guoshoujing Telescope, is a major national scientific facility for astronomical research located in Xinglong, China. Beginning with a pilot survey in 2011, LAMOST has been surveying the night sky for more than 10 years. The LAMOST survey covers various objects in the Universe, from normal stars to peculiar ones, from the Milky Way to other galaxies, and from stellar black holes and their companions to quasars that ignite ancient galaxies. Until the latest data release 8, the LAMOST survey has released spectra for more than 10 million stars, ~220,000 galaxies, and ~71,000 quasars. With this largest celestial spectra database ever constructed, LAMOST has helped astronomers to deepen their understanding of the Universe, especially for our Milky Way galaxy and the millions of stars within it. In this article, we briefly review the characteristics, observations, and scientific achievements of LAMOST. In particular, we show how astrophysical knowledge about the Milky Way has been improved by LAMOST data.

preprint2022arXiv

The Common Envelope Evolution Outcome -- A Case Study on Hot Subdwarf B Stars

Common envelope evolution (CEE) physics plays a fundamental role in the formation of binary systems, such as mergering stellar gravitational wave sources, pulsar binaries and type Ia supernovae. A precisely constrained CEE has become more important in the age of large surveys and gravitational wave detectors. We use an adiabatic mass loss model to explore how the total energy of the donor changes as a function of the remnant mass. This provides a more self-consistent way to calculate the binding energy of the donor. For comparison, we also calculate the binding energy through integrating the total energy from the core to the surface. The outcome of CEE is constrained by total energy conservation at the point at which both component's radii shrink back within their Roche lobes. We apply our results to 142 hot subdwarf binaries. For shorter orbital period sdBs, the binding energy is highly consistent. For longer orbital period sdBs in our samples, the binding energy can differ by up to a factor of 2. The CE efficiency parameter $β_\mathrm{CE}$ becomes smaller than $α_\mathrm{CE}$ for the final orbital period $\log_{10} P_{\mathrm{orb}}/\mathrm{d} > -0.5$. We also find the mass ratios $\log_{10} q$ and CE efficiency parameters $\log_{10} α_{\mathrm{CE}}$ and $\log_{10} β_{\mathrm{CE}}$ linearly correlate in sdBs, similarly to De Marco et al. (2010) for post-AGB binaries.

preprint2022arXiv

TYC~2990-127-1: an Algol-type SB2 binary system of subgiant and red giant with a probable ongoing mass-transfer

We present a study of the spectroscopic binary TYC 2990-127-1 from the LAMOST survey. We use full-spectrum fitting to derive radial velocities and spectral parameters. The high mass ratio indicates that the system underwent mass transfer in the past. We compute the orbital solution and find that it is a very close sub-giant/red giant pair on circular orbit, slightly inclined to the sky-plane. Fitting of the TESS photometrical data confirms this and suggests an inclination of i \sim 39.8 deg. The light curve and spectrum around H_alpha show signs of irregular variability, which supports ongoing mass transfer. The binary evolution simulations suggest that the binary may experience non-conservative mass transfer with accretion efficiency 0.3, and the binary will enter into common envelope phase in the subsequent evolution. The remnant product after the ejection of common envelope may be a detached double helium white dwarf (He WD) or a merger.

preprint2021arXiv

LAMOST Time-Domain Survey: First Results of four $K$2 plates

From Oct. 2019 to Apr. 2020, LAMOST performs a time-domain spectroscopic survey of four $K$2 plates with both low- and med-resolution observations. The low-resolution spectroscopic survey gains 282 exposures ($\approx$46.6 hours) over 25 nights, yielding a total of about 767,000 spectra, and the med-resolution survey takes 177 exposures ($\approx$49.1 hours) over 27 nights, collecting about 478,000 spectra. More than 70%/50% of low-resolution/med-resolution spectra have signal-to-noise ratio higher than 10. We determine stellar parameters (e.g., $T_{\rm eff}$, log$g$, [Fe/H]) and radial velocity (RV) with different methods, including LASP, DD-Payne, and SLAM. In general, these parameter estimations from different methods show good agreement, and the stellar parameter values are consistent with those of APOGEE. We use the $Gaia$ DR2 RV data to calculate a median RV zero point (RVZP) for each spectrograph exposure by exposure, and the RVZP-corrected RVs agree well with the APOGEE data. The stellar evolutionary and spectroscopic masses are estimated based on the stellar parameters, multi-band magnitudes, distances and extinction values. Finally, we construct a binary catalog including about 2700 candidates by analyzing their light curves, fitting the RV data, calculating the binarity parameters from med-resolution spectra, and cross-matching the spatially resolved binary catalog from $Gaia$ EDR3. The LAMOST TD survey is expected to get breakthrough in various scientific topics, such as binary system, stellar activity, and stellar pulsation, etc.

preprint2021arXiv

The Binarity of Early-type Stars from LAMOST Medium-resolution Spectroscopic Survey

Massive binaries play significant roles in many fields. Identification of massive stars, particularly massive binaries, is of great importance. In this paper, by adopting the technique of measuring the equivalent widths of several spectral lines, we identified 9,382 early-type stars from LAMOST medium-resolution survey and divided the sample into four groups, T1 ($\sim$O-B4), T2 ($\sim$B5), T3 ($\sim$B7), and T4 ($\sim$B8-A). The relative radial velocities $RV_{\rm rel}$ were calculated using the Maximum Likelihood Estimation. The stars with significant changes of $RV_{\rm rel}$ and at least larger than 15.57km s$^{-1}$ were identified as spectroscopic binaries. We found that the observed spectroscopic binary fractions for the four groups are $24.6\%\pm0.5\%$, $20.8\%\pm0.6\%$, $13.7\%\pm0.3\%$, and $7.4\%\pm0.3\%$, respectively. Assuming that orbital period ($P$) and mass ratio ($q$) have intrinsic distributions as $f(P) \propto P^π$ (1\textless$P$\textless1000 days) and $f(q) \propto q^κ$ (0.1\textless$q$\textless1), respectively, we conducted a series of Monte-Carlo simulations to correct observational biases for estimating the intrinsic multiplicity properties. The results show that the intrinsic binary fractions for the four groups are 68$\%\pm8\%$, 52$\%\pm3\%$, 44$\%\pm6\%$, and 44$\%\pm6\%$, respectively. The best estimated values for $π$ are -1$\pm0.1$, -1.1$\pm0.05$, -1.1$\pm0.1$, and -0.6$\pm0.05$, respectively. The $κ$ cannot be constrained for groups T1 and T2 and is -2.4$\pm0.3$ for group T3 and -1.6$\pm0.3$ for group T4. We confirmed the relationship of a decreasing trend in binary fractions towards late-type stars. No correlation between the spectral type and the orbital period distribution has been found yet, possibly due to the limitation of observational cadence.

preprint2020arXiv

Adiabatic mass loss in binary stars. III. From the base of the red giant branch to the tip of asymptotic giant branch

The distinguishing feature of the evolution of close binary stars is the role played by the mass exchange between the component stars. Whether the mass transfer is dynamically stable is one of the essential questions in binary evolution. 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 use the adiabatic mass loss model to systematically survey the thresholds for dynamical timescale mass transfer over the entire span of possible donor star evolutionary states. We also simulate mass loss process with isentropic envelopes, the specific entropy of which is fixed to be that at the base of the convective envelope, to artificially mimic the effect of such mass loss in superadiabatic surface convection regions, where the adiabatic approximation fails. We illustrate the general adiabatic response of 3.2 Msun donor stars at different evolutionary stages. We extend our study to a grid of donor stars with different masses (from 0.1 to 100 Msun with Z = 0.02) and at different evolutionary stages. We proceed to present our criteria for dynamically unstable mass transfer in both tabular and graphical forms. For red giant branch and asymptotic giant branch donors in systems with such mass ratios, they may have convective envelopes deep enough to evolve into common envelopes on a thermal timescale, if the donor star overfills its outer Lagrangian radius. Our results show that the red giant branch and asymptotic giant branch stars tend to be more stable than previously believed, and this may be helpful to explain the abundance of observed post-AGB binary stars with an orbital period of around 1000 days.

preprint2020arXiv

Gravitational Wave Radiation of Double Degenerates with Extremely low-mass WD companions

Double Degenerate systems (DDs) are supposed to be significant gravitational wave (GW) sources for future space-based gravitational-wave detectors, e.g., Laser Interferometer Space Antenna (LISA). Recently, one type of DDs with Extremely low-mass WD (ELM WD; $\leq 0.30\; M_\odot$) companions has been largely found in the ELM Survey. They have very short orbital periods and are therefore important sources for LISA detection. Besides, due to the thick envelope of ELM WDs compared with massive WDs (e.g. CO WDs), they are much easier to be found by the combination of electromagnetic (EM) and GW observations. In this paper, we first obtain the population of ELM WDs in DDs with considering the detailed evolutionary tracks of ELM WDs, and then analyse the GW radiation of these systems. We found that about $6\times10^3$ sources could be solely detected by LISA, including $\sim2\times10^3$ chirping sources, and $\sim13$ ($\sim107$) more sources are expected to be detected by both LISA and ELM Survey (Gaia).

preprint2020arXiv

KIC 10736223: An Algol-type eclipsing binary just undergone the rapid mass-transfer stage

This paper reports the discovery of an Algol system KIC 10736223 that just past the rapid mass transfer stage. From the light curve and radial-velocity modelling we find KIC 10736223 to be a detached Algol system with the less-massive secondary nearly filling its Roche lobe. Based on the short-cadence Kepler data, we analyzed intrinsic oscillations of the pulsator and identified six secured independent $δ$ Scuti-type pulsation modes ($f_{1}$, $f_3$, $f_{9}$, $f_{19}$, $f_{42}$, and $f_{48}$). We compute two grids of theoretical models to reproduce the $δ$ Scuti freqiencies, and find fitting results of mass-accreting models meet well with those of single-star evolutionary models. The fundamental parameters of the primary star yielded with asteroseismology are $M$ = $1.57^{+0.05}_{-0.09}$ $M_{\odot}$, $Z$ = 0.009 $\pm$ 0.001, $R$ = $1.484^{+0.016}_{-0.028}$ $R_{\odot}$, $\log g$ = $4.291^{+0.004}_{-0.009}$, $T_{\rm eff}$ = $7748^{+230}_{-378}$ K, $L$ = $7.136^{+1.014}_{-1.519}$ $L_{\odot}$. The asteroseismic parameters match well with the dynamical parameters derived from the binary model. Moreover, our asteroseismic results show that the pulsator is an almost unevolved star with an age between 9.46-11.65 Myr for single-star evolutionary models and 2.67-3.14 Myr for mass-accreting models. Thereofore, KIC 10736223 may be an Algol system that has just undergone the rapid mass-transfer process.

preprint2020arXiv

LAMOST Medium-Resolution Spectroscopic Survey (LAMOST-MRS): Scientific goals and survey plan

Since September 2018, LAMOST starts a new 5-year medium-resolution spectroscopic survey (MRS) using bright/gray nights. We present the scientific goals of LAMOST-MRS and propose a near optimistic strategy of the survey. A complete footprint is also provided. Not only the regular medium-resolution survey, but also a time-domain spectroscopic survey is being conducted since 2018 and will be end in 2023. According to the detailed survey plan, we expect that LAMOST-MRS can observe about 2 million stellar spectra with ~7500 and limiting magnitude of around G=15 mag. Moreover, it will also provide about 200 thousand stars with averagely 60-epoch observations and limiting magnitude of G~14 mag. These high quality spectra will give around 20 elemental abundances, rotational velocities, emission line profiles as well as precise radial velocity with uncertainty less than 1 km/s. With these data, we expect that LAMOST can effectively leverage sciences on stellar physics, e.g. exotic binary stars, detailed observation of many types of variable stars etc., planet host stars, emission nebulae, open clusters, young pre-main-sequence stars etc.

preprint2020arXiv

Simulating Kilonovae in the ΛCDM Universe

Kilonovae are optical flashes produced in the aftermath of neutron star-neutron star mergers (NNMs) or neutron star-black hole mergers (NBMs). In this work, we use the Millennium Simulation, combined with a semi-analytic galaxy formation model--GABE (Galaxy Assembly with Binary Evolution) which adopts binary stellar population synthesis models, to explore the cosmic event rate of kilonovae, and the properties of their host galaxies in a cosmological context. We find that model with supernova kick velocity of 0 km/s fits the observation best, in agreement with the exception of some formation channels of binary neutron star. This indicates that NNMs prefer to originate from binary systems with low kick velocities. With V$_{\rm kick}$=0 km/s, the cosmic event rate of NNMs and NBMs at z=0 are 283 Gpc$^{-3}$yr$^{-1}$ and 91 Gpc$^{-3}$yr$^{-1}$, respectively, marginally consistent with the constraint from LVC GWTC-1. For Milky Way-mass galaxies, we predict the NNM rate is $25.7^{+59.6}_{-7.1}$ Myr$^{-1}$, which is also in good agreement with the observed properties of binary neutron stars in the Milky Way. Taking all the NNMs into account in the history of Milky Way-mass galaxies, we find that the averaged r-process elements yield with A>79 in a NNM and NBM event should be 0.01 M$_{\odot}$ to be consistent with observation. We conclude that NGC 4993, the host galaxy of GW170817, is a typical host galaxy for NNMs. However, generally NNMs and NBMs tend to reside in young, blue, star-forming, late-type galaxies, with stellar mass and gaseous metallicity distribution peaking at $10^{10.65}$ M$_{\odot}$ and 8.72-8.85, respectively. By studying kilonovae host galaxies in the cosmological background, it is promising to constrain model details better when we have more events in the forthcoming future. (abridged)

preprint2020arXiv

The interaction of Type Iax supernova ejecta with a helium companion star

Type Iax supenovae (SNe Iax) are an important sub-luminous class of SNe Ia. However, their progenitors and explosion mechanism have not been understood yet. It has been suggested that SNe Iax may be produced from weak deflagration explosions of Chandrasekhar-mass white dwarfs (WDs) in binary systems with a helium (He) star donor. In such progenitor systems, the stripped He companion material caused by the ejecta-companion interaction is expected to present some He features in their late-time spectra. However, the detection of He lines in late-time spectra of SNe Iax has not yet been successful, which gives an upper limit on the amount of stripped He mass of < 0.002-0.1 solar masses. In this work, we study the interaction between SN Iax ejecta and a He star companion by performing three-dimensional hydrodynamical simulations with a weak pure deflagration explosion model. We find that about 0.004 solar masses of He material can be stripped off from the companion star by SN explosion, which is very close to (or lower than) the observational upper-limit on the total stripped He mass in SNe Iax. We, therefore, conclude that non-detection of He lines in late-time spectra of SNe Iax could be reasonably explained if they are indeed resulted from weak pure deflagration explosions of Chandrasekhar-mass WDs in progenitor systems with a He star donor.

preprint2020arXiv

The thermal equilibrium mass loss model and its applications in binary evolution

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 corresponding objects relative to what is expected from single stellar evolution. What is 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 Msun 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 Msun 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 Msun 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.

preprint2019arXiv

Bayesian discrimination of the panchromatic spectral energy distribution modelings of galaxies

Fitting the multi-wavelength spectral energy distributions (SEDs) of galaxies is a widely used technique to extract information about the physical properties of galaxies. However, a major difficulty lies in the numerous uncertainties regarding almost all ingredients of the SED modeling of galaxies. The Bayesian methods provide a consistent conceptual basis for dealing with the problem of inference with many uncertainties. While the Bayesian parameter estimation method have become quite popular in the field of SED fitting of galaxies, the Bayesian model comparison method, which is based on the same Bayes' rule, is still not widely used in this field. With the application of Bayesian model comparison method in a series of papers, we show that the results obtained with Bayesian model comparison are understandable in the context of stellar/galaxy physics. These results indicate that Bayesian model comparison is a reliable and very powerful method for the SED fitting of galaxies.

preprint2019arXiv

Hot subdwarf B stars with neutron star components II: Binary population synthesis

Context: Subdwarf B stars (sdBs) play a crucial role in stellar evolution, asteroseismology, and far-UV radiation of early-type galaxies, and have been intensively studied with observation and theory. It has theoretically been predicted that sdBs with neutron star (NS) companions exist in the Galaxy, but none have been discovered yet. This remains a puzzle in this field. In a previous study (hereafter Paper I), we have studied the formation channels of sdB+NS binaries from main-sequence (MS) stars plus NS binaries by establishing a model grid, but it is still unclear how these binaries consisting of MS stars and NS binaries came to be in the first place. Aims: We systematically study the formation of sdB+NS binaries from their original zero-age main-sequence progenitors. We bridge the gap left by our previous study in this way. We obtain the statistical population properties of sdB+NS binaries and provide some guidance for observational efforts. Methods: We first used Hurley's rapid binary evolution code BSE to evolve 10^7 primordial binaries to the point where the companions of NS+MS, NS+Hertzsprung gap (HG) star, and NS+Giant Branch (GB) star binaries have just filled their Roche lobes. Next, we injected these binaries into the model grid we developed in Paper I to obtain the properties of the sdB+NS populations. We adopted two prescriptions of NS natal kicks. Different values of common-envelope ejection efficiency were chosen to examine the effect of common-envelope evolution on the results. Conclusions: Most sdB+NS binaries are located in the Galactic disk with small RV semi-amplitudes. SdB+NS binaries with large RV semi-amplitudes are expected to be strong GWR sources, some of which could be detected by LISA in the future.

preprint2019arXiv

The formation of subdwarf A-type stars

Subdwarf A-type stars (sdAs) are objects that have hydrogen-rich spectra with surface gravity similar to that of hot subdwarf stars but effective temperature below the zero-age horizontal branch (ZAHB). They are considered to be metal-poor main sequence (MS) stars or extremely low-mass white dwarfs (ELM WDs). In this work, using the stellar evolution code, Modules for Experiments in Stellar Astrophysics (MESA), we investigate the sdAs formed both by the evolution of (pre-)ELM WDs in double degenerate systems (DDs) and metal-poor main sequence with the single evolution models. We find that both of the evolutionary tracks of ELM WDs and metal-poor MSs can explain the observation properties of sdAs. However, the proportion between these two populations are uncertain. In this work, we adopt the method of binary population synthesis of both ELM WDs in the disk and metal-poor MS in the halo to obtain the populations of ELM WDs and metal-poor MSs at different stellar population ages and calculate their proportions. We find that the proportion of metal-poor MSs to sdA population for a stellar population of 10 Gyr is $\sim$98.5 percent, which is consistent with the conclusion that most sdAs ($>$95 percent) are metal-poor MSs. And the proportion of ELM WDs (metal-poor MSs) to sdA population increases (decreases) from 0.1\% (99.9\%) to 20\% (80\%) with the stellar population ages from 5 to 13.7 Gyr.

preprint2016arXiv

Hot Subdwarf Stars Observed in LAMOST DR1 - Atmospheric parameters from single-lined spectra

We present a catalog of 166 spectroscopically identified hot subdwarf stars from LAMOST DR1, 44 of which show the characteristics of cool companions in their optical spectra. Atmospheric parameters of 122 non-composite spectra subdwarf stars were measured by fitting the profiles of hydrogen (H) and helium (He) lines with synthetic spectra from non-LTE model atmospheres. Most of the sdB stars scatter near the Extreme Horizontal Branch in the $T_{\rm eff}-\log{g}$ diagram and two well defined groups can be outlined. A clustering of He-enriched sdO stars appears near $T_{\rm eff}=45\,000$ K and $\log(g) = 5.8$. The sdB population separates into several nearly parallel sequences in the $T_{\rm eff}-{\rm He}$ abundance diagram with clumps corresponding to those in the $T_{\rm eff}-\log{g}$ diagram. Over $38\,000$ K (sdO) stars show abundance extremes, they are either He-rich or He-deficient and we observe only a few stars in the $ -1 < \log(y) < 0$ abundance range. With increasing temperature these extremes become less prominent and the He abundance approaches to $\log(y)\sim-0.5$. A unique property of our sample is that it covers a large range in apparent magnitudes and galactic latitudes, therefore it contains a mix of stars from different populations and galactic environments. Our results are consistent with the findings of Hirsch (2009) and we conclude that He-rich and He-deficient sdB stars ($\log(y) < 1$) probably origin from different populations. We also find that most sdO and sdB stars lie in a narrow strip in the luminosity and helium abundance plane, which suggests that these atmospheric parameters are correlated.

preprint2016arXiv

Modelling nova populations in galaxies

Theoretical modelling of the evolution of classical and recurrent novae plays an important role in studies of binary evolution, nucleosynthesis and accretion physics. However, from a theoretical perspective the observed statistical properties of novae remain poorly understood. In this paper, we have produced model populations of novae using a hybrid binary population synthesis approach for differing star formation histories (SFHs): a starburst case (elliptical-like galaxies), a constant star formation rate case (spiral-like galaxies) and a composite case (in line with the inferred SFH for M31). We found that the nova rate at 10\;Gyr in an elliptical-like galaxy is $\sim 10-20$ times smaller than a spiral-like galaxy with the same mass. The majority of novae in elliptical-like galaxies at the present epoch are characterized by low mass white dwarfs (WDs), long decay times, relatively faint absolute magnitudes and long recurrence periods. In contrast, the majority of novae in spiral-like galaxies at 10\;Gyr have massive WDs, short decay times, are relatively bright and have short recurrence periods. The mass loss time distribution for novae in our M31-like galaxy is in agreement with observational data for Andromeda. However, it is possible that we underestimate the number of bright novae in our model. This may arise in part due to the present uncertainties in the appropriate bolometric correction for novae.

preprint2016arXiv

The Detailed Science Case for the Maunakea Spectroscopic Explorer: the Composition and Dynamics of the Faint Universe

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

preprint2016arXiv

The X-ray/radio and UV luminosity expected from symbiotic systems as the progenitor of SNe Ia

We carried out a series of binary stellar evolution calculations, in which the effect of tidally enhanced wind on the evolution of WD + RG systems is incorporated. The WDs increase their mass to the Chandrasekhar mass limit, and then explode as SNe Ia. Based on the binary evolution results, we estimated the X-ray/radio (the excess UV) luminosity from the interactions between supernova ejecta and the CSM (the secondary) via some published standard models. We found that the X-ray flux may be high enough to be detected for a nearby SN Ia from a symbiotic system, while the radio flux is more likely to de detected when the companion is an asymptotic giant branch (AGB) star, and for a first giant branch (FGB) companion, the radio flux is generally lower than the detection limit. For two well observed SNe Ia, 2011fe and 2014J, almost all symbiotic systems are excluded by X-ray observations, but WD + FGB systems may not be ruled out by radio observations. The excess UV luminosity resulting from the collision of supernova ejecta with the RG secondary may be high enough to be detected if the secondary fills its Roche lobe at the moment of supernova explosion. The X-ray/radio emissions are more prevalent in SNe Ia from WD + AGB systems, although SNe Ia from such systems are rare. The UV luminosity from the collision of supernova ejecta to RG secondary is high enough, but only 1 in every few hundred SNe Ia manifests the signal from the collision.

preprint2015arXiv

A possible formation channel for blue hook stars in globular cluster

The formation mechanism for blue hook (BHk) stars in globular clusters (GCs) is still unclear. Following one of the possible scenario, named late hot flash scenario, we proposed that tidally enhanced stellar wind in binary evolution may provide the huge mass loss on the red giant branch (RGB) and produce BHk stars. Employing the detailed stellar evolution code, Modules for Experiments in Stellar Astrophysics (MESA), we investigated the contributions of tidally enhanced stellar wind as a possible formation channel for BHk stars in GCs. We evolved the primary stars with different initial orbital periods using the binary module in MESA (version 6208) from zero age main-sequence (ZAMS) to post horizontal branch (HB) stage, and obtained their evolution parameters which are compared with the observation. The results are consistent with observation in the color-magnitude diagram (CMD) and the logg-Teff plane for NGC 2808, which is an example GC hosting BHk stars. However, the helium abundance in the surface for our models is higher than the one obtained in BHk stars. This discrepancy between our models and observation is possibly due to the fact that gravitational settling and radiative levitation which are common processes in hot HB stars are not considered in the models as well as the fact that the flash mixing efficiency may be overestimated in the calculations. Our results suggested that tidally enhanced stellar wind in binary evolution is able to naturally provide the huge mass loss on the RGB needed for late hot flash scenario and it is a possible and reasonable formation channel for BHk stars in GCs.

preprint2015arXiv

A super-Eddington wind scenario for the progenitors of type Ia supernovae: binary population synthesis calculations

The super-Eddington wind scenario has been proposed as an alternative way for producing type Ia supernovae (SNe Ia). The super-Eddington wind can naturally prevent the carbon--oxygen white dwarfs (CO WDs) with high mass-accretion rates from becoming red-giant-like stars. Furthermore, it works in low-metallicity environments, which may explain SNe Ia observed at high redshifts. In this article, we systematically investigated the most prominent single-degenerate WD+MS channel based on the super-Eddington wind scenario. We combined the Eggleton stellar evolution code with a rapid binary population synthesis (BPS) approach to predict SN Ia birthrates for the WD+MS channel by adopting the super-Eddington wind scenario and detailed mass-accumulation efficiencies of H-shell flashes on the WDs. Our BPS calculations found that the estimated SN Ia birthrates for the WD+MS channel are ~0.009-0.315*10^{-3}{yr}^{-1} if we adopt the Eddington accretion rate as the critical accretion rate, which are much lower than that of the observations (<10% of the observed SN Ia birthrates). This indicates that the WD+MS channel only contributes a small proportion of all SNe Ia. The birthrates in this simulation are lower than previous studies, the main reason of which is that new mass-accumulation efficiencies of H-shell flashes are adopted. We also found that the critical mass-accretion rate has a significant influence on the birthrates of SNe Ia. Meanwhile, the results of our BPS calculations are sensitive to the values of the common-envelope ejection efficiency.

preprint2015arXiv

Adiabatic Mass Loss in Binary Stars. II. From Zero-Age Main Sequence to the Base of the Giant Branch

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 of mass 0.10 Msun to 100 Msun from the zero age main sequence to the base of the giant branch, or to central hydrogen exhaustion for lower main sequence stars. For intermediate- and high-mass stars, dynamical mass transfer is preceded by an extended phase of thermal time scale mass transfer as the star is stripped of most of its envelope mass. The critical mass ratio qad 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 q_ad declining with decreasing mass. Our calculated qad agree well with the behavior of timedependent models by Chen & Han (2003) of intermediate-mass stars initiating mass transfer in the Hertzsprung gap. Application of our results to cataclysmic variables, as systems which must be stable against rapid mass transfer, nicely circumscribes the range in qad as a function of orbital period in which they are found. These results are intended to advance the verisimilitude of population synthesis models of close binary evolution.

preprint2015arXiv

Population synthesis of accreting white dwarfs: II. X-ray and UV emission

Accreting white dwarfs (WDs) with non-degenerate companions are expected to emit in soft X-rays and the UV, if accreted H-rich material burns stably. They are an important component of the unresolved emission of elliptical galaxies, and their combined ionizing luminosity may significantly influence the optical line emission from warm ISM. In an earlier paper we modeled populations of accreting WDs, first generating WD with main-sequence, Hertzsprung gap and red giant companions with the population synthesis code \textsc{BSE}, and then following their evolution with a grid of evolutionary tracks computed with \textsc{MESA}. Now we use these results to estimate the soft X-ray (0.3-0.7keV), H- and He II-ionizing luminosities of nuclear burning WDs and the number of super-soft X-ray sources for galaxies with different star formation histories. For the starburst case, these quantities peak at $\sim 1$ Gyr and decline by $\sim 1-3$ orders of magnitude by the age of 10 Gyr. For stellar ages of $\sim$~10 Gyr, predictions of our model are consistent with soft X-ray luminosities observed by Chandra in nearby elliptical galaxies and He II 4686$Å/\rm{H}β$ line ratio measured in stacked SDSS spectra of retired galaxies, the latter characterising the strength and hardness of the UV radiation field. However, the soft X-ray luminosity and He~II~4686$Å/\rm{H}β$ ratio are significantly overpredicted for stellar ages of $\lesssim 4-8$ Gyr. We discuss various possibilities to resolve this discrepancy and tentatively conclude that it may be resolved by a modification of the typically used criteria of dynamically unstable mass loss for giant stars.

preprint2015arXiv

SNe Ia as a cosmological probe

Type Ia supernovae luminosities can be corrected to render them useful as standard candles able to probe the expansion history of the universe. This technique was successful applied to discover the present acceleration of the universe. As the number of SNe Ia observed at high redshift increases and analysis techniques are perfected, people aim to use this technique to probe the equation of state of the dark energy. Nevertheless, the nature of SNe Ia progenitors remains controversial and concerns persist about possible evolution effects that may be larger and harder to characterize than the more obvious statistical uncertainties.

preprint2015arXiv

Super-Eddington wind scenario for the progenitors of type Ia supernovae: Accreting He-rich matter onto white dwarfs

Supernovae of type Ia (SNe Ia) are believed to be thermonuclear explosions of carbon-oxygen white dwarfs (CO WDs). However, the mass accretion process onto CO WDs is still not completely understood. In this paper, we study the accretion of He-rich matter onto CO WDs and explore a scenario in which a strong wind forms on the surface of the WD if the total luminosity exceeds the Eddington limit. Using a stellar evolution code called modules for experiments in stellar astrophysics (MESA), we simulated the He accretion process onto CO WDs for WDs with masses of 0.6-1.35Msun and various accretion rates of 10^{-8}-10^{-5}Msun/yr. If the contribution of the total luminosity is included when determining the Eddington accretion rate, then a super-Eddington wind could be triggered at relatively lower accretion rates than those of previous studies based on steady-state models. The super-Eddington wind can prevent the WDs with high accretion rates from evolving into red-giant-like He stars. We found that the contributions from thermal energy of the WD are non-negligible, judging by our simulations, even though the nuclear burning energy is the dominating source of luminosity. We also provide the limits of the steady He-burning regime in which the WDs do not lose any accreted matter and increase their mass steadily, and calculated the mass retention efficiency during He layer flashes for various WD masses and accretion rates. These obtained results can be used in future binary population synthesis computations.

preprint2014arXiv

A pair of CO + He white dwarfs as the progenitor of 2005E-like supernovae?

Ca-rich transients (CRTs, 2005E-like supernovae) exhibit unusually strong Ca features during their nebular phase, and their distribution in their host galaxies indicates that they belong to a metal-poor old population. We carried out a series of binary population synthesis and presented four different channels of forming CO WD + He WD pairs. We selected the systems fulfilling the constraints of the old population and the birth rate from all the CO WD + He WD pairs by constraining the component mass of the WD pairs. For the four channels, the stable Roche lobe overflow (RLOF) could have a significant influence on the formation of the WD + WD pairs. Based on their position on the $M_{\rm CO}$-$M_{\rm He}$ plane, the mass-transfer between the components for most of the CO WD + He WD pairs is neither always unstable nor always stable. We found that it is necessary that the CO WDs are less massive than 0.6 $M_{\odot}$ and the He WDs are less massive than 0.25 $M_{\odot}$ if CO WD + He WD pairs fulfill the constraints from both the old population and the birth rate of CRTs. However, the He WD mass is lower than, while the total mass of the low-mass WD pairs is larger than{\bf ,} the ejecta mass of the CRTs derived from observations. Our results imply that the CO WDs participate in CRT explosions and at the same time, a bound remnant could be left after the CRT explosion if the low-mass WD pairs are the progenitors of CRTs. Therefore, whether or not the helium detonation on a low-mass CO WD may lead to the second detonation in the center of the CO WDs, and whether or not the binding remnant is left after the thermonuclear explosion should be examined carefully in the future.

preprint2014arXiv

BayeSED: A General Approach to Fitting the Spectral Energy Distribution of Galaxies

We present a newly developed version of BayeSED, a general Bayesian approach to the spectral energy distribution (SED) fitting of galaxies. The new BayeSED code has been systematically tested on a mock sample of galaxies. The comparison between estimated and inputted value of the parameters show that BayeSED can recover the physical parameters of galaxies reasonably well. We then applied BayeSED to interpret the SEDs of a large Ks-selected sample of galaxies in the COSMOS/UltraVISTA field with stellar population synthesis models. With the new BayeSED code, a Bayesian model comparison of stellar population synthesis models has been done for the first time. We found that the model by Bruzual & Charlot (2003), statistically speaking, has larger Bayesian evidence than the model by Maraston (2005) for the Ks-selected sample. Besides, while setting the stellar metallicity as a free parameter obviously increases the Bayesian evidence of both models, varying the IMF has a notable effect only on the Maraston (2005) model. Meanwhile, the physical parameters estimated with BayeSED are found to be generally consistent with those obtained with the popular grid-based FAST code, while the former exhibits more natural distributions. Based on the estimated physical parameters of galaxies in the sample, we qualitatively classified the galaxies in the sample into five populations that may represent galaxies at different evolution stages or in different environments. We conclude that BayeSED could be a reliable and powerful tool for investigating the formation and evolution of galaxies from the rich multi-wavelength observations currently available. A binary version of MPI parallelized BayeSED code is publicly available at https://bitbucket.org/hanyk/bayesed.

preprint2014arXiv

Binary interactions as a possible scenario for the formation of multiple stellar populations in globular clusters

Observations revealed the presence of multiple stellar populations in globular clusters (GCs) that exhibit wide abundance variations and multiple sequences in Hertzsprung-Russell (HR) diagram. We present a scenario for the formation of multiple stellar populations in GCs. In this scenario, initial GCs are single-generation clusters, and our model predicts that the abundance anomalous stars observed in GCs are the merged stars and the accretor stars produced by binary interactions, which are rapidly rotating stars at the moment of their formation and are more massive than normal single stars in the same evolutionary stage. We find that due to their own evolution, these rapidly rotating stars have different surface abundances, effective temperatures and luminosities from normal single stars in the same evolutionary stage. The stellar population with binaries can reproduce two important observational evidences of multiple stellar populations, the Na-O anticorrelation and the multiple sequences in HR diagram. This suggests that binary interactions may be a possible scenario for the formation of multiple stellar populations in GCs.

preprint2014arXiv

Effect of Binary Fraction on Horizontal Branch Morphology under Tidally Enhanced Stellar Wind

Tidally enhanced stellar wind may affect horizontal branch (HB) morphology in globular clusters (GCs) by enhancing the mass loss of primary star during binary evolution. Lei et al. (2013a, 2013b) studied the effect of this kind of wind on HB morphology in details, and their results indicated that binary is a possible secondparameter (2P) candidate in GCs. Binary fraction is an very important fact in the tidally-enhanced-stellar-wind model. In this paper, we studied the effect of binary fraction on HB morphology by removing the effects of metallicity and age. Five different binary fractions (i.e., 10%, 15%, 20%, 30% and 50%) are adopted in our model calculations. The synthetic HB morphologies with different binary fractions are obtained at different metallicities and ages. We found that, due to the great influence of metallicity and age, the effect of binary fraction on HB morphology may be masked by these two parameters. However, when the effects of metallicity and age are removed, the tendency that HB morphologies become bluer with increasing of binary fractions is clearly presented. Furthermore, we compared our results with the observation by Milone et al. (2012). Our results are consistent well with the observation at metalrich and metal-poor GCs. For the GCs with intermediate metallicity, when the effect of age on HB morphology is removed, a weak tendency that HB morphologies become bluer with increasing of binary fractions is presented in all regions of GCs, which is consistent with our results obtained in this metallicity range.

preprint2014arXiv

Next generation population synthesis of accreting white dwarfs: I. Hybrid calculations using BSE + MESA

Accreting, nuclear-burning white dwarfs have been deemed to be candidate progenitors of type Ia supernovae, and to account for supersoft X-ray sources, novae, etc. depending on their accretion rates. We have carried out a binary population synthesis study of their populations using two algorithms. In the first, we use the binary population synthesis code \textsf{BSE} as a baseline for the "rapid" approach commonly used in such studies. In the second, we employ a "hybrid" approach, in which we use \textsf{BSE} to generate a population of white dwarfs (WD) with non-degenerate companions on the verge of filling their Roche lobes. We then follow their mass transfer phase using the detailed stellar evolution code \textsf{MESA}. We investigate the evolution of the number of rapidly accreting white dwarfs (RAWDs) and stably nuclear-burning white dwarfs (SNBWDs), and estimate the type Ia supernovae (SNe Ia) rate produced by "single-degenerate" systems (SD). We find significant differences between the two algorithms in the predicted numbers of SNBWDs at early times, and also in the delay time distribution (DTD) of SD SNe Ia. Such differences in the treatment of mass transfer may partially account for differences in the SNe Ia rate and DTD found by different groups. Adopting 100\% efficiency for helium burning, the rate of SNe Ia produced by the SD-channel in a Milky-way-like galaxy in our calculations is $2.0\times10^{-4}\rm{yr}^{-1}$, more than an order of magnitude below the observationally inferred value. In agreement with previous studies, our calculated SD DTD is inconsistent with observations.

preprint2014arXiv

On the evolution of rotating accreting white dwarfs and type Ia supernovae

The potential importance of the angular momentum which is gained by accreting white dwarfs (WDs) has been increasingly recognized in the context of type Ia supernova (SN Ia) single-degenerate model. The expectation that the spin of the WD can delay the explosion should help the single-degenerate model to be consistent with the observed properties of most SNe Ia, in particular by avoiding hydrogen contamination. In this article, we attempt to study the most prominent single-degenerate supersoft (WD + MS) channel when the rotation of accreting WDs is considered. We present a detailed binary population synthesis study to examine the predicted population of SNe Ia for this channel. For our standard model, we find that 77% of these SNe Ia explode with WD masses which are low enough to be supported by solid-body rotation (1.378-1.5Msun); this is a substantially higher proportion than found by previous work. Only 2% have WD explosion masses >2.0Msun; these require the initial WD mass to be larger than 1.0 Msun. We further discuss the possible origin of the diversity of SNe Ia from the pre- and post- accretion properties of the WDs in this population. We also suggest that some SN Ia progenitors with substantial circumstellar hydrogen, including some apparent type IIn SNe, might be related to WDs which required support from differential rotation to avoid explosion, since these can still be accreting from hydrogen-rich donors with a relatively high mass-transfer rate at the time of the SN explosion.

preprint2014arXiv

Post-merger evolution of carbon-oxygen + helium white dwarf binaries and the origin of R Coronae Borealis and extreme helium stars

Orbital decay by gravitational-wave radiation will cause some close-binary white dwarfs (WDs) to merge within a Hubble time. The results from previous hydrodynamical WD-merger simulations have been used to guide calculations of the post-merger evolution of carbon-oxygen + helium (CO+He) WD binaries. Our models include the formation of a hot corona in addition to a Keplerian disk. We introduce a 'destroyeddisk' model to simulate the effect of direct disk ingestion into the expanding envelope. These calculations indicate significant lifetimes in the domain of the rare R Coronae Borealis (RCB) stars, before a fast evolution through the domain of the hotter extreme helium (EHe) stars. Surface chemistries of the resulting giants are in partial agreement with the observed abundances of RCB and EHe stars. The production of 3He, 18O and 19F are discussed. Evolutionary timescales combined with binary white-dwarf merger rates from binary-star population synthesis are consistent with present-day numbers of RCBs and EHes, provided that the majority come from relatively recent (< 2Gyr) star formation. However, most RCBs should be produced by CO-WD + low-mass HeWD mergers, with the He-WD having a mass in the range 0.20-0.35 Msolar. Whilst, previously, a high He-WD mass (> 0.40 Msolar) was required to match the carbon-rich abundances of RCB stars, the 'destroyed-disk' model yields a high-carbon product with He-WD mass > 0.30 Msolar, in better agreement with population synthesis results.

preprint2014arXiv

The hybrid CONe WD + He star scenario for the progenitors of type Ia supernovae

The hybrid CONe white dwarfs (WDs) have been suggested to be possible progenitors of type Ia supernovae (SNe Ia). In this article, we systematically studied the hybrid CONe WD + He star scenario for the progenitors of SNe Ia, in which a hybrid CONe WD increases its mass to the Chandrasekhar mass limit by accreting He-rich material from a non-degenerate He star. According to a series of detailed binary population synthesis simulations, we obtained the SN Ia birthrates and delay times for this scenario. The SN Ia birthrates for this scenario are ~0.033-0.539*10^(-3)yr^(-1), which roughly accounts for 1-18% of all SNe Ia. The estimated delay times are ~28Myr-178Myr, which are the youngest SNe Ia predicted by any progenitor model so far. We suggest that SNe Ia from this scenario may provide an alternative explanation of type Iax SNe. We also presented some properties of the donors at the point when the WDs reach the Chandrasekhar mass. These properties may be a good starting point for investigating the surviving companions of SNe Ia, and for constraining the progenitor scenario studied in this work.

preprint2013arXiv

A Super-Eddington Wind Scenario for the Progenitors of Type Ia Supernovae

The accretion of hydrogen-rich material onto carbon-oxygen white dwarfs (CO WDs) is crucial for understanding type Ia supernova (SN Ia) from the single-degenerate model, but this process has not been well understood due to the numerical difficulties in treating H and He flashes during the accretion. For the CO WD masses from 0.5 to $1.378\,{M}_\odot$ and accretion rates in the range from $10^{-8}$ to $10^{-5}\,{M}_\odot\,\mbox{yr}^{-1}$, we simulated the accretion of solar-composition material onto CO WDs using the state-of-the-art stellar evolution code of {\sc MESA}. For comparison with the steady-state models (e.g \citet{nskh07}), we firstly ignored the contribution from nuclear burning to the luminosity when determining the Eddington accretion rate and found that the properties of H burning in our accreting CO WD models are similar to those from the steady-state models, except that the critical accretion rates at which the WDs turn into red giants or H-shell flashes occur on their surfaces are slightly higher than those from the steady-state models. However, the super-Eddington wind is triggered at much lower accretion rates, than previously thought, when the contribution of nuclear burning to the total luminosity is included. This super-Eddington wind naturally prevents the CO WDs with high accretion rates from becoming red giants, thus presenting an alternative to the optically thick wind proposed by \cite{hkn96}. Furthermore, the super-Eddington wind works in low-metallicity environments, which may explain SNe Ia observed at high redshifts.

preprint2013arXiv

Double-detonation explosions as progenitors of type Iax supernovae

It has recently been proposed that one sub-class of type Ia supernovae (SNe Ia) is sufficiently both distinct and common to be classified separately from the bulk of SNe Ia, with a suggested class name of "type Iax supernovae" (SNe Iax), after SN 2002cx. We show that the population properties of this class can be understood if the events originate from helium double-detonation sub-Chandrasekhar mass explosions, in which a carbon--oxygen white dwarf (CO WD) accumulates a helium layer from a non-degenerate helium star. We have incorporated detailed binary evolution calculations for the progenitor systems into a binary population synthesis model to obtain birthrates and delay times for such events. The predicted Galactic event rate is $\sim$$0.6-1.8\times 10^{-3}\,{\rm yr}^{-1}$, in good agreement with the measured rates of SNe Iax. In addition, predicted delay times are $\sim$70\,Myr$-$710\,Myr, consistent with the fact that SNe Iax have so far only been discovered in late-type galaxies. Based on the CO WD mass at explosion and previous detonation models, we also estimate the distribution of resulting SN brightness ($-13 \gtrsim M_{\rm bol} \gtrsim -19$ mag), which can reproduce the empirical diversity of SNe Iax. We speculate on why binaries with non-degenerate donor stars might lead to SNe Iax if similar systems with degenerate donors do not and suggest that the higher mass of the helium layer necessary for ignition at the lower accretion rates typically delivered from non-degenerate donors might provide the explanation.

preprint2013arXiv

Fitting Formulae for the Effects of Binary Interactions on Lick Indices and Colours of Stellar Populations

More than about 50% stars are in binaries, but most stellar population studies take single star stellar population (ssSSP) models, which do not take binary interactions into account. In fact, the integrated peculiarities of ssSSPs are various from those of stellar populations with binary interactions (bsSSPs). Therefore, it is necessary to investigate the effects of binary interactions on the Lick indices and colours of populations detailedly. We show some formulae for calculating the difference between the Lick indices and colours of bsSSPs, and those of ssSSPs. Twenty-five Lick indices and 12 colours are studied in the work. The results can be conveniently used for calculating the effects of binary interactions on stellar population studies and for adding the effects of binary interactions into present ssSSP models. The electrical data and fortran procedures of the paper can be obtained on request to the authors.

preprint2013arXiv

Formation of black widows and redbacks -- two distinct populations of eclipsing binary millisecond pulsars

Eclipsing binary millisecond pulsars (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 millisecond pulsars. 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-1.0\;{\rm days}$ their companion masses differ by an order of magnitude. In this paper, we investigate the formation of these systems via 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 millisecond pulsars 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) is responsible for the strong bimodality of these two populations. Finally, we conclude that redback systems do not evolve into black widow systems with time.

preprint2013arXiv

On the minimum mass ratio of W UMa binaries

Using Eggleton's stellar evolution code, we study the minimum mass ratio ($q_{\rm min}$) of W Ursae Majoris (W UMa) binaries that have different primary masses. It is found that the minimum mass ratio of W UMa binaries decreases with increasing mass of the primary if the primary's mass is less than about 1.3$M_{\rm \odot}$, and above this mass the ratio is roughly constant. By comparing the theoretical minimum mass ratio with the observational data, it is found that the existence of low-$q$ systems can be explained by the different structure of the primaries with different masses. This suggests that the dimensionless gyration radius ($k_1^2$) and thus the structure of the primary is very important in determining the minimum mass ratio. In addition, we investigate the mass loss during the merging process of W UMa systems and calculate the rotation velocities of the single stars formed by the merger of W UMa binaries due to tidal instability. It is found that in the case of the conservation of mass and angular momentum, the merged single stars rotate with a equatorial velocity of about $588\sim819$ km s$^{-1}$, which is much larger than their break-up velocities ($v_{\rm b}$). This suggests that the merged stars should extend to a very large radius (3.7$\sim$5.3 times the radii of the primaries) or W UMa systems would lose a large amount of mass (21$\sim$33 per cent of the total mass) during the merging process. If the effect of magnetic braking is considered, the mass loss decreases to be 12$\sim$18 per cent of their total masses. This implies that the significant angular momentum and mass might be lost from W UMa systems in the course of the merging process, and this kind of mass and angular momentum loss might be driven by the release of orbital energy of the secondaries, which is similar to common-envelope evolution.

preprint2013arXiv

Producing type Iax supernovae from a specific class of helium-ignited WD explosions?

It has recently been proposed that one sub-class of type Ia supernovae (SNe Ia) is sufficiently both distinct and common to be classified separately from the bulk of SNe Ia, with a suggested class name of "type Iax supernovae" (SNe Iax), after SN 2002cx. However, their progenitors are still uncertain. We study whether the population properties of this class might be understood if the events originate from a subset of sub-Chandrasekhar mass explosions. In this potential progenitor population, a carbon--oxygen white dwarf (CO WD) accumulates a helium layer from a non-degenerate helium star; ignition of that helium layer then leads to ignition of the CO WD. We incorporated detailed binary evolution calculations for the progenitor systems into a binary population synthesis model to obtain rates and delay times for such events. The predicted Galactic event rate of these explosions is ~1.5\times10^{-3}{yr}^{-1} according to our standard model, in good agreement with the measured rates of SNe Iax. In addition, predicted delay times are ~70Myr-800Myr, consistent with the fact that most of SNe Iax have been discovered in late-type galaxies. If the explosions are assumed to be double-detonations -- following current model expectations -- then based on the CO WD masses at explosion we also estimate the distribution of resulting SN brightness (-13 \gtrsim M_{bol} \gtrsim -19mag), which can reproduce the empirical diversity of SNe Iax. We speculate on why binaries with non-degenerate donor stars might lead to SNe Iax if similar systems with degenerate donors do not. We suggest that the high mass of the helium layer necessary for ignition at the lower accretion rates typically delivered from non-degenerate donors might be necessary to produce SN 2002cx-like characteristics, perhaps even by changing the nature of the CO ignition.

preprint2013arXiv

The binary merger channel for the progenitor of the fastest rotating O-type star VFTS102

VFTS102 has a projected rotational velocity ($>$500\,km s$^{-1}$) and would appear to be the fastest rotating O-type star. We show that its high rotational velocity could be understood within the framework of the binary merger. In the binary merger channel, the progenitor binary of VFTS102 would evolve into contact while two components are still on the main sequence, and then merge into a rapidly rotating single star. Employing Eggleton's stellar evolution code, we performed binary stellar evolution calculations and mapped out the initial parameters of the progenitor of VFTS102 in the orbital period-mass ratio ($P-q$) plane. We found that the progenitor binary of VFTS102 with initial mass ratio $q_0 \la 0.7$ should have an initial orbital period shorter than $3.76-4.25$\,days, while above this mass ratio it should have an initial orbital period shorter than $1.44-1.55$\,days. The progenitor of VFTS102 would evolve into contact during the rapid mass transfer phase or during the subsequent slow mass transfer phase, and might ultimately merge into a rapidly rotating massive star. In addition, we performed Monte Carlo simulations to investigate the binary merger channel. We estimated the fraction of binaries that would merge into single stars and the fraction of single stars that might be produced from the binary merger channel. It is found that about 8.7% of binaries would evolve into contact and merge into rapidly rotating single stars, and about 17.1% of single stars might be produced from the binary merger channel and should have similar properties to VFTS102. This suggests that the binary merger channel might be one of the main channels for the formation of rapidly rotating massive stars like VFTS102.

preprint2013arXiv

The detached-binary channel for the formation of contact binaries

The detached-binary channel is an important channel for the formation of contact binaries, according to which a detached binary might evolve into contact by evolutionary expansion of the components, or angular momentum loss through the effect of magnetic braking (MB). We have carried out a detailed binary population synthesis (BPS) study of this channel, and obtained the parameter regions for detached binaries to evolve into contact. Combining the observations from the $\emph{Kepler}$ satellite with our results, we found that the ratio of the birth rate of the progenitors of contact binaries to that of contact binaries is greater than about 1.2. This suggests that for the detached-binary channel, the progenitors can be sufficient to produce observed contact binaries. In this channel, we find that the distribution of orbital period for contact binaries has a peak at about 0.25 days and a tail extending to longer periods, and the formation timescale of contact binaries has a large range from $\sim1$\,Myr to 15\,Gyr. These results show that the detached-binary channel could explain satisfactorily the main observational characteristics of contact binaries, in particular the distribution of orbital period shown by the $\emph{Kepler}$ observations and the existence of very young contact binaries.

preprint2013arXiv

The impact of Type Ia supernova explosions on helium companions in the Chandrasekhar-mass explosion scenario

In the version of the SD scenario of SNe Ia studied here, a CO WD explodes close to the Chandrasekhar limit after accreting material from a non-degenerate He companion. In the present study, we employ the Stellar GADGET code to perform 3D hydrodynamical simulations of the interaction of the SN Ia ejecta with the He companion taking into account its orbital motion and spin. It is found that only 2%--5% of the initial companion mass are stripped off from the outer layers of He companions due to the SN impact. The dependence of the unbound mass (or the kick velocity) on the orbital separation can be fitted in good approximation by a power law for a given companion model. After the SN impact, the outer layers of a He donor star are significantly enriched with heavy elements from the low-expansion-velocity tail of SN Ia ejecta. The total mass of accumulated SN-ejecta material on the companion surface reaches about > 10e-3 M_sun for different companion models. This enrichment with heavy elements provides a potential way to observationally identify the surviving companion star in SN remnants. Finally, by artificially adjusting the explosion energy of the W7 explosion model, we find that the total accumulation of SN ejecta on the companion surface is also dependent on the explosion energy with a power law relation in good approximation.

preprint2013arXiv

The orbital periods of subdwarf B binaries produced by the first stable Roche overflow channel

Long-orbital-period subdwarf B (sdB) stars with main-sequence companions are believed to be the product of stable Roche Lobe overflow (RLOF), a scenario challenged by recent observations. Here we represent the results of a systematic study of the orbital-period distribution of sdB binaries in this channel using detailed binary evolution calculations. We show that the observed orbital-period distribution of long-period sdB binaries can be well explained by this scenario. Furthermore, we find that, if the progenitors of the sdB stars have initial masses below the helium flash mass, the sdB binaries produced from stable RLOF follow a unique mass -- orbital period relation for a given metallicity $Z$; increasing the orbital period from $\sim 400$ to $\sim 1100$\,d corresponds to increasing the mass of the sdB star from $\sim 0.40$ to $\sim 0.49\,M_\odot$ for $Z=0.02$. We suggest that the longest sdB binaries (with orbital period $> 1100$\,d) could be the result of atmospheric RLOF. The mass -- orbital period relation can be tested observationally if the mass of the sdB star can be determined precisely, e.g.\ from asteroseismology. Using this relation, we revise the orbital period distribution of sdB binaries produced by the first stable RLOF channel for the best fitting model of Han et al (2003), and show that the orbital period has a peak around 830\,d.

preprint2013arXiv

Uncertainties in the calibrations of star formation rate

The calibrations of star formation rate (SFR) are prone to be affected by many factors, such as metallicity, initial mass function (IMF), evolutionary population synthesis (EPS) models and so on. In this paper we will discuss the effects of binary interactions, metallicity, EPS models and IMF on several widely used SFR calibrations based on the EPS models of Yunnan with and without binary interactions, BC03, SB99, PEGASE and POPSTAR. The inclusion of binary interactions makes these SFR conversion coefficients smaller (less than 0.2dex), and these differences increase with metallicity. The differences in the calibration coefficient between SFR and the luminosity of $\rm Hα$ recombination line (C$_{\rm Hα}$) and that between SFR and the ultraviolet (UV) fluxes at 1500 and 2800\,$\rm Å$ (C$_{i, {\rm UV}}$), caused by IMF, are independent of metallicity (0.03-0.33\,dex) except $Δ$C$_{\rm Hα, IMF}$ when using the POPSTAR and $Δ$C$_{i, {\rm UV, IMF}}$ when using the PEGASE models. Moreover, we find that $L_{\rm 2800}$ is not suitable to the linear calibration of SFR at low metallicities. At last, we compare the effects of these several factors on the SFR calibrations considered in this paper. The effects of metallicity/IMF and EPS models on the C$_{\rm Hα}$ and C$_{\rm FIR}$ (the conversion coefficient between SFR and the far-infrared flux) are the largest among these factors, respectively. For the calibration between SFR and C$_{i, {\rm UV}}$, the effects of these several factors are comparable.

preprint2012arXiv

An Algorithm for Preferential Selection of Spectroscopic Targets in LEGUE

We describe a general target selection algorithm that is applicable to any survey in which the number of available candidates is much larger than the number of objects to be observed. This routine aims to achieve a balance between a smoothly-varying, well-understood selection function and the desire to preferentially select certain types of targets. Some target-selection examples are shown that illustrate different possibilities of emphasis functions. Although it is generally applicable, the algorithm was developed specifically for the LAMOST Experiment for Galactic Understanding and Exploration (LEGUE) survey that will be carried out using the Chinese Guo Shou Jing Telescope. In particular, this algorithm was designed for the portion of LEGUE targeting the Galactic halo, in which we attempt to balance a variety of science goals that require stars at fainter magnitudes than can be completely sampled by LAMOST. This algorithm has been implemented for the halo portion of the LAMOST pilot survey, which began in October 2011.

preprint2012arXiv

Decoding spectral energy distributions of dust-obscured starburst-AGN

We present BayeSED, a general purpose tool for doing Bayesian analysis of SEDs by using whatever pre-existing model SED libraries or their linear combinations. The artificial neural networks (ANNs), principal component analysis (PCA) and multimodal nested sampling (MultiNest) techniques are employed to allow a 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 evidences obtained for a pure Starburst, a pure AGN, and a linear combination of Starburst+AGN models show that the Starburst+AGN model have the highest evidence for all galaxies in this sample. The Bayesian evidences for the three models and the estimated contributions of starburst and AGN to infrared luminosity show that HLIRGs can be classified into two groups: one dominated by starburst and the other dominated by AGN. Other parameters and corresponding uncertainties about starburst and AGN are also estimated by using the model with the highest Bayesian evidence. We found that the starburst region of the HLIRGs dominated by starburst tends to be more compact and has a higher fraction of OB star than that of HLIRGs dominated by AGN. Meanwhile, the AGN torus of the HLIRGs dominated by AGN tend to be more dusty than that of HLIRGs dominated by starburst. 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 from decoding their SEDs.

preprint2012arXiv

Evolution of luminosity function and obscuration of AGN: connecting X-ray and Infrared

We present a detailed comparison between the 2-10 keV hard X-ray and infrared (IR) luminosity function (LF) of active galactic nuclei (AGN). The composite X-ray to IR spectral energy distributions (SEDs) of AGN used for connecting the hard X-ray LF (HXLF) and IR LF (IRLF) are modeled with a simple but well tested torus model based on the radiative transfer and photoionization code CLOUDY. Four observational determinations of the evolution of 2-10 keV HXLF and six evolution models of the obscured type-2 AGN fraction ($f_2$) have been considered. The 8.0 and 15 \micron LFs for the total, unobscured type-1 and obscured type-2 AGN are predicted from the HXLFs, and then compared with the measurements currently available. We find that the IRLFs predicted from HXLFs tend to underestimate the number of the most IR-luminous AGN. This is independent of the choices of HXLF and $f_2$, and even more obvious for the HXLFs recently measured. We show that the discrepancy between the HXLFs and IRLFs can be largely resolved when the anticorrelation between the UV to X-ray slope $α_{\mathrm{ox}}$ and UV luminosity $L_{\rm UV}$ is appropriately considered. We also discuss other possible explanations for the discrepancy, such as the missing population of Compton-thick AGN and possible contribution of star-formation in the host to the mid-IR. Meanwhile, we find that the HXLFs and IRLFs of AGN can be more consistent with each other if the obscuration mechanisms of quasars and Seyferts are assumed to be different, corresponding to their different triggering and fueling mechanisms. More accurate measurements of the IRLFs of AGN, especially that determined at smaller redshift bins and more accurately separated to that for type-1 and type-2, are very helpful for clarifying these interesting issues.

preprint2012arXiv

LAMOST Experiment for Galactic Understanding and Exploration (LEGUE) The survey science plan

We describe the current plans for a spectroscopic survey of millions of stars in the Milky Way galaxy using the Guo Shou Jing Telescope (GSJT, formerly the Large Area Multi-Object Spectroscopic Telescope - LAMOST). The survey will obtain spectra for 2.5 million stars brighter than $r<19$ during dark/grey time, and 5 million stars brighter than $r<17$ or $J<16$ on nights that are moonlit or have low transparency. The survey will begin in fall of 2012, and will run for at least four years. The telescope design constrains the optimal declination range for observations to $10^\circ<δ<50^\circ$, and site conditions lead to an emphasis on stars in the direction of the Galactic anticenter. The survey is divided into three parts with different target selection strategies: disk, anticenter, and spheroid. The resulting dataset will be used to study the merger history of the Milky Way, the substructure and evolution of the disks, the nature of the first generation of stars through identification of the lowest metallicity stars, and star formation through study of open clusters and the OB associations. Detailed design of the LEGUE survey will be completed after a review of the results of the pilot survey in summer 2012.

preprint2012arXiv

Progenitors of type Ia supernovae

Type Ia supernovae (SNe Ia) play an important role in astrophysics and are crucial for the studies of stellar evolution, galaxy evolution and cosmology. They are generally thought to be thermonuclear explosions of accreting carbon-oxygen white dwarfs (CO WDs) in close binaries, however, the nature of the mass donor star is still unclear. In this article, we review various progenitor models proposed in the past years and summarize many observational results that can be used to put constraints on the nature of their progenitors. We also discuss the origin of SN Ia diversity and the impacts of SN Ia progenitors on some fields. The currently favourable progenitor model is the single-degenerate (SD) model, in which the WD accretes material from a non-degenerate companion star. This model may explain the similarities of most SNe Ia. It has long been argued that the double-degenerate (DD) model, which involves the merger of two CO WDs, may lead to an accretion-induced collapse rather than a thermonuclear explosion. However, recent observations of a few SNe Ia seem to support the DD model, and this model can produce normal SN Ia explosion under certain conditions. Additionally, the sub-luminous SNe Ia may be explained by the sub-Chandrasekhar mass model. At present, it seems likely that more than one progenitor model, including some variants of the SD and DD models, may be required to explain the observed diversity of SNe Ia.

preprint2012arXiv

Some constraints on the lower mass limit for double-degenerate progenitors of Type Ia supernovae

Recent theoretical and observational studies both argue that the merging of double carbon-oxygen white dwarfs (WDs) is responsible for at least some Type Ia supernovae (SNe Ia). Previous (standard) studies of the anticipated SN birthrate from this channel have assumed that the merger process is conservative and that the primary criterion for explosion is that the merged mass exceeds the Chandrasekhar mass. Han & Webbink (1999) demonstrated that mass transfer and merger in close double WDs will in many cases be non-conservative. Pakmor et al. (2011) further suggested that the merger process should be violent in order to initiate an explosion. We have therefore investigated how the SN Ia birthrate from the double-degenerate (DD) channel is affected by these constraints. Using the binary-star population-synthesis method, we have calculated the DD SN Ia birthrate under conservative and non-conservative approximations, and including lower mass and mass-ratio limits indicated by recent smoothed-particle-hydrodynamic calculations. The predicted DD SN Ia rate is significantly reduced by all of these constraints. With dynamical mass loss alone (violent merger) the birthrate is reduced to 56% of the conservative rate. Requiring that the mass ratio $q>2/3$ further reduces the birthrate to 18% that of the standard assumption. An upper limit of 0.0061 SNuM, or a Galactic rate of $4.6 \times 10^{-4}{\rm yr}^{-1}$, might be realistic.

preprint2012arXiv

The Evolution and Star Formation History of M33

We construct a parameterized model to explore the main properties of the star formation history of M33. We assume that the disk originates and grows by the primordial gas infall and adopt the simple form of gas accretion rate with one free parameter, the infall time-scale. We also include the contribution of gas outflow process. A major update of the model is that we adopt a molecular hydrogen correlated star formation law and calculate the evolution of the atomic and molecular gas separately. Comparisons between the model predictions and the observational data show that the model predictions are very sensitive to the adopted infall time-scale, while the gas outflow process mainly influences the metallicity profile. The model adopting a moderate outflow rate and an inside-out formation scenario can be in good agreement with most of observed constraints of M33 disk. We also compare the model predictions based on the molecular hydrogen correlated star formation law and that based on the Kennicutt star formation law. Our results imply that the molecular hydrogen correlated star formation law should be preferred to describe the evolution of the M33 disk, especially the radial distributions of both the cold gas and the stellar population.

preprint2012arXiv

The LEGUE High Latitude Bright Survey Design for the LAMOST Pilot Survey

We describe the footprint and input catalog for bright nights in the LAMOST Pilot Survey, which began in October 2011. Targets are selected from two stripes in the north and south Galactic Cap regions, centered at $α$= 29$^\circ$, with 10$^\circ$ width in declination, covering right ascension of 135$^\circ-290^\circ$ and -30$^\circ$ to 30$^\circ$ respectively. We selected spectroscopic targets from a combination of the SDSS and 2MASS point source catalogs. The catalog of stars defining the field centers (as required by the Shack-Hartmann wavefront sensor at the center of the LAMOST field) consists of all V < 8m stars from the Hipparcos catalog. We employ a statistical selection algorithm that assigns priorities to targets based on their positions in multidimensional color/magnitude space. This scheme overemphasizes rare objects and de-emphasizes more populated regions of magnitude and color phase space, while ensuring a smooth, well-understood selection function. A demonstration of plate design is presented based on the Shack-Hartmann star catalog and an input catalog that was generated by our target selection routines.

preprint2012arXiv

The LEGUE Input Catalogue for Dark Night Observing in the LAMOST Pilot Survey

We outline the design of the dark nights portion of the LAMOST Pilot Survey, which began observations in October 2011. In particular, we focus on Milky Way stellar candidates that are targeted for the LEGUE (LAMOST Experiment for Galactic Understanding and Exploration) survey. We discuss the regions of sky in which spectroscopic candidates were selected, and the motivations for selecting each of these sky areas. Some limitations due to the unique design of the telescope are discussed, including the requirement that a bright (V < 8) star be placed at the center of each plate for wavefront sensing and active optics corrections. The target selection categories and scientific goals motivating them are briefly discussed, followed by a detailed overview of how these selection functions were realized. We illustrate the difference between the overall input catalog - Sloan Digital Sky Survey (SDSS) photometry - and the final targets selected for LAMOST observation.

preprint2012arXiv

The selection of LEGUE disk targets for LAMOST's pilot survey

We describe the target selection algorithm for the low latitude disk portion of the LAMOST Pilot Survey, which aims to test systems in preparation for the LAMOST spectroscopic survey. We use the PPMXL (Roeser et al. 2010) astrometric catalog, which provides positions, proper motions, B/R/I magnitudes (mostly) from USNO-B (Monet et al. 2003) and J/H/Ks from The Two Micron All Sky Survey (2MASS, see Skrutskie et al. 2006) as well. We chose 8 plates along the Galactic plane, in the region $0^\circ<α<67^\circ$ and $42^\circ<δ<59^\circ$, that cover 22 known open clusters with a range of ages. Adjacent plates may have small overlapping. Each plate covers an area $2.5^\circ$ in radius,with central star (for Shack-Hartmann guider) brighter than $\sim8^{\rm th}$ magnitude. For each plate, we create an input catalog in the magnitude range $11.3<Imag<16.3$ and $Bmag$ available from PPMXL. The stars are selected to satisfy the requirements of the fiber positioning system and have a uniform distribution in the $I$ vs. $B-I$ color-magnitude diagram. Our final input catalog consists of 12,000 objects on each of 8 plates that are observable during the winter observing season in Xinglong Station of the National Astronomical Observatory of China.

preprint2012arXiv

The site conditions of the Guo Shou Jing Telescope

The weather at Xinglong Observing Station, where the Guo Shou Jing Telescope (GSJT) is located, is strongly affected by the monsoon climate in north-east China. The LAMOST survey strategy is constrained by these weather patterns. In this paper, we present a statistics on observing hours from 2004 to 2007, and the sky brightness, seeing, and sky transparency from 1995 to 2011 at the site. We investigate effects of the site conditions on the survey plan. Operable hours each month shows strong correlation with season: on average there are 8 operable hours per night available in December, but only 1-2 hours in July and August. The seeing and the sky transparency also vary with seasons. Although the seeing is worse in windy winters, and the atmospheric extinction is worse in the spring and summer, the site is adequate for the proposed scientific program of LAMOST survey. With a Monte Carlo simulation using historical data on the site condition, we find that the available observation hours constrain the survey footprint from 22h to 16h in right ascension; the sky brightness allows LAMOST to obtain the limit magnitude of V = 19.5mag with S/N = 10.

preprint2011arXiv

Helium star donor channel to type Ia supernovae and their surviving companion stars

Employing Eggleton's stellar evolution code with the optically thick wind assumption, we systematically studied the He star donor channel of type Ia supernovae (SNe Ia), in which a carbon-oxygen white dwarf (WD) accretes material from a He main-sequence star or a He subgiant to increase its mass to the Chandrasekhar mass. We mapped out the initial parameters for producing SNe Ia in the orbital period--secondary mass plane for various WD masses from this channel. Based on a detailed binary population synthesis approach, we find that this channel can produce SNe Ia with short delay times (\sim100Myr) implied by recent observations. We obtained many properties of the surviving companions of this channel after SN explosion, which can be verified by future observations. We also find that the surviving companions from the SN explosion scenario have a high spatial velocity (>400km/s), which could be an alternative origin for hypervelocity stars (HVSs), especially for HVSs such as US 708.

preprint2011arXiv

Testing three derivative methods of stellar population synthesis models

The number of methods used to study the properties of galaxies is increased, and testing these methods is very important. Galactic globular clusters (GCs) provide an excellent medium for such test, because they can be considered as simple stellar populations. We present ages and metallicities for 40 Galactic GCs as determined from three publicly available techniques, including colour, Lick-index and spectrum-fitting methods, based on Bruzual & Charlot evolutionary population synthesis (EPS) models. By comparing with the ages obtained from colour-magnitude diagrams (CMDs) and metallicities obtained from spectra of stars, we are able to estimate the ability of 'these methods on determination of GCs parameters, which is absolutely necessary. As a result, we find that: (i) for the metallicity, our derived metallicities agree with those derived from the spectra of stars, Lick-index method is suitable to study metallicity for the stellar population systems in the range of -1.5=<[Fe/H]=<-0.7 and spectrum- fitting method is suitable to study metallicity for the stellar population systems in the range of -2.3=<[Fe/H]=<-1.5; (ii) for the age, these three methods have difficulties in age determination, our derived ages are smaller (about 2.0 Gyr, on average) than the results of CMDs for all these three methods. We use Vazdekis and Maraston models to analyze whether our results are dependent on EPS models, and find that the tendency of these two models is the same as that of Bruzual & Charlot models. Our results are independent of the EPS models. In addition, our test is based on the old GCs and our conclusions may hold for old stellar population systems. The whole abstract can be found in my PDF version.

preprint2011arXiv

The short-period limit of contact binaries

The stability of mass transfer is important in the formation of contact binaries from detached binaries when the primaries of the initially detached binaries fill their Roche lobes. Using Eggleton's stellar evolution code, we investigate the formation and the short-period limit of contact binaries by considering the effect of the instability of mass transfer. It is found that with decreasing initial primary mass from 0.89M$_{\rm \odot}$ to 0.63M$_{\rm \odot}$, the range of the initial mass ratio decreases for detached binaries that experience stable mass transfer and evolve into contact. If the initial primary mass is less than 0.63M$_{\rm \odot}$, detached binaries would experience dynamically unstable mass transfer when the primaries of detached binaries fill their Roche lobes. These systems would evolve into a common envelope situation and probably then to a complete merger of two components on a quite short timescale. This results in a low mass limit at about 0.63M$_{\rm \odot}$ for the primary mass of contact binaries, which might be a main reason why the period distribution of contact binaries has a short limit of about 0.22 days. By comparing the theoretical period distribution of contact binaries with the observational data, it is found that the observed contact binaries are above the low mass limit for the primary mass of contact binaries and no observed contact binaries are below this limit. This suggests that the short-period limit of contact binaries can be explained by the instability of the mass transfer that occurs when the primaries of the initially detached binaries fill their Roche lobes.

preprint2011arXiv

Tidally enhanced stellar wind: a way to make the symbiotic channel to Type Ia supernova viable

In the symbiotic (or WD+RG) channel of the single degenerate scenario for type Ia supernovae (SNe Ia) the explosions occur a relatively long time after star formation. The birthrate from this channel would be too low to account for all observed SNe Ia were it not for some mechanism to enhance the rate of accretion on to the white dwarf. A tidally enhanced stellar wind, of the type which has been postulated to explain many phenomena related to giant star evolution in binary systems, can do this. Compared to mass-stripping this model extends the space of SNe Ia progenitors to longer orbital periods and hence increases the birthrate to about $0.0069\,{\rm yr^{-1}}$ for the symbiotic channel. Two symbiotic stars, T CrB and RS Oph, considered to be the most likely progenitors of SNe Ia through the symbiotic channel, are well inside the period--companion mass space predicted by our models.

preprint2010arXiv

A likely candidate of type Ia supernova progenitors: the X-ray pulsating companion of the hot subdwarf HD 49798

HD 49798 is a hydrogen depleted subdwarf O6 star and has an X-ray pulsating companion (RX J0648.0-4418). The X-ray pulsating companion is a massive white dwarf. Employing Eggleton's stellar evolution code with the optically thick wind assumption, we find that the hot subdwarf HD 49798 and its X-ray pulsating companion could produce a type Ia supernova (SN Ia) in future evolution. This implies that the binary system is a likely candidate of SN Ia progenitors. We also discussed the possibilities of some other WD + He star systems (e.g. V445 Pup and KPD 1930+2752) for producing SNe Ia.

preprint2010arXiv

Adiabatic Mass Loss in Binary Stars - I. Computational Method

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 time scale mass transfer, and to establish the range of post-common envelope binaries that are allowed energetically.

preprint2010arXiv

Birthrates and delay times of Type Ia supernovae

Type Ia supernovae (SNe Ia) play an important role in diverse areas of astrophysics, from the chemical evolution of galaxies to observational cosmology. However, the nature of the progenitors of SNe Ia is still unclear. In this paper, according to a detailed binary population synthesis study, we obtained SN Ia birthrates and delay times from different progenitor models, and compared them with observations. We find that the Galactic SN Ia birthrate from the double-degenerate (DD) model is close to those inferred from observations, while the birthrate from the single-degenerate (SD) model accounts for only about 1/2-2/3 of the observations. If a single starburst is assumed, the distribution of the delay times of SNe Ia from the SD model is a weak bimodality, where the WD + He channel contributes to the SNe Ia with delay times shorter than 100Myr, and the WD + MS and WD + RG channels to those with age longer than 1Gyr.

preprint2010arXiv

Companion stars of Type Ia supernovae and single low-mass white dwarfs

Recent investigations of the WD + MS channel of Type Ia supernovae (SNe Ia) imply that this channel may be the main contribution to the old population (>1Gyr) of SNe Ia. In the WD + MS channel, the WD could accrete material from a main-sequence or a slightly evolved star until it reaches the Chandrasekhar mass limit. The companions in this channel would survive after SN explosion and show distinguishing properties. In this Letter, based on SN Ia production regions of the WD + MS channel and three formation channels of WD + MS systems, we performed a detailed binary population synthesis study to obtain the properties of the surviving companions. The properties can be verified by future observations. We find that the surviving companions of the old SNe Ia have a low mass, which provides a possible way to explain the formation of the population of single low-mass WDs (<0.45Msun).

preprint2010arXiv

Habitable Zones and UV Habitable Zones around Host Stars

Ultraviolet radiation is a double-edged sword to life. If it is too strong, the terrestrial biological systems will be damaged. And if it is too weak, the synthesis of many biochemical compounds can not go along. We try to obtain the continuous ultraviolet habitable zones, and compare the ultraviolet habitable zones with the habitable zones of host stars. Using the boundary ultraviolet radiation of ultraviolet habitable zone, we calculate the ultraviolet habitable zones of host stars with masses from 0.08 to 4.00 \mo. For the host stars with effective temperatures lower than 4,600 K, the ultraviolet habitable zones are closer than the habitable zones. For the host stars with effective temperatures higher than 7,137 K, the ultraviolet habitable zones are farther than the habitable zones. For hot subdwarf as a host star, the distance of the ultraviolet habitable zone is about ten times more than that of the habitable zone, which is not suitable for life existence.

preprint2010arXiv

Habitable Zones of Host Stars During the Post-MS Phase

A star will become brighter and brighter with stellar evolution, and the distance of its habitable zone will become farther and farther. Some planets outside the habitable zone of a host star during the main sequence phase may enter the habitable zone of the host star during other evolutionary phases. A terrestrial planet within the habitable zone of its host star is generally thought to be suited to life existence. Furthermore, a rocky moon around a giant planet may be also suited to life survive, provided that the planet-moon system is within the habitable zone of its host star. Using Eggleton's code and the boundary flux of habitable zone, we calculate the habitable zone of our Solar after the main sequence phase. It is found that Mars' orbit and Jupiter's orbit will enter the habitable zone of Solar during the subgiant branch phase and the red giant branch phase, respectively. And the orbit of Saturn will enter the habitable zone of Solar during the He-burning phase for about 137 million years. Life is unlikely at any time on Saturn, as it is a giant gaseous planet. However, Titan, the rocky moon of Saturn, may be suitable for biological evolution and become another Earth during that time. For low-mass stars, there are similar habitable zones during the He-burning phase as our Solar, because there are similar core masses and luminosities for these stars during that phase.

preprint2010arXiv

Observations of Doppler Boosting in Kepler Lightcurves

Among the initial results from Kepler were two striking lightcurves, for KOI 74 and KOI 81, in which the relative depths of the primary and secondary eclipses showed that the more compact, less luminous object was hotter than its stellar host. That result became particularly intriguing because a substellar mass had been derived for the secondary in KOI 74, which would make the high temperature challenging to explain; in KOI 81, the mass range for the companion was also reported to be consistent with a substellar object. We re-analyze the Kepler data and demonstrate that both companions are likely to be white dwarfs. We also find that the photometric data for KOI 74 show a modulation in brightness as the more luminous star orbits, due to Doppler boosting. The magnitude of the effect is sufficiently large that we can use it to infer a radial velocity amplitude accurate to 1 km/s. As far as we are aware, this is the first time a radial-velocity curve has been measured photometrically. Combining our velocity amplitude with the inclination and primary mass derived from the eclipses and primary spectral type, we infer a secondary mass of 0.22+/-0.03 Msun. We use our estimates to consider the likely evolutionary paths and mass-transfer episodes of these binary systems.

preprint2010arXiv

On the formation of single and binary helium-rich sdO stars

We propose a formation channel for the previously unexplained helium-rich subdwarf O (He-rich sdO) stars in which post-subdwarf B (sdB) stars (i.e. hybrid COHe white dwarfs) reignite helium burning in a shell after gaining matter from their helium white-dwarf (WD) companions. Such short-period binaries containing post-sdB WDs and helium WDs are predicted by one of the major binary formation channels for sdB stars. In the majority of cases, mass transfer is expected to lead to a dynamically unstable merger event, leaving a single-star remnant. Calculations of the evolution of these stars show that their properties are consistent with the observed He-rich sdO stars. The luminosity of these stars is about an order of magnitude higher than that of canonical sdB stars. We also suggest that binary systems such as PG 1544+488 (Ahmad et al. 2004) and HE 0301-3039 (Lisker et al. 2004), which each contain two hot subdwarfs, could be the outcome of a double-core common-envelope phase. Since this favours intermediate-mass progenitors, this may also explain why the subdwarfs in these systems are He-rich.

preprint2010arXiv

Probability Distribution of Terrestrial Planets in Habitable Zones around Host Stars

With more and more exoplanets being detected, it is paid closer attention to whether there are lives outside solar system. We try to obtain habitable zones and the probability distribution of terrestrial planets in habitable zones around host stars. Using Eggleton's code, we calculate the evolution of stars with masses less than 4.00 \mo. We also use the fitting formulae of stellar luminosity and radius, the boundary flux of habitable zones, the distribution of semimajor axis and mass of planets and the initial mass function of stars. We obtain the luminosity and radius of stars with masses from 0.08 to 4.00 \mo, and calculate the habitable zones of host stars, affected by stellar effective temperature. We achieve the probability distribution of terrestrial planets in habitable zones around host stars. We also calculate that the number of terrestrial planets in habitable zones of host stars is 45.5 billion, and the number of terrestrial planets in habitable zones around K type stars is the most, in the Milky Way.

preprint2010arXiv

Stellar adiabatic mass loss model and applications

Roche-lobe overflow and common envelope evolution are very important in binary evolution, which is believed to be the main evolutionary channel to hot subdwarf stars. The details of these processes are difficult to model, but adiabatic expansion provides an excellent approximation to the structure of a donor star undergoing dynamical time scale mass transfer. We can use this model to study the responses of stars of various masses and evolutionary stages as potential donor stars, with the urgent goal of obtaining more accurate stability criteria for dynamical mass transfer in binary population synthesis studies. As examples, we describe here several models with the initial masses equal to 1 Msun and 10 Msun, and identify potential limitations to the use of our results for giant-branch stars.

preprint2010arXiv

The helium star donor channel for the progenitors of type Ia supernovae and their surviving companion stars

The nature of type Ia supernovae (SNe Ia) is still unclear. Employing Eggleton's stellar evolution code with the optically thick wind assumption, we systematically studied the He star donor channel of SNe Ia, in which a carbon-oxygen white dwarf accretes material from a He main-sequence star or a He subgiant to increase its mass to the Chandrasekhar mass. We mapped out the initial parameters for producing SNe Ia in the orbital period--secondary mass plane for various WD masses from this channel. According to a detailed binary population synthesis approach, we find that this channel can produce SNe Ia with short delay times (~100Myr) implied by recent observations. We obtained many properties of the surviving companions of this channel after SN explosion, which can be verified by future observations. We also find that the surviving companions from the SN explosion scenario have a high spatial velocity (>400km/s), which could be an alternative origin for hypervelocity stars (HVSs), especially for HVSs such as US 708.

preprint2010arXiv

The helium star donor channel for the progenitors of type Ia supernovae with different metallicities

Context: The nature of type Ia supernovae (SNe Ia) is still unclear. Metallicities may have an important effect on their properties. Aims: In this paper, we study the He star donor channel towards SNe Ia comprehensively and systematically at various metallicities. Methods: Employing Eggleton's stellar evolution code with the optically thick wind assumption, we calculated about 10 000 WD + He star systems and obtained SN Ia production regions of the He star donor channel with metallicities Z=0.03, 0.02, 0.004 and 0.0001. According to a detailed binary population synthesis approach, we also obtained SN Ia birthrates at various metallicities. Results: Our study shows that both the initial mass of the He donor star and the initial orbital period for SNe Ia increase with metallicity, while the minimum initial mass of the carbon--oxygen white dwarfs producing SNe Ia decreases with metallicity. For a constant star-formation galaxy, SN Ia birthrates increase with metallicity. If a single starburst is assumed, SNe Ia occur systemically earlier and the peak value of the birthrate is larger for a high Z, and the He star donor channel with different metallicities can produce the young SNe Ia with delay times $\sim$45-220Myr.

preprint2010arXiv

WD+RG systems as the progenitors of Type Ia supernovae

Type Ia supernovae (SNe Ia) play an important role in the study of cosmic evolution, especially in cosmology. There are several progenitor models for SNe Ia proposed in the past years. In this paper, by considering the effect of accretion disk instability on the evolution of white dwarf (WD) binaries, we performed detailed binary evolution calculations for the WD + red-giant (RG) channel of SNe Ia, in which a carbon-oxygen WD accretes material from a RG star to increase its mass to the Chandrasekhar mass limit. According to these calculations, we mapped out the initial and final parameters for SNe Ia in the orbital period--secondary mass plane for various WD masses for this channel. We discussed the influence of the variation of the duty cycle value on the regions for producing SNe Ia. Similar to previous studies, this work also indicates that the long-period dwarf novae offer a possible ways for producing SNe Ia. Meanwhile, we find that the surviving companion stars from this channel have a low mass after SN explosion, which may provide a means for the formation of the population of single low-mass WDs (<0.45Msun).

preprint2009arXiv

Companion stars of Type Ia supernovae and hypervelocity stars

{Context} Recent investigations of the white dwarf (WD) + He star channel of Type Ia supernovae (SNe Ia) imply that this channel can produce SNe Ia with short delay times. The companion stars in this channel would survive and be potentially identifiable. {Aims} In this Letter, we study the properties of the companion stars of this channel at the moment of SN explosion, which can be verified by future observations. {Methods} According to SN Ia production regions of the WD + He star channel and three formation channels of WD + He star systems, we performed a detailed binary population synthesis study to obtain the properties of the surviving companions. {Results} We obtained the distributions of many properties of the companion stars of this channel at the moment of SN explosion. We find that the surviving companion stars have a high spatial velocity (>400 km/s) after SN explosion, which could be an alternative origin for hypervelocity stars (HVSs), especially for HVSs such as US 708.

preprint2009arXiv

Evolving to Type Ia Supernovae with Short Delay Times

The single-degenerate model is currently a favourable progenitor model for Type Ia supernovae (SNe Ia). Recent investigations on the WD + He star channel of the single-degenerate model imply that this channel is noteworthy for producing SNe Ia. In this paper we studied SN Ia birthrates and delay times of this channel via a detailed binary population synthesis approach. We found that the Galactic SN Ia birthrate from the WD + He star channel is $\sim 0.3\times 10^{-3} {\rm yr}^{-1}$ according to our standard model, and that this channel can explain SNe Ia with short delay times ($\sim4.5\times10^7 - 1.4\times10^8$ yr). Meanwhile, these WD + He star systems may be related to the young supersoft X-ray sources prior to SN Ia explosions.

preprint2009arXiv

The energy transfer and its effects on the secondaries in W Ursae Majoris-type contact binaries

The energy transfer of W UMa contact binaries and its effects on the secondary in W UMa contact binaries are investigated. Relations are given between the mass ratio (q) and the relative energy transfer rates, i.e. U1, the ratio of the transferred luminosity to the surface luminosity of the primary and U2, the ratio of the transferred luminosity to the nuclear luminosity of the secondary. The theoretical curves(U1 vs q and U2 vs q) are derived. Although these curves can reflect the distribution of U1 vs q and U2 vs q, some observational systems are significantly deviated from these curves. It is mainly resulted from the difference in the effective temperatures of the components in W UMa systems.The radius and the density of the secondary are related to the relative energy tran-sfer rate U2: the higher is U2, the greater is the expansion and the lower is the density of the secondaries in W UMa systems. The temperature difference of W UMa binary components is correlated with the relative energy transfer rate U1 and decreases with increasing U1. This might suggest that there is a thermal coupling between two components in W UMa contact binaries, and that the classif-ication of W UMa contact binaries into A- or W-types depends on the energy tran-sfer from the primary to the secondary. The temperature difference of W UMa bin-ary components is poorly correlated with the mass of the primary. This suggests that the properties of the common envelope of W UMa contact binaries might not have a significant effect on the energy transfer between two components.

preprint2009arXiv

WD + He star systems as the progenitors of Type Ia supernovae and their surviving companion stars

Employing Eggleton's stellar evolution code with an optically thick wind assumption, we have systematically studied the WD + He star channel of Type Ia supernovae (SNe Ia), in which a carbon-oxygen WD accretes material from a He main-sequence star or a He subgiant to increase its mass to the Chandrasekhar mass. We mapped out the parameter spaces for producing SNe Ia. According to a detailed binary population synthesis approach, we find that the Galactic SN Ia birthrate from this channel is $\sim$$0.3\times 10^{-3} {\rm yr}^{-1}$, and that this channel can produce SNe Ia with short delay times ($\sim45-$140 Myr). We also find that the surviving companion stars in this channel have a high spatial velocity ($>$400 km/s) after SN explosion, which could be an alternative origin for hypervelocity stars (HVSs), especially for HVSs such as US 708.

preprint2006arXiv

Binary Stellar Population Synthesis Study of Elliptical Galaxies

We determine relative stellar ages and metallicities mainly for about 80 elliptical galaxies in low and high density environments via the latest binary stellar population (BSP) synthesis model and test a latest hierarchical formation model of elliptical galaxies which adopted the new $Λ$CDM cosmology for the first time. The stellar ages and metallicities of galaxies are estimated from two high-quality published spectra line indices, i.e. H$β$ and [MgFe]. The results show that elliptical galaxies have stellar populations older than 3.9 Gyr and more metal rich than 0.02. Most of our results are in agreement with predictions of the model: First, elliptical galaxies in denser environment are redder and have older populations than field galaxies. Second, elliptical galaxies with more massive stellar components are redder while have older and more metal rich populations than less massive ones. Third, the most massive galaxies are shown to have the oldest and most metal rich stars. However, some of our results are found to be different with predictions of the galaxy formation model, i.e. the metallicity distributions of low- and high-density elliptical galaxies and the relations relating to cluster-centric distance.