Search

The formation of compact binary systems is largely driven by their evolution through a common envelope (CE) phase, crucial for understanding phenomena such as type Ia supernovae and black hole mergers. Despite their importance, direct observational evidence for CE material has been elusive due to the transient nature of these envelopes. Numerical simulations suggest that some envelope material may persist post-ejection. In this study, we investigate circumstellar material (CSM) surrounding hot subdwarf (sdB) stars, focusing on material ejected during the CE phase of binary evolution. We analyze Ca II K absorption lines in 727 sdB candidates from the LAMOST-LRS survey, selecting 145 stars with strong absorption features, indicating the presence of CSM. We compare the velocities of the Ca II K lines with the systemic velocities of sdB binaries, confirming that the material originates from ejected common-envelope material. The results show that the CSM persists long after the CE event, suggesting the formation of a stable, long-lived circumstellar environment around sdB stars. This study enhances our understanding of the role of CSM in post-CE evolution and provides new insights into the physical processes shaping the evolution of sdB binaries., Comment: 13 pages, 10 figures, Accepted for publication in MNRAS

Linkage between core-collapse supernovae (SNe) and their progenitors is not fully understood and ongoing effort of searching and identifying the progenitors is needed. $\mathrm{SN\,2024abfl}$ is a recent Type II supernova exploded in the nearby star-bursting galaxy $\mathrm{NGC\,2146}$, which is also the host galaxy of $\mathrm{SN\,2018zd}$. From archival Hubble Space Telescope (HST) data, we have found a red source ($\mathrm{m_{F814W} \sim 25}$) near the location (angular distance $\leq 0.2"$) of $\mathrm{SN\,2024abfl}$ before its explosion. With F814W and F606W photometry, we found that the properties of this source matched a typical red supergiant (RSG) moderately reddened by interstellar dust at the distance of the host galaxy. We conclude that the $\mathrm{SN\,2024abfl}$ had an RSG progenitor with initial mass of $\mathrm{10M_{\odot}}$--$\mathrm{16\,M_{\odot}}$., Comment: 11 pages, 6 figures, in review (submitted on Dec 9th, 2024)

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.