5

17

## Changing-look Active Galactic Nuclei from the Dark Energy Spectroscopic Instrument. II. Statistical Properties from the First Data Release

WEI-JIAN GUO , 1 HU ZOU , 1 CLAIRE L. GREENWELL , 2 DAVID M. ALEXANDER , 3, 2 VICTORIA A. FAWCETT , 4 ZHIWEI PAN , 5 MAŁGORZATA SIUDEK , 6 JESSICA NICOLE AGUILAR, 7 STEVEN AHLEN , 8 DAVID BROOKS, 9 TODD CLAYBAUGH, 7 KYLE DAWSON, 10 AXEL DE LA MACORRA , 11 PETER DOEL, 9 ANDREU FONT-RIBERA , 9, 12 ENRIQUE GAZTA˜ NAGA, 13,14,6 SATYA GONTCHO A GONTCHO , 7 GASTON GUTIERREZ, 15 ROBERT KEHOE, 16 THEODORE KISNER , 7 MARTIN LANDRIAU , 7 LAURENT LE GUILLOU , 17 MARC MANERA , 18, 12 AARON MEISNER , 19 RAMON MIQUEL, 20,12 JOHN MOUSTAKAS , 21 FRANCISCO PRADA , 22 GRAZIANO ROSSI, 23 EUSEBIO SANCHEZ , 24 MICHAEL SCHUBNELL, 25,26 DAVID SPRAYBERRY, 19 JIPENG SUI, 27 GREGORY TARL´ E , 26 BENJAMIN ALAN WEAVER, 19 YUN-AO XIAO , 27 AND SIWEI ZOU 28, 29

1 Key Laboratory of Optical Astronomy, National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100012, China

Email:guowj@bao.ac.cn, zouhu@nao.cas.cn

2 Institute for Computational Cosmology, Department of Physics, Durham University, South Road, Durham DH1 3LE, UK

3 Centre for Extragalactic Astronomy, Department of Physics, Durham University, South Road, Durham, DH1 3LE, UK

4 School of Mathematics, Statistics and Physics, Newcastle University, Newcastle, UK

Kavli Institute for Astronomy and Astrophysics at Peking University, PKU, 5 Yiheyuan Road, Haidian District, Beijing 100871, P.R. China

6 Institute of Space Sciences, ICE-CSIC, Campus UAB, Carrer de Can Magrans s/n, 08913 BELlaterra, Barcelona, Spain

7 Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA

8 Physics Dept., Boston University, 590 Commonwealth Avenue, Boston, MA 02215, USA

9 Department of Physics &amp; Astronomy, University College London, Gower Street, London, WC1E 6BT, UK

10 Department of Physics and Astronomy, The University of Utah, 115 South 1400 East, Salt Lake City, UT 84112, USA

11 Instituto de F´ ısica, Universidad Nacional Aut´noma de M´ exico, Cd. de M´ exico C.P . 04510, M´ exico o

12 Institut de F´ ısica d'Altes Energies (IFAE), The Barcelona Institute of Science and Technology, Campus UAB, 08193 BELlaterra Barcelona, Spain 13 Institut d'Estudis Espacials de Catalunya (IEEC), 08034 Barcelona, Spain

14

Institute of Cosmology and Gravitation, University of Portsmouth, Dennis Sciama Building, Portsmouth, PO1 3FX, UK

15 Fermi National Accelerator Laboratory, PO Box 500, Batavia, IL 60510, USA

16 Department of Physics, Southern Methodist University, 3215 Daniel Avenue, Dallas, TX 75275, USA

Sorbonne Universit´ e, CNRS/IN2P3, Laboratoire de Physique Nucl´ eaire et de Hautes Energies (LPNHE), FR-75005 Paris, France

18 Departament de F´ ısica, Serra H´nter, Universitat Aut`noma de Barcelona, 08193 BELlaterra (Barcelona), Spain u o

19 NSF NOIRLab, 950 N. Cherry Ave., Tucson, AZ 85719, USA

20 Instituci´ Catalana de Recerca i Estudis Avanc ¸ats, Passeig de Llu´ ıs Companys, 23, 08010 Barcelona, Spain o

21 Department of Physics and Astronomy, Siena College, 515 Loudon Road, Loudonville, NY 12211, USA

22 Instituto de Astrof´ ısica de Andaluc´ ıa (CSIC), Glorieta de la Astronom´ ıa, s/n, E-18008 Granada, Spain

23 Department of Physics and Astronomy, Sejong University, Seoul, 143-747, Korea

24 CIEMAT, Avenida Complutense 40, E-28040 Madrid, Spain

25 Department of Physics, University of Michigan, Ann Arbor, MI 48109, USA

26 University of Michigan, Ann Arbor, MI 48109, USA

- 27 National Astronomical Observatories, Chinese Academy of Sciences, A20 Datun Rd., Chaoyang District, Beijing, 100012, P.R. China
- 28 Chinese Academy of Sciences South America Center for Astronomy, National Astronomical Observatories, CAS, Beijing 100101, China 29 Department of Astronomy, Tsinghua University, Beijing 100084, China

(Received xxx; Revised xxx; Accepted xxx)

Submitted to the Astrophysical Journal Supplement

## ABSTRACT

We present the identification of changing-look active galactic nuclei (CL-AGNs) from the Dark Energy Spectroscopic Instrument First Data Release and Sloan Digital Sky Survey Data Release 16 at z ≤ 0 9 . . To confirm the CL-AGNs, we utilize spectral flux calibration assessment via an [O III]-based calibration, pseudophotometry examination, and visual inspection. This rigorous selection process allows us to compile a statistical catalog of 561 CL-AGNs, encompassing 527 H β , 149 H α , and 129 Mg II CL behaviors. In this sample, we find

1) a 283:278 ratio of turn-on to turn-off CL-AGNs. 2) the critical value for CL events is confirmed around Eddington ratio λ Edd ∼ 0 01 . . 3) a strong correlation between the change in the luminosity of the broad emission lines (BEL) and variation in the continuum luminosity, with Mg II and H β displaying similar responses during CL phases. 4) the Baldwin-Phillips-Terlevich diagram for CL-AGNs shows no statistically difference from the general AGN catalog. 5) five CL-AGNs are associated with asymmetrical mid-infrared flares, possibly linked to tidal disruption events. Given the large CL-AGNs and the stochastic sampling of spectra, we propose that some CL events are inherently due to typical AGN variability during low accretion rates, particularly for CL events of the singular BEL. Finally, we introduce a Peculiar CL phase, characterized by a gradual decline over decades in the light curve and the complete disappearance of entire BEL in faint spectra, indicative of a real transition in the accretion disk.

Keywords: Accretion (14); Active galaxies (17); Active galactic nuclei (16); Supermassive black holes (1663); Catalogs (205);

## 1. INTRODUCTION

Supermassive black holes (SMBHs) are easily identified when they accrete gas and emit intense radiation (Blandford &amp; Znajek 1977; Alexander &amp; Hickox 2012; Kormendy &amp; Ho 2013). Active galactic nuclei (AGNs) or quasars, which are the extreme cases of accretion, exhibit distinctive fluctuations across various electromagnetic wavelengths, spanning from gamma rays to radio (Rees 1984; Vaughan et al. 2003; Ho 2008; Abdo et al. 2010; Padovani et al. 2017). In the optical and ultraviolet spectrum, the variability of the AGN continuum typically ranges between 10-30%, possibly due to the thermal instability of the accretion disk that primarily generates the continuum radiation (Giveon et al. 1999; Kelly et al. 2009; Zuo et al. 2012). An interesting behavior of AGN variability is the changing-look (CL) phenomenon, whose deeper understanding may require additional physical mechanisms (Shappee et al. 2014; Denney et al. 2014; Ricci &amp; Trakhtenbrot 2023). The variability of CL-AGN are distinct from that of typical AGNs, often described by the damped random walk (DRW) model, which characterize the variability as white noise (Kelly et al. 2009; Kasliwal et al. 2015).

CL-AGN offers a distinct avenue for enhancing our comprehension of AGN populations as dramatic changes in broad emission lines (BEL) carry out the potential link among different types of AGN in their spectra. Specifically, CL-AGNs exhibit BEL either appearing (turn-on) or disappearing (turn-off) over periods spanning several months to years (MacLeod et al. 2016; Yang et al. 2018; Graham et al. 2020; Wolf et al. 2020; Green et al. 2022; Hon et al. 2022; Guo et al. 2024), emanating from the ionized gas within the Broad Line Region (BLR). These variations defy the conventional unified model that classifies AGN types according to the observer's line of sight, where the dusty torus may obscure the presence of BEL due to a high dust covering factor in Type 2 (narrow-line) AGNs (Antonucci 1993; Netzer 2015). An alternative model proposing an evolution sequence among different AGNtypes (Type 1 → 1.2/1.5 → 1.8/1.9 → 2) suggests that the accretion rate may influence this progression, implying that Type 2 AGNs could represent a transient phase akin to 'flameout' Type 1 AGNs where the continuum is temporarily off (Penston &amp; Perez 1984; Trump et al. 2011; Stern &amp; Laor 2012; Elitzur et al. 2014; Runnoe et al. 2016).

Furthermore, during CL events, substantial fluctuations in the optical or mid-infrared bands, often exceeding one magnitude, are also observed (Sheng et al. 2017; MacLeod et al. 2019). These significant variations pose a challenge to the standard accretion disk model, which assumes stable conditions and and typical transition times of several hundred years for viscous timescales (Shakura &amp; Sunyaev 1973; Shapovalova et al. 2010; Runnoe et al. 2016; Gezari et al. 2017; Noda &amp; Done 2018). CL-AGNs also provide an excellent opportunity to explore the relationship between AGNs and their host galaxies, such as co-evolution and feedback, as the spectra exhibit apparent starlight features in the dim state (Gezari et al. 2017; Frederick et al. 2019). The SMBH mass can be determined through the R BLR -L 5100 scaling relation by measuring BEL and continuum in the bright state or M BH -σ ⋆ scaling relation by measuring stellar velocity dispersion in the dim state (Yu et al. 2020; Jin et al. 2022). The population of CL-AGNs would allow for a detailed comparison and correction of the two scaling relations, especially concerning the population of low-mass SMBHs and young host galaxies.

Although the primary cause of the significant variations in both BEL and continuum is generally ascribed to the intrinsic accretion rate reaching a critical threshold (MacLeod et al. 2016, 2019; Green et al. 2022), certain puzzling behaviors hidden in these visual-defined CL-AGNs may have alternative originations. For example, tidal disruption events (TDEs) could contribute to this phenomenon, but their signal might be washed out by the AGN (Chan et al. 2020; Li et al. 2022). Additionally, variations in the dusty covering factor or the obscuration of dusty gas, in conjunction with changes in the accretion rate, could contribute to CL events (Runnoe et al. 2016; Gezari et al. 2017; Zeltyn et al. 2022). However, this phenomena appears to be comparatively

rare amongst typical CLAGN samples. Alternative theoretical frameworks, such as tidal torque generated by a close binary black hole (Wang &amp; Bon 2020), also hold promise in elucidating certain peculiar CL behaviors (Cao et al. 2023; Wu &amp; Gu 2023). To determine the relative contributions of these various processes to the CLAGN phenomenon requires careful and detailed analyses using large statistical samples constructed from a comprehensive spectroscopic dataset.

To date most CL-AGN samples are identified based on their highly unusual variability in optical or mid-IR light curves (Sheng et al. 2017; L´pez-Navas et al. 2022; Wang et al. 2023, 2024b), particularly those associated with turn-on CL-AGN. However, this o method excludes CL-AGNs whose spectra are host-dominant and photometric variability is diluted by the host galaxy. Therefore, a detailed examination of the CL phenomenon, including the dependence on type transition and occurrence rate, still necessitates large-area spectroscopic projects. To enhance our comprehension of the AGN population and the CL physical mechanism, we aim to compile a statistically robust catalog of CL-AGNs following the method outlined Guo et al. (2024) through cross-matching Dark Energy Spectroscopic Instrument (DESI) First Data Release (DR1) and Sloan Digital Sky Survey (SDSS) Data Release 16 (DR16) database.

The paper is structured as follows. Section 2 details the spectral and photometric data. Section 3 describes the selection process for CL-AGNs samples, encompassing spectral fitting and flux calibration assessment. Section 4 presents the measurement of emission lines and physical properties of CL-AGNs. Discussions and summary are found in Section 5 and 6, respectively. Throughout the study, we adopt a Λ CDM cosmology with H 0 = 67 km s -1 Mpc -1 , Ω Λ = 0 68 . , and Ω m = 0 32 . as reported by Planck Collaboration et al. 2020, which is consistent within 2 σ with first cosmological results from DESI (DESI Collaboration et al. 2024a).

## 2. DATA

2.1. Spectroscopic Survey

2.1.1. DESI

DESI, a Stage IV ground-based dark energy experiment, focuses primarily on investigating the nature of dark-energy experiment and exploring cosmological constraints through measurements of baryon acoustic oscillations (BAO) (Levi et al. 2013; DESI Collaboration et al. 2016a,b; Abareshi et al. 2022; DESI Collaboration et al. 2022). The key projects of DESI encompass two-point clustering measurements and validation (DESI Collaboration et al. 2024b,c), BAO measurements from galaxies, quasars, and the Lyman-alpha forest (DESI Collaboration et al. 2024d,e), a comprehensive study of the shapes of galaxies and quasars (DESI Collaboration et al. 2024f), as well as cosmological constraints derived from BAO measurements (DESI Collaboration et al. 2024a), full-shape measurements (DESI Collaboration et al. 2024g), and constraints on primordial non-Gaussianities (DESI Collaboration et al. 2024h).

To achieve these goals, the DESI team utilizes the NOIRLab 4m Mayall telescope located at Kitt Peak and equipped with 5,000 fibers on the focal plate to conduct a comprehensive and extensive optical multi-object spectral survey (Silber et al. 2023; Miller et al. 2023). The first phase of the DESI survey (DESI-I) reaches a depth of 23 magnitudes in the r -band and is projected to amass more than 40 million spectra of galaxies and quasars within the initial five years (DESI Collaboration et al. 2016a). The fibers feed ten three-arm spectrographs (1.5 ′′ optical diameter) and the spectral resolution for the three channels (blue: 3600-5900 ˚ , A green: 5660-7220 ˚ , A and red: 7470-9800 ˚ ) A is approximately λ/ ∆ λ ∼ 2100, λ/ ∆ λ ∼ 3200, and λ/ ∆ λ ∼ 4100, respectively (DESI Collaboration et al. 2016b).

To validate the quality of the spectroscopic data and the data processing/classification pipeline (Guy et al. 2022; Moustakas et al. 2023; Schlafly et al. 2023), the DESI team initiated a Survey Validation (SV) project, which primarily focused on approximately 1% of the Main Survey (DESI Collaboration et al. 2023a,b). Through visual inspection (VI), the data team also validated various preliminary target selections, including bright galaxy, luminous red galaxies, emission line galaxies and quasars (Lan et al. 2023; Raichoor et al. 2023; Zhou et al. 2022; Alexander et al. 2023; Hahn et al. 2022; Chaussidon et al. 2022; Juneau et al. 2024), aimed at improving the performance of the standard spectroscopic classifier (Guy et al. 2022; Bailey et al. 2023).

Upon completion of operations in Year 1 (from May 2021 to June 2022), DESI internally designates the Year 1 data as DR1. The spectral data in the DR1 catalog dataset stems from a combination of various tiles or exposures, some of which involve repeated observations of specific targets. In this work, we adopt the mean value as the spectral Modified Julian Date (MJD) within the DR1 catalog, encompassing over one million quasar and ten million galaxy spectra directly classified by the 'Redrock' spectral template-fitting code, which is a software package used for redshift determination (Brodzeller et al. 2023). The final quasar catalog for DR1 will incorporate additional insights derived from the 'Afterburner' (a broad MgII emission line identifier) and 'QuasarNET' classification methodologies (Farr et al. 2020). Here we only adopt Redrock classification results as our parent sample, but the future quasar catalog would further supplement the quantity of the CL-AGNs. The spectra are corrected for the galactic extinction of our Milky Way by using the extinction curve of Fitzpatrick (1999) with R V = 3 1 . .

SDSS utilized the 2.5 m Sloan Foundation Telescope at Apache Point Observatory to study the large-scale distribution of galaxies and quasars and establish an imaging and spectroscopic legacy for the astronomical community. Observations have been conducted routinely since April 2000 (York et al. 2000; Gunn et al. 2006). SDSS-I and SDSS-II performed a spectroscopic survey using a 640-fiber-fed pair of multiobject double spectrographs with a 3 ′′ optical diameter, covering wavelengths from 3800 ˚ Ato 9200 ˚ Aand achieving a median resolution of approximately λ/ ∆ λ ∼ 2000 (Abazajian et al. 2004; Adelman-McCarthy et al. 2008). For the investigation of dark energy and cosmological parameters through the Baryon Oscillation Spectroscopic Survey (BOSS), which is the primary survey of SDSS-III, SDSS-IV, and SDSS-V, upgraded to 1000-fiber multiobject double spectrographs with a 2 ′′ optical diameter. The resolution improved from λ/ ∆ λ ∼ 1300 at 3600 ˚ A to λ/ ∆ λ ∼ 3000 at 10600 ˚ A at 10600 ˚ A after 2008 (Abazajian et al. 2009; Ahn et al. 2012; Eisenstein et al. 2011; Smee et al. 2013; Alam et al. 2017).

Over two decades, the SDSS has amassed a significant spectral database encompassing millions of quasars and galaxies. This collection is invaluable for investigating the variability of AGNs, particularly in the context of CL-AGNs. Leveraging the Data Release 16 Quasar (DR16Q) catalog for the most recent and comprehensive quasar dataset (Lyke et al. 2020), we utilize DR16 and DR16Q to extract the galaxy and quasar samples respectively (Lyke et al. 2020; Ahumada et al. 2020). The catalogs contain a compilation of over 0.7 million quasar spectra and 4 million galaxy spectra. Within the time domain spectroscopic survey (TDSS) and reverberation mapping (RM) program, DR16 includes repeated observations to capture multiple spectra of stars and quasars exhibiting unusual characteristics (Ahumada et al. 2020). We filter out these duplicate observations between DR16 and DR16Q, retaining only the initial exposure in the quasar catalog.

## 2.2. Image Survey and Light Curve

To delve deeper into the galaxy morphology of CL-AGNs, false-color images, combined from different wavelengths into a single image, of CL-AGNs from the DESI and SDSS datasets were obtained through the DESI Legacy Survey 1 . The DESI Legacy Survey Data Release 10 (LS-DR10) includes green ( g ) at 4671 ˚ , red ( A r ) at 6624 ˚ , near infrared ( ) at 8060 ˚ , infrared A i A ( z ) at 9200 ˚ for a combined ( A griz ) image from the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z -band Legacy Survey, covering a period from 2014 to 2019 with a total sky coverage of 14,000 deg 2 (Flaugher et al. 2015; Dey et al. 2019; Zhou et al. 2022; Schlegel et al. 2023). Photometric magnitudes in the DESI legacy surveys are derived from single-band coadded images, with the MJD being the adopted average value. Between 1998 and 2008, the SDSS legacy survey compiled a 14,500 deg 2 ugriz image (ultraviolet at 3543 ˚ ) using the 2.5 m Sloan Telescope, primarily from SDSS-I and A SDSS-II (Fukugita et al. 1996; Gunn et al. 1998). The images in the five bands were obtained through drift-scan observations with the same MJD.

Given that CL-AGN variations typically occur over more than a decade, it is essential to have long-term photometric data to investigate their characteristic patterns. To achieve this, we have utilized light curves from various sources: Catalina Real-time Transient Survey (CRTS; V -band 2 ; Drake et al. 2009), Pan-STARRS1 (PS1; g - and r -band; Chambers et al. 2016), the Palomar Transient Factory (PTF; g -band; Law et al. 2009), and the Zwicky Transient Facility (ZTF; g - and r -band; Masci et al. 2019) for the optical band, as well as mid-infrared light curves from the Wide-field Infrared Survey Explorer (WISE; W1- and W2-band; Mainzer et al. 2014) and the Near-Earth Object Wide-field Infrared Survey Explorer (NEOWISE; W1- and W2-band).

CRTS is an astronomical survey designed to detect transients using three dedicated telescopes covering approximately 33,000 deg 2 of the sky (Drake et al. 2009). To enhance the depth and sky coverage of the survey, PS1 employs a 1.8-meter telescope to conduct a five-band, time-domain survey, encompassing about 3/4 of the entire sky (Chambers et al. 2016). PTF and ZTF are transient detection systems based on the Palomar Samuel Oschin 48-inch telescope, providing sky coverage of approximately 25,000 to 30,000 deg 2 respectively (Law et al. 2009; Masci et al. 2019). In 2010, WISE operated the all-sky survey with four bands at 3.4, 4.6, 12, and 22 µ m (W1-, W2-, W3-, and W4-band). Since 2014, NEOWISE has started the routine survey (six months) in the W1- and W2-band, which primarily traces the dusty torus in AGN.

## 3. SAMPLE SELECTION

To compile a catalog of CL-AGNs for Mg II, H β , and H α , we are limiting the redshift to z ≤ 0 9 . , which ensures a more reliable calculation of black hole mass through the R BLR -L relationship using H β and H α rather than Mg II and C IV. Additionally, the spectral flux calibration for high redshift quasars poses a challenge. However, we can use [O III]-based calibration and pseudophotometry to address this issue for objects at z ≤ 0 9 . . Figure 1 - 3 showcase representative examples of CL-AGNs featuring

1 https://www.legacysurvey.org

2 The CRTS observations use an unfiltered system, in which the raw data are transformed and calibrated to the V -band based on a few 10's to 100's stars within each frame.

Figure 1. The light curves, images, and spectra for an Mg II CL-AGN. The top-left panel displays the optical light curve over 20 years, with data from SDSS (squares), CRTS (stars), PS1 (triangles), PTF (pentagons), DESI (diamonds), and ZTF (dots) where available. Spectral pseudophotometry is marked with red 'x' markers for DESI and blue 'x' markers for SDSS. The middle-left panel shows the W1 and W2 bands of the MIR light curves from WISE (squares) and NEOWISE (diamonds). The top-right and middle-right panels display false-color images from the DESI Legacy Survey for SDSS and DESI, respectively, with circles indicating the fiber diameter for SDSS and DESI. The bottom-left panel shows smoothed spectra from SDSS and DESI, represented by blue and red lines. The bottom-right panel features the differential spectrum, obtained by subtracting the dim spectrum from the bright one, shown as a grey line.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00402v1-figure-618cb1b7-58a8-424c-82c9-2c745d17ed64.jpg)

Notes: 1) DESI photometry uses the average MJD from the different telescopes from 2013 to 2019. 2) the optical light curves are in the AB magnitude, while the mid-infrared light curves are in the Vega magnitude.

(The complete figure set, consisting of 561 images, is available online.)

MgII, H β , and H α , respectively. We selected our CL-AGN sample following Guo et al. (2024). We summarise the key steps in Figure 4: step (1) spectral integration; step (2) spectral fitting; and step (3) flux calibration assessment, which we describe in the following sub sections.

## 3.1. Parent Sample

We construct the parent sample from DESI DR1 and SDSS DR16/DR16Q catalogs. There are 44,728 AGN-AGN pairs or AGN-galaxy pairs by cross-matching SDSS DR16Q and DESI DR1 ( spectype = galaxy or spectype = QSO ). To obtain galaxy-AGN pairs, we also cross-match SDSS DR16 ( class = galaxy ) with DESI DR1 catalog ( spectype == QSO ). We keep 4,905 non-repeated galaxy-AGN pairs since the results contain duplicate sources between DR16 and DR16Q. Finally, the parent sample before the selection consists of 49,633 pairs. We caution that the final catalog might miss some CL-AGNs because we do not take into account galaxy-galaxy pairs from SDSS DR16 ( class = galaxy ) and DESI DR1 ( spectype = galaxy ), where an AGN might be misclassified as a galaxy.

Figure 2. The light curves, images, and spectra for an H β CL-AGN, J120447.91+170256.8, reported by Wang et al. (2019). The legend is the same as Figure 1.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00402v1-figure-f7569ece-3cae-4f5d-ba44-11e3d09c6dfb.jpg)

Building on previous studies (MacLeod et al. 2016; Guo et al. 2024), we utilize N σ to denote the significance of the variation in the BEL maximum flux and R to represent the total BEL flux change.

$$N _ { \sigma } = ( f _ { \text{bright} } - f _ { \text{dim} } ) / \sqrt { \sigma _ { \text{bright} } ^ { 2 } + \sigma _ { \text{dim} } ^ { 2 } },$$

where f and σ are the spectral flux and variance for a single BEL respectively in erg cm -2 s -1 ˚ A -1 .

$$R = ( F _ { \text{bright} } - F _ { \text{dim} } ) / F _ { \text{dim} },$$

where F bright or F dim is the BEL total flux in the bright or dim state.

The step (1) involves using the spectral integration method to calculate the total flux of the broad emission lines (BELs). We apply two selection criteria: Max(N ) σ &gt; 3 and R &gt; 1 5 . for at least one BEL (Mg,II, H β , or H α ) to perform a quick selection.

## 3.2. Spectral Decomposition

We perform spectral decomposition with a reduced sample of AGNs to more accurately quantify the BEL flux change and produce a final CL-AGN catalog. In addition to determining the luminosity change of Mg II, H β and H α BEL through spectral decomposition, we can calculate the continuum luminosity and BEL to derive further the AGN's physical properties, such as black hole mass and Eddington ratio. Before fitting, we rebin the spectral flux and adjust it to 4 A per pixel to enhance the ˚ signal-to-noise ratio (SNR) per pixel, sacrificing resolution for both SDSS and DESI data.

Figure 3. The light curves, images, and spectra for an H α CL-AGN, J143016.06+230344.5, regarded as SMBH binary by Jiang et al. (2022). The legend is the same as Figure 1.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00402v1-figure-fa35b84f-8320-4d8f-8fd3-987dd3341797.jpg)

We undertake the spectral decomposition using DASpec , which utilizes the Levenberg-Marquardt method to minimize χ 2 and fit all spectral components simultaneously (Du et al. 2018). Consistent with Hu et al. (2015) and Guo et al. (2022), the spectral components include:

- · a set of templates with 640 Myr, 900 Myr, 1.4 Gyr, 2.5 Gyr, 5 Gyr, 11 Gyr stellar-population ages from Bruzual &amp; Charlot (2003) to reproduce the host galaxy starlight;
- · a featureless power law for the quasar continuum;
- · a Fe II pseudo-continuum template from Boroson &amp; Green (1992);
- · a double Gaussian for the broad component of Mg II, H δ , H γ , He II, H β , and H α ;
- · a double Gaussian 3 for the narrow component H β , [O III] λλ 4959, 5007;
- · a single Gaussian for the rest of the forbidden lines.

In DASpec , we fit the entire spectrum to obtain an initial value of the AGN continuum and host galaxy. The different emission line components are added to the decomposition, which would simultaneously be fitted in the results. The Balmer continuum is neglected by skipping the corresponding wavelength since the Balmer jump is insignificant in most CL-AGNs. We choose four potential emission line fitting windows depending on the spectral coverage:

- · 2600A ˚ -3000A ˚ for Mg II λλ 2796, 2803;
- · 4500A ˚ -5200A ˚ for He II λ 4686, H β λ 4861, and [O III] λλ 4959, 5007;
- · 3700A ˚ -4500A ˚ for [O II] λ 3727, [Ne III] λ 3869, H δ λ 4102, and H γ λ 4340;
- · 6200A ˚ -6800A ˚ for [O I] λ 6300, [N II] λ 6548, H α λ 6563, [N II] λ 6583, and [Si II] λλ 6716, 6731.

3 When the [O III] is affected by outflow, a single Gaussian profile may not be sufficient to characterize its complex profile accurately.

Figure 4. The selection flow for the CL-AGNs. The top dash lines represent the cross-matching between DESI and SDSS. The orange diamonds are key target-selection steps, and the grey rectangles are the selected sample. The BEL flux in the first (second) orange diamond is calculated by the spectral integration (fitting).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00402v1-figure-54e33338-33be-48d0-a67e-d496f604465c.jpg)

We tie the profiles of narrow components in each emission line fitting window. We also adopt the minimize χ 2 of the host galaxy template as the best host galaxy in the fitting results. However, we note that the decomposition between the continuum and host, especially for H α CL-AGNs, is less reliable as the flux and power-law index of the continuum are challenging to ascertain when the host galaxy dominates the spectrum. Figure 5 depicts an example of spectral fitting results for J162829.18+432948.5, where the continuum flux in the dim state remains a challenge. Additionally, analyzing the host galaxies of CL-AGNs based on stellar population synthesis would involve considering dust attenuations and star formation history through additional tools, such as FastSpecFit 4 (Moustakas et al. 2023).

Figure 5. An example of our spectral decomposition for SDSS (left panel) and DESI (right panel). The top panel shows the decomposition results, and the bottom panel shows the residuals. The black line represents our best fit. The fitting window is shown in the grey shadow region. The orange and blue lines represent the power-law for the continuum and the host galaxy template, respectively. The broad components and narrow components are marked in different colored lines.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00402v1-figure-3fd06b36-fea8-4cf9-abfa-5386a968ab29.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00402v1-figure-6c162ad4-187c-454b-b9a0-5bef32646ce9.jpg)

Figure 6. The histograms for spectral variation of CL-AGNs and parent samples at different redshift bins. The left panel shows the number of turn-on (grey) and turn-off (pink) patterns for CL-AGNs. The right panel shows the number of turn-bright (grey) and turn-dim (pink) patterns for the parent sample. Note that we directly use the median value of the differential spectrum to determine the turn-bright (DESI-SDSS &gt; 0) or turn-dim (DESI-SDSS &lt; 0).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00402v1-figure-ecf737c5-cd65-47e7-9b82-35ffb5f6ce53.jpg)

Within the BEL flux from spectral fitting in step (2), we utilize R &gt; 1 5 . to select the remaining subsample because the BEL flux might change after subtracting an old-age host galaxy template with absorption features in H β . After re-selection, the catalog contains 924 pairs remaining.

## 3.3. Flux Calibration Assessment

A careful flux calibration assessment is conducted through a [O III]-based calibration, pseudo-photometry examination, and VI because differences in the spectral flux calibration between the SDSS and DESI can be ascribed to the deviations in fiber positioning or fiber drop (Guo et al. 2020).

We adopt a 20% flux difference threshold for [O III] (5007 ˚ ) between DESI and SDSS as a calibration step. A If the flux difference exceeds this threshold, we discard the object instead of scaling their spectra based on the [O III] flux. Typically, the

[O III] flux should be consistent across two epochs due to the narrow line region (NLR), hundreds or thousands of parsecs in size, which should not vary over decadal timescales. However, the difference in fiber diameters between SDSS (2 ′′ or 3 ′′ ) and DESI (1.5 ′′ ) becomes more significant at z ≤ 0 25 . , potentially causing a loss of [O III] flux due to variations in seeing and extended NLR size. Compared to Guo et al. (2024), we chose a 20% tolerance to account for numerous CL-AGNs even at z ≤ 0 1 . , and pseudo-photometry could further verify the spectral flux calibration.

We use pseudo-photometry examination to eliminate spurious CL-AGNs, which pass our selection criteria due to a spectral flux calibration issue. Pseudo-photometry magnitudes are conducted by applying specific photometric filters from CRTS, PS1, or ZTF to the spectral flux from SDSS or DESI. The light curves are combined within 10 days for optical and 180 days for mid-infrared band. Pseudo-photometry magnitudes are given in Figure 1, which is used to compare with the light curve.

VI is also employed to scrutinize CL-AGNs whose spectra fail to pass [O III]-based calibration and pseudo-photometry examination. Objects meeting the following criteria are retained: 1) weak [O III] or low SNR in the [O III] region for accurate fitting in dim spectra and passing the pseudo-photometry test; 2) pseudo-photometry magnitudes which might not be accurate at z ≤ 0 25 . and whose light curve is contaminated by host galaxy contributions. For example, the extended host galaxy would dominate the photometry in ZTF while DESI fiber diameter only covers the central region in Figure 3. After the flux calibration assessment, we create the final catalog comprising 561 CL-AGNs.

## 4. RESULTS

In this sample, 561 CL-AGNs exhibit 527 instances of H β CL behavior, 149 instances of H α CL behavior, and 129 instances of Mg II CL behavior. An example CL-AGN is shown in Figure 1 where a distinct blue core is observed, attributed to the strong AGNcontinuum in the spectrum, in the false-color image of SDSS indicates an AGN system is present in J110455.17+011856.6. However, after J110455.17+011856.6 turns off, the false-color image of DESI only shows an extended galaxy morphology, despite the different exposure times across the griz bands in the DESI image. These false-color changes between the SDSS and DESI images directly demonstrate the 'bluer when brighter' phenomenon under extreme variation, which could provide a new alternative method for identifying CL-AGN candidates.

The turn-on (283) and turn-off (278) patterns are nearly equal in our sample. We categorized CL-AGNs based on their behavior (turn-on vs. turn-off) and redshift bins in Figure 6 to explore whether there is a redshift evolution or dependencies. To illustrate the homogeneity of this sample, we adopt the median value of the differential spectrum as the indicator (turn-bright or turn-dim pattern 5 ) of the parent sample in Figure 6. The ratio bias between turn-bright and turn-dim increases as the redshift bin decreases. One of the reasons may be that many turn-bright AGNs are misclassified to turn-dim as the flux contribution of the host increases. At z = 0 4 . ∼ 0 9 . the ratio bias is approximately unity, suggesting that the parent sample is less biased at these higher redshifts. Whether there is a redshift dependency investigated at z &gt; 0 9 . will be the part of our future work on high-redshift sample. As for lower redshift bins, the reason why the number of turn-on CL-AGNs markedly surpasses turn-off CL-AGNs at z ≤ 0 2 . is unclear and complicated, which might be related to selection effects in our work or target selection bias in DESI or SDSS. For example, the classification pipeline of DESI or SDSS may classify weak AGN into galaxy at low redshift.

## 4.1. Physical Properties

This section provides the physical properties of CL-AGNs, including monochromatic luminosities, black hole masses, and Eddington ratios, which are summarized in Table 1. The monochromatic luminosity L 5100 at 5100 A is derived from spectral ˚ decomposition with a power-law model. Considering the degeneracy between continuum flux and power law index for the hostdominated AGN, the extrapolated monochromatic luminosity L 2500 at 2500 A might supply an additional measurement for the ˚ host-dominated systems.

We utilize the BEL from the bright phase spectra obtained from SDSS or DESI to estimate the black hole mass. Specifically, we measure the full width at half maximum (FWHM) of the H β line ( V H β ) or H α line ( V H α ):

$$M _ { \bullet } = f _ { \text{BLR} } \frac { V _ { \text{BLR} } ^ { 2 } R _ { \text{BLR} } } { G }, \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where G is the gravitational constant and f BLR = 1 12 . is the virial factor from Woo et al. (2015). The BLR radius R BLR is determined using the size-luminosity ( R H β -L 5100 ) scaling relation, which is a well-established method in RM studies:

$$\log ( R _ { \text{H} \beta } / \text{ld} ) = K _ { 1 } + \alpha _ { 1 } \log \ell _ { 4 4 },$$

5 For normal AGN, turn-bright refers to when the average flux of the DESI spectrum is greater than that of the SDSS spectrum, while turn-dim refers to when the average flux of the SDSS spectrum is greater than that of the DESI spectrum.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00402v1-figure-97e405ab-a9b1-445c-a6aa-5ed2db99b522.jpg)

Figure 7. The histograms for monochromatic luminosity at 5100 ˚ A(top-left), black hole mass (top-right), Eddington ratio from 5100 ˚ A(bottomleft) and 2500 ˚ A(bottom-right). The blue represents the bright state while the red represent the dim state of CL AGN. The grey histogram is the properties from SDSS DR14Q of z &lt; 0 9 . from Rakshit et al. (2020). The dashed line is the median value.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00402v1-figure-2c6ae400-f5ee-47fa-9022-79c95d643552.jpg)

•

•

Figure 8. The histograms black hole mass measurements from dim state vs. bright state of the H β (left) and H α (right) if available.

where ℓ 44 = L 5100 / 10 44 erg s -1 and the coefficient values are K 1 = 1 527 . ± 0 031 . and α 1 = 0 533 . ± 0 034 . are compiled by Bentz et al. (2013).

For the broad H α line, we apply a similar size-luminosity relation:

$$\log ( R _ { \text{H} \alpha } / \text{ld} ) = K _ { 2 } + \alpha _ { 2 } \log \ell _ { 4 4 },$$

with coefficients K 2 = 1 59 . ± 0 05 . and α 2 = 0 58 . ± 0 04 . are taken from Cho et al. (2023).

Figure 9. The variation of monochromatic luminosity at 5100 ˚ A(top) and 2500 ˚ A(bottom) vs. the variation of broad Mg II (left), H β (middle), and H α (right) luminosity. The red line is the linear fitting between continuum and different broad BEL.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00402v1-figure-2a4dc1c8-bbe8-4d55-b199-c33573185032.jpg)

With monochromatic luminosity and black hole mass, we can calculate the Eddington ratio, which is a key indicator of the accretion efficiency relative to the CL phenomenon. The Eddington ratio for both λ Edd , 5100 and λ Edd , 2500 is determined using:

$$\lambda _ { \text{Edd} } = L _ { \text{bol} } / L _ { \text{Edd} } = C _ { \text{bol} } \, L _ { \lambda } / L _ { \text{Edd} },$$

where C bol = 9 26 . for L 5100 and C bol = 5 00 . for L 2500 are the bolometric correction factors provided by Shen et al. (2008) and Heckman et al. (1997).

In Figure 7, we present the distributions of 5100 A luminosity ( ˚ L 5100 ), black hole mass ( M · ), and Eddington ratio ( λ Edd 5100 , and λ Edd 2500 , ) for the CL-AGNs in our sample. Generally, the median value of the 5100 A luminosity is ˚ log L 5100 = 43 67 erg s . -1 in the dim state and log L 5100 = 44 10 erg s . -1 in the bright state. For comparison, the median value for the SDSS DR14Q sample ( z ≤ 0 9 . ) is log L 5100 = 44 38 erg s . -1 , indicating a 0.7 dex disparity between CL-AGNs in the dim state and the overall SDSS-identified AGN population.

The black hole mass distribution for CL-AGNs spans 10 6 5 . to 10 9 5 . M ⊙ , showing no significant difference from the SDSS DR16Q sample in the top-right panel of Figure 7. The slight excess in black hole mass for CL-AGNs might be attributed to continuum overestimation for the host-dominated AGN, leading to a larger BLR size. We also compare black hole mass for some CL-AGNs with broad H α and H β lines available in both dim and bright states. In Figure 8, the black hole mass estimates from H β are more reliable than H α , as continuum measurements for H α CL-AGNs remain challenging due to host galaxy contamination.

In the bottom panel of Figure 7, we compare the Eddington ratio for CL-AGNs with the SDSS DR14Q sample ( log λ Edd = -1 06 . ) since Rakshit et al. (2020) gives the spectral decomposition for DR14Q. For CL-AGNs, the median values for our sample are log λ Edd 5100 , = -1 94 . and log λ Edd 2500 , = -2 15 . in the dim state, and log λ Edd 5100 , = -1 51 . and log λ Edd 2500 , = -1 50 . in the bright state. These findings align with previous studies, such as MacLeod et al. (2019), Green et al. (2022), and Wang et al. (2024b), which suggest that CL events typically occur around log λ Edd = -2 . Our results reinforce the view that this Eddington-ratio threshold is a critical driver of CL-AGN events.

4.2. BEL Variablity

We conducted regression analyses to investigate the correlation between the change in continuum luminosity and BEL luminosity for Mg II, H β , and H α . The results in Figure 9 demonstrate a strong correlation for all three BEL, which is consistent with previous studies by MacLeod et al. (2019) and Green et al. (2022), suggesting that the central continuum still photo-ionizes the BEL in CL-AGNs. These correlations have already been used in RM techniques for measuring black hole masses, highlighting the significance of the responsiveness of broad H β or H α .

In the study on Mg II variability, Sun et al. (2015) noted that the change in H β within the SDSS-RM sample exceeded that of MgII by a factor of 1.5. However, our analysis of CL-AGN samples did not reveal any significant discrepancy in the variations of H β and Mg II, which differs slightly from the findings of Green et al. (2022) and Zeltyn et al. (2024), whose samples show less variation in Mg II. Further quantitative investigation into the variability of Mg II in DESI would help us understand it better. Additionally, the R Mg II -L relation remains poorly constrained (Homayouni et al. 2020; Yu et al. 2023). Therefore, given higher variability aids in more precise RM, we propose that CL-AGNs with significant variation in Mg II can be used as RM targets, even though the BEL is nearly losing (turn-off) the central continuum.

## 4.3. TDE Candidates

In Figure 10, five CL-AGNs display an asymetrical mid-infrared flare in their WISE light-curve, characterized by a rapid increase in brightness followed by a gradual decline. Among these, J015804.75-005221.9, also known as PS16dtm, has garnered attention as a TDE occurring within an AGN system. Blanchard et al. (2017) first reported on this event, highlighting the interaction of stellar debris with the preexisting accretion disk as a key feature. Subsequent investigations by Jiang et al. (2017) elucidated the gradual decline in the mid-infrared light curve, attributing it to a high covering factor and structural changes in pre-existing dust induced by intense radiation.

Furthermore, Jiang et al. (2021) identified dozens of TDE candidates or sporadic gas accretion processes from the mid-infrared light curve, including J151345.77+311125.1, J153310.03+272920.3, and J154357.60+100613.6. Interestingly, J115103.77+530140.6 exhibits a mid-infrared outburst akin to the events above. Additionally, its spectrum in the bright state unveils a prominent Fe II feature, indicative of super-Eddington accretion, a trait often associated with TDEs (Gezari 2021). These CL-AGNs might bear significant associations with TDEs, particularly within low-luminosity AGN systems where selfaccretion processes are not predominant, demonstrating the complicated physical mechanism of CL-AGNs.

## 4.4. Recurring CL-AGN Candidates

The study by Zeltyn et al. (2022) describes a transient (J162829.17+432948.5) that experienced a CL turn-off in the broad H β component within one year and then brightened again in several months. This 'V-shape' trend in the light curve (see Figure 2 in Zeltyn et al. 2022) suggests the possibility of recurring CL events over several months or years. In our sample, we have observed that some CL-AGNs display a 'V-shape' or 'reverse V-shape' in their light curves (Figure 11), similar to J162829.17+432948.5. We have coincidentally identified these CL events at the extreme time point in the light curve, where we obtained a spectrum from DESI for these recurring CL-AGN candidates. Furthermore, the latest photometry in the ZTF light curve indicates a return to the bright state for 'V-shape' or a dim state for 'reverse V-shape', almost consistent with the pseudo-photometry magnitude of the SDSS spectrum. Based on the correlation of luminosity variation between BEL and continuum in Section 4.2, it is expected that the luminosity of the BEL should revert to the state observed in the SDSS spectrum. Our follow-up observation in DESI or fortuitous observation in the TDSS project of SDSS would validate the recurring phenomenon (Blanton et al. 2017; Green et al. 2022; Zeltyn et al. 2024). Further analysis of these recurring candidates will enhance our understanding of the underlying mechanisms driving CL phenomena, mainly by providing an accurate timescale instead of an upper-limit transition between two epochs. Whether recurring CL events in the short timescale can be distinguished from variations of typical AGN is also essential for understanding the AGN population, especially whether CL-AGN BELongs to a distinctive species or a different phase in the evolution of all AGN (Green et al. 2022).

Despite the apparent disappearance of BEL in these CL-AGNs, spectral decomposition reveals persistent broad components indicative of the underlying continuity of CL events. Hence, the classification between CL-AGNs and CL candidates based solely on spectral characteristics in Guo et al. (2024) is insufficient considering DESI and SDSS randomly sample the light curve. For example, a CL-AGN in Figure 11 would be classified as a candidate if the best snapshot time (the most extreme time point) in the light curve is missed in the DESI spectrum. Moreover, this also underscores the importance of considering sampling effects when interpreting the ratio of CL-AGNs, as spectroscopic observations may not always coincide with periods of pronounced variability.

Mag

Mag

12

13

14

15

14

15

16

12 5

.

13 0

.

13 5

.

11 5

.

12 0

.

12 5

.

13 0

.

13 5

.

14 0

.

14 5

.

Mag

Mag

Mag allwise-W1

allwise-W2

neowise-W1

neowise-W2

52000

54000

allwise-W1

allwise-W2

neowise-W1

neowise-W2

52000

54000

allwise-W1

allwise-W2

neowise-W1

neowise-W2

52000

54000

allwise-W1

allwise-W2

neowise-W1

neowise-W2

52000

54000

allwise-W1

allwise-W2

neowise-W1

neowise-W2

52000

54000

56000

MJD

56000

MJD

56000

MJD

56000

MJD

56000

58000

58000

58000

58000

58000

60000

60000

60000

60000

60000

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00402v1-figure-714a20ca-98b8-41eb-9de9-6a00dff9bafa.jpg)

16

1

DESI

115103.77+530140.6

-

˚

A

17

2

-

-

10

cm

1

-

∗

λ

f

erg s

1

-

˚

A

17

2

-

-

10

cm

1

-

∗

λ

f

1

-

˚

A

17

2

-

-

10

cm

1

-

∗

λ

f

erg s

1

-

˚

A

17

2

-

-

10

cm

1

-

∗

λ

f

erg s erg s

14

12

10

8

6

4

2

0

64

56

48

40

32

24

16

8

180

160

140

120

100

80

60

40

20

21

18

15

12

9

6

3

SDSS

Mg ii

3000

DESI

SDSS

4000

DESI

SDSS

4000

DESI

SDSS

4000

H

β

4000

5000

6000

Rest-Frame Wavelength( ˚ )

A

151345.77+311125.1

H

5000

6000

7000

8000

9000

β

H

α

Rest-Frame Wavelength( ˚ )

A

153310.03+272920.3

H

5000

6000

7000

8000

9000

β

H

α

Rest-Frame Wavelength( ˚ )

A

154357.60+100613.6

H

5000

6000

7000

8000

β

H

α

MJD

Rest-Frame Wavelength( ˚ )

A

Figure 10. The mid-infrared light curve and spectra for CL-AGNs with asymmetrical flares. The left panel displays the W1-band and W2-band of MIR light curves from WISE (square) and NEOWISE (diamond). The dashed line represents the spectral MJD for DESI and SDSS. The right panel displays the smoothed SDSS and DESI spectra represented by blue and red lines respectively.

Figure 11. The light curves and spectra for recurring CL-AGN candidates can be seen by re-bright or re-dim features in the latest light curve of ZTF, which is close to the pseudo-photometry magnitude of the SDSS spectrum. The top-left panel displays the optical light curve over 20 years from SDSS (square), CRTS (star), PS1 (triangle), PTF (pentagon), DESI (Diamond), and ZTF (dot) when available. The spectral pseudophotometry is marked with red 'x' markers for DESI and blue 'x' for SDSS. The bottom-left panel displays the W1-band and W2-band of mid-infrared light curves from WISE (square) and NEOWISE (diamond). The right panel is smoothed SDSS and DESI represented by blue and red lines respectively.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00402v1-figure-05cd2c2f-fa4a-4428-920d-5aebc22c3e04.jpg)

Figure 12. The BPT diagram for CL-AGNs and SDSS DR14Q. The violet triangles are the measurements of the CL-AGNs in this work, using the emission line flux ratio of the dim spectra. The grey dots are values from SDSS DR14Q at z &lt; 0 9 . from Rakshit et al. (2020). The upper and right panel represent the distribution of log (N II/H α ) and log ([O III]/H β ) respectively.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00402v1-figure-da25eac6-6a8c-4c32-898b-a282b70cffb2.jpg)

We employ the Baldwin-Phillips-Terlevich (BPT) diagram to delineate the dominant energy sources for CL-AGNs. The classification criteria are determined by empirical relations established by Kauffmann et al. (2003) for Composite and Starforming galaxies, Kewley et al. (2001) for Composite and Seyfert/LINERs (Low Ionization Nuclear Emission Line Regions), and Kewley et al. (2006) for Seyfert and LINERs:

$$\log ( [ O \, \i m ] / H \beta ) & = 0. 6 1 / [ \log ( N \, \i m / H \alpha ) - 0. 0 5 ] + 1. 3,$$

$$\log ( [ O \, \i i ] / H \beta ) & = 0. 6 1 / [ \log ( N \, \i i / H \alpha ) - 0. 3 7 ] + 1. 1 9,$$

$$\log ( [ O \, \mathfrak { m } ] / H \beta = 1. 0 1 \log ( N \, \mathfrak { m } / H \alpha ) + 0. 4 8.$$

Figure 12 compares the CL-AGNs from our study with the SDSS DR14Q sample by Rakshit et al. (2020). Our CL-AGNs mainly fall within the Seyfert region of the BPT diagram, showing no significant statistical difference compared to the SDSS DR14Q sample. However, the median value of N II/ H α does exhibit a slight deviation towards LINERs.

This result is not unexpected, given that the BPT diagram generally represents the average AGN activity over several thousand years (Fischer et al. 2013; Xu &amp; Wang 2022), particularly for the slit or fiber spectrum. Furthermore, the CL event might be a flicker for NLR since most AGNs with low Eddington likely experience recurring CL events only the most extreme time point in the light curve. Under this assumption, the activity level of CL-AGNs remains consistent with that of the overall AGN population in the BPT diagram.

However, a theatrical CL event where an AGN completely transforms into a galaxy or vice versa might significantly impact the NLR after several decades or centuries. The utility of Integral Field Unit (IFU) observations, tracing the long-term variability of AGNs (Raimundo et al. 2019), would have a great chance to capture CL event histories in the BPT diagram. Supposing a dramatic CL event lasts for decades, the variations in the BPT emission-line ratios will produce ring-like structures or stratified cones in IFU images, analogous to the appendix Figure B.1 in Husemann et al. (2022).

## 5. DISCUSSION

To identify some of the most intriguing physical mechanisms in CL-AGNs, especially those requiring additional energy sources beyond changes in their own accretion states, we present additional CL-AGN classifications based on the characteristics of the light curves in Section 5.1 and discuss the corresponding physical mechanisms in Section 5.2.

Figure 13. The light curves and spectra for CL-AGNs that show a peculiar changing-look phase. The legend is the same as Figure 11.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00402v1-figure-82afaeae-ced7-404b-8969-b149c22e42ed.jpg)

Before classifying CL-AGNs, it is essential to review the optical variability of quasars, which remains poorly understood (Giveon et al. 1999; Bauer et al. 2009; Ai et al. 2010). Currently, the variability of AGN light curves can be described by various mathematical models, such as the structure-function, DRW, continuous autoregressive moving-average, and coronaheated accretion-disk reprocessing models (Kelly et al. 2009; MacLeod et al. 2010; Kelly et al. 2014; Kasliwal et al. 2015; Sun et al. 2020). The optical variability studies with these models have shown a clear anti-correlation between variability amplitude and the accretion rate or Eddington ratio (Giveon et al. 1999; Kelly et al. 2009; Zuo et al. 2012; Caplar et al. 2017; S´nchez-S´ aez a et al. 2018; Lu et al. 2019). However, it is limited to around 1% of the Eddington ratio, BELow which the host galaxy significantly dilutes the light curve. Therefore, it needs to be clarified whether such anti-correlation still holds BELow 1% of the Eddington ratio.

Ricci &amp; Trakhtenbrot (2023) have proposed a classification scheme for CL-AGNs based on differences in physical mechanisms: 1) changing-obscuration AGNs, attributed to variations in the line-of-sight column density; and 2) changing-state AGNs, triggered by changes in the accretion rate dominated by the innermost regions of the accretion disk. Here, we propose another classification based on the variability of CL-AGNs to determine whether an additional physical mechanism drives the CL instead of either changes to the obscuring column along the line of sight or accretion-disk instabilities (Kawaguchi et al. 1998; Dexter &amp; Agol 2011).

- · Inherent variability. Firstly, the variability of CL-AGNs is well described by a DRW model (Wang et al. 2024b), which cannot be distinguished from Type 1 AGNs except for larger amplitudes. The variability patterns of most of the CL-AGNs in this work also resembled 20-year photometric light curves in SDSS Stripe 82 (Stone et al. 2022). Afterward, the correlation between BEL and continuum (Section 4.2) and the BPT diagram (Section 4.5) for CL-AGNs show nondistinctive features compared to the general AGN population. A significant difference in CL-AGNs is the lower Eddington ratio ( ∼ 1%) compared to Type 1 AGNs, as discussed in Section 4.1. Additionally, through photometrical or repeated spectral selection (MacLeod et al. 2016; Yan et al. 2019; MacLeod et al. 2019; Graham et al. 2020; Green et al. 2022; Guo et al. 2024; Zeltyn et al. 2024; Wang et al. 2024b), the accumulated number of CL-AGNs has significantly increased (over 1,000 cases). It indicates that the peculiarity might be due to the lack of frequently repeated spectral measurements, especially for CL-AGNs with one BEL change. Overall, we proposed that most CL-AGNs inherently BELong to a subclass characterized by a low Eddington ratio, which shows a more considerable amplitude variability. Frequently repeated observation of a large population with a low Eddington ratio will be conducive to restricting our claim. For distinguishing inherent variability, there is no quiescent or nearly quiescent state in their light curves.
- · Fading AGN. Unlike that random variability, we notice that some CL events may relate to the AGN duty cycle, which ultimately turns a type 1 AGN to a type 2 AGN or galaxy in the observed spectrum. Such drastic change is unlikely to be explained by accretion disk instability without additional physical mechanisms, which might stem from fundamental transformations of the accretion disk or the cessation of AGN activity. Among these fading AGNs, we propose a peculiar changing-look (PCL) phase whose light curve is characterized by a slow, decades-long decline in luminosity and the eventual disappearance of all BEL. Figure 13 showcases four candidates with PCL phase in our sample. The behavior of the PCL phase in the light curve is also unexpected in the DRW model since the characteristic time scale for a typical AGNcan be constrained to be about 10 years (Kozłowski 2017; Stone et al. 2022), and the DRW model has less possibility of maintaining regular pattern (or even unchanged after PCL phase) in a type 1 AGN. In addition, we also conjugate a typical AGN light curve from SDSS Stripe 82 with a PCL light curve to illustrate this phase in Figure 14. Such conjugated patterns are similar to the variability of Mrk 590 and Mrk 1018, exhibit decade-long drops, and enter a state of quiescent (Denney et al. 2014; Noda &amp; Done 2018). However, the fading AGN after the PCL phase might remain in a long-term non-accretion state (Case I without observable variability or re-ignite after several years or decades (Case II), which is still unclear. If the PCL phase represents a distinct phenomenon, understanding how long and why they remain inactive is vital for understanding the origin of the AGN variability or AGN duty cycle. Such
- · Transient flares or outbursts. Analogous to TDE candidates shown in Figure 10, some CL-AGNs display minimal variability over long periods, which might be diluted by the host galaxy. However, sudden flares or outbursts in the light curve occur during a CL event. Similar variability patterns have been observed in additional examples from NGC 2617 and six CL LINERs (Shappee et al. 2014; Frederick et al. 2019). These flares are distinct from the typical variation of type 1 AGN, potentially indicating unique accretion processes or transients near SMBHs. Supposing the Eddington ratio and variability amplitude are still anticorrelated BELow 1% of the Eddington ratio, the CL event can manifest itself by drastic continuum variations since the accretion disk instability may reach extreme levels.

Figure 14. The peculiar changing-look phase transfers a powerful AGN into a fading AGN. The green points depict a typical light curve for a type 1 AGN from SDSS Stripe 82 over a span of 20 years (Stone et al. 2022). The red points are the light curve from J122638.66-001114.0. The two arrows indicate the two possible patterns after the peculiar changing-look state: Case I refers to the AGN fading or switching off, while Case II refers to the AGN being re-triggered. The upper-left (right) panel show the spectra of J122638.66-001114.0 from the SDSS (DESI) spectrum to represents the type 1 AGN state (type 2 AGN state or galaxy state).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00402v1-figure-36901982-fb56-40eb-871e-1efa15b1f51b.jpg)

As noted in several previous studies (e.g., LaMassa et al. 2015; Runnoe et al. 2016; Green et al. 2022), changes in the intrinsic accretion rate, rather than cloud obscuration, are likely the primary cause for the majority of CL events. Initial evidence against obscuration suggests that the timescale for a cloud to move outside the BLR would need several decades to produce a CL event (LaMassa et al. 2015; Sheng et al. 2017; Stern et al. 2018), which contradicts the observed CL timescale of several months to years. Furthermore, the deredden test does not support the obscuration scenario, as the extinction from the continuum is not consistent with with the emission-line strengths of the BEL (LaMassa et al. 2015; Ruan et al. 2016). X-ray observations indicate that the dramatic change in the continuum is primarily driven by the central ionization source (Noda &amp; Done 2018), and the critical value for CL (1% Eddington ratio) is generally consistent with a reversal of the UV-to-X-ray spectral index α ox trend (Ruan et al. 2019; Jin et al. 2021; Yang et al. 2023). Currently, accretion-state changes are likely still the dominant reason for the majority of CL-AGNs, directly linking to the inherent variability in our classification. However, obscuration may be due to the average change in the dust covering factor rather than a single cloud, and the extinction curve is model-dependent and varies depending on dust components (Laor &amp; Draine 1993). As demonstrated by Zeltyn et al. (2022), some short-timescale CL events, which appear as extreme points in the light curve, may be influenced by line-of-sight obscuration. Therefore, it is impossible to rule out the possibility of obscuration without analysis of multi-wavelength data, especially in the mid-infrared and X-ray emission, which will be the focus of our future work.

The physical origins for fading AGNs and transient flares might be very complex but related to the change in structure of accretions disks. The first possibility is an inhomogeneous accretion disk or accretion state change, which might be also related to AGN triggering or quenching. When an AGN is in the lower Eddington ratio state ( &lt; 1% ), and the inner part of the accretion disk is unstable, an accretion rate change from temporary gas supply would exceed the variation predicted by the disk instability model and reionize the BLR (Kawaguchi et al. 1998). Wang et al. (2024a) also suggests that the characteristic time of Bondi accretion under such circumstances could align with the CL transition timescale.

In addition, external physical mechanisms beyond the AGN system might affect the accretion rate, such as nuclear supernovae and TDEs. Since recognizing TDEs and supernovae is challenging due to the inherent variation of AGNs, most TDEs have been identified in inactive galaxies. As indicated by TDE candidates in Figure 10, TDEs (or nuclear supernovae) may be responsible for triggering some CL events. Future research will focus on recognizing TDEs or supernovae in AGNs to investigate the possibility of these occurrences.

Stellar-mass binary black hole mergers in AGN and SMBH binaries can also be considered external physical mechanisms that drive CL events. As suggested by Graham et al. (2017) and Wang et al. (2021), stellar-mass binary black hole mergers could produce significant AGN flares and enhance the continuum luminosity, potentially serving as electromagnetic counterparts identified by the Laser Interferometer Gravitational-Wave Observatory (LIGO). Graham et al. (2023) also investigated unusual

AGNflaring events in ZTF as feasible counterparts to gravitational wave events. From the perspective of galactic evolution, AGNs or quasars may predominantly occur in the final stage of galactic mergers where the SMBH binary is often merged (Hopkins et al. 2008). In this process, dramatic changes in the dust covering or accretion rate enhances the activity of AGN and are accompanied by CL events. Wang &amp; Bon (2020) also determined the probability the probability of CL events caused by the orbit change of the SMBH binary. One of our CL-AGN, J143016.06+23044.5, was proposed as a candidate for a binary black hole merger system (Jiang et al. 2022). Hence, CL-AGNs can be regarded as promising candidates for the electromagnetic counterpart of gravitational waves that the Laser Interferometer Space Antenna (LISA) could detect in the future (Amaro-Seoane et al. 2017).

## 6. SUMMARY

In the analysis of DESI DR1 and SDSS DR16 datasets, we created a catalog with 561 CL-AGNs, which yields several key findings:

- · The ratio between turn-on and turn-off events is 283:278 in our sample.
- · The average Eddington ratio for CL events is observed to be around log λ Edd ∼ -2 .
- · A strong correlation is detected between changes in the luminosity of BEL and continuum luminosity variations.
- · Five CL-AGNs exhibit asymmetrical mid-infrared flares, possibly related to TDEs.
- · The BPT diagram reveals no considerable discrepancy between CL-AGNs and SDSS DR14Q.

We also provide a new classification based on the variability of CL-AGNs, which will aid in delving into the physical mechanism behind CL events for forthcoming investigations. We are conducting a comprehensive statistical analysis of X-ray, midinfrared, and radio properties to elucidate the effect of obscuration and change of accretion ratio in CL-AGNs. Additional host galaxy analysis is more critical in CL-AGNs as it sheds light on co-evolutionary processes. We are committed to unraveling the intricate relationship between AGN phenomena and their host environments by acknowledging the indispensable role of host galaxies in shaping AGN activity.

## ACKNOWLEDGEMENTS

We thank Linhua Jiang, Peter Clark, Suijian Xue, Jian-Min Wang, Cheng Cheng, and Zhi-Xiang Zhang for very helpful discussions and comments. We also thank DESI PubBoard Handler (Antonella Palmese) to provide timely help. The work is supported by the funding supports from the National Key R&amp;D Program of China (grant Nos. 2023YFA1607800 and 2022YFA1602902) and Strategic Priority Research Program of the Chinese Academy of Sciences with Grant Nos. XDB0550100 and XDB0550000. The authors also acknowledge the supports from the National Natural Science Foundation of China (NSFC; grant Nos. 12120101003, 12373010, 12173051, and 12233008), the National Key R&amp;D Program of China (2023YFA1607804, 2023YFA1608100, and 2023YFF0714800), Beijing Municipal Natural Science Foundation (grant No. 1222028), and China Manned Space Project with Nos. CMS-CSST-2021-A02, CMS-CSST-2021-A04 and CMS-CSST-2021-A05. SZ acknowledges support from the National Science Foundation of China (no. 12303011) . This work has been supported by the Polish National Agency for Academic Exchange (Bekker grant BPN/BEK/2021/1/00298/DEC/1), the European Union's Horizon 2020 Research and Innovation programme under the Maria Sklodowska-Curie grant agreement (No. 754510).

This material is based upon work supported by the U.S. Department of Energy (DOE), Office of Science, Office of High-Energy Physics, under Contract No. DE-AC02-05CH11231, and by the National Energy Research Scientific Computing Center, a DOE Office of Science User Facility under the same contract. Additional support for DESI was provided by the U.S. National Science Foundation (NSF), Division of Astronomical Sciences under Contract No. AST-0950945 to the NSF's National Optical-Infrared Astronomy Research Laboratory; the Science and Technology Facilities Council of the United Kingdom; the Gordon and Betty Moore Foundation; the Heising-Simons Foundation; the French Alternative Energies and Atomic Energy Commission (CEA); the National Council of Humanities, Science and Technology of Mexico (CONAHCYT); the Ministry of Science, Innovation and Universities of Spain (MICIU/AEI/10.13039/501100011033), and by the DESI Member Institutions: https://www.desi.lbl.gov/ collaborating-institutions.

The DESI Legacy Imaging Surveys consist of three individual and complementary projects: the Dark Energy Camera Legacy Survey (DECaLS), the Beijing-Arizona Sky Survey (BASS), and the Mayall z-band Legacy Survey (MzLS). DECaLS, BASS and MzLS together include data obtained, respectively, at the Blanco telescope, Cerro Tololo Inter-American Observatory, NSF's NOIRLab; the Bok telescope, Steward Observatory, University of Arizona; and the Mayall telescope, Kitt Peak National Observatory, NOIRLab. NOIRLab is operated by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation. Pipeline processing and analyses of the data were supported by

Table 1. Fits Catalog Description and Column Information of CL-AGNs.

| Number   | Column Name    | Format      | Unit                                                              | Description                                                          |
|----------|----------------|-------------|-------------------------------------------------------------------|----------------------------------------------------------------------|
| HDU1     |                |             |                                                                   | Information of DESI spectrum                                         |
| 1        | TARGETID       | int         |                                                                   | Target ID                                                            |
| 2        | R.A.           | float       | degree                                                            | Right Ascension (J2000)                                              |
| 3        | Dec.           | float       | degree                                                            | Declination (J2000)                                                  |
| 4        | Redshift       | float       |                                                                   | Redshift                                                             |
| 5        | SURVEY         | int         |                                                                   | Survey ID                                                            |
| 6        | PROGRAM        | string      |                                                                   | Program ID                                                           |
| 7        | HEALPIX        | int         |                                                                   | HEALPix ID                                                           |
| 8        | MJD            | float       |                                                                   | Mean MJD of the coadded spectra                                      |
| 9        | LOG L5100      | float       | erg s - 1                                                         | Logarithmic Continuum luminosity at 5100A ˚                          |
| 10       | LOG L5100 ERR  | float       | erg s - 1                                                         | Error of LOG L5100                                                   |
| 11       | LOG L1350      | float       | erg s - 1                                                         | Logarithmic Extrapolated luminosity at 2500A ˚                       |
| 12       | LOG L1350 ERR  | float       | erg s - 1                                                         | Error of LOG L2500                                                   |
| 13       | LOG LMgII      | float       | erg s - 1                                                         | Logarithmic luminosity of Mg II broad component                      |
| 14       | LOG LMgII ERR  | float       | erg s - 1                                                         | Error of LOG LMgII                                                   |
| 15       | LOG LHB        | float       | erg s - 1                                                         | Logarithmic luminosity of Hbeta broad component                      |
| 16       | LOG LHB ERR    | float       | erg s - 1                                                         | Error of LOG LHB                                                     |
| 17       | LOG LHA        | float       | erg s - 1                                                         | Logarithmic luminosity of Halpha broad component                     |
| 18       | LOG LHA ERR    | float       | erg s - 1                                                         | Error of LOG LHA                                                     |
| 19       | LOG MBH HB     | float       | M ⊙                                                               | Black hole mass estimation from Hbeta                                |
| 20       | LOG MBH HB ERR | float       | M ⊙                                                               | Error of LOG MBH HB                                                  |
| 21       | LOG MBH HA     | float       | M ⊙                                                               | Black hole mass estimation from Halpha                               |
| 22       | LOG MBH HA ERR | float       | M ⊙                                                               | Error of LOG MBH HA                                                  |
| 23       | WAVELENGTH     | float       | ˚ A                                                               | One-dimensional spectral wavelength                                  |
| 24       | FLUX           | float       | 10 - 17 erg s - 1 cm - 2 ˚ A - 1 10 - 17 erg s - 1 cm - 2 ˚ A - 1 | One-dimensional spectral flux                                        |
| HDU2     |                |             |                                                                   | Information of SDSS spectrum                                         |
| 1        | SDSS NAME      | string      |                                                                   | SDSS object                                                          |
| 2        | R.A.           | float       | degree                                                            | Right Ascension (J2000)                                              |
| 3        | Dec.           | float       | degree                                                            | Declination (J2000)                                                  |
| 4        | Redshift       | float       |                                                                   | Redshift                                                             |
| 5        | PLATE          | int         |                                                                   | Plate ID                                                             |
| 6        | MJD            | int         |                                                                   | MJD of SDSS spectrum                                                 |
| 7        | FIBERID        | int         |                                                                   | Fiber ID                                                             |
| 8        | Transition     | string      |                                                                   | Changing state between SDSS and DESI                                 |
| 9        | LOG L5100      | float       | erg s - 1                                                         | Continuum luminosity at 5100A ˚                                      |
| 10       | LOG L5100 ERR  | float       | erg s - 1                                                         | Error of LOG L5100                                                   |
| 11       | LOG L1350      | float       | erg s - 1                                                         | Extrapolated luminosity at 2500A ˚                                   |
| 12       | LOG L1350 ERR  | float       | erg s - 1                                                         | Error of LOG L2500                                                   |
| 13       | LOG LMgII      | float       | erg s - 1                                                         | Logarithmic luminosity of Mg II broad component                      |
| 14       | LOG LMgII ERR  | float       | erg s - 1                                                         | Error of LOG LMgII                                                   |
| 15       | LOG LHB        | float       | erg s - 1                                                         | Logarithmic luminosity of Hbeta broad component                      |
| 16       | LOG LHB ERR    | float       | erg s - 1                                                         | Error of LOG LHB                                                     |
| 17       | LOG LHA        | float       | erg s - 1                                                         | Logarithmic luminosity of Halpha broad component                     |
| 18       | LOG LHA ERR    | float       | erg s - 1                                                         | Error of LOG LHA                                                     |
| 19       | LOG MBH HB     | float       | M ⊙                                                               | Black hole mass estimation from Hbeta                                |
| 20       | LOG MBH HB ERR | float       | M ⊙                                                               | Error of LOG MBH HB                                                  |
| 21       | LOG MBH HA     | float       | M ⊙                                                               | Black hole mass estimation from Halpha Error of LOG MBH HA           |
| 22       | LOG MBH HA ERR | float       | M ⊙                                                               |                                                                      |
| 23       | WAVELENGTH     | float       | ˚ A                                                               | One-dimensional spectral wavelength                                  |
| 24 25    | FLUX FLUX ERR  | float float | 10 - 17 erg s - 1 cm - 2 ˚ A - 1 10 - 17 erg s - 1 cm - 2 ˚ A - 1 | One-dimensional spectral flux One-dimensional error of spectral flux |

NOTE-This table is available in its entirety in the machine-readable format.

NOIRLab and the Lawrence Berkeley National Laboratory. Legacy Surveys also uses data products from the Near-Earth Object Wide-field Infrared Survey Explorer (NEOWISE), a project of the Jet Propulsion Laboratory/California Institute of Technology, funded by the National Aeronautics and Space Administration. Legacy Surveys was supported by: the Director, Office of Science, Office of High Energy Physics of the U.S. Department of Energy; the National Energy Research Scientific Computing Center, a DOE Office of Science User Facility; the U.S. National Science Foundation, Division of Astronomical Sciences; the National Astronomical Observatories of China, the Chinese Academy of Sciences and the Chinese National Natural Science Foundation. LBNL is managed by the Regents of the University of California under contract to the U.S. Department of Energy. The complete acknowledgments can be found at https://www.legacysurvey.org/.

Any opinions, findings, and conclusions or recommendations expressed in this material are those of the author(s) and do not necessarily reflect the views of the U. S. National Science Foundation, the U. S. Department of Energy, or any of the listed funding agencies.

The authors are honored to be permitted to conduct scientific research on Iolkam Du'ag (Kitt Peak), a mountain with particular significance to the Tohono O'odham Nation.

SDSS is managed by the Astrophysical Research Consortium for the Participating Institutions of the SDSS Collaboration including the Brazilian Participation Group, the Carnegie Institution for Science, Carnegie Mellon University, Center for Astrophysics - Harvard &amp; Smithsonian, the Chilean Participation Group, the French Participation Group, Instituto de Astrof´ ısica de Canarias, The Johns Hopkins University, Kavli Institute for the Physics and Mathematics of the Universe(IPMU)/University of Tokyo, the Korean Participation Group, Lawrence Berkeley National Laboratory, Leibniz Institut f¨rAstrophysik Potsdam (AIP), u Max-Planck-Institut f¨r Astronomie (MPIA Heidelberg), Max-Planck-Institut f¨ ur Astrophysik (MPA Garching), Max-Plancku Institut f¨r ExtraterrestrischePhysik (MPE), National Astronomical Observatories of China, New Mexico State University, New u York University, University of Notre Dame, Observat´ ario Nacional/MCTI, The Ohio State University, Pennsylvania State University, Shanghai Astronomical Observatory, United Kingdom Participation Group, Universidad Nacional Aut´ onoma de M´ exico, University of Arizona, University of Colorado Boulder, University of Oxford, University of Portsmouth, University of Utah, University of Virginia, University of Washington, University of Wisconsin, Vanderbilt University, and Yale University.

The Catalina Sky Survey is funded by the National Aeronautics and Space Administration under Grant No. NNG05GF22G issued through the Science Mission Directorate Near-Earth Objects Observations Program. The CRTS survey is supported by the US National Science Foundation under grants AST-0909182 and AST-1313422. The CRTS survey is supported by the US National Science Foundation under grants AST-0909182 and AST-1313422. The PS1 has been made possible through contributions by the Institute for Astronomy, the University of Hawaii, the Pan-STARRS Project Office, the Max-Planck Society and its participating in- stitutes, the Max Planck Institute for Astronomy, Heidelberg and the Max Planck Institute for Extraterrestrial Physics, Garching, The Johns Hopkins University, Durham University, the University of Edinburgh, Queen's University BELfast, the Harvard-Smithsonian Center for Astrophysics, the Las Cumbres Observatory Global Telescope Network Incorporated, the National Central University of Taiwan, the Space Telescope Science Institute, the National Aeronautics and Space Administration under Grant No. NNX08AR22G issued through the Planetary Science Division of the NASA Science Mission Directorate, the National Science Foundation under Grant No. AST-1238877, the University of Maryland, and Eotvos Lorand University (ELTE). PTF image data are obtained with the Samuel Oschin Telescope and the 60 inch Telescope at the Palomar Observatory as part of the Palomar Transient Factory project, a scientific collaboration between the California Insti- tute of Technology, Columbia University, Las Cumbres Obser- vatory, the Lawrence Berkeley National Laboratory, the National Energy Research Scientific Computing Center, the University of Oxford, and the Weizmann Institute of Science. ZTF is supported by the National Science Foundation under Grant No. AST-2034437 and a collaboration including Caltech, IPAC, the Weizmann Institute for Science, the Oskar Klein Center at Stockholm University, the University of Maryland, Deutsches Elektronen-Synchrotron and Humboldt University, the TANGO Consortium of Taiwan, the University of Wis- consin at Milwaukee, Trinity College Dublin, Lawrence Livermore National Laboratories, and IN2P3, France. Operations are conducted by COO, IPAC, and UW.

## REFERENCES

Abareshi, B., Aguilar, J., Ahlen, S., et al. 2022, AJ, 164, 207,

doi: 10.3847/1538-3881/ac882b

Abazajian, K., Adelman-McCarthy, J. K., Ag¨ ueros, M. A., et al.

2004, AJ, 128, 502, doi: 10.1086/421365

Abazajian, K. N., Adelman-McCarthy, J. K., Ag¨ ueros, M. A., et al.

2009, ApJS, 182, 543, doi: 10.1088/0067-0049/182/2/543

Abdo, A. A., Ackermann, M., Agudo, I., et al. 2010, ApJ, 716, 30, doi: 10.1088/0004-637X/716/1/30

Adelman-McCarthy, J. K., Ag¨ ueros, M. A., Allam, S. S., et al.

2008, ApJS, 175, 297, doi: 10.1086/524984

Ahn, C. P., Alexandroff, R., Allende Prieto, C., et al. 2012, ApJS,

203, 21, doi: 10.1088/0067-0049/203/2/21

- Ahumada, R., Prieto, C. A., Almeida, A., et al. 2020, ApJS, 249, 3, doi: 10.3847/1538-4365/ab929e
- Ai, Y. L., Yuan, W., Zhou, H. Y., et al. 2010, ApJL, 716, L31, doi: 10.1088/2041-8205/716/1/L31
- Alam, S., Ata, M., Bailey, S., et al. 2017, MNRAS, 470, 2617, doi: 10.1093/mnras/stx721
- Alexander, D. M., &amp; Hickox, R. C. 2012, NewAR, 56, 93, doi: 10.1016/j.newar.2011.11.003
- Alexander, D. M., Davis, T. M., Chaussidon, E., et al. 2023, AJ, 165, 124, doi: 10.3847/1538-3881/acacfc

Amaro-Seoane, P., Audley, H., Babak, S., et al. 2017, arXiv

- e-prints, arXiv:1702.00786, doi: 10.48550/arXiv.1702.00786 Antonucci, R. 1993, ARA&amp;A, 31, 473,

doi: 10.1146/annurev.aa.31.090193.002353

Bailey et al. 2023, in preparation

Bauer, A., Baltay, C., Coppi, P., et al. 2009, ApJ, 696, 1241, doi: 10.1088/0004-637X/696/2/1241

- Bentz, M. C., Denney, K. D., Grier, C. J., et al. 2013, ApJ, 767, 149, doi: 10.1088/0004-637X/767/2/149
- Blanchard, P. K., Nicholl, M., Berger, E., et al. 2017, ApJ, 843, 106, doi: 10.3847/1538-4357/aa77f7
- Blandford, R. D., &amp; Znajek, R. L. 1977, MNRAS, 179, 433, doi: 10.1093/mnras/179.3.433
- Blanton, M. R., Bershady, M. A., Abolfathi, B., et al. 2017, AJ, 154, 28, doi: 10.3847/1538-3881/aa7567

Boroson, T. A., &amp; Green, R. F. 1992, ApJS, 80, 109, doi: 10.1086/191661

- Brodzeller, A., Dawson, K., Bailey, S., et al. 2023, arXiv e-prints, arXiv:2305.10426, doi: 10.48550/arXiv.2305.10426
- Bruzual, G., &amp; Charlot, S. 2003, MNRAS, 344, 1000, doi: 10.1046/j.1365-8711.2003.06897.x
- Cao, X., You, B., &amp; Wei, X. 2023, MNRAS, 526, 2331, doi: 10.1093/mnras/stad2877
- Caplar, N., Lilly, S. J., &amp; Trakhtenbrot, B. 2017, ApJ, 834, 111, doi: 10.3847/1538-4357/834/2/111
- Chambers, K. C., Magnier, E. A., Metcalfe, N., et al. 2016, arXiv e-prints, arXiv:1612.05560, doi: 10.48550/arXiv.1612.05560
- Chan, C.-H., Piran, T., &amp; Krolik, J. H. 2020, ApJ, 903, 17, doi: 10.3847/1538-4357/abb776
- Chaussidon, E., Y` eche, C., Palanque-Delabrouille, N., et al. 2022, arXiv e-prints, arXiv:2208.08511.

https://arxiv.org/abs/2208.08511

- Cho, H., Woo, J.-H., Wang, S., et al. 2023, ApJ, 953, 142, doi: 10.3847/1538-4357/ace1e5
- Denney, K. D., De Rosa, G., Croxall, K., et al. 2014, ApJ, 796, 134, doi: 10.1088/0004-637X/796/2/134
- DESI Collaboration, Aghamousa, A., Aguilar, J., et al. 2016a, arXiv e-prints, arXiv:1611.00036.

https://arxiv.org/abs/1611.00036

| -. 2016b, arXiv e-prints, arXiv:1611.00037, doi: 10.48550/arXiv.1611.00037                                                    |
|-------------------------------------------------------------------------------------------------------------------------------|
| DESI Collaboration, Abareshi, B., Aguilar, J., et al. 2022, AJ, 164, 207, doi: 10.3847/1538-3881/ac882b                       |
| DESI Collaboration, Adame, A. G., Aguilar, J., et al. 2023a, arXiv e-prints, arXiv:2306.06308, doi: 10.48550/arXiv.2306.06308 |
| -. 2023b, arXiv e-prints, arXiv:2306.06307,                                                                                   |
| doi: 10.48550/arXiv.2306.06307 -. 2024a, arXiv e-prints, arXiv:2404.03002, doi: 10.48550/arXiv.2404.03002                     |
| -. 2024b, in preparation                                                                                                      |
| -. 2024c, in preparation                                                                                                      |
| -. 2024d, arXiv e-prints, arXiv:2404.03000,                                                                                   |
| doi: 10.48550/arXiv.2404.03000                                                                                                |
| -. 2024e, arXiv e-prints, arXiv:2404.03001,                                                                                   |
| doi: 10.48550/arXiv.2404.03001                                                                                                |
| -. 2024f, in preparation                                                                                                      |
| -. 2024g, in preparation                                                                                                      |
| -. 2024h, in preparation                                                                                                      |
| Dexter, J., &Agol, E. 2011, ApJL, 727, L24, doi: 10.1088/2041-8205/727/1/L24                                                  |
| Dey, A., Schlegel, D. J., Lang, D., et al. 2019, AJ, 157, 168, doi: 10.3847/1538-3881/ab089d                                  |
| Drake, A. J., Djorgovski, S. G., Mahabal, A., et al. 2009, ApJ, 696, 870, doi: 10.1088/0004-637X/696/1/870                    |
| Du, P., Brotherton, M. S., Wang, K., et al. 2018, ApJ, 869, 142, doi: 10.3847/1538-4357/aaed2c                                |
| Eisenstein, D. J., Weinberg, D. H., Agol, E., et al. 2011, AJ, 142, 72, doi: 10.1088/0004-6256/142/3/72                       |
| Elitzur, M., Ho, L. C., &Trump, J. R. 2014, MNRAS, 438,                                                                       |
| 3340, doi: 10.1093/mnras/stt2445                                                                                              |
| T. C., Crenshaw, D. M., Kraemer, S. B., &Schmitt, H. R. 2013, ApJS, 209, 1, doi: 10.1088/0067-0049/209/1/1                    |
| Fischer,                                                                                                                      |
| 10.1086/316293 B., Diehl, H. T., Honscheid, K., et al. 2015, AJ, 150, doi: 10.1088/0004-6256/150/5/150                        |
| 150,                                                                                                                          |
| Flaugher, 883, 31,                                                                                                            |
| M., Ichikawa, T., Gunn, J. E., et al. 1996, AJ, 111, 1748,                                                                    |
| Fukugita, doi: 10.1086/117915                                                                                                 |
| doi: 10.1146/annurev-astro-111720-030029 Gezari, S., Hung, T., Cenko, S. B., et al. 2017, ApJ, 835,                           |
| 144, doi: 10.3847/1538-4357/835/2/144                                                                                         |
| Giveon, U., Maoz, D., Kaspi, S., Netzer, H., &Smith, P. S. 1999, MNRAS, 306, 637, doi: 10.1046/j.1365-8711.1999.02556.x       |
| Graham, M. J., Djorgovski, S. G., Drake, A. J., et al. 2017,                                                                  |
| MNRAS, 470, 4112, doi: 10.1093/mnras/stx1456                                                                                  |
| Fitzpatrick, E. L. 1999, PASP, 111,                                                                                           |
| Farr, J., Font-Ribera, A., doi:                                                                                               |
| Frederick, S., Gezari, S., Graham,                                                                                            |
| 2020, 015,                                                                                                                    |
| &Pontzen,                                                                                                                     |
| A.                                                                                                                            |
| 10.1088/1475-7516/2020/11/015                                                                                                 |
| 2020,                                                                                                                         |
| JCAP,                                                                                                                         |
| Gezari, S. 2021, ARA&A, 59, 21,                                                                                               |
| 63, doi:                                                                                                                      |
| M. J., et al. 2019, ApJ, doi: 10.3847/1538-4357/ab3a38                                                                        |

| Graham, M. J., Ross, N. P., Stern, D., et al. 2020, MNRAS, 491, 4925, doi: 10.1093/mnras/stz3244                          |
|---------------------------------------------------------------------------------------------------------------------------|
| Graham, M. J., McKernan, B., Ford, K. E. S., et al. 2023, ApJ, 942, 99, doi: 10.3847/1538-4357/aca480                     |
| Green, P. J., Pulgarin-Duque, L., Anderson, S. F., et al. 2022, ApJ, 933, 180, doi: 10.3847/1538-4357/ac743f              |
| Gunn, J. E., Carr, M., Rockosi, C., et al. 1998, AJ, 116, 3040, doi: 10.1086/300645                                       |
| Gunn, J. E., Siegmund, W. A., Mannery, E. J., et al. 2006, AJ, 131, 2332, doi: 10.1086/500975                             |
| Guo, H., Peng, J., Zhang, K., et al. 2020, ApJ, 905, 52, doi: 10.3847/1538-4357/abc2ce                                    |
| Guo, W.-J., Li, Y.-R., Zhang, Z.-X., Ho, L. C., &Wang, J.-M. 2022, ApJ, 929, 19, doi: 10.3847/1538-4357/ac4e84            |
| Guo, W.-J., Zou, H., Fawcett, V. A., et al. 2024, ApJS, 270, 26, doi: 10.3847/1538-4365/ad118a                            |
| Guy, J., Bailey, S., Kremin, A., et al. 2022, arXiv e-prints, arXiv:2209.14482. https://arxiv.org/abs/2209.14482          |
| Hahn, C., Wilson, M. J., Ruiz-Macias, O., et al. 2022, arXiv e-prints, arXiv:2208.08512. https://arxiv.org/abs/2208.08512 |
| Heckman, T. M., Gonz´lez-Delgado, a R., Leitherer, C., et al. 1997, ApJ, 482, 114, doi: 10.1086/304139                    |
| Ho, L. C. 2008, ARA&A, 46, 475, doi: 10.1146/annurev.astro.45.051806.110546                                               |
| Homayouni, Y., Trump, J. R., Grier, C. J., et al. 2020, ApJ, 901, 55, doi: 10.3847/1538-4357/ababa9                       |
| Hon, W. J., Wolf, C., Onken, C. A., Webster, R., &Auchettl, K. 2022, MNRAS, 511, 54, doi: 10.1093/mnras/stab3694          |
| Hopkins, P. F., Hernquist, L., Cox, T. J., &Kereˇ, s D. 2008, ApJS, 175, 356, doi: 10.1086/524362                         |
| Hu, C., Du, P., Lu, K.-X., et al. 2015, ApJ, 804, 138, doi: 10.1088/0004-637X/804/2/138                                   |
| Husemann, B., Singha, M., Scharw¨chter, a J., et al. 2022, A&A, 659, A124, doi: 10.1051/0004-6361/202141312               |
| Jiang, N., Wang, T., Yan, L., et al. 2017, ApJ, 850, 63, doi: 10.3847/1538-4357/aa93f5                                    |
| Jiang, N., Wang, T., Dou, L., et al. 2021, ApJS, 252, 32, doi: 10.3847/1538-4365/abd1dc                                   |
| Jiang, N., Yang, H., Wang, T., et al. 2022, arXiv e-prints, arXiv:2201.11633, doi: 10.48550/arXiv.2201.11633              |
| Jin, J.-J., Wu, X.-B., &Feng, X.-T. 2022, ApJ, 926, 184, doi: 10.3847/1538-4357/ac410c                                    |
| Jin, X., Ruan, J. J., Haggard, D., et al. 2021, ApJ, 912, 20, doi: 10.3847/1538-4357/abeb17                               |
| Juneau, S., Canning, R., Alexander, D. M., et al. 2024, arXiv e-prints, arXiv:2404.03621, doi: 10.48550/arXiv.2404.03621  |
| Kasliwal, V. P., Vogeley, M. S., &Richards, G. T. 2015, MNRAS, 451, 4328, doi: 10.1093/mnras/stv1230                      |
| Kauffmann, G., Heckman, T. M., Tremonti, C., et al. 2003,                                                                 |
| MNRAS, 346, 1055, doi: 10.1111/j.1365-2966.2003.07154.x                                                                   |

| Kawaguchi, T., Mineshige, S., Umemura, M., &Turner, E. L. 1998, ApJ, 504, 671, doi: 10.1086/306105                                                                                                  |
|-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| Kelly, B. C., Bechtold, J., &Siemiginowska, A. 2009, ApJ, 698, 895, doi: 10.1088/0004-637X/698/1/895                                                                                                |
| Kelly, B. C., Becker, A. C., Sobolewska, M., Siemiginowska, A., & Uttley, P. 2014, ApJ, 788, 33, doi: 10.1088/0004-637X/788/1/33                                                                    |
| Kewley, L. J., Groves, B., Kauffmann, G., &Heckman, T. 2006, MNRAS, 372, 961, doi: 10.1111/j.1365-2966.2006.10859.x                                                                                 |
| Kewley, L. J., Heisler, C. A., Dopita, M. A., &Lumsden, S. 2001, ApJS, 132, 37, doi: 10.1086/318944                                                                                                 |
| Kormendy, J., &Ho, L. C. 2013, ARA&A, 51, 511, doi: 10.1146/annurev-astro-082708-101811                                                                                                             |
| Kozłowski, S. 2017, A&A, 597, A128, doi: 10.1051/0004-6361/201629890                                                                                                                                |
| LaMassa, S. M., Cales, S., Moran, E. C., et al. 2015, ApJ, 800, 144, doi: 10.1088/0004-637X/800/2/144                                                                                               |
| Lan, T.-W., Tojeiro, R., Armengaud, E., et al. 2023, ApJ, 943, 68, doi: 10.3847/1538-4357/aca5fa                                                                                                    |
| Laor, A., &Draine, B. T. 1993, ApJ, 402, 441, doi: 10.1086/172149                                                                                                                                   |
| Law, N. M., Kulkarni, S. R., Dekany, R. G., et al. 2009, PASP, 121, 1395, doi: 10.1086/648598                                                                                                       |
| Levi, M., Bebek, C., Beers, T., et al. 2013, arXiv e-prints, arXiv:1308.0847, doi: 10.48550/arXiv.1308.0847                                                                                         |
| Li, R., Ho, L. C., Ricci, C., et al. 2022, ApJ, 933, 70, doi: 10.3847/1538-4357/ac714a                                                                                                              |
| L´pez-Navas, o E., Mart´nez-Aldama, ı M. L., Bernal, S., et al. 2022, MNRAS, 513, L57, doi: 10.1093/mnrasl/slac033                                                                                  |
| Lu, K.-X., Huang, Y.-K., Zhang, Z.-X., et al. 2019, ApJ, 877, 23, doi: 10.3847/1538-4357/ab16e8                                                                                                     |
| Lyke, B. W., Higley, A. N., McLane, J. N., et al. 2020, ApJS, 250, 8, doi: 10.3847/1538-4365/aba623                                                                                                 |
| MacLeod, C. L., Ivezi´, c ˇ Z., Kochanek, C. S., et al. 2010, ApJ, 721, 1014, doi: 10.1088/0004-637X/721/2/1014                                                                                     |
| MacLeod, C. L., Ross, N. P., Lawrence, A., et al. 2016, MNRAS, 457, 389, doi: 10.1093/mnras/stv2997                                                                                                 |
| MacLeod, C. L., Green, P. J., Anderson, S. F., et al. 2019, ApJ, 874, 8, doi: 10.3847/1538-4357/ab05e2                                                                                              |
| Mainzer, A., Bauer, J., Cutri, R. M., et al. 2014, ApJ, 792, 30, doi: 10.1088/0004-637X/792/1/30                                                                                                    |
| Masci, F. J., Laher, R. R., Rusholme, B., et al. 2019, PASP, 131, 018003, doi: 10.1088/1538-3873/aae8ac                                                                                             |
| Miller, T. N., Doel, P., Gutierrez, G., et al. 2023, arXiv e-prints, arXiv:2306.06310, doi: 10.48550/arXiv.2306.06310                                                                               |
| Moustakas, J., Scholte, D., Dey, B., &Khederlarian, A. 2023, FastSpecFit: Fast spectral synthesis and emission-line fitting of DESI spectra, Astrophysics Source Code Library, record ascl:2308.005 |
| Netzer, H. 2015, ARA&A, 53, 365,                                                                                                                                                                    |

doi: 10.1146/annurev-astro-082214-122302

| Noda, H., &Done, C. 2018, MNRAS, 480, 3898, doi: 10.1093/mnras/sty2032                                                                                                    |
|---------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| Padovani, P., Alexander, D. M., Assef, R. J., et al. 2017, A&A Rv, 25, 2, doi: 10.1007/s00159-017-0102-9                                                                  |
| Penston, M. V., &Perez, E. 1984, MNRAS, 211, 33P, doi: 10.1093/mnras/211.1.33P                                                                                            |
| Planck Collaboration, Aghanim, N., Akrami, Y., et al. 2020, A&A, 641, A6, doi: 10.1051/0004-6361/201833910                                                                |
| Raichoor, A., Moustakas, J., Newman, J. A., et al. 2023, AJ, 165, 126, doi: 10.3847/1538-3881/acb213                                                                      |
| Raimundo, S. I., Vestergaard, M., Koay, J. Y., et al. 2019, MNRAS, 486, 123, doi: 10.1093/mnras/stz852                                                                    |
| Rakshit, S., Stalin, C. S., &Kotilainen, J. 2020, ApJS, 249, 17, doi: 10.3847/1538-4365/ab99c5                                                                            |
| Rees, M. J. 1984, ARA&A, 22, 471, doi: 10.1146/annurev.aa.22.090184.002351                                                                                                |
| Ricci, C., &Trakhtenbrot, B. 2023, Nature Astronomy, 7, 1282, doi: 10.1038/s41550-023-02108-4                                                                             |
| Ruan, J. J., Anderson, S. F., Eracleous, M., et al. 2019, ApJ, 883, 76, doi: 10.3847/1538-4357/ab3c1a                                                                     |
| Ruan, J. J., Anderson, S. F., Cales, S. L., et al. 2016, ApJ, 826, 188, doi: 10.3847/0004-637X/826/2/188                                                                  |
| Runnoe, J. C., Cales, S., Ruan, J. J., et al. 2016, MNRAS, 455, 1691, doi: 10.1093/mnras/stv2385                                                                          |
| S´nchez-S´ez, a a P., Lira, P., Mej´a-Restrepo, ı J., et al. 2018, ApJ, 864, 87, doi: 10.3847/1538-4357/aad7f9                                                            |
| Schlafly, E. F., Kirkby, D., Schlegel, D. J., et al. 2023, arXiv e-prints, arXiv:2306.06309, doi: 10.48550/arXiv.2306.06309                                               |
| Schlegel et al. 2023, in preparation                                                                                                                                      |
| Shakura, N. I., &Sunyaev, R. A. 1973, A&A, 500, 33 Shapovalova, A. I., Popovi´, c L. ˇ C., Burenkov, A. N., et al. 2010, A&A, 509, A106, doi: 10.1051/0004-6361/200912311 |
| Shappee, B. J., Prieto, J. L., Grupe, D., et al. 2014, ApJ, 788, 48, doi: 10.1088/0004-637X/788/1/48                                                                      |
| Shen, Y., Greene, J. E., Strauss, M. A., Richards, G. T., & Schneider, D. P. 2008, ApJ, 680, 169, doi: 10.1086/587475                                                     |
| Sheng, Z., Wang, T., Jiang, N., et al. 2017, ApJL, 846, L7, doi: 10.3847/2041-8213/aa85de                                                                                 |
| Silber, J. H., Fagrelius, P., Fanning, K., et al. 2023, AJ, 165, 9, doi: 10.3847/1538-3881/ac9ab1                                                                         |
| Smee, S. A., Gunn, J. E., Uomoto, A., et al. 2013, AJ, 146, 32, doi: 10.1088/0004-6256/146/2/32                                                                           |
| Stern, D., McKernan, B., Graham, M. J., et al. 2018, ApJ, 864, 27, doi: 10.3847/1538-4357/aac726                                                                          |
| Stern, J., &Laor, A. 2012, MNRAS, 426, 2703, doi: 10.1111/j.1365-2966.2012.21772.x                                                                                        |
| Stone, Z., Shen, Y., Burke, C. J., et al. 2022, MNRAS, 514, 164, doi: 10.1093/mnras/stac1259                                                                              |

| Sun, M., Trump, J. R., Shen, Y., et al. 2015, ApJ, 811, 42, doi: 10.1088/0004-637X/811/1/42                           |
|-----------------------------------------------------------------------------------------------------------------------|
| Sun, M., Xue, Y., Guo, H., et al. 2020, ApJ, 902, 7, doi: 10.3847/1538-4357/abb1c4                                    |
| Trump, J. R., Impey, C. D., Kelly, B. C., et al. 2011, ApJ, 733, 60, doi: 10.1088/0004-637X/733/1/60                  |
| Vaughan, S., Edelson, R., Warwick, R. S., &Uttley, P. 2003, MNRAS, 345, 1271, doi: 10.1046/j.1365-2966.2003.07042.x   |
| Wang, J., Xu, D. W., Cao, X., et al. 2024a, arXiv e-prints, arXiv:2405.10663, doi: 10.48550/arXiv.2405.10663          |
| Wang, J., Xu, D. W., Wang, Y., et al. 2019, ApJ, 887, 15, doi: 10.3847/1538-4357/ab4d90                               |
| Wang, J., Zheng, W. K., Brink, T. G., et al. 2023, ApJ, 956, 137, doi: 10.3847/1538-4357/acf5e0                       |
| Wang, J.-M., &Bon, E. 2020, A&A, 643, L9, doi: 10.1051/0004-6361/202039368                                            |
| Wang, J.-M., Liu, J.-R., Ho, L. C., Li, Y.-R., &Du, P. 2021, ApJL, 916, L17, doi: 10.3847/2041-8213/ac0b46            |
| Wang, S., Woo, J.-H., Gallo, E., et al. 2024b, arXiv e-prints, arXiv:2402.18131, doi: 10.48550/arXiv.2402.18131       |
| Wolf, C., Golding, J., Hon, W. J., &Onken, C. A. 2020, MNRAS, 499, 1005, doi: 10.1093/mnras/staa2794                  |
| Woo, J.-H., Yoon, Y., Park, S., Park, D., &Kim, S. C. 2015, ApJ, 801, 38, doi: 10.1088/0004-637X/801/1/38             |
| Wu, W.-B., &Gu, W.-M. 2023, ApJ, 958, 146, doi: 10.3847/1538-4357/acf839                                              |
| Xu, X., &Wang, J. 2022, ApJ, 933, 110, doi: 10.3847/1538-4357/ac7222                                                  |
| Yan, L., Wang, T., Jiang, N., et al. 2019, ApJ, 874, 44, doi: 10.3847/1538-4357/ab074b                                |
| Yang, Q., Wu, X.-B., Fan, X., et al. 2018, ApJ, 862, 109, doi: 10.3847/1538-4357/aaca3a                               |
| Yang, Q., Green, P. J., MacLeod, C. L., et al. 2023, arXiv e-prints, arXiv:2303.06733, doi: 10.48550/arXiv.2303.06733 |
| York, D. G., Adelman, J., Anderson, John E., J., et al. 2000, AJ, 120, 1579, doi: 10.1086/301513                      |
| Yu, X., Shi, Y., Chen, Y., et al. 2020, MNRAS, 498, 3985, doi: 10.1093/mnras/staa2627                                 |
| Yu, Z., Martini, P., Penton, A., et al. 2023, MNRAS, 522, 4132, doi: 10.1093/mnras/stad1224                           |
| Zeltyn, G., Trakhtenbrot, B., Eracleous, M., et al. 2022, ApJL, 939, L16, doi: 10.3847/2041-8213/ac9a47               |
| Zhou, R., Dey, B., Newman, J. A., et al. 2022, arXiv e-prints, arXiv:2208.08515. https://arxiv.org/abs/2208.08515     |
| Zuo, W., Wu, X.-B., Liu, Y.-Q., &Jiao, C.-L. 2012, ApJ, 758, 104, doi: 10.1088/0004-637X/758/2/104                    |
Wei-Jian Guo Hu Zou Claire L. Greenwell David M. Alexander Victoria A. Fawcett Zhiwei Pan Malgorzata Siudek Jessica Nicole Aguilar Steven Ahlen David Brooks Todd Claybaugh Kyle Dawson Axel De La Macorra Peter Doel Andreu Font-Ribera Enrique Gaztanaga Satya Gontcho A Gontcho Gaston Gutierrez Robert Kehoe Theodore Kisner Martin Landriau Laurent Le Guillou Marc Manera Aaron Meisner Ramon Mique John Moustakas Francisco Prada Graziano Rossi Eusebio Sanchez Michael Schubnell David Sprayberry Jipeng Sui Gregory Tarle Benjamin Alan Weaver Yun-Ao Xiao Siwei Zou
Changing-look Active Galactic Nuclei from the Dark Energy Spectroscopic Instrument. II. Statistical Properties from the First Data Release
We present the identification of changing-look active galactic nuclei
(CL-AGNs) from the Dark Energy Spectroscopic Instrument First Data Release and
Sloan Digital Sky Survey Data Release 16 at z \leq 0.9. To confirm the CL-AGNs,
we utilize spectral flux calibration assessment via an [O\,{\sc iii}]-based
calibration, pseudo-photometry examination, and visual inspection. This
rigorous selection process allows us to compile a statistical catalog of 561
CL-AGNs, encompassing 527 $\rm H\beta$, 149$\rm H\alpha$, and 129 Mg II CL
behaviors. In this sample, we find 1) a 283:278 ratio of turn-on to turn-off
CL-AGNs. 2) the critical value for CL events is confirmed around Eddington
ratio \sim 0.01. 3) a strong correlation between the change in the luminosity
of the broad emission lines (BEL) and variation in the continuum luminosity,
with Mg II and $\rm H\beta$ displaying similar responses during CL phases. 4)
the Baldwin-Phillips-Terlevich diagram for CL-AGNs shows no statistically
difference from the general AGN catalog. 5) five CL-AGNs are associated with
asymmetrical mid-infrared flares, possibly linked to tidal disruption events.
Given the large CL-AGNs and the stochastic sampling of spectra, we propose that
some CL events are inherently due to typical AGN variability during low
accretion rates, particularly for CL events of the singular BEL. Finally, we
introduce a Peculiar CL phase, characterized by a gradual decline over decades
in the light curve and the complete disappearance of entire BEL in faint
spectra, indicative of a real transition in the accretion disk.