Astrophysics of Galaxies (astro-ph.GA)

Abstract: We select a disk-like galaxy sample with observations of the $HI$, $H_{2}$ and dust from Herschel Reference Survey (HRS), and derive inner HI masses within the optical radius. We find that the inner gas-to-dust ratio is almost independent of gas-phase metallicity, and confirm that the inner gas mass ($HI$+$H_{2}$) shows tighter relationship with dust mass and monochromatic 500 $\mu m$ luminosity than the integral gas mass. It supports that dust is more closely associated with co-spatial cold gas than the overall cold gas. Based on the newly calibrated relationship between inner gas mass and dust mass, we predict dust masses for disk-dominated galaxies from the xCOLD GASS sample. The predicted dust masses show scaling relations consistent with fiducial ones in the literature, supporting their robustness. Additionally, we find that at a given dust mass and star formation rate (SFR), the galactic WISE W3 luminosities show significant dependence on the [NII] luminosity and the stellar mass surface density. Such dependence highlights the caveat of using the W3 luminosity as integral SFR indicator, and is consistent with findings of studies which target star-forming regions in more nearby galaxies and accurately derive dust masses based on mapping-mode spectroscopy.

Abstract: We present an identification of dust-attenuated star-forming galactic-disk substructures in a typical star-forming galaxy (SFG), UDF2, at $z = 2.696$. To date, substructures containing significant buildup of stellar mass and actively forming stars have yet to be found in typical (i.e., main-sequence) SFGs at $z > 2$. This is due to the strong dust attenuation common in massive galaxies at the epoch and the scarcity of high-resolution, high-sensitivity extinction-independent imaging. To search for disk substructures, we subtracted the central stellar-mass disk from the JWST/NIRCam rest-frame 1.2 $\mu$m image ($0.13''$ resolution) and subtracted, in the visibility plane, the central starburst disk from ALMA rest-frame 240 $\mu$m observations ($0.03''$ resolution). The residual images revealed substructures at rest-frame 1.2 $\mu$m co-located with those found at rest-frame 240 $\mu$m, $\simeq 2$ kpc away from the galactic center. The largest substructure contains $\simeq20$% of the total stellar mass and $\simeq5$% of the total SFR of the galaxy. While UDF2 exhibits a kinematically-ordered velocity field of molecular gas consistent with a secularly evolving disk, more sensitive observations are required to characterize the nature and the origin of this substructure (spiral arms, minor merger, or other types of disk instabilities). UDF2 resides in an overdense region ($N \geqslant 4$ massive galaxies within 70 kpc projected distance at $z=2.690-2.697$) and the substructures may be associated with interaction-induced instabilities. Importantly, a statistical sample of such substructures identified with JWST and ALMA could play a key role in bridging the gap between the bulge-forming starburst and the rest of the galaxy.

Abstract: The formation of the first supermassive black holes is expected to have occurred in some most pronounced matter and galaxy overdensities in the early universe. We have conducted a sub-mm wavelength continuum survey of 54 $z\sim6$ quasars using the Submillimeter Common-User Bolometre Array-2 (SCUBA2) on the James Clerk Maxwell Telescope (JCMT) to study the environments around $z \sim 6$ quasars. We identified 170 submillimeter galaxies (SMGs) with above 3.5$\sigma$ detections at 450 or 850 \um\, maps. Their FIR luminosities are 2.2 - 6.4 $\times$ 10$^{12} L_{\odot}$, and star formation rates are $\sim$ 400 - 1200 M$_{\odot}$ yr$^{-1}$. We also calculated the SMGs differential and cumulative number counts in a combined area of $\sim$ 620 arcmin$^2$. To a $4\sigma$ detection (at $\sim$ 5.5 mJy), SMGs overdensity is $0.68^{+0.21}_{-0.19}$($\pm0.19$), exceeding the blank field source counts by a factor of 1.68. We find that 13/54 quasars show overdensities (at $\sim$ 5.5 mJy) of $\delta_{SMG}\sim$ 1.5 - 5.4. The combined area of these 13 quasars exceeds the blank field counts with the overdensity to 5.5 mJy of \dsmg $\sim$ $2.46^{+0.64}_{-0.55}$($\pm0.25$) in the regions of $\sim$ 150 arcmin$^2$. However, the excess is insignificant on the bright end (e.g., 7.5 mJy). We also compare results with previous environmental studies of Lyman alpha emitters (LAEs) and Lyman-Break Galaxies (LBGs) on a similar scale. Our survey presents the first systematic study of the environment of quasars at $z\sim6$. The newly discovered SMGs provide essential candidates for follow-up spectroscopic observations to test whether they reside in the same large-scale structures as the quasars and search for protoclusters at an early epoch.

Abstract: The CO-to-H$_2$ conversion factor ($\alpha_\rm{CO}$) is central to measuring the amount and properties of molecular gas. It is known to vary with environmental conditions, and previous studies have revealed lower $\alpha_\rm{CO}$ in the centers of some barred galaxies on kpc scales. To unveil the physical drivers of such variations, we obtained ALMA Band 3, 6, and 7 observations toward the inner 2 kpc of NGC 3627 and NGC 4321 tracing $^{12}$CO, $^{13}$CO, and C$^{18}$O lines on 100 pc scales. Our multi-line modeling and Bayesian likelihood analysis of these datasets reveal variations of molecular gas density, temperature, optical depth, and velocity dispersion, which are among the key drivers of $\alpha_\rm{CO}$. The central 300 pc nuclei in both galaxies show strong enhancement of temperature $T_\rm{k}>100$ K and density $n_\rm{H_2}>10^3$ cm$^{-3}$. Assuming a CO-to-H$_2$ abundance of $3\times10^{-4}$, we derive 4-15 times lower $\alpha_\rm{CO}$ than the Galactic value across our maps, which agrees well with previous kpc-scale measurements. Combining the results with our previous work on NGC 3351, we find a strong correlation of $\alpha_\rm{CO}$ with low-J $^{12}$CO optical depths ($\tau_\rm{CO}$), as well as an anti-correlation with $T_\rm{k}$. The $\tau_\rm{CO}$ correlation explains most of the $\alpha_\rm{CO}$ variation in the three galaxy centers, whereas changes in $T_\rm{k}$ influence $\alpha_\rm{CO}$ to second order. Overall, the observed line width and $^{12}$CO/$^{13}$CO 2-1 line ratio correlate with $\tau_\rm{CO}$ variation in these centers, and thus they are useful observational indicators for $\alpha_\rm{CO}$ variation. We also test current simulation-based $\alpha_\rm{CO}$ prescriptions and find a systematic overprediction, which likely originates from the mismatch of gas conditions between our data and the simulations.

Abstract: Feedback likely plays a crucial role in resolving discrepancies between observed and theoretical predictions of dwarf galaxy properties. Stellar feedback was once believed to be sufficient to explain these discrepancies, but it has thus far failed to fully reconcile theory and observations. The recent discovery of energetic galaxy-wide outflows in dwarf galaxies hosting Active Galactic Nuclei (AGN) suggests that AGN feedback may have a larger role in the evolution of dwarf galaxies than previously suspected. In order to assess the relative importance of stellar versus AGN feedback in these galaxies, we perform a detailed Keck/KCWI optical integral field spectroscopic study of a sample of low-redshift star-forming (SF) dwarf galaxies that show outflows in ionized gas in their SDSS spectra. We characterize the outflows and compare them to observations of AGN-driven outflows in dwarfs. We find that SF dwarfs have outflow components that have comparable widths (W$_{80}$) to those of outflows in AGN dwarfs, but are much less blue-shifted, indicating that SF dwarfs have significantly slower outflows than their AGN counterparts. The outflows in SF dwarfs are spatially resolved and significantly more extended than those in AGN dwarfs. The mass loss rates, momentum and energy rates of SF-driven outflows are much lower than those of AGN-driven outflows. Our results indicate that AGN feedback in the form of gas outflows may play an important role in dwarf galaxies and should be considered along with SF feedback in models of dwarf galaxy evolution.

Abstract: Weak emission-line quasars (WLQs) are a subset of Type 1 quasars that exhibit extremely weak Ly$\alpha +$N V $\lambda$1240 and/or C IV $\lambda$1549 emission lines. We investigate the relationship between emission-line properties and accretion rate for a sample of 230 `ordinary' Type 1 quasars and 18 WLQs at $z < 0.5$ and $1.5 < z < 3.5$ that have rest-frame ultraviolet and optical spectral measurements. We apply a correction to the H$\beta$-based black-hole mass ($M_{\rm BH}$) estimates of these quasars using the strength of the optical Fe II emission. We confirm previous findings that WLQs' $M_{\rm BH}$ values are overestimated by up to an order of magnitude using the traditional broad emission-line region size-luminosity relation. With this $M_{\rm BH}$ correction, we find a significant correlation between H$\beta$-based Eddington luminosity ratios and a combination of the rest-frame C IV equivalent width and C IV blueshift with respect to the systemic redshift. This correlation holds for both ordinary quasars and WLQs, which suggests that the two-dimensional C IV parameter space can serve as an indicator of accretion rate in all Type 1 quasars across a wide range of spectral properties.

Abstract: AGN have been generally considered to be less frequent in denser environments due to the lower number of galaxy-galaxy interactions and/or the removal of their gas-rich reservoirs by the dense intergalactic medium. However, recent observational and theoretical works suggest that the effect of ram-pressure stripping might reduce the angular momentum of their gas, causing it to infall towards the super massive black hole (SMBH) at their centre, activating the AGN phase. In this work we explore the connection between environment and nuclear activity by evaluating the variation in the incidence of ionized outflows in AGN across different environments. We select a sample of $\sim3300$ optical AGN from the Sloan Digital Sky Survey Data Release 13 that we match with the group catalogue from Lim et al. 2017. We further probe their environment through the projected distance to the central galaxy of the group/cluster and the projected surface density to the 5th neighbour ($\delta_5$). We find that at lower masses ($<10^{10.3}$M$_{\odot}$), the fraction of ionized outflows is significantly lower in satellite ($\sim7$%) than in isolated ($\sim22$%) AGN. The fraction of outflows in all satellite AGN decreases towards closer distances to the central, whereas only the lower-mass ones display a significant decline with $\delta_5$. Although this study does not include AGN in the densest regions of galaxy clusters, our findings suggest that AGN in dense environments accrete less gas than those in the field potentially due to the removal of the gas reservoirs via stripping or starvation, consistent with a negative connection between environment and AGN activity. We propose that the observed change in the incidence of outflows towards denser regions of groups and clusters could contribute to the higher gas metallicities of cluster galaxies compared to field ones, especially at lower masses.

Abstract: We present a detailed study of the column density and covering fraction profiles of C IV absorption around 86 redshift $z \approx 3.3$ Ly$\alpha$ emitters (LAEs) detected in 8 Multi-Unit Spectroscopic Explorer (MUSE) fields of $1'\times 1'$ centered on 8 bright background quasars as part of the MUSEQuBES survey. Using Voigt profile fitting of all the C IV absorbers detected along these 8 sightlines, we generated a ``blind'' absorbers' catalog consisting of 489 C IV absorption components. We cross-matched this blind C IV catalog with the MUSE-detected LAE catalog and found a significant enhancement of C IV components within $\approx \pm$400 $\rm km\, s^{-1}$ of the systemic redshifts of the LAEs. Neither the C IV column density ($N$) nor the Doppler parameter ($b$) of individual C IV components shows any significant anti-correlation with impact parameter ($\rho$) of the LAEs in the 68 percentile range of $90\leq \rho \leq 230$ physical kpc (pkpc). We find a covering fraction of $\approx 60\%$ for a threshold $N$(C IV) of $10^{12.5}\, \rm \rm cm^{-2}$, which is roughly twice as high as in random regions. The C IV covering fraction remains constant at $\approx50\%$ for impact parameters in the range 150--250~pkpc ($\approx 3-6 R_{200}$). Using the covering fraction profile, we constrained the LAE--C IV absorber two-point correlation function, and obtained $r_0 = 3.2~h^{-1}$ comoving Mpc (cMpc) and $\gamma = 1.2$ for a threshold $N$(C IV) of $10^{13.0}\, \rm cm^{-2}$. The C IV covering fraction is found to be enhanced for the LAEs that are part of a ``pair/group'' compared to the isolated ones.

Abstract: The direct comparison between hydrodynamical simulations and observations is needed to improve the physics included in the former and test biases in the latter. Post-processing radiative transfer and synthetic observations are now the standard way to do this. We report on the first application of the \texttt{SKIRT} radiative transfer code to simulations of a star-forming cloud. The synthetic observations are then analyzed following traditional observational workflows. We find that in the early stages of the simulation, stellar radiation is inefficient in heating dust to the temperatures observed in Galactic clouds, thus the addition of an interstellar radiation field is necessary. The spectral energy distribution of the cloud settles rather quickly after $\sim3$ Myr of evolution from the onset of star formation, but its morphology continues to evolve for $\sim8$ Myr due to the expansion of \textsc{Hii} regions and the respective creation of cavities, filaments, and ridges. Modeling synthetic \textit{Herschel} fluxes with 1- or 2-component modified black bodies underestimates total dust masses by a factor of $\sim2$. Spatially-resolved fitting recovers up to about $70\%$ of the intrinsic value. This ``missing mass'' is located in a very cold dust component with temperatures below $10$ K, which does not contribute appreciably to the far-infrared flux. This effect could bias real observations if such dust exists in large amounts. Finally, we tested observational calibrations of the SFR based on infrared fluxes and concluded that they are in agreement when compared to the intrinsic SFR of the simulation averaged over $\sim100$ Myr.