Figure 1: Illustration of the lab frame and center-of-mass frame, drawn to satisfy classical addition of the velocities if the masses are set to 1. Relativistic velocity addition is not so simple.

Consider a particle, with mass
*m _{2}*
and kinetic energy

It is common to call
*m _{2}* and

So the Lorentz invariant quantity for the initial two-particle system is

.

Algebraic manipulation gives the simplified result,

.
Note that if identical particles were to collide with equal and opposite momenta in the lab frame, the portion of the invariant mass available for reactions (the second term in the previous equation) would go like
4
(*mT _{k}*
+

However, I do nuclear physics experiments, where we play different games than the particle physicists. One is to make precise measurments at low energies that are relevant to astrophysics: for reactions in stars, novae, supernovae, and big bang cosmology. Another is to perform experiments with beams of radioactive nuclei, to study exotic nuclei to further our understanding of nuclear structure. So enough of this particle physics talk about big accelerators!

In the center-of-mass frame, the momenta of
*m _{1}* and

Since *s* is just a constant, this equation determines
*p _{cm}*.

*p _{cm}* is the momentum of both particles in the center-of-mass frame. To determine the rapidity of the center-of-mass frame, we consider that a boost to the target with the cm rapidity gives the taget the momentum

A simple fact about the hyperbolic trig functions,

,

allows us to calculate the cm rapidity in a nice closed form:

.

Given the invariant mass and the rapidity in the cm frame, it is simple to calculate the energies of the reaction products as a function of angle. The cm momentum of the reaction products follows from the previous equations, but with
*m _{1}* and

Figure 2: Triangles formed by the reaction momenta and angles in the cm and lab frames. The picture is strictly classical, using Galilean addition of the velocities, and
*v _{3}* and

Given
*E _{k}* and all the masses,

The only undetermined parameter is the cm angle. The game we usually play when doing nuclear physics experiments is to determine the probability for a given reaction to happen, either as a function of the cm energy, or as a function of the cm angle, or both.

All measurements are made with detectors in a laboratory environment, so it is important to convert readily between cm and lab quantities.

Inspecting the simplified figure on the right, all the signs for the boost to the ejectile are positive. The ejectile energy and momenta in the lab are

,

, and

.

There is a minus sign to consider for the boost to the recoil. The recoil energy and momenta in the lab are

,

, and

.

That's really all there is to it. Now, if you were deranged and masochistic, you might want to calculate one of the lab energies as a function of the lab angle. This is possible, but quite ugly. To do it, begin by rearranging the momenta so that the cm angle terms are isolated and square them both, so that addition of the squared sine and cosine terms eliminates the cm angle altogether. The remaining equations are quadratic. The form of the quadratic equation is identical for
*p _{3}* and

The quadratic equation gives each momentum as a function of its angle:

These solutions merit some discussion. The quantity under the radical in the second term of the numerator must be greater than or equal to zero, since the momentum must be real. Depending on the kinematics, this could imply a limit on the lab angle. The condition is,

When this condition applies, the positive sign gives the high-energy reaction products, the positive and negative sign solutions become equal at the maximum lab angle (the term under the radical becomes zero), and the negative sign solution gives the low-energy reaction products. This condition usually applies for one of the reaction products in nuclear physics experiments, because beam energies are generally small compared to nuclear masses on the order of many GeV. The lighter reaction products, on the other hand, often have no maximum lab angle.

If the condition does not apply, so that there is no maximum lab angle for the given product, then the positive sign is correct. To see this, imagine that one of the outgoing particles is a photon. Taking the limit as *m/p' _{cm}* goes to zero, the numerator becomes

.

Since the magnitude of the momentum can never be negative, in this case the positive sign is correct.

For planning and analyzing experiments, it is often useful to convert between center-of-mass and laboratory cross sections. To do this, one must know how the solid angles vary with respect to each other. Using the relativistic transformations to write the laboratory angle as a function of the center-of-mass angle, one possible expression is

.

To confuse persistent readers, the notation has been changed. *p _{cm}* no longer has the apostrophe but it still refers to the center-of-mass momentum of the

Taking the implicit derivative of both sides of this equation, the derivative of the laboratory angle with respect to the center-of-mass angle is

.

Multiply this equation by the ratio of the sines to find

This equation is for the ejectile. The only difference between the ejectile and the recoil is the minus sign, so the result is the same for the recoil with a minus sign in front of the cos(*θ _{cm}*) term.

To make plots of lab energy vs lab angle using these equations, try the Relativistic Kinematics Calculator for Two-Body Nuclear Reactions. Please send any questions, comments, or suggestions to Sky.