Search

We conduct a study of the gas kinematics of two quasar host galaxies at $z\gtrsim6$ traced by the [CII] emission line using ALMA. By combining deep observations at both low and high resolution, we recover the diffuse emission, resolve its structure, and measure the rotation curves from the inner region of the galaxy to its outskirts using DysmalPy and 3DBarolo. Assuming that both galaxies exhibit disk rotation driven by the gravitational potential of the galaxy, we find that the best-fit disk models have a $V_{\rm rot}/\sigma \approx 2$ and inferred circular velocities out to $\sim$6-8 kpc scales, well beyond the likely stellar distribution. We then determine the mass profiles of each component (stars, gas, dark matter) with priors on the baryon and dark matter properties. We find relatively large dark matter fractions within their effective radii ($f_{\rm DM}(R

Star-forming clumps have been found to significantly influence the star formation of gas-rich $z>1$ galaxies. Using public data from JWST/NIRCam (COSMOS-Web) and ALMA (FMOS-COSMOS), we study a sample of 32 massive ($>10^{10.5}\,\rm M_{\odot}$) main-sequence galaxies at $z_{\rm spec}\sim1.5$ with $\sim0.3\,\rm kpc$ resolution. We create composite morphological models consisting of bulge, disk, and clumps to fully 'deconstruct' the galaxy images. With the resulting measurements of the flux and size of these components, we find the following: (I)The combined contribution of clumps is $1-30\%$ towards the net star formation rate (SFR) of the host while contributing $1-20\%$ to its stellar mass. The clumps show a correlation between their stellar mass and SFR, but have an increased specific-SFR (sSFR) relative to the star-forming main sequence, with offsets ranging from $0\lesssim\Delta\log\rm sSFR\lesssim 0.4$. They feature star formation surface densities of $10^{-2}-10^{2}\,\rm M_{\odot}/yr/kpc^{2}$, consistent with values observed in local star-forming and starburst galaxies. (II)The clumps span a large range of characteristic sizes ($r_{e}\sim0.1-1\,\rm kpc$) and stellar masses ($\sim 10^{8.0-9.5}\,\rm M_{\odot}$). We estimate a mass-size relation ($r_{e}\propto\rm M_{\star}^{\,0.52\pm0.07}$) along with a stellar mass function (slope, $\alpha=-1.85\pm 0.19$), both suggesting a hierarchical nature similar to that expected in star-forming regions in local galaxies. (III)Our measurements agree with the properties of stellar clumps in $z\gtrsim1$ lensed systems, bridging the gap between lensed and unlensed studies by detecting structures at sub-kpc scales.(IV)Clumps are found to be preferentially located along spiral features visible primarily in the residual rest-frame near-IR images. In conclusion, we present an observation-based, coherent picture of star-forming clumps at $z>1$., Comment: Published in MNRAS; 22 pages, 24 figures

Recent JWST observations have revealed the prevalence of spiral structures at $z > 1$. Unlike in the local Universe, the origin and the consequence of spirals at this epoch remain unexplored. We use public JWST/NIRCam data from the COSMOS-Web survey to map spiral structures in eight massive ($> 10^{10.5}\,\rm M_{\odot}$) star-forming galaxies at $z_{\rm spec} \sim 1.5$. We present a method for systematically quantifying spiral arms at $z>1$, enabling direct measurements of flux distributions. Using rest-frame near-IR images, we construct morphological models accurately tracing spiral arms. We detect offsets ($\sim 0.2 - 0.8\,\rm kpc$) between the rest-frame optical and near-IR flux distributions across most arms. Drawing parallels to the local Universe, we conclude that these offsets reflect the presence of density waves. For nine out of eighteen arms, the offsets indicate spiral shocks triggered by density waves. Five arms have offsets in the opposite direction and are likely associated with tidal interactions. For the remaining cases with no detected offsets, we suggest that stochastic 'clumpy' star formation is the primary driver of their formation. In conclusion, we find a multi-faceted nature of spiral arms at $z > 1$, similar to that in the local Universe., Comment: Accepted for publication in ApJ Letters; 12 pages, 6 figures

``Little Red Dots'' (LRDs) are an abundant high-redshift population newly discovered by the James Webb Space Telescope (JWST). They are characterized by a red color in the rest-frame optical band, compact morphology, and broad Balmer emission lines (${\rm FWHM} \gtrsim 1000~{\rm km\,s^{-1}}$) that suggest an AGN nature. Using a method of pixel-by-pixel color selection and relaxing the compactness criteria, we identify three of the first dual LRD candidates in the COSMOS-Web survey with projected separations of $0.\!\!^{\prime\prime}2-0.\!\!^{\prime\prime}4$ (1-2 pkpc at their photometric redshifts). A comparison between existing LRD samples and mock data reveals that the projected separations of these dual LRD candidates are unlikely to result from chance projections of objects at different redshifts. In one case (CW-B5-15958), the dual LRD includes two bright sources ($m_{\rm F444W}=24.3$ and $24.8$) with characteristic V-shape spectral energy distribution (SEDs) and photometric redshifts consistent with each other. We find that CW-B5-15958 has a faint off-centered component and a companion galaxy. In the other two dual systems, the brighter LRD exhibits a V-shape SED, while the fainter LRD ($m_{\rm F444W}\gtrsim26$) is undetected in both F115W and F150W. These discoveries suggest that the angular auto-correlation function (ACF) of LRDs exhibits a significant excess ($\sim3\times10^2$ times) on sub-arcsec (kilo-parsec) separations compared to the extrapolation of a power-law ACF of JWST-found AGNs measured over $10^{\prime\prime}-100^{\prime\prime}$. Follow-up spectroscopic confirmation of their redshifts and the construction of a larger sample are essential to advance our understanding of the evolution of supermassive black holes and the importance of mergers in the early universe., Comment: 12 pages, 5 figures, and 1 table. Comments are welcome

We investigate the deconvolved color profiles of 223 disk galaxies at redshifts of $z=1$-3 observed by the James Webb Space Telescope (JWST) as part of the Cosmic Evolution Early Release Science survey (CEERS). The filters were selected to approximate the rest-frame $B-Y$ color, which is used to identify U-shaped color profiles -- those becoming progressively bluer with increasing radius, then turning redder beyond a specific point. We find that 36% of Type II (down-bending) disks exhibit U-shaped color profiles with a minimum at or near the disk break. In contrast, no Type I (single-exponential) disks and only 9% of Type III (up-bending) disks show such a profile. The presence of U-shaped color profiles in Type II disks likely arises from the interplay between a star-formation threshold and spiral- or bar-driven secular radial migration of older stars outward. The fraction of Type II disks exhibiting a U-shaped color profile remains almost consistent across two redshift bins, $z=1$-$2$ and $z=2$-$3$, but is significantly lower than that observed in the local Universe, likely because the secular process of radial migration at high redshift may not have had sufficient time to significantly influence the disk structure. The absence of U-shaped color profiles in Type II disks could point to rapid rather than secular radial star migration potentially caused by violent clump instabilities, transporting both younger and older stars to the outer disk. Our results provide useful constraints on the formation and evolution models of disk galaxies in the early Universe., Comment: 5 pages, 4 figures. Accepted for publication in A&A

A recent study on the spectral energy distribution (SED) of AGN combined unobscured X-ray sources from the eROSITA eFEDS Survey with high quality optical imaging from Subaru's Hyper Suprime-Cam (HSC). The HSC data enabled accurate host galaxy subtraction as well as giving a uniform black hole mass estimator from the stellar mass. The resulting stacked optical/X-ray SEDs for black holes at fixed mass show a dramatic transition, where the dominating disc component in bright AGN evaporates into an X-ray hot plasma below $L/L_{\rm Edd}\sim 0.01$. The models fit to these datasets predicted the largest change in SED in the rest frame UV ($< 3000\,\r{A}$), but this waveband was not included in the original study. Here we use archival $u$-band and UV photometry to extend the SEDs into this range, and confirm the UV is indeed intrinsically faint in AGN below $L/L_{\rm Edd}\sim 0.01$ as predicted. This dramatic drop in UV photo-ionising flux is also seen from its effect on the broad emission lines. We stack the recently released SDSS DR18 optical spectra for this sample, and show that the broad H$\beta$ line disappears along with the UV bright component at $L/L_{\rm Edd}\sim 0.01$. This shows that there is a population of unobscured, X-ray bright, UV faint AGN which lack broad emission lines (true type 2 Seyferts)., Comment: 11 pages, 3 figure, 3 appendices. MNRAS accepted. Comments are very welcome!

The current paradigm for the co-evolution of galaxies and their supermassive black holes postulates that dust-obscured active galactic nuclei (AGNs) represent a transitional phase towards a more luminous and unobscured state. However, our understanding of dusty AGNs and their host galaxies at early cosmic times is inadequate due to observational limitations. Here, we present JWST observations of CID-931, an X-ray-detected AGN at a spectroscopic redshift of $z_{\rm spec}=4.91$. Multiband NIRCam imaging from the COSMOS-Web program reveals an unresolved red core, similar to JWST-discovered dusty AGNs. Strikingly, the red core is surrounded by at least eight massive star-forming clumps spread over $1.\!\!^{\prime\prime}6 \approx 10~{\rm kpc}$, each of which has a stellar mass of $10^9-10^{10}M_\odot$ and $\sim0.1-1~{\rm kpc}$ in radius. The whole system amounts to $10^{11}M_\odot$ in stellar mass, higher than typical star-forming galaxies at the same epoch. In this system, gas inflows and/or complex merger events may trigger clump formation and AGN activity thus leading to the rapid formation of a massive galaxy hosting a supermassive black hole. Future follow-up observations will provide new insights into the evolution of the galaxy-black hole relationship during such transitional phases in the early universe., Comment: 13 pages, 10 figures, submitted to PASJ

Understanding the rapid formation of supermassive black holes (SMBHs) in the early universe requires an understanding of how stellar mass grows in the host galaxies. Here, we perform an analysis of rest-frame optical spectra and imaging from JWST of two quasar host galaxies at z>6 which exhibit Balmer absorption lines. These features in the stellar continuum indicate a lack of young stars, similar to low-redshift post-starburst galaxies whose star formation was recently quenched. We find that the stellar mass (log(M_* / M_sun) > 10.6) of each quasar host grew in a starburst episode at redshift 7 or 8. One of the targets exhibits little ongoing star formation, as evidenced by the photometric signature of the Balmer break and a lack of spatially resolved H-alpha emission, placing it well below the star formation main sequence at z = 6. The other galaxy is transitioning to a quiescent phase; together, the two galaxies represent the most distant massive post-starburst galaxies known. The maturity of these two galaxies is further supported by the stellar velocity dispersions of their host galaxies, placing them slightly above the upper end of the local M_BH - sigma_* relation. The properties of our two post-starburst galaxies, each hosting an active SMBH with log(M_BH / M_sun) > 9, suggests that black holes played a major role in shaping the formation of the first massive galaxies in the Universe., Comment: 24 pages, 7 figures, submitted to a Nature journal

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

We present a study of the molecular gas in five closely-spaced ($R_{\perp}<20$ kpc) dual quasars ($L_{\rm bol}\gtrsim10^{44}~\mathrm{erg~s}^{-1}$) at redshifts $0.4

We introduce an improved method for decomposing the emission of active galactic nuclei (AGN) and their host galaxies using templates from principal component analysis (PCA). This approach integrates prior information from PCA with a penalized pixel fitting mechanism which improves the precision and effectiveness of the decomposition process. Specifically, we have reduced the degeneracy and over-fitting in AGN-host decomposition, particularly for those with low signal-to-noise ratios (SNR), where traditional methods tend to fail. By applying our method to 76,565 SDSS Data Release 16 quasars with $z<0.8$, we achieve a success rate of $\approx$ 94%, thus establishing the largest host-decomposed spectral catalog of quasars to date. Our fitting results consider the impact of the host galaxy on the overestimation of the AGN luminosity and black hole mass ($M_{\rm BH}$). Furthermore, we obtained stellar velocity dispersion ($\sigma_*$) measurements for 4,137 quasars. The slope of the $M_{\rm BH}-\sigma_*$ relation in this subsample is generally consistent with previous quasar studies beyond the local universe. Our method provides a robust and efficient approach to disentangle the AGN and host galaxy components across a wide range of SNRs and redshifts., Comment: 21 pages, 9 figures, 2 tables. Accepted by ApJ

JWST has revealed a population of compact and extremely red galaxies at $z>4$, which likely host active galactic nuclei (AGN). We present a sample of 434 ``little red dots'' (LRDs), selected from the 0.54 deg$^2$ COSMOS-Web survey. We fit galaxy and AGN SED models to derive redshifts and physical properties; the sample spans $z\sim5$-$9$ after removing brown dwarf contaminants. We consider two extreme physical scenarios: either LRDs are all AGN, and their continuum emission is dominated by the accretion disk, or they are all compact star-forming galaxies, and their continuum is dominated by stars. If LRDs are AGN-dominated, our sample exhibits bolometric luminosities $\sim10^{45-47}$ erg\,s$^{-1}$, spanning the gap between JWST AGN in the literature and bright, rare quasars. We derive a bolometric luminosity function (LF) $\sim100$ times the (UV-selected) quasar LF, implying a non-evolving black hole accretion density of $\sim10^{-4}$ M$_\odot$ yr$^{-1}$ Mpc$^{-3}$ from $z\sim2$-$9$. By contrast, if LRDs are dominated by star formation, we derive stellar masses $\sim10^{8.5-10}\,M_\odot$. MIRI/F770W is key to deriving accurate stellar masses; without it, we derive a mass function inconsistent with $\Lambda$CDM. The median stellar mass profile is broadly consistent with the maximal stellar mass surface densities seen in the nearby universe, though the most massive $\sim50$\% of objects exceed this limit, requiring substantial AGN contribution to the continuum. Nevertheless, stacking all available X-ray, mid-IR, far-IR/sub-mm, and radio data yields non-detections. Whether dominated by dusty AGN, compact star-formation, or both, the high masses/luminosities and remarkable abundance of LRDs implies a dominant mode of early galaxy/SMBH growth., Comment: 27 pages, 13 figures. Submitted to ApJ, comments welcome! Data access at https://github.com/hollisakins/akins24_cw

The structure of the accretion flow onto supermassive black holes (SMBH) is not well understood. Standard disc models match to zeroth order in predicting substantial energy dissipation within optically-thick material producing a characteristic strong blue/UV continuum. However they fail at reproducing more detailed comparisons to the observed spectral shapes along with their observed variability. Based on stellar mass black holes within our galaxy, accretion discs should undergo a transition into an X-ray hot, radiatively inefficient flow, below a (mass scaled) luminosity of $\sim 0.02\,L_{\rm{Edd}}$. While this has been seen in limited samples of nearby low-luminosity active galactic nuclei (AGN) and a few rare changing-look AGN, it is not at all clear whether this transition is present in the wider AGN population across cosmic time. A key issue is the difficulty in disentangling a change in spectral state from increased dust obscuration and/or host galaxy contamination, effectively drowning out the AGN emission. Here we use the new eROSITA eFEDS Survey to identify unobscured AGN from their X-ray emission, matched to excellent optical imaging from Subaru's Hyper Suprime-Cam; allowing the subtraction of the host galaxy contamination. The resulting, uncontaminated, AGN spectra reveal a smooth transition from a strongly disc dominated state in bright AGN, to the collapse of the disc into an inefficient X-ray plasma in the low luminosity AGN, with the transition occurring at $\sim 0.02\,L_{\rm{Edd}}$; revealing fundamental aspects of accretion physics in AGN., Comment: 11 pages, 5 figure, 3 appendices - Accepted for publication in MNRAS

We present ALMA [CII] 158 $\mu$m line and underlying far-infrared (FIR) continuum emission observations ($0''.57 \times 0''.46$ resolution) toward a quasar-quasar pair system recently discovered at $z = 6.05$ (Matsuoka et al. 2024). The quasar nuclei (C1 and C2) are faint ($M_{\rm 1450} \gtrsim -23$ mag), but we detect very bright [CII] emission bridging the 12 kpc between the two objects and extending beyond them (total luminosity $L_{\rm [CII]} \simeq 6 \times 10^9~L_\odot$). The [CII]-based total star formation rate of the system is $\sim 550~M_\odot$ yr$^{-1}$ (IR-based dust-obscured SFR is $\sim 100~M_\odot$ yr$^{-1}$), with a [CII]-based total gas mass of $\sim 10^{11}~M_\odot$. The dynamical masses of the two galaxies are large ($\sim 9 \times 10^{10}~M_\odot$ for C1 and $\sim 5 \times 10^{10}~M_\odot$ for C2). There is a smooth velocity gradient in [CII], indicating that these quasars are a tidally interacting system. We identified a dynamically distinct, fast [CII] component around C1: detailed inspection of the line spectrum there reveals the presence of a broad wing component, which we interpret as the indication of fast outflows with a velocity of $\sim 600$ km s$^{-1}$. The expected mass loading factor of the outflows, after accounting for multiphase gas, is $\gtrsim 2-3$, which is intermediate between AGN-driven and starburst-driven outflows. Hydrodynamic simulations in the literature predicted that this pair will evolve to a luminous ($M_{\rm 1450} \lesssim -26$ mag), starbursting ($\gtrsim 1000~M_\odot$ yr$^{-1}$) quasar after coalescence, one of the most extreme populations in the early universe., Comment: 16 pages incl. Appendix, 4 main figures, 1 table, 2 appendix figures. Accepted for publication in the Astrophysical Journal on 2024 June 11

We report the discovery of two quasars at a redshift of $z$ = 6.05, in the process of merging. They were serendipitously discovered from the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The quasars, HSC $J$121503.42$-$014858.7 (C1) and HSC $J$121503.55$-$014859.3 (C2), both have luminous ($>$10$^{43}$ erg s$^{-1}$) Ly$\alpha$ emission with a clear broad component (full width at half maximum $>$1000 km s$^{-1}$). The rest-frame ultraviolet (UV) absolute magnitudes are $M_{1450} = -23.106 \pm 0.017$ (C1) and $-22.662 \pm 0.024$ (C2). Our crude estimates of the black hole masses provide $\log (M_{\rm BH}/M_\odot) = 8.1 \pm 0.3$ in both sources. The two quasars are separated by 12 kpc in projected proper distance, bridged by a structure in the rest-UV light suggesting that they are undergoing a merger. This pair is one of the most distant merging quasars reported to date, providing crucial insight into galaxy and black hole build-up in the hierarchical structure formation scenario. A companion paper will present the gas and dust properties captured by Atacama Large Millimeter/submillimeter Array observations, which provide additional evidence for and detailed measurements of the merger and also demonstrate that the two sources are not gravitationally-lensed images of a single quasar., Comment: Published from ApJL. A companion paper (T. Izumi et al.) presents the associated ALMA observations

We report the identification of a quasar overdensity in the BOSSJ0210 field, dubbed Cosmic Himalayas, consisting of 11 quasars at $z=2.16-2.20$, the densest overdensity of quasars ($17\sigma$) in the $\sim$10,000 deg$^2$ of the Sloan Digital Sky Survey. We present the spatial distributions of galaxies and quasars and an HI absorption map of the intergalactic medium (IGM). On the map of 465 galaxies selected from the MAMMOTH-Subaru survey, we find two galaxy density peaks that do not fall on the quasar overdensity but instead exist at the northwest and southeast sides, approximately 25 $h^{-1}$ comoving-Mpc apart from the quasar overdensity. With a spatial resolution of 15 $h^{-1}$ comoving Mpc in projection, we produce a three-dimensional HI tomography map by the IGM Ly$\alpha$ forest in the spectra of 23 SDSS/eBOSS quasars behind the quasar overdensity. Surprisingly, the quasar overdensity coincides with neither an absorption peak nor a transmission peak of IGM HI but lies near the border separating opaque and transparent volumes, with the more luminous quasars located in an environment with lesser IGM HI. Hence remarkably, the overdensity region traced by the 11 quasars, albeit all in coherently active states, has no clear coincidence with peaks of galaxies or HI absorption densities. Current physical scenarios with mixtures of HI overdensities and quasar photoionization cannot fully interpret the emergence of Cosmic Himalayas, suggesting this peculiar structure is an excellent laboratory to unveil the interplay between galaxies, quasars, and the IGM., Comment: 19 pages, 11 figures, submitted to ApJ, comments are welcome

JWST is revealing a new remarkable population of high-redshift ($z\gtrsim4$), low-luminosity Active Galactic Nuclei (AGNs) in deep surveys and detecting the host galaxy stellar light in the most luminous and massive quasars at $z\sim 6$ for the first time. Latest results claim supermassive black holes (SMBHs) in these systems to be significantly more massive than expected from the local BH mass - stellar mass ($\mathcal{M}_{\rm BH} - \mathcal{M}_\star$) relation and that this is not due to sample selection effects. Through detailed statistical modeling, we demonstrate that the coupled effects of selection biases (i.e., finite detection limit and requirements on detecting broad lines) and measurement uncertainties in $\mathcal{M}_{\rm BH}$ and $\mathcal{M}_\star$ can in fact largely account for the reported offset and flattening in the observed $\mathcal{M}_{\rm BH} - \mathcal{M}_\star$ relation toward the upper envelope of the local relation, even for those at $\mathcal{M}_{\rm BH} < 10^8\,M_{\odot}$. We further investigate the possible evolution of the $\mathcal{M}_{\rm BH} - \mathcal{M}_\star$ relation at $z\gtrsim 4$ with careful treatment of observational biases and consideration of the degeneracy between intrinsic evolution and dispersion in this relation. The bias-corrected intrinsic $\mathcal{M}_{\rm BH} - \mathcal{M}_\star$ relation in the low-mass regime suggests that there might be a large population of low-mass BHs (${\rm log}\,\mathcal{M}_{\rm BH} \lesssim 5$), possibly originating from lighter seeds, remaining undetected or unidentified even in the deepest JWST surveys. These results have important consequences for JWST studies of BH seeding and the coevolution between SMBHs and their host galaxies at the earliest cosmic times., Comment: 21 pages, 9 figures, submitted to ApJ

Resolved stellar morphology of $z>1$ galaxies was inaccessible before JWST. This limitation, due to the impact of dust on rest-frame UV light, had withheld major observational conclusions required to understand the importance of clumps in galaxy evolution. Essentially independent of this issue, we use the rest-frame near-IR for a stellar-mass dependent clump detection method and determine reliable estimations of selection effects. We exploit publicly available JWST/NIRCam and HST/ACS imaging data from CEERS, to create a stellar-mass based picture of clumps in a mass-complete sample of 418 galaxies within a wide wavelength coverage of $0.5-4.6\,\mu$m and a redshift window of $1 < z < 2$. We find that a near-IR detection gives access to a larger, and possibly different, set of clumps within galaxies, with those also detected in UV making up only $28\%$. Whereas, $85\%$ of the UV clumps are found to have a near-IR counterpart. These near-IR clumps closely follow the UVJ classification of their respective host galaxies, with these hosts mainly populating the star-forming regime besides a fraction of them ($16\%$) that can be considered quiescent. The mass of the detected clumps are found to be within the range of $10^{7.5-9.5}\,\rm M_{\odot}$, therefore expected to drive gas into galaxy cores through tidal torques. The clump stellar mass function is found to have a slope of $-1.50 \pm 0.14$, indicating a hierarchical nature similar to that of star-forming regions in the local Universe. Finally, we observe a radial gradient of increasing clump mass towards the centre of galaxies., Comment: Revised arXiv version after moderate revisions and publication in MNRAS. 17 pages, 16 figures

Our knowledge of relations between supermassive black holes and their host galaxies at $z\gtrsim1$ is still limited, even though being actively sought out to $z\sim6$. Here, we use the high resolution and sensitivity of JWST to measure the host galaxy properties for 107 X-ray-selected type-I AGNs at $0.68

We report the identification of 64 new candidates of compact galaxies, potentially hosting faint quasars with bolometric luminosities of $L_\mathrm{bol} = 10^{43}$--10$^{46}$ erg s$^{-1}$, residing in the reionization epoch within the redshift range of $6 \lesssim z \lesssim 8$. These candidates were selected by harnessing the rich multiband datasets provided by the emerging JWST-driven extragalactic surveys, focusing on COSMOS-Web, as well as JADES, UNCOVER, CEERS, and PRIMER. Our search strategy includes two stages: applying stringent photometric cuts to catalog-level data and detailed spectral energy distribution fitting. These techniques effectively isolate the quasar candidates while mitigating contamination from low-redshift interlopers, such as brown dwarfs and nearby galaxies. The selected candidates indicate physical traits compatible with low-luminosity active galactic nuclei, likely hosting $\approx10^5$--$10^7~M_\odot$ supermassive black holes (SMBHs) living in galaxies with stellar masses of $\approx10^8$--$10^{10}~M_\odot$. The SMBHs selected in this study, on average, exhibit an elevated mass compared to their hosts, with the mass ratio distribution slightly higher than those of galaxies in the local Universe. As with other high-$z$ studies, this is at least in part due to the selection method for these quasars. An extensive Monte Carlo analysis provides compelling evidence that heavy black hole seeds from the direct collapse scenario appear to be the preferred pathway to mature this specific subset of SMBHs by $z\approx7$. This work underscores the significance of further spectroscopic observations, as the quasar candidates presented here offer exceptional opportunities to delve into the nature of the earliest galaxies and SMBHs that formed during cosmic infancy., Comment: Accepted for publication in the Astronomy & Astrophysics journal. 19 pages, 10 figures, and 4 tables. We welcome comments from the reader

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

We present observations and analysis of the starburst, PACS-819, at z=1.45 ($M_*=10^{10.7}$ M$_{ \odot}$), using high-resolution ($0^{\prime \prime}.1$; 0.8 kpc) ALMA and multi-wavelength JWST images from the COSMOS-Web program. Dissimilar to HST/ACS images in the rest-frame UV, the redder NIRCam and MIRI images reveal a smooth central mass concentration and spiral-like features, atypical for such an intense starburst. Through dynamical modeling of the CO J=5--4 emission with ALMA, PACS-819 is rotation-dominated thus has a disk-like nature. However, kinematic anomalies in CO and asymmetric features in the bluer JWST bands (e.g., F150W) support a more disturbed nature likely due to interactions. The JWST imaging further enables us to map the distribution of stellar mass and dust attenuation, thus clarifying the relationships between different structural components, not discernable in the previous HST images. The CO J = 5 -- 4 and FIR dust continuum emission are co-spatial with a heavily-obscured starbursting core (<1 kpc) which is partially surrounded by much less obscured star-forming structures including a prominent arc, possibly a tidally-distorted dwarf galaxy, and a clump, either a sign of an ongoing violent disk instability or a recently accreted low-mass satellite. With spatially-resolved maps, we find a high molecular gas fraction in the central area reaching $\sim3$ ($M_{\text{gas}}$/$M_*$) and short depletion times ($M_{\text{gas}}/SFR\sim$ 120 Myrs) across the entire system. These observations provide insights into the complex nature of starbursts in the distant universe and underscore the wealth of complementary information from high-resolution observations with both ALMA and JWST., Comment: 21 pages(including appendix), 13 figures, accepted by ApJ

We take a deep learning-based approach for galaxy merger identification in Subaru HSC-SSP, specifically through the use of deep representation learning and fine-tuning, with the aim of creating a pure and complete merger sample within the HSC-SSP survey. We can use this merger sample to conduct studies on how mergers affect galaxy evolution. We use Zoobot, a deep learning representation learning model pre-trained on citizen science votes on Galaxy Zoo DeCALS images. We fine-tune Zoobot for the purpose of merger classification of images of SDSS and GAMA galaxies in HSC-SSP PDR 3. Fine-tuning is done using 1200 synthetic HSC-SSP images of galaxies from the TNG simulation. We then find merger probabilities on observed HSC images using the fine-tuned model. Using our merger probabilities, we examine the relationship between merger activity and environment. We find that our fine-tuned model returns an accuracy on the synthetic validation data of 76%. This number is comparable to those of previous studies where convolutional neural networks were trained with simulation images, but with our work requiring a far smaller number of training samples. For our synthetic data, our model is able to achieve completeness and precision values of 80%. In addition, our model is able to correctly classify both mergers and non-mergers of diverse morphologies and structures, including those at various stages and mass ratios, while distinguishing between projections and merger pairs. For the relation between galaxy mergers and environment, we find two distinct trends. Using stellar mass overdensity estimates for TNG simulations and observations using SDSS and GAMA, we find that galaxies with higher merger scores favor lower density environments on scales of 0.5 to 8 h^-1 Mpc. However, below these scales in the simulations, we find that galaxies with higher merger scores favor higher density environments., Comment: 36 pages, 15 figures, accepted to Astronomy and Astrophysics

A key question in galaxy evolution has been the importance of the apparent `clumpiness' of high redshift galaxies. Until now, this property has been primarily investigated in rest-frame UV, limiting our understanding of their relevance. Are they short-lived or are associated with more long-lived massive structures that are part of the underlying stellar disks? We use JWST/NIRCam imaging from CEERS to explore the connection between the presence of these `clumps' in a galaxy and its overall stellar morphology, in a mass-complete ($log\,M_{*}/M_{\odot} > 10.0$) sample of galaxies at $1.0 < z < 2.0$. Exploiting the uninterrupted access to rest-frame optical and near-IR light, we simultaneously map the clumps in galactic disks across our wavelength coverage, along with measuring the distribution of stars among their bulges and disks. Firstly, we find that the clumps are not limited to rest-frame UV and optical, but are also apparent in near-IR with $\sim 60\,\%$ spatial overlap. This rest-frame near-IR detection indicates that clumps would also feature in the stellar-mass distribution of the galaxy. A secondary consequence is that these will hence be expected to increase the dynamical friction within galactic disks leading to gas inflow. We find a strong negative correlation between how clumpy a galaxy is and strength of the bulge. This firmly suggests an evolutionary connection, either through clumps driving bulge growth, or the bulge stabilizing the galaxy against clump formation, or a combination of the two. Finally, we find evidence of this correlation differing from rest-frame optical to near-IR, which could suggest a combination of varying formation modes for the clumps., Comment: Resubmitted to ApJ after addressing minor-comments from referee

At fixed galaxy stellar mass, there is a clear observational connection between structural asymmetry and offset from the star forming main sequence, $\Delta$SFMS. Herein, we use the TNG50 simulation to investigate the relative roles of major mergers (stellar mass ratios $\mu\geq0.25$), minor ($0.1 \leq \mu < 0.25$), and mini mergers ($0.01 \leq \mu < 0.1$) in driving this connection amongst star forming galaxies (SFGs). We use dust radiative transfer post-processing with SKIRT to make a large, public collection of synthetic Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) images of simulated TNG galaxies over $0.1\leq z \leq 0.7$ with $\log M_{\star} / \mathrm{M}_{\odot}\geq9$ ($\sim750$k images). Using their instantaneous SFRs, known merger histories/forecasts, and HSC-SSP asymmetries, we show (1) that TNG50 SFGs qualitatively reproduce the observed trend between $\Delta$SFMS and asymmetry and (2) a strikingly similar trend emerges between $\Delta$SFMS and the time-to-coalescence for mini mergers. Controlling for redshift, stellar mass, environment, and gas fraction, we show that individual mini merger events yield small enhancements in SFRs and asymmetries that are sustained on long timescales (at least $\sim3$ Gyr after coalescence, on average) -- in contrast to major/minor merger remnants which peak at much greater amplitudes but are consistent with controls only $\sim1$ Gyr after coalescence. Integrating the boosts in SFRs and asymmetries driven by $\mu\geq0.01$ mergers since $z=0.7$ in TNG50 SFGs, we show that mini mergers are responsible for (i) $55$ per cent of all merger-driven star formation and (ii) $70$ per cent of merger-driven asymmetric structure. Due to their relative frequency and prolonged boost timescales, mini mergers dominate over their minor and major counterparts in driving star formation and asymmetry in SFGs., Comment: 31 pages; 19 figures; MNRAS accepted; Image data available via the TNG website: www.tng-project.org/bottrell23

In this letter, we report the discovery of the highest redshift, heavily obscured, radio-loud AGN candidate selected using JWST NIRCam/MIRI, mid-IR, sub-mm, and radio imaging in the COSMOS-Web field. Using multi-frequency radio observations and mid-IR photometry, we identify a powerful, radio-loud (RL), growing supermassive black hole (SMBH) with significant spectral steepening of the radio SED ($f_{1.28 \mathrm{GHz}} \sim 2$ mJy, $q_{24\mu m} = -1.1$, $\alpha_{1.28-3\mathrm{GHz}}=-1.2$, $\Delta \alpha = -0.4$). In conjunction with ALMA, deep ground-based observations, ancillary space-based data, and the unprecedented resolution and sensitivity of JWST, we find no evidence of AGN contribution to the UV/optical/NIR data and thus infer heavy amounts of obscuration (N$_{\mathrm{H}} > 10^{23}$ cm$^{-2}$). Using the wealth of deep UV to sub-mm photometric data, we report a singular solution photo-z of $z_\mathrm{phot}$ = 7.7$^{+0.4}_{-0.3}$ and estimate an extremely massive host-galaxy ($\log M_{\star} = 11.4 -12\,\mathrm{M}_{\odot}$) hosting a powerful, growing SMBH (L$_{\mathrm{Bol}} = 4-12 \times 10^{46}$ erg s$^{-1}$). This source represents the furthest known obscured RL AGN candidate, and its level of obscuration aligns with the most representative but observationally scarce population of AGN at these epochs., Comment: Accepted to ApJL

We report the discovery of 15 exceptionally luminous $10\lesssim z\lesssim14$ candidate galaxies discovered in the first 0.28 deg$^2$ of JWST/NIRCam imaging from the COSMOS-Web Survey. These sources span rest-frame UV magnitudes of $-20.5>M_{\rm UV}>-22$, and thus constitute the most intrinsically luminous $z\gtrsim10$ candidates identified by JWST to-date. Selected via NIRCam imaging with Hubble ACS/F814W, deep ground-based observations corroborate their detection and help significantly constrain their photometric redshifts. We analyze their spectral energy distributions using multiple open-source codes and evaluate the probability of low-redshift solutions; we conclude that 12/15 (80%) are likely genuine $z\gtrsim10$ sources and 3/15 (20%) likely low-redshift contaminants. Three of our $z\sim12$ candidates push the limits of early stellar mass assembly: they have estimated stellar masses $\sim5\times10^{9}\,M_\odot$, implying an effective stellar baryon fraction of $\epsilon_{\star}\sim0.2-0.5$, where $\epsilon_{\star}\equiv M_{\star}/(f_{b}M_{halo})$. The assembly of such stellar reservoirs is made possible due to rapid, burst-driven star formation on timescales $<$100\,Myr where the star-formation rate may far outpace the growth of the underlying dark matter halos. This is supported by the similar volume densities inferred for $M_\star\sim10^{10}\,M_\odot$ galaxies relative to $M_\star\sim10^{9}\,M_\odot$ -- both about $10^{-6}$ Mpc$^{-3}$ -- implying they live in halos of comparable mass. At such high redshifts, the duty cycle for starbursts would be of order unity, which could cause the observed change in the shape of the UVLF from a double powerlaw to Schechter at $z\approx8$. Spectroscopic redshift confirmation and ensuing constraints of their masses will be critical to understanding how, and if, such early massive galaxies push the limits of galaxy formation in $\Lambda$CDM., Comment: 30 pages, 9 figures; ApJ submitted

We report the identification of 15 galaxy candidates at $z\ge9$ using the initial COSMOS-Web JWST observations over 77 arcmin$^2$ through four NIRCam filters (F115W, F150W, F277W, F444W) with an overlap with MIRI (F770W) of 8.7 arcmin$^2$. We fit the sample using several publicly-available SED fitting and photometric redshift codes and determine their redshifts between $z=9.3$ and $z=10.9$ ($\langle z\rangle=10.0$), UV-magnitudes between M$_{\rm UV}$ = $-$21.2 and $-$19.5 (with $\langle $M$_{\rm UV}\rangle=-20.2$) and rest-frame UV slopes ($\langle \beta\rangle=-2.4$). These galaxies are, on average, more luminous than most $z\ge9$ candidates discovered by JWST so far in the literature, while exhibiting similar blue colors in their rest-frame UV. The rest-frame UV slopes derived from SED-fitting are blue ($\beta\sim$[$-$2.0, $-$2.7]) without reaching extremely blue values as reported in other recent studies at these redshifts. The blue color is consistent with models that suggest the underlying stellar population is not yet fully enriched in metals like similarly luminous galaxies in the lower redshift Universe. The derived stellar masses with $\langle \log_{\rm 10} ($M$_\star/$M$_\odot)\rangle\approx8-9$ are not in tension with the standard $\Lambda$CDM model and our measurement of the volume density of such UV luminous galaxies aligns well with previously measured values presented in the literature at $z\sim9-10$. Our sample of galaxies, although compact, are significantly resolved., Comment: Submitted to ApJ

A large fraction of the accreting supermassive black hole population is shrouded by copious amounts of gas and dust, particularly in the distant ($z\gtrsim1$) Universe. While much of the obscuration is attributed to a parsec-scale torus, there is a known contribution from the larger-scale host galaxy. Using JWST/NIRCam imaging from the COSMOS-Web survey, we probe the galaxy-wide dust distribution in X-ray selected AGN up to $z\sim2$. Here, we focus on a sample of three AGNs with their host galaxies exhibiting prominent dust lanes, potentially due to their edge-on alignment. These represent 27% (3 out of 11 with early NIRCam data) of the heavily obscured ($N_H>10^{23}$ cm$^{-2}$) AGN population. With limited signs of a central AGN in the optical and near-infrared, the NIRCam images are used to produce reddening maps $E(B-V)$ of the host galaxies. We compare the mean central value of $E(B-V)$ to the X-ray obscuring column density along the line-of-sight to the AGN ($N_H\sim10^{23-23.5}$ cm$^{-2}$). We find that the extinction due to the host galaxy is present ($0.6\lesssim E(B-V) \lesssim 0.9$; $1.9 \lesssim A_V \lesssim 2.8$) and significantly contributes to the X-ray obscuration at a level of $N_H\sim10^{22.5}$ cm$^{-2}$ assuming an SMC gas-to-dust ratio which amounts to $\lesssim$30% of the total obscuring column density. These early results, including three additional cases from CEERS, demonstrate the ability to resolve such dust structures with JWST and separate the different circumnuclear and galaxy-scale obscuring structures., Comment: 12 pages, 8 figures, Accepted for publication in ApJL

We present the quasar luminosity function (LF) at $z = 7$, measured with 35 spectroscopically confirmed quasars at $6.55 < z < 7.15$. The sample of 22 quasars from the Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, combined with 13 brighter quasars in the literature, covers an unprecedentedly wide range of rest-frame ultraviolet magnitudes over $-28 < M_{1450} < -23$. We found that the binned LF flattens significantly toward the faint end populated by the SHELLQs quasars. A maximum likelihood fit to a double power-law model has a break magnitude $M^*_{1450} = -25.60^{+0.40}_{-0.30}$, a characteristic density $\Phi^* = 1.35^{+0.47}_{-0.30}$ Gpc$^{-3}$ mag$^{-1}$, and a bright-end slope $\beta = -3.34^{+0.49}_{-0.57}$, when the faint-end slope is fixed to $\alpha = -1.2$ as observed at $z \le 6$. The overall LF shape remains remarkably similar from $z = 4$ to $7$, while the amplitude decreases substantially toward higher redshifts, with a clear indication of an accelerating decline at $z \ge 6$. The estimated ionizing photon density, $10^{48.2 \pm 0.1}$ s$^{-1}$ Mpc$^{-3}$, is less than 1 % of the critical rate to keep the intergalactic medium ionized at $z = 7$, and thus indicates that quasars are not a major contributor to cosmic reionization., Comment: The Astrophysical Journal Letters, in press

How does the host galaxy morphology influence a central quasar or vice versa? We address this question by measuring the asymmetries of 2424 SDSS quasar hosts at $0.245$). This then manifests itself into quasar hosts being more asymmetric, on average, when they harbor a more massive and highly accreting black hole. The merger fraction also positively correlates with $L_{\rm bol}$ and reaches up to 35\% for the most luminous. Compared to control galaxies, quasar hosts are marginally more asymmetric (excess of 0.017 in median at 9.4$\sigma$ level) and the merger fractions are similar ($\sim 16.5\%$). We quantify the dependence of asymmetry on optical band which demonstrates that mergers are more likely to be identified with the bluer bands and the correlation between $L_{\rm bol}$ and asymmetry is also stronger in such bands. We stress that the band dependence, indicative of a changing stellar population, is an important factor in considering the influence of mergers on AGN activity., Comment: 27 pages, 28 figures

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

We investigate the physical properties, such as star-forming activity, disk vs. bulge nature, galaxy size, and obscuration of 3796 X-ray selected AGNs at $0.2

The detection of starlight from the host galaxies of quasars during the reionization epoch ($z>6$) has been elusive, even with deep HST observations. The current highest redshift quasar host detected, at $z=4.5$, required the magnifying effect of a foreground lensing galaxy. Low-luminosity quasars from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) mitigate the challenge of detecting their underlying, previously-undetected host galaxies. Here we report rest-frame optical images and spectroscopy of two HSC-SSP quasars at $z>6$ with JWST. Using NIRCam imaging at 3.6$\mu$m and 1.5$\mu$m and subtracting the light from the unresolved quasars, we find that the host galaxies are massive (stellar masses of $13\times$ and $3.4\times$ $10^{10}$ M$_{\odot}$, respectively), compact, and disk-like. NIRSpec medium-resolution spectroscopy shows stellar absorption lines in the more massive quasar, confirming the detection of the host. Velocity-broadened gas in the vicinity of these quasars enables measurements of their black hole masses ($1.4\times 10^9$ and $2.0\times$ $10^{8}$ M$_{\odot}$, respectively). Their location in the black hole mass - stellar mass plane is consistent with the distribution at low redshift, suggesting that the relation between black holes and their host galaxies was already in place less than a billion years after the Big Bang., Comment: Matched to the published Nature version of the article. 27 pages, 4 main figures, 1 table, 6 supplementary figures, 2 supplementary table

We present the survey design, implementation, and outlook for COSMOS-Web, a 255 hour treasury program conducted by the James Webb Space Telescope in its first cycle of observations. COSMOS-Web is a contiguous 0.54 deg$^2$ NIRCam imaging survey in four filters (F115W, F150W, F277W, and F444W) that will reach 5$\sigma$ point source depths ranging $\sim$27.5-28.2 magnitudes. In parallel, we will obtain 0.19 deg$^2$ of MIRI imaging in one filter (F770W) reaching 5$\sigma$ point source depths of $\sim$25.3-26.0 magnitudes. COSMOS-Web will build on the rich heritage of multiwavelength observations and data products available in the COSMOS field. The design of COSMOS-Web is motivated by three primary science goals: (1) to discover thousands of galaxies in the Epoch of Reionization ($64$ and place constraints on the formation of the Universe's most massive galaxies ($M_\star>10^{10}$\,M$_\odot$), and (3) directly measure the evolution of the stellar mass to halo mass relation using weak gravitational lensing out to $z\sim2.5$ and measure its variance with galaxies' star formation histories and morphologies. In addition, we anticipate COSMOS-Web's legacy value to reach far beyond these scientific goals, touching many other areas of astrophysics, such as the identification of the first direct collapse black hole candidates, ultracool sub-dwarf stars in the Galactic halo, and possibly the identification of $z>10$ pair-instability supernovae. In this paper we provide an overview of the survey's key measurements, specifications, goals, and prospects for new discovery., Comment: 46 pages, 16 figures, ApJ accepted

The conditions under which galactic nuclear regions become active are largely unknown, although it has been hypothesized that secular processes related to galaxy morphology could play a significant role. We investigate this question using optical i-band images of 3096 SDSS quasars and galaxies at 0.3

We measure the host galaxy properties of five quasars with $z\sim 1.6 - 3.5$ selected from SDSS and AEGIS, which fall within the JWST/HST CEERS survey area. A PSF library is constructed based on stars in the full field-of-view of the data and used with the 2-dimensional image modeling tool galight to decompose the quasar and its host with multi-band filters available for HST ACS+WFC3 and JWST NIRCAM (12 filters covering HST F606W to JWST F444W). As demonstrated, JWST provides the first capability to detect quasar hosts at $z>3$ and enables spatially-resolved studies of the underlying stellar populations at $z\sim2$ within morphological structures (spiral arms, bar) not possible with HST. Overall, we find quasar hosts to be disk-like, lack merger signatures, and have sizes generally more compact than typical star-forming galaxies at their respective stellar mass, thus in agreement with results at lower redshifts. The fortuitous face-on orientation of SDSSJ1420+5300A at $z = 1.646$ enables us to find higher star formation and younger ages in the central $2-4$ kpc region relative to the outskirts, which may help explain the relatively compact nature of quasar hosts and pose a challenge to AGN feedback models., Comment: 12 pages, 1 table, 6 figures; ApJL accepted; Comments welcome

We carry out a comparative analysis of the relation between the mass of supermassive black holes (BHs) and the stellar mass of their host galaxies at $0.2

We assess the dark matter halo masses of luminous AGNs over the redshift range 0.2 to 1.2 using galaxy-galaxy lensing based on imaging data from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). We measure the weak lensing signal of a sample of 48907 AGNs constructed using HSC and WISE photometry. %The lensing detection around AGNs has a signal to noise ratio of \textbf{29}. As expected, we find that the lensing mass profile of total AGN sample is consistent with that of massive galaxies ($\rm log(M_{*}/h^{-2}M_\odot)\sim$ 10.61). Surprisingly, the lensing signal remains unchanged when the AGN sample is split into four stellar mass bins of host galaxies. Specifically, we find that the excess surface density (ESD) of AGNs, residing in galaxies with high stellar masses, significantly differs from that of the control sample. We further fit a halo occupation distribution model to the data to infer the posterior distribution of parameters including the average halo mass. We find that the characteristic halo mass of the full AGN population lies near the knee ($\rm log(M_h/h^{-1}M_{\odot})=12.0$) of the stellar-to-halo mass relation (SHMR). Illustrative of the results given above, the halo masses of AGNs residing in host galaxies with high stellar masses (i.e., above the knee of the SHMR) falls below the calibrated SHMR while the halo mass of the low stellar mass sample is more consistent with the established SHMR. These results indicate that massive halos with higher clustering bias tends to suppress AGN activity, probably due to the lack of available gas., Comment: 13 pages,10 figures, The Astrophysical Journal, Volume 977, Issue 1, id.59, 13 pp

The correlation between the mass of supermassive black holes (SMBHs; $\mathcal{M}_{\rm BH}$) and their host galaxies ($\mathcal{M}_\star$) in the reionization epoch provides valuable constraints on their early growth. High-redshift quasars typically have a $\mathcal{M}_{\rm BH}$/$\mathcal{M}_\star$ ratio significantly elevated in comparison to the local value. However, the degree to which this apparent offset is driven by observational biases is unclear for the most distant quasars. To address this issue, we model the sample selection and measurement biases for a compilation of 20 quasars at $z\sim6$ with host properties based on ALMA observations. We find that the observed distribution of quasars in the $\mathcal{M}_{\rm BH} - \mathcal{M}_\star$ plane can be reproduced by assuming that the underlying SMBH population at $z\sim6$ follows the relationship in the local universe. However, a positive or even a negative evolution in $\mathcal{M}_{\rm BH}$/$\mathcal{M}_\star$ can also explain the data, depending on whether the intrinsic scatter evolves and the strength of various systematic uncertainties. To break these degeneracies, an improvement in the accuracy of mass measurements and an expansion of the current sample to lower $\mathcal{M}_{\rm BH}$ limits are needed. Furthermore, assuming a radiative efficiency of 0.1 and quasar duty cycles estimated from the active SMBH fraction, significant outliers in $\mathcal{M}_{\rm BH}$/$\mathcal{M}_\star$ tend to move toward the local relation given their instantaneous BH growth rate and star formation rate. This may provide evidence for a self-regulated SMBH-galaxy coevolution scenario that is in place at $z\sim6$, with AGN feedback being a possible driver., Comment: 10 pages, 5 figures, published in ApJL

We present the spectroscopic discovery of 69 quasars at 5.8 < z < 7.0, drawn from the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) imaging survey data. This is the 16th publication from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, and completes identification of all but the faintest candidates (i.e., i-band dropouts with zAB < 24 and y-band detections, and z-band dropouts with yAB < 24) with Bayesian quasar probability Pq > 0.1 in the HSC-SSP third public data release (PDR3). The sample reported here also includes three quasars with Pq < 0.1 at z ~ 6.6, which we selected in an effort to completely cover the reddest point sources with simple color cuts. The number of high-z quasars discovered in SHELLQs has now grown to 162, including 23 type-II quasar candidates. This paper also presents identification of seven galaxies at 5.6 < z < 6.7, an [O III] emitter at z = 0.954, and 31 Galactic cool stars and brown dwarfs. High-z quasars and galaxies comprise 75 % and 16 % respectively of all the spectroscopic SHELLQs objects that pass our latest selection algorithm with the PDR3 photometry. That is, a total of 91 % of the objects lie at z > 5.6. This demonstrates that the algorithm has very high efficiency, even though we are probing an unprecedentedly low-luminosity population down to M1450 ~ -21 mag., Comment: Accepted for publication in The Astrophysical Journal Supplement Series

Galight is a Python-based open-source package that can be used to perform two-dimensional model fitting of optical and near-infrared images to characterize the light distribution of galaxies with components including a disk, bulge, bar, and quasar. The decomposition of stellar components has been demonstrated in published studies of inactive galaxies and quasar host galaxies observed by the Hubble Space Telescope and Subaru's Hyper Suprime-Cam. Galight utilizes the image modeling capabilities of lenstronomy while redesigning the user interface for the analysis of large samples of extragalactic sources. The package is user-friendly with some automatic features such as determining the cutout size of the modeling frame, searching for PSF-stars in field-of-view, estimating the noise map of the data, identifying all the objects to set the initial model, and associated parameters to fit them simultaneously. These features minimize the manpower and allow the automatic fitting tasks. The software is distributed under the MIT license. The source code, installation guidelines, and example notebooks code can be found at https://galight.readthedocs.io/en/latest/, Comment: 10 pages, 5 figures. This manuscript is posted as a user manual and will be constantly updated in the future with more upcoming features

We present the galaxy size-mass ($R_{e}-M_{\ast}$) distributions using a stellar-mass complete sample of $\sim1.5$ million galaxies, covering $\sim100$ deg$^2$, with $\log(M_{\ast}/M_{\odot})>10.2~(9.2)$ over the redshift range $0.2

We measure the relationship between stellar mass and stellar metallicity, the stellar mass--metallicity relation (MZR), for 1336 star-forming galaxies at $1.6\le z\le3.0$ (=2.2) using rest-frame far-ultraviolet spectra from the zCOSMOS-deep survey. High signal-to-noise composite spectra containing stellar absorption features are fit with population synthesis model spectra of a range of metallicity. We find stellar metallicities, which mostly reflect iron abundances, scaling as $(Z_{Fe,\ast}/Z_{Fe,\odot})=-(0.81\pm0.01)+(0.32+0.03)\log(M_\ast/10^{10}M_\odot)$ across the mass range of $10^9\lesssim M_\ast/M_\odot\lesssim10^{11}$, being $\approx6\times$ lower than seen locally at the same masses. The instantaneous oxygen-to-iron ratio ($\alpha$-enhancement) inferred using the gas-phase oxygen MZRs, is on average found to be [O/Fe]$\approx0.47$, being higher than the local [O/Fe]$\approx0$. The observed changes in [O/Fe] and [Fe/H] are reproduced in simple flow-through gas-regulator models with steady star-formation histories (SFHs) that follow the evolving main sequence. Our models show that the [O/Fe] is determined almost entirely by the instantaneous specific star formation rate alone while being independent of the SFHs, mass, and the gas-regulation characteristics of the systems. We find that the locations of $\sim10^{10}M_\odot$ galaxies at z~2 in the [O/Fe]--metallicity planes are in remarkable agreement with the sequence of low-metallicity thick-disk stars in our Galaxy. This manifests a beautiful concordance between the results of Galactic archaeology and observations of high-redshift Milky Way progenitors. However, there remains a question of how and when the old metal-rich, low-$\alpha$/Fe stars seen in the bulge had formed by z~2 because such a stellar population is not seen in our data and difficult to explain in the context of our models., Comment: Submitted to ApJ; 30 pages, 20 figures

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

We measure the evolution of the $M_{\rm BH}-M_*$ relation using 584 uniformly-selected SDSS quasars at $0.2