\section{Introduction}
\label{intro}
The constructive description of the smoothness of functions in terms of polynomial approximations is a classical problem that started by D. Jackson and S.Bernstein results. In 1984 E. M. Dynkin gave a constructive characterization of holomorphic Besov spaces in simply connected domains in $\mathbb{C}$ with "good" boundary. We continue the research (see \cite{R13,R18_1,R18_2,Sh89}) devoted to the description of spaces of functions of several complex variables and consider spaces of holomorphic functions on strictly pseudoconvex domains, which boundary values are in Besov space. The main obstacle for polynomial approximations in strictly pseudoconvex is that polynomials are not always dense in space of holomorphic functions that are continuous up to the boundary. We restrict our consideration to strictly pseudoconvex Runge domains.
\begin{definition}
The domain $\O\subset \mathds{C}^n $ is Runge domain if for every function $f$ holomorphic in the neighbourhood of $\overline\O$ and every $\varepsilon>0$ there exist a polynomial $P$ for which $\abs{f(z)-P(z)}<\varepsilon,\ z\in\overline\O.$
\end{definition}
The condition that $\O$ is Runge is necessary and sufficient to obtain our results (see \cite{Sh89}).

\section{Main notations and definitions.}
\label{notations}

Let $ \mathds{C}^n $ be the space of $n$ complex variables, $n\geq 1,$ $z =
(z_1,\ldots, z_n),\ z_j = x_j + i y_j;$
\begin{equation*}
\partial_j f =\frac{\partial f}{\partial z_j} = \frac{1}{2}\left( \frac{\partial f}{\partial x_j} - i \frac{\partial f}{\partial y_j}\right), \quad \bar\partial_j f = \frac{\partial f}{\partial\bar{z}_j} = \frac{1}{2}\left( \frac{\partial f}{\partial x_j} + i \frac{\partial f}{\partial y_j}\right);
\end{equation*}

\[
\partial f = \sum\limits_{k=1}^{n} \frac{\partial f}{\partial z_k} dz_k,\quad \bar{\partial} f = \sum\limits_{k=1}^{n} \frac{\partial f}{\partial\bar{z}_k} d\bar{z}_k,\quad df=\partial f+  \bar{\partial} f;
\]

\[
\abs{\bar\partial f} = \abs{\bar\partial_1 f} + \ldots+\abs{\bar\partial_n f}
\]
We use the notation
$
\scp{\eta}{w} = \sum\limits_{k=1}^{n} \eta_k w_k
$
to indicate the action of the differential form $\eta=\sum\limits_{k=1}^{n}\eta_k dz_k$ of type $(1,0)$ on the vector
$w\in \mathds{C}^n .$

Let $\Omega = \left\{z\in\Cn : \rho(z)<0 \right\}$ be a strictly pseudoconvex domain with a $C^2$-smooth defining function $\rho: \mathds{C}^n \to\mathbb{R}.$ We also consider a family of domains $\O_t=\left\{z\in \mathds{C}^n  : \rho(z)<t \right\},$ where t is small real parameter, and a bijection $\Phi_t: {\partial\Omega} \to {\partial\Omega} _t$ given by the exponential map of a normal vector field. This allows us to define the projection $\Psi:\O_\varepsilon\setminus\O_{-\varepsilon}\to {\partial\Omega} $ by $\Psi(\xi)=\Phi_{\rho(\xi)}^{-1}(\xi).$

For $\xi\in {\partial\Omega} _r = \left\{\xi\in \mathds{C}^n  : \rho(\xi)=r
\right\}$ we define the tangent space
\[
 T_\xi^{\mathbb{R}} = \left\{ z\in \mathds{C}^n  : \text{\rm{Re}} \scp{\partial{\rho}(\xi)}{\xi-z} = 0  \right\},
\]
the complex tangent space
\[
T_\xi = \left\{ z\in \mathds{C}^n  : \scp{\partial{\rho}(\xi)}{\xi-z} = 0  \right\}
\]
and the complex normal vector
\begin{equation}
 n(\xi)=\abs{\bar\partial\rho(\xi)}^{-1}\left(\bar\partial_1\rho(\xi),\ldots,\bar\partial_n\rho(\xi)\right).
\end{equation}

We denote the space of holomorphic functions as $H(\O)$ and consider
the Hardy space (see~\cite{S76,FS72})
\[
H^p(\O):=\left\{f\in
H(\O):\ \norm{f}_{H^p(\O)}=\sup\limits_{-\varepsilon<t<0} \norm{f}_{L^p( {\partial\Omega} _t)} <\infty\right\},
\]
where Lebesgue spaces $L^p( {\partial\Omega} _t)$ are defined by induced on the boundary of $\O_t$ Lebesgue measure $d\sigma_t.$ We also denote $d\sigma = d\sigma_0.$ Notice that every function $f\in H^p(\O)$ has nontangential boundary limit $f^*\in L^p( {\partial\Omega} )$ and $\norm{f^*}_{H^p(\O)}\asymp\norm{f^*}_{L^p( {\partial\Omega} )}.$

Throughout this paper we use notations $\lesssim,\ \asymp.$  We write
$f\lesssim g$ if $f\leq c g$ for some constant $c>0,$ that doesn't
depend on main arguments of functions $f$ and $g$ and usually depend
only on the dimension $n$ and the domain $\O.$ Also $f\asymp g$ if $c^{-1}
g\leq f\leq c g$ for some $c>1.$ We denote the
Lebesgue measure in $ \mathds{C}^n $ as~$d\mu.$

\section{Leray-Koppelman formula and pseudoanalytic continuation}
 In the context of the theory of several complex variables there is no canonical reproducing formula, however we can use the Leray theorem that allows us to construct holomorphic reproducing kernels. The following lemma allows us to consider polynomial Leray maps for strictly pseudoconvex Runge domains.

 \begin{lemma}[Shirokov \cite{Sh89}]\label{lm:Sh}
Let $\Omega = \left\{z\in\Cn : \rho(z)<0 \right\}$ be a strictly pseudoconvex Runge domain with the $C^2$-smooth defining function $\rho.$ Then there exist functions $v(\xi,z),\  q(\xi,z),\ w_j(\xi,z),$ constants $\varepsilon,\beta,\delta,s>0$ and a domain $G\subset\mathbb{C}$ such that
\begin{enumerate}
 \item $v(\xi,z),\ q(\xi,z),\ w_j(\xi,z)$ are polynomials in $z$ and continuous in $\xi\in{\overline{\Omega}_\varepsilon\setminus\Omega}.$
 \item $v(\xi,z)= \sum\limits_{j=1}^n w_j(\xi,z)(\xi_j-z_j)=\scp{w(\xi,z)}{\xi-z}.$
 \item $w_j(\xi,z)=\partial_j\rho(\xi)+\sum\limits_{k=1}^n P_{kj}(\xi,z)(\xi_k-z_k)$ where $P_{kj}$ are polynomials in $z.$
 \item For $\xi\in{\overline{\Omega}_\varepsilon\setminus\Omega}$ and $z\in\overline{\O}$
  \begin{align*}
&  \abs{v(\xi,z)}\geq \rho(\xi)-\rho(z)+\beta\abs{\xi-z}^2,\ \abs{\xi-z}<\delta; \\
&   \abs{v(\xi,z)}\geq s>0,\ \abs{\xi-z}>\delta.
  \end{align*}
 \item For every $\xi\in{\overline{\Omega}_\varepsilon\setminus\Omega}$ and $z\in\O$ the point $\lambda(\xi,z)=v(\xi,z)q(\xi,z)$ lies in the domain $G$ and $G$ is simply connected region with $C^2$-smooth boundary, which is tangent to the $y-$axis at the origin.
\end{enumerate}
\end{lemma}
By Leray-Koppelman formula for every $f\in H^1(\O)$  we have
\begin{equation} \label{eq:CLF}
 f(z) = K f(z) =  \int\limits_{ {\partial\Omega} }  \frac{f^*(\xi) \omega(\xi,z)}{\scp{w(\xi,z)}{\xi-z}^n},\ z\in\Omega,
\end{equation}
where $\omega(\xi,z) =\frac{1}{(2\pi i)^n} w(\xi,z)\wedge \left(\bar\partial_\xi w(\xi,z)\right)^{n-1}.$

\begin{definition} The function $f\in C^1(\O_\varepsilon\setminus\overline\O)$ with $\supp{\mathbf{f}}\in\O_\varepsilon$ is a pseudoanalytic continuation of function $f\in H^1(\O)$ if nontangential boundary values of $f$ and ${\mathbf{f}}$ coinside on $ {\partial\Omega} .$
\end{definition}
Applying Stoke's theorem to Leray-Koppelman integral and pseudoanalytic continuation we have (see \cite{R18_1} for details)
\begin{equation} \label{CLF_co}
f(z) = \int\limits_{\O_\varepsilon\setminus\O} \frac{ \bar\partial {\mathbf{f}}(\xi)\wedge \omega(\xi,z)}{\scp{w(\xi,z)}{\xi-z}^n},\ z\in\O,
\end{equation}
since
\[
d_\xi\frac{ \omega(\xi,z)}{\scp{w(\xi,z)}{\xi-z}^n}  = 0,\quad z\in\O,\ \xi\in\Omega_\varepsilon\setminus\Omega.
\]
We denote the kernel by $K(\xi,z)=\scp{w(\xi,z)}{\xi-z}^{-n}.$

The function $d(\xi,z)=\abs{v(\xi,z)}=\abs{\scp{w(\xi,z)}{\xi-z}}$ defines on $ {\partial\Omega} $ a quasimetric with constant $A>0$ such that for $z,z_0,\xi\in {\partial\Omega} $
\begin{equation}\label{eq:QM}
d(z,z_0)\leq A(d(z,\xi)+d(\xi,z_0)),\ v(\xi,z)\leq Av(z,\xi).\end{equation}
If $B(z,\delta)=\{\xi\in {\partial\Omega} :d(\xi,z)<\delta\}$ is a quasiball with respect to $d$ then $\sigma(B(z,\delta))\lesssim \delta^n.$ Adapting ideas of L. Lanzani and E.M. Stein in \cite{LS13} we have the following estimate.
\begin{lemma}
$d(\xi,z)\asymp d(\Psi(\xi),z)+\rho(\xi),\ \xi\in\Omega_\varepsilon\setminus\Omega,\ z\in {\partial\Omega} .$
\end{lemma}

Thus $\mu(V(\xi,\delta))\lesssim\delta^{n+1},$ where $V(\xi,\delta)=\{z\in\Omega_\varepsilon\setminus\Omega : d(\xi,z)<\delta\}$ and analogously to \cite{LS13,R13} we have the following classical estimates.
\begin{lemma}  \label{LerayEst}
 Let $\alpha>0$ and $0<r<\delta<\varepsilon.$ Then
\begin{align*}
& \int\limits_{z\in {\partial\Omega} ,\ d(\xi,z)>\delta} \frac{d\sigma(z)}{d(\xi,z)^{n+\alpha}} \lesssim \delta^{-\alpha},\quad \xi\in\Omega_\varepsilon\setminus\Omega;\\
& \int\limits_{\xi\in {\partial\Omega} _r,\ d(\xi,z)>\delta}
\frac{d\sigma_r(\xi)}{d(\xi,z)^{n+\alpha}} \lesssim
\delta^{-\alpha},\quad z\in\Omega;\\
& \int\limits_{z\in {\partial\Omega} ,\ d(\xi,z)<\delta} \frac{d\sigma(z)}{d(\xi,z)^n} \lesssim 1+ \log{\frac{\delta}{r}},\quad \rho(\xi) = r<\varepsilon;\\
& \int\limits_{\xi\in {\partial\Omega} _r,\ d(\xi,z)<\delta} \frac{d\sigma_r(\xi)}{d(\xi,z)^n} \lesssim 1+ \log{\frac{\delta}{r}},\quad z\in\Omega.
\end{align*}

\end{lemma}

Concluding ideas by N. Shirokov we have the following polynomial approximation of a kernel.

\begin{lemma}[Shirokov \cite{Sh89}] \label{lm:KGlob}
Let $\O$ be a strongly pseudoconvex Runge domain and $\alpha>0.$ Then for every
$m\in\mathbb{N}$ there exists a function $K^{glob}_m(\xi,z)$ which is continuous in $\xi\in\Omega_\varepsilon\setminus\Omega,$  polynomial in $z$ with $\deg K^{glob}_m(\xi,\cdot)\lesssim m$ and satisfies the
following properties:

\begin{equation} \label{eq:CLF_approx1}
 \abs{K(\xi,z) - K^{glob}_m(\xi,z)} \lesssim \frac{1}{m^{\alpha}} \frac{1}{d(\xi,z)^{n+\alpha}},\quad d(\xi,z)\geq \frac{1}{m};
\end{equation}

\begin{equation} \label{eq:CLF_approx2}
 \abs{K^{glob}_m(\xi,z)}\lesssim m^n,\quad d(\xi,z)\leq\frac{1}{m}.
\end{equation}

\end{lemma}

\section{Holomorphic Besov spaces}
The main result of this work is the characterization of holomorphic functions with boundary values in Besov spaces in terms of polynomial approximations. We use definition of Besov spaces in terms of polynomial modulus of smoothness, but it coincides with classical definition in case of smooth domains (see~\cite{R13}).

\begin{definition}
Let $K\subset {\partial\Omega} ,\ 1\leq p\leq\infty.$ The best approximation of a function $f\in L^p(K)$ by polynomials of degree not greater than $m\in\mathbb{N}$ we denote by
\begin{equation} \label{df:Em}
 E_{m}(f,K)_p = \inf_{\deg P\leq m} \norm{f-P}_{L^p(K)}.
\end{equation}
\end{definition}
Notice that
\begin{equation} \label{eq:E1Ep}
 E_m(f,K)_1 \leq \sigma(K)^{1-\frac{1}{p}} E_m(f,K)_p.
\end{equation}

\begin{definition}
We say that a set $F\subset {\partial\Omega} $ is $t-$fat if there exist such $t>0$ and $\xi\in {\partial\Omega} $ that $B(\xi,t/2)\subset F \subset B(\xi,t).$ We say that a family $\{F_j\}$ of disjoint subsets is a $t-$decomposition of $ {\partial\Omega} $ if $ {\partial\Omega} =\bigcup F_j$ and every $F_j$ is $t-$fat.
\end{definition}

\begin{definition}\label{df:om_poly}
Let $f\in L^p( {\partial\Omega} ),\ 1
\leq p\leq\infty.$ We define the polynomial modulus of smoothness of degree $m\geq 0$ as
\begin{equation}
\omega_m(f,h)_p = \sup\norm{ \left\{E_{m-1}(f,F_j)_p\right\}_{j=1}^N}_{l^p},
\end{equation}
where supremum is taken over all $t-$decompositions
$\left\{F_j\right\}$ of $ {\partial\Omega} .$
\end{definition}

\begin{definition} \label{df:Bps_poly}
  Let $1\leq p,q\leq\infty,\ s>0$ and $f\in L^p( {\partial\Omega} ).$ We let $f\in B_{pq}^s( {\partial\Omega} )$ if
        $c_{pq}(f) = \norm{ \omega_m(f,t)_p t^{-s-1/q}}_{L^q(0,\varepsilon)} < \infty,$
        where $m\in\mathbb{N}$ and $m>s.$ The definition does not depend on $m>s.$
\end{definition}
We define holomorphic Besov space as
\[
A^s_{pq}(\O)=\{f\in H^p(\O): f^*\in B^s_{pq}( {\partial\Omega} )\}
\]
with norm $\norm{f}_{A^s_{pq}} = \norm{f^*}_{L^p( {\partial\Omega} )} + c_{pq}(f^*).$
\section{Local approximations on the boundary.}

In this section we study an interpolation construction of local approximation on fat subsets of $ {\partial\Omega} .$

Let $P_0$ be  polynomial in $\mathbb{R}$ of degree $m$ such that
\begin{equation*}
\int\limits_0^1 P_0(x) dx =1,\quad \int\limits_0^1 x^kP_0(x) dx =
0,\ k=1,...,m,
\end{equation*}
and $P_1$ be a polynomial in $\mathbb{C}$ of degree $m$ such that
\begin{equation*}
    \int\limits_{[-1,1]^2} P_1(z) d\mu(z) =1,\quad \int\limits_{[-1,1]^2} z^k \overline{P_1(z)} d\mu(z) = 0,\ k=1,...,m.
\end{equation*}

Without loss of generality we may assume that $\xi=0$ and $T^{\mathbb{R}}_{\xi} = \{z\in \mathds{C}^n  : z_n\in\mathbb{R}\}.$ Then projections $ \text{\rm{pr}}_\xi,\ \pi_\xi$ to the tangent and complex tangent spaces are defined by
\[
 \text{\rm{pr}}_{\xi} (z) = (z',\text{\rm{Re}} w),\ \pi_\xi(z) = (z',0),\ z=(z',w)\in\mathbb{C}^{n-1}\times\mathbb{C}.
\]
For $u\in {\partial\Omega} ,\ h,h_1>0,$ we define a rectangle $\tilde{Q}_u\subset T^{\mathbb{R}}$ by
\[
\tilde{Q}_u= \tilde{Q}_u(h,h_1) =  \text{\rm{pr}}_\xi(u) + \{z\in \mathds{C}^n  : z'\in
[0,h)^{2n-2}\in\mathbb{C}^{n-1},\ w \in[0,h_1)\}.
\]

The projection $ \text{\rm{pr}}_\xi: {\partial\Omega} \to T^{\mathbb{R}}_{\xi}$ is locally invertible and assuming that $h,h_1>0$ and $\abs{\xi-u}$ are small enough we define $Q_u =  \text{\rm{pr}}_\xi^{-1} (\tilde{Q}_u)$ and
\[
u=(u',w_0) = \text{\rm{pr}}_\xi^{-1} (z',\text{\rm{Re}}(u_n)),\ (z',w_1(z')) = \text{\rm{pr}}_\xi^{-1} (z',\text{\rm{Re}}(u_n)+h_1),
\]
where $(z',0)\in\pi_\xi(Q_u).$

Consider
a curve $\tilde{\gamma}_{z'}$ that connects points
$(z',w_0)$ and $(z',w_1(z'))$ in $ {\partial\Omega} $
\[
\tilde{\gamma}_{z'} =  \text{\rm{pr}}_{\xi}^{-1}\{ (z',\text{\rm{Re}}(u_n)+th_1),\ 0\leq t\leq 1  \}.
\]
Finally we denote by ${\gamma}_{z'}$ the curve in $\mathbb{C}$ which is generated by $n$-th coordinate of $\tilde{\gamma}_{z'}$ and connects points $w_0$ and $w_1(z').$

Let $Q'_u = \pi_\xi(Q_u)$ and define a function
\begin{equation*}
P_{Q_u}(z',w) = \frac{1}{h^{2(n-1)}}\bar{P}_{n-1}\left(\frac{z'-u'}{h}\right)
\frac{1}{w_1(z')-w_0}P_0\left(\frac{w-w_0}{w_1(z')-w_0}\right).
\end{equation*}
where $\bar{P}_{n-1}(z') =  \overline{P_1(z_1)\cdot\ldots\cdot P_1(z_{n-1})}.$

Then for every polynomial $T$ with degree not greater $m$ in every variable
\begin{multline*}
\int\limits_{Q'_u} d\mu(z') \int\limits_{\gamma_{z'}} T(z',w)
P_{J_u}(z',w) dw\\ =
 \int\limits_{Q'_u} \frac{d\mu(z')}{h^{2(n-1)}}\bar{P}_{n-1}\left(\frac{z'-u'}{h}\right) \int\limits_{w_0}^{w_1(z')}  P_0\left(\frac{w-w_0}{w_1(z')-w_0}\right)\frac{T(z',w) dw}{w_1(z')-w_0} \\
=\int\limits_{Q'_u}\bar{P}_{n-1}\left(\frac{z'-u'}{h}\right)\frac{d\mu(z')}{h^{2(n-1)}}\int\limits_{0}^1 T(z',w_0+ v (w_1(z')-w_0)) dv \\
=\int\limits_{Q'_u} T(z',w_0)\bar{P}_{n-1}\left(\frac{z'-u'}{h}\right)\frac{d\mu(z')}{h^{2(n-1)}} = T(u',w_0) = T(u).
\end{multline*}

We can replace the integral over a curve $\gamma_{z'}$ that connects points $w_0(z')$ and $w_1(z')$ with the integral over a segment $[w_0(z'),w_1(z')]$ because internal function is holomorphic.

We pass now to the description of operator that provides almost best polynomial approximation
on a $h-$fat set $K\in {\partial\Omega} $ with $B(\xi,h/2)\subset K\subset B(\xi,h).$ Note that $$ \text{\rm{pr}}_\xi(B(\xi,h))\subset [-c\sqrt{h},c\sqrt{h})^{2(n-1)}\times [-ch,ch)=\tilde{Q}\subset T^{\mathbb{R}}_\xi$$ for some $c>0.$ We cover the rectangle $\tilde{Q}$ by a union of $(m+1)^n$ rectangles $\tilde{Q}^j=\tilde{Q}_{u^j}(2c\sqrt{h}/[\sqrt{m+1}],2ch/(m+1)),$ $u_j\in {\partial\Omega} .$

Let $Q^j =  \text{\rm{pr}}_\xi^{-1}(\tilde{Q}^j).$ We define $\P_K(f)$ as a polynomial such that
\begin{equation} \label{P_K}
    \P_K(f)(u^j) = \int_{\pi_\xi(Q^j)} d\mu(z') \int_{\gamma_{z'}}
P_{Q^j}(z',w) f(z',w)dw , \ j=0,\ldots, m^n.
\end{equation}

Then $\P_K(T)(u^j)=T(u^j)$ and $\P_K(T)=T$ for every polynomial $T$ with degree not greater than $m$ in every coordinate, moreover,
\begin{equation} \label{P_est}
  \max\limits_{\dist{z}{K} \leq \lambda h} \abs{\P_K(f)(z)} \leq\frac{c(\lambda,m,\O)}{\sigma(K)}\int_K \abs{f}d\sigma,
\end{equation}
    as  $\max\limits_{B(\xi,\lambda h)} \abs{P_{Q_j}(z)} <\frac{c(\lambda,m,\Omega)}{h^n},$
    where the constant $c$ depends on $\lambda,\ m$ и $\Omega,$ and does not depend on a set $K$ and function $f.$ This implies
\[
  \norm{ f- \P_K f }_{L^p(K)}\leq (c+1) E_m(f,\ K)_p,
\]
and the polynomial $\P_K f$ gives almost best polynomial approximation of function $f$.

\begin{note} \label{rm:part_poly}
 Let $f\in H^p( {\partial\Omega} ),\ 1
 \leq p
 \leq\infty$ and $\{F_j\}$ be a $h$-decomposition of $ {\partial\Omega} .$ Define a function $T_h$ that is equal on $F_j$ to a polynomial $\P_{F_j}f.$ Then
\[
  \omega_m(f,h)_p\asymp \sup \norm{ f-T_h }_{L^p( {\partial\Omega} )},
\]
where the supremum is taken over all such decompositions.
\end{note}

\begin{note}
 The existence and uniqueness of a polynomial $P_K(f)$ needs some explanation. The uniqueness is simple induction by dimension $n.$ For $n=1$ it is obvious. Assume that we know result for dimension $(n-1).$ Let $T$ be a polynomial of degree not greater than $m$ in every variable such that $T(u^j)=0.$ Rearrange points $u^j=u^{(\alpha)}=(u_1^{\alpha_1},\ldots,u_{n-1}^{\alpha_{n-1}},u_n^{\alpha})$ by multiindex $\alpha=(\alpha_1,\ldots,\alpha_n)$ with $\alpha_k=0..m.$ Notice that points $(u_1^{\alpha_1},\ldots,u_{n-1}^{\alpha_{n-1}})$ form a rectangular net in $Q_h\subset T_\xi.$ By induction for every fixed $\alpha_1$ we have $T(u^{\alpha_1}_1,u_2,\ldots,u_n)= 0.$ Thus $T(u_1,\ldots,u_n)\equiv 0$ for every fixed $(u_2,\ldots,u_n)\in\mathbb{C}^{n-1}$ as polynomial of degree not greater than $m$ in $u_1\in\mathbb{C}$ that has $m+1$ roots.

 Proceed to a construction. Let $n_{kj}=\overline{u^k-u^j},$ then the polynomial
 $$\sum\limits_{k=0}^{m} s_k\prod\limits_{j\neq k}\frac{\scp{n_{kj}}{z-u^j}}{\scp{n_{kj}}{u^k-u^j}}$$
 has degree not greater $m$ in every variable and is equal to $s_k$ at $u^k.$
\end{note}

\section{Two methods of pseudoanalytic continuation}
Let ${\mathbf{f}}$ be a pseudoanalytic continuation of function $f\in H^1(\O)$ and $1\leq p\leq\infty.$ We introduce the following important characteristics of function ${\mathbf{f}}$
\begin{equation} \label{eq:SP}
 S_p({\mathbf{f}},r) = \norm{ \bar{\partial} {\mathbf{f}} }_{L^p( {\partial\Omega} _r)},\ r>0.
\end{equation}
In this section we generalize ideas by E.M. Dynkin \cite{D81} to construct pseudoanalytic continuations with some estimates in this value.

\subsection{Pseudoanalytic continuation by global polynomial approximations.\label{Cont_glob}}
By \cite{S76} strict pseudoconvexity of domain $\O$ implies that functions holomorphic in neighbourhood of $\O$ are dense in $H^1(\Omega).$ Also every holomorphic in neighbourhood of $\O$ function can be approximated on $\overline{\O}$ by polynomials since $\O$ is Runge. Thus there exists a polynomial sequence $P_1,P_2,\ldots$ converging to $f^*$ in $L^1( {\partial\Omega} ).$ Let
\begin{equation}
\lambda(z) = \rho(z)^{-1} \abs{P_{2^{m+1}}(z) - P_{2^{m}}(z)},\quad 2^{-m}<\rho(z)\leq 2^{-m+1}.
\end{equation}

\begin{theorem} \label{thm:ContGlob}
 Assume that $\lambda\in L^p(\Omega_\varepsilon\setminus\Omega)$ for some $p\geq 1.$ Then there exist a pseudoanalytic continuation ${\mathbf{f}}$ of the function $f$
 such that
\begin{equation}\label{eq:PAC_la}
\abs{\bar{\partial} \mathbf{f}(z)} \lesssim \lambda(z),\quad
z\in\O_\varepsilon\setminus\overline{\O}.
\end{equation}
\end{theorem}

\begin{proof}
Consider a function $\chi\in C^\infty(0,\infty)$ such that $\chi(t)=1$
for $t\leq \frac{5}{4}$ and $\chi(t)=0$ for $t\geq \frac{7}{4}.$ We let
for $m\in\mathbb{N}$
\[
\mathbf{f}_0(z) =  P_{2^{m}}(z) + \chi(2^m\rho(z)) (P_{2^{m+1}}(z) - P_{2^{m}}(z)),\ 2^{-m}<\rho(z)<2^{-m+1},
\]
and define the continuation of the function $f$ by formula ${\mathbf{f}}=\chi(2\rho(z)/\varepsilon){\mathbf{f}}_0(z).$

Now $\mathbf{f}$ is $C^1$-function on $ \mathds{C}^n \setminus\overline{\Omega}$ and
$\abs{\bar{\partial} \mathbf{f}(z)} \lesssim \lambda(z).$ We define
a function $F_m(z)$ as $F_m(z)= \mathbf{f}(z)$ for $\rho(z)>2^{-m}$
and as $F_m(z) = P_{2^{m+1}}(z)$ for $\rho(z)<2^{-m}.$ The function $F_m$ is smooth and holomorphic in $\Omega_{2^{-m}},$ and
$\abs{\bar\partial F_m(z)}\lesssim\lambda(z)$ for
$z\in \mathds{C}^n \setminus\Omega_{2^{-m}}.$ Thus
 we get
$$P_{2^{m+1}}(z) = F_m(z) = \int\limits_{\O_\varepsilon\setminus\O} \frac{\bar{\partial}F_m(\xi) \wedge\omega(\xi,z)}{\scp{w(\xi,z)}{\xi-z}^n},\ z\in\Omega,$$

We can pass to the limit in this formula by the dominated
convergence theorem; hence, the function~${\mathbf{f}}$ is a pseudoanalytic continuation of the
function~$f$. \qed
\end{proof}

\begin{lemma}\label{lm:Pm}
Let $P_{2^m}$ be a polynomial of degree $2^m$ and $1\leq q\leq\infty.$ Then
\begin{equation}
\norm{P_{2^m}}_{L^q( {\partial\Omega} _r)}\lesssim \norm{P_{2^m}}_{L^q( {\partial\Omega} )},\ 2^{-m}\leq r\leq 2^{-m+1}
\end{equation}
and the constant does not depend on $m$.
\end{lemma}
\begin{proof}
Let $\xi\in {\partial\Omega} $ and $n=n(\xi)$ be a complex normal at this point. For some $\delta>0$ and every $w\in \mathbb{CP}^{n-1}$ a set $$\gamma_w=\{u\in\mathbb{C}: (n+\xi+w)u\in {\partial\Omega} ,\ \abs{u}<\delta \}$$ is a simple non closed $C^2$-smooth curve. Also curves $\tilde{\gamma}_w=\{u(n+\xi+w):u\in\gamma_w,\ |u|<\delta\}$ cover a neighbourhood of $\xi\in {\partial\Omega} $ in a set $ {\partial\Omega} \setminus T_\xi.$

Consider a conformal map $\psi_w:\mathbb{C}\setminus\gamma_w\to\{|v|>1\}$ such that $\psi_w'(\infty)>0.$ The smoothness of $ {\partial\Omega} $ implies that the map $\psi_w$ can be continued to a conformal map of closed domains. The smoothness of $ {\partial\Omega} $ implies that there exists some neighbourhood $V=V_\xi$ of point $\xi$ in $ \mathds{C}^n $ and constant $c_1,c_2>0$ such that
\begin{equation}\label{eq}
\gamma_{w,r}=\{u\in\mathbb{C}: u(n+w)\in {\partial\Omega} _r\cap V\}\subset\{u\in\mathbb{C}: c_1r<\psi_w(u)-1<c_2r\}.
\end{equation}
Thus
\[
\abs{\psi_w(u)}^{2^{m+2}}\asymp (1+2^{-m})^{2^m}\asymp 1,\ 2^{-m}\leq \rho(u(n+w))< 2^{-m+1};
\]
\[
\abs{\psi_w(u)}=1,\ u\in\gamma_w.
\]
Consider a function
\begin{equation*}
H_w(u)=P_{2^m}(u(n+w))\psi_w(u)^{-2^{n+2}},
\end{equation*}
holomorphic in $\mathbb{C}\setminus\gamma_w$ such that $\abs{u}H_w(u) \to 0,\ u\to\infty.$ Then
\begin{equation*}
\sup\limits_{r>0}\norm{H_w}_{L^q(\abs{\psi_w}=1+r)}
\lesssim\norm{H_w}_{L^q(\gamma_w)}
\end{equation*}
and
\begin{equation*}
\sup\limits_{r>0}\norm{H_w}_{L^q(\gamma_{w,r})}
\lesssim\norm{H_w}_{L^q(\gamma_w)}.
\end{equation*}
Then
\begin{equation*}
 \norm{P_{2^m}}_{L^q(\gamma_{w,r})}\lesssim \norm{H_w}_{L^q(\gamma_{w,r})}\lesssim \norm{H_w}_{L^q(\gamma_{w})} =\norm{P_{2^m}}_{L^q(\gamma_{w})}
\end{equation*}
for $2^{-m }\leq r< 2^{-m+1}$
and integrating this estimate by $w\in \mathbb{CP}^{n-1}$ we get
\begin{equation*}
 \norm{P_{2^m}}_{L^q(V\cap {\partial\Omega} _r)}\lesssim \norm{P_{2^m}}_{L^q( {\partial\Omega} \bigcap V)},\ 2^{-m}\leq r<2^{-m+1}.
\end{equation*}

Finally, we choose the finite covering of $ {\partial\Omega} $ that also covers a set $\O_{2^{-m+1}}\setminus\O_{2^{-m}}$ and obtain the desired estimate.

\end{proof}
\begin{corollary}\label{cor:ContGlob}
The continuation ${\mathbf{f}}$ in Theorem~\ref{thm:ContGlob} satisfies an estimate
\begin{equation} \label{eq:SpEn1}
S_p({\mathbf{f}},r) \lesssim 2^m E_{2^m}(f)_p,\  2^{-m}\leq r\leq 2^{-m+1}.
\end{equation}
\end{corollary}

\subsection{Pseudoanalytic continuation by local polynomial approximations}
Let $f\in H^1(\O),\ m>0,\ z\in\Omega_\varepsilon\setminus\Omega.$ Denote $J(z)=B(\Psi(z),\rho(z)/10)$ and
\begin{equation}
E(f,z)=E_m(f,J(z))_1.
\end{equation}

\begin{theorem} \label{thm:LocCont}
 Assume that $p\geq 1$ and
 \begin{equation}
   E(f,z)\rho(z)^{-(n+1)}\in L^p(\Omega_\varepsilon\setminus\Omega).
 \end{equation}
Then $f\in L^p( {\partial\Omega} )$ and there exist such pseudoanalytic continuation ${\mathbf{f}}$ of function $f$ such that
 \begin{equation}
 \abs{\bar\partial{\mathbf{f}}(z)}\lesssim E(f,z)\rho(z)^{-(n+1)}.
 \end{equation}
\end{theorem}

\begin{proof}
First, we will prove that $f\in L^p( {\partial\Omega} ).$ Assume that $J\subset {\partial\Omega} $ is $h-fat$ and let
\[
L(J) =\left\{z\in\Omega_\varepsilon\setminus\Omega: c_1 h<d(\xi,z)<c_2 h\ \text{for some}\ \xi\in J\right\}
\]
where constants $0<c_1<c_2<1$ are chosen such that $J(z)\subset J$ for every $z\in L(J).$
Then $\mu(L(J)) \asymp h^{n+1}$ and $E_m(f,J)_1\leq E(f,z)
$ for $z\in L(J).$

Consider a series of decompositions
$  {\partial\Omega} =\bigcup\limits_{k=1}^{2^{nm}} J_k^m $
into $2^{-m}$-fat sets
assuming that every succeeding decomposition is obtained by dividing every set of previous to  $2^n$ parts.

We define a function $T_m$ on $ {\partial\Omega} $ assuming $T_m(z) = \P_{J_k^m}f(z)$
for $z\in J_k^m,$ where $\P_{J_k^m}$ is a polynomial projector defined in previous section. Then
\[
f = T_{m_0} + \sum\limits_{k=m_0}^{\infty} (T_{m+1}-T_m)
\]
and it is enough to prove that $\sum\limits_{k=m_0}^\infty\norm{T_{m+1}-T_m }_{L^p( {\partial\Omega} )}<\infty. $

For $z\in J^{m+1}_k$ we have
\begin{multline*}
\abs{T_{m+1}(z)-T_m(z)} = \abs{\P_{J_k^{m+1}}f(z) - \P_{J_k^m}f(z) }= \abs{\P_{J_k^{m+1}}\left(f(z) - \P_{J_k^m}f(z)\right) }\\ \lesssim \frac{1}{\sigma(J_k^{m+1})} \int_{J_k^{m+1}} \abs{f- \P_{J_k^m}f} d\sigma
 \lesssim 2^{nm}E(f,J_k^m)_1\\ \lesssim 2^{nm} \inf\{ E(f,z) : z\in L(J_k^m) \}.
\end{multline*}

Thus
\[
\norm{T_{m+1}-T_m}_{L^p(J_k^{m+1})}\lesssim 2^{-m(n+1)(1-1/p)} \norm{E(f,z)\rho(z)^{-(n+1)}}_{L^p(L(J_k^m))}.
\]
Consider $L_m = \{ z\in \mathds{C}^n \setminus\Omega :
c_12^{-m}<\dist{z}{\Omega}<c_22^{-m}\}.$ Then $L(J^k_m)\subset L_m$ and
\[
\norm{T_{m+1}-T_m }_{L^p( {\partial\Omega} )}\lesssim 2^{-m(n+1) \left(1-{1}/{p}\right) } \norm{E(f,z)\rho(z)^{-(n+1)}}_{L^p(L_m)}
\]
Finally
\[
 \sum\limits_{m=m_0}^{\infty} \norm{ T_{m+1}-T_m }_{L^p( {\partial\Omega} )} \leq \norm{E(f,z)\rho(z)^{-(n+1)} }_{L^p(\Omega_\varepsilon\setminus\Omega)}<\infty.
 \]

Now we will construct a continuation of function $f$. Consider Whitney decomposition in $\Omega_\varepsilon\setminus\Omega,$ $\sum\limits_{k=1}^\infty\chi_k =1,\ \chi_k\in C_0^\infty(\Omega_\varepsilon\setminus\Omega),\
\abs{\grad{\chi_k}(z)}\lesssim\rho(z)^{-1}.$ Let
$z_k\in\supp\chi_k$ and $J_k = B(\Psi(z_k),c_1\rho(z)).$  Define
\[
\mathbf{f}(z) = \sum\limits_{k=1}^{\infty} \chi_k(z)\P_{J_k}f(z), \quad z\in \mathds{C}^n \setminus\Omega.
\]
Note that for every polynomial $T(z)$ we have
\[
\mathbf{f}(z) = T(z) + \sum\limits_{k=1}^{\infty} \chi_k(z)(\P_{J_k}f(z) - T(z)), \quad z\in \mathds{C}^n \setminus\Omega
\]
and
\[
\bar{\partial} \mathbf{f}(z) = \sum\limits_{k=1}^{\infty} (\P_{J_k}f(z) - T(z)) \bar{\partial}\chi_k(z), \quad z\in \mathds{C}^n \setminus\Omega.
\]
Let $T = \P_{J(z)}f.$ Then $J_k\subset J(z)$ if
$\chi_k(z)\neq 0.$ Thus by estimate (\ref{P_est}) we get
\[
 \abs{\P_{J_k}f(z) - T(z) }= \abs{ \P_{J_k}(f-T)(z) } \lesssim \sigma(J(z))^{-1} E(f,z)\lesssim \rho(z)^{-n} E(f,z).
 \]
Also $\supp{\mathbf{f}}\subset\O_\varepsilon$ and
\[
 \abs{ \bar{\partial} \mathbf{f}(z) } \lesssim  E(f,z)\rho(z)^{-(n+1)} \in L^p(\Omega_\varepsilon\setminus\Omega).
 \]
This concludes the proof.
\end{proof}

\begin{corollary}\label{cor:LocCont}
 The continuation ${\mathbf{f}}$ of function $f$ in Theorem~\ref{thm:LocCont} satisfies an estimate $S_p({\mathbf{f}},r)\lesssim\omega_m(f,h)_p.$
\end{corollary}
\begin{proof}
Notice that by estimate (\ref{eq:E1Ep}) we have
\[
E(f,z) \leq \sigma(J(z))^{1-\frac{1}{p}}E_m(f,J(z))_p
\lesssim \rho(z)^{\left(1-\frac{1}{p}\right)n}E_m(f,J(z))_p,
\]
 where $J(z)=B(\Psi(z),\rho(z)/10).$ Then
\[
 S_p({\mathbf{f}},r) \lesssim r^{-n/p-1}\norm{ E_m(f,J(z))_p}_{L^p( {\partial\Omega} _r)}.
\]
We will estimate the integral on the right side. Quasiballs $J(z)$ cover $ {\partial\Omega} $ and we can choose a finite covering
$J_k=\{\xi\in {\partial\Omega}  : d(\xi,z_k)<r\}$ such that for every $z\in {\partial\Omega} _r$ the set $J(z)$ is a subset of $J_k$ for some $k$ and every point is covered not more than by $5^{n}$ sets. Consequently
\[
 \norm{E_m(f,J(z))}_{L^p( {\partial\Omega} _r)}\lesssim r^{n/p}\norm{\{ E_m(f,J_k)\}_{k=1}^N}_{l^p} \lesssim r^{n/p} \omega_m(f,10r)_p,
\]
and $S_p({\mathbf{f}},r)\lesssim \omega_m(f,10r)_p/r.$
\end{proof}

\subsection{Construction of local approximation from pseudoanalytic continuation.}
In the proof of Theorem~\ref{thm:ContGlob} we constructed a pseudoanalytic continuation by local approximations of a function. Assume now that we have some pseudoanalytic continuation ${\mathbf{f}}$ of function $f.$ We will construct a local polynomial approximations of $f$ on $ {\partial\Omega} $ in terms of function ${\mathbf{f}}.$

Let $J\subset B(z_0,h/A)$ for some $z_0\in {\partial\Omega} $ and $h>0,$ where $A>0$ is a constant of quasimetric $d$ defined in (\ref{eq:QM}). We consider a polynomial of degree $\leq Cm$
\begin{equation}
P_J(z) = \int\limits_{d(\xi,z_0)>2h} \bar{\partial}{\mathbf{f}}(\xi)\wedge\omega(\xi,z) K^{loc}_m(\xi,z,z_0),
\end{equation}
where the kernel $K^{loc}_m$ is a polynomial in $z$ of degree $\leq Cm$ and is defined by following equation.
\begin{multline} \label{K_loc}
K(\xi,z) = \frac{1}{v(\xi,z)^n} = \frac{1}{\left(v(\xi,z_0)+\Delta\right)^n} =
\frac{1}{v(\xi,z_0)^n}\left( 1 +
\frac{\Delta}{v(\xi,z_0)}
\right)^{-n}\\
=\frac{1}{v(\xi,z_0)^n}\Bigg( \sum\limits_{k=0}^m C_n^k\frac{\Delta^k}{v(\xi,z_0)^k}
+
O\left(\frac{\abs{\Delta}}{d(\xi,z_0)}\right)^{m+1}\Bigg)\\
= K^{loc}_m(\xi,z,z_0) +
O\left(\frac{\abs{\Delta}^{m+1}}{d(\xi,z_0)^{m+n+1}}\right),
\end{multline}
where $\Delta =v(\xi,z)-v(\xi,z_0)=\scp{w(\xi,z)}{\xi-z}-\scp{w(\xi,z_0)}{\xi-z_0}.$

Note that for $d(\xi,z_0)>2h$ and $d(z,z_0)<h/A$ we have $d(\xi,z)\asymp d(\xi,z_0)$ and
\begin{multline*}
\abs{\Delta}\leq\abs{\scp{w(\xi,z)-w(\xi,z_0)}{\xi-z}} + \abs{\scp{w(\xi,z_0)}{z-z_0}}\\
\lesssim\abs{z-z_0}\abs{\xi-z}+\abs{\scp{w(\xi,z_0)-w(z,z_0)}{z-z_0}}
+\abs{\scp{w(z,z_0)}{z-z_0}}\\
\lesssim \abs{z-z_0}\abs{\xi-z}+\abs{w(\xi,z_0)-w(z,z_0)}\abs{z-z_0}+d(z,z_0)\\
\lesssim d(z,z_0)^{1/2}d(\xi,z)^{1/2}+d(z,z_0)\lesssim d(z,z_0)^{1/2}d(\xi,z)^{1/2},
\end{multline*}
because
\begin{multline*}
\abs{w(\xi,z_0)-w(z,z_0)}\lesssim\abs{\partial\rho(\xi)-\partial\rho(z)}\\
+\sum\limits_{k,j=1}^n\left( \abs{P_{jk}(z_0,\xi)}\abs{\xi-z_0}+\abs{P_{jk}(z_0,z)}\abs{z-z_0}\right)\\
\lesssim\abs{\xi-z}+\abs{\xi-z_0}+\abs{z-z_0}\lesssim d(\xi,z)^{1/2}.
\end{multline*}
Consequently,
\begin{equation} \label{eq:ker_cont_m}
\abs{ K(\xi,z)- K^{loc}_m(\xi,z,z_0) }\lesssim
\frac{d(z,z_0)^{\frac{m+1}{2}}}{d(\xi,z)^{n+\frac{m+1}{2}}},\ d(\xi,z_0)
>2h,\ d(z,z_0)<h/A.
\end{equation}
This implies that for $z\in J$
\begin{multline} \label{eq:PolyByCont}
 \abs{ f(z) - P_J(z) }\\
  \lesssim \int\limits_{d(\xi,z_0)<2h}\abs{\bar{\partial}\mathbf{f}(\xi)} \frac{d\mu(\xi)}{d(\xi,z)^n} +
 \int\limits_{d(\xi,z_0)>2h}\abs{\bar{\partial}\mathbf{f}(\xi)} \frac{ h^{\frac{m+1}{2}} d\mu(\xi)}{d(\xi,z)^{n+\frac{m+1}{2}}}.
\end{multline}

\section{Pseudoanalytic continuation of holomorphic Besov functions\label{Cont_Aps}}
To prove two following theorems we use the Hardy inequality. Let $f$ be a positive function on $(0,\infty)$ and for $r\neq1$ define a function $F(x)$ as follows:
\begin{equation*}
    F(x) =  \left\{
  	\begin{array}{ll}
  	  \int\limits_0^{x} f(t) dt, & r>1; \\
 	   \int\limits_x^{\infty} f(t) dt, & r<1.
	\end{array}
\right.
\end{equation*}
Then for $1\leq p<\infty$ we have
\begin{equation} \label{ineq:hardy}
\int\limits_0^{\infty} x^{-r} F^p(x) dx <
\left(\frac{p}{\abs{r-1}}\right)^p \int\limits_0^{\infty} x^{-r}
(xf(x))^p dx.
\end{equation}
The proof and detailed discussion can be found in \cite{HLP52}.

\begin{theorem} \label{thm:Aps_pa}

Let $1\leq p,q\leq\infty,\ s>0$ and $f\in H^1(\Omega).$ Then $f\in
A^s_{pq}(\Omega)$ if and only if there exist a pseudoanalytic continuation ${\mathbf{f}}$ of function $f$ such that
\begin{equation} \label{eq:contcond}
\norm{S_p({\mathbf{f}},r)r^{1-1/q-s}}_{L^q[0,\varepsilon]}<\infty.
\end{equation}

\end{theorem}

\begin{proof} The sufficiency is satisfied due to Corollary~\ref{cor:LocCont} of the Theorem~\ref{thm:LocCont} and the definition of Besov spaces.
To prove necessity assume that there exist a continuation ${\mathbf{f}}$ that satisfies  condition (\ref{eq:contcond}) and that $\supp{{\mathbf{f}}}\subset\O_\varepsilon.$
Then by estimate (\ref{eq:PolyByCont}) we have
\[
	\omega_m(f,c\delta)_p \lesssim \norm{g}_{L^p( {\partial\Omega} )} + \norm{h}_{L^p( {\partial\Omega} )}
\]
for some $c>0,$ where
\begin{align*}
& g(z) = \int\limits_{d(\xi,z)<2\delta} \frac{
\abs{\bar{\partial} {\mathbf{f}}(\xi)} }{ d(\xi,z)^n } d\mu(\xi);\
h(z) = \int\limits_{d(\xi,z)>\delta}\frac{\abs{\bar{\partial} {\mathbf{f}}(\xi)}
\delta^{\frac{m+1}{2}}}{d(\xi,z)^{n+\frac{m+1}{2}}}
d\mu(\xi),\ z\in {\partial\Omega} .
\end{align*}

We introduce the following functions
\begin{align*}
 g_r(z) = \int\limits_{\substack{\xi\in {\partial\Omega} _r\\ d(\xi,z)<2\delta}} \frac{ \abs{\bar{\partial} {\mathbf{f}}(\xi)} }{ d(\xi,z)^n } d\sigma_r(\xi);\
 h_r(z) = \int\limits_{\substack{\xi\in {\partial\Omega} _r\\ d(\xi,z)>\delta}}\frac{\abs{\bar{\partial} {\mathbf{f}}(\xi)}
\delta^{\frac{m+1}{2}}}{d(\xi,z)^{n+\frac{m+1}{2}}}
d\sigma_r(\xi), z\in {\partial\Omega} .
\end{align*}
Then
\[
 \omega_m(f,c\delta)_p \lesssim \int\limits_{0}^{2\delta}\norm{g_r}_{L^p( {\partial\Omega} )} dr + \int\limits_{0}^{1}\norm{h_r}_{L^p( {\partial\Omega} )} dr.
\]
We will prove that
\begin{align*}
& \norm{g_r}_{L^p( {\partial\Omega} )} \lesssim S_p({\mathbf{f}},r) \log{\frac{\delta}{r}},\quad 0<r<2\delta;\\
& \norm{h_r}_{L^p( {\partial\Omega} )} \lesssim S_p({\mathbf{f}},r) ,\quad 0<r<\delta;\quad
\norm{h_r}_{L^p( {\partial\Omega} )} \lesssim S_p({\mathbf{f}},r)
\frac{\delta^{\frac{m+1}{2}}}{r^{\frac{m+1}{2}}} ,\quad r>\delta.
\end{align*}
By Riesz-Thorin interpolation theorem it is enough to prove these estimates for $p=1$ and $p=\infty.$ By Lemma~\ref{LerayEst} for $r>\delta$ we have
\[
\norm{h_r}_{L^1( {\partial\Omega} )} = \int\limits_{ {\partial\Omega} _r} \abs{\bar{\partial} {\mathbf{f}}(\xi)} d\sigma(\xi) \int\limits_{\substack{z\in {\partial\Omega} \\ d(\xi,z)>\delta}}\frac{\delta^{\frac{m+1}{2}}d\sigma(z)}{d(\xi,z)^{n+\frac{m+1}{2}}} \lesssim S_1({\mathbf{f}},r) \frac{\delta^{\frac{m+1}{2}}}{r^{\frac{m+1}{2}}},
\]
 for  $p=1$ and for $p=\infty$
\[
\abs{h_r(z) } \lesssim S_\infty({\mathbf{f}},r) \int\limits_{\substack{z\in {\partial\Omega} \\ d(\xi,z)>\delta}}\frac{\delta^{\frac{m+1}{2}}d\sigma(z)}{d(\xi,z)^{n+\frac{m+1}{2}}} \lesssim  S_\infty({\mathbf{f}},r)   \frac{\delta^{\frac{m+1}{2}}}{r^{\frac{m+1}{2}}},\ z\in {\partial\Omega} ,
\]
which concludes the desired estimate. Second estimate of function $h_r$ and estimate of function $g_r$
can be proven analogously.

Thus for every $\varepsilon>0,\ m\in\mathbb{N}$
\begin{equation} \label{eq:om_constr}
 \omega_m(f,c\delta)_p \lesssim \int\limits_{0}^{2\delta} S_p({\mathbf{f}},r) \frac{\delta^\varepsilon}{r^\varepsilon} dr + \delta^{\frac{m+1}{2}}\int\limits_{\delta}^{1} S_p({\mathbf{f}},r) \frac{dr}{r^{\frac{m+1}{2}}}.
\end{equation}
Assume that $1\leq q<\infty.$ We will use Hardy inequality to prove an estimate
\begin{equation*}
 \int\limits_0^\varepsilon\omega_m(f,\delta)_p^q \delta^{-qs-1}d\delta < \infty.
\end{equation*}
Consider the first term in
(\ref{eq:om_constr}) letting $\varepsilon=\frac{s}{2}$ and
$
F_1(\delta) = \int\limits_0^\delta S_p({\mathbf{f}},r) \frac{dr}{r^\varepsilon}.
$
Then applying Hardy inequality (\ref{ineq:hardy}) with $r=1+sq>1$ we have
\begin{multline*}
\int\limits_0^1 (\delta^\varepsilon F_1(\delta))^q \delta^{-1-sq} d\delta \lesssim \int\limits_0^1 r^{-1-sq+\varepsilon q} \left(S_p({\mathbf{f}},r) r^{-\varepsilon}r \right)^q dr\\ = \int\limits_0^1 \left(\frac{S_p({\mathbf{f}},r)}{r^{s-1}}\right)^q \frac{dr}{r}<\infty.
\end{multline*}
Analogously we can estimate the second term in
(\ref{eq:om_constr}) letting $s<(m+1)/2.$

For $q=\infty$ we have $S_p({\mathbf{f}},r)\lesssim r^{s-1}$ and
\[
\omega_m(f,c\delta) \lesssim \delta^\varepsilon\int_0^{2\delta}r^{s-1-\varepsilon}dr+\delta^{\frac{m+1}{2}}\int_\delta^\varepsilon r^{s-1-\frac{m+1}{2}}dr
\lesssim \delta^s.
\]
This concludes the proof of the theorem.
\end{proof}

\section{Main theorem.}

\begin{theorem} \label{thm:ConstrCond}

Let $1\leq p,q\leq\infty,\ s>0,$ and $f\in H^p( {\partial\Omega} ).$ Then $f\in
A^s_{pq}(\Omega)$ if and only if $\left\{2^{ms}E_m(f)_{p}\right\}_{m=1}^{\infty}\in l^q.$
\end{theorem}

\begin{proof}

Assume that $\left\{2^{ms}E_{2^m}(f)_{p}\right\}_{m=1}^{\infty}\in l^q.$
By Corollary~\ref{cor:ContGlob} of Theorem~\ref{thm:ContGlob} we have a pseudoanalytic continuation such that
$$S_p({\mathbf{f}},r)\lesssim 2^mE_{2^m}(f)_p. $$
Then
\begin{equation}\label{ConstrCond1}
 \norm{{S_p({\mathbf{f}},r)}{r^{-s-1/q+1}}}_{L^q(0,\varepsilon)} \lesssim \norm{ \{ {2^{ms}} E_{2^m}(f)_{p}\}_{m=1}^\infty}_{l^q}<\infty,
\end{equation}
and by Theorem~\ref{thm:Aps_pa} $f\in A^s_{pq}(\Omega).$

Conversely, let $f\in A^s_{pq}(\Omega).$ Then there exist a pseudoanalytic continuation ${\mathbf{f}}$ of function $f$ satisfying estimate (\ref{eq:contcond}). We will prove that
\begin{equation}\label{eq:Em}
 E_m(f)_p \lesssim \int\limits_0^{1/m} S_p({\mathbf{f}},r) \frac{dr}{(mr)^{\varepsilon}} dr +\int\limits_{1/m}^\infty S_p({\mathbf{f}},r) \frac{dr}{(mr)^{\alpha}},
\end{equation}
where the parameter $\alpha$ can be chosen arbitrary great and $\varepsilon>0$ is small enough. We construct polynomials approximating a function $f$ using an approximation of a kernel $K(\xi,z)$ from Lemma~\ref{lm:KGlob}. Let
\[
  P_m(z) = \int_{\CO}
\bar{\partial} {\mathbf{f}}(\xi)\wedge\omega(\xi) K^{glob}_m(\xi,z).
\]
 Then
\begin{multline*}
 \abs{f(z)-P_m(z)} \lesssim \int_{\CO} \abs{\bar{\partial}{\mathbf{f}}(\xi)} \abs{K(\xi,z) - K^{glob}_m(\xi,z) } d\mu(\xi)\\
 \lesssim  \int\limits_{d(\xi,z)<1/m} \frac{\abs{\bar{\partial} {\mathbf{f}}(\xi)}}{d(\xi,z)^n}d\mu(\xi) +  \int\limits_{\substack{d(\xi,z)<1/m,\\ \rho(\xi)<1/m}} \frac{\abs{\bar{\partial}{\mathbf{f}}(\xi)}}{m^\alpha d(\xi,z)^{n+\alpha}}d\mu(\xi)\\ + \int\limits_{\rho(\xi)>1/m}
\frac{\abs{\bar{\partial}{\mathbf{f}}(\xi)}}{m^\alpha d(\xi,z)^{n+\alpha}}d\mu(\xi)
 = U_m(z) + V_m(z) + W_m(z).
\end{multline*}

Using notations of Theorem~\ref{thm:Aps_pa} we get
\[
	U_m(z)\lesssim \int\limits_0^{1/m} dr \int_{ {\partial\Omega} _r} \frac{\abs{\bar{\partial}{\mathbf{f}}(\xi)}d\sigma_r(\xi)}{d(\xi,z)^n} =\int\limits_0^{1/m} g_r(z) dr.
\]
Then $\norm{g_r}_{L^p( {\partial\Omega} )} \lesssim
S_p({\mathbf{f}},r)\log{\frac{2}{mr}},\ 0<r<1/m,$ and by Minkovsky integral inequality
\[
	\norm{U_m}_{L^p( {\partial\Omega} )} \lesssim \int\limits_0^{1/m} S_p({\mathbf{f}},r)\log{\frac{2}{mr}}\ dr\lesssim\int\limits_0^{1/m} S_p({\mathbf{f}},r)\frac{dr}{(mr)^\varepsilon}
\]
for every $\varepsilon>0.$ Analogously
\[
	\norm{V_m}_{L^p( {\partial\Omega} )} \lesssim \int\limits_0^{1/m} S_p({\mathbf{f}},r)\ dr ;\
	\norm{W_m}_{L^p( {\partial\Omega} )} \lesssim \int\limits_{1/m}^\infty S_p({\mathbf{f}},r)
\frac{dr}{(mr)^{\alpha}},
\]
which concludes the proof of (\ref{eq:Em}).

As in the proof of Theorem~\ref{thm:Aps_pa} the final estimate for $1\leq q<\infty$ follows from Hardy inequality (\ref{ineq:hardy}).
Let $F(t) =
\int\limits_0^{t} S_p({\mathbf{f}},r)\frac{dr}{r^\varepsilon}$ and $\varepsilon<s.$ Then by monotonicity of function $F$
\begin{multline*}
 \sum\limits_{m=1}^\infty 2^{msq}\norm{U_m}_{L^p( {\partial\Omega} )}\lesssim \sum\limits_{m=1}^\infty 2^{msq} \left(\int\limits_0^{2^{-m}} S_p({\mathbf{f}},r)\frac{dr}{(2^mr)^\varepsilon}\right)^q\\
 =\sum\limits_{m=1}^\infty 2^{m(s-\varepsilon)q}F(2^{-m}) \lesssim \int\limits_{0}^{\infty} t^{(\varepsilon-s)q-1} F(t)^q dt  \\
\lesssim\int\limits_{0}^{\infty} r^{(\varepsilon-s)q-1} ( S_p({\mathbf{f}},r) r^{1-\varepsilon})^q
dr = \int\limits_{0}^{\infty} S_p({\mathbf{f}},r)^q r^{-q(s-1)-1} dr<\infty.
\end{multline*}
Analougously we get a second estimate for $V_m$ and $W_m.$

For $q=\infty$ we have $S_p({\mathbf{f}},r)\lesssim r^{s-1}$ and $E_m(f)_p\lesssim m^{-s}$ by estimate (\ref{eq:Em}).
This concludes the proof of the theorem.
\end{proof}


\section{Introduction}

The Heath-Jarrow-Morton-Musiela (HJMM) equation is a stochastic partial differential equation that models the evolution of forward rates in a market of zero coupon bonds; we refer to \cite{fillnm} for further details. It has been studied in a series of papers, see, e.g. \cite{Rusinek, Barski}, \cite{Filipovic-Tappe, Positivity} and references therein. The state space, which contains the forward curves, is a separable Hilbert space $H$ consisting of functions $h : \mathbb{R}_+ \rightarrow \mathbb{R}$. In practice, forward curves have the following features:
\begin{itemize}
\item The functions $h \in H$ become flat at the long end.

\item Consequently, the limit $\lim_{x \rightarrow \infty} h(x)$ exists.
\end{itemize}
The second property is taken into account by choosing the Hilbert space
\begin{align*}
L_{\beta}^2 \oplus \mathbb{R},
\end{align*}
where $L_{\beta}^2$ denotes the weighted Lebesgue space
\begin{align}\label{def-L-space}
L_{\beta}^2 := L^2(\mathbb{R}_+,e^{\beta x} dx)
\end{align}
for some constant $\beta > 0$.
Such spaces have been used, e.g., in \cite{Rusinek, Barski}.
As flatness of a function is measured by its derivative, the first property is taken into account by choosing the space
\begin{align}\label{def-H-space}
H_{\gamma} := \{ h : \mathbb{R}_+ \rightarrow \mathbb{R} : h \text{ is absolutely continuous with } \| h \|_{\gamma} < \infty \}
\end{align}
for some constant $\gamma > 0$, where the norm is given by
\begin{align}\label{norm-subspace}
\| h \|_{\gamma} := \bigg( |h(0)|^2 + \int_{\mathbb{R}_+} |h'(x)|^2
e^{\gamma x} dx \bigg)^{1/2}.
\end{align}
Such spaces have been introduced in \cite{fillnm} (even with more general weight functions) and further utilized, e.g., in \cite{Filipovic-Tappe, Positivity}.
Our goal of this note is to show that for all $\gamma > \beta > 0$ we have the compact embedding
\begin{align*}
H_{\gamma} \subset \subset L_{\beta}^2 \oplus \mathbb{R},
\end{align*}
that is, the forward curve spaces used in \cite{fillnm} and forthcoming papers are contained in the forward curve spaces used in \cite{Rusinek}, and the embedding is even compact. Consequently, the embedding operator between these spaces can be approximated by a sequence of finite-rank operators, and hence, when considering the HJMM equation in the state space $H_{\gamma}$, applying these operators its solutions can be approximated by a sequence of finite dimensional processes in the larger state space $L_{\beta}^2 \oplus \mathbb{R}$; we refer to Section~\ref{sec-main} for further details.

The remainder of this note is organized as follows. In Section~\ref{sec-pre} we provide the required preliminaries. In Section~\ref{sec-main} we present the embedding result and its proof, and we outline the described approximation result concerning solutions of the HJMM equation.

\section{Preliminaries and notation}\label{sec-pre}

In this section, we provide the required preliminary results and some basic notation. Concerning the upcoming results about Sobolev spaces and Fourier transforms, we refer to any textbook about functional analysis, such as \cite{Rudin} or \cite{Werner}.

As noted in the introduction, for positive real numbers $\beta,\gamma > 0$ the separable Hilbert spaces $L_{\beta}^2 \oplus \mathbb{R}$ and $H_{\gamma}$ are given by (\ref{def-L-space}) and (\ref{def-H-space}), respectively. These spaces and the forthcoming Sobolev spaces will be regarded as spaces of complex-valued functions. For every $h \in H_{\gamma}$ the limit $h(\infty) := \lim_{x \rightarrow \infty} h(x)$ exists and the subspace
\begin{align*}
H_{\gamma}^0 := \{ h \in H_{\gamma} : h(\infty) = 0 \}
\end{align*}
is a closed subspace of $H_{\gamma}$, see \cite{fillnm}.
For an open set $\Omega \subset \mathbb{R}$ we denote by $W^1(\Omega)$ the Sobolev space
\begin{align*}
W^1(\Omega) := \{ f \in L^2(\Omega) : f' \in L^2(\Omega) \text{ exists} \},
\end{align*}
which, equipped with the inner product
\begin{align}\label{norm-Sobolev}
\langle f,g \rangle_{W^1(\Omega)} = \langle f,g \rangle_{L^2(\Omega)} + \langle f',g' \rangle_{L^2(\Omega)},
\end{align}
is a separable Hilbert space. Here, derivatives are understood as weak derivatives.

For a function $h \in W^1((0,\infty))$ the extension $h \mathbbm{1}_{(0,\infty)} : \mathbb{R} \rightarrow \mathbb{C}$ does, in general, not belong to $W^1(\mathbb{R})$. In the present situation, this technical problem can be resolved as follows.
Let $h : (0,\infty) \rightarrow \mathbb{C}$ be a continuous function such that the limit $h(0) := \lim_{x \rightarrow 0} h(x)$ exists. Then we define the reflection $h^* : \mathbb{R} \rightarrow \mathbb{C}$ as
\begin{align*}
h^*(x) :=
\begin{cases}
h(x), & \text{if $x \geq 0$,}
\\ h(-x), & \text{if $x < 0$.}
\end{cases}
\end{align*}

\begin{lemma}\label{lemma-reflection}
The following statements are true:
\begin{enumerate}
\item For each $h \in W^1((0,\infty))$ we have $h^* \in W^1(\mathbb{R})$.

\item The mapping $W^1((0,\infty)) \rightarrow  W^1(\mathbb{R})$, $h \mapsto h^*$ is a bounded linear operator.

\item For each $h \in  W^1((0,\infty))$ we have
\begin{align*}
\| h \|_{W^1((0,\infty))} &\leq \| h^* \|_{W^1(\mathbb{R})} \leq \sqrt{2} \| h \|_{W^1((0,\infty))},
\\ \| h \|_{L^2((0,\infty))} &\leq \| h^* \|_{L^2(\mathbb{R})} \leq \sqrt{2} \| h \|_{L^2((0,\infty))}.
\end{align*}
\end{enumerate}
\end{lemma}

\begin{proof}
This follows from a straightforward calculation following the proof of \cite[Theorem~8.6]{Brezis}.
\end{proof}

\begin{lemma}\label{lemma-forward-Sobolev}
Let $\gamma > \beta > 0$ be arbitrary. Then the following statements are true:
\begin{enumerate}
\item We have $H_{\gamma}^0 \subset H_{\beta}^0$ and
\begin{align}\label{embedding-same}
\| h \|_{\beta} \leq \| h \|_{\gamma} \quad \text{for all $h \in H_{\gamma}^0$.}
\end{align}

\item We have $H_{\gamma}^0 \subset L_{\beta}^2$ and there is a constant $C_1 = C_1(\beta,\gamma) > 0$ such that
\begin{align}\label{embedding-diff}
\| h \|_{L_{\beta}^2} \leq C_1 \| h \|_{\gamma} \quad \text{for all $h \in H_{\gamma}^0$.}
\end{align}

\item For each $h \in H_{\gamma}^0$ we have
\begin{align*}
h e^{(\beta / 2) \bullet}|_{(0,\infty)} \in W^1((0,\infty)), \quad (h e^{(\beta / 2) \bullet}|_{(0,\infty)})^* \in W^1(\mathbb{R}),
\end{align*}
and there is a constant $C_2 = C_2(\beta,\gamma) > 0$ such that
\begin{align*}
\| (h e^{(\beta / 2) \bullet} |_{(0,\infty)})^* \|_{W^1(\mathbb{R})} &\leq C_2 \| h \|_{\gamma} \quad \text{for all $h \in H_{\gamma}^0$.}
\end{align*}
\end{enumerate}
\end{lemma}

\begin{proof}
The first statement is a direct consequence of the representation of the norm on $H_{\gamma}^0$ given by (\ref{norm-subspace}). Let $h \in H_{\gamma}^0$ be arbitrary.
By the Cauchy-Schwarz inequality we obtain
\begin{align*}
\| h \|_{L_{\beta}^2}^2 &= \int_{\mathbb{R}_+} |h(x)|^2 e^{\beta x} dx =
\int_{\mathbb{R}_+} \bigg( \int_{x}^{\infty} h'(\eta)
e^{(\gamma/2) \eta} e^{-(\gamma/2) \eta} d \eta
\bigg)^2 e^{\beta x} dx
\\ &\leq \int_{\mathbb{R}_+} \bigg( \int_{x}^{\infty} |h'(\eta)|^2 e^{\gamma \eta} d\eta
\bigg) \bigg( \int_{x}^{\infty} e^{- \gamma \eta} d\eta \bigg)
e^{\beta x} dx
\\ &\leq \int_{\mathbb{R}_+} \bigg( \int_{\mathbb{R}_+} |h'(\eta)|^2 e^{\gamma \eta} d\eta \bigg) \frac{1}{\gamma} e^{-\gamma x} e^{\beta x} dx
\\ &\leq \frac{1}{\gamma} \bigg( \int_{\mathbb{R}_+} e^{-(\gamma - \beta)x} dx \bigg) \| h \|_{\gamma}^2 = \frac{1}{\gamma(\gamma - \beta)} \| h \|_{\gamma}^2,
\end{align*}
proving the second statement. Furthermore, by (\ref{embedding-diff}) we have
\begin{align*}
\| h e^{(\beta / 2) \bullet} |_{(0,\infty)} \|_{L^2((0,\infty))}^2
&= \int_{\mathbb{R}_+} |h(x) e^{(\beta / 2)x}|^2 dx = \int_{\mathbb{R}_+} |h(x)|^2 e^{\beta x} dx
\\ &= \| h \|_{L_{\beta}^2}^2 \leq C_1^2 \| h \|_{\gamma}^2,
\end{align*}
and by estimates (\ref{embedding-same}), (\ref{embedding-diff}) we obtain
\begin{align*}
&\| (d/dx) (h e^{(\beta / 2) \bullet} |_{(0,\infty)} ) \|_{L^2((0,\infty))}^2 = \int_{\mathbb{R}_+} \Big| \frac{d}{dx} \big( h(x) e^{(\beta/2) x} \big) \Big|^2 dx
\\ &= \int_{\mathbb{R}_+} \Big| h'(x) e^{(\beta/2) x} + \frac{\beta}{2} h(x) e^{(\beta/2) x} \Big|^2 dx
\\ &\leq 2 \bigg( \int_{\mathbb{R}_+} |h'(x)|^2 e^{\beta x} dx + \frac{\beta^2}{4} \int_{\mathbb{R}_+} |h(x)|^2 e^{\beta x} dx \bigg)
\\ &\leq 2 \| h \|_{\beta}^2 + \frac{\beta^2}{2} \| h \|_{L_{\beta}^2} \leq \bigg( 2 + \frac{\beta^2 C_1^2}{2} \bigg) \| h \|_{\gamma}^2,
\end{align*}
which, together with Lemma~\ref{lemma-reflection}, concludes the proof.
\end{proof}

For $h \in L^1(\mathbb{R})$ the \emph{Fourier transform} $\mathcal{F}h : \mathbb{R} \rightarrow \mathbb{C}$ is defined as
\begin{align}\label{Fourier-transform}
(\mathcal{F} h)(\xi) := \frac{1}{\sqrt{2 \pi}} \int_{\mathbb{R}} h(x) e^{-i \xi x} dx, \quad \xi \in \mathbb{R}.
\end{align}
Recall that $C_0(\mathbb{R})$ denotes the space of all continuous functions vanishing at infinity, which, equipped with the supremum norm, is a Banach space. We have the following result:

\begin{lemma}\label{lemma-Fourier-cont}
The Fourier transform $\mathcal{F} : L^1(\mathbb{R}) \rightarrow C_0(\mathbb{R})$ is a continuous linear operator with $\| \mathcal{F} \| \leq 1 / \sqrt{2 \pi}$.
\end{lemma}

\begin{lemma}\label{lemma-functional}
Let $\gamma > \beta > 0$ be arbitrary. Then the following statements are true:
\begin{enumerate}
\item For each $h \in H_{\gamma}^0$ we have $(h e^{(\beta / 2) \bullet}|_{(0,\infty)})^* \in L^1(\mathbb{R})$ and there is a constant $C_3 = C_3(\beta,\gamma) > 0$ such that
\begin{align*}
\| (h e^{(\beta / 2) \bullet}|_{(0,\infty)})^* \|_{L^1(\mathbb{R})} \leq C_3 \| h \|_{\gamma} \quad \text{for all $h \in H_{\gamma}^0$.}
\end{align*}

\item For each $\xi \in \mathbb{R}$ the mapping
\begin{align*}
H_{\gamma}^0 \rightarrow \mathbb{R}, \quad h \mapsto \mathcal{F}(h e^{(\beta / 2) \bullet}|_{(0,\infty)})^*(\xi)
\end{align*}
is a continuous linear functional.
\end{enumerate}
\end{lemma}

\begin{proof}
We set $\delta := \frac{1}{2}(\beta + \gamma) \in (\beta,\gamma)$.
Let $h \in H_{\gamma}^0$ be arbitrary. By the Cauchy-Schwarz inequality and Lemma~\ref{lemma-forward-Sobolev} we have
\begin{align*}
&\| (h e^{(\beta / 2) \bullet}|_{(0,\infty)})^* \|_{L^1(\mathbb{R})} = 2 \| h e^{(\beta / 2) \bullet} \|_{L^1(\mathbb{R}_+)} = 2 \int_{\mathbb{R}_+} |h(x) e^{(\beta / 2)x}| dx
\\ &= 2 \int_{\mathbb{R}_+} |h(x) | e^{(\delta / 2)x} e^{-((\delta - \beta)/2) x} dx
\\ &\leq 2 \bigg( \int_{\mathbb{R}_+} |h(x)|^2 e^{\delta x} dx \bigg)^{1/2} \bigg( \int_{\mathbb{R}_+} e^{-(\delta - \beta)x} dx \bigg)^{1/2}
\\ &= 2 \sqrt{\frac{1}{\delta - \beta}} \| h \|_{L_{\delta}^2} \leq 2 C_1(\delta,\gamma) \sqrt{\frac{1}{\delta - \beta}} \| h \|_{\gamma},
\end{align*}
showing the first statement. Moreover, we have
\begin{align*}
\| e^{((\beta / 2) - \delta) \bullet} \|_{L_{\delta}^2}^2 = \int_{\mathbb{R}_+} e^{2((\beta / 2) - \delta)x} e^{\delta x} dx = \int_{\mathbb{R}_+} e^{-(\delta - \beta)x} dx = \frac{1}{\delta - \beta},
\end{align*}
showing that $e^{((\beta / 2) - \delta) \bullet} \in L_{\delta}^2$. Let $h \in H_{\gamma}^0$ and $\xi \in \mathbb{R}$ be arbitrary. By Lemma~\ref{lemma-forward-Sobolev} we have $h \in L_{\delta}^2$, and hence
\begin{align*}
&\mathcal{F}(h e^{(\beta / 2) \bullet}|_{(0,\infty)})^*(\xi)
\\ &= \frac{1}{\sqrt{2 \pi}} \bigg( \int_0^{\infty} h(x) e^{(\beta/2) x} e^{- i \xi x} dx + \int_{-\infty}^0 h(-x) e^{-(\beta / 2) x} e^{- i \xi x} dx \bigg)
\\ &= \frac{1}{\sqrt{2 \pi}} \bigg( \int_0^{\infty} h(x) e^{(\beta/2) x} e^{- i \xi x} dx + \int_0^{\infty} h(x) e^{(\beta / 2) x} e^{i \xi x} dx \bigg)
\\ &= \frac{1}{\sqrt{2 \pi}} \big\langle h, e^{((\beta / 2) - \delta) \bullet} \big( e^{-i \xi \bullet} + e^{i \xi \bullet} \big) \big\rangle_{L_{\delta}^2},
\end{align*}
proving the second statement.
\end{proof}

We can also define the Fourier transform on $L^2(\mathbb{R})$ such that
$\mathcal{F} : L^2(\mathbb{R}) \rightarrow L^2(\mathbb{R})$ is a bijection and we have the Plancherel isometry
\begin{align}\label{Plancherel}
\langle \mathcal{F}f,\mathcal{F}g \rangle_{L^2(\mathbb{R})} = \langle f,g \rangle_{L^2(\mathbb{R})} \quad \text{for all $f,g \in L^2(\mathbb{R})$.}
\end{align}
Moreover, the two just reviewed definitions of the Fourier transform coincide on $L^1(\mathbb{R}) \cap L^2(\mathbb{R})$. For each $h \in W^1(\mathbb{R})$ we have
\begin{align}\label{Fourier-diff}
(\mathcal{F} h')(\xi) = i \xi (\mathcal{F} h)(\xi), \quad \xi \in \mathbb{R}.
\end{align}

\begin{lemma}\label{lemma-est-Fourier-diff}
For every $h \in W^1(\mathbb{R})$ we have
\begin{align*}
\| \bullet \mathcal{F} h \|_{L^2(\mathbb{R})} \leq \| h \|_{W^1(\mathbb{R})}.
\end{align*}
\end{lemma}

\begin{proof}
Let $h \in W^1(\mathbb{R})$ be arbitrary. By identity (\ref{Fourier-diff}) and the Plancherel isometry (\ref{Plancherel}) we have
\begin{align*}
\| \bullet \mathcal{F} h \|_{L^2(\mathbb{R})} = \| \mathcal{F} h' \|_{L^2(\mathbb{R})} = \| h' \|_{L^2(\mathbb{R})} \leq \| h \|_{W^1(\mathbb{R})},
\end{align*}
finishing the proof.
\end{proof}

\section{The embedding result and its proof}\label{sec-main}

In this section, we present the compact embedding result and its proof.

\begin{theorem}\label{thm-comp-embedding}
For all $\gamma > \beta > 0$ we have the compact embedding
\begin{align*}
H_{\gamma} \subset \subset L_{\beta}^2 \oplus \mathbb{R}.
\end{align*}
\end{theorem}

\begin{proof}
Noting that $H_{\gamma} \cong H_{\gamma}^0 \oplus \mathbb{R}$, it suffices to prove the compact embedding $H_{\gamma}^0 \subset \subset L_{\beta}^2$.
Let $(h_j)_{j \in \mathbb{N}} \subset H_{\gamma}^0$ be a bounded sequence. Then there exists a subsequence which converges weakly in $H_{\gamma}^0$. Without loss of generality, we may assume that the original sequence $(h_j)_{j \in \mathbb{N}}$ converges weakly in $H_{\gamma}^0$. We shall prove that $(h_j)_{j \in \mathbb{N}}$ is a Cauchy sequence in $L_{\beta}^2$. According to Lemma~\ref{lemma-forward-Sobolev}, the sequence $(g_j)_{j \in \mathbb{N}}$ given by
\begin{align*}
g_j := ( h_j e^{(\beta / 2) \bullet}|_{(0,\infty)} )^*, \quad j \in \mathbb{N}
\end{align*}
is a bounded sequence in $W^1(\mathbb{R})$.
By Lemma~\ref{lemma-reflection} and the Plancherel isometry (\ref{Plancherel}), for all $j,k \in \mathbb{N}$ we get
\begin{align*}
&\| h_k - h_j \|_{L_{\beta}^2}^2 = \| h_k e^{(\beta / 2) \bullet} - h_j e^{(\beta / 2) \bullet} \|_{L^2(\mathbb{R}_+)}^2 \leq \| g_k - g_j \|_{L^2(\mathbb{R})}^2
\\ &= \| \mathcal{F} g_k - \mathcal{F} g_j \|_{L^2(\mathbb{R})}^2
= \int_{\mathbb{R}} | (\mathcal{F} g_k)(x) - (\mathcal{F} g_j)(x) |^2 dx.
\end{align*}
Thus, for every $R > 0$ we obtain the estimate
\begin{equation}\label{comp-conv}
\begin{aligned}
\| h_k - h_j \|_{L_{\beta}^2}^2 &\leq \int_{\{ |x| \leq R \}} | (\mathcal{F} g_k)(x) - \mathcal{F} (g_j)(x) |^2 dx
\\ &\quad + \int_{\{ |x| > R \}} | (\mathcal{F} g_k)(x) - \mathcal{F} (g_j)(x) |^2 dx.
\end{aligned}
\end{equation}
By Lemma~\ref{lemma-est-Fourier-diff},  the sequence $( \bullet \mathcal{F} g_j )_{j \in \mathbb{N}}$ is bounded in $L^2(\mathbb{R})$. Therefore, for an arbitrary $\epsilon > 0$ there exists a real number $R > 0$ such that
\begin{equation}\label{comp-conv-1}
\begin{aligned}
&\int_{\{ |x| > R \}} | ( \mathcal{F} g_k)(x) - (\mathcal{F} g_j)(x) |^2 dx
\\ &\leq \frac{1}{R^2} \int_{\{ |x| > R \}} |x|^2 | (\mathcal{F} g_k)(x) - (\mathcal{F} g_j)(x) |^2 dx < \epsilon \quad \text{for all $j,k \in \mathbb{N}$.}
\end{aligned}
\end{equation}
By Lemma~\ref{lemma-functional}, for each $\xi \in \mathbb{R}$ the mapping
\begin{align*}
H_{\gamma}^0 \rightarrow \mathbb{R}, \quad h \mapsto \mathcal{F}(h e^{(\beta / 2) \bullet}|_{(0,\infty)})^*(\xi)
\end{align*}
is a continuous linear functional. Consequently, since $(h_j)_{j \in \mathbb{N}}$ converges weakly in $H_{\gamma}^0$, for each $\xi \in \mathbb{R}$ the real-valued sequence $( (\mathcal{F} g_j)(\xi) )_{j \in \mathbb{N}}$
is convergent. Moreover, by Lemmas~\ref{lemma-Fourier-cont} and \ref{lemma-functional}, for all $h \in H_{\gamma}^0$ we have the estimate
\begin{align*}
\| \mathcal{F}((h e^{(\beta / 2) \bullet}|_{(0,\infty)})^*) \|_{C_0(\mathbb{R})} \leq \frac{1}{\sqrt{2 \pi}} \| (h e^{(\beta / 2) \bullet}|_{(0,\infty)})^* \|_{L^1(\mathbb{R})}
\leq \frac{C_3}{\sqrt{2 \pi}} \| h \|_{\gamma}.
\end{align*}
Therefore, the sequence $( \mathcal{F} g_j )_{j \in \mathbb{N}}$ is bounded in $C_0(\mathbb{R})$. Using Lebesgue's dominated convergence theorem, we deduce that
\begin{align}\label{comp-conv-2}
\int_{\{ |x| \leq R \}} | (\mathcal{F} g_k)(x) - (\mathcal{F} g_j)(x) |^2 dx \rightarrow 0 \quad \text{for $j,k \rightarrow \infty$.}
\end{align}
Combining (\ref{comp-conv}) together with (\ref{comp-conv-1}) and (\ref{comp-conv-2}) shows that $(h_j)_{j \in \mathbb{N}}$ is a Cauchy sequence in $L_{\beta}^2$, completing the proof.
\end{proof}

\begin{remark}
Note that the proof of Theorem~\ref{thm-comp-embedding} has certain analogies to the proof of the classical Rellich embedding theorem (see, e.g., \cite[Satz~V.2.13]{Werner}), which states the compact embedding $H_0^1(\Omega) \subset \subset L^2(\Omega)$ for an open, bounded subset $\Omega \subset \mathbb{R}^n$. Here $H_0^1(\Omega)$ denotes the Sobolev space $H_0^1(\Omega) = \overline{\mathcal{D}(\Omega)}$, where $\mathcal{D}(\Omega)$ is the space of all $C^{\infty}$-functions on $\Omega$ with compact support, and where the closure is taken with respect to the topology induced by the inner product $\langle \cdot,\cdot \rangle_{W^1}$. Let us briefly describe the analogies and differences between the two results:
\begin{itemize}
\item In the classical Rellich embedding theorem the domain $\Omega$ is assumed to be bounded, whereas in Theorem~\ref{thm-comp-embedding} we have $\Omega = \mathbb{R}_+$. Moreover, we consider weighted function spaces with weight functions of the type $w(x) = e^{\beta x}$ for some constant $\beta > 0$. This requires a careful analysis of the results regarding Fourier transforms which we have adapted to the present situation; see Lemma~\ref{lemma-functional}.

\item $H_{\gamma}$ and $H_0^1(\Omega)$ are different kinds of spaces. While the norm on $H_0^1(\Omega)$ given by (\ref{norm-Sobolev}) involves the $L^2$-norms of a function $h$ and its derivative $h'$, the norm (\ref{norm-subspace}) on $H_{\gamma}$ only involves the $L^2$-norm of the derivative $h'$ and a point evaluation. Therefore, the embedding $H_0^1(\Omega) \subset L^2(\Omega)$ follows right away, whereas we require the assumption $\beta < \gamma$ for the embedding $H_{\gamma}^0 \subset L_{\beta}^2$; see Lemma~\ref{lemma-forward-Sobolev}.

\item The classical Rellich embedding theorem does not need to be true with $H_0^1(\Omega)$ being replaced by $W^1(\Omega)$. The reason behind it is that, in general, it is not possible to extend a function $h \in W^1(\Omega)$ to a function $\tilde{h} \in W^1(\mathbb{R}^n)$, which, however, is crucial in order to apply the results about Fourier transforms. Usually, one assumes that $\Omega$ satisfies a so-called cone condition, see, e.g., \cite{Adams} for further details. In our situation, we have to ensure that every function $h \in H_{\gamma}^0$ can be extended to a function $\tilde{h} \in W^1(\mathbb{R})$, and this is provided by Lemma~\ref{lemma-forward-Sobolev}.
\end{itemize}
\end{remark}

For the rest of this section, we shall describe the announced application regarding the approximation of solutions to semilinear stochastic partial differential equations (SPDEs), which in particular applies to the modeling of interest rates. Consider a SPDE of the form
\begin{align}\label{SPDE}
\left\{
\begin{array}{rcl}
dr_t & = & (A r_t + \alpha(t,r_t)) dt + \sigma(t,r_t)dW_t
\\ && + \int_E \gamma(t,r_{t-},\xi) (\mathfrak{p}(dt,d\xi) - \nu(d\xi)dt) \medskip
\\ r_0 & = & h_0
\end{array}
\right.
\end{align}
on some separable Hilbert space $H_1$ with $A$ denoting the generator of some strongly continuous semigroup on $H_1$, driven by a Wiener process $W$ and a homogeneous Poisson random measure $\mathfrak{p}$ with compensator $dt \otimes \nu(d\xi)$ on some mark space $E$.
We assume that the standard Lipschitz and linear growth conditions are satisfied which ensure for each initial condition $h_0 \in H_1$ the existence of a unique weak solution $r$ to (\ref{SPDE}), that is, for each $\zeta \in \mathcal{D}(A^*)$ we have almost surely
\begin{align*}
\langle \zeta,r_t \rangle &= \langle \zeta,h_0 \rangle_{H_1} + \int_0^t \big( \langle A^* \zeta,r_s \rangle_{H_1} + \langle \zeta,\alpha(s,r_s) \rangle_{H_1} \big) ds + \int_0^t \langle \zeta,\sigma(s,r_s) \rangle_{H_1} dW_s
\\ &\quad + \int_0^t \int_E \langle \zeta,\gamma(s,r_{s-},\xi) \rangle_{H_1} (\mathfrak{p}(ds,d\xi) - \nu(d\xi)ds) \quad \text{for all $t \geq 0$,}
\end{align*}
see, e.g., \cite{SPDE} for further details. Let $H_2$ be a larger separable Hilbert space with compact embedding $H_1 \subset \subset H_2$. By virtue of Theorem~\ref{thm-comp-embedding}, this is in particular satisfied for the forward curve spaces $H_1 = H_{\gamma}$ and $H_2 = L_{\beta}^2 \oplus \mathbb{R}$ for $\gamma > \beta > 0$. If, furthermore, $A = d/dx$ is the differential operator, which is generated by the translation semigroup $(S_t)_{t \geq 0}$ given by $S_t h = h(t + \bullet)$, and $\alpha = \alpha_{\rm HJM}$ is given by the so-called HJM drift condition
\begin{align*}
\alpha_{\rm HJM}(t,h) &= \sum_j \sigma^j(t,h) \int_0^{\bullet} \sigma^j(t,h)(\eta) d\eta
\\ &\quad - \int_E \gamma(t,h,\xi) \bigg[ \exp \bigg( - \int_0^{\bullet} \gamma(t,h,\xi)(\eta)d\eta \bigg) - 1 \bigg] \nu(d\xi),
\end{align*}
then the SPDE (\ref{SPDE}), which in this case becomes the mentioned HJMM equation, describes the evolution of interest rates in an arbitrage free bond market; we refer to \cite{Positivity} for further details.

By virtue of the compact embedding $H_1 \subset \subset H_2$, there exist orthonormal systems $(e_k)_{k \in \mathbb{N}}$ of $H_1$ and $(f_k)_{k \in \mathbb{N}}$ of $H_2$, and a decreasing sequence $(s_k)_{k \in \mathbb{N}} \subset \mathbb{R}_+$ with $s_k \rightarrow 0$ such that
\begin{align*}
h = \sum_{k=1}^{\infty} s_k \langle h,e_k \rangle_{H_1} f_k \quad \text{for all $h \in H_1$,}
\end{align*}
see, e.g., \cite[Satz~VI.3.6]{Werner}. The numbers $s_k$ are the singular numbers of the identity operator ${\rm Id} : H_1 \rightarrow H_2$. Defining the sequence $(T_n)_{n \in \mathbb{N}}$ of finite-rank operators
\begin{align*}
T_n : H_1 \rightarrow F_n, \quad T_n h := \sum_{k=1}^n s_k \langle h,e_k \rangle_{H_1} f_k,
\end{align*}
where $F_n :=  \langle f_1,\ldots,f_n \rangle$,
we even have $T_n \rightarrow {\rm Id}$ with respect to the operator norm
\begin{align*}
\| T \| := \sup_{\| h \|_{H_1} \leq 1} \| Th \|_{H_2},
\end{align*}
see, e.g., \cite[Korollar~VI.3.7]{Werner}. Consequently, denoting by $r$ the weak solution to the SPDE (\ref{SPDE}) for some initial condition $h_0 \in H_1$, the sequence $(T_n(r))_{n \in \mathbb{N}}$ is a sequence of $F_n$-valued stochastic processes, and we have almost surely
\begin{align}\label{norm-conv}
\| T_n(r_t) - r_t \|_{H_2} \leq \| T_n - {\rm Id} \| \, \| r_t \|_{H_1} \rightarrow 0 \quad \text{for all $t \geq 0$,}
\end{align}
showing that the weak solution $r$ -- when considered on the larger state space $H_2$ -- can be approximated by the sequence of finite dimensional processes $(T_n(r))_{n \in \mathbb{N}}$ with distance between $T_n(r)$ and $r$ estimated in terms of the operator norm $\| T_n - {\rm Id} \|$, as shown in (\ref{norm-conv}).
However, the sequence $(T_n(r))_{n \in \mathbb{N}}$ does not need to be a sequence of It\^{o} processes. This issue is addressed by the following result:

\begin{proposition}\label{prop-approx}
Let $(\epsilon_n)_{n \in \mathbb{N}} \subset (0,\infty)$ be an arbitrary decreasing sequence with $\epsilon_n \rightarrow 0$. Then, for every initial condition $h_0 \in H_1$ there exists a sequence $(r^{(n)})_{n \in \mathbb{N}}$ of $F_n$-valued It\^{o} processes such that almost surely
\begin{align}\label{est-epsilon-n}
\| r_t^{(n)} - r_t \|_{H_2} \leq ( \| T_n - {\rm Id} \| + \epsilon_n ) \| r_t \|_{H_1} \rightarrow 0 \quad \text{for all $t \geq 0$,}
\end{align}
where $r$ denotes the weak solution to (\ref{SPDE}).
\end{proposition}

\begin{proof}
According to \cite[Theorems~13.35.c and 13.12]{Rudin}, the domain $\mathcal{D}(A^*)$ is dense in $H_1$. Therefore, for each $n \in \mathbb{N}$ there exist elements $\zeta_1^{(n)},\ldots,\zeta_n^{(n)} \in \mathcal{D}(A^*)$ such that
\begin{align*}
\| \zeta_k^{(n)} - e_k \|_{H_1} < \frac{\epsilon_n}{2^k \cdot s_k} \quad \text{for all $k=1,\ldots,n$,}
\end{align*}
where we use the convention $\frac{x}{0} := \infty$ for $x > 0$.
We define the sequence $(S_n)_{n \in \mathbb{N}}$ of finite-rank operators as
\begin{align*}
S_n : H_1 \rightarrow F_n, \quad S_n h := \sum_{k=1}^n s_k \langle h,\zeta_k^{(n)} \rangle_{H_1} f_k.
\end{align*}
By the geometric series, for all $n \in \mathbb{N}$ we have
\begin{align*}
\| S_n - {\rm Id} \| &\leq \| S_n - T_n \| + \| T_n - {\rm Id} \| \leq \sum_{k=1}^n s_k \| \langle \bullet,\zeta_k^{(n)} - e_k \rangle_{H_1} \| + \| T_n - {\rm Id} \|
\\ &\leq \epsilon_n \sum_{k=1}^n \frac{1}{2^k} + \| T_n - {\rm Id} \| \leq \epsilon_n + \| T_n - {\rm Id} \|.
\end{align*}
For each $n \in \mathbb{N}$ let $r^{(n)}$ be the $F_n$-valued It\^{o} process
\begin{align*}
r_t^{(n)} = h_0^{(n)} + \int_0^t \alpha_s^{(n)} ds + \int_0^t \sigma_s^{(n)} dW_s + \int_0^t \int_E \delta_s^{(n)}(\xi) (\mathfrak{p}(ds,d\xi) - \nu(d\xi,ds)),
\end{align*}
with parameters given by
\begin{align*}
h_0^{(n)} &= \sum_{k=1}^n s_k \langle \zeta_k^{(n)},h_0 \rangle_{H_1} f_k, \quad \alpha_t^{(n)} = \sum_{k=1}^n s_k \big( \langle A^* \zeta_k^{(n)},r_t \rangle_{H_1} + \langle \zeta_k^{(n)},\alpha(t,r_t) \rangle_{H_1} \big) f_k,
\\ \sigma_t^{(n)} &= \sum_{k=1}^n s_k \langle \zeta_k^{(n)},\sigma(t,r_t) \rangle_{H_1} f_k, \quad \delta_t^{(n)}(\xi) = \sum_{k=1}^n s_k \langle \zeta_k^{(n)},\delta(t,r_{t-},\xi) \rangle_{H_1} f_k.
\end{align*}
Since $r$ is a weak solution to (\ref{SPDE}), we obtain almost surely
\begin{align*}
S_n(r_t) &= \sum_{k=1}^n s_k \langle \zeta_k^{(n)},r_t \rangle_{H_1} f_k
\\ &= \sum_{k=1}^n s_k \bigg( \langle \zeta_k^{(n)},h_0 \rangle_{H_1} + \int_0^t \big( \langle A^* \zeta_k^{(n)},r_s \rangle_{H_1} + \langle \zeta_k^{(n)},\alpha(s,r_s) \rangle_{H_1} \big) ds
\\ &\qquad\qquad\,\,\, + \int_0^t \langle \zeta_k^{(n)},\sigma(s,r_s) \rangle_{H_1} dW_s
\\ &\qquad\qquad\,\,\, + \int_0^t \int_E \langle \zeta_k^{(n)},\delta(s,r_{s-},\xi) \rangle_{H_1} (\mathfrak{p}(ds,d\xi) - \nu(d\xi,ds)) \bigg) f_k
\\ &= h_0^{(n)} + \int_0^t \alpha_s^{(n)} ds + \int_0^t \sigma_s^{(n)} dW_s + \int_0^t \int_E \delta_s^{(n)}(\xi) (\mathfrak{p}(ds,d\xi) - \nu(d\xi,ds))
\\ &= r_t^{(n)} \quad \text{for all $t \geq 0$,}
\end{align*}
which finishes the proof.
\end{proof}

We shall conclude this section with further consequences regarding the speed of convergence of the approximations $(r^{(n)})_{n \in \mathbb{N}}$ provided by Proposition~\ref{prop-approx}. Let $h_0 \in H_1$ be an arbitrary initial condition and denote by $r$ the weak solution to (\ref{SPDE}). Furthermore, let $T > 0$ be a finite time horizon. Since
\begin{align*}
\mathbb{E} \bigg[ \sup_{t \in [0,T]} \| r_t \|_{H_1}^2 \bigg] < \infty,
\end{align*}
see, e.g., \cite[Corollary~10.3]{SPDE}, by (\ref{est-epsilon-n}) there exists a constant $K > 0$ such that
\begin{align*}
\mathbb{E} \bigg[ \sup_{t \in [0,T]} \| r_t^{(n)} - r_t \|_{H_2}^2 \bigg]^{1/2} \leq K (\| T_n - {\rm Id} \| + \epsilon_n) \rightarrow 0,
\end{align*}
providing a uniform estimate for the distance of $r^{(n)}$ and $r$ in the mean-square sense. Moreover, considering the pure diffusion case
\begin{align*}
\left\{
\begin{array}{rcl}
dr_t & = & (A r_t + \alpha(t,r_t)) dt + \sigma(t,r_t)dW_t \medskip
\\ r_0 & = & h_0,
\end{array}
\right.
\end{align*}
the sample paths of $r$ are continuous, for every constant $K > \| h_0 \|_{H_1}$ the stopping time
\begin{align*}
\tau := \inf \{ t \geq 0 : \| r_t \| \geq K \}
\end{align*}
is strictly positive, and by (\ref{est-epsilon-n}) for the stopped processes we obtain almost surely
\begin{align}\label{uni-local}
\sup_{t \in \mathbb{R}_+} \| r_{t \wedge \tau}^{(n)} - r_{t \wedge \tau} \|_{H_2} \leq K (\| T_n - {\rm Id} \| + \epsilon_n) \rightarrow 0,
\end{align}
i.e., locally the solution $r$ stays in a bounded subset of $H_{\gamma}$ and we obtain the uniform convergence (\ref{uni-local}).

\section*{Acknowledgement}

The author is grateful to an anonymous referee for valuable comments and suggestions.