EDP Sciences logo
Subscriber Authentication Point
Search
Menu
Home All issues Volume 621 (January 2019) A&A, 621 (2019) A37 Full HTML
Free Access
Issue
A&A
Volume 621, January 2019
Article Number A37
Number of page(s) 12
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/201832682
Published online 03 January 2019

© ESO 2019

1 Introduction

High-mass X-ray binaries (HMXBs) are systems consisting of a massive, early-type star (spectral class O or B) in the process of stellar evolution and a compact object, neutron star, or black hole, revolving around their common center of mass. The compact object accretes mass from the stellar wind of its companion or, if the companion star fills its Roche lobe, the matter flows directly onto the compact object (e.g., Lewin et al. 1997; Lewin & van der Klis 2006; Liu et al. 2006).

Owing to the past supernova explosion of one of its components, possibly combined with close dynamical encounters of the system with other stars in its parent cluster as well, a significant portion of HMXBs possess high (runaway) space velocities (e.g., van Oijen 1989; Kaper et al. 2004). Furthermore, their galactic height distribution also implies their runaway nature (van Oijen 1989).

Traditionally, runaway stars can be revealed via measurements of their proper motions, yielding their tangential velocities. However, this method requires the distance to the object to be determined and the proper motion to be measured with sufficient precision. Determination of the radial velocities via spectroscopic measurements is also a viable method. Nevertheless, these methods are only applicable to a portion of relatively close HMXBs, having sufficiently large proper motion to be reliably measured and/or bright optical counterpart for the spectroscopic studies. With the advent of modern high-energy missions, there has been a sharp increase in the number of confirmed and candidate HMXB systems that are detected at considerable distances, are often highly absorbed, and only have weak infrared counterparts. As they are hard objects to study with the above methods, the kinematics of these objects is poorly known.

Another means of searching for runaway star candidates is through the detection of bow shocks associated with these objects (e.g., van Buren & McCray 1988; van Buren & Noriega-Crespo 1995; Gvaramadze & Bomans 2008; Gvaramadze et al. 2011a; Peri et al. 2012). This method does not rely on the knowledge of the star’s proper motion, precise distance, or the characteristics of the local standard of rest (LSR) of the studied system. Stellar bow shocks are also parsec-scale structures, manifesting strongly in the mid-infrared, so they can be detected at considerable distances due to their size and low levels of extinction at these wavelengths (e.g., Gvaramadze et al. 2010, 2011b). Thus, this method could be especially potent when studying distant HMXBs, as most of them do not have their proper motions and distances determined or they are measured with low significance.

The peculiar radial velocity of HMXBs, their Galactic height distribution, and low OB association memberships provides a strong evidence that HMXBs are a subclass of high-velocity stars (van Oijen 1989). This conclusion was verified by the peculiar tangential velocity measurements by Chevalier & Ilovaisky (1998) and van den Heuvel et al. (2000). Both studies also note a disparity between the peculiar tangential velocities of OB-supergiant X-ray binaries and Be/X-ray binaries; the former group has a mean peculiar tangential velocity ⟨vpec,tan⟩ = 42 ± 14 km s-1, while the latter only has a modest ⟨vpec,tan⟩ = 15 ± 6 km s-1 (van den Heuvel et al. 2000). Despite the considerable difference in the peculiar velocity between the two subgroups, the faster members of the Be/X-ray binary subgroup should have sufficient peculiar velocity to produce an observable bow shock (Meyer et al. 2017).

The first detection of a bow shock associated with a HMXB was reported by Kaper et al. (1997), who discovered a faint arcuate structure around a well-known HMXB Vela X-1 in a Hα image obtained with the 1.54 m Danish telescope at the European Southern Observatory. The first bow shock survey of HMXBs was conducted by Huthoff & Kaper (2002), choosing a sample of 11 HMXBs for which the conditions for bow shock formation seemed favorable: a large peculiar velocity and an OB companion with a strong stellar wind and a high luminosity. Using high-resolution IRAS maps, Huthoff & Kaper (2002) detected an infrared counterpart to the already known bow shock generated by Vela X-1. Another survey was conducted by Gvaramadze et al. (2011c) on the same sample of HMXB systems. They utilized the data from the Spitzer Space Telescope, which were of higher quality but they covered only 5 of 11 HMXB from the original sample. Nevertheless, these authors were able to discover a new bow shock around 4U 1907+09 and provided a more detailed view of Vela X-1 bow shock as well.

We are currently aware of more than a hundred HMXBs or HMXB candidates in the Galaxy and their numbers are increasing (Liu et al. 2006). Together with the advent of space infrared missions such as Spitzer Space Telescope (Werner et al. 2004) and Wide-field Infrared Survey Explorer (WISE; Wright et al. 2010), this presents an opportunity to have a more detailed look at the bow shock occurence around HMXBs. In this paper, we have undertaken a search for bow shocks around the Galactic HMXBs from the fourth edition of the Catalog of High Mass X-ray binaries in the Galaxy (Liu et al. 2006) and around systems listed in Walter et al. (2015), which are not listed in Liu et al. (2006).

The paper is organized as follows. The rationale, observations, and general data processing are described in Sect. 2. A detailed description of the selected targets and supplementary astrometry is given, respectively, in Sects. 3 and 4. We discuss the nature of the emission around the studied objects in Sect. 5. Finally, the summary and conclusions are provided in Sect. 6.

2 Rationale and methods

A number of past successful studies of bow shocks conducted in the mid-infrared (e.g., Peri et al. 2012, 2015; Kobulnicky et al. 2016 and the references therein) and extensive computational work (e.g., Meyer et al. 2016, 2017) suggest that this part of the electromagnetic spectrum is the most suited for bow shock searches.

To search for bow shocks around HMXBs, we used data obtained by Multiband Imaging Photometer for Spitzer (MIPS; Rieke et al. 2004) extracted from the Spitzer Space Telescope archive and data obtained from WISE Image Service. The available Spitzer data cover the fields containing less than a half of catalogued HMXBs, the data from WISE cover all surveyed HMXB systems. To obtain science-grade images, the Spitzer data was processed by Mosaicking and Point source Extraction (MOPEX; Makovoz & Khan 2005). Visual inspection of Spitzer and WISE images was used to search for the extended emission possibly associated with the HMXBs. To search for the possible counterparts of the detected emission, we used the SuperCOSMOS H-alpha Survey (SHS; Parker et al. 2005), data obtained by the Photodetector Array Camera and Spectrometer (PACS; Poglitsch et al. 2010), and the Spectral and Photometric Imaging Receiver (SPIRE; Griffin et al. 2010) on board the Herschel Space Observatory (Pilbratt et al. 2010), and data obtained by the InfraRed Array Camera (IRAC; Fazio et al. 2004) on board the Spitzer Space Telescope.

3 Remarks on individual objects

Owing to the short lifetimes of HMXBs, these objects cannot travel far from the Galactic plane where they formed. This means that they are often projected onto a region with complex diffuse infrared emission and it is often hard to conclude whether the emission is actually physically related to the system or is of a foreground/background origin. However, a physical relation is likely if the diffuse emission is centered at the system or exhibits some degree of symmetry wherein the HMXB lies on or near the axis of symmetry. Below, we briefly describe the systems with such infrared emission.

3.1 γ Cas

γ Cas has attracted a lot of attention since its discovery as the first of what became known as Be stars (Secchi 1866). Despite considerable interest, the origin of X-ray emission from the object is still highly debated. The proposed scenarios include the accretion onto a neutron star or a white dwarf and amagnetic star–disk interaction (see review by Smith et al. 2016 and the references therein). The nature of this system is so peculiar that it has become a prototype of a separate class of X-ray emitters knowns as “γ Cas analogs”. While the recent study by Postnov et al. (2017) reconciles the peculiar X-ray properties of γ Cas by invoking a fast spinning neutron star as a companion, their results are disputed (Smith et al. 2017). Because of its atypical properties and the disputed nature, we consider γ Cas as a HMXB candidate only.

γ Cas is a relatively nearby object. Both Perryman et al. (1997) and Coleiro & Chaty (2013) put it at a similar distance of 0.19 and 0.17 kpc, respectively. Megier et al. (2009) estimated a lower distance of 0.12 kpc. We adopt a distance of 0.15 kpc as a compromise for the subsequent analysis.

WISE 22 μm image (Fig. 1) reveals a faint arc-like nebula to the north of the system. Unfortunately, the system is not covered by the Herschel data and Spitzer data are limited to the first two IRAC channels, taken in the post-cryo mode, which do not show any discernable diffuse emission. The system is slightly offset from the axis of the symmetry of the arc. The system is also embedded in a fainter 22 μm diffuse emission.

3.2 EXO 051910+3737.7 (V420 Aur)

The X-ray source was discovered by Uhuru satellite (Forman et al. 1978) and associated with a Be star V420 Aur by Polcaro et al. (1984). No X-ray pulses have been detected from the system (Liu et al. 2006).

The system does not have a significant HIPPARCOS parallax, Chevalier & Ilovaisky (1998) estimated the distance to be approximately 1.05 kpc using a typical luminosity of a B0 subgiant. Astraatmadja & Bailer-Jones (2016) estimated the distance to be 1.5 ± 0.5 kpc utilizingTycho-Gaia astrometric solution (TGAS; Lindegren et al. 2016) parallaxes and Milky Way prior. A more recent catalog by Bailer-Jones et al. (2018), employing the second Gaia data release (hereafter GDR2; Gaia Collaboration 2018), puts the system at 1.29-0.09+0.11$\mathrm{1.29_{-0.09}^{+0.11}}$ kpc.

The system appears to be a center of an irregular wispy infrared nebula (Fig. 2), readily apparent on WISE 12 μ m, 22 μ m, and MIPS 24 μm images. The central, brighter part of the nebula, is divided into two filaments, aligned in the Galactic north–south direction: one apparently connected with the system and one slightly ahead of it. The nebula also possesses four fainter arms; the most dominant arm extends toward the Galactic south and the fainter three point to the Galactic north, northwest, and northeast.

thumbnail Fig. 1

WISE 22 μm image of the field around γ Cas with 2MASS Ks contours (sensitive to starlight) overlaid on top. The system position is indicated by the cyan circle.
thumbnail Fig. 2

Spitzer MIPS 24 μm image of the field around EXO 051910+3737.7 (V420 Aur) with Spitzer IRAC 3.6 μ m contours (mainly sensitive to starlight) overlaid on top. The system position is indicated by the cyan circle.

3.3 IGR J11435-6109

IGR J11435-6109 is a HMXB discovered by Grebenev et al. (2004) via INTEGRAL mission. The system is an X-ray pulsar with a pulse period of about 162 s and orbital period of 52.5 days (in’t Zand & Heise 2004; Corbet & Remillard 2005). The optical counterpart was determined by Tomsick et al. (2007) and confirmed by Negueruela et al. (2007). Coleiro et al. (2013) refined the spectral classification to B0.5Ve, confirming the system as a Be/X-ray binary.

The system is a distant X-ray binary (>6 kpc; Negueruela et al. 2007). Masetti et al. (2009) estimated the distance to be about 8.6 kpc. Coleiro & Chaty (2013) derived a slightly higher distance, d = 9.8 ± 0.86 kpc. Bailer-Jones et al. (2018) confirmed the distant nature of this HMXB, placing it at 8.59−1.81+2.54$8.59_{-1.81}^{+2.54}$ kpc.

Figure 3 shows the enviroment around the system. IGR J11435-6109 appears to be projected onto the eastern part of a large cavity.

thumbnail Fig. 3

WISE 22 μm image of the field around IGR J11435-6109 with 2MASS K band contours (sensitive to starlight) overlaid on top. The system position is indicated by the cyan circle.

3.4 GX 301-02

GX 301-2 is an obscured HMXB system that is very bright in X-rays, thanks to the slow and very dense stellar wind of its massive, hyper giant companion (Kaper et al. 2006). The system contains a pulsar with a period of around 700 s in 41.5 days, highly eccentric orbit (Sato et al. 1986).

The distance to GX 301-2 was estimated by Kaper et al. (2006) to be in a range of 3–4 kpc. Coleiro & Chaty (2013) estimated the distance to be 3.1 ± 0.64 kpc via a spectral energy distribution procedure. The system is covered in TGAS, however, the measured parallax is not significant (0.34 ± 0.75), so the distance estimates by Astraatmadja & Bailer-Jones (2016) vary significantly based on the prior used. The distance estimate was refined after GDR2; Bailer-Jones et al. (2018) estimated the distance to the system d = 3.53−0.40+0.52$3.53_{-0.40}^{+0.52}$ kpc, confirming the previous estimates.

The system lies in a region of very bright and complex infrared emission (see Fig. 4). This system is embraced by a region of bright 24 μm emission to the Galactic south and a bright ridge of extended emission to the north. The PACS and SPIRE data suggest that the system could be enclosed in a partial far-infrared bubble, open at one side.

thumbnail Fig. 4

Spitzer MIPS 24 μm image of the field around GX 301-2 (BP Cru) with Spitzer IRAC 3.6 μm contours (mainly sensitive to starlight) overlaid on top. The system position is indicated by the cyan circle.

3.5 1H 1255-567

Listed as a HMXB candidate in Liu et al. (2006), but as a Be star in Simbad (Levenhagen & Leister 2006), the system is a part of a visual double with μ1 Cru (Hoffleit & Jaschek 1982).

1H 1255-567 appears to be a close system. Chevalier & Ilovaisky (1998) estimated its distance to be about 111 ± 8 pc. Chen et al. (2012) listed a slightly larger value of 125 ± 5 pc. More recently, Bailer-Jones et al. (2018) estimated the distance to be 112−3+2$112_{-3}^{+2}$ pc.

Figure 5 shows a WISE 22 μm image of the visual double with both members visible in far-infrared emission, enshrouded in a fine extended nebula. Thenebula bears a complex shape and has a system of fine filaments and wisps, emanating from μ1 Cru (the prominent unmarked source to the Galactic south of 1H 1255-567) and 1H 1255-567 as well.

3.6 GX 304-01

The HMXB GX 304-01 (4U 1258-61) was discovered by balloon X-ray observations in 1970 (McClintock et al. 1971). The optical counterpart of the X-ray source was identified by Mason et al. (1978), making V850 Cen the most likely counterpart. A subsequent spectral analysis by Parkes et al. (1980) showed the presence of a strong double-peaked emission in Hα and deep sharp absorption lines indicating a Be star. The most recent measurements give a more precise spectral classification that indicates the secondary is a B0.7 Ve star (Ziolkowski 2002). Using the data from Vela satellite, Priedhorsky & Terrell (1983) found a periodicity in the X-ray flux of 133 days and interpreted this as being the orbital period of a neutron star in an orbit around the Be star, the variability in the X-ray flux being caused by enhanced accretion onto the neutron star as it passes through the periastron. This establishes GX 304-01 as a HMXB with a Be star companion (i.e., a Be/X-ray binary).

Parkes et al. (1980) estimated the distance of the system to be 2.4 ± 0.5 kpc. Coleiro & Chaty (2013) provided an estimate of the distance of 1.3 kpc. Using GDR2, Bailer-Jones et al. (2018) derived a larger distance of 2.01−0.13+0.15$2.01_{-0.13}^{+0.15}$ kpc.

The 22 μm WISE image of the system (see Fig. 6) reveals a fine arcuate nebula that has a noticable dimple at the apex of the arc, close to the HMXB. Unfortunately, the system is not covered by the Spitzer data. The arc is also visible in PACS 70 μm image; no Hα counterpart was detected. The arc exhibits a clear symmetric morphology without any discernable clumps or finer structure. The infrared counterpart of GX 304-01 lies along the axis of symmetry of the arc, making it likely that this emission is associated with the HMXB system. Because of the apex dimple, limited resolution of WISE 22 μm channel, and close proximity of the nebula to the system, the HMXB seems to be immersed in the nebula.

thumbnail Fig. 5

WISE 22 μm image of the field around 1H 1255-567 (μ2 Cru) with 2MASS Ks contours (sensitive to starlight) overlaid on top. The system position is indicated by the cyan circle.
thumbnail Fig. 6

WISE 22 μm image of the field around GX 304-01 (V850 Cen) with 2MASS Ks contours (sensitive to starlight) overlaid on top. The system position is denoted by the cyan circle.

3.7 AX J1639.0-4642

The X-ray source AX J1639.0-4642 was discovered by the Advanced Satellite for Cosmology and Astrophysics (ASCA) observatory (Sugizaki et al. 2001). The system is likely a supergiant HMXB. It is highly absorbed X-ray pulsar with a spin period of 911 s (Bodaghee et al. 2006). The infrared counterpart identification has been problematic. 2MASS J16390535-4642137 was proposed as a possible counterpart (Chaty et al. 2008). Further refinement of the system position by Bodaghee et al. (2012a) suggested another close, faint, and blended star as the most probable counterpart. The precise orbital period of the system is disputed, but is most likely around 4.2 days (e.g., Corbet et al. 2010; Corbet & Krimm 2013). These values of pulse and orbital period place the system among the other wind-fed HMXBs in Corbet’s diagram of pulse versus orbital period (Corbet 1986).

The faintness of the infrared counterpart hampers the distance determination. Assuming a luminosity typical for an accretion-driven X-ray pulsar, Bodaghee et al. (2006) estimated the distance of the source to be approximately 10 kpc.

Figure 7 shows the field around the system. The HMXB is adjacent to a bright 24 μ m nebula, centered on the group of faint infrared stars. The HMXB is immersed into the fainter outskirts of this nebula. The nebula also dominates the field at IRAC and Herschel wavelengths.

thumbnail Fig. 7

Spitzer MIPS 24 μm image of the field around AX J1639.0-4642 with Spitzer IRAC 3.6 μm contours (mainly sensitive to starlight) overlaid on top. The system position is denoted by the cyan circle.

3.8 4U 1700-37

4U 1700-37 is an eclipsing HMXB with a very massive companion of O6Iafcp spectral class (Jones et al. 1973; Sota et al. 2014). The nature of the compact object is unclear. The absence of pulsation suggests the presence of a low-mass black hole, but a neutron starseems more likely (e.g., Reynolds et al. 1999; Boroson et al. 2003).

Ankay et al. (2001) estimated the distance to be about 1.9 kpc. Megier et al. (2009) and Coleiro & Chaty (2013) adopted 2.12 ± 0.34 and 1.8 ± 0.15 kpc, respectively. Depending on the adopted prior, Astraatmadja & Bailer-Jones (2016) derived vastly different distances from the TGAS data (approximately from 0.7 to 3.3 kpc), which also suffer from large uncertainties. Utilizing GDR2, Bailer-Jones et al. (2018) derived a more precise distance estimate of 1.75−0.19+0.24$1.75_{-0.19}^{+0.24}$ kpc.

Figure 8 shows that the system lies on the axis of symmetry of an arcuate cavity that is excavated into a nearby bright cloud of 24 μm.

thumbnail Fig. 8

WISE 22 μm image of the field around 4U 1700-37 (V884 Sco) with WISE 3.4 μm contours (mainly sensitive to starlight) overlaid on top. The system position is denoted by the cyan circle.

3.9 EXO 1722-363

EXO 1722-363 was discovered in the Galactic ridge by EXOSAT observations in 1984 (Warwick et al. 1988). Analysis of Ginga observations by Tawara et al. (1989) suggested the presence of a dense envelope of circumstellar matter around the system. Takeuchi et al. (1990) placed a 9 days lower limit for the orbital period, which was later refined by Corbet et al. (2005) and Thompson & Tomsick (2007) to be 9.74 days. Using INTEGRAL observations, Zurita Heras et al. (2006) found a probable infrared counterpart to the system, 2MASS J17251139-3616575, 1′′ away from the best position of the source as given by INTEGRAL. Subsequent observations of the counterpart by Rahoui et al. (2008) and Mason et al. (2009) supported the earlier assumptions that the system is a HMXB containing an early supergiant B star, which produces strong stellar wind that fuels the accretion onto the neutron star.

The absence of the visual counterpart and the location of the system in the Galactic ridge region make the distance determination difficult. Zurita Heras et al. (2006), considering a typical luminosity of an active accretion pulsar, estimated the distance to the system to be 7 kpc. Thompson & Tomsick (2007) placed the system between 5.3 and 8.7 kpc. Using the broadband data, Rahoui et al. (2008) derived a slightly smaller distance to the system of 6.1 kpc. Distance determinations were also made by Mason et al. (2009), later refined by Mason et al. (2010), placing the system to be 7.1 ≤ d ≤ 7.9 kpc.

The nebula around EXO 1722-363 is visible in both the Spitzer MIPS 24 μ m (see Fig. 9) and 22 μm WISE images. The arc is also visible in the PACS 70 μm image and a gleam of emission possibly associated with the arc is also visible in the IRAC 8 μ m image. No Hα counterpart was detected. The Spitzer MIPS 24 μm image reveals an irregular, curved morphology, where the part of the nebula to the Galactic east of the system appears to be brighter andmore prominent. This seems to be a feature of the structure itself, as our inspection of this area in the different bands does not reveal any possible background or foreground sources that could be responsible for the enhanced emission. The infrared counterpart of EXO 1722-363 appears to lie within the nebula. If we consider the outer boundary of the nebula tracing an arc, the system lies along the approximate symmetry axis of this arc. This, together with the lack of any notable brighter point sources within the nebula, suggests that the structure is likely associated with the HMXB.

thumbnail Fig. 9

Spitzer MIPS 24 μm image of the field around EXO 1722-363 with Spitzer IRAC 3.6 μm contours (mainly sensitive to starlight) overlaid on top. The system position is indicated by the cyan circle.

3.10 XTE J1739-302

XTE J1739-302 (IGR J17391-3021) is a supergiant X-ray transient (SFXT) discovered with the proportional counter array on board the Rossi X-Ray Timing Explorer (RXTE; Smith et al. 1998). The system consists of a neutron star on an eccentric 51.47 days orbit (Drave et al. 2010) around a O8 Iab(f) supergiant companion (Negueruela et al. 2006a, Rahoui et al. 2008).

The source lies in the direction of the Galactic center. Using spectroscopy and photometry of the counterpart, Negueruela et al. (2006a) estimated its distance to be 2.3 ± 0.6 kpc. The subsequent analysis by Rahoui et al. (2008) derived a slightly larger distance of about 2.7 kpc. The distance estimate from GDR2 seems to be less constraining. Bailer-Jones et al. (2018) provided an estimate for the distance of 5.32−2.12+3.92$5.32_{-2.12}^{+3.92}$ kpc. It is therefore questionable whether the HMXB is a foreground source, as suggested by the earlier studies, or is, potentially, a Galactic center object.

Figure 10 shows that the system lies in a region of complicated mid-infrared emission, typical for Galactic center region. The system could be lying on the axis of symmetry of a faint arc of extended emission or a partial bubble.

thumbnail Fig. 10

Spitzer MIPS 24 μm image of the field around XTE J1739-302 with Spitzer IRAC 3.6 μm contours (mainly sensitive to starlight) overlaid on top. The system position is indicated by the cyan circle.

3.11 AX J1841.0-0536

AX J1841.0-0536 was discovered as a variable X-ray source by ASCA (Bamba et al. 2001) while undergoing two bright flares lasting approximately 1 h each. The flaring activity of the system puts it into a SFXT class of HMXBs (e.g., Romano et al. 2011). Chandra observations by Halpern et al. (2004) pinpointed the infrared counterpart. Negueruela et al. (2006b) provided constraints on the spectral type of the secondary, being of a luminous B0-1 type. This is further refined by the infrared spectroscopy done by Nespoli et al. (2008), who derived the spectral type B1Ib for the secondary. The orbital period of the system of 6.45 days was reported by González-Galán (2015).

The distance estimates of the system vary. Nespoli et al. (2008) provided a broad estimate of the distance of 3.2−1.5+2$3.2^{+2}_{-1.5}$ kpc. A more recent analysis by Coleiro & Chaty (2013) put the system further away, at 7.8 ± 0.74 kpc. This is in agreement with the distance estimated from GDR2, d = 7.60−2.16+3.06$7.60_{-2.16}^{+3.06}$ kpc (Bailer-Jones et al. 2018).

Figure 11 shows that the system lies in the center of the region of reduced 24 μ m emission. The cavity does not feature bright rims of mid-IR emission and the evidence of the reduced emission can be seen on the PACS and SPIRE wavebands as well.

thumbnail Fig. 11

Spitzer MIPS 24 μm image of the field around AX J1841.0-0536 with Spitzer IRAC 3.6 μm contours (mainly sensitive to starlight) overlaid on top. The system position is indicated by the cyan circle.

3.12 XTE J1855-026

XTE J1855-026 is a HMXB discovered by RXTE. The system showed pulsations with a period of about 361 s and a flux modulation with a period of 6.1 days, which was interpreted as the orbital period of the system (Corbet et al. 1999). This was confirmed by Negueruela et al. (2008). The optical counterpart of the system was pinpointed by Verrecchia et al. (2002) and its spectral type was determined to be B0Iaep (Negueruela et al. 2008).

Corbet et al. (1999) proposed a distance of approximately 10 kpc. The analysis by Coleiro & Chaty (2013) derived a consistent distance estimate of 10.8 ± 1.0 kpc. The distance given by Bailer-Jones et al. (2018) is slightly smaller, 8.14−1.79+2.58$8.14_{-1.79}^{+2.58}$ kpc.

This system also appears to be associated with a region of reduced 24 μ m emission (Fig. 12). The paucity of emission is not that well pronounced as is the case with the previous system, however, the region appears to be surrounded by a bright partial mid-infrared shell, most noticable to the Galactic south and northwest ofthe system.

thumbnail Fig. 12

WISE 22 μm image of the field around XTE J1855-026 with WISE 3.4 μm contours (mainly sensitive to starlight) overlaid on top. The system position is indicated by the cyan circle.

4 Kinematics of the studied systems

The kinematics of the studied HMXBs is of significant interest, as it can provide constraints for several interesting system characteristics, such as the origin of the system or the properties of the progenitor binary. However, studying it has been difficult owing to considerable distances of these sources and insufficient accuracy of the past astrometric surveys. This resulted in the proper motions measured with a low significance, often varying significantly depending on the survey, making them hard to utilize. Only a few, very close sources did not suffer from these issues. The situation improved somewhat with the advent of the first Gaia data release (GDR1; Gaia Collaboration 2016). While TGAS, part of GDR1 (Lindegren et al. 2016), contained only a small number of HMXBs, re-reductions of past catalogs using GDR1, such as UCAC5 (Zacharias & Finch 2017) and Hot Stuff for One Year (HSOY; Altmann et al. 2017) provided more precise proper motions for many more sources. However, the issue of discrepant proper motions depending on the catalog used still persisted.

We used the parallaxes and proper motions from GDR2 (Gaia Collaboration 2018) to compute peculiar tangential velocities, which are especially interesting for systems associated with an arcuate structure in their vicinity. For our objects, we followed the approach outlined in Luri et al. (2018). The distances and tangential velocities were jointly estimated from the parallaxes and proper motions via Bayesian inferrence and the prior scale lengths are adopted from Bailer-Jones et al. (2018). To obtain the peculiar tangential velocities, we adopted the Galactic constants R0 =8.2 kpc, Θ0 =238 km s−1 and the solar peculiar motion (U⊙, V⊙, W⊙) = (10.0, 11.0, 7.0) km s−1 from Bland-Hawthorn & Gerhard (2016). Table 1 lists the computed peculiar velocity medians with the errors obtained from 68% credibility intervals. The only exception was γ Cas, as it is too bright for Gaia. In this case, we used the proper motion given in HSOY with an adopted distance of 0.15 kpc as discussed above. It must be noted that we computed the errors in the peculiar velocity for γ Cas from the errors of the proper motion measurements only, so they should be considered as low-limit approximations.

Unfortunately, for the rest of the systems not listed in Table 1, the choice of catalogs is limited. Because of the significant distance of these sources or faint optical/infrared counterparts, they are not listed in any proper motion catalogs or they exhibit no measurable proper motion.

A number of the studied objects exhibit increased values of the astrometric flags, such as the astrometric excess noise and goodness-of-fit statistic, possibly indicating problems with the astrometric solution in GDR2. Issues with the astrometry may arise when dealing with the regions with large source densities, such as the Galactic plane and center regions, where HMXBs are predominantly situated. Another caveat is that the astrometric solution in GDR2 does not take the binarity into account, which may have an impact on the astrometry and its quality (Lindegren et al. 2018). To investigate these effects, we queried GDR2 for stars of similar magnitude and color (ΔG = ±0.5 mag, Δ(BP - RP) = ±0.2 mag) within a 10° radius cone centered on each source. For each source we also constructed an equivalent query, extracting the sources from the opposite side of the sky to study the effects of crowding. Comparison of the studied HMXBs with the extracted stars from the source vicinity and the opposite side of the sky showed that, except for AX J1841.0-0536, none of the HMXBs are significant outliers, despite some sources having astrometric flag values indicating problems with the single-star astrometric solution. It is, however, unlikely that the measured parallaxes and proper motions are corrupted by the unmodeled orbital motion. The longest orbital period among our sources is approximately 133 days for GX 304-01 (Priedhorsky & Terrell 1983). Even this orbital period is small compared to the 22-month observing time of GDR2. For such periods, it is expected that the orbital motion should average out and should not significantly impact the parallax and proper motion measurements (Jennings et al. 2018). Interestingly, AX J1841.0-0536, the system with elevated astrometric flags as compared to the extracted sources from the queries, only has an orbital period of 6.45 days (González-Galán 2015).

Table 1

Proper motions and derived peculiar tangential velocities for selected HMXBs.

5 Discussion

The infrared survey of the HMXB sample yielded a variety of extended emission that may suggest the influence of the HMXBs on the surrounding interstellar medium (ISM). This emission, its presence or absence, has important implications not only for the particular system, but also for the HMXB population as a whole.

5.1 Extended emission

The studied objects are projected against complex infrared structures, which are not only arcs bending outward from the HMXB, possibly indicating a bow shock, but also infrared filaments that appear to emanate from the system and emission cavities.

5.1.1 γ Cas arc

Figure 1 shows the arcuate nebulosity in the vicinity of γ Cas. While the infrared arcs around early-type stars are traditionally interpreted as stellar bow shocks, various phenomena can produce morphologies resembling bow shocks, including partial infrared bubbles and H II regions with density gradients. The infrared arcs may also be interpreted as dust waves or bow waves that may form around late-type OB main sequence stars, which may drive weaker stellar winds than expected (the weak wind problem; Ochsendorf et al. 2014a,b). These structures can be attributed to the interaction of radiation pressure from the star with dust carried along by a photoevaporative flow of ionized gas from a dark cloud or inside an evolved interstellar bubble. The structures are very similar in morphology to genuine stellar bow shocks.

The position of the arc is approximately consistent with the position of GCRV 309 E, a H II region, however, we were unable to confirm its presence on the image data in the INT Photometric H-Alpha Survey (IPHAS; Barentsen et al. 2014) Hα images. The absence of other infrared data makes it hard to establish the exact nature of the nebula. However, given the peculiar velocity of the source, its spectral type, a bow or a dust wave interpretation is more likely than a stellar bow shock. Adding to this, the arc approximately points to the infrared pillars, associated with two clouds, IRAS 00560+6037 and IC 59, which may be sources of a material flow toward the system. While the nature of the arc is debatable, it is, most likely, not a stellar bow shock.

5.1.2 EXO 051910+3737.7 nebula

The wispy nebula, fanning approximately in the Galactic north–south direction, may suggest a dusty outflow indicating a systemic mass loss. The nebula bears some resemblance to circumstellar nebulae observed around some B[e] and Be stars in the infrared images (Mayer et al. 2016). The literature on the enviroments around Be stars is sparse, but at least one Be star, SX Aur, also shows the presence of a mid-infrared nebula (Deschamps et al. 2015), however, this nebula is more compact and shows a blob-like morphology instead of fine wispy jets.

There also exists a possibility that the dusty structure is not associated with the system and the system is just passing in its vicinity, heating the dust. However, the position of the system in the central knot of the emission hints at the association. Moreover, the system has a significant peculiar velocity to the Galactic east. If the outflows emanate from the system, it would be expected to observe some bending of the structure toward the west, as the system plows through the surrounding ISM. This bending is apparent in the bright part of the nebula and also in the fainter south wisp.

The nature of the easternmost rim of the nebula is puzzling. It is leading the star in the direction of its peculiar motion, so it cannot be related to the past mass-loss episode. It also exhibits the same degree of bending as the bright bar centered on the system. A bending toward the system due to its systemic peculiar motion would be expected fora stellar bow shock, however, the bar does not seem to be completely detached from the system, which is not expected. The bar may be a part of an infrared filament, locally heated by the system as it passes in its vicinity.

5.1.3 IGR J11435-6109 cavity

The systemlies in a large infrared cavity that is apparent in WISE 12, 22 μm images and in SPIRE wavebands. Coleiro et al. (2013) discussed the cavity briefly and noted that it could be a result of the feedback of the system on the surrounding ISM. The physical relation of the system to the cavity is difficult to prove because of the inherent inaccuracy of the distance determination to both infrared cavities and distant HMXBs. Moreover, the system lies significantly off-center of the cavity. This could be reconciled if the system has a significant peculiar velocity pointing away from the cavity center. The system has proper motion data available but they are of a low quality. However, the peculiar motion component in the Galactic latitude is sufficiently well constrained to suggest that the system is moving to the Galactic south, which is at odds with its supposed origin near the center of the cavity. Considering the points above, it seems that the HMXB and the cavity it is projected against are unrelated.

5.1.4 GX 301-02 nebula

The system is surrounded by an extended infrared emission. Huthoff & Kaper (2002) and Servillat et al. (2014) discussed the possibility of a bow shock or a cavity associated with the system. The peculiar velocity obtained from the proper motion suggests that the system is moving predominantly to the Galactic southeast toward the bright extended emission, so the bright rim of infrared emission to the north of the system is unlikely to be a stellar bow shock. The surrounding emission is most likely related to the infrared bubble [SPK2012] MWP1G300134-001035 to the south of the system. No bow shock is apparent despite the considerable peculiar velocity of the system.

5.1.5 1H 1255-567 nebula

The system is enshrouded in a fine wispy nebula, similar to the wispy nebula associated with EXO 051910+3737.7, but less prominent. Moreover, it is hard to determine if the emission is instead associated with μ1 Cru. The system exhibits no significant peculiar motion, yet the filaments seem to exhibit some curvature, as is evident in Fig. 5. However, also the southern star, μ1 Cru, does not exhibit any significant peculiar systemic motion.

5.1.6 GX 304-01 arc

The HMXB appears to be associated with a mid-infrared smooth arc pointing approximately to the Galactic northeast. The proper motion of the system from the GDR2 data yields a mildly runaway peculiar velocity primarily to the Galactic east, which deviates about 30° from the approximate symmetry axis of the arc. The arcuate mid-infrared nebula upstream of the system could be interpreted as a stellar bow shock. However, there are several problems with this interpretation. The arc does not bear a classical bow shock shape and does not seem to be fully detached from the system. The most puzzling feature of the nebula is its apex dimple, pointing toward the system. Recently, Meyer et al. (2017) conducted a series of simulations of stellar bow shocks of early-type stars in a magnetized medium. Their modeling suggests that in a magnetized ambient medium, the classic bow shock shape gets distorted and compressed, increasing its opening angle and becoming much blunter around the apex, as the stand-off distance decreases significantly. Interestingly, an apex dimple may form. This bow shock morphology change is especially prominent for stars having a modest space velocity, matching that of GX 304-01. This makes the bow shock interpretation appealing, however, the apparent attachment of the emission onto the system still remains an issue despite the expected decreased stand-off distance expected for such system. The mentioned misalignment between the peculiar velocity vector and the symmetry axis of the arc also poses a problem for the bow shock hypothesis, however, this can be reconciled by the presence of a large-scale flow in the ISM. An alternative origin of the emission might be due to the system encountering an infrared filament along its way, locally heating and compressing it. An example of such system is HD 49662, as discussed in Kalas et al. (2002). This system appears to be embedded in an infrared filament, heating it locally. This produces a blister-like infrared emission, visible in WISE 12 and 22 μm images, bulging ahead of the system while the system is embedded in the diffuse emission. However, the peak of the diffuse emission is centered directly on the system, which is not the case for the HMXB. Also, the diffuse emission around GX 304-01 does not bear a blister-like shape. Thus, we classify the emission around GX 304-01 as ambiguous.

5.1.7 Blob near AX J1639.0-4642

AX J1639.0-4642 is adjacent to a blob of strong mid-infrared emission. The emission is very concentrated, however, its outskirts reach the HMXB, which is the most apparent in the longer IRAC wavebands. Several infrared point sources are present in the central part of the blob. This region is designated IRAS 16353-4636 and is a site of star formation (Benaglia et al. 2010). The point sources within the blob make up a protostellar cluster. Benaglia et al. (2010) derived a distance of ~8 kpc to the embedded protocluster. Owing to the inherent evolved nature of HMXBs, and because this protocluster is about 2 kpc closer than the HMXB, this cluster is not related to AX J1639.0-4642.

5.1.8 4U 1700-37 cavity

The system seems to be centered on an arcuate cavity protruding into a nearby mid-infrared cloud to the Galactic southeast. Toalá et al. (2017) suggested that this structure is a stellar bow shock driven by the HMXB. The HMXB has a well-measured proper motion, which yields a runaway peculiar velocity at the adopted distance. The inferred direction of the peculiar motion deviates by about 60° from the approximate axis of symmetry of the cavity. This difference cannot be attributed to the errors in the peculiar velocity, as it is well constrained. Also, considering the magnitude of the peculiar velocity of the system of about 70 km s−1, the peculiar ambient medium velocity would have to be considerably large to produce such deviation. Another point to consider is the absence of the emission enhancement along the boundary of the cavity. The brightness along the rim of the tentative bow shock is practically the same as in the surrounding cloud to the Galactic southeast. This poses problems for the bow shock interpretation. A possible alternative could be a partial cavity, possibly shaped by the feedback from the system, or a chance alignment of a foreground/background structure.

5.1.9 EXO 1722-363 nebula

The system is projected atop a crescent-shaped irregular mid-infrared nebula. Owing to a considerable distance to the system, there is no proper motion information available, so it is not possible to investigate whether there exists a connection with the orientation of the nebula and the systemic peculiar motion. The nebula is not evenly bright; as can be seen in Fig. 9, the east part of the nebula is significantly brighter. The brightening could be related to the systemic motion, however, it is impossible to ascertain given the lack of data. The system may be passing through a larger ISM cloud and shock only a part of it, while heating some of the unshocked material as well, which could explain the fainter filament that is projected downstream of the tentative bow shock. The nebula could be interpreted as a partial wind-blown bubble, however, the system is lying significantly off center. Another possibility is that the filament projected onto the system is a part of an interstellar cloud crossing behind or in front of the HMXB, passing closest to the system on its east-projected side. We classify this structure as ambiguous.

5.1.10 XTE J1739-302 filament

The system is adjacent to a fine mid-infared arc or filament, situated to the Galactic east, and its approximate axis of symmetry is oriented in the same direction. The proper motion data for the system yield a runaway peculiar velocity in the direction approximately to the Galactic northeast, albeit with considerable errors (see Table 1). For this reason, it is not meaningful to investigate the alignment of the peculiar velocity vector with the symmetry axis of the arc. Adding to this, the arcuate filament is not significantly brighter than the surrounding emission to the Galactic south that it seems to be linked with, as would be expected for the heated dust piled at the bow shock front. This implies that the filament cannot be interpreted as a stellar bow shock, and is, most likely, not related to the HMXB system.

5.1.11 AX J1841.0-0536 cavity

The system is projected into a cavity of a reduced 24 μm emission. Taking the peculiar velocity at face value, the system moves rather slowly for a supergiant-hosting HMXB. However, because of its uncertain distance, its peculiar velocity is ill-determined. Interestingly, the system seems to be a part of a small star group, possibly making up a star cluster. The distance to the cavity is not known, as it is not catalogued. However, the presence of a HMXB together with the surrounding stars in its center, make it possible that the cavity is shaped by stellar feedback.

5.1.12 XTE J1855-026 structure

The system seems to be located in a region of reduced 24 μm emission, bracketed by regions of stronger emission to the Galactic south and northwest, possibly making up a partial bubble. The system is a likely runaway owing to its position away from the Galactic plane. The precise kinematics of the system is hard to constrain because of the errors in the parallax measurement, however, none of the bright rims on the either side of the HMXB can be interpreted as stellar bow shocks.

5.2 Paucity of bow shocks around HMXBs

Our search for bow shocks driven by HMXBs yielded only one questionable new detection. Even if we consider the structure around GX 304-01 as a bona fide stellar bow shock, together with the bow shocks associated with Vela X-1 and 4U 1907+09, there are only three known HMXB bow shocks (only about 2% of the HMXBs drive detectable bow shocks). This is considerably lower than the bow shock occurrence around OB runaway stars, although this number has been fluctuating over the years (about 13%, van Buren & Noriega-Crespo 1995; ~40%, Huthoff & Kaper 2002; ~10%, Peri et al. 2012; and ~6%, Peri et al. 2015). This discrepancy can be explained by the hypothesis proposed by Huthoff & Kaper (2002) suggesting different ejection scenarios for HMXBs and OB runaways. OB runaways have predominantly escaped their parent associations or clusters at a relatively early stage, when the cluster was dense and the probability of close encounters and ejections of stars was high. On the other hand, a lot of HMXBs became runaways only after the occurrence of a supernova within the system. Thus, the kinematical age and the distance to their parent clusters should be lower for HMXBs. Therefore, HMXBs are more likely to be still enclosed in the hot and rarefied regions (hot bubbles) that surround OB associations and clusters while the OB runaways have already escaped from these regions. The (isothermal) speed of sound in the ISM increases with temperature, thus, the hotter the medium is, the less likely are runaways moving through it supersonically. Together with the lower ISM density in these regions, this means a smaller chance of bow shock detection. At the time, Huthoff & Kaper (2002) were only able to study the difference in bow shock occurrence between the OB runaways and HMXBs on a sample of 11 HMXB systems. The analysis of the clustering between HMXBs and OB associations by Bodaghee et al. (2012b) also reinforced this notion, suggesting that the massive binaries that are the progenitors of HMXBs tend to remain gravitationally bound to their birth sites until the supernova explosion in the system. Thus, they acquire their runaway velocity only later on, leaving the association or cluster after it has evolved considerably, evacuating cavities in their surroundings, making the formation of an observable bow shock less likely.

Moreover, in addition to bow shocks associated with stars having large peculiar velocities, the young stars in a star cluster can drive significant outflows, causing the local ISM velocities to deviate from the LSR. Thus, bow shocks can also be generated around neighboring stars or stars at cluster outskirts without necessarily having large peculiar velocities (e.g., Povich et al. 2008). These bow shocks are often coined “in situ bow shocks” and their axes of symmetry point predominantly toward the source of outflows. The effects of outflows abates as the cluster ages and thus by the time the first supernova explosions start to occur in a cluster, possibly giving rise to HMXBs or their progenitiors, the outflows are weaker and the ISM around the cluster is more rarified because of feedback effects. Thissuggests that the in situ bow shock generation is also less likely for HMXBs as compared to other OB stars.

All this points to the scenario in which massive binaries, which are progenitors to HMXBs, are not efficiently expelled from star clustersand associations during the early stages of a cluster lifetime. This also implies that the two-step ejection process (Pflamm-Altenburg & Kroupa 2010) is not at work for a considerable number of HMXBs, making it possible to determine their birth sites via the peculiar velocity measurements.

6 Summary and conclusions

We have searched for bow shocks around HMXBs using WISE and MIPS Spitzer Space Telescope archival data. Apart from the already known bow shocks (associated with Vela X-1 and 4U 1907+09), we found only one new structure resembling a bow shock, but even this structure cannot be conclusively interpreted as a bona fide bow shock.

The detection of the bow shock candidate around GX 304-01 suggests that this system possesses at least a moderate runaway velocity. The proper motion measurements also seem to support the runaway interpretation of GX 304-01. The bow shock around GX 304-01, if real, would be the first bow shock detected associated with a Be/X-ray binary.

The relative paucity of bow shocks associated with HMXBs as compared to OB runaway stars supports the hypothesis that HMXBs are kinematically younger objects than OB runaways, meaning that most of them are still moving within hot and tenuous medium near their parent cluster or association and thus unable to form an observable bow shock.

We expect that future releases of proper motion and paralax measurements obtained by the mission Gaia will considerably improve our knowledge of kinematics of these HMXBs. The planned James Webb Space Telescope will also be able to obtain more detailed images of bow shock candidates and nebulosities, allowing us to gain more insight into both the diffuse structures and the HMXBs that are associated with them.

Acknowledgements

The author acknowledges support from grant MUNI/A/0789/2017. This research has made use of the NASA/IPAC Infrared Science Archive, which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. This work is based (in part) on observations made with the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA. This publication makes use of data products from the Wide-field Infrared Survey Explorer, which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, funded by the National Aeronautics and Space Administration. This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC; https://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement. This research has made use of the SIMBAD database and the VizieR catalog access tool, both operated at CDS, Strasbourg, France. This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation. This paper makes use of data obtained as part of the INT Photometric Hα Survey of the Northern Galactic Plane (IPHAS; www.iphas.org) carried out at the Isaac Newton Telescope (INT). The INT is operated on the island of La Palma by the Isaac Newton Group in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias. All IPHAS data are processed by the Cambridge Astronomical Survey Unit, at the Institute of Astronomy in Cambridge. The band-merged DR2 catalog was assembled at the Centre for Astrophysics Research, University of Hertfordshire, supported by STFC grant ST/J001333/1. This research made use of Montage. It is funded by the National Science Foundation under Grant Number ACI-1440620, and was previously funded by the National Aeronautics and Space Administration’s Earth Science Technology Office, Computation Technologies Project, under CooperativeAgreement Number NCC5-626 between NASA and the California Institute of Technology. This research made use of Astropy, a community-developed core Python package for Astronomy (Astropy Collaboration 2018).

References

  1. Altmann, M., Roeser, S., Demleitner, et al. 2017, A&A, 600, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Ankay, A., Kaper, L., de Bruijne, J. H. J., et al. 2001, A&A, 370, 170 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Astraatmadja, T. L., & Bailer-Jones, C. A. L. 2016, ApJ, 833, 119 [NASA ADS] [CrossRef] [Google Scholar]
  4. Astropy Collaboration (Price-Whelan, A. M., et al.) 2018, AJ, 156, 123 [Google Scholar]
  5. Bailer-Jones, C. A. L., Rybizki, J., Fouesneau, M., et al. 2018, AJ, 156, 58 [NASA ADS] [CrossRef] [Google Scholar]
  6. Bamba, A., Yokogawa, J., Ueno, M., et al. 2001, PASJ, 53, 1179 [NASA ADS] [CrossRef] [Google Scholar]
  7. Barentsen, G., Farnhill, H. J., Drew, J. E., et al. 2014, MNRAS, 444, 3230 [NASA ADS] [CrossRef] [Google Scholar]
  8. Benaglia, P., Ribó, M., Combi, J. A., et al. 2010, A&A, 523, A62 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Berger, D. H., & Gies, D. R. 2001, ApJ, 555, 364 [NASA ADS] [CrossRef] [Google Scholar]
  10. Bland-Hawthorn, J., & Gerhard, O. 2016, ARA&A, 54, 529 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Bodaghee, A., Walter, R., Zurita Heras, J. A., et al. 2006, A&A, 447, 1027 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Bodaghee, A., Rahoui, F., Tomsick, J. A., et al. 2012a, ApJ, 751, 113 [NASA ADS] [CrossRef] [Google Scholar]
  13. Bodaghee, A., Tomsick, J. A., Rodriguez, J., et al. 2012b, ApJ, 774, 108 [NASA ADS] [CrossRef] [Google Scholar]
  14. Boroson, B., Vrtilek, S. D., Kallman, T., et al. 2003, ApJ, 592, 516 [NASA ADS] [CrossRef] [Google Scholar]
  15. Chaty, S., Rahoui, F., Foellmi, C., et al. 2008, A&A, 484, 783 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Chen, C. H., Pecaut, M., Mamajek, E. E., et al. 2012, ApJ, 756, 133 [NASA ADS] [CrossRef] [Google Scholar]
  17. Chevalier, C., & Ilovaisky, S.A. 1998, A&A, 330, 201 [NASA ADS] [Google Scholar]
  18. Coleiro, A., & Chaty, S. 2013, ApJ, 764, 185 [NASA ADS] [CrossRef] [Google Scholar]
  19. Coleiro, A., Chaty, S., Zurita Heras, J. A., et al. 2013, A&A, 560, A108 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Corbet, R. H. D. 1986, MNRAS, 220, 1047 [NASA ADS] [CrossRef] [Google Scholar]
  21. Corbet, R. H. D., & Krimm, H. A. 2013, ApJ, 778, 45 [NASA ADS] [CrossRef] [Google Scholar]
  22. Corbet, R. H. D., & Remillard, R. 2005, ATel, 377 [Google Scholar]
  23. Corbet, R. H. D., Marshall, F. E., Peele, A. G., et al. 1999, ApJ, 517, 956 [NASA ADS] [CrossRef] [Google Scholar]
  24. Corbet, R. H. D., Markwardt, C. B., & Swank, J. H. 2005, ApJ, 633, 377 [NASA ADS] [CrossRef] [Google Scholar]
  25. Corbet, R. H. D., Krimm, H. A., Barthelmy, S. D., et al. 2010, ATel, 2570 [Google Scholar]
  26. Deschamps, R., Braun, K., Jorissen, A., et al. 2015, A&A, 577, A55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Drave, S. P., Clark, D. J., Bird, A. J., et al. 2010, MNRAS, 409, 1220 [NASA ADS] [CrossRef] [Google Scholar]
  28. Fazio, G. G., Hora, J. L., Allen, L. E., et al. 2004, ApJS, 154, 10 [NASA ADS] [CrossRef] [Google Scholar]
  29. Forman, W., Jones, C., Cominsky, L., et al. 1978, ApJS, 38, 357 [NASA ADS] [CrossRef] [Google Scholar]
  30. Gaia Collaboration (Brown, A. G. A., et al.) 2016, A&A, 595, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Gaia Collaboration (Brown, A. G. A., et al.) 2018, A&A, 616, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Grebenev, S. A., Ubertini, P., Chenevez, J., et al. 2004, ATel, 350 [Google Scholar]
  33. González-Galán, A. 2015, ArXiv e-prints [arXiv:1503.01087] [Google Scholar]
  34. Griffin, M. J., Abergel, A., Abreu, A. et al. 2010, A&A, 518, L3 [Google Scholar]
  35. Gvaramadze, V. V., & Bomans, D. J. 2008, A&A, 490, 1071 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Gvaramadze, V. V., Kroupa, P., & Pflamm-Altenburg, J. 2010, A&A, 519, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Gvaramadze, V. V., Kniazev, A. Y., Kroupa, P., et al. 2011a, A&A, 535, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Gvaramadze, V. V., Pflamm-Altenburg, J., & Kroupa, P. 2011b, A&A, 525, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Gvaramadze, V. V., Roser, S., Scholz, R. -D., et al. 2011c, A&A, 529, A14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Halpern, J. P., Gotthelf, E. V., Helfand, D. J., et al. 2004, ATel, 289, 1 [NASA ADS] [Google Scholar]
  41. Hoffleit, D., & Jaschek, C. 1982, Bright Star Catalogue (New Haven: Yale University) [Google Scholar]
  42. Huthoff, F., & Kaper, L. 2002, A&A, 383, 999 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. in’t Zand, J., & Heise, J. 2004, ATel, 362 [Google Scholar]
  44. Jennings, R. J., Kaplan, D. L., Chatterjee, S., et al. 2018, ApJ, 864, 26 [NASA ADS] [CrossRef] [Google Scholar]
  45. Jones, C., Forman, W., Tananbaum, H., et al. 1973, ApJ, 181, L43 [NASA ADS] [CrossRef] [Google Scholar]
  46. Kalas, P., Graham, J. R., Beckwith, S. V. W., et al. 2002, ApJ, 567, 999 [NASA ADS] [CrossRef] [Google Scholar]
  47. Kaper, L., van Loon, J.Th., Augusteijn, T., et al. 1997, ApJ, 475, L37 [NASA ADS] [CrossRef] [Google Scholar]
  48. Kaper, L., van der Meer, A., & Tijani, A.H. 2004, The Environment and Evolution of Double and Multiple Stars, Proc. IAU Colloq. 191, 21, 128 [NASA ADS] [Google Scholar]
  49. Kaper, L., van der Meer, A., & Najarro, F. 2006, A&A, 457, 595 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Kobulnicky, H. A., Chick, W. T., Schurhammer, D. P. et al. 2016, ApJS, 2, 18 [NASA ADS] [CrossRef] [Google Scholar]
  51. Levenhagen, R. S., & Leister, N. V. 2006, MNRAS, 371, 252 [NASA ADS] [CrossRef] [Google Scholar]
  52. Lewin, W.H.G., & van der Klis, M. 2006, Compact Stellar X-ray Sources (Cambridge, UK: Cambridge University Press) [CrossRef] [Google Scholar]
  53. Lewin, W.H.G., van Paradijs, J., van den Heuvel, E.P.J., et al. 1997, X-ray Binaries (Cambridge, UK: Cambridge University Press) [Google Scholar]
  54. Lindegren, L., Lammers, U., Bastian, U., et al. 2016, A&A, 595, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  55. Lindegren, L., Hernandez, J., Bombrun, A., et al. 2018, A&A, 616, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. Liu, Q.Z., van Paradijs, J., & van den Heuvel, E.P.J. 2006, A&A, 455, 1165 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Luri, X., Brown, A. G. A., Sarro, L. M., et al. 2018, A&A, 616, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  58. Makovoz, D., & Khan, I. 2005, in Astronomical Data Analysis Software and Systems XIV, eds. P. Shopbell, M. Britton, & R. Ebert, ASP Conf. Ser., 347, 81 [Google Scholar]
  59. Manousakis, A., & Walter, R. 2011, A&A, 526, A62 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  60. Masetti, N., Parisi, P., Palazzi, E., et al. 2009, A&A, 495, 121 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Mason, K. O., Murdin, P. G., Parkes, G. E., et al. 1978, MNRAS, 184, 45P [NASA ADS] [CrossRef] [Google Scholar]
  62. Mason, A. B., Clark, J. S., Norton, A. J., et al. 2009, A&A, 505, 281 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Mason, A. B., Norton, A. J., Clark, J. S., et al. 2010, A&A, 509, A79 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  64. Mayer, A., Deschamps, R., & Jorissen, A. 2016, A&A, 587, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  65. McClintock, J.E., Ricker, G.R., & Lewin, W.H.G. 1971, ApJ, 166, L73 [NASA ADS] [CrossRef] [Google Scholar]
  66. Megier, A., Strobel, A., Galazutdinov, G. A., et al. 2009, A&A, 507, 833 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  67. Meyer, D. M. -A., van Marle, A. -J., Kuiper, R. et al. 2016, MNRAS, 459, 1146 [NASA ADS] [CrossRef] [Google Scholar]
  68. Meyer, D. M. -A., Mignone, A., Kuiper, R. et al. 2017, MNRAS, 464, 3229 [NASA ADS] [CrossRef] [Google Scholar]
  69. Negueruela, I., Smith, D. M., Harrison, T. E., et al. 2006a, ApJ, 638, 982 [NASA ADS] [CrossRef] [Google Scholar]
  70. Negueruela, I., Smith, D. M., Reig, P., et al. 2006b, ESA SP, 604, 165 [NASA ADS] [Google Scholar]
  71. Negueruela, I., Torrejón, J. M., & McBride, V. 2007, ATel, 1239 [Google Scholar]
  72. Negueruela, I., Casares, J., Verrecchia, F., et al. 2008, ATel, 1876 [Google Scholar]
  73. Nespoli, E., Fabregat, J., & Mennickent, R. E. 2008, A&A, 486, 911 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  74. Ochsendorf, B. B., Cox, N. L. J., Krijt, S., et al. 2014a, A&A, 563, A65 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  75. Ochsendorf, B. B., Verdolini, S., Cox, N. L. J., et al. 2014b, A&A, 566, A75 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  76. Parker, Q. A., Phillipps, S., Pierce, M. J., et al. 2005, MNRAS, 362, 689 [NASA ADS] [CrossRef] [Google Scholar]
  77. Parkes, G.E., Murdin, P.G., & Mason, K.O. 1980, MNRAS, 190, 537 [NASA ADS] [CrossRef] [Google Scholar]
  78. Peri, C.S., Benaglia, P., Brookes, D.P., et al. 2012, A&A, 538, A108 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  79. Peri, C. S., Benaglia, P., & Isequilla, N. L. 2015, A&A, 578, A45 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  80. Perryman, M. A. C., Lindegren, L., Kovalevsky, J., et al. 1997, A&A, 323, L49 [NASA ADS] [Google Scholar]
  81. Pflamm-Altenburg, J., & Kroupa, P. 2010, MNRAS, 404, 1564 [NASA ADS] [Google Scholar]
  82. Pilbratt, G.L., Riedinger, J.R., Passvogel, T., et al. 2010, A&A, 518, L1 [CrossRef] [EDP Sciences] [Google Scholar]
  83. Poglitsch, A., Waelkens, C., Geis, N., et al. 2010, A&A, 518, L2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  84. Polcaro, V. F., Bazzano, A., La Padula, C., et al. 1984, A&A, 131, 229 [NASA ADS] [Google Scholar]
  85. Postnov, K., Oskinova, L., & Torrejón, J. M. 2017, MNRAS, 465, L119 [NASA ADS] [CrossRef] [Google Scholar]
  86. Povich, M. S., Benjamin, R. A., Whitney, B. A., et al. 2008, ApJ, 689, 242 [NASA ADS] [CrossRef] [Google Scholar]
  87. Priedhorsky, W.C., & Terrell, J. 1983, ApJ, 273, 709 [NASA ADS] [CrossRef] [Google Scholar]
  88. Rahoui, F., Chaty, S., Lagage, P. O., et al. 2008, A&A, 484, 801 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  89. Reynolds, A. P., Owens, A., Kaper, L., et al. 1999, A&A, 349, 873 [NASA ADS] [Google Scholar]
  90. Rieke, G.H., Young, E.T., Engelbracht, C.W., et al. 2004, ApJS, 154, 25 [NASA ADS] [CrossRef] [Google Scholar]
  91. Romano, P., Mangano, V., Cusumano, G., et al. 2011, MNRAS, 412, L30 [NASA ADS] [Google Scholar]
  92. Sato, N., Nagase, F., Kawai, N., et al. 1986, ApJ, 304, 241 [NASA ADS] [CrossRef] [Google Scholar]
  93. Secchi, A. 1866, Astron. Nachr., 68, 6 [Google Scholar]
  94. Servillat, M., Coleiro, A., Chaty, S., et al. 2014, ApJ, 797, 114 [NASA ADS] [CrossRef] [Google Scholar]
  95. Smith, D. M., Main, D., Marshall, F., et al. 1998, ApJ, 501, L181 [NASA ADS] [CrossRef] [Google Scholar]
  96. Smith, M. A., Lopes de Oliveira, R., & Motch, C. 2016, Adv. Space Res., 58, 782 [NASA ADS] [CrossRef] [Google Scholar]
  97. Smith, M. A., Lopes de Oliveira, R., & Motch, C. 2017, MNRAS, 469, 1502 [NASA ADS] [CrossRef] [Google Scholar]
  98. Sota, A., Maíz Apellániz, J., Morrell, N. I., et al. 2014, ApJS, 211, 10 [NASA ADS] [CrossRef] [Google Scholar]
  99. Sugizaki, M., Mitsuda, K., Kaneda, H., et al. 2001, ApJS, 134, 77 [NASA ADS] [CrossRef] [Google Scholar]
  100. Takeuchi, Y., Koyama, K., & Warwick, R. S. 1990, PASJ, 42, 287 [NASA ADS] [Google Scholar]
  101. Tawara, Y., Yamauchi, S., Awaki, H., et al. 1989, PASJ, 41, 473 [NASA ADS] [Google Scholar]
  102. Thompson, T.W.J., Tomsick, J.A., in ’t Zand, J.J.M., et al. 2007, ApJ, 661, 447 [NASA ADS] [CrossRef] [Google Scholar]
  103. Toalá, J. A., Oskinova, L. M., & Ignace, R. 2017, ApJ, 838, L19 [NASA ADS] [CrossRef] [Google Scholar]
  104. Tomsick, J. A., Chaty, S., Rodriguez, J., et al. 2007, ATel, 1231 [Google Scholar]
  105. van Buren, D., & McCray, R. 1988, ApJ, 329, L93 [NASA ADS] [CrossRef] [Google Scholar]
  106. van Buren, D., Noriega-Crespo, A., & Dgani, R. 1995, AJ, 110, 2914 [NASA ADS] [CrossRef] [Google Scholar]
  107. van den Heuvel, E.P.J., Portegies Zwart, S.F., Bhattacharya, D., et al. 2000, A&A, 364, 563 [NASA ADS] [Google Scholar]
  108. van Oijen, J.G.J. 1989, A&A, 217, 115 [NASA ADS] [Google Scholar]
  109. Verrecchia, F., Negueruela, I., Covino, S., et al. 2002, ATel, 102 [Google Scholar]
  110. Walter, R., Lutovinov, A.A., Bozzo, E., et al. 2015, A&ARv, 23, 2 [NASA ADS] [CrossRef] [Google Scholar]
  111. Warwick, R.S., Norton, A.J., Turner, M.J.L., et al. 1988, MNRAS, 232, 551 [NASA ADS] [CrossRef] [Google Scholar]
  112. Werner, M., Roellig, T., Low, F., et al. 2004, ApJS, 154, 1 [NASA ADS] [CrossRef] [Google Scholar]
  113. Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868 [NASA ADS] [CrossRef] [Google Scholar]
  114. Zacharias, N., Finch, C., & Frouard, J. 2017, AJ, 153, 166 [NASA ADS] [CrossRef] [Google Scholar]
  115. Ziolkowski, J. 2002, Mem. Soc. Astron. It., 73, 1038 [NASA ADS] [Google Scholar]
  116. Zurita Heras, J.A., De Cesare, G., Walter, R., et al. 2006, A&A, 448, 261 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

All Tables

Table 1

Proper motions and derived peculiar tangential velocities for selected HMXBs.

In the text

All Figures

thumbnail Fig. 1

WISE 22 μm image of the field around γ Cas with 2MASS Ks contours (sensitive to starlight) overlaid on top. The system position is indicated by the cyan circle.
In the text
thumbnail Fig. 2

Spitzer MIPS 24 μm image of the field around EXO 051910+3737.7 (V420 Aur) with Spitzer IRAC 3.6 μ m contours (mainly sensitive to starlight) overlaid on top. The system position is indicated by the cyan circle.
In the text
thumbnail Fig. 3

WISE 22 μm image of the field around IGR J11435-6109 with 2MASS K band contours (sensitive to starlight) overlaid on top. The system position is indicated by the cyan circle.
In the text
thumbnail Fig. 4

Spitzer MIPS 24 μm image of the field around GX 301-2 (BP Cru) with Spitzer IRAC 3.6 μm contours (mainly sensitive to starlight) overlaid on top. The system position is indicated by the cyan circle.
In the text
thumbnail Fig. 5

WISE 22 μm image of the field around 1H 1255-567 (μ2 Cru) with 2MASS Ks contours (sensitive to starlight) overlaid on top. The system position is indicated by the cyan circle.
In the text
thumbnail Fig. 6

WISE 22 μm image of the field around GX 304-01 (V850 Cen) with 2MASS Ks contours (sensitive to starlight) overlaid on top. The system position is denoted by the cyan circle.
In the text
thumbnail Fig. 7

Spitzer MIPS 24 μm image of the field around AX J1639.0-4642 with Spitzer IRAC 3.6 μm contours (mainly sensitive to starlight) overlaid on top. The system position is denoted by the cyan circle.
In the text
thumbnail Fig. 8

WISE 22 μm image of the field around 4U 1700-37 (V884 Sco) with WISE 3.4 μm contours (mainly sensitive to starlight) overlaid on top. The system position is denoted by the cyan circle.
In the text
thumbnail Fig. 9

Spitzer MIPS 24 μm image of the field around EXO 1722-363 with Spitzer IRAC 3.6 μm contours (mainly sensitive to starlight) overlaid on top. The system position is indicated by the cyan circle.
In the text
thumbnail Fig. 10

Spitzer MIPS 24 μm image of the field around XTE J1739-302 with Spitzer IRAC 3.6 μm contours (mainly sensitive to starlight) overlaid on top. The system position is indicated by the cyan circle.
In the text
thumbnail Fig. 11

Spitzer MIPS 24 μm image of the field around AX J1841.0-0536 with Spitzer IRAC 3.6 μm contours (mainly sensitive to starlight) overlaid on top. The system position is indicated by the cyan circle.
In the text
thumbnail Fig. 12

WISE 22 μm image of the field around XTE J1855-026 with WISE 3.4 μm contours (mainly sensitive to starlight) overlaid on top. The system position is indicated by the cyan circle.
In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.