A&A 463, 145-152 (2007)
DOI: 10.1051/0004-6361:20066502
O. Tsang - J. G. Kirk
Max-Planck-Institut-für Kernphysik, Saupfercheckweg 1, 69117 Heidelberg, Germany
Received 3 October 2006 / Accepted 15 November 2006
Abstract
Context. The inverse Compton catastrophe is the dramatic rise in the luminosity of inverse-Compton scattered photons predicted to occur when the synchrotron brightness temperature exceeds a threshold value, usually estimated to be 10^{12} K. However, this effect appears to be in contradiction with observation because: (i) the threshold is substantially exceeded by several intra-day variable radio sources, but the inverse Compton emission is not observed, (ii) powerful, extra-galactic radio sources of known angular size do not appear to congregate close to the predicted maximum brightness temperature.
Aims. We re-examine the parameter space available to synchrotron sources using a non-standard electron distribution, in order to see whether the revised threshold temperature is consistent with the data.
Methods. We apply the theory of synchrotron radiation to a population of monoenergetic electrons. The electron distribution and the population of each generation of scattered photons are computed using spatially averaged equations. The results are formulated in terms of the electron Lorentz factors that characterise sources at the threshold temperature and sources in which the particle and magnetic field energy density are in equipartition.
Results. We confirm our previous finding that intrinsic brightness temperatures
K can occur without catastrophic cooling. We show that substantially higher temperatures cannot be achieved either in transitory solutions or in solutions that balance losses with a powerful acceleration mechanism. Depending on the observing frequency, we find strong cooling can set in at a range of threshold temperatures and the imposition of the additional constraint of equipartition between particle and magnetic field energy is not warranted by the data.
Conclusions. Postulating a monoenergetic electron distribution, which approximates one that is truncated below a certain Lorentz factor (
), alleviates several theoretical difficulties associated with the inverse Compton catastrophe, including anomalously high brightness temperatures and the apparent lack of clustering of powerful sources at 10^{12} K.
Key words: galaxies: active - galaxies: high-redshift - galaxies: jets
When interpreting the synchrotron spectrum of powerful extra-galactic radio
sources, it is usual to assume that the underlying electron distribution has a power-law form
.
For a homogeneous source, the spectrum peaks at the point
,
where the optical depth to synchrotron self-absorption
is of the order of unity, and
above this frequency the intensity falls off as
.
Because more electrons become effective at
absorbing the radiation as the frequency decreases, the optically
thick part of the spectrum is not of the Rayleigh-Jeans type, but has
instead
,
independent of the power-law index of the underlying
distribution, (provided q>1/3). Correspondingly, the brightness
temperature peaks at
,
falling off as
to
lower and as
to higher frequencies. In this case,
Kellermann & Pauliny-Toth (1969) found that the
ratio of the luminosity
carried away by inverse Compton
scattered photons to that of synchrotron photons
is
The rapid rise in total luminosity implied by Eq. (1) when exceeds is called the "inverse Compton catastrophe''. The implied energy requirement of a 1 Jy source at z=1is at K, and a prohibitive at K at 1 GHz. A source of 1 GHz photons with magnetic field G, boosts a radio photon into the X-ray band in a single inverse Compton scattering. At K this would imply a nanoJansky X-ray flux, not untypical of strong extragalactic sources. However, at K Eq. (1) predicts an X-ray flux at the milliJansky level, in contradiction with observation (Padovani et al. 2004; Sambruna et al. 2000; Fossati et al. 1998; Tavecchio et al. 2002; Guainazzi et al. 2006).
Despite this, several sources that display intra-day variability (IDV) in their radio emission have an implied brightness temperature that exceeds 10^{12} K by several orders of magnitude, if the observed variability is intrinsic (e.g., Kraus et al. 2003). Even if the variability is caused by scintillation, the implied brightness temperature can still greatly exceed 10^{12} K for some sources (Wagner & Witzel 1995). Currently, the most extreme example is the source PKS 0405-385. This source displays diffractive scintillation (Macquart & de Bruyn 2005), which places an upper limit on its angular size that corresponds to a brightness temperature of 2 10^{14} K. These sources are generally assumed to be beamed, i.e., to be in relativistic motion towards the observer (e.g., Rees 1966; Jones & Burbidge 1973; Singal & Gopal-Krishna 1985). In this case the intrinsic temperature is lower than that deduced for a stationary source by a factor of (for intrinsic variability) or (for scintillating sources). (Here the Doppler factor with the source velocity and the angle between this velocity and the line of sight.) Nevertheless, the observed brightness temperatures are too high to be accounted for by Doppler factors similar to those estimated from observations of superluminal motion (Cohen et al. 2003).
A second problem arises with powerful sources whose angular extent can be measured directly. In an analysis of high brightness temperature radio sources in which Doppler beaming is thought to be absent, Readhead (1994) measured a brightness distribution that cuts off at 10^{11} K; one order of magnitude lower than the inverse Compton limit. This appears consistent with observations of a sample of 48 sources showing superluminal motion (Cohen et al. 2003), in which it was found that the intrinsic brightness temperature cluster around 2 10^{10} K. Readhead (1994) argued that an apparent maximum brightness temperature significantly lower than 10^{12} K could not be caused by catastrophic Compton cooling. Instead, he suggested that sources are driven towards equipartition between their magnetic and particle energy contents. Assuming, in addition, that observations are taken at the peak of the synchrotron spectrum, and that the electron distribution is a power-law, he showed that the expected distribution of brightness temperatures was consistent with that observed.
In this paper, we re-examine these two problems assuming that the source contains a monoenergetic electron distribution instead of the conventional power-law. Although this assumption appears at first sight highly restrictive, the form of the synchrotron emissivity means that under some circumstances such a distribution provides a good approximation to several more commonly encountered cases, including that of a conventional power-law distribution that is truncated to lower energy at a Lorentz factor . Such distributions have been proposed in connection with radio sources for a variety of reasons: the absence of low energy electrons can account for the lack of Faraday depolarisation in parsec-scale emission regions (Wardle 1977; Jones & Odell 1977) and has recently been discussed in connection with statistical trends in the observed distribution of superluminal velocities as a function of observing frequency and redshift (Gopal-Krishna et al. 2004). Also, Blundell et al. (2006) recently examined the radio and X-ray emission from the lobe regions of a giant radio galaxies 6C 0905+3955, and deduced a low energy cutoff of the relativistic particles in the hotspots of .
In Sect. 2 we use standard theory to discuss the general properties of the synchrotron spectra emitted by a homogeneous source. A set of spatially averaged equations describing the evolution of the electron Lorentz factor and both the synchrotron and the associated inverse Compton scattered emission is presented in Sect. 3. Having identified in these equations the threshold for the inverse Compton catastrophe, we discuss the parameter space available to stationary solutions in Sect. 4. Here we confirm the results reported in Kirk & Tsang (2006), indicating that temperatures considerably in excess of 10^{12} K are permitted. We also show that in the case of resolved sources, the onset of catastrophic cooling occurs over a wide range of temperatures, consistent with the observed temperature range. Finally, we address in Sect. 5 the suggestions by Slysh (1992) that extremely high brightness temperatures can be achieved in nonstationary sources either by injecting electrons at high energy, or by balancing their cooling against a powerful acceleration mechanism. A summary of our conclusions is presented in Sect. 6.
We consider a homogeneous source region characterised by a single spatial scale R, that contains monoenergetic electrons and possibly positrons of Lorentz factor and number density immersed in a magnetic field B. Expressions for the synchrotron emissivity and absorption coefficients can be found in many excellent texts (e.g., Rybicki & Lightman (1979, Chap. 6), and Longair (1992, Chap. 18)) and are summarised in our notation in Appendix A.
For any given source there exists a frequency below which absorption is important. Since B and also define a characteristic synchrotron frequency (see Eq. (A.3)), the sources we consider can be divided into two categories: those with weak absorption in which and those with strong absorption . Note that this division is independent of the observing frequency, since it relates only to intrinsic source properties. The synchrotron spectra that emerge in these two cases are quite different, and are illustrated in Fig. 1. A feature they have in common is that the low energy spectrum has the Rayleigh-Jeans form , where is the specific intensity at frequency . This property contrasts with the dependence of at low frequencies of a source containing a power-law distribution of electrons. The reason is that a power-law distribution contains cold (low energy) electrons that contribute to the absorption at low frequencies.
Figure 1: The synchrotron spectra ( upper panel) and brightness temperatures ( lower panel) of sources with monoenergetic electrons in the case of strong (blue) and weak (red) absorption. The green curves show the optically thick ( , see Appendix A) and optically thin ( ) approximation. In the upper panel, is in arbitrary units, and in the lower, the brightness temperature is normalised to the energy of the electron. The abscissa x is the ratio of the frequency to the characteristic synchrotron frequency of the electrons . The blue (red) curves correspond to a source which has an optical depth of unity to synchrotron self-absorption at ( ). For ease of display, the upper panel compares sources with equal flux at high frequency, whereas the lower compares sources with equal flux at low frequency. | |
Open with DEXTER |
The brightness temperature, , where is Boltzmann's constant, is a function of frequency and is also illustrated in Fig. 1. At low frequency, it attains its maximum value roughly in "equilibrium'' with the electrons: , then decreases monotonically to higher frequencies. In the case of weak absorption, for , and then cuts off exponentially as once is exceeded. In strongly absorbed sources, the brightness temperature remains almost constant until the frequency exceeds upon which it falls off as until the source becomes optically thin, after which the exponential cut-off takes over.
Although four parameters (,
,
B and R)
are needed to define a source model, the division between strong and weak
absorption is simple. It occurs at a critical Lorentz factor
given by (see Eq. (A.11))
In his model of high-brightness temperature sources, Slysh (1992) considered the strong absorption case. The most important property of the assumed distribution in this case is the lack of high energy electrons: the addition of a population of cold electrons, which would correspond to a power-law distribution truncated to higher Lorentz factors, would reduce the brightness temperature of the source at x<1 (in Fig. 1) but would not significantly influence this quantity, for x>1.
On the other hand, Crusius-Waetzel (1991) and Protheroe (2003) considered weak absorption, where the key property of the model distribution is the absence of low energy electrons. In this case, the monoenergetic model is a good approximation to a power-law distribution truncated to lower electron energies at . The addition of a high-energy power-law tail affects the spectrum at x>1, but does not change the maximum brightness temperature achieved at . Furthermore, the truncation need not be sharp: provided the opacity at low frequencies is dominated by the contribution of electrons with , the monoenergetic approximation is good. This is the case if, for , the spectrum is sufficiently hard: with . In particular, the low energy tail of a relativistic Maxwellian distribution (q=-2) falls into this category.
In contrast to the pure power-law distribution, where the self-absorption turnover is strongly peaked, the emission of a weakly absorbed source - shown in red in the upper panel of Fig. 1 - is flat over nearly two decades in frequency. It therefore provides a natural explanation of compact flat-spectrum sources, eliminating the need to appeal to a "cosmic conspiracy'' behind the superposition of peaked spectra from different parts of an inhomogeneous source (Marscher 1980).
For the treatment of inverse Compton scattering, it is
necessary to evaluate the energy density in synchrotron photons in a given source. To do this,
must be integrated over angles and
over frequency. The result depends on the geometry and
optical depth as well as the position
within the source. An average value can be estimated by
introducing a geometry dependent factor, ,
defined according to:
(5) |
An approximate, spatially averaged set of equations governing the energy balance of particles and
synchrotron radiation in a source can be found following the
approach of Lightman & Zdziarski (1987) and Mastichiadis & Kirk (1995). In
terms of the time-dependent synchrotron radiation energy density U_{0}(t)
one can write:
The corresponding equation for the particles that takes into account
synchrotron absorption and emission as well as an acceleration term
takes the form
Multiple inverse Compton scatterings can be accounted for as follows: first we label the
photons present in the source according to how many scattering events they
have suffered after production by the synchrotron process.
The energy density of these photons is denoted by U_{i}Thus, i=0 corresponds to photons emitted by the synchrotron process which
have not undergone a scattering, and the corresponding energy density
is governed by Eq. (10).
Assuming the source is optically thin to Thomson
(or Compton) scattering, the dominant loss mechanism for
the energy density of photons
belonging to a given generation
is escape from the source, rather than
conversion to the i+1'th generation. In this case, we can write
for the time-dependence of U_{i}:
If the inverse scattering process proceeds in the Thomson regime
a simple expression can be found for Q_{i}. However, as i increases,
also increases, eventually becoming comparable to the electron
energy when viewed in its rest frame. When this
happens, Klein-Nishina modifications to the Thomson
cross section become important, reducing the value Q_{i}.
We take approximate account of this effect by limiting the number of
scatterings to
,
and using the Thomson approximation
to evaluate Q_{i} for
.
In this case, the average energy of a scattered photon of the i'th generation is
and the rate of such scatterings in unit volume of the source is
.
Therefore
(25) |
For comparison with observations of intra-day variable sources, it is
convenient to formulate the expression for the
specific intensity given in Eq. (A.15) in terms of
quantities accessible to observation. Expressing the result in terms of the
observed (at z=0) brightness temperature, we find, in the case of weak
absorption, and at low frequency (
)
According to Eq. (27), brightness temperatures of K, such as observed in the sources PKS 1519-273 and PKS 0405-385 (Macquart et al. 2000; Rickett et al. 2002) can be understood within a simplified homogeneous synchrotron model in which , implying a relatively modest inverse Compton luminosity, i.e., no catastrophe. Even the extremely compact source J 1819+3845, which has K can be accommodated in a catastrophe-free model provided the Doppler factor is greater than about 15. In each case, a hard spectrum is predicted, extending to Hz. Although the dependence of the brightness temperature on this parameter is quite weak, simultaneous observations in the radio to IR and optical (Ostorero et al. 2006) have the potential to rule out this explanation on a source by source basis.
In a seminal paper, Readhead (1994) discussed the distribution in brightness temperature of a sample of powerful sources whose angular size could either be measured directly, or constrained by interplanetary scintillation. In discussing these objects several simplifications must be made, even within the context of a homogeneous synchrotron model.
Firstly, in the two low frequency samples (81.5 MHz and 430 MHz) considered by Readhead (1994), the emission is thought to be almost isotropic. Doppler boosting is then unimportant and can be neglected. Secondly, these sources are not very compact; their extension on the sky is typically between 0.1 and 1 arcsec. Therefore, for our discussion we fix the linear extent R of the source to 1 kpc, corresponding to an angular size of approximately 0.2 arcsec at redshift z=1. This leaves three parameters needed to specify the source model: the magnetic field strength B, the electron density and the Lorentz factor of the electrons. In order to clarify the physics of a source, we transform from the parameter set to the characteristic Lorentz factors and defined in Eqs. (23) and (24). Our basic parameter set is therefore ). Finally, in order to display on a two-dimensional figure source properties such as brightness temperature and spectral slope at a particular frequency, we consider a slice through this three dimensional parameter space, selecting parameters such that the particle and magnetic energy densities are in equipartition: .
The properties of source models on this slice are shown in the - plane in Fig. 2. This plane can immediately be divided into regions of strong and weak absorption, as defined in Eq. (2). The boundary, drawn as a thick dashed line, represents the locus of the points at which . Weakly absorbed sources lie towards higher and (i.e., the upper-right side) and strongly absorbed sources towards lower and (i.e., the lower-left side). We also show (in white) contours of the magnetic field strength.
Figure 2: The brightness temperature as a function of and assuming equipartition between the magnetic and particle energy densities and a source size 1 kpc. Black contour lines indicate , 10, 11, 12, 13 and 14. The red dot-dashed line is the locus of points at which the characteristic synchrotron frequency of the emitting particles is 81.5 MHz, the yellow short dashed line shows where the source has an optical depth of unity at this frequency. The long dashed line divides regions of strong absorption (to the left) from those of weak absorption (to the right). The diagonal is shown as a dotted line. Contour lines of the magnetic field strength are shown in white, ranging from to 0 (in the bottom right-hand corner). | |
Open with DEXTER |
The remaining source properties depend upon the choice of observing frequency. In Fig. 2 we take this to be 81.5 MHz, corresponding to the low frequency sample discussed by Readhead (1994). In order to determine the spectral slope of a given source, we plot as a yellow short dashed line the locus of points where the observing frequency coincides with the frequency at which the optical depth to absorption is unity, . Sources that lie above this line (on the side of larger ) are optically thin at the chosen observing frequency. In addition, the red dot-dashed line in Fig. 2 gives the locus of points where the observing frequency equals the characteristic frequency of synchrotron radiation . By definition, the intersection point of these lines lies also on the boundary between weak and strong absorption (the long dashed line). The observing frequency lies below on the lower-right side of the dot-dashed line. The (colour) shading gives the intrinsic brightness temperature at the chosen observing frequency.
The two lines (yellow short dashed and red dot-dashed) divide the - -plane in Fig. 2 into four regions with differing spectral properties: in region A, sources have a Rayleigh-Jeans spectrum , in region B, the spectrum is that of low frequency, optically thin synchrotron radiation , in region C, it is close to (see Slysh 1992) and in region D it falls off exponentially . Consequently, flat spectrum sources reside in region B, preferentially close to the yellow short dashed line and in region C, preferentially close to the red dot-dashed line.
Sources that are in equipartition and lie below the threshold of the Compton catastrophe are to be found in the upper left half of Fig. 2, above the dotted line on which . The maximum brightness temperature accessible to these sources occurs close to , and is approximately 10^{12.6} K, in rough agreement with the results of Kellermann & Pauliny-Toth (1969), who, however, did not assume their sources to be in equipartition. The brightest sources are weakly absorbed, (they lie to the right of the long dashed line) and have a magnetic field strength of a few milliGauss. Their optical depth to synchrotron self-absorption lies close to unity at the observation frequency (they lie close to the yellow short dashed line).
Singal & Gopal-Krishna (1985) first discussed the effects of the additional assumption of equipartition on bright sources and used it to estimate Doppler factors for rapidly variable sources. Later, Readhead (1994) introduced the concept of an "equipartition brightness temperature'' to explain the observation that the temperature distribution of resolved sources appears to peak significantly below 10^{12} K. However, the crucial additional assumptions in his treatment is that the source flux is measured at the "synchrotron peak'', and that the electron distribution is a power-law in energy. This implies that the opacity at a given frequency (e.g., at the synchrotron peak) is dominated by those electrons with a corresponding characteristic frequency. In our model, in which the electron distribution is approximated as monoenergetic, these assumptions are roughly equivalent to demanding that the source lies on the red dot-dashed line in Fig. 2 if it is weakly absorbed (i.e., on the boundary of regions B and D), and on the yellow short dashed line if it is strongly absorbed (i.e., on the boundary of regions C and D). This leads to a maximum brightness temperature of a few times 10^{10} K, as found by Readhead (1994). Furthermore, as noted by Readhead (1994), such sources lie far from the threshold temperature, achieved along the dotted line in Fig. 2.
Replacing the assumption that the source flux is measured at the synchrotron peak, by the requirement that its spectrum be flat, i.e., that it lie in region B of Fig. 2, one sees that a wide range of brightness temperatures is available for sources in equipartition, extending up to the threshold temperature found by Kellermann & Pauliny-Toth (1969). Thus, the observed temperature distribution is not explained by the assumption of equipartition.
In order to explain the occurrence of brightness temperatures above 10^{12} K, Slysh (1992) formulated a model involving a monoenergetic electron distribution in a strongly absorbed source, in the sense that , where is defined in Eqs. (2) and (3). He considered two scenarios, (i) a time-dependent one in which electrons were injected at arbitrarily high Lorentz factors and allowed to cool and (ii) one in which a strong continuous re-acceleration of the electrons led to a high brightness temperature equilibrium.
In each case, the assumption that the source is strongly absorbed leads to extreme values of the parameters. For example, in the first scenario in which high energy particles are injected into the source, Slysh (1992) finds that a brightness temperature of 10^{15} K can be sustained over 1 day at an observing frequency of 1 GHz. This is clearly in conflict with our analysis. The electron Lorentz factor required to achieve this temperature is . However, the condition that the source is strongly absorbed, which is used in this model to estimate the cooling rate, combined with the condition GHz required for a flat spectrum, leads to an extremely large Thomson optical depth, , as well as an implausibly low magnetic field 10^{-11} G. The parameter that determines the inverse Compton luminosity is approximately 10^{12}, which implies an extremely large compactness of the inverse Compton radiation from the source. The resulting copious pair production invalidates the analysis and, ultimately, reduces the brightness temperature achievable in the radio range. The same criticism applies also to the second scenario described by Slysh (1992) in which acceleration balances inverse Compton losses to provide a brightness temperature of 10^{14} K at 1 GHz.
In the absence of Klein-Nishina effects on the scattering cross section, we find the time dependence of the particle and photon energies can be described by the three ordinary differential Eqs. (20)-(22). Inspection of these shows that if the threshold temperature is exceeded ( ), the inverse Compton luminosity grows in a timescale of roughly the light-crossing time of the source. Thus, the threshold can only be substantially exceeded if the acceleration process in Eq. (22) operates on a shorter timescale. However, these equations employ a spatial average over the emission region. Although a rapid acceleration rate might be achieved locally in small regions of the source, once an average is taken, no timescale in the system can be shorter than the light-crossing time of the region over which the accelerated particles are distributed. In this case, the threshold temperature cannot be significantly exceeded.
At first sight, Klein-Nishina effects offer a possible escape from this conclusion. If even the first order scattering is suppressed, which requires extremely large Lorentz factors for the electrons ( is needed for Klein-Nishina effects when scattering 10 GHz photons), the strong reduction in the rate of cooling by inverse Compton scattering suggests that higher brightness temperatures might be possible.
This is, however, not the case, because the rate of production of
electron-positron pairs by photon-photon interactions becomes
important. The strength of this effect, which is not included in our
model equations, can be measured in terms of the "compactness parameter''
(see, for
example Mastichiadis & Kirk 1995), defined as
(28) |
We illustrate this in Fig. 3, where we compare two models with the same linear size R (and observing frequency), but different electron densities and different values of B, chosen as follows: For any given set of parameters, R, B and , and observing frequency , the optical depth to synchrotron absorption , as defined in Eq. (A.12), has a single maximum as a function of , located close to the point where equals the characteristic synchrotron frequency. If the source is optically thick to absorption at this point, then , as described in Sect. 2, and the brightness temperature is roughly . If, on the other hand, the source is optically thin at this point, then , but the brightness temperature, given approximately by , decreases to higher , as can be seen from Eq. (A.12). Thus, assuming inverse Compton scattering does not intervene, the maximum brightness temperature is observed at a frequency such that , when , which implies . These conditions are imposed on the parameters of the models presented in Fig. 3. In addition to the source size, chosen to be R=0.01 pc and the observing frequency, set to 1 GHz, this leaves one free parameter, which we choose to be the optical depth to Thomson scattering .
Figure 3: Upper panel: the brightness temperature (black), and the compactness (gray) as functions of time, for two stationary, local sources ( , z=0) with linear size R=0.01 pc, observed at 1 GHz. The Thomson optical depth is (dashed lines) and (solid lines) and the remaining parameters are chosen such that the optical depth to synchrotron self-absorption at (see Eqs. (2) and (3)). A horizontal line is drawn to indicate . Lower panel: the electron Lorentz factor (black dashed) and the optical depth to synchrotron self-absorption (gray dashed) for the case . A horizontal line indicates . | |
Open with DEXTER |
The upper panel in Fig. 3 shows the time-dependence of the brightness temperature found by solving Eqs. (12) and (17) numerically for sources with (dashed black line) and (solid black line), without allowance for Doppler boosting ( ). These sources have and , respectively and, in the absence of inverse Compton cooling, they could potentially achieve brightness temperatures of K and K. In order to do so, rapid acceleration is required, since for these source parameters, inverse Compton cooling leads to a time-asymptotic value of the Lorentz factor that is somewhat lower than for slow acceleration. The exact value of the asymptotic solution depends on the strength of the acceleration. For acceleration on the light-crossing timescale, it corresponds to (see Eq. (22)). In Fig. 3 we choose , which leads to an overshoot that slightly exceeds .
For , the compactness, shown as a function of time by the gray dashed line, remains well below unity, so that the effects of pair production can be neglected. However, this is not the case for . Here, the compactness (solid gray line) rises rapidly, reaching unity at , where 10^{12} K, well below its potential maximum. Thus, the attempt to gain higher brightness temperature by increasing , and, hence, , leads to a breakdown in the model assumptions due to pair production.
The lower panel of Fig. 3 shows the electron Lorentz factor and the optical depth to synchrotron self-absorption as functions of time for the case . The Lorentz factor (black dashed line) overshoots both its time-asymptotic value and . Correspondingly, the optical depth, (shown as the gray dashed line) which initially rises with , reaching unity at goes through a maximum very shortly afterwards. However, the overshoot is not sufficient to push back below unity, and the maximum brightness temperature, which coincides with the maximum Lorentz factor, remains at 10^{12} K, somewhat below the value of K, estimated for large optical depth.
The well-known upper limit on the brightness temperature of a synchrotron source K imposed by the inverse Compton catastrophe, has been reassessed, assuming a monoenergetic electron distribution.
In weakly absorbed sources (see Eq. (2)), this distribution mimics the situation in which the conventional power-law is truncated to lower energies at a Lorentz factor . Using the standard theory of synchrotron emission and self-absorption, we find that, for such sources, the brightness temperature at a frequency of a few GHz can reach approximately 10^{14} K, the precise limit being given in Eq. (27). Physically, this increased limit reflects the absence of cool electrons in monoenergetic distributions and in those that are truncated or hard below a certain Lorentz factor. As a consequence, intra-day variable sources can in principle be understood without recourse to other mechanisms such as unusually large Doppler factors (Rees 1967), coherent emission (e.g., Begelman et al. 2005) or proton synchrotron radiation (Kardashev 2000).
The possibility of exceeding the new limit in a time-dependent solution by balancing losses against a strong acceleration term has been investigated using a set of spatially averaged equations. Provided the acceleration process remains causal i.e., the acceleration time averaged over the source remains longer than the light-crossing time, we find a modest overshoot is possible, but the maximum temperature is still restricted by Eq. (27). In strongly absorbed sources, such as those considered by Slysh (1992), high brightness temperatures cannot be attained in a self-consistent model of the kind we discuss. The underlying reason is that extremely compact sources would be required, in which copious pair-production must be taken into account.
We have examined in detail the parameter space available to homogeneous synchrotron sources of fixed size. In the case of flat spectrum sources, we find that the imposition of the condition of equipartition between the particle and magnetic field energy densities does not result in a lower limit on the brightness temperature than that given by the inverse Compton catastrophe. Suggestions to the contrary (Readhead 1994) are based on the more restrictive twin assumptions that the power-law electron distribution is not truncated within the relevant range, and that the temperature is measured at the point where the optical depth of the source is approximately unity. Consequently, the observed temperature distribution does not support the equipartition hypothesis. We also find that flat spectrum sources close to equipartition can approach the threshold temperature of the inverse Compton catastrophe, in contrast with the finding based on the more restrictive assumptions in Readhead (1994).
A corollary of the theory presented here is that for very bright sources, which are necessarily weakly absorbed, a high degree of intrinsic circular polarisation is predicted (Kirk & Tsang 2006). In addition, the theory can be tested by comparison of the predicted synchrotron spectrum with simultaneous observations of high brightness temperature sources in the radio to infra-red range and comparison of the predicted inverse Compton emission with measurements at MeV to GeV energies.
We consider a region of homogeneous magnetic field B, linear dimension R, (and volume R^{3}) containing monoenergetic electrons/positron of number density
and Lorentz factor .
The volume emissivity for synchrotron radiation, summed over polarisations, is:
(A.10) |
(A.13) |
(A.14) |
= | (A.17) | ||
(A.18) |
To find the energy density
in synchrotron photons in a given source,
must be integrated over angles and over frequency. The result depends on the geometry and optical depth as well as the position within the source. However, an average value can be estimated by introducing a geometry dependent factor
:
(A.20) |
(A.21) |
(A.23) |