Fano Resonance Sensor with Ultra-High Spectral Resolution in a Metallic Waveguide

1. Introduction

Optical sensors based on interferometry represent the gold standard for high-precision metrology, leveraging the principle of superposition to detect minute changes in phase—a quantity far more sensitive than simple light intensity or wavelength shifts[1,2,3]. Traditional, simple interferometers, such as the Mach-Zehnder[4], Michelson[5], and Sagnac[6] configurations, are fundamental but often limited to measurement accuracies bounded by the wavelength of light. To push sensitivity further, state-of-the-art devices employ high-finesse resonance, where light undergoes controlled multiple interference and scattering. Examples include Fiber Bragg Gratings (FBG) [7] and Fabry-Pérot (FP) interferometers [8]. In these devices, increased mirror reflectivity (FP) or strong grating structures (FBG) narrow the spectral transmission bandwidth, leading to high-resolution sensing. The performance limit of these systems is often constrained by the mirror reflectivity and resonator quality (attenuations, flatness, and parallelism), requiring highly complex and expensive engineering. This complexity is exemplified by the Laser Interferometer Gravitational-Wave Observatory (LIGO) project, which successfully measured gravitational strains of the order of 10^-21 using highly refined interferometric techniques[9,10]. While demonstrating extraordinary accuracy, such systems are prohibitively expensive for common applications.

This gap—the need for high performance coupled with low complexity—drives the search for compact, integrated optical sensors capable of ultra-high spectral resolution. We propose a novel approach that achieves extreme sensitivity not through increased reflectivity, but by utilizing Fano resonances [11,12]. Fano resonances offer intrinsically sharp, asymmetric spectral profiles derived from the interference between a localized resonant state (a "bound state") and a continuous spectrum. However, realizing Fano resonances in standard optical systems is challenging because the required bound states, analogous to those in quantum mechanics, do not naturally exist in ordinary optical continuum settings.

To circumvent this fundamental problem, we introduce a new optical architecture: a narrow waveguide with metallic boundaries operating close to its cutoff frequency (similar motivation to work near the cutoff frequency was taken in [13]). By placing a local change in the refractive index (a dielectric impurity) inside this waveguide, we effectively induce a resonant state whose "energy" is lower than the waveguide cutoff. This combination of the local impurity and the strongly confining waveguide creates a scattering environment that yields the necessary Fano interference. A key physical outcome is that at a specific incoming light frequency, the device acts as a perfect reflector, exhibiting exactly zero transmission.

We demonstrate that the spectral bandwidth of this perfect reflection feature can be made exceptionally narrow, simply by reducing the size of the impurity or increasing its proximity to the waveguide boundary. Crucially, this ultra-narrow bandwidth does not affect the depth of the resonance; the point of zero transmission is preserved. This unique configuration allows the device to function as an exceptionally sensitive sensor, capable of measuring light frequency or physical strain with unparalleled precision. Analytically, we employ the Impurity D Function (IDF) [14,15] to model the sub-wavelength dielectric scatterer. To the best of our knowledge, this work represents the first application of the IDF—traditionally used for point potentials in Quantum Mechanics—to the domain of optical systems. Through this analysis, we provide evidence that the spectral resolution of this compact device can exceed 10²⁰, allowing for the measurement of strains smaller than 10^-21.

The remainder of this paper is structured as follows: Section 2 details the theoretical model of the device. Section 3 calculates the device’s transmission using the IDF formalism. Section 4 analyzes the proximity and behavior of the zero-transmission point. Section 5 quantifies the device’s spectral resolution and overall performance, and Section 6 provides a summary of the main results.

2. Theory and Model

The system configuration is illustrated in Figure 1. z-polarized electromagnetic (EM) wave propagates along the x-axis, confined in the y-direction by a metallic waveguide of width w. A local impurity, defined as a perturbation in the medium's dielectric coefficient, is positioned at a distance ε from the metallic boundary.

The electromagnetic field wave equation reads

$$\left({\displaystyle \frac{{\partial }^2}{\partial x^2}}+{\displaystyle \frac{{\partial }^2}{\partial y^2}}\right)E\left(x,y,t\right)-{\displaystyle \frac{n^2\left(x,y,t\right)}{c^2}}{\displaystyle \frac{{\partial }^2}{\partial t^2}}E\left(x,y,t\right)=0$$

(1)

Assuming a monochromatic incident beam with angular frequency $\mathsf{\omega }$, we substitute the ansatz

$$E\left(x,y,t\right)=zE\left(x,y\right)\exp \left(-i\omega t\right),$$

into the time-dependent Maxwell's equations. This substitution simplifies the system into the time-independent Helmholtz equation for the spatial field $E\left(x,y\right)$:

$$\left({\displaystyle \frac{{\partial }^2}{\partial x^2}}+{\displaystyle \frac{{\partial }^2}{\partial y^2}}\right)E\left(x,y\right)+k^2\left(x,y\right)E\left(x,y\right)=0$$

(2)

where $k^2\left(x,y\right)\equiv {\displaystyle \frac{n^2\left(x,y\right)}{c^2}}{\omega }^2$,

and the dielectric coefficient is homogenous inside the waveguide except for the local impurity (see Fig.1), i.e.,

$$n^2\left(r\right)=\left\{\begin{array}{cc}{n}_0^2+\Delta n^2& \left|r-\mathsf{\epsilon }\widehat{y}\right|<R\\ n_0^2& \left|r-\mathsf{\epsilon }\widehat{y}\right|<R,\text{}0<y<w\end{array}\right.$$

(3)

Eq.(2) can be rewritten

$$-\left({\displaystyle \frac{{\partial }^2}{\partial x^2}}+{\displaystyle \frac{{\partial }^2}{\partial y^2}}\right)E\left(x,y,t\right)-\Delta k^2\left(r\right)\left(x,y\right)E\left(x,y,t\right)=k_0^{2}E\left(x,y,t\right)$$

(4)

where

$$\Delta k^2\left(r\right)=\left\{\begin{array}{cc}\Delta n^2{k}_0^2& \left|r-\mathsf{\epsilon }\widehat{y}\right|<R\\ 0& else\end{array}\right..$$

(5)

We assume that the impurity’s dimensions are significantly smaller than the operating wavelength, enabling a point-like approximation, i.e., $R<<\lambda$ or $R{k}_0<<1$, where $k_0\equiv 2\pi n_0/\lambda$, $k_0\equiv 2\pi /{\lambda }_0$. Hence, we can replace $\Delta k^2\left(r-r_0\right)$ with the asymmetric Impurity D Function (IDF) [16,17]

$$D\left(r\right)=-\underset{\mathsf{\rho }\to 0}{\lim }{\displaystyle \frac{2\sqrt{\pi }\delta \left(x\right)\exp \left(-{\left(y/\mathsf{\rho }\right)}^2\right)}{\mathsf{\rho }\ln \left({\mathsf{\rho }}_0/\mathsf{\rho }\right)}}$$

(6)

where (see explanation in Appendix A)

$${\mathsf{\rho }}_0=2Rexp\left({\displaystyle \frac{2}{R^2\Delta k^2}}+\gamma /2-{\displaystyle \frac{1}{2}}\right).$$

(7)

3. Generic Solution

Since the impurity is approximated as a point scatterer, the full field solution can be derived using the Green’s function $G\left(r,r_0\right)$ of the empty waveguide (the unperturbed system), and since the Green function solves

\

$\left({\nabla }^2+k_0^2\right)G\left(r,r_0\right)=-\delta \left(r-r_0\right)$ for $0<y<w$,

the overall solution can be formulated as a sum of the incident wave $E_{in}\left(r\right)$ and the scattered wave, where the latter is determined by the Green’s function of the unperturbed system, i.e., [16,17]

$$E\left(r\right)=E_{in}\left(r\right)-{\displaystyle \frac{E_{in}\left(r_0\right){\displaystyle \int dr\text{'}D\left(r\text{'}\right)}}{1+{\displaystyle \int dr\text{'}D\left(r\text{'}-r_0\right)G\left(r\text{'},r_0\right)}}}G\left(r,r_0\right)$$

(8)

where in case of a waveguide geometry, the Green function is

$$G\left(r,r\text{'}\right)=i{\displaystyle \sum _{n=1}^{\infty }\frac{\sin \left(n\pi y/w\right)\sin \left(n\pi y\text{'}/w\right)}{w\sqrt{k_0^2-{\left(n\pi /w\right)}^2}}\exp \left(i\sqrt{k_0^2-{\left(n\pi /w\right)}^2}\left|x-x\text{'}\right|\right)}.$$

(9)

The incident wave reads

$$E_{in}\left(r\right)=E_0\sin \left(\pi y/w\right)\exp \left(i\sqrt{k_0^2-{\left(\pi /w\right)}^2}x\right).$$

(10)

The overall solution is therefore

$$E_{sc}\left(r\right)=E_0{\displaystyle \sum \sin \left(\frac{\pi qy}{w}\right)\left({\delta }_{1q}\exp \left(i\sqrt{k^2-{\left(\frac{q\pi }{w}\right)}^2}x\right)-A_q\exp \left(i\sqrt{k^2-{\left(\frac{q\pi }{w}\right)}^2}\left|x\right|\right)\right)},$$

(11)

where ${\delta }_{pq}=\left\{\begin{array}{cc}1& p=q\\ 0& p\ne q\end{array}\right.$ is the Kronecker delta function and the coefficient $A_q$ reads (see Appendix B)

$$A_q=\underset{\mathsf{\rho }\to 0}{\lim }{\displaystyle \frac{\sin \left(\pi \mathsf{\epsilon }/w\right)\sin \left(\pi \mathsf{\epsilon }q/w\right)}{iw\sqrt{k^2-{\left(q\pi /w\right)}^2}\left[\frac{\ln \left({\mathsf{\rho }}_0/\mathsf{\rho }\right)}{2\pi }+{\displaystyle \sum _{n=1}^{\infty }\frac{{\sin }^2\left(\pi \mathsf{\epsilon }n/w\right)}{iw\sqrt{k^2-{\left(n\pi /w\right)}^2}}\exp \left(-{\left(\frac{n\pi \mathsf{\rho }}{2w}\right)}^2\right)}\right]}}.$$

(12)

Therefore, beyond the point scatterer, i.e., for $x>0$, the scattered field is

$$E_{sc}\left(r\right)=E_0{\displaystyle \sum \left({\delta }_{1q}-A_q\right)}\sin \left(\pi qy/w\right)\exp \left(i\sqrt{k^2-{\left(q\pi /w\right)}^2}x\right).$$

(13)

In the operational regime where ${\left(\pi /w\right)}^2<k_0^2<{\left(2\pi /w\right)}^2$, only the fundamental mode is propagating. Although the coefficients for the higher-order modes are non-zero, these modes manifest as evanescent waves whose longitudinal power flow vanishes. Consequently, these modes do not contribute to energy transport far from the impurity. The overall transmission amplitude for the fundamental mode is given by the expression:

$$t_{11}=1-\underset{\mathsf{\rho }\to 0}{\lim }{\displaystyle \frac{{\sin }^2\left(\pi \mathsf{\epsilon }/w\right)}{iw\sqrt{k^2-{\left(\frac{\pi }{w}\right)}^2}\left[\frac{\ln \left({\mathsf{\rho }}_0/\mathsf{\rho }\right)}{2\pi }+{\displaystyle \sum _{n=1}^{\infty }\frac{{\sin }^2\left(\pi \mathsf{\epsilon }n/w\right)}{iw\sqrt{k^2-{\left(n\pi /w\right)}^2}}\exp \left(-{\left(\frac{n\pi \mathsf{\rho }}{2w}\right)}^2\right)}\right]}}$$

(14)

In the limit $\mathsf{\rho }\to 0$, the divergence of the summation in the denominator (proportional to $\ln \mathsf{\rho }$) is precisely canceled by the first logarithmic term in the analytic expression, thus ensuring the denominator remains finite in this limit.

To simplify the subsequent analysis, we define the parameter $\overline{\mathsf{\rho }}$ as the analytic continuation of the limiting process:

$$\underset{\mathsf{\rho }\to 0}{\lim }\left\{{\displaystyle \frac{\ln \left(\mathsf{\rho }/\overline{\mathsf{\rho }}\right)}{2\pi }}-{\displaystyle \sum _{n=3}^{\infty }\frac{{\sin }^2\left(\pi \mathsf{\epsilon }n/w\right)}{iw\sqrt{k^2-{\left(n\pi /w\right)}^2}}\exp \left(-{\left(\frac{n\pi \mathsf{\rho }}{2w}\right)}^2\right)}\right\}=0.$$

(15)

Finally, the transmission coefficient can be written

$$t_{11}={\displaystyle \frac{1}{1+\frac{{\sin }^2\left(\pi \mathsf{\epsilon }/w\right)}{iw\sqrt{k^2-{\left(\pi /w\right)}^2}}{\left[\frac{\ln \left({\mathsf{\rho }}_0/\overline{\mathsf{\rho }}\right)}{2\pi }-\frac{{\sin }^2\left(\pi \mathsf{\epsilon }2/w\right)}{w\sqrt{{\left(2\pi /w\right)}^2-k^2}}\right]}^{-1}}}.$$

(16)

4. Analysis–Perfect Reflection Scenario

There exists a specific resonant wavenumber, $k_{*}$ (or optical wavelength), for which the following condition, defining the point of perfect destructive interference, is satisfied:

$${\displaystyle \frac{\ln \left({\mathsf{\rho }}_0/\overline{\mathsf{\rho }}\right)}{2\pi }}={\displaystyle \frac{{\sin }^2\left(\pi \mathsf{\epsilon }2/w\right)}{w\sqrt{{\left(2\pi /w\right)}^2-k_{*}^2}}}$$

(17)

At this precise resonance point, the transmission vanishes exactly ($t_{11}=0$).

Fig. 2 shows the resulting field intensity profile of the exact solution inside the waveguide, specifically evaluated at the resonant wavenumber $k_{*}$.

At the resonant wavenumber, the fundamental (first) mode is perfectly reflected. While higher-order modes are excited, they are evanescent waves that decay exponentially away from the impurity and therefore do not contribute to energy transport in the far field. The coefficients of these non-propagating modes, $A_q$, are

$$A_q={\displaystyle \frac{\sqrt{k_{*}^2-{\left(\pi /w\right)}^2}}{\sqrt{k_{*}^2-{\left(q\pi /w\right)}^2}}}{\displaystyle \frac{\sin \left(\pi \mathsf{\epsilon }q/w\right)}{\sin \left(\pi \mathsf{\epsilon }/w\right)}}$$

(18)

In the regime where the impurity is a surface defect, i.e., $\mathsf{\epsilon }/w<<1$, $k_{*}^2\to {\left(2\pi /w\right)}^2$, in which case

$$A_q\cong \left\{\begin{array}{cc}{\displaystyle \frac{q\sqrt{3}}{\sqrt{4-q^2}}}& q>2\\ {\displaystyle \frac{\sqrt{3}\ln \left({\mathsf{\rho }}_0/\overline{\mathsf{\rho }}\right)}{i{\left(2\pi \mathsf{\epsilon }/w\right)}^2}}& q=2\end{array}\right.$$

(19)

Furthermore, in this surface-defect regime (see Appendix C),

$$\overline{\mathsf{\rho }}\to C_5\mathsf{\epsilon }.$$

(20)

Therefore, the second mode provides the largest contribution to the scattered field. The full field solution can be simplified to a superposition involving only the incident/reflected fundamental mode and this dominant scattered mode:

$$E_{sc}\left(r\right)\cong E_0\left[\sin \left({\displaystyle \frac{\pi y}{w}}\right)\sin \left(\pi \sqrt{3}{\displaystyle \frac{x}{w}}\right)u\left(-x\right)+i{\displaystyle \frac{\sqrt{3}}{d}}\sin \left({\displaystyle \frac{2\pi y}{w}}\right)\exp \left(-2\pi d{\displaystyle \frac{\left|x\right|}{w}}\right)\right]$$

(21)

where

$$d\equiv {\displaystyle \frac{{4\left(\pi \mathsf{\epsilon }/w\right)}^2}{\ln \left({\mathsf{\rho }}_0/C_5\mathsf{\epsilon }\right)}}<<1$$

(22)

and

$u\left(x\right)=\left\{\begin{array}{cc}0& x\le 0\\ 1& x>0\end{array}\right.$ is the Heaviside step function.

Therefore, in this regime, the transmission coefficient of the first mode is

$$t_{11}\cong {\displaystyle \frac{1}{\frac{1}{i2\sqrt{3}}{\left\{\frac{1}{d}-\frac{1}{\sqrt{1-{\left(w/\lambda \right)}^2}}\right\}}^{-1}+1}}.$$

(23)

5. Spectral Resolution

The extreme sensitivity of this device stems from the proximity of the zero-transmission point (${\lambda }_{*}$) and the half-transmission point (${\lambda }_{1/2}$).The exact resonant wavelength, ${\lambda }_{*}$, where the transmission vanishes, is determined by the waveguide width

$${\lambda }_{*}={\displaystyle \frac{w}{\sqrt{1-d^2}}}\cong w\left(1+{\displaystyle \frac{d^2}{2}}+{\displaystyle \frac{3d^4}{8}}\right).$$

(24)

Furthermore, the wavelength at which the transmission coefficient equals one-half is derived using a small parameter expansion:

$${\lambda }_{1/2}\cong w\left(1+{\displaystyle \frac{d^2}{2}}+{\displaystyle \frac{d^3}{2\sqrt{3}}}\right)$$

(25)

Since $d<<1$, the difference between ${\lambda }_{1/2}$ and ${\lambda }_{*}$ is extremely small. The resulting spectral resolution ($\Re$) defined as

$$\Re \equiv {\displaystyle \frac{\lambda }{\Delta \lambda }}={\displaystyle \frac{{\lambda }_{*}}{2\left({\lambda }_{1/2}-{\lambda }_{*}\right)}}\cong {\displaystyle \frac{\sqrt{3}}{d^3}}$$

(26)

This demonstrates that the device's spectral resolution is dramatically enhanced simply by optimizing the parameter d.

Fig. 3 illustrates the transmission coefficient $T_{11}={\left|t_{11}\right|}^2$ as a function of wavelength for two distinct waveguide widths, demonstrating the sensor's response to strain. The observed spectral shift resulting from a simulated strain of $\mathsf{\eta }={10}^{-5}$ is clearly resolvable in the vicinity of the resonant wavelength. In this example, the impurity is set at $\mathsf{\epsilon }/w=0.1$. The overall device sensitivity can be further enhanced by substantially reducing the ratio $\mathsf{\epsilon }/w$, effectively positioning the impurity closer to the metallic boundary.

We identify the half-transmission point ($\lambda ={\lambda }_{1/2}$) as the optimal working point. At this wavelength, the transmission coefficient exhibits a strong linear relationship with respect to both the wavelength and the waveguide width. Placing the device at this point ensures that the change in transmission is directly proportional to the mechanical strain ($\mathsf{\eta }=\Delta w/w$) experienced by the waveguide. In this highly sensitive regime, the transmission is accurately approximated by the linear formula:

$$T_{11}\equiv {\displaystyle \frac{1}{2}}+{\displaystyle \frac{\Delta w}{w}}{\displaystyle \frac{\sqrt{3}}{d^3}}={\displaystyle \frac{1}{2}}+\Re \mathsf{\eta }.$$

(27)

where the spectral resolution $\Re$ serves as the large proportionality constant defining the strain sensitivity.

Beyond its primary function as a standalone gauge, the sensor can be readily integrated as a highly sensitive optical element within established differential interferometric systems. This includes configurations such as Michelson, Mach-Zehnder, or Sagnac setups, thereby facilitating the measurement of relative strains between orthogonal directions.

Figure 4. Spectral resolution vs. proximity to the boundary $\mathsf{\epsilon }$. The different curves represents different values of $\Delta n^2$ (the values shown in the legend). In these cases it was taken that $\mathsf{\epsilon }=R$.

Figure 4. Spectral resolution vs. proximity to the boundary $\mathsf{\epsilon }$. The different curves represents different values of $\Delta n^2$ (the values shown in the legend). In these cases it was taken that $\mathsf{\epsilon }=R$.

According to (7), in the regime of a small impurity, i.e., $R/\lambda <<1$ , $\ln \left({\mathsf{\rho }}_0/C_5\mathsf{\epsilon }\right)~2R^{-2}\Delta k^{-2}$, in which case

$$\Re \equiv {\displaystyle \frac{\lambda }{\Delta \lambda }}~8\sqrt{3}{\left(4{\pi }^2{\displaystyle \frac{\mathsf{\epsilon }R}{\lambda w}}\right)}^{-6}{\left(\Delta n^2\right)}^{-3}$$

(28)

In case $\mathsf{\epsilon }~R$ and $\lambda ~w$

$$\Re ~8\sqrt{3}{\left(2\pi {\displaystyle \frac{\mathsf{\epsilon }}{w}}\right)}^{-12}{\left(\Delta n^2\right)}^{-3}$$

(29)

5.1. Theoretical Resolution and Practical Scaling

Since, in principle, the impurity size can be arbitrarily small and situated arbitrarily close to the boundary, the dimensionless parameter d can be made vanishingly small. Because the spectral resolution $\Re$ scales inversely with the cube of this ratio ($\Re \cong \sqrt{3}/d^3$), the theoretical resolution is potentially enormous.

It must be emphasized that while the interferometer’s resolution $\Re$ is extremely high, the waveguide length must be larger than its width (w) by a factor proportional to $d^{-1}$. To illustrate this scaling advantage, consider the goal of measuring strains on the order of $~{10}^{-18}$: this requires a spectral resolution of $\Re \cong {10}^{18}$. The corresponding required length-to-width ratio would be six orders of magnitude larger than w. Thus, for a $1\mu m$ width waveguide, the necessary length would only be on the order of a meter. Such a compact device, with a spectral resolution of $\Re \cong {10}^{18}$, could theoretically measure strains down to ${10}^{-21}$, assuming current measurement precision (e.g., transmission precision better than 0.1%).

5.2. Construction and Experimental Requirements

In general, the construction imposes few constraints on the local impurity, requiring only that it functions as a subwavelength scatterer whose refractive index exceeds that of the surrounding waveguide medium. This makes the impurity feasible to construct, potentially utilizing existing techniques such as nanofibers [18,19,20,21] situated along the z-direction. However, achieving such high sensitivity experimentally requires stringent controls. The exact implementation must ensure a vibration-free environment, perfectly smooth waveguide boundaries, and the use of low-noise laser sources. The precise engineering methodology required to realize this ultra-sensitive device warrants extended research.

6. Summary and Conclusions

We propose a novel architecture for an ultra-high spectral resolution sensor consisting of a simple, compact metallic-boundary waveguide containing a single point dielectric impurity. A key advantage of this design is its simplicity, as it requires neither complex fine-tuning nor high-finesse mirror systems common in conventional interferometers.

Despite its minimal construction, the device achieves exceptional performance by exploiting Fano resonance, which arises from the quantum-optical interference between the waveguide's continuous modes and a quasi-bound state induced locally by the impurity near the waveguide's cutoff frequency. This mechanism yields a critically sharp spectral feature—a point of perfect reflection (zero transmission)—with an extremely narrow bandwidth.

We demonstrate analytically that the device's spectral resolution ($\Re$) for the case of a small surface defect ($\mathsf{\epsilon }/w<<1$) can be inversely proportional to the 12^th power of the ratio $\mathsf{\epsilon }/w$: $\Re ~{\left(\mathsf{\epsilon }/w\right)}^{-12}$.

This intrinsic high resolution translates directly into extraordinary sensitivity for measuring mechanical strains. When the sensor is calibrated at the 50% transmission working point ($T_{11}\cong 1/2$), the change in transmission becomes linearly proportional to the mechanical strain ($\mathsf{\eta }=\Delta w/w$), with the resolution $\Re$ acting as the proportionality constant: $T_{11}=1/2+\Re \mathsf{\eta }$. Our quantitative analysis shows that for a practical geometric ratio of $\mathsf{\epsilon }/w\cong 1\%$, the theoretical spectral resolution can reach $\Re ~{10}^{20}$. Consequently, this compact, simple device is capable of measuring mechanical strains smaller than $\mathsf{\eta }~{10}^{-21}$, demonstrating performance comparable to large-scale, high-complexity instruments like gravitational wave detectors.