Parametrically driven nonlinear Dirac equation with arbitrary nonlinearity

TL;DR

This paper investigates how the nonlinearity parameter $ppa$ influences the transition between trapped and unbound soliton behavior in a parametrically driven nonlinear Dirac equation, using numerical and variational methods.

Contribution

It extends previous work by analyzing the $ppa$ dependence of soliton dynamics and transition regimes with a variational approximation.

Findings

Transition from trapped to unbound behavior depends on $ppa$, $r$, and $K$.

Higher $ppa$ extends the trapped regime for solitons.

The variational approximation accurately captures the low-order moment dynamics.

Abstract

The damped and parametrically driven nonlinear Dirac equation with arbitrary nonlinearity parameter $\kappa$ is analyzed, when the external force is periodic in space and given by $f(x) =r\cos(K x)$, both numerically and in a variational approximation using five collective coordinates (time dependent shape parameters of the wave function). Our variational approximation satisfies exactly the low-order moment equations. Because of competition between the spatial period of the external force $\lambda=2 \pi/K$, and the soliton width $l_s$, which is a function of the nonlinearity $\kappa$ as well as the initial frequency $\omega_0$ of the solitary wave, there is a transition (at fixed $\omega_0$) from trapped to unbound behavior of the soliton, which depends on the parameters $r$ and $K$ of the external force and the nonlinearity parameter $\kappa$. We previously studied this phenomena when $\kappa=1$ (2019 J. Phys. A: Math. Theor. {\bf 52} 285201) where we showed that for $\lambda \gg l_s$ the soliton oscillates in an effective potential, while for $\lambda \ll l_s$ it moves uniformly as a free particle. In this paper we focus on the $\kappa$ dependence of the transition from oscillatory to particle behavior and explicitly compare the curves of the transition regime found in the collective coordinate approximation as a function of $r$ and $K$ when $\kappa=1/2,1,2$ at fixed value of the frequency $\omega_0$. Since the solitary wave gets narrower for fixed $\omega_0$ as a function of $\kappa$, we expect and indeed find that the regime where the solitary wave is trapped is extended as we increase $\kappa$.