# Spectral decomposition for the Dirac system associated to the DSII equation

## 1 Introduction

Gravity-capillary surface wave packets are described by the Davey–Stewartson (DS) system  which is integrable by inverse scattering tranform in the limit of shallow water . In this paper, we study the focusing DSII equation which can be written in a complex form,

$`iu_t+u_{zz}+u_{\overline{z}\overline{z}}+4(g+\overline{g})u`$ $`=`$ $`0,`$   
$`2g_{\overline{z}}\left(|u|^2\right)_z`$ $`=`$ $`0,`$ (1.1)   
where $`z=x+iy`$, $`\overline{z}=xiy`$, $`u(z,\overline{z},t)`$ and $`g(z,\overline{z},t)`$ are complex functions. This equation appears as the compatibility condition for the two-dimensional Dirac system,

$$\phi _{1\overline{z}}=u\phi _2,\phi _{2z}=\overline{u}\phi _1,$$
(1.2)   
coupled to the equations for the time evolution of the eigenfunctions,

$`i\phi _{1t}+\phi _{1zz}+u\phi _{2\overline{z}}u_{\overline{z}}\phi _2+4g\phi _1`$ $`=`$ $`0,`$   
$`i\phi _{2t}+\phi _{2\overline{z}\overline{z}}+\overline{u}_z\phi _1\overline{u}\phi _{1z}+4\overline{g}\phi _2`$ $`=`$ $`0.`$ (1.3)   
The DSII equation was solved formally through the $`\overline{}`$ problem of complex analysis by Fokas and Ablowitz  and Beals and Coifman . Rigorous results on existence and uniqueness of solutions of the initial-value problem were established under a small-norm assumption . The small-norm assumption was used to eliminate homogeneous solutions of equations of the inverse scattering which correspond to bound states and radially symmetric localized waves (lumps) of the DSII equation. When the potential in the linear system becomes weakly localized (in $`L^2`$ but not in $`L^1`$), homogeneous solutions may exist and the analysis developed in Ref.  is not applicable.

The lump solutions were included formally in , where their weak decay rate was found, $`u\mathrm{O}(R^1)`$ as $`R=\sqrt{x^2+y^2}\mathrm{}`$. This result is only valid for complexified solutions of the DSII equation (when $`|u|^2`$ is considered to be complex). The reality conditions were incorporated in the work of Arkadiev et al.  where lumps were shown to decay like $`u\mathrm{O}(R^2)`$. Multi-lump solutions were expressed as a ratio of two determinants , or, in a special case, as a ratio of two polynomials  but their dynamical role was left out of consideration.

Recently, structural instability of a single lump of the DSII equation was reported by Gadyl’shin and Kiselev . The authors used methods of perturbation theory based on completeness of squared eigenfunctions of the Dirac system . A similar conclusion was announced by Yurov who studied Darboux transformation of the Dirac system .

In this paper, we present an alternative solution of the problem of stability of multi-lump solutions of the DSII equation. The approach generalizes our recent work on spectral decomposition of a linear time-dependent Schrödinger equation with weakly localized (not in $`L^1`$) potentials . We find a new spectral decomposition in terms of single eigenfunctions of the Dirac system. Surprisely enough, the two-component Dirac system in two dimensions has a scalar spectral decomposition. In contrast, we recall that the Dirac system in one dimension (the so-called AKNS system) has a well-known 2 x 2 matrix spectral decomposition .

Using the scalar spectral decomposition, we associate the multi-lump potentials with eigenvalues embedded into a two-dimensional essential spectrum of the Dirac system. Eigenvalues embedded into a one-dimensional essential spectrum occur for instance for the time-dependent Schrödinger problem . They were found to be structurally unstable under a small variation of the potential. Depending on the sign of the variation, they either disappear or become resonant poles in the complex spectral plane which correspond to lump solutions of the KPI equation .

For the Dirac system in two dimensions, the multi-lump potentials and embedded eigenvalues are more exotic. The discrete spectrum of the Dirac system is separated from the continuous spectrum contribution in the sense that the spectral data satisfy certain constraints near the embedded eigenvalues. These constraints are met for special solutions of the DSII equation such as lumps, but may not be satisfied for a generic combination of lumps and radiative waves. As a result, embedded eigenvalues of the Dirac system generally disappear under a local disturbance of the initial data. Physically, this implies that a localized initial data of the DSII equation decays into radiation except for the cases where the data reduce to special solutions such as lumps.

The paper is organized as follows. Elements of inverse scattering for the Dirac system are reviewed in Section 2, where we find that the discrete spectrum of the Dirac system is prescribed by certain constraints on the spectral data. Spectral decomposition is described in Section 3 with the proof of orthogonality and completeness relations through a proper adjoint problem. The perturbation theory for lumps is developed in Section 4 where some of previous results  are recovered. Section 5 contains concluding remarks. Appendix A provides a summary of formulas of the complex $`\overline{}`$-analysis used in proofs of Section 3.

## 2 Spectral Data and Inverse Scattering

Here we review some results on the Dirac system (1.2) and discard henceforth the time dependence of $`u`$, $`g`$ and $`𝝋`$. The potential $`u(z,\overline{z})`$ is assumed to be non-integrable ($`uL^1`$) with the boundary conditions, $`u\mathrm{O}(|z|^2)`$ as $`|z|\mathrm{}`$.

### 2.1 Essential Spectrum of the Dirac System

We define the fundamental matrix solution of Eq. (1.2) in the form ,

$$𝝋=[𝝁(z,\overline{z},k,\overline{k})e^{ikz},𝝌(z,\overline{z},k,\overline{k})e^{i\overline{k}\overline{z}}],$$
(2.1)   
where $`k`$ is a spectral parameter, $`𝝁(z,\overline{z},k,\overline{k})`$ and $`𝝌(z,\overline{z},k,\overline{k})`$ satisfy the system,

$`\mu _{1\overline{z}}`$ $`=`$ $`u\mu _2,\mu _{2z}=ik\mu _2+\overline{u}\mu _1,`$ (2.2)   
$`\chi _{1\overline{z}}`$ $`=`$ $`i\overline{k}\chi _1u\chi _2,\chi _{2z}=\overline{u}\chi _1.`$ (2.3)   
It follows from Eqs. (2.2) and (2.3) that $`𝝁`$ and $`𝝌`$ are related by the symmetry constraint,

$$𝝌(z,\overline{z},k,\overline{k})=𝝈\overline{𝝁}(z,\overline{z},k,\overline{k}),𝝈=\left(\begin{array}{cc}0& 1\\ 1& 0\end{array}\right).$$
(2.4)   
We impose the boundary conditions for $`𝝁(z,\overline{z},k,\overline{k})`$,

$$\underset{|k|\mathrm{}}{lim}𝝁(z,\overline{z},k,\overline{k})=𝐞_1=\left(\begin{array}{cc}1& \\ 0& \end{array}\right).$$
(2.5)   
Solutions of Eq. (2.2) with boundary conditions (2.5) can be expressed through the Green’s functions as Fredholm’s inhomogeneous integral equations ,

$`\mu _1(z,\overline{z},k,\overline{k})`$ $`=`$ $`1{\displaystyle \frac{1}{2\pi i}}{\displaystyle \frac{dz^{}d\overline{z}^{}}{z^{}z}(u\mu _2)(z^{},\overline{z}^{})},`$ (2.6)   
$`\mu _2(z,\overline{z},k,\overline{k})`$ $`=`$ $`{\displaystyle \frac{1}{2\pi i}}{\displaystyle \frac{dz^{}d\overline{z}^{}}{\overline{z}^{}\overline{z}}(\overline{u}\mu _1)(z^{},\overline{z}^{})e^{ik(zz^{})i\overline{k}(\overline{z}\overline{z}^{})}}.`$ (2.7)   
Values of $`k`$ for which the homogeneous system associated to Eqs. (2.6) and (2.7) has bounded solutions are called eigenvalues of the discrete spectrum of the Dirac system. Let us suppose that the homogeneous solutions (eigenvalues) are not supported by the potential $`u(z,\overline{z})`$. We evaluate the departure from analyticity of $`𝝁`$ in the $`k`$ plane by calculating the derivative $`𝝁/\overline{k}`$ directly from the system (2.6)-(2.7)  as

$$\frac{𝝁}{\overline{k}}=b(k,\overline{k})𝐍_\mu (z,\overline{z},k,\overline{k}).$$
(2.8)   
Here $`b(k,\overline{k})`$ is the spectral data,

$$b(k,\overline{k})=\frac{1}{2\pi }𝑑zd\overline{z}\left(\overline{u}\mu _1\right)(z,\overline{z})e^{i(kz+\overline{k}\overline{z})},$$
(2.9)   
and $`𝐍_\mu (z,\overline{z},k,\overline{k})`$ is a solution of Eq. (2.2) which is linearly independent of $`𝝁(z,\overline{z},k,\overline{k})`$ and satisfies the boundary condition,

$$\underset{|k|\mathrm{}}{lim}𝐍_\mu (z,\overline{z},k,\overline{k})e^{i(kz+\overline{k}\overline{z})}=𝐞_2=\left(\begin{array}{cc}0& \\ 1& \end{array}\right).$$
(2.10)   
This solution can be expressed through the Fredholm’s inhomogeneous equations,

$`N_{1\mu }(z,\overline{z},k,\overline{k})`$ $`=`$ $`{\displaystyle \frac{1}{2\pi i}}{\displaystyle \frac{dz^{}d\overline{z}^{}}{z^{}z}(uN_{2\mu })(z^{},\overline{z}^{})},`$ (2.11)   
$`N_{2\mu }(z,\overline{z},k,\overline{k})`$ $`=`$ $`e^{i(kz+\overline{k}\overline{z})}+{\displaystyle \frac{1}{2\pi i}}{\displaystyle \frac{dz^{}d\overline{z}^{}}{\overline{z}^{}\overline{z}}(\overline{u}N_{1\mu })(z^{},\overline{z}^{})e^{ik(zz^{})i\overline{k}(\overline{z}\overline{z}^{})}}.`$ (2.12)   
The following reduction formula  connects $`𝐍_\mu (z,\overline{z},k,\overline{k})`$ and $`𝝁`$$`(z,\overline{z},k,\overline{k})`$,

$$𝐍_\mu (z,\overline{z},k,\overline{k})=𝝌(z,\overline{z},k,\overline{k})e^{i(kz+\overline{k}\overline{z})}=𝝈\overline{𝝁}(z,\overline{z},k,\overline{k})e^{i(kz+\overline{k}\overline{z})}.$$
(2.13)   
If the potential $`u(z,\overline{z})`$ has the boundary values $`u\mathrm{O}(|z|^2)`$ as $`|z|\mathrm{}`$, ($`uL^1`$), then, the integral kernel in Eq. (2.9) is not absolutely integrable, while Eqs. (2.6) and (2.7) are still well-defined. We specify complex integration in the $`z`$-plane of a non-absolutely integrable function $`f(z,\overline{z})`$ according to the formula,

$$𝑑zd\overline{z}f(z,\overline{z})=\underset{R\mathrm{}}{lim}_{|z|R}𝑑zd\overline{z}f(z,\overline{z}).$$
(2.14)   
The same formula is valid for integrating eigenfunctions of the Dirac system in the $`k`$ plane as well. In Section 3, we use (2.14) when computing the inner products and completeness relations for the Dirac system and its adjoint.

### 2.2 The Discrete Spectrum

Suppose here that integral equations (2.6) and (2.7) have homogeneous solutions at an eigenvalue $`k=k_j`$. The discrete spectrum associated to multi-lump potentials was introduced in Refs. . Here we review their approach and give a new result (Proposition 2.1) which clarifies the role of discrete spectrum in the spectral problem (2.2).

For the discrete spectrum associated to the multi-lump potentials, an isolated eigenvalue $`k=k_j`$ has double multiplicity with the corresponding two bound states $`𝚽_j(z,\overline{z})`$ and $`𝚽_j^{}(z,\overline{z})`$ . The bound state $`𝚽_j(z,\overline{z})`$ is a solution of the homogeneous equations,

$`\mathrm{\Phi }_{1j}(z,\overline{z})`$ $`=`$ $`{\displaystyle \frac{1}{2\pi i}}{\displaystyle \frac{dz^{}d\overline{z}^{}}{z^{}z}(u\mathrm{\Phi }_{2j})(z^{},\overline{z}^{})},`$ (2.15)   
$`\mathrm{\Phi }_{2j}(z,\overline{z})`$ $`=`$ $`{\displaystyle \frac{1}{2\pi i}}{\displaystyle \frac{dz^{}d\overline{z}^{}}{\overline{z}^{}\overline{z}}(\overline{u}\mathrm{\Phi }_{1j})(z^{},\overline{z}^{})e^{ik_j(zz^{})i\overline{k}_j(\overline{z}\overline{z}^{})}},`$ (2.16)   
with the boundary conditions as $`|z|\mathrm{}`$,

$$𝚽_j(z,\overline{z})\frac{𝐞_1}{z}.$$
(2.17)   
Equivalently, this boundary condition can be written as renormalization conditions for Eqs. (2.15) and (2.16),

$`{\displaystyle \frac{1}{2\pi i}}{\displaystyle 𝑑z}d\overline{z}(u\mathrm{\Phi }_{2j})(z,\overline{z})`$ $`=`$ $`1,`$ (2.18)   
$`{\displaystyle \frac{1}{2\pi i}}{\displaystyle 𝑑z}d\overline{z}(\overline{u}\mathrm{\Phi }_{1j})(z,\overline{z})e^{i(k_jz+\overline{k}_j\overline{z})}`$ $`=`$ $`0.`$ (2.19)   
The other (degenerate) bound state $`𝚽_j^{}(z,\overline{z})`$ can be expressed in terms of $`𝚽_j(z,\overline{z})`$ using Eq. (2.13),

$$𝚽_j^{}(z,\overline{z})=𝝈\overline{𝚽}_j(z,\overline{z})e^{i(k_jz+\overline{k}_j\overline{z})}.$$
(2.20)   
The behaviour of the eigenfunction $`𝝁(z,\overline{z},k,\overline{k})`$ near the eigenvalue $`k=k_j`$ becomes complicated due to the fact that the double eigenvalue is embedded into the two-dimensional essential spectrum of the Dirac system (2.2). We prove the following result.   

Proposition 2.1. For smooth data $`b(k,\overline{k})C^1`$ at $`kk_j`$, the eigenfunction $`𝝁(z,\overline{z},k,\overline{k})`$ has a pole singularity at $`kk_j`$ only if

$$b_0=\frac{1}{2\pi }𝑑zd\overline{z}\left(\overline{z}\overline{u}\mathrm{\Phi }_{1j}\right)(z,\overline{z})e^{i(k_jz+\overline{k}_j\overline{z})}=0.$$
(2.21)   
Proof. Suppose $`𝝁(z,\overline{z},k,\overline{k})`$ has a pole singularity at $`k=k_j`$. Then, it can be shown from Eq. (2.2) that the meromorphic continuation of $`𝝁(z,\overline{z},k,\overline{k})`$ is given by the limiting relation,

$$\underset{kk_j}{lim}\left[𝝁(z,\overline{z},k,\overline{k})\frac{i𝚽_j(z,\overline{z})}{kk_j}\right]=(z+z_j)𝚽_j(z,\overline{z})+c_j𝚽_j^{}(z,\overline{z}),$$
(2.22)   
where $`z_j`$, $`c_j`$ are some constants. Using Eqs. (2.8), (2.9), and (2.13), we find the differential relation for $`b(k,\overline{k})`$,

$$\frac{b}{\overline{k}}=\frac{b(k,\overline{k})}{2\pi }𝑑zd\overline{z}\left(\overline{u}\overline{\mu }_2\right)(z,\overline{z})\frac{1}{2\pi i}𝑑zd\overline{z}\left(\overline{z}\overline{u}\mu _1\right)(z,\overline{z})e^{i(kz+\overline{k}\overline{z})}.$$

In the limit $`kk_j`$, this equation reduces with the help of Eqs. (2.18) and (2.22) to the form,

$$\frac{b}{\overline{k}}=\frac{b(k,\overline{k})}{\overline{k}\overline{k}_j}\frac{b_0}{kk_j},$$

where $`b_0`$ is given in Eq. (2.21). The reduced equation exhibits the limiting behavior of $`b(k,\overline{k})`$ as $`kk_j`$,

$$b(k,\overline{k})b_0\frac{\overline{k}\overline{k}_j}{kk_j}\mathrm{ln}|\overline{k}\overline{k}_j|.$$
(2.23)   
On the other hand, it follows from Eqs. (2.13), (2.20), and (2.22) that $`𝐍_\mu (z,\overline{z},k,\overline{k})`$ has the limiting behavior,

$$𝐍_\mu (z,\overline{z},k,\overline{k})\frac{i𝚽_j^{}(z,\overline{z})}{\overline{k}\overline{k}_j}.$$
(2.24)   
According to Eqs. (2.23) and (2.24), the right-hand-side of Eq. (2.8) is of order $`\mathrm{O}(b_0|kk_j|^1\mathrm{ln}|kk_j|)`$ as $`kk_j`$. On the other hand, the left-hand-side of Eq. (2.8) must be of order $`O(1)`$ in the limit $`kk_j`$ according to Eq. (2.22). Therefore, the eigenfunction $`𝝁(z,\overline{z},k,\overline{k})`$ has a pole at $`k=k_j`$ only if the constraint $`b_0=0`$ holds. $`\mathrm{}`$   

The limiting relation (2.22) was introduced by Arkadiev et al . However, the authors did not notice that the discrete spectrum is supported only by potentials which satisfy the additional constraint (2.21). In particular, such potentials include the multi-lump solutions for which $`b(k,\overline{k})=0`$ everywhere in the $`k`$-plane.

### 2.3 Expansion Formulas for Inverse Scattering

Combining Eq. (2.8) for the essential spectrum and Eq. (2.22) for the discrete spectrum, we reconstruct the eigenfunction $`𝝁(z,\overline{z},k,\overline{k})`$ ,

$$𝝁(z,\overline{z},k,\overline{k})=𝐞_1+\underset{j=1}{\overset{n}{}}\frac{i𝚽_j(z,\overline{z})}{kk_j}+\frac{1}{2\pi i}\frac{dk^{}d\overline{k}^{}}{k^{}k}b(k^{},\overline{k}^{})𝐍_\mu (z,\overline{z},k^{},\overline{k}^{}),$$
(2.25)   
where $`n`$ is number of distinct eigenvalues $`k_j`$ of double multiplicity. At $`kk_j`$, this system is coupled with the algebraic system for the bound states,

$$(z+z_j)𝚽_j(z,\overline{z})+c_j𝚽_j^{}(z,\overline{z})=𝐞_1+\underset{lj}{}\frac{i𝚽_l(z,\overline{z})}{k_jk_l}+\frac{1}{2\pi i}\frac{dkd\overline{k}}{kk_j}b(k,\overline{k})𝐍_\mu (z,\overline{z},k,\overline{k}).$$
(2.26)   
Expansion (2.25) can be related to the inverse scattering transform for the potential $`u(z,\overline{z})`$ . It follows from Eq. (2.2) that the eigenfunction $`𝝁(z,\overline{z},k,\overline{k})`$ has the asymptotic expansion as $`|k|\mathrm{}`$,

$$𝝁(z,\overline{z},k,\overline{k})=𝐞_1+\frac{1}{ik}𝝁_{\mathrm{}}(z,\overline{z})+\mathrm{O}(|k|^2),$$
(2.27)   
where $`\mu _2\mathrm{}(z,\overline{z})=\overline{u}(z,\overline{z})`$ and

$$\mu _1\mathrm{}(z,\overline{z})=\frac{1}{2\pi i}\frac{dz^{}d\overline{z}^{}}{z^{}z}(|u|^2)(z^{},\overline{z}^{}).$$

We deduce from Eqs. (2.25) and (2.27) that the potential $`\overline{u}(z,\overline{z})`$ is expressed through the eigenfunctions of the Dirac system in the form ,

$$\overline{u}(z,\overline{z})=\underset{j=1}{\overset{n}{}}\mathrm{\Phi }_{2j}(z,\overline{z})\frac{1}{2\pi }𝑑kd\overline{k}b(k,\overline{k})N_{2\mu }(z,\overline{z},k,\overline{k}).$$
(2.28)   
Formulas (2.6) to (2.28) constitute a standard framework for the inverse scattering transform of the DSII equation with a new relation (2.21). The existence and uniqueness of solutions of the Fredholm integral equations (2.6) and (2.7) and the $`\overline{}`$ problem (2.8) and (2.25) were proved in  and  under the small-norm assumption for the potential $`u(x,y)`$,

$$\left(\underset{(x,y)𝐑^2}{sup}|u|(x,y)\right)\left(|u(x,y)|𝑑x𝑑y\right)<\frac{\pi }{8}.$$

In this case $`n=0`$ and $`b(k,\overline{k})0`$. The nonlinear two-dimensional Fourier transform associated to this scheme was discussed in Examples 8-10 of Chapter 7.7 of Ref. . Indeed, the connection formula (2.28) implies that there is a scalar spectral decomposition of $`\overline{u}(z,\overline{z})`$ through $`N_{2\mu }(z,\overline{z},k,\overline{k})`$ for $`n=0`$. In order to close the decomposition, one could use Eqs. (2.9) and (2.13) to construct a “completeness relation” for the expansion of $`\delta (z^{}z)`$ in the form,

$$\delta (z^{}z)=\frac{1}{2\pi ^2i}𝑑kd\overline{k}\overline{N}_{2\mu }(z^{},\overline{z}^{},k,\overline{k})N_{2\mu }(z,\overline{z},k,\overline{k}).$$

However, we show in Proposition 3.4 below that a completeness theorem for Eq. (2.2) is different and is based on the set of eigenfunctions of the adjoint Dirac system.

## 3 Basis for a Scalar Spectral Decomposition

In this section, we specify the adjoint problem for the Dirac system (2.2) and establish orthogonality and completeness relations.

### 3.1 The Adjoint System

The adjoint system for Eq. (2.2) is

$$\mu _{1z}^a=ik\mu _1^au\mu _2^a,\mu _{2\overline{z}}^a=\overline{u}\mu _1^a,$$
(3.1)   
which provides the balance equation,

$$i(k^{}k)\mu _1^a(k^{})\mu _2(k)=\frac{}{z}\left[\mu _1^a(k^{})\mu _2(k)\right]\frac{}{\overline{z}}\left[\mu _2^a(k^{})\mu _1(k)\right].$$
(3.2)   
The system (3.1) admits plane solutions $`𝝁^a(z,\overline{z},k,\overline{k})`$ and oscillatory-type solutions $`𝐍_\mu ^a(z,\overline{z},k,\overline{k})`$ with the boundary conditions,

$`\underset{|k|\mathrm{}}{lim}𝝁^a(z,\overline{z},k,\overline{k})`$ $`=`$ $`𝐞_2,`$ (3.3)   
$`\underset{|k|\mathrm{}}{lim}𝐍_\mu ^a(z,\overline{z},k,\overline{k})e^{i(kz+\overline{k}\overline{z})}`$ $`=`$ $`𝐞_1.`$ (3.4)   
The adjoint eigenfunctions $`𝐍_\mu ^a(z,\overline{z},k,\overline{k})`$ can be expressed through the Green functions,

$`N_{1\mu }^a(z,\overline{z})`$ $`=`$ $`e^{i(kz+\overline{k}\overline{z})}{\displaystyle \frac{1}{2\pi i}}{\displaystyle \frac{dz^{}d\overline{z}^{}}{\overline{z}^{}\overline{z}}(uN_{2\mu }^a)(z^{},\overline{z}^{})e^{ik(zz^{})+i\overline{k}(\overline{z}\overline{z}^{})}},`$ (3.5)   
$`N_{2\mu }^a(z,\overline{z})`$ $`=`$ $`{\displaystyle \frac{1}{2\pi i}}{\displaystyle \frac{dz^{}d\overline{z}^{}}{z^{}z}(\overline{u}N_{1\mu }^a)(z^{},\overline{z}^{})}.`$ (3.6)   
They are related to the adjoint eigenfunctions $`𝝁^a(z,\overline{z},k,\overline{k})`$ by the formula,

$$𝐍_\mu ^a(z,\overline{z},k,\overline{k})=𝝈\overline{𝝁}^a(z,\overline{z},k,\overline{k})e^{i(kz+\overline{k}\overline{z})}.$$
(3.7)   
Using this representation, we prove the following result.   

Lemma 3.1. The spectral data $`b(k,\overline{k})`$ is expressed in terms of the adjoint eigenfunctions as

$$b(k,\overline{k})=\frac{1}{2\pi }𝑑zd\overline{z}(\overline{u}N_{1\mu }^a)(z,\overline{z}).$$
(3.8)   
Proof. Multiplying Eq. (3.5) by $`\overline{u}\mu _1(k)`$, integrating over $`dzd\overline{z}`$ and using Eq. (2.7), we express $`b(k,\overline{k})`$ defined in Eq. (2.9) in the form,

$$b(k,\overline{k})=\frac{1}{2\pi }𝑑zd\overline{z}\left[\overline{u}\mu _1(k)N_{1\mu }^a(k)u\mu _2(k)N_{2\mu }^a(k)\right].$$
(3.9)   
On the other hand, multiplying Eq. (2.6) by $`uN_{1\mu }^a(k)`$, integrating over $`dzd\overline{z}`$, and using Eqs. (3.6) and (3.9), we get Eq. (3.8). $`\mathrm{}`$   

Suppose now that $`k=k_j`$ is an isolated double eigenvalue of Eq. (2.2) with the bound states $`𝚽_j(z,\overline{z})`$ and $`𝚽_j^{}(z,\overline{z})`$ given by Eqs. (2.15) – (2.20). Suppose also that $`k=k_j^a`$ is an eigenvalue of the adjoint system (3.1) with the adjoint bound states $`𝚽_j^a(z,\overline{z})`$ and $`𝚽_j^a(z,\overline{z})`$.   

Lemma 3.2. If $`k_j`$ is a double eigenvalue of the Dirac system (2.2), then $`k_j`$ is also a double eigenvalue of the adjoint system (3.1).

Proof. We use Eq. (3.2) with $`𝝁=𝚽_j(z,\overline{z})`$ and $`𝝁^a=𝚽_j^a(z,\overline{z})`$ at $`k=k_j`$ and $`k^{}=k_j^a`$ and integrate over $`dzd\overline{z}`$ with the help of Eq. (A.3) of Appendix A. The contour contribution of the integral vanishes due to the boundary conditions (2.17) and (3.12) and the resulting expression is

$$(k_j^ak_j)𝑑zd\overline{z}(\mathrm{\Phi }_{1j}^a\mathrm{\Phi }_{2j})(z,\overline{z})=0.$$

The relation $`k_j^a=k_j`$ follows from this formula if the integral is non-zero at $`k_j^a=k_j`$ (which is proved below in Eq. (3.22)). The other possibility is when $`k_j^ak_j`$ but $`\mathrm{\Phi }_{1j}^a`$ is orthogonal to $`\mathrm{\Phi }_{2j}`$. We do not consider such a non-generic situation. The other bound state $`𝚽_j^a`$ at $`k_j^a=k_j`$ can be defined using the symmetry relation (see Eq. (3.15) below). $`\mathrm{}`$   

The adjoint bound state $`𝚽_j^a(z,\overline{z})`$ solves the homogeneous equations,

$`\mathrm{\Phi }_{1j}^a(z,\overline{z})`$ $`=`$ $`{\displaystyle \frac{1}{2\pi i}}{\displaystyle \frac{dz^{}d\overline{z}^{}}{\overline{z}^{}\overline{z}}(u\mathrm{\Phi }_{2j}^a)(z^{},\overline{z}^{})e^{ik_j(zz^{})+i\overline{k}_j(\overline{z}\overline{z}^{})}},`$ (3.10)   
$`\mathrm{\Phi }_{2j}^a(z,\overline{z})`$ $`=`$ $`{\displaystyle \frac{1}{2\pi i}}{\displaystyle \frac{dz^{}d\overline{z}^{}}{z^{}z}(\overline{u}\mathrm{\Phi }_{1j}^a)(z^{},\overline{z}^{})}`$ (3.11)   
with the boundary condition as $`|z|\mathrm{}`$,

$$𝚽_j^a(z,\overline{z})\frac{𝐞_2}{z},$$
(3.12)   
and the normalization conditions,

$`{\displaystyle \frac{1}{2\pi i}}{\displaystyle 𝑑z}d\overline{z}(\overline{u}\mathrm{\Phi }_{1j}^a)(z,\overline{z})`$ $`=`$ $`1,`$ (3.13)   
$`{\displaystyle \frac{1}{2\pi i}}{\displaystyle 𝑑z}d\overline{z}(u\mathrm{\Phi }_{2j}^a)(z,\overline{z})e^{i(k_jz+\overline{k}_j\overline{z})}`$ $`=`$ $`0.`$ (3.14)   
In addition, the bound state $`𝚽_j^a(z,\overline{z})`$ is related to $`𝚽(z,\overline{z})`$ according to the symmetry formula,

$$𝚽_j^a(z,\overline{z})=𝝈\overline{𝚽}_j^a(z,\overline{z})e^{i(k_jz+\overline{k}_j\overline{z})}.$$
(3.15)   
Using Eqs. (3.1)–(3.15), we see that the adjoint eigenfunction $`𝝁^a(z,\overline{z},k,\overline{k})`$ satisfies relations similar to those for $`𝝁(z,\overline{z},k,\overline{k})`$,

$$\frac{𝝁^a}{\overline{k}}=\overline{b}(k,\overline{k})𝐍_\mu ^a(z,\overline{z},k,\overline{k})$$
(3.16)   
and

$$\underset{kk_j}{lim}\left[𝝁^a(z,\overline{z},k,\overline{k})+\frac{i𝚽_j^a(z,\overline{z})}{kk_j}\right]=(z+z_j)𝚽_j^a(z,\overline{z})\overline{c}_j𝚽_j^a(z,\overline{z}).$$
(3.17)   
The expansions for inverse scattering transform of the adjoint eigenfunctions can be found in the form,

$$𝝁^a(z,\overline{z},k,\overline{k})=𝐞_2\underset{j=1}{\overset{n}{}}\frac{i𝚽_j^a(z,\overline{z})}{kk_j}\frac{1}{2\pi i}\frac{dk^{}d\overline{k}^{}}{k^{}k}\overline{b}(k^{},\overline{k}^{})𝐍_\mu ^a(z,\overline{z},k^{},\overline{k}^{})$$
(3.18)   
and

$$(z^{}+z_j)𝚽_j^a(z^{},\overline{z}^{})\overline{c}_j𝚽_j^a(z^{},\overline{z}^{})=𝐞_2\underset{lj}{}\frac{i𝚽_l^a(z^{},\overline{z}^{})}{k_jk_l}\frac{1}{2\pi i}\frac{dkd\overline{k}}{kk_j}\overline{b}(k,\overline{k})𝐍_\mu ^a(z,\overline{z},k,\overline{k}).$$
(3.19)   
### 3.2 Orthogonality and Completeness Relations

Using the Dirac system (2.2) and its adjoint system (3.1), we prove the orthogonality and completeness relations for the set of eigenfunctions $`S=[N_{2\mu }(k,\overline{k}),\{\mathrm{\Phi }_{2j}\}_{j=1}^n]`$ and its adjoint set $`S^a=[N_{1\mu }^a(k,\overline{k}),\{\mathrm{\Phi }_{1j}^a\}_{j=1}^n]`$.   

Proposition 3.3. The eigenfunctions $`N_{2\mu }(z,\overline{z},k,\overline{k})`$ and $`\mathrm{\Phi }_{2j}(z,\overline{z})`$ are orthogonal to the eigenfunctions $`N_{1\mu }^a(z,\overline{z},k,\overline{k})`$ and $`\mathrm{\Phi }_{1j}^a(z,\overline{z})`$ as follows

$$N_{1\mu }^a(k^{})|N_{2\mu }(k)_z=2\pi ^2i\delta (k^{}k),$$
(3.20)   
$$N_{1\mu }^a(k)|\mathrm{\Phi }_{2j}_z=\mathrm{\Phi }_{1j}^a|N_{2\mu }(k)_z=0,$$
(3.21)   
$$\mathrm{\Phi }_{1l}^a|\mathrm{\Phi }_{2j}_z=2\pi i\delta _{jl},$$
(3.22)   
where the inner product is defined as

$$g(k^{})|f(k)_z=𝑑zd\overline{z}g(z,\overline{z},k^{},\overline{k}^{})f(z,\overline{z},k,\overline{k}).$$

Proof. Using Eqs. (2.12) and (3.5), we expand the inner product in Eq. (3.20) as

$$N_{1\mu }^a(k^{})|N_{2\mu }(k)_z=I_0$$

$$+𝑑zd\overline{z}\left(uN_{2\mu }^a\right)(z,k^{})e^{i(k^{}z+\overline{k}^{}\overline{z})}I_1(z)𝑑zd\overline{z}\left(\overline{u}N_{1\mu }\right)(z,k)e^{i(kz+\overline{k}\overline{z})}I_1(z)$$

$$+\frac{1}{2\pi i}𝑑zd\overline{z}\left(uN_{2\mu }^a\right)(z,k^{})e^{i(k^{}z+\overline{k}^{}\overline{z})}\frac{dz^{}d\overline{z}^{}}{\overline{z}^{}\overline{z}}\left(\overline{u}N_{1\mu }\right)(z^{},k)e^{i(kz^{}+\overline{k}\overline{z}^{})}\left[I_1(z)I_1(z^{})\right],$$

where

$$I_0=𝑑zd\overline{z}e^{i(k^{}k)z+i(\overline{k}^{}\overline{k})\overline{z}}=2\pi ^2i\delta (k^{}k)$$
(3.23)   
and

$$I_1(z)=\frac{1}{2\pi i}\frac{dz^{}d\overline{z}^{}}{\overline{z}^{}\overline{z}}e^{i(k^{}k)z^{}+i(\overline{k}^{}\overline{k})\overline{z}^{}}=\frac{1}{i(k^{}k)}e^{i(k^{}k)z+i(\overline{k}^{}\overline{k})\overline{z}}.$$
(3.24)   
The integrals $`I_0`$ and $`I_1(z)`$ are computed in Appendix A. Using these formulas, we find the inner product in Eq. (3.20) in the form,

$$N_{1\mu }^a(k^{})|N_{2\mu }(k)_z=2\pi ^2i\delta (k^{}k)+\frac{1}{i(k^{}k)}(k,k^{}),$$

where the residual term $`(k,k^{})`$ is expressed in the form,

$$(k,k^{})=𝑑zd\overline{z}\left[uN_{2\mu }^a(k^{})N_{2\mu }(k)\overline{u}N_{1\mu }^a(k^{})N_{1\mu }(k)\right],$$

with the help of Eqs. (2.12) and (3.5). We show that $`(k,k^{})=0`$ by multiplying Eq. (3.6) by $`\overline{u}N_{1\mu }(k)`$, integrating over $`dzd\overline{z}`$ and using Eq. (2.11).

The zero inner products in Eqs. (3.21) and (3.22) for $`jl`$ are obtained in a similar way with the help of the Fredholm’s equations for eigenfunctions $`𝚽_j`$, $`𝚽_j^a`$, $`𝐍_\mu `$, and $`𝐍_\mu ^a`$. In order to find the non-zero inner product (3.22) for $`j=l`$ we evaluate the following integral by using the same integral equations,

$`{\displaystyle 𝑑z}d\overline{z}\mathrm{\Phi }_{1j}^a\mu _2(k)`$ $`=`$ $`{\displaystyle \frac{1}{i(k_jk)}}{\displaystyle 𝑑z}d\overline{z}\left[u\mathrm{\Phi }_{2j}^a\mu _2(k)\overline{u}\mathrm{\Phi }_{1j}^a\mu _1(k)\right]`$   
$`=`$ $`{\displaystyle \frac{1}{i(kk_j)}}{\displaystyle 𝑑z}d\overline{z}\overline{u}\mathrm{\Phi }_{1j}^a.`$   
Using Eq. (3.13), the right-hand-side identifies to $`\frac{1}{i(kk_j)}`$. Substituting Eq. (2.25) in the left-hand-side and using the zero inner products (3.21) and (3.22), we find Eq. (3.22) for $`j=l`$. $`\mathrm{}`$   

Proposition 3.4. The eigenfunctions $`N_{2\mu }(z,\overline{z},k,\overline{k})`$ and $`\mathrm{\Phi }_{2j}(z,\overline{z})`$ are complete with respect to the adjoint eigenfunctions $`N_{1\mu }^a(z,\overline{z},k,\overline{k})`$ and $`\mathrm{\Phi }_{1j}^a(z,\overline{z})`$ according to the identity,

$$\delta (z^{}z)=\frac{1}{2\pi ^2i}𝑑kd\overline{k}N_{1\mu }^a(z^{},\overline{z}^{},k,\overline{k})N_{2\mu }(z,\overline{z},k,\overline{k})\frac{1}{\pi }\underset{j=1}{\overset{n}{}}\mathrm{\Phi }_{1j}^a(z^{},\overline{z}^{})\mathrm{\Phi }_{2j}(z,\overline{z}).$$
(3.25)   
Proof. Using the symmetry relations (2.13) and (3.7), we express the integral in Eq. (3.25) as

$$N_{1\mu }^a(z^{})|N_{2\mu }(z)_k=𝑑kd\overline{k}\overline{\mu }_2^a(z^{},\overline{z}^{},k,\overline{k})\overline{\mu }_1(z,\overline{z},k,\overline{k})e^{ik(z^{}z)+i\overline{k}(\overline{z}^{}\overline{z})}.$$
(3.26)   
We use Eqs. (2.25), (2.26), (3.18), and (3.19) and find the pole decomposition for the integrand in Eq. (3.26),

$$\overline{\mu }_2^a(z^{})\overline{\mu }_1(z)=1+\underset{j=1}{\overset{n}{}}\frac{i}{\overline{k}\overline{k}_j}\left[\overline{c}_j\overline{\mathrm{\Phi }}_{2j}^a(z^{},\overline{z}^{})\overline{\mathrm{\Phi }}_{1j}^{}(z,\overline{z})+c_j\overline{\mathrm{\Phi }}_{2j}^a(z^{},\overline{z}^{})\overline{\mathrm{\Phi }}_{1j}(z,\overline{z})\right]$$

$$+\underset{j=1}{\overset{n}{}}\frac{i}{\overline{k}\overline{k}_j}(\overline{z}\overline{z}^{})\overline{\mathrm{\Phi }}_{2j}^a(z^{},\overline{z}^{})\overline{\mathrm{\Phi }}_{1j}(z,\overline{z})+\underset{j=1}{\overset{n}{}}\frac{\overline{\mathrm{\Phi }}_{2j}^a(z^{},\overline{z}^{})\overline{\mathrm{\Phi }}_{1j}(z,\overline{z})}{(\overline{k}\overline{k}_j)^2}$$

$$+\frac{1}{2\pi i}\frac{dk^{}d\overline{k}^{}}{\overline{k}^{}\overline{k}}\left[\overline{\mu }_2^a(z^{})\overline{b}\overline{N}_{1\mu }(z)\overline{\mu }_1(z)b\overline{N}_{2\mu }^a(z^{})\right](k^{},\overline{k}^{}).$$
(3.27)   
We substitute (3.27) into Eq. (3.26) and reduce the integral to the form,

$$N_{1\mu }^a(z^{})|N_{2\mu }(z)_k=I_02\pi \underset{j=1}{\overset{n}{}}\left[\overline{c}_j\overline{\mathrm{\Phi }}_{2j}^a(z^{},\overline{z}^{})\overline{\mathrm{\Phi }}_{1j}^{}(z,\overline{z})+c_j\overline{\mathrm{\Phi }}_{2j}^a(z^{},\overline{z}^{})\overline{\mathrm{\Phi }}_{1j}(z,\overline{z})\right]I_1(k_j)$$

$$+\underset{j=1}{\overset{n}{}}\overline{\mathrm{\Phi }}_{2j}^a(z^{},\overline{z}^{})\overline{\mathrm{\Phi }}_{1j}(z,\overline{z})\left[2\pi (\overline{z}^{}\overline{z})I_1(k_j)+I_2(k_j)\right]$$

$$𝑑kd\overline{k}I_1(k)\left[\overline{\mu }_2^a(z^{})\overline{b}\overline{N}_{1\mu }(z)\overline{\mu }_1(z)b\overline{N}_{2\mu }^a(z^{})\right](k,\overline{k})e^{ik(z^{}z)i\overline{k}(\overline{z}^{}\overline{z})},$$
(3.28)   
where the integrals $`I_0`$ and $`I_1(k)`$ are given in Eqs. (3.23) and (3.24) respectively, with $`z`$ and $`k`$ interchanged, while the integral $`I_2(k_j)`$ is defined by

$$I_2(k_j)=\underset{ϵ0}{lim}_{|kk_j|ϵ}\frac{dkd\overline{k}}{(\overline{k}\overline{k}_j)^2}e^{ik(z^{}z)+i\overline{k}(\overline{z}^{}\overline{z})}.$$
(3.29)   
The integral $`I_2(k)`$ is found in Appendix A in the form, $`I_2(k_j)=2\pi (\overline{z}^{}\overline{z})I_1(k_j)`$, such that the third term in Eq. (3.28) vanishes. In order to express the second term in Eq. (3.28) we use Eqs. (2.26), (3.19), (2.20), and (3.15) and derive the relation,

$$\underset{j=1}{\overset{n}{}}\left[\overline{c}_j\overline{\mathrm{\Phi }}_{2j}^a(z^{},\overline{z}^{})\overline{\mathrm{\Phi }}_{1j}^{}(z,\overline{z})+c_j\overline{\mathrm{\Phi }}_{2j}^a(z^{}\overline{z}^{})\overline{\mathrm{\Phi }}_{1j}(z,\overline{z})\right]e^{ik_j(z^{}z)+i\overline{k}_j(\overline{z}^{}\overline{z})}=$$

$$\underset{j=1}{\overset{n}{}}\left[\overline{c}_j\mathrm{\Phi }_{1j}^a(z^{},\overline{z}^{})\mathrm{\Phi }_{2j}(z,\overline{z})+c_j\mathrm{\Phi }_{1j}^a(z^{},\overline{z}^{})\mathrm{\Phi }_{2j}^{}(z,\overline{z})\right]=$$

$$(z^{}z)\underset{j=1}{\overset{n}{}}\mathrm{\Phi }_{1j}^a(z^{},\overline{z}^{})\mathrm{\Phi }_{2j}(z,\overline{z})+\frac{1}{2\pi i}\underset{j=1}{\overset{n}{}}\frac{dkd\overline{k}}{kk_j}\left[\mathrm{\Phi }_{2j}(z)\overline{b}N_{1\mu }^a(z^{})+\mathrm{\Phi }_{1j}^a(z^{})bN_{2\mu }(z)\right](k,\overline{k}).$$
(3.30)   
Using this expression and Eqs. (3.23), (3.24), and (3.29), we rewrite Eq. (3.28) in the form,

$$N_{1\mu }^a(z^{})|N_{2\mu }(z)_k=2\pi ^2i\delta (z^{}z)2\pi i\underset{j=1}{\overset{n}{}}\mathrm{\Phi }_{1j}^a(z^{},\overline{z}^{})\mathrm{\Phi }_{2j}(z,\overline{z})$$

$$+\frac{i}{z^{}z}𝑑kd\overline{k}\frac{}{\overline{k}}\left[\left(\mu _1^a(z^{})+\underset{j=1}{\overset{n}{}}\frac{i\mathrm{\Phi }_{1j}^a(z^{})}{kk_j}\right)\left(\mu _2(z)\underset{j=1}{\overset{n}{}}\frac{i\mathrm{\Phi }_{2j}(z)}{kk_j}\right)\right].$$

The last integral vanishes according to Eq. (A.3) of Appendix A and the boundary conditions (2.5) and (3.3). $`\mathrm{}`$   

Our main result for the spectral decomposition associated to the Dirac system (2.2) follows from the above orthogonality and completeness relations.   

Proposition 3.5. An arbitrary scalar function $`f(z,\overline{z})`$ satisfying the condition $`f(z,\overline{z})\mathrm{O}(|z|^2)`$ as $`|z|\mathrm{}`$ can be decomposed through the set $`S=[N_{2\mu }(z,\overline{z},k,\overline{k}),\{\mathrm{\Phi }_{2j}(z,\overline{z})\}_{j=1}^n]`$.   

Proof. The spectral decomposition is defined through the orthogonality relations (3.20)–(3.22) as

$$f(z,\overline{z})=𝑑kd\overline{k}\alpha (k,\overline{k})N_{2\mu }(z,\overline{z},k,\overline{k})+\underset{j=1}{\overset{n}{}}\alpha _j\mathrm{\Phi }_{2j}(z,\overline{z}),$$
(3.31)   
where

$$\alpha (k,\overline{k})=\frac{1}{4\pi ^2}N_{1\mu }^a(k)|f_z,\alpha _j=\frac{1}{2\pi i}\mathrm{\Phi }_{1j}^a|f_z.$$
(3.32)   
Provided the condition on $`f(z,\overline{z})`$ is satisfied, we interchange integration with respect to $`dzd\overline{z}`$ and $`dkd\overline{k}`$ and use the completeness formula (3.25). $`\mathrm{}`$   

The spectral decomposition presented here is different from that of Kiselev . In the latter approach, the function $`f(z,\overline{z})`$ is spanned by squared eigenfunctions of the original problem (1.2) defined according to oscillatory-type behaviour at infinity. In our approach, we transformed the system (1.2) to the form (2.2) and defined the oscillatory-type eigenfunctions according to the single eigenfunctions $`𝐍_\mu (z,\overline{z},k,\overline{k})`$. We also notice that the (degenerate) bound states $`𝚽_j^{}(z,\overline{z})`$ are not relevant for the spectral decomposition, although they appear implicitly through the meromorphic contributions of the eigenfunctions $`𝐍_\mu (z,\overline{z},k,\overline{k})`$ at $`k=k_j`$ (see Section 4).

## 4 Perturbation Theory for a Single Lump

We use the scalar spectral decomposition based on Eq. (3.31) and develop a perturbation theory for multi-lump solutions of the DSII equation. We present formulas in the case of a single lump ($`n=1`$), the case of multi-lump potentials can be obtained by summing along the indices $`j`$, $`l`$ occuring in the expressions below.

The single-lump potential $`u(z,\overline{z})`$ has the form ,

$$u(z,\overline{z})=\frac{c_j}{|z+z_j|^2+|c_j|^2}e^{i(k_jz+\overline{k}_j\overline{z})},$$
(4.1)   
where $`c_j`$, $`z_j`$ are complex parameters. The associated bound states follow from Eqs. (2.26) and (3.19) as

$`𝚽_j(z,\overline{z})`$ $`=`$ $`{\displaystyle \frac{1}{|z+z_j|^2+|c_j|^2}}\left[\begin{array}{c}\overline{z}+\overline{z}_j\\ \overline{c}_je^{i(k_jz+\overline{k}_j\overline{z})}\end{array}\right],`$ (4.4)   
$`𝚽_j^a(z,\overline{z})`$ $`=`$ $`{\displaystyle \frac{1}{|z+z_j|^2+|c_j|^2}}\left[\begin{array}{c}c_je^{i(k_jz+\overline{k}_j\overline{z})}\\ \overline{z}+\overline{z}_j\end{array}\right].`$ (4.7)   
We first consider a general perturbation to the single lump subject to the localization condition, $`\mathrm{\Delta }u\mathrm{O}(|z|^2)`$ as $`|z|\mathrm{}`$. We then derive explicit formulas for a special form of the perturbation term $`\mathrm{\Delta }u(z,\overline{z})`$.

### 4.1 General Perturbation of a Single Lump

Suppose the potential is specified as $`u^ϵ=u(z,\overline{z})+ϵ\mathrm{\Delta }u(z,\overline{z})`$, where $`u(z,\overline{z})`$ is given by Eq. (4.1) and $`\mathrm{\Delta }u(z,\overline{z})`$ is a perturbation term. Two bound states $`𝚽_j(z,\overline{z})`$ and $`𝚽_j^{}(z,\overline{z})`$ are supported by a single-lump potential $`u(z,\overline{z})`$ at a single point $`k=k_j`$. The spectral decomposition given by Eq. (3.31) provides a basis for expansion of $`\mu _2^ϵ(z,\overline{z},\kappa ,\overline{\kappa })`$ at $`k=\kappa `$,

$$\mu _2^ϵ(z,\overline{z},\kappa ,\overline{\kappa })=𝑑kd\overline{k}\alpha (k,\overline{k})N_{2\mu }(z,\overline{z},k,\overline{k})+\alpha _j\mathrm{\Phi }_{2j}(z,\overline{z})$$
(4.8)   
where $`\alpha (k,\overline{k})`$ and $`\alpha _j`$ are defined by Eq. (3.32) and depend on the parameter $`\kappa `$. The other component $`\mu _1^ϵ(z,\overline{z},\kappa ,\overline{\kappa })`$ can be expressed from Eq. (2.2) as

$$\mu _1^ϵ(z,\overline{z},\kappa ,\overline{\kappa })=𝑑kd\overline{k}\alpha (k,\overline{k})N_{1\mu }(z,\overline{z},k,\overline{k})+\alpha _j\mathrm{\Phi }_{1j}(z,\overline{z})+ϵ\mathrm{\Delta }\mu _1(z,\overline{z}),$$
(4.9)   
where the remainder term $`\mathrm{\Delta }\mu _1(z,\overline{z})`$ solves the equation,

$$\left(\mathrm{\Delta }\mu _1\right)_{\overline{z}}=\mathrm{\Delta }u\mu _2^ϵ.$$

We write the solution of this equation in the form,

$$\mathrm{\Delta }\mu _1(z,\overline{z})=A\frac{1}{2\pi i}\frac{dz^{}d\overline{z}^{}}{z^{}z}\left(\mathrm{\Delta }u\mu _2^ϵ\right)(z^{},\overline{z}^{}),$$
(4.10)   
subject to the boundary condition as $`|z|\mathrm{}`$,

$$\mathrm{\Delta }\mu _1(z,\overline{z})A+\mathrm{O}(z^1),$$

where $`A`$ is an arbitrary constant. Using the explicit representation (4.10), we transform Eq. (2.2) into the system of integral equations for $`\alpha (k,\overline{k})`$ and $`\alpha _j`$,

$$\alpha (k,\overline{k})=\frac{ϵ}{4\pi ^2i(k\kappa )}\left[𝑑k^{}d\overline{k}^{}K(k,\overline{k},k^{},\overline{k}^{})\alpha (k^{},\overline{k}^{})+K_j(k,\overline{k})\alpha _j+R(k,\overline{k})A\right]+\mathrm{O}(ϵ^2),$$
(4.11)   
$$\alpha _j=\frac{ϵ}{2\pi (k_j\kappa )}\left[𝑑kd\overline{k}P_j(k,\overline{k})\alpha (k,\overline{k})+K_{jl}\alpha _l+R_jA\right]+\mathrm{O}(ϵ^2),$$
(4.12)   
where

$$K(k,\overline{k},k^{},\overline{k}^{})=𝐍_\mu ^a(k)|𝐍_\mu (k^{})_{\mathrm{\Delta }u},K_j(k,\overline{k})=𝐍_\mu ^a(k)|𝚽_j_{\mathrm{\Delta }u},$$

$$P_j(k,\overline{k})=𝚽_j^a|𝐍_\mu (k)_{\mathrm{\Delta }u},K_{jl}=𝚽_j^a|𝚽_l_{\mathrm{\Delta }u},$$

and the scalar product for squared eigenfunction is defined as 

$$𝐟(k)|𝐠(k^{})_h=𝑑zd\overline{z}\left[\overline{h}(z,\overline{z})f_1(z,\overline{z},k,\overline{k})g_1(z,\overline{z},k^{},\overline{k}^{})+h(z,\overline{z})f_2(z,\overline{z},k,\overline{k})g_2(z,\overline{z},k^{},\overline{k}^{})\right].$$

The non-homogeneous terms $`R(k,\overline{k})`$ and $`R_j`$ can be computed exactly as

$`R(k,\overline{k})`$ $`=`$ $`{\displaystyle 𝑑z}d\overline{z}\left(\overline{u}N_{1\mu }^a\right)(z,\overline{z})=2\pi b(k,\overline{k}),`$   
$`R_j`$ $`=`$ $`{\displaystyle 𝑑z}d\overline{z}\left(\overline{u}\mathrm{\Phi }_{1j}^a\right)(z,\overline{z})=2\pi i,`$   
where $`b(k,\overline{k})=0`$ if $`n0`$. We solve the system of equations (4.11) and (4.12) asymptotically for $`\kappa =k_j+ϵ\mathrm{\Delta }\kappa `$ and $`\mathrm{\Delta }\kappa \mathrm{O}(1)`$. The leading order behaviour of the integral kernels follows from the asymptotic representation (2.24) as $`kk_j`$,

$$K(k,\overline{k},k^{},\overline{k}^{})\frac{\overline{K}_{jj}}{(\overline{k}\overline{k}_j)(\overline{k}^{}\overline{k}_j)},K_j(k,\overline{k})\frac{iP_{jj}}{\overline{k}\overline{k}_j},P_j(k,\overline{k})\frac{i\overline{P}_{jj}}{\overline{k}\overline{k}_j},$$
(4.13)   
where

$$\overline{K}_{jj}=𝚽_j^a|𝚽_j^{}_{\mathrm{\Delta }u},P_{jj}=𝚽_j^a|𝚽_j_{\mathrm{\Delta }u},\overline{P}_{jj}=𝚽_j^a|𝚽_j^{}_{\mathrm{\Delta }u}.$$
(4.14)   
Here we have used the symmetry constraints (2.20) and (3.15). The leading order of $`\alpha (k,\overline{k})`$ as $`kk_j`$ follows from Eq. (4.11) as

$$\alpha (k,\overline{k})\frac{ϵ\mathrm{\Delta }\overline{\kappa }\beta _j}{2\pi (k\kappa )(\overline{k}\overline{k}_j)},$$

where $`\beta _j`$ is not yet defined. We use Eq. (A.4) of Appendix A to compute the integral term,

$$\frac{dkd\overline{k}}{(k\kappa )(\overline{k}\overline{k}_j)^2}=\frac{2\pi i}{\overline{k}_j\overline{\kappa }},$$

and reduce the system of integral equations (4.11) and (4.12) to an algebraic system as $`kk_j`$,

$`2\pi \mathrm{\Delta }\kappa \alpha _j`$ $`=`$ $`K_{jj}\alpha _j\overline{P}_{jj}\beta _j2\pi iA,`$ (4.15)   
$`2\pi \mathrm{\Delta }\overline{\kappa }\beta _j`$ $`=`$ $`P_{jj}\alpha _j\overline{K}_{jj}\beta _j.`$ (4.16)   
If $`P_{jj}0`$, the determinant of the above system is strictly positive. Therefore, homogeneous solutions at $`A=0`$ (bound states) are absent for $`ϵ0`$. This result indicates that the double eigenvalue at $`k=k_j`$ disappears under a generic perturbation of the potential $`u(z,\overline{z})`$ with $`P_{jj}0`$ (see also Ref. ).

For $`A0`$, we find inhomogeneous solutions of Eqs. (4.15) and (4.16),

$$\alpha _j=\frac{2\pi iA\left(\overline{K}_{jj}+2\pi \mathrm{\Delta }\overline{\kappa }\right)}{|K_{jj}+2\pi \mathrm{\Delta }\kappa |^2+|P_{jj}|^2},\beta _j=\frac{2\pi iAP_{jj}}{|K_{jj}+2\pi \mathrm{\Delta }\kappa |^2+|P_{jj}|^2}.$$
(4.17)   
The eigenfunction $`𝝁^ϵ(z,\overline{z},\kappa ,\overline{\kappa })`$ given by Eqs. (4.8) and (4.9) satisfies the boundary condition (2.5) if $`A=ϵ^1`$ and has the following asymptotic representation,

$`𝝁^ϵ(z,\overline{z},\kappa ,\overline{\kappa })=𝐞_1+{\displaystyle \frac{2\pi i\left[2\pi (\overline{\kappa }\overline{k}_j)+ϵ\overline{K}_{jj}\right]𝚽_j(z,\overline{z})}{|2\pi (\kappa k_j)+ϵK_{jj}|^2+|ϵP_{jj}|^2}}`$ $``$ $`{\displaystyle \frac{2\pi iϵP_{jj}𝚽_j^{}(z,\overline{z})}{|2\pi (\kappa k_j)+ϵK_{jj}|^2+|ϵP_{jj}|^2}}`$ (4.18)   
$`+`$ $`\mathrm{\Delta }𝝁^ϵ(z,\overline{z}),`$   
where the term $`\mathrm{\Delta }𝝁^ϵ(z,\overline{z})`$ is not singular in the limit $`ϵ0`$ and $`\kappa k_j`$.

In the limit $`ϵ0`$, $`\kappa k_j`$, we find a meromorphic expansion for $`\mu ^ϵ(z,\overline{z},\kappa ,\overline{\kappa })`$ as

$$𝝁^ϵ(z,\overline{z},\kappa ,\overline{\kappa })=𝐞_1+\frac{i𝚽_j(z,\overline{z})}{\kappa k_j}+ϵ\left[\frac{K_{jj}𝚽_j(z,\overline{z})}{2\pi i(\kappa k_j)^2}+\frac{P_{jj}𝚽_j^{}(z,\overline{z})}{2\pi i|\kappa k_j|^2}\right]+\mathrm{O}(ϵ^2).$$
(4.19)   
It is clear that the double pole can be incorporated by shifting the eigenvalue $`k_j`$ to

$$k_j^ϵ=k_j\frac{ϵK_{jj}}{2\pi }.$$

The other double-pole term in the expansion (4.19) has a non-analytic behaviour in the $`k`$-plane and leads to the appearance of the spectral data $`b^ϵ(\kappa ,\overline{\kappa })=ϵ\mathrm{\Delta }b(\kappa ,\overline{\kappa })`$ which measures the departure of $`𝝁^ϵ(z,\overline{z},\kappa ,\overline{\kappa })`$ from analyticity according to Eq. (2.8). We find from Eqs. (2.9) and (4.19) that the spectral data $`\mathrm{\Delta }b(\kappa ,\overline{\kappa })`$ has the following singular behaviour as $`\kappa k_j`$

$$\mathrm{\Delta }b(\kappa ,\overline{\kappa })\frac{P_{jj}}{2\pi |\kappa k_j|^2}.$$
(4.20)   
Thus, if $`P_{jj}0`$ the analyticity of $`𝝁^ϵ(z,\overline{z},\kappa ,\overline{\kappa })`$ is destroyed and the lump disappears. This conclusion as well as the analytical solution (4.17) agree with the results of Gadyl’shin and Kiselev  where the transformation of a single lump into decaying wave packet was also studied.

In the other limit $`ϵ0`$ and $`\kappa k_j^ϵ`$ we find another expansion from Eq. (4.18),

$$𝝁^ϵ(z,\overline{z},\kappa ,\overline{\kappa })=𝐞_1\frac{2\pi i}{ϵ\overline{P}_{jj}}𝚽_j^{}(z,\overline{z})+\mathrm{O}(\kappa k_j^ϵ).$$
(4.21)   
We conclude that the eigenfunction $`𝝁^ϵ(z,\overline{z},\kappa ,\overline{\kappa })`$ is now free of pole singularities . We summarize the main result in the form of a proposition.   

Proposition 4.1. Suppose $`u(z,\overline{z})`$ is given by Eq. (4.1) and $`\mathrm{\Delta }u(z,\overline{z})`$ satisfies the constraint,

$$P_{jj}=𝚽_j^a|𝚽_j_{\mathrm{\Delta }u}0.$$

Then, the potential $`u^ϵ=u(z,\overline{z})+ϵ\mathrm{\Delta }u(z,\overline{z})`$ does not support embedded eigenvalues of the Dirac system (2.2) for $`ϵ0`$.   

### 4.2 Explicit Solution for a Particular Perturbation

Here we specify $`c_j=ce^{i\theta }`$, where $`c`$ and $`\theta `$ are real, and consider a particular perturbation $`\mathrm{\Delta }u(z,\overline{z})`$ to the lump $`u(z,\overline{z})`$ (4.1) in the form,

$$\mathrm{\Delta }u(z,\overline{z})=Q(z,\overline{z})e^{i(k_jz+\overline{k}_j\overline{z}+\theta )},$$

where $`Q(z,\overline{z})`$ is a real function. Using Eqs. (4.14), (4.4), and (4.7), we find explicitly the matrix elements $`K_{jj}`$ and $`P_{jj}`$,

$$K_{jj}=𝑑zd\overline{z}\frac{c(\overline{z}+\overline{z}_j)[Q(z,\overline{z})\overline{Q}(z,\overline{z})]}{[|z+z_j|^2+c^2]^2}=0,$$

$$P_{jj}=𝑑zd\overline{z}\frac{|z+z_j|^2\overline{Q}(z,\overline{z})+c^2Q(z,\overline{z})}{[|z+z_j|^2+c^2]^2}=\frac{1}{2c}𝑑zd\overline{z}\left(u\mathrm{\Delta }\overline{u}+\overline{u}\mathrm{\Delta }u\right).$$

The element $`P_{jj}`$ can be seen as a correction to the field energy,

$$N=\frac{i}{2}𝑑zd\overline{z}|u^ϵ|^2(z,\overline{z})=N_0+iϵcP_{jj}+\mathrm{O}(ϵ^2),$$

where $`N_0=\pi `$ is the energy of the single lump solution (independent of the lump parameters $`k_j`$ and $`c_j`$). Thus, a perturbation which leads to the destruction of a single lump, that is with $`P_{jj}0`$, changes necessarily the value for the lump energy $`N_0`$.

## 5 Concluding Remarks

The main result of our paper is the prediction of structural instability of multi-lump potentials in the Dirac system associated to the DSII equation. The multi-lump potentials correspond to eigenvalues embedded into a two-dimensional continuous spectrum with the spectral data $`b(k,\overline{k})`$ satisfying the additional constraint (2.21). In this case, there is no interaction between lumps and continuous radiation. However, a generic initial perturbation induces coupling between the lumps and radiation and, as a result of their interaction, the embedded eigenvalues disappear. This result indicates that the localized multi-lump solutions decay into continuous wave packets in the nonlinear dynamics of the DSII equation (see also Refs. ).

This scenario is different from the two types of bifurcations of embedded eigenvalues discussed in our previous paper . The type I bifurcation arises from the edge of the essential spectrum when the limiting bounded (non-localized) eigenfunction is transformed into a localized bound state. The type II bifurcation occurs when an embedded eigenvalue splits off the essential spectrum. Both situations persist in the spectral plane when the essential spectrum is one-dimensional and covers either a half-axis or the whole axis. However, in the case of the DSII equation, the essential spectrum is the whole spectral plane and embedded eigenvalue can not split off the essential spectrum. As a result, they disappear due to their structural instability.

## Acknowledgements

We benefited from stimulating discussions with M. Ablowitz, A. Fokas, D. Kaup, O. Kiselev, and A. Yurov. D.P. acknowledges support from a NATO fellowship provided by NSERC and C.S. acknowledges support from NSERC Operating grant OGP0046179.

## Appendix A. Formulas of the $``$-analysis

Here we reproduce some formulas of the complex $`\overline{}`$-analysis  to compute the integrals $`I_0`$, $`I_1(z)`$ and $`I_2(z)`$ defined in Eqs. (3.23), (3.24), and (3.29). We define the complex integration in the $`z`$-plane by

$$𝑑zd\overline{z}f(z,\overline{z})=𝑑\overline{z}dzf(z,\overline{z}),$$

where $`dzd\overline{z}=2idxdy`$. The complex $`\delta (z)`$ distribution is defined by

$$𝑑zd\overline{z}f(z)\delta (zz_0)=2if(z_0),$$
$`(A.1)`$   
where $`\delta (z)=\delta (x)\delta (y)`$. In particular, the $`\delta `$-distribution appears in the $`\overline{}`$-analysis according to the relation ,

$$\frac{}{\overline{z}}\left[\frac{1}{zz_0}\right]=\pi \delta (zz_0).$$
$`(A.2)`$   
Computing the integral $`I_0`$, we get the formula,

$$I_0=𝑑zd\overline{z}e^{ikz+i\overline{k}\overline{z}}=2i𝑑x𝑑ye^{2i\mathrm{Re}(k)x2i\mathrm{Im}(k)y}=2\pi ^2i\delta (k),$$

which proves the identity (3.23).

Using the Green’s theorem , one has the integration identity,

$$_D𝑑zd\overline{z}\left(\frac{f_1}{z}\frac{f_2}{\overline{z}}\right)=_C\left(f_1d\overline{z}+f_2dz\right),$$
$`(A.3)`$   
where $`D`$ is a domain of the complex plane and $`C`$ its boundary. The generalized Cauchy’s formula has the form ,

$$f(z,\overline{z})=\frac{1}{2\pi i}_C\frac{f(z^{},\overline{z}^{})dz^{}}{z^{}z}+\frac{1}{2\pi i}_D\frac{dz^{}d\overline{z}^{}}{z^{}z}\frac{f}{\overline{z}^{}},$$
$`(A.4)`$   
or, equivalently,

$$f(z,\overline{z})=\frac{1}{2\pi i}_C\frac{f(z^{},\overline{z}^{})d\overline{z}^{}}{\overline{z}^{}\overline{z}}+\frac{1}{2\pi i}_D\frac{dz^{}d\overline{z}^{}}{\overline{z}^{}\overline{z}}\frac{f}{z^{}}.$$
$`(A.5)`$   
In order to find the integral $`I_1(z)`$ we use Eq. (A.5) with

$$f(z,\overline{z})=\frac{1}{ik}e^{i(kz+\overline{k}\overline{z})},k0$$

and choose the domain $`D`$ to be a large ball of radius $`R`$ (see Eq. (2.14)). The boundary value integral vanishes since

$$\underset{R\mathrm{}}{lim}_{|z|=R}\frac{d\overline{z}}{\overline{z}}e^{i(kz+\overline{k}\overline{z})}=2\pi i\underset{R\mathrm{}}{lim}J_0(2|k|R)=0,$$
$`(A.6)`$   
where $`J_0(z)`$ is the Bessel function. Equation (A.5) for the function $`f(z,\overline{z})`$ then reduces to Eq. (3.24).

In order to compute the integral $`I_2(z_0)`$, we apply Eq. (A.3) with $`f_1=0`$ and

$$f_2(z,\overline{z})=\frac{1}{\overline{z}\overline{z}_0}e^{i(kz+\overline{k}\overline{z})}.$$

The domain $`D`$ is chosen as above. The boundary value integral vanishes again,

$$\underset{R\mathrm{}}{lim}_{|z|=R}\frac{dz}{\overline{z}}e^{i(kz+\overline{k}\overline{z})}=2\pi i\underset{R\mathrm{}}{lim}J_2(2|k|R)=0,$$
$`(A.7)`$   
where $`J_2(z)`$ is the Bessel function. Equation (A.3) for the function $`f_2(z,\overline{z})`$ reduces to Eq. (3.29).