Patricia Sánchez Blázquez

Using VLT/MUSE data, we study the ionised-gas kinematics in a mass- and volume-limited (M* ≥ 1010M⊙, D ≤ 20 Mpc) sample of 21 nearby galaxies to identify signatures of extended shocks within their inner kiloparsec, which appear as coherent velocity jumps in kinematic maps. By removing angular momentum, shocks in gas cause inflows, which can trigger nuclear star formation and fuel AGN activity. To identify the signatures of extended shocks, we examine residual velocity fields after subtracting a modelled rotating disc, and we study velocity difference between various gas tracers. Combining our kinematic analysis with BPT ionisation diagnostic maps and dust morphology, we find that 12 of 21 galaxies (~57 % show extended shock signatures with velocity jumps consistent with models of bar-driven shocks. This is likely a lower limit, as three additional galaxies (~14 %) exhibit shocks along bars, potentially reaching the nucleus but obscured by AGN outflows. We trace shock signatures inwards close to the resolution limit, which suggests that shocks may be the prevailing mechanism of inflow in the central kpc of galaxies. The only two unbarred galaxies in our sample are also the only systems with unperturbed kinematics and no shocks, strongly linking the perturbed gas dynamics in centres of galaxies to the presence of bars. All galaxies with inner bars show LINER- or Seyfert-like nuclear emission, whereas galaxies without inner bars exhibit all emission types, indicating that regardless of gas supply, inner bars suppress star formation in galactic nuclei.

We investigate the expected uncertainties involved when recovering galaxy cluster mass profiles from upcoming spectroscopic survey data, based on The Three Hundred Project. Using the caustic technique, which leverages galaxy positions and line-of-sight velocities, we assessed the systematic errors introduced by assumptions regarding velocity anisotropy and demonstrated how an iterative correction method can minimize these errors. We also assessed the impact of survey magnitude limits on cluster mass estimates, highlighting potential biases across different observational strategies We focused the analysis on our own CATARSIS survey, which is set to obtain redshift measurements for all galaxies with magnitudes of m_ AB, r < 22 within 2timesR_ 200,c for 16 galaxy clusters with redshifts of 0.14 < z < 0.27. We used the future 8,arcmin2 field-of-view TARSIS integral-field spectrograph of the Calar Alto 3.5-m telescope. These data will enable us to mitigate systematic errors in the determination of density profiles. CATARSIS is aimed at enhancing the precision of mass profile estimates by deepening our understanding of the dynamical states and physical characteristics of galaxy clusters.

Context. The new generation of optical spectrographs, including WEAVE, 4MOST, DESI, and WST, offers huge multiplexing capabilities and excellent spectral resolution. This is an unprecedented opportunity to statistically unveil the details of the star formation histories (SFHs) of galaxies. However, these observations are not easily comparable with the predictions of cosmological simulations. Aims. Our goal is to build a reference framework for comparing spectroscopic observations with cosmological simulations and to test the currently available tools for deriving the stellar population properties of mock galaxies as well as their SFHs. We focus on the observational strategy of the Stellar Population at intermediate redshift Survey (StePS) carried out with the WEAVE instrument. Methods. We created mock datasets of ∼750 galaxies at redshifts of z = 0.3, 0.5, and 0.7 from the TNG50 cosmological simulation. We performed radiative transfer calculations using SKIRT and analyzed the spectra with the pPXF algorithm, treating them as if they were real observations. Results. This work presents the methodology used to generate the mock datasets, providing an initial exploration of stellar population properties (i.e. mass-weighted ages and metallicities) and SFHs of a test sample of three galaxies at z = 0.7 and their descendants at z = 0.5 and 0.3. We show that there is very good agreement between mock WEAVE-like spectra compared to the intrinsic values in TNG50 (average difference of 0.2 ± 0.3 Gyr). We also report that an overall agreement is seen when retrieving the SFHs of galaxies, especially if they form the bulk of their stars on short timescales and at early epochs. While we did identify a tendency to overestimate the weight of old stellar populations in galaxies with complex SFHs, we were able to properly recover the timescales on which galaxies build up 90% of their mass, with almost no difference in the measured and intrinsic cumulative SFHs over the last 4 Gyr. Conclusions. We have released the datasets concurrently with this publication of this paper, which consist of multi-wavelength imaging and spectroscopic data of ∼750 galaxies at redshift z = 0.3, 0.5, and 0.7. This work provides a fundamental bench-test for forthcoming WEAVE observations, providing the community with realistic mock spectra of galaxies that can be used to test currently available tools for deriving first-order stellar populations parameters (i.e. ages and metallicities) as well as more complex diagnostics, such as mass and SFHs.

Context. Gas-phase abundances provide insights into the baryon cycle, with radial gradients and 2D metallicity distributions tracking how metals are built up and redistributed across galaxy disks over cosmic time. Aims. We use a catalog of 22 958 H II regions across 19 nearby spiral galaxies to examine how precisely the radial abundance gradients can be traced when using only the [N II ] λ 5755 electron temperature as a proxy for temperature-based, direct method metallicities. Methods. Using 534 direct detections of the temperature sensitive [N II ] λ 5755 auroral line, we measured gradients in 15 of the galaxies. Leveraging our large catalog of individual H II regions, we carried out a stacking procedure in bins of the H II region [N II ] λ 6583 luminosity and radius to recover stacked radial gradients. Results. We found a good agreement between the metallicity gradients from the stacked spectra and those gradients from individual regions and those from strong-line methods. In addition, particularly in the stacked T _e [N II ] measurements, some galaxies show very low (< 0.05 dex) scatter in metallicities, indicative of a well-mixed ISM. We examined the individual high confidence (S/N > 5) outliers and identified 13 regions across nine galaxies with anomalously low metallicities, although this is not strongly reflected in the strong-line method metallicities. By stacking arm and interarm regions, we found no systematic evidence for offsets in metallicity between these environments, suggesting that enrichment within spiral arms is due to very localized processes. Conclusions. This work demonstrates the potential to systematically exploit the single [N II ] λ 5755 auroral line for detailed gas-phase abundance studies of galaxies. It provides strong validation of previous results, based on the strong-line calibrations, of a well-mixed ISM across typical star-forming spiral galaxies.

We present measurements of the masses associated with ∼18,​​​​​000 H ii regions across 19 nearby star-forming galaxies by combining data from JWST, Hubble Space Telescope, MUSE, Atacama Large Millimeter/submillimeter Array, Very Large Array, and MeerKAT from the multiwavelength PHANGS survey. We report 10 pc-scale measurements of the mass of young stars, ionized gas, and older disk stars coincident with each H ii region, as well as the initial and current mass of molecular gas, atomic gas, and swept-up shell material, estimated from lower-resolution data. We find that the mass of older stars dominates over young stars at ≳10 pc scales, and ionized gas exceeds the stellar mass in most optically bright H ii regions. Combining our mass measurements for a statistically large sample of H ii regions, we derive 10 pc-scale star formation efficiencies of ≈6%–17% for individual H ii regions. Comparing each region’s self-gravity with the ambient interstellar medium (ISM) pressure and total pressure from presupernova stellar feedback, we show that most optically bright H ii regions are overpressured relative to their own self-gravity and the ambient ISM pressure and that they are hence likely expanding into their surroundings. Larger H ii regions in galaxy centers approach dynamical equilibrium. The self-gravity of regions is expected to dominate over presupernova stellar feedback pressure at ≳130 and 60 pc scales in galaxy disks and centers, respectively, but is always subdominant to the ambient ISM pressure on H ii region scales. Our measurements have direct implications for the dynamical evolution of star-forming regions and the efficiency of stellar feedback in ionizing and clearing cold gas.

Context. Stellar carbon abundances are crucial for tracing the star formation history and predicting the near-infrared emission of galaxies. It is still a complex task to derive accurate carbon abundance estimates for a wide variety of stars because it is hard to measure it based on atomic and molecular lines. It therefore remains challenging to include the abundance in stellar population models. Aims. We analyse the carbon abundance determination for the large empirical X-shooter Spectral Library that is commonly used as a benchmark for the development of stellar population models. Methods. We analysed the carbon abundance over strong molecular CH bands in the G-band region. We used the automated spectrum synthesis code GAUGUIN and adopted two different grids of separate reference synthetic spectra, each with the same coverage of the [C/Fe] abundance. We carried out a detailed comparison of the two grids to evaluate the accuracy and the model dependence of the measured [C/Fe] abundances. Results. We obtained a large and precise unbiased [C/Fe] abundance catalogue (~200 stars) from the two theoretical grids, well distributed in the Hertzsprung-Russell diagram and with no a trend with the stellar parameters. We also measured compatible values from each independent CH band, with a high-quality [C/Fe] abundance estimate for both dwarfs and giants. We always observed a dispersed flat trend around [C/Fe] ~ 0.0 dex throughout the entire metallicity regime we covered (–5 < [Fe/H] < +0.5 dex). This agrees with some previous studies. However, we found variations up to |Δ[C/Fe]| ~ 0.8 dex in the [C/Fe] composition of the star depending on the adopted grid. We found no such difference in the α-element measurements. This behaviour implies that the [C/Fe] abundance estimate strongly depends on the model. Conclusions. Potential sources of error might be associated with the use of spectral synthesis methods to derive stellar carbon abundances in the CH 4300 Å band. Intrinsic small differences in the synthetic models for this crowded and blended region may induce a large disparity in the precise abundance estimate for any stellar type, which leads to inadvertently inaccurate carbon measurements.

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.

The PHANGS project is assembling a comprehensive, multi-wavelength dataset of nearby (~5-20 Mpc), massive star-forming galaxies to enable multi-phase, multi-scale investigations into the processes that drive star formation and galaxy evolution. To date, large survey programs have provided molecular gas (CO) cubes with ALMA, optical IFU spectroscopy with VLT/MUSE, high-resolution NUV--optical imaging in five broad-band filters with HST, and infrared imaging in NIRCAM+MIRI filters with JWST. Here, we present PHANGS-HST-Halpha, which has obtained high-resolution (~2-10 pc), narrow-band imaging in the F658N or F657N filters with the HST/WFC3 camera of the warm ionized gas in the first 19 nearby galaxies observed in common by all four of the PHANGS large programs. We summarize our data reduction process, with a detailed discussion of the production of flux-calibrated, Milky Way extinction corrected, continuum-subtracted Halpha maps. PHANGS-MUSE IFU spectroscopy data are used to background subtract the HST-Halpha maps, and to determine the [NII] correction factors for each galaxy. We describe our public data products and highlight a few key science cases enabled by the PHANGS-HST-Halpha observations.