Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “QUASARS”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,153 records · Page 64

STRIDES: Spectroscopic and photometric characterization of the environment and effects of mass along the line of sight to the gravitational lenses DES J0408–5354 and WGD 2038–4008

In time-delay cosmography, three of the key ingredients are (1) determining the velocity dispersion of the lensing galaxy, (2) identifying galaxies and groups along the line of sight with sufficient proximity and mass to be included in the mass model, and (3) estimating the external convergence κ_ext from less massive structures that are not included in the mass model. We present results on all three of these ingredients for two time-delay lensed quad quasar systems, DES J0408–5354 and WGD 2038–4008 . We use the Gemini, Magellan, and VLT telescopes to obtain spectra to both measure the stellar velocity dispersions of the main lensing galaxies and to identify the line-of-sight galaxies in these systems. Next, we identify 10 groups in DES J0408–5354 and two groups in WGD 2038–4008 using a group-finding algorithm. We then identify the most significant galaxy and galaxy-group perturbers using the ‘flexion shift’ criterion. We determine the probability distribution function of the external convergence κ_ext for both of these systems based on our spectroscopy and on the DES-only multiband wide-field observations. Using weighted galaxy counts, calibrated based on the Millennium Simulation, we find that DES J0408–5354 is located in a significantly underdense environment, leading to a tight (width |$\sim 3{{\ \rm per\ cent}}$|⁠), negative-value κ_ext distribution. On the other hand, WGD 2038–4008 is located in an environment of close to unit density, and its low source redshift results in a much tighter κ_ext of |$\sim 1{{\ \rm per\ cent}}$|⁠, as long as no external shear constraints are imposed.

WGD 2038-4008↗

OzDES Reverberation Mapping Programme: the first Mg ii lags from 5 yr of monitoring

ABSTRACT Reverberation mapping is a robust method to measure the masses of supermassive black holes outside of the local Universe. Measurements of the radius–luminosity (R−L) relation using the Mg ii emission line are critical for determining these masses near the peak of quasar activity at z ≈ 1−2, and for calibrating secondary mass estimators based on Mg ii that can be applied to large samples with only single-epoch spectroscopy. We present the first nine Mg ii lags from our 5-yr Australian Dark Energy Survey reverberation mapping programme, which substantially improves the number and quality of Mg ii lag measurements. As the Mg ii feature is somewhat blended with iron emission, we model and subtract both the continuum and iron contamination from the multiepoch spectra before analysing the Mg ii line. We also develop a new method of quantifying correlated spectroscopic calibration errors based on our numerous, contemporaneous observations of F-stars. The lag measurements for seven of our nine sources are consistent with both the H β and Mg ii R−L relations reported by previous studies. Our simulations verify the lag reliability of our nine measurements, and we estimate that the median false positive rate of the lag measurements is $4{{\ \rm per\ cent}}$.

79 ASTRONOMY AND ASTROPHYSICS↗

A comparative study of the physical properties for a representative sample of Narrow and Broad-line Seyfert galaxies

ABSTRACT We present a comparative study of the physical properties of a homogeneous sample of 144 Narrow line Seyfert 1 (NLSy1) and 117 Broad-line Seyfert 1 (BLSy1) galaxies. These two samples are in a similar luminosity and redshift range and have optical spectra available in the 16th data release of Sloan Digital Sky Survey (SDSS-DR16) and X-ray spectra in either XMM-NEWTON or ROSAT. Direct correlation analysis and a principal component analysis (PCA) have been performed using ten observational and physical parameters obtained by fitting the optical spectra and the soft X-ray photon indices as another parameter. We confirm that the established correlations for the general quasar population hold for both types of galaxies in this sample despite significant differences in the physical properties. We characterize the sample also using the line shape parameters, namely the asymmetry and kurtosis indices. We find that the fraction of NLSy1 galaxies showing outflow signatures, characterized by blue asymmetries, is higher by a factor of about 3 compared to the corresponding fraction in BLSy1 galaxies. The presence of high iron content in the broad-line region of NLSy1 galaxies in conjunction with higher Eddington ratios can be the possible reason behind this phenomenon. We also explore the possibility of using asymmetry in the emission lines as a tracer of outflows in the inner regions of Active Galactic Nuclei. The PCA results point to the NLSy1 and BLSy1 galaxies occupying different parameter spaces, which challenges the notion that NLSy1 galaxies are a subclass of BLSy1 galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

Measuring the thermal and ionization state of the low- z IGM using likelihood free inference

ABSTRACT We present a new approach to measure the power-law temperature density relationship $T=T_0 (\rho/ \bar{\rho })^{\gamma -1}$ and the UV background photoionization rate $\Gamma _{{{{\rm H\, {\small I}}}}{}}$ of the intergalactic medium (IGM) based on the Voigt profile decomposition of the Ly α forest into a set of discrete absorption lines with Doppler parameter b and the neutral hydrogen column density $N_{\rm H\, {\small I}}$. Previous work demonstrated that the shape of the $b-N_{{{{\rm H\, {\small I}}}}{}}$ distribution is sensitive to the IGM thermal parameters T0 and γ, whereas our new inference algorithm also takes into account the normalization of the distribution, i.e. the line-density dN/dz, and we demonstrate that precise constraints can also be obtained on $\Gamma _{{{{\rm H\, {\small I}}}}{}}$. We use density-estimation likelihood-free inference (DELFI) to emulate the dependence of the $b-N_{{{{\rm H\, {\small I}}}}{}}$ distribution on IGM parameters trained on an ensemble of 624 nyx hydrodynamical simulations at z = 0.1, which we combine with a Gaussian process emulator of the normalization. To demonstrate the efficacy of this approach, we generate hundreds of realizations of realistic mock HST/COS data sets, each comprising 34 quasar sightlines, and forward model the noise and resolution to match the real data. We use this large ensemble of mocks to extensively test our inference and empirically demonstrate that our posterior distributions are robust. Our analysis shows that by applying our new approach to existing Ly α forest spectra at z ≃ 0.1, one can measure the thermal and ionization state of the IGM with very high precision ($\sigma _{\log T_0} \sim 0.08$ dex, σγ ∼ 0.06, and $\sigma _{\log \Gamma _{{{{\rm H\, {\small I}}}}{}}} \sim 0.07$ dex).

79 ASTRONOMY AND ASTROPHYSICS↗

Efficient analysis routines for single and double peaked Type 2 AGN spectra

Driven by the imminent need to rapidly process and classify millions of AGN spectra drawn from next generation astronomical facilities, we present a spectral fitting routine for Type 2 AGN spectra optimized for high volume processing, using the quasar spectral fitting library (qsfit). We analyse an optically selected sample of 813 luminous Type 2 AGN spectra at z < 0.83 from the Sloan Digital Sky Survey (SDSS) to qualify its performance. We report a median narrow line H α/H β Balmer decrement of 4.5 ± 0.8, alluding to the presence of dust in the narrow line region (NLR). We publish a specialized QSFIT fitting routine for high signal-to-noise ratio spectra and general fitting routine for double peaked Type 2 AGN spectra applied on a subsample of 45 spectra from our parent sample. We report a median red and blue peak velocity separation of 390 ± 60 kms −1 . No trend is found for red or blue peaks to exhibit systematically different luminosity or ionization properties. Emission line diagnostics show that the double peaks in all sources are illuminated by an AGN-powered ionizing continuum. Finally, we examine the morphology of host galaxies of our double peaked sample. We find double peaked Type 2 AGN reside in merging systems at a comparable frequency to single peaked AGN. This suggests that the double peaked AGN phenomenon is likely to have a bi-conical outflow origin in the majority of cases. We publicly release the code used for spectral analysis and produced catalogues used in this work.

79 ASTRONOMY AND ASTROPHYSICS↗

The evolution of the heaviest supermassive black holes in jetted AGNs

ABSTRACT We present the space density evolution, from z = 1.5 up to z = 5.5, of the most massive ($M \ge 10^9\, \mathrm{M}_{\odot }$) black holes hosted in jetted active galactic nuclei (AGNs). The analysis is based on a sample of 380 luminosity-selected (λL1350 ≥ 1046 erg s−1 and $P_{5\, \text{GHz}}\ge 10^{27}$ W Hz−1) flat spectrum radio quasars (FSRQs) obtained from the Cosmic Lens All Sky Survey (CLASS). These sources are known to be face-on jetted AGNs (i.e. blazars) and can be exploited to infer the abundance of all the (misaligned) jetted AGNs, using a geometrical argument. We then compare the space density of the most massive supermassive black holes hosted in jetted AGNs with those present in the total population (mostly composed by non-jetted AGNs). We find that the space density has a peak at z ∼ 3, which is significantly larger than the value observed in the total AGN population with similar optical/UV luminosities (z ∼ 2.2), but not as extreme as the value previously inferred from X-ray-selected blazars (z ≳ 4). The jetted fraction (jetted AGNs/total AGNs) is overall consistent with the estimates in the local Universe (10–20 per cent) and at high redshift, assuming Lorentz bulk factors Γ ≈ 5. Finally, we find a marginal decrease in the jetted fraction at high redshifts (by a factor of ∼2). All these evidences point towards a different evolutionary path in the jetted AGNs compared to the total AGN population.

(galaxies:) quasars: supermassive black holes↗

All-purpose, all-sky photometric redshifts for the Legacy Imaging Surveys Data Release 8

ABSTRACT In this paper, we present photometric redshift (photo-z) estimates for the Dark Energy Spectroscopic Instrument (DESI) Legacy Imaging Surveys, currently the most sensitive optical survey covering the majority of the extragalactic sky. Our photo-z methodology is based on a machine-learning approach, using sparse Gaussian processes augmented with Gaussian mixture models (GMMs) that allow regions of parameter space to be identified and trained separately in a purely data-driven way. The same GMMs are also used to calculate cost-sensitive learning weights that mitigate biases in the spectroscopic training sample. By design, this approach aims to produce reliable and unbiased predictions for all parts of the parameter space present in wide area surveys. Compared to previous literature estimates using the same underlying photometry, our photo-zs are significantly less biased and more accurate at z > 1, with negligible loss in precision or reliability for resolved galaxies at z < 1. Our photo-z estimates offer accurate predictions for rare high-value populations within the parent sample, including optically selected quasars at the highest redshifts (z > 6), as well as X-ray or radio continuum selected populations across a broad range of flux (densities) and redshift. Deriving photo-z estimates for the full Legacy Imaging Surveys Data Release 8, the catalogues provided in this work offer photo-z estimates predicted to be of high quality for ≳9 × 108 galaxies over ${\sim}19\, 400\, \text{deg}^{2}$ and spanning 0 < z ≲ 7, offering one of the most extensive samples of redshift estimates ever produced.

79 ASTRONOMY AND ASTROPHYSICS↗

Measurements of the z > 5 Lyman-α forest flux autocorrelation functions from the extended XQR-30 data set

We present the first observational measurements of the Lyman-α (Ly α) forest flux autocorrelation functions in ten redshift bins from 5.1 ≤ z ≤ 6.0. We use a sample of 35 quasar sightlines at z > 5.7 from the extended XQR-30 data set; these data have signal-to-noise ratios of >20 per spectral pixel. We carefully account for systematic errors in continuum reconstruction, instrumentation, and contamination by damped Ly α systems. With these measurements, we introduce software tools to generate autocorrelation function measurements from any simulation. Our measurements of the smallest bin of the autocorrelation function increase with redshift when normalizing by the mean flux, $\langle{F}\rangle$. This increase may come from decreasing $\langle{F}\rangle$ or increasing mean free path of hydrogen-ionizing photons, λmfp. Recent work has shown that the autocorrelation function from simulations at z > 5 is sensitive to λmfp, a quantity that contains vital information on the ending of reionization. For an initial comparison, we show our autocorrelation measurements with simulation models for recently measured λmfp values and find good agreements. Further work in modelling and understanding the covariance matrices of the data is necessary to get robust measurements of λmfp from this data.

79 ASTRONOMY AND ASTROPHYSICS↗

Measurements of the thermal and ionization state of the intergalactic medium during the cosmic afternoon

We perform the first measurement of the thermal and ionization state of the intergalactic medium (IGM) across 0.9 < z < 1.5 using 301 Ly α absorption lines fitted from 12 archival Hubble Space Telescope Space Telescope Imaging Spectrograph quasar spectra. We employ the machine-learning-based inference method that uses joint Doppler parameter–column density (⁠b-N HI ⁠) distributions obtained from Ly α forest decomposition. Our results show that the Γ HI photoionization rates, ⁠, agree with recent ultraviolet background synthesis models, with log(Γ HI /s -1 ) = $-11.79^{+0.18}_{-0.15}$, $-11.98^{+0.09}_{-0.09}$⁠, and $-12.32^{+0.10}_{-0.12}$⁠, at z = 1.4, 1.2, and 1, respectively. We obtain the IGM temperature at the mean density, T 0 ⁠, and the adiabatic index, γ⁠, as [log(T 0 /K), γ] = $[4.13^{+0.12}_{-0.10}, 1.34^{+0.10}_{-0.15}]$, $[3.79^{+0.11}_{-0.11}, 1.70^{+0.09}_{-0.09}]$, and $[4.12^{+0.15}_{-0.25}, 1.34^{+0.21}_{-0.26}]$ at z = 1.4⁠, 1.2, and 1. Our measurements of T 0 at z = 1.4 and 1.2 are consistent with the trend predicted from previous z < 3 temperature measurements and theoretical expectations, where the IGM cools down after $He\tiny{II}$ reionization in the absence of any non-standard heating. However, our T 0 measurement at z = 1 unexpectedly high IGM temperature. Given the relatively large uncertainty in these measurements, where σ T$_0$ ~ 5000 K, mostly emanating from the limited size of our data set, we cannot conclude whether the IGM cools down as expected. Lastly, we generate mock data sets to test the constraining power of future measurement with larger data sets. The results demonstrate that, with redshift path-length Δz ~ 2 for each redshift bin, three times the current data set, we can constrain the T 0 of IGM within 1500 K, which would be sufficient to constrain the IGM thermal history at z < 1.5 conclusively.

79 ASTRONOMY AND ASTROPHYSICS↗

Stellar population and metal production in AGN discs

ABSTRACT As gravitational wave detections increase the number of observed compact binaries (consisting of neutron stars or blacks), we begin to probe the different conditions producing these binaries. Most studies of compact remnant formation focus either on stellar collapse from the evolution of field binary stars in gas-free environments or on the formation of stars in clusters where dynamical interactions capture the compact objects, forming binaries. But a third scenario exists. In this paper, we study the fate of massive stars formed, accrete gas, and evolve in the dense discs surrounding supermassive black holes. We calculate the explosions produced and compact objects formed by the collapse of these massive stars. Nucleosynthetic yields may provide an ideal, directly observable, diagnostic of the formation and fate of these stars in active galactic nuclei. We present a first study of the explosive yields from these stars, comparing these yields with the observed nucleosynthetic signatures in the discs around supermassive stars with quasars. We show that, even though these stars tend to form black holes, their rapid rotation leads to discs that can eject a considerable amount of iron during the collapse of the star. The nucleosynthetic yields from these stars can produce constraints on the number of systems formed in this manner, but further work is needed to exploit variations from the initial models presented in this paper.

79 ASTRONOMY AND ASTROPHYSICS↗

A New Sample of Gamma-Ray Emitting Jetted Active Galactic Nuclei—Preliminary Results

We are compiling a new list of gamma-ray jetted active galactic nuclei (AGN), starting from the fourth catalog of point sources of the Fermi Large Area Telescope (LAT). Our aim is to prepare a list of jetted AGN with known redshifts and classifications to be used to calibrate jet power. We searched in the available literature for all the published optical spectra and multiwavelength studies useful to characterize the sources. We found new, missed, or even forgotten information leading to a substantial change in the redshift values and classification of many sources. We present here the preliminary results of this analysis and some statistics based on the gamma-ray sources with right ascension within the interval 0 h --12 h (J2000). Although flat-spectrum radio quasars and BL Lac objects are still the dominant populations, there is a significant increase in the number of other objects, such as misaligned AGN, narrow-line Seyfert 1 galaxies, and Seyfert galaxies. We also introduced two new classes of objects: changing-look AGN and ambiguous sources. About one third of the sources remain unclassified.

79 ASTRONOMY AND ASTROPHYSICS↗

Dust Reverberation of 3C 273: Torus Structure and Lag–Luminosity Relation

We monitored the z = 0.158 quasar 3C 273 between 2015 and 2019 in the optical (BVrz) and near-infrared (JHK) with the aim to perform dust reverberation mapping. Accounting for host galaxy and accretion disk contributions, we obtained pure dust light curves in JHK. Cross correlations between the V band and the dust light curves yield an average rest-frame delay for the hot dust of τ {sub cent} ∼ 410 days. This is a factor of two shorter than that expected from the the dust ring radius R {sub x} ∼ 900 lt-day reported from interferometric studies. The dust covering factor (CF) is about 8%, much smaller than that predicted from the half covering angle of 45° found for active galactic nuclei (AGNs). We analyze the asymmetric shape of the correlation functions and explore whether an inclined biconical bowl-shaped dust torus geometry could bring these findings (τ {sub cent}, R {sub x} and CF) into a consistent picture. The hot varying dust emission originates from the edge of the bowl rim with a small covering angle 40° < θ < 45°, and we see only the near side of the biconus. Such a dust gloriole with R {sub x} = 900 ± 200 lt-day and an inclination 12° matches the data remarkably well. Comparing the results of 3C 273 with literature for less luminous AGN, we find a lag–luminosity relation τ ∝ L {sup α} with α = 0.33–0.40, flatter than the widely adopted relation with α ∼ 0.5. We address several explanations for the new lag–luminosity relation.

79 ASTRONOMY AND ASTROPHYSICS↗

ALMA Band-to-band Phase Referencing: Imaging Capabilities on Long Baselines and High Frequencies

High-frequency long-baseline experiments with the Atacama Large Millimeter/submillimeter Array were organized to test the high angular resolution imaging capabilities in the submillimeter wave regime using baselines up to 16 km. Four experiments were conducted, two Band 7 (289 GHz) and two Band 8 (405 GHz) observations. Phase correction using band-to-band (B2B) phase referencing was used with a phase calibrator only [[Formula]] away observed in Band 3 (96 GHz) and Band 4 (135 GHz), respectively. In Band 8, we achieved the highest resolution of 14 × 11 mas. We compared the synthesis images of the target quasar using 20 and 60 s switching cycle times in the phase referencing. In Band 7, the atmosphere had good stability in phase rms (<0.5 rad over 2 minutes), and there was little difference in image coherence between the 20 and 60 s switching cycle times. One Band 8 experiment was conducted under a worse phase rms condition (>1 rad over 2 minutes), which led to a significantly reduced coherence when using the 60 s switching cycle time. One of our four experiments indicates that the residual phase rms error after phase referencing can be reduced to 0.16 rad at 289 GHz in using the 20 s switching cycle time. Such conditions would meet the phase correction requirement of image coherence of >70% in Band 10, assuming a similar phase calibrator separation angle, emphasizing the need for such B2B phase referencing observing at high frequencies.

79 ASTRONOMY AND ASTROPHYSICS↗

On the Accuracy of the ALMA Flux Calibration in the Time Domain and across Spectral Windows

A diverse array of science goals requires accurate flux calibration of observations with the Atacama Large Millimeter/submillimeter array (ALMA); however, this goal remains challenging due to the stochastic time-variability of the “grid” quasars ALMA uses for calibration. In this work, we use 343.5 GHz (Band 7) ALMA Atacama Compact Array observations of four bright and stable young stellar objects over seven epochs to independently assess the accuracy of the ALMA flux calibration and to refine the relative calibration across epochs. The use of these four extra calibrators allows us to achieve an unprecedented relative ALMA calibration accuracy of ~3%. On the other hand, when the observatory calibrator catalog is not up to date, the Band 7 data calibrated by the ALMA pipeline may have a flux calibration poorer than the nominal 10%, which can be exacerbated by weather-related phase decorrelation when self-calibration of the science target is either not possible or not attempted. We also uncover a relative flux calibration uncertainty between spectral windows of 0.8%, implying that measuring spectral indices within a single ALMA band is likely highly uncertain. We thus recommend various methods for science goals requiring high flux accuracy and robust calibration, in particular, the observation of additional calibrators combined with a relative calibration strategy, and observation of solar system objects for high absolute accuracy.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling Stochastic Variability in Multiband Time-series Data

In preparation for the era of time-domain astronomy with upcoming large-scale surveys, we propose a state-space representation of a multivariate damped random walk process as a tool to analyze irregularly-spaced multifilter light curves with heteroscedastic measurement errors. We adopt a computationally efficient and scalable Kalman filtering approach to evaluate the likelihood function, leading to maximum O(k 3 n) complexity, where k is the number of available bands and n is the number of unique observation times across the k bands. This is a significant computational advantage over a commonly used univariate Gaussian process that can stack up all multiband light curves in one vector with maximum O(k 3 n 3 ) complexity. Using such efficient likelihood computation, we provide both maximum likelihood estimates and Bayesian posterior samples of the model parameters. Three numerical illustrations are presented: (i) analyzing simulated five-band light curves for a comparison with independent single-band fits; (ii) analyzing five-band light curves of a quasar obtained from the Sloan Digital Sky Survey Stripe 82 to estimate short-term variability and timescale; (iii) analyzing gravitationally lensed g- and r-band light curves of Q0957+561 to infer the time delay. Two R packages, Rdrw and timedelay, are publicly available to fit the proposed models.

79 ASTRONOMY AND ASTROPHYSICS↗

Alert Classification for the ALeRCE Broker System: The Light Curve Classifier

We present the first version of the Automatic Learning for the Rapid Classification of Events (ALeRCE) broker light curve classifier. ALeRCE is currently processing the Zwicky Transient Facility (ZTF) alert stream, in preparation for the Vera C. Rubin Observatory. The ALeRCE light curve classifier uses variability features computed from the ZTF alert stream and colors obtained from AllWISE and ZTF photometry. We apply a balanced random forest algorithm with a two-level scheme where the top level classifies each source as periodic, stochastic, or transient, and the bottom level further resolves each of these hierarchical classes among 15 total classes. This classifier corresponds to the first attempt to classify multiple classes of stochastic variables (including core- and host-dominated active galactic nuclei, blazars, young stellar objects, and cataclysmic variables) in addition to different classes of periodic and transient sources, using real data. We created a labeled set using various public catalogs (such as the Catalina Surveys and Gaia DR2 variable stars catalogs, and the Million Quasars catalog), and we classify all objects with ≥6 g-band or ≥6 r-band detections in ZTF (868,371 sources as of 2020 June 9), providing updated classifications for sources with new alerts every day. For the top level we obtain macro-averaged precision and recall scores of 0.96 and 0.99, respectively, and for the bottom level we obtain macro-averaged precision and recall scores of 0.57 and 0.76, respectively. Updated classifications from the light curve classifier can be found at the ALeRCE Explorer website (http://alerce.online).

47 OTHER INSTRUMENTATION↗

Functional Principal Component Analysis of Radio–Optical Reference Frame Tie

The Gaia optical reference frame is intrinsically undefined with respect to global orientation and spin, so it needs to be anchored in the radio-based International Celestial Reference Frame (ICRF) to provide a referenced and quasi-inertial celestial coordinate system. The link between the two fundamental frames is realized through two samples of distant extragalactic sources, mostly active galactic nuclei and quasars, but only the smaller sample of radio-loud ICRF sources with optical counterparts is available to determine the mutual orientation. The robustness of this link can be mathematically formulated in the framework of functional principal component analysis using a set of vector spherical harmonics to represent the differences in celestial positions of the common objects. The weakest eigenvectors are computed, which describe the greatest deficiency of the link. The deficient or poorly determined terms are specific vector fields on the sphere that carry the largest errors of absolute astrometry using Gaia in reference to the ICRF. This analysis provides guidelines for the future development of the ICRF maximizing the accuracy of the link over the entire celestial sphere. A measure of the robustness of a least-squares solution, which can be applied to any linear model fitting problem, is introduced to help discriminate between reference frame tie models of different degrees.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Detecting Long-period Variability in the SDSS Stripe 82 Standards Catalog

We report the results of a search for long-period (100 < P < 600 days) periodic variability in the SDSS Stripe 82 standards catalog. The SDSS coverage of Stripe 82 enables such a search because there are on average 20 observations per band in ugriz bands for about one million sources, collected over about 6 yr, with a faint limit of r ~ 22 mag and precisely calibrated 1%–2% photometry. We calculated the periods of variable source candidates in this sample using the Lomb–Scargle periodogram and considered the three highest periodogram peaks in each of the gri filters as relevant. Only those sources with gri periods consistent within 0.1% were later studied. We use the Kuiper statistic to ensure uniform distribution of data points in phased light curves. We present five sources with the spectra consistent with quasar spectra and plausible periodic variability. This SDSS-based search bodes well for future sensitive large-area surveys, such as the Rubin Observatory Legacy Survey of Space and Time, which, due to its larger sky coverage (about a factor of 60) and improved sensitivity (~2 mag), will be more powerful for finding such sources.

79 ASTRONOMY AND ASTROPHYSICS↗