Search

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

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

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

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$.$

The Hot Universe Baryon Surveyor (HUBS) is a proposed space-based X-ray telescope for detecting X-ray emissions from the hot gas content in our universe. With its unprecedented spatially-resolved high-resolution spectroscopy and large field of view, the HUBS mission will be uniquely qualified to measure the physical and chemical properties of the hot gas in the interstellar medium, the circumgalactic medium, the intergalactic medium, and the intracluster medium. These measurements will be valuable for two key scientific goals of HUBS, namely to unravel the AGN and stellar feedback physics that governs the formation and evolution of galaxies, and to probe the baryon budget and multi-phase states from galactic to cosmological scales. In addition to these two goals, the HUBS mission will also help us solve some problems in the fields of galaxy clusters, AGNs, diffuse X-ray backgrounds, supernova remnants, and compact objects. This paper discusses the perspective of advancing these fields using the HUBS telescope., Comment: 52 pages, 22 figures. Accepted for publication in Science China: Physics, Mechanics and Astronomy

The Wide Field Survey Telescope (WFST) is a dedicated photometric surveying facility being built jointly by the University of Science and Technology of China and the Purple Mountain Observatory. It is equipped with a 2.5-meter diameter primary mirror, an active optics system, and a mosaic CCD camera with 0.73 gigapixels on the primary focal plane for high-quality image capture over an FOV of 6.5-square-degree. It is anticipated that WFST will be set up at the Lenghu site in the summer of 2023 and begin to observe the northern sky in four optical bands (u, g, r, and i) with a range of cadences, from hourly/daily in the Deep High-Cadence Survey (DHS) program to semiweekly in the Wide-Field Survey (WFS) program, three months later. During a photometric night, a nominal 30 s exposure in the WFS program will reach a depth of 22.27, 23.32, 22.84, and 22.31 (AB magnitudes) in these four bands, respectively, allowing for the detection of a tremendous amount of transients in the low-z universe and a systematic investigation of the variability of Galactic and extragalactic objects. In the DHS program, intranight 90 s exposures as deep as 23 (u) and 24 mag (g), in combination with target of opportunity follow-ups, will provide a unique opportunity to explore energetic transients in demand for high sensitivities, including the electromagnetic counterparts of gravitational wave events, supernovae within a few hours of their explosions, tidal disruption events and fast, luminous optical transients even beyond a redshift of unity. In addition, the final 6-year co-added images, anticipated to reach g=25.8 mag in WFS or 1.5 mags deeper in DHS, will be of fundamental importance to general Galactic and extragalactic science. The highly uniform legacy surveys of WFST will serve as an indispensable complement to those of LSST that monitor the southern sky., Comment: 48 pages

The Australia Telescope Large Area Survey (ATLAS) and the VLA survey in the XMM-LSS/VIDEO deep field provide deep ($\approx 15$ ${\mu}$Jybeam$^{-1}$) and high-resolution ($\approx$ 4.5--8 arcsec) radio coverage of the three XMM-SERVS fields (W-CDF-S, ELAIS-S1, and XMM-LSS). These data cover a total sky area of 11.3 deg$^2$ and contain $\approx 11000$ radio components. Furthermore, about 3~deg$^2$ of the XMM-LSS field also has deeper MIGHTEE data that achieve a median RMS of 5.6 ${\mu}$Jy beam$^{-1}$ and detect more than 20000 radio sources. We analyze all these radio data and find source counterparts at other wavebands utilizing deep optical and IR surveys. The nature of these radio sources is studied using radio-band properties (spectral slope and morphology), and the IR-radio correlation. %and spectral energy distribution. Radio AGNs are selected and compared with those selected using other methods (e.g. X-ray). We found 1656 new AGNs that were not selected using X-ray and/or MIR methods. We constrain the FIR-to-UV SEDs of radio AGNs using {\sc cigale} and investigate the dependence of radio AGN fraction upon galaxy stellar mass and star-formation rate., Comment: 25 pages, 21+3 figures, 8+1 tables, MNRAS accepted (final version). The catalogues are available from https://personal.psu.edu/wnb3/xmmservs/zhu2023-tables-2023-04-12.tar.zip

We present a routinized and reliable method to obtain source catalogs from the $\textit{Nuclear Spectroscopic Telescope Array}$ ($\textit{NuSTAR}$) extragalactic surveys of the Extended $\textit{Chandra}$ Deep Field-South (E-CDF-S) and $\textit{Chandra}$ Deep Field-North (CDF-N). The $\textit{NuSTAR}$ E-CDF-S survey covers a sky area of $\approx30'\times30'$ to a maximum depth of $\sim$ 230 ks corrected for vignetting in the 3--24 keV band, with a total of 58 sources detected in our E-CDF-S catalog; the $\textit{NuSTAR}$ CDF-N survey covers a sky area of $\approx7'\times10'$ to a maximum depth of $\sim$ 440 ks corrected for vignetting in the 3--24 keV band, with a total of 42 sources detected in our CDF-N catalog that is produced for the first time. We verify the reliability of our two catalogs by crossmatching them with the relevant catalogs from the $\textit{Chandra}$ X-ray observatory, and find that the fluxes of our $\textit{NuSTAR}$ sources are generally consistent with that of their $\textit{Chandra}$ counterparts. Our two catalogs are produced following the exactly same method and made publicly available, thereby providing a uniform platform that facilitates further studies involving these two fields. Our source-detection method provides a systematic approach for source cataloging in other $\textit{NuSTAR}$ extragalactic surveys., Comment: 26 pages, 7 figures, accepted for publication in JUSTC; click "Full text" & "Supplements" for catalogs in https://doi.org/10.52396/JUSTC-2023-0032

We investigate the extreme X-ray variability of a z = 1.608 active galactic nucleus in the 7 Ms Chandra Deep Field-South (XID 403), which showed two significant X-ray brightening events. In the first event, XID 403 brightened by a factor of $>2.5$ in $\lesssim6.1$ rest-frame days in the observed-frame 0.5-5 keV band. The event lasted for $\approx5.0\textrm{-}7.3$ days, and then XID 403 dimmed by a factor of $>6.0$ in $\lesssim6.1$ days. After $\approx1.1\textrm{-}2.5$ years in the rest frame (including long observational gaps), it brightened again with the 0.5-5 keV flux increasing by a factor of $>12.6$. The second event lasted over 251 days and the source remained bright until the end of the 7 Ms exposure. The spectrum is a steep power law (photon index $\Gamma=2.8\pm0.3$) without obscuration during the second outburst, and the rest-frame 2-10 keV luminosity reaches $1.5^{+0.8}_{-0.5}\times10^{43}$ erg s$^{-1}$; there is no significant spectral evolution within this epoch. The infrared-to-UV spectral energy distribution of XID 403 is dominated by the host galaxy. There is no significant optical/UV variability and $R$-band (rest-frame $\approx2500$ $\unicode{xC5}$) brightening contemporaneous with the X-ray brightening. The extreme X-ray variability is likely due to two X-ray unveiling events, where the line of sight to the corona is no longer shielded by high-density gas clumps in a small-scale dust-free absorber. XID 403 is probably a high-redshift analog of local narrow-line Seyfert 1 galaxies, and the X-ray absorber is a powerful accretion-disk wind. On the other hand, we cannot exclude the possibility that XID 403 is an unusual candidate for tidal disruption events., Comment: 23 pages, 9 figures, accepted for publication in ApJ

The shape of the ambient circumnuclear medium (ACM) density profile can probe the history of accretion onto the central supermassive black hole in galaxies and the circumnuclear environment. However, due to the limitation of the instrument resolution, the density profiles of the ACM for most of galaxies remain largely unknown. In this work, we propose a novel method to measure the ACM density profile of active galactic nucleus (AGN) by the equilibrium between the radiation pressure on the warm absorbers (WAs, a type of AGN outflows) and the drag pressure from the ACM. We study the correlation between the outflow velocity and ionization parameter of WAs in each of the five Seyfert 1 galaxies (NGC 3227, NGC 3783, NGC 4051, NGC 4593, and NGC 5548), inferring that the density profile of the ACM is between n\propto r^-1.7 and n \propto r^-2.15 (n is number density and r is distance) from 0.01 pc to pc scales in these five AGNs. Our results indicate that the ACM density profile in Seyfert 1 galaxies is steeper than the prediction by the spherically symmetric Bondi accretion model and the simulated results of the hot accretion flow, but more in line with the prediction by the standard thin disk model., Comment: 8 pages, 3 figures, 2 tables, accepted for publication in The Astrophysical Journal

The Chandra Deep Field-South and North surveys (CDFs) provide unique windows into the cosmic history of X-ray emission from normal (non-active) galaxies. Scaling relations of normal galaxy X-ray luminosity (L_X) with star formation rate (SFR) and stellar mass (M_star) have been used to show that the formation rates of low-mass and high-mass X-ray binaries (LMXBs and HMXBs, respectively) evolve with redshift across z = 0-2 following L_HMXB/SFR ~ 1 + z and L_LMXB/M_star ~ (1 + z)^{2-3}. However, these measurements alone do not directly reveal the physical mechanisms behind the redshift evolution of X-ray binaries (XRBs). We derive star-formation histories for a sample of 344 normal galaxies in the CDFs, using spectral energy distribution (SED) fitting of FUV-to-FIR photometric data, and construct a self-consistent, age-dependent model of the X-ray emission from the galaxies. Our model quantifies how X-ray emission from hot gas and XRB populations vary as functions of host stellar-population age. We find that (1) the ratio L_X/M_star declines by a factor of ~1000 from 0-10 Gyr and (2) the X-ray SED becomes harder with increasing age, consistent with a scenario in which the hot gas contribution to the X-ray SED declines quickly for ages above 10 Myr. When dividing our sample into subsets based on metallicity, we find some indication that L_X/M_star is elevated for low-metallicity galaxies, consistent with recent studies of X-ray scaling relations. However, additional statistical constraints are required to quantify both the age and metallicity dependence of X-ray emission from star-forming galaxies.

The relation between warm absorber (WA) outflows of AGN and nuclear obscuration activities caused by optically-thick clouds (obscurers) crossing the line of sight is unclear. NGC 3227 is a suitable target to study the properties of both WAs and obscurers, because it matches the following selection criteria: WAs in both ultraviolet (UV) and X-rays, suitably variable, bright in UV and X-rays, good archival spectra for comparing with the obscured spectra. To investigate WAs and obscurers of NGC~3227, we used a broadband spectral-energy-distribution model built in our Paper I and the photoionization code of SPEX software to fit archival XMM-Newton and NuSTAR observations in 2006 and 2016. Using unobscured observations, we find four WAs with different ionization states (log$\xi$ [erg cm/s]~-1.0, 2.0, 2.5, 3.0). The highest-ionization WA has a higher hydrogen column density (~$10^{22}$/cm$^2$) than the other three WAs (~$10^{21}$/cm$^2$). Their outflow velocities range from 100 to 1300 km/s, and show a positive correlation with the ionization parameter. These WAs are estimated to be between the outer broad-line-region (BLR) and the narrow line region. Besides, we find an X-ray obscuration event in 2006, which was missed by previous studies. It can be explained by a single obscurer. We also study the previously published obscuration event in 2016, which needs two obscurers in the fit. A high-ionization obscurer (log$\xi$~2.80; covering factor $C_f$~30%) only appears in 2016, which has a high column density (~$10^{23}$/cm$^2$). A low-ionization obscurer (log$\xi$~1.0-1.9; $C_f$~20%-50%) exists in both 2006 and 2016, which has a lower column density (~$10^{22}$/cm$^2$). These obscurers are estimated to be in the BLR by their crossing time of transverse motions. The obscurers and WAs of NGC 3227 have different distances and number densities, which indicate that they might have different origins., Comment: 12 pages, 6 figures, 6 tables; accepted by A&A

Inconsistent conclusions are obtained from recent active galactic nuclei (AGNs) accretion disk inter-band time-lag measurements. While some works show that the measured time lags are significantly larger (by a factor of $\sim 3$) than the theoretical predictions of the Shakura \& Sunyaev disk (SSD) model, others find that the time-lag measurements are consistent with (or only slightly larger than) that of the SSD model. These conflicting observational results might be symptoms of our poor understanding of AGN accretion physics. Here we show that sources with larger-than-expected time lags tend to be less-luminous AGNs. Such a dependence is unexpected if the inter-band time lags are attributed to the light-travel-time delay of the illuminating variable X-ray photons to the static SSD. If, instead, the measured inter-band lags are related not only to the static SSD but also to the outer broad emission-line regions (BLRs; e.g., the blended broad emission lines and/or diffuse continua), our result indicates that the contribution of the non-disk BLR to the observed UV/optical continuum decreases with increasing luminosity ($L$), i.e., an anti-correlation resembling the well-known Baldwin effect. Alternatively, we argue that the observed dependence might be a result of coherent disk thermal fluctuations as the relevant thermal timescale, $\tau_{\mathrm{TH}}\propto L^{0.5}$. With future accurate measurements of inter-band time lags, the above two scenarios can be distinguished by inspecting the dependence of inter-band time lags upon either the BLR components in the variable spectra or the timescales., Comment: 8 pages, 3 figures, 1 table, ApJL accepted

The rest-frame UV/optical variability of the quasars in the Sloan Digital Sky Survey (SDSS) Stripe 82 is used to test the Corona-Heated Accretion-disk Reprocessing (CHAR) model of Sun et al. 2020. We adopt our CHAR model and the observed black-hole masses ($M_{\mathrm{BH}}$) and luminosities ($L$) to generate mock light curves that share the same measurement noise and sampling as the real observations. Without any fine-tuning, our CHAR model can satisfactorily reproduce the observed ensemble structure functions for different $M_{\mathrm{BH}}$, $L$, and rest-frame wavelengths. Our analyses reveal that a luminosity-dependent bolometric correction is disfavored over the constant bolometric correction for UV/optical luminosities. Our work demonstrates the possibility of extracting quasar properties (e.g., the bolometric correction or the dimensionless viscosity parameter) by comparing the physical CHAR model with quasar light curves., Comment: 16 pages, 9 figures, accepted to ApJ

The dusty torus plays a vital role in unifying active galactic nuclei (AGNs). However, the physical structure of the torus remains largely unclear. Here we present a systematical investigation of the torus mid-infrared (MIR) spectroscopic feature, i.e., the 9.7 um silicate line, of $175$ AGNs selected from the Swift/BAT Spectroscopic Survey (BASS). Our sample is constructed to ensure that each of the $175$ AGNs has Spizter/IRS MIR, optical, and X-ray spectroscopic coverage. Therefore, we can simultaneously measure the silicate strength, optical emission lines, and X-ray properties (e.g., the column density and the intrinsic X-ray luminosity). We show that, consistent with previous works, the silicate strength is weakly correlated with the hydrogen column density ($N_\mathrm{H}^\mathrm{X}$), albeit with large scatters. For X-ray unobscured AGNs, the silicate-strength-derived $V$-band extinction and the broad-H$\alpha$-inferred one are both small; however, for X-ray obscured AGNs, the former is much larger than the latter. In addition, we find that the optical type 1 AGNs with strong X-ray absorption on average show significant silicate absorption, indicating that their X-ray absorption might not be caused by dust-free gas in the broad-line region. Our results suggest that the distribution and structure of the obscuring dusty torus are likely to be very complex. We test our results against the smooth and clumpy torus models and find evidence in favor of the clumpy torus model., Comment: 17 pages, 12 figures, Accepted to Research in Astronomy and Astrophysics (a new version with appropriate top margins so that the few bottom lines appear properly)

We investigate the intrinsic spectral energy distribution (SED) of active galactic nuclei (AGNs) at infrared (IR) bands with 42 $z < 0.5$ optically luminous Palomar Green survey quasars through SED decomposition. We decompose the SEDs of the 42 quasars by combining an AGN IR template library Siebenmorgen2015 that covers a wide range of the AGN parameter space with three commonly used galaxy template libraries. We determine the median AGN SED from the best-fitting results. The far-IR (FIR) contribution of our median AGN SED is significantly smaller than that of Symeonidis et al. 2016, but roughly consistent with that of Lyu et al. 2017. The AGN IR SED becomes cooler with increasing bolometric luminosity, which might be due to that more luminous AGNs might have stronger radiative feedback to change torus structures and/or their tori might have higher metallicities. Our conclusions do not depend on the choice of galaxy template libraries. However, since the predicted polycyclic aromatic hydrocarbon (PAH) emission line flux is galaxy template-dependent, cautions should be taken on deriving galaxy FIR contribution from PAH fluxes., Comment: 15 pages, 14 figures, Accepted to ApJ. Tables 2 & 3 and the fitting results can be downloaded from https://www.mysun.work/research/agn-ir-seds/

Active galactic nuclei (AGNs) have long been observed to "twinkle" (i.e., their brightness varies with time) on timescales from days to years in the UV/optical bands. Such AGN UV/optical variability is essential for probing the physics of supermassive black holes (SMBHs), the accretion disk, and the broad-line region. Here we show that the temperature fluctuations of an AGN accretion disk, which is magnetically coupled with the corona, can account for observed high-quality AGN optical light curves. We calculate the temperature fluctuations by considering the gas physics of the accreted matter near the SMBH. We find that the resulting simulated AGN UV/optical light curves share the same statistical properties as the observed ones as long as the dimensionless viscosity parameter $\alpha$, which is widely believed to be controlled by magnetohydrodynamic (MHD) turbulence in the accretion disk, is about $0.01$---$0.2$. Moreover, our model can simultaneously explain the larger-than-expected accretion disk sizes and the dependence of UV/optical variability upon wavelength for NGC 5548. Our model also has the potential to explain some other observational facts of AGN UV/optical variability, including the timescale-dependent bluer-when-brighter color variability and the dependence of UV/optical variability on AGN luminosity and black hole mass. Our results also demonstrate a promising way to infer the black-hole mass, the accretion rate, and the radiative efficiency, thereby facilitating understanding of the gas physics and MHD turbulence near the SMBH and its cosmic mass growth history by fitting the AGN UV/optical light curves in the era of time-domain astronomy., Comment: 25 pages, 25 figures. Accepted to ApJ. The two figsets can be downloaded from https://www.mysun.work/projects/charm/

The spectral energy distribution of blazars around the synchrotron peak can be well described by the log-parabolic model that has three parameters: peak energy ($E_\textrm{p}$), peak luminosity ($L_\textrm{p}$) and the curvature parameter ($b$). It has been suggested that $E_\textrm{p}$ shows relations with $L_\textrm{p}$ and $b$ in several sources, which can be used to constrain the physical properties of the emitting region and/or acceleration processes of the emitting particles. We systematically study the $E_\textrm{p}$-$L_\textrm{p}$ and $E_\textrm{p}$-(1$/b$) relations for 14 BL Lac objects using the 3-25~keV $RXTE$/PCA and 0.3-10~keV $Swift$/XRT data. Most objects (9/14) exhibit positive $E_\textrm{p}$-$L_\textrm{p}$ correlations, three sources show no correlation, and two sources display negative correlations. In addition, most targets (7/14) present no correlation between $E_\textrm{p}$ and 1$/b$, five sources pose negative correlations, and two sources demonstrate positive correlations. 1ES~1959+650 displays two different $E_\textrm{p}$-$L_\textrm{p}$ relations in 2002 and 2016. We also analyze $E_\textrm{p}$-$L_\textrm{p}$ and $E_\textrm{p}$-(1$/b$) relations during flares lasting for several days. The $E_\textrm{p}$-$L_\textrm{p}$ relation does not exhibit significant differences between flares, while the $E_\textrm{p}$-(1$/b$) relation varies from flare to flare. For the total sample, when $L_\textrm{p}$ < $\textrm{10}^\textrm{45}\ \textrm{erg}\ \textrm{s}^\textrm{-1}$, there seems to be a positive $E_\textrm{p}$-$L_\textrm{p}$ correlation. $L_\textrm{p}$ and the slope of $E_\textrm{p}$-$L_\textrm{p}$ relation present an anti-correlation, which indicates that the causes of spectral variations might be different between luminous and faint sources. $E_\textrm{p}$ shows a positive correlation with the black hole mass. We discuss the implications of these results., Comment: 19 pages, 10 figures, 2 tables; Accepted for publication in ApJ

Two bright X-ray transients were reported from the Chandra Deep Field South archival data, namely CDF-S XT1 and XT2. Whereas the nature of the former is not identified, the latter was suggested as an excellent candidate for a rapidly spinning magnetar born from a binary neutron star (BNS) merger. Here we propose a unified model to interpret both transients within the framework of the BNS merger magnetar model. According to our picture, CDF-S XT2 is observed from the "free zone" where the magnetar spindown powered X-ray emission escapes freely, whereas CDF-S XT1 originates from the "trapped zone" where the X-ray emission is initially blocked by the dynamical ejecta and becomes transparent after the ejecta is pushed to a distance where Thomson optical depth drops below unity. We fit the magnetar model to the light curves of both transients and derived consistent parameters for the two events, with magnetic field, initial spin period and X-ray emission efficiency being ($B_p=10^{16}\,G$, $P=1.2\,\rm ms$, $\eta = 0.001$) and ($B_p=10^{15.8}\,G$, $P=4.4\, \rm ms$, $\eta = 0.001$) for XT1 and XT2, respectively. The "isotropic equivalent" ejecta mass of XT1 is $M_{\rm ej} \sim 10^{-3}$ $M_{\odot}$, while it is not constrained for XT2. Our results suggest that more extreme magnetar parameters are required to have XT1 detected from the trapped zone. The model parameters for both events are generally consistent with those derived from SGRB X-ray plateau observations. The host galaxy properties of both transients are also consistent with those of SGRBs. The event rate densities of both XT1 and XT2 are consistent with that of BNS mergers., Comment: 8 pages, 5 figures, accepted for publication in ApJ

We present a detailed X-ray spectral analysis of 1152 AGNs selected in the Chandra Deep Fields (CDFs), in order to identify highly obscured AGNs ($N_{\rm H} > 10^{23}\ \rm cm^{-2}$). By fitting spectra with physical models, 436 (38%) sources with $L_{\rm X} > 10^{42}\ \rm erg\ s^{-1}$ are confirmed to be highly obscured, including 102 Compton-thick (CT) candidates. We propose a new hardness-ratio measure of the obscuration level which can be used to select highly obscured AGN candidates. The completeness and accuracy of applying this method to our AGNs are 88% and 80%, respectively. The observed logN-logS relation favors cosmic X-ray background models that predict moderate (i.e., between optimistic and pessimistic) CT number counts. 19% (6/31) of our highly obscured AGNs that have optical classifications are labeled as broad-line AGNs, suggesting that, at least for part of the AGN population, the heavy X-ray obscuration is largely a line-of-sight effect, i.e., some high-column-density clouds on various scales (but not necessarily a dust-enshrouded torus) along our sightline may obscure the compact X-ray emitter. After correcting for several observational biases, we obtain the intrinsic NH distribution and its evolution. The CT-to-highly-obscured fraction is roughly 52% and is consistent with no evident redshift evolution. We also perform long-term (~17 years in the observed frame) variability analyses for 31 sources with the largest number of counts available. Among them, 17 sources show flux variabilities: 31% (5/17) are caused by the change of NH, 53% (9/17) are caused by the intrinsic luminosity variability, 6% (1/17) are driven by both effects, and 2 are not classified due to large spectral fitting errors., Comment: 32 pages, 21 figures, 9 tables, accepted for publication in ApJ

In this paper we describe the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to accretion flows in the strong field gravity regime around both stellar-mass and supermassive black-holes. eXTP has the unique capability of using advanced 'spectral-timing-polarimetry' techniques to analyze the rapid variations with three orthogonal diagnostics of the flow and its geometry, yielding unprecedented insight into the inner accreting regions, the effects of strong field gravity on the material within them and the powerful outflows which are driven by the accretion process., Comment: Accepted for publication on Sci. China Phys. Mech. Astron. (2019)

Utilizing all the 16-year $RXTE$ observations, we analyze the X-ray spectra of 32 TeV blazars, and perform a systematic investigation of X-ray spectral variability for the 5 brightest sources during their major flares that lasted several days. We obtain photon spectral index ($\alpha$), flux and synchrotron radiation peak energy ($E_{\rm{p}}$) from empirical spectral fitting, and electron spectral index ($p$) from theoretical synchrotron radiation modeling. We find that both $\alpha$ and $p$ generally display a harder-when-brighter trend, confirming the results of many previous works. Furthermore, we confirm and strengthen the result that $p$ must vary in order to explain the observed X-ray spectral variability during flares, which would have useful implications for interpreting the associated higher-energy spectral variability. We see apparent electron spectral hysteresis in many but not all $p$-flux plots that takes a form of "loop" or oblique"8". We obtain a tight $p$-hardness ratio (HR) relation and a tighter $p$-$\alpha$ relation using spectra of flaring periods, both of which are also applicable to stacked data of quiescent periods. We demonstrate that these two empirical relations can be used efficiently to estimate $p$ from HR or $\alpha$ that is readily achieved. Finally, we find that, when considering TeV blazars as a whole, $\alpha$ and X-ray luminosity are positively correlated, $E_{\rm{p}}$ is negatively correlated with $p$ and $\alpha$, and $E_{\rm{p}}$ is positively correlated with HR; all these correlations are in line with the blazar sequence. However, after correcting for the Doppler boosting effect, $\alpha$ and intrinsic X-ray luminosity follow an anti-correlation., Comment: 23 pages, 11 figures, 2 tables. Accepted for publication in ApJ

We explore the evolution of the time variability (in the optical $g$-band and on timescales of weeks to years) of SDSS Stripe 82 quasars along the quasar main sequence. A parent sample of $1004$ quasars within $0.5\leq z \leq 0.89$ are used for our statistical studies, we then make subsamples from our parent sample: a subsample of $246$ quasars with similar luminosities, and a subsample of $399$ quasars with similar Rfe (i.e., the ratio of the equivalent width of FeII within $4435$--$4685\ \mathrm{\AA}$ to that of \Hbeta). We find the variability amplitude decreases with luminosity ($L_{\mathrm{bol}}$). The anti-correlation between the variability amplitude and Rfe is weak but statistically significant. The characteristic timescale, $\tau$, correlates mostly with quasar luminosity, its dependence on Rfe is statistically insignificant. After controlling luminosity and Rfe, the high- and low-FWHM samples have similar structure functions. These results support the framework that Rfe is governed by Eddington ratio and FWHM of Hbeta is mostly determined by orientation. We then provide new empirical relations between variability parameters and quasar properties (i.e., luminosity and Rfe). Our new relations are consistent with the scenario that quasar variability is driven by the thermal fluctuations in the accretion disk, $\tau$ seems to correspond to the thermal timescale. From our new relations, we find the short-term variability is mostly sensitive to $L_{\mathrm{bol}}$. Basing on this, we propose that quasar short-term (a few months) variability might be a new type of "Standard Candle" and can be adopted to probe cosmology., Comment: 13 pages, 8 figures, Accepted for publication in ApJ

Recent multi-band variability studies have revealed that active galactic nucleus (AGN) accretion disc sizes are generally larger than the predictions of the classical thin disc by a factor of $2\sim 3$. This hints at some missing key ingredient in the classical thin disc theory: here, we propose an accretion disc wind. For a given bolometric luminosity, in the outer part of an accretion disc, the effective temperature in the wind case is higher than that in the no-wind one; meanwhile, the radial temperature profile of the wind case is shallower than the no-wind one. In presence of winds, for a given band, blackbody emission from large radii can contribute more to the observed luminosity than the no-wind case. Therefore, the disc sizes of the wind case can be larger than those of the no-wind case. We demonstrate that a model with the accretion rate scaling as $\dot{M}_0 (R/R_{\mathrm{S}})^{\beta}$ (i.e., the accretion rate declines with decreasing radius due to winds) can match both the inter-band time lags and the spectral energy distribution of NGC 5548. Our model can also explain the inter-band time lags of other sources. Therefore, our model can help decipher current and future continuum reverberation mapping observations., Comment: 7 pages, 4 figures, MNRAS in press

NGC 5548 has been intensively monitored by the AGN Space Telescope and Optical Reverberation Mapping collaboration. Approximately after half of the light curves, the correlation between the broad emission lines and the lag-corrected ultraviolet continua becomes weak. This anomalous behavior is accompanied by an increase of soft X-ray emission. We propose a simple model to understand this anomalous behavior, i.e., the corona might fall down, thereby increasing the covering fraction of the inner disk. Therefore, X-ray and extreme ultraviolet emission suffer from spectral variations. The ultraviolet continua variations are driven by both X-ray and extreme ultraviolet variations. Consequently, the spectral variability induced by the falling corona would dilute the correlation between the broad emission lines and the ultraviolet continua. Our model can explain many additional observational facts, including the dependence of the anomalous behavior on velocity and ionization energy. We also show that the time lag and correlation between the X-ray and the ultraviolet variations change as NGC 5548 displays the anomalous behavior. The time lag is dramatically longer than the expectation from disk reprocessing if the anomalous behavior is properly excluded. During the anomalous state, the time lag approaches the light-travel timescale of disk reprocessing albeit with a much weaker correlation. We speculate that the time lag in the normal state is caused by reprocessing of the broad line region gas. As NGC 5548 enters the abnormal state, the contribution of the broad line region gas is smaller; the time lag reflects disk reprocessing. We also discuss alternative scenarios., Comment: 10 pages, 6 figures, Accepted for publication in ApJ

We use the multi-epoch spectra of 362 quasars from the Sloan Digital Sky Survey Reverberation Mapping project to investigate the dependence of the blueshift of CIV relative to MgII on quasar properties. We confirm that high-blueshift sources tend to have low CIV equivalent widths (EWs), and that the low-EW sources span a range of blueshift. Other high-ionization lines, such as HeII, also show similar blueshift properties. The ratio of the line width (measured as both the full-width at half maximum and the velocity dispersion) of CIV to that of MgII increases with blueshift. Quasar variability might enhance the connection between the CIV blueshift and quasar properties (e.g., EW). The variability of the MgII line center (i.e., the wavelength that bisects the cumulative line flux) increases with blueshift. In contrast, the CIV line center shows weaker variability at the extreme blueshifts. Quasars with the high-blueshift CIV lines tend to have less variable continuum emission, when controlling for EW, luminosity, and redshift. Our results support the scenario that high-blueshift sources tend to have large Eddington ratios., Comment: 18 pages, 25 figures, Accepted for publication in ApJ

We present a detailed spectral analysis of the brightest Active Galactic Nuclei (AGN) identified in the 7Ms Chandra Deep Field South (CDF-S) survey over a time span of 16 years. Using a model of an intrinsically absorbed power-law plus reflection, with possible soft excess and narrow Fe K$\alpha$ line, we perform a systematic X-ray spectral analysis, both on the total 7Ms exposure and in four different periods with lengths of 2-21 months. With this approach, we not only present the power-law slopes, column densities $N_H$, observed fluxes, and absorption-corrected 2-10~keV luminosities $L_X$ for our sample of AGNs, but also identify significant spectral variabilities among them on time scales of years. We find that the $N_H$ variabilities can be ascribed to two different types of mechanisms, either flux-driven or flux-independent. We also find that the correlation between the narrow Fe line EW and $N_H$ can be well explained by the continuum suppression with increasing $N_H$. Accounting for the sample incompleteness and bias, we measure the intrinsic distribution of $N_H$ for the CDF-S AGN population and present re-selected subsamples which are complete with respect to $N_H$. The $N_H$-complete subsamples enable us to decouple the dependences of $N_H$ on $L_X$ and on redshift. Combining our data with that from C-COSMOS, we confirm the anti-correlation between the average $N_H$ and $L_X$ of AGN, and find a significant increase of the AGN obscured fraction with redshift at any luminosity. The obscured fraction can be described as $f_{obscured}\thickapprox 0.42\ (1+z)^{0.60}$., Comment: 35 figures, ApJS accepted

Based on high-resolution ultraviolet spectroscopy obtained with FUSE and COS, we present new detections of O VI and N V emission from the black-hole X-ray binary (XRB) system LMC X-3. We also update the ephemeris of the XRB using recent radial velocity measurements obtained with the echelle spectrograph on the Magellan-Clay telescope. We observe significant velocity variability of the UV emission, and we find that the O VI and N V emission velocities follow the optical velocity curve of the XRB. Moreover, the O VI and N V intensities regularly decrease between binary phase = 0.5 and 1.0, which suggests that the source of the UV emission is increasingly occulted as the B star in the XRB moves from superior to inferior conjunction. These trends suggest that illumination of the B-star atmosphere by the intense X-ray emission from the accreting black hole creates a hot spot on one side of the B star, and this hot spot is the origin of the O VI and N V emission. However, the velocity semiamplitude of the ultraviolet emission, K_{UV}~180 km/s, is lower than the optical semiamplitude; this difference could be due to rotation of the B star. If our hypothesis about the origin of the highly ionized emission is correct, then careful analysis of the emission occultation could, in principle, constrain the inclination of the XRB and the mass of the black hole., Comment: To be published in the Astronomical Journal; 10 pages, 9 figures, 3 tables