Examples of non-scattering inhomogeneities

Lucas Chesnel¹, Houssem Haddar¹, Hongjie Li², Jingni Xiao^3
1 IDEFIX, Ensta Paris, Institut Polytechnique de Paris, Palaiseau, France;
² Yau Mathematical Sciences Center, Tsinghua University, Beijing, China;
³ Department of Mathematics, Drexel University, Philadelphia, Pennsylvania, USA.
E-mails: lucas.chesnel@inria.fr, houssem.haddar@inria.fr, hongjieli@tsinghua.edu.cn, jingni.xiao@drexel.edu.
(June 25, 2024)

Abstract. We consider the scattering of waves by a penetrable inclusion embedded in some reference medium. We exhibit examples of materials and geometries for which non-scattering frequencies exist, i.e., for which at some frequencies there are incident fields which produce null scattered fields outside of the inhomogeneity. We show in particular that certain domains with corners or even cusps can support non-scattering frequencies. We relate the latter, for some inclusions, to resonance frequencies for Dirichlet or Neumann cavities. We also find situations where incident non-scattering fields solve the Helmholtz equation in a neighborhood of the inhomogeneity and not in the whole space. Finally, in relation with invisibility, we give examples of inclusions of anisotropic materials which are non-scattering for all real frequencies. We prove that corresponding material indices must have a special structure on the boundary.
Key words. Non-scattering frequencies, transmission eigenvalues, interior transmission problem.
Mathematics Subject Classification. 35J05, 35P25, 35N25

1 Introduction

The so-called transmission eigenvalues play an important role in the study of inverse scattering from inhomogeneous media. They can be helpful in addressing theoretical questions such as uniqueness of the perturbation or in the justification of some reconstruction methods such as the Linear Sampling Method [12, 6]. These special frequencies can also be exploited in imaging algorithms for highly cluttered media [1, 2]. In general, transmission eigenvalues are associated with eigenfunctions that cannot be exactly represented as solutions to the Helmholtz equation in the whole space. In that case, exact non-scattering does not occur. Proving that incident fields always scatter in certain configurations has been the subject of several works and is usually associated with the presence of some geometric singularities such as corners in the boundary of the inhomogeneity. The techniques and results differ according to the considered model for the propagation of waves inside the inclusion, namely

$$\Delta u+k^{2}qu=0$$ (1)

or

$$\mathrm{div}(A\nabla u)+k^{2}qu=0.$$ (2)

Here $k$ denotes the frequency and $A$, $q$ are functions characterizing the physical properties of the materials. We refer the reader respectively to [3, 25, 16, 9, 28, 19] and [4, 11, 29, 10] for studies concerning cases (1) and (2). In particular, it is proven for model (1) that if there is a corner in the boundary of the support of $q-1$, then at any frequency an incident wave will always produce a non zero scattered field outside the inhomogeneity, provided that $q$ is sufficiently regular and different from the background coefficient $1$ near the corner. However, things can be different for model (2) especially when $A$ is different from the background coefficient $\mathrm{Id}$ near the boundary of the support of the inhomogeneity. The main goal of the present article is to provide examples for which non-scattering occur for model (2). We show that this can happen even when the domain has singularities like corners or cusps. This analysis can be helpful in better understanding the optimality of some results concerning the absence of non-scattering frequencies in the literature, and provide some insight on the difference between cases (1) and (2). We also give examples of situations where non-scattering incident fields have some singularities outside the inhomogeneity. This motivates the definition of non-scattering frequencies adopted below.

Let us first describe the scattering problem we are considering. It is associated with the scalar acoustic wave equation in $\mathbb{R}^{d}$, $d\geq 2$. We assume that some inclusion is located in a bounded domain $\Omega$ with Lipschitz boundary $\partial\Omega$ such that $\mathbb{R}^{d}\setminus\overline{\Omega}$ is connected. In accordance with (2), the inclusion is characterized by some material properties $A\in\mathrm{L}^{\infty}(\Omega,\mathbb{R}^{d\times d})$, $q\in\mathrm{L}^{\infty}(\Omega,\mathbb{R})$ such that

$$A=\mathrm{Id},\ q=1\mbox{ in }\mathbb{R}^{d}\setminus\overline{\Omega};\qquad\underset{\mathrm{x}\in\Omega}{\mbox{ess inf }}\underset{\xi\in\mathbb{C}^{d},|\xi|=1}{\mbox{inf }}(\xi\cdot A(\mathrm{x})\overline{\xi})>0;\qquad\underset{\mathrm{x}\in\Omega}{\mbox{ess inf }}q(\mathrm{x})>0.$$ (3)

Below, we refer without distinction to $(A,q;\Omega)$ as the inclusion or the inhomogeneity. Consider some incident field $u_{i}$ satisfying

$$\Delta u_{i}+k^{2}u_{i}=0\quad\mbox{ in }\tilde{\Omega},$$ (4)

where $\tilde{\Omega}$ is a neighborhood of $\Omega$, more precisely a domain such that $\overline{\Omega}\subset\tilde{\Omega}$. Note that we do not assume $u_{i}$ be defined and smooth in the whole $\mathbb{R}^{d}$. This allows us to take into account the possibility of illuminating the inclusion for example by incident fields due to point sources. The scattering of $u_{i}$ due to the inhomogeneity is governed by the following time-harmonic problem where $u_{s}$ denotes the scattered field,

$$\begin{array}[]{|rcll}\mathrm{div}(A\nabla u_{s})+k^{2}qu_{s}&=&-\mathrm{div}((A-\mathrm{Id})\nabla u_{i})-k^{2}(q-1)u_{i}&\mbox{ in }\mathbb{R}^{d}\\[5.0pt] \vrule\lx@intercol\hfil\displaystyle{\displaystyle\lim_{R\to+\infty}\int_{|x|=R}\left|\frac{\partial u_{s}}{\partial r}-iku_{s}\right|^{2}ds(\mathrm{x})=0.}\hfil\lx@intercol\end{array}$$ (5)

The derivatives in the first line of (5) are to be understood in the weak sense where $u_{i}$ is extended by zero outside $\tilde{\Omega}$. The last line of (5), in which $r=|\mathrm{x}|$, is the so-called Sommerfeld radiation condition. The function $u:=u_{i}+u_{s}$ is usually referred to as the total field. For all $k>0$, Problem (5) has a unique solution $u_{s}$ belonging to $\mathrm{H}^{1}(\mathcal{O})$ for all bounded domains $\mathcal{O}\subset\mathbb{R}^{d}$ if $d=2$. When $d=3$, the existence of solutions can be established under the same assumptions. The uniqueness for $d=3$ has been proven with the additional condition that $A$ is Lipschitz continuous in $\overline{\Omega}$, by using the unique continuation principle, see e.g. [6].
We define in the following the so-called non-scattering frequencies.

Definition 1.1.

We say that $k>0$ is a non-scattering frequency if there is $\tilde{\Omega}$ such that $\overline{\Omega}\subset\tilde{\Omega}$ and $u_{i}$ solving (4) for which the corresponding scattered field satisfies $u_{s}=0$ in $\mathbb{R}^{d}\setminus\Omega$. Such $u_{i}$ is referred to as a non-scattering incident field. We denote by $\mathscr{S}_{\mathrm{NS}}$ the set of non-scattering frequencies of (5).

If $u_{i}$ is a non-scattering incident field, then by setting $v=u_{i}$, we see that $(u,v)\in\mathrm{H}^{1}(\Omega)\times\mathrm{H}^{1}(\Omega)$ solves the interior transmission eigenvalue problem

$$\begin{array}[]{|rclll}\mathrm{div}(A\nabla u)+k^{2}qu&=&0&\mbox{ in }\Omega\\% [3.0pt] \Delta v+k^{2}v&=&0&\mbox{ in }\Omega\\[3.0pt] u-v&=&0&\mbox{ on }\partial\Omega\\[3.0pt] \nu\cdot A\nabla u-\nu\cdot\nabla v&=&0&\mbox{ on }\partial\Omega.\end{array}$$ (6)

Here $\nu$ stands for the unit outward normal vector to $\partial\Omega$.

Definition 1.2.

We say that $k>0$ is a transmission eigenvalue if there is a non-zero $(u,v)$ which satisfies (6). We denote by $\mathscr{S}_{\mathrm{TE}}$ the set of transmission eigenvalues.

The above discussion shows that $\mathscr{S}_{\mathrm{NS}}\subset\mathscr{S}_{\mathrm{TE}}$. However the converse does not hold in general. More precisely, we have the following characterization result whose proof is straightforward.

Proposition 1.3.

A real frequency $k\in\mathscr{S}_{\mathrm{TE}}$ is also in $\mathscr{S}_{\mathrm{NS}}$ if and only if there is a non-zero eigenpair $(u,v)$ of (6) associated with $k$ such that $v$ can be extended in $\tilde{\Omega}$, a neighborhood of $\Omega$, as a function $\tilde{v}$ satisfying

$$\Delta\tilde{v}+k^{2}\tilde{v}=0\mbox{ in }\tilde{\Omega}.$$ (7)

Concerning (6), it has first been shown that $\mathscr{S}_{\mathrm{TE}}$ is discrete for rather mild conditions on the coefficients $A$ and $n$ [6]. Under a bit more restrictive hypothesis on $A$ and $q$, people have been able to prove that $\mathscr{S}_{\mathrm{TE}}$ contains a countably infinite sequence of eigenvalues accumulating only at $+\infty$ (see [26, 17, 8, 7, 27]). In contrast, the set $\mathscr{S}_{\mathrm{NS}}$ has been less studied. Certain simple examples of geometries (rectangle, balls) where $\mathscr{S}_{\mathrm{NS}}\neq\emptyset$ have been given in [21, 10]. On the other hand, it has been established that for certain scatterers with non-smooth geometries (having corners in 2D or conical tips/edges in 3D), there holds $\mathscr{S}_{\mathrm{NS}}=\emptyset$ (see [3, 25, 16]).

As indicated above, the main goal of the present article is to provide examples of configurations where $\mathscr{S}_{\mathrm{NS}}\neq\emptyset$. We first consider the case where the material properties have the special form $A=a\,\mathrm{Id}$, $q=a$, with $a>0$ (Section 2). This configuration is interesting for several reasons. For instance, one can completely characterize in this case the set $\mathscr{S}_{\mathrm{TE}}$ as the union of Dirichlet and Neumann cavity eigenvalues. One can therefore exhibit non-scattering frequencies by considering either Dirichlet or Neumann eigenvalues. We show in particular that some polygonal domains (convex or not) possess non-scattering frequencies. We also exhibit examples, some inspired by and some extracted from [15], where non-scattering incident fields are not entire solutions to the Helmholtz equation (i.e. $\tilde{\Omega}\neq\mathbb{R}^{d}$). We consider separately Dirichlet and Neumann eigenvalues. While examples for the first case are already present in the literature, the examples for Neumann eigenvalues are harder to find analytically and seem less known. We propose a method to construct such domains based on analysing the characteristics for the gradient of the solutions to the Helmholtz equation. This allows for instance to prove existence of non-scattering frequencies for certain domains with corners of arbitrary angles and even cusps. We emphasize that on the contrary, non-scattering domains for the Dirichlet eigenvalues can only have corners of particular apertures. In Section 3, we work with anisotropic materials and exhibit configurations where non-scattering occurs at any frequencies. The first examples are inspired by the concept of invisibility by diffeomorphisms [18]. For associated anisotropies we prove in particular that the matrix $A$ has a special structure at irregular points of the boundary as it must satisfy $A\nu=\nu$ (see also [10]). We also provide other examples of anisotropies (not associated with diffeomorphism transformations) for which non-scattering occur at any frequency in the case of regular and non-regular domains.

Note: we say that a function is an entire solution of the Helmholtz equation if it solves the homogeneous Helmholtz equation in $\mathbb{R}^{d}$.

2 Non-scattering in the case $A=a\,\mathrm{Id}$, $q=a$, with $a>0$

In this section, we make the assumption that there holds

$$A=a\,\mathrm{Id}\quad\text{and}\quad q=a\qquad\mbox{ in }\Omega,$$ (8)

where $a>0$ is a constant such that $a\neq 1$. In this particular situation, we first show that $\mathscr{S}_{\mathrm{TE}}$ is formed by the union of the eigenvalues of the Dirichlet and Neumann Laplacians in $\Omega$. More precisely, introduce the problems

$$\begin{array}[]{|rl}\Delta w_{D}+k^{2}w_{D}=0&\mbox{ in }\Omega\\[3.0pt] w_{D}=0&\mbox{ on }\partial\Omega\end{array}\qquad\mbox{ and }\qquad\begin{array}[]{|rl}\Delta w_{N}+k^{2}w_{N}=0&\mbox{ in }\Omega\\[3.0pt] \partial_{\nu}w_{N}=0&\mbox{ on }\partial\Omega.\end{array}$$ (9)

We denote by $\mathscr{S}_{\mathrm{D}}$ (resp. $\mathscr{S}_{\mathrm{N}}$) the values of $k>0$ such that the Dirichlet (resp. Neumann) problem (9) admits non-zero solutions in $\mathrm{H}^{1}(\Omega)$. We have the following result.

Lemma 2.1.

When $A$, $q$ satisfy (8), we have $\mathscr{S}_{\mathrm{TE}}=\mathscr{S}_{\mathrm{D}}\cup\mathscr{S}_{\mathrm{N}}$.

Proof.

Suppose that $(u,v)$ is a transmission eigenpair solving (6). Define $w_{D}:=u-v$ and $w_{N}:=au-v$. Then one readily checks that $w_{D}$, $w_{N}$ solve respectively the Dirichlet and Neumann eigenvalue problems (9). Moreover, at least one of the functions $w_{D}$ and $w_{N}$ must be non-trivial. Hence, we have $\mathscr{S}_{\mathrm{TE}}\subset(\mathscr{S}_{\mathrm{D}}\cup\mathscr{S}_{\mathrm{N}})$.
Conversely, assume that $w_{D}$ is a non-zero solution of the Dirichlet problem appearing in (9). Then $(u,v)=(aw_{D},w_{D})$ is a (non-trivial) eigenpair of (6). Similarly, if $w_{N}$ is a non-zero solution of the Neumann problem appearing in (9), then $(u,v)=(w_{N},w_{N})$ solves (6). This guarantees that $(\mathscr{S}_{\mathrm{D}}\cup\mathscr{S}_{\mathrm{N}})\subset\mathscr{S}_{\mathrm{TE}}$. ∎

Remark 2.2.

Note that the case where $A$, $q$ satisfy (8) is a situation where all the $k^{2}\in\mathbb{C}$ such that (6) admits a non-zero solution are real. It is not known if this can happen with other pairs $A$, $q$.

From Lemma 2.1, when $A$, $q$ satisfy (8), the question of finding non-scattering eigenvalues can be reformulated as “are there eigenfunctions of the Dirichlet or Neumann Laplacian in $\Omega$ that can be extended as solutions to the Helmholtz equation in a neighborhood of $\overline{\Omega}$?”. Some positive answers can be given in three situations:
- For certain geometries $\Omega$, one can compute analytically the eigenfunctions of the Dirichlet/Neumann Laplacians and observe that they are defined in larger domains than $\Omega$;
- For certain $\Omega$, one can use reflections to extend the eigenfunctions of the Dirichlet/Neumann Laplacians in $\Omega$ to larger domains;
- Given a function $u_{i}$ solving the Helmholtz equation (4) in some given domain $\tilde{\Omega}$, one can look for $\Omega$, with $\overline{\Omega}\subset\tilde{\Omega}$, for which $u_{i}$ is an eigenfunction of the Dirichlet or Neumann Laplacian in $\Omega$.
We present corresponding results in the following subsections.

2.1 Domains where analytic expressions can be obtained for the eigenfunctions

Assume here that $\Omega$ coincides with the unit square and consider $k\in\mathscr{S}_{\mathrm{D}}$ as well as $w_{D}$ a corresponding eigenfunction of the Dirichlet problem appearing in (9). In an appropriate set of coordinates, $w_{D}$ writes as a linear combination of the functions

$$\sin(m\pi x)\sin(n\pi y),\qquad m,n\in\mathbb{N}^{\ast}:=\{1,2,3,\dots\}.$$

In the proof of Lemma 2.1, we have seen that $(u,v)=(aw_{D},w_{D})$ constitutes an eigenpair of (6). Clearly $v$ extends as a function solving the homogeneous Helmholtz equation (7) in $\mathbb{R}^{2}$. From Proposition 1.3, we infer that $k\in\mathscr{S}_{\mathrm{NS}}$. Interestingly, as explained in [3, 10], in that case one can exhibit non-scattering incident fields $u_{i}$ which are simple combinations of plane waves. More precisely, $u_{i}$ such that

$$\begin{array}[]{rcl}u_{i}(x,y)&=&4\sin(m\pi x)\sin(n\pi y)\\[3.0pt] &=&e^{i\pi(mx-ny)}+e^{-i\pi(mx-ny)}-e^{i\pi(mx+ny)}-e^{-i\pi(mx+ny)}\end{array}$$

produces a scattered field which is exactly zero outside of $\Omega$. Similarly we establish that $\mathscr{S}_{\mathrm{N}}\subset\mathscr{S}_{\mathrm{NS}}$.

This reasoning can be adapted to deal with other simple domains $\Omega$ where we can use separation of variables. This allows one to state the following result.

Proposition 2.3.

Assume that $\Omega\subset\mathbb{R}^{2}$ is a rectangle, a disk, or an elliptical domain (whose boundary is an ellipse). Then when $A$, $q$ satisfy (8), we have $\mathscr{S}_{\mathrm{NS}}=\mathscr{S}_{\mathrm{TE}}=\mathscr{S}_{\mathrm{D}}\cup\mathscr{S}_{\mathrm{N}}$.

Note that to consider the case of elliptic domains, we work with the elliptic coordinates $(\mu,\phi)$ such that, after a rigid change of variables,

$$x=\alpha\cosh\mu\cos\phi,\qquad y=\alpha\sinh\mu\sin\phi,$$

where $\mu\geq 0$, $\phi\in\mathbb{R}/(2\pi)$ and $(\pm\alpha,0)$ stand for the foci of the ellipse. In particular, the curves $\mu=\mathrm{constant}>0$ are confocal ellipses. Solutions of (9) can be decomposed on functions with separate variables in $\mu$, $\phi$. For the latter functions, we find that the dependences in $\mu$, $\phi$ satisfy respectively a modified Mathieu’s equation and a Mathieu’s equation with periodic boundary conditions, which are both Sturm-Liouville problems. Moreover it is known that the solutions of the Mathieu’s equation which are regular at zero can be extended in $\mathbb{R}$ so that all eigenfunctions of the Dirichlet and Neumann Laplacians in a bounded elliptical domain can be extended as solutions of an Helmholtz equation in $\mathbb{R}^{2}$. For more details, we refer the readers to [14, 23].
Similar to Proposition 2.3, we have the following statement in $\mathbb{R}^{3}$.

Proposition 2.4.

Assume that $\Omega\subset\mathbb{R}^{3}$ is a ball, an ellipsoid, or, in an appropriate system of coordinates, there holds $\Omega=I\times\omega$ where $I$ is a bounded open interval and $\omega\subset\mathbb{R}^{2}$ is a rectangle, a disk, or an elliptical domain. Then when $A$, $q$ satisfy (8), we have $\mathscr{S}_{\mathrm{NS}}=\mathscr{S}_{\mathrm{TE}}=\mathscr{S}_{\mathrm{D}}\cup\mathscr{S}_{\mathrm{N}}$.

Figure 1: Sector (left), spherical wedge (center) and cylindrical wedge (right).

We now turn our attention to other geometries. Consider some constants $\alpha\in(0,2\pi)$ and $\ell>0$. Define in $\mathbb{R}^{2}$ the sector

$$\Omega=\{(r\cos\phi,r\sin\phi)\,|\,r\in(0,\ell),\,\phi\in(0,\alpha)\}$$

and in $\mathbb{R}^{3}$, the spherical wedge

$$\Omega=\{(r\cos\phi\sin\theta,r\sin\phi\sin\theta,r\cos\theta)\,|\,r\in(0,\ell),\theta\in(0,\pi),\phi\in(0,\alpha)\}$$

as well as the cylindrical wedge

$$\Omega=\{(r\cos\phi,r\sin\phi,z)\,|\,r\in(0,\ell),\,\phi\in(0,\alpha),z\in I\}$$

(see Figure 1 for illustrations). Here $I$ is an open bounded interval. For theses domains, depending on the opening angle, there may exist zero or infinitely many non-scattering frequencies.

Proposition 2.5.

With the notation above, assume that $\Omega$ is either a sector in $\mathbb{R}^{2}$, a spherical wedge in $\mathbb{R}^{3}$ or a cylindrical wedge in $\mathbb{R}^{3}$, of angle $\alpha$. When $A$, $q$ satisfy (8), we have:
- if $\alpha\in(0,2\pi)\cap\mathbb{Q}\pi$, then $\mathscr{S}_{\mathrm{NS}}$ contains an unbounded sequence which accumulates only at $+\infty$;
- if $\alpha\in(0,2\pi)\cap(\mathbb{R}\pi\setminus\mathbb{Q}\pi)$, then $\mathscr{S}_{\mathrm{NS}}=\emptyset$.

Remark 2.6.

In fact, we show in the proof that $\mathscr{S}_{\mathrm{NS}}$ coincides exactly with $\mathscr{S}_{\mathrm{TE}}=\mathscr{S}_{\mathrm{D}}\cup\mathscr{S}_{\mathrm{N}}$ if and only if $\alpha\in(0,2\pi)\cap(\pi/\mathbb{N}^{\ast})$.

Remark 2.7.

In the second item of Proposition 2.5, the fact that $\mathscr{S}_{\mathrm{NS}}=\emptyset$ in that domains $\Omega$ is not only due to the value of the opening angle $\alpha$ but also to the shape of the whole $\Omega$. Indeed in §2.3.2, we will exhibit examples of geometries in the plane supporting non-scattering frequencies with corners of arbitrary $\alpha\in(0,2\pi)$.

Proof.

Let $\Omega$ be a sector in $\mathbb{R}^{2}$. The eigenfunctions of the Dirichlet Laplacian in $\Omega$ coincide with the functions $J_{m\pi/\alpha}(\mu_{m}r)\sin(m\pi\phi/\alpha)$, with $m\in\mathbb{N}^{\ast}$ and $\mu_{m}>0$ such that $J_{m\pi/\alpha}(\mu_{m}\ell)=0$. Here $J_{m\pi/\alpha}$ is the Bessel function of the first kind with order $m\pi/\alpha$. As a consequence, by looking at the singularity of Bessel functions at the origin and periodicity of the sin function, we infer that an eigenfunction can be extended as an entire solution of the Helmholtz equation if and only if $\alpha\in(0,2\pi)\cap\pi/\mathbb{N}^{\ast}$. The analysis is similar for eigenfunctions of the Neumann Laplacian.
The cases of spherical and cylindrical wedges in $\mathbb{R}^{3}$ can be dealt similarly by working, respectively, with spherical and cylindrical coordinates. ∎

With this approach, we could also consider 3D conical tip, i.e. domains of the form

$$\Omega=\{(r\cos\theta\cos\phi,r\cos\theta\sin\phi,r\sin\theta)\,|\,0<r<\ell,0<\theta<\theta_{0},0\leq\phi<2\pi\}$$

with $\theta_{0}>0$. In this case, the possibility of extending eigenfunctions of the Dirichlet/Neumann Laplacians in $\Omega$ depends on whether $\cos\theta_{0}$ is a zero of the Legendre polynomial $P^{m}_{n}$ for some integers $n$ and $m\in[-n,n]$. We choose not to elaborate much on this direction.

Let us mention that some of the discussions above for 2D cases can also be found in [20].

2.2 Domains where eigenfunctions can be extended by reflections

Now we present other geometries where some or all the Dirichlet/Neumann eigenfunctions of (9) can be extended, this time by working with reflections. We start with a definition.

Definition 2.8.

Let $\Omega$ be a polygon of $\mathbb{R}^{2}$ (resp. a polyhedron of $\mathbb{R}^{3}$). We say that $\Omega$ is a proper paving unit if the following two conditions hold:
- One can find a combination of successive reflections of $\overline{\Omega}$ with respect to the edges (resp. faces) that allows one to cover a domain $\tilde{\Omega}\supset\overline{\Omega}$ with the overlaps only on the skeleton $\mathcal{S}$ consisting of the edges (resp. faces) of $\partial\Omega$ and their images by the reflections.
- If one assigns a different colour to each of the edges (resp. faces) of $\partial\Omega$, then each element of $\mathcal{S}$ may inherit two or more colours due to reflections. We require that every element of $\mathcal{S}$ has only one colour under reflections.

In Figure 2, we present an example of domain $\Omega$ and covering which does not satisfy the second item of Definition 2.8.

(a) Domain $\Omega$. The vertices of the triangle are such that $\angle A=2\pi/5$ and $\angle B=\angle C=3\pi/10$.

(b) By reflecting $\overline{\Omega}$, we can fill the above pentagon.

(c) The second item of Definition 2.8 is not satisfied because we have two different colours on $[AB]$.

Figure 2: A domain $\Omega$ that is not a proper paving unit.

Proposition 2.9.

Assume that $\Omega\subset\mathbb{R}^{d}$, $d=2,3$, is a proper paving unit. Then when $A$, $q$ satisfy (8), we have $\mathscr{S}_{\mathrm{NS}}=\mathscr{S}_{\mathrm{TE}}=\mathscr{S}_{\mathrm{D}}\cup\mathscr{S}_{\mathrm{N}}$.

Proof.

Let $\Omega$ be a proper paving unit. Consider one particular edge or face of $\partial\Omega$ and introduce some system of coordinates $x^{\prime},x_{d}$ with $x^{\prime}\in\mathbb{R}^{d-1}$, such that this edge or face lies in the hypersurface $\{x_{d}=0\}$. Let $\Omega_{1}$ denote the reflection of $\Omega$ with respect to the line or plane $\{x_{d}=0\}$. Consider $v$ an eigenfunction of the Dirichlet Laplacian in $\Omega$. Classically we can extend it to $\Omega_{1}$ by defining $v(x^{\prime},x_{d})=-v(x^{\prime},-x_{d})$ for all $(x^{\prime},x_{d})\in\Omega_{1}$. Similarly, Neumann eigenfunctions can be extended in $\Omega_{1}$ by defining $v(x^{\prime},x_{d})=v(x^{\prime},-x_{d})$. It can be verified straightforwardly that such extended function $v$ satisfies $\Delta v+k^{2}v=0$ in the whole interior of $\overline{\Omega\cup\Omega_{1}}$. Repeating this reflection argument successively, we can extend Dirichlet/Neumann eigenfunctions to $\tilde{\Omega}$ as solutions to the homogeneous Helmholtz equation. The proof is completed. ∎

From classical results of interior regularity, this shows that in proper paving unit, Dirichlet and Neumann eigenfunctions are smooth up to the boundary because their extensions are real-analytic in $\tilde{\Omega}$.

One can check that triangles that are equilateral (60°-60°-60°), hemiequilateral (30°-60°-90°) or isosceles right (45°-45°-90°) are all proper paving unit of $\mathbb{R}^{2}$. In fact, for these three types of triangles the Dirichlet/Neumann eigenfunctions can be expressed with trigonometric functions as discovered by G. Lamé (see [22]). Besides, observe that if $\omega$ is a proper paving unit of $\mathbb{R}^{2}$ then $\Omega=I\times\omega$ is a proper paving unit of $\mathbb{R}^{3}$ for any bounded interval $I$.

Finally, let us mention that for domains $\omega$, with $\overline{\omega}\subset\tilde{\Omega}$, obtained by “properly” reflecting a given proper paving unit $\Omega$ a finite number of times (see an illustration in 2D with Figure 3), $\mathscr{S}_{\mathrm{NS}}$ contains an unbounded sequence. Indeed, any eigenfunction of the Dirichlet/Neumann Laplacian in $\Omega$ extended to $\tilde{\Omega}$ according to the process described above is clearly an eigenfunction of the Dirichlet/Neumann Laplacian in $\omega$.

Figure 3: Examples of domains for which $\mathscr{S}_{\mathrm{NS}}$ contains an unbounded sequence.

2.3 Combining solutions to the homogeneous Helmholtz equation

In the previous subsections, we presented canonical examples of geometries $\Omega$ where eigenfunctions of the Dirichlet/Neumann Laplacians can be extended as smooth solutions of the homogeneous Helmholtz equation in a neighborhood of $\Omega$. Here, by working directly with solutions to the homogeneous Helmholtz equation, we provide examples of more general domains such that non-scattering frequencies exist. To keep notations simple, we work in $\mathbb{R}^{2}$ only. The ideas can be applied to generate examples in $\mathbb{R}^{3}$.

2.3.1 Extendable Dirichlet eigenfunctions

Consider the function $v_{0}$ such that

$$v_{0}(x,y)=\sin x\sin y.$$

It is a particular solution of the equation

$$\Delta v+2v=0\quad\mbox{ in }\mathbb{R}^{2}.$$ (10)

Now let us, roughly speaking, rotate this function by an angle $\alpha\in(0,2\pi)$ by defining $v_{\alpha}$ such that

$$v_{\alpha}(x,y)=v_{0}(\cos(\alpha)x-\sin(\alpha)y,\sin(\alpha)x+\cos(\alpha)y).$$

Note that $v_{\alpha}$ also solves (10). Then for $\lambda\in\mathbb{R}$, set

$$v^{\lambda}_{\alpha}=v_{0}+\lambda\,v_{\alpha}$$

and define the nodal set $\mathscr{N}^{\lambda}_{\alpha}:=\{(x,y)\in\mathbb{R}^{2}\,|\,v^{\lambda}_{\alpha}(x,y)=0\}$. As classical in literature (see e.g. [13]), we call nodal domains of $v^{\lambda}_{\alpha}$ the maximally connected subsets of $\mathbb{R}^{2}$ for which $v^{\lambda}_{\alpha}$ does not change sign. From Proposition 1.3 and Lemma 2.1, we obtain the following statement.

Proposition 2.10.

Any bounded nodal domain of $v^{\lambda}_{\alpha}$ is a domain where $2\in\mathscr{S}_{\mathrm{NS}}$.

Let us make a few comments concerning this approach. First the choice $k^{2}=2$ is arbitrary here and we could consider any other $k^{2}>0$. On the other hand, we combined only two particular solutions of (10). We could have worked similarly with any other linear combinations of functions satisfying (10). To exhibit different solutions of (10), we can proceed to rotations of $v_{0}$ as we did. We can also translate $v_{0}$ or its rotated versions. This a priori offers a large variety of functions. A natural question then is “how rich is the family of corresponding bounded nodal domains?”. Finally, observe that we prove only the existence of one element in $\mathscr{S}_{\mathrm{NS}}$ and not that $\mathscr{S}_{\mathrm{NS}}$ contains an unbounded sequence as in the statements of the previous subsections.

In Figure 4, we display the nodal sets $\mathscr{N}^{\lambda}_{\alpha}$ for $\alpha=\pi/4$ and $\lambda\in\{-1,-0.2,0.5,2\}$. In each situation we observe that there exist bounded nodal domains. The question of proving the existence of bounded nodal domains would deserve to be studied in more details.

Figure 4: Nodal sets $\mathscr{N}^{\lambda}_{\alpha}$ for $\alpha=\pi/4$ and, from left to right, $\lambda=-1,-0.2,0.5,2$.

Interestingly, this technique can be exploited to exhibit domains $\Omega$ such that eigenfunctions of the Dirichlet Laplacian in $\Omega$ can be extended only to an open neighborhood of $\overline{\Omega}$ but not to the whole $\mathbb{R}^{2}$. This has been found by J. Eckmann and C. Pillet in [15] and we reproduce their work here.

Given $L\in\mathbb{N}^{\ast}$ and $\mu\in\mathbb{R}_{+}:=(0,+\infty)$, we define the function $v$ such that

$$v(r,\theta)=\sum_{l=0}^{L-1}J_{\mu}(k\rho_{l})\cos(\mu\psi_{l}).$$ (11)

Here $(r,\theta)$ stands for the polar coordinates with $r>0$, $\theta\in(-\pi,\pi]$, $J_{\mu}$ is again the Bessel function of the first kind with order $\mu$ and $\rho_{l}>0$, $\psi_{l}\in(-\pi,\pi]$ satisfy the relations

$$\rho_{l}\cos\psi_{l}=a+r\cos(\theta+2\pi l/L),\qquad\qquad\rho_{l}\sin\psi_{l}=r\sin(\theta+2\pi l/L),$$

with $a\in\mathbb{R}$. In particular, the term for $l=0$ in (11) corresponds to the spherical wave function $J_{\mu}(kr)\cos(\mu\theta)$ with the center translated from $(0,0)$ to $(-a,0)$. If we further rotate the initial spherical wave function $J_{\mu}(kr)\cos(\mu\theta)$ clockwise by $2\pi/L$ (resp. $2\pi l/L$), we obtain the term in (11) for $l=1$ (resp. $l$ in general). More generally, we can consider functions $v$ of the form

$$v(r,\theta)=\sum_{l=0}^{L-1}b_{l}J_{\mu_{l}}(k\rho_{l})\cos\big{(}\mu_{l}(\psi_{l}-\phi_{l})\big{)},$$ (12)

with different Bessel orders $\mu_{l}$, rotational parameters $\phi_{l}$, and weights $b_{l}$, for each $l=0,\ldots,L-1$. We can also impose different translational and rotational effects to each term by defining

$$\rho_{l}\cos\psi_{l}=a_{l}+r\cos(\theta+\theta_{l})\qquad\text{and}\qquad\rho_{l}\sin\psi_{l}=r\sin(\theta+\theta_{l})$$

with different $a_{l}$ and $\theta_{l}$ for each $l=0,\ldots,L-1$. Throughout this subsection, we shall let $k$ be the first positive zero of $J_{\mu}$.

The following example, taken from [15], gives a bounded and simple connected domain $\Omega$ with analytic boundary where the Dirichlet Laplacian has an eigenfunction which is extendable in an open neighborhood of $\overline{\Omega}$ but not to the whole $\mathbb{R}^{2}$.

Example 2.11.

Let us consider the function $v$ in (11) with $\mu=3/2$, $L=3$, $a=0.6$. Since $\mu=3/2$, the function $\theta\mapsto\cos(\mu\theta)$ has a derivative which is not continuous at $\theta=\pi$ and then $v$ is real analytic only in $\mathbb{R}^{2}\setminus\Sigma$ with $\Sigma:=\{(r\cos\theta,r\sin\theta)\,|\,r\in[a,\infty),\,\theta=-\pi+2l\pi/L,\,l=0,\ldots,2\}$. In particular, $v$ is real analytic in the open ball $B_{a}$ centred at the origin and of radius $a$. Let us explain how to prove that there is a bounded and simple-connected domain $\Omega$ with real-analytic boundary such that $v$ is real-analytic in (an open neighborhood of) $\Omega$ and $v=0$ on $\partial\Omega$.
First notice that $(x,y)\mapsto J_{3/2}(kr)\cos(3\theta/2)$ is positive in $D_{1}=\{(r\cos\theta,r\sin\theta)\,|\,r\in(0,1),\theta\in(-\pi/3,\pi/3)\}$ and negative in $B_{1}\setminus\overline{D_{1}}$ (see Figure 5(a)). Using this property, we display in Figure 5(b) a yellow region where the three terms in the sum defining $v$ are all positive and a blue region where these three terms are all negative. We note by $\tilde{\Omega}$ the green region in which the sign of $v$ is a priori not clear. Looking at every direction $\theta\in(-\pi,\pi]$, with the intermediate value theorem, we can show that there is $(x_{\theta},y_{\theta})$ in $\tilde{\Omega}$ such that $v(x_{\theta},y_{\theta})=0$. Moreover, by exploiting that $v$ is real analytic in $\tilde{\Omega}$, we infer that there is a closed analytic curve $\Gamma\subset\tilde{\Omega}$ such that $v=0$ on $\Gamma$ (see Figure 5(c) for a plot of $\Gamma$). Denote by $\Omega$ the bounded open set surrounded by $\Gamma$.
Then $v$ is an eigenfunction of the Dirichlet Laplacian in $\Omega$ which can be analytically extended to $\tilde{\Omega}$ (actually to $\mathbb{R}^{2}\setminus\Sigma$) as a solution to the Helmholtz equation. Because of this property, we conclude that we have $\mathscr{S}_{\mathrm{NS}}\neq\emptyset$ for this $\Omega$. Let us emphasize however that $v$ can not be extended in the whole $\mathbb{R}^{2}$ analytically because its derivative with respect to $\theta$ is not continuous on $\Sigma$.

(a) Properties of $(x,y)\mapsto J_{3/2}(kr)\cos(3\theta/2)$: positive in the yellow regions and negative in the blue regions.

(b) Properties of $v$: in the yellow (resp. blue) region, the three terms in the sum defining $v$ are positive (resp. negative).

(c) Curves of constant value of the function $v$.
l
l

Figure 5: Illustration of the properties of the functions involved in Example 2.11.

We can modify Example 2.11 by taking different parameters in the definition of the function $v$ appearing in (11). This allows us to exhibit other analytic domains $\Omega$ where the Dirichlet Laplacian has eigenfunctions which are extendable only in a neighborhood of $\Omega$. In Figure 6, we display certain of these domains.

(a) $L=2$, $\mu=5/2$, $a=0.58$.

(b) $L=3$, $\mu=5/2$, $a=0.6$.

(c) $L=4$, $\mu=5/2$, $a=0.5$

Figure 6: Curves of constant value of the function $v$ for different parameters in (11). In the analytic domains enclosed by the curve $v=0$, the Dirichlet Laplacian has eigenfunctions which are only locally extendable.

Note that if we take the Bessel order $\mu$ in (11) in $\mathbb{N}^{\ast}$, we can find domains $\Omega$ for which the eigenfunctions of the Dirichlet Laplacian can be extended in the whole $\mathbb{R}^{2}$. We present such examples in Figure 7. Let us emphasize that the domains $\Omega$ may have analytic boundaries as in Figure 7(a), or Lipschitz boundaries as in Figures 2–7(c). The latter situation appears when two or more nodal curves of the function $v$ intersect. For example, in Figure 2, both $v$ and $\nabla v$ vanishes at the two points $\pm(a+1,0)$, but $\mathcal{H}v$, the Hessian of $v$, is different from zero at $\pm(a+1,0)$. Hence there are two nodal curves intersecting at $\pm(a+1,0)$. In Figure 7(c), at $(0,-\sqrt{2}/2)$ we have $\partial^{\alpha}v=0$ for all multi-indices $\alpha\in\mathbb{N}^{2}$ with $0\leq|\alpha|\leq 2$ but $\partial_{x}^{3}v\neq 0$. This explains why there are three nodal curves passing through $(0,-\sqrt{2}/2)$.

(a) $L=2$, $\mu=1$, $a=0.55$.
$\sqrt{2}$

(b) $L=2$, $\mu=1$,
$a=\frac{{k_{2}}-{k}}{2k}$ , $\phi_{0}=\pi/2=-\phi_{1}$.

(c) $L=2$, $\mu=1$,
$a=\sqrt{2}/2$, $\phi_{0}=\pi/4$, $\phi_{1}=3\pi/4$.

Figure 7: Curves of constant value of the function $v$ for different parameters in (11). Here we get analytic and Lipschitz domains where the Dirichlet Laplacian has eigenfunctions which extend to solutions of the Helmholtz equation in $\mathbb{R}^{2}$.

Similarly, we can find Lipschitz domains $\Omega$ for which there exist eigenfunctions are the Dirichlet Laplacian which can be extended in a neighborhood of $\overline{\Omega}$ but not to the whole $\mathbb{R}^{2}$ (see Figure 8). We refer the reader to the Appendix for the justification of existence of such nodal curves as shown in Figures 6–8.

(a) $L=2$, $\mu=7/2$, $a=\cos(\pi/7)$, $b_{0}=1=-b_{1}$.

(b) $L=2$, $\mu=5/2$, $a=1$, $\phi_{0}=\pi/5=-\phi_{1}$.

(c) $L=2$, $\mu=3/2$, $a=\sqrt{2}/2$, $\phi_{0}=\pi/{12}=-\phi_{1}$.

Figure 8: Curves of constant value of the function $v$ for different parameters in (11). Here we get Lipschitz domains where the Dirichlet Laplacian has eigenfunctions which are only locally extendable.

We conclude this paragraph with an important result which generalizes the constraint we have in Proposition 2.5 concerning the angle of a sector so that non-scattering frequencies can exist.

$$\begin{array}[]{|l}\mbox{Let $\Omega$ be a domain such that the Dirichlet Laplacian has an eigenfunction $v$ that}\\ \mbox{can be extended into a function $\tilde{v}$ satisfying $\Delta\tilde{v}+k^{2}\tilde{v}=0$ in a neighborhood of $\overline{\Omega}$.}\end{array}$$ (13)

First, for $\Omega$ as in (13), due to the implicit functions theorem, $\partial\Omega$ must be piecewise analytic. Additionally, we have the following statement (see also [13, Chapter V, §16]).

Proposition 2.12.

If $P\in\partial\Omega$ is a corner point of a domain $\Omega$ satisfying (13) with opening angle $\alpha$, then necessarily $\alpha\in(0,2\pi)\cap\mathbb{Q}\pi$.

Proof.

We start by observing that $\partial\Omega$ belongs to the nodal set of $v$, an eigenfunction of the Dirichlet Laplacian in $\Omega$. In an adapted system of coordinates, we can assume that $P$ coincides with the origin $O$. By the assumption that $v$ satisfies the Helmholtz equation in an open ball around $O$, we obtain that $v$ is of the form

$$v(x,y)=c_{N}J_{N}(kr)\cos(N\theta-\theta_{0})+f_{N+1}(x,y),$$

near $O$, where $(r,\theta)$ are the polar coordinates of $(x,y)$. Here $c_{N}\in\mathbb{R}\backslash\{0\}$ and $\theta_{0}\in[0,2\pi)$ are constants depending on $v$, $N\in\mathbb{N}^{\ast}$ is the vanishing order of $v$ at $P$ and $f_{N+1}=O(r^{N+1})$ is a real-analytic function in $(x,y)$. As a consequence, the Taylor series of $v$ near $O$ is of the form

$$v(x,y)=\tilde{c}_{N}r^{N}\cos(N\theta-\theta_{0})+\tilde{f}_{N+1}(x,y),$$

with $\tilde{c}_{N}\neq 0$ and $\tilde{f}_{N+1}$ real-analytic. Hence, there are exactly $N$ analytic nodal curves of $v$ intersecting at $P$, with tangential directions along angles in the set $\{(\pi/2+\kappa\pi+\theta_{0})/N\,|\,\kappa\in\mathbb{Z}\}$. Therefore we see that the angle between two tangential directions is necessarily equal to $\kappa^{\prime}\pi/N$ for some $\kappa^{\prime}\in\mathbb{N}$. ∎

2.3.2 Extendable Neumann eigenfunctions

What we did in §2.3.1 to find domains $\Omega$ where the Dirichlet Laplacian admits extendable eigenfunctions can not be directly applied for Neumann boundary conditions. Instead of looking for nodal curves, we should consider the characteristics associated with the gradient of solutions to the Helmholtz equation. Let us explain in the following.

Start with some non-zero real valued function $v$ satisfying $\Delta v+k^{2}v=0$ in some domain, say $\mathbb{R}^{2}$ for simplicity, for a given $k>0$. Pick some point $(x_{0},y_{0})$ such that $\nabla v(x_{0},y_{0})\neq 0$. Then consider the gradient system

$$\begin{array}[]{|rcl}X^{\prime}(t)&=&\nabla v(X(t))\\[3.0pt] X(0)&=&(x_{0},y_{0}).\end{array}$$ (14)

The Cauchy-Lipschitz theorem ensures that (14) admits a unique maximal solution $X$ defined on an interval $I\subset\mathbb{R}$ containing $0$. Set $\Gamma:=\{X(t)\,|\,t\in I\}\subset\mathbb{R}^{2}$ and introduce $\nu$ a unit normal vector to the smooth curve $\Gamma$ whose orientation is arbitrarily fixed. For $t\in I$, we have

$$\partial_{\nu}v(X(t))=\nu(X(t))\cdot\nabla v(X(t))=\nu(X(t))\cdot X^{\prime}(t)=0$$

because $X^{\prime}(t)$ and $\nu(X(t))$ are respectively tangent and normal to the orbit $t\mapsto X(t)$. Thus $v$ satisfies an homogeneous Neumann boundary condition on $\Gamma$. Therefore if we are able to find a bounded domain $\Omega$ whose boundary coincides with the union of orbits of Problem (14) for different initial conditions $(x_{0},y_{0})$, then this shows that $k$ belongs to the set $\mathscr{S}_{\mathrm{NS}}$ associated with $\Omega$. Let us illustrate the approach with some specific examples.

Example 2.13.

Consider the function

$$v(x,y)=-(\cos x+\cos y)$$ (15)

which satisfies $\Delta v+2v=0$ in $\mathbb{R}^{2}$. For this $v$, (14) reads

$$\begin{array}[]{|rcl}x^{\prime}(t)&=&\sin x(t)\\[3.0pt] y^{\prime}(t)&=&\sin y(t)\\[3.0pt] (x(0),y(0))&=&(x_{0},y_{0}).\end{array}$$ (16)

The stationary points of this system of ordinary differential equations are the $(m\pi,n\pi)$, $m,n\in\mathbb{Z}$. These are constant orbits. Moreover, one can check that

$$\{0\}\times(0,\pi),\quad\{\pi\}\times(0,\pi),\quad(0,\pi)\times\{0\},\quad(0,\pi)\times\{\pi\}$$

are also orbits of (16). Now pick some $(x_{0},y_{0})\in(0,\pi)^{2}$. Since different orbits cannot cross, we know that the one passing through $(x_{0},y_{0})$ stays inside $(0,\pi)^{2}$ and so the associated trajectory is global, i.e. defined for all $t\in\mathbb{R}$. It turns out that we can compute it explicitly and we find

$$x(t)=2\arctan(e^{t}\tan\frac{x_{0}}{2}),\qquad\quad y(t)=2\arctan(e^{t}\tan\frac{y_{0}}{2}).$$

After a few operations, we obtain that the corresponding orbit coincides with the curve

$$\Gamma(\delta):=\{(x,y(x))=2\arctan(\delta\tan(x/2)))\,|\,x\in(0,\pi)\}$$

where $\delta:=\tan(y_{0}/2)\tan(x_{0}/2)>0$ (see Figure 9 left for representations of different $\Gamma(\delta)$). For any $0<\delta_{1}<\delta_{2}$, denote $\Omega(\delta_{1},\delta_{2})$ the domain located between $\Gamma(\delta_{1})$ and $\Gamma(\delta_{2})$. Then $v$ is an eigenfunction of the Neumann Laplacian in $\Omega(\delta_{1},\delta_{2})$ (see Figure 9 right for pictures of some $\Omega(\delta_{1},\delta_{2})$). Since $v$ satisfies $\Delta v+2v=0$ in $\mathbb{R}^{2}$, we deduce that $2$ belongs to $\mathscr{S}_{\mathrm{NS}}$ for this domain.
Note that we also have $2\in\mathscr{S}_{\mathrm{NS}}$ for the domain delimited for example by the curves

$$\{0\}\times(0,\pi),\qquad\{\pi\}\times(0,\pi),\qquad\Gamma(\delta),\qquad\mbox{for any $\delta>0$}.$$

Figure 9: Left: curves $\Gamma(\delta)$ for $\delta=0.5,1,4$. Right: function $v$ defined in (15) in a few domains where it is an eigenfunction of the Neumann Laplacian associated with the eigenvalue $2$.

At this stage, let us make a few comments. First, when one enlarges a bounded domain, due to the min-max principle, the eigenvalues of the Dirichlet Laplacian decrease. This is not true for the eigenvalues of the Neumann Laplacian, roughly speaking because one cannot extend them by zero. This is well illustrated by the above example. We have a family of continously deformed shapes which all have $2$ among the eigenvalues of the Neumann Laplacian.

Second, this example shows that one can find domains where the Neumann Laplacian has eigenfunctions which extend as solutions to the Helmholtz solution in the whole space and which admit corners with arbitrary apertures in $(0,2\pi)$. This is very different from the Dirichlet case where in order to be able to extend the eigenfunctions of the Laplace operator, the aperture of the corners must be in $(0,2\pi)\cap\mathbb{Q}\pi$ (see Proposition 2.12).

The reason we can have such domains $\Omega$ with corners of arbitrary aperture in Example 2.13 is that $\mathcal{H}v(O)$, the Hessian of $v$ at the stationary point $O=(0,0)$, is equal to the identity matrix. As a consequence, the linearized problem of (16) at $O$ simply writes $U^{\prime}(t)=U(t)$ with $U=(U_{x},U_{y})^{\top}$. Solving it, we get $U_{x}(t)=a_{x}\,e^{t}$, $U_{y}(t)=a_{y}\,e^{t}$ for some constants $a_{x}$, $a_{y}$. Thus we see that the curve $\{(U_{x}(t),U_{y}(t)),\,t\in(-\infty;0)\}$ is a line which can take any direction by choosing properly $a_{x}$, $a_{y}$. This explains why the orbits of (16) can leave from $O$ asymptotically with any angle.

Let us consider a second example.

Example 2.14.

Start from the function $v$ such that $v(x,y)=\cos x\cos 2y$. It satisfies $\Delta v+5v=0$ in $\mathbb{R}^{2}$. The set of stationary points associated with (14) is $\{(m\pi,n\pi/2),(\pi/2+m\pi,\pi/4+n\pi/2)\,|\,m,n\in\mathbb{Z}\}$. By computing the Hessian $\mathcal{H}v$, one finds that the $(\pi/2+m\pi,\pi/4+n\pi/2)$ are saddle points. At saddle points, orbits of (14) can meet the equilibrium point at only two particular directions (see Figure 10(b)). On the other hand, for $m,n\in\mathbb{Z}$, we have $\mathcal{H}v(m\pi,n\pi/2)=(-1)^{m+n-1}\mathrm{diag}(1,4)$. Hence the orbits of (14) are all tangential to the horizontal line at $(0,0)$ except one that is tangential to the vertical line. Indeed, for example for $m=n=0$, the linearized problem associated with (14) writes $U^{\prime}(t)=-\mathrm{diag}(1,4)\cdot U(t)$ so that its solution is given by $U_{x}(t)=a_{x}\,e^{-t}$, $U_{y}(t)=a_{y}\,e^{-4t}$ for some constants $a_{x}$, $a_{y}$. Thus we see that except when $a_{x}=0$, the curve $\{(U_{x}(t),U_{y}(t)),\,t\in(0;+\infty)\}$ is tangential to the horizontal line when $t\to+\infty$.
We can make this more explicitly by computing the orbits. First we find that

$$(0,0),\quad(\pi,0),\quad(\pi/2,0),\quad(\pi/2,\pi),\quad(0,\pi)\times\{0\},\quad\{\pi\}\times(0,\pi/2),\quad(0,\pi)\times\{\pi/2\},\quad\{0\}\times(0,\pi/2)$$

are particular orbits of (14). Then we obtain that the (global) trajectory $(x(t),y(t))$ passing through $(x_{0},y_{0})\in(0,\pi)\times(0,\pi/2)$ is such that

$$\sin(2y(t))=\delta\sin^{4}x(t)$$

for $\delta>0$. We show in Figure 10 the normalized gradient flow and some of the “Neumann curves”. The latter enclose domains $\Omega$ for which $v$ is a Neumann eigenfunction so that $5$ belongs to $\mathscr{S}_{\mathrm{NS}}$. We get that $\Omega$ can be of class $\mathscr{C}^{1}$, Lipschitz but may also possess cusps.

(a) Normalized gradient $\nabla v/|\nabla v|$.

(b) Neumann curves.

Figure 10: Case $v(x,y)=\cos x\cos 2y$.

Example 2.15 (Neumann eigenfunctions extendable in a neighborhood of $\overline{\Omega}$ but not to the whole $\mathbb{R}^{d}$).

Let us defined $v$ as in (12) with $L=2$, $\mu=7/2$, $a=\cos(\pi/7)$ and $b_{0}=1=-b_{1}$, with $k$ being the first positive zero of $J_{\mu}$ (see Figure 8(a)). We can calculate the gradient $\nabla v$ explicitly, say, when $|x_{1}|<a$, as

$$\nabla v(x,y)=\sum_{l=0}^{1}(-1)^{l}\Big{(}kJ^{\prime}_{\mu}(k\rho_{l})\cos(\mu\psi_{l})\nabla\rho_{l}-\mu J_{\mu}(k\rho_{l})\sin(\mu\psi_{l})\nabla\psi_{l}\Big{)},$$

with

$$\rho_{l}=\sqrt{\left(a+(-1)^{l}x_{1}\right)^{2}+x_{2}^{2}},\qquad\text{and}\qquad\psi_{l}=(-1)^{l}\arcsin\frac{x_{2}}{\rho_{l}},\qquad l=0,1.$$

It can be shown that $\nabla v=0$ at $(0,\pm\sin(\pi/7))$ and $(\pm t_{0},0)$, where $t_{0}$ is the (only) solution to $J_{\mu}^{\prime}\left(k(a+t_{0})\right)+J_{\mu}^{\prime}\left(k(a-t_{0})\right)=0$ in the interval $(0,a)$. The normalized gradient flow is shown in Figure 11(a). Moreover, by (numerically) solving the system of ODEs (14), we can construct curves on which $\partial_{\nu}v=0$, and some of those curves enclose bounded $\mathscr{C}^{1}$ or Lipschitz domains (with corners of opening angle $\pi/2$), or domains with cusps. We display some of these domains in Figure 11, where the starting points $X(0)=(x_{0},y_{0})$ satisfy $x_{0}=0$ and $y_{0}\in\{\alpha\sin(\pi/7)\,|\,|\alpha|<1\}$. In the numerics, we take $\pm\alpha\in\{0,0.1,0.5,0.8,0.9,0.9999\}$. The curves on the left side of the starting points are obtained by solving (14) and letting $t\to+\infty$. We get $\lim_{t\to+\infty}X(t)=(-t_{0},0)$. Those on the right side are obtained by solving the “backward” system $X^{\prime}(t)=-\nabla v(X(t))$, in which case $\lim_{t\to+\infty}X(t)=(t_{0},0)$. An explicit computation gives $\mathcal{H}v(\pm t_{0},0)=\pm\text{diag}\,(\lambda_{1},\lambda_{2})$ with $0<\lambda_{2}<\lambda_{1}$. As a consequence, though this does not appear clearly in Figure 11(b) because here we do not zoom at $(\pm t_{0},0)$, all the orbits except the horizontal one, are tangential to the vertical axis at $(\pm t_{0},0)$.

(a) Normalized gradient $\nabla v/|\nabla v|$.

(b) Neumann curves.

Figure 11: Case $v$ as in (12) with $L=2$, $\mu=7/2$, $a=\cos(\pi/7)$, $b_{0}=1=-b_{1}$.

3 Non-scattering with anisotropic materials

In this section, we provide examples of anisotropic materials, more precisely materials for which $A$ in (2) is a matrix valued function, that support non-scattering frequencies. For some of them, non-scattering occurs at all $k>0$ so that in particular the set of transmission eigenvalues $\mathscr{S}_{\mathrm{TE}}$ contains $(0,+\infty)$.

3.1 Non-scattering via diffeomorphism

We first look for non-scattering waves and media by using diffeomorphisms following for example [18] (see also [10]). Let $v$ be an entire solution to the Helmholtz equation, i.e. such that

$$\Delta v+k^{2}v=0\qquad\mbox{in $\mathbb{R}^{d}$}.$$

Given a bounded Lipschitz domain $\Omega$ in $\mathbb{R}^{d}$, let $\Psi$ be a $\mathscr{C}^{m}(\overline{\Omega})$, $m\geq 2$, diffeomorphism satisfying $\Psi=\operatorname{Id}$ on $\partial\Omega$. Define

$${A}=\frac{D\Psi D^{T}\Psi}{\left\lvert\det D\Psi\right\rvert}\circ\Psi^{-1}\quad\text{and}\quad{q}=\frac{1}{\left\lvert\det D\Psi\right\rvert}\circ\Psi^{-1}\qquad\mbox{in $\Omega$},$$ (17)

where $D\Psi$ is the Jacobian matrix of $\Psi$. Note in particular that $A$ is symmetric. Then the function $u:=v\circ\Psi^{-1}$ satisfies $u=v$ on $\partial\Omega$ and

$$\begin{split}\mathrm{div}(A\nabla u)+k^{2}qu=0\quad\mbox{in $\Omega$},&\qquad\partial_{\nu}^{{A}}u=\partial_{\nu}v\quad\mbox{on $\partial\Omega$,}\end{split}$$ (18)

with $\partial_{\nu}^{A}\cdot:=\nu\cdot A\nabla\cdot$. As a consequence, we have the following lemma.

Lemma 3.1.

Let $\Omega$ be a bounded Lipschitz domain in $\mathbb{R}^{d}$ and let $\Psi$ be a $\mathscr{C}^{2}(\overline{\Omega})$ diffeomorphism on $\overline{\Omega}$ satisfying $\Psi=\operatorname{Id}$ on $\partial\Omega$. Define $A$, $q$ by (17). Then $A$, $q$ satisfy assumptions (3) and $(A,q;\Omega)$ is non-scattering for any incident waves at any frequencies.

Remark 3.2.

Notice for $A$ defined in (17) we always have $\det A=\left\lvert\det D\Psi\right\rvert^{2-d}\circ\Psi^{-1}$. In particular $\det A=1$ in dimension two and $\det A=q$ in dimension three. This may shed light to the verification of optimality for results on the discreteness of transmission eigenvalues; See also [21, 24].

As observed in [10], given any bounded Lipschitz domain $\Omega\subset\mathbb{R}^{d}$, one can construct infinitely many diffeomorphisms that satisfy the conditions in Lemma 3.1. For instance, let $\Phi:\mathbb{R}^{d}\to\mathbb{R}^{d}$ be a $\mathscr{C}^{2}$ map such that $\Phi=0$ in $\mathbb{R}^{d}\setminus\Omega$. Then the mapping $\Psi$ defined as

$$\Psi=\mathrm{Id}+\varepsilon\Phi.$$ (19)

is a $\mathscr{C}^{2}(\overline{\Omega})$ diffeomorphism for $\varepsilon>0$ sufficiently small that satisfies $\Psi=\operatorname{Id}$ on $\partial\Omega$.

Next, we give two explicit non-scattering examples (where $\varepsilon$ in (19) does not have to be small) with $\Omega$ being a square and a disk, respectively.

Example 3.3 (Geometry with corners).

Let $\Omega=(-1,1)^{2}\subset\mathbb{R}^{2}$. For a fixed constant $\alpha\in(-1/2,1/2)$, define

$$\Psi(x,y)=\left(x+\alpha(1-x^{2})(1-y^{2})\,,\,y\right).$$

It can be verified that $\Psi:\overline{\Omega}\to\overline{\Omega}$ is a $\mathscr{C}^{\infty}$ diffeomorphism with $\Psi|_{\partial\Omega}=\operatorname{Id}$. Then the medium $(A,q;(-1,1)^{2})$ as in (17) is non-scattering for all incident waves at all frequencies. We can express $A$ and $q$ in a more explicit form as

$$A\circ\Psi(x,y)=\frac{1}{1-2\alpha x(1-y^{2})}\left(\begin{array}[]{cc}\left(1-2\alpha x(1-y^{2})\right)^{2}+4\alpha^{2}y^{2}(1-x^{2})^{2}&-2\alpha y(1-x^{2})\\ -2\alpha y(1-x^{2})&1\end{array}\right)$$

and

$$q\circ\Psi=\frac{1}{1-2\alpha x(1-y^{2})}.$$

Example 3.4 (Non-spherically stratified disk).

Let $f\in\mathscr{C}^{2}_{0}(\mathbb{R}^{+})$ satisfy $f(1)=0$. Define $\Psi$, in polar coordinates, as

$$\Psi(r,\theta)=(r,\theta+f(r)).$$

Then $\Psi:\overline{B_{1}}\to\overline{B_{1}}$ is a $\mathscr{C}^{2}$ diffeomorphism with $\Psi|_{\partial\Omega}=\operatorname{Id}$ (this latter property comes from the constraint $f(1)=0$). Moreover, $\Psi$ is orientation- and area-preserving with $\det D\Psi=1$. As a consequence, the inhomogeneity $(A,1;B_{1})$ with $A$ as in (17) is non-scattering for all incident waves at all frequencies. Note that we find

$$A(x,y)=\mathrm{Id}-\frac{f^{\prime}(r)}{r}\left(\begin{array}[]{cc}2xy&y^{2}-x^{2}\\ y^{2}-x^{2}&-2xy\end{array}\right)+\left(f^{\prime}(r)\right)^{2}\left(\begin{array}[]{cc}y^{2}&-xy\\ -xy&x^{2}\end{array}\right),\quad r=\left(x^{2}+y^{2}\right)^{1/2}.$$

On the structure of the anisotropy induced by diffemorphism transforms.

The following two results reveal an interesting geometry property of diffeomorphisms satisfying the conditions in Lemma 3.1.

Lemma 3.5.

Let $\Omega$ and $\Psi$ be as in Lemma 3.1. Then we must have

$$\Big{(}\frac{D^{T}\Psi}{\left\lvert\det D\Psi\right\rvert}-\mathrm{Id}\Big{)}\nu=0\qquad\mbox{on $\partial\Omega$}.$$

Proof.

Let $j\in\{1,\ldots,d\}$ be fixed. Let $v(\mathrm{x})=x_{j}$ and let $u=v\circ\Psi^{-1}$. Then $u$, $v$ satisfy

$$\mathrm{div}(A\nabla u)=0\quad\text{and}\quad\Delta v=0\quad\mbox{in $\Omega$},\qquad u=v\quad\text{and}\quad\partial_{\nu}^{A}u=\partial_{\nu}v\quad\mbox{on $\partial\Omega$},$$

with $A$ defined in (17). Furthermore, by straightforward calculations we have

$$\partial_{\nu}v=\nu^{T}\vec{e}_{j}\qquad\text{and}\qquad\partial_{\nu}^{A}u=\nu^{T}\frac{D\Psi D^{T}\Psi}{\left\lvert\det D\Psi\right\rvert}(D^{T}\Psi^{-1})\,\vec{e}_{j}=\nu^{T}\frac{D\Psi}{\left\lvert\det D\Psi\right\rvert}\,\vec{e}_{j}\qquad\mbox{on $\partial\Omega$}.$$

The proof is complete by taking all $j\in\{1,\ldots,d\}$. ∎

Combining this result with [10, Theorem 2.1], we obtain the following statement when $\Omega$ is not smooth.

Proposition 3.6.

Let $\Omega$ be a $\mathscr{C}^{1,\alpha}$ bounded domain in $\mathbb{R}^{d}$, and let $\Psi$ be a $\mathscr{C}^{m+2}$ (resp. $\mathscr{C}^{\infty}$ or (real) analytic) diffeomorphism on $\overline{\Omega}$ with $m\geq 1$ satisfying $\Psi|_{\partial\Omega}=\operatorname{Id}$. Then

$$\left(D\Psi-\mathrm{Id}\right)\nu=0\qquad\mbox{for\quad$\mathrm{x}_{0}\in\partial\Omega$},$$ (20)

where $\mathrm{x}_{0}\in\partial\Omega$ is such that $\partial\Omega$ is not $\mathscr{C}^{m,\alpha}$ (resp. $\mathscr{C}^{\infty}$ or (real) analytic) in any neighborhood of $\mathrm{x}_{0}$.

Proof.

Let $\mathrm{x}_{0}\in\partial\Omega$ be such that $\partial\Omega$ is not $\mathscr{C}^{m,\alpha}$ (resp. $\mathscr{C}^{\infty}$ or (real) analytic) in any neighborhood of $\mathrm{x}_{0}$. Assume by contradiction that $(D\Psi-\mathrm{Id})\nu\neq 0$ at $\mathrm{x}_{0}$. Let

$$A=\frac{D\Psi D^{T}\Psi}{\left\lvert\det D\Psi\right\rvert}\quad\mbox{ in $\overline{\Omega}$}.$$

Then by Lemma 3.5 we have

$$(A-\mathrm{Id})\,\nu=(D\Psi-\mathrm{Id})\,\nu\qquad\mbox{on $\partial\Omega$},$$

and so $(A-\mathrm{Id})\nu\neq 0$ at $\mathrm{x}_{0}\in\partial\Omega$. Exploiting this, we can find a solution $v$ to the homogeneous Helmholtz equation in $\mathbb{R}^{d}$ such that $\nu^{T}(A-\mathrm{Id})\nabla v\neq 0$ at $\mathrm{x}_{0}$. But from Lemma 3.1 we know that $v$ is non-scattering for $(A,q;\Omega)$ with $A$, $q$ given in (17). Hence by [10, Theorem 2.1] we must have that $\partial\Omega$ is $\mathscr{C}^{m,\alpha}$ (resp. $\mathscr{C}^{\infty}$ or (real) analytic) near $\mathrm{x}_{0}$, which leads to a contradiction. ∎

Remark 3.7.

In particular, we show in the proof that the matrix $A$ defined in (17) must satisfy

$$A\,\nu=\nu$$

at the considered “non-smooth” point $\mathrm{x}_{0}$ of $\partial\Omega$

3.2 Other explicit non-scattering examples

In this section, we are focused on non-scattering examples where the total and incident fields are identical in the whole $\mathbb{R}^{d}$.

Example 3.8.

Let $\Omega=(0,\pi)^{2}$ and

$$A=\left(\begin{matrix}a_{1}&0\\ 0&a_{2}\end{matrix}\right),\qquad\quad q=q_{0}\qquad\mbox{in $\Omega$,}$$ (21)

with $a_{1},a_{2},q_{0}\in\mathbb{R}_{+}$. Assume that these coefficients satisfy one of the following assumptions:
$\bullet$ $a_{1}\neq 1$, $a_{2}\neq 1$ and there exist $m,n\in\mathbb{N}$, with $(m,n)\neq(0,0$), such that $m^{2}(q_{0}-a_{1})=n^{2}(a_{2}-q_{0})$;
$\bullet$ $a_{1}=1$, $a_{2}\neq 1$ and $(q_{0}-1)(a_{2}-1)>0$;
$\bullet$ $a_{1}\neq 1$, $a_{2}=1$ and $(q_{0}-1)(a_{1}-1)>0$.
Then $\mathscr{S}_{\mathrm{NS}}$ contains an unbounded sequence which accumulates at $+\infty$.

Remark 3.9.

It is known for the last two cases in Example 3.8 that $\mathscr{S}_{\mathrm{TE}}$ is at most countable, and hence such $(A,q;\Omega)$ cannot be non-scattering for all frequencies. The same is true for the first case provided $a_{1}q_{0}\neq 1$ or $a_{2}q_{0}\neq 1$ (see [21, 24]); Otherwise if $a_{1}q_{0}=1$ or $a_{2}q_{0}=1$ but $q_{0}\neq 1$, the discreteness of $\mathscr{S}_{\mathrm{TE}}$ is to our best knowledge still an open question.

Proof.

Suppose that $a_{1},a_{2},q_{0}$ satisfy the first set of assumptions. Then it can be verified straightforwardly that for any $\kappa\in\mathbb{N}^{\ast}$, the functions $u$, $v$ such that

$$u(\mathrm{x})=v(\mathrm{x})=\cos(\kappa mx)\cos(\kappa ny)\quad\mbox{in $\Omega$}$$

solve the interior transmission eigenvalue problem (ITEP) (6) with $k=\kappa\sqrt{m^{2}+n^{2}}$.
Now assume that $a_{1},a_{2},q_{0}$ satisfy the second set of assumptions (the third one can be dealt with similarly). Then one can verify that for any $m\in\mathbb{N}^{\ast}$, the functions $u$, $v$ such that

$$u(\mathrm{x})=v(\mathrm{x})=\big{(}c_{1}\cos(bx)+c_{2}\sin(bx)\big{)}\cos(my)\quad\mbox{in $\Omega$},$$

with $b=\sqrt{k^{2}-m^{2}}\in\mathbb{R}_{+}\cup i\mathbb{R}_{+}$ and where $c_{1}$, $c_{2}$ are non zero constants, solve (6) for $k=m\sqrt{(a_{2}-1)/(q_{0}-1)}$. ∎

Remark 3.10.

By scaling and/or rigid change of coordinates, we can adapt Example 3.8 to the case of any rectangle, with properly modified form of $A(\mathrm{x})$ (not necessarily diagonal) and conditions between $A$ and $q$. Analogous arguments and results apply when $\Omega$ is a cuboid in $\mathbb{R}^{3}$.

The case in Example 3.8 but for $q=a_{j}=1$ for $j=1$ or $2$ can be generalized to the following situation where $q\equiv 1$ and $A-\mathrm{Id}$ is rank deficient. In this situation we lost the discreteness of the sets $\mathscr{S}_{\mathrm{NS}}$ and $\mathscr{S}_{\mathrm{TE}}$. A particular case of this example appears in [21, p.1170].

Example 3.11.

Given any bounded Lipschitz domain $\Omega$ in $\mathbb{R}^{d}$, suppose that $A=A(\mathrm{x})$ satisfies

$$A(\mathrm{x})=U\left(\begin{array}[]{cc}1&0\\ 0&A_{n-1}(\mathrm{x})\end{array}\right)U^{T},\qquad\mathrm{x}\in\overline{\Omega},$$

with some $(n-1)\times(n-1)$ symmetric positive definite matrix $A_{n-1}$ and a constant orthogonal matrix $U$. Importantly, note that in this case $1$ is an eigenvalue of $A(\mathrm{x})$ for all $\mathrm{x}\in\overline{\Omega}$. Then any $k>0$ is a non-scattering frequency for $(A,1;\Omega)$.

Proof.

For any non-zero constants $c_{1}$, $c_{2}$, define $v_{0}(\mathrm{x})=v_{0}(x_{1},x^{\prime})=c_{1}\cos(kx_{1})+c_{2}\sin(kx_{1})$. It can be verified straightforwardly that for any bounded Lipschitz domain $\tilde{\Omega}$ and any fixed $k>0$, $(u_{0},v_{0})=(v_{0},v_{0})$ is a pair of eigenfunctions to the ITEP (6) for $(\tilde{A},1;\tilde{\Omega})$ with $\tilde{A}=\mathrm{diag}(1,A_{n-1})$ (note in particular that $(\tilde{A}-I)\nabla v_{0}\equiv 0$). Then consider the change of coordinates $\mathrm{y}=U\mathrm{x}$ for $\mathrm{x}\in\mathbb{R}^{d}$ and define $v$ such that $v(\mathrm{y}):=v_{0}\big{(}U^{-1}\mathrm{y}\big{)}$. The pair $(u,v)=(v,v)$ satisfies (6) for $(A,1;\Omega)$ with $\Omega=\{\mathrm{y}=U\mathrm{x}\,|\,\mathrm{x}\in\tilde{\Omega}\}$. Moreover, $v$ is an entire solution to the Helmholtz equation. The proof is complete. ∎

The following can be viewed as a generalization of Example 3.8 when $q=a_{j}\neq 1$ for $j=1$ or $2$.

Example 3.12.

Let $\Omega=(b_{1},b_{2})\times D\subset\mathbb{R}^{d}$ with $D$ an open interval if $d=2$ or a bounded Lipschitz domain in $\mathbb{R}^{2}$ if $d=3$. Then there are infinitely many non-scattering frequencies for $(A,q;\Omega)$ with

$$A(\mathrm{x})=\left(\begin{array}[]{cc}a_{0}&0\\ 0&A_{n-1}(\mathrm{x})\end{array}\right)\qquad\text{and}\qquad q(\mathrm{x})=a_{0},\qquad\mathrm{x}\in\Omega,$$

for some positive constant $a_{0}$ and some $(n-1)\times(n-1)$ symmetric positive definite matrix $A_{n-1}$.

Proof.

It can be verified straightforwardly that $v(\mathrm{x})=\cos(k(x_{1}-b_{1}))$ is a non-scattering incident field with $u=v$ being the total field for $k=m\pi/|b_{2}-b_{1}|$, $m\in\mathbb{N}^{\ast}$. ∎

Remark 3.13.

We notice from Example 3.11 that if $a_{0}=1$ in Example 3.12, the sets $\mathscr{S}_{\mathrm{NS}}$ and hence $\mathscr{S}_{\mathrm{TE}}$ coincide with $\mathbb{R}^{+}$. However, in some other cases, for example when both $a_{0}$ and the smallest eigenvalue of $A_{n-1}(\mathrm{x})$ are larger than $1$ for all $\mathrm{x}$ in a small neighborhood of $\partial\Omega$, the sets $\mathscr{S}_{\mathrm{TE}}$ and hence $\mathscr{S}_{\mathrm{NS}}$ are discrete [5, 24]. The discreteness of $\mathscr{S}_{\mathrm{TE}}$ or $\mathscr{S}_{\mathrm{NS}}$ in general is unknown for Example 3.12 especially if $\det A=1$ on $\partial\Omega$.

Appendix

The following figures illustrate how to justify the existence of analytic nodal curves as shown in Figures 6–8. Below the corresponding functions are positive in the yellow regions and negative in the blue regions.

As an example of such justification, let us look at the nodal set of Figure 6(a). That is, when $L=2$, $\mu=5/2$ and $a=0.58$ in (11). Based on the “generating function” $J_{\mu}(kr)\cos(\mu\theta)$, we can identify the positive and negative regions of each terms $l=0,\ldots,L-1$ in (11). See Figure 12(a) for an example for the term $l=0$ in (11), where the dashed lines are the nodal curves of the terms $l=1,\ldots,L-1$ in (11). Then we can find certain regions where $v$ in (11) has to be either positive/negative, say, the domains where every term is of the same sign (see Figure 12(b)). Therefore, there must be a nodal set in between the strictly positive and the strictly negative regions. The justification for Figures 6(b), 6(c), and 7(a) is similar.

(a) The term $l=0$ in (11).

(b) $v$ in (11).

(c) The nodal curves of $v$.

Figure 12: Justification for Figure 6(a): $L=2$, $\mu=5/2$, $a=0.58$.

(a) The term $l=0$ in (11).

(b) $v$ in (11).

(c) The nodal curves of $v$.

Figure 13: Justification for Figure 6(b): $L=3$, $\mu=5/2$, $a=0.6$.

(a) The term $l=0$ in (11).

(b) $v$ in (11).

(c) The nodal curves of $v$.

Figure 14: Justification for Figure 6(c): $L=4$, $\mu=5/2$, $a=0.5$.

(a) The term $l=0$ in (11).

(b) $v$ in (11).

(c) The nodal curves of $v$.

Figure 15: Justification for Figure 7(a): $L=2$, $\mu=1$, $a=0.55$.

The situation for the case of Figure 2 is slightly different. First, we can identify the positive and negative regions for each term in (12) as before, and hence some one-sign regions for $v$ in (12) (see Figures 16(a) and 16(b)). In particular, thanks to the symmetry of the two terms in (12), we have that $v=0$ on the axis $\{y=0\}$. Additionally, we observe that $v=0$ at the two points $(0,\pm\sqrt{1-a^{2}})$, where in fact both terms in (12) are zero. Moreover, it can be verified straightforwardly that $\nabla v\neq 0$ at $(0,\pm\sqrt{1-a^{2}})$. So there is exactly one (real-analytic) nodal curve of $v$ passing through the point $(0,\sqrt{1-a^{2}})$. This curve will continue, say, to the right side of $(0,\sqrt{1-a^{2}})$ until it hits some point $(x_{1},0)$ for some $x_{1}\in[1-a,1+a]$. On the other hand, using the fact that $a=(k_{2}/k-1)/2$ with $k_{2}$ the second positive zero of $J_{2}$, it can be verified that $\nabla v=0$ at the two points $\pm(a+1,0)$ and $\nabla v\neq 0$ at any $(x_{1},0)$ with $x_{1}\in(-a-1,a+1)$. Therefore, the nodal curve of $v$ generated from $(0,\sqrt{1-a^{2}})$ must also pass through $(a+1,0)$. In fact, one can also verify that the Hessian of $v$ satisfies $\mathcal{H}v\neq 0$ at $\pm(a+1,0)$. Hence there are exactly two (real-analytic) nodal curves passing at $(a+1,0)$: one is the line $\{y=0\}$ and the other is the one coming from $(0,\sqrt{1-a^{2}})$. Applying analogous arguments for the other three quadrant we can then complete the justification of a closed (real-analytic) nodal curve of $v$.

(a) The term $l=0$ in (12).

(b) $v$ in (12).

(c) The nodal curves of $v$.

Figure 16: Justification for Figure 2: $L=2$, $\mu=1$, $a=\frac{{k_{2}}-{k}}{2k}$, $\phi_{1}=\pi/2=-\phi_{2}$.

In Figure 7(c), we make use of the function $v$ in (12) with $\mu=1$, $L=2$, $a=\sqrt{2}/2$, $\phi_{1}=\pi/4$ and $\phi_{2}=3\pi/4$. The choice of $\phi_{1}$ and $\phi_{2}$ ensures the symmetry of $v$ with respect to the $y$-axis. In particular, $v=0$ on $\{x=0\}$. In addition, thanks to the value of $a$, one can verify that at the point $P_{0}=(0,-\sqrt{2}/2)$ there holds $v=\partial_{x}v=\partial_{y}v=\partial^{2}_{ij}v=0$ for all $i,j\in\{x,y\}$ but $\partial_{x}^{3}v\neq 0$. Hence, there are exactly three nodal curves passing through $(0,-\sqrt{2}/2)$. Furthermore, one can verify straightforwardly that $v=0$ while $\nabla v\neq 0$ at the points $P_{j}$, $j=1,\ldots,9$, where the coordinates of each point can be calculated explicitly in terms of $a$, $k$ and $k_{2}$. Finally, based on the locations of strictly positive and negative regions we can conclude to the existence of real-analytic nodal curves for $v$ that enclose some bounded Lipschitz domains, as shown Figure 17(c).

(a) The term $l=0$ in (12).

(b) $v$ in (12).

(c) The nodal curves of $v$.

Figure 17: Justification for Figure 7(c): $L=2$, $\mu=1$, $a=\sqrt{2}/2$, $\phi_{1}=\pi/4$, $\phi_{2}=3\pi/4$.

The justification of the results appearing in Figure 8 is similar.

(a) The term $l=0$ in (12).

(b) $v$ in (12)

(c) The nodal curves of $v$.

Figure 18: Justification for Figure 8(a): $L=2$, $\mu=7/2$, $a=\cos(\pi/7)$, $b_{0}=1=-b_{1}$.

(a) The term $l=0$ in (12).

(b) $v$ in (12).

(c) The nodal curves of $v$.

Figure 19: Justification for Figure 8(b): $L=2$, $\mu=5/2$, $a=1$, $\phi_{0}=\pi/5=-\phi_{1}$.

(a) The term $l=0$ in (12).

(b) $v$ in (12).

(c) The nodal curves of $v$.

Figure 20: Justification for Figure 8(c): $L=2$, $\mu=3/2$, $a=\sqrt{2}/2$, $\phi_{0}=\pi/12=-\phi_{1}$.