Open cluster kinematics with Gaia DR2 Soubiran, C.; Cantat-Gaudin, T.; Romero-Gómez, M. ...
Astronomy and astrophysics (Berlin),
11/2018, Letnik:
619
Journal Article
Recenzirano
Odprti dostop
Context. Open clusters are very good tracers of the evolution of the Galactic disc. Thanks to Gaia, their kinematics can be investigated with an unprecedented precision and accuracy. Aims. The ...distribution of open clusters in the 6D phase space is revisited with Gaia DR2. Methods. The weighted mean radial velocity of open clusters was determined, using the most probable members available from a previous astrometric investigation that also provided mean parallaxes and proper motions. Those parameters, all derived from Gaia DR2 only, were combined to provide the 6D phase-space information of 861 clusters. The velocity distribution of nearby clusters was investigated, as well as the spatial and velocity distributions of the whole sample as a function of age. A high-quality subsample was used to investigate some possible pairs and groups of clusters sharing the same Galactic position and velocity. Results. For the high-quality sample of 406 clusters, the median uncertainty of the weighted mean radial velocity is 0.5 km s−1. The accuracy, assessed by comparison to ground-based high-resolution spectroscopy, is better than 1 km s−1. Open clusters nicely follow the velocity distribution of field stars in the close solar neighbourhood as previously revealed by Gaia DR2. As expected, the vertical distribution of young clusters is very flat, but the novelty is the high precision to which this can be seen. The dispersion of vertical velocities of young clusters is at the level of 5 km s−1. Clusters older than 1 Gyr span distances to the Galactic plane of up to 1 kpc with a vertical velocity dispersion of 14 km s−1, typical of the thin disc. Five pairs of clusters and one group with five members might be physically related. Other binary candidates that have been identified previously are found to be chance alignments.
We study the spiral arm influence on the solar neighbourhood stellar kinematics. As the nature of the Milky Way (MW) spiral arms is not completely determined, we study two models: the Tight-Winding ...Approximation (TWA) model, which represents a local approximation, and a model with self-consistent material arms named sPiral arms potEntial foRmed by obLAte Spheroids (PERLAS). This is a mass distribution with more abrupt gravitational forces. We perform test particle simulations after tuning the two models to the observational range for the MW spiral arm properties. We find that some of the currently observed MW spiral arm properties are not in obvious agreement with the TWA model. We explore the effects of the arm properties and find that a significant region of the allowed parameter space favours the appearance of kinematic groups. The velocity distribution is mostly sensitive to the relative spiral arm phase and pattern speed. In all cases the arms induce strong kinematic imprints for pattern speeds around 17 km s−1 kpc−1 (close to the 4:1 inner resonance) but no substructure is induced close to corotation. The groups change significantly if one moves only ∼0.6 kpc in galactocentric radius, but ∼2 kpc in azimuth. The appearance time of each group is different, ranging from 0 to more than 1 Gyr. Recent spiral arms can produce strong kinematic structures. The stellar response to the two potential models is significantly different near the Sun, both in density and in kinematics. The PERLAS model triggers more substructure for a larger range of pattern speed values. The kinematic groups can be used to reduce the current uncertainty about the MW spiral structure and to test whether this follows the TWA. However, groups such as the observed ones in the solar vicinity can be reproduced by different parameter combinations. Data from velocity distributions at larger distances are needed for a definitive constraint.
ABSTRACT
We study the vertical perturbations in the galactic disc of the Milky Way-size high-resolution hydrodynamical cosmological simulation named GARROTXA. We detect phase spirals in the vertical ...projection Z − VZ of disc’s stellar particles for the first time in this type of simulations. Qualitatively similar structures were detected in the recent Gaia data, and their origin is still under study. In our model, the spiral-like structures in the phase space are present in a wide range of times and locations across the disc. By accounting for an evolving mix of stellar populations, we observe that, as seen in the data, the phase spirals are better observed in the range of younger-intermediate star particles. We measure the intensity of the spiral with a Fourier decomposition and find that these structures appear stronger near satellite pericentres. Current dynamical models of the phase spiral considering a single perturber required a mass at least of the order of 1010 M⊙, but all three of our satellites have masses of the order of ∼108 M⊙. We suggest that there are other mechanisms at play which appear naturally in our model such as the physics of gas, collective effect of multiple perturbers, and a dynamically cold population that is continuously renovated by the star formation. Complementing collisionless isolated N-body models with the use of fully cosmological simulations with enough resolution can provide new insights into the nature/origin of the phase spiral.
Context.
The so-called action variables are specific functions of the positions and velocities that remain constant along the stellar orbit. The astrometry provided by
Gaia
Early Data Release 3 ...(EDR3), combined with the velocities inferred from the Radial Velocity Spectrograph (RVS) spectra of
Gaia
DR3, allows for the estimation of these actions for the largest volume of stars to date.
Aims.
We explore such actions with the aim of locating structures in the Galactic disc.
Methods.
We computed the actions and the orbital parameters of the
Gaia
DR3 stars, assuming an axisymmetric model for the Milky Way. Using
Gaia
DR3 photometric data, we also selected a subset of giant stars with better astrometry as a control sample.
Results.
We find that the maps of the percentiles of the radial action
J
R
reveal arc-like segments. We found a high
J
R
region centered at
R
≈ 10.5 kpc of 1 kpc width, as well as three arc-shape regions dominated by circular orbits at inner radii. We also identified the spiral arms in the overdensities of the giant population.
Conclusions.
For Galactic coordinates (
X
,
Y
,
Z
), we find good agreement with the literature in the innermost region for the Scutum-Sagittarius spiral arms. At larger radii, the low
J
R
structure tracks the Local arm at negative
X
, while for the Perseus arm, the agreement is restricted to the
X
< 2 kpc region, with a displacement with respect to the literature at more negative longitudes. We detected a high
J
R
area at a Galactocentric radii of ∼10.5 kpc, consistent with some estimations of the Outer Lindblad Resonance location. We conclude that the pattern in the dynamics of the old stars is consistent in several places with the spatial distribution of the spiral arms traced by young populations, with small potential contributions from the moving groups.
Context. The Gaia astrometric sample 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. Aims. 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. Methods. 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. Results. We reveal a sharp yet complex anticentre dominated by Monoceros (MNC) and the Anticentre Stream (ACS) in the north – which we find have intensities 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 the ACS, both of which span more than 100° in longitude, and to provide a high purity sample of giants with which we track MNC down to latitudes as low as ∼5°. Their colour-magnitude diagram is consistent with extended structures at a distance of ∼10−11 kpc that originated in the disc, with a very low ratio of RR Lyrae over M giants, and with kinematics compatible with the rotation curve at those distances or slightly slower. Conclusions. We present a precise characterisation of MNC and the 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 will allow future studies to model their chemo-dynamics and evolution, thus constraining some of the most influential processes that shaped the MW.
Context. The young local associations (YLAs) constitute an excellent sample for the study of a variety of astrophysical topics, especially the star formation process in low-density environments. Data ...from the Gaia mission allows us to undertake studies of the YLAs with unprecedented accuracy. Aims. We determine the dynamical age and place of birth of a set of associations in a uniform and dynamically consistent manner. There are nine YLAs in our sample ϵ Chamaeleontis, TW Hydrae, β Pictoris, Octans, Tucana-Horologium, Columba, Carina, Argus, and AB Doradus. Methods. We designed a method for deriving the dynamical age of the YLAs based on the orbital integration. The method involves a strategy to account for the effect of observational errors. We tested the method using mock YLAs. Finally, we applied it to our set of nine YLAs with astrometry from the first Gaia data release and complementary on-ground radial velocities from the literature. Results. Our orbital analysis yields a first estimate of the dynamical age of 3 3−0+9 $^{+9}_{-0}$−0+9 Myr, 13 13−0+7 $^{+7}_{-0}$−0+7 Myr, and 5 5−0+23 $^{+23}_{-0}$−0+23 Myr for ϵ Chamaeleontis, β Pictoris, and Tucana-Horologium, respectively. For four other associations (Octans, Columba, Carina, and Argus), we provide a lower limit for the dynamical age. Our rigorous error treatment indicates that TW Hydrae and AB Doradus deserve further study. Conclusions. The dynamical ages that we obtain are compatible spectroscopic and isochrone fitting ages obtained elsewhere. From the orbital analysis, we suggest a scenario for these YLAs where there were two episodes of star formation: one ~40 Myr ago in the first quadrant that gave birth to ϵ Chamaeleontis, TW Hydrae, and β Pictoris, and another 5−15 Myr ago close to the Sun that formed Tucana-Horologium, Columba, and Carina. Future Gaia data will provide the necessary accuracy to improve the present results, especially for the controversial age determinations, and additional evidence for the proposed scenario once a complete census of YLAs and better membership can be obtained.
Extensive surveys of star-forming regions with Spitzer have revealed populations of disk-bearing young stellar objects. These have provided crucial constraints, such as the timescale of dispersal of ...protoplanetary disks, obtained by carefully combining infrared data with spectroscopic or X-ray data. While observations in various regions agree with the general trend of decreasing disk fraction with age, the Lupus V and VI regions appeared to have been at odds, having an extremely low disk fraction. Here we show, using the recent Gaia data release 2 (DR2), that these extremely low disk fractions are actually due to a very high contamination by background giants. Out of the 83 candidate young stellar objects (YSOs) in these clouds observed by Gaia, only five have distances of ~150 pc, similar to YSOs in the other Lupus clouds, and have similar proper motions to other members in this star-forming complex. Of these five targets, four have optically thick (Class II) disks. On the one hand, this result resolves the conundrum of the puzzling low disk fraction in these clouds, while, on the other hand, it further clarifies the need to confirm the Spitzer selected diskless population with other tracers, especially in regions at low galactic latitude like Lupus V and VI. The use of Gaia astrometry is now an independent and reliable way to further assess the membership of candidate YSOs in these, and potentially other, star-forming regions.
Context. Recent studies have suggested that moving groups have a dynamic or "resonant" origin. Under this hypothesis, these kinematic structures become a powerful tool for studying the large-scale ...structure and dynamics of the Milky Way. Aims. Here we aim to characterize these structures in the U-V-age-Fe/H space and establish observational constraints that will allow us to study their origin and evolution. Methods. We apply multiscale techniques-wavelet denoising (WD)-to an extensive compendium of more than 24 000 stars in the solar neighbourhood with the best available astrometric, photometric and spectroscopic data. Results. We confirm that the dominant structures in the U-V plane are the branches of Sirius, Coma Berenices, Hyades-Pleiades and Hercules. These branches are nearly equidistant in this kinematic plane and they show a negative slope. The abrupt drops in the velocity density distribution are characterized. We find a certain dependence of these kinematic structures on Galactic position with a significant change of contrast among substructures inside the branches. A large spread of ages is observed for all branches. The Hercules branch is detected in all subsamples with ages older than similar to 2 {\rm\,Gyr} and the set of the other three branches is well established for stars >400 {\rm\,Myr}. The age-metallicity relation of each branch is examined and the relation between kinematics and metallicity is studied. Conclusions. Not all of these observational constraints are successfully explained by the recent models proposed for the formation of such kinematic structures. Simulations incorporating stellar ages and metallicities are essential for future studies. The comparison of the observed and simulated distributions obtained by WD will provide a physical interpretation of the existence of the branches in terms of local or large-scale dynamics.
In this first paper, we simulate the population of disc red clump stars to be observed by Gaia. We generate a set of test particles and we evolve it in a 3D barred Milky Way like galactic potential. ...We assign physical properties of the red clump trace population and a realistic 3D interstellar extinction model. We add Gaia observational constraints and an error model according to the pre-commissioning scientific performance assessments. We present and analyse two mock catalogues, offered to the community, that are an excellent test bed for testing tools being developed for the future scientific exploitation of Gaia data. The first catalogue contains stars up to Gaia
G∼20, while the second is the subset containing Gaia radial velocity data with a maximum error of
$\sigma _{V_{\rm r}}=10$
km s−1. Here, we present first attempts to characterize the density structure of the Galactic bar in the Gaia space of observables. The Gaia large errors in parallax and the high interstellar extinction in the inner parts of the Galactic disc prevent us to model the bar overdensity. This result suggests the need to combine Gaia and IR data to undertake such studies. We find that IR photometric distances for this Gaia sample allow us to recover the Galactic bar orientation angle with an accuracy of ∼5°.
We use self-consistent numerical simulations of the evolution and disruption of the Sun's birth cluster in the Milky Way potential to investigate the present-day phase-space distribution of the Sun's ...siblings. The simulations include the gravitational N-body forces within the cluster and the effects of stellar evolution on the cluster population. In addition, the gravitational forces due to the Milky Way potential are accounted for in a self-consistent manner. Our aim is to understand how the astrometric and radial velocity data from the Gaia mission can be used to pre-select solar sibling candidates. We vary the initial conditions of the Sun's birth cluster, as well as the parameters of the Galactic potential. In particular, we use different configurations and strengths of the bar and spiral arms. We show that the disruption time-scales of the cluster are insensitive to the details of the non-axisymmetric components of the Milky Way model and we make predictions, averaged over the different simulated possibilities, about the number of solar siblings that should appear in surveys such as Gaia or GALAH. We find a large variety of present-day phase-space distributions of solar siblings, which depend on the cluster initial conditions and the Milky Way model parameters. We show that nevertheless robust predictions can be made about the location of the solar siblings in the space of parallaxes (ϖ), proper motions (μ) and radial velocities (V
r). By calculating the ratio of the number of simulated solar siblings to that of the number of stars in a model Galactic disc, we find that this ratio is above 0.5 in the region given by: ϖ ≥ 5 mas, 4 ≤ μ ≤ 6 mas yr−1, and −2 ≤ V
r ≤ 0 km s−1. Selecting stars from this region should increase the probability of success in identifying solar siblings through follow-up observations. However the proposed pre-selection criterion is sensitive to our assumptions, in particular about the Galactic potential. Using a more realistic potential (e.g. including transient spiral structure and molecular clouds) would make the pre-selection of solar sibling candidates based on astrometric and radial velocity data very inefficient. This reinforces the need for large-scale surveys to determine precise astrophysical properties of stars, in particular their ages and chemical abundances, if we want to identify the solar family.