Search

The Taurus region is one of the most extensively studied star-forming regions. Surveys indicate that the young stars in this region are comprised of Young Stellar Objects (YSOs) that cluster in groups associated with the molecular cloud (Grouped Young Stellar Objects, GYSOs), and some older ones that are sparsely distributed throughout the region (Distributed Young Stellar Objects, DYSOs). To bridge the age gap between the GYSOs ($\le$5 Myr) and the DYSOs (10-20 Myr), we conducted a survey to search for new YSOs in this direction. Based on infrared color excesses and Li I absorption lines, we identified 145 new YSOs. Combining these with the previously identified GYSOs and DYSOs, we constructed a sample of 519 YSOs that encompasses the entire region. Subsequently, we calculated the ages of the sample based on their proximity to the local bubble. The age versus Distance to the Local Bubble ($D_{\rm LB}$) relationship for the DYSOs shows a clear trend: the farther they are from the Local Bubble, the younger they are, which is consistent with the supernovae-driven formation scenario of the Local Bubble. The GYSOs also exhibit a mild age versus $D_{\rm LB}$ trend. However, they are significantly younger and are mostly confined to distances of 120 to 220 pc. Considering their distribution in the age versus $D_{\rm LB}$ space is well separated from the older DYSOs, they may also be products of the Local Bubble but formed in more recent and localized events., Comment: 10 pages, Accepted for publication in SCPMA

Based on the LAMOST spectroscopy and TESS time-series photometry, we have obtained a main-sequence star sample of $\delta$ Scuti and $\gamma$ Doradus stars. The sample includes 1534 $\delta$ Sct stars, 367 $\gamma$ Dor stars, 1703 $\delta$ Sct$| \gamma$ Dor stars, 270 $\gamma$ Dor$| \delta$ Sct stars, along with 105 '$\delta$ Sct candidates' and 32 '$\gamma$ Dor candidates'. After correcting for projection effects, we derived the equatorial rotational velocity distribution for $\delta$ Sct and $\gamma$ Dor stars and compared it with that of normal stars. The rotational velocity distributions of $\delta$ Sct and $\gamma$ Dor stars are extremely similar, with the only difference potentially due to the rotational variable stars that have not been completely removed. In contrast, the rotational velocity distribution of normal stars is more dispersed compared to pulsating stars. Additionally, the peak rotational velocity of the pulsating stars is about 10 km s$^{-1}$ higher than that of normal stars. Unlike the normal stars, which show a monotonic increase in peak velocity with mass between 1.8 and 2.5 $M_{\odot}$, the rotational velocity distribution of $\delta$ Sct stars does not exhibit a strong mass dependence. We also found that normal stars accelerate during the late main-sequence evolutionary phase, while $\delta$ Sct stars decelerate. Furthermore, there may still be unclassified stars with diverse rotational properties in the normal star sample compared to the $\delta$ Sct stars, which is likely to be an important contributor to the broader dispersion observed in its rotational velocity distribution. The photometric amplitude in $\delta$ Sct stars is modulated with rotational velocity, with high-amplitude stars typically rotating slowly and low-amplitude stars showing a broad distribution of rotational velocities., Comment: 13 pages, 8 figures, 1 table, APJ accepted

LEIA (Lobster Eye Imager for Astronomy) detected a new X-ray transient on November 7, 2022, identified as a superflare event occurring on a nearby RS CVn-type binary HD 251108. The flux increase was also detected in follow-up observations at X-ray, UV and optical wavelengths. The flare lasted for about 40 days in soft X-ray observations, reaching a peak luminosity of ~1.1 * 10^34 erg/s in 0.5-4.0 keV, which is roughly 60 times the quiescent luminosity. Optical brightening was observed for only one night. The X-ray light curve is well described by a double "FRED" (fast rise and exponential decay) model, attributed to the cooling process of a loop arcade structure formed subsequent to the initial large loop with a half-length of ~1.9 times the radius of the host star. Time-resolved X-ray spectra were fitted with a two-temperature apec model, showing significant evolution of plasma temperature, emission measure, and metal abundance over time. The estimated energy released in the LEIA band is ~3 * 10^39 erg, suggesting this is likely the most energetic X-ray stellar flare with the longest duration detected to date., Comment: submitted to ApJL, 22 pages, 9 figures, 7 tables

The cloud cover and meteorological parameters serve as fundamental criteria for the qualification of an astronomical observatory working in optical and infrared wavelengths. In this paper, we present a systematic assessment of key meteorological parameters at the Lenghu site. The datasets adopted in this study includes the meteorological parameters collected at the local weather stations at the site and in the Lenghu Town, the sky brightness at the local zenith acquired by the Sky Quality Meters and night sky all-sky images from a digital camera, the ERA5 reanalysis database and global climate monitoring data. From 2019 to 2023, the fractional observable time of photometric condition is 69.70%, 74.97%, 70.26%, 74.27% and 65.12%, respectively. The fractional observing time is inversely correlated with surface air temperature, relative humidity, precipitable water vapor, and dew temperature, demonstrating that the observing conditions are influenced by these meteorological parameters. Large-scale air-sea interactions affect the climate at Lenghu site, which in fact delivers a clue to understand the irregularity of 2023. Specifically, precipitable water vapor at Lenghu site is correlated to both the westerly wind index and the summer North Atlantic Oscillation index, the yearly average temperature of Lenghu site is observed to increase significantly during the occurrence of a strong El Ni\~no event and the relative humidity anomaly at Lenghu site is correlated to the Pacific Decadal Oscillation index. The decrease of fractional observing time in 2023 was due to the ongoing strong El Ni\~no event and relevant global climate change. We underscore the substantial role of global climate change in regulating astronomical observing conditions and the necessity for long-term continuous monitoring of the astronomical meteorological parameters at Lenghu site., Comment: accepted for publication in MNRAS

We obtain the largest catalog of multi-mode $\delta$ Sct stars in the northern sky to date using the Zwicky Transient Facility (ZTF) Data Release 20 (DR20). The catalog includes 2254 objects, of which 2181 are new to our study. Among these multi-mode $\delta$ Sct stars, 2142 objects are double-mode $\delta$ Sct, while 109 objects are triple-mode $\delta$ Sct and 3 are quadruple-mode $\delta$ Sct. By analyzing the light curves in the $r$ and $g$ bands of the ZTF, we determine the basic parameters of multi-mode $\delta$ Sct stars, including the periods and amplitudes. Periods are checked by comparison with the OGLE catalog of double-mode $\delta$ Sct stars. On the Petersen diagram, multi-mode $\delta$ Sct stars form six sequences. We find that in Galactic coordinates, the periods of 1O/F double-mode $\delta$ Sct stars at high latitudes are shorter than those of 1O/F double-mode $\delta$ Sct stars in the disk, due to metallicity variations. In the future, our catalog can be used to establish the period--luminosity relation and the period--metallicity relation of double-mode $\delta$ Sct stars, and to study the Galactic structure., Comment: 18 pages, 7 figures, accepted for publication in ApJS

The presence of young stellar populations in the Large Magellanic Cloud cluster NGC 1783 has caught significant attention, with suggestions ranging from it being a genuine secondary stellar generation to a population of blue straggler stars or simply contamination from background stars. Thanks to multi-epoch observations with the Hubble Space Telescope, proper motions for stars within the field of NGC 1783 have been derived, thus allowing accurate cluster membership determination. Here, we report that the younger stars within NGC 1783 indeed belong to the cluster, and their spatial distribution is more extended compared to the bulk of the older stellar population, consistent with previous studies. Through N-body simulations, we demonstrate that the observed characteristics of the younger stars cannot be explained solely by blue straggler stars in the context of the isolated dynamical evolution of NGC 1783. Instead, accretion of the external, low-mass stellar system can better account for both the inverse spatial concentration and the radial velocity isotropy of the younger stars. We propose that NGC 1783 may have accreted external stars from low-mass stellar systems, resulting in a mixture of external younger stars and blue straggler stars from the older bulk population, thereby accounting for the characteristics of the younger sequence., Comment: 11 pages, 6 figures; accepted for publication in ApJ

Classical Cepheids (CCs) are excellent tracers for understanding the structure of the Milky Way disk. The latest Gaia Data Release 3 provides a large number of line-of-sight velocity information for Galactic CCs, offering an opportunity for studying the kinematics of the Milky Way. We determine the three-dimensional velocities of 2057 CCs relative to the Galactic center. From the projections of the 3D velocities onto the XY plane of the Galactic disk, we find that $V_{R}$ and $V_{\phi}$ velocities of the northern and southern warp (directions with highest amplitude) are different. This phenomenon may be related to the warp precession or the asymmetry of the warp structure. By investigating the kinematic warp model, we find that the vertical velocity of CCs is more suitable for constraining the warp precession rate than the line of nodes angles. Our results suggest that CCs at $12-14$ kpc are the best sample for determining the Galactic warp precession rate. Based on the spatial structure parameters of Cepheid warp from Chen et al (arXiv:1902.00998), we determine a warp precession rate of $\omega = 4.9\pm1.6$ km s$^{-1}$ kpc$^{-1}$ at 13 kpc, which supports a low precession rate in the warp model. In the future, more kinematic information on CCs will help to better constrain the structure and evolution of the Milky Way., Comment: 15 pages, 8 figures, accepted for publication in ApJ

RR Lyrae (RR Lyr) stars are a well-known and useful distance indicator for old stellar populations such as globular clusters and dwarf galaxies. Fundamental-mode RR Lyr (RRab) stars are commonly used to measure distances, and the accuracy of the determined distance is strongly constrained by metallicity. Here, we investigate the metallicity dependence in the period-luminosity (PL) relation of double-mode RR Lyr (RRd) stars. We find and establish a linear relation between metallicity and period or period ratio for RRd stars. This relation can predict the metallicity as accurately as the low-resolution spectra. Based on this relation, we establish a metallicity-independent PL relation for RRd stars. Combining the distance of the Large Magellanic Cloud and Gaia parallaxes, we calibrate the zero point of the derived PL relation to an error of 0.022 mag. Using RRd stars, we measure the distances of globular clusters and dwarf galaxies with an accuracy of 2-3% and 1-2%, respectively. In the future, RRd stars could anchor galaxy distances to an accuracy of 1.0% and become an independent distance ladder in the Local Group., Comment: 6 pages, 2 figures, to appear on the Proceedings of IAU Symposium 376, At the cross-roads of astrophysics and cosmology: Period-luminosity relations in the 2020s

RR Lyrae stars are one of the primary distance indicators for old stellar populations such as globular clusters, dwarf galaxies and galaxies. Typically, fundamental-mode RR Lyr stars are used for distance measurements, and their accuracy is strongly limited by the dependence of absolute magnitudes on metallicity, in both the optical and infrared bands. Here, we report the discovery of a period-(period ratio)-metallicity relation for double-mode RR Lyr stars, which can predict metallicity as accurately as the low-resolution spectra. With theoretical and observational evidence, we propose that the period-luminosity relation of double-mode RR Lyr stars is not affected by the metallicity. Combining the Large Magellanic Cloud distance and Gaia parallaxes, we calibrate the zero point of the period-luminosity relation to an error of 0.022 mag, which means that in the best case double-mode RR Lyr stars can anchor galaxy distances to an accuracy of 1.0%. For four globular clusters and two dwarf galaxies, we obtain distances using double-mode RR Lyr stars with a distance accuracy of 2-3% and 1-2%, respectively. With future telescopes such as the China Space Station Telescope and the Vera C. Rubin Observatory, double-mode RR Lyr stars will be established as an independent distance ladder in the near-field universe., Comment: 19 pages, 3 figures

Superluminous supernovae are among the most energetic stellar explosions in the Universe, but their energy sources remain an open question. Here we present long-term observations of one of the closest examples of the hydrogen-poor subclass (SLSNe-I), SN~2017egm, revealing the most complicated known luminosity evolution of SLSNe-I. Three distinct post-peak bumps were recorded in its light curve collected at about $100$--350\,days after maximum brightness, challenging current popular power models such as magnetar, fallback accretion, and interaction between ejecta and a circumstellar shell. However, the complex light curve can be well modelled by successive interactions with multiple circumstellar shells with a total mass of about $6.8$--7.7\,M$_\odot$. In this scenario, large energy deposition from interaction-induced reverse shocks results in ionization of neutral oxygen in the supernova ejecta and hence a much lower nebular-phase line ratio of [O\,\textsc{i}] $\lambda6300$/([Ca\,\textsc{ii}] + [O\,\textsc{ii}]) $\lambda7300$ ($\sim 0.2$) compared with that derived for other superluminous and normal stripped-envelope SNe. The pre-existing multiple shells indicate that the progenitor of SN~2017egm experienced pulsational mass ejections triggered by pair instability within 2 years before explosion, in robust agreement with theoretical predictions for a pre-pulsation helium-core mass of 48--51\,M$_{\odot}$. Finally, this work shows that the final explosion product may be a black hole with about 40\,M$_{\odot}$, and has significant implication for the formation of such heavy black holes that have been recently observed by LIGO-Virgo gravitational wave detectors., Comment: 34 pages, 13 figures, 6 tables

The Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST), also known as the Guoshoujing Telescope, is a major national scientific facility for astronomical research located in Xinglong, China. Beginning with a pilot survey in 2011, LAMOST has been surveying the night sky for more than 10 years. The LAMOST survey covers various objects in the Universe, from normal stars to peculiar ones, from the Milky Way to other galaxies, and from stellar black holes and their companions to quasars that ignite ancient galaxies. Until the latest data release 8, the LAMOST survey has released spectra for more than 10 million stars, ~220,000 galaxies, and ~71,000 quasars. With this largest celestial spectra database ever constructed, LAMOST has helped astronomers to deepen their understanding of the Universe, especially for our Milky Way galaxy and the millions of stars within it. In this article, we briefly review the characteristics, observations, and scientific achievements of LAMOST. In particular, we show how astrophysical knowledge about the Milky Way has been improved by LAMOST data., Comment: 16 pages, 10 figures. Published by The Innovation

The development of large-scale time-domain surveys provides an opportunity to study the physical properties as well as the evolutionary scenario of B-type subdwarfs (sdB) and M-type dwarfs (dM). Here, we obtained 33 sdB+dM eclipsing binaries based on the Zwicky Transient Facility (ZTF) light curves and $Gaia$ early data release 3 (EDR3) parallaxes. By using the PHOEBE code for light curve analysis, we obtain probability distributions for parameters of 29 sdB+dM. $R_1$, $R_2$, and $i$ are well determined, and the average uncertainty of mass ratio $q$ is 0.08. Our parameters are in good agreement with previous works if a typical mass of sdB is assumed. Based on parameters of 29 sdB+dM, we find that both the mass ratio $q$ and the companion's radius $R_2$ decrease with the shortening of the orbital period. For the three sdB+dMs with orbital periods less than 0.075 days, their companions are all brown dwarfs. The masses and radii of the companions satisfy the mass--radius relation for low-mass stars and brown dwarfs. Companions with radii between $0.12R_\odot$ and $0.15R_\odot$ seem to be missing in the observations. As more short-period sdB+dM eclipsing binaries are discovered and classified in the future with ZTF and $Gaia$, we will have more information to constrain the evolutionary ending of sdB+dM., Comment: 16 pages, 10 figures, 5 tables, 2 appendices (14 pages, 27 figures), accepted for publication in Research in Astronomy and Astrophysics

The behavior of the shock wave in the atmosphere of the non-fundamental mode RR Lyrae pulsator remains a mystery. In this work, we firstly report a blueshifted Mg triplet emission in continuous spectroscopic observations for a non-Blazhko RRc pulsator (Catalina-1104058050978) with LAMOST medium resolution spectra. We analyse the photometric observations from Catalina Sky Survey of this RRc pulsator with pre-whitening sequence method and provide the ephemeris and phases. An additional frequency signal with $P_1/P_x = 0.69841$ is detected and discussed. The redshift and radial velocity of the spectra are provided by fitting process with $S\acute{e}rsic$ functions and cross-correlation method. Moreover, we plot the variation of H$\alpha$ and Mg lines in a system comoving with the pulsation. Clear evolution of comoving blueshifted hydrogen and Mg emission is observed, which further confirms the existence of shock waves in RRc pulsators. The shock-triggered emission lasts over $15\%$ of the pulsation cycle, which is much longer than the previous observations.

Angular momentum is a key property regulating star formation and evolution. However, the physics driving the distribution of the stellar rotation rates of early-type main-sequence stars is as yet poorly understood. Using our catalog of 40,034 early-type stars with homogeneous $v\sin i$ parameters, we review the statistical properties of their stellar rotation rates. We discuss the importance of possible contaminants, including binaries and chemically peculiar stars. Upon correction for projection effects and rectification of the error distribution, we derive the distributions of our sample's equatorial rotation velocities, which show a clear dependence on stellar mass. Stars with masses less than $2.5\ {M_\odot}$ exhibit a unimodal distribution, with the peak velocity ratio increasing as stellar mass increases. A bimodal rotation distribution, composed of two branches of slowly and rapidly rotating stars, emerges for more massive stars ($M>2.5\ {M_\odot}$). For stars more massive than $3.0\ {M_\odot}$, the gap between the bifurcated branches becomes prominent. For the first time, we find that metal-poor ([M/H] $< -0.2$ dex) stars only exhibit a single branch of slow rotators, while metal-rich ([M/H] $> 0.2$ dex) stars clearly show two branches. The difference could be attributed to unexpectedly high spin-down rates and/or in part strong magnetic fields in the metal-poor subsample., Comment: 18 pages, 9 figures. Accepted for publication in ApJ; revised references for section 3.2

We derive stellar parameters and abundances (`stellar labels') of 40,034 late-B and A-type main-sequence stars extracted from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope Medium Resolution Survey (LAMOST--MRS). The primary selection of our early-type sample was obtained from LAMOST Data Release 7 based on spectral line indices. We employed the Stellar LAbel Machine (SLAM) to derive their spectroscopic stellar parameters, drawing on Kurucz spectral synthesis models with 6000 K $< T_\mathrm{eff} <$ 15,000 K and $-1$ dex $< \mathrm{[M/H]} <$ 1 dex. For a signal-to-noise ratio of $\sim 60$, the cross-validated scatter is $\sim$75 K, 0.06 dex, 0.05 dex, and $\sim 3.5\,\mathrm{km\,s^{-1}}$ for $T_\mathrm{eff}$, $\log g$, [M/H], and $v\sin i$, respectively. A comparison with objects with prior, known stellar labels shows great consistency for all stellar parameters, except for $\log g$. Although this is an intrinsic caveat that comes from the MRS's narrow wavelength coverage, it only has a minor effect on estimates of the stellar rotation rates because of the decent spectral resolution and the profile-fitting method employed. The masses and ages of our early-type sample stars were inferred from non-rotating stellar evolution models. This paves the way for reviewing the properties of stellar rotation distributions as a function of stellar mass and age., Comment: 18 pages, 10 figures, 1 table. Accepted for publication in ApJS

Hydrogen emissions of RR Lyrae variables are the imprints of shock waves traveling through their atmospheres. We develop a pattern recognition algorithm, which is then applied to single-epoch spectra of SDSS and LAMOST. These two spectroscopic surveys covered $\sim$ 10,000 photometrically confirmed RR Lyrae stars. We discovered in total 127 RR Lyrae stars with blueshifted Balmer emission feature, including 103 fundamental mode (RRab), 20 first-overtone (RRc), 3 double-mode (RRd), and 1 Blazhko type (temporary classification for RR Lyrae stars with strong Blazhko modulation in Catalina sky survey that cannot be characterized) RR Lyrae variable. This forms the largest database to date of the properties of hydrogen emission in RR Lyrae variables. Based on ZTF DR5, we carried out a detailed light-curve analysis for the Blazhko type RR Lyrae star with hydrogen emission of long-term modulations. We characterize the Blazhko type RR Lyrae star as an RRab and point out a possible Blazhko period. Finally, we set up simulations on mock spectra to test the performance of our algorithm and on the real observational strategy to investigate the occurrence of the "first apparition"., Comment: ApJ, published. 20 pages, 14 figures, 2 tables; The final version

We independently determine the zero-point offset of the Gaia early Data Release-3 (EDR3) parallaxes based on $\sim 110,000$ W Ursae Majoris (EW)-type eclipsing binary systems. EWs cover almost the entire sky and are characterized by a relatively complete coverage in magnitude and color. They are an excellent proxy for Galactic main-sequence stars. We derive a $W1$-band Period-Luminosity relation with a distance accuracy of $7.4\%$, which we use to anchor the Gaia parallax zero-point. The final, global parallax offsets are $-28.6\pm0.6$ $\mu$as and $-25.4\pm4.0$ $\mu$as (before correction) and $4.2\pm0.5$ $\mu$as and $4.6\pm3.7$ $\mu$as (after correction) for the five- and six-parameter solutions, respectively. The total systematic uncertainty is $1.8$ $\mu$as. The spatial distribution of the parallax offsets shows that the bias in the corrected Gaia EDR3 parallaxes is less than 10 $\mu$as across $40\%$ of the sky. Only $15\%$ of the sky is characterized by a parallax offset greater than 30 $\mu$as. Thus, we have provided independent evidence that the parallax zero-point correction provided by the Gaia team significantly reduces the prevailing bias. Combined with literature data, we find that the overall Gaia EDR3 parallax offsets for Galactic stars are $[-20, -30]$ $\mu$as and 4-10 $\mu$as, respectively, before and after correction. For specific regions, an additional deviation of about 10 $\mu$as is found., Comment: 7 pages, 4 figures. Accepted for publication in The Astrophysical Journal Letters

The color-magnitude diagrams (CMDs) of young star clusters show that, particularly at ultraviolet wavelengths, their upper main sequences (MSs) bifurcate into a sequence comprising the bulk population and a blue periphery. The spatial distribution of stars is crucial to understand the reasons for these distinct stellar populations. This study uses high-resolution photometric data obtained with the Hubble Space Telescope to study the spatial distributions of the stellar populations in seven Magellanic Cloud star clusters. The cumulative radial number fractions of blue stars within four clusters are strongly anti-correlated with those of the high-mass-ratio binaries in the bifurcated region, with negative Pearson coefficients < -0.7. Those clusters generally are young or in an early dynamical evolutionary stage. In addition, a supporting N-body simulation suggests the increasing percentage of blue-MS stars from the cluster centers to their outskirts may be associated with the dissolution of soft binaries. This study provides a different perspective to explore the MS bimodalities in young clusters and adds extra puzzles. A more comprehensive study combined with detailed simulations is needed in the future., Comment: 19 pages, 9 figures, 3 tables. Accepted for publication in ApJ

The impact of stellar rotation on the morphology of star cluster colour-magnitude diagrams is widely acknowledged. However, the physics driving the distribution of the equatorial rotation velocities of main-sequence turn-off (MSTO) stars is as yet poorly understood. Using Gaia Data Release 2 photometry and new Southern African Large Telescope medium-resolution spectroscopy, we analyse the intermediate-age ($\sim1\,$Gyr-old) Galactic open clusters NGC 3960, NGC 6134 and IC 4756 and develop a novel method to derive their stellar rotation distributions based on SYCLIST stellar rotation models. Combined with literature data for the open clusters NGC 5822 and NGC 2818, we find a tight correlation between the number ratio of slow rotators and the clusters' binary fractions. The blue-main-sequence stars in at least two of our clusters are more centrally concentrated than their red-main-sequence counterparts. The origin of the equatorial stellar rotation distribution and its evolution remains as yet unidentified. However, the observed correlation in our open cluster sample suggests a binary-driven formation mechanism., Comment: 9 pages, 7 figures, 1 table. Accepted for publication in MNRAS

We present a catalog of eclipsing binaries in the northern Galactic Plane from the Kiso Wide-Field Camera Intensive Survey of the Galactic Plane (KISOGP). We visually identified 7055 eclipsing binaries spread across $\sim$330 square degrees, including 4197 W Ursa Majoris/EW-, 1458 $\beta$ Lyrae/EB-, and 1400 Algol/EA-type eclipsing binaries. For all systems, $I$-band light curves were used to obtain accurate system parameters. We derived the distances and extinction values for the EW-type objects from their period--luminosity relation. We also obtained the structure of the thin disk from the distribution of our sample of eclipsing binary systems, combined with those of high-mass star-forming regions and Cepheid tracers. We found that the thin disk is inhomogeneous in number density as a function of Galactic longitude. Using this new set of distance tracers, we constrain the detailed structure of the thin disk. Finally, we report a global parallax zero-point offset of $ \Delta \pi=-42.1\pm1.9\mbox{(stat.)}\pm12.9\mbox{(syst.)}$ $\mu$as between our carefully calibrated EW-type eclipsing binary positions and those provided by Gaia Early Data Release 3. Implementation of the officially recommended parallax zero-point correction results in a significantly reduced offset. Additionally, we provide a photometric characterization of our EW-type eclipsing binaries that can be applied to further analyses., Comment: 17 pages, 9 figures, 2 tables. Accepted for publication in AJ

The study of planet occurrence as a function of stellar mass is important for a better understanding of planet formation. Estimating stellar mass, especially in the red giant regime, is difficult. In particular, stellar masses of a sample of evolved planet-hosting stars based on spectroscopy and grid-based modelling have been put to question over the past decade with claims they were overestimated. Although efforts have been made in the past to reconcile this dispute using asteroseismology, results were inconclusive. In an attempt to resolve this controversy, we study four more evolved planet-hosting stars in this paper using asteroseismology, and we revisit previous results to make an informed study of the whole ensemble in a self-consistent way. For the four new stars, we measure their masses by locating their characteristic oscillation frequency, $\mathrm{\nu}_{\mathrm{max}}$, from their radial velocity time series observed by SONG. For two stars, we are also able to measure the large frequency separation, $\mathrm{\Delta\nu}$, helped by extended SONG single-site and dual-site observations and new TESS observations. We establish the robustness of the $\mathrm{\nu}_{\mathrm{max}}$-only-based results by determining the stellar mass from $\mathrm{\Delta\nu}$, and from both $\mathrm{\Delta\nu}$ and $\mathrm{\nu}_{\mathrm{max}}$. We then compare the seismic masses of the full ensemble of 16 stars with the spectroscopic masses from three different literature sources. We find an offset between the seismic and spectroscopic mass scales that is mass-dependent, suggesting that the previously claimed overestimation of spectroscopic masses only affects stars more massive than about 1.6 M$_\mathrm{\odot}$., Comment: 13 pages, 12 figures, Accepted for publication in the Main Journal of MNRAS

The number of known periodic variables has grown rapidly in recent years. Thanks to its large field of view and faint limiting magnitude, the Zwicky Transient Facility (ZTF) offers a unique opportunity to detect variable stars in the northern sky. Here, we exploit ZTF Data Release 2 (DR2) to search for and classify variables down to r ~ 20.6 mag. We classify 781,602 periodic variables into 11 main types using an improved classification method. Comparison with previously published catalogs shows that 621,702 objects (79.5%) are newly discovered or newly classified, including ~700 Cepheids, ~5000 RR Lyrae stars, ~15,000 Delta Scuti variables, ~350,000 eclipsing binaries, ~100,000 long-period variables, and about 150,000 rotational variables. The typical misclassification rate and period accuracy are on the order of 2% and 99%, respectively. 74% of our variables are located at Galactic latitudes, $|b|<10^\circ$. This large sample of Cepheids, RR Lyrae, Delta Scuti stars, and contact (EW-type) eclipsing binaries is helpful to investigate the Galaxy's disk structure and evolution with an improved completeness, areal coverage, and age resolution. Specifically, the northern warp and the disk's edge at distances of 15--20 kpc are significantly better covered than previously. Among rotational variables, RS Canum Venaticorum and BY Draconis-type variables can be separated easily. Our knowledge of stellar chromospheric activity would benefit greatly from a statistical analysis of these types of variables., Comment: accepted by ApJS. The full catalog and LCs can be accessed at http://variables.cn:88/ztf/ or doi.org/10.5281/zenodo.3886372

We present the physical parameters of 2335 late-type contact binary (CB) systems extracted from the Catalina Sky Survey (CSS). Our sample was selected from the CSS Data Release 1 by strictly limiting the prevailing temperature uncertainties and light-curve fitting residuals, allowing us to almost eliminate any possible contaminants. We developed an automatic Wilson--Devinney-type code to derive the relative properties of CBs based on their light-curve morphology. By adopting the distances derived from CB (orbital) period--luminosity relations (PLRs), combined with the well-defined mass--luminosity relation for the systems' primary stars and assuming solar metallicity, we calculated the objects' masses, radii, and luminosities. Our sample of fully eclipsing CBs contains 1530 W-, 710 A-, and 95 B-type CBs. A comparison with literature data and with the results from different surveys confirms the accuracy and coherence of our measurements. The period distributions of the various CB subtypes are different, hinting at a possible evolutionary sequence. W-type CBs are clearly located in a strip in the total mass versus mass ratio plane, while A-type CBs may exhibit a slightly different dependence. There are no significant differences among the PLRs of A- and W-type CBs, but the PLR zero points are affected by their mass ratios and fill-out factors. Determination of zero-point differences for different types of CBs may help us improve the accuracy of the resulting PLRs. We demonstrate that automated approaches to deriving CB properties could be a powerful tool for application to the much larger CB samples expected to result from future surveys., Comment: 2o pages, 10 figures, 2 tables. Accepted for publication in ApJS

We investigate the influence of Wolf-Rayet (W-R) stars on their surrounding star-forming molecular clouds. We study five regions containing W-R stars in the inner Galactic plane ($l\sim$[14$^\circ$-52$^\circ$]), using multi-wavelength data from near-infrared to radio wavelengths. Analysis of $^{13}$CO line data reveals that these W-R stars have developed gas-deficient cavities in addition to molecular shells with expansion velocities of a few km s$^{-1}$. The pressure owing to stellar winds primarily drives these expanding shells and sweeps up the surrounding matter to distances of a few pc. The column densities of shells are enhanced by a minimum of 14% for one region to a maximum of 88% for another region with respect to the column densities within their central cavities. No active star formation - including molecular condensations, protostars, or ionized gas - is found inside the cavities, whereas such features are observed around the molecular shells. Although the expansion of ionized gas is considered an effective mechanism to trigger star formation, the dynamical ages of the HII regions in our sample are generally not sufficiently long to do so efficiently. Overall, our results hint at the possible importance of negative W-R wind-driven feedback on the gas-deficient cavities, where star formation is quenched as a consequence. In addition, the presence of active star formation around the molecular shells indicates that W-R stars may also assist in accumulating molecular gas, and that they could initiate star formation around those shells., Comment: 24 pages, 15 figures (4 additional figures in Appendix), 3 tables, Accepted for publication in the Astrophysical Journal

Stars spend most of their lifetimes on the `main sequence' (MS) in the Hertzsprung--Russell diagram. The obvious double MSs seen in the equivalent color--magnitude diagrams characteristic of Milky Way open clusters pose a fundamental challenge to our traditional understanding of star clusters as `single stellar populations.' The clear MS bifurcation of early-type stars with masses greater than $\sim1.6 M_\odot$ is thought to result from a range in the stellar rotation rates. However, direct evidence connecting double MSs to stellar rotation properties has yet to emerge. Here, we show through analysis of the projected stellar rotational velocities ($v\sin i$, where $i$ represents the star's inclination angle) that the well-separated double MS in the young, $\sim200Myr$-old Milky Way open cluster NGC 2287 is tightly correlated with a dichotomous distribution of stellar rotation rates. We discuss whether our observations may reflect the effects of tidal locking affecting a fraction of the cluster's member stars in stellar binary systems. We show that the slow rotators could potentially be initially rapidly rotating stars that have been slowed down by tidal locking by a low mass-ratio companion in a cluster containing a large fraction of short-period, low-mass-ratio binaries. This demonstrates that stellar rotation drives the split MSs in young, $\lessapprox 300$Myr-old star clusters. However, special conditions, e.g., as regards the mass-ratio distribution, might be required for this scenario to hold., Comment: 15 pages, 8 figures, 1 table. Accepted for publication in ApJ

We present a detailed analysis of the projected stellar rotational velocities of the well-separated double main sequence (MS) in the young, $\sim200$Myr-old Milky Way open cluster NGC 2287 and suggest that stellar rotation may drive the split MSs in NGC 2287. We find that the observed distribution of projected stellar rotation velocities could result from a dichotomous distribution of stellar rotation rates. We discuss whether our observations may reflect the effects of tidal locking affecting a fraction of the cluster's member stars in stellar binary systems. The slow rotators are likely stars that initially rotated rapidly but subsequently slowed down through tidal locking induced by low-mass-ratio binary systems. However, the cluster may have a much larger population of short-period binaries than is usually seen in the literature, with relatively low secondary masses., Comment: 4 pages, 1 figure. To appear in the Proceedings of the 351 IAU Symposium "Star Clusters: From the Milky Way to the Early Universe"

Combing Gaia DR2 with LAMOST DR5, we spectroscopically identified 924 hot subdwarf stars, among which 32 stars exhibit strong double-lined composite spectra. We measured the effective temperature $T_{\rm eff}$, surface gravity $\log\,g$, helium abundance $y=n{\rm He}/n{\rm H}$, and radial velocities of 892 non-composite spectra hot subdwarf stars by fitting LAMOST observations with Tlusty/Synspec non-LTE synthetic spectra. We outlined four different groups in the $T_{\rm eff}-\log\,g$ diagram with our helium abundance classification scheme and two nearly parallel sequences in the $T_{\rm eff}-\log(y)$ diagram. 3D Galactic space motions and orbits of 747 hot subdwarf stars with $(G_{BP}-G_{RP})_{0}<-0.36$ mag were computed using LAMOST radial velocities and Gaia parallaxes and proper motions. Based on the $U-V$ velocity diagram, $J_{z}-$eccentricity diagram, and Galactic orbits, we derived Galactic population classifications and the fractional distributions of the four hot subdwarf helium groups in the halo, thin disk and thick disk. Comparisons with the predictions of binary population synthesis calculations (Han 2008) suggest that He-rich hot subdwarf stars with $\log(y)\ge0$ are from the double helium white dwarfs merger, He-deficient hot subdwarf stars with $-2.2\le\log(y)<-1$ from the common envelope ejection, and He-deficient hot subdwarf stars with $\log(y)<-2.2$ from the stable Roche lobe overflow channels. The relative number of He-rich hot subdwarf stars with $-1\le\log(y)<0$ and $\log(y)\ge0$ in the halo is more than twice the prediction of Zhang et al.(2017), even more than six times in the thin disk, which implies that the mergers of helium white dwarfs with low mass main sequence stars may not be the main formation channel of He-rich hot subdwarf stars with $-1\le\log(y)<0$, specially in younger environments., Comment: Accepted for publication in the Astrophysical Journal

The origin of extended main sequence turnoffs (eMSTOs) in intermediate-age (1-3 Gyr) clusters is one of the most intriguing questions in current star cluster research. Unlike the split main sequences found in some globular clusters. which are caused by bimodal populations in age and/or chemical abundances, eMSTOs are believed to be owing to stellar rotation. We present a spectroscopic survey of MSTO stars in a nearby, intermediate-age (\unit[0.9]{Gyr}), low mass ($\sim 1.7\times10^3 M_\odot$) Galactic open cluster, NGC 5822. We derive a clean sample of member stars based on \textit{Gaia} proper motions and parallaxes and confirm the existence of an eMSTO. Using medium-resolution ($R\sim 4000$) Southern African Large Telescope (SALT) spectra, we derive the rotational velocities of 24 member stars (representing 20\% completeness around the eMSTO region) and find that the loci of the main sequence stars in the eMSTO region show a clear correlation with the projected rotational velocities in the sense that fast rotators are located on the red side of the eMSTO and slow rotators are found on the blue side. By comparison with a synthetic cluster model, we show that the stellar rotational velocities and the eMSTO of NGC 5822 can be well reproduced and we conclude that stellar rotation is the main cause of the eMSTO in NGC 5822., Comment: 13 pages, 7 figures, 1 table. Accepted for publication in ApJ

Using {\sl Gaia} Data Release 2 photometry, we report the detection of extended main-sequence turnoff (eMSTO) regions in the color--magnitude diagrams (CMDs) of the $\sim 14$ Myr-old double clusters $h$ and $\chi$ Persei (NGC 869 and NGC 884). We find that stars with masses below $\sim$1.3 $M_{\odot}$ in both $h$ and $\chi$ Persei populate narrow main sequences (MSs), while more massive stars define the eMSTO, closely mimicking observations of young Galactic and Magellanic Cloud clusters (with ages older than $\sim$30 Myr). Previous studies based on clusters older than $\sim$30 Myr find that rapidly rotating MS stars are redder than slow rotators of similar luminosity, suggesting that stellar rotation may be the main driver of the eMSTO. By combining photometry and projected rotational velocities from the literature of stars in $h$ and $\chi$ Persei, we find no obvious relation between the rotational velocities and colors of non-emission-line eMSTO stars, in contrast with what is observed in older clusters. Similarly to what is observed in Magellanic Cloud clusters, most of the extremely rapidly rotating stars, identified by their strong H$\alpha$ emission lines, are located in the red part of the eMSTOs. This indicates that stellar rotation plays a role in the color and magnitude distribution of MSTO stars. By comparing the observations with simulated CMDs, we find that a simple population composed of coeval stars that span a wide range of rotation rates is unable to reproduce the color spread of the clusters' MSs. We suggest that variable stars, binary interactions, and stellar rotation affect the eMSTO morphology of these very young clusters., Comment: 14 pages, 12 figures, ApJ accepted

The Milky Way's neutral hydrogen (HI) disk is warped and flared. However, a dearth of accurate HI-based distances has thus far prevented the development of an accurate Galactic disk model. Moreover, the extent to which our Galaxy's stellar and gas disk morphologies are mutually consistent is also unclear. Classical Cepheids, primary distance indicators with distance accuracies of 3-5%, offer a unique opportunity to develop an intuitive and accurate three-dimensional picture. Here, we establish a robust Galactic disk model based on 1339 classical Cepheids. We provide strong evidence that the warp's line of nodes is not oriented in the Galactic Center-Sun direction. Instead, it subtends a mean angle of 17.5 \pm 1 (formal) \pm 3 (systematic) and exhibits a leading spiral pattern. Our Galaxy thus follows Briggs' rule for spiral galaxies, which suggests that the origin of the warp is associated with torques forced by the massive inner disk. The stellar disk traced by Cepheids follows the gas disk in terms of their amplitudes; the stellar disk extends to at least 20 kpc. This morphology provides a crucial, updated map for studies of the kinematics and archaeology of the Galactic disk., Comment: 23 pages, 9 figures, Published in Nature Astronomy

Asteroseismic analysis of solar-like stars allows us to determine physical parameters such as stellar mass, with a higher precision compared to most other methods. Even in a well-studied cluster such as the Hyades, the masses of the red giant stars are not well known, and previous mass estimates are based on model calculations (isochrones). The four known red giants in the Hyades are assumed to be clump (core-helium-burning) stars based on their positions in colour-magnitude diagrams, however asteroseismology offers an opportunity to test this assumption. Using asteroseismic techniques combined with other methods, we aim to derive physical parameters and the evolutionary stage for the planet hosting star epsilon Tau, which is one of the four red giants located in the Hyades. We analysed time-series data from both ground and space to perform the asteroseismic analysis. By combining high signal-to-noise (S/N) radial-velocity data from the ground-based SONG network with continuous space-based data from the revised Kepler mission K2, we derive and characterize 27 individual oscillation modes for epsilon Tau, along with global oscillation parameters such as the large frequency separation and the ratio between the amplitude of the oscillations measured in radial velocity and intensity as a function of frequency. The latter has been measured previously for only two stars, the Sun and Procyon. Combining the seismic analysis with interferometric and spectroscopic measurements, we derive physical parameters for epsilon Tau, and discuss its evolutionary status., Comment: 13 pages, 13 figures, 4 tables, accepted for publication in Astronomy & Astrophysics

We present long-term, multi-color photometric monitoring and spectroscopic observations of the short-period M-type eclipsing binary BX Tri. Six flare events were recorded in 4 years from 2014 to 2017. Three of them were detected on one night within an orbital cycle. The strongest one was identified on December 23, 2014. With the amplitudes of $\Delta$B=0.48 mag, $\Delta$V=0.28 mag, $\Delta$R=0.10 mag, $\Delta$I=0.02 mag, the total energy due to the flare event was measured to be 4.08 ($\pm$0.24) $\times$$10^{34}$ erg, exceeding the superflare level ($10^{34}$). Based on the observations, the evolutionary status of the binary system as well as the long-term orbital period changes were analyzed. It reveals that BX Tri is probably a semi-detached system with the primary component filling its Roche lobe. The extremely high occurrence of flare events of the binary could be related to the rapid mass transfer between components., Comment: 16 pages, 10figures; Accepted for publication in ApJ

Blue straggler stars have been proposeed as powerful indicators to measure the dynamical state of Galactic globular clusters. Here we examine for the first time if this framework of blue straggler stars as dynamical clocks, which was specifically developed for Galactic globular clusters, may also hold for younger globular clusters in the Large Magellanic Cloud. Using Hubble Space Telescope observations, we study seven Large Magellanic Cloud star clusters with ages between ~700 Myr and 7 Gyr. We confirm that our sample clusters are consistent with the empirical correlation previously derived for Galactic globular clusters, within a modest tolerance. We suggest that to further examine if blue straggler stars can measure the dynamical state of their host clusters over an extended range of dynamical ages, more studies of dynamically older Magellanic Cloud clusters are required. We discuss the physical implications of our results in terms of their central, dimensionless King potential, as well as the initial retention fraction of black holes., Comment: 16 pages, 15 figures, accepted by the ApJ

Li et al. (2018) (hereafter L18) detected two apparently distinct populations of blue straggler stars (BSSs) in the young globular cluster NGC 2173, a similar feature as observed in numerous Galactic globular clusters (GCs). Recently, Dalessandro et al. (2018) (D18) compared the observed color--magnitude diagrams (CMDs) of both the cluster and a nearby reference field (although observed with the ACS/WFC instrument). They conclude that the bifurcated pattern of BSSs in NGC 2173 observed by L18 is a field contamination artefact. In this note, we explore the central concentration properties of the removed `field stars' identified by D18. Our purpose is to examine if these `field stars' are spatially homogeneously distributed. Employing a Monte Carlo-based approach, we have carefully studied the probability that any such central concentration may be caused by small number statistics. We find that, in most cases of, the `field stars' removed by D18 exhibit a clear central concentration, which cannot be explained on the basis of small number statistics alone. Therefore, we suggest that D18 may well have overestimated the field contamination level, implying that the bifurcated BSS pattern in NGC 2173 cannot, in fact, be explained by field contamination., Comment: 5 pages, 3 figures, extended version of the RNAAS 2:215, 2018

Although the formation of blue straggler stars (BSSs) is routinely attributed to stellar interactions in binary systems, the relative importance of the direct collision and slow(er) stellar coalescence formation channels is still poorly understood. We selected a sample of 24 Magellanic Cloud star clusters for which multi-passband Hubble Space Telescope images are available to address this outstanding question. We compiled a BSS database, containing both traditional and evolved BSSs. We found a robust correlation between the number of BSSs in a cluster's core and its core mass, $N_{\rm BSS, core} \propto M_{\rm core}^{0.51 \pm 0.07}$, which supports the notion that BSS formation is linked to a population's binary fraction. At low stellar collision rates, the mass-normalised number of BSSs does not appear to depend on the collision rate, which implies that the coalescence-driven BSS formation channel dominates. Comparison with simulations suggests that stellar collisions contribute less than 20% to the total number of BSSs formed., Comment: 7 pages, 2 figures; to appear in: Instability Phenomena and Evolution of the Universe, Byurakan Astrophys. Obs. (Armenia), September 2018, Commun. Byurakan Astrophys. Obs., 65, in press

We model the oscillations of the SONG target $\mu$ Herculis to estimate the parameters of the star. The $\ell$ = 1 mixed modes of $\mu$ Her provide strong constraints on stellar properties. The mass and age given by our asteroseismic modelling are 1.10$^{+0.11}_{-0.06}$ M$_{\odot}$ and 7.55$^{+0.96}_{-0.79}$ Gyr. The initial helium abundance is also constrained at around $Y_{\rm{init}}$ = 0.28, suggesting a ratio in the elements enrichment law ($\Delta Y/\Delta Z$) around 1.3, which is closed to the solar value. The mixing-length parameter converges to about 1.7, which is $\sim$ 10\% lower than the solar value and consistent with the results from hydrodynamic simulations. Our estimates of stellar mass and age agree very well with the previous modelling results with different input physics. Adding asteroseismic information makes these determinations less model-dependent than is typically the case when only surface information is available. Our studies of the model dependence (mass, initial helium and metallicity fractions, and the mixing length parameter) of the age determination indicate that accurate stellar ages ($\lesssim$ 10\%) can be expected from asteroseismic modelling for stars similar to $\mu$ Her. The $\ell$ = 1 bumped modes, which are sensitive to the mean density of the helium core, provide a useful 'clock' that provides additional constraints on its age., Comment: 11 pages, 9 figures

We report the discovery of a new eclipsing EL CVn-type binary, consisting of a $\delta$ Sct-type pulsator and a thermally bloated low-mass pre-He white dwarf (WD). Spectroscopy and time-series $BV$ photometry of V1224 Cas were carried out. The spectroscopy reveals a spectral type of A3 for the star. Light-curve modelling indicate that V1224 Cas is a short-period detached system containing a possible low-mass WD with an effective temperature of about 9516 K. Based on the effective temperature and the surface gravity of the A-type primary star from the spectroscopic results, the absolute parameters of the components were estimated as: $M_{P}$=2.16$\pm$0.22$M_{\sun}$, $R_{P}$=3.54$\pm$0.12$R_{\sun}$, $L_{P}$=55.9$\pm$6.9$L_{\sun}$, $M_{S}$=0.19$\pm$0.02$M_{\sun}$, $R_{S}$=0.97$\pm$0.04$R_{\sun}$, and $L_{S}$=6.9$\pm$0.6$L_{\sun}$. We therefore introduce V1224 Cas as a new EL CVn-type binary candidate. The light curves in both filters all show multi-periodic pulsations, superimposed on binary effects. We performed a preliminary frequency analysis of the light residuals after removing the synthetic eclipsing curve from the original observational data. The results suggest that the rapid light variations among the light curves could be attributed to the $\delta$ Sct-type primary component. We therefore conclude that V1224 Cas is very likely a WD+$\delta$ Sct binary., Comment: Accepted for publication in AJ. 13 pages, 6 figures

We analyze {\sl Hubble Space Telescope} observations of nine Large Magellanic Cloud star clusters with ages of 1--2 Gyr to search for evolved counterparts of blue straggler stars. Near the red clump regions in the clusters' color--magnitude diagrams, we find branches of evolved stars that are much brighter than normal evolved stars. We examine the effects of photometric artifacts, differential reddening, and field contamination. We conclude that these bright evolved stars cannot be explained by any of these effects. Our statistical tests show that the contributions of photometric uncertainties and crowding effects, as well as that owing to differential reddening, to these bright evolved stars are insufficient to fully explain the presence of these bright evolved stars. Based on isochrone fitting, we also ruled out the possibility that these bright evolved stars could be reproduced by an internal chemical abundance spread. The spatial distributions of the bright evolved stars exhibit clear concentrations that cannot be explained by homogeneously distributed field stars. This is further confirmed based on Monte Carlo-based tests. By comparing our observations with stellar evolution models, we find that the masses of most of bright evolved stars do not exceed twice the average mass of normal evolved stars. We suggest that these bright evolved stars are, in fact, evolved blue straggler stars., Comment: 19 pages, 17 figures, AJ accepted

We have compiled the first all-sky mid-infrared variable-star catalog based on Wide-field Infrared Survey Explorer (WISE) five-year survey data. Requiring more than 100 detections for a given object, 50,282 carefully and robustly selected periodic variables are discovered, of which 34,769 (69%) are new. Most are located in the Galactic plane and near the equatorial poles. A method to classify variables based on their mid-infrared light curves is established using known variable types in the General Catalog of Variable Stars. Careful classification of the new variables results in a tally of 21,427 new EW-type eclipsing binaries, 5654 EA-type eclipsing binaries, 1312 Cepheids, and 1231 RR Lyraes. By comparison with known variables available in the literature, we estimate that the misclassification rate is 5% and 10% for short- and long-period variables, respectively. A detailed comparison of the types, periods, and amplitudes with variables in the Catalina catalog shows that the independently obtained classifications parameters are in excellent agreement. This enlarged sample of variable stars will not only be helpful to study Galactic structure and extinction properties, they can also be used to constrain stellar evolution theory and as potential candidates for the James Webb Space Telescope., Comment: 30 pages, 13 figures, accepted by The Astrophysical Journal Supplement Series

We have analyzed populations of blue straggler stars (BSSs) in 24 Magellanic Cloud star clusters using multi-passband Hubble Space Telescope images. We compiled a homogeneous BSS database, containing both traditional and evolved BSSs. We uncovered a sub-linear correlation between the number of BSSs in the cluster cores and the clusters' core masses, characterized by a power-law index of $0.51\pm 0.07$. For low stellar collision rates, the mass-normalized number of BSSs depends only weakly (or perhaps not at all) on the collision rate, implying that the binary-driven BSS formation channel dominates. Comparison with simulations suggests that stellar collisions contribute less than 20\% to the total number of BSSs formed. Further tests, including analysis of the BSS specific frequencies and their population numbers at larger cluster radii, suggest that binary interactions may be their main formation channel, hinting at an anti-correlation between a cluster's binary fraction and its core mass., Comment: 14 pages, 5 figures. Accepted for publication in ApJ

We present {complexity} of the Galactic outer disc {by fitting the stellar volume densities of the red giant branch stars with a two-disc component model}. {The discs are confirmed to} extend to $R\sim19$\,kpc. The radial density profile of the discs shows {two breaks at $R\sim11$ and $\sim14$\,kpc, respectively, which} separate the radial profile into three segments with different scale lengths of $2.12\pm0.26$, $1.18\pm0.08$, and $2.72$\,kpc at $R<11$, $11\leq R\leq14$, and $R>14$\,kpc, respectively. The first {break} is likely due to the sudden drop in the radial profile of the thin disc, which may be an evidence of the radial migration. {Beyond $14$\,kpc, the thick disc becomes prominent and the transition from thin to thick disc leads to the second break.} This implies that the geometrically defined thick disc is more {radially} extended than the thin disc. This is also supported by the larger scale length of the thick disc than that of the thin disc. Meanwhile, {the scale height of the thicker component increases from $0.637_{-0.036}^{+0.056}$ at $R=8$ to $1.284_{-0.079}^{+0.086}$\,kpc at $R=19$\,kpc, showing an intensive flared disc}. Moreover, rich substructures are displayed in the residuals of the stellar density. Among them, the substructures $D14+2.0$ and $O14-1.5$ show a north-south asymmetry, which can be essentially explained by southward shifting of the thick disc. However, no significant overdensity is found for the Monoceros ring. Finally, the thick disc shows a ripple-like feature with unclear origin at $9

Distances and extinction values are usually degenerate. To refine the distance to the general Galactic Center region, a carefully determined extinction law (taking into account the prevailing systematic errors) is urgently needed. We collected data for 55 classical Cepheids projected toward the Galactic Center region to derive the near- to mid-infrared extinction law using three different approaches. The relative extinction values obtained are AJ/AKs = 3.005, AH/AKs = 1.717, A[3.6]/AKs = 0.478, A[4.5]/AKs = 0.341, A[5.8]/AKs = 0.234, A[8.0]/AKs = 0.321, AW1/AKs = 0.506, and AW2/AKs = 0.340. We also calculated the corresponding systematic errors. Compared with previous work, we report an extremely low and steep mid-infrared extinction law. Using a seven-passband 'optimal distance' method, we improve the mean distance precision to our sample of 55 Cepheids to 4%. Based on four confirmed Galactic Center Cepheids, a solar Galactocentric distance of R_0 = 8.10\pm0.19\pm0.22 kpc is determined, featuring an uncertainty that is close to the limiting distance accuracy (2.8%) for Galactic Center Cepheids., Comment: 29 pages, 6 figures, accepted for publication by The Astrophysical Journal

W Ursa Majoris (W UMa)-type contact binary systems (CBs) are useful statistical distance indicators because of their large numbers. Here, we establish (orbital) period-luminosity relations (PLRs) in 12 optical-to-mid-infrared bands (GBV RIJHKsW1W2W3W4) based on 183 nearby W UMa-type CBs with accurate Tycho-Gaia parallaxes. The 1{\sigma} dispersion of the PLRs decreases from optical to near- and mid-infrared wavelengths. The minimum scatter, 0.16 mag, implies that W UMa-type CBs can be used to recover distances to 7% precision. Applying our newly determined PLRs to 19 open clusters containing W UMa-type CBs demonstrates that the PLR and open cluster CB distance scales are mutually consistent to within 1%. Adopting our PLRs as secondary distance indicators, we compiled a catalog of 55,603 CBs candidates, of which 80% have distance estimates based on a combination of optical, near-, and mid-infrared photometry. Using Fourier decomposition, 27,318 high-probability W UMa-type CBs were selected. The resulting 8% distance accuracy implies that our sam- ple encompasses the largest number of objects with accurate distances within a local volume with a radius of 3 kpc available to date. The distribution of W UMa-type CBs in the Galaxy suggests that in different environments, the CB luminosity function may be different: larger numbers of brighter (longer-period) W UMa-type CBs are found in younger environments., Comment: 26 pages, 6 figures, accepted by The Astrophysical Journal

Recent discoveries of bimodal main sequences (MSs) associated with young clusters (with ages $\lesssim 1$ Gyr) in the Magellanic Clouds have drawn a lot of attention. One of the prevailing formation scenarios attributes these split MSs to a bimodal distribution in stellar rotation rates, with most stars belonging to a rapidly rotating population. In this scenario, only a small fraction of stars populating a secondary blue sequence are slowly or non-rotating stars. Here, we focus on the blue MS in the young cluster NGC 1850. We compare the cumulative number fraction of the observed blue-MS stars to that of the high-mass-ratio binary systems at different radii. The cumulative distributions of both populations exhibit a clear anti-correlation, characterized by a highly significant Pearson coefficient of $-0.97$. Our observations are consistent with the possibility that blue-MS stars are low-mass-ratio binaries, and therefore their dynamical disruption is still ongoing. High-mass-ratio binaries, on the other hand, are more centrally concentrated., Comment: 8 pages, 6 figures, accepted to ApJ

We present a three dimensional velocity analysis of Milky Way disk kinematics using LAMOST K giant stars and the GPS1 proper motion catalogue. We find that Galactic disk stars near the anticenter direction (in the range of Galactocentric distance between $R=8$ and $13$ kpc and vertical position between $Z=-2$ and $2$ kpc) exhibit asymmetrical motions in the Galactocentric radial, azimuthal, and vertical components. Radial motions are not zero, thus departing from circularity in the orbits; they increase outwards within $R\lesssim 12$ kpc, show some oscillation in the northern ($0 < Z < 2$ kpc) stars, and have north-south asymmetry in the region corresponding to a well-known nearby northern structure in the velocity field. There is a clear vertical gradient in azimuthal velocity, and also an asymmetry that shifts from a larger azimuthal velocity above the plane near the solar radius to faster rotation below the plane at radii of 11-12 kpc. Stars both above and below the plane at $R\gtrsim 9$ kpc exhibit net upward vertical motions. We discuss some possible mechanisms that might create the asymmetrical motions, such as external perturbations due to dwarf galaxy minor mergers or dark matter sub-halos, warp dynamics, internal processes due to spiral arms or the Galactic bar, and (most likely) a combination of some or all of these components., Comment: 10 pages, 14 figures. MNRAS in press

Stellar fundamental parameters are important in the asteroseismic study of Kepler light curves. However, the most used estimates in the Kepler Input Catalog (KIC) are not accurate enough for hot stars. Using a sample of B stars from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) spectral survey, we confirmed the systematic underestimation in KIC effective temperature and overestimation in KIC surface gravity. The good agreement between LAMOST and other follow-up observations proved the accuracy of effective temperature and surface gravity of B stars derived from LAMOST low-resolution spectra. By searching through LAMOST data, we found four misclassified main-sequence B stars in the Kepler field, which had been previously classified as A-type variables. We present spectroscopic and detailed frequency analysis of these four stars based on LAMOST spectra and Kepler photometry., Comment: 19 pages, 12 figures, 7 tables, to appear in ApJ

Bifurcated patterns of blue straggler stars in their color--magnitude diagrams have atracted significant attention. This type of special (but rare) pattern of two distinct blue straggler sequences is commonly interpreted as evidence of cluster core-collapse-driven stellar collisions as an efficient formation mechanism. Here, we report the detection of a bifurcated blue straggler distribution in a young Large MagellanicCloud cluster, NGC 2173. Because of the cluster's low central stellar number density and its young age, dynamical analysis shows that stellar collisions alone cannot explain the observed blue straggler stars. Therefore, binary evolution is instead the most viable explanation of the origin of these blue straggler stars. However, the reason why binary evolution would render the color--magnitude distribution of blue straggler stars bifurcated remains unclear., Comment: 9 pages, 8 figures, ApJ accepted

Classical Cepheids are well-known and widely used distance indicators. As distance and extinction are usually degenerate, it is important to develop suitable methods to robustly anchor the distance scale. Here, we introduce a near-infrared (near-IR) optimal distance method to determine both the extinction values of and distances to a large sample of 288 Galactic classical Cepheids. The overall uncertainty in the derived distances is less than 4.9%. We compare our newly determined distances to the Cepheids in our sample with previously published distances to the same Cepheids with Hubble Space Telescope parallax measurements and distances based on the IR surface brightness method, Wesenheit functions, and the main-sequence fitting method. The systematic deviations in the distances determined here with respect to those of previous publications is less than 1-2%. Hence, we constructed Galactic mid-IR period-luminosity (PL) relations for classical Cepheids in the four Wide-Field Infrared Survey Explorer (WISE) bands (W1, W2, W3, and W4) and the four Spitzer Space Telescope bands ([3.6], [4.5], [5.8] and [8.0]). Based on our sample of hundreds of Cepheids, the WISE PL relations have been determined for the first time; their dispersion is approximately 0.10 mag. Using the currently most complete sample, our Spitzer PL relations represent a significant improvement in accuracy, especially in the [3.6] band which has the smallest dispersion (0.066 mag). In addition, the average mid-IR extinction curve for Cepheids has been obtained: A_W1/A_Ks=0.560, A_W2/A_Ks=0.479, A_W3/A_Ks=0.507, A_W4/A_Ks=0.406, A_[3.6]/A_Ks=0.481, A_[4.5]/A_Ks=0.469, A_[5.8]/A_Ks=0.427, and A_[8.0]/A_Ks=0.427 mag., Comment: 19 pages, 5 figures, accepted for publication in The Astrophysical Journal

Based on an extensive spectral study of a photometrically confirmed sample of Mira variables, we find a relationship between relative Balmer emission-line strength and spectral temperature of O-rich Mira stars. The $F_{\rm H\delta}/F_{\rm H\gamma}$ flux ratio increases from less than unity to five as stars cool down from M0 to M10, which is likely driven by increasing TiO absorption above the deepest shock-emitting regions. We also discuss the relationship between the equivalent widths of the Balmer emission lines and the photometric luminosity phase of our Mira sample stars. Using our 291 Mira spectra as templates for reference, 191 Mira candidates are newly identified from the LAMOST DR4 catalog. We summarize the criteria adopted to select Mira candidates based on emission-line indices and molecular absorption bands. This enlarged spectral sample of Mira variables has the potential to contribute significantly to our knowledge of the optical properties of Mira stars and will facilitate further studies of these late-type, long-period variables., Comment: 21 pages; ApJS, in press

Numerical simulations were carried out to study the origin of multiple stellar populations in the intermediate-age clusters NGC 411 and NGC 1806 in the Magellanic Clouds. We performed NBODY6++ simulations based on two different formation scenarios, an ad hoc formation model where second-generation (SG) stars are formed inside a cluster of first-generation (FG) stars using the gas accumulated from the external intergalactic medium and a minor merger model of unequal mass ($M_{SG}$/$M_{FG}$ ~5-10%) clusters with an age difference of a few hundred million years. We compared our results such as the radial profile of the SG-to-FG number ratio with observations on the assumption that the SG stars in the observations are composed of cluster members, and confirmed that both the ad hoc formation and merger scenarios reproduce the observed radial trend of the SG-to-FG number ratio which shows less centrally concentrated SG than FG stars. It is difficult to constrain the formation scenario for the multiple populations by only using the spatial distribution of the SG stars. SG stars originating from the merger scenario show a significant velocity anisotropy and rotational features compared to those from the ad hoc formation scenario. Thus, observations aimed at kinematic properties like velocity anisotropy or rotational velocities for SG stars should be obtained to better understand the formation of the multiple populations in these clusters. This is, however, beyond current instrumentation capabilities., Comment: 11 pages, 8 figures, 1 table; accepted for publication in MNRAS