The aim of this PhD Thesis is to characterize a representative sample of Supergiant X-ray Binaries (SGXBs) formed by 4 sources: XTE J1855-026, a classical SGXB with long-term stable X-ray flux; AX J1841.0-0535 and AX J1845.0-0433, two supergiant fast X-ray transients (SFXTs) with the X-ray emission mostly dominated by flaring; and IGR J00370+6122, something in between these 2 sub-groups. The physical processes that produce these observable differences are still a matter of debate. In this PhD Thesis I performed a study of these 4 different systems to provide new data to constrain the models. This study consists of:(i) the determination of the orbital solution,(ii) a systematic study of the wind behavior along the orbit by the measure of Halpha variations,(iii) a model of stellar atmospheres of the donor star,(iv) establish whether there are X-ray flux variations modulated by the orbital period. The study of the wind shows that Halpha variations are dominated by intrinsic wind processes. The stellar atmospheres study shows that the supergiant stars that harbor these binaries have a higher projected rotational velocity, higher He abundance and higher N/C ratio than that of isolated supergiant stars. The results of this study show that the eccentricity of the binary does not have a simple correlation with the differences in the X-ray flux of the different sub-groups. Therefore, the idea of a more complex scenario is consolidated. The discovery of a new type of system, IGR J00370+6122, which properties do not fit in any of the established sub-groups, reinforces the idea of a continuum in the observed properties more than a strict classification of systems. Furthermore, I have developed a pipeline to reduce spectra of the FRODOSpec spectrograph at the Liverpool Telescope optimized for the reduction of the red spectra of the obscured supergiants that we find in these SGXBs.

Non-thermal radiation has been predicted within bow shocks around runaway stars by recent theoretical works. We present X-ray observations towards the runaway stars $\zeta$ Oph (Chandra and Suzaku) and BD+43$^{\circ}$3654 (XMM-Newton) to search for the presence of non-thermal X-ray emission. We found no evidence of non-thermal emission spatially coincident with the bow shocks, nonetheless, diffuse emission is detected in the vicinity of $\zeta$ Oph. After a careful analysis of its spectral characteristics we conclude that this emission has a thermal nature with a plasma temperature of $T \approx 2 \times10^{6}$ K. The cometary shape of this emission seems to be in line with recent predictions of radiation-hydrodynamic models of runaway stars. The case of BD+43$^{\circ}$3654 is puzzling as non-thermal emission has been reported in a previous work for this source.

BD+6073 is the optical counterpart of the X-ray source IGR J00370+6122, a probable accretion-powered X-ray pulsar. The X-ray light curve of this binary system shows clear periodicity at 15.7 d, which has been interpreted as repeated outbursts around the periastron of an eccentric orbit. We obtained high-resolution spectra of BD+6073 at different epochs. We used the FASTWind code to generate a stellar atmosphere model to fit the observed spectrum and obtain physical magnitudes. The synthetic spectrum was used as a template for cross-correlation with the observed spectra to measure radial velocities. The radial velocity curve provided an orbital solution for the system. We have also analysed the RXTE/ASM and Swift/BAT light curves to confirm the stability of the periodicity. BD +6073 is a BN0.7 Ib low-luminosity supergiant located at an approximate distance of 3.1 kpc, in the CasOB4 association. We derive Teff=24000 K and log gc=3.0, and chemical abundances consistent with a moderately high level of evolution. The spectroscopic and evolutionary masses are consistent at the 1 sigma level with a mass of 15 solar masses. The recurrence time of the X-ray flares is the orbital period of the system. The NS is in a high eccentricity (e=0.56) orbit, and the X-ray emission is strongly peaked around orbital phase 0.2, though the observations are consistent with some level of X-ray activity happening at all orbital phases. The X-ray behaviour of IGR J00370+6122 is reminiscent of intermediate SFXTs, though its peak luminosity is rather low. The orbit is somewhat wider than those of classical persistent supergiant X-ray binaries, which, combined with the low luminosity of the mass donor, explains the low X-ray luminosity. IGR J00370+6122 will likely evolve towards a persistent supergiant system, highlighting the evolutionary connection between different classes of wind-accreting X-ray sources.

We aim both to complement the existing data on the spin history of the peculiar accreting X-ray pulsar GX 1+4 with more past and current data from BeppoSAX, INTEGRAL, and Fermi and to interpret the evolution in the framework of accretion theory. We used source light curves obtained from BeppoSAX/WFC and INTEGRAL/ISGRI to derive pulse periods using an epoch-folding analysis. Fermi/GBM data were analyzed by fitting a constant plus a Fourier expansion to background-subtracted rates, and maximizing the Y2 statistic. We completed the sample with hard X-ray light curves from Swift/BAT. The data were checked for correlations between flux and changes of the pulsar spin on different timescales. The spin-down of the pulsar continues with a constant change in frequency, i.e., an apparently accelerating change in the period. Over the past three decades, the pulse period has increased by about ~50%. Short-term fluctuations on top of this long-term trend do show anti-correlation with the source flux. Possible explanations of the observed long-term frequency and its dependence on flux are discussed.

We report on the long-term evolution of the spin period of the symbiotic X-ray pulsar GX 1+4 and a possible interpretation within a model of quasi-spherical accretion. New period measurements from BeppoSAX, INTEGRAL and Fermi observations have been combined with previously published data from four decades of observations. During the 1970s GX 1+4 was spinning up with the fastest rate among the known X-ray pulsars at the time. In the mid 1980s it underwent a change during a period of low X-ray ux and started to spin down with a rate similar in magnitude to the previous spin up rate. The spin period has changed from ~110 s to ~160 s within the last three decades. Our results demonstrate that the overall spin down trend continues and is stronger than ever. We compare the observations with predictions from a model assuming quasi-spherical accretion from the slow wind of the M giant companion.

We report on the long-term evolution of the spin period of the symbiotic X-ray pulsar GX 1+4 and a possible interpretation within a model of quasi-spherical accretion. New period measurements from BeppoSAX/WFC, INTEGRAL/ISGRI and Fermi/GBM observations have been combined with previously published data from four decades of observations. During the 1970's GX 1+4 was spinning up with the fastest rate among the known X-ray pulsars at the time. In the mid 1980's it underwent a change during a period of low X-ray flux and started to spin down with a rate similar in magnitude to the previous spin up rate. The spin period has changed from ~110 s to ~160 s within the last three decades. Our results demonstrate that the overall spin down trend continues and is stronger than ever. We compare the observations with predictions from a model assuming quasi-spherical accretion from the slow wind of the M giant companion.

The properties of wind accretion in symbiotic X-ray binaries (SyXBs) consisting of red-giant and magnetized neutron star (NS) are discussed. The spin-up/spin-down torques applied to NS are derived based on a hydrodynamic theory of quasi-spherical accretion onto magnetized NSs. In this model, a settling subsonic accretion proceeds through a hot shell formed around the NS magnetosphere. The accretion rate onto the NS is determined by the ability of the plasma to enter the magnetosphere.Due to large Reynolds numbers in the shell, the interaction of the rotating magnetosphere with plasma initiates a subsonic turbulence. The convective motions are capable of carrying the angular momentum through the shell. We carry out a population synthesis of SyXBs in the Galaxy with account for the spin evolution of magnetized NS. The Galactic number of SyXBs with bright (M_v<1) low-mass red-giant companion is found to be from \sim 40 to 120, and their birthrate is \sim 5\times 10^{-5}-10^{-4} per year. According to our model, among known SyXBs, Sct X-1 and IRXS J180431.1-273932 are wind-fed accretors. GX 1+4 lies in the transition from the wind-fed SyXBs to SyXBs in which the giants overflow their Roche lobe. The model successfully reproduces very long NS spins (such as in IGR J16358-4724 and 4U 1954+31) without the need to invoke very strong magnetic fields.