ABSTRACT We explore the chemodynamical properties of a sample of barred galaxies in the Auriga magnetohydrodynamical cosmological zoom-in simulations, which form boxy/peanut (b/p) bulges, and compare these to the Milky Way (MW). We show that the Auriga galaxies which best reproduce the chemodynamical properties of stellar populations in the MW bulge have quiescent merger histories since redshift z ∼ 3.5: their last major merger occurs at $t_{\rm lookback}\gt 12\, \rm Gyr$, while subsequent mergers have a stellar mass ratio of ≤1:20, suggesting an upper limit of a few per cent for the mass ratio of the recently proposed Gaia Sausage/Enceladus merger. These Auriga MW-analogues have a negligible fraction of ex-situ stars in the b/p region ($\lt 1{{\ \rm per\ cent}}$), with flattened, thick disc-like metal-poor stellar populations. The average fraction of ex-situ stars in the central regions of all Auriga galaxies with b/p’s is 3 per cent – significantly lower than in those which do not host a b/p or a bar. While the central regions of these barred galaxies contain the oldest populations, they also have stars younger than 5 Gyr (>30 per cent) and exhibit X-shaped age and abundance distributions. Examining the discs in our sample, we find that in some cases a star-forming ring forms around the bar, which alters the metallicity of the inner regions of the galaxy. Further out in the disc, bar-induced resonances lead to metal-rich ridges in the Vϕ − r plane – the longest of which is due to the Outer Lindblad Resonance. Our results suggest the Milky Way has an uncommonly quiet merger history, which leads to an essentially in-situ bulge, and highlight the significant effects the bar can have on the surrounding disc.
ABSTRACT Merging of galaxy clusters are some of the most energetic events in the Universe, and they provide a unique environment to study galaxy evolution. We use a sample of 84 merging and relaxed SPT galaxy clusters candidates, observed with the Dark Energy Camera in the 0.11 < z < 0.88 redshift range, to build colour–magnitude diagrams to characterize the impact of cluster mergers on the galaxy population. We divided the sample between relaxed and disturbed, and in two redshifts bin at z = 0.55. When comparing the high-z to low-z clusters we find the high-z sample is richer in blue galaxies, independently of the cluster dynamical state. In the high-z bin, we find that disturbed clusters exhibit a larger scatter in the red sequence, with wider distribution and an excess of bluer galaxies compared to relaxed clusters, while in the low-z bin we find a complete agreement between the relaxed and disturbed clusters. Our results support the scenario in which massive cluster halos at z < 0.55 galaxies are quenched as satellites of another structure, i.e. outside the cluster, while at z ≥ 0.55 the quenching is dominated by in situ processes.
Galactic bars drive the internal evolution of spiral galaxies, while their formation is tightly coupled to the properties of their host galaxy and dark matter halo. To explore what drives bar formation in the cosmological context and how these structures evolve throughout cosmic history, we use the Auriga suite of magneto-hydrodynamical cosmological zoom-in simulations. We find that bars are robust and long-lived structures, and we recover a decreasing bar fraction with increasing redshift which plateaus around $\sim20\%$ at $z\sim3$. We find that bars which form at low and intermediate redshifts grow longer with time, while bars that form at high redshifts are born `saturated' in length, likely due to their merger-induced formation pathway. This leads to a larger bar-to-disc size ratio at high redshifts as compared to the local Universe. We subsequently examine the multi-dimensional parameter space thought to drive bar formation. We find that barred galaxies tend to have lower Toomre $Q$ values at the time of their formation, while we do not find a difference in the gas fraction of barred and unbarred populations when controlling for stellar mass. Barred galaxies tend to be more baryon-dominated at all redshifts, assembling their stellar mass earlier, while galaxies that are baryon-dominated but that do not host a bar, have a higher ex-situ bulge fraction. We explore the implications of the baryon-dominance of barred galaxies on the Tully-Fisher relation, finding an offset from the unbarred relation; confirming this in observations would serve as additional evidence for dark matter, as this behaviour is not readily explained in modified gravity scenarios.
We study the phase-space behaviour of nearby trajectories in integrable potentials. We show that the separation of nearby orbits initially diverges very fast, mimicking a nearly exponential behaviour, while at late times it grows linearly. This initial exponential phase, known as Miller's instability, is commonly found in N-body simulations, and has been attributed to short-term (microscopic) N-body chaos. However we show here analytically that the initial divergence is simply due to the shape of an orbit in phase-space. This result confirms previous suspicions that this transient phenomenon is not related to an instability in the sense of non-integrable behaviour in the dynamics of N-body systems.
The Milky Way has accreted many ultra-faint dwarf galaxies (UFDs), and stars from these galaxies can be found throughout our Galaxy today. Studying these stars provides insight into galaxy formation and early chemical enrichment, but identifying them is difficult. Clustering stellar dynamics in 4D phase space ($E$, $L_z$, $J_r$, $J_z$) is one method of identifying accreted structure which is currently being utilized in the search for accreted UFDs. We produce 32 simulated stellar halos using particle tagging with the \textit{Caterpillar} simulation suite and thoroughly test the abilities of different clustering algorithms to recover tidally disrupted UFD remnants. We perform over 10,000 clustering runs, testing seven clustering algorithms, roughly twenty hyperparameter choices per algorithm, and six different types of data sets each with up to 32 simulated samples. Of the seven algorithms, HDBSCAN most consistently balances UFD recovery rates and cluster realness rates. We find that even in highly idealized cases, the vast majority of clusters found by clustering algorithms do not correspond to real accreted UFD remnants and we can generally only recover $6\%$ of UFDs remnants at best. These results focus exclusively on groups of stars from UFDs, which have weak dynamic signatures compared to the background of other stars. The recoverable UFD remnants are those that accreted recently, $z_{\text{accretion}}\lesssim 0.5$. Based on these results, we make recommendations to help guide the search for dynamically-linked clusters of UFD stars in observational data. We find that real clusters generally have higher median energy and $J_r$, providing a way to help identify real vs. fake clusters. We also recommend incorporating chemical tagging as a way to improve clustering results.
We analyse the orbital kinematics of the Milky Way (MW) satellite system utilizing the latest systemic proper motions for 38 satellites based on data from Gaia Data Release 2. Combining these data with distance and line-of-sight velocity measurements from the literature, we use a likelihood method to model the velocity anisotropy, β, as a function of Galactocentric distance and compare the MW satellite system with those of simulated MW-mass haloes from the APOSTLE (A Project Of Simulating The Local Environment) and Auriga simulation suites. The anisotropy profile for the MW satellite system increases from β ∼ −2 at r ∼ 20 kpc to β ∼ 0.5 at r ∼ 200 kpc, indicating that satellites closer to the Galactic centre have tangentially biased motions while those farther out have radially biased motions. The motions of satellites around APOSTLE host galaxies are nearly isotropic at all radii, while the β(r) profiles for satellite systems in the Auriga suite, whose host galaxies are substantially more massive in baryons than those in APOSTLE, are more consistent with that of the MW satellite system. This shape of the β(r) profile may be attributed to the central stellar disc preferentially destroying satellites on radial orbits, or intrinsic processes from the formation of the MW system.
We use the Auriga cosmological simulations of Milky Way (MW)-mass galaxies and their surroundings to study the satellite populations of dwarf galaxies in $\Lambda$CDM. As expected from prior work, the number of satellites above a fixed stellar mass is a strong function of the mass of the primary dwarf. For galaxies as luminous as the Large Magellanic Cloud (LMC), and for halos as massive as expected for the LMC (determined by its rotation speed), the simulations predict about 3 satellites with stellar masses exceeding $M_*>10^5\, M_\odot$. If the LMC is on its first pericentric passage, then these satellites should be near the LMC and should have orbital angular momenta roughly coincident with that of the LMC. We use 3D positions and velocities from the 2nd data release of the Gaia mission to revisit which of the "classical" MW dwarf spheroidals could plausibly be LMC satellites. The new proper motions of the Fornax and Carina dwarf spheroidals place them on orbits closely aligned with the orbital plane of the Magellanic Clouds, hinting at a potential Magellanic association. Together with the Small Magellanic Cloud (SMC), this result raises to $3$ the number of LMC satellites with $M_*>10^5\, M_\odot$, as expected from simulations. This also fills the 12-mag luminosity gap between the SMC and the ultra-faints Hyi1, Car2, Hor1, and Car3, the few ultra-faint satellites confirmed to have orbits consistent with a Magellanic origin.
Current studies of large-scale asymmetries (i.e. lopsidedness) in the stellar density distribution of disk galaxies have mainly focused on the local Universe. Recent observations have found a significant fraction (over 60%) of lopsided galaxies at high-redshift ($1.5 < z < 3$), which is significantly larger than the fraction (~30%) observed in the nearby Universe. We aim to understand whether the more widespread lopsidedness at high- than low-redshift can be associated to environmental mechanisms being more effective in producing lopsided perturbations at high-redshift. At each redshift between $0 < z < 2$, we independently select a sample of disk-like galaxies from the IllustrisTNG simulations. We then characterize lopsidedness in the disks of galaxies at each redshift, study the relevant mechanisms generating lopsidedness, as well as the correlation between such perturbation, the local environment and the galaxy internal properties as a function of redshift. Consistent with previous and new observational results, we find that: 1) simulations predict a significant fraction (~60%) of lopsided galaxies at high-redshift ($1.5 < z < 2$), 2) the fraction of lopsided galaxies, as well as the lopsided amplitude, decreases from high- to low-redshift, and 3) there is not a significant dependence of lopsidedness on the local environment, but there is a strong correlation between the lopsided amplitude and basic galaxies' structural properties at all redshift between $0 < z < 2$. This means that, independent of the mechanisms on-setting lopsidedness, galaxies with low central stellar mass density and more extended disks are more susceptible of developing strong lopsidedness. We find that both recent interactions with mass-ratio >1:10 and gas accretion with subsequent star formation can produce lopsided perturbations at all redshift, but they are both significantly more effective at high-redshift.
Lopsided galaxies are late-type galaxies with a non-axisymmetric disc due to an uneven distribution of their stellar mass. Despite being a relatively common perturbation, several questions regarding its origin and the information that can be extracted from them about the evolutionary history of late-type galaxies. The advent of several large multi-band photometric surveys will allow us to statistically analyze this perturbation, with information that was not previously available. Given the strong correlation between lopsidedness and the structural properties of the galaxies, this paper aims to develop a method to automatically classify late-type galaxies between lopsided and symmetric. We seek to explore if an accurate classification can be obtain by only considering their internal properties, without additional information about the environment. We select 8000 late-type galaxies from TNG50. A Fourier decomposition of their stellar mass surface density is used to label galaxies as lopsided and symmetric. We trained a Random Forest classifier to rapidly and automatically identify this type of perturbations, exclusively using galaxies internal properties. We test different algorithms to deal with the imbalance of our data and select the most suitable approach based on the considered metrics. We show that our trained algorithm can provide a very accurate and rapid classification of lopsided galaxies. The excellent results obtained by our classifier strongly supports the hypothesis that lopsidedness is mainly a tracer of galaxies internal structures. We show that similar results can be obtained using observable quantities, readily obtainable from multi-bad photometric surveys. Our results show it allows a rapid and accurate classification of lopsided galaxies, allowing us to explore whether lopsidedness in present-day disc galaxies is connected to galaxies specific evolutionary histories.
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.
'%3e%3cpath fill='%23012169' d='M0 0v30h60V0z'/%3e%3cpath stroke='%23fff' stroke-width='6' d='m0 0 60 30m0-30L0 30'/%3e%3cpath stroke='%23C8102E' stroke-width='4' d='m0 0 60 30m0-30L0 30' clip-path='url(%23b)'/%3e%3cpath stroke='%23fff' stroke-width='10' d='M30 0v30M0 15h60'/%3e%3cpath stroke='%23C8102E' stroke-width='6' d='M30 0v30M0 15h60'/%3e%3c/g%3e%3c/svg%3e)
'/%3e%3cuse xlink:href='%23a' transform='rotate(45 400 250)'/%3e%3cuse xlink:href='%23a' transform='rotate(67.5 400 250)'/%3e%3cpath id='b' fill='%2385340a' d='M404.31 274.41c.453 9.054 5.587 13.063 4.579 21.314 2.213-6.525-3.124-11.583-2.82-21.22m-7.649-23.757 19.487 42.577-16.329-43.887'/%3e%3cuse xlink:href='%23b' transform='rotate(22.5 400 250)'/%3e%3cuse xlink:href='%23b' transform='rotate(45 400 250)'/%3e%3cuse xlink:href='%23b' transform='rotate(67.5 400 250)'/%3e%3c/g%3e%3cuse xlink:href='%23c' transform='rotate(90 400 250)'/%3e%3cuse xlink:href='%23c' transform='rotate(180 400 250)'/%3e%3cuse xlink:href='%23c' transform='rotate(270 400 250)'/%3e%3ccircle cx='400' cy='250' r='27.778' fill='%23f6b40e' stroke='%2385340a' stroke-width='1.5'/%3e%3cpath id='h' fill='%23843511' d='M409.47 244.06c-1.897 0-3.713.822-4.781 2.531 2.136 1.923 6.856 2.132 10.062-.219a7.333 7.333 0 0 0-5.281-2.312zm-.031.438c1.846-.034 3.571.814 3.812 1.656-2.136 2.35-5.55 2.146-7.687.437.935-1.495 2.439-2.067 3.875-2.094z'/%3e%3cuse xlink:href='%23d' transform='matrix(-1 0 0 1 800.25 0)'/%3e%3cuse xlink:href='%23e' transform='matrix(-1 0 0 1 800.25 0)'/%3e%3cuse xlink:href='%23f' transform='translate(18.862)'/%3e%3cuse xlink:href='%23g' transform='matrix(-1 0 0 1 800.25 0)'/%3e%3cpath fill='%2385340a' d='M395.75 253.84c-.913.167-1.563.977-1.563 1.906 0 1.062.878 1.906 1.938 1.906a1.89 1.89 0 0 0 1.563-.812c.739.556 1.764.615 2.312.625.084.002.193 0 .25 0 .548-.01 1.573-.069 2.313-.625.36.516.935.812 1.562.812 1.06 0 1.938-.844 1.938-1.906 0-.929-.65-1.74-1.563-1.906.513.18.844.676.844 1.219a1.28 1.28 0 0 1-1.281 1.281c-.68 0-1.242-.54-1.282-1.219-.208.417-1.034 1.655-2.656 1.719-1.622-.064-2.447-1.302-2.656-1.719-.04.679-.6 1.219-1.281 1.219a1.28 1.28 0 0 1-1.281-1.281c0-.542.33-1.038.843-1.219zm2.09 5.69c-2.138 0-2.983 1.937-4.906 3.219 1.068-.427 1.91-1.27 3.406-2.125 1.496-.855 2.772.187 3.625.187h.031c.853 0 2.13-1.041 3.625-.187 1.497.856 2.369 1.698 3.438 2.125-1.924-1.282-2.8-3.219-4.938-3.219-.426 0-1.271.23-2.125.656h-.031c-.853-.426-1.698-.656-2.125-.656z'/%3e%3cpath fill='%2385340a' d='M397.12 262.06c-.844.037-1.96.207-3.563.688 3.848-.855 4.697.437 6.407.437h.03c1.71 0 2.56-1.292 6.407-.438-4.274-1.282-5.124-.437-6.406-.437h-.031c-.802 0-1.437-.312-2.844-.25z'/%3e%3cpath fill='%2385340a' d='M393.75 262.72c-.248.003-.519.005-.813.031 4.488.428 2.331 3 7.032 3h.03c4.702 0 2.575-2.572 7.063-3-4.7-.426-3.214 2.344-7.062 2.344h-.031c-3.608 0-2.496-2.421-6.22-2.375zm10.1 6.94a3.848 3.848 0 0 0-3.846-3.846 3.848 3.848 0 0 0-3.847 3.846 3.955 3.955 0 0 1 3.847-3.04 3.952 3.952 0 0 1 3.846 3.04z'/%3e%3cpath id='e' fill='%2385340a' d='M382.73 244.02c4.915-4.273 11.11-4.915 14.53-1.709.837 1.121 1.373 2.32 1.593 3.57.43 2.433-.33 5.062-2.236 7.756.215-.001.643.212.856.427 1.697-3.244 2.297-6.577 1.74-9.746a13.815 13.815 0 0 0-.67-2.436c-4.7-3.845-11.11-4.272-15.81 2.138z'/%3e%3cpath id='d' fill='%2385340a' d='M390.42 242.74c2.777 0 3.419.642 4.7 1.71 1.284 1.068 1.924.854 2.137 1.068.213.215 0 .854-.426.64s-1.284-.64-2.564-1.708c-1.283-1.07-2.563-1.069-3.846-1.069-3.846 0-5.983 3.205-6.41 2.991-.426-.214 2.137-3.632 6.41-3.632z'/%3e%3cuse xlink:href='%23h' transform='translate(-19.181)'/%3e%3ccircle id='f' cx='390.54' cy='246.15' r='1.923' fill='%2385340a'/%3e%3cpath id='g' fill='%2385340a' d='M385.29 247.44c3.633 2.778 7.265 2.564 9.402 1.282 2.136-1.282 2.136-1.709 1.71-1.709-.427 0-.853.427-2.564 1.281-1.71.856-4.273.856-8.546-.854z'/%3e%3c/svg%3e)