OSA's Digital Library

Optics Express

Optics Express

  • Editor: Michael Duncan
  • Vol. 11, Iss. 5 — Mar. 10, 2003
  • pp: 452–459
« Show journal navigation

Spatial soliton formation in photonic crystal fibers

Albert Ferrando, Mario Zacarés, Pedro Fernández de Córdoba, Daniele Binosi, and Juan A. Monsoriu  »View Author Affiliations


Optics Express, Vol. 11, Issue 5, pp. 452-459 (2003)
http://dx.doi.org/10.1364/OE.11.000452


View Full Text Article

Acrobat PDF (506 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

We demonstrate the existence of spatial soliton solutions in photonic crystal fibers (PCF’s). These guided localized nonlinear waves appear as a result of the balance between the linear and nonlinear diffraction properties of the inhomogeneous photonic crystal cladding. The spatial soliton is realized self-consistently as the fundamental mode of the effective fiber defined simultaneously by the PCF linear and the self-induced nonlinear refractive indices. It is also shown that the photonic crystal cladding is able to stabilize these solutions, which would be unstable otherwise if the medium was entirely homogeneous.

© 2003 Optical Society of America

1. Introduction

The existence of self-trapped light beam solutions, or nonlinear localized guided waves, in homogeneous media is a classical result in nonlinear optics [1

1. R. Y. Chiao, E. Garmire, and C. H. Townes, “Self-trapping of optical beams,” Phys. Rev. Lett. 13, 479–482 (1964). [CrossRef]

, 2

2. H. A. Haus, “Higher order trapped light beam solutions,” Appl. Phys. Lett. 8, 128–129 (1966). [CrossRef]

]. These rotationally symmetric solutions are, however, unstable and only exist at a critical power. More recently, it has been proven that stable modes can be localized in two-dimensional (2D) nonlinear photonic crystals constituted by periodic arrays of poles of nonlinear material both for infinite-sized [3

3. S. John and N. Aközbek, “Nonlinear optical solitary waves in a photonic band gap,” Phys. Rev. Lett. 71, 1168–1171 (1993). [CrossRef] [PubMed]

, 4

4. S. F. Mingaleev and Y. S. Kivshar, “Self-trapping and stable localized modes in nonlinear photonic crystals,” Phys. Rev. Lett. 86, 5474–5477 (2001). [CrossRef] [PubMed]

] and finite-sized structures [5

5. P. Xie, Z.-Q. Zhang, and X. Zhang, “Gap solitons and soliton trains in finite-sized two-dimensional periodic and quasiperiodic photonic crystals,” Phys. Rev. E 67, 026607-1–026607-5 (2003). [CrossRef]

]. These solutions are non propagating in the direction of invariance of the photonic crystal (axial direction) and, physically, they correspond to the situation of in-plane illumination. Guided localized nonlinear solutions have also been found in 2D linear photonic crystals of dielectric poles with embedded nonlinear arrays [6

6. S. F. Mingaleev, Y. S. Kivshar, and R. A. Sammut, “Long-range interaction and nonlinear localized modes in photonic crystal waveguides,” Phys. Rev. E 62, 5777–5782 (2000). [CrossRef]

] (also for in-plane illumination) and, experimentally, in 1D waveguide arrays, in the form of the so-called discrete solitons propagating along the waveguide direction [7

7. H. S. Eisenberg, Y. Silberberg, R. Morandotti, A. R. Boyd, and J. S. Aitchison, “Discrete spatial optical solitons in waveguide arrays,” Phys. Rev. Lett. 81, 3383–3387 (1998). [CrossRef]

]. Modeling has been performed using an effective “quasi-free photon approximation”[3

3. S. John and N. Aközbek, “Nonlinear optical solitary waves in a photonic band gap,” Phys. Rev. Lett. 71, 1168–1171 (1993). [CrossRef] [PubMed]

] (valid for lowindex contrasts) or a discretized version of the scalar nonlinear equation [4

4. S. F. Mingaleev and Y. S. Kivshar, “Self-trapping and stable localized modes in nonlinear photonic crystals,” Phys. Rev. Lett. 86, 5474–5477 (2001). [CrossRef] [PubMed]

].

2. Description of the method

We consider the guidance problem, at a given frequency, in a triangular PCF which can have a nonlinear response in silica. That is, we solve the following nonlinear equation:

[t2+k02(n02(x)+n22(x)ϕ2)]ϕ=β2ϕ,
(1)

Fig. 1. (a) Schematic representation of the transverse section of a PCF. (b) The box with dimensions D×D corresponds to the unit cell used to implement periodic boundary conditions. In our simulations we have chosen D=7Λ, Λ being the spatial periodicity, or pitch, of the photonic crystal cladding.

3. Results

We take normalized field amplitudes (∫ ϕ*ϕ/A 0=1), which is equivalent to define a power-dependent dimensionless nonlinear coupling γ=Pn2(silica)2/A 0, where P is the total power carried by the unnormalized field and A 0 is a magnitude with dimensions of area characterizing the core size (here we choose A 0=π(Λ/2)2, Λ being the spatial periodicity—or pitch—of the PCF). Solutions of the nonlinear eigenvalue equation are found for different values of the PCF parameters and the nonlinear coupling. In all cases, the solution we find correspond to the eigenvector with highest β2 of the nonlinear eigenvalue problem. The solution can be envisaged as the fundamental mode of the effective fiber generated by the combined effect of the PCF refractive index and the nonlinear index induced by the solution amplitude itself (n 2(x⃗)=n02(x⃗)+n22(x⃗)|ϕsol(x⃗)|2). Since the fundamental mode of the linear theory is also a localized solution, it is convenient to provide an efficient criterion to establish a distinction between localized nonlinear solutions and the fundamental mode of a linear PCF. We use as a measure of distinction the difference, or gap, (normalized to the given vacuum wavenumber k 0) Δ≡(βsolfund)/k 0, where βsol and βfund are the propagation constants of the solution of the nonlinear equation and the fundamental mode of the linear fiber, respectively. The generation of a gap implies that the nonlinear solution has a different shape than the fundamental linear mode.

Although Eq. (1) admits solutions for ordinary-scale PCF’s, their experimental feasibility is uncertain because of the high intensities involved when focusing light in such small cores (a typical core radius of a conventional PCF is approximately 1 µm). These high intensities can produce a silica breakdown when they overcome a critical value. In terms of our dimensionless nonlinear coupling, this critical intensity implies the existence of a maximum value of γ above which the fiber collapses (γmax≈10-3). For these reasons, it is more realistic to consider PCF structures that possess larger cores. One simple way of generating such structures is by magnifying the conventional PCF parameters, i.e., the hole radius a and the pitch Λ, (aMa, ΛMΛ). The core size area A 0 is then scaled by a M 2 factor and γ by 1/M 2. In our simulation, we choose M=10 with respect to an ordinary PCF configuration characterized by Λ=2.3 µm and a ranging from 0.2 to 1.0 microns. That is, we simulate large-scale PCF’s with Λ=23 µm and a=2–10 µm. In this way, we guarantee that all solutions we find will lie below the breakdown intensity threshold.

Using these values, our simulation shows the existence of nonlinear guided solutions in all PCF configurations even for small nonlinear couplings. The shape of the soliton solutions change with γ in the form presented in Fig. 2(a)-(c). As the nonlinear coupling increases, the self-consistent solution is gradually narrower. In terms of the gap function Δ, this fact is reflected in a monotonously increasing behavior of the generated gap with γ, as shown in Fig. 3(a) (solid curves). However, the nonlinear coupling cannot be increased arbitrarily. There is a given value of γ (not to be confused with the maximum γ for silica breakdown), which can depend on a, at which the soliton acquires a critical shape and propagation constant (Fig. 2(d)). At the critical coupling, our calculations recover the rotationally invariant nonlinear solution obtained in an homogeneous medium (also known as Townes soliton) [1

1. R. Y. Chiao, E. Garmire, and C. H. Townes, “Self-trapping of optical beams,” Phys. Rev. Lett. 13, 479–482 (1964). [CrossRef]

]. Above this critical value, our method finds infinitely narrowso lutions (within numerical precision) characterized by an infinite eigenvalue β2 (Fig. 2(e)). The physical interpretation of these solutions is clear. As we increase the nonlinear coupling, the localized mode is more confined and it feels the photonic crystal cladding less and less. At the critical coupling, it stops seeing the cladding completely. Thus, its behavior corresponds to a nonlinear mode in an homogeneous medium. That is why we recover the rotationally invariant solution (no discrete symmetry is left) of an homogeneous medium. This argument also applies above the critical coupling, where we expect to recover the same physics as in an homogeneous medium. It is known that the Townes soliton in an homogeneous medium is unstable under power perturbations. In particular, if power is increased the solution collapses experimenting a process of filamentation induced by a self-focusing instability and generating an infinitely narrow solution with an infinite β2 [12

12. G. Fibich and A. L. Gaeta, “Critical power for self-focusing in bulk media and in hollow waveguides,” Opt. Lett. 25335–337 (2000). [CrossRef]

]. This is the same solution our method detects above the critical coupling (Fig. 2(e)).

It is interesting to notice that we have also simulated with our method the limit where all holes disappear (a→0), which corresponds to the homogeneous medium case. The differences and similarities with respect to the PCF case are enlightening. Below the critical coupling, the method does not detect any nontrivial solution and only finds the trivial zero solution of Eq. (1). This is also consistent with known results for an homogeneous medium, since no stationary nonlinear solutions are found below the critical power that corresponds to the Townes soliton [1

1. R. Y. Chiao, E. Garmire, and C. H. Townes, “Self-trapping of optical beams,” Phys. Rev. Lett. 13, 479–482 (1964). [CrossRef]

, 12

12. G. Fibich and A. L. Gaeta, “Critical power for self-focusing in bulk media and in hollow waveguides,” Opt. Lett. 25335–337 (2000). [CrossRef]

]. At and above the critical coupling, the Townes soliton profile and the infinitely narrow solution are recovered, respectively.

Fig. 2. Intensity distribution of different solutions of Eq. (1) for a PCF with pitch Λ=23µm, radius a=4µm, λ=1.55 µm and different nonlinear couplings: (a), linear mode (γ=0); (b) and (c), spatial soliton solutions for γ=0.0010 and γ=0.0015, respectively;(d), unstable nonlinear solution of the homogeneous medium (γ=γc=0.0017); and (e), self-focusing instability (γ>γc).
Fig. 3. (a) Dependence of the gap function Δ on the nonlinear coupling γ for different hole sizes. As in Fig. 2, Λ=23µm and λ=1.55µm. (b) Diagram of existence of solutions for a nonlinear PCF. The shaded region is the nonlinear soliton phase. The other region corresponds to the homogeneous-instability phase. The inter-phase is given by the γc(a) curve. As before, Λ=23 µm and λ=1.55 µm.

The soliton formation pattern depends on the interplay between the PCF geometry and the nonlinear coupling. It can be more easily visualized using the diagram of existence of solutions in Fig. 3(b). A systematic study of different soliton solutions of Eq. (1) for different values of the hole-radius a and the nonlinear coupling γ allows to establish two different regions, or phases, separated by their corresponding inter-phase curve. Since nonlinearities can dynamically induce a propagation constant gap with respect to the linear regime, the phase diagram in Fig. 3(b) is constructed according to the value of the gap function Δ. A finite non-zero gap characterizes the region where spatial soliton solutions with finite propagation constant exist. The region above the critical coupling corresponds to infinitely narrowso lutions characterized by an infinite gap. As mentioned previously, they are the result of a self-focusing instability acting on the Townes soliton of the homogeneous medium when the critical power is overcome. This fact justifies to call the white area in Fig. 3(b) as the self-focusing instability region. The curve that describes the inter-phase between the two regions is given by the function γc(a). Solutions lying right on the interface curve γc(a) constitute the family of Townes solitons of the homogeneous medium characterized by different beam diameters. PCF’s with larger holes radii a generate more localized solutions and, consequently, also give rise to Townes solitons with smaller beam diameters. Their propagation constant is thus larger and so the gap function is, as reflected in Fig. 3(a). Notice that our numerical simulation provides a nearly vertical inter-phase curve γc(a). This is consistent with known results for an homogeneous medium, where all solutions of the Townes soliton family are characterized by the same critical power, independently of the beam diameter [12

12. G. Fibich and A. L. Gaeta, “Critical power for self-focusing in bulk media and in hollow waveguides,” Opt. Lett. 25335–337 (2000). [CrossRef]

, 13

13. R. W. Boyd, Nonlinear Optics (Academic Press, 1992).

]. In terms of our dimensionless nonlinear coupling constant, γ=Pn2(silica)2/A 0, solutions lying right on the interface curve must share an identical critical γc(a). It is also interesting to point out that in the homogeneous case the diagram of existence would be reduced to the vertical inter-phase curve γc(a) exclusively.

Fig. 4. (2.55 MB) Evolution of the field amplitude in z for a large-scale PCF with Λ=23µm, a=8µm and λ=1.55µm. We show the transient from an initial Gaussian profile towards an asymptotic spatial soliton solution.

Up to now, we have been only concerned with the existence of nonlinear solutions. We have found solutions of the nonlinear eigenvalue equation (1) but their stability has not been proven yet. We can only claim that nonlinear solutions lying on the inter-phase curve γc(a) are unstable since we recognized them as Townes solitons of an homogeneous medium and instability is a known feature of them. Stability is a dynamical issue which implies the study of the evolution equation associated to Eq. (1):

[t2+k02(n02(x)+n22(x)ϕ2)]ϕ=2ϕz2.
(2)

Fig. 5. Typical evolution behavior in z of 〈ϕ|L(ϕ)|ϕ〉 for a large-scale PCF with Λ=23µm, a=8µm and λ=1.55µm.

In such a case, 〈ϕ|L|ϕ〉=β2. In Fig. 5 we show a typical behavior of the evolution of 〈ϕ|L|ϕ〉 (z) where the stabilization mechanism becomes evident: the expectation value evolves (oscillating) until it asymptotically reaches a plateau indicating the existence of a spatial soliton solution. The asymptotic β2 value corresponds to the spatial soliton solution of Eq. (1) for the same PCF parameters with a γ given by the asymptotic power. Notice that power is not conserved because of radiation and decreases during propagation, so that the asymptotic solution possesses a smaller power than any other intermediate solution. It is worth mentioning that we have also checked that these spatial solitons are stable under both small transverse displacements relative to the photonic crystal cladding and launching with small transverse momentum (slight offaxis illumination). Analysis of the evolution of the expectation value 〈ϕ|L|ϕ〉 (z) shows in both cases a similar asymptotic behavior as for the Gaussian profile in Fig. 5: damped oscillations tending towards a spatial soliton plateau.

4. Conclusions

We have demonstrated that large-scale PCF’s can support stable nonlinear localized solutions. They constitute a newclass of spatial soliton solutions characterized by a discrete symmetry and different from the so-called discrete spatial solitons. These new spatial solitons in PCF’s distinguish from discrete spatial solitons in the fact that they cannot be described by the discretized nonlinear Schrödinger equation. Despite both enjoy the sixth-fold symmetry of the photonic crystal, their physical properties are substantially different. We have shown that spatial solitons in PCF’s are closer to nonlinear solutions in a continuous medium. This feature is apparent in the fact that the unstable rotationally invariant solution of an homogeneous medium is recovered for specific large nonlinear couplings and also in the appearance of identical self-focusing instabilities.

On the other hand, our results point out the crucial role played by the photonic crystal cladding in the stabilization mechanism of these newso lutions. Besides the ordinary interplay between diffraction and nonlinearity, the photonic crystal cladding provides an extra and peculiar localization mechanism. In this way, even for small effective nonlinear couplings (or equivalently, for small input powers) stable nonlinear solutions can be generated. Without the existence of the additional confinement provided by the photonic crystal cladding, nonlinearities would not be able to build such solutions in this regime, as known in the homogeneous case. Remarkably enough, the photonic crystal structure does not only help to generate these solutions but it also gives rise to their stabilization.

Finally, we present a proposal that overcomes at first instance the problems derived from collapsing the fiber due to high intensities. Large-scale PCF’s prevent intensities from reaching values above the threshold silica breakdown.

We are thankful to Alex Gaeta for useful discussions. This work was financially supported by the Plan Nacional I+D+I (grant TIC2002-04527-C02-02), Ministerio de Ciencia y Tecnología (Spain) and FEDER funds. M. Zacarés gratefully acknowledges Fundación Ramón Areces grant. D. Binosi gratefully acknowledges a grant, BFM2001-0262, from Ministerio de Educación y Cultura, Spain.

References and links

1.

R. Y. Chiao, E. Garmire, and C. H. Townes, “Self-trapping of optical beams,” Phys. Rev. Lett. 13, 479–482 (1964). [CrossRef]

2.

H. A. Haus, “Higher order trapped light beam solutions,” Appl. Phys. Lett. 8, 128–129 (1966). [CrossRef]

3.

S. John and N. Aközbek, “Nonlinear optical solitary waves in a photonic band gap,” Phys. Rev. Lett. 71, 1168–1171 (1993). [CrossRef] [PubMed]

4.

S. F. Mingaleev and Y. S. Kivshar, “Self-trapping and stable localized modes in nonlinear photonic crystals,” Phys. Rev. Lett. 86, 5474–5477 (2001). [CrossRef] [PubMed]

5.

P. Xie, Z.-Q. Zhang, and X. Zhang, “Gap solitons and soliton trains in finite-sized two-dimensional periodic and quasiperiodic photonic crystals,” Phys. Rev. E 67, 026607-1–026607-5 (2003). [CrossRef]

6.

S. F. Mingaleev, Y. S. Kivshar, and R. A. Sammut, “Long-range interaction and nonlinear localized modes in photonic crystal waveguides,” Phys. Rev. E 62, 5777–5782 (2000). [CrossRef]

7.

H. S. Eisenberg, Y. Silberberg, R. Morandotti, A. R. Boyd, and J. S. Aitchison, “Discrete spatial optical solitons in waveguide arrays,” Phys. Rev. Lett. 81, 3383–3387 (1998). [CrossRef]

8.

A. Ferrando, E. Silvestre, J. J. Miret, P. Andrés, and M. V. Andrés, “Donor and acceptor guided modes in photonic crystal fibers,” Opt. Lett. 25, 1238–1330 (2000). [CrossRef]

9.

A. Ferrando, E. Silvestre, J. J. Miret, P. Andrés, and M. V. Andrés, “Full-vector analysis of a realistic photonic crystal fiber,” Opt. Lett. 24, 276–278 (1999). [CrossRef]

10.

A. Ferrando, E. Silvestre, J. J. Miret, P. Andrés, and M. V. Andrés, “Vector description of higherorder modes in photonic crystal fibers,” J. Opt. Soc. Am. A 17, 1333–1340 (2000). [CrossRef]

11.

A. Ferrando, M. Zacarés, and P. F. de Córdoba, “Ansatz-independent solution of a soliton in a strong dispersion-management system,” Phys. Rev. E 62, 7320–7329 (2000). [CrossRef]

12.

G. Fibich and A. L. Gaeta, “Critical power for self-focusing in bulk media and in hollow waveguides,” Opt. Lett. 25335–337 (2000). [CrossRef]

13.

R. W. Boyd, Nonlinear Optics (Academic Press, 1992).

OCIS Codes
(060.4370) Fiber optics and optical communications : Nonlinear optics, fibers
(190.4370) Nonlinear optics : Nonlinear optics, fibers

ToC Category:
Research Papers

History
Original Manuscript: January 30, 2003
Revised Manuscript: February 26, 2003
Published: March 10, 2003

Citation
Albert Ferrando, Mario Zacares, Pedro Fernandez de Cordoba, Daniele Binosi, and Juan Monsoriu, "Spatial soliton formation in photonic crystal fibers," Opt. Express 11, 452-459 (2003)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-11-5-452


Sort:  Journal  |  Reset  

References

  1. R. Y. Chiao, E. Garmire, and C. H. Townes, �??Self-trapping of optical beams,�?? Phys. Rev. Lett. 13, 479�??482 (1964). [CrossRef]
  2. H. A. Haus, �??Higher order trapped light beam solutions,�?? Appl. Phys. Lett. 8, 128�??129 (1966). [CrossRef]
  3. S. John and N. Akozbek, �??Nonlinear optical solitary waves in a photonic band gap,�?? Phys. Rev. Lett. 71, 1168�??1171 (1993). [CrossRef] [PubMed]
  4. S. F. Mingaleev and Y. S. Kivshar, �??Self-trapping and stable localized modes in nonlinear photonic crystals,�?? Phys. Rev. Lett. 86, 5474�??5477 (2001). [CrossRef] [PubMed]
  5. P. Xie, Z.-Q. Zhang, and X. Zhang, �??Gap solitons and soliton trains in .nite-sized two-dimensional periodic and quasiperiodic photonic crystals,�?? Phys. Rev. E 67, 026607-1�?? 026607-5 (2003). [CrossRef]
  6. S. F.Mingaleev, Y. S. Kivshar, and R. A. Sammut, �??Long-range interaction and nonlinear localized modes in photonic crystal waveguides,�?? Phys. Rev. E 62, 5777�??5782 (2000). [CrossRef]
  7. H. S. Eisenberg, Y. Silberberg, R. Morandotti, A. R. Boyd, and J. S. Aitchison, �??Discrete spatial optical solitons in waveguide arrays,�?? Phys. Rev. Lett. 81, 3383�??3387 (1998). [CrossRef]
  8. A. Ferrando, E. Silvestre, J. J. Miret, P. Andres, and M. V. Andres, �??Donor and acceptor guided modes in photonic crystal fibers,�?? Opt. Lett. 25, 1238�??1330 (2000). [CrossRef]
  9. A. Ferrando, E. Silvestre, J. J. Miret, P. Andres, and M. V. Andres, �??Full-vector analysis of a realistic photonic crystal fiber,�?? Opt. Lett. 24, 276�??278 (1999). [CrossRef]
  10. A. Ferrando, E. Silvestre, J. J. Miret, P. Andres, and M. V. Andres, �??Vector description of higherorder modes in photonic crystal fibers,�?? J. Opt. Soc. Am. A 17, 1333�??1340 (2000). [CrossRef]
  11. A. Ferrando, M. Zacares, and P. F. de Cordoba, �??Ansatz-independent solution of a soliton in a strong dispersion-management system,�?? Phys. Rev. E 62, 7320�??7329 (2000). [CrossRef]
  12. G. Fibich and A. L. Gaeta, �??Critical power for self-focusing in bulk media and in hollow waveguides,�?? Opt. Lett. 25 335�??337 (2000). [CrossRef]

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.

Multimedia

Multimedia FilesRecommended Software
» Media 1: GIF (2611 KB)     

« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited