Abstract

Based on the overlap integral of electromagnetic fields in neighboring cores, a calculating method is proposed for obtaining the coupling coefficient between two adjacent trench-assisted non-identical cores. And a kind of heterogeneous trench-assisted multi-core fiber (Hetero-TA-MCF) with 12 cores is proposed to achieve large effective area (Aeff) and high density of cores. As bending radius becomes larger than 50 mm, the crosstalk value at 1550-nm wavelength of the Hetero-TA-MCF is about −42 dB after 100-km propagation and the Aeff of this Hetero-TA-MCF can reach 100 µm2.

©2012 Optical Society of America

1. Introduction

How to cope with the exponentially increasing demand for transmission capacity per fiber is a hot topic nowadays. As an approach to achieve space division multiplexing (SDM), multi-core fiber (MCF) has been proposed to solve the issue related above [1].

Recently, several kinds of homogeneous MCFs (Homo-MCFs) in which all the cores are identical to each other have been designed and fabricated in order to realize long-haul transmission with low crosstalk [2, 3]. Furthermore, a type of optical fiber called heterogeneous MCF (Hetero-MCF) has been proposed to obtain much lower crosstalk, in which there are not only identical cores but also non-identical cores and the cores are more closely packed in a definite space [4]. On the other hand, a trench-assisted MCF (TA-MCF) that realizes much smaller crosstalk and larger effective area (Aeff) comparing to MCF with step-index profile also has been proposed [5]. It has been proved that the crosstalk between the identical cores will become larger and larger as bending radius (R) increases [6]. However the Hetero-MCF is insensitive to the bending radius after the R reaching a threshold value which was called Rpk in [7]. Moreover, if the cores have slight differences in their core radii and core refractive indices, the maximum power transferred between the cores goes down drastically [4]. Therefore, Hetero-MCF is a good candidate for the research of fiber under the bending condition. In addition, if we want to accommodate more cores inside the fiber, the core pitch between each pair of cores needs to be reduced. But small core pitch will result in a large crosstalk between cores. So in order to lower the crosstalk and meanwhile increase the core number in the fiber, a Hetero-MCF which has an index trench structure around each core (Hetero-TA-MCF) can be a solution.

In this paper, we propose a Hetero-TA-MCF with high density of cores to realize a large-scale SDM transmission. Besides increasing the number of cores, we also aim at achieving low crosstalk and enlarging the Aeff as much as possible. Additionally, when we analyze the crosstalk characteristics between a couple of non-identical cores with trench regions, we should know the coupling coefficient (κ) between them. However, to the best of our knowledge, an analytical method which can be used to figure out the value of κ between two trench-assisted non-identical cores has not been reported yet. Therefore, before proposing a new type of Hetero-TA-MCF and analyzing the characteristics of it, we will introduce an effective method for calculating the κ between two trench-assisted non-identical cores by considering the overlap integral of electromagnetic fields in the neighboring cores. Then we will compare such method with finite element method (FEM) [8] by calculating the crosstalk of Homo-TA-MCF in order to estimate the accuracy degree of it. Based on this calculating method of κ, an optimal design condition for the Hetero-TA-MCF will be given.

2. Calculation of κ between trench-assisted non-identical cores

2.1 Analysis of refractive-index distribution in the coupled region

Figure 1 shows the refractive index profile of two trench-assisted non-identical cores and the part outside the cores. The refractive-index distribution in the entire coupled region can be written as

N2(r,θ)=N12(r,θ)+N22(r,θ)n2(r,θ),
where N1(r, θ) and N2(r, θ) represent the refractive-index distribution of each core with trench structure, and n(r, θ) means the refractive-index distribution outside the cores, which are shown in Fig. 1(a), Fig. 1(b), and Fig. 1(c), respectively. The expression of coupling coefficient between two waveguides is expressed as [9]
κpq=ωε0++(N2Nq2)Ep*Eqdxdy++uz(Ep*×Hp+Ep×Hp*)dxdy,
where ω is an angular frequency of the sinusoidally varying electromagnetic fields, ε0 is the permittivity of the medium, and uz means the outward-directed unit vector. The pair of p and q is either (1, 2) or (2, 1). E and H represent the electric and magnetic fields respectively.

 figure: Fig. 1

Fig. 1 The profile of refractive index in two trench-assisted non-identical cores and the part outside the cores.

Download Full Size | PPT Slide | PDF

As shown in Fig. 2(a) , N2−N22 is zero except the region inside core 1, so the difference of the refractive-index distribution inside the core 1 is n12−ncl2, while according to Fig. 2(b), N2−N12 is zero except the region inside core 2, so the difference of the refractive-index distribution inside the core 2 is n22−ncl2. Furthermore, the denominator of κpq equals 4P, where P means the total power flow [9]. So the expression of κpq can be rewritten as

κpq=ωε04P02π0a1p(np2ncl2)Ep*Eqrdrdθ,
where Ep and Eq represent the amplitude of electric field distribution of core p inside the range of core p, and the amplitude of electric field distribution of core q inside the range of core p, respectively. a1-p denotes the radius of core p.

 figure: Fig. 2

Fig. 2 Difference of the refractive-index distributions. (a) N2N22. (b) N2N12.

Download Full Size | PPT Slide | PDF

2.2 Derivation of the expression of κ between trench-assisted non-identical cores

The electric fields in optical fibers are expressed in cylindrical coordinates as

E˜=E(r,θ)ej(ωtβz).
Substituting Eq. (4) into Maxwell’s equation, we can obtain the wave equation as
2Ezr2+1rEzr+1r22Ezθ2+[k2n(r,θ)2β2]Ez=0.
Here, we define the wave number in the core m (m = 1, 2), cladding and trench along the transversal direction as follows:
ξm=nm2k2βm2,
σm=βm2ncl2k2,
γm=βm2ntr2k2,
where β is the propagation constant and k is the wavenumber in a vacuum. The normalized frequency (V1-m, V2-m), the normalized transverse wave number in core m (U1-m), that in cladding (W1-m) and that in trench (W2-m) can be expressed as follows:
V1m=a1mknm2ncl2,
V2m=a1mkncl2ntr2,
W1m=a1mσm=1.1428V1m0.996,
U1m=a1mξm=V1m2W1m2,
W2m=a1mγm=V2m2+W1m2,
where a1-m is the radius of core m and Eq. (11) is the approximation which is in error by less than 0.2 percent for 1.5 ≤ V1-m ≤ 2.5 [10].

In the core region of fiber, the solutions for Eq. (5) of TM modes are the 0th-order Bessel function J0(ξmrm) and the 0th-order Neumann function N0(ξmrm) [9], respectively. However, N0(ξmrm) diverges infinitely at r = 0. Therefore J0(ξmrm) is the proper solution for the field in the core. In the cladding region of fiber, the solutions for Eq. (5) of TM modes are the modified Bessel function of the first kind I0(σmrm) and modified Bessel functions of the second kind K0(σmrm), respectively. However, I0(σmrm) diverges infinitely at r = ∞. Therefore K0(σmrm) is the proper solution for the field in the cladding. In the same principle, K0(σmrm) is the proper solution for the field in the trench. Nevertheless, in hybrid modes, the solutions for Eq. (5) are given by the product of the nth-order Bessel functions and cos( + ψ). Thus, the z-components of the electric field can be obtained as

Emz={AmJn(ξmrm)cos(nθ+ψ)(inCom)BmKn(σmrm)cos(nθ+ψ)(inICm)CmKn(γmrm)cos(nθ+ψ)(inTrm)DmKn(σmrm)cos(nθ+ψ)(inOC)EmKn(γmrm)cos(nθ+ψ)(inTrm')FmKn(σmrm)cos(nθ+ψ)(inICm'),
the azimuthal dependency of the electric fields in axially symmetric fibers is expressed by cos ( + ψ), where n is an integer and ψ denotes the phase. As shown in Fig. 3 , core m, inner cladding between core m and trench m, trench m, outer cladding outside trench m, trench m’ and inner cladding inside trench m’ are abbreviated as Co m, IC m, Tr m, OC, Tr m’, and IC m’. D is the core pitch, a1-m is the radius of core m, a2-m is the distance from the center of core m to the inner circumference of trench m, and a3-m is the distance from the center of core m to the outer circumference of trench m. R, R1 and R2 denote the distance between the center of core m to the objective point in IC m’, the distance from the center of core m to the outer circumference of trench m’ and the distance from the center of core m to the inner circumference of trench m’, respectively.

 figure: Fig. 3

Fig. 3 The profile of core m with trench structures.

Download Full Size | PPT Slide | PDF

Without any doubt, the electric field should be continuous in each boundary, therefore the solutions in these six parts have the relationship which is shown as follows:

{AmJn(U1m)=BmKn(W1m)BmKn(W1ma2ma1m)=CmKn(W2ma2ma1m)CmKn(W2ma3ma1m)=DmKn(W3ma3ma1m)DmKn(W3mR1a1m)=EmKn(W4mR1a1m)EmKn(W4mR2a1m)=FmKn(W5mR2a1m).
Based on the boundary condition related above, the expression of Dm and Fm can be obtained, which are shown as follows:
Dm=LmAm,
Fm=QmAm,
where
Lm=Jn(U1m)Kn(W1ma2ma1m)Kn(W2ma3ma1m)Kn(W1m)Kn(W2ma2ma1m)Kn(W3ma3ma1m),
Qm=LmKn(W3mR1a1m)Kn(W4mR2a1m)Kn(W4mR1a1m)Kn(W5mR2a1m).
And based on [9], the amplitude coefficient Am of the field is given by
Am=U1mW1mβma1m2V1mJ1(U1m)2Pπε0n1mc,
where c is the velocity of light in a vacuum. Figure 4 illustrates the geometries for the calculation of the coupling coefficient. Setting n = 1 and using the equation of electric fields in [9] by assuming s = s1 = s2 = −1, the electronic fields of the fundamental HE11 mode inside the core p can be express as follows:
Ep={Epx=jApβpa1pU1pJ0(U1pra1p)cosψEpy=jApβpa1pU1pJ0(U1pra1p)sinψEpz=ApJ1(U1pra1p)cos(θ+ψ),
Eq={Eqx=jQqAqβqa1qW1qK0(W1qRa1q)cosψEqy=jQqAqβqa1qW1qK0(W1qRa1q)sinψEqz=QqAqK1(W1qRa1q)cos(Θ+ψ).
where
R=D2+r22DrcosθDrcosθ,
a3p=D2+R122DR1cos(πΘ)DR1cos(πΘ),
r=D2+R22DRcos(πΘ)DRcos(πΘ),
R1(Da3p)(Drcosθ)Dr.
In the same principle,

 figure: Fig. 4

Fig. 4 Geometries for the calculation of the coupling coefficient.

Download Full Size | PPT Slide | PDF

R2(Da2p)(Drcosθ)Dr.

Inside Eq. (3), Ep*Eq can be expressed as

Ep*Eq=QqApAqβpβqa1pa1qU1pW1qJ0(U1pra1p)K0(W1qRa1q)+QqApAqJ1(U1pra1p)K1(W1qRa1q)cos(θ+ψ)cos(Θ+ψ).
Since the second term of the right-hand side of the equation above is sufficiently smaller than the first term, the integration of the first part in the square brackets of Eq. (3) becomes:

S1=02π0a1p(np2ncl2)Ep*Eqrdrdθ=(np2ncl2)LqApAqβpβqa1pa1qU1pW1q×02π0a1pJ0(U1pra1p)K0(W1qRa1q)K1(P1DrcosθDr)K1(Y1DrcosθDr)K1(P2DrcosθDr)K1(Y2DrcosθDr)rdrdθ,

whereP1=W1-q(Da3-p)/a1-q,P2=W2-q(Da3-p)/a1-q,Y1=W2-q(Da2-p)/a1-q,Y2=W1-q(Da2-p)/a1-q. When the argument of the modified Bessel function Kn(z) in Eq. (29) is large, it can be approximated as

Kn(z)π2zexp(z).

Substitution of Eq. (30) into Eq. (29) gives

S1=(np2ncl2)LqApAqβpβqa1pa1qU1pW1qπa1q2W1qDexp(W1qDa1q)×0a1pJ0(U1pra1p)exp[(P2P1+Y2Y1)DDr]rdr02πexp[(W1qa1qP2P1+Y2Y1Dr)rcosθ]dθ.
By using the integral formulas of the Bessel functions [11]:
I0(z)=1π0πexp(zcosθ)dθ,
the coupling coefficient κpq can be expressed as

κpq=k(np2ncl2)W1pU1qLqπa1q2W1qDexp(W1qDa1q)npnqa1pa1qV1pV1qJ1(U1p)J1(U1q)×0a1pJ0(U1pra1p)I0[(W1qa1qP2P1+Y2Y1Dr)r]exp[(P2P1+Y2Y1)DDr]rdr.

2.3 Comparison with finite element method (FEM)

In order to estimate the accuracy degree of this method, we compare it with FEM by calculating the coupling coefficient between adjacent cores in Homo-TA-MCF. Here we use the Homo-TA-7-core fiber as a model to do this comparison. Figure 5 shows index profile of a core with trench and the schematic of the Homo-TA-7-core model. The coiling diameter of the fiber was assumed to be 210 mm, which is the same with the value in [5]. The parameters which we used are summarized in Table 1 .

 figure: Fig. 5

Fig. 5 Schematic of a core with index trench and Homo-TA-7-core model.

Download Full Size | PPT Slide | PDF

Tables Icon

Table 1. Structural Parameters for Calculation

Figure 6 illustrates the length dependence of simulated crosstalk (XT) of step-index MCF (S-MCF) and trench-assisted MCF (TA-MCF) at 1550-nm wavelength. The blue solid line relates the simulation result of S-MCF. On the other hand, the red solid line represents the result of TA-MCF which is calculated by using the above-mentioned method and the green solid line represents the result of TA-MCF that is obtained by using the FEM [8]. The error between crosstalk values which were calculated by these two methods is about −0.4 dB, a sufficiently small value, which proves the feasibility of this analytical method.

 figure: Fig. 6

Fig. 6 Simulated crosstalk at 1550-nm wavelength as function of length.

Download Full Size | PPT Slide | PDF

3. Model of Hetero-TA-MCF with high density of cores

3.1 Design of Hetero-TA-MCF

Figure 7 shows the schematic of Hetero-TA-12-core model. The index profile of the core with trench of this structure is the same with the one shown in Fig. 5. In the simulation of this work, we assumed two sorts of cores which were represented by two different colors — orange and green separately. The calculated wavelength was set to be 1550 nm, while relative refractive-index difference between ncl and ntr2) was assumed to be −0.7%, which is the structural parameter used in the fabricated TA-MCFs [5]. And we define the distance between the outer circumferences of adjacent trenches (Dtr) as a value that is not smaller than 2 µm for making sure the trench not overlap to each other.

 figure: Fig. 7

Fig. 7 Schematic of Hetero-TA-12-core model.

Download Full Size | PPT Slide | PDF

The reason why we design such Hetero-TA-MCF with a ring layout is that the cutoff wavelength of each core will not be very long in this case. As a result, we can ensure a wide wavelength band which is used for the single-mode transmission. The cutoff wavelength related here is defined as the wavelength at which the bending loss of LP11 mode equals to 1 dB/m at the bending radius of 140 mm. If we arrange several layers inside the fiber, the cutoff wavelength of the cores at the center part will have a longer cutoff wavelength than that of the cores in the outer layer. This is because the tight confinement of the LP11 mode of the cores at the center will result in a long cutoff wavelength when index trench structures are deployed [5].

In addition, if we expect the Aeff to reach 110 µm2 or 80 µm2, the outer cladding thickness (OCT) needs to be at least 40 µm [12] or 30 µm [13], respectively, for reducing the micro-bending loss. And if we want to decrease the failure probability of a fiber in order to guarantee the mechanical reliability, the cladding diameter (CD) should not be larger than 200 µm [14]. So based on the linear relationship of Aeff and OCT, the OCT should be at least 37 µm when the Aeff equals 100 µm2. Thus, if we set the maximum of the CD to be around 200 µm, the limit value of Λ should be 33 or so when the Aeff equals 100 µm2. In this case, the core number of 12 here is the limit value for this ring layout due to the required value of CD.

3.2 Crosstalk characteristics

In Hetero-MCFs, there is a threshold value of R which was proposed as Rpk [7]. The crosstalk is degraded at R < Rpk due to the phase-matching between non-identical cores. In this phase-matching region, the bend perturbations are crucial. In the non-phase-matching region of R > Rpk, on the other hand, the crosstalk is dominated by the statistical properties [15]. Therefore, a large effective index difference (Δneff) between cores will be required for pushing the value of Rpk toward sufficiently small range. Figure 8 shows the required Δneff value between two non-identical cores as function of Λ and Rpk. If we try to shift Rpk to R of smaller than 5 cm with Λ being about 33 µm, Δneff should be around 0.001 or larger than it between non-identical cores.

 figure: Fig. 8

Fig. 8 Required Δneff as function of the Λ and Rpk.

Download Full Size | PPT Slide | PDF

Figure 9 illustrates the effective index value (neff) of the fundamental mode at 1550-nm wavelength as function of core radius and core Δ1 in four conditions — (a) r2/r1 = 2.0, W/r1 = 1.0, (b) r2/r1 = 2.0, W/r1 = 1.1, (c) r2/r1 = 2.0, W/r1 = 1.2, (d) r2/r1 = 2.0, W/r1 = 1.3, where neff was simulated by FEM. The black solid lines and the black dashed lines represent the values of neff and Aeff, respectively. The couple of white solid lines correspond to the upper limit of bending loss (BL) of the higher-order mode (HOM) at 1530 nm when R equals 140 mm and the lower limit of bending loss (BL) of the fundamental mode (FM) at 1625 nm when R equals 30 mm. To define the single-mode operation, the bending loss of LP11-like HOM should be > 1 dB/m at R = 140 mm [14] and we assume the limit value of the bending loss of FM to be 0.5 dB/100 turns at R = 30 mm, which is described in ITU-T recommendations G.655 and G.656. Therefore, in order to guarantee the transmission with single-mode operation and low bending loss from C-band to L-band, we only research the field that is surrounded by the couple of white solid lines.

 figure: Fig. 9

Fig. 9 Effective index value of the fundamental mode as function of core radius and core Δ1, where (a) r2/r1 = 2.0, W/r1 = 1.0, (b) r2/r1 = 2.0, W/r1 = 1.1, (c) r2/r1 = 2.0, W/r1 = 1.2, and (d) r2/r1 = 2.0, W/r1 = 1.3.

Download Full Size | PPT Slide | PDF

Here, we fix r2/r1 to be 2.0 and change W/r1 from 1.0 to 1.3 to investigate how the width of trench region influences the crosstalk between the neighboring cores. In the conditions (a), (b), (c), and (d), we do not set the target value of Aeff to be 110 µm2, because the core radii of cores are too large to ensure the required value of CD. So we select four pairs of cores with Aeff equaling 100 µm2 for the simulation in conditions (a), (b), (c), and (d). In Fig. 9(a), one of the cores has a radius of 5.34 µm and Δ1 = 0.304% and another core has a radius of 4.86 µm and Δ1 = 0.242%. In Fig. 9(b), one of the cores has a radius of 5.27 µm and Δ1 = 0.293% and another core has a radius of 4.83 µm and Δ1 = 0.239%. In Fig. 9(c), one of the cores has a radius of 5.18 µm and Δ1 = 0.280% and another core has a radius of 4.77 µm and Δ1 = 0.234%. In Fig. 9(d), one of the cores has a radius of 5.03 µm and Δ1 = 0.263% and another core has a radius of 4.75 µm and Δ1 = 0.240%. Furthermore, in order to make sure the trench not overlap to each other, we defined Dtr to be not smaller than 2 µm. Under this requirement, we set each required Λ. The optical properties of the cores in each condition are summarized in Table 2 .

Tables Icon

Table 2. Optical Properties of the Cores in Different Conditions (1550 nm)

In Table 2, the relative value of core multiplicity factor (RCMF) was used to compare the core density of fibers [12]. And the core multiplicity factor (CMF) is defined as

CMF=NcoreAeffπ(CD/2)2,
where Ncore is a number of core and CD means a cladding diameter. And RCMF is a ratio between CMF of a MCF and a standard single core single mode fiber with Aeff = 80 µm2 at 1.55 µm and CD = 125 µm.

According to the equation of κ which was proposed above, the crosstalk can be calculated by using the novel power coupling coefficient (PCC) definition [15] and the coupled-power theory [3]. The PCC was proposed with correlation length (dc) and dc of 0.05-m was proved to agree well with the measurement in [6]. Therefore, we assume dc to be 0.05 m in this simulation of crosstalk, and the twisting rate is defined to be 5 turns per 100 m.

Figure 10 shows the crosstalk of the Hetero-TA-12-core fiber at 100-km propagation as function of bending radius under the four kinds of conditions mentioned above. We can find that the crosstalk becomes smaller as the width of trench structure increases. But meanwhile the Δneff between the adjacent cores decreases and CD should be larger. Therefore, if we want the Rpk of the designed fiber smaller than 50 mm, the Δneff should be around 0.001. In this case, r2/r1 = 2.0 and W/r1 = 1.0 is the optimal design condition for the Hetero-TA-MCF and under this condition, as R becomes larger than 50 mm, the worst-case crosstalk is around −42 dB after 100-km propagation and it is insensitive to the bending radius.

 figure: Fig. 10

Fig. 10 Crosstalk of Hetero-TA-12-core fiber at 100-km propagation as function of bending radius.

Download Full Size | PPT Slide | PDF

Figure 11 illustrates the crosstalk of the Homo-TA-12-core fiber at 100-km propagation as function of bending radius under same four conditions. The Homo-TA-12-core fiber mentioned here has the same ring structure with the Hetero-TA-12-core fiber and we choose four sorts of cores for the Homo-TA-12-core fiber under these four conditions. In order to compare the crosstalk characteristics of the Homo-TA-12-core fiber with that of the Hetero-TA-12-core fiber, we assume these four kinds of cores to have the same core parameters with the first kind of core in each condition of the Hetero-TA-12-core fiber that we described above. For the condition (a), r1 = 5.34 µm, Δ1 = 0.304%, and Λ = 33 µm. For condition (b), r1 = 5.27 µm, Δ1 = 0.293%, and Λ = 34 µm. For condition (c), r1 = 5.18 µm, Δ1 = 0.280%, and Λ = 34 µm. And for condition (d), r1 = 5.03 µm, Δ1 = 0.263%, and Λ = 35 µm. We can find obviously that the crosstalk of the Homo-TA-MCF become larger and larger as increasing the bending radius. Therefore, we can see the merit of Hetero-TA-MCF clearly from this comparison.

 figure: Fig. 11

Fig. 11 Crosstalk of Homo-TA-12-core fiber at 100-km propagation as function of bending radius.

Download Full Size | PPT Slide | PDF

4. Conclusion

We proposed an analytical method to get the κ value between two adjacent trench-assisted non-identical cores. And we designed a sort of Hetero-TA-12-core fiber with Aeff of 100-µm2 and showed that RCMF can reach 5.8. Based on the calculating method of κ, crosstalk characteristics of the Hetero-TA-MCF were analyzed. The parameters that Λ = 33 µm, Aeff = 100 µm2, r2/r1 = 2.0, and W/r1 = 1.0 were proved to be the best design condition for this model. As R becomes larger than 50 mm, the worst-case crosstalk at 1550-nm wavelength is about −42 dB after 100-km propagation, which also demonstrates that the Hetero-TA-MCF with high density of cores is insensitive to the bending radius.

Acknowledgments

This work was partially supported by the National Institute of Information and Communication Technology (NICT), Japan under “Research on Innovative Optical Fiber Technology”.

References and links

1. T. Morioka, “New generation optical infrastructure technologies: “EXACT initiative” towards 2020 and beyond,” in Proceedings of 14th OptoElectronics and Communications Conference (Institute of Electrical and Electronics Engineers, 2009), paper FT4.

2. K. Takenaga, S. Tanigawa, N. Guan, S. Matsuo, K. Saitoh, and M. Koshiba, “Reduction of crosstalk by quasi-homogeneous solid multi-core fiber,” in Optical Fiber Communication Conference, OSA Technical Digest (CD) (Optical Society of America, 2010), paper OWK7.

3. K. Takenaga, Y. Arakawa, S. Tanigawa, N. Guan, S. Matsuo, K. Saitoh, and M. Koshiba, “An investigation on crosstalk in multi-core fibers by introducing random fluctuation along longitudinal direction,” IEICE Trans. Commun. E 94-B, 409–416 (2011).

4. M. Koshiba, K. Saitoh, and Y. Kokubun, “Heterogeneous multi-core fibers: proposal and design principle,” IEICE Electron. Express 6(2), 98–103 (2009). [CrossRef]  

5. K. Takenaga, Y. Arakawa, S. Tanigawa, N. Guan, S. Matsuo, K. Saitoh, and M. Koshiba, “Reduction of crosstalk by trench-assisted multi-core fiber,” in Optical Fiber Communication Conference, OSA Technical Digest (CD) (Optical Society of America, 2011), paper OWJ4.

6. S. Matsuo, K. Takenaga, Y. Arakawa, Y. Sasaki, S. Tanigawa, K. Saitoh, and M. Koshiba, “Crosstalk behavior of cores in multi-core fiber under bent condition,” IEICE Electron. Express 8(6), 385–390 (2011). [CrossRef]  

7. T. Hayashi, T. Nagashima, O. Shimakawa, T. Sasaki, and E. Sasaoka, “Crosstalk variation of multi-core fiber due to fiber bend,” in Proceedings of 36th European Conference and Exhibition on Optical Communication (Institute of Electrical and Electronics Engineers, 2010), paper We.8.F.6.

8. K. Saitoh and M. Koshiba, “Full-vectorial imaginary-distance beam propagation method based on a finite element scheme: Application to photonic crystal fibers,” IEEE J. Quantum Electron. 38(7), 927–933 (2002). [CrossRef]  

9. K. Okamoto, Fundamentals of Optical Waveguides (Corona Publishing, 1992), Chap. 4.

10. H. D. Rudolph and E. G. Neuman, “Approximations for the eigenvalues of the fundamental mode of a step index glass fiber waveguide,” Nachrichtentech. Elektron. 29, 328–329 (1976).

11. A. W. Snyder and J. D. Love, Optical Waveguide Theory (Chapman & Hall, 1983), Chap. 37.

12. K. Takenaga, Y. Arakawa, Y. Sasaki, S. Tanigawa, S. Matsuo, K. Saitoh, and M. Koshiba, “A large effective area multi-core fiber with an optimized cladding thickness,” Opt. Express 19(26), B543–B550 (2011). [CrossRef]   [PubMed]  

13. S. Matsuo, K. Takenaga, Y. Arakawa, Y. Sasaki, S. Taniagwa, K. Saitoh, and M. Koshiba, “Large-effective-area ten-core fiber with cladding diameter of about 200 μm,” Opt. Lett. 36(23), 4626–4628 (2011). [CrossRef]   [PubMed]  

14. T. Matsui, K. Nakajima, and C. Fukai, “Applicability of photonic crystal fiber with uniform air-hole structure to high-speed and wide-band transmission over conventional telecommunication bands,” J. Lightwave Technol. 27(23), 5410–5416 (2009). [CrossRef]  

15. M. Koshiba, K. Saitoh, K. Takenaga, and S. Matsuo, “Multi-core fiber design and analysis: coupled-mode theory and coupled-power theory,” Opt. Express 19(26), B102–B111 (2011). [CrossRef]   [PubMed]  

References

  • View by:
  • |
  • |
  • |

  1. T. Morioka, “New generation optical infrastructure technologies: “EXACT initiative” towards 2020 and beyond,” in Proceedings of 14th OptoElectronics and Communications Conference (Institute of Electrical and Electronics Engineers, 2009), paper FT4.
  2. K. Takenaga, S. Tanigawa, N. Guan, S. Matsuo, K. Saitoh, and M. Koshiba, “Reduction of crosstalk by quasi-homogeneous solid multi-core fiber,” in Optical Fiber Communication Conference, OSA Technical Digest (CD) (Optical Society of America, 2010), paper OWK7.
  3. K. Takenaga, Y. Arakawa, S. Tanigawa, N. Guan, S. Matsuo, K. Saitoh, and M. Koshiba, “An investigation on crosstalk in multi-core fibers by introducing random fluctuation along longitudinal direction,” IEICE Trans. Commun. E 94-B, 409–416 (2011).
  4. M. Koshiba, K. Saitoh, and Y. Kokubun, “Heterogeneous multi-core fibers: proposal and design principle,” IEICE Electron. Express 6(2), 98–103 (2009).
    [Crossref]
  5. K. Takenaga, Y. Arakawa, S. Tanigawa, N. Guan, S. Matsuo, K. Saitoh, and M. Koshiba, “Reduction of crosstalk by trench-assisted multi-core fiber,” in Optical Fiber Communication Conference, OSA Technical Digest (CD) (Optical Society of America, 2011), paper OWJ4.
  6. S. Matsuo, K. Takenaga, Y. Arakawa, Y. Sasaki, S. Tanigawa, K. Saitoh, and M. Koshiba, “Crosstalk behavior of cores in multi-core fiber under bent condition,” IEICE Electron. Express 8(6), 385–390 (2011).
    [Crossref]
  7. T. Hayashi, T. Nagashima, O. Shimakawa, T. Sasaki, and E. Sasaoka, “Crosstalk variation of multi-core fiber due to fiber bend,” in Proceedings of 36th European Conference and Exhibition on Optical Communication (Institute of Electrical and Electronics Engineers, 2010), paper We.8.F.6.
  8. K. Saitoh and M. Koshiba, “Full-vectorial imaginary-distance beam propagation method based on a finite element scheme: Application to photonic crystal fibers,” IEEE J. Quantum Electron. 38(7), 927–933 (2002).
    [Crossref]
  9. K. Okamoto, Fundamentals of Optical Waveguides (Corona Publishing, 1992), Chap. 4.
  10. H. D. Rudolph and E. G. Neuman, “Approximations for the eigenvalues of the fundamental mode of a step index glass fiber waveguide,” Nachrichtentech. Elektron. 29, 328–329 (1976).
  11. A. W. Snyder and J. D. Love, Optical Waveguide Theory (Chapman & Hall, 1983), Chap. 37.
  12. K. Takenaga, Y. Arakawa, Y. Sasaki, S. Tanigawa, S. Matsuo, K. Saitoh, and M. Koshiba, “A large effective area multi-core fiber with an optimized cladding thickness,” Opt. Express 19(26), B543–B550 (2011).
    [Crossref] [PubMed]
  13. S. Matsuo, K. Takenaga, Y. Arakawa, Y. Sasaki, S. Taniagwa, K. Saitoh, and M. Koshiba, “Large-effective-area ten-core fiber with cladding diameter of about 200 μm,” Opt. Lett. 36(23), 4626–4628 (2011).
    [Crossref] [PubMed]
  14. T. Matsui, K. Nakajima, and C. Fukai, “Applicability of photonic crystal fiber with uniform air-hole structure to high-speed and wide-band transmission over conventional telecommunication bands,” J. Lightwave Technol. 27(23), 5410–5416 (2009).
    [Crossref]
  15. M. Koshiba, K. Saitoh, K. Takenaga, and S. Matsuo, “Multi-core fiber design and analysis: coupled-mode theory and coupled-power theory,” Opt. Express 19(26), B102–B111 (2011).
    [Crossref] [PubMed]

2011 (5)

K. Takenaga, Y. Arakawa, S. Tanigawa, N. Guan, S. Matsuo, K. Saitoh, and M. Koshiba, “An investigation on crosstalk in multi-core fibers by introducing random fluctuation along longitudinal direction,” IEICE Trans. Commun. E 94-B, 409–416 (2011).

S. Matsuo, K. Takenaga, Y. Arakawa, Y. Sasaki, S. Tanigawa, K. Saitoh, and M. Koshiba, “Crosstalk behavior of cores in multi-core fiber under bent condition,” IEICE Electron. Express 8(6), 385–390 (2011).
[Crossref]

K. Takenaga, Y. Arakawa, Y. Sasaki, S. Tanigawa, S. Matsuo, K. Saitoh, and M. Koshiba, “A large effective area multi-core fiber with an optimized cladding thickness,” Opt. Express 19(26), B543–B550 (2011).
[Crossref] [PubMed]

S. Matsuo, K. Takenaga, Y. Arakawa, Y. Sasaki, S. Taniagwa, K. Saitoh, and M. Koshiba, “Large-effective-area ten-core fiber with cladding diameter of about 200 μm,” Opt. Lett. 36(23), 4626–4628 (2011).
[Crossref] [PubMed]

M. Koshiba, K. Saitoh, K. Takenaga, and S. Matsuo, “Multi-core fiber design and analysis: coupled-mode theory and coupled-power theory,” Opt. Express 19(26), B102–B111 (2011).
[Crossref] [PubMed]

2009 (2)

2002 (1)

K. Saitoh and M. Koshiba, “Full-vectorial imaginary-distance beam propagation method based on a finite element scheme: Application to photonic crystal fibers,” IEEE J. Quantum Electron. 38(7), 927–933 (2002).
[Crossref]

1976 (1)

H. D. Rudolph and E. G. Neuman, “Approximations for the eigenvalues of the fundamental mode of a step index glass fiber waveguide,” Nachrichtentech. Elektron. 29, 328–329 (1976).

Arakawa, Y.

K. Takenaga, Y. Arakawa, S. Tanigawa, N. Guan, S. Matsuo, K. Saitoh, and M. Koshiba, “An investigation on crosstalk in multi-core fibers by introducing random fluctuation along longitudinal direction,” IEICE Trans. Commun. E 94-B, 409–416 (2011).

S. Matsuo, K. Takenaga, Y. Arakawa, Y. Sasaki, S. Tanigawa, K. Saitoh, and M. Koshiba, “Crosstalk behavior of cores in multi-core fiber under bent condition,” IEICE Electron. Express 8(6), 385–390 (2011).
[Crossref]

K. Takenaga, Y. Arakawa, Y. Sasaki, S. Tanigawa, S. Matsuo, K. Saitoh, and M. Koshiba, “A large effective area multi-core fiber with an optimized cladding thickness,” Opt. Express 19(26), B543–B550 (2011).
[Crossref] [PubMed]

S. Matsuo, K. Takenaga, Y. Arakawa, Y. Sasaki, S. Taniagwa, K. Saitoh, and M. Koshiba, “Large-effective-area ten-core fiber with cladding diameter of about 200 μm,” Opt. Lett. 36(23), 4626–4628 (2011).
[Crossref] [PubMed]

Fukai, C.

Guan, N.

K. Takenaga, Y. Arakawa, S. Tanigawa, N. Guan, S. Matsuo, K. Saitoh, and M. Koshiba, “An investigation on crosstalk in multi-core fibers by introducing random fluctuation along longitudinal direction,” IEICE Trans. Commun. E 94-B, 409–416 (2011).

Kokubun, Y.

M. Koshiba, K. Saitoh, and Y. Kokubun, “Heterogeneous multi-core fibers: proposal and design principle,” IEICE Electron. Express 6(2), 98–103 (2009).
[Crossref]

Koshiba, M.

K. Takenaga, Y. Arakawa, S. Tanigawa, N. Guan, S. Matsuo, K. Saitoh, and M. Koshiba, “An investigation on crosstalk in multi-core fibers by introducing random fluctuation along longitudinal direction,” IEICE Trans. Commun. E 94-B, 409–416 (2011).

S. Matsuo, K. Takenaga, Y. Arakawa, Y. Sasaki, S. Tanigawa, K. Saitoh, and M. Koshiba, “Crosstalk behavior of cores in multi-core fiber under bent condition,” IEICE Electron. Express 8(6), 385–390 (2011).
[Crossref]

M. Koshiba, K. Saitoh, K. Takenaga, and S. Matsuo, “Multi-core fiber design and analysis: coupled-mode theory and coupled-power theory,” Opt. Express 19(26), B102–B111 (2011).
[Crossref] [PubMed]

S. Matsuo, K. Takenaga, Y. Arakawa, Y. Sasaki, S. Taniagwa, K. Saitoh, and M. Koshiba, “Large-effective-area ten-core fiber with cladding diameter of about 200 μm,” Opt. Lett. 36(23), 4626–4628 (2011).
[Crossref] [PubMed]

K. Takenaga, Y. Arakawa, Y. Sasaki, S. Tanigawa, S. Matsuo, K. Saitoh, and M. Koshiba, “A large effective area multi-core fiber with an optimized cladding thickness,” Opt. Express 19(26), B543–B550 (2011).
[Crossref] [PubMed]

M. Koshiba, K. Saitoh, and Y. Kokubun, “Heterogeneous multi-core fibers: proposal and design principle,” IEICE Electron. Express 6(2), 98–103 (2009).
[Crossref]

K. Saitoh and M. Koshiba, “Full-vectorial imaginary-distance beam propagation method based on a finite element scheme: Application to photonic crystal fibers,” IEEE J. Quantum Electron. 38(7), 927–933 (2002).
[Crossref]

Matsui, T.

Matsuo, S.

S. Matsuo, K. Takenaga, Y. Arakawa, Y. Sasaki, S. Taniagwa, K. Saitoh, and M. Koshiba, “Large-effective-area ten-core fiber with cladding diameter of about 200 μm,” Opt. Lett. 36(23), 4626–4628 (2011).
[Crossref] [PubMed]

K. Takenaga, Y. Arakawa, Y. Sasaki, S. Tanigawa, S. Matsuo, K. Saitoh, and M. Koshiba, “A large effective area multi-core fiber with an optimized cladding thickness,” Opt. Express 19(26), B543–B550 (2011).
[Crossref] [PubMed]

S. Matsuo, K. Takenaga, Y. Arakawa, Y. Sasaki, S. Tanigawa, K. Saitoh, and M. Koshiba, “Crosstalk behavior of cores in multi-core fiber under bent condition,” IEICE Electron. Express 8(6), 385–390 (2011).
[Crossref]

K. Takenaga, Y. Arakawa, S. Tanigawa, N. Guan, S. Matsuo, K. Saitoh, and M. Koshiba, “An investigation on crosstalk in multi-core fibers by introducing random fluctuation along longitudinal direction,” IEICE Trans. Commun. E 94-B, 409–416 (2011).

M. Koshiba, K. Saitoh, K. Takenaga, and S. Matsuo, “Multi-core fiber design and analysis: coupled-mode theory and coupled-power theory,” Opt. Express 19(26), B102–B111 (2011).
[Crossref] [PubMed]

Nakajima, K.

Neuman, E. G.

H. D. Rudolph and E. G. Neuman, “Approximations for the eigenvalues of the fundamental mode of a step index glass fiber waveguide,” Nachrichtentech. Elektron. 29, 328–329 (1976).

Rudolph, H. D.

H. D. Rudolph and E. G. Neuman, “Approximations for the eigenvalues of the fundamental mode of a step index glass fiber waveguide,” Nachrichtentech. Elektron. 29, 328–329 (1976).

Saitoh, K.

S. Matsuo, K. Takenaga, Y. Arakawa, Y. Sasaki, S. Tanigawa, K. Saitoh, and M. Koshiba, “Crosstalk behavior of cores in multi-core fiber under bent condition,” IEICE Electron. Express 8(6), 385–390 (2011).
[Crossref]

K. Takenaga, Y. Arakawa, S. Tanigawa, N. Guan, S. Matsuo, K. Saitoh, and M. Koshiba, “An investigation on crosstalk in multi-core fibers by introducing random fluctuation along longitudinal direction,” IEICE Trans. Commun. E 94-B, 409–416 (2011).

K. Takenaga, Y. Arakawa, Y. Sasaki, S. Tanigawa, S. Matsuo, K. Saitoh, and M. Koshiba, “A large effective area multi-core fiber with an optimized cladding thickness,” Opt. Express 19(26), B543–B550 (2011).
[Crossref] [PubMed]

M. Koshiba, K. Saitoh, K. Takenaga, and S. Matsuo, “Multi-core fiber design and analysis: coupled-mode theory and coupled-power theory,” Opt. Express 19(26), B102–B111 (2011).
[Crossref] [PubMed]

S. Matsuo, K. Takenaga, Y. Arakawa, Y. Sasaki, S. Taniagwa, K. Saitoh, and M. Koshiba, “Large-effective-area ten-core fiber with cladding diameter of about 200 μm,” Opt. Lett. 36(23), 4626–4628 (2011).
[Crossref] [PubMed]

M. Koshiba, K. Saitoh, and Y. Kokubun, “Heterogeneous multi-core fibers: proposal and design principle,” IEICE Electron. Express 6(2), 98–103 (2009).
[Crossref]

K. Saitoh and M. Koshiba, “Full-vectorial imaginary-distance beam propagation method based on a finite element scheme: Application to photonic crystal fibers,” IEEE J. Quantum Electron. 38(7), 927–933 (2002).
[Crossref]

Sasaki, Y.

Takenaga, K.

K. Takenaga, Y. Arakawa, Y. Sasaki, S. Tanigawa, S. Matsuo, K. Saitoh, and M. Koshiba, “A large effective area multi-core fiber with an optimized cladding thickness,” Opt. Express 19(26), B543–B550 (2011).
[Crossref] [PubMed]

S. Matsuo, K. Takenaga, Y. Arakawa, Y. Sasaki, S. Taniagwa, K. Saitoh, and M. Koshiba, “Large-effective-area ten-core fiber with cladding diameter of about 200 μm,” Opt. Lett. 36(23), 4626–4628 (2011).
[Crossref] [PubMed]

M. Koshiba, K. Saitoh, K. Takenaga, and S. Matsuo, “Multi-core fiber design and analysis: coupled-mode theory and coupled-power theory,” Opt. Express 19(26), B102–B111 (2011).
[Crossref] [PubMed]

K. Takenaga, Y. Arakawa, S. Tanigawa, N. Guan, S. Matsuo, K. Saitoh, and M. Koshiba, “An investigation on crosstalk in multi-core fibers by introducing random fluctuation along longitudinal direction,” IEICE Trans. Commun. E 94-B, 409–416 (2011).

S. Matsuo, K. Takenaga, Y. Arakawa, Y. Sasaki, S. Tanigawa, K. Saitoh, and M. Koshiba, “Crosstalk behavior of cores in multi-core fiber under bent condition,” IEICE Electron. Express 8(6), 385–390 (2011).
[Crossref]

Taniagwa, S.

Tanigawa, S.

K. Takenaga, Y. Arakawa, Y. Sasaki, S. Tanigawa, S. Matsuo, K. Saitoh, and M. Koshiba, “A large effective area multi-core fiber with an optimized cladding thickness,” Opt. Express 19(26), B543–B550 (2011).
[Crossref] [PubMed]

K. Takenaga, Y. Arakawa, S. Tanigawa, N. Guan, S. Matsuo, K. Saitoh, and M. Koshiba, “An investigation on crosstalk in multi-core fibers by introducing random fluctuation along longitudinal direction,” IEICE Trans. Commun. E 94-B, 409–416 (2011).

S. Matsuo, K. Takenaga, Y. Arakawa, Y. Sasaki, S. Tanigawa, K. Saitoh, and M. Koshiba, “Crosstalk behavior of cores in multi-core fiber under bent condition,” IEICE Electron. Express 8(6), 385–390 (2011).
[Crossref]

IEEE J. Quantum Electron. (1)

K. Saitoh and M. Koshiba, “Full-vectorial imaginary-distance beam propagation method based on a finite element scheme: Application to photonic crystal fibers,” IEEE J. Quantum Electron. 38(7), 927–933 (2002).
[Crossref]

IEICE Electron. Express (2)

S. Matsuo, K. Takenaga, Y. Arakawa, Y. Sasaki, S. Tanigawa, K. Saitoh, and M. Koshiba, “Crosstalk behavior of cores in multi-core fiber under bent condition,” IEICE Electron. Express 8(6), 385–390 (2011).
[Crossref]

M. Koshiba, K. Saitoh, and Y. Kokubun, “Heterogeneous multi-core fibers: proposal and design principle,” IEICE Electron. Express 6(2), 98–103 (2009).
[Crossref]

IEICE Trans. Commun. E (1)

K. Takenaga, Y. Arakawa, S. Tanigawa, N. Guan, S. Matsuo, K. Saitoh, and M. Koshiba, “An investigation on crosstalk in multi-core fibers by introducing random fluctuation along longitudinal direction,” IEICE Trans. Commun. E 94-B, 409–416 (2011).

J. Lightwave Technol. (1)

Nachrichtentech. Elektron. (1)

H. D. Rudolph and E. G. Neuman, “Approximations for the eigenvalues of the fundamental mode of a step index glass fiber waveguide,” Nachrichtentech. Elektron. 29, 328–329 (1976).

Opt. Express (2)

Opt. Lett. (1)

Other (6)

A. W. Snyder and J. D. Love, Optical Waveguide Theory (Chapman & Hall, 1983), Chap. 37.

K. Takenaga, Y. Arakawa, S. Tanigawa, N. Guan, S. Matsuo, K. Saitoh, and M. Koshiba, “Reduction of crosstalk by trench-assisted multi-core fiber,” in Optical Fiber Communication Conference, OSA Technical Digest (CD) (Optical Society of America, 2011), paper OWJ4.

T. Morioka, “New generation optical infrastructure technologies: “EXACT initiative” towards 2020 and beyond,” in Proceedings of 14th OptoElectronics and Communications Conference (Institute of Electrical and Electronics Engineers, 2009), paper FT4.

K. Takenaga, S. Tanigawa, N. Guan, S. Matsuo, K. Saitoh, and M. Koshiba, “Reduction of crosstalk by quasi-homogeneous solid multi-core fiber,” in Optical Fiber Communication Conference, OSA Technical Digest (CD) (Optical Society of America, 2010), paper OWK7.

T. Hayashi, T. Nagashima, O. Shimakawa, T. Sasaki, and E. Sasaoka, “Crosstalk variation of multi-core fiber due to fiber bend,” in Proceedings of 36th European Conference and Exhibition on Optical Communication (Institute of Electrical and Electronics Engineers, 2010), paper We.8.F.6.

K. Okamoto, Fundamentals of Optical Waveguides (Corona Publishing, 1992), Chap. 4.

Cited By

OSA participates in Crossref's Cited-By Linking service. Citing articles from OSA journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (11)

Fig. 1
Fig. 1 The profile of refractive index in two trench-assisted non-identical cores and the part outside the cores.
Fig. 2
Fig. 2 Difference of the refractive-index distributions. (a) N2N22. (b) N2N12.
Fig. 3
Fig. 3 The profile of core m with trench structures.
Fig. 4
Fig. 4 Geometries for the calculation of the coupling coefficient.
Fig. 5
Fig. 5 Schematic of a core with index trench and Homo-TA-7-core model.
Fig. 6
Fig. 6 Simulated crosstalk at 1550-nm wavelength as function of length.
Fig. 7
Fig. 7 Schematic of Hetero-TA-12-core model.
Fig. 8
Fig. 8 Required Δneff as function of the Λ and Rpk.
Fig. 9
Fig. 9 Effective index value of the fundamental mode as function of core radius and core Δ1, where (a) r2/r1 = 2.0, W/r1 = 1.0, (b) r2/r1 = 2.0, W/r1 = 1.1, (c) r2/r1 = 2.0, W/r1 = 1.2, and (d) r2/r1 = 2.0, W/r1 = 1.3.
Fig. 10
Fig. 10 Crosstalk of Hetero-TA-12-core fiber at 100-km propagation as function of bending radius.
Fig. 11
Fig. 11 Crosstalk of Homo-TA-12-core fiber at 100-km propagation as function of bending radius.

Tables (2)

Tables Icon

Table 1 Structural Parameters for Calculation

Tables Icon

Table 2 Optical Properties of the Cores in Different Conditions (1550 nm)

Equations (34)

Equations on this page are rendered with MathJax. Learn more.

N 2 (r,θ)= N 1 2 (r,θ)+ N 2 2 (r,θ) n 2 (r,θ),
κ pq = ω ε 0 + + ( N 2 N q 2 ) E p * E q dxdy + + u z ( E p * × H p + E p × H p * )dxdy ,
κ pq = ω ε 0 4P 0 2π 0 a 1p ( n p 2 n cl 2 ) E p * E q rdrdθ ,
E ˜ =E(r,θ) e j(ωtβz) .
2 E z r 2 + 1 r E z r + 1 r 2 2 E z θ 2 +[ k 2 n (r,θ) 2 β 2 ] E z =0.
ξ m = n m 2 k 2 β m 2 ,
σ m = β m 2 n cl 2 k 2 ,
γ m = β m 2 n tr 2 k 2 ,
V 1m = a 1m k n m 2 n cl 2 ,
V 2m = a 1m k n cl 2 n tr 2 ,
W 1m = a 1m σ m =1.1428 V 1m 0.996,
U 1m = a 1m ξ m = V 1m 2 W 1m 2 ,
W 2m = a 1m γ m = V 2m 2 + W 1m 2 ,
E m z ={ A m J n ( ξ m r m )cos(nθ+ψ)(inCo m) B m K n ( σ m r m )cos(nθ+ψ)(inIC m) C m K n ( γ m r m )cos(nθ+ψ)(inTrm) D m K n ( σ m r m )cos(nθ+ψ)(inOC) E m K n ( γ m r m )cos(nθ+ψ)(inTrm') F m K n ( σ m r m )cos(nθ+ψ)(inIC m') ,
{ A m J n ( U 1m )= B m K n ( W 1m ) B m K n ( W 1m a 2m a 1m )= C m K n ( W 2m a 2m a 1m ) C m K n ( W 2m a 3m a 1m )= D m K n ( W 3m a 3m a 1m ) D m K n ( W 3m R 1 a 1m )= E m K n ( W 4m R 1 a 1m ) E m K n ( W 4m R 2 a 1m )= F m K n ( W 5m R 2 a 1m ) .
D m = L m A m ,
F m = Q m A m ,
L m = J n ( U 1m ) K n ( W 1m a 2m a 1m ) K n ( W 2m a 3m a 1m ) K n ( W 1m ) K n ( W 2m a 2m a 1m ) K n ( W 3m a 3m a 1m ) ,
Q m = L m K n ( W 3m R 1 a 1m ) K n ( W 4m R 2 a 1m ) K n ( W 4m R 1 a 1m ) K n ( W 5m R 2 a 1m ) .
A m = U 1m W 1m β m a 1m 2 V 1m J 1 ( U 1m ) 2P π ε 0 n 1m c ,
E p ={ E p x =j A p β p a 1p U 1p J 0 ( U 1p r a 1p )cosψ E p y =j A p β p a 1p U 1p J 0 ( U 1p r a 1p )sinψ E p z = A p J 1 ( U 1p r a 1p )cos(θ+ψ) ,
E q ={ E q x =j Q q A q β q a 1q W 1q K 0 ( W 1q R a 1q )cosψ E q y =j Q q A q β q a 1q W 1q K 0 ( W 1q R a 1q )sinψ E q z = Q q A q K 1 ( W 1q R a 1q )cos(Θ+ψ) .
R= D 2 + r 2 2Drcosθ Drcosθ,
a 3p = D 2 + R 1 2 2D R 1 cos(πΘ) D R 1 cos(πΘ),
r= D 2 + R 2 2DRcos(πΘ) DRcos(πΘ),
R 1 (D a 3p )(Drcosθ) Dr .
R 2 (D a 2p )(Drcosθ) Dr .
E p * E q = Q q A p A q β p β q a 1p a 1q U 1p W 1q J 0 ( U 1p r a 1p ) K 0 ( W 1q R a 1q ) + Q q A p A q J 1 ( U 1p r a 1p ) K 1 ( W 1q R a 1q )cos(θ+ψ)cos(Θ+ψ).
S 1 = 0 2π 0 a 1p ( n p 2 n cl 2 ) E p * E q rdrdθ =( n p 2 n cl 2 ) L q A p A q β p β q a 1p a 1q U 1p W 1q × 0 2π 0 a 1p J 0 ( U 1p r a 1p ) K 0 ( W 1q R a 1q ) K 1 ( P 1 Drcosθ Dr ) K 1 ( Y 1 Drcosθ Dr ) K 1 ( P 2 Drcosθ Dr ) K 1 ( Y 2 Drcosθ Dr ) rdrdθ ,
K n (z) π 2z exp(z).
S 1 =( n p 2 n cl 2 ) L q A p A q β p β q a 1p a 1q U 1p W 1q π a 1q 2 W 1q D exp( W 1q D a 1q ) × 0 a 1p J 0 ( U 1p r a 1p )exp[( P 2 P 1 + Y 2 Y 1 ) D Dr ]rdr 0 2π exp[( W 1q a 1q P 2 P 1 + Y 2 Y 1 Dr )rcosθ]dθ .
I 0 (z)= 1 π 0 π exp(zcosθ) dθ,
κ pq = k( n p 2 n cl 2 ) W 1p U 1q L q π a 1q 2 W 1q D exp( W 1q D a 1q ) n p n q a 1p a 1q V 1p V 1q J 1 ( U 1p ) J 1 ( U 1q ) × 0 a 1p J 0 ( U 1p r a 1p ) I 0 [( W 1q a 1q P 2 P 1 + Y 2 Y 1 Dr )r]exp[( P 2 P 1 + Y 2 Y 1 ) D Dr ]rdr.
CMF= N core A eff π (CD/2) 2 ,

Metrics