Search

We present the discoveries of two giant planets orbiting the red giant branch (RGB) star HD 112570 and the red clump (RC) star HD 154391, based on the radial velocity (RV) measurements from Xinglong station and Okayama Astrophysical Observatory (OAO). Spectroscopic and asteroseismic analyses suggest that HD 112570 has a mass of $1.15\pm0.12\,M_{\odot}$, a radius of $9.85\pm0.23\,R_{\odot}$, a metallicity [Fe/H] of $-0.46\pm0.1$ and a ${\rm log}\,g$ of $2.47\pm0.1$. With the joint analysis of RV and Hipparcos-Gaia astrometry, we obtain a dynamical mass of $M_{\rm p}={3.42}_{-0.84}^{+1.4}\ M_{\rm Jup}$, a period of $P={2615}_{-77}^{+85}$ days and a moderate eccentricity of $e={0.20}_{-0.14}^{+0.16}$ for the Jovian planet HD 112570 b. For HD 154391, it has a mass of $2.07\pm0.03\,M_{\odot}$, a radius of $8.56\pm0.05\,R_{\odot}$, a metallicity [Fe/H] of $0.07\pm0.1$ and a ${\rm log}\,g$ of $2.86\pm0.1$. The super-Jupiter HD 154391 b has a mass of $M_{\rm p}={9.1}_{-1.9}^{+2.8}\ M_{\rm Jup}$, a period of $P={5163}_{-57}^{+60}$ days and an eccentricity of $e={0.20}_{-0.04}^{+0.04}$. We found HD 154391 b has one of the longest orbital period among those ever discovered orbiting evolved stars, which may provide a valuable case in our understanding of planetary formation at wider orbits. Moreover, while a mass gap at $4\,M_{\rm Jup}$ seems to be present in the population of giant stars, there appears to be no significant differences in the distribution of metallicity among giant planets with masses above or below this threshold. Finally, The origin of the abnormal accumulation near 2 au for planets around large evolved stars ($R_{\star}>21\,R_{\odot}$), remains unclear., Comment: 26 pages, 18 figures, 8 tables, Accepted for publication in The Astronomical Journal

In N-body simulations, nearly radial mergers can form shell-like overdensities in the sky position and phase space ($r-v_r$) due to the combination of dynamical friction and tidal stripping. The merger event of Gaia-Sausage-Enceladus has provided a unique opportunity to study the shells in the phase space. To search for them, we integrate the orbits of 5949 GSE-related halo K giants from the LAMOST survey and record their positions at all time intervals in $r-v_r$ diagram. After the subtraction of a smoothed background, we find six significant and complete thin chevron-like overdensities. The apocenters $r_\mathrm{apo}$ of stars in the six chevrons are around 6.75, 12.75, 18.75, 25.25, 27.25, and 30.25 kpc. These chevrons reveal the multiple pile-ups of GSE stars at different apocenters. The application of a different Milky Way mass $M_\mathrm{vir}$ will change the opening angles of these chevrons, while leave their apocenters almost unchanged. By comparing with a recent study of the phase space overdensities of local halo stars from Gaia RVS survey, our results are more inclined to a medium $M_\mathrm{vir}$ of $10^{12}\,M_\odot$. The application of a non-axisymmetric Galactic potential with a steadily rotating bar has a blurring effect on the appearance of these chevron-like overdensities, especially for the chevrons with $r_\mathrm{apo} > 20$ kpc., Comment: Accepted by ApJ, 8 figures

We present a catalog of distances for 19544 K giants drawn from LAMOST DR8. Most of them are located in the halo of the Milky Way up to ~120~kpc. There are 15% K giants without SDSS photometry, for which we supplements with Pan-STARRS1 (PS1) photometry calibrated to SDSS photometric system. The possible contamination of the red clumps/horizontal branch are removed according to metallicities and colors before the distance determination. Combining the LAMOST spectroscopic metallicities with the SDSS/PS1 photometry, we estimate the absolute magnitudes in SDSS $r-$band, the distance moduli, and the corresponding uncertainties through an Bayesian approach devised by Xue et al. (2014) for the SEGUE halo K-giants. The typical distance precision is about 11%. The stars in the catalog lie in a region of 4-126 kpc from the Galactic center, of which with 6, 320 stars beyond 20 kpc and 273 stars beyond 50 kpc, forming the largest spectroscopic sample of distant tracers in the Milky Way halo so far., Comment: Accepted for publication in AJ

The halo stars on highly radial orbits should inevitably pass the center regions of the Milky Way. Under the assumption that the stellar halo is in dynamical equilibrium and axisymmetric, we integrate the orbits of $\sim 10,000$ halo K-giants at $5\leq r \leq 50$ kpc cross-matched from LAMOST DR5 and $Gaia$ DR3. By carefully considering the selection function, we construct the stellar halo distribution at the entire regions of $r \leq 50$ kpc. We find that a double-broken power-law function well describes the stellar halo density distribution with shallower slopes in the inner regions and the two breaks at $r=10$ kpc and $r=25$ kpc, respectively. The stellar halo becomes flatter from outer to inner regions but has $q\sim 0.5$ at $r \lesssim 5$ kpc. The stellar halo becomes isotropic with a slight prograde rotation in the inner 5 kpc, and reaches velocity dispersions of $\sim 250\rm \ km\ s^{-1}$. We get a weak negative metallicity gradient of $-0.005$ dex kpc$^{-1}$ at $5\leq r \leq 50$ kpc, while there is an excess of relative metal-rich stars with [Fe/H]$>-1$ in the inner 10 kpc. The halo interlopers at $r \leq 5$ kpc from integration of our sample has a mass of $\sim1.2 \times 10^8\ M_{\odot}$ ($\sim 4.7 \times 10^7\ M_{\odot}$ at [Fe/H]$<-1.5$), which can explain 50-100% of the metal-poor stars with [Fe/H]$<-1.5$ directly observed in the Galactic central regions., Comment: 22 pages, 12 figures, Accepted for publication in AJ

We measure the enclosed Milky Way mass profile to Galactocentric distances of $\sim70$ and $\sim50$ kpc using the smooth, diffuse stellar halo samples of Bird et al. The samples are LAMOST and SDSS/SEGUE K giants (KG) and SDSS/SEGUE blue horizontal branch (BHB) stars with accurate metallicities. The 3D kinematics are available through LAMOST and SDSS/SEGUE distances and radial velocities and {\it Gaia} DR2 proper motions. Two methods are used to estimate the enclosed mass: 3D spherical Jeans equation and Evans et al. tracer mass estimator (TME). We remove substructure via the Xue et al. method based on integrals of motion. We evaluate the uncertainties on our estimates due to random sampling noise, systematic distance errors, the adopted density profile, and non-virialization and non-spherical effects of the halo. The tracer density profile remains a limiting systematic in our mass estimates, although within these limits we find reasonable agreement across the different samples and the methods applied. Out to $\sim70$ and $\sim50$ kpc, the Jeans method yields total enclosed masses of $4.3\pm0.95$ (random) $\pm0.6$ (systematic) $\times10^{11}$ M$_\odot$ and $4.1\pm1.2$ (random) $\pm0.6$ (systematic) $\times10^{11}$ M$_\odot$ for the KG and BHB stars, respectively. For the KG and BHB samples we find a dark matter virial mass of $M_{200}=0.55^{+0.15}_{-0.11}$ (random) $\pm0.083$ (systematic) $\times10^{12}$ M$_\odot$ and $M_{200}=1.00^{+0.67}_{-0.33}$ (random) $\pm0.15$ (systematic) $\times10^{12}$ M$_\odot$, respectively., Comment: Accepted for publication in MNRAS. Comments welcomed

We present a study of the influence of the Gaia-Sausage-Enceladus (GSE) on the density shape of the Galactic stellar halo using 11624 K giants from the LAMOST survey. Every star is assigned a probability of being a member of the GSE based on its spherical velocities and metallicity by a Gaussian Mixture Model. We divide the stellar halo into two parts by the obtained probabilities, of which one is composed of the GSE members and defined as the GSE-related halo, and the other one is referred to as the GSE-removed halo. Using a non-parametric method, the radial number density profiles of the two stellar halos can be well described by a single power law with a variable flattening $q$ ($r = \sqrt{R^2+[(Z/q(r))]^2}, \nu = {\nu_0}r^{-\alpha}$). The index $\alpha$ is $4.92\pm0.12$ for the GSE-related halo and $4.25\pm0.14$ for the GSE-removed halo. Both the two stellar halos are vertically flattened at smaller radii but become more spherical at larger radii. We find that the GSE-related halo is less vertically flattened than the GSE-removed halo, and the difference of $q$ between the two stellar halos ranges from 0.07 to 0.15. However, after the consideration of the bias, it is thought to be within 0.08 at most of the radii. Finally, we compare our results with two Milky Way analogues which experience a significant major merger in the TNG50 simulation. The study of the two analogues also shows that the major merger-related stellar halo has a smaller ellipticity than the major merger-removed stellar halo., Comment: 19 Pages, 14 figures, Accepted by AJ

We explore the contribution of the Gaia-Sausage to the stellar halo of the Milky Way by making use of a Gaussian Mixture model (GMM) and applying it to halo star samples of LAMOST K giants, SEGUE K giants, and SDSS blue horizontal branch stars. The GMM divides the stellar halo into two parts, of which one represents a more metal-rich and highly radially biased component associated with an ancient, head-on collision referred to as the Gaia-Sausage, and the other one is a more metal-poor and isotropic halo. A symmetric bimodal Gaussian is used to describe the distribution of spherical velocity of the Gaia-Sausage, and we find that the mean absolute radial velocity of the two lobes decreases with Galactocentric radius. We find that the Gaia-Sausage contributes about $41\%-74\%$ of the inner (Galactocentric radius $r_\mathrm{gc} < 30$ kpc) stellar halo. The fraction of stars of the Gaia-Sausage starts to decline beyond $r_\mathrm{gc} \sim$ $25-30$ kpc, and the outer halo is found to be significantly less influenced by the Gaia-Sausage than the inner halo. After the removal of halo substructures found by integrals of motion, the contribution of the \textit{Gaia}-Sausage falls slightly within $r_\mathrm{gc} \sim$ 25 kpc, but is still as high as $30\%-63\%$. Finally, we select several possible Sausage-related substructures consisting of stars on highly eccentric orbits. The GMM/Sausage component agrees well with the selected substructure stars in their chemodynamical properties, which increases our confidence in the reliability of the GMM fits., Comment: accepted for publication in ApJ

We examine the vertical structure of Galactic disk kinematics over a Galactocentric radial distance range of $R=5-15$ $\rm{kpc}$ and up to $3$ $\rm{kpc}$ away from the Galactic plane, using the K-type giants surveyed by LAMOST. Based on robust measurements of three-dimensional velocity moments, a wobbly disk is detected in a phenomenological sense. An outflow dominates the radial motion of the inner disk, while in the outer disk there exist alternate outward and inward flows. The vertical bulk velocities is a combination of breathing and bending modes. A contraction-like breathing mode with amplitudes increasing with the distance to the plane and an upward bending mode dominate the vertical motion outside $R_0$, and there are reversed breathing mode and bending mode at $R

In this work, we present a systematic search for stellar groups in the Taurus field by applying DBSCAN algorithm to the data from Gaia DR2. We find 22 groups, consisting of 8 young groups (Groups 1-8) at ages of 2-4Myr and distances of ~130-170pc, 14 old groups (Groups 9-22) at ages of 8-49Myr and distances of ~110-210pc. We characterize the disk properties of group members and find 19 new disk-bearing stars, 8 of which are in the young groups and 11 others belong to the comparatively old groups at ages of 8-11 Myr. We characterize the accretion properties of the group members with H$\alpha$ emission line in their LAMOST spectra, and discover one source in Group 10 at an age of 10 Myr which still shows accretion activity. We investigate the kinematic relations among the old groups, and find that Group 9 is kinematically related to the known Taurus members and exclude any kinematic relations between Groups 10-22 and the known Taurus members., Comment: 36 pages, 17 figures and 2 tables, accepted for publication in ApJs

We characterize the kinematic and chemical properties of 589 Galactic Anticenter Substructure Stars (GASS) with K-/M- giants in Integrals-of-Motion space. These stars likely include members of previously identified substructures such as Monoceros, A13, and the Triangulum-Andromeda cloud (TriAnd). We show that these stars are on nearly circular orbits on both sides of the Galactic plane. We can see velocity($V_{Z}$) gradient along Y-axis especially for the south GASS members. Our GASS members have similar energy and angular momentum distributions to thin disk stars. Their location in [$\alpha$/M] vs. [M/H] space is more metal poor than typical thin disk stars, with [$\alpha$/M] \textbf{lower} than the thick disk. We infer that our GASS members are part of the outer metal-poor disk stars, and the outer-disk extends to 30 kpc. Considering the distance range and $\alpha$-abundance features, GASS could be formed after the thick disk was formed due to the molecular cloud density decreased in the outer disk where the SFR might be less efficient than the inner disk., Comment: 17 pages, 8 figures

We use a distribution function analysis to estimate the mass of the Milky Way out to 100 kpc using a large sample of halo stars. These stars are compiled from the literature, and the vast majority (~98%) have 6D phase-space information. We pay particular attention to systematic effects, such as the dynamical influence of the Large Magellanic Cloud (LMC), and the effect of unrelaxed substructure. The LMC biases the (pre-LMC infall) halo mass estimates towards higher values, while realistic stellar halos from cosmological simulations tend to underestimate the true halo mass. After applying our method to the Milky Way data we find a mass within 100 kpc of M(< 100 kpc) = 6.07 +/- 0.29 (stat.) +/- 1.21 (sys.) x 10^11 M_Sun. For this estimate, we have approximately corrected for the reflex motion induced by the LMC using the Erkal et al. model, which assumes a rigid potential for the LMC and MW. Furthermore, stars that likely belong to the Sagittarius stream are removed, and we include a 5% systematic bias, and a 20% systematic uncertainty based on our tests with cosmological simulations. Assuming the mass-concentration relation for Navarro-Frenk-White haloes, our mass estimate favours a total (pre-LMC infall) Milky Way mass of M_200c = 1.01 +/- 0.24 x 10^12 M_Sun, or (post-LMC infall) mass of M_200c = 1.16 +/- 0.24 x 10^12 M_Sun when a 1.5 x 10^11 M_Sun mass of a rigid LMC is included., Comment: 9 pages, 6 figures. MNRAS accepted

A wealth of recent studies have shown that the LMC is likely massive, with a halo mass $>10^{11} M_\odot$. One consequence of having such a nearby and massive neighbour is that the inner Milky Way is expected to be accelerated with respect to our Galaxy's outskirts (beyond $\sim 30$ kpc). In this work we compile a sample of $\sim 500$ stars with radial velocities in the distant stellar halo, $r_{\rm GC}> 50$ kpc, to test this hypothesis. These stars span a large fraction of the sky and thus give a global view of the stellar halo. We find that stars in the Southern hemisphere are on average blueshifted, while stars in the North are redshifted, consistent with the expected, mostly downwards acceleration of the inner halo due to the LMC. We compare these results with simulations and find the signal is consistent with the infall of a $1.5\times10^{11} M_\odot$ LMC. We cross-match our stellar sample with \textit{Gaia} DR2 and find that the mean proper motions are not yet precise enough to discern the LMC's effect. Our results show that the outer Milky Way is significantly out of equilibrium and that the LMC has a substantial effect on our Galaxy., Comment: 9 pages, 7 figures, submitted to MNRAS. Comments welcome!

We use K-giant stars selected from the LAMOST DR5 to study the variation of the rotational velocity of the galactic halo at different space positions. Modelling the rotational velocity distribution with both the halo and disk components, we find that the rotational velocity of the halo population decreases almost linearly with increasing vertical distance to the galactic disk plane, $Z$, at fixed galactocentric radius, $R$. The samples are separated into two parts with $6

We analyze the anisotropy profile of the Milky Way's smooth, diffuse stellar halo using SDSS/SEGUE blue horizontal branch stars and SDSS/SEGUE and LAMOST K giants. These intrinsically luminous stars allow us to probe the halo to approximately 100 kpc from the Galactic center. Line-of-sight velocities, distances, metallicities, and proper motions are available for all stars via SDSS/SEGUE, LAMOST, and Gaia, and we use these data to construct a full 7D set consisting of positions, space motions, and metallicity. We remove substructure from our samples using integrals of motion based on the method of Xue et al. We find radially dominated kinematic profiles with nearly constant anisotropy within 20 kpc, beyond which the anisotropy profile gently declines although remains radially dominated to the furthest extents of our sample. Independent of star type or substructure removal, the anisotropy depends on metallicity, such that the orbits of the stars become less radial with decreasing metallicity. For $-1.7<$ [Fe/H] $<-1$, the smooth, diffuse halo anisotropy profile begins to decline at Galactocentric distances $\sim20$ kpc, from $\beta\sim0.9$ to 0.7 for K giants and from $\beta\sim0.8$ to 0.1 for blue horizontal branch stars. For [Fe/H] $<-1.7$, the smooth, diffuse halo anisotropy remains constant along all distances with $0.2<\beta<0.7$ depending on the metallicity range probed, although independent on star type. These samples are ideal for estimating the total Galactic mass as they represent the virialized stellar halo system., Comment: 35 pages, 20 figures, 2 tables, accepted by ApJ, comments welcome

We characterize the kinematic and chemical properties of $\sim$3,000 Sagittarius (Sgr) stream stars, including K-giants, M-giants, and BHBs, select from SEGUE-2, LAMOST, and SDSS separately in Integrals-of-Motion space. The orbit of Sgr stream is quite clear from the velocity vector in $X$-$Z$ plane. Stars traced by K-giants and M-giants present the apogalacticon of trailing steam is $\sim$ 100 kpc. The metallicity distributions of Sgr K-, M-giants, and BHBs present that the M-giants are on average the most metal-rich population, followed by K-giants and BHBs. All of the K-, M-giants, and BHBs indicate that the trailing arm is on average more metal-rich than leading arm, and the K-giants show that the Sgr debris is the most metal-poor part. The $\alpha$-abundance of Sgr stars exhibits a similar trend with the Galactic halo stars at lower metallicity ([Fe/H] $<\sim$ $-$1.0 dex), and then evolve down to lower [$\alpha$/Fe] than disk stars at higher metallicity, which is close to the evolution pattern of $\alpha$-element of Milky Way dwarf galaxies. We find $V_Y$ and metallicity of K-giants have gradients along the direction of line-of-sight from the Galactic center in $X$-$Z$ plane, and the K-giants show that $V_Y$ increases with metallicity at [Fe/H] $>\sim-$1.5 dex. After dividing the Sgr stream into bright and faint stream according to their locations in equatorial coordinate, the K-giants and BHBs show that the bright and faint stream present different $V_Y$ and metallicities, the bright stream is on average higher in $V_Y$ and metallicity than the faint stream., Comment: 26 pages, 15 figures

We construct a large halo K-giant sample by combining the positions, distances, radial velocities, and metallicities of over 13,000 LAMOST DR5 halo K giants with the Gaia DR2 proper motions, which covers a Galactocentric distance range of 5-120 kpc. Using a position-velocity clustering estimator (the 6Distance), we statistically quantify the presence of position-velocity substructure at high significance: K giants have more close pairs in position-velocity space than a smooth stellar halo. We find that the amount of substructure in the halo increases with increasing distance and metallicity. With a percolation algorithm named friends-of-friends (FoF) to identify groups, we identify members belonging to Sagittarius (Sgr) Streams, Monoceros Ring, Virgo overdensity, Hercules-Aquila Cloud, Orphan Streams and other unknown substructures and find that the Sgr streams account for a large part of grouped stars beyond 20 kpc and enhance the increase of substructure with distance and metallicity. For the first time, we identify spectroscopic members of Monoceros Ring in the south and north Galactic hemisphere, which presents a rotation of about 185 km s^{-1} and mean metallicity is -0.66 dex.

We study the radial metallicity gradient $\Delta[M/H]/\Delta R_g$ as a function of [Mg/Fe] and $|Z|$ with the help of a guiding radius based on the Apache Point Observatory Galactic Evolution Experiment and Gaia and then analyze the radial migration effect on the radial metallicity gradient and metallicity-rotation gradient between the Galactic thin and thick disks. The derived trend of gradient $\Delta[M/H]/\Delta R_g$ versus [Mg/Fe] shows a transition at [Mg/Fe] $\sim 0.18$ dex, below which the gradient is negative and varies a little as [Mg/Fe] increases; however, it changes sharply in [Mg/Fe] ranges of 0.16-0.18, above which the gradient increases linearly with increasing [Mg/Fe], being a positive value at [Mg/Fe]$\gtrsim 0.22$ dex. These positive gradients in the high-[Mg/Fe] populations are found at $|Z| < 0.8$ kpc, and there are nearly no gradients toward higher $|Z|$. By comparing the metallicity distributions, the radial metallicity gradients $\Delta[M/H]/\Delta R$ and the metallicity-rotation gradients between the total sample and $|R-R_g|<2$ kpc subsample (or $|R-R_g|>2$ kpc subsample), we find that, for the thick disk, blurring flattens the gradient $\Delta[M/H]/\Delta R$ and favors metal-poor high-eccentricity stars. These stars are responsible for the measured positive metallicity-rotation gradient of the thick disk., Comment: Published in ApJ, including 8 pages, 8 figures

The anisotropy parameter $\beta$ characterizes the extent to which orbits in stellar systems are predominantly radial or tangential, and is likely to constrain, for the stellar halo of the Milky Way, scenarios for its formation and evolution. We have measured the anisotropy $\beta$ as a function of Galactocentric radius from $5-100$ kpc for over 8600 metal poor ([Fe/H] $<-1.3$) halo K giants from the LAMOST catalog with line-of-sight velocities and distances, matched to proper motions from the second $Gaia$ data release. We construct full 6-D positions and velocities for the K giants to directly measure the 3 components of the velocity dispersion $(\sigma_r, \sigma_\theta, \sigma_\phi)$ (in spherical coordinates). We find that the orbits in the halo are radial over our full Galactocentric distance range reaching over 100 kpc. The anisotropy remains remarkably unchanged with Galactocentric radius from approximately 5 to 25 kpc, with an amplitude that depends on the metallicity of the stars, dropping from $\beta \approx 0.9$ for $-1.8 \leq$ [Fe/H] $< -1.3$ (for the bulk of the stars) to $\beta \approx 0.6$ for the lowest metallicities ([Fe/H] $< -1.8$). Considering our sample as a whole, $\beta\approx0.8$ and, beyond 25 kpc, the orbits gradually become less radial and anisotropy decreases to $\beta<0.3$ past 100 kpc. Within 8 kpc, $\beta<0.8$. The measurement of anisotropy is affected by substructure and streams, particularly beyond a Galactocentric distance of approximately 25 kpc, where the Sagittarius stream is prominent in the data. These results are complimentary to recent analysis of simulations by Loebman et al. and of SDSS/$Gaia$ DR1 data by Belokurov et al., Comment: 9 pages, 5 figures, published version