Single and Multi-Soliton Solutions for a Spectrally Deformed Set of Maxwell-Bloch Equations

The spectral deformation technique is a tool of applied mathematics, which enables one to generate many new integrable partial differential equations (pde’s) from a given integrable pde. In its main lines, it can be summarized as follows: Each integrable pde has an associated linear overdetermined system (i.e, a zero curvature representation (ZCR) or a Lax pair) which contains a spectral parameter, λ. According to the established methods of soliton theory, this parameter must be a constant number. It was first suggested in [1] that under some restrictive assumptions, λ can be choosen as a function of the independent parameters of the problem, rather than a constant. The main restriction on λ is that it should be the solution of an overdetermined, nonlinear set of pde’s, which was derived in [2]. General solutions for this set of nonlinear pde’s for λ are presently not known. However, many different particular solutions can be found by inspection. Each of these particular solutions generates a new, integrable pde from the given integrable pde. These new pde’s are called the “spectral deformations”, or simply “deformations” of the original pde. For some applications of the spectral deformation technique in nonlinear optics, see [3–5]. In 1994, the spectral deformation technique was applied to the Maxwell-Bloch (MB) equations of nonlinear optics, which model resonant light–matter interactions [2]. The result was three different deformations, corresponding to three different particular solutions for λ. Two of these deformations were subsequently solved by employing the dressing and Darboux transformations [2,6]. The third one is not solved yet. Our aim in this paper is to investigate this third equation, which we will call “the deformed Maxwell–Bloch equation” (DMB). In Section 1, we will summarize the deformation process through which one can generate DMB from MB. Section 2 will develop a Darboux transformation methodology and construct single soliton solutions for DMB. The behaviour of these solitons will be investigated. Section 3 will construct two soliton solutions and investigate soliton collisions.


Introduction
The spectral deformation technique is a tool of applied mathematics, which enables one to generate many new integrable partial differential equations (pde's) from a given integrable pde.In its main lines, it can be summarized as follows: Each integrable pde has an associated linear overdetermined system (i.e, a zero curvature representation (ZCR) or a Lax pair) which contains a spectral parameter, λ.According to the established methods of soliton theory, this parameter must be a constant number.It was first suggested in [1] that under some restrictive assumptions, λ can be choosen as a function of the independent parameters of the problem, rather than a constant.The main restriction on λ is that it should be the solution of an overdetermined, nonlinear set of pde's, which was derived in [2].General solutions for this set of nonlinear pde's for λ are presently not known.However, many different particular solutions can be found by inspection.Each of these particular solutions generates a new, integrable pde from the given integrable pde.These new pde's are called the "spectral deformations", or simply "deformations" of the original pde.For some applications of the spectral deformation technique in nonlinear optics, see [3][4][5].
In 1994, the spectral deformation technique was applied to the Maxwell-Bloch (MB) equations of nonlinear optics, which model resonant light-matter interactions [2].The result was three different deformations, corresponding to three different particular solutions for λ.Two of these deformations were subsequently solved by employing the dressing and Darboux transformations [2,6].The third one is not solved yet.Our aim in this paper is to investigate this third equation, which we will call "the deformed Maxwell-Bloch equation" (DMB).
In Section 1, we will summarize the deformation process through which one can generate DMB from MB. Section 2 will develop a Darboux transformation methodology and construct single soliton solutions for DMB.The behaviour of these solitons will be investigated.Section 3 will construct two soliton solutions and investigate soliton collisions.

Maxwell-Bloch Equations and Spectral Deformation
Maxwell-Bloch (MB) equations can be expressed as where is a constant dependent on the physical properties of the system, with n 0 being the number of atoms per unit volume, p 0 is the polarization and ω 0 is the frequency of the electric field.E denotes the complex envelope of the electric field; ρ is the amount of polarization in the resonant medium, and N is the amount of population inversion between the two energy levels.models frequency shift from the resonance.It is a trivial field as it can be eliminated from the MB equations by rescaling.However, we will not rescale it, as it will become non-trivial when the system is deformed.
It is well-known that MB equations are integrable and possess the following zero curvature representation (ZCR): where and λ is a complex constant.The MB Equations ( 1)-( 4) can be recovered from the ZCR ( 5) and ( 6) by equating the mixed derivatives, i.e by Ψ ζη = Ψ ηζ , which results in where [.,.] denotes the commutator.If we substitute the expressions for U, V and J as given in (7) into the above equations, (8) will yield (1), and ( 4) and (9) will yield ( 3) and ( 2).This shows the equivalence of ZCR ( 5) and ( 6) to the MB equations.Note that Equation (3) differs from the original version derived in [2], Equation ( 9).These two sets of equations are equivalent when the field is real.However, if the field is taken to be complex, they cease to be equivalent.Furthermore, the ZCR given in [2] becomes invalid, which indicates that with a complex field, Equation ( 9) in [2] loses integrability, while (3) continues to be integrable.
Allowing to become complex is important, as the Darboux transformations forces it to become complex during the first and subsequent iterations (see Equation ( 15)).This implies that Darboux transformation cannot be applied directly to Equation ( 9) in [2] and their ZCR, as these cannot handle a complex field.Hence, the modifications resulting in Equation ( 11) is crucial.Furthermore, note that in Equation (11) while is allowed to be complex, only its real part affects the dynamics of the system.Spectral deformation of the MB equations is obtained when λ is considered not as a constant but as a function of the independent variables ζ and η.We are interested in the following specific deformation: where b is a real constant and k is a complex integration constant which is called the 'hidden spectral parameter' in the spectral deformation literature.Enforcing the consistency conditions or, when ( 7) is substituted, These equations will be called the deformed Maxwell-Bloch equations (DMB) in the rest of this article.The physical meaning of Equation ( 11) is not clear, but we consider it interesting anyway to study them from a mathematical point of view.

Single Soliton Solution for DMB
There are many methods to compute soliton solutions of integrable equations and these methods usually work well with the spectrally deformed versions of the same equations.Recent examples that specifically relate to optical solitons are [7][8][9][10][11], which use a plethora of closely related techniques like Darboux transformations, Dressing transformations or Backlund transformations.This article will use Darboux transformations to construct solitons of DMB.
For the seed solution, choose by inspection, where N 0 , 0 are real constants.Substituting these into (7), and solving for Ψ in ( 5)-( 6) with λ as defined in (10) gives The Darboux iteration maps a solution S 0 = {U, V, Ψ(λ)} of the ZCR ( 5) and ( 6) to another solution S 1 = { Ū, V, Ψ(λ)} of the same ZCR.In practice, it is enough to specify a transformation that maps Ψ(λ) to Ψ(λ).Two other transformations which map {U, V} to { Ū, V} follow as an extension to this transformation.The following ansatz is proposed for mapping Ψ(λ) to Ψ(λ): where Q(x, t) is a 2 × 2 matrix whose exact form is to be determined.This ansatz differs from the standard one used commonly in Darboux transformation theory, given as which only generates trivial solutions for DMB.By contrast, the newly proposed transformation (12) generates nontrivial solutions.To specify Q(ζ, η) it is sufficient to fix a zero of Ψ(λ), i.e, choose λ 1 such that Ψ(λ 1 ) = 0 for all ζ and η.Here λ 1 denotes λ evaluated at a hidden spectral parameter k = k 1 .
This implies the form of the matrix Q: The subscript 1 above (as in Λ 1 and Ψ 1 ) denotes that k and ζ 0 in Λ and Ψ are fixed at k = k 1 and After substituting (14) into the the zero curvature representation above and using ( 5)-( 6) to convert the result into an algebraic equation, we arrive at the final result: Carrying out the algebra, we find all the fields (which are elements of the matrices U(λ 1 ) and V(λ 1 )) corresponding to a single-soliton solution: where Here (.), (.) denote the real and imaginary parts of a complex number.Note that due to the term 2b (λ 1 )i the field (λ 1 ) becomes complex after the application of the Darboux transformation.
Figure 1 plots the electric field of the single-soliton solution vs. ζ and η.Contrary to the undeformed MB equation, the soliton of the DMB equation is a transient.It pops out of the background, attains a maximum, and then decays back into the background.This can be seen much better in Figure 2, which indicates that the soliton has a "lifetime", i.e., it has a significant magnitude only for −300 < η < 300.Furthermore, the soliton does not travel at a uniform speed.The soliton stops for some time before it grows or after it decays, which can be observed in Figure 3.One can also ask the question if the DMB soliton behaves more like the MB soliton as b → 0. Figures 4 and 5 demonstrate that this is the case.There is a region around η = 0 in which the soliton travels with approximately constant amplitude and speed, like the MB soliton.As b becomes smaller, this region gets wider.

Discussion
A Darboux-transformation based methodology to construct the soliton solutions of the DMB equation is proposed.These solitons are shown to behave quite differently compared to the well-known MB solitons, as they have a finite lifetime, constantly changing shape and nonuniform speed.It is also shown that as the deformation parameter b goes to zero, DMB solitons behave more and more like the MB soliton.In all other respects, DMB solitons have the usual soliton behaviour, like elastic collisions.

Figure 2 .
Figure 2. Magnitude of the electric field peak vs. η.Parameters are the same with Figure 1.

Figure 3 .
Figure 3. Position of the electric field peak vs. η.Parameters are the same with Figure 1.

Figure 4 .Figure 5 .
Figure 4. Magnitude of the electric field peak vs. η for various values of b, from b = 0.001 to b = 0.0241 with steps ∆b = 0.003.All other parameters are the same with Figure 1.