Abstract

We consider a natural integrable generalization of nonlinear Schrödinger equation to dimensions. By studying the associated spectral operator we discover a rich discrete spectrum associated with regular rationally decaying solutions, the lumps, which display interesting nontrivial dynamics and scattering. Particular interest is placed in the dynamical evolution of the associated pulses. For all cases under study we find that the relevant dynamics corresponds to a central configuration of a certain -body problem.

1. The Physical System

In this paper we consider the system of equations:where is a complex function, depending on three real variables , , . We derive via the inverse scattering method several large classes of solutions that satisfy the boundary conditions (),

Note that in this case the system readsDue to the overall simplicity, elegance, and potential physical interest, the above system is a natural interesting nonlinear -dimensional integrable equation to study, even though it has not been related to any physical situation so far. To further motivate its study we note the following:(1)The more general problem corresponding to the boundary conditionswhere is a given nonnull constant is also solvable and reducible to problem (2) as follows. Let be the unique solution to (1) under the boundary conditions (2). Then one proves easily that also solves (1) and satisfies (4).(2)System (3) arises as the compatibility of a Lax pair (see Section 2 below) and hence it is integrable. We show here that the operator has also a discrete spectrum, corresponding to potentials that solve (1) and (2): that is, it can be written as , where is a regular, weakly decaying function.(3)If and only satisfy and , then (1) readsThis equation was first noticed to be integrable by Fokas [1] and also in Shul’man [2] who also embedded it into a general multiparameter family of equations. Note that in potentials satisfying there is not discrete spectrum (see [3, 4]).(4)The further reduction to the manifold yields the nonlinear Schrödinger (NLS) equation. Hence (1) generalizes NLS to dimensions.(5)The relationship between the previous problem and that corresponding to the boundary conditions (2) and where and , has been investigated in [5]. We note that the conserved quantities depend on the election of .(6)There is an infinite number of conservation laws for this equation (see [5]), which includes, in particular, the functional We note that the “amplitude” of interesting solutions to (3) is naturally expected to be the physical observable but it is not definite positive; hence, the renormalized mass of the field can take both signs.(7)An interesting parabolic real version of (1) (i.e., letting ) is obtained as a particular reduction of the Self-Dual Yang Mills Equations; see [6].

Lump configurations (i.e., regular and rationally decaying solutions) are paradigmatic solutions of integrable equations in dimensions and as such they have been extensively studied in the last years. Their spectral interpretation was first unlocked in relation to KPI [7] (see also [8, 9]). A description of the KP equation, its physical origins, integrability, and soliton solutions can be found in [10]. Subsequently, lumps have been found in other integrable equations like Davey-Stewartdon II (DSII) (see [11]) and the -Toda lattice (see [1214]).

The dynamics of standard lumps on the plane is trivial: the motion is uniform; further, upon interaction, only a parallel shift on the asymptotic trajectory is to be found. Such uniform motion has been considered to rule the dynamics of localized pulses of integrable equations up until the late nineties. Remarkably it was found that KPI possess a new class of localized, real valued solutions that have a nontrivial asymptotic dynamics. Even though the simplest of these solutions has been known for a long time, the interesting scattering properties that they exhibit went unnoticed up until 1995 [15, 16]. The spectral interpretation along with a derivation of the general class of those solutions was unlocked in [17, 18]. Nonstandard lumps were found to be associated with a new discrete spectrum of the time dependent Schrödinger operator corresponding to meromorphic eigenfunctions with poles of higher multiplicity and to what we term nonstandard pole divisors. (A comprehensive account of the spectrum of both KPI and KPII is given in [19].) The extension of these ideas and solutions to DSII equation via spectral analysis of the Dirac operator on the plane has been considered in [20]. We also note that KPI possesses, in addition, other localized, nondecaying solutions like line solitons [21, 22]. The solution of the Cauchy problem in such a background is considered in [14, 23].

In this paper we show that the spectral operator that linearizes (3) possesses a rich discrete spectrum corresponding to smooth, rationally “decaying” lump configurations whose dynamical evolution and scattering are nontrivial. The physical behavior of these solutions is reminiscent to the aforementioned discrete spectrum and nonstandard lumps for KPI. We find that they are associated with higher-order pole meromorphic eigenfunctions of a similar discrete spectrum. Characterization of the discrete spectrum of this operator involves giving the pole multiplicity and an adequate pole divisor, which could be, in the spirit of the ideas of [18, 20], associated with integer winding numbers.

Direct methods to study this class of solutions have also been developed; see [24, 25] for KPI and [26] for DSII. See also [27] for related ideas. These direct algebraic methods, although they do not shed any light on the associated spectral problem, are quite powerful to undertake a general classification of the relevant class of solutions. In this regard note that in [2831] Painleve’s test for (3) was first considered and some special solutions, like line solitons, lumps, and dromion solutions, were found.

The organization of the paper is as follows. A convenient form of the linearizing Lax pair associated with (3) with boundary conditions (2) is introduced in Section 2, following the ideas of [32]. In Section 3 we study the discrete spectrum of meromorphic eigenfunctions. Different relations between the Laurent coefficients (LC) compatible with the Lax pair are established. In Section 4 we study several classes of lump solutions and their dynamics and physical properties. Particular interest is addressed to the problem of determining the number of solitons that ensue from given analytic structure and to determine the motion of these entities. We find that typically lump’s dynamics is a superposition of the “center of mass” motion, which proceeds in a uniform way and individual, lump depending, motion that behaves as with . Interestingly for all cases analyzed the dynamics corresponds to a central configuration of a certain -body problem (a central configuration of the Euler-Lagrange -body problem corresponds to interacting point masses wherein for each body the acceleration vector is directed towards the center of mass and proportional to the distance to the center of mass) where they fall to the mass center. Particularly interesting is the case, similar to the “homothetic Lagrange” solution in the tree body problem, when pulses are located at the vertexes of an equilateral triangle which collapses to the center along a straight line as .

2. Linear Problem

As we have already pointed out system (3) arises as the compatibility of a pair of operators; see [1, 10, 33, 34] for genial ideas. Here we restrict our study to the class of potentials that are nonsingular and satisfy (2). Under these boundary conditions a convenient form of the Lax pair, depending on a complex spectral parameter , is given by the pair of operators:

We assume the existence of an eigenfunction normalized to as . Then, (8) implies that has the asymptotic expansion (with asymptotic coefficients ):

This expression permits recovering the physical amplitude . To determine the phase of the potential consideration of the next order in the expansion is required. We findThus, the potential is determined from

3. Discrete Spectrum and Relationships among Laurent Coefficients

Basic Assumption. The basic assumption of this work is the existence of a solution to (8) and (9), which has a finite or denumerable number of singularities . We suppose that all singularities are higher-order poles and that is analytic away from the singularities. Thus, around any pole , must have a local Laurent expansion:whereare, respectively, the regular and singular parts of the eigenfunction at the pole and , are the Laurent coefficients (LC). Notice that is also termed the singularity principal part or the pole divisor.

Unlike what happens in the regular case, where the eigenfunction is fixed by the corresponding -problem, when singularities exist the inverse problem does not fix uniquely the singular part. Additional information relating different coefficients of the poles divisor is required. In this section we consider examples of meromorphic eigenfunctions related to nonsingular potentials and determine different relationships between coefficients of the poles divisor. We posit the existence of a linear relationship between coefficients of the Laurent expansion in the formwhere is the pole multiplicity and are certain functions to be determined. We say that the integer is the index.

3.1. Simple Poles

We shall first suppose that is an eigenfunction of (8), (9) with a meromorphic dependence on and assume that is a simple pole of . Let be the residue of at . Equation (14) reads

Indeed, letting in (8), (9) we obtain, at dominant orders, that and must satisfy the system of equations:

3.1.1. Simple Poles of Index One

We first consider the simplest linear relationship when Then, from the previous equations must satisfyWe can satisfy these equations by requiring to solve the following system of constant differential equations:

It follows thatTherefore, eigenfunctions of (8) with the structure (16) can be constructed by requiring that solves (17) and that the relationshipholds for some complex constant . A similar relation was first established for KPI in [7]. Whenever such a situation holds we say that the pole divisor is standard. We next study the possibility of having nonstandard pole divisors.

3.1.2. Simple Poles of Index 2

We now show that the analytic structure (16) for eigenfunctions of (8) does not necessarily imply that (24) holds and different linear relationship between the LCs may hold; concretely we consider here the case when are related asIn this case we also need to supplement (17), (18) with higher-order equations:

By insertion of (25) into (26) we obtainThis equation will be satisfied if is again given by (23) and solves By integrationwhere and are complex constants. Equation (30) gives yet another relationship between the three first LCs , , and of simple pole eigenfunctions compatible with the Lax pair.

3.2. Poles of Order Two

We now show how the above ideas can be extended to cover the case when the eigenfunction has a pole of multiplicity at some point : that is, in (14) and , where, to avoid an awkward notation, we write , and also . By letting in (8) we find that the main coefficients must satisfy , and also

Again, the assumed analytic structure does not fix uniquely the way LCs are related. We next consider several possibilities.

3.2.1. Poles of Order Two and Index Two

We first assume that the span of all linear combinations of , and with coefficients dependent on , , is generated by just one of them; concretely we requirewhere , depend on , , . As previously mentioned, from the -equations we find to be given by (23), while the remaining equations yield that must satisfyor, alternatively, , . It follows that (32) is compatible with the given analytic structure provided is given by (23) and for some complex constants ,

3.2.2. Poles of Order Two and Index Three

A different possibility obtains by requiring now In this case, letting in (8) we obtain that

Hence, in view of (36) Requiring independence of , we obtain that , must solveHence, we obtain that and where is a complex constant and an arbitrary function of . The -equation gives Thus, (36) is satisfied provided is given by (23) and

4. Determination of Classes of Potentials

In this section we determine several classes solutions of (1) by considering proper election of meromorphic eigenfunctions and relationships between LCs of the form (15).

4.1. Potentials Corresponding to Simple Poles of Index One

Suppose first that is a meromorphic eigenfunction that has the representation (14) with a finite even number of simple poles. We make the following assumptions (denoted as condition C1):(1)Poles come in pairs , with corresponding residues , :(2)At every pole , the pole divisor is standard and hence satisfies (24):(3)

Suppose that has just two poles , . In this case from the previous assumptions and (24) we must haveSolving this system of linear equations and using (10) we find the potentialwhere is the determinant of the associated matrix:

The field given by (46) will be called the basic lump solution of (3). Obviously, if , where , the solution is regular on the entire plane and decreases rationally.

Consider a new Galilean frame moving with velocities :

One finds also convenient to go to a frame of skew coordinates , where , and are defined in (48) and is a shear transformation along the -axis, where represents a slanting coefficient. In this frame the potential readsThus, in this inertial frame the solution is at rest for all time (objects moving uniformly with velocities in the unprimed frame remain at rest in the transformed primed frame) and strongly localized. However, the maxima structure is richer than what might have been expected. Inspection shows that critical points solve , . With we find that critical points are andwhile for one hasThese points are all candidates to maxima and minima; the number of them will vary depending on the values of the parameters. The parameter space is a two-dimensional plane deprived of the straight line . To describe the situation in the general case we restrict, with no loss of generality, to the first quadrant on the parameter space. The situation varies according to which of the regions do parameters belong. If is in the interior of then two of the four points described above are maxima and the other two are minima. Point is a saddle point between the formers at which .

Particular cases are (previously studied) and . In this case It has two symmetric maxima located at at which while minima are to be found at the mirror images points at which .

When parameters satisfy we obtain a configuration having a maximum at which ; the minimum lump’s amplitude is to be found at the points . At the other boundary point the situation is opposite and is rather reminiscent of dark solitons; it corresponds to a configuration having minimum amplitude at the origin and maximum amplitude at the points . Thus, generically the configuration is multipeaked, localized on an entire region containing several maxima.

4.2. Potentials Corresponding to Index Two

Here we investigate potentials associated with meromorphic eigenfunctions satisfying the following conditions (C2):(1)Poles come in pairs , , with :(2)For the poles are simple with index 2 (i.e., they satisfy (30)). Besides , , are double poles with index 2 in the sense of (32), (34) with constants satisfying , , .

We suppose that has just two poles , with the above properties. In this case one can prove thatwhere is the following definite-positive polynomial:

Note how it depends on complex parameters , , .

We next study the dynamical properties of the field of (56), (57). Clearly is a rational function, regular on the entire plane which decreases as as . Again we find it convenient to use coordinates defined by (48) (see also (49)) in terms of which the “”-function reads (we drop some irrelevant parameters)Note also thatwhere , , , .

Inspection of the solution shows that it corresponds to a coherent two pulsed structures. Each of the pulses has a rich internal structure with several maxima and minima in a similar disposition to that corresponding to the standard lump.

The dynamical behavior of this configuration is different and far more subtle than previously found. It turns out that lumps are not at rest in any inertial frame.

As indicated we consider a system of skew coordinates , defined in (48). We claim that the solutions togive the pulse’s positions as . Indeed, Note that solves (60) if does; thus, the second pulse is the mirror image of the first respect to the origin, and it suffices to describe the motion of the first of them.

We first consider the case for which the dynamics and scattering process are easier to understand. In this case the moving frame (48) and solution are given by

We then find from (60) that asymptotically the first of the pulses has coordinates , where , , and

Thus, in the moving reference frame the trajectories as are straight lines with different slopes. Notice that initially () pulses are located in the second and fourth quadrants and will move to the first and third ones as . For moderate times they collide head-on henceforth undergoing a scattering process. In Figures 1 and 2 we show the collision path and form of the pulses in solution (57) before and after scattering. The scattering angle is easily found to be given byNote that attains a maximum value if (formally) and decreases towards as . Thus, transparent scattering is obtained when while perpendicular scattering corresponds to lumps with .

In the rest at frame the motion is the composition of a uniform motion and the slower one given by (63). Indeed when both coordinate frames are related by a pure Galilean transformation and objects at rest in the primed frame pick a uniform motion with velocities in the original frame. Hence, in the stationary frame the path of the second pulse is no longer the mirror image of the first.

4.3. Potentials Corresponding to Double Poles with Index Three

We now study potentials that correspond to eigenfunctions with double poles with index three. Concretely, assume the following:(1)Poles come in pairs , , with . Hence,(2)The multiplicity of every pole is two and the index three; namely, all pole divisors satisfy (36).

We consider the simplest such potential when .

4.3.1. Pulses of Index Three

We consider in detail the physical properties and dynamics in the case . The potential is given by formula (46) where the tau function isHere and , are complex constants, , are defined in (49), and The position and dynamics of the associated pulses are determined mainly by the constant (alternatively by ). Concretely, we have the following.

Proposition 1. () If the constant the solution is a multipeaked traveling wave of solitonic nature (i.e., stationary in the frame moving with the soliton).
() When the solution is nonstationary. There are exactly three traveling pulses describing “homothetic motion.” Namely, the motion in the frame at rest is a superposition of a “center of mass” uniform motion and an individual and slower motion proportional to .
() In the “center of mass frame” pulses are located at the vertex of a triangle with center of gravity at the origin (see Figure 3 where this situation is shown). The time evolution contracts but does not deform the whole structure. As a result, pulses collapse onto the origin whereupon the configuration regains shape. No deflection angle appears.

To prove the claims we note that the solution reads neater with a skew transformation and a dilatation of coordinates defined by (48) and , . Then dropping some irrelevant constants we have (see (42))

We simply consider here the case when , namely, when the parameters are chosen to satisfyThe -dependence of drops out and hence the entire configuration moves, with respect to the frame at rest, with constant velocity given by . A plot of shows that decreases to zero away from three regions. Thus, the configuration behaves like a multipeaked solution of solitonic nature.

Appendix

Mittag-Leffler Expansions of Eigenfunctions

Recall that we assume that higher-order poles are the only singularities of ; in particular, we exclude the possibility of having essential singularities and/or accumulation points of pole singularities in the finite plane. We study here the possibility of certain Mittag-Leffler expansions of eigenfunctions (see [35] for general ideas).

Theorem A.1. Let be a meromorphic eigenfunction having poles at , with canonical normalization. Suppose that(1)the distance between eigenvalues is bounded below: ,(2) is uniformly bounded away from eigenvalues.Then, for , can be represented as where is the multiplicity of and the pole divisor at .
Reciprocally, if is given by the formal series (A.1), condition is satisfied andThen is a meromorphic eigenfunction satisfying both condition and also along any curve such that . Further,

Proof. The proposition is trivial if the number of pole singularities is finite so we shall assume an infinite number. Note then that implies that they cannot accumulate anywhere on the finite plane. Thus, with no loss of generality we can assume that singularities pile up at infinity but not too quickly. We can order the eigenvalues so that . Next, let set be the set obtained by choosing one eigenvalue from all those that have the same modulus. Let , . It follows that the sequence is strictly increasing: and satisfies , . Further, implies that for all there exists such that for all such that . It follows that there exists such that if , the unit circle of radius ; then uniformly in and : . For fixed the function has poles with multiplicity at and a simple pole at with residue . At the former an easy calculation shows that Res . Thus, Cauchy’s theorem yields thatwhere the sum extends to all eigenvalues satisfying .
Given such a fixed let be the minimum distance between and all : (note that ). We can bound the integrand, uniformly in as follows: where we recall that , say. Further, , . By the theorem of dominated convergence we find that for fixed Sending the result follows.
Reciprocally, if is given by the formal series (A.1) and condition () is satisfied (with , say), we see that exists and, by formally commuting the limit with the sum, find that tends to at infinity. The interchange of limit and sum is validated by noting that if Dist, To prove that the potential is bounded, we use (10):where we have taken the limit under the sum. To validate this note that if Dist, the modulus-maximum principle yields the estimate uniform on : Thus, if the conditions hold we are guaranteed that has an asymptotic expansion , where is given by (A.8). Operating with and since (A.2) validates differentiation under the sum we arrive to Formula (A.3).

The above result can be sharpened when eigenvalues are disposed along a curve.

Proposition A.2. Suppose that the set is contained in a given line for some and satisfies condition . Suppose that the formal series (A.1) satisfies all conditions (A.2) but the third; that is, we drop the requirement . For any let be the cone of complex numbers whose phase satisfies . Then, has, along any curve contained in , an asymptotic expansion , where is given by (A.8). Hence, it is associated with a nonsingular potential given again by (A.3).

Proof. We have the uniform bound:where we use first the modulus-maximum principle for analytic functions and then, with no loss of generality, suppose . For given the sup is attained at the normal projection of onto the ray (note that this estimate can be sharpened but with no further gain). It implies that, along any curve contained in ,

Proposition A.3. Suppose that is given by the formal series (A.1) and conditions () and () are satisfied. Suppose that the set is contained in a given line for some . With no loss of generality order eigenvalues are such that . Further suppose that . Then the formal series (A.1) is convergent to having canonical normalization and is associated with a nonsingular potential .

A proof of this assertion follows by noting that (24) imply that for long , can be expanded asFor both conditions of (A.2) are satisfied. This implies the claim.

Competing Interests

The authors declare herewith that there is no conflict of interests regarding the publication of this article. Our professional judgment concerning the validity of research is not influenced by a secondary interest such as financial gain.

Acknowledgments

The authors acknowledge support from the Spanish Ministerio de Economía y Competitividad under Contracts MTM2012-38445, MAT2013-46308, and Junta de Castilla y Leon SA226U13.