Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Easily fabricated, robust fiber-optic probe for weak fluorescence detection: modeling and initial experimental evaluation

Open Access Open Access

Abstract

Detection of fluorescence from a low quantum yield fluorophore is challenging for a fiber-optic probe, especially when an inexpensive and robust construction is desired. We propose a conceptually straightforward theoretical model to optimize the factors affecting the fluorescence-capture capability of a bifurcated/coaxial fiber-optic probe. Experimentally we verify that such a probe, if optimized, can detect the fluorescence of a polymer fluorophore with a low quantum yield of 0.0065.

©2012 Optical Society of America

1. Introduction

A fiber-optic probe wins over other means for chemical species analysis due to its well-known advantages: it is lightweight, flexible, non-toxic, low-loss, cost-effective and has a remote-monitoring capability. These features make it possible to construct an elegant hand-held probe for on-site chemical species assay in various environments [18].

A bifurcated fiber probe [38] or its improved version, a coaxial fiber probe [3, 4], is often adopted for chemical species assay, whose configuration is shown in Figs. 1(a) and 1(b), respectively. Either probe architectures could be easily constructed without sophisticated facilities. Its tiny probing end packs, for example, two or seven 430 μm fibers and the protecting sleeve into a cylinder approximately 2 mm in diameter. This small dimension makes bifurcated/coaxial fiber-optic probes widely used for chemical and biological species investigations and for in vitro/vivo detection [18]. However, the discussions in the relevant research papers seldom evaluate how well the flat-tipped probe detects the weak signal from a low quantum-yield fluorophore. Instead, a popular solution to cope with the weak fluorescent power is to use a fiber-optic probe with a better ability to capture light, such as one with a beveled probing tip. However, in comparison to the flat-tipped probe, the bevel-tipped end is hard to fabricate and offers poor protection of the sensing area. Given these issues, we investigate the use of the flat-tipped coaxial fiber-optic probe to detect the fluorescent signal of the low quantum-yield fluorophore under the simplest procedures and at the lowest possible cost. The fluorophore is a low band-gap polymer, represented here by P1, which contains a D-π-A system and has a low quantum yield (Q=0.0065).

 figure: Fig. 1

Fig. 1 Bifurcated and co-axial fiber-optic probes. (a) Typical bifurcated fiber-optic probe. (b) Coaxial fiber-optic probe for enhanced light collection. (c) The bifurcated fiber-probe with parameters for theoretical analysis. Only one overlapping volume is available to capture fluorescent power. S1 and S (shaded) are two cross sections of this volume. In contrast to (c), six identical active volumes (not drawn) and thus six active areas are available for a coaxial fiber-optic probe, leading to a six-fold increase of the overall captured fluorescent power.

Download Full Size | PDF

2. Improving the bifurcated/coaxial fiber-optic probe for low quantum-yield fluorescence detection

Detection of weak fluorescence requires an optimized coaxial fiber-probe with enhanced ability to capture light. To date, this subject has been investigated by a number of papers, with references 4 and 5 being perhaps the most rigorous and comprehensive but very complicated. In this paper, we propose a straightforward and intuitive approach to calculating the fluorescent-power capture ability of a flat-tipped bifurcated/coaxial fiber probe as shown in Fig. 1. Figure 1(c) gives detailed information on the two architectures required for the theoretical investigation. The excitation light and the captured signal light are limited respectively by the excitation cone (e-cone) and the receiving cone (r-cone) represented by2θ0. In the most common construction, which has two identical fibers, the angles of these two cones that open to the sample are equal and determined by their numerical aperture NA=nsinθ0, where n is the refractive index (RI) of the sample. The e-cone starts to overlap with the r-cone at the distance H0 from the probe tip. It is obvious that the captured fluorescent light from the sample can only arises from this overlapping or active volume. The shape of this active volume could be understood from the overlapping or active areas of S1 and S, which are formed by two cross sections of the two cones. The excitation power density Pe_d varies from H0 to H, or obeys:

Pe_d=Pe_d(h)const.

Thus the fluorescent powerIfshould be expressed by:

IfH0HPe_d(h)S(h)dh

Obviously, Eq. (2) has to include the factors representing the characteristics of the fluorophore. For a diluted solution, If is related to the excitation power intensityIe, molar absorptivity ε, quantum yield Q, the concentration c and the path length l by:

If=kQIeεcl
where k is a constant coefficient that may be determined by calibration. By applying Eq. (3) to Fig. 1(c), the total fluorescence power captured in the volume S(h)dh at any depth h(H0<hH) could be written as:
dIf=kQIe(h)εcdh=kQ[Pe_d(h)S(h)]εcdh
where

Ie(h)=Pe_d(h)S(h)

Clearly, Eq. (4) is equivalent with Eq. (2) but includes detailed information describing the fluorescence-capable solution.

Equation (5) indicates the overall excitation power incident on the overlapping area S(h)at depth h(H0<hH) while Eq. (4) assumes that within the small thickness dh the value Ie(h) remains unchanged.

To simplify the discussion (a reasonable approach for a diluted solution), we ignore the attenuation of the excitation power when traveling in the solution, and thus:

Pe_d(h)S(h)=Pe_d(0)S(0)

In contrast to the rigorous but complex treatment of the excitation light beam in [4], Eq. (6) assumes a uniform distribution of excitation power across the entire cross section of the light beam, which not only greatly simplifies the derivation, but also makes the concept easier to understand. Any error associated with this approximation is likely due to overestimating the excitation power in the overlapping volume, and could be greatly reduced by including the power penalty factor in the coefficient kin Eq. (4). Critically, since all calculation results given later (Fig. 2 ) are based on ratio to obtain normalized intensities and thus allow the coefficient k to be cancelled out. From this point of view, the results from such an approximation provide a satisfactory degree of accuracy.

 figure: Fig. 2

Fig. 2 Some important conclusions from Eq. (15) are illustrated here. The data points or curve labeled as “This work” are associated with the coaxial fiber-optic probe in our experiment. The shaded areas of each figure indicate the parameters fixed for all curves in that figure. (a) Normalized captured power If at h = 0 vs. cladding thickness when H = 1000 μm. (b) Same as (a) except H = 3000 μm. (c) Impact of NA on If for two fibers, one with r = 100 μm, the other with r = 200 μm. (d) Impact of NA on initial depth H0 defined in Fig. 1(c), which describes the dead space in front of the fibers. (e) If vs the maximum detection depth H for the fiber used in this work and for the standard multimode fiber in fiber-optic communication. (f) Significant impact of fiber core radius r alone on If for fixed t, NA and H.

Download Full Size | PDF

With S(0)=πr2 and S(h)=πR2, where r andR are the radius of the fiber core and the radius of the cross section of the e-cone or r-cone at the depth h(H0<hH), we can write:

Pe_d(h)=Pe_d(0)(rR)2

Any overlapping area S(h) has the same the geometry as S=S(H)in Fig. 1(c). The following expression is not difficult to find:

S(h)=2[R2sec1(Rr)rR2r2]

Replacing the associated terms in Eq. (4) with Eqs. (7) and Eq. (8) yields:

dIf(h)=2×kQPe_d(0)(rR)2[R2sec1(Rr)rR2r2]εcdh

By assuming γ=2×kQεcPe_d(0)and lossless transmission of dIf from depth h=H to h=0, the captured fluorescent power at the entrance of the r-fiber (h=0), represented by If(h=0), could be straightforwardly written from Eq. (9):

If(h=0)=γH0H[(rR)2(R2sec1(Rr)rR2r2)]dh

We change Eq. (10) to the form:

If(h=0)=γr2H0H[sec1(Rr)rR1(rR)2]dh

Figure 1(c) validates the following relations:

R=htgθ0+r,h=(Rr)ctgθ0,

Forh=H0=tctgθ0, Eq. (12) givesR0=r+t, where tis the cladding thickness.

Equation (11) thus becomes:

If(h=0)=γr2R0R[sec1(Rr)rR1(rR)2]d[(Rr)ctgθ0]=γR0R(sec1(Rr)rR1(rR)2)dR
where γ=γr2ctgθ0

Equation (13) is further simplified by assuming u=Rr and w=rR:

If(h=0)=γ{rR0rRrsec1udurrR0rRw1w2d1w}

Assuming If=If(h=0) Eq. (14) yields:

If=γ{[RrarcsecRrln(Rr+(Rr)21][R0rarcsecR0rln(R0r+(R0r)21]}+γ{[1(rR)2ln1+1(rR)2rR][1(rR0)2ln1+1(rR0)2rR0]}
where

γ=2×kQεcPe_d(0)r3ctgθ0

Equation (15) expresses the overall fluorescent power captured by the r-fiber from the overlapping volume between H0andHfor a lossless liquid solution. For a coaxial fiber-optic probe, there are six identical overlapping volumes that contribute to the overall captured fluorescent power. In this case, we can write:

If_coaxial=6If

The calculation of Eqs. (15) and Eq. (16) generates the results presented in Fig. 2, from which we can draw the following conclusions. First, the core radius rgreatly affects If via r3, which is reflected by all curves in Figs. 2(a), 2(b), 2(c), 2(e) and 2(f). Second, from Figs. 2(a) and 2(b), when a greater detecting depth H is considered, which is meaningful for an extremely low or a lossless sample medium, cladding thickness variation affects Ifon a smaller scale. This becomes obvious when a larger core fiber is involved. A core radius as small as 5 μm, as in a single-mode fiber, yields If0 for any cladding thickness. Third, a larger NA could lead to a better capturing capability via ctgθ0 as described by Fig. 2(c), which also indicates a rapid rise of If for a larger core fiber. Fourth, the benefit of a large NA also includes a shorter initial depth H0, as shown in Fig. 2(d). The depth H0 is proportional to the dead space in front of the probe tip. While this might not be a concern for lossless medium such as we assume in this paper for Eqs. (10) and (15), when a lossy medium fills this space – as in fact occurs to most samples – Eqs. (10) and (15) must include the attenuation coefficients for both excitation and fluorescent power, and If<<Ifmay result. Fifth, instead of altering the NA, using a smaller cladding thickness t is another efficient way to reduce the dead space, as indicated in Fig. 2(d). This is perhaps a more meaningful option since all spectrometers accepting fiber-optic connection are designed for NA=0.22, which is the lowest NA found in Figs. 2(c) and 2(d). For this reason, Fig. 2 emphasizes the results for NA=0.22. Sixth, the standard multimode fiber for fiber-optic communication, with its core/cladding diameter of 50 / 125 μm and NA=0.22, is an inexpensive and convenient fiber to use in construction of the bifurcated/coaxial fiber-optic probe. However, the limited capability of this fiber to capture fluorescent power, as evidenced in Fig. 2(e), contrasts with the fluorescence-capture ability of the fiber with core/cladding diameter of 400/430 μm and NA=0.22 used here. Seventh, the high performance of the fiber-optic probe presented in this work is evident in all parts of Fig. 2. For the reasons outlined above, this must be considered the optimal design taking into account all the factors discussed so far.

3. Experiments and discussions

A P1 solution is prepared by dissolving 2.9 mg P1 in 1 ml of CHCL3 and stirring for one hour and further diluted in 10 ml DMF solvent. After filtering, the P1 solution (3 × 10−4 M) is ready. For the purpose of performance comparison, we also prepared a rhodamine 6G (R6G) water solution (29 μg / ml) as the high quantum-yield (Q=0.95) fluorophore.

The coaxial fiber-optic probe shown in Fig. 1(b), with optimized performance as illustrated in Fig. 2, is used for fluorescence measurement. All fibers are identical, having a core/cladding diameter of 400/430 μm and NA of 0.22. The excitation light (633 nm for P1 and 532 nm for R6G) is launched to the e-lead that comprises six e-fibers surrounding one r-fiber in the probing end. The r-fiber constitutes the r-lead and is connected to an Ocean Optics QE65000 spectrometer, which is interrogated by a laptop via USB connection. The Ocean Optics SpectraSuite software package is used to record the fluorescence spectra.

In the first experiment, the fluorescent spectrum is recorded by immersing the fiber probing tip into the R6G solution. Similarly, the second experiment is carried out by immersing the probing end into the P1 solution. The experimental result in Fig. 3(a) shows that the R6G affords a high quality fluorescent spectrum owing to its high quantum yield (Q = 0.95) that requires only 3 ms integration time, with complete elimination of the excitation laser trace even without an emission filter. Yet, despite the very low quantum yield (Q = 0.0065) of the P1 solution, which produces a fluorescent power 150 times weaker than R6G, Fig. 3(b) indicates that the optimized coaxial probe is still able to capture the fluorescent spectrum from P1, albeit with an integration time of 5,000 ms. However, this 5,000 ms integration time leads to a strong excitation-power trace in the captured spectrum in Fig. 3(b), an effect stemming from the limitation of the fiber probe itself. Since a fiber-optic-based system affords many well known advantages but at the cost of an enormous power penalty during both excitation and fluorescent light coupling, a high integration time is needed in order for the QE65000 spectrometer to respond to the weak P1 fluorescence.

 figure: Fig. 3

Fig. 3 Fluorescence spectra of R6G and P1 in DMF captured by the optimized coaxial fiber-optic probe. (a) Strong and high quality fluorescence of R6G without an excitation laser trace; (b) Weak fluorescence of P1 with a strong excitation laser trace. Note the significant differences in their integration times and the quantum yields.

Download Full Size | PDF

4. Conclusions

In conclusion, this paper proposes a straightforward model for calculating and optimizing the ability of a coaxial fiber-optic probe to capture fluorescent light from a lossless liquid. The analytical expression obtained shows how the captured fluorescent power is affected by such important parameters as the molar absorptivity, the sample concentration, the quantum yield, the excitation power intensity, the dimensions of the fiber core and cladding, and the value of NA. The fiber parameters for maximum possible fluorescent power capture, presented here in a series of graphs, could be used to optimize the fiber-optic probe by matching the demands of various measurement environments and a particular spectrometer design. Evaluation of the performance of the coaxial fiber-optic probe optimized on the basis of this model shows a satisfactory quality of the fluorescence spectrum captured from a polymer fluorophore with an extremely low quantum yield of 0.0065. The lowest concentration this probe can cope with is 3 × 10−4 M.

The work presented here paves the way towards using the coaxial fiber probe to detect trace heavy metal by turning on and off fluorescence of this P1 polymer via titration measurement, which will be discussed in our next paper.

Acknowledgments

The authors gratefully acknowledge support for this work from the Natural Sciences and Engineering Research Council of Canada, from the Canada Research Chairs Program and from the Ministère de Développement économique, Innovation et Exportation du Québec.

References and links

1. B. A. Flusberg, E. D. Cocker, W. Piyawattanametha, J. C. Jung, E. L. M. Cheung, and M. J. Schnitzer, “Fiber-optic fluorescence imaging,” Nat. Methods 2(12), 941–950 (2005). [CrossRef]   [PubMed]  

2. M. G. Shim, B. C. Wilson, E. Marple, and M. Wach, “Study of fiber-optic probes for in vivo medical Raman spectroscopy,” Appl. Spectrosc. 53(6), 619–627 (1999). [CrossRef]  

3. L. Wen-xu and C. Jian, “Continuous monitoring of adriamycin in vivo using fiber optic-based fluorescence chemical sensor,” Anal. Chem. 75(6), 1458–1462 (2003). [CrossRef]   [PubMed]  

4. T. F. Coony, H. T. Skinner, and S. M. Angel, “Comparative study of some fiber-optic remote Raman probe designs. Part I: model for liquids and transparent solids,” Appl. Spectrosc. 50(7), 836–848 (1996). [CrossRef]  

5. U. Bunting, F. Lewitzka, and P. Karlitschek, “Mathematical model of a laser-induced fluorescence fiber-optic sensor head for trace detection of pollutants in soil,” Appl. Spectrosc. 53(1), 49–56 (1999). [CrossRef]  

6. K. R. Rogers and E. J. Poziomek, “Fiber optic sensors for environmental monitoring,” Chemosphere 33(6), 1151–1174 (1996). [CrossRef]  

7. J. Ma and W. Bock, “Modeling of photonic crystal fiber with air holes sealed at the fiber end and its application to fluorescent light collection efficiency enhancement,” Opt. Express 13(7), 2385–2393 (2005). [CrossRef]   [PubMed]  

8. X. B. Zhang, C. C. Guo, Z. Z. Li, G. L. Shen, and R. Q. Yu, “An optical fiber chemical sensor for mercury ions based on a porphyrin dimer,” Anal. Chem. 74(4), 821–825 (2002). [CrossRef]   [PubMed]  

Cited By

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

Alert me when this article is cited.


Figures (3)

Fig. 1
Fig. 1 Bifurcated and co-axial fiber-optic probes. (a) Typical bifurcated fiber-optic probe. (b) Coaxial fiber-optic probe for enhanced light collection. (c) The bifurcated fiber-probe with parameters for theoretical analysis. Only one overlapping volume is available to capture fluorescent power. S1 and S (shaded) are two cross sections of this volume. In contrast to (c), six identical active volumes (not drawn) and thus six active areas are available for a coaxial fiber-optic probe, leading to a six-fold increase of the overall captured fluorescent power.
Fig. 2
Fig. 2 Some important conclusions from Eq. (15) are illustrated here. The data points or curve labeled as “This work” are associated with the coaxial fiber-optic probe in our experiment. The shaded areas of each figure indicate the parameters fixed for all curves in that figure. (a) Normalized captured power If at h = 0 vs. cladding thickness when H = 1000 μm. (b) Same as (a) except H = 3000 μm. (c) Impact of NA on If for two fibers, one with r = 100 μm, the other with r = 200 μm. (d) Impact of NA on initial depth H0 defined in Fig. 1(c), which describes the dead space in front of the fibers. (e) If vs the maximum detection depth H for the fiber used in this work and for the standard multimode fiber in fiber-optic communication. (f) Significant impact of fiber core radius r alone on If for fixed t, NA and H.
Fig. 3
Fig. 3 Fluorescence spectra of R6G and P1 in DMF captured by the optimized coaxial fiber-optic probe. (a) Strong and high quality fluorescence of R6G without an excitation laser trace; (b) Weak fluorescence of P1 with a strong excitation laser trace. Note the significant differences in their integration times and the quantum yields.

Equations (17)

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

P e_d = P e_d (h)const.
I f H 0 H P e_d (h)S(h)dh
I f =kQ I e εcl
d I f =kQ I e (h)εcdh=kQ[ P e_d (h)S(h) ]εcdh
I e (h)= P e_d (h)S(h)
P e_d (h)S(h)= P e_d (0)S(0)
P e_d (h)= P e_d (0) ( r R ) 2
S(h)=2[ R 2 sec 1 ( R r )r R 2 r 2 ]
d I f (h)=2×kQ P e_d (0) ( r R ) 2 [ R 2 sec 1 ( R r )r R 2 r 2 ]εcdh
I f (h=0)=γ H 0 H [ ( r R ) 2 ( R 2 sec 1 ( R r )r R 2 r 2 ) ]dh
I f (h=0)=γ r 2 H 0 H [ sec 1 ( R r ) r R 1 ( r R ) 2 ]dh
R=htg θ 0 +r, h=(Rr)ctg θ 0 ,
I f (h=0)=γ r 2 R 0 R [ sec 1 ( R r ) r R 1 ( r R ) 2 ]d[ (Rr)ctg θ 0 ] = γ R 0 R ( sec 1 ( R r ) r R 1 ( r R ) 2 )dR
I f (h=0)= γ { r R 0 r R r sec 1 udur r R 0 r R w 1 w 2 d 1 w }
I f = γ { [ R r arcsec R r ln( R r + ( R r ) 2 1 ][ R 0 r arcsec R 0 r ln( R 0 r + ( R 0 r ) 2 1 ] } + γ { [ 1 ( r R ) 2 ln 1+ 1 ( r R ) 2 r R ][ 1 ( r R 0 ) 2 ln 1+ 1 ( r R 0 ) 2 r R 0 ] }
γ =2×kQεc P e_d (0) r 3 ctg θ 0
I f_coaxial =6 I f
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.