Search

The growing ``Hubble tension'' has prompted the need for precise measurements of cosmological distances. This paper demonstrates a purely geometric approach for determining the distance to extragalactic binaries through a joint analysis of spectroastrometry (SA), radial velocity (RV), and light curve (LC) observations. A parameterized model for the binary system is outlined, and simulated SA, RV, and LC data are computed to infer the probability distribution of model parameters based on the mock data. The impact of data quality and binary parameters on distance uncertainties is comprehensively analyzed, showcasing the method's potential for high-precision distance measurements. For a typical eclipsing binary in the Large Magellanic Cloud (LMC), the distance uncertainty is approximately 6% under reasonable observational conditions. Within a specific range of data quality and input parameters, the distance measurement precision of individual binary star systems is generally better than 10%. As a geometric method based on the simplest dynamics, it is independent of empirical calibration and the systematics caused by model selections can be tested using nearby binaries with known distances. By measuring multiple binary star systems or monitoring one binary system repeatedly, geometric distance measurements of nearby galaxies can be achieved, providing valuable insights into the Hubble tension and advancing our understanding of the universe's structure and evolution., Comment: 16 pages, 9 figures, accepted by ApJ

EL CVn-type systems represent a rare evolutionary stage in binary star evolution, providing ideal laboratories for investigating stable mass transfer processes and the formation of extremely low-mass white dwarfs (ELM WDs). The Transiting Exoplanet Survey Satellite (TESS) has delivered an extensive collection of high-precision time-domain photometric data, which is invaluable for studying EL CVn binaries. In this study, we identified 29 EL CVn systems from the TESS eclipsing binary catalogs (sectors 1-65), 11 of which are newly discovered. These systems consist of smaller, hotter pre-He white dwarfs and A/F main-sequence stars. The orbital periods of these binaries range from 0.64 to 2.5 days. Utilizing TESS light curves, Gaia distances, and multi-band photometric data (e.g., GALEX, 2MASS, WISE, SkyMapper), we modeled the light curves and spectral energy distributions to derive system parameters, including effective temperatures, masses, and radii. These systems were then compared with the white dwarf mass-period relation and the evolutionary tracks of ELM WDs. The comparison reveals that these binaries are consistent with the expected mass-period relation for white dwarfs and align well with the evolutionary tracks on the Teff-logg diagram for ELM WDs. This result suggests that these EL CVn systems likely formed through stable mass transfer processes. We provide a catalog of complete parameters for 29 EL CVn systems identified from the TESS survey. This catalog will serve as an essential resource for studying binary mass transfer, white dwarf formation, and pulsation phenomena in EL CVn-type systems., Comment: 13 pages, 10 figures, 2 tables

Accurate measurements of europium abundances in cool stars are essential for an enhanced understanding of the r-process mechanisms. We measure the abundance of Eu in solar spectra and a sample of metal-poor stars in the Galactic halo and metal-poor disk, with the metallicities ranging from \GG{$-2.4$} to $-0.5$ dex, using non-local thermodynamic equilibrium (NLTE) line formation. We compare these measurements with Galactic Chemical Evolution (GCE) models to \GG{explore the impact of the NLTE corrections on the contribution of r-process site in Galactic chemical evolution. In this work, we use NLTE line formation, as well as one-dimensional (1D) hydrostatic and spatial averages of 3D hydrodynamical ($<$3D$>$) model atmospheres to measure the abundance of Eu based on both the Eu II 4129 \AA\ and Eu II 6645 \AA\ lines for solar spectra and metal-poor stars. We find that \GG{for Eu II 4129 \AA\ line the NLTE modelling leads to higher (0.04 dex) solar Eu abundance in 1D and higher (0.07 dex) in \GG{$<$3D$>$} NLTE while} NLTE modelling leads to higher (0.01 dex) solar Eu abundance in 1D and lower (0.03 dex) in \GG{$<$3D$>$} NLTE for Eu II 6645 \AA\ line. Although the NLTE corrections for the Eu II $\lambda$ 4129 \AA\ and Eu II $\lambda$ 6645 \AA\ lines are opposite, the discrepancy between the abundances derived from these individual lines reduces after applying NLTE corrections, highlighting the critical role of NLTE abundance determinations. By comparing these measurements with Galactic chemical evolution (GCE) models, we find that the \G{amount of NLTE correction does not require significant change of the parameters for Eu production} in the GCE models.

Triple star systems are critical for understanding stellar dynamics and compact objects in astrophysics, yet confirmed hierarchical triples identified via spectroscopy remain limited. In this study, we identified 23 triple systems by cross-matching the Gaia DR3 non-single star catalog with LAMOST DR10 spectroscopic data; 18 of them are new discoveries. For two well-observed triples, we performed radial velocity curve fitting and light curve analysis to determine their orbital parameters, with inner and outer periods of 1.26 days and 656 days for one triple, and 3.42 days and 422 days for the other. We compared the results with other studies. We also analyzed the radial velocities (RVs) of these 23 triples, revealing a range of $V$ from approximately 40~km~s$^{-1}$ to 210~km~s$^{-1}$. Due to spectral resolution and detection limitations, velocity differences below 45~km~s$^{-1}$ in binaries and below 90~km~s$^{-1}$ in the inner binaries of triple systems are challenging to detect. Consequently, our detection range for inner orbital periods is restricted to 0.2--20 days, with the highest efficiency for periods under 10 days. These findings underscore the advantage of spectroscopic observations for identifying triple systems with short inner orbital periods.

We present an analysis of LAMOST J171013.211+532646.04 (hereafter J1710), a binary system comprising a hot subdwarf B star (sdB) and a white dwarf (WD) companion. Multi-epoch spectroscopy reveals an orbital period of 109.20279 minutes, consistent with TESS and ZTF photometric data, marking it as the sixth detached system known to harbor a WD companion with a period less than two hours. J1710 is remarkably close to Earth, situated at a distance of only \(350.68^{+4.20}_{-4.21} \, \mathrm{pc}\), with a GAIA G-band magnitude of 12.59, rendering it conducive for continuous observations. The spectral temperature is around 25164 K, in agreement with SED fitting results (\(25301^{+839}_{-743} \, \mathrm{K}\)). The TESS light curve displays ellipsoidal variation and Doppler beaming without eclipsing features. Through fitting the TESS light curve using the Wilson-Devinney code, we determined the masses for the sdB (\(M_1 = 0.44^{+0.06}_{-0.07} \, M_{\odot}\)) and the compact object (\(M_2 = 0.54^{+0.10}_{-0.07} \, M_{\odot}\)), with the compact object likely being a WD. Furthermore, MESA models suggest that the sdB, with a helium core mass of 0.431 \(M_{\odot}\) and a hydrogen envelope mass of \(1.3 \times 10^{-3}\, M_{\odot}\), is in the early helium main-sequence phase. The MESA binary evolution shows that the J1710 system is expected to evolve into a double white dwarf system, making it an important source of low-frequency gravitational waves., Comment: 11 pages, 8 figures, 2tables, accepted in Astronomy & Astrophysics journal

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 apply multiple-epochs spectral fitting to derive RV and spectral parameters. The analysis of eclipses in CoRoT data show the relative sizes of the stellar components and almost edge-on circular orbit. Combining the spectral and photometrical solutions we estimate masses and radii of the components: $M_{A,B}=2.063\pm0.033({\rm stat.})\pm0.095({\rm syst.}),~0.218\pm0.004({\rm stat.})\pm0.018({\rm syst.})~M_\odot$, $R_{A,B}=2.394\pm0.014,~2.860\pm0.016~R_\odot$. SED analysis and Gaia parallax allow us to get estimation of temperatures ${T_{ eff}}_{A,B}=7624^{+194}_{-174},~5184^{+130}_{-123}$ K and distance $d=907\pm11$ pc. We identify three $\delta$ Scuti type pulsation frequencies in the primary component, while we also suspect TV~Mon having a spot activity in the secondary component. This system experienced intensive mass transfer and mass ratio reversal in the past, but currently shows no signs of mass transfer in the spectra. The low mass component will lose its outer envelope and shrink to the helium white dwarf, the mass and orbital period of which are in good agreement with evolutionary models predictions., Comment: accepted in MNRAS 31.10.2024

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

The properties of H-rich, type II-plateau supernova (SN II-P) progenitors remain uncertain, and this is primarily due to the complexities associated with red supergiant (RSG) wind-mass loss. Recent studies have suggested that the interaction of the ejecta with a standard RSG wind should produce unambiguous signatures in the optical (e.g., a broad, boxy H${\alpha}$ profile) and in the UV (especially Ly ${\alpha}$ and Mg ii ${\lambda}{\lambda}$ 2795, 2802) a few years following the explosion. Such features are expected to be generic in all SNe II-P and can be utilized to constrain RSG winds. Here, we investigate the possibility of detecting late-time (0.3-10 years since explosion) SNe II-P in the NUV with the China Space Station Telescope (CSST). Convolving the existing model spectra of ejecta-wind interactions in SNe II-P with the transmission functions of the CSST, we calculated the associated multiband light curves, in particular, the NUV (255 nm${\sim}$317 nm) band, as well as the $NUV-r$ color. We find that the CSST will be able to detect the NUV radiation associated with ejecta-wind interaction for hundreds SNe II-P out to a few hundred Mpc over its ten-year main sky survey. The CSST will therefore provide a sizable sample of SNe II-P with the NUV signatures of ejecta-wind interaction. This will be helpful for understanding the mass loss history of SN II-P progenitors and their origins., Comment: 10 pages, 7 figures, accepted for publication in Astronomy & Astrophysics

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

Gravitational Waves (GWs) provide a unique way to explore our Universe. The ongoing ground-based detectors, e.g., LIGO, Virgo, and KAGRA, and the upcoming next-generation detectors, e.g., Cosmic Explorer and Einstein Telescope, as well as the future space-borne GW antennas, e.g., LISA, TianQin, and TaiJi, cover a wide range of GW frequencies {from $\sim 10^{-4}\;\rm Hz$ to $\sim 10^3\;\rm Hz$} and almost all types of compact objects in close orbits serve as the potential target sources for these GW detectors. The synergistic multi-band GW and EM observations would allow us to study fundamental physics from stars to cosmology. {The formation of stellar GW sources has been extensively explored in recent years, and progress on physical processes in binary interaction has been made as well. Furthermore, some studies have shown that the progress in binary evolution may significantly affect the properties of the stellar GW sources.}, Comment: 33 pages, 8 figures; Invited Review of Results in Physics; Accepted for publication; References updated

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

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

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

Double white dwarf systems are of great astrophysical importance in the field of gravitational wave and Type Ia supernova. While the binary fraction of CO core white dwarf is about a few percents, the extremely low mass white dwarfs are all thought to be within binary systems. In this work, we report the orbital solution of a double degenerate system: J033847.06+413424.24, an extremely low mass He core white dwarf orbiting a CO core white dwarf. With LAMOST and P200, time domain spectroscopic observations have been made and spectral atmosphere parameters are estimated to be $T_{\rm eff}\sim22500$ K and log $g\sim5.6$ dex. Combining Gaia parallax, 3D extinction, and evolution tracks, we estimate a radius of $\sim0.12$ $R_{\odot}$ and a mass of $\sim0.22$ $M_{\odot}$. With the 37 single exposure spectra, the radial velocities are measured and the orbital parameters are estimated to be $P=0.1253132(1)$ days, $K1=289\pm4$ km/s and $V_{sys}=-41\pm3$ km/s. The radial velocity based system ephemeris is also provided. The light curves from several photometric surveys show no orbital modulation. The orbital solution suggests that the invisible companion has a minimum mass of about 0.60 $M_{\odot}$ and is $\sim0.79$ $M_{\odot}$ for an inclination of $60.0^{\circ}$, indicating most probably a CO core white dwarf. The system is expected to merge in about 1 Gyr. With present period and distance ($\sim596$ pc) it can not irradiate strong enough gravitational wave for LISA. More double degenerate systems are expected to be discovered and parameterized as the LAMOST survey goes on., Comment: 12 pages, 11 figures

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

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

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

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

The Extremely Low Mass White Dwarfs (ELM WDs) and pre-ELM WDs are helium core white dwarfs with mass $<\sim 0.3M_{\odot}$. Evolution simulations show that a lower mass limit for ELM WDs exists at $\approx0.14M_{\odot}$ and no one is proposed by observation to be less massive than that. Here we report the discovery of a binary system, LAMOST J224040.77-020732.8 (J2240 in short), which consists of a very low mass hot star and a compact companion. Multi-epoch spectroscopy shows an orbital period $P_{orb} =$0.219658$\pm0.000002$ days and a radial velocity semi-amplitude $K1=318.5\pm3.3km/s$, which gives the mass function of 0.74$M_{\odot}$, indicating the companion is a compact star. The F-type low resolution spectra illustrate no emission features, and the temperature ($\sim 7400K$) is consistent with that from Spectral Energy Distribution fitting and multi-color light curve solution. The optical light curves, in ZTF g, r and i bands and Catalina V band, show ellipsoidal variability with amplitudes $\sim30\%$, suggesting that the visible component is heavily tidally distorted. Combining the distance from Gaia survey, the ZTF light curves are modeled with Wilson-Devinney code and the result shows that the mass of the visible component is $M1=0.085^{+0.036}_{-0.024}M_{\odot}$, and the mass of the invisible component is $M2=0.98^{+0.16}_{-0.09}M_{\odot}$. The radius of the visible component is $R1=0.29^{+0.04}_{-0.03}R_{\odot}$. The inclination angle is approximately between 60$^{\circ}$ and 90$^{\circ}$. The observations indicate the system is most likely a pre-ELM WD + WD/NS binary, and the mass of pre-ELM is possibly lower than the $0.14M_{\odot}$ theoretical limit., Comment: 18 pages, 10 figures

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

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

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

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.

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^\pi$ (1\textless$P$\textless1000 days) and $f(q) \propto q^\kappa$ (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 $\pi$ are -1$\pm0.1$, -1.1$\pm0.05$, -1.1$\pm0.1$, and -0.6$\pm0.05$, respectively. The $\kappa$ 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.

The maximum mass of neutron stars (NSs) is of great importance for constraining equations of state of NSs and understanding the mass gap between NSs and stellar-mass black holes. NSs in X-ray binaries would increase in mass by accreting material from their companions (known as recycling process), and the uncertainties in the accretion process give challenge to study the NS mass at birth. {In this work, we investigate the NS accreted mass with considering the effect of NS spin evolution and give the maximum accreted mass for NSs in the recycling process. By exploring a series of binary evolution calculations, we obtain the final NS mass and the maximum accreted mass for a given birth mass of NS and a mass transfer efficiency. Our results show that the NSs can accrete relatively more material for binary systems with the donor masses in the range of $1.8\sim 2.4M_\odot$, the NSs accrete relatively more mass when the remnant WD mass is in the range of $\sim 0.25-0.30M_\odot$, and the maximum accreted mass is positively correlated with the initial NS mass. For a $1.4M_\odot$ NS at birth with a moderate mass transfer efficiency of 0.3, the maximum accreted mass could be $0.27M_\odot$. The results can be used to estimate the minimum birth mass for systems with massive NSs in observations., Comment: 12 pages, 6 figures; accepted for publication in ApJ

Based on the two-minutes TESS data, we analyzed intrinsic oscillations of the primary component and identified seven confident independent $\delta$ Scuti frequencies ($f_1$, $f_2$, $f_3$, $f_4$, $f_7$, $f_{11}$, and $f_{12}$). Both of single-star evolutionary models and mass-accreting models are computed to reproduce the $\delta$ Scuti freqiencies. Fitting results of them match well with each other. The stellar parameters of the primary star yielded by asteroseismology are $M$ = $1.92^{+0.10}_{-0.02}$ $M_{\odot}$, $Z$ = 0.011$^{+0.006}_{-0.001}$, $R$ = $2.068^{+0.050}_{-0.007}$ $R_{\odot}$, $\log g$ = $4.090^{+0.010}_{-0.002}$, $T_{\rm eff}$ = $8346^{+244}_{-320}$ K, $L$ = $18.65^{+3.31}_{-2.82}$ $L_{\odot}$, which match well with the dynamic ones by the binary model. Furthermore, our asteroseismic results show that OO Dra is another Algol system that has just undergone the rapid mass-transfer stage. Fitting results of single-star evlutionary models indicate that the pulsator is helium-poor star with an age of 8.22$^{+0.12}_{-1.33}$ Myr, and the further mass-accreting models show that the primary star looks like an almost unevolved star formed by an extremely helium-poor mass accretion in Case A evolutionary scenario., Comment: 25 pages, 8 figures, 7 tables. accepted for publication by ApJ. arXiv admin note: text overlap with arXiv:2004.09916

Most subdwarf B (sdB) + Helium white dwarf (He WD) binaries are believed to be formed from a particular channel. In this channel, the He WDs are produced first from red giants (RGs) with degenerate cores via stable mass transfer and sdB stars are produced from RGs with degenerate cores via common envelope (CE) ejection. They are important for the studies of CE evolution, binary evolution, and binary population synthesis. However, the relation between WD mass and orbital period of sdB + He WD binaries has not been specifically studied. In this paper, we first use a semi-analytic method to follow their formation and find a WD mass and orbital period relation. Then we use a detailed stellar evolution code to model their formation from main-sequence binaries. We find a similar relation between the WD mass and orbital period, which is in broad agreement with observations. For most sdB + He WD systems, if the WD mass (orbital period) can be determined, the orbital period (WD mass) can be inferred with this relation and then the inclination angle can be constrained with the binary mass function. In addition, we can also use this relation to constrain the CE ejection efficiency and find that a relative large CE ejection efficiency is favoured. If both the WD and sdB star masses can be determined, the critical mass ratios of dynamically unstable mass transfer for RG binaries can also be constrained., Comment: 6 pages, 5 figures, accepted for publication in MNRAS

Close-orbit low-mass X-ray binaries (LMXBs), radio binary millisecond pulsars (BMSPs) with extremely low-mass helium WDs (ELM He~WDs) and ultra-compact X-ray binaries (UCXBs) are all part of the same evolutionary sequence. It is therefore of uttermost importance to understand how these populations evolve from one specie to another. Moreover, UCXBs are important gravitational wave (GW) sources and can be detected by future space-borne GW observatories. However, the formation and evolutionary link between these three different populations of neutron star (NS) binaries are not fully understood. In particular, a peculiar fine-tuning problem has previously been demonstrated for the formation of these systems. In this investigation, we test a newly suggested magnetic braking prescription and model the formation and evolution of LMXBs. We compute a grid of binary evolution models and present the initial parameter space of the progenitor binaries which successfully evolve all the way to produce UCXBs. We find that the initial orbital period range of LMXBs, which evolve into detached NS+ELM~He~WD binaries and later UCXBs, becomes significantly wider compared to evolution with a standard magnetic braking prescription, and thus helps to relieve the fine-tuning problem. However, we also find that formation of wide-orbit BMSPs is prohibited for strong versions of this new magnetic braking prescription, which therefore calls for a revision of the prescription. Finally, we present examples of the properties of UCXBs as Galactic GW sources and discuss their detection by the LISA, TianQin and Taiji observatories., Comment: 12 pages, 13 figures, accepted for publication in MNRAS

Binary interactions lead to the formation of intriguing objects, such as compact binaries, supernovae, gamma ray bursts, X-ray binaries, pulsars, novae, cataclysmic variables, hot subdwarf stars, barium stars, and blue stragglers. To study the evolution of binary populations and the consequent formation of these objects, many methods have been developed over the years, of which a robust approach named binary population synthesis (BPS) warrants special attention. This approach has seen widespread use in many areas of astrophysics, including but not limited to analyses of the stellar content of galaxies, research on galactic chemical evolution, and studies concerning star formation and cosmic re-ionization. In this review, we discuss the role of BPS, its general picture, and the various components that comprise it. We pay special attention to the stability criteria for mass transfer in binaries, as this stability largely determines the fate of binary systems. We conclude with our perspectives regarding the future of this field., Comment: 18 pages, 5 figures, accepted publication for RAA

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., Comment: 26 pages, 14 figures, 4 tables; Accepted for publication in ApJ

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., Comment: 25 pages, 10 figues. Submitted to Reseach in Astronomy and Astrophysics