The reflection principle in the control problem of the heat equation

TL;DR

This paper investigates how reflection symmetry in a domain affects the null-controllability of the heat equation, providing new insights and explicit control cost bounds for various geometric domains.

Contribution

It establishes that null-controllability on a symmetric domain implies controllability on its parts with bounded control costs, extending to specific geometric shapes.

Findings

Null-controllability on the whole domain implies controllability on symmetric parts.

Control cost on parts does not exceed that of the whole domain.

Explicit bounds for control costs on triangles and prisms.

Abstract

We consider the control problem for the generalized heat equation for a Schroedinger operator on a domain with a reflection symmetry with respect to a hyperplane. We show that if this system is null-controllable, then so is the system on its respective parts and the corresponding control cost does not exceed the one on the whole domain. As an application, we obtain null-controllability results for the heat equation on half-spaces, orthants, and sectors of angle $\pi/2$. As a byproduct, we also obtain explicit control cost bounds for the heat equation on certain triangles and corresponding prisms in terms of geometric parameters of the control set.

Equations

$${\theta }_1\Vert \xi {\Vert }_{{\mathbb{C}}^d}^2\le ⟨A(x)\xi ,\xi {⟩}_{{\mathbb{C}}^d}\le {\theta }_2\Vert \xi {\Vert }_{{\mathbb{C}}^d}^2\forall \xi \in {\mathbb{C}}^d,\text{ a.e.leavevmodenobreak }x\in \Omega ,$$

$$-\nabla \cdot (A\nabla )+V$$

$${\partial }_{t}u(t)+H_{\Omega }^{\bullet }u(t)={\chi }_{\omega }v(t)\text{ for }0<t<T,u(0)=u_0,$$

$${\mathrm{e}}^{-T{H}_{\Omega }^{\bullet }}{u}_0+{\int }_0^T{\mathrm{e}}^{-(T-s)H_{\Omega }^{\bullet }}{\chi }_{\omega }v(s)ds=0\text{ in }{L}^2(\Omega ).$$

$$C_T:=\underset{\Vert u_0{\Vert }_{L^2(\Omega )}=1}{\sup }\inf \{\Vert v{\Vert }_{L^2((0,T),L^2(\omega ))}:v\text{ satisfies }\text{eqref}eq:nullcontrol\}$$

$$\tilde{A}(x)=\left\{\begin{array}[]{ll}A(x),&x\in\Omega\\ U(A\circ M)(x)U,&x\in M(\Omega)\end{array}\right.,\qquad U:=\operatorname{diag}(-1,1,\ldots,1),$$

$$\tilde{V}(x)=\left\{\begin{array}[]{ll}V(x),&x\in\Omega\\ (V\circ M)(x),&x\in M(\Omega)\end{array}\right..$$

$$-\nabla \cdot (\tilde{A}\nabla )+\tilde{V}$$

$${\partial }_t\tilde{u}(t)+H_{\tilde{\Omega }}^{\bullet }\tilde{u}(t)={\chi }_{\tilde{\omega }}\tilde{v}(t)\text{ for }0<t<T,\tilde{u}(0)={\tilde{u}}_0,$$

$${\partial }_{t}u(t)-{\Delta }_{\Omega }^{\bullet }u(t)={\chi }_{\omega }v(t)\text{ for }0<t<T,u(0)=u_0,$$

$$|S\cap (x+[0,a_1]\times \dots \times [0,a_d])|\ge \gamma \prod _{j=1}^d{a}_j\forall x\in {\mathbb{R}}^d,$$

$$C_T\le \frac{1}{\sqrt{T}}{\left(\frac{K^d}{\gamma }\right)}^{Kd/2}\exp \left(\frac{K\Vert a{\Vert }_1^2{\ln }^2(K^d/\gamma )}{2T}\right),$$

$${\partial }_t\tilde{u}(t)-{\Delta }_{\tilde{\Omega }}^{\bullet }\tilde{u}(t)={\chi }_{\tilde{\omega }}\tilde{v}(t)\text{ for }0<t<T,\tilde{u}(0)={\tilde{u}}_0,$$

$$|\tilde{S}\cap (x+Q_{\tilde{a}})|\ge \frac{\gamma }{2}2a_1\prod _{j=2}^d{a}_j=\tilde{\gamma }|Q_{\tilde{a}}|\forall x\in {\mathbb{R}}^d.$$

$$|\tilde{S}\cap (x+Q_{\tilde{a}})|\ge |S^{\prime }\cap (y+Q_a)|=|S\cap (y+Q_a)|\ge \gamma \prod _{j=1}^d{a}_j=\tilde{\gamma }|Q_{\tilde{a}}|.$$

$$|\tilde{S}\cap (x+Q_{\tilde{a}})|\ge |M(S^{\prime })\cap M(y+Q_a)|=|S\cap (y+Q_a)|\ge \tilde{\gamma }|Q_{\tilde{a}}|,$$

$$C_T\le \frac{1}{\sqrt{T}}{\left(\frac{2K^d}{\gamma }\right)}^{Kd/2}\exp \left(\frac{K(2a_1+a_2+\dots +a_d{)}^2{\ln }^2(2K^d/\gamma )}{2T}\right),$$

$$C_T\le \frac{1}{\sqrt{T}}{\left(\frac{(2K{)}^d}{\gamma }\right)}^{Kd/2}\exp \left(\frac{4K\Vert a{\Vert }_1^2{\ln }^2((2K{)}^d/\gamma )}{2T}\right),$$

$$C_{T}\leq\frac{1}{\sqrt{T}}\Bigl{(}\frac{2^{3n-4}K^{2}}{\tilde{\gamma}}\Bigr{)}^{K}\exp\Bigl{(}\frac{2^{3n-4}K\|{\tilde{a}}\|_{1}^{2}\ln^{2}(2^{3n-4}K^{2}/\tilde{\gamma})}{2T}\Bigr{)},$$

$$\tilde{\gamma }=\frac{\gamma a_1{a}_2}{4(a_1^2+a_2^2)},\tilde{a}=(2\sqrt{a_1^2+a_2^2},2\sqrt{a_1^2+a_2^2}),$$

$$\frac{\tilde{\gamma }{\tilde{a}}_1{\tilde{a}}_2}{4({{\tilde{a}}_1}^2+{{\tilde{a}}_2}^2)}=\frac{\tilde{\gamma }}{8}\text{ and }(2\sqrt{{{\tilde{a}}_1}^2+{{\tilde{a}}_2}^2},2\sqrt{{{\tilde{a}}_1}^2+{{\tilde{a}}_2}^2})=2\sqrt{2}\tilde{a}.$$

$${\mathcal{T}}_L:=\{(x,y)\in {\Lambda }_L^2:y<x\}\text{ and }{\mathcal{P}}_L:={\mathcal{T}}_L\times (0,L).$$

$$C_T\le \frac{1}{\sqrt{T}}{\left(\frac{K^d}{\tilde{\gamma }}\right)}^{Kd/2}\exp \left(\frac{K\Vert \tilde{a}{\Vert }_1^2{\ln }^2(K^d/\tilde{\gamma })}{2T}\right),$$

$${\partial }_{t}u(t)+(-{\Delta }_{\Omega }^{\bullet }+V)u(t)={\chi }_{\omega }v(t)\text{ for }0<t<T,u(0)=u_0,$$

$$S\supset \bigcup _{j\in {\mathbb{Z}}^d}B(z_j,\delta )$$

$$C_T\le {\left(\frac{\delta }{G}\right)}^{-D\cdot (1+G^{4/3}\Vert V{\Vert }_{\infty }^{2/3})}\frac{D}{\sqrt{T}}\exp \left(\frac{D{G}^2{\ln }^2(\delta /G)}{2T}+\Vert V_{-}{\Vert }_{\infty }T\right),$$

$${\partial }_t\tilde{u}(t)+(-{\Delta }_{\tilde{\Omega }}^{\bullet }+\tilde{V})\tilde{u}(t)={\chi }_{\tilde{\omega }}\tilde{v}(t)\text{ for }0<t<T,\tilde{u}(0)={\tilde{u}}_0,$$

$$C_T\le {\left(\frac{\delta }{\tilde{G}}\right)}^{-D\cdot (1+{\tilde{G}}^{4/3}\Vert V{\Vert }_{\infty }^{2/3})}\frac{D}{\sqrt{T}}\exp \left(\frac{D{\tilde{G}}^2{\ln }^2(\delta /\tilde{G})}{2T}+\Vert V_{-}{\Vert }_{\infty }T\right),$$

$$C_T\le {\left(\frac{\delta }{2G}\right)}^{-R\cdot (1+(2G{)}^{4/3}\Vert V{\Vert }_{\infty }^{2/3})}\frac{R}{\sqrt{T}}\exp \left(\frac{R(2G{)}^2{\ln }^2(\delta /(2G))}{2T}+\Vert V_{-}{\Vert }_{\infty }T\right),$$

$${\partial }_{t}u(t)+Hu(t)=Bv(t)\text{ for }0<t<T,u(0)=u_0,$$

Full text

1. Introduction and main result

Let $\Omega\subset\mathbb{R}^{d}$ be open, $d\in\mathbb{N}$, and let $A\colon\Omega\rightarrow\mathbb{R}^{d\times d}$ be measurable with $A(x)$ a symmetric matrix for almost every $x\in\Omega$. Suppose, in addition, that there exist $\theta_{1},\theta_{2}>0$ such that

[TABLE]

and let $V\in L^{\infty}(\Omega)$ be real-valued. We denote by $H_{\Omega}^{D}=H_{\Omega}^{D}(A,V)$ and $H_{\Omega}^{N}=H_{\Omega}^{N}(A,V)$ the Dirichlet and Neumann realizations of the differential expression

[TABLE]

as self-adjoint lower semibounded operators on $L^{2}(\Omega)$ defined via their quadratic forms with form domain $H_{0}^{1}(\Omega)$ and $H^{1}(\Omega)$, respectively. For details of this construction we refer the reader to the discussion in Section 4 below.

Let $T>0$, $\omega\subset\Omega$ be a measurable subset, and $\bullet\in\{D,N\}$. We consider the heat-like system

[TABLE]

with $u_{0}\in L^{2}(\Omega)$ and $v\in L^{2}((0,T),L^{2}(\Omega))$. Here, $\chi_{\omega}$ denotes the characteristic function of $\omega$, and we call $\omega$ a control set for the system (1.2).

System (1.2) is said to be null-controllable in time $T>0$ if for every initial data $u_{0}\in L^{2}(\Omega)$ there exists a control function $v\in L^{2}((0,T),L^{2}(\Omega))$ such that the mild solution of (1.2) satisfies $u(T)=0$, that is,

[TABLE]

In this case, the quantity

[TABLE]

is called control cost.

For bounded domains $\Omega$, there is already a rich literature on null-controllability results for system (1.2) with $H_{\Omega}^{\bullet}$ being the Laplace or a Schrödinger operator with an open or measurable control set $\omega$, see, e.g., [16, 28, 2, 20, 10]. On the other hand, null-controllability of these systems on unbounded domains has been an issue of growing interest only recently, see for example [9, 10, 20] and the references therein. As far as bounds on the control cost are concerned, one classically asks for the type of time dependency, which has been largely explored in [12, 22, 28, 7, 19, 11, 17, 20], see also the references therein. However, the dependency of the control cost on geometric parameters of $\Omega$ and $\omega$ has only recently been addressed in [9, 8, 10, 20], see also [18] for previous results. Especially the work [20] focused on these dependencies, which have been exploited in several asymptotic regimes there in order to discuss homogenization.

The aim of the present paper is to prove null-controllability results for the above situation on some new unbounded domains, such as sectors of certain angles, half-spaces, and orthants, with bounds on the associated control cost that are fully or partially consistent with the known ones on cubes and the full space. As a byproduct, we also obtain explicit control cost bounds in terms of only geometric parameters of the control set for the corresponding system on some bounded domains with flat boundary that were not accessible before. We believe that the lack of the control cost’s dependency on the domain may be an effect of the flatness of its boundary. However, it is not our current scope to explore this matter in depth, but merely to add examples to support this conjecture. All these results are discussed in detail in Section 2 below as applications of a reflection principle which we now describe.

The key idea in our considerations is to relate system (1.2) for certain domains $\Omega$ to a corresponding symmetrized system on a larger domain, where null-controllability results are available. In order to make this precise, let $M\colon\mathbb{R}^{d}\rightarrow\mathbb{R}^{d}$ be the reflection with respect to the first coordinate, that is, $M(x_{1},\ldots,x_{d})=(-x_{1},x_{2},\ldots,x_{d})$. We suppose that $\Omega$ is contained in the half-space $(0,+\infty)\times\mathbb{R}^{d-1}$ such that $\Omega=\tilde{\Omega}\cap((0,+\infty)\times\mathbb{R}^{d-1})$ for some open set $\tilde{\Omega}\subset\mathbb{R}^{d}$ with $\Gamma:=\tilde{\Omega}\cap(\{0\}\times\mathbb{R}^{d-1})\neq\emptyset$ and which is symmetric with respect to the reflection $M$, that is, $M(\tilde{\Omega})=\tilde{\Omega}$. In particular, we have $\tilde{\Omega}=\Omega\cup\Gamma\cup M(\Omega)$, cf. Figure 1.1.

Let $\tilde{A}\colon\tilde{\Omega}\rightarrow\mathbb{R}^{d\times d}$ and $\tilde{V}\colon\tilde{\Omega}\rightarrow\mathbb{R}$ be measurable with

[TABLE]

and

[TABLE]

Observe that by construction $\tilde{A}(x)$ is for almost every $x\in\tilde{\Omega}$ a symmetric matrix satisfying (1.1) with the same constants $\theta_{1},\theta_{2}$.

Let the self-adjoint operators $H_{\tilde{\Omega}}^{\bullet}=H_{\tilde{\Omega}}^{\bullet}(\tilde{A},\tilde{V})$ on $L^{2}(\tilde{\Omega})$ associated with the differential expression

[TABLE]

be defined analogously to $H_{\Omega}^{\bullet}(A,V)$ above. Set $\tilde{\omega}:=\omega\cup M(\omega)\subset\tilde{\Omega}$, and consider the corresponding symmetrized system

[TABLE]

with $\tilde{u}_{0}\in L^{2}(\tilde{\Omega})$, $\tilde{v}\in L^{2}((0,T),L^{2}(\tilde{\Omega}))$, and control set $\tilde{\omega}$. The notions of null-controllability and control cost for this system carry over verbatim.

The main result of the present paper now reads as the following theorem, which can actually be formulated also for much more general situations than the ones we discuss here, cf. Section 3 below.

Theorem 1.1**.**

If the system (1.5) is null-controllable in time $T>0$ with control cost $\tilde{C}_{T}$, then also system (1.2) is null-controllable in time $T>0$ with control cost $C_{T}\leq\tilde{C}_{T}$.

The restriction to the reflection symmetry with respect to the hyperplane $\{0\}\times\mathbb{R}^{d-1}$ in Theorem 1.1 is not essential. Indeed, by rotating the whole system, we can deal with reflection symmetries with respect to any hyperplane in $\mathbb{R}^{d}$. This way, the above theorem allows us to infer null-controllability on the respective parts of a domain with a reflection symmetry if the system on the whole domain is null-controllable.

The rest of the paper is organized as follows. In Section 2 we discuss the null-controllability results and control cost bounds announced above for the case where $A(x)$ is the identity matrix as consequences of Theorem 1.1. Here we will make explicit use of recent results from [10, 20] for cubes and the whole space. The assumption on the essential boundedness of the potential $V$ is tailored towards these applications but it is not essential for the general argument. In Section 3 we prove an abstract result which is not only the core of the proof of Theorem 1.1 but also promises to have a much broader range of application and is therefore of its own interest, see, e.g., Remarks 2.10, 4.2, and 4.3 below. Section 4 deals with the proof of Theorem 1.1 based on the abstract result of the preceding section. Finally, Appendix A provides an integration by parts formula used in Section 4.

2. Applications

We here discuss the applications of Theorem 1.1 mentioned in the previous section for the particular case where each $A(x)$ is the identity matrix, that is, $H_{\Omega}^{\bullet}=-\Delta_{\Omega}^{\bullet}+V$. These applications draw upon null-controllability results from the works [10, 20] on cubes and the full space. As considered there, we also take control sets of the form $\omega=\Omega\cap S$, where $S\subset\mathbb{R}^{d}$ is some measurable set with certain geometric properties, namely a thick set or an equidistributed set. We treat the two types of sets in the subsections below separately.

The general strategy for our applications is to build from the given set $S$ a symmetric set $\tilde{S}$ with respect to the reflection such that $\Omega\cap\tilde{S}=\Omega\cap S$. After having verified corresponding geometric properties of $\tilde{S}$, this new set yields a suitable control set $\tilde{\omega}=\tilde{\Omega}\cap\tilde{S}$ for the symmetrized system on $\tilde{\Omega}$, where a null-controllability result is available. Applying Theorem 1.1 then gives the desired null-controllability result for the original system on $\Omega$. In summary, our main task here is to establish the geometric properties of the symmetrized set $\tilde{S}$ from those of the given set $S$.

Recall that by rotation of the whole system Theorem 1.1 can be applied for a reflection with respect to any hyperplane in $\mathbb{R}^{d}$. Since the Laplacian is rotation invariant, only the potential $V$ then gets rotated, but remains essentially bounded. In this context, we remind the reader that $M\colon\mathbb{R}^{d}\to\mathbb{R}^{d}$ denotes the reflection with respect to the hyperplane $\{0\}\times\mathbb{R}^{d-1}$. Furthermore, for the rest of this section, we use the following notation:

• •

$H_{\theta}:=\{(x_{1},\ldots,x_{d})\in\mathbb{R}^{d}\,\colon\,x_{2}<(\tan\theta)x_{1}\}$ for $0\leq\theta<\pi/2$;

• •

$M_{\theta}\colon\mathbb{R}^{d}\rightarrow\mathbb{R}^{d}$ denotes the reflection with respect to $\partial H_{\theta}$;

• •

$\Lambda_{L}^{d}:=(0,L)^{d}$ for $L>0$.

2.1. Null-controllability from thick sets

Let $V=0$, and consider the system

[TABLE]

with $u_{0}\in L^{2}(\Omega)$ and $v\in L^{2}((0,T),L^{2}(\Omega))$. Moreover, assume that the control set $\omega$ is of the form $\omega=\Omega\cap S$ with some $S\subset\mathbb{R}^{d}$ that is $(\gamma,a)$-thick in the sense of the following definition.

Definition 2.1**.**

A measurable set $S\subset\mathbb{R}^{d}$ is called thick if there exist $\gamma\in(0,1]$ and $a=(a_{1},\ldots,a_{d})\in(0,+\infty)^{d}$ such that

[TABLE]

where $\lvert{\,\cdot\,}\rvert$ denotes the Lebesgue measure. In this case, $S$ is also referred to as $(\gamma,a)$-thick to emphasise the parameters.

The above definition has played a crucial role in a recent development on the null-controllability of the heat equation on $\mathbb{R}^{d}$. In [10], see also [30], it is shown that thickness of $S$ is a necessary and sufficient condition for the heat equation (2.1) on $\Omega=\mathbb{R}^{d}$ to be null-controllable. Moreover, in [10] the authors give an explicit estimate of the control cost in dependence of the thickness parameters of the set $S$ and the time $T$. In the same paper the authors also consider the heat equation (2.1) on the cube $\Omega=\Lambda_{L}^{d}$ with control set $\omega=\Omega\cap S$, where $S$ is $(\gamma,a)$-thick with $a_{j}\leq L$ for all $j=1,\ldots,d$. In this case, they show that null-controllability holds in any time $T>0$ with a bound on the control cost of the same form as for the full space case and independent of the scale $L$.

In [20] these bounds have been strengthened in both time and geometric parameters dependency and the authors have shown that they are close to optimality in certain asymptotic regimes, see [20, Section 4]. For both cases $\Omega\in\{\mathbb{R}^{d},\Lambda_{L}^{d}\}$ they can be written as

[TABLE]

where $K>0$ is a universal constant and $\|{a}\|_{1}=a_{1}+\dots+a_{d}$, see [20, Theorem 3.7].

We show below that system (2.1) on half-spaces, positive orthants, and sectors of angle $\pi/2^{n}$, $n\geq 2$, is null-controllable with a control cost bound fully or partially consistent with (2.2). Moreover, we obtain null-controllability and explicit control cost bounds for the same system on isosceles right-angled triangles and corresponding prisms. We once again emphasise that the results for the listed unbounded domains are new, even though the heat equation on half-spaces and positive orthants with Dirichlet boundary conditions could also be treated with the techniques from [27] (see also Remark 2.6 below), while the novelty of the result on triangles and corresponding prisms is the explicitness of the control cost bound on the model parameters.

In order to enter the setting of the main theorem, let $\tilde{\Omega}\in\{\mathbb{R}^{d},\Lambda_{L}^{d}\}$ and $\tilde{\omega}=\tilde{\Omega}\cap\tilde{S}$ with some $(\tilde{\gamma},\tilde{a})$-thick set $\tilde{S}\subset\mathbb{R}^{d}$, and consider the system

[TABLE]

with $\tilde{u}_{0}\in L^{2}(\tilde{\Omega})$ and $\tilde{v}\in L^{2}((0,T),L^{2}(\tilde{\Omega}))$.

We need the following easy geometric lemma, parts of which are already contained in [10].

Lemma 2.2**.**

Let $S\subset\mathbb{R}^{d}$ be $(\gamma,a)$-thick.

• (a)

Let $S^{\prime}=S\cap((0,+\infty)\times\mathbb{R}^{d-1})$. Then, the set $\tilde{S}=S^{\prime}\cup M(S^{\prime})$ is $(\tilde{\gamma},\tilde{a})$-thick with $\tilde{\gamma}=\gamma/2$ and $\tilde{a}=(2a_{1},a_{2},\ldots,a_{d})$.

• (b)

The set $\tilde{S}=\{(x_{1},\dots,x_{d})\,\colon\,(\lvert{x_{1}}\rvert,\dots,\lvert{x_{d}}\rvert)\in S\}$ is $(\gamma/2^{d},2a)$-thick, where $2a=(2a_{1},\ldots,2a_{d})$.

• (c)

Let $S^{\prime}=S\cap H_{\theta}$. Then, $\tilde{S}=S^{\prime}\cup M_{\theta}(S^{\prime})$ is a $(\tilde{\gamma},\tilde{a})$-thick set with parameters $\tilde{\gamma}=\frac{\gamma a_{1}a_{2}}{4(a_{1}^{2}+a_{2}^{2})}$ and $\tilde{a}=(2\sqrt{a_{1}^{2}+a_{2}^{2}},2\sqrt{a_{1}^{2}+a_{2}^{2}},a_{3},\ldots,a_{d})$.

Proof.

We abbreviate $Q_{a}:=[0,a_{1}]\times\dots\times[0,a_{d}]$, $a=(a_{1},\dots,a_{d})\in(0,+\infty)^{d}$.

For part (a) we need to show that

[TABLE]

Let $x\in\mathbb{R}^{d}$. There is $y\in\mathbb{R}^{d}$ such that $y+Q_{a}\subset(x+Q_{\tilde{a}})\cap([0,+\infty)\times\mathbb{R}^{d-1})$ or $M(y+Q_{a})\subset(x+Q_{\tilde{a}})\cap((-\infty,0])\times\mathbb{R}^{d-1})$, cf. Figure 2.1 (a). In the first case, we have

[TABLE]

In the second case, we have

[TABLE]

which completes the proof of part (a).

For part (b), we observe that $\tilde{S}$ can be obtained from $S\cap[0,+\infty)^{d}$ by successive reflection with respect to all coordinate axes. In this regard, part (b) follows by analogous arguments as in part (a), cf. Figure 2.1 (b); for more details, see also the proof of [10, Theorem 4].

Finally, we prove part (c). Let $x\in\mathbb{R}^{d}$. Then, there exists $y\in\mathbb{R}^{d}$ such that $y+Q_{a}\subset(x+Q_{\tilde{a}})\cap\overline{H_{\theta}}$ or $M_{\theta}(y+Q_{a})\subset(x+Q_{\tilde{a}})\cap(\mathbb{R}^{d}\setminus H_{\theta})$. This follows from the fact that $x+Q_{\tilde{a}}$ contains the cylinder x^{\prime}+\big{(}B_{2}(0,\sqrt{a_{1}^{2}+a_{2}^{2}})\bigtimes_{j=3}^{d}[0,a_{j}]\big{)} for some $x^{\prime}\in\mathbb{R}^{d}$, where $B_{2}(0,r)$ stands for the Euclidean 2-dimensional ball centred at zero with radius $r$, and such a cylinder contains a parellelepiped of type $y+Q_{a}$ as well as its reflection with respect to $M_{\theta}$, cf. Figure 2.1 (c). The claim now follows by analogous calculations as in part (a), taking into account that $\gamma\lvert{Q_{a}}\rvert=\frac{\gamma a_{1}a_{2}}{4(a_{1}^{2}+a_{2}^{2})}\lvert{Q_{\tilde{a}}}\rvert=\tilde{\gamma}\lvert{Q_{\tilde{a}}}\rvert$. ∎

Remark 2.3*.*

The choice of the parameters $\tilde{\gamma}$ and $\tilde{a}$ in part (c) of Lemma 2.2 takes care of every possible orientation of reflected parallelepipeds, which are in general not parallel to coordinate axes (with the exception of $\theta=\pi/4$). In this regard, the angle $\theta=\pi/8$ constitutes a kind of worst case, whereas especially the angle $\theta=\pi/4$ could have been handled in a slightly more efficient way. We refrained from doing so for the sake of simplicity.

Proposition 2.4** (The heat equation on half-spaces).**

The system (2.1) on $\Omega=(0,+\infty)\times\mathbb{R}^{d-1}$ is null-controllable in any time $T>0$ with control cost

[TABLE]

where $K>0$ is the universal constant from (2.2).

Proof.

Choose the $(\gamma/2,(2a_{1},a_{2}\ldots,a_{d}))$-thick set $\tilde{S}$ as in Lemma 2.2 (a). Then, the heat equation (2.3) on $\tilde{\Omega}=\mathbb{R}^{d}$ is null-controllable in any time $T>0$ from the control set $\tilde{\omega}=\tilde{S}$. The associated bound on the control cost $\tilde{C}_{T}$ from (2.2) now reads as in (2.4).

Since $\tilde{\Omega}=\Omega\cup M(\Omega)\cup\Gamma$ with $\Gamma=\{0\}\times\mathbb{R}^{d-1}$ and $\tilde{\omega}=\omega\cup M(\omega)$, the claim follows by Theorem 1.1. ∎

Proposition 2.5** (Heat equation on positive orthants).**

The system (2.1) on $\Omega=(0,+\infty)^{d}$ is null-controllable in any time $T>0$ with control cost

[TABLE]

where $K>0$ is the universal constant from (2.2).

Proof.

Choose the $(\gamma/2^{d},2a)$-thick set $\tilde{S}$ as in Lemma 2.2 (b). Then, the heat equation (2.3) on $\tilde{\Omega}=\mathbb{R}^{d}$ is null-controllable in any time $T>0$ from the control set $\tilde{\omega}=\tilde{S}$. The associated bound on the control cost $\tilde{C}_{T}$ now reads as in (2.5).

Since the set $\tilde{S}$ is symmetric with respect to every coordinate axis, the claim now follows by successively applying Theorem 1.1 for the reflections with respect to the coordinate axes. This step by step leads to null-controllability of the system (2.1) on $\Omega=(0,+\infty)^{j}\times\mathbb{R}^{d-j}$, $j=1,\dots,d$, with control set $\omega=((0,+\infty)^{j}\times\mathbb{R}^{d-j})\cap\tilde{S}$. Each control cost does not exceed the one from the previous step and, thus, the corresponding bound reads as in (2.5). Taking into account that $(0,+\infty)^{d}\cap\tilde{S}=(0,+\infty)^{d}\cap S$, the final step $j=d$ then proves the claim. ∎

Remark 2.6*.*

Control cost estimates on the half-space and the positive orthant can alternatively be obtained by means of the exhaustion approach studied in [27]. In this case, the corresponding bounds will not incorporate any change in the thickness parameters as the above bounds do. However, the considerations from [27] currently allow to apply this approach only in the case of Dirichlet boundary conditions.

Proposition 2.7** (Heat equation on a sector of angle $\pi/2^{n}$, $n\geq 2$).**

The system (2.1) on the sector $\Omega=\{(x,y)\in(0,+\infty)^{2}\,\colon\,y<(\tan(\pi/2^{n}))x\}$ of angle $\pi/2^{n}$, $n\geq 2$, is null-controllable in any time $T>0$ with control cost

[TABLE]

where

[TABLE]

and $K>0$ is the universal constant from (2.2).

Proof.

We prove the claim by induction on $n$. For $n=2$, we choose the $(\tilde{\gamma},\tilde{a})$-thick set $\tilde{S}$ as in Lemma 2.2 (c) with $d=2$ and $\theta=\pi/4$. Then, the heat equation (2.3) on the orthant $\tilde{\Omega}=(0,+\infty)^{2}$ is null-controllable in any time $T>0$, and the bound on the associated control cost from Proposition 2.5 reads as in (2.6) for $n=2$. The claim for $n=2$ then follows by applying Theorem 1.1 with respect to the reflection $M_{\pi/4}$.

Now, suppose that the claim holds for some $n\geq 2$. Choose the $(\tilde{\gamma},\tilde{a})$-thick set $\tilde{S}$ as in Lemma 2.2 (c) with $d=2$ and $\theta=\pi/2^{n+1}$. Since $\tilde{a}=(\tilde{a}_{1},\tilde{a}_{2})$ is a multiple of $(1,1)$, we observe that

[TABLE]

Thus, by the induction hypothesis, the system (2.3) with $\tilde{\Omega}$ being the sector of angle $\pi/2^{n}$ is null-controllable in any time $T$, and the corresponding bound on the control cost reads as in (2.6), but with $\tilde{\gamma}$ and $\tilde{a}$ replaced by $\tilde{\gamma}/8$ and $2\sqrt{2}\tilde{a}$, respectively. The claim for $\Omega$ the sector of angle $\pi/2^{n+1}$ therefore follows by applying Theorem 1.1 with respect to the reflection $M_{\pi/2^{n+1}}$. This concludes the proof. ∎

For the last proposition of this subsection, we suppose that the thickness parameter $a=(a_{1},\dots,a_{d})$ additionally satisfies $2\sqrt{a_{1}^{2}+a_{2}^{2}}\leq L$ and $a_{j}\leq L$ for $j\in\{3,\ldots,d\}$.

Proposition 2.8** (Heat equation on triangles and triangular prisms).**

Let

[TABLE]

Then, the system (2.1) on $\Omega\in\{{\mathcal{T}}_{L},{\mathcal{P}}_{L}\}$ is null-controllable in any time $T>0$ with control cost

[TABLE]

where $\tilde{\gamma}$ and $\tilde{a}$ are as in part (c) of Lemma 2.2 and $K>0$ is the universal constant from (2.2). Here, $d=2$ if $\Omega={\mathcal{T}}_{L}$ and $d=3$ if $\Omega={\mathcal{P}}_{L}$.

Proof.

This is proved analogously to the case $n=2$ in Proposition 2.7 with $\tilde{\Omega}=\Lambda_{L}^{d}$ and the corresponding control cost bound from (2.2). ∎

Remark 2.9*.*

In the recent work [8], it has been shown that the system (2.1) on the strip $\Omega=\Lambda_{L}^{d-1}\times\mathbb{R}$ is null-controllable in any time $T>0$ if and only if $S$ is a thick set (which can be arbitrarily changed outside the strip), and an explicit control cost bound has been provided.

With similar arguments as in the propositions above, one can infer null-controllability results for system (2.1) on the product spaces $\Lambda_{L}^{d-1}\times(0,+\infty)$, ${\mathcal{T}}_{L}\times\mathbb{R}$, and ${\mathcal{T}}_{L}\times(0,+\infty)$ by means of Theorem 1.1.

Remark 2.10*.*

(a) Null-controllability results and corresponding bounds on the control cost for fractional heat equations, that is, system (2.3) with $(-\Delta_{\tilde{\Omega}}^{\bullet})^{\theta}$, $\theta>\frac{1}{2}$, instead of $-\Delta_{\tilde{\Omega}}^{\bullet}$, have also been investigated in [20, Theorem 4.10] on $\tilde{\Omega}\in\{\mathbb{R}^{d},\Lambda_{L}^{d}\}$ with a control cost bound of the form (2.2) but with $\|{a}\|_{1}^{2}\ln^{2}(K^{d}/\gamma)/T$ replaced by $(\|{a}\|_{1}\ln(K^{d}/\gamma))^{2\theta/(2\theta-1)}/T^{1/(2\theta-1)}$.

Although fractional Laplacians are formally not of divergence form as the operators considered in Theorem 1.1, the more general considerations in Section 3 below nevertheless allow us to obtain corresponding results on the domains considered in Propositions 2.4–2.8, cf. Remark 4.2 (a) below.

(b) Recent developments on null-controllability for parabolic equations in Banach spaces [4, 13] together with the abstract considerations in Section 3 also open the way towards similar results for certain strongly elliptic differential operators with constant coefficients on $L^{p}$, $1<p<\infty$; a bound on the associated control cost on $\mathbb{R}^{d}$ of a form close to (2.2) is given in [4, Theorem 3.6], see also [13, Corollary 4.6]. We briefly demonstrate how to infer a corresponding result for the $L^{p}$-Laplacian on the half-space $(0,+\infty)\times\mathbb{R}^{d-1}$ in Remark 4.2 (b) below.

2.2. Null-controllability from equidistributed sets

Let $V\in L^{\infty}(\Omega)$ be real-valued, and consider the heat-like system

[TABLE]

with $u_{0}\in L^{2}(\Omega)$ and $v\in L^{2}((0,T),L^{2}(\Omega))$. Moreover, let the control set $\omega$ be of the form $\omega=\Omega\cap S$ with some $S\subset\mathbb{R}^{d}$ that is $(G,\delta)$-equidistributed in the sense of the following definition.

Definition 2.11**.**

Let $G>0$ and $\delta\in(0,G/2)$. A measurable set $S\subset\mathbb{R}^{d}$ is called $(G,\delta)$-equidistributed if

[TABLE]

for some sequence $(z_{j})_{j\in\mathbb{Z}^{d}}\subset\mathbb{R}^{d}$ with $B(z_{j},\delta)\subset\left(Gj+\Lambda_{G}^{d}\right)$ for all $j\in\mathbb{Z}^{d}$, where $B(z_{j},\delta)$ denotes the open ball in $\mathbb{R}^{d}$ of radius $\delta$ around $z_{j}$.

It is shown in [20] that system (2.8) is null-controllable in any time $T>0$ on domains of the form $\Omega=\bigtimes_{j=1}^{d}(a_{j},b_{j})$ with $a_{j},b_{j}\in\mathbb{R}\cup\{\pm\infty\}$, $a_{j}<b_{j}$, such that $\Lambda_{G}\subset\Omega$. A corresponding bound on the control cost is given by

[TABLE]

where $V_{-}=\max\{-V,0\}$ denotes the negative part of $V$ and $D=D(d)$ depends only on the dimension, cf. [20, Theorem 4.11] and [27, Proposition 2.4 and Corollary 2.5].

Below we extend the above result to sectors of angle $\pi/2^{n}$, $n\geq 2$, isosceles right-angled triangles, and corresponding prisms. A treatment of some related product spaces would also be possible, but we will not do this here for simplicity. In the framework of the main theorem, let $\tilde{\Omega}\in\{\Lambda_{L}^{d},(0,+\infty)^{d}\}$, $L\geq 2G$, and $\tilde{\omega}=\tilde{\Omega}\cap\tilde{S}$ with some $(2G,\delta)$-equidistributed set for $\tilde{S}$. As in the previous subsection, we use the above mentioned results for the system

[TABLE]

with $\tilde{u}_{0}\in L^{2}(\tilde{\Omega})$ and $\tilde{v}\in L^{2}((0,T),L^{2}(\tilde{\Omega}))$ to infer null-controllability and related control cost bounds on the desired domains.

Analogously to Lemma 2.2 in the case of thick sets, we first need the following geometric consideration.

Lemma 2.12**.**

Let $S\subset\mathbb{R}^{d}$ be $(G,\delta)$-equidistributed, and let $S^{\prime}=S\cap H_{\theta}$. Then, the set $\tilde{S}=S^{\prime}\cup M_{\theta}(S^{\prime})$ is $(4G,\delta)$-equidistributed.

Proof.

We need to show that for every cube $\Lambda_{j}:=4Gj+\Lambda_{4G}^{d}$, $j\in\mathbb{Z}^{d}$, there exists $j^{\prime}\in\mathbb{Z}^{d}$ such that $Gj^{\prime}+\Lambda_{G}^{d}\subset\Lambda_{j}\cap\overline{H_{\theta}}$ or $M_{\theta}(Gj^{\prime}+\Lambda_{G}^{d})\subset\Lambda_{j}\cap(\mathbb{R}^{d}\setminus H_{\theta})$.

To this end, we observe that for every $j\in\mathbb{Z}^{d}$ there exists $k\in\mathbb{Z}^{d}$ such that $Gk+\Lambda_{G}^{d}\subset\Lambda_{j}\cap\overline{H_{\theta}}$ or $3Gk+\Lambda_{3G}^{d}\subset\Lambda_{j}\cap(\mathbb{R}^{d}\setminus H_{\theta})$, cf. Figure 2.2 (a). In the first case, we are done with $j^{\prime}=k$. In the second case, the cube $3Gk+\Lambda_{3G}^{d}$ contains a ball of radius $\sqrt{2}G$, which, in turn, contains the reflection of a cube $Q\subset H_{\theta}$ with sides of length $2G$ parallel to coordinate axes, cf. Figure 2.2 (b). This cube $Q$ must contain at least one cube of the form $Gj^{\prime}+\Lambda_{G}^{d}$ with $j^{\prime}\in\mathbb{Z}^{d}$, for which $M_{\theta}(Gj^{\prime}+\Lambda_{G}^{d})\subset 3Gk+\Lambda_{3G}^{d}\subset\Lambda_{j}\cap(\mathbb{R}^{d}\setminus H_{\theta})$. This concludes the proof. ∎

Remark 2.13*.*

If $\theta=\pi/4$, it is possible to slightly strengthen the above lemma. Indeed, in this case the resulting set $\tilde{S}$ is $(2G,\delta)$-equidistributed since the reflection $M_{\pi/4}$ maps cells of the lattice $(G\mathbb{Z})^{d}$ in $\overline{H_{\pi/4}}$ to cells of the lattice in $\mathbb{R}^{d}\setminus H_{\pi/4}$ and vice versa, and each cube $2Gj+\Lambda_{2G}^{d}$, $j\in\mathbb{Z}^{d}$, contains at least one cell of the lattice $(G\mathbb{Z})^{d}$ that belongs either to $\overline{H_{\pi/4}}$ or to $\mathbb{R}^{d}\setminus H_{\pi/4}$. However, for the sake of simplicity we opted for a unified statement valid for every reflection $M_{\theta}$.

Proposition 2.14** (Heat-like equation on a sector of angle $\pi/2^{n}$, $n\geq 2$).**

The system (2.8) on the sector $\Omega=\{(x,y)\in(0,+\infty)^{2}\,\colon\,y<(\tan(\pi/2^{n}))x\}$ of angle $\pi/2^{n}$, $n\geq 2$, is null-controllable in any time $T>0$ with control cost

[TABLE]

where $\tilde{G}=4^{n-1}G$ and $D=D(2)$ is the constant from (2.9) for $d=2$.

Proof.

This is proved by induction on $n$ analogously to Proposition 2.7, the only difference being the use of Lemma 2.12 instead of Lemma 2.2 (c). ∎

Proposition 2.15** (Heat-like equation on triangles and triangular prisms).**

Let $L\geq 2G$, and let ${\mathcal{T}}_{L}$ and ${\mathcal{P}}_{L}$ as in Proposition 2.8. Then, the system (2.8) on $\Omega\in\{{\mathcal{T}}_{L},{\mathcal{P}}_{L}\}$ is null-controllable in any time $T>0$ with control cost

[TABLE]

where $R=\max\{D(2),D(3)\}$ with $D(d)$ from (2.9).

Proof.

Taking into account Remark 2.13, this is proved analogously to the case $n=2$ in Proposition 2.14 with $\tilde{\Omega}=\Lambda_{L}^{d}$. ∎

3. Abstract Result

In this section we prove a general version of Theorem 1.1 for more abstract control systems, which in principle allows one to apply this result, for instance, also to other types of differential operators than discussed here so far, such as fractional Laplacians, second order elliptic operators, and magnetic Schrödinger operators.

Let ${\mathcal{H}}$ and ${\mathcal{U}}$ be Banach spaces, $-H$ an infinitesimal generator of a $C_{0}$-semigroup $(S(t))_{t\geq 0}$ on ${\mathcal{H}}$, $B\colon{\mathcal{U}}\to{\mathcal{H}}$ a bounded linear operator, and $T>0$. The abstract Cauchy problem

[TABLE]

with $u_{0}\in{\mathcal{H}}$ and $v\in L^{1}((0,T),{\mathcal{U}})$ is said to be null-controllable in time $T>0$ with respect to $L^{p}((0,T),{\mathcal{U}})$, $p\in[1,\infty]$, if for every initial datum $u_{0}\in{\mathcal{H}}$ there is a function $v\in L^{p}((0,T),{\mathcal{U}})$ with

[TABLE]

that is, if the mild solution to (3.1) vanishes at time $T$; see, e.g., [21] for the notion of a mild solution to abstract Cauchy problems. The associated control cost in time $T>0$ is then defined as

[TABLE]

The main result of this section is as follows.

Theorem 3.1**.**

Let ${\mathcal{H}},\widetilde{{\mathcal{H}}},{\mathcal{U}},\widetilde{{\mathcal{U}}}$ be Banach spaces, $-H$ and $-\widetilde{H}$ infinitesimal generators of $C_{0}$-semigroups on ${\mathcal{H}}$ and $\widetilde{{\mathcal{H}}}$, respectively, and let $B\colon{\mathcal{U}}\to{\mathcal{H}}$ and $\widetilde{B}\colon\widetilde{{\mathcal{U}}}\to\widetilde{{\mathcal{H}}}$ be bounded linear operators.

Suppose that there is a bounded linear operator $Y\colon\widetilde{{\mathcal{H}}}\to{\mathcal{H}}$ with a bounded right inverse $\hat{Y}\colon{\mathcal{H}}\to\widetilde{{\mathcal{H}}}$ such that

[TABLE]

that is, $Yg\in\operatorname{{\mathcal{D}}}(H)$ and $Y\widetilde{H}g=HYg$ for all $g\in\operatorname{{\mathcal{D}}}(\widetilde{H})$. Furthermore, suppose that for some bounded linear operator $Z\colon\widetilde{{\mathcal{U}}}\to{\mathcal{U}}$ one has

[TABLE]

If the system

[TABLE]

with $\widetilde{u}_{0}\in\widetilde{{\mathcal{H}}}$ and $\widetilde{v}\in L^{1}((0,T),\widetilde{{\mathcal{U}}})$ is null-controllable in time $T>0$ with respect to $L^{p}((0,T),\widetilde{{\mathcal{U}}})$ with control cost $\widetilde{C}_{T,p}>0$, then also the system (3.1) is null-controllable in time $T>0$ with respect to $L^{p}((0,T),{\mathcal{U}})$ with control cost $C_{T,p}>0$ satisfying

[TABLE]

Proof.

Let $(S(t))_{t\geq 0}$ and $(\widetilde{S}(t))_{t\geq 0}$ be the $C_{0}$-semigroups associated to $-H$ and $-\widetilde{H}$, respectively. We first show that

[TABLE]

To this end, let $\widetilde{u}_{0}\in\operatorname{{\mathcal{D}}}(\widetilde{{\mathcal{H}}})$, and define the function $w\colon[0,\infty)\to{\mathcal{H}}$ by $w(t):=Y\widetilde{S}(t)\widetilde{u}_{0}$. Since $\widetilde{S}(t)\widetilde{u}_{0}\in\operatorname{{\mathcal{D}}}(\widetilde{H})$ holds for all $t\geq 0$, we then have $w(t)\in\operatorname{{\mathcal{D}}}(H)$ and, in particular, $w(0)=Y\widetilde{u}_{0}\in\operatorname{{\mathcal{D}}}(H)$. Moreover,

[TABLE]

for $t>0$. By uniqueness of the solution, see, e.g., [21, Theorem 4.1.3], we conclude that $w(t)=S(t)w(0)$, that is, $Y\widetilde{S}(t)\widetilde{u}_{0}=S(t)Y\widetilde{u}_{0}$ for all $t\geq 0$. Since $\operatorname{{\mathcal{D}}}(\widetilde{H})$ is dense in $\widetilde{{\mathcal{H}}}$ and $S(t)$ and $\widetilde{S}(t)$ are bounded operators, this proves (3.6).

Let now $u_{0}\in{\mathcal{H}}$, and set $\widetilde{u}_{0}:=\hat{Y}u_{0}\in\widetilde{{\mathcal{H}}}$. Then, $Y\widetilde{u}_{0}=u_{0}$, and by hypothesis there is $\widetilde{v}\in L^{p}((0,T),\widetilde{{\mathcal{U}}})$ such that

[TABLE]

Define $v\in L^{p}((0,T),{\mathcal{U}})$ by $v(s):=Z\widetilde{v}(s)$. By (3.4) and (3.6), we then conclude that

[TABLE]

so that system (3.1) is indeed null-controllable in time $T>0$ with respect to $L^{p}((0,T),{\mathcal{U}})$.

Finally, using $\|{v(s)}\|_{{\mathcal{U}}}\leq\|{Z}\|_{\widetilde{{\mathcal{U}}}\to{\mathcal{U}}}\|{\widetilde{v}(s)}\|_{\widetilde{{\mathcal{U}}}}$ for $s\in(0,T)$, we observe that $\|{v}\|_{L^{p}((0,T),{\mathcal{U}})}\leq\|{Z}\|_{\widetilde{{\mathcal{U}}}\to{\mathcal{U}}}\|{\widetilde{v}}\|_{L^{p}((0,T),\widetilde{{\mathcal{U}}})}$. Since also $\|{\widetilde{u}_{0}}\|_{\widetilde{{\mathcal{H}}}}\leq\|{\hat{Y}}\|_{{\mathcal{U}}\to\widetilde{{\mathcal{U}}}}\|{u_{0}}\|_{{\mathcal{H}}}$, the claimed estimate for the control cost $C_{T,p}$ is immediate. This completes the proof. ∎

Remark 3.2*.*

Since the extension relation (3.3) is only used to establish the semigroup relation (3.6), the conclusion of Theorem 3.1 still holds if (3.6) is directly taken as an assumption.

We refrained from doing so because the extension relation between the operators gives a nice understanding of the interplay between the systems (3.1) and (3.5).

In the particular case where $H$ and $\widetilde{H}$ are lower semibounded self-adjoint operators on Hilbert spaces ${\mathcal{H}}$ and $\widetilde{{\mathcal{H}}}$, respectively, relation (3.6) for the semigroups can be extended to far more general functions of the operators. The following result formulates this explicitly. It is of particular interest, for instance, when considering fractional Laplacians of the form $(-\Delta)^{\theta}$ with $\theta>1/2$, see Remark 4.2 (a) below. The result itself is probably well known and can be proved, for instance, with a standard reasoning using Stone’s formula. We give a brief variant of this reasoning below for the convenience of the reader.

Lemma 3.3**.**

Let $H$ and $\widetilde{H}$ be lower semibounded self-adjoint operators on Hilbert spaces ${\mathcal{H}}$ and $\widetilde{{\mathcal{H}}}$, respectively, and let $Y\colon\widetilde{{\mathcal{H}}}\to{\mathcal{H}}$ be a bounded linear operator such that $Y\widetilde{H}\subset HY$. Then, the spectral families $\mathsf{E}_{\tilde{H}}$ and $\mathsf{E}_{H}$ for $\tilde{H}$ and $H$, respectively, satisfy

[TABLE]

Moreover, for every Borel measurable function $\phi\colon\mathbb{R}\to\mathbb{R}$ the relation

[TABLE]

holds, where $\phi(H)$ and $\phi(\widetilde{H})$ are defined by functional calculus.

Proof.

It suffices to consider the case where ${\mathcal{H}}$ and $\widetilde{{\mathcal{H}}}$ are complex Hilbert spaces. The case of real Hilbert spaces can be obtained from this by complexification.

In view of the relation $Y\tilde{H}\subset HY$, we have $Y(\widetilde{H}-s\pm\mathrm{i}\varepsilon)\subset(H-s\pm\mathrm{i}\varepsilon)Y$ for every $s\in\mathbb{R}$ and $\varepsilon>0$, so that $Y(\widetilde{H}-s\pm\mathrm{i}\varepsilon)^{-1}g=(H-s\pm\mathrm{i}\varepsilon)^{-1}Yg$ for all $g\in\widetilde{{\mathcal{H}}}$ and, hence,

[TABLE]

Stone’s formula for the spectral families (see, e.g., [31, Satz 8.11]) then implies that

[TABLE]

for all $f\in{\mathcal{H}}$, $g\in\widetilde{{\mathcal{H}}}$, and $\lambda\in\mathbb{R}$, that is, $Y\mathsf{E}_{\widetilde{H}}(\lambda)=\mathsf{E}_{H}(\lambda)Y$ for all $\lambda\in\mathbb{R}$, which proves the first claim of the lemma.

Now, let $\phi\colon\mathbb{R}\to\mathbb{R}$ be a Borel measurable function and $g\in\operatorname{{\mathcal{D}}}(\phi(\widetilde{H}))$. The above then implies that $\langle\mathsf{E}_{H}(\lambda)Yg,Yg\rangle_{{\mathcal{H}}}\leq\|{Y}\|_{\widetilde{{\mathcal{H}}}\to{\mathcal{H}}}^{2}\langle\mathsf{E}_{\widetilde{H}}(\lambda)g,g\rangle_{{\mathcal{H}}}$, so that by functional calculus we conclude that $Yg\in\operatorname{{\mathcal{D}}}(H)$ and

[TABLE]

This proves the second claim and, hence, completes the proof. ∎

Remark 3.4*.*

In the case of complex Hilbert spaces in Lemma 3.3, also complex-valued Borel measurable functions $\phi$ can be considered, obviously, without any change to the proof.

4. Proof of the main theorem

We deduce Theorem 1.1 from the abstract result Theorem 3.1 with the choices ${\mathcal{H}}={\mathcal{U}}=L^{2}(\Omega)$, $\widetilde{{\mathcal{H}}}=\widetilde{{\mathcal{U}}}=L^{2}(\tilde{\Omega})$, $H=H_{\Omega}^{\bullet}$, $\widetilde{H}=H_{\tilde{\Omega}}^{\bullet}$, $B=\chi_{\omega}$, $\widetilde{B}=\chi_{\tilde{\omega}}$, and $p=2$. Set $\lambda^{N}:=1$ and $\lambda^{D}:=-1$, and for $\bullet\in\{D,N\}$ define

[TABLE]

with respect to the orthogonal decomposition $L^{2}(\tilde{\Omega})=L^{2}(\Omega)\oplus L^{2}(M(\Omega))$, where $J\colon L^{2}(\Omega)\rightarrow L^{2}(M(\Omega))$ is defined as $Jf:=f\circ M$.

The adjoint $(X^{\bullet})^{*}$ of $X^{\bullet}$ acts as

[TABLE]

for all $g\in L^{2}(\tilde{\Omega})$. Indeed, for $f\in L^{2}(\Omega)$ and $g\in L^{2}(\tilde{\Omega})$ we have

[TABLE]

by change of variables, so that

[TABLE]

It is then easy to see that

[TABLE]

Moreover, since $\tilde{\omega}=\omega\cup M(\omega)$ by construction, we have $\chi_{\tilde{\omega}}\circ M=\chi_{\tilde{\omega}}$ and, hence,

[TABLE]

Hence, with the choice $Y=Z=(X^{\bullet})^{*}$, condition (3.4) in Theorem 3.1 is satisfied, and we may take $\hat{Y}=X^{\bullet}/2$ as a bounded right inverse of $Y$ with $\|{\hat{Y}}\|\|{Z}\|=1$.

It remains to verify the extension relation (3.3) for the operators $H_{\Omega}^{\bullet}$ and $H_{\tilde{\Omega}}^{\bullet}$, that is, $(X^{\bullet})^{*}H_{\tilde{\Omega}}^{\bullet}\subset H_{\Omega}^{\bullet}(X^{\bullet})^{*}$. To this end, let us first briefly recall how the operators $H_{\Omega}^{\bullet}$ and $H_{\tilde{\Omega}}^{\bullet}$ are constructed. Since it is in general hard to define them directly by their differential expression, especially if the matrix function $A$ has discontinuities, the established approach to do that is via their quadratic forms: Consider the form ${\mathfrak{h}}_{\Omega}^{\bullet}={\mathfrak{a}}_{\Omega}^{\bullet}+V$ with

[TABLE]

where

[TABLE]

and

[TABLE]

The ellipticity condition (1.1) guarantees that this form ${\mathfrak{h}}_{\Omega}^{\bullet}$ is densely defined, lower semibounded, and closed, cf. [26, Example 4.19]. Hence, there is a self-adjoint operator $H_{\Omega}^{\bullet}=H_{\Omega}^{\bullet}(A,V)$ on $L^{2}(\Omega)$ given by

[TABLE]

and

[TABLE]

see, e.g., [15, Theorem VI.2.6], [25, Theorem 10.7], or [24, Theorem VIII.15]. The operator $H_{\tilde{\Omega}}^{\bullet}$ on $L^{2}(\tilde{\Omega})$ is defined completely analogous via the form ${\mathfrak{h}}_{\tilde{\Omega}}^{\bullet}={\mathfrak{a}}_{\tilde{\Omega}}^{\bullet}+\tilde{V}$ with

[TABLE]

and

[TABLE]

The desired extension relation $(X^{\bullet})^{*}H_{\tilde{\Omega}}^{\bullet}\subset H_{\Omega}^{\bullet}(X^{\bullet})^{*}$ can now be formulated in terms of the quadratic forms as

[TABLE]

provided that $X^{\bullet}f\in\operatorname{{\mathcal{D}}}({\mathfrak{h}}_{\tilde{\Omega}}^{\bullet})=\operatorname{{\mathcal{D}}}({\mathfrak{a}}_{\tilde{\Omega}}^{\bullet})$ and $(X^{\bullet})^{*}g\in\operatorname{{\mathcal{D}}}({\mathfrak{h}}_{\Omega}^{\bullet})=\operatorname{{\mathcal{D}}}({\mathfrak{a}}_{\Omega}^{\bullet})$; see the formal proof below.

In order to verify (4.6), we first take care of the mapping properties of $X^{\bullet}$ and $(X^{\bullet})^{*}$.

Lemma 4.1**.**

1. (a)

$(X^{\bullet})^{*}$* maps $\operatorname{{\mathcal{D}}}({\mathfrak{a}}_{\tilde{\Omega}}^{\bullet})$ into $\operatorname{{\mathcal{D}}}({\mathfrak{a}}_{\Omega}^{\bullet})$.* 2. (b)

$X^{\bullet}$* maps $\operatorname{{\mathcal{D}}}({\mathfrak{a}}_{\Omega}^{\bullet})$ into $\operatorname{{\mathcal{D}}}({\mathfrak{a}}_{\tilde{\Omega}}^{\bullet})$, and one has*

[TABLE]

for $f\in\operatorname{{\mathcal{D}}}({\mathfrak{a}}_{\Omega}^{\bullet})$ with $U=\operatorname{diag}(-1,1,\dots,1)$.

Proof.

(a). The case $\bullet=N$ is clear by (4.1). For $\bullet=D$ let first $g\in C_{c}^{\infty}(\tilde{\Omega})$. Then, $g-g\circ M\in C_{c}^{\infty}(\tilde{\Omega})$ and $g-g\circ M=0$ on $\Gamma$, and part (c) of Lemma A.1 in the appendix yields $(X^{D})^{*}g=(g-g\circ M)|_{\Omega}\in H_{0}^{1}(\Omega)$. The case of general $g\in H_{0}^{1}(\tilde{\Omega})$ follows from this by approximation.

(b). Let $f\in\operatorname{{\mathcal{D}}}({\mathfrak{a}}_{\Omega}^{\bullet})$ and $\varphi\in C_{c}^{\infty}(\tilde{\Omega})$, and abbreviate $\alpha_{1}:=-1$ and $\alpha_{k}:=1$ for $k\geq 2$. We then have to show that

[TABLE]

for $k=1,\dots,d$. To this end, observe that $(\partial_{k}\varphi)\circ M=\alpha_{k}\partial_{k}(\varphi\circ M)$, so that

[TABLE]

by change of variables. Now, for each $k\in\{1,\dots,d\}$ and $\bullet\in\{N,D\}$, parts (a) and (b) of Lemma A.1 allow to perform integration by parts without boundary terms on the right-hand side of (4.8). We therefore obtain

[TABLE]

where the last equality follows again by a change of variables. This equality together with (4.8) proves equation (4.7) and, hence, completes the proof. ∎

We are finally in position to prove our main result.

Proof of Theorem 1.1.

We need to verify identity (4.6). To this end, let $f\in\operatorname{{\mathcal{D}}}({\mathfrak{h}}_{\Omega}^{\bullet})=\operatorname{{\mathcal{D}}}({\mathfrak{a}}_{\Omega}^{\bullet})$ and $g\in\operatorname{{\mathcal{D}}}({\mathfrak{h}}_{\tilde{\Omega}}^{\bullet})=\operatorname{{\mathcal{D}}}({\mathfrak{a}}_{\tilde{\Omega}}^{\bullet})$. Since $U(\nabla g)\circ M=\nabla(g\circ M)$ and, by definition, $\tilde{A}(x)=UA(M(x))U$ for $x\in M(\Omega)$, we obtain from part (b) of Lemma 4.1 with change of variables that

[TABLE]

Thus,

[TABLE]

where we have taken into account (4.1) and part (a) of Lemma 4.1. Since also $\tilde{V}\circ M=\tilde{V}$ and, therefore, $(X^{\bullet})^{*}\tilde{V}=V(X^{\bullet})^{*}$, this proves (4.6). In turn, for $g\in\operatorname{{\mathcal{D}}}(H_{\tilde{\Omega}}^{\bullet})\subset\operatorname{{\mathcal{D}}}({\mathfrak{h}}_{\tilde{\Omega}}^{\bullet})$ we conclude that

[TABLE]

for all $f\in\operatorname{{\mathcal{D}}}({\mathfrak{h}}_{\Omega}^{\bullet})$, so that $(X^{\bullet})^{*}g\in\operatorname{{\mathcal{D}}}(H_{\Omega}^{\bullet})$ and $H_{\Omega}^{\bullet}(X^{\bullet})^{*}g=(X^{\bullet})^{*}H_{\tilde{\Omega}}^{\bullet}g$, which proves the desired extension relation (3.3). Since also (3.4) is satisfied by (4.3), applying Theorem 3.1 with $\hat{Y}=X^{\bullet}/2$ completes the proof. ∎

Remark 4.2*.*

(a) The above proof in particular shows the extension relation $(X^{\bullet})^{*}(-\Delta_{\tilde{\Omega}}^{\bullet})\subset(-\Delta_{\Omega}^{\bullet})(X^{\bullet})^{*}$. By Lemma 3.3 with the particular choice $\phi=(\cdot)^{\theta}\mathbbm{1}_{[0,\infty)}$, $\theta>1/2$, we then also have $(X^{\bullet})^{*}(-\Delta_{\tilde{\Omega}}^{\bullet})^{\theta}\subset(-\Delta_{\Omega}^{\bullet})^{\theta}(X^{\bullet})^{*}$. This allows us to apply Theorem 3.1 also in this situation, which gives an analogue of Theorem 1.1 for such fractional Laplacians, as claimed in Remark 2.10 (a).

(b) A modification of the above allows to consider also certain operators on $L^{p}$, $p\in(1,\infty)$. We demonstrate this briefly for the pure Laplacian on $L^{p}(\Omega)$, $\Omega=(0,\infty)\times\mathbb{R}^{d-1}$: Set $\tilde{\Omega}=\mathbb{R}^{d}$ and let $-\Delta_{\Omega,p}^{\bullet}$ and $-\Delta_{\tilde{\Omega},p}$ be the realizations on the differential expression $-\Delta$ on $L^{p}(\Omega)$ and $L^{p}(\tilde{\Omega})$, respectively, as closed sectorial operators, see, e.g., [6, Theorem 5.4 and Corollary 6.11]. These agree with $-\Delta_{\Omega}^{\bullet}=-\Delta_{\Omega,2}^{\bullet}$ and $-\Delta_{\tilde{\Omega}}=-\Delta_{\tilde{\Omega},2}$, respectively, on common elements of their domains. Consider now the bounded operator $Y^{\bullet}_{p}\colon L^{p}(\tilde{\Omega})\rightarrow L^{p}({\Omega})$, $Y_{p}^{\bullet}g=(g+\lambda^{\bullet}g\circ M)|_{\Omega}$, with $\|{Y_{p}^{\bullet}}\|_{L^{p}\rightarrow L^{p}}=2^{1/q}$, $1/p+1/q=1$, a right inverse of which is given by $\hat{Y}^{\bullet}_{p}\colon L^{p}({\Omega})\rightarrow L^{p}(\tilde{\Omega})$, $\hat{Y}^{\bullet}_{p}f=2^{-1}(f\oplus\lambda^{\bullet}f\circ M)$, with $\|{\hat{Y}_{p}^{\bullet}}\|_{L^{p}\rightarrow L^{p}}=2^{1/p-1}$. Let $g\in\mathcal{S}(\mathbb{R}^{d})$ be a Schwartz function. We then have $Y_{p}^{\bullet}g=(X^{\bullet})^{*}g\in\operatorname{{\mathcal{D}}}(-\Delta_{\Omega,p}^{\bullet})\cap\operatorname{{\mathcal{D}}}(\Delta_{\Omega,2}^{\bullet})$ and

[TABLE]

The extension relation $Y^{\bullet}_{p}(-\Delta_{\tilde{\Omega},p})\subset(-\Delta_{\Omega,p})Y^{\bullet}_{p}$ is then obtained by approximation since $\mathcal{S}(\mathbb{R}^{d})$ is an operator core for $-\Delta_{\tilde{\Omega},p}$. Theorem 3.1 can then be applied with $Y=Z=Y^{\bullet}_{p}$ and $\hat{Y}=\hat{Y}^{\bullet}_{p}$.

Remark 4.3*.*

In the above considerations, the operator $X^{\bullet}$ models the extension of a function in $L^{2}(\Omega)$ to a function in $L^{2}(\tilde{\Omega})$ by reflection with respect to one hyperplane. However, in view of the general form of Theorem 3.1, also other forms of extensions are feasible. For instance, one could consider to prolong functions on a sector of angle $\pi/(2n)$ in the plane to functions on an orthant by successive reflections with respect to different lines. In a similar way, in the work [23] the author considers the prolongation of functions on the equilateral triangle with corners $(0,0)$, $(0,1)$, $(0,1/\sqrt{3})$ to functions on the rectangle $(0,\sqrt{3})\times(0,1)$. One then has to prove results analogous to Lemma 4.1 and relation (4.9), allowing to infer null-controllability results on such sectors and the equilateral triangle from those on the orthant and the rectangle, respectively. In fact, in case of the equilateral triangle, these analogous results follow to some extend from the considerations in [23] already.

Appendix A An integration by parts formula

The following result is probably folklore, but in lack of a suitable reference we give an elementary proof using only integration by parts on hypercubes. We emphasize that we are not assuming any boundary regularity for the set $\tilde{\Omega}$.

Lemma A.1**.**

Let $\tilde{\Omega}\subset\mathbb{R}^{d}$ be an open set with $\Gamma:=\tilde{\Omega}\cap(\{0\}\times\mathbb{R}^{d-1})\neq\emptyset$, and set $\Omega:=\tilde{\Omega}\cap((0,+\infty)\times\mathbb{R}^{d-1})$. Let $f\in H^{1}(\Omega)$ and $\phi\in C_{c}^{\infty}(\tilde{\Omega})$.

1. (a)

One has

[TABLE]

for $k\in\{2,\dots,d\}$. 2. (b)

If, in addition, $f\in H_{0}^{1}(\Omega)$ or $\phi|_{\Gamma}=0$, then

[TABLE] 3. (c)

If $\phi|_{\Gamma}=0$, then $\phi|_{\Omega}\in H_{0}^{1}(\Omega)$.

Proof.

Since $\operatorname{supp}\phi\subset\tilde{\Omega}$ is compact, there is $\varepsilon>0$ with

[TABLE]

where the distance is taken with respect to the maximum norm. We now cover $\operatorname{supp}\phi\cap([0,+\infty)\times\mathbb{R}^{d-1})$ with suitable cubes from an equidistant lattice of mesh size $\varepsilon$. More precisely, we choose a finite subset $J\subset\mathbb{N}\times\mathbb{Z}^{d-1}$ with $\mathbb{N}=\{1,2,\dots\}$ such that the compact set

[TABLE]

is a subset of $\tilde{\Omega}\cap([0,+\infty)\times\mathbb{R}^{d-1})$ containing $\operatorname{supp}\phi\cap([0,+\infty)\times\mathbb{R}^{d-1})$ and such that $\operatorname{supp}\phi\cap((0,+\infty)\times\mathbb{R}^{d-1})$ lies in the interior $Q^{\circ}$ of $Q$. Clearly, we have $Q\setminus(\{0\}\times\mathbb{R}^{d-1})\subset\Omega$.

An integration by parts on each cube $\Lambda=\Lambda_{\varepsilon}(j)\subset\Omega$, $j\in J$, yields

[TABLE]

where $\Gamma_{k}^{+}$ and $\Gamma_{k}^{-}$ denote the two opposite faces of the cube $\Lambda$ (relatively open in the corresponding hyperplane) with outward unit normal in positive and negative direction with respect to the $k$-th coordinate axis, respectively. Now, each of the faces $\Gamma_{k}^{+}$ and $\Gamma_{k}^{-}$ either corresponds to the face of exactly one other cube $\Lambda_{\varepsilon}(l)$, $l\in J\setminus\{j\}$, but then with the opposite direction of the corresponding outward unit normal, or the face belongs to the boundary of $Q$. In the latter case, the integral over the face vanishes unless the face belongs to the hyperplane $\{0\}\times\mathbb{R}^{d-1}$, since $\operatorname{supp}\phi\cap((0,+\infty)\times\mathbb{R}^{d-1})$ lies in the interior of $Q$. This means that after summing over all cubes $\Lambda_{\varepsilon}(j)$, $j\in J$, in (A.1) only boundary integrals for faces in the hyperplane $\{0\}\times\mathbb{R}^{d-1}$ can survive, and these faces have an outward unit normal in the negative direction in the first coordinate ($k=1$).

Thus, after summing in (A.1) over all cubes $\Lambda_{\varepsilon}(j)$, $j\in J$, in both cases (a) and (b) we have no remaining boundary integrals, so that

[TABLE]

where for $k=1$ we may approximate $f\in H_{0}^{1}(\Omega)$ with functions from $C_{c}^{\infty}(\Omega)$ and then take the limit. Since $\operatorname{supp}\phi\cap((0,+\infty)\times\mathbb{R}^{d-1})\subset Q^{\circ}\subset\Omega$, this proves (a) and (b).

It remains to show that $\phi|_{\Omega}\in H_{0}^{1}(\Omega)$ if $\phi|_{\Gamma}=0$. For that, we choose an open cube $\tilde{\Lambda}\subset\mathbb{R}^{d}$ with $\operatorname{supp}\phi\subset\tilde{\Lambda}$ and $\Lambda^{\prime}:=\tilde{\Lambda}\cap((0,+\infty)\times\mathbb{R}^{d-1})\neq\emptyset$. Extending $\phi$ by zero on $\tilde{\Lambda}\setminus\tilde{\Omega}$, and taking into account that $\phi|_{\Gamma}=0$ by hypothesis, we have $\phi\in C_{c}^{\infty}(\tilde{\Lambda})$ and $\phi=0$ on $\partial\Lambda^{\prime}$. Since $\Lambda^{\prime}$ is convex and therefore has a Lipschitz boundary, see, e.g., [14, Corollary 1.2.2.3], this implies that $\phi|_{\Lambda^{\prime}}\in H_{0}^{1}(\Lambda^{\prime})$, see, e.g., [14, Corollary 1.5.1.6] or [1, Lemma A6.10]. Thus, there exists a sequence $(\varphi_{k})$ in $C_{c}^{\infty}(\Lambda^{\prime})$ such that $\varphi_{k}\to\phi|_{\Lambda^{\prime}}$ in $H^{1}(\Lambda^{\prime})$ as $k\to\infty$. We choose a smooth cutoff function $\eta\in C_{c}^{\infty}(\tilde{\Omega})$ with $\eta=1$ on $\operatorname{supp}\phi$. Since $(\eta\phi)|_{\Omega}=\phi|_{\Omega}$, it is straightforward to verify that $(\eta\varphi_{k})|_{\Omega}\to\phi|_{\Omega}$ in $H^{1}(\Omega)$ as $k\to\infty$. Since also $(\eta\varphi_{k})|_{\Omega}\in C_{c}^{\infty}(\Omega)$, we conclude that $\phi|_{\Omega}\in H_{0}^{1}(\Omega)$. This proves (c) and, hence, completes the proof of the lemma. ∎

Acknowledgements

The authors are grateful to Christoph Schumacher for an inspiring discussion and to Christian Seifert for a helpful correspondence. The authors also thank the anonymous referee of an earlier version of this manuscript for a useful comment that helped to improve this manuscript.