Search

Dual Active Galactic Nuclei (dual AGNs), a phase in some galaxy mergers during which both central supermassive black holes (SMBHs) are active, are expected to be a key observable stage leading up to SMBH mergers. Constraining the population of dual AGNs in both the nearby and high-z universe has proven to be elusive until very recently. We present a multi-wavelength follow-up campaign to confirm the nature of a sample of 20 candidate dual AGNs at cosmic noon (z~2) from the VODKA sample. Through a combination of Hubble Space Telescope (HST) and Very Large Array (VLA) imaging, we refute the possibility of gravitational lensing in all but one target. We find evidence of dual AGNs in four systems, while seven exhibit single AGN in galaxy pairs, either through strong radio emission or ancillary emission line data. The remaining systems are either confirmed as quasar-star superpositions (six) or non-lensed pairs (two) that require further investigations to establish AGN activity. Among the systems with radio detections, we find a variety of radio spectral slopes and UV/optical colors suggesting that our sample contains a range of AGN properties, from obscured radio-quiet objects to those with powerful synchrotron-emitting jets. This study presents one of the largest dedicated multi-wavelength follow-up campaigns to date searching for dual AGNs at high redshift. We confirm several of the highest-z systems at small physical separations, thus representing some of the most evolved dual AGN systems at the epoch of peak quasar activity known to date., Comment: 26 pages, 8 figures, submitted to ApJ, includes appendix of 9 additional figures and 18 tables

We present Gemini/GMOS and HST/STIS optical spectra for 27 dual quasar candidates selected based on their variability-induced astrometric noise or double detections in Gaia (the VODKA project). From this follow-up, we spectroscopically identify 10 star superpositions and 8 dual/lensed quasars. Among the remaining targets, 2 are likely dual/lensed quasars based on additional radio imaging, while the rest are quasars with unknown companions. Notably, WISE J1649+0812 is a newly confirmed dual quasar with a projected separation of 5 kpc at $z=1.39$ and a significant velocity offset of 183$\pm$76 km/s, highlighting the utility of narrow emission lines in identifying genuine dual quasars. Without prior photometric or spectroscopic selection, we find the star contamination rate to be 37-63%, while the dual/lensed quasar fraction is $\gtrsim$ 30% in the follow-up VODKA sample. However, when combined with existing unresolved spectra and spatially-resolved two-band color cuts, the dual/lensed quasar fraction can be increased to $\gtrsim$ 67%. High signal-to-noise ratio spectra ($\gtrsim$ 20 per spectral element) with adequate spectral resolution ($R \gtrsim$ 1000) are essential for identifying faint absorption lines in foreground stars and detecting dual quasars through velocity offsets., Comment: 14 pages, 6 figures, 4 tables, submitted to ApJ, comments are welcome

We present JWST MIRI/NIRSpec observations of the extremely red quasar SDSS J165202.64+172852.3 at z~3, one of the most luminous quasars known to date, driving powerful outflows and hosting a clumpy starburst, amidst several interacting companions. We estimate the black hole (BH) mass of the system based on the broad H$\alpha$ and H$\beta$ lines, as well as the Pa$\beta$ emission in the IR and MgII in the UV. We recover a very broad range of mass estimates, with constraints ranging between log $M_{\rm BH}$=9 and 10.1, which is exacerbated if imposing a uniform BLR geometry at all wavelengths. Several factors may contribute to the large spread: measurement uncertainties (insufficient sensitivity to detect the broadest component of the faint Pa$\beta$ line, spectral blending, ambiguities in the broad/narrow component distinction), lack of virial equilibrium (in a system characterised by powerful outflows and rapid accretion), and uncertainties on the luminosity-inferred size of the broad line region, a.o. given central dust obscuration. We constrain the stellar mass via SED fitting, suggesting the host to be extremely massive at $10^{12.8\pm 0.5} M_\odot$ - ~2 dex above the characteristic mass of the Schechter fit to the z=3 stellar mass function. Notably, J1652's central BH might be interpreted as being either undermassive, overmassive, or in line with the BH mass-stellar mass relation, depending on the choice of assumptions. The recovered Eddington ratio varies accordingly, but exceeds 10% in any case. We put our results into context by providing an extensive overview and discussion of recent literature results and their associated assumptions. Our findings provide an important demonstration of the uncertainties inherent in virial BH mass estimates, which are of particular relevance in the JWST era given the growing number of studies on rapidly accreting quasars at high redshift., Comment: 15 pages, 7 figures, submitted to A&A, comments welcome

The search for dual supermassive black holes (SMBHs) is of immense interest in modern astrophysics. Galaxy mergers may be an important route to fuel and to produce SMBH pairs. Actively accreting SMBH pairs can be observed as a dual quasar, which are vital probes of SMBH growth. Gaia observations have enabled a novel technique to systematically search for such dual quasars at previously unreachable sub-kpc scales, based on the small jitters of the light centroid as the two quasars vary stochastically. Here we present the first detailed study of a 0.46'', 3.8 kpc separation, VODKA-selected dual quasar, J0749+2255, at $z=2.17$ using JWST/NIRSpec integral field unit spectroscopy. This is one of the most distant, small separation dual quasars identified today. Dual quasars at cosmic noon are not well characterized. We detect the faint ionized gas of the host galaxy, best traced by the narrow \ha\ emission. Line ratio diagnostics show a mix of ionization from the two quasars and intense star formation. The spatially-resolved spectra of the two quasars suggest that they have very similar black hole properties (two $M_{BH}\sim 10^9\ \textrm{M}_{\odot}$ with large Eddington ratio reaching $L/L_{Edd}\sim0.2$) hinting at the possible synchronized growth and accretion from the same gas supply. Surprisingly, the ionized gas kinematics suggest an extended, rotating disk rather than a disturbed system that would be expected in a major gas-rich galaxy merger. While it is unclear if J0749+2255 is representative of the dual quasar evolution, the observations with JWST revealed a major puzzle. It would be interesting to see what observations of other dual quasars will show., Comment: 21 pages, 9 figures, submitted

Dual quasars, two active supermassive black holes at galactic scales, represent crucial objects for studying the impact of galaxy mergers and quasar activity on the star formation rate (SFR) within their host galaxies, particularly at cosmic noon when SFR peaks. We present JWST/MIRI mid-infrared integral field spectroscopy of J074922.96+225511.7, a dual quasar with a projected separation of 3.8 kilo-parsec at a redshift $z$ of 2.17. We detect spatially extended [Fe II] 5.34$\rm \mu$m and polycyclic aromatic hydrocarbon (PAH) 3.3$\mu$m emissions from the star formation activity in its host galaxy. We derive the SFR of 10$^{3.0\pm0.2}$ M$_{\odot}$ yr$^{-1}$ using PAH 3.3$\mu$m, which is five times higher than that derived from the cutoff luminosity of the infrared luminosity function for galaxies at $z\sim2$. While the SFR of J0749+2255 agrees with that of star-forming galaxies of comparable stellar mass at the same redshifts, its molecular gas content falls short of expectations based on the molecular Kennicutt-Schmidt law. This discrepancy may result from molecular gas depletion due to the longer elevated stage of star formation, even after the molecular gas reservoir is depleted. We do not observe any quasar-driven outflow that impacts PAH and [Fe II] in the host galaxy based on the spatially resolved maps. From the expected flux in PAH-based star formation, the [Fe II] line likely originates from the star-forming regions in the host galaxy. Our study highlights the stardust nature of J0749+2255, indicating a potential connection between the dual quasar phase and intense star formation activities., Comment: 12 pages, 7 figures, 1 table, submitted to ApJ, comments are welcome

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/

We study dual AGN host galaxy morphologies at $z=2$ using the ASTRID simulation, selecting black hole (BH) pairs with small separation ($\Delta r<30\rm{kpc}$), high mass ($M_{\text{BH,12}}>10^7M_\odot$), and luminosity ($L_{\text{bol,12}}>10^{43}\rm{erg/s}$). We kinematically decompose (using MORDOR) $\sim1000$ dual AGN hosts into standard components - a `disk' (thin and thick disk, pseudo-bulge) and 'bulge' (bulge and halo) and define disk-dominated galaxies by the disk-to-total $D/T\geq0.5$. In ASTRID, $60.9\pm2.1\%$ of dual AGN hosts (independent of separation) are disk-dominated, with the $D/T$ distribution peaking at $\sim0.7$. Notably, hosts of BH pairs have similar morphologies (most either both disk or bulge-dominated). In dual-AGN hosts, the $D/T$ increases from $\sim17\% $ at $M_{\rm *}\sim 10^{9} M_{\odot}$ to $ 64\% $ for $M_{\rm *} \sim 10^{11.5} M_{\odot}$, and the pseudo-bulge is the dominant component of the disk fraction at the high mass end. Moreover, dual AGN hosts exhibit a higher fraction of disk/large pseudo-bulge than single-AGN hosts. The Disk-to-Total ratio is approximately constant with BH mass or AGN luminosity. We also create mock images of dual AGN host galaxies, employing morphological fitting software Statmorph to calculate morphological parameters and compare them with our kinematic decomposition results. Around $83.3\pm2.4\%$ of galaxies display disk-like profiles, of which $\sim60.7\pm2.2\%$ are kinematically confirmed as disks. Se\'rsic indices and half-mass radii of dual AGN host galaxies align with observational measurements from HST at $z\sim2$. Around $34\%$ are identified as mergers from the $\text{Gini}-M_{20}$ relation. We find two dual AGN hosted by galaxies that exhibit disk-like se\'rsic index $n_{12}<1$ and $(D/T)_{12}>0.5$, which are in remarkable agreement with properties of recently discovered dual quasars in disk galaxies at $z\sim 2$., Comment: 15 pages, 12 figures, submitted to the Open Journal of Astrophysics

Lightcurves of many classical novae deviate from the canonical "fast rise - smooth decline" pattern and display complex variability behavior. We present the first TESS-space-photometry-based investigation of this phenomenon. We use TESS Sector 41 full-frame images to extract a lightcurve of the slow Galactic nova V606 Vul that erupted nine days prior to the start of the TESS observations. The lightcurve covers the first of two major peaks of V606 Vul that was reached 19 days after the start of the eruption. The nova reached its brightest visual magnitude V=9.9 in its second peak 64 days after the eruption onset, following the completion of Sector 41 observations. To increase the confidence level of the extracted lightcurve, we performed the analysis using four different codes implementing the aperture photometry (Lightkurve, VaST) and image subtraction (TESSreduce, tequila_shots) and find good agreement between them. We performed ground-based photometric and spectroscopic monitoring to complement the TESS data. The TESS lightcurve reveals two features: periodic variations (0.12771d, 0.01mag average peak-to-peak amplitude) that disappeared when the source was within 1mag of peak optical brightness and a series of isolated mini-flares (with peak-to-peak amplitudes of up to 0.5mag) appearing at seemingly random times. We interpret the periodic variations as the result of azimuthal asymmetry of the photosphere engulfing the nova-hosting binary that was distorted by and rotating with the binary. Whereas we use spectra to associate the two major peaks in the nova lightcurve with distinct episodes of mass ejection, the origin of mini-flares remains elusive., Comment: 31 pages, 10 figures, submitted to ApJ; comments welcome

Dual and off-nucleus active supermassive black holes are expected to be common in the hierarchical structure formation paradigm, but their identification at parsec scales remains a challenge due to strict angular resolution requirements. We conduct a systematic study using the Very Long Baseline Array (VLBA) to examine 23 radio-bright candidate dual and off-nucleus quasars. The targets are selected by a novel astrometric technique ("varstrometry") from Gaia, aiming to identify dual or off-nucleus quasars at (sub)kilo-parsec scales. Among these quasars, 8 exhibit either multiple radio components or significant (>3$\sigma$) positional offsets between the VLBA and Gaia positions. The radio emission from the three candidates which exhibit multiple radio components is likely to originate from small-scale jets based on their morphology. Among the remaining five candidates with significant VLBA-Gaia offsets, three are identified as potential dual quasars at parsec scales, one is likely attributed to small-scale jets, and the origin of the last candidate remains unclear. We explore alternative explanations for the observed VLBA-Gaia offsets. We find no evidence for optical jets at kilo-parsec scales, nor any contamination to Gaia astrometric noise from the host galaxy; misaligned coordinate systems are unlikely to account for our offsets. Our study highlights the promise of the varstrometry technique in discovering candidate dual or off-nucleus quasars and emphasizes the need for further confirmation and investigation to validate and understand these intriguing candidates., Comment: 17 pages, 8 figures, 4 tables, accepted to ApJ

The hierarchical structure formation paradigm predicts the formation of pairs of supermassive black holes in merging galaxies. When both (or one) members of the SMBH pair are unobscured AGNs, the system can be identified as a dual (or offset) AGN. Quantifying the abundance of these AGN pairs as functions of separation, redshift and host properties is crucial to understanding SMBH formation and AGN fueling in the broad context of galaxy formation. The High Latitude Wide Area Survey with Roman, with its unprecedented combination of sensitivity, spatial resolution, area and NIR wavelength coverage, will revolutionize the study of galactic-scale environments of SMBH pairs. This white paper summarizes the science opportunities and technical requirements on the discovery and characterization of SMBH pairs down to galactic scales (i.e., less than tens of kpc) over broad ranges of redshift (11E42 erg/s)., Comment: Roman Core Community Survey White Paper, focusing on the High Latitude Wide Area Survey

We report Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) deep IR (F160W) imaging of SDSS J1608+2716. This system, located at a redshift of z=2.575, was recently reported as a triple quasar candidate with subarcsecond separations ($\sim0.25''$) based on selection from Gaia astrometry and follow-up Keck adaptive optics-assisted integral field unit spectroscopy. Our new HST deep IR imaging reveals the presence of a fourth point-like component located $\sim0.9''$ away from the triple system. Additionally, we detect an edge-on disk galaxy located in between the four point sources. The entire system exhibits a characteristic cusp structure in the context of strong gravitational lensing, and the observed image configuration can be successfully reproduced using a lens model based on a singular isothermal ellipsoid mass profile. These findings indicate that this system is a quadruply lensed quasar. Our results highlight the challenges associated with identifying dual/multiple quasars on $\sim$kpc scales at high redshifts, and emphasize the crucial role of deep, high-resolution IR imaging in robustly confirming such systems., Comment: 10 pages, 3 figures, accepted by ApJL

Dual quasars at small physical separations are an important precursor phase of galaxy mergers, ultimately leading to the coalescence of the two supermassive black holes. Starting from a sample of dual/lensed quasar candidates discovered using astrometric jitter in Gaia data, we present a pilot case study of one of the most promising yet puzzling candidate dual quasars at cosmic noon (z$\sim$1.8). Using multi-wavelength imaging and spectroscopy from X-ray to radio, we test whether the J0823+2418 system is two individual quasars in a bound pair at separation$\sim$0.64'', or instead a single quasar being gravitationally lensed by a foreground galaxy. We find consistent flux ratios ($\sim$1.25-1.45) between the two sources in optical, NIR, UV, and radio, and thus similar spectral energy distributions, suggesting a strong lensing scenario. However, differences in the radio spectral index, as well as changing X-ray fluxes, hint at either a dual quasar with otherwise nearly identical properties, or perhaps lensing-based time lag of $\sim$3 days paired with intrinsic variability. We find with lens mass modeling that the relative NIR positions and magnitudes of the two quasars and a marginally detected central galaxy are consistent with strong lensing. Archival SDSS spectra likewise suggest a foreground absorber via Mg II absorption lines. We conclude that J0823+2418 is likely a lensed quasar, and therefore that the VODKA sample contains a population of these lensed systems (perhaps as high as 50%) as well as dual quasars., Comment: 25 pages, 14 figures, submitted to ApJ

Galaxy mergers produce pairs of supermassive black holes (SMBHs), which may be witnessed as dual quasars if both SMBHs are rapidly accreting. The kiloparsec (kpc)-scale separation represents a physical regime sufficiently close for merger-induced effects to be important yet wide enough to be directly resolvable with the facilities currently available. Whereas many kpc-scale dual active galactic nuclei--the low-luminosity counterparts of quasars--have been observed in low-redshift mergers, no unambiguous dual quasar is known at cosmic noon ($z\approx2$), the peak of global star formation and quasar activity. Here we report multiwavelength observations of SDSS J0749+2255 as a kpc-scale, dual-quasar system hosted by a galaxy merger at cosmic noon ($z=2.17$). We discover extended host galaxies associated with the much brighter compact quasar nuclei (separated by 0.46" or 3.8 kpc) and low-surface-brightness tidal features as evidence for galactic interactions. Unlike its low-redshift and low-luminosity counterparts, SDSS J0749+2255 is hosted by massive compact disk-dominated galaxies. The apparent lack of stellar bulges and the fact that SDSS J0749+2255 already follows the local SMBH mass--host stellar mass relation, suggest that at least some SMBHs may have formed before their host stellar bulges. While still at kpc-scale separations where the host-galaxy gravitational potential dominates, the two SMBHs may evolve into a gravitationally bound binary system in around 0.22 Gyr., Comment: 88 pages, 18 figures, 8 tables; published

The statistics of galactic-scale quasar pairs can elucidate our understanding of the dynamical evolution of supermassive black hole (SMBH) pairs, the duty cycles of quasar activity in mergers, or even the nature of dark matter, but have been challenging to measure at cosmic noon, the prime epoch of massive galaxy and SMBH formation. Here we measure a double quasar fraction of $\sim 6.2\pm0.5\times 10^{-4}$ integrated over $\sim 0.3-3$ arcsec separations (projected physical separations of $\sim 3-30\,{\rm kpc}$ at $z\sim 2$) in luminous ($L_{\rm bol}>10^{45.8}\,{\rm erg\,s^{-1}}$) unobscured quasars at $1.5

We study the optical $gri$ photometric variability of a sample of 190 quasars within the SDSS Stripe 82 region that have long-term photometric coverage during $\sim 1998-2020$ with SDSS, PanSTARRS-1, the Dark Energy Survey, and dedicated follow-up monitoring with Blanco 4m/DECam. With on average $\sim 200$ nightly epochs per quasar per filter band, we improve the parameter constraints from a Damped Random Walk (DRW) model fit to the light curves over previous studies with 10-15 yr baselines and $\lesssim 100$ epochs. We find that the average damping timescale $\tau_{\rm DRW}$ continues to rise with increased baseline, reaching a median value of $\sim 750$ days ($g$ band) in the rest-frame of these quasars using the 20-yr light curves. Some quasars may have gradual, long-term trends in their light curves, suggesting that either the DRW fit requires very long baselines to converge, or that the underlying variability is more complex than a single DRW process for these quasars. Using a subset of quasars with better-constrained $\tau_{\rm DRW}$ (less than 20\% of the baseline), we confirm a weak wavelength dependence of $\tau_{\rm DRW}\propto \lambda^{0.51\pm0.20}$. We further quantify optical variability of these quasars over days to decades timescales using structure function (SF) and power spectrum density (PSD) analyses. The SF and PSD measurements qualitatively confirm the measured (hundreds of days) damping timescales from the DRW fits. However, the ensemble PSD is steeper than that of a DRW on timescales less than $\sim$ a month for these luminous quasars, and this second break point correlates with the longer DRW damping timescale., Comment: 20 pages, 24 figures

We present a sample of 706, $z < 1.5$ active galactic nuclei (AGNs) selected from optical photometric variability in three of the Dark Energy Survey (DES) deep fields (E2, C3, and X3) over an area of 4.64 deg$^2$. We construct light curves using difference imaging aperture photometry for resolved sources and non-difference imaging PSF photometry for unresolved sources, respectively, and characterize the variability significance. Our DES light curves have a mean cadence of 7 days, a 6 year baseline, and a single-epoch imaging depth of up to $g \sim 24.5$. Using spectral energy distribution (SED) fitting, we find 26 out of total 706 variable galaxies are consistent with dwarf galaxies with a reliable stellar mass estimate ($M_{\ast}<10^{9.5}\ M_\odot$; median photometric redshift of 0.9). We were able to constrain rapid characteristic variability timescales ($\sim$ weeks) using the DES light curves in 15 dwarf AGN candidates (a subset of our variable AGN candidates) at a median photometric redshift of 0.4. This rapid variability is consistent with their low black hole masses. We confirm the low-mass AGN nature of one source with a high S/N optical spectrum. We publish our catalog, optical light curves, and supplementary data, such as X-ray properties and optical spectra, when available. We measure a variable AGN fraction versus stellar mass and compare to results from a forward model. This work demonstrates the feasibility of optical variability to identify AGNs with lower black hole masses in deep fields, which may be more "pristine" analogs of supermassive black hole seeds., Comment: 21 pages, 22 figures incl. 3 appendices; accepted for publication in MNRAS

We present results on a multi-wavelength analysis of SDSS J025214.67-002813.7, a system which has been previously classified as a binary AGN candidate based on periodic signals detected in the optical light curves. We use available radio-X-ray observations of the system to investigate the true accretion nature. Analyzing new observations from XMM-Newton and NuSTAR, we characterize the X-ray emission and search for evidence of circumbinary accretion. Although the 0.5-10 keV spectrum shows evidence of an additional soft emission component, possibly due to extended emission from hot nuclear gas, we find the spectral shape consistent with a single AGN. Compiling a full multi-wavelength SED, we also search for signs of circumbinary accretion, such as a "notch" in the continuum due to the presence of minidisks. We find that the radio-optical emission agrees with the SED of a standard, radio-quiet, AGN, however there is a large deficit in emission blueward of ~1400 A. Although this deficit in emission can plausibly be attributed to a binary AGN system, we find that the SED of SDSS J0252-0028 is better explained by emission from a reddened, single AGN. However, future studies on the expected hard X-ray emission associated with binary AGN (especially in the unequal-mass regime), will allow for more rigorous analyses of the binary AGN hypothesis., Comment: Submitted to ApJ. 11 pages, 3 figures