Decay of solutions of nonlinear Dirac equations

Basic Information

• Paper ID: 2503.05410
• Title: Decay of solutions of nonlinear Dirac equations
• Authors: Sebastian Herr, Christopher Maulén, Claudio Muñoz
• Classification: math.AP (Analysis of PDEs)
• Publication Date: October 13, 2025 (arXiv version v4)
• Paper Link: https://arxiv.org/abs/2503.05410

Abstract

This paper investigates the long-time behavior of solutions to a broad class of nonlinear Dirac-type equations, encompassing both small and large solutions. The results are categorized according to the 1D massless and massive cases, the 3D general case, and the n-dimensional general case. In the 1D massless case, it is proven that any globally defined solution converges to zero as time approaches infinity, with convergence occurring within spatial regions expanding at the rate $t \log^{-2} t$. This result requires no assumptions on the smallness of initial data or specific powers of the nonlinearity, thereby excluding the existence of standing wave breathers or solitary wave structures in this case. In the 1D massive case, solitary waves are known to exist. By introducing novel virial identities adapted to the special algebra of Dirac operators, it is shown that under certain "holomorphic" odd nonlinearities, globally defined small odd solutions decay to zero on compact spatial sets as time approaches infinity.

Research Background and Motivation

Problem Background

1. Importance of Dirac equations: Nonlinear Dirac equations are among the most important models in relativistic quantum mechanics, describing self-interactions of Dirac fermions with significant applications in high-energy physics and condensed matter physics.
2. Necessity of studying long-time behavior: Understanding the long-time asymptotic behavior of solutions to nonlinear partial differential equations is a central problem in mathematical physics, crucial for predicting the evolution of physical systems.
3. Limitations of existing research:
   • The energy functional lacks definite sign, making traditional energy methods difficult to apply directly
   • The Hamiltonian is unbounded both above and below, reflecting the existence of particles and antiparticles
   • Lack of effective decay estimation techniques applicable to Dirac equations
4. Existence of breathers and solitary waves: Whether Dirac equations admit non-decaying structures analogous to breather solutions in nonlinear Schrödinger equations remains an important open problem.

Core Contributions

1. Universal decay in 1D massless case: Proves that any global solution converges to zero in the region $I(t) = (-\frac{|t|}{\log^2|t|}, \frac{|t|}{\log^2|t|})$ without restrictions on initial data size or nonlinearity powers.
2. Novel virial identity technique: Introduces specialized virial identities for Dirac equations for the first time, avoiding traditional "squaring" techniques and performing estimates directly within the Dirac framework.
3. Conditional decay in 1D massive case: Under "harmonic" odd nonlinearities satisfying Cauchy-Riemann-type conditions, proves local decay of small odd solutions.
4. Extension to 3D case: Under $H^1$ boundedness conditions without parity assumptions, establishes decay properties for an important class of Dirac models.
5. $L^2$ decay in exterior light cone regions for higher dimensions: For n-dimensional cases, proves $L^2$ decay in the "exterior light cone" region $I_j^b(t) = \{x \in \mathbb{R}^n | |x_j| \geq (1+b)t\}$.

Detailed Methodology

Problem Formulation

Study the long-time asymptotic behavior of nonlinear Dirac equations: $i\partial_t\psi = H\psi + (m + V(x,\psi))\beta\psi$ where $H = -i\alpha \cdot \nabla$ is the Hamiltonian and $V$ is a (possibly) nonlinear potential function.

Core Technical Methods

1. Construction of Virial Identities

For n-dimensional Dirac equations, construct modified charge functionals: $I = \frac{1}{\mu(t)}\int \phi\left(\frac{x+\rho(t)}{\lambda(t)}\right)\psi^*(t,x)\psi(t,x)dx$

Through careful selection of weight functions $\phi$ and parameters $\mu(t), \lambda(t), \rho(t)$, obtain the key virial identity:

Theorem 3.1: For $\psi = u_1 + iu_2$, $\frac{d}{dt}I = -\frac{\mu'(t)}{\mu^2(t)}\int \phi\left(\frac{x+\rho(t)}{\lambda(t)}\right)[u_1^*u_1 + u_2^*u_2]$$+ \text{other terms involving gradients and mass}$

2. Special Treatment of 1D Case

For 1D systems: $i(\partial_t u + \partial_x u) + mv = \partial_u W_1(u,\bar{u},v,\bar{v})$$i(\partial_t v - \partial_x v) + mu = \partial_v W_2(u,\bar{u},v,\bar{v})$

Construct functionals $K = \int \phi(\frac{x}{\lambda(t)})(|u|^2 + |v|^2)$ and $J = \int \phi(\frac{x}{\lambda(t)})(|u|^2 - |v|^2)$.

Lemma 3.2: In the massless case ($m=0$): $\frac{d}{dt}K = \int \phi'(|u|^2 - |v|^2)$

3. Treatment of 3D Radial Case

For 3D radial solutions, use partial wave decomposition to reduce the 4×4 complex system to a 2×2 system: $i\partial_t\phi_1 = \left(\partial_r + \frac{2}{r}\right)\phi_2 + m\phi_1 - W_1$$i\partial_t\phi_2 = -\partial_r\phi_1 - m\phi_2 + W_2$

Construct four virial functionals $K_1, \tilde{K}_1, K_2, \tilde{K}_2$, selecting special weights $\phi = \frac{r^{3/2}}{1+r}$.

Technical Innovations

1. Direct Dirac estimates: Avoid traditional Klein-Gordon reduction techniques, working directly within the Dirac framework, which is technically more challenging but more fundamental.
2. Careful selection of weight functions:
   • 1D: $\phi(x) = \tanh(x)$ for the massless case
   • 3D: $\phi(r) = \frac{r^{3/2}}{1+r}$ balancing various nonlinear terms
3. Introduction of harmonic conditions: For the 1D massive case, introduce Cauchy-Riemann-type conditions: $\partial_a W_2 + \partial_c W_1 = \partial_b W_2 - \partial_d W_1 = 0$ etc., ensuring the preservation of odd solutions.

Main Theorems and Results

Theorem 1.1 (1D Massless Decay)

For any global solution $(u,v)$ of the massless 1D Dirac equation, $\lim_{t \to \infty}\|(u,v)(t)\|_{L^2(I(t))} = 0$ where $I(t) = (-\frac{|t|}{\log^2|t|}, \frac{|t|}{\log^2|t|})$.

Theorem 1.3 (1D Massive Conditional Decay)

For odd nonlinearities satisfying harmonic conditions, small global odd solutions satisfy: $\lim_{t \to \infty}\|(\psi_1,\psi_2)(t)\|_{(L^2 \cap L^\infty)(I)} = 0$ for any compact interval $I$.

Theorem 1.4 (3D Radial Decay)

For 3D radial Dirac equations, under the condition $\|W_j\| \leq C|\phi|^p$ ($p \geq 3$), small global solutions satisfy: $\lim_{t \to \infty}\|\phi(t)\|_{L^2(B(0,R))} = 0$ for any $R > 0$.

Theorem 1.5 (High-Dimensional Exterior Light Cone Decay)

For n-dimensional Dirac equations, any global $L^2$ solution decays in the exterior light cone region: $\lim_{t \to \infty}\int_{I_j^b(t)}\psi^*(t,x)\psi(t,x)dx = 0$

Proof Techniques and Key Estimates

1. Integrability Estimates

Through virial identities obtain key integrability: $\int_0^\infty \frac{1}{\lambda(t)}\int \text{sech}^2\left(\frac{x}{\lambda(t)}\right)(|u|^2 + |v|^2)dx dt < \infty$

2. Sequential Decay

Utilize integrability to establish the existence of a sequence $t_n \to \infty$ such that corresponding $L^2$ norms tend to zero.

3. Continuity Argument

Lift sequential decay to continuous decay through precise virial estimates.

4. Treatment of Nonlinear Terms

• 1D case: Use harmonic conditions to eliminate mixed terms
• 3D case: Employ Sobolev embedding and special properties of radial functions

Related Work

Linear Theory

• D'Ancona & Fanelli et al. established dispersive estimates for magnetic Dirac equations
• Kopylova proved weighted energy decay for 1D Dirac equations

Nonlinear Theory

• Escobedo & Vega: Well-posedness in 3D massive case
• Candy & Herr: Bilinear Fourier restriction methods
• Pelinovsky: Scattering theory for 1D case

Solitary Wave Theory

• Known existence of standing waves in Soler models
• Variational methods by Merle et al. for construction
• Spectral stability studies by Boussaïd & Comech

Conclusions and Discussion

Main Conclusions

1. Universality in massless case: The 1D massless Dirac equation admits no non-decaying localized solutions, contrasting with Klein-Gordon equations.
2. Effectiveness of virial method: The newly introduced virial identities provide powerful tools for studying Dirac equations.
3. Dimensional effects: Decay properties improve in higher dimensions, consistent with general dispersive theory expectations.

Limitations

1. Restrictions in 1D massive case: Requires special harmonic conditions, excluding classical integrable models.
2. Smallness assumption in 3D case: Solutions must be uniformly bounded and sufficiently small in $H^1 \cap L^\infty$.
3. Radial symmetry: 3D results apply only to partial wave subspaces.

Future Directions

1. Asymptotic stability: Study asymptotic stability of Dirac solitary waves using virial methods.
2. 2D case: Address more complex 2D Dirac systems.
3. General non-radial case: Extend to non-radial 3D situations.

In-Depth Evaluation

Strengths

1. Technical innovation: First systematic virial theory for Dirac equations, avoiding technical difficulties of traditional Klein-Gordon reduction.
2. Strength of results: The 1D massless result requires no smallness assumptions, which is quite rare in nonlinear PDEs.
3. Generality of methods: Virial techniques applicable to a broad class of Dirac-type models.
4. Mathematical rigor: Proofs are complete and rigorous with careful handling of technical details.

Weaknesses

1. Scope of applications: Both 1D massive and 3D cases require additional technical assumptions, limiting applicability.
2. Physical interpretation: Lack of clear physical motivation for technical assumptions such as harmonic conditions.
3. Optimality: Unclear whether decay regions and decay rates are optimal.

Impact

1. Theoretical contribution: Provides new analytical tools for nonlinear Dirac equation theory.
2. Methodological value: Virial techniques may apply to other relativistic equations.
3. Physical significance: Excludes breather existence in certain cases, relevant to understanding relativistic quantum field theory.

Application Scenarios

• Long-time dynamics in relativistic quantum mechanics
• Dirac material research in condensed matter physics
• Dispersive equation theory in mathematical physics
• Solitary wave stability analysis

References

The paper contains 49 references covering well-posedness theory for Dirac equations, solitary wave theory, dispersive estimates, and other related fields, providing a solid theoretical foundation for this research.