Search

Using multifrequency observations, from radio to γ-rays of the blazar Mrk 501, we constructed their corresponding light curves and built periodograms using RobPer and Lomb–Scargle algorithms. Long-term variability was also studied using the power density spectrum and the detrended function analysis. Using the software VARTOOLS Version 1.40, we also computed the analysis of variance, box-least squares and discrete fourier transform. The result of these techniques showed an achromatic periodicity ≲229d. This, combined with the result of pink-color noise in the spectra, led us to propose that the periodicity was produced via a secondary eclipsing supermassive binary black hole orbiting the primary one locked inside the central engine of Mrk 501. We built a relativistic eclipsing model of this phenomenon using Jacobi elliptical functions, finding a periodic relativistic eclipse occurring every ∼224d in all the studied wavebands. This implies that the frequency of the emitted gravitational waves falls slightly above 0.1 mHz, well within the operational range of the upcoming LISA space-based interferometer, and as such, these gravitational waves must be considered as a prime science target for future LISA observations.

The [C ii ] 158 μ m emission line and the underlying far-infrared (FIR) dust continuum are important tracers for studying star formation and kinematic properties of early galaxies. We present a survey of the [C ii ] emission lines and FIR continua of 31 luminous quasars at z > 6.5 using the Atacama Large Millimeter Array (ALMA) and the NOrthern Extended Millimeter Array at sub-arcsec resolution. This survey more than doubles the number of quasars with [C ii ] and FIR observations at these redshifts and enables statistical studies of quasar host galaxies deep into the epoch of reionization. We detect [C ii ] emission in 27 quasar hosts with a luminosity range of L _[C _II _] = (0.3–5.5) × 10 ^9 L _⊙ and detect the FIR continuum of 28 quasar hosts with a luminosity range of L _FIR = (0.5–13.0) × 10 ^12 L _⊙ . Both L _[C _II _] and L _FIR are correlated ( ρ ≃ 0.4) with the quasar bolometric luminosity, albeit with substantial scatter. The quasar hosts detected by ALMA are clearly resolved with a median diameter of ∼5 kpc. About 40% of the quasar host galaxies show a velocity gradient in [C ii ] emission, while the rest show either dispersion-dominated or disturbed kinematics. Basic estimates of the dynamical masses of the rotation-dominated host galaxies yield M _dyn = (0.1–7.5) × 10 ^11 M _⊙ . Considering our findings alongside those of literature studies, we found that the ratio between M _BH and M _dyn is about 10 times higher than that of local M _BH – M _dyn relation on average but with substantial scatter (the ratio difference ranging from ∼0.6 to 60) and large uncertainties.

Theoretical models predict that z ≳ 6 quasars are hosted in the most massive halos of the underlying dark matter distribution and thus would be immersed in protoclusters of galaxies. However, observations report inconclusive results. We investigate the 1.1 proper-Mpc ^2 environment of the z = 7.54 luminous quasar ULAS J1342+0928. We search for Lyman-break galaxy (LBG) candidates using deep imaging from the Hubble Space Telescope (HST) in the Advanced Camera for Surveys (ACS)/F814W, Wide Field Camera 3 (WFC3)/F105W/F125W bands, and Spitzer/Infrared Array Camera at 3.6 and 4.5 μ m. We report a ${z}_{\mathrm{phot}}={7.69}_{-0.23}^{+0.33}$ LBG with mag _F125W = 26.41 at 223 projected proper kpc (pkpc) from the quasar. We find no HST counterpart to one [C ii ] emitter previously found with the Atacama Large millimeter/submillimeter Array (ALMA) at 27 projected pkpc and z _[C II] =7.5341 ± 0.0009 (Venemans et al. 2020 ). We estimate the completeness of our LBG candidates using results from Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey/GOODS deep blank field searches sharing a similar filter setup. We find that >50% of the z ∼ 7.5 Lyman-break galaxies (LBGs) with mag _F125W > 25.5 are missed due to the absence of a filter redward of the Lyman break in F105W, hindering the UV color accuracy of the candidates. We conduct a QSO-LBG clustering analysis revealing a low LBG excess of ${0.46}_{-0.08}^{+1.52}$ in this quasar field, consistent with an average or low-density field. Consequently, this result does not present strong evidence of an LBG overdensity around ULAS J1342+0928. Furthermore, we identify two LBG candidates with a z _phot matching a confirmed z = 6.84 absorber along the line of sight to the quasar. All these galaxy candidates are excellent targets for follow-up observations with JWST and/or ALMA to confirm their redshift and physical properties.

We present results from the Quasar hosts Unveiled by high Angular Resolution Techniques survey studying the circumgalactic medium (CGM) by observing rest-frame UV emission lines Ly α , C iv , and He ii around two radio-loud quasars, 3C 9 ( z = 2.02) and 4C 05.84 ( z = 2.32), using the Keck Cosmic Web Imager. We detect large-scale Ly α nebulae around both quasars with projected diameters ∼100 kpc, with spatially resolved, embedded 15–30 kpc He ii and C iv nebulae around both quasars as well as kinematically distinct He ii and C iv nebulae at a physical separation of ∼15 kpc from both quasars. Observations of H α , H β , and [O iii ] emission using Keck MOSFIRE spectroscopically confirm that the Ly α nebulae extend to companion galaxies and that these quasars are in a protogroup/protocluster environment. We confirm that the He ii and C iv emission is kinematically and spatially coincident with the companion galaxies. We estimate the virial masses of the companion galaxies, their metallicities, and star formation rates, and investigate the sources of ionization. We measure the dynamical mass of the host dark matter halos and estimate that the dark matter halos of these systems will grow to a mass of 2 × 10 ^14 M _⊙ (3C 9) and 2 × 10 ^13 M _⊙ (4C 05.84) by z = 0. The combined CGM and companion galaxies observations indicate Ly α substructure can indicate the presence of companion galaxies in the CGM.

We present Keck Cosmic Web Imager Ly α integral field spectroscopy of the fields surrounding 14 damped Ly α absorbers (DLAs) at z ≈ 2. Of these 14 DLAs, nine have high metallicities ([M/H] > − 0.3), and four of those nine feature a CO-emitting galaxy at an impact parameter ≲30 kpc. Our search reaches median Ly α line flux sensitivities of ∼2 × 10 ^−17 erg s ^−1 cm ^−2 over apertures of ∼6 kpc and out to impact parameters of ∼50 kpc. We recover the Ly α flux of three known Ly α -emitting H i -selected galaxies in our sample. In addition, we find two Ly α emitters at impact parameters of ≈50–70 kpc from the high-metallicity DLA at z ≈ 1.96 toward QSO B0551-366. This field also contains a massive CO-emitting galaxy at an impact parameter of ≈15 kpc. Apart from the field with QSO B0551-366, we do not detect significant Ly α emission in any of the remaining eight high-metallicity DLA fields. Considering the depth of our observations and our ability to recover previously known Ly α emitters, we conclude that H i -selected galaxies associated with high-metallicity DLAs at z ≈ 2 are dusty and therefore might feature low Ly α escape fractions. Our results indicate that complementary approaches—using Ly α , CO, H α , and [C ii ] 158 μ m emission—are necessary to identify the wide range of galaxy types associated with z ≈ 2 DLAs.

Purpose: Estimate of the environment density of giant (with the linear size of about megaparsec) radio structures for galaxies and quasars with steep low-frequency spectra taken from the UTR-2 catalogue. Study of the cosmological evolution of environment density of giant radio sources. Determination of dependence of contribution of radio lobes into the emission of giant sources with respect to their environment density. Design/methodology/approach: We use the sample of sources from the UTR-2 catalogue of extragalactic sources to estimate the environment density for giant sources with steep low-frequency spectra. The selection criteria for the examined objects are the following: 1) the spectral index value is equal or larger than 1; 2) the fl ux density of emission at the frequency of 25 MHz is larger than 10 Jy; 3) the sample sources are optically identifi ed. The value of environment density of examined sources is obtained with the assumption of equality of source jet luminosity (at the synchrotron mechanism of radio emission) and its corresponding kinetic luminosity. The analysis of the estimates of environment densities is made for different classes of the sample objects (for galaxies and quasars with linear steep spectra and with break steep spectra). Findings: The estimates of environment density have been derived for giant radio structures formed by the jets of sources with steep spectrum from the UTR-2 catalogue. On the average, the environment density for the quasar structure (~ 10-28 g/sm3) is lesser than the one for the galaxies (~ 10-27 g/sm3 to ~ 10-26 g/sm3). The larger jet environment density is typical for the galaxies and quasars with the break steep spectra than for those with the linear steep spectra. The inverse power relation of the jet environment density and the source redshift (the cosmological evolution of the jet environment density) has been derived. The contribution of jet-related radio lobes into the emission of sources displays the inverse power relation for the environment density of the corresponding radio structures. Conclusions: The mean values of obtained estimates of environment density of giant jets of radio sources with steep low-frequency spectra indicate the lesser environment density of quasar jets than that for the galaxy jets. Giant radio sources with steep low-frequency spectrum (especially, with break steep spectrum) reveal considerable evolution of environment density of jets. The larger contribution of radio lobes (jets) into the emission of sources corresponds to the lesser environment density of sources taken from the UTR-2 catalogue. It can be due to propagation of jets (surrounded by radio lobes) from powerful radio sources to distances of about megaparsec, until the balance of source’s environment density and extragalactic environment density is reached.

We search for quasi-stellar objects (QSOs) in a wide area of the south ecliptic pole (SEP) field, which has been and will continue to be intensively explored through various space missions. For this purpose, we obtain deep broadband optical images of the SEP field covering an area of ∼14.5 × 14.5 deg ^2 with the Korea Microlensing Telescope Network (KMTNet). The 5 σ detection limits for point sources in the BVRI bands are estimated to be ∼22.59, 22.60, 22.98, and 21.85 mag, respectively. Utilizing data from the Wide-field Infrared Survey Explorer, unobscured QSO candidates are selected among the optically pointlike sources using mid-infrared (MIR) and optical–MIR colors. To refine our selection further and eliminate any contamination not adequately removed by the color-based selection, we perform spectral energy distribution fitting with archival photometric data ranging from optical to MIR. As a result, we identify a total of 2383 unobscured QSO candidates in the SEP field. We also apply a similar method to the north ecliptic pole field using Pan-STARRS data and obtain a similar result of identifying 2427 candidates. The differential number count per area of our QSO candidates is in good agreement with those measured from spectroscopically confirmed ones in other fields. Finally, we compare the results with the literature and discuss how this work will impact future studies, especially upcoming space missions.

We use high-resolution zoom-in cosmological simulations to model outflow triggered by radiation and thermal drivers around the central mass accumulation during direct collapse within the dark matter (DM) halo. The maximal resolution is 1.3 × 10 ^−5 pc, and no restrictions are put on the geometry of the inflow/outflow. The central mass is considered prior to the formation of the supermassive black hole seed at a redshift of z ∼ 15.9 and can constitute either a supermassive star (SMS) of ∼10 ^5 M _⊙ surrounded by a growing accretion disk or a self-gravitating disk. The radiation transfer is modeled using the ray-tracing algorithm. Due to the high accretion rate of ∼1 M _⊙ yr ^−1 determined by the DM halo, accretion is mildly supercritical, resulting in mildly supercritical luminosity that has only a limited effect on the accretion rate, with a duty cycle of ∼0.9. We observe a fast development of hot cavities, which quickly extend into polar funnels and expand dense shells. Within the funnels, fast winds, ∼10 ^3 km s ^−1 , are mass-loaded by the accreting gas. We follow the expanding shells to ∼1 pc, when the shell velocity remains substantially (∼5 times) above the escape speed. The ionization cones formed by the central UV/X-ray completely ionize the cavities. Extrapolating the outflow properties shows that the halo material outside the shell will have difficulty stopping it. We therefore conclude that the expanding wind-driven shell will break out of the central parsec and will reach the halo virial radius. Finally, the anisotropic accretion flow on subparsec scales will attenuate the UV/soft X-rays on the H _2 . Hence, the formation of funnels and powerful outflows around, e.g., SMSs can have interesting observational corollaries.

We report XMM-Newton observations of two examples of an unclassified type of X-ray weak quasars from the Pu et al. survey of X-ray weak quasars in the Chandra archive, SDSS J083116.62+321329.6 at z = 1.797 and SDSS J142339.87+042041.1 at z = 1.702. They do not belong to the known populations of X-ray weak quasars that show broad absorption lines, weak ultraviolet (UV) broad emission lines, or red optical/UV continua. Instead, they display typical quasar UV spectra and spectral energy distributions. In the XMM-Newton observations, both quasars show nominal levels of X-ray emission with typical quasar X-ray spectral shapes (power-law photon indices of ${1.99}_{-0.23}^{+0.27}$ and ${1.86}_{-0.14}^{+0.15}$ ), displaying strong X-ray variability compared to the archival Chandra data (variability factors of ${4.0}_{-1.4}^{+1.6}$ and ${9.0}_{-3.8}^{+7.4}$ in terms of the 2 keV flux density). Simultaneous optical (rest-frame UV) spectra indicate no strong variability compared to the archival spectra. Long-term optical/UV and infrared light curves do not show any substantial variability either. We consider that the X-ray weakness observed in the Chandra data is due to X-ray obscuration from a small-scale dust-free absorber, likely related to accretion-disk winds. Such X-ray weak/absorbed states are probably rare in typical quasars, and thus both targets recovered to X-ray nominal-strength states in the XMM-Newton observations.

We present a long-term and intraday variability study of optical multiwaveband ( U BV RI ) data from the blazar AO 0235+164 collected by various telescopes for ∼44 yr (1975–2019). The blazar was found to be significantly variable over the years in all wave bands with a variation of about 6 mag between its low and active states. The variations in the different wave bands are highly correlated without any time lag. We did not observe any significant trend in color variation with time, but we observed a bluer-when-brighter trend between the B − I color index and the R magnitude. Optical BV R -band spectral energy distributions always show a convex shape. Significant intraday variability was frequently seen in the quasi-simultaneous observations of AO 0235+164 made on 22 nights in the R and V bands by the CASLEO and CAHA telescopes during 1999–2019. We also estimated a central supermassive black hole mass of 7.9 × 10 ^7 M _⊙ by analyzing the broad Mg ii emission line in AO 0235+164's spectrum. We briefly explore the probable physical scenarios responsible for the observed variability.

The radiative and jet power in active galactic nuclei is generated by accretion of material on to supermassive galactic-centre black holes. For quasars, where the radiative power is by definition very high, objects with high radio luminosities form 10 per cent of the population, although it is not clear whether this is a stable phase. Traditionally, quasars with high radio luminosities have been thought to present jets with edge-brightened morphology (Fanaroff-Riley II—FR II) due to the limitations of previous radio surveys (i.e., FRIs were not observed as part of the quasar population). The LOw Frequency ARray (LOFAR) Two-metre Sky Survey (LoTSS) with its unprecedented sensitivity and resolution covering wide sky areas has enabled the first systematic selection and investigation of quasars with core-brightened morphology (Fanaroff-Riley I—FR). We carried out a Very Large Array (VLA) snapshot survey to reveal inner structures of jets in selected quasar candidates; 15 (25 per cent) out of 60 sources show clear inner jet structures that are diagnostic of FRI jets and 13 quasars (∼22 per cent) show extended structures similar to those of FRI jets. Black hole masses and Eddington ratios do not show a clear difference between FRI and FRII quasars. FRII quasars tend to have higher jet powers than FRI quasars. Our results show that the occurrence of FRI jets in powerful radiatively efficient systems is not common, probably mainly due to two factors: galaxy environment and jet power.

A non-negligible fraction of quasars are red at optical wavelengths, indicating (in the majority of cases) that the accretion disc is obscured by a column of dust which extinguishes the shorter-wavelength blue emission. In this paper, we summarize recent work by our group, where we find fundamental differences in the radio properties of SDSS optically-selected red quasars. We also present new analyses, using a consistent color-selected quasar parent sample matched to four radio surveys (FIRST, VLA Stripe 82, VLA COSMOS 3 GHz, and LoTSS DR1) across a frequency range 144 MHz–3 GHz and four orders of magnitude in radio flux. We show that red quasars have enhanced small-scale radio emission (∼kpc) that peaks around the radio-quiet threshold (defined as the ratio of 1.4 GHz luminosity to 6 μm luminosity) across the four radio samples. Exploring the potential mechanisms behind this enhancement, we rule out star-formation and propose either small-scale synchrotron jets, frustrated jets, or dusty winds interacting with the interstellar medium; the latter two scenarios would provide a more direct connection between opacity (dust; gas) and the production of the radio emission. In our future study, using new multi-band uGMRT data, we aim to robustly distinguish between these scenarios.