Source author record

Ranjeev Misra

Ranjeev Misra appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

53works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

53 published item(s)

preprint2024arXiv

Black hole spin measurements in LMC X-1 and Cyg X-1 are highly model-dependent

The black-hole spin parameter, $a_*$, was measured to be close to its maximum value of 1 in many accreting X-ray binaries. In particular, $a_*\gtrsim 0.9$ was found in a number of studies of LMC X-1. These measurements were claimed to take into account both statistical and systematic uncertainties. We perform new measurements using a recent simultaneous observation by NICER and NuSTAR, providing a data set of high quality. We use the disk continuum method together with improved models for coronal Comptonization. With the standard relativistic disk model and optically thin Comptonization, we obtain values of $a_*$ similar to those obtained before. We then consider modifications to the standard model. Using a color correction of 2, we find $a_*\approx 0.64$--0.84. We then consider disks with dissipation in surface layers. To account for that, we assume the standard disk is covered by a warm and optically thick Comptonizing layer. Our model with the lowest $χ^2$ yields then $a_*\approx 0.40^{+0.41}_{-0.32}$. In order to test the presence of such effects in other sources, we also study an X-ray observation of Cyg X-1 by Suzaku in the soft state. We confirm the previous findings of $a_*>0.99$ using the standard model, but then we find a weakly constrained $a_*\approx 0.82^{+0.16}_{-0.74}$ when including an optically thick Comptonizing layer. We conclude that determinations of the spin using the continuum method can be highly sensitive to the assumptions about the disk structure.

preprint2022arXiv

AstroSat Detection of a Quasi-periodic Oscillation at ~ 42 Hz in Cygnus X-2

We report the results of AstroSat observations of Cygnus X-2 during February 2016. The source's power density spectrum generated using LAXPC data revealed the presence of a prominent Quasi-periodic Oscillation (QPO) at $\sim42$ Hz with broadband continuum noise at lower frequencies at $\sim10$ Hz. The large effective area of LAXPC at $\gtrsim$30 keV allowed for an unprecedented study of the energy dependence of the QPO and the broad noise continuum. The fractional r.m.s increases with energy, and its shape is similar for both the QPO and the continuum noise, suggesting a common radiative origin. However, while the QPO exhibits hard time lags, with the high energy photons lagging the low ones by a few milliseconds, the continuum noise shows the opposite behavior. The photon spectrum from SXT and LAXPC in $0.7-30$ keV band comprises the soft component from a disc and a hard Comptonized component from a hot corona. While the energy dependence of the r.m.s shows that the QPO and the continuum noise variability are dominated by the Comptonized component, the change in sign of the time-lag suggests that the dynamic origin of the QPO may be in the disk while the noise continuum may originate from the corona.

preprint2022arXiv

Broadband Spectral properties of MAXI J1348--630 using AstroSat Observations

We present broadband X-ray spectral analysis of the black hole X-ray binary MAXI J1348--630, performed using five AstroSat observations. The source was in the soft spectral state for the first three and in the hard state for the last two observations. The three soft state spectra were modelled using a relativistic thin accretion disc with reflection features and thermal Comptonization. Joint fitting of the soft state spectra constrained the spin parameter of the black hole $a_*$ $>$ 0.97 and the disc inclination angle $i$ = 32.9$\substack{+4.1 \\ -0.6}$ degrees. The bright and faint hard states had bolometric flux a factor of $\sim 6$ and $\sim 10$ less than that of the soft state. Their spectra were fitted using the same model except that the inner disc radius was not assumed to be at the last stable orbit. However, the estimated values do not indicate large truncation radii and the inferred accretion rate in the disc was an order of magnitude lower than that of the soft state. Along with earlier reported temporal analysis, AstroSat data provides a comprehensive picture of the evolution of the source.

preprint2022arXiv

Correlations between X-ray spectral parameters of Mkn 421 using long-term Swift-XRT data

We have performed a detailed analysis of the X-ray spectra of the blazar Mkn 421 using Swift-XRT observations taken between 2005 and 2020, to quantify the correlations between spectral parameters for different models. In an earlier work, it has been shown that such spectral parameter correlations obtained from a single short flare of duration $\sim$ 5-days of Mkn 421, can be used to distinguish spectrally degenerate models and provide estimates of physical quantities. In this work, we show that the results from the long-term spectral parameter correlations are consistent with those obtained from the single flare. In particular, that the observed spectral curvature is due to maximum cutoff energy in the particle distribution is ruled out. Instead, models where the curvature is due to the energy dependence of escape or acceleration time-scale of the particles are favored. The estimated values of the physical parameters for these models are similar to the ones obtained from the single flare analysis and are somewhat incompatible with the physical assumption of the models, suggesting that more complex physical models are required. The consistency of the results obtained from the long and short-term evolution of the source, underlines the reliability of the technique to use spectral parameter correlations to distinguish physical models.

preprint2022arXiv

Detection of low-frequency breaks in power density spectrum of GX 339-4 in faint low/hard state observations using AstroSat data

We present the spectro-timing analysis of the black hole binary system GX 339-4 using AstroSat data extracted at the beginning of outbursts in 2017 and 2019. The joint spectral fitting of LAXPC and SXT data revealed that the source was in a faint low/hard state for both 2017 and 2019 observations with a nearly equal photon index of $\sim$1.57 and $\sim$1.58 and with Eddington ratio, $L/L_{Edd}$, of 0.0011 and 0.0046 respectively. The addition of a reflection component into the spectral modelling improved the fit ($Δχ^2$ $\approx~6$ for 2017 and $Δχ^2$ $\approx~7$ for 2019), pointing towards the presence of weak reflection features arising due to irradiation of the accretion disk. The power density spectrum (PDS) consisted of strong band-limited noise with a break at low frequencies, described with a combination of few zero-centered Lorentzian. The fitting revealed a low-frequency break at $\sim$6~mHz for 2017 and at $\sim$11~mHz for 2019 observation, whose detection is validated by results from independent detectors (LAXPCs and SXT). The break frequency is roughly consistent with results obtained from earlier observations that showed an evolution of the frequency with flux, which is in accordance with the truncated disk model. Associating the break frequency with viscous time scale of the accretion disk, we estimated a truncation radius of $\sim$93 and $\sim$61 gravitational radius for 2017 and 2019 observation, respectively.

preprint2022arXiv

On the energy dependence of the QPO phenomenon in the black hole system MAXI J1535-571

Previous analysis of AstroSat observations of the black hole system MAXI J1535-571, have revealed the presence of a strong Quasi-Periodic Oscillation (QPO) whose frequency is correlated with the high energy spectral index. Here, we fit the spectra as emitted from a truncated disc with an inner hot corona, study the QPO frequency dependence on other spectral parameters and model the energy dependent r.m.s and time-lag of the QPO to identify the physical spectral parameters whose variation are responsible for the QPO. The QPO frequency is found to also correlate with the scattering fraction (i.e. the fraction of the soft photons Comptonized) and its dependence on the accretion rate and inner disc radii is consistent with it being the dynamical frequency. The time-lag between the hard and soft photons is negative for QPO frequency > 2.2 Hz and is positive for lesser values, making this the second black hole system to show this behaviour after GRS 1915+105. Modelling the energy dependent time-lag and r.m.s requires correlated variation of the accretion rate, inner disc radii and the coronal heating rate, with the latter having a time-lag compared to the other two for QPO frequencies less than < 2.2 Hz and which changes sign (i.e. the coronal heating variation precedes the accretion rate one) for higher values. The implications of the results are discussed.

preprint2022arXiv

Probing the shot behaviour in Cygnus X-1 using simultaneous AstroSat-NICER observation

We analyse the aperiodic flaring features, also known as shots, observed in Cyg X-1 in the 0.1-80 keV energy band using a 6.39 ks simultaneous observation with AstroSat and NICER. We detect 49 simultaneous shots in the soft and hard X-ray bands with NICER and AstroSat-LAXPC, respectively. We observe the shot profile for the first time in soft X-rays (0.1-3 keV), which shows a spectral peak at $\sim$2 keV. Using time-averaged spectroscopy, we measured the truncation of the inner accretion disk at $6.7\pm0.2$ gravitational radii. The shot-phase resolved spectroscopy allowed us to identify the origin of some of the brightest aperiodic peaks in the soft X-rays. We find that the accretion rate is consistent with a constant during the shots while the inner edge of the accretion disk moves inwards/outwards as these shots rise/decay. We discuss the possible mechanisms causing the swing in the inner radius.

preprint2022arXiv

Spectral properties of soft X-ray transient MAXI J0637$-$430 using AstroSat

Soft X-ray transients are systems that are detected when they go into an outburst, wherein their X-ray luminosity increases several orders of magnitude. These outbursts are markers of the poorly understood change in the spectral state of these systems from low/hard state to high/soft state. We report the spectral properties of one such soft X-ray transient: MAXI J0637$-$430, with data from the \textit{SXT} and \textit{LAXPC} instruments on-board \textit{AstroSat} mission. The source was observed for a total of $\sim$ 60 ks over two observations on 8$^{th}$ and 21$^{st}$ November, 2019 soon after its discovery. Flux resolved spectral analysis of the source indicates the presence of a multi-colour blackbody component arising from the accretion disk and a thermal Comptonization component. The stable low temperature ($\sim$ 0.55 $keV$) of the blackbody component, points to a cool accretion disk with an inner disk radius of the order of a few hundred $km$. In addition, we report the presence of a relativistically broadened Gaussian line at 6.4 $keV$. The disk dominated flux and photon power law index of $\gtrapprox 2$ and a constant inner disk radius indicate the source to be in the soft state. From the study we conclude that MAXI J0637$-$430 is a strong black hole X-ray binary candidate.

preprint2022arXiv

Time-resolved spectroscopy on the heartbeat state of GRS 1915+105 using AstroSat

AstroSat spectra of the black hole system GRS 1915+105 during the heartbeat state (with a varying oscillation period from 150 to 100 secs) were analysed using a truncated relativistic disc model along with a Comptonization component. Spectra were fitted for segments of length ~ 24 secs. The oscillation can be described as coordinated variations of the accretion rate, Comptonised flux, and the inner disc radius, with the latter ranging from 1.235-5 gravitational radii. Comparison with results from the $χ$ and Intermediate states shows that while the accretion rate and the high energy photon index were similar, the inner disc radius and the fraction of Comptonised photons were larger for these states than for the heartbeat one. The coronal efficiency $η\equiv L_{ac}/\dot M c^2$, where $L_{ac}$ is the radiative luminosity generated in the corona is found to be approximately $\propto \dot M^{-2/3}$ for all the observations. The efficiency decreases with inner radii for the heartbeat state but has similar values for the $χ$ and Intermediate states where the inner radii is larger. The implications of these results are discussed.

preprint2021arXiv

AstroSat/LAXPC detection of milli-second phenomena in 4U 1728-34

The low mass X-ray binary 4U 1728-24 was observed with AstroSat/LAXPC on 8th March 2016. Data from a randomly chosen one orbit of over 3 ks was analyzed for detection of rapid intensity variations. We found that the source intensity was nearly steady but towards the end of the observation a typical Type-1 burst was detected. Dynamical power spectrum of the data in the 3 to 20 keV band, reveals presence of a kHz Quasi-Periodic Oscillation (QPO) whose frequency drifted from around 815 Hz at the beginning of the observation to about 850 Hz just before the burst. The QPO is also detected in the 10 to 20 keV band, which was not obtainable by earlier RXTE observations of this source. Even for such a short observation with a drifting QPO frequency, the time-lag between the 5 to 10 and 10 to 20 keV bands can be constrained to be less than 100 microseconds. The Type-1 burst that lasted for about 20 secs had a typical profile. During the first four seconds dynamic power spectra reveal a burst oscillation whose frequency increased from 361.5 to 363.5 Hz. This is consistent with the earlier results obtained with RXTE/PCA, showing the same spin frequency of the neutron star. The present results demonstrate the capability of LAXPC instrument for detecting millisecond variability even from short observations. After RXTE ceased operation, LAXPC on AstroSat is the only instrument at present with capability of detecting kHz QPOs and other kind of rapid variations from 3 keV to 20 keV and possibly at higher energies also.

preprint2021arXiv

Testing evolution of LFQPOs with mass accretion rate in GRS 1915+105 with Insight-HXMT

Using the Insight-HXMT observations of GRS 1915+105 when it exhibits low frequency quasi-periodic oscillations (QPOs), we measure the evolution of the QPO frequency along with disk inner radius and mass accretion rate. We find a tight positive correlation between the QPO frequency and mass accretion rate. Our results extend the finding of previous work with AstroSat to a larger range of accretion rate with independent instruments and observations. Treating the QPO frequency of GRS 1915+105 as the relativistic dynamic frequency of a truncated disk, we are able to confirm the high spin nature of the black hole in GRS 1915+105. We also address the potential of our finding to test general relativity in the future.

preprint2020arXiv

Broadband X-ray Spectral Study of Ultra-luminous X-ray Source M81 X--6

Ultra-luminous X-ray sources (ULXs) are a class of extra-galactic, point-like, off-nuclear X-ray sources with X-ray luminosity from $\sim 10^{39}$~erg s$^{-1}$ to $10^{41}$~erg s$^{-1}$. We investigated the temporal and broadband X-ray spectral properties of the ULX M81~X--6 using simultaneous \textit{Suzaku} and \textit{NuSTAR} observations. To understand the nature of the source, we searched for pulsating signals from the source using the \textit{NuSTAR} observation. However, we failed to identify any strong pulsating signals from the source. Alternatively, the broadband spectral modelling with accreting magnetic neutron star continuum model provides a statistically acceptable fit, and the inferred spectral parameters and X-ray colours are consistent with other pulsating ULXs. Thus, our analysis suggests that M81~X--6 is another candidate ULX pulsar.

preprint2020arXiv

Correlating nonlinear time series and spectral properties of IGR J17091-3624: Is it similar to GRS 1915+105?

We explore the nonlinear properties of IGR J17091-3624 in the line of the underlying behaviour for GRS 1915+105, following correlation integral method. We find that while the latter is known to reveal the combination of fractal (or even chaotic) and stochastic behaviours depending on its temporal class, the former mostly shows stochastic behaviour. Therefore, although several observations argue IGR J17091-3624 to be similar to GRS 1915+105 with common temporal classes between them, underlying nonlinear time series analyses offer differently. Nevertheless, the Poisson noise to $rms$ variation value for IGR J17091-3624 turns out to be high, arguing them to be Poisson noise dominated and hence may plausibly lead to suppression of the nonlinear character, if any. Indeed, it is a very faint source compared to GRS 1915+105. However, by increasing time bin some of the temporal classes of IGR J17091-3624 show deviation from stochasticity, indicating plausibility of higher fractal dimension. Along with spectral analysis, overall IGR J17091-3624 argues to reveal three different accretion classes: slim, Keplerian and advective accretion discs.

preprint2020arXiv

Correlation between relativistic reflection fraction and photon index in NuSTAR sample of Seyfert 1 AGN

The primary X-ray emission from AGN, described by a power law, irradiates the accretion disc producing reflection features in the spectrum. The reflection features arising from the inner regions of the disc can be significantly modified by the relativistic effects near the black hole. We investigate the relationship between the relativistic reflection fraction $R_f$, defined as the ratio of the coronal intensity that illuminates the accretion disc to the coronal intensity observed directly, and the hard X-ray photon index $Γ$ of a NuSTAR sample of Seyfert 1 galaxies. The X-ray spectra are modelled using relxill code which helps to directly obtain the reflection fraction of a relativistically smeared reflection component. The parameter $R_f$ depends on the amount of Comptonised X-ray emission intercepted by the inner accretion disc. We found a positive correlation between $Γ$ and $R_f$ in our sample. Seed photons from a larger area of an accretion disc entering the corona will result in increased cooling of the coronal plasma, giving rise to steeper X-ray spectrum. The corona irradiating the larger area of the disc will result in higher reflection fraction. Thus, the observed $R_f-Γ$ relation is most likely related to the variations in the disc-corona geometry of AGN.

preprint2020arXiv

Early Warning Signals Indicate a Critical Transition in Betelgeuse

Critical transitions occur in complex dynamical systems, when the system dynamics undergoes a regime shift. These can often occur with little change in the mean amplitude of system response prior to the actual time of transition. The recent dimming and brightening event in Betelgeuse occured as a sudden shift in the brightness and has been the subject of much debate. Internal changes or an external dust cloud have been suggested as reasons for this change in variability. We examine whether the dimming and brightening event of 2019-20 could be due to a critical transition in the pulsation dynamics of Betelgeuse, by studying the characteristics of the light curve prior to transition. We calculate the quantifiers hypothesised to rise prior to a critical transition for the light curve of Betelgeuse upto the dimming event of 2019-20. These include the autocorrelation at lag-1, variance and the spectral coefficient calculated from detrended fluctation analysis (DFA), apart from two measures that quantify the recurrence properties of the light curve. Significant rises are confirmed using the Mann-Kendall trend test. We see a significant increase in all quantifiers (p < 0.05) prior to the dimming event of 2019-20. This suggests that the event was a critical transition related to the underlying nonlinear dynamics of the star. Together with results that suggests minimal change in $T_{eff}$ and infra-red flux, a critical transition in the pulsation dynamics could be a possible reason for the unprecedented dimming of Betelgeuse. The rise in the studied quantifiers prior to the dimming event, supports this possibility.

preprint2020arXiv

Gamma-ray Emission Properties of Four Bright Fermi-LAT AGNs: Implications on Emission Processes

The X-ray, Ultraviolet, Optical emission from radio-quiet AGNs, black hole binaries, and other compact sources, in general, follow a lognormal flux distribution, a linear rms-flux relation, and a (broken) power-law power spectral densities (PSDs). These characteristics are normally attributed to the multiplicative combination of fluctuations in the accretion disk. Similar features have been inferred for some well-observed blazars in different energy bands, but a systematic study over a long duration is still missing. Using a continuous gamma-ray light curves over 3-days cadence from August 2008 - October 2015, we present the first systematic study of these features in four sources: the FR I radio galaxy NGC 1275 and three blazars- Mrk 421, B2 1520+31 and PKS 1510-089. For all, except Mrk 421, the flux spans $\gtrsim$ 2 orders of magnitude. For blazars, a log-normal profile describes the flux histograms better compared to a Gaussian, while none is favored for NGC 1275, the only non-blazar source, suggesting a complex distribution. Regardless of flux histogram profile, the rms-flux relation is linear for all with PSDs being consistent with a power-law shot noise spectrum despite hints of breaks. The inferred results are consistent with the properties of unresolved magnetic reconnection sites, as inferred in the X-ray emission from the whole Solar disk and the statistical characteristics of magnetic reconnection based minijets-in-a-jet model. The results, thus, suggest a strong jet-accretion-disk coupling with energy input from the central source being distributed over a wide range in time and energy by the reconnection process depending on the geometry and local physical conditions.

preprint2020arXiv

Identification of QPO frequency of GRS 1915+105 as the relativistic dynamic frequency of a truncated accretion disk

We have analyzed AstroSat observations of the galactic micro-quasar system GRS 1915+105, when the system exhibited C-type Quasi-periodic Oscillations (QPOs) in the frequency range of 3.4-5.4 Hz. The broad band spectra (1-50 keV) obtained from simultaneous LAXPC and SXT can be well described by a dominant relativistic truncated accretion disk along with thermal Comptonization and reflection. We find that while the QPO frequency depends on the inner radii with a large scatter, a much tighter correlation is obtained when both the inner radii and accretion rate of the disk are taken into account. In fact, the frequency varies just as the dynamic frequency (i.e. the inverse of the sound crossing time) as predicted decades ago by the relativistic standard accretion disk theory for a black hole with spin parameter of ~0.9. We show that this identification has been possible due to the simultaneous broad band spectral coverage with temporal information as obtained from AstroSat.

preprint2020arXiv

NuSTAR observation of Ark 564 reveals the variation of coronal temperature with flux

The hard X-ray spectral index of some AGN has been observed to steepen with the source flux. This has been interpreted in a Comptonization scenario, where an increase in the soft flux decreases the temperature of the corona, leading to steepening of the photon index. However, the variation of the coronal temperature with flux has been difficult to measure due to the presence of complex reflection component in the hard X-rays and the lack of high-quality data at that energy band. Recently, a 200 ks {\it NuSTAR} observation of Ark 564 in 3--50 keV band, revealed the presence of one of the coolest coronae with temperature $kT_e \sim 15$ keV in the time-averaged spectrum. Here, we re-analyse the data and examined the spectra in four flux levels. Our analysis shows that the coronal temperature decreased from $\sim 17$ to $\sim 14$ keV as the flux increased. The high energy photon index $Γ\sim 2.3$ varied by less than $0.1$, implying that the optical depth of the corona increased by about 10\% as the flux increased. This first reporting of coronal temperature variation with flux shows that further long observation by {\it NuSTAR} of this and other sources would shed light on the geometry and dynamics of the inner regions of the accretion flow.

preprint2020arXiv

On the determination of log-normal flux distributions for astrophysical systems

Determining whether the flux distribution of an Astrophysical source is a Gaussian or a log-normal, provides key insight into the nature of its variability. For lightcurves of moderate length ($< 10^3$), a useful first analysis is to test the Gaussianity of the flux and logarithm of the flux, by estimating the skewness and applying the Anderson-Darling (AD) method. We perform extensive simulations of lightcurves with different lengths, variability, Gaussian measurement errors and power spectrum index $β$ (i.e. $P(f) \propto f^{-β}$), to provide a prescriptionand guidelines for reliable use of these two tests. We present empirical fits for the expected standard deviation of skewness and tabulated AD test critical values for $β= 0.5$ and $1.0$, which differ from the values given in the literature which are for white noise ($β= 0$). Moreover, we show that for white noise, for most practical situations, these tests are meaningless, since binning in time alters the flux distribution. For $β\gtrsim 1.5$, the skewness variance does not decrease with length and hence the tests are not reliable. Thus, such tests can be applied only to systems with $β\gtrsim 0.5$ and $β\lesssim 1.0$. As an example of the prescription given in this work, we reconfirm that the Fermi data of the blazar, 3FGL\,J0730.2-1141, shows that its $γ$ ray flux is consistent with a log-normal distribution and not with a Gaussian one.

preprint2020arXiv

Unveiling the temporal properties of MAXI J1820+070 through AstroSat observations

We present here the results of the first broadband simultaneous spectral and temporal studies of the newly detected black hole binary MAXI J1820+070 as seen by SXT and LAXPC on-board AstroSat.The observed combined spectra in the energy range 0.7-80 keV were well modeled using disk black-body emission, thermal Comptonization and a reflection component. The spectral analysis revealed that the source was in its hard spectral state (Gamma=1.61) with a cool disk (kTin=0.22 keV). We report the energy dependent time-lag and root mean squared (rms) variability at different frequencies in the energy range 3-80 keV using LAXPC data. We also modeled the flux variability using a single zone stochastic propagation model to quantify the observed energy dependence of time-lag and fractional rms variability and then compared the results with that of Cygnus X-1. Additionally, we confirm the detection of a quasi-periodic oscillation with the centroid frequency at 47.7 mHz.

preprint2019arXiv

AstroSat observation of GX 5-1: Spectral and timing evolution

We report on the first analysis of AstroSat observation of the Z-source GX 5- 1 on February 26-27, 2017. The hardness-intensity plot reveals that the source traced out the horizontal and normal branches. The 0.8-20 keV spectra from simultaneous SXT and LAXPC data at different locations of the hardness-intensity plot can be well described by a disk emission and a thermal Comptonized component. The ratio of the disk flux to the total i.e. the disk flux ratio increases monotonically along the horizontal to the normal one. Thus, the difference between the normal and horizontal branches is that in the normal branch, the disk dominates the flux while in the horizontal one it is the Comptonized component which dominates. The disk flux scales with the inner disk temperature as T_{in}^{5.5} and not as T_{in}{4} suggesting that either the inner radii changes dramatically or that the disk is irradiated by the thermal component changing its hardness factor. The power spectra reveal a Quasi Periodic Oscillation whose frequency changes from \sim 30 Hz to 50 Hz. The frequency is found to correlate well with the disk flux ratio. In the 3-20 keV LAXPC band the r.m.s of the QPO increases with energy (r.m.s \prop E0.8), while the harder X-ray seems to lag the soft ones with a time-delay of a milliseconds. The results suggest that the spectral properties of the source are characterized by the disk flux ratio and that the QPO has its origin in the corona producing the thermal Comptonized component.

preprint2019arXiv

Study of long-term flux and photon Index distributions of blazars using RXTE observations

We present a detailed study of flux and index distributions of three blazars (one FSRQ and two BL Lacs) by using 16 years of Rossi X-ray Timing Explorer archival data. The three blazars were chosen such that their flux and index distributions have sufficient number of data points ($\geq$90) with relatively less uncertainty ($\bar{σ_{err}^{2}}/σ^{2}$<0.2) in light curves. Anderson-Darling (AD) test and histogram fitting shows that flux distribution of FSRQ 3C273 is log-normal, while its photon index distribution is Gaussian. This result is consistent with linear Gaussian perturbation in the particle acceleration time-scale, which produces log-normal distribution in flux. However, for two BL Lacs viz. Mkn501 and Mkn421, AD test shows that their flux distributions are neither Gaussian nor log-normal, and their index distributions are non-normal. The histogram fitting of Mkn501 and Mkn421, suggests that their flux distributions are more likely to be a bi-modal, and their index distributions are double Gaussian. Since, Sinha et al. (2018) had shown that Gaussian distribution of index produces a log-normal distribution in flux, double Gaussian distribution of index in Mkn501 and Mkn421 indicates that their flux distributions are probably double log-normal. Observation of double log-normal flux distribution with double Gaussian distribution in index reaffirms two flux states hypothesis. Further, the difference observed in the flux distribution of FSRQ (3C273) and BL Lacs (Mkn501 and Mkn421) at X-rays, suggest that the low energy emitting electrons have a single log-normal flux distribution while the high energy ones have a double log-normal flux distribution.

preprint2016arXiv

A model for the energy-dependent time-lag and rms of the heartbeat oscillations in GRS 1915+105

Energy dependent phase lags reveal crucial information about the causal relation between various spectral components and about the nature of the accretion geometry around the compact objects. The time-lag and the fractional root mean square (rms) spectra of GRS 1915+105 in its heartbeat oscillation class/$ρ$ state show peculiar behaviour at the fundamental and harmonic frequencies where the lags at the fundamental show a turn around at $\sim$ 10 keV while the lags at the harmonic do not show any turn around at least till $\sim$ 20 keV. The magnitude of lags are of the order of few seconds and hence cannot be attributed to the light travel time effects or Comptonization delays. The continuum X-ray spectra can roughly be described by a disk blackbody and a hard X-ray power-law component and from phase resolved spectroscopy it has been shown that the inner disk radius varies during the oscillation. Here, we propose that there is a delayed response of the inner disk radius (DROID) to the accretion rate such that $r_{in}(t)\propto \dot{m}^β(t-τ_d)$. The fluctuating accretion rate drives the oscillations of the inner radius after a time delay $τ_d$ while the power-law component responds immediately. We show that in such a scenario a pure sinusoidal oscillation of the accretion rate can explain not only the shape and magnitude of energy dependent rms and time-lag spectra at the fundamental but also the next harmonic with just four free parameters.

preprint2016arXiv

A Monte Carlo estimate for the fraction of thermal Comptonized photons that impinge back on the soft source in neutron star LMXBs

In earlier works, it was shown that the energy dependent soft time lags observed in kHz QPOs of neutron star low mass X-ray binaries (LMXBs) can be explained as being due to Comptonization lags provided a significant fraction ($η\sim 0.2$ - $0.8$) of the Comptonized photons impinge back into the soft photon source. Here we use a Monte Carlo scheme to verify if such a fraction is viable or not. In particular we consider three different Comptonizing medium geometries: (i) a spherical shell, (ii) a boundary layer like torus and (iii) a corona on top of an accretion disk. Two set of spectral parameters corresponding to the 'hot' and 'cold' seed photon models were explored. The general result of the study is that for a wide range of sizes, the fraction lies within $η\sim 0.3$ - $0.7$, and hence compatible with the range required to explain the soft time lags. Since there is a large uncertainty in the range, we cannot concretely rule out any of the geometries or spectral models, but the analysis suggests that a boundary layer type geometry with a 'cold' seed spectral model is favoured over an accretion corona model. Better quality data will allow one to constrain the geometry more rigorously. Our results emphasise that there is significant heating of the soft photon source by the Comptonized photons and hence this effect needs to be taken into account for any detailed study of these sources.

preprint2016arXiv

AstroSat/LAXPC reveals the high energy variability of GRS 1915+105 in the chi class

We present the first quick look analysis of data from nine {\it AstroSat}'s LAXPC observations of GRS 1915+105 during March 2016 when the source had the characteristics of being in Radio-quiet $χ$ class. We find that a simple empirical model of a disk blackbody emission, with Comptonization and a broad Gaussian Iron line can fit the time averaged 3--80 keV spectrum with a systematic uncertainty of 1.5\% and a background flux uncertainty of 4\%. A simple deadtime-corrected Poisson noise level spectrum matches well with the observed high frequency power spectra till 50 kHz and as expected the data show no significant high frequency ($> 20$ Hz) features. Energy dependent power spectra reveal a strong low frequency (2 - 8 Hz) Quasi-periodic oscillation (LFQPO) and its harmonic along with broad band noise. The QPO frequency changes rapidly with flux (nearly 4 Hz in ~ 5 hours). With increasing QPO frequency, an excess noise component appears significantly in the high energy regime (> 8 keV). At the QPO frequencies, the time-lag as a function of energy has a non-monotonic behavior such that the lags decrease with energy till about 15 -20 keV and then increase for higher energies. These first look results benchmark the performance of LAXPC at high energies and confirms that its data can be used for more sophisticated analysis such as flux or frequency-resolved spectro-timing studies.

preprint2016arXiv

Broadband X-ray emission and the reality of the broad iron line from the Neutron Star - White Dwarf X-ray binary 4U 1820-30

Broad relativistic iron lines from neutron star X-ray binaries are important probes of the inner accretion disk. The X-ray reflection features can be weakened due to strong magnetic fields or very low iron abundances such as is possible in X-ray binaries with low mass, first generation stars as companions. Here we investigate the reality of the broad iron line detected earlier from the neutron star low mass X-ray binary 4U~1820--30 with a degenerate helium dwarf companion. We perform a comprehensive, systematic broadband spectral study of the atoll source using \suzaku{} and simultaneous \nustar{} \& \swift{} observations. We have used different continuum models involving accretion disk emission, thermal blackbody and thermal Comptonization of either disk or blackbody photons. The \suzaku{} data show positive and negative residuals in the region of iron K band. These features are well described by two absorption edges at $7.67\pm0.14\kev$ and $6.93\pm0.07\kev$ or partial covering photoionized absorption or by blurred reflection. Though, the simultaneous \swift{} and \nustar{} data do not clearly reveal the emission or absorption features, the data are consistent with the presence of either absorption or emission features. Thus, the absorption based models provide an alternative to the broad iron line or reflection model. The absorption features may arise in winds from the inner accretion disk. The broadband spectra appear to disfavour continuum models in which the blackbody emission from the neutron star surface provides the seed photons for thermal Comptonization. Our results suggest emission from a thin accretion disk ($kT_{disk} \sim 1\kev$), Comptonization of disk photons in a boundary layer most likely covering a large fraction of the neutron star surface and innermost parts of the accretion disk, and blackbody emission ($kT_{bb} \sim 2\kev$) from the polar regions.

preprint2016arXiv

Constraining the size of the Comptonizing medium by modeling the energy dependent time-lags of kHz QPOs of Neutron star system

In earlier works, we had shown that the observed soft lags and r.m.s versus energy of the lower kHz QPO of neutron star binaries can be explained in the framework of a thermal Comptonization model. It was also shown that such an interpretation can provide estimates of the size and geometry of the Comptonizing medium. Here we study the dependence of these estimates on the time-averaged spectral model assumed and on the frequency of the QPO. We use the high quality time lag and r.m.s obtained during March 3rd 1996 observation of 4U 1608-52 by RXTE as well as other observations of the source at different QPO frequencies where a single time-lag between two broad energy bands have been reported. We compare the results obtained when assuming that the time-averaged spectra are represented by the spectrally degenerate "hot" and "cold" seed photon spectral models. We find that for the "hot" seed photon model the medium size is in the range of 0.3-2.0 kms and the size decreases with increasing QPO frequency. On the other hand for the "cold" seed photon model the range for the sizes are much larger 0.5-20 kms and hence perhaps show no variation with QPO frequency. Our results emphasis the need for broad band spectral information combined with high frequency timing to lift this degeneracy. We further show that the r.m.s as a function of energy for the upper kHz QPO is similar to the lower one and indeed we find that the driver for this QPO should be temperature variations of the corona identical to the lower kHz QPO. However, the time lag reported for the upper kHz QPO is hard, which if confirmed, would challenge the simple Comptonization model presented here. It would perhaps imply that reverberation lags are also important and/or the dominate spectral component is not a single temperature medium but a multi-temperature complex one.

preprint2016arXiv

Evidence for two Lognormal States in Multi-wavelength Flux Variation of FSRQ PKS 1510-089

We present a systematic characterization of multi-wavelength emission from blazar PKS 1510-089 using well-sampled data at infrared(IR)-optical, X-ray and $γ$-ray energies. The resulting flux distributions, except at X-rays, show two distinct lognormal profiles corresponding to a high and a low flux level. The dispersions exhibit energy dependent behavior except for the LAT $γ$-ray and optical B-band. During the low level flux states, it is higher towards the peak of the spectral energy distribution, with $γ$-ray being intrinsically more variable followed by IR and then optical, consistent with mainly being a result of varying bulk Lorentz factor. On the other hand, the dispersions during the high state are similar in all bands expect optical B-band, where thermal emission still dominates. The centers of distributions are a factor of $\sim 4$ apart, consistent with anticipation from studies of extragalactic $γ$-ray background with the high state showing a relatively harder mean spectral index compared to the low state.

preprint2016arXiv

Multi-wavelength Temporal Variability of the Blazar 3C 454.3 during 2014 Activity Phase

We present a multi-wavelength temporal analysis of the blazar 3C 454.3 during the high $γ$-ray active period from May-December, 2014. Except for X-rays, the period is well sampled at near-infrared (NIR)-optical by the \emph{SMARTS} facility and the source is detected continuously on daily timescale in the \emph{Fermi}-LAT $γ$-ray band. The source exhibits diverse levels of variability with many flaring/active states in the continuously sampled $γ$-ray light curve which are also reflected in the NIR-optical light curves and the sparsely sampled X-ray light curve by the \emph{Swift}-XRT. Multi-band correlation analysis of this continuous segment during different activity periods shows a change of state from no lags between IR and $γ$-ray, optical and $γ$-ray, and IR and optical to a state where $γ$-ray lags the IR/optical by $\sim$3 days. The results are consistent with the previous studies of the same during various $γ$-ray flaring and active episodes of the source. This consistency, in turn, suggests an extended localized emission region with almost similar conditions during various $γ$-ray activity states. On the other hand, the delay of $γ$-ray with respect to IR/optical and a trend similar to IR/optical in X-rays along with strong broadband correlations favor magnetic field related origin with X-ray and $γ$-ray being inverse Comptonized of IR/optical photons and external radiation field, respectively.

preprint2016arXiv

NuSTAR and Swift joint view of neutron star X-ray binary 4U~1728--34 : Disc reflection in the island and lower banana states

We analyze two simultaneous \nustar{} and \swift{} data of the Atoll type neutron star (NS) X-ray binary 4U~1728--34 observed on 1 and 3 October, 2013. Based on power density spectra, hardness ratio and colour-colour diagram, we infer that the first and second observations belong to the island state and lower banana state, respectively. During island state, a low-frequency quasi-periodic oscillation (LFQPO) at $\sim4.53$ Hz is observed along with four type-I X-ray bursts ($\sim 10-15$ s duration each). The X-ray luminosity of the source during the island and lower banana states are found to be $L_{X}= 1.1$ and $1.6\times10^{37}$ erg s$^{-1}$, respectively which correspond to $\sim 6\%$ and $\sim 9\%$ of the Eddington luminosity. The combined burst spectrum is well represented by a blackbody with a characteristic temperature of $2.22\pm0.05\kev$ and the blackbody radius of $9.87\pm1.86$ km which is consistent with the typical radius of the neutron star. The persistent energy spectra from both observations in the energy band $1-79\kev$ are well described with thermal emission from the NS surface ($kT_{bb}\simeq1-2.5\kev$), Comptonized emission of thermal seed photons from the hot boundary layer/corona and the strong reflection component from the accretion disc. We detect a broad Iron line in the $5-8 \kev$ band and reflection hump in the $15-30 \kev$ band. These features are well modelled by the {\tt relxill} reflection model. From joint spectral fitting, we constrain the inclination angle of the binary system and inner disc radius to be $22\degree-40\degree$ and $(2.8-4.3)\times R_{ISCO}$, respectively. From inner disc radius, we estimate the magnetic field to be $(3.3-6.5)\times10^{8}$ Gauss.

preprint2016arXiv

On the reality of broad iron L lines from the narrow line Seyfert 1 galaxies 1H0707-495 and IRAS 13224-3809

We performed time resolved spectroscopy of 1H0707-495 and IRAS 13224-3809 using long XMM-Newton observations. These are strongly variable narrow line Seyfert 1 galaxies and show broad features around 1 keV that has been interpreted as relativistically broad Fe L$α$ lines. Such features are not clearly observed in other AGN despite sometimes having high iron abundance required by the best fitted blurred reflection models. Given the importance of these lines, we explore the possibility if rapid variability of spectral parameters may introduce broad bumps/dips artificially in the time averaged spectrum, which may then be mistaken as broadened lines. We tested this hypothesis by performing time resolved spectroscopy using long (> 100 ks) XMM-Newton observations and by dividing it into segments with typical exposure of few ks. We extracted spectra from each such segment and modelled using a two component phenomenological model consisting of a power law to represent hard component and a black body to represent the soft emission. As expected both the sources showed variations in the spectral parameters. Using these variation trends, we simulated model spectra for each segment and then co-added to get a combined simulated spectrum. In the simulated spectra, we found no broad features below 1 keV and in particular no deviation near 0.9 keV as seen in the real average spectra. This implies that the broad Fe L? line that is seen in the spectra of these sources is not an artifact of the variation of spectral components and hence providing evidence that the line is indeed genuine.

preprint2016arXiv

UV and X-Ray Variability of the Narrow-Line Seyfert 1 Galaxy Ark 564

We analyse eight XMM-Newton observations of the bright Narrow-Line Seyfert 1 galaxy Arakelian 564 (Ark 564). These observations, separated by ~6 days, allow us to look for correlations between the simultaneous UV emission (from the Optical Monitor) with not only the X-ray flux but also with the different X-ray spectral parameters. The X-ray spectra from all the observations are found to be adequately fitted by a double Comptonization model where the soft excess and the hard X-ray power law are represented by thermal Comptonization in a low temperature plasma and hot corona, respectively. Apart from the fluxes of each component, the hard X-ray power law index is found to be variable. These results suggest that the variability is associated with changes in the geometry of the inner region. The UV emission is found to be variable and well correlated with the high energy index while the correlations with the fluxes of each component are found to be weaker. Using viscous time-scale arguments we rule out the possibility that the UV variation is due to fluctuating accretion rate in the outer disc. If the UV variation is driven by X-ray reprocessing, then our results indicate that the strength of the X-ray reprocessing depends more on the geometry of the X-ray producing inner region rather than on the X-ray luminosity alone.

preprint2016arXiv

X-ray/UV variability and the origin of soft X-ray excess emission from II Zw 177

We study X-ray and UV emission from the narrow-line Seyfert 1 galaxy II~Zw~177 using a $137\ks$ long and another $13\ks$ short \xmm{} observation performed in 2012 and 2001, respectively. Both observations show soft X-ray excess emission contributing $76.9\pm4.9\%$ in 2012 and $58.8\pm10.2\%$ in 2001 in the $0.3-2\kev$ band. We find that both blurred reflection from an ionized disc and Comptonized disc emission describe the observed soft excess well. Time-resolved spectroscopy on scales of $\sim20\ks$ reveals strong correlation between the soft excess and the powerlaw components. The fractional variability amplitude $F_{var}$ derived from EPIC-pn lightcurves at different energy bands is nearly constant ($F_{var} \sim20\%$). This is in contrast to other AGNs where the lack of short term variation in soft X-ray excess emission has been attributed to intense light bending in the framework of the "lamppost" model. Thus, the variations in powerlaw emission are most likely intrinsic to corona rather than just due to the changes of height of compact corona. The variable UV emission ($F_{var} \sim 1\%$) is uncorrelated to any of the X-ray components on short timescales suggesting that the UV emission is not dominated by the reprocessed emission. The gradual observed decline in the UV emission in 2012 may be related to the secular decline due to the changes in the accretion rate. In this case, the short term X-ray variability is not due to the changes in the seed photons but intrinsic to the hot corona.

preprint2015arXiv

Autoencoding Time Series for Visualisation

We present an algorithm for the visualisation of time series. To that end we employ echo state networks to convert time series into a suitable vector representation which is capable of capturing the latent dynamics of the time series. Subsequently, the obtained vector representations are put through an autoencoder and the visualisation is constructed using the activations of the bottleneck. The crux of the work lies with defining an objective function that quantifies the reconstruction error of these representations in a principled manner. We demonstrate the method on synthetic and real data.

preprint2015arXiv

Constraining distance and inclination angle of V4641 Sgr using ${\it SWIFT}$ and ${\it NuSTAR}$ observations during low soft spectral state

We present results from ${\it NuSTAR}$ and ${\it SWIFT}$/XRT joint spectral analysis of V4641 Sgr during a disk dominated or soft state as well as a powerlaw dominated or hard state. The soft state spectrum is well modeled by a relativistically blurred disk emission, a powerlaw, a broad Iron line, two narrow emission lines and two edges. The Markov Chain Monte Carlo simulation technique and the relativistic effects seen in the disk and broad Iron line allow us to self-consistently constrain the inner disk radius, disk inclination angle and distance to the source at 2.43$^{+0.39}_{-0.17}$ R$_g$ (GM/c$^2$), 69.5$^{+12.8}_{-4.2}$ degrees and 10.8$^{+1.6}_{-2.5}$ kpc respectively. For the hard state, the spectrum is a power-law with a weakly broad Iron line and an edge. The distance estimate gives a measure of the Eddington fraction, $L_{2.0 - 80.0 keV}/L_{Edd}$, to be $\sim$1.3 $\times$ 10$^{-2}$ and $\sim$1.9 $\times$ 10$^{-3}$ for the soft and hard states respectively. Unlike many other typical black hole systems which are always in a hard state at such low Eddington fraction, V4641 Sgr shows a soft, disk dominated state. The soft state spectrum shows narrow emission lines at $\sim 6.95$ and $\sim 8.31$ keV which can be identified as being due to emission from highly ionized Iron and Nickel in an X$-$ray irradiated wind respectively. If not due to instrumental effect or calibration error, this would be the first detection of a Ni fluorescent line in a black hole X$-$ray binary.

preprint2015arXiv

Investigating the connection between Quasi Periodic Oscillations and spectral components with NuSTAR data of GRS 1915+105

The low frequency quasi periodic oscillations (QPOs) are commonly observed during hard states of black hole binaries. Several studies have established various observational/empirical correlations between spectral parameters and QPO properties, indicating a close link between the two. However, the exact mechanism of generation of QPO is not yet well understood. In this paper, we present our attempts to comprehend the connection between the spectral components and the low frequency QPO observed in GRS 1915+105 using the data from NuSTAR. Detailed spectral modeling as well as the presence of the low frequency QPO and its energy dependence during this observation have been reported by Miller et al. (2013) and Zhang et al. (2015) respectively. We investigate the compatibility of the spectral model and energy dependence of the QPO by simulating light curves in various energy bands for small variation of the spectral parameters. The basic concept here is to establish connection, if any, between the QPO and the variation of either a spectral component or a specific parameter, which in turn can shed some light on the origin of the QPO. We begin with the best fit spectral model of Miller et al. (2013) and simulate the light curve by varying the spectral parameter at frequencies close to the observed QPO frequency in order to generate the simulated QPO. Further we simulate similar light curves in various energy bands in order to reproduce the observed energy dependence of RMS amplitude of the QPO. We find that the observed trend of increasing RMS amplitude with energy can be reproduced qualitatively if the spectral index is assumed to be varying with the phases of the QPO. Variation of any other spectral parameter does not reproduce the observed energy dependence.

preprint2015arXiv

Study of underlying particle spectrum during huge X-ray flare of Mkn 421 in April 2013

Context: In April 2013, the nearby (z=0.031) TeV blazar, Mkn 421, showed one of the largest flares in X-rays since the past decade. Aim: To study all multiwavelength data available during MJD 56392 to 56403, with special emphasis on X-ray data, and understand the underlying particle energy distribution. Methods: We study the correlations between the UV and gamma bands with the X-ray band using the z-transformed discrete correlation function. We model the underlying particle spectrum with a single population of electrons emitting synchrotron radiation, and do a statistical fitting of the simultaneous, time-resolved data from the Swift-XRT and the NuSTAR. Results: There was rapid flux variability in the X-ray band, with a minimum doubling timescale of $1.69 \pm 0.13$ hrs. There were no corresponding flares in UV and gamma bands. The variability in UV and gamma rays are relatively modest with $ \sim 8 \% $ and $\sim 16 \% $ respectively, and no significant correlation was found with the X-ray light curve. The observed X-ray spectrum shows clear curvature which can be fit by a log parabolic spectral form. This is best explained to originate from a log parabolic electron spectrum. However, a broken power law or a power law with an exponentially falling electron distribution cannot be ruled out either. Moreover, the excellent broadband spectrum from $0.3-79$ keV allows us to make predictions of the UV flux. We find that this prediction is compatible with the observed flux during the low state in X-rays. However, during the X-ray flares, the predicted flux is a factor of $2-50$ smaller than the observed one. This suggests that the X-ray flares are plausibly caused by a separate population which does not contribute significantly to the radiation at lower energies. Alternatively, the underlying particle spectrum can be much more complex than the ones explored in this work.

preprint2015arXiv

The Dust Content and Radiation Fields of Sample of Galaxies in the ELAIS-N1 Field

The Mid-IR colors ($F_{8}/F_{24}$) of galaxies together with their IR-UV luminosity correlations can be used to get some insight into the relative abundance of the different dust grain populations present in them. The ELAIS-N1 field contains thousands of galaxies which do not have optical spectra but have been observed in the Mid-IR by {\it Spitzer} and UV by {\it GALEX} making it ideal for these studies. As part of this work we have selected a sample of galaxies from the ELAIS-N1 field which have photometric observations in the MIR and UV as well as photometric redshifts from the SDSS database. We put the constraint that the redshifts are $\le$ 0.1, thereby giving us a total of 309 galaxies. We find that the majority of the galaxies in the sample are PAH dominated due to their high MIR flux ratio. We also find a reasonable correlation between the Mid-IR and the UV luminosities out of which the Mid-IR emission from PAHs at 8 $μ$m is marginally better correlated than the 24 $μ$m VSG emission with the UV luminosities. However, if we divide the sample based on their $F_{8}/F_{24}$ ratios which is also an indicator of metallicity, the MIR-UV correlation seems to increase with the $F_{8}/F_{24}$ ratio. But the MIR-UV correlations are not very different for the PAHs and the VSG population within the individual metallicity groups.

preprint2015arXiv

The Effect of X-ray Irradiation on the Time Dependent Behaviour of Accretion Disks with Stochastic Perturbations

The UV emission from X-ray binaries, is more likely to be produced by reprocessing of X-rays by the outer regions of an accretion disk. The structure of the outer disk may be altered due to the presence of X-ray irradiation and we discuss the physical regimes where this may occur and point out certain X-ray binaries where this effect may be important. The long term X-ray variability of these sources is believed to be due to stochastic fluctuations in the outer disk, which propagate inwards giving rise to accretion rate variation in the X-ray producing inner regions. The X-ray variability will induce structural variations in the outer disk which in turn may affect the inner accretion rate. To understand the qualitative behaviour of the disk in such a scenario, we adopt simplistic assumptions that the disk is fully ionised and is not warped. We develop and use a time dependent global hydrodynamical code to study the effect of a sinusoidal accretion rate perturbation introduced at a specific radius. The response of the disk, especially the inner accretion rate, to such perturbations at different radii and with different time periods is shown. While we didn't find any oscillatory or limit cycle behaviour, our results show that irradiation enhances the X-ray variability at time-scales corresponding to the viscous time-scales of the irradiated disk.

preprint2015arXiv

The influence of soft spectral components on the structure and stability of warm absorbers in AGN

The radiation from the central regions of active galactic nuclei, including that from the accretion disk surrounding the black hole, is likely to peak in the extreme ultraviolet $\sim 13 -100$ eV. However, due to Galactic absorption, we are limited to constrain the physical properties, i.e. the black hole mass and the accretion rate, from what observations we have below $\sim 10$ eV or above $\sim 100$ eV. In this paper we predict the thermal and ionization states of warm absorbers as a function of the shape of the unobservable continuum. In particular we model an accretion disk at $kT_{in} \sim 10$ eV and a {\it soft excess} at $kT_{se} \sim 150$ eV. The warm absorber, which is the highly ionized gas along the line of sight to the continuum, shows signatures in the $\sim 0.3 - 2$ keV energy range consisting of numerous absorption lines and edges of various ions, some of the prominent ones being H- and He-like oxygen, neon, magnesium and silicon. We find that the properties of the warm absorber are significantly influenced by the changes in the temperature of the accretion disk, as well as by the strength of the {\it soft excess}, as they affect the optical depth particularly for iron and oxygen. These trends may help develop a method of characterising the shape of the unobservable continuum and the occurrence of warm absorbers.

preprint2014arXiv

Energy Dependent Time Delays of kHz Oscillations due to Thermal Comptonization

We study the energy dependent photon variability from a thermal Comptonizing plasma that is oscillating at kHz frequencies. In particular, we solve the linearised time dependent Kompaneets equation and consider the oscillatory perturbation to be either in the soft photon source or in the heating rate of the plasma. For each case, we self consistently consider the energy balance of the plasma and the soft photon source. The model incorporates the possibility of a fraction of the Comptonized photons impinging back into the soft photon source. We find that when the oscillation is due to the soft photon source, the variation of the fractional root mean sqaure (r.m.s) is nearly constant with energy and the time-lags are hard. However, for the case when the oscillation is due to variation in the heating rate of the corona, and when a significant fraction of the photons impinge back into the soft photon source, the r.m.s increases with energy and the time lags are soft. As an example, we compare the results with the $\sim 850$ Hz oscillation observed on March 3, 1996 for 4U 1608-52 and show that both the observed soft time-lags as well as the r.m.s versus energy can be well described by such a model where the size of the Comptonizing plasma is $\sim 1$ km. Thus, modelling of the time lags as due to Comptonization delays, can provide tight constraints on the size and geometry of the system. Detailed analysis would require well constrained spectral parameters.

preprint2014arXiv

Estimation of the Extragalactic Background Light using TeV Observations of BL~Lacs

The very high energy (VHE) gamma ray spectral index of high energy peaked blazars correlates strongly with its corresponding redshift whereas no such correlation is observed in the X-ray or the GeV bands. We attribute this correlation to a result of photon-photon absorption of TeV photons with the extragalactic background light (EBL) and utilizing this, we compute the allowed flux range for the EBL, which is independent of previous estimates. The observed VHE spectrum of the sources in our sample can be well approximated by a power-law, and if the de-absorbed spectrum is also assumed to be a power law, then we show that the spectral shape of EBL will be $εn(ε) \sim k log(\fracε{ε_p}) $. We estimate the range of values for the parameters defining the EBL spectrum, $k$ and $ε_p$, such that the correlation of the intrinsic VHE spectrum with redshift is nullified. The estimated EBL depends only on the observed correlation and the assumption of a power law source spectrum. Specifically, it does not depend on the spectral modeling or radiative mechanism of the sources, nor does it depend on any theoretical shape of the EBL spectrum obtained through cosmological calculations. The estimated EBL spectrum is consistent with the upper and lower limits imposed by different observations. Moreover, it also agrees closely with the theoretical estimates obtained through cosmological evolution models.

preprint2013arXiv

Alternating lags of QPO harmonics : A Generic model and its application to the 67 millihertz QPO of GRS 1915+105

A generic model for alternating lags in QPO harmonics is presented where variations in the photon spectrum are caused by oscillations in two parameters that characterize the spectrum. It is further assumed that variations in one of the parameters is linearly driven by variations in the other after a time delay $t_d$. It is shown that alternating lags will be observed for a range of $t_d$ values. A phenomenological model based on this generic one is developed which can explain the amplitude and phase lag variation with energy of the fundamental and the next three harmonics of the 67 mHz QPO observed in GRS 1915+105. The phenomenological model also predicts the variation of the Bicoherence phase with energy, which can be checked by further analysis of the observational data.

preprint2013arXiv

Comparison of time/phase lags in the hard state and plateau state of GRS 1915+105

We investigate the complex behavior of energy- and frequency-dependent time/phase lags in the plateau state and the radio-quiet hard state of GRS 1915+105. In our timing analysis, we find that when the source is faint in the radio, QPOs are observed above 2 Hz and typically exhibit soft lags (soft photons lag hard photons), whereas QPOs in the radio-bright plateau state are found below 2.2 Hz and consistently show hard lags. The phase lag at the QPO frequency is strongly anti-correlated with the QPO frequency, changing sign at 2.2 Hz. However, the phase lag at the frequency of the first harmonic is positive and nearly independent of frequency at at ~0.172 rad, regardless of the radio emission. The lag-energy dependence at the first harmonic is also independent of radio flux. However, the lags at the QPO frequency are negative at all energies during the radio-quiet state, but lags at the QPO frequency during the plateau state are positive at all energies and show a 'reflection-type' evolution of the lag-energy spectra with respect to the radio-quiet state. The lag-energy dependence is roughly logarithmic, but there is some evidence for a break around 4-6 keV. Finally, the Fourier frequency-dependent phase lag spectra are fairly flat during the plateau state, but increase from negative to positive during the radio-quiet state. We discuss the implications of our results in the light of some generic models.

preprint2013arXiv

Interpreting the large amplitude X-ray variation of GRS 1915+105 and IGR J17091-3624 as modulations of an accretion disc

Using the flux resolved spectroscopy for the first time, we analyse the RXTE/PCA data of the black hole X-ray binaries GRS 1915+105 and IGR J17091-3624, when both sources show large amplitude, quasi-regular oscillations in 2.0-60.0 keV X-ray light curves. For different observations, we extract spectra during the peak (spectrally soft) and dip (spectrally hard) intervals of the oscillation, and find that their spectra are phenomenologically complex, requiring at least two distinct spectral components. Besides a thermal Comptonization component, we find that the disc emission is better modelled by an index-free multicolour disc blackbody component (p-free disc model) rather than that from a standard accretion disc. While the peak and dip spectra are complex, remarkably, their difference spectra constructed by treating dip spectra as the background spectra of the peak spectra, can be modelled as a single p-free disc component. Moreover, the variability at different time-scales and energy bands of the peak flux level is always greater than or equal to the variability of the dip flux level, which strengthens the possibility that the peak flux level may be due to an independent spectral component added to the dip one. We also verify it using joint spectral analysis of peak and dip spectra with a variable emission component. Our result substantiates the standard hypothesis that the oscillations are due to the limit cycle behaviour of an unstable radiation pressure dominated inner disc. However, in this interpretation, the flux variation of the unstable disc can be several order of magnitudes as expected from some theoretical simulations and need not be fine tuned to match the factor ~10 variation seen between the peak and dip levels.

preprint2013arXiv

Properties of unique hard X-ray dips observed from GRS 1915+105 and IGR J17091-3624 and their implications

We report a comprehensive study on spectral and timing properties of hard X-ray dips uniquely observed in some so-called variability classes of the micro-quasars GRS 1915+105 and IGR J17091-3624. These dips are characterized by a sudden decline in the 2.0-60.0 keV X-ray intensity by a factor of 4-12 simultaneous with the increase in hardness ratio by a factor of 2-4. Using 31 observations of GRS 1915+105 with RXTE/PCA, we show that different behavior are observed in different types of variability classes, and we find that a dichotomy is observed between classes with abrupt transitions vs those with smoother evolution. For example, both energy-lag spectra and frequency-lag spectra of hard X-ray dips in classes with abrupt transitions and shorter dip intervals show hard-lag (hard photons lag soft photons), while both lag spectra during hard dips in classes with smoother evolution and longer dip intervals show soft-lag. Both lag time-scales are of the order of 100-600 msec. We also show that timing and spectral properties of hard X-ray dips observed in light curves of IGR J17091-3624 during its 2011 outburst are consistent with the properties of the abrupt transitions in GRS 1915+105 rather than smooth evolution. A global correlation between the X-ray intensity cycle time and hard dip time is observed for both abrupt and smooth transition which may be due to two distinct physical processes whose time-scales are eventually correlated. We discuss implications of our results in the light of some generic models.

preprint2013arXiv

The kilo-second variability of X-ray sources in nearby galaxies

Chandra observations of 17 nearby galaxies were analysed and 166 bright sources with X-ray counts > 100, were chosen for temporal analysis. Fractional root mean square variability amplitudes were estimated for lightcurves binned at ~ 4 ksec and of length ~ 40 ksec. While there are nine ultra-luminous X-ray sources (ULXs) with unabsorbed luminosity (in 0.3-8.0 keV band) L > 10^39 erg/s in the sample for which the fractional r.m.s variability is constrained to be < 10%, only two of them show variability. One of the variable ULXs exhibits a secular transition and has a ultra-soft spectrum with temperature ~ 0.3 keV while the other is a rapidly varying source in NGC 0628, which has been previously compared to the Galactic micro-quasar GRS1915+105. These results seem to indicate that ULXs are typically not highly variable in ksec time-scales, except for some ultra-soft ones. Among the relatively low luminosity sources (L ~ 10^38 erg/s) we find five of them to be variable. Apart from an earlier known source in NGC 1569, we identify a source in NGC 2403, which exhibits persistent high amplitude fluctuations.The variability of the sources in general, do not seem to be correlated with hardness,which indicates that they may not be due to variations in any absorbing material, but instead could reflect inner accretion disk instabilities.

preprint2013arXiv

The uncorrelated long term gamma-ray and X-ray variability of blazars and its implications on disk-jet coupling

We examine the long term (~10 years) gamma-ray variability of blazars observed by EGRET and Fermi and find that for six sources the average flux varied by more than an order of magnitude. For two of these sources (PKS 0208-512 and PKS 0528+134), there were extensive observations (at various observing periods) by EGRET. Hence these dramatic variations are not due to a single short time-scale flare, but reflect long term changes in the average flux. Over the last twenty years, these two sources were also the target of several X-ray observatories (e.g. ROSAT, ASCA, RXTE, BeppoSAX, Chandra, Suzaku, XMM-Newton and Swift). While the ratios of the average gamma-ray fluxes between EGRET and Fermi observations are 22.9 +/- 1.9 and 12.6 +/- 1.5, their estimated 2-10 keV X-ray flux do not show such dramatic variations. The X-ray emission from such flat spectrum radio quasars (FSRQs) are believed to be due to synchrotron self Compton, while gamma-rays originate from inverse Comptonization of external soft photons from an accretion disk and/or broad line region. We argue that in this scenario, the only explanation for the uncorrelated variability is that there was an order of magnitude decrease in the external soft photons, while the jet parameters remained more or less constant. This result indicates that perhaps the accretion and jet processes are not tightly coupled in these sources.

preprint2012arXiv

Energy-dependent orbital modulation of X-rays and constraints on emission of the jet in Cyg X-3

We study orbital modulation of X-rays from Cyg X-3, using data from Swift, INTEGRAL and RXTE. Using the wealth of the presently available data and an improved averaging method, we obtain energy-dependent folded and averaged light curves with unprecedented accuracy. We find that above ~5 keV, the modulation depth decreases with the increasing energy, which is consistent with the modulation being caused by both bound-free absorption and Compton scattering in the stellar wind of the donor, with minima corresponding to the highest optical depth, which occurs around the superior conjunction. We find a decrease of the depth below ~3 keV, which appears to be due to re-emission of the absorbed continuum by the wind in soft X-ray lines. Based on the shape of the folded light curves, any X-ray contribution from the jet in Cyg X-3, which emits gamma-rays detected at energies >0.1 GeV in soft spectral states, is found to be minor up to ~100 keV. This implies the presence of a rather sharp low-energy break in the jet MeV-range spectrum. We also calculate phase-resolved RXTE X-ray spectra, and show the difference between the spectra corresponding to phases around the superior and inferior conjunctions can indeed be accounted for by a combined effect of bound-free absorption in an ionized medium and Compton scattering.

preprint2012arXiv

Stability of black hole accretion disks

We discuss the issues of stability of accretion disks that may undergo the limit-cycle oscillations due to the two main types of thermal-viscous instabilities. These are induced either by the domination of radiation pressure in the innermost regions close to the central black hole, or by the partial ionization of hydrogen in the zone of appropriate temperatures. These physical processes may lead to the intermittent activity in AGN on timescales between hundreds and millions of years. We list a number of observational facts that support the idea of the cyclic activity in high accretion rate sources. We conclude however that the observed features of quasars may provide only indirect signatures of the underlying instabilities. Also, the support from the sources with stellar mass black holes, whose variability timescales are observationally feasible, is limited to a few cases of the microquasars. Therefore we consider a number of plausible mechanisms of stabilization of the limit cycle oscillations in high accretion rate accretion disks. The newly found is the stabilizing effect of the stochastic viscosity fluctuations.

preprint2012arXiv

Stabilization of radiation pressure dominated accretion disks by viscous fluctuations

The standard thin accretion disk model has been successfully used to explain the soft X-ray spectra of Galactic black hole systems and perhaps the UV emission of Active Galactic Nuclei. However, radiation pressure dominated disks are known to be viscously unstable and should produce large amplitude oscillations that are typically not observed. Instead, these sources exhibit stochastic variability which may naturally arise due to viscous fluctuations in a turbulent disk. Here we investigate whether these aperiodic viscous fluctuations can stabilize the inner radiation pressure dominated disks and hence maybe the answer to a forty year old problem in accretion disk theory.

preprint2011arXiv

Soft time lags in the X-ray emission of Mrk 1040

Temporal analysis of X-ray binaries and Active Galactic Nuclei have shown that hard X-rays react to variation of soft ones after a time delay. The opposite trend, or soft lag, has only been seen in a few rare Quasi-periodic Oscillations in X-ray binaries and recently for the AGN, 1H 0707-495, on short timescales of ~ 10^3 secs. Here, we report analysis of a XMM-Newton observation of Mrk~1040, which reveals that on the dominant variability timescale of ~ 10^4 secs, the source seems to exhibit soft lags. If the lags are frequency independent, they could be due to reverberation effects of a relativistically blurred reflection component responding to a varying continuum. Alternatively, they could be due to Comptonization delays in the case when high energy photons impinge back on the soft photon source. Both models can be verified and their parameters tightly constrained, because they will need to predict the photon spectrum, the r.m.s variability and time lag as a function of energy. A successful application of either model will provide unprecedented information on the radiative process, geometry and more importantly the size of the system, which in turn may provide stringent test of strong general relativistic effects.

preprint2007arXiv

Statistical Evidence for Three classes of Gamma-ray Bursts

Two different multivariate clustering techniques, the K-means partitioning method and the Dirichlet process of mixture modeling, have been applied to the BATSE Gamma-ray burst (GRB) catalog, to obtain the optimum number of coherent groups. In the standard paradigm, GRB are classified in only two groups, the long and short bursts. However, for both the clustering techniques, the optimal number of classes was found to be three, a result which is consistent with previous statistical analysis. In this classification, the long bursts are further divided into two groups which are primarily differentiated by their total fluence and duration and hence are named low and high fluence GRB. Analysis of GRB with known red-shifts and spectral parameters suggests that low fluence GRB have nearly constant isotropic energy output of 10^{52} ergs while for the high fluence ones, the energy output ranges from 10^{52} to 10^{54} ergs. It is speculated that the three kinds of GRBs reflect three different origins: mergers of neutron star systems, mergers between white dwarfs and neutron stars, and collapse of massive stars.