劉慶敏，侯尚林，雷景麗

（蘭州理工大學(xué)理學(xué)院，蘭州 730050）

0 Introduction

Surface Plasmon Resonance（SPR）is a prominent optical phenomenon that arises as to the extent of energy transferring from photons to surface plasmon waves under appropriate conditions［1］.In the past few years，this optical effect，owing to its high sensitivity，real-time detection，and anti-interference，has already been extensively investigated and applied in the fields of medical treatment［2］，environment monitoring［3］，biomedical sensing［4］and so on.The earliest SPR sensing devices were based on prism［5］，such as the Kretschmann-Raether prism［6］，whose sensitivity reached up to 10-6RIU.Even though the prism can achieve high sensitivity，it has the disadvantage of costly integration，limited mechanical reliability and difficulties in mass production，which were limited to use in remote monitoring.Optical fiber-based SPR sensors have been exploited in recent years due to their advantages of easy miniaturization，low manufacturing cost and in-situ monitoring［7］，etc.Especially，the sensors based on Photonic Crystal Fiber（PCF）can be easy phase-matching and encapsulation because of the minimized size，tunable geometric parameters and flexible design［8］，so it becomes a promising candidate for assembling compact and ultra-sensitive SPR sensors with fast，sustainable and real-time monitoring.

As is well known that the sensing performance of Photonic Crystal Fiber-surface Plasmon Resonance（PCF-SPR）sensors depends on the metallic materials and the structure of PCF.Among the various plasmonic materials，gold［9］and silver［10］are mostly used due to their relatively low loss in the visible and near-infrared region，which makes it more effective to monitor the high energy transfer on the surface of plasmonic materials.To date，much effort has been devoted to the structure design of PCF to enhance the sensing performance of PCF-SPR sensors［11］，and considerable variety of PCF-SPR structures have been reported，which can be divided into two categories，internal metal-coated structures using selectively coating or nanowires filled，and external metal-coated structures such as D-shaped，slotted and exposed core PCFs［12］.Among them，D-shaped PCF has attracted many interests compared with other PCF structures，because it can not only provide a good solution to the difficulties in a uniform coating which means the precise control of the thickness and roughness of the metal film on the outer surface of the optical fiber，but also solve the problem of filling analytes to Nano-size holes［13］.For example，AN Guowen，et al［14］proposed a D-shaped PCF-SPR sensor with a triangular lattice and four large-size channels，which can obtain the maximum wavelength sensitivity of 10 493 nm/RIU at 1.38.SINGH S and PRAJAPATI Y［15］presented a D-shape PCF-SPR sensor with gold-graphene layers on the surface.The sensitivity is 33 500 nm/RIU and the effective RI resolution is as high as 2.98×10-5RIU for analyte RIs between 1.32 and 1.40.It can be seen from these reported literatures that the most reported Dshaped SPR sensors exhibit good sensing performance when the Refractive Index（RI）of the analyte is higher than 1.30.

However，with the development of biomedicine and material chemistry，the detection for low refractive index analyte has been paid more attention.For example，the RI of sevoflurane in the drug is around 1.27［16］，the refractive index of some fluorine-containing organics compounds in water pollution is about 1.30［17］.Hence，it is essential for us to exploit sensors with high sensitivity for detecting low RI.In recent years，some research work about the RI of analytes lower than 1.30 has been conducted.For example，WANG Famei et al［18］studied a D-shaped sensor with two parallel to detect low RI ranging from 1.27 to 1.32 with the maximum spectral sensitivity up to 13 500 nm/RIU and the resonance wavelength ranging between 720 nm and 1 680 nm.LIU Chao et al［19］reported a mid-infrared PCF-SPR sensor with RI detection range of 1.23～1.29 and average sensitivity of 5 500 nm/RIU.CHEN Xin et al［20］proposed a D-type PCF-SPR sensor which operates in the mid-infrared band from ranging 2 530 nm to 2 750 nm with RI of 1.26～1.29 and the maximum sensitivity of 11 055 nm/RIU.In most of these studies，the reports of PCF-SPR sensor with operating wavelength in midinfrared range are relatively few.In fact，SPR sensors operating in the mid-infrared wavelength have unique advantages compared to the visible and near-infrared bands［21］.For example，it is conducive to improve the detection accuracy and sensitivity of the sensor to large samples as living cells［22］，and avoid the light damage to the biological samples［23］.Therefore，it is necessary to exploit the PCF-SPR sensor with higher sensitivity in mid-infrared waveband.

Herein，a novel D-shaped PCF-SPR sensor for detecting low RI analytes is proposed and analyzed.The open-ring channel coated with gold film can not only reduce the coating area，but also simplify the fabrication process.Numerical results show that a maximum spectral sensitivity of 38 800 nm/RIU can be obtained when analyte RI is 1.30 with a high resolution of 2.37×10-6RIU.The operating wavelength of this fiber sensor includes the near-and mid-infrared regions ranging from 2 020 nm to 3 036 nm.This proposed sensor is of significance in environmental engineering，biosensors and healthcare.

1 Design and analysis of the model

The cross-section of the proposed D-shaped SPR sensor is illustrated in Fig.1.This structure contains three layers of air rings.The radii of air holes in the first-layer and third-layer arer1andr3，respectively.While the second-layer air ring consists of air holes with two different radii，r2andrs.The refractive index of air is fixed atnair=1 and the radius of the cladding isR.A thin gold film with thicknesstgis deposited on the inner surface of the micro-opening analyte channel on the upper side，the radius and the center location of the channel arersand 2.5×Λ-1.25×rs，respectively.The first layer of air holes is composed of regular octagonal air holes with a radius ofr1and a distance ofΛfrom the fiber core，the second layer of air holes consists of a large air hole with radius ofrsand ten small air holes with radius ofr2.The large air hole is located directly below the fiber core and the distance from the fiber core is 2.5×Λ-1.25×rs，while the small air holes with radius ofr2are arranged in two ways.One is distributed on both sides of the large air hole and the distance from the fiber core is 2.5×Λ-1.25×rs，it is composed of four air holes with a radius ofr2which are 30 degrees from the vertical center line.The remaining six air holes are left-right symmetrical structure，in which the abscissa of the three air holes on the right isthe distance between the three adjacent air holes on the right isΛ，and the remaining three air holes on the left are symmetrically distributed with the three air holes on the right.The outermost air holes are obtained by rotating the air holes with a radius ofr3with an abscissa of 2.5Λand an ordinate of 0.5Λfor three times，and each rotation angle is 30°.Then four air holes will be obtained with a radius ofr3，while the four air holes in the outermost layer on the left are symmetrically distributed with the four air holes with a radius ofr3on the right.The air hole structure of the proposed sensor is symmetrical so as to enlarge the energy leakage from the fiber core to enhance the coupling between optical and plasmonic modes.

The D-shaped sensors can be fabricated by stack-and-draw method［24］and side-polishing technique［25］，the gold film can be coated on the inner surface of the micro-opening channel by chemical vapor method［26］.The fiber material is fused silica and the RI is determined by the Sellmeier equation［27］

whereλrepresents the operating wavelength.Furthermore，the relative dielectric constant of gold can be demonstrated by Drude-Lorentz model，which is expressed as follows［28］

whereε∞=9.75 is the permittivity of gold at high frequencies，ωp=1.36×1016is the plasma frequency of gold，andωc=1.45×1014is the scattering frequency of electrons，andωis the angular frequency of the incident light.

In order to prevent the light reflects from the PCF boundary，an artificial boundary condition known as the Perfectly Matched Layer（PML）［29］is added to absorb the radiation energy at the outermost layer.SPR sensor is based on the coupling between the core mode and the plasmonic mode，a loss peak occurs at resonance wavelength.The coupling strength between the fiber core mode and the plasmonic mode is described by the confinement loss of the core mode，which is expressed as［30］

whereαlossrepresents the confinement loss，Im(neff) is the imaginary part of the complex effective refractive index.

2 Results and discussion

Fig.2 shows the relationship between the real part of the effective refractive index of core mode and plasmonic mode and the variation of confinement loss at a wavelength ranging from 1 800 nm to 2 300 nm at RI=1.2，in which the green curve stands for core mode，while the black curve represents the plasmonic mode.Meanwhile，the relationship between the confinement loss of the core mode and wavelength is represented by the red curve.As shown in Fig.2，when the real part of the effective refractive index of the core mode and the plasmonic mode are equal，a strong SPR effect can be observed with a prominent peak in the loss spectrum.In other words，a maximum energy transfers from the core mode to the plasmonic mode at the resonance wavelength.

Fig.3 shows the energy transfer process between core and plasmonic mode when the wavelength changes from 1 800 nm to 2 300 nm.As it can be seen from the figure，when the wavelength is away from the optimal phase matching point，most of the energy is concentrated in the core，while only a small portion of the energy leaks into the sensing interface.Obviously，the maximum energy transfers at the resonant wavelength，which is caused by the phase matching condition of core mode and plasmonic mode.

In order to optimize the sensing performance，the air holes radii，the pitch and gold film thickness were investigated separately.The refractive index of the analyte is fixed at 1.2 and the structural parameters of the PCF are calculated by keeping the rest of the structural parameters constant.The initial structural parametersr1，r2，r3，rsof the PCF are 0.5 μm，0.6 μm，0.7 μm and 1 μm，respectively.The innermost air hole pitch is 2.4 μm and the thickness of the gold film is 50 nm.It is worth noting that the variation of air hole diameter affects the coupling intensity and phase matching condition of core mode and plasmonic mode，thereby also affecting the energy conversion between them.

It can be seen from Fig.4（a）that the resonance peak gradually decreases from 229.43 dB/cm to 203.49 dB/cm asr1changes from 0.5 μm to 0.54 μm，and the resonance wavelength shifts from 2 070 nm to a longer wavelength of 2 150 nm.This can be attributed to the fact that more energy is confined to the core whenr1changes from 0.5 μm to 0.54 μm，which affects the coupling between the core mode and the plasmonic mode.Fig.4（b）describes the confinement loss decrease when radiusr2changes from 0.6 μm to 0.8 μm，and the corresponding blue shift occurs with the resonance peaks moving toward a shorter wavelength from 2 100 nm to 1 650 nm over the process.The reason is that the increase ofr2will increase the refractive index difference between the plasmonic mode and core mode，which will affect the coupling between them.Therefore，with the increase ofr2，the shorter wavelength can excite the plasmonic mode，resulting in the phenomenon of wavelength blue shift in the loss spectrum.

Fig.4（c）describes the relationship between the confinement loss and the radiusr3.As can be seen from the figure，the variation ofr3has almost no effect on the resonance wavelength and the confinement loss，which means the air holer3has a slight impact on the coupling between the core mode and the plasmonic mode.So the precision requirement of the radiusr3is low and therefore the fabrication difficulty of the sensor can be greatly reduced.Fig.4（d）shows the effect of the different open-ring channel radiirson the confinement loss of the core mode.It can be seen that the confinement loss reduces from 229.43 dB/cm to 54.21 dB/cm asrsincreases slightly from 1 μm to 1.2 μm and the resonance wavelength decreases from 2 070 nm to 1 550 nm.Therefore，the performance of the proposed sensor can be optimized by adjusting the size of the open-ring channel.In this case，we setr1=0.5 μm，r2=0.6 μm，r2=0.7 μm，andrs=1 μm，respectively，for the proposed sensor.

Fig.5（a）depicts the variation of the confinement loss of the proposed sensor with different pitches ranging from 2.4 μm to 2.56 μm.It can be seen that whenΛincreases，the resonance wavelength gradually decreasing from 2 070 nm to 1 870 nm with the resonance peak decreases from 229.40 dB/cm to 96.76 dB/cm.The reason is that the change ofΛhas an impact on the refractive index of core mode and plasmonic mode，which in turn affects the phase matching condition and energy coupling between them.Meanwhile，if the pitch is too large，it will prevent the concentration of the energy in the core.Totally，the pitch of 3 μm is selected to prevent the gap from increasing.

The thickness of gold film plays a vital role in the sensing performance.If the gold film is too thick，the electric field can not penetrate the gold film，which will reduce the sensitivity of the proposed sensor.While if the gold film is too thin，the plasmonic wave will be strongly suppressed due to radiation damping.Therefore，the thickness of gold film can significantly affect the coupling between the core mode and the plasmonic mode.Fig.5（b）exhibits the loss spectra variation as the thickness of gold film ranges from 40 nm to 60 nm.As shown in Fig.5（b），the reduction of the confinement loss and the blue shift appears as the thickness of the gold film increasing.A similar phenomenon was also noted in the paper［31］.Therefore，the sensitivity of the sensor can be improved by optimizing the thickness of the gold film.

The amplitude sensitivity is a critical parameter for evaluating the quality of the proposed sensor，and it is defined by the following equation［32］

where ?α(λ，na)represents the difference of the confinement loss of the adjacent analyte RI，α(λ，na)denotes the overall confinement loss，and ?narepresents difference of the analyte RI.

The resonance wavelength is used to detect the refractive index of the analyte，which means the analyte has a significant effect on the coupling intensity between the core mode and the plasmonic mode.The wavelength dependence of the loss spectra and the amplitude sensitivity of the sensor for analytes RI ranging from 1.18 to 1.30 has been shown in Fig.6（a）and Fig.6（b），respectively.The resonance wavelength shifts to a longer wavelength and the peak of the confinement loss also increases gradually when the analyte RI increases from 1.18 to 1.30.This can be attributed to the incomplete coupling between the core mode and the plasmonic mode.As shown in Fig.6（a），with the RI increasing，the confinement loss increases from 209.98 dB/cm to 579.23 dB/cm.We can notice that the maximum amplitude sensitivity is 68.13 RIU-1from Fig.6（b）.

Fig.7 illustrates the fitted curve of the resonance wavelength with the analyte RI changing from 1.18 to 1.30.The inserted table shows the results of the polynomial fitting，the slope of the fitted curve represents the sensitivity of the proposed sensor.The adjusted R-square value of the fitting curve is 0.995 08，which means a high agreement of fitting.It also can be seen that the resonant point moves to longer wavelengths as the analyte RI increases，with a maximum spectral sensitivity of 388 00 nm/RIU at the analyte RI of 1.30.

The spectral sensitivity has equal significance to the amplitude sensitivity at estimating the performance of the PCF sensor.The spectral sensitivity is defined by the following equation［33］

where Δλpeakrepresents the difference of the adjacent resonance wavelength.

In addition to the spectral and amplitude sensitivities，the resolution is a critical parameter of the proposed sensor.The resolution is expressed as［34］

where Δλminis the minimum spectral resolution.A high resolution of 2.37×10-6RIU can be gained at 1.30 in the proposed structure as shown in Table 1.

Table 1 shows the performance of the proposed sensor compared to previous sensors.It can be seen that the sensor proposed in this work has the advantage of higher refractive index detection.Meanwhile，the detection range of the proposed sensor is wider than that of the previous sensor.

Table 1 Comparison of the proposed sensor with previous designs

3 Conclusion

In this work，a novel D-shaped PCF sensor based on SPR effect is proposed and numerically investigated.Unlike conventional D-shaped structures，the SPR effect is excited by coating in the open-ring channel with a thin gold film，which can enhance the resonance effect and improve the spectral sensitivity of proposed sensor.In addition，the usage of this channel reduces the coating area of gold film.This makes the fabrication process more simplified and cost-effective.Simulation results show that the sensor can detect low refractive indexes ranging from 1.18 to 1.30 with the highest spectral sensitivity and resolution of 38 800 nm/RIU and 2.37×10-6RIU.Due to the low refractive index sensing range and the simple compact fiber design，the proposed SPR sensor could be a competitive candidate in low refractive index detection.