\documentclass[CRMATH,Unicode,biblatex,XML]{cedram}

\TopicFR{\'Equations aux dérivées partielles}
\TopicEN{Partial differential equations}

\addbibresource{CRMATH_Colin-de-Verdire_20240031.bib}

\usepackage{mathtools}
\usepackage{xcolor}

\newcommand{\R}{{\mathbb{R}}}
\newcommand{\N}{{\mathbb{N}}}
\newcommand{\ha}{\frac{1}{2}}
\newcommand{\ra}{\to}
\newcommand{\ga}{\alpha}
\newcommand{\gb}{\beta}
\newcommand{\gl}{\lambda}
\newcommand\OPD{pseudo-differential operator}

\newcommand{\dd}{\mathrm{d}}
\newcommand{\irm}{\mathrm{i}}
\DeclareMathOperator{\expe}{e}


                                   
%\def\YCV{Y. Colin de Verdi{\`e}re}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\makeatletter
\let\save@mathaccent\mathaccent
\newcommand*\if@single[3]{%
 \setbox0\hbox{${\mathaccent"0362{#1}}^H$}%
 \setbox2\hbox{${\mathaccent"0362{\kern0pt#1}}^H$}%
 \ifdim\ht0=\ht2 #3\else #2\fi
 }
%The bar will be moved to the right by a half of \macc@kerna, which is computed by amsmath:
\newcommand*\rel@kern[1]{\kern#1\dimexpr\macc@kerna}
%If there's a superscript following the bar, then no negative kern may follow the bar;
%an additional {} makes sure that the superscript is high enough in this case:
\newcommand*\widebar[1]{\@ifnextchar^{{\wide@bar{#1}{0}}}{\wide@bar{#1}{1}}}
%Use a separate algorithm for single symbols:
\newcommand*\wide@bar[2]{\if@single{#1}{\wide@bar@{#1}{#2}{1}}{\wide@bar@{#1}{#2}{2}}}
\newcommand*\wide@bar@[3]{%
 \begingroup
 \def\mathaccent##1##2{%
%Enable nesting of accents:
 \let\mathaccent\save@mathaccent
%If there's more than a single symbol, use the first character instead (see below):
 \if#32 \let\macc@nucleus\first@char \fi
%Determine the italic correction:
 \setbox\z@\hbox{$\macc@style{\macc@nucleus}_{}$}%
 \setbox\tw@\hbox{$\macc@style{\macc@nucleus}{}_{}$}%
 \dimen@\wd\tw@
 \advance\dimen@-\wd\z@
%Now \dimen@ is the italic correction of the symbol.
 \divide\dimen@ 3
 \@tempdima\wd\tw@
 \advance\@tempdima-\scriptspace
%Now \@tempdima is the width of the symbol.
 \divide\@tempdima 10
 \advance\dimen@-\@tempdima
%Now \dimen@ = (italic correction / 3) - (Breite / 10)
 \ifdim\dimen@>\z@ \dimen@0pt\fi
%The bar will be shortened in the case \dimen@<0 !
 \rel@kern{0.6}\kern-\dimen@
 \if#31
 \overline{\rel@kern{-0.6}\kern\dimen@\macc@nucleus\rel@kern{0.4}\kern\dimen@}%
 \advance\dimen@0.4\dimexpr\macc@kerna
%Place the combined final kern (-\dimen@) if it is >0 or if a superscript follows:
 \let\final@kern#2%
 \ifdim\dimen@<\z@ \let\final@kern1\fi
 \if\final@kern1 \kern-\dimen@\fi
 \else
 \overline{\rel@kern{-0.6}\kern\dimen@#1}%
 \fi
 }%
 \macc@depth\@ne
 \let\math@bgroup\@empty \let\math@egroup\macc@set@skewchar
 \mathsurround\z@ \frozen@everymath{\mathgroup\macc@group\relax}%
 \macc@set@skewchar\relax
 \let\mathaccentV\macc@nested@a
%The following initialises \macc@kerna and calls \mathaccent:
 \if#31
 \macc@nested@a\relax111{#1}%
 \else
%If the argument consists of more than one symbol, and if the first token is
%a letter, use that letter for the computations:
 \def\gobble@till@marker##1\endmarker{}%
 \futurelet\first@char\gobble@till@marker#1\endmarker
 \ifcat\noexpand\first@char A\else
 \def\first@char{}%
 \fi
 \macc@nested@a\relax111{\first@char}%
 \fi
 \endgroup
}
\makeatother

\let\oldbar\bar
\renewcommand*{\bar}[1]{\mathchoice{\widebar{#1}}{\widebar{#1}}{\widebar{#1}}{\oldbar{#1}}}

\let\oldtilde\tilde
\renewcommand*{\tilde}[1]{\mathchoice{\widetilde{#1}}{\widetilde{#1}}{\oldtilde{#1}}{\oldtilde{#1}}}
%\let\tilde\widetilde

%\let\hat\widehat
\let\oldhat\hat
\renewcommand*{\hat}[1]{\mathchoice{\widehat{#1}}{\widehat{#1}}{\oldhat{#1}}{\oldhat{#1}}}
%\let\tilde\widetilde

\renewcommand*{\to}{\mathchoice{\longrightarrow}{\rightarrow}{\rightarrow}{\rightarrow}}
\let\oldmapsto\mapsto
\renewcommand*{\mapsto}{\mathchoice{\longmapsto}{\oldmapsto}{\oldmapsto}{\oldmapsto}}

\newcommand\blfootnote[1]{%
  \begingroup
  \renewcommand\thefootnote{}\footnote{#1}%
  \addtocounter{footnote}{-1}%
  \endgroup
}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\newcommand*{\mk}{\mkern -1mu}
\newcommand*{\Mk}{\mkern -2mu}
\newcommand*{\mK}{\mkern 1mu}
\newcommand*{\MK}{\mkern 2mu}

\hypersetup{urlcolor=purple, linkcolor=blue, citecolor=red}


\newcommand*{\romanenumi}{\renewcommand*{\theenumi}{\roman{enumi}}}
\newcommand*{\Romanenumi}{\renewcommand*{\theenumi}{\Roman{enumi}}}
\newcommand*{\alphenumi}{\renewcommand*{\theenumi}{\alph{enumi}}}
\newcommand*{\Alphenumi}{\renewcommand*{\theenumi}{\Alph{enumi}}}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
																		
\def\dagger{**}

\author[L. Boutet de Monvel]{\firstname{Louis} \lastname{Boutet de Monvel}\dead} 

\author[Y. Colin de Verdière]{\firstname{Yves} \lastname{Colin de Verdière}
\IsCorresp}

\address{Institut Fourier,  Unit{\'e} mixte  de recherche CNRS-UJF 5582,  BP 74, 38402-Saint Martin d'H\`eres Cedex, France}

\email{yves.colin-de-verdiere@univ-grenoble-alpes.fr}


\dedicatory{Dedicated to the memory of Louis Boutet de Monvel.}


\keywords{\kwd{Pseudo-differential operators} \kwd{Dirichlet-to-Neumann map} \kwd{Poisson operator}}


\altkeywords{\kwd{Opérateurs pseudos-différentiels} \kwd{application Dirichlet-to-Neumann} \kwd{Opérateurs de Poisson}}


\subjclass{35S05, 35S15}

\title{The action of pseudo-differential operators on functions harmonic outside a smooth hyper-surface}

\alttitle{L'action d'opérateurs pseudo-différentiels sur les fonctions harmoniques en dehors d'une hypersurface lisse}

\begin{abstract}
The goal of this note is to describe the action of  pseudo-differential operators on the space of square integrable  functions which are harmonic outside a smooth closed hyper-surface  of a compact Riemannian manifold.
\end{abstract}

\begin{altabstract}
Dans cette note, nous étudions un opérateur du type
$B=\mathcal{P} ^* A \mathcal{P}$ où $A$ est un opérateur pseudo-différentiel et $\mathcal{P}$ l'opérateur de Poisson bilatéral associé à une hypersurface. Nous montrons que, sous certaines conditions, B est un opérateur pseudo-différentiel sur cette hypersurface dont nous calculons le symbole principal.
\end{altabstract}


\begin{document}




\maketitle




The goal of this note is to describe the action of \OPD s
 on the space $\mathcal{H}$  of $L^2$ functions 
which are harmonic  outside a smooth closed 
hyper-surface $Z$ of a compact  Riemannian
manifold without boundary  $(X,g)$ and whose traces from  both sides of $Z$
co\"incide. We will
 represent these $L^2$  harmonic functions as harmonic extensions
 of functions  in the Sobolev space  $H^{-1/2}(Z) $ by
 a Poisson operator $\mathcal{P}$.
The main result says that,  if $A$ is a \OPD ~of degree $d<3$ on $X$,  the operator
\[
 B=\mathcal{P}^\star  \circ A \circ \mathcal{P} ~
\]
is a \OPD~ on $Z$  of  degree $d-1 $  whose principal symbol of degree
 $d-1$  can be computed by integration of the principal symbol of $A$
 on
 the co-normal bundle of $Z$.

 These ``bilateral'' extensions are  simpler (at least
 for the Laplace operator) 
 than the ``unilateral'' ones whose study is the theory of
 \OPD s on manifolds with boundary (see \cite {Cal,LBM1_1, LBM1_2,LBM2,Grub,Taylor}).

\section{Symbols}
The following classes of symbols are defined in the books
\cite[Section~7.1]{Grub} and in \cite[Section~18.1]{Ho}. A \emph{symbol of degree}
$d$ on $U_x\times \R^n_\xi$ where $U$ is an open set in $\R^N$ is a
 smooth
 complex valued function $a(x,\xi)$ on $U\times \R^n$ which satisfies
 the following estimates:
for any multi-indices $(\ga,\gb)$, there exists a constant
 $C_{\ga,\gb}$ so that \[ | D^\ga _x D_\xi ^\gb a (x;\xi)|
 \leq C_{\ga,\gb} (1+\| \xi \|)^{d-|\gb|}~.\]
The symbol $a$ is called \emph{classical} if $a $ admits
 an expansion $a \sim \sum_{l=0}^\infty  a_{d-l}$ where  $a_j$ 
 is homogeneous of degree $j$ ($j$ an integer) 
for  $\xi\in \R^n$ large enough; more precisely,
 for any $J\in \N$,  $ a-\sum_{j=0}^J a_{d-j} $  is a symbol of degree $d-J-1$.

We will need the 
\begin{lemm}\label{lemm:intsymb} 
If $a (x;\xi,\eta )$ is a  symbol of degree $d<-1$ defined
 on $U_x\times \left( \R^{n}_{\xi}\times \R_\eta \right)$,
  $ b(x;\xi )=\int _\R a(x;\xi,\eta ) \,\dd \eta$ is  a  symbol
 of degree $d+1$ defined on $U_x\times \R^{n}_\xi$. Moreover, 
 if $a$ is classical, $b$ is also classical and the homogeneous
 components of $b$ 
are given for $l\leq d+1$, by
 $ b_l(x;\xi )=\int _\R a_{l-1}(x;\xi,\eta ) \dd \eta$
\end{lemm}

\section{A general reduction Theorem for \OPD s}
We  choose local coordinates in some neighborhood of a point
 in $Z$ denoted $x=(z,y)\in \R^{d-1} \times \R$,  so that  $Z=\{ y=0
 \}$.
  We denote by $(\Omega _j ,~ j=1,\ldots, N)$  a finite cover of $Z$ by
 such charts and denote by $\Omega _0$ an open set disjoint from $Z$
 so that $X=\bigcup_{j=0}^N \Omega _j $. We  choose the charts $\Omega
 _j$
 so that the densities $|\dd z|$ and $|\dd x| $ are the Lebesgue measures.

If $X$ is a smooth manifold, we denote by  $\mathcal{D'} (X)$ the space
 of generalized functions  on $X$ of which the space of
 smooth functions on $X$ 
is a dense subspace. We assume that $X$ and $Z$
 are equipped with smooth densities $|\dd x| $ and $|\dd z|$.
 This allows to identify generalized functions  with  Schwartz
 distributions, i.e. linear functionals 
  on test functions; this duality extending the $L^2$ product
 is denoted by  $\langle \;\mid \;\rangle $. 
We introduce the extension operator
 $\mathcal{E}:\mathcal{D'}  (Z) \ra \mathcal{D'}(X)$ sending the distribution
  $f $ to the distribution  $f\delta (y=0)$ defined 
\[
\langle f\delta (y=0)\mid \phi(z,y)\rangle = \langle f | \phi(z,0)
\rangle \]
and its adjoint, the trace $\mathcal{ T }:C^\infty (X)\ra C^\infty (Z)$
 defined by $\phi \ra \phi_{|Z}$.
Let  $A$ be  a \OPD ~ on $X$: let us call $A_j$ the restriction
 of $A$ to test functions compactly supported in $\Omega _j$.  We will
 work
 with one of the $A_j$'s given by the following ``quantization'' rule 
\[ 
A_ju(z,y)=\frac{1}{(2\pi)^{d}}\int_{ \R^{2d}}
  \expe^{\irm\left( \langle z-z'|\zeta\rangle  + (y-y')\eta \right)}
a_j(z,y;\zeta,\eta ) u(z',y') \,\dd z' \,\dd y' \,\dd \zeta \dd \eta ~.
\]
So we have formally, using the facts that the densities
 on $X$ and $Z$ are given by the Lebesgue measures in these local coordinates: 
\[
\mathcal{T}\circ A_j\circ\mathcal{E}v(z)= \frac{1}{(2\pi)^{d}}
\int_{\R^{2d-1}} \expe^{\irm \langle z-z'|\zeta\rangle    }
a_j(z,0;\zeta,\eta ) v(z') \,\dd z'  \,\dd \zeta \,\dd \eta ~,
\]
which we can rewrite
\[\mathcal{T}\circ A_j\circ\mathcal{E}v(z)=\frac{1}{(2\pi)^{d-1}}\int_{\R^{2d}}  
 \expe^{\irm\langle z-z'|\zeta\rangle  } b_j(z;\zeta ) v(z') \, \dd z'  \, \dd \zeta ~,
 \]
with
\begin{equation} \label{equ:symbol}
 b_j(z;\zeta ) =\frac{1}{2\pi}\int_\R  a_j(z,0;\zeta,\eta ) \, \dd \eta ~.
\end{equation}

We have the
\begin{theo}\label{theo:main}
 If $A$ is a \OPD~ on $X$ of degree $m<-1$ whose full symbol in the
 chart $\Omega _j$ is $a_j$, 
then  the operator $\mathcal{T}\circ A\circ\mathcal{E}$ is a \OPD~ on $Z$
  of degree $m+1$ whose  symbol is given in the charts $\Omega _j\cap
  Z$
 by Equation (\ref{equ:symbol}).
 \end{theo} 
 This is proved by looking at the actions on test functions compactly
 supported in the chart $\Omega _j,~j\geq 1$:  then we use Lemma \ref{lemm:intsymb}. 
 \begin{rema}
The principal symbol can be described in a more intrinsic way:
 let $z\in Z$ be given, from the  smooth densities on $T_zX$ and
 on $T_zZ$ given by $|\dd x|$ and $|\dd z|$, we get, using the Liouville
  densities, densities on the dual bundles  $T_z^\star Z$ and
  $T_z^\star X$.
 Let us denote by $\Omega^1(E)$ the
1-dimensional space of densities on the vector space $E$. From the
exact
 sequence
\[
0 \ra N^\star _z Z\ra T^\star _zX \ra T^\star _z Z \ra 0~,
\]
we deduce
\[ 
\Omega ^1 (T^\star X) \equiv \Omega ^1 (N^\star Z)\otimes 
 \Omega ^1 (T^\star Z )~
 \]
and a canonical density $\dd m (z)$  in  $\Omega ^1 (N_z^\star Z)$.
 The principal symbol of $B=\mathcal{T}\circ A\circ\mathcal{E}$ is given
 in coordinates  by $b(z,\zeta)=(1/2\pi)\int_{N_z^\star Z} a(z;\zeta,\eta )\dd m(\eta) $.
\end{rema}

\section{ The ``bilateral''  Dirichlet-to-Neumann operator}
 We will assume that the local coordinates $x=(z,y)$  along $Z$ are
 chosen
 so that $g(z,0)=  h(\dd z) +\dd y^2$ and the Riemannian volume along $Z$ is
 $|\dd x|_g=|\dd z|_h \,|\dd y|$. We will choose the associated densities on $X$
 and $Z$.
 We will denote by $\Delta _g$ the Laplace--Beltrami operator on
 $(X,g)$
 as defined by Riemannian geometers (i.e. with a minus sign in front
 of
 the second order derivatives).

 If $f $ is given on $Z$, let us denote by $\mathcal{DN}(f)$
  minus  the sum of the interior  normal derivatives on both
 sides of $Z$ of  the harmonic extension $F$  of $f$;
this always  makes sense, even if the normal bundle of $Z$ is not orientable.
 We have the
\begin{lemm} \label{lemm:harm-ext}
 The distributional Laplacian of the harmonic extension  $F$
 of a smooth function $f$ on $Z$  is 
$\Delta_g F=\mathcal{E}(\mathcal{DN}(f))$.
\end{lemm}
\begin{proof} The proof is a simple application of the Green's formula:
 by definition of the action of the Laplacian on distributions, if
 $\phi $ is a test function on $X$, $ \langle \Delta_g F\mid \phi \rangle
 \coloneqq 
 \langle  F \mid \Delta_g \phi \rangle $.
We can compute the righthandside  integral as an integral on
$X\setminus Z$
 using  Green's formula.
\[
 \int _{X\setminus Z}(F\Delta_g \phi -\phi \Delta _g F)|\dd x|_g=
 \int _Z (F\delta \phi - \phi \delta F )|\dd z|_h
 \]
where $\delta $ is   the sum of the interior normal derivatives from
 both sides of $Z$. Using the fact that $\Delta _g F =0 $ in
 $X\setminus Z$
 and $\delta \phi=0$, we get the result.
\end{proof} 

Denoting by $\Delta_g^{-1}$ the ``quasi-inverse''  of $\Delta_g$
defined
 by $\Delta_g^{-1}\phi_j =\gl_j^{-1} \phi_j$ for the eigenfunctions
 $\phi_j$ of $\Delta _g$  with non-zero eigenvalue $\gl_j$ and 
 $\Delta_g^{-1}1=0$, we have $f=\left(\mathcal{T}\circ\Delta_g
   ^{-1}\circ\mathcal{E}\right) \circ \mathcal{DN}(f)~({\rm mod~
   constants})$.
 By Theorem \ref{theo:main}, the operator $B=\mathcal{T}\circ\Delta_g
 ^{-1}
\circ\mathcal{E}$ is an elliptic self-adjoint \OPD~ on $Z$.
 The operator $\mathcal{DN}$ is a right inverse of $B$  modulo smoothing
 operators and hence also a left  inverse modulo smoothing operators.
 So that $\mathcal{DN}=B^{-1}$ is an elliptic self-adjoint of principal
 symbol the inverse of
\[
 \frac{1}{2\pi}\int_\R  (\|\zeta\|_h^2+\eta^2)^{-1} \dd \eta = 
\frac{1}{2\| \zeta \|_h} ~,
\]
namely $2\| \zeta \|_h$. Hence
\begin{theo}\label{theo:dn} 
 The bilateral Dirichlet-to-Neumann operator $\mathcal{DN}$ is   a self-adjoint
  elliptic  \OPD~  of  degree $1$ on 
 $L^2(Z,|dz|)$ and of principal symbol $2\| \zeta \|_h$.
 The kernel of  $\mathcal{DN}$ is the space of constant functions.  
\end{theo}
The full symbol of  $\mathcal{DN}$ can be computed in
 a similar way from the full symbol of the resolvent $\Delta_g^{-1}$ along $Z$.



\section{The Poisson operator}

Let $A$ be an \OPD~ on $X$ of principal symbol $a$. We are
 interested to the restriction to the space $\mathcal{H}$ of 
 the quadratic form $Q_A(F)=\langle AF |F \rangle $ associated to $A
 $.
 We will parametrize $\mathcal{H}$ as harmonic extensions  of functions
 which are in $H^{-\ha}(Z)$ by the so-called Poisson operator denoted
 by $\mathcal{P}$; the pull-back $R_A$  of $Q_A$ on $L^2(Z)$ is defined by
\[
 R_A (f)=\langle A\mathcal{P}f \mid \mathcal{P}f \rangle =
 \langle \mathcal{P}^\star A\mathcal{P}f \mid f \rangle~.
\]
The goal of this section is to compute the operator
 $B=\mathcal{P}^\star A  \mathcal{P}$ associated to the quadratic form~$R_A$.

From Lemma \ref{lemm:harm-ext}, we have,   modulo smoothing operators,
 \[
  \mathcal{P}= \Delta_g^{-1}\circ \mathcal{E}\circ \mathcal{DN }~.
  \]
Hence
\[ 
B= \mathcal{DN }\circ \left[ \mathcal{T}\circ \left(\Delta_g^{-1}\circ A 
\circ\Delta_g^{-1}\right)\circ \mathcal{E}\right] \circ \mathcal{DN }~.
\]
   
The operator  $ \Delta_g^{-1}\circ A \circ\Delta_g^{-1}$
 is a \OPD~ of principal symbol $a/(\| \zeta\|_h^2 +\eta^2)^2$ near
 $Z$.
 

  Applying Theorem \ref{theo:main} to the inner bracket 
 and Theorem \ref{theo:dn}, we get the:
\begin{theo} \label{theo:poisson}
If $A$ is a \OPD ~ of degree $d<3$ on $X$ and $\mathcal{P}$
the Poisson operator associated to $Z$, the operator
 $B=\mathcal{P}^\star A  \mathcal{P}$ is a \OPD ~ of degree $d-1$ on $Z$
of principal symbol
\[
 b(z,\zeta)=\frac{2}{\pi} \| \zeta \|_h^2 \int _\R
 \frac{a(z,0;\zeta, \eta)}{(\| \zeta \|_h^2+ \eta ^2)^2}\dd \eta ~.
\]
\end{theo}
\begin{rema} Note that if $A$ is a pseudo-differential
 operator without the transmission property, 
the operator $A\circ \mathcal{P}$ may be ill-behaved and have disagreeable
 singularities along $Z$; however $\mathcal{P}^*A\mathcal{P}$ is always a
 good pseudo-differential operator on $Z$.
\end{rema}

\section*{Note}
This note was written with Louis in 2012. We  had the project  to publish  it as an Appendix to a work still in progress with Gregory Berkolaiko.  Finally, I decided to publish it independently and to dedicate it to the memory of Louis.

\section*{Declaration of interests}
The authors do not work for, advise, own shares in, or receive funds from
any organization that could benefit from this article, and have declared no affiliations other than their research organizations.
\printbibliography

\end{document}
 


\bibliographystyle{plain}
\begin{thebibliography}{99}

\bibitem{Cal} Alberto  Calder\'on. Boundary value problems for
  elliptic
 equations.  \emph{ Outlines of the joint Soviet-American Symposium
 on Partial Differential Equations, Novosibirsk, } pp 303--304  (1963).


\bibitem{LBM1} Louis Boutet de Monvel.
 Comportement d'un op\'erateur
 pseudo-diff\'erentiel sur une vari\'et\'e  \`a
 bord I, II. \emph{J. Anal. Math.} {\bf 17}:241--304 (1966).

\bibitem{LBM2} Louis Boutet de Monvel. Boundary problems
 for pseudo-differential operators.
\emph{Acta Mathematica} {\bf 126}:11--51 (1971).


\bibitem{Grub} Gerd Grubb. Distributions and Operators. \emph{Springer } (2008).

\bibitem{Ho} Lars H\"ormander. The Analysis of Linear Partial
 Differential Operators III.
Pseudo-Differential Operators. \emph{Grundlehren der Mathematischen
 Wissenschaften (Springer)}
{\bf 274} (1985).

\bibitem{Taylor} Michael E. Taylor. Partial Differential Equations II,
 Qualitative study of Linear Equations.
\emph{Springer, Applied Math. Sciences} {\bf 116} (1996).


\end{thebibliography}


\end{document}
