Spatially-resolved emission line kinematics are invaluable to investigating fundamental galaxy properties and have become increasingly accessible for galaxies at $z\gtrsim0.5$ through sensitive near-infrared imaging spectroscopy and millimeter interferometry. Kinematic modeling is at the core of the analysis and interpretation of such data sets, which at high-z present challenges due to lower signal-to-noise ratio (S/N) and resolution compared to data of local galaxies. We present and test the 3D fitting functionality of DysmalPy, examining how well it recovers intrinsic disk rotation velocity and velocity dispersion, using a large suite of axisymmetric models, covering a range of galaxy properties and observational parameters typical of $z\sim1$-$3$ star-forming galaxies. We also compare DysmalPy's recovery performance to that of two other commonly used codes, GalPak3D and 3DBarolo, which we use in turn to create additional sets of models to benchmark DysmalPy. Over the ranges of S/N, resolution, mass, and velocity dispersion explored, the rotation velocity is accurately recovered by all tools. The velocity dispersion is recovered well at high S/N, but the impact of methodology differences is more apparent. In particular, template differences for parametric tools and S/N sensitivity for the non-parametric tool can lead to differences up to a factor of 2. Our tests highlight and the importance of deep, high-resolution data and the need for careful consideration of: (1) the choice of priors (parametric approaches), (2) the masking (all approaches) and, more generally, evaluating the suitability of each approach to the specific data at hand. This paper accompanies the public release of DysmalPy.
We present an analysis of millimeter CO observations to search and quantify signatures of molecular gas outflows. We exploit the large sample of $0.5 < z < 2.6$ galaxies observed as part of the PHIBSS1/2 surveys with the IRAM Plateau de Bure interferometer, focusing on the 154 typical massive star-forming galaxies with CO detections (mainly CO(3-2), but including also CO(2-1) and CO(6-5)) at signal-to-noise (SNR) > 1.5 and available properties (stellar mass, star formation rate, size) from ancillary data. None of the individual spectra exhibit a compelling signature of CO outflow emission even at high SNR > 7. To search for fainter outflow signatures, we carry out an analysis of stacked spectra, including the full sample, as well as subsets, split in terms of stellar mass, redshift, inclination, offset in star formation rate (SFR) from the main sequence, and AGN activity. None of the physically motivated subsamples show any outflow signature. We report a tentative detection in a subset statistically designed to maximize outflow signatures. We derive upper limits on molecular gas outflow rate and mass loading factors $\eta$ based on our results and find $\eta \leq$ 2.2-35.4, depending on the subsample. Much deeper CO data and observations of alternative tracers are needed to decisively constrain the importance of cold molecular gas component of outflows relative to other gas phases.
As we learn more about the multi-scale interstellar medium (ISM) of our Galaxy, we develop a greater understanding for the complex relationships between the large-scale diffuse gas and dust in Giant Molecular Clouds (GMCs), how it moves, how it is affected by the nearby massive stars, and which portions of those GMCs eventually collapse into star forming regions. The complex interactions of those gas, dust and stellar populations form what has come to be known as the ecology of our Galaxy. Because we are deeply embedded in the plane of our Galaxy, it takes up a significant fraction of the sky, with complex dust lanes scattered throughout the optically recognizable bands of the Milky Way. These bands become bright at (sub-)millimetre wavelengths, where we can study dust thermal emission and the chemical and kinematic signatures of the gas. To properly study such large-scale environments, requires deep, large area surveys that are not possible with current facilities. Moreover, where stars form, so too do planetary systems, growing from the dust and gas in circumstellar discs, to planets and planetesimal belts. Understanding the evolution of these belts requires deep imaging capable of studying belts around young stellar objects to Kuiper belt analogues around the nearest stars. Here we present a plan for observing the Galactic Plane and circumstellar environments to quantify the physical structure, the magnetic fields, the dynamics, chemistry, star formation, and planetary system evolution of the galaxy in which we live with AtLAST; a concept for a new, 50m single-dish sub-mm telescope with a large field of view which is the only type of facility that will allow us to observe our Galaxy deeply and widely enough to make a leap forward in our understanding of our local ecology.
We develop a new framework to constrain the source redshift. The method jointly accounts for the detection/non-detection of spectral lines and the prior information from the photometric redshift and total infrared luminosity from spectral energy distribution analysis. The method uses the estimated total infrared luminosity to predict the line fluxes at given redshifts and generates model spectra. The redshift-dependent spectral models are then compared with the observed spectra to find the redshift. Results. We apply the aforementioned joint redshift analysis method to four high-z dusty star-forming galaxy candidates selected from the NIKA2 observations of the HLSJ091828.6+514223 (HLS) field, and further observed by NOEMA with blind spectral scans. These sources only have SPIRE/Herschel photometry as ancillary data. They were selected because of very faint or no SPIRE counterparts, as to bias the sample towards the highest redshift candidates. The method finds the spectroscopic redshift of 4 in the 5 NOEMA-counterpart detected sources, with z>3. Based on these measurements, we derive the CO/[CI] lines and millimeter continuum fluxes from the NOEMA data and study their ISM and star-formation properties. We find cold dust temperatures in some of the HLS sources compared to the general population of sub-millimeter galaxies, which might be related to the bias introduced by the SPIRE-dropout selection. Our sources, but one, have short gas depletion time of a few hundred Myrs, which is typical among high-z sub-millimeter galaxies. The only exception shows a longer gas depletion time, up to a few Gyrs, comparable to that of main-sequence galaxies at the same redshift. Furthermore, we identify a possible over-density of dusty star-forming galaxies at z=5.2, traced by two sources in our sample, as well as the lensed galaxy HLSJ091828.6+514223. (abridged)
We aim to understand the physical mechanisms that drive star formation in a sample of mass-complete (>10$^{9.5}M_{\odot}$) star-forming galaxies (SFGs) at 1.2 $\leq z$ < 1.6. We selected SFGs from the COSMOS2020 catalog and applied a $uv$-domain stacking analysis to their archival Atacama Large Millimeter/submillimeter Array (ALMA) data. Our stacking analysis provides precise measurements of the mean molecular gas mass and size of SFGs. We also applied an image-domain stacking analysis on their \textit{HST} $i$-band and UltraVISTA $J$- and $K_{\rm s}$-band images. Correcting these rest-frame optical sizes using the $R_{\rm half-stellar-light}$-to-$R_{\rm half-stellar-mass}$ conversion at rest 5,000 angstrom, we obtain the stellar mass size of MS galaxies. Across the MS (-0.2 < $\Delta$MS < 0.2), the mean molecular gas fraction of SFGs increases by a factor of $\sim$1.4, while their mean molecular gas depletion time decreases by a factor of $\sim$1.8. The scatter of the MS could thus be caused by variations in both the star formation efficiency and molecular gas fraction of SFGs. The majority of the SFGs lying on the MS have $R_{\rm FIR}$ $\approx$ $R_{\rm stellar}$. Their central regions are subject to large dust attenuation. Starbursts (SBs, $\Delta$MS>0.7) have a mean molecular gas fraction $\sim$2.1 times larger and mean molecular gas depletion time $\sim$3.3 times shorter than MS galaxies. Additionally, they have more compact star-forming regions ($\sim$2.5~kpc for MS galaxies vs. $\sim$1.4~kpc for SBs) and systematically disturbed rest-frame optical morphologies, which is consistent with their association with major-mergers. SBs and MS galaxies follow the same relation between their molecular gas mass and star formation rate surface densities with a slope of $\sim1.1-1.2$, that is, the so-called KS relation.
In the framework of a systematic ALMA study of IR-selected main-sequence and starburst galaxies at z~1-1.7 at typical ~1" resolution, we report on the effects of mid-IR- and X-ray-detected active galactic nuclei (AGN) on the reservoirs and excitation of molecular gas in a sample of 55 objects. We find detectable nuclear activity in ~30% of the sample. The presence of dusty tori influences the IR SED of galaxies, as highlighted by the strong correlation among the AGN contribution to the total IR luminosity budget (fAGN = LIR,AGN/LIR), its hard X-ray emission, and the Rayleigh-Jeans to mid-IR (S1.2mm/S24um) observed color, with consequences on the empirical SFR estimates. Nevertheless, we find only marginal effects of AGN on the CO (J=2,4,5,7) or neutral carbon ([CI](1-0), [CI](2-1)) line luminosities and on the derived molecular gas excitation as gauged by line ratios and the full SLEDs. The [CI] and CO emission up to J=5,7 thus primarily traces the properties of the host in typical IR luminous galaxies. However, we highlight the existence of a large variety of line luminosities and ratios despite the homogeneous selection. In particular, we find a sparse group of AGN-dominated sources with the highest LIR,AGN/LIR,SFR ratios, >3, that are more luminous in CO(5-4) than what is predicted by the L'CO(5-4)-LIR,SFR relation, which might be the result of the nuclear activity. For the general population, our findings translate into AGN having minimal effects on quantities such as gas and dust fractions and SFEs. If anything, we find hints of a marginal tendency of AGN hosts to be compact at far-IR wavelengths and to display 1.8x larger dust optical depths. In general, this is consistent with a marginal impact of the nuclear activity on the gas reservoirs and star formation in average star-forming AGN hosts with LIR>5e11 Lsun, typically underrepresented in surveys of quasars and SMGs.
Abstract We present 5.5 GHz observations with the Very Large Array of a sample of nearby galaxies with energetic nuclear outbursts at mid-infrared (MIR) bands. These observations reach a uniform depth down to a median rms of ∼10 μ Jy, representing one of the most sensitive searches for radio emission associated with nuclear transients. We detect radio emission in 12 out of 16 galaxies at a level of >5 σ , corresponding to a detection rate of 75%. Such a high detection is remarkably different from previous similar searches in stellar tidal disruption events. The radio emission is compact and not resolved for the majority of sources on scales of ≲0.″5 (<0.9 kpc at z < 0.1). We find that the possibility of the star formation contributing to the radio emission is low, but an active galactic nucleus (AGN) origin remains a plausible scenario, especially for sources that show evidence of AGN activity in their optical spectra. If the detections could represent radio emission associated with a nuclear transient phenomenon such as a jet or outflow, we could use the blast wave model by analogy with the gamma-ray burst afterglows to describe the evolution of radio light curves. In this context, the observations are consistent with a decelerating jet with an energy of ∼10 51–52 erg viewed at 30°–60° off-axis at later times, suggesting that powerful jets may be ubiquitous among MIR-burst galaxies. Future continuous monitoring observations will be crucial to decipher the origin of radio emission through detections of potential flux and spectral evolution. Our results highlight the importance of radio observations to constrain the nature of nuclear MIR outbursts in galaxies.
Our knowledge of relations between supermassive black holes and their host galaxies at $z\gtrsim1$ is still limited, even though being actively sought out to $z\sim6$. Here, we use the high resolution and sensitivity of JWST to measure the host galaxy properties for 61 X-ray-selected type-I AGNs at $0.7
The NOrthern Extended Millimeter Array (NOEMA) formIng Cluster survEy (NICE) is a NOEMA large programme targeting 69 massive galaxy group candidates at z > 2 over six deep fields with a total area of 46 deg 2 . Here we report the spectroscopic confirmation of eight massive galaxy groups at redshifts 1.65 ≤ z ≤ 3.61 in the Cosmic Evolution Survey (COSMOS) field. Homogeneously selected as significant overdensities of red IRAC sources that have red Herschel colours, four groups in this sample are confirmed by CO and [CI] line detections of multiple sources with NOEMA 3 mm observations, three are confirmed with Atacama Large Millimeter Array (ALMA) observations, and one is confirmed by H α emission from Subaru/FMOS spectroscopy. Using rich ancillary data in the far-infrared and sub-millimetre, we constructed the integrated far-infrared spectral energy distributions for the eight groups, obtaining a total infrared star formation rate (SFR) of 260–1300 M ⊙ yr −1 . We adopted six methods for estimating the dark matter masses of the eight groups, including stellar mass to halo mass relations, overdensity with galaxy bias, and NFW profile fitting to radial stellar mass densities. We find that the radial stellar mass densities of the eight groups are consistent with a NFW profile, supporting the idea that they are collapsed structures hosted by a single dark matter halo. The best halo mass estimates are log( M h /M ⊙ ) = 12.8 − 13.7 with a general uncertainty of 0.3 dex. Based on the halo mass estimates, we derived baryonic accretion rates (BARs) of (1 − 8)×10 3 M ⊙ /yr for this sample. Together with massive groups in the literature, we find a quasi-linear correlation between the integrated SFR/BAR ratio and the theoretical halo mass limit for cold streams, M stream / M h , with SFR/BAR = 10 −0.46 ± 0.22 ( M stream / M h ) 0.71 ± 0.16 with a scatter of 0.40 dex. Furthermore, we compared the halo masses and the stellar masses with simulations, and find that the halo masses of all structures are consistent with those of progenitors of M h ( z = 0) > 10 14 M ⊙ galaxy clusters, and that the most massive central galaxies have stellar masses consistent with those of the brightest cluster galaxy progenitors in the TNG300 simulation. Above all, the results strongly suggest that these massive structures are in the process of forming massive galaxy clusters via baryonic and dark matter accretion.
We present new CO ($J=5-4$ and $7-6$) and [CI] ($^3P_2\,-\, ^3P_1$ and $^3P_1\,-\, ^3P_0$) emission line observations of the star-forming galaxy D49 at the massive end of the Main Sequence at $z=3$. We incorporate previous CO ($J=3-2$) and optical-to-millimetre continuum observations to fit its spectral energy distribution (SED). Our results hint at high-$J$ CO luminosities exceeding the expected location on the empirical correlations with the infrared luminosity. [CI] emission fully consistent with the literature trends is found. We do not retrieve any signatures of a bright active galactic nucleus that could boost the $J=5-4,\,7-6$ lines in either the infrared or X-ray bands, but warm photon-dominated regions, shocks or turbulence could in principle do so. We suggest that mechanical heating could be a favourable mechanism able to enhance the gas emission at fixed infrared luminosity in D49 and other main-sequence star-forming galaxies at high redshift, but further investigation is necessary to confirm this explanation. We derive molecular gas masses from dust, CO, and [CI] that all agree within the uncertainties. Given its large star formation rate (SFR) $\sim 500~M_\odot~{\rm yr}^{-1}$ and stellar mass $>10^{11.5}~M_\odot$, the short depletion time scale of $<0.3$ Gyr might indicate that D49 is experiencing its last growth spurt and will soon transit to quiescence.
'/%3e%3cuse xlink:href='%23a' transform='rotate(23.036 .093 25.536)'/%3e%3cuse xlink:href='%23a' transform='rotate(45.87 1.273 16.18)'/%3e%3cuse xlink:href='%23a' transform='rotate(69.945 .996 12.078)'/%3e%3cuse xlink:href='%23a' transform='rotate(20.66 -19.689 31.932)'/%3e%3c/svg%3e)
