Search

The outshining light from active galactic nuclei (AGNs) poses significant challenges in studying the properties of AGN host galaxies. To address this issue, we propose a novel approach which combines image decomposition and spectral energy distribution (SED) decomposition to constrain properties of AGN host galaxies. Image decomposition allows us to disentangle optical flux into AGN and stellar components, thereby providing additional constraints on the SED models to derive more refined stellar mass. To test the viability of this approach, we obtained a sample of 24 X-ray selected type-I AGNs with redshifts ranging from 0.73 to 2.47. We estimated the stellar masses for our sample and found that our results are generally consistent with earlier estimates based on different methods. Through examining the posterior distribution of stellar masses, we find that our method could derive better constrained results compared to previous SED decomposition methods. With the derived stellar masses, we further studied the $M_{\rm BH}-M_\star$ relation of our sample, finding a higher intrinsic scatter in the correlation for our entire sample compared to the local quiescent correlation, which could be caused by a few black hole monsters in our sample. We propose that based on our method, future works could extend to larger samples of high-redshift AGN host galaxies, thereby enhancing our understanding of their properties., Comment: 18 pages, 7 figures, submitted to AAS journal

Tidal disruption events (TDEs) could be an important growth channel for massive black holes in dwarf galaxies. Theoretical work suggests that the observed active galactic nuclei (AGNs) in dwarf galaxies are predominantly TDE-powered. To assess this claim, we perform variability analyses on the dwarf-hosted AGNs detected in the $7$ Ms Chandra Deep Field-South (CDF-S) survey, with observations spanning $\approx 16$ years. Based on the spectral energy distribution (SED) modeling with X-CIGALE, we select AGNs hosted by dwarf galaxies (stellar mass below $10^{10}\ M_\odot$). We focus on X-ray sources with full-band detections, leading to a sample of $78$ AGNs (0.122 $\leq$ $z$ $\leq$ 3.515). We fit the X-ray light curves with a canonical TDE model of $t^{-5/3}$ and a constant model. If the former outperforms the latter in fitting quality for a source, we consider the source as a potential TDE. We identify five potential TDEs, constituting a small fraction of our sample. Using true- and false-positive rates obtained from fitting models to simulated light curves, we perform Bayesian analysis to obtain the posterior of the TDE fraction for our sample. The posterior peaks close to zero ($2.56\%$), and we obtain a $2$-$\sigma$ upper limit of $9.80\%$. Therefore, our result indicates that the observed AGNs in dwarf galaxies are not predominantly powered by TDEs., Comment: 16 pages, 10 figures, accepted for publication in ApJ

Thanks to the rapidly increasing time-domain facilities, we are entering a golden era of research on gamma-ray bursts (GRBs). In this Letter, we report our observations of GRB 240529A with the Burst Optical Observer and Transient Exploring System, the 1.5-meter telescope at Observatorio Sierra Nevada, the 2.5-meter Wide Field Survey Telescope of China, the Large Binocular Telescope, and the Telescopio Nazionale Galileo. The prompt emission of GRB 240529A shows two comparable energetic episodes separated by a quiescence time of roughly 400 s. Combining all available data on the GRB Coordinates Network, we reveal the simultaneous apparent X-ray plateau and optical re-brightening around $10^3-10^4$ s after the burst. Rather than the energy injection from the magnetar as widely invoked for similar GRBs, the multi-wavelength emissions could be better explained as two shocks launched from the central engine separately. The optical peak time and our numerical modeling suggest that the initial bulk Lorentz factor of the later shock is roughly 50, which indicates that the later jet should be accretion-driven and have a higher mass loading than a typical one. The quiescence time between the two prompt emission episodes may be caused by the transition between different accretion states of a central magnetar or black hole, or the fall-back accretion process. A sample of similar bursts with multiple emission episodes in the prompt phase and sufficient follow-up could help to probe the underlying physics of GRB central engines., Comment: Resubmitted to ApJL after addressing the referee's comments; comments are welcome

Supermassive black holes (SMBHs) can grow through both accretion and mergers. It is still unclear how SMBHs evolve under these two channels from high redshifts to the SMBH population we observe in the local universe. Observations can directly constrain the accretion channel but cannot effectively constrain mergers yet, while cosmological simulations provide galaxy merger information but can hardly return accretion properties consistent with observations. In this work, we combine the observed accretion channel and the simulated merger channel, taking advantage of both observations and cosmological simulations, to depict a realistic evolution pattern of the SMBH population. With this methodology, we can derive the scaling relation between the black-hole mass ($M_\mathrm{BH}$) and host-galaxy stellar mass ($M_\star$) and the local black-hole mass function (BHMF). Our scaling relation is lower than those based on dynamically measured $M_\mathrm{BH}$, supporting the claim that dynamically measured SMBH samples may be biased. We show that the scaling relation has little redshift evolution. The BHMF steadily increases from $z=4$ to $z=1$ and remains largely unchanged from $z=1$ to $z=0$. The overall SMBH growth is generally dominated by the accretion channel, with possible exceptions at high mass ($M_\mathrm{BH}\gtrsim10^{8}~M_\odot$ or $M_\star\gtrsim10^{11}~M_\odot$) and low redshift ($z\lesssim1$). We also predict that around 25% of the total SMBH mass budget in the local universe may be locked within long-lived, wandering SMBHs, and the wandering mass fraction and wandering SMBH counts increase with $M_\star$., Comment: 29 pages, 22 figures, 5 tables, accepted for publication in ApJ

The characteristic timescale at which the variability of active galactic nuclei (AGNs) turns from red noise to white noise can probe the accretion physics around supermassive black holes (SMBHs). A number of works have studied the characteristic timescale of quasars and obtained quite different scaling relations between the timescale and quasar physical properties. One possible reason for the discrepancies is that the characteristic timescale can be easily underestimated if the light curves are not long enough. In this work, we construct well-defined AGN samples to observationally test the relationships between the characteristic timescale and AGN properties obtained by previous works. Our samples eliminate the effects of insufficient light-curve lengths. We confirm that the timescale predictions \citep{Zhou2024} of the Corona Heated Accretion disk Reprocessing model are consistent with our timescale measurements. The timescale predictions by empirically relations \citep[e.g.,][]{Kelly2009} are systematically smaller than our measured ones. Our results provide further evidence that AGN variability is driven by thermal fluctuations in SMBH accretion disks. Future flagship time-domain surveys can critically test our conclusions and reveal the physical nature of AGN variability., Comment: 20 pages, 7 figures, Accepted to ApJ

Resolving the inner structures of active galactic nuclei (AGNs) provides the "standard ruler" to measure the parallax distances of the Universe and a powerful way to weigh supermassive black holes (SMBHs). Thanks to time-domain observations, it is possible to use the reverberation mapping (RM) technique to measure time delays between different light curves that probe the structures of the SMBH accretion disks and broad line regions (BLRs), which are otherwise often too compact to be spatially resolved for most AGNs. Despite decades of RM studies, the critical physical process that controls the structures of SMBH accretion disk and BLR and their temporal evolution remains unclear. Here we report the variation of the SMBH accretion disk structure of NGC 4151 in response to changes in luminosity within 6 years. In the high-flux state, the time delays measured from our continuum RM with high-cadence (2 days) spectroscopy are 3.8 times larger than that in the low-flux state and 15 times longer than the classical standard thin disk (SSD) prediction. This result provides the first piece of direct evidence that the SMBH disk structure "breathes" in highly-variable AGN manifestations. The time-delay change severely challenges the popular X-ray reprocessing of the SSD model, with or without BLR contributions. More importantly, the continuum time delays can be comparable with the time delay between the broad Hb line and the nearby optical continuum, and the latter is commonly used to calculate the BLR sizes. Hence, the BLR sizes are significantly underestimated if the continuum time delays are not properly considered. This underestimation introduces up to 0.3 dex systematic uncertainties on RM SMBH masses and BLR parallax distances. Our findings underscore that simultaneous continuum and BLR RM studies are vital for better deciphering the SMBH mass growth and the cosmological expansion history., Comment: 18 pages, 8 figures, comments welcome

The emission of Anomalous X-ray Pulsars (AXPs) and Soft Gamma-Ray Repeaters (SGRs) is believed to be powered by the dissipation of their strong magnetic fields, which coined the name `magnetar'. By combining timing and energy observational results, the magnetar model can be easily appreciated. From a timing perspective, the magnetic field strengths of AXPs and SGRs, calculated assuming dipole radiation, are extremely strong. From an energy perspective, the X-ray/soft gamma-ray luminosities of AXPs and SGRs are larger than their rotational energy loss rates (i.e., $ L_{\rm X}>\dot E_{\rm rot}$). It is thus reasonable to assume that the high-energy radiation comes from magnetic energy decay, and the magnetar model has been extensively discussed (or accepted). However, we argue that: (1) calculating magnetic fields by assuming that rotational energy loss is dominated by dipole radiation (i.e., $\dot{E}_{\rm rot}\simeq\dot{E}_{\mu}$)) may be controversial, and we suggest that the energies carried by outflowing particles should also be considered; and (2) the fact that X-ray luminosity is greater than the rotational energy loss rate does not necessarily mean that the emission energy comes from the magnetic field decaying, which requires further observational testing. Furthermore, some observational facts conflict with the `magnetar' model, such as observations of anti-magnetars, high magnetic field pulsars, and radio and X-ray observations of AXPs/SGRs. Therefore, we propose a crusted strange star model as an alternative, which can explain many more observational facts of AXPs/SGRs., Comment: 12 pages, 3 figures, accepted for publication in JUSTC

We introduce our project, AGNSTRONG (Active Galactic Nuclei and STaR fOrmation in Nearby Galaxies). Our research goals encompass investigating the kinematic properties of ionized and molecular gas outflows, understanding the impact of AGN feedback, and exploring the coevolution dynamics between AGN strength activity and star formation activity. We aim to conduct a thorough analysis to determine whether there is an increase or suppression in SFRs among targets with and without powerful relativistic jets. Our sample consists of 35 nearby AGNs with and without powerful relativistic jet detections. Utilizing sub-millimeter (sub-mm) continuum observations at 450 {\mu}m and 850 {\mu}m from SCUBA-2 at the James Clerk Maxwell Telescope, we determine star-formation rates (SFRs) for our sources using spectral energy distribution (SED) fitting models. Additionally, we employ high-quality, spatially resolved spectra from UV-optical to near-infrared bands obtained with the Double Spectrograph and Triple Spectrograph mounted on the 200-inch Hale telescope at Palomar Observatory to study their multiphase gas outflow properties. This paper presents an overview of our sample selection methodology, research strategy, and initial results of our project. We find that the SFRs determined without including the sub-mm data in the SED fitting are overestimated by approximately 0.08 dex compared to those estimated with the inclusion of sub-mm data. Additionally, we compare the estimated SFRs in our work with those traced by the 4000{\AA} break, as provided by the MPA-JHU catalog. We find that our determined SFRs are systematically higher than those traced by the 4000{\AA} break. Finally, we outline our future research plans., Comment: 13 pages, accepted for publication in The Astronomical Journal

Stellar black holes (sBHs) are widely believed to exist in the accretion disks of active galactic nuclei (AGNs). Previous studies often focus on the transient emission produced by embedded sBHs. Here, we explore the possible observational consequences of an AGN accretion disk that contains a population of accreting sBHs. Embedded accreting sBHs change the effective temperature distribution of the AGN accretion disk by heating gas in the outer regions. Two possible observational consequences are presented. First, the spectral energy distribution has a turnover feature at $\sim 4700\ \textrm{\AA}$ when the supermassive black hole (SMBH) mass is $\sim 10^8\ M_{\odot}$, which can help explain the observed shallow spectral shape at wavelengths $>5000\ \textrm{\AA}$ for the Sloan Digital Sky Survey quasar composite spectrum. Second, the half-light radius of a given relatively long wavelength is significantly larger than for an AGN disk without sBHs, which can be tested by microlensing observations. With appropriate sBH distributions, the model can be reconciled with quasar microlensing disk sizes. We propose that the half-light radius-wavelength relation can be utilized to investigate the distributions of embedded sBHs in AGN accretion disks., Comment: 11 pages, 5 figures, Accepted to ApJL

The UV/optical light curves of Active Galactic Nuclei (AGNs) are commonly described by the Damped Random Walk (DRW) model. However, the physical interpretation of the damping timescale, a key parameter in the DRW model, remains unclear. Particularly, recent observations indicate a weak dependence of the damping timescale upon both wavelength and accretion rate, clearly being inconsistent with the accretion-disk theory. In this study, we investigate the damping timescale in the framework of the Corona Heated Accretion disk Reprocessing (CHAR) model, a physical model that describes AGN variability. We find that while the CHAR model can reproduce the observed power spectral densities of the 20-year light curves for 190 sources from \cite{Stone2022}, the observed damping timescale, as well as its weak dependence on wavelength, can also be well recovered through fitting the mock light curves with DRW. We further demonstrate that such weak dependence is artificial due to the effect of inadequate durations of light curves, which leads to best-fitting damping timescales lower than the intrinsic ones. After eliminating this effect, the CHAR model indeed yields a strong dependence of the intrinsic damping timescale on the bolometric luminosity and rest-frame wavelength. Our results highlight the demand for sufficiently long light curves in AGN variability studies and important applications of the CHAR model in such studies., Comment: 19 pages, 16 figures, accepted to ApJ

High-cadence, multiwavelength observations have continuously revealed the diversity of tidal disruption events (TDEs), thus greatly advancing our knowledge and understanding of TDEs. In this work, we conducted an intensive optical-UV and X-ray follow-up campaign of TDE AT2023lli, and found a remarkable month-long bump in its UV/optical light curve nearly two months prior to maximum brightness. The bump represents the longest separation time from the main peak among known TDEs to date. The main UV/optical outburst declines as $t^{-4.10}$, making it one of the fastest decaying optically selected TDEs. Furthermore, we detected sporadic X-ray emission 30 days after the UV/optical peak, accompanied by a reduction in the period of inactivity. It is proposed that the UV/optical bump could be caused by the self-intersection of the stream debris, whereas the primary peak is generated by the reprocessed emission of the accretion process. In addition, our results suggest that episodic X-ray radiation during the initial phase of decline may be due to the patched obscurer surrounding the accretion disk, a phenomenon associated with the inhomogeneous reprocessing process. The double TDE scenario, in which two stars are disrupted in sequence, is also a possible explanation for producing the observed early bump and main peak. We anticipate that the multicolor light curves of TDEs, especially in the very early stages, and the underlying physics can be better understood in the near future with the assistance of dedicated surveys such as the deep high-cadence survey of the 2.5-meter Wide Field Survey Telescope (WFST)., Comment: 14 pages, 8 figures,accepted for publication by ApJL

Highly radio-loud quasars (HRLQs; $\log R>2.5$) at $z\gtrsim 4$ show apparent enhanced X-ray emission compared to matched HRLQs at lower redshifts, perhaps due to a redshift-dependent fractional contribution to the X-ray luminosity from inverse-Compton scattering of cosmic microwave background photons (IC/CMB). Using new {\it Chandra} observations and archival X-ray data, we investigate this phenomenon with an optically flux-limited sample of 41 HRLQs at $z = 4$--5.5 all with sensitive X-ray coverage, the largest sample utilized to date by a wide margin. X-ray enhancements are assessed using X-ray-to-optical flux ratios and spectral energy distributions. We confirm the presence of X-ray enhancements at a 4.9--5.3$\sigma$ significance level, finding that the median factor of enhancement is $\approx 1.8$ at our sample median redshift of $z\approx 4.4$. Under a fractional IC/CMB model, the expected enhancement at lower redshifts is modest; e.g., $\approx 4$% at $z\approx 1.5$. We also investigate a sample of seven radio-loud quasars (RLQs; $\log R>1$) at even higher redshifts of $z=5.6$--6.8, using new and archival X-ray data. These RLQs also show evidence for X-ray enhancements by a median factor of $\approx 2.7$ at a 3.7--4.9$\sigma$ significance level. The X-ray spectral and other properties of these $z=5.6$--6.8 RLQs, however, pose challenges for a straightforward fractional IC/CMB interpretation of their enhancements., Comment: 26 pages, 12 figures, accepted for publication in MNRAS

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

Quasi-periodic oscillations (QPOs) have been widely observed in black hole X-ray binaries (BHBs), which often exhibit significant X-ray variations. Extensive research has explored the long-term evolution of the properties of QPOs in BHBs. In contrast, such evolution in active galactic nuclei (AGNs) has remained largely unexplored due to limited observational data. By using the 10 new XMM-Newton observations for the narrow-line Seyfert 1 galaxy RE J1034+396 from publicly available data, we analyze the characteristics of its X-ray QPOs and examine their long-term evolution. The hard-band (1--4 keV) QPOs are found in all 10 observations and the frequency of these QPOs evolves ranging at $(2.47\text{--}2.83)\times10^{-4}\rm\ Hz$. Furthermore, QPO signals in the soft (0.3--1 keV) and hard bands exhibit strong coherence, although, at times, the variations in the soft band lead those in the hard band (the hard-lag mode), while at other times, it is the reverse (the soft-lag mode). The observations presented here serendipitously captured two ongoing lag reversals within about two weeks, which are first seen in RE J1034+396 and also among all AGNs. A transition in QPO frequency also takes place within a two-week timeframe, two weeks prior to its corresponding lag reversal, indicating a possible coherence between the transitions of QPO frequency and lag mode with delay. The diagram of time lag versus QPO frequency clearly evidences this interconnected evolution with hysteresis, which is, for the first time, observed among AGNs., Comment: 11 pages, 4 figures, accepted for publication in ApJ Letters

"Changing-look" active galactic nuclei (CL-AGNs) challenge our basic ideas about the physics of accretion flows and circumnuclear gas around supermassive black holes. Using first-year Sloan Digital Sky Survey V (SDSS-V) repeated spectroscopy of nearly 29,000 previously known AGNs, combined with dedicated follow-up spectroscopy, and publicly available optical light curves, we have identified 116 CL-AGNs where (at least) one broad emission line has essentially (dis-)appeared, as well as 88 other extremely variable systems. Our CL-AGN sample, with 107 newly identified cases, is the largest reported to date, and includes $\sim0.4\%$ of the AGNs reobserved in first-year SDSS-V operations. Among our CL-AGNs, 67% exhibit dimming while 33% exhibit brightening. Our sample probes extreme AGN spectral variability on months to decades timescales, including some cases of recurring transitions on surprisingly short timescales ($\lesssim 2$ months in the rest frame). We find that CL events are preferentially found in lower-Eddington-ratio ($f_{Edd}$) systems: Our CL-AGNs have a $f_{Edd}$ distribution that significantly differs from that of a carefully constructed, redshift- and luminosity-matched control sample (Anderson-Darling test yielding $p_{\rm AD}\approx 6\times10^{-5}$; median $f_{Edd}\approx0.025$ vs. $0.043$). This preference for low $f_{Edd}$ strengthens previous findings of higher CL-AGN incidence at lower $f_{Edd}$, found in smaller samples. Finally, we show that the broad MgII emission line in our CL-AGN sample tends to vary significantly less than the broad H$\beta$ emission line. Our large CL-AGN sample demonstrates the advantages and challenges in using multi-epoch spectroscopy from large surveys to study extreme AGN variability and physics., Comment: Submitted to ApJ. Full tables and figure-sets will be published upon acceptance, and can be made available upon request$.$