Astrophysics of Galaxies (astro-ph.GA)

Abstract: The number density of extragalactic 21-cm radio sources as a function of their spectral line-widths -- the HI width function (HIWF) -- is a tracer of the dark matter halo mass function. The ALFALFA 21-cm survey measured the HIWF in northern and southern Galactic fields finding a systematically higher number density in the north; an asymmetry which is in tension with $\Lambda$ cold dark matter models which predicts the HIWF should be identical everywhere if sampled in sufficiently large volumes. We use the Sibelius-DARK N-body simulation and semi-analytical galaxy formation model GALFORM to create mock ALFALFA surveys to investigate survey systematics. We find the asymmetry has two origins: the sensitivity of the survey is different in the two fields, and the algorithm used for completeness corrections does not fully account for biases arising from spatial galaxy clustering. Once survey systematics are corrected, cosmological models can be tested against the HIWF.

Abstract: Observed scaling relations in galaxies between baryons and dark matter global properties are key to shed light on the process of galaxy formation and on the nature of dark matter. Here, we study the scaling relation between the neutral hydrogen (HI) and dark matter mass in isolated rotationally-supported disk galaxies at low redshift. We first show that state-of-the-art galaxy formation simulations predict that the HI-to-dark halo mass ratio decreases with stellar mass for the most massive disk galaxies. We then infer dark matter halo masses from high-quality rotation curve data for isolated disk galaxies in the local Universe, and report on the actual universality of the HI-to-dark halo mass ratio for these observed galaxies. This scaling relation holds for disks spanning a range of 4 orders of magnitude in stellar mass and 3 orders of magnitude in surface brightness. Accounting for the diversity of rotation curve shapes in our observational fits decreases the scatter of the HI-to-dark halo mass ratio while keeping it constant. This finding extends the previously reported discrepancy for the stellar-to-halo mass relation of massive disk galaxies within galaxy formation simulations to the realm of neutral atomic gas. Our result reveals that isolated galaxies with regularly rotating extended HI disks are surprisingly self-similar up to high masses, which hints at mass-independent self-regulation mechanisms that have yet to be fully understood.

Abstract: Quasars experiencing strong lensing offer unique viewpoints on subjects like the cosmic expansion rate, the dark matter profile within the foreground deflectors, and the quasar host galaxies. Unfortunately, identifying them in astronomical images is challenging since they are overwhelmed by the abundance of non-lenses. To address this, we have developed a novel approach by ensembling cutting-edge convolutional networks (CNNs) -- i.e., ResNet, Inception, NASNet, MobileNet, EfficientNet, and RegNet -- along with vision transformers (ViTs) trained on realistic galaxy-quasar lens simulations based on the Hyper Suprime-Cam (HSC) multiband images. While the individual model exhibits remarkable performance when evaluated against the test dataset, achieving an area under the receiver operating characteristic curve of $>$97.4% and a median false positive rate of 3.1%, it struggles to generalize in real data, indicated by numerous spurious sources picked by each classifier. A significant improvement is achieved by averaging these CNNs and ViTs, resulting in the impurities being downsized by factors up to 40. Subsequently, combining the HSC images with the UKIRT, VISTA, and unWISE data, we retrieve approximately 60 million sources as parent samples and reduce this to 892,609 after employing a photometry preselection to discover $z>1.5$ lensed quasars with Einstein radii of $\theta_\mathrm{E}<5$ arcsec. Afterward, the ensemble classifier indicates 3991 sources with a high probability of being lenses, for which we visually inspect, yielding 161 prevailing candidates awaiting spectroscopic confirmation. These outcomes suggest that automated deep learning pipelines hold great potential in effectively detecting strong lenses in vast datasets with minimal manual visual inspection involved.

Abstract: Context: The hierarchical model of galaxy formation suggests that galaxies are continuously growing. However, our position inside the Milky Way prevents us from studying the disk edge. Truncations are low surface brightness features located in the disk outskirts of external galaxies. They indicate where the disk brightness abruptly drops and their location is thought to change dynamically. In previous analyses of Milky Way-like galaxies, truncations were detected up to 3 kpc above the mid-plane but whether they remain present beyond that height remains unclear. Aims: Our goal is to determine whether truncations can be detected above 3 kpc height in the Milky Way-like galaxy NGC 4565, thus establishing the actual disk thickness. We also aim to study how the truncation relates to disk properties such as star formation activity or the warp. Methods: We perform a vertical study of the disk of NGC 4565 edge in unprecedented detail. We explore the truncation radius at different heights above/below the disk mid-plane (0<z<8 kpc) and at different wavelengths. We use new ultra-deep optical data ($\mu_{g,\rm{lim}}=30.5$ mag arcsec$^{-2}$; $3 \sigma$ within $10 \times 10$ arcsec$^{2}$ boxes) in the $g$, $r$ and $i$ broad bands, along with near- and far-ultraviolet, H$\alpha$, and \ion{H}{i} observations. Results: We detect the truncation up to 4 kpc in the $g$, $r$ and $i$ ultra-deep bands which is 1 kpc higher than in any previous study for any galaxy. The radial position of the truncation remains constant up to 3 kpc while higher up it is located at a smaller radius. This result is independent of the wavelength but is affected by the presence of the warp. Conclusions: We propose an inside-out growth scenario for the formation of the disk of NGC 4565. Our results point towards the truncation feature being linked to a star-forming threshold and to the onset of the disk warp.

Abstract: Many $z=1.5$ galaxies with a stellar mass ($M_{\star}$) $\geq 10^{10}\,\mathrm{M}_\odot$ are already quenched in both galaxy clusters ($>50$ per cent) and the field ($>20$ per cent), with clusters having a higher quenched fraction at all stellar masses compared to the field. A puzzling issue is that these massive quenched galaxies have stellar populations of similar age in both clusters and the field. This suggests that, despite the higher quenched fraction in clusters, the dominant quenching mechanism for massive galaxies is similar in both environments. In this work, we use data from the cosmological hydrodynamic simulations Hydrangea and EAGLE to test whether the excess quenched fraction of massive galaxies in $z = 1.5$ clusters results from fundamental differences in their halo properties compared to the field. We find that (i) at $10^{10} \leq$ $M_{\star}/\,\mathrm{M}_\odot\leq 10^{11}$, quenched fractions in the redshift range $1.5 < z < 3.5$ are consistently higher for galaxies with higher peak maximum circular velocity of the dark matter halo ($v_{\mathrm{max, peak}}$), and (ii) the distribution of $v_{\mathrm{max, peak}}$ is strongly biased towards higher values for cluster satellites compared to the field centrals. Due to this difference in the halo properties of cluster and field galaxies, secular processes alone may account for (most of) the environmental excess of massive quenched galaxies in high-redshift (proto) clusters. Taken at face value, our results challenge a fundamental assumption of popular quenching models, namely that clusters are assembled from an unbiased subset of infalling field galaxies. If confirmed, this would imply that such models must necessarily fail at high redshift, as indicated by recent observations.

Abstract: Extreme emission line galaxies (EELGs) are considered the primary contributor to cosmic reionization and are valuable laboratories to study the astrophysics of massive stars. It is strongly expected that Roman's High Latitude Wide Area Survey (HLWAS) will find many strongly gravitationally lensed [O III] emitters at Cosmic Noon (1 < z < 2.8). Roman imaging and grism spectroscopy alone will simultaneously confirm these strong lens systems and probe their interstellar medium (ISM) and stellar properties on small scales ($\lesssim$ 100 pc). Moreover, these observations will synergize with ground-based and space-based follow-up observations of the discovered lensed [O III] emitters in multi-wavelength analyses of their properties (e.g., massive stars and possible escape of ionizing radiation), spatially resolved on the scales of individual star cluster complexes. Only Roman can uniquely sample a large number of lensed [O III] emitters to study the small scale (~ 100 pc) ISM and stellar properties of these extreme emission line galaxies, detailing the key physics of massive stars and the ISM that govern cosmic reionization.

Abstract: The Astro2020 Decadal Survey has identified the mapping of the circumgalactic medium (CGM, gaseous plasma around galaxies) as a key objective. We explore the prospects for characterizing the CGM in and around nearby galaxy halos with future large grasp X-ray microcalorimeters. We create realistic mock observations from hydrodynamical simulations (EAGLE, IllustrisTNG, and Simba) that demonstrate a wide range of potential measurements, which will address the open questions in galaxy formation and evolution. By including all background and foreground components in our mock observations, we show why it is impossible to perform these measurements with current instruments, such as X-ray CCDs, and only microcalorimeters will allow us to distinguish the faint CGM emission from the bright Milky Way (MW) foreground emission lines. We find that individual halos of MW mass can, on average, be traced out to large radii, around R500, and for larger galaxies even out to R200, using the OVII, OVIII, or FeXVII emission lines. Furthermore, we show that emission line ratios for individual halos can reveal the radial temperature structure. Substructure measurements show that it will be possible to relate azimuthal variations to the feedback mode of the galaxy. We demonstrate the ability to construct temperature, velocity, and abundance ratio maps from spectral fitting for individual galaxy halos, which reveal rotation features, AGN outbursts, and enrichment.

Abstract: According to the current concordance cosmological model, the dark matter (DM) particles are collision-less and produce self-gravitating structures with a central cusp which, generally, is not observed. The observed density tends to a central plateau or core, explained within the cosmological model through the gravitational feedback of baryons on DM. This mechanism becomes inefficient when decreasing the galaxy stellar mass so that in the low-mass regime (Mstar << 10**6 Msun) the energy provided by the baryons is insufficient to modify cusps into cores. Thus, if cores exist in these galaxies they have to reflect departures from the collision-less nature of DM. Measuring the DM mass distribution in these faint galaxies is extremely challenging, however, their stellar mass distribution can be characterized through deep photometry. Here we provide a way of using only the stellar mass distribution to constrain the underlying DM distribution. The so-called Eddington inversion method allows us to discard pairs of stellar distributions and DM potentials requiring (unphysical) negative distribution functions in the phase space. In particular, cored stellar density profiles are incompatible with the Navarro, Frenk, and White (NFW) potential expected from collision-less DM if the velocity distribution is isotropic and the system spherically symmetric. Through a case-by-case analysis, we are able to relax these assumptions to consider anisotropic velocity distributions and systems which do not have exact cores. In general, stellar distributions with radially biased orbits are difficult to reconcile with NFW-like potentials, and cores in the baryon distribution tend to require cores in the DM distribution.

Abstract: We present the first results from the Revealing Low-Luminosity Active Galactic Nuclei (ReveaLLAGN) survey, a JWST survey of seven nearby LLAGN. We focus on two observations with the Mid-Infrared Instrument's (MIRI) Medium Resolution Spectrograph (MRS) of the nuclei of NGC 1052 and Sombrero (NGC 4594 / M104). We also compare these data to public JWST data of a higher-luminosity AGN, NGC 7319. JWST clearly resolves the AGN component even in Sombrero, the faintest target in our survey; the AGN components have very red spectra. We find that the emission-line widths in both NGC 1052 and Sombrero increase with increasing ionization potential, with FWHM > 1000 km/s for lines with ionization potential > 50 eV. These lines are also significantly blue-shifted in both LLAGN. The high ionization potential lines in NGC 7319 show neither broad widths or significant blue shifts. Many of the lower ionization potential emission lines in Sombrero show significant blue wings extending > 1000 km/s. These features and the emission-line maps in both galaxies are consistent with outflows along the jet direction. Sombrero has the lowest luminosity high-ionization potential lines ([Ne V] and [O IV]) ever measured in the mid-IR, but the relative strengths of these lines are consistent with higher luminosity AGN. On the other hand, the [Ne V] emission is much weaker relative to the [Ne III}] and [Ne II] lines of higher-luminosity AGN. These initial results show the great promise that JWST holds for identifying and studying the physical nature of LLAGN.

Abstract: The hot, X-ray-emitting phase of the circumgalactic medium in galaxies is believed to be the reservoir of baryons from which gas flows onto the central galaxy and into which feedback from AGN and stars inject mass, momentum, energy, and metals. These effects shape the velocity fields of the hot gas, which can be observed by X-ray IFUs via the Doppler shifting and broadening of emission lines. In this work, we analyze the gas kinematics of the hot circumgalactic medium of Milky Way-mass disk galaxies from the TNG50 simulation with synthetic observations to determine how future instruments can probe this velocity structure. We find that the hot phase is often characterized by outflows outward from the disk driven by feedback processes, radial inflows near the galactic plane, and rotation, though in other cases the velocity field is more disorganized and turbulent. With a spectral resolution of $\sim$1 eV, fast and hot outflows ($\sim$200-500 km s$^{-1}$) can be measured, depending on the orientation of the galaxy on the sky. The rotation velocity of the hot phase ($\sim$100-200 km s$^{-1}$) can be measured using line shifts in edge-on galaxies, and is slower than that of colder gas phases but similar to stellar rotation velocities. By contrast, the slow inflows ($\sim$50-100 km s$^{-1}$) are difficult to measure in projection with these other components. We find that the velocity measured is sensitive to which emission lines are used. Measuring these flows will help constrain theories of how the gas in these galaxies forms and evolves.

Abstract: We report new observations of the cool diffuse gas around 29, $2.3<z<6.3$ galaxies, using deep JWST/NIRCam slitless grism spectroscopy around the sightline to the quasar J0100+2802. The galaxies span a stellar mass range of $7.1 \leq \log M_{*}/M_{sun} \leq 10.7$, and star-formation rates of $-0.1 < \log \; SFR/M_{sun}yr^{-1} \; <2.3$. We find galaxies for seven MgII absorption systems within 300 kpc of the quasar sightline. The MgII radial absorption profile falls off sharply with radii, with most of the absorption extending out to 2-3$R_{200}$ of the host galaxies. Six out of seven MgII absorption systems are detected around galaxies with $\log M_{*}/M_{sun} >$9. MgII absorption kinematics are shifted from the systemic redshift of host galaxies with a median absolute velocity of 135 km/s and standard deviation of 85 km/s. The high kinematic offset and large radial separation ($R> 1.3 R_{200}$), suggest that five out of the seven MgII absorption systems are gravitationally not bound to the galaxies. In contrast, most cool circumgalactic media at $z<1$ are gravitationally bound. The high incidence of unbound MgII gas in this work suggests that towards the end of reionization, galaxy halos are in a state of remarkable disequilibrium, and are highly efficient in enriching the intergalactic medium. Two strongest MgII absorption systems are detected at $z\sim$ 4.22 and 4.5, the former associated with a merging galaxy system and the latter associated with three kinematically close galaxies. Both these galaxies reside in local galaxy over-densities, indicating the presence of cool MgII absorption in two "proto-groups" at $z>4$.

Abstract: We look for simulated star-forming linear wakes such as the one recently discovered by van Dokkum et al. (2023) in the cosmological hydrodynamical simulation ASTRID. Amongst the runaway black holes in ASTRID, none are able to produce clear star-forming wakes. Meanwhile, fly-by encounters, typically involving a compact galaxy (with a central black hole) and a star-forming galaxy (with a duo of black holes) reproduce remarkably well many of the key properties (its length and linearity; recent star formation, etc.) of the observed star-forming linear feature. We predict the feature to persist for approximately 100 Myr in such a system and hence constitute a rare event. The feature contains a partly stripped galaxy (with $M_{\rm gal}=10^9 \sim 10^{10}M_\odot$) and a dual BH system ($M_{\rm BH}=10^5 \sim 10^7\,M_\odot$) in its brightest knot. X-ray emission from AGN in the knot should be detectable in such systems. After $100\sim 200\,{\rm Myrs}$ from the first fly-by, the galaxies merge leaving behind a triple black hole system in a (still) actively star-forming early-type remnant of mass $\sim 5\times 10^{10}\,M_\odot$. Follow-up JWST observations may be key for revealing the nature of these linear features by potentially detecting the older stellar populations constituting the bright knot. Confirmation of such detections may therefore help discriminate a fly-by encounter from a massive BH wake to reveal the origin of such features.

Abstract: We derive predictions from state-of-the-art cosmological galaxy simulations for the spatial distribution of the hot circumgalactic medium (CGM, ${\rm [0.1-1]R_{200c}}$) through its emission lines in the X-ray soft band ($[0.3-1.3]$ keV). In particular, we compare IllustrisTNG, EAGLE, and SIMBA and focus on galaxies with stellar mass $10^{10-11.6}\, \MSUN$ at $z=0$. The three simulation models return significantly different surface brightness radial profiles of prominent emission lines from ionized metals such as OVII(f), OVIII, and FeXVII as a function of galaxy mass. Likewise, the three simulations predict varying azimuthal distributions of line emission with respect to the galactic stellar planes, with IllustrisTNG predicting the strongest angular modulation of CGM physical properties at radial range ${\gtrsim0.3-0.5\,R_{200c}}$. This anisotropic signal is more prominent for higher-energy lines, where it can manifest as X-ray eROSITA-like bubbles. Despite different models of stellar and supermassive black hole (SMBH) feedback, the three simulations consistently predict a dichotomy between star-forming and quiescent galaxies at the Milky-Way and Andromeda mass range, where the former are X-ray brighter than the latter. This is a signature of SMBH-driven outflows, which are responsible for quenching star formation. Finally, we explore the prospect of testing these predictions with a microcalorimeter-based X-ray mission concept with a large field-of-view. Such a mission would probe the extended hot CGM via soft X-ray line emission, determine the physical properties of the CGM, including temperature, from the measurement of line ratios, and provide critical constraints on the efficiency and impact of SMBH feedback on the CGM.

Abstract: The formation of molecular clouds out of HI gas is the first step toward star formation. Its metallicity dependence plays a key role to determine star formation through the cosmic history. Previous theoretical studies with detailed chemical networks calculate thermal equilibrium states and/or thermal evolution under one-zone collapsing background. The molecular cloud formation in reality, however, involves supersonic flows, and thus resolving the cloud internal turbulence/density structure in three dimension is still essential. We here perform magnetohydrodynamics simulations of 20 km s^-1 converging flows of Warm Neutral Medium (WNM) with 1 uG mean magnetic field in the metallicity range from the Solar (1.0 Zsun) to 0.2 Zsun environment. The Cold Neutral Medium (CNM) clumps form faster with higher metallicity due to more efficient cooling. Meanwhile, their mass functions commonly follow dn/dm proportional to m^-1.7 at three cooling times regardless of the metallicity. Their total turbulence power also commonly shows the Kolmogorov spectrum with its 80 percent in the solenoidal mode, while the CNM volume alone indicates the transition towards the Larson's law. These similarities measured at the same time in the unit of the cooling time suggest that the molecular cloud formation directly from the WNM alone requires a longer physical time in a lower metallicity environment in the 1.0-0.2 Zsun range. To explain the rapid formation of molecular clouds and subsequent massive star formation possibly within 10 Myr as observed in the Large/Small Magellanic Clouds (LMC/SMC), the HI gas already contains CNM volume instead of pure WNM.

Abstract: Using the SWIFT simulation code we study different forms of active galactic nuclei (AGN) feedback in idealized galaxy groups and clusters. We first present a physically motivated model of black hole (BH) spin evolution and a numerical implementation of thermal isotropic feedback (representing the effects of energy-driven winds) and collimated kinetic jets that they launch at different accretion rates. We find that kinetic jet feedback is more efficient at quenching star formation in the brightest cluster galaxies (BCGs) than thermal isotropic feedback, while simultaneously yielding cooler cores in the intracluster medium (ICM). A hybrid model with both types of AGN feedback yields moderate star formation rates, while having the coolest cores. We then consider a simplified implementation of AGN feedback by fixing the feedback efficiencies and the jet direction, finding that the same general conclusions hold. We vary the feedback energetics (the kick velocity and the heating temperature), the fixed efficiencies and the type of energy (kinetic versus thermal) in both the isotropic and the jet case. The isotropic case is largely insensitive to these variations. In particular, we highlight that kinetic isotropic feedback (used e.g. in IllustrisTNG) is similar in its effects to its thermal counterpart (used e.g. in EAGLE). On the other hand, jet feedback must be kinetic in order to be efficient at quenching. We also find that it is much more sensitive to the choice of energy per feedback event (the jet velocity), as well as the efficiency. The former indicates that jet velocities need to be carefully chosen in cosmological simulations, while the latter motivates the use of BH spin evolution models.