Galaxy formation and evolution

Abstract:

We present a comprehensive study of the star-forming main sequence (SFMS) and its scatter at redshifts , using Near Infrared Camera (NIRCam) photometry from the JADES (JWST Advanced Deep Extragalactic Survey) survey in the Great Observatories Origins Deep Survey (GOODS) South (GOODS-S) and North (GOODS-N) fields. Our analysis is based on a sample of galaxies that is stellar mass complete down to . The redshift evolution of the SFMS at an averaging time-scale of 10 Myr follows a relation, quantified by the specific star formation rates (sSFR), of with , in good agreement with theoretical predictions and the specific mass accretion rate of dark matter haloes. We find that the SFMS normalization varies in a complex way with the SFR averaging time-scale, reflecting the combined effects of bursty star formation and rising star formation histories (SFHs). We quantify the scatter of the SFMS, revealing that it decreases with longer SFR averaging time-scales, from at 10 Myr to at 100 Myr, indicating that shorter term fluctuations dominate the scatter, although long-term variations in star formation activity are also present. Our findings suggest that bursty SFHs are more pronounced at lower stellar masses. Furthermore, we explore the implications of our results for the observed overabundance of UV-bright galaxies at , concluding that additional mechanisms, such as top-heavy initial mass functions, increased star formation efficiencies, or increased burstiness in star formation are needed to explain these observations. Finally, we emphasize the importance of accurate stellar mass completeness limits when fitting the SFMS, especially for galaxies with bursty SFHs.

Abstract:

We utilize James Webb Space Telescope (JWST) Mid Infrared Instrument (MIRI) integral field unit observations to investigate the behavior and excitation of H2 in the nearby Seyfert galaxies NGC 3081 and NGC 5506, both part of the Galactic Activity, Torus, and Outflow Survey (or GATOS). We compare population levels of the S(1) to S(8) rotational H2 emission lines visible to JWST/MIRI spectroscopy to models assuming local thermodynamic equilibrium (LTE), in order to estimate the column density and thermal scaling of the molecular gas. For the nuclear regions, we incorporate Very Large Telescope Spectrograph for INtegral Field Observations in the Near Infrared (or VLT/SINFONI) K-band observations to estimate population levels for available rovibrational H2 emission lines, and compare the resultant population curves to non-LTE radiative transfer models and shock modeling. We report a differing set of prominent active galactic nuclei (AGN)-driven excitation mechanisms between the two galaxies. For NGC 3081, we find that a non-LTE radiative transfer environment is adequate to explain observations of the nuclear region, indicating that the primary mode in which the AGN transfers excitation energy is likely irradiation. We estimate the extent of AGN photoionization along the ionization bicone to be ≈330 pc. In contrast, for NGC 5506, we find a shock scenario to be a more plausible excitation mechanism, a conclusion bolstered by an observed spatial correlation between higher-energy rotational H2 and [Fe II]5.34μm emission. In addition, we identify potential nuclear H2 outflows resulting from an interaction between the ionization bicone and the rotational disk. By isolating the outflowing component of the H2 emission, we estimate the warm molecular mass outflow rate to be 0.07 M⊙ yr−1.

Abstract:

To maximize the amount of information extracted from cosmological datasets, simulations that accurately represent these observations are necessary. However, traditional simulations that evolve particles under gravity by estimating particle-particle interactions (𝑁-body simulations) are computationally expensive and prohibitive to scale to the large volumes and resolutions necessary for the upcoming datasets. Moreover, modeling the distribution of galaxies typically involves identifying virialized dark matter halos, which is also a time- and memory-consuming process for large 𝑁-body simulations, further exacerbating the computational cost. In this study, we introduce CHARM, a novel method for creating mock halo catalogs by matching the spatial, mass, and velocity statistics of halos directly from the large-scale distribution of the dark matter density field. We develop multistage neural spline flow-based networks to learn this mapping at redshift 𝑧 =0.5 directly with computationally cheaper low-resolution particle mesh simulations instead of relying on the high-resolution 𝑁-body simulations. We show that the mock halo catalogs and painted galaxy catalogs have the same statistical properties as obtained from 𝑁-body simulations in both real space and redshift space. Finally, we use these mock catalogs for cosmological inference using redshift-space galaxy power spectrum, bispectrum, and wavelet-based statistics using simulation-based inference, performing the first inference with accelerated forward model simulations and finding unbiased cosmological constraints with well-calibrated posteriors.

Abstract:

The distribution of molecular gas on small scales regulates star formation and the growth of supermassive black holes in galaxy centers. Yet, the role of active galactic nuclei (AGN) feedback in shaping this distribution remains poorly constrained. We investigate how AGNs influence the small-scale structure of molecular gas in galaxy centers by measuring the clumpiness of CO($3-2$) emission observed with the Atacama Large Millimeter/submillimeter Array (ALMA) in the nuclear regions ($50-200$ pc from the AGNs) of 16 nearby Seyfert galaxies from the Galaxy Activity, Torus, and Outflow Survey (GATOS). To quantify clumpiness we applied three different methods: (1) the median of the pixel-by-pixel contrast between the original and smoothed maps; (2) the ratio of the total excess flux to the total flux, after subtracting the background smoothed emission; and (3) the fraction of total flux coming from clumpy regions, interpreted as the mass fraction in clumps. We find a negative correlation between molecular gas clumpiness and AGN X-ray luminosity (L_ X ), suggesting that higher AGN activity is associated with smoother gas distributions. All methods reveal a turnover in this relation around L_ X erg s^-1, possibly indicating a threshold above which AGN feedback becomes efficient at dispersing dense molecular structures and suppressing future star formation. Our findings provide new observational evidence that AGN feedback can smooth out dense gas structures in galaxy centers.