Test particle simulations of Milky Way type galactic disks are being conducted to map the evolution of the stellar large scale kinematic response to the bar and spiral structure. Second and third order moments of the velocity distribution function prove to be good indicators of both, the velocity ellipsoid misalignment near the arms and the bar, and the degree of kinematic substructure in the UV plane, that is, the presence of moving groups. A large scale analysis all through the galactic disk allows to establish the kinematic behavior near resonances and the correlation between the kinematic parameters and properties of the non-axisymmetric components, such as its transient nature or its mass overdensity. N-body simulaions are being run in order to test these results in self-consistent models.
We study the components of cool and warm/hot gas in the circumgalactic medium (CGM) of simulated galaxies and address the relative production of OVI by photoionization versus collisional ionization, as a function of halo mass, redshift, and distance from the galaxy halo center. This is done utilizing two different suites of zoom-in hydro-cosmological simulations, VELA (6 halos; $z>1$) and NIHAO (18 halos; to $z=0$), which provide a broad theoretical basis because they use different codes and physical recipes for star formation and feedback. In all halos studied in this work, we find that collisional ionization by thermal electrons dominates at high redshift, while photoionization of cool or warm gas by the metagalactic radiation takes over near $z\sim2$. In halos of $\sim 10^{12}M_{\odot}$ and above, collisions become important again at $z 3\times10^{11}~M_{\odot}$, at $z\sim 0$ most of the photoionized OVI is in a warm, not cool, gas phase ($T\lesssim 3\times 10^5$~K). We also find that collisions are dominant in the central regions of halos, while photoionization is more significant at the outskirts, around $R_{\textrm v}$, even in massive halos. This too may be explained by the presence of warm gas or, in lower mass halos, by cool gas inflows.
Stellar tidal streams are a key tracer of galaxy evolution and have the potential to provide an indirect means for tracing dark matter. For the Local Group, many diffuse substructures have been identified and their link to galaxy evolution has been traced. However, the Local Group does not offer a statistically significant sample of stellar tidal streams. Thus, an analysis of a larger sample beyond the Local Group is required to better probe the frequency and characteristics of these streams to verify whether these properties are in agreement with the predictions of the Lambda CDM model and its implementation in cosmological simulations, taking into account the impact of the baryonic physics modelling. The main scope of the Stellar Stream Legacy Survey is to obtain a statistically significant sample of stellar streams in the local Universe to be able to trace and study minor mergers and their contribution to galaxy evolution with respect to the Lambda CDM theory. For that purpose, we are carrying out the first systematic survey of faint stellar debris from tidally disrupted dwarf satellites around nearby galaxies up to a distance of 100 Mpc. In this paper, we present a catalogue with the results of the first harvest of stellar tidal streams found by visual inspection in deep images of sim 700 galaxies from the Dark Energy Survey (DES). We also include, for the first time, a photometric characterisation of the streams obtained by measuring their surface brightnesses and colours. We found a total of 63 streams in our sample at distances between 40 and 100 Mpc, including 58 that were not previously reported. We measured their average surface brightness for the $g$ band, the $r$ band and the $z$ band, to be 28.35pm 0.20, 27.81pm 0.13, and 27.62pm 0.09 mag arcsec$^ $, respectively. By applying a statistical analysis to our findings, we obtained a stream detection frequency of 9.1<!PCT!> pm 1.1<!PCT!> for the given surface brightness limit of the DES image sample, in agreement with previous studies. We identified stream progenitors in 5--14<!PCT!> of our stream sample, depending on the confidence level. The first catalogue of streams in the local Universe presented here will be complemented by future stream surveys within the Stellar Stream Legacy Survey and can be exploited in studies pertaining to galaxy evolution and cosmological models. In this work, we have learnt that the faintest measured stream surface brightness can be significantly brighter than the surface brightness limit of an image measured at the pixel level (in our case up to sim 1 mag arcsec$^ $ for the $r$ band) mainly due to correlated noise present in the images.
Context: The Large Magellanic Cloud (LMC) internal kinematics have been studied in unprecedented depth thanks to the excellent quality of the Gaia mission data, revealing the disc's non-axisymmetric structure. Aims: We want to constrain the LMC bar pattern speed using the astrometric and spectroscopic data from the Gaia mission. Methods: We apply three methods to evaluate the bar pattern speed: it is measured through the Tremaine-Weinberg (TW) method, the Dehnen method and a bisymmetric velocity (BV) model. The methods provide additional information on the bar properties such as the corotation radius and the bar length and strength. The validity of the methods is tested with numerical simulations. Results: A wide range of pattern speeds are inferred by the TW method, owing to a strong dependency on the orientation of the galaxy frame and the viewing angle of the bar perturbation. The simulated bar pattern speeds (corotation radii, respectively) are well recovered by the Dehnen method (BV model). Applied to the LMC data, the Dehnen method finds a pattern speed Omega_p = -1.0 +/- 0.5 km s-1 kpc-1, thus corresponding to a bar which barely rotates, slightly counter-rotating with respect to the LMC disc. The BV method finds a LMC bar corotation radius of Rc = 4.20 +/- 0.25 kpc, corresponding to a pattern speed Omega_p = 18.5^{+1.2}_{-1.1} km s-1 kpc-1. Conclusions: It is not possible to decide which global value best represents an LMC bar pattern speed with the TW method, due to the strong variation with the orientation of the reference frame. The non-rotating bar from the Dehnen method would be at odds with the structure and kinematics of the LMC disc. The BV method result is consistent with previous estimates and gives a bar corotation-to-length ratio of 1.8 +/- 0.1, which makes the LMC hosting a slow bar.
ABSTRACT We introduce a new set of simulations of Milky Way (MW)-sized galaxies using the AMR code ART + hydrodynamics in a Λ cold dark matter cosmogony. The simulation series is called GARROTXA and it follows the formation of a halo/galaxy from z = 60 to z = 0. The final virial mass of the system is ∼7.4 × 10 11 M ⊙ . Our results are as follows. (a) Contrary to many previous studies, the circular velocity curve shows no central peak and overall agrees with recent MW observations. (b) Other quantities, such as (6 × 10 10 M ⊙ ) and R d (2.56 kpc), fall well inside the observational MW range. (c) We measure the disk-to-total ratio kinematically and find that D/T = 0.42. (d) The cold-gas fraction and star formation rate at z = 0, on the other hand, fall short of the values estimated for the MW. As a first scientific exploitation of the simulation series, we study the spatial distribution of hot X-ray luminous gas. We have found that most of this X-ray emitting gas is in a halo-like distribution accounting for an important fraction but not all of the missing baryons. An important amount of hot gas is also present in filaments. In all our models there is not a massive disk-like hot-gas distribution dominating the column density. Our analysis of hot-gas mock observations reveals that the homogeneity assumption leads to an overestimation of the total mass by factors of 3–5 or to an underestimation by factors of 0.7–0.1, depending on the used observational method. Finally, we confirm a clear correlation between the total hot-gas mass and the dark matter halo mass of galactic systems.
Only recently, complex models that include the global dynamics from dwarf satellite galaxies, dark matter halo structure, gas infalls, and stellar disk in a cosmological context became available to study the dynamics of disk galaxies such as the Milky Way (MW). We use a MW model from a high-resolution hydrodynamical cosmological simulation named GARROTXA to establish the relationship between the vertical disturbances seen in its galactic disk and multiple perturbations, from the dark matter halo, satellites and gas. We calculate the bending modes in the galactic disk in the last 6 Gyr of evolution. To quantify the impact of dark matter and gas we compute the vertical acceleration exerted by these components onto the disk and compare them with the bending behavior with Fourier analysis. We find complex bending patterns at different radii and times, such as an inner retrograde mode with high frequency, as well as an outer slower retrograde mode excited at different times. The amplitudes of these bending modes are highest during the early stages of the thin disk formation and reach up to 8.5 km s-1 in the late disk evolution. We find that the infall of satellite galaxies leads to a tilt of the disk, and produces anisotropic gas accretion with subsequent star formation events, and supernovae, creating significant vertical accelerations onto the disk plane. The misalignment between the disk and the inner stellar/dark matter triaxial structure, formed during the ancient assembly of the galaxy, creates a strong vertical acceleration on the stars. We conclude that several agents trigger the bending of the stellar disk and its phase spirals in this simulation, including satellite galaxies, dark sub-halos, misaligned gaseous structures, and the inner dark matter profile, which coexist and influence each other, making it challenging to establish direct causality.
In this fourth paper from the AGORA Collaboration, we study the evolution down to redshift $z=2$ and below of a set of cosmological zoom-in simulations of a Milky Way mass galaxy by eight of the leading hydrodynamic simulation codes. We also compare this CosmoRun suite of simulations with dark matter-only simulations by the same eight codes. We analyze general properties of the halo and galaxy at $z=4$ and 3, and before the last major merger, focusing on the formation of well-defined rotationally-supported disks, the mass-metallicity relation, the specific star formation rate, the gas metallicity gradients, and the non-axisymmetric structures in the stellar disks. Codes generally converge well to the stellar-to-halo mass ratios predicted by semi-analytic models at $z\sim$2. We see that almost all the hydro codes develop rotationally-supported structures at low redshifts. Most agree within 0.5 dex with the observed MZR at high and intermediate redshifts, and reproduce the gas metallicity gradients obtained from analytical models and low-redshift observations. We confirm that the inter-code differences in the halo assembly history reported in the first paper of the collaboration also exist in CosmoRun, making the code-to-code comparison more difficult. We show that such differences are mainly due to variations in code-dependent parameters that control the time-stepping strategy of the gravity solver. We find that variations in the early stellar feedback can also result in differences in the timing of the low-redshift mergers. All the simulation data down to $z=2$ and the auxiliary data will be made publicly available.
The nature of the few known solar-mass stars simultaneously containing debris disks and planets remains an open question. A number of works have shown that this property appears to be independent of planetary masses as well as of stellar age, but possible correlations with stellar kinematics and metallicity have not been investigated. In this paper, we show that the majority of known stars containing both debris disks and planets belong to the metal-enriched Galactic thin disk. The few exceptions are stars that seem to be born in the star formation peak occurring in times of thick disk formation (i.e., HD 10700, HD 20794, and HD 40307), that is, between 11 and 8 Gyr. The mass of the dusty disk of these three old stars measured at 70 $\mu$m is very small - in fact, it is lower than that of the Kuiper belt of our Solar system by several orders of magnitude. These results are not surprising, as they remain within the values expected for the stellar disk evolution of such primitive stars. In parallel, we found another six thick-disk stars containing only debris disks or planets. These results enable us to establish a correlation between stellar metallicity and the mass of the dust disk modulated by the different formation epochs of the thick and thin Galactic disks.
We model the effects of the spiral arms of the Milky Way on the disk stellar kinematics in the Gaia observable space. We also estimate the Gaia capabilities of detecting the predicted signatures. We use both controlled orbital integrations in analytic potentials and self-consistent simulations. We introduce a new strategy, which consists of comparing the stellar kinematics of symmetric Galactic longitudes (+l and -l), in particular the median transverse velocity (from parallaxes and proper motions). This approach does not require the assumption of an axisymmetric model. The typical differences between the transverse velocity in symmetric longitudes in the models are of ~2 km/s, but can be larger than 10 km/s for certain longitudes and distances. The kinematic differences for +l and -l show trends that depend on the properties of spiral arms. Thus, this method can be used to quantify the importance of the effects of spiral arms on the orbits of stars in the different disk regions, and to constrain the location of the arms, main resonances and, thus, pattern speed. Moreover, the method allows us to test the dynamical nature of the spiral structure (e.g. grand design versus transient multiple arms). We estimate the number of stars of each spectral type that Gaia will observe in representative Galactic longitudes, their errors in distance and transverse velocity, and the error in computing the median velocity as a function of distance. We will be able to measure the median transverse velocity with Gaia data, with precision smaller than ~1 km/s up to distances of ~4-6 kpc for certain giant stars, and up to ~2-4 kpc and better kinematic precision (<= 0.5 km\s) for certain sub-giants and dwarfs. These are enough to measure the typical signatures seen in the models. We also show that a similar strategy can be used with line-of-sight velocities from other upcoming spectroscopic surveys.(Abridged)
Context. Complex models recently became available for studying the dynamics of disk galaxies such as the Milky Way (MW). These models include the global dynamics from dwarf satellite galaxies, dark matter halo structure, gas infall, and stellar disks in a cosmological context. Aims. We use a MW model from a suite of high-resolution hydrodynamical cosmological simulations named GARROTXA to establish the relationship between the vertical disturbances seen in its galactic disk and multiple perturbations from the dark matter halo, satellites, and gas. Methods. We calculated the bending modes in the galactic disk in the last 6 Gyr of evolution. We computed the vertical acceleration exerted by dark matter and gas in order to quantify the impact of these components on the disk, and compared this with the bending behavior with Fourier analysis. Results. We find complex bending patterns at different radii and times, such as an inner retrograde mode with high frequency and an outer slower retrograde mode excited at different times. The amplitudes of these bending modes are highest during the early stages of formation of the thin disk (20 km s −1 ) and reach up to 8.5 km s −1 in the late disk evolution. We find that the infall of satellite galaxies leads to a tilt of the disk, and produces strong anisotropic gas accretion with a misalignment of 8° with subsequent star formation events and supernovae, creating significant vertical accelerations on the disk plane. The misalignment between the disk and the inner stellar and dark matter triaxial structure, which formed during the ancient assembly of the galaxy, also leads to a strong vertical acceleration of the stars. We also find dark matter subhalos that temporally coincide with the appearance of bending waves in certain periods. Conclusions. We conclude that several agents trigger the bending of the stellar disk and its phase spirals in this simulation, including satellite galaxies, dark subhalos, misaligned gaseous structures, and the inner dark matter profile. These phenomena coexist and influence each other, sometimes making it challenging to establish direct causality.