\documentclass[reqno]{amsart}
\usepackage{hyperref}

\AtBeginDocument{{\noindent\small
\emph{Electronic Journal of Differential Equations},
Vol. 2012 (2012), No. 189, pp. 1--30.\newline
ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu
\newline ftp ejde.math.txstate.edu}
\thanks{\copyright 2012 Texas State University - San Marcos.}
\vspace{9mm}}

\begin{document}
\title[\hfilneg EJDE-2012/189\hfil Generators with interior degeneracy]
{Generators with interior degeneracy on spaces of $L^2$ type}

\author[G. Fragnelli, G. R. Goldstein, J. A. Goldstein, S. Romanelli 
\hfil EJDE-2012/189\hfilneg]
{Genni Fragnelli, Gis\`ele Ruiz Goldstein \\
 Jerome A. Goldstein, Silvia Romanelli}  % in alphabetical order

\address{Genni Fragnelli \newline 
Dipartimento di Matematica\\
Universit\`{a} di Bari "Aldo Moro"\\
Via E. Orabona 4\\ 70125 Bari, Italy}
\email{genni.fragnelli@uniba.it}

\address{Gis\`ele Ruiz Goldstein \newline
Department of Mathematical Sciences\\
University of Memphis, Memphis, TN 38152, USA}
\email{ggoldste@memphis.edu}

\address{Jerome A. Goldstein \newline
Department of Mathematical Sciences \\
University of Memphis, Memphis, TN 38152, USA}
\email{jgoldste@memphis.edu}

\address{Silvia Romanelli \newline
Dipartimento di Matematica\\
Universit\`a di Bari ``Aldo Moro''\\
Via E. Orabona 4\\
70125 Bari, Italy}
\email{silvia.romanelli@uniba.it}

\thanks{Submitted September 6, 2012. Published October 31, 2012.}
\subjclass[2000]{47D06, 35K65, 47B25, 47N20}
\keywords{Degenerate operators in divergence and non divergence form;
\hfill\break\indent
linear and semiliner parabolic equations; linear hyperbolic equations}

\begin{abstract}
 We consider operators in divergence and in nondivergence form with
 degeneracy at the interior of the space domain. Characterizing the
 domain of the operators, we prove that they generate positive
 analytic semigroups on spaces of $L^2$ type. Finally, some
 applications to linear and semilinear parabolic evolution problems
 and to linear hyperbolic ones are presented.
\end{abstract}

\maketitle
\numberwithin{equation}{section}
\newtheorem{theorem}{Theorem}[section]
\newtheorem{lemma}[theorem]{Lemma}
\newtheorem{proposition}[theorem]{Proposition}
\newtheorem{corollary}[theorem]{Corollary}
\newtheorem{remark}[theorem]{Remark}
\newtheorem{definition}[theorem]{Definition}
\newtheorem{hypothesis}[theorem]{Hypothesis}
\allowdisplaybreaks

\def\R{{\mathbb{R}}}
\def\N{{\mathbb{N}}}

\section{Introduction}

This article is concerned with the generation property of a second
order ordinary differential degenerate operator in divergence or in
nondivergence form under Dirichlet boundary conditions in the real
setting. In particular, we consider both the operators
 $A_1u := (au')'$ and $A_2u:= au''$ with a suitable domain, where the function
$a$ vanishes at an interior point of the domain. Degenerate
parabolic operators  naturally arise in many problems. For some
examples involving degeneracy, let us recall some applications
arising in aeronautics (the Crocco equation, see, e.g., \cite{cmv}),
in physics (boundary layer models, see, e.g., \cite{br}), in
genetics (Wright-Fisher and Fleming-Viot models, see, e.g.,
\cite{mr,s}) and in mathematical finance
(Black-Merton-Scholes models, see, e.g., \cite{egg,gmr,mro}).
 Moreover, degenerate operators have been extensively
studied since Feller's investigations in \cite{F1,F2}, whose
main motivation was the probabilistic interest of the associated
parabolic equation for transition probabilities. After that, the
degenerate operator $A_1u$ or $A_2u$ has been studied under
different boundary conditions, see, for example,
\cite{cmp,cfr,crv, fggr,frt,f,gl,mp,t1,t2,v}. In
particular, \cite{frt,mp,t1,t2} develop a
functional analytic approach to the construction of Feller
semigroups generated by degenerate elliptic operators with Wentzell
boundary conditions. In \cite{gl}, the authors consider degenerate
operators with boundary conditions of Dirichlet, Neumann, periodic,
or nonlinear Robin type. In \cite{s2}, A. Stahel proves that the
Dirichlet parabolic problem associated to a degenerate operator in
divergence form with degeneracy at the boundary of the domain, has a
unique solution under suitable assumptions on the degenerate
function, using an approximation technique.
In \cite{acf, cfr,f}, the authors consider the degenerate
operator in divergence and in non divergence form with Dirichlet or
Neumann boundary conditions, giving more importance to
controllability problems of the associated parabolic evolution
equations. However, all the previous papers deal with a degenerate
operator with degeneracy at the boundary of the domain. For example,
as $a$, one can consider the double power function
\[
a(x)= x^k(1-x)^\alpha, \quad x \in [0,1],
\]
where $k$ and $\alpha$ are positive constants.

To the best of our knowledge, Stahel's paper \cite{s1} is the first
treating a problem with a degeneracy which may be interior. In
particular, Stahel considers a parabolic problem in $\R^N$ with
Dirichlet, Neumann or mixed boundary conditions, associated with a
$N\times N$ matrix $a$, which is positive definite and symmetric,
but whose smallest eigenvalue might converge to $0$ as the space
variable approaches a singular set contained in the closure of the
space domain. In this case, he proves that the corresponding
abstract Cauchy problem has a solution, provided that
$\underline{a}^{-1} \in L^q(\Omega, \R)$ for some $q >1$, where
\[
 \underline{a}(x):= \min \{a(x)\xi\cdot \xi: \|\xi \|=1\}.
\]
In the present paper we generalize his assumption when $N=1$ and the
degeneracy is interior. More
precisely, we shall admit two types of degeneracy for $a$, namely
weak and strong degeneracy according to the following definitions:

\begin{definition}\label{def1}\rm
  The operators $A_1u:=(au')'$ and $A_2u= au''$ are weakly degenerate if there exists
$x_0 \in (0,1)$ such that $a(x_0)=0$, $a>0$ on
$[0, 1]\setminus \{x_0\}$, $a\in C[0,1]$ and $\frac{1}{a} \in
L^1(0,1)$.
\end{definition}

For example, as $a$, we can consider $a(x)= |x- x_0|^{\alpha}, \;
0<\alpha<1$.

\begin{definition}\label{def2} \rm
The operators $A_1u :=(au')'$ and $A_2u= au''$ are strongly
degenerate if there exists $x_0 \in (0,1)$ such that $a(x_0)=0$,
$a>0$ on $[0, 1]\setminus \{x_0\}$, $a\in W^{1, \infty}(0,1)$ and
$\frac{1}{a} \not \in L^1(0,1)$.
\end{definition}

For example, as $a$, one can consider $a(x)= |x- x_0|^{\alpha}, \;\alpha\ge1$.

We remark that, while  in \cite{s1} only the existence of a solution
for the parabolic problem is considered, here we analyze in detail
the underlying degenerate operator in the spaces $L^2 (0,1)$ with or
without weight, proving that under suitable assumptions it is
nonpositive and selfadjoint, hence it generates an analytic
semigroup which is positivity preserving (see Theorems
\ref{prop}-\ref{prop-sur-AS}, \ref{theorem_analytic_nondivergence} -
\ref{teorema_derivata1}).

Moreover, in the strongly degenerate case, we are able to
characterize the domain of the operator both in divergence and in
non divergence cases (see Propositions \ref{domain}, \ref{domain4})
under, possibly, additional assumptions on the function $a$. We
point out that the generation property of the operator in
nondivergence form cannot be deduced by the generation property of
the operator in divergence form without additional assumptions on
the function $a$. In fact, the operator
\[
A_2u= au''
\]
can be recast using divergence form as follows
\begin{equation}\label{pri0}
A_2u=A_1u-a'u'
\end{equation}
at the price of adding the drift term $-a'u'$ to the divergence form
operator $A_1$ defined in Definition \ref{def1} or Definition
\ref{def2}. Such an addition has major consequences. For example, as
described in \cite{bfr}, degenerate operators of the form
\eqref{pri0}, where $a$ degenerates at the boundary of the domain,
generate a strongly continuous semigroup in $L^2(0,1)$ under the
structural assumption
\begin{equation}\label{ip}
|a'(x)|\le C\sqrt{a(x)}
\end{equation}
for any $x \in [0,1]$, where $C$  is a positive constant. Hence, it
is natural to expect that a similar result holds also if $a$ has an
interior degeneracy. Similar considerations also hold interchanging
the role of divergence and nondivergence operators. This is proved
in Section \ref{non_divergence}, where, however, we establish the
generation property of $A_2u= au''$ on a suitable weighted space
without any additional technical assumption on $a$. The generator
property of $A_1u= (au')'$ is proved without any further assumptions
on $a$ in Section \ref{divergence_form}.


The paper is organized in the following way. In Section
\ref{divergence_form} we consider the degenerate operator in
divergence form and we prove that it is nonpositive and selfadjoint
on $L^2(0,1).$ In Section \ref{non_divergence} we prove the same
result for the operator in nondivergence form on the space
$L^2_{1/a}(0,1)$. Moreover, under hypothesis \eqref{ip}, we
prove that both degenerate operators in nondivergence and in
divergence form generate analytic semigroups on $L^2(0,1)$ and
$L^2_{1/a}(0,1)$, respectively (see Theorems
\ref{teorema_derivata} and \ref{teorema_derivata1}). As a
consequence of the results proved in Sections \ref{divergence_form}
and \ref{non_divergence}, in Section \ref{applications} we obtain
existence results for linear and semilinear parabolic evolution
problems and linear hyperbolic evolution problems associated with
the operators under consideration. Finally, for the reader
convenience, in the Appendix we give some compactness theorems which
are crucial for our proofs in the semilinear parabolic cases. In
this paper we will always consider spaces of real valued functions.


\section{Divergence form}\label{divergence_form}
In this section we consider the operator in divergence form, that is
$A_1u=(au')'$, and we distinguish two cases: the weakly degenerate
case and the strongly one.

\subsection{Weakly degenerate operator}

Throughout this subsection we assume that the operator is weakly
degenerate.
To prove that $A_1$, with a suitable domain, generates a
strongly continuous semigroup, we introduce, as in \cite{acf}, the
following weighted space:
\begin{align*}
 H^1_a (0,1):=\big\{& u \in L^2(0,1) : u \text{ absolutely
continuous in } [0,1],\\
 &\sqrt{a} u' \in  L^2(0,1) \text{ and } u(0)=u(1)=0 \big\}
\end{align*}
with the norm
\[
\|u\|_{H^1_a(0,1)}:= \left(\|u\|^2_{L^2(0,1)} +
\|\sqrt{a}u'\|^2_{L^2(0,1)}\right)^{1/2}
\]
 and consider
\[
H^2_a(0,1) := \{ u \in H^1_a(0,1) : au' \in H^1(0,1)\}.
\]
Then, define the operator $A_1$ by
\[
 D(A_1)= H^2_a(0,1),
\]
 and for any $ u \in D(A_1)$,
 \[
 A_1u:=(au')'.
\]
We prove the following Green's formula:

\begin{lemma}\label{green}
For all $(u,v)\in H^2_a(0,1)\times H^1_a(0,1)$ one has
\begin{equation}\label{greenformula}
\int_0^1(au')' v dx= - \int_0^1 au'v' dx.
\end{equation}
\end{lemma}

\begin{proof}
Let $(u,v)\in H^2_a(0,1)\times H^1_a(0,1)$.
For any sufficiently small $\delta >0$ one has
\begin{equation}\label{integration}
\begin{aligned}
\int_0^1 (au')' v dx
&= \int_0^{x_0-\delta} (au')' v dx+
\int_{x_0-\delta}^{x_0+ \delta}(au')' v dx + \int_{x_0+
\delta}^1(au')' v dx \\
&=(au'v)(x_0- \delta) - (au'v)(0) \\
&\quad -\int_0^{x_0-\delta}au'v'dx +\int_{x_0- \delta}^{x_0+\delta}(au')'v dx  \\
&\quad+ (au'v)(1) - (au'v)(x_0+ \delta) - \int_{x_0+ \delta}^1 au'v' dx\\
&= (au'v)(x_0- \delta)  -\int_0^{x_0-\delta}au'v'dx +
\int_{x_0- \delta}^{x_0+\delta}(au')'v dx  \\
&\quad - (au'v)(x_0+ \delta) - \int_{x_0+ \delta}^1 au'v' dx,
\end{aligned}
\end{equation}
since $au' \in H^1(0,1)$ and $v(0)=v(1)=0$. Now, we prove that
\[
\lim_{\delta \to 0} \int_0^{x_0- \delta}au'v'  dx = \int_0^{x_0}au'v'  dx, \quad
\lim_{\delta \to 0} \int_{x_0+ \delta}^1au'v'  dx = \int_{x_0}^1au'v'  dx
\]
and
\begin{equation}\label{delta2}
\lim_{\delta \to 0} \int_{x_0- \delta}^{x_0+ \delta}(au')'v dx =0.
\end{equation}
Toward this end, observe that
\begin{equation}\label{delta1'}
 \int_0^{x_0- \delta}au'v' dx =
\int_0^{x_0}au'v' dx - \int_{x_0- \delta}^{x_0}au'v'
dx\end{equation}
and
\begin{equation}\label{delta1''}
\int_{x_0+ \delta}^1au'v' dx= \int_{x_0}^1au'v' dx -
\int_{x_0}^{x_0+ \delta}au'v' dx.
\end{equation}
Moreover, $(au')'v$ and  $au'v' \in L^1(0,1)$. Thus, for any
$\epsilon> 0$, by the absolute continuity of the integral, there
exists $\delta:= \delta(\epsilon) >0 $ such that
\begin{gather*}
\big| \int_{x_0- \delta}^{x_0}au'v'  dx \big|
\le \big| \int_{x_0- \delta}^{x_0}|au'v'  |dx \big| < \epsilon,
\\
\big|\int_{x_0- \delta}^{x_0+ \delta}(au')'v dx\big|
\le\big|\int_{x_0- \delta}^{x_0+ \delta}|(au')'v| dx\big|< \epsilon,
\\
\big|\int_{x_0}^{x_0+ \delta}au'v'  dx \big|
\le \big| \int_{x_0}^{x_0+ \delta}|au'v'  |dx \big| < \epsilon.
\end{gather*}
Now, take such a $\delta$ in \eqref{integration}. Thus, $\epsilon$
being arbitrary,
\[
\lim_{\delta \to 0} \int_{x_0- \delta}^{x_0}au'v'  dx =
\lim_{\delta \to 0} \int_{x_0- \delta}^{x_0+ \delta}(au')'v
dx = \lim_{\delta \to 0}\int_{x_0}^{x_0+ \delta}au'v'  dx
=0.
\]
The previous equalities and \eqref{delta1'}, \eqref{delta1''} imply
\[
\lim_{\delta \to 0} \int_0^{x_0- \delta}au'v'  dx =
\int_0^{x_0}au'v'  dx\quad \text{and} \quad \lim_{\delta \to
0} \int_{x_0+ \delta}^1au'v'  dx = \int_{x_0}^1au'v'  dx.
\]
 To obtain the desired result it is sufficient to prove that
\begin{equation}\label{delta3}
\lim_{\delta \to 0} (au'v)(x_0- \delta) =\lim_{\delta
\to 0} (au'v)(x_0+ \delta).
\end{equation}
 Since $au' \in H^1(0,1)$ and $v \in H^1_a(0,1)$,
\begin{equation}\label{equality}
\lim_{\delta \to 0} (au'v)(x_0- \delta) =(au'v)(x_0)
=\lim_{\delta \to 0} (au'v)(x_0+ \delta).\end{equation}
Thus, by \eqref{integration} -- \eqref{delta2} and \eqref{delta3},
it follows that
\[
\int_0^1(au')' v dx= - \int_0^1 au'v' dx. \qedhere
\]
\end{proof}
As a consequence of the previous lemma one has the next result.

\begin{theorem}\label{prop}
The operator $A_1: D(A_1) \to L^2(0,1)$ is
nonpositive and self-adjoint on $L^2(0,1)$. Moreover, the semigroup
$\{T(t)= e^{tA_1}: t \ge 0\}$ generated by $A_1$ is positivity
preserving.
\end{theorem}

\begin{proof}
Observe that $D(A_1)$ is dense in $L^2(0,1)$. To show that
$A_1$ is nonpositive and self-adjoint it suffices to prove that
$A_1$ is symmetric, nonpositive and $(I-A_1)(D(A_1))=L^2(0,1)$ (see
e.g. \cite[Theorem B.14]{abhn} or \cite{g}).

\noindent\textbf{$\boldsymbol A_1$  is symmetric.}
Indeed, for any $u,v\in D(A_1)$, one has
\[
\langle v, A_1u\rangle_{L^2(0,1)}= \int_0^1 v(au')' dx
=-\int_0^1 au'v' dx= \int_0^1(av')'u dx= \langle A_1v, u
\rangle_{L^2(0,1)}.
\]


\noindent\textbf{$\boldsymbol A_1$ is nonpositive.}
By \eqref{greenformula}, it follows that, for any $u\in D(A_1)$
\[
\langle A_1u, u \rangle_{L^2(0,1)} = \int_0^1 (au')' u dx =
-\int_0^1 a (u')^2 dx \le 0.
\]

\noindent \textbf{$\boldsymbol{ I-A_1}$ is surjective.}
First of all, observe that $H^1_a(0,1)$ equipped with the inner
product
\[
(u, v)_1:= \int_0^1 (u v + au'v') dx,
\]
for any $u, v \in H^1_a(0,1)$, is a Hilbert space. Moreover,
\begin{equation}\label{tripla}
H^1_a(0,1) \hookrightarrow L^2(0,1) \hookrightarrow (H^1_a(0,1))^*,
\end{equation}
where $(H^1_a(0,1))^*$ is the dual space of $H^1_a(0,1)$ with
respect to $L^2(0,1)$ (cf. \eqref{tripla}). Indeed, the continuous
embedding of $H^1_a(0,1)$ in $L^2(0,1)$ is readily seen. In
addition, for any $f, \varphi\in L^2(0,1)$
\[
|\langle f, \varphi\rangle_{L^2(0,1)}| =\left|\int_0^1 f\varphi dx
\right|\le \|f\|_{L^2(0,1)} \|\varphi\|_{L^2(0,1)} \le
\|f\|_{H^1_a(0,1)}\|\varphi\|_{L^2(0,1)}.
\]
Hence,
 $L^2(0,1) \hookrightarrow (H^1_a(0,1))^*$. Then,
$(H^1_a(0,1))^*$ is the completion of $L^2(0,1)$ with respect to the
norm of $(H^1_a(0,1))^*$. Now, for
 $f \in L^2(0,1)$,
consider the functional $F: H^1_a(0,1) \to \R$ defined as
$ F(v) := \int_0^1 f v dx$. Since $H^1_a(0,1)
\hookrightarrow L^2(0,1)$, we have that $F \in (H^1_a(0,1))^*$. As a
consequence, by Riesz's Theorem, there exists a unique $u \in
H^1_a(0,1)$ such that for all $v \in H^1_a(0,1)$
\begin{equation}\label{distributional}
(u, v)_1= \int_0^1 f v dx.
\end{equation}
In particular, since $C_c^\infty(0,1) \subset H^1_a(0,1)$,
\eqref{distributional} holds for all $v \in C_c^\infty(0,1)$; i.e.,
\[
\int_0^1au'v'dx = \int_0^1 (f-u)v dx, \quad \text{ for all }   v
\in C_c^\infty(0,1).
\]
Thus, the distributional derivative of $au'$ is a function in
$L^2(0,1)$, that is $au' \in H^1(0,1)$ (recall that $\sqrt{a}u' \in
L^2(0,1)$) and $(au')' = u-f$ a.e. in $(0,1)$. Then $u \in
H^2_a(0,1)$ and, by \eqref{distributional} and Lemma \ref{green}, we
have
\[
\int_0^1(u -(au')'- f)v dx =0.
\]
Consequently,
\begin{equation}\label{adjoint1}
u \in D(A_1) \quad \text{ and } \quad u- A_1u=f.
\end{equation}

As an immediate consequence of the Stone-von Neumann spectral
theorem and functional calculus associated with the spectral
theorem, one has that the operator $(A_1, D(A_1))$ generates a
cosine family and an analytic semigroup of angle $\frac{\pi}{2}$ on
$L^2(0,1)$. Positivity preserving follows as a consequence of the
positive minimum principle.
\end{proof}

\subsection{Strongly degenerate operator}

In this subsection we assume that the operator is strongly
degenerate. Following \cite{acf},  we introduce the weighted space
\[
\begin{aligned}
H^1_a(0,1):=\big\{& u \in L^2(0,1) : u \text{ locally absolutely
continuous in } [0,x_0) \cup (x_0,1], \\ &\sqrt{a} u' \in  L^2(0,1)
\text{ and } u(0)= u(1)=0 \big\}
\end{aligned}
\]
with the norm
\[
\|u\|_{H^1_a(0,1)}:= \left(\|u\|^2_{L^2(0,1)} +
\|\sqrt{a}u'\|^2_{L^2(0,1)}\right)^{1/2}.
\]
Define the operator $A_1$ by
 \[
D(A_1)= H^2_a(0,1),
\]
and for any  $ u \in D(A_1),$
\[
 A_1u:=(au')',
\]
where $H^2_a(0,1)$ is defined as before. Since in this case a
function $u \in H^2_a(0,1)$ is locally absolutely continuous in
$[0,1]\setminus\{x_0\}$ and not necessarily absolutely continuous in
$[0,1]$ as for the weakly degenerate case, equality \eqref{equality}
is not true a priori. Thus, we have to prove again the Green
formula. To do this, an idea is to characterize the domain of $A_1$.
The next results hold:

\begin{proposition}\label{characterization}
Let
\begin{align*}
X:=\Big\{& u \in L^2(0,1): u \text{ locally
absolutely continuous in } [0,1]\setminus \{x_0\},  \\
&\sqrt{a}u' \in L^2(0,1), au \text{ is continuous at } x_0 \;\text{and} \\
& (au)(x_0)=0=u(0)=u(1)\Big\}.
\end{align*}
Then
$H^1_a(0,1)=X$.
\end{proposition}

\begin{proof}
 Obviously, $X \subseteq H^1_a$. Now we take $u \in
H^1_a$, and we prove that $(au)(x_0)=0$, that is $u \in X$. Toward
this end, observe that $au \in H^1_0(0,1)$. Indeed, using the
assumptions on $u$, one has that $au \in L^2(0,1)$ and $(au)(0)=
(au)(1)=0$. Moreover, since $a \in W^{1, \infty}(0,1)$, $(au)' = a'
u+ au' \in L^2(0,1)$.  Thus $au \in H^1_0(0,1)\subset C[0,1]$. This
implies that there exists $\lim_{x \to x_0}(au)(x)=(au)(x_0)
=L \in \R$. If $L \neq 0$, then there exists $C>0$ such that
\[
|(au)(x)| \ge C
\]
for all $x$ in a neighborhood of $x_0$, $x\neq x_0$. Thus, setting
$C_1:= \frac{C^2}{\max_{[0,1]}a(x)} > 0$, it follows
that
\[
|u^2(x)| \ge \frac{C^2}{a^2(x)}\ge \frac{C_1}{a(x)},
\]
for all $x$ in a neighborhood of $x_0$, $x\neq x_0$. But, since the
operator is strongly degenerate, $ \frac{1}{a} \not \in
L^1(0,1)$ thus $u \not \in L^2(0,1)$. Hence $L=0$, that is
$(au)(x_0)=0$.
\end{proof}

Using the previous result, one can prove the following
characterization.

\begin{proposition}\label{domain}
Let
\begin{align*} 
D:=\big\{& u \in L^2(0,1) : u \text{ locally
absolutely continuous in } [0,1]\setminus \{x_0\}, \\
& au \in H^1_0(0,1), a u' \in H^1(0,1) \text{ and } (au)(x_0)=(au')(x_0)=0\big\}.
\end{align*}
Then
$ D(A_1)=D$.
\end{proposition}

To prove Proposition \ref{domain}, the following lemma is
crucial:

\begin{lemma} For all $u \in D$ we have that
\[
|a(x)u(x)| \le \|(au)'\|_{L^2(0,1)}\sqrt{|x-x_0|},
\]
and
\begin{equation}\label{domain2}
|a(x)u'(x)| \le \|(au')'\|_{L^2(0,1)}\sqrt{|x-x_0|},
\end{equation}
for all $x \in [0,1]$.
\end{lemma}

\begin{proof}
Let $u \in D$. Since $(au)(x_0)=0$, then
\[
|(au)(x)|= \big|\int_{x_0}^x (au)'(s)ds\big|
\le \|(au)'\|_{L^2(0,1)}\sqrt{|x-x_0|},
\]
for all $ x \in [0,1]$. Analogously, using the fact that
$(au')(x_0)=0$,
\[
|(au')(x) |=\big|\int_{x_0}^x(au')'(s) ds\big| \le
\|(au')'\|_{L^2(0,1)}\sqrt{|x-x_0|},
\]
for all $x \in [0,1]$.
\end{proof}

\begin{proof}[Proof of Proposition \ref{domain}]
 Let us prove that $D=D(A_1)$.

{$\boldsymbol{D\subseteq D(A_1):}$} Let $u \in D$. It is sufficient
to prove that $\sqrt{a}u' \in L^2(0,1)$. Since $au' \in H^1(0,1)$
and $u(1)=0$ (recall that $a >0$ in $[0,1]\setminus \{x_0\}$), for
$x \in (x_0,1]$ we have
\[
\int_x^1[(au')'u](s)ds= [au'u]_x^1-\int_x^1(a(u')^2)(s)ds =
-[au'u](x) -\int_x^1(a(u')^2)(s)ds.
\]
Thus
\[
(au'u)(x)=-\int_x^1[(au')'u](s)ds-\int_x^1(a(u')^2)(s)ds.
\]
Since $u \in D$, $(au')'u\in L^1(0,1)$. Hence, there exists
\[
\lim_{x \to x_0^+}(au'u)(x) =L\in[-\infty , +\infty).
\]
 If $L\neq 0$, there exists $C>0$ such that
\[
|(au'u)(x)|\ge C
\]
for all $x$ in a right neighborhood of $x_0$, $x \neq x_0$. Thus, by
\eqref{domain2}, there exists $C_1>0$ such that
\[
|u(x)| \ge \frac{C}{|(au')(x)|} \ge\frac{C_1}{\sqrt{(x-x_0)}},
\]
for all $x$ in a right neighborhood of $x_0$, $x\neq x_0$. This
implies that $u \not\in L^2(0,1)$. Hence $L=0$ and
\begin{equation}\label{int}
\int_{x_0}^1[(au')'u](s)ds=-\int_{x_0}^1(a(u')^2)(s)ds.
\end{equation}
 If $x\in
[0, x_0)$, proceeding as before and using the condition $u(0)=0$, it
follows that:
\[
(au'u)(x)=\int_0^x[(au')'u](s)ds+\int_0^x(a(u')^2)(s)ds
\]
and there exists
\[
 \lim_{x \to x_0^-}(au'u)(x) =L\in(-\infty ,
+\infty].
\]
\noindent As before, if $L\neq 0$, there exist two positive
constants $C'$ and $C_1'$ such that
\[
|u(x)| \ge \frac{C'}{|(au')(x)|} \ge\frac{C_1'}{\sqrt{(x_0-x)}},
\]
for all $x$ in a left neighborhood of $x_0$, $x \neq x_0$. This
implies that $u \not\in L^2(0,1)$. Hence $L=0$ and
\begin{equation}\label{int1}
\int_0^{x_0}[(au')'u](s)ds=-\int_0^{x_0}(a(u')^2)(s)ds.
\end{equation}
By \eqref{int} and \eqref{int1}, it follows that
\[
\int_0^1[(au')'u](s)ds=-\int_0^1(a(u')^2)(s)ds.
\]
 Since $(au')'u\in L^1(0,1)$, then $\sqrt{a}u' \in
L^2(0,1)$. Hence, $D \subseteq D(A_1)$.

{$\boldsymbol {D(A_1) \subseteq D:}$} Let $u \in D(A_1)$. As in the
proof of Proposition \ref{characterization}, we can prove that $au
\in H^1_0(0,1)$. Moreover, by Proposition \ref{characterization},
$(au)(x_0)=0$. Thus, it is sufficient to prove that $(au')(x_0)=0$.
Toward this end, observe that, since $au' \in H^1(0,1)$, there
exists $L\in \R$ such that $ \lim_{x \to
x_0}(au')(x)=(au')(x_0)= L$. If $L \neq 0$, there exists $C>0$ such
that
\[
|(au')(x)| \ge C,
\]
for all $x$ in a neighborhood of $x_0$, $x \neq x_0$. Thus
\[
|(a(u')^2)(x)| \ge \frac{C^2}{a(x)},
\]
for all $x$ in a neighborhood of $x_0$, $ x \neq x_0$. This implies
that $\sqrt{a}u' \not \in L^2(0,1)$. Hence $L=0$, that is
$(au')(x_0)=0$.
\end{proof}

We point out the fact that the condition $\frac{1}{a}
\not \in L^1(0,1)$ is crucial to prove the previous
characterizations. Clearly this condition is not satisfied if the
operator is weakly degenerate.

 As for the weakly degenerate case and using the previous characterization, we can prove the
following Green's formula.

\begin{lemma}\label{greenS}
For all $(u,v)\in H^2_a(0,1)\times H^1_a(0,1)$ one has
\[
\int_0^1(au')' v dx= - \int_0^1 au'v' dx.
\]
\end{lemma}

\begin{proof}
Let $(u,v)\in H^2_a(0,1)\times H^1_a(0,1)$.
As for the weak case, one can prove that, for any $\delta >0$:
\begin{equation}\label{integrationS}
\begin{aligned} \int_0^1 (au')' v dx
&= (au'v)(x_0- \delta) -\int_0^{x_0-\delta}au'v'dx\\
&\quad + \int_{x_0-\delta}^{x_0+\delta}(au')'v dx  - (au'v)(x_0+ \delta) - \int_{x_0+
\delta}^1 au'v' dx.
\end{aligned}
\end{equation}
Moreover,
\begin{equation}\label{delta1S}
\lim_{\delta \to 0} \int_0^{x_0- \delta}au'v'  dx =
\int_0^{x_0}au'v'  dx, \quad \lim_{\delta \to 0} \int_{x_0+
\delta}^1au'v'  dx = \int_{x_0}^1au'v' dx
\end{equation}
and
\begin{equation}\label{delta2S}
\lim_{\delta \to 0} \int_{x_0- \delta}^{x_0+ \delta}(au')'v
dx =0.
\end{equation}
 To obtain the
desired result it is sufficient to prove that
\begin{equation}\label{delta3S}
\lim_{\delta \to 0} (au'v)(x_0- \delta) =\lim_{\delta
\to 0} (au'v)(x_0+ \delta).
\end{equation}
First of all, observe that
\[
(au'v)(x_0-\delta)= \int_0^{x_0-\delta}((au')'v)(s)ds +
\int_0^{x_0-\delta}(au'v')(s)ds
\]
and
\[
(au'v)(x_0+\delta)= -\int_{x_0+\delta}^1((au')'v)(s)ds -
\int_{x_0+\delta}^1(au'v')(s)ds.
\]
 Since $(au')',\; v \in L^2(0,1)$ and
$\sqrt{a}u', \; \sqrt{a}v' \in L^2(0,1)$, by H\"{o}lder's
inequality, $(au')'v \in L^1(0,1)$ and $au'v'\in L^1(0,1)$. Thus,
there exist $L_1, \; L_2 \in \R$ such that
\begin{align*}
\lim_{\delta \to 0}(au'v)(x_0-\delta)&= \lim_{\delta
\to 0}\int_0^{x_0-\delta}((au')'v)(s)ds +\lim_{\delta
\to 0} \int_0^{x_0-\delta}(au'v')(s)ds\\
& =\int_0^{x_0}((au')'v)(s)ds + \int_0^{x_0}(au'v')(s)ds =L_1
\end{align*}
and
\begin{align*}
\lim_{\delta \to 0}(au'v)(x_0+\delta)&= -\lim_{\delta
\to 0}\int_{x_0+\delta}^1((au')'v)(s)ds -\lim_{\delta
\to 0} \int_{x_0+\delta}^1(au'v')(s)ds
\\&=-\int_{x_0}^1((au')'v)(s)ds-\int_{x_0}^1(au'v')(s)ds =L_2.
\end{align*}
If $L_1 \neq 0$, then there exists $C>0$ such that
\[
|(au'v)(x)| \ge C
\]
for all $x$ in a left neighborhood of $x_0$, $x \neq x_0$. Thus, by
\eqref{domain2},
\[
|v(x)| \ge \frac{C}{|(au')(x)|}\ge \frac{C_1}{\sqrt{x_0-x}}
\]
for all $x$ in a left neighborhood of $x_0$, $x \neq x_0$, and for a
suitable positive constant $C_1$. This implies that $v \not \in
L^2(0,1)$. Hence $L_1=0$. Analogously, one can prove that $L_2=0$.
Thus \eqref{delta3S} holds. In particular,
\[
\lim_{\delta \to 0} (au'v)(x_0- \delta) =\lim_{\delta
\to 0} (au'v)(x_0+ \delta)=0. \qedhere
\]
\end{proof}

As for the weakly degenerate case, one has the next result.

\begin{theorem}\label{prop-sur-AS}
The operator $A_1: D(A_1) \to L^2(0,1)$ is
self-adjoint and nonpositive on $L^2(0,1)$. Moreover, the semigroup
$\{T(t)= e^{tA_1}: t \ge 0\}$ generated by $A_1$ is positivity
preserving.
\end{theorem}


\section{Non divergence form}\label{non_divergence}

Now, we consider the operator $A_2u=au''$ and we distinguish again
between the weakly and the strongly degenerate cases.

\subsection{Weakly degenerate operator}\label{weakly_nondivergence}

Throughout this subsection we consider the weakly degenerate
operator and, as in \cite{cfr}, we consider the following Hilbert
spaces:
\begin{gather*}
L^2_{1/a}(0,1) :=\big\{ u \in L^2(0,1) :
 \int_0^1 \frac{u^2}{a} dx <\infty \big\},
\\
H^1_{1/a}(0,1) :=L^2_{1/a}(0,1)\cap H^1_0(0,1)
\end{gather*}
with the norms
\[
\|u\|_{L^2_{1/a}(0,1)}^2:= \int_0^1 \frac{u^2}{a} dx,
\]
and
\[
\|u\|_{H^1_{1/a}(0,1)} :=
\Big(\|u\|_{L^2_{1/a}(0,1)}^2 + \|u'\|^2_{L^2(0,1)}\Big)^{1/2},
\]
respectively. Using the
previous spaces, we define the operator $A_2$ by
$ D(A_2)= H^2_{1/a} (0,1)$,
 and for any $u \in D(A_2) $,
\[
A_2u:=au'',
\]
where
\[
H^2_{1/a}(0,1) :=\big\{ u \in H^1_{1/a}(0,1)\cap
W^{2,1}_{\text{loc}}(0,1) : au'' \in L^2_{1/a}(0,1)\big\}.
\]
 The following characterization is immediate.

\begin{corollary}
The spaces $H^1_{1/a}(0,1)$ and $H^1_0(0,1)$ coincide
algebraically. Moreover the two norms are equivalent.
\end{corollary}

\begin{proof}
Clearly $H^1_{1/a}(0,1) \subseteq H^1_0(0,1)$.
Now,  if $u \in H^1_0(0,1)$ then
\[
\int_0^1 \frac{u^2}{a} dx \le \max_{[0,1]}u^2\int_0^1 \frac{1}{a} dx
\in \R,
\]
using the fact that $ \frac{1}{a} \in L^1(0,1)$, that
is $u \in L^2_{1/a}(0,1)$. This implies $u\in
H^1_{1/a}(0,1)$.

Moreover, using the embedding of $H^1_0(0,1)$ in $C[0,1]$, one has,
for all  $u \in H^1_0(0,1)$,
\[
\int_0^1 \frac{u^2}{a} dx \le \|u\|_{C[0,1]}^2\|\frac{1}{a}
\|_{L^1(0,1)} \le C\|u\|_{H^1_0(0,1)}^2,
\]
for a positive constant $C$. Thus, for all  $u \in H^1_0(0,1)$,
\[
\|u\|_{H^1_{1/a}(0,1)} \le (C+1)\|u\|_{H^1_0(0,1)} \le
(C+1)\|u\|_{H^1_{1/a}(0,1)},
\]
for a positive constant $C$.
\end{proof}
As a consequence of the previous corollary one has that
$C_c^\infty(0,1)$ is dense in $H^1_{1/a}(0,1)$.


 Also for the weakly degenerate case, in order to prove that the operator generates
a cosine family and hence an analytic semigroup on $L^2_{1/a}(0,1)$,
 we have to use a Green formula similar to the one stated in Lemma
 \ref{green} or Lemma \ref{greenS}. However, we have to
prove it again, since we cannot adapt the proof of these Lemma. In
fact, to proceed as before, we would need $u', v' \in H^1(0,1)$, but
this is not the case.

The following Green's formula holds.

\begin{lemma}\label{green2}
 For all $(u,v)\in H^2_{1/a}(0,1)\times H^1_{1/a}(0,1)$
 one has
\begin{equation}\label{greenformula2}
\int_0^1u'' v dx= - \int_0^1 u'v' dx.
\end{equation}
\end{lemma}

\begin{proof}
First, we prove that  $H^1_c(0,1):=\big\{v\in
H^1(0,1):\operatorname{supp}\{v\}\subset (0,1)\setminus \{x_0\}\big\}$ is
dense in $H^1_{1/a}(0,1)$.
Indeed, if we consider the sequence $(v_n)_{n \ge
\max\left\{\frac{4}{x_0}, \frac{4}{1-x_0}\right\}}$, where $v_n:
=\xi_nv$ for a fixed function $v\in H^1_{1/a}(0,1)$ and
\[
\xi_n(x):=\begin{cases}
0,&x\in [0,1/n]\cup [ x_0-\frac{1}{n}, x_0+\frac{1}{n}]\cup [1-1/n,1],
\\
1,&x\in   [2/n, x_0 -2/n ] \cup [ x_0+2/n, 1-2/n],
\\
nx-1, &x\in  ( 1/n ,2/n ), \\
n(x_0-x) -1, & x \in  (x_0-2/n, x_0 -1/n ),\\
n(x-x_0)-1, & x \in (x_0+ 1/n, x_0+ 2/n),
\\
n(1-x)-1,  &x\in  (1-2/n ,1-1/n ),
\end{cases}
\]
it is easy to see that $v_n\to v$ in
$L^2_{1/a}(0,1)$. Moreover, one has that
\begin{equation}\label{h1c}
\begin{aligned}
&\int_0^1((v_n- v)')^2dx\\
&\le 2\int_0^1(1-\xi_n)^2(v')^2dx +2\int_0^1(\xi_n')^2v^2dx
\\
& = 2\int_0^1(1-\xi_n)^2(v')^2dx \\
&\quad +2n^2 \Big( \int_{1/n}^{2/n} v^2dx
 + \int_{x_0-\frac{2}{n}}^{x_0-\frac{1}{n}}  v^2dx
 + \int_{x_0+\frac{1}{n}}^{x_0+\frac{2}{n}}  v^2dx
 + \int_{1-\frac{2}{n}}^{1-\frac{1}{n}} v^2dx \Big).
\end{aligned}
\end{equation}
Obviously, the first term in the last member of \eqref{h1c}
converges to zero. Furthermore, since $v\in H^1_0(0,1)$, by
H\"older's inequality, one has that
$$
v^2(x)\le x\int_0^x (v')^2(y)dy,
$$
for all $ x \in [0,1]$. Therefore,
\begin{equation}\nonumber
n^2\int_{1/n}^{2/n}v^2dx \le
n^2\int_{0}^{2/n}(v')^2dx
\int_{1/n}^{2/n}x\,dx  \to   0,\quad \text{as } n\to\infty.
\end{equation}
 Since the remaining terms in \eqref{h1c} can be similarly estimated, our claim is proved.

Now, set $\varPhi(v):=\int_0^1(u'v)'dx$, with $u\in H^2_{1/a}(0,1)$. It follows that
\[
\varPhi(v) =0,
\]
 for all $v \in
H^1_c(0,1)$. In fact, let $v \in H^1_c(0,1)$. Let $\delta >0$ be
such that $v(x_0+ \delta)= v(x_0- \delta)=0$. Then, there holds
\begin{equation}\label{integration2}
\begin{aligned}
\int_0^1 u'' v dx
&= \int_0^{x_0-\delta} u''v dx+
\int_{x_0-\delta}^{x_0+ \delta}u'' v dx + \int_{x_0+
\delta}^1u'' v dx \\
&=(u'v)(x_0- \delta) - (u'v)(0) -\int_0^{x_0-\delta}u'v'dx +
\int_{x_0- \delta}^{x_0+\delta}u''v dx  \\
&\quad + (u'v)(1) -
(u'v)(x_0+ \delta) - \int_{x_0+ \delta}^1 u'v' dx\\
&= -\int_0^{x_0-\delta}u'v'dx + \int_{x_0- \delta}^{x_0+\delta}u''v
dx  - \int_{x_0+ \delta}^1 u'v' dx,
\end{aligned}
\end{equation}
since $v \in H^1_c(0,1)$. Now we prove that
\[
\lim_{\delta \to 0} \int_0^{x_0- \delta}u'v' dx =
\int_0^{x_0}u'v'  dx, \quad \lim_{\delta \to 0} \int_{x_0+
\delta}^1u'v'  dx = \int_{x_0}^1u'v'  dx
\]
and
\begin{equation}\label{delta2_2}
\lim_{\delta \to 0} \int_{x_0- \delta}^{x_0+ \delta}u''v dx
=0.
\end{equation}
Toward this end, observe that
\begin{equation}\label{delta1'_2}
 \int_0^{x_0- \delta}u'v' dx =
\int_0^{x_0}u'v' dx - \int_{x_0- \delta}^{x_0}u'v' dx\end{equation}
and
\begin{equation}\label{delta1''_2}
\int_{x_0+ \delta}^1u'v' dx= \int_{x_0}^1u'v' dx - \int_{x_0}^{x_0+
\delta}u'v' dx.
\end{equation}
Moreover, using the H\"{o}lder inequality, one can prove that $u''v$
and $u'v' \in L^1(0,1)$. Thus, for any $\epsilon> 0$, by the
absolute continuity of the integral, there exists
$\delta:=\delta(\epsilon)>0 $ such that
\begin{gather*}
\big| \int_{x_0- \delta}^{x_0}u'v'  dx \big|
\le \big| \int_{x_0-\delta}^{x_0}|u'v'  |dx \big| < \epsilon,
\\
\big|\int_{x_0- \delta}^{x_0+ \delta}u''v dx\big|
\le\big|\int_{x_0- \delta}^{x_0+ \delta}|u''v| dx\big|< \epsilon,
\\
\big|\int_{x_0}^{x_0+ \delta}u'v'  dx \big|
\le \big| \int_{x_0}^{x_0+ \delta}|u'v'  |dx \big| < \epsilon.
\end{gather*}
Now, take such a $\delta$ in \eqref{integration2}. Thus, $\epsilon$
being arbitrary,
\[
\lim_{\delta \to 0} \int_{x_0- \delta}^{x_0}u'v'  dx =
\lim_{\delta \to 0} \int_{x_0- \delta}^{x_0+ \delta}u''v dx
= \lim_{\delta \to 0}\int_{x_0}^{x_0+ \delta}u'v'  dx =0.
\]
The previous equalities and \eqref{delta1'_2}, \eqref{delta1''_2}
imply
\[
\lim_{\delta \to 0} \int_0^{x_0- \delta}u'v'  dx =
\int_0^{x_0}u'v'  dx\quad \text{and} \quad \lim_{\delta \to
0} \int_{x_0+ \delta}^1u'v'  dx = \int_{x_0}^1u'v' dx.
\]
Thus, by the previous equalities and by \eqref{integration2}--
\eqref{delta2_2}, it follows that
\[
\int_0^1u'' v dx= - \int_0^1 u'v' dx\quad \text{if and only if}\quad
\varPhi(v)=\int_0^1(u'v)'dx =0,
\]
for all $v \in H^1_c(0,1)$.  Then, $\varPhi$ is a bounded linear
functional on $H^1_{1/a}(0,1)$ such that $\varPhi=0$ on
$H^1_c(0,1)$. Thus, $\varPhi=0$ on $H^1_{1/a}(0,1)$, that
is, \eqref{greenformula2} holds.

\end{proof}
As a consequence of the previous lemma one has the next proposition,
whose proof is similar to the proof of Theorem \ref{prop}.

\begin{theorem}\label{theorem_analytic_nondivergence}
The operator $A_2: D(A_2) \to L^2_{1/a}(0,1)$ is
self-adjoint and nonpositive on $L^2_{1/a}(0,1)$.  Moreover,
the semigroup $\{T(t)= e^{tA_2}: t \ge 0\}$ generated by $A_2$ is
positivity preserving.
\end{theorem}


\subsection{Strongly degenerate operator}

 Assume that the operator
$A_2$ is strongly degenerate and consider, as in \cite{cfr}, the
spaces introduced in Section \ref{weakly_nondivergence}. Observe
that also in this case the space $C_c^\infty(0,1)$ is dense in
$H^1_{1/a}(0,1)$, since it is clearly dense in $H^1_0(0,1)$
and dense in $L^2_\frac{1}{a}(0,1)$. Then, the conclusions of Lemma
\ref{green2}  and Theorem \ref{theorem_analytic_nondivergence} hold
in this case.

\begin{theorem}\label{theorem_analytic_nondivergence_strongly}
The operator $A_2: D(A_2) \to L^2_{1/a}(0,1)$ is
self-adjoint and nonpositive on $L^2_{1/a}(0,1)$. Moreover,
the semigroup $\{T(t)= e^{tA_2}: t \ge 0\}$ generated by $A_2$ is
positivity preserving.
\end{theorem}

 Moreover in this case, under an additional assumption on the function $a$, one can
characterize the spaces
$H^1_{1/a}(0,1)$ and $H^2_{1/a}(0,1)$. We point out
the fact that in non divergence form, the characterization of the
domain of the operator is not important to prove the Green formula
as in divergence form.

From now on, we make the following assumption on $a$.

\begin{hypothesis} \label{Ass02} \rm
Assume that there exists a positive constant $K$ such that
$\frac{1}{a(x)} \le \frac{K}{|x-x_0|^2}$, for all $x
\in [0,1]\setminus\{x_0\}$
 (e.g. $a(x)=|x- x_0|^{K}, \; 1 \le K \le 2$).
\end{hypothesis}

\begin{proposition}\label{domain3}
Let
\[
\begin{aligned} X:=\{ u \in H^1_{1/a}(0,1) : u(x_0)=0\}.
\end{aligned}
\]
 If Hypothesis \ref{Ass02} is satisfied, then
\[
H^1_{1/a}(0,1)=X
\]
and the norm $\|u\|_{H^1_{\frac{1}{a}(0,1)}}$ and
$\big(\int_0^1(u')^2dx\big)^{1/2}$ are equivalent.
\end{proposition}

To prove Proposition \ref{domain3} the following lemma is crucial.

\begin{lemma}\label{hardy}
Assume that Hypothesis \ref{Ass02} is satisfied. Then, there exists
a positive constant $C$ such that
\[
\int_0^1 v^2\frac{1}{a} dx \le C\int_0^1 (v')^2dx,
\]
for all $v \in \, X$.
\end{lemma}

\begin{proof} Let $v\in X$.
By assumption, there exists $K>0$ such that
$\frac{1}{a(x)} \le \frac{K}{|x-x_0|^2}$, for all $x
\in [0,1]\setminus \{x_0\}$.
 Then, for  a suitable  $\varepsilon>0$ and using the assumption on $a$
and the Hardy inequality, one has
 \begin{align*}
 \int_0^1 v^2\frac{1}{a} dx &=\int_{0}^{x_0-\epsilon}v^2\frac{1}{a} dx
 + \int_{x_0-\epsilon}^{x_0+ \epsilon} v^2\frac{1}{a} dx +
 \int_{x_0+
 \epsilon}^{1} v^2\frac{1}{a} dx \\
 &\le
 \frac{1}{\min_{[0, x_0-\epsilon]}a(x)}\int_0^{x_0-\epsilon} v^2dx
 + \int_{x_0-\epsilon}^{x_0+\epsilon} v^2\frac{1}{a}dx\\
& \quad +
 \frac{1}{\min_{[x_0+\epsilon,1]}a(x)}\int_{x_0 +\epsilon}^1 v^2dx \\
& \le  \frac{1}{\min_{[0, x_0-\epsilon]\cup [x_0+\epsilon,1]}a(x)}\int_0^1 v^2 dx
 + \int_{x_0-\epsilon}^{x_0}v^2\frac{1}{a} dx  + \int_{x_0}^{x_0+ \epsilon}v^2\frac{1}{a} dx\\
& \le  \frac{1}{\min_{[0, x_0-\epsilon]\cup [x_0+\epsilon,1]}a(x)}\int_0^1 v^2 dx
 + K\int_{x_0-\epsilon}^{x_0}v^2\frac{1}{|x-x_0|^2} dx \\
&\quad +K\int_{x_0}^{x_0+ \epsilon}v^2\frac{1}{|x-x_0|^2} dx \\
& \le  \frac{1}{\min_{[0, x_0-\epsilon]\cup [x_0+\epsilon,1]}a(x)}\int_0^1 v^2 dx\\
&\quad + C_H \int_{x_0-\epsilon}^{x_0}  (v')^2 dx
 +C_H \int_{x_0}^{x_0+ \epsilon}  (v')^2dx \\
& \le  C \Big( \int_0^{1}v^2 dx  + \int_0^1 (v')^2 dx\Big),
 \end{align*}
for a positive constant $C$. Here $C_H$ is the Hardy constant.
 By Poincar\'e's inequality, it follows that
 \[
 \int_0^1 v^2\frac{1}{a} dx\le C \int_0^1 (v')^2dx
 \]
 for a suitable constant $C$.
\end{proof}

\begin{proof}[Proof of Proposition \ref{domain3}]
Obviously $X \subseteq H^1_{1/a}(0,1)$. Now, take $u \in
H^1_{1/a}(0,1)$ and prove that $u(x_0)=0$, that is $u \in
X$. Since $u \in H^1_0(0,1)$, then there exists
\[
\lim_{x \to x_0}u(x) =u(x_0)=L\in\R.
\]
\noindent If $L\neq 0$, then there exists $C>0$ such that
\[
|u(x)|\ge C,
\]
for all $x$ in a neighborhood of $x_0$, $ x \neq x_0$. Thus,
\[
\frac{u^2(x)}{a(x)} \ge \frac{C^2}{a(x)},
\]
for all $x$ in a neighborhood of $x_0$, $x \neq x_0$. Since the
operator is strongly degenerate, $ \frac{1}{a} \not\in
L^1(0,1)$, then $u \not\in L^2_{1/a}(0,1)$. Hence $L=0$.

Now, we prove that the two norms are equivalent. Take $u \in X$.
Then, by Lemma \ref{hardy}, there exists a positive constant $C$
such that
\[
\|u'\|^2_{L^2(0,1)} \le \|u\|^2_{H^1_{1/a}(0,1)}\le
C\|u'\|^2_{L^2(0,1)}.
\]
Hence the two norms are equivalent.
\end{proof}

 An immediate consequence is the following result.

\begin{proposition}\label{domain4}
Let
\[
D:=\big\{u \in H^1_{1/a}(0,1) : au'' \in
 L^2_{1/a}(0,1),\; au' \in H^1(0,1),\; 
 u(x_0)=(au')(x_0)=0\big\}.
\]
 If Hypothesis \ref{Ass02} is satisfied, then
$H^2_{1/a}(0,1)=D$.
\end{proposition}

\begin{proof}
Obviously, $D \subseteq H^2_{1/a}(0,1)$. Now, we take $u \in
H^2_{1/a}(0,1)$ and we prove that $u \in D$.  By Proposition
\ref{domain3}, $u(x_0)=0$. Thus, it is sufficient to prove that $au'
\in H^1(0,1)$ and $(au')(x_0)=0$. Since $u \in H^1_0(0,1)$ and $a
\in W^{1, \infty}(0,1)$, then $u' \in L^2(0,1)$, $\sqrt{a}u' \in
L^2(0,1)$ and $au' \in L^2(0,1)$. Moreover $(au')'=a'u'+ au'' \in
L^2(0,1)$ (recall that $au''\in L^2_{1/a}(0,1) \subset
L^2(0,1)$). Thus $au' \in H^1(0,1)$. This implies that there exists
$\lim_{x \to x_0}(au')(x) = (au')(x_0)=L \in \R$. If $L \neq
0$, then there exists $C>0$ such that
\[
|(au')(x)| \ge C
\]
for all $x$ in a neighborhood of $x_0$, $x \neq x_0$. Thus,
\[
|a(u')^{2}(x)| \ge \frac{C^2}{a(x)},
\]
for all $x$ in a neighborhood of $x_0$, $ x \neq x_0$. But
$ \frac{1}{a} \not \in L^1(0,1)$, thus $\sqrt{a}u' \not
\in L^2(0,1)$. Hence $L=0$, that is $(au')(x_0)=0$.
\end{proof}

We point out the fact that also in non divergence form the condition
$ \frac{1}{a} \not \in L^1(0,1)$ is crucial to
characterize the domain of the strongly degenerate operator.

\subsection{The operator in non divergence form in the space $\bf{L^2(0,1)}$}

In the previous subsections we have seen that without any additional
assumptions on the function $a$, the operator in non divergence form
generates a cosine family, and hence an analytic semigroup, on the
space $L^2_{1/a}(0,1)$.  In this subsection we will prove
that this operator generates an analytic semigroup also on
$L^2(0,1)$ under a suitable assumption on $a$. In particular we make
the following hypothesis:

\begin{hypothesis}\label{derivata}\rm
\begin{equation}\label{WD}
\begin{aligned}
&\text{In the  \emph{weakly degenerate case}, assume further  that } a \in W^{1,\infty} (0,1) \\
&\text{ and there exists } C \ge 0 \text{ such that } |a'(x)| \le C \sqrt{a(x)} \text{ a.e. }
x \in (0,1).
\end{aligned}
\end{equation}

\begin{equation}\label{SD}
\begin{aligned}
&\text{In the  \emph{strongly degenerate case}, assume further  that there exists }\\
& C \ge 0 \text{ such that } |a'(x)| \le C \sqrt{a(x)} \text{ a.e. } x \in (0,1).
\end{aligned}
\end{equation}
\end{hypothesis}

Now, define the operator $A$ by $ D(A)= H^2_{a} (0,1)$,
 and for any $ u \in D(A)$,
\[
Au:=au'',
\]
where $H^2_a(0,1)$ is the Sobolev space introduced in Section
\ref{divergence_form}. Observe that if $a$ satisfies Hypothesis
\ref{derivata}, then
\[
au'' \in L^2(0,1) \text{ if and only if } (au')' \in L^2(0,1),
\]
for all $u \in D(A)$.

\begin{theorem}\label{teorema_derivata}
Assume that the operator $A$ is weakly
or strongly degenerate. If $a$ satisfies Hypothesis \ref{derivata},
then $A$ generates an analytic semigroup on $L^2(0,1)$.
\end{theorem}

\begin{proof}
Let $Bu:= a'u'$, for all $u \in D(A)$. Then $Au= A_1u-Bu$, where, we recall,
$A_1u:=(au')'$. Therefore, one has that $B$ is $A_1 -$ bounded with
$A_1 -$ bound $b_0=0$, where, we recall,
\begin{equation}\label{bound}
\begin{aligned}
b_0 := \inf\big\{& b \ge 0: \text{ there exists } c \in \R_+ \text{ such
that }\\
& \|Bu\|_{L^2(0,1)}\le b\|A_1 u\|_{L^2(0,1)} +c \|u\|_{L^2(0,1)}\big\}
\end{aligned}
\end{equation}
(see, for example, \cite[Definition III.2.1]{en}). Indeed, using the
assumption on $a$, one has
\begin{align*}
\|Bu\|_{L^2(0,1)}^2
&= \int_0^1(a')^2(u')^2dx \le C\int_0^1au'u' dx\\
&= -C\int_0^1 (au')'udx  \le \epsilon \frac{K}{2}\int_0^1
((au')')^2dx + \frac{K}{\epsilon}\int_0^1 u^2dx\\
&=\epsilon \frac{K}{2}\|A_1 u\|_{L^2(0,1)}^2 +\frac{K}{\epsilon}\|u\|_{L^2(0,1)}^2
\end{align*}
for all $u \in D(A_1)$, for all $\epsilon >0$ and for a positive
constant $K$. Hence $B$ is a Kato perturbation of $A_1$.
 The conclusion follows by \cite[Theorem III.2.10]{en} and Theorems \ref{prop} and \ref{prop-sur-AS}.
\end{proof}

\subsection{The operator in divergence form in the space $\bf{L^2_{1/a}(0,1)}$}

As in the previous subsection, interchanging the role of the
divergence and nondivergence operators, we can prove that the
operator in divergence form generates an analytic semigroup also on
$L^2_{1/a}(0,1)$ under Hypothesis \ref{derivata}. Indeed,
define the operator $A$ by
$ D(A)= H^2_{1/a} (0,1)$,
 and for any $u \in D(A)$,
\[
Au:=(au')',
\]
where $H^2_{1/a}(0,1)$ is the Sobolev space introduced in
Section \ref{non_divergence}. Also in this case we have that if $a$
satisfies Hypothesis \ref{derivata}, then
\[
au'' \in L^2_{1/a}(0,1) \text{ if and only if } (au')' \in L^2_{1/a}(0,1),
\]
for all $u \in D(A)$. The following theorem holds.

\begin{theorem}\label{teorema_derivata1}
Assume that the operator $A$ is weakly
or strongly degenerate. If $a$ satisfies Hypothesis \ref{derivata},
then $A$ generates an analytic semigroup on
$L^2_{1/a}(0,1)$.
\end{theorem}

\begin{proof}
Let $Bu:= a'u'$, for all $u \in D(A)$. Then $Au= A_2u+Bu$, where, we recall,
$A_2u:=au''$. As before, $B$ is $A_2 -$ bounded with $A_2 -$ bound
$b_0=0$, where $b_0$ is defined in \eqref{bound}. Indeed, using the
assumption on $a$, one has
\begin{align*}
\|Bu\|_{L^2_{1/a}(0,1)}^2 &=
\int_0^1\frac{(a')^2(u')^2}{a}dx \le C\int_0^1u'u' dx\\&= -C\int_0^1
u''udx =-C\int_0^1\sqrt{a}u''\frac{u}{\sqrt{a}}dx \\&\le \epsilon
\frac{K}{2}\int_0^1
a(u'')^2dx + \frac{K}{\epsilon}\int_0^1 \frac{u^2}{a}dx\\
&=\epsilon \frac{K}{2}\|A_2 u\|_{L^2_{1/a}(0,1)}^2
+\frac{K}{\epsilon}\|u\|_{L^2_{1/a}(0,1)}^2
\end{align*}
for all $u \in D(A_2)$, for all $\epsilon >0$ and for a positive
constant $K$. Hence $B$ is a Kato perturbation of $A_2$.
 The conclusion follows by \cite[Theorem III.2.10]{en} and Theorems
\ref{theorem_analytic_nondivergence} and
\ref{theorem_analytic_nondivergence_strongly}.
\end{proof}

\section{Applications}\label{applications}

\subsection{Linear problems}

As an application of the previous theorems, consider the linear
operator $B(t)$ defined as
\[
B(t) u := -b(t,\cdot)\frac{\partial{u}}{\partial x}- c(t,\cdot) u,
\]
where $b, c\in L^\infty(\R_+\times(0,1))$ and there exists a
positive constant $C$ such that $|b(t,x)| \le C \sqrt{a(x)}$, a.e.
$x \in (0,1)$. Then $B(t)$ is an $A_1 -$ bounded operator on
$L^2(0,1)$ or an $A_2-$ bounded operator on
$L^2_{1/a}(0,1)$. Observe that if $b \equiv 0$ and
$c(t,\cdot) \le 0$, the operator $A_i-B(t)$ with domain $D(A_i)$,
$i=1,2$, is still selfadjoint and nonpositive for each $t$.
Moreover, $D(A_i-B(t))$ is independent of $t$. Thus using evolution
operator theory (see e.g. \cite{g}, pp. 140-147), we can prove that
the problem
\begin{equation}\label{parabolic}
\begin{gathered}
{\frac{\partial u}{\partial t}} - Au + b(t,x)\frac{\partial{u}}{\partial x}+ c(t,x)u =h(t,x),
\quad  (t,x) \in \R_+ \times (0,1),\\
u(t,0)=u(t,1)=0, \quad  t \ge 0, \\
u(0,x)=u_0(x), \quad x \in (0,1),
\end{gathered}
\end{equation}
is wellposed in the sense of evolution operator theory, provided
that
 $A:=A_1$ or $A:=A_2$.
In addition, setting $Q_T := (0,T) \times (0,1)$ for a fixed $T>0$,
the next results follow by Theorems \ref{prop}, \ref{prop-sur-AS},
\ref{theorem_analytic_nondivergence},
\ref{theorem_analytic_nondivergence_strongly}.

\begin{theorem}\label{th-parabolic}
Assume that the operator $A=A_1$ is
weakly or strongly degenerate. If $b(\cdot, x), \, c(\cdot, x)\in
C^1(\R_+)$ for all $x \in [0,1]$ and there exists a positive
constant $C$ such that $|b(t,x)| \le C \sqrt{a(x)}$, a.e. $x \in
(0,1)$, then, for all $h \in  L^2(Q_T)$ and $u_0 \in L^2(0,1)$,
there exists a unique weak solution  $u \in C([0,T];
L^2(0,1)) \cap L^2 (0,T; H^1_a(0,1))$ of \eqref{parabolic} and
\begin{equation}\label{stima}
\sup_{t \in [0,T]} \|u(t)\|^2_{L^2(0,1)}+\int_0^T\|u(t)\|^2_{H^1_a
(0,1)} dt \le C_T(\|u_0\|^2_{L^2(0,1)}+\|h\|^2_{L^2(Q_T)}),
\end{equation}
for a positive constant $C_T$.
 Moreover, if $u_0 \in D(A_1)$,
then
\begin{equation}\label{regularity}
u \in H^1(0,T; L^2(0,1)) \cap L^2(0,T; H^2_a(0,1))\cap
C([0,T]; H^1_a(0,1)),
\end{equation}
 and there exists a positive constant $C$ such
that
\begin{equation}\label{stima1}
\begin{aligned}
&\sup_{t \in [0,T]}\left(\|u(t)\|^2_{H^1_a(0,1)} \right)+ \int_0^{T}
\Big(\big\|\frac{\partial u}{\partial t}\big\|^2_{L^2(0,1)} +
\big\|\frac{\partial }{\partial x}\big(a\frac{\partial u}{\partial
x}\big)\big\|^2_{L^2(0,1)}\Big)dt\\
&\le C \left(\|u_0\|^2_{H^1_a(0,1)} + \|h\|^2_{L^2(Q_T)}\right).
\end{aligned}
\end{equation}
\end{theorem}

\begin{remark}\label{remark} \rm
Actually, since $D(A_1)$ is dense in $H^1_a(0,1)$, one can prove
that \eqref{regularity} and \eqref{stima1} also hold if $u_0 \in
H^1_a(0,1)$.
\end{remark}

\begin{theorem}\label{theorem_nondivergence}
Assume that the operator $A=A_2$ is weakly or strongly degenerate.
If $b(\cdot, x), c(\cdot,x) \in C^1(\R_+)$ for all $x \in [0,1]$ and
there exists a positive constant $C$ such that $|b(t,x)| \le C
\sqrt{a(x)}$, a.e. $x \in (0,1)$, then, for all $h \in
L^2_{1/a}(Q_T)$ and $u_0 \in L^2_{1/a}(0,1)$,
there exists a unique weak solution $u \in C\big([0,T];
L^2_{1/a}(0,1)\big) \cap L^2 \big(0,T;
H^1_{1/a}(0,1)\big)$ of \eqref{parabolic} and
\begin{equation}\label{stima2}
\sup_{t \in [0,T]}
\|u(t)\|^2_{L^2_{1/a}(0,1)}+\int_0^T\|u(t)\|^2_{H^1_{1/a}
(0,1)} dt \le
C_T\Big(\|u_0\|^2_{L^2_{1/a}(0,1)}+\|h\|^2_{L^2_{1/a}(Q_T)}\Big),
\end{equation}
for a positive constant $C_T$. Moreover, if $u_0 \in D(A_2)$, then
\begin{equation}\label{regularity1}
u\in H^1\Big(0,T; L^2_{1/a}(0,1)\Big) 
\cap L^2\Big(0,T;H^2_{1/a}(0,1)\Big)\cap C\Big([0,T];
H^1_{1/a}(0,1)\Big),
\end{equation}
 and there exists a positive constant
$C$ such that
\begin{equation}\label{stima3}
\begin{aligned}
&\sup_{t \in [0,T]}\Big(\|u(t)\|^2_{H^1_{1/a}(0,1)}\Big)
+ \int_0^{T} \Big(\big\|\frac{\partial u}{\partial
t}\big\|^2_{L^2_{1/a}(0,1)} + \big\|a\frac{\partial^2
u}{\partial x^2}\big\|^2_{L^2_{1/a}(0,1)}\Big)dt\\
&\le C \Big(\|u_0\|^2_{H^1_{1/a}(0,1)} + \|h\|^2_{L^2(Q_T)}\Big).
\end{aligned}
\end{equation}
\end{theorem}

\begin{remark}\label{remark1} \rm
Also in this case \eqref{regularity1} and \eqref{stima3} hold if
$u_0\in H^1_{1/a}(0,1)$,  since $D(A_2)$ is dense in  $H^1_{1/a}(0,1)$.
\end{remark}

 By Theorem \ref{teorema_derivata}, one has the following result.

\begin{theorem}\label{th-parabolic_nondiv}
Assume that the operator $A=A_2$ is
weakly or strongly degenerate. If $a$ satisfies Hypothesis
\ref{derivata}, $b(\cdot, x), \, c(\cdot,x) \in C^1(\R_+)$ for all
$x \in [0,1]$ and there exists a positive constant $C$ such that
$|b(t,x)| \le C \sqrt{a(x)}$, a.e. $x \in (0,1)$,  then, for all $h
\in L^2(Q_T)$ and $u_0 \in L^2(0,1)$, there exists a unique weak
solution  $u \in  C([0,T]; L^2(0,1)) \cap L^2 (0,T;
H^1_a(0,1))$ of \eqref{parabolic}. If $u_0 \in D(A_1)$, then
$$
u \in H^1\big(0,T; L^2(0,1)\big) \cap L^2\big(0,T;
H^2_a(0,1)\big)\cap  C\big([0,T]; H^1_a(0,1)\big).
$$
Moreover, \eqref{stima} and \eqref{stima1} hold.
\end{theorem}

By Theorem \ref{teorema_derivata1}, one has the following result.

\begin{theorem}\label{th-parabolic_div} Assume that the operator $A=A_1$ is
weakly or strongly degenerate. If $a$ satisfies Hypothesis
\ref{derivata}, $b(\cdot, x), \, c(\cdot,x) \in C^1(\R_+)$ for all
$x \in [0,1]$ and there exists a positive constant $C$ such that
$|b(t,x)| \le C \sqrt{a(x)}$, a.e. $x \in (0,1)$,  then, for all $h
\in L^2_{1/a}(Q_T)$ and $u_0 \in L^2_{1/a}(0,1)$,
there exists a unique weak solution $u \in C\left([0,T];
L^2_{1/a}(0,1)\right) \cap L^2\left (0,T;
H^1_{1/a}(0,1)\right)$ of \eqref{parabolic}. If $u_0 \in
D(A_2)$, then
$$
u \in H^1\big(0,T; L^2_{1/a}(0,1)\big) \cap
 L^2\big(0,T;H^2_{1/a}(0,1)\big)\cap C\big([0,T];H^1_{1/a}(0,1)\big).
$$
Moreover, \eqref{stima2} and \eqref{stima3} hold.
\end{theorem}


Let \[B(t)u:=-c(t, \cdot) u,\] where $c \in
L^\infty_{\text{loc}}(\R_+ \times[0,1])$ and $c(\cdot, x) \in
C^1(\R_+)$ for each $x \in [0,1]$. Consider the nonautonomous wave
equation
\begin{equation}\label{J}
\begin{gathered}
{  \frac{\partial ^2 u}{\partial t^2}} - Au  +c(t,x){ \frac{\partial u}{\partial t}} =h(t,x), \quad
  (t,x) \in \R_+ \times (0,1),\\
u(t,0)=u(t,1)=0, \quad t \ge 0, \\
u(0,x)=u_0(x), \quad  x \in (0,1), \\
{ \frac{\partial u}{\partial t}}(0,x)=u_1(x), \quad  x \in (0,1).
\end{gathered}
\end{equation}
We take $A=A_j$ acting in $H_j$, for $j=1,2$, where
\[
H_1 =L^2(0,1), \quad H_2 = L^2_{1/a}(0,1).
\]
In either case $A_j= A_j^* \le 0$. Let $D_j$ be the completion of
$D(A_j)$ in the norm
\[
\|f\|_{D_j}= \|(-A_j)^{1/2}f\|_{H_j}.
\]
Then \eqref{J} can be rewritten as
\begin{gather*}
\frac{d}{dt} U(t) = GU(t) +P(t)U(t) + H(t), \\
U(0)=U_0=\begin{pmatrix}u_0\\u_1\end{pmatrix},
\end{gather*}
where
\[
G = \begin{pmatrix} 0&I\\A_j &0\end{pmatrix},\quad
P(t)= \begin{pmatrix}0&0\\ 0& - c(t, \cdot)\end{pmatrix}, \quad
H(t)= \begin{pmatrix}0\\h(t, \cdot) \end{pmatrix},
\]
for $j=1, 2$.
 By \cite[Section II.7.1]{g}, on the completion $K_j$ of
 $
 D_j \oplus H_j$ in the norm
 \[
\big\| \begin{pmatrix} v \\w\end{pmatrix} \big\|_{K_j} =
\left(\|(-A_j)^{1/2}v\|_{H_j}^2 + \|w\|^2_{H_j}
\right)^{1/2},
 \]
 $G$ is skewadjoint and generates a ($C_0$) unitary group. Moreover, each $P(t)$ is
bounded and $P\in C^1(\R_+, K_j)$. Finally, we assume
\[
h \in L^\infty_{\text{loc}} (\R_+ \times [0,1])
\]
and $h(\cdot, x) \in C^1(\R_+)$ for each $x \in [0,1]$. Then it
follows that
\[
H\in C^1(\R_+, K_j).
\]
The wellposedness of \eqref{J} now follows immediately by a simple
modification of \cite[Theorem II.13.9]{g}.


\subsection{Semilinear problems}

In this subsection we extend the existence results obtained in the
previous theorems to semilinear degenerate parabolic systems of the
type:
\begin{equation} \label{nl}
\begin{gathered}
u_t - Au  + f(t,x,u,u_x)  =h(t,x), \quad (t,x) \in (0,T) \times (0,1), \\
u(t,1)=u(t,0)=0, \quad t \in (0,T),\\
u(0,x)=u_0(x) , \quad  x \in (0,1),
\end{gathered}
\end{equation}
where $Au:=A_1u$ or $Au:=A_2u$,
$u_t(t,x):= \frac{\partial u (t,x)}{\partial t}$ and
$u_x(t,x):=  \frac{\partial u (t,x)}{\partial x}$. In particular,
recalling that $Q_T=(0,T)\times(0,1)$, we give the following
definitions:

\begin{definition} \label{def_div} \rm
Assume that $Au:=A_1u$, $u_0 \in L^2(0,1)$ and $h \in L^2(Q_T)$.
A function $u$ is said to be a solution of \eqref{nl} if
\[
u \in C(0, T; L^2(0,1)) \cap L^2(0, T; H^1_a(0,1))
\]
and satisfies
\begin{align*}
&\int_0^1 u(T,x)\varphi(T,x) dx - \int_0^1 u_0(x) \varphi(0,x) dx -
\int_0^T \int_0^1 \varphi_t (t,x)u(t,x)\,dx\,dt\\
& = - \int_0^T\int_0^1 a  u_x
\varphi_x\,dx\,dt + \int_0^T\int_0^1(-f(t,x,u,u_x)+ h(t,x)) \varphi(t,x)\,dx\,dt,
\end{align*}
for all $\varphi \in H^1(0, T; L^2(0,1)) \cap L^2(0, T;H^1_a(0,1))$.
\end{definition}

\begin{definition} \label{def_nondiv}\rm
 Assume that $Au:=A_2u$, $u_0 \in L^2_{1/a}(0,1)$ and $h\in L^2_{1/a}(Q_T)$.
A function $u$ is said to be a solution of \eqref{nl} if
 \[
u \in C(0, T; L^2_{1/a}(0,1)) \cap L^2(0, T;
H^1_{1/a}(0,1))
\]
and satisfies
\begin{align*}
&\int_0^1 \frac{ u(T,x)\varphi(T,x)}{a(x)} dx - \int_0^1
\frac{u_0(x) \varphi(0,x)}{a(x)} dx - \int_0^T \int_0^1
\frac{\varphi_t (t,x)u(t,x)}{a(x)}\,dx\,dt\\
& = - \int_0^T\int_0^1 u_x
\varphi_x\,dx\,dt + \int_0^T\int_0^1(-f(t,x,u,u_x)+ h(t,x))
\frac{\varphi(t,x) }{a(x)}\,dx\,dt,
\end{align*}
for all $\varphi \in H^1(0, T; L^2_{1/a}(0,1)) \cap L^2(0, T; H^1_{1/a}(0,1))$.
\end{definition}

On the functions $a$ and $f$ we make the following assumptions:

\begin{hypothesis}\label{5.0}  \rm
There exists $x_0 \in (0,1)$ such that $a(x_0)=0$, $a>0$ on $[0, 1]\setminus
\{x_0\}$, $a\in  C[0,1]$ and $\frac{1}{a} \in L^1(0,1)$.
\end{hypothesis}

\begin{hypothesis}\label{5.1} \rm
Let $f : [0,T] \times [0,1]\times \mathbb R \to \mathbb R$ be such
that
\begin{gather*}
\forall (q,p) \in \mathbb R^2, \quad (t,x) \mapsto
f(t,x,q,p)\text{ is measurable; }
\\
\text{for a.e. }  (t,x)  \in (0,T) \times  (0,1), \quad f(t, x,0,0)=0;
\end{gather*}
 $f_p(t,x,q,p)$ exists, $f_p$ is a Carath\'{e}odory function; i.e.,
\begin{gather*}
\forall (q,p) \in \mathbb R^2, \quad (t,x) \mapsto f_p(t,x,,q,p)  \text{ is measurable and}\\
\text{for a.e. }  (t,x) \in (0,T) \times (0,1),  \quad (q,p)
\mapsto f_p(t,x,q,p) \text{ is continuous},
\end{gather*}
and  there exists $L >0$ such that $\text{for a.e.} (t,x) \in (0, T) \times (0,1)$ and
for any $(q, p) \in \mathbb{R} \times \mathbb{R} $,
\begin{equation}\label{ipo.probabilistica}
 |f_p(t,x,q,p)| \le L \sqrt{a(x)};
\end{equation}
$f_q(t,x,q,p)$ exists, $f_q$ is a Carath\'eodory function and
 there exists $C >0$ such that $\text{for a.e.} (t,x) \in (0, T) \times (0,1)$ and for all
$(q, p) \in \mathbb{R} \times \mathbb{R} $
\begin{equation}
\label{f2}  |f_q(t,x,q,p)| \le C.
\end{equation}
\end{hypothesis}
However, to prove that \eqref{nl} has a solution in the non
divergence case, i.e. $Au:= A_2u$, it is sufficient to substitute
\eqref{ipo.probabilistica} with the more general condition: there
exists $ L >0 $ such that $\text{ for a.e. } (t,x) \in (0, T) \times
(0,1)$ and $\forall \, (q, p) \in \mathbb{R} \times \mathbb{R} $
\begin{equation}\label{ipo}
|f_p(t,x,q,p)| \le L.
\end{equation}

 Moreover,  to obtain the desired result, we
introduce on the spaces $H^2_a(0,1)$ and $H^2_{1/a}(0,1)$
the scalar products
\[
\int_0^1 u(x)v(x)dx
+\int_0^1a(x)u'(x)v'(x)dx+\int_0^1(a(x)u'(x))'(a(x)v'(x))'dx,
\]
for all $u,v \in H^2_a(0,1)$, and
\[
\int_0^1\frac{u(x)v(x)}{a(x)}dx+\int_0^1u'(x)v'(x)dx+\int_0^1au''(x)v''(x)dx,
\]
for all $u, v\in H^2_{1/a}(0,1)$, respectively. The previous
scalar products induce on $H^2_a(0,1)$ and $H^2_{1/a}(0,1)$
the following two norms
\begin{gather*}
\|u\|_{H^2_a(0,1)} := \Big(\|u\|^2_{H^1_a(0,1)} +
\|(au')'\|^2_{L^2(0,1)}\Big)^{1/2},
\\
\|u\|_{H^2_{1/a}(0,1)} :=\Big(\|u\|_{H^1_{1/a}(0,1)}^2 +
\|au''\|^2_{L^2_{1/a}(0,1)}\Big)^{1/2}.
\end{gather*}
Clearly, $H^2_a(0,1)$ and $H^2_{1/a}(0,1)$ are Hilbert
spaces.

As a first step, we study  \eqref{nl} with  $u_0 \in H^1_a(0,1)$, if
$Au:= A_1u$, or $u_0 \in H^1_{1/a}(0,1)$, if $Au:= A_2u$.

To prove the existence results we will use, as in \cite{acf} or in \cite{cf}, a
fixed  point method. To this aim, we rewrite, first of all, the function $f$
in the following way $f(t,x,u, u_x) = b(t,x, u)u_x +  c(t,x,u)u$, where
\[
       b(t,x,u) := \int_0^1 f_p(t,x,\lambda u,\lambda u_x)d\lambda,
       \quad
      c(t,x,u) := \int_0^1 f_q(t,x,\lambda u,\lambda
       u_x)d\lambda.
\]
In fact
\begin{align*}
       f(t,x, u, u_x) &= \int_0^1\frac{d}{d\lambda}f(t,x,\lambda
       u,\lambda u_x)d\lambda \\
&= \int_0^1 f_q(t,x,\lambda u,\lambda
       u_x)ud\lambda  + \int_0^1 f_p(t,x,\lambda u,\lambda
       u_x)u_xd\lambda.
\end{align*}
Using the fact that $f_p$ and $f_q$ are Carath\'eodory functions
and the Lebesgue Theorem, we can prove the following properties:

\begin{proposition}\label{prop1}
For the functions $b$ and $c$ one has the following properties:
\begin{itemize}
\item $b(t,x,u(t,x))$ and $c(t,x, u(t,x))$ belong to
       $L^{\infty}((0, T) \times (0,1))$;

\item $|b(t,x,u)| \le L
       \sqrt{a(x)}$;

\item if $\lim_{k \to + \infty}v_k = v$ in $X$, then
\[
\lim_{k \to + \infty} \frac{b(t,x; v_k)}{\sqrt{a(x)}} =\frac{b(t,x;
v)}{\sqrt{a(x)}}, \quad \text{a.e.},
\]
and
\[
\lim_{k \to + \infty}c(t,x; v_k) =c(t,x; v), \quad \text{a.e.}
\]
\end{itemize}
Here
 $$
X:=C(0, T; L^2(0,1)) \cap L^2(0, T; H^1_a(0,1))
$$
or
$$
X:= C(0, T; L^2_{1/a}(0,1)) \cap L^2(0, T; H^1_{1/a}(0,1)),
$$
and $L$ is the same constant of \eqref{ipo.probabilistica}.
\end{proposition}

\begin{theorem}\label{semilinear1}
 Assume that Hypotheses \ref{5.0}
and \ref{5.1} are satisfied. Then, for all $h \in L^2(Q_T)$, the
problem
\begin{equation}\label{sc}
\begin{gathered}
u_t -(au_x)_x + f(t,x,u, u_x)  = h(t,x), \quad (t,x) \in (0, T) \times (0,1), \\
u(t,1)= u(t,0)= 0, \quad t \in (0, T),\\
u(0,x)= u_0(x) \in H^1_a(0,1), \quad  x \in (0,1),
\end{gathered}
\end{equation}
has a solution $u \in H^1(0,T; L^2(0,1)) \cap L^2(0,T; H^2_a(0,1))$.
\end{theorem}

\begin{proof}
Let $X:= C(0, T; L^2(0,1)) \cap L^2(0, T; H^1_a(0,1))$ and, for any
$(t,x) \in (0, T) \times (0,1)$, set $b^{v}(t,x) := b(t,x,
       v(t,x))$ and $c^{v}(t,x) := c(t,x, v(t,x))$ for any $v \in X$.  Then, consider
the function
\[
\mathcal{T}: v \in X \mapsto u^{v} \in X, \] where $u^{v}$ is the
unique solution of
       \begin{equation}\label{linear-v}
       \begin{gathered}
u_t - \left(a(x)u_x \right) _x
 + b^{v}(t,x)u_x + c^{v}(t,x)u
 =h(t,x),\\
u(t,1) = u(t,0)=0,\\
u(0,x)=u_0(x) \in H^1_a(0,1).
       \end{gathered}
       \end{equation}
 By Theorem \ref{parabolic}, problem \eqref{linear-v}  has a unique weak solution $u \in X$.  Hence, if we
       prove that ${\mathcal{T}}$ has a fixed point $u^{
       v}$, i.e. ${\mathcal{T}}(u^{v})=u^{v}
       $, then $u^{v}$ is solution of \eqref{sc}.

To prove that
       $\mathcal{T}$ has a fixed point, by the Schauder's Theorem, it is sufficient to
       prove that
       \begin{enumerate}
       \item
       $\mathcal{T}: B_X \to B_X$,
       \item
       $\mathcal{T}$ is a compact function,
       \item
       $\mathcal{T}$ is a continuous function.
       \end{enumerate}
       Here $B_X:= \{v \in X: \|v\|_X \le R\}$,
$R:= C_T (\|u_0\|^2_{L^2(0,1)}+
       \|h\|^2_{L^2((0, T)\times (0,1))})$ ($C_T$ is the same constant of
       Theorem \ref{parabolic}) and
       $\|v\|_X:= \sup_{t \in
       [0, T]}\left(\|u(t)\|^2_{L^2(0,1)} \right)+ \int_0^{T}
       \|\sqrt{a}u_x\|^2_{L^2(0,1)}dt$.

The first item is a consequence of Theorem \ref{parabolic}. Indeed,
one has that $\mathcal{T}: X  \to B_X$ and in particular $\mathcal{T}:\;B_X \to B_X$. Moreover,
       it is easy to see that item (2)
is a simple consequence
       of the compactness Theorem
       \ref{compact_generale} below. This theorem is also useful for the
       proof of item (3). Indeed,
let $v_k \in X$ be such
       that $v_k \to v$ in $X$, as ${k \to +\infty}$. We
       want to prove that $u^k:=u^{v_k} \to u^{v}$ in
       $X$, as ${k \to +\infty}$. Here $u^k$ and
       $u^{v}$ are the solutions of \eqref{sc} associated to
       $v_k$ and $v$, respectively. By Remark \ref{remark},
       $u^k \in B_Y,$ where $Y:=   H^1(0,T; L^2(0,1)) \cap L^2(0,T; H^2_a(0,1))$.
       Proceeding as in Theorem \ref{compact1} below, one has that, up to
       subsequence, $u^k$ converges weakly to some $\bar{u}$ in $Y$
       and, thanks to Theorem \ref{compact_generale} below, strongly
       in $X$. Now, we prove that $\bar{u}=u^v$.
      Multiplying the equation
      \[
u^k_t - \left(a(x)u^k_x \right) _x
 + b^{v^k}(t,x)u^k_x + c^{v^k}(t,x)u^k
 =h(t,x),\\
      \]
      by a test
      function $\varphi \in H^1(0, T; L^2(0,1)) \cap L^2(0, T;
H^1_a(0,1))$ and integrating over $Q_T$ (we recall that
$Q_T:=(0,T)\times(0,1)$), we have:
\begin{align*}
&\int_0^1 u^k(T,x)\varphi(T,x) dx - \int_0^1 u_0(x) \varphi(0,x) dx
- \int_0^T \int_0^1 \varphi_t (t,x)u^k(t,x)\,dx\,dt\\
&= - \int_0^T\int_0^1 a(x)  u^k_x(t,x)\varphi_x(t,x)\,dx\,dt
 - \int_0^T\int_0^1b^{v^k}(t,x)u^k_x(t,x)\varphi(t,x) \,dx\,dt\\
&\quad -\int_0^T\int_0^1 c^{v^k}(t,x)u^k(t,x) \varphi(t,x) \,dx\,dt
+ \int_0^T\int_0^1 h(t,x)\varphi(t,x) \,dx\,dt.
\end{align*}
Since $u^k$ converges strongly to $\bar{u}$ in $X$, it is immediate
to prove
\begin{gather*}
\lim_{k\to+\infty}\int_0^1 u^k(T,x)\varphi(T,x)dx =
\int_0^1\bar{u}(T,x)\varphi(T,x) dx,
\\
\lim_{k\to+\infty}\int_0^T \int_0^1 \varphi_t
(t,x)u^k(t,x)\,dx\,dt=\int_0^T \int_0^1 \varphi_t (t,x)\bar{u}(t,x)\,dx\,dt
\\
\lim_{k\to+\infty}\int_0^T\int_0^1 a  u^k_x
\varphi_x\,dx\,dt=\int_0^T\int_0^1 a  \bar{u}_x \varphi_x\,dx\,dt
\end{gather*}
for all $\varphi \in H^1(0, T; L^2(0,1)) \cap L^2(0, T;
H^1_a(0,1))$. Moreover
\begin{equation}\label{ck}
\lim_{k\to+\infty}\int_0^T\int_0^1 c^{v^k}(t,x)u^k(t,x)
\varphi(t,x) \,dx\,dt=\int_0^T\int_0^1 c^{v}(t,x)\bar{u}(t,x)
\varphi(t,x) \,dx\,dt.
\end{equation}
Indeed,
\begin{align*}
&\big|\int_0^T\int_0^1(c^{v^k}(t,x)u^k(t,x)-c^{v}(t,x)\bar{u}(t,x))\varphi(t,x)
\,dx\,dt \big| \\
&\le \int_0^T\int_0^1 \big|
(u^k(t,x)-\bar{u}(t,x))c^{v^k}(t,x)\varphi(t,x)\big|\,dx\,dt\\
&\quad +\int_0^T\int_0^1 \big|c^{v^k}(t,x)-c^{v}(t,x)\big|
|\bar{u}(t,x)\varphi(t,x)|\,dx\,dt\\
&\le C\|u^k-\bar{u}\|_X\|\varphi\|_{L^2(0,1)}+\int_0^T\int_0^1
\big| c^{v^k}(t,x)-c^{v}(t,x)\big||\bar{u}(t,x)\varphi(t,x)|\,dx\,dt,
\end{align*}
where $C$ is the constant of \eqref{f2}. By the Lebesgue Theorem and
by Proposition \ref{prop1}, it follows
\[
\lim_{k\to+\infty}\int_0^T\int_0^1 \big|
c^{v^k}(t,x)-c^{v}(t,x)\big||\bar{u}(t,x)\varphi(t,x)|\,dx\,dt=0.
\]
Thus, since $u^k$ converges strongly to $\bar{u}$ in $X$, \eqref{ck}
holds.

Finally,
\begin{equation}\label{bk}
\lim_{k\to+\infty}\int_0^T\int_0^1b^{v^k}(t,x)u^k_x(t,x)
\varphi(t,x) \,dx\,dt=\int_0^T\int_0^1b^{v}(t,x)\bar{u}_x(t,x)
\varphi(t,x) \,dx\,dt.
\end{equation}
Indeed,
\begin{equation}\label{bk1}
\begin{aligned}
&\big|\int_0^T\int_0^1(b^{v^k}(t,x)u_x^k(t,x)-b^{v}(t,x)\bar{u}_x(t,x))\varphi(t,x)
\,dx\,dt \big| \\
&\le \int_0^T\int_0^1 \big|(u^k_x(t,x)-\bar{u}_x(t,x))b^{v^k}(t,x)\varphi(t,x)\big|\,dx\,dt\\
&\quad +\int_0^T\int_0^1 \big|b^{v^k}(t,x)-b^{v}(t,x)\big| |\bar{u}_x(t,x)\varphi(t,x)|\,dx\,dt\\
&\le L\|u^k-\bar{u}\|_X\|\varphi\|_{L^2(0,1)}\\
&\quad +\int_0^T\int_0^1 \frac{\left|b^{v^k}(t,x)-b^{v}(t,x)\right|}{\sqrt{a(x)}}
|\sqrt{a(x)}\bar{u}_x(t,x)\varphi(t,x)|\,dx\,dt,
\end{aligned}
\end{equation}
where $L$ is the constant of \eqref{ipo.probabilistica}. The
strongly convergence of $u^k$ to $\bar{u}$ in $X$, Proposition
\ref{prop1}, the Lebesgue Theorem and \eqref{bk1} imply that
\eqref{bk} holds. Thus $u^{ v}= \bar{u}$, that is $\bar{u}\in Y$  is the solution of
\eqref{sc} associated to $v$.

Hence, ${\mathcal{T}}$ has a fixed point  $u^{
       v}\in Y$ and, in particular, $u^{v}$ solves \eqref{sc}.
\end{proof}

To prove that the existence result holds also if the initial data
$u_0$ is in $L^2(0,1)$, we observe that \eqref{ipo.probabilistica}
implies for a.e. $(t,x) \in (0, T)\times (0,1)$ and for all $(u,p,q)
\in \mathbb{R}^3,$
\begin{equation}\label{f7}
 |f(t,x,u,p) - f(t,x, u,q)| \le L\sqrt{a(x)} |p-q|.
\end{equation}
Moreover, by \eqref{f2} and \eqref{f7}, it follows that for a.e.
$(t,x) \in (0, T) \times (0,1)$ and for all
$(u,v,p,q) \in \mathbb{R}^4$,
\begin{equation}\label{1'}
|(f(t,x,u,p) - f(t,x, v, q))(u-v)| \le M[|u-v|^2 + \sqrt{a(x)}
|p-q||u-v|],
\end{equation}
for some positive constant $M$.

Hence, using the previous estimates and Theorem \ref{semilinear1},
one can prove, as in \cite{cf}, the following result.

\begin{theorem}\label{semilinear2}
 Assume that Hypotheses \ref{5.0}
and \ref{5.1} are satisfied. Then, for all $h \in L^2(Q_T)$, the
problem
\begin{equation}\label{sc1}
\begin{gathered}
u_t -(au_x)_x + f(t,x,u, u_x)  = h(t,x), \quad (t,x) \in (0, T) \times (0,1), \\
u(t,0)= 0, \quad  t \in (0, T),\\
u(0,x)= u_0(x) \in L^2(0,1), \quad  x \in (0,1),
\end{gathered}
\end{equation}
has a solution $u \in C([0,T]; L^2(0,1)) \cap L^2 (0,T;H^1_a(0,1))$.
\end{theorem}

\begin{proof}
Let $(u^j_0)_j \in H^1_a(0,1)$ be such that $\lim_{j \to +
\infty}\|u_0^j - u_0\|_{L^2(0,1)} = 0$. Denote by $u^j$ the
solutions of \eqref{sc1} with respect to $u^j_0$. By Theorem
\ref{semilinear1}, $u^j \in H^1(0,T; L^2(0,1)) \cap L^2(0,T;
H^2_a(0,1))$. Then $(u^j)_j$ is a Cauchy sequence in $X:=C([0,T];
L^2(0,1)) \cap L^2 (0,T; H^1_a(0,1))$. In fact $u^j - u^i$ solves
the system
\begin{gather*}
(u^j - u^i)_t -(a(u^j - u^i)_x)_x + f(t,x,u^j, u^j_x)
- f(t,x, u^i, u^i_x)  = 0, \\
(u^j - u^i)(t,1)=  (u^j - u^i)(t,0) =0, \\
(u^j - u^i)(0,x)= (u^j_0 - u^i_0)(x),
\end{gather*}
where $(t,x) \in (0, T) \times (0,1) $. Multiplying
\[
(u^j - u^i)_t -(a(u^j - u^i)_x)_x + f(t,x,u^j, u^j_x)- f(t,x, u^i,
u^i_x)  = 0
\]
by $u^j- u^i$ and integrating over $(0,1)$, one has, using
\eqref{1'},
\[
\frac{1}{2}\frac{d}{dt}\int_0^1 |u^j - u^i|^2dx + \int_0^1a|u^j_x-
u^i_x|^2dx \le \int_0^1 M [ |u^j- u^i|^2 + \sqrt{a}|u^j_x -
u^i_x||u^j - u^i|]dx.
\]
Integrating over $(0, t)$:
\begin{align*}
&\frac{1}{2}\|(u^j - u^i)(t)\|^2_{L^2(0,1)} +
\int_0^t\int_0^1a|(u^j- u^i)_x|^2\,dx\,ds \\
&\le \frac{1}{2}\|u^j_0 - u^i_0\|^2_{L^2(0,1)}
 + M\int_0^t\int_0^1 |u^j- u^i|^2\,dx\,ds + \frac{\epsilon
M}{2}\int_0^t\int_0^1a|u^j_x-u^i_x|^2\,dx\,ds \\
&\quad + \frac{M}{2\epsilon}\int_0^t\int_0^1|u^j- u^i|^2\,dx\,ds.
\end{align*}
Thus
\begin{equation}\label{stima?}
\begin{aligned}
&\frac{1}{2}\|(u^j-u^i)(t)\|^2_{L^2(0,1)} + \big(1- \frac{\epsilon
M}{2} \big)\int_0^t\int_0^1a|u^j_x -u^i_x|^2 \,dx\,ds \\
&\le \frac{1}{2}\|u^j_0- u^i_0\|^2_{L^2(0,1)} +
M_{\epsilon}\int_0^t\int_0^1|u^j- u^i|^2\,dx\,ds,
\end{aligned}
\end{equation}
where $M_{\epsilon}:=\frac{M(1+2\epsilon)}{2\epsilon}$. By
Gronwall's Lemma
\[
\|(u^j- u^i)(t)\|^2_{L^2(0,1)} \le e^{M_{\epsilon}t}\|u^j_0 -
u^i_0\|^2_{L^2(0,1)},
\]
and
\begin{equation}\label{stima??}
\sup_{t \in [0, T]}\|(u^j- u^i)(t)\|^2_{L^2(0,1)} \le
e^{M_{\epsilon}T}\|u^j_0 - u^i_0\|^2_{L^2(0,1)}.
\end{equation}
This implies that $(u^j)_j$ is a Cauchy sequence in $C(0, T;L^2(0,1))$.
Moreover, by \eqref{stima?}, one has
\[
\big(1- \frac{\epsilon M}{2} \big)\int_0^t\int_0^1a|u^j_x
-u^i_x|^2\,dx\,ds \le \frac{1}{2}\|u^j_0- u^i_0\|^2_{L^2(0,1)} +
M_{\epsilon}\int_0^t\int_0^1|u^j- u^i|^2\,dx\,ds.
\]
Using \eqref{stima??}, it follows
\begin{align*}
\int_0^t\|\sqrt{a}(u^j_x -u^i_x)\|^2_{L^2(0,1)}ds & \le M_{\epsilon,
T}(\|u^j_0- u^i_0\|^2_{L^2(0,1)}
+ \sup_{t \in [0, T]}\|u^j- u^i\|^2_{L^2(0,1)} )\\
&\le M_{\epsilon, T}\|u^j_0- u^i_0\|^2_{L^2(0,1)}.
\end{align*}
Thus $(u^j)_j$ is a Cauchy sequence also in $L^2(0, T; H^1_a(0,1))$.
Then there exists $\bar{u} \in X$ such that
\[
\lim_{j \to + \infty}\|u^j - \bar{u}\|_X= 0.
\]
Proceeding as in the proof of Theorem \ref{semilinear1}, one can
prove that $\bar{u}$ is a solution of \eqref{sc1}.
\end{proof}

Analogously, recalling Definition \ref{def_nondiv}, one can prove
the following results:
\begin{theorem}\label{semilinear3} Assume that Hypotheses \ref{5.0}
and \ref{5.1} are satisfied. Then, for all $h \in
L^2_{1/a}(Q_T)$, the problem
\begin{gather*}
u_t -au_{xx} + f(t,x,u, u_x)  = h(t,x), \quad (t,x) \in (0, T) \times (0,1), \\
u(t,0)= 0, \quad  t \in (0, T),\\
u(0,x)= u_0(x) \in H^1_{1/a}(0,1), \quad  x \in (0,1),
\end{gather*}
has a solution $u \in H^1(0,T; L^2_{1/a}(0,1)) \cap L^2(0,T;
H^2_{1/a}(0,1))$.
\end{theorem}

\begin{theorem}\label{semilinear4}
Assume that Hypotheses \ref{5.0}
and \ref{5.1} are satisfied. Then, for all $h \in
L^2_{1/a}(Q_T)$, the problem
\begin{gather*}
u_t -au_{xx} + f(t,x,u, u_x)  = h(t,x), \quad (t,x) \in (0, T) \times (0,1), \\
u(t,0)= 0, \quad t \in (0, T),\\
u(0,x)= u_0(x) \in L^2_{1/a}(0,1), \quad x \in (0,1),
\end{gather*}
has a solution $u \in C([0,T]; L^2_{1/a}(0,1)) \cap L^2
(0,T; H^1_{1/a}(0,1))$.
\end{theorem}

We recall that Theorems \ref{semilinear3} and \ref{semilinear4}
still hold if we substitute \eqref{ipo.probabilistica} with
\eqref{ipo}. In this case \eqref{f7} becomes for a.e.
$(t,x) \in (0,T)\times (0,1)$ and for all $(u,p,q) \in \mathbb{R}^3$,
\[
|f(t,x,u,p) - f(t,x,u, q)| \le L|p-q|.
\]

\section{Appendix}
In this section we will give some compactness theorems that we have
used in the previous section.

\begin{theorem}\label{compact}
Assume that the function $a$ satisfies Hypothesis \ref{5.0}. Then, the space $H^1_a(0,1)$ is
compactly imbedded in $L^2(0,1)$.
\end{theorem}

\begin{proof}
Clearly, $H^1_a(0,1)$ is continuously imbedded in $L^{2}(0,1)$.
Now, let $\epsilon >0$ and $M>0$. We want to prove that there exists
$\delta >0$ such that for all $u \in H^{1}_{a}(0,1)$ with $\|u\|_{H^1_a(0,1)}^2 \le M$ and
for all $|h| < \delta$ it results
\begin{gather}\label{comp1}
\int_{\delta}^{1-\delta}|u(x+h) - u(x)|^{2}dx <\epsilon,\\
\label{comp2} \int_{1-\delta}^{1}|u(x)|^{2}dx
+\int_{0}^{\delta}|u(x)|^{2}dx < \epsilon.
\end{gather}
Hence, let $u \in H^{1}_{a}(0,1)$ such that
 $\|u\|_{H^1_a(0,1)}^2 \le M$ and take $\delta >0$ such that
\begin{equation}\label{g}
\delta < g(\epsilon):=\min \Big\{\frac{x_0}{2}, \frac{1+x_0}{2},
\sqrt{\frac{\epsilon}{2M \max_{[0, \frac{x_0}{2}]}\frac{1}{a}}},
\sqrt{\frac{\epsilon}{2M \max_{[\frac{1+x_0}{2},
1]}\frac{1}{a}}}\Big\}.
\end{equation}
Then,
\begin{align*}
\int_{0}^{\delta}|u(x)|^{2}dx
&\le \int_{0}^{\delta}\Big|\int_0^x\frac{1}{\sqrt{a(y)}}\sqrt{a(y)}u'(y)dy
\Big|^2 dx
\le M \int_0^\delta \delta\max_{[0, \frac{x_0}{2}]}\frac{1}{a}dx \\
&< M \delta^2\max_{[0,\frac{x_0}{2}]}\frac{1}{a}<\frac{\epsilon}{2}.
\end{align*}
Analogously,
\[
\int_{1-\delta}^{1}|u(x)|^{2}dx < \frac{\epsilon}{2}.
\]
Now, let $h$ be such that $|h| < \delta$ and, for simplicity, assume
$h>0$ (the case $h<0$ can be treated in the same way). Then
\[
|u(x+h) - u(x)|^{2}\le
\|u\|^{2}_{H^1_a(0,1)}\int_{x}^{x+h}\frac{dy}{a(y)}.
\]
Integrating over $(\delta, 1-\delta)$, it results
\begin{align*}
\int_{\delta}^{1-\delta}|u(x+h) - u(x)|^{2}dx
& \le \|u\|^{2}_{H^1_a(0,1)}\int_{\delta}^{1-\delta}dx\int_{x}^{x+\delta}
\frac{dy}{a(y)}\\
&\le M\int_{\delta}^{1}\frac{dy}{a(y)}\int_{y-\delta}^{y}dx\le
M\delta \left\| \frac{1}{a}\right\|_{L^1(0,1)}.
\end{align*}
Moreover, since $\lim_{\delta \to 0}M\delta \|\frac{1}{a}\|_{L^1(0,1)} =0$,
there exists $\eta (\epsilon)>0$ such that if $\delta < \eta(\epsilon)$, then
$M\delta \left\| \frac{1}{a}\right\|_{L^1(0,1)} <\epsilon$. Thus, taking
$\delta < \min\{g(\epsilon), \eta(\epsilon)\}$, \eqref{comp1} and \eqref{comp2} are verified and
the thesis follows (see, e.g., \cite[Chapter IV]{b}).
\end{proof}


Using Theorem \ref{compact} one can prove the next theorem.

\begin{theorem}\label{compact1}
Assume that the function $a$ satisfies Hypothesis \ref{5.0}. Then,
the space $H^2_a(0,1)$ is compactly imbedded in $H^1_a(0,1)$.
\end{theorem}

\begin{proof}
Take $(u_n)_n \in \overline{B}_{H^2_a(0,1)}$. Here $B_{H^2_a(0,1)}$
denotes the unit ball of $H^2_a(0,1)$. Since $H^2_a(0,1)$ is
reflexive, then, up to subsequence, there exists $u \in H^2_a(0,1)$
such that $u_n$ converges weakly to $u$ in $H^2_a(0,1)$. In
particular, $u_n$ converges weakly to $u$ in $H^1_a(0,1)$ and in
$L^2(0,1)$. But, since by the previous Theorem $H^1_a(0,1)$ is
compactly imbedded in $L^2(0,1)$, then, up to subsequence, there
exists $v \in L^2(0,1)$ such that $u_n$ converges strongly to $v$ in
$L^2(0,1)$. Thus $u_n$ converges weakly to $v$ in $L^2(0,1)$. By
uniqueness $v\equiv u$. Then we can conclude that the sequence $u_n$
converges strongly to $u$ in $L^2(0,1)$.

Moreover, one has
\[
\|\sqrt{a}u_{n}' - \sqrt{a}u' \|_{L^2(0,1)} \to 0, \quad \text{as }
n \to + \infty.
\]
Indeed, using the H\"older inequality, one has
\begin{align*}
\|\sqrt{a}(u_{n} - u)'\|^2_{L^2(0,1)}
&= \int_0^1 a(u_{n} - u)'(u_{n} - u)'dx \\
&=- \int_0^1 (a(u_{n} - u)')'(u_n-u)dx \\
&\le \|(a(u_n-u)')'\|_{L^2(0,1)}\|u_n -u\|_{L^2(0,1)} \to 0,
\end{align*}
as $ n \to + \infty$. Hence the sequence $u_n$ converges strongly to
$u$  in $H^1_a(0,1)$.
\end{proof}

For the proof of Theorem \ref{compact_generale} below we will use
Aubin's Theorem, that we give here for the reader's convenience.

\begin{theorem}[{\cite[Chapter 5]{an} or \cite[Theorem 5.1]{l}}]\label{compact_th}
Let $X_0$, $X_1$ and $X_2$ be three Banach spaces such that $X_0
\subset X_1 \subset X_2$, $X_0$, $ X_2$ are reflexive spaces and the
injection of $X_0$ into $X_1$ is compact. Let $r_0$, $r_1 \in (1,
+\infty)$ and $a,b \in \mathbb{R}$, $a<b$. Then the space
\[
L^{r_0}(a,b; X_0) \cap W^{1,r_1}(a,b; X_2)
\]
is compactly imbedded in $L^{r_0}(a,b; X_1)$.
\end{theorem}

Now we are ready to prove the next compactness Theorem.

\begin{theorem}\label{compact_generale} Assume that function $a$
satisfies Hypothesis \ref{5.0}. Then, the space
 $H^1(0, T;L^2(0,1)) \cap L^2(0, T; H^2_a(0,1))$ is
compactly imbedded in $L^2(0, T; H^1_a(0,1))$ $ \cap$ $ C(0,T;
L^2(0,1))$.
\end{theorem}

\begin{proof}
Using Aubin's Theorem \ref{compact_th} with $r_0= r_1=2$, $X_0 =
H^2_a(0,1)$, $X_1 = H^1_a(0,1)$, $X_2 = L^2(0,1)$, $a=0$ and $b=T$,
one has
\[
H^1(0, T;L^2(0,1)) \cap L^2(0, T; H^2_a(0,1))
\hookrightarrow\hookrightarrow L^2(0, T; H^1_a(0,1)).
\]
Moreover, since $H^1(0, T;L^2(0,1)) \cap L^2(0, T;
H^2_a(0,1)))\hookrightarrow H^1(0, T; L^2(0,1))$ and $H^1(0, T;
L^2(0,1))\hookrightarrow\hookrightarrow C(0, T; L^2(0,1) )$, then
\[
H^1(0, T;L^2(0,1)) \cap L^2(0, T; H^2_a(0,1))\hookrightarrow
\hookrightarrow C(0, T; L^2(0,1) ).
\]
Thus
 \[
 H^1(0, T;L^2(0,1)) \cap L^2(0, T;
H^2_a(0,1))\hookrightarrow\hookrightarrow L^2(0, T; H^1_a(0,1)) \cap
C(0, T;L^2(0,1) ).
\]
\end{proof}

Analogously to Theorem \ref{compact_generale}, one can obtain the
following compactness result.

\begin{theorem}\label{compact_generale1}
Assume that the function $a$ satisfies Hypothesis \ref{5.0}. Then,
the space $H^1(0, T;L^2_{1/a}(0,1)) \cap L^2(0, T;H^2_{1/a}(0,1))$
is compactly imbedded in the space
$L^2(0, T;H^1_{1/a}(0,1))\cap C(0,T; L^2_{1/a}(0,1))$.
\end{theorem}

To prove Theorem \ref{compact_generale1}, it is sufficient to prove
the analogous results of Theorems \ref{compact} and \ref{compact1}.
In particular, we have:
\begin{theorem}\label{compact2} Assume that the function $a$
satisfies Hypothesis \ref{5.0}. Then, the space
$H^1_{1/a}(0,1)$ is compactly imbedded in
$L^2_{1/a}(0,1)$.
\end{theorem}

\begin{proof}
Clearly, $H^1_{1/a}(0,1)$ is continuously imbedded in
$L^{2}_{1/a}(0,1)$.
Now, let $\epsilon >0$ and $M>0$. We want to prove that there exists
$\delta >0$ such that for all $u \in H^{1}_{1/a}(0,1)$ with $\|u\|_{H^1_{1/a}(0,1)}^2 \le M$
and for all $|h| < \delta$
it holds that
\begin{gather}\label{comp1bis}
\int_{\delta}^{1-\delta}\frac{|u(x+h) - u(x)|^{2}}{a(x)}dx <\epsilon,\\
\label{comp2bis} \int_{1-\delta}^{1}\frac{|u(x)|^{2}}{a(x)}dx
+\int_{0}^{\delta}\frac{|u(x)|^{2}}{a(x)}dx < \epsilon.
\end{gather}
Hence, let $u \in H^{1}_{1/a}(0,1)$ such that
$\|u\|_{H^1_{1/a}(0,1)}^2 \le M$ and take $\delta >0$ such
that $ \delta < g(\epsilon)$, where $g(\epsilon)$ is defined in
\eqref{g}. Then,
\[
\int_{0}^{\delta}\frac{|u(x)|^{2}}{a(x)}dx \le \max_{[0,
\frac{x_0}{2}]}\frac{1}{a} \int_{0}^{\delta}\Big|\int_0^xu'(y)dy
\Big|^2 dx< M \delta^2\max_{[0,
\frac{x_0}{2}]}\frac{1}{a}<\frac{\epsilon}{2}.
\]
Analogously,
\[
\int_{1-\delta}^{1}\frac{|u(x)|^{2}}{a(x)}dx < \frac{\epsilon}{2}.
\]
Now, let $h$ be such that $|h| < \delta$ and, for simplicity, assume
$h>0$ (the case $h<0$ can be treated in the same way). Then
\[
\int_{\delta}^{1-\delta}\frac{|u(x+h) - u(x)|^{2}}{a(x)}dx =
\int_{\delta}^{1-\delta} \frac{1}{a(x)}\\Big|\int_x^{x+h}u'(y)dy
\Big|^2dx\le M\delta \| \frac{1}{a}\|_{L^1(0,1)}.
\]
Moreover, since $\lim_{\delta \to 0}M\delta \|\frac{1}{a}\|_{L^1(0,1)} =0$, there exists
$\eta (\epsilon) >0$ such that if $\delta < \eta(\epsilon)$, then
$M\delta \left\| \frac{1}{a}\right\|_{L^1(0,1)} <\epsilon$. Thus, taking $\delta < \min\{g(\epsilon),
\eta(\epsilon)\}$, \eqref{comp1bis} and \eqref{comp2bis} are
verified and the thesis follows (see, e.g., \cite[Chapter IV]{b}).
\end{proof}

Using Theorem \ref{compact2} one can prove the next theorem.

\begin{theorem}\label{compact3} Assume that the function $a$
satisfies Hypothesis \ref{5.0}. Then, the space
$H^2_{1/a}(0,1)$ is compactly imbedded in
$H^1_{1/a}(0,1)$.
\end{theorem}

\begin{proof}
Take $(u_n)_n \in \overline{B}_{H^2_{1/a}(0,1)}$. Here
$B_{H^2_{1/a}(0,1)}$ denotes the unit ball of
$H^2_{1/a}(0,1)$. As before, we can prove that, up to
subsequence, there exists $u \in H^2_{1/a}(0,1)$ such that
$u_n$ converges strongly to $u$ in $L^2_{1/a}(0,1)$.

Moreover, one has
\[
\|u_{n}' - u' \|_{L^2(0,1)} \to 0, \quad \text{as } n \to + \infty.
\]
Indeed, using the H\"older inequality, one has
\begin{align*}
\|(u_{n} - u)'\|^2_{L^2(0,1)}
&= \int_0^1 (u_{n} - u)'(u_{n} - u)'dx \\
&=- \int_0^1 (u_{n} - u)''(u_n-u)dx =- \int_0^1 \sqrt{a}(u_n-u)''\frac{u_n-u}{\sqrt{a}}dx\\
&\le \|a(u_n-u)''\|_{L^2(0,1)}\|u_n -u\|_{L^2_{1/a}(0,1)}
\to 0,
\end{align*}
as $ n \to + \infty$. Hence the sequence $u_n$ converges strongly to
$u$  in $H^1_{1/a}(0,1)$.
\end{proof}

\begin{thebibliography}{99}

\bibitem{acf}
   F. Alabau-Boussouira, P. Cannarsa, G. Fragnelli;
\emph{Carleman   estimates for degenerate parabolic operators with applications to
   null controllability}, J. Evol. Equ. \textbf{6} (2006), 161-204.

\bibitem{an} S. Ani\c{t}a;
\emph{Analysis and Control of Age-Dependent Population Dynamics},
Mathematical Modelling: Theory and Applications, 11. Kluwer Academic
Publishers, Dordrecht, 2000.

\bibitem{abhn} W. Arendt, C. J. K. Batty, M. Hieber, F. Neubrander;
\emph{ Vector-valued Laplace Transforms and Cauchy Problems}, Monographs in
Mathematics \textbf{96} (2001), Birkh\"{a}user Verlag, Basel.

\bibitem{bfr} V. Barbu, A. Favini, S. Romanelli;
\emph{Degenerate evolution equations and regularity of their associated semigroups}, Funkcial.
Ekvac. \textbf{39 } (1996), 421--448.

\bibitem{b} H. Brezis;
\emph{Analyse fonctionnelle}, Collection of
Applied Mathematics for the Master's Degree, Masson, Paris, 1983.

\bibitem{br} J. M. Buchot, J. P. Raymond;
\emph{A linearized model for boundary
layer equations}, in "Optimal control of complex structures"
(Oberwolfach, 2000), Internat. Ser. Numer. Math. \textbf{139}
(2002), Birkh\"{a}user, Basel, 31--42.

\bibitem{cmp} M. Campiti, G. Metafune, D. Pallara;
\emph{Degenerate self-adjoint
evolution equations on the unit interval}, Semigroup Forum
\textbf{57} (1998), 1--36.

\bibitem{cf} P. Cannarsa, G. Fragnelli;
\emph{Null controllability of semilinear degenerate parabolic equations in bounded domains},
 Electron. J. Differential Equations \textbf{2006} (2006), 1--20.

\bibitem{cfr}  P. Cannarsa, G. Fragnelli, D. Rocchetti;
\emph{Controllability results for a class of one-dimensional degenerate parabolic
          problems in nondivergence form}, J. Evol. Equ. \textbf{8} (2008), 583--616.

\bibitem{cmv} P. Cannarsa, P. Martinez, J. Vancostenoble;
\emph{Null controllability of the degenerate heat equations}, Adv. Differential
Equations \textbf{10} (2005), 153--190.

\bibitem{crv} P. Cannarsa, D. Rocchetti, J. Vancostenoble;
\emph{Generator of analytic semi-groups in $L^2$ for a
class of second order degenerate elliptic operators},  Control
Cybernet. \textbf{37} (2008), 831--878.

\bibitem{egg} H. Emamirad, G. R. Goldstein, J. A. Goldstein;
\emph{Chaotic solution for the Black-Scholes equation}, Proc. Amer. Math. Soc. \textbf{140}
(2012), 2043--2052.

\bibitem{en} K.J. Engel, R. Nagel;
\emph{One-Parameter Semigroups for  Linear Evolution Equations}, Springer-Verlag,
 New York, Berlin, Heidelberg, 1999.

\bibitem{fggr} A. Favini, G. R. Goldstein, J. A. Goldstein, S. Romanelli;
\emph{Degenerate Second Order Differential Operators Generating Analytic Semigroups in $L^p$
and $W^{1,p}$}, Math. Nachr. \textbf{238} (2002), 78--102.

\bibitem{F1} W. Feller;
\emph{The parabolic differential equations and the
associated semigroups of transformations}, Ann. of Math. \textbf{55}
(1952), 468--519.

\bibitem{F2} W. Feller;
 \emph{Diffusion processes in one dimension}, Trans. Am. Math. Soc. \textbf{97} (1954),
1--31.

\bibitem{f} G. Fragnelli;
\emph{Null controllability of degenerate parabolic
equations in non divergence form via Carleman estimates}, to appear
in Discrete Contin. Dyn. Syst. -
                         S.
\bibitem{g} J. A. Goldstein;
\emph{Semigroups of Linear Operators and Applications}, Oxford Univ.
Press, Oxford, New York, 1985.

\bibitem{gl} J. A. Goldstein, A. Y. Lin;
\emph{An $L^ p$ semigroup approach to degenerate parabolic boundary value
problems}, Ann. Mat. Pura Appl. \textbf{159} (1991), 211--227.

\bibitem{gmr} J. A. Goldstein, R. M. Mininni, S. Romanelli;
\emph{A new explicit formula for the solution of the Black-Merton-Scholes equation}, Infinite
dimensional stochastic analysis. Quantum Probab. White Noise Anal.,
World Sci. Publ., Hackensack, NJ, \textbf{22} (2008), 226--235.

\bibitem{l} J. L. Lions;
\emph{Quelques m\'{e}thodes de r\'{e}solution des probl\`{e}mes aux limites non lin\'{e}aires},
 Dunod-Gauthier-Villars, 1969.

 \bibitem{mp} G. Metafune, D. Pallara;
\emph{Trace formulas for some singular differential operators and applications},
Math. Nachr. \textbf{211} (2000), 127--157.

\bibitem{mro} R. M. Mininni, S. Romanelli;
\emph{Martingale estimating functions for Feller diffusion processes
generated by degenerate elliptic operators}, J. Concr. Appl. Math.
\textbf{1} (2003), 191--216.

\bibitem{mr} D. Mugnolo, A. Rhandi;
\emph{On the domain of a Fleming-Viot type operator on an $L^p$ space with invariant
measure}, Note di Matematica \textbf{31} (2011), 139--148.

\bibitem{s2} A. Stahel;
\emph{Existence and uniqueness of solutions to degenerate semilinear parabolic
equations}, J. Differential Equations \textbf{93} (1991), 311-331.

\bibitem{s1} A. Stahel; \emph{Degenerate semilinear parabolic
equations}, Differential Integral Equations \textbf{5} (1992),
683--691.
 \bibitem{s} N. Shimakura;
\emph{Partial Differential Operators of
   elliptic type}, Translations of Mathematical Monographs \textbf{99} (1992),
   American Mathematical Society, Providence, RI.

\bibitem{t1} K. Taira;
\emph{On the existence of Feller semigroups with boundary conditions,
I}, Mem. Amer. Math. Soc.\textbf{ 99 }(1992), 8--73.

\bibitem{t2} K. Taira;
\emph{On the existence of Feller semigroups with boundary conditions,
II}, J. Funct. Anal. \textbf{129} (1995), 108--131.

\bibitem{frt} K. Taira, A. Favini, S. Romanelli;
\emph{Feller semigroups and degenerate elliptic operators with Wentzell boundary
conditions}, Studia Mathematica \textbf{145} (2001), 17--53.

\bibitem{v} V. Vespri;
\emph{Analytic semigroups, degenerate elliptic operators and applications to nonlinear Cauchy
problems}, Ann. Mat. Pura Appl. (IV) CLV (1989), 353--388.

 \end{thebibliography}

\end{document}
