Exact zero transmission during the Fano resonance phenomenon in non symmetric waveguides

TL;DR

This paper proves that in non-symmetric waveguides, there is always a real frequency at which the transmission coefficient is zero, meaning complete backscattering occurs, regardless of symmetry assumptions.

Contribution

It demonstrates that zero transmission can occur in non-symmetric waveguides, extending previous results that required symmetry for such phenomena.

Findings

Existence of a real frequency with zero transmission in non-symmetric waveguides

Numerical validation of the theoretical results

Energy is fully backscattered at the zero transmission frequency

Abstract

We investigate a time-harmonic wave problem in a waveguide. We work at low frequency so that only one mode can propagate. It is known that the scattering matrix exhibits a rapid variation for real frequencies in a vicinity of a complex resonance located close to the real axis. This is the so-called Fano resonance phenomenon. And when the geometry presents certain properties of symmetry, there are two different real frequencies such that we have either $R=0$ or $T=0$, where $R$ and $T$ denote the reflection and transmission coefficients. In this work, we prove that without the assumption of symmetry of the geometry, quite surprisingly, there is always one real frequency for which we have $T=0$. In this situation, all the energy sent in the waveguide is backscattered. However in general, we do not have $R=0$ in the process. We provide numerical results to illustrate our theorems.

Figures

Equations

$$\{(x,y)\in \mathbb{R}\times (0;1)\}$$

$$\begin{array}[]{|rcll}\Delta u+\lambda u&=&0&\mbox{ in }\Omega\\[3.0pt] \partial_{\nu}u&=&0&\mbox{ on }\partial\Omega.\end{array}$$

$$w_{\pm }(x,y)=e^{\pm ix\sqrt{\lambda }}$$

$$u_{+}=\begin{array}[]{|ll}w_{+}+R_{+}\,w_{-}+\dots,&\mbox{ for }x<-d\\ \hskip 31.2982ptT\ w_{+}+\dots,&\mbox{ for }x>d,\end{array}\qquad\qquad u_{-}=\begin{array}[]{|ll}\hskip 28.45274ptT\ w_{-}+\dots,&\mbox{ for }x<-d\\ w_{-}+R_{-}\,w_{+}+\dots,&\mbox{ for }x>d.\end{array}$$

$$\mathfrak{s}:=\left(\begin{array}[]{cc}R_{+}&T\\ T&R_{-}\\ \end{array}\right)\in\mathbb{C}^{2\times 2}.$$

$$u_{\mathrm{tr}}=K{e}^{-x\sqrt{{\pi }^2-{\lambda }^0}}\cos (\pi y)+{\tilde{u}}_{\mathrm{tr}}\text{mbox}forx\ge d,$$

$${\lambda }^{\epsilon }={\lambda }^0+\epsilon {\lambda }^{\prime }$$

$$n(s)=\epsilon H(s).$$

$$\begin{array}[]{|rcll}\Delta u^{\varepsilon}+\lambda^{\varepsilon}u^{\varepsilon}&=&0&\mbox{ in }\Omega^{\varepsilon}\\[3.0pt] \partial_{\nu^{\varepsilon}}u^{\varepsilon}&=&0&\mbox{ on }\partial\Omega^{\varepsilon},\end{array}$$

$$\mathfrak{s}(\epsilon ,\lambda ),T(\epsilon ,\lambda ),R_{+}(\epsilon ,\lambda ),R_{-}(\epsilon ,\lambda )$$

$${\mathfrak{s}}^0:=\mathfrak{s}(0,{\lambda }^0),T^0:=T(0,{\lambda }^0),R_{+}^0:=R_{+}(0,{\lambda }^0),R_{-}^0:=R_{-}(0,{\lambda }^0).$$

$$\kappa (H)$$

$$\alpha$$

$$\beta (H)$$

$$\underset{\epsilon \to 0}{\lim }\mathfrak{s}(\epsilon ,{\lambda }^0+\epsilon {\lambda }^{\prime })={\mathfrak{s}}^0.$$

$$\underset{\epsilon \to 0}{\lim }\mathfrak{s}(\epsilon ,{\lambda }^0+\epsilon \kappa (H)+{\epsilon }^2\mu )={\mathfrak{s}}^0+\frac{{\tau }^{\top }\tau }{i\tilde{\mu }-|\tau {|}^2/2},$$

$$\begin{array}[]{lcl}\mathfrak{s}^{\varepsilon}(\mu)&:=&\mathfrak{s}(\varepsilon,\lambda^{0}+\varepsilon\kappa(H)+\varepsilon^{2}\mu)\\[2.0pt] T^{\varepsilon}(\mu)&:=&T(\varepsilon,\lambda^{0}+\varepsilon\kappa(H)+\varepsilon^{2}\mu)\\[2.0pt] R_{+}^{\varepsilon}(\mu)&:=&R_{+}(\varepsilon,\lambda^{0}+\varepsilon\kappa(H)+\varepsilon^{2}\mu)\\[2.0pt] R_{-}^{\varepsilon}(\mu)&:=&R_{-}(\varepsilon,\lambda^{0}+\varepsilon\kappa(H)+\varepsilon^{2}\mu).\end{array}$$

$$|{\mathfrak{s}}^{\epsilon }(\mu )-s^{\mathrm{asy}}(\mu )|\le C\epsilon \text{mbox}with{s}^{\mathrm{asy}}(\mu )={\mathfrak{s}}^0+\frac{{\tau }^{\top }\tau }{i\tilde{\mu }-|\tau {|}^2/2},$$

$$|T^{\epsilon }(\mu )-T^{\mathrm{asy}}(\mu )|\le C\epsilon \text{mbox}with{T}^{\mathrm{asy}}(\mu )=T^0+\frac{{\tau }_1{\tau }_2}{i\tilde{\mu }-(|{\tau }_1{|}^2+|{\tau }_2{|}^2)/2}.$$

$$\{T^{\mathrm{asy}}(\mu ),\mu \in \mathbb{R}\}={\mathcal{C}}^{\mathrm{asy}}\setminus \{T^0\}$$

$$\frac{|{\tau }_1{|}^2+|{\tau }_2{|}^2}{2}=\Re e\left(\frac{{\tau }_1{\tau }_2}{T^0}\right).$$

$$\overline{U}{\mathfrak{s}}^0=U.$$

$$R_{+}^{0}\,\overline{u_{+}}+T^{0}\,\overline{u_{-}}=\begin{array}[]{|ll}R_{+}^{0}\,(w_{-}+\overline{R_{+}^{0}}\,w_{+})+T^{0}\,(\overline{T^{0}}\,w_{+})+\dots,&\mbox{ for }x<-d\\ R_{+}^{0}\,(\overline{T^{0}}\,w_{-})+T^{0}\,(w_{+}+\overline{R_{-}^{0}}\,w_{-})+\dots,&\mbox{ for }x>d.\end{array}$$

$$\overline{\tau }{\mathfrak{s}}^0=\tau .$$

$$\overline{\tau }{\mathfrak{s}}^0=(\kappa (H)\overline{\alpha }-\overline{\beta (H)})\overline{{\mathfrak{s}}^0}{\mathfrak{s}}^0=\kappa (H)\overline{\alpha }-\overline{\beta (H)}.$$

$$\begin{array}[]{|ll}R_{+}^{0}\,\overline{\tau_{1}}+T^{0}\,\overline{\tau_{2}}=a\\[4.0pt] T^{0}\,\overline{\tau_{1}}+R_{-}^{0}\,\overline{\tau_{2}}=b.\end{array}$$

$$T^0\overline{{\tau }_1}-\frac{\overline{R_{+}^0}{T}^0}{\overline{T^0}}\overline{{\tau }_2}={\tau }_2.$$

$$R_{+}^0=\frac{{\tau }_1-T^0\overline{{\tau }_2}}{\overline{{\tau }_1}}.$$

$$T^0(|{\tau }_1{|}^2+|{\tau }_2{|}^2)-\frac{T^0}{\overline{T^0}}\overline{{\tau }_1{\tau }_2}={\tau }_1{\tau }_2\iff |{\tau }_1{|}^2+|{\tau }_2{|}^2=2\Re e\left(\frac{{\tau }_1{\tau }_2}{T^0}\right).$$

$$C^{\epsilon }=\{T^{\epsilon }(\mu ),\mu \in I^{\epsilon }\}$$

$$-R_{+}^{\epsilon }(\mu )/\overline{R_{-}^{\epsilon }(\mu )}=T^{\epsilon }(\mu )/\overline{T^{\epsilon }(\mu )}\forall \mu \in I^{\epsilon }.$$

Full text

**Exact zero transmission during the Fano resonance

phenomenon in non symmetric waveguides**

Lucas Chesnel1, Sergei A. Nazarov2

1 INRIA/Centre de math matiques appliqu es, École Polytechnique, Universit Paris-Saclay, Route de Saclay, 91128 Palaiseau, France;

2 St. Petersburg State University, Universitetskaya naberezhnaya, 7-9, 199034, St. Petersburg, Russia;

E-mails: lucas.chesnel@inria.fr, srgnazarov@yahoo.co.uk, s.nazarov@spbu.ru

(March 13, 2024)

Abstract. We investigate a time-harmonic wave problem in a waveguide. We work at low frequency so that only one mode can propagate. It is known that the scattering matrix exhibits a rapid variation for real frequencies in a vicinity of a complex resonance located close to the real axis. This is the so-called Fano resonance phenomenon. And when the geometry presents certain properties of symmetry, there are two different real frequencies such that we have either $R=0$ or $T=0$, where $R$ and $T$ denote the reflection and transmission coefficients. In this work, we prove that without the assumption of symmetry of the geometry, quite surprisingly, there is always one real frequency for which we have $T=0$. In this situation, all the energy sent in the waveguide is backscattered. However in general, we do not have $R=0$ in the process. We provide numerical results to illustrate our theorems.

Key words. Waveguides, Fano resonance, zero transmission, scattering matrix.

1 Introduction

The Fano resonance is a universal phenomenon in physics which appears in many areas. For a general presentation, we refer the reader to [18] for the seminal paper and to [27, 26] for recent reviews. In this work, we consider its expression on a model problem of propagation of time-harmonic waves in a waveguide which is unbounded in one direction. This problem appears naturally for instance in acoustics, in water-waves theory or in electromagnetism. In this context, the Fano resonance mechanism can be described as follows. Assume that the Neumann Laplacian (for the problem we consider below) has a real eigenvalue $\lambda^{0}$ embedded in the continuous spectrum. In this case, the corresponding eigenfunctions are the so-called trapped modes which are exponentially decaying at infinity. Then perturbing slightly the setting, for example the geometry or the material index, in general this real eigenvalue will turn into a complex resonance [2, 45, 36]. And for real spectral parameters $\lambda$ (proportional to the square of the frequency) varying in a neighbourhood of $\lambda^{0}$, the scattering matrix will exhibit a rapid variation. This variation is even quicker as the imaginary part of the complex resonance is small. When $\lambda^{0}$ is between the first and the second thresholds in the continuous spectrum, so that only two conjugated waves can propagate in the waveguide, the symmetric scattering matrix is composed of two reflection coefficients $R$, $\tilde{R}$ and one transmission coefficient $T$ (see the notation in (3)). In this case, under certain properties of symmetry of the configuration, one can show that the scattering coefficients take zero values for some real $\lambda$ around $\lambda^{0}$. Such particular values for $R$, $\tilde{R}$ are studied in particular in the context of Perfect Transmission Resonances (PTRs), see e.g. [39, 38, 24, 44, 28]. For the presentation of simple models in optics explaining the Fano resonance phenomenon, we refer the reader to [16, 17]. For more mathematical approaches, one can consult [41, 40, 42, 1, 7]. For computations of complex resonances and numerical investigations of the Fano resonance phenomenon in waveguides, we refer the reader to [12, 6, 13, 19, 21, 20]. For results concerning the existence of trapped modes associated with eigenvalues embedded in the continuous spectrum, see e.g. [43, 14, 15, 11, 25, 29, 34, 35]. Finally, note that another approach to get rigorously a zero transmission coefficient can be found in [8, 9, 10].

The goal of this note is to show that without assumption of symmetry of the configuration, the transmission coefficient $T$ still takes the zero value throughout the Fano resonance phenomenon. This was intuited in [23] using a continuation idea from a symmetric setting. In the present work, we prove rigorously the result using a different approach which does not require to start from a symmetric setting. The outline of the article is as follows. First we present the setting in Section 2. Then we perturb the geometry and the frequency of the configuration supporting trapped modes via a small parameter $\varepsilon>0$ and we recall the results of [7] providing an asymptotic expansion of the scattering matrix with respect to $\varepsilon$ tending to zero. In Section 4, we show that miraculously (we have no physical explanation for that), the main asymptotic term in the expansion of the transmission coefficient passes through zero for real $\lambda$ around $\lambda^{0}$. Then in Section 5, working as in [10], we demonstrate that the unitary structure of the scattering matrix is enough to deduce that the transmission coefficient itself passes through zero for real $\lambda$ around $\lambda^{0}$. We provide some numerical results to illustrate this analysis in Section 6. Finally, we give short concluding remarks. The main result of this work is Theorem 5.1.

2 Setting

Let $\Omega\subset\mathbb{R}^{2}$ be a domain, that is a connected open set, with Lipschitz boundary $\partial\Omega$ which coincides with the reference strip

[TABLE]

for $|x|\geq d$ where $d>0$ is fixed (see Figure 1). We assume that the propagation of time-harmonic waves in $\Omega$ is governed by the Helmholtz equation with Neumann boundary conditions

[TABLE]

In this problem, $u$ is the quantity of interest (acoustic pressure, velocity potential, component of the electromagnetic field,…), $\Delta$ denotes the 2D Laplace operator, $\lambda$ is a parameter which is proportional to the square of the frequency and $\nu$ stands for the normal unit vector to $\partial\Omega$ directed to the exterior of $\Omega$. Note that from time to time, abusively we will call $\lambda$ the frequency. We emphasize that we consider an academic $\mathrm{2D}$ problem only to simplify the presentation. Other configurations can be dealt with in a completely similar way. In particular, the analysis is the same in higher dimension and in waveguides for which the two unbounded branches are not aligned. Moreover, we can also impose Dirichlet or periodic boundary conditions in (1) to study quantum waveguides or gratings. For $\lambda\in(0;\pi^{2})$, only the plane waves $w_{\pm}$ defined by

[TABLE]

can propagate in $\Omega$. For $\lambda\in(0;\pi^{2})$, the problem (1) has solutions $u_{\pm}$ admitting the decompositions

[TABLE]

Here $R_{\pm}\in\mathbb{C}$ are reflection coefficients and $T\in\mathbb{C}$, which is the same both for $u_{+}$ and $u_{-}$ due to the reciprocity relation, is the transmission coefficient. Moreover, the dots stand for remainders in $\mathrm{H}^{1}(\Omega)$ which decay as $O(e^{-|x|\sqrt{\pi^{2}-\lambda}})$ when $|x|\to+\infty$. Physically, $u_{+}$ (resp. $u_{-}$) models the scattering of the incident rightgoing wave $w_{+}$ (resp. leftgoing wave $w_{-}$) by the perturbation in the geometry with respect to the reference strip $\mathbb{R}\times(0;1)$. We define the scattering matrix

[TABLE]

It is a classical exercise to show that $\mathfrak{s}$ is unitary ($\mathfrak{s}\overline{\mathfrak{s}}^{\top}=\mathrm{Id}$) and symmetric ($\mathfrak{s}=\mathfrak{s}^{\top}$). The functions $u_{\pm}$ are uniquely defined if and only if trapped modes (non-zero solutions of (1) which are in $\mathrm{L}^{2}(\Omega)$) do not exist at the chosen $\lambda$. If trapped modes exist, we define uniquely $u_{\pm}$ as the functions admitting the expansions (3) and which are orthogonal to the linear space of trapped modes (which is of finite dimension) in $\mathrm{L}^{2}(\Omega)$.

We assume that the geometry $\Omega$ is such that $\lambda=\lambda^{0}\in(0;\pi^{2})$ is a simple eigenvalue of the Neumann Laplacian. In other words, we assume there is a non zero $u_{\mathrm{tr}}\in\mathrm{L}^{2}(\Omega)$ satisfying $\Delta u_{\mathrm{tr}}+\lambda^{0}u_{\mathrm{tr}}=0$ in $\Omega$, $\partial_{\nu}u_{\mathrm{tr}}=0$ on $\partial\Omega$ and that any $\mathrm{L}^{2}$ solution of (1) is proportional to $u_{\mathrm{tr}}$. Note that since the continuous spectrum of the Neumann Laplacian in $\Omega$ is $\sigma_{c}=[0;+\infty)$, the eigenvalue is embedded in $\sigma_{c}$. To set ideas, we impose that $\|u_{\mathrm{tr}}\|_{\mathrm{L}^{2}(\Omega)}=1$. Using decomposition in Fourier series, we obtain the expansion

[TABLE]

where $K$ is a constant and $\tilde{u}_{\mathrm{tr}}$ is a remainder which decays as $O(e^{-x\sqrt{4\pi^{2}-\lambda^{0}}})$ when $x\to+\infty$. We assume that $u_{\mathrm{tr}}$ has a slow decay as $x\to+\infty$, i.e. $K\neq 0$. In case $K=0$, the analysis below must be adapted but can be done. Without lost of generality, we can impose that $K>0$. Note that the choice of making an assumption on the decay of $u_{\mathrm{tr}}$ as $x\to+\infty$ is arbitrary. Considering the change $x\mapsto-x$, the analysis below can be developed completely similarly imposing the behaviour as $x\to-\infty$.

3 Perturbation of the frequency and of the geometry

Now, we perturb slightly the original setting supporting trapped modes. First, the spectral parameter $\lambda^{0}$ is changed for

[TABLE]

where $\lambda^{\prime}\in\mathbb{R}$ is given and $\varepsilon>0$ is small. Second, we make a perturbation of amplitude $\varepsilon$ of the geometry to change $\Omega$ into some new waveguide $\Omega^{\varepsilon}$. More precisely, consider $\gamma\subset\partial\Omega$ a smooth arc. In a neighbourhood $\mathscr{V}$ of $\gamma$, we introduce natural curvilinear coordinates $(n,s)$ where $n$ is the oriented distance to $\gamma$ such that $n>0$ outside $\Omega$ and $s$ is the arc length on $\gamma$. Additionally, let $H\in\mathscr{C}^{\infty}_{0}(\gamma)$ be a smooth profile function which vanishes in a neighbourhood of the two endpoints of $\gamma$. Outside $\mathscr{V}$, we assume that $\partial\Omega^{\varepsilon}$ coincides with $\partial\Omega$ and inside $\mathscr{V}$, $\partial\Omega^{\varepsilon}$ is defined by the equation

[TABLE]

In other words if $\gamma$ is parametrized as $\gamma=\{P(s)\in\mathbb{R}^{2}\,|\,s\in I\}$ where $I$ is a given interval of $\mathbb{R}$, then $\gamma^{\varepsilon}:=\{P(s)+\varepsilon H(s)\nu(s)\,|\,s\in I\}$. Here $\nu(s)$ is the unit vector normal to $\gamma$ at point $P(s)$ directed to the exterior of $\Omega$. Finally we consider the perturbed problem

[TABLE]

where $\nu^{\varepsilon}$ stands for the normal unit vector to $\partial\Omega^{\varepsilon}$ directed to the exterior of $\Omega^{\varepsilon}$. We denote by

[TABLE]

the scattering parameters introduced in the previous section in the geometry $\Omega^{\varepsilon}$ at frequency $\lambda$. And for short, we set

[TABLE]

To recall the Theorem 5.1 of [7] describing the behaviour of the scattering matrix $\mathfrak{s}(\varepsilon,\lambda^{0}+\varepsilon\lambda^{\prime})$ as $\varepsilon$ goes to zero, and which will be the basis of our analysis below, we need to introduce a few quantities. Set $U:=(u_{+},u_{-})$ where $u_{\pm}$ are the functions introduced in (3) for $\lambda=\lambda^{0}$. Set also

[TABLE]

Theorem 3.1**.**

$\star$* Assume that $\lambda^{\prime}\neq\kappa(H)$. Then we have*

[TABLE]

$\star$* Assume that $H$ is such that $\kappa(H)\alpha\neq\beta(H)\in\mathbb{C}\times\mathbb{C}$. Then we have*

[TABLE]

with $\tau:=(\kappa(H)\alpha-\beta(H))\,\mathfrak{s}$ and $\tilde{\mu}:=A\mu+B$ for some unimportant real constants $A$, $B$ with $A\neq 0$. We emphasize that $A$, $B$ are independent of $\varepsilon$, $\mu$.

Let us comment this result. The methodology to prove it is the following. First, we compute an asymptotic expansion of an auxiliary object called the augmented scattering matrix, which has been introduced in [37, 22] and [30, 32] as $\varepsilon\to 0$. The essential property is that this augmented scattering matrix considered as a function of $(\varepsilon,\lambda)$ is smooth at $(0,\lambda^{0})$. The procedure and the proof of error estimates are detailed in [31, 32, 33]. Then using the relation existing between the usual scattering matrix and the augmented scattering matrix, we can get the statement of the theorem.

As explained in [7], Theorem 3.1 shows that the scattering matrix $\mathfrak{s}(\cdot,\cdot)$ is not continuous at the point $(0,\lambda^{0})$ (setting where trapped modes exist). Indeed, the function $\mathfrak{s}(\cdot,\cdot)$ valued on different parabolic paths $\{(\varepsilon,\lambda^{0}+\varepsilon\kappa(H)+\varepsilon^{2}\mu),\,\varepsilon\in(0;\varepsilon_{0})\}$ (see Figure 3) have different limits when $\varepsilon$ tends to zero. And for $\varepsilon_{0}\neq 0$ small fixed, the usual scattering matrix $\lambda\mapsto\mathfrak{s}(\varepsilon_{0},\lambda)$ exhibits a quick change in a neighbourhood of $\lambda^{0}+\varepsilon_{0}\kappa(H)$. Indeed, the map $\mu\mapsto\mathfrak{s}(\varepsilon_{0},\lambda^{0}+\varepsilon_{0}\kappa(H)+\varepsilon_{0}^{2}\mu)$ has a large variation for $\mu\in[-C\varepsilon_{0}^{-1};C\varepsilon_{0}^{-1}]$ for some arbitrary $C>0$ (which is only a small change for $\lambda^{\varepsilon_{0}}$). Said differently, a change of order $\varepsilon$ of the frequency leads to a change of order one of the scattering matrix. This is nothing but the the Fano resonance phenomenon. For a given $C>0$, outside an interval of length $C\varepsilon_{0}$ centred at $\lambda^{0}+\varepsilon_{0}\kappa(H)$, $\mathfrak{s}(\varepsilon_{0},\cdot)$ is approximately equal to $\mathfrak{s}^{0}$.

Remark 3.1**.**

When $H$ is such that $\kappa(H)\alpha=\beta(H)\in\mathbb{C}\times\mathbb{C}$, in general a fast Fano resonance phenomenon appears. More precisely, for a given $\varepsilon_{0}\neq 0$ small, the variation of $\mathfrak{s}(\varepsilon_{0},\cdot)$ of order one occurs on a range of frequencies of length $O(\varepsilon_{0}^{2})$ (instead of $O(\varepsilon_{0})$ when $\beta(H)\neq\kappa(H)t$). We write “in general” because we can also show that for well-chosen geometric perturbations, obtained solving a fixed-point problem, no Fano resonance phenomenon happens and the real eigenvalue embedded in the continuous spectrum keeps this property instead of becoming a complex resonance. In particular, this latter result allows one to construct non symmetric waveguides with eigenvalues embedded in the continuous spectrum (see [32, 33]).

From now, we denote by $\tau_{1}$, $\tau_{2}\in\mathbb{C}$ the two components of $\tau$, so that $\tau=(\tau_{1},\tau_{2})$, and we set

[TABLE]

With this notation, the analysis developed in [7] provides the estimate

[TABLE]

where in (12), for any compact set $I\subset\mathbb{R}$, the constant $C>0$ can be chosen independent of $\mu\in I$. In particular, we have

[TABLE]

In order to prove that we have $T^{\varepsilon}(\mu)=0$ for some $\mu\in\mathbb{R}$ for $\varepsilon$ small enough, we first show that the map $\mu\mapsto T^{\mathrm{asy}}(\mu)$ vanishes in $\mathbb{R}$. This is the object of the next section.

4 Asymptotic behaviour of the transmission coefficient

Proposition 4.1**.**

Assume that $T^{0}\neq 0$. Then we have

[TABLE]

where $\mathscr{C}^{\mathrm{asy}}$ is a circle passing through $T^{0}$ and zero.

Proof.

Classical results concerning the Möbius transform guarantee that $\{T^{\mathrm{asy}}(\mu),\,\mu\in\mathbb{R}\}$ coincides with $\mathscr{C}^{\mathrm{asy}}\setminus\{T^{0}\}$ where $\mathscr{C}^{\mathrm{asy}}$ is a circle passing through $T^{0}$. Let us show that $\mathscr{C}^{\mathrm{asy}}$ also passes through zero. From (12), one finds that $T^{\mathrm{asy}}(\mu)=0$ for some $\mu\in\mathbb{R}$ if and only if there holds

[TABLE]

In order to establish (21), we need to derive some relations between $T^{0}$ and $\tau=(\tau_{1},\tau_{2})$. To proceed, first we notice that $U=(u_{+},u_{-})$ satisfies

[TABLE]

Indeed, the first component of $\overline{U}\mathfrak{s}^{0}$ is equal to $R_{+}^{0}\,\overline{u_{+}}+T^{0}\,\overline{u_{-}}$ and using (3), one finds that this function admits the expansion

[TABLE]

From the unitarity of $\mathfrak{s}^{0}$, we infer that $R_{+}^{0}\,\overline{u_{+}}+T^{0}\,\overline{u_{-}}$ has the same expansion as $u_{+}$ at infinity. Using that $u_{\pm}$ are orthogonal to $u_{\mathrm{tr}}$ in $\mathrm{L}^{2}(\Omega)$, we deduce that $R_{+}^{0}\,\overline{u_{+}}+T^{0}\,\overline{u_{-}}=u_{+}$. Similarly, we show that $T^{0}\,\overline{u_{+}}+R_{-}^{0}\,\overline{u_{-}}=u_{-}$, which allows us to conclude to (14). Now we exploit (14) to establish the identity

[TABLE]

From the expressions (8)-(10) of $\kappa(H)$, $\alpha$, $\beta(H)$ and the properties of $\mathfrak{s}^{0}$, we obtain

[TABLE]

Then replacing $U$ by $\overline{U}\mathfrak{s}^{0}$ (identity (14)) in $\kappa(H)\overline{\alpha}-\overline{\beta(H)}$, we get $\overline{\tau}\mathfrak{s}^{0}=(\kappa(H)\alpha-\beta(H))\mathfrak{s}^{0}=\tau$. This proves (15) or equivalently

[TABLE]

Finally, we use (16) to establish (21). The unitarity of $\mathfrak{s}^{0}$ imposes $R_{-}^{0}=-\overline{R_{+}^{0}}T^{0}/\overline{T^{0}}$. Inserting this relation in the second line of (16) gives

[TABLE]

The first line of (16) implies

[TABLE]

Inserting (18) in (17) and multiplying by $\tau_{1}$ leads to

[TABLE]

This is identity (21). ∎

Remark 4.1**.**

The reason why $\mathscr{C}^{\mathrm{asy}}$ passes through zero is quite mysterious. When $\Omega$, $\Omega^{\varepsilon}$ are symmetric with respect to the $(Oy)$ axis, this can be shown quite simply working with half-waveguides problems (see e.g. [7]). But without assumption of symmetry, we cannot provide a physical interpretation of this fact.

Denote $\mu_{\star}$ the value of $\mu$ such that $T^{\mathrm{asy}}(\mu_{\star})=0$ and for $\varepsilon>0$, define the interval $I^{\varepsilon}=(\mu_{\star}-\sqrt{\varepsilon};\mu_{\star}+\sqrt{\varepsilon})$. From (12), for $\varepsilon>0$ small, we know that the curve

[TABLE]

passes close to zero. It remains to show that $C^{\varepsilon}$ passes exactly through zero for $\varepsilon$ small enough.

5 Exact zero transmission

Now, we state and prove the main result of the article. Its proof relies on Proposition 4.1 and an argument presented in [10] (see also [23]).

Theorem 5.1**.**

Assume that $T^{0}\neq 0$. Then there is $\varepsilon_{0}>0$ such that for all $\varepsilon\in(0;\varepsilon_{0}]$, there exists $\mu\in\mathbb{R}$ such that $T^{\varepsilon}(\mu)=0$.

Proof.

Let us first give the general idea of the proof. Assume by contradiction that for all $\varepsilon>0$, $\mu\mapsto T^{\varepsilon}(\mu)$ does not pass through zero in $I^{\varepsilon}$. Since $\mathfrak{s}^{\varepsilon}(\mu)$ is unitary, there holds $R_{+}^{\varepsilon}(\mu)\,\overline{T^{\varepsilon}(\mu)}+T^{\varepsilon}(\mu)\,\overline{R_{-}^{\varepsilon}(\mu)}=0$ and so

[TABLE]

But if $\mu\mapsto T^{\varepsilon}(\mu)$ does not pass through zero on $I^{\varepsilon}$, using Proposition 4.1 one can verify that the point $T^{\varepsilon}(\mu)/\overline{T^{\varepsilon}(\mu)}=e^{2i\mathrm{arg}(T^{\varepsilon}(\mu))}$ must run rapidly on the unit circle for $\mu\in I^{\varepsilon}$ as $\varepsilon\to 0$. On the other hand, $R_{+}^{\varepsilon}(\mu)/\overline{R_{-}^{\varepsilon}(\mu)}$ tends to a constant in $I^{\varepsilon}$ as $\varepsilon\to 0$. This way we obtain a contradiction. We emphasize that the unitary structure of $\mathfrak{s}^{\varepsilon}(\mu)$ is the key ingredient of this step of the proof. Now we make this discussion more rigorous.

Since the circle $\mathscr{C}^{\mathrm{asy}}$ passes through zero, there is $\eta\in(-\pi/2;\pi/2]$ such that $\mathscr{C}^{\mathrm{asy}}$ is tangent to the line $\{\rho\,e^{i\eta}\in\mathbb{C},\ \rho\in\mathbb{R}\}$. Define the quadrants

[TABLE]

see Figure 4. The graph of the map $\mu\mapsto T^{\mathrm{asy}}(\mu)$ crosses both quadrants $Q_{1}$ and $Q_{2}$ in $I_{\varepsilon}$. On the other hand, we have $|T^{\varepsilon}(\mu)-T^{\mathrm{asy}}(\mu)|\leq C\varepsilon$ where $C>0$ is independent of $\mu\in I_{\varepsilon}$ for all $\varepsilon\in(0;\varepsilon_{0}]$. As a consequence, there is $\varepsilon_{0}$ such that for all $\varepsilon\in(0;\varepsilon_{0}]$, the graph of the map $\mu\mapsto T^{\varepsilon}(\mu)$ intersects both $Q_{1}$ and $Q_{2}$ on $I_{\varepsilon}$.

If $\mu\mapsto T^{\varepsilon}(\mu)$ does not vanish in $I^{\varepsilon}$, since $\mu\mapsto T^{\varepsilon}(\mu)$ is continuous, we deduce that for all $\varepsilon\in(0;\varepsilon_{0}]$, there are $a_{\varepsilon},\,b_{\varepsilon}\in I^{\varepsilon}$ such that $T^{\varepsilon}(a_{\varepsilon})=t_{\varepsilon}\,e^{i(\eta-\pi/4)}$ and $T^{\varepsilon}(b_{\varepsilon})=\tilde{t}_{\varepsilon}\,e^{i(\eta+\pi/4)}$, with $t_{\varepsilon}$, $\tilde{t}_{\varepsilon}\in\mathbb{R}\setminus\{0\}$. Taking successively $\mu=a_{\varepsilon}$, $\mu=b_{\varepsilon}$ in the relation preceding (19), we obtain

[TABLE]

Introduce the functions $R_{\pm}^{\mathrm{asy}}$ such that

[TABLE]

From (11), we know that there is $\varepsilon_{0}>0$ such that, for all $\varepsilon\in(0;\varepsilon_{0}]$, we have

[TABLE]

where for $z_{0}\in\mathbb{C}$, $B(z_{0},r)$ denotes the open disk of $\mathbb{C}$ of radius $r>0$ centred at $z_{0}$. From (20), we deduce that we must have both

[TABLE]

This is impossible for $\varepsilon$ small enough because $|R_{-}^{\mathrm{asy}}(\mu_{\star})|=1$ (remember that $T^{\mathrm{asy}}(\mu_{\star})=0$). Thus, we deduce that for all $\varepsilon\in(0;\varepsilon_{0}]$, $\mu\mapsto T^{\varepsilon}(\mu)$ cancels in $I^{\varepsilon}$. ∎

Concerning the zeros of $\mu\mapsto R_{+}^{\varepsilon}(\mu)$, we can make the following comments. When $\varepsilon$ tends to zero, from (11), we know that the curve $\{R_{+}^{\varepsilon}(\mu),\,\mu\in\mathbb{R}\}$ gets closer and closer to $\{R_{+}^{\mathrm{asy}}(\mu),\,\mu\in\mathbb{R}\}$. The set $\{R_{+}^{\mathrm{asy}}(\mu),\,\mu\in\mathbb{R}\}$ is a circle. It passes through zero if and only if we have

[TABLE]

Dividing the first line of (16) by $R_{+}^{0}$ and computing the square of the modulus, we obtain the identity

[TABLE]

Using the above equality, we obtain that (21) is satisfied if and only if there holds

[TABLE]

As a consequence, if $|a|\neq|b|$, for $\varepsilon$ small enough, $\mu\mapsto R_{+}^{\varepsilon}(\mu)$ does not pass through zero. Using the definition of $\tau$ in Theorem 3.1, we observe that we have $|a|=|b|$ if $\Omega$ and $H$ are symmetric with respect to the $(Oy)$ axis. However, surely it is not necessary to consider symmetric geometries to have (22). But we emphasize that if (22) holds in a non symmetric setting, then we cannot work as in the proof of Theorem 5.1 to get exactly $R_{+}^{\varepsilon}(\mu)=0$ for some $\mu\in\mathbb{R}$. Everything lies in the fact that the identity (19) cannot be exploited similarly for the reflection and the transmission coefficients. Therefore, a priori nothing guarantees that exact zero reflection occurs during the Fano resonance phenomenon in a non symmetric waveguide, even when (22) is satisfied.

6 Numerical results

In this section, we illustrate the results obtained above. In the first series of experiments, we work in the geometry

[TABLE]

pictured in Figure 5 left. In $\Omega:=\Omega^{0}$ the obstacle is symmetric with respect to the line $\mathbb{R}\times\{1/2\}$. According to the results of the literature (see e.g. [15]), we know that there are trapped modes for certain real frequencies in this geometry. Using Perfectly Matched Layers [4, 3, 5], we find that they exist for $\sqrt{\lambda^{0}}\approx 1.9939$. Figure 5 right represents such a trapped mode in $\Omega$.

The domain $\Omega^{\varepsilon}$ is obtained from $\Omega$ by shifting by $\varepsilon$ the obstacle along the $(Oy)$ axis. Admittedly, this kind of perturbation is not exactly the one considered in (6). However, since there exists an almost identical mapping from $\Omega$ to $\Omega^{\varepsilon}$, results are similar. We emphasize that for $\varepsilon>0$, $\Omega^{\varepsilon}$ has no symmetry property. In Figure 6, we display the values of the complex scattering coefficients $R_{+}(\varepsilon,\lambda)$, $T(\varepsilon,\lambda)$ appearing in the decomposition (3) of $u_{+}$ for $\varepsilon=0.05$ and for $\sqrt{\lambda}\in(1.97;2.03)$ (note that this interval contains the value $\sqrt{\lambda^{0}}$). To proceed, we use a $\mathrm{P2}$ finite element method in a truncated geometry. On the artificial boundary created by the truncation, a Dirichlet-to-Neumann operator with $20$ terms serves as a transparent condition. As expected, we observe that $\lambda\mapsto T(\varepsilon,\lambda)$ passes through zero.

In Figure 7, we display the curves $\lambda\mapsto|T(\varepsilon,\lambda)|$ for several $\varepsilon$ and a range of values of $\lambda$. The right picture is a zoom of the left picture around $\lambda^{0}$. As expected we observe that for the different $\varepsilon$, we have $T(\varepsilon,\lambda)=0$ for one $\lambda$ close to $\lambda^{0}$. We also note that the smaller $\varepsilon>0$ is, the faster the Fano resonance phenomenon occurs. This is also expected. Finally, in Figure 8, we display the real part of $u_{+}$ (see (3)) in $\Omega^{\varepsilon}$ for $\varepsilon=0.05$ and $\sqrt{\lambda}=2.0072$. In this setting, there holds $T(\varepsilon,\lambda)\approx 0$. And indeed we observe that the incident rightgoing wave $w_{+}$ is completely backscattered, this is the mirror effect.

In the second series of experiments, we work in the geometry of Figure 10. Using Perfectly Matched Layers, we find a complex resonance $\lambda_{c}$ such that $\sqrt{\lambda_{c}}\approx 2.49-0.15i$. In Figure 9, we display the values of the complex scattering coefficients $R_{+}(\lambda)$, $T(\lambda)$ appearing in the decomposition (3) of $u_{+}$ for $\sqrt{\lambda}\in(2.1;2.8)$ (note that this interval contains the value $\Re e\,\sqrt{\lambda_{c}}$). Though this experiment does not strictly enter the framework presented in this note (we do not start from a situation where trapped modes exist), we observe that the curve $\lambda\mapsto T(\lambda)$ passes through zero for $\lambda$ in a neighbourhood of $\Re e\,\lambda_{c}$. In Figure 10, we display the real part of $u_{+}$ for $\sqrt{\lambda}=2.4016$. In this setting, we have $T(\lambda)\approx 0$.

7 Concluding remarks

In this note, we proved that during the Fano resonance phenomenon in monomode regime, without assumption of symmetry of the geometry, the transmission coefficient passes through zero. Physically, when the transmission coefficient is null, the energy of an incident wave propagating through the structure is completely backscattered. As already mentioned, everything presented here is also valid in higher dimension and with Dirichlet or periodic boundary conditions instead of Neumann ones. We considered a geometrical perturbation of the walls of the waveguide. We could also have worked with a penetrable inclusion placed in the waveguide. Then perturbing the material parameter, we would have obtained similar results. Importantly, the above analysis applies only in monomode regime, that is for our geometry when $\lambda^{0}$ belongs to $(0;\pi^{2})$. It is not clear what happens in multimodal regime ($\lambda^{0}>\pi^{2})$. Moreover, we assumed that $\lambda^{0}$ is a simple eigenvalue embedded in the continuous spectrum of the Neumann Laplacian. When $\lambda^{0}$ is not simple, the analysis has to be done.

Acknowledgments

The research of S.A. Nazarov was supported by the grant No. 17-11-01003 of the Russian Science Foundation.