Search

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.

Context. Open clusters (OCs) are valuable probes of stellar population characteristics. Their age and metallicity provide insights into the chemical enrichment history of the Milky Way. By studying the metallicity of OCs, we can explore the spatial distribution of composition across the Galaxy and understand stellar birth radii through chemical tagging. However, inferring the original positions of OCs remains a challenge. Aims. This study investigates the distribution of metallicity in the solar neighborhood using data from Gaia DR3 and LAMOST spectra. By measuring accurate ages and metallicities, we aim to derive birth radii and understand stellar migration patterns. Methods. We selected 1131 OCs within 3 kpc of the Sun from Gaia DR3 and LAMOST DR8 low-resolution spectra (R=1800). To correct the LAMOST data, we incorporated high-resolution spectra from GALAH DR3 (R=28000) using an artificial neural network. The average metallicity of the OCs was derived from reliable [Fe/H] values of their members. We examined the metallicity distribution across the Galaxy and calculated birth radii based on age and metallicity. Results. The correction method reduces the systematic offset in LAMOST data. We found a metallicity gradient as a function of Galactocentric distance and guiding radii. Comparisons with chemo-dynamic simulations show that observed metallicity values are slightly lower than predicted when uncertainties are ignored, but the metallicity gradients align with previous studies. We also inferred that many OCs near the Sun likely originated from the outer Galactic disk.

Based on astrometric measurements and spectral analysis from $Gaia$ DR3, two quiescent black hole (BH) binaries, $Gaia$ BH1 and BH2, have been identified. Their origins remain controversial, particularly for $Gaia$ BH1. By considering a rapidly rotating ($\omega/\omega_{\rm crit} = 0.8$) and strongly magnetized ($B_{\rm 0} = 5000$ G) merger product, we find that, at typical Galactic metallicity, the merger product can undergo efficient chemically homogeneous evolution (CHE). This results in the merger product having a significantly smaller radius during its evolution compared to that of a normally evolving massive star. Under the condition that the initial triple stability is satisfied, we use the Multiple Stellar Evolution (MSE) code and the MESA code to identify an initial hierarchical triple that can evolve into $Gaia$ BH1. It initially consists of three stars with masses of 9.03 $M_{\odot}$, 3.12 $M_{\odot}$, and 1 $M_{\odot}$, with inner and outer orbital periods of 2.21 days and 121.92 days, and inner and outer eccentricities of 0.41 and 0.45, respectively. This triple initially experiences triple evolution dynamics instability (TEDI) followed by Roche lobe overflow (RLOF). During RLOF, the inner orbit shrinks, and tidal effects gradually suppress the TEDI. Eventually, the inner binary undergoes a merger through contact (or collision). Finally, using models of rapidly rotating and strongly magnetic stars, along with standard core-collapse supernova (SN) or failed supernova (FSN) models, we find that a PMB consisting of an 12.11 $M_{\odot}$ merger product and a 1 $M_{\odot}$ companion star (originally an outer tertiary) can avoid RLOF. After a SN or FSN with a low ejected mass of $\sim$0.22 $M_{\odot}$ and a low kick velocity ($46^{+25}_{-33}$ ${\rm km/s}$ or $9^{+16}_{-8}$ ${\rm km/s}$), the PMB can form $Gaia$ BH1 in the Galactic disk., Comment: 11 pages, 3 figures, Accepted to APJL

The source of the Galactic Lithium (Li) has long been a puzzle. With the discovery of Li in novae, extensive research has been conducted. However, there still exists a significant disparity between the observed abundance of lithium in novae and the existing theoretical predictions. Using the Modules for Experiments in Stellar Astrophysics (MESA), we simulate the evolution of nova with element diffusion and appropriately increased the amount of 3^He in the mixtures. Element diffusion enhances the transport efficiency between the nuclear reaction zone and the convective region on the surface of the white dwarf during nova eruptions, which results in more 7^Be to be transmitted to the white dwarf surface and ultimately ejected. Compared to the previous predictions, the abundance of 7^Be in novae simulated in our model significantly increases. And the result is able to explain almost all observed novae. Using the method of population synthesis, we calculate Li yield in the Galaxy. We find that the Galactic occurrence rate of nova is about 130 yr^{-1}, and about 110M Li produced by nova eruption is ejected into the interstellar medium (ISM). About 73\% of Li in the Galactic ISM originates from novae, and approximately 15\%-20\% of the entire Galaxy. It means that novae are the important source of Li in the Galactic., Comment: 12 pages, 4 figures. Accepted for publication in Astrophysical Journal

Wolf-Rayet stars (WRs) are very important massive stars. However, their origin and the observed binary fraction within the entire WR population are still debated. We investigate some possible merger channels for the formation of WRs, including main sequence (MS)/ Hertzsprung Gap (HG) + MS, He + HG/ Giant Branch (GB). We find that many products produced via binary merger can evolve into WRs, the MS/ HG + MS merger channel can explain WRs with luminosities higher than $\sim 10^{5.4}$\,L$_{\odot}$, while the He + HG/ GB merger channel can explain low-luminosity WRs in the range of $10^{4.7}$\,L$_{\odot}$\,$\sim$\,$10^{5.5}$\,L$_{\odot}$. In the population synthesis analysis of WRs, we assume an initial binary fraction ($f_{\rm ini,bin}$) of 50\% and 100\% for massive stars. We also assume that MS/ HG + MS merger products are non-rotating or rapidly rotating ($\omega/\omega_{\rm crit}=0.8$). In different cases, the calculated single fractions of WRs range from $22.2\%$ to $60.6\%$ in the Milky Way (MW) and from $8.3\%$ to $70.9\%$ in the Large Magellanic Cloud (LMC). The current observations fall within the range of our calculations. When the merger product of MS/HG+MS rotates rapidly, we estimate that there are approximately 1015 to 1396 WRs in the MW and 128 to 204 WRs in the LMC. Our model also roughly reproduces the observed single-peak luminosity distribution of WRs in the MW. However, the weak bimodal luminosity distribution observed in the LMC is not reproduced in our model. We assess that this may be due to the model underestimating the mass-loss rate in the LMC. In conclusion, we consider that the binary merger is significant formation channel for WR formation, and can explain the observed high fraction of the single WRs in the total population., Comment: 14 pages, 8 figures , Accepted to 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

Neutrino losses play a crucial role in the evolution of massive stars. We study the neutrino luminosity of stars ranging from 20 to 90 M_{\odot} from Zero Age Main Sequence (ZAMS) to Fe Core Collapse (FeCC) with different rotation and metallicity in a neutrino Hertzsprung-Russell diagram. In our simulations, we consider {\omega}/{\omega}crit = 0 and 0.7 to represent non-rotation and high rotation, respectively, and set the metallicities to 0.014, 0.001, and 0.0001. During hydrogen burning stages, neutrino luminosity primarily originates from CNO cycle, and increases with higher stellar mass while decreasing with increasing metallicity. For the high metallicity models (Z = 0.014) during the helium burning stage, the reduction of the hydrogen envelope caused by a larger mass loss rate leads to a gradual decrease in neutrino luminosity. The rapid rotation results in extra mixing inside massive stars, which increases the neutrino luminosity during main sequence (MS), while decreases the neutrino luminosity during helium burning phase. Simultaneously, the rapid rotation also increases CO core mass, which enhances the neutrino luminosity during C and O burning phase. We also investigate the effect of neutrino magnetic moment (NMM) on the massive stars. We find that the energy loss caused by the NMM does not have effects on the evolutionary destiny of massive stars, and it does not significant change the compactness at the time of Fe core collapse., Comment: This paper has been accepted by Monthly Notices of the Royal Astronomical Society on October 3, 2023

Symbiotic stars are interacting binary systems, making them valuable for studying various astronomical phenomena, such as stellar evolution, mass transfer, and accretion processes. Despite recent progress in the discovery of symbiotic stars, a significant discrepancy between the observed population of symbiotic stars and the number predicted by theoretical models. To bridge this gap, this study utilized machine learning techniques to efficiently identify new symbiotic stars candidates. Three algorithms (XGBoost, LightGBM, and Decision Tree) were applied to a dataset of 198 confirmed symbiotic stars and the resulting model was then used to analyze data from the LAMOST survey, leading to the identification of 11,709 potential symbiotic stars candidates. Out of the these potential symbiotic stars candidates listed in the catalog, 15 have spectra available in the SDSS survey. Among these 15 candidates, two candidates, namely V* V603 Ori and V* GN Tau, have been confirmed as symbiotic stars. The remaining 11 candidates have been classified as accreting-only symbiotic star candidates. The other two candidates, one of which has been identified as a galaxy by both SDSS and LAMOST surveys, and the other identified as a quasar by SDSS survey and as a galaxy by LAMOST survey., Comment: 18 pages, 6 figures, accepted 11-Jul-2023 by RAA

The actinides, such as the uranium (U) element, are typically synthesized through the rapid neutron-capture process (r-process), which can occur in core-collapse supernovae or double neutron star mergers. There exist nine r-process giant stars exhibiting conspicuousUabundances, commonly referred to as U-rich giants. However, the origins of these U-rich giants remain ambiguous. We propose an alternative formation scenario for these U-rich giants whereby a red giant (RG) engulfs an Earth-like planet. To approximate the process of a RG engulfing an Earth-like planet, we employ an accretion model wherein the RG assimilates materials from said planet. Our findings demonstrate that this engulfment event can considerably enhance the presence of heavy elements originating from Earth-like planets on the surfaces of very metal-poor stars (Z = 0.00001), while its impact on solar-metallicity stars is comparatively modest. Importantly, the structural and evolutionary properties of both very metalpoor and solar-metallicity stars remain largely unaffected. Notably, our engulfment model effectively accounts for the observed U abundances in known U-rich giants. Furthermore, the evolutionary trajectories of U abundances on the surfaces of RGs subsequent to the engulfment of Earth-like planets encompass all known U-rich giants. Therefore, it is plausible that U-rich giants are formed when a RG engulfs an Earth-like planet., Comment: 9 pages, 8 figures, accepted 2023 July 10 by MNRAS

Failed supernovae (FSN) are a possible channel for the formation of heavy stellar-mass black holes ($M_{ BH}>\sim 30$ M$_\odot$). However, the effects of metallicity, rotation and magnetic field on the islands of explodabilty of massive stars are not clear. Here, we simulate the stellar structure and evolution in the mass range between 6 and 55 $M_{\odot}$ with different initial rotational velocities, metallicities, and magnetic fields from zero-age main sequence (ZAMS) to pre-collapse. We find that the rapid rotating stars can remain lower $\rm ^{12}C$ mass fraction at the time of C ignition, which allows the transition, from convective carbon burning to radiative burning, to occur at lower $M_{\rm ZAMS}$ than those from stars without rotation. However, the rapid rotation is unfavorable for FSN occurring but is conducive to long gamma-ray bursts (lGRBs) because it results in the specific angular momentum in the CO core is greater than the last stable orbit at core collapse. The increasing metallicity does not affect FSN islands, but high metallicity inhibits rotational mixing and is unfavorable for producing lGRBs. A magnetic field can constrain the mass-loss rate even for rapid rotating stars, resulting in higher mass at pre-collapse. The magnetic braking triggered by the magnetic field can reduce the rotation velocity for high-metallicity models, which decreases the specific angular momentum in the CO core and is favorable for FSN occurring. We suggest that the heavy-mass black holes detected by LIGO may originate from rapidly rotating massive stars with strong magnetic fields, rather than those with very low metallicity., Comment: 16 pages, 11 figures

Observations have shown that the effective temperature of hydrogen-free Wolf-Rayet (WR) stars is considerably lower than that of the standard model, which means that the radius of the observed H-free WR stars is several times larger than that estimated by the standard model. The envelope inflation structure (EIS) caused by the radiation luminosity being close to the Eddington luminosity in the iron opacity peak region of H-free WR stars may be the key to resolve the radius problem of H-free WR stars. We try to explain the H-free WR stars observed in the Milk Way (MW) and the Large Magellanic Cloud (LMC) by the He stars. Using the Modules for Experiments in Stellar Astrophysics code, we compute the evolution of He stars with and without MLT++ prescriptions and discuss their effects on the EIS. We have calculated the evolution of He stars using a new mass-loss rate formula and three different relative rotational velocity and compared our results with observations on Hertzsprung-Russell diagrams. The low luminosity (log$(L/L_{{\odot}})\leq5.2$) H-free WR stars in the MW and the LMC can be explained by the helium giant phase in low-mass He stars, the high $X_{C}$ and $X_{O}$ in WC stars can only evolve through low-mass He stars with a rapid rotation. High-mass He stars with the EIS can explain H-free WR stars with a luminosity exceeding $10^{5.7} L_{{\odot}}$ and an effective temperature above $10^{4.7}$ K in the MW. They can also explain H-free WR stars on the right-hand side of the He zero-age main sequence in the LMC. High-mass stars with the EIS evolve into WO stars at the final evolution stage, and the shorter lifetime fraction is consistent with the small number of observed WO stars., Comment: 9 pages, 7 figures 1 tables, Accepted to A&A

Subdwarf A-type stars (sdAs) are objects that have hydrogen-rich spectra with surface gravity similar to that of hot subdwarf stars but effective temperature below the zero-age horizontal branch (ZAHB). They are considered to be metal-poor main sequence (MS) stars or extremely low-mass white dwarfs (ELM WDs). In this work, using the stellar evolution code, Modules for Experiments in Stellar Astrophysics (MESA), we investigate the sdAs formed both by the evolution of (pre-)ELM WDs in double degenerate systems (DDs) and metal-poor main sequence with the single evolution models. We find that both of the evolutionary tracks of ELM WDs and metal-poor MSs can explain the observation properties of sdAs. However, the proportion between these two populations are uncertain. In this work, we adopt the method of binary population synthesis of both ELM WDs in the disk and metal-poor MS in the halo to obtain the populations of ELM WDs and metal-poor MSs at different stellar population ages and calculate their proportions. We find that the proportion of metal-poor MSs to sdA population for a stellar population of 10 Gyr is $\sim$98.5 percent, which is consistent with the conclusion that most sdAs ($>$95 percent) are metal-poor MSs. And the proportion of ELM WDs (metal-poor MSs) to sdA population increases (decreases) from 0.1\% (99.9\%) to 20\% (80\%) with the stellar population ages from 5 to 13.7 Gyr., Comment: 9 pages, 6 figures, 1 table,Accepter for publication in ApJ

We report the distances of molecular clouds at high Galactic latitudes (|b|>10$^\circ$) derived from parallax and G band extinction (A$_{\rm G}$) measurements in the second Gaia data release, Gaia DR2. Aided by Bayesian analyses, we determined distances by identifying the breakpoint in the extinction A$_{\rm G}$ towards molecular clouds and using the extinction A$_{\rm G}$ of Gaia stars around molecular clouds to confirm the breakpoint. We use nearby star-forming regions, such as Orion, Taurus, Cepheus, and Perseus, whose distances are well-known to examine the reliability of our method. By comparing with previous results, we found that the molecular cloud distances derived from this method are reliable. The systematic error in the distances is approximately 5%. In total, 52 molecular clouds have their distances well determined, most of which are at high Galactic latitudes, and we provide reliable distances for 13 molecular clouds for the first time., Comment: Accepted by A&A. 13 pages, 9 figures

Context. Quasars are spatially stationary, and they are essential objects in astrometry when defining reference frames. However, the census of quasars is far from complete. Mid-infared colors can be used to find quasar candidates because AGNs show a peculiar appearance in mid-infrared color, but these methods are incapable of separating quasars from AGNs. Aims. The aim of our study is to use astrometric and mid-infrared methods to select quasars and get a reliable quasar candidates catalog. Methods. We used a near-zero proper motion criterion in conjuction with WISE (all-sky Wide-field Infrared Survey Explorer) [W1-W2] color to select quasar candidates. The [W1-W2] color criterion is defined by the linear boundary of two samples: LAMOST DR5 quasars, which serve as the quasar sample, and LAMOST DR5 stars and galaxies, which serve as the non-quasar sample. The contamination and completeness are evaluated. Results. We present a catalog of 662 753 quasar candidates, with a completeness of about 75% and a reliability of 77.2%., Comment: 7 pages, 7 figures, accepted by A&A

Gaia data release 2 (DR2) provides the best non-rotating optical frame aligned with the radio frame (ICRF) thanks to the inclusion of about half-million quasars in the 5-parameter astrometric solution. Given their crucial diagnostic role for characterizing the properties of the celestial reference frame, we aim to make an independent assessment of the astrometry of quasars in DR2. We cross-match with Gaia DR2 the quasars from LQAC3, SDSS and LAMOST, obtaining 208743 new sources (denominated as KQCG). With the quasars already identified in DR2, we estimate the global offsets of parallaxes and proper motions of different quasar subsets to check their astrometric consistency. The features of the proper motion field are analyzed by means of vectorial spherical harmonics (VSH); the scalar field of parallaxes is expanded into spherical harmonics to investigate their spatial correlation. We find a bias of $\sim$ $-0.030$ $mas$ in KQCG parallaxes, and a bias of $-9.1$ $\mu as/yr$ in $\mu_{\alpha\ast}$; a $\sim$ +10 $\mu as/yr$ bias in $\mu_{\delta}$ consistently found in the entire quasar sample. The results of the VSH analysis on different subsets indicate good agreement between them. The proper motion field exhibits a ((-6,-5,-5) $\pm$ 1) $\mu as/yr$ rotation in the northern hemisphere and a rotation of ((0,+1,+3) $\pm$ 1 ) $\mu as /yr$ in the southern one. The spherical harmonics expansion of the parallax field reveals an angular scale of systematics of $\approx$ 18 degrees with a RMS amplitude of 13 $\mu$as. The positional comparison shows that the axes of the current Gaia Celestial Reference Frame and the ICRF2 are aligned within 30 $\mu$as.