This letter deals with the Darboux representation of multisoliton solutions of the nonlinear Schrödinger equation. As carriers of information, these multisolitonic signals offer a promising solution to the problem of nonlinear signal distortions in fiber optic channels 
. In any nonlinear Fourier transform (NFT) based transmission methodology seeking to modulate the discrete spectrum of the multisolitons, the unbounded support of such signals (as well as its Fourier spectrum) presents some challenges in achieving the best possible spectral efficiency[2, 3] which forms part of the motivation for this work.
The Darboux transformation has proven to be an extremely powerful tool in handling the discrete part of the nonlinear Fourier spectrum. The rational structure of the associated Darboux matrix was recently exploited to obtain fast inverse NFT algorithms in [4, 5]. In the particular case of -soliton solutions, the rational structure of the Darboux matrix facilitates the exact solution of the Zakharov-Shabat problem  for (doubly-) truncated version of the signal via the solution of an associated Riemann-Hilbert problem . In [4, 7], an exact method for computing the energy content of the “tails” of -soliton solutions was reported which was again based on the Darboux transformation. This method was further used in  to establish the sufficient conditions for either one-sided or compact support of the signals resulting from “addition” of boundstates. Given that the complexity of computing the Darboux matrix coefficients is
, these methods turn out to be extremely efficient eliminating any need for heuristic approaches based on the asymptotic expansions (with respect to the windowing parameter, say,). The present work, therefore, tries to further reinforce the idea that the Darboux representation can potentially facilitate a number of design and signal processing aspects of -soliton solutions. For the general case, when the reflection coefficient is bandlimited, the work presented in  may allow us to compute the Jost solutions of the seed potential with extremely high accuracy.
The present work is also motivated by the fact that the recent attempts [2, 3, 10]111The readers are warned that the representation of the norming constants used in these papers do not follow the standard convention and are incorrect in many cases. The relationship , where in the tuple comprising the discrete eigenvalue and the corresponding norming constant, does not hold when does not have an analytic continuation in . In fact, from  it is known that the remains invariant when boundstates are added to an arbitrary profile so that even if exists. towards optimizing the generated multisolitonic signals are either based on brute-force methods or asymptotic expansions with respect to the windowing parameter. These methods have serious drawbacks either because they do not scale well in complexity when the number of boundstates are only moderately high or because they are not reliable in the absence of a prior knowledge of the goodness of the approximations made. In this letter, we present some new identities that can potentially make the optimization problem amenable to some of the powerful optimization procedures available in the literature (see  and the references therein) at the same time completely circumventing the need for any heuristics. Given that the independent variables (i.e. discrete eigenvalues and the norming constants) in the optimization procedure are complex in nature, the framework based on Wirtinger calculus presented in  appears to be more convenient.
Ii Temporal Width
The temporal width a -soliton solution can be defined via the -norm of the profile which is also related to the energy of the pulse. Let the energy content of the “tails” of the profile, denoted by , be defined by
characterizes the total energy in the tails which is a fraction of the total energy of a -soliton solution. The total energy is given by . Now, if the tolerance for the fraction of total energy in the tails is , then effective the support, , of the profile must be chosen such that . Therefore, the effective temporal width of the -soliton solution can then be defined as which is a function of .
For the specific case of -soliton solutions, an exact recipe for computing was reported in [4, 7] which we summarize briefly as follows: Let be the matrix form of the Jost solutions. Then, from the standard theory of scattering transforms , it is known that, for ,
Now, the -soliton potentials along with their Jost solutions can be computed quite easily using the Darboux transformation (DT) [4, 7]. Let be the discrete spectrum of a -soliton potential. The seed solution here corresponds to the null potential; therefore, . The augmented matrix Jost solution can be obtained from the seed solution using the Darboux matrix as for . The Darboux transformation can be implemented as a recursive scheme. Let us define the successive discrete spectra such that for where are distinct elements of . The Darboux matrix of degree can be factorized into Darboux matrices of degree one as
where are the successive Darboux matrices of degree one with the convention that is the bound state being added to the seed potential whose discrete spectra is . Note that the Darboux matrices of degree one can be stated as
for and the successive Jost solutions, , needed in this ratio are computed as
The potential is given by
so that with . Next, our objective is to compute the derivatives of with respect to the discrete spectra of -soliton in order to facilitate optimization procedures that are based on gradients. In the following, we use the notation
For real-valued function , it is known that . In Wirtinger calculus, the complex gradient is defined as
Here, we have described the Darboux transformation only for the case of -soliton potentials; however, the recipe can be easily adapted to the case of arbitrary seed potentials. Note that this would require explicit knowledge of and .
Ii-a Derivatives with respect to norming constants
Note that and are independent of ; therefore,
By direct calculation, we have
where we have used to the Wronskian relation
Using the identity (12), it is straightforward to work out:
Note that is the last eigenvalue to be added using the DT iterations. Given that there is no restriction on the order in which the eigenvalues can be added, we can always choose to be added last. This would determine using DT iterations for arbitrary . Thus, the complexity of computing derivatives works out to be .
Before we conclude this discussion, let us examine the case of multisoliton solutions when is large. In this limit, we have
Thus, the stationary condition translates into . Therefore, asymptotically, minimizes the energy in the tails. This result can be easily verified from (9) which in the limit of large gives
Ii-B Derivatives with respect to discrete eigenvalues
Let and define so that
The ratio can also be computed in terms of the modified Jost solutions on account of the fact that falls out of the equation while taking the ratio:
This gives us the opportunity to compute the derivatives with respect to recursively:
Using the notation for the Wronskian of scalar functions, let us introduce
By direct calculation, we have
Note that and are independent of ; therefore, using the above identity, it is straightforward to obtain
Following as in the case of norming constants, we can always choose to be added last so that can be determined using DT iterations for arbitrary . Thus, the complexity of computing derivatives again works out to be .
Iii Spectral Width
Consider the Fourier spectrum of the multisoliton potential denoted by . Let us observe that the following quantities can be expressed entirely in terms of the discrete eigenvalues:
with . These quantities do not evolve as the pulse propagates along the fiber. From 
, the varianceis given by
This quantity characterizes the width of the Fourier spectrum. The biquadratic integral above cannot be computed exactly in general, however, can be computed in a straightforward manner: From (8), we have , we have which yields where (correcting a typographical error in )
Note that this inequality holds irrespective of how the pulse evolves as it propagates along the fiber.
Iv-a One-sided effective support
Let us consider the case where we want to introduce a boundstate with eigenvalue to any arbitrary profile such that the energy content of the tail is . The problem is to determine the norming constant which minimizes . To this end, setting , we have
It is easy to verify that the first choice corresponds to maximum which leaves us with so that , i.e., no part of the soliton’s energy goes into the tail . By a recursive argument, the conclusion holds for any number of boundstates provided .
Iv-B Adding a boundstate to a symmetric profile
Let us consider the case where we want to introduce a boundstate with eigenvalue to any arbitrary seed profile. The energy content of the tail of the seed profile is . The problem is to determine the norming constant which minimizes . To this end, setting , we have
For the sake of simplicity, we assume that the seed profile is symmetric. Further, we also assume that so that
with . In the following, we set . Then, using the symmetry relations, we obtain
Physically, is related to the translation of the profile; therefore, it is easy to conclude, for a symmetrical profile, that the extrema is obtained for . Putting
in (29) and using , we have . The solution of this equation works out to be
From the relations
In order to show that , consider
which shows that . Therefore, the extremal points for are given by
Iv-B1 Symmetric -soliton
The general result derived above can be applied to a symmetric -soliton potential. Let us assume that the seed potential is a symmetric -soliton potential with the discrete spectrum given by . The boundstate being introduced is characterized by . Expression for the Jost solutions can be obtained from (4) which leads to so that . It can be directly verified that corresponds to the minima of for all as follows: Given the symmetric nature of the profile, it suffices to find the minima of which reads as
It is straightforward to show that , and . Now, from
it follows that ; therefore, the minima of occurs at or .
-  S. K. Turitsyn, J. E. Prilepsky, S. T. Le, S. Wahls, L. L. Frumin, M. Kamalian, and S. A. Derevyanko, “Nonlinear Fourier transform for optical data processing and transmission: advances and perspectives,” Optica, vol. 4, no. 3, pp. 307–322, Mar 2017.
-  A. Span, V. Aref, H. Bülow, and S. Ten Brink, “On time-bandwidth product of multi-soliton pulses,” in Information Theory (ISIT), 2017 IEEE International Symposium on. IEEE, 2017, pp. 61–65.
-  ——, “Optimization of multi-soliton joint phase modulation for reducing the time-bandwidth product,” in Photonic Networks; 19th ITG-Symposium. VDE, 2018, pp. 1–8.
-  V. Vaibhav, “Fast inverse nonlinear Fourier transformation using exponential one-step methods: Darboux transformation,” Phys. Rev. E, vol. 96, p. 063302, 2017.
-  ——, “Fast inverse nonlinear Fourier transform,” Phys. Rev. E, vol. 98, p. 013304, 2018.
-  M. J. Ablowitz, D. J. Kaup, A. C. Newell, and H. Segur, “The inverse scattering transform - Fourier analysis for nonlinear problems,” Stud. Appl. Math., vol. 53, no. 4, pp. 249–315, 1974.
-  V. Vaibhav, “Exact solution of the Zakharov–Shabat scattering problem for doubly-truncated multisoliton potentials,” Commun. Nonlinear Sci. Numer. Simul., vol. 61, pp. 22–36, 2018.
-  ——, “Nonlinear Fourier transform of time-limited and one-sided signals,” J. Phys. A: Math. Theor., vol. 51, no. 42, p. 425201, 2018.
-  ——. (2018) Nonlinearly bandlimited signals. [Online]. Available: http://arxiv.org/abs/1811.06338
-  A. Span, V. Aref, H. Bülow, and S. Ten Brink. (2018) Time-bandwidth product perspective for multi-soliton phase modulation. [Online]. Available: http://arxiv.org/abs/1812.04443
-  J. Lin, “Evolution of the scattering data under the classical Darboux transform for SU(2) soliton systems,” Acta Mathematicae Applicatae Sinica, vol. 6, no. 4, pp. 308–316, 1990.
-  L. Sorber, M. van Barel, and L. De Lathauwer, “Unconstrained optimization of real functions in complex variables,” SIAM J. Optim., vol. 22, no. 3, pp. 879–898, 2012.