Ortho-positronium annihilation: steps toward computing the first order radiative corrections

The 0.2% experimental accuracy of the 1968 Beers and Hughes measurement of the annihilation lifetime of ortho-positronium motivates the attempt to compute the first order quantum electrodynamic corrections to this lifetime. The theoretical problems arising in this computation are here studied in detail up to the point of preparing the necessary computer programs and using them to carry out some of the less demanding steps -- but the computation has not yet been completed. Analytic evaluation of the contributing Feynman diagrams is superior to numerical evaluation, and for this process can be carried out with the aid of the Reduce algebra manipulation computer program.

The relation of the positronium decay rate to the electronpositron annihilation-in-flight amplitude is derived in detail, and it is shown that at threshold annihilation-in-flight, Coulomb divergences appear while infrared divergences vanish. The threshold Coulomb divergences in the amplitude cancel against like divergences in the modulating continuum wave function.

Using the lowest order diagrams of electron-positron annihilation into three photons as a test case, various pitfalls of computer algebraic manipulation are discussed along with ways of avoiding them. The computer manipulation of artificial polynomial expressions is preferable to the direct treatment of rational expressions, even though redundant variables may have to be introduced.

Special properties of the contributing Feynman diagrams are discussed, including the need to restore gauge invariance to the sum of the virtual photon-photon scattering box diagrams by means of a finite subtraction.

A systematic approach to the Feynman-Brown method of Decomposition of single loop diagram integrals with spin-related tensor numerators is developed in detail. This approach allows the Feynman-Brown method to be straightforwardly programmed in the Reduce algebra manipulation language.

The fundamental integrals needed in the wake of the application of the Feynman-Brown decomposition are exhibited and the methods which were used to evaluate them -- primarily dis persion techniques are briefly discussed.

Finally, it is pointed out that while the techniques discussed have permitted the computation of a fair number of the simpler integrals and diagrams contributing to the first order correction of the ortho-positronium annihilation rate, further progress with the more complicated diagrams and with the evaluation of traces is heavily contingent on obtaining access to adequate computer time and core capacity.

Plasma loop and strapping field dynamics: reproducing solar eruptions in the laboratory

Coronal mass ejections (CMEs) are dramatic eruptions of large, plasma structures from the Sun. These eruptions are important because they can harm astronauts, damage electrical infrastructure, and cause auroras. A mysterious feature of these eruptions is that plasma-filled solar flux tubes first evolve slowly, but then suddenly erupt. One model, torus instability, predicts an explosive-like transition from slow expansion to fast acceleration, if the spatial decay of the ambient magnetic field exceeds a threshold.

We create arched, plasma filled, magnetic flux ropes similar to CMEs. Small, independently-powered auxiliary coils placed inside the vacuum chamber produce magnetic fields above the decay threshold that are strong enough to act on the plasma. When the strapping field is not too strong and not too weak, expansion force build up while the flux rope is in the strapping field region. When the flux rope moves to a critical height, the plasma accelerates quickly, corresponding to the observed slow-rise to fast-acceleration of most solar eruptions. This behavior is in agreement with the predictions of torus instability.

Historically, eruptions have been separated into gradual CMEs and impulsive CMEs, depending on the acceleration profile. Recent numerical studies question this separation. One study varies the strapping field profile to produce gradual eruptions and impulsive eruptions, while another study varies the temporal profile of the voltage applied to the flux tube footpoints to produce the two eruption types. Our experiment reproduced these different eruptions by changing the strapping field magnitude, and the temporal profile of the current trace. This suggests that the same physics underlies both types of CME and that the separation between impulsive and gradual classes of eruption is artificial.