## Abstract

The coupling properties of multi-core fibers are analyzed using the multipole method and coupled mode theory in order to gain insights into the performance of image fibers that are commonly used in flexible endoscopes. It is explained that coherent fiber bundles with high core density are able to transport images because nonuniformity in the pixel size reduces the inter-core coupling that causes crosstalk. The wavelength, average core size and separation, and degree of core size variation determine the strength of coupling between adjacent cores, such that fibers with a smaller core size and separation at longer wavelengths require more nonuniformity in order for reliable image transmission. Guidelines are given for assessing the performance of image fibers in a particular system.

© 2007 Optical Society of America

## 1. Introduction

A fiber bundle contains thousands of high index cores in a common lower index cladding. The core sizes and separations are on the order of a few microns, thus differing from a simple bundle of step-index fibers. The index contrast is increased from that of standard single mode fiber in order to more tightly confine the light and reduce crosstalk between cores. In coherent fiber bundles, or image fibers, the input and the output are spatially correlated such that an image transmitted through the fiber can be reconstructed from either end face in an identical manner. Image fibers are used primarily for confocal and multiphoton endoscopic imaging [1–10] but have also been employed in areas such as optical interconnects [11–13], sensing [14], and optical coherence tomography [15]. In imaging systems, the resolution of the transmitted image is limited by the pixilated nature of the fiber bundle as well as by the degree of coupling between cores.

Applications for image fibers motivate the production of higher resolution or higher information density fibers with smaller cores that are more closely packed. Decreasing the fiber dimensions, however, will increase the strength of coupling or crosstalk between cores, resulting in blurred images and a lower signal to background ratio. Available image fibers manufactured by Fujikura and Sumitomo, for example, already have dimensions well below the acceptable theoretical limit for crosstalk. Numerical simulations of these fibers reveal short coupling lengths relative to the typical length of fiber utilized. That is, the strength of the coupling between adjacent cores is predicted to be such that power will begin to disperse amongst the cores after a propagation distance of much less than a meter. Yet these image degrading effects are not dramatic in practice. A closer look at these fibers reveals a significant degree of nonuniformity in the size and shape of individual pixels or cores. In this paper, we will describe how the introduction of nonuniformities disrupts coupling by creating a mode mismatch between adjacent cores. In certain fiber structures, despite high core density, even a small mismatch is sufficient to cause the cores to de-couple. We further verified this behavior using coupled mode theory, suggesting methods for assessing the performance of a particular image fiber and giving guidelines for choosing an image fiber for a particular wavelength. The authors anticipate that a greater understanding of coupling in multi-core fibers and the effects of nonuniformities on crosstalk will not only motivate the production of more reliable image fibers but also facilitate the more effective use of existing fibers.

## 2. Methods

The numerical simulations in this paper are based on two image fibers manufactured by Fujikura, distributed by Myriad Fiber Imaging Technology, Inc. The fiber specifications are given in Table 1 and, unless otherwise noted, the values are those given by the distributor. The samples are approximately one foot in length. These two fiber types were chosen because their dimensions are typical of imaging applications and of current research [3, 4, 10]. SEM imaging was performed on each fiber in order to verify their dimensions and the pixel size for the 350S fiber was corrected. The final column in Table 1 provides an estimate for the variation in pixel diameter taken from the standard deviation of elliptical fits to the cores in the SEM images using an image processing program. The nonuniformity in the pixel cores is evidenced by the SEM image of Fig. 1. The degree of nonuniformity in the core shape and size observed in these fiber appears to be typical of other commercially available fiber bundles, such as Sumitomo, and has been noted previously [9, 16–21].

Because the coupling behavior between the thousands of cores of an image fiber would be extremely complicated to model, smaller systems of two and seven cores are examined in detail in order to make informed assumptions about the larger system. The nature of the coupling between cores can be understood by calculating the power or energy distribution in each core at different points along the fiber length. In a multi-core system, the power oscillates between cores as a function of the propagation distance, z. The percentage of power transferred between cores is referred to as the coupling efficiency, related to crosstalk, and the rate of power transfer is described by the coupling or beat length. The coupling length and the coupling efficiency are used to quantify the strength of coupling and depend on the interactions between cores or on the beating of the modes of the fiber. The power can be determined from a normal mode expansion of the field solution in the fiber given a Gaussian single core input field. The method used here employs the multipole method to solve for the fiber modes and follows Ref. [22].

When solving for the modes of multi-core fibers, the mode group with the largest propagation constants contains 2N modes, where N is the number of cores in the system and the factor of two is consistent with two polarization states. The modes in this group will be referred to as the fundamental modes because the distribution of energy within the individual cores for each of the modes is approximately Gaussian with azimuthal symmetry, as characteristic of a fundamental mode. The normal mode expansion is truncated after these 2N lowest order modes since they represent the most significant contribution to the sum due to their large overlap with the Gaussian input field. Our modeling showed that ignoring the higher order modes results in negligible error because the modal amplitudes, a_{j}, become extremely small.

A two-core fiber has nondegenerate mode solutions of even and odd nature for each polarization, analogous behavior to a coupled oscillator system. The energy coupled into each of these modes at the fiber endface, represented by the coefficients a_{j}, will propagate at slightly different speeds along the fiber due to the splitting in the mode effective indices. The resulting effect on the power in each core is the “beating” shown in Fig. 2, where the power in each core, as determined from a normal mode expansion, is plotted over z for two-core fibers with dimensions taken from the two image fiber samples of Table 1. The energy distribution in the cores at different points along the z-axis is included below Fig. 2(a). The input field is an x-polarized Gaussian, E_{y} = 0, that is centered on the left core with a spot diameter (1/e^{2} intensity) equal to the radius of the core. As shown in Fig. 2, the coupling efficiency is 100% when the two cores have identical diameters, indicating that all the power coupled into the fiber oscillates back and forth between the two cores. The coupling length is defined in Fig. 2(b) as the distance after which all power, when incident initially on one of the cores, will be coupled into the second core. For the case when the cores are identical, the coupling length for each polarization can be calculated directly from the beat frequency of the modes in the following manner:

where e and o indicate the even and odd modes. This distance is 0.25 m for fiber 500N and 1.4 cm for fiber 350S. Because practical applications with fiber bundles use lengths much longer than these distances, coupling between two cores should be readily observed in endoscopic systems.

In an image fiber, a single core has not just one but six nearest neighbor cores and the impact of these additional cores is to increase the rate at which power is transferred between the cores. In the plot of Fig. 3, the power in each core of a seven-core fiber is plotted versus the propagation distance when light is incident on the central core. The core dimensions correspond to the FIGH-10-350S fiber and the wavelength is 600 nm. Along the length of the fiber, the system oscillates between a state where all the power is in the central core (the initial state) and a state where one-seventh of the power is in each of the outer cores and the central core. The slightly differing behavior between cores in the outer ring, visible in Fig. 3 at the point where the power in the outer cores is labeled, is a result of the linear polarization of the input field and the six-fold symmetry of the seven-core system. Note that the coupling length, or the distance to the first minimum of the power in the input core, is reduced from that of a two-core system [compare with Fig. 2(b)], providing further evidence that severe crosstalk is predicted to greatly reduce the practical use of these fibers.

Numerical simulations obviously demonstrate a discrepancy between theoretical predictions of strong crosstalk and the experimentally demonstrated successful image transmission of the image fibers in Table 1. Considering the case of only two cores, in order for at least 90% of the propagating power to remain in the initial core, simulations predict that the fiber FIGH-10-500N could be no longer than 10 cm while a length less than 0.5 cm would be necessary for fiber FIGH-10-350S. In addition, the trend demonstrated by Fig. 3 indicates that this length will decrease when more neighboring cores are considered. Typical flexible endoscopes utilize fibers on the order of a meter long; thus, core coupling or crosstalk should produce detrimental effects on the transmitted images. Yet, experimental results with similar sized fibers [8–10] do not show the significant blurring predicted here. This obvious contradiction between demonstrated fact and numerical results indicates that the simulated system does not accurately represent the fiber used in the lab. As is evident in the images of Fig. 1, available image fibers are not composed of identical cores. Numerical simulations lead to drastically different conclusions when nonuniformity is introduced into the cross-sections of the simulated multi-core fibers [22, 23].

## 3. Numerical analysis

Although both the shape and size of the image fiber cores are quite irregular, only variations in the diameters of adjacent circular cores will be examined here. The effects of other types of nonuniformity, such as core ellipticity and irregular core separation will be briefly discussed later in the paper. In order to numerically study random variations in the core diameter, each core in a multi-core fiber is randomly assigned a diameter from a Gaussian distribution with an average value of d_{o} (taken from Table 1) and a standard deviation of δd. The degree of nonuniformity in the fiber is quantified by the ratio of the standard deviation to the mean of this distribution (δd/d_{o}) times 100, referred to as a percentage of variation. If the average value used for the diameter is 2.0 μm, a variation of 1% indicates a standard deviation of 20 nm. The core separation is the average value from Table 1.

Returning first to the simple system of only two cores, asymmetry resulting from random variations in the core diameters will alter the coupling properties of a fiber differently depending on the average size of the cores, the wavelength, and the degree of nonuniformity. In general, both the coupling efficiency and the coupling length decrease in magnitude from their value when the cores are identical. Figure 4 shows an example of how the power in each core is altered from that of Fig. 2(b) when the core diameters of a FIGH-10-350S two-core fiber are no longer equal. Both types of fibers were analyzed at three different wavelengths, 600 nm, 980 nm and 1300 nm (material dispersion was neglected), for several degrees of variation. Because 100% power transfer may not occur, the coupling length is defined as the distance to the first power minimum of the incident core and the efficiency is the maximum power transferred out of the incident core. In Fig. 5, the average coupling efficiency is shown to decrease as the degree of variation increases and as the wavelength decreases for both types of fibers at three different wavelengths. Each marker represents the average efficiency for a set of 99 two-core fibers whose core diameters were randomly assigned from a Gaussian distribution of values. Notice the FIGH-10-500N image fiber is more sensitive to nonuniformity; smaller variations produce a more drastic reduction in the coupling efficiency for the FIGH-10-500N fiber as compared to the FIGH-10-350S type. When the wavelength is 600 nm, a variation of 1% in a FIGH-10-500N fiber results in an average efficiency of approximately 2e-3, indicating essentially independent core propagation or greatly reduced crosstalk. A variation of at least 10% would be needed in the FIGH-10-350S type fiber in order to obtain this same condition. If the coupling efficiency is small, very little power oscillates between the cores of a fiber and the coupling length ceases to be a relevant estimate for crosstalk or the strength of coupling. The introduction of nonuniformities in core size can therefore reduce the crosstalk in cores that would otherwise be strongly coupled. The conclusions made from studying two-core fibers are applicable to multi-core fibers because they correctly predict the response of larger seven-core systems when the core diameters are no longer equal. For a situation where the diameters of seven cores have been randomly selected with a percentage variation of 2% and one outer core happens to have a diameter nearly identical to that of the central core, the power evolution in all seven cores is shown in Fig. 6(a). The potentially complicated behavior of a seven-core system is essentially reduced to that of a two-core system. The two cores with almost identical diameters—the difference is less than a tenth of a nanometer—transfer power in a manner similar to the two-core system of Fig. 2, while the remaining cores exhibit behavior similar to Fig 4, participating very little in the power exchange. Additional simulations reveal that when the diameter of a second core in the outer ring is forced to match that of the already participating outer core, the system behaves like the respective three-core system, depending on the orientation of the three similar cores, either in a line, Fig. 6(b), or a triangle, Fig. 6(c).

The average coupling behavior of a seven-core system with randomly varying core diameters is similar to that of the two-core fibers in Fig. 5. When the diameters of the cores are no longer identical, the maximum amount of power remaining in the central core is, in general, increased. In Fig. 7, the average coupling efficiency of the central core is plotted versus the percentage variation for type FIGH-10-350S seven-core fibers at λ = 600 nm. Notice that the average efficiency in Fig. 7, for each degree of variation shown, is approximately six times the corresponding value in Fig. 5(a) for the same wavelength. Figure 5(a), therefore, provides, in the weak coupling regime, an average efficiency for each two-core interaction in a many core system. When a single core is surrounded by six cores, and on average the interactions with each of these cores have a coupling efficiency of 10%, then approximately 60% of the power will oscillate back and forth between the center core and the outer cores while the remaining 40% will be retained in the central core. Obviously, this argument is no longer relevant when the coupling is strong; for fibers with strongly coupled cores, the coupling behavior is similar to that of the symmetric structure in Fig. 3. Thus, as in the two-core system, a large percentage of variation will increase the power that remains in the input core; however, this amount will be decreased by the interactions with each neighboring core.

A more practical way to view the data in Figs. 5 and 7 is as a function of the diameter difference between the input core and the other cores in the fiber. In this way, the coupling properties of two cores can be related directly to their size mismatch. These relationships are plotted in Fig. 8 for the same fibers used for the data in Fig. 5. Each marker represents a different randomly generated fiber, including all the degrees of variation. As the difference between the diameter of the input core and that of the other core, or Δd, increases, both the coupling length and the efficiency decrease. This decrease occurs more or less quickly depending on the wavelength and the average core size. For example, the efficiency decreases more dramatically at shorter wavelengths and for fiber type FIGH-10-500N, indicating that sensitivity to nonuniformity increases as the wavelength decreases and as the average core size increases. In addition, the coupling length is nearly wavelength independent when the mismatch is large. The spread in the data is due to the fact that the diameters of both cores are varied randomly about an average value; therefore, each Δd does not refer to a unique structure and can be associated with a different coupling length and coupling efficiency.

Similar trends appear when the efficiency and the coupling length of each of the outer cores of a seven-core fiber are plotted versus their difference in diameter from the central core, as shown in Fig. 9 for the 350S fiber at a wavelength of 600 nm. For these fibers, the central core is kept at d_{0}, the average value, and the diameters of the outer cores are randomly assigned according to a percentage variation of up to 4%, thus reducing the spread in the data points from the two-core data. The close resemblance of Figs. 8(a) and 8(b) to Figs. 9(a) and 9(b) indicate that the behavior of a seven-core fiber can be predicted by studying two-core fibers. Provided an estimate can be made for Δd in a particular fiber, understanding these relationships could lead to a method for predicting the coupling behavior of image fibers with variations.

## 4. Coupled mode theory

Coupled mode theory (CMT) is commonly used for obtaining approximate analytical solutions to systems of coupled waveguides, such as multi-core fibers. In CMT, each core is solved for independently and interactions with neighboring cores are treated as a perturbation. The total field solution of a system of two waveguides, label them *a* and *b*, is approximated as a linear combination of the mode fields of the individual waveguides [24, 25], as shown in Eq. (2).

$$\overrightarrow{H}\left(x,y,z\right)=a\left(z\right){\overrightarrow{h}}^{a}\left(x,y\right)+b\left(z\right){\overrightarrow{h}}^{b}\left(x,y\right)$$

The power as a function of z is calculated from Eq. (2) in the standard way as shown in Eq. (3).

$${C}_{\mathit{ab}}=\frac{1}{2}\int \int {\overrightarrow{e}}^{a}\left(x,y\right)\times {\overrightarrow{h}}^{*b}\left(x,y\right)\xb7\hat{z}\mathit{dxdy}$$

The amplitude coefficients of the modes of the two cores obey the coupled amplitude equations of Eq. (4) as derived from the general reciprocity relation [24, 26].

$$\frac{d}{\mathit{dz}}\mathit{b}\left(z\right)=i{K}_{\mathit{ba}}a\left(z\right)+{\mathit{i\beta}}_{\mathit{b}}b\left(z\right)$$

In the weak coupling regime, or conventional CMT, the cross field power (the terms in Eq. (3) containing the cross field integrals, *C _{ab}* and

*C*) is ignored. Although this is strictly only true when

_{ba}*K*=

_{ab}*K*

_{ba}^{*}and the cores are identical, this approximation is valid when

*C*and

_{ab}*C*are very small and

_{ba}*K*≈

_{ab}*K*

_{ba}^{*}. In this case,

*β*and

_{a}*β*are the propagation constants for the individual, independent waveguides and the coupling coefficients,

_{b}*K*and

_{ba}*K*, are proportional to the overlap integral of the mode fields of the individual waveguides in each core [24, 27]. Improved CMT considers the situation where

_{ab}*K*and

_{ab}*K*are very different, therefore second order terms are added to the constants in Eq. (4) and the cross field power is included in order to conserve power [24, 27]. There has been much discussion over the last two decades concerning the limitations of conventional CMT and the regimes for which ICMT is necessary [27–29]. At this time, the authors continue to ignore the cross field power as the simplicity of conventional CMT has its advantages. The validity of these assumptions will be assessed later in the paper.

_{ba}Solving the system of equations in Eq. (4) and applying the initial condition that light is incident on waveguide *a* at z = 0, or *a*(0) = 1 and *b*(0) = 0, produces the expressions for the power given in Eq. (5) [27]. The notation used here is that of Snyder [27] and differs from that used by the authors in previous work [22] because the definitions of the constants in Eq. (5) are better suited to the arguments presented in this paper.

$${P}_{b}\left(z\right)={\mid b\left(z\right)\mid}^{2}=\frac{{\mathit{K}}_{\mathit{ab}}}{{K}_{\mathit{ba}}}{F}^{2}{\mathrm{sin}}^{2}{\beta}_{d}z$$

$$F=\frac{1}{\sqrt{1+\frac{{\left({\beta}_{a}-{\beta}_{b}\right)}^{2}}{4{\mathit{K}}_{\mathit{ab}}{\mathit{K}}_{\mathit{ba}}}}}$$

$${\beta}_{\mathit{d}}=\frac{\sqrt{{K}_{\mathit{ab}}{\mathit{K}}_{\mathit{ba}}}}{F}$$

The maximum power transferred out of waveguide *a*, or the coupling efficiency, is *F*
^{2} and *β _{d}* is related to the beat frequency or the rate of power transfer. The coupling length is the distance, L

_{c}, at which the power in waveguide

*b*has oscillated to its first maximum, ie.

*β*/

_{d}*L*=π/2 . As the mismatch between the modes of the individual cores,

_{c}*β*–

_{a}*β*or

_{b}*Δβ*, increases, these relationships also dictate that the coupling efficiency and the coupling length will decrease, as observed previously.

The expressions in Eq. (5) can now be applied to the data presented in Fig. 8. Specifically, the trend in the coupling efficiency in Figs. 8(a) and 8(c) should be described by *F*
^{2}, or

and the coupling length is π/(2*β _{d}*) or

The mismatch, or *Δβ*, is then considered to be a function of *Δd* and the product *K _{ab}K_{ba}* is treated as a constant which is typically a valid assumption [30]. Because

*Δβ*is determined from the propagation constants of the modes of the individual waveguides, this value can be approximated by calculating

*β*over a range of core diameters for a standard step index fiber [31, 32]. For the degree of variation considered here,

*β*is an approximately linear function of diameter, and

*Δβ*can be written as the slope times

*Δd*. Equations (6) and (7) are then rewritten as the following:

$${L}_{c}=\frac{b}{\sqrt{1+a{\left(\Delta d\right)}^{2}}}$$

where the constant *a* is equal to slope^{2}/(4*K _{ab}*

*K*) and

_{ba}*b*is the coupling length for the system when the two cores are identical. The fitting of Eq. (8) to the data in Fig. 8 is shown for fiber FIGH-10-500N in Fig. 10. Table 2 displays the product

*K*

_{ab}*K*as solved for from the fit parameter

_{ba}*a*for the four plots in Fig. 8. The values for

*K*

_{ab}*K*derived from the two fit parameters in the coupling length plots are approximately equal; therefore, only

_{ba}*K*

_{ab}*K*as derived from the

_{ba}*a*parameter are shown in Table 2 for Figs. 8(b) and 8(d). The final column of Table 2 gives the value for

*K*

_{ab}*K*that is calculated from a two-core system with identical cores by setting Eq. (1) equal to Eq. (7) for Δβ = 0 and solving for

_{ba}*K*

_{ab}*K*.

_{ba}The values for *K _{ab}*

*K*in Table 2 are very similar demonstrating that simple conventional CMT is adequate for predicting the coupling behavior for image fibers with core diameter mismatch and that the product

_{ba}*K*

_{ab}*K*is essentially constant for the degree of nonuniformity studied here. In addition, these fits are consistent for the long wavelength data as well, indicating that the expressions in Eq. (6) and Eq. (7), derived using conventional CMT, retain accuracy even in a stronger coupling regime. Therefore, these equations can be used for predicting the coupling behavior of nonuniform two-core fibers, which we have shown are a good model for multi-core fibers such as fiber bundles.

_{ba}The quality of images transmitted by an image fiber depends on the amount of crosstalk between cores, or the coupling efficiency; however, because crosstalk can be reduced through the introduction of mismatch between adjacent cores, the performance of a particular fiber depends on the amount of nonuniformity, the wavelength, and the average core size. In order to better understand the interplay of these different parameters and the sensitivity of a particular fiber to core diameter mismatch we compare the two parameters that determine the efficiency, *K _{ab}*

*K*and

_{ba}*μβ*

^{2}/4, see Eq. (6). The relative magnitude of these values determines how sensitive a system is to asymmetry, and also the average amount of crosstalk in a fiber. Because

*K*

_{ab}*K*is essentially constant over variations of approximately 10%, it is therefore a property of the symmetric system and can be calculated reliably from Eq. (1) for any two-core fiber. The parameter

_{ba}*Δβ*can be estimated from the degree of nonuniformity in a particular image fiber. If a particular two-core fiber is weakly coupled with a small value for

*K*

_{ab}*K*, a small

_{ba}*Δβ*caused from nonuniformities can significantly decrease the coupling efficiency; this fiber would be considered as very sensitive. In >Fig. 11, these two parameters are plotted versus wavelength for the fiber type FIGH-10-350S in Fig. 11(a), and type FIGH-10-500N in Fig. 11(b). The solid lines indicate

*K*

_{ab}*K*, as calculated from Eq. (1). The non-solid lines are

_{ba}*Δβ*

^{2}/4 for different degrees of asymmetry, where Δ

*β*=

*β*(

^{b}*λ*)-

*λ*(

^{a}*λ*) . For simplicity, the core labeled,

*a*, maintained a diameter of d

_{0}, the average core diameter given in Table 1, while the diameter of core

*b*was increased by 1%, 4% and 10%. The propagation constant,

*β*, was calculated over wavelength for each core separately and subtracted in order to obtain

*Δβ*.

The utility of these plots can be understood by examining, for example, the point where a solid line crosses a non-sold line. At the wavelength of this crossing point, the predicted efficiency, from Eq. (6), will be 50%. Above this wavelength, the efficiency will increase and at lower wavelengths it will decrease. As the wavelength and *K _{ab}*

*K*decrease, a fiber becomes more sensitive to nonuniformities, approaching a region of independent core propagation. For example, when the wavelength and value for

_{ba}*Δβ*fall in the gray shaded region, the two cores will have less than 1.67% coupling efficiency. For the opposite trend in wavelength and

*K*

_{ab}*K*, a fiber enters a strongly coupled regime where reducing crosstalk becomes more and more difficult. Notice also that fiber type FIGH-10-500N has a larger wavelength range with low efficiency coupling demonstrating again that this fiber is more sensitive to nonuniformity. Due to the scale invariance of Maxwell’s equations, the results in Fig. 11 can be applied to situations not within the range shown and to fibers of different core sizes, if the wavelength and fiber dimensions are scaled appropriately. These plots can also be generated for any two-core fiber given a step-index fiber mode solver and the average coupling length when the cores are identical.

_{ba}## 5. Discussion

The dramatic impact of nonuniformity on the coupling properties of multi-core fibers can be understood more intuitively by examining the changes that occur in the modal fields when variation in core size is present. When the core diameters of a seven-core fiber differ significantly, the modes of the multi-core system decouple into those of the individual cores. The existence of localized modes in a disordered waveguide system such as a fiber bundle has also been observed by others [33]. Figure 12 shows the seven unique energy distributions of the fourteen nondegenerate modes of a decoupled seven-core system, which can be compared to the modes of the identical core fiber shown in Fig. 3. Light incident on a single core of the fiber in Fig. 12 will couple almost completely into just one mode of the system and therefore not mix with other cores or modes along the length of the fiber. The field expansion will contain essentially one mode and the energy distribution is no longer a strong function of z. The cores of this image fiber behave independently of their neighbors because differing adjacent core diameters have created a mode mismatch that will inhibit coupling.

In addition to variations in size, the cores of image fibers also vary in shape and separation. The impact of these nonuniformities can be estimated to be less than that of variation in core size. The increase in *Δβ* caused by core ellipticity can be estimated from studies on the geometric birefringence of elliptical core fibers [34–36]. For example, if the variation in core diameter is 10%, a maximum birefringence can be estimated by defining the major and minor axes of an elliptical core as 10% more and 10% less than the average core diameter, resulting in an ellipticity of approximately 0.182. This shape distortion would cause a *Δβ* on the order of 10^{-4} [34–36]. Comparing the square of this value to the non-solid lines in Fig. 11, it is evident that ignoring the effect of nonuniformity in core shape is appropriate due to the relatively small size of this parameter compared to the *Δβ* caused by variation in the core diameters. Ellipticity results in a minimal amount of mode mismatch because the overall core area is not altered as dramatically as when the core radius is changed; the propagation constant is effected more by a change in the modal area than by a change in the modal shape. Irregularity in the core spacing will affect the coupling length rather than the coupling efficiency when the core diameters are the same. For non-identical cores, changing the separation will make the two-core system more or less sensitive to core diameter mismatch, depending on whether the cores become closer or farther apart. Because an increase in the core separation has the opposite consequence from decreasing the core separation, the net effect is expected to be small. Core diameter mismatch is, therefore, the dominant effect in these multi-core fibers.

Our numerical and analytical assessments of multi-core fibers provide a method for analyzing the crosstalk in a particular fiber bundle. The plots of Fig. 11, based on the parameters in Eq. (6), can be used to estimate the coupling efficiency of two core interactions given a certain degree of nonuniformity. Based on arguments present in this paper, this twocore analysis can be applied to cores with six neighboring cores. In this manner, the reliability of a pixel in an image fiber can be estimated give a certain acceptable amount of crosstalk and an average degree of nonuniformity. For example, if the pixel size of the FIGH-10-350S fiber type is preferred for an experiment, yet a coupling efficiency of no more than 0.1 can be tolerated at λ = 780 nm—in other words, 90% of the signal power is expected to be retained in each core—the average coupling efficiency for each core pair should then be approximately 0.0167, when six near neighbor interactions are considered. The value for *Δβ* that would produce this efficiency would need to be at least sixty times larger than *K _{ab}*

*K*(using Eq. (6) ); the shaded region in Fig. 11 indicates the values for

_{ba}*Δβ*that satisfy this condition. Referring to Fig. 11, a diameter difference of at least 4% would result in a difference of this order of magnitude. Since this type of image fiber has approximately a 10% variation in core diameter, as measured by SEM (Table 1), and given a Gaussian distribution, there is a 70% chance that two cores will have a diameter difference of exactly 4% or higher. Therefore, approximately 30% of the pixels will be limited by crosstalk given this situation and tolerance. Due to the higher sensitively of the FIGH-10-500N fiber, better performance would be anticipated given the same efficiency restrictions. This example illustrates the fact that the coupling will typically not be reduced below a desired limit in one hundred percent of the cores in an image fiber; some degree of error is inevitable. Table 3 provides an estimate of the pixel accuracy of the two fibers studied given two different values for the total percentage of power coupled out of the input core that can be tolerated in an experiment. The estimates also assume the fiber exhibits 10% variation in the core size and that the wavelength is 780 nm. A percentage of 95% indicates that around 500 of the total ten thousand cores will exhibit stronger coupling than desired. In fact, strong core coupling in these image fibers has been demonstrated experimentally and is described elsewhere [37]. These effects will be manifested by an overall reduction in the image resolution and the reliability of each pixel value.

## 6. Conclusion

Theory and simulation indicate that due to the small core size and separation of current image fibers, crosstalk between pixels should make the transmission of images through flexible endoscopes severely blurred. However, nonuniformity in the core size is shown to reduce the efficiency of the coupling between adjacent cores, substantially suppressing this crosstalk. The manufacture of functional image fibers with closer and smaller pixels requires the pixels to be nonuniform in order to enhance mismatch between neighboring cores. The degree to which two cores will interact or couple depends on the wavelength, the degree of variation, and the average core size. Coupled mode theory was used to generate expressions that can predict the amount of nonuniformity that will reduce crosstalk to acceptable levels in addition to assessing the pixel reliability of available image fibers.

## Acknowledgments

The authors would like to acknowledge the Cornell Nanofabrication Facility (CNF) staff for assisting with SEM images, Jim McDonald at Myriad Fiber Imaging Tech., Inc., for the fiber samples and information about the fibers, Larry Chen and other members of the Xu group for support and assistance, and finally the University of Sydney through its ARC Centre of Excellence - Ultrahigh-bandwidth Devices for Optical Systems, School of Physics as the exclusive Licensor of the CUDOS MOF UTILITIES Software used by the authors as a mode solver for multi-core fibers.

## References and links

**1. **A. F. Gmitro and D. Aziz, “Confocal microscopy through a fiber-optic imaging bundle,” Opt. Let. **18**, 565 (1993). [CrossRef]

**2. **A. F. Gmitro, A. R. Rouse, and A. Kano, “In vivo fluorescence confocal microendoscopy,” in *Biomedical Imaging, 2002. Proceedings 2002 IEEE International Symposium on*, pp. 277–280, (2002).

**3. **Y. S. Sabharwal, A. R. Rouse, L. Donaldson, M. F. Hopkins, and A. F. Gmitro, “Slit-scanning Confocal Microendoscope for high resolution in Vivo Imaging,” Appl. Opt. **38**, 7133–7144 (1999). [CrossRef]

**4. **J. Knittel, L. Schnieder, G. Buess, B. Messerschmidt, and T. Possner, “Endoscope-compatible confocal microscope using a gradient index-lens system,” Opt. Commun. **188**, 267–273 (2001). [CrossRef]

**5. **K.-B. Sung, C. Liang, M. Descour, T. Collier, M. Follen, and R. Richards-Kortum, “Fiber-optic confocal reflectance microscope with miniature objective for in vivo imaging of human tissues,” Biomedical Engineering, IEEE Transactions on **49**, 1168–1172 (2002). [CrossRef]

**6. **V. Dubaj, A. Mazzolini, A. Wood, and M. Harris, “Optic fibre bundle contact imaging probe employing a laser scanning confocal microscope,” J. Microsc. **207**, 108–117 (2002). [CrossRef] [PubMed]

**7. **W. Göbel, J. N. D. Kerr, A. Nimmerjahn, and F. Helmchen, “Miniaturized two-photon microscope based on a flexible coherent fiber bundle and a gradient-index lens objective,” Opt. Lett. **29**, 2521–2523 (2004). [CrossRef] [PubMed]

**8. **R. Juskaitis, T. Wilson, and T. F. Watson, “Real-time white light reflection confocal microscopy using a fibre-optic bundle,” Scanning **19**, 15–19 (1997). [CrossRef]

**9. **G. L. Goualher, A. Perchant, M. Genet, C. Cavé, B. Viellerobe, F. Berier, B. Abrat, and N. Ayache, “Towards optical biopsies with an integrated fibered Confocal Fluorescence Microscope,” Lecture notes in Computer Science **3217**, 761–768 (2004). [CrossRef]

**10. **E. Laemmel, M. Genet, G. Le Goualher, A. Perchant, J-F. Le Gargasson, and E. Vicaut, “Fibered Confocal Fluorescence Microscopy (Cell-viZio) facilitates extended imaging in the field of Microcirculation,” J. Vasc. Res. **41**, 400–411 (2004). [CrossRef] [PubMed]

**11. **Y. Li, T. Wang, H. Kosaka, and S. K. Kasahara, “Fiber-image-guide-based bit-parallel optical interconnects,” Appl. Opt. **35**, 6920–6933 (1996). [CrossRef] [PubMed]

**12. **H. Kosaka, M. Kajita, Y. Li, and Y. Sugimoto, “A two-dimensional optical parallel transmission using a vertical-cavity surface-emitting laser array module and an image fiber,” Photon. Technol. Lett. **9**, 253–255 (1997). [CrossRef]

**13. **R. K. Kostuk and J. Carriere, “Interconnect characteristics of fiber image guides,” Appl. Opt. **40**, 2428–2434 (2001). [CrossRef]

**14. **C. Amatore, A. Chovin, P. Garrigue, L. Servant, N. Sojic, S. Szunerits, and L. Thouin, “Remote fluorescence imaging of dynamic concentration profiles with micrometer resolution using a coherent optical fiber bundle,” Anal. Chem. **76**, 7202–7210 (2004). [CrossRef] [PubMed]

**15. **T. Xie, D. Mukai, S. Guo, M. Brenner, and Z. Chen, “Fiber-optic-bundle-based optical coherence tomography,” Opt. Lett. **30**, 1803–1805 (2005). [CrossRef] [PubMed]

**16. **A. Komiyama and M. Hashimoto, “A new class of crosstalk in image fibers,” Opt. Commun. **107**, 49–53 (1994). [CrossRef]

**17. **M. Nakamura, T. Otsubo, and K. Kitayama, “Skew characteristics of image fiber for high-speed 2-D parallel optical data link,” J. Lightwave Technol. **18**, 1214–1219 (2000). [CrossRef]

**18. **R. Conde, C. Depeursinge, B. Gisin, N. Gisin, and B. Groebli, “Refractive index profile and geometry measurements in multicore fibres,” Pure Appl. Opt. **5**, 269–274 (1996). [CrossRef]

**19. **K. Nishioka, K. Ono, and M. Shiraiwa, “Image Fiber,” in www.freepatentsonline.com (Olympus Optical Co., Ltd., Tokyo, Japan, USA, 1995).

**20. **K. Ono, M. Shiraiwa, and K. Nishioka, “Image fiber and method of fabricating the same,” in www.freepatentsonline.com (Olympus Optical Co., Ltd., Tokyo, Japan, USA, 2000).

**21. **S. Kumar, U. H. Manyam, and V. Srikant, “Optical fibers having cores with different propagation constants, and methods of manufacturing same,” in www.freepatentsonline.com (Corning Incorporated, Corning, NY, USA, 2003).

**22. **K. L. Reichenbach and C. Xu, “Independent core propagation in two-core photonic crystal fibers resulting from structural nonuniformities,” Opt. Express **13**, 10336–10348 (2005). [CrossRef] [PubMed]

**23. **K. L. Reichenbach and C. Xu, “The effects of randomly occurring nonuniformities on propagation in photonic crystal fibers,” Opt. Express **13**, 2799–2807 (2005). [CrossRef] [PubMed]

**24. **S. L. Chuang, *Physics of Optoelectronics Devices* (Wiley-Interscience, 1995).

**25. **A. W. Snyder, “Coupled-Mode Theory for Optical Fibers,” J. Opt. Soc. Am. A **62**, 1267 (1972). [CrossRef]

**26. **S.-L. Chuang, “A coupled-mode theory for multiwaveguide systems satisfying the reciprocity theorem and power conservation,” J. Lightwave Technol. **5**, 174–183 (1987). [CrossRef]

**27. **A. W. Snyder and A. Ankiewicz, “Optical fiber couplers-optimum solution for unequal cores,” J. Lightwave Technol. **6**, 463–474 (1988). [CrossRef]

**28. **H. A. Haus, W.-P. Huang, and A. W. Snyder, “Coupled-mode formulations,” Opt. Lett. **14**, 1222 (1989). [CrossRef] [PubMed]

**29. **W.-P. Huang, “Coupled-mode theory for couple optical wavedguides: an overview,” J. Opt. Soc. Am A **11**, 963- (1994). [CrossRef]

**30. **E. Marcatili, “Improved coupled-mode equations for dielectric guides,” IEEE J. Quantum Electron **22**, 988–993 (1986). [CrossRef]

**31. **A. Snyder and J. Love, *Optical Waveguide Theory* (Kluwer, London, 1983).

**32. **A. Yariv, *Optical Electronics in Modern Communications* (Oxford University Press, Inc., New York, 1997).

**33. **A. Komiyama, “Localization of mode waves in a disordered multi-waveguide system,” Opt. Commun. **151**, 25–30 (1998). [CrossRef]

**34. **A. Kumar, R. K. Varshney, and K. Thyagarajan, “Birefringence calculations in elliptical-core optical fibres,” Electron. Lett. **20**, 112–113 (1984). [CrossRef]

**35. **R. Dyott, J. Cozens, and D. Morris, “Preservation of polarization in optical-fiber waveguides with elliptical cores,” Electron. Lett. **15**, 380–382 (1979). [CrossRef]

**36. **K. Okamoto, T. Hosaka, and Y. Sasaki, “Linearly single polarization fibers with zero polarization mode dispersion,” IEEE J. Quantum Electron. **18**, 496–503 (1982). [CrossRef]

**37. **K. L. Reichenbach and C. Xu have submitted for publication a manuscript entitled “Analysis and measurement of light propagation in coherent fiber bundles”.