Search

The NOEMA formIng Cluster survEy (NICE) is a large program targeting 69 massive galaxy group candidates at $z>2$ in six deep fields. We report spectroscopic confirmation of eight groups at $1.65\leq z\leq3.61$ in COSMOS. Homogeneously selected as significant overdensities of red IRAC sources with red Herschel colors, four groups are confirmed by CO and [CI] with NOEMA 3mm observations, three are confirmed with ALMA, and one is confirmed by H$\alpha$ from Subaru/FMOS. We constructed the integrated FIR SEDs for the eight groups, obtaining total IR SFR $=260-1300~{\rm M_\odot}$~yr$^{-1}$. We adopted six methods to estimate the dark matter masses, including stellar mass to halo mass relations, overdensity with galaxy bias, and NFW profile fitting to radial stellar mass density. We found the radial stellar mass density are consistent with a NFW profile, supporting that they are collapsed structures hosted by a single dark matter halo. The best halo mass estimates are $\log(M_{\rm h}/{\rm M_\odot})=12.8-13.7$ with uncertainty of 0.3 dex. From halo mass estimates, we derive baryonic accretion rate ${\rm BAR}=(1-8)\times10^{3}\,{\rm M_{\odot}/yr}$ for this sample. We find a quasi-linear correlation between the integrated SFR/BAR and the theoretical halo mass limit for cold streams, $M_{\rm stream}/M_{\rm h}$, with ${\rm SFR/BAR}=10^{-0.46\pm0.22}\left({M_{\rm stream}/M_{\rm h}}\right)^{0.71\pm0.16}$ with a scatter of $0.40\,{\rm dex}$. Further, we compare halo masses and stellar masses with simulations, and find all structures are consistent with being progenitors of $M_{\rm h}(z=0)>10^{14}\,{\rm M_{\odot}}$ galaxy clusters, and the most massive central galaxies have stellar masses consistent with brightest cluster galaxies (BCGs) progenitors in the TNG300 simulation. The results strongly suggest these structures are forming massive galaxy clusters via baryonic and dark matter accretion., Comment: 44 pages (27pp appendix), 32 figures, 18 tables, accepted for publication in A&A

The study of distant galaxy groups and clusters at the peak epoch of star formation is limited by the lack of a statistically and homogeneously selected and spectroscopically confirmed sample. Recent discoveries of concentrated starburst activities in cluster cores have opened a new window to hunt for these structures based on their integrated IR luminosities. Hereby we carry out the large NOEMA (NOrthern Extended Millimeter Array) program targeting a statistical sample of infrared-luminous sources associated with overdensities of massive galaxies at z>2, the Noema formIng Cluster survEy (NICE). We present the first result from the ongoing NICE survey, a compact group at z=3.95 in the Lockman Hole field (LH-SBC3), confirmed via four massive (M_star>10^10.5M_sun) galaxies detected in CO(4-3) and [CI](1-0) lines. The four CO-detected members of LH-SBC3 are distributed over a 180 kpc physical scale, and the entire structure has an estimated halo mass of ~10^13Msun and total star formation rate (SFR) of ~4000Msun/yr. In addition, the most massive galaxy hosts a radio-loud AGN with L_1.4GHz, rest = 3.0*10^25W/Hz. The discovery of LH-SBC3 demonstrates the feasibility of our method to efficiently identify high-z compact groups or forming cluster cores. The existence of these starbursting cluster cores up to z~4 provides critical insights into the mass assembly history of the central massive galaxies in clusters., Comment: 10 pages, 8 figures, published by A&A

Molecular gas plays a critical role in explaining the quiescence of star formation (SF) in massive isolated spiral galaxies, which could be a result of either the low molecular gas content and/or the low SF efficiency. We present IRAM 30m observations of the CO lines in the Sombrero galaxy (NGC~4594), the most massive spiral at $d\lesssim30\rm~Mpc$. We detect at least one of the three CO lines covered by our observations in all 13 observed positions located at the galactic nucleus and along a $\sim25\rm~kpc$-diameter dusty ring. The total extrapolated molecular gas mass of the galaxy is $M_{\rm H_2}\approx4\times10^{8}\rm~M_\odot$. The measured maximum CO gas rotation velocity of $\approx379\rm~km~s^{-1}$ suggests that NGC~4594 locates in a dark matter halo with a mass $M_{\rm200}\gtrsim10^{13}\rm~M_\odot$. Comparing to other galaxy samples, NGC~4594 is extremely gas poor and SF inactive, but the SF efficiency is apparently not inconsistent with that predicted by the Kennicutt-Schmidt law, so there is no evidence of enhanced SF quenching in this extremely massive spiral with a huge bulge. We also calculate the predicted gas supply rate from various sources to replenish the cold gas consumed in SF, and find that the galaxy must experienced a starburst stage at high redshift, then the leftover or recycled gas provides SF fuels to maintain the gradual growth of the galactic disk at a gentle rate., Comment: 21 pages, 13 figures, accepted for publication in MNRAS

Both Galactic and extragalactic studies on star formation suggest that stars form directly from dense molecular gas. To trace such high volume density gas, HCN and HCO+ J=1-0 have been widely used for their high dipole moments, relatively high abundances, and often being the strongest lines after CO. However, HCN and HCO+ J=1-0 emission could be arguably dominated by the gas components at low volume densities. HCN J=2-1 and HCO+ J=2-1, with more suitable critical densities and excitation requirements, would trace typical dense gas closely related to star formation. Here we report new observations of HCN J=2-1 and HCO+ J=2-1 towards 17 nearby infrared-bright galaxies with the APEX 12-m telescope. The correlation slopes between luminosities of HCN J=2-1, and HCO+ J=2-1 and total infrared emission are 1.03 +- 0.05 and 1.00 +- 0.05, respectively. The correlations of their surface densities, normalised with the area of radio/sub-millimeter continuum, show even tighter relations (Slopes: 0.99 +- 0.03 and 1.02 +- 0.03). The eight AGN-dominated galaxies show no significant difference from the eleven star-formation dominated galaxies in above relations. The average HCN/HCO+ ratios are 1.15 +- 0.26 and 0.98 +- 0.42 for AGN-dominated and star-formation dominated galaxies, respectively, without obvious dependencies on infrared luminosity, dust temperature, or infrared pumping. The Magellanic Clouds roughly follow the same correlations, expanding to eight orders of magnitude. On the other hand, ultra-luminous infrared galaxies with active galactic nucleus (AGN) systematically lay above the correlations, indicating potential biases introduced by AGNs., Comment: 28 pages, 12 figures, 6 tables

Using NOEMA and ALMA 3mm line scans, we measure spectroscopic redshifts of six new dusty galaxies at 3.56, based on template IR spectrum energy distribution (SED) from known submillimeter galaxies at z=4--6. Dust SED analyses explain the photo-z overestimate from seemingly cold dust temperatures (Td) and steep Rayleigh-Jeans (RJ) slopes, providing additional examples of cold dusty galaxies impacted by the Cosmic Microwave Background (CMB). We therefore study the general properties of the enlarged sample of 10 ``cold" dusty galaxies over 3.5~4, together with the severe impact of the CMB on their RJ observables, paving the way for the diagnostics of optically thick dust in the early universe. Conventional gas mass estimates based on RJ dust continuum luminosities implicitly assume an optically thin case, overestimating gas masses by a factor of 2--3 on average in compact dusty star-forming galaxies., Comment: Published version in A&A

We present well-resolved near-IR and sub-mm analysis of the three highly star-forming massive ($>10^{11}\,\rm M_{\odot}$) galaxies within the core of the RO-1001 galaxy group at $\rm z=2.91$. Each of them displays kpc-scale compact star-bursting cores with properties consistent with forming galaxy bulges, embedded at the center of extended, massive stellar disks. Surprisingly, the stellar disks are unambiguously both quiescent, and severely lopsided. Therefore, `outside-in' quenching is ongoing in the three group galaxies. We propose an overall scenario in which the strong mass lopsidedness in the disks (ranging from factors of 1.6 to $>$3), likely generated under the effects of accreted gas and clumps, is responsible for their star-formation suppression, while funnelling gas into the nuclei and thus creating the central starbursts. The lopsided side of the disks marks the location of accretion streams impact, with additional matter components (dust and stars) detected in their close proximity directly tracing the inflow direction. The interaction with the accreted clumps, which can be regarded as minor-mergers, leads the major axes of the three galaxies to be closely aligned with the outer Lyman-$\alpha$-emitting feeding filaments. These results provide the first observational evidence of the impact of cold accretion streams on the formation and evolution of the galaxies they feed. In the current phase, this is taking the form of the rapid buildup of bulges under the effects of accretion, while still preserving massive quiescent and lopsided stellar disks at least until encountering a violent major-merger., Comment: Accepted for publication in A&A

Submillimetre galaxies represent a rapid growth phase of both star formation and massive galaxies. Mapping SMGs in galaxy protoclusters provides key insights into where and how these extreme starbursts take place in connections with the assembly of the large-scale structure in the early Universe. We search for SMGs at 850$\,\mu m$ using JCMT/SCUBA-2 in two massive protoclusters at $z=2.24$, BOSS1244 and BOSS1542, and detect 43 and 54 sources with $S_{850}>4\,$mJy at the $4\sigma$ level within an effective area of 264$\,$arcmin$^2$, respectively. We construct the intrinsic number counts and find that the abundance of SMGs is $2.0\pm0.3$ and $2.1\pm0.2$ times that of the general fields, confirming that BOSS1244 and BOSS1542 contain a higher fraction of dusty galaxies with strongly enhanced star formation. The volume densities of the SMGs are estimated to be $\sim15-$30 times the average, significantly higher than the overdensity factor ($\sim 6$) traced by H$\alpha$ emission-line galaxies (HAEs). More importantly, we discover a prominent offset between the spatial distributions of the two populations in these two protoclusters -- SMGs are mostly located around the high-density regions of HAEs, and few are seen inside these regions. This finding may have revealed for the first time the occurrence of violent star formation enhancement in the outskirts of the HAE density peaks, likely driven by the boosting of gas supplies and/or starburst triggering events. Meanwhile, the lack of SMGs inside the most overdense regions at $z\sim2$ implies a transition to the environment disfavouring extreme starbursts., Comment: 16pages, 11figures, accepted for publication in MNRAS

While the slope of the dust attenuation curve ($\delta$) is found to correlate with effective dust attenuation ($A_V$) as obtained through spectral energy distribution (SED) fitting, it remains unknown how the fitting degeneracies shape this relation. We examine the degeneracy effects by fitting SEDs of a sample of local star-forming galaxies (SFGs) selected from the Galaxy And Mass Assembly survey, in conjunction with mock galaxy SEDs of known attenuation parameters. A well-designed declining starburst star formation history is adopted to generate model SED templates with intrinsic UV slope ($\beta_0$) spanning over a reasonably wide range. The best-fitting $\beta_0$ for our sample SFGs shows a wide coverage, dramatically differing from the limited range of $\beta_0<-2.2$ for a starburst of constant star formation. Our results show that strong degeneracies between $\beta_0$, $\delta$, and $A_V$ in the SED fitting induce systematic biases leading to a false $A_V$--$\delta$ correlation. Our simulation tests reveal that this relationship can be well reproduced even when a flat $A_V$--$\delta$ relation is taken to build the input model galaxy SEDs. The variations in best-fitting $\delta$ are dominated by the fitting errors. We show that assuming a starburst with constant star formation in SED fitting will result in a steeper attenuation curve, smaller degeneracy errors, and a stronger $A_V$--$\delta$ relation. Our findings confirm that the $A_V$--$\delta$ relation obtained through SED fitting is likely driven by the systematic biases induced by the fitting degeneracies between $\beta_0$, $\delta$, and $A_V$., Comment: 21 pages, 13 figures, accepted for publication in the MNRAS, Comments welcome!

We present a new study of the merger dynamics of Abell~1775 by analyzing the high-quality Chandra and XMM-Newton archival data. We confirm/identify an arc-shaped edge (i.e., the head) at $\sim48$~kpc west of the X-ray peak, a split cold gas tail that extends eastward to $\sim163$~kpc, and a plume of spiral-like X-ray excess (within about $81-324$~kpc northeast of the cluster core) that connects to the end of the tail. The head, across which the projected gas temperature rises outward from $3.39_{-0.18}^{+0.28}$~keV to $5.30_{-0.43}^{+0.54}$~keV, is found to be a cold front with a Mach number of $\mathcal{M}\sim0.79$. Along the surfaces of the cold front and tail, typical KHI features (noses and wings, etc.) are found and are used to constrain the upper limit of the magnetic field ($\sim11.2~\mu$G) and the viscosity suppression factor ($\sim0.01$). Combining optical and radio evidence we propose a two-body merger (instead of systematic motion in a large-scale gas environment) scenario and have carried out idealized hydrodynamic simulations to verify it. We find that the observed X-ray emission and temperature distributions can be best reproduced with a merger mass ratio of 5 after the first pericentric passage. The NAT radio galaxy is thus more likely to be a single galaxy falling into the cluster center at a relative velocity of 2800~$\rm km~s^{-1}$, a speed constrained by its radio morphology. The infalling subcluster is expected to have a relatively low gas content, because only a gas-poor subcluster can cause central-only disturbances as observed in such an off-axis merger., Comment: 12 figures, 5 tables, accepted for publication in the ApJ

To study the high-transition dense-gas tracers and their relationships to the star formation of the inner $\sim$ 2 kpc circumnuclear region of NGC253, we present HCN $J=4-3$ and HCO$^+ J=4-3$ maps obtained with the James Clerk Maxwell Telescope (JCMT). With the spatially resolved data, we compute the concentration indices $r_{90}/r_{50}$ for the different tracers. HCN and HCO$^+$ 4-3 emission features tend to be centrally concentrated, which is in contrast to the shallower distribution of CO 1-0 and the stellar component. The dense-gas fraction ($f_\text{dense}$, traced by the velocity-integrated-intensity ratios of HCN/CO and HCO$^+$/CO) and the ratio $R_\text{31}$ (CO 3-2/1-0) decline towards larger galactocentric distances, but increase with higher SFR surface density. The radial variation and the large scatter of $f_\text{dense}$ and $R_\text{31}$ imply distinct physical conditions in different regions of the galactic disc. The relationships of $f_\text{dense}$ versus $\Sigma_\text{stellar}$, and SFE$_\text{dense}$ versus $\Sigma_\text{stellar}$ are explored. SFE$_\text{dense}$ increases with higher $\Sigma_\text{stellar}$ in this galaxy, which is inconsistent with previous work that used HCN 1-0 data. This implies that existing stellar components might have different effects on the high-$J$ HCN and HCO$^+$ than their low-$J$ emission. We also find that SFE$_\text{dense}$ seems to be decreasing with higher $f_\text{dense}$, which is consistent with previous works, and it suggests that the ability of the dense gas to form stars diminishes when the average density of the gas increases. This is expected in a scenario where only the regions with high-density contrast collapse and form stars., Comment: accepted to MNRAS

We present maps of atomic carbon [CI](1-0) and [CI](2-1) at a linear resolution ~1kpc scale for a sample of one HII, six LINER, three Seyfert and five starburst galaxies observed with Herschel. We compare spatial distributions of two [CI] lines with that of CO(1-0) emission, and find that both [CI] lines distribute similarly to CO(1-0) emission in most galaxies. We present luminosity ratio maps of L'_[CI](1-0)/L'_CO(1-0), L'_[CI](2-1)/L'_CO(1-0), L'_[CI](2-1)/L'_[CI](1-0) (hereafter R_[CI]) and f_70/f_160. L'_[CI](2-1)/L'_CO(1-0), R_[CI] and f_70/f_160 are centrally peaked in starbursts; whereas remain relatively constant in LINERs, indicating that star-forming activity can enhance carbon emission, especially for [CI](2-1). We explore the correlations between the luminosities of CO(1-0) and [CI] lines, and find that L'_CO(1-0) correlates tightly and almost linearly with both L'_[CI](1-0) and L'_[CI](2-1), suggesting that [CI] lines, similar as CO(1-0), can trace total molecular gas in our resolved galaxies on kpc scales. We investigate the dependence of L'_[CI](1-0)/L'_CO(1-0), L'_[CI](2-1)/L'_CO(1-0) and [CI] excitation temperature T_ex on dust temperature T_dust, and find non-correlation, a weak and modest correlation, respectively. The ratio of L'_[CI](1-0)/L'_CO(1-0) stays smooth distribution in most galaxies, indicating that the conversion factor of [CI](1-0) luminosity to H_2 mass (X_[CI](1-0)) changes with CO(1-0) conversion factor (\alpha_CO) proportionally. Under optically thin and LTE assumptions, we derive a galaxy-wide average carbon excitation temperature T_ex ~ 19.7 \pm 0.5K and an average neutral carbon abundance X[CI]/X[H_2] ~2.5 \pm 1.0 * 10^{-5} in our resolved sample, which is comparable to the usually adopted value of 3*10^{-5}, but ~3 times lower than the carbon abundance in local (U)LIRGs. We conclude that the carbon abundance varies in different galaxy types., Comment: 15 pages, 5 figures, accepted for publication in ApJ

We report the discovery of a remarkable concentration of massive galaxies with extended X-ray emission at $z_{spec} = 2.506$, which contains 11 massive ($M_{*} \gtrsim 10^{11} M_{\odot}$) galaxies in the central 80kpc region (11.6$\sigma$ overdensity). We have spectroscopically confirmed 17 member galaxies with 11 from CO and the remaining ones from $H\alpha$. The X-ray luminosity, stellar mass content and velocity dispersion all point to a collapsed, cluster-sized dark matter halo with mass $M_{200c} = 10^{13.9\pm0.2} M_{\odot}$, making it the most distant X-ray-detected cluster known to date. Unlike other clusters discovered so far, this structure is dominated by star-forming galaxies (SFGs) in the core with only 2 out of the 11 massive galaxies classified as quiescent. The star formation rate (SFR) in the 80kpc core reaches $\sim$3400 $M_{\odot}$ yr$^{-1}$ with a gas depletion time of $\sim 200$ Myr, suggesting that we caught this cluster in rapid build-up of a dense core. The high SFR is driven by both a high abundance of SFGs and a higher starburst fraction ($\sim25\%$, compared to 3\%-5\% in the field). The presence of both a collapsed, cluster-sized halo and a predominant population of massive SFGs suggests that this structure could represent an important transition phase between protoclusters and mature clusters. It provides evidence that the main phase of massive galaxy passivization will take place after galaxies accrete onto the cluster, providing new insights into massive cluster formation at early epochs. The large integrated stellar mass at such high redshift challenges our understanding of massive cluster formation., Comment: matched to the published version by ApJ

We measure the dust and gas content of the three sub-millimeter galaxies (SMGs) in the GN20 proto-cluster at $z=4.05$ using new IRAM Plateau de Bure interferometer (PdBI) CO(4-3) and 1.2-3.3 mm continuum observations. All these three SMGs are heavily dust obscured, with UV-based star formation rate (SFR) estimates significantly smaller than the ones derived from the bolometric infrared (IR), consistent with the spatial offsets revealed by HST and CO imaging. Based also on evaluations of the specific SFR, CO-to-H$_2$ conversion factor and gas depletion timescale, we classify all the three galaxies as starbursts (SBs), although with a lower confidence for GN20.2b that might be a later stage merging event. We place our measurements in the context of the evolutionary properties of main sequence (MS) and SB galaxies. ULIRGs have 3-5 times larger $L'_{\rm CO}/M_{\rm dust}$ and $M_{\rm dust}/M_\star$ ratios than $z=0$ MS galaxies, but by $z\sim2$ the difference appears to be blurred, probably due to differential metallicity evolution. SB galaxies appear to slowly evolve in their $L'_{\rm CO}/M_{\rm dust}$ and $M_{\rm dust}/M_\star$ ratios all the way to $z>6$ (consistent with rapid enrichment of SB events), while MS galaxies rapidly increase in $M_{\rm dust}/M_\star$ from $z=0$ to 2 (due to gas fraction increase, compensated by a decrease of metallicities). While no IR/submm continuum detection is available for indisputably normal massive galaxies at $z>2.5$, we show that if metallicity indeed decrease rapidly for these systems at $z>3$ as claimed in the literature, we should expect a strong decrease of their $M_{\rm dust}/M_\star$, consistent with recent PdBI and ALMA upper limits. We conclude that the $M_{\rm dust}/M_\star$ ratio could be a powerful tool for distinguishing starbursts from normal galaxies at $z>4$., Comment: 17 pages, 9 figures. Updated to A&A accepted version

We present deep IRAM Plateau de Bure Interferometer (PdBI) observations, searching for CO-emission toward two massive, non-lensed Lyman break galaxies (LBGs) at z=3.216 and 4.058. With one low significance CO detection (3.5 sigma) and one sensitive upper limit, we find that the CO lines are >~ 3-4 times weaker than expected based on the relation between IR and CO luminosities followed by similarly, massive galaxies at z=0-2.5. This is consistent with a scenario in which these galaxies have low metallicity, causing an increased CO-to-H_2 conversion factor, i.e., weaker CO-emission for a given molecular (H_2) mass. The required metallicities at z>3 are lower than predicted by the fundamental metallicity relation (FMR) at these redshifts, consistent with independent evidence. Unless our galaxies are atypical in this respect, detecting molecular gas in normal galaxies at z>3 may thus remain challenging even with ALMA., Comment: 9 pages, 3 figures, accepted for publication in Astrophysical Journal Letters

We present simultaneous observations of CO,13CO and C18O J=1-0 emission in 11 nearby (cz<1000 km/s) bright infrared galaxies. Both 12CO and 13CO are detected in the centers of all galaxies, except for 13CO in NGC 3031. We have also detected C18O, CS J=2-1, and HCO+ J=1-0 emission in the nuclear regions of M82 and M51. These are the first systematical extragalactic detections of 12CO and its isotopes from the PMO 14m telescope. We have conducted half-beam spacing mapping of M82 over an area of 4'*2.5' and major axis mapping of NGC 3627, NGC 3628, NGC 4631, and M51. The radial distributions of 12CO and 13CO in NGC 3627, NGC 3628, and M51 can be well fitted by an exponential profile. The 12CO/13CO intensity ratio,R,decreases monotonically with galactocentric radius in all mapped sources. The average R in the center and disk of the galaxies are 9.9+/-3.0 and 5.6+/-1.9 respectively, much lower than the peculiar R(~24) found in the center of M82. The intensity ratios of 13CO/C18O, 13CO/HCO+ and 13CO/CS (either ours or literature data) show little variations with galactocentric radius, in sharp contrast with the greatly varied R. This supports the notion that the observed gradient in R could be the results of the variations of the physical conditions across the disks. The H_2 column density derived from C18O shows that the Galactic standard conversion factor (X-factor) overestimates the amount of the molecular gas in M82 by a factor of ~2.5. These observations suggest that the X-factor in active star-forming regions (i.e., nuclear regions) should be lower than that in normal star-forming disks, and the gradient in R can be used to trace the variations of the X-factor., Comment: 27 pages, 7 figures, accepted by RAA