This thesis presents the results of a comprehensive multi-wavelength study of the high mass X-ray binary 2S 0114+650. This enigmatic source has previously been proposed to be the first in a new class of super-slow X-ray pulsars, containing a neutron star revolving once every 2.7 h. The 11.6 d orbital period of this system has been well established in both X-ray and optical wavelengths. During the initial stages of the research presented in this thesis we discovered an additional 30.7 d ???super-orbital??? modulation in the X-ray emission from this source. While similar periodicities seen in other X-ray binaries are commonly attributed to the precession of a warped accretion disc, the observational evidence suggests the absence of such a disc in 2S 0114+650. The purpose of this project is thus to determine the nature of the super-orbital modulation and to better constrain the astrophysical parameters of this system. To investigate the long-term variability we analysed ~8.5 yr of archived data from the Rossi X-ray Timing Explorer space telescope. The problem of the spurious ~24 h periods in this data was solved as a by-product of these studies. Follow-up pointed observations were obtained with this satellite in order to examine the spectral and temporal behaviour over the spin, orbital and super-orbital timescales. Independent confirmation of the super-orbital modulation was performed using ~2 yr of data from the INTEGRAL satellite obtained during a long-term monitoring campaign of the Cassiopeia region. The evolution of the spin, orbital and super-orbital periods over ~10 yr was examined using archived data from the Rossi X-ray Timing Explorer satellite. Radio observations were performed with the Giant Meterwave Radio Telescope to search for any radio emission associated with this source and to determine whether it is variable over the known periodicities. Near infrared observations were performed with the Mt Abu telescope to determine whether a Be star nature can be ruled out for the optical component. Finally, a statistical analysis of the properties of the confirmed super-orbital X-ray binaries was performed in order to search for commonalities between these systems.

This thesis presents the results of a comprehensive multi-wavelength study of the high mass X-ray binary 2S 0114+650. This enigmatic source has previously been proposed to be the first in a new class of super-slow X-ray pulsars, containing a neutron star revolving once every 2.7 h. The 11.6 d orbital period of this system has been well established in both X-ray and optical wavelengths. During the initial stages of the research presented in this thesis we discovered an additional 30.7 d ???super-orbital??? modulation in the X-ray emission from this source. While similar periodicities seen in other X-ray binaries are commonly attributed to the precession of a warped accretion disc, the observational evidence suggests the absence of such a disc in 2S 0114+650. The purpose of this project is thus to determine the nature of the super-orbital modulation and to better constrain the astrophysical parameters of this system. To investigate the long-term variability we analysed ~8.5 yr of archived data from the Rossi X-ray Timing Explorer space telescope. The problem of the spurious ~24 h periods in this data was solved as a by-product of these studies. Follow-up pointed observations were obtained with this satellite in order to examine the spectral and temporal behaviour over the spin, orbital and super-orbital timescales. Independent confirmation of the super-orbital modulation was performed using ~2 yr of data from the INTEGRAL satellite obtained during a long-term monitoring campaign of the Cassiopeia region. The evolution of the spin, orbital and super-orbital periods over ~10 yr was examined using archived data from the Rossi X-ray Timing Explorer satellite. Radio observations were performed with the Giant Meterwave Radio Telescope to search for any radio emission associated with this source and to determine whether it is variable over the known periodicities. Near infrared observations were performed with the Mt Abu telescope to determine whether a Be star nature can be ruled out for the optical component. Finally, a statistical analysis of the properties of the confirmed super-orbital X-ray binaries was performed in order to search for commonalities between these systems.

More than half a decade of X-ray astronomy with various balloon borne and space orbiting X-ray instruments, have led to discoveries and detailed studies of X-ray binaries. An important property of X-ray binaries is intensity variations of different magnitudes in a wide range of timescales from milliseconds (quasi-periodic oscillations, millisecond pulsations), to a few weeks (orbital and super-orbital modulations) or longer (outbursts etc). In this thesis, different types of variabilities of X-ray binaries are considered in X-ray binary population studies and to investigate certain aspects of some individual systems. In Chapter 1, we provide an introduction to various types of variabilities seen in different classes of X-ray binaries. We mention in detail the various periodic and aperiodic variabilities seen in X-ray binaries. In Chapter 2, we describe, in some detail, the various X-ray all sky monitors and X-ray observatories, data from which has been utilized in the work carried out in this thesis. We also describe the various data analysis techniques that we have used. The rest of the thesis is divided into two major sections: Variability studies of indi-vidual systems and X-ray binary population studies Variability studies of individual systems In Chapter 3, we report results from an investigation of energy resolved orbital in-tensity pro les and from exhaustive orbital phase resolved spectroscopic measurements of GX 301{2 with MAXI{GSC . The orbital variation of the spectral parameters, es-pecially the relation between the equivalent width of Fe line and the column density of absorbing matter are then utilized to examine the models for the mode of accretion onto the neutron star in GX 301{2: circumstellar disk model by Pravdo & Ghosh (2001), and the accretion stream model by Leahy & Kostka (2008). A very large equivalent width of the iron line along with a small value of the column density in the orbital phase range 0.10-0.30 after the periastron passage indicates an asymmetry in the distribution of the matter around the neutron star, strongly favoring the accretion stream model by Leahy & Kostka (2008). Presence of an eclipse in an X-ray binary can be useful in determining orbital param-eters like inclination and in estimating the orbital evolution by eclipse timing method, which is reported in Chapter 4. For the HMXB system IGR J16393{4643, we found a short eclipse in the Swift{BAT light-curve and utilized it to constrain the orbital in-clination of the system. We have also studied, for the rst time, broad-band pulsation and spectral characteristic of the system with a Suzaku observation, showing sub-orbital intensity variations. For another eclipsing and non-pulsing HMXB 4U 1700{37, the orbital evolution is studied using mid-eclipse times from observations with narrow eld instruments as well as from long term light-curves of X-ray all sky monitors. The orbital period decay rate is estimated to be 5 10 7 /yr, an order slower than a previous measurement by Rubin et al.(1996). Since no pulsations are detected in this system, it is difficult to estimate its orbital parameters, especially its eccentricity. Using mid-eclipse times from 10 years of Swift{BAT data, we have independently constrained the eccentricity of the binary system. X-ray binary population studies In Chapter 5, we report results from an analysis of the 16 years light-curves of X-ray binaries in 2-10 keV energy band of RXTE{ASM , used to construct the differential and integral probability distributions of count-rates. These distributions are then employed to construct multiple snapshots of X-ray binary luminosity functions of the Milky Way instead of averaging the luminosities, an improvement over previous analysis by Grimm et al. (2002). We found that the averaged luminosities of highly variable X-ray binaries do not represent their true positions in XLFs and the variability of X-ray binaries do indeed signi cantly affect the luminosity functions. In Chapter 6, the measurements of the averaged spectra of X-ray binaries using MAXI{GSC data are reported and are used for constructing the composite X-ray spec-trum. These composite X-ray binary spectra are useful in estimating the contribution of X-ray binaries in extra-galactic SEDs constructed from the simultaneous Chandra / XMM{Newton and NuSTAR observations of these galaxies. These SEDs will also serve as a useful input in estimating the contribution of X-ray binary heating at high redshift IGM during the Epoch of Re-ionization. In Chapter 7, we summarize the main conclusions of the work carried out in this thesis and discuss some future prospects related to this thesis.