Search

We present a comprehensive study of the star formation histories of massive-quenched galaxies at $z=3$ in 3 semi-analytic models (SHARK, GAEA, GALFORM) and 3 cosmological hydrodynamical simulations (EAGLE, Illustris-TNG, Simba). We study the predicted number density and stellar mass function of massive-quenched galaxies, their formation and quenching timescales and star-formation properties of their progenitors. Predictions are disparate in all these diagnostics, for instance: (i) some simulations reproduce the observed number density of very massive-quenched galaxies ($>10^{11}\rm M_{\odot}$) but underpredict the high density of intermediate-mass ones, while others fit well the lower masses but underpredict the higher ones; (ii) In most simulations, except for GAEA and EAGLE, most massive-quenched galaxies had starburst periods, with the most intense ones happening at $4300\rm M_{\odot}\,yr^{-1}$; (iii) quenching timescales are in the range $\approx 20-150$~Myr depending on the simulation; among other differences. These disparate predictions can be tied to the adopted Active Galactic Nuclei (AGN) feedback model. For instance, the explicit black-hole (BH) mass dependence to trigger the "radio mode" in Illustris-TNG and Simba makes it difficult to produce quenched galaxies with intermediate stellar masses, also leading to higher baryon collapse efficiencies ($\approx 15-30$%); while the strong bolometric luminosity dependence of the AGN outflow rate in GAEA leads to BHs of modest mass quenching galaxies. Current observations are unable to distinguish between these different predictions due to the small sample sizes. However, these predictions are testable with current facilities and upcoming observations, allowing a "true physics experiment" to be carried out., Comment: Submitted for publication in MNRAS. 32 pages (21 of main body and 11 of appendices). Comments welcome!

Schwarzschild dynamical models are now regularly employed in large surveys of galaxies in the local and distant Universe to derive information on galaxies' intrinsic properties such as their orbital structure and their (dark matter and stellar) mass distribution. Comparing the internal orbital structures and mass distributions of galaxies in the distant Universe with simulations is key to understanding what physical processes are responsible for shaping galaxy properties. However it is first crucial to understand whether observationally derived properties are directly comparable with intrinsic ones in simulations. To assess this, we build Schwarzschild dynamical models for MUSE-like IFS cubes (constructed to be like those obtained by the MAGPI survey) of 75 galaxies at z ~ 0.3 from the Eagle simulations. We compare the true particle-derived properties with the galaxies' model-derived properties. In general, we find that the models can recover the true galaxy properties qualitatively well, with the exception of the enclosed dark matter, where we find a median offset of 48%, which is due to the assumed NFW profile not being able to reproduce the dark matter distribution in the inner region of the galaxies. We then compare our model-derived properties with Schwarzschild models-derived properties of observed MAGPI galaxies and find good agreement between MAGPI and Eagle: the majority of our galaxies (57%) have non-oblate shapes within 1 effective radius. More triaxial galaxies show higher fractions of hot orbits in their inner regions and tend to be more radially anisotropic., Comment: 22 pages. Accepted for publication in MNRAS

We measure the ionised gas velocity dispersions of star-forming galaxies in the MAGPI survey ($z\sim0.3$) and compare them with galaxies in the SAMI ($z\sim0.05$) and KROSS ($z\sim1$) surveys to investigate how the ionised gas velocity dispersion evolves. For the first time, we use a consistent method that forward models galaxy kinematics from $z=0$ to $z=1$. This method accounts for spatial substructure in emission line flux and beam smearing. We investigate the correlation between gas velocity dispersion and galaxy properties to understand the mechanisms that drive gas turbulence. We find that in both MAGPI and SAMI galaxies, the gas velocity dispersion more strongly correlates with the star-formation rate surface density ($\Sigma_{\rm SFR}$) than with a variety of other physical properties, and the average gas velocity dispersion is similar, at the same $\Sigma_{\rm SFR}$, for SAMI, MAGPI and KROSS galaxies. The results indicate that mechanisms related to $\Sigma_{\rm SFR}$ could be the dominant driver of gas turbulence from $z\sim1$ to $z\sim0$, for example, stellar feedback and/or gravitational instability. The gas velocity dispersion of MAGPI galaxies is also correlated with the non-rotational motion of the gas, illustrating that in addition to star-formation feedback, gas transportation and accretion may also contribute to the gas velocity dispersion for galaxies at $z\sim 0.3$. KROSS galaxies only have a moderate correlation between gas velocity dispersion and $\Sigma_{\rm SFR}$ and a higher scatter of gas velocity dispersion with respect to $\Sigma_{\rm SFR}$, in agreement with the suggestion that other mechanisms, such as gas transportation and accretion, are relatively more important at higher redshift galaxies., Comment: 15 pages, 10 figures, accepted for publication in MNRAS

We use the `Middle Ages Galaxy Properties with Integral field spectroscopy' (MAGPI) survey to investigate whether galaxies have evolved in the distribution of their stellar angular momentum in the past 3-4 Gyr, as probed by the observational proxy for spin, $\lambda_{R}$. We use 2D stellar kinematics to measure $\lambda_{R}$ along with detailed photometric models to estimate galaxy ellipticity. The combination of these measurements quantifies the kinematic classes of `fast rotators' and the rarer `slow rotators', which show no regular rotation in their line-of-sight velocity fields. We compare 51 MAGPI galaxies with $\log_{10} (M_{\star}/\mathrm{M}_\odot) > 10$ to carefully drawn samples of MaNGA galaxies in the local Universe, selected to represent possible descendants of the MAGPI progenitors. The EAGLE simulations are used to identify possible evolutionary pathways between the two samples, explicitly accounting for progenitor bias in our results and the varied evolutionary pathways a galaxy might take between the two epochs. We find that the occurrence of slow rotating galaxies is unchanged between the MAGPI ($z \sim 0.3$) and MaNGA ($z \sim 0$) samples, suggesting the massive slow rotator population was already in place $\sim 4$ Gyr ago and has not accumulated since. There is a hint of the MAGPI sample having an excess of high $\lambda_{R}$ galaxies compared to the MaNGA sample, corresponding to more ordered rotation, but statistically the samples are not significantly different. The large-scale stellar kinematics, as quantified through the $\lambda_{R}$ parameter, of galaxies at $z \sim 0.3$ have already evolved into the diversity of structures seen today in the local Universe., Comment: 10 pages, 7 figures, accepted for publication in MNRAS

Using adaptive optics with the Multi-Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT), the Middle Ages Galaxy Properties with Integral Field Spectroscopy (MAGPI) survey allows us to study the spatially resolved Universe at a crucial time of ~4 Gyr ago ($z$ ~ 0.3) when simulations predict the greatest diversity in evolutionary pathways for galaxies. We investigate the radial trends in the star formation (SF) activity and luminosity-weighted stellar ages as a function of offset from the star-forming main sequence (SFMS) for a total of 294 galaxies. Using both H$\alpha$ emission and the 4000 Angstrom break (i.e., D4000) as star formation rate (SFR) tracers, we find overall flat radial profiles for galaxies lying on and above the SFMS, suggestive of physical processes that enhance/regulate SF throughout the entire galaxy disc. However, for galaxies lying below the SFMS, we find positive gradients in SF suggestive of inside-out quenching. Placing our results in context with results from other redshift regimes suggests an evolution in radial trends at $z$ ~ 0.3 for SF galaxies above the SFMS, from uniformly enhanced SF at $z$ ~ 1 and $z$ ~ 0.3 to centrally enhanced SF at $z$ ~ 0 (when averaged over a wide range of mass). We also capture higher local SFRs for galaxies below the SFMS compared to that of $z$ ~ 0, which can be explained by a larger population of quenched satellites in the local Universe and/or different treatments of limitations set by the D4000-sSFR relation., Comment: 19 pages, 15 figures, 4 tables, accepted for publication in MNRAS

The merger history of a galaxy is thought to be one of the major factors determining its internal dynamics, with galaxies having undergone different types or mergers (e.g. dry, minor or major mergers) predicted to show different dynamical properties. We study the instantaneous orbital distribution of galaxies in the Eagle simulation, colouring the orbits of the stellar particles by their stellar age, in order to understand whether stars form in particular orbits (e.g. in a thin or thick disc). We first show that Eagle reproduces well the observed stellar mass fractions in different stellar orbital families as a function of stellar mass and spin parameter at z = 0. We find that the youngest stars reside in a thin disc component that can extend to the very inner regions of galaxies, and that older stars have warmer orbits, with the oldest ones showing orbits consistent with both hot and counter-rotating classifications, which is consistent with the trend found in the Milky-Way and other disc galaxies. We also show that counter-rotating orbits trace galaxy mergers - in particular dry mergers, and that in the absence of mergers, counter-rotating orbits can also be born from highly misaligned gas accretion that leads to star formation., Comment: Submitted to MNRAS, 21 pages

We use a sample of $z=0$ galaxies visually classified as slow rotators (SRs) in the EAGLE hydrodynamical simulations to explore the effect of galaxy mergers on their formation, characterise their intrinsic galaxy properties, and study the connection between quenching and kinematic transformation. SRs that have had major or minor mergers (mass ratios $\ge 0.3$ and $0.1-0.3$, respectively) tend to have a higher triaxiality parameter and ex-situ stellar fractions than those that had exclusively very minor mergers or formed in the absence of mergers ("no-merger" SRs). No-merger SRs are more compact, have lower black hole-to-stellar mass ratios and quenched later than other SRs, leaving imprints on their $z=0$ chemical composition. For the vast majority of SRs we find that quenching, driven by active galactic nuclei feedback, precedes kinematic transformation, except for satellite SRs, in which these processes happen in tandem. However, in $\approx 50$\% of these satellites, satellite-satellite mergers are responsible for their SR fate, while environment (i.e. tidal field and interactions with the central) can account for the transformation in the rest. By splitting SRs into kinematic sub-classes, we find that flat SRs prefer major mergers; round SRs prefer minor or very minor mergers; prolate SRs prefer gas-poor mergers. Flat and prolate SRs are more common among satellites hosted by massive halos ($>10^{13.6}\,\rm M_{\odot}$) and centrals of high masses ($M_{\star} > 10^{10.5}\, \rm M_{\odot}$). Although EAGLE galaxies display kinematic properties that broadly agree with observations, there are areas of disagreement, such as inverted stellar age and velocity dispersion profiles. We discuss these and how upcoming simulations can solve them., Comment: Accepted for publications in MNRAS. 22 pages and 20 figures. Main difference with previous version is additional panel in Fig. 7 showing Sersic index vs stellar mass for different merger histories

Large galaxy samples from multi-object IFS surveys now allow for a statistical analysis of the z~0 galaxy population using resolved kinematics. However, the improvement in number statistics comes at a cost, with multi-object IFS survey more severely impacted by the effect of seeing and lower S/N. We present an analysis of ~1800 galaxies from the SAMI Galaxy Survey and investigate the spread and overlap in the kinematic distributions of the spin parameter proxy $\lambda_{Re}$ as a function of stellar mass and ellipticity. For SAMI data, the distributions of galaxies identified as regular and non-regular rotators with \textsc{kinemetry} show considerable overlap in the $\lambda_{Re}$-$\varepsilon_e$ diagram. In contrast, visually classified galaxies (obvious and non-obvious rotators) are better separated in $\lambda_{Re}$ space, with less overlap of both distributions. Then, we use a Bayesian mixture model to analyse the observed $\lambda_{Re}$-$\log(M_*/M_{\odot})$ distribution. Below $\log(M_{\star}/M_{\odot})\sim10.5$, a single beta distribution is sufficient to fit the complete $\lambda_{Re}$ distribution, whereas a second beta distribution is required above $\log(M_{\star}/M_{\odot})\sim10.5$ to account for a population of low-$\lambda_{Re}$ galaxies. While the Bayesian mixture model presents the cleanest separation of the two kinematic populations, we find the unique information provided by visual classification of kinematic maps should not be disregarded in future studies. Applied to mock-observations from different cosmological simulations, the mixture model also predicts bimodal $\lambda_{Re}$ distributions, albeit with different positions of the $\lambda_{Re}$ peaks. Our analysis validates the conclusions from previous smaller IFS surveys, but also demonstrates the importance of using kinematic selection criteria that are dictated by the quality of the observed or simulated data., Comment: 30 pages and 17 figures, accepted for publication in MNRAS. Abstract abridged for Arxiv. The key figures of the paper are: 3, 7, 8, and 11