| Outburst Modelling in Astrophysical Jets |
High-resolution radio observations of astrophysical jets show moving structures or blobs thought to be shock waves propagating down the jet. Such disturbed outflows can be modelled by hydrodynamical simulations, but we present here a different approach, which is intermediate between observations and simulations.
The method consists in a decomposition of the observed total flux variations into a series of self-similar outbursts evolving with time according to a physical shock model. The decomposition uses as many observational constraints as possible from flux measurements at different epochs and wavelengths. The outburst's properties are derived by fitting the model parameters simultaneously to all observations.
This original approach aims at testing the validity of the model and at constraining the physical properties of the jet. It has been applied to both galactic jets from black-hole binaries (micro-quasars) and extra-galactic jets from distant powerful quasars and blazars. The results have been the subject of several publications listed below. We present here a brief description of the model, the method and some results illustrated by a selection of interesting complementary material - mainly animations - which cannot be included in scientific publications.
Figure: Animation of simulated high-resolution jet observations of the blazar 3C 279 for the 1990's decade. The simulation is based on the results of the decomposition of the optical-to-radio light-curves of the source. It shows the evolution of blobs at an observing frequency of 22 GHz and with a resolution (FWHM) of 0.15 milli-arcsecond (mas). We assume that all blobs are moving with a constant (for simplicity) apparent superluminal velocity of 5 times the speed of light. The color-coding and contours represent the observed flux on a square-root scale, chosen to give more emphasis to low-flux features. Lindfors et al. (2006)
The physical model is a generalization of the
shock-in-jet model of Marscher & Gear (1985). The basic principle of this model
is as follows. A shock wave is propagating along the jet. Electrons crossing the
shock front are accelerated. They emit synchrotron radiation in the shocked
medium behind the shock front. The emitted synchrotron spectrum (blue curve)
evolves following 3 stages depending on the dominant emission-loss mechanism of
the electrons: 1. Compton losses, 2. synchrotron losses and 3. adiabatic
(non-radiative) losses. The self-absorption turnover of the spectrum follows a
characteristic 3-stage pattern (red curve). While the turnover frequency is
steadily evolving from high to low freqencies, the turnover flux is first rising
during the Compton-loss stage, stays more or less constant during the
synchrotron-loss stage and decays during final adiabatic-loss stage. According
to the model, the optically thin spectral slope is expected to be steeper
(-s/2) during the first two radiative-loss stages than during the final
decay (-(s-1)/2), where s is the index of the electron energy
distribution N(E) ~ E-s.
The model presented above defines fully the evolution
of a synchrotron outburst resulting from a shock wave at any frequency
(wavelength) and at any time. The evolution of the spectrum with time shown
above defines therefore also the shape of the light-curve (flux evolution with
time) at any given frequency. For observational reasons (many observations at
some fixed wavelengths), it is easier to identify and characterize the outbursts
in the time-domain (light-curves) than in the frequency-domain (spectra). This
is done by a simultaneous fit of the light-curves available at as many different
frequencies as possible to have the best constraints on the model. In practice,
many iterative fits on subsets of all parameters are needed to achieve a good
match to the observations. The result is a best fit decomposition of the
light-curves into a series of self-similar outbursts as shown in the figure. The
fit parameters define both the average properties of a typical outburst and the
specific characteristics of each individual outburst. They give valuable clues
on the physical properties of the jet.
Figure: Long-term (20 years) radio-to-sub-millimetre light-curves of the quasar 3C 273 with the best fit decomposition into self-similar outbursts. The total model light-curve (grey line) is the sum of the flux of the individual outbursts (colored areas), the global decay of the outbursts peaking before 1979 (light grey), a possible constant emission from the underlying jet (medium grey) and the steady (fixed) emission from the outer jet extending on kiloparsec-scales (dark grey). Türler et al. (2000); Türler & Lindfors (2006)
The decomposition of the lightcurves into separate
model outbursts allows to produce simulated jet images at any frequency and
angular resolution. The model flux of each component is known at any frequency
and time and therefore we just have to assume a Gaussian beam and a constant
apparent speed of each knot in the jet to construct the simulated jet images at
any time. The link between the outbursts in the lightcurve and the
conrresponding structures moving in the jet is shown by the use of different
colors in the animated figure. A nicer view of the simulated jet animation of 3C
279 is shown at the top of this page.
Figure: Animation showing the link between moving jet structures and outbursts in the lightcurve of 3C 273 at 22 GHz. The simulated jet image (contours) at the top is derived from the outbursts identified in the lightcurve. Colors are used to relate the outbursts to the moving jet structures. The steady core emission is shown in grey. The apparent speed of the knots in the jet is assumed to be the same for all outbursts and to be constant in time.
The method first used for quasars like 3C 273 with a single observable jet was developped further to include as well the emission of the counter-jet. This was needed to model the outbursts of the galactic black hole binary GRS 1915+105. For this microquasar, we assumed twin outbursts propagating in two opposite jets. The decomposition of the very well sampled lightcurves enabled to derive a mildly relativistic jet speed of 0.6c during the observations of 15 May 1997 corresponding to a "baby-jet" phase (Türler et al. 2004).
Figure: Light-curves of the micro-quasar GRS 1915+105 taken on 15 May 1997 with the best fit decomposition into self-similar outbursts. The total model light-curve (red curve) with its distinct outbursts (red dotted lines) is the sum of the light-curves from the approaching jet (blue) and the receiding jet (green) plus a constant flux in the infrared at 2.2 microns. The grey line in the upper panel is the X-ray light-curve in the 2-60 keV observed by RXTE. Arrows show the onset time of the individal outbursts. Türler et al. (2004)
An interesting area of development of the modelling presented here is the inclusion of the associated synchrotron self-Compton (SSC) emission in the X-ray and gamma-ray range. This is motivated by growing evidence that X-ray and gamma-ray flaring is related to synchrotron outbursts observed at lower frequencies (Lindfors et al. 2006, and references therein). A full self-consistent treatment of the SSC mechanism requires numerical calculations. It is hower possible to do some simple modelling of the X-ray and gamma-ray emission with basic assumptions. In particular, the shape of the Compton spectrum was simply taken to be that of the synchrotron spectrum by squared, thus appearing scaled-up by a factor of two on logarithmic-scales (Lindfors et al. 2005).
Figure: Animation of model SSC spectra based on the outbursts identified and fitted to match the optical-to-radio light-curves of 3C 279. The shape of the outburst spectrum at X-ray and gamma-ray energies is simply assumed to be self-similar at any time to the synchrotron spectrum of each outburst (black lines). The blue spectra represent emission from the quiescent jet and the magenta curve the decaying combined spectrum of outbursts flaring before the start of the modelling in 1979. The overall spectrum is shown by the red curve and is compared to data points corresponding to the average observed spectrum of 3C 279 with error bars showing the typical dispersion due to the source variability. Lindfors et al. (2005)
In our last paper on Cyg X-3 (Lindfors et al. 2007) we included the effect of the modifications to the Marscher & Gear (1985) model proposed by Björnsson & Aslaksen (2000) concerning the initial Compton-loss stage. In this paper and also in Türler & Lindfors (2007) we derive the physical properties of the jet (magnetic field, electron energy density, etc.) at the point where a typical outburst makes its transition to the adiabatic loss stage.
In the future, we plan to further investigate the possibility discussed in Türler & Lindfors (2007) that during the early evolution of the outbursts, the observed spectral turnover might not be due to synchrotron self-absorption, but to a minimal energy of the electron spectrum. On a somewhat longer timescale, we wish to make significant progress in the treatment of the associated synchrotron self-Compton (SSC) emission based - at least in part - on the analytic description proposed by Björnsson & Aslaksen (2000).