No Arabic abstract
The present paper is the culminating one of a series aimed to contribute to the understanding of the kinematic structures of the solar neighbourhood (SN), explaining the origin of the Local Arm and relating the moving groups with the spiral-arms resonances in the disk. With a model for the Galactic potential, with the Sun inside the spiral corotation resonance (CR), we integrate the 2D orbits of test particles distributed in birthplaces along the main spiral arms, the Local Arm, and in the axisymmetric disk. A comparison of the resulting U-V plane of the SN with that provided by Gaia DR2 confirms our previous conclusion that the moving groups of Coma Berenices, Pleiades, and Hyades are associated with the CR, and that the Hercules stream is formed by the bulk of high-order inner Lindblad resonances. The kinematic structures result from stellar orbits trapped by the spiral resonances in a timespan of ~ 1 Gyr, indicating the long-living nature of the spiral structure and challenging recent arguments in favor of short-lived structures originated from incomplete phase mixing in the Galactic disk. As a byproduct, our simulations give some insight into the birthplaces of the stars presently located in the SN; the majority of stars of the main moving groups and the Hercules stream were likely born in the Local Arm, while stars of the Sirius group possibly originated from the outer segment of the Sagittarius-Carina arm. We also propose the spiral resonances as the dynamical origin for the diagonal ridges in the Galactic distribution of rotation velocities.
Ever since a thick disk was proposed to explain the vertical distribution of the Milky Way disk stars, its origin has been a recurrent question. We aim to answer this question by inspecting 19 disk galaxies with stellar mass greater than $10^{10},rm M_odot$ in recent cosmological high-resolution zoom-in simulations: Galactica and NewHorizon. The thin and thick disks are reasonably reproduced by the simulations with scale heights and luminosity ratios as observed. We then spatially classify the thin and thick disks and find that the thick disk stars are older, metal-poorer, kinematically-hotter, and higher in accreted star fraction, while both disks are dominated by the stars formed in situ. Half of the in-situ stars in the thick disks are formed before the galaxies develop their disks, and the rest are formed in spatially and kinematically thinner disks and then thickened with time by heating. However, the 19 galaxies have various properties and evolutionary routes, highlighting the need for statistically-large samples to draw general conclusions. We conclude from our simulations that the thin and thick disk components are not entirely distinct in terms of formation processes, but rather markers of the evolution of galactic disks. Moreover, as the combined result of the thickening of the existing disk stars and the continued formation of young thin-disk stars, the vertical distribution of stars does not change much after the disks settle, pointing to the modulation of both orbital diffusion and star formation by the same confounding factor: the proximity of galaxies to marginal stability.
Using the Gaia data release 2 (DR2), we analyzed the distribution of stars in the close vicinity of the Sun in the full 3D position-velocity space. We have found no evidence of incomplete phase mixing in the vertical direction of the disk, which could be originated by some external events. We show that the vertical phase space spiral $Z$-$V_z$ is produced by the well-known moving groups (MGs), mainly by Coma-Berenices, Pleiades-Hyades and Sirius, when the statistical characteristics (mean, median, or mode) of the azimuthal velocity $V_varphi$ are used to analyze the distribution in the vertical position-velocity plane. This result does not invoke external perturbations and is independent on the internal dynamical mechanisms that originate the MGs. Our conclusions counterbalance current arguments in favor of short-lived (between 300 and 900 Myr) structures in the solar neighborhood. Contrarily, they support the hypothesis of a longer formation time scale (around a few Gyr) for the MGs.
We investigate the properties and evolution of star particles in two simulations of isolated spiral galaxies, and two galaxies from cosmological simulations. Unlike previous numerical work, where typically each star particle represents one `cluster, for the isolated galaxies we are able to model features we term `clusters with groups of particles. We compute the spatial distribution of stars with different ages, and cluster mass distributions, comparing our findings with observations including the recent LEGUS survey. We find that spiral structure tends to be present in older (100s Myrs) stars and clusters in the simulations compared to the observations. This likely reflects differences in the numbers of stars or clusters, the strength of spiral arms, and whether the clusters are allowed to evolve. Where we model clusters with multiple particles, we are able to study their evolution. The evolution of simulated clusters tends to follow that of their natal gas clouds. Massive, dense, long-lived clouds host massive clusters, whilst short-lived clouds host smaller clusters which readily disperse. Most clusters appear to disperse fairly quickly, in basic agreement with observational findings. We note that embedded clusters may be less inclined to disperse in simulations in a galactic environment with continuous accretion of gas onto the clouds than isolated clouds and correspondingly, massive young clusters which are no longer associated with gas tend not to occur in the simulations. Caveats of our models include that the cluster densities are lower than realistic clusters, and the simplistic implementation of stellar feedback.
Tidal disruption events (TDEs) probe properties of supermassive black holes (SMBHs), their accretion disks, and the surrounding nuclear stellar cluster. Light curves of TDEs are related to orbital properties of stars falling SMBHs. We study the origin, density, and velocity distributions of bound and unbound stars in the nuclear star cluster, which are causing TDEs as a function of their orbital eccentricity $e$ and energy $E$. These quantities determine near the SMBH the ratio of the orbits pericenter to tidal disruption radii (denoted as penetration factor, $beta$). We develop an analytical model for the density and velocity distribution of such stars in the cluster, which agrees well with N-body experiments. Our model extends classical models of angular momentum diffusion in the loss cone. We also derive an analytical model for three characteristic eccentricities in the loss cone: the minimum and maximum value for given $beta$, respectively, and $e_{rm lcb}$, which represents the orbital eccentricity defining the boundary between empty and full loss cone regimes. With N-body experiments, we show that stars causing TDEs are distributed between these eccentricity limits on the $e-beta$ plane. Moreover, we find most of the bound stars between $e_{rm lcb}$ and $e=1$ (i.e., the full loss cone regime), whereas the remaining bound stars are originating from the empty loss cone regime. This is consistent with the loss cone theory. We propose that the $e-beta$ distribution of stars in a star cluster or galactic nucleus can be a good tool to diagnose whether the stars can cause TDEs.
Using the astrometry from the ESAs Gaia mission, previous works have shown that the Milky Way stellar halo is dominated by metal-rich stars on highly eccentric orbits. To shed light on the nature of this prominent halo component, we have analysed 28 Galaxy analogues in the Auriga suite of cosmological hydrodynamics zoom-in simulations. Some three quarters of the Auriga galaxies contain significant components with high radial velocity anisotropy, beta > 0.6. However, only in one third of the hosts do the high-beta stars contribute significantly to the accreted stellar halo overall, similar to what is observed in the Milky Way. For this particular subset we reveal the origin of the dominant stellar halo component with high metallicity, [Fe/H]~-1, and high orbital anisotropy, beta>0.8, by tracing their stars back to the epoch of accretion. It appears that, typically, these stars come from a single dwarf galaxy with a stellar mass of order of 10^9-10^10 Msol that merged around 6-10 Gyr ago, causing a sharp increase in the halo mass. Our study therefore establishes a firm link between the excess of radially anisotropic stellar debris in the Milky Way halo and an ancient head-on collision between the young Milky Way and a massive dwarf galaxy