\section{Introduction}

We consider the stochastic reaction-diffusion equation
\begin{align}\label{rd}
\frac{\partial X_t(\xi)}{\partial t}
=\Delta X_t(\xi) +f(X_t(\xi)) + G \frac{\partial W_t(\xi)}{\partial t},
\quad (t, \xi) \in \mathbb R_+ \times \mathscr O.
\end{align}
Here the homogeneous Dirichlet boundary condition on the bounded, open subset $\mathscr O$ of $\mathbb R^d$ is considered, the initial value $X_0=x$ vanishes on the boundary $\partial \mathscr O$ of $\mathscr O$, the nonlinear drift function $f$ has polynomial growth and satisfies a certain dissipative condition, the driving process $(W_t)_{t \ge 0}$ is a $U$-valued cylindrical Wiener process, and the linear operator $G$ is a densely defined closed linear operator on $U$ which could be unbounded (see Section \ref{sec2} for more details).

Let $x \in L^p=L^p(\mathscr O)$, $p \ge 2$, and denote by $(X_t^x)_{t \ge 0}$ a mild solution of Eq. \eqref{rd} with the initial datum $X_0=x$.
Then $(X^x_t)_{t \ge 0}$ is a Markov process and it generates a Markov semigroup $(P_t)_{t \ge 0}$ defined as
\begin{align} \label{pt}
P_t \phi(x):=\mathbb E \phi(X^x_t), \quad
t \ge 0, \ x \in L^p, \ \phi \in \mathcal B_b(L^p).
\end{align}
Our main concern in this paper is to investigate Harnack inequalities and ergodicity of the Markov semigroup $(P_t)_{t \ge 0}$ defined in \eqref{pt} in the Banach space, $L^p$, for all $p \ge 2$.

The stochastic reaction-diffusion equation \eqref{rd} has numerous applications in material sciences and chemical kinetics \cite{ET89}.
When $f(\xi)=\xi-\xi^3$, $\xi \in \mathbb R$, Eq. \eqref{rd} is also called the stochastic Allen--Cahn equation or the stochastic Ginzburg--Landau equation.
It is widely used in many fields, for example, the random interface models and stochastic mean curvature flow \cite{Fun16}.
There are many interesting and important properties for the solution of Eq. \eqref{rd}, which have been investigated in Hilbert settings.
For example, the existence of invariant measures and ergodicity are studied in \cite{BS20, Hai02(PTRF), Kaw05(PA), WXX17(BER)},
the large deviation principles are investigated in \cite{BBP17(ARMA), WRD12(JDE), XZ18(BER)},
and sharp interface limits are derived in \cite{KORV07(CPAM), Web10(CPAM),  Yip98(ARMA)}.
See also \cite{BGJK17, CH19(SINUM), LQ20(IMA), LQ20(SPDE)} and references therein for analysis in the numerical aspect.

In contrast to SPDEs in Hilbert spaces, only a few papers are treating the regularity, such as invariant measures and the ergodicity, of the Markov semigroup $(P_t)_{t \ge 0}$ for SPDEs even with Lipschitz coefficients in Banach spaces.
The authors in \cite{BLS10(JEE)} studied invariant measures for SPDEs in martingale-type (M-type) 2 Banach spaces, under Lipschitz and dissipativity conditions, driven by regular noise.
For white-noise driven stochastic heat equation (Eq. \eqref{see} with Lipschitz coefficients), \cite{BR16(DCDS)} showed the uniqueness of the invariant measure, if it exists, on $L^p(0, 1)$ with $p>4$; the case for $p \in (2,4]$ is unknown.
Recently, their method was extended in \cite{BK18(FS)} to an SPDE, arisen in stochastic finance, in a weighted $L^p$-space.
See also \cite{Nee01(OTAA)} for the uniqueness of the invariant measure, if it exists, of the Ornstein--Uhlenbeck process (with form \eqref{df-wa}) using a pure analytical method.

For SPDEs with non-Lipschitz coefficients, \cite{BG99(SPA)} obtained the existence of an invariant measure for Eq. \eqref{rd} in the space of continuous functions under the martingale solution framework; the uniqueness of the invariant measure was derived in \cite{Cer03(PTRF), Cer05(PTRF)} by taking advantage of the fact that a polynomial is uniformly continuous on bounded subsets of continuous functions.
Recently, \cite{KN13(SPDE)} showed the existence of a unique invariant measure on the space of continuous complex functions for the stochastic complex Ginzburg--Landau (Eq. \eqref{rd} with $f(u)=- \mathbf i |u|^2 u$, where $\mathbf i=\sqrt{-1}$), relying on some good estimates of the solution in the Hilbert--Sobolev spaces $\dot H^\beta$ with $\beta>d/2$, so that the noise is spatially regular enough.

To show the uniqueness of the invariant measure, these authors mainly formulated a Bismut--Elworthy--Li formula for the derivative $D P_t$ to get a gradient estimate, which shows the strong Feller property of $P_t$.
Then the uniqueness of the invariant measure follows immediately by Khas'minskii and Doob theorems, provided the irreducibility holds.
The difficulties for the study of the uniqueness of an invariant measure for SPDEs in Banach settings arise, mainly because the tools frequently used in the Hilbert space framework cannot be extended in a straightforward way to the Banach space settings \cite{BR16(DCDS)}.

In the past decade, Wang-type dimension-free inequalities have been a new and efficient tool to study diffusion semigroups.
They were first introduced in \cite{Wan97(PTRF)} for elliptic diffusion semigroups on non-compact Riemannian manifolds and in \cite{Wan10(JMPA)} for heat semigroups on manifolds with boundary.
Roughly speaking, such inequality for the Markov semigroup $(P_t)_{t \ge 0}$ in a Banach space $E$ is formulated as
\begin{align} \label{har}
\Phi(P_t \phi(x)) & \le P_t (\Phi(\phi)(y)) \exp \Psi(t, x, y),
\quad t>0, \ x, y \in E, \ \phi \in \mathcal B^+_b(E),
\end{align}
where $\Phi: [0, \infty) \rightarrow [0, \infty)$ is convex, $\Psi$ is non-negative on $[0, \infty) \times E \times E$ with $\Psi(t, x, x)=0$ for all $t>0$ and $x \in E$, and $\mathcal B^+_b(E)$ denotes the family of all measurable and bounded, non-negative functions on $E$.

There are two frequently used choices of $\Phi$.
One is given by a power function $\Phi(\xi)=\xi^\mathbf s$, $\xi \ge 0$, for some $\mathbf s>1$,
then \eqref{har} is reduces to
\begin{align} \label{har-1}
(P_t \phi(x))^\mathbf s & \le P_t \phi^\mathbf s(y)  \exp \Psi(t, x, y),
\quad t>0, \ x, y \in E, \ \phi \in \mathcal B^+_b(E).
\end{align}
Another is given by $\Phi(\xi)=e^\xi$, $\xi \in \mathbb R$, in which one may use $\log \phi$ to replace $\phi$, so that \eqref{har} becomes
\begin{align} \label{har-2}
P_t \log \phi(x) & \le \log P_t \phi(y)+\Psi(t, x, y),
\quad t>0, \ x, y \in E, \ \phi \in \mathcal B^+_b(E).
\end{align}

The above inequalities \eqref{har-1} and \eqref{har-2} are called power-Harnack inequality and log-Harnack inequality, respectively.
Both inequalities have been investigated extensively and applied to SODEs and SPDEs via coupling by the change of measure, see, e.g., \cite{Kaw05(PA), Wan07(AOP), WZ13(JMPA), WZ14(SPA), Zha10(PA)}, the monograph \cite{Wan13}, and references therein.
Besides the gradient estimate which yields strong Feller property, these Harnack inequalities also have a lot of other applications.
For example, they are used to study the contractivity of the Markov semigroup $(P_t)_{t \ge 0}$  in \cite{DRW09(JFA), Wan11(AOP), Wan17(JFA)} and to derive almost surely strictly positivity of the solution for an SPDE in \cite{Wan18(PTRF)}.

For SPDEs with polynomial growth drift driven by rough noise in Hilbert settings, we are only aware \cite{Kaw05(PA), Xie19(JDE)} investigating Harnack inequalities in a weighted $L^2$-space and a subspace of $L^2$ consisting of all non-negative functions, respectively.
When the noise is of trace-class, i.e., $G$ appearing in Eq. \eqref{rd} is a Hilbert--Schmidt operator, then the variational solution theory can be used and multiplicative noise can also be considered, as the solution is a semi-martingale so that It\^o formula can be applied; see, e.g., \cite{HLL20(SPL), Liu09(JEE)}.
In the rougher white noise case, i.e., $G$ coincides with the identity operator in $L^2$, the variational solution would not exist; one needs to adopt the mild solution theory instead.
Generally, the mild solution is not a semi-martingale and thus It\^o formula is not available.
We also note that to derive Harnack inequalities for an SPDE with white noise, \cite{WZ14(SPA)} used finite-dimensional approximations to get a sequence of SODEs such that the arguments developed in \cite{Wan11(AOP)} for SODEs can be applied.

\iffalse
As it is known that a polynomial is no longer continuous on bounded subsets in $(L^p)_{p \ge 2}$-spaces.
Also, to handle the $(L^p)_{p \ge 2}$ case, the finite-dimensional approximations argument seems to be impossible.

To overcome the difficulty of polynomial growth, our main idea is to use the negative high order term in the dissipative condition \eqref{f},  to control other terms.
In particular, when $q>2$, i.e., the drift function $f$ has a super-linear growth, we derive a (\eqref{est-tig1} in Theorem \ref{tm-erg}) and a uniform exponential moments' estimation $\sup_{t \ge 0} \mathbb E \exp(\epsilon \|\cdot\|_p^p)$ with some positive constant $\epsilon$ for the solution (see \eqref{exp} in Lemma \ref{lm-exp}), which will be useful to establish the existence of an invariant measure and to derive the Gaussian concentration property for the invariant measure, if it exists, and contractivity property of $(P_t)_{t \ge 0}$.
\fi

Our main idea to derive the Harnack inequalities \eqref{har-log} and \eqref{har-pow} in the first main result, Theorem \ref{main1}, for the Markov semigroup $(P_t)_{t \ge 0}$ defined by \eqref{pt} of Eq. \eqref{rd} on $(L^p)_{p \ge 2}$-spaces, under the dissipativity condition \eqref{f} and the polynomial growth condition \eqref{f+}, are the construction of a coupling (see \eqref{rd-y}) of the change of measure and a uniform pathwise estimate (see \eqref{x-y}) for this coupling.
As by-products, a gradient estimate and the uniqueness of the invariant measure for $(P_t)_{t \ge 0}$, if it exists, follow immediately.
To the best of our knowledge, the Harnack inequalities \eqref{har-log} and \eqref{har-pow} are the first two Harnack inequalities for SPDEs in Banach settings.

To show the existence of an invariant measure for Eq. \eqref{rd} with super-linear growth and without strong dissipativity, e.g., $q>2$ and $\lambda$ defined in \eqref{lam} is non-positive, we derive a uniform estimate of $\mu_n(\|\cdot\|_{q+p-2}^{q+p-2})$ (see \eqref{est-tig1}) for a sequence of probability measures $(\mu_n)_{n \in \mathbb N_+}$ (defined in \eqref{mun}).
However, this is not strong enough to conclude the tightness of $(\mu_n)_{n \in \mathbb N_+}$, since the embedding $L^{q+p-2} \subset L^p$ is not compact.
To overcome this difficulty, we utilize a compact embedding (see \eqref{emb}) by a Sobolev--Slobodeckii space.
This forces us to derive a uniform estimate  of $\mu_n(\|\cdot\|_{\beta, p})$ (see \eqref{est-tig2}, with the $\|\cdot\|_{\beta, p}$-norm defined in \eqref{sob-slo}), where the aforementioned uniform estimate of $\mu_n(\|\cdot\|_{q+p-2}^{q+p-2})$ plays a key role.
Then the tightness of $(\mu_n)_{n \in \mathbb N_+}$ follows and we get the existence of an invariant measure for $(P_t)_{t \ge 0}$.
In combination with the uniqueness result, we obtain the existence of a unique and thus ergodic invariant measure for $(P_t)_{t \ge 0}$ in the second main result, Theorem \ref{main2}.

\iffalse
Finally, we establish a uniform exponential moments' estimation of $\sup_{t \ge 0} \mathbb E \exp(\epsilon \|\cdot\|_p^p)$ with some positive constant
$\epsilon$ for the solution (see \eqref{exp} in Lemma \ref{lm-exp}).
This shows the Gaussian concentration property for the distribution of the solution, which will be used to derive the contractivity property of $(P_t)_{t \ge 0}$.
It is worthy pointed out that such exponential moments' estimation in existing literature has the form $\sup_{t \ge 0} \mathbb E \exp(\epsilon \|\cdot\|_H^2)$ on the underlying Hilbert space $(H, \|\cdot\|_H)$, while our estimate includes a higher index $p$ when $p>2$.
This enables us to obtain the supercontractivity of $(P_t)_{t \ge 0}$ in the last main result, Theorem \ref{main3}.
\fi

The rest of the paper is organized as follows.
Some preliminaries, assumptions, and main results are given in the next section.
We derive a uniform pathwise estimate to get the existence of a unique global solution to Eq. \eqref{rd} in Section \ref{sec3}.
In another part of Section \ref{sec3}, we construct the coupling and derive a uniform pathwise estimation for this coupling.
These estimations ensure the well-posedness of the coupling and will be used in the last section to derive Harnack inequalities and to prove the main results, Theorems \ref{main1} and \ref{main2}.

\section{Preliminaries and Main Results}
\label{sec2}

Let $\mathscr O \subset \mathbb R^d$, $d \ge 1$, be an open, bounded Lipschitz domain.
Throughout $p \ge 2$ is a fixed constant.
Denote by $L^p=L^p(\mathscr O)$ the usual Lebesgue space on $\mathscr O$ with norm $\|\cdot\|_p$, $\mathcal B_b(L^p)$ the class of bounded measurable functions on $L^p$, and $\mathcal B^+_b(L^p)$ the set of positive functions in $\mathcal B_b(L^p)$.
For a function $\phi \in \mathcal B_b(L^p)$, define
\begin{align*}
\|\phi\|_\infty=\sup_{x \in L^p} |\phi(x)|,  \quad
\|\nabla \phi\|_\infty=\sup_{x \in L^p} |\nabla \phi|(x),
\end{align*}
where $|\nabla \phi|(x)=\limsup_{y \rightarrow x} |\phi(y)-\phi(x)|/\|y-x\|$, $x \in L^p$.
In particular, when $p=2$, $L^2$ is a Hilbert space with the norm $\|\cdot\|:=\|\cdot\|_2$ and the inner product $(\cdot, \cdot)$.
Let $(W_t)_{t \ge 0}$ be a $U:=L^2$-valued cylindrical Wiener process with respect to a complete filtered probability space $(\Omega, \mathscr F, (\mathscr F_t)_{t \ge 0},\mathbb P)$ satisfying the usual condition,
i.e., there exists an orthonormal basis $(e_k)_{k=1}^\infty$ of $L^2$ and a family of independent standard real-valued Brownian motions $(\beta_k)_{k=1}^\infty $ such that
\begin{align} \label{df-W}
W_t=\sum_{k=1}^\infty e_k\beta_k(t),\quad t \ge 0.
\end{align}

Define by $F$ the Nemytskii operator associated with $f$, i.e.,
\begin{align} \label{df-F}
F(x)(\xi)=f(x(\xi)), \quad x \in L^p, \ \xi \in \mathscr O.
\end{align}
The precise assumptions on $f$ and $F$ are given in Assumption \ref{ap-f} and Remark \ref{rk-F}.
Then Eq. \eqref{rd} can be rewritten as the stochastic evolution equation
\begin{align}\label{see}
{\rm d}X_t &=(A X_t+F(X_t)) {\rm d}t+ G{\rm d}W_t,
\end{align}
with the initial datum $X_0=x \in L^p$, where $A$ is the Dirichlet Laplacian operator on $L^p$, $F$ is the Nemytskii operator defined in \eqref{df-F} associated with $f$, and $W$ is a $U$-valued cylindrical Wiener process given in \eqref{df-W}.

It is well-known that the Dirichlet Laplacian operator $A$ in Eq. \eqref{see} generates an analytic $\mathcal C_0$-semigroup in $L^p$, denoted by $(S^p_t)_{t \ge 0}$, for each $p \ge 2$; see, e.g., \cite[Theorem 1.4.1]{Dav89}.
These semigroups are consistent, in the sense that $S^{p_1}_t x=S^{p_2}_t x$, for all $t \ge 0$, $x \in L^{p_1} \cap L^{p_2}$, and $p_1, p_2 \ge 2$.
Then we shall denote all $(S^p_t)_{t \ge 0}$, $p \ge 2$, by $(S_t)_{t \ge 0}$, if there is no confusion.
In Section \ref{sec4}, we also need the Sobolev--Slobodeckij spaces $W^{\beta,p}$ and $W_0^{\beta,p}:=\{\phi \in W^{\beta,p}:\ \phi|_{\partial \mathscr O}=0\}$, with $\beta\in (0,1)$, whose norm is defined by
\begin{align} \label{sob-slo}
\|\phi\|_{\beta,p}
:=\bigg(\|\phi\|_p^p
+\int_{\mathscr O}\int_{\mathscr O}
\frac{|\phi(\xi)-\phi(\eta)|^p}{|\xi-\eta|^{d+\beta p}}
{\rm d}\xi {\rm d}\eta \bigg)^\frac1p.
\end{align}
It is known that the following compact embedding holds true (see, e.g., \cite[Theorem 2.3.4]{Agr15}):
\begin{align} \label{emb}
W_0^{\beta, p} \subset L^p,
\quad \beta \in (0, d/p).
\end{align}
Moreover, $(S_t)_{t \ge 0}$ satisfies the following ultracontractivity (see, e.g., \cite[Section 2.1]{Cer03(PTRF)}):
\begin{align} \label{ult}
\|S_t u\|_{\beta, \mathbf r} \le C e^{-\lambda_1 t} t^{-(\frac\beta2+\frac{d(\mathbf r-\mathbf s)}{2\mathbf r \mathbf s})} \|u\|_{\mathbf s},
\quad t>0, \ u \in L^\mathbf s,
\end{align}
for all $\beta \in (0,1)$ and $1 \le \mathbf s \le \mathbf r \le \infty$, where $\lambda_1>0$ is the first eigenvalue of $A$.
For convenience, here and what follows, we frequently use the generic constant $C$, which may be different in each appearance.
When $p=2$, the following Poincar\'e inequality holds:
\begin{align} \label{poin}
\|\nabla u\|_2 \ge \lambda_1 \|u\|_2, \quad u \in H_0^1,
\end{align}
where $H_0^1=H_0^1(\mathscr O)$ denotes the space of weakly differentiable functions, vanishing on the boundary $\partial \mathscr O$, whose derivatives belong to $L^2$.

\subsection{Main assumptions and results}

Let us give the following assumptions on the data of Eq. \eqref{rd}.
We begin with the conditions on the drift function $f$.

\begin{ap} \label{ap-f}
There exist constants $L_f \in \mathbb R$, $\theta, L'_f>0$, and $q \ge 2$ such that for all $\xi, \eta \in \mathbb R$,
\begin{align}
& (f(\xi)-f(\eta)) (\xi-\eta) \le L_f |\xi-\eta|^2 - \theta |\xi-\eta|^q, \label{f} \\
& |f(\xi)-f(\eta)| \le L'_f (1+|\xi|^{q-2}+|\eta|^{q-2})|\xi-\eta|.  \label{f+}
\end{align}
\end{ap}

\begin{rk} \label{rk-f}
A motivating example of $f$ such that Assumption \ref{ap-f} holds true is a polynomial of odd order $q-1$ with a negative leading coefficient (for the stochastic Allen--Cahn equation, $q= 4$), perturbed with a Lipschitz continuous function; see, e.g., \cite[Example 7.8]{PZ14}.
\end{rk}

\begin{rk} \label{rk-F}
It follows from \eqref{f}-\eqref{f+} that the Nemytskii operator $F$, defined in \eqref{df-F}, associated with $f$  is a well-defined, continuous operator from $L^q$ to $L^{q'}$, with $q'=q/(q-1)$.
Moreover,
\iffalse
\begin{align} \label{F}
& \<F(u)-F(v), u-v\>
 \le L_f \|u-v\|^2-\theta \|u-v\|_q^q,
\end{align}
\fi
\begin{align*}
& \<F(u)-F(v), u-v\>
 \le L_f \|u-v\|^2-\theta \|u-v\|_q^q,
 \quad \forall~u, v\in L^q,
\end{align*}
where $\<\cdot, \cdot\>$ is the dualization between $L^{q'}$ and $L^q$ with respect to $L^2$.
\end{rk}

To perform the assumption of the noise part, as we consider the Banach spaces $(L^p)_{p \ge 2}$, let us first recall the required materials of stochastic calculus in Banach spaces, especially the M-type 2 spaces and the $\gamma$-radonifying operators.
It is known that the stochastic calculations in Banach space depend heavily on the geometric structure of the underlying spaces.

We first recall the definitions of the M-type for a Banach space.
Let $E$ be a Banach space and $(\epsilon_n)_{n\in \mathbb N_+}$ be a Rademacher sequence in a probability space $(\Omega',\mathscr F',\mathbb P')$, i.e., a sequence of independent random variables taking the values $\pm 1$ with the equal probability $1/2$.
Then $E$ is called of M-type $2$ if there exists a constant $\tau^M \ge 1$ such that
\begin{align*}
\|f_N\|_{L^2( \Omega'; E)}
\le \tau^M \Big(\|f_0\|_{L^2(\Omega';E)}^2
+\sum_{n=1}^N \|f_n-f_{n-1}\|_{L^2(\Omega';E)}^2 \Big)^\frac12
\end{align*}
for all $E$-valued square integrable martingales $\{f_n\}_{n=0}^N$.

For stochastic calculus in Banach spaces, the so-called $\gamma$-radonifying operators play the important roles instead of Hilbert--Schmidt operators in Hilbert settings.
Let $(\gamma_n)_{n \in \mathbb N_+}$ be a sequence of independent
$\mathcal N(0,1)$-random variables in a probability space $(\Omega',\mathscr F',\mathbb P')$.
Denote by $\mathcal L(U, E)$ the space of linear operators from $U$ to $E$.
An operator $R\in \mathcal L(U, E)$ is called $\gamma$-radonifying if there exists an orthonormal basis $(h_n)_{n\in \mathbb N_+}$ of $U$ such that the Gaussian series $\sum_{n\in \mathbb N_+}\gamma_nRh_n$ converges in $L^2(\Omega';E)$.
In this situation, it is known that the number
\begin{align*}
\|R\|_{\gamma(U, E)}:=\Big\|\sum_{n\in \mathbb N_+}\gamma_nRh_n\Big\|_{L^2(\Omega';E)}
\end{align*}
does not depend on the sequence $(\gamma_n)_{n \in \mathbb N_+}$ and the basis $(h_n)_{n\in \mathbb N_+}$, and
it defines a norm on the space $\gamma(U, E)$ of all $\gamma$-radonifying operators from $U$ to $E$.
In particular, if $E$ reduces to a Hilbert space, then $\gamma(U, E)$ coincides with the space of all Hilbert--Schmidt operators from $U$ to $E$.

Let $T>0$ and $(E, \|\cdot\|_E)$ be an M-type 2 space.
For any $\gamma(U, E)$-valued adapted process $\Phi \in L^\mathbf s(\Omega; L^2(0,T;\gamma(U, E)))$ with $\mathbf s \ge 2$, the following one-sided Burkholder inequality for the $E$-valued stochastic integral $\int_{0}^t \Phi_r{\rm d}W_r$ holds for some constant $C=C(\mathbf s)$ (see, e.g., \cite[Theorem 2.4]{Brz97(SSR)} or \cite[(3)]{HHL19(JDE)}):
\begin{align}\label{bdg}
\mathbb E \sup_{t\in [0,T]} \Big\|\int_{0}^t \Phi_r{\rm d}W_r\Big\|_E^\mathbf s
\le C \mathbb E \Big(\int_0^T \|\Phi\|^2_{\gamma(U, E)} {\rm d} t\Big)^\frac \mathbf s 2.
\end{align}
Coming back to our case, it is known that $(L^p)_{p \ge 2}$ are M-type 2 spaces.
For more details about definitions and properties of M-type 2 spaces and
$\gamma$-radonifying operators, we refer to \cite{Brz97(SSR)}.

With these preliminaries, now we can give our assumptions on $G$ appearing in Eq. \eqref{see}.
Let $\mathcal L_S(U)$ be the set of all densely defined closed linear operators $(L, {\rm Dom}(L))$ on $U$ such that for every $s>0$, $S_s L$ extends to a unique operator in $\gamma(U, L^p)$, which is again denoted by $S_s L$.
Assume that $G \in L_S(U)$ and denote by $W_A$ the stochastic convolution, also known as the Ornstein--Uhlenbeck process, associated with Eq. \eqref{see}, i.e.,
\begin{align} \label{df-wa}
W_A(t)=\int_0^t S_{t-r} G {\rm d}W_r, \quad t \ge 0.
\end{align}
It is clear that $W_A$ is the mild solution of the linear equation
\begin{align} \label{eq-wa}
{\rm d}Z_t &=A Z_t {\rm d}t+G {\rm d}W_t,
\quad Z_0=0,
\end{align}
i.e., Eq. \eqref{see} with $F=0$, with vanishing initial datum.
We always assume that, for each fixed $T>0$,
$\int_0^T \|S_t G\|^2_{\gamma(U, L^p)} {\rm d}t<\infty$,
so it follows from the Burkholder inequality \eqref{bdg} that
\begin{align*}
\mathbb E \|W_A(T)\|_p^2
\le C \int_0^T \|S_t G\|^2_{\gamma(U, L^p)} {\rm d}t <\infty,
\end{align*}
which shows that $W_A(T)$ possesses the bounded second moment in $L^p$.
Moreover, we perform the following stronger assumption, as we shall handle the polynomial drift function $f$ which satisfies Assumption \ref{ap-f}.

\begin{ap} \label{ap-wa}
$W_A$ has a continuous version in $L^{q+p-2}$ such that
\begin{align} \label{wa1}
\sup_{t \ge 0} \mathbb E \|W_A(t)\|_{q+p-2}^{q+p-2} <\infty.
\end{align}
\end{ap}

The above Assumption \ref{ap-wa} on the Ornstein--Uhlenbeck process $W_A$ \eqref{df-wa} is necessary for the study of the well-posedness, Harnack inequalities, and ergodicity in Sections \ref{sec3} and \ref{sec4}.
To indicate that Assumption \ref{ap-wa} is natural, we remark that the condition \eqref{wa1} is valid in various of applications, even when $G$ is an unbounded operator.

\begin{ex} \label{ex-wa}
Consider 1D Eq. \eqref{see} with $\mathscr O=(0,1)$, $U=L^2(0,1)$ (with a uniformly bounded orthonormal basis $(e_k=\sqrt2 \sin (k\pi x))_{k \in \mathbb N_+}$), and $G: =(-\Delta)^{\theta/2}$ for some $\theta<1/2$.
In particular, $\theta=0$ corresponds to white noise and $G$ is an unbounded operator when $\theta \in (0, 1/2)$.

It is clear that $G$ is a densely defined closed linear operator on $U$.
To show $G \in L_S(U)$, we note that for $\mathbf r \ge 2$, $L^\mathbf r$ is a Banach function space with finite cotype, so an operator $\Phi\in \gamma(U, L^\mathbf r)$ if and only if $(\sum_{k=1}^\infty (\Phi e_k)^2)^{1/2}$ belongs to $L^\mathbf r$ and there exist a constant  $C>0$ such that (see \cite[Lemma 2.1]{NVW08(JFA)})
\begin{align} \label{gamma}
\frac1C \Big\|\sum_{k=1}^\infty (\Phi e_k)^2\Big\|_{\mathbf r/2}
\le \|\Phi\|^2_{\gamma(U, L^\mathbf r)}
&\le C \Big\|\sum_{k=1}^\infty (\Phi e_k)^2\Big\|_{\mathbf r/2},
\quad \Phi\in \gamma(U, L^\mathbf r).
\end{align}
It follows from the above estimates with $\Phi=S_r G$ for $r>0$ that
\begin{align*}
\|S_r G\|^2_{\gamma(U, L^\mathbf r)}
\le C \Big\|\sum_{k=1}^\infty (S_r G e_k)^2\Big\|_{\mathbf r/2}
\le C \sum_{k=1}^\infty e^{-2 \lambda_k r} \lambda_k^\theta,
\end{align*}
which is convergent if and only if $\theta<1/2$,
where $\lambda_k=\pi^2 k^2$, $k \in \mathbb N_+$.
This shows $S_r G \in \gamma(U, L^p)$ for every $r>0$ and thus $G \in L_S(U)$.

Finally, one can use the Burkholder inequality \eqref{bdg} and the above equivalence relation \eqref{gamma} to show that
$W_A$ belongs to $\mathcal C([0, \infty); L^{\mathbf s}(\Omega; L^\mathbf r))$ for any $\mathbf s, \mathbf r \ge 2$, following an argument used in \cite[(2.17) in Lemma 2.2]{LQ20(IMA)}.
Indeed, for any $t \ge 0$,
\begin{align*}
\sup_{t \ge 0} \mathbb E \|W_A(t)\|_{\mathbf r}^{\mathbf s}
& \le C(\mathbf s) \sup_{t \ge 0} \Big(\int_0^t \|S_r (-\Delta)^\frac\theta2\|^2_{\gamma(L^2, L^\mathbf r)} {\rm d}r\Big)^\frac \mathbf s 2 \\
& \le C(\mathbf s) \Big(\int_0^\infty \Big\|\sum_{k=1}^\infty e^{-2 \lambda_k r} \lambda_k^\theta e_k^2\Big\|_{\mathbf r/2} {\rm d}t\Big)^\frac \mathbf s 2 \\
& \le  C(\mathbf s) \Big(\sum_{k=1}^\infty \lambda_k^\theta \|e_k\|^2_\mathbf r
\Big(\int_0^\infty e^{-2 \lambda_k r} {\rm d}t \Big) \Big)^\frac \mathbf s 2 \\
& \le  C(\mathbf s) \Big(\sum_{k=1}^\infty \frac1{2 \lambda_k^{1-\theta}} \Big)^\frac \mathbf s 2  <\infty.
\end{align*}
This shows \eqref{wa1} with $\mathbf s=\mathbf r=q+p-2$.
\end{ex}

In the derivation of the existence of an invariant measure in Theorem \ref{tm-erg} when $q>2$, we need the following Sobolev regularity of the Ornstein--Uhlenbeck process $(W_A(t))_{t \ge 0}$.

\begin{ap} \label{ap-wa+}
There exists a $\beta_0>0$ such that
\begin{align} \label{wa+}
\sup_{t \ge 0} \mathbb E \|W_A(t)\|_{\beta_0, p}<\infty.
\end{align}
\end{ap}

\begin{ex} \label{ex-wa+}
As in Example \ref{ex-wa},
\begin{align*}
\sup_{t \ge 0}\mathbb E \|W_A(t)\|_{\beta_0, \mathbf r}^{\mathbf s}
&\le C(\mathbf s) \Big(\int_0^\infty \|(-A)^{\frac{\beta_0+\theta}2} S_r G\|^2_{\gamma(L^2, L^\mathbf r)} {\rm d}r\Big)^\frac \mathbf s 2 \\
& \le  C(\mathbf s) \Big(\sum_{k=1}^\infty \frac1{\lambda^{1-(\beta_0+\theta)}_k} \Big)^\frac \mathbf s 2,
\end{align*}
which is convergent if and only if $\beta_0+\theta<1/2$.
As $\theta<1/2$ in Example \ref{ex-wa}, this shows \eqref{wa+} with $\mathbf r=p$ and $\mathbf s=1$ for any $\beta_0 \in (0, 1/2-\theta)$.
\end{ex}

To derive Wang-type Harnack inequalities and the ergodicity for $(P_t)_{t \ge 0}$ in Section \ref{sec4}, we need the following standard elliptic condition.

\begin{ap}\label{ap-ell}
$G G^*$ is invertible, with inverse $(G G^*)^{-1}$, such that $G^{-1}:=G^* (G G^*)^{-1}$ is a bounded linear operator on $U$:
\begin{align} \label{G-1}
\|G^{-1}\|_\infty:=\sup_{x \in U:~ x \neq 0} \frac{\| G^{-1} x\|}{\|x\|}<\infty.
\end{align}
\end{ap}

\iffalse
\begin{rk}
In general, there exist examples such that $G G^*$ is invertible but $G$ is not. For instance, let $p=2$, $U=L^2$ with orthonormal basis $(e_k)_{k \in \mathbb N_+}$, and $G e_1=0$ and $G e_k=e_{k-1}$, $k \ge 2$, so that $G$ is not invertible.
But $G^* e_k=e_{k+1}$, $k \ge 1$, and thus $G G^*={\rm Id}$, which is invertible.
\end{rk}
\fi

\begin{rk}  \label{rk-G-1}
Let $\theta \in [0, 1/2)$.
Then the operator $G: =(-\Delta)^{\theta/2}$ defined in Example \ref{ex-wa} is invertible with bounded inverse $G^{-1}: =(-\Delta)^{-\theta/2}$.
\end{rk}

\subsection{Main results}

Now we are in the position to present our main results.
Let us first recall some definitions.
The Markov semigroup $(P_t)_{t \ge 0}$ defined in \eqref{pt} is called of strong Feller if $P_T \phi \in \mathcal C_b(L^p)$ for any $T>0$ and $\phi \in \mathcal B_b(L^p)$.
A probability measure $\mu$ on $L^p$ is said to be an invariant measure of
$(P_t)_{t \ge 0}$ or of Eq. \eqref{see}, if
\begin{align*}
\int_{L^p} P_t \phi(x) \mu({\rm d}x)=\int_{L^p} \phi(x) \mu({\rm d}x),
\quad \phi \in \mathcal B_b(L^p), \ t \ge 0.
\end{align*}
{\color{blue} df of ergodicity}

\iffalse
If $(P_t)_{t \ge 0}$ has an invariant measure $\mu$, then $(P_t)_{t \ge 0}$ is called hypercontractive if $\|P_t\|_{2 \rightarrow 4}=1$ for some $t>0$,
where $\|\cdot\|_{\mathbf r \rightarrow \mathbf s}$ denotes the operator norm from $L^\mathbf r(\mu)$ to $L^\mathbf s(\mu)$: $\|P_t\|_{\mathbf r \rightarrow \mathbf s}=\sup\{\|P_t \phi\|_{L^\mathbf s(\mu)}: \ \mu(\phi^\mathbf r) \le 1 \}$, and $(P_t)_{t \ge 0}$ is called supercontractive if $\|P_t\|_{2 \rightarrow 4}<\infty$ for all $t>0$.
\fi

Our first main result is the following log-Harnack inequality and power-Harnack inequality, from which the uniqueness of the invariant measure, if it exists, for $(P_t)_{t \ge 0}$ follows.

\begin{tm} \label{main1}
Let Assumptions \ref{ap-f}, \ref{ap-wa}, and \ref{ap-ell} hold.
For any $T > 0$, $\mathbf s >1$, $x, y \in L^p$, and $\phi \in \mathcal B^+_b(L^p)$,
\begin{align}
& P_T \log \phi(y) \le \log P_T \phi(x)
+\frac{\lambda \|G^{-1}\|_\infty^2}{e^{2 \lambda T}-1} \|x-y\|_p^2, \label{har-log} \\
& (P_T \phi(y))^\mathbf s \le P_T \phi^\mathbf s(x)
\exp\Big(\frac{\mathbf s \lambda \|G^{-1}\|_\infty^2 \|x-y\|_p^2}{(\mathbf s-1) (e^{2 \lambda T}-1)}\Big).  \label{har-pow}
\end{align}
Consequently, $(P_t)_{t \ge 0}$ has at most one invariant measure.
\end{tm}

The next main result is the existence of an invariant measure for $(P_t)_{t \ge 0}$.
In combination with the uniqueness result in Theorem \ref{main1}, $(P_t)_{t \ge 0}$ possesses exactly one ergodic invariant measure.

\begin{tm} \label{main2}
Let Assumptions \ref{ap-f} and \ref{ap-wa} hold.
Assume that $q>2$ such that
\begin{align} \label{d}
d<\frac{2p(q+p-2)}{(p-1)(q-2)},
\end{align}
and Assumption \ref{ap-wa+} holds, or $\lambda$ defined by
\begin{align} \label{lam}
\lambda:= -L_f + \theta \chi_{q=2, p \neq 2}
+ \lambda_1 \chi_{q \neq 2, p=2}
+(\lambda_1+\theta) \chi_{q=p=2},
\end{align}
is positive.
Then $(P_t)_{t \ge 0}$ has an invariant measure $\mu$ with full support on $L^p$ such that $\mu( \|\cdot\|_{q+p-2}^{q+p-2})<\infty$.
Assume furthermore that Assumption \ref{ap-ell} holds, then $\mu$ is the unique ergodic invariant measure of $(P_t)_{t \ge 0}$.
\end{tm}

\iffalse
The last main result is the contractivity of $(P_t)_{t \ge 0}$.

\begin{tm} \label{main3}
Let Assumptions \ref{ap-f} and \ref{ap-ell}, and the condition \eqref{wa2} hold for a sufficiently large $\epsilon_0>0$.
\begin{enumerate}
\item
Assume that $\lambda$ defined in \eqref{lam} is positive, then $(P_t)_{t \ge 0}$ is hypercontractive.
Consequently, $(P_t)_{t \ge 0}$ is compact on $L^2(\mu)$ for large $t>0$, and there exist constants $C, \lambda_0 >0$ such that for all $t \ge 0$,
\begin{align}
& \mu(P_t \phi \log P_t \phi) \le C e^{-\lambda_0 t} \mu(\phi \log \phi),
\quad \ f\in \mathcal B_b^+(L^p), \ \mu(\phi)=1,  \label{ent} \\
& \|P_t \phi-\mu(\phi)\|_{L^2(\mu)} \le C e^{-\lambda_0 t}   \|\phi-\mu(\phi)\|_{L^2(\mu)},
\quad \phi \in L^2(\mu). \label{spe-gap}
\end{align}

\item
Let $q>2$ and Assumptions \ref{ap-f} and \ref{ap-ell} hold.
Assume that \eqref{wa2} holds for all $\epsilon_0>0$ if $p=2$ or holds for a sufficiently large $\epsilon_0>0$ if $p>2$ , then $(P_t)_{t \ge 0}$ is supercontractive.
\end{enumerate}
\end{tm}
\fi

\section{Coupling and Moments' Estimations}
\label{sec3}

The main aims of this section are to show the existence of a unique global solution to Eq. \eqref{rd} and to construct a well-defined coupling process for this solution process.
We also derive several uniform a priori estimates on moments of these two processes, which will be used in Section \ref{sec4} to derive Wang-type Harnack inequalities and the ergodicity of $(P_t)_{t \ge 0}$.

\subsection{Well-posedness and moments' estimations}

Let us first recall that an $L^p$-valued process $(X_t)_{t \in [0, T]}$ is called a mild solution of Eq. \eqref{see} with the initial datum $X_0=x$ if $\mathbb P$-a.s.
\begin{align}\label{mild}
X_t &=S_t x+ \int_0^t S_{t-r} F(X_r) {\rm d}r+W_A(t),
\quad t \in [0, T].
\end{align}
From Remark \ref{rk-F}, the deterministic convolution in Eq. \eqref{mild} makes sense.
Define $Z=X-W_A$.
It is clear that $X$ is a mild solution Eq. \eqref{see} if and only if $Z$ is a mild solution of the random PDE
\begin{align} \label{Z}
\partial_t Z_t =\Delta Z_t+ F(Z_t+W_A(t)),
\quad Z_0=x.
\end{align}

The following results show the existence of a unique mild solution of Eq. \eqref{see} which is a Markov process and depends on the initial data continuously in pathwise sense.

\begin{lm} \label{lm-lq}
Let $T>0$, $x \in L^p$, and Assumptions \ref{ap-f} and \ref{ap-wa} hold.
Eq. \eqref{rd} with the initial datum $X_0=x$ possesses a unique mild solution, in $\mathcal C([0,T]; L^p) \cap L^{q+p-2}(0, T; L^{q+p-2})$ $\mathbb P$-a.s., which is a Markov process.
Moreover, there exists a constant $\lambda$ defined in \eqref{lam} such that
\begin{align} \label{con}
& \|X_t^x-X_t^y\|_p \le e^{- \lambda t}  \|x-y\|_p,
\quad t \ge 0, \ x, y \in L^p.
\end{align}
\end{lm}

\begin{proof}
From \eqref{f+}, $f$ is locally Lipschitz continuous, so it is clear that both Eq. \eqref{Z} with $Z=X-W_A$ and Eq. \eqref{see} exist local solutions on $[0, T_0)$ for some $T_0 \in (0, T]$.
To extend this local solution to the whole time interval $[0, T]$, it suffices to give a priori uniform estimations for $Z$ and $X$.

As Eq. \eqref{Z} is a pathwise random PDE, we test $p |Z_t|^{p-2} Z_t$ on this equation with $t \in [0, T_0)$ and use the conditions \eqref{f}-\eqref{f+} and Young inequality to obtain
\begin{align*}
& \partial_t \|Z_t\|_p^p
+ p(p-1) \int_{\mathscr O} |Z_t|^{p-2} |\nabla Z_t|^2 {\rm d}\xi  \\
& =  p\< |Z_t|^{p-2}Z_t, F(Z_t+W_A(t))-F(W_A(t))\>
+p \<|Z_t|^{p-2}Z_t, F(W_A(t)) \> \\
& \le  p L_f \|Z_t\|_p^p
- p \theta \|Z_t\|_{q+p-2}^{q+p-2}
+p \<|Z_t|^{p-2}Z_t, F(W_A(t)) \> \nonumber \\
& \le  p L_f \|Z_t\|_p^p
- p \theta_1 \|Z_t\|_{q+p-2}^{q+p-2}
+C \|F(W_A(t))\|_\frac{q+p-2}{q-1}^\frac{q+p-2}{q-1}\nonumber \\
& \le  C (1+\|W_A(t)\|_{q+p-2}^{q+p-2})
+p L_f \|Z_t\|_p^p  - p \theta_1 \|Z_t\|_{q+p-2}^{q+p-2},
\end{align*}
where $\theta_1$ could be chosen as any positive number which is smaller than $\theta$.
This yields
\begin{align} \label{Zp}
& \|Z_t\|_p^p
+ p \theta_1 \int_0^t \|Z_r\|_{q+p-2}^{q+p-2} {\rm d} r
+ p(p-1) \int_0^t \int_{\mathscr O} |Z_t|^{p-2} |\nabla Z_t|^2 {\rm d}\xi  {\rm d} r  \\
& \le  \|x\|_p^p+C \int_0^t (1+\|W_A(r)\|_{q+p-2}^{q+p-2}) {\rm d} r
+p L_f \int_0^t \|Z_r\|_p^p  {\rm d} r. \nonumber
\end{align}
Using Gronwall inequality, we obtain
\begin{align*}
& \|Z_t\|_p^p
+ p \theta_1 \int_0^t \|Z_r\|_{q+p-2}^{q+p-2} {\rm d} r
+ p(p-1) \int_0^t \int_{\mathscr O} |Z_t|^{p-2} |\nabla Z_t|^2 {\rm d}\xi  {\rm d} r  \\
& \le e^{p L_f t} \Big(\|x\|_p^p+C \int_0^t (1+\|W_A(r)\|_{q+p-2}^{q+p-2}) {\rm d} r \Big).
\end{align*}

The above uniform estimate, in combination with the condition \eqref{wa1} of $W_A$ in Assumption \ref{ap-wa}, implies the global existence of a mild solution $Z$ to Eq. \eqref{Z} on $[0, T]$ in $\mathcal C([0,T]; L^p) \cap L^{q+p-2}(0, T; L^{q+p-2})$ $\mathbb P$-a.s.
Taking into account the relation $X=Z+W_A$ and the condition \eqref{wa1}, we obtain a global mild solution $X$ to Eq. \eqref{see}.

To show the continuous dependence \eqref{con}, let us note that
\begin{align*}
\partial_t (X_t^x-X_t^y) &=A (X_t^x-Y_t^y)+F(X_t^x)-F(Y_t^y).
\end{align*}
Testing $p |X_t^x-X_t^y|^{p-2} (X_t^x-X_t^y)$ on the above equation, using integration by parts formula, and applying the condition \eqref{f}, we obtain
\iffalse
\begin{align*}
& \partial_t \|X_t^x-X_t^y\|_p^p
+ p(p-1) \int_{\mathscr O} |X_t^x-X_t^y|^{p-2} |\nabla (X_t^x-X_t^y)|^2 {\rm d} \xi \\
& \le p L_f \|X_t^x-X_t^y\|_p^p - p \theta \|X_t^x-X_t^y\|_{q+p-2}^{q+p-2},
\end{align*}
\fi
\begin{align} \label{con-xy}
& \|X_t^x-X_t^y\|_p^p
+p \theta \int_0^t \|X_r^x-X_r^y\|_{q+p-2}^{q+p-2} {\rm d}r \nonumber \\
&\quad + p(p-1) \int_0^t \int_{\mathscr O} |X_r^x-X_r^y|^{p-2} |\nabla (X_r^x-X_r^y)|^2 {\rm d} \xi {\rm d}r  \nonumber \\
& \le \|x-y\|_p^p + p L_f \int_0^t \|X_r^x-X_r^y\|_p^p {\rm d}r.
\end{align}
We conclude \eqref{con} with $\lambda= -L_f$ by Gronwall inequality.
When $q=2$, then \eqref{con} holds with $\lambda=\theta-L_f$, as one can subtract the first integral on the left-hand side of the above inequality; while $p=2$, then using the Poincar\'e inequality \eqref{poin} yields \eqref{con} with $\lambda=\lambda_1-L_f$.
Similarly, \eqref{con} holds with $\lambda=\lambda_1+\theta-L_f$ when $q=p=2$.
These statements show \eqref{con} with $\lambda$ given by \eqref{lam}.

The pathwise continuous dependence clearly implies the uniqueness of the solution to Eq. \eqref{see}.
One can also show the Markov property for this solution using a standard method, see, e.g., \cite[Theorem 9.21]{PZ14}.
This completes the proof.
\end{proof}

\begin{rk}
The pathwise estimate \eqref{con} immediately yields the following estimate between any two solutions in $\mathbf r$-Wasserstein distance for any $\mathbf r \ge 1$:
\begin{align*}
{\rm W}_{\mathbf r}(\mu^1_t, \mu^2_t)
: =\inf (\mathbb E \|X_t^{x_1}-X_t^{x_2}\|^{\mathbf r})^\frac1\mathbf r
\le e^{- \lambda t} {\rm W}_{\mathbf r}(\mu^1_0, \mu^2_0),
\quad t \ge 0,
\end{align*}
where $(\mu^i_0)_{i=1}^2$ are two measures on $L^p$, $(X^{x_i}_t)_{i=1}^2$ are the solutions to Eq. \eqref{see} starting from $(x_i)_{i=1}^2$ of laws $(\mu^i_0)_{i=1}^2$, and the infimum runs over all random variables $(X^{x_i}_t)_{i=1}^2$ with laws $(\mu^i_t)_{i=1}^2$, $t \ge 0$, respectively.
Similar contraction-type estimate in $2$-Wasserstein distance on $\mathbb R^d$ had been investigated in \cite{BGG12(JFA)}.
\end{rk}

\iffalse
Under the stronger condition \eqref{wa2} instead of \eqref{wa1}, we have the following exponential moments' estimation, also known as exponential integrability, which is related to Fernique thoerem and logarithmic Sobolev inequality \cite{AMS94(JFA), BG99(JFA)}.
It will be used to show the Gaussian concentration property for the invariant measure of $(P_t)_{t \ge 0}$, if it exists, in Section \ref{sec4}.

{\color{red}
\begin{lm} \label{lm-exp}
Let $x \in L^p$ and Assumptions \ref{ap-f} and \ref{ap-wa} hold for some $\epsilon_0>0$.
Assume that $q>2$ or $\lambda$ defined in \eqref{lam} is positive, then there exists a constant $\epsilon_1>0$ such that
\begin{align} \label{exp}
\sup_{t \ge 0} \mathbb E \exp(\epsilon_1 \|X_t^x\|_p^2)<\infty.
\end{align}
If \eqref{wa2} is valid for all $\epsilon_0>0$, then \eqref{exp} holds for all $\epsilon_1>0$.
\end{lm}

\begin{proof}
By \eqref{Zp} and the chain rule, we have
\begin{align*}
& \partial_t \exp( \epsilon \|Z_t\|_p^p )
=\epsilon \exp( \epsilon \|Z_t\|_p^p) \partial_t \|Z_t\|_p^p \\
& \le \epsilon \exp( \epsilon \|Z_t\|_p^p)
\Big( - p(p-1) \int_{\mathscr O} |Z_t|^{p-2} |\nabla Z_t|^2 {\rm d}\xi \\
& \quad + C (1+\|W_A(t)\|_{q+p-2}^{q+p-2})
+p L_f \|Z_t\|_p^p - p\theta_1 \|Z_t\|_{q+p-2}^{q+p-2} \Big),
\end{align*}
where $\theta_1$ could be chosen as any positive number smaller than $\theta$.
When $q >2$, by the continuous embedding $L^{q+p-2} \subset L^p$ and Young inequality, we have
\begin{align*}
- p(p-1) \int_{\mathscr O} |Z_t|^{p-2} |\nabla Z_t|^2 {\rm d}\xi
+ p L_f \|Z_t\|_p^p - p \theta_1 \|Z_t\|_{q+p-2}^{q+p-2}
\le C - \theta_2 \|Z_t\|_p^p,
\end{align*}
for some positive constant $\theta_2$.
This is also true when $q=2$ and $L_f<\theta$, with $\theta_2=\theta_1-L_f>0$, as one can choose $\theta_1$ as close as possible to $\theta$.
When $q=p=2$, the Poincar\'e inequality yields that
\begin{align*}
& - p(p-1) \int_{\mathscr O} |Z_t|^{p-2} |\nabla Z_t|^2 {\rm d}\xi
+ p L_f \|Z_t\|_p^p - p \theta_1 \|Z_t\|_{q+p-2}^{q+p-2} \\
& \le -2(\lambda_1+\theta_1-L_f) \|Z_t\|_p^p.
\end{align*}
Therefore, when $q>2$ or $\lambda >0$, there exists $\theta_3>0$ such that
\begin{align*}
\partial_t \exp( \epsilon \|Z_t\|_p^p )
\le \epsilon \exp( \epsilon \|Z_t\|_p^p)
[ C (1+\|W_A(t)\|_{q+p-2}^{q+p-2})- \theta_3 \|Z_t\|_p^p].
\end{align*}
Using the elementary inequality: for any constants $\epsilon, C_1, C_2>0$, there exists a constant $C_3:=\exp((\epsilon C_1+1)/C_2)$ such that
\begin{align*}
\epsilon e^{\epsilon x} (C_1-C_2 x) \le C_3-e^{\epsilon x},
\quad  x \ge 0,
\end{align*}
\iffalse
where we choose $C_3 \ge e^\frac{\epsilon C_1}{C_2}$ if $x \le \frac{C_1}{C_2}$, $C_3>0$ if $x>\frac{C_1+\epsilon}{C_2}$, and  $C_3  \ge e^\frac{\epsilon C_1+1}{C_2}$ if
$x \in (\frac{C_1}{C_2}, \frac{C_1+\epsilon}{C_2})$,
\fi
we obtain
\begin{align*}
\partial_t \exp( \epsilon \|Z_t\|_p^p)
\le C_{\epsilon, \theta} \exp(C_{\epsilon, \theta} \|W_A(t)\|_{q+p-2}^{q+p-2})-\exp( \epsilon \|Z_t\|_p^p),
\end{align*}
and thus
\begin{align*}
\exp( \epsilon \|Z_t\|_p^p)
\le C_{\epsilon, \theta} \exp(C_{\epsilon, \theta} \|W_A(t)\|_{q+p-2}^{q+p-2}).
\end{align*}

It follows from the relation $X=Z+W_A$ and the above inequality that
\begin{align*}
\mathbb E \exp( \epsilon \|X_t\|_p^p)
& \le \sqrt{\mathbb E \exp( \epsilon 2^p \|Z_t\|_p^p)}
\cdot \sqrt{\mathbb E \exp( \epsilon 2^p \|W_A(t)\|_p^p)} \nonumber  \\
& \le C_{\epsilon, \beta, p} \mathbb E \exp(C_{\epsilon, \beta, p} \|W_A(t)\|_{q+p-2}^{q+p-2}).
\end{align*}
Using the condition \eqref{wa2}, we get \eqref{exp}.
If \eqref{wa2} is valid for all $\epsilon_0>0$, then the above estimate yields  \eqref{exp} for any $\epsilon_1>0$.
\end{proof}
}
\fi

\iffalse
We note that the restrictions on the numbers $q$ and $p$, the dissipative constants $L_f$ and $\theta$ in \eqref{f}-\eqref{f+}, and the Poincar\'e constant $\lambda_1$ in \eqref{poin}, appearing in Lemma \ref{lm-exp}, will be frequently used in the rest of the present paper.

\begin{rk} \label{rk-pq}
When $q>2$, i.e., the dissipative drift function $f$ in \eqref{f} and \eqref{f+} has super-linear growth, the Gaussian concentration property \eqref{exp} for Eq. \eqref{see} is always valid.
Otherwise, one needs to ensure that the constant $\lambda$, given in \eqref{lam}, is positive.
This corrosponding the cases $q=2$, $p>2$, and $\theta>L_f$, or $q=2$, $p=2$, and $\lambda_1+\theta>L_f$.
In particular, in the latter case with $p=2$, one can utilize the dissipativity of the Dirichlet Laplacian operator in the Hilbert space $L^2$ to reduce the requirement of the dissipative constants in \eqref{f}: $L_f<\lambda_1+\theta$.
If $p \neq 2$, the term $\int_0^t \int_{\mathscr O} |Z_t|^{p-2} |\nabla Z_t|^2 {\rm d}\xi  {\rm d} r$, involving the Dirichlet Laplacian operator, has no dissipativity, so that the stronger restriction $L_f < \theta$ has to be imposed.
\end{rk}
\fi

\subsection{Construction of coupling and moments' estimations}

Let $T>0$ be fixed throughout the rest of Section \ref{sec3} and set
\begin{align} \label{gt}
\gamma_t:=\int_0^{T-t} e^{2 \lambda r} {\rm d}r
=\frac{e^{2 \lambda (T-t)}-1}{2 \lambda},
\quad t \in [0, T],
\end{align}
where $\lambda$ is given in \eqref{lam}.
For convention, if $\lambda=0$, we set $\frac{e^{2 \lambda t}-1}{2 \lambda}:=t$ for $t \in [0, T]$.
Then $\gamma$ is smooth, strictly positive, and strictly decreasing on $[0,T)$ (with $\gamma_T=0$) such that
\begin{align} \label{gt=}
\gamma_t'+2 \lambda \gamma_t+1=0.
\end{align}
Moreover, the integral of $\gamma^{-1}$ on $[0, T)$ diverges:
\begin{align} \label{gt>}
\int_0^T \frac1{\gamma_t} {\rm d}t=\infty.
\end{align}

Now we can define the coupling $Y$ of $X$ as the mild solution of the coupling equation
\begin{align} \label{rd-y}
{\rm d}Y_t =( AY_t+F(Y_t)+\gamma_t^{-1} (X_t-Y_t) ){\rm d}t+G {\rm d}W_t,
\end{align}
with an initial datum $Y_0=y \in L^p$.
Since the additional drift term $\gamma_t^{-1} (X_t-Y_t)$ is Lipschitz continuous for each fixed $t \in [0, T)$ and $\omega \in \Omega$, one can use similar arguments in Lemma \ref{lm-lq} to show that $Y$ is a well-defined continuous process on $[0, T)$.

\begin{rk} \label{rk-Y}
As $\gamma^{-1}$ is continuous and thus integrable on $[0, T_0] \subset [0, T)$ for any $T_0 \in (0, T)$, one can use the arguments in Lemma \ref{lm-lq} to extend the local solution to $[0, T)$.
However, it is difficult to get a uniform a priori estimation, following the idea in Lemma \ref{lm-lq}, to conclude the well-posedness of $Y$ at $T$ as $\gamma_T^{-1}$ is singular satisfying \eqref{gt>}.
\end{rk}

For each $s \in [0, T)$, we set
\begin{align} \label{v}
v_s:=\frac{G_s^{-1} (X_s-Y_s)}{\gamma_t}, \quad
\widetilde{W}_s:=W_s+\int_0^s v_r {\rm d}r,
\end{align}
and define
\begin{align} \label{M}
M_s: & =\exp\Big( -\int_0^s (v_r, {\rm d}W_r)
-\frac12 \int_0^s \|v_r\|^2 {\rm d}r\Big).
\end{align}
From \eqref{v} and \eqref{M}, $M$ can also be rewritten as
\begin{align} \label{M+}
M_s: & =\exp\Big( -\int_0^s (v_r, {\rm d}{\widetilde W}_r)
+\frac12 \int_0^s \|v_r\|^2 {\rm d}r\Big), \quad s \in [0, T).
\end{align}

It is clear that ${\mathbb Q}_s:=M_s \mathbb P$ is a probability measure.
By the representation \eqref{v} and the non-degenerate condition \eqref{G-1} in Assumption \ref{ap-ell}, we have
\begin{align} \label{est-v}
\frac12 \int_0^s \|v_r\|^2 {\rm d}r
\le \frac{ \|G^{-1}\|_\infty^2}2 \int_0^s \frac{\|X_r-Y_r\|_p^2}{\gamma_r^2} {\rm d}r.
\end{align}
We first show that for any $s \in (0, T)$, $(\widetilde W_t)_{t \in [0, s]}$ is a $U$-valued cylindrical Wiener process under the probability measure ${\mathbb Q}_s$ through Girsanov theorem ensured by the Novikov condition \eqref{nov}.

\begin{lm} \label{lm-R-}
Let Assumptions \ref{ap-f}, \ref{ap-wa}, and \ref{ap-ell} hold.
For any $s \in (0, T)$, $(\widetilde W_t)_{t \in [0, s]}$ is a $U$-valued cylindrical Wiener process under ${\mathbb Q}_s$.
\end{lm}

\begin{proof}
It follows from Eq. \eqref{see} and Eq. \eqref{rd-y} on $[0, s]$ that
\begin{align} \label{x-y-p}
\partial_t (X_t-Y_t) &=A (X_t-Y_t)+F(X_t)-F(Y_t)
- \gamma_t^{-1} (X_t-Y_t), \quad \mathbb P\text{-a.s.}
\end{align}
As in the proof of the inequality \eqref{con-xy}, we test $p |X_t-Y_t|^{p-2} (X_t-Y_t)$ on the above equation, use integration by parts formula, and apply the condition \eqref{f} to obtain
\begin{align*}
& \partial_t \|X_t-Y_t\|_p^p
+ p(p-1) \int_{\mathscr O} |X_t-Y_t|^{p-2} |\nabla (X_t-Y_t)|^2 {\rm d} \xi
\\
& \le p L_f \|X_t-Y_t\|_p^p - p \theta \|X_t-Y_t\|_{q+p-2}^{q+p-2}
- p\gamma_t^{-1} \|X_t-Y_t\|_p^p.
\end{align*}
It follows from the chain rule that
\begin{align*}
& \partial_t \|X_t-Y_t\|_p^2
=\frac2p \|X_t-Y_t\|_p^{2-p}\partial_t \|X_t-Y_t\|_p^p \\
& \le - 2(p-1) \|X_t-Y_t\|_p^{2-p} \int_{\mathscr O} |X_t-Y_t|^{p-2} |\nabla (X_t-Y_t)|^2 {\rm d} \xi \\
& +2 L_f \|X_t-Y_t\|^2_p
- 2 \theta \|X_t-Y_t\|_p^{2-p} \|X_t-Y_t\|_{q+p-2}^{q+p-2}
- 2 \gamma_t^{-1}\|X_t-Y_t\|^2_p,
\end{align*}
and thus
\begin{align*}
\partial_t \|X_t-Y_t\|_p^2
\le -2 \lambda \|X_t-Y_t\|^2_p  - 2 \gamma_t^{-1}\|X_t-Y_t\|^2_p.
\end{align*}
The product rule of differentiation and the equality \eqref{gt=} yield that
\begin{align*}
\partial_t (\gamma_t^{-1} \|X_t-Y_t\|_p^2 )
& =\gamma_t^{-1} \partial_t \|X_t-Y_t\|_p^2
-\gamma_t^{-2} \gamma'_t \|X_t-Y_t\|_p^2 \\
& \le -\gamma_t^{-2} (\gamma'_t+2 \lambda \gamma_t+2) \|X_t-Y_t\|_p^2 \\
& = - \gamma_t^{-2} \|X_t-Y_t\|_p^2.
\end{align*}
Integrating on both sides from $0$ to $s$, we obtain
\begin{align} \label{x-y}
\frac{\|X_s-Y_s\|_p^2}{\gamma_s}
+\int_0^s \frac{\|X_t-Y_t\|_p^2}{\gamma_t^2} {\rm d}t
\le \frac{\|x-y\|_p^2}{\gamma_0}, \quad \mathbb P\text{-a.s.}
\end{align}
This pathwise estimate, in combination with the estimate \eqref{est-v}, particularly implies the Novikov condition
\begin{align} \label{nov}
\mathbb E \exp \Big(\frac12 \int_0^s \|v_t\|^2 {\rm d}t \Big)
\le \exp \Big(\frac{\|G^{-1}\|_\infty^2 \|x-y\|_p^2}{2 \gamma_0}\Big)
<\infty.
\end{align}
By Girsanov theorem, $(\widetilde W_t)_{t \in [0, s]}$ is a cylindrical Wiener process under the probability measure ${\mathbb Q}_s$.
\end{proof}

Next, we will give two uniform moments' estimations \eqref{est-R} and \eqref{est-R+} for certain functionals of $(M_s)$ for all $s \in [0, T)$.
The first moments' estimation \eqref{est-R} will indicate that $(M_s)_{s \in [0, T]}$ defined in \eqref{M+} is indeed a uniformly integrable martingale.

\begin{lm} \label{lm-R}
Let Assumptions \ref{ap-f}, \ref{ap-wa}, and \ref{ap-ell} hold.
Then
\begin{align} \label{est-R}
& \sup_{s \in [0, T)} \mathbb E[M_s \log M_s]
\le \frac{\lambda \|G^{-1}\|_\infty^2}{e^{2 \lambda T}-1}\|x-y\|_p^2,
\quad x, y \in L^p.
\end{align}
Consequently, $M_T:=\lim_{s \uparrow T} M_s$ exists and $(M_s)_{s \in [0, T]}$ is a martingale.
\end{lm}

\begin{proof}
Let $s \in [0, T)$ be fixed.
By the construction \eqref{v}, we can rewrite Eq. \eqref{see} and Eq. \eqref{rd-y} on $[0, s]$ as
\begin{align}
{\rm d}X_t &=( AX_t+F(X_t)-\gamma_t^{-1}(X_t-Y_t) ) {\rm d}t +G {\rm d}{\widetilde W}_t, \label{rd-x+}   \\
{\rm d}Y_t &=(AY_t+F(Y_t)) {\rm d}t+G {\rm d}{\widetilde W}_t, \label{rd-y+}
\end{align}
with initial values $X_0=x$ and $Y_0=y$, respectively.
Then Eq. \eqref{x-y-p} about $X-Y$ also holds ${\mathbb Q}_s$-a.s.
Therefore, the pathwise estimate \eqref{x-y} is valid ${\mathbb Q}_s$-a.s., which in combination with the equality \eqref{M+} and the estimate \eqref{est-v} implies that
\begin{align} \label{log-m}
\log M_s & =-\int_0^s (v_r, {\rm d}{\widetilde W}_r)
+\frac12 \int_0^s \|v_r\|^2 {\rm d}r  \nonumber \\
& \le -\int_0^s (v_r, {\rm d}{\widetilde W}_r)
+\frac{\|G^{-1}\|_\infty^2 \|x-y\|_p^2}{2 \gamma_0}, \quad {\mathbb Q}_s\text{-a.s.}
\end{align}
Taking into account the fact in Lemma \ref{lm-R-} that $(\widetilde W_t)_{t \in [0, s]}$ is a $U$-valued cylindrical Wiener process under ${\mathbb Q}_s$, we arrive at
\begin{align*}
& \mathbb E[M_s \log M_s]
=\mathbb E_s \log M_s
\le \frac{\|G^{-1}\|_\infty^2 \|x-y\|_p^2}{2 \gamma_0},
\end{align*}
and thus obtain \eqref{est-R}, noting that $\gamma_0$ is given in \eqref{gt} with $t=0$, where $\mathbb E_s$ denotes the expectation with respect to ${\mathbb Q}_s$.
By the martingale convergence theorem, $M_T:=\lim_{s \uparrow T} M_s$ exists and $(M_t)_{t \in [0, T]}$ is a well-defined martingale.
\end{proof}

Lemmas \ref{lm-R-} and \ref{lm-R} ensures that $(\widetilde W_t)_{t \in [0, T]}$ is a $U$-valued cylindrical Wiener proces under the probability measure ${\mathbb Q}:=M_T \mathbb P$ and
\begin{align} \label{est-R-T}
& \sup_{s \in [0, T]} \mathbb E[M_s \log M_s]
\le \frac{\lambda \|G^{-1}\|_\infty^2}{e^{2 \lambda T}-1}\|x-y\|_p^2,
\quad x, y \in L^p.
\end{align}
Then by \eqref{see} and \eqref{rd-y+}, the coupling $(X_t, Y_t)$ is well constructed under ${\mathbb Q}$ for $t \in [0, T]$.

\begin{lm} \label{lm-R+}
Let Assumptions \ref{ap-f}, \ref{ap-wa}, and \ref{ap-ell} hold.
For any $x, y \in L^p$ and $\mathbf s>1$,
\begin{align} \label{est-R+}
& \sup_{s \in [0, T]} \mathbb E M_s^{\frac \mathbf s{\mathbf s-1}}
\le \exp\Big(\frac{\mathbf s \lambda \|G^{-1}\|_\infty^2 \|x-y\|_p^2}{(\mathbf s-1)^2 (e^{2 \lambda T}-1)}\Big).
\end{align}
\end{lm}

\begin{proof}
Let $s \in [0, T]$.
Denote by $v^\mathbf s_r:=-\frac1{\mathbf s-1} v_r$ for $r \in [0, s] \subset [0, T]$.
The representation \eqref{M+} and the pathwise estimate \eqref{x-y} with $\gamma_0$ given in \eqref{gt} yield that
\begin{align*}
M_s^{\frac1{\mathbf s-1}}
& =\exp \Big( -\frac1{\mathbf s-1}\int_0^s \<v_r, {\rm d} {\widetilde W}_r\> + \frac1{2(\mathbf s-1)} \int_0^s \|v_r\|^2 {\rm d}r \Big)  \\
& =\exp \Big( \int_0^s \<v^\mathbf s_r, {\rm d}{\widetilde W}_r \>
- \frac12 \int_0^s \|v^\mathbf s_r\|^2 {\rm d}r \Big) \\
& \qquad \times \exp \Big(\frac\mathbf s{2(\mathbf s-1)^2} \int_0^s \| v_r \|^2 {\rm d}r \Big) \\
& \le {\widetilde M}_s \exp\Big(\frac{\mathbf s \lambda \|G^{-1}\|_\infty^2 \|x-y\|_p^2}{(\mathbf s-1)^2 (e^{2 \lambda T}-1)}\Big), \quad {\mathbb Q}_s\text{-a.s.}
\end{align*}
where ${\widetilde M}_s:=\exp( \int_0^s \<v^\mathbf s_r, {\rm d}{\widetilde W}_r \>
- \frac12 \int_0^s \|v^\mathbf s_r\|^2 {\rm d}r)$.
It follows that
\begin{align*}
\mathbb E M_s^{\frac \mathbf s{\mathbf s-1}}
=\mathbb E_s M_s^{\frac1{\mathbf s-1}}
\le \exp\Big(\frac{\mathbf s \lambda \|G^{-1}\|_\infty^2 \|x-y\|_p^2}{(\mathbf s-1)^2 (e^{2 \lambda T}-1)}\Big) \mathbb E_s {\widetilde M}_s.
\end{align*}
Taking into account the fact that $(\widetilde W_t)_{t \in [0, T]}$ is a $U$-valued cylindrical Wiener process under ${\mathbb Q}$, which shows that $({\widetilde M}_t)_{t \in [0, T]}$ is a martingale under ${\mathbb Q}$, we have
$\mathbb E_s {\widetilde M}_s=\mathbb E_s {\widetilde M}_0=1$
and obtain \eqref{est-R+}.
\end{proof}

\section{Harnack Inequalities and Ergodicity}
\label{sec4}

In the last section, we derive Harnack inequalities and the ergodicity for the Markov semigroup $(P_t)_{t \ge 0}$, in the following three parts.
At the first two parts, we give the proof of our main results, Theorems \ref{main1} and \ref{main2}, espectively.
Other applications, including several estimates for the density of $(P_t)_{t \ge 0}$,  are also derived.

\subsection{Harnack inequalities}

We begin with the following Harnack inequalities.

\begin{tm} \label{tm-har}
Let Assumptions \ref{ap-f}, \ref{ap-wa}, and \ref{ap-ell} hold.
Then \eqref{har-log} and \eqref{har-pow} hold for any $T > 0$, $\mathbf s >1$, $x, y \in L^p$, and $\phi \in \mathcal B^+_b(L^p)$.
\end{tm}

\begin{proof}
We first show that $X_T=Y_T$ ${\mathbb Q}$-a.s.
From Lemma \ref{lm-R}, $(M_t)_{t \in [0, T]}$ is a uniformly integrable martingale and $({\widetilde W}_t)_{[0, T]}$ is a cylindrical Wiener process under the probability measure ${\mathbb Q}$.
So $Y_t$ can be solved up to time $T$.
Let
\begin{align*}
\tau:=\inf\{t \in [0, T]: \ X_t=Y_t\}
\quad \text{with} \quad \inf \emptyset : =\infty.
\end{align*}

Suppose that for some $\omega \in \Omega$ such that $\tau(\omega)>T$, then the continuity of the process $X-Y$, in Lemma \ref{lm-lq} and Remark \ref{rk-Y}, yields
\begin{align*}
\inf_{t \in [0, T]} \|X_t-Y_t\|_p^2(\omega)>0.
\end{align*}
By the divergence relation \eqref{gt>},
\begin{align*}
\int_0^T \frac{\|X_t-Y_t\|_p^2(\omega)}{\gamma_t^2} {\rm d}t=\infty
\end{align*}
holds on the set $(\tau>T):=\{\omega: \tau(\omega)>T\}$.
But according to the pathwise estimate \eqref{x-y} which holds ${\mathbb Q}$-a.s.,
\begin{align*}
\mathbb E_{{\mathbb Q}} \int_0^T \frac{\|X_t-Y_t\|_p^2(\omega)}{\gamma_t^2} {\rm d}t
\le \frac{\|x-y\|_p^2}{\gamma_0}<\infty,
\end{align*}
where $\mathbb E_{{\mathbb Q}}$ denotes the expectation with respect to ${\mathbb Q}$.
It follows from the above two estimates that ${\mathbb Q}(\tau>T)=0$, i.e., $\tau \le T$ ${\mathbb Q}$-a.s.
By the definition of $\tau$, $X_T=Y_T$ ${\mathbb Q}$-a.s.

Therefore, we get a coupling $(X,Y)$ by the change of measure, with changed probability ${\mathbb Q}= M_T \mathbb P$, such that $X_T=Y_T$ ${\mathbb Q}$-a.s.
Consequently, the inequalities \eqref{har-log} and \eqref{har-pow} follow from the following known inequalities (see, e.g., \cite[Theorem 1.1.1]{Wan13}):
\begin{align*}
P_T \log \phi(y) & \le \log P_T \phi(x)+\mathbb E [M_t \log M_t], \\
(P_T \phi(y))^\mathbf s & \le (P_T \phi^\mathbf s(x)) (\mathbb E M_t^{\frac \mathbf s{\mathbf s-1}})^{\mathbf s-1},
\end{align*}
and the estimations \eqref{est-R} and \eqref{est-R+}, in Lemmas \ref{lm-R} and \ref{lm-R+}, respectively.
\end{proof}

The log-Harnack inequality \eqref{har-log} imply the following gradient estimate and regularity properties for the Markov semigroup $(P_t)_{t \ge 0}$.

\begin{cor} \label{cor-fel}
Let Assumptions \ref{ap-f}, \ref{ap-wa}, and \ref{ap-ell} hold.
For any $T>0$ and $\phi \in \mathcal B_b(L^p)$,
\begin{align}\label{est-gra}
\|D P_T \phi\| \le \sqrt{\frac{2 \lambda \|G^{-1}\|_\infty^2}{e^{2 \lambda T}-1}} \sqrt{P_T \phi^2-(P_T \phi)^2}.
\end{align}
Consequently, $(P_t)_{t \ge 0}$ is strong Feller and has at most one invariant measure, and if it has one, the density of $(P_t)_{t \ge 0}$ with respect to the invariant measure is strictly positive.
\end{cor}

\begin{proof}
The gradient estimate \eqref{est-gra} and the uniqueness of the invariant measure for $(P_t)_{t \ge 0}$ with a strictly positive density, if it exists, are direct consequence of the log-Harnack inequality \eqref{har-log}, see Proposition 1.3.8 and Theorem 1.4.1 in \cite{Wan13}, respectively.
Finally, the strong Feller property of $(P_t)_{t \ge 0}$ follows easily from the gradient estimate \eqref{est-gra}.
\end{proof}

\begin{rk}
The uniqueness of the invariant measure, if it exists, for 1D stochastic heat equation (Eq. \eqref{see} with Lipschitz coefficients) driven by white noise on $L^p(0, 1)$ with $p>4$ was shown in \cite{BR16(DCDS)}.
So Corollary \ref{cor-fel} can be seen as filling the gap for $p \in (2, 4]$ in the additive white noise case.
\end{rk}

\begin{rk}
Under the conditions in Corollary \ref{cor-fel}, there exist constants $C, T_0 >0$ such that
\begin{align} \label{tv}
\|\mathcal L(X_t^x)-\mathcal L(X_t^y)\|_{\rm TV} \le C e^{-\lambda t} \|x-y\|_p,
\quad t \ge T_0,  \ x, y \in L^p,
\end{align}
where $\mathcal L(X)$ denotes the distribution of $X$ on $L^p$, $\lambda$ is given in \eqref{lam}, and $\|\cdot \|_{\rm TV}$ denotes the total variation norm betweem two signed measures, i.e.,
$\|\mu-\nu\|_{\rm TV}:=\sup_{\|\phi\|_\infty \le 1}|\int_{L^p} \phi {\rm d}\mu-\int_{L^p} \phi {\rm d}\nu|$ for two signed measures $\mu$ and $\nu$.
\end{rk}

\iffalse
\begin{proof}
Let $T_0$ be any positive number.
By the gradient estimate \eqref{est-gra}, we have
\begin{align*}
|P_{T_0} f(x)-P_{T_0} f(y)|^2
& \le \|D P_{T_0} f\|^2 \|x-y\|_p^2
\le \frac{2 \lambda \|G^{-1}\|_\infty^2}{e^{2 \lambda {T_0}}-1} \|f\|_\infty^2 \|x-y\|_p^2.
\end{align*}
Combining this with the Markov property, we obtain
\begin{align*}
& \|\mathcal L(X_{t+t_0}^x)-\mathcal L(X_{t+t_0}^y)\|_{\rm TV} \\
& \le \sup_{\|f\|_\infty \le 1} \mathbb E |P_{T_0}f(X_t^x)-P_{T_0}f(X_t^y)| \\
& \le \sqrt{\frac{2 \lambda \|G^{-1}\|_\infty^2}{e^{2 \lambda {T_0}}-1} } \mathbb E \|X_t^x-X_t^y\|_p,
\quad t \ge 0.
\end{align*}
Then \eqref{tv} with $C=\sqrt{\frac{2 \lambda \|G^{-1}\|_\infty^2}{e^{2 \lambda {T_0}}-1} }$ follows from \eqref{con}.
\end{proof}
\fi

\begin{proof}[Proof of Theorem \ref{main1}]
Theorem \ref{main1} follows from Theorem \ref{tm-har} and Corollary \ref{cor-fel}.
\end{proof}

\subsection{Ergodicity}

In this part, we show the existence of an invariant measure for the Markov semigroup $(P_t)_{t \ge 0}$.
In combination with the uniqueness of the invariant measure, as shown in Corollary \ref{cor-fel}, we derive the existence of a unique and thus ergodic invariant measure.
We also note that \cite[Theorem 6.1]{BG99(SPA)} used the factorization approach to obtain the existence of an invariant measure in $L^p$ with $p \ge 2$ under the martingale solution framework.

\begin{tm} \label{tm-erg}
Let Assumptions \ref{ap-f}, \ref{ap-wa}, and \ref{ap-wa+} hold.
Assume that $q>2$ and \eqref{d} holds.
Then $(P_t)_{t \ge 0}$ has an invariant measure $\mu$ with full support on $L^p$ such that $\mu( \|\cdot\|_{q+p-2}^{q+p-2})<\infty$.
Assume furthermore that Assumption \ref{ap-ell} holds, then $\mu$ is the unique and thus ergodic invariant measure of $(P_t)_{t \ge 0}$.
\end{tm}

\begin{proof}
The uniqueness of the invariant measure and the strong Feller Markov property for $(P_t)_{t \ge 0}$ have been shown in Corollary \ref{cor-fel}.
Thus, to show the existence of an invariant measure, by Krylov--Bogoliubov theorem, it suffices to verify the tightness of the sequence of probability measures $(\mu_n)_{n \in\mathbb N_+}$ defined by
\begin{align} \label{mun}
\mu_n:=\frac1n \int_0^n \delta_0 P_t {\rm d}t, \quad n \in \mathbb N_+,
\end{align}
where $\delta_0 P_t$ is the distribution of $X_t^0$, the solution of Eq. \eqref{see} with the initial datum $X_0=0$.

It follows from the relation $X=Z+W_A$, the estimate \eqref{Zp} with $x=0$, and Young inequality that
\begin{align} \label{est-tig}
& \|X^0_t\|_p^p
\le 2^{p-1} \|Z^0_t\|_p^p+2^{p-1}\|W_A(t)\|_p^p  \\
& \le  C \int_0^t (1+\|W_A(r)\|_{q+p-2}^{q+p-2}) {\rm d} r
+2^{p-1} p L_f \int_0^t \|Z^0_r\|_p^p  {\rm d} r  \nonumber \\
& \quad - 2^{p-1} p \theta_1 \int_0^t \|Z^0_r\|_{q+p-2}^{q+p-2} {\rm d} r +2^{p-1} \|W_A(t)\|_p^p \nonumber \\
& \le  C \int_0^t (1+\|W_A(r)\|_{q+p-2}^{q+p-2}) {\rm d} r
- \theta_4 \int_0^t \|Z^0_r\|_{q+p-2}^{q+p-2} {\rm d} r
 +2^{p-1}\|W_A(t)\|_p^p \nonumber \\
 & \le  C \int_0^t (1+\|W_A(r)\|_{q+p-2}^{q+p-2}) {\rm d} r
- \theta_5 \int_0^t \|X^0_r\|_{q+p-2}^{q+p-2} {\rm d} r
 +2^{p-1}\|W_A(t)\|_p^p,  \nonumber
\end{align}
for some constants $\theta_4, \theta_5>0$, where we have used the elementary inequality $|\xi-\eta|^\mathbf r \ge 2^{1-r} \xi^r-\eta^r$ for $\xi, \eta \ge 0$ and $\mathbf r \ge 1$, in the last inequality.
Then we have
\begin{align*}
& \theta_5 \int_0^t \|X^0_r\|_{q+p-2}^{q+p-2} {\rm d} r
\le C \int_0^t (1+\|W_A(r)\|_{q+p-2}^{q+p-2}) {\rm d} r
+ 2^{p-1} \|W_A(t)\|_p^p.
\end{align*}
The above estimate, in combination with the condition \eqref{wa1}, yields that there exists a constant $C$ such that for all $n \ge 1$,
\begin{align} \label{est-tig1}
& \mu_n(\|\cdot\|_{q+p-2}^{q+p-2})
=\frac1n \int_0^n \mathbb E \|X^0_r\|_{q+p-2}^{q+p-2} {\rm d} r
\nonumber \\
& \le \frac C{\theta_5}
\Big( 1+\frac{\mathbb E \|W_A(n)\|_p^p}{n}+
\frac1n \int_0^n \mathbb E \|W_A(r)\|_{q+p-2}^{q+p-2} {\rm d} r \Big)
\le C.
\end{align}

It follows from the ultracontractivity \eqref{ult}, with $\mathbf r=p$ and $\mathbf s=\frac{q+p-2}{q-1}$, and Young convolution inequality that
\begin{align*}
& \int_0^n \Big\| \int_0^t S_{t-r} F(X^0_r) {\rm d}r \Big\|_{\beta, p} {\rm d}t   \\
 & \le C \int_0^n \int_0^t e^{-\lambda_1 (t-r)}
  (t-r)^{-\alpha}(1+\|X^0_r\|^{q-1}_{q+p-2}) {\rm d}r {\rm d}t  \\
   & \le C \Big(\int_0^n e^{-\lambda_1 t} t^{-\alpha}{\rm d}t \Big)
\Big(\int_0^n  (1+\|X^0_t\|^{q-1}_{q+p-2}) {\rm d}t \Big),
\end{align*}
where $\alpha=\frac\beta 2+\frac{d(p-1)(q-2)}{2p(q+p-2)} \in (0, 1)$ provided that $\beta>0$ is sufficiently small, since $d<\frac{2p(q+p-2)}{(p-1)(q-2)}$.
The fact that
\begin{align*}
& \sup_{n \ge 1}\Big(\int_0^n e^{-\lambda_1 t} t^{-\alpha}{\rm d}t \Big)
\le \int_0^\infty e^{-\lambda_1 t} t^{-\alpha}{\rm d}t
<\infty,
\end{align*}
for all $\lambda_1>0$ and $\alpha \in (0, 1)$, and Young inequality imply that
\begin{align*}
& \int_0^n \Big\| \int_0^t S_{t-r} F(X^0_r) {\rm d}r \Big\|_{\beta, p} {\rm d}t
\le C \int_0^n  (1+\|X^0_t\|^{q+p-2}_{q+p-2}) {\rm d}t.
\end{align*}
By Fubini theorem, the estimate \eqref{est-tig1}, and the condition \eqref{wa+}, we arrive at
\begin{align} \label{est-tig2}
& \mu_n(\|\cdot\|_{\beta,p})
=\frac1n \int_0^n \mathbb E \|X^0_r\|_{\beta,p} {\rm d} r \\
& \le \frac1n \mathbb E \int_0^n \Big\| \int_0^t S_{t-r} F(X^0_r) {\rm d}r \Big\|_{\beta, p} {\rm d} r
+\frac1n \int_0^n \mathbb E \|W_A(t)\|_{\beta,p} {\rm d} r \nonumber \\
& \le \frac Cn \int_0^n  (1+ \mathbb E \|X^0_t\|^{q+p-2}_{q+p-2}) {\rm d}t
+\frac1n \int_0^n \mathbb E \|W_A(t)\|_{\beta,p} {\rm d} r
\le C<\infty, \nonumber
\end{align}
for all $n \ge 1$ and $\beta<(1-\frac{d(p-1)(q-2)}{2p(q+p-2)}) \wedge \beta_0$.
For any fixed $p \ge 2$, we take $\beta<(1-\frac{d(p-1)(q-2)}{2p(q+p-2)}) \wedge \beta_0 \wedge \frac d p$, so that the embedding $W_0^{\beta,p} \subset L^p$ in \eqref{emb} is compact.
Consequently, the above estimate \eqref{est-tig2} shows that
$\{ u \in L^p: \ \|u\|_{\beta, p} \le N \}$ is relatively compact in $L^p$ for any $N>0$, and thus $(\mu_n)_{n \in \mathbb N_+}$ is tight.
This shows the existence of an invariant measure, denoted by $\mu$, of $(P_t)_{t \ge 0}$.
\iffalse
By the Prohorov theorem, there exists a probability measure $\mu$ and a subsequence  $(\mu_{n_k})_{k \in \mathbb N_+}$ such that  $\mu_{n_k} \rightarrow \mu$ weakly as $k \rightarrow \infty$.
One can easily show that $\mu$ is an invariant probability measure for
$(P_t)_{t \ge 0}$.
\fi

To show that the invariant measure $\mu$ has full support on $L^p$, let us choose $\mathbf s=2$, $\phi =\chi_\Gamma$, in \eqref{har-pow}, with $\Gamma$ being a Borel set in $L^p$, and get
\begin{align*}
& (P_T \chi_\Gamma(x))^2 \int_{L^p} \exp\Big(-\frac{2 \lambda \|G^{-1}\|_\infty^2}{e^{2 \lambda T}-1} \|x-y\|_p^2\Big) \mu({\rm d} y) \\
& \le \int_H P_T \chi_\Gamma(y) \mu({\rm d} y)
=\int_H \chi_\Gamma(y) \mu({\rm d} y)=\mu(\Gamma),
\quad T>0, \ x \in L^p.
\end{align*}
This shows that the transition kernel of $(P_t)_{t \ge 0}$ is absolutely continuous with respect to $\mu$ so that it has a density $p_T(x, y)$.
Suppose that ${\rm supp} \ \mu\neq L^p$, then there exist $x_0 \in L^p$ and $r>0$ such that $\mu(B(x_0, r))=0$, where $B(x_0, r)$ is a ball in $L^p$ with radius $r$ and center $x_0$.
Then $p_T(x_0, B(x_0, r))=0$ and
$\mathbb P(\|X_T^{x_0}-x_0\|_p \le r)=0$ for all $T>0$.
This contradicts with the fact that $X_T^{x_0}$ is a continuous process on $L^p$ as shown in Lemma \ref{lm-lq}.

Similarly to \eqref{est-tig1}, we have (with $n=1$ and $X_0=x$)
\begin{align*}
\int_0^1 P_t \|\cdot\|_{q+p-2}^{q+p-2}(x) {\rm d}t
=\int_0^1 \mathbb E \|X^x_t\|_{q+p-2}^{q+p-2} {\rm d} t
\le C (1+\|x\|_p^p).
\end{align*}
Integrating on $L^p$ with respect to the invariant measure $\mu$ and using Fubini theorem, we obtain
\begin{align*}
\mu( \|\cdot\|_{q+p-2}^{q+p-2})
=\int_0^1 \int_H P_t \|\cdot\|_{q+p-2}^{q+p-2}(x)\mu({\rm d}x) {\rm d}t
\le C (1+\mu(\|\cdot\|_p^p))<\infty.
\end{align*}
This shows that $\mu( \|\cdot\|_{q+p-2}^{q+p-2})<\infty$ and completes the proof.
\end{proof}

\begin{rk}
In the case $q>2=p$, the condition \eqref{d} is equivalent to $d<4+8/(q-2)$, which will be always valid in $d=1,2,3$-dimensional cases.
\end{rk}

\begin{rk} \label{rk-erg}
Under Assumptions \ref{ap-f} and \ref{ap-wa}, if $\lambda$ defined in \eqref{lam} is positive, one can use the standard remote control method to show that $(P_t)_{t \ge 0}$ has an invariant measure $\mu$, once we derive similar estimates as \eqref{con} and \eqref{Zp}, without the restriction \eqref{d}; see, e.g., \cite[Theorem 6.3.2]{PZ96}.
In this case, $\mu$ also has full support on $L^p$ such that $\mu( \|\cdot\|_{p+q-2}^{p+q-2})<\infty$ holds.
\iffalse
Indeed, one can introduce another $U$-valued cylindrical Wiener process
$(V_t)_{t \ge 0}$, independent of $(W_t)_{t \ge 0}$, and define
${\widehat W}_t = W_t$ when $t \ge 0$,
${\widehat W}_t = V_t$ when $t \le 0$, and
${\widehat \mathscr F}_t=\sigma({\widehat W}_s, \ s \le t)$ for $t \in \mathbb R$.
Denote by $X_{t,s}^x$ the solution of
\begin{align}\label{see+}
{\rm d}X_t &=(A X_t+F(X_t)) {\rm d}t+ G{\rm d}{\widehat W}_t,
\quad t \ge s,
\end{align}
with $X_s=x$.
Then similar idea deriving the estimate \eqref{est-tig} yields that
\begin{align*}
\|X^x_{t, s} \|_p^p
& \le  \|x\|_p^p + C \int_s^t (1+\|{\widehat W}_A(r,s)\|_p^p) {\rm d} r
- \theta_5 \int_s^t \|X^0_r\|_p^p {\rm d} r \\
& \quad  +2^{p-1}\|{\widehat W}_A(t,s)\|_p^p,  \nonumber
\end{align*}
where ${\widehat W}_A(t,s)=\int_s^t S_{r-s} {\rm d}{\widehat W}_s$.
This, in combination with \eqref{wa1}, implies the following moments' estimation for $X_{t,s}^x$:
\begin{align*}
& \mathbb E \|X^x_{t, s} \|_p^p
\le  (C+\|x\|_p^p) e^{-\theta_5 (t-s)},   \quad t \ge s.
\end{align*}
Using the same arguments to derive the estimates \eqref{Zp} and \eqref{est-tig}, we obtain
\begin{align*}
\mathbb E \|X_{0, -t_1}^x-X_{0, -t_2}^x\|_p^p
\le  e^{-\theta_6 t_1} \mathbb E \|X^x_{-t_1, -t_2} \|_p^p
\le  (C+\|x\|_p^p) e^{-\theta_7 t_2},
\end{align*}
for constants $\theta_6, \theta_7>0$ and all $t >-t_1>-t_2>0$.
It follows that the sequence of random variables $(X_{0, -t}^x)_{t \ge 0}$ is a Cauchy sequence in $L^p(\Omega; L^p)$ as $t \rightarrow \infty$, and therefore it is convergent to a random variable $\eta \in L^p(\Omega; L^p)$, which is independent of $x$.
It is not difficult to show that $\eta$ is an invariant measure of $(P_t)_{t \ge 0}$.
This completes the proof of the existence of an invariant measure when $q>2$ and $\lambda>0$.
\fi
\end{rk}

\iffalse
\begin{rk}
As usual, the existence of an invariant measure $\mu$ of $(P_t)_{t \ge 0}$ is also hold without the non-degenerate Assumption \ref{ap-ell}.
Indeed, following the proof in Theorem \ref{tm-erg}, one just need to show the Feller property of $(P_t)_{t \ge 0}$, which is ensured by the pathwise continuity \eqref{con}.
\end{rk}
\fi

\begin{rk}
Under the conditions of Theorem \ref{tm-erg} or Remark \ref{rk-erg}, $(P_t)_{t \ge 0}$ has a unique invariant measure $\mu$ with full support on $L^p$, which shows that $(P_t)_{t \ge 0}$ is irreducible, i.e., $P_T\chi_{\Gamma}(x)>0$ for arbitrary non-empty open set $\Gamma \subset L^p$, $x \in L^p$, and $T>0$.
Indeed, the power-Harnack inequality \eqref{har-pow} with $f=\chi_\Gamma$ yields that
\begin{align} \label{irr}
 (P_T \chi_\Gamma(y))^\mathbf s \le P_T \chi_\Gamma(x)
\exp\Big(\frac{\mathbf s \lambda \|G^{-1}\|_\infty^2 \|x-y\|_p^2}{(\mathbf s-1) (e^{2 \lambda T}-1)}\Big),
\quad y \in L^p.
\end{align}
The facts that $\mu$ is $P_T$-invariant and has full support on $L^p$ imply
\begin{align*}
\int_{L^p} P_T \chi_\Gamma(y) \mu({\rm d}y)
\le \int_{L^p} \chi_\Gamma(y) \mu({\rm d}y)
=\mu(\Gamma)>0,
\end{align*}
which shows that there is a $y \in L^p$ such that $P_T \chi_\Gamma(y)>0$.
Then \eqref{irr} yields that $P_T \chi_\Gamma(x)>0$ for all $x \in L^p$, so that the  irreducibility holds.
\end{rk}

\begin{proof}[Proof of Theorem \ref{main2}]
Theorem \ref{main2} follows from Theorem \ref{tm-erg} and Remark \ref{rk-erg}.
\end{proof}

\subsection{Estimates of density}

Finally, we use the Harnack inequalities \eqref{har-log} and \eqref{har-pow} to derive an estimate of the density, denoted by $p_T(x, y)$, with respect to the invariant measure $\mu$ of $(P_t)_{t \ge 0}$.

\begin{cor}  \label{cor-den}
Let Assumptions \ref{ap-f} and \ref{ap-wa} hold.
Assume that $q>2$ such that \eqref{d} hold and Assumption \ref{ap-wa+} hold, or $\lambda$ defined in \eqref{lam} is positive.
Then for all $T>0$, $x \in L^p$, and $\mathbf s >1$,
\begin{align}\label{den}
& \|p_T (x, \cdot)\|_{L^\mathbf s(\mu)}
\le \Big( \int_H \exp\Big(-\frac{\mathbf s \lambda \|G^{-1}\|_\infty^2 \|x-y\|_p^2}{e^{2 \lambda T}-1} \Big) \mu({\rm d}y)  \Big)^{-\frac{\mathbf s-1}\mathbf s}.
\end{align}
\end{cor}

\begin{proof}
Using \eqref{har-pow} with $\mathbf s$ replaced by $\frac\mathbf s{\mathbf s-1}$, and noting that $\mu$ is $(P_t)_{t \ge 0}$-invariant, we obtain
\begin{align*}
& \|p_T (x, \cdot)\|_{L^\mathbf s(\mu)} \\
& =\sup \{\< p_T (x, \cdot), \phi\>_\mu: \ \phi \in B_b^+(L^p), ~ \mu(\phi^\frac\mathbf s{\mathbf s-1}) \le 1 \} \\
& =\sup \{P_T \phi(x) : \ \phi \in B_b^+(L^p), \ \mu(\phi^\frac\mathbf s{\mathbf s-1}) \le 1\} \\
& \le \Big( \int_H \phi^\frac\mathbf s{\mathbf s-1}(y) \mu({\rm d}y) \Big)^\frac{\mathbf s-1}\mathbf s
\Big( \int_H \exp\Big(-\frac{\mathbf s \lambda \|G^{-1}\|_\infty^2 \|x-y\|_p^2}{e^{2 \lambda T}-1} \Big) \mu({\rm d}y) \Big)^\frac{\mathbf s-1}\mathbf s \\
& = \mu(\phi^\frac\mathbf s{\mathbf s-1})
\Big( \int_H \exp\Big(-\frac{\mathbf s \lambda \|G^{-1}\|_\infty^2 \|x-y\|_p^2}{e^{2 \lambda T}-1} \Big) \mu({\rm d}y) \Big)^\frac{\mathbf s-1}\mathbf s.
\end{align*}
This shows the density estimate \eqref{den}.
\end{proof}

\begin{rk}
According to \cite[Proposition 2.4]{Wan10(JMPA)},
the log-Harnack inequality \eqref{har-log} and the power-Harnack inequality \eqref{har-pow} are equivalent to the following two heat kernel inequalities, respectively, provided $P_T$ have a strictly positive density $p_T(x, y)$ with respect to a Radon measure of $P_T$:
\begin{align*}
\int_{L^p} p_T(x, z) \log \frac{p_T(x, z)}{p_T(y, z)} \mu({\rm d}z)
& \le \frac{\lambda \|G^{-1}\|_\infty^2}{e^{2 \lambda T}-1} \|x-y\|_p^2,  \\
\int_{L^p} p_T(x, z) \Big(\frac{p_T(x, z)}{p_T(y, z)}\Big)^\frac1{\mathbf s-1} \mu({\rm d}z)
& \le \frac{\mathbf s \lambda\|G^{-1}\|_\infty^2}{(\mathbf s-1)^2 (1-e^{-2 \lambda T)}}\|x-y\|_p^2.
\end{align*}
Under the conditions in Corollary \ref{cor-den}, $(P_t)_{t \ge 0}$ has a unique invariant measure $\mu$ such that $p_T(x, y)$ is strictly positive.
Then the above two heat kernel inequalities are direct consequence of Theorem \ref{tm-har}, Theorem \ref{tm-erg}, and Remark \ref{rk-erg}.
\end{rk}

\iffalse
Finally, we study the contractivity property of $(P_t)_{t \ge 0}$.

\begin{cor} \label{cor-hyp}
Let Assumptions \ref{ap-f} and \ref{ap-ell}, and the condition \eqref{wa2} hold for a sufficiently large $\epsilon_0>0$.
Assume that $\lambda$ defined in \eqref{lam} is positive, then $(P_t)_{t \ge 0}$ is hypercontractive.
Consequently, $(P_t)_{t \ge 0}$ is compact in $L^2(\mu)$ for large $t>0$, and there exist constants $C, \lambda_0 >0$ such that \eqref{ent} and \eqref{spe-gap} hold for all $t \ge 0$.
\end{cor}

\begin{proof}
According to \cite[Theorem 2.1]{Wan17(JFA)}, to show the hypercontractivity for some $t>0$ and compactness in $L^2(\mu)$ for large $t>0$ of $(P_t)_{t \ge 0}$ such that \eqref{ent} and \eqref{spe-gap} hold, it suffices to show that there exists a measurable function $\phi: [0, \infty) \rightarrow (0, \infty)$ with $\lim_{t \uparrow \infty}\phi(t)=0$ and constants $T_0, C_0, \epsilon>0$ such that the following conditions hold:
\begin{align*}
\|X_t^x-X_t^y\|_p \le \phi(t) \|x-y\|_p,
& \quad t \ge 0, \ x, y \in L^p; \\
(P_{T_0} f(y))^2 \le P_{T_0} f^2(x)
\exp(C_0 \|x-y\|_p^2), & \quad x, y \in L^p, \ \phi \in \mathcal B_b(L^p); \\
\mu(\exp(\epsilon \|\cdot\|_p^2))<\infty. &
\end{align*}

The first coupling property follows from \eqref{con} with $\phi(t)=e^{- \lambda t}$.
The second condition is ensured by the power-Harnack inequality \eqref{har-pow} with $\mathbf s=2$ and $C_0=2 \lambda \|G^{-1}\|_\infty^2 (e^{2 \lambda T_0}-1)^{-1}$ for any $T_0>0$ and $\phi \in \mathcal B^+_b(L^p)$ (and,  clearly, for all $\phi \in \mathcal B_b(L^p)$).
Finally, by the dominated convergence theorem and the fact $\mu_t^0 \rightarrow \mu$ weakly as $t \rightarrow \infty$, we have
\begin{align} \label{mu-e}
\mu(\exp(\epsilon \|\cdot\|_p^2))
& =\lim_{N \rightarrow \infty} \mu(N \wedge \exp(\epsilon \|\cdot\|_p^2)) \nonumber \\
& =\lim_{N \rightarrow \infty} \lim_{t \rightarrow \infty}
\mu_t^0(N \wedge \exp(\epsilon \|\cdot\|_p^2))  \nonumber  \\
& =\lim_{N \rightarrow \infty} \lim_{t \rightarrow \infty} \mathbb E [N \wedge \exp(\epsilon \|X_t^0\|_p^2)]  \nonumber \\
& \le \sup_{t \ge 0} \mathbb E \exp(\epsilon \|X_t^0\|_p^2).
\end{align}
Then the last concentration property follows from \eqref{mu-e} and \eqref{exp}.
\end{proof}

In the case $q>2$, even $\lambda$ defined in \eqref{lam} is non-positive, the existence of an invariant measure for $(P_t)_{t \ge 0}$ has been shown in Theorem \ref{tm-erg}.
In this case, we have the following supercontractivity of $(P_t)_{t \ge 0}$.

\begin{cor} \label{cor-sup}
Let $q>2$ and Assumptions \ref{ap-f} and \ref{ap-ell} hold.
Assume that \eqref{wa2} holds for all $\epsilon_0>0$, then $(P_t)_{t \ge 0}$ is supercontractive.
\end{cor}

\begin{proof}
The power-Harnack inequality \eqref{har-pow} with $\mathbf s=2$ yields that
\begin{align*}
& (P_T \phi(y))^2 \le P_T \phi^2(x)
\exp\Big(\frac{2 \lambda \|G^{-1}\|_\infty^2}{e^{2 \lambda T}-1}
\|x-y\|_p^2\Big), \quad T>0, \ \phi \in \mathcal B_b(L^p).
\end{align*}
Then for any $\phi$ with $\mu(\phi^2)=1$, we have
\begin{align*}
& (P_T \phi(y))^2  \int_{L^p} \exp\Big(-\frac{2 \lambda \|G^{-1}\|_\infty^2}{e^{2 \lambda T}-1} \|x-y\|_p^2 \Big) \mu({\rm d}x) \\
& \le \int_{L^p} P_T \phi^2(x) \mu({\rm d}x)
=\int_{L^p} f^2(x) \mu({\rm d}y)
=\mu(\phi^2)=1.
\end{align*}
It is clear that
\begin{align*}
& \int_{L^p} \exp\Big(-\frac{2 \lambda \|G^{-1}\|_\infty^2}{e^{2 \lambda T}-1} \|x-y\|_p^2 \Big) \mu({\rm d}x)  \\
& \ge \int_{B(0,1)} \exp\Big(-\frac{2 \lambda \|G^{-1}\|_\infty^2}{e^{2 \lambda T}-1} (1+\|y\|_p)^2 \Big) \mu({\rm d}x) \\
& \ge \mu(B(0,1)) \exp\Big(-\frac{2 \lambda \|G^{-1}\|_\infty^2}{e^{2 \lambda T}-1} (1+\|y\|_p)^2 \Big).
\end{align*}
Then we arrive at
\begin{align*}
P_T \phi(y)
& \le \frac{\exp(T^{-1} \|G^{-1}\|_\infty^2)}{\sqrt{\mu(B(0,1))}} \exp(T^{-1} \|G^{-1}\|_\infty^2 \|y\|_p^2).
\end{align*}
This estimate, in combination with \eqref{mu-e}, yields
\begin{align*}
\|P_T\|_{2 \rightarrow 4}
&
\le \frac{\exp(T^{-1} \|G^{-1}\|_\infty^2)}{\sqrt{\mu(B(0,1))}} \sqrt[4]{\mu(\exp(4 T^{-1} \|G^{-1}\|_\infty^2 \|\cdot\|_p^2)} \\
& \le \frac{\exp(T^{-1} \|G^{-1}\|_\infty^2)}{\sqrt{\mu(B(0,1))}} \sqrt[4]{\sup_{t \ge 0} \mathbb E \exp(4 T^{-1} \|G^{-1}\|_\infty^2 \|X_t^0\|_p^2)},
\end{align*}
Since \eqref{wa2} holds for any $\epsilon_0>0$, by Lemma \ref{lm-exp}, the estimate \eqref{exp} is valid with any $\epsilon_1>0$, and thus $\|P_T\|_{2 \rightarrow 4}<\infty$ for any $T>0$, which shows the supercontractivity of $(P_t)_{t \ge 0}$.
\iffalse
To show the ultracontractivity of $(P_t)_{t \ge 0}$ when $q>2$, according to \cite[Corollary 1.3]{Wan11(AOP)}) and \cite[Theorem 2.3]{RW03(FM)}, it suffices to show that $\|P_t \exp(\epsilon \|\cdot\|_p^2)\|_\infty$ is bounded for any $t, \epsilon>0$.
This condition implies the tightness of $(\frac1n \int_0^n \nu P_t {\rm d}t)_{n \in \mathbb N_+}$, where $\nu$ is a probability measure.
Hence it implies the existence of an invariant probability measure according to a standard argument.

This is ensured by the exponential moments' estimation \eqref{exp}.
Indeed, when $p=2$, then for any $t, \epsilon>0$,
\begin{align*}
P_t \exp(\epsilon \|\cdot\|_2^2)
=\mathbb E \exp(\epsilon \|X_t^x\|_2^2)<\infty.
\end{align*}
When $p>2$, there exists a sufficiently small $\epsilon_0>0$ and a constant $C_{\epsilon_0}>0$ such that
\begin{align*}
P_t \exp(\epsilon \|\cdot\|_p^2)
=C_{\epsilon_0}\mathbb E \exp(\epsilon_0 \|X_t^x\|_p^p)<\infty.
\end{align*}

When $q>2$, Theorem \ref{tm-erg} shows the existence of the unique invariant measure $\mu$.

\begin{align*}
P_t \exp(\epsilon \|\cdot\|_2^2)
=\mathbb E \exp(\epsilon \|X_t^x\|_2^2)<\infty.
\end{align*}
When $p>2$, there exists a sufficiently small $\epsilon_0>0$ and a constant $C_{\epsilon_0}>0$ such that
\begin{align*}
P_t \exp(\epsilon \|\cdot\|_p^2)
=C_{\epsilon_0}\mathbb E \exp(\epsilon_0 \|X_t^x\|_p^p)<\infty.
\end{align*}

which is bounded for all $T>0$, by \eqref{exp}.

\begin{align*}
\|P_{2T}\|_{2 \rightarrow \infty}
& \le \|P_T P_T \phi\|_\infty
\le \sup_{y \in L^p} P_T |P_T \phi(y)| \\
& \le e^{C_1/T}\sup_{y \in L^p} P_T \exp(C_2 T^{-1} \|y\|_p^2)
\end{align*}
\fi
\end{proof}

\begin{proof}[Proof of Theorem \ref{main3}]
Theorem \ref{main3} follows from Corollary \ref{cor-hyp} and Corollary \ref{cor-sup}.
\end{proof}
\fi

\bibliographystyle{amsalpha}