Source author record

Xuefei Chen

Xuefei Chen 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

29works
6topics
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

29 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.

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

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

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

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

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

Mass transfer of low-mass binaries and chemical anomalies among unevolved stars in globular clusters

While it is well known that mass transfer in binaries can pollute the surfaces of the accretors, it is still unclear whether this mechanism can reproduce the observed chemical inhomogeneities in globular clusters. We study the surface abundances of the accretors in low-mass binaries, as a first step towards understanding whether mass transfer in low-mass binaries is one of the potential origins of the aforementioned abundance anomalies in globular clusters. We use the Modules for Experiments in Stellar Astrophysics code to calculate binary evolutionary models with different initial donor masses between 0.9 and 1.9 $ \rm {M}_\odot$ for an initial metallicity of Z = 0.0034. The results show that in some low-mass binary systems, the accretors exhibit peculiar chemical patterns when they are still unevolved stars, e.g. C and O depletion; Na and N enhancement; constant Mg, Al and C+N+O. The abundance patterns of the accretors are significantly different from their initial abundances (or that of normal single stars), and can match the observed populations. These abundance patterns strongly depend not only on the initial parameters of binaries (donor mass, mass ratio, and orbital period), but also on the assumptions regarding mass-transfer efficiency and angular momentum loss. These results support the hypothesis that mass transfer in low-mass binaries is at least partly responsible for the unevolved anomalous stars in globular clusters. More work on binary evolutionary models and binary population synthesis is required to fully evaluate the contribution of this scenario.

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.

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.

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.

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

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

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.

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.

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

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

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.

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.