Synchrotron Sources

Minimum Energy and Equipartition

The existence of a synchrotron source implies the presence of relativistic electrons with some energy density $U_{\rm e}$ and a magnetic field whose energy density is $U_{\rm B} = B^2 / (8 \pi)$.  What is the minimum energy in relativistic particles and magnetic fields required to produce a synchrotron source of a given luminosity?

To estimate $U_{\rm e}$, we assume a power-law electron energy distribution
$$ N(E) \approx K E^{-\delta} $$ spanning the energy range $E_{\rm min}$ to $E_{\rm max}$ that is needed to produce synchrotron radiation over any frequency range $\nu_{\rm min}$ to $\nu_{\rm max}$. Then $$U_{\rm e} = \int_{E_{\rm min}}^{E_{\rm max}} E N(E)\, d E$$ For a given synchrotron luminosity
$$L = \int_{\nu_{\rm min}}^{\nu_{\rm max}} L_\nu d \nu ,$$ $$ {U_{\rm e} \over L} \propto { \int_{E_{\rm min}}^{E_{\rm max}} E N(E)\, dE \over -\int_{E_{\rm min}}^{E_{\rm max}} (dE/dt) N(E)\, dE}$$ Substituting $N(E) = K E^{-\delta}$ and the synchrotron power emitted per electron $(-dE/dt) \propto B^2 E^2$ gives
$${U_{\rm e} \over L} \propto { K \int_{E_{\rm min}}^{E_{\rm max}} E^{1 - \delta} \, dE \over K B^2 \int_{E_{\rm min}}^{E_{\rm max}} E^{2-\delta} \, dE}$$ $${U_{\rm e} \over L} \propto {E^{2 - \delta} \vert_{E_{\min}}^{E_{\rm max}} \over B^2 E^{3 - \delta} \vert_{E_{\rm min}}^{E_{\rm max}} }$$

Since electrons with energy $E$ emit most of the radiation seen at frequency $\nu \propto E^2 B$, the electron energy needed to produce radiation at frequency $\nu$ scales as
$$E \propto B^{-1/2}$$ If we consider the energy content of only those electrons that emit in a fixed frequency range (e.g., from $\nu_{\rm min} \sim 10^7$ Hz to $\nu_{\rm max} \sim 10^{11}$ Hz), then the energy limits $E_{\rm min}$ and $E_{\rm max}$ are both proportional to $B^{-1/2}$ and
$${U_{\rm e} \over L } \propto { (B^{-1/2})^{2 - \delta} \over B^2 (B^{-1/2})^{3 - \delta} } = {B^{-1 + \delta/2} \over B^2 B^{-3/2 + \delta/2}} = B^{-3/2}$$ We conclude that
$$\bbox[border:3px blue solid,7pt]{U_{\rm e} \propto B^{-3/2}}\rlap{\quad \rm {(5D1)}} $$
and we already know that $$U_{\rm B} \propto B^2~.$$

The "invisible" cosmic-ray protons and heavier ions emit negligible synchrotron power but they still contribute to the total cosmic-ray particle energy. If we call the ion/electron energy ratio $\eta$, then the total energy density in cosmic rays is $(1 + \eta) U_{\rm e}$.  The total energy density $U$ of both cosmic rays and magnetic fields is
$$\bbox[border:3px blue solid,7pt]{U = (1 + \eta) U_{\rm e} + U_{\rm B}}\rlap{\quad \rm {(5D2)}}$$
We cannot measure $\eta$ directly in distant radio sources, but cosmic rays collected near the Earth have $\eta \approx 40$. 

The greatly differing dependences of $U_{\rm e}$ and $U_{\rm B}$ on $B$ means that the total (cosmic ray plus magnetic) energy density $U(B)$ has a fairly sharp minimum near the point at which $(1 + \eta) U_{\rm e} \approx U_{\rm B}$.



For a source of a given synchrotron luminosity, the particle energy density $U_{\rm E} \equiv (1 + \eta)U_{\rm e}$ is proportional to $B^{-3/2}$ and the magnetic energy density $U_{\rm B}$ is proportional to $B^2$.  The total energy density $U = U_{\rm E} + U_{\rm B}$ has a fairly sharp minimum near equipartition of the particle and magnetic energy densities ($U_{\rm E} \approx U_{\rm B}$).

The minimum of the total energy density $U$ occurs at
$${d U \over d B} = {d [(1 + \eta) U_{\rm e} + U_{\rm B}] \over d B} = 0$$ First we evaluate the derivative of the electron energy density $U_{\rm e} \propto B^{-3/2}$.
$${d U_{\rm e} \over d B} \cdot U_{\rm e}^{-1} = -\biggl( {3 \over 2} \biggr) B^{-5/2} B^{3 / 2} = -{3 \over 2 B}$$ so
$${dU_{\rm e} \over d B} = - {3 U_{\rm e} \over 2 B}$$ Next we evaluate the derivative of the magnetic-field energy density $U_{\rm B} \propto B^2$.
$${d U_{\rm B} \over d B} \cdot U_{\rm B}^{-1} = { 2 B \over B^2 } = {2 \over B}$$ so
$${dU_{\rm B} \over d B} = { 2 U_{\rm B} \over B}$$ Inserting these results into the minimum-energy equation gives
$${d [(1 + \eta) U_{\rm e}] \over d B} + {d U_{\rm B} \over d B} = 0 = -{ 3 (1 + \eta) U_{\rm e} \over 2 B} + {2 U_{\rm B} \over B}$$

At minimum energy, the ratio of cosmic-ray particle energy to field energy is
$$\bbox[border:3px blue solid,7pt]{{{\rm particle~energy} \over {\rm field~energy}} = {(1 + \eta) U_{\rm e} \over U_{\rm B}} = {4 \over 3}}\rlap{\quad \rm {(5D3)}}$$ This ratio is nearly unity. Thus minimum energy implies (near) equipartition of energy: the total cosmic-ray energy density (including the energy of the nonradiating ions) $(1 + \eta) U_{\rm e}$ is nearly equal to the total magnetic energy density $U_{\rm B}$. We don't really know if equipartition exists in most sources, but radio astronomers often assume so, for several reasons:

(1) It is physically plausible—systems with interacting components often tend toward equipartition.

(2) Large and luminous extragalactic radio sources such as Cyg A have enormous energy requirements even near equipartition; the problem of explaining the large energy is even worse otherwise.

(3) It eliminates an unknown parameter and permits estimates of the relativistic particle energies and the magnetic field strengths of radio sources.

Getting the actual numerical values of the particle and magnetic-field energies from the synchrotron emission coefficient is a straightforward by tedious algebraic chore (see Rohlfs & Wilson Section 9.10).  Below are the results (from Pacholczyk's Radio Astrophysics, p. 171).  The functions $c_{12}$ and $c_{13}$ in these equations absorb the integrations from frequency $\nu_{\rm min}$ to $\nu_{\rm max}$ and the physical constants in Gaussian cgs units.


For a spherical radio source with radius $R$ and magnetic field strength $B$, the total magnetic energy is
$$E_{\rm B} = U_{\rm B} V = {B^2 \over 8 \pi} {4 \pi R^3 \over 3} = {B^2 R^3 \over 6}~.$$ The minimum-energy magnetic field is
$$\bbox[border:3px blue solid,7pt]{B_{\min} = [4.5 (1 +\eta) c_{12} L]^{2/7} R^{-6/7} {\rm ~Gauss}}\rlap{\quad \rm {(5D4)}}$$ and the corresponding energy in relativistic particles is
$$\bbox[border:3px blue solid,7pt]{E_{\rm min}{\rm (total)} = c_{13} [(1+ \eta) L]^{4/7} R^{9/7} {\rm ~ergs}}\rlap{\quad \rm {(5D5)}}$$ where the radio luminosity $L$ is conventionally integrated over the observable frequency range conventionally taken to be $\nu_{\rm min} = 10^7$ Hz to $\nu_{\rm max}= 10^{11}$ Hz,
$$L = \int_{\nu_{\rm min}}^{\nu_{\rm max}} L_\nu d \nu~, \qquad L_\nu = 4 \pi D^2 S_\nu~,$$ $D$ is the source distance, $S_\nu$ is its flux density at frequency $\nu$, and
$$ 1 + \eta \equiv {{\rm energy~in~all~relativistic~particles} \over {\rm energy~in~relativistic~electrons}}~.$$

The minimum total energy in relativistic particles and fields occurs when
$$ {(1+\eta) U_{\rm e} \over U_{\rm B}} \approx {4 \over 1 + \delta} ,$$ where
$$\delta = 2 \alpha + 1 \sim 2.4~.$$

The synchrotron lifetime of a source is defined as the ratio of electron energy to the energy loss rate $L$ from synchrotron radiation:
$$\tau \equiv {E_{\rm e} \over L}~. $$ $$\bbox[border:3px blue solid,7pt]{\tau \approx c_{12} B_\bot^{-3/2}}\rlap{\quad \rm {(5D6)}}$$

The functions $c_{12}$ and $c_{13}$ in Gaussian cgs units are plotted below.



Plots of $c_{12}$ in Gaussian cgs units as a function of [negative] spectral index $\alpha \equiv - d \log S / d  \log \nu$ for $\nu_{\rm min} = 10^6$ Hz (dashed curves) and $10^7$ Hz (solid curves) and $\nu_{\rm max} = 10^{10}$ Hz, $10^{11}$ Hz, and $10^{12}$ Hz.



Plots of $c_{13}$ in Gaussian cgs units as a function of [negative] spectral index $\alpha \equiv - d \log S / d  \log \nu$ for $\nu_{\rm min} = 10^6$ Hz (dashed curves) and $10^7$ Hz (solid curves) and $\nu_{\rm max} = 10^{10}$ Hz, $10^{11}$ Hz, and $10^{12}$ Hz.

Example: What are the minimum-energy magnetic field strength and the minimum total energy of Cygnus A, a luminous double radio source (see the VLA image below) at distance $D \approx 230$ Mpc (for $H_0 = 75$ km s$^{-1}$ Mpc$^{-1}$)? The lobe radii are $R \approx 30$ kpc and the total flux density of Cyg A is
$$S_\nu \approx 2000 {\rm ~Jy~} \biggl( {\nu \over {\rm GHz} } \biggr)^{-0.8}$$

VLA image of Cyg A

First we convert the data from "astronomical" units to Gaussian cgs units:
$$R = 30 {\rm ~kpc} \times {10^3 {\rm ~pc} \over {\rm kpc}} \times {3.09 \times 10^{18} {\rm ~cm} \over {\rm pc}} \approx 9.0 \times 10^{22} {\rm ~cm}$$ $$S_\nu = 2000 {\rm ~Jy} \biggl( {10^{-23} {\rm ~erg~s}^{-1} {\rm ~Hz}^{-1} {\rm ~cm}^{-2} \over {\rm Jy}} \biggr) \biggl( {\nu \over 10^9 {\rm ~Hz}}\biggr)^{-0.8}$$ $$S_\nu = {3.17 \times 10^{-13} {\rm ~erg~s}^{-1} {\rm ~Hz}^{-1} \over {\rm cm}^2} \times\biggl( {\nu \over {\rm Hz}} \biggr)^{-0.8}$$ $$D = 230 {\rm ~Mpc} \times {10^6 {\rm ~pc} \over {\rm Mpc}} \times { 3.09 \times 10^{18} {\rm ~cm} \over {\rm pc}} \approx 7.1 \times 10^{26} {\rm ~cm}$$

The spectral luminosity of Cyg A is
$$L_\nu \approx 4 \pi D^2 S_\nu = 4 \pi (7.1 \times 10^{26} {\rm ~cm})^2 \times {3.17 \times 10^{-13} {\rm ~erg~s}^{-1} {\rm ~Hz}^{-1} \over {\rm cm}^2} \times \biggl({\nu \over {\rm Hz}}\biggr)^{-0.8}$$ $$L_\nu \approx 2.0 \times 10^{42} {\rm ~erg~s}^{-1} {\rm ~Hz}^{-1} \times \biggl({\nu \over {\rm Hz}}\biggr)^{-0.8}$$

The total radio luminosity of Cyg A in the frequency range $10^7$ Hz to $10^{11}$ Hz is
$$L = \int_{10^7 {\rm ~Hz}}^{10^{11} {\rm ~Hz}} L_\nu d \nu$$ $$L \approx 2.0 \times 10^{42} {\rm ~erg~s}^{-1} {\rm ~Hz}^{-1} \, \biggl( { \nu^{0.2} \over 0.2} \biggr) \bigg|_{10^7 {\rm ~Hz}}^{10^{11} {\rm ~Hz}}$$ $$L \approx 2.0 \times 10^{42} {\rm ~erg~s}^{-1} \, \biggl[ { (10^{11})^{0.2} - (10^7)^{0.2} \over 0.2} \biggr]$$ $$L \approx 1.33 \times 10^{45} {\rm ~erg~s}^{-1}$$

In units of the bolometric solar luminosity $L_\odot \approx 3.83 \times 10^{33}$ erg s$^{-1}$, the radio luminosity of Cyg A is
$${L \over L_\odot} \approx {1.33 \times 10^{45} {\rm ~erg~s}^{-1} \over 3.83 \times 10^{33} {\rm ~erg~s}^{-1}} \approx 3.5 \times 10^{11}$$ Thus the radio power from Cyg A exceeds the bolometric output from a galaxy of stars similar to our Galaxy.  The energy appears to originate in a compact object at the center of the host galaxy.  How massive must this compact object be to produce such a luminous source?


Eddington Limit

What is the maximum luminosity of an astronomical object of total mass $M$ in a steady state? The outward radiation pressure cannot exceed gravity. Otherwise, radiation pressure would expel the outer layers of a star in the form of a wind, or accretion onto a compact object would be disrupted. Even if the atmosphere or infalling material is ionized hydrogen, the free electrons will Thomson-scatter outflowing radiation. Each electron being blown away by radiation pressure will drag along one proton ($m_{\rm p} \gg m_{\rm e}$) to maintain charge neutrality. Balancing the radiation and gravitational forces on each electron/proton pair at distance $r$ from the accreting object gives the Eddington Luminosity:
$${L_{\rm E} \over 4 \pi r^2 c} \sigma_{\rm T} = {GM (m_{\rm p} + m_{\rm e}) \over r^2} \approx {GM m_{\rm p} \over r^2}$$ Note that the distance $r$ drops out and
$$L_{\rm E} \approx { 4 \pi G M m_{\rm p} c \over \sigma_{\rm T}}$$ depends only on the mass $M$ of the source.  In cgs units
$$L_{\rm E} ({\rm erg~s}^{-1}) = {4 \pi \times 6.67 \times 10^{-8} {\rm ~dyne~cm}^2 {\rm ~g}^{-2} \times M \times 1.66 \times 10^{-24} {\rm ~g} \times 3 \times 10^{10} {\rm ~cm~s}^{-1} \over 6.65 \times 10^{-25} {\rm ~cm}^2}$$ $$L_{\rm E} ({\rm erg~s}^{-1}) = 6.28 \times 10^4\, M ({\rm g})$$ Normalized to "solar" units $L_\odot \approx 3.83 \times 10^{33}$ erg s$^{-1}$ and $M_\odot \approx 1.99 \times 10^{33}$ g,
$$ \biggl( { L_{\rm E} \over L_\odot} \biggr) \approx {6.28 \times 10^4 \times 1.99 \times 10^{33}{\rm ~g} \over 3.83 \times 10^{33} {\rm ~erg~s}^{-1} } \biggl( { M \over M_\odot} \biggr)$$ $$\bbox[border:3px blue solid,7pt]{\biggl( { L_{\rm E} \over L_\odot } \biggr) \approx 3.3 \times 10^4 \biggl( { M \over M_\odot} \biggr)}\rlap{\quad \rm {(5D7)}}$$

As the mass of a main-sequence star approaches $M \approx 100 M_\odot$, its luminosity approaches its Eddington luminosity. Very massive stars often have radiation-driven winds, and stable stars more massive than $100 M_\odot$ may not be possible.


Example: What does the Eddington limit give for the minimum mass for a source of luminosity $L \approx 3.5 \times 10^{11} L_\odot$? We make the assumption that the average luminosity of the central mass was at least this much at some times, so
$$\biggl( { M \over M_\odot} \biggl) \geq {3.5 \times 10^{11} \over 3.3 \times 10^4} \approx 10^7$$ Note that the Eddington mass limit depends only on the instantaneous power emitted by the source, not on the total energy of the source, the source age, or any other indicator of its history.



Temperatures of Eddington-limited accretion disks near black holes

The Eddington limit is directly applicable to quasars and other sources having having luminous accretion disks surrounding supermassive black holes. If the gravitational energy released by accretion is thermalized, the hottest and brightest material will be concentrated just outside the innermost stable orbit surrounding the central rotating black hole. The radius of this orbit is
$$r = 3 r_{\rm g} = 3 \times {2 G M \over c^2}~,$$ where $r_g = 2 G M / c^2$ is the gravitational radius or Schwarzschild radius. If the compact object is accreting enough matter to approach its Eddington luminosity, the combination of luminosity and radius determines the blackbody temperature $T$ of the inner accretion disk.
$$L \approx 4 \pi r^2 \sigma T^4 \approx L_{\rm E}$$ $$ 4 \pi \biggl( {6 G M \over c^2} \biggr)^2 \sigma T^4 \approx 4 \pi \biggl( {G M m_{\rm p} c \over \sigma_{\rm T} } \biggr)$$ $$T^4 \approx \biggl( {m_{\rm p} c^5 \over 36 G \sigma \sigma_{\rm T}} \biggr) M^{-1}$$ The more massive the black hole, the cooler.  Inserting cgs values for the constants gives
$$T^4 \approx \biggl( {1.66 \times 10^{-24} { \rm ~g}\, (3 \times 10^{10} {\rm ~cm~s}^{-1})^5 \over 36 \times 6.67 \times 10^{-8} {\rm ~erg~cm}^{-2}\, {\rm s}^{-1}\, {\rm K}^{-4} \times 6.65 \times 10^{-25} {\rm ~cm}^2} \biggr) M^{-1}$$ $$T^4 \approx 4.46 \times 10^{62} {\rm ~g~K}^4 \, M^{-1}$$ In units of $M_{\odot} \approx 1.99 \times 10^{33}$ g,
$$\biggl( {T \over {\rm K}} \biggr) \approx 2.2 \times 10^7 \biggl( {M_\odot \over M} \biggr)^{1/4}$$


Example: The spectrum of the quasar 3C 273 is the superposition of a power law from synchrotron radiation and a thermal "big blue bump" that peaks at ultraviolet wavelengths.

big blue bump of 3C 273

The spectrum of 3C273 (Malkan, M. A., & Sargent, W. L. W. 1982, ApJ, 254, 22).


If $M \sim 10^9 M_\odot$ then $T \sim 10^5$ K, in agreement with the observed thermal spectrum and accounting for the strong emission-line spectrum of ionized hydrogen. Many quasars have similar blue bumps and appear to be accreting at rates approaching the Eddington limit. Note that black holes with masses of only a few $M_\odot$ accreting at the Eddington limit will have much higher temperatures $T \sim 10^7$ K and be strong thermal X-ray sources.



Returning to Cyg A, we estimate the magnetic field strength $B_{\rm min}$ that minimizes the total energy in relativistic particles and magnetic fields. We approximate Cyg A by two spherical lobes of radius $R \approx 30$ kpc and luminosity $L/2$, where $L$ is the luminosity of the whole source.
$$B_{\rm min} \approx [4.5 (1 + \eta) c_{12} (L/2)]^{2/7} R^{-6/7}$$ $$B_{\rm min} \approx (4.5 \times 3.9 \times 10^7 \times 1.33 \times 10^{45} {\rm ~erg~s}^{-1} / 2)^{2/7} (9 \times 10^{22} {\rm ~cm} )^{-6/7} (1 +\eta)^{2/7}$$ The value of $\eta$ is poorly constrained in extragalactic radio sources such as Cyg A. The cosmic rays accelerated by a supermassive black hole might be primarily electrons and positrons. Electrons and positrons have the same mass and charge (except for sign), so they are equally efficient at emitting synchrotron radiation and $\eta \approx 1$. If electrons and protons are acclerated to the same velocities (same $\gamma$), then the protons carry $m_{\rm p} / m_{\rm e} \sim 2 \times 10^3$ as much energy but emit almost nothing and $\eta \sim 2 \times 10^3$. Fortunately, $B_{\rm min} \propto \eta^{2/7}$ is only weakly dependent on $\eta$. Varying $\eta$ from 1 to $2 \times 10^3$ only changes $(1 +\eta)^{2/7}$ from about 1 to 9.
$$B_{\rm min} \approx 1.45 \times 10^{15} \times 2.1 \times 10^{-20} \times { \rm (1~to~9)~Gauss}$$ $$B_{\rm min} \approx ({\rm 30~to~300}) \times 10^{-6} \sim 10^{-4}{\rm ~Gauss}$$

The minimum total energy of Cyg A is twice the energy of each lobe: $$E_{\rm min} \approx 2 {\rm (lobes)} \times c_{13} [(1+\eta) L]^{4/7} R^{9/7}$$ $$E_{\rm min} \approx 2 \times 2.0 \times 10^4 \biggl( {1.33 \times 10^{45} {\rm ~erg~s}^{-1 } \over 2} \biggr)^{4/7} (9 \times 10^{22} {\rm ~cm})^{9/7} \times (1+\eta)^{4/7}~,$$ where $(1+\eta)^{4/7}$ is in the range of about 1 to 80.
$$E_{\rm min} \approx 4 \times 10^4 \times 4.1 \times 10^{25} \times 3.26 \times 10^{29} \times ({\rm 1~to~80}) {\rm ~ergs}$$ $$E_{\rm min} \approx 5.4 \times 10^{59} \times ({\rm 1 ~to~80}) {\rm ~ergs} \sim 5 \times 10^{60} {\rm ~ergs}$$

This enormous energy implies an independent lower limit to the mass of the central object powering the radio source. If mass could be converted to energy with 100% efficiency, the minimum mass needed to produce $E_{\rm min}$ would be
$$M \geq {E_{\rm min} \over c^2} \approx {5 \times 10^{60} {\rm ~ergs} \over (3 \times 10^{10} {\rm ~cm~s}^{-1})^2 } \approx 6 \times 10^{39} {\rm ~g}$$ $$M \geq 6 \times 1 0^{39} {\rm ~g} \biggl( { M_\odot \over 1.99 \times 10^{33} {\rm ~g}} \biggr) \approx 3 \times 10^6 M_\odot$$ This is a very conservative lower limit. Nuclear fusion can only convert mass to energy with about 1% efficiency, so $M > 3 \times 10^8 M_\odot$ if the energy source were nuclear fusion. Accretion onto a spinning black hole can result in efficiencies up to $(1 - 3^{-1/2}) \approx 0.4$ in theory, so it is consistent with $M > 10^7 M_\odot$. Many authors assume that mass is converted to energy with about 10% efficiency; this yields $M > 3 \times 10^7 M_\odot$. The small size of the radio core implied by Very Long Baseline Interferometery (VLBI) and variability of the core flux on time scales of months to years combined with the large minimum masses estimated from the Eddington limit and the total energy of the radio lobes together imply that the compact, massive object powering the radio source is a supermassive black hole. The adjective supermassive is used for black holes much more massive than the most massive stars, $\sim 100 M_\odot$.

A lower limit to the age of the radio source Cyg A is the synchrotron lifetime of the relativistic electrons estimated by taking the ratio of the electron energy to the observed synchrotron luminosity:
$$\tau \geq {E_{\rm min}/(1+\eta) \over L} \approx { 5.4 \times 10^{59} {\rm ~erg} \,(1+\eta)^{4/7} \over 1.33 \times 10^{45} {\rm ~erg~s}^{-1} \,(1+\eta) } \approx 4 \times 10^{14} {\rm ~s} \times \eta^{-3/7} \sim 10^{14} {\rm ~s} \sim 3 \times 10^6 {\rm ~yr}$$

Since each electron radiates energy at a rate proportional to $E^2$ and the critical frequency is proportional to $E^2$, the most energetic electrons emitting at the highest frequencies have the shortest lifetimes. The rapid depletion of high-energy electrons causes the emitted radio spectrum to steepen at high frequencies.


Baars et al. cyg A spectrum plot

The radio spectrum of Cyg A (and Cas A, Vir A) from Baars, J. W. M. et al. 1977, A&A, 61, 99.  Note the spectral steepening above $\nu \sim 10^3$ MHz.


Suppose that new relativistic electrons are continuously injected with a power-law energy distribution $N(E) \propto E^{-\delta_0}$ into a radio source. After a long time, electrons emitting at frequencies higher than $\nu$ will have been depleted by radiative losses $\propto E^2$, so these high-energy electrons will eventually have an energy distribution $N(E) \propto E^{-(\delta_0 + 1)}$. Consequently, the (negative) spectral index will be $\alpha_0 = (\delta_0 - 1) / 2$ at low frequencies and approach $\alpha = (\delta_0 + 1 - 1) / 2 = (\alpha_0 + 1/2)$ at higher frequencies; that is, the high-frequency spectrum steepens by $\Delta \alpha = 1/2$.

If the observed cutoff frequency $\nu$ is very high, the synchrotron lifetime of electrons with $\nu_{\rm c} \sim \nu$ may be less than the time needed for new relativistic electrons to travel from the core to the emitting feature in a jet or lobe. This implies in situ acceleration—something outside the core (e.g., shocks in the jet) must replenish the supply of relativistic electrons.


M87 jet from HST


Optical synchrotron emission in the radio jet of Virgo A = M87. Image credit



Synchrotron Self-Absorption

For every emission process there is an associated absorption process. The emitting particles in a source in local thermodynamic equilibrium (LTE) have a Maxwellian energy distribution, and such a source is called a thermal source. If the temperature of the particles is $T$, the source cannot have a brightness temperature greater than $T$. If the energy distribution of relativistic electrons in a synchrotron source were a (relativistic) Maxwellian, those electrons would have a characteristic temperature $T \sim E / 3 k$, and synchrotron self-absorption would prevent the brightness temperature from exceeding $T$. Astrophysical synchrotron sources are often called nonthermal sources because the energy distribution of the relativistic electrons is a power law and there is no single electron temperature $T$. However, self-absorption occurs regardless of the energy distribution.

In the approximation that electrons with energy $E = \gamma m_{\rm e} c^2$ in a magnetic field of strength $B$ emit only at the critical frequency
$$\nu_{\rm c} \sim {\gamma^2 e B \over 2 \pi m _{\rm e} c}~,$$ the Lorentz factor $\gamma$ of electrons emitting at frequency $ \nu$ is $$\gamma \approx \biggl( {2 \pi m_{\rm e} c \nu \over e B} \biggr)^{1/2 }~.$$ Since only electrons of one particular energy contribute to the emission and absorption at any one frequency in that approximation, the other electrons could have a relativistic Maxwellian energy distribution to match without changing the resulting emission and absorption at that frequency. Thus we expect that a sufficiently bright synchrotron source will be optically thick, and the brightness temperature at any frequency cannot exceed the effective temperature of those electrons emitting at that frequency.

In an ultrarelativistic gas, the ratio of specific heats at constant pressure and at constant volume is $c_{\rm p} / c_{\rm v} = 4/3$, not the nonrelativistic $5/3$, so the relation between electron energy $E$ and temperature $T_{\rm e}$ is
$$E = 3 k T_{\rm e} {\rm ~,~not~} (3/2) k T_{\rm e}~.$$ Thus the effective temperature of a relativistic electron is
$$T_{\rm e} = { E \over 3 k} = {\gamma m_{\rm e} c^2 \over 3 k}$$ Eliminating $\gamma$ in favor of $\nu$ gives the effective temperature of those electrons accounting for most of the radiation at frequency $\nu$: $$T_{\rm e} \approx \biggl( { 2 \pi m_{\rm e} c \nu \over e B} \biggr)^{1/2} {m_{\rm e} c^2 \over 3 k }$$ Numerically,
$$\bbox[border:3px blue solid,7pt]{\biggl( { T_{\rm e} \over {\rm K} } \biggr) \approx 1.18 \times 10^6 \biggl( { \nu \over {\rm Hz}} \biggr)^{1/2} \biggl( { B \over {\rm Gauss }} \biggr)^{-1/2}}\rlap{\quad \rm {(5D8)}}$$


Example: What is the effective temperature of the relativistic electrons emitting synchrotron radiation at $\nu = 0.1 {\rm ~GHz} = 10^8$ Hz if $B = 100\,\mu$Gauss $= 10^{-4}$ Gauss? $$\biggl( {T_{\rm e} \over {\rm K}} \biggr) \approx 1.18 \times 10^6 \times (10^8)^{1/2} (10^{-4})^{-1/2} \approx 10^{12}$$


At a sufficiently low frequency $\nu$, the brightness temperature $T_{\rm b}$ of any synchrotron source will approach the effective electron temperature $T_{\rm e}$ at that frequency and the source will become opaque. Starting with the definition of $T_{\rm b}$: $$I_\nu = {2 k T_{\rm b} \nu^2 \over c^2}$$ and setting $T_{\rm b} \approx T_{\rm e}$ gives
$$I_\nu \approx { 2 k T_{\rm e} \nu^2 \over c^2} \propto \nu^{1/2} \nu^2 B^{-1/2}$$ Then, because flux density is proportional to $I _\nu$ for a source subtending a given solid angle $\Omega$, the spectrum of a synchrotron self-absorbed and spatially homogeneous source is a power law of slope $5/2$:
$$\bbox[border:3px blue solid,7pt]{S(\nu) \propto \nu^{5/2}}\rlap{\quad \rm {(5D9)}}$$ independent of the slope $\delta$ of the electron-energy spectrum. The flux density of an opaque but truly thermal source (e.g., an HII region) is proportional to $\nu^2$; the extra $\nu^{1/2}$ for synchrotron radiation comes from the fact that $T_{\rm e} \propto \nu^{1/2}$.



The spectrum of a homogeneous cylindrical synchrotron source is a power law with slope $-\alpha = 5/2$ at low frequencies where $\tau \gg 1$.  Astrophysical sources are inhomogeneous, so their actual low-frequency spectral slopes are smaller than 5/2 and their spectral peaks are broader.  The optical depth of the source is $\tau = 1$ at $\nu = \nu_1$.


We can invert Equation 5D8 and use $T_{\rm e}\approx T_{\rm b}$ to estimate the magnetic field strength in a self-absorbed source whose brightness temperature has been measured.
$$\bbox[border:3px blue solid,7pt]{\biggl( { B \over {\rm Gauss}} \biggr) \approx 1.4 \times 10^{12} \biggl( {\nu \over {\rm Hz}} \biggr) \biggl( { T_{\rm b} \over {\rm K} } \biggr)^{-2}}\rlap{\quad \rm {(5D10)}}$$


Example: If a self-absorbed radio source is observed to have $T_{\rm b} \approx 10^{11}$ K at $\nu = 1$ GHz, the magnetic field strength is
$$ \biggl( { B \over {\rm Gauss}} \biggr) \approx 1.4 \times 10^{12} \times 10^9 \times (10^{ 11})^{-2} \approx 0.1$$


The spectra of "real" radio sources reflect the idealized spectra of uniform sources, but they are more complex because real sources have nonuniform magnetic fields and electron energy distributions in geometrically complex structures. Representative spectra of powerful radio galaxies and quasars are illustrated below.


Spectra of 3C sources GEGRA2, Fig 13.3
Spectra of radio galaxies and quasars.  The radio source 3C 84 in the nearby galaxy NGC 1275 contains a very compact nuclear component that is opaque below about 20 GHz.  The radio galaxy 3C 123 is transparent at all plotted frequencies, and its spectrum steepens above a few GHz.  The quasar 3C 48 is synchrotron self-absorbed only below 100 MHz, while the quasar 3C 454.3 contains structures that become opaque at widely differing frequences.  [Kellermann, K. I., & Owen, F. N. 1988, in Galactic and Extragalactic Radio Astronomy, eds. G. L. Verschuur & K. I. Kellermann (Springer Verlag)]