Search

During the last decade, studies about highly attenuated and massive red star-forming galaxies (RedSFGs) at $z \sim 4$ have suggested that they could constitute a crucial population for unraveling the mechanisms driving the transition from vigorous star formation to quiescence at high redshifts. Since such a transition seems to be linked to a morphological transformation, studying the morphological properties of these RedSFGs is essential to our understanding of galaxy evolution. To this end, we are using JWST/NIRCam images from the CEERS survey to assemble a mass-complete sample of 188 massive galaxies at $z=3-4$, for which we perform resolved-SED fit. After classifying galaxies into typical blue SFGs (BlueSFGs), RedSFGs and quiescent galaxies (QGs), we compare the morphologies of each population in terms of stellar mass density, SFR density, sSFR, dust-attenuation and mass-weighted age. We find that RedSFGs and QGs present similar stellar surface density profiles and that RedSFGs manifest a dust attenuation concentration significantly higher than that of BlueSFGs at all masses. This indicates that to become quiescent, a BlueSFG must transit through a major compaction phase once it has become sufficiently massive. At the same time, we find RedSFGs and QGs to account for more than $50\%$ of galaxies with ${\rm log}(M_\ast/M_\odot)> 10.4$ at this redshift. This transition mass corresponds to the "critical mass" delineating the bimodality between BlueSFGs and QGs in the local Universe. We then conclude that there is a bimodality between extended BlueSFGs and compact, highly attenuated RedSFGs that have undergone a major gas compaction phase enabling the latter to build a massive bulb in situ. There is evidence that this early-stage separation is at the origin of the local bimodality between BlueSFGs and QGs, which we refer to as a "primeval bimodality".

The majority of stars in today's Universe reside within spheroids, which are bulges of spiral galaxies and elliptical galaxies. Their formation is still an unsolved problem. Infrared/submm-bright galaxies at high redshifts have long been suspected to be related to spheroids formation. Proving this connection has been hampered so far by heavy dust obscuration when focusing on their stellar emission or by methodologies and limited signal-to-noise ratios when looking at submm wavelengths. Here we show that spheroids are directly generated by star formation within the cores of highly luminous starburst galaxies in the distant Universe. This follows from the ALMA submillimeter surface brightness profiles which deviate significantly from those of exponential disks, and from the skewed-high axis-ratio distribution. The majority of these galaxies are fully triaxial rather than flat disks: the ratio of the shortest to the longest of their three axes is half, on average, and increases with spatial compactness. These observations, supported by simulations, reveal a cosmologically relevant pathway for in-situ spheroid formation through starbursts likely preferentially triggered by interactions (and mergers) acting on galaxies fed by non-co-planar gas accretion streams., Comment: Nature in press

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 stellar structures of star-forming galaxies (SFGs) undergo significant size growth during their mass assembly and must pass through a compaction phase as they evolve into quiescent galaxies (QGs). To shed light on the mechanisms behind this structural evolution, we study the morphology of the star-forming components of 665 SFGs at 010^{10.5}M_\odot$). The compact star-forming components of these galaxies are also concentrated and slightly spheroidal, suggesting that this compaction phase can build dense bulge in-situ. Finally, we identify a third population of SFGs (19%), with both compact stellar and star-forming components. The density of their stellar cores resemble those of QGs and are compatible with being the descendants of EC galaxy. Overall, the structural evolution of SFGs is mainly dominated by a secular inside-out growth, which can, however, be interrupted by violent compaction phase(s) that can build dominant stellar bulges like those in massive SFGs or QGs., Comment: 24 pages, 17 figures, submitted to A&A, comments are welcome

[Abridged] In recent years, conflicting results have provided an uncertain view of the dust-attenuated properties of $z>4$ star-forming galaxies (SFGs). To solve this, we used the deepest data publicly available in COSMOS to build a mass-complete ($>10^{9.5}\,M_{\odot}$) sample of SFGs at $4M_\ast)$ converges at $M_\ast<10^{9}\,M_\odot$ and is dominated by SFGs with $M_\ast\sim10^{9.5-10.5}\,M_\odot$. The fraction of the cosmic SFRD that is attenuated by dust, ${\rm SFRD}_{\rm IR}(>M_\ast)/ {\rm SFRD}(>M_\ast)$, is $90\pm4\%$ for $M_\ast\,=\,10^{10}\,M_\odot$, $68\pm10\%$ for $M_\ast=10^{8.9}\,M_\odot$ (i.e., $0.03\times M^\star$; $M^\star$ being the characteristic stellar mass of SFGs) and this value converges to $60\pm10\%$ for $M_\ast=10^{8}\,M_\odot$. Even at this early epoch, the fraction of the cosmic SFRD that is attenuated by dust remains thus significant., Comment: Accepted for publication in A&A; 14 pages and 8 figures

Age-dating and weighting stellar populations in galaxies are essential steps to study galaxy formation through cosmic times. Evolutionary population synthesis models with different input physics are used towards this aim. In particular, the contribution from the thermally pulsing asymptotic giant branch (TP-AGB) stellar phase, which peaks for intermediate-age 0.6-2 Gyr systems, has been debated upon for decades. Here we report the detection of strong cool-star signatures in the rest-frame near-infrared spectra of three young (~1 Gyr), massive (~10^10 Msun) quiescent galaxies at large look-back time, z=1-2, using JWST/NIRSpec. The co-existence of oxygen- and carbon-type absorption features, spectral edges and features from rare species such as Vanadium, and possibly Zirconium, reveal a strong contribution from TP-AGB stars. Population synthesis models with significant TP-AGB contribution reproduce the observations better than those with weak TP-AGB, which are commonly used. These findings call for revisions of published stellar population fitting results, pointing to lower masses and younger ages, with additional implications on cosmic dust production and chemical enrichment. New generations of improved models are needed, informed by these and future observations., Comment: Nature Astronomy, in press

Modern (sub)millimeter interferometers, such as ALMA and NOEMA, offer high angular resolution and unprecedented sensitivity. This provides the possibility to characterize the morphology of the gas and dust in distant galaxies. To assess the capabilities of current softwares in recovering morphologies and surface brightness profiles in interferometric observations, we test the performance of the Spergel model for fitting in the $uv$-plane, which has been recently implemented in the IRAM software GILDAS (uv$\_$fit). Spergel profiles provide an alternative to the Sersic profile, with the advantage of having an analytical Fourier transform, making them ideal to model visibilities in the $uv$-plane. We provide an approximate conversion between Spergel index and Sersic index, which depends on the ratio of the galaxy size to the angular resolution of the data. We show through extensive simulations that Spergel modeling in the $uv$-plane is a more reliable method for parameter estimation than modeling in the image-plane, as it returns parameters that are less affected by systematic biases and results in a higher effective signal-to-noise ratio (S/N). The better performance in the $uv$-plane is likely driven by the difficulty of accounting for correlated signal in interferometric images. Even in the $uv$-plane, the integrated source flux needs to be at least 50 times larger than the noise per beam to enable a reasonably good measurement of a Spergel index. We characterise the performance of Spergel model fitting in detail by showing that parameters biases are generally low (< 10%) and that uncertainties returned by uv$\_$fit are reliable within a factor of two. Finally, we showcase the power of Spergel fitting by re-examining two claims of extended halos around galaxies from the literature, showing that galaxies and halos can be successfully fitted simultaneously with a single Spergel model., Comment: 23 pages, 15 figures, accepted for publication in A&A

Understanding the gas content in galaxies, its consumption and replenishment, remains pivotal in our comprehension of the evolution of the Universe. Numerous studies have addressed this, utilizing various observational tools and analytical methods. These include examining low-transition $^{12}$CO millimeter rotational lines and exploring the far-infrared and the (sub-)millimeter emission of galaxies. With the capabilities of present-day facilities, much of this research has been centered on relatively bright galaxies. We aim at exploring the gas reservoirs of a more general type of galaxy population at $1.0\leq z\leq 3.0$. We stack ALMA 1.1 mm data to measure the gas content of a mass-complete sample down to $\sim10^{8.6}$ M$_\odot$ at $z=1$ ($\sim10^{9.2}$ M$_\odot$ at $z=3$), extracted from the HST/CANDELS sample in GOODS-S. The sample is composed of 5,530 on average blue ($\sim0.12$ mag, $\sim0.81$ mag), star-forming main sequence objects ($\Delta$MS$\sim-0.03$). We report measurements at $10^{10-11}$ M$_\odot$ and upper limits for the gas fractions at $10^{8-10}$ M$_\odot$. At $10^{10-11}$ M$_\odot$, our f$_{\mathrm{gas}}$, ranging from 0.32 to 0.48, agree well with other studies based on mass-complete samples down to $10^{10}$ M$_\odot$, and are lower than expected according to other works more biased to individual detections. At $10^{9-10}$ M$_\odot$, we obtain 3$\sigma$ upper limits for f$_{\mathrm{gas}}$ ranging from 0.69 to 0.77. These upper limits are on the level of the extrapolations of scaling relations based on mass-complete samples down to $10^{10}$ M$_\odot$. As such, it suggests that the gas content of low-mass galaxies is at most what is extrapolated from literature scaling relations. The comparison of our results with previous works reflects how the inclusion of bluer, less obscured, and more MS-like objects progressively pushes the gas level to lower values., Comment: 18 pages, 8 figures, 4 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

Recent JWST observations have revealed an unexpected abundance of massive galaxy candidates in the early Universe, extending further in redshift and to lower luminosity than what had previously been found by sub-millimeter surveys. These JWST candidates have been interpreted as challenging the $\Lambda$CDM cosmology, but, so far, they have mostly relied only on rest-frame ultraviolet data and lacked spectroscopic confirmation of their redshifts. Here we report a systematic study of 36 massive dust-obscured galaxies with spectroscopic redshifts between $z_{\rm spec}=5-9$ from the JWST FRESCO survey. We find no tension with the $\Lambda$CDM model in our sample. However, three ultra-massive galaxies (log$M_{\star}/M_{\odot}$ $\gtrsim11.0$) require an exceptional fraction of 50% of baryons converted into stars -- two to three times higher than even the most efficient galaxies at later epochs. The contribution from an active nucleus is unlikely because of their extended emission. Ultra-massive galaxies account for as much as 17% of the total cosmic star formation rate density at $z\sim5-6$., Comment: Nature in press. Updated to the accepted version. 24 pages, 4 main figures, 7 supplementary figures, 3 supplementary tables

We investigate the morphology and resolved physical properties of a sample of 22 IR-selected DSFG at cosmic noon using the JWST/NIRCam images obtained in the EGS field for the CEERS survey. The resolution of the NIRCam images allowed to spatially resolve these galaxies up to 4.4um and identify their bulge even when extinguished by dust. The goal of this study is to obtain a better understanding of the formation and evolution of FIR-bright galaxies by spatially resolving their properties using JWST in order to look through the dust and bridge the gap between the compact FIR sources and the larger optical SFG. Based on RGB images from the NIRCam filters, we divided each galaxy into several uniformly colored regions, fitted their respective SEDs, and measured physical properties. After classifying each region as SF or quiescent, we assigned galaxies to three classes, depending on whether active SF is located in the core, in the disk or in both. We find (i) that galaxies at a higher z tend to have a fragmented disk with a low core mass fraction. They are at an early stage of bulge formation. When moving toward a lower z, the core mass fraction increases, and the bulge growth is associated with a stabilization of the disk: the NIRCam data clearly point toward bulge formation in preexisting disks. (ii) Lopsidedness is a common feature of DSFGs. It could have a major impact on their evolution; (iii) 23% of galaxies have a SF core embedded in a quiescent disk. They seem to be undergoing outside-in quenching, often facilitated by their strong lopsidedness inducing instabilities. (iv) We show that half of our galaxies with SF concentrated in their core are good SMG counterpart candidates, demonstrating that compact SMGs are usually surrounded by a larger, less obscured disk. (v) Finally, we found surprising evidence for clump-like substructures being quiescent or residing in quiescent regions., Comment: 23 pages, 24 figures, Accepted by A&A

We present the results of a systematic study of the rest-frame optical morphology of quiescent galaxies at $z \geq 3$ using the Near-Infrared Camera (NIRCam) onboard $JWST$. Based on a sample selected by $UVJ$ color or $NUVUVJ$ color, we focus on 26 quiescent galaxies with $9.8<\log{(M_\star/M_\odot)}<11.4$ at $2.810.3$ implies that our size - stellar mass relations are below those at lower redshifts, with the amplitude of $\sim0.6\, {\rm kpc}$ at $M_\star = 5\times 10^{10}\, M_\odot$. This value agrees with the extrapolation from the size evolution of quiescent galaxies at $z<3$ in the literature, implying that the size of quiescent galaxies increases monotonically from $z\sim3-5$. Our sample is mainly composed of galaxies with bulge-like structures according to their median S\'ersic index and axis ratio of $n\sim3-4$ and $q\sim0.6-0.8$, respectively. On the other hand, there is a trend of increasing fraction of galaxies with low S\'ersic index, suggesting $3

Surface densities of gas, dust and stars provide a window into the physics of star-formation that, until the advent of high-resolution far-infrared/sub-millimeter observations, has been historically difficult to assess amongst dusty galaxies. To study the link between infrared (IR) surface densities and dust properties, we leverage the Atacama Large Millimetre/Submillimetre Array (ALMA) archive to measure the extent of cold dust emission in 15 $z\sim2$ IR selected galaxies selected on the basis of having available mid-IR spectroscopy from Spitzer. We use the mid-IR spectra to constrain the relative balance between dust heating from star-formation and active galactic nuclei (AGN), and to measure emission from Polycylic Aromatic Hydrocarbons (PAHs) -- small dust grains that play a key role in the photoelectric heating of gas. In general, we find that dust-obscured star-formation at high IR surface densities exhibits similar properties at low- and high-redshift, namely: local luminous IR galaxies have comparable PAH luminosity to total dust mass ratios as high-$z$ galaxies, and star-formation at $z\sim0-2$ is more efficient at high IR surface densities despite the fact that our sample of high$-z$ galaxies are closer to the main-sequence than local luminous IR galaxies. High star-formation efficiencies are coincident with a decline in the PAH/IR luminosity ratio reminiscent of the deficit observed in far-infrared fine-structure lines. Changes in the gas and dust conditions arising from high star-formation surface densities might help drive the star-formation efficiency up. This could help explain high efficiencies needed to reconcile star-formation and gas volume densities in dusty galaxies at cosmic noon., Comment: 18 pages, 10 figures, accepted to ApJ

[Abridged] We combined HST images from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey with JWST images from the Cosmic Evolution Early Release Science (CEERS) survey to measure the stellar and dust-obscured star formation distributions of a mass-complete ($>10^{10}M_\odot$) sample of 69 star-forming galaxies (SFGs) at $0.175$; 35 galaxies). At lower $S/N$, simulations show that simultaneous measurements of the size and S\'ersic index become unreliable. We extended our study to fainter sources ($S/N>10$; 69 galaxies) by fixing their S\'ersic index to unity. The S\'ersic index of bright galaxies ($S/N>75$) has a median value of 0.7, which, together with their axis ratio distribution, suggests a disk-like morphology in the rest-MIR. Galaxies above the main sequence (MS; i.e., starbursts) have rest-MIR sizes that are a factor 2 smaller than their rest-optical sizes. The median rest-optical to rest-MIR size ratio of MS galaxies increases with stellar mass, from 1.1 at $10^{9.8}M_\odot$ to 1.6 at $10^{11}M_\odot$. This mass-dependent trend resembles the one found in the literature between the rest-optical and rest-near-infrared sizes of SFGs, suggesting that it is due to radial color gradients affecting rest-optical sizes and that the sizes of the stellar and star-forming components of SFGs are, on average, consistent at all masses. There is, however, a small population of SFGs (15%) with a compact star-forming component embedded in a larger stellar structure. This could be the missing link between galaxies with an extended stellar component and those with a compact stellar component, the so-called blue nuggets., Comment: Published in A&A. 18 pages, 14 figures

We study a new population of extremely red objects (EROs) recently discovered by JWST based on their NIRCam colors F277W$-$F444W $>1.5$ mag. We find 37 EROs in the CEERS field with F444W $<28$ mag and photometric redshifts between $57$ by up to a factor $\sim$60. Similarly, if they are OSOs with luminosities in the L$_{\rm bol}>10^{46-47}$ erg s$^{-1}$ range, their number would exceed that of bright blue QSOs by more than two orders of magnitude. Additional photometry at mid-IR wavelengths will reveal the true nature of the red continuum emission in these EROs and will place this puzzling population in the right context of galaxy evolution., Comment: 26 pages, 10 figures, submitted to ApJ

We present NOEMA CO (2-1) line and ALMA 870 $\mu$m continuum observations of a main-sequence galaxy at $z=1.38$. The galaxy was initially selected as a "gas-giant", based on the gas mass derived from sub-mm continuum (log$(M_{\rm gas}/M_{\odot})=11.20\pm0.20$), however the gas mass derived from CO (2-1) luminosity brings down the gas mass to a value consistent with typical star-forming galaxies at that redshift (log$(M_{\rm gas}/M_{\odot})=10.84\pm0.03$). Despite that the dust-to-stellar mass ratio remains elevated above the scaling relations by a factor of 5. We explore the potential physical picture and consider an underestimated stellar mass and optically thick dust as possible causes. Based on the updated gas-to-stellar mass ratio we rule out the former, and while the latter can contribute to the dust mass overestimate it is still not sufficient to explain the observed physical picture. Instead, possible explanations include enhanced HI reservoirs, CO-dark H$_2$ gas, an unusually high metallicity, or the presence of an optically dark, dusty contaminant. Using the ALMA data at 870 $\mu$m coupled with $HST$/ACS imaging, we find extended morphology in dust continuum and clumpy star-formation in rest-frame UV in this galaxy, and a tentative $\sim 10$ kpc dusty arm is found bridging the galaxy center and a clump in F814W image. The galaxy shows levels of dust obscuration similar to the so-called $HST$-dark galaxies at higher redshifts, and would fall into the optically faint/dark $JWST$ color-color selection at $z>2$. It is therefore possible that our object could serve as low-$z$ analog of the $HST$-dark populations. This galaxy serves as a caveat to the gas masses based on the continuum alone, with a larger sample required to unveil the full picture., Comment: 9 pages, 5 figures, 1 table. Submitted to A&A Letters

We present a search for extremely red, dust-obscured, $z>7$ galaxies with $\textit{JWST}$/NIRCam+MIRI imaging over the first 20 arcmin$^2$ of publicly-available Cycle 1 data from the COSMOS-Web, CEERS, and PRIMER surveys. Based on their red color in F277W$-$F444W ($\sim 2.5$ mag) and detection in MIRI/F770W ($\sim 25$ mag), we identify two galaxies$\unicode{x2014}$COS-z8M1 and CEERS-z7M1$\unicode{x2014}$which have best-fit photometric redshifts of $z=8.5^{+0.3}_{-0.4}$ and $z=7.6^{+0.1}_{-0.1}$, respectively. We perform SED fitting with a variety of codes (including BAGPIPES, PROSPECTOR, BEAGLE, and CIGALE), and find a $>95\%$ probability that these indeed lie at $z>7$. Both sources are compact ($R_{\rm eff} \lesssim 200$ pc), highly obscured ($A_V \sim 1.5$$\unicode{x2013}$$2.5$), and, at our best-fit redshift estimates, likely have strong [OIII]+H$\beta$ emission contributing to their $4.4\,\mu$m photometry. We estimate stellar masses of $\sim 10^{10}~M_\odot$ for both sources; by virtue of detection in MIRI at $7.7\,\mu$m, these measurements are robust to the inclusion of bright emission lines, for example, from an AGN. We identify a marginal (2.9$\sigma$) ALMA detection at 2 mm within $0.5''$ of COS-z8M1, which if real, would suggest a remarkably high IR luminosity of $\sim 10^{12} L_\odot$. These two galaxies, if confirmed at $z\sim 8$, would be extreme in their stellar and dust masses, and may be representative of a substantial population of modestly dust-obscured galaxies at cosmic dawn., Comment: 23 pages, 8 figures, submitted to ApJL, comments welcome

In recent years, observations have uncovered a population of massive galaxies that are invisible or very faint in deep optical/near-infrared (near-IR) surveys but brighter at longer wavelengths. However, the nature of these optically dark or faint galaxies (OFGs; one of several names given to these objects) is highly uncertain. In this work, we investigate the drivers of dust attenuation in the JWST era. In particular, we study the role of stellar mass, size, and orientation in obscuring star-forming galaxies (SFGs) at $3 < z < 7.5$, focusing on the question of why OFGs and similar galaxies are so faint at optical/near-IR wavelengths. We find that stellar mass is the primary proxy for dust attenuation, among the properties studied. Effective radius and axis ratio do not show a clear link with dust attenuation, with the effect of orientation being close to random. However, there is a subset of highly dust attenuated ($A_V > 1$, typically) SFGs, of which OFGs are a specific case. For this subset, we find that the key distinctive feature is their compact size (for massive systems with $\log (M_{*}/M_{\odot}) > 10$); OFGs exhibit a 30% smaller effective radius than the average SFG at the same stellar mass and redshift. On the contrary, OFGs do not exhibit a preference for low axis ratios (i.e., edge-on disks). The results in this work show that stellar mass is the primary proxy for dust attenuation and compact stellar light profiles behind the thick dust columns obscuring typical massive SFGs., Comment: Accepted by A&A. 13 pages, 9 figures

We present the results of a systematic search for candidate quiescent galaxies in the distant Universe in eleven $JWST$ fields with publicly available observations collected during the first three months of operations and covering an effective sky area of $\sim145$ arcmin$^2$. We homogeneously reduce the new $JWST$ data and combine them with existing observations from the $Hubble\,Space\,Telescope$. We select a robust sample of $\sim80$ candidate quiescent and quenching galaxies at $3 < z < 5$ using two methods: (1) based on their rest-frame $UVJ$ colors, and (2) a novel quantitative approach based on Gaussian Mixture Modeling of the $NUV-U$, $U-V$, and $V-J$ rest-frame color space, which is more sensitive to recently quenched objects. We measure comoving number densities of massive ($M_\star\geq 10^{10.6} M_\odot$) quiescent galaxies consistent with previous estimates relying on ground-based observations, after homogenizing the results in the literature with our mass and redshift intervals. However, we find significant field-to-field variations of the number densities up to a factor of $2-3$, highlighting the effect of cosmic variance and suggesting the presence of overdensities of red quiescent galaxies at $z>3$, as it could be expected for highly clustered massive systems. Importantly, $JWST$ enables the robust identification of quenching/quiescent galaxy candidates at lower masses and higher redshifts than before, challenging standard formation scenarios. All data products, including the literature compilation, are made publicly available., Comment: 15 pages, 6 Figures + Appendix. Accepted for publication in ApJ on Feb, 9. Data release: - Reduced HST+JWST mosaics + photometric catalogs and Eazy-py modeling: https://erda.ku.dk/archives/7166d013c1ca1371aac3c57b9e73190d/published-archive.html - Supplementary material and tables: https://zenodo.org/record/7614908#.Y-4ZruzMLmE - MAST: https://doi.org/10.17909/g3nt-a370 - See also Gould et al. 2023

We present CEERS JWST/NIRCam imaging of a massive galaxy group at z=1.85, to explore the early JWST view on massive group formation in the distant Universe. The group contains >16 members (including 6 spectros. confirmations) down to log10(Mstar/Msun)=8.5, including the brightest group galaxy (BGG) in the process of actively assembling at this redshift. The BGG is comprised of multiple merging components extending ~3.6" (30kpc) across the sky. The BGG contributes 69% of the group's total galactic stellar mass, with one of the merging components containing 76% of the total mass of the BGG and a SFR>1810Msun/yr. Most importantly, we detect intra-halo light (IHL) in several HST and JWST/NIRCam bands, allowing us to construct a state-of-the-art rest-frame UV-NIR Spectral Energy Distribution of the IHL for the first time at this high redshift. This allows stellar population characterisation of both the IHL and member galaxies, as well as the morphology distribution of group galaxies vs. their star-formation activity when coupled with Herschel data. We create a stacked image of the IHL, giving us a sensitivity to extended emission of 28.5 mag/arcsec2 at rest-frame 1um. We find that the IHL is extremely dust poor (Av~0), containing an evolved stellar population of log10(t50/yr)=8.8, corresponding to a formation epoch for 50% of the stellar material 0.63Gyr before z=1.85. There is no evidence of ongoing star-formation in the IHL. The IHL in this group at z=1.85 contributes ~10% of the total stellar mass, comparable with what is observed in local clusters. This suggests that the evolution of the IHL fraction is more self-similar with redshift than predicted by some models, challenging our understanding of IHL formation during the assembly of high-redshift clusters. JWST is unveiling a new side of group formation at this redshift, which will evolve into Virgo-like structures in the local Universe., Comment: 14 pages + appendix, 11 figures, 4 tables. Accepted to A&A on 15th May 2023

Using the novel $JWST$/NIRCam observations in the Abell 2744 field, we present a first spatially resolved overview of an $HST$-dark galaxy, spectroscopically confirmed at $z=2.58$ with magnification $\mu\approx1.9$. While being largely invisible at $\sim$1 $\mu$m with NIRCam, except for sparse clumpy sub-structures, the object is well-detected and resolved in the long-wavelength bands with a spiral shape clearly visible in F277W. By combining ancillary ALMA and $Herschel$ data, we infer that this object is an edge-on dusty spiral with an intrinsic stellar mass log$(M_*/M_\odot)\sim11.3$ and a dust-obscured SFR$\sim 300~M_\odot$~yr$^{-1}$. A massive quiescent galaxy (log$(M_*/M_\odot)\sim10.8$) with tidal features lies 2\farcs{0} away ($r$$\sim$9 kpc), at a consistent redshift as inferred by $JWST$ photometry, indicating a potential major merger. The dusty spiral lies on the main-sequence of star formation, and shows high dust attenuation in the optical ($34$ would be only detectable in F356W and F444W in UNCOVER-like survey, and become totally $JWST$-dark at $z\sim6$. This suggests that detecting highly attenuated galaxies in the Epoch of Reionization might be a challenging task for $JWST$., Comment: 15 pages, 5 figures, 1 table. Accepted to ApJL

We report the discovery of a compact group of galaxies, CGG-z5, at z~5.2 in the EGS field covered by the JWST/CEERS survey. CGG-z5 was selected as the highest overdensity of galaxies at z>2 in recent JWST public surveys and it consists of six candidate members lying within a projected area of $1.5"\times3"$ (10$\times$20~kpc$^2$). All group members are HST/F435W and HST/F606W dropouts while securely detected in the JWST/NIRCam bands, yielding a narrow range of robust photometric redshifts $5.011$) at z~1. This implies that CGG-z5 could be a "proto-massive galaxy" captured during a short-lived phase of massive galaxy formation., Comment: Published in A&A Letter

The new capabilities that JWST offers in the near- and mid-infrared (IR) are used to investigate in unprecedented detail the nature of optical/near-IR faint, mid-IR bright sources, HST-dark galaxies among them. We gather JWST data from the CEERS survey in the EGS, jointly with HST data, and analyze spatially resolved optical-to-mid-IR spectral energy distributions (SEDs) to estimate both photometric redshifts in 2 dimensions and stellar populations properties in a pixel-by-pixel basis. We select 138 galaxies with F150W-F356W>1.5 mag, F356W<27.5 mag. The nature of these sources is threefold: (1) 71% are dusty star-forming galaxies at 2100 Gyr^-1); (2) 18% are quiescent/dormant (i.e., subject to reignition and rejuvenation) galaxies at 3

Our current understanding of the cosmic star formation history at z>3 is primarily based on UV-selected galaxies (i.e., LBGs). Recent studies of H-dropouts have revealed that we may be missing a large proportion of star formation that is taking place in massive galaxies at z>3. In this work, we extend the H-dropout criterion to lower masses to select optically dark/faint galaxies (OFGs), in order to complete the census between LBGs and H-dropouts. Our criterion (H> 26.5 mag & [4.5] < 25 mag) combined with a de-blending technique is designed to select not only extremely dust-obscured massive galaxies but also normal star-forming galaxies. In total, we identified 27 OFGs at z_phot > 3 (z_med=4.1) in the GOODS-ALMA field, covering a wide distribution of stellar masses with log($M_{\star}$/$M_{\odot}$) = 9.4-11.1. We find that up to 75% of the OFGs with log($M_{\star}$/$M_{\odot}$) = 9.5-10.5 were neglected by previous LBGs and H-dropout selection techniques. After performing stacking analyses, the OFGs exhibit shorter gas depletion timescales, slightly lower gas fractions, and lower dust temperatures than typical star-forming galaxies. Their SFR_tot (SFR_ IR+SFR_UV) is much larger than SFR_UVcorr (corrected for dust extinction), with SFR_tot/SFR_UVcorr = $8\pm1$, suggesting the presence of hidden dust regions in the OFGs that absorb all UV photons. The average dust size measured by a circular Gaussian model fit is R_e(1.13 mm)=1.01$\pm$0.05 kpc. We find that the cosmic SFRD at z>3 contributed by massive OFGs is at least two orders of magnitude higher than the one contributed by equivalently massive LBGs. Finally, we calculate the combined contribution of OFGs and LBGs to the cosmic SFRD at z=4-5 to be 4 $\times$ 10$^{-2}$ $M_{\odot}$ yr$^{-1}$Mpc$^{-3}$, which is about 0.15 dex (43%) higher than the SFRD derived from UV-selected samples alone at the same redshift., Comment: 26 pages, 14 figures, 8 tables. Accepted for publication in A&A

We report a galaxy group candidate HPC1001 at $z\approx3.7$ in the COSMOS field. This structure was selected as a high galaxy overdensity at $z>3$ in the COSMOS2020 catalog. It contains ten candidate members, of which eight are assembled in a $10''\times10''$ area with the highest sky density among known protoclusters and groups at $z>3$. Four out of ten sources were also detected at 1.2$~$mm with Atacama Large Millimeter Array continuum observations. Photometric redshifts, measured by four independent methods, fall within a narrow range of $3.5

Lyman Break Galaxy (LBG) candidates at z>10 are rapidly being identified in JWST/NIRCam observations. Due to the (redshifted) break produced by neutral hydrogen absorption of rest-frame UV photons, these sources are expected to drop out in the bluer filters while being well detected in redder filters. However, here we show that dust-enshrouded star-forming galaxies at lower redshifts (z<7) may also mimic the near-infrared (near-IR) colors of z>10 LBGs, representing potential contaminants in LBG candidate samples. First, we analyze CEERS-DSFG-1, a NIRCam dropout undetected in the F115W and F150W filters but detected at longer wavelengths. Combining the JWST data with (sub)millimeter constraints, including deep NOEMA interferometric observations, we show that this source is a dusty star-forming galaxy (DSFG) at z~5.1. We also present a tentative 2.6sigma SCUBA-2 detection at 850um around a recently identified z~16 LBG candidate in the same field and show that, if the emission is real and associated with this candidate, the available photometry is consistent with a z~5 dusty galaxy with strong nebular emission lines despite its blue near-IR colors. Further observations on this candidate are imperative to mitigate the low confidence of this tentative submillimeter emission and its positional uncertainty. Our analysis shows that robust (sub)millimeter detections of NIRCam dropout galaxies likely imply z=4-6 redshift solutions, where the observed near-IR break would be the result of a strong rest-frame optical Balmer break combined with high dust attenuation and strong nebular line emission, rather than the rest-frame UV Lyman break. This provides evidence that DSFGs may contaminate searches for ultra high-redshift LBG candidates from JWST observations., Comment: Published in The Astrophysical Journal Letters (updated to match the published version)

We report the discovery of a candidate galaxy with a photo-z of z~12 in the first epoch of the JWST Cosmic Evolution Early Release Science (CEERS) Survey. Following conservative selection criteria we identify a source with a robust z_phot = 11.8^+0.3_-0.2 (1-sigma uncertainty) with m_F200W=27.3, and >7-sigma detections in five filters. The source is not detected at lambda < 1.4um in deep imaging from both HST and JWST, and has faint ~3-sigma detections in JWST F150W and HST F160W, which signal a Ly-alpha break near the red edge of both filters, implying z~12. This object (Maisie's Galaxy) exhibits F115W-F200W > 1.9 mag (2-sigma lower limit) with a blue continuum slope, resulting in 99.6% of the photo-z PDF favoring z > 11. All data quality images show no artifacts at the candidate's position, and independent analyses consistently find a strong preference for z > 11. Its colors are inconsistent with Galactic stars, and it is resolved (r_h = 340 +/- 14 pc). Maisie's Galaxy has log M*/Msol ~ 8.5 and is highly star-forming (log sSFR ~ -8.2 yr^-1), with a blue rest-UV color (beta ~ -2.5) indicating little dust though not extremely low metallicity. While the presence of this source is in tension with most predictions, it agrees with empirical extrapolations assuming UV luminosity functions which smoothly decline with increasing redshift. Should followup spectroscopy validate this redshift, our Universe was already aglow with galaxies less than 400 Myr after the Big Bang., Comment: 17 pages, 6 figures, 3 tables, ApJL in press. Summary of changes from original submission: Improvements in astrometry generated a weak detection in F150W that reduces the photo-z to 11.8 but does not increase the likelihood of lower-z solutions. A full discussion of changes from the original version is available at: https://web.corral.tacc.utexas.edu/ceersdata/papers/Maisie_update.pdf

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

One of the most prominent features of galaxy clusters is the presence of a dominant population of massive ellipticals in their cores. Stellar archaeology suggests that these gigantic beasts assembled most of their stars in the early Universe via starbursts. However, the role of dense environments and their detailed physical mechanisms in triggering starburst activities remain unknown. Here we report spatially resolved Atacama Large Millimeter/submillimeter Array (ALMA) observations of the CO $J= 3-2$ emission line, with a resolution of about 2.5 kiloparsecs, toward a forming galaxy cluster core with starburst galaxies at $z=2.51$. In contrast to starburst galaxies in the field often associated with galaxy mergers or highly turbulent gaseous disks, our observations show that the two starbursts in the cluster exhibit dynamically cold (rotation-dominated) gas-rich disks. Their gas disks have extremely low velocity dispersion ($\sigma_{\mathrm{0}} \sim 20-30$ km s$^{-1}$), which is three times lower than their field counterparts at similar redshifts. The high gas fraction and suppressed velocity dispersion yield gravitationally unstable gas disks, which enables highly efficient star formation. The suppressed velocity dispersion, likely induced by the accretion of corotating and coplanar cold gas, might serve as an essential avenue to trigger starbursts in massive halos at high redshifts., Comment: 19 pages, 14 figures, 3 tables. Accepted for publication in A&A

We analyse the far-infrared properties of $\sim$ 5,000 star-forming galaxies at $z<4.5$, drawn from the deepest, super-deblended catalogues in the GOODS-N and COSMOS fields. We develop a novel panchromatic SED fitting algorithm, $\texttt{Stardust}$, that models the emission from stars, AGN, and infrared dust emission, without relying on energy balance assumptions. Our code provides robust estimates of the UV-optical and FIR physical parameters, such as the stellar mass ($M_*$), dust mass ($M_{\rm dust}$), infrared luminosities ($L_{\rm IR}$) arising from AGN and star formation activity, and the average intensity of the interstellar radiation field ($\langle U \rangle$). Through a set of simulations we quantify the completeness of our data in terms of $M_{\rm dust}$, $L_{\rm IR}$ and $\langle U \rangle$, and subsequently characterise the distribution and evolution of these parameters with redshift. We focus on the dust-to-stellar mass ratio ($f_{\rm dust}$), which we parametrise as a function of cosmic age, stellar mass, and specific star formation rate. The $f_{\rm dust}$ is found to increase by a factor of 10 from $z=0$ to $z=2$ and appears to remain flat at higher$-z$, mirroring the evolution of the gas fraction. We also find a growing fraction of warm to cold dust with increasing distance from the main sequence, indicative of more intense interstellar radiation fields, higher star formation efficiencies and more compact star forming regions for starburst galaxies. Finally, we construct the dust mass functions (DMF) of star-forming galaxies up to $z=1$ by transforming the stellar mass function to DMF through the scaling relations derived here. The evolution of $f_{\rm dust}$ and the recovered DMFs are in good agreement with the theoretical predictions of the Horizon-AGN and IllustrisTNG simulations., Comment: 39 pages, 21 figures, 3 tables. Accepted to ApJ July 10, 2021

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., Comment: 16 pages, 10 figures + Appendix, accepted in A&A on August 28th, 2021. The data compilation will be available on Vizier. (Minor) updates for the AGN population with respect to the data release of Valentino et al. 2020, A&A, 641, 155. For an early access, please contact the corresponding author. Abstract slightly modified to adjust to arXiv's requirements

Deep ALMA and HST observations reveal the presence of a quenched massive galaxy within the $z=2.91$ galaxy group RO-1001. With a mass-weighted stellar age of $1.6 \pm 0.4 \,$Gyr this galaxy is one of the oldest known at $z\sim3$, implying that most of its $10^{11}\rm \, M_{\odot}$ of stars were rapidly formed at $z>6$--8. This is a unique example of the predominantly passive evolution of a galaxy over at least $3

We study the interstellar medium (ISM) properties as a function of the molecular gas size for 77 infrared-selected galaxies at $z \sim 1.3$. Molecular gas sizes are measured on ALMA images that combine CO(2-1), CO(5-4) and underlying continuum observations, and include CO(4-3), CO(7-6)+[CI]($^3 P_2-^3P_1$), [CI]($^3 P_1-^3P_0$) observations for a subset of the sample. The $\gtrsim 46 \%$ of our galaxies have a compact molecular gas reservoir, and lie below the optical disks mass-size relation. Compact galaxies on and above the main sequence have higher CO excitation and star formation efficiency than galaxies with extended molecular gas reservoirs, as traced by CO(5-4)/CO(2-1) and CO(2-1)/$L_{\rm IR, SF}$ ratios. Average CO+[CI] spectral line energy distributions indicate higher excitation in compacts relative to extended sources. Using CO(2-1) and dust masses as molecular gas mass tracers, and conversion factors tailored to their ISM conditions, we measure lower gas fractions in compact main-sequence galaxies compared to extended sources. We suggest that the sub-millimetre compactness, defined as the ratio between the molecular gas and the stellar size, is an unavoidable information to be used with the main sequence offset to describe the ISM properties of galaxies, at least above $M_{\star} \geqslant 10^{10.6}$ M$_{\odot}$, where our observations fully probe the main sequence scatter. Our results are consistent with mergers driving the gas in the nuclear regions, enhancing the CO excitation and star formation efficiency. Compact main-sequence galaxies are consistent with being an early post-starburst population following a merger-driven starburst episode, stressing the important role of mergers in the evolution of massive galaxies., Comment: Accepted on MNRAS on October 11th 2021. Added Figure 10, 11, 12 following referee's comments. Abstract slightly modified to fit to arXiv's requirements

We report the detection of multiple faint radio sources, that we identify as AGN-jets, within CLJ1449+0856 at z=2 using 3 GHz VLA observations. We study the effects of radio-jet based kinetic feedback at high redshifts, which has been found to be crucial in low redshift clusters to explain the observed thermodynamic properties of their ICM. We investigate this interaction at an epoch featuring high levels of AGN activity and a transitional phase of ICM in regards to the likelihood of residual cold-gas accretion. We measure a total flux of $\rm 30.6 \pm 3.3~\mu Jy$ from the 6 detected jets. Their power contribution is estimated to be $1.2 ~(\pm 0.6)~ \times 10^{44} ~\rm ergs~ s^{-1}$, although this value could be up to $4.7 ~ \times 10^{44} ~\rm ergs~ s^{-1}$. This is a factor $\sim 0.25 - 1.0$ of the previously estimated instantaneous energy injection into the ICM of CLJ1449+0856 from AGN outflows and star formation, that have already been found to be sufficient in globally offsetting the cooling flows in the cluster core. In line with the already detected abundance of star formation, this mode of feedback being distributed over multiple sites, contrary to a single central source observed at low redshifts, points to accretion of gas into the cluster centre. This also suggests a 'steady state' of the cluster featuring non cool-core like behaviour. Finally, we also examine the TIR-radio luminosity ratio for the known sample of galaxies within the cluster core and find that dense environments do not have any serious consequence on the compliance of galaxies to the IR-radio correlation.

We examine the Fundamental Plane (FP) and mass-to-light ratio ($M/L$) scaling relations using the largest sample of massive quiescent galaxies at $1.511.26$, subset of 8 quiescent galaxies at $z>2$, from Stockmann et al. (2020), we show that they cannot passively evolve to the local Coma cluster relation alone and must undergo significant structural evolution to mimic the sizes of local massive galaxies. The evolution of the FP and $M/L$ scaling relations, from $z=2$ to present-day, for this subset are consistent with passive aging of the stellar population and minor merger structural evolution into the most massive galaxies in the Coma cluster and other massive elliptical galaxies from the MASSIVE Survey. Modeling the luminosity evolution from minor merger added stellar populations favors a history of merging with "dry" quiescent galaxies., Comment: 17 pages, 6 figures

We present the detection of multiple CO line transitions with ALMA in a few tens of infrared-selected galaxies on and above the main sequence at z=1.1-1.7. We reliably detected the emission of CO(5-4), CO(2-1), and CO(7-6)+[CI](2-1) in 50, 33, and 13 galaxies, respectively, and we complemented this information with available CO(4-3) and [CI](1-0) fluxes for part of the sample, and modeling of the optical-to-mm SEDs. We retrieve a quasi-linear relation between LIR and CO(5-4) or CO(7-6) for main-sequence galaxies and starbursts, corroborating the hypothesis that these transitions can be used as SFR tracers. We find the CO excitation to steadily increase as a function of the star formation efficiency, the mean intensity of the radiation field warming the dust, the surface density of SFR, and, less distinctly, with the distance from the main sequence. This adds to the tentative evidence for higher excitation of the CO+[CI] SLED of starbursts relative to that for main-sequence objects, where the dust opacities play a minor role in shaping the high-J CO transitions in our sample. However, the distinction between the average SLED of upper main-sequence and starburst galaxies is blurred, driven by a wide variety of intrinsic shapes. LVG radiative transfer modeling demonstrates the existence of a highly excited component that elevates the CO SLED of high-redshift main-sequence and starbursting galaxies above the typical values observed in the disk of the Milky Way. This excited component is dense and it encloses ~50% of the total molecular gas mass in main-sequence objects. We interpret the observed trends involving the CO excitation as mainly driven by a combination of large SFRs and compact sizes, naturally connected with enhanced dense molecular gas fractions and higher dust and gas temperatures, due to increasing UV radiation fields, cosmic ray rates, and dust/gas coupling. [Abridged], Comment: 18 pages, 9 figures + Appendix, accepted in A&A. The data compilation will be released as a .fits table on the journal's website. For an early access, please contact the corresponding author. The plots of spectra and far-infrared SEDs will be available on the journal's website and temporarily at this link: https://bit.ly/spectra_and_seds_valentino_2020c

As one of the prime contributors to the interstellar medium energy budget, magnetic fields naturally play a part in shaping the evolution of galaxies. Galactic magnetic fields can originate from strong primordial magnetic fields provided these latter remain below current observational upper limits. To understand how such magnetic fields would affect the global morphological and dynamical properties of galaxies, we use a suite of high-resolution constrained transport magneto-hydrodynamic cosmological zoom simulations where we vary the initial magnetic field strength and configuration along with the prescription for stellar feedback. We find that strong primordial magnetic fields delay the onset of star formation and drain the rotational support of the galaxy, diminishing the radial size of the galactic disk and driving a higher amount of gas towards the centre. This is also reflected in mock UVJ observations by an increase in the light profile concentration of the galaxy. We explore the possible mechanisms behind such a reduction in angular momentum, focusing on magnetic braking. Finally, noticing that the effects of primordial magnetic fields are amplified in the presence of stellar feedback, we briefly discuss whether the changes we measure would also be expected for galactic magnetic fields of non-primordial origin., Comment: 23 pages, 15 figures, accepted for publication in MNRAS

We present a detailed analysis of a large sample of spectroscopically confirmed ultra-massive quiescent galaxies (${\rm{log}}(M_{\ast}/M_{\odot})\sim11.5$) at $z\gtrsim2$. This sample comprises 15 galaxies selected in the COSMOS and UDS fields by their bright K-band magnitudes and followed up with VLT/X-shooter spectroscopy and HST/WFC3 $H_{F160W}$ imaging. These observations allow us to unambiguously confirm their redshifts ascertain their quiescent nature and stellar ages, and to reliably assess their internal kinematics and effective radii. We find that these galaxies are compact, consistent with the high mass end of the mass-size relation for quiescent galaxies at $z=2$. Moreover, the distribution of the measured stellar velocity dispersions of the sample is consistent with the most massive local early-type galaxies from the MASSIVE Survey showing that evolution in these galaxies, is dominated by changes in size. The HST images reveal, as surprisingly high, that $40\ \%$ of the sample have tidal features suggestive of mergers and companions in close proximity, including three galaxies experiencing ongoing major mergers. The absence of velocity dispersion evolution from $z=2$ to $0$, coupled with a doubling of the stellar mass, with a factor of four size increase and the observed disturbed stellar morphologies support dry minor mergers as the primary drivers of the evolution of the massive quiescent galaxies over the last 10 billion years., Comment: 30 pages, 10 figures, accepted in ApJ

We present the first stellar velocity dispersion measurement of a massive quenching galaxy at z=4.01. The galaxy is first identified as a massive z>~4 galaxy with suppressed star formation from photometric redshifts based on deep multi-band data in the UKIDSS Ultra Deep Survey field. A follow-up spectroscopic observation with MOSFIRE on Keck revealed strong multiple absorption features, which are identified as Balmer absorption lines, giving a secure redshift of z=4.01. Thanks to the high S/N of the spectrum, we are able to estimate the stellar velocity dispersion, sigma=268+/-59 km/s. This velocity dispersion is consistent with that of massive galaxies today, implying no significant evolution in stellar velocity dispersion over the last 12 Gyr. Based on an upper limit on its physical size from deep optical images (r_eff<1.3 kpc), we find that its dynamical mass is consistent with the stellar mass inferred from photometry. Furthermore, the galaxy is located on the mass fundamental plane extrapolated from lower redshift galaxies. Combining all these results, we find that the velocity dispersion does not significantly evolve with redshift, although the size and mass of massive quenched galaxies do. This suggests that the mass in the core of massive galaxies does not evolve significantly, while most of the mass growth occurs in the outskirts of the galaxies, which also increases the size. This picture is consistent with a two-phase formation scenario in which mass and size growth is due to accretion in the outskirts of galaxies via mergers., Comment: Published in the Astrophysical Journal letters. Fixed an error in dynamical mass

We report two secure ($z=3.775, 4.012$) and one tentative ($z\approx3.767$) spectroscopic confirmations of massive and quiescent galaxies through $K$-band observations with Keck/MOSFIRE and VLT/X-Shooter. The stellar continuum emission, the absence of strong nebular emission lines and the lack of significant far-infrared detections confirm the passive nature of these objects, disfavoring the alternative solution of low-redshift dusty star-forming interlopers. We derive stellar masses of $\mathrm{log}(M_{\star}/M_\odot)\sim11$ and ongoing star formation rates placing these galaxies $\gtrsim 1-2$ dex below the main sequence at their redshifts. The adopted parametrization of the star formation history suggests that these sources experienced a strong ($\langle \rm SFR \rangle \sim 1200-3500\,M_\odot\,\mathrm{yr}^{-1}$) and short ($\sim 50$ Myr) burst of star formation, peaking $\sim 150-500$ Myr before the time of observation, all properties reminiscent of the characteristics of sub-millimeter galaxies (SMGs) at $z>4$. We investigate this connection by comparing the comoving number densities and the properties of these two populations. We find a fair agreement only with the deepest sub-mm surveys detecting not only the most extreme starbursts, but also more normal galaxies. We support these findings by further exploring the Illustris-TNG cosmological simulation, retrieving populations of both fully quenched massive galaxies at $z\sim3-4$ and SMGs at $z\sim4-5$, with number densities and properties in agreement with the observations at $z\sim3$, but in increasing tension at higher redshift. Nevertheless, as suggested by the observations, not all the progenitors of quiescent galaxies at these redshifts shine as bright SMGs in their past and, similarly, not all bright SMGs quench by $z\sim3$, both fractions depending on the threshold assumed to define the SMGs themselves., Comment: 20 pages (+appendix), 11 figures, accepted for publication in the Astrophysical Journal (ApJ). This work comes with a companion letter by Tanaka, Valentino, Toft et al. 2019, ApJL, 885, L34

We present the rest-frame optical sizes of massive quiescent galaxies (QGs) at $z\sim4$ measured at $K'$-band with the Infrared Camera and Spectrograph (IRCS) and AO188 on the Subaru telescope. Based on a deep multi-wavelength catalog in the Subaru XMM-Newton Deep Survey Field (SXDS), covering a wide wavelength range from the $u$-band to the IRAC $8.0\mu m$ over 0.7 deg$^2$, we evaluate photometric redshift to identify massive ($M_{\star}\sim10^{11}\ M_\odot$) galaxies with suppressed star formation. These galaxies show a prominent 4000$\rm \AA$ break feature at $z\sim4$, suggestive of an evolved stellar population. We then conduct follow-up $K'$-band imaging with adaptive optics for the five brightest galaxies ($K_{AB,total}=22.5\sim23.4$). Compared to lower redshift ones, QGs at $z\sim4$ have smaller physical sizes of effective radii $r_{eff}=0.2$ to $1.8$ kpc. The mean size measured by stacking the four brightest objects is $r_{eff}=0.7\rm\ kpc$. This is the first measurement of the rest-frame optical sizes of QGs at $z\sim4$. We evaluate the robustness of our size measurements using simulations and find that our size estimates are reasonably accurate with an expected systematic bias of $\sim0.2$ kpc. If we account for the stellar mass evolution, massive QGs at $z\sim4$ are likely to evolve into the most massive galaxies today. We find their size evolution with cosmic time in a form of $\log(r_e/{\rm kpc})= -0.44+1.77 \log(t/\rm Gyr)$. Their size growth is proportional to the square of stellar mass, indicating the size-stellar mass growth driven by minor dry mergers., Comment: 15 pages, 11 figures, ApJ accepted

The Hubble Ultra Deep field (HUDF) is the deepest region ever observed with the Hubble Space Telescope. With the main objective of unveiling the nature of galaxies up to $z \sim 7-8$, the observing and reduction strategy have focused on the properties of small and unresolved objects, rather than the outskirts of the largest objects, which are usually over-subtracted. We aim to create a new set of WFC3/IR mosaics of the HUDF using novel techniques to preserve the properties of the low surface brightness regions. We created ABYSS: a pipeline that optimises the estimate and modelling of low-level systematic effects to obtain a robust background subtraction. We have improved four key points in the reduction: 1) creation of new absolute sky flat fields, 2) extended persistence models, 3) dedicated sky background subtraction and 4) robust co-adding. The new mosaics successfully recover the low surface brightness structure removed on the previous HUDF published reductions. The amount of light recovered with a mean surface brightness dimmer than $\overline{\mu}=26$ mar arcsec$^{-2}$ is equivalent to a m=19 mag source when compared to the XDF and a m=20 mag compared to the HUDF12. We present a set of techniques to reduce ultra-deep images ($\mu>32.5$ mag arcsec$^{-2}$, $3\sigma$ in $10\times10$ arcsec boxes), that successfully allow to detect the low surface brightness structure of extended sources on ultra deep surveys. The developed procedures are applicable to HST, JWST, EUCLID and many other space and ground-based observatories. We will make the final ABYSS WFC3/IR HUDF mosaics publicly available at http://www.iac.es/proyecto/abyss/., Comment: Published in Astronomy & Astrophysics

We present morphological properties of dusty star-forming galaxies at z=1-3 determined with high-resolution (FWHM~0"19) Atacama Large Milllimeter/submilimeter Array (ALMA) 1-mm band maps of our ASAGAO survey covering a 26-arcmin^2 area in GOODS-S. In conjunction with the ALMA archival data, the present sample consists of 42 ALMA sources with a wide rest-frame far-infrared (FIR) luminosity L_FIR range of ~10^11-10^13 Lo. To obtain an average rest-frame FIR profile, we perform individual measurements and careful stacking of the ALMA sources using the uv-visibility method that includes positional-uncertainty and smoothing-effect evaluations through Monte-Carlo simulations. We find that the dusty star-forming galaxies have the average FIR-wavelength Sersic index and effective radius of n_FIR=1.2+/-0.2 and R_e,FIR=1.0-1.3 kpc, respectively, additionally with a point source at the center, indicative of the existence of AGN. The average FIR profile agrees with a morphology of an exponential-disk clearly distinguished from a spheroidal profile (Sersic index of 4). We also examine the rest-frame optical Sersic index n_opt and effective radius R_e,opt with the deep Hubble Space Telescope (HST) images. Interestingly, we obtain n_opt=0.9+/-0.3 (~n_FIR) and R_e,opt=3.2+/-0.6 kpc (>R_e,FIR), suggesting that the FIR-emitting disk is embedded within a larger stellar disk. The rest-frame UV and FIR data of HST and ALMA provide us a radial surface density profile of the total star-formation rate (SFR), where the FIR SFR dominates over the UV SFR at the center. Under the simple assumption of a constant SFR, a compact stellar distribution found in z~1-2 compact quiescent galaxies (cQGs) is well reproduced, while a spheroidal stellar morphology of cQGs (n_opt=4) cannot, suggestive of other important mechanisms such as dynamical dissipation., Comment: 13 pages, 6 figures, ApJ in press

At redshift z = 2, when the Universe was just three billion years old, half of the most massive galaxies were extremely compact and had already exhausted their fuel for star formation(1-4). It is believed that they were formed in intense nuclear starbursts and that they ultimately grew into the most massive local elliptical galaxies seen today, through mergers with minor companions(5,6), but validating this picture requires higher-resolution observations of their centres than is currently possible. Magnification from gravitational lensing offers an opportunity to resolve the inner regions of galaxies(7). Here we report an analysis of the stellar populations and kinematics of a lensed z = 2.1478 compact galaxy, which surprisingly turns out to be a fast-spinning, rotationally supported disk galaxy. Its stars must have formed in a disk, rather than in a merger-driven nuclear starburst(8). The galaxy was probably fed by streams of cold gas, which were able to penetrate the hot halo gas until they were cut off by shock heating from the dark matter halo(9). This result confirms previous indirect indications(10-13) that the first galaxies to cease star formation must have gone through major changes not just in their structure, but also in their kinematics, to evolve into present-day elliptical galaxies., Comment: Nature (June 22, 2017)

Several studies of nearby active galactic nuclei (AGN) have shown that the soft X-ray emission presents a size and morphology that resembles that of the narrow-line region (NLR) traced by [O III]. Since the NLR is mainly constituted by gas photoionised by the AGN, it seems logical to assume that this is also the primary source of the soft X-ray emission. However, these results are based on individual sources or small samples, particularly focused on type-2 Seyfert galaxies. Very little has been said concerning other types of AGN. The purpose of this work is to compare the circumnuclear morphologies of soft X-ray and [O III] images to test whether they match in different optical classes of AGN. Our sample is composed of 27 AGN: nine type-1 Seyferts, 10 type-2 Seyferts, and eight low ionisation nuclear emission-line regions (LINERs). We find a good match in 100% of the type-2 Seyferts in our sample. This correspondence is less frequent in type-1 Seyferts (22%) and it is not seen in LINERs. The good resemblance in type-2 Seyferts constitutes an evidence for a common physical origin. We argue that the lack of correspondence in type-1 Seyferts might be due to the line of sight perpendicular to the accretion disk. Based on the morphologies of the eight LINERs in our sample, we discard a common origin for the soft X-ray and [O III] emissions in these objects. Regarding the X-ray properties, both high column density and hard X-ray luminosity are associated with matched morphologies., Comment: Accepted by MNRAS. 20 pages, 9 figures

Studies measuring the star formation rate density, luminosity function, and properties of star-forming galaxies are numerous. However, it exists a gap at $0.5