Luiz A. Silva-Lima

Stellar populations in the central regions of disc galaxies provide key constraints on their formation and secular evolution. We investigate the properties of young stellar populations in the centres of nearby barred galaxies from the TIMER survey, using four stellar population models, two based on empirical stellar libraries (MILES and X-Shooter) and two on theoretical libraries (Starburst99 and HR-pyPopStar), complemented by independent gas-phase metallicities from emission-line diagnostics. The inferred stellar metallicities depend strongly on the adopted model. Measured as the 5th–95th percentile spread across the three galaxies analysed in detail, the empirical libraries recover variations of 1.11 and 0.98 dex, whereas the theoretical ones return 0.40 and 0.49 dex; the gas in the same regions spans 0.49 dex. Gas metallicities are indeed lower in the central star-forming regions, but the amplitude of that depression (about 0.3 dex in NGC 1097) exceeds what bar-driven inflow of disc gas can account for, given typical abundance gradients in literature. The age–metallicity relations consistently reveal two distinct components; an older population consistent with a self-enrichment scenario, and a younger population formed from gas inflows with a different chemical enrichment history. The empirical libraries produce a monotonic decline of metallicities from old to young populations, which would require a sustained inflow of metal-poor gas outpacing its own enrichment. The theoretical libraries instead recover a young component whose metallicity departs from the older populations only at the most recent ages and by a few tenths of a dex, compatible with a single dilution episode. Our results demonstrate that model choice critically affects inferred metallicities, with models based on empirical stellar libraries recovering age–metallicity relations that demand physical conditions the gas measurements cannot sustain.

The vertical evolution of galactic discs is governed by the sub-structures within them. Several of these features, including bulges and kinematically distinct discs, are best studied in edge-on galaxies, as the viewing angle allows the easier separation of component light. For this work, we examined the diversity of kinematic sub-structure present in the first 12 galaxies observed from the GECKOS survey, a VLT/MUSE large programme providing a systematic study of 36 edge-on Milky Way-mass disc galaxies. Employing the NGIST analysis pipeline, we derived the mean luminosity-weighted line-of-sight stellar velocity (V_⋆), velocity dispersion (σ_⋆), skew (h₃), and kurtosis (h₄) for the sample, and examined 2D maps and 1D line profiles. Common clear kinematic signatures were observed: all galaxies display h₃ – V_⋆ sign mismatches in the outer disc regions consistent with a (quasi-)axisymmetric, rotating disc of stars. After scrutinising visual morphologies, we found that the majority of this sample (8/12) possess boxy-peanut bulges and host the corresponding kinematic structure predicted for stellar bars viewed in projection. Inferences were made on the bar viewing angle with respect to the line of sight from the strength of these kinematic indicators; we found one galaxy whose bar is close to side-on with respect to the observer, and two that are close to end-on. Four galaxies exhibit strong evidence for the presence of nuclear discs, including central h₃–V_⋆ profile anti-correlations, croissant-shaped central depressions in σ_⋆ maps, strong gradients in h₃, and positive h₄ plateaus over the expected nuclear disc extent. The strength of the h₃ feature corresponds to the size of the nuclear disc, measured from the h₃ turnover radius, taking into account geometric effects. We can explain the features within the kinematic maps of the four unbarred galaxies via disc structure(s) alone. We do not find any need to invoke the existence of dispersion-dominated bulges in any of the sample galaxies. Obtaining the specialised data products for this paper and the broader GECKOS survey required significant development of existing integral field spectroscopic (IFS) analysis tools. Therefore, we also present the NGIST pipeline: a modern, sophisticated, and easy-to-use pipeline for the analysis of galaxy IFS data, and the key tool employed by the GECKOS survey for producing value-added data products. We conclude that the variety of kinematic sub-structures seen in GECKOS galaxies requires a contemporary view of galaxy morphology, expanding on the traditional view of galaxy structure, and uniting the kinematic complexity observed in the Milky Way with the extragalactic.

Galaxy evolution is driven by spatially distributed processes with varying time-scales. Integral field spectroscopy provides spatially resolved information about these processes. Nevertheless, disentangling these processes, which are related to both the underlying stellar populations and the interstellar medium can be challenging. We present a case study on NGC 613, observed with MUSE (Multi-Unit Spectroscopic Explorer) for the TIMER (Time Inference with MUSE in Extragalactic Rings) project, a local barred galaxy, which shows several gas ionization mechanisms and is rich in both large and inner-scale stellar structures. We develop a set of steps to overcome fundamental problems in the modelling of emission lines with multiple components, together with the characterization of the stellar populations. That results in the disentanglement of the gas ionization mechanisms and kinematics, along with an optimal parametrization for star formation history recovery. Our analysis reveals evidence of gas inflows, which are associated with the bar dust lanes traced with Hubble Space Telescope. In addition, we show the gas kinematics in a central biconical outflow, which is aligned with a radio jet observed with Very Large Array. The emission line provides estimates of electron density, gas-phase metallicity, and the mass outflow rate, allowing us to distinguish intertwined ionization mechanisms and to identify a part of the multiphase gas cycle in NGC 613. It traces the gas kinematics from the bar lanes to inner scale gas reservoirs, where it can eventually trigger star formation or AGN activity, as observed in the outflow.

Galaxy evolution is driven by spatially distributed processes with varying timescales. Integral field spectroscopy provides spatially-resolved information about these processes. Nevertheless, disentangling these processes, which are related to both the underlying stellar populations and the interstellar medium can be challenging. We present a case study on NGC~613, observed with MUSE (Multi-Unit Spectroscopic Explorer) for the TIMER (Time Inference with MUSE in Extragalactic Rings) project, a local barred galaxy, which shows several gas ionisation mechanisms and is rich in both large and inner-scale stellar structures. We develop a set of steps to overcome fundamental problems in the modelling of emission lines with multiple components, together with the characterisation of the stellar populations. That results in the disentanglement of the gas ionisation mechanisms and kinematics, along with an optimal parametrisation for star formation history recovery. Our analysis reveals evidence of gas inflows, which are associated with the bar dust lanes traced with \textit{Hubble} Space Telescope (HST). In addition, we show the gas kinematics in a central biconical outflow, which is aligned with a radio jet observed with Very Large Array (VLA). The emission line provides estimates of electron density, gas-phase metallicity, and the mass outflow rate, allowing us to distinguish intertwined ionisation mechanisms and to identify a part of the multiphase gas cycle in NGC 613. It traces the gas kinematics from the bar lanes to inner scale gas reservoirs, where it can eventually trigger star formation or AGN activity, as observed in the outflow.

Bars drive gas inflow. As the gas flows inward, shocks and shear occur along the bar dust lanes. Such shocks and shear can affect the star formation (SF) and change the gas properties. For four barred galaxies, we present Hα velocity gradient maps that highlight bar-driven shocks and shear using data from the PHANGS-MUSE and PHANGS-ALMA surveys, which allow us to study bar kinematics in unprecedented detail. Velocity gradients are enhanced along the bar dust lanes, where shocks and shear are shown to occur in numerical simulations. Velocity gradient maps also efficiently pick up H ii regions that are expanding or moving relative to the surroundings. We put pseudo-slits on the regions where velocity gradients are enhanced and find that Hα and CO velocities jump up to ∼170 km s⁻¹, even after removing the effects of circular motions due to the galaxy rotation. Enhanced velocity gradients either coincide with the peak of CO intensity along the bar dust lanes or are slightly offset from CO intensity peaks, depending on the objects. Using the Baldwin–Philips–Terlevich BPT diagnostic, we identify the source of ionization on each spaxel and find that SF is inhibited in the high-velocity gradient regions of the bar, and the majority of those regions are classified as a low-ionization nuclear emission-line region (LINER) or composite. This implies that SF is inhibited where bar-driven shear and shocks are strong. Our results are consistent with the results from the numerical simulations that show SF is inhibited in the bar where the shear force is strong.

The high luminosity displayed by an active galactic nucleus (AGN) requires that gas be transported to the centre of the galaxy by some mechanism. Bar-driven processes are often pointed out in this context and a number of studies have addressed the bar–AGN connection, but with conflicting results. Some of the inconsistencies can be explained by the different spatial- and timescales involved in bar-driven gas inflows, accretion by the central black hole, and AGN emission. However, the discrepant results could also be due to sample biases, because both the AGN activity determination and the bar detection are influenced by the method employed. We revisit the bar–AGN connection in a sample of galaxies from SDSS, looking for evidence of the influence of bars on AGN activity. We determine AGN activity by emission line diagnostics and the properties of the bar were previously estimated withBUDDA, which performs 2D bulge–bar–disk decomposition. Before comparing active and inactive galaxies, we made a careful selection of the sample to minimise selection biases. We created control samples by matching them with the AGN sample using propensity score matching. This technique offers an analytical approach for creating control samples given some object parameters. We find that AGN are preferentially found in barred galaxies and that the accretion rate is higher in barred galaxies, but only when differentM–σrelations are used to estimate the black hole massM_•in barred and unbarred galaxies (from the central velocity dispersionσ). On the other hand, we find no correlation between activity level and bar strength. Altogether, our results strengthen theoretical predictions that the bar is an important mechanism in disc galaxies, creating a gas reservoir to feed AGN, but they also indicate that other mechanisms can play a major role, particularly at scales ≲100 pc.