A Be-type star with a black-hole companion

Stellar-mass black holes have all been discovered through X-ray emission, which arises from the accretion of gas from their binary companions (this gas is either stripped from low-mass stars or supplied as winds from massive ones). Binary evolution models also predict the existence of black holes accreting from the equatorial envelope of rapidly spinning Be-type stars (stars of the Be type are hot blue irregular variables showing characteristic spectral emission lines of hydrogen). Of the ~80 Be X-ray binaries known in the Galaxy, however, only pulsating neutron stars have been found as companions. A black hole was formally allowed as a solution for the companion to the Be star MWC 656 (also known as HD 215227), although that was based on a single radial velocity curve of the Be star, a mistaken spectral classification and rough estimates of the inclination angle. Here we report observations of an accretion disk line mirroring the orbit of the Be star. This, together with an improved radial velocity curve of the Be star through fitting sharp Fe II profiles from the equatorial disk, and a refined Be classification (to that of a B1.5-B2 III star), reveals a black hole of 3.8 to 6.9 solar masses orbiting MWC 656, the candidate counterpart of the gamma-ray source AGL J2241+4454. The black hole is X-ray quiescent and fed by a radiatively inefficient accretion flow giving a luminosity less than 1.6 x 10-7 times the Eddington luminosity. This implies that Be binaries with black-hole companions are difficult to detect by conventional X-ray surveys.

On the binary nature of the γ-ray sources AGL J2241+4454 (= MWC 656) and HESS J0632+057 (= MWC 148)

We present optical spectroscopy of MWC 656 and MWC 148, the proposed optical counterparts of the gamma-ray sources AGL J2241+4454 and HESS J0632+0 57, respectively. The main parameters of the Halpha emission line (EW, FWHM and centroid velocity) in these stars are modulated on the proposed orbital periods of 60.37 and 321 days, respectively. These modulations are likely produced by the resonant interaction of the Be discs with compact stars in eccentric orbits. We also present radial velocity curves of the optical stars folded on the above periods and obtain the first orbital elements of the two gamma-ray sources thus confirming their binary nature. Our orbital solution support eccentricities e~0.4 and 0.83+-0.08 for MWC 656 and MWC 148, respectively. Further, our orbital elements imply that the X-ray outbursts in HESS J0632+057/MWC 148 are delayed ~0.3 orbital phases after periastron passage, similarly to the case of LS I +61 303. In addition, the optical photometric light curve maxima in AGL J2241+4454/MWC 656 occur ~0.25 phases passed periastron, similar to what is seen in LS I +61 303. We also find that the orbital eccentricity is correlated with orbital period for the known gamma-ray binaries. This is explained by the fact that small stellar separations are required for the efficient triggering of VHE radiation. Another correlation between the EW of Halpha and orbital period is also observed, similarly to the case of Be/X-ray binaries. These correlations are useful to provide estimates of the key orbital parameters Porb and e from the Halpha line in future Be gamma-ray binary candidates.

Gamma-ray emission from massive young stellar objects

Context.Massive stars form in dense and massive molecular cores. The exact formation mechanism is unclear, but it is possible that some massive stars are formed by processes similar to those that produce the low-mass stars, with accretion/ejection phenomena occurring at some point of the evolution of the protostar. This picture seems to be supported by the detection of a collimated stellar wind emanating from the massive protostar IRAS 16547-4247. A triple radio source is associated with the protostar: a compact core and two radio lobes. The emission of the southern lobe is clearly non-thermal. Such emission is interpreted as synchrotron radiation produced by relativistic electrons locally accelerated at the termination point of a thermal jet. Since the ambient medium is determined by the properties of the molecular cloud in which the whole system is embedded, we can expect high densities of particles and infrared photons. Because of the confirmed presence of relativistic electrons, inverse Compton and relativistic Bremsstrahlung interactions are unavoidable. Aims.We aim to make quantitative predictions of the spectral energy distribution of the non-thermal spots generated by massive young stellar objects, with emphasis on the particular case of IRAS 16547-4247. Methods.We study the high-energy emission generated by the relativistic electrons which produce the non-thermal radio source in IRAS 16547-4247. We also study the result of proton acceleration at the terminal shock of the thermal jet and make estimates of the secondary gamma rays and electron-positron pairs produced by pion decay. Results.We present spectral energy distributions for the southern lobe of IRAS 16547-4247, for a variety of conditions. We show that high-energy emission might be detectable from this object in the gamma-ray domain. The source may also be detectable in X-rays through long exposures with current X-ray instruments. Conclusions.Gamma-ray telescopes such as GLAST, and even ground-based Cherenkov arrays of new generation can be used to study non-thermal processes occurring during the formation of massive stars.

Photometric and Halpha observations of LSI+61·303: detection of a 26 day V and JHK band modulation

We present new optical and infrared photometric observations and high resolution H α spectra of the periodic radio star LSI+61◦303. The optical photometric data set covers the time interval 1985-1993 and amounts to about a hundred nights. A period of ∼26 days is found in the V band. The infrared data also present evidence for a similar periodicity, but with higher amplitude of variation ((0.m 2). The spectroscopic observations include 16 intermediate and high dispersion spectra of LSI+61◦303 collected between January 1989 and February 1993. The H α emission line profile and its variations are analyzed. Several emission line parameters -- among them the H α EW and the width of the H α red hump -- change strongly at or close to radio maximum, and may exhibit periodic variability. We also observe a significant change in the peak separation. The H α profile of LSI+61◦303 does not seem peculiar for a Be star. However, several of the observed variations of the H α profile can probably be associated with the presence of the compact, secondary star.

Evidence of Halpha periodicities in LSI+61·303

We present the results of analyzing H$\alpha$ spectra of the radio emitting X-ray binary LS I+61303. For the first time, the same 26.5 d radio period is clearly detected in the H$\alpha$ emission line. Moreover, the equivalent width and the peak separation of the H$\alpha$ emission line seem also to vary over a time scale of 1600 days. This points towards the $\sim4$ yr modulation, detected in the radio outburst amplitude, being probably a result of variations in the mass loss rate of the Be star and/or density variability in the circumstellar disk. In addition, the dependence of the peak separation from the equivalent width informs us that the LS I+61303 circumstellar disk is among the densest of Be-stars.

Leptonic secondary emission in a hadronic microquasar model

Context.It has been proposed that the origin of the very high-energy photons emitted from high-mass X-ray binaries with jet-like features, so-called microquasars (MQs), is related to hadronic interactions between relativistic protons in the jet and cold protons of the stellar wind. Leptonic secondary emission should be calculated in a complete hadronic model that includes the effects of pairs from charged pion decays inside the jets and the emission from pairs generated by gamma-ray absorption in the photosphere of the system. Aims.We aim at predicting the broadband spectrum from a general hadronic microquasar model, taking into account the emission from secondaries created by charged pion decay inside the jet. Methods.The particle energy distribution for secondary leptons injected along the jets is consistently derived taking the energy losses into account. The spectral energy distribution resulting from these leptons is calculated after assuming different values of the magnetic field inside the jets. We also compute the spectrum of the gamma-rays produced by neutral pion-decay and processed by electromagnetic cascades under the stellar photon field. Results.We show that the secondary emission can dominate the spectral energy distribution at low energies (~1 MeV). At high energies, the production spectrum can be significantly distorted by the effect of electromagnetic cascades. These effects are phase-dependent, and some variability modulated by the orbital period is predicted.

The star forming region Monoceros R2 as a gamma-ray source

After the release of the gamma-ray source catalog produced by the Fermi satellite during its first two years of operation, a significant fraction of sources still remain unassociated at lower energies. In addition to well-known high-energy emitters (pulsars, blazars, supernova remnants, etc.), theoretical expectations predict new classes of gamma-ray sources. In particular, gamma-ray emission could be associated with some of the early phases of stellar evolution, but this interesting possibility is still poorly understood. Aims: The aim of this paper is to assess the possibility of the Fermi gamma-ray source 2FGL J0607.5-0618c being associated with the massive star forming region Monoceros R2. Methods: A multi-wavelength analysis of the Monoceros R2 region is carried out using archival data at radio, infrared, X-ray, and gamma-ray wavelengths. The resulting observational properties are used to estimate the physical parameters needed to test the different physical scenarios. Results: We confirm the 2FGL J0607.5-0618c detection with improved confidence over the Fermi two-year catalog. We find that a combined effect of the multiple young stellar objects in Monoceros R2 is a viable picture for the nature of the source.

A microquasar model applied to unidentified gamma-ray sources

Among unidentified gamma-ray sources in the galactic plane, there are some that present significant variability and have been proposed to be high-mass microquasars. To deepen the study of the possible association between variable low galactic latitude gamma-ray sources and microquasars, we have applied a leptonic jet model based on the microquasar scenario that reproduces the gamma-ray spectrum of three unidentified gamma-ray sources, 3EG J1735-1500, 3EG J1828+0142 and GRO J1411-64, and is consistent with the observational constraints at lower energies. We conclude that if these sources were generated by microquasars, the particle acceleration processes could not be as efficient as in other objects of this type that present harder gamma-ray spectra. Moreover, the dominant mechanism of high-energy emission should be synchrotron self-Compton (SSC) scattering, and the radio jets may only be observed at low frequencies. For each particular case, further predictions of jet physical conditions and variability generation mechanisms have been made in the context of the model. Although there might be other candidates able to explain the emission coming from these sources, microquasars cannot be excluded as counterparts. Observations performed by the next generation of gamma-ray instruments, like GLAST, are required to test the proposed model.

Radio continuum and near-infrared study of the MGRO J2019+37 region

MGRO J2019+37 is an unidentified extended source of very high energy gamma-rays originally reported by the Milagro Collaboration as the brightest TeV source in the Cygnus region. Its extended emission could be powered by either a single or several sources. The GeV pulsar AGL J2020.5+3653 , discovered by AGILE and associated with PSR J2021+3651 , could contribute to the emission from MGRO J2019+37 . Aims. Our aim is to identify radio and near-infrared sources in the field of the extended TeV source MGRO J2019+37 , and study potential counterparts to explain its emission. Methods. We surveyed a region of about 6 square degrees with the Giant Metrewave Radio Telescope (GMRT) at the frequency 610 MHz. We also observed the central square degree of this survey in the near-infrared -band using the 3.5 m telescope in Calar Alto. Archival X-ray observations of some specific fields are included. VLBI observations of an interesting radio source were performed. We explored possible scenarios to produce the multi-TeV emission from MGRO J2019+37 and studied which of the sources could be the main particle accelerator. Results. We present a catalogue of 362 radio sources detected with the GMRT in the field of MGRO J2019+37 , and the results of a cross-correlation of this catalog with one obtained at near-infrared wavelengths, which contains ~3105 sources, as well as with available X-ray observations of the region. Some peculiar sources inside the ~1° uncertainty region of the TeV emission from MGRO J2019+37 are discussed in detail, including the pulsar PSR J2021+3651 and its pulsar wind nebula PWN G75.2+0.1 , two new radio-jet sources, the H II region Sh 2-104 containing two star clusters, and the radio source NVSS J202032+363158 . We also find that the hadronic scenario is the most likely in case of a single accelerator, and discuss the possible contribution from the sources mentioned above. Conclusions. Although the radio and GeV pulsar PSR J2021+3651 / AGL J2020.5+3653 and its associated pulsar wind nebula PWN G75.2+0.1 can contribute to the emission from MGRO J2019+37 , extrapolation of the GeV spectrum does not explain the detected multi-TeV flux. Other sources discussed here could contribute to the emission of the Milagro source.

Radio continuum and near-ingrared study of the MGRO J2019+37 region

Context. MGRO J2019+37 is an unidentified extended source of very high energy gamma-rays originally reported by the Milagro Collaboration as the brightest TeV source in the Cygnus region. Its extended emission could be powered by either a single or several sources. The GeV pulsar AGL J2020.5+3653, discovered by AGILE and associated with PSR J2021+3651, could contribute to the emission from MGRO J2019+37. Our aim is to identify radio and near-infrared sources in the field of the extended TeV source MGRO J2019+37, and study potential counterparts to explain its emission. Methods: We surveyed a region of about 6 square degrees with the Giant Metrewave Radio Telescope (GMRT) at the frequency 610 MHz. We also observed the central square degree of this survey in the near-infrared Ks-band using the 3.5 m telescope in Calar Alto. Archival X-ray observations of some specific fields are included. VLBI observations of an interesting radio source were performed. We explored possible scenarios to produce the multi-TeV emission from MGRO J2019+37 and studied which of the sources could be the main particle accelerator. Results: We present a catalogue of 362 radio sources detected with the GMRT in the field of MGRO J2019+37, and the results of a cross-correlation of this catalog with one obtained at near-infrared wavelengths, which contains ∼3 × 105 sources, as well as with available X-ray observations of the region. Some peculiar sources inside the ∼1◦ uncertainty region of the TeV emission from MGRO J2019+37 are discussed in detail, including the pulsar PSR J2021+3651 and its pulsar wind nebula PWN G75.2+0.1, two new radio-jet sources, the Hii region Sh 2-104 containing two star clusters, and the radio source NVSS J202032+363158. We also find that the hadronic scenario is the most likely in case of a single accelerator, and discuss the possible contribution from the sources mentioned above. Conclusions: Although the radio and GeV pulsar PSR J2021+3651 / AGL J2020.5+3653 and its associated pulsar wind nebula PWN G75.2+0.1 can contribute to the emission from MGRO J2019+37, extrapolation of the GeV spectrum does not explain the detected multi-TeV flux. Other sources discussed here could contribute to the emission of the Milagro source

The changing milliarcsecond radio morphology of the gamma-ray binary LS 5039

LS 5039 is one of the few TeV emitting X-ray binaries detected so far. The powering source of its multiwavelength emission can be accretion in a microquasar scenario or wind interaction in a young nonaccreting pulsar scenario. Aims.To present new high-resolution radio images and compare them with the expected behavior in the different scenarios. Methods.We analyze Very Long Baseline Array (VLBA) radio observations that provide morphological and astrometric information at milliarcsecond scales. Results.We detect a changing morphology between two images obtained five days apart. In both runs there is a core component with a constant flux density, and an elongated emission with a position angle (PA) that changes by 12 $\pm$ $3\degr$ between both runs. The source is nearly symmetric in the first run and asymmetric in the second one. The astrometric results are not conclusive. Conclusions.A simple and shockless microquasar scenario cannot easily explain the observed changes in morphology. An interpretation within the young nonaccreting pulsar scenario requires the inclination of the binary system to be very close to the upper limit imposed by the absence of X-ray eclipses.

Classical line shapes based on analytical solutions of bimolecular trajectories in collision induced emission. II. Reactive collisions

The classical theory of collision induced emission (CIE) from pairs of dissimilar rare gas atoms was developed in Paper I [D. Reguera and G. Birnbaum, J. Chem. Phys. 125, 184304 (2006)] from a knowledge of the straight line collision trajectory and the assumption that the magnitude of the dipole could be represented by an exponential function of the inter-nuclear distance. This theory is extended here to deal with other functional forms of the induced dipole as revealed by ab initio calculations. Accurate analytical expression for the CIE can be obtained by least square fitting of the ab initio values of the dipole as a function of inter-atomic separation using a sum of exponentials and then proceeding as in Paper I. However, we also show how the multi-exponential fit can be replaced by a simpler fit using only two analytic functions. Our analysis is applied to the polar molecules HF and HBr. Unlike the rare gas atoms considered previously, these atomic pairs form stable bound diatomic molecules. We show that, interestingly, the spectra of these reactive molecules are characterized by the presence of multiple peaks. We also discuss the CIE arising from half collisions in excited electronic states, which in principle could be probed in photo-dissociation experiments.

High-energy gamma-ray emission from microquasars: LS 5039 and LS I+61303

The possible associations between the microquasars LS 5039 and LS I +61 303 and the EGRET sources 3EG J1824-1514 and 3EG J0241+6103 suggest that microquasars could also be sources of high-energy gamma-rays. In this work, we present a detailed numerical inverse Compton (IC) model, based on a microquasar scenario, that reproduces the high-energy gamma-ray spectra and variability observed by EGRET for the mentioned sources. Our model considers a population of relativistic electrons entrained in a cylindrical inhomogeneous jet that interact through IC scattering with both the radiation and the magnetic fields.

Exploring the high frequency emission of radio loud X-ray binaries

We report millimetre-wave continuum observations of the X-ray binaries Cygnus X-3, SS 433, LSI+61 303, Cygnus X-1 and GRS 1915+105. The observations were carried out with the IRAM 30 m-antenna at 250 GHz (1.25 mm) from 1998 March 14 to March 20. These millimetre measurements are complemented with centimetre observations from the Ryle Telescope, at 15 GHz (2.0 cm) and from the Green Bank Interferometer at 2.25 and 8.3 GHz (13 and 3.6 cm). Both Cygnus X-3 and SS 433 underwent moderate flaring events during our observations, whose main spectral evolution properties are described and interpreted. A significant spectral steepening was observed in both sources during the flare decay, that is likely to be caused by adiabatic expansion, inverse Compton and synchrotron losses. Finally, we also report 250 GHz upper limits for three additional undetected X-ray binary stars: LSI+65 010, LSI+61 235 and X Per.

Exploring the connection between the stellar wind and the non-thermal emission in LS 5039

Context.LS 5039 has been observed with several X-ray instruments so far showing quite steady emission in the long term and no signatures of accretion disk. The source also presents X-ray variability at orbital timescales in flux and photon index. The system harbors an O-type main sequence star with moderate mass-loss. At present, the link between the X-rays and the stellar wind is unclear. Aims.We study the X-ray fluxes, spectra, and absorption properties of LS 5039 at apastron and periastron passages during an epoch of enhanced stellar mass-loss, and the long term evolution of the latter in connection with the X-ray fluxes. Methods.New XMM-Newton observations were performed around periastron and apastron passages in September 2005, when the stellar wind activity was apparently higher. April 2005 Chandra observations on LS 5039 were revisited. Moreover, a compilation of H EW data obtained since 1992, from which the stellar mass-loss evolution can be approximately inferred, was carried out. Results.XMM-Newton observations show higher and harder emission around apastron than around periastron. No signatures of thermal emission or a reflection iron line indicating the presence of an accretion disk are found in the spectrum, and the hydrogen column density () is compatible with being the same in both observations and consistent with the interstellar value. 2005 Chandra observations show a hard X-ray spectrum, and possibly high fluxes, although pileup effects preclude conclusive results from being obtained. The H EW shows yearly variations of 10%, and does not seem to be correlated with X-ray fluxes obtained at similar phases, unlike what is expected in the wind accretion scenario. Conclusions.2005 XMM-Newton and Chandra observations are consistent with 2003 RXTE/PCA results, namely moderate flux and spectral variability at different orbital phases. The constancy of the seems to imply that either the X-ray emitter is located at 1012 cm from the compact object, or the density in the system is 3 to 27 times smaller than that predicted by a spherical symmetric wind model. We suggest that the multiwavelength non-thermal emission of LS 5039 is related to the observed extended radio jets and is unlikely to be produced inside the binary system.