Researcher profile

Zhanwen Han

Zhanwen Han contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
28works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

28 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&#39; 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&#39;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.

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.