Combining the precise parallaxes and optical photometry delivered by Gaia’s second data release with the photometric catalogues of Pan-STARRS1, 2MASS, and AllWISE, we derived Bayesian stellar ...parameters, distances, and extinctions for 265 million of the 285 million objects brighter than G = 18. Because of the wide wavelength range used, our results substantially improve the accuracy and precision of previous extinction and effective temperature estimates. After cleaning our results for both unreliable input and output data, we retain 137 million stars, for which we achieve a median precision of 5% in distance, 0.20 mag in V-band extinction, and 245 K in effective temperature for G ≤ 14, degrading towards fainter magnitudes (12%, 0.20 mag, and 245 K at G = 16; 16%, 0.23 mag, and 260 K at G = 17, respectively). We find a very good agreement with the asteroseismic surface gravities and distances of 7000 stars in the Kepler, K2-C3, and K2-C6 fields, with stellar parameters from the APOGEE survey, and with distances to star clusters. Our results are available through the ADQL query interface of the Gaia mirror at the Leibniz-Institut für Astrophysik Potsdam (gaia.aip.de) and as binary tables at data.aip.de. As a first application, we provide distance- and extinction-corrected colour-magnitude diagrams, extinction maps as a function of distance, and extensive density maps. These demonstrate the potential of our value-added dataset for mapping the three-dimensional structure of our Galaxy. In particular, we see a clear manifestation of the Galactic bar in the stellar density distributions, an observation that can almost be considered direct imaging of the Galactic bar.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
We measure the escape speed curve of the Milky Way based on the analysis of the velocity distribution of ~2850 counter-rotating halo stars from the Gaia Data Release 2. The distances were estimated ...through the StarHorse code, and only stars with distance errors smaller than 10% were used in the study. The escape speed curve is measured at Galactocentric radii ranging from ~5 kpc to ~10.5 kpc. The local Galactic escape at the Sun’s position is estimated to be ve(r⊙) = 580 ± 63 km s−1, and it rises towards the Galactic centre. Defined as the minimum speed required to reach three virial radii, our estimate of the escape speed as a function of radius implies for a Navarro–Frenk–White profile and local circular velocity of 240 km s−1 a dark matter mass M200 = 1.28−0.50+0.68 × 1012 M⊙ M 200 = 1.28 − 0.50 + 0.68 × 10 12 M ⊙ $ M_{200}=1.28_{-0.50}^{+0.68}\times10^{12}M_\odot $ and a high concentration c200 = 11.09−1.79+2.94 c 200 = 11.09 − 1.79 + 2.94 $ c_{200}=11.09_{-1.79}^{+2.94} $ . Assuming the mass-concentration relation of ΛCDM, we obtain M200 = 1.55−0.51+0.64 × 1012 M⊙ M 200 = 1.55 − 0.51 + 0.64 × 10 12 M ⊙ $ M_{200}=1.55_{-0.51}^{+0.64}\times10^{12}M_\odot $ and c200 = 7.93−0.27+0.33 c 200 = 7.93 − 0.27 + 0.33 $ c_{200}=7.93_{-0.27}^{+0.33} $ for a local circular velocity of 228 km s−1.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
The second data release of the Gaia mission has revealed a very rich structure in local velocity space. In terms of in-plane motions, this rich structure is also seen as multiple ridges in the ...actions of the axisymmetric background potential of the Galaxy. These ridges are probably related to a combination of effects from ongoing phase-mixing and resonances from the spiral arms and the bar. We have recently developed a method for capturing the behaviour of the stellar phase-space distribution function at a resonance by re-expressing it in terms of a new set of canonical actions and angles variables valid in the resonant region. Here, by properly treating the distribution function at resonances, and by using a realistic model for a slowly rotating large Galactic bar with pattern speed Ωb = 39 km s−1 kpc−1, we show that no fewer than six ridges in local action space can be related to resonances with the bar. Two of these ridges at low angular momentum correspond to the corotation resonance, and can be associated with the Hercules moving group in local velocity space. Another ridge at high angular momentum corresponds to the outer Lindblad resonance, and can tentatively be associated with the velocity structure seen as an arch at high azimuthal velocities in Gaia data. The other ridges are associated with the 3:1, 4:1, and 6:1 resonances. The last can be associated with the so-called “horn” of the local velocity distribution. While it is clear that effects from spiral arms and incomplete phase-mixing related to external perturbations also play a role in shaping the complex kinematics revealed by Gaia data, the present work demonstrates that, contrary to common misconceptions, the bar alone can create multiple prominent ridges in velocity and action space.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
ABSTRACT
We study the late-time evolution of the central regions of two Milky Way (MW)-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, Rb, measurements fluctuate on a dynamical time-scale by up to 100 per cent, depending on the spiral structure strength and measurement threshold. The bar amplitude oscillates by about 15 per cent, correlating with Rb. The Tremaine–Weinberg method estimates of the bars’ instantaneous pattern speeds show variations around the mean of up to $\sim \!20{{\ \rm per\ cent}}$, typically anticorrelating 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 per cent of bar measurements in MW-mass external galaxies, the bar lengths of SBab-type galaxies are overestimated by $\sim \!15{{\ \rm per\ cent}}$ and those of SBbc types by $\sim \!55{{\ \rm per\ cent}}$. Consequently, bars longer than their corotation radius reported in the literature, dubbed ‘ultrafast bars’, may simply correspond to the largest biases. Given that the Scutum–Centaurus arm is likely connected to the near half of the MW bar, recent direct measurements may be overestimating its length by 1–1.5 kpc, while its present pattern speed may be 5–10 $\rm km\ s^{-1}\ kpc^{-1}$ smaller than its time-averaged value.
We use data from the Radial Velocity Experiment (RAVE) and the Tycho-Gaia astrometric solution (TGAS) catalogue to compute the velocity fields yielded by the radial (V-R), azimuthal (V-phi), and ...vertical (V-z) components of associated Galactocentric velocity. We search in particular for variation in all three velocity components with distance above and below the disc midplane, as well as how each component of V-z (line-of-sight and tangential velocity projections) modifies the obtained vertical structure. To study the dependence of velocity on proper motion and distance, we use two main samples: a RAVE sample including proper motions from the Tycho-2, PPMXL, and UCAC4 catalogues, and a RAVE-TGAS sample with inferred distances and proper motions from the TGAS and UCAC5 catalogues. In both samples, we identify asymmetries in V-R and V-z. Below the plane, we find the largest radial gradient to be partial derivative V-R/partial derivative R = -7.01 +/- 0.61 km s(-1) kpc(-1), in agreement with recent studies. Above the plane, we find a similar gradient with partial derivative V-R/partial derivative R = -9.42 +/- 1.77 km s(-1) kpc(-1). By comparing our results with previous studies, we find that the structure in V-z is strongly dependent on the adopted proper motions. Using the Galaxia Milky Way model, we demonstrate that distance uncertainties can create artificial wave-like patterns. In contrast to previous suggestions of a breathing mode seen in RAVE data, our results support a combination of bending and breathing modes, likely generated by a combination of external or internal and external mechanisms.
Aims.
We infer the gravitational potential of the Galactic disk by analysing the phase-space densities of 120 stellar samples in 40 spatially separate sub-regions of the solar neighbourhood, using
...Gaia
’s second data release (DR2), in order to quantify spatially dependent systematic effects that bias this type of measurement.
Methods.
The gravitational potential was inferred under the assumption of a steady state in the framework of a Bayesian hierarchical model. We performed a joint fit of our stellar tracers’ three-dimensional velocity distribution, while fully accounting for the astrometric uncertainties of all stars as well as dust extinction, and we also masked angular areas of known open clusters. The inferred gravitational potential is compared, post-inference, to a model for the baryonic matter and halo dark matter components.
Results.
We see an unexpected but clear trend for all 40 spatially separate sub-regions: Compared to the potential derived from the baryonic model, the inferred gravitational potential is significantly steeper close to the Galactic mid-plane (≲60 pc), but flattens such that the two agree well at greater distances (∼400 pc). The inferred potential implies a total matter density distribution that is highly concentrated to the Galactic mid-plane and decays quickly with height. We see a dependence on the Galactic radius that is consistent with a disk scale length of a few kiloparsecs. Apart from this, there are discrepancies between stellar samples, implying spatially dependent systematic effects which are, at least in part, explained by substructures in the phase-space distributions.
Conclusions.
In terms of the inferred matter density distribution, the very low matter density that is inferred at greater heights (≳300 pc) is inconsistent with the observed scale height and matter distribution of the stellar disk, which cannot be explained by a misunderstood density of cold gas or other hidden mass. Our interpretation is that these results must be biased by a time-varying phase-space structure, possibly a breathing mode, that is large enough to affect all stellar samples in the same manner.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
Context.
The structure and dynamics of the central bar of the Milky Way (MW) are still under debate whilst being fundamental ingredients for the evolution of our Galaxy. The recent
Gaia
DR3 offers an ...unprecedented detailed view of the 6D phase space of the MW, allowing for a better understanding of the complex imprints of the bar on the phase space.
Aims.
We aim to identify and characterise the dynamical moving groups across the MW disc, and use their large-scale distribution to help constrain the properties of the Galactic bar.
Methods.
We used 1D wavelet transforms of the azimuthal velocity (
V
ϕ
) distribution in bins of radial velocity to robustly detect the kinematic substructure in the
Gaia
DR3 catalogue. We then connected these structures across the disc to measure the azimuthal (
ϕ
) and radial (
R
) gradients of
V
ϕ
of the moving groups. We simulated thousands of perturbed distribution functions using backward integration, sweeping a large portion of parameter space of feasible Galaxy models that include a bar, in order to compare them with the data and to explore and quantify the degeneracies.
Results.
The radial gradient of the Hercules moving group (∂
V
ϕ
/∂
R
= 28.1 ± 2.8 km s
−1
kpc
−1
) cannot be reproduced by our simple models of the Galaxy that show much larger slopes both for a fast and a slow bar. This suggests the need for more complex dynamics (e.g. a different bar potential, spiral arms, a slowing bar, a complex circular velocity curve, external perturbations, etc.). We measured an azimuthal gradient for Hercules of ∂
V
ϕ
/∂
ϕ
= −0.63 ± 0.13 km s
−1
deg
−1
and find that it is compatible with both the slow and fast bar models. Our analysis points out that in using this type of analysis, at least two moving groups are needed to start breaking the degeneracies.
Conclusions.
We conclude that it is not sufficient for a model to replicate the local velocity distribution; it must also capture its larger-scale variations. The accurate quantification of the gradients, especially in the azimuthal direction, will be key for the understanding of the dynamics governing the disc.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
The second data release of the Gaia mission has revealed, in stellar velocity and action space, multiple ridges, the exact origin of which is still debated. Recently, we demonstrated that a large ...Galactic bar with pattern speed 39 km s−1 kpc−1 creates most of the observed ridges. Among these ridges, the Hercules moving group would then be associated with orbits trapped at the co-rotation resonance of the bar. Here we show that a distinctive prediction of such a model is that the angular momentum of Hercules at the Sun’s radius must significantly decrease with increasing Galactocentric azimuth (i.e. when getting closer to the major axis of the bar). We show that this dependence of the angular momentum of trapped orbits on the azimuth on the other hand does not happen close to the outer Lindblad resonance of a faster bar, unless the orbital distribution is still far from phase-mixed, namely for a bar perturbation younger than ∼2 Gyr. Using Gaia DR2 and Bayesian distances from the StarHorse code, and tracing the average Galactocentric radial velocity as a function of angular momentum and azimuth, we show that the Hercules angular momentum changes significantly with azimuth as expected for the co-rotation resonance of a dynamically old large bar.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
Kinematics with Gaia DR2: the force of a dwarf Carrillo, I; Minchev, I; Steinmetz, M ...
Monthly notices of the Royal Astronomical Society,
11/2019, Volume:
490, Issue:
1
Journal Article
Peer reviewed
Open access
ABSTRACT
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ϕ, VR, and Vz) with Galactic radius, azimuth, and distance from the disc mid-plane. We confirm recent work showing that bulk vertical motions in the R–z plane are consistent with a combination of breathing and bending modes. In the x–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 VR gradient was identified in the inner disc, transitioning smoothly from 16 km s−1 kpc−1 at an azimuth of 30° < ϕ < 45° ahead of the Sun-Galactic centre line to −16 km s−1 kpc−1 at an azimuth of −45° < ϕ < −30° 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–2 Gyr.
Context.
Data-driven methods play an increasingly important role in the field of astrophysics. In the context of large spectroscopic surveys of stars, data-driven methods are key in deducing physical ...parameters for millions of spectra in a short time. Convolutional neural networks (CNNs) enable us to connect observables (e.g. spectra, stellar magnitudes) to physical properties (atmospheric parameters, chemical abundances, or labels in general).
Aims.
We test whether it is possible to transfer the labels derived from a high-resolution stellar survey to intermediate-resolution spectra of another survey by using a CNN.
Methods.
We trained a CNN, adopting stellar atmospheric parameters and chemical abundances from APOGEE DR16 (resolution
R
= 22 500) data as training set labels. As input, we used parts of the intermediate-resolution RAVE DR6 spectra (
R
∼ 7500) overlapping with the APOGEE DR16 data as well as broad-band ALL_WISE and 2MASS photometry, together with
Gaia
DR2 photometry and parallaxes.
Results.
We derived precise atmospheric parameters
T
eff
, log(
g
), and M/H, along with the chemical abundances of Fe/H,
α
/M, Mg/Fe, Si/Fe, Al/Fe, and Ni/Fe for 420 165 RAVE spectra. The precision typically amounts to 60 K in
T
eff
, 0.06 in log(
g
) and 0.02−0.04 dex for individual chemical abundances. Incorporating photometry and astrometry as additional constraints substantially improves the results in terms of the accuracy and precision of the derived labels, as long as we operate in those parts of the parameter space that are well-covered by the training sample. Scientific validation confirms the robustness of the CNN results. We provide a catalogue of CNN-trained atmospheric parameters and abundances along with their uncertainties for 420 165 stars in the RAVE survey.
Conclusions.
CNN-based methods provide a powerful way to combine spectroscopic, photometric, and astrometric data without the need to apply any priors in the form of stellar evolutionary models. The developed procedure can extend the scientific output of RAVE spectra beyond DR6 to ongoing and planned surveys such as
Gaia
RVS, 4MOST, and WEAVE. We call on the community to place a particular collective emphasis and on efforts to create unbiased training samples for such future spectroscopic surveys.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK