Materials, Photonic Devices, and Sensors

Trapping efficiency of fluorescent optical fibers

[+] Author Affiliations
Jonathan D. Weiss

JSA Photonics, LLC, 439 Tayne Lane, P.O. Box 2930, Corrales, New Mexico 87048, United States

Opt. Eng. 54(2), 027101 (Feb 17, 2015). doi:10.1117/1.OE.54.2.027101
History: Received November 21, 2014; Accepted January 27, 2015
Text Size: A A A

Open Access Open Access

Abstract.  Fluorescent optical fibers are used in nuclear detection and other forms of fiber-optic sensors. The trapping efficiency of a fluorescent optical fiber is defined by the optical energy trapped (or guided) by the fiber divided by the total energy emitted within it by the fluorescers that dope the fiber core. This characteristic is clearly important in determining the size of signals from these devices. A calculation of the trapping efficiency has been performed under the assumption that the fluorescence radiation is emitted isotropically by the individual fluorescers that are uniformly distributed throughout the core and are equally likely to be excited by particles or shorter-wavelength light. At the price of increased complexity, nothing in the analysis precludes the lifting of these restrictions. What is included in this analysis is the contribution of skew rays, which, to the author’s knowledge, is not presented elsewhere. A very simple expression for the trapping efficiency as a function of the cladding-to-core index ratio is derived. Also important in determining signal size is the transmission loss of the fluorescence radiation to either end of the fiber from the point of its generation. However, as it is a separate matter, it is not discussed here.

Figures in this Article

Plastic fluorescent optical fibers have been widely discussed in detail as nuclear particle detectors as a result of the scintillation radiation they produce in their interior when bombarded with such particles. Some of the light travels to either end of the fiber, where it generates an electrical signal in a photodetector. Four examples of such discussions are in Refs. 1234, along with a textbook on radiation detection.5 Fluorescent fibers are also used in at least one optical sensor that is unrelated to particle detection.6 In this case, the fluorescence at a given wavelength is generated by a pump light of a shorter wavelength operating continuously. Saint Gobain7 and Nanoptics8 are two companies that manufacture fluorescent fibers.

The fraction of the fluorescence power generated within the fiber that is trapped or guided by it is the subject of this paper. Its practical importance lies in its contribution to the size of the signal in the above applications, along with the contribution from the attenuation of the guided power between the points of generation and detection. Scattering or absorption of the fluorescence radiation by the fluorescers themselves adds to the attenuation caused by the bulk material of the fiber. Although the wavelength- and fluorescer-dependent loss of the signal during propagation will not be considered here, as it is a completely separate matter, measurements of transmission loss are presented in Refs. 3, 4, and 7.

Before presenting the calculation of the trapped power, we state our assumptions that form the basis for that calculation. The first is that the fluorescers are uniformly distributed throughout the fiber core and that the distribution can be characterized by a volume density, ρ, of fluorescers. The second is that the fluorescers emit radiation isotropically into 4π steradians. The third is that they are equally likely to emit, which means that the stimulus is uniform throughout the fiber cross-section, at least on a time-average basis. The stimulus, I, can be thought of as an optical power per unit area or the number of particles crossing a unit area per unit time. If the probability of emission per unit stimulus and per unit solid angle is p, then the probability of emission per unit solid angle is pI. The last is that the fibers are highly multimode. In other words, they have a large enough diameter and numerical aperture that a ray-optics treatment is completely adequate.

The first and fourth assumptions are eminently reasonable. The second assumption could be a bit of an oversimplification because the fluorescers may emit dipole radiation, which is not isotropic. On the other hand, if they are randomly oriented in a small macroscopic volume, their angular average could be close to isotropic. To determine the details of emission, one would have to examine how fluorescent dye molecules are incorporated into a polymer matrix, such as with respect to their alignment, and how those molecules respond to various stimuli, I. The third assumption is one of linearity, which is valid for small enough stimuli. Are the stimuli encountered in practice likely to be small enough for linearity to hold? Incorporating the answers to these questions into a calculation of the trapping efficiency is outside the scope of this paper and is almost certain to require a numerical approach. However, doing so may provide an illuminating comparison to the results obtained in this self-contained, analytic treatment.

Based on these assumptions, we can readily see that the total fluorescent power per unit thickness of fiber, PT, is Display Formula

PT=(πa2)(4π)(ρ)(pI).(1)

In this equation, a is the fiber radius. When considering the trapping efficiency, the explanation provided by, for example, Saint Gobain, is as shown in Fig. 1.

Graphic Jump Location
Fig. 1
F1 :

Emission of guided radiation from an on-axis fluorescer within a fluorescent optical fiber.

A particle excites a fluorescer located along the fiber axis that consequently emits over an angle θ, such that the radiation is guided. The condition for guidance or total internal reflection (TIR) is that cosθ>ncl/nco, where the two refractive indices refer to the core and cladding (1.6 and 1.49, for the Saint Gobain fiber). Then the fraction of the total power emitted by these axial fluorescers (Pf) that is guided can easily be shown to be Display Formula

Pf=1cosθ=1ncl/nco,(2)
where θ is the maximum angle of emission that leads to guidance. For these core and cladding indices, Eq. (2) yields 6.88% and a maximum emission angle of 21.4deg. The percentage is the total for both forward and backward propagation of the fluorescence radiation. The result is, however, a simplification because it leaves out the off-axis sources of fluorescence. These off-axis sources produce skew rays, which corkscrew their way along the fiber while never crossing the axis, as opposed to the meridional rays that always cross the axis. The on-axis sources only produce meridional rays. Both of these are depicted in Fig. 2.

Graphic Jump Location
Fig. 2
F2 :

Skew and meridional ray paths within an optical fiber.

It is mentioned in 7 that 7% of the emission originating near the core-cladding interface will be guided, and passing consideration is given to them in 5. To the author’s knowledge, however, no quantitative analysis of the trapping efficiency has been published that includes skew rays. We now turn our attention to them.

We should first note that because of the continuous distribution of fluorescers within the fiber core and their continuous emission over angles, all but an infinitesimal fraction of the fluorescent radiation is in the form of skew rays. What is critical about skew rays in the context of trapping efficiency is the fact that their angle of incidence at the core-cladding interface depends on their azimuth, as well as their elevation (π/2θ). Their azimuthal angle is φ in Fig. 3, where, in a perspective view, a skew ray is shown emerging from a fluorescer at a distance r from the center of the fiber, along an arbitrary radius.

Graphic Jump Location
Fig. 3
F3 :

Geometry of a skew ray.

The angle α is that between the horizontal projection of the skew ray and the radius of the fiber. The angle of incidence is the slant angle β, which is larger than it would be if φ were 0 or π radians. If it were either, then the skew ray would be a meridional ray. Thus, skew rays may be incident at an angle greater than the critical angle over part, all, or none of the range of φ, depending on r and θ. This added TIR is the basis of the comment in 7 as mentioned above. Taking projections of the skew ray onto the horizontal plane and the radius, we obtain: Display Formula

cosβ=sinθcosα.(3)

Based on the law of sines, Display Formula

sinα=rsinφa.(4)

Thus, Display Formula

cosβ=sinθ1(rsinφ/a)2.(5)

It is clear from this equation that the angle of incidence is smallest when (r/a) or sinφ, or both are zero. It then increases as the product of these two factors increases, reaching 90 deg when the product is unity. This implies that if θ is small enough (θ21.4deg for Saint Gobain fibers) to initially allow for TIR, then TIR will exist throughout the entire range of φ and (r/a). The behavior of the angle of incidence is shown in Figs. 4 and Fig. 5, which are graphs of β versus φ.

Graphic Jump Location
Fig. 4
F4 :

Angle of incidence, β, versus azimuthal angle, φ, for various values of (r/a) and for θ=32.1deg.

Graphic Jump Location
Fig. 5
F5 :

Angle of incidence, β, versus azimuthal angle, φ, for various values of θ and for (r/a)=0.5.

In Fig. 4, the value of θ for all curves is 32.1 deg, or 1.5 times the largest angle that allows for TIR when (r/a)sinφ is zero. Thus, β< the critical angle for some range of φ. For the smaller values of (r/a), this is true over the entire range of φ. However, for (r/a)=0.75 and 0.95, respectively, TIR does occur for φ>76 and 45deg.

In Fig. 5, the value of (r/a)=0.5 for all curves, but θ is allowed to vary upward in the immediate vicinity of 21.4 deg. At 21.4 deg, TIR occurs for all φ. As θ increases, the crossover to TIR occurs at progressively larger values of φ. At θ=22.6deg, for example, TIR only occurs for φ>40deg. At θ=24.9deg, no TIR occurs for any value of φ.

Before presenting the calculation of the trapping efficiency based on the above considerations, we should note only for completeness that very technically, skew rays that would not be guided if they were meridional rays are not guided in a round fiber, even if they satisfy the condition for TIR at the point of incidence. The reason for this is because for skew rays, the boundary is curved in their plane of incidence and Fresnel’s laws of reflection are valid only for a planar boundary. The result is that although these rays are predicted to experience TIR, they lose energy by tunneling out of the core. For a radius of curvature of the boundary that is much larger than the wavelength of light, which is certainly the case here, the loss is extremely gradual with the propagation distance and is only of consequence over many kilometers of fiber length, which is certainly not the case here. It should also be noted that when these rays are predicted to refract out of the core, they do so even in the case of a curved boundary. The curved boundary is of no practical consequence here and is being ignored. The whole notion of tunneling or weakly leaky rays was developed by Snyder and coworkers in the 1970s and is discussed in 9 and the references contained therein.

As stated earlier, the trapping efficiency is the fraction of the total power emitted in 4π steradians by all fluorescers per unit length that is guided by the fiber. We break up the power that is guided into those rays that are emitted into an angle 0θθo (referred to here as narrow-angle rays) and those for which (π/2)θθo (or wide-angle rays). For Saint Gobain fibers, θo21.4deg. The contribution of the narrow-angle rays to the total guided optical power, PN, is Display Formula

PN=(πa2)(2)(2π)(ρ)(pI)(1nclnco).(6)

The numerical prefactors involve the area of the fiber, the propagation forward and backward, and the integral over φ. To obtain the fraction of the total power, we divide PN by PT in Eq. (1). The additional contribution is from wide-angle rays and for which ββcrit, the critical angle over some range of φ and (r/a). We now consider that contribution, denoted by PW. We note that Display Formula

cosβcrit=sinθo=NA/nco,(7)
where the numerical aperture of the fiber, NA, is {nco2ncl2}(1/2). Since ββcrit, cosβcosβcrit. This condition and the one on wide-angle rays, sinθo>NA/nco, form the basis for calculating the contribution to the guided power from these rays. The calculation will involve φ up to π/2, even though it extends to 2π. By symmetry, the result can be multiplied by 4. It is then multiplied by 2 because both forward and backward propagation are being included. Finally, it is multiplied by 2π to include the angular contribution to the circular area of the fiber. The prefactor is then 16πρpI. This prefactor is then multiplied by a triple integral involving the solid angle over which ray emission takes place and the range of the radial position of the fluorescer. The resulting expression for PW is Display Formula
Pw=16πρpIθLθUsinθdθφLφUdφ(ra)L(ra)U(a2)(ra)d(ra).(8)

We should note that the order of the integration over φ and (r/a) can be reversed, such that the integral over φ is first. The result is unchanged. The values of θU and φU are π/2 for reasons discussed above. The value of (r/a)U is obviously 1. The value of θL is sin1(NA/nco) because we are considering only wide-angle rays. The lower limit of (r/a), (r/a)L, is determined from Eq. (7), and β>βcrit. The result is that Display Formula

(ra)1(NAncosinθ)21sinφ=(ra)L.(9)

Since (r/a)1, Display Formula

sinφ1(NAncosinθ)2=sinφL.(10)

With these relationships, the contribution to the trapped energy from wide-angle rays is, including the prefactors, Display Formula

PW=8π2a2ρpI[(π2)(nclnco)θLθUsinθsin11(NAncosinθ)2dθ(NAnco)θLθU1(NAncosinθ)2dθ].(11)

Let us define the two integrals in Eq. (11) as I1 and I2. Then we can rewrite the equation as Display Formula

PW=8π2a2ρpI[(π2)(nclnco)I1(NAnco)I2].(12)

By combining a certain amount of ingenuity with standard integration techniques, one can prove,10 for these particular limits of integration, that both integrals in Eq. (12) are equal to Display Formula

I1=I2=(π2)[11(ncl/nco)2].(13)

By plugging this expression into Eq. (12), it readily follows that Display Formula

PW=4π2a2ρpI(nclnco)[1(nclnco)].(14)

In addition, by combining this equation with Eq. (6), we obtain Display Formula

PW=PN(nclnco).(15)

This equation tells us that as (ncl/nco) diminishes, the relative contribution from the wide-angle rays also diminishes. This trend occurs because, as the ratio of the two indices drops, the narrow angle widens and the wide angle narrows. For the Saint Gobain fiber, where the core and cladding indices are 1.6 and 1.49, respectively, the contribution to the fluorescence power from the wide-angle rays is 93% of that from the narrow-angle rays. Even at a cladding index=1.0, the relative contribution is 62.5%.

The trapping efficiency (TE), as described at the beginning of this section, can be expressed as Display Formula

TE=(PW+PN)/4π2a2ρpl.(16)

It then follows, using Eqs. (6), (15), and (16), that the TE is Display Formula

TE=1(ncl/nco)2.(17)

We note that, within this treatment, the TE is independent of the volume density of fluorescers. This simple expression is plotted in Fig. 6. In this figure, the core index is the same as that of the Saint Gobain fiber, 1.6. The maximum cladding index is also the same as the Saint Gobain fiber, 1.49, and the minimum index is that of an air cladding, 1. The TE of the Saint Gobain fiber is 13%, while for an air cladding, it is 61%. The indices between these two are not included just to exhibit a trend; down to 1.3, they can be achieved with various degrees of fluorination of conventional polymers, for which ncl/nco=0.8125. Indices <1.3 to values <1.1 can be achieved by creating nanopores (i.e., small empty spaces) in a material such as silica.11

Graphic Jump Location
Fig. 6
F6 :

Trapping efficiency (TE) versus ncl/nco [TE=1(ncl/nco)2].

The TE of a fluorescent optical fiber has been calculated, assuming a uniform distribution of fluorescers over the cross-sectional area of the fiber, a uniform excitation of the fluorescers, and a uniform emission of fluorescence radiation over a solid angle. These assumptions are reasonable but could be relaxed. The calculation includes the contribution from rays over the full range of azimuth and those that are guided over only part of that range. The rays in the latter category make a substantial contribution to the total, but less than half, which diminishes along with the ratio of the cladding-to-core index. Furthermore, a simple analytic expression for the TE is presented.

Bross  A. D., “Scintillating plastic optical fiber radiation detectors in high-energy particle physics,” Proc. SPIE. 1592, , 122 –132 (1991). 0277-786X CrossRef
Thevenin  J. C.  et al., “Scintillating and fluorescent plastic optical fibers for sensors applications,” Proc. SPIE. 0514, , 133 –142 (1984). 0277-786X CrossRef
Amos  N. A., , Bross  A. D., and Lundin  M. C., “Optical attenuation length measurements of scintillating fibers,” Instrum. Methods Phys. Res. A. 297, (3 ), 396 –403 (1990). 0168-9002 CrossRef
Abbot  B.  et al., “Scintillating fibers and waveguides for tracking applications,” IEEE Trans. Nucl. Sci.. 38, (2 ), 441 –445 (1991). 0018-9499 CrossRef
Knoll  G. F., “Scintillation detection principles,” Chapter 8 in Radiation Detection and Measurement. , 4th ed., pp. 261 –270,  Wiley ,  Hoboken, NJ  (2010).
Weiss  J. D., “A fluorescent long-line fiber-optic position sensor,” Sens. Mag.. 22, (3 ), 28 –32 (2005).CrossRef
www.nanoptics.com (12  February  2015).
Snyder  A. W., and Love  J. D., “Leaky rays,” Chapter 7 in Optical Waveguide Theory. , 1st ed., pp. 134 –153,  Chapman and Hall ,  London, England  (1983).
Mercer  J. O., National Council of Teachers of Mathematics, private communication (2014).
“Low & ultra-low refractive index polymers,” (17  January  2011), http://omnexus.specialchem.com/tech-library/article/low-ultra-low-refractive-index-polymers (12  February  2015).

Jonathan D. Weiss is the chief technology officer of JSA Photonics LLC, a small high-technology company in Corrales, New Mexico. The company specializes in optical sensors. Before coming to JSA, he worked as a physicist at Sandia National Laboratories, where much of his work was in the area of optical sensors. He holds seventeen patents for such sensors, several of which have been licensed to private industry.

© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.

Citation

Jonathan D. Weiss
"Trapping efficiency of fluorescent optical fibers", Opt. Eng. 54(2), 027101 (Feb 17, 2015). ; http://dx.doi.org/10.1117/1.OE.54.2.027101


Figures

Graphic Jump Location
Fig. 1
F1 :

Emission of guided radiation from an on-axis fluorescer within a fluorescent optical fiber.

Graphic Jump Location
Fig. 2
F2 :

Skew and meridional ray paths within an optical fiber.

Graphic Jump Location
Fig. 3
F3 :

Geometry of a skew ray.

Graphic Jump Location
Fig. 5
F5 :

Angle of incidence, β, versus azimuthal angle, φ, for various values of θ and for (r/a)=0.5.

Graphic Jump Location
Fig. 4
F4 :

Angle of incidence, β, versus azimuthal angle, φ, for various values of (r/a) and for θ=32.1deg.

Graphic Jump Location
Fig. 6
F6 :

Trapping efficiency (TE) versus ncl/nco [TE=1(ncl/nco)2].

Tables

References

Bross  A. D., “Scintillating plastic optical fiber radiation detectors in high-energy particle physics,” Proc. SPIE. 1592, , 122 –132 (1991). 0277-786X CrossRef
Thevenin  J. C.  et al., “Scintillating and fluorescent plastic optical fibers for sensors applications,” Proc. SPIE. 0514, , 133 –142 (1984). 0277-786X CrossRef
Amos  N. A., , Bross  A. D., and Lundin  M. C., “Optical attenuation length measurements of scintillating fibers,” Instrum. Methods Phys. Res. A. 297, (3 ), 396 –403 (1990). 0168-9002 CrossRef
Abbot  B.  et al., “Scintillating fibers and waveguides for tracking applications,” IEEE Trans. Nucl. Sci.. 38, (2 ), 441 –445 (1991). 0018-9499 CrossRef
Knoll  G. F., “Scintillation detection principles,” Chapter 8 in Radiation Detection and Measurement. , 4th ed., pp. 261 –270,  Wiley ,  Hoboken, NJ  (2010).
Weiss  J. D., “A fluorescent long-line fiber-optic position sensor,” Sens. Mag.. 22, (3 ), 28 –32 (2005).CrossRef
www.nanoptics.com (12  February  2015).
Snyder  A. W., and Love  J. D., “Leaky rays,” Chapter 7 in Optical Waveguide Theory. , 1st ed., pp. 134 –153,  Chapman and Hall ,  London, England  (1983).
Mercer  J. O., National Council of Teachers of Mathematics, private communication (2014).
“Low & ultra-low refractive index polymers,” (17  January  2011), http://omnexus.specialchem.com/tech-library/article/low-ultra-low-refractive-index-polymers (12  February  2015).

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging & repositioning the boxes below.

Related Book Chapters

Topic Collections

PubMed Articles
Advertisement
  • Don't have an account?
  • Subscribe to the SPIE Digital Library
  • Create a FREE account to sign up for Digital Library content alerts and gain access to institutional subscriptions remotely.
Access This Article
Sign in or Create a personal account to Buy this article ($20 for members, $25 for non-members).
Access This Proceeding
Sign in or Create a personal account to Buy this article ($15 for members, $18 for non-members).
Access This Chapter

Access to SPIE eBooks is limited to subscribing institutions and is not available as part of a personal subscription. Print or electronic versions of individual SPIE books may be purchased via SPIE.org.