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 $\beta$ (i.e. $P(f) \propto f^{-\beta}$), 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 $\beta = 0.5$ and $1.0$, which differ from the values given in the literature which are for white noise ($\beta = 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 $\beta \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 $\beta\gtrsim 0.5$ and $\beta \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 $\gamma$ ray flux is consistent with a log-normal distribution and not with a Gaussian one.

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.

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 $\Gamma$ 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 $\Gamma$ 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-\Gamma$ relation is most likely related to the variations in the disc-corona geometry of AGN.

Results on timing and spectral properties of the Be/X-ray binary pulsar 3A 0726-260 (4U 0728-25) are presented. The binary was observed on 2016 May 6-7 with the Large Area X-ray Proportional Counter (LAXPC) and Soft X-ray Telescope (SXT) instruments onboard the AstroSat satellite. During this observation the source was in non-flaring persistent state at a flux level of $\sim$ 8.6 $\pm$ 0.3 $\times$10$^{-11}$ ergs cm$^{-2}$ sec$^{-1}$ in 0.4-20 keV. Strong X-ray pulsations with a period of 103.144$\pm$0.001 seconds are detected in 0.3-7 keV with the SXT and in 3-40 keV with the LAXPC. The pulse profile is energy dependent, and there is an indication that the pulse shape changes from a broad single pulse to a double pulse at higher energy. At energies above 20 keV, we report the first time detection of pulsation period 103.145$\pm$0.001 seconds and the double peaked pulse profile from the source. The energy spectrum of the source is derived from the combined analysis of the SXT and LAXPC spectral data in 0.4-20 keV. The best spectral fit is obtained by a power law model with a photon index (1.7$\pm$0.03) with high energy spectral cut-off at 12.9 $\pm$ 0.7 keV. A broad Iron line at $\sim$ 6.3 keV is detected in the energy spectrum. We briefly discuss the implications of these results.

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.

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.

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.

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

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.

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.

We report the analysis of simultaneous multi-wavelength data of the high energy peaked blazar RGB J0710+591 from the LAXPC, SXT and UVIT instruments on-board AstroSat. The wide band X-ray spectrum (0.35 -- 30 keV) is modelled as synchrotron emission from a non-thermal distribution of high energy electrons. The spectrum is unusually curved, with a curvature parameter $\beta_p \sim 6.4$ for a log parabola particle distribution, or a high energy spectral index $p_2 > 4.5$ for a broken power-law distribution. The spectrum shows more curvature than an earlier quasi-simultaneous analysis of Swift-XRT/NuSTAR data where the parameters were $\beta_p \sim 2.2$ or $p_2 \sim 4$. It has long been known that a power-law electron distribution can be produced from a region where particles are accelerated under Fermi process and the radiative losses in acceleration site decide the maximum attainable Lorentz factor, $\gamma_{max}$. Consequently, this quantity decides the energy at which the spectrum curves steeply. We show that such a distribution provides a more natural explanation for the AstroSat data as well as the earlier XRT/NuSTAR observation, making this as the first well constrained determination of the photon energy corresponding to $\gamma_{max}$. This in turn provides an estimate of the acceleration time-scale as a function of magnetic field and Doppler factor. The UVIT observations are consistent with earlier optical/UV measurements and reconfirm that they plausibly correspond to a different radiative component than the one responsible for the X-ray emission.

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{\sigma_{err}^{2}}/\sigma^{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.

We present the X-ray timing and spectral analysis of the new Galactic X-ray transient Swift J1658.2-4242 observed with LAXPC and SXT instruments onboard {\it Astrosat}. We detect prominent C-type quasi-periodic oscillations (QPOs) of frequencies varying from $\sim 1.5$ Hz to $\sim 6.6$ Hz along with distinct 2nd harmonics and sub-harmonics. The QPO detected at $\sim 1.56$ Hz drifts to a higher centroid frequency of $\sim 1.74$ in the course of the observation, while the QPO detected at $\sim 6.6$\,Hz disappeared during hard flarings. The fractional rms at the QPO and the sub-harmonic frequencies increases with photon energy, while at the 2nd harmonic frequencies the rms seems to be constant. In addition, we have observed soft time lag at QPO and sub-harmonic frequencies up to a time scale of $\sim 35$ ms, however, at the 2nd harmonic frequencies there is weak/zero time lag. We attempt spectral modeling of the broadband data in the 0.7--25 keV band using the doubly absorbed disk plus thermal Comptonization model. Based on the spectral and timing properties, we identified the source to be in the hard intermediate state of black hole X-ray binaries. To quantitatively fit the energy and frequency-dependent fractional rms and time lag, we use a single zone fluctuation propagation model and discuss our results in the context of that model.

We report on simultaneous sub-second optical and X-ray timing observations of the low mass X-ray binary black hole candidate MAXI J1820+070. The bright 2018 outburst rise allowed simultaneous photometry in five optical bands ($ugriz_s$) with HiPERCAM/GTC (Optical) at frame rates over 100 Hz, together with NICER/ISS observations (X-rays). Intense (factor of two) red flaring activity in the optical is seen over a broad range of timescales down to $\sim$10 ms. Cross-correlating the bands reveals a prominent anti-correlation on timescales of $\sim$seconds, and a narrow sub-second correlation at a lag of $\approx$+165 ms (optical lagging X-rays). This lag increases with optical wavelength, and is approximately constant over Fourier frequencies of $\sim$0.3-10 Hz. These features are consistent with an origin in the inner accretion flow and jet base within $\sim$5000 Gravitational radii. An additional $\sim$+5 s lag feature may be ascribable to disc reprocessing. MAXI J1820+070 is the third black hole transient to display a clear $\sim$0.1s optical lag, which may be common feature in such objects. The sub-second lag $variation$ with wavelength is novel, and may allow constraints on internal shock jet stratification models.

Close binary stars are binary stars where the component stars are close enough such that they can exchange mass and/or energy. They are subdivided into semi-detached, overcontact or ellipsoidal binary stars. A challenging problem in the context of close binary stars, is their classification into these subclasses, based solely on their light curves. Conventionally, this is done by observing subtle features in the light curves like the depths of adjacent minima, which is tedious when dealing with large datasets. In this work we suggest the use of machine learning algorithms applied to measures of recurrence networks and nonlinear time series analysis to differentiate between classes of close binary stars. We show that overcontact binary stars occupy a region different from semi-detached and ellipsoidal binary stars in a plane of characteristic path length(CPL) and average clustering coefficient(CC), computed from their recurrence networks. We use standard clustering algorithms and report that the clusters formed corresponds to the standard classes with a high degree of accuracy.

We report the results of the analysis of an AstroSat observation of the Black Hole candidate MAXI J1535-571 during its Hard Intermediate state. We studied the evolution of the spectral and timing parameters of the source during the observation. The observation covered a period of $\sim$5 days and consisted of 66 continuous segments, corresponding to individual spacecraft orbits. Each segment was analysed independently. The source count rate increased roughly linearly by $\sim$30 %. We modelled the spectra as a combination of radiation from a thermal disk component and a power-law. The timing analysis revealed the presence of strong Quasi Periodic Oscillations with centroid frequency $\nu_{\rm{QPO}}$ fluctuating in the range 1.7-3.0 Hz. We found a tight correlation between the QPO centroid frequency $\nu_{\rm{QPO}}$ and the power-law spectral index $\Gamma$, while $\nu_{\rm{QPO}}$ appeared not to be correlated with the linearly-increasing flux itself. We discuss the implications of these results on physical models of accretion.

We present rapid, multiwavelength photometry of the low-mass X-ray binary Swift J1357.2-0933 during its 2017 outburst. Using several sets of quasi-simultaneous ULTRACAM/NTT (optical), NuSTAR (X-ray), XRT/Swift (X-ray), SALT (optical) and ATCA (radio) observations taken during outburst decline, we confirm the frequent optical dipping that has previously been noted both in outburst and in quiescence. We also find: 1) that the dip frequency decreases as the outburst decays, similar to what was seen in the previous outburst, 2) that the dips produce a shape similar to that in binary systems with partial disc occultations, 3) that the source becomes significantly bluer during these dips, indicating an unusual geometry compared to other LMXB dippers, and 4) that dip superposition analysis confirms the lack of an X-ray response to the optical dips. These very unusual properties appear to be unique to Swift J1357.2-0933, and are likely the result of a high binary inclination, as inferred from features such as its very low outburst X-ray luminosity. From this analysis as well as X-ray/optical timing correlations, we suggest a model with multi-component emission/absorption features with differing colours. This could include the possible presence of a sporadically occulted jet base and a recessed disc. This source still hosts many puzzling features, with consequences for the very faint X-ray transients population.

We report the results from analysis of six observations of Cygnus X-1 by Large Area X-ray Proportional Counters (LAXPC) and Soft X-ray Telescope (SXT) on-board AstroSat, when the source was in the hard spectral state as revealed by the broad band spectra. The spectra obtained from all the observations can be described by a single temperature Comptonizing region with disk and reflection components. The event mode data from LAXPC provides unprecedented energy dependent fractional root mean square (rms) and time-lag at different frequencies which we fit with empirical functions. We invoke a fluctuation propagation model for a simple geometry of a truncated disk with a hot inner region. Unlike other propagation models, the hard X-ray emission (> 4 keV) is assumed to be from the hot inner disk by a single temperature thermal Comptonization process. The fluctuations first cause a variation in the temperature of the truncated disk and then the temperature of the inner disk after a frequency dependent time delay. We find that the model can explain the energy dependent rms and time-lag at different frequencies.

Overcontact binary stars are systems of two stars where the component stars are in contact with each other. This implies that they share a common envelope of gas. In this work we seek signatures of nonlinearity and chaos in these stars by using time series analysis techniques. We use three main techniques, namely the correlation dimension,f (\alpha) spectrum and the bicoherence. The former two are calculated from the reconstructed dynamics, while the latter is calculated from the Fourier transforms of the time series of intensity variations(light curves) of these stars. Our dataset consists of data from 463 overcontact binary stars in the Kepler field of view [1]. Our analysis indicates nonlinearity and signatures of chaos in almost all the light curves. We also explore whether the underlying nonlinear properties of the stars are related to their physical properties like fill-out-factor, a measure of the extend of contact between the components of an overcontact binary system . We observe that significant correlations exist between the fill out factor and the nonlinear quantifiers. This correlation is more pronounced in specific subcategories constructed based on the mass ratios and effective temperatures of the binaries. The correlations observed can be indicative of variations in the nonlinear properties of the star as it ages. We believe that this study relating nonlinear and astrophysical properties of binary stars is the first of its kind and is an important starting point for such studies in other astrophysical objects displaying nonlinear dynamical behaviour.

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.

The Be X-ray Binary 4U 0115+63 was observed by Large Area X-ray Proportional Counter (LAXPC) instrument on AstroSat on 2015 October 24 during the peak of a giant Type II outburst. Prominent intensity oscillations at ~ 1 and ~ 2 mHz frequency were detected during the outburst. Nuclear Spectroscopic Telescope Array (NuSTAR) observations made during the same outburst also show mHz quasi periodic oscillations (QPOs). Details of the oscillations and their characteristics deduced from LAXPC/AstroSat and NuSTAR observations are reported in this paper. Analysis of the archival Rossi X-ray Timing Explorer (RXTE)/Proportional Counter Array (PCA) data during 2001-11 also show presence of mHz QPOs during some of the outbursts and details of these QPOs are also reported. Possible models to explain the origin of the mHz oscillations are examined. Similar QPOs, albeit at higher frequencies, have been reported from other neutron star and black hole sources and both may have a common origin. Current models to explain the instability in the inner accretion disk causing the intense oscillations are discussed.

We present a detailed temporal and spectral study of the blazar 3C\,279 using multi-wavelength observations from Swift-XRT, Swift-UVOT and Fermi-LAT during a flare in 2018 January. The temporal analysis of $\gamma$-ray light curve indicates a lag of $\sim 1$ d between the 0.1--3 GeV and 3--500 GeV emission. Additionally, the $\gamma$-ray light curve shows asymmetry with slow rise--fast decay in energy band 0.1--3 GeV and fast rise--slow decay in the 3--500 GeV band. We interpret this asymmetry as a result of shift in the Compton spectral peak. This inference is further supported by the correlation studies between the flux and the parameters of the log-parabola fit to the source spectra in the energy range 0.1--500 GeV. We found that the flux correlates well with the peak spectral energy and the log-parabola fit parameters show a hard index with large curvature at high flux states. Interestingly, the hardest index with large curvature was synchronous with a very high energy flare detected by H.E.S.S. Our study of the spectral behavior of the source suggests that $\gamma$-ray emission is most likely to be associated with the Compton up-scattering of IR photons from the dusty environment. Moreover, the fit parameters indicate the increase in bulk Lorentz factor of emission region to be a dominant cause for the flux enhancement.

We present results from an observation of the Low Mass X-ray Binary 4U 1636-536 obtained with the LAXPC instrument aboard AstroSat. The observations of 4U 1636-536 made during the performance verification phase of AstroSat showed seven thermonuclear X-ray bursts in a total exposure of ~ 65 ks over a period of about two consecutive days. Moreover, the light curve of 4U 1636-536 revealed the presence of a rare triplet of X-ray bursts, having a wait time of about 5.5 minutes between second and the third bursts. We also present results from time-resolved spectroscopy performed during these seven X-ray bursts. In addition, we have also detected a transient Quasi-periodic oscillation (QPO) at ~ 5 Hz. However, we did not find any evidence of kilo-hertz QPOs and/or X-ray burst oscillations, perhaps due to the hard spectral state of the source during this observation.

In this work, we present a $\sim$90 ks continuous monitoring of the Galactic micro-quasar GRS 1915+105 with AstroSat when the source undergoes a major transition from a non-variable, $\chi$ class (similar to radio-quiet $\chi$ class) to a structured, large amplitude, periodic heartbeat state (similar to $\rho$ class). We show that such transition takes place via an intermediate state when the large-amplitude, irregular variability of the order of hundreds of seconds in the soft X-ray band turned into a 100-150 sec regular, structured, nearly periodic flares. The properties of a strong low-frequency quasi-periodic oscillations (LF QPO) in the frequency range 3-5 Hz also evolve marginally during these variability transitions. We also study time-lag and rms spectra at the QPO and harmonic component and the dynamic power spectra. We note few important differences between the heartbeat state and the $\rho$ class. Interestingly, the time-averaged LF QPO properties in the hard X-ray band is relatively stable in three states when compared to the significant evolution observed in the slow variability properties at mHz frequencies. Such relative stability of LF QPOs implies the inner disk-corona coupled accretion flow which determines the LF QPO properties, may be uninterrupted by the launch of long, large-amplitude flares.

We present X-ray spectral variability studies of the ultra-luminous X-ray source (ULX) M81 X--6 using {\it Suzaku} and {\it XMM-Newton} observations performed during 2001--2015. The spectra were first fitted by a standard multi-temperature disk and a thermal Comptonization component which revealed spectral variability where the primary distinction is the change in the optical depth of the Comptonizing component, similar to what has been observed for other ULXs. We also fitted the spectra with a general relativistic accretion disk emission and a power-law component and found that it can reproduce a large part but not all of the spectral variability of the source. The parameters for the black hole mass and spin were found to be degenerate, but the high spin and larger mass ($20-100\,\rm M_\odot$) solutions provided near Eddington accretion rates consistent with the assumptions of the model. The spectral variation is found to be driven by accretion rate changes leading to three different spectral classes. Thus, our results suggest the possibility of a dominant relativistic disk emission component for some of the spectral states of the source.

The X-ray spectral curvature of blazars is traditionally explained by an empirical log-parabola function characterized by three parameters, namely the flux, curvature and spectral index at a given energy. Since their exact relationship with the underlying physical quantities is unclear, interpreting the physical scenario of the source through these parameters is difficult. To attain an insight on the X-ray spectral shape, we perform a detailed study of the X-ray spectra of the blazar MKN 421, using an analytical model where the electron diffusion from the particle acceleration site is energy-dependent. The resultant synchrotron spectrum is again determined by three parameters, namely, the energy index of the escape time scale, the quantity connecting the electron energy to the observed photon energy and the normalization. The X-ray observations of MKN 421, during July 2012 - April 2013 by NuSTAR and Swift-XRT are investigated using this model and we find a significant correlation between model parameters and the observational quantities. Additionally, a strong anti-correlation is found between the fit parameters defining the spectral shape, which was not evident from earlier studies using empirical models. This indicates the flux variations in MKN 421 and possibly other blazars, may arise from a definite physical process that needs to be further investigated.

The observed log-normal flux distributions in the high energy emission from blazars have been interpreted as being due to variability stemming from non-linear multiplicative processes generated dynamically from the accretion disc. On the other hand, rapid minute scale variations in the flux point to a compact emitting region inside the jet, probably disconnected from the disc. In this work, we show that linear Gaussian variations of the intrinsic particle acceleration or escape time-scales can produce distinct non-Gaussian flux distributions, including log-normal ones. Moreover, the spectral index distributions can provide confirming evidence for the origin of the variability. Thus, modelling of the flux and index distributions can lead to quantitative identification of the micro-physical origin of the variability in these sources. As an example, we model the X-ray flux and index distribution of Mkn 421 obtained from over 9 years of MAXI observations and show that the variability in the X-ray emission is driven by Gaussian fluctuations of the particle acceleration process rather than that of the escape rate.

We present the first results from a detailed spectral-timing analysis of a long ($\sim$130 ks) XMM-Newton observation and quasi-simultaneous NuSTAR and Swift observations of the highly-accreting narrow-line Seyfert 1 galaxy Mrk 1044. The broadband (0.3$-$50 keV) spectrum reveals the presence of a strong soft X-ray excess emission below $\sim$1.5 keV, iron K$_{\alpha}$ emission complex at $\sim$6$-$7 keV and a `Compton hump' at $\sim$15$-$30 keV. We find that the relativistic reflection from a high-density accretion disc with a broken power-law emissivity profile can simultaneously explain the soft X-ray excess, highly ionized broad iron line and the Compton hump. At low frequencies ($[2-6]\times10^{-5}$ Hz), the power-law continuum dominated 1.5$-$5 keV band lags behind the reflection dominated 0.3$-$1 keV band, which is explained with a combination of propagation fluctuation and Comptonization processes, while at higher frequencies ($[1-2]\times10^{-4}$ Hz), we detect a soft lag which is interpreted as a signature of X-ray reverberation from the accretion disc. The fractional root-mean-squared (rms) variability of the source decreases with energy and is well described by two variable components: a less variable relativistic disc reflection and a more variable direct coronal emission. Our combined spectral-timing analyses suggest that the observed broadband X-ray variability of Mrk~1044 is mainly driven by variations in the location or geometry of the optically thin, hot corona.

We present the results of the $\gamma$-ray flux distribution study on the brightest blazars which are observed by the \emph{Fermi}-LAT. We selected 50 brightest blazars based on the maximum number of detection reported in the LAT third AGN catalog. We performed standard unbinned maximum likelihood analysis on the LAT data during the period between August 2008 and December 2016, in order to obtain the average monthly flux. After quality cuts, blazars for which at least 90\% of the total flux was survived were selected for the further study, and this includes 19 FSRQs and 19 BL Lacs. The Anderson-Darling and $\chi^2$ tests suggest that the integrated monthly flux follow a log-normal distribution for all sources, except for three FSRQs for which neither a normal nor a log-normal distribution was preferred. A double log-normal flux distribution tendency were observed in these sources, though it has to be confirmed with improved statistics. We also found that, the standard deviation of the log-normal flux distribution increases with the mean spectral index of the blazar, and can be fitted with a line of slope 0.24$\pm$0.04. We repeat our study on three additional brightest unclassified blazars to identify their flux distribution properties. Based on the features of their log-normal flux distribution, we infer these unclassified blazars may be closely associated with FSRQs. We also highlight that considering the log-normal behavior of the flux distribution of blazars, averaging their long term flux in linear scale can largely under estimate the nominal flux and this discrepancy can propagate down to the estimation of source parameters through spectral modeling.

It has recently been shown that the persistent emission of a neutron star low-mass X-ray binary (LMXB) evolves during a thermonuclear (type-I) X-ray burst. The reason of this evolution, however, is not securely known. This uncertainty can introduce significant systematics in the neutron star radius measurement using burst spectra, particularly if an unknown but significant fraction of the burst emission, which is reprocessed, contributes to the changes in the persistent emission during the burst. Here, by analyzing individual burst data of AstroSat/LAXPC from the neutron star LMXB 4U 1728--34 in the soft state, we show that the burst emission is not significantly reprocessed by a corona covering the neutron star. Rather, our analysis suggests that the burst emission enhances the accretion disk emission, possibly by increasing the accretion rate via disk. This enhanced disk emission, which is Comptonized by a corona covering the disk, can explain an increased persistent emission observed during the burst. This finding provides an understanding of persistent emission components, and their interaction with the thermonuclear burst emission. Furthermore, since burst photons are not significantly reprocessed, non-burst and burst emissions can be reliably separated, which is required to reduce systematic uncertainties in the stellar radius measurement.

The broad band optical/UV and X-ray spectra of blazars have been often modeled as synchrotron component arising from a broken power-law distribution of electrons. A broken power-law distribution is expected, since the high energy electrons undergo radiative losses effectively. The change in the energy index should then be approximately 1 and corresponds to a spectral index difference of 0.5. However, one of the long outstanding problems has been that the observed index change is significantly different. On the other hand, recent high quality observations of blazars suggest that their local spectra may not be a power-law, instead have a slight curvature and often represented by a log parabola model. Using XMM-Newton observations spanning over 12 years for the BL Lac PKS 2155-304, we show that the optical/UV and X-ray spectra can be well represented by a broken log parabola model. Further, we show that such a spectrum can indicate the energy dependence of the electron escape timescale from the main acceleration zone. This novel approach, besides addressing the observed difference in the photon spectral indices, also tries to explain the spectral turn over in far UV/soft X-rays as a result of the radiative losses.

One of the popular models for the low/hard state of Black Hole Binaries is that the standard accretion disk is truncated and the hot inner region produces via Comptonization, the hard X-ray flux. This is supported by the value of the high energy photon index, which is often found to be small $\sim$ 1.7 ($<$ 2) implying that the hot medium is seed photons starved. On the other hand, the suggestive presence of a broad relativistic Fe line during the hard state would suggest that the accretion disk is not truncated but extends all the way to the inner most stable circle orbit. In such a case, it is a puzzle why the hot medium would remain photon starved. The broad Fe line should be accompanied by a broad smeared reflection hump at $\sim$ 30 keV and it may be that this additional component makes the spectrum hard and the intrinsic photon index is larger, i.e. $>$ 2. This would mean that the medium is not photon deficient, reconciling the presence of a broad Fe line in the observed hard state. To test this hypothesis, we have analyzed the RXTE observations of GX 339-4 from the four outbursts during 2002-2011 and identify the observations when the system was in the hard state and showed a broad Fe line. We have then attempted to fit these observations with models, which include smeared reflection to understand whether the intrinsic photon index can indeed be large. We find that, while for some observations the inclusion of reflection does increase the photon index, there are hard state observations with broad Fe line that have photon indices less than 2.

We present X-ray spectral and timing behavior of Cyg X-3 as observed by AstroSat during the onset of a giant radio flare on 01-02 April 2017. Within a time-scale of few hours, the source shows a transition from the hypersoft state (HPS) to a more luminous state (we termed as the very high state) which coincides with the time of the steep rise in radio flux density by an order of magnitude. Modeling the SXT and LAXPC spectra jointly in 0.5-70.0 keV, we found that the first few hours of the observation is dominated by the HPS with no significant counts above 17 keV. Later, an additional flat powerlaw component suddenly appeared in the spectra which extends to very high energies with the powerlaw photon index of 1.49 +/- 0.04. Such a flat powerlaw component has never been reported from Cyg X-3. Interestingly the fitted powerlaw model in 25-70 keV, when extrapolated to the radio frequency, predicts the radio flux density consistent with the trend measured from RATAN-600 telescope at 11.2 GHz. This provides a direct evidence of the synchrotron origin of flat X-ray powerlaw component and the most extensive monitoring of the broadband X-ray behavior at the moment of decoupling the giant radio jet base from the compact object in Cyg X-3. Using SXT and LAXPC observations, we determine the giant flare ejection time as MJD 57845.34 +/- 0.08 when 11.2 GHz radio flux density increases from ~100 to ~478 mJy.

The broadband spectral energy distribution(SED) of blazars is generally interpreted as radiation arising from synchrotron and inverse Compton mechanisms.Traditionally, the underlying source parameters responsible for these emission processes,like particle energy density, magnetic field, etc., are obtained through simple visual reproduction of the observed fluxes. However, this procedure is incapable of providing the confidence range on the estimated parameters. In this work, we propose an efficient algorithm to perform a statistical fit of the observed broadband spectrum of blazars using different emission models. Moreover, in this work we use the the observable quantities as the fit parameters, rather than the direct source parameters which govern the resultant SED. This significantly improves the convergence time and eliminates the uncertainty regarding the initial guess parameters. This approach also has an added advantage of identifying the degenerate parameters, which can be removed by including more observable information and/or additional constraints. A computer code developed based on this algorithm is implemented as an user-defined routine in the standard X-ray spectral fitting package, XSPEC. Further, we demonstrate the efficacy of the algorithm by fitting the well sampled SED of the blazar, 3C 279, during its gamma ray flare in 2014.

The synchrotron peak of the X-ray bright High Energy Peaked Blazar (HBL) PKS 2155$-$304 occurs in the UV-EUV region and hence its X-ray emission (0.6--10 keV) lies mostly in the falling part of the synchrotron hump. We aim to study the X-ray emission of PKS 2155$-$304 during different intensity states in 2009$-$2014 using XMM$-$Newton satellite. We studied the spectral curvature of all of the observations to provide crucial information on the energy distribution of the non-thermal particles. Most of the observations show curvature or deviation from a single power-law and can be well modeled by a log parabola model. In some of the observations, we find spectral flattening after 6 keV. In order to find the possible origin of the X-ray excess, we built the Multi-band Spectral Energy distribution (SED). We find that the X-ray excess in PKS 2155--304 is difficult to fit in the one zone model but, could be easily reconciled in the spine/layer jet structure. The hard X-ray excess can be explained by the inverse Comptonization of the synchrotron photons (from the layer) by the spine electrons.

We present detailed broadband UV/optical to X-ray spectral variability of the Seyfert 1 galaxy 1H 0419-577 using six XMM-Newton observations performed during 2002-2003. These observations covered a large amplitude variability event in which the soft X-ray (0.3-2 keV) count rate increased by a factor of ~4 in six months. The X-ray spectra during the variability are well described by a model consisting of a primary power law, blurred and distant reflection. The 2-10 keV power-law flux varied by a factor ~7 while the 0.3-2 keV soft X-ray excess flux derived from the blurred reflection component varied only by a factor of ~2. The variability event was also observed in the optical and UV bands but the variability amplitudes were only at the 6-10% level. The variations in the optical and UV bands appear to follow the variations in the X-ray band. During the rising phase, the optical bands appear to lag behind the UV band but during the declining phase, the optical bands appear to lead the UV band. Such behavior is not expected in the reprocessing models where the optical/UV emission is the result of reprocessing of X-ray emission in the accretion disc. The delayed contribution of the broad emission lines in the UV band or the changes in the accretion disc/corona geometry combined with X-ray reprocessing may give rise to the observed behavior of the variations.

We present a statistical characterization of the $\gamma$-ray emission from the four \emph{Fermi}-LAT sources: FR I radio galaxy NGC 1275, BL Lac Mrk 421, FSRQs B2 1520+31 and PKS 1510-089 detected almost continuously over a time integration of 3-days between August 2008 - October 2015. The observed flux variation is large, spanning $\gtrsim 2$ orders of magnitude between the extremes except for Mrk~421. We compute the flux distributions and compare with Gaussian and lognormal ones. We find that the 3 blazars have distribution consistent with a lognormal, suggesting that the variability is of a non-linear, multiplicative nature. This is further supported by the computation of the flux-rms relation, which is observed to be linear for the 3 blazars. However, for NGC 1275, the distribution does not seem to be represented either by a lognormal or a Gaussian, while its flux-rms relation is still found to be linear. We also compute the power spectra, which suggest the presence of a break, but are consistent with typical scale-free power-law shot noise. The results are broadly consistent with the statistical properties of the magnetic reconnection powered minijets-in-a-jet model. We discuss other possible scenarios and implications of these observations on jet processes and connections with the central engine.

We investigate the long-term spectral variability in the ultra-luminous X-ray source Holmberg IX X--1. By analyzing the data from eight {\it Suzaku} and 13 {\it XMM-Newton} observations conducted between 2001 and 2015, we perform a detailed spectral modeling for all spectra with simple models and complex physical models. We find that the spectra can be well explained by a disc plus thermal Comptonization model. Applying this model, we unveil correlations between the X-ray luminosity ($L_{\rm X}$) and the spectral parameters. Among the correlations, a particular one is the statistically significant positive correlation between $L_{\rm X}$ and the photon index ($\Gamma$), while at the high luminosities of $> 2\times10^{40}\,{\rm~erg\ s}^{-1}$, the source becomes marginally hard and that results a change in the slope of the $\Gamma - L_{\rm X}$ correlation. Similar variability behavior is observed in the optical depth of the source around $L_{\rm X} \sim 2\times10^{40}\,{\rm~erg\ s}^{-1}$ as the source becomes more optically thick. We consider the scenario that a corona covers the inner part of the disc, and the correlations can be explained as to be driven by the variability of seed photons from the disc input into the corona. On the basis of the disc-corona model, we discuss the physical processes that are possibly indicated by the variability of the spectral parameters. Our analysis reveals the complex variability behavior of Holmberg IX X--1 and the variability mechanism is likely related to the geometry of the X-ray emitting regions.

We present here results from the X-ray timing and spectral analysis of the X-ray binary Cyg X-3 using observations from Large Area X-ray Proportional Counter (LAXPC) on-board AstroSat. Consecutive lightcurves observed over a period of one year show the binary orbital period of 17253.56 +/- 0.19 sec. Another low-amplitude, slow periodicity of the order of 35.8 +/- 1.4 days is observed which may be due to the orbital precession as suggested earlier by Molteni et al. (1980). During the rising binary phase, power density spectra from different observations during flaring hard X-ray state show quasi-periodic oscillations (QPOs) at ~5-8 mHz, ~12-14 mHz, ~18-24 mHz frequencies at the minimum confidence of 99%. However, during the consecutive binary decay phase, no QPO is detected up to 2-sigma significance. Energy-dependent time-lag spectra show soft lag (soft photons lag hard photons) at the mHz QPO frequency and the fractional rms of the QPO increases with the photon energy. During the binary motion, the observation of mHz QPOs during the rising phase of the flaring hard state may be linked to the increase in the supply of the accreting material in the disk and corona via stellar wind from the companion star. During the decay phase, the compact source moves in the outer wind region causing the decrease in the supply of material for accretion. This may cause weakening of the mHz QPOs below the detection limit. This is also consistent with the preliminary analysis of the orbital phase-resolved energy spectra presented in this paper.

In the unified scheme of active galactic nuclei, a dusty torus absorbs and then reprocesses a fraction of the intrinsic luminosity which is emitted at longer wavelengths. Thus, subject to radiative transfer corrections, the fraction of the sky covered by the torus as seen from the central source (known as the covering factor $f_c$) can be estimated from the ratio of the infrared to the bolometric luminosities of the source as $f_c=L_{\rm torus}/L_{\rm Bol}$. However, the uncertainty in determining $L_{\rm Bol}$ has made the estimation of covering factors by this technique difficult, especially for AGN in the local Universe where the peak of the observed SEDs lies in the UV (ultraviolet). Here, we determine the covering factors of an X-ray/optically selected sample of 51 type~1 AGN. The bolometric luminosities of these sources are derived using a self-consistent, energy-conserving model that estimates the contribution in the unobservable far-UV region, using multi-frequency data obtained from SDSS, \textit{XMM-Newton}, \textit{WISE}, 2MASS and UKIDSS. We derive a mean value of $f_c\sim$0.30 with a dispersion of 0.17. Sample correlations, combined with simulations, show that $f_c$ is more strongly anti-correlated with $\lambda_{\rm Edd}$ than with $L_{\rm Bol}$. This points to large-scale torus geometry changes associated with the Eddington-dependent accretion flow, rather than a receding torus, with its inner sublimation radius determined solely by heating from the central source. Furthermore, we do not see any significant change in the distribution of $f_c$ for sub-samples of radio-loud sources or Narrow Line Seyfert~1 galaxies (NLS1s), though these sub-samples are small.

The observed variability of X-ray binaries over a wide range of time-scales can be understood in the framework of a stochastic propagation model, where viscous fluctuations at different radii induce accretion rate variability that propagate inwards to the X-ray producing region. The scenario successfully explains the power spectra, the linear rms-flux relation as well as the time-lag between different energy photons. The predictions of this model have been obtained using approximate analytical solutions or empirically motivated models which take into account the effect of these propagating variability on the radiative process of complex accretion flows. Here, we study the variation of the accretion rate due to such viscous fluctuations using a hydro-dynamical code for the standard geometrically thin, gas pressure dominated $\alpha$-disc with a zero torque boundary condition. Our results confirm earlier findings that the time-lag between a perturbation and the resultant inner accretion rate variation depends on the frequency (or time-period) of the perturbation. Here we have quantified that the time-lag $t_{lag} \propto f^{-0.54}$, for time-periods less than the viscous time-scale of the perturbation radius and is nearly constant otherwise. This, coupled with radiative process would produce the observed frequency dependent time-lag between different energy bands. We also confirm that if there are random Gaussian fluctuations of the $\alpha$-parameter at different radii, the resultant inner accretion rate has a power spectrum which is a power-law.

Recurrence networks and the associated statistical measures have become important tools in the analysis of time series data. In this work, we test how effective the recurrence network measures are in analyzing real world data involving two main types of noise, white noise and colored noise. We use two prominent network measures as discriminating statistic for hypothesis testing using surrogate data for a specific null hypothesis that the data is derived from a linear stochastic process. We show that the characteristic path length is especially efficient as a discriminating measure with the conclusions reasonably accurate even with limited number of data points in the time series. We also highlight an additional advantage of the network approach in identifying the dimensionality of the system underlying the time series through a convergence measure derived from the probability distribution of the local clustering coefficients. As examples of real world data, we use the light curves from a prominent black hole system and show that a combined analysis using three primary network measures can provide vital information regarding the nature of temporal variability of light curves from different spectroscopic classes.

Deriving meaningful information from observational data is often restricted by many limiting factors, the most important of which is the presence of noise. In this work, we present the use of the bicoherence function to extract information about the underlying nonlinearity from noisy time series. We show that a system evolving in the presence of noise which has its dynamical state concealed from quantifiers like the power spectrum and correlation dimension D2, can be revealed using the bicoherence function. We define an index called main peak bicoherence function as the bicoherence associated with the maximal power spectral peak. We show that this index is extremely useful while dealing with quasi-periodic data as it can distinguish strange non chaos from quasi periodicity even with added noise. We demonstrate this in a real world scenario, by taking the bicoherence of variable stars showing period doubling and strange non-chaotic behavior. Our results indicate that bicoherence analysis can also bypass the method of surrogate analysis using Fourier phase randomization, used to differentiate linear stochastic processes from non linear ones, in conventional methods involving measures like D2.

We present the calibration and background model for the Large Area X-ray Proportional Counter (LAXPC) detectors on-board AstroSat. LAXPC instrument has three nominally identical detectors to achieve large collecting area. These detectors are independent of each other and in the event analysis mode, they record the arrival time and energy of each photon that is detected. The detectors have a time-resolution of 10 $\mu$s and a dead-time of about 42 $\mu$s. This makes LAXPC ideal for timing studies. The energy resolution and peak channel to energy mapping were obtained from calibration on ground using radioactive sources coupled with GEANT4 simulations of the detectors. The response matrix was further refined from observations of the Crab X-ray source after launch. At around 20 keV the energy resolution of detector is about 10--15\%, while the combined effective area of the 3 detectors is about 6000 cm$^2$.

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.

Large Area X-ray Propositional Counter (LAXPC) instrument on AstroSat is aimed at providing high time resolution X-ray observations in 3 to 80 keV energy band with moderate energy resolution. To achieve large collecting area, a cluster of three co-aligned identical LAXPC detectors, is used to realize an effective area in access of about 6000 cm2 at 15 keV. The large detection volume of the LAXPC detectors, filled with xenon gas at about 2 atmosphere pressure, results in detection efficiency greater than 50%, above 30 keV. In this article, we present salient features of the LAXPC detectors, their testing and characterization in the laboratory prior to launch and calibration in the orbit. Some preliminary results on timing and spectral characteristics of a few X-ray binaries and other type of sources, are briefly discussed to demonstrate that the LAXPC instrument is performing as planned in the orbit.

Large Area X-ray Proportional Counter (LAXPC) is one of the major AstroSat payloads. LAXPC instrument will provide high time resolution X-ray observations in 3 to 80 keV energy band with moderate energy resolution. A cluster of three co-aligned identical LAXPC detectors is used in AstroSat to provide large collection area of more than 6000 cm2 . The large detection volume (15 cm depth) filled with xenon gas at about 2 atmosphere pressure, results in detection efficiency greater than 50%, above 30 keV. With its broad energy range and fine time resolution (10 microsecond), LAXPC instrument is well suited for timing and spectral studies of a wide variety of known and transient X-ray sources in the sky. We have done extensive calibration of all LAXPC detectors using radioactive sources as well as GEANT4 simulation of LAXPC detectors. We describe in brief some of the results obtained during the payload verification phase along with LXAPC capabilities.

We perform a detailed spectral study of a recent flaring activity from the Flat Spectrum Radio Quasar (FSRQ), 3C\,454.3, observed simultaneously in optical, UV, X-ray and $\gamma$-ray energies during 16 to 28 August, 2015. The source reached its peak $\gamma$-ray flux of $(1.9\pm0.2)\times\,10^{-05} \; {\rm ph\,cm^{-2}\,s^{-1}}$ on 22 August. The time averaged broadband spectral energy distribution (SED) is obtained for three time periods, namely "flaring state"; covering the peak $\gamma$-ray flux, "post flaring state"; immediately following the peak flare and "quiescent state"; separated from the flaring event and following the post flaring state. The SED corresponding to the flaring state is investigated using different emission models involving synchrotron, synchrotron self Compton (SSC) and external Compton (EC) mechanisms. Our study suggests that the X-ray and $\gamma$-ray emission from 3C\,454.3 cannot be attributed to a single emission mechanism and instead, one needs to consider both SSC and EC mechanisms. Moreover, the target photon energy responsible for the EC process corresponds to an equivalent temperature of 564 K, suggesting that the flare location lies beyond the broad line emitting region of the FSRQ. SED fitting of the other two flux states further supports these inferences.

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 secondsdynamic power spectra reveal a burst oscillation whose frequency increased from 361.5 to 363.5Hz. 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.

Temporal analysis of radiation from Astrophysical sources like Active Galactic Nuclei, X-ray Binaries and Gamma-ray bursts provide information on the geometry and sizes of the emitting regions. Establishing that two light-curves in different energy bands are correlated and measuring the phase and time-lag between them is an important and frequently used temporal diagnostic. In the present work, expressions to estimate the errors on the cross-correlation, phase and time-lag between two light-curves are presented and the same have been verified using simulations. Earlier estimates depended upon numerically expensive simulations or on dividing the light-curves in large number of segments to find the variance. The estimates presented here allow for analysis of light-curves with relatively small (~1000) number of points, as well as to obtain information on the longest time-scales available. The error estimation is verified using simulations of light-curves derived from both white and 1/f stochastic processes with measurement errors. As a demonstration, we apply this technique to the XMM-Newton light-curves of the Active Galactic Nucleus, Akn 564.

We present results from simultaneous optical (SALT) and X-ray (Swift and INTEGRAL) observations of GS 1354-64/BW Cir during the 2015 hard state outburst. During the rising phase, optical/X-ray time series show a strong anti-correlation with X-ray photons lagging optical. Optical and X-ray power spectra show quasi-periodic oscillations at a frequency of ~18 mHz with a confidence level of at least 99%. Simultaneous fitting of Swift/XRT and INTEGRAL spectra in the range 0.5-1000 keV shows non-thermal, power-law dominated (> 90%) spectra with a hard power-law index of 1.48 +/- 0.03, inner disc temperature of 0.12 +/- 0.01 keV and inner disc radius of ~3000 km. All evidence is consistent with cyclo-synchrotron radiation in a non-thermal, hot electron cloud extending to ~100 Schwarzschild radii being a major physical process for the origin of optical photons. At outburst peak about one month later, when the X-ray flux rises and the optical drops, the apparent features in the optical/X-ray correlation vanish and the optical auto correlation widens. Although ~0.19 Hz QPO is observed from the X-ray power spectra, the optical variability is dominated by the broadband noise, and the inner disc temperature increases. These results support a change in the dominant optical emission source between outburst rise and peak, consistent with a weakening of hot flow as the disc moves in.

The observations of the prompt emission of gamma ray bursts (GRB) by GLAST Burst Monitor (GBM), on board Fermi Gamma-ray Space Telescope, suggest the presence of a significant thermal spectral component, whose origin is not well understood. Recently, it has been shown that for long duration GRBs, the spectral width as defined as the logarithm of the ratio of the energies at which the spectrum falls to half its peak value, lie in the range of 0.84-1.3 with a median value of 1.07. Thus, while most of the GRB spectra are found to be too narrow to be explained by synchrotron emission from an electron distribution, they are also significantly broader than a blackbody spectrum whose width should be 0.54. Here, we consider the possibility that an intrinsic thermal spectrum from a fire-ball like model, may be observed to be broadened if the system undergoes a rapid temperature evolution. We construct a toy-model to show that for bursts with durations in the range 5-70 seconds, the widths of their 1 second time-averaged spectra can be at the most $\lesssim 0.557$. Thus, while rapid temperature variation can broaden the detected spectral shape, the observed median value of $\sim 1.07$ requires that there must be significant sub-photospheric emission and/or an anisotropic explosion to explain the broadening for most GRB spectra.

We report the first analysis of data from AstroSat/LAXPC observations of Cygnus X-1 in January 2016. LAXPC spectra reveals that the source was in the canonical hard state, represented by a prominent thermal Comptonization component having a photon index of 1.8 and high temperature kT of electron > 60 keV along with weak reflection and possible disk emission. The power spectrum can be characterized by two broad lorentzian functions centered at 0.4 and 3 Hz. The r.m.s of the low frequency component decreases from 15% at around 4 keV to 10% at around 50 keV, while that of the high frequency one varies less rapidly from 13.5% to 11.5% in the same energy range. The time lag between the hard (20 to 40 keV) and soft (5 to 10 keV) bands varies in a step-like manner being nearly constant at 50 Milli-seconds from 0.3 to 0.9 Hz, decreasing to 8 Milli-seconds from 2 to 5 Hz and finally dropping to 2 Milli-seconds for higher frequencies. The time lags increase with energy for both the low and high frequency components. The event mode LAXPC data allows for flux resolved spectral analysis on a time-scale of 1 second, which clearly shows that the photon index increased from 1.72 to 1.80 as the flux increased by nearly a factor of two. We discuss the results in the framework of the fluctuation propagation model.

We study multi-wavelength variability of a bare Seyfert 1 galaxy Fairall~9 using \swift{} monitoring observations consisting of $165$ usable pointings spanning nearly two years and covering six UV/optical bands and X-rays. Fairall~9 is highly variable in all bands though the variability amplitude decreases from X-ray to optical bands. The variations in the X-ray and UV/optical bands are strongly correlated. Our reverberation mapping analysis using the {\tt JAVALIN} tool shows that the variation in the UV/optical bands lag behind the variations in the X-ray band by $\sim 2-10{\rm~days}$. These lag measurements strongly suggest that the optical/UV variations are mainly caused by variations in the X-rays, and the origin of most of the optical/UV emission is X-ray reprocessing. The observed lags are found to vary as $\tau\propto\lambda^{1.36\pm0.13}$, consistent with the prediction, $\tau\propto\lambda^{4/3}$, for X-ray reprocessing in a standard accretion disc. However, the predicted lags for an standard accretion disc with X-ray reprocessing using black hole mass ($M_{BH} \sim 2.6\times10^{8}~M_{\odot}$) estimated from the reverberation mapping of broad emission lines and accretion rate relative to the Eddington rate ($\dot{m_E} =0.02$) are shorter than the observed lags. These observations suggest that accretion disc in Fairall~9 is larger than that predicted by the standard disc model, and confirm similar findings in a few other Seyfert 1 galaxies such as NGC~5548.

We present a multi-wavelength temporal analysis of the blazar 3C 454.3 during the high $\gamma$-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 $\gamma$-ray band. The source exhibits diverse levels of variability with many flaring/active states in the continuously sampled $\gamma$-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 $\gamma$-ray, optical and $\gamma$-ray, and IR and optical to a state where $\gamma$-ray lags the IR/optical by $\sim$3 days. The results are consistent with the previous studies of the same during various $\gamma$-ray flaring and active episodes of the source. This consistency, in turn, suggests an extended localized emission region with almost similar conditions during various $\gamma$-ray activity states. On the other hand, the delay of $\gamma$-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 $\gamma$-ray being inverse Comptonized of IR/optical photons and external radiation field, respectively.

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

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

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.

We perform X-ray/ultraviolet (UV) spectral and X-ray variability studies of the radio-loud narrow-line Seyfert 1 (NLS1) galaxy RX J1633.3+4719 using XMM-Newton and Suzaku observations from 2011 and 2012. The 0.3-10 keV spectra consist of an ultrasoft component described by an accretion disc blackbody (kT_in = 39.6^{+11.2}_{-5.5} eV) and a power law due to the thermal Comptonization ({\Gamma} = 1.96^{+0.24}_{-0.31}) of the disc emission. The disc temperature inferred from the soft excess is at least a factor of 2 lower than that found for the canonical soft excess emission from radio-quiet NLS1s. The UV spectrum is described by a power law with photon index 3.05^{+0.56}_{-0.33}. The observed UV emission is too strong to arise from the accretion disc or the host galaxy, but can be attributed to a jet. The X-ray emission from RX J1633.3+4719 is variable with fractional variability amplitude $F_{\rm var}$=13.5$\pm1.0$ per cent. In contrast to radio-quiet active galactic nuclei (AGN), X-ray emission from the source becomes harder with increasing flux. The fractional rms variability increases with energy and the rms spectrum is well described by a constant disc component and a variable power-law continuum with the normalization and photon index being anticorrelated. Such spectral variability cannot be caused by variations in the absorption and must be intrinsic to the hot corona. Our finding of possible evidence for emission from the inner accretion disc, jet and hot corona from RX J1633.3+4719 in the optical to X-ray bands makes this object an ideal target to probe the disc-jet connection in AGN.

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.

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 $\gamma$-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 $\gamma$-ray and optical B-band. During the low level flux states, it is higher towards the peak of the spectral energy distribution, with $\gamma$-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 $\gamma$-ray background with the high state showing a relatively harder mean spectral index compared to the low state.

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 ($\eta \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 $\eta \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.

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.

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.

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/$\rho$ 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}^\beta (t-\tau_d)$. The fluctuating accretion rate drives the oscillations of the inner radius after a time delay $\tau_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.

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.

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.

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.

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.

Observational data, especially astrophysical data, is often limited by gaps in data that arises due to lack of observations for a variety of reasons. Such inadvertent gaps are usually smoothed over using interpolation techniques. However the smoothing techniques can introduce artificial effects, especially when non-linear analysis is undertaken. We investigate how gaps can affect the computed values of correlation dimension of the system, without using any interpolation. For this we introduce gaps artificially in synthetic data derived from standard chaotic systems, like the R{\"o}ssler and Lorenz, with frequency of occurrence and size of missing data drawn from two Gaussian distributions. Then we study the changes in correlation dimension with change in the distributions of position and size of gaps. We find that for a considerable range of mean gap frequency and size, the value of correlation dimension is not significantly affected, indicating that in such specific cases, the calculated values can still be reliable and acceptable. Thus our study introduces a method of checking the reliability of computed correlation dimension values by calculating the distribution of gaps with respect to its size and position. This is illustrated for the data from light curves of three variable stars, R Scuti, U Monocerotis and SU Tauri. We also demonstrate how a cubic spline interpolation can cause a time series of Gaussian noise with missing data to be misinterpreted as being chaotic in origin. This is demonstrated for the non chaotic light curve of variable star SS Cygni, which gives a saturated D$_{2}$ value, when interpolated using a cubic spline. In addition we also find that a careful choice of binning, in addition to reducing noise, can help in shifting the gap distribution to the reliable range for D$_2$ values.

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.

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 $\mu$m is marginally better correlated than the 24 $\mu$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.

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.

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.

We present a comprehensive flux resolved spectral analysis of the bright Narrow line Seyfert I AGNs, Mrk~335 and Ark~564 using observations by XMM-Newton satellite. The mean and the flux resolved spectra are fitted by an empirical model consisting of two Comptonization components, one for the low energy soft excess and the other for the high energy power-law. A broad Iron line and a couple of low energies edges are required to explain the spectra. For Mrk~335, the 0.3 - 10 keV luminosity relative to the Eddington value, L{$_{X}$}/L$_{Edd}$, varied from 0.002 to 0.06. The index variation can be empirically described as $\Gamma$ = 0.6 log$_{10}$ L{$_{X}$}/L$_{Edd}$ + 3.0 for $0.005 < L{_{X}}/L_{Edd} < 0.04$. At $ L_{{X}}/L_{Edd} \sim 0.04$ the spectral index changes and then continues to follow $\Gamma$ = 0.6 log$_{10}$ L$_{{X}}$/L$_{Edd}$ + 2.7, i.e. on a parallel track. We confirm that the result is independent of the specific spectral model used by fitting the data in the 3 - 10 keV band by only a power-law and an Iron line. For Ark~564, the index variation can be empirically described as $\Gamma$ = 0.2 log$_{10}$ L$_{{X}}$/L$_{Edd}$ + 2.7 with a significantly large scatter as compared to Mrk~335. Our results indicate that for Mrk~335, there may be accretion disk geometry changes which lead to different parallel tracks. These changes could be related to structural changes in the corona or enhanced reflection at high flux levels. There does not seem to be any homogeneous or universal relationship for the X-ray index and luminosity for different AGNs or even for the same AGN.

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.

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 $\epsilon n(\epsilon) \sim k log(\frac{\epsilon}{\epsilon_p}) $. We estimate the range of values for the parameters defining the EBL spectrum, $k$ and $\epsilon_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.

We present long term optical variability studies of bright X-ray sources in four nearby elliptical galaxies with {\it Chandra} Advanced CCD Imaging Spectrometer array (ACIS-S) and {\it Hubble Space Telescope (HST)} Advanced Camera for Surveys observations. Out of the 46 bright (X-ray counts $> 60$) sources that are in the common field of view of the {\it Chandra} and {\it HST} observations, 34 of them have potential optical counterparts, while the rest of them are optically dark. After taking into account of systematic errors, estimated using the field optical sources as reference, we find that four of the X-ray sources (three in NGC1399 and one in NGC1427) have variable optical counterparts at a high significance level. The X-ray luminosities of these source are $\sim 10^{38}$ $\rm ergs~s^{-1}$ and are also variable on similar time-scales. The optical variability implies that the optical emission is associated with the X-ray source itself rather than being the integrated light from a host globular cluster. For one source the change in optical magnitude is $> 0.3$, which is one of the highest reported for this class of X-ray sources and this suggests that the optical variability is induced by the X-ray activity. However, the optically variable sources in NGC1399 have been reported to have blue colours ($g - z > 1$). All four sources have been detected in the infra-red (IR) by {\it Spitzer} as point sources, and their ratio of $5.8$ to $3.6 \mu m$ flux are $> 0.63$ indicating that their IR spectra are like those of Active Galactic Nuclei (AGN). While spectroscopic confirmation is required, it is likely that all four sources are background AGNs. We find none of the X-ray sources having optical/IR colours different from AGNs, to be optically variable.

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.

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.

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.

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.

We make a detailed time resolved spectroscopy of bright long gamma ray bursts (GRBs) which show significant GeV emissions (GRB 080916C, GRB 090902B, and GRB 090926A). In addition to the standard Band model, we also use a model consisting of a blackbody and a power-law to fit the spectra. We find that for the latter model there are indications for an additional soft component in the spectra. While previous studies have shown that such models are required for GRB 090902B, here we find that a composite spectral model consisting of two black bodies and a power law adequately fit the data of all the three bright GRBs. We investigate the evolution of the spectral parameters and find several generic interesting features for all three GRBs, like a) temperatures of the black bodies are strongly correlated to each other, b) flux in the black body components are strongly correlated to each other, c) the temperatures of the black body trace the profile of the individual pulses of the GRBs, and d) the characteristics of the power law component like the spectral index and the delayed onset bear a close similarity to the emission characteristics in the GeV regions. We discuss the implications of these results to the possibility of identifying the radiation mechanisms during the prompt emission of GRBs.

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.

We analyze Suzaku and XMM-Newton data of the highly variable Seyfert 2, IRAS 18325-5926. The spectra of the source are well modeled as a primary component described as an absorbed power law and a secondary power law component which is consistent with being scattered emission from an on-axis extended highly ionized medium. We show that while the primary component varies on a wide range of timescales from $10^{4} - 10^{8}$ s, the scattered emission is variable only on timescales longer than $10^{5}$ s. This implies that the extent of the scattering medium is greater than $10^{16}$ cm. The ratio of the scattered to primary flux ($\sim 0.03$) implies a column density for the scattering medium to be $\sim 10^{23}$ cm$^{-2}$. We argue that for such a medium to be highly ionized it must be located less than $10^{17}$ cm from the X-ray source. Thus we localize the position and extent of scattering region to be $\sim$ a few $\times 10^{16}$ cm, with an average particle density of $\sim 10^{6}$ cm$^{-3}$. We consider the physical interpretation of these results and as an aside, we confirm the presence of a broad Iron line emission in both the {\it XMM-Newton} and {\it Suzaku} observations.

It is known that the IR and X-ray luminosities of Active Galactic Nuclei (AGN) are correlated with L_{IR} ~ L_{X}. Moreover, the IR flux ratio between the 5.8 and 3.6 \mu m bands is a good distinguishing characteristic of AGN or AGN-like behaviour. On the other hand, Galactic X-ray binaries (GXB) are under-luminous in the IR with L_{IR} << L_{X}. Since Ultra-luminous X-ray sources in nearby galaxies may be an intermediate class between AGN and GXB, it is interesting to study if their IR properties indicate which kind of objects they resemble. We use Spitzer IRAC images to identify mid-IR counterparts of bright X-ray sources, detected by Chandra in the elliptical galaxy NGC 1399. We find that for sources with AGN-like IR flux ratios, the IR luminosity strongly correlates with that in X-rays, L_{IR} ~ L_{X}, while for the others, there is no correlation between the two. Some of the former objects may be background AGN. If they are not strongly contaminated by background AGN, this result extends the IR-X-ray luminosity correlation down to L_{X} ~ 10^{39} ergs/s. We calculate their g-z colours and find that the bright X-ray sources with IR counterparts are typically blue in optical color. This is in contrast to typical X-ray sources, without IR counterparts which have predominantly red optical counterparts. We highlight the need for IR or optical spectra of these sources to distinguish background AGN and unveil the effect of the X-ray emission on the different environments of these systems.

We report on our search for very-long-term variability (weeks to years) in X-ray binaries (XRBs) in the giant elliptical galaxy M87. We have used archival Chandra imaging observations to characterise the long-term variability of 8 of the brightest members of the XRB population in M87. The peak brightness of some of the sources exceeded the ultra luminous X-ray source (ULX) threshold luminosity of ~ 10^{39} erg/s, and one source could exhibit dips or eclipses. We show that for one source, if it has similar modulation amplitude as in SS433, then period recoverability analysis on the current data would detect periodic modulations, but only for a narrow range of periods less than 120 days. We conclude that a dedicated monitoring campaign, with appropriately defined sampling, is essential if we are to investigate properly the nature of the long-term modulations such as those seen in Galactic sources.

We present broadband X-ray spectral study of two ultraluminous X-ray sources (ULXs), M81 X-6 and Holmberg IX X-1 based on Suzaku & XMM-Newton observations. We perform joint broadband spectral analysis of the brightest sources in the field, i.e. the two ULXs and the active galactic nucleus (AGN) in M81, and demonstrate that the X-ray spectra of the ULXs cut off at energies more than about 3keV with negligible contribution at high energies in the Suzaku HXD/PIN band. The 90% upper limit on the 10-30keV band luminosity of an underlying broadband power-law component is 3.5e38 ergs/s for M81 X-6 and 1.2e39ergs/s for Holmberg IX X-1. These limits are more than an order of magnitude lower than the bolometric (0.1-30 keV) luminosity of 6.8e39ergs/s for M81 X-6 and 1.9e40 ergs/s for Holmberg IX X-1. Our results confirm earlier indications of spectral cut-offs inferred from the XMM-Newton observations of bright ULXs and show that there isn't an additional high energy power-law component contributing significantly to the X-ray emission. The spectral form of the two ULXs are very different from those of Galactic black hole X-ray binaries (BHBs) or AGNs. This implies that the ULXs are neither simply scaled-up versions of stellar mass BHBs or scaled-down versions of AGNs.

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.

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.

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.

We analyse a large number ($> 500$) pointed RXTE observations of Cyg X-1 and model the spectrum of each one. A subset of the observations for which there is simultaneous reliable measure of the hardness ratio by the All Sky Monitor, shows that the sample covers nearly all the spectral shapes of Cyg X-1. The relative strength, width of the Iron line and the reflection parameter are in general correlated with the high energy photon spectral index $\Gamma$. This is broadly consistent with a geometry where for the hard state (low $\Gamma \sim 1.7$) there is a hot inner Comptonizing region surrounded by a truncated cold disk. The inner edge of the disk moves inwards as the source becomes softer till finally in the soft state (high $\Gamma > 2.2$) the disk fills the inner region and active regions above the disk produce the Comptonized component. However, the reflection parameter shows non-monotonic behaviour near the transition region ($\Gamma \sim 2$), suggestive of a more complex geometry or physical state of the reflector. Additionally, the inner disk temperature, during the hard state, is on the average higher than in the soft one, albeit with large scatter. These inconsistencies could be due to limitations in the data and the empirical model used to fit them. The flux of each spectral component is well correlated with $\Gamma$ which shows that unlike some other black hole systems, Cyg X-1 does not show any hysteresis behaviour. In the soft state, the flux of the Comptonized component is always similar to the disk one, which confirms that the ultra-soft state (seen in other brighter black hole systems) is not exhibited by Cyg X-1. The rapid variation of the Compton Amplification factor with $\Gamma$, naturally explains the absence of spectra with $\Gamma < 1.6$, despite a large number having $\Gamma \sim 1.65$.

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.

Origin of extragalactic gamma-ray background (EGRB) is still a matter of debate. EGRB can either have truly diffuse or unresolved discrete point sources origin. Majority of the Fermi and EGRET detected identified sources are blazar. So, they are expected to be a significant contributors to the EGRB. In order to estimate their contribution to the EGRB one needs to construct their luminosity functions. Here, we examine and construct the evolution and luminosity function of blazars using first LAT AGN catalog. We consider both pure luminosity and pure density evolution models for flat spectrum radio quasars (FSRQs). We also describe the methodology to estimate the contribution from misaligned blazars to the EGRB.

We undertake a spectral study of a sample of bright X-ray sources taken from six XMM-Newton fields at high galactic latitudes, where AGN are the most populous class. These six fields were chosen such that the observation had an exposure time more than 60 ksec, had data from the EPIC-pn detector in the full-Frame mode and lying at high galactic latitude $|b| > 25^o$ . The analysis started by fitting the spectra of all sources with an absorbed power-law model, and then we fitted all the spectra with an absorbed power-law with a low energy black-body component model.The sources for which we added a black body gave an F-test probability of 0.01 or less (i.e. at 99% confidence level), were recognized as sources that display soft excess. We perform a comparative analysis of soft excess spectral parameters with respect to the underlying power-law one for sources that satisfy this criterion. Those sources, that do not show evidence for a soft excess, based on the F-test probability at a 99% confidence level, were also fitted with the absorbed power-law with a low energy black-body component model with the black-body temperature fixed at 0.1 and 0.2 keV. We establish upper limits on the soft excess flux for those sources at these two temperatures. Finally we have made use of Aladdin interactive sky atlas and matching with NASA/IPAC Extragalactic Database (NED) to identify the X-ray sources in our sample. For those sources which are identified in the NED catalogue, we make a comparative study of the soft excess phenomenon for different types of systems.

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.

We present the results of a broadband spectral study of the transient Low Mass X-ray Binary Aql X-1 observed by Suzaku and Rossi X-ray Timing Explorer satellites. The source was observed during its 2007 outburst in the High/Soft (Banana) state and in the Low/Hard (Extreme Island) state. Both the Banana state and the Extreme Island state spectra are best described by a two component model consisting of a soft multi-colour blackbody emission likely originating from the accretion disk and a harder Comptonized emission from the boundary layer. Evidence for a hard tail (extending to ~50 keV) is found during the Banana state; this further (transient) component, accounting for atleast ~1.5% of the source luminosity, is modeled by a power-law. Aql X-1 is the second Atoll source after GX 13+1 to show a high energy tail. The presence of a weak but broad Fe line provides further support for a standard accretion disk extending nearly to the neutron star surface. The input photons for the Comptonizing boundary layer could either be the disk photons or the hidden surface of the star or both. The luminosity of the boundary layer is similar to the disk luminosity in the banana state and is about six times larger in the extreme island state. The temperature of the Comptonizing boundary layer changes from ~2 keV in the banana state to ~20 keV in the extreme island state.

Estimation of the black hole mass in bright X-ray sources of nearby galaxies is crucial to the understanding of these systems and their formation. However, the present allowed black hole mass range spans five order of magnitude (10Msun < M < 10^5 Msun) with the upper limit obtained from dynamical friction arguments. We show that the absence of a detectable optical counterpart for some of these sources, can provide a much more stringent upper limit. The argument is based only on the assumption that the outer regions of their accretion disks is a standard one. Moreover, such optically dark X-ray sources cannot be foreground stars or background active galactic nuclei, and hence must be accreting systems residing within their host galaxies. As a demonstration we search for candidates among the point-like X-ray sources detected with Chandra in thirteen nearby elliptical galaxies. We use a novel technique to search for faint optical counterparts in the HST images whereby we subtract the bright galaxy light based on isophotal modeling of the surface brightness. We show that for six sources with no detectable optical emission at the 3$-$sigma level, their black hole masses M_{BH} < 5000Msun. In particular, an ultra-luminous X-ray source (ULX) in NGC 4486 has M_{BH} < 1244 Msun. We discuss the potential of this method to provide stringent constraints on the black hole masses, and the implications on the physical nature of these sources.

We consider implications of a possible presence of a Thomson-thick, low-temperature, plasma cloud surrounding the compact object in the binary system Cyg X-3. The presence of such a cloud was earlier inferred from the energy-independent orbital modulation of the X-ray flux and the lack of high frequencies in its power spectra. Here, we study the effect of Compton scattering by the cloud on the X-ray energy and power spectra, concentrating on the hard spectral state. The process reduces the energy of the high-energy break/cut-off in the energy spectra, which allows us to determine the Thomson optical depth. This, together with the observed cut-off in the power spectrum, determines the size of the plasma to be 2x10^9 cm. At this size, the cloud will be in thermal equilibrium in the photon field of the X-ray source, which yields the cloud temperature of 3 keV, which refines the determination of the Thomson optical depth to 7. At these parameters, thermal bremsstrahlung emission of the cloud becomes important as well. The physical origin of the cloud is likely to be collision of the very strong stellar wind of the companion Wolf-Rayet star with a small accretion disc formed by the wind accretion. Our model thus explains the peculiar X-ray energy and power spectra of Cyg X-3.

GRS 1915+105 is a prominent black hole system exhibiting variability over a wide range of time scales and its observed light curves have been classified into 12 temporal states. Here we undertake a complete analysis of these light curves from all the states using various quantifiers from nonlinear time series analysis, such as, the correlation dimension (D_2), the correlation entropy (K_2), singular value decomposition (SVD) and the multifractal spectrum ($f(\alpha)$ spectrum). An important aspect of our analysis is that, for estimating these quantifiers, we use algorithmic schemes which we have proposed recently and tested successfully on synthetic as well as practical time series from various fields. Though the schemes are based on the conventional delay embedding technique, they are automated so that the above quantitative measures can be computed using conditions prescribed by the algorithm and without any intermediate subjective analysis. We show that nearly half of the 12 temporal states exhibit deviation from randomness and their complex temporal behavior could be approximated by a few (3 or 4) coupled ordinary nonlinear differential equations. These results could be important for a better understanding of the processes that generate the light curves and hence for modelling the temporal behavior of such complex systems. To our knowledge, this is the first complete analysis of an astrophysical object (let alone a black hole system) using various techniques from nonlinear dynamics.

(Abridged) Stellar mass black hole X-ray binaries exhibit X-ray spectral states which also have distinct and characteristic temporal properties. These states are believed to correspond to different accretion disk geometries. We present analysis of two long XMM-Newton observations of the Ultra-Luminous X-ray source (ULX) NGC 1313 X-1, which reveal that the system was in two different spectral/temporal states. With a count rate of 1.5 counts/s and a fractional variability amplitude of ~15%, the ULX was in a high flux and strongly variable state in March 2006. In October 2006, the count rate of the ULX had reduced by a factor of ~2 and the spectral shape was distinctly different with the presence of a soft component. No strong variability was detected during this low flux state with an upper limit on the amplitude < 3%. Moreover, the spectral properties of the two states implies that the accretion disk geometry was different for them. The low flux state is consistent with a model where a standard accretion disk is truncated at a radius of ~17 Schwarzschild radius around a ~200 Msun black hole. The inner hot region Comptonizes photons from the outer disk to give the primary spectral component. The spectrum of the high flux state is not compatible with such a geometry. Instead, it is consistent with a model where a hot corona covers a cold accretion disk and Comptonizes the disk photons. The variability as a function of energy is also shown to be consistent with the corona model. Despite these broad analogies with Galactic black hole systems, the spectral nature of the ULX is distinct in having a colder Comptonizing temperature (~2 keV) and higher optical depth (~15) than what is observed for the Galactic ones.

We present the results obtained by a detailed study of the extragalactic Z source, LMC X-2, using broad band Suzaku data and a large ($ \sim 750$ ksec) data set obtained with the proportional counter array (PCA) onboard RXTE. The PCA data allows for studying the complete spectral evolution along the horizontal, normal and flaring branches of the Z-track. Comparison with previous study show that the details of spectral evolution (like variation of Comptonizing electron temperature), is similar to that of GX 17+2 but unlike that of Cyg X-2 and GX 349+2. This suggests that Z sources are heterogeneous group with perhaps LMC X-2 and GX 17+2 being member of a subclass. However non monotonic evolution of the Compton y-parameter seems to be generic to all sources. The broad band {\it Suzaku} data reveals that the additional soft component of the source modelled as a disk blackbody emission is strongly preferred over one where it is taken to be a blackbody spectrum. This component as well as the temperature of seed photons do not vary when source goes into a flaring mode and the entire variation can be ascribed to the Comptonizing cloud. The bolometric unabsorbed luminosity of the source is well constrained to be $ \sim 2.23 \times 10^{38}$ ergs/sec which if the source is Eddington limited implies a neutron star mass of 1.6 M$_\odot$. We discuss the implications of these results.

We investigate the X-ray variability characteristics of hard X-ray selected AGNs (based on Swift/BAT data) in the soft X-ray band using the RXTE/ASM data. The uncertainties involved in the individual dwell measurements of ASM are critically examined and a method is developed to combine a large number of dwells with appropriate error propagation to derive long duration flux measurements (greater than 10 days). We also provide a general prescription to estimate the errors in variability derived from rms values from unequally spaced data. Though the derived variability for individual sources are not of very high significance, we find that, in general, the soft X-ray variability is higher than those in hard X-rays and the variability strengths decrease with energy for the diverse classes of AGN. We also examine the strength of variability as a function of the break time scale in the power density spectrum (derived from the estimated mass and bolometric luminosity of the sources) and find that the data are consistent with the idea of higher variability at time scales longer than the break time scale.

We report variability of the X-ray source, X-7, in NGC 6946, during a 60 ksec Chandra observation when the count rate decreased by a factor of ~1.5 in ~5000 secs. Spectral fitting of the high and low count rate segments of the light curve reveal that the simplest and most probable interpretation is that the X-ray spectra are due to disk black body emission with an absorbing hydrogen column density equal to the Galactic value of 2.1 X 10^{21} cm^{-2}. During the variation, the inner disk temperature decreased from ~0.29 to ~0.26 keV while the inner disk radius remained constant at ~6 X 10^8 cm. This translates into a luminosity variation from 3.8 to 2.8 X 10^{39} ergs cm^{-2} sec^{-1} and a black hole mass of ~400 solar masses. More complicated models like assuming intrinsic absorption and/or the addition of a power-law component imply a higher luminosity and a larger black hole mass. Even if the emission is beamed by a factor of ~5, the size of the emitting region would be > 2.7 X 10^8 cm implying a black hole mass > 180 solar masses. Thus, these spectral results provide strong evidence that the mass of the black hole in this source is definitely > 100 solar masses and more probably ~400 solar masses.

We develop a model of an accretion disc in which the variability induced at a given radius is governed by a damped harmonic oscillator at the corresponding epicyclic frequency. That variability induces both linear and non-linear responses in the locally emitted radiation. The total observed variability of a source is the sum of these contributions over the disc radius weighted by the energy dissipation rate at each radius. It is shown that this simple model, which effectively has only three parameters including the normalization, can explain the range of the power spectra observed from Cyg X-1 in the soft state. Although a degeneracy between the black hole mass and the strength of the damping does not allow a unique determination of the mass, we can still constrain it to <16--20 solar masses. We also show that our model preserves the observed linear rms-flux relationship even in the presence of the non-linear flux response.

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.

Data from {\it Chandra} observations of thirty nearby galaxies were analyzed and 365 X-ray point sources were chosen whose spectra were not contaminated by excessive diffuse emission and not affected by photon pile up. The spectra of these sources were fitted using two spectral models (an absorbed power-law and a disk blackbody) to ascertain the dependence of estimated parameters on the spectral model used. It was found that the cumulative luminosity function depends on the choice of the spectral model, especially for luminosities $> 10^{40}$ ergs/s. In accordance with previous results, a large number ($\sim 80$) of the sources have luminosities $> 10^{39}$ ergs/s (Ultra-Luminous X-ray sources) with indistinguishable average spectral parameters (inner disk temperature $\sim 1$ keV and/or photon index $\Gamma \sim 2$) with those of the lower luminosities ones. After considering foreground stars and known background AGN,we identify four sources whose minimum luminosity exceed $10^{40}$ ergs/s, and call them Extremely Luminous X-ray sources (ELX). The spectra of these sources are in general better represented by the disk black body model than the power-law one. These ELX can be grouped into two distinct spectral classes. Two of them have an inner disk temperature of $< 0.5$ keV and hence are called ``supersoft'' ELX, while the other two have temperatures $\gtrsim 1.3$ keV and are called ``hard'' ELX. The estimated inner disk temperatures of the supersoft ELX are compatible with the hypothesis that they harbor intermediate size black holes, which are accreting at $\sim 0.5$ times their Eddington Luminosity. The radiative mechanism for hard ELX, seems to be Inverse Comptonization, which in contrast to standard black holes systems, is probably saturated.

Emission spectra of hot accretion disks characteristic of advection dominated accretion flow (ADAF) models are investigated for comparison with the brightest ultra-luminous source, X-1, in the galaxy M82. If the spectral state of the source is similar to the low luminosity hard state of stellar mass black holes in our Galaxy, a fit to the {\it Chandra} X-ray spectrum and constraints from the radio and infrared upper limits, require a black hole mass in the range of $9 \times 10^4 - 5 \times 10^5 \msun$. Lower black hole masses ($\la 10^4 \msun$) are possible if M82 X-1 corresponds to the high luminosity hard state of Galactic black hole X-ray binary sources. Both of these spectrally degenerate hot accretion disk solutions lead to an intermediate mass black hole interpretation for M82 X-1. Since these solutions have different spectral variability with X-ray luminosity and predict different infrared emission, they can be distinguished by future off axis {\it Chandra} observations or simultaneous sensitive infrared detections.

In the inner regions of an accretion disk around a black hole, relativistic protons can interact with ambient matter to produce electrons, positrons and $\gamma$-rays. The resultant steady state electron and positron particle distributions are self-consistently computed taking into account Coulomb and Compton cooling, $e^-e^+$ pair production (due to $\gamma-\gamma$ annihilation) and pair annihilation. While earlier works used the diffusion approximation to obtain the particle distributions, here we solve a more general integro-differential equation that correctly takes into account the large change in particle energy that occur when the leptons Compton scatter off hard X-rays. Thus this formalism can also be applied to the hard state of black hole systems, where the dominant ambient photons are hard X-rays. The corresponding photon energy spectrum is calculated and compared with broadband data of black hole binaries in different spectral states. The results indicate that the $\gamma$-ray spectra ($E > 0.8$ MeV) of both the soft and hard spectral states and the entire hard X-ray/$\gamma$-ray spectrum of the ultra-soft state, could be due to $p-p$ interactions. These results are consistent with the hypothesis that there always exists in these systems a $\gamma$-ray spectral component due to $p-p$ interactions which can contribute between 0.5 to 10% of the total bolometric luminosty. The model predicts that {\it GLAST} would be able to detect black hole binaries and provide evidence for the presence of non-thermal protons which in turn would give insight into the energy dissipation process and jet formation in these systems.

We present an adaptation of the standard Grassberger-Proccacia (GP) algorithm for estimating the Correlation Dimension of a time series in a non subjective manner. The validity and accuracy of this approach is tested using different types of time series, such as, those from standard chaotic systems, pure white and colored noise and chaotic systems added with noise. The effectiveness of the scheme in analysing noisy time series, particularly those involving colored noise, is investigated. An interesting result we have obtained is that, for the same percentage of noise addition, data with colored noise is more distinguishable from the corresponding surrogates, than data with white noise. As examples for real life applications, analysis of data from an astrophysical X-ray object and human brain EEG, are presented.

Using non-linear time series analysis, along with surrogate data analysis, it is shown that the various types of long term variability exhibited by the black hole system GRS 1915+105, can be explained in terms of a deterministic non-linear system with some inherent stochastic noise. Evidence is provided for a non-linear limit cycle origin of one of the low frequency QPO detected in the source, while some other types of variability could be due to an underlying low dimensional chaotic system. These results imply that the partial differential equations which govern the magneto-hydrodynamic flow of the inner accretion disk, can be approximated by a small number ($\approx 3 -5$) of non-linear but {\it ordinary} differential equations. While this analysis does not reveal the exact nature of these approximate equations, they may be obtained in the future, after results of magneto-hydrodynamic simulation of realistic accretion disks become available.

The results of a spectral analysis, using {\it XMM-Newton} and {\it Chandra} data of the brightest ultra luminous X-ray source in the nearby galaxy M82, are presented. The spectrum of M82 X-1, was found to be unusually hard (photon spectral index $\Gamma \approx 1$) with a sharp cutoff at $\approx 6$ keV. Disk black body emission model requires a nonphysically high temperature. Instead, the spectrum is better described, with a lower reduced $\chi^2$, as emission due to nearly saturated Comptonization of photons in an optically thick ($\tau \approx 10-30$, depending on the geometry) plasma having a temperature $kT \approx 2$ keV. This is in contrast to the high energy spectra of other black hole systems, which are relatively steeper ($\Gamma > 1.5$) and hence are modeled as un-saturated thermal and/or non-thermal Comptonization of soft photons, in an optically thin ($\tau \approx 1$) high temperature plasma. An Iron line emission which is marginally resolved ($\sigma \sim 0.2$ keV) is required to fit the data. We argue that the standard geometries for the X-ray producing region, which are an optically thin inner disk or an uniform/patchy corona on top of a cold disk, are not applicable to this source. Alternatively, the geometry of the X-ray producing region could be a large sphere surrounding a cold accretion disk or an optically thick inner disk region which cools by bremsstrahlung self-Comptonization. For the latter scenario, such an inner disk region, whose effective optical depth to absorption is less than unity, is expected in the standard accretion disk theory for near Eddington accretion rates.

The dynamics of relativistically moving blobs ejected out of a central AGN, are considered. It is assumed that the collision between two blobs are completely inelastic, such that the bulk kinetic energy lost in the collision is used to energize electrons to relativistic energies via acceleration in internal shocks which are formed due to the collision. These high energy electrons which are produced on a time-scale corresponding to the collision time-scale, cool by radiative losses due to synchrotron and Inverse Compton processes. The model is applied to the radio/X-ray knots of several AGN. For three of these sources we have analyzed long (> 40 ksec) Chandra observations and report on constrains on the X-ray spectral indices. In the framework of this model the AGN are inferred to sporadically eject relativistic blobs on time-scales ranging from 10^{11} to 10^{12} secs for different sources. It is shown that the collision time-scales can be longer than the age of the knot and, hence, non-thermal electrons are continuously being injected into the system. This continuous injection, in contrast to an instantaneous one time injection, gives rise to a characteristic spectral break, rather than a high-energy cutoff in the spectrum.

Nearly $1.32 \times 10^6$ sec of non-zero count rate data were obtained from five years of {\rm RXTE} observations of the atoll source 4U 1636-53. The data was divided into 10309 segments of 128 sec each. The histogram of the number of segments as a function of count rate shows that the system can naturally be classified into four flux states. For each segment an automated search for kHz QPO was undertaken and the histograms of the maximum Leahy normalised power $P_{max}$ for the four flux states were created. These were fitted by probability distribution functions under the assumption that the intrinsic amplitude of the QPO has a normal distribution. The best fit distribution functions, showed that the fraction $f$ of the time the system exhibits a kHz QPO, decreases from near unity for the lowest flux state ($\approx 1400$ c/s) to zero for the highest one ($> 2550$ c/s) while the average amplitude remains nearly constant. Based on the best fit probability distribution functions, a threshold $P_{thre}$ was defined, such that 95% of the segments with $P_{max} > P_{thre}$ correspond to a real QPO signal. For segments selected by this criterion, the frequency versus count rate plot reproduced the earlier known ``parallel track'' variation with the tracks coinciding roughly with the flux states. It is shown that the gaps between the tracks are not caused by uneven sampling, but rather the QPO phenomenon is absent or weak when the system's flux level is intermediate between two flux states.

A modified non-linear time series analysis technique, which computes the correlation dimension $D_2$, is used to analyze the X-ray light curves of the black hole system GRS 1915+105 in all twelve temporal classes. For four of these temporal classes $D_2 $ saturates to $\approx 4-5$ which indicates that the underlying dynamical mechanism is a low dimensional chaotic system. Of the other eight classes, three show stochastic behavior while five show deviation from randomness. The light curves for four classes which depict chaotic behavior have the smallest ratio of the expected Poisson noise to the variability ($ < 0.05$) while those for the three classes which depict stochastic behavior is the highest ($ > 0.2$). This suggests that the temporal behavior of the black hole system is governed by a low dimensional chaotic system, whose nature is detectable only when the Poisson fluctuations are much smaller than the variability.

In the inner regions of an accretion disk around a black hole, the presence of non-thermal protons would via p-p collisions produce electrons, positrons and $\gamma$-rays. For such a scenario, the steady state electron-positron distribution is computed taking into account Compton cooling, $e^{-} e^{+}$ pair production (due to $\gamma-\gamma$ interactions) and pair annihilation. The resultant spectra has a broad feature around $1 - 10$ MeV which may be tested with observations by INTEGRAL and GLAST. Using the available OSSE data for the black hole system, GRS 1915+115 we illustrate that an upper limit on the non-thermal proton fraction can be obtained, which may put constraints on the acceleration process active in this system.

We consider the evolution of a spherically expanding plasma cloud, where there is continuous injection of non-thermal electrons. We compute the time dependent electron distribution and resultant photon spectra taking into account synchrotron, adiabatic and inverse Compton cooling. This model is different from previous works where, instead of a continuous injection of particles, a short injection period was assumed. We apply this model to the radio/optical knots in the large scale jets of AGN, detected in X-rays by {\it Chandra} and find that the overall broadband spectral features can be reproduced. It is shown that for some sources, constraints on the X-ray spectral index (by a longer {\it Chandra} observation) will be able to differentiate between the different models. This in turn will put a strong constraint on the acceleration mechanism active in these sources.

The surface flux (and the corresponding observed flux) of a geometrically thick ``funnel'' shaped disk is computed taking into account the radiation impinging on the surface from other parts of the disk. It is found that the ratio of the maximum apparent luminosity to the real luminosity of the disk is only a factor $\approx 5$ even when the opening angle of the disk is small ($\approx 15^o$). Thus, geometrically beamed emission from ``funnel'' shaped sub-Eddington disks around stellar mass black holes, cannot explain the Ultra-Luminous X-ray sources detected in nearby galaxies.

Two pseudo-Newtonian potentials, which approximate the angular and epicyclic frequencies of the relativistic accretion disk around rotating (and counter rotating) compact objects, are presented. One of them, the Logarithmically Modified Potential, is a better approximation for the frequencies while the other, the Second-order Expanded potential, also reproduces the specific energy for circular orbits in close agreement with the General Relativistic values. These potentials may be included in time dependent hydrodynamical simulations to study the temporal behavior of such accretion disks.

The vertical structure of optically thick accretion disks is investigated in the two-zone approximation. The disk is divided into an underlying disk and a corona, where the latter is defined as the upper surface layers for which the effective optical depth is unity. It is found that a significant part of the accretion flow (or dissipation rate) can take place in the corona if the scale height of the magnetic field is larger than that of the disk. The presence of such a dissipative corona leads to a modification in the topology of local disk solutions. For example, these solutions are found from local stability analysis to be both secularly and thermally stable, for accretion rates which are a factor $\approx$ four higher than those inferred from the stability of standard disk solutions. Thus, the applicability of optically thick disks with dissipative coronas are not as restrictive as disks without such coronas and can provide an attractive explanation for the origin of the soft spectral component observed in black hole X-ray binary systems.

We compute spectra from accretion disks around rapidly rotating neutron stars. The full effect of general relativity is considered for the structure calculation of the stars. We take into account the Doppler shift, gravitational redshift and light-bending effects in order to compute the observed spectra. To facilitate direct comparison with observations, a simple empirical function is presented which describes the numerically computed spectra well. This function can in principle be used to distinguish between the Newtonian spectra and the relativistic spectra. We also discuss the possibility of constraining neutron star's equation of state using our spectral models.

An empirical Compton up-scattering model is described which reproduces both the fractional amplitude (RMS) vs. energy and the soft time lags in the 830 Hz QPO observed in 4U1608-52 on Mar. 3, 1996. A combination of two coherent variations in the coronal and soft photon temperatures (with their relative contributions determined by enforcing energy conservation) gives rise to the QPO's energy dependent characteristics. All input parameters to the model, save a characteristic plasma size and the fraction of Comptonized photons impinging on the soft photon source, are derived from the time-averaged photon energy spectrum of the same observation. Fits to the fractional RMS and phase lag data for this kilohertz QPO imply that the spatial extent of the plasma is in the range from 4 to 15 km.

The effect of a possible hydrodynamical wind on the nature of hot accretion disk solutions is studied. It is found that the advection dominated branch, in the presence of a wind, maintains the self-similar solution for the disk structure with the temperature, $\theta \propto 1/r$, optical depth $\tau \propto r^{P-1/2}$ and accretion rate $\dot m \propto r^P$. Based on global solutions, cooling due to wind energy loss and advection are found to be equally important. For a wide range of viscosity and wind parameters, the temperature is about one-tenth of the virial value and $P \approx 0.9$, independent of the mass accretion rate and radius. In the context of cooling by unsaturated Comptonization of soft photons, solutions also exist where radiative cooling, advection, and wind cooling are important. In this case, wind regulated solutions are possible. Here, the radial dependence of the critical mass accretion rate above which solutions do not exist is unchanged from those solutions without winds. The wind/advection dominated solutions are locally unstable to a new type of instability called ``wind driven insbility'' where the presence of a wind causes the disk to be unstable to long wavelength perturbations of the surface density. The growth rate of this instability is inversely proportional to the ratio of the radiative cooling to the gravitational energy dissipation rates and it can grow on a timescale much longer than the viscous timescale in the disk for sufficiently small radiative cooling efficiencies.

We calculate the temperature profiles of (thin) accretion disks around rapidly rotating neutron stars (with low surface magnetic fields), taking into account the full effects of general relativity. We then consider a model for the spectrum of the X-ray emission from the disk, parameterized by the mass accretion rate, the color temperature and the rotation rate of the neutron star. We derive constraints on these parameters for the X-ray source Cygnus X-2 using the estimates of the maximum temperature in the disk along with the disk and boundary layer luminosities, using the spectrum inferred from the EXOSAT data. Our calculations suggest that the neutron star in Cygnus X-2 rotates close to the centrifugal mass-shed limit. Possible constraints on the neutron star equation of state are also discussed.

General relativistic spectra from accretion disks around rotating neutron stars in the appropriate space-time geometry for several different equation of state, spin rates and mass of the compact object have been computed. The analysis involves the computation of the relativistically corrected radial temperature profiles and the effect of Doppler and gravitational red-shifts on the spectra. Light bending effects have been omitted for simplicity. The relativistic spectrum is compared with the Newtonian one and it is shown that the difference between the two is primarily due to the different radial temperature profile for the relativistic and Newtonian disk solutions. To facilitate direct comparison with observations, a simple empirical function has been presented which describes the numerically computed relativistic spectra well. This empirical function (which has three parameters including normalization) also describes the Newtonian spectrum adequately. Thus the function can in principle be used to distinguish between the two. In particular, the best fit value of one of the parameters ($\beta$-parameter) $\approx 0.4$ for the Newtonian case, while it ranges from 0.1 to 0.35 for relativistic case depending upon the inclination angle, EOS, spin rate and mass of the neutron star. Constraining this parameter by fits to future observational data of X-ray binaries will indicate the effect of strong gravity in the observed spectrum.

Hydrodynamical winds from a spherical two-temperature plasma surrounding a compact object are constructed. The mass-loss rate is computed as a function of electron temperature, optical depth and luminosity of the sphere, the values of which can be constrained by the fitting of the spectral energy distributions for known X-ray binary systems. The sensitive dependence of the mass loss rate with these parameters leads to the identification of two distinct regions in the parameter space separating wind-dominated from non wind dominated systems. A critical optical depth, tau_c, as a function of luminosity and electron temperature, is defined which differentiates these two regions. Systems with optical depths significantly smaller than tau_c are wind-dominated. The results are applied to black hole candidate X-ray binary systems in the hard spectral state (Cyg X-1, GX 339-4 and Nova Muscae), and it is found that the inferred optical depth (tau) is similar to tau_c suggesting that they are wind regulated systems. On the other hand, for X-ray binary systems containing a neutron star (e.g., Cyg X-2) tau is much larger than tau_c indicating the absence of significant hydrodynamical winds.

We show that waves propagating in a transition disk can explain the short term temporal behavior of Cygnus X-1. In the transition disk model the spectrum is produced by saturated Comptonization within the inner region of the accretion disk where the temperature varies rapidly with radius. Recently, the spectrum from such a disk has been shown to fit the average broad band spectrum of this source better than that predicted by the soft-photon Comptonization model. Here, we consider a simple model where waves are propagating cylindrically symmetrical ly in the transition disk with a uniform propagation speed ($c_p$). We show that this model can qualitatively explain (a) the variation of the power spectral density (PSD) with energy, (b) the hard lags as a function of frequency and (c) the hard lags as a function of energy for various frequencies. Thus the transition disk model can explain the average spectrum and the short term temporal behavior of Cygnus X-1.

We reanalyze the ASCA and BeppoSAX data of MCG-6-30-15, using a double zone model for the Iron line profile. In this model, the X-ray source is located around $\approx 10$ Schwarzschild radius and the regions interior and exterior to the X-ray source produce the line emission. We find that this model fits the data with similar reduced $\chi^2$ as the standard single zone model. The best fit inclination angle of the source ($i \approx 10^o$) for the medium intensity ASCA data set is compatible with that determined by earlier modeling of optical lines. The observed variability of the line profile with intensity can be explained as variations of the X-ray source size. That several AGN with broad lines have the peak centroid near 6.4 keV can be explained within the framework of this model under certain conditions. We also show that the simultaneous broad band observations of this source by BeppoSAX rules out the Comptonization model which was an alternative to the standard inner disk one. We thereby strengthen the case that the line broadening occurs due to the strong gravitational influence of a Black Hole.

The broad Iron line in MCG-6-30-15 is fitted to the Comptonization model where line broadening occurs due to Compton down-scattering in a highly ionized optically thick cloud. These results are compared to the disk line model where the broadening is due to Gravitational/Doppler effects in the vicinity of a black hole. We find that both models fit the data well and it is not possible to differentiate between them by fitting only the ASCA data. The best fit temperature and optical depth of the cloud are found to be kT = 0.54 keV and $\tau = 4.0$ from the Comptonization model. This model further suggests that while the temperature can be assumed to be constant, the optical depth varies during the observation period. We emphasis an earlier conclusion that simultaneous broad band data ($3 - 50$ keV) can rule out (or confirm) the Comptonization model.

We model the UV/optical spectrum of the black hole binary Nova Muscae as a sum of black body emissions from the outer region of an accretion disk. We show for self-consistency that scattering effects in this region are not important. The black hole mass ($M \approx 6 M_\odot$), the inclination angle ($\mu \approx 0.5$) and the distance to the source ($D \approx 5$ kpc) have been constrained by optical observations during quiescence (Orosz et al. 1996). Using these values we find that the accretion rate during the peak was ${\dot M} \approx 8 \times 10^{19}$ g sec$^{-1}$ and subsequently decayed exponentially. We define a radiative fraction ($f$) to be the ratio of the X-ray energy luminosity to the total gravitational power dissipated for a keplerian accretion disk. We find that $f \approx 0.1$ and remains nearly constant during the Ultra-soft and Soft spectral states. Thus for these states, the inner region of the accretion disk is advection dominated. $f$ probably increased to $\approx 0.5$ during the Hard state and finally decreased to $\approx 0.03$ as the source returned to quiescence.

We show that the Iron K emission line from MCG-6-30-15 could be broadened due to Comptonization by a surrounding highly ionized cloud with radius $\sim 10^{14}$ cms. We calculate the temperature of the cloud to be $\sim0.21\kev$, provided a reasonable estimate of the UV flux is made. The X-ray/$\gamma$-ray emission observed from the source is compatible with this model. Such a cloud should be highly ionized and strong absorption edges are not expected from the source (Fabian et al 1995). For a $\onlyten{6}\msol$ black hole the size of the could corresponds to about 300 Schwarzschild radius. The intrinsic line could then be emitted far from the black hole and gravitational red-shift and Doppler effects would be negligible. If the black hole mass is much larger than $\onlyten{6}\msol$, gravitational/Doppler red-shifts would also contribute significantly to the broadening. We argue that the broad red wing observed in the source does not by itself imply emission from regions close ($R < 5 r_s$) to the black hole. However, Comptonization cannot produce a double peak. The presence of such a feature is a clear sign of inner disk emission influenced by gravitational and Doppler effects, perhaps broadened by the Comptonization. We note that simultaneous broad band (2-100 keV) study of this source can also reveal (or rule out) the presence of such a Comptonizing cloud.

We have used the transition disk model to fit the simultaneous broad band ($2-500$ keV) spectrum of Cygnus X-1 from OSSE and Ginga observations. In this model, the spectrum is produced by saturated Comptonization within the inner region of the accretion disk, where the temperature varies rapidly with radius. In an earlier attempt, we demonstrated the viability of this model by fitting the data from EXOSAT, XMPC balloon and OSSE observations, though these were not made simultaneously. Since the source is known to be variable, however, the results of this fit were not conclusive. In addition, since only once set of observations was used, the good agreement with the data could have been a chance occurrence. Here, we improve considerably upon our earlier analysis by considering four sets of simultaneous observations of Cygnus X-1, using an empirical model to obtain the disk temperature profile. The vertical structure is then obtained using this profile and we show that the analysis is self- consistent. We demonstrate conclusively that the transition disk spectrum is a better fit to the observations than that predicted by the soft photon Comptonization model. Since the temperature profile is obtained by fitting the data, the unknown viscosity mechanism need not be specified. The disk structure can then be used to infer the viscosity parameter $\alpha$, which appears to vary with radius and luminosity. This behavior can be understood if $\alpha$ depends intrinsically on the local parameters such as density, height and temperature. However, due to uncertainties in the radiative transfer, quantitative statements regarding the variation of $\alpha$ cannot yet be made.

The long-term transitions of the black hole candidate Cygnus X-1 (between the states gamma_1, gamma_2, and gamma_3) include the occasional appearance of a strong ~ MeV bump (gamma_1), whose strength appears to be anti-correlated with the continuum flux (~ 400 keV) due to the Compton upscattering of cold disk photons by the inner, hot corona. We develop a self-consistent disk picture that accounts naturally for these transitions and their corresponding spectral variations. We argue that the bump is due to the self-Comptonization of bremsstrahlung photons emitted predominantly near the plane of the corona itself. Our results suggest that a decrease by a factor of approx 2 in the viscosity parameter alpha is responsible for quenching this bump and driving the system to the gamma_2 state, whereas a transition from gamma_2 to gamma_3 appears to be induced by an increase of about 25 % in the accretion rate Mdot. In view of the fact that most of the transitions observed in this source seem to be of the gamma_2 to gamma_3 variety, we conclude that much of the long term gamma-ray spectral variability in Cygnus X-1 is due to these small fluctuations in Mdot. The unusual appearance of the gamma_1 state apparently reflects a change in the dissipative processes within the disk.