Quasi-periodic two-scale homogenisation and effective spatial dispersion in high-contrast media

TL;DR

This paper extends two-scale homogenisation theory to high-contrast media with disconnected stiff components, revealing complex spectral behavior and spatial dispersion effects due to quasimomentum dependence.

Contribution

It introduces an extended two-scale convergence framework to accurately capture spectral properties in high-contrast media with disconnected stiff phases.

Findings

Persistent waves of all periods are shown to exist asymptotically.

Homogenised equations exhibit non-trivial quasimomentum dependence.

Spectral limits reveal rich dispersion phenomena.

Abstract

The convergence of spectra via two-scale convergence for double-porosity models is well known. A crucial assumption in these works is that the stiff component of the body forms a connected set. We show that under a relaxation of this assumption the (periodic) two-scale limit of the operator is insufficient to capture the full asymptotic spectral properties of high-contrast periodic media. Asymptotically, waves of all periods (or quasi-momenta) are shown to persist and an appropriate extension of the notion of two-scale convergence is introduced. As a result, homogenised limit equations with none trivial quasimomentum dependence are found as resolvent limits of the original operator family, resulting in limiting spectral behaviour with a rich dependence on quasimomenta.

Figures

Equations

$$-{\rm div}\big{(}a_{\varepsilon}(\tfrac{x}{\varepsilon})\boldsymbol{\nabla}u\big{)}=\omega^{2}u\qquad\text{in $\Omega$}$$

$$a_{\epsilon }(y)=\{\begin{array}{cc}{a}_{1\epsilon }(y),& y\in Q_1,\\ a_{0\epsilon }(y),& y\in Q_0,\end{array}$$

$$a_{1\epsilon }=a_1,a_{0\epsilon }=a_0,\text{for bounded elliptic a0,a1.}$$

$$a_{1\epsilon }=a_1,a_{0\epsilon }={\epsilon }^2{a}_0.$$

$$\mathbf{u}=(u_1,u_2,u_3)=\sum _{i=1}^3{u}_i{\mathbf{e}}_i.$$

$$C_{\mathbf{\theta }}^{\infty }(\square )=\{\varphi |\varphi =e^{\mathrm{i}\mathbf{\theta }\cdot \mathrm{y}}\psi ,\psi \in C_{\#}^{\infty }(\square )\}.$$

$$H_{\mathbf{\theta }}^1(\square ):=\{e^{\mathrm{i}\mathbf{\theta }\cdot \mathrm{y}}{u}_{\#}|u_{\#}\in H_{\#}^1(\square )\}.$$

$$A_{\varepsilon}u:=-{\rm div}\big{(}a_{\varepsilon}(\tfrac{\cdot}{\varepsilon})\boldsymbol{\nabla}u\big{)}$$

$${\nu }_{1}I\le a_{\epsilon }\le {\nu }_{2}I.$$

$$\sigma (A_{\epsilon })=\bigcup _{\mathbf{\Theta }\in {\left[0,\frac{2\pi }{\epsilon }\right)}^d}\sigma (A_{\epsilon }(\mathbf{\Theta }))$$

$$\displaystyle-{\rm div}\big{(}a_{\varepsilon}(\tfrac{x}{\varepsilon})\boldsymbol{\nabla}e^{\rm i{\boldsymbol{\Theta}}\cdot x}u\big{)}$$

$$\displaystyle-{\rm div}\big{(}a_{\varepsilon}(\tfrac{x}{\varepsilon})\boldsymbol{\nabla}e^{\rm i{\boldsymbol{\Theta}}\cdot x}u_{\varepsilon}\big{)}$$

$$\displaystyle-{\rm div}\big{(}\varepsilon^{-2}a_{\varepsilon}(y)\boldsymbol{\nabla}w_{\varepsilon}\big{)}$$

$$w_{\epsilon }(y+z)=e^{\mathrm{i}\mathbf{\theta }\cdot \mathrm{z}}{w}_{\epsilon }(y),y\in \square ,z\in {\mathbb{Z}}^d,$$

$$H_{\mathbf{\theta }}^1(\square ):=\{e^{\mathrm{i}\mathbf{\theta }\cdot \mathrm{y}}{u}_{\#}|u_{\#}\in H_{\#}^1(\square )\}.$$

$$C_{\mathbf{\theta }}^{\infty }(\square )=\{\varphi |\varphi =e^{\mathrm{i}\mathbf{\theta }\cdot \mathrm{y}}\psi ,\psi \in C_{\#}^{\infty }(\square )\}.$$

$$\underset{\epsilon \to 0}{\lim }{\int }_{\Omega }{\left|\varphi (x)\phi \left(\frac{x}{\epsilon }\right)\right|}^2\mathrm{d}x={\int }_{\Omega }{\int }_{\square }{\left|\varphi (x)\phi (y)\right|}^2,$$

$$\underset{\epsilon \to 0}{\lim }{\int }_{\Omega }{u}_{\epsilon }(x)\overline{\varphi (x)\phi \left(\frac{x}{\epsilon }\right)}\mathrm{d}x={\int }_{\Omega }{\int }_{\square }u(x,y)\overline{\varphi (x)\phi (y)}\mathrm{d}y\mathrm{d}x,\forall \varphi \in C_0^{\infty }(\Omega ),\forall \phi \in C_{\mathbf{\theta }}^{\infty }(\square )$$

$${\int }_{\Omega }{u}_{\epsilon }(x)\overline{\varphi (x)\phi (\frac{x}{\epsilon })}\mathrm{d}x$$

$$\displaystyle\int_{\Omega}\int_{\square}{\rm exp}\big{(}\rm i{\boldsymbol{\theta}}\cdot y\big{)}\widetilde{u}(x,y)\overline{\phi(x)\varphi(y)}\,{\rm d}y{\rm d}x$$

$$\underset{\epsilon \to 0}{\lim }{\int }_{\Omega }{u}_{\epsilon }(x)\overline{\varphi (x)\psi \left(\frac{x}{\epsilon }\right)}\mathrm{d}x={\int }_{\Omega }{\int }_{\square }u(x,y)$$

$$\exp (-\mathrm{i}\mathbf{\theta }\cdot \frac{\mathrm{x}}{\epsilon }){\mathrm{u}}_{\epsilon }(\mathrm{x})\rightharpoonup {\int }_{\square }\exp (-\mathrm{i}\mathbf{\theta }\cdot \mathrm{y})\mathrm{u}(\mathrm{x},\mathrm{y})\mathrm{d}\mathrm{y}\text{weakly in L2(Ω).}$$

$$\underset{\epsilon \to 0}{\mathrm{liminf}}{\int }_{\Omega }|u_{\epsilon }(x){|}^2\mathrm{d}x\ge {\int }_{\Omega }{\int }_{\square }|u(x,y){|}^2\mathrm{d}y\mathrm{d}x.$$

$$\underset{\epsilon }{\sup }||u_{\epsilon }|{|}_{L^2(\Omega )}<\infty ,$$

$$u_{\epsilon }\stackrel{2-\mathbf{\theta }}{\rightharpoonup }u,$$

$$\epsilon u_{\epsilon }\stackrel{2-\mathbf{\theta }}{\rightharpoonup }0.$$

$${\int }_{\Omega }\epsilon \nabla u_{\epsilon }(x)\cdot \overline{\varphi (x)\mathbf{\Psi }(\frac{x}{\epsilon })}\mathrm{d}x$$

$$=-{\int }_{\Omega }{u}_{\epsilon }(x)\overline{\epsilon {\nabla }_x\varphi (x)\cdot \mathbf{\Psi }(\frac{x}{\epsilon })}\mathrm{d}x-{\int }_{\Omega }{u}_{\epsilon }(x)\overline{\varphi (x){\mathrm{div}}_y\mathbf{\Psi }(\frac{x}{\epsilon })}\mathrm{d}x,$$

$${\int }_{\Omega }{\int }_{\square }\mathbf{\chi }(x,y)\cdot \overline{\varphi (x)\mathbf{\Psi }(y)}\mathrm{d}y\mathrm{d}x=-{\int }_{\Omega }{\int }_{\square }u(x,y)\overline{\varphi (x){\mathrm{div}}_y\mathbf{\Psi }(y)}\mathrm{d}y\mathrm{d}x.$$

$${\int }_{\square }\mathbf{\chi }(x,y)\cdot \overline{\mathbf{\Psi }(y)}\mathrm{d}y=-{\int }_{\square }u(x,y)\overline{{\mathrm{div}}_y\mathbf{\Psi }(y)}\mathrm{d}y,\forall \mathbf{\Psi }\in C_{\mathbf{\theta }}^{\infty }(\square ;{\mathbb{C}}^d).$$

Full text

1 Introduction

The model problem to study time-harmonic waves, with frequency $\omega$, in media with microstructure is

[TABLE]

where the wave $u$ represents the information being propagated, such as pressure in acoustics, deformation in elasticity or electromagnetic fields in electromagnetism111In elasticity and electromagnetism the wave equation describes certain polarised waves: e.g. Shear polarised wave in elasticity or Transverse Electric and Transverse Magnetic polarised waves for the Maxwell system.. The microstructured nature of the media is characterised by periodic coefficients a_{\varepsilon}\222The implied non-trivial dependence of $a$ on $\varepsilon$ is deliberate and, as we shall see, important. :

[TABLE]

where $a_{0\varepsilon}$, $a_{1\varepsilon}$ are (the square-root of) the wave speeds of the individual constitutive material components, see figure 1.

The parameter $\varepsilon$ represents the ratio between the size of the microstructure and the observable length scale, and is typically taken to be small. From the point of view of applications, it is important to study the asymptotic behaviour of these waves in the limit of vanishing $\varepsilon$.

A classical approximation, provided by the homogenisation theorem333Also called the long-wavelength or quasi-static approximation depending on the community., states that for fixed frequency $\omega$ the microstructured media can be approximated by an ‘effective’ homogeneous media whose wave speed $a^{\rm hom}$ is constant and determined directly from the ‘local periodic’ behaviour of the problem. The intuition behind why the homogenisation theorem holds is that the ‘wavelength’ of $u$ is long with respect to the microstructure: variations in $u$ appear over much longer distances than the media’s period. Mathematically, this is ensured by assuming that $a_{\varepsilon}$ are taken to be uniformly bounded and elliptic with respect to $\varepsilon$, for example

[TABLE]

It has been known for some time now that interesting effects appear when the above elliptic conditions are not uniform. This happens for example in so-called high-contrast media. In the context of waves, high-contrast media of particular interest are the so-called double porosity models which admit the ‘critical’ scaling:

[TABLE]

Physically, this critical scaling corresponds to the wavelength of $u$ remaining ‘long’ within the media $Q_{1}$ but in media $Q_{0}$ the wavelength is at the ‘resonant’ scale, i.e. of the same order as the size of the microstructure. Thus violating the underlying intuition for the long-wavelength approximation.

The mathematical analysis of high-contrast problems has given rise to rigorous descriptions of various scale-interaction phenomena such as memory effects and other non-local effects (e.g. [17, 28, 5, 7, 8, 9, 14]). Within the context of wave propagation, an important feature of high-contrast problems is that they contain spectral gaps (cf. [18, 30, 31]): frequencies at which no wave can propagate through the underlying medium. Such gaps are important from the point of view of wave-guiding applications such as photonic crystal fibres.

An important initial work in the study of the spectrum of high-contrast elliptic operators was undertaken by V. V. Zhikov [30, 31]. Therein, the homogenisation theory for double porosity-type problems was developed within the framework of the so-called two-scale convergence of G. Nguetseng-G. Allaire [27, 1]. Using this theory, Zhikov derived two-scale limit spectral equations that contain a non-trivial coupling between micro- and macro-scales. Such a coupling leads to an eigenvalue problem with a highly non-linear spectral dependence, described by a function $\beta$. The convergence of spectra (in the appropriate sense) was proved and, by doing so, demonstrates that this $\beta$ function provides an explicit description of the asymptotic structure of the spectrum. Such an explicit description of the limit spectral behaviour via two-scale homogenisation has made way for mathematical studies of high-contrast media as wave-guides: in [20] using multi-scale asymptotics and supplemented with analysis based on two-scale convergence in [10].

Moreover, the Zhikov $\beta$ function was later independently discovered by G. Bouchitté and D. Felbacq [6] in the specific context of TM-polarised electromagnetic wave propagation in a dilute dielectric two-dimensional photonic crystal fibre; therein the authors made the interesting interpretation of the $\beta$-function playing the role of effective negative magnetism. Later, in the context of elasticity, a matrix analogue of the $\beta$ function is derived and plays the role of frequency-dependent effective density [3, 4, 32]. Such works demonstrate that the unusual phenomena observed in high-contrast media can be described by non-standard constitutive laws provided via two-scale homogenisation.

The idea that high-contrast media can result in the appearance of non-standard constitutive laws and give rise to composite media with complex wave phenomena near micro-resonances has prompted a recent energetic pursuit of such laws in the contexts of elasticity (e.g. [24]) and electromagnetism (e.g. [22, 12]). Applications can be found in areas such as cloaking (e.g. [25, 26]). It was shown in the work of V. P. Smyshlyaev [29], building on related ideas in [14], that the two-scale homogenised limit of various anisotropic elastic media contain not only the temporal non-locality (as described by the Zhikov $\beta$ function) but also exhibit spatial non-locality. The presence of which leads to the phenomena of ‘directional’ localisation: the number of admissible propagating wave modes depends not only on the frequency but on the direction of propagation. Such a feature is important for cloaking applications. These motivating works have led to recent systematic study containing rigorous asymptotic and spectral analysis of general mathematical constructions containing ‘high-contrasts’ [21]. Analysis based on the work [21] has led to the demonstration that the two-scale convergence is insufficient to fully study the spectrum of general high-contrast problems, see [15, Chapter 5],[11]. The reason for this inconsistency is due to the presence of quasi-periodic micro-oscillations that persist at leading-order in general high-contrast media.

In this work we appropriately develop homogenisation theory to study quasi-periodic micro-oscillations. This is achieved by extending the two-scale convergence framework to admit such oscillations. We explain the (lack of the) role these micro-oscillations in the numerous previous works on high-contrast problems. Then, we apply this theory in the spectral analysis of a novel class of high-contrast media. In particular, we shall show that by relaxing the geometric assumptions in the double-porosity model leads to multi-scale homogenised models that contain a new feature: the effective wave speed depends on the quasi-momenta in a highly discontinuous fashion. Specifically, the non-standard constitutive equations for such high-contrast media exhibit spatial dispersion. The presence of this novel feature is related to the contribution of the quasi-periodic waves on the microscale.

Notations

We end the introduction with some words on the notation used in this article.

Vectors and vector-valued functions are represented by lower-case boldface symbols with the exception of the co-ordinate points. $\{{\bf{e}}_{1},{\bf{e}}_{2},{\bf{e}}_{3}\}$ denotes the Euclidean basis in $\mathbb{R}^{3}$. For a vector ${\bf{u}}\in\mathbb{R}^{3}$, we denote by $u_{i}$ its component with respect to ${\bf{e}}_{i}$, and write

[TABLE]

Points in $\mathbb{R}^{3}$ will be denoted by the symbol $x$ and points in the unit cell $\square:=[0,1)^{3}$ will be denoted by $y$. The notation $\partial_{i}$ will be used to denote partial differentiation with respect to the $i$-th coordinate variable, and we shall replace the suffix $i$ with $x_{i}$ or $y_{i}$ when we wish to emphasis the macroscopic or microscopic variable. Similarly, the notion ${\rm div}_{x}$, ${\rm div}_{y}$, $\boldsymbol{\nabla}_{x}$, and $\boldsymbol{\nabla}_{y}$, are used for the divergence or gradient of a function in terms of $x$ or $y$.

Throughout $\Omega$ is a domain in $\mathbb{R}^{d}$, $d\geq 1$, $\square:=(0,1)^{d}$ and ${{\boldsymbol{\theta}}}\in[0,2\pi)^{d}$. All of the functions, even if real-valued, are considered to take values in the complex field.

The space $C^{\infty}_{\#}(\square)$ denotes the usual space of smooth $\square$-periodic functions. Whereas, $C^{\infty}_{{\boldsymbol{\theta}}}(\square)$ shall denote the space of smooth functions $\varphi(y)$ whose functions and derivatives are ${{\boldsymbol{\theta}}}$-quasi-periodic with respect to $y$: $\varphi(y+{\bf{e}}_{j})=\mathrm{exp}(\rm i{{\boldsymbol{\theta}}}_{j})\varphi(y)$ for each $y\in\square$ and each Euclidean basis vector ${\bf{e}}_{j}$, $j=1,\ldots,d$. Equivalently,

[TABLE]

Note that $C^{\infty}_{\boldsymbol{0}}(\square)=C^{\infty}_{\#}(\square)$ and use the latter to avoid confusion with the notation for the space compactly supported smooth functions.

The Sobolev space $H^{1}_{\#}(\square)$ is the usual Sobolev space of $H^{1}$ $\square$-periodic functions. Whereas $H^{1}_{\boldsymbol{\theta}}(\square)$ is defined as the closure of $C^{\infty}_{\boldsymbol{\theta}}(\square)$ with respect to the $H^{1}$ norm, or equivalently as

[TABLE]

Also, we note $H^{1}_{\boldsymbol{0}}(\square)=H^{1}_{\#}(\square)$ and, in this situation, we use the latter to avoid confusion with the Sobolev space of zero trace $H^{1}$ functions.

For subsets $\{S_{\varepsilon}\}_{\varepsilon}$ and $S$ of $\mathbb{R}^{d}$ we say that $S_{\varepsilon}$ converges to $S$ in the Hausdorff sense if the following conditions hold:

For every $\lambda_{\varepsilon}\in S_{\varepsilon}$ such that $\lambda_{\varepsilon}$ converges to some $\lambda_{0}$, then $\lambda_{0}\in S$. 2. 2.

For every $\lambda_{0}\in S$ there exists$\lambda_{\varepsilon}\in S_{\varepsilon}$ such that $\lim_{\varepsilon}\lambda_{\varepsilon}=\lambda_{0}$.

We shall use the notation $\lim_{\varepsilon}S_{\varepsilon}=S$ when a sequence of sets $S_{\varepsilon}$ Hausdorff converges to $S$.

The Einstein summation convention will not be used in this article, that is we do not sum with respect to repeated indices.

2 Quasi-periodic two-scale convergence

In this section we introduce an appropriate notion of convergence that will account for the presence of microscopic oscillations that are quasi-periodic in nature. This convergence will turn out to be a natural extension of the standard (periodic) two-scale convergence introduced by G. Nguetseng[27]-G. Allaire [1]. In particular, we aim to use this extended notion of two-scale convergence to study the spectral properties of operator families in homogenisation theory in a similar vein to that first introduced by V. V. Zhikov [30, 31].

2.1 Motivation

We shall motivate the notion of quasi-periodic two-scale convergence here. This motivation is based on the principle goal of characterising the spectral asymptotics of high-contrast elliptic operators.

Let $\varepsilon$ be a sequence of positive real numbers with limit zero. Consider the differential operator $A_{\varepsilon}:\mathcal{D}(A_{\varepsilon})\subset L^{2}(\mathbb{R}^{d})\rightarrow L^{2}(\mathbb{R}^{d})$ whose action is given by

[TABLE]

and domain $\mathcal{D}(A_{\varepsilon})$ consists of $u$ for which $A_{\varepsilon}u\in L^{2}(\mathbb{R}^{d})$. Here $a_{\varepsilon}$ are $\square$-periodic measurable functions, that are bounded and elliptic: $\exists\nu_{1},\nu_{2}>0$ such that

[TABLE]

In this article, we focus on $a_{\varepsilon}$ that are uniformly bounded, i.e. $\nu_{2}$ is independent of $\varepsilon$, but $a_{\varepsilon}$ may degenerate in the sense that $\nu_{1}=\nu_{1}(\varepsilon)$ with $\lim_{\varepsilon}\nu_{1}(\varepsilon)\geq 0$. We are interested in analysing the structure of the spectrum $\sigma(A_{\varepsilon})$ of $A_{\varepsilon}$ in the limit of $\varepsilon$. The strategy of the study is to establish the existence of some operator $A_{0}$ such that $\sigma(A_{\varepsilon})$ Hausdorff converges to $\sigma(A_{0})$; i.e. the following conditions hold:

For every $\lambda_{\varepsilon}\in\sigma(A_{\varepsilon})$ such that $\lambda_{\varepsilon}$ converges to some $\lambda_{0}$ we deduce that $\lambda_{0}\in\sigma(A_{0})$. 2. 2.

For every $\lambda_{0}\in\sigma(A_{0})$ we find $\lambda_{\varepsilon}\in\sigma(A_{\varepsilon})$ such that $\lim_{\varepsilon}\lambda_{\varepsilon}=\lambda_{0}$.

A crucial question is how to determine the operator $A_{0}$. For example, in classical and semi-classical high-contrast problems, $A_{0}$ turns out the be the strong two-scale resolvent homogenised limit of $A_{\varepsilon}$, cf. [30, 31, 16]. To develop some intuition on what to expect in the general case, let us recall an important result from the spectral theory of elliptic operators with $\varepsilon\square$-periodic coefficients: the Floquet-Bloch decomposition (see for example [23] for more details). This result states that the following characterisation of $\sigma(A_{\varepsilon})$ holds:

[TABLE]

where $A_{\varepsilon}({\boldsymbol{\Theta}}):\mathcal{D}(A_{\varepsilon}({\boldsymbol{\Theta}}))\subset L^{2}(\varepsilon\square)\rightarrow L^{2}(\varepsilon\square)$, describe a family of densely defined self-adjoint operators with compact resolvent given by the action that $A_{\varepsilon}({\boldsymbol{\Theta}})u=f\in L^{2}(\square)$ if $u\in H^{1}_{\#}(\varepsilon\square)$ solves

[TABLE]

Taking the above into consideration we see that $\lambda_{\varepsilon}\in\sigma(A_{\varepsilon})$ if, and only if, there exists ${\boldsymbol{\Theta}}_{\varepsilon}\in\left[0,\tfrac{2\pi}{\varepsilon}\right)^{d}$ and non-trivial $u_{\varepsilon}\in H^{1}_{\#}(\varepsilon\square)$ such that

[TABLE]

By a change of variables $y=x/\varepsilon$ and ${\boldsymbol{\theta}}=\varepsilon{\boldsymbol{\Theta}}$, we see that $w_{\varepsilon}(y):=e^{\rm i{{\boldsymbol{\theta}}}\cdot y}u_{\varepsilon}(\varepsilon y)$ solves

[TABLE]

and $w_{\varepsilon}$ belongs to the space of $H^{1}(\square)$ functions that satisfy the condition

[TABLE]

for some ${\boldsymbol{\theta}}\in[0,2\pi)^{d}$. This condition is typically referred to as the Bloch or quasi-periodic condition and ${\boldsymbol{\theta}}$ is known as the quasi-momentum. Note that ${\boldsymbol{\theta}}=\bf 0$ is the usual periodicity condition. The Sobolev space of $H^{1}(\square)$ ${\boldsymbol{\theta}}$-quasi-periodic functions coincides with $H^{1}_{\boldsymbol{\theta}}(\square)$ which, we recall, is be defined as

[TABLE]

The general principle to observe here is that if we wish to study the asymptotic behaviour of the spectrum $\sigma(A_{\varepsilon})$ we need to keep track of the eigenfunctions that are ${\boldsymbol{\theta}}$-quasi-periodic on micro-scale $y:=x/\varepsilon$, for all ${\boldsymbol{\theta}}\in[0,2\pi)^{d}$. The notion of quasi-periodic two-scale convergence, introduced in Section 2.2 below, performs such a task.

We note here that in the case of the whole space, as discussed above, one need not refer to a notion of quasi-periodic two-scale convergence to study the asymptotics of the spectrum; one may study the norm-resolvent limits of the operators $A_{\varepsilon}(\theta)$ to study spectral asymptotics, cf. [18] where the point-wise (in $\theta$) limits or [13] where the uniform limits were considered in the double-porosity setting. That being said, for boundary-value problems, the Bloch decomposition does not hold; nevertheless, the whole space (or Bloch) spectrum is expected to contribute asymptotically to the bounded domain spectrum and the precluding discussion is still relevant. It is this setting that the quasi-periodic two-scale convergence will be particularly useful.

Finally, we comment that the above discussion leads to the natural question: why in previous cases considered was it sufficient to consider the standard (periodic) two-scale limit of $A_{\varepsilon}$ to ensure spectral convergence? Or put another way, when in the asymptotic limit of $\varepsilon$ do we need consider all quasi-periodicity and not just ${\boldsymbol{\theta}}=\bf 0$. This shall be explained in Section 2.3.

2.2 Definition and basic properties

This section is dedicated to the introduction of the notion quasi-periodic two-scale convergence and an exposition of results that are appropriate to homogenisation theory.

Recall $C^{\infty}_{\#}(\square)$ denotes the usual space of smooth $\square$-periodic functions and

[TABLE]

The following mean-value property will be important: For every $\varphi\in C^{\infty}_{{\boldsymbol{\theta}}}(\square)$ and every $\phi\in C^{\infty}_{0}(\Omega)$ the following convergence

[TABLE]

holds. This fact follows by noting that the assertion holds for elements in $C^{\infty}_{\#}(\square)$, see for example [1, Lemma 1.3], and observing that multiplication by $\mathrm{exp}(-\rm i{{\boldsymbol{\theta}}}\cdot y)$ defines an isomorphism between $C^{\infty}_{{\boldsymbol{\theta}}}(\square)$ and $C^{\infty}_{\#}(\square)$ that preserves absolute value. Indeed, $\varphi$ belongs to $C^{\infty}_{{\boldsymbol{\theta}}}(\square)$ if, and only if, $\mathrm{exp}(-\rm i{{\boldsymbol{\theta}}}\cdot y)\varphi$ belongs to $C^{\infty}_{\#}(\square)$ and $|\varphi|=|\mathrm{exp}(-\rm i{{\boldsymbol{\theta}}}\cdot y)\varphi|$ in $\square$.

We remark here that because $C^{\infty}_{\boldsymbol{\theta}}(\square)$ is isomorphic to $C^{\infty}_{\#}(\square)=C^{\infty}_{\boldsymbol{0}}(\square)$, with isomorphism {\rm exp}\big{(}-\rm i{\boldsymbol{\theta}}\cdot y\big{)}, then the results presented in this section444The results in this section can be established by first principles making no reference to such an isomorphism. are immediately established for each ${\boldsymbol{\theta}}\in[0,2\pi)^{3}$ if proved for ${\boldsymbol{\theta}}=\boldsymbol{0}$. We shall demonstrate this with the first result of the section and omit the remaining proofs which follow in a similar manner.

Definition 1**.**

Let $u_{\varepsilon}\in L^{2}(\Omega)$ be a bounded sequence and $u\in L^{2}(\Omega\times\square)$. Then, we say $u_{\varepsilon}$ (weakly) ${{\boldsymbol{\theta}}}$-quasi-periodic two-scale converges to $u$, denoted by $u_{\varepsilon}\overset{2-{{\boldsymbol{\theta}}}}{\rightharpoonup}u$, if the following convergence

[TABLE]

holds.

Remark 1*.*

Notice that for ${{\boldsymbol{\theta}}}=0$, this is the standard notion of two-scale convergence.

The next important result states that bounded sequences in $L^{2}(\Omega)$ are relatively compact with respect to quasi-periodic two-scale convergence.

Proposition 1**.**

If $u_{\varepsilon}$ is bounded in $L^{2}(\Omega)$ then, up to a subsequence, $u_{\varepsilon}$ weakly ${{\boldsymbol{\theta}}}$-quasi-periodic two-scale converges to some $u\in L^{2}(\Omega\times\square)$.

Proof.

The result has been established previously for the case ${\boldsymbol{\theta}}=\boldsymbol{0}$, see for example [27, 1, 30]. Let us consider ${\boldsymbol{\theta}}\neq\boldsymbol{0}$. Note that the function \widetilde{u}_{\varepsilon}={\rm exp}\big{(}-\rm i{\boldsymbol{\theta}}\cdot y\big{)}u_{\varepsilon} is bounded in $L^{2}$ and therefore by the assertion for ${\boldsymbol{\theta}}=\boldsymbol{0}$, up to a discarded subsequence, $\widetilde{u}_{\varepsilon}$ ($\boldsymbol{0}$-quasi-periodically) two-scale converges to some $\widetilde{u}\in L^{2}(\Omega\times Q)$. Now, the result follows from this fact and noting that for fixed $\varphi\in C^{\infty}_{\boldsymbol{\theta}}(\square)$ one has

[TABLE]

and that {\rm exp}\big{(}-\rm i{\boldsymbol{\theta}}\cdot y\big{)} is a smooth periodic function. Hence u_{\varepsilon}\overset{2-{{\boldsymbol{\theta}}}}{\rightharpoonup}{\rm exp}\big{(}\rm i{\boldsymbol{\theta}}\cdot y\big{)}\widetilde{u}. ∎

An important result from the point of view of homogenisation theory is that the test functions $\varphi$ in (2.3) can be taken to be quasi-periodic elements of $L^{2}(\square)$, i.e. the following result holds.

Proposition 2**.**

If $u_{\varepsilon}\in L^{2}(\Omega)$ ${{\boldsymbol{\theta}}}$-quasi-periodic two-scale converges to $u\in L^{2}(\Omega\times\square)$, then the following convergence

[TABLE]

holds for all $\phi\in C^{\infty}_{0}(\Omega)$, and for all $\psi\in L^{2}(\square)$ such that $\psi(y+{\bf{e}}_{j})=\mathrm{exp}(\rm i{{\boldsymbol{\theta}}}_{j})\psi(y)$ for almost every $y\in\square$ and $j=1,\ldots,d$.

Remark 2*.*

If $\Omega$ is a bounded domain, as in this article, then additionally the test functions $\phi$ can be taken to be elements of $C(\overline{\Omega})$.

The following results are of interest.

Proposition 3**.**

For $u_{\varepsilon}\in L^{2}(\Omega)$ ${{\boldsymbol{\theta}}}$-quasi-periodic two-scale converging to $u\in L^{2}(\Omega\times\square)$ one has that

[TABLE] 2. 2.

For $u_{\varepsilon}\in L^{2}(\Omega)$ ${{\boldsymbol{\theta}}}$-quasi-periodic two-scale converging to $u\in L^{2}(\Omega\times\square)$ then

[TABLE]

A result of particular interest in high-contrast homogenisation problems is the following.

Proposition 4**.**

Let $u_{\varepsilon}\in H^{1}(\Omega)$ satisfy

[TABLE]

Then, there exists $u\in L^{2}(\Omega;H^{1}_{\boldsymbol{\theta}}(\square))$ such that, up to a subsequence, the following convergences hold:

[TABLE]

Recall here that $H^{1}_{\boldsymbol{\theta}}(\square)$ is given by (1.1).

Proof.

Let $\phi$ and $\boldsymbol{\Psi}$ denote respectively fixed arbitrary elements of $C^{\infty}_{0}(\Omega)$ and $C^{\infty}_{{\boldsymbol{\theta}}}(\square;\mathbb{C}^{d})$. By Proposition 1, there exists $u\in L^{2}(\Omega\times\square)$ and $\boldsymbol{\chi}\in L^{2}(\Omega\times\square;\mathbb{C}^{d})$ such that, up to a discarded subsequence, the following convergences hold:

[TABLE]

Note that, since $u_{\varepsilon}$ is bounded in $L^{2}(\Omega)$, then $\varepsilon u_{\varepsilon}$ strongly converges to zero in $L^{2}(\Omega)$ and from Proposition 3 part 2. we conclude that

[TABLE]

Let us prove that $u\in L^{2}(\Omega;H^{1}_{{\boldsymbol{\theta}}}(Q))$. Using the convergences (2.4) and (2.5) we pass to the limit in the identity

[TABLE]

to deduce that

[TABLE]

Therefore, for almost every $x$, the functions $\boldsymbol{\chi}(x,\cdot)$ and $u(x,\cdot)$ are related by the identity

[TABLE]

It is clear that C^{\infty}_{0}\big{(}(0,1)^{d}\big{)}\subset C^{\infty}_{{\boldsymbol{\theta}}}(\square) and so $u\in H^{1}(\square)$ with $\boldsymbol{\nabla}_{y}u=\boldsymbol{\chi}$. It remains to show $u$ belongs to $H^{1}_{\boldsymbol{\theta}}(\square)$. This follows from noting that after performing integration by parts in the above identity we arrive at

[TABLE]

Setting $\boldsymbol{\Psi}=e^{\rm i{\boldsymbol{\theta}}\cdot y}\boldsymbol{\Psi}_{\#}$ above, for arbitrary smooth $\square$-periodic $\boldsymbol{\Psi}_{\#}$, demonstrates that $u_{\#}(x,\cdot):=e^{-\rm i{\boldsymbol{\theta}}\cdot y}u(x,\cdot)$ is an element of $H^{1}(\square)$ that satisfies periodic boundary conditions with respect to $y$. That is, $u_{\#}(x,\cdot)\in H^{1}_{\#}(\square)$ and so (see definition (1.1)) $u(x,\cdot)\in H^{1}_{\boldsymbol{\theta}}(\square)$. ∎

We end this section with a result that is illuminating when it comes to studying the convergence of spectra for parameter-dependent operator families. It readily provides a one-sided justification for the Hausdorff convergence of the high-contrast spectra to the spectrum associated to quasi-periodic two-scale limits. This result is based on the following definition, which extends the notion of strong resolvent two-scale convergence first introduced by V. V. Zhikov in [30, 31].

Definition 2**.**

Fix ${{\boldsymbol{\theta}}}\in[0,2\pi)^{d}$, and let $A_{\varepsilon}$ and $A$ be non-negative self-adjoint operators respectively defined in $L^{2}(\Omega)$ and $H$ a closed subset of $L^{2}(\Omega\times\square)$. We say that $A_{\varepsilon}$ strong resolvent ${{\boldsymbol{\theta}}}$-quasi-periodic two-scale converges to $A$ if for every $f_{\varepsilon}(x)\in L^{2}(\Omega)$ that ${{\boldsymbol{\theta}}}$-quasi-periodic two-scale converges to $f(x,y)\in L^{2}(\Omega\times\square)$, the following convergence

[TABLE]

holds. Here, $P$ is the orthogonal projection onto $H$ in $L^{2}(\Omega\times\square)$.

Here, we state an important consequence of such resolvent convergence. The proof, omitted here, follows standard spectral theoretic arguments, see for example [30].

Proposition 5**.**

If $A_{\varepsilon}$ strong resolvent ${{\boldsymbol{\theta}}}$-quasi-periodic two-scale converges to $A$ then the spectrum $\sigma(A)$ of $A$ is related to the spectrum $\sigma(A_{\varepsilon})$ of $A_{\varepsilon}$ in the following sense:

For every $\lambda\in\sigma(A)$ there exists $\lambda_{\varepsilon}\in\sigma(A_{\varepsilon})$ such that $\lambda_{\varepsilon}$ converges to $\lambda$ as $\varepsilon$ tends to zero.

2.3 On the relevance of quasi-periodic two-scale convergence in spectral asymptotics

Proposition 5 informs us that, in principle, one should consider all strong quasi-periodic two-scale limits of an operator $A_{\varepsilon}$ to fully characterise its limit spectrum (in the Hausdorff sense). Yet, clearly this is not always the case: such a notion of convergence has not appeared previously, nor was it needed, to study the spectral asymptotics of classical and particular double-porosity operators. The reason for this shall be elucidated here. Moreover, at the end of this section we shall argue when quasi-periodic convergence is necessary via a model problem that we later study in detail in this article.

2.3.1 Classical homogenisation

Consider the resolvent problem: For fixed $f\in L^{2}(\Omega)$ find $u_{\varepsilon}\in H^{1}_{0}(\Omega)$ such that

[TABLE]

where the symmetric matrix-valued function $a$ is $\square$-periodic, elliptic and bounded: $\exists\nu>0$ such that

[TABLE]

The following homogenisation theorem is classical.

Theorem 1** (Classical homogenisation theorem).**

Let $\varepsilon$ be a sequence with limit [math], and $f_{\varepsilon}\in L^{2}(\Omega)$ a sequence such that $f_{\varepsilon}$ weakly converges in $L^{2}(\Omega)$ to some $f_{0}$ as $\varepsilon$ tends to zero. Then $u_{\varepsilon}\in H^{1}_{0}(\Omega)$ the solution to (2.6), for $f=f_{\varepsilon}$, converges weakly in $H^{1}_{0}(\Omega)$ (and strongly in $L^{2}(\Omega)$) to $u_{0}\in H^{1}_{0}(\Omega)$ the solution to

[TABLE]

Here $a^{\rm hom}$ is the constant symmetric homogenised matrix determined by $a$:

[TABLE]

It is well-known that the homogenisation theorem implies the Hausdorff convergence of spectra (cf. [2, Section 2]):

[TABLE]

Let us study the quasi-periodic two-scale limits of $u_{\varepsilon}$.

Proposition 6**.**

Fix ${\boldsymbol{\theta}}\in(0,2\pi)^{d}$ and consider $f_{\varepsilon}\in L^{2}(\Omega)$ such that $f_{\varepsilon}\overset{2-{{\boldsymbol{\theta}}}}{\rightharpoonup}f_{0}$. Then, $u_{\varepsilon}\in H^{1}_{0}(\Omega)$ ${\boldsymbol{\theta}}$-quasi-periodic two-scale converges to zero; that is $u_{\varepsilon}\overset{2-{{\boldsymbol{\theta}}}}{\rightharpoonup}0$.

Remark 3*.*

Proposition 6 informs us that for the classical resolvent problem (2.6), the non-zero quasi-periodic micro-oscillations at leading order do not contribute to the spectral asymptotics. So one need only study the ${\boldsymbol{\theta}}=0$ quasi-periodic oscillations, i.e. the standard two-scale limit. It is well-known that the (periodic) two-scale limit coincides with the classical limit provided by Theorem 1, see for example [1, 30]. 2. 2.

The part of the spectrum corresponding the ${\boldsymbol{\theta}}$-quasi-periodic micro-oscillations, for ${\boldsymbol{\theta}}\neq\bf 0$, actually resides in an $\varepsilon^{-2}$ neighbourhood of infinity; this can be formally seen from the considerations of Section 2.1: for $a_{\varepsilon}$ independent of $\varepsilon$, the eigenvalues $\lambda_{\varepsilon}$ in (2.1) are clearly of the order $\varepsilon^{-2}$. To study such ‘high-frequency’ spectrum one can consider the re-scaled operator $\varepsilon^{2}A_{\varepsilon}$, that is consider coefficients of the form $a_{\varepsilon}=\varepsilon^{2}a$. The precise study of such high-frequency spectra was performed in [2] for a broader class of moderately contrasting locally periodic coefficients. Therein, the authors provide a rigorous description of the high-frequency spectral asymptotics in terms of non-trivial quasi-momenta ${\boldsymbol{\theta}}$. This was done by introducing an appropriate notion of “Bloch wave homogenisation”. For the reduced setting of (globally) periodic coefficients, the Bloch-wave operator-limits determined therin can readily be shown to be equivalent to the ${\boldsymbol{\theta}}$-quasi-periodic two-scale limits.

Proof of Proposition 6.

The sequence $f_{\varepsilon}$ weakly converges, cf. Proposition 3 part 1., and so is bounded. Multiplying (2.6) (for $f=f_{\varepsilon}$) and integrating over $\Omega$, and using the ellipticity of $a$, produces the a-priori bound

[TABLE]

Applying Proposition 4, we deduce that there exists $u\in L^{2}(\Omega;H^{1}_{\boldsymbol{\theta}}(\square))$ such that, up to a subsequence, the following convergences hold:

[TABLE]

Let us show $u=0$: $\boldsymbol{\nabla}u_{\varepsilon}$ is a bounded sequence and so $\varepsilon\boldsymbol{\nabla}u_{\varepsilon}$ strongly converges to zero in $L^{2}$. Therefore, by Proposition 3 part 2., we deduce $\boldsymbol{\nabla}_{y}u=0$. As $\square$ is connected it follows that $u$ is constant. Yet $u_{\varepsilon}\in H^{1}_{\boldsymbol{\theta}}(\square)$ and there are no non-trivial constant ${\boldsymbol{\theta}}$-quasi-periodic functions for ${\boldsymbol{\theta}}\neq 0$, see (1.1). Hence, $u=0$. ∎

2.3.2 Double-porosity model

Consider the resolvent problem: For fixed $f\in L^{2}(\Omega)$ find $u_{\varepsilon}\in H^{1}_{0}(\Omega)$ such that

[TABLE]

Here

[TABLE]

where $Q_{0}$ is a smooth compactly contained subset of $\square$ such that, for $Q_{1}:=\square\backslash\overline{Q_{0}}$, the periodic extension

[TABLE]

forms a connected set in $\mathbb{R}^{d}$. The functions $a_{i}$, $i=0,1$ are taken to be real-valued, elliptic and bounded on $Q_{i}$. The following homogenisation theorem is established in [30, Theorem 5.1].

Theorem 2**.**

Suppose $f=f_{\varepsilon}$ in the right-hand side of (2.7) two-scale converges to some $f_{0}$, that is $f_{\varepsilon}\overset{2-{{\boldsymbol{\theta}}}}{\rightharpoonup}f_{0}$ for ${\boldsymbol{\theta}}=0$. Then, the sequence of solutions $u_{\varepsilon}$ two-scale converges to $u_{0}(x,y)=u(x)+v(x,y)$, where $(u,v)$ belongs to

[TABLE]

and uniquely solves

[TABLE]

Here, $a^{\rm hom}_{dp}$ is the constant symmetric and positive homogenised matrix for perforated domains determined by $a_{1}$:

[TABLE]

This result informs us that $A_{\varepsilon}$ strongly two-scale converges to the operator $A_{\boldsymbol{0}}$, defined in $L^{2}(\Omega\times\square)$, associated to the above two-scale limit resolvent problem. Therefore, by Proposition 5 for ${\boldsymbol{\theta}}=\boldsymbol{0}$, the lower-semicontinuity of the spectral convergence is ensured. In fact, Zhikov proved in [30, Theorem 8.1], under the condition that $F_{1}$ is connected in $\mathbb{R}^{d}$, the stronger result

[TABLE]

Let us determine the strong resolvent quasi-periodic two-scale limits of $A_{\varepsilon}$.

Proposition 7**.**

Fix ${\boldsymbol{\theta}}\in(0,2\pi)^{d}$. Suppose $f=f_{\varepsilon}$ in the right-hand side of (2.7) such that $f_{\varepsilon}\overset{2-{{\boldsymbol{\theta}}}}{\rightharpoonup}f_{0}$ to some $f_{0}\in L^{2}(\Omega\times\square)$. Then, the sequence of solutions $u_{\varepsilon}$ ${\boldsymbol{\theta}}$-quasi-periodically two-scale converges to $v_{0}(x,y)\in L^{2}(\Omega;H^{1}_{0}(Q_{0}))$ the solution to

[TABLE]

Remark 4*.*

1. Note that $A_{\boldsymbol{\theta}}=\mathcal{A}$ is independent of ${\boldsymbol{\theta}}$ for ${\boldsymbol{\theta}}\neq\boldsymbol{0}$, and its spectrum is the point spectrum given by the operator whose action is $u\mapsto-{\rm div}(a_{0}\boldsymbol{\nabla}u)$ with domain $\{u\in H^{1}_{0}(Q_{0})\,|\,-{\rm div}(a_{0}\boldsymbol{\nabla}u)\in L^{2}(Q_{0})\}$.

2. It is easy to see $\mathcal{A}\subset A_{\boldsymbol{0}}$ (by noting that setting $u=\phi=0$ in (2.8) gives (2.9)) and so

[TABLE]

The set

[TABLE]

is the limit spectrum arriving from quasi-periodic micro-oscillations.

The restriction of the limit spectrum $\sigma(A_{\boldsymbol{0}})$ to $\bigcup_{{\boldsymbol{\theta}}\neq\boldsymbol{0}}\sigma(A_{{\boldsymbol{\theta}}})$ is achieved by considering the purely macro-scopic component $u(x)$ (of eigenfunctions) to be zero. For this reason, we coin this spectrum to be pure Bloch spectrum. In the simplified setting of double-porosity the pure Bloch spectrum is point spectrum (due to the fact $A_{\boldsymbol{\theta}}=\mathcal{A}$ is independent of ${\boldsymbol{\theta}}$ for ${\boldsymbol{\theta}}\neq\boldsymbol{0}$). In general, we expect this spectrum to have band-gap structure, and the gaps have only contracted to points here due to the geometric constraint that $F_{1}$ is connected in $\mathbb{R}^{d}$. This expectation is verified in Section 6.

3. Even though the strong resolvent ${\boldsymbol{\theta}}$-quasi-periodic limit of $A_{\varepsilon}$ exists, it has trivial dependence on ${\boldsymbol{\theta}}$, ${\boldsymbol{\theta}}\neq\boldsymbol{0}$ and more importantly is a restriction of the two-scale limit $A_{\boldsymbol{0}}$. Hence, one need only consider $A_{\boldsymbol{0}}$, and this explains why in this setting one is to expect Zhikov’s result $\lim_{\varepsilon}\sigma(A_{\varepsilon})=\sigma(A_{\boldsymbol{0}})$. In general, the limit $A_{\boldsymbol{0}}$ will not be sufficient to capture the full spectral asymptotics.

Proof of Proposition 7.

Let us consider $a_{1}$ (respect. $a_{0}$) to be extended by zero into $Q_{0}$ (respect $Q_{1}$), and consider $\nu>0$ to be the constant such that

[TABLE]

The solution $u_{\varepsilon}$ solves

[TABLE]

Setting $\phi=u_{\varepsilon}$ in the above variational problem and using the fact that $a_{1}+a_{0}\geq\nu$, we deduce the a-priori bound (for $\varepsilon\leq 1$)

[TABLE]

Additionally, we have the bound

[TABLE]

Indeed, $a_{1}\geq 0$ and

[TABLE]

By Proposition 4 it follows that, up to a discarded subsequence,

[TABLE]

for some $v_{0}\in L^{2}(\Omega;H^{1}_{\boldsymbol{\theta}}(Q))$.

Let us show that $v_{0}\in L^{2}(\Omega;H^{1}_{0}(Q_{0}))$. By Proposition 2 it follows that

[TABLE]

Yet $\sqrt{a_{1}}(\tfrac{x}{\varepsilon})\varepsilon\boldsymbol{\nabla}u_{\varepsilon}$ strongly converges to zero in $L^{2}(\Omega)$. Therefore $\sqrt{a_{1}}(y)\boldsymbol{\nabla}_{y}v_{0}=0$, which is equivalent to $\boldsymbol{\nabla}_{y}v_{0}=0$ on $Q_{1}$ (recall $a_{1}$ is positive on $Q_{1}$ and zero on $Q_{0}$). As $Q_{1}$ is connected then $v_{0}$ is constant in $Q_{1}$. Now, since the periodic extension $F_{1}=\bigcup_{z\in\mathbb{Z}^{d}}(Q_{1}+z)$ forms a connected set then $v_{0}$ is constant in $F_{1}$. Yet, $v_{0}\in H^{1}_{\boldsymbol{\theta}}(\square)$ for ${\boldsymbol{\theta}}\neq\boldsymbol{0}$ and consequently this constant is zero, cf (1.1).

It remains to prove $v_{0}$ solves (2.9). This can easily be deduced by passing the the ${\boldsymbol{\theta}}$-quasi-periodic limit in (2.10) for test functions $\phi(x)=\psi(x)\varphi(\tfrac{x}{\varepsilon})$, $\psi\in H^{1}_{0}(\Omega)$, $\varphi\in H^{1}_{0}(Q_{0})$ and using convergences (2.11). ∎

2.3.3 An example with non-trivial quasi-periodic limits

Let us provide an example which demonstrates that in general the family $A_{{\boldsymbol{\theta}}}$, of strong resolvent ${\boldsymbol{\theta}}$-quasi-periodic limits to $A_{\varepsilon}$, are not restrictions of $A_{\boldsymbol{0}}$.

Suppose, we consider (2.7) for coefficients

[TABLE]

and $a_{i}$ are real-valued $\square$-functions such that $a_{1}\geq 0$ and $a_{1}+a_{0}\geq\nu>0$. Let $f_{\varepsilon}$ be a bounded sequence and $u_{\varepsilon}$ solve (2.7) for $f=f_{\varepsilon}$. Arguing as in the proof of Proposition 7, we see that $u_{\varepsilon}$ to solution to (2.7) will satisfy the a-priori bounds

[TABLE]

In particular, Proposition 4 informs us that up to a subsequence

[TABLE]

for some $u_{0}\in L^{2}(\Omega;H^{1}_{\boldsymbol{\theta}}(Q))$. Moreover, by an application of Proposition 2 we deduce that

[TABLE]

Denoting by $V_{\boldsymbol{\theta}}$ the closed linear subspace of $H^{1}_{\boldsymbol{\theta}}(\square)$ given by555The space $V_{\boldsymbol{0}}$ was first introduced in [29] and coined the space of microscopic oscillations. We have appropriately extended this notion to ${\boldsymbol{\theta}}$-quasi-periodic oscillations here, ${\boldsymbol{\theta}}\neq\boldsymbol{0}$.

[TABLE]

Suppose we show an example where $V_{{\boldsymbol{\theta}}}$ is not a subset of $V_{\boldsymbol{0}}$ for some non-trivial $V_{\boldsymbol{\theta}}$, ${\boldsymbol{\theta}}\neq\boldsymbol{0}$. Then, for such examples we should not expect that $A_{\boldsymbol{0}}$ is an extension of $A_{{\boldsymbol{\theta}}}$, nor should we expect $\sigma(A_{\boldsymbol{\theta}})\subset\sigma(A_{\boldsymbol{0}})$. Let us provide such an example. The conjectures (stated immediately above) based on this example will be proved rigorously in the remainder of the article.

Suppose $Q_{1}$ is the cylinderical domain

[TABLE]

and

[TABLE]

Notice that $F_{1}:=\bigcup_{z\in\mathbb{Z}^{3}}(Q_{1}+z)$ the periodic extension of $Q_{1}$ into $\mathbb{R}^{d}$ consists of infinitely many mutually disjoint cylinders $C_{l}:=\mathbb{R}\times[\tfrac{1}{4}+l_{1},\tfrac{3}{4}+l_{2}]^{2}$, $\forall l\in\mathbb{Z}^{2}$. That is, the assumptions of Subsection 2.3.2, and in particular [30, 31], do not hold.

Now $v\in V_{\boldsymbol{\theta}}$ if, and only if, $v\in H^{1}_{\boldsymbol{\theta}}(\square)$ with $v=c$ for some $c\in\mathbb{C}$ on $Q_{1}$. This implies, cf (1.1), that

[TABLE]

(in the sense of trace). Then, for $(y_{2},y_{3})\in(\tfrac{1}{4},\tfrac{3}{4})^{2}$ we arrive at the condition

[TABLE]

Therefore, if $\theta_{1}\neq 0$, then the above condition only holds if $c=0$. That is $v$ must necessarily be zero on $Q_{1}$. On the other hand, if $\theta_{1}=0$ then any $H^{1}_{\boldsymbol{\theta}}(\square)$ function that is constant on $Q_{1}$ belongs to $V_{\boldsymbol{\theta}}$. In particular we see that $V_{\boldsymbol{\theta}}$ does not belong to $V_{\boldsymbol{0}}$ for all ${\boldsymbol{\theta}}\in(0,2\pi)^{3}$.

Remark 5*.*

Note that if $Q_{1}$ contains a connected subset which joins two opposite faces of the square $\square$ then the space $V_{{\boldsymbol{\theta}}}$ non-trivially depends on ${\boldsymbol{\theta}}$. Consequently, non-trivial limit Bloch spectrum is expected for ${\boldsymbol{\theta}}$ aligned orthogonally to these faces.

The remainder of the article is dedicated to determining the strong quasi-periodic two-scale limits of $A_{\varepsilon}$ for such fibre-like inclusions. Moreover, we demonstrate that indeed $A_{\boldsymbol{\theta}}$ are not restrictions of $A_{\boldsymbol{0}}$ and that $\sigma(A_{\boldsymbol{\theta}})$ form non-trivial subsets of $\lim_{\varepsilon}\sigma(A_{\varepsilon})$.

3 Problem formulation and Homogenisation

In this article we are concerned with the asymptotic analysis of the resolvent problem

[TABLE]

where $\varepsilon<1$ is a small parameter, $\Omega$ is a smooth open bounded star-shaped domain666All the results and proofs follow through in an identical manner for the case where $\Omega$ is the whole space. and $f_{\varepsilon}\in L^{2}(\Omega)$ known. The coefficient $a_{\varepsilon}$ is given by

[TABLE]

and the regions $Q_{0}$ and $Q_{1}$ are described as follows (cf. Figure 2).

Geometric assumptions. For $j=1,2,3$ we consider smooth domains $S_{j}$ compactly contained in $(0,1)^{2}$ that have mutually disjoint closures. We denote by $C_{j}$ the cylinder aligned to the $j$-th co-ordinate axis with cross-section $S_{j}$, i.e. C_{1}:=\{y\in(0,1)^{3}\,\big{|}\,y\in(0,1)\times S_{1}\}, C_{2}:=\{y\in(0,1)^{3}\,\big{|}\,y=(z_{2},z_{3},z_{1}),z\in(0,1)\times S_{2}\} and C_{3}:=\{y\in(0,1)^{3}\,\big{|}\,y=(z_{3},z_{1},z_{2}),z\in(0,1)\times S_{3}\}.

Then, for a given non-empty subset $\mathcal{I}$ of $\{1,2,3\}$, we consider $Q_{1}=\cup_{i\in\mathcal{I}}C_{i}$. We denote by $\Gamma_{i}=\partial C_{i}\backslash\partial Q$ and $\Gamma=\bigcup_{i\in\mathcal{I}}\Gamma_{i}=\partial Q_{1}\backslash\partial Q$.

Under this geometric assumption we determine for each ${\boldsymbol{\theta}}\in[0,2\pi)^{3}$ the strong resolvent ${\boldsymbol{\theta}}$-quasi-periodic two-scale limit, cf. Section 2, of the self-adjoint operator $A_{\varepsilon}$ associated to resolvent problem (3.1). That is, for a fixed ${{\boldsymbol{\theta}}}\in[0,2\pi)^{3}$ and a given bounded sequence $f_{\varepsilon}\in L^{2}(\Omega)$ such that $f_{\varepsilon}\overset{2-{{\boldsymbol{\theta}}}}{\rightharpoonup}f$, i.e.

[TABLE]

we aim to determine the ${\boldsymbol{\theta}}$-quasi-periodic two-scale limit behaviour of the solution $u_{\varepsilon}\in H^{1}_{0}(\Omega)$ to

[TABLE]

As $f_{\varepsilon}$ is bounded in $L^{2}(\Omega)$, upon setting $\phi=u_{\varepsilon}$ in (3.3) we deduce that the sequences

[TABLE]

are bounded. Let us describe the ${\boldsymbol{\theta}}$-quasi-periodic two-scale limit, referring to Section 4 for the details.

The limit of $u_{\varepsilon}(x)$ will be a function $u(x,y)$, of two variables $x\in\Omega$, $y\in Q$, that is ${\boldsymbol{\theta}}$-quasi-periodic with respect to the second variable $y$, cf Proposition 4. Furthermore, due to the fact that in each cylinder $C_{i}$, $i\in\mathcal{I}$, the gradient of $u_{\varepsilon}$ is bounded, the limit $u$ necessarily belongs to the (Bochner) space $L^{2}(\Omega;V_{\boldsymbol{\theta}})$ where

[TABLE]

It follows from this (see (1.1)) that $u$ is non-zero in cylinder $C_{i}$ if and only if the $i$-th component $\theta_{i}$ of ${\boldsymbol{\theta}}$ is zero. If $\theta_{i}=0$, then we determine that $u_{i}$ is not only non-trivial but it is actually more regular in the $x_{i}$-th coordinate direction: $\partial_{x_{i}}u_{i}\in L^{2}(\Omega)$.

More precisely, for $\mathcal{I}^{\boldsymbol{\theta}}$ the subset of indexes $\mathcal{I}\subseteq\{1,2,3\}$ given by $\mathcal{I}^{\boldsymbol{\theta}}:=\{i\in\mathcal{I}\,|\,\theta_{i}=0\}$, we denote by $\mathbb{C}^{\boldsymbol{\theta}}$ the closed subspace of $\mathbb{C}^{3}$ spanned by $\{{\bf{e}}_{i}\}_{i\in\mathcal{I}^{\boldsymbol{\theta}}}$,777Note that $\mathbb{C}^{\boldsymbol{\theta}}$ is either the whole space, a plane or a line in $\mathbb{C}^{3}$. and show that the function $u$ belongs to the set

[TABLE]

which is clearly a Hilbert space when endowed with the inner product

[TABLE]

Here, $\boldsymbol{\nu}$ is the outer unit normal to $\partial\Omega$.

For each fibre $C_{i}$ there corresponds an effective constant material parameter $a^{\rm hom}_{i}>0$ given by

[TABLE]

where N^{(i)}\in H^{1}_{\#_{i}}(C_{i}):=\{u\in H^{1}(C_{i})\,|\,u\text{ is 1$-periodic in the variable$y_{i}}\} is the unique solution to the cell problem

[TABLE]

Then, for each ${\boldsymbol{\theta}}\in[0,2\pi)^{3}$, the ${\boldsymbol{\theta}}$-quasi-periodic two-scale limit problem is formulated as follows: For $f\in L^{2}(\Omega\times Q)$ find $u\in U_{\boldsymbol{\theta}}$ such that

[TABLE]

As $a^{\rm hom}_{i}$ are positive numbers and $a_{0}^{-1}\in L^{\infty}(Q_{0})$ it follows that the left-hand side of the above problem defines an equivalent inner product on the space $U_{\boldsymbol{\theta}}$, and consequently the existence and uniqueness of solutions $u$ to (3.9) are ensured by the Riesz representation theorem.

Setting $\phi=0$ on $Q_{1}$ in (3.9) gives the equation

[TABLE]

and a subsequent integration by parts in (3.9) leads to the variational formula

[TABLE]

for all $\boldsymbol{\phi}\in L^{2}(\Omega,\mathbb{C}^{\boldsymbol{\theta}})$, such that $\partial_{i}\phi_{i}\in L^{2}(\Omega)$. For each fixed $j\in\mathcal{I}^{\boldsymbol{\theta}}$ we set $\phi_{j}=\phi$, $\phi\in C^{\infty}_{0}(\Omega)$ and $\phi_{i}=0$ for $i\neq j$, above. This leads to the ${\boldsymbol{\theta}}$-quasi-periodic two-scale homogenised system of equations.

[TABLE]

Here

[TABLE]

for ${\bf{n}}$ the outer unit normal of $\Gamma_{j}=\partial C_{j}\backslash\partial Q$. We now state the main result of the article.

Theorem 3**.**

Consider $f_{\varepsilon}\in L^{2}(\Omega)$, $f\in L^{2}(\Omega\times Q)$ such that $f_{\varepsilon}\overset{2-{{\boldsymbol{\theta}}}}{\rightharpoonup}f$, and $u_{\varepsilon}$ the solution to (3.3). Then $u_{\varepsilon}$ converges, up to some subsequence, in the ${{\boldsymbol{\theta}}}$-quasi-periodic sense to $u\in U_{\boldsymbol{\theta}}$ the unique solution to (3.9), equivalently (3.10).

An immediate consequence of Theorem 3 is that for each ${\boldsymbol{\theta}}\in[0,2\pi)^{3}$, the operator $A_{\varepsilon}$ strong resolvent ${\boldsymbol{\theta}}$-quasi-periodically two-scale converges to the operator $A^{\rm hom}_{\boldsymbol{\theta}}$ associated to problem (3.10), see Definition 2 in Section 2. Consequently, Proposition 5 informs us that the lower semi-continuity of the spectra in the Hausdorff sense is ensured:

[TABLE]

The structure of the limit spectrum $\bigcup\limits_{{\boldsymbol{\theta}}\in[0,2\pi)^{3}}\sigma(A^{\rm hom}_{\boldsymbol{\theta}})$ is analysed in Section 6 and described in Proposition 10.

Remark 6*.*

A seperate issue, not explored here, is the so-called spectral completeness statement, i.e. the question of whether or not the remaining criterion for Hausdorff convergence of spectra is satisfied: does it follow that

[TABLE]

In general this will not be true due to the presence of the boundary, and the fact that $Q_{0}$ intersects the boundary. This leads to the expectation that there exists non-trivial spectrum due to surface waves asymptotically localised near the boundary, cf. [2] for analoguous results in the context of classical locally periodic media. For the case of $\Omega$ being the torus or the whole space the above assertion is expected to hold and will be explored in future works.

4 Proof of the homogenisation theorem

This section is dedicated to the proof of Theorem 3. To do this, we shall develop an appropriate quasi-periodic two-scale variation of a powerful method first introduced in [21] in the context of standard (periodic) two-scale convergence, i.e. ${{\boldsymbol{\theta}}}$-quasi-periodic two-scale convergence for ${{\boldsymbol{\theta}}}=0$. In what follows $\phi$, $\varphi$ and $\boldsymbol{\Phi}$ will denote respectively fixed arbitrary elements of $C^{\infty}_{0}(\Omega)$, $C^{\infty}_{{\boldsymbol{\theta}}}(Q)$ and $C^{\infty}_{{\boldsymbol{\theta}}}(Q;\mathbb{C}^{3})$.

4.1 Technical preliminaries

The following results will be of importance in the proof of the homogenisation theorem.

Lemma 1**.**

Let $B$ be the closure of a smooth domain and let $B_{1}$ be a smooth bounded domain such that $B\subset B_{1}$ and $A=B_{1}\backslash B$ is a connected set. Then every u\in H^{1}\big{(}(0,1)\times A\big{)} can be extended to $(0,1)\times B_{1}$ as a function \widetilde{u}\in H^{1}\big{(}(0,1)\times B_{1}\big{)} such that

[TABLE]

where $c$ does not depend on $u\in H^{1}((0,1)\times A)$.

Proof.

Suppose u\in H^{1}\big{(}(0,1)\times A\big{)}. Then, by Fubini’s theorem, for almost every $x_{1}\in(0,1)$ the function $u(x_{1},\cdot)$ belongs to $H^{1}(A)$ and let $Eu(x_{1},\cdot)$ be the Sobolev extension of $u(x_{1},\cdot)$ into $B_{1}$ given in [19, Lemma 3.2, pg. 88]. In particular, one has

[TABLE]

where $c$ does not depend on $u$ nor $x_{1}$. Here, $\boldsymbol{\nabla}^{\prime}$ denotes the gradient vector $(0,\partial_{x_{2}},\partial_{x_{3}})$.

Consider $\widetilde{u}$ given by

[TABLE]

Then $\boldsymbol{\nabla}^{\prime}\widetilde{u(x_{1},\cdot)}=\boldsymbol{\nabla}^{\prime}Eu(x_{1},\cdot)$ and from assertion (4.2) it follows that

[TABLE]

To prove (4.1), it remains to demonstrate that \partial_{x_{1}}\widetilde{u}\in L^{2}\big{(}(0,1)\times B_{1}\big{)} and

[TABLE]

For each $t\in\mathbb{R}$, the difference quotient is given by

[TABLE]

where we have extended $\widetilde{u}$ trivially by zero into $\mathbb{R}\backslash(0,1)$. Notice that $D_{t}\widetilde{u}=ED_{t}u$, i.e. the extension into $B_{1}$ of the function $\tfrac{u(x_{1}+t,\cdot)-u(x_{1},\cdot)}{t}=D_{t}u(x_{1}+t,\cdot)$, and consequently

[TABLE]

Since $D_{t}u$ converges strongly in L^{2}\big{(}(0,1)\times A\big{)} to $\partial_{x_{1}}u$, it follows that $D_{t}\widetilde{u}$ is a Cauchy sequence in L^{2}\big{(}(0,1)\times B_{1}\big{)} and this limit can be identified, using the fact that $(D_{t}\widetilde{u},\phi)_{L^{2}}=(\widetilde{u},D_{-t}\phi)_{L^{2}}$ for \phi\in C^{\infty}_{0}\big{(}(0,1)\times B_{1}\big{)}, as $\partial_{x_{1}}\widetilde{u}$. Furthermore, passing to the limit in the above inequality yields (4.4). ∎

Proposition 8**.**

Fix ${{\boldsymbol{\theta}}}\in[0,2\pi)^{3}$. There exists a constant $C=C_{{\boldsymbol{\theta}}}>0$ such that

[TABLE]

Here V_{{\boldsymbol{\theta}}}={\rm ker}\sqrt{a_{1}}\boldsymbol{\nabla}_{{\boldsymbol{\theta}}}=\{u\in H^{1}_{{\boldsymbol{\theta}}}(Q)\,\big{|}\,\sqrt{a_{1}}\boldsymbol{\nabla}u\equiv 0\}.

Proof.

For each $i\in\mathcal{I}$, let $S_{i}$ be the cross-section of the cylinder $C_{i}$. Since $S_{i}$ are compactly contained in $(0,1)^{2}$ and have mutually disjoint closures then there exists open $A_{i}$ such that $\overline{{S_{i}}}\subset A_{i}\subset(0,1)^{2}$ and $A_{i}$ are mutually disjoint. Let $\chi_{i}\in C^{\infty}_{0}(A_{i})$ be smooth cut-off functions that are identity on $S_{i}$, we extend $\chi_{i}$ by zero to $(0,1)^{2}$.

Now using Lemma 1, let $\widetilde{u_{i}}$ be the extension of $u|_{C_{i}}\in H^{1}(C_{i})$ to $H^{1}(D_{i})$, where $D_{i}$ is the cylinder whose axis is parallel to $y_{i}$ with cross-section $A_{i}$. Note that since $u$ is $\theta_{i}$-quasi-periodic in the variable $y_{i}$, then the extension will be also, see (4.3). By Lemma 1 it follows that

[TABLE]

For $\theta_{i}\neq 0$, the following Poincaré inequality

[TABLE]

holds888This follows from noting the lower bound on the spectrum of the laplacian on the space of $H^{1}{\boldsymbol{\theta}}(D_{i})$ functions that are $\theta_{i}$-quasi-periodic in direction $y_{i}$.. For $\theta_{i}=0$, one has

[TABLE]

for some $C>0$. Here $\langle\widetilde{u_{i}}\rangle:=\tfrac{1}{D_{i}}\int_{D_{i}}\widetilde{u_{i}}$.

Recalling, $\mathcal{I}^{{\boldsymbol{\theta}}}=\{i\in\mathcal{I}\,|\,\theta_{i}=0\}$, we set $\widetilde{u}=\sum_{i\in\mathcal{I}^{\boldsymbol{\theta}}}\chi_{i}(\widetilde{u_{i}}-\langle\widetilde{u_{i}}\rangle)+\sum_{i\in\mathcal{I}\backslash\mathcal{I}^{\boldsymbol{\theta}}}\chi_{i}\widetilde{u_{i}}$, here $\chi_{i}$ are taken to be constant in the variable $y_{i}$ and as above the complementary directions. It follows that $\widetilde{u}\in H^{1}_{\boldsymbol{\theta}}(Q)$ and $u-\widetilde{u}\in V_{\boldsymbol{\theta}}$. Note that, by construction and inequalities (4.6), and (4.7), one has

[TABLE]

Now, the positivity of $a_{1}$ on $Q_{1}$ and (4.5) imply that the element $v:=u-\widetilde{u}$ of $V_{\boldsymbol{\theta}}$ is such that

[TABLE]

and the result follows.

∎

4.2 Proof of Theorem 3

Consider the sequence $u_{\varepsilon}$ of solutions to (3.3), i.e.

[TABLE]

for $f_{\varepsilon}\overset{2-{{\boldsymbol{\theta}}}}{\rightharpoonup}f$, and recall, cf. (3.4), that

[TABLE]

Consequently, Proposition 4 informs us that a subsequence of $u_{\varepsilon}$ ${\boldsymbol{\theta}}$-quasi-periodic two-scale converges to some $u\in L^{2}(\Omega;H^{1}_{\boldsymbol{\theta}}(Q))$, and moreover $\varepsilon\boldsymbol{\nabla}u_{\varepsilon}\overset{2-{{\boldsymbol{\theta}}}}{\rightharpoonup}\boldsymbol{\nabla}_{y}u$. Let us study the structure of this limit $u$ further.

We begin by introducing the densely defined unbounded linear operator $\sqrt{a_{1}}\boldsymbol{\nabla}_{{\boldsymbol{\theta}}}:H^{1}_{{\boldsymbol{\theta}}}(Q)\subset L^{2}(Q)\rightarrow L^{2}(Q;\mathbb{C}^{3})$ which is given by the action

[TABLE]

We now argue that a generalised Weyl’s decomposition holds, which was first introduced and proved for the case ${{\boldsymbol{\theta}}}=0$ in [21].

Lemma 2**.**

Let $(\sqrt{a_{1}}\boldsymbol{\nabla}_{\boldsymbol{\theta}})^{*}$ denote the adjoint of $\sqrt{a_{1}}\boldsymbol{\nabla}_{\boldsymbol{\theta}}$. Then, the orthogonal decomposition

[TABLE]

holds.

Remark 7*.*

Lemma 2 is a generalisation of the well-known fact that (periodic) divergence-free vector fields are mutually orthogonal to gradients of (periodic) potentials in $L^{2}$. In fact, this classical result can be deduced from the above lemma by (formally)999In fact, as expected the proof of this statement for $I$ is much easier as $I$ is positive where as $\sqrt{a_{1}}$ is non-negative. setting $\sqrt{a_{1}}=I$ on $\square$.

Proof of Lemma 2.

By the Banach closed ranged theorem, this result will follow if we demonstrate that the range of $\sqrt{a_{1}}\boldsymbol{\nabla}_{{\boldsymbol{\theta}}}$ is closed, and this fact is implied by Proposition 8.

Indeed, suppose $u_{n}\in{\rm Ran}(\sqrt{a_{1}}\boldsymbol{\nabla}_{{\boldsymbol{\theta}}})$ converges strongly in $L^{2}(Q;\mathbb{C}^{3})$ to some $u$ as $n\rightarrow\infty$, i.e. there exists $w_{n}\in H^{1}_{{\boldsymbol{\theta}}}(Q)$ such that $\sqrt{a_{1}}\boldsymbol{\nabla}w_{n}$ converges strongly in $L^{2}(Q;\mathbb{C}^{3})$ to $u$. By Proposition 8, the sequence $w_{n}^{\perp}$, where $w_{n}^{\perp}$ denotes the orthogonal projection of $w_{n}$ onto the orthogonal complement $V^{\perp}_{{\boldsymbol{\theta}}}$ of $V_{{\boldsymbol{\theta}}}$ in $H^{1}_{{\boldsymbol{\theta}}}(Q)$, is a Cauchy sequence in $H^{1}_{{\boldsymbol{\theta}}}(Q)$ and therefore converges, up to some subsequence, to $w\in H^{1}_{{\boldsymbol{\theta}}}(Q)$. In particular, $\sqrt{a_{1}}\boldsymbol{\nabla}w_{n}=\sqrt{a_{1}}\boldsymbol{\nabla}w^{\perp}_{n}$ converges strongly in $L^{2}(Q)$ to $\sqrt{a_{1}}\boldsymbol{\nabla}w$ and, consequently $u=\sqrt{a_{1}}\boldsymbol{\nabla}w$. Hence, the range of $\sqrt{a_{1}}\boldsymbol{\nabla}_{{\boldsymbol{\theta}}}$ is closed. ∎

Let us now describe $u$ in detail.

Lemma 3**.**

*The function $u$ belongs to the Bochner space $L^{2}(\Omega;V_{\boldsymbol{\theta}})$. *

Proof.

Recall that

[TABLE]

and so we aim to show that $\sqrt{a_{1}}\boldsymbol{\nabla}_{{\boldsymbol{\theta}}}u=0$.

On the one hand we deduce from (4.9) and (2.2) that

[TABLE]

Yet, on the other hand, Proposition 2 and the assertion $\varepsilon\boldsymbol{\nabla}u_{\varepsilon}\overset{2-{{\boldsymbol{\theta}}}}{\rightharpoonup}\boldsymbol{\nabla}_{y}u$ imply

[TABLE]

Therefore, as finite sums of $\phi(x)\boldsymbol{\Phi}(y)$ are dense in $L^{2}(\Omega\times Q;\mathbb{C}^{3})$ it follows that $a_{1}\boldsymbol{\nabla}_{y}u=0$ and since $\sqrt{a_{1}}^{-1}\in L^{\infty}(Q)$ we find that $u\in L^{2}(\Omega;V_{{\boldsymbol{\theta}}})$. ∎

The following result is of fundamental importance in characterising the (${\boldsymbol{\theta}}$-quasi-periodic) limit of the flux $a_{1}(\tfrac{\cdot}{\varepsilon})\boldsymbol{\nabla}u_{\varepsilon}$ in terms of the limit $u$ of the function $u_{\varepsilon}$. Put another way, this identity is crucial for determining the homogenised coefficients.

Lemma 4**.**

There exists $\boldsymbol{\xi}\in L^{2}(\Omega\times Q;\mathbb{C}^{3})$ such that, up to a subsequence, $\sqrt{a_{1}}(\tfrac{\cdot}{\varepsilon})\boldsymbol{\nabla}u_{\varepsilon}\overset{2-{{\boldsymbol{\theta}}}}{\rightharpoonup}\boldsymbol{\xi}$. Moreover, $\boldsymbol{\xi}$ belongs to the Bochner space L^{2}\big{(}\Omega;{\rm ker}\big{(}(\sqrt{a_{1}}\boldsymbol{\nabla}_{{\boldsymbol{\theta}}})^{*}\big{)}\big{)} and the pair $(u,\boldsymbol{\xi})$ satisfies the identity

[TABLE]

Proof.

By Proposition 1 and (4.9) there exists $\boldsymbol{\xi}\in L^{2}(\Omega\times Q;\mathbb{C}^{3})$ such that, up to a subsequence that we discard, one has

[TABLE]

To prove \boldsymbol{\xi}\in L^{2}\big{(}\Omega;{\rm ker}\big{(}(\sqrt{a_{1}}\boldsymbol{\nabla}_{{\boldsymbol{\theta}}})^{*}\big{)}, we take in (4.8) test functions of the form $\varepsilon\phi(x)\varphi(\tfrac{x}{\varepsilon})$, $\phi\in C^{\infty}_{0}(\Omega)$, $\varphi\in C^{\infty}_{{\boldsymbol{\theta}}}(Q)$, and use (4.9), (4.11) to pass to the limit in $\varepsilon$ and deduce that

[TABLE]

Therefore, for almost every $x\in\Omega$ one has

[TABLE]

and, hence by Lemma 2 it follows that \boldsymbol{\xi}(x,y)\in L^{2}\big{(}\Omega;{\rm ker}\big{(}(\sqrt{a_{1}}\boldsymbol{\nabla}_{{\boldsymbol{\theta}}})^{*}\big{)}.

Let us now prove assertion (4.10). Henceforth, we consider \Psi\in{\rm ker}\big{(}(\sqrt{a_{1}}\boldsymbol{\nabla}_{{\boldsymbol{\theta}}})^{*}\big{)} to be ${{\boldsymbol{\theta}}}$-quasi-periodically extended to $\mathbb{R}^{3}$. We shall prove below the following “integration by parts” formula:

[TABLE]

Using Proposition 2, (4.11) and the convergence $u_{\varepsilon}\overset{2-{{\boldsymbol{\theta}}}}{\rightharpoonup}u$, we pass to the limit in the above formula to readily arrive at (4.10).

To prove (4.12), it is sufficient to prove the following: for every $w\in H^{1}(\mathbb{R}^{3})$ one has

[TABLE]

Indeed, (4.12) follows from utilising (4.13) and the following facts: for $\phi\in C^{\infty}(\Omega)$ then $u_{\varepsilon}\phi$ belongs to $H^{1}_{0}(\Omega)$, as $u_{\varepsilon}\in H^{1}_{0}(\Omega)$, and can be trivially extended to $H^{1}(\mathbb{R}^{3})$, and that

[TABLE]

Let us now prove (4.13). For $Q^{(z)}_{\varepsilon}=\prod_{i=1}^{3}\varepsilon(z_{i},z_{i}+1)$, $z\in\mathbb{Z}^{3}$, it follows that

[TABLE]

where the last equality comes from the change of variables $x=\varepsilon(y+z)$ and recalling that $a_{1}(y)$ is periodic and $\Psi$ is ${{\boldsymbol{\theta}}}$-quasi-periodic. By noting, for $w\in H^{1}(\mathbb{R}^{3})$, that

[TABLE]

is an element of $H^{1}_{{\boldsymbol{\theta}}}(Q)$, and that

[TABLE]

the identity (4.13) follows. ∎

We are now ready to describe the properties of the macroscopic part of $u$ and express the flux $\boldsymbol{\xi}$ in terms of $u$.

Lemma 5**.**

Let $(u,\boldsymbol{\xi})$, $u\in L^{2}(\Omega;V_{\boldsymbol{\theta}})$ and \boldsymbol{\xi}\in L^{2}(\Omega;{\rm ker}\big{(}(\sqrt{a_{1}}\boldsymbol{\nabla}_{\boldsymbol{\theta}})^{*}\big{)}, be a pair which satisfies the identity (4.10). Then, $u\in U_{\boldsymbol{\theta}}$, see (3.6). That is, for every for $i\in\mathcal{I}^{{\boldsymbol{\theta}}}=\{i\in\mathcal{I}\,|\,\theta_{i}=0\}$, one has $u=u_{i}$ on $\Omega\times C_{i}$, where $\partial_{i}u_{i}\in L^{2}(\Omega)$ with $u_{i}\nu_{i}=0$ on $\partial\Omega$, for $\boldsymbol{\nu}$ the outer unit normal to $\partial\Omega$. Furthermore,

[TABLE]

Here, $N^{(i)}$ solve (3.8).

The following result immediately follows from the above lemma.

Proposition 9**.**

For every for $i\in\mathcal{I}^{{\boldsymbol{\theta}}}$, one has

[TABLE]

Here, $a^{\rm hom}$ are given by (3.7), i.e.

[TABLE]

for N^{(i)}\in H^{1}_{\#_{i}}(C_{i})=\{u\in H^{1}(C_{i})\,|\,u\text{ is 1$-periodic in the variable$y_{i}}\} is the unique solution to the cell problem

[TABLE]

Proof.

Equation (4.14) implies

[TABLE]

For each $j\in\{1,2,3\}\backslash\{i\}$, we set $\phi=y_{j}$ in (4.15) to determine that

[TABLE]

Hence, it follows that

[TABLE]

for almost every $x\in\Omega$. Finally, from (4.15) it follows that

[TABLE]

Then, the positivity of $a^{\rm hom}_{i}$ can be seen by the inequality

[TABLE]

and noting that the right-hand side of this inequality can not be zero for this would contradict the periodicity of $N^{(i)}$ in the $y_{i}$ variable. ∎

Proof of Lemma 5.

As $u\in L^{2}(\Omega;V_{{\boldsymbol{\theta}}})$, see Lemma 3, then $u$ is constant in each fibre $C_{i}$, $i\in\mathcal{I}$. Now if ${\theta}_{i}\neq 0$ then $u$ is necessarily zero in $C_{i}$. On the other hand, if ${\theta}_{i}=0$, i.e. $i\in\mathcal{I}^{{\boldsymbol{\theta}}}$, then $u(x,y)=u_{i}(x)$ for $x\in\Omega$, $y\in C_{i}$. That is, $u=u_{i}$ on $\Omega\times C_{i}$ for some ${\bf{u}}\in L^{2}(\Omega,\mathbb{C}^{\boldsymbol{\theta}})$, where we recall that $\mathbb{C}^{\boldsymbol{\theta}}$ is the closed subspace of $\mathbb{C}^{3}$ spanned by $\{{\bf{e}}_{i}\}_{i\in\mathcal{I}^{\boldsymbol{\theta}}}$.

Let us now demonstrate that $u$ belongs to the Hilbert space $U_{\boldsymbol{\theta}}.$ By substituting $u=u_{i}$ on $\Omega\times C_{i}$, $i\in\mathcal{I}$, into (4.10), we deduce that

[TABLE]

For fixed $j\in\mathcal{I}^{{\boldsymbol{\theta}}}$, we will show directly below that there exists a function \boldsymbol{\Psi}^{(j)}{\ker}\big{(}(\sqrt{a_{1}}\boldsymbol{\nabla}_{\boldsymbol{\theta}})^{*}\big{)} such that

[TABLE]

Therefore

[TABLE]

and consequently substituting $\boldsymbol{\Psi}^{(j)}$ into (4.16) gives

[TABLE]

That is, $\partial_{x_{j}}u_{j}(x)=\int_{Q}\boldsymbol{\xi}(x,y)\cdot\overline{\boldsymbol{\Psi}^{(j)}}\,{\rm d}y\in L^{2}(\Omega)$ and $u_{j}\nu_{j}=0$ on $\partial\Omega$ where $\boldsymbol{\nu}$ is the outer unit normal to $\partial\Omega$, i.e. $u\in U_{\boldsymbol{\theta}}$ if (4.17) holds.

To show (4.17), we note that under the geometric assumptions on cylinders $C_{i}$, $i\in\mathcal{I}$, there exists a function

[TABLE]

Then, for each $j\in\mathcal{I}^{\boldsymbol{\theta}}=\{i\in\mathcal{I}\,|\,\theta_{i}=0\}$, the function $\boldsymbol{\Psi}^{(j)}=\tfrac{1}{|C_{j}|\sqrt{a_{1}}}\chi_{j}{\bf{e}}_{j}$ clearly satisfies (4.17). Furthermore, $\boldsymbol{\Psi}^{(j)}$ belongs to {\ker}\big{(}(\sqrt{a_{1}}\boldsymbol{\nabla})^{*}\big{)}: Indeed, as ${\theta}_{j}=0$, an element $\phi\in H^{1}_{{\boldsymbol{\theta}}}(Q)$ is $1$-periodic in the variable $y_{j}$, and it follows

[TABLE]

Therefore, (4.17) holds.

Let us now demonstrate (4.14). For $i\in\mathcal{I}^{{\boldsymbol{\theta}}}$, and almost every $x\in\Omega$, notice that

[TABLE]

and, by the geometric assumption of the cylinders, we can extend $y_{j}$ into $Q$ such that the extensions are elements of $H^{1}_{{\boldsymbol{\theta}}}(Q)$ and equal to zero on $C_{j}$. Therefore, it follows that

[TABLE]

Consequently, (4.16) takes the form

[TABLE]

Integrating by parts above, which is permissible since $\partial_{i}u_{i}\in L^{2}(\Omega)$, we deduce that

[TABLE]

That is, for almost every $x$, $\boldsymbol{\xi}(x,\cdot)$ is the projection onto \ker\big{(}(\sqrt{a_{1}}\boldsymbol{\nabla}_{{\boldsymbol{\theta}}})^{*}\big{)} of the function

[TABLE]

For each $i\in\mathcal{I}^{{\boldsymbol{\theta}}}$, let $\chi_{i}$ given by (4.18), and we introduce $\widetilde{N^{i}}\in H^{1}(Q)$ the extension into $Q$, given by Lemma 1, of the function $N^{(i)}\in H^{1}_{\#_{i}}(C_{i})$ that solves (4.15). It follows that $\sum_{i\in\mathcal{I}^{{\boldsymbol{\theta}}}}\chi_{i}\widetilde{N^{(i)}}$ belongs to $H^{1}_{{\boldsymbol{\theta}}}(Q)$ and

[TABLE]

for all $\phi\in H^{1}_{{\boldsymbol{\theta}}}(Q)$. That is, $\sqrt{a_{1}}\chi_{i}[\boldsymbol{\nabla}\widetilde{N^{(i)}}+{\bf{e}}_{i}]$ belongs to \ker\big{(}(\sqrt{a_{1}}\boldsymbol{\nabla}_{{\boldsymbol{\theta}}})^{*}\big{)}. Obviously

[TABLE]

and

[TABLE]

since $\chi_{i}$ is piece-wise constant on $C$. Consequently, as $\boldsymbol{\xi}(x,\cdot)$ is the projection of $w(x,\cdot)$ onto \ker\big{(}(\sqrt{a_{1}}\boldsymbol{\nabla})^{*}\big{)}, we have

[TABLE]

Hence, (4.14) holds and the proof is complete.

∎

We now conclude with the proof of Theorem 3. That is, we show that $u$ solves (3.9). We being by stating that under the assumption that $\Omega$ is star-shaped, standard pull-back and mollification type arguments prove that functions smooth in $x$ are dense in the Hilbert space $U_{\boldsymbol{\theta}}$. Therefore, it is sufficient to show (3.9) holds for such test functions $\phi$. Let us take such a $\phi$ and consider the test functions $\phi_{\varepsilon}(x)=\phi(x,\tfrac{x}{\varepsilon})$, $x\in\Omega$ in (4.8). Utilising the convergences

[TABLE]

we pass to the limit $\varepsilon\rightarrow 0$ in (4.8) to deduce that

[TABLE]

Then, as $\phi=\phi_{i}$ on $\Omega\times C_{i}$, with $\phi\neq 0$ only if $i\in\mathcal{I}^{\boldsymbol{\theta}}$, Proposition 9 implies that

[TABLE]

and (3.9) follows.

5 Quasi-periodic two-scale limit operator

For ${\boldsymbol{\theta}}\in[0,2\pi)^{3}$, we consider the subspace $H$ which is the closure of $U_{\boldsymbol{\theta}}$ in $L^{2}(\Omega\times Q)$, i.e.

[TABLE]

Indeed, for $f\in H$, we have $f=f_{i}$ on $\Omega\times C_{i}$, $i\in\mathcal{I}$, and consequently we deduce that

[TABLE]

and therefore, $H$ is closed in $L^{2}(\mathbb{R}^{3};L^{2}(Q))$. It is also straightforward to show that $U_{\boldsymbol{\theta}}$ is dense in $H$. Defining on $U_{\boldsymbol{\theta}}$ the form

[TABLE]

we find that, since $a^{\rm hom}_{i}$ are positive constants and $a^{-1}_{0}\in L^{\infty}(Q_{0})$, $Q_{\boldsymbol{\theta}}$ is closed when considered as a form on $H$. Setting $A^{\rm hom}_{\boldsymbol{\theta}}:D(A^{\rm hom}_{\boldsymbol{\theta}})\subset H\rightarrow H$ to be the unbounded self-adjoint operator generated by $Q_{\boldsymbol{\theta}}$, for $f\in L^{2}(\mathbb{R}^{3};Q)$ the ${\boldsymbol{\theta}}$-quasi-periodic two-scale homogenised limit problem (3.9) takes the form $A^{\rm hom}_{\boldsymbol{\theta}}u=P_{\boldsymbol{\theta}}f$. Here, $P_{\boldsymbol{\theta}}:L^{2}(\mathbb{R}^{3};L^{2}(Q))\rightarrow H$ is the orthogonal projection given by

[TABLE]

An immediate consequence of Theorem 3 is that for each ${\boldsymbol{\theta}}\in[0,2\pi)^{3}$, the operator $A_{\varepsilon}$ strong ${\boldsymbol{\theta}}$-quasi-periodic two-scale resolvent converges to $A^{\rm hom}_{\boldsymbol{\theta}}$, see Section 2 definition 2.

5.1 Spatial operators

Introducing the notation

[TABLE]

we consider the Hilbert space

[TABLE]

endowed with the inner product

[TABLE]

and the following bilinear form defined on $H$:

[TABLE]

Note that for ${\boldsymbol{\theta}}$ such that $\mathcal{I}^{\boldsymbol{\theta}}=\{i\in\mathcal{I}\,|\,\theta_{i}=0\}=\emptyset$ then $H_{\boldsymbol{\theta}}$ is zero and for such ${\boldsymbol{\theta}}$ we define our ‘spatial’ operator $A_{\boldsymbol{\theta}}$ to be the zero map. Otherwise, $\alpha_{\boldsymbol{\theta}}$ is a positive form on $H_{\boldsymbol{\theta}}$ and therefore has a positive self-adjoint operator $A_{\boldsymbol{\theta}}$, densely defined in $L^{2}(\Omega;\mathbb{C}^{\boldsymbol{\theta}})$, associated with the form. The space $H$ is compactly embedded101010This follows from an application of Vitali’s theorem, which is permissible by noting that since $u_{i}$ has an $L^{2}$ weak derivative in the $x_{i}$-th direction one can use the fundamental theorem of calculus to prove that any bounded sequence in $H$ is 2-equi-integrable. into $L^{2}$, and consequently the spatial operator $A_{\boldsymbol{\theta}}$ has compact resolvent and therefore its spectrum is discrete.

5.2 Pure Bloch operators

Consider the space

[TABLE]

which is a closed subspace of $H^{1}_{{{\boldsymbol{\theta}}}}(Q)$, and therefore is a Hilbert space when equipped with standard $H^{1}_{{{\boldsymbol{\theta}}}}(Q)$ norm. Define the sesquilinear form

[TABLE]

Since $a_{0}$ is positive and bounded on $Q_{0}$, and elements of $\mathcal{V}_{{\boldsymbol{\theta}}}$ have zero trace on the part of the boundary $\Gamma=\partial Q_{1}$, then by Poincaré’s inequality the form $\beta_{{\boldsymbol{\theta}}}$ is (uniformly in ${{\boldsymbol{\theta}}}$) coercive and bounded on $\mathcal{V}_{{\boldsymbol{\theta}}}$, i.e. there exists $c_{1}$ and $c_{2}$ independent of ${{\boldsymbol{\theta}}}$ such that

[TABLE]

for all $u,v\in\mathcal{V}_{{\boldsymbol{\theta}}}$. This implies that for every $f\in L^{2}(Q_{0})$ there exists a unique solution $u\in\mathcal{V}_{{\boldsymbol{\theta}}}$ such that

[TABLE]

Consequently, the unbounded self-adjoint linear operator $B_{{\boldsymbol{\theta}}}$, defined in $L^{2}(Q_{0})$, given by $B_{{\boldsymbol{\theta}}}u=f$, is positive and, moreover, by the Rellich embedding theorem has compact resolvent. Therefore the spectrum of $B_{{\boldsymbol{\theta}}}$ is discrete, and we order the eigenvalues in accordance with the min-max principle. These eigenvalues can be shown to be continuous functions of ${{\boldsymbol{\theta}}}$, in fact the following result holds.

Lemma 6**.**

For each $n\in\mathbb{N}$, let $\mu^{(n)}_{{\boldsymbol{\theta}}}$ denote the $n$-th eigenvalue of $B_{{\boldsymbol{\theta}}}$ as ordered according to the min-max principle, i.e.

[TABLE]

where $v\perp v_{i}$ is shorthand for $v$ is orthogonal to $v_{i}$ in $L^{2}(Q_{0})$. Then, for each $n\in\mathbb{N}$ the function $\lambda_{n}({{\boldsymbol{\theta}}}):=\mu^{(n)}_{{\boldsymbol{\theta}}}$ is Lipschitz continuous, that is there exists a $C_{n}>0$ such that

[TABLE]

The proof relies on an important observation is that the spaces $\mathcal{V}_{{\boldsymbol{\theta}}}$, ${{\boldsymbol{\theta}}}\in[0,2\pi)^{3}$, are mutually isomorphic. Indeed, if ${{\boldsymbol{\theta}}},{{\boldsymbol{\theta}}}^{\prime}\in[0,2\pi)^{3}$ then it is clear that the isometric mapping $\mathcal{U}({{\boldsymbol{\theta}}},{{\boldsymbol{\theta}}}^{\prime}):L^{2}(Q)\rightarrow L^{2}(Q)$ defined as multiplication by the function \mathrm{exp}\big{(}\rm i({{\boldsymbol{\theta}}}^{\prime}-{{\boldsymbol{\theta}}})\cdot y\big{)} defines an isomorphism between $\mathcal{V}_{{{\boldsymbol{\theta}}}}$ and $\mathcal{V}_{{{\boldsymbol{\theta}}}^{\prime}}$.

Proof.

Let $v$ be $L^{2}(Q_{0})$-normalised element of $\mathcal{V}_{{{\boldsymbol{\theta}}}}$ and consider v^{\prime}:=\mathcal{U}({{\boldsymbol{\theta}}},{{\boldsymbol{\theta}}}^{\prime})v=\mathrm{exp}\big{(}\rm i({{\boldsymbol{\theta}}}^{\prime}-{{\boldsymbol{\theta}}})\cdot y\big{)}v. Then, $v^{\prime}$ is an $L^{2}(Q_{0})$-normalised element of $\mathcal{V}_{{{\boldsymbol{\theta}}}^{\prime}}$ and the following identity

[TABLE]

holds. Therefore, one has

[TABLE]

Consequently, as the isometric mapping $\mathcal{U}({{\boldsymbol{\theta}}},{{\boldsymbol{\theta}}}^{\prime}):L^{2}(Q_{0})\rightarrow L^{2}(Q_{0})$ is an isomorphism between $\mathcal{V}_{{\boldsymbol{\theta}}}$ and $\mathcal{V}_{{{\boldsymbol{\theta}}}^{\prime}}$, the above inequality and the min-max formula (5.2) implies that

[TABLE]

Now, if we consider the self-adjoint Dirichlet operator in $L^{2}(Q_{0})$ associated with the form

[TABLE]

then, since $H^{1}_{0}(Q_{0})$ is embedded in $\mathcal{V}_{{\boldsymbol{\theta}}}$ for all ${{\boldsymbol{\theta}}}$, one has

[TABLE]

Here $\mu_{n}$ is the $n$-th eigenvalue111111The spectrum of $B_{D}$ is discrete, which again is a consequence of the Rellich theorem. of the operator $B_{D}$, defined in a similar manner as $B_{{\boldsymbol{\theta}}}$ above. Hence, we deduce from (5.3) that $\lambda_{n}({{\boldsymbol{\theta}}})$ is Lipschitz continuous with a Lipschitz constant bounded from above by $2||a_{0}||^{1/2}_{L^{\infty}(Q_{0})}\mu_{n}.$ ∎

6 Quasi-periodic two-scale limit spectrum

In this section we study the spectrum

[TABLE]

In particular we shall characterise the spectrum in terms of the spatial and pure Bloch operators introduced in Section 5. This leads to an appropriate analogue of the Zhikov $\beta$ function, cf. [30].

Let us fix ${\boldsymbol{\theta}}\in[0,2\pi)^{3}$ and suppose that $\lambda_{{\boldsymbol{\theta}}}$ is in the spectrum of $A_{{\boldsymbol{\theta}}}^{\rm hom}$. Then, there exists an eigenfunction $u_{{\boldsymbol{\theta}}}\in V_{{\boldsymbol{\theta}}}$ that solves the spectral problem

[TABLE]

Here, we recall that

[TABLE]

There are two subcases to study: when ${{\boldsymbol{\theta}}}\in\cup_{i\in\mathcal{I}}\Pi_{i}$, for $\Pi_{i}:=\{{{\boldsymbol{\theta}}}\in[0,2\pi^{2})^{3}\,|\,{\boldsymbol{\theta}}\cdot{\bf{e}}_{i}=0\}$, and {{\boldsymbol{\theta}}}\in[0,2\pi)^{3}\backslash\big{(}\cup_{i\in\mathcal{I}}\Pi_{i}\big{)}.

Pure Bloch spectrum. If {{\boldsymbol{\theta}}}\in[0,2\pi)^{3}\backslash\big{(}\cup_{i\in\mathcal{I}}\Pi_{i}\big{)}, then $\lambda_{\boldsymbol{\theta}},u_{{\boldsymbol{\theta}}}$ solves the problem

[TABLE]

Therefore, setting $u_{{\boldsymbol{\theta}}}(x,y)=\phi(x)v_{{\boldsymbol{\theta}}}(y)$ for a sufficiently arbitrary $\phi$, we find that $v_{{\boldsymbol{\theta}}}$ solves

[TABLE]

Therefore, the spectrum of $A_{{\boldsymbol{\theta}}}^{\rm hom}$ for {{\boldsymbol{\theta}}}\in[0,2\pi)^{3}\backslash\big{(}\cup_{i\in\mathcal{I}}\Pi_{i}\big{)} consists of eigenvalues of infinite multiplicity, and these eigenvalues coincide with the eigenvalues the pure Bloch operator $B_{{\boldsymbol{\theta}}}$ introduced in Section 5.2. Lemma 6 implies that these eigenvalues are continuous with respect to ${{\boldsymbol{\theta}}}$, and by continuously extending ${{\boldsymbol{\theta}}}$ from [0,2\pi)^{3}\backslash\big{(}\cup_{i\in\mathcal{I}}\Pi_{i}\big{)} to $[0,2\pi)^{3}$ we deduce that

[TABLE]

It is for this reason that we call $\bigcup_{{{\boldsymbol{\theta}}}\in[0,2\pi)^{3}}\sigma(B_{{\boldsymbol{\theta}}})$ the pure Bloch spectrum of $A^{\rm hom}_{{\boldsymbol{\theta}}}$.

Spatial spectrum. Let us now suppose that ${{\boldsymbol{\theta}}}\in\cup_{i\in\mathcal{I}}\Pi_{i}$ and $\lambda_{\boldsymbol{\theta}}\in\sigma(A_{\boldsymbol{\theta}})$ is not a pure Bloch eigenvalue, i.e. $\lambda_{\boldsymbol{\theta}}\notin\bigcup_{{{\boldsymbol{\theta}}}\in[0,2\pi)^{3}}\sigma(B_{{\boldsymbol{\theta}}})$. Introducing, for $i\in\mathcal{I}^{\boldsymbol{\theta}}$ the functions $b^{(i)}_{{\boldsymbol{\theta}}}\in H^{1}_{{\boldsymbol{\theta}}}(Q)$ that satisfy

[TABLE]

we represent $u_{{\boldsymbol{\theta}}}$ as follows

[TABLE]

and substitute this representation into (6.1) to deduce that $v_{{\boldsymbol{\theta}}}(x,y)\in\mathcal{V}_{{\boldsymbol{\theta}}}$, see (5.1), solves

[TABLE]

Denoting respectively by $\mu^{(m)}_{{\boldsymbol{\theta}}}$ and $v^{(m)}_{{\boldsymbol{\theta}}}$ the $m$-th eigenvalue and orthonormal eigenfunction of $B_{{\boldsymbol{\theta}}}$, we perform a spectral decomposition of $v_{\boldsymbol{\theta}}$ and $b^{(i)}_{\boldsymbol{\theta}}$ to conclude that

[TABLE]

for some $c_{m}({\boldsymbol{\theta}},x)$ and constants $b^{(i)}_{m}({\boldsymbol{\theta}})=\int_{Q_{0}}b^{(i)}_{\boldsymbol{\theta}}\overline{v_{\boldsymbol{\theta}}^{(m)}}.$ Substituting the spectral representations into (6.5) gives

[TABLE]

Therefore, $u_{{\boldsymbol{\theta}}}$ admits the form

[TABLE]

Consequently, we calculate

[TABLE]

Recalling that $b^{(i)}_{\boldsymbol{\theta}}$ solves (6.4), $v_{\boldsymbol{\theta}}^{(m)}$ solves (6.3), and utilising Green’s identity, we deduce that

[TABLE]

Therefore

[TABLE]

and for each $i\in\mathcal{I}^{{\boldsymbol{\theta}}}$, $u_{i}$ solves the problem

[TABLE]

for

[TABLE]

Hence, we have demonstrated the following.

Proposition 10**.**

The spectrum of $\bigcup\limits_{{\boldsymbol{\theta}}\in[0,2\pi)^{3}}\sigma(A_{\boldsymbol{\theta}}^{\rm hom})$ is the union of the following two sets:

• •

The pure Bloch spectrum:

[TABLE]

where $\mu^{(m)}_{\boldsymbol{\theta}}$ are the eigenvalues of $B_{\boldsymbol{\theta}}$ ordered according to the min-max principle.

• •

The spatial spectrum: $\{\lambda\in\ [0,\infty)\,|\,$ $\gamma({\boldsymbol{\theta}})(\lambda)\in\sigma(A_{\boldsymbol{\theta}})\}$, where $A_{\boldsymbol{\theta}}$ is an operator with compact resolvent. Here $\gamma:\mathbb{R}^{3}\rightarrow\mathbb{S}^{3}$ is for each ${\boldsymbol{\theta}}$ a (possibly) sign-indefinite symmetric matrix defined by setting for $i\notin\mathcal{I}^{\boldsymbol{\theta}}$, $\gamma_{ij}({\boldsymbol{\theta}})=0$ for all $j$, and $\gamma_{ij}({\boldsymbol{\theta}})=\beta^{(ij)}_{\boldsymbol{\theta}}$ otherwise.

Acknowledgements

This work was performed under the financial support of the Engineering and Physical Sciences Research Council Grant EP/M017281/1: “Operator asymptotics, a new approach to length-scale interactions in metamaterials.”