In this note, we exhibit a three dimensional structure that permits to guide waves. This structure is obtained by a geometrical perturbation of a 3D periodic domain that consists of a three dimensional grating of equi-spaced thin pipes oriented along three orthogonal directions. Homogeneous Neumann boundary conditions are imposed on the boundary of the domain. The diameter of the section of the pipes, of order , is supposed to be small. We prove that, for small enough, shrinking the section of one line of the grating by a factor of () creates guided modes that propagate along the perturbed line. Our result relies on the asymptotic analysis (with respect to ) of the spectrum of the Laplace-Neumann operator in this structure. Indeed, as tends to , the domain tends to a periodic graph, and the spectrum of the associated limit operator can be computed explicitly.
Keywords : guided waves, periodic media, spectral theory.
AMS codes : 78M35, 35J05, 58C40
Existence of guided waves due to a lineic perturbation of a 3D periodic medium
Bérangère Delourme, Patrick Joly, Elizaveta Vasilevskaya
: Université Paris 13, Sorbonne Paris Cité, LAGA, UMR 7539, 93430 Villetaneuse, France
: POEMS, UMR 7231 CNRS/ENSTA/Inria, ENSTA ParisTech, 828 boulevard des Maréchaux, 91762 Palaiseau Cedex, France
1 Statement of the problem
Let , and be three Lipschitz bounded domains of of same area () containing the origin , let be a parameter (that is going to be small), and let , and be three positive real numbers. We denote by , the standard basis of . For any , we consider the three dimensional domain defined by
which is an unbounded cylinder of constant cross section . It is infinite along the direction (invariant with respect to ) and contains the point . Similarly, for any , we define the domains
and we consider the periodic domain given by
The domain is a three dimensional grating of equi-spaced parallel pipes (of constant cross section) oriented along the three orthogonal directions , and . It is -periodic with respect to . Moreover, the points , , belong to .
In order to create guided modes, we introduce a linear defect (see --) in the periodic structure by modifying the section size of one pipe of the grating (it is conjectured that guided modes cannot appear in the purely periodic structure, see  for the proof in the case of a symmetric medium). More precisely, we assume that the domain is replaced with the domain
where is a positive parameter. In other words, we enlarge () or shrink () the section of one pipe of the domain by a factor (see Fig 3). The corresponding perturbed domain is denoted by . Its precise definition is given by
is still -periodic with respect to . However, the presence of the perturbed pipe breaks the periodicity with respect to and . We emphasize that the domain (as well as ) tends to a 3D periodic graph as tends to .
We look for guided modes, i.e. solutions to the wave equation in , satisfying homogeneous Neumann boundary conditions on (see  for the investigation of the Dirichlet case), that propagate along the defect pipe (i.e. in the direction) but stay confined in the transversal directions. More precisely, denoting by the restriction of the domain to the band ,
we search solutions of the form , where is a real parameter and is an -periodic function in . In fact, it is easily seen that the -quasiperiodic fonction is an eigenfunction of the operator
To study the spectral properties of , we investigate its (formal) limit as tends to . The operator is defined on the limit graph (see Fig. 3) and its spectrum can be explicitly computed. In particular, its spectrum has infinitely many gaps (Lemma 2.1), i.e. open intervals such that the intersection of with the spectrum is reduced to . Moreover, for , there is at least one eigenvalue in each gap (Lemma 2.5). Since, in addition, for sufficiently small, the spectrum of is close to the spectrum of , the existence of guided modes is guaranteed (Theorem 3.1).
2 The spectrum of the limit operator
2.1 Definition of the limit operator
The limit operator is defined on the infinite periodic graph obtained as the limit of as tends to : is made of the vertices connected by the edges
It is -periodic with respect to and -periodic with respect to (see Fig. 3).
For any function defined on , we denote by (resp. ) its value at the vertex (resp. ). The restriction of to the edge (resp. and ) is denoted by (resp. and ).
The definition of also requires the introduction of the function spaces and defined as
and, for any , is the weight coefficient equal to for and otherwise.
The unbounded limit operator in has domain
and is defined by
The functions of are continuous on and quasi-periodic. Moreover, they satisfy the Kirchhoff conditions (8) that enforce the weighted sum of the outward derivatives of to vanish at each vertex (). We point out that the perturbation, which results from a geometrical modification of the domain for the problem (4), is taken into account at the limit by means of the Kirchhoff condition (8) at the vertex (). The formal derivation of the limit model can be found in . It is easily verified that the operator is self-adjoint (for the weighted scalar product associated with (6)), see also . The objective of the following two sections is to study the spectrum of .
2.2 Characterization and properties of the essential spectrum of
By a compact perturbation argument, one can prove that , where is the purely periodic operator corresponding to for . The computation of its spectrum relies on the Floquet-Bloch theory (see ). More precisely, we can prove that if and only if either and or and there exists such that
Based on the previous characterization, we prove that the operator has a countable infinity of gaps that can be separated into three categories (see  for the proof):
The following properties hold :
For any , the operator has infinitely many gaps whose ends tend to infinity.
Let and .
If an interval is a spectral gap of , then, one of the following possibilities holds:
, , and there is a unique .
, and .
, and .
2.3 Computation of the discrete spectrum
Let us now determine the discrete spectrum of . If is an eigenvalue of , then the corresponding eigenfunction satisfies the linear differential equation on each edge of the graph . Solving explicitly this equation (on each edge), taking into account the quasi-periodicity of and the Kirchhoff conditions (8), we can replace the eigenvalue problem with a set of finite differences equations for :
Assume that . is an eigenfunction of if and only if belongs to and satisfies
where we have defined
As well-known, finite difference schemes may be investigated using the discrete Fourier transform
where is an isometry between and . This, together with Lemma 2.2, provides the following characterization for the discrete spectrum of :
Assume that . is an eigenfunction of if and only if the discrete Fourier transform of belongs to and satisfies
Under the assumption of Lemma 2.3, (10) can be written as . It follows that does not belong to if and only if, for any , does not vanish. As a consequence, as soon as , the function is continuous and bounded. Then, the inverse discrete Fourier transform can be applied to (13) to obtain
Writing the previous relation for yields the following criterion of existence of an eigenvalue:
Assume that and that . Then, is an eigenvalue of if and only if
The study of the behavior of the function leads to the existence of at least one eigenvalue in each gap of as soon as , the minimal number of eigenvalues in each gap depending on the type of gaps (cf. Lemma. 2.1-3 for the classification):
For , the operator has no eigenvalue. For , let be a spectral gap of the operator :
If is a gap of type (i), then has at least two eigenvalues and that satisfy (see Lemma. 2.1-3 for the definition of ).
If is a gap of type (ii) or (iii), then has at least one eigenvalue such that .
The sketch of the proof of the previous lemma is the following, a complete proof being available in  (Theorem 5.2.1): First, one can verify that in any gap, which, together with (14) proves that has no eigenvalue for . Then, if is a gap of type (i), one can show that
By continuity of inside the gap, (a) directly results from the intermediate value theorem and (14). If is a gap of type (ii), the intermediate value theorem also permits us to conclude since
A similar argument works for a gap of type (iii).
3 Guided modes for the operator : an asymptotic result
Finally, thanks to the general result  (Theorem 2.13 convergence of the spectrum of toward the spectrum of ), we can prove the following result of existence of eigenvalue for the operator :
Let , be a spectral gap of the operator and be a (simple) eigenvalue of this operator. Then, there exists such that if the operator has an eigenvalue inside a spectral gap . Moreover,
- Ammari, H. and Santosa, F., Guided waves in a photonic bandgap structure with a line defect, SIAM Journal on Applied Mathematics, Vol. 64, Num. 6, p 2018-2033, 2004
- Fliss, Sonia, A Dirichlet-to-Neumann approach for the exact computation of guided modes in photonic crystal waveguides, SIAM Journal on Scientific Computing, Vol. 35, Num. 2, p B438–B461, 2013
- Kuchment P., Quantum graphs: an introduction and a brief survey, Analysis on graphs and its applications,. Ser. Proc. Sympos. Pure Math., Vol. 77, p. 291–312, Amer. Math. Soc., Providence, RI, 2008
- Friedlander, L., On the spectrum of a class of second order periodic elliptic differential operators, Communications in Mathematical Physics, Vol. 229, p. 49–55, 2002
- Kuchment, P. and Ong, B. S. , On guided waves in photonic crystal waveguides, Waves in periodic and random media, Vol. 339, p. 105–115, 2003
- Kuchment, P. and Zeng, H., Convergence of spectra of mesoscopic systems collapsing onto a graph, Journal of Mathematical Analysis and Applications, Vol. 258, Num. 2, p. 671-700, 2001
- Nazarov, S.A., On the spectrum of the Laplace operator on the infinite Dirichlet ladder, St. Petersburg Mathematical Journal, Vol. 23, p. 1023–1045, 2012
- Post O. , Spectral convergence of quasi-one-dimensional spaces, Ann. Henri Poincaré, Vol. 7, Num. 5, p. 933-977, 2006
- Reed, M. and Simon, B., Methods of modern mathematical physics v. I-IV, Academic Press,1972-1978
- Vasilevskaya, E., Open periodic waveguides. Theory and application, PhD thesis, 2016