# Heat Equation with Inverse-Square Potential of Bridging Type Across Two Half-Lines

Matteo Gallone, Alessandro Michelangeli, Eugenio Pozzoli

► To cite this version:

Matteo Gallone, Alessandro Michelangeli, Eugenio Pozzoli. Heat Equation with Inverse-Square Potential of Bridging Type Across Two Half-Lines. Qualitative Properties of Dispersive PDEs, INdAM, 2021, Bari, Italy. pp.141-164, 10.1007/978-981-19-6434-3_7. hal-04309754

HAL Id: hal-04309754

https://hal.science/hal-04309754v1

Submitted on 3 Nov 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Heat equation with inverse-square potential of bridging type across two half-lines

Matteo Gallone*

Alessandro Michelangeli†

Eugenio Pozzoli‡

December 3, 2021

## Abstract

The heat equation with inverse square potential on both half-lines of $\mathbb{R}$ is discussed in the presence of bridging boundary conditions at the origin. The problem is the lowest energy (zero-momentum) mode of the transmission of the heat flow across a Grushin-type cylinder, a generalisation of an almost Riemannian structure with compact singularity set. This and related models are reviewed, and the issue is posed of the analysis of the dispersive properties for the heat kernel generated by the underlying positive self-adjoint operator. Numerical integration is shown that provides a first insight and relevant qualitative features of the solution at later times.

## 1 Introduction: the bridging heat equation in 1D.

For fixed $\alpha \in [0, 1)$ we discuss in this note the following initial value problem in the unknown $u \equiv u(t, x)$, with $t \ge 0$ and $x \in \mathbb{R} \setminus \{0\}$:

$$ egin{cases} 
\frac{\partial u}{\partial t} - \frac{\partial^2 u}{\partial x^2} + \frac{\alpha(\alpha + 2)}{4x^2} u = 0 \\ 
u_0(t)^- = u_0(t)^+ \quad \text{where } u_0^{\pm}(t) := \lim_{x \to 0^{\pm}} |x|^{\frac{\alpha}{2}} u(t, x), \\ 
u_1(t)^- = -u_1(t)^+ \quad \text{where } u_1^{\pm}(t) := \lim_{x \to 0^{\pm}} |x|^{-(1+\frac{\alpha}{2})} (u(t, x) - |x|^{-\frac{\alpha}{2}} u_0^{\pm}(t)), \\ 
u(0, x) = \varphi(x) \quad \text{where } \varphi \in L^2(\mathbb{R}), 
\end{cases} 
\quad (1) $$

seeking for solutions $u$ that, for (almost every) $t$, belong to $L^2(\mathbb{R})$. In the above formulation (1) the existence of the limits $u_0^{\pm}$ and $u_1^{\pm}$ is part of the problem. We shall also consider the special case where the initial datum itself satisfies the very boundary conditions required at later times.

We are in particular concerned with the well-posedness of the problem and the dispersive properties of the solution(s).

At the same time, in this note we review the origin and meaning of the problem (1) in the context of geometric quantum confinement or transmission across the metric's singularity for a particle constrained on a degenerate Riemannian manifold and only subject to the geometry of the constraining manifold, thus "free to evolve" over that manifold in analogy to a classical particle moving along geodesics.

<footnote>*Department of Mathematics, Milan University, via Cesare Saldini 50, 20133 Milan (Italy) e-mail: matteo.gallone@unimi.it</footnote>

<footnote>†Institute for Applied Mathematics, and Hausdorff Center for Mathematics, University of Bonn Endenicher Allee 60, D-53115 Bonn (Germany) e-mail: michelangeli@iam.uni-bonn.de</footnote>

<footnote>‡Institut de Mathématiques de Bourgogne, UMR 5584, CNRS, Université Bourgogne Franche-Comté, F-21000 Dijon (France) e-mail: eugenio.pozzoli@u-bourgogne.fr.</footnote>
The latter viewpoint is attracting an increasing amount of interest in recent years, making then natural to investigate the time-dependent equations arising in such contexts. Ours here is a ‘pilot’ analysis of a more systematic study that unfolds ahead of us concerning dispersive and Strichartz estimates, and it has therefore the purpose of some propaganda and overview of the state of the art and on the future perspectives. Besides, here we only deal with the heat evolution, and not the Schrödinger evolution, as we shall comment in due time.

Prior to outlining the geometric background, let us comment on the structure of the problem (1). The considered PDE is a heat-type equation governed by the second order, elliptic (Schrödinger) differential operator

$$ -\frac{d^2}{dx^2} + \frac{C_\alpha}{x^2}, \qquad C_\alpha := \frac{\alpha(\alpha + 2)}{4x^2} \qquad (2) $$

(the precise meaning of the parameter $\alpha$ and its presence through the coefficient $C_\alpha$ will be clear after discussing the parent geometric model). As such, the complete description of square-integrable solutions to the associated heat equation is achieved through a standard PDE analysis, once certain features of (2) are known as a linear operator on $L^2(\mathbb{R})$. For concreteness, a limit-circle/limit-point argument [23, Theorems X.11] shows that when $C_\alpha \ge \frac{3}{4}$, i.e., $\alpha \in (-\infty, -3] \cup [1, +\infty)$, the linear operator (2) minimally defined on smooth functions compactly supported away from $x = 0$ is actually essentially self-adjoint on $L^2(\mathbb{R})$. Denoting its closure with $A$, one concludes that $A$ is a self-adjoint operator with strictly positive spectrum and domain $\mathcal{D}(A)$ that explicitly, when $C_\alpha > \frac{3}{4}$, is the Sobolev space $H_0^2(\mathbb{R})$. As a straightforward consequence of the abstract theory of differential equation on Hilbert space [24, Proposition 6.6], one then concludes that the heat equation $\frac{d}{dt}u = -Au$ with initial datum $\varphi \in L^2(\mathbb{R})$ admits a unique solution in $C^1(\mathbb{R}_+, L^2(\mathbb{R}_x))$, with $u(t, \cdot) \in \mathcal{D}(A)$ at ever later $t > 0$, explicitly given by $u(t, x) = (e^{-tA}\varphi)(x)$.

In fact, it is worth recalling that the inverse square potential differential operator (2) is greatly studied and deeply understood from many standpoints, in particular, both as far radial space-time (Strichartz) estimates are concerned both in the linear and non-linear Schrödinger evolution (see, e.g., [9, 19] and references therein), and as a Bessel operator on the $L^2$-space of the half-line (see, e.g., [11] and its precursors in that prolific research line).

Yet, in addition to the differential side, the problem (1) prescribes the solutions $u$ to satisfy certain boundary conditions at $x = 0$. The first one, $g_0^+(t) = g_0^-(t)$ can be interpreted as the continuity of the function, up to the weight $|x|^{\frac{\alpha}{2}}$ that allows $u$ to have some degree of singularity at the origin; analogously, the condition $g_1^+(t) = -g_1^-(t)$ links the right and left derivative at zero, up to certain weights, and taken directionally from each side. In the regime $C_\alpha > \frac{3}{4}$ such conditions are obviously redundant, but when $C_\alpha \le \frac{3}{4}$ an ad hoc analysis is needed to recognise that the prescribed behaviour at the origin expresses another condition of self-adjointness and positivity. Such an analysis has been carried out in several recent works [7, 17, 18, 20, 15] and is concisely reviewed in Section 2. The net result, for the sake of the present discussion, is the following.

**Theorem 1** ([18]). Let $\alpha \in [0, 1)$ and let $C_\alpha$ be given by (2).

(i) The space

$$ \mathcal{D} := \left\{ g \in L^2(\mathbb{R}) \left| \left( -\frac{d^2}{dx^2} + \frac{C_\alpha}{x^2} \right) g \in L^2(\mathbb{R}) \right. \right\} $$

is a dense subspace of $L^2(\mathbb{R})$ and for every $g \in \mathcal{D}$ the following limits exist and are finite:

$$ \begin{aligned} g_0^{\pm} &= \lim_{x \to 0^{\pm}} |x|^{\frac{\alpha}{2}} g^{\pm}(x), \\ g_1^{\pm} &= \lim_{x \to 0^{\pm}} |x|^{-(1+\frac{\alpha}{2})} (g^{\pm}(x) - g_0^{\pm}|x|^{-\frac{\alpha}{2}}). \end{aligned} \qquad (3) $$
<figure><img src="image_2.png" /><figcaption>Figure 1: Manifold $M_\alpha$ for different $\alpha$: $\alpha < 0$ (left), $\alpha = 0$ (center) and $\alpha > 0$ (right).</figcaption></figure>

(ii) The operator

$$ \begin{aligned} \mathcal{D}(A_{\alpha}^{\mathrm{B}}) &= \{g \in \mathcal{D} | g_{0}^{+} = g_{0}^{-}, g_{1}^{+} = -g_{1}^{-}\}, \\ A_{\alpha}^{\mathrm{B}}g &= -g'' + C_{\alpha}|x|^{-2}g \end{aligned} \qquad (4) $$

is self-adjoint on $L^2(\mathbb{R})$ and non-negative. Its spectrum is $[0, +\infty)$, and is all essential and absolutely continuous.

In Theorem 1 we only consider the regime $\alpha \in [0, 1)$. The remaining regime $\alpha \in [-3, 0)$ is simply less relevant from the viewpoint of the underlying geometric model, as will be argued in Section 2. And, as discussed above, when $\alpha \in (-\infty, -3) \cup (1, +\infty)$ one applies standard limit-point/limit-circle considerations.

In view of Theorem 1, the initial value problem (1) is immediately interpreted as the problem for the one-dimensional heat equation governed by the positive and self-adjoint operator $A_\alpha^B$, and therefore it admits unique solution $u(t) = e^{-tA_\alpha^B}\varphi$, again by abstract facts of differential equations on Hilbert space [24, Proposition 6.6].

The well-posedness of (1) is therefore fully controlled in $C^1(\mathbb{R}_t^+, L^2(\mathbb{R}_x))$ with $u(t, \cdot) \in \mathcal{D}(A_\alpha^B)$ at every $t > 0$.

The superscript ‘B’ in $A_\alpha^B$ is to refer to certain ‘bridging’ features of optimal transmission across the origin, allowing in a precise sense complete communication between the right and left half-line, induced by $A_\alpha^B$, as compared to a whole family of similar transmission protocols.

Indeed, in Section 2 it will be recapped how initial value problems like (1), and their counterparts with the Schrödinger equation, describe the heat-type or Schrödinger-type propagation over a particular almost-Riemannian structure, customarily referred to as a ‘Grushin cylinder’, constituted by an infinite two-dimensional cylinder with a non-flat metric that becomes suitably singular on a given orthogonal section. Depending on the magnitude of the metric’s singularity, which is quantified by the parameter $\alpha$, the transmission is either inhibited, so that a function initially supported on one half-cylinder remains confined in that half at later times, or on the opposite it is allowed, through a precise set of boundary conditions between the two halves. This can be qualitatively visualised as in Figure 1 with cylinders that shrink to one point or get flattened in correspondence of the given singular section. One-dimensional problems like (1) emerge for the evolution on the lowest-energy mode, which corresponds to functions on the cylinder that are constant along the compact variable. It turns out that in certain regimes of metric’s singularity (and $\alpha \in [0, 1)$ is the physically most significant regime) an infinity of transmission protocols emerge, each characterised by suitable boundary conditions of self-adjointness, and each yielding a heat-type or Schrödinger-type equation. Among them, the ‘bridging’ protocol described by the operator $A_\alpha^B$ given by (4) displays distinguished features of optimal transmission. In practice – see equations (30)-(31) below – the heat (and Schrödinger) equation of ‘bridging’ type between the two half-lines is the one that
describes a crossing at $x = 0$ 'without spatial filter' (continuity of the function) and 'without energy filter' (the fraction of transmitted flux does not depend on the incident energy).

The bridging protocol $A_\alpha^B$ was first identified in the recent work [7]. The comparison analysis of the bridging protocol with respect to the whole family of the other physically meaningful ones was subsequently analysed in [17, 15]. Section 3 reports on the recent literature of closely related models and results. Moreover, in [7, 6] the bridging and some other protocols were analysed from the perspective of the conservation of the total heat, or equivalently, the perspective of infinite lifespan of the stochastic processes generated by such operators (stochastic completeness).

It should be then sufficiently clear at this point that the initial value problem (1) describes the low-energy transmission of bridging-type between the two halves of a Grushin cylinder with metric singularity at zero. In view of that, beside the already guaranteed well-posedness, *dispersive properties* of the solution come to have great relevance in connection with the underlying physical transmission protocol, in particular $L^p-L^q$ estimates, smoothing estimates, and space-time (Strichartz) estimates for the heat semi-group associated with $A_\alpha^B$.

These are the analogue of the well-known estimates for the classical heat equation $(\frac{\partial}{\partial t} - \frac{\partial^2}{\partial x^2})u(t,x) = 0$, $u(0,x) = \varphi(x)$, for which one has [25, Section 2.2]

$$ 
\begin{align} 
\|u(t, \cdot)\|_{L^p} &\lesssim t^{-\frac{1}{2}(\frac{1}{r}-\frac{1}{p})} \|\varphi\|_{L^r} &\quad 1 \le r \le p \le \infty, \nonumber \\
\left\| \frac{\partial}{\partial x} u(t, \cdot) \right\|_{L^p} &\lesssim t^{-\frac{1}{2}(\frac{1}{r}-\frac{1}{p}+1)} \|\varphi\|_{L^r} &\quad 1 \le r \le p \le \infty, \nonumber \\
\|u\|_{L^q(\mathbb{R}_t^+, L_x^p)} &\lesssim \|\varphi\|_{L^2} &\quad 2 \le p < \infty, \quad q = \frac{4p}{p-2}, \tag{5} \\
\left\| \frac{\partial}{\partial x} u \right\|_{L^2(\mathbb{R}_t^+, L_x^2)} &\lesssim \|\varphi\|_{L^2}. \nonumber
\end{align} 
$$

Establishing the analogue of (5) for the heat type semi-group $\exp(-tA_\alpha^B)$, and in fact eventually also for the Schrödinger unitary group $\exp(-itA_\alpha^B)$, as well as for the corresponding semi-groups and groups induced by other transmission protocols, and more generally on other geometries and classes of almost-Riemannian structures beyond the Grushin cylinder, appear to be one of most relevant challenges in this field, with abstract interest per se and impact on applications, of quantum control concern in the first place.

In this respect we intend with this note to promote the above questions and advertise them for future investigations, in particular posing them in the rigorous context of geometric confinement and transmission protocols. Such an overview is given, as mentioned, in Sections 2 and 3.

Last, in Section 4 we present a glance at numerical computations of the solution $u$ of (1) when the initial datum is well localised on one half-line. The numerical evidence is strong on the dynamical formation of the bridging boundary conditions at $x = 0$ at later times, and on a general behaviour that is qualitatively comparable to the classical heat propagation.

In fact, at present no analytic computation is available of the heat (and, in the future, the Schrödinger) propagator generated by $A_\alpha^B$, and numerics is a first natural approach to infer meaningful properties to be rigorously proved in forthcoming investigations.

## 2 A concise review of geometric confinement and transmission protocols in a Grushin cylinder

We have already anticipated that the problem (1) provides the one-dimensional description for the heat flow across the singularity of a Grushin cylinder, and let us give in this Section a concise overview of the problem from that perspective.
Grushin cylinders are Riemannian manifolds $M_\alpha \equiv (M, g_\alpha)$, with parameter $\alpha \in \mathbb{R}$, where

$$ M^{\pm} := \mathbb{R}_{x}^{\pm} \times \mathbb{S}_{y}^{1}, \quad \mathcal{Z} := \{0\} \times \mathbb{S}_{y}^{1}, \quad M := M^{+} \cup M^{-} \quad (6) $$

and with degenerate Riemannian metric

$$ g_{\alpha} := dx \otimes dx + |x|^{-2\alpha} dy \otimes dy. \quad (7) $$

Thus, $M_\alpha$ is a two-dimensional manifold built upon the cylinder $\mathbb{R} \times \mathbb{S}^1$, with singularity locus $\mathcal{Z}$ and incomplete Riemannian metric both on the right and the left half-cylinder $\mathbb{R}^{\pm} \times \mathbb{S}^1$ meaning that geodesics cross smoothly the singularity $\mathcal{Z}$ at finite times). The values $\alpha = -1, \alpha = 0$, and $\alpha = 1$ select, respectively, the flat cone, the Euclidean cylinder, and the standard 'Grushin cylinder' [10, Chapter 11]: in the latter case one has an 'almost-Riemannian structure' on $\mathbb{R} \times \mathbb{S}^1 = M^+ \cup \mathcal{Z} \cup M^-$ in the rigorous sense of [2, Sec. 1] or [22, Sect. 7.1]. Actually, $g_\alpha$ is defined as the unique metric for which the distribution of vector fields globally defined on $\mathbb{R} \times \mathbb{S}^1$ as

$$ X_1(x, y) := \frac{\partial}{\partial x}, \quad X_2^{(\alpha)}(x, y) := |x|^{\alpha} \frac{\partial}{\partial y} \quad (8) $$

is an orthonormal frame at every $(x, y) \in \mathbb{R} \times \mathbb{S}^1$: in this regard, the Grushin cylinder ($\alpha = 1$) is a two-dimensional almost-Riemannian manifold of step two, meaning that

$$ \operatorname{span}\{X_1, X_2^{(1)}, [X_1, X_2^{(1)}]\}|_{(x,y)} = \mathbb{R}^2 \quad \forall (x, y) \in \mathbb{R} \times \mathbb{S}^1, $$

where $[X_1, X_2^\alpha]$ denotes the Lie brackets of vector fields. In fact, $M_\alpha$ is a hyperbolic manifold whenever $\alpha > 0$, with Gaussian (sectional) curvature

$$ K_{\alpha}(x, y) = - \frac{\alpha(\alpha + 1)}{x^2}. \quad (9) $$

To each $M_\alpha$ one naturally associates the Riemannian volume form

$$ \mu_{\alpha} := \operatorname{vol}_{g_{\alpha}} = \sqrt{\det g_{\alpha}} \, dx \wedge dy = |x|^{-\alpha} \, dx \wedge dy, \quad (10) $$

the Hilbert space

$$ \mathcal{H}_{\alpha} := L^2(M, d\mu_{\alpha}), \quad (11) $$

understood as the completion of $C_c^\infty(M)$ with respect to the scalar product

$$ \langle \psi, \varphi \rangle_{\alpha} := \iint_{\mathbb{R} \times \mathbb{S}^1} \overline{\psi(x, y)} \varphi(x, y) \frac{1}{|x|^\alpha} \, dx \, dy, \quad (12) $$

and the (Riemannian) Laplace-Beltrami operator $\Delta_{\mu_\alpha} := \operatorname{div}_{\mu_\alpha} \circ \nabla$ acting on functions over $M_\alpha$. A standard computation (see, e.g. [17, Sect. 2]) yields explicitly

$$ \Delta_{\mu_{\alpha}} = \frac{\partial^2}{\partial x^2} + |x|^{2\alpha} \frac{\partial^2}{\partial y^2} - \frac{\alpha}{|x|} \frac{\partial}{\partial x}. \quad (13) $$

As a linear operator on $\mathcal{H}_\alpha$ one minimally defines $\Delta_{\mu_\alpha}$ on the dense subspace of $L^2(M, d\mu_\alpha)$-functions that are smooth and compactly supported away from the metric's singularity locus $\mathcal{Z}$, thus introducing

$$ H_{\alpha} := -\Delta_{\mu_{\alpha}}, \quad \mathcal{D}(H_{\alpha}) := C_{c}^{\infty}(M). \quad (14) $$

The Green identity implies that $H_\alpha$ is symmetric and non-negative. One has the following classification.
**Theorem 2** ([5, 7, 17]). (i) If $\alpha \in (-\infty, -3] \cup [1, +\infty)$, then the operator $H_\alpha$ is essentially self-adjoint.

(ii) If $\alpha \in (-3, -1]$, then $H_\alpha$ is not essentially self-adjoint with deficiency index 2.

(iii) If $\alpha \in (-1, +1)$, then $H_\alpha$ is not essentially self-adjoint and it has infinite deficiency index.

The same holds, separately, for the symmetric operators $H_\alpha^\pm$ minimally defined on the $L^2$-space of each half-cylinder.

With respect to the Hilbert space orthogonal decomposition

$$ \mathcal{H}_{\alpha} = L^{2}(M, d\mu_{\alpha}) \cong L^{2}(M^{-}, d\mu_{\alpha} \oplus L^{2}(M^{+}, d\mu_{\alpha})) \quad (15) $$

the operator $H_\alpha$ is reduced as $H_\alpha = H_\alpha^- \oplus H_\alpha^+$, and therefore in the regime of essential self-adjointness its closure is the reduced, non-negative, self-adjoint operator $\overline{H_\alpha} = \overline{H_\alpha^-} \oplus \overline{H_\alpha^+}$. This implies that both the Schrödinger equation $\partial_t u = -i\overline{H_\alpha}u$ and the heat equation $\partial_t u = -\overline{H_\alpha}u$ decompose to uncoupled equations on each half-cylinder, or, better to say, group and semi-group decompose, respectively, as $e^{-it\overline{H_\alpha}} = e^{-it\overline{H_\alpha^-}} \oplus e^{-it\overline{H_\alpha^+}}$ and $e^{-t\overline{H_\alpha}} = e^{-t\overline{H_\alpha^-}} \oplus e^{-t\overline{H_\alpha^+}}$, with the consequence that an initial datum supported, say, only on one half, keeps evolving in that half at each later time. This phenomenon is customarily referred to as 'heat-geometric confinement' and 'quantum-geometric confinement', respectively, to emphasise the sole effect of the geometry (meaning that we are not considering any potential energy on the manifold, but only the kinetic one), with no coupling boundary conditions – hence no interaction – declared at $\mathcal{Z}$. Quantum-mechanically, in this regime of the Grushin metric, a quantum particle constrained on $M_\alpha$ and left 'free' to evolve only under the effect of the underlying geometry, never happen to cross the singularity locus $\mathcal{Z}$.

The scenario becomes much more diversified when $H_\alpha$ is not essentially self-adjoint and therefore admits non-trivial self-adjoint extensions. Our regime of interest includes $\alpha \in (0, 1)$, the sub-case of greatest relevance because it corresponds to an actual local singularity (and not vanishing) of the metric $g_\alpha$, and for the purposes of this note we shall only consider such $\alpha$'s. Qualitatively analogous results can be established in the remaining non-self-adjoint regime $\alpha \in (-3, 0)$.

There is in fact a giant family of inequivalent self-adjoint realisations of $H_\alpha$ when $\alpha \in [0, 1)$, as the deficiency index is infinite. Each one is characterised by boundary conditions at $\mathcal{Z}$ that prescribe a one-sided or two-sided interaction with the boundary, or more generally a protocol of left↔right transmission. Such a family include physically unstable realisations (those that are not lower semi-bounded), as well as a huge amount of unphysical realisations, such as those with non-local boundary conditions at $\mathcal{Z}$.

An extensive and fairly explicit classification of physical extensions of $H_\alpha$ was recently completed in [18].

**Theorem 3** ([18]). Let $\alpha \in [0, 1)$. $H_\alpha$ defined in (14) admits, among others, the following families of self-adjoint extensions with respect to $L^2(M, d\mu_\alpha)$:

* Friedrichs extension: $H_{\alpha,F}$;

* Family $I_R$: $\{H_{\alpha,R}^{[\gamma]} | \gamma \in \mathbb{R}\}$;

* Family $I_L$: $\{H_{\alpha,L}^{[\gamma]} | \gamma \in \mathbb{R}\}$;

* Family $\Pi_\alpha$ with $a \in \mathbb{C}$: $\{H_{\alpha,a}^{[\gamma]} | \gamma \in \mathbb{R}\}$;

* Family III: $\{H_\alpha^{[\Gamma]} | \Gamma = (\gamma_1, \gamma_2, \gamma_3, \gamma_4) \in \mathbb{R}^4\}$.

Each member of any such family acts precisely as the differential operator $-\Delta_{\mu_\alpha}$ on a domain of functions $f \in L^2(M, d\mu_\alpha)$ satisfying the following properties.
(i) *Integrability and regularity:*

$$ \sum_{\pm} \iint_{\mathbb{R}_x^{\pm} \times \mathbb{S}_y^1} |(\Delta_{\mu_\alpha} f)(x, y)|^2 d\mu_\alpha(x, y) < +\infty. \quad (16) $$

(ii) *Boundary condition: The limits*

$$ f_0^{\pm}(y) = \lim_{x \to 0^{\pm}} f(x, y) \quad (17) $$

$$ f_1^{\pm}(y) = \pm(1+\alpha)^{-1} \lim_{x \to 0^{\pm}} \left( \frac{1}{|x|^\alpha} \frac{\partial f(x, y)}{\partial x} \right) \quad (18) $$

exist and are finite for almost every $y \in \mathbb{S}^1$, and depending on the considered type of extension, and for almost every $y \in \mathbb{R}$,

$$ f_0^{\pm}(y) = 0 \quad \text{if } f \in \mathcal{D}(H_{\alpha,F}), \quad (19) $$

$$ \begin{cases} f_0^-(y) = 0 \\ f_1^+(y) = \gamma f_0^+(y) \end{cases} \quad \text{if } f \in \mathcal{D}(H_{\alpha,R}^{[\gamma]}), \quad (20) $$

$$ \begin{cases} f_1^-(y) = \gamma f_0^-(y) \\ f_0^+(y) = 0 \end{cases} \quad \text{if } f \in \mathcal{D}(H_{\alpha,L}^{[\gamma]}), \quad (21) $$

$$ \begin{cases} f_0^+(y) = a f_0^-(y) \\ f_1^-(y) + \bar{a} f_1^+(y) = \gamma f_0^-(y) \end{cases} \quad \text{if } f \in \mathcal{D}(H_{\alpha,a}^{[\gamma]}), \quad (22) $$

$$ \begin{cases} f_1^-(y) = \gamma_1 f_0^-(y) + (\gamma_2 + i\gamma_3) f_0^+(y) \\ f_1^+(y) = (\gamma_2 - i\gamma_3) f_0^-(y) + \gamma_4 f_0^+(y) \end{cases} \quad \text{if } f \in \mathcal{D}(H_{\alpha}^{[\Gamma]}). \quad (23) $$

One can further select those extensions that are non-negative and then induce a heat type flow.

**Theorem 4** ([15]). • The Friedrichs extension $H_{\alpha,F}$ is non-negative.

• Extensions in the family $I_R$, $I_L$, and $\Pi_a$, $a \in \mathbb{C}$, are non-negative if and only if $\gamma \ge 0$.

• Extensions in the family III are non-negative if and only if so is the matrix

$$ \tilde{\Gamma} := \begin{pmatrix} \gamma_1 & \gamma_2 + i\gamma_3 \\ \gamma_2 - i\gamma_3 & \gamma_4 \end{pmatrix}, $$

i.e., if and only if $\gamma_1 + \gamma_4 > 0$ and $\gamma_1\gamma_4 \ge \gamma_2^2 + \gamma_3^2$.

A customary *quantum-mechanical* quantification of the transmission modelled by each extension is the fraction of Schrödinger flux that gets transmitted vs reflected when a beam of particles are shot free from infinity towards $\mathcal{Z}$. This analysis, albeit in a *Schrödinger equation* framework, elucidates the qualitative properties of the crossing at $x = 0$ and was recently done in [15]. Intuitively speaking, far away from $\mathcal{Z}$ the metric tends to become flat and the action $-\Delta_{\mu_\alpha}$ of each self-adjoint ‘free Hamiltonian’ tends to resemble that of the free Laplacian $-\Delta$, plus the correction due to the $(|x|^{-1}\partial_x)$-term, on wave functions $f(x, y)$ that are constant in $y$. This suggests that at very large distances a quantum particle evolves free from the effects of the underlying geometry, and one can speak of scattering states of energy $E > 0$. The precise shape of the wave function $f_{\text{scatt}}$ of such a scattering state can be easily guessed to be of the form

$$ f_{\text{scatt}}(x, y) \sim |x|^{\frac{\alpha}{2}} e^{\pm i x \sqrt{E}} \quad \text{as } |x| \to +\infty. \quad (24) $$
Indeed, $-\Delta_{\mu_a} f_{\text{scatt}} \sim E f_{\text{scatt}} + \frac{\alpha(2+\alpha)}{4|x|^2} f_{\text{scatt}}$, that is, up to a very small $O(|x|^{-2})$-correction, $f_{\text{scatt}}$ is a generalised eigenfunction of $-\Delta_{\mu_a}$ with eigenvalue $E$. All this can be fully justified on rigorous grounds [15] and leads naturally to the definition of the 'transmission coefficient' and 'reflection coefficient' for the scattering, namely the spatial density of the transmitted flux and the reflected flux, normalised with respect to the density of the incident flux. Obviously, no scattering across the singularity occurs for Friedrichs, or type-I<sub>R</sub>, or type-I<sub>L</sub> quantum protocols, whereas in type-II<sub>a</sub> scattering one obtains the following (analogous conclusions can be made for type-III scattering).

**Theorem 5** ([15]). Let $\alpha \in [0, 1)$, $a \in \mathbb{C}$, $\gamma \in \mathbb{R}$. The transmission coefficient $T_{\alpha,a,\gamma}(E)$ and the reflection coefficient $R_{\alpha,a,\gamma}(E)$ at given energy $E > 0$ for the Schrödinger transmission protocol governed by $H_{\alpha,a}^{[\gamma]}$ are given by

$$ T_{\alpha,a,\gamma}(E) = \left| \frac{E^{\frac{1+\alpha}{2}} (1 + e^{i\pi\alpha}) \Gamma(\frac{1-\alpha}{2}) \bar{a}}{E^{\frac{1+\alpha}{2}} \Gamma(\frac{1-\alpha}{2}) (1 + |a|^2) + i\gamma 2^{1+\alpha} e^{i\frac{\pi}{2}\alpha} \Gamma(\frac{3+\alpha}{2})} \right|^2, \quad (25) $$

$$ R_{\alpha,a,\gamma}(E) = \left| \frac{E^{\frac{1+\alpha}{2}} \Gamma(\frac{1-\alpha}{2}) (1 - |a|^2 e^{i\pi\alpha}) + i\gamma 2^{1+\alpha} e^{i\frac{\pi}{2}\alpha} \Gamma(\frac{3+\alpha}{2})}{E^{\frac{1+\alpha}{2}} \Gamma(\frac{1-\alpha}{2}) (1 + |a|^2) + i\gamma 2^{1+\alpha} e^{i\frac{\pi}{2}\alpha} \Gamma(\frac{3+\alpha}{2})} \right|^2. $$

They satisfy

$$ T_{\alpha,a,\gamma}(E) + R_{\alpha,a,\gamma}(E) = 1, \quad (26) $$

and when $\gamma = 0$ they are independent of $E$. The scattering is reflection-less ($R_{\alpha,a,\gamma}(E) = 0$) when

$$ E = \left( \frac{2^{1+\alpha} \gamma \Gamma(\frac{3+\alpha}{2}) \sin \frac{\pi}{2}\alpha}{\Gamma(\frac{1-\alpha}{2})(1 - \cos \pi\alpha)} \right)^{\frac{2}{1+\alpha}}, \quad (27) $$

provided that $\alpha \in (0, 1)$, $|a| = 1$, and $\gamma > 0$. In the high energy limit the scattering is independent of the extension parameter $\gamma$ and one has

$$ \lim_{E \to +\infty} T_{\alpha,a,\gamma}(E) = \frac{2|a|^2(1 + \cos \pi\alpha)}{(1 + |a|^2)^2}, \quad (28) $$

$$ \lim_{E \to +\infty} R_{\alpha,a,\gamma}(E) = \frac{1 + |a|^4 - 2|a|^2 \cos \pi\alpha}{(1 + |a|^2)^2}, $$

whereas in the low energy limit, for $\gamma \neq 0$,

$$ \lim_{E \downarrow 0} T_{\alpha,a,\gamma}(E) = 0, \quad (29) $$

$$ \lim_{E \downarrow 0} R_{\alpha,a,\gamma}(E) = 1. $$

Upon inspection of the boundary conditions (19)-(23) one sees that the type-II<sub>a</sub> extension $H_{\alpha,a}^{[\gamma]}$ with $a = 1$ and $\gamma = 0$ imposes the local behaviour

$$ \lim_{x \to 0^{-}} f(x, y) = \lim_{x \to 0^{+}} f(x, y) \\ \lim_{x \to 0^{-}} \left( \frac{1}{|x|^{\alpha}} \frac{\partial f(x, y)}{\partial x} \right) = \lim_{x \to 0^{+}} \left( \frac{1}{|x|^{\alpha}} \frac{\partial f(x, y)}{\partial x} \right), \quad (30) $$

namely the distinguished feature of having a domain of functions that are continuous across the Grushin singularity, together with their weighted derivative. $H_{\alpha,1}^{[0]}$ is called the 'bridging extension' of $H_{\alpha}$, and for it we shall simply write $H_{\alpha}^{\text{B}}$. In view of the results reviewed so far, the transmission modelled by the bridging extension
• has no spatial filter in the sense of (30) (in fact, all type-II$_{a}$ protocols with a = 1 impose local continuity at Z ; quantum-mechanically this is interpreted as a lack of jump in the particle’s probability density from one side to the other of the singularity),

• and has no energy filter in the Schrödinger scattering, indeed, H B α and all type-II$_{a}$ protocols with γ = 0 induce a scattering where the fraction of transmitted and reflected flux does not depend on the incident energy (see (25) above),

$$
T_{\alpha}^{B} := T_{\alpha,1,0}(E) = \frac{1}{2} (1 + \cos \pi\alpha), \qquad (31)
$$

$$
R_{\alpha}^{B} := R_{\alpha, 1, 0}(E) = \frac{1}{2} (1 - \cos \pi\alpha), \qquad (31)
$$

meaning that the singularity does not act as a filter in the energy.

The overall picture surveyed so far poses naturally the question of the analysis of the heat type flow generated by the positive and self-adjoint realisations of the Laplace-Beltrami operator on H α (Theorem 4), as well as the Schrödinger type flow generated by self-adjoint realisations (Theorem 3), let alone the study of non-linear heat and Schrödinger equations on M α with linear part given by a self-adjoint Laplace-Beltrami operator. This appears to be a completely uncharted territory of notable relevance in abstract terms and for applications. The gap between such future goals and the current knowledge is a lack of informative characterisation of the heat and Schrödinger propagator’s kernel.

To complete the present review, let us make the connection explicit between the two-dimensional heat type equation induced by H B α and the one-dimensional problem (1).

This is done [7, 18] by means of the canonical Hilbert space unitary isomorphism H α ≅ → H , where

$$
\mathcal{H} := \mathcal{F}_2 U_\alpha L^2(M, d\mu_\alpha) \cong \ell^2(\mathbb{Z}, L^2(\mathbb{R}, dx)) \cong \mathcal{H}^- \oplus \mathcal{H}^+ \cong \bigoplus_{k \in \mathbb{Z}} \mathfrak{h}, \quad (32)
$$

$$
\mathfrak{h} := L^2(\mathbb{R}^{-}, dx) \oplus L^2(\mathbb{R}^{+}, dx) \cong L^2(\mathbb{R}, dx),
$$

recalling that

$$
\mathcal{H}_{\alpha} \cong L^2(M^{-}, d\mu_{\alpha}) \oplus L^2(M^{+}, d\mu_{\alpha}), \quad (33)
$$

and where the unitary transformations U$_{α}$ := U − α ⊕ U + α and F$_{2}$ := F − 2 ⊕ F + 2 are defined, respectively, as

$$
U_{\alpha}^{\pm} : L^2(\mathbb{R}^{\pm} \times \mathbb{S}^1, |x|^{-\alpha} dxdy) \xrightarrow{\cong} L^2(\mathbb{R}^{\pm} \times \mathbb{S}^1, dxdy) \quad (34)
$$

$$
f \mapsto \phi := |x|^{-\frac{\alpha}{2}} f,
$$

and

$$
\begin{align}
\mathcal{F}_2^{\pm} : L^2(\mathbb{R}^{\pm} \times \mathbb{S}^1, dxdy) &\xrightarrow{\cong} L^2(\mathbb{R}^{\pm}, dx) \otimes \ell^2(\mathbb{Z}), \nonumber \\
\phi &\mapsto \psi \equiv (\psi_k)_{k \in \mathbb{Z}}, \tag{35} \\
e_k(y) := \frac{e^{iky}}{\sqrt{2\pi}}, \quad \psi_k(x) := \int_0^{2\pi} e_k(y) \phi(x, y) dy, \quad x \in \mathbb{R}^{\pm} \nonumber
\end{align}
$$

(thus, φ ( x, y ) = ∑ k ∈ Z ψ$_{k}$ ( x ) e$_{k}$ ( y ) in the L $^{2}$-convergent sense). This provides, up to isomorphism, the orthogonal sum decomposition of the Hilbert space of interest into identical ‘bilateral’ fibres h = L $^{2}$( R $^{−}$, dx ) ⊕ L $^{2}$( R $^{+}$, dx ) ≅ L $^{2}$( R , dx ). The decomposition is discrete, as a consequence of having taken the Fourier transform F$_{2}$ only in the compact variable y .
**Theorem 6** ([18]). Let $\alpha \in [0, 1)$. Through the isomorphism (32) the self-adjoint bridging operator $H_{\alpha}^{B}$ on $\mathcal{H}_{\alpha} = L^{2}(M, d\mu_{\alpha})$ is unitarily equivalent to the self-adjoint operator $\mathcal{H}_{\alpha}^{B}$ on $\mathcal{H} \cong \ell^{2}(\mathbb{Z}, L^{2}(\mathbb{R}))$, namely

$$H_{\alpha}^{B} = (U_{\alpha})^{-1}(\mathcal{F}_{2})^{-1}\mathcal{H}_{\alpha}^{B} \mathcal{F}_{2} U_{\alpha}, \qquad (36)$$

where

$$\mathcal{H}_{\alpha}^{B} = \bigoplus_{k \in \mathbb{Z}} A_{\alpha}(k) \qquad (37)$$

and each $A_{\alpha}(k)$ is the self-adjoint operator on $L^{2}(\mathbb{R})$ given by

$$\mathcal{D}(A_{\alpha}(k)) = \left\{ \begin{array}{ll} g = g^{-} \oplus g^{+}, & g^{\pm} \in L^{2}(\mathbb{R}^{\pm}, \mathrm{d}x) \text{ such that} \\ \left(-\frac{d^2}{dx^2} + k^2|x|^{2\alpha} + \frac{\alpha(2+\alpha)}{4x^2}\right)g^{\pm} \in L^{2}(\mathbb{R}^{\pm}, \mathrm{d}x) \\ g_0^- = g_0^+, \quad g_1^- = -g_1^+ \end{array} \right\}, \qquad (38)$$

$$A_{\alpha}(k)g = \bigoplus_{\pm} \left( -\frac{d^2}{dx^2} + k^2|x|^{2\alpha} + \frac{\alpha(2+\alpha)}{4x^2} \right) g^{\pm},$$

where $g_0^{\pm}, g_1^{\pm} \in \mathbb{C}$ are the existing and finite limits

$$\begin{aligned} g_0^{\pm} &= \lim_{x \to 0^{\pm}} |x|^{\frac{\alpha}{2}} g(x) \\ g_1^{\pm} &= \lim_{x \to 0^{\pm}} |x|^{-(1+\frac{\alpha}{2})} (g(x) - g_0^{\pm} |x|^{-\frac{\alpha}{2}}). \end{aligned} \qquad (39)$$

In Theorem 6 the existence and finiteness of the limits (39) is guaranteed by the distributional constraint $(-\frac{d^2}{dx^2} + k^2|x|^{2\alpha} + \frac{\alpha(2+\alpha)}{4x^2})g^{\pm} \in L^2(\mathbb{R}^{\pm}, \mathrm{d}x)$. A completely analogous unitary equivalence and fibred decomposition like (36)-(37) holds for all other self-adjoint realisations of the Laplace-Beltrami operator on Grushin cylinder, as classified in Theorem 3 [18].

Each $A_{\alpha}(k)$ is the k-th transversal momentum mode of the operator $H_{\alpha}^{B}$ on cylinder, in the sense of the isomorphism (32), namely with respect to the momentum conjugate to the y-variable. By compactness, these are discrete modes and, as seen from (38), the boundary condition at $x = 0$ has the same form ($g_0^- = g_0^+, g_1^- = -g_1^+$) in each mode, and moreover it does not couple distinct modes. Because of this structure, the bridging operator $H_{\alpha}^{B}$ is said to be 'uniformly fibred', and in fact all other extensions classified in Theorem 3 are uniformly fibred too [18]. Uniformly fibred extensions generate a heat or Schrödinger flow that is reduced into the discrete modes $k$.

A careful spectral analysis [15] shows that for each (uniformly fibred) extension from Theorem 3, the transversal momentum modes are energetically increasingly ordered in the sense of increasing $|k|$, meaning in particular that the zero-th mode is the lowest energy one, and that for the bridging operator all modes have only non-negative, essential, absolutely continuous spectrum.

Comparing (38) with (4) one recognises that $A_{\alpha}(0) = A_{\alpha}^{B}$. This and the considerations made in Section 1 finally show that the heat flow generated by the bridging operator $H_{\alpha}^{B}$ starting with a function $f_{\text{iniz}}$ on the cylinder which belongs to the zero-th transversal momentum mode and therefore is constant in y, say, $f_{\text{iniz}}(x, y) = \varphi(x)$, produces at times $t > 0$ and evoluted function

$$f(t; x, y) = u(t, x) \qquad (40)$$

(still belonging to the zero mode) where $u$ solves the one-dimensional initial value problem (1) with initial datum $\varphi$.
### 3 Related settings: Grushin planes and almost Riemannian manifolds

The subject of geometric quantum confinement away from the metric's singularity, and transmission across it, for quantum particles or for the heat flow on degenerate Riemannian manifolds is experiencing a fast growth in the recent years. Such themes are particularly active with reference to Grushin structures on cylinder, cone, and plane [5, 8, 7, 17, 20, 6, 4, 3], as well as, more generally, on two-dimensional orientable compact almost-Riemannian manifolds of step two [5], d-dimensional regular almost-Riemannian and sub-Riemannian manifolds [22, 14].

Of significant relevance is the counterpart model to the Grushin-type cylinder, but in the lack of compact variable. This leads to related almost Riemannian structures called 'Grushin-type planes'. In complete analogy to Section 2, these are Riemannian manifolds $M_\alpha \equiv (M, g_\alpha)$, for some $\alpha \in \mathbb{R}$, where now

$$ M^{\pm} := \mathbb{R}_{x}^{\pm} \times \mathbb{R}_{y}, \quad \mathcal{Z} := \{0\} \times \mathbb{R}_{y}, \quad M := M^{+} \cup M^{-} \quad (41) $$

and again with degenerate Riemannian metric

$$ g_{\alpha} := dx \otimes dx + |x|^{-2\alpha} dy \otimes dy. \quad (42) $$

The standard 'Grushin plane' corresponds to $\alpha = 1$. Also for a Grushin-type plane one builds the Hilbert space $\mathcal{H}_\alpha$, defined as in (11), and the Laplace-Beltrami differential operator $\Delta_{\mu_\alpha} := \text{div}_{\mu_\alpha} \circ \nabla$, explicitly given again by the analogue of (13), and minimally realised as the analogue of (14) on smooth functions compactly supported within each open half-plane. This yields the densely defined, non-negative, symmetric operator $H_\alpha$, and poses the problem of self-adjointness of $H_\alpha$, in order to analyse the generated heat or Schrödinger flow.

**Theorem 7** ([14, 21, 17, 20]). (i) If $\alpha \in [-1, 1)$, then $H_\alpha$ is not essentially self-adjoint in $\mathcal{H}_\alpha$ and has infinite deficiency index.

(ii) If $\alpha \in (-\infty, -1) \cup [1, +\infty)$, then $H_\alpha$ is essentially self-adjoint and therefore the Grushin-type plane $M_\alpha$ induces geometric quantum confinement.

The above regime of essential self-adjointness was implicitly established in [14] as an adaptation of the previous perturbative analysis [22] devised for the compactified version of the manifold; the complete identification of essential self-adjointness and lack thereof was subsequently obtained in [21, 17, 20] within a non-perturbative, novel scheme of constant-fibre direct integral decomposition of the Hilbert space $\mathcal{H}_\alpha = L^2(M, d\mu_\alpha)$ that generalises the direct integral decomposition (32)-(35) one performs in the compact case. This replaces uniformly fibred extensions on cylinder of the form (37) discussed above, namely,

$$ \bigoplus_{k \in \mathbb{Z}} A_{\alpha}(k), $$

$A_\alpha(k)$ acting self-adjointly on the fibre Hilbert space $\mathfrak{h} = L^2(\mathbb{R})$, with uniformly fibred direct integral extensions

$$ \int_{\mathbb{R}} A_{\alpha}(\xi) d\xi, $$

where the fibre operator $A_\alpha(\xi)$ on $\mathfrak{h}$ now depends on the continuous Fourier mode $\xi$, dual to the non-compact variable $y$.

It is worth observing that the regime of self-adjointness for $\alpha$-Grushin cylinders and planes differ when $\alpha \in (-3, -1)$ (compare Theorems 2 and 7). This is due to the different nature of the direct sum and direct integral decompositions: indeed, when $\alpha \in (-3, -1)$, the only Fourier mode that is not self-adjoint is the
zero-th one, which brings a non-trivial contribution to the sum, but not to the integral. As a consequence, when $\alpha \in (-3, -1)$ the zero mode of a generic function $\psi \equiv \psi(x, y)$ hitting $\mathcal{Z}$ in the cylinder, namely the average on $\mathbb{S}^1$

$$ \psi_0(x) = \frac{1}{\sqrt{2\pi}} \int_{\mathbb{S}^1} \psi(x, y) \, dy, $$

does cross the singularity, whereas the zero mode of $\psi$ on the plane, namely

$$ \psi_0(x) = \frac{1}{\sqrt{2\pi}} \int_{\mathbb{R}} \psi(x, y) \, dy, $$

does not cross the singularity. The case $\alpha = -1$ is different as well between cylinder and plane: indeed, the non-self-adjoint Fourier modes are $\xi \in (-1, 1)$ for the plane, and $k = 0$ for the cylinder, thus yielding deficiency index of $H_\alpha$ equal to infinity for the plane, and equal to 2 for the cylinder.

As a matter of fact, the lack of compactness makes the systematic identification of non-trivial self-adjoint extensions of $H_\alpha$ considerably harder and so far no explicit classification is available that mirrors Theorem 3 for the plane.

Beside the above concrete cylindrical and planar settings, the deep connection between geometry and self-adjointness is investigated for the problem of geometric confinement on more general almost-Riemannian structures. This includes ‘two-step two-dimensional almost-Riemannian structures’, characterised by an orthonormal frame for the metric in the vicinity of the singularity locus $\mathcal{Z}$ of the form [2]

$$ X_1(x, y) = \frac{\partial}{\partial x}, \quad X_2(x, y) = x e^{\phi(x, y)} \frac{\partial}{\partial y} \qquad (43) $$

(to be compared to (8) with $\alpha = 1$). The essential self-adjointness of the corresponding minimally defined Laplace-Beltrami in the case of compactified $\mathcal{Z}$ was established in [5].

From a related perspective, the already observed circumstance that Grushin-type cylinders or planes are, classically, geodesically incomplete, but can induce, quantum-mechanically, geometric confinement (a condition that occurs more generally for regular almost-Riemannian manifold with compact singular set), poses an intriguing problem as far as semi-classical analysis is concerned. Indeed, reinstating Planck's constant in the Schrödinger equation

$$ i\partial_t\psi + \varepsilon^2\Delta_{\mu\alpha}\psi = 0, \quad \varepsilon > 0 \qquad (44) $$

(in the regime of $\alpha$ in which the minimally defined $\Delta_{\mu\alpha}$ is unambiguously realised self-adjointly), semi-classics show, informally speaking, that as $\varepsilon \downarrow 0$ solutions get concentrated and evolves around geodesics. Therefore, the above-mentioned classical/quantum discrepancy makes the semi-classical analysis necessarily brake down in the limit.

Such a discordance between classical and quantum picture can be at least partially resolved by appealing to different quantisation procedures on the considered Riemannian manifold, in practice considering corrections of the Laplace-Beltrami operator that have a suitable interpretation of free kinetic energy, much in the original spirit of [12]. Most of coordinate-invariant quantisation procedures (including path integral quantisation, covariant Weyl quantisation, geometric quantisation, and finite-dimensional approximation to Wiener measures) modify $\Delta_{\mu\alpha}$ with a term that depends on the scalar curvature $R_\alpha$ (which, in two dimensions, is twice the Gaussian curvature $K_\alpha$). This produces a replacement in (44), in two dimensions, of $-\Delta_{\mu\alpha}$ with the ‘curvature Laplacian’

$$ - \Delta_{\mu\alpha} + cK_{\alpha} \qquad (45) $$

for suitable $c \ge 0$. In the recent work [4] it was indeed shown, for generic two-step two-dimensional almost-Riemannian manifolds with compact singular set, that irrespective of $c \in (0, \frac{1}{2})$ the above correction
washes essential self-adjointness out, yielding a quantum picture where the Schrödinger particle does reach the singularity much as the classical particle does. (At the expenses of some further technicalities, the whole regime $c > 0$ can be covered as well.) For concreteness, in the Grushin cylinder the effect of the curvature correction is evidently understood as a compensation between $K = -\frac{2}{x^2}$ (see (9) above) and the singular term $\frac{3}{4x^2}$ of the (unitary equivalent) Laplace-Beltrami operator. Still, the classical/quantum discrepancy discussed so far remains unexplained in more general settings.

Concerning, instead, the heat flow, a satisfactory interpretation of the heat-confinement in the Grushin cylinder is known in terms of Brownian motions [6] and random walks [1]: roughly speaking, random particles are lost in the infinite area accumulated along $\mathcal{Z}$: the latter, in practice, acts as a barrier. Clearly, whereas curvature Laplacians are meaningful in the above context of inducing a non-confining (transmitting) Schrödinger flow on two-step two-dimensional almost-Riemannian manifolds (including the Grushin cylinders), thus making quantum and classical picture more alike and well connected by semi-classics, this has no direct meaning instead in application to the heat flow on Riemannian or almost Riemannian manifolds. Indeed, as long as one regards the heat equation on manifold as a limit of a space-time discretised random walk, the stochastic process' generator is the Laplace-Beltrami operator.

Generalisations of [5] have been established in [14, 22] from two-step two-dimensional almost-Riemannian structures to any dimensions, any step, and even to sub-Riemannian geometries, provided that certain geometrical assumptions on the singular set are taken. The main difficulty is the treatment of the 'tangency' (or 'characteristic') points: these are points belonging to the singularity of the metric structure where the vector distribution is tangent to the singularity. They are never present in Grushin cylinder or two-step almost-Riemannian structures, but may appear for example in three-step structures, such as, for instance,

$$X_1(x, y) = \frac{\partial}{\partial x}, \quad X_2(x, y) = (y - x^2)\frac{\partial}{\partial y}, \quad (x, y) \in \mathbb{R}^2, \quad (46)$$

where the singularity is the parabola $y = x^2$ and the origin $(0, 0)$ is a tangency point. Virtually nothing is known on the heat or the quantum confinement on such singular structures, including the simplest example (46) (see [13] for further remarks). First preliminary results in this respect were recently obtained in [3], where the interpretation of almost-Riemannian structures as special Lie manifolds permits to study some closure properties of singular perturbations of the Laplace-Beltrami operator even in the presence of tangency points. This opens new perspectives of treating several types of different singularities in sub-Riemannian geometry within the same unifying theory.

## 4 A numerical glance at the bridging heat evolution

In this final Section we present and comment on qualitative features of the solution to the one-dimensional problem (1), obtained by numerical integration, also in comparison with the initial value problem for the classical heat equation on $\mathbb{R}$.

As already argued, it is the determination of the (integral kernel of) the heat propagator $\exp(-tA_\alpha^B)$, $t > 0$, to be hard analytically, and this is due to the presence of boundary conditions for the solution at $x = 0$ and any positive time.

Numerics then represent a first, valuable way to access relevant aspects of the transmission of the heat flow between positive and negative half-line with bridging boundary conditions, and one may envisage that a systematic comparison will be launched numerically between analogous heat flows with different transmission protocols among those surveyed in Section 2. Ours, here, is only an initial numerical glance at the bridging heat evolution to provide some insight and anticipate future investigations.

Our numerical approach consists in approximating the solution $u = e^{-tA_\alpha^B}\varphi$ to the problem (1) by means of an approximated version of both the spatial convolution integral between propagator's kernel
and $\varphi$, and the complex line integral that turns the resolvent of $A_{\alpha}^{B}$ into its semi-group.
More precisely, let us write

$$u(t, x) = (e^{-tA_{\alpha}^{B}} \varphi)(x) = \int_{\mathbb{R}} \mathcal{K}_{\alpha}^{B}(t; x, y) \varphi(x) \, dx \quad (47)$$

where $\mathcal{K}_{\alpha}^{B}(\cdot, \cdot)$ is the integral kernel of the bridging heat propagator. In turn, let us exploit the relation

$$\begin{aligned} e^{-tA_{\alpha}^{B}} &= \mathcal{L}^{-1}((A_{\alpha}^{B} - (\cdot)1)^{-1})(t) \\ &= \frac{1}{2\pi i} \int_{\Gamma} e^{-zt}((A_{\alpha}^{B} - z1)^{-1}) \, dz, \quad t > 0, \end{aligned} \quad (48)$$

as an identity between bounded operators on $L^2(\mathbb{R})$ and with the integral understood in the Riemann sense in the strong operator topology, $\Gamma$ being a straight line in $\mathbb{C}$ orthogonal to the real axis in the open left half-plane, and $\mathcal{L}^{-1}$ denoting the inverse Laplace transform (the non-negativity of $A_{\alpha}^{B}$ has led here to the choice $\Re z < 0$). (48) connects the resolvent of $A_{\alpha}^{B}$ at the complex point $z$ with the semi-group at time $t > 0$, and in terms of the integral kernels $(A_{\alpha}^{B} - z1)^{-1}(x, y)$ of the resolvent and $\mathcal{K}_{\alpha}^{B}(t; x, y)$ of the propagator it reads

$$\mathcal{K}_{\alpha}^{B}(t; x, y) = \frac{1}{2\pi i} \int_{-1-i\cdot\infty}^{-1+i\cdot\infty} e^{-zt}((A_{\alpha}^{B} - z1)^{-1})(x, y) \, dz, \quad t > 0. \quad (49)$$

The combinations of (47) and (49) produces the solution $u$ and the two integrations contained therein may be computed numerically with standard packages.

Of course, for (47) and (49) to be implementable one needs to know the (integral kernel of) the resolvent $(A_{\alpha}^{B} - z1)^{-1}$. This is not so hard knowledge to achieve from the underlying structure (4) of the operator $A_{\alpha}^{B}$, once it is interpreted as a self-adjoint extension of the differential operator (2) minimally defined on smooth functions compactly supported on $\mathbb{R}$ away from the origin. For this status of extension operator, one can appeal to the general Krečn-Višik-Birman theory of self-adjoint extensions of lower semi-bounded and densely defined symmetric operators on Hilbert space [16], and obtain $(A_{\alpha}^{B} - z1)^{-1}$ fairly explicitly.

The net result of this computation gives the following expression for the integral kernel of $(A_{\alpha}^{B} - z1)^{-1}$. With respect to the canonical decomposition

$$L^2(\mathbb{R}, dx) \stackrel{\sim}{\longrightarrow} L^2(\mathbb{R}^+, dx) \oplus L^2(\mathbb{R}^-, dx), \quad u \longmapsto \begin{pmatrix} u^+ \\ u^- \end{pmatrix} \quad (50)$$

(that is, $u^{\pm}(x) := u(x)$ for $x \ge 0$), consider the unitary transformation

$$\begin{gathered} U : L^2(\mathbb{R}^+, dx) \oplus L^2(\mathbb{R}^-, dx) \stackrel{\sim}{\longrightarrow} L^2(\mathbb{R}^+, dx) \oplus L^2(\mathbb{R}^-, dx) \\ U \begin{pmatrix} u^+ \\ u^- \end{pmatrix}(x) = \begin{pmatrix} u^+(x) \\ u^-(x) \end{pmatrix}, \quad x > 0, \end{gathered} \quad (51)$$

and set $\mathcal{R}_{\alpha}^{B}(z) := U(A_{\alpha}^{B} - z1)^{-1}U^{-1}$. Then

$$(A_{\alpha}^{B} - z1)^{-1} = U^{-1} \mathcal{R}_{\alpha}^{B}(z) U \quad (52)$$

and the integral kernel of $\mathcal{R}_{\alpha}^{B}(z)$ is given by

$$\mathcal{R}_{\alpha}^{B}(z)(x, y) = \mathcal{G}_{\alpha, z}(x, y) \begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix} - \frac{i\pi}{8} \cos\left(\frac{\pi\alpha}{2}\right) e^{i\frac{\pi\alpha}{2}} \begin{pmatrix} 1 & 1 \\ 1 & 1 \end{pmatrix} P_{\alpha, z}(x) P_{\alpha, z}(y), \quad (53)$$
$$x > 0, y > 0,$$
<figure><img src="image_3.png" /><figcaption>Figure 2: Solution $u(t,x)$ to the heat-bridging initial value problem (1) with Gaussian initial datum $\varphi(x) = e^{-(x-2)^2}$ (red curve). Plot of $|u(t, \cdot)|$ at $t = 0.5$ (magenta dotted line) and $t = 2$ (blue dotted line).</figcaption></figure>

where, in terms of the Bessel functions of first and second kind $J_\nu$ and $Y_\nu$,

$$ 
\begin{aligned}
P_{\alpha,z}(x) &= \sqrt{x} J_{\frac{1+\alpha}{2}}(x\sqrt{z}) + i\sqrt{x} Y_{\frac{1+\alpha}{2}}(x\sqrt{z}) & (\Im\sqrt{z} > 0), \\
Q_{\alpha,z}(x) &= 2\sqrt{x} J_{\frac{1+\alpha}{2}}(x\sqrt{z}) &
\end{aligned}
\qquad (54) 
$$

and

$$ \mathcal{G}_{\alpha,z}(x,y) = -\frac{i\pi}{4} \begin{cases} P_{\alpha,z}(x) Q_{\alpha,z}(y), & \text{if } 0 < y < x, \\ Q_{\alpha,z}(x) P_{\alpha,z}(y), & \text{if } 0 < x < y. \end{cases} \qquad (55) $$

When formulas (47), (49), (52), and (53) are implemented numerically we obtain a scenario exemplified in Figures 2, 3, 4, 5 below.

For concreteness, the bridging-heat evolution is considered, namely the solution $u \equiv u(t,x)$ to (1), of an initial datum $\varphi$ essentially supported on the right half-line. An initial Gaussian is seen to evolve at later times with the typical heat-flow flattening of the solution, with the immediate formation of the characteristic bridging behaviour at $x = 0$ (Fig. 2).

Notably, if $\varphi$ is additionally shot with an initial non-zero momentum towards the singularity, its evolution displays an oscillation given by the superposition of an in-going wave and a component that bounces backwards (Fig. 3), as compared with the evolution at the same time of the same Gaussian with no initial momentum.

It is also pretty transparent that the bridging-heat flow has a regularising effect at every $t > 0$, as observed with the evolution of an initial step function. (Fig. 4).

We have also further evidence of a qualitatively similar behaviour of the free heat flow and the bridging-heat flow, but of course for the characteristic boundary condition of bridging type at the origin (Fig. 5).

Whereas, as said, this provides only a first glance at the qualitative properties of the bridging-heat evolution on two connected half-lines, the evidences collected here are encouraging and further corroborate the quest for the analytic identification of counterpart $L^p-L^q$ estimates, smoothing estimates, and space-time (Strichartz) estimates for the bridging-heat flow, as compared to (5) for the classical heat flow.
<figure><img src="image_4.png" /><figcaption>Figure 3: Comparison at time $t = 1.5$ between the solution to the heat-bridging initial value problem (1) with zero-momentum Gaussian initial datum $\varphi(x) = e^{-(x-2)^2}$ (magenta dotted curve) and with non-zero momentum Gaussian initial datum $\varphi_2(x) = e^{-3ix}e^{-(x-2)^2}$ towards left. Both plots are of $|u(t, \cdot)|$. The evolution of the Gaussian initially shot towards left displays in-going + outgoing oscillation.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>Figure 4: Solution $u(t, x)$ to the heat-bridging initial value problem (1) with initial datum $\varphi(x)$ given by the characteristic function of the interval $[\frac{1}{2}, \frac{5}{2}]$ (red curve). Plot of $|u(t, \cdot)|$ at $t = 0.5$ (blue dotted line)</figcaption></figure>

## Acknowledgement

This work is partially supported by the QUACO project (grant ANR-17-CE-40-0007-01), the EIPHI Graduate School (ANR-17-EURE-0002), the CONSTAT project (funded by the Conseil Régional de Bourgogne Franche Comté and the European Union through the PO FEDER Bourgogne 2014/2020 programmes), the MIUR-PRIN 2017 project MaQuMA cod. 2017ASFLJR, the INdAM – Italian National Institute for Higher Mathematics, and the Alexander von Humboldt Foundation.
<figure><img src="image_6.png" /><figcaption>Figure 5: Comparison at time $t = 1.5$ between the solution to the heat-bridging initial value problem (1) with Gaussian initial datum $\varphi(x) = e^{-(x-2)^2}$ (blue dotted curve) and the solution to the ordinary heat equation on $\mathbb{R}$ (green curve).</figcaption></figure>

References

[1] A. AGRACHEV, U. BOSCAIN, R. NEEL, AND L. RIZZI, *Intrinsic random walks in Riemannian and sub-Riemannian geometry via volume sampling*, ESAIM Control Optim. Calc. Var., 24 (2018), pp. 1075–1105.

[2] A. AGRACHEV, U. BOSCAIN, AND M. SIGALOTTI, *A Gauss-Bonnet-like formula on two-dimensional almost-Riemannian manifolds*, Discrete Contin. Dyn. Syst., 20 (2008), pp. 801–822.

[3] I. BESCHASTNYI, *Closure of the Laplace-Beltrami operator on 2D almost-Riemannian manifolds and semi-Fredholm properties of differential operators on Lie manifolds*, arXiv:2104.07745 (2021).

[4] I. BESCHASTNYI, U. BOSCAIN, AND E. POZZOLI, *Quantum Confinement for the Curvature Laplacian* –$\Delta + cK$ on 2D-Almost-Riemannian Manifolds, Potential Analysis, (2021).

[5] U. BOSCAIN AND C. LAURENT, *The Laplace-Beltrami operator in almost-Riemannian geometry*, Ann. Inst. Fourier (Grenoble), 63 (2013), pp. 1739–1770.

[6] U. BOSCAIN AND R. W. NEEL, *Extensions of Brownian motion to a family of Grushin-type singularities*, Electron. Commun. Probab., 25 (2020), pp. Paper No. 29, 12.

[7] U. BOSCAIN AND D. PRANDI, *Self-adjoint extensions and stochastic completeness of the Laplace-Beltrami operator on conic and anticonic surfaces*, J. Differential Equations, 260 (2016), pp. 3234–3269.

[8] U. BOSCAIN, D. PRANDI, AND M. SERI, *Spectral analysis and the Aharonov-Bohm effect on certain almost-Riemannian manifolds*, Comm. Partial Differential Equations, 41 (2016), pp. 32–50.

[9] N. BURQ, F. PLANCHON, J. G. STALKER, AND A. S. TAHVILDAR-ZADEH, *Strichartz estimates for the wave and Schrödinger equations with the inverse-square potential*, J. Funct. Anal., 203 (2003), pp. 519–549.
[10] O. CALIN AND D.-C. CHANG, *Sub-Riemannian geometry*, vol. 126 of Encyclopedia of Mathematics and its Applications, Cambridge University Press, Cambridge, 2009. General theory and examples.

[11] J. DEREZIŃSKI AND V. GEORGESCU, *On the domains of Bessel operators*, Ann. Henri Poincaré, 22 (2021), pp. 3291–3309.

[12] B. S. DEWITT, *Dynamical Theory in Curved Spaces. I. A Review of the Classical and Quantum Action Principles*, Rev. Mod. Phys., 29 (1957), pp. 377–397.

[13] V. FRANCESCHI, D. PRANDI, AND L. RIZZI, *Recent results on the essential self-adjointness of sub-Laplacians, with some remarks on the presence of characteristic points*, Séminaire de théorie spectrale et géométrie, 33 (2015-2016), pp. 1–15.

[14] ——, *On the essential self-adjointness of singular sub-Laplacians*, Potential Anal., 53 (2020), pp. 89–112.

[15] M. GALLONE AND A. MICHELANGELI, *Quantum particle across Grushin singularity*, Journal of Physics A: Mathematical and Theoretical, (2021).

[16] M. GALLONE, A. MICHELANGELI, AND A. OTTOLINI, *Kreïn-Višik-Birman self-adjoint extension theory revisited*, in Mathematical Challenges of Zero Range Physics, A. Michelangeli, ed., INdAM-Springer series, Vol. 42, Springer International Publishing, 2020, pp. 239–304.

[17] M. GALLONE, A. MICHELANGELI, AND E. POZZOLI, *On geometric quantum confinement in Grushin-type manifolds*, Z. Angew. Math. Phys., 70 (2019), pp. Art. 158, 17.

[18] ——, *Geometric confinement and dynamical transmission of a quantum particle in Grushin cylinder*, arXiv:2003.07128 (2020).

[19] A. MICHELANGELI, *Global well-posedness of the magnetic Hartree equation with non-Strichartz external fields*, Nonlinearity, 28 (2015), pp. 2743–2765.

[20] E. POZZOLI, *Quantum Confinement in $\alpha$-Grushin planes*, in Mathematical Challenges of Zero-Range Physics, A. Michelangeli, ed., Springer INdAM Series, Springer International Publishing, 2021, pp. 229–237.

[21] ——, *Models of quantum confinement and perturbative methods for point interactions*, Master Thesis (2018).

[22] D. PRANDI, L. RIZZI, AND M. SERI, *Quantum confinement on non-complete Riemannian manifolds*, J. Spectr. Theory, 8 (2018), pp. 1221–1280.

[23] M. REED AND B. SIMON, *Methods of modern mathematical physics. II. Fourier analysis, self-adjointness*, Academic Press [Harcourt Brace Jovanovich, Publishers], New York-London, 1975.

[24] K. SCHMÜDGEN, *Unbounded self-adjoint operators on Hilbert space*, vol. 265 of Graduate Texts in Mathematics, Springer, Dordrecht, 2012.

[25] B. WANG, Z. HUO, C. HAO, AND Z. GUO, *Harmonic analysis method for nonlinear evolution equations. I*, World Scientific Publishing Co. Pte. Ltd., Hackensack, NJ, 2011.