Search

Quasar feedback may play a key role in the evolution of massive galaxies. The dust-reddened quasar, F2M110648.35$+$480712 at $z = 0.4352$ is one of the few cases at its redshift that exhibits powerful quasar feedback through bipolar outflows. Our new observation with the integral field unit mode of Near-infrared Spectrograph onboard JWST opens a new window to examine this spectacular outflow through Pa$\alpha$ emission line with $\sim$3$\times$ better spatial resolution than previous work. The morphology and kinematics of the Pa$\alpha$ nebula confirm the existence of a bipolar outflow extending on a scale of $\sim$17$\times$14 kpc and with a velocity reaching $\sim$1100 km s$^{-1}$. The higher spatial resolution of our new observation leads to more reliable measurements of outflow kinematics. Considering only the spatially resolved outflow and assuming an electron density of 100 cm$^{-2}$, the mass, momentum and kinetic energy outflow rates are $\sim$50-210 M$_{\odot}$ yr$^{-1}$, $\sim$0.3-1.7$\times$10$^{36}$ dynes ($\sim$14-78\% of the quasar photon momentum flux) and $\sim$0.16-1.27$\times$10$^{44}$ erg s$^{-1}$ ($\sim$0.02-0.20\% of the quasar bolometric luminosity), respectively. The local instantaneous outflow rates generally decrease radially. We infer that the quasar is powerful enough to drive the outflow, while stellar processes cannot be overlooked as a contributing energy source. The mass outflow rate is $\sim$0.4-1.5 times the star formation rate, and the ratio of kinetic energy outflow rate to the quasar bolometric luminosity is comparable to the minimum value required for negative quasar feedback in simulations. This outflow may help regulate the star formation activity within the system to some extent., Comment: 14 pages, 6 figures, 1 table, ApJ in review

JWST has revealed an abundance of compact, red objects at $z\approx5-8$ dubbed "little red dots" (LRDs), whose SEDs display a faint blue UV continuum followed by a steep rise in the optical. Despite extensive study of their characteristic V-shaped SEDs, the nature of LRDs remains unknown. We present a new analysis of the NIRSpec/PRISM spectrum of A2744-QSO1, a triply imaged LRD at $z=7.04$ from the UNCOVER survey. The spectrum shows a strong Balmer break and broad Balmer emission lines, both of which are difficult to explain with models invoking exclusively AGN or stellar contributions. Our fiducial model decomposes the spectrum into a post-starburst galaxy dominating the UV-optical continuum and a reddened AGN being sub-dominant at all wavelength and contributing at $\sim20\%$ level. However, our most credible model infers a stellar mass of $M_\star\approx 4\times10^9\,\mathrm{M_\odot}$ within a radius of $r_\mathrm{e}<30\,$pc, driving its central density to the highest among observations to date. This high central density could be explained if A2744-QSO-1 is the early-forming core of a modern-day massive elliptical galaxy that later puffed up via the inside-out growth channel. The models also necessitate an unusually steep dust law to preserve the strong break strength, though this steepness may be explained by a deficit of large dust grains. It is also probable that these challenges reflect our ignorance of A2744-QSO1's true nature. Future variability and reverberation mapping studies could help disentangle the galaxy and AGN contribution to the continuum, and deeper redder observations could also unveil the dust properties in LRDs., Comment: 19 pages, 7 figures, 2 tables. Submitted to ApJ. Comments are welcome

JWST has revealed diverse new populations of high-redshift ($z\sim4-11$) AGN and extreme star-forming galaxies that challenge current models. In this paper, we use rest-frame UV emission-line diagnostics to identify AGN candidates and other exceptional ionizing sources, complementing previous studies predominantly focused on broad-line AGN. In this paper, we use rest-frame UV emission-line diagnostics to identify AGN candidates and other exceptional ionizing sources, complementing previous studies predominantly focused on broad-line AGN. From a parent sample of 205 $\mathrm{z_{spec}}>3$ UNCOVER galaxies with NIRSpec/PRISM follow-up, we identify 12 C IV, He II, and C III] emitters. Leveraging the combined rest-optical and UV coverage of PRISM, we limit the emission-line model space using the sample's [O III]/H$\beta$ distribution, significantly decreasing the overlap between AGN and star-formation models in the UV diagnostics. We then find that the five He II emitters are the strongest AGN candidates, with further support from two [Ne V] detections and one X-ray detection from Chandra. We cannot robustly quantify the AGN fraction in this sample, but we note that close to 20% of $\mathrm{M_{*}>2\times10^{9}\,M_{\odot}}$ parent sample galaxies are AGN candidates. The lower-mass line emitters, which are consistent with both AGN and star-forming photoionization models, have more compact sizes and higher specific star formation rates than the parent sample. Higher-resolution and deeper data on these UV line emitters should provide much stronger constraints on the obscured AGN fraction at $z > 3$., Comment: Submitted to ApJ

Active galactic nuclei (AGN) are the signposts of black hole growth, and likely play an important role in galaxy evolution. An outstanding question is whether AGN of different spectral types indicate different evolutionary stages in the coevolution of black holes and galaxies. We present the angular correlation function between an AGN sample selected from the Hyper Suprime Camera Subaru Strategic Program (HSC-SSP) optical + Wide-field Infrared Survey Explorer (WISE) mid-IR photometry, and a luminous red galaxy (LRG) sample from HSC-SSP. We investigate AGN clustering strength as a function of their luminosity and spectral features across three independent HSC fields totaling $\sim600\,{\rm deg^{2}}$, for $z\in0.6-1.2$ and AGN with $L_{6\mu m}>3\times10^{44}{\rm\,erg\,s^{-1}}$. There are $\sim28,500$ AGN and $\sim1.5$ million LRGs in our primary analysis. We determine the inferred average halo mass for the full AGN sample ($M_h \approx 10^{12.9}h^{-1}M_\odot$), and note that it does not evolve significantly as a function of redshift (over this narrow range) or luminosity. We find that, on average, unobscured AGN ($M_h \approx10^{13.3}h^{-1}M_\odot$) occupy $\sim4.5\times$ more massive halos than obscured AGN ($M_h \approx10^{12.6}h^{-1}M_\odot$), at $5\sigma$ statistical significance using 1-D uncertainties, and at $3\sigma$ using the full covariance matrix, suggesting a physical difference between unobscured and obscured AGN, beyond the line-of-sight viewing angle. Furthermore, we find evidence for a halo mass dependence on reddening level within the Type I AGN population, which could support the existence of a previously claimed dust-obscured phase in AGN-host galaxy coevolution. However, we also find that even quite small systematic shifts in the redshift distributions of the AGN sample could plausibly explain current and previously observed differences in $M_{h}$., Comment: 31 pages, 14 figures, submitted to ApJ

We present the design and observations of low resolution JWST/NIRSpec PRISM spectroscopy from the Ultradeep NIRSpec and NIRCam ObserVations before the Epoch of Reionization (UNCOVER) Cycle 1 JWST Treasury program. Targets are selected using JWST/NIRCam photometry from UNCOVER and other programs, and cover a wide range of categories and redshifts to ensure the legacy value of the survey. These categories include the first galaxies at $z\gtrsim10$, faint galaxies during the Epoch of Reionization ($z\sim6-8$), high redshift AGN ($z\gtrsim6$), Population III star candidates, distant quiescent and dusty galaxies ($1\lesssim z \lesssim 6$), and filler galaxies sampling redshift--color--magnitude space from $z\sim 0.1-13$. Seven NIRSpec MSA masks across the extended Abell 2744 cluster were observed, along with NIRCam parallel imaging in 8 filters (F090W, F115W, F150W, F200W, F277W, F356W, F410M, F444W, F480M) over a total area of ~26 arcmin$^2$, overlapping existing HST coverage from programs including the Hubble Frontier Fields and BUFFALO. We successfully observed 553 objects down to $m_{\mathrm{F444W}}\sim30\mathrm{AB}$, and by leveraging mask overlaps, we reach total on-target exposure times ranging from 2.4-16.7h. We demonstrate the success rate and distribution of confirmed redshifts, and also highlight the rich information revealed by these ultradeep spectra for a subset of our targets. An updated lens model of Abell 2744 is also presented, including 14 additional spectroscopic redshifts and finding a total cluster mass of $M_{\mathrm{SL}}=(2.1\pm0.3)\times10^{15}\,\mathrm{M}_{\odot}$. We publicly release reduced 1D and 2D spectra for all objects observed in Summer 2023 along with a spectroscopic redshift catalog and the updated lens model of the cluster (https://jwst-uncover.github.io/DR4.html)., Comment: 19 pages, 9 figures, 4 tables, submitted to ApJ, comments welcome! Data available at: https://jwst-uncover.github.io/DR4.html (v2: figure format correction)

Supermassive black holes (SMBHs) are found in the centers of massive galaxies, and galaxy mergers should eventually lead to SMBH mergers. Quasar activity has long been associated with galaxy mergers, so here we investigate if supermassive black hole binaries (SMBHBs) are preferentially found in quasars. Our multimessenger investigation folds together a gravitational wave background signal from NANOGrav, a sample of periodic AGN candidates from the Catalina Real-Time Transient Survey, and a quasar mass function, to estimate an upper limit on the fraction of quasars which could host a SMBHB. We find at 95\% confidence that quasars are at most seven times as likely to host a SMBHB as a random galaxy. Quasars should therefore be prioritized as targets for SMBHB searches in pulsar timing arrays., Comment: 9 pages, 5 figures, 1 appendix

In this paper, we describe the "Medium Bands, Mega Science" JWST Cycle 2 survey (JWST-GO-4111) and demonstrate the power of these data to reveal both the spatially-integrated and spatially-resolved properties of galaxies from the local universe to the era of cosmic dawn. Executed in November 2023, MegaScience obtained ~30 arcmin^2 of deep multiband NIRCam imaging centered on the z~0.3 Abell 2744 cluster, including eleven medium-band filters and the two shortest-wavelength broad-band filters, F070W and F090W. Together, MegaScience and the UNCOVER Cycle 1 treasury program provide a complete set of deep (~28-30 mag) images in all NIRCam medium- and broad-band filters. This unique dataset allows us to precisely constrain photometric redshifts, map stellar populations and dust attenuation for large samples of distant galaxies, and examine the connection between galaxy structures and formation histories. MegaScience also includes ~17 arcmin^2 of NIRISS parallel imaging in two broad-band and four medium-band filters from 0.9-4.8um, expanding the footprint where robust spectral energy distribution (SED) fitting is possible. We provide example SEDs and multi-band cutouts at a variety of redshifts, and use a catalog of JWST spectroscopic redshifts to show that MegaScience improves both the scatter and catastrophic outlier rate of photometric redshifts by factors of 2-3. Additionally, we demonstrate the spatially-resolved science enabled by MegaScience by presenting maps of the [OIII] line emission and continuum emission in three spectroscopically-confirmed z>6 galaxies. We show that line emission in reionization-era galaxies can be clumpy, extended, and spatially offset from continuum emission, implying that galaxy assembly histories are complex even at these early epochs. We publicly release fully reduced mosaics and photometric catalogs for both the NIRCam primary and NIRISS parallel fields., Comment: 21 pages, 10 figures. Fully reduced imaging, photometric catalogs, and photometric redshift fits publicly available at https://jwst-uncover.github.io/megascience/

The current and upcoming large data volume galaxy surveys require the use of machine learning techniques to maximize their scientific return. This study explores the use of Self-Organizing Maps (SOMs) to estimate galaxy parameters with a focus on handling cases of missing data and providing realistic probability distribution functions for the parameters. We train a SOM with a simulated mass-limited lightcone assuming a ugrizYJHKs+IRAC dataset, mimicking the Hyper Suprime-Cam (HSC) Deep joint dataset. For parameter estimation, we derive SOM likelihood surfaces considering photometric errors to derive total (statistical and systematic) uncertainties. We explore the effects of missing data including which bands are particular critical to the accuracy of the derived parameters. We demonstrate that the parameter recovery is significantly better when the missing bands are "filled-in" rather than if they are completely omitted. We propose a practical method for such recovery of missing data., Comment: 22 pages, 12 figures, 1 table, accepted for publication in The Astronomical Journal

The JWST discovery of ``little red dots'' (LRDs) is reshaping our picture of the early Universe, yet the physical mechanisms driving their compact size and UV-optical colors remain elusive. Here we report an unusually bright LRD ($z=3.1$) observed as part of the RUBIES program. This LRD exhibits broad emission lines (FWHM $\sim4000$km/s), a blue UV continuum, a clear Balmer break and a red continuum sampled out to rest 4 $\mu$m with MIRI. We develop a new joint galaxy and AGN model within the Prospector Bayesian inference framework and perform spectrophotometric modeling using NIRCam, MIRI, and NIRSpec/Prism observations. Our fiducial model reveals a $M_*\sim 10^9M_\odot$ galaxy alongside a dust-reddened AGN driving the optical emission. Explaining the rest-frame optical color as a reddened AGN requires $A_{\rm v}\gtrsim4$, suggesting that a great majority of the accretion disk energy is re-radiated as dust emission. Yet despite clear AGN signatures, we find a surprising lack of hot torus emission, which implies that either the dust emission in this object must be cold, or the red continuum must instead be driven by a massive, evolved stellar population of the host galaxy -- seemingly inconsistent with the high EW broad lines (H$\alpha$ EW $\sim800$\AA). The widths and luminosities of Pa$\beta$, Pa$\delta$, Pa$\gamma$, and H$\alpha$ imply a modest black hole mass of $M_{\rm BH}\sim10^8M_\odot$. Additionally, we identify a narrow blue-shifted HeI absorption in G395M spectra, signaling an ionized outflow with kinetic energy up to $\sim1$\% the luminosity of the AGN. The low redshift of RUBIES-BLAGN-1 combined with the depth and richness of the JWST imaging and spectroscopic observations provide a unique opportunity to build a physical model for these so-far mysterious LRDs, which may prove to be a crucial phase in the early formation of massive galaxies and their supermassive black holes., Comment: 20 pages, 8 figures, 4 tables

We report the spectroscopic confirmation of a massive ($\log(M_\star/M_\odot)=10.34 \pm_{0.07}^{0.06}$), HST-dark ($m_\mathrm{F150W} - m_\mathrm{F444W} = 3.6$) quiescent galaxy at $z_{spec}=3.97$ in the UNCOVER survey. NIRSpec/PRISM spectroscopy and a non-detection in deep ALMA imaging surprisingly reveals that the galaxy is consistent with a low ($<$10 $M_\odot \ \mathrm{yr^{-1}}$) star formation rate despite evidence for moderate dust attenuation. The F444W image is well modeled with a two component \sersic fit that favors a compact, $r_e\sim200$ pc, $n\sim2.9$ component and a more extended, $r_e\sim1.6$ kpc, $n\sim1.7$ component. The galaxy exhibits strong color gradients: the inner regions are significantly redder than the outskirts. Spectral energy distribution models that reproduce both the red colors and low star formation rate in the center of UNCOVER 18407 require both significant ($A_v\sim1.4$ mag) dust attenuation and a stellar mass-weighted age of 900 Myr, implying 50\% of the stars in the core already formed by $z=7.5$. Using spatially resolved annular mass-to-light measurements enabled by the galaxy's moderate magnification ($\mu=2.12\pm_{0.01}^{0.05}$) to reconstruct a radial mass profile from the best-fitting two-component \sersic model, we infer a total mass-weighted $r_\mathrm{eff} = 0.72 \pm_{0.11}^{0.15}$ kpc and log$(\Sigma_\mathrm{1 kpc} \ [\mathrm{M_\odot/kpc^2}]) = 9.61 \pm_{0.10}^{0.08}$. The early formation of a dense, low star formation rate, and dusty core embedded in a less attenuated stellar envelope suggests an evolutionary link between the earliest-forming massive galaxies and their elliptical descendants. Furthermore, the disparity between the global, integrated dust properties and the spatially resolved gradients highlights the importance of accounting for radially varying stellar populations when characterizing the early growth of galaxy structure., Comment: 17 pages, 9 figures, 2 tables. Resubmitted to ApJ after response to referee and update to include new medium band imaging from the JWST MEGASCIENCE program. Comments welcome!

We present the result of a spectroscopic campaign targeting Active Galactic Nucleus (AGN) candidates selected using a novel unsupervised machine-learning (ML) algorithm trained on optical and mid-infrared (mid-IR) photometry. AGN candidates are chosen without incorporating prior AGN selection criteria and are fainter, redder, and more numerous, $\sim$340 AGN deg$^{-2}$, than comparable photometric and spectroscopic samples. In this work we obtain 178 rest-optical spectra from two candidate ML-identified AGN classes with the Hectospec spectrograph on the MMT Observatory. We find that our first ML-identified group, is dominated by Type I AGNs (85%) with a $<3$% contamination rate from non-AGNs. Our second ML-identified group is comprised mostly of Type II AGNs (65%) with a moderate contamination rate of 15% primarily from star-forming galaxies. Our spectroscopic analyses suggest that the classes recover more obscured AGNs, confirming that ML techniques are effective at recovering large populations of AGNs at high levels of extinction. We demonstrate the efficacy of pairing existing WISE data with large-area and deep optical/near-infrared photometric surveys to select large populations of AGNs and recover obscured SMBH growth. This approach is well suited to upcoming photometric surveys, such as Euclid, Rubin, and Roman., Comment: 23 pages, 16 figures. Accepted for publication in The Astronomical Journal

Drawing from the Chandra archive and using a carefully selected set of nearby dwarf galaxies, we present a calibrated high-mass X-ray binary (HMXB) luminosity function in the low-mass galaxy regime and search for an already hinted at dependence on metallicity. Our study introduces a new sample of local dwarf galaxies (D < 12.5 Mpc and M* < 5 x 10^9 M_sun), expanding the specific star-formation rates (sSFR) and gas-phase metallicities probed in previous investigations. Our analysis of the observed X-ray luminosity function indicates a shallower power-law slope for the dwarf galaxy HMXB population. In our study, we focus on dwarf galaxies that are more representative in terms of sSFR compared to prior work. In this regime, the HMXB luminosity function exhibits significant stochastic sampling at high luminosities. This likely accounts for the pronounced scatter observed in the galaxy-integrated HMXB population's Lx/SFR versus metallicity for our galaxy sample. Our calibration is necessary to understand the AGN content of low mass galaxies identified in current and future X-ray survey fields and has implications for binary population synthesis models, as well as X-ray driven cosmic heating in the early universe., Comment: 19 pages, 7 figures

With its sensitivity in the rest-frame optical, the James Webb Space Telescope (JWST) has uncovered active galactic nuclei (AGN), comprising both intrinsically faint and heavily reddened sources, well into the first billion years of the Universe, at $z \sim 4-11$. In this work, we revisit the AGN contribution to reionization given the high number densities associated with these objects. We use the DELPHI semi-analytic model, base-lined against the latest high-redshift datasets from the JWST and the Atacama Large millimetre Array (ALMA) to model early star forming galaxies and AGN. We calculate the escape fractions of ionizing radiation from both star formation and AGN and include the impact of reionization feeback in suppressing the baryonic content of low-mass galaxies in ionized regions. This model is validated against the key observables for star forming galaxy, AGN and reionization. In our {\it fiducial} model, reionization reaches its mid-point at $z \sim 6.9$ and ends by $z \sim 5.9$. Low stellar mass ($M_*\leq 10^9M_\odot$) star forming galaxies are found to be the key drivers of the reionization process, providing about $77\%$ of the total photon budget. Despite their high numbers, high accretion rates and higher escape fractions compared to star forming galaxies at $z \sim 5$, AGN only provide about $23\%$ of the total reionization budget which is dominated by black holes in high stellar mass systems (with $M_* \geq 10^9M_\odot$). This is because AGN number densities become relevant only at $z \leq 7$ - as a result, AGN contribute as much as galaxies as late as $z \sim 6.2$, when reionization is already in its end stages. Finally, we find that even contrasting models of the AGN ionizing photon escape fraction (increasing or decreasing with stellar mass) do not qualitatively change our results., Comment: Submitted to A&A; comments welcome

Over the past year, JWST has uncovered galaxies at record-breaking distances up to $z \sim 13$. The JWST UNCOVER (ultra-deep NIRSpec and NIRcam observations before the epoch of reionization) program has obtained ultra-deep multiwavelength NIRCam imaging of the massive galaxy cluster Abell 2744 over $\sim 45$ arcmin$^{2}$ down to $\sim 29.5$ AB mag. Here, we present a robust ultraviolet (UV) luminosity function derived through lensing clusters at $99$ in recent JWST studies. In particular, a double power-law (DPL) describes better the bright-end of the luminosity function compared to the classical Schechter form. The number density of these bright galaxies is 10-100 times larger than theoretical predictions and previous findings based on Hubble Space Telescope (HST) observations. Additionally, we measure a star formation rate density of $\rho_{\rm SFR} = 10^{-2.64}$ M$_{\odot}$~yr$^{-1}$~Mpc$^{-3}$ at these redshifts, which is 4 to 10 times higher than galaxy formation models that assume a constant star formation efficiency. Future wide-area surveys and accurate modeling of lensing-assisted observations will reliably constrain both the bright and the dim end of the UV luminosity function at $z>9$, which will provide key benchmarks for galaxy formation models., Comment: Published in MNRAS

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

JWST is revealing a new population of dust-reddened broad-line active galactic nuclei (AGN) at redshifts $z\gtrsim5$. Here we present deep NIRSpec/Prism spectroscopy from the Cycle 1 Treasury program UNCOVER of 15 AGN candidates selected to be compact, with red continua in the rest-frame optical but with blue slopes in the UV. From NIRCam photometry alone, they could have been dominated by dusty star formation or AGN. Here we show that the majority of the compact red sources in UNCOVER are dust-reddened AGN: $60\%$ show definitive evidence for broad-line H$\alpha$ with FWHM$\, >2000$ km/s, for $20\%$ current data are inconclusive, and $20\%$ are brown dwarf stars. We propose an updated photometric criterion to select red $z>5$ AGN that excludes brown dwarfs and is expected to yield $>80\%$ AGN. Remarkably, among all $z_{\rm phot}>5$ galaxies with F277W$-$F444W$>1$ in UNCOVER at least $33\%$ are AGN regardless of compactness, climbing to at least $80\%$ AGN for sources with F277W$-$F444W$>1.6$. The confirmed AGN have black hole masses of $10^7-10^9$ M$_{\odot}$. While their UV-luminosities ($-16>M_{\rm UV}>-20$ AB mag) are low compared to UV-selected AGN at these epochs, consistent with percent-level scattered AGN light or low levels of unobscured star formation, the inferred bolometric luminosities are typical of $10^7-10^9$ M$_{\odot}$ black holes radiating at $\sim 10-40\%$ of Eddington. The number densities are surprisingly high at $\sim10^{-5}$ Mpc$^{-3}$ mag$^{-1}$, 100 times more common than the faintest UV-selected quasars, while accounting for $\sim1\%$ of the UV-selected galaxies. While their UV-faintness suggest they may not contribute strongly to reionization, their ubiquity poses challenges to models of black hole growth., Comment: 23 pages, 9 figures, 5 tables, submitted to ApJ

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

Deep observations with JWST have revealed an emerging population of red point-like sources that could provide a link between the postulated supermassive black hole seeds and observed quasars. In this work we present a JWST/NIRSpec spectrum from the JWST Cycle 1 UNCOVER Treasury survey, of a massive accreting black hole at $z=8.50$, displaying a clear broad-line component as inferred from the H$\beta$ line with FWHM = $3439\pm413$ km s$^{-1}$, typical of the broad line region of an active galactic nucleus (AGN). The AGN nature of this object is further supported by high ionization, as inferred from emission lines, and a point-source morphology. We compute the black hole mass of log$_{10}(M_{\rm BH}/M_\odot)=8.17\pm0.42$, and a bolometric luminosity of $L_{\rm bol}\sim6.6\times10^{45}$ erg s$^{-1}$. These values imply that our object is accreting at $\sim 40\%$ of the Eddington limit. Detailed modeling of the spectral energy distribution in the optical and near-infrared, together with constraints from ALMA, indicate an upper limit on the stellar mass of log$_{10}(M_{\rm *}/M_\odot)<8.7$, which would lead to an unprecedented ratio of black hole to host mass of at least $\sim 30 \%$. This is orders of magnitude higher compared to the local QSOs, but is consistent with recent AGN studies at high redshift with JWST. This finding suggests that a non-negligible fraction of supermassive black holes either started out from massive seeds and/or grew at a super-Eddington rate at high redshift. Given the predicted number densities of high-$z$ faint AGN, future NIRSpec observations of larger samples will allow us to further investigate the galaxy-black hole co-evolution in the early Universe., Comment: 14 pages, 6 figures, 2 tables. Accepted for a publication in ApJL

We present JWST NIRSpec prism spectroscopy of gravitationally lensed galaxies at $z\gtrsim9$ found behind the massive galaxy cluster Abell 2744 in the UNCOVER Cycle 1 Treasury Program. We confirm the source redshift via emission lines and/or the Ly$\alpha$ break feature for ten galaxies at z=8.50-13.08 down to $M_{\rm UV}=-17.3$. We achieve a high confirmation rate of 100\% for $z>9$ candidates reported in Atek et al. (2023). Using six sources with multiple emission line detections, we find that the offset of the redshift estimates between the lines and the Ly$\alpha$ break alone with prism can be as large as $\pm0.2$, raising caution in designing future follow-up spectroscopy for the break-only sources. With spec-$z$ confirmed sources in UNCOVER and the literature, we derive lower limits on the rest-frame ultraviolet (UV) luminosity function (LF) at $z\simeq9$-12 and find these lower limits to be consistent with recent photometric measurements. We identify at least two unambiguous and several possible active galactic nucleus (AGN) systems based on X-ray emission, broad line (BL) H$\beta$, high ionization line (e.g., NIV]1487, CIV1549) detections, and excess in UVLF. This requires the AGN LFs at $z\simeq$ 9-10 to be comparable or even higher than the X-ray AGN LF estimated at $z\sim6$ and indicates a plausible cause of the high abundance of $z>9$ galaxies claimed in recent photometric studies may be AGNs. One UV-luminous source is confirmed at the same redshift as a dusty BL AGN at $z=8.50$ with a physical separation of 380 kpc in the source plane. These two sources show blueward Ly$\alpha$ line or continuum emission, suggesting that they reside in the same ionized bubble with a radius of $7.69\pm0.18$ pMpc. Our results imply that AGNs have a non-negligible contribution to cosmic reionization., Comment: 27 pages, 11 figures, 4 tables, submitted to ApJ (See also arXiv:2308.11610)

Early JWST observations have uncovered a new population of red sources that might represent a previously overlooked phase of supermassive black hole growth (Kocevski et al. 2023; Matthee et al. 2023, Labb\'e et al. 2023). One of the most intriguing examples is an extremely red, point-like object that was found to be triply-imaged by the strong lensing (SL) cluster Abell 2744 (Furtak et al. 2023). Here we present deep JWST/NIRSpec observations of this object, Abell2744-QSO1. The spectroscopy confirms that the three images are of the same object, and that it is a highly reddened ($A_V\simeq3$) broad emission-line Active Galactic Nucleus (AGN) at a redshift of $z_{\mathrm{spec}}=7.0451\pm0.0005$. From the width of H$\beta$ ($\mathrm{FWHM}=2800\pm250\,\frac{\mathrm{km}}{\mathrm{s}}$) we derive a black hole mass of $M_{\mathrm{BH}}=4_{-1}^{+2}\times10^7\,\mathrm{M}_{\odot}$. We infer a very high ratio of black hole to galaxy mass of at least 3%, an order of magnitude more than is seen in local galaxies (Bennert et al. 2011), and possibly as high as 100%. The lack of strong metal lines in the spectrum together with the high bolometric luminosity ($L_{\mathrm{bol}}=(1.1\pm0.3)\times10^{45}\,\frac{\mathrm{erg}}{\mathrm{s}}$) indicate that we are seeing the black hole in a phase of rapid growth, accreting at 30% of the Eddington limit. The rapid growth and high black hole to galaxy mass ratio of A2744-QSO1 suggest that it may represent the missing link between black hole seeds (Volonteri et al. 2021) and the first luminous quasars (Fan et al. 2022)., Comment: Published in Nature. Updated to the accepted version

Observations of high-redshift galaxies provide a critical direct test to the theories of early galaxy formation, yet to date, only three have been spectroscopically confirmed at $z>12$. Due to strong gravitational lensing over a wide area, the galaxy cluster field A2744 is ideal for searching for the earliest galaxies. Here we present JWST/NIRSpec observations of two galaxies: a robust detection at $z_{\rm spec} = 12.393^{+0.004}_{-0.001}$, and a plausible candidate at $z_{\rm spec} = 13.079^{+0.013}_{-0.001}$. The galaxies are discovered in JWST/NIRCam imaging and their distances are inferred with JWST/NIRSpec spectroscopy, all from the JWST Cycle 1 UNCOVER Treasury survey. Detailed stellar population modeling using JWST NIRCam and NIRSpec data corroborates the primeval characteristics of these galaxies: low mass ($\sim 10^8~{\rm M_\odot}$), young, rapidly-assembling, metal-poor, and star-forming. Interestingly, both galaxies are spatially resolved, having lensing-corrected rest-UV effective radii on the order of 300-400 pc, which are notably larger than other spectroscopically confirmed systems at similar redshifts. The observed dynamic range of $z \gtrsim 10$ sizes spans over 1 order of magnitude, implying a significant scatter in the size-mass relation at early times. Deep into the epoch of reionization, these discoveries elucidate the emergence of the first galaxies., Comment: ApJL in press; 16 pages, 6 figures, 2 tables

The James Webb Space Telescope is now detecting early black holes (BHs) as they transition from "seeds" to supermassive BHs. Recently Bogdan et al. (2023) reported the detection of an X-ray luminous supermassive BH, UHZ-1, with a photometric redshift at $z > 10$. Such an extreme source at this very high redshift provides new insights on seeding and growth models for BHs given the short time available for formation and growth. Harnessing the exquisite sensitivity of JWST/NIRSpec, here we report the spectroscopic confirmation of UHZ-1 at $z = 10.073 \pm 0.002$. We find that the NIRSpec/Prism spectrum is typical of recently discovered z~10 galaxies, characterized primarily by star-formation features. We see no clear evidence of the powerful X-ray source in the rest-frame UV/optical spectrum, which may suggest heavy obscuration of the central BH, in line with the Compton-thick column density measured in the X-rays. We perform a stellar population fit simultaneously to the new NIRSpec spectroscopy and previously available photometry. The fit yields a stellar mass estimate for the host galaxy that is significantly better constrained than prior photometric estimates ($M_*\sim 1.4^{+0.3}_{-0.4} \times 10^8 M_\odot$). Given the predicted BH mass ($M_{\rm BH}\sim10^7-10^8 M_\odot$), the resulting ratio of $M_{\rm BH}/M_*$ remains two to three orders of magnitude higher than local values, thus lending support to the heavy seeding channel for the formation of supermassive BHs within the first billion years of cosmic evolution., Comment: 10 pages, 4 figures, ApJL in-press. Updated to accepted version incl. X-ray image

The recent Chandra-JWST discovery of a quasar in the z = 10.1 galaxy UHZ1 reveals that accreting supermassive black holes (SMBHs) were already in place 470 million years after the Big Bang. The Chandra X-ray source detected in UHZ1 is a Compton-thick quasar with a bolometric luminosity of $L_{\rm bol}\sim5\times10^{45}\ \rm{erg\ s^{-1}},$ which corresponds to an estimated BH mass of $\sim4\times10^7 \ \rm{M_{\odot}}$ assuming accretion at the Eddington rate. JWST NIRCAM and NIRSpec data yield a stellar mass estimate for UHZ1 comparable to its BH mass. These characteristics are in excellent agreement with prior theoretical predictions for a unique class of transient, high-redshift objects, Over-massive Black Hole Galaxies [OBGs] by Natarajan et al. that harbor a heavy initial black hole seed that likely formed from the direct collapse of the gas. Based on the excellent agreement between the observed multi-wavelength properties of UHZ1 with theoretical model template predictions, suggests that UHZ1 is the first detected OBG candidate. Our assertion rests on multiple lines of concordant evidence between model predictions and the following observed properties of UHZ1: its X-ray detection and the estimated ratio of the X-ray flux to the IR flux that is consistent with theoretical expectations for a heavy initial BH seed; its high measured redshift of z = 10.1, as predicted for the transient OBG stage (9 < z< 12); the amplitude and shape of the detected JWST Spectral Energy Distribution (SED) between 1 - 5 microns, which is in very good agreement with simulated template SEDs for OBGs; and the extended JWST morphology of UHZ1 that is suggestive of a recent merger, also expected for the formation of transient OBGs. As the first OBG candidate, UHZ1 provides compelling evidence for the formation of heavy initial seeds from direct collapse in the early Universe., Comment: 8 pages, 4 figures, submitted to ApJ, under review

We study the broad-band spectral energy distribution (SED) of the prototypical low-mass active galactic nucleus (AGN) in NGC 4395. We jointly model the optical through mid-infrared SED with a combination of galaxy and AGN light, and find that on arcsecond scales, the AGN dominates at most wavelengths. However, there is still some ambiguity about emission from the galaxy, owing partially to the strong short-term variability of the black hole. We investigate the use of smooth and clumpy-torus models in order to disentangle the nuclear infrared emission, as well as exploring the use of poloidal wind emission to account for the blue spectral slope observed in the near-IR. Even when simultaneously fitting the full optical-IR spectral range, we find that degeneracies still remain in the best-fit models. We conclude that high spatial resolution and wider wavelength coverage with the James Webb Space Telescope is needed to understand the mid-infrared emission in this complex highly-variable object, which is the best nearby example to provide a blueprint to finding other low-mass AGN via their mid-infrared emission in the future., Comment: 17 pages, 8 figures

Quasar-driven galactic outflows are a major driver of the evolution of massive galaxies. We report observations of a powerful galactic-scale outflow in a $z=3$ extremely red, intrinsically luminous ($L_{\rm bol}\simeq 5\times 10^{47}$erg s$^{-1}$) quasar SDSSJ1652+1728 with the Near Infrared Spectrograph (NIRSpec) on board JWST. We analyze the kinematics of rest-frame optical emission lines and identify the quasar-driven outflow extending out to $\sim 10$ kpc from the quasar with a velocity offset of ($v_{r}=\pm 500$ km s$^{-1}$) and high velocity dispersion (FWHM$=700-2400$ km s$^{-1}$). Due to JWST's unprecedented surface brightness sensitivity in the near-infrared -- we unambiguously show that the powerful high velocity outflow in an extremely red quasar (ERQ) encompasses a large swath of the host galaxy's interstellar medium (ISM). Using the kinematics and dynamics of optical emission lines, we estimate the mass outflow rate -- in the warm ionized phase alone -- to be at least $2300\pm1400$ $M_{\odot}$ yr$^{-1}$. We measure a momentum flux ratio between the outflow and the quasar accretion disk of $\sim$1 on kpc scale, indicating that the outflow was likely driven in a relatively high ($>10^{23}$cm$^{-2}$) column density environment through radiation pressure on dust grains. We find a coupling efficiency between the bolometric luminosity of the quasar and the outflow of 0.1$\%$, matching the theoretical prediction of the minimum coupling efficiency necessary for negative quasar feedback. The outflow has sufficient energetics to drive the observed turbulence seen in shocked regions of the quasar host galaxy, likely directly responsible for prolonging the time it takes for gas to cool efficiently., Comment: 14 pages, 7 figures, 1 table, submitted for publication in ApJ

The James Webb Space Telescope (JWST) is revolutionizing our knowledge of $z>5$ galaxies and their actively accreting black holes. Using the JWST Cycle 1 Treasury program Ultradeep NIRSpec and NIRCam ObserVations before the Epoch of Reionization (UNCOVER) in the lensing field Abell 2744, we report the identification of a sample of little red dots at $3 < z_{\rm{phot}} < 7$ that likely contain highly-reddened accreting supermassive black holes. Using a NIRCam-only selection to F444W$<27.7$ mag, we find 26 sources over the $\sim45$ arcmin$^{2}$ field that are blue in F115W$-$F200W$\sim0$ (or $\beta_{\rm UV}\sim-2.0$ for $f_{\lambda} \propto \lambda^\beta$), red in F200W$-$F444W = $1-4$ ($\beta_{\rm opt} \sim +2.0$), and are dominated by a point-source like central component. Of the 20 sources with deep ALMA 1.2-mm coverage, none are detected individually or in a stack. For the majority of the sample, SED fits to the JWST+ALMA observations prefer models with hot dust rather than obscured star-formation to reproduce the red NIRCam colors and ALMA 1.2-mm non-detections. While compact dusty star formation can not be ruled out, the combination of extremely small sizes ($\langle r_e \rangle\approx50$ pc after correction for magnification), red rest-frame optical slopes, and hot dust can by explained by reddened broad-line active galactic nuclei (AGNs). Our targets have faint $M_{\rm 1450} \approx -14\ \, {\rm to} -18$ mag but inferred bolometric luminosities of $L_{\rm bol} = 10^{43}-10^{46}$ erg/s, reflecting their obscured nature. If the candidates are confirmed as AGNs with upcoming UNCOVER spectroscopy, then we have found an abundant population of reddened luminous AGN that are at least ten times more numerous than UV-luminous AGN at the same intrinsic bolometric luminosity., Comment: submitted to ApJ

We present the results of a search for high-redshift ($z>9$) galaxy candidates in the JWST UNCOVER survey, using deep NIRCam and NIRISS imaging in 7 bands over $\sim45$ arcmin$^2$ and ancillary HST observations. The NIRCam observations reach a $5-\sigma$ limiting magnitude of $\sim 29.2$ AB. The identification of high$-z$ candidates relies on a combination of a dropout selection and photometric redshifts. We find 16 candidates at $99$. A few galaxies at $z\sim9-10$ appear to show a clear Balmer break between the F356W and F444W/F410M bands, which helps constrain their stellar mass. We estimate blue UV continuum slopes between $\beta=-1.8$ and $-2.3$, typical for early galaxies at $z>9$ but not as extreme as the bluest recently discovered sources. We also find evidence for a rapid redshift-evolution of the mass-luminosity relation and a redshift-evolution of the UV continuum slope for a given range of intrinsic magnitude, in line with theoretical predictions. These findings suggest that deeper JWST observations are needed to reach the fainter galaxy population at those early epochs, and follow-up spectroscopy will help better constrain the physical properties and star formation histories of a larger sample of galaxies., Comment: Submitted to MNRAS

Quasar feedback may regulate the growth of supermassive black holes, quench coeval star formation, and impact galaxy morphology and the circumgalactic medium. However, direct evidence for quasar feedback in action at the epoch of peak black hole accretion at z ~ 2 remains elusive. A good case in point is the z = 1.6 quasar WISEA J100211.29+013706.7 (XID 2028) where past analyses of the same ground-based data have come to different conclusions. Here we revisit this object with the integral field unit of the Near Infrared Spectrograph (NIRSpec) on board the James Webb Space Telescope (JWST) as part of Early Release Science program Q3D. The excellent angular resolution and sensitivity of the JWST data reveal new morphological and kinematic sub-structures in the outflowing gas plume. An analysis of the emission line ratios indicates that photoionization by the central quasar dominates the ionization state of the gas with no obvious sign for a major contribution from hot young stars anywhere in the host galaxy. Rest-frame near-ultraviolet emission aligned along the wide-angle cone of outflowing gas is interpreted as a scattering cone. The outflow has cleared a channel in the dusty host galaxy through which some of the quasar ionizing radiation is able to escape and heat the surrounding interstellar and circumgalactic media. The warm ionized outflow is not powerful enough to impact the host galaxy via mechanical feedback, but radiative feedback by the AGN, aided by the outflow, may help explain the unusually small molecular gas mass fraction in the galaxy host., Comment: 17 pages, 9 figures, accepted for publication in The Astrophysical Journal

Massive galaxies formed most actively at redshifts $z=1-3$ during the period known as `cosmic noon.' Here we present an emission-line study of an extremely red quasar SDSSJ165202.64+172852.3 host galaxy at $z=2.94$, based on observations with the Near Infrared Spectrograph (NIRSpec) integral field unit (IFU) on board JWST. We use standard emission-line diagnostic ratios to map the sources of gas ionization across the host and a swarm of companion galaxies. The quasar dominates the photoionization, but we also discover shock-excited regions orthogonal to the ionization cone and the quasar-driven outflow. These shocks could be merger-induced or -- more likely, given the presence of a powerful galactic-scale quasar outflow -- these are signatures of wide-angle outflows that can reach parts of the galaxy that are not directly illuminated by the quasar. Finally, the kinematically narrow emission associated with the host galaxy presents as a collection of 1 kpc-scale clumps forming stars at a rate of at least 200 $M_{\odot}$ yr$^{-1}$. The ISM within these clumps shows high electron densities, reaching up to 3,000 cm$^{-3}$ with metallicities ranging from half to a third solar with a positive metallicity gradient and V band extinctions up to 3 magnitudes. The star formation conditions are far more extreme in these regions than in local star-forming galaxies but consistent with that of massive galaxies at cosmic noon. JWST observations reveal an archetypical rapidly forming massive galaxy undergoing a merger, a clumpy starburst, an episode of obscured near-Eddington quasar activity, and an extremely powerful quasar outflow simultaneously., Comment: 19 pages, 8 figures. Accepted for publication in ApJ

We present visual classifications of merger-induced tidal disturbances in 143 $\rm{M}_* \sim 10^{11}\rm{M}_\odot$ post-starburst galaxies at z$\sim$0.7 identified in the SQuIGG$\vec{L}$E Sample. This sample spectroscopically selects galaxies from the Sloan Digital Sky Survey that have stopped their primary epoch of star formation within the past $\sim$500 Myrs. Visual classifications are performed on Hyper Suprime Cam (HSC) i-band imaging. We compare to a control sample of mass- and redshift-matched star-forming and quiescent galaxies from the Large Early Galaxy Census and find that post-starburst galaxies are more likely to be classified as disturbed than either category. This corresponds to a factor of $3.6^{+2.9}_{-1.3}$ times the disturbance rate of older quiescent galaxies and $2.1^{+1.9}_{-.73}$ times the disturbance rate of star-forming galaxies. Assuming tidal features persist for $\lesssim500$ Myr, this suggests merging is coincident with quenching in a significant fraction of these post-starbursts. Galaxies with tidal disturbances are younger on average than undisturbed post-starburst galaxies in our sample, suggesting tidal features from a major merger may have faded over time. This may be exacerbated by the fact that, on average, the undisturbed subset is fainter, rendering low surface brightness tidal features harder to identify. However, the presence of ten young ($\lesssim150$ Myr since quenching) undisturbed galaxies suggests that major mergers are not the only fast physical mechanism that shut down the primary epoch of star formation in massive galaxies at intermediate redshift., Comment: 13 pages, 8 figures, accepted to ApJ

Quantifying the connection between galaxies and their host dark matter halos has been key for testing cosmological models on various scales. Below $M_\star \sim 10^9\,M_\odot$, such studies have primarily relied on the satellite galaxy population orbiting the Milky Way. Here we present new constraints on the connection between satellite galaxies and their host dark matter subhalos using the largest sample of satellite galaxies in the Local Volume ($D \lesssim 12\,\mathrm{Mpc}$) to date. We use $250$ confirmed and $71$ candidate dwarf satellites around 27 Milky Way (MW)-like hosts from the Exploration of Local VolumE Satellites (ELVES) Survey and use the semi-analytical SatGen model for predicting the population of dark matter subhalos expected in the same volume. Through a Bayesian model comparison of the observed and the forward-modeled satellite stellar mass functions (SSMF), we infer the satellite stellar-to-halo mass relation. We find that the observed SSMF is best reproduced when subhalos at the low mass end are populated by a relation of the form $M_\star \propto M^\alpha_\mathrm{peak}$, with a moderate slope of $\alpha_\mathrm{const}=2.10 \pm 0.01$ and a low scatter, constant as a function of the peak halo mass, of $\sigma_\mathrm{const}=0.06^{+0.07}_{-0.05}$. A model with a steeper slope ($\alpha_\mathrm{grow}=2.39 \pm 0.06$) and a scatter that grows with decreasing $M_\mathrm{peak}$ is also consistent with the observed SSMF but is not required. Our new model for the satellite-subhalo connection, based on hundreds of Local Volume satellite galaxies, is in line with what was previously derived using only the Milky Way satellites., Comment: Accepted for publication in ApJ. Figure 8 shows the key result -- the Satellite Stellar to Halo Mass relation obtained in this work, in comparison to previous studies

Extremely red quasars, with bolometric luminosities exceeding $10^{47}$ erg s$^{-1}$, are a fascinating high-redshift population that is absent in the local universe. They are the best candidates for supermassive black holes accreting at rates at or above the Eddington limit, and they are associated with the most rapid and powerful outflows of ionized gas known to date. They are also hosted by massive galaxies. Here we present the first integral field unit (IFU) observations of a high-redshift quasar obtained by the Near Infrared Spectrograph (NIRSpec) on board the James Webb Space Telescope (JWST), which targeted SDSSJ165202.64+172852.3, an extremely red quasar at $z=2.94$. JWST observations reveal extended ionized gas - as traced by [OIII]$\lambda$5007\AA - in the host galaxy of the quasar, its outflow, and the circumgalactic medium. The complex morphology and kinematics imply that the quasar resides in a very dense environment with several interacting companion galaxies within projected distances of 10-15 kpc. The high density of the environment and the large velocities of the companion galaxies suggest that this system may represent the core of a forming cluster of galaxies. The system is a good candidate for a merger of two or more dark matter halos, each with a mass of a few $10^{13}$ M$_\odot$ and traces potentially one of the densest knots at $z\sim3$., Comment: ApJL in press, 11 pages, please also see associated ESA / STScI press release from Oct 20th 2022

The James Webb Space Telescope (JWST) will have the sensitivity to detect early low-mass black holes as they transition from "seeds" to supermassive black holes (BHs). Based on the JAGUAR mock catalog of galaxies, we present a clean color selection that takes advantage of the unique UV slope of accreting supermassive black holes with relatively low mass and high accretion rates. We show for a range of host stellar masses that those galaxies hosting $\sim 10^{6} M_{\odot}$ BHs radiating at >10% of their Eddington luminosity separate in color space from inactive systems. Here we propose a set of 3-band, 2-color selection boxes (with 90% completeness; 90% purity; balanced purity/completeness) with JWST/NIRCam to identify the most promising growing BH candidates at $z \sim 7-10$., Comment: Updated to accepted version. 11 pages, 4 figures, 1 table. Published in ApJL

High-accuracy black hole (BH) masses require excellent spatial resolution that is only achievable for galaxies within ~100 Mpc using present-day technology. At larger distances, BH masses are often estimated with single-epoch scaling relations for active galactic nuclei. This method requires only luminosity and the velocity dispersion of the broad line region (BLR) to calculate a virial product, and an additional virial factor, $f$, to determine BH mass. The accuracy of these single-epoch masses, however, is unknown, and there are few empirical constraints on the variance of $f$ between objects. We attempt to calibrate single-epoch BH masses using spectropolarimetric measurements of nine megamaser galaxies from which we measure the velocity distribution of the BLR. We do not find strong evidence for a correlation between the virial products used for single-epoch masses and dynamical mass, neither for the megamaser sample alone or when combined with dynamical masses from reverberation mapping modeling. Furthermore, we find evidence that the virial parameter $f$ varies between objects, but we do not find strong evidence for a correlation with other observable parameters such as luminosity or broad line width. Although we cannot definitively rule out the existence of any correlation between dynamical mass and virial product, we find tension between allowed $f$ values for masers and those widely used in the literature. We conclude that the single-epoch method requires further investigation if it is to be used successfully to infer BH masses., Comment: 27 pages, 17 figures, accepted by ApJ

We present the statistical redshift distribution of a large sample of low surface brightness (LSB) galaxies identified in the first 200 deg$^2$ of the Hyper Suprime-Cam Strategic Survey Program. Through cross-correlation with the NASA-SDSS Atlas, we find that the majority of objects lie within z<0.15 or ~ 500 Mpc, yielding a mass range of $M_*$ ~ $10^7-10^9$ M_sun and size range of $r_{\rm{eff,g}}$ ~ 1-8 kpc. We find a peak in the distance distribution within 100 Mpc, corresponding mostly to ~ $10^7$ M_sun galaxies that fall on the known mass-size relation. There is also a tail in the redshift distribution out to z~0.15, comprising more massive ($M_*=10^8-10^9$ M_sun) galaxies at the larger end of our size range. We see tentative evidence that at the higher-mass end ($M_* > 10^8$ M_sun) the LSB galaxies do not form a smooth extension of the mass-size relation of higher surface brightness galaxies, perhaps suggesting that the LSB galaxy population is distinct in its formation path., Comment: 15 pages, 6 figures, resubmitted to ApJ after review

We present structural measurements of 145 spectroscopically selected intermediate-redshift (z$\sim$0.7), massive ($M_\star \sim 10^{11} \ M_\odot$) post-starburst galaxies from the SQuIGG$\vec{L}$E Sample measured using wide-depth Hyper Suprime-Cam i-band imaging. This deep imaging allows us to probe the sizes and structures of these galaxies, which we compare to a control sample of star forming and quiescent galaxies drawn from the LEGA-C Survey. We find that post-starburst galaxies systematically lie $\sim0.1$ dex below the quiescent mass-size (half-light radius) relation, with a scatter of $\sim0.2$ dex. This finding is bolstered by non-parametric measures, such as the Gini coefficient and the concentration, which also reveal these galaxies to have more compact light profiles than both quiescent and star-forming populations at similar mass and redshift. The sizes of post-starburst galaxies show either negative or no correlation with the time since quenching, such that more recently quenched galaxies are larger or similarly sized. This empirical finding disfavors the formation of post-starburst galaxies via a purely central burst of star formation that simultaneously shrinks the galaxy and shuts off star formation. We show that the central densities of post-starburst and quiescent galaxies at this epoch are very similar, in contrast with their effective radii. The structural properties of z$\sim$0.7 post-starburst galaxies match those of quiescent galaxies that formed in the early universe, suggesting that rapid quenching in the present epoch is driven by a similar mechanism to the one at high redshift., Comment: 28 pages, 19 figures, 2 appendices, updated to accepted ApJ version

There is growing evidence for physical influence between supermassive black holes and their host galaxies. We present a case study of nearby galaxy NGC 7582, for which we find evidence that galactic substructure plays an important role in affecting the collimation of ionized outflows as well as contributing to the heavy active galactic nucleus (AGN) obscuration. This result contrasts with a simple, small-scale AGN torus model, according to which AGN wind collimation may take place inside the torus itself, at subparsec scale. Using 3D spectroscopy with the MUSE instrument, we probe the kinematics of the stellar and ionized gas components as well as the ionization state of the gas from a combination of emission line ratios. We report for the first time a kinematically distinct core (KDC) in NGC 7582, on a scale of ~600pc. This KDC coincides spatially with dust lanes and starbursting complexes previously observed. We interpret it as a circumnuclear ring of stars and dusty, gas-rich material. We obtain a clear view of the outflowing cones over kpc scales, and demonstrate that they are predominantly photoionized by the central engine. We detect the back cone (behind the galaxy), and confirm previous results of a large nuclear obscuration of both the stellar continuum and HII regions. While we tentatively associate the presence of the KDC to a large-scale bar and/or a minor galaxy merger, we stress the importance of gaining a better understanding of the role of galaxy substructure in controlling the fueling, feedback and obscuration of AGN., Comment: 31 pages, 16 figures (including two animations), accepted for publication in the Astrophysical Journal. For a preview of the animations, see https://noirlab.edu/science/videos/rotating-gas-disk-of-the-ngc-7582-galaxy for Figure 15, and see https://noirlab.edu/science/videos/black-hole-wind-from-the-ngc-7582-galaxy for Figure 16

In the coming decade, a new generation of massively multiplexed spectroscopic surveys, such as PFS, WAVES, and MOONS, will probe galaxies in the distant universe in vastly greater numbers than was previously possible. In this work, we generate mock catalogs for each of these three planned surveys to help quantify and optimize their scientific output. To assign photometry into the UniverseMachine empirical model, we develop the Calibrating Light: Illuminating Mocks By Empirical Relations (CLIMBER) procedure using UltraVISTA photometry. Using the published empirical selection functions for each aforementioned survey, we quantify the mass completeness of each survey. We compare different targeting strategies by varying the area and targeting completeness, and quantify how these survey parameters affect the uncertainty of the two-point correlation function. We demonstrate that the PFS and MOONS measurements will be primarily dominated by cosmic variance, not shot noise, motivating the need for increasingly large survey areas. On the other hand, the WAVES survey, which covers a much larger area, will strike a good balance between cosmic variance and shot noise. For a fixed number of targets, a 5% increased survey area (and $\sim$5% decreased completeness) would decrease the uncertainty of the correlation function at intermediate scales by 0.15%, 1.2%, and 1.1% for our WAVES, PFS, and MOONS samples, respectively. Meanwhile, for a fixed survey area, 5% increased targeting completeness improves the same constraints by 0.7%, 0.25%, and 0.1%. All of the utilities used to construct our mock catalogs and many of the catalogs themselves are publicly available.

We describe the SQuIGGLE survey of intermediate-redshift post-starburst galaxies. We leverage the large sky coverage of the SDSS to select ~1300 recently-quenched galaxies at 0.575% of SQuIGGLE galaxies formed at least a quarter of their total stellar mass in the recent burst, which ended just ~200Myr before observation. We find that SQuIGGLE galaxies are on average younger and more burst-dominated than most other z<~1 post-starburst samples. This large sample of bright post-starburst galaxies at intermediate redshift opens a wide range of studies into the quenching process. In particular, the full SQuIGGLE survey will investigate the molecular gas reservoirs, morphologies, kinematics, resolved stellar populations, AGN incidence, and infrared properties of this unique sample of galaxies in order to place definitive constraints on the quenching process., Comment: 23 pages, 16 figures, accepted to ApJ

The nanohertz gravitational wave background (GWB) is believed to be dominated by GW emission from supermassive black hole binaries (SMBHBs). Observations of several dual active galactic nuclei (AGN) strongly suggest a link between AGN and SMBHBs, given that these dual AGN systems will eventually form bound binary pairs. Here we develop an exploratory SMBHB population model based on empirically constrained quasar populations, allowing us to decompose the GWB amplitude into an underlying distribution of SMBH masses, SMBHB number density, and volume enclosing the GWB. Our approach also allows us to self-consistently predict the number of local SMBHB systems from the GWB amplitude. Interestingly, we find the local number density of SMBHBs implied by the common-process signal in the NANOGrav 12.5-yr dataset to be roughly five times larger than previously predicted by other models. We also find that at most $\sim 25 \%$ of SMBHBs can be associated with quasars. Furthermore, our quasar-based approach predicts $\gtrsim 95\%$ of the GWB signal comes from $z \lesssim 2.5$, and that SMBHBs contributing to the GWB have masses $\gtrsim 10^8 M_\odot$. We also explore how different empirical galaxy-black hole scaling relations affect the local number density of GW sources, and find that relations predicting more massive black holes decrease the local number density of SMBHBs. Overall, our results point to the important role that a measurement of the GWB will play in directly constraining the cosmic population of SMBHBs, as well as their connections to quasars and galaxy mergers., Comment: 12 pages, 6 figures, 1 appendix, accepted for publication in ApJ

In this study, we investigate the X-ray properties of WISE J090924.01+000211.1 (WISEJ0909+0002), an extremely luminous infrared (IR) galaxy (ELIRG) at $z_{\rm spec}$= 1.871 in the eROSITA final equatorial depth survey (eFEDS). WISEJ0909+0002 is a WISE 22 $\mu$m source, located in the GAMA-09 field, detected by eROSITA during the performance and verification phase. The corresponding optical spectrum indicates that this object is a type-1 active galactic nucleus (AGN). Observations from eROSITA combined with Chandra and XMM-Newton archival data indicate a very luminous ($L$ (2--10 keV) = ($2.1 \pm 0.2) \times 10^{45}$ erg s$^{-1}$) unobscured AGN with a power-law photon index of $\Gamma$ = 1.73$_{-0.15}^{+0.16}$, and an absorption hydrogen column density of $\log\,(N_{\rm H}/{\rm cm}^{-2}) < 21.0$. The IR luminosity was estimated to be $L_{\rm IR}$ = (1.79 $\pm$ 0.09) $\times 10^{14}\, L_{\odot}$ from spectral energy distribution modeling based on 22 photometric data (X-ray to far-IR) with X-CIGALE, which confirmed that WISEJ0909+0002 is an ELIRG. A remarkably high $L_{\rm IR}$ despite very low $N_{\rm H}$ would indicate that we are witnessing a short-lived phase in which hydrogen gas along the line of sight is blown outwards, whereas warm and hot dust heated by AGNs still exist. As a consequence of eROSITA all-sky survey, $6.8_{-5.6}^{+16}\times 10^2$ such X-ray bright ELIRGs are expected to be discovered in the entire extragalactic sky ($|b| > 10^\circ$). This can potentially be the key population to constrain the bright-end of IR luminosity functions., Comment: 10 pages, 5 figures, and 3 tables, accepted for publication in A&A Letters (special Issue: First science highlights from SRG/eROSITA)

Classifying catalog objects as stars, galaxies, or AGN is a crucial part of any statistical study of galaxies. We describe our pipeline for binary (star/galaxy) and multiclass (star/galaxy/Type I AGN/Type II AGN) classification developed for the very deep CLAUDS+HSC-SSP $u^*grizy$ dataset. Our method uses the XGBoost implementation of Gradient Boosted Trees (GBT) to train ensembles of models which take photometry, colours, maximum surface brightnesses, and effective radii from all available bands as input, and output the probability that an object belongs to each of the classes under consideration. At $i_{AB}<25$ our binary star/galaxy model has AUC=0.9974 and at the threshold that maximizes our sample's weighted F1 score, selects a sample of galaxies with 99.7% purity and 99.8% completeness. We test the model's ability to generalize to objects fainter than those seen during training and find that extrapolation of ~1-2 magnitudes is reasonable for most applications provided that the galaxies in the training sample are representative of the range of redshifts and colours of the galaxies in the target sample. We also perform an exploratory analysis of the method's ability to identify AGN using a small x-ray selected sample and find that it holds promise for classifying type I AGN, although it performs less well for type II AGN. Our results demonstrate that GBTs provide a flexible, robust and efficient method for performing classification of catalog objects in large astronomical imaging surveys., Comment: Accepted for publication in MNRAS

We report on deep XMM-Newton and NuSTAR observations of the high redshift, z=2.94, extremely red quasar (ERQ), SDSS J165202.60+172852.4, with known galactic ionized outflows detected via spatially-resolved [OIII] emission lines. X-ray observations allow us to directly probe the accretion disk luminosity and the geometry and scale of the circumnuclear obscuration. We fit the spectra from the XMM-Newton/EPIC and NuSTAR detectors with a physically motivated torus model and constrain the source to exhibit a near Compton-thick column density of $N_H=(1.02^{+0.76}_{-0.41})\times10^{24}\textrm{ cm}^{-2}$, a near edge-on geometry with the line-of-sight inclination angle of $\theta_i=85^{\circ}$, and a scattering fraction of $f_{sc}\sim 3$ %. The absorption-corrected, intrinsic 2-10 keV X-ray luminosity of $L_{\textrm{2-10}}=(1.4^{+1}_{-1}) \times10^{45} \textrm{ erg s}^{-1}$ reveals a powerful quasar that is not intrinsically X-ray weak, consistent with observed trends in other ERQs. We also estimate the physical properties of the obscuration, although highly uncertain: the warm ionized scattering density of $n_e \sim 7.5\times(10^2-10^3)\textrm{ cm}^{-3}$ and the obscuration mass of $M_{obsc} \sim 1.7\times(10^4-10^6) M_{\odot}$. As previously suggested with shallower X-ray observations, optical and infrared selection of ERQ has proved effective in finding obscured quasars with powerful outflow signatures. Our observations provide an in-depth view into the X-ray properties of ERQs and support the conclusions of severely photon-limited studies of obscured quasar populations at high redshifts., Comment: Accepted for publication in MNRAS. 12 pages, 6 Figures, 4 Tables

We have performed a search over 3440 deg$^2$ of Epoch 1 (2017-2019) of the Very Large Array Sky Survey to identify unobscured quasars in the optical ($0.2 < z < 3.2$) and obscured active galactic nuclei (AGN) in the infrared that have brightened dramatically in the radio over the past one to two decades. These sources would have been previously classified as "radio-quiet" quasars based on upper limits from the Faint Images of the Radio Sky at Twenty Centimeters survey (1993-2011), but they are now consistent with "radio-loud" quasars ($L_{\rm 3\,GHz} = 10^{40 - 42} \,\, {\rm erg} \,{\rm s}^{-1}$). A quasi-simultaneous, multiband ($\sim1-18$ GHz) follow-up study of 14 sources with the VLA has revealed compact sources ($<0.1^{\prime \prime}$ or $<1$ kpc) with peaked radio spectral shapes. The high-amplitude variability over decadal timescales at 1.5 GHz (100% to $>$2500%), but roughly steady fluxes over a few months at 3 GHz, are inconsistent with extrinsic variability due to propagation effects, thus favoring an intrinsic origin. We conclude that our sources are powerful quasars hosting compact/young jets. This challenges the generally accepted idea that "radio-loudness" is a property of the quasar/AGN population that remains fixed on human timescales. Our study suggests that frequent episodes of short-lived AGN jets that do not necessarily grow to large scales may be common at high redshift. We speculate that intermittent but powerful jets on subgalactic scales could interact with the interstellar medium, possibly driving feedback capable of influencing galaxy evolution., Comment: 28 pages, 9 figures, accepted to ApJ

We increase the sample of ultra diffuse galaxies (UDGs) in lower density environments with characterized globular cluster (GC) populations using new Hubble Space Telescope observations of nine UDGs in group environments. While the bulk of our UDGs have GC abundances consistent with normal dwarf galaxies, two of these UDGs have excess GC populations. These two UDGs both have GC luminosity functions consistent with higher surface brightness galaxies and cluster UDGs. We then combine our nine objects with previous studies to create a catalog of UDGs with analyzed GC populations that spans a uniquely diverse range of environments. We use this catalog to examine broader trends in the GC populations of low stellar mass galaxies. The highest GC abundances are found in cluster UDGs, but whether cluster UDGs are actually more extreme requires study of many more UDGs in groups. We find a possible positive correlation between GC abundance and stellar mass, and between GC abundance and galaxy size at fixed stellar mass. However, we see no significant stellar-mass galaxy-size relation, over our limited stellar mass range. We consider possible origins of the correlation between GC abundance and galaxy size, including the possibility that these two galaxy properties are both dependent on the galaxy dark matter halo, or that they are related through baryonic processes like internal feedback., Comment: 22 pages, 12 figures, submitted to ApJ, posted after referee's comments

Though smooth, extended spheroidal stellar outskirts have long been observed around nearby dwarf galaxies, it is unclear whether dwarfs generically host an extended stellar halo. We use imaging from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) to measure the shapes of dwarf galaxies out to four effective radii for a sample of dwarfs at 0.005 >~ 0.75 (where the ellipsoid is characterized by three principle semi-axes constrained by C<=B<=A). At M_star>10^8.5 M_sun, the galaxies grow from thick disk-like near their centers towards the spheroidal extreme at four effective radii. We also see that although blue dwarfs are, on average, characterized by thinner discs than red dwarfs, both blue and red dwarfs grow more spheroidal as a function of radius. This relation also holds true for a comparison between field and satellite dwarfs. This uniform trend towards relatively spheroidal shapes as a function of radius is consistent with an in-situ formation mechanism for stellar outskirts around low-mass galaxies, in agreement with proposed models where star formation feedback produces round stellar outskirts around dwarfs., Comment: 23 pages, 15 figures. Submitted to ApJ

Active galactic nuclei (AGNs) are key for understanding the coevolution of galaxies and supermassive black holes (SMBHs). AGN activity is thought to affect the properties of their host galaxies, via a process called "AGN feedback", which drives the co-evolution. From a parent sample of 1151 z < 1 type-1 quasars from the Sloan Digital Sky Survey quasar catalog, we detected host galaxies of 862 of them in the high-quality grizy images of the Subaru Hyper Suprime-Cam (HSC) survey. The unprecedented combination of the survey area and depth allows us to perform a statistical analysis of the quasar host galaxies, with small sample variance. We fit the radial image profile of each quasar as a linear combination of the point spread function and the Sersic function, decomposing the images into the quasar nucleus and the host galaxy components. We found that the host galaxies are massive, with stellar mass Mstar > 10^(10) Msun, and are mainly located on the green valley. This trend is consistent with a scenario in which star formation of the host galaxies is suppressed by AGN feedback, that is, AGN activity may be responsible for the transition of these galaxies from the blue cloud to the red sequence. We also investigated the SMBH mass to stellar mass relation of the z < 1 quasars, and found a consistent slope with the local relation, while the SMBHs may be slightly undermassive. However, the above results are subject to our sample selection, which biases against host galaxies with low masses and/or large quasar-to-host flux ratios., Comment: Accepted for publication in PASJ

We study the incidence of nuclear activity in a large sample of massive post-starburst galaxies at z~0.7 selected from the Sloan Digital Sky Survey, and identify active galactic nuclei based on radio continuum and optical emission lines. Over our mass range of 10^10.6-10^11.5 Msun, the incidence of radio activity is weakly dependent on stellar mass and independent of stellar age, while radio luminosity depends strongly on stellar mass. Optical nuclear activity incidence depends most strongly on the Dn4000 line index, a proxy for stellar age, with an active fraction that is ~ten times higher in the youngest versus oldest post-starburst galaxies. Since a similar trend is seen between age and molecular gas fractions, we argue that, like in local galaxies, the age trend reflects a peak in available fueling rather than feedback from the central black hole on the surrounding galaxy., Comment: 10 pages, 4 figures, posted after review by ApJL