1. Introduction
Power scaling of the fiber laser technology has been quite impressive in the recent years. In pulsed regime, multi-megawatt peak power has been obtained for nanosecond pulses in large core fibers operating beyond the single mode regime. Peak power is limited in singlemode fibers by either power handling of the core or nonlinear effects. Special structures have been proposed to decrease the effective non-linear gain in order to scale to higher power. For instance, W like waveguide structures can be designed to reduce the overlap between signal and Raman Stokes wavelength [
1
J. Kim, P. Dupriez, C. Codemard, J. Nilsson, and J. K. Sahu, “Suppression of Stimulated Raman Scattering in a high power Yb-doped Yb-doped fiber amplifier using a W-type core with fundamental mode cut-off,” Opt. Express
14(12), 5103–5112 (2006). [CrossRef]
[PubMed]
]. This solution increases the Raman non-linear threshold but does not address other non-linearities. Moreover the damage threshold is unaffected by this technique. Another approach is then to increase the mode field diameter (MFD) while decreasing the numerical aperture (NA). Standard large mode area (LMA) single mode fibers have a diameter typically limited to 30 µm with 0.06 NA for single mode operation. Photonic crystal fibers (PCF) can be designed to operate in the endlessly single mode regime but are limited by bending losses in the same range of core diameter [
2
J. C. Baggett, T. M. Monro, K. Furusawa, V. Finazzi, and D. Richardson, “Understanding bending losses in holey optical fibers,” Opt. Commun.
227(4-6), 317–335 (2003). [CrossRef]
]. Recently, leakage channel fibers have been introduced to address these limitations [
3
W. S. Wong, X. Peng, J. M. McLaughlin, and L. Dong, “Breaking the limit of maximum effective area for robust single-mode propagation in optical fibers,” Opt. Lett.
30(21), 2855–2857 (2005). [CrossRef]
[PubMed]
]. These fibers are formed by a core surrounded by a single ring of large holes. These fibers require the core refractive index to be equal to the silica index. For fibers operating at eyesafe wavelengths (1.5 and 2µm) this is typically more difficult than for 1 µm amplification because of high dopant index. In 1.5µm Erbium-Ytterbium fibers, a large amount of phosphorous is used which increases the refractive index. In the same manner, in 2 µm Thulium fibers high efficiency can be obtained at high doping level. But the core refractive index is increased accordingly. Pedestal fibers allow to compensate for this larger NA but their NA is still higher than 0.1 [
4
K. Tankala and B. Samson., “New developments in High Power Eye-Safe LMA fibers,” Proc. SPIE
6102, 610206 (2006). [CrossRef]
]. There are thus still many efforts to propose special fiber designs with large mode effective area operating in the single mode regime or close to this regime and with low bending losses.
Recently, such a new type of large mode fiber has been proposed. Its core is composed of high index filaments surrounded by undoped silica (
Fig. 1
). An Erbium-Ytterbium doped 37-filaments fiber has been prepared and tested in MOPA configuration [
5
G. Canat, S. Jetschke, L. Lombard, S. Unger, P. Bourdon, J. Kirchhof, A. Dolfi, V. Jolivet, and O. Vasseur, “Multifilament-core fibers for high energy pulse amplification at 1.5 µm with excellent beam quality,” Opt. Lett.
33(22), 2701–2073 (2008). [CrossRef]
[PubMed]
]. The fiber core size was 24 x 32 µm. The amplifier operated in a quasi-single mode regime with M
2 = 1.3. More recently, a passive 19-filaments fiber has been prepared and tested in the single mode regime [
6
M. M. Vogel, M. Abdou-Ahmed, A. Voss, and T. Graf, “Very-large-mode-area, single-mode multicore fiber,” Opt. Lett.
34(18), 2876–2878 (2009). [CrossRef]
[PubMed]
]. Both experiments show that the MFC fibers exhibit low bending losses at bending radius down to 15 cm. This fiber structure is thus promising. However, the source of these good performances remained unclear. These multifilament fibers can be seen as a version of the multicore fibers with core radius a≤λ. In these fibers, there is a strong coupling between the filaments which enables only a few supermodes to propagate in the core. This situation is different of usual multicore fibers in which the coupling between the cores is low: usual multicore fibers with N cores guide N supermodes which can be studied using the coupled mode theory.
Fig. 1 Cleave of a 28 x 31 µm multifilament core fiber with N = 37 filaments.
In this paper, we propose a theoretical analysis of the multifilament core fibers. Our theory explains the most significant properties of the MFC. We first developed a finite element modelling (FEM) of the MFC fiber that supports the good experimental results. This model is presented in part 2. This method provides little insight into the fiber properties and is computationally demanding. In part 3, we show that the equivalent step index fiber theory (ESI) fails to explain the fiber performances. In part 4, we model the MFC fiber by an improved equivalent step index fiber with a core refractive index dependent on the wavelength. The equivalent core index is equal to the so-called fundamental space filling mode (FSM), similar to the one used in photonic crystal fibers (PCF) [
7
T. A. Birks, J. C. Knight, and P. S. Russell, “Endlessly single-mode photonic crystal fiber,” Opt. Lett.
22(13), 961–963 (1997). [CrossRef]
[PubMed]
]. This semi-analytical model will be referred as the FSM based step index method (ESI-FSM). It is shown to agree with the FEM computations. Part 5 uses this ESI-FSM model to predict the fiber performances.
2. The multifilament core (MFC) fiber
The MFC structure can be described as a microstructured core surrounded by a silica cladding with refractive index n
clad = 1.444 (
Fig. 2
). The core has a periodic structure composed of Erbium-Ytterbium doped filaments with radius a and refractive index n
1 surrounded by a fluorine doped glass with refractive index n
2. The fiber that has been tested has a hexagonal structure with pitch Λ = 5.1 µm. Through this paper we will consider the structure of
Fig. 2 with a = 0.9 µm and Λ = 5.1 µm as well as n
1 = 1.455 and n
2 = 1.443 corresponding to Erbium-Ytterbium filaments and fluoride doped glass. The fiber on
Fig. 2 has 37 filaments. On an hexagonal lattice, N filaments can be packed into P rings in a regular MFC if N = 1 + 3P(P + 1).
Fig. 2 The multifilament fiber core with hexagonal lattice. The high index filament (red disks) radius is a and the lattice pitch is Λ.
Modelling can be performed using vectorial finite elements method. Using the commercial software COMSOL Multiphysics, we were able to compute the variations of the fundamental mode effective index
and effective area with wavelength (
Fig. 3
). This method will be used as reference through this paper. At 1.5 µm, the effective area reaches more than 800 µm
2.
Fig. 3 Evolution of the effective index and effective area of the fundamental mode as a function of wavelength computed with the finite element method.
The bending losses α for the power propagating in the fundamental mode in a fiber with bending radius R can be computed by surrounding the fiber by a perfectly matched cladding layer (PML). The bending of the fiber was taken into account by transforming the index profile n(r) to n
bend(r) with [
8
D. Marcuse, “Influence of curvature on the losses of doubly clad fibers,” Appl. Opt.
21(23), 4208–4213 (1982). [CrossRef]
[PubMed]
], [
9
R. T. Schermer and J. H. Cole, “Improved bend loss formula verified for Optical Fiber by simulation and experiment,” IEEE J. Quantum Electron.
43(10), 899–909 (2007). [CrossRef]
]:
where ξ = 1.27 accounts for the bend induced stress in the fiber.
The propagation constant β of the fundamental mode is complex: its imaginary part accounts for attenuation. The propagation losses are computed using the following equation:
As it was pointed out the multifilament core fibers exhibit low sensitivity to bending considering the very large area of the fundamental mode (
Fig. 4
). At 1545 nm for instance, the effective area is A
eff = 855 µm
2 and the critical bending radius where the bending loss reaches 1 dB/m is around 20 cm.
Fig. 4 Modelling of the bending loss as a function of the bending radius for the MFC fiber using the finite element method (circles) and the equivalent step index fiber from moments (squares).
3. The equivalent step index fiber (ESI) using the theory of moments
To draw a comparison, we can make calculations on an equivalent step index (ESI) fiber. The constant refractive index in the equivalent core will be an averaged version of the MFC fiber refractive index profile. At first sight, the best equivalent step index fiber can be designed using the two first moments of the refractive index distribution [
10
R. J. Black and C. Pask, “Developments in the theory of equivalent-step-index fibers,” J. Opt. Soc. Am. A
1(11), 1129–1131 (1984). [CrossRef]
]. We define the moments G
0 and G
2p of the refractive index profile n(r,θ) by
It can be shown that in the limit λ→∞, two fibers have the same effective index of the fundamental mode if they have the same moments G
p of all orders [
10
R. J. Black and C. Pask, “Developments in the theory of equivalent-step-index fibers,” J. Opt. Soc. Am. A
1(11), 1129–1131 (1984). [CrossRef]
]. In the long wavelength limit, this equivalence is exact. For a step index fiber with core radius a and core index n
core it implies:
We can thus compare the MFC fiber with a step index fiber with core radius a
step and numerical aperture NA
step given by
and
For the MFC of
Fig. 2, G
0 = 0.138 µm
2 and G
2 = 17.1 µm
4. We thus get
= 15.74 µm and
= 0.034. The corresponding equivalent core index is
= 1.444386. In this ESI fiber, we can compute the effective index of the fundamental mode using the characteristic equation for step index fibers (
Fig. 5
) [
11A. W. Snyder, and J. D. Love, “Optical Waveguide Theory”, Chapman and Hall, New York, 1983.
]. We observe that the effective index predicted by the equivalent step index fiber is always far below the FEM effective index. At 1545 nm for instance, the effective index using the ESI fiber is
= 1.444177 whereas
= 1.444361 computed using the FEM solver. At large wavelengths, the field only sees the averaged index profile and both values converge.
Fig. 5 Variation of the effective index of the fundamental mode in the MFC and in the equivalent step index fiber (ESI). The dashed line corresponds to λ = 1545nm where the bending loss modeling has been performed.
The bending loss of the fundamental mode in a bent fiber has an exponential dependence on the bending radius normalized by a. critical bending radius R
crit defined by [
12
D. Marcuse, “Curvature loss formula for optical fibers,” J. Opt. Soc. Am.
66(3), 216–219 (1976). [CrossRef]
]:
with
.
This definition is quite general and can be used for arbitrary refractive index profiles with cylindrical symmetry. The much lower value of
predicted by the ESI theory leads to a value of
w lower than in the real MFC using
. This explains why we cannot compare the bending loss of the MFC and the bending loss of a step index with the same averaged core refractive index.
Figure 4 shows that the bending losses computed with this equivalent fiber are overestimated. The equivalent fiber built using the moments is an insufficient first order approximation of the real MFC profile. The real fiber structure has large index fluctuations due to the high index filament. For wavelengths with realistic values, the effective index of the core is thus wavelength dependent: at short wavelength, the light is confined into the high index filament whereas at long wavelength, the light penetration into both glasses is similar. This phenomenon is similar to the dependence of the effective cladding index in PCF on wavelength. We will thus adapt the effective index model developed for the PCF analysis with a wavelength-dependent cladding index to the MFC with a wavelength-dependent core index.
4. The equivalent step index fiber based on the fundamental space filling mode (ESI-FSM) calculation
Due to the total internal reflection mechanism, a mode will propagate in the MFC if it can travel with a propagation constant β such that the effective index n
eff = β/k is larger than n
clad. If the radius of the MFC core is increased to infinity, the light will only see a periodic composite medium with unit cell identical to the unit cell of the MFC. It will travel with an effective propagation constant β
FSM which is the fundamental space filling mode (FSM) constant [
13
M. Midrio, M. Singh, and C. Someda, “The space filling mode of holey fibers: An analytical Vectorial solution,” J. Lightwave Technol.
18(7), 1031–1037 (2000). [CrossRef]
]. In this ESI-FSM model, the equivalent core index is defined by
= β
FSM/k. For a core of finite size, we have
. The value of β
FSM can be determined by analyzing only the elementary cell of the lattice which is a hexagonal cell with a filament at the center and edge size
(
Fig. 6
). As noticed by Midrio [
13
M. Midrio, M. Singh, and C. Someda, “The space filling mode of holey fibers: An analytical Vectorial solution,” J. Lightwave Technol.
18(7), 1031–1037 (2000). [CrossRef]
], the periodicity of the lattice and the hexagonal geometry impose symmetry properties to the electric and magnetic field at the point P on the symmetry axes which can only be satisfied if the longitudinal components of the fields satisfy E
z(P) = H
z(P) = 0.
Fig. 6 The unit cell of the infinite hexagonal lattice in the fundamental space filling mode theory.
The geometry of the problem can be simplified by considering that the cell is circular with radius R. Equating the filling ratio in both geometries yields
In this cylindrical geometry, the FSM can be computed using a vectorial semi-analytical technique developed for PCF (same as in reference [
13
M. Midrio, M. Singh, and C. Someda, “The space filling mode of holey fibers: An analytical Vectorial solution,” J. Lightwave Technol.
18(7), 1031–1037 (2000). [CrossRef]
]). However, the unit cell is different: in PCF fibers a low index central air hole is surrounded by silica, whereas in the MFC a high index filament is surrounded by a low index glass. In both cases, the circular approximation imposes the boundary condition E
z(R,θ) = H
z(R, θ) = 0.
Following conventional fiber optics, we can write the electric and magnetic longitudinal components using
Eq. (11) and
(12).
In the filament (r<a):
In the low index glass (r>a):
where
are the polar coordinates,
,
,
l is an integer and A, B, C, D are constants.
Function P
l is determined from the boundary condition at r = R:
We obtain the following characteristic equation for the hybrid modes guided by the filaments.
with
where the primes denote derivatives with respect to the argument and I, J and K denote the Bessel functions. The characteristic
Eq. (14) can be solved for
. Its highest value is obtained for
l = 1 corresponding to the HE modes.
is defined to be this highest value. It decreases with wavelength (
Fig. 7
). In the limit of large λ,
converges to the geometric average index value:
with f the filling ratio
. In the short wavelength region, we get the filament refractive index value n
1. The large variation of the equivalent refractive index introduces a large core index dispersion which can be interpreted as waveguide dispersion in the core elementary cell.
Fig. 7 Variation of the core equivalent index from FSM as function of the normalized wavelength. The dashed line corresponds to the equivalent step index computed by the ESI-FSM theory.
Having computed the equivalent step index
, we need to define an equivalent core radius
. We choose
in order to have the same core area in the cylindrical equivalent core and in the hexagonal MFC. For a N filaments MFC fiber, we have:
For the fiber of
Fig. 2, we get
= 16.28 µm and at 1545 nm
= 1.444582. We can define an effective normalized frequency by the equation
Contrary to PCF where the variation of n
clad can compensate for the λ
−1 behaviour, the variation of
in MFC enhances the V variation with λ. We check that the
definition using (17) is correct by comparing the cut-off wavelengths of the LP
11 mode for the MFC and for the equivalent step index fiber for V
eff = 2.405.
5. MFC fiber properties from the effective index model
We can thus model the MFC as a step index fiber with core index n
core =
and core radius
surrounded by silica with refractive index n
clad. We can then compute the value of
using the characteristic equation of ESI-FSM fiber and compare it to the value
computed using the Finite Element Method on the MFC structure.
Figure 8
shows that an excellent agreement between both values of the effective index is obtained.
Fig. 8 Comparison of the effective index of the fundamental mode computed using the finite elements method (FEM) (squares) and the Effective Step Index Approximation (ESI-FSM) (disks) for a = 0.9 µm and Λ = 4.5 µm, 5.1 µm and 6.5 µm.
This ESI-FSM theory can predict any property of the MFC which is only related to
or
.
Equation (18) can be used to determine the single mode regime.
Figure 9
depicts the V
eff value as a function of the fiber design parameters a/Λ and λ/Λ. In step index fibers, the useful operation region is limited at high V values by the single mode limit at V = 2.405 and at the low V values by the bending losses at around V = 1.5. The equivalence between both fibers implies that the useful operating region can be extracted from our theory.
Fig. 9 Variation of the effective normalized frequency V
eff with the normalized wavelength and the a/Λ parameter. The single mode region is in the direction of the arrows. The black circle shows the fiber described in
Fig. 2.
The fundamental mode radius of this equivalent step index fiber is given by the Marcuse relation:
where V
eff is the (effective) normalized frequency given by
Eq. (18).
If we use this relation to estimate the effective area of the fundamental mode A
eff = πω
2
eqwe find results very closed to the finite elements method (
Fig. 10(a)
). The difference tends to increase at short and long wavelengths corresponding to very large and very small V
eq but it is smaller than 10%. We can expect deviation to the Marcuse relation at large V
eff as the mode shape will differ more and more from a gaussian shape.
Fig. 10 Comparison of the fundamental mode effective area (a) and bending loss (b) computed at λ = 1545 nm using the Marcuse equations and the effective index computed with the ESI-FSM compared (line) to numerical calculations using COMSOL (squares).
Bending loss for the MFC can be computed using the Marcuse equation and the value of
[
9
R. T. Schermer and J. H. Cole, “Improved bend loss formula verified for Optical Fiber by simulation and experiment,” IEEE J. Quantum Electron.
43(10), 899–909 (2007). [CrossRef]
]:
where
and
.
Figure 10(b) shows once again an excellent agreement between numerical and semi-analytical calculations.
6. Conclusion
We have proposed a semi-analytical analysis of the multifilament core fiber using the equivalent step index technique based on fundamental space mode (ESI-FSM) technique. The fundamental space filling mode was computed on a circular cell providing an equivalent step index model for the MFC much more accurate than modeling based on the moment theory. It allows to predict the fiber properties using semi-analytical expressions. Using this technique we were able to determine the useful design domains of the MFC from the effective normalized frequency value. The low bending loss and large effective area of the fundamental mode were also recovered.
The main interest of the MFC fiber structure is its emulation of a low NA step index fiber using large NA doped filaments for eyesafe wavelength amplification. This theory will be useful to optimize the MFC fiber design.