Soliton collision and Raman gain regimes in continuous-wave pumped supercontinuum generation

doi:10.1364/OE.14.009391

« Show journal navigation

Soliton collision and Raman gain regimes in continuous-wave pumped supercontinuum generation

Michael H. Frosz, Ole Bang, and Anders Bjarklev »View Author Affiliations

Optics Express, Vol. 14, Issue 20, pp. 9391-9407 (2006)
http://dx.doi.org/10.1364/OE.14.009391

View Full Text Article

Acrobat PDF (482 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Full Text
Article Info
References (43)
Cited By
Figures (10)
Multimedia (1)
Metrics

Abstract

We numerically investigate supercontinuum generation using continuous-wave pumping. It is found that energy transfer during collision of solitons plays an important role. The relative influence of Raman gain on spectral broadening is shown to depend on the width of the calculation time window. Our results indicate that increasing the spectral linewidth of the pump can decrease the supercontinuum spectral width. Using a fiber with smaller dispersion at the pump wavelength reduces the required fiber length by decreasing the temporal width of the solitons formed from modulation instability. This also reduces the sensitivity to the pump spectral linewidth.

1. Introduction

Spatially coherent light with an extremely broad optical spectrum (>100 nm) can be produced by pumping light with a narrow spectrum and high intensity through a nonlinear medium [1

R. R. Alfano, ed., The Supercontinuum Laser Source: Fundamentals With Updated References , 2nd ed. (Springer-Verlag, New York, 2005). ISBN 0-387-24504-9.

]. The interaction between the high intensity light and the nonlinear medium initiates several mechanisms, which broaden the spectrum of the light into a so-called supercontinuum (SC) as it propagates through the medium. Photonic crystal fibers (PCFs) are often used as the nonlinear medium, not only because the silica core confining the light can be made very small, thus increasing the intensity of the light, but also because the dispersion properties of these fibers are highly controllable. The zero dispersion wavelength (ZDW) can, e.g., be shifted to the visible region of the optical spectrum [2

J. K. Ranka, R. S. Windeler, and A. J. Stentz, “Visible continuum generation in air-silica microstructure optical fibers with anomalous dispersion at 800 nm,” Opt. Lett. 25, 25–27 (2000). [CrossRef]

], or there can even be two ZDWs in the optical spectrum [3

K. M. Hilligsøe, T. V. Andersen, H. N. Paulsen, C. K. Nielsen, K. Mølmer, S. Keiding, R. Kristiansen, K. P. Hansen, and J. J. Larsen, “Supercontinuum generation in a photonic crystal fiber with two zero dispersion wavelengths,” Opt. Express 12, 1045–1054 (2004). http://www.opticsexpress.org/abstract.cfm?URI=OPEX-12-6-1045. [CrossRef] [PubMed]

–5

M. H. Frosz, P. Falk, and O. Bang, “The role of the second zero-dispersion wavelength in generation of supercontinua and bright-bright soliton-pairs across the zero-dispersion wavelength,” Opt. Express 13, 6181–6192 (2005). http://www.opticsexpress.org/abstract.cfm?URI=OPEX-13-16-6181. [CrossRef] [PubMed]

]. To some degree, similar dispersion engineering can be performed by tapering standard telecommunications fiber [6

T. A. Birks, W. J. Wadsworth, and P. St. J. Russell, “Supercontinuum generation in tapered fibers,” Opt. Lett. 25, 1415–1417 (2000). [CrossRef]

]. By choosing a PCF with a suitable dispersion-profile, the broadening mechanisms can be made more efficient for generating light in particular wavelength regions [2

Pump light with high intensity has traditionally been achieved using femtosecond lasers, offering pulses with several kilowatts of peak power [2

, 4

G. Genty, M. Lehtonen, H. Ludvigsen, and M. Kaivola, “Enhanced bandwidth of supercontinuum generated in microstructured fibers,” Opt. Express 12, 3471–3480 (2004). http://www.opticsexpress.org/abstract.cfm?URI=OPEX-12-15-3471. [CrossRef] [PubMed]

, 7

J. M. Dudley, L. Provino, N. Grossard, H. Maillotte, R. S. Windeler, B. J. Eggleton, and S. Coen, “Supercontinuum generation in air-silica microstructured fibers with nanosecond and femtosecond pulse pumping,” J. Opt. Soc. Am. B 19, 765–771 (2002). [CrossRef]

]. Femtosecond lasers are typically complex and bulky; therefore more compact picosecond [8

S. Coen, A. H. L. Chau, R. Leonhardt, J. D. Harvey, J. C. Knight, W. J. Wadsworth, and P. St. J. Russell, “Supercontinuum generation by stimulated Raman scattering and parametric four-wave mixing in photonic crystal fibers,” J. Opt. Soc. Am. B 19, 753–764 (2002). [CrossRef]

] and even nanosecond [7

] pump sources have also been used. Nanosecond pumped SC generation technology has reached a point of maturity where commercial sources are available (see, e.g., www.koheras.com, SuperK source).

Advances in fiber laser technology in recent years have also led to continuous-wave (CW) fiber lasers with sufficiently high power (~0.1–15 W) to be used for SC generation [9

M. Prabhu, N. S. Kim, and K. Ueda, “Ultra-Broadband CW Supercontinuum Generation Centered at 1483.4 nm from Brillouin/Raman Fiber Laser,” Jpn. J. Appl. Phys. 39, L291–L293 (2000). http://dx.doi.org/10.1143/JJAP.39.L291. [CrossRef]

–11

A. V. Avdokhin, S. V. Popov, and J. R. Taylor, “Continuous-wave, high-power, Raman continuum generation in holey fibers,” Opt. Lett. 28, 1353–1355 (2003). [CrossRef] [PubMed]

]. It has been demonstrated that a CW-pumped SC source provides sufficiently high output power and broad bandwidth to be used for ultrahigh resolution optical coherence tomography (OCT) [12

P.-L. Hsiung, Y. Chen, T. Ko, J. G. Fujimoto, C. J. S. de Matos, S. V. Popov, J. R. Taylor, and V. P. Gapontsev, “Optical coherence tomography using a continuous-wave, high-power, Raman continuum light source,” Opt. Express 12, 5287–5295 (2004). http://www.opticsinfobase.org/abstract.cfm?URI=oe-12-22-5287. [CrossRef] [PubMed]

]. Experimental characterization of the noise of the CW pumped SC showed that the relative intensity noise (RIN) can be 30–50 dB lower than that achieved using a femtosecond pump [13

C. J. S. de Matos, S. V. Popov, and J. R. Taylor, “Temporal and noise characteristics of continuous-wave-pumped continuum generation in holey fibers around 1300 nm,” Appl. Phys. Lett. 85, 2706–2708 (2004). [CrossRef]

]. Other experiments with a CW pumped SC found RIN values comparable to those with a femtosecond pump [14

J. H. Lee, Y.-G. Han, and S. B. Lee, “Experimental study on seed light source coherence dependence of continuous-wave supercontinuum performance,” Opt. Express 14, 3443–3452 (2006). http://www.opticsinfobase.org/abstract.cfm?URI=oe-14-8-3443. [CrossRef] [PubMed]

CW pumping for SC generation has been investigated both by pumping in the anomalous-dispersion region (ADR) [10

J. W. Nicholson, A. K. Abeeluck, C. Headley, M. F. Yan, and C. G. Jørgensen, “Pulsed and continuous-wave supercontinuum generation in highly nonlinear, dispersion-shifted fibers,” Appl. Phys. B 77, 211–218 (2003). http://dx.doi.org/10.1007/s00340-003-1201-z. [CrossRef]

,11

A. V. Avdokhin, S. V. Popov, and J. R. Taylor, “Continuous-wave, high-power, Raman continuum generation in holey fibers,” Opt. Lett. 28, 1353–1355 (2003). [CrossRef] [PubMed]

,13

,15

A. Mussot, E. Lantz, H. Maillotte, T. Sylvestre, C. Finot, and S. Pitois, “Spectral broadening of a partially coherent CW laser beam in single-mode optical fibers,” Opt. Express 12, 2838–2843 (2004). http://www.opticsexpress.org/abstract.cfm?URI=OPEX-12-13-2838. [CrossRef] [PubMed]

–17

F. Vanholsbeeck, S. Martin-Lopez, M. González-Herráez, and S. Coen, “The role of pump incoherence in continuous-wave supercontinuum generation,” Opt. Express 13, 6615–6623 (2005). http://www.opticsexpress.org/abstract.cfm?URI=OPEX-13-17-6615. [CrossRef] [PubMed]

] and the normal-dispersion region (NDR) [18

A. K. Abeeluck and C. Headley, “Continuous-wave pumping in the anomalous- and normal-dispersion regimes of nonlinear fibers for supercontinuum generation,” Opt. Lett. 30, 61–63 (2005). [CrossRef] [PubMed]

,19

B. Barviau, S. Randoux, and P. Suret, “Spectral broadening of a multimode continous-wave optical field propagating in the normal dispersion regime of a fiber,” Opt. Lett. 31, 1696–1698 (2006). [CrossRef] [PubMed]

]. In the present work we focus on pumping in the ADR. The spectral broadening achieved usingCW pumping in the ADR has previously been explained as being seeded by modulation instability (MI), which causes breakup of the CW into solitons. The solitons are then red-shifted due to the soliton self-frequency shift (SSFS) [13

, 15

–17

]. During the red-shift the solitons may also transfer energy to blue-shifted dispersive waves. Thus, after MI-induced breakup into solitons, the spectral broadening is essentially caused by the same mechanisms as when using fs pumping [4

Solitons generated directly from MI can be too temporally broad to undergo any significant red-shift [20

E. A. Golovchenko, P. V. Mamyshev, A. N. Pilipetskii, and E. M. Dianov, “Numerical analysis of the Raman spectrum evolution and soliton pulse generation in single-mode fibers,” J. Opt. Soc. Am. B 8, 1626–1632 (1991). http://www.opticsinfobase.org/abstract.cfm?URI=josab-8-8-1626. [CrossRef]

]. It was shown by Islam et al. that a significant red-shift can still occur because MI-generated solitons undergo collisions, in which they transfer energy between each other [21

M. N. Islam, G. Sucha, I. Bar-Joseph, M. Wegener, J. P. Gordon, and D. S. Chemla, “Femtosecond distributed soliton spectrum in fibers,” J. Opt. Soc. Am. B 6, 1149–1158 (1989). http://www.opticsinfobase.org/abstract.cfm?URI=josab-6-6-1149. [CrossRef]

]. It is then possible for one of the solitons to acquire sufficient energy to undergo a large red-shift. In this work, we examine the conditions for proper modeling of the soliton collisions and their effect on the spectral broadening. The influence of the CW pump spectral linewidth and the fiber dispersion on the spectral broadening is also investigated.

The paper is organized as follows. In Section 2 we provide an overviewof previous numerical modeling of CW pumped SC generation, and justify our choice of how to model the CW input. Section 3 then describes how we simulate the propagation of the partially coherent CW beam along the fiber, and the implemented phase noise model. In Section 4 we describe the primary mechanisms involved in the spectral broadening during propagation along the fiber. It is shown that soliton collisions play an important role in CW-pumped supercontinuum generation for the parameters used here. It is also shown how the influence of soliton collisions is practically absent, if the time window of the calculations is too narrow. We then present results showing that the output spectrum can differ significantly depending on the spectral linewidth of the CW pump laser. Finally, we show that SC generation in a fiber with smaller dispersion at the pump wavelength is less sensitive to the spectral linewidth of the pump laser.

2. Previous numerical modeling

One important aspect of modeling the nonlinear propagation of a CW beam in an optical fiber is how to model the CW input.

Kobtsev and Smirnov used the one photon per mode approach [16

S. M. Kobtsev and S. V. Smirnov, “Modelling of high-power supercontinuum generation in highly nonlinear, dispersion shifted fibers at CW pump,” Opt. Express 13, 6912–6918 (2005). http://www.opticsexpress.org/abstract.cfm?URI=OPEX-13-18-6912. [CrossRef] [PubMed]

]. This method has previously been used to model SC generation using picosecond pulses [8

] and consists of phenomenologically adding one photon with random phase to each frequency bin of the input field. It was shown by Smith in 1972 that this fictitious injection of photons to the input leads to the same output power at the Raman Stokes wavelength (downshifted from the pump by 13.2 THz) as Raman amplification of spontaneous emission along the fiber length [22

R. G. Smith, “Optical power handling capacity of low loss optical fibers as determined by stimulated Raman and Brillouin scattering,” Appl. Opt. 11, 2489–2494 (1972). [CrossRef] [PubMed]

]. However, as we show in Fig. 1(right) this approach suffers mainly from two drawbacks: (1) There is no inherent spectral linewidth, and (2) the power spectral density at the Raman Stokes wavelength can be significantly underestimated. As mentioned in Subsection 3.2 it is reasonable to assume a Lorentzian spectrum for a partially coherent CW field. We show in Fig. 1(right) that a Lorentzian spectrum for a CW field with an average power of 10W, center wavelength of 1064 nm and a full width at half-maximum of 30 GHz (~0.1 nm) has significantly more spectral power at the Raman Stokes wavelength 1116 nm, than a CW field containing one photon per mode. The one photon per mode method is therefore unsuited to investigate the influence of the pump laser spectral linewidth and may underestimate the Raman Stokes power spectral density at the input. This is critical when Raman scattering is dominant because the input Raman Stokes power acts as a seed for the exponential gain along the fiber.

Vanholsbeeck et al. also take a phenomenological approach by taking the input spectrum of the field to be

\tilde{E} (ν) = \sqrt{S_{m} (ν)} \times ϕ (ν)

, where ϕ(ν) is a random spectral phase and S _m(ν) is the measured power spectral density of the pump laser to be modelled [17

]. This approach mathematically assures that the power spectrum of the input field is as desired, |Ẽ(ν)|²=S _m(ν). But the choice of randomly imposing a spectral phase has no physical justification. The strong temporal fluctuations of the input field E(t) obtained by inverse Fourier transform of Ẽ(ν) are therefore not necessarily present in a real partially coherent CW beam, as also pointed out by the authors in Ref. [17

]. Recently, it was pointed out by Barviau et al. that modifying the statistics of the spectral phase ϕ(ν) has a significant influence on the spectral evolution along the fiber [19

]. One must therefore be careful when using arbitrary statistical properties for the spectral phase.

Cavalcanti et al. used a time-independent average field on which noise appears as small fluctuations in amplitude and phase to study the influence of MI on the statistical properties of a partially coherent beam [23

S. B. Cavalcanti, G. P. Agrawal, and M. Yu, “Noise amplification in dispersive nonlinear media,” Phys. Rev. A 51, 4086–4092 (1995). http://dx.doi.org/10.1103/PhysRevA.51.4086. [CrossRef] [PubMed]

]. Mussot et al. disregarded amplitude fluctuations and included only phase noise to investigate the spectral broadening mechanism [15

]. As shown in Subsection 3.2, the phase noise can be directly related to the measurable linewidth of a CW laser. Furthermore, the underlying phase-diffusion model is widely accepted as an appropriate model for lasers operating far above threshold (see e.g. Refs. [15

, 23

], and Refs. therein). Since this method involves both a physically justified choice of phase noise statistics and includes a finite spectral linewidth, it is physically reasonable and used in the work presented here. We can therefore investigate the influence of the spectral linewidth of the CW pump on the SC generation.

Fig. 1. Left: Calculated dispersion profiles for fibers considered here. Right: Comparison of power spectral density S(λ) for two noise models: CW with one photon per mode (blue, solid) and CW with phase noise corresponding to FWHM linewidth of 30 GHz (~0.1 nm), Lorentzian shaped (green, dashed). Average power of input field is 10W. Red dotted curve: resulting spectrum when the CW is replaced by a super-Gaussian pulse of 30 ps width and 10 W peak power (the spectrum is scaled so integrating over the spectrum gives 10 W to allow comparison with the CW spectrum).

3. Theory

3.1. Propagation equation

We model the propagation of the CW field using the generalized nonlinear Schrödinger equation (NLSE) [24

G. P. Agrawal, Nonlinear Fiber Optics , 3rd ed. (Academic Press, San Diego, CA, USA, 2001).

, 25

P. L. François, “Nonlinear propagation of ultrashort pulses in optical fibers: total field formulation in the frequency domain,” J. Opt. Soc. Am. B 8, 276–293 (1991). http://www.opticsinfobase.org/abstract.cfm?URI=josab-8-2-276. [CrossRef]

]

\frac{\partial E}{\partial z} = i \sum_{m \geq 2} \frac{i^{m} β_{m}}{m!} \frac{\partial^{m} E}{\partial T^{m}} - \frac{α}{2} E + i γ (ω) [1 + \frac{i}{ω_{0}} \frac{\partial}{\partial T}] [E (z, T) \int_{- \infty}^{\infty} R (T') {∣ E (z, T - T') ∣}^{2} d T'],

(1)

where E(z,T) is the field envelope in a retarded time frame T=t-β1z moving with the group velocity 1/β1 of the pump, along the fiber axis z. ω₀ is the carrier angular frequency of the CW input. γ(ω)=n ₂ω₀/[cA_eff (ω)] is the nonlinear parameter, where n ₂=2.6×10^-20 m²/W is the nonlinear-index coefficient for silica, c is the speed of light in vacuum, and A _eff is the effective core area [24

G. P. Agrawal, Nonlinear Fiber Optics , 3rd ed. (Academic Press, San Diego, CA, USA, 2001).

]. The dispersion parameters β _m are obtained from a polynomial fit to the dispersion profile calculated using fully-vectorial plane-wave expansions [26

S. G. Johnson and J. D. Joannopoulos, “Block-Iterative Frequency-Domain Methods for Maxwell’s Equations in a Planewave Basis,” Opt. Express 8, 173–190 (2001). http://www.opticsexpress.org/abstract.cfm?URI=OPEX-8-3-173. [CrossRef] [PubMed]

Dispersion profiles for the PCFs considered in this paper are shown in Fig. 1(left). The PCF with pitch Λ=1.72 µm and relative hole size d/Λ=0.65 corresponds to a fiber made by Crystal Fiber A/S. It was used in Refs. [11

A. V. Avdokhin, S. V. Popov, and J. R. Taylor, “Continuous-wave, high-power, Raman continuum generation in holey fibers,” Opt. Lett. 28, 1353–1355 (2003). [CrossRef] [PubMed]

, 12

] [note that both papers state that d=0.65 µm, but this does not correspond with the reported core size of d _core=2.3 µm (d _core≈2Λ-d) and is believed to be a misprint [27

K. P. Hansen, Crystal Fibre A/S, Personal communication.

]]. To investigate the influence of the dispersion profile, we have also considered a PCF with a relative hole size d/Λ=0.378. The transverse field distribution E(x,y,ω) of the fundamental mode was also calculated with fully-vectorial plane-wave expansions. From this A _eff(ω) was calculated using the more general definition suitable for fibers, where some of the field energy may reside in the air-holes [28

J. Lægsgaard, N. A. Mortensen, and A. Bjarklev, “Mode areas and field-energy distribution in honeycomb photonic bandgap fibers,” J. Opt. Soc. Am. B 20, 2037–2045 (2003). [CrossRef]

]. R(t) is the Raman response function [24

G. P. Agrawal, Nonlinear Fiber Optics , 3rd ed. (Academic Press, San Diego, CA, USA, 2001).

, 29

K. J. Blow and D. Wood, “Theoretical description of transient stimulated Raman scattering in optical fibers,” IEEE J. Quantum Electron. 25, 2665–2673 (1989). [CrossRef]

]

R (t) = (1 - f_{R}) δ (t) + f_{R} h (t) = (1 - f_{R}) δ (t) + f_{R} \frac{τ_{1}^{2} + τ_{2}^{2}}{τ_{1} τ_{2}^{2}} \exp (- t ⁄ τ_{2}) \sin (t ⁄ τ_{1}) Θ (t)

(2)

where f _R=0.18 is the fractional contribution of the delayed Raman response, τ₁=12.2 fs, and τ₂=32 fs. Θ(t) is the Heaviside step function.

The propagation Eq. (1) is solved using the split-step Fourier method [24

G. P. Agrawal, Nonlinear Fiber Optics , 3rd ed. (Academic Press, San Diego, CA, USA, 2001).

]. The relative change in photon number (a measure of the numerical error that is ideally zero in the absence of loss [29

K. J. Blow and D. Wood, “Theoretical description of transient stimulated Raman scattering in optical fibers,” IEEE J. Quantum Electron. 25, 2665–2673 (1989). [CrossRef]

]) was less than 5.5% for all results presented here.

3.2. Noise model

We model the input field envelope as [15

, 23

]

E (0, T) = \sqrt{P_{0}} \exp [i δ ϕ (T)]

(3)

where δϕ is a small fluctuation with zero ensemble average, 〈δϕ〉=0. P ₀ is the power of the CW input field. We will model δϕ as a Gaussian random process.

In the following, we determine how to properly relate δϕ to a measurable quantity, namely the full width at half-maximum (FWHM) power spectral linewidth Δν _FWHM of the CW input. The random phase fluctuations can be seen as arising from a random fluctuation ν_R of the CW frequency ν₀, so that the instantaneous frequency ν_i is [30

J. W. Goodman, Statistical Optics (John Wiley & Sons Inc., 2000). ISBN 0471399167.

]

ν_{i} = ν_{0} + \frac{1}{2 π} \frac{d (δ ϕ)}{d t} = ν_{0} + ν_{R} (t) .

(4)

It can be assumed in most cases that ν _R (t) is a zero mean, statistically stationary fluctuation [30

J. W. Goodman, Statistical Optics (John Wiley & Sons Inc., 2000). ISBN 0471399167.

]. Here, as in Ref. [23

], we model ν _R (t) as Gaussian white noise with zero mean and variance

σ_{ν R}^{2}

. The phase fluctuation δϕ(t) is then found from [30

J. W. Goodman, Statistical Optics (John Wiley & Sons Inc., 2000). ISBN 0471399167.

]

δ ϕ (t) = 2 π \int_{- \infty}^{t} ν_{R} (ξ) d ξ,

(5)

and is a statistically nonstationary random process. In the following, we determine the relation between

σ_{ν R}^{2}

and Δν_FWHM.

It can be shown that the structure function Dδϕ(τ)=〈[δϕ(t+τ)-δϕ(t)]²〉 is given by [30

J. W. Goodman, Statistical Optics (John Wiley & Sons Inc., 2000). ISBN 0471399167.

]

D_{δ ϕ} (τ) = 8 π^{2} τ \int_{0}^{\infty} Γ_{ν_{R}} (η) d η,

(6)

where Γ_νR is the autocorrelation function of ν _R (t). Using the Wiener-Khinchin theorem [30

J. W. Goodman, Statistical Optics (John Wiley & Sons Inc., 2000). ISBN 0471399167.

]

Γ_{ν_{R}} (τ) = \int_{- \infty}^{\infty} S_{ν_{R}} (ν) \exp (- i 2 π ν τ) d ν,

(7)

where S_νR (ν) is the power spectral density of ν_R (t), we obtain for white noise (i.e. Sν_R frequency independent) over a limited bandwidth B:

Γ_{ν_{R}} (τ) = S_{ν_{R}} \int_{- B ⁄ 2}^{B ⁄ 2} \exp (- i 2 π ν τ) d ν .

(8)

The bandwidth limit is introduced because the split-step Fourier method implies both a finite time- and frequency-domain in the calculations. Since

σ_{ν R}^{2}

=Γν_R (0) [30

J. W. Goodman, Statistical Optics (John Wiley & Sons Inc., 2000). ISBN 0471399167.

], we have

σ_{ν_{R}}^{2} = Γ_{ν_{R}} (0) = S_{ν_{R}} \int_{- B ⁄ 2}^{B ⁄ 2} d ν = S_{ν_{R}} B \Rightarrow S_{ν_{R}} = \frac{σ_{ν_{R}}^{2}}{B}

(9)

Inserting this in Eqs. (6) and (8), we obtain

D_{δ ϕ} (τ) = 8 π^{2} τ \frac{σ_{ν_{R}}^{2}}{B} \int_{0}^{\infty} \int_{- B ⁄ 2}^{B ⁄ 2} \exp (- i 2 π ν η) d ν d η = 4 π^{2} τ \frac{σ_{ν_{R}}^{2}}{B} .

(10)

The structure function Dδφ (τ) is related to Δν_FWHM by [23

] D_δϕ(τ)=2πΔν_FWHMτ, which together with Eq. (10) finally gives the following relation between σ ² ν_R and Δν_FWHM:

σ_{ν_{R}}^{2} = \frac{Δ ν_{FWHM} B}{2 π}

(11)

In summary, we model ν_R(t) as Gaussian white noise, with zero mean and variance

σ_{ν R}^{2}

. Equation (5) then gives the phase fluctuation δϕ(t), which is inserted in Eq. (3) to obtain the input field for the propagation simulation. Due to the assumptions of the statistical properties of the phase noise, the spectrum of the input field has a Lorentzian line shape [23

], with a FWHM given by Δν_FWHM.

3.3. Numerical considerations

An ideal CW field has an infinite extent in time, but with the split-step Fourier method we can only model the propagation of an input field with finite temporal width. Furthermore, the method implies that the input field is periodic [24

G. P. Agrawal, Nonlinear Fiber Optics , 3rd ed. (Academic Press, San Diego, CA, USA, 2001).

]. This means that energy going out through one side of the time window will reappear from the other side of the time window.

Due to the phase noise, the phase δϕ(0) of the input field at one edge of the time window, will differ from the phase δϕ(T _max) at the other edge of the time window. There will therefore be a large phase discontinuity in the periodic input field (Fig. 2, bottom). Unfortunately, this phase discontinuity leads to large intensity fluctuations at the edges of the time window, since phase noise is converted to intensity noise upon propagation in a dispersive medium (see, e.g., Ref. [31

W. K. Marshall, B. Crosignani, and A. Yariv, “Laser phase noise to intensity noise conversion by lowest-order group-velocity dispersion in optical fiber:exact theory,” Opt. Lett. 25, 165–167 (2000). http://www.opticsinfobase.org/abstract.cfm?URI=ol-25-3-165. [CrossRef]

]).

To avoid this problem in the numerical simulations, one can artificially replace the CW field with a broad Gaussian [23

] or super-Gaussian [10

] pulse. Here, we multiply our input field [Eq. (3)] with a super-Gaussian pulse and obtain

E (0, T) = \sqrt{P_{0}} \exp [i δ ϕ (T)] \exp [- \frac{1}{2} {(\frac{T}{T_{sG}})}^{2 m}],

(12)

where we use m=10, and T _sG is the 1/e intensity half-width of the super-Gaussian pulse, see Fig. 2(top). We show in Fig. 1(right) that multiplication with the super-Gaussian pulse practically does not affect the spectrum of the CW input field.

There are now two important deviations of the model from a real CW field: (1) T _sG is finite, meaning that due to dispersion some of the pulse energy, e.g. in the form of a soliton with a carrier frequency slightly offset from the center frequency of theCW-field, can move beyond the edges of the super-Gaussian pulse. This soliton will then no longer interact with the remaining super-Gaussian pulse, unless it exits one side of the time window and meets with the other edge of the super-Gaussian pulse. (2) The time window is finite; as explained before this means that pulse energy leaving one side of the time window reappears from the other side of the time window, as was also described in Ref. [17

Fig. 2. An illustration of the modeled periodic super-Gaussian input pulse. Top: power variation in time, bottom: phase fluctuation in time. In this example the width of the time window is 59.0 ps, the width of the super-Gaussian pulse is 2T _sG=30 ps, and Δν_FWHM=265 GHz.

Instead of an infinite CW field, we are thus actually modeling a CW laser which is turned on and off periodically. Each time the laser is turned on, the phase noise fluctuations are the same as before (see Fig. 2). The influence of these deviations decreases as the widths of the super-Gaussian pulse and the time window are increased. In the following Section, we examine the effects of these limitations closer.

It should be noted that an optical spectrum analyzer typically has a minimum integration time on the order of milliseconds. One should therefore in principle simulate a quasi-CW pulse with this duration to obtain sufficient statistical information to compare with an experimental measurement. Unfortunately, the calculation time of the numerical simulations based on the fast Fourier transform scales with N log₂(N) (see e.g. [32

W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical Recipes in C++: The Art of Scientific Computing , 2nd ed. (Cambridge University Press, Cambridge, 2002). http://www.nr.com.

]), where N is the number of computational points and is proportional to the width of the time window for a fixed temporal resolution. To decrease calculation time it is therefore common to simulate several shorter pulses, each only differing by the random initial conditions, and then average over the ensembles [16

, 17

, 23

4. Numerical results

In the following Subsection we first show that there are mainly three physical mechanisms responsible for the spectral broadening observed in this work. Although the mechanisms are well known, the role of soliton collisions in CW pumped SC generation, first described by Islam et al. for a quasi-CW pump [21

], is often overlooked (e.g. Refs. [16

, 17

]). In Subsection 4.2 we consider pumping in the fiber with relative hole size d/Λ=0.65 and find that, depending on pump spectral linewidth and the width of the time window, the red-shift of solitons can be suppressed. Subsection 4.3 then considers pumping in the fiber with relative hole size d/Λ=0.378. In this case it turns out that many solitons are formed, and that the spectral broadening is not so sensitive to the pump spectral linewidth.

4.1. Physical mechanisms

For the cases investigated in this work we have found three dominant physical mechanisms responsible for spectral broadening. One mechanism is Raman gain, which occurs when some of the power of a strong optical field (the pump) is converted to a longer wavelength (the Raman Stokes wave) [24

G. P. Agrawal, Nonlinear Fiber Optics , 3rd ed. (Academic Press, San Diego, CA, USA, 2001).

]. The energy difference between a pump photon and a Raman Stokes photon is transferred to molecular vibrations in the propagation medium. In silica, the material typically used for PCFs, most of the Raman Stokes photons are downshifted from the pump frequency by 13.2 THz. A weak optical field propagating at a frequency upshifted from the pump by 13.2 THz (the anti-Raman Stokes wave) is exponentially damped due to energy transfer to the pump. The second mechanism is MI which causes a CW pump to break up into a periodic pulse train [24

G. P. Agrawal, Nonlinear Fiber Optics , 3rd ed. (Academic Press, San Diego, CA, USA, 2001).

]. The pulse train can then evolve into a series of solitons.

The third mechanism is collision of solitons followed by soliton red-shift [21

]. In the context of optical fibers, a soliton is a pulse representing a solution to the unperturbed NLSE [24

G. P. Agrawal, Nonlinear Fiber Optics , 3rd ed. (Academic Press, San Diego, CA, USA, 2001).

, 33

V. E. Zakharov and A. B. Shabat, “Exact theory of 2-dimensional self-focusing and one-dimenstional self-modulation of waves in nonlinear media,” Sov. Phys. JETP 34, 62 (1972).

\frac{\partial E}{\partial z} = \frac{- i β_{2}}{2} \frac{\partial^{2} E}{\partial T^{2}} + i γ E {∣ E ∣}^{2} .

(13)

The unperturbed soliton propagates without changing its shape, amplitude or width even after collision with a soliton of different carrier frequency (and therefore different group velocity) [34

H. A. Haus and W. S. Wong, “Solitons in optical communications,” Rev. Mod. Phys. 68, 423–444 (1996). http://link.aps.org/abstract/RMP/v68/p423. [CrossRef]

]. Perturbations such as third-order dispersion (TOD), self-steepening, and delayed Raman response [all included in Eq. (1)] lead to, in the case of TOD, emission of radiation and position shifts during collision [35

A. Peleg, M. Chertkov, and I. Gabitov, “Interchannel interaction of optical solitons,” Phys. Rev. E 68, 026605 (2003). http://link.aps.org/abstract/PRE/v68/e026605. [CrossRef]

]. Two solitons colliding in the presence of delayed Raman response transfer energy from the soliton with the highest carrier frequency to the other soliton; this has been shown by both numerical simulations [21

,36

S. Chi and S. Wen, “Raman cross talk of soliton collision in a lossless fiber,” Opt. Lett. 14, 1216–1218 (1989). [CrossRef] [PubMed]

] and theoretical analysis [21

,37

B. A. Malomed, “Soliton-collision problem in the nonlinear Schrödinger equation with a nonlinear damping term,” Phys. Rev. A 44, 1412–1414 (1991). http://link.aps.org/abstract/PRA/v44/p1412. [CrossRef] [PubMed]

,38

Y. Chung and A. Peleg, “Strongly non-Gaussian statistics of optical soliton parameters due to collisions in the presence of delayed Raman response,” Nonlinearity 18, 1555–1574 (2005). http://stacks.iop.org/Non/18/1555. [CrossRef]

]. In general, solitons can transfer energy during collisions in non-integrable models; there is some indication that energy is preferentially transferred to the soliton with largest amplitude [39

O. Bang and M. Peyrard, “Generation of high-energy localized vibrational modes in nonlinear Klein-Gordon lattices,” Phys. Rev. E 53, 4143–4152 (1996). http://link.aps.org/abstract/PRE/v53/p4143. [CrossRef]

]. A soliton propagating in a medium with delayed Raman response is downshifted in frequency (red-shift) during propagation due to the SSFS [24

G. P. Agrawal, Nonlinear Fiber Optics , 3rd ed. (Academic Press, San Diego, CA, USA, 2001).

, 40

J. K. Lucek and K. J. Blow, “Soliton self-frequency shift in telecommunications fiber,” Phys. Rev. A 45, 6666–6674 (1992). http://dx.doi.org/10.1103/PhysRevA.45.6666. [CrossRef] [PubMed]

]. The red-shift is enhanced during collision with another soliton [36

S. Chi and S. Wen, “Raman cross talk of soliton collision in a lossless fiber,” Opt. Lett. 14, 1216–1218 (1989). [CrossRef] [PubMed]

]. It is expected that delayed Raman response is the perturbation with the most significant effect on soliton collisions [38

, 41

Y. Kodama and K. Nozaki, “Soliton interaction in optical fibers,” Opt. Lett. 12, 1038–1040 (1987). [CrossRef] [PubMed]

All of these physical mechanisms are studied in the spectrogram movie shown in Fig. 3. Initially, the CW pump breaks up into a periodic pulse train after just a few meters of propagation due to MI. To verify that the breakup is caused by MI we have in Fig. 4(left) plotted |P̃(ν)|², where P̃(ν) is the Fourier transform of the pulse power P(t)=|E(t)|². MI provides maximum gain for oscillations with a frequency given by [24

G. P. Agrawal, Nonlinear Fiber Optics , 3rd ed. (Academic Press, San Diego, CA, USA, 2001).

]

ν_{max} = \sqrt{\frac{γ P_{0}}{2 π^{2} ∣ β_{2} ∣}},

(14)

which we have indicated in Fig. 4(left) by vertical lines for the two fibers considered here. It is seen that oscillations at the MI frequency are significant for both fibers and are therefore the main cause of pulse train formation.

As seen from the spectrogram movie, Fig. 3, the periodic pulse train evolves into solitons. The number of solitons can be estimated as the ratio between the super-Gaussian pulse width and the period T _m of the periodic pulse train: 2T _sG/T _m=2T _sG×ν_max=30 ps×0.75 THz≈23. This corresponds well with the number of solitons in the top right of Fig. 3. The solitons have slightly different group velocities for two reasons. First, the carrier frequency of the CW pump fluctuates in time due to the phase noise. Second, they have slightly different temporal widths, which means that the rate of SSFS is also different. As mentioned previously, SSFS causes the carrier frequencyν₀ of a soliton to red-shift as it propagates along the fiber. The rate of red-shift dν₀/dz can be expressed analytically in two limits. For a soliton temporal width T ₀ much longer than 2π/Ω_R≈76 fs, where Ω_R=2π×13.2 THz is the angular frequency shift for maximum Raman gain in silica [24

G. P. Agrawal, Nonlinear Fiber Optics , 3rd ed. (Academic Press, San Diego, CA, USA, 2001).

, 42

J. P. Gordon, “Theory of the soliton self-frequency shift,” Opt. Lett. 11, 662–664 (1986). [CrossRef] [PubMed]

, 43

J. Herrmann and A. Nazarkin, “Soliton self-frequency shift for pulses with a duration less than the period of molecular oscillations,” Opt. Lett. 19, 2065–2067 (1994). [CrossRef] [PubMed]

]

\frac{d ν_{0}}{d z} = - \frac{8 ∣ β_{2} ∣ T_{R}}{2 π 15 T_{0}^{4}}, T_{0} ≫ 76 fs,

(15)

where T _R≈3 fs is related to the slope of the Raman gain spectrum [24

G. P. Agrawal, Nonlinear Fiber Optics , 3rd ed. (Academic Press, San Diego, CA, USA, 2001).

]. For a soliton duration shorter than 76 fs the red-shift rate is given by [43

J. Herrmann and A. Nazarkin, “Soliton self-frequency shift for pulses with a duration less than the period of molecular oscillations,” Opt. Lett. 19, 2065–2067 (1994). [CrossRef] [PubMed]

]

\frac{d ν_{0}}{d z} = - \frac{0.09 ∣ β_{2} ∣ Ω_{R}^{2}}{2 π T_{0}}, T_{0} \underset{\approx}{<} 76 fs .

(16)

Fig. 3. Spectrogram movie for one simulation of propagation in the d/Λ=0.65 fiber. Input super-Gaussian pulse width 2T _sG=30 ps, Δν_FWHM=30 GHz, time window T _max=59 ps. Note that the color scale changes during propagation, as some of the solitons acquire higher peak power. The horizontal white line on the bottom Figs. indicates the Raman Stokes wavelength at 1116 nm. pitch1p72 dL0p65 30GHz.avi (1.93 MB).

Fig. 4. Left: Frequency spectrum of pulse power P(t) at z=9 m calculated for one simulation for the d/Λ=0.65 fiber (blue, solid) and one simulation for the d/Λ=0.378 fiber (green, dashed). The vertical lines indicate the frequency with maximum MI gain in the corresponding fibers. Right: Maximum peak power along the fiber length for the same simulation as shown in Fig. 3. The horizontal black line at 152 Windicates the minimum peak power required for a soliton to red-shift to 1116 nm after 74 m of propagation.

Because of the slightly different group velocities, and perhaps to some degree also by soliton attraction/repulsion [24

G. P. Agrawal, Nonlinear Fiber Optics , 3rd ed. (Academic Press, San Diego, CA, USA, 2001).

], the solitons undergo collisions as they propagate. As stated above, each collision leads to a fractional amount of energy transfer to the most red-shifted soliton. A soliton that gains energy from a collision increases its peak power and decreases its temporal width to maintain its fundamental soliton shape [24

G. P. Agrawal, Nonlinear Fiber Optics , 3rd ed. (Academic Press, San Diego, CA, USA, 2001).

]. The decrease in temporal width leads to an increased red-shift rate, which in turn modifies the group velocity [21

]. This soliton will therefore also have a greater probability of colliding with other solitons. The mechanism is therefore self-amplifying. It is seen in Fig. 3 how one particular soliton gradually gains energy and obtains such a short temporal width that it red-shifts significantly faster than the other solitons. Figure 4(right) shows the maximum peak power of the propagated field along the fiber, for the same simulation as in Fig. 3. It is known from Fig. 3 that the soliton formation from MI is complete after less than ~10 m of propagation. It is seen from Fig. 4 (right) that the maximum peak power continues to increase along the fibre after a propagation length of z=10 m because of the energy transfer during soliton collisions.

The temporal width of the MI-generated solitons can be estimated from the MI oscillation period T _m=1/ν_max. Thus, the duration of the generated solitons is 1/(0.75 THz)≈1.3 ps and 1/(2.5 THz)≈0.4 ps, in the d/Λ=0.65 and the d/Λ=0.378 fiber, respectively. Since the red-shift rate for such soliton durations is proportional to

T_{0}^{- 4}

[Eq. (15)], even a small decrease in pulse width could lead to a significant increase in red-shift rate.

To clarify that energy transfer during soliton collisions is necessary for the large observed red-shift, and that the resulting spectrum can not be explained by MI and SSFS alone, the following analysis is performed. The most red-shifted soliton in Fig. 3 has shifted to 1161 nm at z=74 m. This corresponds to a shift of Δν₀=-23.5 THz from the input pump over Δz=74 m. An estimate for the red-shift using Eq. (15) gives Δν₀≈-8|β ₂|T _RΔz/(2π15

T_{0}^{4}

)≈-3×10^-4 THz, where the soliton width T ₀ was approximated by the pulse train period T _MI=1.3 ps. β ₂(λ) changes only slightly from -4.38×10^-26 s²/m at λ0=1064 nm to -6.14×10^-26 s²/m at λ=1161 nm; β ₂(λ) is therefore approximated by β ₂(λ ₀)=β ₂. If there is no energy transfer during soliton collisions, the change in T ₀ during red-shift can also be neglected. It is thus clear that the solitons generated directly from MI, in the cases investigated here, are too temporally long to make a notable red-shift. This is in accordance with the finding by Golovchenko et al., who found that for ν_max≪ν_R=13.2 THz the soliton red-shift is insignificant [20

To support this conclusion further, Fig. 5(left) shows the red-shift rate dν₀/dz in the two regimes where Eqs. (15–16) are valid. The red-shift rate has been scaled by -1/|β ₂| for generality. The observed shift of -23.5 THz corresponds to a mean scaled red-shift rate of (-1/β̄ ₂)Δν₀/Δz≈7.3×10³⁶ s^-3. It is seen from Fig. 5(left) that the scaled red-shift rate for a soliton width T ₀ ≈ 1 ps is four orders of magnitude smaller than the observed mean red-shift rate. This shows that even though there is a

T_{0}^{- 4}

dependence, it is not enough for the MI-generated soliton width to change ‘slightly’, if the soliton is to make the observed red-shift. The soliton width has to be decreased by an order of magnitude, as seen from Fig. 5(left). This is achieved through the energy transfer during soliton collisions.

Fig. 5. Left: Red-shift rate dν₀/dz scaled by -1/|β ₂| in the two regimes where analytical expressions are available: T ₀~76 fs [Eq. (16)] (blue, solid) and T ₀ 76 fs [Eq. (15)] (green, dashed). Right: Illustration of a super-Gaussian pulse with width 2T _sG=15 ps in a time domain of T _max=29.5 ps (top, case A), 2T _sG=15 ps and T _max=59 ps (middle, case B), and 2T _sG=30 ps and T _max=59 ps (bottom, case C).

We note that one can also see soliton formation from MI followed by the emergence of a few quickly red-shifting high-power solitons in the simulations performed in Refs. [16

, 17

]. Estimates of the red-shift rate for the parameters used in those papers indicate that solitons formed directly from MI could not explain the observed red-shift nor the short temporal width of the observed solitons. We therefore believe that these solitons were also the result of energy transfer during collisions. The authors seem to have overlooked this possibly because they did not consider in detail the formation of the high-power solitons using e.g. spectrograms.

The above conclusion on the necessity of soliton collisions to explain the observed red-shift of a few individual solitons is based on the parameters used in this work and in Refs. [16

,17

]. It should be emphasized that, as shown by Eq. (14), the temporal width of MI-generated solitons can be made sufficiently small to obtain a significant red-shift without soliton collisions by increasing γ and/or P ₀, and/or by decreasing |β ₂|. It was, e.g., shown by Golovchenko et al. that for sufficiently large pump power, the MI-generated solitons are short enough for all of them to red-shift [20

4.2. Formation of few or no red-shifting solitons

We have found that for the fiber with Λ=1.72 µm and d/Λ=0.65, only few or even no red-shifting solitons are formed. We treat the cases of a CW pump with 30 GHz or 265 GHz linewidth separately.

4.2.1. 30 GHz linewidth

First, we consider case A illustrated in the top of Fig. 5(right): a super-Gaussian pulse width of 15 ps, and a time window of T _max=29.5 ps. 2¹⁴ points were used with a time resolution of Δt=1.8 fs. This results in a frequency window of 1/Δt≈556 THz, and a frequency resolution of 1/T _max≈34 GHz. We performed 10 simulations with different seeds for the random number generator used to simulate the phase noise. As shown in Fig. 6(left), all the resulting spectra have a peak at 1116 nm and a dip at 1013 nm. These features both correspond very well with a frequency shift from the pump by 13.2 THz, which is the peak frequency shift for Raman gain [24

G. P. Agrawal, Nonlinear Fiber Optics , 3rd ed. (Academic Press, San Diego, CA, USA, 2001).

]. Thus, Raman gain is the dominant physical mechanism responsible for the resulting spectrum in this case.

Fig. 6. Spectra calculated for 10 ensembles differing by the random number seed for the phase noise. Λ=1.72 µm, d/Λ=0.65, and the fibre length L=80 m. Linewidth is 30 GHz, input super-Gaussian pulse width 2T _sG=15 ps. Left: 2¹⁴ points used, T _max=29.5 ps, smoothed over 8 points. Right: 2¹⁵ points used, T _max=59 ps, smoothed over 16 points.

We now consider case B illustrated in the middle of Fig. 5(right) by doubling the width of the time window to 59 ps, but keeping the width of the input super-Gaussian pulse the same as in case A. Since we keep the time resolution constant, we now use 2¹⁵ points. This means that the frequency resolution is halved from 34 GHz to 17 GHz. We assume that since the frequency resolution is already high, this has negligible effect on the results. As seen in Fig. 6(right) the peak and the dip resulting from Raman gain are now less prominent, indicating that less energy is transferred via the Raman effect. Also, in 3 of the 10 simulations a red-shifting soliton is seen with a center wavelength larger than 1100 nm.

The difference in resulting spectra between cases A and B can be explained as follows. We have calculated the difference in β ₁=1/v_g between the pump (1064 nm) and the Raman Stokes wave (1116.3 nm) to be 4.0 ps/m (for the d/Λ=0.65 fiber). This means that pulse energy transferred from the pump to the Raman Stokes wave will separate from the pump in the time domain at a rate of 4.0 ps for each meter of propagation. The Raman Stokes wave is thus spread out over the time window after just a few meters of propagation. This can be seen in the bottom of Fig. 3. The transfer of energy from the pump to the Raman Stokes wave requires a temporal overlap. When we increase the time window from case A to case B, the Raman Stokes wave is spread over a larger time domain, thus decreasing the temporal overlap with the pump. The transfer of energy from the Raman anti-Stokes wave to the pump, and from the pump to the Raman Stokes wave therefore becomes less efficient for case B. Since less energy is removed from the pump wavelength there remains more energy in the solitons undergoing collisions. There is therefore a greater probability of a high-energy soliton being formed, with such a short temporal width that it red-shifts far-away from the pump.

Note that we simulated a super-Gaussian input pulse with the same width in cases A and B, and yet the resulting spectra are very different. This important result shows that the choice of time window width requires careful consideration.

To consider case C (bottom of Fig. 5, right) we keep the same time window as in case B, but now double the width of the super-Gaussian input pulse to 2T _sG=30 ps. The ratio of super-Gaussian pulse width and time window width, 2T _sG/T _max, is now the same as in case A. From the physical explanation above, the Raman energy transfer from e.g. the pump to the Raman Stokes wave should be as effective as in case A. This is because the temporal overlap between the Raman Stokes wave and the pump is the same in both cases A and C, if we assume that the Raman Stokes wave is quickly almost uniformly spread over the entire time window. Indeed, as seen in Fig. 7(left), the Raman anti-Stokes dip and the Raman Stokes peak have practically the same power spectral density as for case A [Fig. 6(left)]. It is also seen in Fig. 3 how the Raman Stokes wave is uniformly spread across the time window. In case A the Raman energy transfer from the pump was so severe that no red-shifting solitons were observed. However, the doubling of the super-Gaussian input pulse width, compared to case A, means that a larger number of colliding solitons will be formed from MI. This should result in a greater probability for the formation of a red-shifting soliton. As expected, it is seen in Fig. 7(left) that in 4 out of the 10 ensembles a soliton was observed to red-shift beyond 1100 nm.

Fig. 7. Spectra calculated for 10 ensembles differing by the random number seed for the phase noise generation. All for Λ=1.72 µm, d/Λ=0.65, and a fibre length of 80 m. 215 points used, 2T _sG=30 ps, T _max=59 ps. Left: Spectral linewidth Δν_FWHM=30 GHz. Right: Δν_FWHM=265 GHz. Smoothed over 16 points. The thin black line indicates the input spectrum for one of the simulations.

4.2.2. 265 GHz linewidth

We again use the computational parameters for case C (2T _sG=30 ps, T _max=59 ps, 2¹⁵ points) but now simulate a spectral FWHM linewidth of Δν_FWHM=265 GHz. The resulting spectra are seen in Fig. 7(right).

The only simulation difference between Fig. 7(right) and Fig. 7(left) is the spectral linewidth of the quasi-CW input. None of the simulations with 265 GHz linewidth show distinct red-shifting solitons, whereas the simulations with 30 GHz linewidth showed a red-shifting soliton in 4 out of 10 simulations.

We suggest two possible explanations for this difference. (1) The larger pump linewidth leads to a larger spread in center frequencies of solitons as they form from the pump. Since the energy transfer between two colliding solitons decreases with increasing separation in carrier frequency [36

S. Chi and S. Wen, “Raman cross talk of soliton collision in a lossless fiber,” Opt. Lett. 14, 1216–1218 (1989). [CrossRef] [PubMed]

], this leads to less overall energy transfer during soliton collisions. This could then hinder the buildup of a quickly red-shifting soliton. (2) As the linewidth of the Lorentzian power spectrum is increased, the pump power spectral density S(λ=1064 nm) is decreased, see Fig. 8(left). At the same time, there is more spectral power S(λ=1116 nm) available at the Raman Stokes wavelength. The reduction in power spectral density at the pump wavelength could also hinder the buildup of a quickly red-shifting soliton. In Fig. 8(right) we have plotted the power within a bandwidth of ~4 nm centered at the pump wavelength and the Raman Stokes wavelength, respectively, during propagation along the fiber. It is seen that for a linewidth of 265 GHz there is slightly less power at the pump wavelength, and significantly more power at the Raman Stokes wavelength, compared to the simulation with a linewidth of 30 GHz.

We have not been able to determine the relative importance of these two mechanisms in hindering the buildup of a quickly red-shifting soliton.

The simulations in Figs. 6(right) and 7(left) showed a quickly red-shifting soliton in only some of the simulations. The solitons also did not undergo the same amount of red-shift. For example, one simulation in Fig. 7(left) resulted in a soliton at ~1170 nm and another simulation had a soliton at ~1270 nm. We show in Fig. 9 how averaging over the ensembles smooths the peaks from individual solitons. Experimental spectra were presented for similar pumping conditions in Refs. [11

A. V. Avdokhin, S. V. Popov, and J. R. Taylor, “Continuous-wave, high-power, Raman continuum generation in holey fibers,” Opt. Lett. 28, 1353–1355 (2003). [CrossRef] [PubMed]

, 12

]. The experimental spectra showed a large pump residual at 1064 nm, a Raman Stokes peak at 1116 nm about 11 dB below the pump residual, and a relatively flat continuum spanning from ~1150 to ~1350 nm about 14 dB below the pump residual. We compare this with the ensemble average in Fig. 9(right), corresponding to case C, because this should be a better approximation to a real CW input than cases A and B. In our simulations the pump residual and Raman Stokes peak are at the same spectral power level. This could be caused by the fact that our simulations assume a polarized input pump and only consider propagation in one polarization axis, whereas the experiments were done with a non-polarized pump laser. Instead of a flat continuum from 1150 nm to 1350 nm our simulations show three distinct peaks in the same wavelength region. This occurs because we have only averaged over 10 ensembles; we know that the red-shift of each soliton is determined by the random initial conditions, so increasing the number of averaged ensembles would smooth the soliton spectra into a continuum [21

]. Ideally, for comparison with the measurement of an optical spectrum analyzer with 1 ms integration time, one would need to average over ~10⁷ ensembles, if each ensemble corresponds to a quasi-CW pulse of 30 ps duration. In practice, averaging over as few as 100 ensembles is found to give reasonable agreement with experimental measurements [16

, 17

Fig. 8. Left: Lorentzian shaped power spectrum for 30 GHz (solid, blue) and 265 GHz linewidth (dashed, green). Right: Power evolution along the fibre at the pump wavelength and the Raman Stokes wavelength of 1116 nm. One simulation for a linewidth of 30 GHz (solid, blue) and one simulation for a 265 GHz linewidth (dashed, green).

4.3. Formation of several red-shifting solitons

We now consider the fiber with the structural parameters Λ=1.72 µm and d/Λ=0.378. This fiber has zero-dispersionwavelengths (ZDWs) at ~1025 and 1281 nm, as seen from Fig. 1(left). It is also seen that the dispersion at the pump wavelength is ~20 times smaller for this fiber, compared to the fiber with d/Λ=0.65. The effective area A _eff is ~2 times larger, and using Eq. (14) this results in a MI maximum gain frequency of ν _max=2.5 THz, as was also indicated in Fig. 4(left). This means that approximately 3 times more solitons should be formed in this fiber than in the fiber with d/Λ=0.65. This is also what we have observed in spectrograms. The larger number of solitons leads to a larger number of soliton collisions which in turn increases the probability of the creation of a quickly red-shifting soliton. Furthermore, the solitons are not only more abundant, but also have a shorter temporal width giving them a larger red-shift rate (Fig. 5, right). The resulting spectra are seen in Fig. 10.

Fig. 9. Left: Ensemble average for the same parameters as in Fig. 6(right). Right: Ensemble average for the same parameters as in Fig. 7(left). The averaged input spectrum is indicated as a thin black line.

Fig. 10. Calculated spectra at z=64 m for the Λ=1.72 µm, d/Λ=0.378 fiber, super-Gaussian input pulse width 2T _sG=30 ps, time window T _max=59 ps, 2¹⁵ points. Left: 30 GHz linewidth. Right: 265 GHz linewidth. 10 ensembles shown in each case. The thin black line indicates the input spectrum for one of the ensembles.

It is seen that the large number of red-shifting solitons results in a relatively smooth spectral broadening on the red side of the pump. The spectrum has achieved the large broadening at a shorter fiber distance (z=64 m) than for the d/Λ=0.65 µm fiber. We also observe the generation of dispersive waves in the normal dispersion region [4

, 5

] at ~857 nm and ~1376 nm.

The formation of many solitons due to the lower dispersion also means that the input linewidth is less critical than for the d/Λ=0.65 fiber. As is seen in Fig. 10, the simulations with 30 GHz and 265 GHz linewidth show much more resemblance to each other than the corresponding simulations for the d/Λ=0.65 fiber (Fig. 7).

5. Conclusion

Recently, it was claimed that solitons formed directly from modulation instability could red-shift far away from the CW pump [16

, 17

]. We found, for the cases investigated here and in Refs. [16

, 17

], that the solitons require a shorter temporal width to make a significant red-shift, and showed that this can be achieved by energy transfer during soliton collisions.

It was shown how the time window width used in the calculations can affect the resulting spectrum and must be chosen carefully.

We demonstrated that increasing the spectral linewidth of the CW pump laser can hinder the formation of quickly red-shifting solitons. This indicates that narrow linewidth pump lasers should be selected for generating a broader supercontinuum. The number of solitons formed from a quasi-CW pulse and their red-shift rate can be increased by choosing a fiber with lower dispersion at the pump wavelength. This reduces the required fiber length and allows the use of a pump laser with broader linewidth while still obtaining quickly red-shifting solitons.

Acknowledgments

The authors would like to thank Jes Broeng and Kim P. Hansen from Crystal Fibre A/S, and Peter E. Andersen, Steen G. Hanson, and Lars Thrane from Risø National Laboratory for helpful discussions. This work was partly supported by the EU project NANO UB-SOURCES, contract no. 017128.

References and links

1.	R. R. Alfano, ed., The Supercontinuum Laser Source: Fundamentals With Updated References , 2nd ed. (Springer-Verlag, New York, 2005). ISBN 0-387-24504-9.
2.	J. K. Ranka, R. S. Windeler, and A. J. Stentz, “Visible continuum generation in air-silica microstructure optical fibers with anomalous dispersion at 800 nm,” Opt. Lett. 25, 25–27 (2000). [CrossRef]
3.	K. M. Hilligsøe, T. V. Andersen, H. N. Paulsen, C. K. Nielsen, K. Mølmer, S. Keiding, R. Kristiansen, K. P. Hansen, and J. J. Larsen, “Supercontinuum generation in a photonic crystal fiber with two zero dispersion wavelengths,” Opt. Express 12, 1045–1054 (2004). http://www.opticsexpress.org/abstract.cfm?URI=OPEX-12-6-1045. [CrossRef] [PubMed]
4.	G. Genty, M. Lehtonen, H. Ludvigsen, and M. Kaivola, “Enhanced bandwidth of supercontinuum generated in microstructured fibers,” Opt. Express 12, 3471–3480 (2004). http://www.opticsexpress.org/abstract.cfm?URI=OPEX-12-15-3471. [CrossRef] [PubMed]
5.	M. H. Frosz, P. Falk, and O. Bang, “The role of the second zero-dispersion wavelength in generation of supercontinua and bright-bright soliton-pairs across the zero-dispersion wavelength,” Opt. Express 13, 6181–6192 (2005). http://www.opticsexpress.org/abstract.cfm?URI=OPEX-13-16-6181. [CrossRef] [PubMed]
6.	T. A. Birks, W. J. Wadsworth, and P. St. J. Russell, “Supercontinuum generation in tapered fibers,” Opt. Lett. 25, 1415–1417 (2000). [CrossRef]
7.	J. M. Dudley, L. Provino, N. Grossard, H. Maillotte, R. S. Windeler, B. J. Eggleton, and S. Coen, “Supercontinuum generation in air-silica microstructured fibers with nanosecond and femtosecond pulse pumping,” J. Opt. Soc. Am. B 19, 765–771 (2002). [CrossRef]
8.	S. Coen, A. H. L. Chau, R. Leonhardt, J. D. Harvey, J. C. Knight, W. J. Wadsworth, and P. St. J. Russell, “Supercontinuum generation by stimulated Raman scattering and parametric four-wave mixing in photonic crystal fibers,” J. Opt. Soc. Am. B 19, 753–764 (2002). [CrossRef]
9.	M. Prabhu, N. S. Kim, and K. Ueda, “Ultra-Broadband CW Supercontinuum Generation Centered at 1483.4 nm from Brillouin/Raman Fiber Laser,” Jpn. J. Appl. Phys. 39, L291–L293 (2000). http://dx.doi.org/10.1143/JJAP.39.L291. [CrossRef]
10.	J. W. Nicholson, A. K. Abeeluck, C. Headley, M. F. Yan, and C. G. Jørgensen, “Pulsed and continuous-wave supercontinuum generation in highly nonlinear, dispersion-shifted fibers,” Appl. Phys. B 77, 211–218 (2003). http://dx.doi.org/10.1007/s00340-003-1201-z. [CrossRef]
11.	A. V. Avdokhin, S. V. Popov, and J. R. Taylor, “Continuous-wave, high-power, Raman continuum generation in holey fibers,” Opt. Lett. 28, 1353–1355 (2003). [CrossRef] [PubMed]
12.	P.-L. Hsiung, Y. Chen, T. Ko, J. G. Fujimoto, C. J. S. de Matos, S. V. Popov, J. R. Taylor, and V. P. Gapontsev, “Optical coherence tomography using a continuous-wave, high-power, Raman continuum light source,” Opt. Express 12, 5287–5295 (2004). http://www.opticsinfobase.org/abstract.cfm?URI=oe-12-22-5287. [CrossRef] [PubMed]
13.	C. J. S. de Matos, S. V. Popov, and J. R. Taylor, “Temporal and noise characteristics of continuous-wave-pumped continuum generation in holey fibers around 1300 nm,” Appl. Phys. Lett. 85, 2706–2708 (2004). [CrossRef]
14.	J. H. Lee, Y.-G. Han, and S. B. Lee, “Experimental study on seed light source coherence dependence of continuous-wave supercontinuum performance,” Opt. Express 14, 3443–3452 (2006). http://www.opticsinfobase.org/abstract.cfm?URI=oe-14-8-3443. [CrossRef] [PubMed]
15.	A. Mussot, E. Lantz, H. Maillotte, T. Sylvestre, C. Finot, and S. Pitois, “Spectral broadening of a partially coherent CW laser beam in single-mode optical fibers,” Opt. Express 12, 2838–2843 (2004). http://www.opticsexpress.org/abstract.cfm?URI=OPEX-12-13-2838. [CrossRef] [PubMed]
16.	S. M. Kobtsev and S. V. Smirnov, “Modelling of high-power supercontinuum generation in highly nonlinear, dispersion shifted fibers at CW pump,” Opt. Express 13, 6912–6918 (2005). http://www.opticsexpress.org/abstract.cfm?URI=OPEX-13-18-6912. [CrossRef] [PubMed]
17.	F. Vanholsbeeck, S. Martin-Lopez, M. González-Herráez, and S. Coen, “The role of pump incoherence in continuous-wave supercontinuum generation,” Opt. Express 13, 6615–6623 (2005). http://www.opticsexpress.org/abstract.cfm?URI=OPEX-13-17-6615. [CrossRef] [PubMed]
18.	A. K. Abeeluck and C. Headley, “Continuous-wave pumping in the anomalous- and normal-dispersion regimes of nonlinear fibers for supercontinuum generation,” Opt. Lett. 30, 61–63 (2005). [CrossRef] [PubMed]
19.	B. Barviau, S. Randoux, and P. Suret, “Spectral broadening of a multimode continous-wave optical field propagating in the normal dispersion regime of a fiber,” Opt. Lett. 31, 1696–1698 (2006). [CrossRef] [PubMed]
20.	E. A. Golovchenko, P. V. Mamyshev, A. N. Pilipetskii, and E. M. Dianov, “Numerical analysis of the Raman spectrum evolution and soliton pulse generation in single-mode fibers,” J. Opt. Soc. Am. B 8, 1626–1632 (1991). http://www.opticsinfobase.org/abstract.cfm?URI=josab-8-8-1626. [CrossRef]
21.	M. N. Islam, G. Sucha, I. Bar-Joseph, M. Wegener, J. P. Gordon, and D. S. Chemla, “Femtosecond distributed soliton spectrum in fibers,” J. Opt. Soc. Am. B 6, 1149–1158 (1989). http://www.opticsinfobase.org/abstract.cfm?URI=josab-6-6-1149. [CrossRef]
22.	R. G. Smith, “Optical power handling capacity of low loss optical fibers as determined by stimulated Raman and Brillouin scattering,” Appl. Opt. 11, 2489–2494 (1972). [CrossRef] [PubMed]
23.	S. B. Cavalcanti, G. P. Agrawal, and M. Yu, “Noise amplification in dispersive nonlinear media,” Phys. Rev. A 51, 4086–4092 (1995). http://dx.doi.org/10.1103/PhysRevA.51.4086. [CrossRef] [PubMed]
24.	G. P. Agrawal, Nonlinear Fiber Optics , 3rd ed. (Academic Press, San Diego, CA, USA, 2001).
25.	P. L. François, “Nonlinear propagation of ultrashort pulses in optical fibers: total field formulation in the frequency domain,” J. Opt. Soc. Am. B 8, 276–293 (1991). http://www.opticsinfobase.org/abstract.cfm?URI=josab-8-2-276. [CrossRef]
26.	S. G. Johnson and J. D. Joannopoulos, “Block-Iterative Frequency-Domain Methods for Maxwell’s Equations in a Planewave Basis,” Opt. Express 8, 173–190 (2001). http://www.opticsexpress.org/abstract.cfm?URI=OPEX-8-3-173. [CrossRef] [PubMed]
27.	K. P. Hansen, Crystal Fibre A/S, Personal communication.
28.	J. Lægsgaard, N. A. Mortensen, and A. Bjarklev, “Mode areas and field-energy distribution in honeycomb photonic bandgap fibers,” J. Opt. Soc. Am. B 20, 2037–2045 (2003). [CrossRef]
29.	K. J. Blow and D. Wood, “Theoretical description of transient stimulated Raman scattering in optical fibers,” IEEE J. Quantum Electron. 25, 2665–2673 (1989). [CrossRef]
30.	J. W. Goodman, Statistical Optics (John Wiley & Sons Inc., 2000). ISBN 0471399167.
31.	W. K. Marshall, B. Crosignani, and A. Yariv, “Laser phase noise to intensity noise conversion by lowest-order group-velocity dispersion in optical fiber:exact theory,” Opt. Lett. 25, 165–167 (2000). http://www.opticsinfobase.org/abstract.cfm?URI=ol-25-3-165. [CrossRef]
32.	W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical Recipes in C++: The Art of Scientific Computing , 2nd ed. (Cambridge University Press, Cambridge, 2002). http://www.nr.com.
33.	V. E. Zakharov and A. B. Shabat, “Exact theory of 2-dimensional self-focusing and one-dimenstional self-modulation of waves in nonlinear media,” Sov. Phys. JETP 34, 62 (1972).
34.	H. A. Haus and W. S. Wong, “Solitons in optical communications,” Rev. Mod. Phys. 68, 423–444 (1996). http://link.aps.org/abstract/RMP/v68/p423. [CrossRef]
35.	A. Peleg, M. Chertkov, and I. Gabitov, “Interchannel interaction of optical solitons,” Phys. Rev. E 68, 026605 (2003). http://link.aps.org/abstract/PRE/v68/e026605. [CrossRef]
36.	S. Chi and S. Wen, “Raman cross talk of soliton collision in a lossless fiber,” Opt. Lett. 14, 1216–1218 (1989). [CrossRef] [PubMed]
37.	B. A. Malomed, “Soliton-collision problem in the nonlinear Schrödinger equation with a nonlinear damping term,” Phys. Rev. A 44, 1412–1414 (1991). http://link.aps.org/abstract/PRA/v44/p1412. [CrossRef] [PubMed]
38.	Y. Chung and A. Peleg, “Strongly non-Gaussian statistics of optical soliton parameters due to collisions in the presence of delayed Raman response,” Nonlinearity 18, 1555–1574 (2005). http://stacks.iop.org/Non/18/1555. [CrossRef]
39.	O. Bang and M. Peyrard, “Generation of high-energy localized vibrational modes in nonlinear Klein-Gordon lattices,” Phys. Rev. E 53, 4143–4152 (1996). http://link.aps.org/abstract/PRE/v53/p4143. [CrossRef]
40.	J. K. Lucek and K. J. Blow, “Soliton self-frequency shift in telecommunications fiber,” Phys. Rev. A 45, 6666–6674 (1992). http://dx.doi.org/10.1103/PhysRevA.45.6666. [CrossRef] [PubMed]
41.	Y. Kodama and K. Nozaki, “Soliton interaction in optical fibers,” Opt. Lett. 12, 1038–1040 (1987). [CrossRef] [PubMed]
42.	J. P. Gordon, “Theory of the soliton self-frequency shift,” Opt. Lett. 11, 662–664 (1986). [CrossRef] [PubMed]
43.	J. Herrmann and A. Nazarkin, “Soliton self-frequency shift for pulses with a duration less than the period of molecular oscillations,” Opt. Lett. 19, 2065–2067 (1994). [CrossRef] [PubMed]

OCIS Codes
(190.4370) Nonlinear optics : Nonlinear optics, fibers
(190.5530) Nonlinear optics : Pulse propagation and temporal solitons
(190.5650) Nonlinear optics : Raman effect

ToC Category:
Nonlinear Optics

History
Original Manuscript: July 5, 2006
Revised Manuscript: September 13, 2006
Manuscript Accepted: September 13, 2006
Published: October 2, 2006

Citation
Michael H. Frosz, Ole Bang, and Anders Bjarklev, "Soliton collision and Raman gain regimes in continuous-wave pumped supercontinuum generation," Opt. Express 14, 9391-9407 (2006)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-14-20-9391

Sort: Author | Year | Journal | Reset

References

R. R. Alfano, ed., The Supercontinuum Laser Source: Fundamentals With Updated References, 2nd ed. (Springer-Verlag, New York, 2005).
J. K. Ranka, R. S. Windeler, and A. J. Stentz, "Visible continuum generation in air-silica microstructure optical fibers with anomalous dispersion at 800 nm," Opt. Lett. 25, 25-27 (2000). [CrossRef]
K. M. Hilligsøe, T. V. Andersen, H. N. Paulsen, C. K. Nielsen, K. Mølmer, S. Keiding,R. Kristiansen, K. P. Hansen, and J. J. Larsen, "Supercontinuum generation in a photoniccrystal fiber with two zero dispersion wavelengths," Opt. Express 12, 1045-1054 (2004). [CrossRef] [PubMed]
G. Genty, M. Lehtonen, H. Ludvigsen, and M. Kaivola, "Enhanced bandwidth of supercontinuum generated in microstructured fibers," Opt. Express 12, 3471-3480 (2004). [CrossRef] [PubMed]
M. H. Frosz, P. Falk, and O. Bang, "The role of the second zero-dispersion wavelength in generation of supercontinua and bright-bright soliton-pairs across the zero-dispersion wavelength," Opt. Express 13, 6181-6192 (2005). [CrossRef] [PubMed]
T. A. Birks, W. J. Wadsworth, and P. St. J. Russell, "Supercontinuum generation in tapered fibers," Opt. Lett. 25, 1415-1417 (2000). [CrossRef]
J. M. Dudley, L. Provino, N. Grossard, H. Maillotte, R. S. Windeler, B. J. Eggleton, and S. Coen, "Supercontinuum generation in air-silica microstructured fibers with nanosecond and femtosecond pulse pumping," J. Opt. Soc. Am. B 19, 765-771 (2002). [CrossRef]
S. Coen, A. H. L. Chau, R. Leonhardt, J. D. Harvey, J. C. Knight, W. J. Wadsworth, and P. St. J. Russell, "Supercontinuum generation by stimulated Raman scattering and parametric four-wave mixing in photonic crystal fibers," J. Opt. Soc. Am. B 19, 753-764 (2002). [CrossRef]
M. Prabhu, N. S. Kim, and K. Ueda, "Ultra-Broadband CW Supercontinuum Generation Centered at 1483.4 nm from Brillouin/Raman Fiber Laser," Jpn. J. Appl. Phys. 39, L291-L293 (2000). http://dx.doi.org/10.1143/JJAP.39.L291. [CrossRef]
J. W. Nicholson, A. K. Abeeluck, C. Headley, M. F. Yan, and C. G. Jørgensen, "Pulsed and continuous-wave supercontinuum generation in highly nonlinear, dispersion-shifted fibers," Appl. Phys. B 77, 211-218 (2003). http://dx.doi.org/10.1007/s00340-003-1201-z. [CrossRef]
A. V. Avdokhin, S. V. Popov, and J. R. Taylor, "Continuous-wave, high-power, Raman continuum generation in holey fibers," Opt. Lett. 28, 1353-1355 (2003). [CrossRef] [PubMed]
P.-L. Hsiung, Y. Chen, T. Ko, J. G. Fujimoto, C. J. S. de Matos, S. V. Popov, J. R. Taylor, and V. P. Gapontsev, "Optical coherence tomography using a continuous-wave, high-power, Raman continuum light source," Opt. Express 12, 5287-5295 (2004). http://www.opticsinfobase.org/abstract.cfm?URI=oe-12-22-5287. [CrossRef] [PubMed]
C. J. S. de Matos, S. V. Popov, and J. R. Taylor, "Temporal and noise characteristics of continuous-wave-pumped continuum generation in holey fibers around 1300 nm," Appl. Phys. Lett. 85, 2706-2708 (2004). [CrossRef]
J. H. Lee, Y.-G. Han, and S. B. Lee, "Experimental study on seed light source coherence dependence of continuous-wave supercontinuum performance," Opt. Express 14, 3443-3452 (2006). http://www.opticsinfobase.org/abstract.cfm?URI=oe-14-8-3443. [CrossRef] [PubMed]
A. Mussot, E. Lantz, H. Maillotte, T. Sylvestre, C. Finot, and S. Pitois, "Spectral broadening of a partially coherent CW laser beam in single-mode optical fibers," Opt. Express 12, 2838-2843 (2004). http://www.opticsexpress.org/abstract.cfm?URI=OPEX-12-13-2838. [CrossRef] [PubMed]
S. M. Kobtsev and S. V. Smirnov, "Modelling of high-power supercontinuum generation in highly nonlinear, dispersion shifted fibers at CW pump," Opt. Express 13, 6912-6918 (2005). [CrossRef] [PubMed]
F. Vanholsbeeck, S. Martin-Lopez, M. González-Herráez, and S. Coen, "The role of pump incoherence in continuous-wave supercontinuum generation," Opt. Express 13, 6615-6623 (2005). [CrossRef] [PubMed]
A. K. Abeeluck and C. Headley, "Continuous-wave pumping in the anomalous- and normal-dispersion regimes of nonlinear fibers for supercontinuum generation," Opt. Lett. 30, 61-63 (2005). [CrossRef] [PubMed]
B. Barviau, S. Randoux, and P. Suret, "Spectral broadening of a multimode continous-wave optical field propagating in the normal dispersion regime of a fiber," Opt. Lett. 31, 1696-1698 (2006). [CrossRef] [PubMed]
E. A. Golovchenko, P. V. Mamyshev, A. N. Pilipetskii, and E. M. Dianov, "Numerical analysis of the Raman spectrum evolution and soliton pulse generation in single-mode fibers," J. Opt. Soc. Am. B 8, 1626-1632 (1991). http://www.opticsinfobase.org/abstract.cfm?URI=josab-8-8-1626. [CrossRef]
M. N. Islam, G. Sucha, I. Bar-Joseph, M. Wegener, J. P. Gordon, and D. S. Chemla, "Femtosecond distributed soliton spectrum in fibers," J. Opt. Soc. Am. B 6, 1149-1158 (1989). http://www.opticsinfobase.org/abstract.cfm?URI=josab-6-6-1149. [CrossRef]
R. G. Smith, "Optical power handling capacity of low loss optical fibers as determined by stimulated Raman and Brillouin scattering," Appl. Opt. 11, 2489-2494 (1972). [CrossRef] [PubMed]
S. B. Cavalcanti, G. P. Agrawal, and M. Yu, "Noise amplification in dispersive nonlinear media," Phys. Rev. A 51, 4086-4092 (1995). http://dx.doi.org/10.1103/PhysRevA.51.4086. [CrossRef] [PubMed]
G. P. Agrawal, Nonlinear Fiber Optics, 3rd ed. (Academic Press, San Diego, CA, USA, 2001).
P. L. François, "Nonlinear propagation of ultrashort pulses in optical fibers: total field formulation in the frequency domain," J. Opt. Soc. Am. B 8, 276-293 (1991). http://www.opticsinfobase.org/abstract.cfm?URI=josab-8-2-276. [CrossRef]
S. G. Johnson and J. D. Joannopoulos, "Block-Iterative Frequency-Domain Methods for Maxwell’s Equations in a Planewave Basis," Opt. Express 8, 173-190 (2001). [CrossRef] [PubMed]
K. P. Hansen, Crystal Fibre A/S, Personal communication.
J. Lægsgaard, N. A. Mortensen, and A. Bjarklev, "Mode areas and field-energy distribution in honeycomb photonic bandgap fibers," J. Opt. Soc. Am. B 20, 2037-2045 (2003). [CrossRef]
K. J. Blow and D. Wood, "Theoretical description of transient stimulated Raman scattering in optical fibers," IEEE J. Quantum Electron. 25, 2665-2673 (1989). [CrossRef]
J. W. Goodman, Statistical Optics (John Wiley & Sons Inc., 2000). ISBN 0471399167.
W. K. Marshall, B. Crosignani, and A. Yariv, "Laser phase noise to intensity noise conversion by lowest-order group-velocity dispersion in optical fiber:exact theory," Opt. Lett. 25, 165-167 (2000). http://www.opticsinfobase.org/abstract.cfm?URI=ol-25-3-165. [CrossRef]
W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical Recipes in C++: The Art of Scientific Computing, 2nd ed. (Cambridge University Press, Cambridge, 2002). http://www.nr.com.
V. E. Zakharov and A. B. Shabat, "Exact theory of 2-dimensional self-focusing and one-dimenstional selfmodulation of waves in nonlinear media," Sov. Phys. JETP 34, 62 (1972).
H. A. Haus and W. S. Wong, "Solitons in optical communications," Rev. Mod. Phys. 68, 423-444 (1996). http://link.aps.org/abstract/RMP/v68/p423. [CrossRef]
A. Peleg, M. Chertkov, and I. Gabitov, "Interchannel interaction of optical solitons," Phys. Rev. E 68, 026605 (2003). http://link.aps.org/abstract/PRE/v68/e026605. [CrossRef]
S. Chi and S. Wen, "Raman cross talk of soliton collision in a lossless fiber," Opt. Lett. 14, 1216-1218 (1989). [CrossRef] [PubMed]
B. A. Malomed, "Soliton-collision problem in the nonlinear Schr¨odinger equation with a nonlinear damping term," Phys. Rev. A 44, 1412-1414 (1991). http://link.aps.org/abstract/PRA/v44/p1412. [CrossRef] [PubMed]
Y. Chung and A. Peleg, "Strongly non-Gaussian statistics of optical soliton parameters due to collisions in the presence of delayed Raman response," Nonlinearity 18, 1555-1574 (2005). http://stacks.iop.org/Non/18/1555. [CrossRef]
O. Bang and M. Peyrard, "Generation of high-energy localized vibrational modes in nonlinear Klein-Gordon lattices," Phys. Rev. E 53, 4143-4152 (1996). http://link.aps.org/abstract/PRE/v53/p4143. [CrossRef]
J. K. Lucek and K. J. Blow, "Soliton self-frequency shift in telecommunications fiber," Phys. Rev. A 45, 6666-6674 (1992). http://dx.doi.org/10.1103/PhysRevA.45.6666. [CrossRef] [PubMed]
Y. Kodama and K. Nozaki, "Soliton interaction in optical fibers," Opt. Lett. 12, 1038-1040 (1987). [CrossRef] [PubMed]
J. P. Gordon, "Theory of the soliton self-frequency shift," Opt. Lett. 11, 662-664 (1986). [CrossRef] [PubMed]
J. Herrmann and A. Nazarkin, "Soliton self-frequency shift for pulses with a duration less than the period of molecular oscillations," Opt. Lett. 19, 2065-2067 (1994). [CrossRef] [PubMed]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Click here to see a list of articles that cite this paper