Free Access
Volume 543, July 2012
Article Number A56
Number of page(s) 10
Section Interstellar and circumstellar matter
Published online 27 June 2012

© ESO, 2012

1. Introduction

Runaway stars have high peculiar velocities, v   ~ 30 km s-1 (e.g. Gies & Bolton 1986; Tetzlaff et al. 2011). These stars move supersonically through the interstellar medium (ISM) and are ejected from their birth associations by gravitational kicks (e.g. Fujii & Zwartz 2011). Runaway stars with strong winds – spectral types OB and Wolf-Rayet (WR) – can produce stellar bowshocks when interacting with the surrounding gas (e.g. Van Buren et al. 1995). Bowshocks develop as arc-shaped structures, with bows pointing ahead of the stars (in the same direction as the stellar velocity).

There are several cataloged objects of this type (e.g. Van Buren et al. 1995; Noriega-Crespo et al. 1997; Kobulnicky et al. 2009; Gvaramadze et al. 2011; Peri et al. 2011). Not all runaway stars produce a bowshock. Van Buren et al. (1995) concluded that only 30% of the cataloged runaway stars present bowshock-like structures. More recently, Peri et al. (2011) found out that only  ~10% of early-type runaway stars develop an observable bowshock in the infrared (IR). Some physical situations exist in which a bowshock does not form. If the star is moving in a low-density, hot plasma then it might be moving subsonically and a shock will not occur. Also, a bowshock might not form if the star has a weak wind or if it has a high space velocity (see Comerón & Kaper 1998; Huthoff & Kaper 2002).

Bowshocks occur around many classes of astrophysical sources: pulsars, cataclysmic variables, colliding wind binaries, cometary H II regions, and even in groups and clusters of galaxies. A lot of work has been done on bowshock modeling (e.g. Van Buren & McCray 1988; Van Buren et al. 1990; Bandiera 1993; Van Buren 1993; Brighenti & Dércole 1995; Wilkin 1996; Comerón 1997; Chen & Huang 1997; Comerón & Kaper 1998; Wilkin 2000; Wareing et al. 2007). Related physical problems such as solving the astrophysical blunt-body problem are treated in e.g. Cantó & Raga (1998) and more recently in Schulreich & Breitschwerdt (2011); the physics of mixing layers is discussed in e.g. Baranov et al. (1976), Cantó & Raga (1991) and Raga et al. (1995).

The supersonic stellar wind sweeps the ISM material and piles it up in the bowshock. The stellar and shock-excited radiation heats this swept-up material. The dust, in turn, re-radiates the energy as mid-to-far IR flux (Van Buren & McCray 1988).

Relativistic particles can be accelerated at strong shocks producing non-thermal emission (Drury 1983). Benaglia et al. (2010) have reported non-thermal radio emission from the bowshock of the runaway star BD +43°3654. This emission is thought to be synchrotron radiation generated by the interaction of relativistic electrons with the magnetic field of the source.

In this work we present a model for the radiative non-thermal emission that takes place in the bowshocks of runaway stars. We compute the emission produced by the relativistic particles accelerated at the shock, and the thermal emission caused by free-free mechanism (Bremsstrahlung). Finally, we apply this model to the well-known stellar bowshock from the O9.5V star ζ Oph. The bowshock shape is also computed, following the analytical method developed by Wilkin (2000).

Our paper is organized as follows. In the next section we present the radiative model developed to compute the non-thermal emission and the absorption that can take place in these types of sources. We also present a brief discussion on particle acceleration in shocks. In Sect. 3 we apply the model to a particular source, the bowshock from ζ Oph. We present the best-fit shell shape and the computed spectral energy distribution (SED). Finally, in Sect. 4 we discuss the detectability of the non-thermal radiation and offer our conclusions.

2. Radiative model

Bowshocks from runaway stars produce thermal emission by radiative heating of the swept-up dust by the stellar radiation field. The infrared signal of the heated dusty bowshock is strongest in the far infrared (e.g. Van Buren & McCray 1988; Kobulnicky et al. 2009). The dust temperature can be understood in terms of dust models. For a typical bowshock of a runaway star, a fraction  ~ 10-2 of the star bolometric luminosity is emitted in the infrared. In this paper we do not calculate this radiation, and leave it for a future work (del Valle et al., in prep.).

The shocked ISM can also produce thermal emission through free-free interactions (Bremsstrahlung). This emission peaks at energies  ~ 1 eV; we calculate it in Sect. 2.5.2.

The interactions of locally accelerated relativistic particles with the matter, radiation and magnetic fields in the system produce non-thermal radiation. We calculate this emission below.

2.1. Shocks and particle acceleration

Shocks transfer kinetic energy to non-thermal particles through particle acceleration. The acceleration mechanism is diffusive shock acceleration, the so-called Fermi I mechanism (e.g. Bell 1978). In this mechanism particles are accelerated by successive bouncing across the shock, gaining energy in each cross. The particle deflection is mediated by magnetic field irregularities. These irregularities are expected from turbulent and magnetic instabilities. For the mechanism to operate is necessary that the particles in the shocked medium (downstream) can effectively diffuse and reach the shock.

The particle energy gain in each cycle (upstream-downstream-upstream) is ΔE/E  ∝  (vs/c), where vs is the shock velocity; and after k cycles the particle energy is , where Ei is the initial energy. The time that a particle requires to accelerate up to an energy E is given by (1)Here, B is the magnetic field in the acceleration region, and η is the acceleration efficiency (Drury 1983): (2)D is the diffusion coefficient, and rg = E/(eB) is the particle gyroradius. In the Bohm limit DB = rgc/3. Fast shocks are more efficient accelerators than slow shocks. The resulting spectrum of injected particles by this mechanism is a power-law, i.e. Q(E) ∝ E − α (e.g. Protheroe 1999), with α ~ 2 − 2.2.

thumbnail Fig. 1

Scheme of the system of shocks (not to scale). The acceleration region is indicated by dots.

The collision of the supersonic stellar wind with the ISM around a runaway star results in a system of two shocks (see Fig. 1). In the steady state, mass and momentum are conserved and a flow is established between the two shocks that carries away the mass and momentum deposited by the colliding fluids. In the reference frame of the star, the ISM can be considered as an incoming wind of parallel streamlines and the stellar wind as a radial outflow (e.g. Wilkin 1996). One of the shocks is a forward shock that propagates in the same direction as the stellar wind with velocity vfs  ~  V  ; and the other is a reverse shock that propagates in the opposite direction with vrs  ~  VW. The stellar wind can be considered as a continuous power source, therefore both shocks reach a steady state. If radiative losses are inefficient within the shock discontinuity, the shock is adiabatic; otherwise, the shock is radiative.

The stellar wind is divergent because its ram pressure decreases with distance, while the ram pressure of the ISM is constant. The point where the ram pressure of the wind and the ISM balance, i.e. , is called the standoff point, where ρw = w/4πR2Vw. This point of the shell defines the standoff radius R0: (3)The physical conditions in bowshocks from runaway stars produce an adiabatic and fast reverse shock at the shocked stellar wind; and the forward shock is radiative and slow (e.g. Van Buren 1993). The radiative shock compresses the material, and as the temperature decreases, the density grows.

Kis et al. (2004) demonstrated that upstream particles in the Earth bowshock undergo diffusive transport into the upstream region, a direct evidence of Fermi I acceleration (e.g. Burgess 2007). Benaglia et al. (2010) reported the detection of non-thermal radio emission from the bowshock of the runaway star BD+43°3654. The non-thermal radiation is expected from synchrotron emission generated by relativistic electrons accelerated at either the forward shock or the reverse shock.

We considered that an initial supra-thermal population of relativistic particles, electrons and protons, are accelerated in the adiabatic shock (reverse shock). This shock is faster than the forward shock, so in the reverse shock the relativistic particles are accelerated more efficiently by the Fermi mechanism (see Eq. (2)). All equations given above are valid for planar shocks. We considered an acceleration region where the shock is nearly flat (see Fig. 1). The width of the shocked stellar wind, Δ, can be estimated as Δ ~ M-2R0, where M is the Mach number of the shocked wind. The acceleration region is assumed as a small region near the bowshock apex, of scale length  ~ Δ; the acceleration volume is given in Table 1.

The available power in the system is the kinetic power from the stellar wind: (4)To estimate the magnetic field in the flow, we consider that the magnetic energy density is in subequipartition with respect to the kinetic energy LT, by a 0.1 factor, i.e.: (5)where A is the area of a sphere of radius R0. This condition ensures that the flow is matter-dominated, i.e. compressible, allowing shocks to develop.

The kinetic power available in the acceleration region is L = fLT, where f is the ratio of the volume of a sphere of radius R0 and the volume of the acceleration region. A small fraction of this kinetic power goes into relativistic particles. Lrel = qrelL. We adopted a standard fraction qrel = 0.1 (e.g. Protheroe 1999). We took into account both hadronic and leptonic content in the relativistic power, Lrel = Lp + aLe. The ratio of relativistic protons to electrons, a, is unknown. We considered two cases a = 1 (equal energy density in both species) and a = 100 (as observed in Galactic cosmic rays, Ginzburg & Syrovatskii 1964).

Here we discuss two different types of massive star: an O4I and an O9I star. The values for the parameters involved are given in Table 1.

Table 1

Parameters for the different types of stars.

2.2. Non-thermal radiative losses

The electrons lose energy mainly by inverse Compton (IC) scattering, synchrotron radiation, and relativistic Bremsstrahlung.

The synchrotron cooling time rate is (6)where σT is the Thomson cross section and UB is the magnetic energy density.

The IC losses can be calculated from (Blumenthal & Gould 1970) (7)where ϵ and ϵ1 are the incident photon and scattered photon energy, respectively, and (8)with (9)where b = 4ϵγ/mc2 and q = ϵ1/ [b(Ee − ϵ1)] . This expression takes into account losses in the Klein-Nishina regimen. The target photon fields are the radiation fields of the acceleration region. The fields considered here are the stellar radiation field at a distance R0 from the star – assumed as a black body at temperature T  –, and the IR radiation from the heated dust, also considered as a black body at TIR. Dust grains are heated by starlight, and cool by radiating in the infrared. To estimate TIR we used a simplified dust model based on Draine & Lee (1984): (10)(e.g. Van Buren & McCray 1988). Here aμm ~ 0.2 μm is the dust grain radius. More detailed and complex dust emission models are beyond the scope of this work. For more details see e.g. Draine & Li (2007) and Draine (2011).

The relativistic Bremsstrahlung losses are calculated considering a complete ionized plasma, using (Berezinskii et al. 1990) (11)where n is the density of target ions in the acceleration region (the shocked stellar wind in this case). The density of the shocked wind according to the Rankine-Hugoniot equations for adiabatic shocks (e.g. Landau & Lifshitz 1959) is 4nw, where nw is the wind density.

Protons lose energy through proton-proton inelastic collisions with the shocked wind material. These interactions produce neutral and charged pions; the former decay and produced gamma rays (e.g. Vila & Aharonian 2010). The latter decay into secondary electrons/positrons and neutrinos. The losses produced by proton-proton interactions occur on time scales of (12)where n is the density of target protons,  ~ 4nw, and Kpp is the ineslasticity (~0.5). The cross section can be approximated by (Kelner et al. 2006) (13)where L = ln(Ep/1TeV); and Eth = 280 MeV.

thumbnail Fig. 2

Acceleration and cooling time scales for electrons and protons for an O4I star (up), and for an O9I star (down). Left panels are for electrons and right panels are for protons. The non-radiative losess produced by convection are also shown.

Photomeson production is irrelevant at the energies considered in this paper.

In Fig. 2 we show the cooling rates for electrons and protons in the acceleration region for an O4I and an O9I star. In the O4I star system the IC scattering of IR photons dominates the radiative losses. For the O9I star, the IC scattering of IR photons and synchrotron radiation prevail among the radiative losses.

The relativistic particles can also suffer from non-radiative losses due to escape from the acceleration region. Particles can be convected away by the stellar wind on a time tconv ~ Δ/Vw. This non-radiative loss dominates the proton cooling rate in both types of runaway stars (O4I and O9I).

The minimum kinetic energy for each particle is considered to be on the order of the rest mass energy. The maximum energy for the electrons and protons is obtained equating the smallest cooling rate to the acceleration rate, given by Eq. (1). In the O4I system the electrons reach energies  ~10 TeV and the protons can reach energies  ~102 TeV. For O9I electrons reach energies  ~TeV, and protons energies  ~10 TeV. In both cases the Hillas criterion is satisfied: (14)where rg is the maximum particle gyroradius in the available space, i.e.  ~Δ.

2.3. Particle distributions

To calculate the steady state particle distributions N(E) for electrons and protons, we solved the transport equation in steady state (Ginzburg & Syrovatskii 1964): (15)where tesc is the convection time. We assumed that all physical properties in the acceleration region are homogeneous.

Q(E) is the injection function, a power-law in the energy of the particles, as expected from diffusive shock acceleration: (16)The normalization constant Q0 for each type of particles is obtained from Le,p as (17)Here V is the volume of the acceleration region.

The exact solution to the equation is also a power-law: (18)with (19)

thumbnail Fig. 3

Non-thermal emission from the shocked stellar wind for an O4I star (left), and for an O9I star (right).

2.4. Non-thermal radiative processes

As mentioned before, the interactions of the relativistic particles with the various fields in the system produce non-thermal emission.

We considered synchrotron emission, inverse Compton scattering with the stellar photon field and the IR emission, and relativistic Bremsstrahlung from electrons. For protons, inelastic collisions between the relativistic protons with the shocked wind material were also calculated.

The synchrotron emission is given by (20)where (21)κSSA is the synchrotron self-absorption (SSA) factor: (22)here τSSA(Eγ) is the SSA optical depth (see Rybicki & Lightman 1979).

The IC emission is calculated using (23)The spectrum of scattered photons is (24)with f(q) given by Eq. (9).

The relativistic Bremsstrahlung luminosity is (25)where (26)and (27)To compute the gamma-ray emission produced by neutral pion decays for Ep < 0.1 TeV, we used the following expression: (28)with Emin = Eγ + mπc4/4Eπ. The π0-emissivity is given by (Aharonian & Atoyan 2000) (29)with κπ ~ 0.17 (Gaisser 1990). For Ep < 0.1 TeV down to the threshold energy, the following replacement is necessary (30)Here is the total number of π0 created per p − p collision, and Ekin = Ep − mpc2 is the proton kinetic energy.

The gamma-ray luminosity in the range 0.1 TeV  ≤ Ep ≤ 105 TeV can be obtained from (Kelner et al. 2006): (31)with a function of Eγ and Ep (see Vila & Aharonian 2009, and references therein).

Secondary electron-positron pairs are created through the p − p inelastic collisions (see Orellana et al. 2007). These particles can also produce non-thermal emission, which we took into account. For the relevant formulae see Orellana et al.’s paper.

Figure 3 shows the computed SED for both types of stars. We calculated the contribution of the secondary pairs, to synchrotron emission and IC luminosity, only for the O9I star (because we assumed a proton dominated content in this case). The non-thermal luminosity of the O4I star is higher, mainly because the available power in the system is greater. The IC emission produced in the interaction with the dust photons dominates both SEDs at high energies. At lower energies synchrotron radiation dominates, but is relatively more important in the O9I system.

The IC emission from the dust photons, in general, increases as the temperature of the dust also increases. This temperature depends, as a first aproximation, on L   and R0 (see Eq. (10)). As argued in Sect. 2.1, B must be at subequipartition values. This parameter not only affects the synchrotron emission, but the acceleration rate (see Eq. (1)). A higher value of B gives a lower acceleration rate and the maximum particle energies change (decrease).

2.5. Emission of the shocked ISM

2.5.1. Non-thermal emission

Some relativistic particles can diffuse from the acceleration region to the shocked ISM before losing most of their energy. The diffusion time can be estimated as Δ2/2DB, where DB is the diffusion coefficient in the Bohm limit as before. Particles that diffuse to the denser shocked ISM (see Fig. 1) can interact there through p − p inelastic collisions and relativistic Bremsstrahlung. The shocked ISM is contained within a very thin layer.

Because of the radiative losses the temperature decreases and the density increases in the shocked ISM; both quantities are described by a temperature profile Tfs(x) and a density profile nfs(x). The temperature profile is given by (Zhekov & Palla 2007)

(32)where x is the distance from the shock; K and , are the corresponding values of temperature and pressure for the adiabatic case.

The density profile is (33)with C0 a normalization constant such that na ≡ nfs(Ta), where Ta is the ambient temperature. We estimated the width of the shocked region as xc, where T(xc) = Ta.

The particle energy distribution of the diffused protons is obtained solving Eq. (15), considering a 3D delta function injection (see Aharonian & Atoyan 1996; Bosch-Ramon et al. 2004): (34)where Np is the initial energy distribution of the particles accelerated at the reverse shock. The general solution is (35)where τpp  ~  yr. Here we consider an average value for the shocked density ; t is the propagation time after the injection into the shocked ISM, t ~ cxc. Rdiff is the so-called diffusion radius. For t ≪ τpp, , with D the diffusion coefficient, we considered D = DB, see Sect. 2.1.

The particle energy distributions of the diffused electrons can be approximated by (36)where . The interactions of diffused particles with the matter of the shocked ISM – through p − p interactions for protons and relativistic Bremsstrahlung for electrons – produce non-thermal emission.

For the two stars considered here, only the protons from the O9I system can diffuse into the shocked ISM. Figure 4 shows this contribution to the non-thermal SED. This emission peaks at E ~ 100 MeV, with L ~ 1028 erg s-1. From Fig. 3 it can be seen that this contribution is negligible compared to the total IC radiation.

2.5.2. Thermal emission

thumbnail Fig. 4

γ-ray emission from the shocked ISM produced by p − p interactions from an O9I star.

thumbnail Fig. 5

Thermal emission from the shocked ISM produced by free-free mechanism for an O4I star and for an O9I star.

The shocked ISM produces thermal emission through free-free mechanism (thermal Bremsstrahlung). For completeness we computed this contribution to the total SED.

We calculated the thermal emission by integrating the emissivity ϵ along the shocked ISM, considering the temperature and density profiles given by Eqs. (32) and (33). The emissivity is given by e.g. (Lang 1999) (37)where g(ν,T) is the free-free Gaunt factor given by  ~0.54ln [5 × 107(T3/2/ν)] ; ni and ne are the ion and electron densities respectively, we assumed ni = ne = nfs. The emission volume is a spherical wedge of width xc and radius Δ, see Fig. 1. In Fig. 5 the result of this contribution is shown.

For the O4I star the free-free emission is negligible (see Fig. 3). For the O9I case, at energies around  ~ 1 eV this thermal contribution dominates over synchrotron radiation.

2.6. Absorption

Gamma rays can be absorbed in the acceleration region immediately after they are created, by photon-photon annihilation. All radiation fields in the acceleration region, thermal and non-thermal, provide target photons for pair creation.

The differential opacity for a gamma ray traveling in the direction eγ due to photons of an energy ϵ in the direction eph is


where dΩ is the solid angle of the emitting surface and nϵ is the radiation density. The photon annihilation cross-section is (Gould & Schréder 1967) (39)where β = (1 − 1/s)1/2, and s = Eγϵ(1 − eγeph)/(mec2)2. Here, Eγ and ϵ are the energies of the gamma ray and the target photon, respectively. The threshold energy is given by (40)The opacity is (41)Here, u = cosϑ, ϑ is the angle between the momenta of the colliding photons, and l is the photon path across the target radiation field (here l = Δ). The target photon fields are those generated within the acceleration region, the IR emission from the heated dust, and the star radiation field at a distance d = R0 − Δ.

The optical depth is a trajectory integral on which the angular dependence has a very significant effect, and consequently the absorption depends strongly on the line of sight (e.g. Romero et al. 2010). Photons traveling from the acceleration region toward an observer placed at A, see Fig. 6, are not absorbed by the stellar photon field. On the other hand, photons traveling toward an observer placed at B strongly interact with the stellar photon field. This interaction depends on the distance d between the star and the photon path, and therefore depends on the angle α.

Photons of lower energies are absorbed by matter through photoionization. The target material can be the shocked ISM and the material along the line of sight corresponding to each particular source. The optical depth τγH can be approximated as (42)Here NH is the column density corresponding to the target density. The cross section σγN was taken from Ryter (1996). For more details see Reynoso et al. (2011) and references therein.

Below we apply the radiative model described so far to the specific case of ζ Oph.

thumbnail Fig. 6

Diagram of the path followed by a gamma ray traveling from the acceleration region toward observers placed in front and behind the runaway star (not to scale).

3. Application to ζ Oph

The star ζ Oph (HD 149757) is one of the brightest massive stars in the northern hemisphere and has been intensively studied. This star has spectral type O9.5V and is a well-known runaway star. It is rapidly rotating with almost break-up velocity, with vsin(i) ~ 400−500 km s-1 (Walker et al. 1979; Repolust et al. 2004).

ζ Oph bowshock has been observed by the InfraRed Astronomical Satellite IRAS (Van Buren & McCray 1988; Noriega-Crespo et al. 1997) and more recently by the Wide-field Infrared Survey Explorer WISE1. The WISE image shows a very clear, regular structure. This source is quite nearby, located at  ~222 pc from the Earth (Megier et al. 2009).

Table 2

Parameters for ζ Oph.

3.1. Bowshock shape

We adopted a mass loss rate for ζ Oph of  ~10-7 M yr-1 (Fullerton et al. 2006) and terminal wind velocity Vw = 1550 km s-1 (e.g. Hubrig et al. 2011). For the ISM density we adopted a value of na ~ 10 cm-3; since ζ Oph is embedded in a cirrus cloud region, the density and ambient temperature might be higher than average (e.g. Vidal et al. 2011). These parameters yield R0 ~ 0.3 pc, which agrees well with the upper limit measured by Peri et al. (2011). The list of the values of the main parameters adopted in our calculations (we follow Marcolino et al. 2009) are in Table 2.

To compute the bowshock shape we used the analytical method developed by Wilkin (2000). We assumed that the angular dependence of the wind momentum flux is equatorial (c2 =  −1 and λ  = 30°, see Sect. 4.3 of Wilkin 2000). Figure 7 shows the 3D computed bowshock shape.

To compare the 2D image observed by WISE we rotated the 3D coordinate system by three angles. We defined the coordinates xp ≡ E and zp ≡ N (see Fig. 8) to describe the WISE image. As usual, the origin is placed at the position of the star. From the WISE image of the bowshock of ζ Oph it is evident that the midpoint from the star to the bowshock (see Fig. 4 from Peri et al. 2011) is in the direction of the star velocity. Figure 8 shows the best fit of ζ Oph bowshock.

Our result shows that a simple model can provide a good representation of the observational data, assuming that the IR image is a good tracer of the actual hydrodynamic bowshock shape. The differences between the observed and the intrinsic shape of the bowshock depend on many factors involving radiative transfer, cooling time, dust characteristics, dust distribution, and so on.

thumbnail Fig. 7

ζ Oph computed bowshock.

thumbnail Fig. 8

Top: image of ζ Oph from WISE ( Dots pointed by arrows correspond to the projection to the 3D structure. Bottom: ζ Oph projected bowshock.

thumbnail Fig. 9

Electron and proton losses, acceleration rates, diffusion and convection time scales – defined in Sect. 2.5 – for ζ Oph.

thumbnail Fig. 10

Internal opacity.

3.2. Spectral energy distribution

In Fig. 9 we show the radiative losses, the acceleration rates and the diffusion and convection times (see Sect. 2.5) for electrons and protons. The maximum energy is  ~TeV for both species of particles. These values are in accordance with the Hillas criterion (see Sect. 2.1).

The internal photon-photon optical depth in the bowshock is shown in Fig. 10. It is negligible in the energy ranges of interest. The photoelectrical absorption is also negligible due to the small amount of material that photons cross on their way to the observer. The external absorption is also negligible given the relative positions of the bowshock, the star and the observer.

The expected non-thermal luminosity of the source is weak. However, since ζ Oph is very nearby, the bowshock might be detectable at gamma-ray and X-ray wavelengths through long exposures, under the assumptions we made. We remark that the sensitivity shown in Fig. 11 for MAGIC is for 50 h of exposure over the source, and for Fermi it is for one year of integration. For these types of sources an instrument like CTA gives the best chance of detection. CTA might become a unique tool to explore the high-energy radiation produced by runaway massive stars and the population of relativistic particles generated in them.

thumbnail Fig. 11

Computed SED for ζ Oph bowshock, at d ~ 222 pc. The sensitivity for CTA, Fermi, MAGIC, XMM-Newton and VLA. IRAS data are also shown.

4. Discussion and conclusions

Bowshocks of runaway massive stars are natural candidates for particle acceleration. The different types of massive stars have different energetics (depending on the wind parameters). The available power for particle acceleration also depends on the distance between the star and the acceleration region. The different available powers produce different non-thermal fluxes. Under the assumptions we made the asymmetries that might arise in runaway bowshocks do not produce a difference in the emitted spectrum. The obtained SEDs depend essentially on the particular assumptions made for the particle acceleration, the magnetic field, and the dust emission.

The emission might be detectable at several wavelengths, provided that the source is close enough and long exposure times are used – a good candidate is ζ Oph –. The synchrotron emission expected at radio wavelengths might be detectable, as in the case of BD +43°3654. The undetectability can establish constraints on parameters such as the magnetic field in the shocked wind. Stellar bowshocks might also be detectable at X-ray wavelengths, although no runaway bowshock has been observed at these energies so far. Finally, a system like ζ Oph might be detectable at γ-rays by the future ground-based detector CTA, as well as by the Fermi satellite. The energy range between 1 GeV − 1 TeV offers the best prospects for the study of runaway stars as non-thermal emitters. Our work suggests that bowshocks of runaway stars might constitute a new class of high-energy sources to be explored in the near future.

Figure 11 shows the SED for the emission from the bowshock of ζ Oph, along with the sensitivity of the gamma-ray detectors CTA (Cherenkov Telescope Array-forthcoming-), MAGIC and Fermi, the X-ray satellite XMM-Newton (theoretical upper limit from Hasinger et al. 2001), and VLA (upper limit from the NVSS survey – Condon, Cotton, Greisen et al. 1998 –, angular resolution is given by Peri et al. 2011). For completeness the IR IRAS data are also shown (Van Buren & MacCray 1988).


We are grateful to Paula Benaglia, Anabella Araudo, Florencia Vieyro and Cintia Peri for insightful discussions. This publication used 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 is supported by PIP 0078 (CONICET) and PICT 2007-00848, Préstamo BID (ANPCyT). G.E.R. received additional support from the Spanish Ministerio de Inovación y Tecnología under grant AYA 2010-21782-c03-01.


  1. Aharonian, F. A., & Atoyan, A. M. 1996, A&A, 309, 917 [NASA ADS] [Google Scholar]
  2. Aharonian, F. A., & Atoyan, A. M. 2000, A&A, 362, 937 [NASA ADS] [Google Scholar]
  3. Bandiera, R. 1993, A&A, 276, 648 [NASA ADS] [Google Scholar]
  4. Baranov, V. B., Krasnobaev, K. V., & Ruderman, M. S. 1976, Ap&SS, 41, 481 [NASA ADS] [CrossRef] [Google Scholar]
  5. Bell, A. R. 1978, MNRAS, 182, 147 [NASA ADS] [CrossRef] [Google Scholar]
  6. Benaglia, P., Romero, G. E., Martí, J., Peri, C. S., & Araudo, A. T. 2010, A&A, 517, L10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Berezinskii, V. S., Buland, V., Dogid, V. A., et al. 1990, Astrophysics of Cosmic Rays (Amsterdam: North-Holland) [Google Scholar]
  8. Blumenthal, G. R., & Gould, R. J. 1970, Rev. Mod. Phys., 42, 237 [NASA ADS] [CrossRef] [Google Scholar]
  9. Bosch-Ramon, V., Aharonian, F. A., & Paredes, J. M. 2005, A&A, 432, 609 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Brighenti, F., & D’Ercole, A. 1995, MNRAS, 277, 53 [NASA ADS] [Google Scholar]
  11. Burgess, D. 2007, Lect. Notes Phys. 725, 161 [Google Scholar]
  12. Cantó, J., & Raga, A. C. 1991, ApJ, 372, 646 [NASA ADS] [CrossRef] [Google Scholar]
  13. Cantó, J., & Raga, A. C. 1998, MNRAS, 297, 383 [NASA ADS] [CrossRef] [Google Scholar]
  14. Chen, Y., & Huang, J. H. 1997, ApJ, 479, L23 [NASA ADS] [CrossRef] [Google Scholar]
  15. Comerón, E. 1997, A&A, 326, 1195 [NASA ADS] [Google Scholar]
  16. Comerón, F., & Kaper, L. 1998, A&A, 338, 273 [NASA ADS] [Google Scholar]
  17. Condon, J. J., Cotton, W. D., Greisen, E. W., et al. 1998, AJ, 115, 1693 [NASA ADS] [CrossRef] [Google Scholar]
  18. Draine, B. T. 2011, Physics of the Interstellar and Intergalactic Medium (Woodstock: Princeton University Press) [Google Scholar]
  19. Draine, B. T., & Lee, H. M. 1984, ApJ, 285, 89 [Google Scholar]
  20. Draine, B. T., & Li, A. 2007, ApJ, 657, 810 [NASA ADS] [CrossRef] [Google Scholar]
  21. Drury, L. O’ C. 1983, Rep. Prog. Phys., 46, 973 [NASA ADS] [CrossRef] [Google Scholar]
  22. Fujii, M. S., & Zwart, S. P. 2011, Science, 334, 380 [Google Scholar]
  23. Fullerton, A. W., Massa, D. L., & Prinja, R. K. 2006, ApJ, 637, 1025 [NASA ADS] [CrossRef] [Google Scholar]
  24. Gaisser, T. K. 1990, Cosmic Rays and Particle Physics, ed. T. K. Gaisser [Google Scholar]
  25. Gies, D., & Bolton, C. T. 1986, ApJS, 61, 419 [NASA ADS] [CrossRef] [Google Scholar]
  26. Ginzburg, L. V., & Syrovatskii, S. I. 1964, The Origin of Cosmic Rays (Oxford: Pergamon Press) [Google Scholar]
  27. Gould, R. J., & Schréreder, G. P. 1967, Phys. Rev., 155, 1404 [NASA ADS] [CrossRef] [Google Scholar]
  28. Gvaramadze, V. V., Kniazev, A. Y., Kroupa, P., & Oh, S. 2011, A&A, 535, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Hasinger, G., Altieri, B., Arnaud, M., et al. 2011, A&A, 365, L45 [Google Scholar]
  30. Hubrig, S., Oskinova, L. M., & Schöller, M. 2011, Astron. Nachr., 332, 147 [NASA ADS] [CrossRef] [Google Scholar]
  31. Kelner, S. R., Aharonian, F. A., & Bugayov, V. V. 2006, Phys. Rev. D, 74, 034018 [NASA ADS] [CrossRef] [Google Scholar]
  32. Kis, A., Scholer, M., Klecker, B., et al. 2004, Geophys. Res. Lett. 31, 208001, 183 [Google Scholar]
  33. Kobulnicky, H. A., Gilbert, I. J., & Kiminki, D. C. 2010, ApJ, 710, 549 [NASA ADS] [CrossRef] [Google Scholar]
  34. Landau, L. D., & Lifshit, E. M. 1959, Course of Theoretical Physics (Oxford: Pergamon Press) [Google Scholar]
  35. Lang, K. R. 1999, Astrophysical Formulae (Berlin: Springer-Verlag) [Google Scholar]
  36. Marcolino, W. L. F., Bouret, J. C., Martins, F., et al. 2009, A&A, 498, 837 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Martins, F., Schaerer, D., & Hillier, D. J. 2005, A&A, 436, 1049 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Megier, A., Strobel, A., Galazutdinov, G. A., et al. 2009, A&A, 507 833 [Google Scholar]
  39. Noriega-Crespo, A., Van Buren, D., & Dgani, R. 1997, AJ, 113, 780 [NASA ADS] [CrossRef] [Google Scholar]
  40. Orellana, M., Bordas, P., Bosch-Ramon, V., Romero, G. E., & Paredes, J. M. 2007, A&A, 476, 9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Peri, C. S., Benaglia, P., Brookes, D. P., Stevens, I. R., & Isequilla, N. 2011, A&A, 538, A108 [Google Scholar]
  42. Protheroe, R. J. 1999, in Topics in Cosmic-Ray Astrophysics, Horizons in World Physics, ed. M. A. DuVernois (New York: Nova Science Publishers) [Google Scholar]
  43. Raga, A. C., Cabrit, S., & Cató 1995, MNRAS, 273, 422 [NASA ADS] [CrossRef] [Google Scholar]
  44. Repolust, T., Puls, J., & Herrero, A. 2004, A&A, 415, 349 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Reynoso, M. M., Medina, M. C., & Romero, G. E. 2011, A&A, 531, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Romero, G. E., & Vila, G. S. 2008, A&A, 485, 623 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Romero, G. E., del Valle, M. V., & Orellana, M. 2010, A&A, 518, A12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. Rybicki, G. B., & Lightman, A. P. 1979, Radiative Processes in Astrophysics (New York: John Wiley & Sons) [Google Scholar]
  49. Ryter, C. E. 1996, Ap&SS, 236, 285 [NASA ADS] [CrossRef] [Google Scholar]
  50. Schulreich, M. M., & Breitschwerdt, D. 2011, A&A, 531, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Tetzlaff, N., Neuhäuser, R., & Hohle, M. M. 2011, MNRAS, 410, 190 [NASA ADS] [CrossRef] [Google Scholar]
  52. Van Buren, D. 1993, ASP Conf. Ser., 35, 315 [NASA ADS] [Google Scholar]
  53. Van Buren, D., & McCray, R. 1988, ApJ, 329, L93 [NASA ADS] [CrossRef] [Google Scholar]
  54. Van Buren, D., Mac Low, M.-M., Wood, D. O. S., & Churchwell, E. 1990, ApJ, 353, 570 [NASA ADS] [CrossRef] [Google Scholar]
  55. Van Buren, D., Noriega-Crespo, A., & Dgani, R. 1995, AJ, 110, 2914 [NASA ADS] [CrossRef] [Google Scholar]
  56. Vidal, M., Casassus, S., Dickinson, C., et al. 2011, MNRAS, 414, 2424 [NASA ADS] [CrossRef] [Google Scholar]
  57. Vila, G. S., & Aharonian, F. 2009, in Compact Objects and their Emission, ed. G. E. Romero, & P. Benaglia (La Plata: Paideia) [Google Scholar]
  58. Walker, G. A. H., Yang, S., & Fahlman, G. G. 1979, ApJ, 233, 199 [NASA ADS] [CrossRef] [Google Scholar]
  59. Wareing, C. J., Zijlstra, A. A., & O’Brien, T. J. 2007, MNRAS, 382, 1233 [NASA ADS] [CrossRef] [Google Scholar]
  60. Wilkin, F. P. 1996, ApJ, 459, L31 [NASA ADS] [CrossRef] [Google Scholar]
  61. Wilkin, F. P. 2000, ApJ, 532, 400 [NASA ADS] [CrossRef] [Google Scholar]
  62. Zhekov, S. A., & Palla, F. 2007, MNRAS, 382, 1124 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1

Parameters for the different types of stars.

Table 2

Parameters for ζ Oph.

All Figures

thumbnail Fig. 1

Scheme of the system of shocks (not to scale). The acceleration region is indicated by dots.

In the text
thumbnail Fig. 2

Acceleration and cooling time scales for electrons and protons for an O4I star (up), and for an O9I star (down). Left panels are for electrons and right panels are for protons. The non-radiative losess produced by convection are also shown.

In the text
thumbnail Fig. 3

Non-thermal emission from the shocked stellar wind for an O4I star (left), and for an O9I star (right).

In the text
thumbnail Fig. 4

γ-ray emission from the shocked ISM produced by p − p interactions from an O9I star.

In the text
thumbnail Fig. 5

Thermal emission from the shocked ISM produced by free-free mechanism for an O4I star and for an O9I star.

In the text
thumbnail Fig. 6

Diagram of the path followed by a gamma ray traveling from the acceleration region toward observers placed in front and behind the runaway star (not to scale).

In the text
thumbnail Fig. 7

ζ Oph computed bowshock.

In the text
thumbnail Fig. 8

Top: image of ζ Oph from WISE ( Dots pointed by arrows correspond to the projection to the 3D structure. Bottom: ζ Oph projected bowshock.

In the text
thumbnail Fig. 9

Electron and proton losses, acceleration rates, diffusion and convection time scales – defined in Sect. 2.5 – for ζ Oph.

In the text
thumbnail Fig. 10

Internal opacity.

In the text
thumbnail Fig. 11

Computed SED for ζ Oph bowshock, at d ~ 222 pc. The sensitivity for CTA, Fermi, MAGIC, XMM-Newton and VLA. IRAS data are also shown.

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.