Search

This work considers the impact of resolution in the construction of mock observations of simulated galaxies. In particular, when building mock integral field spectroscopic observations from galaxy formation models in cosmological simulations, we investigate the possible systematics that may arise given the assumption that all galaxies above some stellar mass limit will provide unbiased and meaningful observable stellar kinematics. We build a catalogue of N-body simulations to sample the range of stellar particle resolutions within the EAGLE Ref0050N0752 simulation box and examine how their observable kinematics vary relative to a higher-resolution N-body control. We use these models to compile a table of the minimum number of particles-per-pixel to reach a given uncertainty in the fitted line-of-sight velocity distribution parameters. Further, we introduce a Voronoi-binning module to the mock observation code, SimSpin, in order to meet these minimum numbers. Using EAGLE, we show the impact of this shot noise on the observed spin-ellipticity plane and the recovery of this space when observations are binned with increasing numbers of particles. In conclusion, we advise binning mock images to meet at least 200 particles-per-pixel to avoid systematically under-estimating the velocity dispersion along a given line-of-sight. We demonstrate that this is important for comparing galaxies extracted from the same simulation, as well as between simulations of varying mass resolution and observations of real galaxies., Comment: 19 pages, 11 figures (including 3 pages of appendix with 2 figures). Accepted 2024 November 3 for publication in MNRAS. Received 2024 October 28; in original form 2024 August 9

We present spectroscopic properties of 22 Lyman-alpha emitters(LAEs) at z=5.5-6.6 with Lyman-alpha(Lya) luminosity log($L_{Lya}$[$ergs^{-1}$])=42.4-43.5, obtained using VLT/MUSE as part of the Middle Ages Galaxy Properties with Integral Field Spectroscopy(MAGPI) survey. Additionally, we incorporate broad-band photometric data from the Subaru Hyper Suprime-Cam(HSC) for 17 LAEs in our sample. The HSC-y band magnitudes show that our LAEs are UV-bright, with rest-frame absolute UV magnitudes -19.7 < $M_{UV}$ < -23.3. We find that the Lya line width increases with luminosity, and this trend becomes more prominent at z > 6 where Lya lines become significantly broadened (> 260 $kms^{-1}$) at luminosities log($L_{Lya}$[$ergs^{-1}$]) > 43. This broadening is consistent with previous studies, suggesting that these sources are located inside larger ionized bubbles. We observe a slightly elevated ionizing photon production efficiency estimated for LAEs at z > 6, indicating that younger galaxies could be producing more ionizing photons per UV luminosity. A tentative anti-correlation between ionizing photon production efficiency and Lya rest-frame equivalent width is noticed, which could indicate a time delay between production and escape of ionizing photon primarily due to supernovae activity. Furthermore, we find a positive correlation between bubble radius and Lya line width, which again suggests that large ionized bubbles are created around these LAEs, allowing them to self-shield from the scattering effects of the intergalactic medium (IGM). We also detect two closely separated LAEs at z=6.046 (projected spatial separation is 15.92 kpc). The size of their respective bubbles suggests that they likely sit inside a common large ionized region. Such a closely-separated LAE pair increases the size of ionized bubble, potentially allowing a boosted transmission of Lya through neutral IGM. (Abridged), Comment: 14 pages, 6 figures, 2 tables, Accepted for publication in PASA

Studies of the internal mass structure of galaxies have observed a `conspiracy' between the dark matter and stellar components, with total (stars $+$ dark) density profiles showing remarkable regularity and low intrinsic scatter across various samples of galaxies at different redshifts. Such homogeneity suggests the dark and stellar components must somehow compensate for each other in order to produce such regular mass structures. We test the conspiracy using a sample of 22 galaxies from the `Middle Ages Galaxy Properties with Integral field spectroscopy' (MAGPI) Survey that targets massive galaxies at $ z \sim 0.3$. We use resolved, 2D stellar kinematics with the Schwarzschild orbit-based modelling technique to recover intrinsic mass structures, shapes, and dark matter fractions. This work is the first implementation of the Schwarzschild modelling method on a sample of galaxies at a cosmologically significant redshift. We find that the variability of structure for combined mass (baryonic and dark) density profiles is greater than that of the stellar components alone. Furthermore, we find no significant correlation between enclosed dark matter fractions at the half-light radius and the stellar mass density structure. Rather, the total density profile slope, $\gamma_{\mathrm{tot}}$, strongly correlates with the dark matter fraction within the half-light radius, as $\gamma_{\mathrm{tot}} = (1.3 \pm 0.2) f_{\mathrm{DM}} - (2.44 \pm 0.04)$. Our results refute the bulge-halo conspiracy and suggest that stochastic processes dominate in the assembly of structure for massive galaxies., Comment: 21 pages, 18 figures

Star formation rates (SFRs), gas-phase metallicities, and stellar masses are crucial for studying galaxy evolution. The different relations resulting from these properties give insights into the complex interplay of gas inside galaxies and their evolutionary trajectory and current characteristics. We aim to characterize these relations at $z\sim 0.3$, corresponding to a 3-4 Gyr lookback time. We utilized optical integral field spectroscopy of 65 emission-line galaxies from the MAGPI survey at a redshift of $0.28

We present a study of kinematic asymmetries from the integral field spectroscopic surveys MAGPI and SAMI. By comparing the asymmetries in the ionsied gas and stars, we aim to disentangle the physical processes that contribute to kinematic disturbances. We normalise deviations from circular motion by $S_{05}$, allowing us to study kinematic asymmetries in the stars and gas, regardless of kinematic temperature. We find a similar distribution of stellar asymmetries in galaxies where we do and do not detect ionised gas, suggesting that whatever is driving the stellar asymmetries does not always lead to gas removal. In both MAGPI and SAMI, we find an anti-correlation between stellar asymmetry and stellar mass, that is absent in the gas asymmetries. After stellar mass and mean-stellar-age matching distributions, we find that at all stellar masses, MAGPI galaxies display larger stellar asymmetry compared to SAMI galaxies. In both MAGPI and SAMI galaxies, we find that star-forming galaxies with old mean-stellar-ages typically have larger asymmetries in their gas compared to their stars, whereas galaxies with young mean-stellar-ages have larger asymmetries in their stars compared to their gas. We suggest that this results from continuous, clumpy accretion of gas., Comment: 12 pages, 10 figures

Since the turn of the century, astronomers have been exploiting the rich information afforded by combining stellar kinematic maps and imaging in an attempt to recover the intrinsic, three-dimensional (3D) shape of a galaxy. A common intrinsic shape recovery method relies on an expected monotonic relationship between the intrinsic misalignment of the kinematic and morphological axes and the triaxiality parameter. Recent studies have, however, cast doubt about underlying assumptions relating shape and intrinsic kinematic misalignment. In this work, we aim to recover the 3D shape of individual galaxies using their projected stellar kinematic and flux distributions using a supervised machine learning approach with mixture density network (MDN). Using a mock dataset of the EAGLE hydrodynamical cosmological simulation, we train the MDN model for a carefully selected set of common kinematic and photometric parameters. Compared to previous methods, we demonstrate potential improvements achieved with the MDN model to retrieve the 3D galaxy shape along with the uncertainties, especially for prolate and triaxial systems. We make specific recommendations for recovering galaxy intrinsic shapes relevant for current and future integral field spectroscopic galaxy surveys., Comment: Accepted for publication in PASA. 18 pages, 12 figures, 2 tables

We report the discovery of three double-peaked Lyman-$\alpha$ emitters (LAEs) exhibiting strong blue peak emission at 2.9 $\lesssim z \lesssim$ 4.8, in the VLT/MUSE data obtained as part of the Middle Ages Galaxy Properties with Integral Field Spectroscopy (MAGPI) survey. These strong blue peak systems provide a unique window into the scattering of Lyman-$\alpha$ photons by neutral hydrogen (HI), suggesting gas inflows along the line-of-sight and low HI column density. Two of them at $z=2.9$ and $z=3.6$ are spatially extended halos with their core regions clearly exhibiting stronger blue peak emissions than the red peak. However, spatial variations in the peak ratio and peak separation are evident over $25\times 26$ kpc ($z=2.9$) and $19\times28$ kpc ($z=3.6$) regions in these extended halos. Notably, these systems do not fall in the regime of Lyman-$\alpha$ blobs or nebulae. To the best of our knowledge, such a Lyman-$\alpha$ halo with a dominant blue core has not been observed previously. In contrast, the LAE at $z\sim4.8$ is a compact system spanning a $9\times9$ kpc region and stands as the highest-redshift strong blue peak emitter ever detected. The peak separation of the bright cores in these three systems ranges from $\Delta_{\mathrm{peak}}\sim370$ to $660$ km/s. The observed overall trend of decreasing peak separation with increasing radius is supposed to be controlled by HI column density and gas covering fraction. Based on various estimations, in contrast to the compact LAE, our halos are found to be good candidates for LyC leakers. These findings shed light on the complex interplay between Lyman-$\alpha$ emission, gas kinematics, and ionising radiation properties, offering valuable insights into the evolution and nature of high-redshift galaxies., Comment: 2 Figures, 1 Table, accepted for A&A Letters

In this work, we present a methodology and a corresponding code-base for constructing mock integral field spectrograph (IFS) observations of simulated galaxies in a consistent and reproducible way. Such methods are necessary to improve the collaboration and comparison of observation and theory results, and accelerate our understanding of how the kinematics of galaxies evolve over time. This code, SimSpin, is an open-source package written in R, but also with an API interface such that the code can be interacted with in any coding language. Documentation and individual examples can be found at the open-source website connected to the online repository. SimSpin is already being utilised by international IFS collaborations, including SAMI and MAGPI, for generating comparable data sets from a diverse suite of cosmological hydrodynamical simulations., Comment: 20 pages, 15 figures, 2 tables. Accepted for publication in PASA. 30/08/23

We investigate the impact of environment on the internal mass distribution of galaxies using the Middle Ages Galaxy Properties with Integral field spectroscopy (MAGPI) survey. We use 2D resolved stellar kinematics to construct Jeans dynamical models for galaxies at mean redshift $z \sim 0.3$, corresponding to a lookback time of $3-4$ Gyr. The internal mass distribution for each galaxy is parameterised by the combined mass density slope $\gamma$ (baryons $+$ dark matter), which is the logarithmic change of density with radius. We use a MAGPI sample of 28 galaxies from low-to-mid density environments and compare to density slopes derived from galaxies in the high density Frontier Fields clusters in the redshift range $0.29

We use comparisons between the SAMI Galaxy Survey and equilibrium galaxy models to infer the importance of disc fading in the transition of spirals into lenticular (S0) galaxies. The local S0 population has both higher photometric concentration and lower stellar spin than spiral galaxies of comparable mass and we test whether this separation can be accounted for by passive aging alone. We construct a suite of dynamically self--consistent galaxy models, with a bulge, disc and halo using the GalactICS code. The dispersion-dominated bulge is given a uniformly old stellar population, while the disc is given a current star formation rate putting it on the main sequence, followed by sudden instantaneous quenching. We then generate mock observables (r-band images, stellar velocity and dispersion maps) as a function of time since quenching for a range of bulge/total (B/T) mass ratios. The disc fading leads to a decline in measured spin as the bulge contribution becomes more dominant, and also leads to increased concentration. However, the quantitative changes observed after 5 Gyr of disc fading cannot account for all of the observed difference. We see similar results if we instead subdivide our SAMI Galaxy Survey sample by star formation (relative to the main sequence). We use EAGLE simulations to also take into account progenitor bias, using size evolution to infer quenching time. The EAGLE simulations suggest that the progenitors of current passive galaxies typically have slightly higher spin than present day star-forming disc galaxies of the same mass. As a result, progenitor bias moves the data further from the disc fading model scenario, implying that intrinsic dynamical evolution must be important in the transition from star-forming discs to passive discs., Comment: Accepted for publication in MNRAS

We present an overview of the Middle Ages Galaxy Properties with Integral Field Spectroscopy (MAGPI) survey, a Large Program on ESO/VLT. MAGPI is designed to study the physical drivers of galaxy transformation at a lookback time of 3-4 Gyr, during which the dynamical, morphological, and chemical properties of galaxies are predicted to evolve significantly. The survey uses new medium-deep adaptive optics aided MUSE observations of fields selected from the GAMA survey, providing a wealth of publicly available ancillary multi-wavelength data. With these data, MAGPI will map the kinematic and chemical properties of stars and ionised gas for a sample of 60 massive (> 7 x 10^10 M_Sun) central galaxies at 0.25 < z <0.35 in a representative range of environments (isolated, groups and clusters). The spatial resolution delivered by MUSE with Ground Layer Adaptive Optics (GLAO, 0.6-0.8 arcsec FWHM) will facilitate a direct comparison with Integral Field Spectroscopy surveys of the nearby Universe, such as SAMI and MaNGA, and at higher redshifts using adaptive optics, e.g. SINS. In addition to the primary (central) galaxy sample, MAGPI will deliver resolved and unresolved spectra for as many as 150 satellite galaxies at 0.25 < z <0.35, as well as hundreds of emission-line sources at z < 6. This paper outlines the science goals, survey design, and observing strategy of MAGPI. We also present a first look at the MAGPI data, and the theoretical framework to which MAGPI data will be compared using the current generation of cosmological hydrodynamical simulations including EAGLE, Magneticum, HORIZON-AGN, and Illustris-TNG. Our results show that cosmological hydrodynamical simulations make discrepant predictions in the spatially resolved properties of galaxies at z ~ 0.3. MAGPI observations will place new constraints and allow for tangible improvements in galaxy formation theory., Comment: 31 pages, 13 figures, PASA accepted

Observers experience a series of limitations when measuring galaxy kinematics, such as variable seeing conditions and aperture size. These effects can be reduced using empirical corrections, but these equations are usually applicable within a restrictive set of boundary conditions (e.g. S\'ersic indices within a given range) which can lead to biases when trying to compare measurements made across a full kinematic survey. In this work, we present new corrections for two widely used kinematic parameters, $\lambda_R$ and $V/\sigma$, that are applicable across a broad range of galaxy shapes, measurement radii and ellipticities. We take a series of mock observations of N-body galaxy models and use these to quantify the relationship between the observed kinematic parameters, structural properties and different seeing conditions. Derived corrections are then tested using the full catalogue of galaxies, including hydro-dynamic models from the EAGLE simulation. Our correction is most effective for regularly-rotating systems, yet the kinematic parameters of all galaxies -- fast, slow and irregularly rotating systems -- are recovered successfully. We find that $\lambda_R$ is more easily corrected than $V/\sigma$, with relative deviations of 0.02 and 0.06 dex respectively. The relationship between $\lambda_R$ and $V/\sigma$, as described by the parameter $\kappa$, also has a minor dependence on seeing conditions. These corrections will be particularly useful for stellar kinematic measurements in current and future integral field spectroscopic (IFS) surveys of galaxies., Comment: Accepted for publication in MNRAS 2020 June 20. 29 pages, 26 figures, including further supplementary material included online in Appendix D

We present SimSpin, a new, public, software framework for generating integral field spectroscopy (IFS) data cubes from N-body/hydrodynamical simulations of galaxies, which can be compared directly with observational datasets. SimSpin provides a consistent method for studying a galaxy's stellar component. It can be used to explore how observationally inferred measurements of kinematics, such as the spin parameter $\lambda_R$, are impacted by the effects of, for example, inclination, seeing conditions, distance, etc. SimSpin is written in R and has been designed to be highly modular, flexible, and extensible. It is already being used by the astrophysics community to generate IFS-like cubes and FITS files for direct comparison of simulations to observations. In this paper, we explain the conceptual framework of SimSpin; how it is implemented in R; and we demonstrate SimSpin's current capabilities, providing as an example a brief investigation of how numerical resolution affects how reliably we can recover the intrinsic stellar kinematics of a simulated galaxy., Comment: 14 pages, 10 figures. Accepted for publication in PASA

A primary goal of integral field spectroscopic (IFS) surveys is to provide a statistical census of galaxies classified by their internal kinematics. As a result, the observational spin parameter, $\lambda_R$, has become one of the most popular methods of quantifying the relative importance of velocity dispersion and rotation in supporting a galaxy's inner structure. The goal of this paper is to examine the relationship between the observationally deduced $\lambda_R$ and one of the most commonly used theoretical spin parameters in the literature, the Bullock et al. (2001) $\lambda'$. Using a set of $N$-body realisations of galaxies from which we construct mock IFS observations, we measure $\lambda_R$ as an observer would, incorporating the effects of beam smearing and seeing conditions. Assuming parameters typical of current IFS surveys, we confirm that there are strong positive correlations between $\lambda_R$ and measurement radius, and strong negative correlations between $\lambda_R$ and size of the PSF, for late-type galaxies; these biases can be reduced using a recently proposed empirical correction. Once observational biases are corrected for, we find that $\lambda_R$ provides a good approximation to $\sim \sqrt{3}/2 \; \lambda'(\rm R_{\rm eff})$, where $\lambda'$ is evaluated for the galactic stellar component within 1 R$_{\rm eff}$., Comment: Accepted for publication in MNRAS. 15 pages, 14 figures