Search

We aim to study the phase spiral in the Milky Way (MW) with Gaia DR3. We used an edge detection algorithm to find the border of the phase spiral, allowing us to robustly quantify its shape at different positions and for different selections. We calculated the time of onset of the phase-mixing by determining the different turns of the phase spiral and using the vertical frequencies from commonly used MW potential models. We find that the phase spiral extends down to $-1.2$ kpc in height below the plane (about 3 to 5 scale heights of the thin disc) and beyond $\pm 50$ km/s in $V_Z$. We see a secondary branch mostly at positive vertical velocities when coloured by azimuthal velocity and in the counts projection. We also find complex variations of the phase spirals with angular momentum and azimuth. All these possibly provide evidence of multiple perturbations (from different times or from different perturbers) and/or of the complexity of the phase mixing process. We detect the phase spiral from 6 to 11 kpc from the Galactic centre and find signatures of vertical asymmetries 1-2 kpc beyond this range. We measure small but clear variations with azimuth. When we determine the phase mixing times from the phase spiral at different angular momenta and using the different spiral turns (at different $Z$) we obtain inconsistent times with systematic differences (times increasing with $|L_Z|$ and with $|Z|$). Our determinations are mostly in the range of [0.3-0.9] Gyr, with an average of 0.5 Gyr. The inconsistencies do not change when using different usual potential models, different stellar distances or frequencies for different kinetic temperatures. They could stem from the inconsistency of potential models with the true MW, and from too simple modelling, in particular neglecting self-gravity, not considering the multiple perturbations and the interference with other processes., Comment: version after proofs corrections

The density of dark matter near the Sun is important for experiments hunting for dark matter particles in the laboratory, and for constraining the local shape of the Milky Way's dark matter halo. Estimates to date have typically assumed that the Milky Way's stellar disc is axisymmetric and in a steady-state. Yet the Milky Way disc is neither, exhibiting prominent spiral arms and a bar, and vertical and radial oscillations. We assess the impact of these assumptions on determinations of the local dark matter density by applying a free-form, steady-state, Jeans method to two different N-body simulations of Milky Way-like galaxies. In one, the galaxy has experienced an ancient major merger, similar to the hypothesized Gaia-Sausage-Enceladus; in the other, the galaxy is perturbed more recently by the repeated passage and slow merger of a Sagittarius-like dwarf galaxy. We assess the impact of each of the terms in the Jeans-Poisson equations on our ability to correctly extract the local dark matter density from the simulated data. We find that common approximations employed in the literature - axisymmetry and a locally flat rotation curve - can lead to significant systematic errors of up to a factor ~1.5 in the recovered surface mass density ~2kpc above the disc plane, implying a fractional error on the local dark matter density of order unity. However, once we add in the tilt term and the rotation curve term in our models, we obtain an unbiased estimate, consistent with the true value within our 95% confidence intervals for realistic 20% uncertainties on the baryonic surface density of the disc. Other terms - the axial tilt, 2:nd Poisson and time dependent terms - contribute less than 10% to the local dark matter density (given current data) and can be safely neglected for now. In the future, as more data become available, these terms will need to be included in the analysis., Comment: 15 pages, 16 figues. Submitted to MNRAS. Comments welcome!

The astrometric sample of Gaia allows us to study the outermost Galactic disc, the halo and their interface. It is precisely at the very edge of the disc where the effects of external perturbations are expected to be the most noticeable. Our goal is to detect the kinematic substructure present in the halo and at the edge of the Milky Way (MW) disc, and provide observational constraints on their phase-space distribution. We download, one HEALpix at a time, the proper motion histogram of distant stars, to which we apply a Wavelet Transformation to reveal the significant overdensities. We then analyse the large coherent structures that appear in the sky. We reveal a sharp yet complex anticentre dominated by Monoceros (MNC) and the Anticentre Stream (ACS) in the north, which we find with an intensity comparable to the Magellanic clouds and the Sagittarius stream, and by MNC south and TriAnd at negative latitudes. Our method allows us to perform a morphological analysis of MNC and ACS, both spanning more than 100$^\circ$ in longitude, and to provide a high purity sample of giants with which we track MNC down to latitudes as low as $\sim$5$^\circ$. Their colour-magnitude diagram is consistent with extended structures at a distance of $\sim$10-11 kpc originated in the disc, with a very low ratio of RR Lyrae over M giants, and kinematics compatible with the rotation curve at those distances or only slightly slower. We present a precise characterisation of MNC and ACS, two previously known structures that our method reveals naturally, allowing us to detect them without limiting ourselves to a particular stellar type and, for the first time, using only kinematics. Our results allow future studies to model their chemo-dynamics and evolution, thus constraining some of the most influential processes that shaped the MW., Comment: Submitted to Astronomy and Astrophysics. Online material available upon acceptance

The recent discovery of a spiral feature in $Z-V_Z$ phase plane in the solar neighborhood implies that the Galactic disk has been remarkably affected by a dwarf galaxy passing through it some hundreds of millions of years ago. Using 429,500 LAMOST K giants stars, we show that the spiral feature exits not only in the solar vicinity; it also extends to about 15 kpc from the Galactic center, and then disappears beyond this radius. Moreover, we find that when the spiral features in a plot of $V_\phi$ as a function of position in the $Z-V_Z$ plane, at various Galactocentric radii, are re-mapped to $R-Z$ plane, the spiral can explain well the observed asymmetric velocity substructures. This is evidence that the phase spiral features are the same as the bulk motions found in previous as well as this work. Test-particle simulations and N-body simulations show that an encounter with a dwarf galaxy a few hundred million years ago will induce a perturbation in the Galactic disk. In addition, we find that the last impact of Sgr dSph can also contribute to the flare. As a consequence of the encounter, the distribution function of disk stars at a large range of radii is imprinted by the gravitational perturbation., Comment: 31pages, 23figures, ApJ accepted

We study the late-time evolution of the central regions of two Milky Way-like simulations of galaxies formed in a cosmological context, one hosting a fast bar and the other a slow one. We find that bar length, R_b, measurements fluctuate on a dynamical timescale by up to 100%, depending on the spiral structure strength and measurement threshold. The bar amplitude oscillates by about 15%, correlating with R_b. The Tremaine-Weinberg-method estimates of the bars' instantaneous pattern speeds show variations around the mean of up to ~20%, typically anti-correlating with the bar length and strength. Through power spectrum analyses, we establish that these bar pulsations, with a period in the range ~60-200 Myr, result from its interaction with multiple spiral modes, which are coupled with the bar. Because of the presence of odd spiral modes, the two bar halves typically do not connect at exactly the same time to a spiral arm, and their individual lengths can be significantly offset. We estimated that in about 50% of bar measurements in Milky Way-mass external galaxies, the bar lengths of SBab type galaxies are overestimated by ~15% and those of SBbc types by ~55%. Consequently, bars longer than their corotation radius reported in the literature, dubbed "ultra-fast bars", may simply correspond to the largest biases. Given that the Scutum-Centaurus arm is likely connected to the near half of the Milky Way bar, recent direct measurements may be overestimating its length by 1-1.5 kpc, while its present pattern speed may be 5-10 km/s/kpc smaller than its time-averaged value., Comment: 25 p., 19 figures, Accepted to MNRAS

We use Gaia DR2 astrometric and line-of-sight velocity information combined with two sets of distances obtained with a Bayesian inference method to study the 3D velocity distribution in the Milky Way disc. We search for variations in all Galactocentric cylindrical velocity components ($V_{\phi}$, $V_R$ and $V_z$) with Galactic radius, azimuth, and distance from the disc mid-plane. We confirm recent work showing that bulk vertical motions in the $R\text{-}z$ plane are consistent with a combination of breathing and bending modes. In the $x\text{-}y$ plane, we show that, although the amplitudes change, the structure produced by these modes is mostly invariant as a function of distance from the plane. Comparing to two different Galactic disc models, we demonstrate that the observed patterns can drastically change in short time intervals, showing the complexity of understanding the origin of vertical perturbations. A strong radial $V_R$ gradient was identified in the inner disc, transitioning smoothly from $16$ km s$^{-1}$ kpc$^{-1}$ at an azimuth of $30^\circ<\phi<45^\circ$ ahead of the Sun-Galactic centre line, to $-16$ km s$^{-1}$ kpc$^{-1}$ at an azimuth of $-45^\circ<\phi<-30^\circ$ lagging the solar azimuth. We use a simulation with no significant recent mergers to show that exactly the opposite trend is expected from a barred potential, but overestimated distances can flip this trend to match the data. Alternatively, using an $N$-body simulation of the Sagittarius dwarf-Milky Way interaction, we demonstrate that a major recent perturbation is necessary to reproduce the observations. Such an impact may have strongly perturbed the existing bar or even triggered its formation in the last $1\text{-}2$ Gyr., Comment: 17 p., 13 fig., accepted MNRAS

We use astrometry, broad-band photometry and variability information from the Data Release 2 of ESA's Gaia mission (GDR2) to identify members of the Orphan Stream (OS) across the whole sky. The stream is traced above and below the celestial equator and in both Galactic hemispheres, thus increasing its visible length to ~ 210 degrees equivalent to ~150 kpc in physical extent. Taking advantage of the large number of RR Lyrae stars in the OS, we extract accurate distances and proper motions across the entire stretch of the tidal debris studied. As delineated by the GDR2 RR Lyrae, the stream exhibits two prominent twists in its shape on the sky which are accompanied by changes in the tangential motion. We complement the RR Lyrae maps with those created using GDR2 Red Giants and the DECam Legacy Survey Main Sequence Turn-Off stars. The behavior of the OS track on the sky is consistent across all three tracers employed. We detect a strong non-zero motion in the across-stream direction for a substantial portion of the stream. Such a misalignment between the debris track and the streaming velocity cannot be reproduced in a static gravitational potential and signals an interaction with a massive perturber., Comment: 18 pages; accepted to MNRAS

In a companion paper by Koposov et al., RR Lyrae from \textit{Gaia} Data Release 2 are used to demonstrate that stars in the Orphan stream have velocity vectors significantly misaligned with the stream track, suggesting that it has received a large gravitational perturbation from a satellite of the Milky Way. We argue that such a mismatch cannot arise due to any realistic static Milky Way potential and then explore the perturbative effects of the Large Magellanic Cloud (LMC). We find that the LMC can produce precisely the observed motion-track mismatch and we therefore use the Orphan stream to measure the mass of the Cloud. We simultaneously fit the Milky Way and LMC potentials and infer that a total LMC mass of $1.38^{+0.27}_{-0.24} \times10^{11}\,\rm{M_\odot}$ is required to bend the Orphan Stream, showing for the first time that the LMC has a large and measurable effect on structures orbiting the Milky Way. This has far-reaching consequences for any technique which assumes that tracers are orbiting a static Milky Way. Furthermore, we measure the Milky Way mass within 50 kpc to be $3.80^{+0.14}_{-0.11}\times10^{11} M_\odot$. Finally, we use these results to predict that, due to the reflex motion of the Milky Way in response to the LMC, the outskirts of the Milky Way's stellar halo should exhibit a bulk, upwards motion., Comment: 17 pages, 11 figures. Updated to version accepted to MNRAS after minor revision

We build a sample of distant (D > 80 kpc) stellar halo stars with measured radial velocities. Faint ($20 < g <22$) candidate blue horizontal branch (BHB) stars were selected using the deep, but wide, multi-epoch Sloan Digital Sky Survey photometry. Follow-up spectroscopy for these A-type stars was performed using the VLT-FORS2 instrument. We classify stars according to their Balmer line profiles, and find 7 are bona fide BHB stars and 31 are blue stragglers (BS). Owing to the magnitude range of our sample, even the intrinsically fainter BS stars can reach out to D ~ 90 kpc. We complement this sample of A-type stars with intrinsically brighter, intermediate-age, asymptotic giant branch stars. A set of 4 distant cool carbon stars is compiled from the literature and we perform spectroscopic follow-up on a further 4 N-type carbon stars using the WHT-ISIS instrument. Altogether, this provides us with the largest sample to date of individual star tracers out to r ~ 150 kpc. We find that the radial velocity dispersion of these tracers falls rapidly at large distances and is surprisingly cold (sigma_r ~ 50-60 km/s) between 100-150 kpc. Relating the measured radial velocities to the mass of the Milky Way requires knowledge of the (unknown) tracer density profile and anisotropy at these distances. Nonetheless, by assuming the stellar halo stars between 50-150 kpc have a moderate density fall-off (with power-law slope alpha < 5) and are on radial orbits (sigma^2_t/sigma^2_r < 1), we infer that the mass within 150 kpc is less than 10^12 M_solar and suggest it probably lies in the range (5-10) x 10^11 M_solar. We discuss the implications of such a low mass for the Milky Way., Comment: 15 pages, 12 figures. Accepted for publication in MNRAS