Optimum Power in a Multi-Span DWDM System Limited by Non-Linear Effects

Jose Alfredo Alvarez-Chavez^{1,2}^{*},
Rafael Sanchez-Lara^{3},
Jaime Rafael Ek-Ek^{2},
Herman Leonard Offerhaus^{1},
Manuel May-Alarcon^{3},
Victor Golikov^{3}

Show more

1. Introduction

Long distance high capacity optical transmission is achieved using optical fibre links in backbone of terrestrial and transoceanic communication systems. Linear and non-linear effects limit the capacity. Among the main non-linear effects are Four-Wave Mixing (FWM), Stimulated Raman Scattering (SRS), Stimulated Brillouin Scattering (SBS), Self-Phase Modulation (SPM), and Cross-Phase Modulation (XPM). Furthermore, Amplified Spontaneous Emission (ASE) noise from Erbium-Doped Fibres Amplifiers (EDFAs) is present. Dense Wavelength Division Multiplexing (DWDM) technology is used to make maximum use of the total bandwidth of single mode optical fibres [1] [2] . The effect of dispersion and nonlinearities limit the performance for transmission over long distances at data rates greater than 2.5 Gbits/s. EDFAs are used to compensate fibre losses and to increase the transmission distance, causing at the same time, an increase of ASE noise and nonlinearities. Different kinds of fibres have been designed to reduce the effect of dispersion and nonlinearities [3] [4] [5] . There are numerous studies on nonlinear impairments showing the importance and complexity of such phenomena for DWDM [6] - [11] .

In this paper, design parameters for multi-span DWDM system limited by FWM, SRS and ASE noise are proposed. The dispersion parameter limit is also calculated, and the optimum power transmission per channel is evaluated for different number of channels or inter-channel spacing. Additionally, our analysis shows different spectral regions where each nonlinearity is dominant, and how to find intersections where general spectral limits can be obtained, taking into account the combined effects of FWM, SRS, and ASE noise.

2. Model

The multi-spam DWDM model used is show in Figure 1.

Where N is the number of channels, L_{a} is the optical amplifier spacing, L is the system length, M is the number of segments, A_{1} and A_{m}_{-1} are the first and last EDFA respectively, f_{1} and f_{N} are the first and last frequency channel.

A total spectral bandwidth $W={f}_{N}-{f}_{1}=3.75$ THz (30 nm) within C-band is considered. The WDM density parameter ∆f is a function of the number of channels N and can be expressed as:

$\Delta f={f}_{i}-{f}_{i-1}=W/\left(N-1\right)$ (1)

The central frequency f is given by:

$f=\frac{{f}_{1}+{f}_{N}}{2}=193\text{\hspace{0.17em}}\text{THz}\left[1550\text{\hspace{0.17em}}\text{nm}\right]$ (2)

C bands central frequency is a well-known value. At 1550 nm its equivalent to 193 THz approximately (f = c/wavelength).

The FWM signal is located at a frequency ${f}_{n}={f}_{i}+{f}_{j}-{f}_{k}\left(i,j\ne k\right)$ and it is given by the mixing of three signal frequencies ${f}_{i},{f}_{j}$ and ${f}_{k}$ co-propagating through a single-mode fibre, which generates a new wave frequency (FWM signal) by means of nonlinear interaction. The power of the FWM signal at the ${f}_{n}$ frequency is given by [10] [12] :

${P}_{ijk}={k}^{2}\ast {P}_{i}\ast {P}_{j}\ast {P}_{k}\ast {\eta}_{ijk}\ast {d}_{ijk}^{2}{\left(\frac{M\ast {L}_{eff}}{{A}_{eff}}\right)}^{2}{\text{e}}^{-a\ast {L}_{a}}$ (3)

Figure 1. DWDM architecture.

$k=\frac{32\ast {\text{\pi}}^{3}\ast {\chi}^{\left(3\right)}}{{\eta}^{2}\ast c\ast \lambda}$ (4)

where,
$\lambda $ is wavelength, c is the light speed in vacuum,
$\eta $ is the core refractive index, a is the linear loss coefficient,
${P}_{i,j,k}$ are launched input power per channel,
${L}_{eff}=\left(1-\mathrm{exp}\left(-a\ast {L}_{a}\right)\right)/a$ is the fibre effective length,
${A}_{eff}$ is the effective area of the guided mode, d is the degeneracy factor (
$d=3$ for
$i=j$ ,
$d=6$ for
$i\ne j$ ),
${\chi}^{\left(3\right)}$ is the third-order nonlinear electric susceptibility (6 × 10^{−14} m^{3} W^{−1}) and
${\eta}_{ijk}$ is the FWM efficiency, given by:

${\eta}_{ijk}=\left(\frac{{a}^{2}}{{a}^{2}+\Delta {\beta}_{ijk}^{2}}\right)\left(1+\left(\frac{4{\text{e}}^{-a\ast {L}_{a}}}{{\left(1-{\text{e}}^{-a\ast {L}_{a}}\right)}^{2}}\right){\mathrm{sin}}^{2}\left(\frac{\Delta {\beta}_{ijk}\ast {L}_{a}}{2}\right)\right)$ (5)

where $\Delta {\beta}_{ijk}$ represents the phase mismatch which may be expressed in terms of signal frequency differences and fibre dispersion as:

$\Delta {\beta}_{ijk}=\left(\frac{2\ast \text{\pi}\ast {\lambda}^{2}}{c}\right)\left(\left|{f}_{i}-{f}_{k}\right|\left|{f}_{j}-{f}_{k}\right|\right)\left(D+\frac{\text{d}D}{\text{d}\lambda}\left(\frac{{\lambda}^{2}}{2\ast c}\right)\left(\left({f}_{i}-{f}_{k}\right)+\left({f}_{j}-{f}_{k}\right)\right)\right)$ (6)

here, D is the fibre chromatic dispersion coefficient which is wavelength dependent. In a multi-wavelength system with N wavelength channels of equal channel spacing, the total FWM power generated at the frequency f_{n} may be expressed as a summation:

${P}_{n}=\underset{{f}_{n}={f}_{i}+{f}_{j}-{f}_{k}}{{\displaystyle \sum}}{P}_{ijk}={k}^{2}\ast {P}_{s}\ast {\text{e}}^{-a\ast {L}_{a}}{\left(\frac{M\ast {L}_{eff}}{{A}_{eff}}\right)}^{2}\underset{{f}_{n}={f}_{i}+{f}_{j}-{f}_{k}}{{\displaystyle \sum}}{\eta}_{ijk}\ast {d}_{ijk}^{2}$ (7)

where the summation is made over all the relevant combinations which satisfy the relation: ${f}_{n}={f}_{i}+{f}_{j}-{f}_{k}\left(i,j\ne k\right)$ .

We can define the term Y indicating the maximum values of the summation ${\eta}_{ijk}{d}_{ijk}^{2}$ for overall channels as follow:

$Y=MAX\left({{\displaystyle \sum}}^{\text{}}\text{\hspace{0.05em}}{\eta}_{ijk}\ast {d}_{ijk}^{2}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}i,j,k=1,2,\cdots ,N$ (8)

In a system with equal power per channel $\left({P}_{i}={P}_{j}={P}_{k}={P}_{s}\right)$ , the FWM power on the worst affected channel ( $n=w$ ) is given by:

${P}_{w}={k}^{2}\ast {P}_{s}{\left(\frac{M\ast {L}_{eff}}{{A}_{eff}}\right)}^{2}Y\ast {\text{e}}^{-a\cdot {L}_{a}}$ (9)

For the case of a system with different power levels, Equation (9) would only affect the Y term by adding each different power per channel. In such case Y would be:

$Y=MAX\left({{\displaystyle \sum}}^{\text{}}{P}_{ijk}\ast {\eta}_{ijk}\ast {d}_{ijk}^{2}\right)$ (10)

Nevertheless, for practical purposes all channels are considered to have equal power.

FWM can be considered as an interference signal, and supposing that the required optical signal to interference ratio is 20 dB (typical value) at the receiver to ensure no appreciable performance penalty [1] , the maximum launched power per channel ${P}_{m}$ can be obtained considering the worst affected channel and it can be approximated as follows:

$\frac{{P}_{m}\ast {\text{e}}^{-a\ast {L}_{a}}}{{P}_{w}}=100$ (11)

${P}_{m}=\frac{{A}_{eff}}{M\ast {L}_{eff}}\sqrt{\frac{1}{100\ast {k}^{2}\ast Y}}$ (12)

The effect of ASE noise due to inline amplifiers is estimated next. The total amount of ASE noise power at the receiver increases with the number of amplifiers. Therefore, in long distance amplified systems the ASE noise dominates over other noise sources (shot noise and thermal noise), and the Bit Error Rate (BER) performance is mainly determined by optical signal-to-noise ratio (OSNR). The SNR can be estimated from a simple Gaussian approximation in [13] :

$SNR=\frac{{Q}^{2}+Q\sqrt{\frac{{B}_{o}}{{B}_{e}}}}{\sqrt{\frac{{B}_{o}}{{B}_{e}}}}$ (13)

where ${B}_{o}$ is the bandwidth of the optical filter, ${B}_{e}$ is the electrical bandwidth of the receiver and Q a parameter related to BER given by:

$BER=\frac{1}{2}erfc\left(\frac{Q}{\sqrt{2}}\right)\approx \frac{{e}^{-\frac{{Q}^{2}}{2}}}{Q\sqrt{2\ast \pi}}$ (14)

The filter bandwidth should be large enough to pass the entire frequency contents of the selected channel but, at the same time, small enough to block the neighbouring channels and noise, typically
${B}_{o}~B$ and it represents the minimum optical bandwidth to avoid blocking signal [1] and
${B}_{e}<B$ , in our simulation we take
${B}_{o}~2B$ and
${B}_{e}~B/2$ to provide sufficient margin. To achieve a BER of 1 × 10^{−12}, Q = 7 and SNR = 15.75 (11.97 dB) from Equations (8) (9). Assuming an SNR = 20 dB (typical value), the total ASE noise at the detector can be expressed as [1] [14] :

${P}_{sp}=2\ast {n}_{sp}\ast \left(G-1\right)\ast h\ast f\ast {B}_{o}\ast M$ (15)

where h is Plancks constant (6.63 × 10^{−34} J・s), f is the centre frequency, G is the gain of the amplifier,
${n}_{sp}$ is the population inversion parameter and 2n_{sp} is the amplifier noise factor and its minimum value is 2 for an ideal amplifier [1] . Considering equal amplifier spacing (M = L/L_{a}) and that amplifier gain G compensates for fibre loss as:

$G=E\ast {\text{e}}^{a\ast {L}_{a}}$ (16)

where E is a factor accounting for possible losses introduced by fibre splices and other optical components inserted in the transmission line, we can assume the value of 1 for practical purposes. The minimum power per channel to ensure the required SNR (20 dB) can then be obtained as:

$P=200\ast {n}_{sp}\ast \left(G-1\right)\ast h\ast f\ast {B}_{o}\ast L/{L}_{a}$ (17)

The FWM effect imposes an upper limit and ASE imposes, on the contrary, a minimum limit to the power per channel value when the transmission distance is increased. Therefore, the intersection between P and P_{m} provides the maximum value of transmission distance (L = L_{max}) and P_{o}__{fwm} = P = P_{m} is the optimum transmission power as shown below:

${L}_{\text{max\_}fwm}={\left\{\frac{{A}_{eff}\ast {L}_{a}^{2}}{200\ast h\ast f\ast {n}_{sp}\ast {B}_{o}\ast \left(E\ast {\text{e}}^{a\ast {L}_{a}}-1\right)\ast {L}_{eff}}\ast \left[\sqrt{\frac{1}{100\ast {k}^{2}\ast Y}}\right]\right\}}^{1/2}$ (18)

${P}_{o\text{\_}fwm}={\left\{200\ast h\ast f\ast {n}_{sp}\ast {B}_{o}\ast \left(E{\text{e}}^{a{L}_{a}}-1\right)\left(1/{L}_{eff}\right)\left[\sqrt{\frac{1}{100\ast {k}^{2}\ast Y}}\right]\right\}}^{1/2}$ (19)

Equations (13) and (14) are valid only for a systems limited by FWM and ASE noise.

Next it is important to calculate the dispersion limit. The limit of transmission distance for Intensity Modulation with Direct Detection (IM/DD) systems due to dispersion effect is given by the following expression [14] :

${L}_{D}=\frac{c}{2\ast {B}^{2}\ast {\lambda}^{2}\ast \left|D\right|}$ (20)

The effect of SRS-ASE can be considered in a similar way. In this case, an upper and lower limit can be obtained for the power per channel due to SRS-ASE. To do this, a worst-case neglecting walk-off effect and dispersion are considered. The maximum power to ensure a SNR degradation of less than 0.5 dB in the worst channel is given by [15] :

${P}_{\text{max\_}srs}=\frac{8.7\ast {10}^{12}}{N\left(N-1\right)M\ast \Delta f\ast {L}_{eff}}$ (21)

The intersection between P_{max_srs} and Equation (12) gives a maximum value of transmission distance considering SRS and ASE noise:

${L}_{\text{max}\_srs}=\sqrt{\frac{8.7\ast {10}^{12}\ast {L}_{A}^{2}}{200\ast h\ast f\ast {n}_{sp}\ast {B}_{o}\ast N\left(N-1\right)\left(E\ast {\text{e}}^{\alpha \ast {L}_{a}}-1\right)\Delta f\ast {L}_{eff}}}$ (22)

So, the optimum transmission power for a system limited only by SRS and ASE noise is given by:

${P}_{o\_srs}=41.7\ast {10}^{6}\ast \sqrt[]{\frac{h\ast f\ast {n}_{sp}\ast {B}_{o}\left(E\ast {\text{e}}^{a\ast {L}_{a}}-1\right)}{N\ast W\ast {L}_{eff}}}$ (23)

The intersection between P_{max_srs} and Equation (21) gives a maximum value of transmission distance considering SRS and ASE noise:

Equations (13) and (17) can be used to evaluate the optimal length for multi-span DWDM systems limited by FWM and SRS separately. Then, we can analyse overlapping regions where both nonlinearities fix the same optimal length. This intersection can provide us with a more general limit for the maximum length in the DWDM system considering the combined effects of FWM, SRS, and ASE respectively. With this procedure one can obtain the spectral limits of a DWDM system without considering the simultaneous analysis of both nonlinearities [16] . In this case, we have to be clear take that we are not interested to study the gradual degradation process through the fibre due to FWM an SRS, but rather the final degradation due to both nonlinearities in the worst channel.

3. Simulation Results and Discussion

The amplifier spacing (L_{a}) can be calculated according to typical values of receiver sensitivity as: S = −24 dBm, α = 0.25 dB/km (a in dB), P_{T} = 1 mW.

A value of L_{a} = 75 km (
${L}_{a}<\left({P}_{T}-S\right)/\alpha $ ) and a bit rate B = 2.5 Gbps per channel is considered in all our simulations.

Figure 2 shows the FWM power per channel (Equation (4)) evaluated for different values of D, with 31 channels and 1 nm channel spacing. As we can see, the FWM power depends on the location of the zero dispersion wavelength. At D = 0 ps/((km・nm)) the maximum FWM power in the central channel is obtained due to the maximum efficiency ${\eta}_{ijk}$ is obtained too, it agrees with [12] . At $D\ne 0$ the FWM power decreases. If the zero dispersion wavelength is far of the WDM signal bandwidth, the right or left extreme channels suffer most (depending on the sign of D).

Figure 3 and Figure 4 show the maximum length considering the combined effects of SRS-ASE, FWM-ASE and Dispersion in a separated form for a 30- and 70-channel system, respectively. Three intervals or domains can be identified by the intersection of the curves, SRS-ASE, FWM-ASE and Dispersion domain.

The evolution of these intersection can be approximated through the following expressions:

FWM-ASE interval domination.

$D<8.9\ast {10}^{-5}\ast {N}^{2}+0.011\ast N+1.1$ (24)

Figure 5 shows the evolution of the optimal power for different intervals in agreement with Equation (19). As it can be seen, they are higher than the power threshold for SRS and SBS (Stimulated Brillouin Scattering) [17] [18] [19] [20] [21] .

Table 1 is obtained via Equations (13) (14) and (19). It shows the maximum value of optimum power (P_{o_fwm}) and maximum length (L_{max}) in a FWM-ASE domain. P_{o_fwm} and L_{max} decrease depending of the interval of D, and value of D.

For a DWDM system with N = 100 channels (∆f = 37.8 GHz), the maximum value of P_{o_fwm} is 0.15 mW with a value of D = 3.09 ps/km・nm, and it can decrease depending on D.

P_{o_fwm} varies in a range of 0.78 - 0.15 mW for a range N of 1.21 - 3.09 ps/km・nm up to N = 100 channels.

SRS-ASE interval domination.

$8.9\ast {10}^{-5}\ast {N}^{2}+0.011\ast N+1.1<D<-0.00066\ast {N}^{2}+0.18\ast N+3.8$ (25)

Figure 6 shows a curve of the evolution of the optimal power for different number of channels in agreement with (20).

Table 2 is obtained using Equations (17) (18) and (20). It shows the maximum value of optimum power (P_{o_srs}) and maximum length (L_{max}) in a SRS-ASE domain. P_{o_srs} and L_{max} decrease depending on the value of the design parameters (N, D, ∆f).

For a DWDM system of N = 100 channels (∆f = 37.8 GHz), the maximum value of P_{o_srs} is 0.005 mW with a value between D = 3.09 to 15.2 ps/km・nm, and it is constant in all interval of D. The interval of D is higher than FWM-ASE case, but the level of optimum power is lower.

Finally, the interval dominated by the effect of linear dispersion D is given by:

$D>-0.00066\ast {N}^{2}+0.18\ast N+3.8$ (26)

Table 1. P_{o_fwm} (mW) and L_{max} (km) in FWM-ASE domain, for different design parameters N (channels), ∆f (GHz), D (ps/km・nm).

Table 2. P_{o_srs} (mW) and L_{max} (km) in SRS-ASE domain, for different design parameters N (channels), ∆f (GHz), D (ps/km・nm).

Figure 2. FWM power per channel for positive and negative values of D.

Figure 3. Limitation due to FWM-SRS-ASE-Dispersion for N = 20 channels.

Figure 4. Limitation due to FWM-SRS-ASE-Dispersion for N = 70 channels.

Figure 5. Optimum power as a function of D, for a system limited by FWM y ASE, for 30, 60 and 100 channels.

Figure 6. Optimum power as a function of N, for a system limited by SRS and ASE.

4. Conclusion

We have explored the limitations of a multi-span DWDM system imposed by FWM, ASE, dispersion and SRS for up to 100 channels. We have determined the conditions under which DWDM systems will be limited mainly by FWM and ASE, or SRS and ASE separately and altogether. In particular, for the FWM-ASE nonlinearities, the optimum power transmission per channel was calculated. The effect of SBS can be neglected because the power needed for these nonlinear effects to appear above their threshold is not reached and Brillouin gain bandwidth is too short. The optimum link length that can be achieved for these systems was calculated from the channel number or inter-channel spacing and the total dispersion based on the crossover of the predicted FWM, SRS and dispersion limit. In particular, for channel spacing higher than 125 GHz, the SRS-ASE nonlinearity becomes dominant in comparison to the FWM-ASE effect, but at the same time, according to the proposed methodology, one can adjust the design parameters in order to obtain a same effect of the FWM and SRS nonlinearities including the ASE noise respectively. Our results are of great interest for the design of novel optimized multi-span DWDM systems.

Acknowledgements

The authors are grateful towards UNACAR, IPN-CIITEC and CONACYT from Mexico, and Optical Sciences Group at Twente University in the Netherlands.

References

[1] Agrawal, G.P. (2010) Fiber-Optic Communication Systems. John Wiley & Sons.

https://doi.org/10.1002/9780470918524

[2] Rostami, A., Rahbari, J. and Andalib, A. (2013) Investigation of Power Penalty in WDM Systems for Dispersion Managed Fibers. Optik, 124, 2072-2075.

https://doi.org/10.1016/j.ijleo.2012.06.065

[3] Karfaa, Y.M., et al. (2007) Effects of Four-Wave Mixing Crosstalk in WDM Networks on the Transmitted Optical Frequencies and Wavelengths of Channels for Various Fiber Types. Asia-Pacific Conference on Applied Electromagnetics APACE, Melaka, 1-5.

https://doi.org/10.1109/APACE.2007.4603976

[4] Djordjevic, I.B. (2001) Transmission Limitation of WDM Transmission Systems with Dispersion Compensated Links in the Presence of Fiber Nonlinearities. Conference on Telecommunications in Modern Satellite, Cable and Broadcasting Service TELSIKS, 5th International, Nis, 19-21.

[5] Rana, Md.M. and Khan, M.A.G. (2012) Dense Wavelength Division Multiplexing Optical Network System. International Journal of Electrical and Computer Engineering, 2, 203-206

[6] Kowalewski, M., et al. (2000) Nonlinear Interactions in Wavelength-Multiplexed Optical Fibre Telecommunications Systems. Journal of Optics A, 2, 319-326.

https://doi.org/10.1088/1464-4258/2/4/314

[7] Gurmeet Kaur, M.I. Singh, M. and Patterh, S. (2009) Theoretical Investigations of the Combined Effect of Fibre Nonlinearities, Amplifier and Receiver Noise in a DWDM Transmission System. Journal of Optical and Fiber Communications Research, 6, 1-10.

https://doi.org/10.1007/s10297-009-9002-8

[8] Wu, M. and Way, W.I. (2004) Fiber Nonlinearity Limitations in Ultra-Dense WDM Systems. IEEE Journal of Lightwave Technology, 22, 1483-1497.

https://doi.org/10.1109/JLT.2004.829222

[9] Louchet, H., et al. (2005) Analytical Model for the Design of Multispan DWDM Transmission Systems. IEEE Photonics Technology Letters, 17, 247-249.

https://doi.org/10.1109/LPT.2004.838142

[10] Maeda, M.W., et al. (1990) The Effect of Four-Wave Mixing in Fibers on Optical Frequency-Division Multiplexed Systems. IEEE Journal of Lightwave Technology, 8, 1402-1408.

https://doi.org/10.1109/50.59171

[11] Ellis, D., Mccarthy, M.E., Alkhateeb, M.A.Z., Sorokina, M. and Doran, N.J. (2017) Performance Limits in Optical Communications Due to Fiber Nonlinearity. Advances in Optics and Photonics, 9, 429-503.

https://doi.org/10.1364/AOP.9.000429

[12] Inoue, K. (1992) Four Wave Mixing in an Optical Fiber in the Zero-Dispersion Wavelength Region. IEEE Journal of Lightwave Technology, 10, 1533-1561.

https://doi.org/10.1109/50.184893

[13] Marcuse, D. (1990) Derivation of Analytical Expressions for the Bit-Error Probability in Lightwave Systems with Optical Amplifiers. IEEE Journal of Lightwave Technology, 8, 1816-1823.

[14] Henry, P.S. (1985) Lightwave Primer. IEEE Journal of Quantum Electronics, 21, 1862-1879.

https://doi.org/10.1109/JQE.1985.1072601

[15] Chraplyvy, A.R. and Tkach, R.W. (1993) What Is the Actual Capacity of Single-Mode Fibers in Amplified Lightwave Systems. IEEE Photonics Technology Letters, 5, 666-668.

https://doi.org/10.1109/68.219704

[16] Vanholsbeeck, F., Coen, S., et al. (2005) Coupled-Mode Analysis of Stimulated Raman Scattering and Four Wave Mixing in Wavelength-Division Multiplexed Systems. Optics Communications, 250, 191-201.

https://doi.org/10.1016/j.optcom.2005.02.011

[17] Sanchez-Lara, R., Alvarez-Chavez, J.A., et al. (2015) Threshold and Maximum Power Evolution of Stimulated Brillouin Scattering and Rayleigh Backscattering in a SM Fiber-Segment. Laser Physics, 25, Article ID: 035103.

https://doi.org/10.1088/1054-660X/25/3/035103

[18] Hekmat, M.J., et al. (2013) Study of the Stimulated Brillouin Scattering Power Threshold in High Power Double-Clad Fiber Lasers. Laser Physics, 23, Article ID: 025104.

https://doi.org/10.1088/1054-660X/23/2/025104

[19] Kovalev, V.I. (2014) Peculiarities of the Single and Repetitive Pulse Threshold Pump Power for Stimulated Brillouin Scattering in Long Optical Fibres. Laser Physics, 24, Article ID: 085407.

https://doi.org/10.1088/1054-660X/24/8/085407

[20] Tithi, F.H., et al. (2012) Compensation of SBS Using Cross-Phase Modulation in WDM Transmission System. 7th International Conference on Electrical & Computer Engineering, Dhaka, 39-42.

https://doi.org/10.1109/ICECE.2012.6471479

[21] Sanchez-Lara, R., Alvarez-Chavez, J.A., et al. (2014) Cascaded Raman Generation of the N’th Stokes in Single Mode-Fibers. Optical Engineering, 53, Article ID: 075102.