X-ray Flux Variability Research Articles

Characterization of the Western Pictor A Hotspot in Hard X-Rays with NuSTAR

The origin of X-ray emission from the resolved kiloparsec-scale jets and hotspots of many active galactic nuclei remains uncertain, particularly where the X-ray emission is separate from the radio-optical synchrotron component. Possible explanations include synchrotron emission from a second electron population and external Compton or synchrotron self-Compton processes—alternatives which imply very different physical conditions within the jet. Until recently, X-ray studies of resolved jets and hotspots have been restricted to below ∼10 keV, often showing a hard spectral index indicating a spectral peak beyond this energy range. Here we present NuSTAR observations of the nearby powerful radio galaxy Pictor A, in which we clearly detect the western hotspot at approximately 4′ from the host galaxy, the most significant detection of hotspot emission above 10 keV to date. The NuSTAR spectrum is best fit by a single power law of index Γ = 2.03 ± 0.04; an exponential cutoff gives a 1σ lower limit on the cutoff energy of 40.7 keV. We confirm previous findings of variations in the soft X-ray flux detected by Chandra over the 2000 to 2015 period, at a significance of 6.5σ. This rises to >8σ in the common 3–8 keV band using the combined 22 yr span of Chandra and NuSTAR observations. The variability of the western Pictor A hotspot strongly confirms the previously argued synchrotron nature of the X-ray emission for the hotspot, while the lower bound to the spectral cutoff energy implies electron energies in the hotspot reach up to at least a few TeV.

Intensive Broadband Reverberation Mapping of Fairall 9 with 1.8 yr of Daily Swift Monitoring

We present 1.8 yr of near-daily Swift monitoring of the bright, strongly variable Type 1 active galactic nucleus (AGN) Fairall 9. Totaling 575 successful visits, this is the largest such campaign reported to date. Variations within the UV/optical are well correlated, with longer wavelengths lagging shorter wavelengths in the direction predicted by thin-disk/lamppost models. The correlations are improved by “detrending,” subtracting a second-order polynomial fit to the UV/optical light curves to remove long-term trends that are not of interest to this study. Extensive testing indicates detrending with higher-order polynomials removes too much intrinsic variability signal on reverberation timescales. These data provide the clearest detection to date of interband lags within the UV, indicating that neither emission from a large disk nor diffuse continuum emission from the broad-line region (BLR) can independently explain the full observed lag spectrum. The observed X-ray flux variations are poorly correlated with those in the UV/optical. Further, subdivision of the data into four ∼160 days light curves shows that the UV/optical lag spectrum is highly stable throughout the four periods, but the X-ray to UV lags are unstable, significantly changing magnitude and even direction from one period to the next. This indicates the X-ray to UV relationship is more complex than predicted by the simple reprocessing model often adopted for AGN. A “bowl” model (lamppost irradiation and blackbody reprocessing on a disk with a steep rim) fit suggests the disk thickens at a distance (∼10 lt-day) and temperature (∼8000 K) consistent with the inner edge of the BLR.

Chandra Study of the Proper Motion of HST-1 in the Jet of M87

The radio galaxy M87 is well known for its jet, which features a series of bright knots observable from radio to X-ray wavelengths. We analyze the X-ray image and flux variability of the knot HST-1 in the jet. Our analysis includes all 112 available Chandra ACIS-S observations from 2000 to 2021, with a total exposure time of ∼884 ks. We use deconvolved images to study the brightness profile of the X-ray jet and measure the relative separation between the core and HST-1. From 2003 to 2005 (which coincides with a bright flare from HST-1), we find a correlation between the flux of HST-1 and its offset from the core. In subsequent data, we find a steady increase in this offset, which implies a bulk superluminal motion for HST-1 of 6.6 ± 0.9 c (2.0 ± 0.3 pc yr−1), in keeping with prior results. We discuss models for the flux–offset correlation that feature either two or four emission regions separated by tens of parsecs. We attribute these results to moving shocks in the jet, which allow us to measure the internal structure of the jet.

Ensemble X-ray variability of optically selected QSOs: dependence on black hole mass and Eddington ratio

ABSTRACT Although flux variability is one of the defining properties of accretion flows on to supermassive black holes, its dependence on physical parameters such as the mass of the compact object and the Eddington ratio remains under discussion. In this paper, we address this issue using the structure function statistic to measure the variability at X-ray wavelengths of a sample of optically selected QSOs with available black hole masses and Eddington ratios. We present a new Bayesian methodology for estimating the structure function tailored to the Poisson nature of the X-ray data. This is applied to 15 548 SDSS DRQ16 QSOs with repeat observations in the XMM–Newton archive and/or the SRG/eROSITA All Sky Survey. The X-ray structure function monotonically increases to time intervals of about 10–15 yrs, consistent with scenarios in which instabilities of the accretion disc contribute to the X-ray variability on long time-scales. Additionally, there is evidence that the amplitude of the stochastic X-ray flux variations rises with decreasing black hole mass and Eddington ratio. This finding imposes stringent constraints on empirical models of Active Galactic Nuclei variability derived from local samples, emphasizing the significance of high-redshift population studies for comprehending the stochastic flux variations in active black holes.

Flickering pulsations in bright X-ray pulsars: the evidence of gravitationally lensed and eclipsed accretion column

ABSTRACT It is expected that the extreme mass accretion rate onto strongly magnetized neutron stars results in the appearance of accretion columns above the stellar surface. For a distant observer, rotation of a star results in periodic variations of X-ray flux. Because the mass accretion rate fluctuates around the average value, the pulse profiles are not stable and demonstrate fluctuations as well. In the case of bright X-ray pulsars, however, pulse fluctuations are not solely attributed to variations in the mass accretion rate. They are also influenced by the variable height of the columns, which is dependent on the mass accretion rate. This study delves into the process of pulse profile formation in bright X-ray pulsars, taking into account stochastic fluctuations in the mass accretion rate, the corresponding variations in accretion column geometry, and gravitational bending. Our analysis reveals that potential eclipses of accretion columns by a neutron star during their spin period should manifest specific features in pulse profile variability. Applying a novel pulse profile analysis technique, we successfully detect these features in the bright X-ray transient V 0332+53 at luminosities $\gtrsim 2\times 10^{38}\, {\rm erg\ \rm s^{-1}}$. This detection serves as compelling evidence for the eclipse of an accretion column by a neutron star. Detection of the eclipse places constraints on the relation between neutron star mass, radius, and accretion column height. Specifically, we can establish an upper limit on the accretion column height, which is crucial for refining theoretical models of extreme accretion.

Broad-band noise and quasi-periodic oscillation characteristics of the X-ray pulsar RX J0440.9+4431

ABSTRACT We present a comprehensive timing analysis on the Be/X-ray binary pulsar RX J0440.9+4431 using observations from NICER and Insight-HXMT during the 2022–2023 outburst. The power density spectrum (PDS) of RX J0440.9+4431 exhibits typical aperiodic variability in X-ray flux across a wide frequency range. During a supercritical accretion state, we detect quasi-periodic oscillations (QPOs) at 0.2–0.5 Hz in the light curves of five pulses for RX J0440.9+4431. The observed QPOs manifest during flares, while the flares appear at the peaks of the pulse profiles on a time-scale of seconds and are primarily caused by an increase in hard photons. These flares can be explained by increased material ingestion in the accretion column at a fixed phase, primarily generating hard photons. Alternatively, an increase in accretion rate, independent of phase, may result in highly beamed hard photons within the accretion column, causing the flares. We argue the origin of QPOs to instabilities within the accretion flow. Additionally, we find that the break frequencies in the noise power spectra align well with $\propto L_{\mathrm{x}}^{3 / 7}$ across three orders of magnitude in the luminosity, which points to a relatively strong magnetic field in RX J0440.9+4431, estimated to be ~1013 G.

Evidence for a dynamic corona in the short-term time lags of black hole X-ray binary MAXI J1820+070

ABSTRACT In X-ray observations of hard state black hole X-ray binaries (BHXRBs), rapid variations in accretion disc and coronal power-law emission are correlated and show Fourier-frequency-dependent time lags. On short ($\sim$0.1 s) time-scales, these lags are thought to be due to reverberation and therefore may depend strongly on the geometry of the corona. Low-frequency quasi-periodic oscillations (QPOs) are variations in X-ray flux that have been suggested to arise because of geometric changes in the corona, possibly due to general relativistic Lense–Thirring precession. Therefore, one might expect the short-term time lags to vary on the QPO time-scale. We performed novel spectral-timing analyses on Neutron Star Interior Composition ExploreR observations of the BHXRB MAXI J1820+070 during the hard state of its outburst in 2018 to investigate how the short-term time lags between a disc-dominated and a coronal power-law-dominated energy band vary on different time-scales. Our method can distinguish between variability due to the QPO and broad-band noise, and we find a linear correlation between the power-law flux and lag amplitude that is strongest at the QPO frequency. We also introduce a new method to resolve the QPO signal and determine the QPO phase dependence of the flux and lag variations, finding that both are very similar. Our results are consistent with a geometric origin of QPOs, but also provide evidence for a dynamic corona with a geometry varying in a similar way over a broad range of time-scales, not just the QPO time-scale.

Three cases of optical periodic modulation in Active Galactic Nuclei

ABSTRACT We report on the case of optical periodic modulation discovered in two Active Galactic Nuclei (AGN) and one candidate AGN. Analysing the archival optical data obtained from large transient surveys, namely the Catalina Real-Transient Survey (CRTS) and the Zwicky Transient Facility (ZTF), we find periodicities of 2169.7, 2103.1, and 1462.6 d in sources J0122 + 1032, J1007 + 1248 (or PG 1004 + 1248), and J2131 − 1127, respectively. The optical spectra of the first two indicate that the first is likely a blazar and the second a type 1 Seyfert galaxy, and while no spectroscopic information is available for the third one, its overall properties suggest that it is likely an AGN. In addition, mid-infrared (MIR) light-curve data of the three sources, taken by the Wide-field Infrared Survey Explorer (WISE), are also analysed. The light curves show significant variations, but not appearing related to the optical periodicities. Based on the widely discussed supermassive black hole binary (SMBHB) scenario, we discuss the origin of the optical modulation. Two possible interesting features, an additional 162-d short optical periodicity in J2131 − 1127 and the consistency of the X-ray flux variations of J1007 + 1248 with its optical periodicity, are also discussed within the SMBHB scenario.

THE SOLAR FLARE IMPACT ON THE LOWER IONOSPHERE AND GEOMAGNETIC FIELD

The geophysical testing ground at Kamchatka provides wide possibilities to monitor the electromagnetic and ionospheric phenomena at the Far East. Here we demonstrate the facilities of this multi-instrument observational complex for the study of the solar flare impact on the lower ionosphere and geomagnetic field. We have used data of the VLF amplitude/phase monitoring along several radio paths, augmented by the geomagnetic observations. A typical VLF amplitude/phase response “replicates” the X-ray flux variations, however the magnitude of the response depends on a preexisting ionization level. The geomagnetic field, besides long-lasting (~30-40 min) magnetic bay, sometimes demonstrates the oscillatory (with quasi-periods ~ 4 mins) response to an impulsive X-ray/UV emission burst, whereas these oscillations are absent in the VLF data. We suggest that some ionospheric heavily damped acoustic-gravity or MHD modes may be responsible for this quasi-periodic response. Solar flare was found to enhance the intensity of magnetospheric Pc3 (~ 50 s) pulsations.

A Multiwavelength Study of the Ongoing Changing-look AGN AT2021fxu

We present multiwavelength analyses of an active optical transient AT2021fxu which shows the appearance of previously absent broad emission lines in a recent optical spectrum, suggesting a Changing-look Active Galactic Nucleus (AGN). During the spectral transition, the brightness in the individual photometric bands increased up to ≈0.6 in the optical bands and up to ≈1.1 mag in the UV bands. The brightening was accompanied by a blueward shift of the optical spectrum. AT2021fxu shows high X-ray (0.3–10 keV) flux variability before and after the outburst, with the average X-ray flux increasing by a factor of ≈2 post-outburst. However, the X-ray spectral shape remains roughly the same, with no significant change in the line-of-sight column density. AT2021fxu’s overall properties are consistent with an accretion-rate-driven transition from a Type-II to a Type-1 AGN.

Black hole mass estimation using X-ray variability measurements in Seyfert galaxies

Aims. Our objective is to critically assess the X-ray flux variability as a tool for measuring the black hole (BH) mass in active galactic nuclei (AGN). We aim to establish a prescription for estimating BH masses based on measurements of the normalised excess variance from X-ray data. We discuss the minimum requirements in terms of the light-curve duration and X-ray signal-to-noise ratio (S/N) to enable a reliable determination that is comparable to what can be derived from the continuum and emission line reverberation studies. Methods. We used the light curves of local Seyfert from the Nuclear Spectroscopic Telescope Array hard X-ray mission (NuSTAR), to compute the normalised excess variance (σNXV2) in the 3–10 and 10–20 keV bands, thus extending the analysis to an energy band higher than 10 keV. The excess variance measurements were then combined with independent BH mass estimates from the literature to establish the MBH vs. σNXV2 relation for different samples and weigh its accuracy in terms of the light-curve duration and X-ray S/N. Results. We find that it is possible to accurately measure the BH mass in AGN using excess variance measurements in the 3–10 and the 10–20 keV bands, however, strong quality requirements should be applied. The minimum necessary S/N and duration of the light curves used to compute the excess variance ought to be ∼3 and ∼80 − 100 ks, respectively. We provide a linear relationship between the normalised excess variance and the black hole mass that can be used to estimate the latter, with an average uncertainty of the order of 0.4 − 0.25 dex (depending on the adopted light-curve segment duration). In general, BH mass estimates from 3–10 keV and 10–20 keV band light curves are expected to be similar. The 10–20 keV band is preferred for sources that are heavily absorbed and the 3–10 keV band is preferred for sources that may be dominated by the X-ray reflection component at energies above 10 keV.

Sensitive Chandra coverage of a representative sample of weak-line quasars: revealing the full range of X-ray properties

ABSTRACT We present deeper Chandra observations for weak-line quasars (WLQs) in a representative sample that previously had limited X-ray constraints, and perform X-ray photometric analyses to reveal the full range of X-ray properties of WLQs. Only 5 of the 32 WLQs included in this representative sample remain X-ray undetected after these observations, and a stacking analysis shows that these 5 have an average X-ray weakness factor of >85. One of the WLQs in the sample that was known to have extreme X-ray variability, SDSS J1539+3954, exhibited dramatic X-ray variability again: it changed from an X-ray normal state to an X-ray weak state within ≈3 months in the rest frame. This short time-scale for an X-ray flux variation by a factor of ≳9 further supports the thick disc and outflow (TDO) model proposed to explain the X-ray and multiwavelength properties of WLQs. The overall distribution of the X-ray-to-optical properties of WLQs suggests that the TDO has an average covering factor of the X-ray emitting region of ∼0.5, and the column density of the TDO can range from $N_{\rm H}\, \sim 10^{23-24}~{\rm cm}^{-2}$ to $N_{\rm H}\, \gtrsim 10^{24}~{\rm cm}^{-2}$, which leads to different levels of absorption and Compton reflection (and/or scattering) among WLQs.

Blazar OJ 287 after First VHE Activity: Tracking the Reemergence of the HBL-like Component in 2020

Abstract We report the reemergence of a new broadband emission through a detailed and systematic study of the multiwavelength spectral and temporal behavior of OJ 287 after its first-ever reported very high energy activity in 2017 to date, which includes the second-highest X-ray flux of the source. The source shows high optical to X-ray flux variations, accompanied mainly by strong spectral changes. The optical to X-ray flux variations are correlated and simultaneous except for two durations when they are anticorrelated. The flux variations, however, are anticorrelated with the X-ray spectral state while correlated with optical–UV (ultraviolet). Weekly binned Fermi-LAT data around the duration of the highest X-ray activity show a few detections with a log-parabola model but none with a power law, yet the extracted LAT spectral energy distribution of the high-activity duration for both the models is similar and shows a hardening above 1 GeV. Further, near-infrared data indicate strong spectral change, resembling a thermal component. Overall, the combined optical to gamma-ray broadband spectrum establishes the observed variations to a new high-energy-peaked broadband emission component, similar to the one seen during the highest reported X-ray flux state of the source in 2017. The observed activities indicate some peculiar features that seem to be characteristic of this emission component, while its appearance a few years around the claimed ∼12 yr optical outbursts strongly indicates a connection between the two.

Light-curve fingerprints: an automated approach to the extraction of X-ray variability patterns with feature aggregation – an example application to GRS 1915+105

ABSTRACT Time series data mining is an important field of research in the era of ‘Big Data’. Next generation astronomical surveys will generate data at unprecedented rates, creating the need for automated methods of data analysis. We propose a method of light-curve characterization that employs a pipeline consisting of a neural network with a long-short term memory variational autoencoder architecture and a Gaussian mixture model. The pipeline performs extraction and aggregation of features from light-curve segments into feature vectors of fixed length that we refer to as light-curve ‘fingerprints’. This representation can be readily used as input of down-stream machine learning algorithms. We demonstrate the proposed method on a data set of Rossi X-ray Timing Explorer observations of the Galactic black hole X-ray binary GRS 1915+105, which was chosen because of its observed complex X-ray variability. We find that the proposed method can generate a representation that characterizes the observations and reflects the presence of distinct classes of GRS 1915+105 X-ray flux variability. We find that this representation can be used to perform efficient classification of light curves. We also present how the representation can be used to quantify the similarity of different light curves, highlighting the problem of the popular classification system of GRS 1915+105 observations, which does not account for intermediate class behaviour.

Spectro-temporal Studies of Rapid Transition of the Quasi-periodic Oscillations in the Black Hole Source H1743-322

Abstract An appearance or disappearance of quasi-periodic oscillations (QPOs) associated with the variation of X-ray flux can be used to decipher the accretion–ejection mechanism of black hole X-ray sources. We searched for and studied such rapid transitions in H1743-322 using archival data from the Rossi X-ray Timing Explorer and found eight such events, where a QPO vanishes suddenly along with the variation of X-ray flux. The appearance/disappearance of QPOs was associated with four events exhibiting type-B QPOs at ∼4.5 Hz: one was a type-A QPO at ν ∼ 3.5 Hz, and the other three were connected to type-C QPOs at ∼9.5 Hz. Spectral studies of the data revealed that the inner disk radius remained at the same location around 2–9 r g , depending on the model used, but power-law indices were varying, indicating that either a corona or a jet is responsible for the events. The probable ejection radii of coronas were estimated to be around 4.2–15.4 r g based on the plasma ejection model. Our X-ray and quasi-simultaneous radio correlation studies suggest that the type-B QPOs are probably related to the precession of a weak jet, though a small and weak corona is present at its base, and the type-C QPOs are associated with the base of a relatively strong jet, which is acting like a corona.

A Time Resolved View of the X-Ray Spectral Variability of the Seyfert-1 Galaxy ESO 511-G030

Abstract The origin of the Soft X-ray excess in the spectra of several nearby active galactic nucleus is still highly debated. In this work we report the monotonic spectral hardening of a nearby bare Seyfert galaxy ESO 511-G030 over a time period of ∼110 ks. We found that both the soft (0.3–2 keV) and the hard X-ray (2–10 keV) bands have shown interesting variability during the 110 ks observation. The soft X-ray flux shows a sharp drop in the last ∼50 ks, by ∼30%, while the hard X-ray flux initially increased and then decreased. The overall effect of these variations is a monotonic increase in the hardness ratio. We conclude that: (A) The soft X-ray flux variations may be related to the changes in the primary hard X-ray. (B) The soft X-ray flux decreases by ∼36% in a timespan of ∼50 ks, indicating that its origin could be from a region as close as <0.5 light-day from the supermassive black hole. A detailed hardness ratio resolved spectral analysis with state-of-the-art reflection models will be reported in a future paper.

A machine learning approach for classification of accretion states of black hole binaries

ABSTRACT In this paper, we employ Machine Learning algorithms on multimission observations for the classification of accretion states of outbursting black hole X-ray binaries for the first time. Archival data from RXTE, Swift, MAXI, and AstroSat observatories are used to generate the hardness intensity diagrams (HIDs) for outbursts of the sources XTE J1859+226 (1999 outburst), GX 339−4 (2002, 2004, 2007, and 2010 outbursts), IGR J17091−3624 (2016 outburst), and MAXI J1535−571 (2017 outburst). Based on variation of X-ray flux, hardness ratios, presence of various types of quasi-periodic oscillations (QPOs), photon indices, and disc temperature, we apply clustering algorithms like K-Means clustering and Hierarchical clustering to classify the accretion states (clusters) of each outburst. As multiple parameters are involved in the classification process, we show that clustering algorithms club together the observations of similar characteristics more efficiently than the ‘standard’ method of classification. We also infer that K-Means clustering provides more reliable results than Hierarchical clustering. We demonstrate the importance of the classification based on machine learning by comparing it with results from ‘standard’ classification.

Analysis of Global Ionospheric Response to Solar Flares Based on Total Electron Content and Very Low Frequency Signals

The response of the global ionosphere to solar flares is an important topic in the field of space weather. The global ionospheric response to solar flares is comprehensively analyzed from the perspectives of total electron content (TEC) and very low frequency (VLF) signals by using solar flare data on the eruption days of X-class flares from 2006 to 2019, including flare level, flare duration, geographical location, and local time. In addition, the relationship between X-ray flux and VLF phase variation is studied through correlation analysis. The concepts of disturbance intensity (DI) and disturbance angle are defined, and a DI evaluation model is established to help detect the difference in sensitivity to solar flares between TEC and VLF signals. Results show the following. (1) The higher the flare level and the longer the duration, the greater the disturbance to the ionosphere. Simultaneously, a phenomenon exists in which the disturbance caused by low-level and long-lasting flares is greater than that caused by high-level flares. (2) VLF phase variation and flare level exhibit a good correlation, and they are also closely related to geographical location and local time. The disturbance degree of a station facing the sun is more evident than that of a station facing away from the sun. The TEC disturbance of stations in the morning (local time) is more obvious than that in the afternoon, and disturbance increases along the direction of Earth’s rotation. (3) When DI is at the same level, the lowest flare level that causes TEC response is higher than the lowest flare level that causes VLF signal response. The disturbance angle of TEC is unevenly distributed within the interval [0°, 90°], and that of VLF signals is more than 85°. The sensitivity of VLF signals to flare response is considerably higher than that of TEC, and the difference between the two even stride across DI level.

UV/Optical Disk Thermal Reverberation in Active Galactic Nuclei: An In-depth Study with an Analytic Prescription for Time-lag Spectra

Abstract Several active galactic nuclei show correlated variations in the UV/optical range, with time delays increasing at longer wavelengths. Thermal reprocessing of the X-rays illuminating the accretion disk has been proposed as a viable explanation. In this scenario, the variable X-ray flux irradiating the accretion disk is partially reflected in X-rays and partially absorbed, thermalized, and reemitted with some delay by the accretion disk at longer wavelengths. We investigate this scenario assuming an X-ray pointlike source illuminating a standard Novikov–Thorne accretion disk around a rotating black hole. We consider all special and general relativistic effects to determine the incident X-ray flux on the disk and in propagating light from the source to the disk and to the observer. We also compute the disk reflection flux, taking into consideration the disk ionization. We investigate the dependence of the disk response function and time lags on various physical parameters, such as the black hole mass and spin; X-ray corona height, luminosity, and photon index; accretion rate; inclination; and inner/outer disk radii. We find it is important to consider relativistic effects and the disk ionization in estimating the disk response. We also find a strong nonlinearity between the X-ray luminosity and the disk response. We present an analytic function for the time-lag dependence on wavelength, which can be used to fit observed time-lag spectra. We also estimate the fraction of the reverberation signal with respect to the total flux, and we suggest possible explanations for the lack of X-ray–UV/optical correlated variations in a few sources.

The Kinematic Evolution of Erupting Structures in Confined Solar Flares

Abstract In this Letter, we study the kinematic properties of ascending hot blobs associated with confined flares. Taking advantage of high-cadence extreme-ultraviolet images provided by the Atmospheric Imaging Assembly on board the Solar Dynamics Observatory, we find that for the 26 events selected here, the hot blobs are first impulsively accelerated outward, but then quickly slow down to motionlessness. Their velocity evolution is basically synchronous with the temporal variation of the Geostationary Operational Environmental Satellite soft X-ray flux of the associated flares, except that the velocity peak precedes the soft X-ray peak by minutes. Moreover, the duration of the acceleration phase of the erupting blobs is moderately correlated with that of the flare rise phase. For nine of the 26 cases, the erupting blobs even appear minutes prior to the onset of the associated flares. Our results show that a fraction of confined flares also involve the eruption of a magnetic flux rope, which sometimes is formed and heated prior to the flare onset. We suggest that the initiation and development of these confined flares are similar to that of eruptive ones, and the main difference may lie in the background field constraint, which is stronger for the former than for the latter.