1. Get photon density of upscattered photons at the base from the model (do not consider synchrotron because the synchrotron

spectrum from the thermal population of electrons/positrons cuts off sharply at energies well below 511 keV, thus they cannot produce any pairs). The figure below shows the photon spectrum for the best fits of J1118 and GX339 reported in the paper.

The relevant model parameters in these two cases were:

J1118 | GX339 | |

$D_{\rm src}$ [kpc] | 1.7 | 6 |

$r_0/c$ [s] | 5e-4 | 1e-3 |

$T_{\rm thermal\ particles}$ [K] | 4.0e+10 | 6.75e+10 |

$n_0$ [#/$cm^{-3}$ ] | 0.6e+14 | 1.8e+14 |

2. Estimate the pair annihilation rate for a given fit.

3. Estimate the pair production rate for a given fit.

4. Results and interpretation: The total rates of pp, pa [Unit: #/cm^3/s] and the fractional rates of pp ($\dot n_{pp}/n_0$) and pa ($\dot n_{pa}/n_0$), calculated using the above procedure gives the following numbers. The fractional rates $\dot n/n_0$ give an estimate of the rate at which particles are gained/lost at the jet-base.

J1118 | GX339 | |

Pair production rate $\dot n_{pp}$ [#/cm^3/s] | $\sim 10^7$ | $\sim 10^{17}$ |

Pair annihilation rate $\dot n_{pa}$ [#/cm^3/s] | $\sim 10^{11}$ | $\sim 10^{12}$ |

Fractional pair production rate $\dot n_{pp}/n_0$ [1/s] | $\sim 10^{-6}$ | $\sim 10^3$ |

Fractional pair annihilation rate $\dot n_{pa}/n_0$ [1/s] | $\sim 10^{-2}$ | $\sim 10^{-2}$ |

For J1118 $\dot n_{pp} \ll \dot n_{pa}$, suggesting that any pair that is created will be quickly annihilated. Also the fact that

$\dot n_{pa}/n_0$ is small implies that the total number density is not affected very much by pair processes. Thus for J1118 (I think) we are safe.

However for GX339 $\dot n_{pp} \gg \dot n_{pa}$ suggesting that pp dominates over pa. Also the fact that $\dot n_{pa}/n_0 \gg 1$ implies that particles produced by pp will quickly dominate the overall number density at the base. Even taking into account the fact that the residence time of the plasma near the jet base $\tau_r \sim R/(\beta_jc) \sim 10^{-3}$s, it looks like pp cannot be ignored for GX339. While I am exploring the model parameter space to see if we can find another minima with lower pp rate, it would be good to be sure that the way these numbers we calculated are reasonable.

—-

**Added on 2009 Apr 02**

Taking the values of the model parameters from 'fit 1' of the submitted paper, froze everything but T_e and r_0. The plots below (click on the plot to see full size) show the variation of number density, pair production rate and pair annihilation rate at the base of the jet.

Here's the Dove et al. paper.

—-