Clumpy stellar winds and high-energy emission in high-mass binaries hosting a young pulsar

# Clumpy stellar winds and high-energy emission in high-mass binaries hosting a young pulsar

V. Bosch-Ramon 1 Departament d’Astronomia i Meteorologia, Institut de Ciències del Cosmos (ICC), Universitat de Barcelona (IEEC-UB), Martí i Franquès 1, 08028 Barcelona, Spain, e-mail: vbosch@am.ub.es 1
Received ¡date¿ / Accepted ¡date¿
###### Key Words.:
Hydrodynamics – X-rays: binaries – Stars: winds, outflows – Radiation mechanisms: nonthermal – Gamma rays: stars
offprints: V. Bosch-Ramon,
###### Abstract

Context:High-mass binaries hosting young pulsars can be powerful gamma-ray emitters. The stellar wind of the massive star in the system, which interacts with the pulsar wind, is expected to be clumpy. Since the high-energy emission comes from the interaction of the two winds, the presence of clumps can affect the spectrum and variability of this radiation.

Aims:We look for the main effects of the presence of clumps in the stellar wind on the two-wind interaction region and on the non-thermal radiation that originates there.

Methods:A simple analytical model for the two-wind interaction dynamics was developed. The model accounts for the lifetime of clumps under the pulsar-wind impact. This time plays a very important role with regard to the evolution of the clump, the magnetic field in the clump-pulsar wind interaction region, and the non-radiative and radiative cooling of the non-thermal particles. We also computed the high-energy emission produced at the interaction of long-living clumps with the pulsar wind.

Results:For reasonable parameters, the clumps will induce small variability on the X-ray and gamma-ray radiation. Sporadically, large clumps can reach closer to the pulsar increasing the magnetic field, triggering synchrotron X-ray flares and weakening other emission components like inverse Compton. The reduction of the emitter size induced by clumps also makes non-radiative losses faster. Stellar wind clumps can also enhance instability development and matter entrainment in the shocked pulsar wind when it leaves the binary. Growth limitations of the clumps from the wind acceleration region may imply that a different origin for the largest clumps is required. The large-scale wind structures behind the observed discrete absorption components in the ultraviolet may be the source of these large clumps.

Conclusions:The presence of structure in the stellar wind can produce substantial energy-dependent variability and thus should not be neglected when studying the broadband emission from high-mass binaries hosting young pulsars.

## 1 Introduction

Several binary systems in the Galaxy are powerful gamma-ray emitters (see Paredes et al., 2013; Dubus, 2013, and references therein). Among these binaries, PSR B125963(/LS 2883) is a system formed by a late O star (Negueruela et al., 2011) and a powerful, relatively young, non-accreting pulsar (Johnston et al., 1992). The remaining gamma-ray emitting binaries, with the exception of the X-ray binaries Cygnus X-1 (to be confirmed) and Cygnus X-3, are objects of unclear nature that may also belong to the same class as PSR B125963.

Young-pulsar high-mass binaries are not expected to be very numerous in the Galaxy because of the short duration of the pulsar non-accreting phase compared with the whole life-span of massive X-ray binaries with neutron stars. In addition to PSR B125963, only a few less powerful systems lacking non-thermal radiation have been found so far (see Table 4 in Papitto et al. 2012). Assuming a lifetime for the non-accreting pulsar of a few times  yr (the age of PSR B125963), and the birth rates for high-mass binaries hosting a neutron star from Portegies Zwart & Yungelson (1999), one might expect of those systems in the Galaxy. Despite their relative scarcity, however, given the complex high-energy phenomena that arise from the conjunction of the orbital motion, the presence of a massive star, and the interaction of the stellar and the pulsar winds, these objects are important for the study of pulsar and stellar wind physics, particle acceleration, supersonic and relativistic fluid dynamics, emission mechanisms, and radiation reprocessing.

Important physical information can be extracted from young-pulsar high-mass binaries through observations in different energy bands and spectral, variability and morphological studies, although this requires a proper understanding of the main factors that come into play. In this regard, since the interaction region between the stellar and the pulsar winds is a natural site for non-thermal radiation (e.g. Maraschi & Treves, 1981; Tavani & Arons, 1997; Kirk et al., 1999; Dubus, 2006; Khangulyan et al., 2007; Neronov & Chernyakova, 2007, etc.), the stellar wind properties are to be properly characterized for a sensible modelling of the dynamics of the interacting flows and their emission.

Concerning the non-thermal emission from the two-wind interaction region, it is expected that the gamma-ray production will be mostly affected by dynamical processes taking place close to or within the binary (e.g. Khangulyan et al. 2007; see also Zabalza et al. 2013), whereas observations show (see Moldón et al., 2011) that radio emission is produced much farther from the system (see also Dubus, 2006; Bosch-Ramon, 2010). Therefore, the structure of the stellar wind and its impact on the non-thermal processes should be better probed at high energies at which the radiation from the two-wind interaction region can present the fastest and sharpest variations.

The structure of massive star winds is characterized by the presence of density and velocity irregularities or clumps (e.g. Runacres & Owocki, 2002). In the past, numerical and analytic studies of the influence of stellar wind clumping on thermal and non-thermal radiation have been performed for massive star binaries (e.g. Pittard, 2007), high-mass microquasars (e.g. Owocki et al., 2009; Araudo et al., 2009; Perucho & Bosch-Ramon, 2012), and LS I +61 303 (e.g. Zdziarski et al., 2010), a high-mass binary candidate that may host a non-accreting pulsar. So far, however, a specific investigation of the characteristic timescales and broadband flux variations introduced by the presence of stellar wind clumps in young pulsar high-mass binaries has not been done. In what follows, we present a simplified treatment of this problem to evaluate the observational impact of the arrival of clumps to the two-wind interaction region at X- and gamma rays. This study can provide a framework for more detailed calculations.

## 2 The clumpy stellar wind

The winds of massive stars are produced through radiation pressure exerted via the interaction of ultraviolet (UV) photons with electrons bound to atoms in the stellar atmosphere (Lucy & Solomon, 1970). Since the acceleration of the wind is a highly unstable process, it leads to strong irregularities in density and velocity (Runacres & Owocki, 2002). These irregularities, often referred to as clumps, have an impact on the observables of massive star winds, in particular affecting the inferred values of the mass-loss rates, which are crucial for the understanding of stellar evolution (e.g. Hirschi, 2008). The characteristics of clumps are still unclear, but it is thought that they carry most of the wind mass, but fill only a small fraction () of the volume in the stellar surroundings, yielding a contrast of in density with respect to its average value. The clump sizes are expected to be much smaller than the star (Owocki & Cohen, 2006), although larger clumps may also be present in the wind like the clumps used, for instance, to explain the X-ray flares in supergiant fast X-ray transients (e.g. Walter & Zurita Heras, 2007), or the large scale structures behind the observed discrete absorption components in the UV in early-type stars (DACs) (e.g. Cranmer & Owocki, 1996; Lobel & Blomme, 2008). The clumps should present a distribution of masses and sizes if they are the product of non-linear processes at the wind base (Moffat, 2008). In this investigation, clumps with different masses and sizes are empirically considered. The values discussed are representative of the majority of the clumps in the wind, i.e. the average clump case, or exceptional, i.e. the arrival of a rare and large matter condensation.

## 3 The two-wind interaction region

### 3.1 Physical framework

Stellar and pulsar winds collide forming two shocks and a region in between filled with shocked wind material. These two shocks contain a contact discontinuity that is located where the ram pressures of the two winds balance each other. The two-wind interaction structure becomes bow-shaped towards the wind with the lower momentum flux. In the top-left corner of Fig. 1, a representation of the colliding two-wind region is shown.

The shocked stellar and pulsar winds flow in contact but at very different speeds, suffering re-acceleration (Bogovalov et al., 2008), the Kelvin-Helmholtz instability, and mixing, as they leave the binary system (Lamberts et al., 2013). All these phenomena, and the orbital motion, play a major role in shaping the interaction structure at larger scales (Okazaki et al., 2011; Lamberts et al., 2012; Bosch-Ramon et al., 2012a), and possibly the highest energy emission as well (see Zabalza et al., 2013). In addition, a clumpy stellar wind can have consequences on the shocked flow dynamics and the radiation formation processes that might be significant and should not be overlooked.

Here, we restrict our study to the region confined within the binary system. In Fig. 1, a schematic of the considered scenario is presented. For simplicity, in what follows we consider the two-wind interaction region to be a circular slab with radius and thickness . The solid angle of the clump-pulsar wind interaction region, as seen from the pulsar, is . The quantities , , and are the distance between the two-wind interaction region and the pulsar, the clump radius, and the thickness of the shocked two-wind zone, respectively. We neglect here the bow shape of the two-wind and clump-pulsar wind shocks, as well as the shocked material flowing beyond the interaction slab. Under all these approximations, the escape timescale of the material in the shocked wind zones, similar to the local adiabatic cooling time, can be taken as and in the two-wind and the clump-pulsar wind shocks, respectively. The flow velocity will be taken as the sound speed, similar to the stellar wind (and clump) velocity , in the shocked stellar wind, and the relativistic sound speed in the shocked pulsar wind. In our study, the wind and its clumps always interact with the pulsar wind well beyond the pulsar magnetosphere.

The properties of the two-wind and the clump-pulsar wind interaction regions are assumed to be uniform, which will greatly simplify the radiation treatment in Sect. 4. Another simplification of the emitting regions done here is to consider them to be point-like (see Sect. 4.2). Finally, we assume that clumps arrive to the two-wind interaction region moving along the line joining the two stars. All these approximations will allow us to obtain order-of-magnitude estimates for the relevant timescales and flux variations at different energy bands, capturing the main features of the influence of the clump on the radiation. For accurate predictions, however, detailed studies of the dynamics and radiation should be combined, accounting also for the shocked flow evolution outside the binary system.

### 3.2 Relevant spatial and timescales of the problem

The observational impact of a clumpy stellar wind strongly depends on the density contrast in the two-wind interaction region, which gives the clump enhancement in momentum flux. The sound speed of the shocked stellar wind when it is smooth is at the region where its ram pressure is balanced by the pulsar wind, whereas the sound speed of a shocked clump at the same location gets reduced by because of the higher density (see e.g. Araudo et al., 2009, and references therein), so the time required to shock the clump when it meets the interaction region is

 tc≈f−1/2Rc/Vw. (1)

For a clump radius such that the clumps will be destroyed and deflected by the shocked wind medium well before crossing it, being quickly integrated and carried away from the binary (Pittard, 2007). Otherwise, for

 Rc≳f1/2ΔR, (2)

the clumps will not be smoothed out within the shocked two-wind zone, but will penetrate into the unshocked pulsar wind. Depending on and the average clump size, the two-wind interaction region as a whole will be smooth or very inhomogeneous. In the latter case, the shape and location of the two-wind interaction region will vary in time.

Under typical pulsar-to-star wind momentum-flux ratios, say (see Sect. 1 in Bosch-Ramon et al., 2012a), the total luminosity injected by the pulsar wind within the binary will be , times larger than that of the star. This suggests that the pulsar wind shock should largely dominate the non-thermal output from the two-wind interaction region. In addition, determines the magnetic field through for a dominant toroidal -component, and the cooling time for escape / adiabatic (non-radiative) cooling. Thus, clump-induced variations in the pulsar wind shock could strongly affect the properties of the non-thermal emission, namely the broadband spectrum and variability.

Considering the solid angle covered by the clump as seen from the pulsar, the luminosity of a clump-pulsar wind shock will be smaller by than the whole two-wind interaction region, where , and is the distance the clump can cover within the unshocked pulsar wind:

 X∼tcVw∼f−1/2Rc. (3)

To study the case , one can write

 (Rc/R′)2∼(RcR−X)2∼(RcR−f−1/2Rc)2 (4)

to order 1 if

 Rc∼R/(1+f−1/2)∼f1/2R. (5)

Equation (5) provides a stricter condition than that derived from Eq. (2) by a factor of several, given the typical thickness of the two-wind interaction region right on the line joining the two stars (see e.g. Bogovalov et al., 2008). Clumps fulfilling this condition will be able to penetrate up to the location at which the pulsar wind ram pressure equals that of the clump. This yields for these long-lived clumps

 R′∼f1/2R∼Rc, (6)

which means that most of the pulsar wind will terminate in the clump. At that point, the clump will stop its advance and dissipate under the impact of the pulsar wind. This process is discussed in Sect. 5.

### 3.3 Clump size

As described in Sect. 2, the stellar wind clumps form through complex non-linear, radiation-hydrodynamical processes in the stellar atmosphere. Although empirical and theoretical clump characterization is still under development, clumps can be considered dense, changing structures embedded in a hotter and very dilute medium. Clumps therefore expand and evolve in the stellar wind in a complex manner. In particular, and thereby the density of the clumps, as well as , seem to be functions of (e.g. Puls et al., 2006), not only because clumps expand as they propagate away from the star, but also because of collisions between clumps that produce their merging or fragmentation (e.g. Owocki & Cohen, 2006).

In this work we are interested in those known binaries that host or may host a non-accreting pulsar (see Paredes et al., 2013; Dubus, 2013, and references therein). These binaries are relatively compact, with periods ranging from days to years, and their orbital semi-major axis typically going from a few to a few cm. To investigate the impact of a clumpy wind, we will normalize , where is the orbital separation distance and the stellar radius, fixed to  cm. The parameter will be normalized at .

The pulsar-to-star wind momentum flux ratio is also important here. Typically, (see Sect. 3.2), so we can write

 R=√ηd/(1+√η)≈η1/2d≈3η1/2−1d1R∗, (7)

where . For the error of the approximation is and roughly decreases as for smaller -values. Equations (6) and (7) provide a minimum radius for the clump to deeply penetrate into the unshocked pulsar wind:

 Rc∼f1/2−1η1/2−1d1R∗. (8)

The possibility of having such large structures from instabilities in the wind acceleration region deserves more careful considerations. Assuming clump-free expansion , Eq. (8) implies clump sizes at the wind base of . This value seems too large if one takes the Sobolev length as the typical clump size at the wind base, of , where is the stellar wind ion thermal speed (see e.g. Owocki & Cohen, 2006). The typical or average clump size at the two-wind interaction region then becomes , hence to fulfill Eq. (8) would require . This would only be possible if were much smaller than the values typically considered, and/or were very small, which would imply a very weak pulsar wind and therefore a faint non-thermal source. Furthermore, clumps originating in the stellar wind acceleration region could be disrupted by Rayleigh-Taylor or thin-shell instabilities, limiting the clump growth when travelling away from the star (Dessart & Owocki, 2003, 2005). It is worth noting that large clump-like structures may still form in the stellar wind for reasons different from instabilities occurring during the line-driven acceleration of the wind. In particular, the stellar rotation, magnetic field inhomogeneities in the stellar surface, or non-radial pulsations, may lead to clump-like, large-scale structures in the wind of size , like those inferred from the observed DACs in the UV. It is unclear whether the density contrast behind these observational UV features could be as large as (Cranmer & Owocki, 1996; Lobel & Blomme, 2008). Nevertheless, the complexity of the stellar wind dynamics prevents us from discarding this possibility, as well as that of sporadic large and/or dense clumps of wind-acceleration region origin.

### 3.4 Clump-pulsar wind interaction rates

An estimate of how often clumps fulfilling Eq. (8) reach the two-wind interaction region can be calculated regardless of the clump origin. One can determine the duty cycle () of these events by deriving their frequency () and multiplying it by their duration , and also by adopting an empirical clump mass distribution without specifying the physical origin of those clumps.

To find , one needs first the number density distribution of clumps with radius

 dNdRcdV=n(Rc). (9)

Assuming that all clumps are spherical and have the same density, the distribution function should meet the following normalization condition:

 4π3∞∫0R3cn(Rc)dRc=f. (10)

Adopting for simplicity for , and otherwise, Eq. (10) yields

 n(Rc)=3f(α−4)4πRα−4c,0R−αc, (11)

and the number density of clumps exceeding the limit given by Eq. (8) is

 Υ=∞∫f1/2Rn(Rc)dRc=3f(α−4)4π(α−1)Rα−4c,0(f1/2R)1−α. (12)

We note that the derivation is valid for , since otherwise the integration in Eq. (10) would diverge. This -range also means that the wind mass is mostly in the smallest clumps, of size .

A large enough clump can strongly affect the emission from the two-wind interaction region if it propagates in the direction of the pulsar within a solid angle . Therefore, the frequency of arrival of large clumps to the vicinity of the pulsar can be estimated as

 ˙N=ΔΩd2VwΥ, (13)

which allows the derivation of :

 DC=tc˙N=f−1/2RcΔΩd2Υ. (14)

We note that this is only approximately true because is a function of clump size and is an average. Taking an empirical mass power-law index , , and  cm s (typical for massive star winds; Muijres et al. 2012), one can obtain and the typical time between large clump arrivals  s.

Clumps fulfilling Eq. (5) are expected to be rare, implying that their interaction with the pulsar wind should appear as a distinguishable event. Nevertheless, albeit unlikely, it cannot be discarded that clumps fulfilling Eqs. (2) or even (5) could have , in which case the simultaneous arrival of many clumps would smooth out the associated variability (see also Owocki et al., 2009).

## 4 Radiation from clumpy stellar/pulsar wind collisions

### 4.1 Main characteristics of the clump-pulsar wind shock radiation

In the clump-pulsar wind region two shocks form. One in the pulsar wind, and one in the clump. Non-thermal high-energy radiation, of a synchrotron and inverse Compton (IC) nature, is expected to originate in the shocked pulsar wind. Relativistic bremsstrahlung is negligible in this light flow. The shock in the much denser clump has a luminosity smaller by than that of the pulsar wind shock. This is much smaller than unity, and implies that particle acceleration within the clump will be rather inefficient. The highest energy electrons may still penetrate in the clump from the shocked pulsar wind through diffusion and radiate there, but to characterize this process, fast radiative and non-radiative losses and diffusive transport should be carefully accounted for. Protons of pulsar wind origin have not been considered in this work.

#### 4.1.1 Dynamics and particle cooling

Two extreme cases illustrate the impact of a clump interacting with the two-wind interaction region: (i) and , or (ii) .

For case (i), will grow by a factor of in the clump-pulsar wind interaction region. Given the relation between and the synchrotron cooling rate , the cooling rate will increase by . Assuming an injection non-thermal luminosity in the clump-pulsar wind interaction region, where is the total non-thermal luminosity, the maximum variation of the synchrotron emission will be . For , , , and (so ), the expected fluctuations will be at the level of %. The variations in IC and non-radiative cooling will be similar or even smaller.

As becomes (case ii), the situation changes quickly. From that point on, a substantial fraction of the pulsar wind luminosity crosses the clump-pulsar wind shock, and synchrotron cooling is now enhanced by a factor of (as ). On the other hand, IC cooling only changes slightly for . Non-radiative cooling grows by . As the loss balance changes, radiation is re-distributed in photon energy.

#### 4.1.2 Non-thermal emission

The rise time of the lightcurve of the emission produced during a clump-pulsar wind interaction is linked to the non-radiative timescale (see Sect. 3.1):

to the acceleration timescale (, ; with the electron energy and the acceleration efficiency):

 tacc≈1ETeVξ−10.1B−1Gs; (16)

to the synchrotron timescale:

 tsync≈400B−1GE−1TeVs; (17)

to the IC timescale:

 tIC≈160u−12E−1GeVs (18)

in the Thomson regime (, ), and

 tIC≈100E0.7TeVT4.5u−12s, (19)

in the Klein-Nishina regime () (for radiation and acceleration timescales, see Khangulyan et al., 2008); and to the clump expansion timescale:

 ∼tc∼104f−1/2−1(Rc/R∗)s. (20)

Since the particle energy gain and loss timescales will be generally much shorter than , the accelerated particles will reach the steady state, or an adiabatic succession of quasi-steady states, before the clump gets destroyed by the pulsar wind. The lightcurve will therefore be strongly linked to the clump dynamics.

Clump expansion is rapid and occurs simultaneously with clump braking and deflection by the pulsar wind. As illustrated in Fig. 1, the clump ends up as a tail of shocked material that will join the tail of the whole two-wind interaction region. For an accurate modelling of the lightcurve, the complexity of the clump destruction deserves a detailed treatment that is beyond the scope of this work. Nevertheless, we note that for the case , the emission from the clump-pulsar wind interaction can be additionally enhanced by a factor of a few during the clump expansion with a sharpening of the lightcurve (see e.g. Barkov et al., 2010; Bosch-Ramon et al., 2012b, and Sect. 5, for the cloud radius evolution).

The balance between the different cooling channels changes with the clump arrival, which, putting aside non-radiative cooling, diverts radiation energy from the highest energies (IC) to lower ones (synchrotron). The synchrotron and the IC luminosity can be approximately described by

before the clump arrival. After the arrival of a clump with , the clump-pulsar wind interaction region emits synchrotron radiation as

and IC radiation as

The cooling balance shows that, when a large clump reaches the pulsar vicinity, the synchrotron emission grows in the clump-pulsar wind interaction region unless it was the dominant cooling channel before the clump arrival. On the other hand, the IC emission always diminishes. Non-radiative losses, which grow as well but more slowly than synchrotron radiation, moderate when they are dominant the -dependence of the two radiation channels.

Inverse Compton scattering of electrons with their own synchrotron photons may also become important as the emitter size decreases from to . The same applies to gamma-ray absorption through pair creation in the synchrotron field. One can estimate the influence of size reduction comparing IC scattering with stellar photons and synchrotron self-Compton (SSC) in the Thomson regime approximation, the former being dominant unless , in which case the internal photon energy density can overcome that of the star. This requirement can be rewritten as

 Lsd>1036η−1f−1L∗38erg s−1, (24)

or

 ˙MVwc>1037f−1L∗38erg s−1, (25)

where is the wind mass-loss rate and , and may be met in sources such as PSR B125963. A similar condition applies for significant pair creation in the synchrotron field given that the pair-creation cross section is similar to that of IC scattering in the Klein-Nishina regime. The stellar and the synchrotron spectra are however different, and thus the transition between the Thomson and the Klein-Nishina regime (e.g. Blumenthal & Gould, 1970) occurs differently in both target fields. This introduces some uncertainty in the estimated constraint, but given that the synchrotron photons in the region are more energetic than the stellar photons, the constraint would actually be stricter.

### 4.2 Semi-analytic modelling of an illustrative case

In this section, we present the results of semi-analytical calculations of the synchrotron, synchrotron self-Compton, stellar photon IC emission, and external (binary) and internal (emitter) pair creation and its subsequent pair radiation. For these calculations, we focused on the case with , , , , all this yielding . We also took , and fixed to  erg s, taking it as equal to the injected non-thermal power (i.e. ). We also adopted two magnetic-to-total energy density ratios in the shocked pulsar wind of and 1, giving an initial magnetic field in the two-wind interaction region of and 3 G. Therefore, after the clump arrival, becomes and 10 G in the clump-pulsar wind interaction region; these two magnetic field values were adopted as fiducial values for a kinetic and a magnetic pulsar wind, respectively. We took an acceleration efficiency in the clump-pulsar wind shock of to obtain the maximum particle energy , with an injection energy distribution . The emitter was assumed to be homogeneous (one zone), with size and and a non-radiative timescale and before and after the clump arrival. The stellar luminosity and temperature were taken  erg s and  K, respectively, with the star being considered a black body and point-like.

To avoid the complexities of IC scattering and gamma-ray absorption from the orbital motion and the emitter-observer geometry relations (e.g. Khangulyan et al., 2008), not relevant at this stage, we adopted an isotropic stellar field. Moreover, we computed the secondary emission in the one-zone approximation. For the pairs created within the system, we took an emitter of size with a radiation energy density fixed to its value at a distance from the star. The -value in the system was taken 10 times below equipartition with the stellar photon field. However, as expected (see Sect. 4.1.2), both internal and external absorption yielded a relatively small number of pairs. Therefore, secondary emission, as well as SSC, well below the external IC component, are not considered further.

We present in Figs. 2 and 3 the broadband spectral energy distribution (SED) of the synchrotron and stellar IC radiation with gamma-ray absorption by the stellar photons. The pre- and post-clump SEDs are presented together, each figure corresponding to a different magnetic field: (Fig. 2) and 1 (Fig. 3). Results are similar for both cases, although the radiative cooling regime starts at higher energies for the lower case. Below the highest energies of the synchrotron and IC components and for the adopted parameters, non-radiative cooling shapes the SED, so the whole synchrotron luminosity grows and the IC emission decreases by the same factor after the clump arrival. At high photon energies of both radiation components, synchrotron cooling shapes the SED before and after the clump arrival in both cases, and therefore the synchrotron luminosity does not change, whereas the IC radiation decreases (hard to see because of the steepness of the SED). We note that for a fixed energy band, synchrotron emission is actually (for a spectral index , where ), so when balanced by non-radiative cooling one gets an X-ray flux enhancement , i.e. a factor of for . This does not occur in IC emission. Another spectral feature is that the SED tends to soften for larger magnetic fields.

For the studied case, non-radiative losses reduce the synchrotron enhancement in X-rays from to a factor of a few. The reduction of IC emission at GeV energies is also moderated by the dominant non-radiative losses. Given that non-radiative losses are , and synchrotron and IC losses , for more compact binaries radiation cooling gets more prominent, and the contrast between luminosities before and after the clump arrival becomes higher at lower energies. For wide systems, on the other hand, non-radiative losses can be dominant in the whole particle energy range.

## 5 Discussion and summary

The structure of stellar winds in massive binaries hosting a powerful pulsar can have significant effects in the non-thermal lightcurve and spectrum. The two dominant factors are the increase in , which directly affects the synchrotron emission and indirectly, the IC emission, and the shorter non-radiative cooling time. Nevertheless, the synchrotron losses grow faster, , than the non-radiative losses, . Inverse Compton cooling is only slightly affected for .

An open issue is whether clumps large and/or dense enough to penetrate deeply in the pulsar wind region can form in the stellar wind. Instabilities in the stellar wind acceleration region could be a possibility, but the related clumps may be too small if their growth is limited by clump disruption beyond some size. The large-scale stellar wind structures behind DACs may be another possibility, although it is unclear whether the density contrast in this case is enough to fulfill the pulsar wind deep-penetration condition.

Regardless of the physical origin of the clumps, if the number density distribution of clumps with radius can be approximated by an empirical power-law at least up to , some clumps can sporadically collide with the two-wind interaction region inducing variations in the emission by a factor of several on timescales of minutes to hours (more on clump expansion below). The duty cycles strongly depend on the structure of the stellar wind, but duty cycles of % seem feasible.

There are other issues that deserve discussion in the studied scenario. We comment on a number of them in what follows: (i) the clump expansion; (ii) the Doppler boosting of the clump radiation when the shocked clump material approaches relativistic speeds; (iii) the effect of clump fragmentation; (iv) the magnetic field; (v) the clump thermal radiation; and (vi) observations.

Clump expansion (i) in the vicinity of the pulsar can be very fast. Using Eqs. (9)-(13) in Barkov et al. (2010) and accounting for the radial geometry of the pulsar wind, we calculated the cloud expansion until the cloud advance is stopped by the pulsar wind, which result is shown in Fig. 4. As seen in the figure, reaches the size of the whole two-wind interaction region around the point where the clump stops. The result is consistent with the shocked clump material eventually forming its own bow-shaped structure. The calculation also shows that the clump expansion accelerates, and its last stages may take minutes after an initial phase lasting about one hour, under the adopted parameter values. These last stages of the expansion may correspond to a sharpening of the lightcurve. The value of found in the calculation is slightly different to that derived analytically adopting the same parameters. However, the calculation of the cloud evolution has been largely simplified, as instabilities and clump fragmentation have been neglected. More detailed numerical simulations of a relativistic outflow impacting a spherical cloud (e.g. Bosch-Ramon et al., 2012b) show that the clump expansion rate, albeit high, is somewhat slower than the one computed here, and more in line with the analytical estimates.

Regarding Doppler boosting (ii), the clump expansion and propagation velocities do not become relativistic before dissipation, although the shocked pulsar wind circumventing the clump should have an accelerating relativistic speed (Bogovalov et al., 2008). Finer calculations should account for this effect, which can introduce line-of-sight dependencies, i.e. an additional variability factor.

The last stages of the clump-pulsar wind interaction may be accompanied by features in the lightcurve and the spectrum linked to cloud fragmentation (iii). Once disrupted, the clump material becomes a slow and dense inhomogeneous flow that joins the shocked two-wind structure and is advected by it. This will excite shocks, turbulence and two-wind mixing in the flow (as suggested, e.g., in Zdziarski et al. 2010; Lamberts et al. 2012), which are in fact already present when a smooth stellar wind interacts with the relativistic pulsar wind (Bosch-Ramon et al., 2012a).

The clump propagation through the shocked stellar and pulsar wind can strongly interfere in the magnetic field advected downstream (iv). The -lines will be warped and stretched by the clump, which will lead to -growth, instabilities, turbulence, and reconnection in the plasma trailing the clump (e.g. Jones et al., 1996); all these effects are potentially relevant for the non-thermal activity around the clump.

Modest thermal emission (v) in the form of soft X-ray lines and thermal Bremsstrahlung is expected from the clump because of the impact of the pulsar wind. Initially, the shocked clump produces thermal photons times cooler than those of the shocked smooth wind, which mainly emits in the soft X-ray band (see e.g. Zabalza et al., 2011). As clump heating and expansion proceed, the thermal photon energy reaches that of the shocked smooth wind. The clump gets then hotter but also significantly diluted, quenching its thermal radiation.

Certain observational features (vi) detected in binaries potentially hosting a pulsar may be explained by invoking the presence of clumps, for instance short flares of scales of seconds to hours found in the X-ray lightcurves of LS 5039 and LS I +61 303 (e.g. Bosch-Ramon et al., 2005; Paredes et al., 2007; Smith et al., 2009; Li et al., 2011). The observational hints of stellar wind inhomogeneities, manifesting themselves mostly through short X-ray variability, and the complex evolution of the wind structures under the impact of the pulsar wind, call for coupled dynamics-radiation studies and their comparison with observations. Besides improving our understanding of the non-thermal processes within the binary, these studies can also probe the stellar wind properties in high-mass binaries hosting young pulsars.

We want to emphasize that direct attempts to simulate the processes considered approximately in this work are very necessary, since as noted an analytic treatment does not allow for precise quantitative predictions and may miss subtle effects of major impact. Despite our conclusions pointing to a smooth effect of moderate clumping on the two-wind post-shock regions (see also Pittard, 2007), the relativistic nature of the pulsar wind (about the related complexities, see e.g. the discussion in Bosch-Ramon et al., 2012a), and the range of possible interaction parameters of the sources of interest, open the door to unexpected effects even when only modest size clumps could be present in the stellar wind. In addition, the presence of additional larger inhomogeneities in the stellar wind (like those behind DACs) could provide a strong source of variability, and also introduce substantial changes not only in the size, but also in the orientation and the overall structure of the two-wind interaction region. Finally, it is noteworthy that PSR B125963, the only confirmed high-mass binary hosting a powerful pulsar so far, has a circumstellar disc, which adds further complications to the evolution of the two-wind interaction structure (see e.g. Okazaki et al., 2011; Takata et al., 2012). In short, numerical three-dimensional relativistic hydrodynamical simulations of a stellar wind interacting with a pulsar wind, accounting for both modest and large wind structures, are needed for a more accurate assessment of the required clump size to have a significant observational impact. This should be done at least on a spatial scale , and the inclusion of a magnetic field would be useful. The required size of the grid plus a minimum resolution of a few cells per mean a rather heavy calculation, although box-type simulations on spatial scales could already be quite informative. Detailed clump time-arrival calculations, as done in Oskinova et al. (2012), would be also informative to interpret lightcurves.

###### Acknowledgements.
We want to thank Stan Owocki for a constructive and really useful review. We also thank Dmitry Khangulyan for his useful comments on the manuscript. We acknowledge support by DGI of the Spanish Ministerio de Economía y Competitividad (MINECO) under grants AYA2010-21782-C03-01 and FPA2010-22056-C06-02. V.B-R. acknowledges financial support from MINECO through a Ramón y Cajal fellowship. This research has been supported by the Marie Curie Career Integration Grant 321520.

## References

• Araudo et al. (2009) Araudo, A. T., Bosch-Ramon, V., & Romero, G. E. 2009, A&A, 503, 673
• Barkov et al. (2010) Barkov, M. V., Aharonian, F. A., & Bosch-Ramon, V. 2010, ApJ, 724, 1517
• Blumenthal & Gould (1970) Blumenthal, G. R. & Gould, R. J. 1970, Reviews of Modern Physics, 42, 237
• Bogovalov et al. (2008) Bogovalov, S. V., Khangulyan, D. V., Koldoba, A. V., Ustyugova, G. V., & Aharonian, F. A. 2008, MNRAS, 387, 63
• Bosch-Ramon (2010) Bosch-Ramon, V. 2010, in 25th Texas Symposium on Relativistic Astrophysics
• Bosch-Ramon et al. (2012a) Bosch-Ramon, V., Barkov, M. V., Khangulyan, D., & Perucho, M. 2012a, A&A, 544, A59
• Bosch-Ramon et al. (2005) Bosch-Ramon, V., Paredes, J. M., Ribó, M., et al. 2005, ApJ, 628, 388
• Bosch-Ramon et al. (2012b) Bosch-Ramon, V., Perucho, M., & Barkov, M. V. 2012b, A&A, 539, A69
• Cranmer & Owocki (1996) Cranmer, S. R. & Owocki, S. P. 1996, ApJ, 462, 469
• Dessart & Owocki (2003) Dessart, L. & Owocki, S. P. 2003, A&A, 406, L1
• Dessart & Owocki (2005) Dessart, L. & Owocki, S. P. 2005, A&A, 437, 657
• Dubus (2006) Dubus, G. 2006, A&A, 456, 801
• Dubus (2013) Dubus, G. 2013, ArXiv e-prints
• Hirschi (2008) Hirschi, R. 2008, in Clumping in Hot-Star Winds, ed. W.-R. Hamann, A. Feldmeier, & L. M. Oskinova, 9
• Johnston et al. (1992) Johnston, S., Manchester, R. N., Lyne, A. G., et al. 1992, ApJ, 387, L37
• Jones et al. (1996) Jones, T. W., Ryu, D., & Tregillis, I. L. 1996, ApJ, 473, 365
• Khangulyan et al. (2008) Khangulyan, D., Aharonian, F., & Bosch-Ramon, V. 2008, MNRAS, 383, 467
• Khangulyan et al. (2007) Khangulyan, D., Hnatic, S., Aharonian, F., & Bogovalov, S. 2007, MNRAS, 380, 320
• Kirk et al. (1999) Kirk, J. G., Ball, L., & Skjaeraasen, O. 1999, Astroparticle Physics, 10, 31
• Lamberts et al. (2012) Lamberts, A., Dubus, G., Lesur, G., & Fromang, S. 2012, A&A, 546, A60
• Lamberts et al. (2013) Lamberts, A., Fromang, S., Dubus, G., & Teyssier, R. 2013, ArXiv e-prints
• Li et al. (2011) Li, J., Torres, D. F., Zhang, S., et al. 2011, ApJ, 733, 89
• Lobel & Blomme (2008) Lobel, A. & Blomme, R. 2008, ApJ, 678, 408
• Lucy & Solomon (1970) Lucy, L. B. & Solomon, P. M. 1970, ApJ, 159, 879
• Maraschi & Treves (1981) Maraschi, L. & Treves, A. 1981, MNRAS, 194, 1P
• Moffat (2008) Moffat, A. F. J. 2008, in Clumping in Hot-Star Winds, ed. W.-R. Hamann, A. Feldmeier, & L. M. Oskinova, 17
• Moldón et al. (2011) Moldón, J., Johnston, S., Ribó, M., Paredes, J. M., & Deller, A. T. 2011, ApJ, 732, L10
• Muijres et al. (2012) Muijres, L. E., Vink, J. S., de Koter, A., Müller, P. E., & Langer, N. 2012, A&A, 537, A37
• Negueruela et al. (2011) Negueruela, I., Ribó, M., Herrero, A., et al. 2011, ApJ, 732, L11
• Neronov & Chernyakova (2007) Neronov, A. & Chernyakova, M. 2007, Ap&SS, 309, 253
• Okazaki et al. (2011) Okazaki, A. T., Nagataki, S., Naito, T., et al. 2011, PASJ, 63, 893
• Oskinova et al. (2012) Oskinova, L. M., Feldmeier, A., & Kretschmar, P. 2012, MNRAS, 421, 2820
• Owocki & Cohen (2006) Owocki, S. P. & Cohen, D. H. 2006, ApJ, 648, 565
• Owocki et al. (2009) Owocki, S. P., Romero, G. E., Townsend, R. H. D., & Araudo, A. T. 2009, ApJ, 696, 690
• Papitto et al. (2012) Papitto, A., Torres, D. F., & Rea, N. 2012, ApJ, 756, 188
• Paredes et al. (2013) Paredes, J. M., Bednarek, W., Bordas, P., et al. 2013, Astroparticle Physics, 43, 301
• Paredes et al. (2007) Paredes, J. M., Ribó, M., Bosch-Ramon, V., et al. 2007, ApJ, 664, L39
• Perucho & Bosch-Ramon (2012) Perucho, M. & Bosch-Ramon, V. 2012, A&A, 539, A57
• Pittard (2007) Pittard, J. M. 2007, ApJ, 660, L141
• Portegies Zwart & Yungelson (1999) Portegies Zwart, S. F. & Yungelson, L. R. 1999, MNRAS, 309, 26
• Puls et al. (2006) Puls, J., Markova, N., Scuderi, S., et al. 2006, A&A, 454, 625
• Runacres & Owocki (2002) Runacres, M. C. & Owocki, S. P. 2002, A&A, 381, 1015
• Smith et al. (2009) Smith, A., Kaaret, P., Holder, J., et al. 2009, ApJ, 693, 1621
• Takata et al. (2012) Takata, J., Okazaki, A. T., Nagataki, S., et al. 2012, ApJ, 750, 70
• Tavani & Arons (1997) Tavani, M. & Arons, J. 1997, ApJ, 477, 439
• Walter & Zurita Heras (2007) Walter, R. & Zurita Heras, J. 2007, A&A, 476, 335
• Zabalza et al. (2013) Zabalza, V., Bosch-Ramon, V., Aharonian, F., & Khangulyan, D. 2013, A&A, 551, A17
• Zabalza et al. (2011) Zabalza, V., Bosch-Ramon, V., & Paredes, J. M. 2011, ApJ, 743, 7
• Zdziarski et al. (2010) Zdziarski, A. A., Neronov, A., & Chernyakova, M. 2010, MNRAS, 403, 1873
You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters