A dynamical mechanism for the origin of nuclear rings

We develop a dynamical theory for the origin of nuclear rings in barred galaxies. In analogy with the standard theory of accretion discs, our theory is based on shear viscous forces among nested annuli of gas. However, the fact that gas follows non circular orbits in an external barred potential has profound consequences: it creates a region of reverse shear in which it is energetically favourable to form a stable ring which does not spread despite dissipation. Our theory allows us to approximately predict the size of the ring given the underlying gravitational potential. The size of the ring is loosely related to the location of the Inner Lindblad Resonance in the epicyclic approximation, but the predicted location is more accurate and is also valid for strongly barred potentials. By comparing analytical predictions with the results of hydrodynamical simulations, we find that our theory provides a viable mechanism for ring formation if the effective sound speed of the gas is low (\cs\lesssim1\kms), but that nuclear spirals/shocks created by pressure destroy the ring when the sound speed is high (\cs\simeq10\kms). We conclude that whether this mechanism for ring formation is relevant for real galaxies ultimately depends on the effective equation of state of the ISM. Promising confirmation comes from simulations in which the ISM is modelled using state-of-the-art cooling functions coupled to live chemical networks, but more tests are needed regarding the role of turbulence driven by stellar feedback. If the mechanism is relevant in real galaxies, it could provide a powerful tool to constrain the gravitational potential, in particular the bar pattern speed.Comment: Accepted for publication in MNRA

The dual origin of the Galactic thick disc and halo from the gas-rich Gaia-Enceladus Sausage merger

We analyse a set of cosmological magnetohydrodynamic simulations of the formation of Milky Way-mass galaxies identified to have a prominent radially anisotropic stellar halo component similar to the so-called 'Gaia Sausage' found in the Gaia data. We examine the effects of the progenitor of the Sausage (the Gaia-Enceladus Sausage, GES) on the formation of major galactic components analogous to the Galactic thick disc and inner stellar halo. We find that the GES merger is likely to have been gas-rich and contribute 10-50 per cent of gas to a merger-induced centrally concentrated starburst that results in the rapid formation of a compact, rotationally supported thick disc that occupies the typical chemical thick disc region of chemical abundance space. We find evidence that gas-rich mergers heated the proto-disc of the Galaxy, scattering stars on to less-circular orbits such that their rotation velocity and metallicity positively correlate, thus contributing an additional component that connects the Galactic thick disc to the inner stellar halo. We demonstrate that the level of kinematic heating of the proto-galaxy correlates with the kinematic state of the population before the merger, the progenitor mass, and orbital eccentricity of the merger. Furthermore, we show that the mass and time of the merger can be accurately inferred from local stars on counter-rotating orbits

Local variations of the stellar velocity ellipsoid - II. The effect of the bar in the inner regions of Auriga galaxies

Theoretical works have shown that off-plane motions of bars can heat stars in the vertical direction during buckling but is not clear how do they affect the rest of components of the stellar velocity ellipsoid (SVE). We study the 2D spatial distribution of the vertical, σz, azimuthal, σφ, and radial, σr velocity dispersions in the inner regions of Auriga galaxies, a set of high-resolution magneto-hydrodynamical cosmological zoom-in simulations, to unveil the influence of the bar on the stellar kinematics. σz and σφ maps exhibit non-axisymmetric features that closely match the bar light distribution with low-σ regions along the bar major axis and high values in the perpendicular direction. On the other hand, σr velocity dispersion maps present more axisymmetric distributions. We show that isophotal profile differences best capture the impact of the bar on the three SVE components providing strong correlations with bar morphology proxies although there is no relation with individual σ. Time evolution analysis shows that these differences are a consequence of the bar formation and that they tightly coevolve with the strength of the bar. We discuss the presence of different behaviours of σz and its connection with observations. This work helps us understand the intrinsic σ distribution and motivates the use of isophotal profiles as a mean to quantify the effect of bars

An ever-present Gaia snail shell triggered by a dark matter wake

We utilize a novel numerical technique to model star formation in cosmological simulations of galaxy formation – called superstars – to simulate a Milky Way-like galaxy with ≳108 star particles to study the formation and evolution of out-of-equilibrium stellar disc structures in a full cosmological setting. In the plane defined by the coordinate and velocity perpendicular to the mid-plane [vertical phase space, (Z, VZ)], stars in solar-like volumes at late times exhibit clear spirals qualitatively similar in shape and amplitude to the Gaia ‘snail shell’ phase spiral. We show that the phase spiral forms at a lookback time of ∼6 Gyr during the pericentric passage of an ∼1010 M. satellite on a polar orbit. This satellite stimulates the formation of a resonant wake in the dark matter halo while losing mass at a rate of ∼0.5–1 dex per orbit loop. The peak magnitude of the wake-induced gravitational torque at the solar radius is ∼8 times that from the satellite, and triggers the formation of a disc warp that wraps up into a vertical phase spiral over time. As the wake decays, the phase spiral propagates several gigayears to present day and can be described as ‘ever-present’ once stable disc evolution is established. These results suggest an alternative scenario to explain the Gaia phase spiral, which does not rely on a perturbation from bar buckling or a recent direct hit from a satellite

The prevalence of pseudo-bulges in the Auriga simulations

We study the galactic bulges in the Auriga simulations, a suite of 30 cosmological magnetohydrodynamical zoom-in simulations of late-type galaxies in Milky Way sized dark matter haloes performed with the moving-mesh code AREPO. We aim to characterize bulge formation mechanisms in this large suite of galaxies simulated at high resolution in a fully cosmological context. The bulges of the Auriga galaxies show a large variety in their shapes, sizes, and formation histories. According to observational classification criteria, such as Sérsic index and degree of ordered rotation, the majority of the Auriga bulges can be classified as pseudo-bulges, while some of them can be seen as composite bulges with a classical component; however, none can be classified as a classical bulge. Auriga bulges show mostly an in situ origin, 21 per cent of them with a negligible accreted fraction (facc < 0.01). In general, their in situ component was centrally formed, with ∼75 per cent of the bulges forming most of their stars inside the bulge region at z = 0. Part of their in situ mass growth is rapid and is associated with the effects of mergers, while another part is more secular in origin. In 90 per cent of the Auriga bulges, the accreted bulge component originates from less than four satellites. We investigate the relation between the accreted stellar haloes and the bulges of the Auriga simulations. The total bulge mass shows no correlation with the accreted stellar halo mass, as in observations. However, the accreted mass of bulges tends to correlate with their respective accreted stellar halo mass

The effects of AGN feedback on the structural and dynamical properties of Milky Way-mass galaxies in cosmological simulations

Feedback from active galactic nuclei (AGNs) has become established as a fundamental process in the evolution of the most massive galaxies. Its impact on Milky Way (MW)-mass systems, however, remains comparatively unexplored. In this work, we use the auriga simulations to probe the impact of AGN feedback on the dynamical and structural properties of galaxies, focusing on the bar, bulge, and disc. We analyse three galaxies - two strongly and one unbarred/weakly barred - using three setups: (i) the fiducial auriga model, which includes both radio and quasar mode feedback, (ii) a setup with no radio mode, and (iii) one with neither the radio nor the quasar mode. When removing the radio mode, gas in the circumgalactic medium cools more efficiently and subsequently settles in an extended disc, with little effect on the inner disc. Contrary to previous studies, we find that although the removal of the quasar mode results in more massive central components, these are in the form of compact discs, rather than spheroidal bulges. Therefore, galaxies without quasar mode feedback are more baryon-dominated and thus prone to forming stronger and shorter bars, which reveals an anticorrelation between the ejective nature of AGN feedback and bar strength. Hence, we report that the effect of AGN feedback (i.e. ejective or preventive) can significantly alter the dynamical properties of MW-like galaxies. Therefore, the observed dynamical and structural properties of MW-mass galaxies can be used as additional constraints for calibrating the efficiency of AGN feedback models