Abstract

An all-solid dispersion-compensating photonic crystal fiber based on mode coupling mechanism in dual-concentric core has been proposed. The mode coupling characteristics, dispersion, confinement loss of the fiber, and the influence on dispersion of some structure parameters are simulated by full-vector finite element method. By using the relationship between phase matching wavelengths and coupling strength with the change of fiber microstructure parameters, an all-solid dual-concentric-core dispersion-compensating photonic crystal fiber is presented. The structure parameters on dispersion characteristic are investigated. The results demonstrate that the proposed fiber has a large negative dispersion value 8465 ps/(nm·km) at 1550 nm. The effective mode area and the splicing loss to the standard single mode fiber are 12.8 μm2 and 1.89 dB at 1550 nm, respectively. At 1550 nm, the confinement loss is less than 1 × 10−3 dB/km and the bending loss with 2 cm bending diameter is less than 1 × 10−2 dB/km.

1. Introduction

In modern high-speed optical fiber communication system, the cumulative dispersion in transmission process is the main factor limiting the communication rate [1]. Therefore, the dispersion compensation for high-speed optical fiber communication system has become a hot research topic. Photonic crystal fiber (PCF) cladding contains a lot of micro air holes with periodic or disordered distribution. Because of this special structure, this kind of fiber has a high degree of design freedom, and it can flexibly adjust the dispersion characteristics of fiber [2]. It allows to arbitrarily shape the dispersion in a wide spectral range [3].

Birks et al. got the dispersion compensation photonic crystal fiber (DCPCF) for the negative dispersion value as −2000 ps/(nm·km) at 1550 nm by using control effect of large air hole and small fiber core on light to increase the influence of waveguide dispersion [4]. Gerome et al. proposed the structure of dual-concentric-core dispersion-compensating photonic crystal fiber (DCC-DCPCF) [5] and used the mode coupling mechanism of dual-concentric-core (DCC) to obtain the large negative dispersion. After then, the design and realization of dispersion compensation characteristics by using microstructural fiber are mostly based on the mode coupling mechanism of DCC [69]; the transmission mode between the inner and outer core is regulated to achieve the conversion of light energy between the inner and outer core, so as to flexibly adjust the dispersion compensation characteristics of the optical fiber. All the fibers mentioned above are air holes or liquid-filled air holes. A liquid-filled hybrid structure of DCC-DCPCF with the negative dispersion value being −40400 ps/(nm · km) at 1550 nm and an inner core high-doping DCC-DCPCF with the negative dispersion value being −51625 ps/(nm · km) at 1550 nm were reported by Hsu et al. [10, 11]. Partha et al. designed a DCC-DCPCF with liquid-filled air holes—the negative dispersion value has reached −52100 ps/(nm · km) at 1550 nm—and analyzed the influence of the material on the phase matching wavelength (λp) and dispersion value at different temperatures [12]. Han et al. developed a negative dispersion value of −41800 ps/(nm · km) at 1550 nm in a DCC-DCPCF by using two methods to expand the first-layer air hole [13]. DCC-DCPCF introduced above are air hole silica or liquid-filled air holes structure, which have disadvantages of complex drawing technique and not easy to be coupled with standard single-mode fiber (SSMF).

Liquid-filled air holes fiber belongs to the fiber postprocessing technology; it is very difficult to achieve liquid filling due to the DCC-DCPCF has a small structure, and this kind of fiber is also not conducive to the application of practical engineering. An all-solid PCF with low dispersion slope and an all-solid PCF with a wide wavelength bands’ dispersion compensation were reported by Xu et al. [14, 15]. Mariusz et al. used an all-solid PCF with flattened normal dispersion to improve spectral flatness of supercontinuum [16]. Wang et al. developed an all-solid dual-concentric-core microstructure fiber for ultra-broadband dispersion compensation, the designed fiber has very wide compensation bandwidth of up to 190 nm. [17]. The all-solid PCF consists of low Ge-doped inner core and F-doped concentric ring inner cladding.

In this paper, based on the mode coupling mechanism between the traditional air hole silica DCC-DCPCF, the mode coupling characteristics of the DCC all-solid photonic crystal fiber are studied by the full-vector finite element method. The influence of the different structure parameters on the mode coupling characteristics is analyzed. The structure parameters on dispersion characteristic are investigated. Based on that, we optimize the fiber structure parameters reasonably and design an all-solid DCC-DCPCF. In addition, the all-solid structure provides a more friendly process for both fabricating and splicing; the designed all-solid DCC-DCPCF is easily fabricated and can easily access the existing optical network systems.

2. Fundamental Theory

A full-vector finite element method with the perfectly matched layer boundary conditions is used to analyze the mode coupling mechanism of all-solid DCC-DCPCF [18].

Starting with Maxwell’s curl equations, the vector equation for the electric field vector E can be derived aswhere the subscript i is section number, k0 = /λ is the wave number in the vacuum, and λ is operation wavelength.

The curvilinear hybrid edge/nodal elements based on linear tangential and quadratic normal vector basis functions are adopted to realize the computational window divisions and perfectly matched layer is incorporated as the boundary condition to absorb waves out of the computational window. Applying the finite element procedure to equation (1), the following eigenvalue equation is defined aswhere [A] and [B] are the global finite element matrices. The eigenvector {H} and the eigenvalue provide, respectively, the full vector magnetic field distribution on the cross section of PCF and the effective index of the mode.

The total dispersion D (λ) of the fiber consists of two parts: the material dispersion Dm (λ) and the waveguide dispersion Dw (λ); that is [19],

The waveguide dispersion can be expressed as [20]where c is the velocity of light in vacuum, λ is the wavelength, and Re [neff] is the real part of the effective refractive index of the fundamental mode.

The material dispersion is defined as [21]where nm is the refractive index of the material, and it is obtained by Sellmeier’s equation at different wavelengths.

The effective mode field area of the fiber can be expressed as [22]where E is an electric field vector and S represents the mode field area of the fiber.

The curved PCF can be replaced by a straight fiber with equivalent refractive index as [23]where R is the curvature diameter and is the original refractive index profile of the PCF.

The confinement loss is obtained from the imaginary part of neff by

The splice loss between the proposed DCC-DCPCF and the SSMF is theoretically evaluated according to the following definition [24]:where ω1 and ω2 represent the mode field diameter of the DCC-DCPCF and SSMF.

3. Structural Design and Simulation Analysis

3.1. All-Solid DCC-DCPCF

The transverse section and the equivalent refractive index of the all-solid DCC-DCPCF are shown in Figure 1. The base material is pure silica, and the refractive index n0 = 1.45. The fiber consists of a high-germanium-doped silica column as an inner core and a low-doped fluoride or boron silica column as the outer core or cladding and the refractive index of the doped rod is controlled within ±0.05. It can be seen that the cladding structure contains three kinds of different size and different refractive index of medium cylinders, whose center medium cylinder is the inner core; its refractive index is n1 and diameter is d1; the outer core is the minimum media cylinder; its refraction index and diameter is n2 and d2, respectively; gray medium column is cladding and the diameter is d3. The cladding refraction index contents n3 < n2< n0< n1. The column spacing between any two silica columns is Λ, the number of low-doped silica columns between the inner core and outer core is N, and the number of layers after the outer core is No. The silica columns of PCF have a hexagonal distribution. The first ring low-doped silica columns has 6 columns, the second ring low-doped silica columns has 12 columns, and the Nth ring low-doped silica columns has 6 ∗ N columns. As for outer core, the low-doped silica columns number of the Noth is 6 ∗ (No + N). For every two adjacent layers in outer cladding, there are six more columns in the outer layer than in the inner layer. There are no air holes in all-solid DCC-DCPCF, which make it easier to fabricate such fiber than conventional air-hole PCF. As long as the diameter of columns is not low than 0.10 μm, it can be fabricated.

3.2. Mode Coupling Characteristics of the All-Solid DCC-DCPCF

To study the mode coupling characteristics of the all-solid DCC-DCPCF, a set of initial parameters is given: n1 = 1.480, n2 = 1.445, n3 = 1.42, d1 = 1.30 μm, d2 = 0.51 μm, d3 = 2.40 μm  = 2.90 μm, N = 2, and No = 2. The relationship between the wavelength and the effective refractive index of the inner core fundamental mode and outscore second-order mode is calculated by the full-vector finite element method.

When the optical fibers are designed to have only inner core or outer core, the relationship between the wavelength and the effective refractive index of the inner core or outer core fundamental mode is shown in Figure 2. We can see that the effective refractive index of the inner core and outer core fundamental mode is approximately linear with the wavelength when it only has inner core or outer core; the effective refractive index of the inner core fundamental mode is greater than the outer core at the short wave, while it is the opposite at the long wave, the two effective refractive index curves have intersection points at a fixed wavelength. But if the structure of designed fiber is DCC, the effective refractive index of the inner core and outer core mode is abrupt in a certain wavelength, and the distribution of the mode field is also changed, which have the two modes of the outer core and inner core distribution, respectively.

According to the coupled mode theory, the core mode of DCC fiber affects each other by evanescent wave. The DCC fiber forms a super mode with guide characteristics by controlling of mode coupling mechanism between the inner core and outer core. The inner core mode and the outer core mode in Figure 2 are called fundamental supermodel and second-order supermodel, respectively. At the short wavelength, the fundamental supermodel is transmitted in the inner core, and the second-order supermodel is transmitted in the outer core. Moreover, most of the energy of the fundamental supermodel and the second-order supermodel is gathered in the inner core and the outer core, respectively. With the increase of wavelength, the effective refractive index of the two modes has a sudden change in a certain wavelength. This wavelength is called λp. At the λp, the energy of the fundamental supermodel is gradually converted to the outer core, while the energy distribution of the second-order supermodel converts to the inner core gradually. After λp, the energy of the fundamental supermodel is more and more transferred from the inner core to the outer core, and the energy of the second-order supermodel is more and more transferred from the outer core to the inner core.

The no-cross effect in the coupling process makes the effective refractive index curve of the fundamental supermodel and the second-order supermodel rapidly change, and the two modes produce large negative dispersion and positive dispersion, respectively, at λp, as shown in Figure 3. Owing to the fact that SSMF has certain positive dispersion in common communication band, the large negative dispersion produced by the fundamental supermodel can compensate the cumulative dispersion in the transmission process. Due to the large loss of the second-order supermodel, it will decay continuously during transmission, so we do not need to consider the influence of second-order supermodel on communication.

3.3. Simulation Analysis

According to the refractive index distribution of the fiber transverse section, a change on the diameter of the medium cylinder is the equivalent to a change on the refractive index of the medium cylinder. The influence of the diameter of the medium cylinder on mode coupling characteristics is analyzed. When the diameter of inner core medium cylinder is changed, the relationship between the effective refractive index of the fundamental supermodel and the wavelength is shown in Figure 4. We can see when d1 is changed, the effective refractive index of the fundamental supermodel at long wavelength is almost unchanged; the effective refractive index at short wavelength increases with the increase of d1. At d1 = 0.90 μm, the effective refractive index of the fundamental supermodel has a distinct change-point, which is the phase matching point. With the increase of d1, the abrupt change of the effective refractive index of the fundamental supermodel gradually becomes slow, but the mutation position can still be seen.

Shown in the insets of Figure 4(a), the fundamental supermodel is transmitted in the inner core at the short wavelength; most of the energy of the fundamental supermodel is gathered in the inner core for d1 = 0.90 μm. At the long wavelength, the fundamental supermodel is transmitted in the outer core; most of the energy of the fundamental supermodel is gathered in the outer core; that is, most of the energy is transferred from the inner core to the outer core. With d1 = 1.70 μm, most of the energy of the fundamental supermodel is also gathered in the inner core at the short wavelength, but at the long wavelength, the fundamental supermodel still has a large part of the energy gathered in the inner core, which means the energy transfer from the inner core to the outer core is slow. Thus, when the effective refractive index of the fundamental supermodel changes drastically from short wavelength to long wavelength, there is an obvious phase matching point. The energy of the fundamental supermodel is transferred rapidly from the inner core to the outer core near the phase matching point; that is, the supermodel between the inner core and outer core has a strongly coupling. When the effective refractive index of the fundamental supermodel changes slowly from short wavelength to long wavelength, there is no obvious phase matching point. The energy of the fundamental supermodel is transferred slowly from the inner core to the outer core near the phase matching point that means the supermodel between the inner core and outer core has a weak coupling. When the equivalent refractive index difference between the inner core and outer core increases, the effective index difference of the fundamental supermodel in the inner core and outer core is also increased. When the effective refractive index difference is small, the coupling process of the fundamental supermodel from inner core to outer core is fast and the coupling strength is large, producing a larger negative dispersion value. When the effective refractive index difference is large, the coupling process is slow and the coupling strength is weak, so a smaller negative dispersion value is generated. Hence, the coupling strength between the inner core and outer core can be expressed by the absolute value of the negative dispersion value, and the greater the absolute value is, the greater the coupling strength is. The variation of λp and coupling strength with d1 is shown in Figure 4(b). λp produces red shift and the coupling strength of the fundamental supermodel decreases gradually with the increasing of d1, which is consistent with the theoretical analysis.

Variation of the effective refractive index of the fundamental supermodel with wavelength when d2 is changed can be seen in Figure 5(a). The effect of d2 change and d1 change is exactly the opposite. The effective refractive index of the fundamental supermodel at the short wavelength is basically unchanged, but the effective refractive index at the long wavelength is increased with the decrease of d2. In fact, the equivalent refractive index of outer core becomes larger and the equivalent refractive index difference between inner core and outer core decreases with the increase of d2. It is equivalent to the equivalent refractive index difference of the fundamental supermodel transferred from the inner core to the outer core becomes smaller, which leads to a rapid coupling process of the fundamental supermodel from the inner core to the outer core and a larger coupling strength. The variation of λp and coupling strength with d2 is presented in Figure 5(b). Figure 5(b) shows that λp produces red shift and the coupling strength of the fundamental supermodel decreases little by little with the increasing of d2.

Variation of the effective refractive index of the fundamental supermodel with wavelength when d3 is changed is analyzed in Figure 6(a). Decreasing d3 increases the equivalent refractive index of the cladding, which is equal to the increase of the effective refractive index of the fundamental supermodel at the short wavelength and the long wavelength with the decrease of d3, but the increase at the short wavelength is less than that at the long wavelength. This result is different from that which only changes the core area effective refractive index of designed fiber. At the same time, we can see that the growth amplitude of the fundamental super mode’s equivalent refractive index increases with the decrease of d3 at short wavelength and long wavelength. So, the change of d3 has an influence on the change of the refractive index of inner core and outer core, and the changes of each other are different. Finally, the variation of λp and coupling strength of the fundamental supermodel is modulated by the refractive index change of the inner core and outer core. The variation of λp and coupling strength with d3 is expressed in Figure 6(b). It can be seen in Figure 6(b) that λp is red shifted first and then blue shifted, but the coupling strength of the fundamental supermodel has been increased, which further shows that this result is adjusted by the refractive index change of the inner core and outer core.

Variation of the effective refractive index of the fundamental supermodel with wavelength when is changed is discussed in Figure 7(a). The effective refractive index of the fundamental supermodel at the short wavelength and the long wavelength increases with the increase of , but the increase at the short wavelength is less than that at the long wavelength. Meanwhile, the growth amplitude of the equivalent refractive index of the fundamental supermodel at the short wavelength and the long wavelength decreases with the increase of , respectively. Changing is equal to changing the diameter of the cladding and core area, further changing the equivalent refractive index distribution in the cladding and core areas. The changing of equivalent refractive index of cladding also has an influence on the equivalent refractive index distribution in the core areas. The final variation of λp and coupling strength of the fundamental supermodel is modulated by the variation of the equivalent refractive index of the cladding and core area. The variation of λp and coupling strength with is shown in Figure 7(b). λp is red shifted and the coupling strength of the fundamental supermodel gradually weakened with the increase of . Because of the limitation of the structure of the fiber, has a small change, so λp and the coupling strength of the fundamental supermodel produce a single trend under the combined action of the equivalent refractive index of core area and cladding.

The fabrication difficulty of the photonic crystal fiber increases with the increase of the number of layers; we study the number of layers N = 1–3. Variation of the effective refractive index of the fundamental supermodel with wavelength when N is changed is shown in Figure 8(a). When the number of layers N varies from 1 to 3, the effective refractive index of the inner core basically remains unchanged, while the refractive index of the outer core increases gradually. But from 2 to 3 layers, the increase of the effective refractive index decreases obviously. The variation of λp and coupling strength with N is shown in Figure 8(b). λp is red shifted and the coupling strength becomes larger with the increase of the number of layers N, while λp is less red shifted and the coupling strength has a great change with N changed from 2 to 3. Therefore, considering the fabrication difficulty of the photonic crystal fiber and the change of the coupling strength with N, the number of layers N = 3 is set in the future research.

3.4. The Result Optimization of the All-Solid DCC-DCPCF

According to the above study of the influence of the fiber structure parameters on λp and the mode coupling strength, an all-solid DCC-DCPCF is designed for n1 = 1.480, n2 = 1.445, n3 = 1.430, d1 = 1.30 μm, d2 = 0.51 μm, d3 = 3.40 μm, Λ = 3.50 μm, N = 3, and No = 3. The variation of the dispersion with wavelength is shown in Figure 9. λp of the designed fiber locates at 1550 nm, and the negative dispersion value reached −8465 ps/(nm · km). At 1550 nm, is 12.8 μm2 and is 1.89 dB. The confinement loss and bending loss of the fiber can be controlled by increasing the number of cladding layers. The confinement loss and bending loss with 2 cm bending radius are less than 1 × 10−3 dB/km and 1 × 10−2 dB/km at 1550 nm, respectively, when No = 6. The field is mostly confined to the inner core for the optimized DCC-DCPCF. The proposed all-solid DCC-DCPCF not only has large negative dispersion at 1550 nm and higher compensation ratio, but also has the advantages of low confinement loss and easy fabrication.

4. Conclusions

In summary, an all-solid DCC-DCPCF with large negative dispersion at 1550 nm has been presented. The influence of the different structure parameters on the mode coupling characteristics is analyzed. The results demonstrate that the proposed all-solid DCC-DCPCF has a large negative dispersion value and low confinement loss. By optimizing the fiber structure parameters, the dispersion value at 1550 nm of the proposed all-solid DCC-DCPCF researches −8465 ps/(nm · km), and the compensation ratio of the proposed fiber to SSMF is 500. The proposed all-solid DCC-DCPCF has important guiding significance for the industrialization of dispersion compensation.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (nos. 61405173 and 61405172) and the Hebei Provincial Natural Science Foundation of China (no. F2018203346).