Multicomplex solitons
- DOI
- 10.1080/14029251.2020.1683963How to use a DOI?
- Keywords
- nonlinear wave equations; solitons; quaternions; coquaternions; octonions
- Abstract
We discuss integrable extensions of real nonlinear wave equations with multi-soliton solutions, to their bicomplex, quaternionic, coquaternionic and octonionic versions. In particular, we investigate these variants for the local and nonlocal Korteweg-de Vries equation and elaborate on how multi-soliton solutions with various types of novel qualitative behaviour can be constructed. Corresponding to the different multicomplex units in these extensions, real, hyperbolic or imaginary, the wave equations and their solutions exhibit multiple versions of antilinear or π«π―-symmetries. Utilizing these symmetries forces certain components of the conserved quantities to vanish, so that one may enforce them to be real. We find that symmetrizing the noncommutative equations is equivalent to imposing a π«π―-symmetry for a newly defined imaginary unit from combinations of imaginary and hyperbolic units in the canonical representation.
- Copyright
- Β© 2020 The Authors. Published by Atlantis and Taylor & Francis
- Open Access
- This is an open access article distributed under the CC BY-NC 4.0 license (http://creativecommons.org/licenses/by-nc/4.0/).
1. Introduction
The underlying mathematical structure of quantum mechanics, a Hilbert space over the field of complex numbers, can be generalized and modified in various different ways. One may for instance re-define the inner product of the Hilbert space or alter, typically enlarge, the field over which this space is defined. The first approach has been pursued successfully since around twenty years [8], when it was first realized that the modification of the inner product allows to include non-Hermitian Hamiltonians into the framework of a quantum mechanical theory. When these non-Hermitian Hamiltonians are π«π―-symmetric/quasi-Hermitian [7, 28, 32] they possess real eigenvalues when their eigenfunctions are also π«π―-symmetric or pairs of complex conjugate eigenvalues when the latter is not the case. The reality of the spectrum might only hold in some domain of the coupling constant, but break down at what is usually referred to as an exceptional point when at least two eigenvalues coalesce. Higher order exceptional points may occur for larger degeneracies. In order to unravel the structure of the neighbourhood of these points one can make use of the second possibility of generalizations of standard quantum mechanics and change the type of fields over which the Hilbert space is defined. This view helps to understand the bifurcation structure at these points and has been recently investigated for the analytically continued Gross-Pitaevskii equation with bicomplex interaction terms [16, 18, 21]. In a similar spirit, systems with finite dimensional Hilbert spaces have been formulated over Galois fields [37]. Hyperbolic extensions of the complex Hilbert space have been studied in [38]. The standard SchrΓΆdinger equation was bicomplexified in [3] and further studied in [4β6, 35]. Quaternionic and coquaternionic quantum mechanics and quantum field theory have been studied for a long time, see e.g. [2, 19, 20], mainly motivated by the fact that they may be related to various groups and algebras that play a central role in physics, such as SO(3), the Lorentz group, the Clifford algebra or the conformal group. Recently it was suggested that they [9] provide a unifying framework for complexified classical and quantum mechanics. Octonionic Hilbert spaces have been utilized for instance in the study of quark structures [22].
Drawing on various relations between the quantum mechanical setting and classical integrable nonlinear systems that possess soliton solutions, such as the formal identification of the L operator in a Lax pair as a Hamiltonian, many of the above possibilities can also been explored in the latter context. Most direct are the analogues of the field extensions. Previously we demonstrated [12, 13] that one may consistently extend real classical integrable nonlinear systems to the complex domain by maintaining the reality of the energy. Here we go further and investigate multicomplex versions of these type of nonlinear equations. We demonstrate how these equations can be solved in several multicomplex settings and study some of the properties of the solutions. We explore three different possibilities to construct solutions that are not available in a real setting, i) using multicomplex shifts in a real solutions, ii) exploiting the complex representations by defining a new imaginary unit in terms of multicomplex ones and iii) exploiting the idempotent representation. We take π«π―-symmetry as a guiding principle to select out physically meaningful solutions with real conserved quantities, notably real energies. We clarify the roles played by the different types of π«π―-symmetries. For the noncommutative versions, that is quaternionic, coquaternionic and octonionic, we find that imposing certain π«π―-symmetries corresponds to symmetrizing the noncommutative terms in the nonlinear differential equations.
Our manuscript is organized as follows: In section 2 we discuss the construction of bicomplex multi-solitons for the standard Korteweg de-Vries (KdV) equation and its nonlocal variant. We present two different types of construction schemes leading to solutions with different types of π«π―-symmetries. We demonstrate that the conserved quantities constructed from these solutions, in particular the energy, are real. In section 3, 4 and 5 we discuss solution procedures for noncommutative versions of the KdV equation in quaternionic, coquaternionic and octonionic form, respectively. Our conclusions are stated in section 6.
2. Bicomplex solitons
2.1. Bicomplex numbers and functions
We start by briefly recalling some key properties of bicomplex numbers and functions to settle our notations and conventions. Denoting the field of complex numbers with imaginary unit Δ± as
The canonical basis is spanned by the units β, Δ±, π₯, k, involving the two imaginary units Δ± and π₯ with Δ±2 = π₯2 = β1, so that the representations in equations (2.2) and (2.3) naturally prompt the notion to view these numbers as a doubling of the complex numbers. The real unit β and the hyperbolic unit k = Δ±π₯ square to 1, β2 = k2 = 1. The multiplication of these units is commutative with further products in the Cayley multiplication table being βΔ± = Δ±, β π₯ = π₯, βk = k, Δ±k = βπ₯, π₯k = βΔ±. The idempotent representation (2.5) is an orthogonal decomposition obtained by using the orthogonal idempotents
Arithmetic operations are most elegantly and efficiently carried out in the idempotent representation (2.5). For the composition of two arbitrary numbers na and nb we have
The hyperbolic numbers (or split-complex numbers) π» = {a1β + a4k | a1, a4 β β} are an important special case of πΉ obtained in the absence of the imaginary units Δ± and π₯, or when taking a2 = a3 = 0.
The same arithmetic rules as in (2.9) then apply to bicomplex functions. In what follows we are most interested in functions depending on two real variables x and t of the form f(x, t) = βp(x, t) + Δ±q(x, t) + π₯r(x, t) + ks(x, t) β πΉ involving four real fields p(x, t), q(x, t), r(x, t), s(x, t) β β. Having kept the functional variables real, we also keep our differential real, so that we can differentiate f (x, t) componentwise as βx f(x, t) = ββx p(x, t) + Δ±βxq(x, t) + π₯βxr(x, t) + kβxs(x, t) and similarly for βt f(x, t). For further properties of bicomplex numbers and functions, such as for instance computing norms, see for instance [17, 26, 29, 30].
2.2. π«π―-symmetric bicomplex functions and conserved quantities
As there are two different imaginary units, there are three different types of conjugations for bicomplex numbers, corresponding to conjugating only Δ±, only π₯ or conjugating both Δ± and π₯ simultaneously. This is reflected in different symmetries that leave the Cayley multiplication table invariant. As a consequence we also have three different types of bicomplex π«π―-symmetries, acting as
Decomposing a density function for any conserved quantity as
2.3. The bicomplex Korteweg-de Vries equation
Using the multiplication law (2.9) for bicomplex functions, the KdV equation for a bicomplex field in the canonical form
We recall that we keep here our space and time variables, x and t, to be both real so that also the corresponding derivatives βx and βt are not bicomplexified.
When acting on the component functions the π«π―-symmetries (2.10)β(2.12) are implemented in (2.16) as
We observe that (2.16) allows for a scaling of space by the hyperbolic unit k as x β k x, leading to a new type of KdV-equation with u β h
Next we consider various solutions to these different versions of the bicomplex KdV-equation, discuss how they may be constructed and their key properties.
2.3.1. One-soliton solutions with broken π«π―-symmetry
We start from the well known bright one-soliton solution of the real KdV equation (2.16)
Noting that the complex solution uiΞΈ,Ξ±(x, t) studied in [13], can be expressed as uiΞΈ,Ξ±(x, t) = pa,ΞΈ;Ξ±(x β a/Ξ±, t)+iqa,ΞΈ;Ξ±(x β a/Ξ±, t), we can also expand the bicomplex solution (2.25) in terms of the complex solution as
In figure 2.3.1 we depict the canonical components of this solution at different times. We observe in all of them that the one-soliton solution is split into two separate one-soliton-like components moving parallel to each other with the same speed. The real p-component can be viewed as the sum of two bright solitons and the hyperbolic s-component is the sum of a bright and a dark soliton. This effect is the results of the decomposition of each of the components into a sum of the functions pa,b;Ξ± or qa,b;Ξ±, as defined in (2.26), at different values of a, b, but the same value of Ξ±. Since a and b control the amplitude and distance, whereas regulates the speed, the constituents travel at the same speed. We recall that this type of behaviour of degenerate solitons can neither be created from a real nor a complex two-soliton solution [11, 15]. So this is a novel type of phenomenon for solitons previously not observed.
In general, the solution (2.24) is not π«π―-symmetric with regard to any of the possibilities defined above. It becomes π«π―Δ±π₯-symmetric when Ο = Ο = 0, π«π―Δ±k-symmetric when Ο = Ο = Ο = 0 and π«π―jk-symmetric when Ο = Ο = ΞΈ = 0.
A solution to the new KdV equation (2.23) is constructed as
In figure 2.3.1 we depict the canonical component functions of this solution. We observe that the one-soliton solution is split into two one-soliton-like structures that scatter head-on with each other. The real p-component consists of a head-on scattering of two bright solitons and hyperbolic the s-component is a head-on collision of a bright and a dark soliton. Given that uΟ,ΞΈ,Ο,Ο;Ξ± and hΟ,ΞΈ,Ο,Ο;Ξ±(x, t) differ in the way that one of its constituent functions is time-reversed this is to be expected.
2.3.2. π«π―ij-symmetric one-soliton solution
An interesting solution can be constructed when we start with a complex π«π―Δ±k and a complex π«π―π₯k symmetric solution to assemble the linear decomposition of an overall π«π―Δ±π₯-symmetric solution with different velocities. Taking in the decomposition (2.17) v(x, t) = uΔ±ΞΈ,Ξ±(x, t) and w(x, t) = uπ₯Ο,Ξ²(x, t), we can build the bicomplex KdV-solution in the idempotent representation
The expanded version in the canonical representation becomes in this case
2.3.3. Multi-soliton solutions
The most compact way to express the N-soliton solution for the real KdV equation in the form (2.16) is
We could now take the shifts ΞΌ1, ΞΌ2, . . ., ΞΌn β πΉ and expand (2.33) into its canonical components to obtain the N-soliton solution for the bicomplex equation. Alternatively we may also construct N-solitons in the idempotent basis in analogy to (2.32). We demonstrate here the latter approach for the two-soliton. From (2.33) we observe that the second derivative will not alter the linear bicomplex decomposition and it is therefore useful to introduce the quantity w(x, t) as u = wx. Thus a complex one-soliton solution can be obtained from
Noting that
Using (2.37) to define the two complex quantities
Then by construction uΞΈ1,ΞΈ2,ΞΈ3,ΞΈ4,Ο1,Ο2,Ο3,Ο4;Ξ±1,Ξ±2,Ξ²1,Ξ²2 = (w2)x is a bicomplex two-soliton solution with four speed parameters. In a similar fashion we can proceed to construct N-soliton for N > 2.
2.3.4. Real and hyperbolic conserved quantities
Next we compute the first conserved quantities the mass m, the momentum p and the energy E, see e.g. [12, 13]
Decomposing the relevant densities into the canonical basis, u as in (2.15), u2 as
These values are the same as those found in [13] for the complex solitons. Given that the π«π―-symmetries are all broken this is surprising at first sight. However, considering the representation (2.28) this is easily understood when using the result of [13]. Then m(uΟ,ΞΈ,Ο,Ο;Ξ±) is simply β/2(2Ξ± + 2Ξ±) + π₯/2(2Ξ± β 2Ξ±) = 2Ξ±β. We can argue similarly for the other conserved quantities.
For the π«π―ij-symmetric solution ûθ,Ο;Ξ±,Ξ² we obtain the following hyperbolic values for the conserved quantities
The values become real and coincide with the expressions (2.45)β(2.47) when we sum up the contributions from the real and hyperbolic component or in the degenerate case when we take the limit Ξ² β Ξ±.
2.4. The bicomplex Alice and Bob KdV equation
Various nonlocal versions of nonlinear wave equations that have been overlooked previously have attracted considerable attention recently. In reference to standard scenarios in quantum cryptography some of them are also often referred to as Alice and Bob systems. These variants of the nonlinear SchrΓΆdinger or Hirota equation [1, 10, 27, 33] arise from an alternative choice in the compatibility condition of the two AKNS-equations. For the KdV equation (2.16) they can be constructed [23β25] by choosing u(x, t) = 1/2[a(x, t) + b(x, t)], with the constraint π«π― a(x, t) = a(βx, βt) = b(x, t), thus converting it into an equation that can be decomposed into two equations, the Alice and Bob KdV (ABKdV) equation
In a similar way as the two AKNS-equations can be made compatible by a suitable transformation map, these two equations are converted into each other by a π«π―-transformation, i.e. π«π― (2.51) β‘ (2.52). Evidently the decomposition is not unique and one may also add and subtract a constrained function of a and b or consider different types of maps to relate the equation.
The bicomplex version of the Alice and Bob system (2.51), (2.52) is obtained by taking a, b β πΉ. In the canonical basis we use the conventions u(x, t) = βp(x, t) + Δ±q(x, t) + π₯r(x, t) + ks(x, t),
A real solution to the ABKdV equations (2.51) and (2.52) that sums up to the standard one-soliton solution (2.24) is found asx
The functions bΟ,ΞΈ,Ο,Ο;Ξ±, or equivalently the individual components
We may also proceed as in subsection 2.3.2 and construct a solution in the idempotent representation. Keeping the parameter real, a solution based on the idempotent decomposition is
Once more, the functions bΟ,ΞΈ,Ο,Ο;Ξ± or
3. Quaternionic solitons
3.1. Quaternionic numbers and functions
The quaternions in the canonical basis are defined as the set of elements
The multiplication of the basis {β,Δ±, π₯, k} is noncommutative with β denoting the real unit element, β2 = 1 and Δ±, π₯, k its three imaginary units with Δ±2 = π₯2 = k2 = β1. The remaining multiplication rules are Δ±π₯ = β π₯Δ± = k, π₯k = βkπ₯ = Δ± and kΔ± = βΔ±k = π₯. The multiplication table remains invariant under the symmetries π«π―Δ±π₯, π«π―Δ±k and π«π―π₯k. Using these rules for the basis, two quaternions in the canonical basis na = a1β + a2Δ± + a3 π₯ + a4k β π and nb = b1β + b2Δ± + b3π₯ + b4k β π are multiplied as
There are various representations for quaternions, see e.g. [31], of which the complex form will be especially useful for what follows. With the help of (3.2) one easily verifies that
3.2. The quaternionic Korteweg-de Vries equation
Applying now the multiplication law (3.2) to quaternionic functions, the KdV equation for a quaternionic field of the form u(x, t) = βp(x, t) + Δ±q(x, t) + π₯r(x, t) + ks(x, t) β π can also be viewed as a set of coupled equations for the four real fields p(x, t), q(x, t), r(x, t), s(x, t) β β
Notice that when comparing the bicomplex KdV equation (2.16) and the quaternionic KdV equation (3.5) only the signs of the penultimate terms in all four equations have changed. This means that also (3.5) is invariant under the π«π―Δ±π₯-symmetry. Alternatively, we may consider here the aforementioned symmetry
When eliminating these terms from (3.5) the remaining set of equations is π«π―Δ±π₯k-symmetric, which appears to be a rather strong imposition. However, the equations without these terms emerge quite naturally when keeping in mind that the product of functions in (3.5) is noncommutative so that one should symmetrize products and replace 6uux β 3uux + 3uxu. This process corresponds precisely to imposing the constraints (3.7).
3.3. π«π―Δ±π₯k-symmetric N-soliton solutions
Due to the noncommutative nature of the quaternions it appears difficult at first sight to find solutions to the quaternionic KdV equation. However, using the complex representation (3.4), and imposing the π«π―Δ±π₯k-symmetric, we may resort to our previous analysis on complex solitons. Following [13] and considering the shifted solution (2.24) in the complex space β(ΞΎ) yields the solution
This solution becomes π«π―Δ±π₯k-symmetric when we carry out a shift in x or t to eliminate the real part of the shift. Reading off the functions p(x, t), q(x, t), r(x, t), s(x, t) from (3.9), it is also obvious that the constraints (3.7) are indeed satisfied. Thus the real β-component is a one-solitonic structure similar to the real part of a complex soliton and the remaining component consists of the imaginary parts of a complex soliton with overall different amplitudes. It is clear that the conserved quantities constructed from this solution must be real, which follows by using the same argument as for the imaginary part in the complex case [13] separately for each of the Δ±, π₯, k-components. By considering all functions to be in β(ΞΎ), it is also clear that multi-soliton solutions can be constructed in analogy to the complex case β(Δ±) treated in [13] with a subsequent expansion into canonical components.
Since the quaternionic algebra does not contain any idempotents, a construction similar to the one carried out in subsection 2.3.2 does not seem to be possible for quaternions. However, we can use (2.36) for two complex solutions
A coquaternionic two-soliton solution to (3.5) is then obtained from (3.10) as u(2) = (w2)x.
4. Coquaternionic solitons
4.1. Coquaternionic numbers and functions
The coquaternions or often also referred to as split-quaternions in the canonical basis are defined as the set of elements
The multiplication of the basis {β,Δ±, π₯, k} is noncommutative with a real unit element β, β2 = 1, two hyperbolic unit elements π₯, k, π₯2 = k2 = 1, and one imaginary unit Δ±2 = β1. The remaining multiplication rules are Δ±π₯ = βπ₯Δ± = k, π₯k = βk π₯= βΔ± and kΔ± = βΔ±k = π₯. The multiplication table remains invariant under the symmetries π«π―Δ±π₯, π«π―Δ±k and π«π―π₯k. Using these rules for the basis, two coquaternions in the canonical basis na = a1β + a2Δ±+ a3π₯ + a4k β π and nb = b1β + b2Δ±+ b3π₯ + b4k β π are multiplied as
There are various coquaternionic representations for numbers and functions. Similar as a quaternion one can formally view a coquaternion, n1 β π, as an element in the complex numbers
4.2. The coquaternionic Korteweg-de Vries equation
Applying now the multiplication law (4.2) to coquaternionic functions, the KdV equation for a quaternionic field of the form u(x, t) = βp(x, t) + Δ±q(x, t) + π₯r(x, t) + ks(x, t) β π can also be viewed as a set of coupled equations for the four real fields p(x, t), q(x, t), r(x, t), s(x, t) β β. The symmetric coquaternionic KdV equation then becomes
Notice that the last three equations of the coupled equation in (4.6) are identical to the symmetric quaternionic KdV equation (3.5) with constraints (3.7).
4.3. π«π―Δ±π₯k-symmetric N-soliton solutions
Using the representation (4.3) we proceed as in subsection 3.3 and consider the shifted solution (2.24) in the complex space β(ΞΆ)
5. Octonionic solitons
We finish our discussion with a comment on the construction of octonionic solitons. Octonions or Cayley numbers are extensions of the quaternions with a doubling of the dimensions. In the canonical basis they can be represented as
The multiplication of the units is defined by noting that each of the seven quadruplets (e0, e1, e2, e3), (e0, e1, e4, e5), (e0, e1, e7, e6), (e0, e2, e4, e6), (e0, e2, e5, e7), (e0, e3, e4, e7) and (e0, e3, e6, e5), constitutes a canonical basis for the quaternions in one-to-one correspondence with (β,Δ±,π₯,k). Hence the octonions have one real unit, 7 imaginary units and the multiplication of two octonions is noncommutative. Similarly as for quaternions and coquaternions we can view an octonion na β π as a complex number
In order to obtain a π«π―o-symmetry we require a π«π―e1e2e3e4e5e6e7-symmetry in the canonical basis.
5.1. The octonionic Korteweg-de Vries equation
Taking now an octonionic field to be of the form u(x, t) = p(x, t)e0 + q(x, t)e1 + r(x, t)e2 + s(x, t)e3 + t(x, t)e4 + v(x, t)e5 + w(x, t)e6 + z(x, t)e7 β π the symmetric octonionic KdV equation, in this form of (4.6) becomes a set of eight coupled equations
5.2. π«π―e1e2e3e4e5e6e7-symmetric N-soliton solutions
Using the representation (5.2) we proceed as in subsection 3.3 and consider the shifted solution (2.24) in the complex space β(o)
6. Conclusions
We have shown that the bicomplex, quaternionic, coquaternionic and octonionic versions of the KdV equation admit multi-soliton solutions. Using the standard folklore we assume that the existence of such type of solutions implies certain integrability of these equations, which we did not formally prove. The bicomplex versions, local and nonlocal, display a particularly rich structure with the two types of solutions found to exhibit very different types of qualitative behaviour. Especially interesting is the solution in the idempotent representation that decomposes a N-soliton into a 2N-solitonic structure. Each one-soliton constituent of the N-soliton has two contributions that even involve two independent speed parameters. Unlike as for the real and complex solitons, where the degeneracy poses a nontrivial technical problem [11, 15], here these parameters can be trivially set to be equal.
For all noncommuative versions of the KdV equation, i.e. quaternionic, coquaternionic and octonionic, we found multi-soliton solutions based on complex representation in which the imaginary unit is built from specific combinations of the imaginary and hyperbolic units. Interestingly in all cases we observe that the π«π―-symmetry needed to ensure that the newly defined imaginary unit can also be used as a π«π―-symmetry imposes constraints that are equivalent to the constraints needed to obtain the symmetric KdV equation from the nonsymmetric one.
Naturally it would be interesting to extend the analysis presented here to other types of nonlinear integrable systems. A more challenging extension is to multi-complexify also the variables x and t which then also impacts on the definition of the derivatives with respect to these variables.
Acknowledgments
JC is supported by a City, University of London Research Fellowship.
References
Cite this article
TY - JOUR AU - Julia Cen AU - Andreas Fring PY - 2019 DA - 2019/10/25 TI - Multicomplex solitons JO - Journal of Nonlinear Mathematical Physics SP - 17 EP - 35 VL - 27 IS - 1 SN - 1776-0852 UR - https://doi.org/10.1080/14029251.2020.1683963 DO - 10.1080/14029251.2020.1683963 ID - Cen2019 ER -