Hard X-ray Variability Research Articles

Investigating the interplay between the coronal properties and the hard X-ray variability of active galactic nuclei with NuSTAR

Investigating the interplay between the coronal properties and the hard X-ray variability of active galactic nuclei with <i>NuSTAR</i>

Long-term hard X-ray variability properties of Swift-BAT blazars

ABSTRACT We present results from the first dedicated study in the time domain of the hard X-ray variability behaviour of blazars on long time-scales based on ∼13 yr of continuous hard X-ray data in the 14–195 keV band. We use monthly binned data from the recent 157-month Swift-BAT (Burst Alert Telescope) catalogue to characterize the hard X-ray variability of 127 blazars and search for potential differences between the variability of BL Lacertae objects (BL Lacs) and flat-spectrum radio quasars (FSRQs). A significant portion of the blazars in the sample (∼37 per cent) do not show statistically significant hard X-ray variability on monthly time-scales, which is deeply at odds with previous studies that show that blazars are highly variable in the X-rays and other energy bands on a wide range of time-scales. We also find that, on average, the FSRQs and BL Lacs for which we do detect variability exhibit similar flux variability; this suggests that the variability in these FSRQs is not necessarily driven by variations in the source function of scattered external radiation arriving from extended regions, and that it is instead possibly driven by processes that lead to variations in particle injection. In addition, only five blazars in our sample show significant spectral variability in the long-term light curves. For three blazars, we find that a power law that changes slope on monthly time-scales is sufficient to characterize the variable hard X-ray spectrum, suggesting that, at least for some bright blazars, the long-term spectra in the hard X-rays may be described in a relatively simple fashion.

Achromatic rapid flares in hard X-rays in the γ-ray binary LS I + 61 303

ABSTRACT We report on the presence of very rapid hard X-ray variability in the γ-ray binary LS I + 61 303. The results were obtained by analysing Nuclear Spectroscopic Telescope Array (NuSTAR) data, which show two achromatic strong flares on ks time-scales before apastron. The Swift/Burst Alert Telescope orbital X-ray light curve is also presented, and the NuSTAR data are put in the context of the system orbit. The spectrum and estimated physical conditions of the emitting region indicate that the radiation is synchrotron emission from relativistic electrons, likely produced in a shocked pulsar wind. The achromaticity suggests that losses are dominated by escape or adiabatic cooling in a relativistic flow, and the overall behaviour in hard X-rays can be explained by abrupt changes in the size of the emitting region and/or its motion relative to the line of sight, with Doppler boosting potentially being a prominent effect. The rapid changes of the emitter could be the result of different situations such as quick changes in the intra-binary shock, variations in the re-accelerated shocked pulsar wind outside the binary, or strong fluctuations in the location and size of the Coriolis shock region. Although future multiwavelength observations are needed to further constrain the physical properties of the high-energy emitter, this work already provides important insight into the complex dynamics and radiation processes in LS I + 61 303.

Correlated short time-scale hard-soft X-ray variability of the blazars Mrk 421 and 1ES 1959+650 using AstroSat

ABSTRACT We study simultaneous soft (0.7–7 keV) and hard (7–20 keV) X-ray light curves at a total of eight epochs during 2016–2019 of two TeV blazars Mrk 421 and 1ES 1959+650 observed by the SXT and LAXPC instruments on-board AstroSat. The light curves are 45–450 ks long and may be sampled with time bins as short as 600–800 s with high signal-to-noise ratio. The blazars show a harder when brighter trend at all epochs. Discrete cross-correlation functions indicate that the hard and soft X-ray variability are strongly correlated. The time lag is consistent with zero in some epochs, and indicates hard or soft lag of a few hours in the rest. In the leptonic model of blazar emission, soft lag may be due to slower radiative cooling of lower energy electrons while hard lag may be caused by gradual acceleration of the high energy electrons emitting at the hard X-ray band. Assuming the above scenario and the value of the Doppler factor (δ) to be 10–20, the hard and soft lags may be used to estimate the magnetic field to be ∼0.1 Gauss and the acceleration parameter to be ∼104 in the emission region. Due to the availability of the high time resolution (∼ minutes to hours) light curves from AstroSat, the value of the illusive acceleration parameter could be estimated, which provides a stringent constraint on the theories of particle acceleration in blazar jets.

Unusual Hard X-Ray Flares Caught in NICER Monitoring of the Binary Supermassive Black Hole Candidate AT2019cuk/Tick Tock/SDSS J1430+2303

The nuclear transient AT2019cuk/Tick Tock/SDSS J1430+2303 has been suggested to harbor a supermassive black hole (SMBH) binary near coalescence. We report results from high-cadence NICER X-ray monitoring with multiple visits per day from 2022 January to August, as well as continued optical monitoring during the same time period. We find no evidence of periodic/quasiperiodic modulation in the X-ray, UV, or optical bands; however, we do observe exotic hard X-ray variability that is unusual for typical active galactic nuclei (AGN). The most striking feature of the NICER light curve is repetitive hard (2–4 keV) X-ray flares that result in distinctly harder X-ray spectra compared to the nonflaring data. In its nonflaring state, AT2019cuk looks like a relatively standard AGN, but it presents the first case of day-long, hard X-ray flares in a changing-look AGN. We consider a few different models for the driving mechanism of these hard X-ray flares, including (1) corona/jet variability driven by increased magnetic activity, (2) variable obscuration, and (3) self-lensing from the potential secondary SMBH. We prefer the variable corona model, as the obscuration model requires rather contrived timescales and the self-lensing model is difficult to reconcile with a lack of clear periodicity in the flares. These findings illustrate how important high-cadence X-ray monitoring is to our understanding of the rapid variability of the X-ray corona and necessitate further high-cadence, multiwavelength monitoring of changing-look AGN like AT2019cuk to probe the corona-jet connection.

Investigating non-linear and stochastic hard X-ray variability of active galactic nuclei using recurrence analysis

ABSTRACT We present results of recurrence analysis of 46 active galactic nuclei (AGNs) using light curves from the 157-month catalogue of the Swift Burst Alert Telescope (BAT) in the 14–150 keV band. We generate recurrence plots and compute recurrence plot metrics for each object. We use the surrogate data method to compare all derived recurrence-based quantities to three sets of stochastic light curves with identical power spectrum, flux distribution, or both, in order to determine the presence of determinism, non-linearity, entropy, and non-stationarity. We compare these quantities with known physical characteristics of each system, such as black hole mass, Eddington ratio, and bolometric luminosity, radio loudness, obscuration, and spectroscopic type. We find that almost all AGNs in this sample exhibit substantial higher order modes of variability than is contained in the power spectrum, with approximately half exhibiting non-linear or non-stationary behaviour. We find that Type 2 AGNs are more likely to contain deterministic variability than Type 1 AGNs, while the same distinction is not found between obscured and unobscured AGNs. The complexity of variability among Type 1 AGNs is anticorrelated with Eddington ratio, while no relationship is found among Type 2 AGNs. The connections between the recurrence properties and AGN class suggest that hard X-ray emission is a probe of distinct accretion processes among classes of AGNs, which supports interpretations of changing-look AGNs and challenges the traditional unification model that classifies AGNs only on viewing angle.

Daily variability at milli-arcsecond scales in the radio-quiet NLSy1 Mrk 110

ABSTRACT The origin of radio emission in the majority of active galactic nuclei is still poorly understood. Various competing mechanisms are likely involved in the production of radio emission and precise diagnostic tools are needed to disentangle them, of which variability is among the most powerful. For the first time, we show evidence for significant radio variability at 5 GHz at milli-arcsecond scales on days to weeks time-scales in the highly accreting and extremely radio-quiet (RQ) Narrow Line Seyfert 1 (NLSy1) Mrk 110. The simultaneous Swift/XRT light curve indicates stronger soft than hard X-ray variability. The short-term radio variability suggests that the GHz emitting region has a size smaller than ∼180 Schwarzschild radii. The high brightness temperature and the radio and X-ray variability rule out a star-formation and a disc wind origin. Synchrotron emission from a low-power jet and/or an outflowing corona is then favoured.

A Broadband X-Ray View of the Precessing Accretion Disk and Pre-eclipse Dip in the Pulsar Her X-1 with NuSTAR and XMM-Newton

Abstract We present a broadband X-ray timing study of the variations in pulse behavior with the superorbital cycle in the low-mass X-ray binary Her X-1. This source shows a 35 day superorbital modulation in X-ray flux that is likely caused by occultation by a warped, precessing accretion disk. Our data set consists of four joint XMM-Newton and NuSTAR observations of Her X-1 which sample a complete superorbital cycle. We focus our analysis on the first and fourth observations, which occur during the bright “main-on” phase, because these observations have strongly detected pulsations. We added an archival XMM-Newton observation during the “short-on” phase of the superorbital cycle because our observations at that phase are lower in signal to noise. We find that the energy-resolved pulse profiles show the same shape at similar superorbital phases, and the profiles are consistent with expectations from a precessing disk. We demonstrate that a simple precessing accretion disk model is sufficient to reproduce the observed pulse profiles. The results of this model suggest that the similarities in the observed pulse profiles are due to reprocessing by a precessing disk that has returned to its original precession phase. We determine that the broadband spectrum is well fit by an absorbed power law with a soft blackbody component and show that the spectral continuum also exhibits a dependence on the superorbital cycle. We also present a brief analysis of the energy-resolved light curves of a pre-eclipse dip, which shows soft X-ray absorption and hard X-ray variability during the dip.

Broadband X-ray observations of four gamma-ray narrow-line Seyfert 1 galaxies

Narrow-line Seyfert 1 galaxies (NLS1s) is one of the few classes of active galactic nuclei harboring powerful relativistic jets and detected in γ-rays. NLS1s are well-known X-ray sources. While in non-jetted sources the origin of this X-ray emission may be a hot corona surrounding the accretion disk, in jetted objects, especially beamed ones, the contribution of corona and relativistic jet is difficult to disentangle without a proper sampling of the hard X-ray emission. For this reason, we observed with NuSTAR the first four NLS1s detected at high energy γ-rays. These data, along with XMM-Newton and Swift/XRT observations, confirmed that X-rays originate both in the jet and in the accretion disk corona. Time variability in hard X-rays furthermore suggests that, as observed in flat-spectrum radio quasars, the dissipation region during flares could change its position from source to source, and it can be located both inside and outside the broad-line region. We find that jetted NLS1s, and other blazars as well, seem not to follow the classical fundamental plane of black hole (BH) activity, which therefore should be used as a BH mass estimator in blazars with extreme care only. Our results strengthen the idea according to which γ-NLS1s are smaller and younger version of flat-spectrum radio quasars, in which both a Seyfert and a blazar component co-exist.

Hard X-ray properties of NuSTAR blazars

Context. Investigation of the hard X-ray emission properties of blazars is key to the understanding of the central engine of the sources and associated jet process. In particular, simultaneous spectral and timing analyses of the intraday hard X-ray observations provide us a means to peer into the compact innermost blazar regions that are not accessible to our current instruments. Aims. The primary objective of the work is to associate the observed hard X-ray variability properties in blazars with their flux and spectral states, thereby, based on the correlation among these states, extract the details about the emission regions and processes occurring near the central engine. Methods. We carried out timing, spectral, and cross-correlation analysis of 31 NuSTAR observations of 13 blazars. We investigated the spectral shapes of the sources using single power-law, broken power-law, and log-parabola models. We also studied the co-relation between the soft and hard emission using z-transformed discrete correlation function. In addition, we attempted to constrain the smallest emission regions using minimum variability timescales derived from the light curves. Results. We found that, for most of the sources, the hard X-ray emission can be well represented by the log-parabola model and that the spectral slopes for different blazar subclasses are consistent with the so-called blazar sequence. We also report the steepest spectra (Γ ∼ 3) in the BL Lacertae PKS 2155–304 and the hardest spectra (Γ ∼ 1.4) in the flat-spectrum radio quasar PKS 2149–306. In addition, we noted a close connection between the flux and spectral slope within the source subclass in the sense that high flux and/or flux states tend to be harder in spectra. In BL Lacertae objects, assuming particle acceleration by diffusive shocks and synchrotron cooling as the dominant processes governing the observed flux variability, we constrain the magnetic field of the emission region to be a few Gauss; whereas in flat-spectrum radio quasars, using external Compton models, we estimate the energy of the lower end of the injected electrons to be a few Lorentz factors.

Nuclear absorption and emission in the AGN merger NGC 6240 : the hard X-ray view

We present the analysis of four NuSTAR observations of the luminous infrared galaxy merger NGC 6240 , hosting a close pair of highly obscured active galactic nuclei (AGNs). Over a period of about 2 yr, the source exhibits hard X-ray variability of the order of 20 per cent, peaking around 20 keV. When the two AGNs are resolved with Chandra , column densities in the range of NH ∼ 1–2 × 1024 cm−2 are estimated for both of them. The exact values are hard to determine, as they appear to depend on aspects that are sometimes overlooked in Compton-thick objects, such as the covering factor of the absorber, iron abundance and the contamination in the Fe K band from foreground hot-gas emission. Nearly spherical covering and slightly subsolar iron abundance are preferred in this case. While the southern nucleus is suggested to be intrinsically more powerful, as also implied by the mid-IR and 2–10 keV brightness ratios, solutions involving a similar X-ray luminosity of the two AGNs cannot be ruled out. The observed variability is rather limited compared to the one revealed by the Swift /BAT light curve, and it can be fully explained by changes in the continuum flux from the two AGNs, without requiring significant column density variations. NGC 6240 is hereby confirmed to represent a unique opportunity to investigate the X-ray (and broad-band) properties of massive galaxy mergers, which were much more frequent in the early Universe.

Detection of the high-energy cut-off from the Seyfert 1.5 galaxy NGC 5273

We perform the NuSTAR and Swift/XRT joint energy spectral fitting of simultaneous observations from the broad-line Seyfert 1.5 galaxy NGC 5273. When fitted with the combination of an exponential cut-off power-law and a reflection model, a high energy cut-off is detected at 143$^{+96}_{-40}$ keV with 2-sigma significance. Existence of such cut-off is also consistent with the observed Comptonizing electron temperature when fitted with a Comptonization model independently. We observe a moderate hard X-ray variability of the source over the time-scale of ~12 years using INTEGRAL/ISGRI observations in the energy range of 20-100 keV. When the hard band count rate (6-20 keV) is plotted against the soft band count rate (3-6 keV), a hard offset is observed. Our results indicate that the cut-off energy may not correlate with the coronal X-ray luminosity in a simple manner. Similarities in parameters that describe coronal properties indicate that the coronal structure of NGC 5273 may be similar to that of the broad-line radio galaxy 3C 390.3 and another galaxy MCG-5-23-16 where the coronal plasma is dominated by electrons, rather than electron-positron pairs. Therefore, the coronal cooling is equally efficient to the heating mechanism keeping the cut-off energy at low even at the low accretion rate.

Hard X-ray variability of V404 Cygni during the 2015 outburst

Aims. Hard X-ray spectra of black hole binaries (BHB) are produced by Comptonization of soft seed photons by hot electrons near the black hole. The slope of the resulting ray spectra is governed by two main parameters: the electron temperature and the optical depth of the emitting plasma. Given the extreme brightness of V404~Cyg during the 2015 outburst, we aim to constrain the source spectral properties using an unprecedented time resolution in hard X-rays, and to monitor their evolution over the outburst. Methods. We have extracted and analyzed 602 X-ray spectra of V404 Cyg obtained by the IBIS/ISGRI instrument on-board INTEGRAL during the 2015 June outburst, using effective integration times ranging between 8 and 176000 seconds. We have fit the resulting spectra in the 20-200 keV energy range. Results. We find that while the light curve and soft X-ray spectra of V404 Cyg are remarkably different from those of other BHBs, the spectral evolution of V404~Cyg in hard X-rays and the relations between the spectral parameters are consistent with those observed in other BHBs. We identify a hard branch} where the electron temperature is anti-correlated with the hard X-ray flux, and a soft flaring branch where the relation reverses. In addition, we find that during long X-ray plateaus detected at intermediate fluxes, the thermal Comptonization models fail to describe the spectra. However, the statistics improve if we allow the absorption to vary freely in the fits to these spectra. Conclusions. We conclude that the hard branch in V404 Cyg is analogous to the canonical hard state of BHBs. V404~Cyg never seems to enter the canonical soft state, although the soft flaring branch bears resemblance to the BHB intermediate state and ultra-luminous state. The X-ray plateaus are likely the result of absorption by a Compton-thick outflow which reduces the observed flux by a factor of about 10.[abridged]

HIGH ENERGY SPECTRAL EVOLUTION OF V404 CYGNI DURING THE 2015 JUNE OUTBURST AS OBSERVED BY INTEGRAL

The black hole binary GS 2023+338 exhibited an unprecedently bright outburst on June 2015. Since June 17th, the high energy instruments on board INTEGRAL detected an extremely variable emission during both bright and low luminosity phases, with dramatic variations of the hardness ratio on time scales of ~seconds. The analysis of the IBIS and SPI data reveals the presence of hard spectra in the brightest phases, compatible with thermal Comptonization with temperature kTe ~ 40 keV. The seed photons temperature is best fit by kT0 ~ 7 keV, that is too high to be compatible with blackbody emission from the disk. This result is consistent with the seed photons being provided by a different source, that we hypothesize to be a synchrotron driven component in the jet. During the brightest phase of flares, the hardness shows a complex pattern of correlation with flux, with a maximum energy released in the range 40-100 keV. The hard X-ray variability for E > 50 keV is correlated with flux variations in the softer band, showing that the overall source variability cannot originate entirely from absorption, but at least part of it is due to the central accreting source.

FIRSTNuSTAROBSERVATIONS OF MRK 501 WITHIN A RADIO TO TeV MULTI-INSTRUMENT CAMPAIGN

We report on simultaneous broadband observations of the TeV-emitting blazar Markarian 501 between 1 April and 10 August 2013, including the first detailed characterization of the synchrotron peak with Swift and NuSTAR. During the campaign, the nearby BL Lac object was observed in both a quiescent and an elevated state. The broadband campaign includes observations with NuSTAR, MAGIC, VERITAS, the Fermi Large Area Telescope (LAT), Swift X-ray Telescope and UV Optical Telescope, various ground-based optical instruments, including the GASP-WEBT program, as well as radio observations by OVRO, Mets\"ahovi and the F-Gamma consortium. Some of the MAGIC observations were affected by a sand layer from the Saharan desert, and had to be corrected using event-by-event corrections derived with a LIDAR (LIght Detection And Ranging) facility. This is the first time that LIDAR information is used to produce a physics result with Cherenkov Telescope data taken during adverse atmospheric conditions, and hence sets a precedent for the current and future ground-based gamma-ray instruments. The NuSTAR instrument provides unprecedented sensitivity in hard X-rays, showing the source to display a spectral energy distribution between 3 and 79 keV consistent with a log-parabolic spectrum and hard X-ray variability on hour timescales. None (of the four extended NuSTAR observations) shows evidence of the onset of inverse-Compton emission at hard X-ray energies. We apply a single-zone equilibrium synchrotron self-Compton model to five simultaneous broadband spectral energy distributions. We find that the synchrotron self-Compton model can reproduce the observed broadband states through a decrease in the magnetic field strength coinciding with an increase in the luminosity and hardness of the relativistic leptons responsible for the high-energy emission.

THE VIOLENT HARD X-RAY VARIABILITY OF MRK 421 OBSERVED BYNuSTARIN 2013 APRIL

The well studied blazar Markarian 421 (Mrk 421, $z$=0.031) was the subject of an intensive multi-wavelength campaign when it flared in 2013 April. The recorded X-ray and very high energy (VHE, E$>$100 GeV) $\gamma$-ray fluxes are the highest ever measured from this object. At the peak of the activity, it was monitored by the hard X-ray focusing telescope {\it Nuclear Spectroscopic Telescope Array} ({\it NuSTAR}) and {\it Swift} X-Ray Telescope (XRT). In this work, we present a detailed variability analysis of {\it NuSTAR} and {\it Swift}-XRT observations of Mrk 421 during this flaring episode. We obtained the shortest flux doubling time of 14.01$\pm$5.03 minutes, which is the shortest hard X-ray (3$-$79 keV) variability ever recorded from Mrk 421 and is on the order of the light crossing time of the black hole's event horizon. A pattern of extremely fast variability events superposed on slowly varying flares is found in most of the {\it NuSTAR} observations. We suggest that these peculiar variability patterns may be explained by magnetic energy dissipation and reconnection in a fast moving compact emission region within the jet. Based on the fast variability, we derive a lower limit on the magnetic field strength of $B \ge 0.73 \delta_1^{-2/3} \, \nu_{19}^{1/3}$~G, where $\delta_1$ is the Doppler factor in units of 10, and $\nu_{19}$ is the characteristic X-ray synchrotron frequency in units of $10^{19}$~Hz.

Soft X-Ray Pulsations in Solar Flares

The soft X-ray emissions of solar flares come mainly from the bright coronal loops at the highest temperatures normally achieved in the flare process. Their ubiquity has led to their use as a standard measure of flare occurrence and energy, although the overwhelming bulk of the total flare energy goes elsewhere. Recently Dolla et al. (2012) noted quasi-periodic pulsations (QPP) in the soft X-ray signature of the X-class flare SOL2011-02-15, as observed by the standard photometric data from the GOES (Geostationary Operational Environmental Satellite) spacecraft. In this paper we analyze the suitability of the GOES data for this kind of analysis and find them to be generally valuable after September, 2010 (GOES-15). We then extend the Dolla et al. result to a complete list of X-class flares from Cycle 24, and show that most of them (80%) display QPPs in the impulsive phase. The pulsations show up cleanly in both channels of the GOES data, making use of time-series of irradiance differences (the digital time derivative on the 2-s sampling). We deploy different techniques to characterize the periodicity of GOES pulsations, considering the red-noise properties of the flare signals, finding a range of chracteristic time scales of the QPPs for each event, but usually with no strong signature of a single period dominating in the power spectrum. The QPP may also appear on somewhat longer time scales during the later gradual phase, possibly with a greater tendency towards coherence, but the sampling noise in GOES difference data for large irradiance values (X-class flares) makes these more uncertain. We show that there is minimal phase difference between the differenced GOES energy channels, or between them and the hard X-ray variations on short time scales. During the impulsive phase the footpoints of the newly-forming flare loops may also contribute to the observed soft X-ray variations.

Origin of the X-ray off-states in Vela X-1

$Context$: Vela X-1 is the prototype of the classical sgHMXB systems. Recent continuous and long monitoring campaigns revealed a large hard X-rays variability amplitude with strong flares and off-states. This activity has been interpreted invoking clumpy stellar winds and/or magnetic gating mechanisms. $Aims$: We are probing if the observed behaviour could be explained by unstable hydrodynamic flows close to the neutron star rather than the more exotic phenomena. $Methods$: We have used the hydrodynamic code VH-1 to simulate the flow of the stellar wind with high temporal resolution and to compare the predicted accretion rate with the observed light-curves. $Results$: The simulation results are similar to the observed variability. Off-states are predicted with a duration of 5 to 120 minutes corresponding to transient low density bubbles forming around the neutron star. Oscillations of the accretion rate with a typical period of $\sim$ 6800 sec are generated in our simulations and observed. They correspond to the complex motion of a bow shock, moving either towards or away from the neutron star. Flares are also produced by the simulations up to a level of $10^{37}$ erg/s. $Conclusions$: We have qualitatively reproduced the hard X-ray variations observed in Vela X-1 with hydrodynamic instabilities predicted by a simple model. More sophisticated phenomena, such as clumpy winds or the magnetic gating mechanism, are not excluded but not required to explain the basic phenomenology.

THE FIRST HARD X-RAY POWER SPECTRAL DENSITY FUNCTIONS OF ACTIVE GALACTIC NUCLEUS

We present results of our Power Spectral Density (PSD) analysis of 30 AGN using the 58 month light curves from Swift's Burst Alert Telescope (BAT) in the 14-150 keV band. PSDs were fit using a Monte Carlo based algorithm to take into account windowing effects and measurement error. All but one source were found to be fit very well using an unbroken power law with a slope of ~-1, consistent at low frequencies with previous studies in the 2-10 keV band, with no evidence of a break in the PSD. For 5 of the highest S/N sources we tested the energy dependence of the PSD and found no significant difference in the PSD at different energies. Unlike previous studies of X-ray variability in AGN, we do not find any significant correlations between the hard X-ray variability and different properties of the AGN including luminosity and black hole mass. The lack of break frequencies and correlations seem to indicate that AGN are similar to the high state of Galactic Black Holes.

FAST TIMING ANALYSIS OF CYGNUS X-1 USING THE SPECTROMETER ON THEINTERNATIONAL GAMMA-RAY ASTROPHYSICS LABORATORY

For the very first time, we report the high frequency analysis of Cyg X-1 up to hard X-ray using SPI on-board INTEGRAL. After analyzing the possible contribution fromthe background, and using INTEGRAL archive from March 2005 to May 2008, Power Density Spectra (PDS) were obtained until 130 keV. First, we show that their overall shape is very similar to that observed at lower energies, as they are well described by sets of Lorentzians. The strength of this fast variability (up to 40 Hz) does not drop at high energy since we show that it remains at ~25% rms, even in the highest energy bands. Second, the hard X-ray variability patterns of Cyg X-1 are state dependent: the softer the spectrum (or the lower the hardness ratio), the lower the total fractional variability and the higher the typical frequencies observed. The strength of total variability as a function of energy and state is then investigated. By comparison with simultaneous and published RXTE/PCA data, we showed that in the hard state, it remains quite constant in the 2-130 keV energy range. In our softer state, it is also flat until 50 keV and may increase at higher energy. The implications of this behavior on the models are then discussed.