Pulse Propagation Models with Bands of Forbidden Frequencies or Forbidden Wavenumbers: a Consequence of Abandoning the Slowly Varying Envelope Approximation and Taking into Account Higher-order Dispersion

We study linear and nonlinear pulse propagation models whose linear dispersion relations present bands of forbidden frequencies or forbidden wavenumbers. These bands are due to the interplay between higher-order dispersion and one of the terms (a second-order derivative with respect to the propagation direction) which appears when we abandon the slowly varying envelope approximation. We show that as a consequence of these forbidden bands, narrow pulses radiate in a novel and peculiar way. We also show that the nonlinear equations studied in this paper have exact soliton-like solutions of different forms, some of them being embedded solitons. The solutions obtained (of the linear as well as the nonlinear equations) are interesting since several arguments suggest that the Cauchy problems for these equations are ill-posed, and therefore the specification of the initial conditions is a delicate issue. It is also shown that some of these equations are related to elliptic curves, thus suggesting that these equations might be related to other fields where these curves appear, such as the theory of modular forms and Weierstrass ℘ functions, or the design of cryptographic protocols.


Introduction
The nonlinear Schrödinger (NLS) equation: plays a central role in the study of light pulses propagating in optical fibers.In this equation z represents the distance along an optical fiber, t is the so-called retarded time, u(z, t) is the envelope of the electric field of a laser beam, ε 2 is a real constant whose value depends on the laser´s frequency, and γ 1 depends on the characteristics of the fiber (and also on the frequency of the light).It is worth mentioning that Equation ( 1) is sometimes referred to as the temporal NLS equation, to distinguish it from the spatial NLS equation, whose physical meaning is different, even though it has the same form as Equation ( 1), but with a spatial variable x (corresponding to a transversal coordinate) instead of the retarded time t.
The effects of including these terms in generalized NLS equations have been thoroughly studied in the literature.Moreover, an additional approximation introduced in the deduction of the NLS equation is the slowly-varying-envelope approximation (SVEA).This approximation is adequate when the complex amplitude of an optical pulse varies slowly in space and time.However, when we deal with optical pulses whose widths only contain a few cycles of the carrier wave, it is necessary to improve the SVEA.Several alternatives have been proposed to improve this approximation.As early as 1985, Christodoulides and Joseph studied an extended NLS equation which contained higher-order dispersion and the additional terms u zz , i(|u| 2 u) τ and (|u| 2 u) ττ , where τ was the standard (laboratory) time [28], and they found an interesting exact soliton solution by an algebraic tour de force.In more recent times completely different models have been proposed to describe ultra-short optical pulses, and the famous Kortweg-de Vries (KdV), modified KdV (mKdV) and sine-Gordon (sG) equations have been found useful for this purpose, as well as a combination of the last two of these equations [29,30].Another two equations which have been useful to describe ultra-short spatiotemporal pulses are a two-dimensional sG equation [31] and a cubic generalized Kadomtsev-Petviashvili (cgKP) equation [32,33].It is worth observing that in all these models (KdV, mKdV, sG, mKdV-sG, 2D-sG and cgKP) the temporal variable which appears in the equations is a delayed time (not the laboratory time), similar to the retarded time that appears in the standard NLS equation.
In the context of spatial solitons the paraxial approximation plays a role analogous to the SVEA in temporal solitons, and the effects of nonparaxiality have also been studied.In particular, the following non-paraxial extensions of the spatial NLS equation have been studied [34][35][36]: The study of Equation (2) addresses a basic question which might be of interest to any reader interested in optical solitons: what is the effect of introducing the non-paraxial term u zz in the standard NLS equation?On the other hand, the study of Equation (3) considers an equally interesting question, albeit a slightly more specialized one: what is the effect of taking into account simultaneously nonparaxiality and higher-order nonlinearities?This question suggests that in the field of temporal solitons it would be interesting to investigate the effect of taking into account simultaneously non-SVEA terms, higher-order dispersion and higher-order nonlinearities such as |u| 4 u.A study of this type was carried out in Ref. [28], where a model with these characteristics was proposed (except that the nonlinearity |u| 4 u was not considered).Complex models of this type, which include many different terms, may provide adequate descriptions for the behavior of very short pulses.However, in such models it might be difficult to appreciate the individual effect of each of the terms which have been taken into account.
In the present communication we are interested in studying approximate linear and nonlinear pulse propagation models which incorporate higher-order dispersion, higher-order nonlinearities and one of the terms that appear when we drop the SVEA (a second-order derivative with respect to the propagation direction).When we abandon the SVEA, the terms u zz , i(|u| 2 u) τ and (|u| 2 u) ττ must be introduced in the NLS equation (τ being the standard laboratory time).However, if we use the retarded time to describe the propagation of the optical pulses, in addition to the term u zz , it is also necessary to introduce a mixed derivative of the form u zt in the resulting equation [37].To study the complete equation containing these two derivatives (u zz and u zt ), the nonlinear terms i(|u| 2 u) t and (|u| 2 u) tt , an additional nonlinearity of the form |u| 4 u, and higher-order dispersive terms (−iu ttt and/or u 4t ) might be interesting and important, but it is not the objective of the present work.In this communication we are only interested in studying simplified models where the terms u zt , i(|u| 2 u) t and (|u| 2 u) tt have been discarded.It is clear that these models do not pretend to provide a quantitative accurate description of a particular real system, but we will see that the study of these simplified models reveals that the simultaneous presence of the non-SVEA term u zz and higher-order dispersive terms such as −iu ttt or u 4t may produce interesting results.In particular, we will show that the combined effects of these terms generate bands of forbidden frequencies, which may be related to the gaps which have been experimentally observed in the spectral profiles of very narrow pulses propagating in photonic crystal fibers [38].Moreover, we will see that just the mathematical structure of the equations considered in this paper is interesting by itself, since we will show that these equations are related to elliptic curves, thus suggesting that a relationship may exist between the study of optical solitons and other fields where these curves play an important role, such as the abstruse theory of modular forms and Weierstrass ℘ functions, or the design of cryptographic protocols.
We will begin this communication by studying the linear equations: iu z + c 0 u zz + c 2 u tt + c 4 u 4t = 0 (5) and afterwards we will focus our attention on the nonlinear equations: The study of the linear Equations ( 4) and (5) will reveal how the interplay between the non-SVE term u zz and the higher-order dispersive terms generates bands of forbidden frequencies or forbidden wavenumbers.Then, we will show that the nonlinear Equations ( 6)-( 8) have exact solitons of different types, some of them being embedded solitons.The structure of this paper is the following.In Section 2 we will show that the linear dispersion relations of Equations ( 4) and ( 5) have highly unusual forms.We shall see that the dispersion relation of Equation ( 4) is an elliptic curve, while that of Equation ( 5) presents bands of forbidden frequencies, or forbidden wavenumbers.As the occurrence of elliptic curves in this context opens the possibility of relating Equation (4)-or its solutions-to other areas where elliptic curves play an important role, in Appendix A we describe in more detail the basic characteristics of these curves, and in Appendix B we briefly discuss the relationship of these curves with modular forms and Fermat´s last theorem.In Section 3 we study the linear Equations ( 4) and (5).We will begin this section by paying attention to the fact that several results indicate that the Cauchy problems for these equations are ill-posed, and consequently the specification of initial conditions is a delicate issue.Then we present different solutions of these equations.These solutions will show that narrow pulses that evolve according to Equations ( 4) and ( 5) emit radiation in a completely novel and unexpected way, never observed in other models of optical pulses.In Section 4 exact soliton solutions of Equations ( 6)-( 8) are presented.It is shown that Equation (6) has two different types of solitons (embedded solitons and solitons with a nonlinear frequency shift), Equation (7) has moving solitons of different heights and different velocities, and the solitons of Equation ( 8) are given by squared hyperbolic secants, thus proving that the interplay between the non-SVEA term u zz and higher-order dispersive and nonlinear terms permits the existence of different types of solitons.Finally Section 5 contains the conclusions of the paper.

Dispersion Relation of Equation (4)
Let us begin by paying attention to the form of the dispersion relation of Equation (4).If we substitute the tentative solution: in Equation ( 4), the following dispersion relation is easily found: and this equation can be rescaled by defining k = c 1/2 0 k and ω = c 1/3 3 ω, thus obtaining (after supressing the tildes): where

3
. This equation describes an elliptic curve (if a 2 = 0 and a 3 = 0), which is an interesting result since elliptic curves have important roles in areas which might seem completely unrelated to the study of light pulses in optical fibers, such as the abstruse theory of modular forms [39], or the development of cryptographic protocols [40].For the sake of simplicity in the following we will also say that Equation (10) itself is an elliptic curve, although, being rigorous, we should say that Equation ( 10) is an elliptic curve up to scaling.In other words, it is an elliptic curve distorted by a non-uniform scaling (i.e., a scale transformation which applies different scale factors in the k and ω axes).Now let us pay more attention to the shape of the curves described by the dispersion relation (10).From Equation (10) it follows that the wavenumber k is given by the following function of the frequency: and the form of this function is quite different from the dispersion relations found in other optical systems.If we choose, for example, the coefficients c 0 = 1/30, c 2 = 1/2 and c 3 = 1/15, the form of the curve k(ω) is shown in Figure 1.
Appl.Sci.2017, 7, 340 4 of 32 (embedded solitons and solitons with a nonlinear frequency shift), Equation (7) has moving solitons of different heights and different velocities, and the solitons of Equation ( 8) are given by squared hyperbolic secants, thus proving that the interplay between the non-SVEA term and higher-order dispersive and nonlinear terms permits the existence of different types of solitons.Finally Section 5 contains the conclusions of the paper.

Dispersion Relation of Equation (4)
Let us begin by paying attention to the form of the dispersion relation of Equation (4).If we substitute the tentative solution: in Equation ( 4), the following dispersion relation is easily found: (10) and this equation can be rescaled by defining / and / , thus obtaining (after supressing the tildes): (11) where / and / .This equation describes an elliptic curve (if 0 and 0), which is an interesting result since elliptic curves have important roles in areas which might seem completely unrelated to the study of light pulses in optical fibers, such as the abstruse theory of modular forms [39], or the development of cryptographic protocols [40].For the sake of simplicity in the following we will also say that Equation (10) itself is an elliptic curve, although, being rigorous, we should say that Equation ( 10) is an elliptic curve up to scaling.In other words, it is an elliptic curve distorted by a non-uniform scaling (i.e., a scale transformation which applies different scale factors in the k and ω axes).Now let us pay more attention to the shape of the curves described by the dispersion relation (10).From Equation (10) it follows that the wavenumber k is given by the following function of the frequency: and the form of this function is quite different from the dispersion relations found in other optical systems.If we choose, for example, the coefficients 1/30, 1/2 and 1/15, the form of the curve is shown in Figure 1.We can see that the curve is symmetrical with respect to a horizontal line, but completely asymmetrical with respect to the k-axis, which is an unusual behavior.However, the behavior of We can see that the curve is symmetrical with respect to a horizontal line, but completely asymmetrical with respect to the k-axis, which is an unusual behavior.However, the behavior of k(ω) can be more bizarre.If we choose the coefficients c 0 = 1/40, c 2 = 1/2 and c 3 = 1/50, k(ω) has the form shown in Figure 2.  We can see that in this case the dispersion relation exhibits two bands of forbidden frequencies, thus implying that small-amplitude linear waves cannot propagate with frequencies in the ranges: where: are the frequencies which make the radicand that appears in Equation ( 12) equal to zero.This radicand becomes negative for frequencies in the intervals shown in (13), and therefore becomes complex for these frequencies.In the next section we will investigate the consequences of the existence of these bands of forbidden frequencies.The fact that the linear dispersion relation of Equation ( 4) is an elliptic curve suggests that this equation might be related to other subjects where these curves also appear.In particular, a relationship may exist between optical solitons and the theory of modular forms, which played a central role in the proof of Fermat's last theorem.For this reason in Appendix A we explain more precisely what an elliptic curve is, and in Appendix B we briefly review what was the role played by elliptic curves and modular forms in the proof of Fermat's last theorem.
To close this sub-section it is worth observing that if we introduced the mixed derivative in Equation ( 4), the resulting dispersion relation would still be an elliptic curve.

Dispersion Relation of Equation (5)
Now let us direct our attention to Equation (5).If we substitute the plane wave (9) into Equation ( 5) we arrive at the following dispersion relation: Let us now study the form of this dispersion relation in more detail.To determine the range of wavenumbers permitted by this relation it is convenient to write it in the form: This equation implies that every real value of k will be permitted (i.e., every real value of k corresponds to a positive ) if and only if: (17) In the opposite case, when: (18) We can see that in this case the dispersion relation exhibits two bands of forbidden frequencies, thus implying that small-amplitude linear waves cannot propagate with frequencies in the ranges: where: are the frequencies which make the radicand that appears in Equation ( 12) equal to zero.This radicand becomes negative for frequencies in the intervals shown in (13), and therefore k(ω) becomes complex for these frequencies.In the next section we will investigate the consequences of the existence of these bands of forbidden frequencies.The fact that the linear dispersion relation of Equation ( 4) is an elliptic curve suggests that this equation might be related to other subjects where these curves also appear.In particular, a relationship may exist between optical solitons and the theory of modular forms, which played a central role in the proof of Fermat's last theorem.For this reason in Appendix A we explain more precisely what an elliptic curve is, and in Appendix B we briefly review what was the role played by elliptic curves and modular forms in the proof of Fermat's last theorem.
To close this sub-section it is worth observing that if we introduced the mixed derivative u zt in Equation ( 4), the resulting dispersion relation would still be an elliptic curve.

Dispersion Relation of Equation (5)
Now let us direct our attention to Equation (5).If we substitute the plane wave (9) into Equation ( 5) we arrive at the following dispersion relation: Let us now study the form of this dispersion relation in more detail.To determine the range of wavenumbers permitted by this relation it is convenient to write it in the form: This equation implies that every real value of k will be permitted (i.e., every real value of k corresponds to a positive ω 2 ) if and only if: In the opposite case, when: the dispersion relation contains a band of forbidden wavenumbers.When the condition ( 18) is satisfied, the analysis of the function defined by Equation (16) shows that ω 2 turns out to be negative if k is in the interval: where: Therefore, when the condition (18) holds, the inequalities ( 19) define a band of forbidden wavenumbers.
In Figure 3 we can see the shape of the dispersion relation k(ω) defined by Equation ( 15) when c 0 = 1/30, c 2 = 1/2 and c 4 = 1/80.As these values satisfy the condition (18), the gap of forbidden wavenumbers seen in the figure is explained.The dashed horizontal lines indicate the wavenumbers C 1 and C 2 of the soliton solutions of Equation ( 6) (with γ 1 = 5 and γ 2 = 1) that will be determined in Section 4.
Appl.Sci.2017, 7, 340 6 of 32 the dispersion relation contains a band of forbidden wavenumbers.When the condition ( 18) is satisfied, the analysis of the function defined by Equation (16) shows that turns out to be negative if k is in the interval: (19) where: Therefore, when the condition (18) holds, the inequalities ( 19) define a band of forbidden wavenumbers.
In Figure 3 we can see the shape of the dispersion relation defined by Equation ( 15) when 1/30, 1/2 and 1/80.As these values satisfy the condition (18), the gap of forbidden wavenumbers seen in the figure is explained.The dashed horizontal lines indicate the wavenumbers and of the soliton solutions of Equation ( 6) (with 5 and 1) that will be determined in Section 4.  5) and ( 6), defined by Equation ( 15), with 1/30, 1/2 and 1/80, which satisfy the condition (18).The dashed lines , indicate the wavenumbers and of the solution of Equation ( 6) given by Equations ( 52)-( 55) with 5 and 1.
Now let us study in more detail what happens when the condition (17) holds.In this case the linear dispersion relation presents two bands of forbidden frequencies, which is a novel situation in the study of optical solitons.To convince ourselves that such forbidden frequency gaps really exist, it is convenient to write the linear dispersion relation (15) in the form: This equation implies that if holds, there are frequencies for which k turns out to be complex.Analyzing the radicand which appears on the right-hand-side of Equation ( 21) we find that these forbidden frequencies are in the intervals: and (22) and the frequencies and are defined by: where the minus sign corresponds to , and the plus sign defines .In Figure 4 we can see the shape of the dispersion relation (21) in the particular case when 1/15 , 1/2 and 1/256.These coefficients satisfy the condition (17), and therefore there are two bands of forbidden frequencies, as we can see in the figure.The consequences of the presence of these forbidden frequency gaps will be investigated in the following section.The dashed horizontal lines  5) and ( 6), defined by Equation (15), with c 0 = 1/30, c 2 = 1/2 and c 4 = 1/80, which satisfy the condition (18).The dashed lines k = C 1,2 indicate the wavenumbers C 1 and C 2 of the solution of Equation (6) given by Equations ( 52)-( 55) with γ 1 = 5 and Now let us study in more detail what happens when the condition (17) holds.In this case the linear dispersion relation presents two bands of forbidden frequencies, which is a novel situation in the study of optical solitons.To convince ourselves that such forbidden frequency gaps really exist, it is convenient to write the linear dispersion relation (15) in the form: This equation implies that if c 0 c 2 2 > c 4 holds, there are frequencies for which k turns out to be complex.Analyzing the radicand which appears on the right-hand-side of Equation ( 21) we find that these forbidden frequencies are in the intervals: −ω 2 < ω < −ω 1 and ω 1 < ω < ω 2 (22) and the frequencies ω 1 and ω 2 are defined by: where the minus sign corresponds to ω 1 , and the plus sign defines ω 2 .In Figure 4 we can see the shape of the dispersion relation (21) in the particular case when c 0 = 1/15, c 2 = 1/2 and c 4 = 1/256.These coefficients satisfy the condition (17), and therefore there are two bands of forbidden frequencies, as we can see in the figure.The consequences of the presence of these forbidden frequency gaps will be investigated in the following section.The dashed horizontal lines shown in this figure indicate the values of the wavenumbers of the soliton solutions of Equation ( 6) (with γ 1 = 5 and γ 2 = 1) that will be found in Section 4.   5) and ( 6), given by Equation ( 21), with 1/15, 1/2 and 1/256, which satisfy the condition (17).The dashed lines , indicate the wavenumbers and of the solution of Equation ( 6) given by Equations ( 52)-( 55) with 5 and 1.
We now know that the dispersion relation (15) has interesting forms: when there is a band of forbidden wavenumbers (as in Figure 3), and when there are two bands of forbidden frequencies (as in Figure 4).The only case that we have not examined is when the coefficients , and are on the surface defined by the equation: (24) Only in this case (intermediate between the inequalities (17) and ( 18)) does the dispersion relation (15) have no gaps.An example of this case is obtained when 1/20, 1/2 and 1/80.With these coefficients the shape of the dispersion relation is shown in Figure 5.We have thus seen that the dispersion relations of Equations ( 4) and ( 5) have interesting and unusual forms.Now we will investigate the behavior of solitary waves that obey these equations.

The Issue of the Initial Conditions
In the study of optical solitons many variants of the nonlinear Shrödinger (NLS) equation have been considered in the past in order to describe the propagation of light pulses in optical fibers under different conditions.However, different as they are, most of these variants share a common denominator: they are first-order equations with respect to the evolution variable z.This characteristic implies that in all these cases the initial condition 0, is sufficient to determine particular  5) and ( 6), given by Equation ( 21), with c 0 = 1/15, c 2 = 1/2 and c 4 = 1/256, which satisfy the condition (17).The dashed lines k = C 1,2 indicate the wavenumbers C 1 and C 2 of the solution of Equation ( 6) given by Equations ( 52)-( 55) with γ 1 = 5 and γ 2 = 1.
We now know that the dispersion relation (15) has interesting forms: when c 0 c 2 2 < c 4 there is a band of forbidden wavenumbers (as in Figure 3), and when c 0 c 2 2 > c 4 there are two bands of forbidden frequencies (as in Figure 4).The only case that we have not examined is when the coefficients c 0 , c 2 and c 4 are on the surface defined by the equation: Only in this case (intermediate between the inequalities ( 17) and ( 18)) does the dispersion relation ( 15) have no gaps.An example of this case is obtained when c 0 = 1/20, c 2 = 1/2 and c 4 = 1/80.With these coefficients the shape of the dispersion relation is shown in Figure 5.   5) and ( 6), given by Equation ( 21), with 1/15, 1/2 and 1/256, which satisfy the condition (17).The dashed lines , indicate the wavenumbers and of the solution of Equation ( 6) given by Equations ( 52)-( 55) with 5 and 1.
We now know that the dispersion relation (15) has interesting forms: when there is a band of forbidden wavenumbers (as in Figure 3), and when there are two bands of forbidden frequencies (as in Figure 4).The only case that we have not examined is when the coefficients , and are on the surface defined by the equation: (24) Only in this case (intermediate between the inequalities ( 17) and ( 18)) does the dispersion relation (15) have no gaps.An example of this case is obtained when 1/20, 1/2 and 1/80.With these coefficients the shape of the dispersion relation is shown in Figure 5.We have thus seen that the dispersion relations of Equations ( 4) and ( 5) have interesting and unusual forms.Now we will investigate the behavior of solitary waves that obey these equations.

The Issue of the Initial Conditions
In the study of optical solitons many variants of the nonlinear Shrödinger (NLS) equation have been considered in the past in order to describe the propagation of light pulses in optical fibers We have thus seen that the dispersion relations of Equations ( 4) and ( 5) have interesting and unusual forms.Now we will investigate the behavior of solitary waves that obey these equations.

The Issue of the Initial Conditions
In the study of optical solitons many variants of the nonlinear Shrödinger (NLS) equation have been considered in the past in order to describe the propagation of light pulses in optical fibers under different conditions.However, different as they are, most of these variants share a common denominator: they are first-order equations with respect to the evolution variable z.This characteristic implies that in all these cases the initial condition u(z = 0, t) is sufficient to determine particular solutions of the corresponding equations.This fact seems to be consistent with physical reality: if we have an optical fiber (with well-defined characteristics: shape, composition, refraction index, type of cladding, etc.), and a laser beam (of known frequency) is sent into the fiber, the light intensity (as a function of time) at the beginning of the fiber (i.e., at z = 0) is all we need to determine completely the behavior of light along the fiber.
On the other hand, the situation with Equations ( 4)-( 8) is completely different.These equations are of second order in z, and this characteristic has two important consequences: With these equations we require two initial conditions, u(z = 0, t) and u z (z = 0, t), in order to determine a particular solution.This requirement is somewhat perplexing since physical intuition suggests that u(z = 0, t) should be enough to determine the evolution of light along the fiber.Therefore, it seems as if we were in the presence of a contradiction: while mathematics indicates that Equations ( 4)-( 8) require two initial conditions to define a particular solution, physics suggests that the knowledge of u(0, t) should be enough to define completely the behavior of the system.

CONSEQUENCE 2:
Although formal proofs are still lacking, several arguments suggest that the Cauchy problems for Equations ( 4)-( 8) might be ill-posed (in certain regions of the space of coefficients c n ).The following results permit us to understand why this ill-posedness is to be expected: (i) The presence of the non-paraxial term u zz in the equation: makes the Cauchy problem for this equation ill-posed [21,41].(ii) The presence of higher-order nonlinearities in equations of the form: makes their Cauchy problems ill-posed if m ≥ 4 and λ > 0 [42].(iii) The presence of higher-order derivative u 4t in the equation: makes the Cauchy problem for this equation linearly ill-posed [43].(iv) The solutions of the linear Equations ( 4) and ( 5) and the linear parts of Equations ( 6)-( 8) contain decaying as well as growing modes of the form exp(±σz − iωt), as we shall see in Sections 3.2 and 3.3.The existence of these modes makes the corresponding Cauchy problems linearly ill-posed [43].(v) Attemps of solving initial value problems for Equations ( 4)-( 8) by finite differences in t, and Runge-Kutta algorithms in z, encounter difficulties that are typical of ill-posed problems.
At first sight these two consequences seem to be unrelated.However, as we explain in the following, there is a connection between them.Let us begin by paying attention to Consequence 2. If the Cauchy problems associated to Equations ( 4)-( 8) are indeed ill-posed, then (as Hadamard first observed [44]) it would be necessary to impose certain compatibility conditions among the initial conditions, for these problems to have global solutions (i.e., solutions defined for all t).Consequently, these compatibility conditions will establish a link between the two initial conditions u(z = 0, t) and u z (z = 0, t), thus implying that, essentially, the initial condition u(z = 0, t) does indeed determine the future of the system, as physical intuition suggested.In this form, the compatibility conditions required by the ill-posedness of these problems (i.e., by Consequence 2), solve the apparent contradiction contained in Consequence 1 (i.e., the contrast between the mathematical requirement of two initial conditions, and the physical suggestion that only one condition should be necessary).
In the following two sub-sections we will examine different solutions of Equations ( 4) and ( 5), corresponding to different initial conditions u(z = 0, t).We will see that the behavior of these solutions is particularly interesting if the initial condition is a narrow solitary wave.

Solutions of Equation (4)
Let us obtain the solution of Equation ( 4) corresponding to a given initial condition u(0, t).We will see that we do not have to specify the value of u z (0, t).
Let us begin by calculating the Fourier transform (FT) of Equation ( 4) (with respect to t): In this equation we have defined: The solution of Equation ( 28) can be obtained immediately: where λ a and λ b are the following functions of ω: where the plus sign in front of the radical corresponds to λ a , the minus sign corresponds to λ b , and the radicand S(ω) is given by: From Equation ( 30) it follows that: and from Equations ( 30) and (33) we can obtain the system: and this system gives us the values of the coefficients a(ω) and b(ω) in terms of the functions U(0, ω) and U z (0, ω): Now we can observe something interesting.Let us consider, for example, the particular case when c 0 = 1/40, c 2 = 1/2 and c 3 = 1/50.In this case the radicand S(ω) (given in (32)) has the form shown in Figure 6.

0,
and this system gives us the values of the coefficients and in terms of the functions 0, and 0, : Now we can observe something interesting.Let us consider, for example, the particular case when 1/40, 1/2 and 1/50.In this case the radicand (given in ( 32)) has the form shown in Figure 6.We can see in this figure that S(ω) > 0 if ω < ω 1 or ω 2 < ω < ω 3 , and S(ω) < 0 if ω 1 < ω < ω 2 or ω > ω 3 , where the frequencies ω n are the values shown in (14).Taking into account the signs of S(ω), we can write λ a and λ b in the following way: We can see that λ a has a positive real part when ω < ω 1 and ω 2 < ω < ω 3 , and therefore for frequencies within these two intervals the first term on the right-hand-side (rhs) of Equation ( 30) grows exponentially in z.As this growth is completely unphysical (as it would imply a fictitious energy growth), the coefficient a(ω) has to be equal to zero, and consequently from Equation (36) it follows that there is an unexpected relationship between U z (0, ω) and U(0, ω): thus implying that: Therefore, as far as Equation ( 4) is concerned, the perplexing situation associated to the Consequence 1 (mentioned in the previous sub-section) is completely clarified: physically realistic solutions of Equation ( 4) are indeed completely determined by the initial condition u(0, t), as physics suggested.Once u(0, t) is chosen, the function u z (0, t) is determined by Equation (41).Now let us calculate the solution of Equation ( 4) corresponding to a soliton-like initial condition of the form: u(0, t) = A sech(Bt) (42) If we take the coefficients c 0 = 1/40, c 2 = 1/2 and c 3 = 1/50 in Equation ( 4), and the parameters A = 1 and B = 1.5 in Equation ( 42), the evolution of the pulse can be seen in Figure 7, which shows the form of |u(z, t)| for z = 0.5, 1.0 and 1.5.We do not see anything special in this figure, which is somewhat disappointing, since we expected an unusual behavior due to the awkward form of the dispersion relation k(ω) shown in Figure 2.
Therefore, as far as Equation ( 4) is concerned, the perplexing situation associated to the Consequence 1 (mentioned in the previous sub-section) is completely clarified: physically realistic solutions of Equation ( 4) are indeed completely determined by the initial condition 0, , as physics suggested.Once 0, is chosen, the function 0, is determined by Equation (41).Now let us calculate the solution of Equation ( 4) corresponding to a soliton-like initial condition of the form: If we take the coefficients 1/40 , 1/2 and 1/50 in Equation ( 4), and the parameters 1 and 1.5 in Equation ( 42), the evolution of the pulse can be seen in Figure 7, which shows the form of | , | for 0.5, 1.0 and 1.5.We do not see anything special in this figure, which is somewhat disappointing, since we expected an unusual behavior due to the awkward form of the dispersion relation shown in Figure 2.  If we now obtain the solution of Equation ( 4) (with the same coefficients indicated above) and using again an initial condition of the form (42), but with A = 1/2 and B = 14, the result is more interesting.In this case small ripples are generated on the leading edge of the pulse (i.e., the left hand side of the pulse), as shown in Figure 8.If we now obtain the solution of Equation ( 4) (with the same coefficients indicated above) and using again an initial condition of the form (42), but with 1/2 and 14, the result is more interesting.In this case small ripples are generated on the leading edge of the pulse (i.e., the left hand side of the pulse), as shown in Figure 8.As the pulse advances along the fiber these ripples form a well-defined radiation wavetrain which propagates to the left of the pulse, as shown in Figures 9 and 10.As the pulse advances along the fiber these ripples form a well-defined radiation wavetrain which propagates to the left of the pulse, as shown in Figures 9 and 10.As the pulse advances along the fiber these ripples form a well-defined radiation wavetrain which propagates to the left of the pulse, as shown in Figures 9 and 10.At first sight, the radiation seen in Figure 9 might seem similar to the small-amplitude radiation waves emitted by solitary pulses which obey equations such as the following [1,2,11]:  As the pulse advances along the fiber these ripples form a well-defined radiation wavetrain which propagates to the left of the pulse, as shown in Figures 9 and 10.At first sight, the radiation seen in Figure 9 might seem similar to the small-amplitude radiation waves emitted by solitary pulses which obey equations such as the following [1,2,11]: At first sight, the radiation seen in Figure 9 might seem similar to the small-amplitude radiation waves emitted by solitary pulses which obey equations such as the following [1,2,11]: However, Figures 9 and 10 show that the radiation emitted by the solution of Equation ( 4) possesses a unique characteristic, that no other model of optical pulses has ever predicted: the emission of radiation stops abruptly at a certain point along the retarded-time axis, and this point moves away from the pulse's peak at a constant velocity.To understand clearly what this means we should remember that the independent variable "t" which appears in Equations ( 1), ( 4)-( 8) and ( 26), and appears at the horizontal axes in Figures 8-10, is not the laboratory time, but the so-called retarded time.Therefore, the graphs shown in Figures 8-10 display the time-dependence of the square of the light intensity as measured by observers placed at different points (i.e., different values of z) along the optical fiber, but each of these observers has shifted his/her time origin in such a way that the time t = 0 corresponds to the instant when the pulse's maximum passes through the point z.Consequently, a graph such as the one shown in Figure 10 tells us that the radiation arrives at the observer placed at the position z before the arrival of the pulse.This implies that the radiation advances faster than the pulse.Moreover, it should be noticed that the point where we see an abrupt interruption of the radiation wave corresponds to the endpoint of the radiation wavetrain, since the first radiation waves that were emitted have already moved far away towards the left of the graph.
If we obtain numerical results similar to those shown in Figures 9 and 10, but for z = 2 and z = 2.5, we can see that if the observer advances a distance ∆z = 0.5 along the optical fiber, the endpoint of the radiation wavetrain recedes to the left of the retarded time axis approximately ∆t ≈ −5, thus implying that this endpoint moves along the time axis with an "inverse velocity" ∆t/∆z ≈ −10.It might be worth emphasizing that in the study of light pulses along optical fibers the evolution variable is not the time, but the distance along the fiber, and therefore the "inverse velocity" ∆t/∆z is more adequate to describe the movement of radiation fronts (or pulses) such as those shown in Figures 9 and 10, when they travel along the retarded time axis.
To understand why the end point of the radiation wave moves to the left of the retarded time axis with an inverse velocity ∆t/∆z ≈ −10 we must remember that the radiation emitted by optical pulses is usually the result of a resonance between these pulses and the small-amplitude radiation waves capable of traveling along the fiber.This resonance occurs when an optical pulse has a wavenumber k p that is contained in the range of wavenumbers permitted by the dispersion relation k(ω) of the system.In the case of Equation ( 4) we know that the dispersion relation is given by Equation ( 12), and the form of k(ω) is shown in Figure 2 for the coefficients c 0 = 1/40, c 2 = 1/2 and c 3 = 1/50.What we do not know is the wavenumber corresponding to the solution of Equation ( 4) whose profiles are shown in Figures 8-10.However, we can determine this wavenumber from the graphs of the real and imaginary parts of u(z, t = 0) shown in Figures 11 and 12.
the position z before the arrival of the pulse.This implies that the radiation advances faster than the pulse.Moreover, it should be noticed that the point where we see an abrupt interruption of the radiation wave corresponds to the endpoint of the radiation wavetrain, since the first radiation waves that were emitted have already moved far away towards the left of the graph.
If we obtain numerical results similar to those shown in Figures 9 and 10, but for 2 and 2.5, we can see that if the observer advances a distance Δ 0.5 along the optical fiber, the endpoint of the radiation wavetrain recedes to the left of the retarded time axis approximately Δ 5, thus implying that this endpoint moves along the time axis with an "inverse velocity" Δ /Δ 10.It might be worth emphasizing that in the study of light pulses along optical fibers the evolution variable is not the time, but the distance along the fiber, and therefore the "inverse velocity" Δ /Δ is more adequate to describe the movement of radiation fronts (or pulses) such as those shown in Figures 9 and 10, when they travel along the retarded time axis.
To understand why the end point of the radiation wave moves to the left of the retarded time axis with an inverse velocity Δ /Δ 10 we must remember that the radiation emitted by optical pulses is usually the result of a resonance between these pulses and the small-amplitude radiation waves capable of traveling along the fiber.This resonance occurs when an optical pulse has a wavenumber that is contained in the range of wavenumbers permitted by the dispersion relation of the system.In the case of Equation ( 4) we know that the dispersion relation is given by Equation ( 12), and the form of is shown in Figure 2 for the coefficients 1/40, 1/2 and 1/50.What we do not know is the wavenumber corresponding to the solution of Equation ( 4) whose profiles are shown in Figures 8-10.However, we can determine this wavenumber from the graphs of the real and imaginary parts of , 0 shown in Figures 11 and 12.   From these figures it follows that k p = −40, and therefore a horizontal line placed at k = −40 intersects the dispersion relation k(ω) shown in Figure 2 at two points: (k, ω) = (−40, 0) and (k, ω) = (−40, 25).The first of these points does not correspond to a traveling wave (as its frequency is zero), but the second one corresponds to a linear wave with a resonant frequency ω = 25, phase velocity k/ω = −1.6 and group velocity dk/dω = −12.5 (i.e., in fact we are talking of inverse velocities, ∆t/∆z, as mentioned in the previous paragraph).We can see, therefore, that the endpoint of the radiation wavetrain that we see in Figures 9 and 10 moves to the left of the retarded time axis at the group velocity corresponding to the resonant frequency ω = 25.
We now understand the origin of the velocity of the endpoint of the radiation wavetrain shown in Figures 9 and 10.However, we still do not know why this endpoint exists.In other words, we still do not know why the pulse stops emitting resonant radiation at a certain instant.
The clue to understand why the emission of radiation ends abruptly lies in the contrast between Figures 7 and 10. Figure 7 shows no radiation at all, while in Figure 10 the radiation is quite conspicuous.The reason for this difference must be necessarily in the initial conditions used to obtain the solutions of Equation ( 4) which led to these two figures.These initial conditions are shown in Figure 13, and the absolute values of their Fourier transforms are shown in Figures 14 and 15.We can see that the initial condition which lead to Figure 10 is much narrower than the initial condition corresponding to Figure 7, but the graphs which really suggest the answer to the presence of radiation in Figure 10, and its absence in Figure 7, are the Fourier transforms (FT) shown in Figures 14 and 15.
The FT shown in Figure 14 shows that the spectrum corresponding to the initial condition which lead to Figure 7 is essentially contained in the interval −5 < ω < 5, and it is practically zero for 14, which is the second interval where the dispersion relation k(ω) is well defined (remember, from Figure 2, that k(ω) is defined in two intervals: ω 1 < ω < ω 2 and ω > ω 3 ) Therefore, it is impossible for the initial condition which lead to Figure 7 to resonate with the linear waves with frequencies ω > ω 3 , since this initial condition has no frequency components in this interval.On the contrary, the spectrum shown in Figure 15 shows that the initial condition which leads to Figure 10 contains significant frequency components with ω > ω 3 .Consequently, when this initial pulse starts traveling along the fiber, it resonates with the linear wave with wavenumber k p = −40 (the pulse's wavenumber) and resonant frequency ω = 25, and thus the emission of radiation begins.However, the pulse's width will begin to grow (and its spectrum will become narrower) as a result of the dispersive terms c 2 u tt and ic 3 u ttt which appear in Equation ( 4).As the dispersion of the pulse will continue indefinitely (since Equation (4) does not contain any nonlinear term which could cancel this dispersion), a moment will necessarily arrive when the pulse will be too wide, and its spectrum will be too narrow to have frequency components with ω > ω 3 .Moreover, the components of the pulse with frequencies in the interval ω 2 < ω < ω 3 will die quickly, since the dispersion relation (shown in Figure 2) does not permit the propagation of linear waves with these frequencies.Therefore, a moment will arrive when the pulse´s spectrum will be essentially contained in the interval ω 1 < ω < ω 2 (where the dispersion relation allows the propagation of linear waves), and then the pulse will no longer be able to resonate with the wave with ω = 25 (the resonant frequency), and the emission of radiation will stop.So, the interruption of the radiation is due to the fact that the pulse's spectrum is confined to live in the interval ω 1 < ω < ω 2 (far from the resonant frequency ω = 25), and the presence of the forbidden frequency gap ω 2 < ω < ω 3 seen in Figure 2 strongly contributes to this confinement.
essentially contained in the interval (where the dispersion relation allows the propagation of linear waves), and then the pulse will no longer be able to resonate with the wave with 25 (the resonant frequency), and the emission of radiation will stop.So, the interruption of the radiation is due to the fact that the pulse's spectrum is confined to live in the interval (far from the resonant frequency 25), and the presence of the forbidden frequency gap seen in Figure 2 strongly contributes to this confinement.waves with these frequencies.Therefore, a moment will arrive when the pulse´s spectrum will be essentially contained in the interval (where the dispersion relation allows the propagation of linear waves), and then the pulse will no longer be able to resonate with the wave with 25 (the resonant frequency), and the emission of radiation will stop.So, the interruption of the radiation is due to the fact that the pulse's spectrum is confined to live in the interval (far from the resonant frequency 25), and the presence of the forbidden frequency gap seen in Figure 2 strongly contributes to this confinement.quickly, since the dispersion relation (shown in Figure 2) does not permit the propagation of linear waves with these frequencies.Therefore, a moment will arrive when the pulse´s spectrum will be essentially contained in the interval (where the dispersion relation allows the propagation of linear waves), and then the pulse will no longer be able to resonate with the wave with 25 (the resonant frequency), and the emission of radiation will stop.So, the interruption of the radiation is due to the fact that the pulse's spectrum is confined to live in the interval (far from the resonant frequency 25), and the presence of the forbidden frequency gap seen in Figure 2 strongly contributes to this confinement.We have thus seen that a very original consequence of having a dispersion relation given by an elliptic curve such as that shown on Figure 2 is the abrupt interruption of the resonant radiation emitted by a narrow pulse which obeys the linear Equation (4).Therefore, this novel radiation process is the consequence of the simultaneous presence of the non-slowly-varying-envelope term u zz and the third-order temporal derivative.
To close this sub-section it is worth observing that no qualitative changes are to be expected if we introduced the mixed derivative u zt in Equation ( 4).The equation thus obtained could be solved by Fourier transforms, just as we solved Equation (4).The form of the Equation (30) would still be valid, although the form of the functions λ a,b (ω) and S(ω) would be slightly different.In this case (when we take into account a term of the form cu zt in Equation ( 4)), instead of Equations ( 31) and ( 32) we would have: It should be emphasized that the key ingredient which produces the peculiar radiation process shown in Figures 9 and 10 is the presence of gaps of forbidden frequencies in the dispersion relation of Equation (4).It is interesting to investigate if these gaps survive if we introduce a term of the form cu zt in Equation ( 4).Therefore, let us consider the equation: We can see that the finite gap of forbidden positive frequencies survives in Figures 16 and 17, but it disappears in Figure 18.Therefore, the presence of the term cu zt in Equation ( 47) modifies the gaps of forbidden frequencies in a significant way.In fact, the gap of forbidden positive frequencies will disappear if the coefficients c, c 0 , c 2 and c 3 are such that the equation S(ω) = 0 has only one real root (S(ω) being the polynomial shown in Equation ( 46)).In particular, for pulses propagating in a typical optical fiber with a group velocity v g ∼ = 0.68c L (c L being the light speed in vacuum) the band of forbidden frequencies may not exist for processes whose characteristic dimensions are Z ∼ = 1 mm and T ∼ = 1 ps, as for these values δ ≈ 10, and Figures 16-18 show that the band of forbidden frequencies does not exist if δ ≥ 3.2 (and c 0 , c 2 and c 3 have the values used to obtain Figures 16-18).However, the value of δ might be smaller in special fibers where v g is closer to c L (as in photonic fibers with air holes), and we consider processes where Z/T is also closer to c L .In such conditions the band of forbidden frequencies might exist.In fact, the generation of a band of forbidden frequencies has already been observed by Fang et al. [38] during the propagation of very short optical pulses in photonic crystal fibers, a result that Fang et al. considered "the most distinctive feature" of this process.
It might be interesting to see how pulses evolve when the width of the band of forbidden frequencies shrinks, as in the dispersion relation shown by the dashed curve in Figure 17.In Figure 19 we can see the modulus of the solution (at z = 1.5) of Equation ( 47) corresponding to a narrow initial pulse of the form (42), with A = 1/2 and B = 14, and the coefficients used to obtain the dashed curve in Figure 17.We can see that the graph shown in Figure 19 is similar to that shown in Figure 10, but the amplitude of the radiation is significantly higher.However, the radiation also stops at a definite time (as in Figure 10), due to the presence of the band of forbidden frequencies shown in Figure 17.
If we now calculate the solution of Equation ( 47) corresponding to the broader initial condition used to obtain Figure 7, and the same coefficients used in Figure 19, we obtain Figure 20.This figure shows that the pulse does not emit any radiation at all, as expected, since the initial condition is now a broad pulse, and its Fourier transform is too narrow to have frequency components which may resonate with the small-amplitude linear waves corresponding to the right branch of the dispersion relation seen in Figure 17. Figure 20 also shows that the position of the pulse has been shifted to the We can see that the graph shown in Figure 19 is similar to that shown in Figure 10, but the amplitude of the radiation is significantly higher.However, the radiation also stops at a definite time (as in Figure 10), due to the presence of the band of forbidden frequencies shown in Figure 17.
If we now calculate the solution of Equation ( 47) corresponding to the broader initial condition used to obtain Figure 7, and the same coefficients used in Figure 19, we obtain Figure 20.This figure shows that the pulse does not emit any radiation at all, as expected, since the initial condition is now a broad pulse, and its Fourier transform is too narrow to have frequency components which may resonate with the small-amplitude linear waves corresponding to the right branch of the dispersion relation seen in Figure 17. Figure 20 also shows that the position of the pulse has been shifted to the left, a result that is due to the presence of the term cu zt Equation ( 47), and the negative value of the coefficient c (Figure 7 shows that the pulse does not move to the left when the term cu zt is absent).We can see that the graph shown in Figure 19 is similar to that shown in Figure 10, but the amplitude of the radiation is significantly higher.However, the radiation also stops at a definite time (as in Figure 10), due to the presence of the band of forbidden frequencies shown in Figure 17.
If we now calculate the solution of Equation ( 47) corresponding to the broader initial condition used to obtain Figure 7, and the same coefficients used in Figure 19, we obtain Figure 20.This figure shows that the pulse does not emit any radiation at all, as expected, since the initial condition is now a broad pulse, and its Fourier transform is too narrow to have frequency components which may resonate with the small-amplitude linear waves corresponding to the right branch of the dispersion relation seen in Figure 17. Figure 20 also shows that the position of the pulse has been shifted to the left, a result that is due to the presence of the term in Equation ( 47), and the negative value of the coefficient (Figure 7 shows that the pulse does not move to the left when the term is absent).

Solutions of Equation (5)
Now let us investigate how the solutions of Equation ( 5) behave.As we saw in Figure 4, when the coefficients , and satisfy the condition (17), the dispersion relation presents two finite bands of forbidden frequencies (defined in (22)).Therefore, we expect that narrow solitons will exhibit a behavior similar to that of the solutions of the Equation ( 4).More precisely, we expect that

Solutions of Equation (5)
Now let us investigate how the solutions of Equation ( 5) behave.As we saw in Figure 4, when the coefficients c 0 , c 2 and c 4 satisfy the condition (17), the dispersion relation presents two finite bands of forbidden frequencies (defined in (22)).Therefore, we expect that narrow solitons will exhibit a behavior similar to that of the solutions of the Equation ( 4).More precisely, we expect that if take an initial condition whose Fourier transform (FT) is wide enough to contain significant frequency components in the intervals ω < −ω 2 and ω > ω 2 , the pulse will start radiating, but the radiation will stop at some time, when the pulse widens and its FT becomes too narrow.To verify this expectation we obtained the numerical solution of Equation ( 5) with coefficients c 0 = 1/15, c 2 = 1/2 and c 4 = 1/256, and using an initial condition u(0, t) of the form (42) with A = 1 and B = 7.As in the case of Equation ( 4), we can obtain the solution of Equation ( 5) by calculating the Fourier transform (FT) of this equation, and the FT of the solution also has the form shown in Equation ( 30), with λ a,b (ω) given by Equation ( 31), but with the function S(ω) which appears within the square root now defined as follows: As this function is positive on two regions of the ω axis, the function U(z, ω) (given in Equation ( 30)) would grow exponentially in z, unless we impose the restriction a(ω) = 0.This restriction implies that Equation ( 41) is also valid in this case (but using Equation (48) in the expression (31) which defines the function λ b (ω)).Therefore, to solve Equation ( 5) we used u(0, t) as defined in Equation ( 42), and we calculated u z (0, t) with Equation (41).In Figure 21 we can see the profile of |u(z = 1.5, t)| corresponding to the initial condition (42) with A = 1 and B = 7.As expected, we see that the pulse starts emitting radiation, but afterwards the radiation stops.The only important difference in comparison with the behavior of the solutions of Equation ( 4) is that the solution of Equation ( 5) emits radiation in both directions, due to the fact that its dispersion relation, shown in Figure 4, is now symmetrical with respect to the wavenumber axis.
Equation ( 30)) would grow exponentially in z, unless we impose the restriction 0. This restriction implies that Equation ( 41) is also valid in this case (but using Equation ( 48) in the expression (31) which defines the function ).Therefore, to solve Equation ( 5) we used 0, as defined in Equation ( 42), and we calculated 0, with Equation (41).In Figure 21 we can see the profile of | 1.5, | corresponding to the initial condition ( 42) with 1 and 7.
As expected, we see that the pulse starts emitting radiation, but afterwards the radiation stops.
The only important difference in comparison with the behavior of the solutions of Equation ( 4) is that the solution of Equation ( 5) emits radiation in both directions, due to the fact that its dispersion relation, shown in Figure 4, is now symmetrical with respect to the wavenumber axis.On the other hand, in Figure 22 we can see the profile of another solution of Equation ( 5) corresponding to an initial condition of the form ( 42), but now using 1 and 2 (and the same coefficients , and ).In this case no radiation is observed because the initial pulse was very wide, and its FT was very narrow, and therefore the initial FT has no significant frequency components in the intervals and (where the resonant frequencies are located).On the other hand, in Figure 22 we can see the profile of another solution of Equation ( 5) corresponding to an initial condition of the form ( 42), but now using A = 1 and B = 2 (and the same coefficients c 0 , c 2 and c 4 ).In this case no radiation is observed because the initial pulse was very wide, and its FT was very narrow, and therefore the initial FT has no significant frequency components in the intervals ω < −ω 2 and ω > ω 2 (where the resonant frequencies are located).We have thus seen that the behavior of the radiating pulses of Equations ( 4) and ( 5) is similar, except that the solutions of Equation ( 5) emit radiation in both directions, and those of Equation ( 4) only radiate to the left.However, Equations ( 4) and ( 5) differ in a significant aspect.In the case of Equation ( 4) the function defined by Equation ( 32) is always positive in a portion of the real axis, and consequently is always necessary to impose the restriction 0 which leads to the relationship between 0, and 0, given by Equation (41).On the other hand, in the case of Equation ( 5), depending on the values of the coefficients , and , the function defined by Equation ( 48) may be positive over certain frequency intervals, and in these cases the restriction 0 is necessary, and the Equation ( 41) again defines a necessary relation between 0, and 0, .However, for certain coefficients, the function might be always negative, and then the values of the functions , would be imaginary for all frequencies.In these cases there is no reason to impose the restriction 0, and with the disappearance of this restriction also the relations ( 40) and ( 41) disappear.Therefore, in the case of Equation ( 5), when and are both imaginary, we face again the question: how do we specify the initial condition 0, ?In the We have thus seen that the behavior of the radiating pulses of Equations ( 4) and ( 5) is similar, except that the solutions of Equation ( 5) emit radiation in both directions, and those of Equation ( 4) only radiate to the left.However, Equations ( 4) and ( 5) differ in a significant aspect.In the case of Equation ( 4) the function S(ω) defined by Equation ( 32) is always positive in a portion of the real axis, and consequently is always necessary to impose the restriction a(ω) = 0 which leads to the relationship between u(0, t) and u z (0, t) given by Equation (41).On the other hand, in the case of Equation ( 5), depending on the values of the coefficients c 0 , c 2 and c 4 , the function S(ω) defined by Equation ( 48) may be positive over certain frequency intervals, and in these cases the restriction a(ω) = 0 is necessary, and the Equation ( 41) again defines a necessary relation between u(0, t) and u z (0, t).However, for certain coefficients, the function S(ω) might be always negative, and then the values of the functions λ a,b (ω) would be imaginary for all frequencies.In these cases there is no reason to impose the restriction a(ω) = 0, and with the disappearance of this restriction also the relations ( 40) and ( 41) disappear.Therefore, in the case of Equation ( 5), when λ a and λ b are both imaginary, we face again the question: how do we specify the initial condition u z (0, t)?In the following paragraph we address this question.
When a(ω) = 0 the Equations ( 40) and ( 41) no longer hold, but we can obtain an approximate relation between u z (0, t) and u(0, t) by considering what happens when we send a light pulse at the beginning of an optical fiber.As the pulse enters into the fiber, the wavelength of the carrier wave changes (since the light's velocity is different inside the fiber), and this change introduces a factor exp(i∆k z) in the mathematical form of the pulse.In other words, at the beginning of the fiber (inside the fiber, near the point z = 0) the form of the light pulse will be: where f (t) defines the temporal profile of the pulse.Consequently, from Equation ( 49) it follows that: and this equation, whose FT is similar to Equation ( 40), implies that the form of the function u(0, t) determines the form of u z (0, t).It should be noticed, however, that we do not know the value of the parameter ∆k.This parameter depends on the exact form of the refractive index of the fiber as a function of the light intensity.However, if we approximate the wavelength of the light inside the fiber in the form λ ∼ = λ 0 /n 0 , where λ 0 is the wavelength outside the fiber, and n 0 is the intensity-independent part of the refractive index, then we can approximate ∆k in the form: Therefore, when there is no justification for using Equation (41) to calculate u z (0, t), this function might be calculated by means of the approximate Equations ( 50) and (51).In fact, in the following section we will see that the exact soliton solutions of Equation ( 6) are indeed consistent with a linear relationship between u z (0, t) and u(0, t), such as that given by Equation (50).
To close this section, it is important to emphasize that the peculiar radiation behavior observed in Figures 9, 10, 19 and 21 is a direct consequence of the presence of bands of forbidden frequencies in the dispersion relations of Equations ( 4), ( 5) and (47), and these forbidden bands are the result of the interplay between the non-SVEA tem u zz and higher-order dispersive terms such as −iu ttt or u 4t .Therefore, the principal lesson to be learnt from the study of Equations ( 4), ( 5) and ( 47) is the following: when a very short pulse is launched along an optical fiber, some of the frequencies which are contained in the initial pulse may not be permitted to travel along the fiber, since bands of forbidden frequencies may appear in the dispersion relation of the equation which controls the propagation of the pulse.Although these forbidden bands may be cancelled by the influence of the term cu zt (which we did not include in Equations ( 4) and ( 5)), a band of this type has already been observed in a photonic crystal fiber, as mentioned before [38].

The Nonlinear Equations (6)-(8)
In this section we will show that the nonlinear Equations ( 6)-( 8) have exact soliton solutions of different forms, some of them being embedded solitons (as we shall explain in the following).

The Solitons of Equation (6)
Equation ( 6) is an interesting generalization of Equation ( 5) because it accepts soliton solutions of two different types.The first of these solutions is the following: where A, B and C are the following constants: Equation ( 55) defines two values for C, and consequently Equation ( 6) has two solutions of the form (52) with different wavenumbers C 1 and C 2 corresponding, respectively, to the signs "plus" and "minus" in front of the radical that appears in Equation (55).
It is worth observing that the solution defined by Equations ( 52)-( 55) is the particular solution of Equation ( 6) corresponding to the initial conditions: with A, B and C given by Equations ( 53)- (55).We can see that Equation (57) has the same form as Equation ( 50), but with a particular value of the wavenumber.Therefore, the use of initial conditions of the form (50) in order to look for particular solutions of Equations ( 5) and ( 6) seems justified.
We must now remember that in Section 2 we saw that if the inequality (17) holds, the range of wavenumbers permitted by the linear dispersion relation of Equation ( 6) (given by Equation (15) or, alternatively, by Equation ( 21)), covers the entire real axis, and consequently, in this case (when c 0 c 2 2 > c 4 ) the soliton wavenumbers C 1 and C 2 are obviously contained within the range of wavenumbers permitted by the linear dispersion relation of the system.This characteristic implies that these solitons are not standard ones (since standard solitons have wavenumbers lying outside the linear spectrum of the system), but they are embedded solitons.It is worth remembering that prior to 1997 it was believed that soliton wavenumbers must necessarily lie outside the linear spectra of the systems.Otherwise (it was believed) the solitons would resonate with the small-amplitude linear waves capable of propagating in the system, and the solitons would emit resonant radiation.However, in 1997 it was discovered that exact radiationless soliton solutions with wavenumbers contained within the range of the linear dispersion relation of the system may indeed exist [11], and the term "embedded soliton" was coined two years later to distinguish these peculiar solitons from the standard ones [45,46].
In Figure 4 we showed the form of the dispersion relation, given by Equation (21), in the particular case when c 0 = 1/15, c 2 = 1/2 and c 4 = 1/256.This figure also shows the position of the soliton wavenumbers, C 1 and C 2 , given by Equation ( 55) with γ 1 = 5 and γ 2 = 1.If we look at this figure an interesting question arises: how is the frequency spectrum of the exact soliton solution (52) in comparison to the dispersion relation?To answer this question, in Figure 23 we have superimposed the Fourier transform (FT) of the exact soliton (52) and the scaled dispersion relation k(ω)/8.6287.The scaling factor (8.6287) −1 was introduced for the amplitude of the FT and the scaled dispersion relation to have the same size at ω = 0.This figure shows that the amplitude of the spectrum of the exact soliton solution is completely insignificant in the range of frequencies where the lateral branches of the dispersion relation are placed.Therefore, even though the two soliton solutions of Equation ( 6) have wavenumbers C 1 and C 2 that intersect these branches (as the dashed lines show in Figure 4), and this fact favors a resonance between the solitons and the radiation modes, the soliton spectrum is practically zero at the frequencies where k(ω) = C 1,2 , thus implying that the solitons do not have frequency components capable of resonating with the radiation modes satisfying the resonance condition k(ω) = C 1,2 .
this figure an interesting question arises: how is the frequency spectrum of the exact soliton solution (52) in comparison to the dispersion relation?To answer this question, in Figure 23 we have superimposed the Fourier transform (FT) of the exact soliton (52) and the scaled dispersion relation /8.6287.The scaling factor (8.6287) −1 was introduced for the amplitude of the FT and the scaled dispersion relation to have the same size at 0. This figure shows that the amplitude of the spectrum of the exact soliton solution is completely insignificant in the range of frequencies where the lateral branches of the dispersion relation are placed.Therefore, even though the two soliton solutions of Equation ( 6) have wavenumbers and that intersect these branches (as the dashed lines show in Figure 4), and this fact favors a resonance between the solitons and the radiation modes, the soliton spectrum is practically zero at the frequencies where , , thus implying that the solitons do not have frequency components capable of resonating with the radiation modes satisfying the resonance condition , .given by Equation ( 21) and 8.6287.
The dashed line shows the Fourier transform of the initial form of the soliton solution ( 52), with A and B given by Equations ( 53) and ( 54).
This qualitative argument is helpful to understand why these embedded solitons do not resonate with the radiation modes, in spite of the fact of having wavenumbers contained in the range of the dispersion relation.However, if desired, it is also possible to use the procedure shown in Refs.[18] and [20] to construct a rigorous proof (not just a qualitative one) which shows that the absence of resonances between the solitons and the radiation modes is the consequence of a delicate balance between the linear and the nonlinear terms of Equation ( 6).This procedure consists in taking the FT of Equation ( 6), using a function of the form (52) to calculate the FT of the nonlinear terms.In this way, an expression for the FT of , is obtained, which contains a quotient of two polynomials in .The polynomial in the numerator depends on the nonlinear coefficients, and that in the denominator depends on the coefficients of the linear terms.Moreover, the roots of the polynomial that appears in the denominator determines the resonant frequencies.The next step is to prove that if the values of the parameters A, B and C which appear in (52) have the values given by the Equations ( 53)- (55), then the two polynomials cancel out, and the denominator disappears.With the disappearance of this denominator, the possibility of a resonance also disappears, thus explaining why the solution given by Equations ( 52)-( 55) does not resonate with the radiation modes.This qualitative argument is helpful to understand why these embedded solitons do not resonate with the radiation modes, in spite of the fact of having wavenumbers contained in the range of the dispersion relation.However, if desired, it is also possible to use the procedure shown in Refs.[18,20] to construct a rigorous proof (not just a qualitative one) which shows that the absence of resonances between the solitons and the radiation modes is the consequence of a delicate balance between the linear and the nonlinear terms of Equation ( 6).This procedure consists in taking the FT of Equation ( 6), using a function of the form (52) to calculate the FT of the nonlinear terms.In this way, an expression for the FT of u(z, t) is obtained, which contains a quotient of two polynomials in ω.The polynomial in the numerator depends on the nonlinear coefficients, and that in the denominator depends on the coefficients of the linear terms.Moreover, the roots of the polynomial that appears in the denominator determines the resonant frequencies.The next step is to prove that if the values of the parameters A, B and C which appear in (52) have the values given by the Equations ( 53)-( 55), then the two polynomials cancel out, and the denominator disappears.With the disappearance of this denominator, the possibility of a resonance also disappears, thus explaining why the solution given by Equations ( 52)-( 55) does not resonate with the radiation modes.
The above paragraph shows that the soliton solutions of Equation ( 6) (defined by Equations ( 52)-( 55)) are embedded solitons if the inequality (17) holds.On the other hand, when the inequality (18) holds, the linear dispersion relation (21) presents a band of forbidden wavenumbers (defined by (19) and ( 20)), and in this case it is not evident if the soliton wavenumbers C 1 and C 2 are located within the forbidden band (in which case the solitons would be standard), or they lie outside this band (and they are contained within the linear spectrum of the system), in which case they would be embedded solitons.A careful comparison of the soliton wavenumbers C 1 and C 2 given by Equation (55), and the boundaries k 1 and k 2 of the band of forbidden frequencies (given by Equation ( 20)), show that: These inequalities imply that the soliton wavenumbers lie outside the forbidden band k 1 < k < k 2 , and consequently there are contained within the range of wavenumbers permitted by dispersion relation, thus implying that also case (when the inequality (18) holds) the solitons defined by Equations ( 52)-( 55) are embedded solitons.In Figure 3 we showed the dispersion relation k(ω) defined by Equation ( 15) when c 0 = 1/30, c 2 = 1/2 and c 4 = 1/80, and the dashed horizontal lines show the position of the soliton wavenumbers C 1 and C 2 (corresponding to γ 1 = 5 and γ 2 = 1), evidencing that these wavenumbers lie outside the forbidden band.Now let us consider the second type of soliton solutions of Equation (6).Direct substitution shows that if the coefficients of the equation satisfy the following relation: then Equation ( 6) has an exact solution of the following form: where: and P is the solution of: which, if c 0 = 0, implies that: The solution given by Equations ( 60)-( 65) is an interesting one, as it contains an unusual nonlinear frequency shift.This type of frequency shift was also found in [17], and if c 0 = 0 our Equations (59)-(64) coincide with the results found in this reference.
It should be noticed that Equation (6) does not accepts simultaneously the solutions ( 52) and (60).The solution (52) exists for certain values of the coefficients c 0 , c 2 , c 4 , γ 1 and γ 2 , and the solution (60) exists for other values of these coefficients.For example, if c 2 and c 4 are both positive, the solution (52) only exists if γ 1 is also positive (as a consequence of Equation ( 53)), while the solution (60) will only exist if γ 1 is negative (as a consequence of Equation ( 61)).
The existence of two different types of solitons in Equation ( 6) suggests that this equation might have additional exact analytical solutions.It would be interesting to investigate if the new methods described in Refs.[47][48][49] and the references therein, might be successful in finding new solutions for Equation (6).

The Solitons of Equation (7)
Direct substitution shows that Equation (7) has moving soliton solutions of the following form: where r, q, V, R and W are constants whose values can be calculated as explained in the following.First we calculate r: r = − c 3 4c 4 (67) and then we define the following parameters: Then we can obtain the possible values of q from the following equation: ( and with each of the possible values of q we can calculate V, R and W with the following expressions: As Equation ( 71) is a sixth-order equation for q, we may have up to six different values for q, and consequently Equation ( 7) might have six different soliton solutions for a given set of coefficients {c 0 , c 2 , c 3 , c 4 , γ 1 , γ 2 }.This result suggests that Equation (7) might have a continuous family of soliton solutions, in which case the solitons defined by Equations (66)-(74) would only be particular elements of this family that can be expressed in terms of hyperbolic secants.Now let us investigate if the solitons defined by Equations (66)-(74) are standard or embedded.Apparently, in order to determine if the soliton (66) is embedded, we would only need to investigate if the soliton's wavenumber q (defined by Equation (71)) is contained in the range of wavenumbers permitted by the linear dispersion relation corresponding to Equation (7), which is given by: However, this is not true.When we deal with a moving soliton (i.e., a soliton which moves along the retarded time axis, as the soliton in Equation (66)), the determination if the soliton is embedded is more involved.As explained in Ref. [18], in this case we must investigate if the intrinsic soliton's wavenumber defined as: is contained in the range of a modified dispersion relation, which is obtained by substituting a linear wave written in the form u = exp[kz − ω(t − Vz)] into the linear part of the equation under study.
In the case of Equation ( 7), this modified dispersion relation takes the form: This equation permits us to write k as function of ω, but now something unexpected occurs.As Equation (77) involves V, and V depends on q (as shown in Equation ( 72)), the dispersion relation (77) will depend on the value of q.In other words: we will have a different dispersion relation for each of the different real values of q defined by Equation (71).Consider, for example, the following set of coefficients: Substituting these values in Equations ( 67)-( 70) we obtain the values of Y, M and L.Then, substituting these values into Equation (71), and solving this equation, we find that there are two complex values for q, and the following four real values: Each of these values defines a dispersion relation given by Equation (77).The form of these dispersion relations can be seen in Figure 24.The dashed horizontal line which appears in each of the four figures indicates the position of the corresponding intrinsic soliton's wavenumber Q i = q i + V(q i )r.The values of these Q i are the following: In each of the four graphs shown in Figure 24 the dashed horizontal line intersects the dispersion relation k(ω), thus implying that the four solitons corresponding to the coefficients shown in (78) are embedded solitons.
of coefficients: 0.025, 0.5, 0.07, 0.02, 5, 1 Substituting these values in Equations ( 67)-( 70) we obtain the values of Y, M and L.Then, substituting these values into Equation (71), and solving this equation, we find that there are two complex values for , and the following four real values: In each of the four graphs shown in Figure 24 the dashed horizontal line intersects the dispersion relation , thus implying that the four solitons corresponding to the coefficients shown in (78) are embedded solitons.
We may wonder if all the solitons of Equation ( 7) are embedded.It is difficult to answer such a question.We have examined the solitons corresponding to other coefficients, different from those shown in (78), and in all the cases that we have studied the solitons turned out to be embedded.However, it is difficult to prove that in general, for every set of coefficients which lead to real values of the parameters M, L, q, R and W (defined by Equations ( 69)-( 74)), the corresponding solitons are all embedded.In fact, it might be impossible to generate such a proof, since we would need to compare the values of the intrinsic solitons' wavenumbers (which depend on the values of the solutions of Equation ( 71)) with the boundaries of the forbidden bands that appear in the dispersion relations, and it is impossible to obtain analytical expressions for in closed form, since Equation ( 71) is a sixth-order equation.

The Solitons of Equation (8)
Equation ( 8) is worth being studied because, in spite of being a particular case of Equation ( 6), it has soliton solutions which are different from those of Equation (6).Direct substitution shows that Equation (8) has exact soliton solutions of the following form: where the parameters D, E and H are defined by following equations: We may wonder if all the solitons of Equation ( 7) are embedded.It is difficult to answer such a question.We have examined the solitons corresponding to other coefficients, different from those shown in (78), and in all the cases that we have studied the solitons turned out to be embedded.However, it is difficult to prove that in general, for every set of coefficients which lead to real values of the parameters M, L, q, R and W (defined by Equations ( 69)-( 74)), the corresponding solitons are all embedded.In fact, it might be impossible to generate such a proof, since we would need to compare the values of the intrinsic solitons' wavenumbers Q i (which depend on the values of the solutions q i of Equation ( 71)) with the boundaries of the forbidden bands that appear in the dispersion relations, and it is impossible to obtain analytical expressions for q i in closed form, since Equation ( 71) is a sixth-order equation.

The Solitons of Equation (8)
Equation ( 8) is worth being studied because, in spite of being a particular case of Equation ( 6), it has soliton solutions which are different from those of Equation (6).Direct substitution shows that Equation (8) has exact soliton solutions of the following form: u(z, t) = D sech 2 (Et) exp(iH 1,2 z) (81) where the parameters D, E and H are defined by following equations: Optical solitons with profiles given by squared hyperbolic secants are not unknown in optics [50,51], but they are far less common than the typical sech-type solitons.Consequently, it is interesting that the solitons of Equation ( 8) have the form given by Equation (81).Moreover, it is also worth observing that even though Equation ( 6) reduces to Equation (8) when γ 2 = 0, neither of the two soliton solutions of Equation ( 6) (given either by Equations ( 52)- (55) or by Equations ( 60)-( 65)) reduces to the soliton solution of Equation (8) given by Equations ( 81)-(84) in the limit when γ 2 → 0 .This is not an unusual situation.It occurs in many systems with embedded solitons.For example, the equation: reduces to the standard NLS equation when ε 2 → 0 and γ 2 → 0 , but it is known that the solitons of Equation (85) do not reduce to NLS solitons in this limit [11].In a similar way, Equation (8) reduces to the NLS equation when c 0 = c 4 = 0, but the solution defined by Equations ( 81)-(84) does not reduce to a NLS soliton in this case.In fact, as c 4 → 0 the soliton's height (given by the parameter D defined by Equation ( 82)) tends to infinity, and the soliton's width (given by the parameter E −1 defined by Equation ( 83)) tends to zero.Are the soliton solutions of Equation ( 8) embedded solitons?This is a question which deserves some attention.In the case when the inequality (17) holds, it is obvious that the solitons defined by Equations ( 81)-(84) are embedded, since in this case the wavenumbers permitted by the linear dispersion relation (15) cover the entire real axis (as explained in Section 2), and therefore the solitons' wavenumbers H 1,2 will necessarily be contained in the range of this dispersion relation.On the other hand, when the inequality (18) holds, the range of the dispersion relation will contain a band of forbidden wavenumbers, and the boundaries of this band are the values k 1 and k 2 given by Equation (20).In this case it is not evident if the solitons defined by Equations (81)-(84) are embedded or not.However, if we choose values of c 0 and c 4 which lead to real values of H 1,2 (when substituted in Equation (84)), an algebraic exercise shows that: Therefore, the soliton solutions of Equation ( 8) are interesting because, depending on the values of c 0 , c 2 and c 4 , they may be embedded solitons or standard ones.end of Section 3.2), that these forbidden bands may be cancelled by the effect of the term cu zt , but there may exist systems where such bands appear.The experimental observation by Fang et al. [38] of a gap of forbidden frequencies in the spectral profile of very short optical pulses propagating in a photonic crystal fiber may be an example of a process of this type.
The main results of this communication are, therefore, the following: • the discovery of a relationship between Equations ( 4), ( 47), (A12) and (A13) and elliptic curves, • the discovery that the interplay between u zz and higher-order dispersive terms generates bands of forbidden frequencies or forbidden wavenumbers in the dispersion relations of Equations ( 4)-( 8) and ( 47), • the discovery that the Cauchy problems associated to Equations ( 4)-( 8) are probably ill-posed, • the discovery that short pulses who evolve according to Equations ( 4), ( 5) and ( 47) radiate in a novel and peculiar way, the discovery of four different types of soliton solutions in Equations ( 6)-( 8), the discovery that some of the solitons of Equation ( 6)-( 8) are embedded solitons, • the discovery that Equations (A12) and (A13) might be related to modular forms or Weierstrass ℘ functions, • the discovery that we might find a relation between the equations studied in this paper and cryptography.
Several of these results suggest further topics of research, and consequently we hope that this article encourages the research along these lines.
In 1975 Yves Hellegouarch discovered that the elliptic curve [56]: where a, b and p > 2 are integers, would have rather unusual properties if a p + b p were also a pth power, i.e., if there existed an integer c such that: STEP 3: Frey's proposition.
In 1984 Gerhard Frey announced in a conference in Oberwolfach, Germany, that if there existed integers a, b, c and p > 2 satisfying (A16), then the elliptic curve (A15) would not have a modular form associated with it [57].In other words: if Fermat's last theorem were false, then the Taniyama-Shimura would also be false, and this implied that: If the Taniyama-Shimura is true, then Fermat's last theorem is also true!Therefore, Frey had discovered that one way to prove Fermat's last theorem was to prove the Taniyama-Shimura conjecture.However, a problem remained: Frey's procedure to "prove" that the elliptic curve (A15) had not a corresponding modular form, contained an error.This problem was solved in 1986 by Ken Ribet, who presented a correct proof of Frey's proposition [57].Therefore, the way to prove Fermat's last theorem was clearly indicated: it was necessary to prove the Taniyama-Shimura conjecture.STEP 4: Wiles' proof.
In 1993 Andrew Wiles announced that he had proved a restricted form of the Taniyama-Shimura conjecture that was sufficient to prove Fermat's last theorem.It turned out that Wiles' proof had a mistake, but in 1994, in collaboration with Richard Taylor, the mistake was corrected, and Fermat's last theorem was finally proved.The full Taniyama-Shimura conjecture (now known as modularity theorem) was proved in 2001 by Christophe Breuil, Brian Conrad, Fred Diamond and Richard Taylor, thus completing the intellectual adventure started by Taniyama and Shimura 46 years before [58].
has the form shown in Figure2.

with 5 and 1 )
figure indicate the values of the wavenumbers of the soliton solutions of Equation (6) (that will be found in Section 4.
Appl.Sci.2017, 7, 340 7 of 32 shown in this figure indicate the values of the wavenumbers of the soliton solutions of Equation (6) (with 5 and 1) that will be found in Section 4.

Figure 11 .
Figure 11.Real part of , 0 , where , is the solution of Equation (4) with the same coefficients and the same initial condition used to obtain Figures 8-10.Figure 11.Real part of u(z, t = 0), where u(z, t) is the solution of Equation (4) with the same coefficients and the same initial condition used to obtain Figures 8-10.

Figure 11 . 32 Figure 12 .
Figure 12.Imaginary part of , 0 , where , is the solution of Equation (4) with the same coefficients and the same initial condition used to obtain Figures 8-10.From these figures it follows that 40, and therefore a horizontal line placed at 40 intersects the dispersion relation shown in Figure 2 at two points: , 40, 0 and , 40, 25 .The first of these points does not correspond to a traveling wave (as its frequency is zero), but the second one corresponds to a linear wave with a resonant frequency 25, phase velocity / 1.6 and group velocity / 12.5 (i.e., in fact we are talking of inverse velocities, Δ /Δ , as mentioned in the previous paragraph).We can see, therefore, that the endpoint

Figure 12 .
Figure 12.Imaginary part of u(z, t = 0), where u(z, t) is the solution of Equation (4) with the same coefficients and the same initial condition used to obtain Figures 8-10.

Figure 14 .
Figure 14.Fourier transform of the curve shown with the continuous line in Figure 13.

Figure 15 .
Figure 15.Fourier transform of the curve shown with the dashed line in Figure 13.

Figure 14 .
Figure 14.Fourier transform of the curve shown with the continuous line in Figure 13.

Figure 15 .
Figure 15.Fourier transform of the curve shown with the dashed line in Figure 13.

Figure 14 .
Figure 14.Fourier transform of the curve shown with the continuous line in Figure 13.

Figure 14 .
Figure 14.Fourier transform of the curve shown with the continuous line in Figure 13.

Figure 15 .
Figure 15.Fourier transform of the curve shown with the dashed line in Figure 13.Figure 15.Fourier transform of the curve shown with the dashed line in Figure 13.

Figure 15 .
Figure 15.Fourier transform of the curve shown with the dashed line in Figure 13.Figure 15.Fourier transform of the curve shown with the dashed line in Figure 13.
) and let us take the same coefficients c 0 = 1/40, c 2 = 1/2 and c 3 = 1/50 used to obtain Figures 7-12.Concerning the coefficient c, one of the reviewers that revised this paper observed that in typical fibers c < 0 and the absolute value of the ratio c/c 0 is approximately equal to the dimensionless number δ = 2v −1 g Z/T, where Z and T are characteristic values of length and time, respectively, and v −1 g is the reciprocal group velocity.In Figures 16-19 the dashed curves show the form of the dispersion relation of Equation (47) for three different values of c (−c 0 , −2.8 c 0 and −3.2 c 0 ), and the continuous curve shows the dispersion relation when the term cu zt is absent.

Figure 24 .
Figure 24.Dispersion relations of the form (77), corresponding to the four values of shown in (79), which are the real roots of Equation (71) with the coefficients and shown in (78).The dashed lines indicate the values of the corresponding intrinsic solitons' wavenumbers given in (80).

Figure 24 .
Figure 24.Dispersion relations of the form (77), corresponding to the four values of q i shown in (79), which are the real roots of Equation (71) with the coefficients c i and γ i shown in (78).The dashed lines indicate the values of the corresponding intrinsic solitons' wavenumbers Q i given in (80).
Each of these values defines a dispersion relation given by Equation (77).The form of these dispersion relations can be seen in Figure24.The dashed horizontal line which appears in each of the four figures indicates the position of the corresponding intrinsic soliton's wavenumber .