Search

The shallow potential wells of star-forming dwarf galaxies make their surrounding circumgalactic and intergalactic medium (CGM/IGM) sensitive laboratories for studying the inflows and outflows thought to regulate galaxy evolution. We present new absorption-line measurements in quasar sightlines probing within projected distances of $<300$ kpc from 91 star-forming field dwarf galaxies with a median stellar mass of $\log{M_\star/\rm{M_\odot}} \approx 8.3$ at $0.077 < z < 0.73$ from the Cosmic Ultraviolet Baryon Survey (CUBS). In this redshift range, the CUBS quasar spectra cover a suite of transitions including H I, low and intermediate metal ions (e.g., C II, Si II, C III, and Si III), and highly ionized O VI. This CUBS-Dwarfs survey enables constraints with samples 9$\times$ larger than past dwarf CGM/IGM studies with similar ionic coverage. We find that low and intermediate ionization metal absorption is rare around dwarf galaxies, consistent with previous surveys of local dwarfs. In contrast, highly ionized O VI is commonly observed in sightlines that pass within the virial radius of a dwarf, and O VI detection rates are non-negligible at projected distances of 1$-$2$\times$ the virial radius. Based on these measurements, we estimate that the O VI-bearing phase of the CGM/IGM accounts for a dominant share of the metal budget of dwarf galaxies. The absorption kinematics suggest that a relatively modest fraction of the O VI-bearing gas is formally unbound. Together, these results imply that low-mass systems at $z\lesssim 1$ effectively retain a substantial fraction of their metals within the nearby CGM and IGM., Comment: Submitted to ApJ. 18 pages, 6 figures, 2 tables

Continuum reverberation mapping (CRM) of active galactic nuclei (AGN) monitors multiwavelength variability signatures to constrain accretion disk structure and supermassive black hole (SMBH) properties. The upcoming Vera Rubin Observatory's Legacy Survey of Space and Time (LSST) will survey tens of millions of AGN over the next decade, with thousands of AGN monitored with almost daily cadence in the deep drilling fields. However, existing CRM methodologies often require long computation time and are not designed to handle such large amount of data. In this paper, we present a fast and flexible inference framework for CRM using simulation-based inference (SBI) with deep learning to estimate SMBH properties from AGN light curves. We use a long-short-term-memory (LSTM) summary network to reduce the high-dimensionality of the light curve data, and then use a neural density estimator to estimate the posterior of SMBH parameters. Using simulated light curves, we find SBI can produce more accurate SMBH parameter estimation with $10^3-10^5$ times speed up in inference efficiency compared to traditional methods. The SBI framework is particularly suitable for wide-field RM surveys as the light curves will have identical observing patterns, which can be incorporated into the SBI simulation. We explore the performance of our SBI model on light curves with irregular-sampled, realistic observing cadence and alternative variability characteristics to demonstrate the flexibility and limitation of the SBI framework., Comment: 14 pages, 8 figures, submitted to ApJ. Comments welcome!

Over three decades of reverberation mapping (RM) studies on local broad-line active galactic nuclei (AGNs) have measured reliable black-hole (BH) masses for $> 100$ AGNs. These RM measurements reveal a significant correlation between the Balmer broad-line region size and the AGN optical luminosity (the $R-L$ relation). Recent RM studies for AGN samples with more diverse BH accretion parameters (e.g., mass and Eddington ratio) reveal a substantial intrinsic dispersion around the average $R-L$ relation, suggesting variations in the overall spectral energy distribution shape as functions of accretion parameters. Here we perform a detailed photoionization investigation of expected broad-line properties as functions of accretion parameters, using the latest models for the AGN continuum implemented in {\tt qsosed}. We compare theoretical predictions with observations of a sample of 67 $z\lesssim0.5$ reverberation-mapped AGNs with both rest-frame optical and UV spectra in the moderate-accretion regime (Eddington ratio $\lambda_{\rm Edd}\equiv L/L_{\rm Edd}<0.5$). The UV/optical line strengths and their dependences on accretion parameters can be reasonably well reproduced by the locally-optimally-emitting cloud (LOC) photoionization models. We provide quantitative recipes that use optical/UV line flux ratios to infer the ionizing continuum, which is not directly observable. In addition, photoionization models with universal values of ionization parameter ($\log U_{\rm H}=-2$) and hydrogen density ($\log n({\rm H})=12$) can qualitatively reproduce the observed global $R-L$ relation for the current AGN sample. However, such models fail to reproduce the observed trend of decreasing BLR size with $L/L_{\rm Edd}$ at fixed optical luminosity, which may imply that the gas density increases with the accretion rate., Comment: 23 pages, 15 figures. Submitted to ApJ. Comments welcome!

We report the discovery of large ionized, [O II] emitting circumgalactic nebulae around the majority of thirty UV luminous quasars at $z=0.4-1.4$ observed with deep, wide-field integral field spectroscopy (IFS) with the Multi-Unit Spectroscopy Explorer (MUSE) by the Cosmic Ultraviolet Baryon Survey (CUBS) and MUSE Quasar Blind Emitters Survey (MUSEQuBES). Among the 30 quasars, seven (23%) exhibit [O II] emitting nebulae with major axis sizes greater than 100 kpc, twenty greater than 50 kpc (67%), and 27 (90%) greater than 20 kpc. Such large, optically emitting nebulae indicate that cool, dense, and metal-enriched circumgalactic gas is common in the halos of luminous quasars at intermediate redshift. Several of the largest nebulae exhibit morphologies that suggest interaction-related origins. We detect no correlation between the sizes and cosmological dimming corrected surface brightnesses of the nebulae and quasar redshift, luminosity, black hole mass, or radio-loudness, but find a tentative correlation between the nebulae and rest-frame [O II] equivalent width in the quasar spectra. This potential trend suggests a relationship between ISM content and gas reservoirs on CGM scales. The [O II]-emitting nebulae around the $z\approx1$ quasars are smaller and less common than Ly$\alpha$ nebulae around $z\approx3$ quasars. These smaller sizes can be explained if the outer regions of the Ly$\alpha$ halos arise from scattering in more neutral gas, by evolution in the cool CGM content of quasar host halos, by lower-than-expected metallicities on $\gtrsim50$ kpc scales around $z\approx1$ quasars, or by changes in quasar episodic lifetimes between $z=3$ and $1$., Comment: 18 pages, 5 figures, 2 tables. Accepted for publication in the Astrophysical Journal

We investigate the group-scale environment of 15 luminous quasars (luminosity $L_{\rm 3000}>10^{46}$ erg s$^{-1}$) from the Cosmic Ultraviolet Baryon Survey (CUBS) at redshift $z\approx1$. Using the Multi Unit Spectroscopic Explorer (MUSE) integral field spectrograph on the Very Large Telescope (VLT), we conduct a deep galaxy redshift survey in the CUBS quasar fields to identify group members and measure the physical properties of individual galaxies and galaxy groups. We find that the CUBS quasars reside in diverse environments. The majority (11 out of 15) of the CUBS quasars reside in overdense environments with typical halo masses exceeding $10^{13}{\rm M}_{\odot}$, while the remaining quasars reside in moderate-size galaxy groups. No correlation is observed between overdensity and redshift, black hole (BH) mass, or luminosity. Radio-loud quasars (5 out of 15 CUBS quasars) are more likely to be in overdense environments than their radio-quiet counterparts in the sample, consistent with the mean trends from previous statistical observations and clustering analyses. Nonetheless, we also observe radio-loud quasars in moderate groups and radio-quiet quasars in overdense environments, indicating a large scatter in the connection between radio properties and environment. We find that the most UV luminous quasars might be outliers in the stellar mass-to-halo mass relations or may represent departures from the standard single-epoch BH relations., Comment: 18 pages, 8 figures, accepted for publication in ApJ

This paper presents a newly established sample of 103 unique galaxies or galaxy groups at $0.4\lesssim z\lesssim 0.7$ from the Cosmic Ultraviolet Baryon Survey (CUBS) for studying the warm-hot circumgalactic medium (CGM) probed by both O VI and Ne VIII absorption. The galaxies and associated neighbors are identified at $< 1$ physical Mpc from the sightlines toward 15 CUBS QSOs at $z_{\rm QSO}\gtrsim 0.8$. A total of 30 galaxies or galaxy groups exhibit associated O VI $\lambda\lambda$ 1031, 1037 doublet absorption within a line-of-sight velocity interval of $\pm250$ km/s, while the rest show no trace of O VI to a detection limit of $\log N_{\rm OVI}/{\rm cm^{-2}}\approx13.7$. Meanwhile, only five galaxies or galaxy groups exhibit the Ne VIII $\lambda\lambda$ 770,780 doublet absorption, down to a limiting column density of $\log N_{\rm NeVIII}/{\rm cm^{-2}}\approx14.0$. These O VI- and Ne VIII-bearing halos reside in different galaxy environments with stellar masses ranging from $\log M_{\rm star}/M_\odot \approx 8$ to $\approx11.5$. The warm-hot CGM around galaxies of different stellar masses and star formation rates exhibits different spatial profiles and kinematics. In particular, star-forming galaxies with $\log M_{\rm star}/M_\odot\approx9-11$ show a significant concentration of metal-enriched warm-hot CGM within the virial radius, while massive quiescent galaxies exhibit flatter radial profiles of both column densities and covering fractions. In addition, the velocity dispersion of O VI absorption is broad with $\sigma_v > 40$ km/s for galaxies of $\log M_{\rm star}/M_\odot>9$ within the virial radius, suggesting a more dynamic warm-hot halo around these galaxies. Finally, the warm-hot CGM probed by O VI and Ne VIII is suggested to be the dominant phase in sub-$L^*$ galaxies with $\log M_{\rm star}/M_\odot\approx9-10$ based on their high ionization fractions in the CGM., Comment: Submitted to ApJ after addressing the referee's comments; 28 pages, 16 figures

High-quality Extragalactic Legacy-field Monitoring (HELM) is a long-term observing program that photometrically monitors several well-studied extragalactic legacy fields with the Dark Energy Camera (DECam) imager on the CTIO 4m Blanco telescope. Since Feb 2019, HELM has been monitoring regions within COSMOS, XMM-LSS, CDF-S, S-CVZ, ELAIS-S1, and SDSS Stripe 82 with few-day cadences in the $(u)gri(z)$ bands, over a collective sky area of $\sim 38$ deg${\rm ^2}$. The main science goal of HELM is to provide high-quality optical light curves for a large sample of active galactic nuclei (AGNs), and to build decades-long time baselines when combining past and future optical light curves in these legacy fields. These optical images and light curves will facilitate the measurements of AGN reverberation mapping lags, as well as studies of AGN variability and its dependences on accretion properties. In addition, the time-resolved and coadded DECam photometry will enable a broad range of science applications from galaxy evolution to time-domain science. We describe the design and implementation of the program and present the first data release that includes source catalogs and the first $\sim 3.5$ years of light curves during 2019A--2022A., Comment: 11 pages, 8 figures, 4 tables, submitted to ApJS. Median source catalogs and light curves of individual objects are publicly available at https://ariel.astro.illinois.edu/helm/

In this paper, we describe the capabilities and constraints of Large Language Models (LLMs) within disparate academic disciplines, aiming to delineate their strengths and limitations with precision. We examine how LLMs augment scientific inquiry, offering concrete examples such as accelerating literature review by summarizing vast numbers of publications, enhancing code development through automated syntax correction, and refining the scientific writing process. Simultaneously, we articulate the challenges LLMs face, including their reliance on extensive and sometimes biased datasets, and the potential ethical dilemmas stemming from their use. Our critical discussion extends to the varying impacts of LLMs across fields, from the natural sciences, where they help model complex biological sequences, to the social sciences, where they can parse large-scale qualitative data. We conclude by offering a nuanced perspective on how LLMs can be both a boon and a boundary to scientific progress.

Turbulent motions in the circumgalactic medium (CGM) play a critical role in regulating the evolution of galaxies, yet their detailed characterization remains elusive. Using two-dimensional velocity maps constructed from spatially-extended [OII] and [OIII] emission, Chen et al. (2023b) measured the velocity structure functions (VSFs) of four quasar nebulae at $z\approx\!0.5$--1.1. One of these exhibits a spectacular Kolmogorov relation. Here we carry out an ensemble study using an expanded sample incorporating four new nebulae from three additional QSO fields. The VSFs measured for all eight nebulae are best explained by subsonic turbulence revealed by the line-emitting gas, which in turn strongly suggests that the cool gas ($T\!\sim\!10^4$ K) is dynamically coupled to the hot ambient medium. Previous work demonstrates that the largest nebulae in our sample reside in group environments with clear signs of tidal interactions, suggesting that environmental effects are vital in seeding and enhancing turbulence within the gaseous halos, ultimately promoting the formation of the extended nebulae. No discernible differences are observed in the VSF properties between radio-loud and radio-quiet QSO fields. We estimate the turbulent heating rate per unit volume, $Q_{\rm turb}$, in the QSO nebulae to be $\sim 10^{-26}$--$10^{-22}$ erg cm$^{-3}$ s$^{-1}$ for the cool phase and $\sim 10^{-28}$--$10^{-25}$ erg cm$^{-3}$ s$^{-1}$ for the hot phase. This range aligns with measurements in the intracluster medium and star-forming molecular clouds but is $\sim10^3$ times higher than the $Q_{\rm turb}$ observed inside cool gas clumps on scales $\lesssim1$ kpc using absorption-line techniques. We discuss the prospect of bridging the gap between emission and absorption studies by pushing the emission-based VSF measurements to below $\approx\!10$ kpc., Comment: 23 pages; 7 figures, and 4 tables in main text; 9 figures in Appendix; accepted by ApJ. Comments welcome

We present the first comprehensive study of a giant, $\approx \! \! 70$ kpc-scale nebula around a radio-quiet quasar at $z<1$. The analysis is based on deep integral field spectroscopy with MUSE of the field of HE$\,$0238$-$1904, a luminous quasar at $z=0.6282$. The nebula emits strongly in $\mathrm{[O \, II]}$, $\rm H \beta$, and $\mathrm{[O \, III]}$, and the quasar resides in an unusually overdense environment for a radio-quiet system. The environment likely consists of two groups which may be merging, and in total have an estimated dynamical mass of $M_{\rm dyn}\approx 4\times 10^{13}$ to $10^{14}\ {\rm M_\odot}$. The nebula exhibits largely quiescent kinematics and irregular morphology. The nebula may arise primarily through interaction-related stripping of circumgalactic and interstellar medium (CGM/ISM) of group members, with some potential contributions from quasar outflows. The simultaneous presence of the giant nebula and a radio-quiet quasar in a rich environment suggests a correlation between such circum-quasar nebulae and environmental effects. This possibility can be tested with larger samples. The upper limits on the electron number density implied by the $\mathrm{[O \, II]}$ doublet ratio range from $\log(n_{\rm e, \, [O \, II]} / \mathrm{cm^{-3}}) < 1.2$ to $2.8$. However, assuming a constant quasar luminosity and negligible projection effects, the densities implied from the measured line ratios between different ions (e.g., $\mathrm{[O\,II]}$, $\mathrm{[O\,III]}$, and $\mathrm{[Ne\,V]}$) and photoionization simulations are often $10{-}400$ times larger. This large discrepancy can be explained by quasar variability on a timescale of $\approx 10^4{-}10^5$ years., Comment: 19 pages, 9 figures, 3 tables; Accepted by MNRAS

This paper presents a new sample of 19 unique galaxies and galaxy groups at $z\approx1$ from the CUBS program, which is designated as the CUBSz1 sample. In this CUBSz1 sample, nine galaxies or galaxy groups show absorption features, while ten systems do not have detectable absorption with 2-$\sigma$ upper limits of log$N$(HeI)/cm$^{-2}\lesssim 13.5$ and log$N$(OV)/cm$^{-2}\lesssim 13.3$. Environmental properties of the galaxies, including galaxy overdensities, the total stellar mass and gravitational potential summed over all nearby neighbors, and the presence of local ionizing sources, are found to have a significant impact on the observed CGM absorption properties. Specifically, massive galaxies and galaxies in overdense regions exhibit a higher rate of incidence of absorption. At the same time, the observed CGM absorption properties in galaxy groups appear to be driven by the galaxy closest to the QSO sightline, rather than by the most massive galaxy or by mass-weighted properties. We introduce a total projected gravitational potential $\psi$, defined as $-\psi/G =\sum M_{{\rm halo}}/d_{{\rm proj}}$ summed over all group members, to characterize the overall galaxy environment. This projected gravitational potential correlates linearly with the maximum density detected in each sightline, consistent with higher-pressure gas being confined in deeper gravitational potential wells. In addition, we find that the radial profile of cool gas density exhibits a general decline from the inner regions to the outskirts, being in pressure balance with the hot halo. Finally, we note that the ionizing flux from nearby galaxies can generate an elevated $N$(HI)/$N$(HeI) ratio, which in turn provides a unique diagnostic of possible local sources contributing to the ionizing radiation field., Comment: 17 pages, 11 figures and 1 table in the main text; 40 figures and 1 table in the appendix; accepted for publication in MNRAS

We present the final data from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project, a precursor to the SDSS-V Black Hole Mapper Reverberation Mapping program. This data set includes 11-year photometric and 7-year spectroscopic light curves for 849 broad-line quasars over a redshift range of 0.1=0.62+-0.07 for the line dispersion measured from the RMS spectrum. The intrinsic scatter of individual virial factors is 0.31+-0.07 dex, indicating a factor of two systematic uncertainty in RM black hole masses. Our lag measurements reveal significant R-L relations for Hbeta and MgII at high redshift, consistent with the latest measurements based on heterogeneous samples. While we are unable to robustly constrain the slope of the R-L relation for CIV given the limited dynamical range in luminosity, we found substantially larger scatter in CIV lags at fixed L1350. Using the SDSS-RM lag sample, we derive improved single-epoch (SE) mass recipes for Hbeta, MgII and CIV, which are consistent with their respective RM masses as well as between the SE recipes from two different lines, over the luminosity range probed by our sample. The new Hbeta and MgII recipes are approximately unbiased estimators at given RM masses, but there are systematic biases in the CIV recipe. The intrinsic scatter of SE masses around RM masses is ~0.45 dex for Hbeta and MgII, increasing to ~0.58 dex for CIV., Comment: Replaced with accepted version (ApJS in press). All measurements remain unchanged from the previous version. 38 pages. Data products available at https://ariel.astro.illinois.edu/sdssrm/final result/ and ftp://quasar.astro.illinois.edu/public/sdssrm/final_result/

We present a high-cadence multi-epoch analysis of dramatic variability of three broad emission lines (MgII, H$\beta$, and H$\alpha$) in the spectra of the luminous quasar ($\lambda L_{\lambda}$(5100\r{A}) = $4.7 \times 10^{44}$ erg s$^{-1}$) SDSS J141041.25+531849.0 at $z = 0.359$ with 127 spectroscopic epochs over 9 years of monitoring (2013-2022). We observe anti-correlations between the broad emission-line widths and flux in all three emission lines, indicating that all three broad emission lines "breathe" in response to stochastic continuum variations. We also observe dramatic radial velocity shifts in all three broad emission lines, ranging from $\Delta{v}$ $\sim$400 km s$^{-1}$ to $\sim$800 km s$^{-1}$, that vary over the course of the monitoring period. Our preferred explanation for the broad-line variability is complex kinematics in the broad-line region gas. We suggest a model for the broad-line variability that includes a combination of gas inflow with a radial gradient, an azimuthal asymmetry (e.g., a hot spot), superimposed on the stochastic flux-driven changes to the optimal emission region ("line breathing"). Similar instances of line-profile variability due to complex gas kinematics around quasars are likely to represent an important source of false positives in radial velocity searches for binary black holes, which typically lack the kind of high-cadence data we analyze here. The long-duration, wide-field, and many-epoch spectroscopic monitoring of SDSS-V BHM-RM provides an excellent opportunity for identifying and characterizing broad emission-line variability, and the inferred nature of the inner gas environment, of luminous quasars.

We measure the correlation between black-hole mass $M_{\rm BH}$ and host stellar mass $M_*$ for a sample of 38 broad-line quasars at $0.2\lesssim z\lesssim 0.8$ (median redshift $z_{\rm med}=0.5$). The black-hole masses are derived from a dedicated reverberation mapping program for distant quasars, and the stellar masses are estimated from two-band optical+IR HST imaging. Most of these quasars are well centered within $\lesssim 1$kpc from the host galaxy centroid, with only a few cases in merging/disturbed systems showing larger spatial offsets. Our sample spans two orders of magnitude in stellar mass ($\sim 10^9-10^{11}\,M_\odot$) and black-hole mass ($\sim 10^7-10^9\,M_\odot$), and reveals a significant correlation between the two quantities. We find a best-fit intrinsic (i.e., selection effects corrected) $M_{\rm BH}-M_{\rm *,host}$ relation of $\log (M_{\rm BH}/M_{\rm \odot})=7.01_{-0.33}^{+0.23} + 1.74_{-0.64}^{+0.64}\log (M_{\rm *,host}/10^{10}M_{\rm \odot})$, with an intrinsic scatter of $0.47_{-0.17}^{+0.24}$dex. Decomposing our quasar hosts into bulges and disks, there is a similar $M_{\rm BH}-M_{\rm *,bulge}$ relation with a slightly larger scatter, likely caused by systematic uncertainties in the bulge-disk decomposition. The $M_{\rm BH}-M_{\rm *,host}$ relation at $z_{\rm med}=0.5$ is similar to that in local quiescent galaxies, with negligible evolution over the redshift range probed by our sample. With direct black-hole masses from reverberation mapping and a large dynamical range of the sample, selection biases do not appear to affect our conclusions significantly. Our results, along with other samples in the literature, suggest that the locally-measured black-hole mass$-$host stellar mass relation is already in place at $z\sim 1$., Comment: 23 pages, 10 figures (Fig 9 is the key figure). Submitted to ApJ. The full figure set and ancillary data products can be found at ftp://quasar.astro.illinois.edu/public/sdssrm/paper_data/Li_2023_HST_host

We present the first empirical constraints on the turbulent velocity field of the diffuse circumgalactic medium around four luminous QSOs at $z\!\approx\!0.5$--1.1. Spatially extended nebulae of $\approx\!50$--100 physical kpc in diameter centered on the QSOs are revealed in [OII]$\lambda\lambda\,3727,3729$ and/or [OIII]$\lambda\,5008$ emission lines in integral field spectroscopic observations obtained using MUSE on the VLT. We measure the second- and third-order velocity structure functions (VSFs) over a range of scales, from $\lesssim\!5$ kpc to $\approx\!20$--50 kpc, to quantify the turbulent energy transfer between different scales in these nebulae. While no constraints on the energy injection and dissipation scales can be obtained from the current data, we show that robust constraints on the power-law slope of the VSFs can be determined after accounting for the effects of atmospheric seeing, spatial smoothing, and large-scale bulk flows. Out of the four QSO nebulae studied, one exhibits VSFs in spectacular agreement with the Kolmogorov law, expected for isotropic, homogeneous, and incompressible turbulent flows. The other three fields exhibit a shallower decline in the VSFs from large to small scales. However, with a limited dynamic range in the spatial scales in seeing-limited data, no constraints can be obtained for the VSF slopes of these three nebulae. For the QSO nebula consistent with the Kolmogorov law, we determine a turbulence energy cascade rate of $\approx\!0.2$ cm$^{2}$ s$^{-3}$. We discuss the implication of the observed VSFs in the context of QSO feeding and feedback in the circumgalactic medium., Comment: 20 pages, 6 figures in main text + 9 figures in appendix, 4 tables; accepted by MNRAS

We report the discovery of giant (50-100 kpc) [O II] emitting nebulae with the Multi-Unit Spectroscopic Explorer (MUSE) in the field of TXS 0206-048, a luminous quasar at z=1.13. Down-the-barrel UV spectra of the quasar show absorption at velocities coincident with those of the extended nebulae, enabling new insights into inflows and outflows around the quasar host. One nebula exhibits a filamentary morphology extending over 120 kpc from the halo toward the quasar and intersecting with another nebula surrounding the quasar host with a radius of 50 kpc. This is the longest cool filament observed to-date and arises at higher redshift and in a less massive system than those in cool-core clusters. The filamentary nebula has line-of-sight velocities >300 km/s from nearby galaxies but matches that of the nebula surrounding the quasar host where they intersect, consistent with accretion of cool inter- or circum-galactic medium or cooling hot halo gas. The kinematics of the nebulae surrounding the quasar host are unusual and complex, with redshifted and blueshifted spiral-like structures. The emission velocities at 5-10 kpc from the quasar match those of inflowing absorbing gas observed in UV spectra of the quasar. Together, the extended nebulae and associated redshifted absorption represent a compelling case of cool, filamentary gas accretion from halo scales into the extended interstellar medium and toward the nucleus of a massive quasar host. The inflow rate implied by the combined emission and absorption constraints is well below levels required to sustain the quasar's radiative luminosity, suggesting anisotropic or variable accretion., Comment: 10 pages (main text), 3 figures, 2 tables, and one appendix figure.. Accepted to ApJ Letters

Dual supermassive black holes (SMBHs) at $\sim$kpc scales are the progenitor population of SMBH mergers and play an important role in understanding the pairing and dynamical evolution of massive black holes in galaxy mergers. Because of the stringent resolution requirement and the apparent rareness of these small-separation pairs, there are scarce observational constraints on this population, with few confirmed dual SMBHs at $<10$kpc separations at $z>1$. Here we present results from a pilot search for kpc-scale dual quasars selected with Gaia Data release 2 (DR2) astrometry and followed up with Hubble Space Telescope (HST) Wide Field Camera 3 dual-band (F475W and F814W) snapshot imaging. Our targets are quasars primarily selected with the varstrometry technique, i.e., light centroid jitter caused by asynchronous variability from both members in an unresolved quasar pair, supplemented by sub-arcsec pairs already resolved by Gaia DR2. We find an overall high fraction of HST-resolved pairs among the varstrometry-selected quasars (unresolved in Gaia DR2), $\sim 30-50\%$, increasing toward high redshift ($\sim 60-80\%$ at $z>1.5$). We discuss the nature of the 43 resolved sub-arcsec pairs based on HST and supplementary data. A substantial fraction ($\sim 40\%$) of these pairs are likely physical quasar pairs or gravitationally lensed quasars. We also discover a triple quasar candidate and a quadruply lensed quasar, which is among the smallest-separation quadruple lenses. These results provide important guidelines to improve varstrometry selection and follow-up confirmation of $\sim$kpc-scale dual SMBHs at high redshift., Comment: 22 pages, 10 figures, 3 tables, accepted to ApJ

Galaxy mergers occur frequently in the early universe and bring multiple supermassive black holes (SMBHs) into the nucleus, where they may eventually coalesce. Identifying post-merger-scale (i.e., <~a few kpc) dual SMBHs is a critical pathway to understanding their dynamical evolution and successive mergers. While serendipitously discovering kpc-scale dual SMBHs at z<1 is possible, such systems are elusive at z>2, but critical to constraining the progenitors of SMBH mergers. The redshift z~2 also marks the epoch of peak activity of luminous quasars, hence probing this spatial regime at high redshift is of particular significance in understanding the evolution of quasars. However, given stringent resolution requirements, there is currently no confirmed <10 kpc physical SMBH pair at z>2. Here we report two sub-arcsec double quasars at z>2 discovered from a targeted search with a novel astrometric technique, demonstrating a high success rate (~50%) in this systematic approach. These high-redshift double quasars could be the long-sought kpc-scale dual SMBHs, or sub-arcsec gravitationally-lensed quasar images. One of these double quasars (at z=2.95) was spatially resolved with optical spectroscopy, and slightly favors the scenario of a physical quasar pair with a projected separation of 3.5 kpc (0.46"). Follow-up observations of double quasars discovered by this targeted approach will be able to provide the first observational constraints on kpc-scale dual SMBHs at z>2., Comment: Published in Nature Astronomy (April 2021)

We present accretion-disk structure measurements from UV-optical reverberation mapping observations of a sample of eight quasars at 0.24

It is well known that reverberation mapping of active galactic nuclei (AGN) reveals a relationship between AGN luminosity and the size of the broad-line region, and that use of this relationship, combined with the Doppler width of the broad emission line, enables an estimate of the mass of the black hole at the center of the active nucleus based on a single spectrum. An unresolved key issue is the choice of parameter used to characterize the line width, either FWHM or line dispersion (the square root of the second moment of the line profile). We argue here that use of FWHM introduces a bias, stretching the mass scale such that high masses are overestimated and low masses are underestimated. Here we investigate estimation of black hole masses in AGNs based on individual or "single epoch" observations, with a particular emphasis in comparing mass estimates based on line dispersion and FWHM. We confirm the recent findings that, in addition to luminosity and line width, a third parameter is required to obtain accurate masses and that parameter seems to be Eddington ratio. We present simplified empirical formulae for estimating black hole masses from the Hbeta 4861 A and C IV 1549 A emission lines. While the AGN continuum luminosity at 5100 A is usually used to predict the Hbeta reverberation lag, we show that the luminosity of the Hbeta broad component can be used instead without any loss of precision, thus eliminating the difficulty of accurately accounting for the host-galaxy contribution to the observed luminosity., Comment: 41 pages, 13 figures, accepted for publication in ApJ

Quasar broad emission lines are largely powered by photoionization from the accretion continuum. Increased central luminosity will enhance line emissivity in more distant clouds, leading to increased average distance of the broad-line-emitting clouds and decreased averaged line width, known as the broad-line region (BLR) "breathing". However, different lines breathe differently, and some high-ionization lines, such as C IV, can even show "anti-breathing" where the line broadens when luminosity increases. Using multi-year photometric and spectroscopic monitoring data from the Sloan Digital Sky Survey Reverberation Mapping project, we quantify the breathing effect ($\Delta$log W=$\alpha\Delta$log L) of broad H$\alpha$, H$\beta$, Mg II, C IV, and C III] for statistical quasar samples over $z\approx 0.1-2.5$. We found that H$\beta$ displays the most consistent normal breathing expected from the virial relation ($\alpha\sim-0.25$), Mg II and H$\alpha$ on average show no breathing ($\alpha\sim 0$), and C IV (and similarly C III] and Si IV mostly shows anti-breathing ($\alpha>0$). The anti-breathing of C IV can be well understood by the presence of a non-varying core component in addition to a reverberating broad-base component, consistent with earlier findings. The deviation from canonical breathing introduces extra scatter (a luminosity-dependent bias) in single-epoch virial BH mass estimates due to intrinsic quasar variability, which underlies the long argued caveats of C IV single-epoch masses. Using the line dispersion instead of FWHM leads to less, albeit still substantial, deviations from canonical breathing in most cases. Our results strengthen the need for reverberation mapping to provide reliable quasar BH masses, and quantify the level of variability-induced bias in single-epoch BH masses based on various lines., Comment: Accepted for publication in ApJ. 34 pages, 26 figures

We present the results of a pilot Hubble Space Telescope (HST) imaging study of the host galaxies of ten quasars from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. Probing more than an order of magnitude in BH and stellar masses, our sample is the first statistical sample to study the BH-host correlations beyond z>0.3 with reliable BH masses from reverberation mapping rather than from single-epoch spectroscopy. We perform image decomposition in two HST bands (UVIS-F606W and IR-F110W) to measure host colors and estimate stellar masses using empirical relations between broad-band colors and the mass-to-light ratio. The stellar masses of our targets are mostly dominated by a bulge component. The BH masses and stellar masses of our sample broadly follow the same correlations found for local RM AGN and quiescent bulge-dominant galaxies, with no strong evidence of evolution in the M_BH-M_*,bulge relation to z~0.6. We further compare the host light fraction from HST imaging decomposition to that estimated from spectral decomposition. We found a good correlation between the host fractions derived with both methods. However, the host fraction derived from spectral decomposition is systematically smaller than that from imaging decomposition by ~30%, indicating different systematics in both approaches. This study paves the way for upcoming more ambitious host galaxy studies of quasars with direct RM-based BH masses at high redshift., Comment: 19 pages, 8 figures, submitted to ApJ

We present reverberation mapping results for the MgII 2800 A broad emission line in a sample of 193 quasars at 0.35

Results from a few decades of reverberation mapping (RM) studies have revealed a correlation between the radius of the broad-line emitting region (BLR) and the continuum luminosity of active galactic nuclei. This "radius-luminosity" relation enables survey-scale black-hole mass estimates across cosmic time, using relatively inexpensive single-epoch spectroscopy, rather than intensive RM time monitoring. However, recent results from newer reverberation mapping campaigns challenge this widely used paradigm, reporting quasar BLR sizes that differ significantly from the previously established radius-luminosity relation. Using simulations of the radius--luminosity relation with the observational parameters of SDSS-RM, we find that this difference is not likely due to observational biases. Instead, it appears that previous RM samples were biased to a subset of quasar properties, and the broader parameter space occupied by the SDSS-RM quasar sample has a genuinely wider range of BLR sizes. We examine the correlation between the deviations from the radius-luminosity relation and several quasar parameters; the most significant correlations indicate that the deviations depend on UV/optical SED and the relative amount of ionizing radiation. Our results indicate that single-epoch black-hole mass estimates that do not account for the diversity of quasars in the radius-luminosity relation could be overestimated by an average of ~0.3 dex., Comment: 11 pages, 12 figures. Submitted to ApJ. Current version includes revisions after referee report

We investigate the performance of different methodologies that measure the time lag between broad-line and continuum variations in reverberation mapping data using simulated light curves that probe a range of cadence, time baseline, and signal-to-noise ratio in the flux measurements. We compare three widely-adopted lag measuring methods: the Interpolated Cross-Correlation Function (ICCF), the z-transformed Discrete Correlation Function (ZDCF) and the MCMC code JAVELIN, for mock data with qualities typical of multi-object spectroscopic reverberation mapping (MOS-RM) surveys that simultaneously monitor hundreds of quasars. We quantify the overall lag detection efficiency, the rate of false detections, and the quality of lag measurements for each of these methods and under different survey designs (e.g., observing cadence and depth) using mock quasar light curves. Overall JAVELIN and ICCF outperform ZDCF in essentially all tests performed. Compared with ICCF, JAVELIN produces higher quality lag measurements, is capable of measuring more lags with timescales shorter than the observing cadence, is less susceptible to seasonal gaps and S/N degradation in the light curves, and produces more accurate lag uncertainties. We measure the Hbeta broad-line region size-luminosity (R-L) relation with each method using the simulated light curves to assess the impact of selection effects of the design of MOS-RM surveys. The slope of the R-L relation measured by JAVELIN is the least biased among the three methods, and is consistent across different survey designs. These results demonstrate a clear preference for JAVELIN over the other two non-parametric methods for MOS-RM programs, particularly in the regime of limited light curve quality as expected from most MOS-RM programs., Comment: 22 pages, 34 figures, Accepted for publication in the Astrophysical Journal

We present reverberation-mapping lags and black-hole mass measurements using the CIV 1549 broad emission line from a sample of 349 quasars monitored as a part of the Sloan Digital Sky Survey Reverberation Mapping Project. Our data span four years of spectroscopic and photometric monitoring for a total baseline of 1300 days. We report significant time delays between the continuum and the CIV 1549 emission line in 52 quasars, with an estimated false-positive detection rate of 10%. Our analysis of marginal lag measurements indicates that there are on the order of 100 additional lags that should be recoverable by adding more years of data from the program. We use our measurements to calculate black-hole masses and fit an updated CIV radius-luminosity relationship. Our results significantly increase the sample of quasars with CIV RM results, with the quasars spanning two orders of magnitude in luminosity toward the high-luminosity end of the CIV radius-luminosity relation. In addition, these quasars are located at among the highest redshifts (z~1.4-2.8) of quasars with black hole masses measured with reverberation mapping. This work constitutes the first large sample of CIV reverberation-mapping measurements in more than a dozen quasars, demonstrating the utility of multi-object reverberation mapping campaigns., Comment: 42 pages (incl. 21 pages of tables): 12 figures, 5 tables. Submitted to AAS journals. Figures 8 and 12 are sets that will be provided in the online published version of the article, and the tables will also be provided fully online

The width of the broad emission lines in quasars is commonly characterized either by the full-width-at-half-maximum (FWHM) or the square root of the second moment of the line profile ($\sigma_{\rm line}$), and used as an indicator of the virial velocity of the broad-line region (BLR) in the estimation of black hole (BH) mass. We measure FWHM and $\sigma_{\rm line}$ for H$\alpha$, H$\beta$ and Mg II broad lines in both the mean and root-mean-square (rms) spectra of a large sample of quasars from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. We introduce a new quantitative recipe to measure $\sigma_{\rm line}$ that is reproducible, less susceptible to noise and blending in the wings, and scales with the intrinsic width of the line. We compare the four definitions of line width (FWHM and $\sigma_{\rm line}$ in mean and rms spectra, respectively) for each of the three broad lines and among different lines. There are strong correlations among different width definitions for each line, providing justification for using the line width measured in single-epoch spectroscopy as a virial velocity indicator. There are also strong correlations among different lines, suggesting alternative lines to H$\beta$ can be used to estimate virial BH masses. We further investigate the correlations between virial BH masses using different line width definitions and the stellar velocity dispersion of the host galaxies, and the dependence of line shape (characterized by the ratio FWHM/$\sigma_{\rm line}$) on physical properties of the quasar. Our results provide further evidence that FWHM is more sensitive to the orientation of a flattened BLR geometry than $\sigma_{\rm line}$, but the overall comparison between the virial BH mass and host stellar velocity dispersion does not provide conclusive evidence that one particular width definition is significantly better than the others., Comment: 14 pages, 9 figures. Accepted by ApJ

The ubiquitous variability of quasars across a wide range of wavelengths and timescales encodes critical information about the structure and dynamics of the circumnuclear emitting regions that are too small to be directly resolved, as well as detailed underlying physics of accretion and feedback processes in these active supermassive black holes. We emphasize the importance of studying quasar variability with time-domain spectroscopy, focusing on two science cases: (1) reverberation mapping (RM) to measure the broad-line region sizes and black hole masses in distant quasars; (2) spectroscopic follow-up of extreme variability quasars that dramatically change their continuum and broad-line flux within several years. We highlight the need for dedicated optical-infrared spectroscopic survey facilities in the coming decades to accompany wide-area time-domain imaging surveys, including: (1) the next phase of the Sloan Digital Sky Survey (SDSS-V; ~2020-2025), an all-sky, time-domain multi-object spectroscopic survey with 2.5m-class telescopes; (2) the planned Maunakea Spectroscopic Explorer, a dedicated 10m-class spectroscopic survey telescope with a 1.5 sq. deg field-of-view and multiplex of thousands of fibers in both optical and near-IR (J+H) to begin operations in 2029; (3) the Time-domain Spectroscopic Observatory (TSO), a proposed Probe-class ~1.3m telescope at L2, with imaging and spectroscopy (R=200, 1800) in 4 bands (0.3 - 5 micron) and rapid slew capability to 90% of sky, which will extend the coverage of Hbeta to z=8., Comment: Science White Paper for the US Astro 2020 Decadal Survey

We present a detailed characterization of the 849 broad-line quasars from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. Our quasar sample covers a redshift range of 0.1

We present accretion-disk structure measurements from continuum lags in the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. Lags are measured using the \texttt{JAVELIN} software from the first-year SDSS-RM $g$ and $i$ photometry, resulting in well-defined lags for 95 quasars, 33 of which have lag SNR $>$ 2$\sigma$. We also estimate lags using the \texttt{ICCF} software and find consistent results, though with larger uncertainties. Accretion-disk structure is fit using a Markov Chain Monte Carlo approach, parameterizing the measured continuum lags as a function of disk size normalization, wavelength, black hole mass, and luminosity. In contrast with previous observations, our best-fit disk sizes and color profiles are consistent (within 1.5~$\sigma$) with the \citet{SS73} analytic solution. We also find that more massive quasars have larger accretion disks, similarly consistent with the analytic accretion-disk model. The data are inconclusive on a correlation between disk size and continuum luminosity, with results that are consistent with both no correlation and with the \citet{SS73} expectation. The continuum lag fits have a large excess dispersion, indicating that our measured lag errors are underestimated and/or our best-fit model may be missing the effects of orientation, spin, and/or radiative efficiency. We demonstrate that fitting disk parameters using only the highest-SNR lag measurements biases best-fit disk sizes to be larger than the disk sizes recovered using a Bayesian approach on the full sample of well-defined lags., Comment: Accepted for publication in ApJ. 20 pages and 15 figures. Complete version of table 1 and full figure set (similar to figure 4) are provided through http://www.phys.uconn.edu/~jtrump/sdssrmdisk_table1.txt http://www.phys.uconn.edu/~jtrump/sdssrmdisk_figset.zip

We present composite broad-line region (BLR) reverberation-mapping lag measurements for \halpha, \hbeta, \HeII\,$\lambda4686$ and \MgII\ for a sample of 144, $z\lesssim 1$ quasars from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. Using only the 32-epoch spectroscopic light curves in the first 6-month season of SDSS-RM observations, we compile correlation-function measurements for individual objects and then coadd them to allow the measurement of the average lags for our sample at mean redshifts of $0.4$ (for \halpha) and $\sim 0.65$ (for the other lines). At similar quasar luminosities and redshifts, the sample-averaged lag decreases in the order of \MgII, \halpha, \hbeta\ and \HeII. This decrease in lags is accompanied by an increase in the mean line width of the four lines, and is roughly consistent with the virialized motion for BLR gas in photoionization equilibrium. These are among the first RM measurements of stratified BLR structure at $z>0.3$. Dividing our sample by luminosity, \halpha\ shows clear evidence of increasing lags with luminosity, consistent with the expectation from the measured BLR size-luminosity relation based on \hbeta. The other three lines do not show a clear luminosity trend in their average lags due to the limited dynamic range of luminosity probed and the poor average correlation signals in the divided samples, a situation that will be improved with the incorporation of additional photometric and spectroscopic data from SDSS-RM. We discuss the utility and caveats of composite-lag measurements for large statistical quasar samples with reverberation-mapping data., Comment: 12 pages, 10 figures, published in ApJ

We present reverberation mapping results from the first year of combined spectroscopic and photometric observations of the Sloan Digital Sky Survey Reverberation Mapping Project. We successfully recover reverberation time delays between the $g+i$-band emission and the broad H$\beta$ emission line for a total of 44 quasars, and for the broad H$\alpha$ emission line in 18 quasars. Time delays are computed using the JAVELIN and CREAM software and the traditional interpolated cross-correlation function (ICCF): Using well defined criteria, we report measurements of 32 H$\beta$ and 13 H$\alpha$ lags with JAVELIN, 42 H$\beta$ and 17 H$\alpha$ lags with CREAM, and 16 H$\beta$ and 8 H$\alpha$ lags with the ICCF. Lag values are generally consistent among the three methods, though we typically measure smaller uncertainties with JAVELIN and CREAM than with the ICCF, given the more physically motivated light curve interpolation and more robust statistical modeling of the former two methods. The median redshift of our H$\beta$-detected sample of quasars is 0.53, significantly higher than that of the previous reverberation-mapping sample. We find that in most objects, the time delay of the H$\alpha$ emission is consistent with or slightly longer than that of H$\beta$. We measure black hole masses using our measured time delays and line widths for these quasars. These black hole mass measurements are mostly consistent with expectations based on the local M-sigma relationship, and are also consistent with single-epoch black hole mass measurements. This work increases the current sample size of reverberation-mapped active galaxies by about two-thirds and represents the first large sample of reverberation mapping observations beyond the local universe (z < 0.3)., Comment: 24 pages, 14 figures, published in ApJ, Volume 851, 1. Figures 6 and 7 are sets that are provided in the online published version of the article, and Table 2 is also fully available online. This version contains updates in Tables 4 and 5, with updated Figures 12, 13, and 14, corresponding to an erratum issued. None of the text required changes; only table values

The recent rapid progress in observations of circumstellar disks and extrasolar planets has reinforced the importance of understanding an intimate coupling between star and planet formation. Under such a circumstance, it may be invaluable to attempt to specify when and how planet formation begins in star-forming regions and to identify what physical processes/quantities are the most significant to make a link between star and planet formation. To this end, we have recently developed a couple of projects. These include an observational project about dust growth in Class 0 YSOs and a theoretical modeling project of the HL Tauri disk. For the first project, we utilize the archive data of radio interferometric observations, and examine whether dust growth, a first step of planet formation, occurs in Class 0 YSOs. We find that while our observational results can be reproduced by the presence of large ($\sim$ mm) dust grains for some of YSOs under the single-component modified blackbody formalism, an interpretation of no dust growth would be possible when a more detailed model is used. For the second project, we consider an origin of the disk configuration around HL Tauri, focusing on magnetic fields. We find that magnetically induced disk winds may play an important role in the HL Tauri disk. The combination of these attempts may enable us to move towards a comprehensive understanding of how star and planet formation are intimately coupled with each other., Comment: Invited talk; 6 pages, 2 figures, 1 table; to appear in proceedings of SFDE conference, eds. D. Johnstone, T. Hoang, F. Nakamura, Q. N. Luong, and J. T. Tranh Van

We are motivated by the recent measurements of dust opacity indices beta around young stellar objects (YSOs), which suggest that efficient grain growth may have occurred earlier than the Class I stage. The present work makes use of abundant archival interferometric observations at submillimeter,millimeter, and centimeter wavelength bands to examine grain growth signatures in the dense inner regions (<1000 AU) of nine Class 0 YSOs. A systematic data analysis is performed to derive dust temperatures, optical depths, and dust opacity indices based on single-component modified black body fittings to the spectral energy distributions (SEDs). The fitted dust opacity indices (beta) are in a wide range of 0.3 to 2.0 when single-component SED fitting is adopted. Four out of the nine observed sources show beta lower than 1.7, the typical value of the interstellar dust. Low dust opacity index (or spectral index) values may be explained by the effect of dust grain growth, which makes beta<1.7. Alternatively, the very small observed values of beta may be interpreted by the presence of deeply embedded hot inner disks, which only significantly contribute to the observed fluxes at long wavelength bands. This possibility can be tested by the higher angular resolution imaging observations of ALMA, or more detailed sampling of SEDs in the millimeter and centimeter bands. The beta values of the remaining five sources are close to or consistent with 1.7, indicating that grain growth would start to significantly reduce the values of beta no earlier than the late-Class 0 stage for these YSOs., Comment: 13 pages, 5 figures, 5 tables, accepted for publication in ApJ

Dust trapping accelerates the coagulation of dust particles, and, thus, it represents an initial step toward the formation of planetesimals. We report H-band (1.6 microns) linear polarimetric observations and 0.87 mm interferometric continuum observations toward a transitional disk around LkH(alpha) 330. As a result, a pair of spiral arms were detected in the H-band emission, and an asymmetric (potentially arm-like) structure was detected in the 0.87 mm continuum emission. We discuss the origin of the spiral arm and the asymmetric structure and suggest that a massive unseen planet is the most plausible explanation. The possibility of dust trapping and grain growth causing the asymmetric structure was also investigated through the opacity index (beta) by plotting the observed spectral energy distribution slope between 0.87 mm from our Submillimeter Array observation and1.3 mm from literature. The results imply that grains are indistinguishable from interstellar medium-like dust in the east side (beta = 2.0 +/- 0.5) but are much smaller in the west side beta = 0.7+0.5 -0.4, indicating differential dust size distribution between the two sides of the disk. Combining the results of near-infrared and submillimeter observations, we conjecture that the spiral arms exist at the upper surface and an asymmetric structure resides in the disk interior. Future observations at centimeter wavelengths and differential polarization imaging in other bands (Y-K) with extreme AO imagers are required to understand how large dust grains form and to further explore the dust distribution in the disk.