Quasi-integrability in the modified defocusing non-linear Schrödinger model and dark solitons

The concept of quasi-integrability has been examined in the context of deformations of the defocusing non-linear Schrödinger model (NLS). Our results show that the quasi-integrability concept, recently discussed in the context of deformations of the sine-Gordon, Bullough-Dodd and focusing NLS models, holds for the modified defocusing NLS model with dark soliton solutions and it exhibits the new feature of an infinite sequence of alternating conserved and asymptotically conserved charges. For the special case of two dark soliton solutions, where the field components are eigenstates of a space-reflection symmetry, the first four and the sequence of even order charges are exactly conserved in the scattering process of the solitons. Such results are obtained through analytical and numerical methods, and employ adaptations of algebraic techniques used in integrable field theories. We perform extensive numerical simulations and consider the scattering of dark solitons for the cubic-quintic NLS model with potential V=ηI2−∈6I3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ V=\eta {I}^2-\frac{\in }{6}{I}^3 $$\end{document} and the saturable type potential satisfying , with a deformation parameter ϵ ∈ and I = |ψ|2. The issue of the renormalization of the charges and anomalies, and their (quasi)conservation laws are properly addressed. The saturable NLS supports elastic scattering of two soliton solutions for a wide range of values of {η, ϵ, q}. Our results may find potential applications in several areas of non-linear science, such as the Bose-Einstein condensation.


Introduction
The theory of soliton collisions in the integrable models is well known in the current literature. In particular, the solitons emerge with their velocities and shapes completely unaltered after collision between them, the only noticed outcome being their phase shifts. However, certain non-linear field theory models with relevant physical applications possess soliton-like solutions and it is difficult to know a priori if they are in fact true solitons. The so-called quasi-integrability concept has recently been put forward in the context of certain deformations of the integrable sine-Gordon (SG), Bullough-Dodd (BD) and focusing non-linear Schrödinger (NLS) models [1][2][3][4]. The idea is that many non-integrable theories posses solitary wave solutions that behave similar to solitons, i.e. the scattering of such JHEP03(2016)005 solitons basically preserve their shapes and velocities. So, certain deformations of the SG, BD and focusing NLS models were considered as quasi-integrable theories possessing infinite number of charges that are asymptotically conserved. This means that during the collision of two soliton-like solutions such charges do vary in time and sometimes significantly, however their values remain the same in the past (before the collision) and in the future (after the collision). The deformations of the topological SG and BD models studied in [1][2][3] furnish kink-type solitary waves presenting soliton-like behaviour. Whereas, the deformed focusing NLS model considered in [4] supports bright solitary waves presenting soliton-like behaviour. The both numerical and analytical techniques were used to describe the concept of quasi-integrability.
Regarding the analytical calculations of solitary wave collisions in deformed NLS models some results have been obtained for the cases of small perturbations of the integrable NLS model [5][6][7][8]. For nearly integrable models it has been implemented a perturbation theory allied to the the inverse scattering transform (IST) method [9]. These methods have been applied mainly to the focusing NLS solitons that decay at infinity, i.e. the bright solitons. On the other hand, the integrable defocusing NLS model supports dark solitons which are intensity dips sitting on a continuous wave background (cwb) with a phase change across their intensity minimum. In order to apply to the dark soliton perturbation and the non-vanishing boundary conditions (nvbc), inherent to this type of solution, various improvements have been made to the calculations and methods [10][11][12] based on the so-called complete set of squared Jost solution (eigenstates of the linearised NLS operator). The implementation of direct methods, however, consider small perturbations around the integrable NLS model [8,13]. Moreover, in [8] it has been argued that the squared Jost functions associated with the dark soliton might be an insufficient basis, so rendering problematic the issue of the completeness of the basis for the solution space of a deformed NLS equation. The solitary wave collisions of certain deformed NLS models, beyond small perturbations, have been considered in the recent literature; e.g. in [14] it has been computed the spatial shifts of fast dark-dark soliton collisions and in [15] some properties of slow dark soliton collisions have been studied.
It would be interesting to understand the issue of quasi-integrability in the other variant of the NLS model, i.e. the deformed defocusing NLS model, since such theories also appear in many areas of non-linear science, condensed matter physics, plasma physics and, in particular, in the study of Bose-Einstein condensates. In this context, our aim is to predict the results of solitary wave collisions and test the quasi-integrability concept in deformed defocusing NLS models with nvbc. However, the non-vanishing boundary condition inherent to the dark solitons, which may change when the deformations are present, introduces some complications when applying the techniques developed for bright solitons in order to study the quasi-integrability concept. In the analytic treatments the nvbc introduces the issue of the renormalization of the charges and anomalies in the expressions of the quasi-conservation laws which must be properly addressed. This is closely related to the fact that the dark soliton solution is a composite object, i.e. it incorporates the continuous wave background (cwb) plus the dark soliton itself. In addition, to simulate the time dependence of field configurations for computing dark solitons properties the choice of the efficient and accurate numerical methods, the so-called time-splitting cosine pseudo-JHEP03(2016)005 spectral finite difference (TSCP) and the time-splitting finite difference with transformation (TSFD-T) methods [16,17] will be made in order to control the highly oscillatory phase background. In fact, these methods allowed us to improve in several orders of magnitude the accuracy in the computation of the charges and anomalies presented in [4].
The deformed focusing NLS with bright soliton solutions and the structures responsible for the phenomenon of quasi-integrability have been discussed in [4]. In this context it has been shown that this model possesses an infinite number of asymptotically conserved charges. An explanation found so far for this behaviour of the charges is that some special soliton type solutions are eigenstates of a space-time parity transformation. On the other hand, the deformed defocusing NLS with dark soliton solutions also presents the above remarkable properties, however in this case we can say even more. As we will show, there are some soliton like solutions which present a special space-reflection symmetry for any time, so this property implies the exact conservation of the sequence of the even order charges, which in the absence of this special symmetry would be conserved only asymptotically, as in the models studied in [1][2][3][4]. Here we will show the connection between the space-reflection parity and the exactly conserved charges which involves an interplay between the spacereflection parity and internal transformations in the affine Kac-Moody algebra underlying the anomalous Lax equation. However, it seems to be that such parity property is a sufficient but not a necessary condition in order to have the sequence of exactly conserved charges. In fact, as we will show by numerical simulations, there are certain soliton like configurations without this symmetry which also exhibit such conserved charges. So, these properties constitute the distinguishing new features associated to the deformed defocusing NLS with dark soliton solutions, as compared to the previous quasi-integrable models [1][2][3][4].
The paper is organized as follows. In the next section we discuss the deformed defocusing NLS model with non-vanishing boundary conditions. In section 3, we discuss the concept of quasi-integrability for deformed defocusing NLS, based on an anomalous zero curvature condition. In 3.1 and 3.2 we discuss the relationships between the space-time parity and asymptotically conserved charges, and the space-reflection parity and the exactly conserved charges, respectively. In section 4 we perform the renormalizations of the charges and anomalies relevant for dark soliton solutions. In sections 5.1 and 5.2 we discuss the space-time and space-reflection symmetries of dark soliton solutions of the integrable NLS model, respectively. In 6 we show the vanishing of the anomalies associated to solitary waves. In 6.1, by using general symmetry considerations, we show that the anomalies β (n) r vanish when evaluated on a general solitary wave of the deformed NLS model. In 6.2 we show, by direct computation, the vanishing of the first non-trivial anomaly β (4) r when evaluated on a general solitary wave of the deformed NLS. In section 7 we present the results of our numerical simulations which allowed us to compute and study various properties of dark soliton scattering of the full equations of motion for two deformations of the NLS model and several values of the deformation parameter . In 7.1 we present the results of the simulations of the first deformation with V = η|ψ| 4 − 6 |ψ| 6 , and in 7.2 the second deformation with V [I] = 2ηI − I q 1+I q . In 8 we present some conclusions and discussions. The appendix A presents the eqs. of motion in the R and ϕ parametrizations, as well as some useful identities and B presents relevant expressions corresponding to the gauge transformation.

JHEP03(2016)005 2 Deformations of defocusing NLS
We will consider non-relativistic models of a complex scalar field in (1 + 1)−dimensions with Lagrangian density given by whereψ stands for the complex conjugation of ψ and V : IR + → IR. The Euler-Lagrange equation following from (2.1) becomes The model (2.2) defines the deformed NLS model and it supports dark soliton type solutions in analytical form for some special functions V [I]. The potential V [I] = ηI 2 , (η > 0), corresponds to the integrable NLS model and supports N-dark soliton solutions. The potential V [I] = ηI 2 − I 3 /6 defines the non-integrable cubic-quintic NLS model (CQNLS) which possesses analytical bright and dark type solitons [18,19]. In [18,20] the bright solitary waves of the cubic-quintic focusing NLS have been regarded as quasi-solitons presenting partially inelastic collisions in certain region of parameter space. Among the models with saturable non-linearities [9], the case V [I] = 1 2 ρ s (I + ρ 2 s I+ρs ) also exhibits analytical dark solitons [21]. The deformed NLS model with V [I] = 2ηI − I q 1+I q will deserve a careful consideration below. It passes the Painlevé test for arbitrary positive integers q ∈ Z Z + and = 1 [22]. However, its Lax pair formulation and analytical solutions for a general set {η, , q}, to the best of our knowledge, are not known in the literature.
The qualitative properties of travelling waves of the NLS model for general nonlinearities and non-vanishing boundary conditions have been studied in [23]. Exact analytical dark-soliton solutions of the deformed NLS model with arbitrary potential V are not available, and one can resort to numerical simulations to obtain such solutions. However, dark-soliton solutions can be presented in an implicit form as a curve in a complex plane (for the complex variable Z = A + iB, where ψ ≡ (A + iB)e (kx+wt) ) (see [9] and references therein) or in a small-amplitude approximation where the deformed NLS model is reduced to the Korteweg-de Vries equation [24].
Among the possible deformations of the NLS model the case V = 2η 2+ |ψ| 2 2+ , ∈ IR, for η < 0 has recently been considered in [4] in order to study the concept of quasiintegrability for bright soliton collisions. An analytical solution with vanishing boundary condition (bright soliton) for this potential is well known in the literature. However, an analytical dark soliton solution, as far as we know, is not known in the literature.
In this paper we will study analytically and numerically some deformations of the NLS model of the type (2.2). The first deformation we will consider in our study is defined by the non-integrable cubic-quintic NLS model dy 2 ). This deformation possesses analytical solitary waves as we will present below.
The second deformation considers a saturable non-linearity and it is defined by the equation (2.4) We do not know any analytical solutions for a general set of parameters of this model, however dark solitary wave solutions will be obtained numerically.
The both deformations reproduce the integrable defocusing NLS model in the limit → 0.

Quasi-integrability of deformed NLS
Next we undertake an analytical study of the properties of the deformed NLS model (2.2), in order to do this we will use the well known techniques from the integrable field theories. We follow the developments and notations put forward in [1][2][3][4] on quasi-integrability. Then, one considers an anomalous zero curvature representation of the deformed NLS model (2.2) with the connection given by where the above Lax potentials are based on a sl(2) loop algebra with commutation relations One can verify that the curvature of the connection (3.1) is given by Note that when the equation of motion (2.2) and its complex conjugate, i.e −i∂ tψ + ∂ 2 xψ −ψ δV δ|ψ| 2 = 0, are satisfied the terms proportional to the Lie algebra generators T 0 ± vanish. In addition, the quantity X vanishes for the usual non-linear Schrödinger potential Then, the curvature vanishes for the NLS model making this theory an integrable field theory. For the generalizations of this theory, i.e. deformations of the potential (3.6), the JHEP03(2016)005 above curvature does not vanish and the model is regarded as non-integrable [4]. Next, we use some algebraic techniques borrowed from the theory of integrable models in order to analyse the equation (3.4) corresponding to a model with non-vanishing anomaly X (3.5), or equivalently F xt in (3.3). Let us denote and parametrize γ = i √ ηe iα , γ = −i √ ηe −iα , η ≡γγ. Substituting the parametrization (3.7) into (2.2) one gets the system of eqs. of motion We specialize the parametrizations for the case of the defocusing NLS, i.e. the case η > 0. In addition, define a new basis of the sl(2) loop algebra as The connection in the new basis is obtained through the gauge transformation So, inserting the connection (3.1) into the above expression one gets So, considering the fields which satisfy the equation (2.2) and its complex conjugate, i.e. the terms in the directions of T 0 ± vanish in (3.4), the curvature becomes (3.14) Next, in order to construct the quasi-conservation laws we perform a further gauge transformation

JHEP03(2016)005
The parameters ζ (−n) i of the transformation will be determined below so that the a x component of the new connection lies in the abelian sub-algebra generated by the elements b n . Under the gauge transformation (3.15) the curvature (3.14) transforms as (3.17) The sl(2) loop algebra is furnished with an integer gradation associated to the grading operator d defined by . . . where As usual we will use the generator b 1 as a semi simple element in the sense that its adjoint action decomposes the sl(2) loop algebra G into its kernel and image as follows such that Ker and Im have no elements in common. According to the commutation relations in (3.11) one has that the b n 's form a basis of Ker, and F n i , i = 1, 2, a basis of Im. Next, one can make all the components of a x in the subspace Im to vanish by choosing conveniently Notice that in this process it has not been used neither the equation of motion (2.2) nor its complex conjugate. The components a (3,−n) x will depend on the real fields R and ∂ x ϕ, and their x-derivatives, but not on the potential V . In the appendix B we provide the first components of a x .
Next, the component a t of the connection does not get rotated onto the subspace Ker even when the equation of motion (2.2) and its complex conjugate are used. Then, one has Since a x lies in Ker one has that the term [a t , a x ] in (3.17) has components only in Im. Let us denote then, we find that In the appendix B we present the first quantities α (3,−n) . Next, let us discuss the nvbc associated to dark soliton solutions. These solutions present a difficulty in the sense that the field ψ does not vanish at x → ±∞. In general, the behaviour of the solutions are such that where θ 0 is a constant phase. The connection (3.1) is not suitable to construct the charges because However, making the gauge transformation (3.12) and the reparametrization (3.7) the new connection A µ depends upon the modulus R and its derivatives, and only on the derivatives of the phase ϕ. Therefore, one has that for dark solitons (3.28) In addition the abelianized potential a µ in (3.15) also satisfy that condition, i.e. in the appendix B and one notes that g depends on R and its derivatives and only on the derivatives of the phase ϕ.
Then, in the defocusing NLS case with dark soliton solutions one has that the a t component of the connection satisfies a non-vanishing boundary condition such that a (3.31) Thus, the non-vanishing of the quantity X given in (3.14) and the anomalies β n above imply the non-conservation of the charges. Therefore, the charges and anomalies in (3.31) are valid for the deformed NLS model with nvbc and dark soliton solutions.

Space-time parity and asymptotically conserved charges
Next, let us discuss a special symmetry which plays an important role in the study of quasi-integrable theories. For certain solutions of the theory (2.1) the charges Q (n) satisfy the so-called mirror symmetry. This symmetry is realized on a special subset of solutions of the deformed NLS model (2.2). For every solution belonging to this subset one can find a point (x ρ , t ρ ) in space-time and define a parity transformation around this point as such that the fields R and ϕ transform as In order to realize this symmetry at the level of the anomalous zero-curvature equation let us combine the above parity transformation with the order two Z Z 2 automorphism of the sl(2) loop algebra given by This property can be used to show that the group element g which enter in the gauge transformation (3.15) is even, i.e. Ω(g) = g . For more details of this demonstration we refer to [4].
Next, let us note that the X in the integrand of the anomaly β n in (3.31) is a x−derivative of a functional of R. Since R is assumed to satisfy (3.33) one has that the time-integrated anomaly vanishes where the integration is meant to be performed on any rectangle with center in (x ρ , t ρ ), and the parameters t o , x o which define the integration intervals are any given fixed values of the shifted space-time coordinates t and x, respectively, defined in (3.32).
A remarkable consequence of the vanishing of the integrated anomaly is that the charges (3.31) satisfy a mirror time-symmetry around the point t ρ This relationship shows that the scattering of two-soliton solutions may present an infinite number of charges which are asymptotically conserved. For some very special two-bright soliton configurations in the context of focusing deformed NLS it has been verified that the charges are asymptotically conserved, i.e. even though they vary in time during collision their values in the far past and the far future are the same [4]. As we will show below, for the case of defocusing NLS it is a remarkable fact that for any two-dark soliton configuration it is possible to find a point in space-time (x ρ , t ρ ), such that the fields R and ϕ always satisfy the space-time parity symmetry (3.32)-(3.33) without any additional restriction on the soliton parameters. Next we will summarize the main steps in the construction of solutions in perturbation theory, for a detailed account see [4]. One can choose as a zero order solution a solution of the NLS equation, such that it satisfies (3.33), i.e. it exhibits only the pair R + 0 , ϕ − 0 and nothing of the pair R − 0 , ϕ + 0 . So, let us expand the fields in powers of the deformation parameter

JHEP03(2016)005
Then the equations of motion for the components R + 1 , ϕ − 1 of the first order solution satisfy non-homogeneous equations and so they can never vanish. However, the pair R − 1 , ϕ + 1 satisfy homogeneous equations and so they can vanish. In addition, if the pair of components without definite parity (R 1 , ϕ 1 ) is a solution, so is the pair Then, one can always choose the components of the first order solution whith definite parity. For this choice one can show that the second order solution has also similar properties, i.e. the components R + 2 , ϕ − 2 satisfy non-homogeneous equations and the components R − 1 , ϕ + 1 satisfy homogeneous equations. This enables one to choose the second order solution to be the pair R + 2 , ϕ − 2 with definite parity, and this process repeats in all orders. Summarizing, we may conclude that 3. Then, one has an infinite set of asymptotically conserved charges, i.e.
Therefore, the values of the charges in the infinite past, before the collision of the solitons, are the same as in the infinite future, after the collision. Theories possessing such properties are dubbed as quasi-integrable theories. In particular, the deformed NLS models (2.3) and (2.4) can be shown to satisfy the requirements (3.44) and (3.45) for some field configurations, and then they belong to the class of quasi-integrable theories, as we will discuss in this work.

Space-reflection symmetry and conserved charges
Next, we consider some special solutions which exhibit a space-reflection symmetry such that the fields R and ϕ transform as As we will see below, in the special case of 2-dark solitons moving in opposite directions and equal velocities, such that they undergo a head-on collision, we have a space-reflection JHEP03(2016)005 symmetry. The implication of this additional symmetry of the fields under space-reflection (for any shifted time), i.e. R and ϕ being even fields, on the behaviour of the quantities α (3,−n) deserves a further analysis. We will show that some α (3,−n) 's possess an even parity under P x for any shifted time. Since X is odd the anomaly β n in (3.31) will vanish identically providing an exact conservation equation for the charges dQ (n) dt = 0, n = 0, 2, 4, . . .. Therefore, these type of charges will not vary during collision, implying their exact conservation in the whole process of scattering.
In order to realize the symmetry (3.48)-(3.49) at the level of the anomalous zero curvature equations (3.14) and (3.17) we introduce a transformation Σ in the internal space along with the space-reflection symmetry P x . So, let us consider another order two Z Z 2 automorphism Σ of the sl(2) loop algebra given by Note that the connection component A x in (3.13) does not have a definite space-reflection P x parity. In addition, it doest not lie in an eigenspace of the automorphism Σ of the loop algebra. So, consider the joint action of the space and internal transformations as Then, the x−component of the connection (3.13) transforms as Actually, each term of A x in (3.13) possesses this property. Next, we analyse the transformation property of the new connection a x in (3.19). Thus, one can have The term A (0) x in the third eq. of (3.19) as defined in (3.20) is odd under Ω, then it follows from the third equation in (3.19) that The F (−n) have been selected to rotate a x into the kernel of the adjoint action of b 1 and the r.h.s. of (3.54) is in the image of the adjoint action, so for consistency the both sides of (3.54) must vanish Therefore, F (−1) is even under Ω. Taking into account this result we see from the fourth equation in (3.19) that Then, following similar arguments we conclude that

JHEP03(2016)005
and so that F (−2) is even under Ω as well. Furthermore, from the fifth equation in (3.19) one has that (1 + Ω)a , and following similar arguments as above we find that (1 + Ω)a (−2) By repeating the same arguments as above one can show that all F (−n) are even under Ω. This property can be used to show that the group element g which enter in the gauge transformation (3.15) is even, i.e. Ω(g) = g . (3.59) Since A x and ∂ x are odd under Ω, and since g is even one can show that a x in (3.15) is odd under Ω. In addition, since the F (−n) are even under Ω and the generators satisfy (3.50), it follows from (3.16) that Ω(ζ , a = 1, 2. One can verify these results by inspecting the explicit expressions for the fields ζ −n a provided in the appendix B. In order to find the parity of the factor α (3,−n) in the integrand of the anomaly β n in (3.31) let us consider some properties of the sl(2) loop algebra. The Killing form is given by It is realized by Tr( ) = 1 2πi dλ λ tr( ), where tr is the usual finite matrix trace, and T n i = λT i , i = 3, ±. So, from (3.24) one has that where the invariance property of the Killing form under Σ has been used, and the fact that the b n transform as in (3.50). Then, taking into account (3.59) one has that and so one can see that the α (3,−n) 's with n = 0, 2, 4, . . . are even under P x . On the other hand, since we have assumed that R is even under P x and that X given in (3.14) is a x−derivative of a functional of R, one has that X is odd under P x , i.e. P x (X) = −X and so that Therefore, the anomalous conservation laws (3.31) imply Consequently, the even order charges are exactly conserved. Notice that it holds for any potential V which depends only on the modulus |ψ|. The only requirement is that the field components satisfy (3.49).

JHEP03(2016)005
Next let us analyse the solutions which admit the system of equations (3.8)-(3.9) such that the field components R, ϕ satisfy the space-reflection symmetry (3.48)-(3.49). We perform this construction in perturbation theory around solutions of the NLS model, so we expand the equations of motion and the solutions into power series in , as The deformed potential and its gradient can been expanded as and An interesting solution is the one satisfying R = R + , ϕ = ϕ + , i.e. the odd components vanish R − = 0, ϕ − = 0. We will show that this property is maintained in all orders of perturbation theory. So, substituting the solution Then, the solution of (3.70)-(3.71) with vanishing odd components will satisfy Now, assuming (R 0 , ϕ 0 ) = R + 0 , ϕ + 0 one has that the first order terms O( ) of the even and odd parts of the eqs. of motion (3.8)-(3.9) satisfy .

JHEP03(2016)005
Some comments are in order here. First, the pair of fields R + 1 and ϕ + 1 satisfies the linear system of equations (3.73) and (3.75) with variable coefficients, similarly the pair R − 1 and ϕ − 1 satisfies the linear system of equations (3.74) and (3.76). The couple of first order fields (R + 1 , ϕ + 1 ) and (R − 1 , ϕ − 1 ) in the eqs. above form two decoupled system of equations, respectively. Second, the equations of motion for the even components R + 1 , ϕ + 1 , (3.73) and (3.75), satisfy a non-homogeneous coupled system of equations (the non-homogeneous term being the last term in (3.75 ) and so they can never vanish.
However, the odd part components R − 1 , ϕ − 1 , (3.74) and (3.76), satisfy a homogeneous system of equations and so they can vanish. Third, if the pair of components without definite parity (R 1 , ϕ 1 ) is a solution, so is the pair Then, one can always choose a first order solution which is even under the space-reflection parity, i.e.
Choosing the zeroth (3.72) and first order (3.77) solutions one has that the second order terms O( 2 ) of (3.8)-(3.9) splits under P x as .

JHEP03(2016)005
Then one can choose again the second order solution to be even, and this process repeats in all orders O( n ). By repeating this procedure, order by order, one can construct a perturbative solution which satisfies (3.48)-(3.49), and so has charges satisfying (3.64). Therefore, the theory (2.2) possesses a subset of solutions such that the charges Q (n) of order n = 0, 1, 2, 3, 4, 6, 8, . . ., are exactly conserved. In fact, the first four charges n = 0, 1, 2, 3 will be related below to the phase jump, normalization, momentum and energy of the deformed NLS model. The set of even order conservation laws (3.64) is extended to include the odd cases n = 1, 3, as we will show in the next section. Let us summarize our results. The deformed NLS model presents an infinite number of asymptotically conserved charges as in (3.41) for solitons satisfying the space-time symmetry (3.32)-(3.33). In addition, for solitons satisfying the same space-time symmetry, as well as the space-reflection symmetry (3.48)-(3.49) one can say more. In this case, the sequence of the even order charges become inded exactly conserved (3.64). So, the model supports an infinite number of alternating conserved and asymptotically conserved charges associated to soliton solutions satisfying the both space-time and space-reflection symmetries.

Renormalized charges and anomalies
The difficulty with dark solitons is the fact that the charges and anomalies (3.31) diverge when evaluated on the dark-solitons. However, the divergence comes from the vacuum behavior of the solution at x → ±∞, i.e. the cwb solution given in (3.26). So, one can renormalize them by substructing the vacuum contribution (which is infinite). Using the notation ψ = √ Re iϕ/2 , we consider the following boundary condition The nvbc (4.1) is suitable for soliton solutions with the continuous wave background such that R = |ψ 0 | and ϕ/2 = kx + wt + x 0 .
The vacuum solution is given in (3.26) and we have So, we define the renormalized charges and anomalies and they satisfy dQ Let us evaluate the first charges and anomalies in (3.31).

JHEP03(2016)005
The anomaly does not get renormalized and it vanishes The charge Q r is associated to the phase difference (or the phase jump) of the solutions.

The Q
(1) and This charge is related to the space translation symmetry of the model and using the expression (3.7) it can be written as (4.14) Since the solution ψ incorporates the cwb solutions along with the dark soliton itself, it is a composite object. The momentum associated to the dark soliton itself must be renormalized even for k = 0, and it is provided by [9] where we have used the eq. of motion for R given in the appendix. Therefore (3.31) implies We then define the fake charge as Discarding the total derivative term ∂ in (4.16) and using the relationship and so the renormalized charge becomes (4.24) So, the renormalized Hamiltonian of the usual NLS model (V = η|ψ| 4 ) can be written as Therefore, the NLS Hamiltonian can be written as a linear combination H NLS = Q r . Notice that the expressions (4.15) and (4.25) evaluated for a NLS 1-dark soliton solution will provide a standard relationship ∂H ∂P = v as in classical mechanics, where v is the velocity of the soliton. This indicates that the dark soliton can be considered as an effective particle [9].

JHEP03(2016)005
Therefore we have that the renormalized charge and anomaly are Note that β the anomaly β (4) does not get renormalized. So, consider the renormalized quasi-conservation law This is the first non-trivial non-vanishing anomaly associated to an asymptotically conserved charge for solutions satisfying the symmetry (3.33). Moreover, one can have vanishing anomaly β (4) r associated to the exactly conserved charge Q (4) r for solutions satisfying the symmetry (3.49). In addition, we will show below that this anomaly vanishes for a general solitary wave solution of the deformed NLS.
Discarding a total derivative associated to the first term in (4.32) (i.e. the term (V − 2η) 2η∂ x R 3 must vanish once integrated in the whole line taking into account the b.c. (4.1)) the anomaly β (−4) r can be rewritten in a form amenable to numerical calculations. Then, discarding the 'surface' term, one is left with the expression For later purposes we defined the integrand function γ(x, t). The time integrated anomaly becomes One notices that in the limit → 0 the anomaly β (4) r vanishes identically since V [R] → 2η in this limit. In addition, this anomaly vanishes when |ψ| = |ψ 0 |, i.e. when evaluated on the cwb. We will compute numerically the anomaly β Next we discuss the both space-time and space-reflection symmetries in the general one and two dark soliton solutions of the integrable defocusing NLS model. So, consider the defocusing NLS equation Here we consider the Hirota method and construct the one and two dark solitons. Let us introduce the Hirota tau functions as [25] where h(x, t) = |ψ 0 | exp [i(kx + wt + x 0 )] and |ψ 0 |, k, w, x 0 are real constants. The tau function G is a complex function, whereas F is real. Notice that for F = G = 1 in (5.2) one has the solution which is a continuous wave background (cwb) solution of (5.1) provided that w = −k 2 − 2η|ψ 0 | 2 . Substituting the expression (5.2) into (5.1) one gets The 1-dark soliton solution is given by where y is a complex variable. The functions θ,θ are defined as where the bar stands for complex conjugation. The parameters satisfy the relationships Putting these tau functions into the NLS function ψ in (5.2) and making use of the identity This solution possesses three arbitrary real parameters, say |ψ 0 |, k and the phase associated to y. In fact, from the relationships (5.9) one can get the remaining ones w, p in terms JHEP03(2016)005 of the set {|ψ 0 |, k, y}. From the third relationship in (5.9) one gets the condition (assume k > 0, p = p R + i p I ) |p I + k| ≤ √ η |ψ 0 |. (5.11) The intensity function |ψ| moves at the velocity v = −2p I , which is the velocity of the dark soliton. The dark soliton approaches constant amplitude |ψ 0 | as |x| → ±∞. As x varies from −∞ to +∞ the soliton acquires a 2δ phase, where this quantity is the phase of the parameter y y = e 2iδ . (5.12) We restrict −π/2 < δ < π/2. Moreover, at the center of the soliton, i.e. for θ(x, t) + θ(x, t) + ρ ≡ 0, one has that the intensity becomes This center intensity is lower than the asymptotic amplitude |ψ 0 | and this property characterizes a dark soliton. Notice that this center intensity is controlled by the parameter δ; i.e. this parameter defines the "darkness" of the soliton. From (5.11) one can see that the velocity v = −2p I of the dark soliton satisfy the where v s is the sound speed. Therefore the speed of the dark soliton must satisfy The 2-dark soliton solution is given by where r is a real constant, and y j (j = 1, 2) are constant phases. The functions θ j ,θ j are defined as where the bar stands for complex conjugation. The parameters satisfy the following relationships For later purposes let us record the next relationship between the soliton parameters in the above 2-dark soliton solution where v j = −2p jI , a j = 2p jR , p j ≡ p jR + ip jI , j = 1, 2. The v j are the soliton velocities and the a j are the parameters associated to the width (∼ 1 |a j | ) of each soliton.

Space-time parity transformation
Let us verify the parity properties of the above solutions under the transformation 1-dark soliton. Let us consider the 1-dark soliton in (5.10) rewritten as Writing the function as ψ = √ Re i ϕ 2 such that ϕ = 2(kx + wt + x 0 ) + 2φ + 2δ, φ ≡ arctan tan δ tanh a x + b t , (5.24) we can show that P : R → R, ϕ → −ϕ + const.
The last expression can be rewritten as Notice that the parameters x ρ , t ρ and Λ 1 , Λ 2 are real. The phase φ of N in (5.29) can be written as Therefore it can directly be verified that the squared absolute value R = |ψ| 2 and phase ϕ = 2(kx + wt It must be emphasized that the above parity symmetry is satisfied for any set of values {δ 1 , δ 2 } characterizing the 'darkness' of each one-dark soliton component. In fact, there is no further restriction on the parameter space of this two-dark soliton solution in order to satisfy the parity property. Therefore, for this type of solutions it is expected that the integrated anomalies to vanish and the associated charges to be asymptotically conserved. This behaviour is in contradistinction from the one described in the two-bright analytical soliton solution of the focusing NLS model, in which only a subset of the general solution satisfies that symmetry (i.e. the solutions such that the 2-soliton parameter c satisfies c = n π 2 , n ∈ Z Z in the eq. (4.23) of [1][2][3]). In [1][2][3] the initial conditions for the simulations have used c = 0 in order to show the vanishing of the integrated anomalies associated to asymptotically conserved charges for the collision of bright solitons in the deformed focusing NLS model. The c = 0 cases exhibit non-vanishing integrated anomalies for the soliton collisions.
The general form of the one and two dark soliton solutions as presented in (5.23) and (5.27), respectively, have been derived in [25,26]. In the first reference a hybrid of the dressing transformation and tau function methods has been used, whereas in the second one they have been derived by the KP-hierarchy reduction method. The ansatz (5.2) for the general N-soliton solution in terms of the tau functions G and F follows the construction presented in [25].

Space-reflection parity transformation
It is interesting to see the behaviour of the fields R and φ for the two-soliton solutions possessing a space-reflection symmetric amplitude R(− x, t) = R( x, t) and phase φ(− x, t) = φ( x, t) for any given shifted time t. In fact, when the continuous wave background is at rest, i.e. k = 0, one has the parameter relationships δ 2 = δ 1 , p 2R = −p 1R , p 2I = −p 1I , which satisfy the eq. (5.21) for two-soliton solutions. These parameters correspond to two JHEP03(2016)005 φ( x, t) = arctan sin 2δ 1 sinh t 1 e −∆/2 cosh x 1 + cos 2δ 1 cosh t 1 . (5.40) In figure 1 we plot the functions R and φ for three successive times. Note that they are symmetric under the reflection x → − x for any value of the variable t = t. There is another solution with similar behaviour for the set of parameters {k = 0, Therefore, it is clear that the above solution satisfies This type of solution will be associated to the even order exactly conserved charges as discussed in section 3.2. In our numerical simulations we will verify the properties for this kind of solutions of deformed NLS models, as well as for configurations without this type of symmetry.

JHEP03(2016)005 6 Solitary waves and anomalies
Next we discuss the solitary wave solutions of the model (2.2) and their associated symmetries and vanishing anomalies. Below we discuss the cubic-quintic NLS solitary wave as a particular example of solution possessing the space-time parity symmetry P . In section 6.1 we present a general computation for any anomaly β (n) r of the deformed NLS (2.2). We will show that in the co-moving frame of the solitary wave the anomalies β (n) r vanish. Once the anomaly is computed in the co-moving frame one performs a Galilean boost in order to compute for travelling waves. So, the computation is performed for a general deformed NLS required that it possesses a solitary wave solution. The cubic-quintic explicit solution is merely used to motivate the discussion since it possesses the symmetry P . In section 6.2 it is computed the particular anomaly β (4) r , through another method, for solitary waves of deformed NLS (2.2). In this case the computation process considers the relationship between the solitary wave components R and ϕ arising from the equations of motion in appendix A.
We will consider translationally invariant solutions of (2.2) of the form satisfying the following non-vanishing boundary condition (nvbc) The nvbc (6.2) is suitable for wave solutions in the form of localized 'dark' pulses created on the continuous wave background (cwb) at rest is the phase shift experienced by the cwb with intensity |ψ 0 | 2 , in the non-linear medium with intensity |ψ| 2 and varying according to the law − ∂V [|ψ| 2 ] ∂|ψ| 2 . Moreover, we will consider a family of equations such that small oscillations around |ψ| = |ψ 0 | of (2.2) defines a hyperbolic system with sound speed Then the dark solitons of the model (2.2) must propagate with the velocity v such that |v| < v s . As discussed above, some deformed NLS models support analytical dark solitary waves. In particular, the cubic-quintic non-integrable NLS system (2.3) possesses solitary waves and they are obtained considering the situations with > 0 and < 0 separately. So, the system (2.3) possesses dark soliton solutions of the form (6.1) provided by [19,27] The signs ± in Φ ± and Θ ± correspond to the case > 0 and < 0, respectively. The parameter ξ 1 is always positive and ξ 2 becomes positive (negative) for > 0 ( < 0). For > 0 they satisfy ξ 1 < |ψ 0 | 2 < ξ 2 , whereas for < 0 it holds ξ 2 < ξ 1 < |ψ 0 | 2 . Notice that k ± characterizes the inverse soliton width and √ ξ 1 is the minimum intensity (dip) of the dark soliton. The maximum intensity of the dark soliton approaches Φ 2 (±∞) = |ψ 0 |. Moreover, for v = 0 the parameter ξ 1 vanishes and the dark soliton becomes a black soliton.
The space-time parity transformation (5.22) can be performed defining In fact, under (5.22) the solution (6.5)-(6.6) transforms as Therefore, it can be verified that the squared absolute value R = |ψ| 2 = Φ 2 ± (z) and phase ϕ = 2(wt + δ) + 2Θ ± of ψ = √ Re i ϕ 2 satisfy P : R → R, ϕ → −ϕ + const. the anomaly β (n) does not get renormalized. In order to compute the general renormalized anomaly β (n) r associated to a solitary wave of the deformed NLS model it is convenient to consider the co-moving coordinate of the soliton. In this frame, the soliton is stationary. In this frame consider the cwb solution of the deformed NLS model as a stationary oscillatory phase as in (6.3). In addition, a stationary soliton can be transformed into a non-stationary form through the Galilean symmetry. In fact, assume that ψ(x, t) is a solution of the deformed NLS equation (2.2) and v a constant, then is also a solution. So, in the co-moving coordinate of the solitary wave one can write the solitary wave solutions Φ 2 ( x) and Θ( x) as dependents only on the variable x. As an example, we can rewrite the solitary wave in (6.5)-(6.6) as

JHEP03(2016)005
Then, the space-time transformation (3.32) applied to this form of solution furnishes the next symmetry transformation for R( x) = |ψ| 2 = Φ 2 ± ( x) and phase ϕ( So, for solutions satisfying this type of space-time symmetry one has that the factor α (3,−n) in the anomaly β n in (3.31) is even according to (3.38). The first expressions of α (3,−n) provided in the appendix B exhibit this property. Therefore, the associated charges Q (n) satisfy a mirror type symmetry (3.40) and an asymptotically conservation law (3.41). Moreover, in this case we can also explore the space-reflection symmetry of the anomalies which can be useful in order to estimate their values, by using general symmetry considerations. So, note that the time dependence of the solution ψ(x, t) in (6.1) for functions of the form Φ 2 ( x) and Θ( x) appears only in the cwb phase (6.3). In addition, the components α (3,−n) are polynomials in R( x) and ∂ x ϕ( x), and their x−derivatives (see the first five expressions in the appendix B). So, the space-time parity symmetry can be used in order to see the properties of these expressions. In fact, the parities of R and ∂ x ϕ of this special soliton solution can be obtained by taking t = 0 in this frame. Therefore, using the result (3.38) one has that the α (3,−n) are even, whereas X is odd under the above space-reflection, then one has that the anomalies β (n) vanish for this type of solitary wave. Therefore, in the co-moving frame one has dQ (n) dt = 0, n = 0, 1, 2, 3, . . . (6.18) Let us discuss this result for travelling waves. For travelling waves the anomaly associated to a solitary wave can be written as x ϕ(z), . . .]; z = x − vt. (6.19) In fact, the integrand of these expressions are some polynomials containing the functions R(z) and ∂ x ϕ(z), and their x−derivatives. Next, we can rewrite the integration (6.19) by making a coordinate transformation x → z, t → t. Furthermore, one can relabel z → x, and the outcome would be an equivalent integral expression to the one in the co-moving frame we have already analized above. Therefore, also for the travelling waves the anomalies vanish and the associated charges will be conserved.  in (4.33). However, we perform this computation using the relationships regarding this type of solutions presented in appendix A for general deformations of the NLS model of type (2.2). One

JHEP03(2016)005
proceeds by discarding some 'surface terms' and then rewriting the improper anomaly integral lim L→∞ [ L −L (. . .)dx] as a definite integral in the field space (with R as an integration variable) of type ± |ψ 0 | 2 R min (. . .)dR, where the limits of integration R min and |ψ 0 | 2 characterize the minimum intensity dip of the dark soliton and the cwb amplitude, respectively. Let us make a coordinate transformation x → z, t → t and rewrite (4.33) as Due to (A.1) one can write ϕ (z) = 2Θ (z) and considering (A.3) one notices that Θ (z) depends only on R. Moreover, using (A.5) one has that the term inside parenthesis can be written as . Therefore, all the terms in the integrand above, except the last one, can be rewritten in the form K[R]R , which can be written as a total derivative in z, then this 'surface term' does not contribute to the anomaly for the boundary condition (4.1) with k = 0. Regarding the last term, we can rewrite the integration in field space as follows

Simulations
Firstly, we have solved numerically the ordinary NLS model (5.1) without deformation ( = 0). We develop this procedure for a general initial condition with a cwb, for the one-dark and two-dark soliton solutions, respectively. Then, it is convenient to rewrite the dark soliton (5.10) as where α ≡ 2p R , v ≡ −2p I . Notice that the dark soliton (7.1) is characterized by the width 1/α and the velocity v. The background density becomes |ψ 0 | = (α 2 + v 2 ) 1/2 / 4|η|.

JHEP03(2016)005
In order to study the interaction of 2-dark solitons we will choose the initial data (t = 0) as (consider η = 1/4) The two solitons are initially centered at ±x 0 (x 0 > 0), the soliton centered initially at −x 0 (t = 0) moves to the right with velocity v 1 (v 1 > 0), whereas −v 2 (v 2 > 0) is the velocity of the soliton initially (t = 0) centered at x 0 and travels to the left, α > 0 is a constant characterizing the width of the solitons. In the numerical simulation we will consider initially well-separated solitons, i.e. the parameter 2x 0 is chosen to be several times the width of the solitons. The equation ( .1) is satisfied for each time step. In our numerical simulations we will use the so-called time-splitting cosine pseudo-spectral finite difference (TSCP) and the time-splitting finite difference with transformation (TSFD-T) methods [16,17] for k = 0 and k = 0, respectively. Our numerical simulations reproduce the main properties already known for dark soliton interactions in the integrable defocusing NLS model. Next, we summarize the main properties of the interaction of two dark solitons: a). The interaction is repulsive, and the initial solitons recover completely their forms and velocities after the collision.
b). For a given asymptotic value |ψ 0 | there exists a critical velocity v c = v s /2, where v s is the sound speed (5.14) which defines two types of solitons characterized by their behaviour under the collision: when 0 < v < v c one has the low speed solitons which are reflected by each other, and when v c < v < v s they are the high speed solitons which are transmitted through each other [28]. c). After collision the two solitons recover completely their properties, i.e. the NLS model is integrable, then all anomalies vanish (all charges are conserved).
The numerical method TSFD is specialized in order to compute the soliton type solutions of the deformed NLS model (2.2). The first deformed model (2.3) possesses a solitary wave presented in (6.5)-(6.6). Then, we can take two one dark solitons some distance apart as the initial condition for our numerical simulations. For the second deformed model (2.4) = 0, q ∈ Z Z + we do not know any analytical expression for its solitary wave; however one can generate numerically a solitary wave, as we will explain below. We will study the dark soliton collisions numerically for the both deformed NLS models and compare the outcomes with the properties summarized above for the integrable NLS model. The collision of dark solitons in the cubic-quintic NLS equation (2.3) can been simulated numerically by considering the initial condition ψ(x) defined as where two solitary wave solutions (6.5)-(6.6) of the cubic-quintic NLS model have been located some distance apart. The function ψ is continuous at x = 0, i.e. it has continuous absolute value |ψ| = Φ(x 0 ) and phase equal to Θ(x 0 ). The parameters {ξ 1 , ξ 2 , r, k ± } are defined in (6.5)-(6.6), so the ± signs correspond to the relevant sign of the deformation parameter . The two solitons are initially centered at ±x 0 (x 0 > 0), the soliton centered initially at −x 0 (t = 0) moves to the right with velocity v > 0, whereas the soliton initially (t = 0) centered at x 0 travels to the left with velocity (−v). Notice that the direction of motion of the soliton is related to the sign of the phase slope in (7.3). In addition, we will consider initially well-separated solitons, i.e. the parameter 2x 0 is chosen to be several times the width of the solitons (∼ 1 k ± ) and 2x 0 < 2L. So, the initial condition uses two of its static analytic solitary waves which are stitched together at the middle point, and then we allow the scattering of them, absorbing the radiation at the edges of the grid. It amounts to maintain the nvbc at the edges of the grid for each time step of the numerical simulation. For solitons colliding with different velocities (different amplitudes) we must modify conveniently this initial condition.
In the figures 2-10 the results for the simulations of soliton collisions in the cubicquintic NLS model (7.3) are presented. The relevant anomaly β  Since the two-soliton solutions admit in general the symmetry (3.33), as presented in section 3.1 through perturbation theory on the deformation parameter, we expect their anomalies to vary during collision and the charges to be asymptotically conserved. However, for a special soliton solutions with space-reflection symmetry, the anomaly β (4) belongs to the sequence of even order charges and it must vanish during the whole collision process. In fact, this behaviour is observed through our numerical simulations for various soliton JHEP03(2016)005    in (4.30). Therefore, the charges must be exactly conserved, within numerical accuracy. So, our numerical simulations suggest that the space-reflection P x symmetry is a sufficient but not a necessary condition for the vanishing of the anomaly β                       Notice that the charge Q (4) in (4.30)-(4.31) is conserved to within an accuracy of ≈ 10 −11 in the second modified NLS model even for deformation parameters of the order of | | ≈ 0.1, whereas in the first modification the CQNLS model this error is of the order of ≈ 10 −5 for | | ≈ 0.01. In our extensive simulations we noticed that in general the vanishing of the anomalies occur with small numerical errors provided that the parameter value is small and the solitons are fast. Then, we are quite confident of our simulations in the second modified NLS model in which the anomalies vanish within an error of less than 10 −10 for | | < 1.5 and n ≥ 2 in the whole process of the collisions of two fast and/or JHEP03(2016)005      slow solitons. The main conclusion that we can make from our simulations for several twosoliton collisions is that the charge Q (4) is not only asymptotically conserved, but exactly conserved during the whole collision process, within the numerical accuracy. The results of our simulations of the second deformation of the NLS model (2.4) and the fact that it possesses a submodel (with = 1, q ∈ Z Z + ) which passes the Painlevé test [22] suggest that it might be an exactly integrable model for some set of the parameters {|ψ 0 |, η, , q}. To the best of our knowledge its Lax pair formulation and analytic soliton solutions are not known in the literature and the search for them are worth attempting. However, the existence of two-solitons is not indicative of integrability. In fact, some nonlinear models present elastic scattering of two solitons but they do not posses three-soliton JHEP03(2016)005   solutions. Nevertheless, according to the ref. [29] the above deformed NLS model (2.4) may be regarded as a partially integrable model in the sense that it possesses two-soliton solutions, even though it has been obtained numerically. Moreover, the existence of threesoliton solutions in partially integrable models is quite demanding and it is very rare. In fact, in the Hirota's bilinear approach the existence of three-solitons is a strict integrability criterion [29]. So, it would be interesting to tackle this problem in the lines discussed above for some deformations of the NLS model, since the presence of three-soliton combinations of any type may indicate integrability of the model.

Discussions and some conclusions
In this paper we have studied some properties of the deformed non-linear Schrödinger model defined by the eq. of motion (2.2), with the potential involving a deformation parameter and the non-vanishing boundary condition (3.26). We have shown that such deformed NLS theories possess some sectors where the solitary wave solutions exhibit properties very similar to solitons in exactly integrable field theories. They possess an infinite number of exactly conserved quantities for one-soliton type solutions, i.e. for dark solitary waves travelling with a constant speed. For solutions possessing a special space-time parity sym-JHEP03(2016)005   metry, using analytical methods borrowed from integrable theories, it can be shown that those quantities are asymptotically conserved for two-soliton type solutions. This means that these quantities can vary in time during the collision process of two one-solitons but return, in the distant future (after the collision), to the values they had in the distant past (before the collision). Furthermore, for solutions possessing the space-time parity symmetry and in addition, a special space-reflection parity symmetry we can say even more. Using an algebraic technique which combines the space-reflection and an order two Z Z 2 automorphism of the affine sl(2) loop algebra, in the context of an anomalous zero-curvature (Lax) equation, where the Lax potentials live on the loop sl(2) affine Kac-Moody algebra, we have shown that those infinite number of charges split into two sets, one of them contains a sequence of even order exactly conserved charges and the other one comprises asymptotically conserved charges. We have verified this property, analytically and numerically, for the head-on collision of two solitons travelling in opposite directions with equal velocities.
The mechanism behind the exact conservation of the sequence of even order charges is not well understood yet. As in all of the examples where the asymptotically conserved charges have been observed so far, associated to a space-time parity symmetry [1][2][3][4], the two-soliton type solutions associated to a sequence of exactly conserved charges present JHEP03(2016)005   special properties under a space-reflection parity transformation. The only explanation we have found, so far, for such behaviour of the even order charges, is that those special solitonlike solutions transform in a special way under a space-reflection parity transformation for a fixed time, where the point in x−coordinate around which space is reversed depends upon the parameters of the particular solution under consideration. In order to compute the vanishing of the anomaly for a general solitary wave solution we have used the properties of the solutions under the Galilei transformation and spacetime transformations in the co-moving coordinate of the soliton to prove that the quantity α (3,−n) satisfies the property (3.38) under the parity transformation for t = 0. In addition, we have shown that by considering the zeroth order solution in the deformation parameter to be a two-soliton solution of the integrable NLS model, satisfying (5.41), one can always construct by a power series expansion on , a solution of the deformed model satisfying the parity property (5.41), which implies the exact conservation of the even order charges.   We have confirmed the predictions of our analytical calculations through an efficient and accurate numerical method, the so-called time-splitting cosine pseudo-spectral finite difference method, which is appropriate for the case in which the cwb is at rest. The numerical code was devised for two important deformations of the NLS. The first deformation is related to the cubic-quintic non-integrable NLS model and its numerical simulation considers as an initial condition two of its solitary waves located some distance apart. The second deformation considers a saturable type potential and its solitary wave is generated numerically by taking as a seed solution the analytic solitary wave of the cubic-quintic NLS model.

JHEP03(2016)005
For the above two deformations we have computed the first non-trivial anomaly β (4) of the Q (4) r charge's non-conservation law. We have verified that this anomaly vanishes, and consequently the exact conservation of the charge Q       Decoupling the last two equations one can write a differential equation for R as one can integrate once the eq. (A.5) to get