🌀 nonlinear sciences

Dirac Equations With All Possible Combinations of Scalar-Scalar, Vector-Vector and Pseudoscalar-Pseudoscalar Interactions

This paper presents exact solitary wave solutions for the nonlinear Dirac equation in 1+1 dimensions with all possible combinations of scalar, vector, and pseudoscalar interactions at nonlinearity parameter $\kappa=1$ , demonstrating that these models reduce to known forms, exhibit diverse wave profiles and frequency constraints, and universally yield bound states where the energy-to-charge ratio is independent of the coupling constant and less than the particle mass.

Original author: Avinash Khare, Fred Cooper, John F. Dawson, Avadh Saxena

Published 2026-11-02

📖 5 min read🧠 Deep dive

CC BY 4.0

Original author: Avinash Khare, Fred Cooper, John F. Dawson, Avadh Saxena

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the author. For technical accuracy, refer to the original paper. Read full disclaimer

In the subatomic world, particles like electrons are not merely tiny billiard balls; they are described by a set of rules that blend the motion of objects with the strange behavior of waves. This description, known as the Dirac equation, is a cornerstone of modern physics, explaining how these particles move at speeds close to light and how they possess an intrinsic property called spin. For decades, physicists have used this equation to understand everything from the structure of atoms to the behavior of light in special materials. However, the standard version of this equation assumes particles do not interact with one another in complex ways. To model real-world scenarios where particles influence each other—such as in certain exotic materials or theoretical models of how matter is held together—scientists must turn to a more complex version called the nonlinear Dirac equation. This advanced form allows for interactions where the particles' own presence alters the field they move through, creating self-sustaining waves that can hold their shape as they travel.

A team of researchers has now taken a significant step in understanding these complex interactions by solving a specific, highly intricate version of this equation. They focused on a scenario where the strength of the interaction is set to a particular value, a condition that corresponds to the most general type of four-particle interaction allowed by the laws of physics. In this setting, the particles can interact through three distinct mechanisms: a scalar interaction that affects their mass-like properties, a vector interaction that relates to their flow and momentum, and a pseudoscalar interaction, which is a more subtle, mirror-image version of the scalar type. The researchers asked a fundamental question: if you mix these three types of interactions in every possible combination, do stable, self-contained waves still exist, and what do they look like?

The answer they found is both surprising and unifying. Despite the apparent complexity of mixing these three different forces, the researchers discovered that at this specific interaction strength, every possible combination of these forces can be mathematically reduced to a simpler mix of just the scalar and vector types. This means that the universe of possibilities, which initially seemed to contain ten distinct and complicated models, actually collapses into a smaller set of known behaviors. It is as if a complex, multi-colored prism, when viewed from a specific angle, reveals that it is made of only two primary colors. This reduction allowed the team to apply existing knowledge to the new, more complex scenarios, confirming that stable, self-sustaining waves do indeed exist in all ten combinations.

These waves, known as solitary waves or solitons, are remarkable because they do not spread out or dissipate over time. Instead, they travel as a single, coherent packet of energy. The researchers calculated the energy and the total charge of these waves for every possible combination of interactions. They found that while the specific amount of energy and charge depends on the strength of the forces involved, the ratio between the two remains constant regardless of that strength. More importantly, they confirmed that in every case, the energy of the wave is less than what would be required to create a free particle of the same mass. This proves that these waves are bound states; they are tightly held together, unable to break apart into individual particles without an external push.

The shape of these waves also revealed a rich variety of behaviors depending on the specific mix of forces and the frequency of the wave. In some cases, the wave appears as a single, smooth hill of probability. In other cases, the wave splits into a double-humped shape, resembling two peaks side by side. The researchers mapped out exactly when this transition occurs. They found that for certain combinations of forces, the wave starts as a double-humped shape at lower frequencies and smoothly transforms into a single-humped shape as the frequency increases. There is a precise tipping point where the wave becomes flat-topped before settling into its single peak. This transition is not random; it is dictated by the specific parameters of the interaction, allowing scientists to predict the exact shape of the wave based on the forces at play.

The study also identified strict limits on where these waves can exist. For some combinations of forces, stable waves can form across the entire range of possible frequencies. For others, the waves can only exist if the frequency is high enough, falling within a specific window between a minimum value and the maximum possible limit. If the frequency drops below this minimum, the wave cannot sustain itself and disappears. This finding is crucial because it defines the boundaries of stability for these exotic states of matter. The researchers also noted that the stability of these waves, which is essential for them to persist in the real world, is linked to the slope of their energy-to-charge ratio. Since this slope is always positive in their models, the waves are predicted to be stable against small disturbances, at least in the non-relativistic regime.

While the mathematical work is complete for this specific interaction strength, the researchers acknowledge that the story is not finished. They point out that understanding whether these waves remain stable when subjected to external fields or when they collide with one another remains a difficult challenge. Furthermore, they have only explored the case where the interaction strength is set to one. They suspect that if this strength is changed to other values, the behavior of the waves could become even more complex, potentially allowing for waves with three, four, or more humps. They also wonder if these theoretical models have a direct connection to physical phenomena observed in the laboratory, such as in the behavior of super-cooled gases known as Bose-Einstein condensates. For now, this work provides a complete and unified map of the landscape for these specific interactions, showing that even in a world of complex forces, nature often finds a way to simplify the rules.

Technical Summary: Nonlinear Dirac Equations With All Possible Combinations of Scalar-Scalar, Vector-Vector and Pseudoscalar-Pseudoscalar Interactions

Problem Statement
The paper addresses the search for exact solitary wave solutions to the nonlinear Dirac (NLD) equation in 1+1 dimensions. While previous studies have extensively analyzed NLD models with scalar-scalar (SS) and vector-vector (VV) interactions, and more recently mixed SS-VV models, the most general NLD models can include an admixture of pseudoscalar-pseudoscalar (PP) interactions. In 1+1 dimensions, pseudovector-pseudovector interactions are equivalent to VV, and tensor interactions relate to pseudoscalar ones, making the inclusion of SS, VV, and PP the most general case. The specific problem tackled is to obtain exact solutions of the form $\psi(x, t) = e^{-i\omega t}\psi(x)$ for all possible combinations of these three interaction types at the nonlinearity parameter $\kappa = 1$ , which corresponds to the most general 4-Fermion interactions. The study focuses on the regime $0 < \omega < m$ , where $\omega$ is the frequency and $m$ is the mass.

Methodology
The authors utilize a Lagrangian density formulation for the SS+VV+PP model, characterized by three parameters $q, p, r$ governing the strength of the respective interactions. The field equation is derived and solved using a specific representation of the $\gamma$ matrices. The authors seek stationary solutions in the rest frame, parameterizing the spinor components using an amplitude $R(x)$ and a phase angle $\theta(x)$ .

Key steps in the methodology include:

Reduction to Ordinary Differential Equations: By applying current conservation and stress-energy tensor conservation, the authors derive a differential equation for the phase $\theta(x)$ and an algebraic relation for the amplitude squared $R^2$ .
Transformation of Interactions: A central methodological insight is the demonstration that at $\kappa = 1$ , the PP interaction term in the denominator of the $R^2$ expression can be algebraically rewritten as a linear combination of SS and VV terms. This allows the reduction of any of the ten possible SS-VV-PP combinations into equivalent models containing only SS and VV interactions with modified coupling parameters ( $\tilde{q}, \tilde{p}$ ).
Analysis of Solutions: For each of the ten combinations, the authors analyze the existence conditions for the charge $Q$ and energy $E$ . They determine the range of $\omega$ for which solutions exist, specifically looking for minimum frequency thresholds ( $\omega_{min}$ ) where the denominator of the integral for $Q$ does not vanish.
Hump Analysis: The second derivative of $R^2$ at $x=0$ is calculated to distinguish between single-hump and double-hump solitary wave solutions. A critical frequency $\omega_0$ is identified where the transition from double-hump to single-hump occurs.
Stability Indicators: The ratio of energy to charge ( $E/Q$ ) is computed. The authors note that for $\kappa=1$ , $E/Q < m$ (indicating bound states) and the slope of $E/Q$ with respect to $\omega$ is positive, which is associated with linear stability in the non-relativistic regime.

Key Contributions and Results
The paper establishes that at $\kappa = 1$ , the most general NLD model with arbitrary SS, VV, and PP interactions is mathematically equivalent to a suitable combination of SS and VV interactions. Consequently, the known properties of SS-VV models apply to the general case. The specific results for the ten possible combinations are:

Equivalence to Known Models: Two combinations are equivalent to $SS - (1/\tilde{p})VV$ , and two are equivalent to $VV - (1/\tilde{q})SS$ . The remaining six combinations reduce to various forms of $SS \pm VV$ depending on the specific values of the parameters $p, q, r$ .
Existence of Solutions:
- For some combinations (e.g., pure SS, pure VV, or specific mixed cases), solutions exist for the entire range $0 < \omega < m$ .
- For other combinations (specifically those equivalent to $SS - VV$ or $VV - SS$ with certain parameter constraints), solutions only exist for $\omega_{min} < \omega < m$ , where $\omega_{min}$ is determined by the coupling parameters.
Solitary Wave Profiles:
- Depending on the parameters and the value of $\omega$ , solitary waves can exhibit either a single-hump or a double-hump profile.
- A critical frequency $\omega_0$ is derived for each case, marking the transition where the solitary wave becomes flat-topped before splitting into a single hump as $\omega$ increases.
Charge and Energy: Explicit analytical expressions for the charge $Q$ $Q$ and energy $E$ $E$ are derived for all ten cases.
- While $E$ and $Q$ individually depend on the coupling constant $g$ , their ratio $E/Q$ is independent of $g$ .
- The ratio satisfies $E/Q > \omega$ and, crucially, $E/Q < m$ for all cases at $\kappa=1$ . This confirms that all solitary waves in these models are bound states.
- The positive slope of $E/Q$ with respect to $\omega$ suggests linear stability in the non-relativistic regime.

Significance and Claims
The paper claims that its primary significance lies in providing the first systematic classification and exact solution set for NLD models containing arbitrary admixtures of SS, VV, and PP interactions at $\kappa = 1$ . By proving the reducibility of these general models to the well-studied SS-VV sector, the authors conclude that the property of being a bound state ( $E/Q < m$ ) holds universally for these interactions at this nonlinearity parameter.

The authors modestly frame the work as a "first step" in this direction. They explicitly state that while they have determined the existence and profile of these solutions, the problem of demonstrating full spectral stability remains difficult and unresolved. They identify several open questions for future work, including:

Investigating the link between spectral stability and the single/double-hump nature of the solutions.
Examining the behavior of these bound states under external fields.
Finding physical relevance for these mixed-interaction models in systems like Bose-Einstein condensates.
Extending the analysis to arbitrary integer values of $\kappa$ (noting differences between odd and even $\kappa$ regarding parity symmetry and the potential for multi-humped solutions), which is addressed in a companion paper.

The paper does not propose new experimental setups or claim immediate physical applications but rather provides the theoretical foundation and exact solutions necessary to explore such applications.

More like this