Information about the kinematics plotter

kinematics calculator and plotter | github repo for kinematics software

Show historical background and plot example

Two-body reaction kinematics

Initial state

Consider a particle with mass m1 and kinetic energy T1 incident on a stationary particle with mass m2 . We can write a reaction of these particles to create two new particles as m1+ m2 m3+ m4 . The energy of the projectile is E1 = m1 + T1 . The momentum of the projectile is p1 = E12 - m12 = T1 (T1 +2 m1) .

A Lorentz invariant for the initial projectile m1 and target m2 system is Etot2- ptot2 , with Etot= ( m1 + T1 + m2 ) and ptot= p1 . Algebraic manipulation gives the simplified Lorentz invariant s = ( m1 + m2 )2 + 2 m2 T1 For identical particles with opposing momenta, the second term would be 2mT 4(m+T) T . So a pair of accelerators gets an additional factor of quadratic gain in energy available for reactions, which is why the largest accelerators like CERN have collided opposing beams of particles.

Center of momentum (CM) reference frame

In the center of momentum (CM) frame, the momenta of m1 and m2 are equal and opposite. The invariant quantity s does not change. In the CM frame it can be expressed as s = ( m12 + pcm2 + m22 + pcm2 ) 2 Since s is an invariant constant, this equation defines the CM momentum. pcm = ( s - m12 - m22 ) 2 - 4 m12 m22 4 s This is the momentum of both particles in the CM frame. To determine the rapidity χ in the CM frame, we consider that a boost to the target with the CM rapidity gives the target the CM momentum: ( m2 , 0 ) ( E2 , p2 ) = ( m2 cosh(χ) , m2 sinh(χ) ) A simple fact about hyperbolic trigonometric functions, coshχ + sinhχ = eχ , allows us to calculate the CM rapidity in closed form: χ = ln ( pcm + m22 + pcm2 m2 )

So far we have found the invariant s and CM rapidity χ from the given initial values of the projectile mass, target mass and projectile kinetic energy. To compute the properties of the reaction products, we must also know the ejectile and recoil masses, m3 and m4. For clarity, we will refer to the CM momenta of the products as q , since it will generally be different from pcm. The relationship between the product masses, the Lorentz invariant and the CM momentum is q = ( s - m32 - m42 ) 2 - 4 m32 m42 4 s

In principle, this is all that is needed to calculate the energies and momenta of the products as a function of either lab angle or CM angle. These relationships are fully determined by conservation of energy and momentum. The quantity of interest in the relevant nuclear and particle physics experiments is the probability for a certain set of reaction products to be produced as a function of the initial energy ( Ecm or s, determined by T1 ) and angle (often θ3cm , inferred from the measured angle in the laboratory of θ3).

Outgoing product states

In the CM frame, the reaction products each have the CM momentum. Their CM energies are E3cm = q2 + m32 and E4cm = q2 + m42 . The corresponding lab energies and momenta are found by applying a boost with the CM rapidity χ to these values as a function of CM angle. The ejectile energy, parallel momentum and perpendicular momentum are given by E3 = E3cmcoshχ + q sinhχ cosθ3cm p3 cosθ3 = q coshχ cosθ3cm + E3cmsinhχ p3 sinθ3 = qsinθ3cm The same values for the recoil are given by the same boost with the sign changed on q: E4 = E4cmcoshχ - q sinhχ cosθ3cm p4 cosθ4 = -q coshχ cosθ3cm + E4cmsinhχ p3 sinθ4 = qsinθ3cm

These equations fully determine the outgoing energies as a function of angle. It is also possible to calculate the lab energy as a function of the lab angle for the ejectile and recoil. To do this, isolate the terms with the CM angle and CM momentum, then square and add them to eliminate the CM angle by using the trigonometric identity sin2θ + cos2θ =1 . Further rearrangement and application of the hyperbolic trigonometric identity cosh2χ - sinh2χ =1 yield quadratic equations for both the ejectile and recoil momomenta as a function of their lab angles: p32 ( 1+ sinh2χ sin2θ3 ) - 2 p3 E3cm coshχ cosθ3 + m32 sinh2χ - q2 =0 p42 ( 1+ sinh2χ sin2θ4 ) - 2 p4 E4cm coshχ cosθ4 + m42 sinh2χ - q2 =0

The solution to the quadratic equation gives p3 = E3cm sinhχ cosθ3 ± coshχ q2 - m32 sinh2χ sin2θ3 1 + sinh2χ sin2θ3 The same solution applies with the 3 subscripts changed to 4. 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. If m3sinhχ > q , then this implies a maximum lab angle for the ejectile of θ3 = arcsin ( q m3sinhχ ) The same condition and solution apply to the recoil as well, using its mass instead.

If there is a maximum scattering angle, then both signs apply in the equation for momentum as a function of angle. If there is no maximum angle, then the positive sign is correct. To see this, consider the case where one of the final state partcles is a photon. The numerator in the equation for momentum as a function of angle in the massless limit becomes qsinhχcosθ ± qcoshχ . The magnitude of the momentum must not be negative, so the positive sign is correct.

Limit for elastic scattering of target

Consider the limiting form of the equations for the recoil at very small values of θ3cm. For elastic scattering, E4cm=m2coshχ. For infinitesimal values of the CM angle, we write its cosine as 1-ε2/2 and its sine as ε. To leading order, we find E4 = m2 ( 1+sinh2χ ε2/2 ) p4 cosθ4 = m2 sinhχ coshχ ε2/2 p4 sinθ4 = m2 sinhχ ε To leading order, the kinetic energy takes the form of the non-relativistic limit and it is completely explained by the perpendicular component of the momentum. The component of momentum parallel to the incident momentum is only proportional to the infinitesimal quantity squared and we can neglect it. The target's momentum is perpendicular to the incident projectile and its limiting angle is 90 degrees, in the case of elastic scattering.

Limiting angle for elastic scattering of a heavy particle from a light particle

Consider one more elastic scattering limit. Assume that the projectile is a heavy particle and the target is a light particle. Just to be specific, let the heavy particle be a proton and the light particle be an electron. In the case of elastic scattering, the initial particles are the same as the final particles, so for our example we have q = me sinhχ . The formula for the maximum scattering angle becomes θ3 = arcsin ( m4 m3 ) θp = me mp The portion written in terms of the electron and proton masses has used the small angle approximation, but this is quite valid for the mass ratio of approximately 1/2000. This result is completely general for any case of a heavy particle scattering off a light particle. Just to give one more concrete example, the maximum scattering angle of a deuteron incident on a proton is the inverse sine of 1/2, or 30 degrees.

Conversion between laboratory and CM frames

For planning and analyzing experiments, it is often useful to convert between CM 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 CM angle, one possible expression is tanθ = qsinθcm qcosθcm coshχ + Ecmsinhχ Take the implicit derivative of both sides with respect to angle to find dθ dθcm = q p ( cosθcm cosθ + sinθcm sinθ coshχ ) Finally, multiply by the ratio of sines to find dΩ dΩcm = q2 p2 ( cosθcm cosθ + sinθcm sinθ coshχ ) The above equation was derived for the ejectile. The only difference for the recoil is that the first term in parentheses has a minus sign.