arxiv: v1 [cs.it] 1 Oct 2018

Transcription

1 The loglog growth of channel capacity for nondispersive nonlinear optical fiber channel in intermediate power range. Extension of the model. A. V. Reznichenko, 1,, A. I. Chernykh, 3,, S. V. Smirnov,, and I. S. Terekhov 1,, 1 Budker Institute of Nuclear Physics of Siberian Branch Russian Academy of Sciences, Novosibirsk, 39 Russia Novosibirsk State University, Novosibirsk, 39 Russia 3 Institute of Automation and Electrometry of the Siberian Branch of the Russian Academy of Sciences, Novosibirsk, 39 Russia arxiv:11.513v1 cs.it] 1 Oct 1 In our previous paper 1] we considered the optical channel modelled by the nonlinear Sсhrödinger equation with zero dispersion and additive Gaussian noise. We found per-sample channel capacity for this model. In the present paper we extend per-sample channel model by introducing the initial signal dependence on time and the output signal detection procedure. The proposed model is a closer approximation of the realistic communications link than the per-sample model where there is no dependence of the initial signal on time. For the proposed model we found the correlators of the output signal both analytically and numerically. Using these correlators we built the conditional probability density function. Then we calculated an entropy of the output signal, a conditional entropy, and the mutual information. Maximizing the mutual information we found the optimal input signal distribution, channel capacity, and their dependence on the shape of the initial signal in the time domain for the intermediate power range. PACS numbers: 9.7.-a, 5.1.Gg I. INTRODUCTION. Nonlinear communication channels have received a lot of attention in last twenty years due to the development of fiber optical communication systems. In these communication systems the Kerr nonlinearity in the optical fiber becomes important when one increases the power of transmitted signal. The problem of capacity finding was considered analytically and numerically in a series of papers, see e.g. 1 11] and references therein. In spite of a lot of publications this problem has not been solved for the case of arbitrary Kerr nonlinearity and the second dispersion parameter of an optical fiber. The nondispersive model is much simpler than the case with nonzero dispersion but it catches-up the main features connected with nonlinearity. Also this model is more convenient for the understanding of dependence of the capacity on the channel nonlinearity. So the analytical form of the conditional probability density function PY X], i.e., the probability density function PDF to receive the output signal Y if the input signal is X, for nondispersive per-sample channel was first obtained in Ref.1]. The upper bound for the capacity at very large input signal power for the model was obtained in Refs. 13, 1]. The capacity of the channel was found in Ref. 1] in the intermediate power range implying both large signal-to-noise ratio SNR and the condition for the next-to-leading corrections in the noise power to be small, see Eq. 3 in Ref. 1]. Electronic address: a.v.reznichenko@inp.nsk.su Electronic address:chernykh@iae.nsk.su Electronic address:smirnov@lab.nsu.ru Electronic address:i.s.terekhov@gmail.com The per-sample model assumes that the input signal does not depend on time. In realistic communication channel the transmitted signal does depend on time. In the recent paper 15] the influence of the receiver, signal, and noise bandwidth on the autocorrelation function and the capacity was discussed within the filter-and-sample model for the channel with zero dispersion. In our opinion, one of the important result of the paper 15] is understanding that the conditional PDF depends significantly on the properties of the receiver. In this paper we consider the nondispersive channel in the intermediate power range in the case where the initial signal depends on time and has the bandwidth much less than the noise bandwidth. We also introduce the detection procedure which takes into account the time resolution characteristics of the detector and we demonstrate the influence of the detector and the noise bandwidth on statistical properties of the channel. Therefore this paper is the generalization of the previous results of Refs.1, ] for the per-sample model to the time-dependant signal. The paper is organized in the following way. In Sec. II we present the model of the signal propagation, the input signal and the receiver model. In Sec. III we obtain the conditional probability density function for the introduced model. In Sec. IV we present numerical results for the correlators and compare these results with analytical ones. And in Sec. V we calculate the optimal input signal distribution and the channel capacity in the intermediate power range. In the Conclusion we discuss our results.

2 II. MODEL OF THE SIGNAL PROPAGATION AND DETECTION In our model the propagation of the signalψz,t is described by the stochastic nonlinear Sсhrödinger equation NLSE with zero dispersion: z ψ iγ ψ ψ = ηz,t, 1 where γ is the Kerr nonlinearity coefficient, the function ψz, t obeys the input and output conditions: ψz =,t = Xt and ψz = L,t = Yt, respectively, L is the length of the signal propagation, and ηz,t is an additive complex noise with zero mean ηz,t η =, and the correlation function in the frequency domain: W ηz,ω ηz,ω η = πqδω ω θ ω δz z, where the bar means complex conjugation, and Q is a power of the noise per unit length and per unit frequency, θω is the Heaviside theta-function, δω is the Dirac delta-function, W is the bandwidth of the noise. The noise ηz,ω is not white due to limited bandwidth. In the time domain this correlator has the form ηz,t ηz,t η = Q W π sinc W t t δz z.3 One can see that if the time difference t t = nπ/w then the correlator 3 is equal to zero, here n is integer. Thus we can solve equation 1 independently for parameters t j = j for different integer j, where = π/w is the time grid spacing. Therefore instead of the continuous time model 1 we will consider the following discrete model: z ψz,t j iγ ψz,t j ψz,t j = ηz,t j for any time moment t j. It means that we obtain the set of independent time channels since the noise in these moments is not correlated. We present the input and output conditions in the discrete form as well: ψz =,t j = Xt j and ψz = L,t j = Yt j. Note that the solutionφz,t j of the equation with zero noise which obeys the input condition Φz =,t j = Xt j has the form: Φz,t j = Xt j e iγz Xtj. 5 Below we assume that the frequency bandwidth W of the noise is much broader than the frequency bandwidth W of the input signal Xt and the frequency bandwidth W of the function Φz = L,t. In our model the input signal Xt has the form: Xt = N k= N C k ft k, where C k are complex random coefficients with some probability density function P X {C}], {C} = {C N,...,C N }; the pulse envelope ft is the real function which is normalized as dt f t = 1. The pulse envelope ft has the following properties: the overlapping of the functions ft k and ft m for k m is negligible: dtft kft m. It means that the function ft has almost the finite support /, /], and the input signal Xt is defined on the interval T = N + 1. Thus the frequency support of the function Xt is infinite. But we imply that W Xω dω W Xω dω, where Xω is the Fourier transformation of Xt. The last relation means that W 1. In our consideration the average input signal power P is fixed: N P = k= N d C k P X {C}] dt T Xt, 7 where d C k = drec k dimc k, and the input signal probability density function P X {C}] is normalized as follows: N k= N d C k P X {C}] = 1. Using the properties of the function ft k we can rewrite equation 7: where P = P m X C m] = d C m P m X C m] C m, 9 N k= N,k m d C k P X {C}], 1 and we imply that the distribution P m X C m] does not depend on m. Let us describe the output signal detection procedure. Our detector recovers the information which is carried by the coefficients {C k }. First, the detector receives the signal ψz = L,t j at the discrete time moments t j = j, here j = M,...,M 1, where M = T/ N. It means that the time resolution of the detector coincides with the time discretization. Since 1/ W our detector can completely recover the input signal in noiseless case. Second, the detector removes the nonlinear phase to obtain the recovered input signal Xt in the following form Xt j = ψz = L,t j e iγl ψz=l,tj. 11 And finally, using Xt detector recovers the coefficients

3 3 C k by projecting on the basis functions ft k : C k = 1 M 1 j= M dtft k Xt ft j k Xt j. 1 One can check that in the case of zero noise Xt = Xt and C k = C k. III. STATISTICS OF C k In the previous paper 1] we obtained the conditional probability function PY X] for the case where input X and output Y signals do not depend on time persample conditional PDF. In the previous section we extend our model 1] by including detector procedure and time dependence of the input signal Xt. Our goal is to obtain conditional probability function P{ C} {C}], i.e., the probability to detect the set of coefficients { C} if the transmitted set is {C}. Using the function P{ C} {C}] we can calculate the probability density function P out { C}] as N P out { C}] = k= N d C k P{ C} {C}]P X {C}]. 13 Since the propagation of the signal in the different time moments t j is independent, and noise is not correlated, the conditional probability function PYt Xt], i.e., the probability density to obtain the output signal Yt for the given input signal Xt, can be presented in the factorized form: PYt Xt] = M 1 j= M P j Y j X j ], 1 where X j = Xt j, Y j = Yt j, and P j Y j X j ] is persample conditional PDF obtained in Ref. 1]. The function P j Y j X j ] in the leading and next-to-leading order in parameter Q can be deduced from the results of Ref. 1], where we have to replace parameter Q by Q/ : { exp 1+µ j /3x j µ jx j y j +yj QL1+µ j P j Y j X j ] = /3 πql 1+µ j /3 5 µ j /3y j } 1 µ j /ρ j 151+µ µ j 15+µ j /3 j x j µ j 3 {µ j µ 135QLρ j 1+µ j /3 j +15µ j x +5 3 j 3µ + j +55µ j 9 x j y j + µ j µ j +117µ j 5 x j y j 3 5µ j +33µ j +3 y 3 j}, 15 here ρ j = X j, X j = ρ j e iφ j, µj = γlρ j, and x j = Re Y j e iφ j iµ j ρj ], yj = Im Y j e iφ j iµ j ρj ]. The expression 15 was obtained in Ref. 1] on the condition that the average input signal power P lies in the intermediate power range: QL P / QL 3 γ, 1 where P = π dρρ 3 Pρ], Pρ] is the distribution function of the quantity ρ, see Ref. 1]. Therefore, our consideration is restricted by the condition 1. The factorization of PYt Xt] in the form 1 means that there are M independent sub-channels. Note that, the signal Xt is completely defined by N+1 coefficients C k, i.e., there are only N +1 independent X j, but all M quantities Y j are independent. However, our detector reduces the function Yt to the set of N+1 coefficients{ C k } by the procedure 11 and 1. Therefore we have to reduce the function PYt Xt] to the function P{ C k } {C k }] by integrating over M N 1 redundant degrees of freedom. Using the conditional PDF PYt Xt] in the form 1 one can calculate all correlators of the coefficients C k : C k1, C k1 Ck, C k1... C kn. Here C k1... C kn = M 1 d Y j PYt Xt] C k1... C kn, 17 j= M where d Y j = drey j dimy j, and C k is defined in equation 1, and in the discrete form it reads: C k = M 1 j= M ft j k Y j e iγl Yj. 1

4 To recover the function P{ C k } {C k }] in the leading approximation in parameter Q it is necessary to know only three correlators: C k, C k Cm, C k Cm. After substitution of Eqs. 1, 15, and 1 to Eq. 17 and performing the integration we obtain in the leading order in the noise parameter Q: C k = C k ic kql γ Cm C m Cn C n iδ m,n C mql γ 1 iγl C k n, 19 3 = n in 3 γl C m, Note that for the first correlator C k C k is proportional to QL/ = QLW /π, i.e., it is proportional to the total noise power. Whereas the correlators and 1 are proportional to QL/ and do not depend on the discretization parameter only in leading order in parameter Q and depend on the parameter in higher order corrections in parameter Q, see Appendix A. Using the correlators 19 1 we obtain the conditional PDF P C C] in the leading order in parameter Q: Cm C m Cn C n = δ m,n QL 1+ n 3 γ L C m, 1 N P C C] = P m C m C m ], 3 m= N where δ m,n is Kronecker symbol and n s = T / / dt f s t. where P m C m C m ] πql 1+ξ µ m/3 exp 1+n µ m/3 x m +x m y m µ m n +y m QL1+ξ µ m /3 ]. Here we have introduced the notations: x m = Re y m = Im e iφm { Cm C m + ic m γql 1 iγl C m n } ], 5 3 e iφm { Cm C m + ic m γql 1 iγl C m n } ], 3 φ m = argc m, µ m = γl C m, 7 ξ = n 3n. The parameter ξ obeys inequality ξ > n > due to Cauchy-Schwarz-Buniakowski inequality. Note that the function P C C] has the factorized form 3 only in the leading approximation in the parameter Q. The Eq. 3 means that we have N + 1 independent information channels, and the channel corresponding to the time slotmis described by the functionp m C m C m ]. The function P m C m C m ] obeys the normalization condition d Cm P m C m C m ] = 1. 9 Since there are N + 1 independent channels, we can choose the input signal distribution P X {C m }] in the factorized form: P X {C}] = N k= N P k X C k], 3 and we can consider only one channel, say m-th channel. For this channel we can calculate the probability distribution function of the coefficients C m : P m out C m ] = d C m P m C m C m ]P m X C m]. 31 We imply that the function P m X C m] is a smooth function that changes on a scale C m P which is much

5 5 greater than QL/ : P QL/ QL/. 3 In other words, the signal power is much greater than the noise power in the channel. The variation scale of the function P m C m C m ] in the variable C m is of order of QL/ therefore we can use Laplace s method 1] for the calculation of the integral 31. Performing the integration in the leading order in parameter Q we obtain P m out C m ] P m X C m ], 33 for details see Appendix C in Ref. 1]. The result 33 implicates that the statistics of the coefficients C m coincides with the statistics of the coefficients C m. IV. NUMERICAL CALCULATIONS OF THE CORRELATORS In order to verify analytical results we performed numerical simulations of pulse propagation through nonlinear nondispersive optical fiber and calculated correlators 19,, and 1. For these purposes we solve numerically Eq. 1 for fixed input signal Xt and for different realizations of the noise ηz, t. Then we numerically perform the detection procedure described by Eqs. 11, 1. Finally, we average the coefficients C k and their quadratic combinations over noise realizations. In our simulations we use two numerical methods of the solution of Eq. 1: the split-step Fourier method and Runge-Kutta method of the fourth order. The results are presented in the following subsections. We have checked that the numerical results do not depend on the numerical method and these results are consistent with analytical ones for different realizations of the form ft of the input pulse. For numerical simulation we choose the following realistic channel parameters. The duration of one pulse is = 1 1 sec; fiber length is equal to L = km; Kerr nonlinearity parameter is γ = 1.5 km W 1. A. Split-step Fourier method Equation 1 was integrated numerically over z from up to communication line length L using split-step Fourier method 17], 1]: ψz +h,t = ψz,texp iγ ψz,t h + ˆF {δq h }, 3 where ψz,t stands for numerical solution of 1, h is a step size of z-mesh, ˆF denotes discrete inverse Fourier transform. The quantity δq h stands for the noise addition per step h which is made in frequency domain according to hq ψz,ω j ψz,ω j + T ηx +iη Y, 35 where j =,...,M 1 stands for index of ω-mesh, M is the number of t- and ω-mesh points, T is the total width of t-mesh, we choose T =, see Eq. 3 below; η X and η Y are independent standard Gauss random numbers with zero mean and σ = 1, additive noise level is Q = 1 1 W/km Hz. The input signal for z = has the form ψz =,t = Xt = C k ft k, 3 k=1 here we use the pulse envelope of the Gaussian form: ft = exp t, 37 T 1 π T 1 where T 1 = /1 = 1 11 sec stands for the characteristic time scale of the function ft. Pulse intersection is negligible. For such pulses coefficients n s defined in Eq. are n = T/T1 π 3.99, n = T/T1 π 1.3, 3 n = T/T13 π 9.79, ξ 5.. π In the numerical simulation we vary the average power 1 k=1 C k of the input signal from.177 mw up to.3 mw. It corresponds to the variation of the peak power C k f from.1 mw up to 5 mw. Simulations are performed for different t-meshes different grid spacing, i.e., for different noise bandwidths and fixed noise parameter Q. These meshes differ from each other by time grid spacing = T/M: 1 = sec, = sec and 3 = sec. These grid spacings determine the widths of conjugated ω-meshes: 1/ 1 = 1. THz, 1/ = 5.1 THz and 1/ 3 =.5 THz. For each average power of the signal and each mesh step size we simulate propagation of the signal for different realizations of the noise and then average obtained results for correlators over realizations. The total number of noise realizations for fixed Xt, see Eq. 3, is determined by the necessary statistic relative error and is chosen as This number of the realizations corresponds to the statistic relative error for correlators and 1 on the level of.% since the total number of pulses is 5. 1 = We performed simulations on z-meshes with different number of points 1,,, and checked out that the results do not depend on step size h. In Figs. 1 5 the numerical and analytical results for correlators 19 1 are presented for different time grid spacing as a function of input signal power. In Fig. 5 the results are presented for the grid spacings 1 and 3 because the results for 1 and almost coincide. One can see that numerical and analytical results are in a good agreement up to 3 mw at least. However the difference between numerical and analytical results for the smallest time grid spacing 1 is maximal. Decreasing of the parameter means the increasing of the spectral bandwidth of the noise. This increasing results in

6 1 3 Re Ck Ck /Ck C k, mw] Re Cm Cm C m, mw] Figure 1: The real part of the relative difference of the coefficient C k and the correlator 19 in units 1 3 as a function of input signal power C k for ft from Eq. 37. The noise power parameter is Q = 1 1 W/km Hz. Dashed doted, dashed, and solid lines correspond to analytic representation 19 for time grid spacings 1,, 3, respectively. Circles, squares, and diamonds correspond to numerical results for time grid spacings 1,, 3, respectively. Figure 3: The real part of the correlator multiplied by 1 as a function of input signal power C m for ft from Eq. 37. The noise power parameter is Q = 1 1 W/km Hz. Dashed doted, dashed, and solid lines correspond to analytic representation with NLOcorrections A1 for time grid spacings 1,, 3, respectively. Circles, squares, and diamonds correspond to numerical results for time grid spacings 1,, 3, respectively. /Ck Ck Ck 1 3 Im C k, mw] Figure : The imaginary part of the relative difference of the coefficient C k and the correlator 19 in units 1 3 as a function of input signal power C k for ft from Eq. 37. The noise power parameter is Q = 1 1 W/km Hz. Dashed doted, dashed, and solid lines correspond to analytic representation 19 for time grid spacings 1,, 3, respectively. Circles, squares, and diamonds correspond to numerical results for time grid spacings 1,, 3, respectively. 1Im Cm Cm C m, mw] Figure : The imaginary part of the correlator multiplied by 1 as a function of input signal power C m for ft from Eq. 37. The noise power parameter is Q = 1 1 W/km Hz. Dashed doted, dashed, and solid lines correspond to analytic representation with NLOcorrections A1 for time grid spacings 1,, 3, respectively. Circles, squares, and diamonds correspond to numerical results for time grid spacings 1,, 3, respectively. the growth of the total noise power received by detector. Note that the analytical expressions for correlators were obtained using the conditional PDF PYt Xt] in the form 15. This form was derived in the approximation of large signal-to-noise ratio SNR = P /QL. Decreasing parameter we diminish the parameter SNR and, as a consequence, the accuracy of our approximation. The difference between numerical and analytical results can be explained by taking into account the nextto-leading order NLO corrections in noise power parameter QL/. Analytical results in Figs. 3 5 are shown with taking into account both leading order results 1 and NLO corrections presented in Appendix A, see Eqs. A1 and A3. B. Runge-Kutta method For the equation 1 the time t is the incoming parameter. Thus the simulation consists in the solution of ordinary differential equation with various initial conditions determined by the real pulse shape ft, the amplitude C m, and independent random noise functions ηz,t. In the second method we used pulse envelopes of the form f n t = A n cos n πt/ 3 for n =,, t /, /]: A = 3 and A = We choose the time discretization parameter = /. The random noise was realized as the tele-

7 7 Cm Cm C m, mw] 1 3 Re Ck Ck /Ck C k, mw] Figure 5: The correlator 1 as a function of input signal power C m for ft from Eq. 37. The noise power parameter is Q = 1 1 W/km Hz. Dashed doted, dashed, and solid lines correspond to analytic representation 1 with NLO-corrections A3 for time grid spacings 1, 3, respectively. Circles and squares correspond to numerical results for time grid spacings 1, 3, respectively. Figure : The real part of the relative difference of the coefficient C k and the correlator 19 in units 1 3 as a function of input signal power C k for f t, see black solid line, and for f t, see black dashed line. The noise power parameter is Q = W/km Hz. Circles and rectangles correspond to numerical results with statistic error on the level of three standard deviations for the functions f and f, respectively. graph process with the step of the length z = 1 L and of the random height with zero average and with the dispersion σ =.3 1 W/km both for real and imaginary parts. The noise power parameter reads as Q = σ z W/km Hz and it is almost six times greater than that in the previous method. We independently control this parameter Q by using the leading order contribution to the correlator 1 numerically simulated for γ =. The noise η is constant within the step. Within the step the equation 1 was solved by the Runge-Kutta method of the fourth order with the step h = z /5. The recovered input signal Xtj was calculated using Eq. 11 at the equidistant points t j. The coefficients C k were calculated using Eq. 1. The average 19 and correlators,1 were calculated over 13 values of C k which were found for various noise realizations. To control the accuracy of the method we solved the equation 1 with zero noise from z = to z = L with the step h and then we performed the backward propagation from z = L to z = with the found solution as the initial condition. In the procedure the input signal was recovered with the relative precision equal to 1. The analytical results in comparison with the numerical results are presented in Figs. 1 for different pulse shapes and average power. The numerical results are presented with statistic errors on the level of three standard deviations. One can see that numerical and analytical results are in a good agreement as well. /Ck Ck Ck 1 3 Im C k, mw] Figure 7: The imaginary part of the relative difference of the coefficient C k and the correlator 19 in units 1 3 as a function of input signal power C k for f t and f t, see black solid line. The noise power parameter is Q = W/km Hz. Circles and rectangles correspond to numerical results with statistic error on the level of three standard deviations for the functions f and f, respectively. V. ENTROPIES AND MUTUAL INFORMATION Now we proceed to the calculation of the output signal entropy H C m ] = conditional entropy H C m C m ] = d Cm P m out C m ]logp m out C m ], 39 d Cm d C m P m C m C m ] P m X C m]logp m C m C m ],

8 Re Cm Cm C m, mw] Cm Cm C m, mw] Figure : The real part of the correlator multiplied by 1 as a function of input signal power C m for f t, see black solid line, and for f t, see black dashed line. Solid and dashed lines correspond to the real part of leading order contribution with the next-to-leading order corrections, see Eq. A1. The noise power parameter is Q = W/km Hz. Circles and rectangles correspond to numerical results with statistic error on the level of three standard deviations for the functions f and f, respectively. 1Im Cm Cm C m, mw] Figure 9: The imaginary part of the correlator multiplied by 1 as a function of input signal power C m for f t, see black solid line, and for f t, see black dashed line. Solid and dashed lines correspond to the imaginary part of leading order contribution with the next-to-leading order corrections, see Eq. A1. The noise power parameter is Q = W/km Hz. Circles and rectangles correspond to numerical results with statistic error on the level of three standard deviations for the functions f and f, respectively. and the mutual information I P m X = H C m ] H C m C m ]. 1 Our calculations of the entropies 39,, and the mutual information 1 are similar to calculations of the entropies and the mutual information for per-sample channel, see Sec.III and Sec.IV of Ref. 1]. Therefore we will not repeat the similar calculations here and present only Figure 1: Correlator 1 as a function of input signal power C m for f t, see black solid line, and for f t, see black dashed line. Solid and dashed lines correspond to leading order contribution 1 with the next-to-leading order corrections, see Eq. A3. The noise power parameter is Q = W/km Hz. Circles and rectangles correspond to numerical results with statistic error on the level of three standard deviations for the functions f and f, respectively. the final results: H C m ] = HC m ] = d C m P m X C m]logp m X C m], H C m C m ] = 1+log π QL ] + 1 d C m P m X C m]log 1+ξ γ L C m 3 ]. 3 To calculate the optimal input signal distribution P m opt C m] we calculate the mutual information substituting Eqs. and 3 to Eq. 1 then we variate the mutual information over P m X C m] with taking into account the normalization condition and the fixed average power 9. Assuming the variation of the mutual information to be zero, we obtain the equation for the optimal input signal distribution P m opt C m ]. We solve the equation and obtain for details of the similar calculations for per-sample channel see the Sec.III of Ref. 1]: opt C m] = N 1+ξ γ L C m /3, P m e λ Cm where parameters N = N P,ξγ and λ = λ P,ξγ are functions of the power P and modified nonlinearity parameter ξγ by virtue of the relations compare with Eqs. and 7 of Ref. 1]: d C m P m opt C m] = dρ πn ρe λρ = 1, 5 1+ξ γ L ρ /3

9 9 Capacity nat/symb.] PmW] Figure 11: Shannon capacity and the mutual information for the parameters Q = 1 1 W/km Hz; L = I P m opt km; γ = 1.5 km W 1 ; = 1 1 sec, and for the Gaussian shape 37 of ft. The black dotted line corresponds to the Shannon limit log, the black solid line P QL corresponds to I m P, see Eq. 7, the black dashed dotted opt line corresponds to the asymptotics 9 for large γlp. P = d C m P m opt C m ] C m = dρ πn ρ 3 e λρ 1+ξ γ L ρ /3. The capacity of one channel m, i.e., the mutual information calculated using the optimal input signal distribution reads C = I P m opt PT = log +Pλ logpn ].7 πeql One can see that the first term in the right-hand side of Eq. 7 corresponds to the Shannon s result 19] for the linear channel at large signal-to-noise ratio, the second and third terms are related with the nonlinearity impact. The result 7 is similar to that obtained for the persample model in Ref. 1] but with modification of the Kerr nonlinearity parameter γ for per-sample model to parameter ξγ for the present model, where ξ = n 3n. There is no simple analytical form for N and λ, see the Secs. III and IV of Ref. 1], therefore we present below the analytical results for the asymptotics of the mutual information for small and large dimensionless nonlinearity parameter ξγlp and the numerical calculations in Fig. 11. Performing the substitution γ ξγ in the results of the Sec.III and Sec IV of Ref. 1] we arrive at following asymptotics of the mutual information for small and large dimensionless nonlinearity parameter γlp: I P opt X X] log PT ξ γ L P, QL 3 for ξγlp 1, and I P opt X X] = loglog 1 log BξγLP/ 3 BξγLP/ 3 log QL ξγe/ 3 + log log BξγLP/ 3 +1 loglog BξγLP/ 3 ] log BξγLP/ 3, 9 for logξγlp 1 and P /QL 3 ξ γ. Here B = e γe, γ E.577 is the Euler constant. Note that the asymptotics 9 is obtained with accuracy 1/log ξγlp, see the Sec. IV of Ref. 1]. VI. CONCLUSION In the present paper we use results obtained in Ref. 1] for per-sample model to calculate the informational characteristics of the channel where the input signal Xt depends on time, see Eq.. For this channel the information is carried by coefficients C k. In the process of the signal propagation the input signal is transformed by the Kerr nonlinearity and the noise in the channel. To recover the transmitted information we introduce the detection procedure which removes the nonlinearity effects, see Eq. 11, and then projects Xt on the basis functions, see Eq. 1, to obtain the coefficients C k. Using the conditional probability density function for persample model obtained in Ref. 1] we calculate the correlators of the coefficients C k, see Eqs We demonstrate that these correlators depend on the noise bandwidth parameter. We also perform the numerical calculations of these correlators using two different methods and show that the numerical and analytical results are in agreement. Using obtained results for correlators we find the conditional probability density function P{ C k } {C k }] in the leading and next-to-leading orders in parameter QL/ P. Then we calculate the informational entropies and the mutual information for the channel in leading order in the parameterql/ P. We perform variation of the mutual information over the input signal distribution function and obtain the optimal input signal distribution function which maximizes the mutual information. We calculate the channel capacity in the leading order in parameter QL/ P and demonstrate that the capacity depends on the pulse envelope through one parameter ξ, see Eq.. The capacity grows as log log P for sufficiently large average power P: ξγql 1 P /QL 3 ξ γ. Note that the same asymptotics was obtained for per-sample model, therefore taking into account the time dependance of the pulse envelope does not change the asymptotics behavior and modifies only the nonlinearity parameter γ to ξγ.

10 1 Acknowledgments Ref. ]: Acknowledgment All authors would like to thank the Russian Science Foundation RSF, grant No Also A.V. would like to thank the Russian Foundation for Basic Research RFBR, grant No Also I.S. would like to thank the Russian Science Foundation RSF, grant No , and Ministry of Education and Science of the Russian Federation 1.Y Appendix A: Correlators and 1 with NLO corrections Let us present the correlator with next-to-leading NLO corrections in the noise power. Cm C m Cn C n = δ m,n QL Cm C m Cm γ { C m Cm 1+ n 9 γ L C m +i n 3 γl C m } = QL γ δ m,n Cm n ] 3 γl C m in + QL γ T C m 9n + n 3 γ L C m + ] i 5n 15 γl C m. A1 Here we have used the relation Cm C m Cm C m =, Cm C m Cm C m Cm C m A the result 19 for Cm C m and the calculation of Cm C m Cm C m on the base of next-toleading order result for PY X] in Ref. ]. In a similar manner it is easy to calculate the following corrections to correlator 1 from the results obtained in 1Im Cm Cm C m, mw] Cm C m Cn C n = δ m,n Cm C m Cm C m QL { } γ C m 1+ n 9 γ L C m QL δ m,n 1+ n ] 3 γ L C m + QL γ T C m n n 9 γ L C m = ].A3 Note that these NLO results A1 and A3 contain the time discretization parameter related with the noise bandwidth W = π/. The relative importance of the NLO corrections in correlators A1 and A3 is governed by the dimensionless parameter QL γl γlp, i.e., it increases linearly for large and increasing P. To demonstrate the importance of these corrections for our numerical results we present the Fig. 1 where for the noise power parameter Q = W/km Hz the imaginary part of the leading order contribution and the next-to-leading order corrections A1 are presented together with the numerical results Runge-Kutta method for the envelope form f t = 3 cos πt/. One can see that our calculations, i.e., Eq. and formulae of the Sec. V based on the leading order results 19 1 are in a good agreement with the numerical calculations up to the average power of order of mw for given noise and channel parameters. Figure 1: The imaginary part of the correlator multiplied by 1 as a function of input signal power C m for f t = 3 cos πt/ in the leading order, see black dashed dotted line, and with the next-to-leading order cor-

11 11 1] I. S. Terekhov, A. V. Reznichenko, Ya. A. Kharkov, and S. K. Turitsyn, Phys. Rev. E 95, ] A. A. Panarin, A. V. Reznichenko, I. S. Terekhov, Phys. Rev. E 95, ] P. P. Mitra and J. B. Stark, Nature 11, ] E. E. Narimanov and P. Mitra, J. Lightwave Technol., 53. 5] J. M. Kahn and K.-P. Ho, IEEE. J. Sel. Topics Quant. Electron. 1, 59. ] R.-J. Essiambre, G. J. Foschini, G. Kramer, and P. J. Winzer, Phys. Rev. Lett. 11, ] R.-J. Essiambre, G. Kramer, P. J. Winzer, G. J. Foschini, and B. Goebel, J. of Lightwave Technol., 1. ] R. Killey and C. Behrens, J. Mod. Opt. 5, ] E. Agrell, A. Alvarado, G. Durisi, and M. Karlsson, IEEE/OSA J. Lightwave Technol. 3, 1. 1] M. A. Sorokina and S. K. Turitsyn, Nat. Comm. 5, ] I. S. Terekhov, A. V. Reznichenko, and S. K. Turitsyn, Phys. Rev. E, 9, ] A. Mecozzi, J. Lightwave Technol. 1, ] K.S. Turitsyn, S.A. Derevyanko, I.V. Yurkevich, and S.K. Turitsyn, Phys. Rev. Lett. 91, ] M. I. Yousefi and F. R. Kschischang, IEEE Trans. on Information Theory 57, ] G. Kramer, Submitted to the IEEE Transactions on Information Theory, 1 arxiv:175.5v. 1] M. A. Lavrentiev and B. V. Shabat, Method of Complex Function Theory. Nauka, Moscow, 197, in Russian; M. Lavrentiev and B. Chabot, Methodes de la Theorie des fonctions d une variable complexe Mir, Moscou, 1977, in French. 17] G. P. Agrawal, Nonlinear Fiber Optics Academic press, th edition, 7. 1] R. H. Hardin and F. D. Tappert, SIAM Rev. Chronicle, 15, ] C. Shannon, Bell Syst. Tech. J. 7, ; 7, 3 19.