New Trends in the Simulation of Nanosplasmonic Optical D-Type Fiber Sensors

This article presents a review of the numerical techniques employed in simulating plasmonic optical sensors based on metal-dielectric nanostructures, including examples, ranging from conventional D-type fiber sensors, to those based on photonic crystal D-type fibers and incorporating metamaterials, nanowires, among other new materials and components, results and applications. We start from the fundamental physical processes, such as optical and plasmonic mode coupling, and discuss the implementation of the numerical model, optical response customization and their impact in sensor performance. Finally, we examine future perspectives.


Introduction
Nowadays, fiber optical sensors integrate information and monitoring systems, as well as different types of devices, with applications from infrastructure and industrial monitoring to environmental control, largely due to their high sensitivity and immunity to electromagnetic and radiofrequency interference [1]. A considerable part of these sensors are refractive index sensors that have been functionalized to respond indirectly to a multitude of parameters, from pH, temperature, humidity, to carbon dioxide (CO 2 ) concentration and many other physical, chemical and biological parameters. This technology has benefited from the exploration of surface plasmon resonance (SPR) that occurs at the interface between metal and dielectric films deposited on the surface of optical fibers, reaching sensitivities in excess of 10 −7 RIU [2].
Still, the recent development of sophisticated fabrication techniques capable of producing optical fibers, incorporating structures that combine dielectric, metallic and semi-conducting materials, with sizes below the micron, reaching in some cases even a few tens of nanometers, has paved the way to consider the development of optical fiber sensors that explore nanophotonic phenomena as sensing principles. In optical structures with these sizes, light can no longer be described by ordinary optical rays or even planewaves, as the strongly localized boundary conditions, imposed at the interface between different materials, force light into optical modes that are strongly dependent of the geometry and material optical properties. This dependency can be explored to customize the optical response and optimize the performance of the sensors, however, a development approach based on their actual construction and testing is very expensive and time consuming. Instead, using numerical models, to simulate the sensor response, constitutes an important design tool to obtain a first and expedite approximation of the performance of a sensor design, that can provide insight on the optical phenomena involved, as well as, a significant reduction of the development and optimization costs.
Typically, the propagation of light in optical structures much larger than the wavelength can be described using ray-tracing models [3,4], whereas the performance of interferometric sensors or those based on overlapping of thin metal-dielectric films usually is calculated using simulations based on the propagation of plane waves, such as the expansion and propagation method (MEP) [5], the method for multilayer structure transfer matrix modeling [6,7], the optical fiber multilayer cylindrical structure [8], the rigorous coupled wave analysis (RCWA) [9] and sometimes wave-mode coupling [10,11]. However, when considering optical structures with sizes of the order or even below the wavelength, and with complex geometries, these methods fail short. Optical devices and nanostructured materials such as photonic crystal fibers (PCF) [12], metamaterials [13], metallic nanowires [14][15][16], etc. require a distinct approach that considers the strongly localized optical modes supported and their mutual influence, that can adequately calculate the spatial distribution of optical power and the phase effects [17], using typically meshing methods. These methods involve the sampling of the electromagnetic field at a specific number of nodes and establishing a relation between them.
One such approach is the finite-difference time domain (FDTD) method [17]. It consists of solving partial differential equations (PDE) in the time domain, being easier to understand and it is a more commonly used method to solve basic problems. Although FDTD provides a good description of the propagation of light in optical structures, it requires a thin regular mesh and consumes large amount of computer memory [18].
Alternatively, the finite element method (FEM) is preferable for waveguiding structures [19,20], such as those characteristic of optical fiber sensors, as it allows to determine the optical modes (sometimes referred as spatial or normal modes) with a relatively low computer memory requirements by considering an irregular mesh that accommodates the specific geometry of the device and their structural components. The FEM is usually used when it is necessary to vary the resolution by which the field is described over different regions of the simulation box. The drawback in using the FEM is that the process to reduce the PDEs for the electromagnetic field into an eigenvalue problem is typically more complex than using FDTD. However, combined with the facility of automated meshing provided by simulation platforms, such as COMSOL or ANSYS, the FEM is becoming the preferred simulation method.
This article reviews some of the theoretical and numerical foundations and the numerical techniques employed in simulating plasmonic optical fiber sensors based on metal-dielectric nanostructures. We describe from conventional multilayer sensors based on D-type fibers, to those based on photonic crystal fibers and incorporating metamaterials, nanowires, among other new materials and components.
The article is organized as follows: after the introduction and in the second section we review some of the fundamental physical principles that underline the operation of this type of sensors, such as the equations for optical and plasmonic modes and their coupling. In Section 3 we present some of the numerical models used to describe the sensors and how to compute their optical response and sensor performance. In particular, we focus on FEM. Section 4 provides some examples of application of these methods. Finally, Section 5 discusses future perspectives and delivers some concluding remarks.

The Fundamentals of Nanoplasmonics
The propagation of light in a nonmagnetic optical medium is governed by the Maxwell equations: where E(r, t) is the electric, D(r, t) is the electric displacement, B(r, t) is the magnetic induction and µ 0 is the magnetic vacuum permeability. Combining these two equations, using continuous harmonic waves of the form E(r, t) = E(r, ω) exp(iωt) and the constitutive relation D(r, ω) = 0 n 2 (r, ω)E(r, ω), yields the wave equation: with r = (x, y, z) and where µ 0 0 ω 2 = k 2 0 , ω is the angular frequency, 0 is the vacuum permittivity and n(r, ω) is the refractive index, which is a complex quantity and can be decomposed into the real and imaginary parts, say n(r, ω) = n r (r, ω) + in i (r, ω). For simplicity, the optical axis is assumed to be the z direction.
This paper considers optical fibers capable of guiding optical and plasmonic excitations characterized by a constant profile of refractive index along the optical axis, i.e., n(r, ω) = n(r ⊥ , ω), and for which it is possible to apply separation of variable with E(r, ω) = E(r ⊥ , ω) exp(−iβz) to obtain: If the optical structures are domains with homogeneous material with well-defined and uniform refractive index n(ω) then, Equation (3) becomes a simpler eigen-problem for each domain m: with a common eigenvalue β, a fixed value k 0 and subjected to the adequate boundary conditions: wheren is the unit normal pointing from domain with subscript (+) into that with subscript (−). The boundary conditions state imposes that every component of the magnetic field and the tangential components of the electric field is continuous across the interface while, the normal component of the electric field times has a discontinuity proportional to the induced surface charge density Σ at the boundary. The combination between the eigenvalue problems in each domain and the boundary conditions usually forces the spectrum of eigenvalues to be discrete. It is common to define the effective refractive index for each mode n e f f according to β = n e f f k 0 , which can also have real and imaginary parts. While the real part of the effective refractive index provides the phase velocity of the guided mode along the fiber, the imaginary part characterizes the losses, typically the physical quantity used to inspect the properties of the environment (such as the refractive index of the analyte). The loss curve can be obtained from the imaginary part of the parameter β or alternatively of the effective refractive index of the mode n e f f (ω) as α(ω) = 2Jm(β) = 2Jm n e f f k 0 .
Multiplying Equation (3) by E * (r ⊥ , ω) and integrating along the transverse direction, one obtains: which implies that the effective refractive index of each mode is in some sense the average of the refractive index over the domains of the optical structure. The contribution of each domain to n e f f is proportional to the relative optical power of the mode in that domain, and therefore for some modes, it is possible to consider only the refractive index structure where the mode is localized. This simplification is convenient to obtain analytical (and sometimes numerical) approximations of the guided modes. For example, the core of an optical fiber is usually a cylinder of radius R of a transparent dielectric with larger refractive index than the surroundings. If their refractive indexes can be assumed to have only a real part, then the guided modes are given in terms of Bessel and von Newmann functions of first kind: where the modal parameters are defined as: It results that U is real, while W is imaginary, and the mode is strongly localized in the region r ⊥ < R with an evanescent part in the cladding. Replacing the core with a metallic wire supporting plasmonic surface waves and if the refractive index in the metal is approximately imaginary, then the mode is evanescent in both metal and cladding: where I is the modified Bessel function.
Complex optical structures can usually support several localized modes that partially overlap, resulting in the coupling between them. For example, consider the localized modes calculated for a core near a metallic nanowire (or any other two structures capable of supporting localized modes): where n 1 (r ⊥ , ω) describes the refractive index structure of the core neglecting the presence of the wires and n 2 (r ⊥ , ω) describes the refractive index structure of the wire neglecting the presence of the core. Notice that the total refractive index distribution combining both the core and the wire satisfies: where ∆n 1 (r ⊥ , ω) and ∆n 2 (r ⊥ , ω) are corrections to the refractive index distributions of n 1 (r ⊥ , ω) and n 2 (r ⊥ , ω) due to the presence of the core and the wire, respectively. Then, replacing in Equation (3) a tentative solution as a linear composition of the modes supported by each of the two structures, say E(r ⊥ , ω) = c 1 E 1 (r ⊥ , ω) + c 2 E 2 (r ⊥ , ω), one obtains: For simplicity, we have assumed that the two localized modes are orthogonal and normalized, in other words that they satisfy E * i (r ⊥ , ω)E j (r ⊥ , ω) d 2 r ⊥ = δ ij , which is not always exactly the case, but constitutes a good approximation for the majority of situations. Then, defining the following integrals: multiplying the previous equation either by E * 1 (r ⊥ , ω) or by E * 2 (r ⊥ , ω) and integrating over the transverse direction results in the following eigen-problem: with eigenvalues are β 2 ± = β 2 ± ∆β 4 + k 4 0 g 211 g 122 and eigenvectors: The eigenvectors constitute the modes of the structure that combines both a core and a nanowire, and result from the hybridization of the localized modes supported by both substructures. Obviously, this analysis is quite general and can be applied to any type of substructures beside a core and a nanowire, but is an excellent model to explore the effects of mode hybridization. For example, if the overlap integrals g ijk are small, then β + ≈ β 2 and β − ≈ β 2 and the two modes do not couple with each other, unless β 1 ≈ β 2 ≈ β. In other words, when the two modes are degenerated in β. In this latter case, corresponding to a crossing point in the dispersion curve of the two localized modes, the coupling and exchange energy between the two localized modes is extremely strong even if the overlap is small. Then, the resonant coupling between the two modes introduces a gap of 2 k 2 0 (g 111 + g 222 ) in the values of β since: In some sense, this resonant coupling resembles the hybridization of atomic orbitals that occurs during the formation of a chemical bond, resulting in the formation of a bonding molecular orbital with a decrease of energy and an anti-bonding molecular orbital with an increase of energy relative to the energy of the original orbitals. The processes of mode hybridization are extremely important for optical sensing as they generate coupling resonances at very specific wavelengths, when the dispersion curves of different localized modes cross and become degenerated. If one of the modes is plasmonic in nature and characterized by strong losses, the coupling at that specific wavelength allows to transfer power from the optical into the plasmonic mode, which is then dissipated, resulting in a peak in the loss curve of the optical mode. The key for optical sensing using this type of processes is to promote the dependency of the properties (usually the dispersion curve) of one of the modes on a specific parameter of an analyte (typically, the refractive index). As that property of the analyte changes, it alters the dispersion curve of the mode and shifts the wavelength of the crossing point and thus the wavelength of the peak in the loss curve.
This simple example of a two-mode hybridization clarifies at a fundamental level the optical processes that underline most sensors based on localized optical and plasmonic modes in micro and nanostructured optical fibers. However, these principles are not limited to just two modes hybridization but can be extended to multiple modes hybridization that, in principle, result in much more critical resonance conditions and increased sensor sensitivities. The following section will elaborate on the implementation of numerical methods targeted at capturing and evaluating these processes, as well as engineering them to improve sensor performance.

Simulation Models
As discussed in the previous section, the key for the analysis of the performance of optical fiber sensors with micro-and nanoscale metal-dielectric structures, capable of supporting optical and plasmonic modes, is the analysis of the eigenproblem associated with Equation (3). If the structure of the sensor is composed only by cylindrical metal-dielectric components, such as optical waveguides (cores) and nanowires, for which there are good approximate analytical solutions, it is possible to expedite the calculation of the mode hybridization using simple numerical models [21,22]. Other methods, such as the expansion and propagation method [5] and the rigorous coupled wave analysis [9,23,24] can provide improved results to the mode hybridization. On the other hand, when considering devices composed of multiple layers of different materials, there are specific methods capable of addressing these problems, which include the method for multilayer structure transfer matrix modeling [6,[25][26][27][28][29] and the optical fiber multilayer cylindrical structure [8]. Yet these methods have some difficulties for devices with nanostructured irregularities, or more complex and intricated structures. In these cases, one can resort to meshing based methods, such finite difference time domain (FDTD) [30][31][32], spectral methods and finite element methods (FEM) [33][34][35][36][37][38][39][40].
As referred earlier, FDTD consists of solving partial differential equations in the time domain, being easier to understand and it is a more commonly used method to solve basic problems. FDTD samples the electromagnetic field in the simulation domain at points in a regular mesh and approximates the spatial derivatives in Equation (3) using finite differences. For example, the gradient in two dimensions of a function φ can be approximated using symmetric first order finite differences as: where h is the spacing between the mesh points and the subscripts i and j refer to the order of the mesh point, along each dimension of the grid. In very complex geometries, it is required substructures with different sizes, being necessary a dense mesh of sampling points, which requires large computing memory and often results in long simulation times. An alternative are the spectral methods, and particularly the Fourier spectral methods, which consider the decomposition of the solution as a sum over a set of base functions, say: which after substitution in Equation (3) also results in a set of algebraic eigenproblems, that can be solved numerically. In some way, these methods extend on the method of mode hybridization described in the previous section and that corresponds to a simple case with only two base functions. In general, the adopted basis is composed of plane waves to take advantage of highly efficient Fast Fourier Transforms implementations available, which reduce memory consumption and computing time. However, such approach is better suited for problems with cyclic or periodic boundary conditions, such as photonic crystals and metamaterials. Finally, FEM considers the sampling of the field at the nodes of an irregular mesh based on the construction of sub-elements within the structure to be studied. For each structure, the calculations of the field equations are represented discreetly in a system of algebraic equations and are solved by an eigenvalue method. FEM is usually used when it is necessary to analyze the field with different resolutions in different regions of the simulation domain, and thus provides a compromise between accuracy, resolution and consumption of computer memory resources. The drawback is that the quality of the solutions obtained by the FEM methods is strongly dependent of the characteristics of the mesh, and finding the best mesh for each specific problem is not a trivial problem, even using the current automatic meshing provided by commercial simulation tools such as COMSOL Multiphysics (1) Reduction of the numerical error due to rounding and respective propagation during the calculation, which is obtained for meshes with fewer elements and, consequently reduces the number of calculations in the simulations. (2) Reduction of methodological error associated with the finite element method, which is achieved using more samples of the field, thus requiring grids with more elements. However, there are drawbacks to choosing meshes with many elements as this increases the complexity of simulations, requiring more computational resources and more time to perform calculations.
The balance between these two sources of error (numerical and methodological) is achieved by selecting grids where there is a denser sampling of the field (smaller finite elements) in areas where the field has larger variations and lowering the sampling in areas, where the field is quasi constant or zero. When comparing the different methods available to compute the modes of an optical structure, clearly the FEM method presents the best advantages and will be used in the next section to analyze the performance of several examples of optical fiber sensors.

Examples of New Nanoplasmonic Optical Fiber Sensors
The previous sections introduced the physical fundaments and compared the main numerical methods behind the modeling of nanostructured optical fiber sensors. In particular, we have shown that the key objective is to determine the optical modes supported by each structure and characterize them in terms of their optical properties, namely the effective refractive index, which allows to determine the losses in each section of the sensors, among other things. Here we have focused on sensors with a uniform refractive index distribution along the entire length of the fiber and whose operation strategy consists in determining the changes in the loss spectrum (i.e., as a function of wavelength) produced by the analyte or changes in the environmental refractive index, however the same principles and methods can be applied to other designs and operation approaches.
Moreover, we have established the importance of the localized modes and how they can couple in determining the features of the loss spectrum of the sensors. In particular, we have shown that the strong hybridization of localized modes when their dispersion curves cross, can introduce important features (such as peaks and deeps) in the loss spectrum that are strongly localized in wavelength.
As a result, one of the key elements of sensor design is to engineer the localized modes not only to promote this hybridization but also, and more importantly, to introduce a strong dependency of some of the localized modes in the parameters that the sensor aims to monitor in order to affect the resonant coupling. The critical nature of the resonant coupling amplifies the impact of the change of external parameters on each localized mode and expresses it, for example, in the loss spectrum of the sensor. The importance of numerical methods in the engendering these sensors results from the facility that it allows in the testing of different designs and the optimization of their performance, as illustrated by the examples depicted in the remaining of this section.

D-Type Fiber Configurations (Conventional and PCF) with Metallic Films
Here we compare the impact of engineering a D-type photonic crystal structure in the sensor performance with a more conventional optical structure based on a D-type fiber. In both configurations the flat surface, has been covered with a thin metal layer ( Figure 1). The key idea is to explore the role of the photonic crystal nanostructure in customizing the optical properties of the plane wave plasmon modes to ultimately promote enhancement of the sensor performance.

D-Type Fiber Configurations (Conventional and PCF) with Metallic Films
Here we compare the impact of engineering a D-type photonic crystal structure in the sensor performance with a more conventional optical structure based on a D-type fiber. In both configurations the flat surface, has been covered with a thin metal layer (Figure 1). The key idea is to explore the role of the photonic crystal nanostructure in customizing the optical properties of the plane wave plasmon modes to ultimately promote enhancement of the sensor performance. Using the Sellmeier equation [25] allows one to calculate the refractive index. The fibers have a D-type profile with a metallic layer deposited in the flat surface, whose the refractive index is obtained from the Drude model. It is assumed the analytic medium to be studied, with a refractive index of next, fills the space outside the fiber. The distance between the gold layer and the center of the fiber is d -residual cladding, the gold layer thickness is dm, the holes' diameter is dhole, and the holes' separation is Λ -pitch [19].
The conventional D-type fiber is based on a singlemode configuration (rco = 4 μm), being the edge of the fiber placed at d = 4.05 μm from the center of core, the index refraction of the core is nco = 1.476 and for the cladding is ncl = 1.452, for a wavelength of 900 nm. In the case of the PCF D-type fiber, the configuration is dhole = 1.61 μm, Λ = 2.3 μm and d = 3.2 μm. For both fibers, the metal thickness is dm = 45 nm, and the sensor length is L = 1 mm. Figure 2 represents the loss as function of wavelength for both configurations (conventional and PCF D-type fibers) for external refractive indices of next= {1.36, 1.38, 1.39}.  Figure 1b shows the configuration of PCF D-type optical fiber plasmonic sensor, composed by a core with refractive index (n co ) and surrounded by an array of dielectric structures (the holes in the microstructured fiber) with refractive index equal to the air [41].
Using the Sellmeier equation [25] allows one to calculate the refractive index. The fibers have a D-type profile with a metallic layer deposited in the flat surface, whose the refractive index is obtained from the Drude model. It is assumed the analytic medium to be studied, with a refractive index of n ext , fills the space outside the fiber. The distance between the gold layer and the center of the fiber is d -residual cladding, the gold layer thickness is d m , the holes' diameter is d hole , and the holes' separation is Λ -pitch [19].
The conventional D-type fiber is based on a singlemode configuration (r co = 4 µm), being the edge of the fiber placed at d = 4.05 µm from the center of core, the index refraction of the core is n co = 1.476 and for the cladding is n cl = 1.452, for a wavelength of 900 nm. In the case of the PCF D-type fiber, the configuration is d hole = 1.61 µm, Λ = 2.3 µm and d = 3.2 µm. For both fibers, the metal thickness is d m 45 nm, and the sensor length is L = 1 mm. Figure 2 represents the loss as function of wavelength for both configurations (conventional and PCF D-type fibers) for external refractive indices of n ext = {1.36, 1.38, 1.39}. To determine the optical performance of the fiber sensors, the attenuation (loss) and the distribution of the Poynting vector in each fiber are observed. The configurations under comparison aim to increase the light interaction with the external medium (permitting more light to escape in this direction). From Figure 2 losses are higher for the PCF D-type fiber configuration, indicating a stronger plasmonic effect due to the higher energy transfer. Figure 3 illustrates the normalized Poynting vector amplitude 1D across the fiber core, for an external refractive index next = 1.39 and for different wavelengths, considering (a) the conventional and (b) the PCF D-type fiber structures. The results indicate the PCF D-type fiber configuration allows a higher interaction of the light with the external medium, around a two times fold increase in the intensity of the Poynting vector. To determine the optical performance of the fiber sensors, the attenuation (loss) and the distribution of the Poynting vector in each fiber are observed. The configurations under comparison aim to increase the light interaction with the external medium (permitting more light to escape in this direction). From Figure 2 losses are higher for the PCF D-type fiber configuration, indicating a stronger plasmonic effect due to the higher energy transfer.  To determine the optical performance of the fiber sensors, the attenuation (loss) and the distribution of the Poynting vector in each fiber are observed. The configurations under comparison aim to increase the light interaction with the external medium (permitting more light to escape in this direction). From Figure 2 losses are higher for the PCF D-type fiber configuration, indicating a stronger plasmonic effect due to the higher energy transfer. Figure 3 illustrates the normalized Poynting vector amplitude 1D across the fiber core, for an external refractive index next = 1.39 and for different wavelengths, considering (a) the conventional and (b) the PCF D-type fiber structures. The results indicate the PCF D-type fiber configuration allows a higher interaction of the light with the external medium, around a two times fold increase in the intensity of the Poynting vector.  Table 1 compares both sensing configurations. For the conventional D-type fiber configuration, one metallic layer thickness is considered, dm = 45 nm, while for the PCF D-type fiber configuration two metallic layer thicknesses are considered, dm = 45 nm and dm = 65 nm, respectively. The resolution (R) and the sensitivity (S) were obtained considering an experimental spectral variation detection of 0.1 nm and using the data of Figure 2 Table 1 compares both sensing configurations. For the conventional D-type fiber configuration, one metallic layer thickness is considered, d m = 45 nm, while for the PCF D-type fiber configuration two metallic layer thicknesses are considered, d m = 45 nm and d m = 65 nm, respectively. The resolution (R) and the sensitivity (S) were obtained considering an experimental spectral variation detection of 0.1 nm and using the data of Figure 2 and the following equations: where λ min is the minimum wavelength value between two spectral lines that can be detected experimentally and λ peak is the wavelength shift of the resonance peak obtained from the simulations for different external refractive indices values (n ext ). The PCF D-type fiber sensor, with d m = 45 nm, presents a better performance, enhancing the sensitivity from 2.8 × 10 3 to 10.2 × 10 3 nm/RIU and the resolution from 3.6 × 10 −5 to 9.8 × 10 −6 RIU. The relevance of engineering the structure is confirmed by the results. The optimization of the available light energy near the metallic film, enhances the plasmonic effect, which in turn improves the sensitivity and resolution of the SPR sensor.

D-Type Fiber Configuration with Metallic Nanowire
Leaving the more conventional plasmonic sensing configuration that integrates a metallic film, addressed in the previous section, which supports plane wave plasmon modes, in this work we draw our attention to the insertion of a metallic nanowire that allows a sensor with strong localized plasmonic modes.
The existence of strong localized optical plasmonic modes introduce strong peaks in the loss spectrum and is a relevant element of the fiber sensor. These localized optical plasmonic modes are strongly guided and this light can be taken from the sensing section of the sensor, which results in a good coupling with conventional optical fibers. When the fiber sensor is integrated in a remote sensing system or in a sensing network this is of special importance. In order to emphasize the role of the plasmonic modes, it is assumed that guided modes result from conventional step index cores. Figure 4 [42] illustrates a novel design of a plasmonic D-type fiber sensor that integrates a gold nanowire mounted on the flat surface, substituting the conventional metallic film. As referred, the optical fiber has a conventional step index core (n co ), made of silica doped with germanium (16%) and enclosed by a cladding of pure silica (n cl ). The Drude model was used to calculate the gold effective refractive index (n m ). In order to improve optical coupling to both the external medium (n ext ) and to the core, as well as to introduce mechanical resistance to the plasmonic sensing configuration, the metallic nanowire was partially embedded in the fiber surface. The gold wire radius was r m = 300 nm, the core radius r co = 1 µm and the distance between the fiber core axis and the metal axis was d = 2 µm.
The PCF D-type fiber sensor, with dm = 45 nm, presents a better performance, enhancing the sensitivity from 2.8 × 10 3 to 10.2 × 10 3 nm/RIU and the resolution from 3.6 × 10 −5 to 9.8 × 10 −6 RIU. The relevance of engineering the structure is confirmed by the results. The optimization of the available light energy near the metallic film, enhances the plasmonic effect, which in turn improves the sensitivity and resolution of the SPR sensor.

D-Type Fiber Configuration with Metallic Nanowire
Leaving the more conventional plasmonic sensing configuration that integrates a metallic film, addressed in the previous section, which supports plane wave plasmon modes, in this work we draw our attention to the insertion of a metallic nanowire that allows a sensor with strong localized plasmonic modes.
The existence of strong localized optical plasmonic modes introduce strong peaks in the loss spectrum and is a relevant element of the fiber sensor. These localized optical plasmonic modes are strongly guided and this light can be taken from the sensing section of the sensor, which results in a good coupling with conventional optical fibers. When the fiber sensor is integrated in a remote sensing system or in a sensing network this is of special importance. In order to emphasize the role of the plasmonic modes, it is assumed that guided modes result from conventional step index cores. Figure 4 [42] illustrates a novel design of a plasmonic D-type fiber sensor that integrates a gold nanowire mounted on the flat surface, substituting the conventional metallic film. As referred, the optical fiber has a conventional step index core (nco), made of silica doped with germanium (16%) and enclosed by a cladding of pure silica (ncl). The Drude model was used to calculate the gold effective refractive index (nm). In order to improve optical coupling to both the external medium (next) and to the core, as well as to introduce mechanical resistance to the plasmonic sensing configuration, the metallic nanowire was partially embedded in the fiber surface. The gold wire radius was rm = 300 nm, the core radius rco = 1 μm and the distance between the fiber core axis and the metal axis was d = 2 μm. The graph in Figure 5a illustrates the real part of the refractive index for the fundamental mode of the optical fiber and for the plasmon modes of the wire (m = 0, 1, 2). The intersection of these The graph in Figure 5a illustrates the real part of the refractive index for the fundamental mode of the optical fiber and for the plasmon modes of the wire (m = 0, 1, 2). The intersection of these dispersion curves are indicated by squared dots and correspond, in wavelength, to the peaks of the curve of the imaginary part of the effective refractive index (solid blue line). These peaks result from a resonate coupling and are denoted as a supermode SM1, circle dot "A" and supermode SM2, circle dot "B". Distribution of the intensity of the light in the metal wire for the plasmon modes of the wire are the detailed in the inset of Figure 5a. As referred, SM1 results from a high plasmonic interaction, noticeable at 900 nm, where the losses reach a high value (point A of Figure 5a) where the supermode presents a fairly strong dipole (upper image of Figure 5b).
After hybridization with the fiber core mode, the light intensity distribution for the supermodes SM1 and SM2 are shown in Figure 5b [43], revealing the transition between the two supermodes. A solid dipolar moment mode (peak A) and a solid quadrupole moment mode (peak B) are observable. Because of the higher losses and potentially higher sensitivity to external refractive index we turn our attention to peak A (supermode SM1). curve of the imaginary part of the effective refractive index (solid blue line). These peaks result from a resonate coupling and are denoted as a supermode SM1, circle dot "A" and supermode SM2, circle dot "B". Distribution of the intensity of the light in the metal wire for the plasmon modes of the wire are the detailed in the inset of Figure 5a. As referred, SM1 results from a high plasmonic interaction, noticeable at 900 nm, where the losses reach a high value (point A of Figure 5a) where the supermode presents a fairly strong dipole (upper image of Figure 5b). After hybridization with the fiber core mode, the light intensity distribution for the supermodes SM1 and SM2 are shown in Figure 5b [43], revealing the transition between the two supermodes. A solid dipolar moment mode (peak A) and a solid quadrupole moment mode (peak B) are observable. Because of the higher losses and potentially higher sensitivity to external refractive index we turn our attention to peak A (supermode SM1).

SM1 (A)
We have carried out a comparison on the performance of this sensor with a conventional D-type fiber based design, also with d = 2 μm and rco = 1 μm, but with a metal film with 45 nm of thickness, mounted on the flat surface [1,44], in order to assess the advantages of this metallic nanowire configuration.
Therefore, the comparison results are presented in Figure 6a. They show that the metal nanowire supports two resonant peaks, that can be strongly wavelength shifted due to the alteration of the external refractive index under measurement, while the conventional configuration with a metal layer (film) it supports only a single resonant peak.
In Figure 6a, for the nanowire structure and considering the peak in the right-supermode SM1, for higher external refractive indexes, a fairly large wavelength shift is noticed, as well as an increase in the losses. The nanowire design structure in both graphs (Figure 6a,b) demonstrates higher losses, indicating a stronger coupling between the external medium and the plasmon and guided modes, than the conventional design structure with a metallic film. Table 2 presents the results of sensitivity and resolution for various for ranges of next, for the two configurations under comparison. We have carried out a comparison on the performance of this sensor with a conventional D-type fiber based design, also with d = 2 µm and r co = 1 µm, but with a metal film with 45 nm of thickness, mounted on the flat surface [1,44], in order to assess the advantages of this metallic nanowire configuration. Therefore, the comparison results are presented in Figure 6a. They show that the metal nanowire supports two resonant peaks, that can be strongly wavelength shifted due to the alteration of the external refractive index under measurement, while the conventional configuration with a metal layer (film) it supports only a single resonant peak.   In Figure 6a, for the nanowire structure and considering the peak in the right-supermode SM1, for higher external refractive indexes, a fairly large wavelength shift is noticed, as well as an increase in the losses. The nanowire design structure in both graphs (Figure 6a,b) demonstrates higher losses, indicating a stronger coupling between the external medium and the plasmon and guided modes, than the conventional design structure with a metallic film. Table 2 presents the results of sensitivity and resolution for various for ranges of n ext , for the two configurations under comparison. The metallic nanowire design structure presented higher performance in comparison with the conventional metallic film design structure, both in sensitivity and resolution, for all ranges of n ext tested. As an example of performance improvement of using the plasmonic metallic nanowire structure, for n ext = [1.38, 1.40], the metallic nanowire sensor has a sensitivity of 8437 nm/RIU, while the conventional metal film sensor has a sensitivity of 3200 nm/RIU.

D-Type Fiber Configuration with Metamaterial
Localized modes can be found in other materials beyond metallic nanowires. For instance, some forms of metamaterials are made of compactly packed metal-dielectric structures with dimensions lower than the wavelength [13]. In this situation the aggregation of localized modes is so solid that they hybridized collectively, whose modes spread over the whole material. Consequently, the metamaterial acts as a bulk material, exhibiting characteristics that result from the nanoscale structure, and therefore can be tailored to stimulate the coupling with guided modes. Truly, the relevant feature of metamaterials is the capacity of exhibiting unusual optical characteristics that are not encountered in natural materials. For instance, metamaterials may present refractive indexes that stimulate the required boundary conditions to create surface plasmons at determined wavelengths or having optical properties that depend on the presence of certain environmental chemical species [16].
For the simulation modelling, the easiest strategy is to compute the optical properties of an ad-hoc metamaterial, using either numerical models [14], direct experimental measurement results [16] or analytical effective medium theory based models [13]. This permits to place the metamaterial as a bulk material in the numerical model. After that one can find the optimal parameters for the metamaterial that optimize sensor performance, such as size and shape, their relative position and composition of the nanostructures.
In this work, it is considered a D-type PCF refractive index sensor based on SPR with a planar structure, where the conventional metal film is replaced by a metamaterial. The performance analysis and a comparison with the conventional structure is conducted. The metamaterial film is composed of thin layers of alumina and silver, as shown in Figure 7 [45].
The relative thickness of the alumina and silver layers define the metamaterial optical properties. Using an effective medium approach, it is possible to calculate the effective refractive index of the metamaterial. For that, the most commonly used models are the Maxwell-Garnett Theory and the Bruggeman effective medium theory [46]. using the Bruggeman effective medium theory model. The metamaterial films with fill ratios from 100% (the case of a pure silver film) until 50%, display a negative dielectric constant for the vast majority of the spectral range considered. Such condition is necessary to promote the existence of surface plasmons. For a fill ration of 50%, the formed metamaterial approaches the epsilon-near-zero regime, limiting the validation of the effective medium approach. 12 metamaterial that optimize sensor performance, such as size and shape, their relative position and composition of the nanostructures.
In this work, it is considered a D-type PCF refractive index sensor based on SPR with a planar structure, where the conventional metal film is replaced by a metamaterial. The performance analysis and a comparison with the conventional structure is conducted. The metamaterial film is composed of thin layers of alumina and silver, as shown in Figure 7 [45]. The relative thickness of the alumina and silver layers define the metamaterial optical properties. Using an effective medium approach, it is possible to calculate the effective refractive index of the metamaterial. For that, the most commonly used models are the Maxwell-Garnett Theory and the Bruggeman effective medium theory [46]. Figure 8 illustrates the real (Figure 8a) and imaginary (Figure 8b) parts of the dielectric, for several filling ratios concentrations of Al2O3 and Ag, calculated using the Bruggeman effective medium theory model. The metamaterial films with fill ratios from 100% (the case of a pure silver film) until 50%, display a negative dielectric constant for the vast majority of the spectral range considered. Such condition is necessary to promote the existence of surface plasmons. For a fill ration of 50%, the formed metamaterial approaches the epsilon-near-zero regime, limiting the validation of the effective medium approach.  In Figure 9, we compare results of the loss as function of wavelength for the metamaterial (80% of Ag and 20% of Al2O3) and for pure silver (only Ag), when subject to different external refractive indices. The peaks of loss, distinctive characteristic of SPR are intensely affected by the external refractive index when the metamaterial layer is considered, while for the case of pure silver they are not vividly affected by the external refractive index change. In Figure 9, we compare results of the loss as function of wavelength for the metamaterial (80% of Ag and 20% of Al 2 O 3 ) and for pure silver (only Ag), when subject to different external refractive indices. The peaks of loss, distinctive characteristic of SPR are intensely affected by the external refractive index when the metamaterial layer is considered, while for the case of pure silver they are not vividly affected by the external refractive index change. Figure 8. Complex dielectric (a) real part and (b) imaginary part of the metamaterial for concentrations of: 50% (small dashed curve), 60% (dash-double dotted curve), 70% (dash-dotted curve), 80% (dotted curve), 90% (dashed curve), and Ag of 100% (solid curve, corresponding to pure silver).
In Figure 9, we compare results of the loss as function of wavelength for the metamaterial (80% of Ag and 20% of Al2O3) and for pure silver (only Ag), when subject to different external refractive indices. The peaks of loss, distinctive characteristic of SPR are intensely affected by the external refractive index when the metamaterial layer is considered, while for the case of pure silver they are not vividly affected by the external refractive index change. From the light distribution present in Figure 9b, it is possible to state that light is less concentrated in the core for the metamaterial film, which permits the formation of a stronger plasmon, promising in a clear way better results when the silver film is substituted by the metamaterial film. Table 3 presents, for different external refractive indices, respectively, the sensor From the light distribution present in Figure 9b, it is possible to state that light is less concentrated in the core for the metamaterial film, which permits the formation of a stronger plasmon, promising in a clear way better results when the silver film is substituted by the metamaterial film. Table 3 presents, for different external refractive indices, respectively, the sensor sensitivity and sensor resolution for pure silver and for the metamaterial with a concentration of 80%Ag and 20%Al 2 O 3 . Table 3. Sensor sensitivity and sensor resolution for the D-type SPR refractive index optical fiber sensor with two different films, silver and the metamaterial with 80%Ag and 20%Al 2 O 3 , respectively. The results show clearly the possibility of performance improvement, in sensitivity and resolution, of SPR optical fiber sensors when using metamaterials. The metamaterial is able to act as an artificial metal, whose optical properties are not found in natural materials. Moreover, the metamaterial can be tailored to optimize the sensor performance characteristics.

D-Type Fiber Configuration with Internal and External Metallic Nanowires for Simultaneous Measurements
Evidently, the guided mode can couple to more than just a single localized plasmonic mode. Therefore, a plasmonic fiber sensor can integrate structures whose respective localized plasmonic modes are sensitive to distinct analytes or parameters, or the distinct analytes or parameters can couple with the guided mode at distinct wavelengths, permitting the development of sensors capable of detecting multiparameters. Due to the intrinsic nature of the plasmonic resonant coupling mechanism, it is important to assure the coupling between the localized plasmonic modes and the guided mode occurs at a suitably wavelength separation, in order to minimize any possible cross-talk and guarantee independent detection.
However, it is important to state that it is possible to couple the guided mode and different localized plasmon modes at the same wavelength (a form of degenerated mode coupling or hybridization) to enhance the sensing mechanism. In this situation, the individual plasmon modes are hybridized and form cooperative plasmon modes with different optical properties that couple with guided mode. This cooperative plasmon modes can reinforce their dependence on the external parameters and augment the sensor performance. One possible example is the configuration shown in sub-section "D-type fiber configuration with metallic nanowire" where the single nanowire would be accompanied by several nanowires, placed adjacently to the original nanowire.
In this sub-section it is analyzed a sensing configuration composed by an engineered D-type fiber with a set of metallic nanowires, where the generated surface plasmon resonances allow to measure refractive index and temperature simultaneously [47]. Figure 10 presents the configuration of the sensor which is based on a D-type fiber with three gold nanowires mounted in the fiber cladding, around the core and one gold nanowire implanted on the flat surface. This latter nanowire is partially embedded on the fiber cladding flat surface to enhance mechanical resistance and to allow optical coupling to the external medium, and therefore exposed to the external refractive index. The three inner metal nanowires are not exposed to the external refractive index and therefore their localized plasmon modes will be sensitive only to temperature.
The important geometrical parameters are the following: distance between the external metallic nanowire center and the fiber core center (d = 2 µm), the core radius (r co = 1 µm), the internal metallic nanowires radii (r I = 500 nm) and the metallic external nanowire radius (r E = 300 nm). Figure 10. Configuration of the proposed SPR D-type fiber with one external and three internal metallic nanowires for simultaneous measurement. Figure 11a illustrates the guided mode loss, per unit of length, versus wavelength for distinct external refractive index values, revealing three main peaks. The first peak at low wavelengths is extremely weak (λP1, present at 700 nm with next = 1.38) and it is not useful for sensing. The higher peak is the second one (λP2, present at 900 nm, with next = 1.38) and shifts towards longer wavelengths when the external refractive index increases, and finally the third peak (λP3, present at 950 nm for next = 1.38) remains practically unchanged when the external refractive index changes. Therefore, this third peak is the adequate peak for temperature measurements. This third peak is weaker in magnitude than the second one, a property that can be tailored either by altering the number of internal nanowires, or by changing the distance between the core and the internal nanowires. Another option is to lower the second peak intensity by increasing the distance between the external wire and the core. Figure 11b illustrates the loss versus wavelength for different temperature values. The left peak corresponds to the second peak (λP2) present in Figure 11a, while the right peak corresponds to the third peak (λP3) present in Figure 11a. The wavelength of λP2 changes with the external refractive  Figure 11a illustrates the guided mode loss, per unit of length, versus wavelength for distinct external refractive index values, revealing three main peaks. The first peak at low wavelengths is extremely weak (λP 1 , present at 700 nm with n ext = 1.38) and it is not useful for sensing. The higher peak is the second one (λP 2 , present at 900 nm, with n ext = 1.38) and shifts towards longer wavelengths when the external refractive index increases, and finally the third peak (λP 3 , present at 950 nm for n ext = 1.38) remains practically unchanged when the external refractive index changes. Therefore, this third peak is the adequate peak for temperature measurements. This third peak is weaker in magnitude than the second one, a property that can be tailored either by altering the number of internal nanowires, or by changing the distance between the core and the internal nanowires. Another option is to lower the second peak intensity by increasing the distance between the external wire and the core. Figure 11b illustrates the loss versus wavelength for different temperature values. The left peak corresponds to the second peak (λP 2 ) present in Figure 11a, while the right peak corresponds to the third peak (λP 3 ) present in Figure 11a. The wavelength of λP 2 changes with the external refractive index and remains practically unchanged with the temperature. On the contrary, λP 3 does not depend on the external refractive index, and for higher temperatures shows a shift toward longer wavelengths.

15
peak is the second one (λP2, present at 900 nm, with next = 1.38) and shifts towards longer wavelengths when the external refractive index increases, and finally the third peak (λP3, present at 950 nm for next = 1.38) remains practically unchanged when the external refractive index changes. Therefore, this third peak is the adequate peak for temperature measurements. This third peak is weaker in magnitude than the second one, a property that can be tailored either by altering the number of internal nanowires, or by changing the distance between the core and the internal nanowires. Another option is to lower the second peak intensity by increasing the distance between the external wire and the core. Figure 11b illustrates the loss versus wavelength for different temperature values. The left peak corresponds to the second peak (λP2) present in Figure 11a, while the right peak corresponds to the third peak (λP3) present in Figure 11a. The wavelength of λP2 changes with the external refractive index and remains practically unchanged with the temperature. On the contrary, λP3 does not depend on the external refractive index, and for higher temperatures shows a shift toward longer wavelengths.  In Table 4 it is presented the sensitivity coefficients of the for refractive index and temperature wavelength peaks, clearly indicating the capacity of this D-type fiber sensor configuration, with gold nanowires, to simultaneously measure refractive index and temperature, based on localized plasmons. This work has proved the possibility to engineer the fiber structure with metal wires to generate localized plasmon modes in the wires, that result in two distinct measurement peaks, capable of simultaneously measuring two distinct physical parameters. In principle, this concept can be further extended to measure even more parameters.

Future Perspectives
The continuous development of fabrication techniques allowed to conceive a new generation of nanostructured optical fiber sensors, that explore a multitude of new degrees of freedom and control of the operation of these devices. The size, shape and material constitution fix the optical characteristics of the localized modes supported by the metal-dielectric structures, while their relative position controls how they couple and hybridize to determine the optical response of the sensor to external parameters or analytes. In this article, we have shown a few examples of how to explore this tremendous potential using numerical tools, but much of the possibilities and combinations have only now begun to be explored and much work is still to be performed.
It has been shown how photonic crystal fibers can be used to generate well localized and guided optical modes that carry light into the sensor and provide a readout channel. In general, for sensors based on the measurement of losses, we are looking for guided modes with very low losses and therefore optical confinement is of crucial importance. Then, at specific wavelengths it is possible to introduce loss peaks, via the resonant coupling with other modes, either localized plasmonic modes or collective modes that are responsive to the parameters being monitored. The challenge here is to engineer these modes to simultaneously create a strong dependency on the external parameters while still satisfying the necessary conditions for a resonant coupling with the localized guided modes. Clearly, the range of structural parameters of the sensor is so wide and the costs of fabrication so large that producing and testing each design is unpractical, leaving numerical solutions as the main tool for sensor performance optimization.
In this review paper, the focus was on a specific subset of nanostructured fiber optical sensors, namely those that have a metal-dielectric nanostructured cross-section uniform throughout the length of the sensing section of the optical fiber and that operate via the monitoring of the optical losses. In these cases, the problem can be reduced to the analysis of a two-dimensional eigen-value problem associated with the modes supported by the device, and ultimately in the determination of the imaginary part of the effective refractive index for each mode. Clearly, the FEM constitutes the ideal method to solve numerically this problem by balancing accuracy with economy of computer resources, specifically computer memory and run time. This approach can be extended to nanostructures much more complex than those presented here, but they constitute examples of the phenomena at hand and how to employ them in sensor design using numerical tools. Specifically, it was exemplified the use of nanowires, metamaterials and photonic crystal fibers, but much more is possible.
These basic principles described in this paper can be adapted to consider more complex situations. One can conceive three dimensional structures that introduce non-uniform distributions of refractive index throughout the length of the optical fiber. These problems will be computationally much more intensive, especially if these structures do not have a repetitive pattern that allow to divide the fiber into cell structures and employ FEM with periodic boundary conditions, much alike the treatment given to study of crystals. Then, it will be necessary to combine beam propagation methods to fully account for the evolution of the optical signal as it propagates through the fiber.
On the other hand, the micro and nanoscale of these structures forces the discretization of the optical modes supported and in some sense these modes become quantized, of which the similarity with the hybridization of atomic orbitals during the formation of a chemical bond constitutes an intuitive clue. In this case, the system is not truly quantum but behaves much like a quantum system. This similarity suggests that it might be possible to combine these structures with truly quantum effects, such mode squeezing, to enhance sensing sensitivity. This would require a completely new approach to the experimental methods used to operate the sensor and necessarily in the modeling tools employed.
Also, it is important to notice that plasmons are inherently nonlinear excitations and as we consider structures with sizes and separations of just a few tens of nanometers these effects might be noticeable, even at low optical powers. Capturing these effects and employing them in optical sensing will require yet a new challenge into numerical modeling. To provide an example, if instead of using the refractive index to model the properties of the metal, we use a fluid description for the motion of the free or conduction electron in the material then, we would obtain necessarily a system of coupled differential equations that are inherently nonlinear.
In all these situations, whether we consider simpler or more complex nanostructured optical sensors, or even quantum and nonlinear effects, the exploration in design and engineering of these sensors requires numerical simulation tools of ever-increasing sophistication. It constitutes a challenge for both those who advance new sensors and to those who develop these simulation tools.