Search

Double helium white dwarfs (He WDs) are one type of gravitational wave source and are greatly important in the studies of binary interaction, particularly in the common envelope (CE) ejection physics. Most double He WDs with mass ratios of q~1 are formed through a particular channel. In this channel, one He WD is initially produced from a red giant (RG) with a degenerate core via stable Roche lobe overflow, and another He WD is formed from an RG with a degenerate core via CE ejection. They may have significant implications for the binary evolution processes yet have not received specific studies, especially for the CE phase. This paper adopts a semi-analytic method and a detailed stellar evolution simulation to model the formation of double He WDs. We find that most double He WDs show mass ratios being slightly greater than 1, and their orbital periods and mass ratios relation are broadly consistent with observations. There is also a relation between the mass ratios and the progenitors' masses of the He WDs produced via CE ejection for double He WDs with determined WD masses. Based on this relation, the mass of the He WD progenitor can be inferred from the mass ratio. Then, the CE ejection efficiency can be constrained with the orbital period. In addition, we constrain the CE ejection efficiency for two double He WDs, J1005-2249 and WD0957-666. The results show that the CE ejection efficiencies increase with the WD progenitor masses., Comment: 12 pages,10 figures, Accepted for publication in ApJ

We present a study of the detached eclipsing binary TV~Mon using spectra from the LAMOST medium-resolution survey and ASAS-SN, CoRoT photometry. We applied multiple-epochs spectral fitting to derive RV and spectral parameters. The analysis of eclipses in CoRoT data told us relative sizes of the stellar components and almost edge-on circular orbit. Combining spectral and photometrical solution we estimated masses and radii of the components: $M_{A,B}=2.063\pm0.033,~0.218\pm0.004~M_\odot$, $R_{A,B}=2.427\pm0.014,~2.901\pm0.016~R_\odot$. SED analysis and Gaia parallax allowed us to get estimation of temperatures $T_{A,B}=7624^{+194}_{-174},~5184^{+130}_{-123}$ K and distance $d=907\pm11$ pc. We identified three $\delta$ Scuti type pulsation frequencies in primary component, while we also suspect TV~Mon having a long period variability with period $P_{\rm long}\sim128$ days and spot activity in secondary component. This system experienced intensive mass transfer and mass ratio reversal in the past, currently showing no signs of mass transfer in the spectra. The low mass component will lose its outer envelope and shrink to the helium white dwarf, which mass and orbital period are in good agreement with evolutionary models predictions., Comment: submitted in MNRAS

Binary stars are responsible for many unusual astrophysical phenomena, including some important explosive cosmic events. The stability criteria for rapid mass transfer and common-envelope evolution are fundamental to binary star evolution. They determine the mass, mass ratio, and orbital distribution of systems such as X-ray binaries and merging gravitational-wave sources. We use our adiabatic mass-loss model to systematically survey metal-poor and solar-metallicity donor thresholds for dynamical timescale mass transfer. The critical mass ratios qad are systematically explored, and the impact of metallicity and nonconservative mass transfer are studied. For metal-poor radiative-envelope donors, qad are smaller than those for solar-metallicity stars at the same evolutionary stage. However, qad do the opposite for convective-envelope donors. Nonconservative mass transfer significantly decreases qad for massive donors. This is because it matters how conservative mass transfer is during the thermal timescale phase immediately preceding a delayed dynamical mass transfer. We apply our theoretical predictions to observed high-mass X-ray binaries that have overfilled their Roche lobes and find a good agreement with their mass ratios. Our results can be applied to study individual binary objects or large samples of binary objects with binary population synthesis codes., Comment: Submitted to ApJ. Comments are welcome

The measurement of atmospheric parameters is fundamental for scientific research using stellar spectra. The Chinese Space Station Telescope (CSST), scheduled to be launched in 2024, will provide researchers with hundreds of millions of slitless spectra for stars during a 10 yr survey. And machine learning has unparalleled efficiency in processing large amounts of data compared to manual processing. Here we studied the stellar parameters of early-type stars (effective temperature Teff more than 15,000 K) based on the design indicators of the CSST slitless spectrum and the machine learning algorithm, Stellar LAbel Machine. We used the Potsdam Wolf-Rayet (POWR) synthetic spectra library for cross validation. Then we tested the reliability of machine learning results by using the Next Generation Spectrum Library (NGSL) from Hubble Space Telescope observation data. We use the spectra with the impact of interstellar extinction (AV = 0, 0.5, 1, 1.5, 2 mag) and radial velocity (RV = -50, -30, 0, 30, 50 km s-1) from the POWR library as the test set. When RV = 0 km s-1 and AV = 0 mag, the average value and standard deviation for 3 wavelength ranges (2550-4050 Ang (R = 287); 4050-6300 Ang (R = 232); 6300-10000 Ang (R = 207)) are -66 K, 550 K, and 356 K for Teff, and 0.004 c.g.s, -0.024 c.g.s, and 0.01 c.g.s for log g. When using the observed data from NGSL as the testing samples, the deviation of Teff is less than 5%, and the deviation of log g is less than 11%. In addition, we also test the influence of shifting of spectra on the parameters accuracy. The deviation of Teff for the case with a shift of 5 Ang and 10 Ang are 3.6% and 4.3%, respectively; the deviation of log g are 4.2% and 5.1%. These results demonstrate that we can obtain relatively accurate stellar parameters of a population of early-type stars with the CSST slitless spectra and a machine-learning method.

Type Ia Supernovae (SNe Ia) discovered at redshift $z\lesssim2.5$ are presumed to be produced from Population (Pop) I/II stars. {In this work, we investigate the production of SNe Ia from Pop III binaries in the cosmological framework. We derive the SN Ia rate as a function of redshift under a theoretical context for the production of first generation stars and emanate the likelihood of their detection by the James Webb Space Telescope (JWST).} {Assuming the initial stellar mass function (IMF) favors low-mass stars as from recent numerical simulations, we found Pop III stars may give rise to a considerable amount of SNe Ia at high redshift and Pop III stars may even be the dominant SN Ia producer at z $\gtrsim 6$.} {In an optimistic scenario, we expect $\sim 1(2)$ SNe Ia from Pop III stars at $z\approx 4(5)$ for a survey of area $300 \;\rm arcmin^2$ during a $3\;\rm yr$ period with JWST. The same survey may record more than $\sim 400$ SNe Ia at lower redshift ($z\lesssim 2.5$) but with only about one of them from Pop III progenitors. There will be $\sim 6$ Pop III SNe Ia in the same field of view at redshifts of $5-10$.} Observational constraints on SN Ia rates at the redshift range of $5-10$ can place crucial constraints on the IMF of Pop III stars., Comment: 19 pages, 7 Figures. Minor updates to match the published version

The unstable mass transfer situation in binary systems will asymptotically cause the adiabatic expansion of the donor star and finally lead to the common envelope phase. This process could happen in helium binary systems once the helium donor star fills its Roche-lobe. We have calculated the adiabatic mass loss model of naked helium stars with a mass range of 0.35\,$M_{\odot}$ to 10\,$M_{\odot}$, and every mass sequence evolved from the He-ZAMS to the cooling track of white dwarf or carbon ignition. In consideration of the influence of stellar wind, massive helium stars are not considered in this paper. Comparing stellar radius with the evolution of the Roche-lobe under the assumption of conservative mass transfer, we give the critical mass ratio $q_{\textrm{crit}}=M_{\textrm{He}}/M_{\textrm{accretor}}$ as the binary stability criteria of low and intermediate-mass helium binary stars. On He-MS, the result shows $1.0

Low mass-ratio (q) contact binary systems are progenitors of stellar mergers such as blue straggles (BS) or fast-rotating FK Com stars. In this study, we present the first light curve analysis of two newly identified low mass-ratio contact binary systems, TIC 55007847 and TIC 63597006, that are identified from TESS. Both stars are classified as A-subtype contact binaries. We obtained the precise orbit periods for the two objects by using the O-C method, i.e. P=0.6117108 d for TIC 55007847 and P=0.7008995 d for TIC 63597006, respectively, and found an obvious periodic signal in the O-C curve of TIC 63597006. We suggest that the periodic signal comes from a third body. We further use the Markov Chain Monte Carlo (MCMC) method with PHOEBE to derive the photometric solutions for the two binaries. The photometric solution for this object shows that the contribution of the third body is about 6%. Our analysis revealed that TIC 55007847 has an extremely low mass ratio of q=0.08. By calculating the ratio of spin angular momentum to the orbital angular momentum Js/Jo, we found that TIC 55007847 is very close to the instability threshold with Js/Jo = 0.31, indicating that it may merge into a single, fast-rotating star in the future. For TIC 63597006, q=0.14 and Js/Jo=0.15. This object is in a relatively stable evolutionary status at present.

Runaway stars are OB-type stars ejected from their birthplace with large peculiar velocities. The leading hypothesis addressed in their formation includes the supernova ejection mechanism and the dynamic ejection scenario. Identification of runaway populations is the first step to investigating their formation and evolution. Here we present our work of searching for Galactic runaway candidate stars from the LAMOST Medium-Resolution Survey DR8 database. After studying the kinematic properties for a collection of 4,432 early-type stars, predominantly B-type stars, using the radial velocity measurements from LAMOST DR8 and astrometric solutions made by Gaia DR3, we identified 229 runaway candidate stars. They span a wide distribution in projected rotational velocities. We investigated the Galactic spatial distribution of the runaway population and noticed that most of them likely reside within the Galactic thin disk. Based upon analyzing the Doppler shifts of the candidate stars, we found two binary runaway candidates displaying velocity variation with estimated orbital periods of 40 and 61 days.

Although supernovae is a well-known endpoint of an accreting white dwarf, alternative theoretical possibilities has been discussing broadly, such as the accretion-induced collapse (AIC) event as the endpoint of oxygen-neon (ONe) white dwarfs, either accreting up to or merging to excess the Chandrasekhar limit (the maximum mass of a stable white dwarf). AIC is an important channel to form neutron stars, especially for those unusual systems, which are hardly produced by core-collapse supernovae. However, the observational evidences for this theoretical predicted event and its progenitor are all very limited. In all of the known progenitors, white dwarfs increase in mass by accretion. Here, we report the discovery of an intriguing binary system Lan 11, consisted of a stripped core-helium-burning hot subdwarf and an unseen compact object of 1.08 to 1.35 $M_{\odot}$. Our binary population synthesis calculations, along with the absence of detection from the deep radio observations of the Five-hundred-meter Aperture Spherical Radio Telescope, strongly suggest that the latter is an ONe white dwarf. The total mass of this binary is 1.67 to 1.92 $M_{\odot}$}, significantly excessing the Chandrasekhar limit. The reproduction of its evolutionary history indicates that the unique system has undergone two phases of common envelope ejections, implying a born nature of this massive ONe white dwarf rather than an accretion growth from its companion. These results, together with short orbital period of this binary (3.65 hours), suggest that this system will merge in 500-540 Myr, largely triggering an AIC event, although the possibility of type Ia supernova cannot be fully ruled out. This finding greatly provides valuable constraints on our understanding of stellar endpoints, whatever leading to an AIC or a supernova., Comment: 25 pages, 14 figures

The development of spectroscopic survey telescopes like Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST), Apache Point Observatory Galactic Evolution Experiment and Sloan Digital Sky Survey has opened up unprecedented opportunities for stellar classification. Specific types of stars, such as early-type emission-line stars and those with stellar winds, can be distinguished by the profiles of their spectral lines. In this paper, we introduce a method based on derivative spectroscopy (DS) designed to detect signals within complex backgrounds and provide a preliminary estimation of curve profiles. This method exhibits a unique advantage in identifying weak signals and unusual spectral line profiles when compared to other popular line detection methods. We validated our approach using synthesis spectra, demonstrating that DS can detect emission signals three times fainter than Gaussian fitting. Furthermore, we applied our method to 579,680 co-added spectra from LAMOST Medium-Resolution Spectroscopic Survey, identifying 16,629 spectra with emission peaks around the H{\alpha} line from 10,963 stars. These spectra were classified into three distinct morphological groups, resulting in nine subclasses as follows. (1) Emission peak above the pseudo-continuum line (single peak, double peaks, emission peak situated within an absorption line, P Cygni profile, Inverse P Cygni profile); (2) Emission peak below the pseudo-continuum line (sharp emission peak, double absorption peaks, emission peak shifted to one side of the absorption line); (3) Emission peak between the pseudo-continuum line.

The hot subdwarf O/B stars (sdO/Bs) are known as extreme horizontal branch stars, which is of great importance in stellar evolution theory. The sdO/Bs are generally thought to have a helium-burning core and a thin hydrogen envelope $(M_{\rm env }<0.02M_\odot)$. In the canonical binary evolution scenario, sdO/Bs are considered to be the stripped cores of red giants. However, such a scenario cannot explain the recently discovered sdO/B binary, SMSS J1920, where the strong Ca H$\&$K lines in the spectrum are found. It suggests that this binary is likely originated from the recent ejection of common envelope (CE). In this {work}, we proposed a new formation channel of massive sdO/Bs, namely sdO/Bs produced from a CE ejection process with an asymptotic giant branch (AGB) star (hereafter AGB CE channel). We constructed the evolutionary model of sdO/Bs and successfully explained most of the important observed parameters of the sdO/B star in SMSS J1920, including the evolutionary age, sdO/B mass, effective temperature, surface gravity and surface helium abundance. The minimum sdO/B mass produced from the AGB CE channel is about $0.48M_\odot$. The evolutionary tracks in $\log T_{\rm eff}-\log g$ plane {may explain a fraction of the observational samples} with high-$\log T_{\rm eff}$ and low-$\log g$. Considering wind mass-loss of sdO/Bs, the model could produce helium-rich hot subdwarfs with $\log (n_{\rm He}/n_{\rm H})\gtrsim-1$., Comment: 9 pages, 5 figures, accepted for publication in ApJ

Millisecond pulsars (MSPs) are a kind of radio pulsars with short spin periods, playing a key role in many aspects of stellar astrophysics. In recent years, some more MSPs with wide orbits ($>30\,\rm d$) have been discovered, but their origin is still highly unclear. In the present work, according to an adiabatic power-law assumption for the mass-transfer process, we carried out a large number of complete binary evolution computations for the formation of MSPs with wide orbits through the iron core-collapse supernova (CCSN) channel, in which a neutron star (NS) originating from a CCSN accretes matter from a red-giant (RG) star and spun up to millisecond periods. We found that this channel can form the observed MSPs with wide orbits in the range of $30-1200\,{\rm d}$, in which the WD companions have masses in the range of $0.28-0.55\,\rm M_{\odot}$. We also found that almost all the observed MSPs can be reproduced by this channel in the WD companion mass versus orbital period diagram. We estimate that the Galactic numbers of the resulting MSPs from the CCSN channel are in the range of $\sim 4.8-8.5\times10^{5}$. Compared with the accretion-induced collapse channel, the CCSN channel provides a main way to produce MSPs with wide orbits., Comment: 19 pages, 8 figures, 3 tables. Accepted for publication in MNRAS

Binary evolution theory predicts that the second common envelope (CE) ejection can produce low-mass (0.32-0.36 Msun) subdwarf B (sdB) stars inside ultrashort-orbital-period binary systems, as their helium cores are ignited under nondegenerate conditions. With the orbital decay driven by gravitational-wave (GW) radiation, the minimum orbital periods of detached sdB binaries could be as short as ~20 minutes. However, only four sdB binaries with orbital periods below an hour have been reported so far, while none of them has an orbital period approaching the above theoretical limit. Here we report the discovery of a 20.5-minute-orbital-period ellipsoidal binary, TMTS J052610.43+593445.1, in which the visible star is being tidally deformed by an invisible carbon-oxygen white dwarf (WD) companion. The visible component is inferred to be an sdB star with a mass of ~0.33 Msun, approaching that of helium-ignition limit, although a He-core WD cannot be completely ruled out. In particular, the radius of this low-mass sdB star is only 0.066 Rsun, about seven Earth radii, possibly representing the most compact nondegenerate star ever known. Such a system provides a key clue to map the binary evolution scheme from the second CE ejection to the formation of AM CVn stars having a helium-star donor, and it will also serve as a crucial verification binary of space-borne GW detectors in the future., Comment: 24 pages, 11 figures, 1 table, published on Nature Astronomy, URL: https://www.nature.com/articles/s41550-023-02188-2

Thanks to the mid-infrared capability offered by Spitzer, systematic searches of dust in SNe have been carried out over the past decade. Studies have revealed the presence of a substantial amount of dust over a broad range of SN subtypes. How normal SNe present mid-IR excess at later time and turn out to be dusty SNe can be affected by several factors, such as mass-loss history and envelope structure of progenitors and their explosion environment. All these can be combined and related to their environmental properties. A systematic analysis of SNe that exploded under a dusty environment could be of critical importance to measure the properties of the dust-veiled exploding stars, and whether such an intense dust production process is associated with the local environment. In this work, we firstly use the IFS data to study the environmental properties of dusty SNe compared to those of normal ones, and analyze correlations between the environmental properties and their dust parameters. We find that dusty SNe have a larger proportion located at higher SFR regions compared to the normal types. The occurrence of dusty SNe is less dependent on metallicity, with the oxygen abundance spanning from subsolar to oversolar metallicity. We also find the host extinction of dusty SNe scatters a lot, with about 40% of dusty SN located at extremely low extinction environments, and another 30% of them with considerably high host extinction of E(B-V)>0.6 mag., Comment: 20 pages, 15 figures, submitted to ApJ

The common envelope evolution (CEE) is vital in forming short orbital period compact binaries. It covers many objects, such as double compact merging binaries, type Ia supernovae progenitors, binary pulsars, and X-ray binaries. Knowledge about the common envelope (CE) eject efficiency still needs to be improved, though progress has been made recently. Short orbital period hot subdwarf B star plus white dwarf binaries are the most straightforward samples to constrain CEE physics. We apply the known orbital period-white dwarf relation to constrain the sdB progenitor of seven sdB+WD binaries with a known inclination angle. The average value of the CE efficiency parameter is 0.32, which is consistent with previous studies. However, the CE efficiency might not be a constant but is a function of the initial mass ratio based on well-constrained sdB progenitor mass and evolutionary stage. Our results can be used as physical inputs for binary population synthesis simulations on related objects. A similar method can also be applied to study other short orbital period WD binaries., Comment: 16 pages, 12 figures, 2 tables, accepted for publication in The Astrophysical Journal

Binary stars are as common as single stars. Binary stars are of immense importance to astrophysicists because that they allow us to determine the masses of the stars independent of their distances. They are the cornerstone of the understanding of stellar evolutionary theory and play an essential role in cosmic distance measurement, galactic evolution, nucleosynthesis and the formation of important objects such as cataclysmic variable stars, X-ray binaries, Type Ia supernovae, and gravitational wave-producing double compact objects. In this article, we review the significant theoretical and observational progresses in addressing binary stars in the new millennium. Increasing large survey projects have led to the discovery of enormous numbers of binary stars, which enables us to conduct statistical studies of binary populations, and therefore provide unprecedented insight into the stellar and binary evolution physics. Meanwhile, the rapid development of theoretical concepts and numerical approaches for binary evolution have made a substantial progress on the alleviation of some long-standing binary-related problems such as the stability of mass transfer and common envelope evolution. Nevertheless, it remains a challenge to have a full understanding of fundamental problems of stellar and binary astrophysics. The upcoming massive survey projects and increasingly sophisticated computational methods will lead to future progress., Comment: An invited review published in Progress in Particle and Nuclear Physics; An open access to the published version, see https://doi.org/10.1016/j.ppnp.2023.104083

Semi-detached binaries are in the stage of mass transfer and play a crucial role in studying mass transfer physics between interacting binaries. Large-scale time-domain surveys provide massive light curves of binary systems, while Gaia offers high-precision astrometric data. In this paper, we develop, validate, and apply a pipeline that combines the MCMC method with a forward model and DBSCAN clustering to search for semi-detached binary and estimate its inclination, relative radius, mass ratio, and temperature ratio using light curve. We train our model on the mock light curves from PHOEBE, which provides broad coverage of light curve simulations for semi-detached binaries. Applying our pipeline to TESS sectors 1-26, we have identified 77 semi-detached binary candidates. Utilizing the distance from Gaia, we determine their masses and radii with median fractional uncertainties of ~26% and ~7%, respectively. With the added 77 candidates, the catalog of semi-detached binaries with orbital parameters has been expanded by approximately 20%. The comparison and statistical results show that our semi-detached binary candidates align well with the compiled samples and the PARSEC model in Teff-L and M-R relations. Combined with the literature samples, comparative analysis with stability criteria for conserved mass transfer indicates that ~97.4% of samples are undergoing nuclear-timescale mass transfer, and two samples (GO Cyg and TIC 454222105) are located within the limits of stability criteria for dynamical- and thermal-timescale mass transfer, which are currently undergoing thermal-timescale mass transfer. Additionally, one system (IR Lyn) is very close to the upper limit of delayed dynamical-timescale mass transfer.

Hierarchical triple-star systems consists of three components organised into an inner binary ($M_{1}$,$M_{2}$) and a more distant outer tertiary ($M_{3}$) star. The LAMOST Medium-Resolution Spectroscopic Survey (LAMOST-MRS) has offered a great sample for the study of triple system populations. We used the Peak Amplitude Ratio (PAR) method to obtain the mass ratio ($q_\mathrm{{in}}$, $q_\mathrm{{out}}$) of a triple system from its normalised spectrum. By calculating Cross-Correlation Function (CCF), we determined the correlation between the mass ratio $q_\mathrm{{out}}$ ($M_{3}$/($M_{1}$+$M_{2}$)) and the amplitude ratio ($A_{3}$/($A_{1}$+$A_{2}$)). We derived $q_\mathrm{{in}}$ of $0.5-1.0$ and $q_\mathrm{{out}}$ between 0.2 and 0.8. By fitting a power-law function of the corrected $q_\mathrm{{in}}$ distribution, the $\gamma_\mathrm{{in}}$ are estimated to be $-0.654\pm2.915$, $4.304\pm1.125$ and $11.371\pm1.309$ for A, F and G type stars. The derived $\gamma_\mathrm{{in}}$-values increase as the mass decrease, indicating that less massive stars are more likely to have companion stars with similar masses. By fitting a power-law function of the corrected $q_\mathrm{{out}}$ distribution, the ${\gamma_\mathrm{{out}}}$ are estimated to be $-2.016\pm0.172$, $-1.962\pm0.853$ and $-1.238\pm0.141$ for G, F and A type stars, respectively. The ${\gamma_\mathrm{{out}}}$-values show a trend of growth toward lower primary star masses.

We use an updated method for the detection of double-lined spectroscopic binaries (SB2) using $v \sin{i}$ values from spectral fits. The method is applied to all spectra from LAMOST-MRS. Using this method, we detect 12426 SB2 candidates, where 4321 are already known and 8105 are new discoveries. We check their spectra manually to minimise possible false positives. We also detect several cases of contamination of the spectra by solar light. Additionally, for candidates with multiple observations we compute mass ratios with systemic velocities and determine Keplerian orbits. We present an updated catalogue of all SB2 candidates together with additional information for some of them in separate data tables., Comment: accepted in MNRAS 2023-10-16. arXiv admin note: substantial text overlap with arXiv:2207.06996

We selected three double-lined spectroscopic binary systems which have extreme mass ratios, if measured using the Wilson method. We analysed medium resolution spectroscopic observations and space-based photometry and find that all these systems are not SB2, but rather triple systems and a chance alignment of another star with SB1 that have an unseen component. Therefore suspicious mass ratios determined by the Wilson method for some double-lined spectroscopic binary systems aren't correct as these systems are more complex., Comment: 10 pages, accepted in MNRAS

The forthcoming CSST wide-field multiband imaging survey will produce seven-band photometric spectral energy distributions (SEDs) for billions of galaxies. The effective extraction of astronomical information from these massive datasets of SEDs relies on the techniques of both SED synthesis (or modeling) and analysis (or fitting). We evaluate the performance of the latest version of BayeSED code combined with SED models with increasing complexity for simultaneously determining the photometric redshifts and stellar population parameters of galaxies in this survey. By using an empirical statistics-based mock galaxy sample without SED modeling errors, we show finding that the random observational errors in photometries are more important sources of errors than the parameter degeneracies and Bayesian analysis method and tool. By using a Horizon-AGN hydrodynamical simulation-based mock galaxy sample with SED modeling errors about the star formation histories (SFHs) and dust attenuation laws (DALs), the simple typical assumptions lead to significantly worse parameter estimation with CSST photometries only. The SED models with more flexible (or complicated) forms of SFH/DAL do not necessarily lead to better estimation of redshift and stellar population parameters. We discuss the selection of the best SED model by means of Bayesian model comparison in different surveys. Our results reveal that the Bayesian model comparison with Bayesian evidence may favor SED models with different complexities when using photometries from different surveys. Meanwhile, the SED model with the largest Bayesian evidence tends to give the best performance of parameter estimation, which is more clear for photometries with larger discriminative power., Comment: Accepted for publication in ApJS (49 pages, 23 figures, 5 tables). Comments are welcome! The new version of BayeSED code, documents, and the scripts used for the performance tests presented in this work will be publicly available at https://bitbucket.org/hanyk/bayesed/, https://bayesed.readthedocs.io/, and https://github.com/hanyk/BayeSED-performance-test/, respectively

Wide binaries play a crucial role in analyzing the birth environment of stars and the dynamical evolution of clusters. When wide binaries are located at greater distances, their companions may overlap in the observed images, becoming indistinguishable and resulting in unresolved wide binaries, which are difficult to detect using traditional methods. Utilizing deep learning, we present a method to identify unresolved wide binaries by analyzing the point-spread function (PSF) morphology of telescopes. Our trained model demonstrates exceptional performance in differentiating between single stars and unresolved binaries with separations ranging from 0.1 to 2 physical pixels, where the PSF FWHM is ~2 pixels, achieving an accuracy of 97.2% for simulated data from the Chinese Space Station Telescope. We subsequently tested our method on photometric data of NGC 6121 observed by the Hubble Space Telescope. The trained model attained an accuracy of 96.5% and identified 18 wide binary candidates with separations between 7 and 140 au. The majority of these wide binary candidates are situated outside the core radius of NGC 6121, suggesting that they are likely first-generation stars, which is in general agreement with the results of Monte Carlo simulations. Our PSF-based method shows great promise in detecting unresolved wide binaries and is well suited for observations from space-based telescopes with stable PSF. In the future, we aim to apply our PSF-based method to next-generation surveys such as the China Space Station Optical Survey, where a larger-field-of-view telescope will be capable of identifying a greater number of such wide binaries., Comment: 17 pages, 12 figures. Accepted for publication in ApJS

Double white dwarf (WD) merger process and their post-merger evolution are important in many fields of astronomy, such as supernovae, gamma-ray bursts, gravitational waves, etc. The evolutionary outcomes of double ultra-massive WD merger remnants are still a subject of debate, though the general consensus is that the merger remnant will collapse to form a neutron star. In this work, we investigate the evolution of a 2.20Msun merger remnant stemmed from the coalescence of double 1.10Msun ONe WDs. We find that the remnant ignites off-centre neon burning at the position near the surface of primary WD soon after the merger, resulting in the stable inwardly propagating oxygen/neon (O/Ne) flame. The final outcomes of the merger remnant are sensitive to the effect of convective boundary mixing. If the mixing cannot stall the O/Ne flame, the flame will reach the centre within 20 years, leading to the formation of super Chandrasekhar mass silicon core, and its final fate probably be neutron star (NS) through iron-core-collapse supernova. In contrast, if the convective mixing is effective enough to prevent the O/Ne flame from reaching the centre, the merger remnant will undergo electron capture supernova to form an ONeFe WD. Meanwhile, we find that the wind mass loss process may hardly alter the final fate of the remnant due to its fast evolution. Our results imply that the coalescence of double ONe WDs can form short lived giant like object, but the final outcomes (NS or ONeFe WD) are influenced by the uncertain convective mixing in O/Ne flame., Comment: 17 pages, 7 figures, a revised version after referee's comments. Comments are welcome

We present the detection of false positive double-lined spectroscopic binaries candidates (SB2) using medium-resolution survey (MRS) spectra from the one time-domain field of LAMOST data release 10 (DR10). The secondary component in all these binaries has near zero radial velocity and solar-like spectral lines. Highly likely this is light from the semi-transparent clouds illuminated by the full Moon. However we also suspect that partially this contamination can be caused by a solar light reflected from the surface of low-orbital artificial satellites launched in the beginning of 2022. We found several possible contaminant candidates using archival orbital data. We propose measures to reduce risk of such contamination for the future observations and methods to find it in archived ones., Comment: accepted in MNRAS letters 2023 July 07

SNe Ia play a key role in the fields of astrophysics and cosmology. It is widely accepted that SNe Ia arise from thermonuclear explosions of WDs in binaries. However, there is no consensus on the fundamental aspects of the nature of SN Ia progenitors and their explosion mechanism. This fundamentally flaws our understanding of these important astrophysical objects. We outline the diversity of SNe Ia and the proposed progenitor models and explosion mechanisms. We discuss the recent theoretical and observational progress in addressing the SN Ia progenitor and explosion mechanism in terms of the observables at various stages of the explosion, including rates and delay times, pre-explosion companion stars, ejecta-companion interaction, early excess emission, early radio/X-ray emission from CSM interaction, surviving companions, late-time spectra and photometry, polarization signals, and SNR properties, etc. Despite the efforts from both the theoretical and observational side, the questions of how the WDs reach an explosive state and what progenitor systems are more likely to produce SNe Ia remain open. No single published model is able to consistently explain all observational features and the full diversity of SNe Ia. This may indicate that either a new progenitor paradigm or the improvement of current models is needed if all SNe Ia arise from the same origin. An alternative scenario is that different progenitor channels and explosion mechanisms contribute to SNe Ia. In the next decade, the ongoing campaigns with the JWST, Gaia and the ZTF, and upcoming extensive projects with the LSST and the SKA will allow us to conduct not only studies of individual SNe Ia in unprecedented detail but also systematic investigations for different subclasses of SNe Ia. This will advance theory and observations of SNe Ia sufficiently far to gain a deeper understanding of their origin and explosion mechanism., Comment: An invited review (accepted for publication in RAA)

The formation path to ultra-compact X-ray binaries (UCXBs) with black hole (BH) accretors is still unclear. In the classical formation scenario, it is difficult to eject the massive envelope of the progenitor star of the BH via common envelope process. Given that some neutron stars (NSs) in binary systems evidently have birth masses close to $\sim 2.0\;M_\odot$, we explore here the possibility that BH-UCXBs may form via accretion-induced collapse (AIC) of accreting NSs, assuming that these previously evolved in LMXBs to masses all the way up to the maximum limit of a NS. We demonstrate this formation path by modelling a few cases of NS-UCXBs with initial NS masses close to the maximum mass of a NS that evolve into BH-UCXBs after the NS accretes material from its He~WD companion. We follow the evolution of the post-AIC BH-UCXB and, based on simple arguments, we anticipate that there is about one BH-UCXB with an AIC origin and a He~WD donor within the current sample of known UCXBs and that 2--5 such BH-UCXBs may be detected in gravitational waves by LISA. In addition, we find that the X-ray luminosity of NS-UCXBs near their orbital period minimum exceeds $\sim 10^{39}\;{\rm erg\;s^{-1}}$ and thus such systems may appear as ultraluminous X-ray sources., Comment: ApJ accepted, 12 pages, 3 figures

Close binary interactions may play a critical role in the formation of the rapidly rotating Be stars. Mass transfer can result in a mass gainer star spun up by the accretion of mass and angular momentum, while the mass donor is stripped of its envelope to form a hot and faint helium star. FUV spectroscopy has led to the detection of about 20 such binary Be+sdO systems. Here we report on a three-year program of high quality spectroscopy designed to determine the orbital periods and physical properties of five Be binary systems. These binaries are long orbital period systems with $P =$ 95 to 237 days and with small semi-amplitude $K_1<11$ km s$^{-1}$. We combined the Be star velocities with prior sdO measurements to obtain mass ratios. A Doppler tomography algorithm shows the presence of the He II $\lambda 4686$ line in the faint spectrum of the hot companion in four of the targets. We discuss the observed line variability and show evidence of phased-locked variations in the emission profiles of HD 157832, suggesting a possible disk spiral density wave due to the presence of the companion star. The stripped companions in HD 113120 and HD 137387 may have a mass larger than the 1.4 $M_\odot$ indicating that they could be progenitors of Type Ib and Ic supernovae., Comment: 50 pages, 23 figures, 16 tables

Although multidimensional simulations have investigated the processes of double WD mergers, post-merger evolution only focused on the carbon-oxygen (CO) WD or helium (He) WD merger remnants. In this work, we investigate for the first time the evolution of the remnants stemmed from the merger of oxygen-neon (ONe) WDs with CO WDs. Our simulation results indicate that the merger remnants can evolve to hydrogen- and helium-deficient giants with maximum radius of about 300Rsun. Our models show evidence that merger remnants more massive than 1.95Msun can ignite Ne before significant mass-loss ensues, and they thus would become electron-capture supernovae (ECSNe). However, remnants with initial masses less than 1.90Msun will experience further core contraction and longer evolutionary time before reaching at the conditions for Ne-burning. Therefore their fates are more dependent on mass-loss rates due to stellar winds, and thus more uncertain. Relatively high mass-loss rates would cause such remnants to end their lives as ONe WDs. Our evolutionary models can naturally explain the observational properties of the double WD merger remnant IRAS 00500+6713 (J005311). As previously suggested in the literature, we propose and justify that J005311 may be the remnant from the coalescence of an ONe WD and an CO WD. We deduce that the final outcome of J005311 would be a massive ONe WD rather than a supernova explosion. Our investigations may be able to provide possible constraints on the wind mass-loss properties of the giants which have CO-dominant envelopes., Comment: 22 pages, 10 figures, 8 tables, accepted for publication in ApJL

The stability criteria of rapid mass transfer and common envelope evolution are fundamental in binary star evolution. They determine the mass, mass ratio and orbital distribution of many important systems, such as X-ray binaries, Type Ia supernovae and merging gravitational wave sources. We use our adiabatic mass-loss model to systematically survey the intermediate-mass stars' thresholds for dynamical-timescale mass transfer. The impact of metallicity on the stellar responses and critical mass ratios is explored. Both tables ($Z=0.001$) and fitting formula ($Z=0.001$ and $Z=0.02$) of critical mass ratios of intermediate-mass stars are provided. An application of our results to intermediate-mass X-ray binaries (IMXBs) is discussed. We find that the predicted upper limit to mass ratios, as a function of orbital period, is consistent with the observed IMXBs that undergo thermal or nuclear timescale mass transfer. According to the observed peak X-ray luminosity $L_\mathrm{X}$, we predict the range of $L_\mathrm{X}$ for IMXBs as a function of the donor mass and the mass transfer timescale., Comment: 16 pages, 19 figures, 3 tables, accepted for publication in The Astrophysical Journal

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 $\alpha_{CE}$. While the results of the polytropic model strongly depend on the adopted value of $\alpha_{ 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., Comment: 28 pages, 17 figures; A&A version

We present a study of eclipsing binary J064726.39+223431.6 using spectra from the LAMOST-MRS and TESS photometry. We use full-spectrum fitting to derive radial velocities and spectral parameters: ${T_{\rm eff}}_{A,B}=6177,\,5820$ K, $v \sin{i}_{A,B}=59,\,50~\kms$ and ${\rm [Fe/H]}_{A,B}=-0.19$ dex. The orbital solution and light curve analysis suggest that it is a close pair of fast rotating stars on circular orbit. We measure their masses to be $M_{A,B}=1.307\pm0.007,\, 1.129\pm0.005\,M_\odot$ and their radii to be $R_{A,B}=1.405\pm0.052,\, 1.219\pm0.060\,R_\odot$ resulting in surface gravities of $\log{g}_{A,B}=4.259\pm0.033,\,4.319\pm0.042$ (cgs). Theoretical models cannot match all of these properties, predicting significantly higher ${T_{\rm eff}}$ for a given metallicity. Derived age of the system 1.56 Gyr indicates that both components are younger than Sun, which contradicts to much longer orbit's circularisation timescale., Comment: submitted to MNRAS. arXiv admin note: text overlap with arXiv:2204.10680

Massive binary stars play a crucial role in many astrophysical fields. Investigating the statistical properties of massive binary stars is essential to trace the formation of massive stars and constrain the evolution of stellar populations. However, no consensus has been achieved on the statistical properties of massive binary stars, mainly due to the lack of a large and homogeneous sample of spectroscopic observations. We study the intrinsic binary fraction $f_{\rm b}^{\rm in}$ and distributions of mass ratio $f(q)$ and orbital period $f(P)$ of early-type stars (comprised of O-, B-, and A-type stars) and investigate their dependences on effective temperature $T_{\rm eff}$, stellar metallicity [M/H], and the projection velocity $v\sin{i}$, based on the homogeneous spectroscopic sample from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) Data Release Eight (DR8). We found that $f_{\rm b}^{\rm in}$ increases with increasing $T_\mathrm{eff}$. The binary fraction is positively correlated with metallicity for spectra in the sample. Over all the $v\sin{i}$ values we considered, the $f_{\rm b}^{\rm in}$ have constant values of $\sim$50\%. It seems that the binary population is relatively evenly distributed over a wide range of $v\sin{i}$ values, while the whole sample shows that most of the stars are concentrated at low values of $v\sin{i}$ (probably from strong wind and magnetic braking of single massive stars) and at high values of $v\sin{i}$ (likely from the merging of binary stars). Stellar evolution and binary interaction may be partly responsible for this.There are no correlations found between $\pi$($\gamma$) and $T_{\rm eff}$, nor for $\pi$($\gamma$) and [M/H]. The uncertainties of the distribution decrease toward a larger sample size with higher observational cadence.

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}$., Comment: 13 pages, 15 figures, to be published in MNRAS

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., Comment: 10 pages, 8 figures, accepted for publication in ApJ

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., Comment: 12 pages, 6 figures, accepted for publication in Research in Astronomy and Astrophysics

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., Comment: 13 pages, 7 figures, accepted for publication in ApJ

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 $\beta_\mathrm{CE}$ becomes smaller than $\alpha_\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} \alpha_{\mathrm{CE}}$ and $\log_{10} \beta_{\mathrm{CE}}$ linearly correlate in sdBs, similarly to De Marco et al. (2010) for post-AGB binaries., Comment: 18 pages, 17 figures, 3 tables, Accepted for publication in ApJ, minor updates for the final proof

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 $\gamma$ 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 $\gamma$-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., Comment: 18 pages, 12 figures, accepted for publication in The Astrophysical Journal

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., Comment: Accepted in MNRAS 22.04.2022. Received 21.04.2022; in original form 17.11.2021

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., Comment: 13 pages, 5 figures, 2 tables, accepted for publication in MNRAS, comments welcome!

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., Comment: Accepted by ApJ, 9 pages, 7 figures

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., Comment: 34 pages, 12 figures, 11 tables

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., Comment: 16 pages, 10 figures. Published by The Innovation

With the third observing run the number of gravitational wave emitting events has increased significantly. The data of the recorded events is inspected to search for overall properties on the population. The properties of subpopulations are determined and compared to predictions from simulations. It appears that the most outstanding systems follow linear relations in the parameter space of the total binary and the chirp mass. Those relations are too tight to have a stochastic origin and are supported by at least five independent events each. The origin of the correlations is still open to be confirmed, while possible sources, ranging from instrumental artefacts to unknown physics, are discussed and partly excluded. Depending on the relation's source the two events (GW190814 and GW200210_092254) having a smaller mass component between 2.5 and 2.9 Msun may reveal in a very different light., Comment: more checks on the relations added; 6+8 pages, 2+10 figures, 2+6 tables, submitted to A&A on Dec. 8 2021

Double white dwarf (DWD) binaries are important for studies of common-envelope (CE) evolution, Type Ia supernova progenitors and Galactic sources of low-frequency gravitational waves (GWs). PTF J0533+0209 is a DWD system with a short orbital period of Porb ~ 20 min and thus a so-called LISA verification source. The formation of this system and other DWDs is still under debate. In this paper, we discuss the possible formation scenarios of this binary and argue that it is not likely to have formed through CE evolution. Applying a new magnetic braking prescription, we use the MESA code to model the formation of this system through stable mass transfer. We find a model which can well reproduce the observed WD masses and orbital period but not the effective temperature and hydrogen abundance of the low-mass He WD component. We discuss the possibility of using H flashes to mitigate this discrepancy. Finally, we discuss the future evolution of this system into a AM CVn binary such as those that will be detected by space-borne GW observatories like LISA, TianQin and Taiji., Comment: Accepted by ApJ, 11 pages, 4 figures

Massive stars play key roles in many astrophysical processes. Deriving atmospheric parameters of massive stars is important to understand their physical properties and thus are key inputs to trace their evolution. Here we report our work on adopting the data-driven technique Stellar LAbel Machine ({\tt SLAM}) with the non-LTE TLUSTY synthetic spectra as the training dataset to estimate the stellar parameters of LAMOST optical spectra for early-type stars. We apply two consistency tests to verify this machine learning method and compare stellar labels given by {\tt SLAM} with that in literature for several objects having high-resolution spectra. We provide the stellar labels of effective temperature ($T_\mathrm{eff}$), surface gravity ($\log{g}$), metallicity ([M/H]), and projected rotational velocity ($v\sin{i}$) for 3,931 and 578 early-type stars from LAMOST Low-Resolution Survey (LAMOST-LRS) and Medium-Resolution Survey (LAMOST-MRS), respectively. To estimate the average statistical uncertainties of our results, we calculated the standard deviation between the predicted stellar label and the pre-labeled published values from the high-resolution spectra. The uncertainties of the four parameters are $\sigma(T_\mathrm{eff}) = 2,185 $K, $\sigma(\log{g}) = 0.29$ dex, and $\sigma(v\sin{i}) = 11\, \rm km\,s^{-1}$ for MRS, and $\sigma(T_\mathrm{eff}) = 1,642 $K, $\sigma(\log{g}) = 0.25$ dex, and $\sigma(v\sin{i}) = 42\, \rm km\,s^{-1}$ for LRS spectra, respectively. We notice that parameters of $T_\mathrm{eff}$, $\log{g}$ and [M/H] can be better constrained using LRS spectra rather than using MRS spectra, most likely due to their broad wavelength coverage, while $v\sin{i}$ is constrained better by MRS spectra than by LRS spectra, probably due to the relatively accurate line profiles of MRS spectra.