## Abstract

Thin glass fibers imbedded into a glass cladding of slightly lower refractive index represent a promising medium for optical communication. This article presents simple formulas and functions for the fiber parameters as a help for practical design work. It considers the propagation constant, mode delay, the cladding field depth, and the power distribution in the fiber cross section. Plots vs frequency of these parameters are given for 70 modes

© 1971 Optical Society of America

## I. Introduction

Recently, glass fibers have been produced that permit the transmission of optical signals over several kilometers. In general, these fibers support many modes, which propagate at different velocities. Since this causes signal distortion over long distances, fibers that transmit only a limited number of modes are of special interest. A fiber waveguide consists of a thin central glass core surrounded by a glass cladding of slightly lower refractive index. Most modes can be suppressed by making the core thin and the index different between core and cladding small. Typically, a difference of a few parts in a thousand is feasible. This avoids propagation of most modes. The modes that do propagate are weakly guided, but in general the guidance is sufficient to negotiate bends with radii of tens of centimeters.

Maxwell’s equations have exact solutions for the dielectric cylinder, but even with the simplifying assumption that the cladding be infinitely thick these solutions are too complicated to be evaluated without computer. Recent efforts in simplifying the theory for weakly guided modes had promising results, but in the region of interest, they did not lead to the kind of simple formulas one would wish to have for fiber design work. The following paper is aimed at such formulas and functions. It is meant as a help for engineering applications directed toward fiber communication systems. Most results are valid for all frequencies and propagation conditions—even at cutoff—with an accuracy of the order of the index difference between core and cladding.

## II. Mode Parameters and Characteristic Equation

Consider the cylindric core of radius a depicted in Fig. 1(a). The refractive index of the core is nc,. Let the cladding material of index n extend to infinity. We shall use both artesian coordinates (x,y) and cylindrical coordinates (r,ϕ). The propagation constant β of any mode of this fiber is limited within the interval nckβnk, where k = 2π/λ is the wavenumber in free space. If we define parameters

$u=a(k2nc2-β2)12$
$w=a(β2-k2n2)12,$

the mode field can be expressed by Bessel function J(ur/a) inside the core and modified Hankel function K(wr/a) outside the core. The quadratic summation

$v2=u2+w2$

$v=ak(nc2-n2)12,$

which can be considered as a normalized frequency. By matching the fields at the core–cladding interface, we obtain characteristic functions u(v) or w(v) for every mode; the propagation constant and all other parameters of interest can be derived from these functions.

For weak guidance, we have

$Δ=(nc-n)/n≪1.$

In this case, we can construct modes whose transverse field is essentially polarized in one direction. This can be deduced from the results of Ref. , but in order to elucidate the approximation involved, let us try a direct derivation. We postulate transverse field components

$Ey=Hx {Z0/ncZ0/n}=El {Jl(ur/a)/Jl(u)Kl(wr/a)/Kl(w)} coslϕ.$

Here, as in the following, the upper line holds for the core and the lower line for the cladding; Z0 is the plane wave impedance in vacuum, and El the electrical field strength at the interface. Figure 1 (b)–(e) illustrate the case l = 1. Since we have the freedom of choosing sin or cos in Eq. (6) and two orthogonal states of polarization, we can construct a set of four modes for every l as long as l > 0. For l = 0, we have only a set of two modes polarized orthogonally with respect to each other.

The longitudinal components can be obtained from the equations

$Ez=iZ0k {1/nc21/n2} ∂Hx∂y,$
$and Hz=(i/kZ0)(∂Ey/∂x).$

By introducing Eq. (6), we therefore have

$Ez=-iEl2ka{uncJl+1(ur/a)Jl(u) sin(l+1)ϕ+uncJl-1(ur/a)Jl(u) sin(l-1)ϕwnKl+1(wr/a)Kl(w) sin(l+1)ϕ-wnKl-1(wr/a)Kl(w) sin(l-1)ϕ},$
$Hz=-iEl2kZ0a{uJl+1(ur/a)Jl(u) cos(l+1)ϕ-uJl-1(ur/a)Jl(u) cos(l-1)ϕwKl+1(wr/a)Kl(w) cos(l+1)ϕ+wKl-1(wr/a)Kl(w) cos(l-1)ϕ}.$

For small Δ, the longitudinal components [Eqs. (8a), (8b)] are small compared to the transverse components. The factors involved are u/ak and w/ka which because of Eqs. (1) and (2) are both of the order $Δ12$. Repeated differentiation of Eqs. (8a) and (b) leads to transverse components which are not identical with the postulated field [Eq. (6)] but small of order Δ compared to it. We shall neglect these fields in the following. It is this approximation that determines the accuracy of our assumption of linearly polarized modes.

To match the fields at the interface let us write Eq. (6) in terms of cylindrical components. We then have

$Eϕ=12 El{Jl(ur/a)/Jl(u)Kl(wr/a)/Kl(w)} [cos(l+1)ϕ+cos(l-1)ϕ],$
$Hϕ=-12 ElZ0{ncJl(ur/a)/Jl(u)nKl(wr/a)/Kl(w)}× [sin(l+1)ϕ-sin(l-1)ϕ].$

If we set nc = n in Eqs. (8) and (9) and use the recurrence relations for Jl and Kl, we can match all tangential field components at the interface by the one equation

$u[Jl-1(u)/Jl(u)]=-w[(Kl-1(w)/Kl(w)].$

This is the characteristic equation for the linearly polarized (LP) modes. Setting w = 0 yields the cutoff values Jl−1(u) = 0. For l = 0, this includes the roots of the Bessel function J−1(u) = −J1(u), which we shall count so as to include J1(0) = 0 as the first root. We then obtain the cutoff values indicated in Fig. 2 for LP0m and LP1m. In the limit of w → ∞, we have Jl(u) = 0. Thus, the solutions for u are between the zeros of Jl−1(u) and Jl(u). Every solution is associated with one set of modes designated LPlm. For l ≥ 1, each set comprises four modes.

The accuracy of the characteristic equation can be improved if we retain n and nc as different in Eqs. (8) and (9). In this case, however, terms with (l + 1)ϕ and (l − 1)ϕ satisfy two different characteristic equations:

$(u/nc) [Jl±1(u)/Jl(u)]=±(w/n) [Kl±1(w)/Kl(w)].$

By using the recurrence relations for Jl and Kl, one can easily show that these two equations converge into Eq. (10) for nc = n. For ncn, this degeneracy ceases to exist; each mode LPlm breaks up into modes with terms (l + 1)ϕ, which can be identified as HEl+1,m, and modes with terms (l − 1)ϕ which form EHl−1,m or TEm and TMm., This association is indicated in the lower half of Fig. 2 for the cases l = 0 or 1. A more rigorous proof of the previous results is given in the Appendix, where Eqs. (10) and (11) are derived directly from the exact characteristic equation. The following calculations are based on Eq. (10), which is found to be sufficiently accurate for most practical applications.

## III. Approximate Analytic Solution

By using Eq. (3) and differentiating both sides of Eq. (10) with respect to v, one can write the characteristic Eq. (10) in the form

$du/dv=(u/v)[1-κl(w)],$
$where κl(w)=Kl2(w)/Kl-1(w)Kl+1(w).$

For large w, we have κl ≈ 1 − (1/v). This can be used to solve Eq. (12) for large v. Unfortunately, parameters of interest like the propagation constant or the field depth in the cladding depend on the difference v2u2. As this difference becomes small in the region of interest, the relative error introduced by the above approximation becomes intolerably large. To be useful, an approximation of u must improve toward smaller v.

With this in mind, we replace Eq. (13) by

$κl≈1-(w2+l2+1)-12,$

which not only approximates Eq. (13) for large w, but provides a reasonable fit throughout. We now use Eq. (3) and replace w2 by v2u2. Yet since u stays in the narrow region between successive roots of adjacent Bessel functions, we may write

$w≈(v2-uc2)12,$

replacing u by its cutoff value uc. For the mode LPlm, uc is the mth root of Jl−1(u). The approximations (14) and (15) are satisfactory for all modes except LP01 = HE11, whose mode parameters u, v, and w all approach zero simultaneously.

If we exclude HE11, we can now use Eqs. (14), (15), and the boundary value uc at cutoff to solve Eq. (12). The result is

$u(v)=uc exp[arcsin(s/uc)-arcsin(s/v)]/s,$
$with s=(uc2-l2-1)12.$

In the case of the HE11 mode, a more careful approximation is necessary, although the basic approach is similar to the one outlined previously. Without going into detail, we list the result

$u(v)=(1+√2)v/[1+(4+v4)14] for HE11.$

We mentioned earlier that u is bound between successive zeros of the Bessel functions Jl−1 and Jl. We thus know the asymptotic value u for v → ∞ exactly. Equations (16) and (18) approximate these values within an error of 2%; this is a good indication of the accuracy of Eqs. (16) and (18). For vs (far enough from cutoff), we can reduce Eqs. (16) and (18) to

$u(ν)=u∞[1-(1/v)]$

for all modes, using the mth root of Jl(u) for u.

## IV. Propagation Constant and Mode Delay

With the help of u(v), we can calculate the propagation constant β from Eq. (1). In order to make the results independent of particular fiber configurations, however, we shall not plot β directly but the ratio

$b(v)=1-(u2/v2)=[(β2/k2)-n2]/(nc2-n2),$

which, for small index difference, reduces to

$b≈[(β/k)-n]/(nc-n).$

From this and Eq. (5) we obtain the propagation constant

$β=nk(bΔ+1)=nk[1+Δ-Δ(u2/v2)].$

Since β and b are proportional, the quantity b can be understood as a normalized propagation constant. Figure 3 shows b(v) for 18 LP modes. A comparison with exact computer solutions showed deviations that were too small to be displayed in Fig. 3.

Direct detection of intensity-modulated light signals recognizes dispersion effects only in the envelope of the light signal. This envelope is influenced by the group delay

$τgr=(L/c)(dβ/dk).$

Here c is the vacuum velocity of light and L the length of the fiber. When we differentiate Eq. (22), we must consider the k dependence of n, Δ, and b. Yet, if the dispersion of the core and the cladding glass is approximately the same, Δ is independent of k. Moreover, for all glasses, kdn/dkn. If we ignore products of Δ with (k/n) (dn/dk), we obtain

$τgr=Lc {[d(nk)/dk]+nΔ[d(vb)/dv]}.$

The first part of Eq. (24) characterizes the material dispersion, which is the same for all modes. The second part, which represents the group delay on account of waveguide dispersion, is governed by the derivative d(vb)/dv. Because of Eqs. (12) and (20), this derivative can be expressed by

$d(vb)/dv=1-(u/v)2(1-2κ),$

with κ from Eq. (13). This function is plotted in Fig. 4. Far from cutoff, it approaches unity for all modes. At cutoff, d(vb)/dv = 2κ(0). This results in cutoff values d(vb)/dv = 0 for l = 0, 1, and 2[1 − (1/l)] for l ≥ 2. As Fig. 4 shows, the mode of largest order has the largest group delay. The difference between this and the slowest mode is approximately 1 − (2/l). For large v, we have lmaxv. Thus we obtain a group spread [(1 − (2/v)](n1n2)L/c for a fiber that propagates many modes.

## V. Power Flow and Power Density

The Poynting vector in axial direction can be calculated from the cross product of the transverse fields given in Eq. (6). Integration over the cross section of core and cladding leads to tabulated integrals,; the results are

$Pcore=[1+(w2/u2)(1/κ)] (πa2/2)(Z0/nc)El2$
$and Pclad=[(1/κ)-1] (πa2/2)(Z0/n)El2$

for the power flow in core and cladding, respectively. If we ignore the small difference between nc and n, the total power in a certain mode becomes

$P=Pcore+Pclad=(v2/u2)(1/κ)(πa2/2)(Z0/n)El2.$

Practical fibers have small heat and scattering losses which cause significant attenuation over long distances. In general, these losses are attributable to certain parts of the fiber and proportional to the power propagating in this part. For considerations of this kind, it is convenient to use the power fractions

$Pcore/P=1-(u2/v2)(1-κ)$
$and Pclad/P=(u2/v2)(1-κ),$

which are plotted in Fig. 5. As expected, the mode power is concentrated in the core far away from cutoff. As cutoff is approached, the power of low order modes (l = 0,1) withdraws into the cladding, whereas modes with l ≥ 2 maintain a fixed ratio of l − 1 between the power in core and cladding at cutoff.

The power density is related to the mode power P by

$p(r)=κu2v2 2Pπa2 {Jl2(ur/a)/Jl2(u)Kl2(wr/a)/Kl2(w)} cos2lϕ.$

By averaging over ϕ at r = a, we obtain the mean density

$p¯(r)=κ(u2/v2)Pπa2 {Jl2(ur/a)/Jl2(u)Kl2(wr/a)/Kl2(w)}.$

At the core–cladding interface, we have r = a and

$p¯(a)=κ(u2/v2)(P/πa2).$

The normalized density πa2 $P¯$(a)/P is plotted in Fig. 6. For modes of order l = 0,1 this density approaches zero both at cutoff and far away from it, having a maximum in between. Modes with l ≥ 2 have $P¯$(a) = [1 − (1/l)]P/πa2 at cutoff.

For ra/w, we can replace the K functions in Eq. (31) by their approximation for large argument and obtain

$p¯(r)≈κ(u2/v2)(P/πar) exp[-2w(r-a)/a] for r≫a,$

as long as w is not too small. The power density decreases exponentially with the distance from the interface. The parameter w is plotted in Fig. 7. It decreases sharply as cutoff is approached and is zero at cutoff. For sufficiently small w we may set u = v and replace the K functions in Eq. (31) by their approximation for small argument, obtaining

$p¯(r)≈κl(P/πa2)(a/r)l for r>a, w=0.$

This function describes the cutoff power distribution in the cladding. It decreases with the distance from the axis for all but the lowest azimuthal order, whose cladding field is independent of the radius. These results are of course based on our theoretical model of a core imbedded into an unbounded cladding.

## VI. Approximations for Multimode Fibers

Fibers that transmit a large number of modes are of particular interest in connection with incoherent light sources, as, for example, light-emitting diodes. This is because the amount of light that the fiber accepts from this source increases with the number of modes it transmits. Clearly, fibers with a large mode volume neither can nor need be evaluated in as much detail as was done in previous sections. A much simpler mode picture is required. To find this, let us introduce a somewhat simpler way of counting the modes.

Figure 8 illustrates the front end of a fiber and the cone of light that this fiber accepts. The cone is limited by those rays that, after entering the front end of the core, are totally reflected at the interface between core and cladding. These rays form an angle

$θ≈sinθ=(nc2-n2)12$

with the axis. Consider now the free-space modes which can enter the front area πa2 of the core. There are pairs of modes that are perpendicularly polarized with respect to each other. Each pair occupies a cone of solid angle πδ2, where

$δ=λ/πa.$

The total number of free-space modes accepted by the fiber is consequently

$N≈2(θ/δ)2$

or, because of Eq. (4),

$N≈v2/2.$

This is also the number of modes transmitted by the fiber. We find this confirmed if we count all cutoff values uc < v, bearing in mind that the lowest value represents two modes and all others four. Equation (36) permits us to label the fiber modes in the sequence of their cutoff values. Vice versa, we can predict the cutoff of the νth mode approximately at

$uc≈(2ν)12.$

This rule, of course, holds only for large ν.

Let us now use this counting method to describe fibers with large mode volume. We restrict ourselves to operation far from cutoff, ignoring, among the large number of modes, those few that are appreciably close to cutoff. We thus set κ = 1. We replace u by uc, bearing in mind that the parameter u never departs much from its cutoff value uc. Because of Eq. (20), the propagation parameter b then becomes

$b=1-(uc2/v2)=1-(ν/N).$

Equation (25) can be simplified accordingly to yield the group delay parameter

$d(vb)/dv=1+(uc2/v2)=1+(ν/N).$

We can estimate the power density at the interface by converting Eq. (32). This yields

$p¯(a)=(P/πa2)(uc2/v2)=(P/πa2)(ν/N).$

The power flow in the cladding can be obtained from Eq. (30) by using the approximation $1-κ≈(v2-uc2)-12$, which is obtained from Eq. (14). Equation (30) then reads

$Pclad≈Puc2/v2(v2-uc2)12=Pν/N(2N-2ν)12.$

The incoherent source, in general, excites every fiber mode with the same amount of power. We find the average power distribution in a fiber under these conditions by averaging Eqs. (40) and (41) over all N modes, treating the mode index ν as a continuous variable. We then obtain for the density at the interface

$[p¯(a)/P]tot=1πa2N∫0NνNdν=12πa2.$

It is interesting to note that this quantity is independent of N or the parameter v. The average cladding power follows from an integration of Eq. (41). The result is

$(Pclad/P)tot=1N∫0NνdνN(2N-2ν)12=43N-12.$

The power flow in the cladding decreases proportional to $N-12$ or, because of Eq. (36), proportional to 1/v.

As an example, consider a fiber with a core index of 1.5, Δ = 0.003, and a core radius of 25 μ. This results in v = 20.3 at 0.9-μ wavelength. The group delay between the fastest and the slowest mode is 13.5 nsec after 1 km. The fiber transmits 206 modes. Roughly 10% of the power propagates in the cladding.

## VII. Conclusions

The theory of dielectric waveguides greatly simplifies if the difference between the refractive indices of core and cladding is small. In this case, linearly polarized modes with simple fields and a simple characteristic equation can be defined. Approximate analytic solutions can be derived which are exact at cutoff, showing a maximum relative error of up to 2% for very high frequencies. The mode parameters that follow from these solutions have sufficient accuracy for most practical applications. We consider the propagation constant (Fig. 3), mode delay (Fig. 4), the cladding field depth (Fig. 7), and the power distribution in the fiber cross section (Figs. 5 and 6). Far from cutoff, the group velocity of all modes is smaller than the plane wave velocity in the core and decreases as cutoff is approached. For certain modes, however, this trend reverses shortly before cutoff is reached. The fraction of mode power that propagates in the cladding is the larger the closer the mode to cutoff. For most modes, however, it is just a small fraction of the total power even at cutoff. The modes with the two lowest azimuthal orders are an exception. Their cutoff distributions are characterized by plane wave fields in the cladding which contain practically all the mode power. In multimode fibers with incoherent input, the power is in general equally distributed among all modes. In this case, we find that the power density at the interface is independent of the number of modes, but the power flow in the cladding decreases proportional to the root of that number.

## Appendix: Another Derivation of the Simplified Characteristic Equations

As mentioned earlier, the problem of the dielectric cylinder with sharp index step can be solved exactly; using the nomenclature defined in Sec. II, one can write the exact characteristic equation in the form

$(Q-D-2Δ{[(l±1)/ω2]± [Kl(ω)/ωKl±1(ω)]})(Q-D)=Q2[1-2Δ(u2/v2)],$

where

$Q=(l±1)(v2/u2w2),$
$D=[Jl(u)/uJl±1(u)]∓[Kl(w)/wKl±1(w)],$
$and 2Δ=(nc2-n2)/nc2.$

The upper sign holds for HEl+1 modes and the lower sign for EHl−1, TM, and TE. Equation (A4) agrees with Eq. (5) in the case of small index differences ncn. If Δ is set to zero in Eq. (A1), we find D = 0; Eq. (A3) then becomes the simplified characteristic Eq. (10). For small Δ, D is also small. Let us now simplify Eq. (A1) to the extent that we retain terms linear in Δ or D. This results in

$D=Δ{Q(u2/v2)-[(l±1)/w2]∓[Kl(w)/wKl±1(w)]}$

and with Eq. (A2)

$D=∓Δ[Kl(w)/wKl±1(w)].$

By introducing this into Eq. (A3) and inserting Eq. (5), we find

$(u/nc) [Jl±1(u)/Jl(u)]=±(w/n)[Kl±1(w)/Kl(w)].$

This is exactly the characteristic Eq. (11).

I am grateful to E. A. J. Marcatili for helpful suggestions.

## Figures Fig. 1 Sketch of the fiber cross section and the four possible distributions of LP11. Fig. 2 The regions of the parameter u for modes of order l = 0,1. Fig. 3 Normalized propagation parameter b = (β/kn)/(ncn) as a function of the normalized frequency v. Fig. 4 Normalized group delay d(vb)/dv as a function of v. Fig. 5 Portion of the mode power which propagates in the cladding plotted vs v. Fig. 6 Normalized power density at the core–cladding interface plotted vs v. Fig. 7 Cladding parameter w plotted vs v. Fig. 8 Sketch of the fiber front face and the cone of light that the fiber accepts.

1. F. P. Kapron, D. B. Keck, and R. D. Maurer, Appl. Phys. Lett. 17, 423 (1970) [CrossRef]  .

2. E. Snitzer, J. Opt. Soc. Am. 51, 491 (1961) [CrossRef]  .

3. K. C. Kao and G. A. Hockham, Proc. IEE 113, 1151 (1966).

4. E. A. J. Marcatili, Bell Syst. Tech. J. 48, 2103 (1969).

5. A. W. Snyder, IEEE Trans. Microwave Theory Techniques MTT-17, 1130 (1969) [CrossRef]  .

6. S. A. Schelkunoff, Electromagnetic Waves (Van Nostrand, New York, 1943), p. 94.

7. A. W. Snyder, IEEE Trans. Microwave Theory Techniques MTT-17, 1138 (1969) [CrossRef]  .

8. G. Toraldo di Francia, J. Opt. Soc. Am. 59, 799 (1969) [CrossRef]   [PubMed]  .

9. P. Moon and D. E. Spencer, Field Theory Handbook (Springer, Berlin, 1961), p. 192.

### References

• View by:
• |
• |
• |

1. F. P. Kapron, D. B. Keck, R. D. Maurer, Appl. Phys. Lett. 17, 423 (1970).
[Crossref]
2. E. Snitzer, J. Opt. Soc. Am. 51, 491 (1961).
[Crossref]
3. K. C. Kao, G. A. Hockham, Proc. IEE 113, 1151 (1966).
4. E. A. J. Marcatili, Bell Syst. Tech. J. 48, 2103 (1969).
5. A. W. Snyder, IEEE Trans. Microwave Theory Techniques MTT-17, 1130 (1969).
[Crossref]
6. S. A. Schelkunoff, Electromagnetic Waves (Van Nostrand, New York, 1943), p. 94.
7. A. W. Snyder, IEEE Trans. Microwave Theory Techniques MTT-17, 1138 (1969).
[Crossref]
8. G. Toraldo di Francia, J. Opt. Soc. Am. 59, 799 (1969).
[Crossref] [PubMed]
9. P. Moon, D. E. Spencer, Field Theory Handbook (Springer, Berlin, 1961), p. 192.

#### 1970 (1)

F. P. Kapron, D. B. Keck, R. D. Maurer, Appl. Phys. Lett. 17, 423 (1970).
[Crossref]

#### 1969 (4)

E. A. J. Marcatili, Bell Syst. Tech. J. 48, 2103 (1969).

A. W. Snyder, IEEE Trans. Microwave Theory Techniques MTT-17, 1130 (1969).
[Crossref]

A. W. Snyder, IEEE Trans. Microwave Theory Techniques MTT-17, 1138 (1969).
[Crossref]

#### 1966 (1)

K. C. Kao, G. A. Hockham, Proc. IEE 113, 1151 (1966).

#### Hockham, G. A.

K. C. Kao, G. A. Hockham, Proc. IEE 113, 1151 (1966).

#### Kao, K. C.

K. C. Kao, G. A. Hockham, Proc. IEE 113, 1151 (1966).

#### Kapron, F. P.

F. P. Kapron, D. B. Keck, R. D. Maurer, Appl. Phys. Lett. 17, 423 (1970).
[Crossref]

#### Keck, D. B.

F. P. Kapron, D. B. Keck, R. D. Maurer, Appl. Phys. Lett. 17, 423 (1970).
[Crossref]

#### Marcatili, E. A. J.

E. A. J. Marcatili, Bell Syst. Tech. J. 48, 2103 (1969).

#### Maurer, R. D.

F. P. Kapron, D. B. Keck, R. D. Maurer, Appl. Phys. Lett. 17, 423 (1970).
[Crossref]

#### Moon, P.

P. Moon, D. E. Spencer, Field Theory Handbook (Springer, Berlin, 1961), p. 192.

#### Schelkunoff, S. A.

S. A. Schelkunoff, Electromagnetic Waves (Van Nostrand, New York, 1943), p. 94.

#### Snyder, A. W.

A. W. Snyder, IEEE Trans. Microwave Theory Techniques MTT-17, 1138 (1969).
[Crossref]

A. W. Snyder, IEEE Trans. Microwave Theory Techniques MTT-17, 1130 (1969).
[Crossref]

#### Spencer, D. E.

P. Moon, D. E. Spencer, Field Theory Handbook (Springer, Berlin, 1961), p. 192.

#### Appl. Phys. Lett. (1)

F. P. Kapron, D. B. Keck, R. D. Maurer, Appl. Phys. Lett. 17, 423 (1970).
[Crossref]

#### Bell Syst. Tech. J. (1)

E. A. J. Marcatili, Bell Syst. Tech. J. 48, 2103 (1969).

#### IEEE Trans. Microwave Theory Techniques (2)

A. W. Snyder, IEEE Trans. Microwave Theory Techniques MTT-17, 1130 (1969).
[Crossref]

A. W. Snyder, IEEE Trans. Microwave Theory Techniques MTT-17, 1138 (1969).
[Crossref]

#### Proc. IEE (1)

K. C. Kao, G. A. Hockham, Proc. IEE 113, 1151 (1966).

#### Other (2)

S. A. Schelkunoff, Electromagnetic Waves (Van Nostrand, New York, 1943), p. 94.

P. Moon, D. E. Spencer, Field Theory Handbook (Springer, Berlin, 1961), p. 192.

### Cited By

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

### Figures (8)

Fig. 1 Sketch of the fiber cross section and the four possible distributions of LP11.
Fig. 2 The regions of the parameter u for modes of order l = 0,1.
Fig. 3 Normalized propagation parameter b = (β/kn)/(ncn) as a function of the normalized frequency v.
Fig. 4 Normalized group delay d(vb)/dv as a function of v.
Fig. 5 Portion of the mode power which propagates in the cladding plotted vs v.
Fig. 6 Normalized power density at the core–cladding interface plotted vs v.
Fig. 8 Sketch of the fiber front face and the cone of light that the fiber accepts.

### Equations (56)

$u = a ( k 2 n c 2 - β 2 ) 1 2$
$w = a ( β 2 - k 2 n 2 ) 1 2 ,$
$v 2 = u 2 + w 2$
$v = a k ( n c 2 - n 2 ) 1 2 ,$
$Δ = ( n c - n ) / n ≪ 1.$
$E y = H x { Z 0 / n c Z 0 / n } = E l { J l ( u r / a ) / J l ( u ) K l ( w r / a ) / K l ( w ) } cos l ϕ .$
$E z = i Z 0 k { 1 / n c 2 1 / n 2 } ∂ H x ∂ y ,$
$and H z = ( i / k Z 0 ) ( ∂ E y / ∂ x ) .$
$E z = - i E l 2 k a { u n c J l + 1 ( u r / a ) J l ( u ) sin ( l + 1 ) ϕ + u n c J l - 1 ( u r / a ) J l ( u ) sin ( l - 1 ) ϕ w n K l + 1 ( w r / a ) K l ( w ) sin ( l + 1 ) ϕ - w n K l - 1 ( w r / a ) K l ( w ) sin ( l - 1 ) ϕ } ,$
$H z = - i E l 2 k Z 0 a { u J l + 1 ( u r / a ) J l ( u ) cos ( l + 1 ) ϕ - u J l - 1 ( u r / a ) J l ( u ) cos ( l - 1 ) ϕ w K l + 1 ( w r / a ) K l ( w ) cos ( l + 1 ) ϕ + w K l - 1 ( w r / a ) K l ( w ) cos ( l - 1 ) ϕ } .$
$E ϕ = 1 2 E l { J l ( u r / a ) / J l ( u ) K l ( w r / a ) / K l ( w ) } [ cos ( l + 1 ) ϕ + cos ( l - 1 ) ϕ ] ,$
$H ϕ = - 1 2 E l Z 0 { n c J l ( u r / a ) / J l ( u ) n K l ( w r / a ) / K l ( w ) } × [ sin ( l + 1 ) ϕ - sin ( l - 1 ) ϕ ] .$
$u [ J l - 1 ( u ) / J l ( u ) ] = - w [ ( K l - 1 ( w ) / K l ( w ) ] .$
$( u / n c ) [ J l ± 1 ( u ) / J l ( u ) ] = ± ( w / n ) [ K l ± 1 ( w ) / K l ( w ) ] .$
$d u / d v = ( u / v ) [ 1 - κ l ( w ) ] ,$
$where κ l ( w ) = K l 2 ( w ) / K l - 1 ( w ) K l + 1 ( w ) .$
$κ l ≈ 1 - ( w 2 + l 2 + 1 ) - 1 2 ,$
$w ≈ ( v 2 - u c 2 ) 1 2 ,$
$u ( v ) = u c exp [ arcsin ( s / u c ) - arcsin ( s / v ) ] / s ,$
$with s = ( u c 2 - l 2 - 1 ) 1 2 .$
$u ( v ) = ( 1 + √ 2 ) v / [ 1 + ( 4 + v 4 ) 1 4 ] for HE 11 .$
$u ( ν ) = u ∞ [ 1 - ( 1 / v ) ]$
$b ( v ) = 1 - ( u 2 / v 2 ) = [ ( β 2 / k 2 ) - n 2 ] / ( n c 2 - n 2 ) ,$
$b ≈ [ ( β / k ) - n ] / ( n c - n ) .$
$β = n k ( b Δ + 1 ) = n k [ 1 + Δ - Δ ( u 2 / v 2 ) ] .$
$τ g r = ( L / c ) ( d β / d k ) .$
$τ g r = L c { [ d ( n k ) / d k ] + n Δ [ d ( v b ) / d v ] } .$
$d ( v b ) / d v = 1 - ( u / v ) 2 ( 1 - 2 κ ) ,$
$P core = [ 1 + ( w 2 / u 2 ) ( 1 / κ ) ] ( π a 2 / 2 ) ( Z 0 / n c ) E l 2$
$and P clad = [ ( 1 / κ ) - 1 ] ( π a 2 / 2 ) ( Z 0 / n ) E l 2$
$P = P core + P clad = ( v 2 / u 2 ) ( 1 / κ ) ( π a 2 / 2 ) ( Z 0 / n ) E l 2 .$
$P core / P = 1 - ( u 2 / v 2 ) ( 1 - κ )$
$and P clad / P = ( u 2 / v 2 ) ( 1 - κ ) ,$
$p ( r ) = κ u 2 v 2 2 P π a 2 { J l 2 ( u r / a ) / J l 2 ( u ) K l 2 ( w r / a ) / K l 2 ( w ) } cos 2 l ϕ .$
$p ¯ ( r ) = κ ( u 2 / v 2 ) P π a 2 { J l 2 ( u r / a ) / J l 2 ( u ) K l 2 ( w r / a ) / K l 2 ( w ) } .$
$p ¯ ( a ) = κ ( u 2 / v 2 ) ( P / π a 2 ) .$
$p ¯ ( r ) ≈ κ ( u 2 / v 2 ) ( P / π a r ) exp [ - 2 w ( r - a ) / a ] for r ≫ a ,$
$p ¯ ( r ) ≈ κ l ( P / π a 2 ) ( a / r ) l for r > a , w = 0.$
$θ ≈ sin θ = ( n c 2 - n 2 ) 1 2$
$δ = λ / π a .$
$N ≈ 2 ( θ / δ ) 2$
$N ≈ v 2 / 2.$
$u c ≈ ( 2 ν ) 1 2 .$
$b = 1 - ( u c 2 / v 2 ) = 1 - ( ν / N ) .$
$d ( v b ) / d v = 1 + ( u c 2 / v 2 ) = 1 + ( ν / N ) .$
$p ¯ ( a ) = ( P / π a 2 ) ( u c 2 / v 2 ) = ( P / π a 2 ) ( ν / N ) .$
$P clad ≈ P u c 2 / v 2 ( v 2 - u c 2 ) 1 2 = P ν / N ( 2 N - 2 ν ) 1 2 .$
$[ p ¯ ( a ) / P ] tot = 1 π a 2 N ∫ 0 N ν N d ν = 1 2 π a 2 .$
$( P clad / P ) tot = 1 N ∫ 0 N ν d ν N ( 2 N - 2 ν ) 1 2 = 4 3 N - 1 2 .$
$( Q - D - 2 Δ { [ ( l ± 1 ) / ω 2 ] ± [ K l ( ω ) / ω K l ± 1 ( ω ) ] } ) ( Q - D ) = Q 2 [ 1 - 2 Δ ( u 2 / v 2 ) ] ,$
$Q = ( l ± 1 ) ( v 2 / u 2 w 2 ) ,$
$D = [ J l ( u ) / u J l ± 1 ( u ) ] ∓ [ K l ( w ) / w K l ± 1 ( w ) ] ,$
$and 2 Δ = ( n c 2 - n 2 ) / n c 2 .$
$D = Δ { Q ( u 2 / v 2 ) - [ ( l ± 1 ) / w 2 ] ∓ [ K l ( w ) / w K l ± 1 ( w ) ] }$
$D = ∓ Δ [ K l ( w ) / w K l ± 1 ( w ) ] .$
$( u / n c ) [ J l ± 1 ( u ) / J l ( u ) ] = ± ( w / n ) [ K l ± 1 ( w ) / K l ( w ) ] .$