Threshold dynamics of stochastic SIRSW infectious disease model with multiparameter perturbation

1.   Introduction

Integrable evolution equations are very common and play a key role in nonlinear science, as shown in studies like ^[1,2,3]. Specifically, the Korteweg-de Vries (KdV) equation, which comes from fluid dynamics, explains how waves with slight dispersion and small amplitude evolve in media with quadratic nonlinearity ^[2,4]. Additionally, the nonlinear Schrödinger (NLS) equation describes the behavior of waves that are slightly nonlinear and quasi-monochromatic in media with cubic nonlinearity ^[3,5,6]. Moreover, the derivative nonlinear Schrödinger (DNLS) equation is used to describe Alfvén waves in plasmas with low-$\beta$ values, as well as large-amplitude waves that run parallel to the magnetic field ^[7,8,9]. These equations are mathematically integrable, meaning they can be solved exactly.

Recently, Ablowitz and Musslimani introduced the following parity time (PT)-symmetric nonlocal NLS equation and derived its explicit solutions^[10]:

where $\varepsilon = \pm 1$ denotes the focusing $(+)$ and defocusing $(-)$ nonlinearity, and the asterisk denotes the complex conjugate. The PT symmetry implies that the self-induced potential $V(x, t) = u(x, t) u^{*}(-x, t)$ fulfills the relation $V^{*}(-x, t) = V(x, t)$, while the nonlocality indicates that the value of the potential $V(x, t)$ at $x$ requires the information on $u(x, t)$ at $x$, as well as at $-x$^[11,12,13]. Since Eq (1.1) exhibits properties akin to those of the standard NLS equation, such as integrability in the Lax and inverse scattering senses, a significant amount of theoretical research has been conducted^[14,15,16], and the study on the equations has extended to higher-order and coupled ones^[17,18,19]. Furthermore, there has been experimental evidence observed in the fields of optics and materials science^[20,21,22].

Subsequently, some other nonlocal integrable equations have been proposed^[23,24,25]. One of these nonlocal equations is the nonlocal DNLS equation^[8,9],

which is obtained from the standard DNLS equation by replacing $|u|^2 u$ with $u(x, t)^2 u^{*}(-x, t)$, and the conserved density $V(x, t) = u(x, t) u^{*}(-x, t)$ satisfies the relation $V^{*}(-x, t) = V(x, t)$. In ^[8], hierarchies of nonlocal DNLS-type equations have been derived utilizing Lie algebra splitting, with the specific nonlocal DNLS equation being termed as nonlocal DNLSI. Zhou^[9] formulates the 2n-fold Darboux transformation (DT) and derives global bounded solutions from zero seed solutions for Eq (1.2). In this paper, we concentrate on how to construct arbitrary order DT of Eq (1.2), what happens if the eigenvalues degenerate, and how to construct solutions from different seeds.

This paper is organized as follows: In Section 2, we give a detailed derivation of the Wronskian-type determinant representation of arbitrary order DT. Moreover, different choice of the eigenvalues may lead different Darboux matrix, even in the same order. If there are duplicate eigenvalues, a generalized form of the DT may prove useful. In Section 3, we obtain solutions of Eq (1.2) from zero and non-zero seed solutions. Periodic solutions and rogue waves are constructed despite of the existence of singularity. Finally, conclusions and discussions are given in Section 4.

2.   Darboux transformation

In this section, we will construct the DT of the coupled nonlocal DNLS equations

which is the compatibility condition of the following Lax pair ^[9]:

where $\Phi = (f(x, t), g(x, t))^{T}$,

with

and Pauli matrix

Under the reduction condition $v(x, t) = -\varepsilon \bar{u}^{*}(x, t) = -\varepsilon u^{*}(-x, t)$, Eq (2.1) will be reduced to the nonlocal DNLS Eq (1.2). Here, $\bar{f}(x, t)$ means $f(-x, t)$. Unlike the local DNLS equation, there are four linear independence solutions of (2.2) associate with one complex eigenvalue. Specifically, if $\left(\begin{array}{cc} f \\ g \\ \end{array} \right)$ are solutions of (2.2) with $\lambda$, then $\left(\begin{array}{cc} f\\ -g \\ \end{array} \right)$, $\left(\begin{array}{cc} \varepsilon \bar{g}^* \\ \bar{f}^* \\ \end{array} \right)$, and $\left(\begin{array}{cc} \varepsilon \bar{g}^* \\ -\bar{f}^* \\ \end{array} \right)$ are solutions of (2.2) with $-\lambda$, $\lambda^*$, and $-\lambda^*$, respectively.

2.1.   The first-order DT

We can suppose $\Phi_{[1]} = T\Phi$ as the DT, where

is the first-order DT matrix, and there exist $U_{[1]}$ and $V_{[1]}$ possessing the same form as $U$ and $V$ with $u$ and $v$ replaced by the new potentials $u_{[1]}$ and $v_{[1]}$, which satisfy Eq (2.2). Therefore, $T[\lambda]$ must satisfy conditions

We set $\Phi_1 = \left(\begin{array}{cc} f_1 \\ g_1 \\ \end{array} \right)$ are solutions of (2.2) with $\lambda = \lambda_1$. Comparing the coefficients of $\lambda^j, (j = 3, 2, 1, 0)$ in Eq (2.8), we obtain $b_1(x, t) = c_1(x, t) = a_0(x, t) = d_0(x, t) = 0$, $b_0(x, t) = b_0, c_0(x, t) = c_0$ and $a_1(x, t) = \frac{g_1}{ f_1}, d_1(x, t) = \frac{f_1}{ g_1}$ by setting $b_0 = c_0 = -\lambda_1$, so we have

where $\rho_1 = \frac{ g_1}{f_1}$.

Theorem 1. Suppose $\lambda_1\neq 0$, which is real when $\varepsilon = 1$ or purely imaginary when $\varepsilon = -1$. Define $\Phi_{[1]}$ and $\rho_1$ as above such that $\rho_1 \bar{\rho_1}^* = 1$. Then Eq (2.9) are new solutions of (1.2), and the corresponding Darboux matrix is

Proof. We only need to prove $v_{[1]} = -\varepsilon \bar{u}_{[1]}{}^*$, which can be easily verified by using (2.9) and $\rho_1 \bar{\rho_1}^*$$\geq$$0.$

□

As $\lambda_1$ will be real or purely imaginary, one-fold DT only allows two eigenvalues $\lambda_1$ and $-\lambda_1$, so one-fold DT is degenerated in a sense. If we choose $\lambda_1$, $-\lambda_1$, $\lambda_1{}^*$, and $-\lambda_1{}^*$ as four different eigenvalues, the real and imaginary part of $\lambda_1$ are both not equal to zero and we need to construct multi-fold DT.

2.2.   The second-order DT

In this subsection, we will construct the two-fold DT of (1.2). The multi-fold DT can be iterated by the one-fold DT, so we set second-order Darboux matrix as

We set $\Phi_k = \left(\begin{array}{cc} f_k \\ g_k \\ \end{array} \right)$ are solutions of (2.2) with $\lambda = \lambda_k$, $\rho_k = \frac{ g_k}{f_k} (k = 1, 2)$, and then $T_2$ maps $(u, v, \Phi)$ to $(u_{[2]}, v_{[2]}, \Phi_{[2]})$. By following the previous steps, we obtain

with

If $\lambda_2 = \lambda_1{}^*$ and $\lambda_1{}^* \neq \pm \lambda_1$, we can obtain Theorem 1 in ^[9]. Moreover, if $\lambda_1, \lambda_2$ are both real or purely imaginary, we can also obtain new solutions of (1.2).

2.3.   The multi-fold DT

From above results, we know the DT will be iterated two times if $\lambda$ is a complex eigenvalue and one times if $\lambda$ is real or purely imaginary. Therefore, if we choose $\lambda_k (k = 0, 1, 2, \cdot\cdot\cdot, n_1)$ as complex eigenvalues and $\lambda_j (j = 0, 1, \cdot\cdot\cdot, n_2)$ as real $(\varepsilon = 1)$ or purely imaginary $(\varepsilon = -1)$ eigenvalues, then we can construct the $(2n_1+n_2)$-fold DT of (1.2). For example, we can construct third-order DT by choosing $(n_1, n_2) = (1, 1)$ or $(0, 3)$ and fourth-order DT by choosing $(n_1, n_2) = (2, 0)$, $(1, 2)$, or $(0, 4)$. Even though the $2n$-fold DT have been derived in ^[9] by choosing $n_1 = n, n_2 = 0$, we construct the Wronskian-type determinant representation of the Darboux matrix.

• $2n$-order DT

We set $2n$-order Darboux matrix as

and $\Phi_k = \left(\begin{array}{cc} f_k \\ g_k \\ \end{array} \right)$ are solutions of (2.2) with $\lambda = \lambda_k$, $\rho_k = \frac{ g_k}{f_k} \, (k = 1, 2, \cdots, 2n)$, and then $T_{2n}$ maps $(u, v, \Phi)$ to $(u_{[2n]}, v_{[2n]}, \Phi_{[2n]})$, where $a_{2k}, d_{2k}\, (k = 0, 1, \cdots, n)$, $b_{2k-1}, c_{2k-1} (k = 1, 2, \cdots, n)$ are $4n+2$ unknown functions. Solving $4n$ equations $T_{2n}(\lambda_k)\Phi_{[2n]} = 0 \, (k = 1, 2, \cdots, 2n)$ by setting $d_{2n} = \frac{1}{a_{2n}}$ and $d_0 = a_0$, we can uniquely determine the $4n+2$ unknown functions, and obtain

with

and

where

and other matrix blocks with the following blocks have been similarly defined.

As $D_1 = \prod_{j = 1}^{2n} \lambda_j D_4$, $D_2 = \prod_{j = 1}^{2n} \lambda_j D_3$, and $a_0 = \prod_{j = 1}^{2n} \lambda_j$, then $d_{2n} = \frac{1}{a_{2n}}$ will be satisfied automatically.

If we choose $\lambda_{n+j} = \lambda_j^*$ and $\Phi_{n+j} = \left(\begin{array}{cc} \varepsilon \bar{g_j}^{*} \\ \bar{f_j}^{*}\\ \end{array} \right)\, (j = 1, 2, \cdots, n)$, then $a_0 = \prod_{j = 1}^{n} |\lambda_j|^2$, $D_1 = \begin{vmatrix} F^{(2, 4, \cdots, 2n)}_{(1, 2, \cdots, n)} & G^{(1, 3, \cdots, 2n-1)}_{(1, 2, \cdots, n)} \\ \varepsilon \bar{G}^{*(2, 4, \cdots, 2n)}_{(1, 2, \cdots, n)} & \bar{F}^{*(1, 3, \cdots, 2n-1)}_{(1, 2, \cdots, n)} \end{vmatrix}\,$ and so are $D_2, \cdots, D_6$, where $\bar{G}^{*(2, 4, \cdots, 2n)}_{(1, 2, \cdots, n)}$ means doing complex conjugate and space-reverse operation to $G^{(2, 4, \cdots, 2n)}_{(1, 2, \cdots, n)}$.

Theorem 2. Suppose $\lambda_j \, (j = 1, 2, \cdots, n)$ are complex, which are all not real or purely imaginary. Then Eq (2.15) are new solutions of (1.2). Especially, when $n = 1$, (2.15) become (2.12).

Proof. We only need to prove $v_{[2n]} = -\varepsilon \bar{u}_{[2n]}{}^*$, which can be easily verified by using $\bar{D_3}^* = D_4,$ and $\bar{D_5}^* = \varepsilon D_6.$ □

Remark 1. It is a little different if we choose $\lambda_{k+j} = \lambda_j^*$, $\Phi_{k+j} = \left(\begin{array}{cc} \varepsilon \bar{g_j}^{*} \\ \bar{f_j}^{*}\\ \end{array} \right)\, (j = 1, 2, \cdots, k)$ and $\lambda_{s}\, (s = 2k+1, \cdots, 2n)$ are all real $(e.g., \varepsilon = 1)$, here $a_0 = \prod_{j = 1}^{k} |\lambda_j|^2 \prod_{s = 2k+1}^{2n} \lambda_s$, $D_1 = \begin{vmatrix} F^{(2, 4, \cdots, 2n)}_{(1, 2, \cdots, k)} & G^{(1, 3, \cdots, 2n-1)}_{(1, 2, \cdots, k)} \\ \varepsilon \bar{G}^{*(2, 4, \cdots, 2n)}_{(1, 2, \cdots, k)} & \bar{F}^{*(1, 3, \cdots, 2n-1)}_{(1, 2, \cdots, k)}\\ F^{(2, 4, \cdots, 2n)}_{(2k+1, \cdots, 2n)} & G^{(1, 3, \cdots, 2n-1)}_{(2k+1, \cdots, 2n)} \end{vmatrix}\,$ and so are $D_2, \cdots, D_6$. Equation (2.15) is new solutions of (1.2) with constraints $\rho_s \bar{\rho_s}^* = 1\, (s = 2k+1, \cdots, 2n)$.

• $(2n-1)$-order DT

We set $(2n-1)$-order Darboux matrix as

and $\Phi_k = \left(\begin{array}{cc} f_k \\ g_k \\ \end{array} \right)$ are solutions of (2.2) with $\lambda = \lambda_k$, $\rho_k = \frac{ g_k}{f_k} \, (k = 1, 2, \cdots, 2n-1)$. Then $T_{2n-1}$ maps $(u, v, \Phi)$ to $(u_{[2n-1]}, v_{[2n-1]}, \Phi_{[2n-1]})$, where $a_{2k-1}, d_{2k-1}\, (k = 1, 2, \cdots, n)$, $b_{2k}, c_{2k} \, (k = 0, 1, \cdots, n-1)$ are $4n$ unknown functions. Solving $4n-2$ equations $T_{2n-1}(\lambda_k)\Phi_{[2n-1]} = 0 \, (k = 1, 2, \cdots, 2n-1)$ by setting $d_{2n-1} = \frac{1}{a_{2n-1}}$ and $c_0 = b_0$, we can uniquely determine the $4n$ unknown functions, and obtain

with

where

Noticing $D_1 = \prod_{j = 1}^{2n-1} \lambda_j D_4$, $D_2 = \prod_{j = 1}^{2n-1} \lambda_j D_3$ and $b_0 = -\prod_{j = 1}^{2n-1} \lambda_j$, then $d_{2n} = \frac{1}{a_{2n}}$ will be satisfied automatically.

If we choose $\lambda_{k+j} = \lambda_j^*$ and $\Phi_{k+j} = \left(\begin{array}{cc} \varepsilon \bar{g_j}^{*} \\ \bar{f_j}^{*}\\ \end{array} \right)\, (j = 1, 2, \cdots, k)$, then $b_0 = - \prod_{j = 1}^{k} |\lambda_j|^2 \prod_{s = 2k+1}^{2n-1} \lambda_{s}$, $D_1 = \begin{vmatrix} F^{(1, 3, \cdots, 2n-1)}_{(1, 2, \cdots, k)} & G^{(2, 4, \cdots, 2n-2)}_{(1, 2, \cdots, k)} \\ \varepsilon \bar{G}^{*(1, 3, \cdots, 2n-1)}_{(1, 2, \cdots, k)} & \bar{F}^{*(2, 4, \cdots, 2n-2)}_{(1, 2, \cdots, k)}\\ F^{(1, 3, \cdots, 2n-1)}_{(2k+1, 2k+2, \cdots, 2n-1)} & G^{(2, 4, \cdots, 2n-2)}_{(2k+1, 2k+2, \cdots, 2n-1)} \end{vmatrix}\,$ and so are $D_2, \cdots, D_6$.

Theorem 3. Suppose $\lambda_j \, (j = 1, 2, \cdots, k)$ are complex, and $\lambda_s \, (s = 2k+1, 2k+2, \cdots, 2n-1)$ are real when $\varepsilon = 1$ or purely imaginary when $\varepsilon = -1$, then Eq (2.21) is new solutions of (1.2) with constraints $\rho_s \bar{\rho_s}^* = 1\, (s = 2k+1, \cdots, 2n-1)$. Especially, when $n = 1$, (2.21) become (2.9).

Proof. Similarly, we only need to prove $v_{[2n-1]} = -\varepsilon \bar{u}_{[2n-1]}{}^*$, which can be easily verified by using $\bar{D_3}^* = (-1)^{(n^2-1)} \varepsilon^{n} D_4,$ and $\bar{D_5}^* = (-1)^{(n^2-1)} \varepsilon^{n-1} D_6.$ □

Until now, we have constructed the arbitrary DT of (1.2) with Wronskian-type determinant representation. Compared with ^[9], we have not only constructed the DT of arbitrary order, but also adopted a simpler and more feasible way.

3.   Solutions of Eq (1.2)

In this section, we will construct several explicit solutions of Eq (1.2).

3.1.   Solutions by applying the first-order DT

• If we set seed solution $u = v = 0$, then solution of (2.2) with $\lambda = \lambda_k$ is $\Phi_k = \left(\begin{array}{cc} f_k \\ g_k \\ \end{array} \right) = \left(\begin{array}{cc} \alpha_k e^{\theta_k} \\ \beta_k e^{-\theta_k} \\ \end{array} \right)$, where $\theta_k = \lambda_k^2 x-2 i \lambda_k^4 t \, (k = 1, 2, \cdots)$. As $\lambda_1$ is real when $\varepsilon = 1$ or purely imaginary when $\varepsilon = -1$ in Theorem 1, background of the new solution of (1.2) will go to infinity, which appears meaningless.

• If we set seed solution $u = \rho e^{-\varepsilon \rho ^2 x +i \phi }$, the background will also go to infinity. So, we take $u = \rho e^{i \phi }$ as the seed solution eclectically to avoid infinity of the background. Solution of (2.2) with $\lambda = \lambda_k$ is $\Phi_k = \left(\begin{array}{cc} f_k \\ g_k \end{array} \right)$, where

with $\theta_k = \sqrt{\lambda_k^4-\lambda_k^2 \rho ^2 \varepsilon } \left[x-i t \left(2 \lambda_k^2+\rho ^2 \varepsilon \right)\right]$, $\rho_k = \frac{ g_k}{f_k}$, $\alpha_k$ and $\beta_k \, (k = 1, 2, \cdots)$ are none-zero complex constants.

Using Theorem 1, we obtain

under the constraint

which means $|\lambda _1| > |\rho|$. This solution has been shown in Figure 1, which displays the periodic one-soliton with $\varepsilon = 1, \rho = 4, \lambda_1 = 5, \alpha _1 = 2, \beta _1 = 1$, $\phi = \pi$, and $\varepsilon = -1, \rho = 4, \lambda_1 = -5i, \alpha _1 = 1/2, \beta _1 = 1$, $\phi = 0$. However, there are periodic singularities in any case for $\varepsilon$. The reason for the existence of the singularity may be that $\theta_k$ in solution (3.1) is too complex. When $\varepsilon$ takes the values of -1 and 1, respectively, eigenvalue $\lambda$ will differ due to the constraints, which subsequently causes variations in other parameters. Additionally, we found that the positions of the singular points have changed.

Figure 1.  The periodic one-soliton from (3.2) with (a) $\varepsilon = 1, \rho = 4, \lambda_1 = 5, \alpha _1 = 2, \beta _1 = 1$, $\phi = \pi$, and (b) $\varepsilon = -1, \rho = 4, \lambda_1 = -5i, \alpha _1 = 1/2, \beta _1 = 1$, $\phi = 0$.

DownLoad: Full-Size Img PowerPoint

When the eigenvalue $\lambda_1 = \pm \rho \, (\varepsilon = 1)$ or $\lambda_1 = \pm i \rho \, (\varepsilon = -1)$, the new solution will be complex constant and trivial, so we should resolve (2.2). Now, the solution of (2.2) $(\varepsilon = 1)$ is $\Phi_1 = \left(\begin{array}{cc} f_1 \\ g_1 \end{array} \right)$, where

with $\alpha _k \, (k = 1, 2)$ are none-zero complex constants. We can obtain a new fundamental rational solution,

The graph of this solution is displayed in Figure 2(a) with $\rho = 1, \alpha _1 = 1, \alpha _2 = 2$ and $\phi = \pi$. It is seen that this is a rogue wave, similar to Eq (40) in Yang's article^[26], which blows up to infinity at $x = -1/2$ and finite time $t = 0$. Furthermore, rising from a constant background, rogue wave (3.5) can be considered as resulting from the degradation of periodic soliton (3.2). If we set $\varepsilon = -1$, then new rogue wave will blow up to infinity at $x = 1/2$ and $t = 0$, which is shown in Figure 2(b). Similarly, variations in the value of parameter $\varepsilon$ only affects the position of solitons or rogue waves, without altering the type or mechanical properties of the solutions. Meanwhile, parameter $\phi$ do not affect the modulus of the solutions.

Figure 2.  The one rogue wave from (3.5) with (a) $\varepsilon = 1, \rho = 1, \lambda_1 = -1, \alpha _1 = 1, \alpha _2 = 2$, $\phi = \pi$, and (b) $\varepsilon = -1, \rho = 1, \lambda_1 = -i, \alpha _1 = -2, \alpha _2 = 4$, $\phi = 0$.

DownLoad: Full-Size Img PowerPoint

3.2.   Solutions by applying the second-order DT

• When we set seed solution $u = v = 0$, the background of new solutions will not go to infinity if and only if $\lambda^2$ is pure imaginary, which means $4 arg\lambda/ \pi$ is an odd integer ^[9]. Hence, we can obtain a new soliton solution with $\lambda_1 = a(1-i)$,

This solution does not have singularity and global if $|\alpha _1| \neq |\beta_1|$, which is periodic in both $x$ and $t$. The graph is displayed in Figure 3 with $\varepsilon = -1, a = 1, \alpha _1 = 1$ and $\beta _1 = 2$, from which we can see that the norm of the solution is periodic in $x$ but independent of $t$. However, if $\lambda_1$ and $\lambda_2$ are both real or purely imaginary, the background of new solution will go to infinity.

Figure 3.  Left panel is the periodic two solitons from (3.6) with $\varepsilon = -1, \lambda_1 = 1-i, \alpha _1 = 1$ and $\beta _1 = 2$; Right panel is the corresponding density plot.

DownLoad: Full-Size Img PowerPoint

• If we set seed solution $u = \rho e^{-\varepsilon \rho ^2 x +i \phi }$, the background will go to infinity. So, we take $u = \rho e^{i \phi }$ as the seed solution, and the solution of (2.2) with $\lambda = \lambda_k$ is (3.1). Using Eq (2.12), we can obtain new solutions of Eq (1.2) and display it in Figure 4. Figure 4 show the propagation of two solitons, the interaction of which seems elastic. However, the solution blows up to infinity in the interaction region at one fixed point. The reason for the existence of singularity may be due to nonlocal effects. Figure 5(a) shows the interaction of two periodic soliton with periodic singularities, as long as $\lambda_1, \lambda_2$ are both real numbers. This solution seems to be a superposition and interaction of two solutions (3.2). Furthermore, if $\lambda_1 = \pm \rho$, then we can see that the interaction of periodic soliton and rogue wave, which has been shown in Figure 5(b), will degenerate into the seed solution if $\lambda_2 = \pm \rho$.

Figure 4.  Left panel is the interaction of two solitons from (2.12) with $\varepsilon = 1, \rho = 1, \lambda_1 = (1+i\sqrt{3})/2, \alpha _1 = 1, \beta_1 = i, \phi = 0$; Right panel is the corresponding density plot.

DownLoad: Full-Size Img PowerPoint

Figure 5.  (a) Inteaction of two periodic solitons from (2.12) with $\varepsilon = 1, \rho = 1, \lambda_1 = 2$, $\lambda_2 = 3, \alpha _1 = \sqrt{\frac{2\sqrt{3}+3}{2\sqrt{3}-3}}, \alpha _2 = \sqrt{\frac{6\sqrt{2}+8}{6\sqrt{2}-8}}, \; \beta_1 = \beta_2 = 1, \; \phi = 0$; (b) Interaction of one periodic soliton and rogue wave from (2.12) with $\varepsilon = 1, \rho = 1, \lambda_1 = -1, \lambda_2 = 2, \alpha _1 = 1+i, \alpha _2 = \sqrt{\frac{2\sqrt{3}+3}{2\sqrt{3}-3}}, \beta_1 = 2, \beta_2 = 1, \phi = 0$.

DownLoad: Full-Size Img PowerPoint

3.3.   Solutions by applying the multi-fold DT

• Considering the third-order DT, three eigenvalues can be one complex and one real (or pure imaginary), or three real. New solutions derived from seeds $u = \rho e^{i \phi }$ would be the superposition of previous solutions (2.9) and (2.12). For example, if we choose three eigenvalues as one complex and one real, new solutions is displayed in Figure 6(a). We can see Figure 6(a) is a superposition and interaction of Figure 1(a) and Figure 4. There is one periodic soliton and two elastic-interacted solitons. If the real eigenvalue equals $\pm \rho$, then the solution will be different from Figure 6(a), which is shown in Figure 6(b). The phases of the solitons undergo changes, and two of them remain parallel both before and after the collision.

Figure 6.  Interaction of three solitons from (2.21) with (a) $\varepsilon = 1, \rho = 1, \lambda_1 = (1+i \sqrt{3})/2, \lambda_2 = 2, \alpha _1 = -1, \alpha _2 = \sqrt{\frac{2\sqrt{3}+3}{2\sqrt{3}-3}}, \beta_1 = \beta_2 = 1, \phi = 0$; (b) $\varepsilon = 1, \rho = \lambda_1 = 1, \lambda_2 = (1+i \sqrt{3})/2, \alpha _1 = \alpha _2 = \beta_1 = \beta_2 = 1, \phi = 0$; (c) $\varepsilon = 1, \rho = 1, \lambda_1 = 2, \lambda_2 = 3$, $\lambda_3 = 4$, $\alpha _1 = \sqrt{\frac{2\sqrt{3}+3}{2\sqrt{3}-3}}, \alpha _2 = \sqrt{\frac{6\sqrt{2}+8}{6\sqrt{2}-8}}, \alpha _3 = \sqrt{\frac{4\sqrt{15}+15}{4\sqrt{15}-15}}, \beta_1 = \beta_2 = \beta_3 = 1, \phi = 0$; (d) $\varepsilon = 1, \rho = 1$, $\lambda_1 = -1, \lambda_2 = 2, \lambda_3 = 2, \alpha _1 = 1, \alpha _2 = \sqrt{\frac{2\sqrt{3}+3}{2\sqrt{3}-3}}, \alpha _3 = \sqrt{\frac{6\sqrt{2}+8}{6\sqrt{2}-8}}, \beta_1 = 2, \beta_2 = \beta_3 = 1, \phi = 0$.

DownLoad: Full-Size Img PowerPoint

If three eigenvalues are all selected real, then the new solution is displayed in Figure 6(c), which are three times superposition and interaction of Figure 1(a). Further more, if one of three real eigenvalues equals $\pm\rho$, then there will be a rogue wave, which has been shown in Figure 6(d).

• Considering the fourth-order DT, we choose zero seed solutions and two complex eigenvalues $\lambda_1$, $\lambda_2$, where $\lambda_k^2 (k = 1, 2)$ is pure imaginary to avoid infinity of the background. If $\lambda_2 \neq \lambda_1$, we can obtain periodic solitons, which have been shown in Figure 7. If $\lambda_2 = \lambda_1$, then solution of (2.2) with $\lambda = \lambda_2$, says $\Phi_2$, will linearly dependent on $\Phi_1$. Following the idea of generalized DT in ^[7], we set $\Phi_1^{[1]} = \Phi_1(\lambda_1+\delta), \delta\ll1$, and obtain

which is another solution of (2.2) with $\lambda = \lambda_1$ but linearly independent on $\Phi_1$. Using (2.15), we obtain new solutions of (1.2),

with $\varepsilon = 1, \alpha _1 = \beta _1 = 1$. This solution is a mixed-type of exponential and rational functions, whose norm is displayed in Figure 8. From it we can see there are three solitons, one $x$-periodic soliton and two elastic-interacted solitons, and no rogue waves are presented obviously. However, the interaction blows up to infinity at two fixed points. As it has been seen that higher order solution can be the superposition and interaction of the lower order solution, we do not show the solutions derived from four real eigenvalues or one complex and two real eigenvalues. At the end, we display new periodic solitons in Figure 9 obtained by using six-order DT from zero seeds.

Figure 7.  Left panel is the periodic four solitons from (2.15) with $\varepsilon = 1, \lambda_1 = 1+i$, $\lambda_2$$=$$2 (1-i)$, $\alpha _1 = \alpha_2 = 1/5$ and $\beta _1 = \beta_2 = 1$; Right panel is the corresponding density plot.

DownLoad: Full-Size Img PowerPoint

Figure 8.  Left panel is interaction of three solitons from (3.8) with $\varepsilon = 1, \lambda_1 = \lambda_2 = 1+i, \alpha _1 = \beta _1 = 1$; Right panel is the corresponding density plot.

DownLoad: Full-Size Img PowerPoint

Figure 9.  Left panel is interaction of periodic multi-solitons from (2.15) (n = 3) with $\varepsilon = 1$, $\lambda_1 = 1+i, \lambda_2 = 2 (1-i), \lambda_3 = 3 (1+i), \alpha _1 = 0.1, \alpha _2 = 0.2, \alpha _3 = 0.3, \beta_1 = \beta_2 = \beta_3 = 1$; Right panel is the corresponding density plot.

DownLoad: Full-Size Img PowerPoint

4.   Conclusions

In this paper, the nonlocal DNLS Eq (1.2) has been analytically studied. Firstly, using the Lax pair of the coupled Eq (2.1), from which the reduction conditions lead to Eq (1.2), we have constructed arbitrary order DTs in a more streamlined manner. From zero and non-zero seed solutions, new Wronskian-type determinant solutions are obtained. Specifically, we have obtained new solutions by applying the first four orders DT.

Starting from zero seeds, due to the nonlocality of the equation, odd-order DTs will cause the background of the solutions to diverge towards infinity. For even-order DTs, the background of the solutions will not go to infinity unless the square of eigenvalues $\lambda^2$ are pure imaginary, in order to maintain the reduction conditions, which could increase the difficulty in selecting parameter values. Soliton solutions without singularities have been obtained for specific parameter values, as demonstrated in Figures 3, 7, and 9. These new solutions are all periodic and bounded, similar to the results reported in ^[9]. Besides, when one eigenvalue tends to another one, generalized DT can be constructed. Taking the solitons in Figure 7 as an example, we can obtain new solutions which are mixed-type of exponential and rational functions with $\lambda_2 = \lambda_1$, which can be seen Figure 8. Although the solution in Figure 7 is bounded with $\lambda_2 \neq \lambda_1$, there are two singularities, causing the solitons to blow up to infinity.

Starting from non-zero seeds, we have obtained various solutions including the periodic solutions and rogue waves. If the eigenvalues are real or pure imaginary numbers, there will be constraints on the parameters. However, these solutions inherently possess singularities. Nevertheless, we have constructed soliton solutions for the first few orders. Figure 1 displays the periodic one-soliton, which degrades into a rogue wave (see Figure 2) when the magnitude of the eigenvalue is identical to the amplitude of the seed solution. Similarly, periodic two- and three-soliton solutions can also degrade or partially degrade into rogue waves, which can be observed in Figures 5, 6(c), and 6(d). When the eigenvalues are taken as complex numbers, the interaction of two solitons has been shown in Figure 4. Figure 6(a) shows a superposition and interaction of Figures 1(a) and 4. When the real eigenvalue equals to $\pm \rho$, the periodic soliton will disappear, as shown in Figure 6(b).

In addition, the value of $\varepsilon$ does not affect much except the location of the singularities, which is different from the nonlocal NLS Eq (1.1). We hope that the work presented in this paper can provide a valuable reference and assistance for the study of nonlocal integrable systems.

Use of Generative-AI tools declaration

The author declares he has not used Artificial Intelligence (AI) tools in the creation of this article.

Acknowledgments

Thanks to the valuable suggestions of the editor and reviewers. Work of this paper was supported by State Key Laboratory of Heavy Oil Processing (SKLHOP2024115808).

Conflict of interest

The author declares no conflict of interest.