Source author record

P. Reig

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

39works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

39 published item(s)

preprint2022arXiv

Average bolometric corrections and optical to X-ray flux measurements as a function of accretion rate for X-ray binaries

In this paper we use an RXTE library of spectral models from 10 black-hole and 9 pulsar X-ray binaries, as well as model spectra available in the literature from 13 extra-galactic Ultra-luminous X-ray sources (ULXs). We compute average bolometric corrections (BC=$\mathrm{L_{band}/L_{bol}}$) for our sample as a function of different accretion rates. We notice the same behaviour between black-hole and pulsars BCs only when ULX pulsars are included. These measurements provide a picture of the energetics of the accretion flow for an X-ray binary based solely on its observed luminosity in a given band. Moreover it can be a powerful tool in X-ray binary population synthesis models. Furthermore we calculate the X-ray (2-10 keV) to optical (V-band) flux ratios at different Eddington ratios for the black-hole X-ray binaries in our sample. This provides a metric of the maximum contribution of the disk to the optical emission of a binary system and better constraints on its nature (donor type etc). We find that the optical to X-ray flux ratio shows very little variation as a function of accretion rate, but testing for different disk geometries scenarios we find that the optical contribution of the disk increases as the $p$ value decreases ($T(r)\sim r^{-p}$). Moreover observational data are in agreement with a thicker disk scenario ($p<0.65$), which could also possibly explain the lack of observed high-inclination systems.

preprint2021arXiv

A new automated tool for the spectral classification of OB stars

(abridged) We develop a tool for the automated spectral classification of OB stars according to their sub-types. We use the regular Random Forest (RF) algorithm, the Probabilistic RF (PRF), and we introduce the KDE-RF method which is a combination of the Kernel-Density Estimation and the RF algorithm. We train the algorithms on the Equivalent Width (EW) of characteristic absorption lines (features) measured in high-quality spectra from large Galactic (LAMOST,GOSSS) and extragalactic surveys (2dF,VFTS) with available spectral-types and luminosity classes. We find that the overall accuracy score is $\sim$70% with similar results across all approaches. We show that the full set of 17 spectral lines is needed to reach the maximum performance per spectral class. We apply our model in other observational data sets providing examples of potential application of our classifier on real science cases. We find that it performs well for both single massive stars and for the companion massive stars in Be X-ray binaries. In addition, we propose a reduced 10-features scheme that can be applied to large data sets with lower S/N. The similarity in the performances of our models indicates the robustness and the reliability of the RF algorithm when it is used for the spectral classification of early-type stars. The score of $\sim$70% is high if we consider (a) the complexity of such multi-class classification problems, (b) the intrinsic scatter of the EW distributions within the examined spectral classes, and (c) the diversity of the training set since we use data obtained from different surveys with different observing strategies. In addition, the approach presented in this work, is applicable to data of different quality and of different format (e.g.,absolute or normalized flux) while our classifier is agnostic to the Luminosity Class of a star and, as much as possible, metallicity independent.

preprint2021arXiv

Illumination of the accretion disk in black hole binaries: An extended jet as the primary source of hard X-rays

The models that seek to explain the reflection spectrum in black hole binaries usually invoke a point-like primary source of hard X-rays. This source illuminates the accretion disk and gives rise to the discrete (lines) and continuum-reflected components. The main goal of this work is to investigate whether the extended, mildly relativistic jet that is present in black hole binaries in the hard and hard-intermediate states is the hard X-ray source that illuminates the accretion disk. We use a Monte Carlo code that simulates the process of inverse Compton scattering in a mildly relativistic jet. Blackbody photons from the thin accretion disk are injected at the base of the jet and interact with the energetic electrons that move outward. Despite the fact that the jet moves away from the disk at a mildly relativistic speed, we find that approximately $15-20$\% of the input soft photons are scattered back toward the accretion disk. The vast majority of the Comptonized, back-scattered photons escape very close to the black hole ($h\lesssim 6 r_g$, where $r_g$ is the gravitational radius), but a non-negligible amount escape at a wide range of heights. At high heights, $h\sim 500-2000\,r_g$, the distribution falls off rapidly. The high-height cutoff strongly depends on the width of the jet at its base and is almost insensitive to the optical depth. The disk illumination spectrum is softer than the direct jet spectrum of the radiation that escapes in directions that do not encounter the disk. We conclude that an extended jet is an excellent candidate source of hard photons in reflection models.

preprint2021arXiv

X-ray spectral-timing variability of 1A 0535+262 during the 2020 giant outburst

The Be/X-ray binary 1A 0535+262 underwent a giant X-ray outburst in November 2020, peaking at $\sim1\times 10^{38}$ erg s$^{-1}$ (1--100 keV, 1.8 kpc), the brightest outburst recorded for this source so far. Our goal is to search for patterns of correlated spectral and timing behavior that can be used to characterize the accretion states in hard X-ray transient pulsars. We have studied the evolution of the spectral continuum emission using hardness-intensity diagrams and the aperiodic variability of the source. The hardness-intensity diagram displays three distinct branches that can be identified with different accretion regimes. The characteristic frequency of the noise components correlates with the luminosity. Our observations cover the highest end of this correlation, at luminosities not previously sampled. We have found evidence for a flattening of the correlation at those high luminosities, which might indicate that the accretion disk reached the closest distance from the neutron star surface during the peak of the outburst. We also find evidence for hysteresis in the spectral and timing parameters: at the same luminosity level, the spectrum is harder and the characteristic noise frequency larger during the rise than during the decay of the outburst.

preprint2020arXiv

Inclination Dependence of The Time-Lag -- Photon-Index Correlation in BHXRBs and its Explanation with a Simple Jet Model

Recently, we reported an observational correlation between a) the time-lag of the hard (9 - 15 keV) with respect to the soft (2 - 5 keV) X-ray photons in black-hole X-ray binaries (BHXRBs) and b) the power-law photon index $Γ$ of the X-ray spectrum. This was physically explained with a simple jet model, i.e., a model where the Comptonization (the Compton upscattering of soft photons) happens in the jet. Here, we report the inclination dependence of this correlation, which we also explain with our jet model. Photons that emerge at different polar angles from the jet axis have different spectra and different time-lags. Because of this, we can explain quantitatively the type-B QPOs of GX 339-4 as resulting from a precessing jet.

preprint2020arXiv

Optical counterpart to Swift J0243.6+6124

Swift J0243.6+6124 was discovered during a giant X-ray outburst in October 2017. While there are numerous studies in the X-ray band, very little is known about the optical counterpart. We have performed an spectral and photometric analysis of the optical counterpart of this intriguing source. We find that the optical counterpart to Swift J0243.6+6124 is a V = 12.9, O9.5Ve star, located at a distance of $\sim5$ kpc. The optical extinction in the direction of the source is $A_V=3.6$ mag. The rotational velocity of the O-type star is 210 km s$^{-1}$. The long-term optical variability agrees with the growth and subsequent dissipation of the Be circumstellar disk after the giant X-ray outburst. The optical and X-ray luminosity are strongly correlated during the outburst, suggesting a common origin. We did not detect short-term periodic variability that could be associated with nonradial pulsations from the Be star photosphere.

preprint2019arXiv

Full orbital solution for the binary system in the northern Galactic disc microlensing event Gaia16aye

Gaia16aye was a binary microlensing event discovered in the direction towards the northern Galactic disc and was one of the first microlensing events detected and alerted to by the Gaia space mission. Its light curve exhibited five distinct brightening episodes, reaching up to I=12 mag, and it was covered in great detail with almost 25,000 data points gathered by a network of telescopes. We present the photometric and spectroscopic follow-up covering 500 days of the event evolution. We employed a full Keplerian binary orbit microlensing model combined with the motion of Earth and Gaia around the Sun to reproduce the complex light curve. The photometric data allowed us to solve the microlensing event entirely and to derive the complete and unique set of orbital parameters of the binary lensing system. We also report on the detection of the first-ever microlensing space-parallax between the Earth and Gaia located at L2. The properties of the binary system were derived from microlensing parameters, and we found that the system is composed of two main-sequence stars with masses 0.57$\pm$0.05 $M_\odot$ and 0.36$\pm$0.03 $M_\odot$ at 780 pc, with an orbital period of 2.88 years and an eccentricity of 0.30. We also predict the astrometric microlensing signal for this binary lens as it will be seen by Gaia as well as the radial velocity curve for the binary system. Events such as Gaia16aye indicate the potential for the microlensing method of probing the mass function of dark objects, including black holes, in directions other than that of the Galactic bulge. This case also emphasises the importance of long-term time-domain coordinated observations that can be made with a network of heterogeneous telescopes.

preprint2016arXiv

Accretion regimes in the X-ray pulsar 4U 1901+03

The source 4U 1901+03 is a high-mass X-ray pulsar than went into outburst in 2003. Observation performed with the Rossi X-ray Timing Explorer showed spectral and timing variability, including the detection of flares, quasi-periodic oscillations, complex changes in the pulse profiles, and pulse phase dependent spectral variability. We re-analysed the data covering the 2003 X-ray outburst and focused on several aspects of the variability that have not been discussed so far. These are the 10 keV feature and the X-ray spectral states and their association with accretion regimes, including the transit to the propeller state at the end of the outburst. We find that 4U 1901+03 went through three accretion regimes over the course of the X-ray outburst. At the peak of the outburst and for a very short time, the X-ray flux may have overcome the critical limit that marks the formation of a radiative shock at a certain distance above the neutron star surface. Most of the time, however, the source is in the subcritical regime. Only at the end of the outburst, when the luminosity decreased below ~10^{36} (d/10 kpc)^2 erg/s, did the source enter the propeller regime. Evidence for the existence of these regimes comes from the pulse profiles, the shape of the hardness-intensity diagram, and the correlation of various spectral parameters with the flux. The 10 keV feature appears to strongly depend on the X-ray flux and on the pulse phase, which opens the possibility to interpret this feature as a cyclotron line.

preprint2016arXiv

Long-term optical variability of high-mass X-ray binaries. II. Spectroscopy

We present the results of our monitoring program to study the long-term variability of the Halpha line in high-mass X-ray binaries. We have carried out the most complete optical spectroscopic study of the global properties of high-mass X-ray binaries so far with the analysis of more than 1100 spectra of 20 sources. Our aim is to characterise the optical variability timescales and study the interaction between the neutron star and the accreting material. Our results can be summarised as follows: i) we find that Be/X-ray binaries with narrow orbits are more variable than systems with long orbital periods, ii) we show that a Keplerian distribution of the gas particles provides a good description of the disks in Be/X-ray binaries, as it does in classical Be stars, iii) a decrease in the Halpha equivalent width is generally observed after major X-ray outbursts, iv) we confirm that the Halpha equivalent width correlates with disk radius, v) while systems with supergiant companions display, multi-structured profiles, most of the Be/X-ray binaries show at some epoch double-peak asymmetric profiles, indicating that density inhomogeneities is a common property in the disk of Be/X-ray binaries, vi) the profile variability (V/R ratio) timescales are shorter and the Halpha equivalent width are smaller in Be/X-ray binaries than in isolated Be stars, and vii) we provide new evidence that the disk in Be/X-ray binaries is on average denser than in classical Be stars.

preprint2016arXiv

On the origin of the hard X-ray tail in neutron-star X-ray binaries

Neutron star X-ray binaries emit a compact, optically thick, relativistic radio jet during low-luminosity, usually hard states, as Galactic black-hole X-ray binaries do. When radio emission is bright, a hard power-law tail without evidence for an exponential cutoff is observed in most systems. We have developed a jet model that explains many spectral and timing properties of black-hole binaries in the states where a jet is present. Our goal is to investigate whether our jet model can reproduce the hard tail, with the correct range of photon index and the absence of a high-energy cutoff, in neutron-star X-ray binaries. We have performed Monte Carlo simulations of the Compton upscattering of soft, accretion-disk or boundary layer photons, in the jet and computed the emergent energy spectra, as well as the time lag of hard photons with respect to softer ones as a function of Fourier frequency. We demonstrate that our jet model explains the observed power-law distribution with photon index in the range 1.8-3. With an appropriate choice of the parameters, the cutoff expected from Comptonization is shifted to energies above ~300 keV, producing a pure power law without any evidence for a rollover, in agreement with the observations. Our results reinforce the idea that the link between the outflow (jet) and inflow (disk) in X-ray binaries does not depend on the nature of the compact object, but on the process of accretion. Furthermore, we address the differences of jets in black-hole and neutron-star X-ray binaries and predict that the break frequency in the spectral energy distribution of neutron-star X-ray binaries, as a class, will be lower than that of black-hole binaries.

preprint2016arXiv

Optical linear polarization of 74 white dwarfs with the RoboPol polarimeter

We present the first linear polarimetric survey of white dwarfs (WDs). Our sample consists of WDs of DA and DC spectral types in the SDSS r magnitude range from 13 to 17. We performed polarimetric observations with the RoboPol polarimeter attached to the 1.3-m telescope at the Skinakas Observatory. We have 74 WDs in our sample, of which almost all are low polarized WDs with polarization degree (PD) smaller than 1%, while only 2 have PD higher than 1%. There is an evidence that on average the isolated WDs of DC type have higher PD (with median PD of 0.78%) than the isolated DA type WDs (with median PD of 0.36%). On the other hand, the median PD of isolated DA type WDs is almost the same, i.e. 0.36% as the median PD of DA type white dwarfs in binary systems with red dwarfs (dM type), i.e. 0.33%. This shows, as expected, that there is no contribution to the PD from the companion if the WD companion is the red dwarf, which is the most common situation for WDs binary systems. We do not find differences in the polarization degree between magnetic and non-magnetic WDs. Because 97% of WDs in our sample have PD lower than 1%, they can be used as faint zero--polarized standard star in the magnitude range from 13 up to 17 of SDSS r filter. They cover the Northern sky between 13 hour to 23 hour in right ascension and from -11 degrees to 78 degrees in declination. Additionally, we found that for low extinction values (< 0.04) the best model that describes the dependence of PD on E(B-V) is given by the equation: PD_{max, ISM}[%] = 0.65 x E(B-V)^{0.12}.

preprint2016arXiv

Optical polarization map of the Polaris Flare with RoboPol

The stages before the formation of stars in molecular clouds are poorly understood. Insights can be gained by studying the properties of quiescent clouds, such as their magnetic field structure. The plane-of-the-sky orientation of the field can be traced by polarized starlight. We present the first extended, wide-field ($\sim$10 $\rm deg^2$) map of the Polaris Flare cloud in dust-absorption induced optical polarization of background stars, using the RoboPol polarimeter at the Skinakas Observatory. This is the first application of the wide-field imaging capabilities of RoboPol. The data were taken in the R-band and analysed with the automated reduction pipeline of the instrument. We present in detail optimizations in the reduction pipeline specific to wide-field observations. Our analysis resulted in reliable measurements of 641 stars with median fractional linear polarization 1.3%. The projected magnetic field shows a large scale ordered pattern. At high longitudes it appears to align with faint striations seen in the Herschel-SPIRE map of dust emission (250 $μm$), while in the central 4-5 deg$^2$ it shows an eddy-like feature. The overall polarization pattern we obtain is in good agreement with large scale measurements by Planck of the dust emission polarization in the same area of the sky.

preprint2016arXiv

RoboPol: Do optical polarization rotations occur in all blazars?

We present a new set of optical polarization plane rotations in blazars, observed during the third year of operation of RoboPol. The entire set of rotation events discovered during three years of observations is analysed with the aim of determining whether these events are inherent in all blazars. It is found that the frequency of the polarization plane rotations varies widely among blazars. This variation cannot be explained either by a difference in the relativistic boosting or by selection effects caused by a difference in the average fractional polarization. We conclude that the rotations are characteristic of a subset of blazars and that they occur as a consequence of their intrinsic properties.

preprint2016arXiv

RoboPol: First season rotations of optical polarization plane in blazars

We present first results on polarization swings in optical emission of blazars obtained by RoboPol, a monitoring program of an unbiased sample of gamma-ray bright blazars specially designed for effective detection of such events. A possible connection of polarization swing events with periods of high activity in gamma rays is investigated using the dataset obtained during the first season of operation. It was found that the brightest gamma-ray flares tend to be located closer in time to rotation events, which may be an indication of two separate mechanisms responsible for the rotations. Blazars with detected rotations have significantly larger amplitude and faster variations of polarization angle in optical than blazars without rotations. Our simulations show that the full set of observed rotations is not a likely outcome (probability $\le 1.5 \times 10^{-2}$) of a random walk of the polarization vector simulated by a multicell model. Furthermore, it is highly unlikely ($\sim 5 \times 10^{-5}$) that none of our rotations is physically connected with an increase in gamma-ray activity.

preprint2016arXiv

RoboPol: optical polarization-plane rotations and flaring activity in blazars

We present measurements of rotations of the optical polarization of blazars during the second year of operation of RoboPol, a monitoring programme of an unbiased sample of gamma-ray bright blazars specially designed for effective detection of such events, and we analyse the large set of rotation events discovered in two years of observation. We investigate patterns of variability in the polarization parameters and total flux density during the rotation events and compare them to the behaviour in a non-rotating state. We have searched for possible correlations between average parameters of the polarization-plane rotations and average parameters of polarization, with the following results: (1) there is no statistical association of the rotations with contemporaneous optical flares; (2) the average fractional polarization during the rotations tends to be lower than that in a non-rotating state; (3) the average fractional polarization during rotations is correlated with the rotation rate of the polarization plane in the jet rest frame; (4) it is likely that distributions of amplitudes and durations of the rotations have physical upper bounds, so arbitrarily long rotations are not realised in nature.

preprint2016arXiv

RoboPol: The optical polarization of gamma-ray--loud and gamma-ray--quiet blazars

We present average R-band optopolarimetric data, as well as variability parameters, from the first and second RoboPol observing season. We investigate whether gamma- ray--loud and gamma-ray--quiet blazars exhibit systematic differences in their optical polarization properties. We find that gamma-ray--loud blazars have a systematically higher polarization fraction (0.092) than gamma-ray--quiet blazars (0.031), with the hypothesis of the two samples being drawn from the same distribution of polarization fractions being rejected at the 3σ level. We have not found any evidence that this discrepancy is related to differences in the redshift distribution, rest-frame R-band lu- minosity density, or the source classification. The median polarization fraction versus synchrotron-peak-frequency plot shows an envelope implying that high synchrotron- peaked sources have a smaller range of median polarization fractions concentrated around lower values. Our gamma-ray--quiet sources show similar median polarization fractions although they are all low synchrotron-peaked. We also find that the random- ness of the polarization angle depends on the synchrotron peak frequency. For high synchrotron-peaked sources it tends to concentrate around preferred directions while for low synchrotron-peaked sources it is more variable and less likely to have a pre- ferred direction. We propose a scenario which mediates efficient particle acceleration in shocks and increases the helical B-field component immediately downstream of the shock.

preprint2014arXiv

Discovery of X-ray pulsations in the Be/X-ray binary IGR J21343+4738

We report on the discovery of X-ray pulsations in the Be/X-ray binary IGR J21343+4738 during an XMM-Newton observation. We obtained a barycentric corrected pulse period of 320.35+-0.06 seconds. The pulse profile displays a peak at low energy that flattens at high energy. The pulse fraction is 45+-3$% and independent of energy within the statistical uncertainties. The 0.2-12 keV spectrum is well fit by a two component model consisting of a blackbody with kT=0.11+-0.01 keV and a power law with photon index Gamma=1.02+-0.07. Both components are affected by photoelectric absorption with a equivalent hydrogen column density NH=(1.08+-0.15)x 10^{22} cm^{-2} The observed unabsorbed flux is 1.4x10^{-11} erg cm^{-2} s^{-1} in the 0.2-12 keV energy band. Despite the fact that the Be star's circumstellar disc has almost vanished, accretion continues to be the main source of high energy radiation. We argue that the observed X-ray luminosity (LX~10^{35} erg s^{-1}) may result from accretion via a low-velocity equatorial wind from the optical companion.

preprint2014arXiv

Early-time polarized optical light curve of GRB 131030A

We report the polarized optical light curve of a gamma-ray burst afterglow obtained using the RoboPol instrument. Observations began 655 seconds after the initial burst of gamma-rays from GRB131030A, and continued uninterrupted for 2 hours. The afterglow displayed a low, constant fractional linear polarization of $p = (2.1 \pm 1.6)\,\%$ throughout, which is similar to the interstellar polarization measured on nearby stars. The optical brightness decay is consistent with a forward-shock propagating in a medium of constant density, and the low polarization fraction indicates a disordered magnetic field in the shock front. This supports the idea that the magnetic field is amplified by plasma instabilities on the shock front. These plasma instabilities produce strong magnetic fields with random directions on scales much smaller than the total observable region of the shock, and the resulting randomly-oriented polarization vectors sum to produce a low net polarization over the total observable region of the shock.

preprint2014arXiv

Long-term variability of high-mass X-ray binaries. I.Photometry

We present photometric observations of the field around the optical counterparts of high-mass X-ray binaries. Our aim is to study the long-term photometric variability in correlation with their X-ray activity and derive a set of secondary standard stars that can be used for time series analysis. We find that the donors in Be/X-ray binaries exhibit larger amplitude changes in the magnitudes and colours than those hosting a supergiant companion. The amplitude of variability increases with wavelength in Be/X-ray binaries and remains fairly constant in supergiant systems. When time scales of years are considered, a good correlation between the X-ray and optical variability is observed. The X-rays cease when optical brightness decreases. These results reflect the fact that the circumstellar disk in Be/X-ray binaries is the main source of both optical and X-ray variability. We also derive the colour excess, E(B-V), selecting data at times when the contribution of the circumstellar disk was supposed to be at minimum, and we revisit the distance estimates.

preprint2014arXiv

The high optical polarization in the Be/X-ray binary EXO 2030+375

Polarization in classical Be stars results from Thomson scattering of the unpolarized light from the Be star in the circumstellar disc. Theory and observations agree that the maximum degree of polarization from isolated Be stars is < 4%. We report on the first optical polarimetric observations of the Be/X-ray binary EXO\,2030+375. We find that the optical (R band) light is strongly linearly polarized with a degee of polarization of 19%, the highest ever measured either in a classical or Be/X-ray binary. We argue that the interstellar medium cannot account for this high polarization degree and that a substantial amount must be intrinsic to the source. We propose that it may result from the alignment of non-spherical ferromagnetic grains in the Be star disc due to the strong neutron star magnetic field.

preprint2014arXiv

The Large Observatory For x-ray Timing

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

preprint2014arXiv

The quiescent state of the accreting X-ray pulsar SAX J2103.5+4545

We present an X-ray timing and spectral analysis of the Be/X-ray binary SAX J2103.5+4545 at a time when the Be star's circumstellar disk had disappeared and thus the main reservoir of material available for accretion had extinguished. In this very low optical state, pulsed X-ray emission was detected at a level of L_X~10^{33} erg/s. This is the lowest luminosity at which pulsations have ever been detected in an accreting pulsar. The derived spin period is 351.13 s, consistent with previous observations. The source continues its overall long-term spin-up, which reduced the spin period by 7.5 s since its discovery in 1997. The X-ray emission is consistent with a purely thermal spectrum, represented by a blackbody with kT=1 keV. We discuss possible scenarios to explain the observed quiescent luminosity and conclude that the most likely mechanism is direct emission resulting from the cooling of the polar caps, heated either during the most recent outburst or via intermittent accretion in quiescence.

preprint2014arXiv

The RoboPol optical polarization survey of gamma-ray - loud blazars

We present first results from RoboPol, a novel-design optical polarimeter operating at the Skinakas Observatory in Crete. The data, taken during the May - June 2013 commissioning of the instrument, constitute a single-epoch linear polarization survey of a sample of gamma-ray - loud blazars, defined according to unbiased and objective selection criteria, easily reproducible in simulations, as well as a comparison sample of, otherwise similar, gamma-ray - quiet blazars. As such, the results of this survey are appropriate for both phenomenological population studies and for tests of theoretical population models. We have measured polarization fractions as low as $0.015$ down to $R$ magnitude of 17 and as low as $0.035$ down to 18 magnitude. The hypothesis that the polarization fractions of gamma-ray - loud and gamma-ray - quiet blazars are drawn from the same distribution is rejected at the $10^{-3}$ level. We therefore conclude that gamma-ray - loud and gamma-ray - quiet sources have different optical polarization properties. This is the first time this statistical difference is demonstrated in optical wavelengths. The polarization fraction distributions of both samples are well-described by exponential distributions with averages of $\langle p \rangle =6.4 ^{+0.9}_{-0.8}\times 10^{-2}$ for gamma-ray--loud blazars, and $\langle p \rangle =3.2 ^{+2.0}_{-1.1}\times 10^{-2}$ for gamma-ray--quiet blazars. The most probable value for the difference of the means is $3.4^{+1.5}_{-2.0}\times 10^{-2}$. The distribution of polarization angles is statistically consistent with being uniform.

preprint2014arXiv

The RoboPol Pipeline and Control System

We describe the data reduction pipeline and control system for the RoboPol project. The RoboPol project is monitoring the optical $R$-band magnitude and linear polarization of a large sample of active galactic nuclei that is dominated by blazars. The pipeline calibrates and reduces each exposure frame, producing a measurement of the magnitude and linear polarization of every source in the $13'\times 13'$ field of view. The control system combines a dynamic scheduler, real-time data reduction, and telescope automation to allow high-efficiency unassisted observations.

preprint2013arXiv

Disc-loss episode in the Be shell optical counterpart to the high-mass X-ray binary IGR J21343+4738

The main goal of this work is to determine the properties of the optical counterpart to the INTEGRAL source IGR J21343+4738, and study its long-term optical variability. We present optical photometric BVRI and spectroscopic observations covering the wavelength band 4000-7500 A. We find that the optical counterpart to IGR J21343+4738 is a V=14.1 B1IVe shell star located at a distance of ~8.5 kpc. The Halpha line changed from an absorption dominated profile to an emission dominated profile, and then back again into absorption. In addition, fast V/R asymmetries were observed once the disc developed. Although the Balmer lines are the most strongly affected by shell absorption, we find that shell characteristics are also observed in He I lines. The optical spectral variability of IGR J21343+4738 is attributed to the formation of an equatorial disc around the Be star and the development of an enhanced density perturbation that revolves inside the disc. We have witnessed the formation and dissipation of the circumstellar disc. The strong shell profile of the Halpha and He I lines and the fact that no transition from shell phase to a pure emission phase is seen imply that we are seeing the system near edge-on.

preprint2013arXiv

Evidence for a change in the radiation mechanism in the hard state of GRO J1655-40. Hysteresis in the broad-band noise components

We have analysed RXTE archival data to study the aperiodic variability of the black-hole binary GRO J1655-40 during the hard state of the 2005 outburst. We computed the 0.008-64 Hz power spectral density during the rise and decay of the 2005 outburst, and we found that they were reasonably well modelled by the sum of two, broad Lorenztian functions in most cases (plus a narrow QPO) which correspond to three different variability components. Our aim is to study the evolution of the timing properties of the source during the outburst, by studying the correlation between the characteristics of the broad-band noise components in the power spectra and the source luminosity. Our results suggest that the whole power spectrum shifts to high (low) frequencies as the source luminosity increases (decreases), in agreement with previous studies of other black hole binaries. However, we also detect a strong "hysteresis" pattern in the "frequency-luminosity" plots, and show that the "critical" luminosity limit, above which the timing properties of the source change, is different during the rise and the decay phase of the outburst. We discuss the general implications of these results in the context of the truncated disc model.

preprint2013arXiv

Spectral states in Be/X-ray pulsars

In the last quarter of a century, a unified characterization of the spectral evolution of low-mass X-ray binaries, both containing a neutron star and a black hole, was possible. In this context, the notion of source states characterizing the X-ray emission from black-hole binaries and neutron-star low-mass X-ray binaries revealed to be a very useful tool to disentangle the complex spectral and aperiodic phenomenology displayed by those classes of accreting objects. Be/X-ray binaries constitute another major class of transient accreting binaries, for which very little work has been done on the correlated timing and spectral variability. Especially, no definition of source states exists for this class, in spite of their highly variable X-ray emission. When active, Be/X-ray binaries are among the brightest objects in the X-ray sky and are characterized by dramatic variability in brightness on timescales ranging from seconds to years. It is then worth it to ask whether a definition of spectral states is possible for these systems. In this work, we try to address such a question, investigating whether accreting X-ray pulsars display source states and characterizing those states through their spectral properties. Our results show that Be/X-ray pulsars trace two different branches in their hardness-intensity diagram: the horizontal branch, a low-intensity state, and the diagonal branch, a high-intensity state that only appears when the X-ray luminosity exceeds a critical limit. We propose that the two branches are the phenomenological signature of two different accretion modes -- in agreement with recently proposed models -- depending on whether the luminosity of the source is above or below a critical value.

preprint2012arXiv

Accreting magnetars: a new type of high-mass X-ray binaries?

The discovery of very slow pulsations (Pspin=5560s) has solved the long-standing question of the nature of the compact object in the high-mass X-ray binary 4U 2206+54 but has posed new ones. According to spin evolutionary models in close binary systems, such slow pulsations require a neutron star magnetic field strength larger that the quantum critical value, suggesting the presence of a magnetar. We present the first XMM-Newton observations of 4U 2206+54 and investigate its spin evolution. We find that the observed spin-down rate agrees with the magnetar scenario. We analyse ISGRI/INTEGRAL observations of 4U 2206+54 to search for the previously suggested cyclotron resonance scattering feature at ~30 keV. We do not find a clear indication of the presence of the line, although certain spectra display shallow dips, not always at 30 keV. The association of these dips with a cyclotron line is very dubious because of its apparent transient nature. We also investigate the energy spectrum of 4U 2206+54 in the energy range 0.3-10 keV with unprecedented detail and report for the first time the detection of very weak 6.5 keV fluorescence iron lines. The photoelectric absorption is consistent with the interstellar value, indicating very small amount of local matter, which would explain the weakness of the florescence lines. The lack of matter locally to the source may be the consequence of the relatively large orbital separation of the two components of the binary. The wind would be too tenuous in the vicinity of the neutron star.

preprint2012arXiv

LOFT: the Large Observatory For X-ray Timing

The LOFT mission concept is one of four candidates selected by ESA for the M3 launch opportunity as Medium Size missions of the Cosmic Vision programme. The launch window is currently planned for between 2022 and 2024. LOFT is designed to exploit the diagnostics of rapid X-ray flux and spectral variability that directly probe the motion of matter down to distances very close to black holes and neutron stars, as well as the physical state of ultra-dense matter. These primary science goals will be addressed by a payload composed of a Large Area Detector (LAD) and a Wide Field Monitor (WFM). The LAD is a collimated (<1 degree field of view) experiment operating in the energy range 2-50 keV, with a 10 m^2 peak effective area and an energy resolution of 260 eV at 6 keV. The WFM will operate in the same energy range as the LAD, enabling simultaneous monitoring of a few-steradian wide field of view, with an angular resolution of <5 arcmin. The LAD and WFM experiments will allow us to investigate variability from submillisecond QPO's to year-long transient outbursts. In this paper we report the current status of the project.

preprint2012arXiv

New insights into the Be/X-ray binary system MXB 0656-072

The X-ray transient MXB 0656-072 is a poorly studied member of high-mass X-ray binaries. Based on the transient nature of the X-ray emission, the detection of pulsations, and the early-type companion, it has been classified as a Be X-ray binary (Be/XRB). However, the flaring activity covering a large fraction of a giant outburst is somehow peculiar. Our goal is to investigate the multiwavelength variability of the high-mass X-ray binary MXB 0656-072. We carried out optical spectroscopy and analysed all RXTE archive data, performing a detailed X-ray-colour, spectral, and timing analysis of both normal (type-I) and giant (type-II) outbursts from MXB 0656-072 This is the first detailed analysis of the optical counterpart in the classification region (4000-5000 A). From the strength and ratio of the elements and ions, we derive an O9.5Ve spectral type, in agreement with previous classification. This confirms its Be nature. The characterisation of the Be/XRB system relies on Balmer lines in emission in the optical spectra, long-term X-ray variability, and the orbital period vs. spin period and EW(Hα) relation. The peculiar feature that distinguishes the type-II outburst is flaring activity, which occurs during the whole outburst peak, before a smoother decay. We interpret it in terms of magneto-hydrodynamic instability. Colour and spectral analysis reveal a hardening of the spectrum as the flux increases. We explored the aperiodic X-ray variability of the system for the first time, finding a correlation of the central frequency and rms of the main timing component with luminosity, which extends up to a "saturation" flux of 1E-8 erg/cm^2/s. A correlation between timing and spectral parameters was also found, pointing to an interconnection between the two physical regions responsible for both phenomenologies.

preprint2011arXiv

Multi-frequency observations of SWIFT J1626.6-5156

SWIFT J1626.6-5156 is an X-ray pulsar that was discovered in December 2005 during an X-ray outburst. Although the X-ray data suggest that the system is a high-mass X-ray binary, very little information exists on the nature of the optical counterpart. We have analysed all RXTE observations since its discovery, archived optical spectroscopic and photometric data and obtained for the first time near-IR spectra. The K-band spectrum shows HeI 20581 A and HI 21660 A (Brackett-gamma) in emission, which confine the spectral type of the companion to be earlier than B2.5. The H-band spectrum exhibits the HI Br-18-11 recombination series in emission. The most prominent feature of the optical band spectrum is the strong emission of the Balmer line Halpha. The 4000-5000 A spectrum contains HeII and numerous HeI ines in absorption, indicating an early B-type star. The source shows three consecutive stages characterised by different types of variability in the X-ray band: a smooth decay after the peak of a large outburst, large-amplitude flaring variability (reminiscent of type I oytbursts) and quiescence. We observed that the spectrum becomes softer as the flux decreases and that this is a common characteristic of the X-ray emission for all observing epochs. An emission line feature at ~6.5 keV is also always present. The X-ray/optical/IR continuum and spectral features are typical of an accreting X-ray pulsar with an early-type donor. The long-term X-ray variability is also characteristic of hard X-ray transients. We conclude that SWIFT J1626.6-5156 is a Be/X-ray binary with a B0Ve companion located at a distance of ~10 kpc.

preprint2010arXiv

Non-radial pulsations in the Be/X binaries 4U0115+63 and SAXJ2103.5+4545

The discovery of non-radial pulsations (NRP) in the Be/X binaries of the Magellanic Clouds (MC, eg. Fabrycky 2005, Coe et al. 2005, Schmidtke & Cowley 2005) provided a new approach to understand these complex systems, and, at the same time, favoured the synergy between two different fields: stellar pulsations and X-ray binaries. This breakthrough was possible thanks to the MACHO and OGLE surveys. However, in our Galaxy, only two Be/X have been reported to show NRP: GROJ2058+42 (Kiziloglu et al. 2007) and LSI+61 235 (Sarty et al. 2009). Our objective is to study the short-term variability of Galactic Be/X binaries, compare them to the Be/X of the MC and to the isolated Galactic Be observed with CoRoT and Kepler. We present preliminary results of two Be/X stars, namely 4U0115+63 and SAXJ2103.5+4545 showing multiperiodicity and periodicity respectively, most probably produced by non-radial pulsations.

preprint2010arXiv

The optical counterpart to IGR J06074+2205: a Be/X-ray binary showing disc loss and V/R variability

Present X-ray missions are regularly discovering new X/gamma-ray sources. The identification of the counterparts of these high-energy sources at other wavelengths is important to determine their nature. In particular, optical observations are an essential tool in the study of X-ray binary populations in our Galaxy. The main goal of this work is to determine the properties of the optical counterpart to the INTEGRAL source IGR J06074+2205, and study its long-term optical variability. Although its nature as a high-mass X-ray binary has been suggested, little is known about its physical parameters. We have been monitoring IGR J06074+2205 since 2006 in the optical band. We present optical photometric BVRI and spectroscopic observations covering the wavelength band 4000-7000 A. The blue spectra allow us to determine the spectral type and luminosity class of the optical companion; the red spectra, together with the photometric magnitudes, were used to derive the colour excess E(B-V) and estimate the distance. We have carried out the first detailed optical study of the massive component in the high-mass X-ray binary IGR J06074+2205. We find that the optical counterpart is a V=12.3 B0.5Ve star located at a distance of ~4.5 kpc. The monitoring of the Halpha line reveals V/R variability and an overall decline of its equivalent width. The Halpha line has been seen to revert from an emission to an absorption profile. We attribute this variability to global changes in the structure of the Be star's circumstellar disc which eventually led to the complete loss of the disc. The density perturbation that gives rise to the V/R variability vanishes when the disc becomes too small.

preprint2010arXiv

The scaling of X-ray variability with luminosity in Ultra-luminous X-ray sources

We investigated the relationship between the X-ray variability amplitude and X-ray luminosity for a sample of 14 bright Ultra-luminous X-ray sources (ULXs) with XMM-Newton/EPIC data, and compare it with the well established similar relationship for Active Galactic Nuclei (AGN). We computed the normalised excess variance in the 2-10 keV light curves of these objects and their 2-10 keV band intrinsic luminosity. We also determined model "variability-luminosity" relationships for AGN, under several assumptions regarding their power-spectral shape. We compared these model predictions at low luminosities with the ULX data. The variability amplitude of the ULXs is significantly smaller than that expected from a simple extrapolation of the AGN "variability-luminosity" relationship at low luminosities. We also find evidence for an anti-correlation between the variability amplitude and L(2-10 keV) for ULXs. The shape of this relationship is consistent with the AGN data but only if the ULXs data are shifted by four orders of magnitudes in luminosity. Most (but not all) of the ULXs could be "scaled-down" version of AGN if we assume that: i) their black hole mass and accretion rate are of the order of ~(2.5-30)x 10E+03 Msolar and ~ 1-80 % of the Eddington limit, and ii) their Power Spectral Density has a doubly broken power-law shape. This PDS shape and accretion rate is consistent with Galactic black hole systems operating in their so-called "low-hard" and "very-high" states.

preprint2009arXiv

Characterizing black hole variability with nonlinear methods: the case of the X-ray Nova 4U 1543-47

We investigate the possible nonlinear variability properties of the black hole X-ray nova 4U1543-47 to complement the temporal studies based on linear techniques, and to search for signs of nonlinearity in Galactic black hole (GBH) light curves. First, we apply the weighted scaling index method (WSIM) to characterize the X-ray variability properties of 4U1543-47 in different spectral states during the 2002 outburst. Second, we use surrogate data to investigate whether the variability is nonlinear in any of the different spectral states. The main findings can be summarized as follows. The mean weighted scaling index appears to be able to parametrize uniquely the temporal variability properties of this GBH: the 3 different spectral states of the 2002 outburst of 4U1543-47 are characterized by different and well constrained values. The search for nonlinearity reveals that the variability is linear in all light curves with the notable exception of the very high state. Our results imply that we can use the WSIM to assign a single number, namely the mean weighted scaling index, to a light curve, and in this way discriminate among the different spectral states of a source. The detection of nonlinearity in the VHS, that is characterized by the presence of most prominent QPOs, suggests that intrinsically linear models which have been proposed to account for the low frequency QPOs in GBHs may be ruled out (abridged).

preprint2009arXiv

Correlated optical/X-ray variability in the high-mass X-ray binary SAX J2103.5+4545

SAX J2103.5+4545 is the Be/X-ray binary with the shortest orbital period. It shows extended bright and faint X-ray states that last for a few hundred days. The main objective of this work is to investigate the relationship between the X-ray and optical variability and to characterise the spectral and timing properties of the bright and faint states. We have found a correlation between the spectral and temporal parameters that fit the energy and power spectra. Softer energy spectra correspond to softer power spectra. That is to say, when the energy spectrum is soft the power at high frequencies is suppressed. We also present the results of our monitoring of the Halpha line of the optical counterpart since its discovery in 2003. There is a correlation between the strength and shape of the Halpha line, originated in the circumstellar envelope of the massive companion and the X-ray emission from the vicinity of the neutron star. Halpha emission, indicative of an equatorial disc around the B-type star, is detected whenever the source is bright in X-rays. When the disc is absent, the X-ray emission decreases significantly. The long-term variability of SAX J2103.5+4545 is characterised by fast episodes of disc loss and subsequent reformation. The time scales for the loss and reformation of the disc (about 2 years) are the fastest among Be/X-ray binaries.

preprint2009arXiv

Spectroscopy of the bright optical counterparts of X-ray sources in the direction of M 31. II

A recent survey of the Local Group spiral galaxy M 31 with XMM-Newton yielded a large number of X-ray sources. This is the second in a series of papers with the aim of identifying the optical counterparts of these X-ray sources. We have obtained optical spectra for 21 bright optical counterparts of 20 X-ray sources in the direction of M 31, using the 1.3-m Skinakas telescope in Crete, Greece. For 17 of the 20 X-ray sources, we have identified the optical counterpart as a normal late type star (of type F or later) in the foreground (i.e. in the Milky Way). For two more sources there were two possible optical counterparts in each case, while two more objects have X-ray properties that are not compatible with the spectral characteristics of late type non-flaring stars.

preprint2006arXiv

Fourier resolved spectroscopy of 4U 1543-47 during the 2002 outburst

We have obtained Fourier-resolved spectra of the black-hole binary 4U 1543-47 in the canonical states (high/soft, very high, intermediate and low/hard) observed in this source during the decay of an outburst that took place in 2002. Our objective is to investigate the variability of the spectral components generally used to describe the energy spectra of black-hole systems, namely a disk component, a power-law component attributed to Comptonization by a hot corona and the contribution of the iron line due to reprocessing of the high energy (E > 7 keV) radiation. We find that i) the disk component is not variable on time scales shorter than ~100 seconds, ii) the reprocessing emission as manifest by the variability of the Fe Kalpha line responds to the primary radiation variations down to time scales of ~70 ms in the high and very-high states, but longer than 2 s in the low state, iii) the low-frequency QPOs are associated with variations of the X-ray power law spectral component and not to the disk component and iv) the spectra corresponding to the highest Fourier frequency are the hardest (show the flatter spectra) at a given spectral state. These results question models that explain the observed power spectra as due to modulations of the accretion rate alone, as such models do not provide any apparent reason for a Fourier frequency dependence of the power law spectral indices.