An Analytic Prescription for $t$-channel Singularities

Basic Information

• Paper ID: 2505.10890
• Title: An Analytic Prescription for $t$-channel Singularities
• Authors: Kento Asai, Nagisa Hiroshima, Joe Sato, Ryusei Sato, Masaki J. S. Yang
• Classification: hep-ph (High Energy Physics - Phenomenology)
• Publication Date: November 21, 2025 (v3)
• Paper Link: https://arxiv.org/abs/2505.10890
• Institutions: Yukawa Institute for Theoretical Physics (Kyoto University), Institute for Cosmic Ray Research (University of Tokyo), Yokohama National University, RIKEN, Saitama University

Abstract

Singularities in $t$-channel scattering processes represent divergences in scattering amplitudes that occur when stable particles become on-shell states during $t$-channel propagation. Such situations arise both in collider experiments and in cosmological particle production backgrounds. Currently, no universally applicable treatment scheme exists. This paper proposes a novel analytic method to identify and remove divergence sources. The scheme is completely analytical and expected to have broad applications, providing valuable tools for this research field.

Research Background and Motivation

Core Problem

In quantum field theory scattering processes, singularities arise when propagators become on-shell, causing scattering amplitudes to diverge. According to Feynman diagram topology, singularities fall into two categories:

1. $s$-channel singularities: Extensively studied in baryogenesis and leptogenesis, typically regularized through finite decay widths of unstable particles (narrow width approximation)
2. $t$-channel singularities: When intermediate states are stable particles (e.g., neutrinos), regularization via decay width is impossible—this is the core focus of this work

Importance of the Problem

$t$-channel singularities appear in multiple physical systems:

• Collider experiments: Such as $\mu^+\mu^- \to W^+e^-\bar{\nu}_e$ in muon colliders, where propagating $\nu_\mu$ can become a real particle
• Cosmological particle production: In the early universe when multiple unstable particles share the same daughter particles, $t$-channel singularity effects become more pronounced
• Double-counting problem: Scattering processes actually contain contributions from decay and inverse decay, requiring correct separation to avoid repetitive counting in Boltzmann equations

Limitations of Existing Methods

Several regularization schemes proposed in the literature have limitations:

1. Finite beam width effects: Only applicable to collider experiments, cannot be applied to cosmological systems
2. Imaginary parts of initial particle masses: Violates energy-momentum conservation at vertices
3. Background particle interactions: Only applicable at high energies (finite temperature effects), cannot study low-temperature cosmology
4. Numerical methods: The authors' previous work 23 used numerical methods to identify $1/\epsilon$ terms, but with high computational cost

Research Motivation

The goal of this work is to develop a completely analytical, universally applicable scheme for handling $t$-channel singularities, significantly reducing computational cost and enabling exploration of vast parameter spaces in cosmological systems with $t$-channel singularities.

Core Contributions

1. Proposes a completely analytical $t$-channel singularity treatment scheme: Through variable substitution and asymptotic expansion, analytically identifies and removes divergent terms
2. Clear physical interpretation: Explicitly associates divergent terms $1/\epsilon$ with decay/inverse-decay processes, achieving rigorous separation of scattering and decay contributions
3. Avoids double-counting: Correctly handles particle production processes in Boltzmann equations without introducing non-physical parameters
4. Demonstrates practical application: Applies the scheme to the $U(1)_{L_\mu-L_\tau}\times U(1)_L$ model, finding that scattering processes can dominate Majoron production at high temperatures
5. Broad applicability: The scheme applies to both collider physics and cosmology, providing tools for multiple research fields

Detailed Methodology

Task Definition

Given a scattering process containing a $t$-channel propagator, when the intermediate state particle satisfies the on-shell condition:

• Input: Squared scattering amplitude $|M^t|^2$ containing terms like $\frac{1}{((p_3-p_2)^2-m^2)^2+\epsilon^2}$
• Output: The "pure scattering" part after removing on-shell contributions, for use in cosmological evolution or collider cross-section calculations
• Constraints: Preserve energy-momentum conservation, introduce no non-physical parameters, complete analyticity

Theoretical Foundation: Field Theory Interpretation

Using the simple scalar field model in Figure 1 as an example, the cross section for scattering process $\phi_3\phi_4 \to \phi_1\phi_2$ can be decomposed as:

$\sigma = \sigma_{\text{off-shell}} + \sigma_{\text{on-shell}}$

where the on-shell part can be explicitly written as:

$\sigma_{\text{on-shell}} = \frac{1}{\epsilon}E \cdot \Gamma(\phi_3 \to \Phi\phi_1) \cdot \sigma(\phi_4\Phi \to \phi_2)$

This shows that the $1/\epsilon$ divergent term corresponds to sequential decay and inverse-decay processes (Figure 2), not true scattering.

Analytical Treatment Scheme

Core Idea

For integrals containing $t$-channel propagators:

$\int_a^b dX \frac{1}{X^2+\epsilon^2}f(p_1,p_2,p_3,p_4)$

where $X=(p_3-p_2)^2-m^2$, singularities appear when $a<0<b$.

Variable Substitution

Introduce new variable $X=\epsilon\bar{X}$, the integral becomes:

$\int_{a/\epsilon}^{b/\epsilon} d(\epsilon\bar{X})\frac{1}{\epsilon^2}\frac{1}{\bar{X}^2+1}g(p_1,p_2,p_3,p_4)$

Key observation: The function $g$ as a Lorentz invariant can be expressed as a power series in $\bar{X}$. Extracting the zeroth-order term (divergence source):

$\int_{a/\epsilon}^{b/\epsilon} d(\epsilon\bar{X})\frac{1}{\epsilon^2}\frac{1}{\bar{X}^2+1} = \frac{1}{\epsilon}\left[\arctan\left(\frac{b}{\epsilon}\right)-\arctan\left(\frac{a}{\epsilon}\right)\right]$

Asymptotic Expansion

Using standard asymptotic expansion:

$\arctan\frac{x}{\epsilon} = \text{sgn}(x)\frac{\pi}{2} - \frac{\epsilon}{x} + O(\epsilon^3)$

Case-by-Case Analysis

Based on the signs of integration interval endpoints:

1. $ab>0$ (no singularity): $\frac{1}{\epsilon}\left(\pm\frac{\pi}{2}-\frac{\epsilon}{b} \mp \frac{\pi}{2}+\frac{\epsilon}{a}\right) = \frac{1}{a}-\frac{1}{b}+O(\epsilon^2)$ No divergence, no special treatment needed.
2. $ab<0$ (contains singularity): $\frac{1}{\epsilon}\left(-\frac{\pi}{2}-\frac{\epsilon}{b}-\frac{\pi}{2}+\frac{\epsilon}{a}\right) = -\frac{\pi}{\epsilon}+\frac{1}{a}-\frac{1}{b}+O(\epsilon^2)$ The first term $-\pi/\epsilon$ corresponds to decay/inverse-decay and should be removed. Retain: $\frac{1}{a}-\frac{1}{b}$
3. $ab=0$ (boundary case): If $b=0$: $\frac{1}{\epsilon}\left(0-\frac{\pi}{2}+\frac{\epsilon}{a}\right) = -\frac{\pi}{2\epsilon}+\frac{1}{a}+O(\epsilon^2)$ Remove $-\pi/(2\epsilon)$, retain $\frac{1}{a}$

Technical Innovations

1. Complete analyticity: No need for numerical integration to identify divergent terms; directly calculate using analytical formulas
2. Physical transparency: Clearly associate $1/\epsilon$ terms with physical processes (decay/inverse-decay)
3. Computational efficiency: Significantly reduced computational cost compared to numerical methods
4. Universality: Independent of specific models, applicable to any process with $t$-channel singularities
5. No artificial parameters: No introduction of beam width, temperature cutoff, or other non-physical parameters

Experimental Setup

Application Model: $U(1)_{L_\mu-L_\tau}\times U(1)_L$ Model

This model contains:

• $Z'$ boson: Corresponding to $U(1)_{L_\mu-L_\tau}$ gauge symmetry
• Majoron $\phi$: Pseudo-Nambu-Goldstone boson from global lepton number $U(1)_L$ symmetry breaking

Lagrangian: $\mathcal{L} = \mathcal{L}_{\text{SM}} + \mathcal{L}_{Z'} + \mathcal{L}_\phi$

where: $\mathcal{L}_{Z'} = -\frac{1}{4}Z'_{\rho\sigma}Z'^{\rho\sigma} + \frac{1}{2}m_{Z'}^2Z'_\rho Z'^\rho + g_{Z'}Z'_\rho J^\rho_{Z'} + \epsilon eZ'_\mu J^\mu_{\text{EM}}$

$\mathcal{L}_\phi = -\frac{1}{2}m_\phi^2\phi^2 + (h_{\alpha\beta}\bar{\nu}_\alpha\nu_\beta^c\phi + \text{H.c.})$

Research Process

Examines Compton-type scattering process $Z'\nu_\alpha \leftrightarrow \phi\bar{\nu}_\beta$ (Figure 3), where:

• $s$-channel: No singularities
• $t$-channel: Intermediate neutrino can become on-shell, producing singularities

Physical Quantity Calculations

Squared Scattering Amplitude

Complete expression given in equation (24), containing three terms:

• First term: $s$-channel contribution
• Second term: $t$-channel contribution (contains singularity)
• Third term: Interference term

The $t$-channel term's denominator contains $(p_3-p_2)^4$ and $(p_3-p_2)^2$, which are singularity sources.

Integral Treatment

Introduce variable $x=(p_3-p_2)^2$, integration interval is $[m_{Z'}^2-2|p_2|\sqrt{s}, m_{Z'}^2-2|p_2|\frac{m_{Z'}^2}{\sqrt{s}}]$

Applying the scheme in this paper, the $t$-channel integral $I_t(\epsilon)$ (equations 37-38) decomposes into:

• Divergent term: $\frac{1}{\epsilon}\tan^{-1}\frac{x}{\epsilon}$ term
• Finite term: Logarithmic and algebraic terms

After removing divergences, obtain $I_0$ (equation 39), which is the pure scattering contribution.

Cosmological Application

Boltzmann Equation

Majoron number density evolution: $\frac{dn_a}{dt} = -3Hn_a + \frac{\delta n_a}{\delta t}\bigg|_S + \frac{\delta n_a}{\delta t}\bigg|_{\text{ID}}$

where scattering and inverse-decay terms are calculated separately, using singularity-removed amplitudes to avoid double-counting.

Interaction Rate

Define thermal equilibrium interaction rate: $\langle\Gamma_{S,\text{ID}}\rangle = \frac{1}{n_{\phi,\text{eq}}}\frac{\delta n_\phi}{\delta t}\bigg|_{S,\text{ID}}$

Determine whether processes are in thermal equilibrium via Gamow criterion $\langle\Gamma\rangle/H>1$.

Parameter Settings

• Majoron mass: $m_\phi = 0.1$ MeV
• $Z'$ boson mass: $m_{Z'} = 13$ MeV
• Majoron-neutrino coupling: $h_{\alpha\beta} = 10^{-8}, 10^{-9}$ (assuming all couplings identical)
• Temperature range: $10^{-3}m_\phi < T < m_\phi$
• Discussion limited to region after QCD phase transition

Experimental Results

Main Results

Figure 4 shows comparison of production efficiency for scattering and inverse-decay processes:

Key Findings

1. Scattering dominates at high temperature: At $T \sim m_\phi$ (normalized temperature $\sim 1$), scattering rate can exceed inverse-decay rate
2. Coupling dependence:
   • For $h=10^{-8}$: Scattering rate $\langle\Gamma_S\rangle/H \sim 10^{-1}$, satisfies Gamow criterion ($>1$), process in thermal equilibrium
   • For $h=10^{-9}$: Scattering rate $\sim 10^{-3}$, does not satisfy criterion
3. Low-temperature behavior: As temperature decreases, inverse-decay gradually dominates

Physical Interpretation

• At high temperatures, phase space factors make scattering processes more efficient
• Different energy dependence of scattering cross section versus inverse-decay leads to different temperature-dependent behavior

Important Conclusions

Corrections to Previous Work

Previous assumption: Majoron primarily produced through inverse-decay, scattering contribution negligible

This work's results:

• Scattering contribution non-negligible when $h \gtrsim 10^{-8}$
• Neglecting scattering only valid for $h < 10^{-9}$
• This threshold stronger than previous work 27,28

Cosmological Significance

• Majoron production efficiency underestimated
• Constraints on dark matter candidates require reassessment
• Affects cosmological history of lepton number violation processes

Method Verification

Successfully achieved through analytical scheme:

1. Removed $t$-channel singularity divergences
2. Separated scattering and decay contributions
3. Obtained physically reasonable results (finite and causality-preserving)
4. Revealed previously overlooked important physical effects

Related Work

$s$-channel Singularity Treatment

• Breit-Wigner resonance 4: Classical treatment scheme
• Narrow width approximation 7,8: Regularization via $1/(p^2-m^2+im\Gamma)$
• Leptogenesis 6: Study of right-handed neutrinos as intermediate states

History of $t$-channel Singularities

• First identified 9: Peierls (1961) in hadronic process $\pi N^* \to \pi N^*$
• General proof 10: Brayshaw et al. (1978) proved universality when meson satisfies on-shell condition and process decomposes into two subprocesses

Existing Regularization Schemes

Collider Physics

• Finite beam effects 11-15: Melnikov-Serbo scheme, regularization via beam transverse size
• Phase space truncation 14,15: Dams-Kleiss treatment method
• Complex momentum method 16,17: Nowakowski-Pilaftsis and Karamitros-Pilaftsis schemes

Cosmological Applications

• Quasiparticle picture 18,19: Grzadkowski et al. and Iglicki, imaginary part from background gas interactions
• Virtual mass scheme 20: Ginzburg's method, but violates energy conservation
• Neglect approximation 21,22: Escudero et al., assumes scattering as subdominant contribution

This Work's Innovation

Compared to existing work:

1. Completely analytical: First completely analytical universal scheme
2. No external parameters: No introduction of beam width, temperature cutoff, etc.
3. Physical transparency: Clearly separates physical processes
4. Broad applicability: Simultaneously applicable to collider and cosmology
5. Computational efficiency: Significantly reduces numerical computation cost

Conclusions and Discussion

Main Conclusions

1. Establishment of analytical scheme: Proposes completely analytical method for handling $t$-channel singularities through variable substitution $X=\epsilon\bar{X}$ and asymptotic expansion, analytically identifying and removing $1/\epsilon$ divergent terms
2. Clear physical picture: Divergent term $1/\epsilon$ explicitly corresponds to real particle propagation (decay + inverse-decay sequence); after removal, obtains pure scattering contribution
3. Successful practical application: In $U(1)_{L_\mu-L_\tau}\times U(1)_L$ model, finds scattering processes dominate Majoron production at high temperature, correcting previous assumption of neglecting scattering
4. Strengthened constraints: Coupling threshold where scattering contribution becomes non-negligible is $h \sim 10^{-9}$, stronger than previous work

Limitations

Theoretical Limitations

1. Leading-order approximation: Scheme based on tree-level calculations; loop correction effects not considered
2. Finite temperature effects: High-temperature quantum field theory modifications (e.g., thermal masses) not included; authors acknowledge more refined treatment needed at high temperatures 18,19
3. Stable particle assumption: Scheme specifically targets stable intermediate states; transition region for extremely narrow-width particles not discussed in detail

Application Limitations

1. Parameter space exploration: Complete parameter space scan and phenomenological conclusions for $U(1)_{L_\mu-L_\tau}\times U(1)_L$ model left for future work
2. Other model verification: Although scheme is universal, demonstrated in only one specific model
3. Collider applications: Mentioned applicability to muon colliders but no specific calculations provided

Technical Details

1. Numerical implementation: Although analytical, practical calculation still requires numerical phase space integration
2. Higher-order terms: Contributions from $O(\epsilon^2)$ and higher orders may be non-zero in some cases
3. Multiple singularities: Cases with simultaneously occurring multiple $t$-channel singularities not discussed

Future Directions

Short-term Applications

1. Complete study of $U(1)_{L_\mu-L_\tau}\times U(1)_L$ model: Systematic parameter space scan, providing cosmological and astrophysical constraints
2. Other beyond-standard-model applications: Application to other models containing Majorons or light scalars
3. Muon collider physics: Detailed calculations of processes like $\mu^+\mu^- \to W^+e^-\bar{\nu}_e$, comparison with finite beam scheme

Medium-term Development

1. Finite temperature corrections: Systematically include thermal effects at high temperatures, quantitative comparison with refs. 18,19
2. Loop corrections: Study single-loop correction effects on $t$-channel singularities
3. Multi-particle processes: Generalization to $2\to 3$ or more complex processes

Long-term Perspectives

1. Automated tools: Develop software packages for automatic identification and treatment of $t$-channel singularities
2. Quantum field theory foundations: Deep understanding from renormalization group and effective field theory perspectives
3. Non-equilibrium field theory: Rigorous treatment in Schwinger-Keldysh formalism

In-Depth Evaluation

Strengths

Method Innovation (★★★★★)

1. Breakthrough progress: First to achieve complete analytical treatment of $t$-channel singularities, solving long-standing technical challenges
2. Mathematical elegance: Variable substitution $X=\epsilon\bar{X}$ is concise and effective, transforming complex limit problems into standard asymptotic expansion
3. Deep physical insight: Explicitly establishes correspondence between $1/\epsilon$ terms and decay/inverse-decay, not merely mathematical technique

Experimental Sufficiency (★★★★☆)

1. Practical application verification: Demonstrates scheme feasibility through $U(1)_{L_\mu-L_\tau}\times U(1)_L$ model
2. Reasonable physical results: Obtained scattering and inverse-decay rates consistent in magnitude and temperature dependence with physical intuition
3. New physics discovery: Reveals previously overlooked scattering-dominated region

Shortcoming: Demonstration in only one model; verification in more models would strengthen persuasiveness

Result Convincingness (★★★★★)

1. Resolves double-counting: Rigorously avoids repetitive counting in Boltzmann equations, crucial for cosmological calculations
2. Quantitative improvement: Provides specific coupling threshold ($h \sim 10^{-9}$), directly applicable to model constraints
3. Clear comparison with previous work: Explicitly indicates applicable ranges of previous assumptions

Writing Clarity (★★★★★)

1. Rigorous logic: From physical motivation → field theory interpretation → analytical scheme → practical application, well-structured
2. Complete technical details: Key formula derivations detailed (equations 1-18), reproducible
3. Intuitive physical picture: Figures 1-3 effectively aid understanding

Shortcomings

Method Limitations

1. Missing temperature effects: Acknowledges need for finite temperature field theory corrections at high temperatures, but provides no specific scheme or quantitative comparison
2. Limited applicability range: Focuses on stable intermediate states; transition region for narrow-width particles ($\Gamma \ll m$ but $\Gamma \neq 0$) insufficiently discussed
3. Higher-order term treatment: $O(\epsilon^2)$ terms neglected, but may be important in some precision calculations

Experimental Design

1. Single parameter choice: Figure 4 shows only $m_\phi=0.1$ MeV, $m_{Z'}=13$ MeV case, lacking parameter scan
2. Simplified coupling assumption: Assumes all $h_{\alpha\beta}$ identical; actual models may have hierarchical structure
3. Missing error analysis: No discussion of numerical integration precision and $\epsilon$ value effects

Analysis Depth

1. Collider applications undeveloped: Although mentioning muon collider applicability, provides no specific calculations or comparison with beam width scheme
2. Incomplete phenomenological conclusions: Complete cosmological constraints and observational predictions for $U(1)_{L_\mu-L_\tau}\times U(1)_L$ model left for future
3. Relationship with finite-temperature scheme: No quantitative comparison of which temperature region this scheme agrees with refs. 18,19

Impact

Contribution to Field (★★★★★)

1. Solves long-standing problem: $t$-channel singularities, identified since 1961, lacked universal treatment; this work fills the gap
2. Methodological significance: Provides directly applicable tool, will promote research in related fields
3. Interdisciplinary value: Simultaneously applicable to particle physics (colliders) and cosmology

Practical Value (★★★★☆)

1. Computational efficiency improvement: Significantly reduces cost compared to numerical methods, enabling large-scale parameter scans
2. Directly usable: Equations (13-18) provide analytical formulas for different cases, directly programmable
3. Avoids artificial parameters: No need for beam width, temperature cutoff, or other difficult-to-determine quantities

Potential limitation: High-temperature region still requires finite-temperature corrections

Reproducibility (★★★★★)

1. Complete formulas: Key derivation steps detailed, especially equations (1-18)
2. Clear parameters: Figure 4 calculation parameters clearly provided
3. Data availability: Authors commit to "Plot data available upon request"

Applicable Scenarios

Ideal Application Scenarios

1. Low-temperature cosmology: Early universe with $T \ll m_W$, production of light particles (Majorons, axions, dark photons)
2. Neutrino physics: Processes involving neutrinos as intermediate states, such as lepton-number-violating scattering
3. Beyond-standard-model searches: Models requiring precise $t$-channel contribution calculations

Applicable Conditions

• Intermediate state particles stable or extremely long-lived ($c\tau \gg$ system scale)
• Leading-order perturbation theory applicable
• Not too high temperature ($T \lesssim$ GeV, avoiding strong thermal effects)

Inapplicable Scenarios

• High-temperature plasma ($T \sim$ 100 GeV and above), requiring complete finite-temperature field theory
• Strong coupling systems where perturbation theory fails
• Intermediate states with significant width ($\Gamma/m \gtrsim 0.1$), where traditional methods suffice

Potential Extensions

1. Automated implementation: Integration into FeynCalc, MadGraph and similar tools
2. Machine learning assistance: Neural networks to accelerate phase space integration
3. Gravitational wave cosmology: Application to particle physics processes in gravitational wave production
4. Dark matter direct detection: Precise scattering cross-section calculations including $t$-channel contributions

Technical Highlights

Mathematical Technique Sophistication

The variable substitution $X=\epsilon\bar{X}$ appears simple yet is ingenious:

• Transforms singular limit $\epsilon \to 0$ into asymptotic behavior $\bar{X} \to \infty$
• Exploits known asymptotic expansion of $\arctan$ function, avoiding complex regularization
• Naturally separates $1/\epsilon$ term with clear physical meaning

Physical Intuition Manifestation

Equation (6) derivation demonstrates deep physical understanding: $\sigma_{\text{on-shell}} = \frac{1}{\epsilon}E \cdot \Gamma(\phi_3 \to \Phi\phi_1) \cdot \sigma(\phi_4\Phi \to \phi_2)$

This is not merely a mathematical result but reveals the essence of $t$-channel singularities: on-shell propagators correspond to real particle creation and annihilation.

Practical Considerations

The case-by-case analysis in equations (13-18) covers all possible scenarios ($ab>0, ab<0, ab=0$), demonstrating scheme completeness and practicality.

Selected References

Foundational Theory

• 1 Lehmann, Symanzik, Zimmermann (1955): LSZ reduction formula
• 2 Cutkosky (1960): Cutting rules
• 3 Mandelstam (1958): Mandelstam variables

$s$-channel Singularities

• 4 Breit-Wigner (1936): Resonance theory
• 7,8 Kolb-Wolfram, Giudice et al.: Narrow width approximation

$t$-channel Singularity History

• 9 Peierls (1961): First identification
• 10 Brayshaw et al. (1978): General proof
• 11-15 Melnikov-Serbo, Dams-Kleiss: Collider schemes

Cosmological Applications

• 18,19 Grzadkowski et al., Iglicki: Finite-temperature schemes
• 21,22 Escudero et al.: Scattering-neglecting studies
• 23 Authors' previous work: Numerical scheme

Model-Related

• 24-25 $U(1)_{L_\mu-L_\tau}\times U(1)_L$ model
• 27,28 Escudero-Witte, Escudero: Majoron constraints

Overall Assessment

This is a high-quality theoretical physics methodology paper solving long-standing technical challenges in particle physics and cosmology. The method shows strong innovation, mathematical rigor, and clear physical picture, with broad application prospects. The main contribution is providing the first completely analytical $t$-channel singularity treatment scheme, significantly improving calculation efficiency and reliability for related computations.

Recommendation Index: 9/10

Target Audience: Researchers working on beyond-standard-model physics, early universe cosmology, and collider physics. This method is expected to become a standard tool in related fields.