Search

We present an analysis of the structure, kinematics, and chemo-dynamical properties of the Milky Way bulge using data from OGLE, APOGEE, and Gaia. Firstly, we identified 2,156 ab-type RR Lyrae stars (RRabs) from OGLE-IV, then through their apocenters derived from orbital integration, we distinguished three populations: the central bulge RRabs, the inner bulge RRabs and halo interlopers. Inner bulge RRabs kinematically align with the Galactic bar, while central bulge RRabs show slower rotation with lower velocity dispersion, which do not trace the bar. Higher velocity dispersion stars were identified as halo interlopers. Then, orbital analysis of 28,188 APOGEE Red Clump and Red Giant Branch stars revealed kinematic properties consistent with RRabs, and chemical abundance distribution displayed a bimodal stellar density pattern, suggesting complex star evolution histories and also slightly different star formation histories for the inner bulge and central bulge. The differences in the density distribution on the $|\mathrm{Z}|_{\text{max}}$-eccentricity plane for the central bulge, inner bulge, and halo regions are clearly detected. Finally, the chemodynamical analysis of 301,485 Gaia DR3 red giants without orbital integration indicated that metal-rich bulge stars form a bar-like structure, while metal-poor bulge stars are dominated by velocity dispersion. It is found that the classification of bulge stars based on orbital parameters, rather than solely on metallicity, provides a more accurate population separation. Our results also support that secular evolution of the Galactic disk is the primary origin of the bulge, and boxy/peanut (B-P) bulge population might be more dominant than X-shape bulge population., Comment: 24 pages, 20 figures

The Wide-field Spectroscopic Telescope (WST) is proposed as a new facility dedicated to the efficient delivery of spectroscopic surveys. This white paper summarises the initial concept as well as the corresponding science cases. WST will feature simultaneous operation of a large field-of-view (3 sq. degree), a high multiplex (20,000) multi-object spectrograph (MOS) and a giant 3x3 sq. arcmin integral field spectrograph (IFS). In scientific capability these requirements place WST far ahead of existing and planned facilities. Given the current investment in deep imaging surveys and noting the diagnostic power of spectroscopy, WST will fill a crucial gap in astronomical capability and work synergistically with future ground and space-based facilities. This white paper shows that WST can address outstanding scientific questions in the areas of cosmology; galaxy assembly, evolution, and enrichment, including our own Milky Way; origin of stars and planets; time domain and multi-messenger astrophysics. WST's uniquely rich dataset will deliver unforeseen discoveries in many of these areas. The WST Science Team (already including more than 500 scientists worldwide) is open to the all astronomical community. To register in the WST Science Team please visit https://www.wstelescope.com/for-scientists/participate, Comment: 194 pages, 66 figures. Comments are welcome (wstelescope@gmail.com)

We present a novel, deep-learning based method -- dubbed Galactic-Seismology Substructures and Streams Hunter, or GS$^{3}$ Hunter for short, to search for substructures and streams in stellar kinematics data. GS$^{3}$ Hunter relies on a combined application of Siamese Neural Networks to transform the phase space information and the K-means algorithm for the clustering. As a validation test, we apply GS$^{3}$ Hunter to a subset of the Feedback in Realistic Environments (FIRE) cosmological simulations. The stellar streams and substructures thus identified are in good agreement with corresponding results reported earlier by the FIRE team. In the same vein, we apply our method to a subset of local halo stars from the Gaia Early Data Release 3 and GALAH DR3 datasets, and recover several, previously known dynamical groups, such as Thamnos 1+2, Hot Thick Disk, ED-1, L-RL3, Helmi 1+2, and Gaia-Sausage-Enceladus, Sequoia, VRM, Cronus, Nereus. Finally, we apply our method without fine-tuning to a subset of K-giant stars located in the inner halo region, obtained from the LAMOST Data Release 5 (DR5) dataset. We recover three, previously known structures (Sagittarius, Hercules-Aquila Cloud, and the Virgo Overdensity), but we also discover a number of new substructures. We anticipate that GS$^{3}$ Hunter will become a useful tool for the community dedicated to the search of stellar streams and structures in the Milky Way (MW) and the Local group, thus helping advance our understanding of the stellar inner and outer halos, and of the assembly and tidal stripping history in and around the MW., Comment: 41 pages, 34 figures, 3 tables. Accepted for Publication in the Astrophysical Journal

The unexpected discovery of hot Jupiters challenged the classical theory of planet formation inspired by our solar system. Until now, the origin and evolution of hot Jupiters are still uncertain. Determining their age distribution and temporal evolution can provide more clues into the mechanism of their formation and subsequent evolution. Using a sample of 383 giant planets around Sun-like stars collected from the kinematic catalogs of the Planets Across Space and Time (PAST) project, we find that hot Jupiters are preferentially hosted by relatively younger stars in the Galactic thin disk. We subsequently find that the frequency of hot Jupiters declines with age. In contrast, the frequency of warm/cold Jupiters shows no significant dependence on age. Such a trend is expected from the tidal evolution of hot Jupiters' orbits, and our result offers supporting evidence using a large sample. We also perform a joint analysis on the planet frequencies in the stellar age-metallicity plane. The result suggests that the frequencies of hot Jupiters and warm/cold Jupiters, after removing the age dependence are both correlated with stellar metallicities. Moreover, we show that the above correlations can explain the bulk of the discrepancy in hot Jupiter frequencies inferred from the transit and radial velocity (RV) surveys, given that RV targets tend to be more metal-rich and younger than transits., Comment: Published in PNAS; 7 pages, 5 figures in the main text; 17 pages, 29 figures, 5 tables in the supporting information

In this study we estimate the stellar ages of LAMOST DR8 Red Giant Branch (RGB) stars based on the Gradient Boosting Decision Tree algorithm (GBDT). We used 2,643 RGB stars extracted from the APOKASC-2 astero-seismological catalog as training data-set. After selecting the parameterses ([$\alpha$/Fe], [C/Fe], T$_{eff}$, [N/Fe], [C/H], log g) highly correlated with age using GBDT, we apply the same GBDT method to the new catalog of more than 590,000 stars classified as RGB stars. The test data-set shows that the median relative error is around 11.6$\%$ for the method. We also compare the predicted ages of RGB stars with other studies (e.g., based on APOGEE), and find systematic differences. The final uncertainty is about 15 to 30$\%$ compared to open clusters' ages. Then we present the spatial distribution of the RGB sample having an age determination, which could recreate the expected result, and discuss systematic biases. All these diagnostics show that one can apply the GBDT method to other stellar samples to estimate atmospheric parameters and age., Comment: 17 pages, 17 figures, Accepted for publication in The Astrophysical Journal

The application of asymmetric drift (AD) tomography across different populations provides valuable insights into the kinematics, dynamics, and rotation curves of the Galactic disk. By leveraging common stars identified in both the LAMOST and Gaia surveys, alongside Gaia DR3's circular velocity curve, we conducted a qualitative exploration of asymmetric drift distributions within the Galactic disk spanning distances from 8 to 16 kpc. In the R-Z plane, we observed that the asymmetric drift is minimal near the mid-plane of the Galactic disk and gradually increases with vertical distance, resulting in a distinctive ``horn" shape. Additionally, our analysis revealed that populations with higher [$\alpha$/Fe] ratios exhibit greater asymmetric drift compared to those with lower [$\alpha$/Fe] ratios. Specifically, we found the asymmetric drift around the solar location to be approximately 6 km s$^{-1}$, with a median value of 16 km s$^{-1}$ across the entire sample. Notably, the median asymmetric drift in the northern region of the Galactic disk (20 km s$^{-1}$) surpasses that in the southern region (13 km s$^{-1}$), with errors remaining within 2 km s$^{-1}$. Furthermore, our investigation into mono-age stellar populations unveiled that older stellar populations tend to exhibit larger asymmetric drift and velocity dispersion, aligning closely with predictions from previous numerical models. Finally, based on chemical compositions, we observed that the median asymmetric drift of the thick disk significantly exceeds that of the thin disk and found that star formation within the thick disk primarily occurred earlier than 8-10 billion years, whereas the thin disk's predominant star formation period spanned 6-8 billion years ago., Comment: 14 pages and 17 figures. The Open Journal of Astrophysics. AD analysis is the key target and the time scale picture of the milky way disk are also revealed here in another point of view, see the Fig. 9

We present the Galactic disk vertical velocity analysis using OB type stars (OB), Red Clump stars (RC), and Main-Sequence-Turn-Off stars (MSTO) with different average age populations crossed matched with LAMOST DR5 and Gaia DR3. We reveal the vertical velocities of the three populations varies clearly with the Galactocentric distance ($R$) and the younger stellar population has stronger increasing trend in general. The bending and breathing modes indicated by the vertical motions are dependent on the populations and they are varying with spatial locations. These vertical motions may be due to the Galactic warp, or minor mergers, or non-equilibrium of the disk. Assuming the warp is the dominant component, we find that the warp amplitude ($\gamma$, $Z_\omega$) for OB (younger population) is larger than that for RC (medium population) and the later one is also larger than that for MSTO (older population), which is in agreement with other independent analyses of stellar density distribution, and supports the warp is long-lived, non-steady structure and has time evolution. This conclusion is robust whether or not the line-of-nodes $\phi_w$ is fixed or as a free parameter (with $\phi_w$ is around 3$-$8.5$^{\circ}$ as best fit). Furthermore, we find that warp is lopsided with asymmetries along azimuthal angle ($\phi$)., Comment: 15 pages, 13 figures. Accepted for publication in The Astrophysical Journal

We apply a statistical deconvolution of the parallax errors based on Lucy's inversion method (LIM) to the Gaia-DR3 sources to measure their three dimensional velocity components in the range of Galactocentric distances $R$ between 8 kpc and 30 kpc with their corresponding errors and root mean square values. We find results that are consistent with those obtained by applying LIM to the Gaia-DR2 sources, and we conclude that the method gives convergent and more accurate results by improving the statistics of the data-set and lowering observational errors. The kinematic maps reconstructed with LIM up to $R \approx 30$ kpc show that the Milky Way is characterized by asymmetrical motions with significant gradients in all velocity components. Furthermore, we determine the Galaxy rotation curve $V_C(R)$ up to $\approx 27.5$ kpc with the cylindrical Jeans equation assuming an axisymmetric gravitational potential. We find that $V_C(R)$ is significantly declining up to the largest radius investigated. Finally, we also measure $V_C(R)$ at different vertical heights, showing that, for $R <15$ kpc, there is a marked dependence on $Z$, whereas at larger $R$ the dependence on $Z$ is negligible., Comment: 14 pages, 10 figures, 1 table. Accepted for publication in The Astrophysical Journal

We investigate the kinematic-chemical distribution of Red Giant Branch (RGB) stars from the LAMOST survey crossed matched with Gaia DR2 proper motions, and present time tagging for the well-known ridge structures (diagonal distributions for $V_R$ in the $R$, $V_\phi$ plane) in the range of Galactocentric distance $R$ = 8 to 15 kpc. We detect six ridge structures, including five ridges apparent in the radial velocity distribution and three ridges apparent in the vertical velocity, the sensitive time of which to the perturbations are from young population (0$-$3 Gyr) to old population (9$-$14 Gyr). Based on an analysis of the evolution of angular momentum distribution, we find that four ridges are relatively stationary, while another is evolving with time, which is confirmed by the difference analysis at different populations and supporting that there might be two kinds of dynamical origins. Furthermore, ridge features are also vividly present in the chemical properties ([Fe/H], [$\alpha$/Fe]). The comparison between the north and south hemispheres of the Galaxy does show some differences and the ridge features are asymmetrical. Moreover, we find that diagonal ridge structures may affect the shape of the rotation curve, which is manifested as fluctuations and undulations on top of a smooth profile. Finally we speculate that the bar dynamics should be not enough to explain all ridge properties including the break feature in the $V_Z$-$L_Z$ plane., Comment: 17 pages, 14 figures. Accepted publication in The Astronomical Journal

We present a catalog of 948,216 stars with mass label and a catalog of 163,105 red clump (RC) stars with mass and age labels simultaneously. The training dataset is cross matched from the LAMOST (The Large Sky Area Multi-Object Fiber Spectroscopic Telescope) DR5 and high resolution asteroseismology data, mass and age are predicted by random forest method or convex hull algorithm. The stellar parameters with high correlation with mass and age are extracted and the test dataset shows that the median relative error of the prediction model for the mass of large sample is 3\% and meanwhile, the mass and age of red clump stars are 4\% and 7\%. We also compare the predicted age of red clump stars with the recent works and find that the final uncertainty of the RC sample could reach 18\% for age and 9\% for mass, in the meantime, final precision of the mass for large sample with different type of stars could reach 13\% without considering systematics, all these are implying that this method could be widely used in the future. Moreover, we explore the performance of different machine learning methods for our sample, including bayesian linear regression (BYS), gradient boosting decision Tree (GBDT), multilayer perceptron (MLP), multiple linear regression (MLR), random forest (RF) and support vector regression (SVR). Finally we find that the performance of nonlinear model is generally better than that of linear model, and the GBDT and RF methods are relatively better., Comment: 19 pages, 23 figures, Accepted for publication in the ApJS

The infall of the Sagittarius (Sgr) Dwarf Spheroidal Galaxy in the Milky Way halo is an unique opportunity to understand how the different components of a dwarf galaxy could be tidally removed. In this work, we reconstruct the Sgr core morphology and kinematics on the basis of a model that has already successfully reproduced the Sgr stream. Here we use a very high resolution model that almost resolves individual stars in the Sgr core. It reproduces most of the observed morphology and kinematic properties, without specific fine-tuning. We also show that the dark matter may have been almost entirely stripped by Milky Way tides after two passages at pericenter. Finally the model predicts that the Sgr core will be fully disrupted within the next 2 Gyr., Comment: 9 pages, 5 figures,1 table. Accepted for publication in ApJL

It is a challenge to reproduce the full 6D space-phase properties of Sagittarius (Sgr) dwarf galaxy and its Stream simultaneously. Using N-body simulations with a Milky Way mass of 5.2$\times10^{11}$ M$_{\odot}$ and a ``scaling down'' Sgr mass of 9.3$\times10^{8}$ M$_{\odot}$, from a qualitative point of view, we have been able to reproduce well all 3D spatial features of Sgr stream, including its core, leading and trailing arms, and their associated bifurcations, moreover, the overall trend of the reported 3D kinematics properties of the Sgr stream have also been reproduced without fine tuning. Furthermore, we also find that our model fails in reproducing the exact behaviours of the line-of-sight velocity and angular-energy distributions. It let us to suggest that significant further progress might be achievable after introducing a major component in the Sgr progenitor, which is the gas that dominates all Irregular dwarf galaxies in the Sgr mass range and can slow down the radial velocity of Sgr before its removal, if gas can not solve this problem then we will consider a non-spherical Milky Way halo with hot gas, LMC, etc. As the first step for us to towards the complete understanding of the Sgr system, this progress is also advancing our understanding of the bifurcations, the generation of which might be due to the MW shocks at each pericenter passage, and also be linked to the orientation and disk-shape in the initial conditions., Comment: 9 pages. We have ultimately decided to defer the modelling of the full Sagittarius dwarf galaxy with gas to a later contribution, which will also be part of the China-France collaboration on the Local Group. So please note that this paper II (purely N-body) has to be delayed to contribute later (Paper I: ApJL, 940, L3). The project (Full modelling of the MW+Sgr) is still ongoing and stay tuned

The Kepler telescope has discovered over 4,000 planets (candidates) by searching ? 200,000 stars over a wide range of distance (order of kpc) in our Galaxy. Characterizing the kinematic properties (e.g., Galactic component membership and kinematic age) of these Kepler targets (including the planet (candidate) hosts) is the first step towards studying Kepler planets in the Galactic context, which will reveal fresh insights into planet formation and evolution. In this paper, the second part of the Planets Across the Space and Time (PAST) series, by combining the data from LAMOST and Gaia and then applying the revised kinematic methods from PAST I, we present a catalog of kinematic properties(i.e., Galactic positions, velocities, and the relative membership probabilities among the thin disk, thick disk, Hercules stream, and the halo) as well as other basic stellar parameters for 35,835 Kepler stars. Further analyses of the LAMOST-Gaia-Kepler catalog demonstrate that our derived kinematic age reveals the expected stellar activity-age trend. Furthermore, we find that the fraction of thin(thick) disk stars increases (decreases) with the transiting planet multiplicity (Np = 0, 1, 2 and 3+) and the kinematic age decreases with Np, which could be a consequence of the dynamical evolution of planetary architecture with age. The LAMOST-Gaia-Kepler catalog will be useful for future studies on the correlations between the exoplanet distributions and the stellar Galactic environments as well as ages., Comment: Accepted for publication in Astronomical Journal, 19 pages, 9 figures, 5 tables. The full data of the LAMOST-Gaia-Kepler kinematic catalog (Table 3) will be published with the journal. If you need to use the data, please contact us via the email (jwxie@nju.edu.cn, dg1826003@smail.nju.edu.cn)

Over 4,000 exoplanets have been identified and thousands of candidates are to be confirmed. The relations between the characteristics of these planetary systems and the kinematics, Galactic components, and ages of their host stars have yet to be well explored. Aiming to addressing these questions, we conduct a research project, dubbed as PAST (Planets Across Space and Time). To do this, one of the key steps is to accurately characterize the planet host stars. In this paper, the Paper I of the PAST series, we revisit the kinematic method for classification of Galactic components and extend the applicable range of velocity ellipsoid from about 100 pc to 1, 500 pc from the sun in order to cover most known planet hosts. Furthermore, we revisit the Age-Velocity dispersion Relation (AVR), which allows us to derive kinematic age with a typical uncertainty of 10-20% for an ensemble of stars. Applying the above revised methods, we present a catalog of kinematic properties (i.e. Galactic positions, velocities, the relative membership probabilities among the thin disk, thick disk, Hercules stream, and the halo) as well as other basic stellar parameters for 2,174 host stars of 2,872 planets by combining data from Gaia, LAMOST, APOGEE, RAVE, and the NASA exoplanet archive. The revised kinematic method and AVR as well as the stellar catalog of kinematic properties and ages lay foundation for future studies on exoplanets from two dimensions of space and time in the Galactic context., Comment: 26 pages, 20 figure, 10 tables, Accepted for publication in ApJ. The full data of the planet kinematic catalog (Table 10) will be published with the journal. If you need to use the data, please contact us via the email (jwxie@nju.edu.cn, dg1826003@smail.nju.edu.cn)

We present analysis of the spatial density structure for the outer disk from 8$-$14 \,kpc with the LAMOST DR5 13534 OB-type stars and observe similar flaring on north and south sides of the disk implying that the flaring structure is symmetrical about the Galactic plane, for which the scale height at different Galactocentric distance is from 0.14 to 0.5 \,kpc. By using the average slope to characterize the flaring strength we find that the thickness of the OB stellar disk is similar but flaring is slightly stronger compared to the thin disk as traced by red giant branch stars, possibly implying that secular evolution is not the main contributor to the flaring but perturbation scenarios such as interactions with passing dwarf galaxies should be more possible. When comparing the scale height of OB stellar disk of the north and south sides with the gas disk, the former one is slightly thicker than the later one by $\approx$ 33 and 9 \,pc, meaning that one could tentatively use young OB-type stars to trace the gas properties. Meanwhile, we unravel that the radial scale length of the young OB stellar disk is 1.17 $\pm$ 0.05 \,kpc, which is shorter than that of the gas disk, confirming that the gas disk is more extended than stellar disk. What is more, by considering the mid-plane displacements ($Z_{0}$) in our density model we find that almost all of $Z_{0}$ are within 100 \,pc with the increasing trend as Galactocentric distance increases., Comment: 14 pages and 8 figures, Accepted by APJ for publication

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