Breit–Wheeler Pair Production Cross Section Calculator

Breit–Wheeler Pair Creation in Photon–Photon Collisions

The direct conversion of light into matter is one of the most striking predictions of quantum electrodynamics. In the Breit–Wheeler process, two photons collide and produce an electron–positron pair, written as γγ → e⁺e⁻. Although the idea is conceptually simple, the conditions required are demanding: the photons must carry enough combined energy in the center-of-mass frame to create the rest mass of both particles. That threshold requirement is why the process matters so much in high-energy astrophysics, laser-plasma physics, and studies of the early universe. This calculator gives a quick numerical estimate of the unpolarized Breit–Wheeler cross section for a head-on collision using two photon energies entered in MeV.

Cross section is a probability-like measure of how likely an interaction is to occur. A larger cross section means the process is more effective when beams or radiation fields overlap. For photon–photon pair production, the cross section is zero below threshold, rises rapidly once pair creation becomes kinematically allowed, reaches a maximum somewhat above threshold, and then gradually decreases at very high energies. That shape is important in practice. It helps determine how far gamma rays can travel through ambient radiation fields, how efficiently intense laser setups may generate pairs, and where attenuation becomes significant in astrophysical environments.

In natural units where ℏ = c = 1, the center-of-mass energy squared for two photons of energies E₁ and E₂ colliding head-on is s = 4E₁E₂. Pair production is possible only if s exceeds 4mₑ², where mₑ is the electron mass. Equivalently, the product of the two photon energies must satisfy E₁E₂ ≥ mₑ². Using mₑ = 0.511 MeV, the threshold product is about 0.261 MeV². If the product is smaller than that value, the calculator reports that the reaction is below threshold and sets the cross section to zero.

Introduction

This page is designed for readers who want both a working calculator and a clear explanation of what the result means. You do not need to derive the full quantum electrodynamics result to use the tool, but it helps to know the physical picture. Two photons have no rest mass, yet if they collide with enough energy they can create massive particles. In this case the final state is an electron and a positron. Because the initial particles are photons, there is no material target in the usual sense; the interaction is controlled entirely by the collision geometry and the available invariant energy.

The calculator assumes the standard head-on configuration, which is the simplest and most common form used for quick estimates. In that limit, the energy product E₁E₂ directly determines whether the threshold is crossed and how large the cross section becomes. This makes the tool useful for classroom exploration, back-of-the-envelope astrophysical estimates, and preliminary planning for more detailed simulations. It is especially helpful when you want to see how changing one photon energy can compensate for changing the other, since the threshold depends on their product rather than on either energy alone.

Because the result is reported in barns, it is easy to compare with familiar cross-section scales from nuclear and particle physics. One barn equals 10⁻²⁸ m². The Breit–Wheeler cross section is typically of order one barn near its peak, which is large enough to matter in dense radiation fields but still small enough that ordinary laboratory observation is challenging unless photon densities are very high. That balance between conceptual simplicity and experimental difficulty is one reason the process remains so interesting.

How to Use

Using the calculator is straightforward. Enter the energy of the first photon in the field labeled E₁ and the energy of the second photon in the field labeled E₂. Both values should be given in MeV. Then press the compute button. The script multiplies the two energies, checks whether the threshold condition is satisfied, and if so evaluates the Breit–Wheeler formula. The result area then displays the dimensionless speed parameter β and the cross section σ in barns.

When interpreting the inputs, remember that the calculator is built for head-on photon collisions. If your physical situation involves a different collision angle, the true invariant energy would be different from the simple head-on expression used here. For many educational and quick-estimate purposes, however, the head-on case is exactly the right place to start because it gives the maximum center-of-mass energy for fixed photon energies.

A few practical tips help avoid confusion. First, use positive energies only. Second, keep the units consistent: the form expects MeV, not eV, keV, or GeV. Third, if you are exploring threshold behavior, try values whose product is just below and just above 0.261 MeV². You will see the result switch from zero to a finite cross section, which is a good way to build intuition for the kinematics. Finally, note that two very different energy pairs can produce the same answer if their product E₁E₂ is the same.

Formula

The exact Breit–Wheeler cross section for unpolarized photons in the head-on case is written in terms of the lepton speed parameter β. This parameter is defined by the available invariant energy and ranges from 0 at threshold to values approaching 1 far above threshold. Physically, β represents the speed of the produced electron and positron in the center-of-mass frame, measured in units of the speed of light.

The calculator uses the following expression:

The script computes β from the threshold relation

using the implementation form β = √(1 - mₑ²/(E₁E₂)) for head-on collisions. The classical electron radius is taken as rₑ ≈ 2.8179 × 10⁻¹⁵ m. After evaluating the formula in square meters, the code converts the result to barns by dividing by 10⁻²⁸.

Near threshold, β is small, so the cross section starts at zero because there is almost no available phase space for the outgoing particles. As the energy increases, β grows and the cross section rises. At still higher energies, the logarithmic growth inside the bracket is outweighed by the overall suppression, so the cross section eventually decreases. This non-monotonic behavior is a hallmark of the Breit–Wheeler process and explains why the largest interaction probability occurs not at threshold and not at arbitrarily high energy, but in an intermediate range above threshold.

Worked Example

Suppose the first photon has energy E₁ = 2 MeV and the second has energy E₂ = 1 MeV. Their product is 2 MeV², which is well above the threshold value of about 0.261 MeV². That means pair production is allowed. The calculator then evaluates β = √(1 - 0.511²/2), which gives a value close to 0.93. Substituting that into the Breit–Wheeler expression produces a cross section of order a barn. The exact displayed value depends on the numerical evaluation in the script, but the important point is that this energy pair lies in the regime where the process is quite efficient compared with values just above threshold.

Now compare that with E₁ = 1 MeV and E₂ = 0.2 MeV. The product is only 0.2 MeV², which is below threshold. In that case the calculator returns “Below threshold: σ = 0.” This is not a numerical artifact or a rounding issue; it is a direct consequence of energy conservation. There simply is not enough invariant energy to create the rest mass of the electron–positron pair.

You can also test the symmetry of the result. If you enter E₁ = 10 MeV and E₂ = 0.2 MeV, the product is again 2 MeV², so the cross section is the same as for 2 MeV and 1 MeV in the head-on approximation. This is a useful reminder that the interaction depends on the invariant combination of the two photon energies, not on which photon is labeled first or second.

Interpreting the Result

The output contains two pieces of information. The first is β, the speed parameter of the produced leptons in the center-of-mass frame. Values near zero mean the system is just above threshold, so the electron and positron emerge slowly in that frame. Values near one indicate a collision far above threshold, where the outgoing particles are highly relativistic. The second quantity is the cross section σ in barns. This tells you the effective interaction strength for the chosen energy pair under the assumptions of the model.

A result near zero can mean one of two things. Either the energies are below threshold, in which case pair production cannot occur at all, or the energies are so close to threshold that the process is allowed but still strongly suppressed. A moderate positive value, often around one barn, indicates that the chosen energies are in the most favorable range for pair creation. If you continue increasing the energies while keeping the collision head-on, you will eventually see the cross section decrease again even though β keeps increasing. That decline is expected from the analytic form of the Breit–Wheeler formula.

In astrophysical applications, the cross section is usually not the final quantity of interest by itself. Instead, it enters an integral over a distribution of target photons. Even so, a single-energy calculation is extremely useful because it shows where the kernel of that integral is large or small. In laboratory contexts, the same result can help estimate whether a proposed photon source is likely to produce an observable number of pairs when combined with another beam or radiation field.

Representative Values

The table below lists cross sections for several representative energy pairs. It illustrates the general trend: zero below threshold, a rapid rise above threshold, a broad maximum, and then a gradual decline at high energy.

These values are meant as orientation points rather than a substitute for calculation. If you are studying a specific problem, use the form below with your own energies. The examples are still useful because they show the scale of the answer and confirm that the cross section does not simply increase forever with energy.

Limitations and Assumptions

This calculator intentionally uses a simplified but standard setup. The most important assumption is that the photons collide head-on. In a more general geometry, the invariant energy depends on the collision angle, and the threshold condition changes accordingly. If the photons are not counter-propagating, the same pair of energies may produce a smaller center-of-mass energy than the calculator assumes. As a result, the tool is best understood as a head-on estimate rather than a universal photon–photon collision solver.

The formula implemented here is the tree-level Breit–Wheeler result for unpolarized photons. It does not include polarization-dependent effects, higher-order quantum corrections, beam spread, finite pulse duration, angular distributions of the outgoing leptons, or environmental effects such as external fields and plasma backgrounds. Those refinements matter in precision studies and in strongly non-linear laser regimes, but they are beyond the scope of a compact educational calculator.

Another limitation is numerical rather than physical. The script uses ordinary double-precision JavaScript arithmetic, which is more than adequate for typical educational and exploratory use. However, if you are performing high-precision research calculations, integrating over broad photon spectra, or comparing with detailed Monte Carlo simulations, you should treat this page as a quick estimator and not as a replacement for a dedicated computational framework.

Even with those limitations, the calculator remains valuable because it captures the essential threshold physics and the correct qualitative energy dependence of the Breit–Wheeler cross section. It is a practical way to build intuition about when light can turn into matter and how strongly that conversion proceeds once the kinematic barrier is crossed.