## Abstract

We present a method for designing optical fibers that support field-flattened, ring-like higher order modes, and show that the effective and group indices of its modes can be tuned by adjusting the widths of the guide’s field-flattened layers or the average index of certain groups of layers. The approach provides a path to fibers that have simultaneously large mode areas and large separations between the propagation constants of their modes.

© 2013 OSA

## 1. Introduction

Most nonlinear propagation artifacts in glass waveguides can be reduced by spreading the power they carry over a large area. Many telecommunications and laser applications, however, require the power to be confined to a single transverse spatial mode. Unfortunately, as a mode’s area increases, its propagation constant approaches those of its neighboring modes, making it susceptible to power cross-coupling and potentially degrading the mode’s spatial or temporal fidelity.

Optical fibers that propagate power in a high-order mode [1,2] offer a path to simultaneously increasing the effective area [3] of a mode and the spacing between the desired mode’s propagation constant and those of its neighbors. Unfortunately, the high-order modes of a step index fiber can have hotspots – regions in their transverse profiles where the local irradiance significantly exceeds the average value – which may make them more susceptible to damage or nonlinear artifacts than modes whose power is relatively uniformly distributed, such as the fundamental.

Optical fibers having a flattened fundamental [4–8] are also attractive, as they spread the propagating power uniformly, and in an amplifier allow for uniform and efficient extraction of energy from the gain medium. Like all waveguides, though, they are bound by a mode size-spacing tradeoff, and we show below that in this regard they are only moderately better than more economically-manufactured conventional guides.

We present here a design methodology that combines the benefits of the two waveguides described above, enabling the construction of flattened high-order modes. Specifically, we provide rules for creating structures that support field-flattened segments, stitching structures that interconnect these segments, and terminating regions that bind (impedance match) the resulting engineered mode to the cladding. The effective index of these modes remains fixed as the number of rings and sizes of the rings are varied, while the effective indices of all other modes, and the group index of the flattened mode, change.

We also present examples of fibers that support field-flattened, ring-like propagation modes and compare their characteristics to the characteristics of high-order modes of step-index fibers, showing that the former are more compact and can have larger intermodal spacings. Though not discussed here, we expect that selective-doping will further discriminate the preferred flattened mode from other modes, as suggested by others for fibers having selectively-doped flattened fundamental modes [9].

## 2. Definitions

In the step-like structures of the following designs, the field’s continuity is enforced between steps by matching the field and its radial derivative across the interfaces. The modes of the guides are analyzed by the transfer matrices of Appendix II and by a separate two-dimensional mode solver that finds the eigenmodes of the scalar Helmholtz equation.

Most of the underlying mathematics and physics are considered in the Appendices. Appendix I presents Bessel solutions to the equation governing axially-symmetric waveguides such as a conventional telecom fiber; its results can be used to determine the refractive indices and thickness of the layers that comprise the flattened, stitching, and termination groups defined below. Appendix II presents transfer matrices that can also be used to determine layer indices and thicknesses, and to determine the properties of all bound modes of the fiber. Appendix III presents closed-form solutions to the mode normalization integral. Appendix IV defines several mode size-spacing products and shows that for a given waveguide these products are fixed, a consequence of the radiance theorem.

#### 2.1 Scaled quantities

A characteristic numerical aperture of the fiber, *NA _{flat}*, is defined as:

*n*is the refractive index of the cladding and

_{clad}*n*is the index of the layer or layers over which the field is to be flattened. The scaled radial coordinate, v, is defined as:

_{flat}*λ*is the wavelength of the guided light and

*r*is the radial coordinate. The scaled refractive index profile,

*η*(v), is defined as:

For the flattened waveguides described here, *n _{flat}* is usually chosen to be the minimum refractive index that can be well controlled. For silica fibers, the flattened layer might be lightly doped with an index-raising dopant such as germanium or doped with a rare-earth element along with index-raising and lowering dopants. Alternatively,

*n*might be pure silica and the cladding might be lightly doped with an index depressing agent such as fluorine; in this case, the dopant only needs to extend to the penetration depth of the desired mode.

_{flat}A layer group’s area-averaged index, 〈*η*〉, is defined as:

*η*and

_{i}*A*represent the scaled index and cross-sectional area of the

_{i}*i*

^{th}layer of the group. In the layer groups defined below, we sometimes constrain this value; 〈

*η*〉 sometimes alters the number of allowed modes or the guide’s intermodal spacings.

Several of the examples that follow list a mode’s scaled effective area and illustrate its scaled field. The scaled area is defined such that the physical area, *A _{eff}*, is given by Eq. (57):

*ψ*, is given from Eq. (50):

*P*

_{0}is the power carried by the mode.

In the following examples, *η* is assumed to range between ± 10, which is achievable for germanium and fluorine-doped silica provided *NA _{flat}* is on the order of 0.06. In silica, other dopants might extend this range moderately, or in phosphate glasses or holey structures, various dopants or air holes can extend this range significantly. Moreover, in holey fibers

*NA*might be controlled to a much smaller value, which would proportionally extend the range of

_{flat}*η*. A larger range of indices is generally advantageous, as it reduces the portion of the guide devoted to the stitching and matching groups described below.

#### 2.2 Flattened layers

A flattened layer is one in which the field does not vary with radius; that is, one where:

*n*=

_{flat}*n*) and the azimuthal order,

_{eff}*l*, must be equal to zero. Furthermore, it is necessary that a flattened layer be joined to appropriate stitching or termination groups, as defined below.

#### 2.3 Stitching groups

A stitching group is a layer or group of layers in which the field’s slope is zero at both endpoints (to match that of the adjacent flattened region) and is predominantly nonzero between those points, usually crossing zero one or more times. This can be accomplished in different ways to produce a variety of mode shapes; several examples are presented here.

Figure 1, Fig. 2, and Fig. 3 illustrate stitching groups that might form a portion of a guide that supports a flattened mode. In the figures, *η _{flat}* is 1 (from Eq. (3) since

*n*(v) =

*n*), the minimum and maximum values of

_{flat}*η*are assumed to fall between ± 10, and the left edge of each group starts at v

_{0}= 0.5

*π*, an arbitrarily chosen value. The thicknesses of the layers that comprise the groups were determined numerically from Bessel solutions to the wave equation, as outlined in Appendix I.

### 2.3.1 Half-wave stitching

Figure 1 illustrates three half-wave stitching groups, that is, three groups that cause the field’s polarity to change sign an odd number of times.

Figure 1(a) shows half-wave stitching accomplished with a single layer. The field changes by a factor of −0.78 as determined by its Bessel solution’s behavior. Simulations show that for a single layer, as the left side of the group is placed at higher values of v_{0}, the ratio of the magnitude of the fields approaches unity and:

*η*is the layer’s scaled index, the numeral one arises from the assumption that the layer is surrounded by field-flattened layers having

*η*= 1, and

*m*is an odd integer. This is also the condition for single layer, half-wave stitching in a one-dimensional slab waveguide (in that case, independent of v

_{0}).

In Fig. 1(b) a second layer is added to make the magnitude of the field to the right of the group the same as the magnitude to its left. We mention without illustration that if the sequence of the layers in Fig. 1(b) is reversed – that is, if the higher index layer is place to the right of the lower index layer – the field on that group’s right can be made an even smaller fraction of the field on its left, when compared to the single layer example of Fig. 1(a).

Figure 1(c) illustrates an evanescent half-wave stitching group, a term that here refers to groups having at least one layer in which the field is the sum of exponentially growing and decaying functions. The thicknesses of the layers that comprise the group are adjusted to also make the 〈*η*〉 = 1 for the group (see Eq. (4)) and to make *ψ* = −1 and *ψ*´ = 0 on the group’s right edge.

### 2.3.2 Full-wave stitching

Figure 2 illustrates three full-wave stitching groups, that is, three groups that cause the field’s polarity to change sign an even number of times.

Figure 2(a) shows half-wave stitching accomplished with a single layer. The field changes by a factor of 0.66 due to its Bessel solution’s behavior. As v_{0} is increased an equation similar to Eq. (8) holds, but whose right-hand side is proportional to an even multiple of *π*.

Figure 2(b) illustrates a two layer full-wave group that returns the field’s magnitude and polarity to their original values. The thickness of the group’s first layer is chosen to make the field zero at the right boundary of the first layer. The thickness and index of the second layer are determined numerically to make *ψ* = 1 and *ψ*´ = 0 on the group’s right edge.

Figure 2(c) illustrates a five-layer evanescent full-wave stitching group. The thickness of the first two layers and a portion of the thickness of the third layer are chosen so that the slope is returned to zero, the field is changed by a factor of −0.707 (*ψ*^{2} drops by a factor of two) within the third layer; we also require that, for the group, 〈*η*〉 = 1 (see Eq. (4)). The thicknesses of the second portion of the third layer and of the remaining two layers are determined in the same fashion, but now with the constraint that *ψ* = 1 and *ψ*´ = 0 on the group’s right edge.

### 2.3.3 Fractional-wave stitching

Figure 3 illustrates three fractional wave stitching groups, that is, three groups that return the field’s slope to zero without allowing the field’s polarity change to sign.

Figure 3(a) illustrates a central stitching layer. The central index is lower than the cladding’s and the field consequently grows exponentially with position; the field on-axis is not zero, here it is 2% of the field at the layer’s edge, and hence it is not classified as a half-wave group. Simulations show that layers such as this can efficiently disrupt the properties of a guide’s non-flattened mode, though their disadvantage is that they carry very little power. Note that the central index of Fig. 3(a) could be made higher than the cladding’s index, resulting in a field similar to that in Fig. 1(a) or Fig. 2(a).

Figure 3(b) illustrates a three layer stitching group in which the field dips but does not pass through zero, creating an effect similar to those created by other structures such as those suggested in [10] (though the index dips of those structures do not return the field or slope to their original values). Simulations suggest that groups such as those in Fig. 3(b) may be difficult to manufacture since their behavior varies strongly with the layers’ thicknesses.

Figure 3(c) illustrates a three layer stitching group in which the field’s magnitude rises within the group. The resulting hotspot may be advantageous for applications where field effects are to be enhanced, but problematic for most other high power laser applications. Like the example of Fig. 3(b), simulations suggest that such a group may be difficult to manufacture.

#### 2.4 Termination groups

A termination group is a layer or group of layers placed between one region of a guide, here most often a region in which the slope of the desired mode’s field is zero, and the guide’s cladding. The indices and thicknesses of the layers that comprise the group are chosen to force the cladding’s exponentially-growing term to zero, and to thus bind the mode to the guide. Termination (binding) is analogous to impedance matching.

The examples of this and the following section give the flattened mode’s scaled effective area and illustrate its scaled field, quantities defined by Eq. (5) and Eq. (6). For example and comparison, consider a step-index fiber that supports the LP_{01} mode and is at the cusp of supporting the LP_{02} mode, that is, v = 1.23*π.* It can be shown that its fundamental mode has a scaled effective area of 37.5; therefore, if the guide’s design operates at *λ* = 1 μm and its core has a numerical aperture of 0.06, its effective area will be 260 μm^{2}. It can also be shown that this mode has a scaled peak field of 0.219 = 1/√20.8. If the fiber carries 1 kW of power its peak field will be 2.61 W^{1/2}/μm and its peak irradiance will be (2.61 W^{1/2}/μm)^{2} = 6.8 W/μm^{2}. Note that the peak irradiance is 1.8 times higher than the simple ratio of the power to the effective area (37.5 ÷ 20.8). For flattened modes, this ratio is closer to unity, typically 1.15.

Figure 4 illustrates three termination groups applied to three flattened waveguides. In the figure, *η _{flat}* is 1 (from Eq. (3) since

*n*(v) =

*n*) and the minimum and maximum values of

_{flat}*η*are limited to ± 10. The thickness of the flattened layer is chosen so that each guide is on the cusp of allowing one axially-symmetric mode beyond the flattened mode. The thicknesses of the layers that comprise the groups were determined numerically from Bessel solutions to the wave equation, applying the constraints listed for each example.

Figure 4(a) illustrates a single-layer termination group. Note that the field extends relatively far into the cladding; at the cladding interface the field is 93% of its value in the flattened region and 21% of the mode’s power is guided in the cladding.

Since the effective index of the guide’s flattened mode is predetermined (because *n _{eff}* =

*n*), the mode’s decay constant in the cladding is fixed and consequently the field in the cladding can only be reduced by reducing the field at the cladding interface – the purpose of the additional layers in Fig. 4(b) and Fig. 4(c).

_{flat}Figure 4(b) illustrates a two-layer termination group, similar to those described in [8]. In this group, the group-averaged scaled index, Eq. (3), serves as an additional constraint; simulations show that it strongly affects the field at the cladding interface. In the example, the layers’ thicknesses are varied to make the field at the cladding boundary 50% of the field in the flattened layer (this occurs with the group’s average index, Eq. (4), set to 〈*η*〉 = 0.7), and to match the field’s slope at the cladding interface. Roughly 7% of the mode’s power is guided in the cladding.

Figure 4(c) illustrates a three-layer termination group. The field is set to zero at the interface between the first and second layer, the local minima in the second layer is 50% of the field in the flattened layer, and the group-averaged index, Eq. (4), is set to 〈*η*〉 = 0.7. The field at the cladding interface is −3% of the field in the flattened region, and 0.04% of the mode’s power is guided in the cladding, though now a significant power-fraction is guided by the termination group.

Termination groups of the type shown in Fig. 4(c) enhance the mode’s confinement but also allow at least one additional axially-symmetric mode, plus the asymmetric modes that may accompany it. Relative to the desired mode, the additional modes can have very different propagation constants, very different transverse power distributions, or both; thus, they may not readily couple to the desired mode and may not be problematic.

## 3. Example waveguides

Waveguides that propagate a flattened high order mode are created by interleaving flattening layers with stitching groups, typically starting from the inside of the guide and working outward, then binding the mode to the cladding with a termination group.

Table 1 lists designs for three waveguides; A and B both support a three-ringed, flattened mode, and C supports several higher-order modes. A and B each have three flattened layers (*i*, *v* and *ix*), two three-layer half-wave stitching groups similar to those illustrated in Fig. 1(c) (*ii-iv* and *vi-viii*), and a two-layer termination group similar to the one in Fig. 4(b) (*x*-*xi*). Surrounding these layers is the cladding having *η* = 0.

In Design A the flattened layers have equal cross-sectional areas, both stitching groups have 〈*η*〉 = 3.0, and the termination group has 〈*η*〉 = 0.7. In Design B the flattened layers have equal widths, both stitching groups have 〈*η*〉 = 2.4, and the termination group has 〈*η*〉 = 0.7.

We compare the flattened LP_{03} modes of Designs A and B to the LP_{03} mode of a few-mode step index design, Design C. Design C is similar to the high-order mode fibers reported by others [2], but has a smaller v-number to make its mode count similar to those of A and B.

Figure 5 shows line-outs of the field distributions for the three designs. For Design A, the scaled area is 140 and the scaled peak field is 1/√122; for Design B the values are 150 and 1/√134; and for Design C the values are 140 and 1/√30.8. Thus for equivalent areas, the peak intensity of the modes in Designs A and B would be one-fourth the peak intensity of the mode of Design C, greatly reducing the damage threshold. (The large disparity between the two measures of mode size for C – 140 for its effective area vs. 30.8 for the reciprocal of its peak irradiance, a ratio of 4.5 – is due to its central hotspot.)

Figure 6 shows the transverse field distributions of the LP_{03} and LP_{13} modes of the three designs; when bent, the LP_{03}’s will morph toward their respective LP_{13}’s. Note that the power is more compactly packed in the flattened modes than in the step-index mode.

Note, too, that the inner rings of LP_{13} modes of the flattened designs have essentially the same diameter as the inner rings of their corresponding LP_{03} modes. The inner ring of the LP_{13} mode for the step-index design, though, has a substantially larger diameter than its corresponding LP_{03} mode. This suggests the latter’s mode will experience a larger shift in its centroid when that fiber is bent. The design of the high-order mode fiber in [2] has a central spike in its index profile, perhaps to keep its mode centered.

Figure 7 compares the size-spacing products, Θ* _{eff}* (the phase index-area spacing – essentially the radiance), defined by Eq. (59) in Appendix IV, for the modes of the three designs. The size-spacing products are an invariant of a design. Larger values are likely preferable, since they imply that larger-sized modes may be fabricated while keeping the intermodal spacing constant, and thus keeping the likelihood of intermodal coupling constant. Bear in mind that the effective area term in the Θ

*equation is the same for all of a design’s modes; for each design here, it is chosen to be the area of the design’s LP*

_{eff}_{03}mode.

For Designs A and B, the spacing between the Θ* _{eff}*’s for the three highest-order symmetric modes, the LP

_{02}, LP

_{03}(flattened mode) and LP

_{04}(on the cusp of existence), have been made equal by choosing an appropriate thickness for the flattened layers and by choosing an appropriate value of 〈

*η*〉 (Eq. (3)) for each design’s stitching groups.

For A and B, the size-spacing differential for the axially-symmetric modes is 2.5 times larger than it is for Design C, and three times larger than for the designs in Fig. 4. This implies that for the same manufacturing tolerances, the three-ringed flattened design can have 2.5 times the area of C, or three times the area of the designs in Fig. 4.

Note that the effective index spectra of A and B are strongly affected by the relative widths of the flattened layers; a relatively large spacing has been created between the LP_{12} and LP_{22} modes of B (red arrow in Fig. 7(b)).

Figure 8 compares the size-spacing products, Θ* _{g}* (group index-area spacing), defined by Eq. (62) in Appendix IV, for the modes of the three designs. The size-spacing products are an invariant of a design. Larger values are likely preferable, since they imply that larger-sized modes may be fabricated while the keeping the intermodal spacing constant, and thus keeping the likelihood of intermodal coupling constant. Keep in mind that the effective area term in Θ

*equation is the same for all of a design’s modes; for each design, it is chosen to be the area of the design’s LP*

_{g}_{03}mode.

Note that the group index spacings of the two flattened designs, A and B, are significantly larger than those of the step-index design, C; the larger spacings may help reduce linear and nonlinear modal coupling in pulsed laser applications. Simulations show that the group delay spectra of A and B are strongly affected by the relative widths of the flattened layers. Note that a local maxima has been created for the LP_{12} mode of B (red arrow in Fig. 8(b)), and that in A and C the flattened mode is the slowest axially-symmetric mode, while in B it is the fastest of all modes.

## 4. Discussion

The design philosophy presented here is atypical – it begins with the desired mode’s shape and then constructs a waveguide that will allow it. Flattening layers are interleaved with stitching groups and a termination group binds the flattened mode to the guide; the latter is analogous to impedance matching. For axially-symmetric waveguides, the thicknesses or indices of the layers that comprise the stitching groups must be changed when the group’s radial placement is changed; the examples presented here should be considered starting points for user-specific designs.

The high-order flattened modes allow two size-spacing invariants – one relating to the phase index spacing, one relating to the group index spacing – to be tailored. In particular, we have shown that the effective index (phase index) spacing of the guide’s axially-symmetric modes can be increased substantially, and show that this spacing grows in proportion to the number of rings added to the structure.

Note that the flattened modes do not suffer potentially problematic hotspots, they inherently pack the propagated power into a compact cross-section, and they may reduce a mode’s susceptibility to some artifacts such as nonlinear self-focusing. In an amplifier, they allow power to be extracted uniformly and efficiently across the mode’s cross section. Furthermore, in amplifier applications the stitching and termination groups would not likely be doped with rare-earth ions, allowing for better control of their indices, and since the field of the flattened mode is near-zero in those regions, avoiding leaving regions of unsaturated gain that might contribute to noise or amplification of undesired modes.

Here we have qualitatively considered the bending properties of the flattened high-order modes by inspecting the transverse structure of the neighboring mode that they would couple to, and find that the flattened modes will stay well-centered. Quantitative bend-loss studies are in progress.

Comparisons to the high-order modes of a step-index fiber are complicated by the fact that the effective area, as conventionally defined, does not account for hotspots in a mode’s peak irradiance. We have used the effective area metric here though in some applications it may give an overly optimistic representation of the performance of high order step-index modes. Despite this (sometimes) lenient metric, the high-order mode of the step-index example fiber is less attractive than the flattened modes described here in terms of intermodal spacing, peak irradiance, and the compactness. While increasing the v-number of the step-index design would improve the intermodal spacing, it would also increase its mode count, accentuate its central hotspot, and further reduce its mode’s packing density.

In principle, flattened high-order modes could be manufactured with conventional telecom techniques such as modified chemical vapor deposition and outside vapor deposition, but the tighter manufacturing tolerances allowed by holey-fiber construction techniques may, however, be preferable or necessary.

## Appendix I: Bessel solutions

Consider the equation that governs the radially-varying portion of the field in an axially symmetric waveguide such as a conventional telecom optical fiber [11]:

*ψ*represents the field of a guided mode,

*l*is the azimuthal order,

*n*(

*r*) is the index at radial coordinate

*r*,

*n*is the effective index (propagation constant) of the mode, and

_{eff}*λ*is the vacuum wavelength of the guided light. In the discussion that follows, we assume that the radial index profile varies in discreet steps, or layers.

Define the dimensionless and scaled variables:

*n*is the refractive index of the layer or layers in which the field will ultimately be flattened (in the method prescribed in this paper,

_{flat}*n*is chosen before the waveguide is designed). In these terms, the wave equation becomes:

_{flat}For weak waveguides, the field and its radial derivative are continuous across the step-like boundaries between layers. Since the radial derivative is continuous, so is the quantity:

To determine the field distribution of the modes of a complex waveguide, we track *ψ* and *ζ*; we begin by determining analytic solutions for the field in layers whose index is greater than, less than, and equal to the propagation constant. Each analytic solution has two unknown constants, which can be determined by the boundary conditions.

Begin by considering layers that are neither the inner-most layer, here referred to as the “core,” nor the outermost layer, referred to as the “cladding.” The cladding is presumed to extend to infinity.

In layers where *η* > *η _{eff}* (

*n*>

*n*), the solution to the wave equation is:

_{eff}*J*and

_{l}*Y*are oscillatory Bessel functions,

_{l}*A*and

*B*are unknown constants, and:

If *ψ* and *ζ* are known at some position *v*_{1}, such as at one of the layer’s boundaries, then *A* and *B* can be expressed:

*A*and

*B*were determined with the help of the identity [12]:

*η*<

*η*(

_{eff}*n*<

*n*) the solution to the wave equation is:

_{eff}*I*and

_{l}*K*are exponentially growing and decaying modified Bessel functions and

_{l}*A*and

*B*are unknown constants. If

*ψ*and

*ζ*are known at some position

*v*

_{1}, such as at one of the layer’s boundaries, then

*A*and

*B*can be expressed:

*A*and

*B*we used the identity:

*η*=

*η*(

_{eff}*n*=

*n*) the wave equation reduces to:

_{eff}*l*≠ 0 the solution is:

*A*and

*B*become:

*l*= 0 the solution is:

*A*and

*B*become:

*B*to zero (from Eq. (34), this is equivalent to making the field’s slope zero); thus a necessary condition is that

*n*=

*n*. Comparing Eq. (29) and Eq. (32) we see that the field can only be flattened if, in addition to

_{eff}*n*=

*n*, the azimuthal order,

_{eff}*l*, is also zero.

Now consider the inner-most layer, the core, and the outer-most layer, the cladding. In these, only a single Bessel solution is allowed. In the core the solutions are:

## Appendix II: Transfer matrices

The solutions for the constants *A* and *B* can be substituted into the original expressions for *ψ* and the corresponding expressions for *ζ* to obtain transfer matrices, **M**, that relate *ψ* and *ζ* at position *v*_{2} to their known values at position *v*_{1}:

*x*

_{1}is the quantity

*x*, defined by Eq. (16), evaluated at position v

_{1}and index

*η*

_{12}(the index between v

_{1}and v

_{2}), and

*x*

_{2}is

*x*evaluated at v

_{2}and index

*η*

_{12}.

The determinant of each matrix is unity, but they are not orthogonal. Their inverses are found by exchanging their diagonal elements and changing the signs of their off-diagonal elements.

In layers where *η* > *η _{eff}* (

*n*>

*n*):

_{eff}*η*<

*η*:

_{eff}*η*=

*η*and

_{eff}*l*≠ 0:

*η*=

*η*and

_{eff}*l*= 0:

*is (from Eq. (39)):*

_{clad}*x*is the term

_{clad}*x*, as defined by Eq. (16), evaluated at position v

*and index*

_{clad}*η*= 0. Note that the Bessel derivates can be calculated from Eq. (27). Ω

_{clad}*is similarly calculated from Eq. (35), Eq. (36), Eq. (37), or Eq. (38) at the core’s boundary.*

_{core}The matrix **M** is the product of the matrices that represent the layers between the core and cladding; it takes advantage of the fact that *ψ* and *ζ* are continuous across layer boundaries. For a given waveguide, the propagation constant *η _{eff}* is determined iteratively – that is, by varying its value until the transfer matrix solution is satisfied.

In the above, (*const*) refers to a multiplicative constant related to the total power carried by a mode, as discussed in the following Appendix.

## Appendix III: Mode normalization

This appendix gives closed-form solutions for the mode normalization integral, and defines scaled fields.

Mode normalization involves choosing the (*const*) term of Eq. (46) to make the power carried by a mode equal to some preselected value, *P*_{0}:

*ψ*such that:

_{scaled}The integration is typically performed numerically, though with the expressions that follow, which we believe are novel, it can be calculated analytically. The solutions were obtained by integrating the above expression by parts twice and taking advantage of the fact that the bound modes’ fields satisfy the original wave equation, Eq. (13).

For *η* ≠ *η _{eff}* (

*n*≠

*n*):

_{eff}*η*=

*η*and

_{eff}*l*= 0:

*η*=

*η*and

_{eff}*l*= 1:

*η*=

*η*and

_{eff}*l*≥ 2:

These are the indefinite solutions to the integrals; the contribution from an individual layer is found by evaluating its solution (depending on its index relative to the propagation constant) at the its boundaries, and subtracting one from the other. The full integral (from zero to infinity) is found by summing the individual contributions. Note that, for any waveguide design, the right-most terms of the piece-wise integrals contributes the following series to the full integral:

*ψ*are continuous across interfaces, this reduces to $\pi {\left[{\text{v}}^{2}{\psi}^{2}\right]}_{\text{\hspace{0.17em}}0}^{\text{\hspace{0.17em}}\infty}$, which is zero for all bound modes. Thus, while the right-most terms contribute to the piece-wise integrals, they do not contribute to the full integral.

The closed form solutions can also be used to quickly calculate the group index of a mode via Eq. (60).

## Appendix IV: Size-spacing products

This appendix defines several mode size-spacing products and shows that for a given waveguide design, these are fixed. It refers to scaled terms defined in Appendix I.

Once the scaled index profile (Eq. (11)) is specified, the scaled propagation constants, Eq. (12), and the shapes of the allowed modes are completely determined, as implied by the form of the scaled wave equation, Eq. (13). To relate scaled quantities to those that can be measured in a laboratory, begin by noting that the effective mode area can be written:

*, is also fixed:*

_{eff}*is fixed, if a mode’s size is increased, its effective index necessarily approaches the cladding index. Since this holds for all modes, it follows that as a desired mode’s size is increased, the effective indices of all modes necessarily approach each other.*

_{eff}The effective index is the phase index of the mode. When evaluating pulse propagation effects, the group index, *n _{g}*, is also important. Using an integral form of the group index [13] it can be shown that:

*, it can be shown that the following quantity is also fixed for each mode of a waveguide:*

_{eff}*n*is the product of a mode’s phase and group indices. Like Θ

_{eff}n_{g}*, this is a strict invariant of a design (within the strictures of the weak-guiding approximation), but unfortunately the separations between the Θ*

_{eff,}*’s are not obvious indicators of the separations between the group indices. The following term is more transparent:*

_{eff,g}*depends on*

_{g}*A*it is not a true invariant of the guide. However, if the mode’s area is sufficiently large (usually the case for high power laser applications) then the third term can be neglected and since Θ

_{eff}*and Θ*

_{eff}_{eff,}

*are true invariants, then Θ*

_{g}*is approximately invariant.*

_{g}## Acknowledgments

This work performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344.

## References and links

**1. **J. M. Fini and S. Ramachandran, “Natural bend-distortion immunity of higher-order-mode large-mode-area fibers,” Opt. Lett. **32**(7), 748–750 (2007). [CrossRef] [PubMed]

**2. **S. Ramachandran, J. M. Fini, M. Mermelstein, J. W. Nicholson, S. Ghalmi, and M. F. Yan, “Ultra-large effective-area, higher-order mode fibers: a new strategy for high-power lasers,” Laser Photonics Rev. **2**(6), 429–448 (2008). [CrossRef]

**3. **R. H. Stolen and C. Lin, “Self-phase-modulation in silica optical fibers,” Phys. Rev. A **17**(4), 1448–1453 (1978). [CrossRef]

**4. **A. K. Ghatak, I. C. Goyal, and R. Jindal, “Design of a waveguide refractive index profile to obtain a flat modal field,” Proc. SPIE **3666**, 40–44 (1999). [CrossRef]

**5. **J. W. Dawson, R. Beach, I. Jovanovic, B. Wattellier, Z. Liao, S. A. Payne, and C. P. J. Barty, “Large flattened mode optical fiber for reduction of non-linear effects in optical fiber lasers,” Proc. SPIE **5335**, 132–139 (2004). [CrossRef]

**6. **W. Torruellas, Y. Chen, B. McIntosh, J. Farroni, K. Tankala, S. Webster, D. Hagan, M. J. Soileau, M. Messerly, and J. Dawson. “High peak power ytterbium-doped fiber amplifiers,” Proc SPIE **6102**, 61020N (2006).

**7. **B. Ward, C. Robin, and M. Culpepper, “Photonic crystal fiber designs for power scaling of single-polarization amplifiers,” Proc. SPIE **6453**, 645307, 645307-9 (2007). [CrossRef]

**8. **C. Zhao, Z. Tang, Y. Ye, L. Shen, and D. Fan, “Design guidelines and characteristics for a kind of four-layer large flattened mode fibers,” Optik **119**(15), 749–754 (2008). [CrossRef]

**9. **N. N. Elkin, A. P. Napartovich, V. N. Troshchieva, and D. V. Vysotsky, “Modeling of large flattened mode area fiber lasers,” Laser Phys. **20**(2), 304–310 (2010). [CrossRef]

**10. **D. Jain, C. Baskiotis, and J. K. Sahu, “Mode area scaling with multi-trench rod-type fibers,” Opt. Express **21**(2), 1448–1455 (2013). [CrossRef] [PubMed]

**11. **A. Yariv, *Optical Electronics, 3rd Edition*, (Holt, Rinehart and Winston, 1985).

**12. **M. Abramowitz and I. Stegun, *Handbook of Mathematical Functions* (Dover, 1972).

**13. **A. W. Snyder and J. D. Love, *Optical Waveguide Theory* p.644 (Chapman and Hall Ltd, 1983).