Theoretical Design of A High Sensitivity SPR-based Optical Fiber Pressure Sensor
Theoretical Design of A High Sensitivity SPR-based Optical Fiber Pressure Sensor
fully edited. Content may change prior to final publication. Citation information: DOI
10.1109/JLT.2015.2477353, Journal of Lightwave Technology
JLT-17721-2015 1
Copyright (c) 2015 IEEE. Personal use of this material is permitted. However, permission to use this material for any other
purposes must be obtained from the IEEE by sending a request to pubs-permissions@ieee.org.
0733-8724 (c) 2015 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission. See
http://www.ieee.org/publications_standards/publications/rights/index.html for more information.
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI
10.1109/JLT.2015.2477353, Journal of Lightwave Technology
JLT-17721-2015 2
where Rcl, ncl, Rcr, and ncr are the radius and refractive index of
the cladding and core fiber, respectively.
C. Numerical Simulation
As already described, the formation of the SPR will only
happen since the necessary conditions of energy transfer given
by the fiber cladding, the metal layer and the silicon rubber
medium of the block are satisfied. For these three layers, a
transfer matrix formalism was used in the simulation of the
system and the reflectance was calculated based on an N-layer
model [14], [16], [17]. Maxwell equations will give the
solution that relates the electromagnetic amplitude quantities,
electric (E) and magnetic (H) field vectors, by the formula:
= [M ]× N
E0 E (6)
H0 HN
− i sin δ k
M 12 N −1 cos δ k
[M ] =
M 11
=∏
(7)
ηk
M 22 k =1
21 cos δ k
M
Fig. 2. Schematic representation of the U-shaped SPR optical fiber probe − i η k sin δ k
inside the rubber block. All parameters here presented are described in the
text. with δk the phase vector and ηk the optical admittance of the kth
the normal of the incident plane is lower than the critical layer which, for p-polarized light, are dependent of the
angle. Otherwise, the light ray will continue its path inside the dielectric function εk, the thickness dk of the layer, the incident
core. All rays capable of propagate inside the cladding angle θ0 and wavelength λ, defined by:
medium will then encounter the outer radius metal layer with
an incident angle ϕ with respect to the normal of incident 2πd k
δk = (ε k − n0 sin 2 θ 0 )1 / 2 (8)
plane. This incident light will have all the requirements to λ
create plasmonic wave excitation for certain wavelengths
while the other rays will reflect and continue their paths εk
ηk = . (9)
through the cladding medium, directing up to the fiber core. (ε k − n0 sin 2 θ 0 )1 / 2
Once it reaches the cladding-core boundary, Snell refraction
phenomena will happen, propagating through the core until it The p-polarized reflectance of the structure is finally given
reaches the inner core-cladding boundary and, when possible, in term of Fresnel reflection coefficient (frp) that is related to
the inner radius metal layer. Therefore, the incident rays will the elements of M by:
have an angle δ with respect to the normal of incident plane.
Again, plasmonic wave excitation for certain wavelengths will 2
2 ( M 11 + M 12 ⋅ η N )η 0 − ( M 21 + M 22 ⋅ η N ) . (10)
occur. Non-absorbed light will then be reflected and Frp = f rp =
propagated through a new cycle until the bending region is ( M 11 + M 12 ⋅ η N )η 0 + ( M 21 + M 22 ⋅ η N )
over.
The reflection angles in the outer (φ) and the inner (δ)
surfaces of the U-region can be calculated with the application In this proposed structure it is required two different
of the sine rule and Snell law considering the dependence of reflection angles (φ and δ in (4) and (5)) corresponding to the
both parameters on h and θ. inner and outer metal interface. For each angle and, by
association with each interface, it will also have two different
calculated reflectance values in (10), namely Fφ and Fδ where
n R + Rcl − Rcr + h θ0, in (8), will be respectively φ and δ.
ϕ = sin −1 cr ⋅ sin(θ ) (4)
ncl R + 2 Rcl For the dependency of the refractive index of the silica fiber
core, in this case a step-index multimode fiber structure, with
Copyright (c) 2015 IEEE. Personal use of this material is permitted. However, permission to use this material for any other
purposes must be obtained from the IEEE by sending a request to pubs-permissions@ieee.org.
0733-8724 (c) 2015 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission. See
http://www.ieee.org/publications_standards/publications/rights/index.html for more information.
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI
10.1109/JLT.2015.2477353, Journal of Lightwave Technology
JLT-17721-2015 4
the wavelength, it was used the Sellmeier dispersion relation, completes a full reflection cycle, as seen in Fig. 2. These
given by the following expression [18]: quantities will have dependence on the morphology of the
fiber and its bent radius:
a1 λ 2 a 2 λ2 a3 λ2 (11)
ncr = 1 + + + R + Rcl − Rcr + h
λ2 − b12 λ2 − b22 λ 2 − b32 d1 = θ − sin −1 sin(θ ) × (R + Rcl ) (15)
R + Rcl + Rcr
where a1, a2, a3, b1, b2 and b3 are the Sellmeier coefficients
with the respective values of a1=0.6961663, a2=0.4079426, n R + Rcl − Rcr + h
d 2 = 2 × sin −1 cr ⋅ sin(θ )
a3=0.8974794, b1=0.068404×(T/300)2, b2=0,1162414×(T/300)2 ncl R + Rcl + Rcr (16)
and b3= 9.896161, with the temperature T in Kelvin and the
n R + Rcl − Rcr + h
wavelength λ in micrometer. The refractive index of the fiber − sin −1 cr ⋅ sin(θ ) × (R + Rcl )
cladding is directly derived from its fixed numerical aperture, ncl R + 2 Rcl
NA.
The dependence of the refractive index of the metal layers R + Rcl − Rcr + h
with the wavelength was determined using the experimental d 3 = sin −1 sin(θ )
data obtained by Johnson and Christy [19]. The values were R + Rcl − Rcr (17)
converted to the dielectric function and an interpolation was R + Rcl − Rcr + h
applied between the existing data. − sin −1 sin(θ ) × (R + Rcl )
Assuming a collimated light source launched at one end of R + R cl + R cr
the optical fiber and focused with a microscopic objective to
the fiber end-face, where the output is collected at the other n R + Rcl − Rcr + h
end of the fiber with a detection system, the power distribution d 4 = 2 × sin −1 cr ⋅ sin(θ )
of light can be expressed as [20]: ncl R (18)
n R + Rcl − Rcr + h
n 2 sin θ cosθ − sin −1 cr ⋅ sin(θ ) × (R + Rcl )
P(θ )dθ ∝ 0 2 dθ (12) ncl R + Rcl − Rcr
(1 − n0 cos 2 θ ) 2
With the detection of the transmitted power resonant dip
where θ is the angle of the ray with the normal of the core– shifts by a spectral analyzer, it is possible to determine the
cladding interface. The range of θ values possible for the fiber bending variation (Rv) and, therefore, the applied
launched rays at the input of the fiber is limited in its lower pressure (Pr) to the block. Knowing the Young modulus (Y) of
value depending on the fiber core and cladding refractive the silicon rubber and the initial bending radius (R0) of the
indices. This dependence is called critical angle, given by the sensor, we have the following relation:
relation: θcr=sin-1(ncl/ncr). The maximum value of θ is π/2. The
h value is related to the distance between the core-cladding R − R0 .
interface and the spatial point where the ray reaches the Pr = εY with ε = v (19)
boundary of the straight and the U-shaped region. The values R0
of h can range from 0 to the fiber core diameter (2Rcr). The
contribution of all rays and multi-reflections are necessary to The sensitivity (S) of the sensor will be determined by the
determine the total transmitted power: correlation of the resonance wavelength shift δλres to
increments of pressure δPr defined by:
π
2 Rcr n02 sin θ cos θ
∫ ∫θ
0
2
[ Fϕ (θ , h) × Fδ (θ , h)] N R (θ ,h )
(1 − n02 cos 2 θ ) 2
dθ dh
δλres
PT = S=
cr
π
(20)
2 Rcr n02 sin θ cos θ δPr
∫ ∫θ
0
2
cr (1 − n02 cos 2 θ ) 2
dθ dh
Copyright (c) 2015 IEEE. Personal use of this material is permitted. However, permission to use this material for any other
purposes must be obtained from the IEEE by sending a request to pubs-permissions@ieee.org.
0733-8724 (c) 2015 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission. See
http://www.ieee.org/publications_standards/publications/rights/index.html for more information.
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI
10.1109/JLT.2015.2477353, Journal of Lightwave Technology
JLT-17721-2015 5
Copyright (c) 2015 IEEE. Personal use of this material is permitted. However, permission to use this material for any other
purposes must be obtained from the IEEE by sending a request to pubs-permissions@ieee.org.
0733-8724 (c) 2015 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission. See
http://www.ieee.org/publications_standards/publications/rights/index.html for more information.
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI
10.1109/JLT.2015.2477353, Journal of Lightwave Technology
JLT-17721-2015 6
REFERENCES
[1] H. Fu, M. T. C. Wu, L. Zhang, K. Cheng, H. Tam, B. Guan and C. Lu,
“High pressure sensor based on photonic crystal fiber for downhole
application,” Appl. Opt., vol. 49, no. 14, pp. 2639-2643, May 2010.
[2] D. Yun, L. Haijun, D. A. Olson and Y. Miao, “Hybrid Miniature Fabry–
Perot Sensor with Dual Optical Cavities for Simultaneous Pressure and
Temperature Measurements,” J. Lightw. Technol., vol. 32, no. 8, pp.
1585-1593, 2014.
[3] A. Sapienza, R. Rocha, J.P. Faria, F.S. Ferrari and A. Paterno, “Surface
plasma analysis on a palladium cylindrical shell covering weakly guided
silica optical fiber” Microwave and Optoelectronics Conference, 2007.
IMOC 2007. SBMO/IEEE MTT-S International, pp.260-264, 2007
[4] S. K. Srivastava and B. D. Gupta, “Fiber Optic Plasmonic Sensors: Past,
Present and Future,” Open Opt. J., vol. 7, pp. 58-83, Dec. 2013.
[5] M. Consales, A. Ricciardi, A. Crescitelli, E. Esposito, A. Cutolo and A.
Cusano, “Lab-on-Fiber Technology: Toward Multifunctional Optical
Nanoprobes,” ACS Nano, vol. 6, no. 4, pp. 3163-3170, 2012.
[6] H. P. Ho and C. L. Wong, “Application of spectral surface plasmon
resonance for gas pressure sensing,” in Proc. IEEE Electron Devices
Meeting (IEEE, 2002), pp. 73-77, 2002.
[7] C. H. Wong, H. Ho, K. Chan, S. Wu and C. Lin, “Application of spectral
surface plasmon resonance to gas pressure sensing,” Opt. Eng., vol. 44,
no. 12, 124403, 2005.
[8] N. Wu, X. Zou, Y. Tian, J. Fitek, M. Maffeo, C. Niezrecki, J. Chen and
X. Wang, “An ultra-fast fiber optic pressure sensor for blast event
measurements,” Meas. Sci. Technol., vol. 23, no. 5, 055102, 2012.
[9] W. Zhang, L. Liu, F. Li and Y. Liu, “Fiber Bragg grating pressure sensor
with enhanced sensitivity,” Chin Opt Lett, vol. 5, no. 9, pp. 507-508,
2007.
[10] R. Nieuwland, L. Cheng, P. Harmsma, M. Lemmen, R. Oostenbrink and
J. Schreuder, “A fiber Bragg based semi distributed pressure sensor
system for in-vivo vascular applications,” in Advanced Photonics, OSA
Technical Digest, SeW1C.1, 2014.
[11] K. Bhowmik, G. Rajan, E. Ambikairajah and G. Peng, “Hydrostatic
Fig. 4. a) SPR wavelength shift with applied pressure for the three metal thin
pressure sensitivity of standard polymer fibre Bragg gratings and etched
films. b) Sensitivity variation with pressure for the three metals.
polymer fibre Bragg gratings,” in Proc. SPIE 9157, 91573G, 2014.
[12] H. Mi, Y. Wang, P. Jin and L. Lei, “Design of a ultrahigh-sensitivity
range that the sensor here proposed can achieve (0.12 MPa). SPR-based optical fiber pressure sensor,” Optik, vol. 124, no. 21, pp.
The reported optical sensors have an operation range up to 5248–5250, 2013.
0.58 MPa for FP interferometer and up to 1 MPa for the Bragg [13] R. K. Verma and B. D. Gupta, “Theoretical modelling of a bi-
dimensional U-shaped surface plasmon resonance based fibre optic
grating, significantly higher than the proposed sensor, being sensor for sensitivity enhancement,” J. Phys. D: Appl. Phys., vol. 41,
this parameter one of its main disadvantages in comparison. no.9, 095106, 2008.
[14] D. Duarte, N. Alberto, L. Bilro and R. Nogueira, “Theoretical modeling
IV. CONCLUSION of an U-shaped SPR fiber sensor in 1550-nm spectral range for sensing
applications,” in Proc. SPIE 9286, 928640, 2014.
A novel theoretical design of a high sensitivity SPR-based [15] S. A. Maier, Plasmonics: Fundamentals and Applications, Springer,
optical fiber pressure sensor was presented. This design uses a 2007.
[16] M. Kanso, S. Cuenot and G. Louarn, “Sensitivity of optical fiber sensor
new U-shaped configuration embedded in a silicon rubber based on surface plasmon resonance: modeling and experiments,”
block and the advantages of this solution compared with PLASMONICS, vol. 3, no. 2-3, pp. 49-57, 2008.
others already proposed in the literature is the maintaining of [17] M. Born and E. Wolf, Principles of Optics: Electromagnetic Theory of
the fiber cladding and thus the easier of the overall production Propagation, Interference and Diffraction of Light, Cambridge
University Press, 1999.
process. This sensor is also designed to work in the optical [18] A. K. Ghatak and K. Thyagarajan, Introduction to Fiber Optics,
communication wavelengths, so the already available Cambridge University Press, 1998.
components and devices operating at these bands can be used. [19] P. B. Johnson and R. W. Christy, “Optical Constants of the Noble
The simulation results showed that for gold thin film layer, Metals,” Phys. Rev. B, vol. 6 no. 12, pp. 4370-4379, 1972.
[20] Y. Xu, N. B. Jones, J. C. Fothergill and C. D. Hanning, “Analytical
theoretical sensitivities higher than 1700 nm/MPa can be estimates of the characteristics of surface plasmon resonance fibre-optic
achieved, for pressure values below 0.03 MPa, corresponding sensors,” J. Mod. Opt., vol. 47, no. 6, pp. 1099-1110, 2000.
to a resolution of 5.9x10-4 kPa. These sensitivities decrease [21] B. D. Gupta and R. K. Verma, “Surface Plasmon Resonance-Based
with increasing pressure until 0.12 MPa. A similar behavior is Fiber Optic Sensors: Principle, Probe Designs, and Some Applications,”
J. Sens., vol. 2009, 979761, 2009.
Copyright (c) 2015 IEEE. Personal use of this material is permitted. However, permission to use this material for any other
purposes must be obtained from the IEEE by sending a request to pubs-permissions@ieee.org.
0733-8724 (c) 2015 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission. See
http://www.ieee.org/publications_standards/publications/rights/index.html for more information.
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI
10.1109/JLT.2015.2477353, Journal of Lightwave Technology
JLT-17721-2015 7
Copyright (c) 2015 IEEE. Personal use of this material is permitted. However, permission to use this material for any other
purposes must be obtained from the IEEE by sending a request to pubs-permissions@ieee.org.
0733-8724 (c) 2015 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission. See
http://www.ieee.org/publications_standards/publications/rights/index.html for more information.