An Integrated Detection Based on a Multi-Parameter Plasmonic Optical Fiber Sensor

In this paper, a multi-parameter integrated detection photonic crystal fiber (PCF) sensor based on surface plasmon resonance (SPR) is proposed for its application in detecting temperature, magnetic field, and refractive index. The air holes on both sides of the fiber core were coated with gold film and introduced to the temperature-sensitive medium (PDMS) and magnetic fluid (MF), detecting temperature and magnetic field, respectively. The graphene layer is also presented on the gold film of the D-type side polished surface to improve the sensor sensitivity. The sensor’s critical parameters’ influence on its performance is investigated using a mode solver based on the finite element method (FEM). Simulation results show when the samples refractive index (RI) detection is a range of 1.36~1.43, magnetic field detection is a range of 20~550 Oe, and the temperature detection is a range of 5~55 °C; the maximum sensor’s sensitivity obtains 76,000 nm/RIU, magnetic field intensity sensitivity produces 164.06 pm/Oe, and temperature sensitivity obtains −5001.31 pm/°C.


Introduction
Compared with traditional optical fiber, photonic crystal fiber (PCF) has attracted more and more researchers' attention due to its high birefringence, low loss, high nonlinearity, and flexible structure. Photonic crystal fiber (PCF) introduces a stomatal filling structure into the cladding. It fills the pores with different sensitive materials to show photonic crystal fiber sensors' other sensing characteristics in different environments [1][2][3]. Surface plasmon resonance (SPR) is a physical phenomenon occurring on dielectric and metal surfaces [4][5][6].When the evanescent wave generated by the total internal reflection of incident light appears at the interface of the metal medium and the surface plasmonic wave meeting the phase-matching condition, resonance will occur, resulting in a large amount of incident light energy being coupled to the surface plasmon wave, the incident light energy will sharply decrease, and show a loss resonance peak in the loss spectrum. Because SPR can effectively detect small changes in the refractive index of the surrounding environment, in 1968, Otto combined a prism with a metal film to design the first SPR sensor structure [7], which attracted much attention because of its high sensitivity. However, the design of the Otto structure increases the complexity of the overall system of the sensor. Therefore, in 1971, Kretschmann improved the Otto structure [8] and significantly reduced the structural by finite element method. The multi-parameter sensor proposed in this paper is shown in Figure 1. Firstly, the gold film was coated on the PCF cladding polishing plane to form the SPR sensing channel. When the refractive index of the measured liquid changed, the wavelength shift of the SPR loss peak was analyzed to realize the sensor's refractive index sensing measurement. Then, the sensitivity of the sensor was improved by introducing a graphene layer on the gold film. The air hole on the left side of the fiber core was filled with magnetic fluid Fe 2 O 3 to form channel 1, and that on the right side of the center was filled with the temperature-sensitive medium PDMS to form channel 2. Channel 1 and channel 2 were symmetrical, and their inner walls were coated with a gold film. They used the magnetic fluid and the temperature-sensitive medium to respond differently to magnetic field changes and temperature changes. The peak values of SPR loss and the wavelength shift of the two channels' resonance wavelength were different. Then, through the sensitivity matrix method, the simultaneous sensing measurement of temperature and a magnetic field was realized; Finally, a D-type photonic crystal fiber SPR sensor was designed and acknowledged, which integrated multi-parameter detection of temperature, magnetic field, and refractive index. This paper analyzes the modal characteristics of an SPR D-type PCF sensor with multi-parameter integrated detection of temperature, magnetic field intensity, and refractive index by finite element method. The multi-parameter sensor proposed in this paper is shown in Figure 1. Firstly, the gold film was coated on the PCF cladding polishing plane to form the SPR sensing channel. When the refractive index of the measured liquid changed, the wavelength shift of the SPR loss peak was analyzed to realize the sensor's refractive index sensing measurement. Then, the sensitivity of the sensor was improved by introducing a graphene layer on the gold film. The air hole on the left side of the fiber core was filled with magnetic fluid Fe2O3 to form channel 1, and that on the right side of the center was filled with the temperature-sensitive medium PDMS to form channel 2. Channel 1 and channel 2 were symmetrical, and their inner walls were coated with a gold film. They used the magnetic fluid and the temperature-sensitive medium to respond differently to magnetic field changes and temperature changes. The peak values of SPR loss and the wavelength shift of the two channels' resonance wavelength were different. Then, through the sensitivity matrix method, the simultaneous sensing measurement of temperature and a magnetic field was realized; Finally, a D-type photonic crystal fiber SPR sensor was designed and acknowledged, which integrated multi-parameter detection of temperature, magnetic field, and refractive index. As shown in Figure 1, the air holes in the optical fiber are composed of 19 air holes arranged in parallel in four layers, in which the air holes in the first layer are arranged in a horizontal straight line with the fiber core and two-channel, and two air holes are distributed on the left and right; there are six air holes in the second layer, five in the third layer, and four in the fourth layer. The outermost air hole center of these four layers of air holes is connected to form half of a regular hexagon. By introducing air holes into the fiber cladding to create a solid photonic crystal fiber, the light wave can be better constrained in the core to gather high field strength in the center of the fiber and finally enhance the nonlinear effect of the fiber.
In the PCF sensor designed in this paper, the substrate material of optical fiber was molten SiO2, and the relationship between refractive index (RI) n and wavelength λ can be defined according to the Sellmeier Equation [17]. The expression is as follows: Figure 1. Schematic diagram of a plasmonic optical fiber sensor: (a) its three-dimensional structure; (b) its two-dimensional cross-section structure.
As shown in Figure 1, the air holes in the optical fiber are composed of 19 air holes arranged in parallel in four layers, in which the air holes in the first layer are arranged in a horizontal straight line with the fiber core and two-channel, and two air holes are distributed on the left and right; there are six air holes in the second layer, five in the third layer, and four in the fourth layer. The outermost air hole center of these four layers of air holes is connected to form half of a regular hexagon. By introducing air holes into the fiber cladding to create a solid photonic crystal fiber, the light wave can be better constrained in the core to gather high field strength in the center of the fiber and finally enhance the nonlinear effect of the fiber.
In the PCF sensor designed in this paper, the substrate material of optical fiber was molten SiO 2 , and the relationship between refractive index (RI) n and wavelength λ can be defined according to the Sellmeier Equation [17]. The expression is as follows: The relative dielectric constant of gold film Au is defined by the Drude model [18]. For the graphene layer, the refractive index expression [19] is as follows: Generally speaking, the thickness of the monolayer graphene is 0.34 nm. For Equation (2), C 1 is constant. Its value is 5.446 µm −1 , and λ is the vacuum wavelength. In this design, the total thickness of graphene was D = 0.34 × N, where N is the number of graphene layers.
Magnetic fluid (MF) is a kind of colloid solution. The changes of external temperature and magnetic field intensity will affect the refractive index of MF [20][21][22]. The relationship between them can be described by the Langevin Function [23]. The magnetic fluid material selected in this paper was water-based Fe 3 O 4 . When the volume ratio of Fe 3 O 4 to water is 3%, the refractive index of MF is 1.3592, and the thermo-optic coefficient and magnetooptic coefficient are −2.4 × 10 −4 / • C and 4.98 × 10 −5 /Oe, respectively. The function [24] of the magnetic fluid designed in this paper can be described as follows: PDMS (Polydimethylsiloxane) is a new type of polymer material. In existing temperature sensors, ethanol is mostly used as a sensitive material to detect the temperature. Compared with the PDMS used in this paper, ethanol is toxic, easy to leak and volatile, and PDMS has good mechanical properties, easy processing and higher thermooptic coefficient.When the temperature is 20 • C, the thermo-optic coefficient of PDMS is −4.5 × 10 −4 / • C, and the relationship between its refractive index (n PDMS ) and temperature (T) [25] can be expressed as follows: When the transmission constant of the SPP mode and the fundamental mode satisfies the matching condition, the two modes are coupled to excite the SPR. However, the effective refractive index of the SPP mode is greatly affected by the change of refractive index of the liquid to be measured, and the effective refractive index of the core mode is weakly affected, which leads to the resonance wavelength shifts when the two modes meet the wave vector matching. Therefore, we can detect the change in the effective refractive index of the liquid by measuring the offset of the resonant wavelength. The sensitivity S(λ) is one of the important parameters to evaluate the performance of the PCF sensor, and its definition [26] is as follows: Among them, ∆λ peak is the offset of the resonant wavelength, and ∆n s is the change of the refractive index of the liquid to be measured.
However, for PCF sensors, when the sensitivity is too high, it may lead to an increase in internal loss, resulting in too high a full width at half maxima (FWMH) value of the resonant peak, resulting in a decrease in the resolution of the sensor. Therefore, we define the quality factor figure of merit (FOM) to measure the characteristics of the transmission curve, which is as follows: Among them, sensor resolution R = ∆n s × ∆λ min /∆λ peak . In this design, channel 1 and channel 2 were filled with different sensitive materials. When the external magnetic field or temperature changes, the refractive index of the sensitive materials in the two channels also changes; thus, the resonance coupling phenomenon between the SPP mode and the fiber core mode changes, resulting in the resonance wavelength shift; that is, integrated detection of magnetic field and temperature can be realized by measuring the resonant wavelength offset of channel 1 and channel 2, respectively. In the wavelength measurement method, the calculation formula of sensitivity is as follows: K ch (T) = ∆λ ch /∆T K ch (H) = ∆λ ch /∆H (8) Therefore, the relationship between resonant wavelength offset, temperature change, and magnetic field intensity change is expressed by a two-parameter sensitivity matrix as follows: Among them, ∆λ i,j , i = 1,2, j = T, H denotes the shift of the resonant wavelength caused by the change of temperature or magnetic field intensity in two channels. ∆T represents the change in external temperature, and ∆H represents the change in the intensity of the external magnetic field, K chi (j), i = 1, 2, j = T, H, denotes the temperature sensitivity and magnetic field sensitivity of the two channels, respectively.
The sensitivity matrix in Equation (9) is inverted. That is, the sensing matrix for calculating the change of temperature and magnetic field can be obtained, and the expression is as follows: From Equation (10), we can calculate the change of temperature and magnetic field by measuring the resonant wavelength's offset and realizing the simultaneous detection of temperature and magnetic field.

Results and Discussion
In this paper, using the full vector finite element method, the solution region was divided into finite elements. The process of mathematical approximation numerically simulated the physical structure designed in this paper. In the simulation, the free triangle meshing method was adopted. The frame was divided into 74,803 domain units, of which the minimum element mass was 0.4486. The numerical analysis of a surface plasmon resonance D-type PCF sensor integrated with temperature, magnetic field intensity, and refractive index was carried out under the boundary condition of a perfectly matched layer (PML). By changing the test environment and using wavelength modulation, the drift of the resonant wavelength was measured, and the sensor function of multi-parameter integrated detection was realized.

Sensor Structure Optimization
To obtain better sensor performance, then we optimized the structure of the sensor. First of all, under the condition of keeping other parameters unchanged, we used the transition boundary conditions to optimize the structure of photonic crystal fiber, and finally, used OriginPro software to process the data. Next, the sensor's design was optimized, including the air hole spacing Λ, the side polishing depth d h , the air hole diameter d, and the diameter d ch1 and d ch2 of channel 1 and channel 2.
First of all, as can be seen from Figure 2a, the air hole spacing Λ changed from 15 µm to 17 µm and the loss decreased gradually. This is because the increase in the air hole spacing meant that the lattice period of the PCF became larger. The fiber had a larger core size, which made the core more restrained to light, more energy was constrained in the core, and the loss gradually became smaller. However, the change of air hole spacing did not affect the resonant wavelength of the sensor, but it caused the FWHM to increase at first and then decrease. Therefore, after a comprehensive analysis, we chose the air hole spacing Λ = 16 µm.
Then, the side depth d h was optimized, and the simulation results are shown in Figure 2b. As shown in Figure 2b, when the side polishing depth d h changed from 0.6 Λ to 0.8 Λ, the loss peak value decreased gradually, and the resonance wavelength was red-shifted. This is because, with the increase in d h , the distance between the core base mode and the SPW mode increases, resulting in a decrease in the coupling strength with the SPW mode and a reduction in the energy transferred from the core to the SPW mode, thus showing a peak decrease in the loss peak spectrum. Although the change of d h will lead to the shift of resonance wavelength, the change curve of resonance wavelength is the same in different refractive index environments; therefore, the change of d h does not affect the sensitivity of the sensor, but it gradually reduces the FWHM of the sensor. Therefore, based on the above analysis, we finally chose the side depth d h = 0.7 Λ.
Sensors 2021, 21, x FOR PEER REVIEW 6 of 14 core, and the loss gradually became smaller. However, the change of air hole spacing did not affect the resonant wavelength of the sensor, but it caused the FWHM to increase at first and then decrease. Therefore, after a comprehensive analysis, we chose the air hole spacing Λ = 16 µm.  Figure 2b, when the side polishing depth dh changed from 0.6 Λ to 0.8 Λ, the loss peak value decreased gradually, and the resonance wavelength was redshifted. This is because, with the increase in dh, the distance between the core base mode and the SPW mode increases, resulting in a decrease in the coupling strength with the SPW mode and a reduction in the energy transferred from the core to the SPW mode, thus showing a peak decrease in the loss peak spectrum. Although the change of dh will lead to the shift of resonance wavelength, the change curve of resonance wavelength is the same in different refractive index environments; therefore, the change of dh does not affect the sensitivity of the sensor, but it gradually reduces the FWHM of the sensor. Therefore, based on the above analysis, we finally chose the side depth dh = 0.7 Λ.
Then, we optimized the air hole diameter of the sensor, and the simulation results are shown in Figure 2c. The simulation results show that when the air holes diameter d changed from 0.3 Λ to 0.5 Λ, the loss of the sensor increased gradually. The reason is opposite to the reasons for the change in Figure 2a; the increase in the diameter of the air hole will lead to the smaller size of the fiber core, which will lead to the deterioration of the constraining ability of the fiber core to light. More energy will be transferred from the fiber core to SPW mode during resonance, and the loss will gradually increase.
Similarly, we can observe that the size of the air hole has no effect on the position of the resonant wavelength because the change of the air hole size has little impact on the core mode, so the phase-matching condition of SPR between the core mode and the SPW mode will not be affected. Therefore, we chose the diameter of the air hole as 0.5 Λ. Moreover, because the air hole spacing Λ, the side polishing depth dh, and the air hole diameter d do not affect the sensor's sensitivity, the proposed sensor in this paper had a good fabrication tolerance.
Next, we continued to optimize the diameter of channel 1 and channel 2 of the sensor, and the simulation results are shown in Figure 3. From Figure 3a,d, we can see that when the diameter dch1 and dch2 increase, the loss of channel 1 and channel 2 increase gradually but have little effect on the RI. This is because when the diameter of channel 1 or channel 2 become larger, the defect hole diameter increases, making the distribution of the refractive index of the fiber more asymmetric and the birefringence of PCF stronger. It will also  Then, we optimized the air hole diameter of the sensor, and the simulation results are shown in Figure 2c. The simulation results show that when the air holes diameter d changed from 0.3 Λ to 0.5 Λ, the loss of the sensor increased gradually. The reason is opposite to the reasons for the change in Figure 2a; the increase in the diameter of the air hole will lead to the smaller size of the fiber core, which will lead to the deterioration of the constraining ability of the fiber core to light. More energy will be transferred from the fiber core to SPW mode during resonance, and the loss will gradually increase.
Similarly, we can observe that the size of the air hole has no effect on the position of the resonant wavelength because the change of the air hole size has little impact on the core mode, so the phase-matching condition of SPR between the core mode and the SPW mode will not be affected. Therefore, we chose the diameter of the air hole as 0.5 Λ. Moreover, because the air hole spacing Λ, the side polishing depth d h , and the air hole diameter d do not affect the sensor's sensitivity, the proposed sensor in this paper had a good fabrication tolerance.
Next, we continued to optimize the diameter of channel 1 and channel 2 of the sensor, and the simulation results are shown in Figure 3. From Figure 3a,d, we can see that when the diameter d ch1 and d ch2 increase, the loss of channel 1 and channel 2 increase gradually but have little effect on the RI. This is because when the diameter of channel 1 or channel 2 become larger, the defect hole diameter increases, making the distribution of the refractive index of the fiber more asymmetric and the birefringence of PCF stronger. It will also weaken the energy constraint ability of the fiber core in the x-axis direction, and more x-direction means more fundamental mode energy, which is coupled to channel 1 and channel 2, which increases the loss of channel 1 and channel 2. Simultaneously, the increase in channel 1 or channel 2 will improve the magnetic fluid content, and PDMS filled in channel 1 and channel 2 coupled more core energy to the channels on both sides. However, because the refractive index of MF is much smaller than that of the core, the energy coupled to channel 1 will be less than that of channel 2. At the same time, we can also find that the size changes of channel 1 and channel 2 do not affect the RI, which indicates that the size changes of the two channels have no effect on the base mode energy in the y-direction; when the diameter of the two channels increases, the metal film coated on the two channels will not only increase the contact area but also reduce the sensing distance from the fiber core diameter. However, the refractive index of MF is much smaller than that of the fiber core, so the shift of the resonant wavelength is not apparent, and the resonant wavelength is unchanged. However, the effect of channel 2 is more prominent. As shown in Figure 3d, the resonance wavelength will find a redshift, but under different diameters, the shifting trend of the resonance wavelength is the same; it has no effect on the sensor's sensitivity. Additionally, from Figure 3b, it can be found that the increase in channel 1 size will make the FWHM decrease at first and then increase; Figure 3e shows that with the rise of channel 2 diameter, the FWHM will increase. When the magnetic field intensity changes, we can find from Figure 3c that the resonant wavelength of channel 1 changes, while RI and channel 2 are not affected. However, when the temperature changes, although channel 1 is mainly responsible for magnetic field detection because its filled magnetic fluid also has a temperature sensitivity coefficient, it will also have an impact on channel 1, as shown in Figure 3f. Similarly, temperature changes do not affect RI. This can reflect the multiparameter integrated detection sensor we designed. There is no cross-influence among the parameters. Because d ch1 and d ch2 do not affect the sensor's sensitivity, the sensor proposed in this paper has a good fabrication tolerance. nel 1 and channel 2 coupled more core energy to the channels on both sides. However, because the refractive index of MF is much smaller than that of the core, the energy coupled to channel 1 will be less than that of channel 2. At the same time, we can also find that the size changes of channel 1 and channel 2 do not affect the RI, which indicates that the size changes of the two channels have no effect on the base mode energy in the ydirection; when the diameter of the two channels increases, the metal film coated on the two channels will not only increase the contact area but also reduce the sensing distance from the fiber core diameter. However, the refractive index of MF is much smaller than that of the fiber core, so the shift of the resonant wavelength is not apparent, and the resonant wavelength is unchanged. However, the effect of channel 2 is more prominent. As shown in Figure 3d, the resonance wavelength will find a redshift, but under different diameters, the shifting trend of the resonance wavelength is the same; it has no effect on the sensor's sensitivity. Additionally, from Figure 3b, it can be found that the increase in channel 1 size will make the FWHM decrease at first and then increase; Figure 3e shows that with the rise of channel 2 diameter, the FWHM will increase. When the magnetic field intensity changes, we can find from Figure 3c that the resonant wavelength of channel 1 changes, while RI and channel 2 are not affected. However, when the temperature changes, although channel 1 is mainly responsible for magnetic field detection because its filled magnetic fluid also has a temperature sensitivity coefficient, it will also have an impact on channel 1, as shown in Figure 3f. Similarly, temperature changes do not affect RI. This can reflect the multi-parameter integrated detection sensor we designed. There is no cross-influence among the parameters. Because dch1 and dch2 do not affect the sensor's sensitivity, the sensor proposed in this paper has a good fabrication tolerance.

Graphene Layer Number
Then, we used the transition boundary conditions to optimize the number of graphene layers while keeping other parameters unchanged. Finally, we used OriginPro software to process the data, and the results are shown in Figure 4. As shown in Figure 4a,b, in the same external environment, with the increase in graphene thickness, the resonance wavelength increases, the loss decreases, and the resonance curve widens, and FWHM increases. This is mainly because the effective refractive index of the metal interface increases with the increase in graphene thickness, so a higher resonance wavelength is needed to meet the resonance condition. At this time, the damping loss of graphene increases. Therefore, when SPR occurs, the core mode's energy coupled to the SPP film is reduced, resulting in a reduced loss. From Figure 4c, we can see that when the refractive index of the external environment is in the range of 1.36 to 1.42, the average resonance wavelength shift increases with the increase in graphene thickness, so the average sensitivity of the sensor increases. When the refractive index of the external environment is in the range of 1.36 to 1.43, the average resonant wavelength shift fluctuates slightly with the rise of graphene thickness, reaches the maximum at N = 3, and then decreases gradually. This is because graphene is a kind of lossy medium with real and imaginary parts. When the number of layers of graphene (N) increases, the imaginary part of graphene also increases exponentially, which increases the carrier concentration and leads to the increase in the chemical potential (µ c ) of graphene, and the corresponding resonance wavelength is red-shifted [27]. After a comprehensive analysis, we chose the layer number of graphene N = 3.
Sensors 2021, 21, x FOR PEER REVIEW 9 of 14 the refractive index changes, channel 1 and channel 2 will not change, indicating that the RI channel will not have cross-effects with the other two channels.

Gold Film Thickness
The thickness of the gold film has a significant influence on the PCF sensor's sensing characteristics based on SPR. This is mainly because the energy of the evanescent wave excited by the fiber core attenuates exponentially during propagation, so the thickness of the gold film is smaller than the penetration depth of an evanescent wave d p . Among them, the formula of penetration depth d p [28] is as follows: Among them, n 1 is the refractive index of the fiber core, n aq is the refractive index of the surrounding medium, and θ int is the incident angle. This paper is mainly designed in the visible and near-mid-infrared frequency band. As can be determined from Equation (11), Sensors 2021, 21, 803 9 of 14 the thickness of the metal film of this design is at least less than 60.5 nm. This design optimizes the thickness of a gold film from 25~60 nm, which corresponds with the propagation of an evanescent wave.
We simulated the gold film thickness of the side polishing plane, channel 1, and channel 2. Firstly, the structure of the gold film thickness t coated on the side polishing plane was optimized, and the result is shown in Figure 5. As can be seen from Figure 5a, when the thickness t of the gold film becomes more extensive, the loss of the sensor decreases gradually, the resonance wavelength is red-shifted, and the FWHM decreases slowly. This is mainly because the thicker the gold film is, the more the energy of evanescent wave attenuates and less energy is coupled with the SPW mode, so the loss becomes smaller. When the thickness of the gold film increases, the effective refractive index of SPW increases, but the core remains unchanged, which the phase matching point of SPW mode and the core fundamental mode redshifts when SPR occurs, that is, the resonance wavelength redshifts. In this case, FWHM decreases with the increase in the gold film thickness because when the gold film is relatively thin, the evanescent field's energy focused on the film is very weak, and the excitation effect of SPR is not good, thus making the FWHM larger. Besides, it can be found that when the thickness of the gold film changes or when the refractive index changes, channel 1 and channel 2 will not change, indicating that the RI channel will not have cross-effects with the other two channels.
RI channel will not have cross-effects with the other two channels.   When the refractive index n s changes successively from 1.36 to 1.43 in turn by step length 0.01, and the thickness of the gold film varies from 40 nm to 60 nm in turn by step length 5 nm, from Figure 5b, we know that the average resonant wavelength shift and maximum resonant wavelength shift gradually increase. Based on the above analysis, the gold film thickness t = 55 nm was selected in this paper. In this case, the average resonant wavelength offset is 175.71 nm, and the maximum resonant wavelength offset is 760 nm. According to Formula (7), the average sensitivity of the sensor designed in this paper is 17,571 nm/RIU. Maximum sensitivity can reach 76,000 nm/RIU. When the minimum resolution of the spectrometer is ∆λ min = 0.1 nm, the average minimum spectral resolution is 5.69 × 10 −6 (RIU). From Figure 5c, the minimum FWHM of RI can reach 26.86 nm at t = 55 nm. According to Formula (6), the FOM corresponding to the average sensitivity and the maximum sensitivity are 654.17 RIU −1 and 2829.49 RIU −1 .
Next, the structure of the gold film thickness t 1 of channel 1 was optimized, and the result is shown in Figure 6. It can be seen from Figure 6a,b, when the thickness of the gold film t 1 becomes larger, channel 1 affects the x-direction fundamental mode of the fiber core, so that the loss of the sensor decreases gradually, the resonance wavelength is redshifted, and the FWHM decreases slowly. However, the increase in t 1 does not affect the symmetrical channel 2 in the x-direction. This is because the change of the thickness of the gold film in the two channels has little effect on the energy constraint ability of the core in the x-axis direction. With the increase in magnetic field intensity, the magnetic particles of MF gather, and the effective refractive index increases, which leads to the change of phase matching point during resonance coupling, the red shift of resonance wavelength, the better effect of SPR excitation, and the decrease in FWHM. At this point, from Figure 6a,b, we can see that channel 2 and the RI channel have not changed, indicating that channel 1 will not have a cross effect with the RI channel and channel 2. Under different magnetic field ranges, when the gold film thickness t 1 changes sequentially from 25 nm to 50 nm in turn by step length 5 nm, the magnetic field sensitivity tends to increase by observing the Figure 6c. Because when the gold film is too thin, the SPR excitation effect is not good, the FWHM is too high, the gold film is too thick, and the loss peak is too low to be observed. With comprehensive consideration, the designers chose channel 1 gold film thickness t 1 = 40 nm. maximum resonant wavelength shift gradually increase. Based on the above analysis, the gold film thickness t = 55 nm was selected in this paper. In this case, the average resonant wavelength offset is 175.71 nm, and the maximum resonant wavelength offset is 760 nm. According to Formula (7), the average sensitivity of the sensor designed in this paper is 17,571 nm/RIU. Maximum sensitivity can reach 76,000 nm/RIU. When the minimum resolution of the spectrometer is Δλmin = 0.1 nm, the average minimum spectral resolution is 5.69 × 10 −6 (RIU). From Figure 5c, the minimum FWHM of RI can reach 26.86 nm at t = 55 nm. According to Formula (6), the FOM corresponding to the average sensitivity and the maximum sensitivity are 654.17 RIU −1 and 2829.49 RIU −1 .
Next, the structure of the gold film thickness t1 of channel 1 was optimized, and the result is shown in Figure 6. It can be seen from Figure 6a,b, when the thickness of the gold film t1 becomes larger, channel 1 affects the x-direction fundamental mode of the fiber core, so that the loss of the sensor decreases gradually, the resonance wavelength is redshifted, and the FWHM decreases slowly. However, the increase in t1 does not affect the symmetrical channel 2 in the x-direction. This is because the change of the thickness of the gold film in the two channels has little effect on the energy constraint ability of the core in the x-axis direction. With the increase in magnetic field intensity, the magnetic particles of MF gather, and the effective refractive index increases, which leads to the change of phase matching point during resonance coupling, the red shift of resonance wavelength, the better effect of SPR excitation, and the decrease in FWHM. At this point, from Figure  6a,b, we can see that channel 2 and the RI channel have not changed, indicating that channel 1 will not have a cross effect with the RI channel and channel 2. Under different magnetic field ranges, when the gold film thickness t1 changes sequentially from 25 nm to 50 nm in turn by step length 5 nm, the magnetic field sensitivity tends to increase by observing the Figure 6c. Because when the gold film is too thin, the SPR excitation effect is not good, the FWHM is too high, the gold film is too thick, and the loss peak is too low to be observed. With comprehensive consideration, the designers chose channel 1 gold film thickness t1 = 40 nm. When the magnetic field intensity varies at 20~550 Oe, Figure 6d shows that the resonant wavelength shift increases with the magnetic field intensity offset. At this time, through the linear fitting of OriginPro software, we can obtain ∆λ 1H = 0.15631∆H − 4.89956, the sensitivity of magnetic field intensity was 156.31 pm/Oe, and the degree of linearity was 0.99027; the wavelength resolution of the sensor was assumed to be ∆λmin = 0.1 nm, and the magnetic field resolution was 0.64 Oe [29]. Channel 2 was not affected by the change of magnetic field intensity. As shown in Figure 6e, when t 1 = 40 nm, the FWHM of channel 1 decreased with the increase in magnetic field intensity, and the minimum can reach 34.97 nm. In this case, according to Equation (6), FOM = 4.47 × 10 −3 Oe −1 .
Finally, the structure of the gold film thickness t 2 of channel 2 was optimized, and the result is shown in Figure 7. As shown in Figure 7a,b, when t 2 increases, channel 2 influences the primary mode in the x-direction of the fiber core, making the loss of channel 2 gradually decrease, the resonance wavelength red-shift, and a gradual increase in FWHM. Simultaneously, the rise of t 2 does not affect the symmetrical channel 1 in the x-direction. Figure 7a,b shows that channel 2 and RI channels have not changed, indicating that channel 1 will not have cross effects with RI channels and channel 2. However, the thickness of the gold film t 2 will affect the distance between the RI channel's resonant wavelengths and channel 2. When the space is too close, it will lead to the superposition of the two resonant wavelengths. With comprehensive consideration, the designers chose the channel 1 gold film thickness t 2 = 45 nm.
(e) The FWHM of channel 1 varies with the intensity of the magnetic field for t1 = 40 nm.
When the magnetic field intensity varies at 20~550 Oe, Figure 6d shows that the resonant wavelength shift increases with the magnetic field intensity offset. At this time, through the linear fitting of OriginPro software, we can obtain Δλ1H = 0.15631ΔH − 4.89956, the sensitivity of magnetic field intensity was 156.31 pm/Oe, and the degree of linearity was 0.99027; the wavelength resolution of the sensor was assumed to be Δλmin = 0.1 nm, and the magnetic field resolution was 0.64 Oe [29]. Channel 2 was not affected by the change of magnetic field intensity. As shown in Figure 6e, when t1 = 40 nm, the FWHM of channel 1 decreased with the increase in magnetic field intensity, and the minimum can reach 34.97 nm. In this case, according to Equation (6), FOM = 4.47 × 10 −3 Oe −1 .
Finally, the structure of the gold film thickness t2 of channel 2 was optimized, and the result is shown in Figure 7. As shown in Figure 7a,b, when t2 increases, channel 2 influences the primary mode in the x-direction of the fiber core, making the loss of channel 2 gradually decrease, the resonance wavelength red-shift, and a gradual increase in FWHM. Simultaneously, the rise of t2 does not affect the symmetrical channel 1 in the x-direction. Figure 7a,b shows that channel 2 and RI channels have not changed, indicating that channel 1 will not have cross effects with RI channels and channel 2. However, the thickness of the gold film t2 will affect the distance between the RI channel's resonant wavelengths and channel 2. When the space is too close, it will lead to the superposition of the two resonant wavelengths. With comprehensive consideration, the designers chose the channel 1 gold film thickness t2 = 45 nm.  The temperature T is changed sequentially from 5 • C to 55 • C, with a step size of 5 • C for t 2 = 45 nm. By observing Figure 7c, it can be found that both channel 1 and channel 2 have a blue shift in resonance wavelength and a gradual decrease in the loss with the increase in temperature. This is because both MF and PDMS of channel 1 and channel 2 have negative temperature coefficients. Figure 7d shows that the resonant wavelength offset of the two channels decreases with the increase in the temperature shift; after the linear fitting of OriginPro software, the linear fitting expression of channel 1 could be obtained as ∆λ 1T = −0.50909∆T + 1.45455, the temperature sensitivity was −509.09 pm/ • C, the degree of linearity was 0.95455, the linear fitting expression of channel 2 was ∆λ 2T = −5.00131∆T − 4.73636, the temperature sensitivity was 5001.31 pm/ • C, and the degree of linearity was 0.99275. The wavelength resolution of the sensor was assumed to be ∆λmin = 0.1 nm, and the temperature resolutions of channel 1 and channel 2 were 0.196 • C and 0.02 • C, respectively. From Figure 7e, when t 2 = 45 nm, the FWHM of channel 1 increased with the increase in temperature, the minimum FWHM was 40.43 nm, and the FWHM of channel 2 decreased with the rise of temperature, and the minimum could reach 38.84 nm. According to Equation (6), the FOM of channel 1 and channel 2 are −1.26 × 10 −2 • C −1 and −0.129 • C −1 , respectively.
Through the above analysis, we obtain the temperature sensitivity of the sensor channel 1 and channel 2 as −509.09 pm/ • C and −5001.31 pm/ • C, and magnetic field intensity sensitivity as 156.31 pm/Oe and 0 pm/Oe, respectively. According to Equation (9) and Equation (10) To summarize, the plasmonic optical fiber sensor proposed in this paper, which integrates temperature, magnetic field intensity, and refractive index multi-parameter detection, can not only detect these three parameters at the same time, but also has high refractive index sensitivity, temperature sensitivity, and magnetic field sensitivity.

Conclusions
In summary, a surface plasmon resonance D-type photonic crystal fiber (PCF) sensor integrated with multi-parameter detection of temperature, magnetic field intensity, and refractive index has been proposed. The theoretical model was simulated and analyzed by the finite element method. The D-type photonic crystal fiber SPR sensor for multiparameter integrated detection was realized using the magneto-optic effect of magnetic fluid and the temperature-sensitive development of SPR and PDMS. The advantages of this design are: it can effectively avoid the cross-sensitivity among the three parameters, realize the integrated detection of multi-parameters, overcome the limitation of single-parameter measurement of traditional sensors, and provide an idea for the field of multi-function sensor detection. The graphene layer was introduced to improve the RI channel's sensitivity so that the average sensitivity of the designed RI was 17,571 nm/RIU, and the maximum sensitivity was 76,000 nm/RIU, which was higher than that of some other published sensors [30][31][32]. This design has the advantages of a simple structure and principle, and high magnetic field sensitivity and temperature sensitivity in visible and near-mid-infrared bands. This design also has good production tolerance. Therefore, the sensor designed in this paper is a promising tool in the field of multi-sensor detection and biosensor.