## Designing the properties of dispersion-flattened photonic crystal fibers

Optics Express, Vol. 9, Issue 13, pp. 687-697 (2001)

http://dx.doi.org/10.1364/OE.9.000687

Acrobat PDF (738 KB)

### Abstract

We present a systematic study of group-velocity-dispersion properties in photonic crystal fibers (PCF’s). This analysis includes a thorough description of the dependence of the fiber geometrical dispersion on the structural parameters of a PCF. The interplay between material dispersion and geometrical dispersion allows us to established a well-defined procedure to design specific predetermined dispersion profiles. We focus on flattened, or even ultraflattened, dispersion behaviors both in the telecommunication window (around 1.55 *µ*m) and in the *Ti-Za* laser wavelength range (around 0.8 *µ*m). We show the different possibilities of obtaining normal, anomalous, and zero dispersion curves in the above frequency domains and discuss the limits for the existence of the above dispersion profiles.

© Optical Society of America

## 1 Introduction

*a*and the lattice period, or pitch, Λ of a 2D triangular photonic crystal cladding) to generate an enormous variety of different configurations.

1. 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]

*µ*m where found [2

2. D. Mogilevtsev, T. A. Birks, and P. S. J. Russell, “Dispersion of photonic crystal fibers,” Opt. Lett. **23**, 1662–1664 (1998). [CrossRef]

3. M. J. Gander, R. McBride, J. D. C. Jones, D. Mogilevtsev, T. A. Birks, J. C. Knight, and P. S. J. Russell, “Experimental measurement of group velocity dispersion in photonic crystal fibers,” Electron. Lett. **35**, 63–64 (1999). [CrossRef]

4. P. J. Bennet, T. M. Monro, and D. J. Richardson, “Toward practical holey fiber technology: fabrication, splicing, modeling, and fabrication,” Opt. Lett. **24**, 1203–1205 (1999). [CrossRef]

5. A. Ferrando, E. Silvestre, J. J. Miret, J. A. Monsoriu, M. V. Andrés, and P. S. J. Russell, “Designing a photonic crystal fibre with flattened chromatic dispersion,” Electron. Lett. **24**, 325–327 (1999). [CrossRef]

6. J. Broeng, D. Mogilevtsev, S. E. Barkou, and A. Bjarklev, “Photonic crys tal fibers : a new clas s of optical waveguides,” Opt. Fib. Tech. **5**, 305–330 (1999). [CrossRef]

7. A. Ferrando, E. Silvestre, J. J. Miret, and P. Andrés, “Nearly zero ultraflattened dispersion in photonic crystal fibers,” Opt. Lett. **25**, 790–792 (2000). [CrossRef]

7. A. Ferrando, E. Silvestre, J. J. Miret, and P. Andrés, “Nearly zero ultraflattened dispersion in photonic crystal fibers,” Opt. Lett. **25**, 790–792 (2000). [CrossRef]

*D*is proportional to the second derivative of the modal effective index with respect to the wavelength λ. For this reason, it is very sensitive to the precision with which the dependence of

*n*

_{eff}with λ is calculated. On the other hand, the calculation method has to provide a reasonable computational time for evaluating the dispersion curves of different configurations. Ideally, it should also be flexible to accommodate divers geometric proposals (such as different photonic-crystal lattice geometries, non-perfect lattices and holes, asymmetries, different materials, and so on) with a little effort. Here, we use a full-vector modal method developed by our own group to describe electromagnetic propagation in general systems with translational invariance [8

8. E. Silvestre, M. V. Andrés, and P. Andrés, “Biorthonormal-basis method for the vector description of optical-fiber modes,” J. Lightwave Technol. **16**, 923–928 (1998). [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 higher-order modes in photonic crystal fibers,” J. Opt. Soc. Am. A **17**, 1333–1340 (2000). [CrossRef]

## 2 Designing Procedure

*a*and the lattice period, or pitch, Λ. However, in the design procedure we have recognized that it is more convenient to consider two alternate dimensionless parameters instead. First, we consider the so-called filling fraction

*f*, defined as

*f*=(4π/3)(

*a*/Λ)

^{2}, that involves the dimensionless ratio

*a*/Λ and provides the proportion of air with respect to silica in the photonic crystal structure. A change in

*f*produces a variation of the amount of air in the structure, as shown in Fig. 1(a). The second parameter we take into account is the magnification

*M*, which simply consists in a simultaneous scale transformation of both

*a*and Λ in the same amount, as shown in Fig. 1(b). In order to define

*M*operatively, it is convenient to select a reference value of the pitch (in our case, w e choose Λ=2.3

*µ*m). The magnification

*M*has also an appealing practical interest. In the pulling process during the fabrication of the fiber,

*M*is the parameter that can be controlled in a natural way. This is so because, under optimal conditions, the pulling process should preserve the proportions of the original structure.

*n*

_{eff}=

*β*[λ,

*n*

_{m}(λ)]/

*k*

_{0},

*β*is the propagation constant,

*k*

_{0}=2

*π*/λ is the free-space w ave number, and

*n*

_{m}=

*n*

_{m}(λ) is the chromatic dispersion of the material, silica in this case. According to the above equation, there are two different sources of obtaining non-zero dispersion due to the existence of two different types of dependence of

*β*on λ. One of them is originated by the explicit dependence of the propagation constant of the mode on λ and it occurs even if the material is, or it can be considered, non-dispersive (

*n*

_{m}(λ)=const). Since the dispersion generated in this way is not produced by the chromatic dispersion of the material but by the geometry of the PCF refractive index distribution that determines the dispersion relation of the guided mode,

*β*=

*β*[λ,

*n*

_{m}(λ)=const], we call it geometrical dispersion. Its definition is, accordingly, the same as in Eq. (1) but supplemented with the condition that the material is non-dispersive;

*β*on λ through the chromatic dispersion of the material,

*n*

_{m}=

*n*

_{m}(λ). Consequently, this type of dispersion is called material dispersion,

*D*

_{m}, and we calculate it as in Eq. (1) by substituting

*n*

_{eff}(λ) by

*n*

_{m}(λ).

*D*by a sum of the geometrical and material dispersion [11];

*D*is written in this way. The virtue of Eq. (2) is that permits to split both sources of dispersion into two different terms explicitly.

*n*

_{m}(λ) is an input of the problem and consequently, so is

*D*

_{m}. All the design power is stored in the geometrical dispersion, In this sense it is very important to recognize the following fact. The effective refractive index of a guided mode

*n*

_{eff}, for the calculation of which we assume no material dispersion,

*explicitly*depends on the photonic crystal cladding parameters,

*a*and Λ, and the wavelength λ. Inasmuch as

*n*

_{eff}is a dimensionless function, this dependence can only occur through dimensionless ratios of these three parameters. For our discussion, it is convenient to take as independent parameters

*a*/λ and Λ/λ, so that

*n*

_{eff}=

*n*

_{eff}(

*a*/λ,Λ/λ). This property determines the dependence of

*D*

_{g}on

*M*completely. According to the definition of the geometrical dispersion, it is clear that under a scale transformation of λ, we obtain

*f*and setting

*M*=1, or equivalently, fixing

*a*and Λ=2.3

*µ*m) to analytically obtain all the family of dispersion curves parametrized by

*M*, as shown in Fig. 2(a). The linear part of these curves modifies its slope and it is simulteanously shitfed when

*M*is changed.

*D*

_{g}with the filling fraction

*f*. In practice, the only way to determine this dependence is by calculating

*D*

_{g}numerically. We thus start with a reference configuration, e.g.

*M*=1 (i.e., with Λ=2.3

*µ*m) and evaluate the geometrical dispersion curves for different filling fractions f simply by changing

*a*. The result is represented in Fig. 2(b). The remarkable feature of these curves is that, besides they are shifted, the slope of their linear part is approximately preserved when the filling fraction is changed. This property will show to be very helpful in the design process.

*D*

_{g}, the sign-changed material dispersion -

*D*

_{m}, and the total dispersion

*D*, are represented in blue, black and red, respectively. According to Eq. (4), the red curve corresponding to total dispersion is obtained by subtracting the values of the black curve from the blue one.

7. A. Ferrando, E. Silvestre, J. J. Miret, and P. Andrés, “Nearly zero ultraflattened dispersion in photonic crystal fibers,” Opt. Lett. **25**, 790–792 (2000). [CrossRef]

*D*

_{g}. The sign-changed material dispersion -

*D*

_{m}is a smooth curve in most of the infrared region, so that it can be well approximated by a linear function around different λ’s belonging to this region over pretty wide intervals. It is clear, in view of Fig. 3, that in the region of λ’s in which the linear part of both the black and blue curves can be set parallel, the total dispersion will achieve an ideally perfect flattened behavior.

*D*

_{m}curve is fixed, we perform a scale transformation of

*D*

_{g}parametrized by the magnification

*M*in such a way that provides an scaled

*D*

_{g}curve having a linear region with the same given slope. If the wavelength region (centered at the specified λ) where

*D*

_{m}behaves linearly overlaps the wavelength region of linear behavior of

*D*

_{g}, we will obtain an ultraflattened total dispersion curve in the overlapping wavelength range.

*M*. However, it remains still one degree of freedom to play with, the filling fraction

*f*. As shown in Fig. 2(b), note that a change in

*f*does not alter the value of the slope of the linear part of

*D*

_{g}. Therefore, if we proceed to change the value of

*f*preserving the value of

*M*obtained above, the difference between both curves will change, and also the overlapping range, but the parallelism condition will remain unaltered. This means that one simultaneously modifies the value of the total dispersion

*D*and the width of the wavelength window where ultraflattened behavior occurs by acting on the filling fraction

*f*. Since these two properties are not independent, one has the choice to select

*f*either to obtain a desired value of

*D*or to maximize the range of ultraflattened dispersion operation. In both cases, all possible configurations will provide ultraflattened dispersion profiles.

*M*and

*f*providing a previously established ultraflattened dispersion behavior at a reference wavelength. The design procedure does not supply us with the exact parameters because it is based on an approximate expression for

*D*(see Eq. (2)). If one is interested in a more accurate evaluation of

*M*and

*f*, one can undertake a second fine-tuning search of the actual parameters starting from their approximate values. This search is performed by calculating the real dispersion exactly, that is, by including the chromatic dispersion of the material explicitly in the determination of the dispersion relation of the mode (see Eq. (1)). Since the starting values are pretty close to the real ones, the fine-tuning process converges rapidly.

## 3 Some Specific Designs

*µ*m. Using the procedure described in the previous section, it is possible to systematically obtain ultraflattened dispersion configurations using 1.55

*µ*m as the reference wavelength providing the magnification

*M*. These configurations can be designed to show positive, negative or nearly-zero

*D*just by properly adjusting the value of the filling fraction

*f*. However, despite the dispersion profiles obtained this way are certainly ultraflattened, they are not necessarily optimized. This means that they will not provide the widest windows of ultraflattened behavior. Nevertheless, they constitute a good starting point to search configurations with optimal width at a given

*D*.

*D*. We would like to point out that it is easier to achieve broader wavelength windows for ultraflattened curves with positive dispersion. This feature can be understood by analyzing the way they are obtained in the design process. If we look at the material dispersion in Fig. 3, we see that it is precisely in this region where -

*D*

_{m}behaves more smoothly, thus providing a larger wavelength interval of linear behavior. On the other hand, we need

*D*

_{g}curves with the same slope as the -

*D*

_{m}curve at the reference wavelength but, and this is the crucial point, that simultaneously remain above it. This means we need to use curves corresponding to higher filling fraction in Fig. 2(b) to achieve increasingly higher positive values of

*D*. This is so because the larger

*f*is, the more shifted to higher λ’s these curves are, and, consequently, the higher the total dispersion becomes. However, this is not the only effect that occurs as

*f*increases. It turns out that the linear part of these curves also increases when shifted to higher λ’s, as depicted in Fig. 2(b). As a consequence, the overlapping region of linear behavior of the -

*D*

_{m}and

*D*

_{g}curves potentially becomes larger as

*D*increases. The final outcome is a wider ultraflattened dispersion profiles for higher positive

*D*’s, as shown in Fig. 4. An analogous, but opposite, argument can be used to show why the wavelength windows of ultraflattened dispersion behavior get reduced as we search configurations with increasingly negative dispersion, as it is apparent in Fig. 4 as well. In any case, notice the remarkably wide windows of ultraflattened dispersion behavior. These windows are defined according to the usual criterion used to define flattened dispersion, namely, that the maximum of the dispersion variation in the given window has to be lesser than a fixed small amount (in our case, 2psnm

^{-1}km

^{-1}).

*D*

_{m}curve behavior that becomes less and less smooth as we decrease λ due to its growing curvature. The slope of this curve increases quickly yielding smaller wavelength regions where the overlap of the intervals of linear behavior of -

*D*

_{m}and

*D*

_{g}can occur. So, it is clear that there exist a limit for ultraflattened behavior as we move toward shorter wavelengths. If λ is small enough, the behavior of the material dispersion is so badly represented by a linear approximation that there is not even the possibility to work with the same design procedure explained in the previous section to obtain ultraflattened dispersion profiles. This is the case of the

*Ti-Za*wavelength window centered at 0.8

*µ*m.

*D*

_{m}for silica at 0.8

*µ*m is approximately 120 psnm

^{-1}km

^{-1}, a high value to compensate with

*D*

_{g}if one is looking for positive or nearly zero total dispersion. The second one is that, in this wavelength range, the curvature of the geometrical and material curves, unlike in the 1.55

*µ*m region, has always opposite signs. This fact is clearly appreciated in Fig. 2, where the curvature of Dg is negative around 0.8

*µ*m in all cases, whereas, according to Fig. 3, the curvature of -

*D*

_{m}remains positive even for values of λ beyond the zero material dispersion point at 1.3

*µ*m. The issue now is not to play with the slope of

*D*

_{g}, as before, but to be able to achieve values of the geometrical dispersion large enough to compensate for the high value of the sign-changed material dispersion at this wavelength. Our strategy will consist in finding configurations whose geometrical dispersion around 0.8

*µ*m exceeds the value of -

*D*

_{m}in the same wavelength region, as shown in Fig. 5. Because of the opposite sign of the curvature of these two curves, it is granted that the profile of the total dispersion will have the form shown by the red curve in Fig. 5. It will include one point of zero third-order dispersion, located at the wavelength for which the difference between the

*D*

_{g}and -

*D*

_{m}curves reaches a maximum, or, equivalently, at the point for which the negative slopes of both curves are equal. This type of behavior has been already proven to exist in PCF’s, although in a different wavelength window and only for nearly-zero dispersion [5

5. A. Ferrando, E. Silvestre, J. J. Miret, J. A. Monsoriu, M. V. Andrés, and P. S. J. Russell, “Designing a photonic crystal fibre with flattened chromatic dispersion,” Electron. Lett. **24**, 325–327 (1999). [CrossRef]

*M*=1, or Λ=2.3

*µ*m, in Fig. 2 (a)), whose maximum occurs close to 0.8

*µ*m. At this wavelength, the geometrical dispersion has a small value (

*D*

_{g}≈ 25 psnm

^{-1}km

^{-1}) as compared to that of -

*D*

_{m}at the same wavelength. The properties depicted in Fig. 2 for

*D*

_{g}will guide us in increasing the values of the geometrical dispersion. According to Fig. 5, we will focus on the region of these curves that have negative slope and near maximum. It is clear from Fig. 2(b) that, as we increase

*f*, the maximum of the

*D*

_{g}curve moves upwards and simultaneously is shifted to the right. Despite that we are able to increase the value of the maximum of the

*D*

_{g}curve, this maximum moves away from the 0.8

*µ*m window. We can relocate the

*D*

_{g}curve in such a way the region near maximum moves back to the desired window by acting nowon the magnification

*M*. By reducing

*M*, we simultaneously displace the maximum to shorter wavelengths and increase its value, as depicted in Fig. 2(a). The global effect on the value of the dispersion of this twofold operation is additive, so that we can considerably increase the value of the geometrical dispersion in the 0.8

*µ*m window by a suitable selection first of

*f*(increase) and then of

*M*(decrease). The high value of -

*D*

_{m}at 0.8

*µ*m can be in fact overcome, as shown by the positive dispersion curve in Fig. 6.

*µ*m (Fig. 5). Another complementary question we can formulate now is whether it is possible to obtain also flattened dispersion configurations in windows centered at 1.5

*µ*m. The answer to this question is certainly yes. It is not surprising because, according to what it has been discussed before for the 0.8

*µ*m window, the requirement for flattened dispersion behavior is less demanding than that necessary to achieve an ultraflattened dispersion profile. Therefore, it is possible to analyze the interesting technological issue of the tunableness of flattened dispersion in the telecommunication window centered at 1.5

*µ*m. In another words, we search all type of flattened dispersion configurations that can have positive, negative, or nearly-zero dispersion centered at this wavelength. In this way, we will be able to identify the PCF’s geometrical parameters that permit to tailor an specified flattened dispersion profile in this window. This is precisely the analysis depicted in Fig. 7. We can appreciate here that there is a considerable range of tunableness in the dispersion of a PCF in this window. This range extends from configurations owning positive dispersion (up to +45 ps nm

^{-1}km

^{-1}) to configurations with similar but negative value of dispersion (up to -43 psnm

^{-1}km

^{-1}). In all cases, the wavelength extension of the flattened dispersion behavior is around or above 200 nm.

## 4 Conclusions

*Ti-Za*windows. We have focused on a specially technologically interesting type of dispersion profiles, namely, that corresponding to configurations owning constant dispersion over wide wavelength windows. In this direction, we have formulated two different design procedures depending on whether we are interested to achieve flattened dispersion (one point of zero third-order dispersion) or ultraflattened dispersion (one point of zero fourth-order dispersion). In all cases, the basic triangular geometry of the PCF cladding has proven to be very rich in yielding configurations covering a large dispersion spectrum, ranging from rather large positive values to equally large negative values of dispersion. The final conclusion of our systematic analysis on PCF dispersion is that, despite the enormous size of the parameter space, a good comprehension of the guidance mechanism permits an extreme simplification of the search procedure, resulting in the implementation of a successful design prescription to achieve constant dispersion configurations.

## References and links

1. | 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. |

2. | D. Mogilevtsev, T. A. Birks, and P. S. J. Russell, “Dispersion of photonic crystal fibers,” Opt. Lett. |

3. | M. J. Gander, R. McBride, J. D. C. Jones, D. Mogilevtsev, T. A. Birks, J. C. Knight, and P. S. J. Russell, “Experimental measurement of group velocity dispersion in photonic crystal fibers,” Electron. Lett. |

4. | P. J. Bennet, T. M. Monro, and D. J. Richardson, “Toward practical holey fiber technology: fabrication, splicing, modeling, and fabrication,” Opt. Lett. |

5. | A. Ferrando, E. Silvestre, J. J. Miret, J. A. Monsoriu, M. V. Andrés, and P. S. J. Russell, “Designing a photonic crystal fibre with flattened chromatic dispersion,” Electron. Lett. |

6. | J. Broeng, D. Mogilevtsev, S. E. Barkou, and A. Bjarklev, “Photonic crys tal fibers : a new clas s of optical waveguides,” Opt. Fib. Tech. |

7. | A. Ferrando, E. Silvestre, J. J. Miret, and P. Andrés, “Nearly zero ultraflattened dispersion in photonic crystal fibers,” Opt. Lett. |

8. | E. Silvestre, M. V. Andrés, and P. Andrés, “Biorthonormal-basis method for the vector description of optical-fiber modes,” J. Lightwave Technol. |

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. |

10. | A. Ferrando, E. Silvestre, J. J. Miret, P. Andrés, and M. V. Andrés, “Vector description of higher-order modes in photonic crystal fibers,” J. Opt. Soc. Am. A |

11. | D. Davidson, Optical-Fiber Transmission (E. E. Bert Basch, ed., Howard W. Sams & Co, 1987). |

**OCIS Codes**

(060.2270) Fiber optics and optical communications : Fiber characterization

(060.2280) Fiber optics and optical communications : Fiber design and fabrication

(060.2430) Fiber optics and optical communications : Fibers, single-mode

**ToC Category:**

Focus Issue: Photonic crystal fiber

**History**

Original Manuscript: November 5, 2001

Published: December 17, 2001

**Citation**

Albert Ferrando, Enrique Silvestre, Pedro Andres, Juan Miret, and Miguel Andres, "Designing the properties of dispersion-flattened photonic crystal fibers," Opt. Express **9**, 687-697 (2001)

http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-9-13-687

Sort: Journal | Reset

### References

- 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]
- D. Mogilevtsev, T. A. Birks, and P. S. J. Russell, "Dispersion of photonic crystal fibers," Opt. Lett. 23, 1662-1664 (1998). [CrossRef]
- M. J. Gander, R. McBride, J. D. C. Jones, D. Mogilevtsev, T. A. Birks, J. C. Knight, and P. S. J. Russell, "Experimental measurement of group velocity dispersion in photonic crystal fibers," Electron. Lett. 35, 63-64 (1999). [CrossRef]
- P. J. Bennet, T. M. Monro, and D. J. Richardson, "Toward practical holey fiber technology: fabrication, splicing, modeling, and fabrication," Opt. Lett. 24, 1203-1205 (1999). [CrossRef]
- A. Ferrando, E. Silvestre, J. J. Miret, J. A. Monsoriu, M. V. Andres, and P. S. J. Russell, "Designing a photonic crystal fibre with flattened chromatic dispersion," Electron. Lett. 24, 325-327 (1999). [CrossRef]
- J. Broeng, D. Mogilevtsev, S. E. Barkou, and A. Bjarklev, "Photonic crystal fibers: a new class of optical waveguides," Opt. Fib. Tech. 5, 305-330 (1999). [CrossRef]
- A. Ferrando, E. Silvestre, J. J. Miret, and P. Andres, "Nearly zero ultraflattened dispersion in photonic crystal fibers," Opt. Lett. 25, 790-792 (2000). [CrossRef]
- E. Silvestre, M. V. Andres, and P. Andres, "Biorthonormal-basis method for the vector description of optical-fiber modes," J. Lightwave Technol. 16, 923-928 (1998). [CrossRef]
- 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]
- A. Ferrando, E. Silvestre, J. J. Miret, P. Andres, and M. V. Andres, "Vector description of higher-order modes in photonic crystal fibers," J. Opt. Soc. Am. A 17, 1333-1340 (2000). [CrossRef]
- D. Davidson, Optical-Fiber Transmission (E. E. Bert Basch, ed., Howard W. Sams & Co, 1987).

## 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.

« Previous Article | Next Article »

OSA is a member of CrossRef.