Remembering John Brown

From RHESSI Wiki

Revision as of 20:28, 1 January 2020 by Amackinnon (Talk | contribs)
Jump to: navigation, search

John C Brown (1947-2019)

John Brown passed away unexpectedly on 16 November 2019, a great sadness to all who knew him and a huge loss particularly to the Glasgow Astronomy and Astrophysics group. His many contributions to solar physics were highly influential and we were prompted to think back over some of the big themes of his research career. In this nugget we look at one of those big themes, inverse problems in astrophysics. This particular topic spans the whole of John's research career and combines a technical side with a more philosophical approach to the character of data and the limits on what one can deduce from them.


Brown(1971)

John's first paper was published in 1971. As of 22 November 2019, ADS says it has been cited 889 times. It introduced several ideas that became key to the study of flares and their X-radiation. A previous nugget dealt with the thick target model. Another new idea was the possibility of "inverting" observed X-ray spectra to deduce the energy distribution of the emitting electrons. Electrons are accelerated to 10s of keV energy in the flare. The hard X-rays they emit via bremsstrahlung yield information on their number and energies.

The bremsstrahlung cross-section dσ/dε gives the probability that an electron of energy E will emit a photon of energy ε. Integrating over all electrons and over the whole of the source volume, the flux at Earth is I(ε):

Equation 1

Failed to parse (PNG conversion failed; check for correct installation of latex, dvips, gs, and convert): I(\epsilon) = \frac{1}{4\pi d^2} \int_\epsilon^\infty \frac{mathrm{d}\sigma}{\mathrm{d}\epsilon} v F(E) \mathrm{d}E

Here d= 1 AU, v is the speed of an electron of energy E and Failed to parse (syntax error): F(E) = \int n({\bf r}) \int f(E,\Omega,{\bf r}) \mathrm{d}^2\Omega \mathrm{d}^3{\bf r} <\math> i.e., the angle-averaged, density-weighted electron distribution. With limiting assumptions about the character of the source (thick or thin target), ''F(E)'' can be simply reduced to the flux of electrons being injected to the X-ray emission region - information that should be intimately related to the behaviour and physical character of the energy release process. The key task, then, is to deduce ''F(E)'' from observed spectra. In the simplest approach one assumes a parametric form, e.g. a power-law in energy ''E<sup>-&delta;'', and adjusts the parameters (e.g. &delta;) to obtain the best fit. John explored the opposite process, direct deduction of ''F(E)'' from the observed properties of ''I(&epsilon;)''. Instead of susbstituting assumed forms of ''FE()'' to see what comes out, He spotted that one could go a long way in this direction by using the Born approximation bremsstrahlung cross-section due to Bethe and Heitler, in its non-relativistic limit: <blah> Here <math>J(\epsilon) = K \int_\epsilon^\infty blah blah \mathrm{d}E

Here J is the energy flux detected at Earth (erg keV-1 cm-2 s-1) and N is blah blah, after an appropriate average over the source volume. The probability of an energetic particle emitting a photon in a close encounter with another charged particle is given by the cross-section. John used the Born approximation cross-section differential in photon energy due to Bethe and Heitler, in its non-relativistic limit. This is accurate enough in the relevant energy regime while also being simple enough to allow a largely analytical discussion. Then there is a solution for F(E) in terms of I(ε) and its derivatives. 

Failed to parse (PNG conversion failed; check for correct installation of latex, dvips, gs, and convert): \frac{\mathrm{d}\sigma}{\mathrm{d}\epsilon} = \frac{\sigma_0}{\epsilon E}\ln\{\frac{1+\sqrt{1-\epsilon/E}}{1-\sqrt{1-\epsilon/E}}\}

This was both reasonably accurate for the relevant energy range, and simple enough algebraically to allow a largely analytical discussion.

With some algebra, this expression can be rewritten Failed to parse (PNG conversion failed; check for correct installation of latex, dvips, gs, and convert): J(\epsilon) \sim \int_\epsilon^\infty F(E) \frac{N(E)}{\sqrt{E-\epsilon}} \mathrm{d}E

This is Abel's integral equation, which has a solution for N(E) in terms of J(ε) and its derivatives:


What a discovery! Before this people had tried to fit observed hard X-ray spectra by assuming parametric forms for N(E), e.g. a power-law in electron energy E and looking for the value of δ that gave the best fit to the observed spectrum. Now we realise that the electron distribution in the source can be deduced directly, in a non-parametric way - if the observations giving the spectrum J(ε) are sufficiently free from noise, and sufficiently well resolved. The scintillators in use in the 1970s had a typical energy resolution δε/ε of about 10%, inadequate for implementation of this approach but the possibility might be realisable with other sorts of X-ray spectrometer.

If we use more precise expressions for the cross-section the analytical inversion of Abel's integral equation is no longer applicable but a numerical (matrix) process can still be constructed, deducing a vector of values (N(E1),N(E2),N(E3)...N(EM)) from a histogram of values (J1, J2,...,JL).

Ill-conditioned problems

A worse difficulty lay just below the surface.

Personal tools
Namespaces
Variants
Actions
Navigation
Toolbox