We analyse a set of cosmological magneto-hydrodynamic 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$\%$ 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 onto 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.

The magnetic fields of spiral galaxies are so strong that they cannot be primordial. Their typical values are over one billion times higher than any value predicted for the early Universe. Explaining this immense growth and incorporating it in galaxy evolution theories is one of the long-standing challenges in astrophysics. So far, the most successful theory for the sustained growth of the galactic magnetic field is the alpha-omega dynamo. This theory predicts a characteristic dipolar or quadrupolar morphology for the galactic magnetic field, which has been observed in external galaxies. However, so far, there has been no direct demonstration of a mean-field dynamo operating in direct, multi-physics simulations of spiral galaxies. We do so in this work. We employ numerical models of isolated, star-forming spiral galaxies that include a magnetized gaseous disk, a dark matter halo, stars, and stellar feedback. Naturally, the resulting magnetic field has a complex morphology that includes a strong random component. Using a smoothing of the magnetic field on small scales, we are able to separate the mean from the turbulent component and analyze them individually. We find that a mean-field dynamo naturally occurs as a result of the dynamical evolution of the galaxy and amplifies the magnetic field by an order of magnitude over half a Gyr. Despite the highly dynamical nature of these models, the morphology of the mean component of the field is identical to analytical predictions. This result underlines the importance of the mean-field dynamo in galactic evolution. Moreover, by demonstrating the natural growth of the magnetic field in a complex galactic environment, it brings us a step closer to understanding the cosmic origin of magnetic fields.

We explore the chemodynamical properties of a sample of barred galaxies in the Auriga magneto-hydrodynamical 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\sim3.5$: their last major merger occurs at $t_{\rm lookback}>12\,\rm Gyrs$, while subsequent mergers have a stellar mass ratio of $\leq$1:20, suggesting an upper limit of a few percent 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 ($<1\%$), 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% -- 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 5Gyrs (>30%) 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_{\phi}-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.

Stellar populations in barred galaxies save an imprint of the influence of the bar on the host galaxy's evolution. We present a detailed analysis of star formation histories (SFHs) and chemical enrichment of stellar populations in nine nearby barred galaxies from the TIMER project. We use integral field observations with the MUSE instrument to derive unprecedented spatially resolved maps of stellar ages, metallicities, [Mg/Fe] abundances and SFHs, as well as H$\alpha$ as a tracer of ongoing star formation. We find a characteristic V-shaped signature in the SFH perpendicular to the bar major axis which supports the scenario where intermediate age stars ($\sim 2$-$6\ \mathrm{Gyr}$) are trapped on more elongated orbits shaping a thinner part of the bar, while older stars ($> 8\ \mathrm{Gyr}$) are trapped on less elongated orbits shaping a rounder and thicker part of the bar. We compare our data to state-of-the-art cosmological magneto-hydrodynamical simulations of barred galaxies and show that such V-shaped SFHs arise naturally due to the dynamical influence of the bar on stellar populations with different ages and kinematic properties. Additionally, we find an excess of very young stars ($< 2\ \mathrm{Gyr}$) on the edges of the bars, predominantly on the leading side, confirming typical star formation patterns in bars. Furthermore, mass-weighted age and metallicity gradients are slightly shallower along the bar than in the disc likely due to orbital mixing in the bar. Finally, we find that bars are mostly more metal-rich and less [Mg/Fe]-enhanced than the surrounding discs. We interpret this as a signature that the bar quenches star formation in the inner region of discs, usually referred to as star formation deserts. We discuss these results and their implications on two different scenarios of bar formation and evolution.

Boxy, peanut- or X-shaped "bulges" are observed in a large fraction of barred galaxies viewed in, or close to, edge-on projection, as well as in the Milky Way. They are the product of dynamical instabilities occurring in stellar bars, which cause the latter to buckle and thicken vertically. Recent studies have found nearby galaxies that harbour two such features arising at different radial scales, in a nested configuration. In this paper we explore the formation of such double peanuts, using a collisionless N-body simulation of a pure disc evolving in isolation within a live dark matter halo, which we analyse in a completely analogous way to observations of real galaxies. In the simulation we find a stable double configuration consisting of two X/peanut structures associated to the same galactic bar - rotating with the same pattern speed - but with different morphology, formation time, and evolution. The inner, conventional peanut-shaped structure forms early via the buckling of the bar, and experiences little evolution once it stabilises. This feature is consistent in terms of size, strength and morphology, with peanut structures observed in nearby galaxies. The outer structure, however, displays a strong X, or "bow-tie", morphology. It forms just after the inner peanut, and gradually extends in time (within 1 to 1.5 Gyr) to almost the end of the bar, a radial scale where ansae occur. We conclude that, although both structures form, and are dynamically coupled to, the same bar, they are supported by inherently different mechanisms.

In this work we analyse the structural and photometric properties of 21 barred simulated galaxies from the Auriga Project. These consist of Milky Way-mass magneto-hydrodynamical simulations in a $\Lambda$CDM cosmological context. In order to compare with observations, we generate synthetic SDSS-like broad-band images from the numerical data at z = 0 with different inclinations (from face-on to edge-on). Ellipse fits are used to determine the bar lengths, and 2D bulge/disc/bar decompositions with galfit are also performed, modelling the bar component with the modified Ferrer profile. We find a wide range of bar sizes and luminosities in the sample, and their structural parameters are in good agreement with the observations. All bulges present low S\'ersic indexes, and are classified as pseudobulges. In regard to the discs, the same breaks in the surface brightness profiles observed in real galaxies are found, and the radii at which these take place are in agreement with the observations. Also, from edge-on unsharp-masked images at z = 0, boxy or peanut-shaped (B/P) structures are clearly identified in the inner part of 4 bars, and also 2 more bars are found in buckling phase. The sizes of the B/P match fairly well with those obtained from observations. We thus conclude that the observed photometric and structural properties of galaxies with bars, which are the main drivers of secular evolution, can be developed in present state-of-the-art $\Lambda$CDM cosmological simulations.

Using a set of 15 high-resolution magnetohydrodynamic cosmological simulations of Milky Way formation, we investigate the origin of the baryonic material found in stars at redshift zero. We find that roughly half of this material originates from subhalo/satellite systems and half is smoothly accreted from the Inter-Galactic Medium (IGM). About $90 \%$ of all material has been ejected and re-accreted in galactic winds at least once. The vast majority of smoothly accreted gas enters into a galactic fountain that extends to a median galactocentric distance of $\sim 20$ kpc with a median recycling timescale of $\sim 500$ Myr. We demonstrate that, in most cases, galactic fountains acquire angular momentum via mixing of low-angular momentum, wind-recycled gas with high-angular momentum gas in the Circum-Galactic Medium (CGM). Prograde mergers boost this activity by helping to align the disc and CGM rotation axes, whereas retrograde mergers cause the fountain to lose angular momentum. Fountain flows that promote angular momentum growth are conducive to smooth evolution on tracks quasi-parallel to the disc sequence of the stellar mass-specific angular momentum plane, whereas retrograde minor mergers, major mergers and bar-driven secular evolution move galaxies towards the bulge-sequence. Finally, we demonstrate that fountain flows act to flatten and narrow the radial metallicity gradient and metallicity dispersion of disc stars, respectively. Thus, the evolution of galactic fountains depends strongly on the cosmological merger history and is crucial for the chemo-dynamical evolution of Milky Way-sized disc galaxies.

Stellar feedback plays a significant role in modulating star formation, redistributing metals, and shaping the baryonic and dark structure of galaxies -- however, the efficiency of its energy deposition to the interstellar medium is challenging to constrain observationally. Here we leverage HST and ALMA imaging of a molecular gas and dust shell ($M_{H2} \sim 2\times 10^{5} ~{\rm M}_{\odot}$) in an outflow from the nuclear star forming ring of the galaxy NGC 3351, to serve as a boundary condition for a dynamical and energetic analysis of the outflowing ionised gas seen in our MUSE TIMER survey. We use \texttt{STARBURST99} models and prescriptions for feedback from simulations to demonstrate that the observed star formation energetics can reproduce the ionised and molecular gas dynamics -- provided a dominant component of the momentum injection comes from direct photon pressure from young stars, on top of supernovae, photoionisation heating and stellar winds. The mechanical energy budget from these sources is comparable to low luminosity AGN, suggesting that stellar feedback can be a relevant driver of bulk gas motions in galaxy centres - although here $\lesssim 10^{-3}$ of the ionized gas mass is escaping the galaxy. We test several scenarios for the survival/formation of the cold gas in the outflow, including in-situ condensation and cooling. Interestingly, the geometry of the molecular gas shell, observed magnetic field strengths and emission line diagnostics are consistent with a scenario where magnetic field lines aided survival of the dusty ISM as it was initially launched (with mass loading factor $\lesssim 1$) from the ring by stellar feedback. This system's unique feedback driven morphology can hopefully serve as a useful litmus test for feedback prescriptions in magnetohydrodynamical galaxy simulations.

We explore the origin of phase-space substructures revealed by the second Gaia data release in the disc of the Milky Way, such as the ridges in the $V_{\phi}$-$r$ plane, the undulations in the $V_{\phi}$-$r$-$V_r$ space and the streams in the $V_{\phi}$-$V_r$ plane. We use a collisionless N-body simulation with co-spatial thin and thick discs, along with orbit integration, to study the orbital structure close to the Outer Lindblad Resonance (OLR) of the bar. We find that a prominent, long-lived ridge is formed in the $V_{\phi}$-$r$ plane due to the OLR which translates to streams in the $V_{\phi}$-$V_r$ plane and examine which closed periodic and trapped librating orbits are responsible for these features. We find that orbits which carry out small librations around the $x_1(1)$ family are preferentially found at negative $V_r$, giving rise to a `horn'-like feature, while orbits with larger libration amplitudes, trapped around the $x_1(2)$ and $x_1(1)$ families, constitute the positive $V_r$ substructure, i.e. the Hercules-like feature. This changing libration amplitude of orbits will translate to a changing ratio of thin/thick disc stars, which could have implications on the metallicity distribution in this plane. We find that a scenario in which the Sun is placed close to the OLR gives rise to a strong asymmetry in $V_r$ in the $V_{\phi}$-$V_r$ plane (i.e. Hercules vs. `the horn') and subsequently to undulations in the $V_{\phi}$-$r$-$V_r$ space. We also explore a scenario in which the Sun is placed closer to the bar corotation and find that the bar perturbation $alone$ cannot give rise to the these features.

We study the galactic bulges in the Auriga simulations, a suite of thirty cosmological magneto-hydrodynamical 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 Sersic 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 percent of them with a negligible accreted fraction (facc < 0.01). In general,their in-situ component was centrally formed, with 75 percent 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 percent 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.

We present a convenient, all-in-one framework for the scientific analysis of fully reduced, (integral-field) spectroscopic data. The GIST pipeline (Galaxy IFU Spectroscopy Tool) is entirely written in Python3 and conducts all steps from the preparation of input data, over the scientific analysis to the production of publication-quality plots. In its basic setup, it extracts stellar kinematics, performs an emission-line analysis and derives stellar population properties from full spectral fitting as well as via the measurement of absorption line-strength indices by exploiting the well-known pPXF and GandALF routines, where the latter has now been implemented in Python. The pipeline is not specific to any instrument or analysis technique and provides easy means of modification and further development, as of its modular code architecture. An elaborate, Python-native parallelisation is implemented and tested on various machines. The software further features a dedicated visualization routine with a sophisticated graphical user interface. This allows an easy, fully-interactive plotting of all measurements, spectra, fits, and residuals, as well as star formation histories and the weight distribution of the models. The pipeline has successfully been applied to both low and high-redshift data from MUSE, PPAK (CALIFA), and SINFONI, as well as to simulated data for HARMONI@ELT and WEAVE and is currently being used by the TIMER, Fornax3D, and PHANGS collaborations. We demonstrate its capabilities by applying it to MUSE TIMER observations of NGC 1433.

Action space synthesizes the orbital information of stars and is well-suited to analyse the rich amount of kinematic disk substructure in the Gaia DR2 radial velocity sample (RVS). In this work, we revisit one of the strongest perturbers in the Milky Way (MW) disk: the m=2 bar. We investigate how its resonances affect the actions of individual test particle stars, i.e., (JR,Lz,Jz) estimated in an axisymmetric MW potential. We confirm that the stars' behaviour is well approximated by scattering and oscillation along a slope Delta JR / Delta Lz = l/m centered on the l:m resonance lines. The Outer Lindblad Resonance (OLR, l=+1,m=2) creates signatures in the stellar action space that can be used to identify the Galactic bar's OLR in the Gaia DR2 RVS data: (a) The JR dependence of the oscillation causes an overdensity ridge (underdensity region) at Lz larger (smaller) than the resonance line in the (Lz,JR) plane. (b) For the first time, we demonstrate that the OLR is expected to cause a gradient in average Jz with Lz across the resonance. (c) We show that the change of predominantly outward to inward motions at the OLR occurs along the resonance line in action space. The latter signature allows us to identify three candidates for the bar's OLR - and therefore its pattern speed Omega_bar - in the Gaia data within 3 kpc from the Sun: 1.85 Omega0, 1.2 Omega0, and 1.6 Omega0 (with ~0.1 Omega0 uncertainty). This demonstrates that (i) the local Gaia action data is consistent with both the short-fast and long-slow bar models in the literature, and that (ii) axisymmetrically estimated actions are a powerful diagnostic even in non-axisymmetric systems.

The formation of two stellar bars within a galaxy has proved challenging for numerical studies. It is yet not clear whether the inner bar is born via a star formation process promoted by gas inflow along the outer bar, or whether it is dynamically assembled from instabilities in a small-scale stellar disc. Observational constraints to these scenarios are scarce. We present a thorough study of the stellar content of two double-barred galaxies observed by the MUSE TIMER project, NGC 1291 and NGC 5850, combined with a two-dimensional multi-component photometric decomposition performed on the 3.6{\mu}m images from S4G. Our analysis confirms the presence of {\sigma}-hollows appearing in the stellar velocity dispersion distribution at the ends of the inner bars. Both galaxies host inner discs matching in size with the inner bars, suggestive of a dynamical formation for the inner bars from small-scale discs. The analysis of the star formation histories for the structural components shaping the galaxies provides constraints on the epoch of dynamical assembly of the inner bars, which took place >6.5 Gyr ago for NGC 1291 and >4.5 Gyr ago for NGC 5850. This implies that inner bars are long-lived structures.

In the MW bulge, metal-rich stars form a strong bar and are more peanut-shaped than metal-poor stars. It has been recently claimed that this behavior is driven by the initial in-plane radial velocity dispersion of these populations, rather than by their initial vertical random motions. This has led to the suggestion that a thick disc is not necessary to explain the characteristics of the MW bulge. We rediscuss this issue by analyzing two dissipationless N-body simulations of boxy/peanut (b/p)-shaped bulges formed from composite stellar discs, made of kinematically cold and hot stellar populations, and we conclude that initial vertical random motions are as important as in-plane random motions in determining the relative contribution of cold and hot disc populations with height above the plane, the metallicity and age trends. Previous statements emphasizing the dominant role of in-plane motions in determining these trends are not confirmed. However, differences exist in the morphology and strength of the resulting b/p-shaped bulges: a model where disc populations have initially only different in-plane random motions, but similar thickness, results into a b/p bulge where all populations have a similar peanut shape, independently on their initial kinematics, or metallicity. We discuss the reasons behind these differences, and also predict the signatures that these two extreme initial conditions would leave on the vertical age and metallicity gradients of disc stars, outside the bulge region. We conclude that a metal-poor, kinematically (radial and vertical) hot component, that is a thick disc, is necessary in the MW before bar formation, supporting the scenario traced in previous works. [abridged]

Double bars are thought to be important features for secular evolution in the central regions of galaxies. However, observational evidence about their origin and evolution is still scarce. We report on the discovery of the first Box-Peanut (B/P) structure in an inner bar detected in the face-on galaxy NGC 1291. We use the integral field data obtained from the MUSE spectrograph within the TIMER project. The B/P structure is detected as bi-symmetric minima of the $h_4$ moment of the line-of-sight velocity distribution along the major axis of the inner bar, as expected from numerical simulations. Our observations demonstrate that inner bars can follow a similar evolutionary path as outer bars, undergoing buckling instabilities. They also suggest that inner bars are long-lived structures, thus imposing tight constraints to their possible formation mechanisms

The Time Inference with MUSE in Extragalactic Rings (TIMER) project is a survey with the VLT-MUSE integral-field spectrograph of 24 nearby barred galaxies with prominent central structures (e.g., nuclear rings or inner discs). The main goals of the project are: (i) estimating the cosmic epoch when discs of galaxies settle, leading to the formation of bars; (ii) testing the hypothesis whereby discs in more massive galaxies are assembled first; and (iii) characterising the history of external gas accretion in disc galaxies. We present details on the sample selection, observations, data reduction, and derivation of high-level data products, including stellar kinematics, ages and metallicities. We also derive star formation histories and physical properties and kinematics of ionised gas. We illustrate how this dataset can be used for a plethora of scientific applications, e.g., stellar feedback, outflows, nuclear and primary bars, stellar migration and chemical enrichment, and the gaseous and stellar dynamics of nuclear spiral arms, barlenses, box/peanuts and bulges. Amongst our first results - based on a few selected galaxies -, we show that the dynamics of nuclear rings and inner discs is consistent with the picture in which they are formed by bars, that the central few hundred parsecs in massive disc galaxies tend to show a pronounced peak in stellar metallicity, and that nuclear rings can efficiently prevent star formation in this region. Finally, we present evidence that star-bursting nuclear rings can be fed with low-metallicity gas from low-mass companions.

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.

While the number of stars in the Galactic bulge with detailed chemical abundance measurements is increasingly rapidly, the inner Galactic bulge ( |b| < 2$^\circ$) remains poorly studied, due to heavy interstellar absorption and photometric crowding. We have carried out a high-resolution IR spectroscopic study of 72 M giants in the inner bulge using the CRIRES (ESO/VLT) facility. Our spectra cover the wavelength range of 2.0818 - 2.1444 $\mu$m with the resolution of R$\sim$50,000 and have signal-to-noise ratio of 50 - 100. Our stars are located along the bulge minor axis at l = 0$^\circ$, b = $\pm$0$^\circ$, $\pm$1$^\circ$, $\pm$2$^\circ$ and +3$^\circ$, selected with the aim of investigating any North-South asymmetries in the metallicity distribution function and composition, and comparing them to the outer bulge fields. Our sample was analysed in a homogeneous way using the most current K-band line list. We clearly detect a bimodal MDF with a metal-rich peak at $\rm \sim +0.3\,dex$ and a metal-poor peak at $\rm \sim -0.5\,dex$. Only a single star is found to exceed $\rm [Fe/H]=+0.5\,dex$. They show a symmetric behaviour along the $\pm$1$^\circ$, $\pm$2$^\circ$ fields. The Galactic Center field reveals in contrast a mainly metal-rich population with a mean metallicity of $\rm +0.3\,dex$. We derived $\rm [Mg/Fe]$ and $\rm [Si/Fe]$ abundances which are consistent with trends from the outer bulge, with a gradually decreasing trend with increasing metallicity. We confirm for the supersolar metallicity stars the decreasing trend in \mgfe\ and \sife\ as expected from chemical evolution models. With the caveat of a relatively small sample, we do not find significant differences in the chemical abundances between the Northern and the Southern fields, hence the evidence is consistent with symmetry in chemistry between North and South.

There is a long-standing debate on the origin of the metal-poor stellar populations of the Milky Way (MW) bulge, with the two leading scenarios being that these populations are either i) part of a classical metal-poor spheroid or ii) the same population as the chemically defined thick disc seen at the Solar neighbourhood. Here we test whether the latter scenario can reproduce the observed chemical properties of the MW bulge. To do so we compare an N-body simulation of a composite (thin+thick) stellar disc -- which evolves secularly to form a bar and a boxy/peanut (b/p) bulge -- to data from APOGEE DR13. This model, in which the thick disc is massive and centrally concentrated, can reproduce the morphology of the metal-rich and metal-poor stellar populations in the bulge, as well as the mean metallicity and [$\alpha$/Fe] maps as obtained from the APOGEE data. It also reproduces the trends, in both longitude and latitude, of the bulge metallicity distribution function (MDF). Additionally, we show that the model predicts small but measurable azimuthal metallicity variations in the inner disc due to the differential mapping of the thin and thick disc in the bar. We therefore see that the chemo-morphological relations of stellar populations in the MW bulge are naturally reproduced by mapping the thin and thick discs of the inner MW into a b/p.

Although there is consensus that metal-rich stars in the Milky Way bulge are formed via secular evolution of the thin disc, the origin of their metal-poor counterparts is still under debate. Two different origins have been invoked for metal-poor stars: they might be classical bulge stars or stars formed via internal evolution of a massive thick disc. We use N-body simulations to calculate the kinematic signature given by the difference in the mean Galactocentric radial velocity ($\Delta\rm V_{\rm GC}$) between metal-rich stars ([Fe/H] $\ge$ 0) and moderately metal-poor stars (-1.0 $\le$ [Fe/H] $<$ 0) in two models, one containing a thin disc and a small classical bulge (B/D=0.1), and the other containing a thin disc and a massive centrally concentrated thick disc. We reasonably assume that thin-disk stars in each model may be considered as a proxy of metal-rich stars. Similarly, bulge stars and thick-disc stars may be considered as a proxy of metal-poor stars. We calculate $\Delta\rm V_{\rm GC}$ at different latitudes ($b=0^\circ$, $-2^\circ$, $-4^\circ$,$-6^\circ$, $-8^\circ$ and $-10^\circ$) and longitudes ($l=0^\circ$, $\pm5^\circ$, $\pm10^\circ$ and $\pm15^\circ$) and show that the $\Delta\rm V_{\rm GC}$ trends predicted by the two models are different. We compare the predicted results with ARGOS data and APOGEE DR13 data and show that moderately metal-poor stars are well reproduced with the co-spatial stellar discs model, which has a massive thick disc. Our results give more evidence against the scenario that most of the metal-poor stars are classical bulge stars. If classical bulge stars exists, most of them probably have metallicities [Fe/H] $<$ -1 dex, and their contribution to the mass of the bulge should be a small percentage of the total bulge mass.

We show that the bulge and the disk of the Milky Way (MW) at R$\lesssim$7~kpc are well described by a unique chemical evolution and a two-phase star-formation history (SFH). We argue that the populations within this inner disk, not the entire disk, are the same, and that the outer Lindblad resonance (OLR) of the bar plays a key role in explaining this uniformity. In our model of a two-phase star formation history, the metallicity, [$\alpha$/Fe] and [$\alpha$/H] distributions, and age-metallicity relation are all compatible with the observations of both the inner disk and bulge. The dip at [Fe/H]$\sim$0 dex seen in the metallicity distributions of the bulge and inner disk reflects the quenching episode in the SFH of the inner MW at age $\sim$8 Gyr, and the common evolution of the bulge and inner disk stars. We show that at z$\le$1.5, when the MW was starting to quench, transitioning between the end of the $\alpha$-enhanced thick disk formation to the start of the thin disk, and yet was still gas rich, the gas accretion rate could not have been significant. The [$\alpha$/Fe] abundance ratio before and after this quenching phase would be different, which is not observed. The present analysis suggests that Milky Way history, and in particular at the transition from the thick to the thin disk -- the epoch of the quenching -- must have been driven by a decrease of the star formation efficiency. We argue that the decline in the intensity of gas accretion, the formation of the bar, and the quenching of the SFR at the same epoch may be causally connected thus explaining their temporal coincidence (abridged).

We examine the metallicity trends in the Milky Way (MW) bulge - using APOGEE DR13 data - and explore their origin by comparing two N-body models of isolated galaxies which develop a bar and a boxy/peanut (b/p) bulge. Both models have been proposed as scenarios for reconciling a disc origin of the MW bulge with a negative vertical metallicity gradient. The first model is a superposition of co-spatial, i.e. overlapping, disc populations with different scaleheights, kinematics and metallicities. In this model the thick, metal-poor, and centrally concentrated disc populations contribute significantly to the stellar mass budget in the inner galaxy. The second model is a single disc with an initial steep radial metallicity gradient, which is mapped by the bar into the b/p bulge in such a way that the vertical metallicity gradient of the MW bulge is reproduced -- as shown already in previous works in the literature. However, as we show here, the latter model does not reproduce the positive longitudinal metallicity gradient of the inner disc, nor the metal-poor innermost regions of the Bulge seen in the data. On the other hand, the model with co-spatial thin and thick disc populations reproduces all the aforementioned trends. We therefore see that it is possible to reconcile a (primarily) disc origin for the MW bulge with the observed trends in metallicity by mapping the inner thin and thick discs of the MW into a b/p. For this scenario to reproduce the observations, the $\alpha$-enhanced, metal-poor, thick disc populations must have a significant mass contribution in the inner regions -- as has been suggested for the Milky Way.

We explore trends in the morphology and line-of-sight (los) velocity of stellar populations in the inner regions of disc galaxies, using N-body simulations with both a thin (kinematically cold) and a thick (kinematically hot) disc which form a bar and boxy/peanut (b/p) bulge. The bar in the thin disc component is $\sim$50\% stronger than the thick disc bar and is more elongated, with an axis ratio almost half that of the thick disc bar. The thin disc b/p bulge has a pronounced X-shape, while the thick disc b/p is weaker with a rather boxy shape. This leads to the signature of the b/p bulge in the thick disc to be weaker and further away from the plane than in the thin disc. Regarding the kinematics, we find that the los velocity of thick disc stars in the outer parts of the b/p bulge can be \emph{larger} than that of thin disc stars, by up to 40\% and 20\% for side-on and Milky Way-like orientations of the bar respectively. This is due to the different orbits followed by thin and thick disc stars in the bar-b/p region, which are affected by the fact that: i) thin disc stars are trapped more efficiently in the bar - b/p instability and thus lose more angular momentum than their thick disc counterparts and ii) thick disc stars have large radial excursions and therefore stars from large radii with high angular momenta can be found in the bar region. We also find that the difference between the los velocities of the thin and thick disc in the b/p bulge ($\Delta v_{los}$) correlates with the initial difference between the radial velocity dispersions of the two discs ($\Delta \sigma$) . We therefore conclude that stars in the bar - b/p bulge will have considerably different morphologies and kinematics depending on the kinematic properties of the disc population they originate from.

We present a pilot study on the nearby massive barred galaxy NGC 1291, in which we use dynamical modelling to constrain the disc mass-to-light ratio (M/L), thus breaking the degeneracy between the baryonic and dark matter in its central regions. We use the gas, specifically the morphology of the dust lanes on the leading side of the bar, as a tracer of the underlying gravitational potential. We run a large number of hydrodynamic gas response simulations, in potentials obtained directly from near-infrared images of the galaxy, which have three free parameters: the M/L, the bar pattern speed and the height function. We explore the three-dimensional parameter space, by comparing the morphology of the shocks created in the gas response simulations with those of the observed dust lanes, and find the best-fitting models; these suggest that the M/L of NGC 1291 agrees with that predicted by stellar population synthesis models in the near-infrared ($\approx$0.6\,$M_{\odot}/L_{\odot}$), which leads to a borderline maximum disc for this galaxy. Furthermore, we find that the bar rotates fast, with a corotation radius which is $\leq$ 1.4 times the bar length.

In this letter we investigate the effect of boxy/peanut (b/p) bulges on bar-induced gas inflow to the central kiloparsec, which plays a crucial role on the evolution of disc galaxies. We carry out hydrodynamic gas response simulations in realistic barred galaxy potentials, including or not the geometry of a b/p bulge, to investigate the amount of gas inflow induced in the different models. We find that b/p bulges can reduce the gas inflow rate to the central kiloparsec by more than an order of magnitude, which leads to a reduction in the amount of gas available in the central regions. We also investigate the effect of the dark matter halo concentration on these results, and find that for maximal discs, the effect of b/p bulges on gas inflow remains significant. The reduced amount of gas reaching the central regions due to the presence of b/p bulges could have significant repercussions on the formation of discy- (pseudo-) bulges, on the amount of nuclear star formation and feedback, on the fuel reservoir for AGN activity, and on the overall secular evolution of galaxies.

We present a pilot study on the nearby massive galaxy NGC 1291, in which we aim to constrain the dark matter in the inner regions, by obtaining a dynamical determination of the disc mass-to-light ratio (M/L). To this aim, we model the bar-induced dust lanes in the galaxy, using hydrodynamic gas response simulations. The models have three free parameters, the M/L of the disc, the bar pattern speed and the disc height function. We explore the parameter space to find the best fit models, i.e. those in which the morphology of the shocks in the gas simulations matches the observed dust lanes. The best-fit models suggest that the M/L of NGC 1291 agrees with that predicted by stellar population synthesis models in the near-infrared ($\approx$0.6\,$M_{\odot}/L_{\odot}$), which leads to a borderline maximum disc for this galaxy. The bar rotates fast, with corotation radius $\leq$ 1.4 times the bar length. Additionally, we find that the height function has a significant effect on the results, and can bias them towards lower or higher M/L.

We examine the effects that the modelling of a Boxy/Peanut (B/P) bulge will have on the estimates of the stellar gravitational potential, forces, orbital structure and bar strength of barred galaxies. We present a method for obtaining the potential of disc galaxies from surface density images, assuming a vertical density distribution (height function), which is let to vary with position, thus enabling it to represent the geometry of a B/P. We construct a B/P height function after the results from a high-resolution, N-body+SPH simulation of an isolated galaxy and compare the resulting dynamical model to those obtained with the commonly used, position-independent "flat" height functions. We show that methods that do not allow for a B/P can induce errors in the forces in the bar region of up to 40% and demonstrate that this has a significant impact on the orbital structure of the model, which in turn determines its kinematics and morphology. Furthermore, we show that the bar strength is reduced in the presence of a B/P. We conclude that neglecting the vertical extent of a B/P can introduce considerable errors in the dynamical modelling. We also examine the errors introduced in the model due to uncertainties in the parameters of the B/P and show that even for generous but realistic values of the uncertainties, the error will be noticeably less than that of not modelling a B/P bulge at all.

In this paper we present an improved fitting formula for the dark matter bispectrum motivated by the previous phenomenological approach of Scoccimarro & Couchman (2001). We use a set of LCDM simulations to calibrate the fitting parameters in the k-range of 0.03 h/Mpc<k<0.4 h/Mpc and in the redshift range of 0<z<1.5. This new proposed fit describes well the BAO-features although it was not designed to. The deviation between the simulations output and our analytic prediction is typically less than 5% and in the worst case is never above 10%. We envision that this new analytic fitting formula will be very useful in providing reliable predictions for the non-linear dark matter bispectrum for LCDM models.