Adriano Poci

The central regions of disc galaxies host a rich variety of stellar structures: nuclear discs, bars, bulges, and boxy-peanut bulges. These components are often difficult to disentangle, both photometrically and kinematically, particularly in star-forming galaxies where dust obscuration and complex stellar motions complicate interpretation. In this work, we used data from the GECKOS-MUSE survey to investigate the impact of dust on axisymmetric Jeans Anisotropic Multi-Gaussian Expansion (JAM) models and assess their ability to recover kinematic structures in edge-on disc galaxies. We constructed JAM models for a sample of seven edge-on ( i ⪆ 85°) galaxies that span a range of star formation rates, dust content, and kinematic complexity. We find that when dust is appropriately masked, the disc regions of each galaxy are fit to χ _reduced ² ≤ 5. We analysed 2D residual velocity fields to identify signatures of non-axisymmetric structure. We find that derived dynamical masses are constant within 10% for each galaxy across all dust masking levels. In NGC 3957, a barred boxy galaxy in our sample, we identified velocity residuals that persist even under aggressive dust masking, aligned with bar orbits and supported by photometric bar signatures. We extended this analysis to reveal a bar in IC 1711 and a possible side-on bar in NGC 0522. Our results highlight both the capabilities and limitations of JAM in dusty, edge-on systems and attempt to link residual velocities to known non-axisymmetric kinematic structure.

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.

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 ∼ 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, γ_tot, strongly correlates with the dark matter fraction within the half-light radius, as γ_tot = (1.3 ± 0.2)f_DM − (2.44 ± 0.04). Our results refute the bulge-halo conspiracy and suggest that stochastic processes dominate in the assembly of structure for massive galaxies.