Search

We present results from Atacama Large Millimeter/submillimeter Array (ALMA) observations of CS from the nearby galaxy NGC 1068 ($\sim14$ Mpc). This Seyfert 2 barred galaxy possesses a circumnuclear disc (CND, $r\sim200$ pc) and a starburst ring (SB ring, $r\sim1.3$ kpc). These high-resolution maps ($\sim0.5$", $\sim35$ pc) allow us to analyse specific sub-regions in the galaxy and investigate differences in line intensity ratios and physical conditions, particularly those between the CND and SB ring. Local thermodynamic equilibrium (LTE) analysis of the gas is used to calculate CS densities in each sub-region, followed by non-LTE analysis conducted using the radiative transfer code RADEX to fit observations and constrain gas temperature, CS column density and hydrogen density. Finally, the chemical code UCLCHEM is used to reconstruct the gas, allowing an insight into its origin and chemical history. The density of hydrogen in the CND is found to be $\geq10^5$ cm$^{-2}$, although exact values vary, reaching $10^6$ cm$^{-2}$ at the AGN. The conditions in the two arms of the SB ring appear similar to one another, though the density found ($\sim10^4$ cm$^{-2}$) is lower than in the CND. The temperature in the CND increases from east to west, and is also overall greater than found in the SB ring. These modelling methods indicate the requirement for multi-phase gas components in order to fit the observed emission over the galaxy. A larger number of high resolution transitions across the SLED may allow for further constraining of the conditions, particularly in the SB ring., Comment: 32 pages, 19 figures, 6 tables, accepted for publication in MNRAS

Scaling laws of dust, HI gas and metal mass with stellar mass, specific star formation rate and metallicity are crucial to our understanding of the buildup of galaxies through their enrichment with metals and dust. In this work, we analyse how the dust and metal content varies with specific gas mass ($M_{\text{HI}}$/$M_{\star}$) across a diverse sample of 423 nearby galaxies. The observed trends are interpreted with a set of Dust and Element evolUtion modelS (DEUS) - incluidng stellar dust production, grain growth, and dust destruction - within a Bayesian framework to enable a rigorous search of the multi-dimensional parameter space. We find that these scaling laws for galaxies with $-1.0\lesssim \log M_{\text{HI}}$/$M_{\star}\lesssim0$ can be reproduced using closed-box models with high fractions (37-89$\%$) of supernova dust surviving a reverse shock, relatively low grain growth efficiencies ($\epsilon$=30-40), and long dus lifetimes (1-2\,Gyr). The models have present-day dust masses with similar contributions from stellar sources (50-80\,$\%$) and grain growth (20-50\,$\%$). Over the entire lifetime of these galaxies, the contribution from stardust ($>$90\,$\%$) outweighs the fraction of dust grown in the interstellar medium ($<$10$\%$). Our results provide an alternative for the chemical evolution models that require extremely low supernova dust production efficiencies and short grain growth timescales to reproduce local scaling laws, and could help solving the conundrum on whether or not grains can grow efficiently in the interstellar medium., Comment: 56 pages, 30 figures, accepted for publication in MNRAS

We use HI and H2 global gas measurements of galaxies from xGASS and xCOLD GASS to investigate quenching paths of galaxies below the star formation main sequence (SFMS). We show that the population of galaxies below the SFMS is not a 1:1 match with the population of galaxies below the HI and H2 gas fraction scaling relations. Some galaxies in the transition zone (TZ) 1-sigma below the SFMS can be as HI-rich as those in the SFMS, and have on average longer gas depletion timescales. We find evidence for environmental quenching of satellites, but central galaxies in the TZ defy simple quenching pathways. Some of these so-called "quenched" galaxies may still have significant gas reservoirs and be unlikely to deplete them anytime soon. As such, a correct model of galaxy quenching cannot be inferred with SFR (or other optical observables) alone, but must include observations of the cold gas. We also find that internal structure (particularly, the spatial distribution of old and young stellar populations) plays a significant role in regulating the star formation of gas-rich isolated TZ galaxies, suggesting the importance of bulges in their evolution., Comment: 15 pages, 11 figures, accepted for publication in MNRAS

We determine the cosmic abundance of molecular hydrogen (H2) in the local universe from the xCOLD GASS survey. To constrain the H2 mass function at low masses and correct for the effect of the lower stellar mass limit of 10^9 Msun in the xCOLD GASS survey, we use an empirical approach based on an observed scaling relation between star formation rate and gas mass. We also constrain the HI and HI+H2 mass functions using the xGASS survey, and compare it to the HI mass function from the ALFALFA survey. We find the cosmic abundance of molecular gas in the local Universe to be Omega_H2=(5.34+/-0.47)x10^-5 h^-1. Molecular gas accounts for 19.6 +/- 3.9% of the total abundance of cold gas, Omega_HI+H2=(4.66+/-0.70)x10^-4 h^-1. Galaxies with stellar masses in excess of 10^9 Msun account for 89% of the molecular gas in the local Universe, while in comparison such galaxies only contain 73% of the cold atomic gas as traced by the HI 21cm line. The xCOLD GASS CO, molecular gas and cold gas mass functions and Omega_H2 measurements provide constraints for models of galaxy evolution and help to anchor blind ALMA and NOEMA surveys attempting to determine the abundance of molecular gas at high redshifts., Comment: 8 pages, 6 figures, submitted to MNRAS

Outflows form an integral component in regulating the gas cycling in and out of galaxies, although their impact on the galaxy hosts is still poorly understood. Here we present an analysis of 405 high mass (log M$_{*}$/M$_{\odot}\geqslant10$), star-forming galaxies (excluding AGN) with low inclinations at $z\sim$0, using stacking techniques of the NaD $\lambda\lambda$5889,5895 A neutral gas tracer in IFU observations from the MaNGA DR15 survey. We detect outflows in the central regions of 78/405 galaxies and determine their extent and power through the construction of stacked annuli. We find outflows are most powerful in central regions and extend out to $\sim$1R$_{e}$, with declining mass outflow rates and loading factors as a function of radius. The stacking of spaxels over key galaxy quantities reveals outflow detections in regions of high $\Sigma_{\text{SFR}}$ ($\gtrsim$0.01 M$_{\odot}$yr$^{-1}$kpc$^{-2}$) and $\Sigma_{M_{*}}$ ($\gtrsim$10$^{7}$ M$_{\odot}$kpc$^{-2}$) along the resolved main sequence. Clear correlations with $\Sigma_{\text{SFR}}$ suggest it is the main regulator of outflows, with a critical threshold of $\sim$0.01 M$_{\odot}$yr$^{-1}$kpc$^{-2}$ needed to escape the weight of the disk and launch them. Furthermore, measurements of the H$\delta$ and D$_{n}$4000 indices reveal virtually identical star formation histories between galaxies with outflows and those without. Finally, through stacking of HI 21 cm observations for a subset of our sample, we find outflow galaxies show reduced HI gas fractions at central velocities compared to their non-detection control counterparts, suggestive of some removal of HI gas, likely in the central regions of the galaxies, but not enough to completely quench the host., Comment: 20 pages, 11 figures. Accepted for publication in MNRAS

We present a comparison of galaxy atomic and molecular gas properties in three recent cosmological hydrodynamic simulations, Simba, EAGLE, and Illustris-TNG, versus observations from $z\sim 0-2$. These simulations all rely on similar sub-resolution prescriptions to model cold interstellar gas which they cannot represent directly, and qualitatively reproduce the observed $z\approx 0$ HI and H$_2$ mass functions (HIMF, H2MF), CO(1-0) luminosity functions (COLF), and gas scaling relations versus stellar mass, specific star formation rate, and stellar surface density $\mu_*$, with some quantitative differences. To compare to the COLF, we apply an H$_2$-to-CO conversion factor to the simulated galaxies based on their average molecular surface density and metallicity, yielding substantial variations in $\alpha_{\rm CO}$ and significant differences between models. Using this, predicted $z=0$ COLFs agree better with data than predicted H2MFs. Out to $z\sim 2$, EAGLE's and Simba's HIMF and COLF strongly increase, while TNG's HIMF declines and COLF evolves slowly. EAGLE and Simba reproduce high $L_{\rm CO1-0}$ galaxies at $z\sim 1-2$ as observed, owing partly to a median $\alpha_{\rm CO}(z=2)\sim 1$ versus $\alpha_{\rm CO}(z=0)\sim 3$. Examining \HI, H$_2$, and CO scaling relations, their trends with $M_*$ are broadly reproduced in all models, but EAGLE yields too little HI in green valley galaxies, TNG and Simba overproduce cold gas in massive galaxies, and Simba overproduces molecular gas in small systems. Using Simba variants that exclude individual AGN feedback modules, we find that Simba's AGN jet feedback is primarily responsible by lowering cold gas contents from $z\sim 1\to0$ by suppressing cold gas in $M_*> 10^{10}{\rm M}_\odot$ galaxies, while X-ray feedback suppresses the formation of high-$\mu_*$ systems., Comment: 22 pages, MNRAS accepted

We introduce a method to estimate the HI mass within the optical radius of disk galaxies from integrated HI spectra, with an uncertainty of 0.09 dex. We use these estimates to study how inner HI fuels star formation in late-type disk galaxies. We find that star formation rate (SFR) at a given stellar mass ($M_*$) is well correlated with the inner HI surface density ($\Sigma_{\rm HI,in}$) and inner HI mass-to-stellar mass ratio. For the massive ($M_*>10^{10} M_{\odot}$) disk galaxies, higher SFR at a given stellar mass is also related with higher efficiency of converting inner HI to molecular gas, but no such correlation is found for the total HI mass. The highest $\Sigma_{\rm HI,in}$ and the fastest depletion of the total neutral gas within the optical disks are found in the most compact and star-forming disk galaxies at a given stellar mass. These results highlight the important role of inner HI as an intermediate step of fueling star formation in disk galaxies., Comment: 13 pages, 12 figures, accepted for publication at ApJ

Feedback from AGN is thought to be key in shaping the life cycle of their host galaxies by regulating star-formation activity. Therefore, to understand the impact of AGN on star formation, it is essential to trace the molecular gas out of which stars form. In this paper we present the first systematic study of the CO properties of AGN hosts at z~2 for a sample of 27 X-ray selected AGN spanning two orders of magnitude in AGN bolometric luminosity (Lbol= 10^44.7-10^46.9 erg/s) by using ALMA Band 3 observations of the CO(3-2) transition (~1" angular resolution). To search for evidence of AGN feedback on the CO properties of the host galaxies, we compared our AGN with a sample of inactive (i.e., non-AGN) galaxies from the PHIBSS survey with similar redshift, stellar masses, and SFRs. We used the same CO transition as a consistent proxy for the gas mass for the two samples in order to avoid systematics involved when assuming conversion factors. By adopting a Bayesian approach to take upper limits into account, we analyzed CO luminosities as a function of stellar masses and SFRs, as well as the ratio LCO(3-2)/M* (proxy for the gas fraction). The two samples show statistically consistent trends in the LCO(3-2)-Lfir and LCO(3-2)-M* planes. However, there are indications that AGN feature lower CO(3-2) luminosities (0.4-0.7 dex) than inactive galaxies at the 2-3sigma level when we focus on the subset of parameters where the results are better constrained and on the distribution of the mean LCO(3-2)/M*. Therefore, even by conservatively assuming the same excitation factor r31, we would find lower molecular gas masses in AGN, and assuming higher r31 would exacerbate this difference. We interpret our result as a hint of the potential effect of AGN activity (e.g., radiation and outflows), which may be able to heat, excite, dissociate, and/or deplete the gas reservoir of the host galaxies. (abridged), Comment: Key figures: 2, 3, 4. 22 pages, 9 figures (5 in appendices), 3 tables. Accepted for publication in A&A

We perform a joint-analysis of high spatial resolution molecular gas and star-formation rate (SFR) maps in main-sequence star-forming galaxies experiencing galactic-scale outflows of ionised gas. Our aim is to understand the mechanism that determines which galaxies are able to launch these intense winds. We observed CO(1-0) at 1" resolution with ALMA in 16 edge-on galaxies, which also have 2" spatial resolution optical integral field observations from the SAMI Galaxy Survey. Half the galaxies in the sample were previously identified as harbouring intense and large-scale outflows of ionised gas ("outflow-types"), the rest serve as control galaxies. The dataset is complemented by integrated CO(1-0) observations from the IRAM 30-m telescope to probe the total molecular gas reservoirs. We find that the galaxies powering outflows do not possess significantly different global gas fractions or star-formation efficiencies when compared with a control sample. However, the ALMA maps reveal that the molecular gas in the outflow-type galaxies is distributed more centrally than in the control galaxies. For our outflow-type objects, molecular gas and star-formation is largely confined within their inner effective radius ($\rm r_{eff}$), whereas in the control sample the distribution is more diffuse, extending far beyond $\rm r_{eff}$. We infer that outflows in normal star-forming galaxies may be caused by dynamical mechanisms that drive molecular gas into their central regions, which can result in locally-enhanced gas surface density and star-formation., Comment: 21 pages

We present the host galaxy molecular gas properties of a sample of 213 nearby (0.0110^44 erg/s) increases by ~10-100 between a molecular gas mass of 10^8.7 Msun and 10^10.2 Msun. Higher Eddington ratio AGN galaxies tend to have higher molecular gas masses and gas fractions. Higher column density AGN galaxies (Log NH>23.4) are associated with lower depletion timescales and may prefer hosts with more gas centrally concentrated in the bulge that may be more prone to quenching than galaxy wide molecular gas. The significant average link of host galaxy molecular gas supply to SMBH growth may naturally lead to the general correlations found between SMBHs and their host galaxies, such as the correlations between SMBH mass and bulge properties and the redshift evolution of star formation and SMBH growth., Comment: 53 pages, 37 figures, accepted in ApJ