Soliton dynamics and stability in the ABS spinor model with a PT-symmetric periodic potential

Basic Information

• Paper ID: 2511.05167
• Title: Soliton dynamics and stability in the ABS spinor model with a PT-symmetric periodic potential
• Authors: Franz G. Mertens, Bernardo Sánchez-Rey, Niurka R. Quintero
• Classification: nlin.PS (Nonlinear Sciences - Pattern Formation and Solitons)
• Publication Date: November 7, 2025
• Paper Link: https://arxiv.org/abs/2511.05167

Abstract

This study investigates the effects of introducing PT-symmetric complex potentials on soliton dynamics in the ABS model, a specific family of the (1+1)-dimensional cubic Dirac equation. By exploiting the Lagrangian formalism inherent to nonlinear Dirac equations, potentials are introduced whose imaginary parts relate to gain and loss, manifesting as spatially periodic damping acting simultaneously on both spinor components. Through the development of collective coordinate theory, the position, rapidity, momentum, frequency, and phase of moving solitons are all treated as functions of time. Treating the complex potential as a perturbation, the consistency between numerical solutions of the collective coordinate equations of motion and simulations of the nonlinear Dirac equation is verified. The primary effect of the imaginary potential is to induce oscillations in charge and energy with the same frequency and phase as the momentum. Long-lived solitons are found even under conditions of large charge and energy oscillations. Furthermore, empirical stability criteria previously successfully applied in nonlinear Schrödinger equations are extended to nonlinear Dirac equations.

Research Background and Motivation

Problem Background

1. Rise of PT-symmetric systems: Since the pioneering work of Bender and Boettcher, PT-symmetric non-Hermitian Hamiltonians have attracted widespread attention due to their real eigenvalue properties
2. Optical realization: Optical lattices provide excellent testing platforms for PT-symmetric systems, as Schrödinger equations with complex potentials are formally identical to paraxial optical wave equations
3. Soliton stability: PT-symmetric nonlinear optical lattices support soliton solutions that exhibit stability over broad parameter ranges

Research Motivation

1. Theoretical extension: Extending the concept of PT symmetry from nonlinear Schrödinger equations to nonlinear Dirac equations
2. Novel implementation: Proposing a new method to introduce PT-symmetric terms in the ABS model by exploiting its Lagrangian formalism
3. Stability analysis: Verifying and extending empirical stability criteria to the nonlinear Dirac equation domain

Limitations of Existing Methods

1. Limited conventional approaches for implementing PT-symmetric Dirac equations
2. Lack of systematic theoretical analysis of soliton dynamics under complex potentials
3. Applicability of stability criteria in nonlinear Dirac equations remains unverified

Core Contributions

1. Novel PT-symmetric implementation: Proposes a new method for introducing complex potentials in the ABS model through Lagrangian density and dissipation functions
2. Collective coordinate theory: Develops variational theory incorporating five collective coordinates (position, rapidity, momentum, frequency, phase)
3. Analytical approximate solutions: Obtains analytical approximate expressions for charge and energy in the non-relativistic limit
4. Stability criterion verification: First verification of empirical stability criteria in nonlinear Dirac equations
5. Long-lived soliton discovery: Discovers long-lived solitons even under conditions of large charge and energy oscillations

Methodology Details

Problem Definition

Investigating the dynamical behavior of ABS solitons under PT-symmetric complex potential $V(x) + iW(x)$, where:

• $V(x) = -V_0 \cos(kx)$ (real part, even function)
• $W(x) = -W_0 \sin(lx)$ (imaginary part, odd function)

ABS Model Architecture

The fundamental equations of the ABS model are:

i(u_t - u_x) + v + u*v² = 0
i(v_t + v_x) + u + v*u² = 0

where $u(x,t)$ and $v(x,t)$ are the two components of the Dirac spinor.

Complex Potential Introduction Method

Real Part Introduction

Real potential is introduced through gauge-covariant derivatives $i\partial_μΨ → (i\partial_μ - eA_μ)Ψ$:

L₃ = -V(x)[|u(x,t)|² + |v(x,t)|²]

Imaginary Part Introduction

Imaginary potential is introduced through the dissipation function:

F = -iW(x)[uu*_t - u*u_t + vv*_t - v*v_t]

Collective Coordinate Theory

Trial Wave Function

Adopting the same functional form as exact solutions, but with time-dependent parameters:

u(z) = e^(-β/2) a(z) e^(iθ(z)) e^(-iφ+ipz/cosh(β))
v(z) = -e^(β/2) a(z) e^(-iθ(z)) e^(-iφ+ipz/cosh(β))

where $z = \cosh[β(t)](x - q)$

Equations of Motion

Five equations of motion for collective coordinates are obtained through Lagrange-Euler equations:

1. Frequency equation: $\frac{dQ}{dω}ω̇ = -2W_0 I_9(β,ω)\sin(lq)$
2. Velocity-rapidity relation: $q̇ = \tanh(β) - \frac{2W_0\cos(lq)}{\cosh(β)Q(ω)}I_8(β,ω)$
3. Momentum equation: $Ṗ = -\frac{∂U}{∂q} - 2W_0\sin(lq)[pI_9(β,ω) + \sinh(β)I_7(β,ω)]$

Experimental Setup

Numerical Methods

1. Spatial discretization: $\Delta x = 0.02$, system length far exceeding soliton width
2. Time integration: Fifth-order Runge-Kutta-Verner algorithm with adaptive step size
3. Spatial derivatives: Spectral method computation
4. Boundary conditions: Periodic boundary conditions

Parameter Settings

• Frequency: $ω^{(0)} = 0.9, 0.74$
• Wave number: $k = π/32$, $l = k, 2k, k/2$
• Potential amplitude: $V_0 = 0.001$, $W_0 ∈ [0, 0.01]$
• Initial position: $q^{(0)} = 0$
• Initial phase: $φ^{(0)} = 0$

Evaluation Metrics

1. Soliton position: $q(t) = \frac{\int dx\,x\,ρ(x,t)}{\int dx\,ρ(x,t)}$
2. Momentum: $P = \int dx\,\frac{i}{2}[u*_x u - u_x u* + v*_x v - v_x v*]$
3. Charge: $Q = \int dx\,ρ(x,t) = \int dx[|u|² + |v|²]$
4. Energy: $E = \int dx\,T^{00}(x,t)$

Experimental Results

Main Results

Critical Velocity Analysis

The critical velocity $v_c$ exhibits non-monotonic behavior with respect to $W_0$, reaching a minimum at $W_0 = 0.0035$.

Soliton Dynamics Characteristics ($l = k$ case)

1. Oscillatory motion ($q̇^{(0)} < v_c$):
   • Momentum, charge, and energy oscillate in phase
   • Frequency given by zeroth-order approximation: $Ω_{osc} = \frac{π}{4}\frac{kq̇^{(0)}m_1}{K(1/m_1)}$
2. Unbounded motion ($q̇^{(0)} > v_c$):
   • Soliton continues moving while undergoing periodic oscillations
   • Frequency: $Ω_{unb} = \frac{π}{2}\frac{kq̇^{(0)}}{K(m_1)}$

Analytical Approximation Verification

In the non-relativistic limit, first-order correction results:

Q^{(1)}(t) = (2/kV₀)[P^{(0)}(t) - P^{(0)}]
E^{(1)}(t) ≈ ω^{(0)}Q^{(1)}(t)

Show excellent agreement with numerical simulations.

Stability Analysis

Empirical Stability Criterion

When the slope of normalized momentum $\tilde{P}(q̇) = P/Q$ is negative, i.e.:

∂P̃/∂q̇ < 0

the soliton becomes unstable.

Stability Regions

• Stable region: Parameter regions where $\tilde{P}(q̇)$ slope is positive
• Unstable region: Blue shaded regions in Figure 1 where $\tilde{P}(q̇)$ exhibits negative slope branches

Different Wave Number Cases ($l ≠ k$)

$l = 2k$ Case

• Charge and energy oscillations exhibit dual-frequency characteristics
• $\tilde{P}(q̇)$ curve maintains positive slope, soliton remains stable long-term

$l = k/2$ Case

• Oscillations exhibit single frequency
• $\tilde{P}(q̇)$ curve develops a cusp at certain point
• Despite large charge and energy oscillation amplitudes (dropping to 25% of initial values), soliton maintains long lifetime

Related Work

PT-symmetric System Development

1. Quantum mechanics foundation: Bender and Boettcher's pioneering work established PT-symmetric non-Hermitian Hamiltonian theory
2. Optical realization: Implementation of PT-symmetric systems through refractive index distribution and gain-loss element design
3. Soliton research: Discovery and stability studies of soliton solutions in PT-symmetric nonlinear optical lattices

Nonlinear Dirac Equations

1. Exact solutions: Exact soliton solutions in Thirring model, Gross-Neveu model, and ABS model
2. External field effects: Soliton dynamics under real potential fields
3. Parametric driving: Research on parametrically driven nonlinear Dirac equations

Stability Criteria

1. NLS equations: Establishment and application of empirical stability criteria in nonlinear Schrödinger equations
2. Extended applications: Successful applications in parametrically driven and generalized NLS equations
3. NLD equations: First verification of this criterion in nonlinear Dirac equations in this work

Conclusions and Discussion

Main Conclusions

1. Complex potential effects: Imaginary potential primarily induces charge and energy oscillations with the same frequency and phase as momentum
2. Long-lived solitons: Solitons persist long-term even under large-amplitude oscillations
3. Stability criterion: Empirical stability criterion is verified for the first time in nonlinear Dirac equations
4. Theory verification: Collective coordinate theory shows excellent agreement with numerical simulations

Limitations

1. Small parameter approximation: Analytical results valid only for small $W_0$ and non-relativistic limit
2. Specific potential forms: Only periodic combinations of cosine and sine potentials considered
3. Parameter range: Stability analysis limited to specific parameter regions

Future Directions

1. Complex potential forms: Exploring other complex potential forms, such as double PT-symmetric potentials
2. Reverse configurations: Investigating cases where imaginary potentials have opposite distributions in two spinor components
3. Higher-dimensional extension: Extending theory to higher dimensions
4. Experimental realization: Exploring implementation possibilities in actual physical systems

In-Depth Evaluation

Strengths

1. Theoretical innovation: Proposes novel method for introducing PT-symmetric complex potentials in ABS model through dissipation functions
2. Systematic analysis: Establishes complete collective coordinate theory framework with five time-dependent parameters
3. Analytical-numerical combination: Obtains analytical approximations in non-relativistic limit verified through numerical simulations
4. Stability verification: First verification of empirical stability criteria in nonlinear Dirac equations
5. Rich physical phenomena: Discovers diverse dynamical behaviors including in-phase oscillations and multi-frequency oscillations

Weaknesses

1. Parameter restrictions: Applicability of analytical results limited to small parameters and non-relativistic approximation
2. Single potential form: Only specific periodic potential function forms considered
3. Stability analysis depth: Physical understanding of instability mechanisms requires further investigation
4. Experimental verification absent: Lacks experimental verification in actual physical systems

Impact

1. Theoretical contribution: Provides new theoretical framework for PT-symmetric nonlinear Dirac equation research
2. Methodological value: Application of collective coordinate method in complex potential systems has universal significance
3. Stability theory: Extends theoretical tools for soliton stability analysis
4. Cross-disciplinary applications: Results potentially applicable in optics, condensed matter physics, and related fields

Applicable Scenarios

1. Theoretical physics research: PT-symmetric systems, soliton theory, nonlinear dynamics research
2. Optical applications: Nonlinear optics, optical soliton transmission, optical lattice design
3. Condensed matter physics: Spintronics, topological insulators, and related fields
4. Numerical methods: Provides reference for numerical studies of complex nonlinear systems

References

The paper cites 40 important references covering multiple aspects including PT-symmetric theory foundations, optical implementations, nonlinear Dirac equations, and soliton stability, reflecting the comprehensiveness and cutting-edge nature of the research.