<figure><img src="image_1.png" /></figure>

HAL
open science

# Optimizing Evanescent Efficiency of Chalcogenide Tapered Fiber

Xudong Zhao, Ni Yao, Xianghua Zhang, Lei Zhang, Guangming Tao, Zijian

Li, Quan Liu, Xiujian Zhao, Yinsheng Xu

► To cite this version:

Xudong Zhao, Ni Yao, Xianghua Zhang, Lei Zhang, Guangming Tao, et al. Optimizing Evanescent Efficiency of Chalcogenide Tapered Fiber. Materials, 2022, 15 (11), pp.3834. 10.3390/ma15113834. hal-03715188

HAL Id: hal-03715188

https://hal.science/hal-03715188

Submitted on 6 Jul 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Optimizing Evanescent Efficiency of Chalcogenide Tapered Fiber

Xudong Zhao<sup>1, ●</sup>, Ni Yao<sup>2,*</sup>, Xianghua Zhang<sup>1,3, ●</sup>, Lei Zhang<sup>2,4</sup>, Guangming Tao<sup>5</sup>, Zijian Li<sup>1</sup>, Quan Liu<sup>1</sup>,
Xiujian Zhao<sup>1</sup> and Yinsheng Xu<sup>1,*</sup>

**Abstract:** Evanescent wave absorption-based mid-infrared chalcogenide fiber sensors have prominent advantages in multicomponent liquid and gas detection. In this work, a new approach of tapered-fiber geometry optimization was proposed, and the evanescent efficiency was also theoretically calculated to evaluate sensing performance. The influence of fiber geometry (waist radius ($R_w$), taper length ($L_t$), waist deformation) on the mode distribution, light transmittance ($T$), evanescent proportion ($T_O$) and evanescent efficiency ($\tau$) is discussed. Remarkably, the calculated results show that the evanescent efficiency can be over 10% via optimizing the waist radius and taper length. Generally, a better sensing performance based on tapered fiber can be achieved if the proportion of the $LP_{11}$-like mode becomes higher or $R_w$ becomes smaller. Furthermore, the radius of the waist boundary ($R_L$) was introduced to analyze the waist deformation. Mode proportion is almost unchanged as the $R_L$ increases, while $\tau$ is halved. In addition, the larger the micro taper is, the easier the taper process is. Herein, a longer waist can be obtained, resulting in larger sensing area which increases sensitivity greatly.

**Keywords:** evanescent wave; tapered fiber; sensor; evanescent wave efficiency

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Zhao, X.; Yao, N.; Zhang, X.; Zhang, L.; Tao, G.; Li, Z.; Liu, Q.; Zhao, X.; Xu, Y. Optimizing Evanescent Efficiency of Chalcogenide Tapered Fiber.
*Materials* **2022**, *15*, 3834. https://doi.org/10.3390/ma15113834

Academic Editor: George Kenanakis

Received: 2 May 2022

Accepted: 24 May 2022

Published: 27 May 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2022 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

## 1. Introduction

Evanescent wave absorption-based fiber sensors have widespread applications ranging from biosensing to gas sensing. To date, various types of fiber sensors based on evanescent wave have been proposed, among which chalcogenide fiber possess unique advantages for its wide transmission window which covers the mid-infrared molecular fingerprint. In recent years, researchers from University of Rennes 1 [1–3], Ningbo University [4–9], Zhejiang University [10,11], and other institutions [12–24] have carried out extensive research on the mid-infrared evanescent wave absorption sensing technology, which proves that this technique can realize the real-time in situ qualitative and quantitative analysis of a variety of organic compounds. Dai et al. developed a gas sensor based on a four-hole suspended core $As_2S_3$ optical fiber, the sensitivity of which is less than 100 ppm for methane, and response time is estimated to be less than 20 s [8]. Romanova et al. prepared the $Ge_{26}As_{17}Se_{25}Te_{32}$ fiber loop for the detection of an anti-gel additive in diesel fuel. The absorption coefficients of the solution are well approximated by a linear function [13].

Evanescent wave absorption sensing, based on the resonance absorption between the organic groups and the electromagnetic wave of a specific wavelength, can realize the molecule identification and concentration detection. As shown in Figure 1, the light

<footnote n="1">¹ State Key Laboratory of Silicate Materials for Architectures, Wuhan University of Technology, Wuhan 430070, China; 291035@whut.edu.cn (X.Z.); xzhang@univ-rennes1.fr (X.Z.); zijian_li@whut.edu.cn (Z.L.); quan_liu@whut.edu.cn (Q.L.); opluse@whut.edu.cn (X.Z.)</footnote>

<footnote n="2">² Research Center for Intelligent Sensing, Zhejiang Laboratory, Hangzhou 311121, China; zhang_lei@zju.edu.cn</footnote>

<footnote n="3">³ Laboratoire des Verres et Céramiques, UMR-CNRS 6226, Sciences Chimiques de Rennes,
Université de Rennes 1, 35042 Rennes, France</footnote>

<footnote n="4">⁴ State Key Laboratory of Modern Optical Instrumentation, College of Optical Science and Engineering, Zhejiang University, Hangzhou 310027, China</footnote>

<footnote n="5">5 Wuhan National Laboratory for Optoelectronics, Huazhong University of Science and Technology, Wuhan 430074, China; tao@hust.edu.cn</footnote>

<footnote>* Correspondence: yaoni20@zhejianglab.com (N.Y.); xuyinsheng@whut.edu.cn (Y.X.)</footnote>
is propagated through the fiber via successive total internal reflection, among which $E = A\exp(-\kappa z)\exp[i(k_x x - \omega t)]$ is the formula of the evanescent wave. The $A\exp(-\kappa z)$ and $\exp[i(k_x x - \omega t)]$ are the amplitude and phase of the evanescent wave, respectively. When the amplitude of the evanescent wave decays to $1/e$ of its maximum value, this position is defined as penetration depth, that is, $z = 1/\kappa$. The imaginary part of the refractive index is introduced to characterize the light absorption ability of external medium. It should be noted that the molecule identification and concentration detection are realized by analyzing the infrared spectra, including the intrinsic absorption and intensity of molecules.

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Mathematical model for numerical simulation.</figcaption></figure>

How to improve the sensitivity of the fiber sensor has always been a research hotspot. The plasmonic-based fibers have high sensitivity in the visible region by using the surface plasmon resonance (SPR) effect of coating nanoparticles [22–24]. However, the SPR effect is not suitable for improving the sensitivity of infrared fibers. Geometry optimization is another way to improve the sensitivity of the fiber sensor. To date, varieties of the fiber geometry, such as tapered fiber [25], U-shaped fiber [26] and D-shaped fiber [27], have been widely studied, for these structures can strengthen the evanescent field greatly. Among them, the tapered fiber adjusts the mode distribution by reducing the core diameter to realize more evanescent wave energy existence in the boundary region. In the past years, waveguiding properties of tapered fiber have been extensively investigated theoretically, but mostly on adiabatic tapered silica fiber, which cannot transmit mid-infrared light [28–30]. Due to the great difference in optical properties between chalcogenide fiber and silica fiber, it is essential to carry out simulations for chalcogenide tapered fiber. To the best of our knowledge, no one has ever systematically studied the waveguiding properties of chalcogenide tapered fiber theoretically for sensing applications.

Here, we systematically investigate the evanescent field in different geometries of tapered fiber by using the finite-element-method-based full-vector solver COMSOL. Specifically, the electromagnetic wave frequency domain (ewfd) module is used for calculation. A general method for geometry design of chalcogenide tapered fiber is presented in this paper, which could be used for the tapered fiber optimization.
## 2. Basic Model

The mathematical model for our numerical simulation is illustrated in Figure 1, in which a tapered As₂Se₃ fiber with a radius of R<sub>w</sub> is embedded in glucose solution. For simplicity, the radius (R<sub>o</sub>) of the untapered fiber is fixed at 100 µm, and the adjustable radius (R<sub>w</sub>) of the waist and the length (L<sub>t</sub>) of the cone are schematically shown in Figure 1.

To achieve high sensitivity for glucose, the input light wavelength for detection is assumed to be 7 µm, which corresponds to the position of the characteristic peak of glucose, of which the absorption coefficient is large, and the loss of this chalcogenide fiber at 7 µm is small. Another reason for this selection is that the cut-off radius corresponding to the same mode is different for different wavelengths. Among them, the shorter the wavelength, the smaller the corresponding mode cut-off radius and the smaller the radius of the optimized structure. The refractive indices of As₂Se₃ and water at 7 µm are 2.70 and 1.33, respectively.

As shown in Figure 2, a potential fiber-puller system, which consists of an electrical heater, two motorized translation stages, a temperature controller, a CCD and a computer, can be used to draw the fiber taper. To achieve the flexible fiber geometry design, the modified resistance heating platform with different electrical heaters is controlled by a temperature controller. The chalcogenide fiber is fixed in fiber clamps and preheated before being pulled by the stepper motor controller. The diameter of the chalcogenide fiber is monitored in situ by CCD linked with a computer. Remarkably, the parameters chosen in this work have the potential to be realized by our homemade tapering platform.

<figure><img src="image_6.png" /><figcaption>Figure 2. Experimental setup for the chalcogenide fiber tapering.</figcaption></figure>

The different geometries are investigated to evaluate the transmission losses, and the transmittance (T) can be described as follows:

$$T = \sum (S_{21})^2 \qquad (1)$$

where S<sub>21</sub> is the energy-output ratio of different modes. Typically, for the adiabatic taper, T is approximately equal to 1.

In order to ensure the accuracy of the calculation, it is necessary to ensure the accuracy of the energy at the incident end (P<sub>in</sub>) before calculating the evanescent efficiency. According to the wave optics, the energy at the incident end can be obtained by taking the surface integral of the Poynting vector. However, since this is a two-dimensional model, it is necessary to determine the height of the two-dimensional model first, as shown in the following formula:

$$P_{\text{in}} = \int_{l} S_{\text{av}} dl \cdot H \qquad (2)$$
where *l* is the incident end of the fiber, *S*<sub>av</sub> is the Poynting vector and *H* refers to the height of the model. Significantly, the height is also applicable to the calculation of the following formula. In this model, *H* is calculated to be 1.939.

The profile of the power distribution is essential to evanescent-wave-based optical sensing. Taking the plane integral of the Poynting vector, Tong et al. calculated the fractional power inside the core ($T_1$) under single mode condition and took it as the evaluation criteria of the fiber sensing performance as shown in Figure 3 [31]. Here, $A_1$ and $B_1$ refer to the planes of core and cladding, respectively. However, the evanescent field power changes periodically in the direction of light under the multimode condition due to mode interference. The formula calculating the evanescent field ($T_o$) can be optimized by taking the volume integral of the average power. These two calculation methods are only applicable to adiabatic taper due to that the structural energy loss is negligible. To obtain arbitrary structural loss, the optimized evanescent proportion ($T_o$) and evanescent efficiency ($\tau$) are calculated as follows:

$$T_1 = \frac{\iint_{A_1} S_{av} dB}{\iint_{A_1+B_1} S_{av} dA} \qquad (3)$$

$$T_O = \frac{\iiint_{A_2} W_{av} dB}{\iiint_{A_2+B_2} W_{av} dA} \qquad (4)$$

$$\tau = P \cdot T_O = \left[ \sum (S_{21})^2 \right] \cdot \frac{\iiint_{A_2} W_{av} dB}{\iiint_{A_2+B_2} W_{av} dA} \qquad (5)$$

where $W_{av}$ is the average energy density, and $A_2$ and $B_2$ refer to the volume of liquid and core marked in Figure 3. More formula details can be obtained from the Supplementary Materials.

<figure><img src="image_7.png" /><figcaption>Figure 3. Calculation of the fractional power inside the fiber. $A_1$ and $B_1$ refer to the plane of core and cladding, respectively. $A_2$ and $B_2$ refer to the volume of liquid and core, respectively.</figcaption></figure>

### 3. Effect of Waist Radius

The effect of $R_w$ is first analyzed assuming a fixed $L_t$ with 0.4 mm. The mode numbers of the fiber can accommodate are directly related to the normalized frequency:

$$V = \kappa_0 R_w \sqrt{n_1^2 - n_2^2} \qquad (6)$$

where $\kappa_0$ is the vacuum wave number, $R_w$ is the radius of the fiber and $n_1$ and $n_2$ are the core and cladding indices of the fiber, respectively. The electric field distributions for radius with 1 µm, 2 µm, 3 µm and 5 µm at 7 µm wavelength are shown in Figure 4. Generally, energy from fundamental mode transfer to high-order modes occurs during the transition from core to cladding because of the non-adiabaticity. The mode capacity decreases as the radius decreases further. A noticeable leakage is observed in the transition region between cone and waist when $R_w$ with 1 µm can only accommodate the fundamental mode, as shown in Figure 4a. More details are shown in Figures S1 and S2. The calculation results of the light transmittance ($T$), $T_o$, and the $\tau$ are 42.91%, 7.85% and 3.87%, respectively. Remarkably, the $\tau$ is much smaller than $T_o$ due to the non-adiabatic taper. Figure 4b shows that the electric field distribution inside the waist changes periodically along the light
propagation direction field when $R_w$ is 2 µm. The position of the electromagnetic vibration peaks indicates that the $LP_{11}$-like mode appears. Figure 4c demonstrates that the electric field distributions become more complex when $R_w$ are 3 µm and 5 µm, respectively. It can be seen from Table 1 that $T$ increases significantly from 42.91% to 98.72% and $T_o$ decreases gradually from 7.85% to 1.67% with the increase in $R_w$. Meanwhile, $\tau$ increases first and then decreases, indicating that some interesting phenomena occurred when $R_w$ changed from 1 µm to 3 µm.

<figure><img src="image_8.png" /><figcaption>**Figure 4.** Schematic of the modal evolution in the transition region with (a) $R_w$ = 1 µm; (b) $R_w$ = 2 µm; (c) $R_w$ = 3 µm; (d) $R_w$ = 5 µm; (e) calculated modes proportion; (f) light transmittance, evanescent power ratio, evanescent efficiency in the waist with different $R_w$.</figcaption></figure>

<table><caption>**Table 1.** Calculated light transmittance $T$, evanescent power ratio $T_o$, evanescent efficiency $\tau$.</caption><thead><tr><th>R<sub>w</sub> (µm)</th><th>Mode Evolution</th><th>T (%)</th><th>T<sub>o</sub> (%)</th><th>τ (%)</th></tr></thead><tbody><tr><td>1</td><td>single mode → multimode → single mode</td><td>42.91</td><td>7.85</td><td>3.87</td></tr><tr><td>2</td><td>single mode → multimode → multimode</td><td>76.62</td><td>7.09</td><td>5.43</td></tr><tr><td>3</td><td>single mode → multimode → multimode</td><td>83.78</td><td>5.01</td><td>4.20</td></tr><tr><td>5</td><td>single mode → multimode → multimode</td><td>98.72</td><td>1.67</td><td>1.65</td></tr></tbody></table>

Figure 4e shows the power ratio of different modes of $R_w$ from 0.8 µm to 9 µm. With the increasing of $R_w$, the transmittance increases as a result of the increasing $V$, and $\tau$ is lower than $T_o$ on the whole due to the optical loss caused by the non-adiabatic taper. When $R_w$ is below 1.5 µm, only the fundamental mode can propagate in the waist. With $R_w$ gradually decreasing, the electric field confinement of the waist gradually weakens, resulting in a larger evanescent field. Some peaks from the line of $T_o$ and $\tau$ appeared when $R_w$ was above 1.5 µm, corresponding to the cut-off radius of the higher-order mode. As $R_w$ increases further, the proportion of modes in the waist remains the same, while the
confinement of the waist on the electromagnetic field strengthens further, and $T_o$ and $\tau$ start to decrease again. At the same time, as $T$ approaches 1, the difference between $T_o$ and $\tau$ gradually disappears (Figure 4f), which shows evanescent efficiency is suitable for both adiabatic and non-adiabatic taper.

The simulation shows that $T_o$ and $\tau$ do not decrease monotonically with the $R_w$ in the case of non-adiabatic taper. By increasing the taper length ($L_t$), the non-adiabatic taper will gradually convert to an adiabatic taper. Next, the influence of $L_t$ on light transmittance, evanescent field proportion and evanescent efficiency is analyzed.

### 4. Effect of the Length of Taper

Firstly, $R_w$ is assumed to be fixed at 1.45 µm, on which condition the waist is single-mode fiber. Parametric scanning is applied to investigate the effect of $L_t$ on waveguiding properties by increasing $L_t$ from 0.3 mm to 0.9 mm with a 0.01 mm step. Figure 5a shows the electromagnetic field distribution in the waist with different $L_t$, indicating that more energy is conserved when the cone is longer. Figure 5c demonstrates that more energy is coupled to the fundamental mode instead of the leaky mode with the increase in $L_t$. The calculated $T_o$ shows a downward trend due to the increasing fundamental mode. Meanwhile, $\tau$ shows an upward trend due to the significant increase in $T$.

<figure><img src="image_9.png" /><figcaption>**Figure 5.** Calculated electric field in the waist as (a) $R_w$ = 1.45 µm and (b) 2 µm; calculated modes proportion, light transmittance, evanescent power ratio, evanescent efficiency with different $L_t$ as $R_w$ fixed to (c) 1.45 µm and (d) 2 µm; calculated (e) modes proportion; (f) light transmittance, evanescent power ratio, evanescent efficiency in the waist as $R_w$ is fixed at 4 µm.</figcaption></figure>

Secondly, $R_w$ is assumed to be fixed at 2 µm, on which condition the $LP_{11}$-like mode can be transmitted in the waist. Figure 5b shows the electromagnetic field distribution in the waist. The energy of the $LP_{11}$-like mode center at the boundary of the waist benefits the interaction between light and molecular groups. The calculated power proportion of different modes as a function of $L_t$ is shown in Figure 5d. It is noticed that more energy couples to the fundamental mode as $L_t$ increases, while the proportion of the $LP_{11}$-like mode increases first and then decreases. The $T$, $T_o$ and $\tau$ in the waist as a function of $L_t$ are also presented in Figure 5d, respectively.

Finally, $R_w$ is assumed to be fixed at 4 µm, on which condition three kinds of modes can transmit in the waist, as Figure 5e,f show. The $T_o$ and $\tau$ are quite small because the light field is well restricted in the waist with a larger $R_w$, indicating a poor sensing performance.
**5. Optimization of Multimode Tapered Fiber**

In this section, different $R_w$s which are close to the cut-off radius of the LP$_{11}$-like mode are chosen for the $L_t$ optimization calculation. Figure 6a–c give the calculated power proportion of LP$_{01}$, LP$_{11}$ and light transmittance, respectively. With the increase in $L_t$, the proportion of LP$_{01}$ increases, and the proportion of LP$_{11}$ increases first and then decreases, which is consistent with the above results. Comparing Figure 6a,b,d,e, when $R_w$ increases from 1.52 µm to 1.64 µm, the fundamental mode proportion stays unchangeable, and more LP$_{11}$ transmits in the waist, while evanescent proportion decreases gradually. Notably, evanescent efficiency increases first and then decreases, which can be seen clearly in Figure 6f. Obviously, maximum evanescent efficiency exists at the place where the LP$_{11}$ mode proportion does not reach its extreme. For the fiber with an original diameter of 100 µm in this simulation, when $R_w$ and $L_t$ are about 1.54 µm and 0.42 mm, respectively, the fiber sensor obtains the highest evanescent efficiency of 15%. Meanwhile, we can control the $R_w$ in the range of 1.54~1.6 µm and the $L_t$ in the range of 0.35~0.7 mm by the taper experiment, which is shown in Figure 2 to obtain a fiber sensor with evanescent efficiency higher than 10%.

<figure><img src="image_10.png" /><figcaption>**Figure 6.** Calculated (a) fundamental mode proportion; (b) LP$_{11}$-like mode proportion; (c) light transmittance; (d) evanescent proportion; (e) evanescent efficiency with different $L_t$ as $R_w$ changes from 1.5 to 1.61 µm; (f) calculated evanescent efficiency with different $R_w$ as $L_t$ changes from 0.38 to 0.46 mm.</figcaption></figure>

By numerical simulation, we obtained the optimized geometry of tapered microfiber, which is very important for obtaining a high-precision and highly repeatable fiber sensor. It is noteworthy that this method can also be applied to the size design of U-shaped fiber, D-shaped fiber and other fibers with special shapes.

**6. Analysis of Micro Deformation in Waist**

The optimal size of the tapered As$_2$Se$_3$ fiber was achieved in Section 5, which could be realized by the fiber-puller system mentioned in Section 2. However, for the demand of liquid pool detection, in which condition a long waist for sensing is necessary, it is too hard to realize the geometry mentioned above. Thus, this geometry is only expected to be applied to single droplet detection, as the waist length is small. Therefore, we make a further calculation for the case that there is a micro cone in the waist.

The waist region of the fiber with the micro taper is a part with changing radius, and the change is quite slow, so it can be approximated as a straight line with a small slope.
Because the decreasing of fiber diameter will cause the change in mode, the light field in the micro cone is actually a process of coupling from high-order mode to low-order mode. Figure 7a gives the mathematical model for this section. The $R_L$ increases from 1.56 µm to 3.5 µm to show the micro deformation of the fiber waist. Figure 7b,c give the mode proportion at $R_w$ and the calculated $T$, $T_o$ and $\tau$ in the waist, respectively. We find that with the increase in $R_L$, the distribution of mode at $R_w$ changes little, while the calculated $T_o$ and $\tau$ decrease gradually. In other words, when the deformation is very small in the waist, the mode distribution is mainly determined by $R_w$. This is because, on the one hand, $R_w$ is fixed, while the mode field distribution in the optical fiber is mainly related to its radius. On the other hand, the taper in the waist area is very small, resulting in the mode conversion efficiency close to 1. Therefore, the mode proportion at $R_w$ is almost unchanged at a different $R_L$. Meanwhile, as the radius changes from $R_L$ to $R_w$, the equivalent radius of the waist area is greater than $R_w$. The larger the radius is, the smaller the $\tau$ is (as shown in Figure 4). Although the evanescent efficiency decreases with $R_L$, it can still be maintained above 6%.

<figure><img src="image_11.png" /><figcaption>Figure 7. (a) Mathematical model; (b) calculated modes proportion; and (c) light transmittance, evanescent proportion and evanescent efficiency with different $R_L$.</figcaption></figure>

In some cases, low evanescent efficiency does not mean poor sensing performance. As shown in Figure 1, raising the length of the sensing region can also increase the light absorbed by the molecular groups but needs more liquids. Although the deformation decreases the evanescent efficiency, it also reduces the requirement of precision in the fiber taper, which enables the manufacturing of fiber with a longer waist. Deformation is inevitable during the taper processing, but if we control it properly, the sensitivity can be improved.

As shown in Figure 2, different kinds of heaters will provide different temperature fields which could flexibly control the taper geometry. In addition, according to the literature [11], optimizing the taper drawing method can realize the long waist area. In practical terms, the fragility of the tapered fiber can be overcome by encapsulation with the UV curing adhesive in a liquid cell. However, we assume the loss the tapered fiber can be neglected in this simulation work. In practical tapered fiber, as the mode of light cannot be confined in the fiber when the radius is too small, the loss increased dramatically with the decrease in the waist radius, which limits the application of the tapered fiber with very small radius.
## 7. Conclusions

We systematically calculated the mode distribution in chalcogenide tapered fiber with flexible size design and defined the evanescent efficiency characterizing its sensing performance. Generally, $R_w$ determines the mode capacity and the ability to limit light, and $L_t$ determines the mode's proportion and light transmittance. By optimizing $R_w$ and $L_t$, evanescent efficiency over 10% can be obtained for droplet liquid detection. For the liquid detection, although a high-precision taper cannot be achieved due to the inevitable micro deformation, we can lengthen the waist to obtain a better sensing performance. All in all, we give the general geometry design strategy of chalcogenide micro-nano fiber for sensing applications. While the tapered fiber has high sensitivity in liquid detection, there are still many challenges we have to face. The loss reduction and device encapsulation work can further promote the application of the chalcogenide fiber sensor. The experimental work on the taper fiber sensor is in progress.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ma15113834/s1, Figure S1: Transmission of light in a non-adiabatic taper; Figure S2: The variation in electric field vibration from fundamental mode to higher-order mode and to fundamental mode in tapered fiber; 40–470 represent the different positions from incident end to exit end.

**Author Contributions:** Conceptualization, X.Z. (Xudong Zhao) and Q.L.; Methodology, X.Z. (Xudong Zhao); and Z.L.; Software, L.Z. and G.T.; Formal Analysis, X.Z. (Xianghua Zhang) and X.Z. (Xiujian Zhao); Investigation, Y.X., G.T. and X.Z. (Xiujian Zhao); Data Curation, L.Z., Y.X. and X.Z. (Xianghua Zhang); Writing—Original Draft Preparation, X.Z. (Xudong Zhao), Z.L. and N.Y.; Writing—Review and Editing, Q.L., N.Y. and Y.X. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work is financially supported by the National Natural Science Foundation of China (61975156), Fundamental Research Funds for the Central Universities (203134001, 215201006), Hubei Natural Science Foundation (2020CFB641).

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The raw data required to reproduce these results cannot be shared at this time as the data also form part of an ongoing study.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Wu, Z.; Xu, Y.; Qi, D.; Nie, Q.; Zhang, X. Progress in preparation and applications of Te-As-Se chalcogenide glasses and fibers. *Infrared Phys. Technol.* **2019**, *102*, 102981–102991. [CrossRef]

2. Bensaid, S.; Kachenoura, A.; Costet, N.; Bensalah, K.; Tariel, H.; Senhadji, L. Noninvasive detection of bladder cancer using mid-infrared spectra classification. *Expert Syst. Appl.* **2017**, *89*, 333–342. [CrossRef]

3. Starecki, F.; Braud, A.; Doualan, J.-L.; Ari, J.; Boussard-Plédel, C.; Michel, K.; Nazabal, V.; Camy, P. All-optical carbon dioxide remote sensing using rare earth doped chalcogenide fibers. *Opt. Laser. Eng.* **2019**, *122*, 328–334. [CrossRef]

4. Yang, C.; Wang, X.; Su, J.; Dai, S. Spectroscopy analysis of mixed organic liquid detection with Ge₂₀Se₆₀Te₂₀ glass-tapered fiber. *J. Non Cryst. Solids* **2018**, *500*, 377–381. [CrossRef]

5. Wang, X.; Su, J.; Wang, Y.; Yang, C.; Dai, S.; Zhang, P. High-sensitivity sensing in bare Ge-Sb-Se chalcogenide tapered fiber with optimal structure parameters. *J. Non Cryst. Solids* **2021**, *559*, 120686–120691. [CrossRef]

6. Su, J.; Dai, S.; Gan, N. Optimized Ge-As-Se-Te chalcogenide glass fiber sensor with polydopamine-coated tapered zone for the highly sensitive detection of p-xylene in waters. *Opt. Express* **2020**, *28*, 184–193. [CrossRef] [PubMed]

7. Wang, M.; Yang, F.; Dai, S.; Cao, Z.; Su, J.; Ding, S.; Zhang, P. Effect of the Geometries of Ge-Sb-Se Chalcogenide Glass Tapered Fiber on the Sensitivity of Evanescent Wave Sensors. *J. Lightwave Technol.* **2021**, *39*, 4828–4836. [CrossRef]

8. Wang, L.; Ma, W.; Zhang, P.; Zhu, L.; Yang, D.; Wang, X.; Dai, S. Mid-Infrared Gas Detection Using a Chalcogenide Suspended-Core Fiber. *J. Lightwave Technol.* **2019**, *37*, 5193–5198. [CrossRef]

9. Yang, D.; Zhang, P.; Zeng, J.; Wang, X.; Song, B.; Dai, S. SRI-Immune Highly Sensitive Temperature Sensor of Long-Period Fiber Gratings in Ge-Sb-Se Chalcogenide Fibers. *J. Lightwave Technol.* **2017**, *35*, 3974–3979. [CrossRef]
10. Xie, Y.; Cai, D.; Pan, J.; Zhou, N.; Gao, Y.; Jin, Y.; Jiang, X.; Qiu, J.; Wang, P.; Guo, X.; et al. Batch Fabrication of High-Quality Infrared Chalcogenide Microsphere Resonators. *Small* **2021**, *17*, 2100140–2100146. [CrossRef]

11. Xie, Y.; Cai, D.; Wu, H.; Pan, J.; Zhou, N.; Xin, C.; Yu, S.; Wang, P.; Jiang, X.; Qiu, J.; et al. Mid-infrared chalcogenide microfiber knot resonators. *Photon. Res.* **2020**, *8*, 616–621. [CrossRef]

12. Jiang, X.; Jha, A. Engineering of a Ge–Te–Se glass fibre evanescent wave spectroscopic (FEWS) mid-IR chemical sensor for the analysis of food and pharmaceutical products. *Sens. Actuators B Chem.* **2015**, *206*, 159–169. [CrossRef]

13. Romanova, E.A.; Korsakova, S.V.; Rozhnev, A.G.; Velmuzhov, A.P.; Kotereva, T.V.; Sukhanov, M.V.; Shiryaev, V.S. Chalcogenide fiber loop probe for the mid-IR spectroscopy of oil products. *Opt. Express* **2020**, *28*, 5267–5272. [CrossRef] [PubMed]

14. Zhang, K.; Alamgir, I.; Rochette, M. Midinfrared Compatible Tunable Bandpass Filter Based on Multimode Interference in Chalcogenide Fiber. *J. Lightwave Technol.* **2020**, *38*, 857–863. [CrossRef]

15. Markos, C.; Bang, O. Nonlinear Label-Free Biosensing With High Sensitivity Using As₂S₃ Chalcogenide Tapered Fiber. *J. Lightwave Technol.* **2015**, *33*, 2892–2898. [CrossRef]

16. Bednarska, K.; Sobotka, P.; Woliński, T.R.; Zakręcka, O.; Pomianek, W.; Nocoń, A.; Lesiak, P. Hybrid Fiber Optic Sensor Systems in Structural Health Monitoring in Aircraft Structures. *Materials* **2020**, *13*, 2249. [CrossRef] [PubMed]

17. Marro Bellot, C.; Sangermano, M.; Olivero, M.; Salvo, M. Optical Fiber Sensors for the Detection of Hydrochloric Acid and Sea Water in Epoxy and Glass Fiber-Reinforced Polymer Composites. *Materials* **2019**, *12*, 379. [CrossRef] [PubMed]

18. Ermatov, T.; Skibina, J.S.; Tuchin, V.V.; Gorin, D.A. Functionalized Microstructured Optical Fibers: Materials, Methods, Applications. *Materials* **2020**, *13*, 921. [CrossRef]

19. Yu, H.; Gao, W.; Jiang, X.; Guo, H.; Jiang, S.; Zheng, Y. On-line Writing of Fiber Bragg Grating Array on a Two-mode Optical Fiber for Sensing Applications. *Materials* **2019**, *12*, 1263. [CrossRef] [PubMed]

20. Sheng, J.; Peng, Q.; Sun, W.; Wang, Y.; Yang, K.; Feng, Z.; Wang, J.; Wu, G.; Wang, J.; Jiao, K.; et al. Influence of extrusion on the properties of chalcogenide glasses and fibers. *Opt. Commun.* **2022**, *513*, 128091. [CrossRef]

21. Shiryaev, V.S.; Sukhanov, M.V.; Velmuzhov, A.P.; Karaksina, E.V.; Kotereva, T.V.; Snopatin, G.E.; Denker, B.I.; Galagan, B.I.; Sverchkov, S.E.; Koltashev, V.V.; et al. Core-clad terbium doped chalcogenide glass fiber with laser action at 5.38 µm. *J. Non-Cryst. Solids* **2021**, *567*, 120939. [CrossRef]

22. Esfahani Monfared, Y. Overview of Recent Advances in the Design of Plasmonic Fiber-Optic Biosensors. *Biosensors* **2020**, *10*, 77. [CrossRef] [PubMed]

23. Cennamo, N.; Arcadio, F.; Zeni, L.; Catalano, E.; Del Prete, D.; Buonanno, G.; Minardo, A. The Role of Tapered Light-Diffusing Fibers in Plasmonic Sensor Configurations. *Sensors* **2021**, *21*, 6333. [CrossRef]

24. Carvalho, W.O.F.; Mejía-Salazar, J.R. Surface Plasmon Resonances in Sierpinski-Like Photonic Crystal Fibers: Polarization Filters and Sensing Applications. *Molecules* **2020**, *25*, 4654. [CrossRef] [PubMed]

25. Tong, L. Micro/Nanofibre Optical Sensors: Challenges and Prospects. *Sensors* **2018**, *18*, 903. [CrossRef] [PubMed]

26. Zhong, N.; Xin, X.; Liu, H.; Yu, X.; Chang, H.; Tang, B.; Zhong, D.; Zhao, M.; Zhang, H.; Zhao, J. Plastic optical fiber sensor for temperature-independent high-sensitivity detection of humidity. *Appl. Opt.* **2020**, *59*, 5708–5713. [CrossRef]

27. Xin, X.; Zhong, N.; Liao, Q.; Cen, Y.; Wu, R.; Wang, Z. High-sensitivity four-layer polymer fiber-optic evanescent wave sensor. *Biosens. Bioelectron.* **2017**, *91*, 623–628. [CrossRef]

28. Zhang, L.; Tang, Y.; Tong, L. Micro-/Nanofiber Optics: Merging Photonics and Material Science on Nanoscale for Advanced Sensing Technology. *iScience* **2020**, *23*, 100810–100832. [CrossRef]

29. Ravets, S.; Hoffman, J.E.; Kordell, P.R.; Wong-Campos, J.D.; Rolston, S.L.; Orozco, L.A. Intermodal energy transfer in a tapered optical fiber: Optimizing transmission. *J. Opt. Soc. Am. A* **2013**, *30*, 2361–2371. [CrossRef]

30. Brambilla, G. Optical fibre nanowires and microwires: A review. *J. Opt.* **2010**, *12*, 043001–043019. [CrossRef]

31. Tong, L.; Lou, J.; Mazur, E. Single-mode guiding properties of subwavelength-diameter silica and silicon wire waveguides. *Opt. Express* **2004**, *12*, 1025–1035. [CrossRef] [PubMed]