Abstract
We treat the linear heat equation in a periodic waveguide \(\Pi \subset {{\mathbb {R}}}^d\), with a regular enough boundary, by using the Floquet transform methods. Applying the Floquet transform \({{\textsf{F}}}\) to the equation yields a heat equation with mixed boundary conditions on the periodic cell \(\varpi \) of \(\Pi \), and we analyse the connection between the solutions of the two problems. The considerations involve a description of the spectral projections onto subspaces \({{\mathcal {H}}}_S \subset L^2(\Pi )\) corresponding certain spectral components. We also show that the translated Wannier functions form an orthonormal basis in \({{\mathcal {H}}}_S\).
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
We consider the boundary-initial value problem for the classical linear heat equation with an unknown function u depending on place x and time t,
in a periodic, regular enough (see below) x-domain \(\Pi \subset {{\mathbb {R}}}^d\), \(d=3,4,5\ldots \), which is contained in a cylinder with a bounded cross-section and axis in the \(x_1\)-direction. Here g is the given initial data, on which we assume that
so that there holds \(g \in L^\infty (\Pi ) \cap L^1(\Pi ) \subset L^2(\Pi )\). Moreover, \(\partial _t = \partial / \partial t\) is the time-derivative, \(\Delta = \nabla \cdot \nabla \) and \(\nabla \) is the gradient in the x-variable, \(\partial _\nu \) is the outward normal derivative on the surface of a given x-domain and \(\iota =0\) or 1 so that in (1.3) we have the homogeneous Dirichlet or Neumann condition in the x-domain.
The main aim of this paper is to analyse the unique, classical solution \(u = u(x,t)\) of the problem (1.1)–(1.3) by Floquet transform methods and in particular to relate the spectral representation of the solution with the eigenfunctions of the so called model problem (3.4)–(3.6) (see below), which is a spectral elliptic problem on the bounded periodic cell of \(\Pi \).
We start by recalling some well-known facts on the solution u. First, all partial derivatives up to the first order in t and second order in x exist and are continuous functions of all variables. Moreover, u can be written by using the spectral representation of the elliptic (unbounded, self-adjoint, positive) operator \({{\mathcal {T}}}\), which is defined in a subdomain of \(L^2(\Pi )\) ([5, 40], Theorem VIII.6; [41], Theorem 13.30):
Here \(\sigma = \sigma ({{\mathcal {T}}}) \subset [0, \infty )\) is the spectrum of \({{\mathcal {T}}}\), \( (E_\lambda )_{\lambda \in [0,\infty )} \) is the resolution of the identity (see (3.14) of [45]), or the spectral family, associated with the spectral measure \(P(\Omega )\) of \({{\mathcal {T}}}\), where \(\Omega \subset \sigma \) is measurable (see the spectral theorem, e.g. Theorem 3.6. in [45], or Theorem VIII.6 of [40]), i.e. \(E_\lambda = P((-\infty , \lambda ])\). Formula (1.5) shows, among other things that we have
The elliptic problem related to (1.1)–(1.3) is the Laplace-Dirichlet or Neumann spectral problem (2.6)–(2.7) (see below) which is a standard spectral problem for the operator \({{\mathcal {T}}}\) in \(L^2(\Pi )\). Our assumptions will include that the spectrum \(\sigma \) coincides with the essential \(\sigma _\textrm{ess} = \sigma _\textrm{ess} ({{\mathcal {T}}})\) or continuous spectra of \({{\mathcal {T}}}\) and has band-gap structure (cf. Theorem 3.2), which is a consequence of the Floquet transform theory. The application of the Floquet transform \({{\textsf{F}}}\) to the elliptic problem (2.6)–(2.7) leads to a spectral problem which depends on the Floquet parameter \(\eta \in [-\pi ,\pi ]\), in a bounded periodic cell. The problem has a discrete spectrum consisting of a sequence of eigenvalues \(\Lambda _k(\eta )\), \(k \in {{\mathbb {N}}}= \{1,2,3,\ldots \}\), corresponding to an orthonormal sequence of eigenfunctions \(V_k(\cdot ;\eta )\). The general spectral band formula, see (3.18) below, shows that the spectral bands of \(\sigma \) emerge from the unions of the sets of eigenvalues \(\Lambda _k(\eta )\).
The Floquet transform \({{\textsf{F}}}\) is applied in Sect. 5 also directly to the original initial-boundary value problem (1.1)–(1.3), and the resulting description of the solution u is given in formula (5.11) of Proposition 5.1. In Sect. 6 we describe the heat semigroup and the spectral measure and projections of the operator \({{\mathcal {T}}}\) in terms of the eigenfunctions \(V_k\) of the model problem, see Theorems 6.3 and 6.5.
In Sect. 6 we consider a spectral subspace \({{\mathcal {H}}}_S\) corresponding to a disjoint component of \(\sigma \). Since the related spectral bands do not overlap in this case, the presentation of the spectral projection corresponding to \({{\mathcal {H}}}_S\) in Theorem 6.5 can be simplified, see Proposition 7.1. Moreover, in Theorem 7.2 an orthonormal basis of \({{\mathcal {H}}}_S\): it is formed by certain inverse Floquet transforms \({{\textsf{F}}}^{-1}V_k\) of the eigenfunctions and their translates (these are analogous to the so called Wannier functions, see the next paragraph). The consequences of these considerations to long time decay estimates of the parabolic problem are straightforward, as mentioned at the end of the paper. On the other hand, the exponentially decaying factor \(e^{-\lambda t}\), cf. (1.5), which is characteristic for the solutions of the heat equation, is essential for our calculations and treatment of the spectral projections.
As for literature, we mention that band-gap spectra in geometric situations similar to ours have been considered e.g. in the papers [3, 4, 36, 37]. However, the authors do not know papers dealing with the consequences of these studies to parabolic initial-boundary value problems. (Floquet-theory in the setting of parabolic equations has been studied e.g. in [19, 23, 24, 30] and others, but these works concern periodicity with respect to the time-variable, which is unrelated to periodicity in the x-domain studied in the present work.) We also mention the following references as an incomplete list of other literature on the subject of spectral gaps in elliptic boundary value problems of various types: [2, 7, 8, 10,11,12,13, 15, 17, 18, 31,32,33,34, 47, 48]. For presentations of the theory of Floquet transform, see [24, 26, 35]. In the setting of the time-independent Schrödinger equation with periodic potential, the functions corresponding to \({{\textsf{F}}}^{-1} V_k\) are known as Wannier functions. For example the papers [22, 25, 38, 39] contain studies of the question, whether the spectral subspace corresponding to a spectral band of the Schrödinger equation has an orthonormal basis spanned by translates of exponentially decreasing Wannier functions. See also [25] for more references especially to the literature on molecular chemistry.
As for the notation used in this article, B(a, r) or \(B_d(a,r)\) denote a Euclidean ball in \({{\mathbb {R}}}^d\) with center a and radius \(r > 0\). For a measurable set \(A \subset {{\mathbb {R}}}^d\), |A| means the volume, and if \(x \in {{\mathbb {R}}}^d\), then dist (x, A) denotes the Euclidean distance of x from the set A. By C, \(C'\), \(C_1\) etc. (respectively, \(C_k\), \(C_\eta '\) etc.) we denote generic positive constants which do not depend on the functions or variables in the calculations in question (resp. depend only on k, \(\eta \)) and the exact value of which may change from place to place. We write \({{\mathbb {N}}}= \{1,2,3,\ldots \}\), \({{\mathbb {Z}}}= \{0, \pm 1, \pm 2,\}\), \({{\mathbb {R}}}_+ =(0,\infty )\). If \(x \in {{\mathbb {R}}}\), [x] stands for the largest integer not larger than x. Given \(x=(x_1,x') \in {{\mathbb {R}}}\times {{\mathbb {R}}}^{d-1}\), we write \(x^{\textsf{tr}}= ( x_1-[x_1], x')\). For \(j \in {{\mathbb {N}}}\) we denote \(\partial _j:= \partial / \partial {x_j}\). Moreover, given a domain \(\Omega \subset {{\mathbb {R}}}^d\), we denote by \(\Vert u \Vert _\Omega \) and \((u,v)_{\Omega }\) the norm and scalar product of the space \(L^2 (\Omega ) \) with the complex scalar field. The notation for the other \(L^p(\Omega )\)-spaces is standard. In general, given a Banach space X, its norm is denoted by \(\Vert \cdot \Vert _X\). The standard Sobolev-Hilbert space of order \(k \in {{\mathbb {N}}}\) is denoted by \(H^k(\Omega )\), and \(H_0^k(\Omega )\) is its closed subspace spanned by the compactly supported infinitely smooth test functions, i.e., elements of the space \(C_0^\infty (\Omega )\). Moreover, \(C_B(\Omega )\) (respectively, \(C_B^1(\Omega )\)) denotes the space of bounded, continuous functions on \(\Omega \) (resp. continuously differentiable functions on \(\Omega \) with bounded partial derivatives of the first order). The Hilbert space of square summable sequences of complex numbers is denoted by \(\ell ^2\), or by \(\ell ^2({{\mathbb {N}}})\) or \( \ell ^2({{\mathbb {Z}}})\), if it is necessary to indicate the index set.
We will need to pay attention to some details of Banach-space valued integrals. For the basic theory of vector valued \(L^2\)-spaces and the Bochner-integral we refer to [20] (Fig. 1).
2 Preliminaries
Let us proceed with the details of the geometric assumptions. By \(\varpi \subset {{\mathbb {R}}}^d\) we denote the primary periodic cell
which is a bounded domain in \({{\mathbb {R}}}^d\) with boundary \(\partial \varpi \). We write \(T^\pm = \{ x \in {{\mathbb {R}}}^d \,: \, x_1 = \pm 1/2 \}\) and \(P^{\pm } = (\pm 1/2, 0, \ldots ,0 ) \in T^\pm \) and assume that there holds
for some open set \(\omega \subset {{\mathbb {R}}}^{d-1}\) containing the origin of \((0, \ldots , 0)\) of \({{\mathbb {R}}}^{d-1}\) (hence, \(P^\pm \subset \partial \varpi \)). Furthermore, we denote the translates of \(\varpi \) by
and define the main periodic domain (waveguide) \(\Pi \) as the interior of the union
Assumption 1
We require that the boundary \(\partial \Pi \) is of smoothness \(C^\kappa \), where \(\kappa \in {{\mathbb {N}}}\) is the smallest integer such that
We will consider the following spectral Laplace-Dirichlet or Neumann problem for an unknown function v on the domain \(\Pi \), related with the initial-boundary value problem (1.1)–(1.3):
The variational formulation of this problem, see [28, 29], reads in the Neumann case \(\iota =1\) as
and in the Dirichlet case \(\iota =0\) the space \(H^1(\Pi )\) is replaced by \(H_0^1(\Pi )\). Also, the negative Laplacian \(-\Delta \) defines an unbounded self–adjoint positive operator \({{\mathcal {T}}}: {{\mathcal {D}}}\rightarrow L^2 (\Pi ) \) with the domain of definition \({{\mathcal {D}}}= H_0^1(\Pi ) \cap H^2(\Pi ) \) in the Dirichlet case or \({{\mathcal {D}}}= H_\bullet ^1(\Pi ) \cap H^2(\Pi ) \) in the Neumann case, where
By positivity, the spectrum \(\sigma = \sigma ({{\mathcal {T}}})\) (the same as the spectrum of the problem (2.6)–(2.7)) lies in \(\overline{{{\mathbb {R}}}}_+ = [0, \infty )\), and since the embedding \(H^1 ( \Pi ) \subset L^2(\Pi ) \) is not compact in the unbounded domain \(\Pi \), the essential spectrum \(\sigma _\textrm{ess}= \sigma _\textrm{ess} ({{\mathcal {T}}})\) is not empty (cf. [5, Th.10.1.5]). Actually, it is a well known consequence of the Floquet transform theory, that the spectrum has the so called band-gap structure. We will recall this in Sect. 3, and we later use this to derive consequences on the long-time asymptotics of the solutions of the problem (1.1)–(1.3).
Actually, we pose the following
Assumption 2
The spectrum \(\sigma \) is absolutely continuous.
This assumption in particular implies that \(\sigma ({{\mathcal {T}}}) = \sigma _\textrm{ess} ({{\mathcal {T}}})\) and that, for the problem (2.6)–(2.7), there does not exist eigenvalues of finite or infinite multipliticity (with eigenfunctions belonging to \(L^2(\Pi )\). It is generally believed that indeed \(\sigma \) is absolutely continuous in the situation of this paper, although this has not yet been proved. For example, by the result of [14], the absolute continuity holds, if the boundary \(\partial \Pi \) is smooth and in addition \(\varpi \) and thus \(\Pi \) are symmetric under the reflection \(x_1 \mapsto - x_1\). Partial positive results on the absolute continuity can also be found e.g. [9, 24, 27, 42,43,44]; see also the references on p. 383 of [26].
3 Floquet transform and band–gap spectrum of the elliptic problem
In this section we recall the basics of the Floquet transform theory and apply it in order to turn problem (2.6)–(2.7) into a model problem in the periodic cell \(\varpi \). The well-known Theorem 3.2 recalls how the band-gap structure of the spectrum of \({{\mathcal {T}}}\) emerges from the spectra of the model problems.
We recall the definition of the Floquet transform (see for example [16, 24, 26, 35]):
where \(v \in L^2(\Pi )\), \(x =(x_1,x') \in \Pi \) on the left, \(\eta \in [-\pi , \pi ]\), and \(x \in \varpi \) on the right. This operator is a linear isometry from \(L^2 ( \Pi ) \) onto \(L^2 ( -\pi ,\pi ; L^2 (\varpi ) )\), which also preserves the inner product; the inner product of \(L^2 ( -\pi ,\pi ; L^2 (\varpi ) ) \ni f,g \) is defined by
Moreover, it is a (Banach-space) isomorphism from \(H^1 ( \Pi ) \) onto \(L^2 ( -\pi ,\pi ; H^1_{\eta } (\varpi ) )\) (see e.g. [26, Th. 4.2, Th. 4.8], [35, § 3.4], [31, Cor. 3.4.3]). Here \(L^2 ( -\pi ,\pi ; L^2 (\varpi ) )\) consists of (complex) \(L^2 (\varpi ) \)–valued functions on the interval \([-\pi ,\pi ]\), which are \(L^2\)-integrable in the sense of Bochner, and the space \(L^2 ( -\pi ,\pi ; H^1_{\eta } (\varpi ) )\) is the closed subspace of \(L^2 ( -\pi ,\pi ; H^1 (\varpi ) )\) consisting of such vector valued functions f that \(f(\cdot ; \eta ) \in H_\eta ^1 (\varpi )\) for almost all \(\eta \). Moreover, \(H_\eta ^1 (\varpi )\) is the subspace of \(H^1 (\varpi )\) consisting of functions that satisfy, in the sense of traces, the quasiperiodic boundary condition which is the first equation in (3.6), below. The inverse transform has the formula
where \( x\in \Pi \), \( g \in L^2 ( -\pi ,\pi ; L^2 (\varpi ) )\), and \(x^{\textsf{tr}}= ( x_1-[x_1], x')\) so that \(x^{\textsf{tr}}\in \varpi \) for all \(x \in \Pi \).
Using the Floquet transform, problem (2.6)–(2.7) turns into the following model spectral problem in the periodic cell \(\varpi \), (2.1),
where \(\Lambda = \Lambda (\eta ) \) is just a new notation for the spectral parameter. The variational formulation of the problem (3.4)–(3.6) in the case \(\iota =1\) is to find \(V = V(\cdot ; \eta ) \in H_{\eta }^1 (\varpi )\), \(V \not = 0\), and \(\Lambda \) such that
for all \(U \in H_{\eta }^1 (\varpi )\); the inner product is that of \(L^2 (\varpi ) \). The case \(\iota =0\) is similar, except that in (3.7) the space \(H_\eta ^1(\varpi )\) is replaced by its closed subspace \(H_{0,\eta }^1(\varpi )\), which consists of functions vanishing on the boundary \(\partial \varpi \cap \partial \Pi \). In a standard way we also define, for every \(\eta \in [-\pi ,\pi ]\), an unbounded, self-adjoint, positive operator \(S= S(\eta ): {{\mathcal {D}}}_\eta \rightarrow L^2(\varpi )\) such that problem (3.4)–(3.6) is equivalent to the equation \(S(\eta ) V = \Lambda V\). We will not need to specify the domains \({{\mathcal {D}}}_\eta \subset L^2(\varpi )\) precisely, but in the Neumann case \(\iota =1\) we will need the sesquilinear, lower semibounded form associated with \(S(\eta )\), which is the mapping
see e.g. [5], Section 10.1–2. Since the domain \(\varpi \) is bounded, the embedding \(H_{\eta }^1 (\varpi ) \subset L^2(\varpi )\) is compact and the spectra of the problem in both cases \(\iota =0,1\) are discrete, see [5], Th. 10.1.5. The numbers
form the sequence of eigenvalues (counting multiplicities) for the model problem, and the corresponding eigenvectors are denoted by \(V_k = V_k(\cdot ; \eta ) \in H^1_{\eta } (\varpi )\), \(k \in {{\mathbb {N}}}\). We require the orthogonality and normalization conditions
so that for any fixed \(\eta \in [-\pi ,\pi ]\) the functions \(V_k(\cdot ; \eta )\), \(k \in {{\mathbb {N}}}\), form an orthonormal basis of \(L^2(\varpi )\).
Remark 3.1
Let us derive some Sobolev- and sup-norm estimates for the eigenfunctions \(V_k\). Let \(B \subset (-1/2, 1/2) \times {{\mathbb {R}}}^{d-1} \subset {{\mathbb {R}}}^d\) be a ball such that \(B \cap \varpi \not = \emptyset \) and \(B \cap T^\pm = \emptyset \). Then, an eigenfunction \(V_k\) satisfies
Also, according to Assumption 1, the boundary \(\Gamma \) is of smoothness \(C^\kappa \), \(\kappa > d/2 + 3/2\), hence, we can apply the local elliptic estimate in Theorem 15.2 in [1] to get
for all \(\ell = 2,3, \ldots , \kappa \). Thus, (3.10) and (3.13) directly yield a bound for the \(H^2\)-norm of \(V_k\), which can be used a second time with \(\ell =4\) to obtain an estimate for the \(H^4\)-norm of \(V_k\). Repeating this procedure, if necessary, one finally obtains
for some integer \(d^*\) depending on d (to be precise, \(d^* = [d/4] +1\)). Then, the Sobolev embedding \(H^\kappa (\Pi \cap B) \hookrightarrow C_B^1 (\Pi \cap B)\) (see 2.5) gives us
Note that the constants C both in (3.14) and (3.15) can be chosen independently of \(\eta \) or k.
The estimate (3.15) can also be shown for any ball \({\widetilde{B}} \subset (0,1) \times {{\mathbb {R}}}^{d-1} \subset {{\mathbb {R}}}^d\) with radius, say, at most 1/4 and \({\widetilde{B}} \cap \varpi \not = \emptyset \) and \({\widetilde{B}} \cap T^+ \not = \emptyset \), by using the quasiperiodicity. Namely, we define the function
and observe (using (3.6) among other things) that this is an eigenfunction with eigenvalue \(\Lambda _k(\eta )\) of a problem similar to (3.4)–(3.6), but in the cell \({\widetilde{\varpi }} = \Pi \cap (0,1) \times {{\mathbb {R}}}^{d-1} \) instead of \(\varpi \). The previous argument applies and yields the same Hölder estimate for \({\widetilde{V}}_k\) and thus for \(V_k\) in the ball \({\widetilde{B}}\).
The balls with \(B \cap T^- \not = \emptyset \) are treated in the same way. Since \(\varpi \) can be covered by finitely many balls B and \({\widetilde{B}}\), we can combine these estimates and deduce that every eigenfunction \(V_k(\cdot ; \eta )\) belongs to \(H^\kappa (\varpi )\) and \(C^1_B (\varpi )\), and
where \(d^* \) is as in (3.14) and the constant C in particular does not depend k or \(\eta \).
For the following fact, see [35, Th. 3.4.6], [31, Th. 2.1], and [24, 26] in the case of the operators in the entire Euclidean space.
Theorem 3.2
A number \(\lambda \) belongs to the resolvent set or to the discrete spectrum of \({{\mathcal {T}}}\), if and only if it does not coincide with \(\Lambda _k(\eta )\) for any \(\eta \in [-\pi ,\pi ]\) and k.
Hence, \(\sigma _\textrm{ess} ({{\mathcal {T}}})\) and thus also \(\sigma ({{\mathcal {T}}})\) (see Assumption 2) get the band–gap structure, namely
where the kth spectral band is defined by
It is known that the functions \({{\mathbb {R}}}\ni \eta \mapsto \Lambda _k(\eta )\) are continuous and even piecewise real analytic (see [21, Ch.7]), which implies that every \(\Upsilon _k\) is a closed interval, and \(2 \pi \)–periodic: replacing \(\eta \) by \(\eta + 2 \pi \ell \), the function \(e^{- i 2 \pi \ell z} V_k\) is again an eigenfunction with eigenvalue \(\Lambda _k(\eta )\), for all \(\ell \in {{\mathbb {Z}}}\). We denote the endpoints of the spectral bands \(\Upsilon _{k}\), \(k \in {{\mathbb {N}}}\), by
Remark 3.3
Our Assumption 2 implies that none of the functions \( \eta \mapsto \Lambda _k(\eta )\) is constant (even on a subset of \([-\pi ,\pi ]\) with positive measure). Namely, if \(\Lambda _k(\eta ) = r\) for some number \(r \in {{\mathbb {R}}}\) on a set of positive measure, then r is an eigenvalue of \({{\mathcal {T}}}\), which we are assuming not to happen. See Section 6.3. of [26].
4 Lower bounds for the eigenvalues of the model problem.
The aim of this section is to prove a lower bound for the eigenvalues (3.9) of the model problem. For the purposes of this paper, there is no need to approach optimal results, namely, quite rough estimates will suffice to guarantee that the series (5.8) converges well for all \(t > 0\). The result we need is contained in the following lemma. In the Dirichlet case, the lower bound can be obtained by comparing the eigenvalues with those of the corresponding problem on a cube and using the max-min-principle. In the Neumann case it is possible to give a simpler proof with the help of associated semibounded forms and a comparison with the eigenvalues of the pure Neumann problem in \(\varpi \).
Lemma 4.1
There is a constant \(\gamma >0\) such that for all \(k \in {{\mathbb {N}}}\), \(\eta \in [-\pi ,\pi ]\), we have the lower bound
for the eigenvalues (3.9).
Proof. We first consider the Dirichlet case \(\iota =0\) and choose a large enough number \(M> 0\) such that \(\varpi \subset (-1/ 2, 1/2) \times (-M,M) ^{d-1} =: Q_M \). We fix k and \(\eta \) and consider the following eigenvalue problem,
The variational formulation of this problem is composed in the same way as for the problem (3.4)–(3.6) and it amounts to finding \(V\in H_{0,\eta }^1(Q_M)\) and \(\Lambda \) such that
for all \(U \in H_{0,\eta }^1(Q_M)\); this Sobolev-type space is defined in the same way as in Sect. 3. Also, the solution of the problem consists of an increasing, positive, unbounded sequence of eigenvalues \({\widetilde{\Lambda }}_k = {\widetilde{\Lambda }}_k ( \eta )\) together with eigenfunctions, which are of the form
where \(n = (n_1, \ldots ,n_d)\) is running over all multi-indices with \(n_1 \ge 0\) and \(n_j> 0\) for \(j=2,\ldots ,d\). The function (4.6) corresponds to the eigenvalue \( (\eta + 2 \pi n_1)^2 + \pi ^2 (4\,M)^{-2} \sum _{j=2}^d n_j^2 \), hence, given \(k \in {{\mathbb {N}}}\), there are less than \(C_M k^{d/2}\) such eigenvalues belonging to the interval [0, k], where \(C_M > 0\) is a large enough constant. We thus obtain the lower bound
for the kth eigenvalue of the problem (4.2)–(4.4).
On the other hand, the eigenvalue \({\widetilde{\Lambda }}_k (\eta )\) is obtained from the max-min-principle ([5, 40], Theorem 10.2.2]) as
where in the supremum, \(H^{(k)}\) is running over all \(k-\)1-codimensional subspaces of \(H_{0,\eta }^1(Q_M)\). We now fix \(k \in {{\mathbb {N}}}\) and extend all eigenfunctions \(V_j(\cdot ;\eta )\), \(j =1, \ldots ,k\), as 0 to the domain \(Q_M \) and denote the extensions by \({\widetilde{V}}_j( \cdot ;\eta )\). These functions belong to \(H_{0,\eta }^1(Q_M)\) and still form an orthonormal sequence in \(L^2(Q_M)\), by (3.10), and also an orthogonal sequence in \(H_{0,\eta }^1(Q_M)\), by (3.7). Thus, given an arbitrary \(k-1\)-codimensional subspace \(H^{(k)}\) of \(H_{0,\eta }^1(\varpi )\), there are coefficients \(\varrho _j \in {{\mathbb {C}}}\) such that
Given an arbitrary \(\varepsilon >0\) we choose a \(k-1\)-codimensional subspace \(\textbf{H}\) of \(H_{0,\eta }^1(Q_M)\) such that
and choose the coefficients \(\varrho _j\) such that (4.9) holds for \(\textbf{H}\). Noting that \((V_j, V_\ell )_\varpi = \Lambda _j \delta _{j, \ell }\) by (3.7), (3.10), we obtain from (4.8), (3.9), (4.10)
This and (4.7) imply (4.1) in the Dirichlet case for some constant \(\gamma > 0\), since \(\varepsilon > 0\) was arbitrary.
In the Neumann case \(\iota =1\) we consider the pure Neumann problem in the periodic cell,
Again, this problem has an unbounded, non-negative sequence of eigenvalues \(0= {\widehat{\Lambda }}_1 < {\widehat{\Lambda }}_2 \le {\widehat{\Lambda }}_3 \le \cdots \), for which there holds
for example just by the classical Weyl estimate [46]. Moreover, the standard problem operator of (4.12), denoted here by \(S_N\), is associated with the sesquilinear, lower semibounded form
Now, comparing with (3.8), there holds \(S(\eta ) > S_N\) in the sense of the discussion preceding Theorem 10.2.4. of [5]. This reference and (4.13) imply that we have
This proves the lemma.
5 Solution of the initial-boundary value problem via Floquet transform
Applying the Floquet transform formally with respect to the \(x_1\)-variable to the problem (1.1)–(1.3) leads to the following system for the function \(U = {{\textsf{F}}}u\):
where g is the Cauchy data in (1.2) and \(\iota =0\) (Dirichlet case) or 1 (Neumann case). Note that the assumption (1.4) and the very definition of the Floquet transform, (3.1) yield
For a fixed \(\eta \), (5.1)–(5.4) is an initial-boundary value problem for the heat equation in a bounded x-domain, although the mixed boundary conditions (5.3)–(5.4) for the unknown function U are somewhat unconventional. Yet, the corresponding spectral elliptic problem is nothing but (3.4)–(3.6), associated with the positive, self-adjoint operator \({{\mathcal {S}}}= {{\mathcal {S}}}( \eta )\) with domain \(D_\eta \subset L^2(\varpi )\). The operators \({{\mathcal {S}}}\) are in particular maximal monotone, hence the Hille-Yosida-theorem ( [6], Theorem 7.4) implies that problem (5.1)–(5.4) has a unique solution
and we moreover have
for every \(\eta \).
We will show shortly that the following series coincides with (5.6):
Here,
are the coordinates of the initial data \(G(\cdot ;\eta )\) in the orthonormal basis \(\big (V_k(\cdot ;\eta ) \big )_{k=1}^\infty \) of \(L^2(\varpi )\). Thus, for all \(\eta \) and \(t \ge 0\), (5.8) converges in \(L^2(\varpi )\), and the sequence \(\big ( G_k(\eta ) \big )_{k=1}^\infty \) belongs to \(\ell ^2\). Note that due to (5.5) and the normalization (3.10), there exists a constant \(C> 0\) such that
for all \(k \in {{\mathbb {N}}}\) and \(\eta \in [-\pi ,\pi ]\).
Moreover, if \(t > 0\), then Lemma 4.1 implies that the series (5.8) together with its termwise t-derivative converge uniformly with respect to t and x. Thus, termwise differentiation yields the function \(\partial _t U\). Similarly, termwise differentiation with \(\Delta \) yields \(-\Delta U\), since the terms of the series contain eigenfunctions of \(\Delta \) and the resulting coefficients \(\Lambda _k(\eta ) e^{- \Lambda _k(\eta ) t}\) are still small enough. We obtain that the series (5.8) is a solution of (5.1).
As for the boundary conditions (5.3)–(5.4), we obtain by (3.17) and Lemma 4.1 that the series (5.8) converges absolutely also in \(H^\kappa (\varpi )\) for every \(t > 0\), hence, these boundary conditions are satisfied by the sum of the series, since this is true for all terms.
Finally, (5.2) holds since due to the orthonormality of the functions \(V_k\) we have for U of the series (5.8)
as a limit in the space \(L^2 (\varpi )\), for every \(\eta \). Since the solution of (5.1)–(5.4) is unique, (5.6) and (5.8) coincide.
We next verify that \({{\textsf{F}}}^{-1} U \) is the unique solution of the original initial-boundary problem (1.1)–(1.3). Indeed, by (5.6), (5.7) and (5.5), the function \(\partial _t U(\cdot , t;\eta )\) (respectively, due to (5.1), also \(\Delta U(\cdot ,t;\eta )\) ) belong to \(L^2(\varpi )\) for every \(t \ge 0\) and \(\eta \), with an \(\eta \)-uniform bound for its norm in \(L^2(\varpi )\). By the Banach-space valued Lebesgue dominated convergence theorem, e.g. Proposition 1.2.5 in [20], and just the integral definition of the inverse Floquet transform, this implies that the derivative \(\partial _t\) (respectively, \(\Delta \)) and \({{\textsf{F}}}^{-1}\) can be commuted. We obtain that \({{\textsf{F}}}^{-1} U \) satisfies (1.1). As a consequence of (5.3), it also satisfies the boundary condition (1.2). Finally, by (5.6), \(U(,t;\eta ) \rightarrow G(\cdot ;\eta )\) in \(L^2(\varpi )\) as \(t \rightarrow 0^+\) for every \(\eta \), hence \({{\textsf{F}}}^{-1} U (\cdot ,t) \rightarrow g\) again by the vector valued Lebesgue theorem. Hence, also the initial condition (1.3) holds, which completes the proof of the claim.
Proposition 5.1
The unique solution u of the problem (1.1)–(1.3) can be written for \(t > 0\) as the series
which converges absolutely in \(L^2(\Pi )\). Here, \(G_k\) is as in (5.9) and \(V_k\) is the kth eigenfunction in (3.10).
Proof. As a consequence of (3.10) and Lemma 4.1, the series (5.8) converges in \(L^2(\varpi )\) absolutely and uniformly with respect to \(\eta \), hence it converges absolutely in \(L^2(-\pi , \pi ;L^2(\varpi ))\). The absolute convergence of the series (5.11) in \(L^2(\Pi )\) thus follows from the continuity of the operator \({{\textsf{F}}}^{-1}\) in the relevant Hilbert spaces. Hence, (5.11) coincides with \({{\textsf{F}}}^{-1} U\), which by the argument above is the unique solution of (1.1)–(1.3). \(\Box \)
6 Spectral projections: general formulas
In this section we present the main general results concerning the presentation of the spectral projections \(P(\Omega ) \) of the operator \({{\mathcal {T}}}\), see (1.5) in terms of the eigenfunctions \(V_k\) of the \(\eta \)-dependent model problem. The main technical issues consist of handling the overlapping of the spectral bands and clarifying the relations of the Floquet and spectral parameters especially as regards to integration by substitution.
We recall that the solution of the problem (1.1)–(1.3) can be written as
where \((E_\lambda )_{\lambda \in [0,\infty )}\) is the resolution of the identity (see (3.14) of [45]) associated with the spectral measure \(P(\Omega \)) of the operator \({{\mathcal {T}}}\) (Theorem 3.6. in [45], or Theorem VIII.6 of [40]) so that there holds \(E_\lambda = P( (-\infty , \lambda ])\).
Given \(k \in {{\mathbb {N}}}\), we decompose the interval \([-\pi , \pi ]\) into finitely many subintervals \(J_{k,q}\), \(q = 1, \ldots , q_k\), such that the mapping \(\Lambda _k: \eta \mapsto \Lambda _k(\eta )\) is monotone, i.e. a bijection, on each of them and denote the inverse mappings \(\eta _{k,q}: I_{k,q} \rightarrow J_{k,q}\), where \(I_{k,q} = \Lambda _k(J_{k,q}) \subset \Upsilon _k\). This is possible because of Remark 3.3. Obviously, there holds the change of variables formula
where f is a scalar valued Lebesgue-integrable function or more generally a Banach-space valued Bochner-integrable function on \([-\pi ,\pi ]\). Note that here indeed k is arbitrary.
For every \(\lambda \in \sigma \) we denote by \({{\mathbb {K}}}_\lambda \subset {{\mathbb {N}}}\times {{\mathbb {N}}}\) the (finite) set of all pairs (k, q) such that \(\lambda \in I_{k,q} \subset \Upsilon _k\). The spectrum \(\sigma \) can be decomposed into a disjoint union of countably many subintervals such that the set-valued mapping \(\lambda \mapsto {{\mathbb {K}}}_\lambda \) is constant on these subintervals. If \(m \in {{\mathbb {N}}}\) is arbitrary, then Lemma 4.1 implies that there are at most m indices k such that \((k,q) \in {{\mathbb {K}}}_\lambda \) for some \(\lambda \le C_\varpi m^{1/(dd^*)}\).
Formula (6.3) in the next definition looks complicated, but note that it simplifies a lot by applying the change of variables (6.2); see e.g. the second identity in (6.7).
Definition 6.1
For \(\lambda \in \sigma \) and all pairs \((k,q) \in {{\mathbb {K}}}_\lambda \) we denote
where \( x \in \Pi \), \(V_k\) is the eigenfunction (3.7) of the model problem and \(G_k\) is the inner product (5.9) for the function \(g \in L^2(\Pi )\). For all \(M> 0\) we also denote by \(F_M = F_M(\cdot ; \lambda ,k,q)\) the restriction of F to the set \(\Pi ^M:= \{x \in \Pi \,: \, |x| \le M \}\).
Note that usually F does not belong to \(L^2(\Pi )\) or \(L^1(\Pi )\), since its modulus is a periodic function. This fact complicates some of following calculations and the formulation of the results. Consequently, we need to formulate the following technical observation.
Lemma 6.2
For all \(M> 0\) the function
is Bochner-integrable on \(\sigma \cap I\) (as a function \(\sigma \cap I \rightarrow L^2(\Pi ^M)\)), for every bounded interval I. Moreover, the limit
converges in \(L^2(\Pi )\).
Before proceeding with the proof, we remark that the sets \({{\mathbb {K}}}_\lambda \) are finite and depend piecewise constantly on \(\lambda \in I\), hence, there are altogether only finitely many terms in the sums in (6.5). Thus, due to (5.10), (3.17) and the boundedness of the derivatives \(\eta _{k,q}'\), both \(\lambda \)-integrals in (6.5) converge for every fixed \(x \in \Pi \).
Proof. Let \(M> 0\) and a bounded \(\lambda \)-interval I be given. Denoting by \({{\mathbb {K}}}(I)\) the finite set of all k such that the pair (k, q) appears in (6.4) for some \(q= 1, \ldots , q_k\), we have by (3.10), (5.10), (6.3) and the definition of the notation \(x^{\textsf{tr}}\)
This proves the Bochner-integrability.
Moreover, for a given \(M> 0\) and for almost all \(x \in \Pi ^M\) (x fixed for a moment), the change of the integration variable (6.2) yields
where \(N_k\) is a subset of the interval \([-\pi ,\pi ]\) (as \(\Upsilon _k\) may not be completely contained in I) and \(\chi _k\) is its characteristic function. If \(x \in \Pi {\setminus } \Pi ^M\), then the expression in (6.7), i.e., \(G_M(x)\), is null, due the very definition of \(F_M\).
Now, the function
belongs to \(L^2(-\pi ,\pi ; L^2(\varpi ))\), due to (3.10), (5.10), hence, \(G:= {{\textsf{F}}}^{-1} H\) belongs to \(L^2(\Pi )\), and by (6.7) we have \(G_M(x) = G(x) \) for \(x \in \Pi ^M\). Since \(G_M(x)= 0\) for x outside \(\Pi ^M\), we have \(G_M \rightarrow G\) in \(L^2(\Pi )\). Finally, the same calculation as in (6.7) shows that G coincides with the right-hand side of (6.5). \(\Box \)
In principle, the proof of the next theorem follows by changing the integration variables, but it is technically complicated by the problems mentioned in the beginning of this section.
Theorem 6.3
If \(g \in L^2(\Pi )\) is a function satisfying (1.4) and \(t>0\), then the limit
converges in the space \(L^2(\Pi )\) and it coincides with
Proof. Fix \(t > 0\), and then let \(\varepsilon > 0\) be arbitrary. Fix \(k_0 = k_0(t) \in {{\mathbb {N}}}\) such that
where the constant \(\gamma \) is as in Lemma 4.1, and then, let \(\lambda _0\) be any number such that
where \(M_{k_0}\) is the upper end of the spectral band \(\Upsilon _{k_0}\), see (3.20). Thus, for every \(k \le k_0\), all bands \(\Upsilon _k \) and intervals \( I_{k,q} \subset \Upsilon _k \) with \(k \le k_0\), \(1=1, \ldots ,q_k\), are contained in \([0,\lambda _0]\).
Also, for \(k > k_0\) it may still happen that \([0,\lambda _0] \cap \Upsilon _k \not = \emptyset \), which requires the next considerations. First, we denote by \(K(\lambda _0)\) the set of those indices \(k > k_0\) such that \(\Upsilon _k \cap [0, \lambda _0] \not = \emptyset \). Moreover, for all \(\lambda \le \lambda _0\), we denote
where \({{\mathbb {K}}}_\lambda \) was defined before Definition 6.1; consequently, given \(\lambda \le \lambda _0\), the set \(R_\lambda \) consists of all indices (k, q) such that \(\lambda \in I_{k,q} \subset \Upsilon _k\) and \(k > k_0\). Note that again, the set valued mapping \(\lambda \mapsto R_\lambda \) is piecewise constant. Finally, for every \(k \in K(\lambda _0)\) and \(\lambda \in [0, \lambda _0]\) we also denote
Note that also the set valued mapping \(\lambda \rightarrow Q_{k,\lambda }\) is piecewise constant. The lower bound (4.1) implies that there exists a constant \(C>0\) (depending on \(\varpi \) only) such that for
we have \(\Upsilon _k \cap [0,\lambda _0] = \emptyset \), hence, \(k < C \lambda _0^{d/2 }\) holds for every \(k \in K(\lambda _0)\) with the same constant C as in (6.14).
By applying formula (6.2) and the definition of the sets \({{\mathbb {K}}}_\lambda \) and \(R_\lambda \) to all \(\eta \)-integrals, we get
Next, we derive an estimate for
For a given \(k \in K(\lambda _0)\), a change of the integration variable yields (cf. 6.8, 6.3)
where \(N_k\) is a subset of the interval \([-\pi ,\pi ]\) (as \(\Upsilon _k\) is not completely contained in \([0,\lambda _0]\), see above). Using (3.10), (5.10), (4.1) we can estimate
Due to the isometric property of the Floquet transform we obtain from (6.17), (6.18) and (6.10)
since \(k > k_0\) for \(k \in K(\lambda _0)\).
Taking into account (6.15), (6.19) yields
for all large enough \(k_0\) and \(\lambda _0\) as in (6.11), hence, in view of the absolute convergence of the series (5.11) in \(L^2(\Pi )\), the claims (6.8) and (6.9) are proved. \(\Box \)
Remark 6.4
In the proof of the next theorem we will need the trivial observation that the “scalar” versions of Lemma 6.2 and Theorem 6.3 also hold: given an arbitrary function \(h \in L^2(\Pi )\), the function
is Lebesgue-integrable for every \(M > 0\), and formulas (6.5), (6.8) holds, when the integrands are replaced on both sides with their \(L^2(\Pi )\)-inner products with h.
The following is a consequence of the previous theorem. Recall that \(P(\Omega )\) denotes the spectral measure of the operator \({{\mathcal {T}}}\). Obviously, formula (6.22) determines the spectral measure in terms of the eigenfunctions of the model problem (5.1)–(5.4).
Theorem 6.5
For every set \(I = [a,b] \cap \sigma \subset {{\mathbb {R}}}_0^+\) with \(0 \le a < b\) there holds the formula
where P is the spectral measure of the operator \({{\mathcal {T}}}\) in (1.5), \(g \in L^2(\Pi )\), F is the expression (6.3) and the integral is defined by (6.5).
Proof. Functions g satisfying (1.4) form a dense subspace in \(L^2(\Pi )\). Moreover, in addition to the orthogonal projection P(I), also the sum on the right hand side of (6.22) defines a bounded operator \(L^2(\Pi ) \rightarrow L^2(\Pi )\),
To see this, the right hand side of (6.22) or (6.23) equals, by using the substitution (6.2) as in (6.17), the expression
where \(\chi _k\) is the characteristic function of some subset of the interval \([-\pi ,\pi ]\) and the number of the terms is bounded by a constant not depending on g. Now, for all g in the unit ball B of \(L^2(\Pi )\), the definition (5.9) implies that \(|G_k(\eta )| \le 1\), hence, the norm of \(\chi _k G_k V_k\) in \(L^2(-\pi ,\pi ;L^2(\varpi ))\) is bounded by a constant independent of \(g \in B\). Using the isometric property of \({{\textsf{F}}}^{-1}\) shows that the \(L^2(\Pi )\)-norm of the right hand side of (6.24) is bounded by a constant independent of \(g \in B\).
Hence, we may assume that (1.4) holds for g and thus the result of Theorem 6.3 can be used. Let \((e_n)_{n=1}^\infty \) be an orthonormal basis of the space \(L^2(\Pi )\). Taking the inner product of (6.8), which is the same as (6.9), with any vector \(e_n\) yields the identity
(Namely, on the left-hand side of (6.22), the definition of the spectral family or the resolution of the identity \((E_\lambda )\) includes that the function \(\lambda \mapsto (E_\lambda g, e_n)_{\Pi }\) is of bounded variation on every closed interval and thus the Stieltjes integral is well defined. On the right-hand side, to commute the integral and inner product, one needs to check all the steps of the definition of the integral and also use Remark 6.4.)
We interpret (6.25) as the equality of two standard scalar Laplace-Stieltjes transforms. Hence, the measures \(( dE_\lambda g, e_n)_{\Pi }\) and \(\big ( \sum _{(k,q) \in {{\mathbb {K}}}_\lambda } F(x, \lambda , k,q), e_n \big )_{\Pi } \, d \lambda \) must be identical, and an integration with respect to \(\lambda \) over the set I yields
for all \(n\in {{\mathbb {N}}}\). This yields (6.22). \(\Box \)
7 Spectral projections of disjoint components
In this final section we consider the situation where the spectrum \(\sigma = \sigma _\textrm{ess} \subset {{\mathbb {R}}}^+\) is not connected and thus contains a bounded disjoint component S. This may happen for example in the case when the adjacent periodic cells are connected with thin ligaments, see e.g. [36]. We will represent the spectral projection onto the subspace \(H \subset L^2(\Pi )\), corresponding to the component S, by constructing an orthonormal basis of H consisting of Wannier functions and their translates. The results, of course, yield relevant information on the large time behavior of the solution of the heat equation (1.1)–(1.3). The emphasis here is in the simplicity of the presentation, and it should be mentioned that analogous results in other contexts are known to experts at least, see e.g. [25].
We assume that for some finite set \(K \subset {{\mathbb {N}}}\) of indices,
forms a disjoint, connected subset of the spectrum of \({{\mathcal {T}}}\), or, more generally, S of formula (7.1) is a finite union of such disjoint components of spectrum. We denote
is the spectral projection corresponding to the component S.
Proposition 7.1
The spectral projection P(S) can be written as
where \(G_k(\eta )\) and \(V_k\) are as in (5.9) and (3.10), respectively. Consequently, if \(g \in {{\mathcal {H}}}_S\), then the solution of the initial-boundary value problem (1.1)–(1.3) can be written as
Proof. From (6.22) we obtain
and the same change of integration variables as in the proof of (6.15) shows that
note that an analogue of the index set \(R_\lambda \) of formula (6.15) does not exist now, since we are assuming that \(\Upsilon _k\) is separated from the other spectral bands \(\Upsilon _k\) with \(k \notin K\). \(\Box \)
Next we single out an orthonormal basis for the space \({{\mathcal {H}}}_S\). In the context of self-adjoint elliptic operators with periodic coefficients in \({{\mathbb {R}}}^d\), like the Schrödinger operators with periodic potentials, the functions (7.7) are called Wannier functions. See [25], Definition 3.1. We give an elementary proof for our case.
Theorem 7.2
An orthonormal basis of the space \({{\mathcal {H}}}_S\) is formed by the functions
with \(k \in K\) and all translates \({{\textsf{F}}}^{-1} V_k \circ \tau _m\), \(m \in {{\mathbb {Z}}}\).
Here, \(f \circ \tau _m (x) = f (x_1 - m, x')\) for a function f of \(x \in \Pi \).
Proof. Note that \({{\textsf{F}}}^{-1} V_k(\cdot ;\eta ) \circ \tau _m = {{\textsf{F}}}^{-1} (e^{- im \eta } V_k(\cdot ;\eta ) )\). Denoting \(h:= {{\textsf{F}}}^{-1} V_k(\cdot ;\eta ) \circ \tau _m\) we thus have for every \(j \in K\), by (3.10),
so that \(P(S) h = h\) follows from the representation of the spectral projection in (7.3). Moreover, if \(k, j \in K\) and \(m, n \in {{\mathbb {Z}}}\) are arbitrary, we have
It remains to show that the functions span the subspace \({{\mathcal {H}}}_S\). Given \(h \in L^2(\Pi )\), we first note that both functions \({{\textsf{F}}}h (x;\eta ) \) and \(V_k (x;\eta )\) are \(2 \pi \)-periodic with respect to \(\eta \), hence the inner product \(G_k(\eta )\), defined as in (5.9) with h instead of g, is a continuous, periodic function on \([-\pi ,\pi ]\). We expand it as a Fourier series
and note that
The series converges in \(L^2(\Pi )\), due to (7.10) and the observation in (7.9) that the functions \({{\textsf{F}}}^{-1} V_k \circ \tau _{-n} \), \(k \in K\), \(n \in {{\mathbb {Z}}}\), form an orthonormal sequence. \(\Box \)
Let us finally consider the situation where the disjoint subset S of the spectrum in (7.1) consist of all bands \(\Upsilon _k\) with \(k=1,2, \ldots , N\) for some \(N > 1\), i.e., the set K in (7.1) equals the sequence \((1,2, \ldots , N)\). It is now plain from Theorem 7.2 that we have the decay rate
\(m_{N +1}\) as in (3.20), for all initial data g satisfying the orthogonality condition
Data availability statement
The manuscript does not have associated data.
References
S. Agmon, A. Douglis, L. Nirenberg, Estimates near the boundary for solutions of elliptic differential equations satisfying general boundary conditions. I. Comm. Pure Appl. Math. 12 (1959), 623–727.
F.L.Bakharev, S.A Nazarov, Gaps in the spectrum of a waveguide composed of domains with different limiting dimensions, Siberian Math. J. 56,4 (2015), 575–592.
F.L. Bakharev, K. Ruotsalainen, J. Taskinen, Spectral gaps for the linear surface wave model in periodic channels. Quaterly J. Mechanics Appl. Math. 67,3 (2014), 343–362.
F.L. Bakharev, J. Taskinen, Bands in the spectrum of a periodic elastic waveguide. Zeitschrift Angew. Math. Phys. 68, 5 (2017).
M.S. Birman, M.Z. Solomyak, Spectral theory of self–adjoint operators in Hilber space. Reidel Publ. Company, Dordrecht (1986).
H. Brezis, Functional analysis, Sobolev spaces and partial differential equations. Universitext, Springer, New York (2011).
G. Cardone, V. Minutolo, S.A. Nazarov, Gaps in the essential spectrum of periodic elastic waveguides, Zeitschrift für Angewandte Mathematik und Mechanik (ZAMM), in print.
G. Cardone, S.A. Nazarov, C. Perugia, A gap in the continuous spectrum of a cylindrical waveguide with a periodic perturbation of the surface, Math. Nachr., to appear.
V.I. Derguzov, On the discreteness of the spectrum of the periodic boundary-value problem associated with the study of periodic waveguides, Siberian Math. J. 21 (1980), 664–672.
F. Ferraresso, J. Taskinen, Singular perturbation Dirichlet problem in a double-periodic perforated plane. Ann. Univ. Ferrara 61,2 (2015), 277-290.
A. Figotin, P. Kuchment, Band-gap structure of spectra of periodic dielectric and acoustic media: I. Scalar model, SIAM J. Appl. Math. 56 (1996), 68-88; Band-gap structure of spectra of periodic dielectric and acoustic media: II. Two-dimensional photonic crystals, SIAM J. Appl. Math. 56 (1996), 1561-1620.
N. Filonov, Gaps in the spectrum of the Maxwell operator with periodic coefficients, Comm. Math. Physics. 240, 1-2 (2003), 161–170.
L. Friedlander, On the density of states of periodic media in the large coupling limit, Comm. Partial Diff. Equations 27 (2002), 355–380 .
L. Friedlander, Absolute continuity of the spectra of periodic waveguides, in Waves in periodic and random media (Peter Kuchment, editor) AMS Contemporary Mathematics 339 (2003), 37–42.
L. Friedlander, M. Solomyak, On the spectrum of narrow periodic waveguides, Russ. J. Math. Phys. 15,2 (2008), 238–242.
I.M. Gelfand, Expansion in characteristic functions of an equation with periodic coefficients (in Russian), Dokl.Akad.Nauk SSSR 73 (1950), 1117–1120.
E.L. Green, Spectral theory of Laplace-Beltrami operators with periodic metrics, J. Differential Equations 133 (1997), 15–29.
R. Hempel, K. Lineau, Spectral properties of the periodic media in large coupling limit, Comm. Partial Diff. Equations 25 (2000), 1445–1470.
J. Húska, Exponential separation and principal Floquet bundles for linear parabolic equations on general bounded domains: nondivergence case, Trans. Amer. Math. Soc. 360, 9 (2008), 4639–4679.
T. Hytönen, J. van Neerven, M. Veraar, L. Weis, Analysis in Banach Spaces, Ergebnisse der Mathematik und ihrer Grenzgebiete. 3. Folge / A Series of Modern Surveys in Mathematics (MATHE3, volume 63), Springer, Berlin (2016).
T. Kato, Perturbation theory for linear operators, 2nd edition, Springer-Verlag, Berlin–New York (1976).
W. Kohn, Analytic properties of Bloch waves and Wannier functions, Phys. Rev. 115 (1959), 809–21.
P. Kuchment, On the Floquet theory for parabolic and elliptic boundary value problems in a cylinder, Dokl. Akad. Nauk SSSR 258, 2 (1981), 296–299.
P. Kuchment, Floquet theory for partial differential equations, Birchäuser, Basel (1993).
P. Kuchment, Tight frames of exponentially decaying Wannier functions, J. Phys. A: Math. Theor. 42, 025203 (2009).
P. Kuchment, An overview of periodic elliptic operators, Bull. Amer. Math. Soc. 53,3 (2016), 343–414.
P. Kuchment, J. Taskinen, Dispersion relations and spectra of periodically perforated structures. Pure Appl. Functional Anal. 7, 2 (2022), 721–731.
O.A. Ladyzhenskaya, Boundary value problems of mathematical physics. Springer Verlag , New York (1985).
J.L. Lions, E. Magenes, Non–homogeneus boundary value problems and applications (in French), Dunod, Paris (1968). (English transl.: Springe–Verlag, Berlin–Heidelberg–New York (1972).)
A.I. Miloslavskii, On the Floquet theory for parabolic equations, Teor. Funktsii Funktsional. Anal. i Prilozhen 31(1979), 98–102 .
S.A. Nazarov, Properties of spectra of boundary value problems in cylindrical and quasicylindrical domains, Sobolev Spaces in Mathematics, vol. II (Maz’ya V., Ed.) International Mathematical Series 9 (2008), 261–309 .
S.A. Nazarov, A gap in the continuous spectrum of an elastic waveguide, C. R. Mecanique 336 (2008), 751–756.
S.A. Nazarov, Gap detection in the spectrum of an elastic periodic waveguide with a free surface, Zh. Vychisl. Mat. i Mat. Fiz. 49,2 (2009), 323–333. (English transl.: Comput. Math. and Math. Physics. 49,2 (2009), 332-343 .)
S.A Nazarov, Opening a gap in the essential spectrum of the elasticity problem in a periodic half-layer, Algebra i analiz 21,2 (2009), 166–202. (English transl.: St. Petersburg Math. J. 21,2 (2009).)
S.A Nazarov, B.A Plamenevskii, Elliptic problems in domains with piecewise smooth boundaries, Walter be Gruyter, Berlin, New York (1994).
S.A. Nazarov, K. Ruotsalainen, J. Taskinen , Essential spectrum of a periodic elastic waveguide may contain arbitrarily many gaps. Applicable Analysis 89 (2010), 109–124.
S.A. Nazarov, J. Taskinen , Spectral gaps for periodic piezoelectric waveguides. Z. Angew. Math. Phys. 66 (2015), 3017–3047.
G. Nenciu, Existence of the exponentially localized Wannier functions, Comm. Math. Phys. 91 (1983), 81–5.
G. Panati, Triviality of Bloch and Bloch-Dirac bundles, Ann. Henri Poincaré 8 (2007), 995-1011.
M. Reed, B.Simon, Methods of modern mathematical physics I, IV Academic Press San Diego (1980, 1978)
W. Rudin, Functional analysis, McGraw-Hill, New York (1973)
A.V. Sobolev, J. Walthoe, Absolute continuity in periodic waveguides, Proc. London Math. Soc. 85, 3 (2002), 717–741.
T.A. Suslina, On the absence of eigenvalues of a periodic matrix Schröodinger operator in a layer, Russ. J. Math. Phys. 8, 4 (2001), 463–486.
T.A. Suslina, R.G. Shterenberg, Absolute continuity of the spectrum of the magnetic Schrödinger operator with a metric in a two-dimensional periodic waveguide, St. Petersburg Math. J. 14,2 (2003), 305–343.
G. Teschl, Mathematical methods in quantum mechanics with applications to Schrödinger operators, second ed. Graduate studies in mathematics 157, American Mathematical Society Providence, Rhode Island.
H. Weyl, Das asymptotische Verteilungsgesetz der Eigenwerte linearer partieller Differentialgleichungen, Math. Ann. 71 (1911), 441–479.
K. Yoshitomi, Band gap of the spectrum in periodically curved quantum waveguides, J.Diff. Equations 142(1998), 123–166 .
V. Zhikov, Gaps in the spectrum of some elliptic operators in divergent form with periodic coefficients, Algebra i Analiz 16,5 (2004), 34–58. (English transl.: St. Petersburg Math. J. 16,5 (2005), 773–790.)
Funding
Open Access funding provided by University of Helsinki (including Helsinki University Central Hospital).
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors declare that they have no conflict of interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The first named author was supported by the Magnus Ehrnrooth Foundation. The second named author was partially supported by the Väisälä Foundation of the Finnish Academy of Science and Letters.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Rosenberg, M., Taskinen, J. Some aspects of the Floquet theory for the heat equation in a periodic domain. J. Evol. Equ. 24, 23 (2024). https://doi.org/10.1007/s00028-024-00951-0
Accepted:
Published:
DOI: https://doi.org/10.1007/s00028-024-00951-0