Bethe-Heitler Pair Production

• Bethe-Heitler pair production is a QED process where high-energy photons convert into electron-positron pairs in the Coulomb field of nuclei or protons, defined by distinct threshold energies and screening effects.
• It is modeled using Born-level approximations, Monte Carlo algorithms, and particle-in-cell simulations to capture cross section scaling, interference, and various radiative corrections.
• The process finds applications in astrophysical jet modeling, gamma-ray burst studies, and all-optical positron-source design, highlighting its relevance across high-energy physics and plasma environments.

to=arxiv_search 彩神争霸破解_code: {"query":"Bethe-Heitler pair production", "max_results": 10} to=arxiv_search 东臣_code: {"all_fields":"Bethe-Heitler pair production","start":0} Bethe-Heitler pair production usually denotes the conversion of a high-energy photon into an electron-positron pair in the Coulomb field of a nucleus, $\gamma + (Z) \rightarrow e^+ + e^- + (Z)$. In several astrophysical papers, the same name also denotes the photopair channel $p+\gamma \rightarrow p+e^+ + e^-$. Across these usages, the process functions as an important QED mechanism, a Monte Carlo module in particle-in-cell simulations, a source of secondary pairs in relativistic plasmas, and an irreducible Standard-Model background in fixed-target searches for light mediators (Martinez et al., 2019, Petropoulou et al., 2014, Pustyntsev et al., 18 May 2026).

1. Canonical process and kinematic domain

In the nuclear-field formulation, the process is treated in the rest frame of a heavy nucleus of charge $Ze$. The threshold condition is that the photon energy in that frame satisfy $\omega \ge 2 m_e c^2$. Energy conservation may be written schematically as $\omega + M \simeq E_+ + E_- + M'$, while momentum conservation requires the nuclear recoil momentum even when the recoil energy is negligible. A standard approximation is to treat the nucleus as an external field, neglect nuclear recoil, and work at leading order in the Born approximation (Karavola et al., 2024).

A representative Born-level amplitude in the heavy-nucleus rest frame has the form

$$M^\mu \approx e Z\, \bar u(p_-)\gamma^\mu v(p_+) / q^2,$$

with $q^\mu \approx p_+ + p_- - k$ the momentum transfer to the nucleus. In this formulation the cross section scales with $Z^2$, and screening enters through the atomic form factor $F(q)$ (Martinez et al., 2019).

In hadronic astrophysics, the same mechanism is recast as proton-photon photopair production,

$p+\gamma \rightarrow p+e^+ + e^-.$

The threshold in the proton rest frame is $p+\gamma \rightarrow p+e^+ + e^-$0, and one frequently uses the head-on condition

$p+\gamma \rightarrow p+e^+ + e^-$1

for a proton of Lorentz factor $p+\gamma \rightarrow p+e^+ + e^-$2 interacting with a target photon of energy $p+\gamma \rightarrow p+e^+ + e^-$3. A more exact threshold quoted for the proton rest frame is

$p+\gamma \rightarrow p+e^+ + e^-$4

with the lab-frame condition $p+\gamma \rightarrow p+e^+ + e^-$5 (Petropoulou et al., 2014, Bégué et al., 2021, Romero et al., 3 Sep 2025).

A common misconception is that “Bethe-Heitler pair production” always refers to a single kinematic setup. The literature here shows two recurrent realizations: photon conversion in an external Coulomb field and photopair production off a proton. The underlying threshold structure and pair-production kernel are related, but the target, observable, and approximation scheme differ.

2. Born-level cross sections and screening structure

For an unpolarized photon of energy $p+\gamma \rightarrow p+e^+ + e^-$6 converting in the Coulomb field of a nucleus of charge $p+\gamma \rightarrow p+e^+ + e^-$7, one form of the Born-level differential cross section is

$p+\gamma \rightarrow p+e^+ + e^-$8

where $p+\gamma \rightarrow p+e^+ + e^-$9. In the leading-logarithm approximation with screening cutoff $Ze$0,

$Ze$1

These expressions make explicit the $Ze$2 scaling and the sensitivity to the infrared cutoff provided by screening (Zhu, 2019).

For ultra-relativistic photons and leptons in partially ionized plasmas, a screened form used in particle-in-cell work is

$Ze$3

with $Ze$4, $Ze$5, and $Ze$6. The screening corrections are encoded in

$Ze$7

$Ze$8

Here $Ze$9 is the atomic form factor (Martinez et al., 2019).

The same PIC-oriented treatment models the atomic Coulomb potential in an arbitrarily ionized plasma as a Thomas-Fermi-plus-Debye screened Yukawa potential,

$\omega \ge 2 m_e c^2$0

with

$\omega \ge 2 m_e c^2$1

The corresponding form factor obeys

$\omega \ge 2 m_e c^2$2

This construction reduces to the usual TF-only result at $\omega \ge 2 m_e c^2$3 and to Debye-only screening at full ionization (Martinez et al., 2019).

These formulations show that the “Bethe-Heitler cross section” is not a single immutable expression. Screening, ionization state, recoil assumptions, and asymptotic limits alter both the analytic form and the computational strategy.

3. Angular momentum transfer, polarization, and strong-field generalizations

A twisted-state treatment of $\omega \ge 2 m_e c^2$4 makes the angular-momentum content explicit. For a circularly polarized plane-wave photon with helicity $\omega \ge 2 m_e c^2$5, and twisted final leptons labeled by orbital index $\omega \ge 2 m_e c^2$6 and spin projection $\omega \ge 2 m_e c^2$7, the amplitude obeys the selection rule

$\omega \ge 2 m_e c^2$8

With $\omega \ge 2 m_e c^2$9 and $\omega + M \simeq E_+ + E_- + M'$0, the averages satisfy

$\omega + M \simeq E_+ + E_- + M'$1

The reported result is that the spin angular momentum of the incident photon is converted not only into spin angular momentum of the produced pair but also into orbital angular momentum, and that the average OAM gained by the leptons surpasses the average SAM while their orientations coincide. For $\omega + M \simeq E_+ + E_- + M'$2 on a Cu nucleus, the dominant channel has $\omega + M \simeq E_+ + E_- + M'$3 and $\omega + M \simeq E_+ + E_- + M'$4; numerically $\omega + M \simeq E_+ + E_- + M'$5 and $\omega + M \simeq E_+ + E_- + M'$6 for each lepton. Raising $\omega + M \simeq E_+ + E_- + M'$7 to $\omega + M \simeq E_+ + E_- + M'$8–$\omega + M \simeq E_+ + E_- + M'$9 or selecting larger opening angles increases the width of the $$M^\mu \approx e Z\, \bar u(p_-)\gamma^\mu v(p_+) / q^2,$$0-spectrum and further enhances $$M^\mu \approx e Z\, \bar u(p_-)\gamma^\mu v(p_+) / q^2,$$1 (Lei et al., 2023).

In intense laser fields, Bethe-Heitler pair creation becomes a nonlinear process of the form

$$M^\mu \approx e Z\, \bar u(p_-)\gamma^\mu v(p_+) / q^2,$$2

A Volkov-state $$M^\mu \approx e Z\, \bar u(p_-)\gamma^\mu v(p_+) / q^2,$$3-matrix in a bichromatic field leads to rates summed over photon numbers $$M^\mu \approx e Z\, \bar u(p_-)\gamma^\mu v(p_+) / q^2,$$4, with energy-momentum conservation enforced by

$$M^\mu \approx e Z\, \bar u(p_-)\gamma^\mu v(p_+) / q^2,$$5

When the two laser frequencies are commensurable, distinct absorption pathways can reach the same final energy and interfere. The interference terms are phase sensitive, modulate the total rate $$M^\mu \approx e Z\, \bar u(p_-)\gamma^\mu v(p_+) / q^2,$$6, and visibly reshape the angular distributions. In the weak-intensity regime, the partial rates scale as powers of the intensity parameters $$M^\mu \approx e Z\, \bar u(p_-)\gamma^\mu v(p_+) / q^2,$$7, and the relative phase can switch interference from destructive to constructive (Augustin et al., 2013, Augustin et al., 2013).

A separate nonperturbative analysis considers a relativistic proton beam colliding with a bichromatic laser composed of a strong low-frequency mode and a weak high-frequency assisting mode. In the pure strong-field limit one recovers

$$M^\mu \approx e Z\, \bar u(p_-)\gamma^\mu v(p_+) / q^2,$$8

whereas the assisting mode lowers the effective tunneling barrier and yields

$$M^\mu \approx e Z\, \bar u(p_-)\gamma^\mu v(p_+) / q^2,$$9

Numerically, the singly assisted channel dominates up to $q^\mu \approx p_+ + p_- - k$0 and provides an enhancement by up to two orders of magnitude around $q^\mu \approx p_+ + p_- - k$1–$q^\mu \approx p_+ + p_- - k$2. The assisted spectra are broader in angle and energy, and the lab-frame emission is compressed into a few-degree cone by the Lorentz boost (Augustin et al., 2014).

4. Radiative corrections and alternative high-energy formulas

For the reaction $q^\mu \approx p_+ + p_- - k$3, the Born-level Bethe-Heitler amplitude consists of the two standard tree diagrams in which the lepton pair is emitted from the photon line and the proton couples through its electromagnetic vertex. In the soft-photon approximation, the first-order QED correction factorizes on the Born cross section,

$q^\mu \approx p_+ + p_- - k$4

The virtual and real soft terms are infrared divergent separately but cancel in the sum. For the kinematics at MAMI, the reported corrections are of the percent level for muons and of order $q^\mu \approx p_+ + p_- - k$5 for electrons. In the proposed lepton-universality observable, these corrections shift the cross-section ratio by $q^\mu \approx p_+ + p_- - k$6–$q^\mu \approx p_+ + p_- - k$7, which is larger than the target experimental precision $q^\mu \approx p_+ + p_- - k$8 (Heller et al., 2018).

A distinct proposal is the “improved” static-limit formula for pair creation in a very thin, fully ionized gas treated as infinitely massive and completely unscreened out to large impact parameters. In that limit the differential cross section is written as

$q^\mu \approx p_+ + p_- - k$9

to be compared with the standard screened Bethe-Heitler form proportional to $Z^2$0. The quoted ratio is

$Z^2$1

For $Z^2$2 and $Z^2$3, the reported enhancement is several thousand compared to the Bethe-Heitler prediction (Zhu, 2019).

This suggests that the status of recoil and screening assumptions is not merely technical. Within the literature summarized here, standard screened Born formulas, soft-photon radiative corrections, and static-limit alternatives define materially different asymptotic regimes rather than small perturbations of a unique closed form.

5. Monte Carlo sampling, PIC integration, and positron-source design

In weighted particle-in-cell simulations, Bethe-Heitler events can be implemented as local pairwise interactions between photon and ion macro-particles. A representative algorithm is: pair the photon with a randomly chosen ion macro-particle of weight $Z^2$4; compute local densities $Z^2$5; evaluate the total cross section $Z^2$6 in the ion rest frame; draw an interaction with probability

$Z^2$7

if an event occurs, sample the positron Lorentz factor from the normalized differential cross section, create $Z^2$8 and $Z^2$9 macro-particles along the photon direction, remove the photon, and Lorentz-transform back to the simulation frame. Energy-momentum conservation is enforced by replacing a photon of energy $F(q)$0 with a pair of total energy $F(q)$1, neglecting recoil. The Monte Carlo timestep must satisfy $F(q)$2 (Martinez et al., 2019).

One-dimensional PIC studies of fast-electron transport through micrometric Cu foils show that the Bethe-Heitler yield $F(q)$3 grows roughly as $F(q)$4 for $F(q)$5, and as $F(q)$6 for larger $F(q)$7. In the example $F(q)$8, $F(q)$9, the PIC yields are $p+\gamma \rightarrow p+e^+ + e^-.$0 and $p+\gamma \rightarrow p+e^+ + e^-.$1 per kJ, whereas the corresponding Myatt 0D values are $p+\gamma \rightarrow p+e^+ + e^-.$2 and $p+\gamma \rightarrow p+e^+ + e^-.$3. The analysis attributes the discrepancy to electron energy losses from plasma expansion, with ion acceleration draining $p+\gamma \rightarrow p+e^+ + e^-.$4–$p+\gamma \rightarrow p+e^+ + e^-.$5 of the electron energy into ion kinetic energy over $p+\gamma \rightarrow p+e^+ + e^-.$6 (Martinez et al., 2019).

A recent compression of the same PIC workflow replaces pre-calculated Bethe-Heitler tables by a neural-network surrogate. The network is trained to predict Bethe-Heitler pair production cross sections for atomic numbers $p+\gamma \rightarrow p+e^+ + e^-.$7–$p+\gamma \rightarrow p+e^+ + e^-.$8 and photon energies between $p+\gamma \rightarrow p+e^+ + e^-.$9 and $p+\gamma \rightarrow p+e^+ + e^-$00 in the PIC code OSIRIS. The reported model is as accurate as pre-calculated tables and requires a hundred times less memory to store. It is first validated against a theoretical estimate in a simplified context and then shown to have similar performance in a typical relativistic laser-plasma interaction scenario (Amaro et al., 2024).

In all-optical positron-source design, Bethe-Heitler pairs generated in a foil can be trapped and accelerated by direct laser acceleration in a plasma channel. One reported configuration achieves an 8-fold increase in positron retention compared to previous studies. For a $p+\gamma \rightarrow p+e^+ + e^-$01 laser with $p+\gamma \rightarrow p+e^+ + e^-$02, $p+\gamma \rightarrow p+e^+ + e^-$03, a $p+\gamma \rightarrow p+e^+ + e^-$04 Al foil at $p+\gamma \rightarrow p+e^+ + e^-$05, the reported numbers are $p+\gamma \rightarrow p+e^+ + e^-$06 and $p+\gamma \rightarrow p+e^+ + e^-$07 after $p+\gamma \rightarrow p+e^+ + e^-$08 propagation; the accelerated spectrum is quasi-thermal up to $p+\gamma \rightarrow p+e^+ + e^-$09, the mean energy gain is $p+\gamma \rightarrow p+e^+ + e^-$10, and the beam divergence narrows from $p+\gamma \rightarrow p+e^+ + e^-$11 to $p+\gamma \rightarrow p+e^+ + e^-$12 (Gamiz et al., 2024).

6. Astrophysical and experimental roles

Domain	Bethe-Heitler role	Reported signature
Blazars	$p+\gamma \rightarrow p+e^+ + e^-$13 photopair injection	“pe bump” at $p+\gamma \rightarrow p+e^+ + e^-$14–$p+\gamma \rightarrow p+e^+ + e^-$15
GRBs	Secondary-pair synchrotron in proton-synchrotron models	subdominant tail to few tens or hundreds of MeV, or $p+\gamma \rightarrow p+e^+ + e^-$16 regime
ULXs	Secondary pairs in super-Eddington funnels	$p+\gamma \rightarrow p+e^+ + e^-$17–$p+\gamma \rightarrow p+e^+ + e^-$18 emission
M87* flares	BH-triggered electromagnetic cascade	UV-X-ray hump and suppressed GeV tail
MAGIX@MESA	Dominant Standard-Model background and BSM production mode	reach to $p+\gamma \rightarrow p+e^+ + e^-$19 couplings

In leptohadronic blazar models, the process $p+\gamma \rightarrow p+e^+ + e^-$20 injects broad secondary-pair distributions whose synchrotron emission fills the gap between the low- and high-energy humps of the spectral energy distribution. One study predicts a “pe bump” in the $p+\gamma \rightarrow p+e^+ + e^-$21–$p+\gamma \rightarrow p+e^+ + e^-$22 range, alongside $p+\gamma \rightarrow p+e^+ + e^-$23–$p+\gamma \rightarrow p+e^+ + e^-$24 neutrino emission, and derives an approximate analytical expression for the photopair loss rate that is accurate to a few percent for $p+\gamma \rightarrow p+e^+ + e^-$25 and $p+\gamma \rightarrow p+e^+ + e^-$26. A later study provides an empirical monoenergetic kernel for the injected pair spectrum that reproduces the shape of $p+\gamma \rightarrow p+e^+ + e^-$27 and $p+\gamma \rightarrow p+e^+ + e^-$28 to better than $p+\gamma \rightarrow p+e^+ + e^-$29–$p+\gamma \rightarrow p+e^+ + e^-$30 over five orders of magnitude in $p+\gamma \rightarrow p+e^+ + e^-$31, for $p+\gamma \rightarrow p+e^+ + e^-$32, and concludes that $p+\gamma \rightarrow p+e^+ + e^-$33-rays in low- and intermediate-peaked blazars may arise from Bethe-Heitler pairs in regions of the jet with typical transverse size $p+\gamma \rightarrow p+e^+ + e^-$34 and co-moving magnetic field $p+\gamma \rightarrow p+e^+ + e^-$35–$p+\gamma \rightarrow p+e^+ + e^-$36 (Petropoulou et al., 2014, Karavola et al., 2024).

In proton-synchrotron models of gamma-ray bursts, the process is labeled BeHe and constrains the allowed parameter space. Two regimes are reported. At high bulk Lorentz factor, large radius, and low luminosity, proton synchrotron emission dominates and a subdominant power law from BeHe pairs extends up to few tens or hundreds of MeV. At low bulk Lorentz factor, small radius, and high luminosity, BeHe cooling dominates and the spectrum becomes a single power law with spectral index $p+\gamma \rightarrow p+e^+ + e^-$37 across the entire GBM/Swift window, which is stated to be incompatible with observations (Bégué et al., 2021).

In ultraluminous X-ray sources, relativistic protons interacting with ambient photons in the funnel produce secondary pairs through the Bethe-Heitler channel. Using self-similar accretion-disk models with strong winds, one study reports nonthermal radiation in the $p+\gamma \rightarrow p+e^+ + e^-$38–$p+\gamma \rightarrow p+e^+ + e^-$39 band with luminosities from $p+\gamma \rightarrow p+e^+ + e^-$40 up to $p+\gamma \rightarrow p+e^+ + e^-$41. In the modeled ULX funnels, Bethe-Heitler pairs outnumber primary electrons by $p+\gamma \rightarrow p+e^+ + e^-$42, and the predicted MeV excess is identified as a possible diagnostic of relativistic protons in super-Eddington accretion flows (Romero et al., 3 Sep 2025).

For reconnection-driven flares in M87*, disk photons act as targets for Bethe-Heitler pair production by accelerated protons. The resulting pairs emit very high-energy synchrotron photons, $p+\gamma \rightarrow p+e^+ + e^-$43, which are then attenuated by the disk photon field and feed an electromagnetic cascade. In the baseline model summarized in the literature, the effective fraction of proton power lost to Bethe-Heitler pairs is $p+\gamma \rightarrow p+e^+ + e^-$44, the synchrotron-cascade hump peaks at $p+\gamma \rightarrow p+e^+ + e^-$45–$p+\gamma \rightarrow p+e^+ + e^-$46 between $p+\gamma \rightarrow p+e^+ + e^-$47 and $p+\gamma \rightarrow p+e^+ + e^-$48, and the escaping $p+\gamma \rightarrow p+e^+ + e^-$49 luminosity is suppressed by a factor $p+\gamma \rightarrow p+e^+ + e^-$50 relative to pure proton-synchrotron expectations (Petropoulou et al., 26 Feb 2026).

In fixed-target spectroscopy, Bethe-Heitler production is both signal and background. For MAGIX@MESA, the dominant Standard-Model background to resonant $p+\gamma \rightarrow p+e^+ + e^-$51 signals in the few-to-hundred-MeV regime is Bethe-Heitler pair production off the nuclear Coulomb field in

$p+\gamma \rightarrow p+e^+ + e^-$52

With $p+\gamma \rightarrow p+e^+ + e^-$53 and $p+\gamma \rightarrow p+e^+ + e^-$54, $p+\gamma \rightarrow p+e^+ + e^-$55, and two weeks of continuous running corresponding to $p+\gamma \rightarrow p+e^+ + e^-$56, the projected reach extends to mediator-electron couplings down to $p+\gamma \rightarrow p+e^+ + e^-$57. Specific numbers reported are $p+\gamma \rightarrow p+e^+ + e^-$58 for a vector mediator near $p+\gamma \rightarrow p+e^+ + e^-$59 in Setup I at $p+\gamma \rightarrow p+e^+ + e^-$60, $p+\gamma \rightarrow p+e^+ + e^-$61 in Setup I at $p+\gamma \rightarrow p+e^+ + e^-$62, and a best reach $p+\gamma \rightarrow p+e^+ + e^-$63 near $p+\gamma \rightarrow p+e^+ + e^-$64 in Setup II at $p+\gamma \rightarrow p+e^+ + e^-$65 (Pustyntsev et al., 18 May 2026).

A plausible implication of these diverse applications is that Bethe-Heitler pair production is best understood not as a narrow textbook process but as a family of closely related QED channels whose practical meaning depends on the target field, screening model, radiation environment, and observable of interest.