Skip to main content

Boundary value problems for modified Dirac operators in Clifford analysis

Abstract

In this paper, we discuss two kinds of Riemann type boundary value problems for the operator \(\widetilde{D}_{\lambda}\), where λ is a complex number. Furthermore, we establish the Almansi type expansion for the operator \(\widetilde{D}_{\lambda}^{k}\), where \(k\in\mathbf{N}\). As applications of the expansion, we investigate the Riemann type boundary value problem and the generalized Riquier problem for the operator \(\widetilde{D}_{\lambda}^{k}\).

1 Introduction

The uniqueness and existence theorems for the solutions of boundary value problems for systems of partial differential equations are sufficiently well known. Such problems have remarkable applications in mathematical physics, the mechanics of deformable bodies, electromagnetism, relativistic quantum mechanics, and some of their natural generalizations. Almost all such problems can be set in the context of Clifford analysis (see [1, 2]). Clifford analysis is centered around the concept of monogenic functions, i.e. null solutions of a first order vector valued rotation invariant differential operator called the Dirac operator which factorizes the Laplace operator (see [3, 4]). As to the mathematical study of boundary value problems in Clifford analysis, there are several different approaches known in the literature. Without claiming completeness, we mention some of them. First of all, we have the approach originating with Bernstein, whose approach is to translate boundary value problems to the corresponding singular integral equations, then use the properties of the Fredholm operator to discuss the solvability of singular integral equations (see [5]). Another important approach is based on complex analysis. In this case, first we use analytic function theory to solve these kinds of boundary value problems, then we use the results of boundary value problems to solve singular integral equations (see [6, 7]). The advantage of this method is that the explicit representation of solutions can be obtained, but in the higher dimensional space there still exist many obstacles to generalize this method. In this paper, we continue to use the method in [6, 7] to solve boundary value problems for the modified Dirac operators.

The paper is organized as follows. In Section 2, we review some results on the theory of Clifford analysis. In Section 3, applying the Plemelj formula for the modified Dirac operator [6], we consider Riemann type boundary value problems for the operator \(\widetilde{D}_{\lambda}\). In Section 4, using the Euler operator in Clifford analysis, we obtain the Almansi type expansion for the operator \(\widetilde{D}_{\lambda}^{k}\). In Section 5, as applications of the expansion, we investigate the Riemann type boundary value problem and the generalized Riquier problem for the operator \(\widetilde{D}_{\lambda}^{k}\).

2 Preliminaries

2.1 Clifford analysis

Let \(\mathbf{R}_{0,m}\) be the real associative Clifford algebra generated by \(\{e_{1}, e_{2}, \ldots, e_{m}\}\), where the basic vectors \(e_{1}, e_{2}, \ldots, e_{m}\) satisfy the relations \(e_{i}e_{j}+e_{j}e_{i}=-2\delta_{i,j}\), \(i,j=1,\ldots,m\). Let \(\varepsilon_{i}=-e_{1}e_{i}\), \(i=1,\ldots,m\), then the universal Clifford algebra \(\mathbf{R}_{0, m-1}\) for \(\mathbf{R}^{m-1}\) is generated by \(\{\varepsilon_{1}, \varepsilon_{2}, \ldots, \varepsilon_{m}\}\), where the vectors \(\varepsilon_{1}, \varepsilon_{2}, \ldots, \varepsilon_{m}\) satisfy the following relations:

$$\left \{ \textstyle\begin{array}{@{}l} \varepsilon_{1}\varepsilon_{i}=\varepsilon_{i}\varepsilon_{1},\quad i=1,\ldots,m,\\ \varepsilon_{i}\varepsilon_{j}+\varepsilon_{j}\varepsilon_{i}=-2\delta_{i,j}, \quad i,j=2,\ldots,m. \end{array}\displaystyle \right . $$

Each of the elements in \(\mathbf{R}_{0, m-1}\) may be written as \(a=\sum_{A} a_{A}\varepsilon_{A}\), where \(a_{A}\) are real numbers and \(\varepsilon_{A}=\varepsilon_{\alpha _{1}}\varepsilon_{\alpha_{2}}\cdots\varepsilon_{\alpha_{h}}\) with \(A=\{\alpha_{1},\ldots,\alpha_{h}\}\subset\{2,\ldots,m\}\). We define the norm of a as \(|a|= (\sum_{A} |a_{A}|^{2} )^{\frac{1}{2}}\). If there exists \(b\in\mathbf{R}_{0,m-1}\) such that \(ab=ba=\varepsilon _{1}\), then b is called the inverse of a, which is denoted as \(a^{-1}\).

A typical element of \(\mathbf{R}^{m}\) is denoted by \(x=x_{1}\varepsilon _{1}+x_{2}\varepsilon_{2}+\cdots+x_{m}\varepsilon_{m}\) with \(x_{i}\in \mathbf{R}\). We define \(\overline{x}=x_{1}\varepsilon_{1}-x_{2}\varepsilon_{2}-\cdots -x_{m}\varepsilon_{m}\), then \(x\overline{x}=\overline{x}x=|x|^{2}\). Obviously, for \(x\neq0\), we have \(x^{-1}=\frac{\overline{x}}{|x|^{2}}\).

One of the main aims of Clifford analysis is to construct a first order operator, the so-called Dirac operator, factorizing the Laplace operator and to study the function-theoretical properties of the null solutions of this operator. When working over \(\mathbf{R}^{m}\), this Dirac operator is defined by

$$ D=\sum_{i=1}^{m}{e}_{i} \partial_{{ x}_{i}}. $$
(1)

Then the modified Dirac operator is defined as

$$ \widetilde{D}=\sum_{i=1}^{m} \varepsilon_{i}\partial_{x_{i}}. $$
(2)

When studying the modified Dirac operator in this setting, we consider functions f which are e.g. elements of spaces such as \(C^{k}(\Omega)\otimes\mathbf{R}_{0,m-1}\) with Ω some open domain in \(\mathbf{R}^{m}\). This means that f can be written as

$$ f=\sum_{A} f_{A}(x){ \varepsilon}_{A} $$
(3)

with \(f_{A}(x)\in C^{k}(\Omega)\). Denote by \(|f|= (\sum_{A} |f_{A}(x)| )^{\frac{1}{2}}\) the norm of \(f\in C^{k}(\Omega)\otimes\mathbf{R}_{0,m-1}\).

3 Boundary value problems for the operator \(\widetilde{D}_{\lambda}\)

3.1 Riemann type problem for the operator \(\widetilde{D}_{\lambda}\)

Let

$$ E(x)=\frac{1}{\omega_{m}}\frac{\overline{x}}{|x|}, $$
(4)

where \(\omega_{m}=\frac{2\pi^{\frac{m}{2}}}{\Gamma(\frac{m}{2})}\) is the surface area of the unit sphere in \(\mathbf{R}^{m}\). Then \(E(x)\) satisfies the equation \(\widetilde{D}f=0\).

Let f be a Hölder continuous function on Ω and take its Cauchy transform

$$ f(x)=\frac{1}{A_{m}}\int_{\partial\Omega}E(y-x)\,d \sigma_{y}f(y),\quad x\in\mathbf{R}^{m}\setminus{\partial \Omega}. $$
(5)

Then \(f(x)\) satisfies the equation \(\widetilde{D}f=0\) in \(\mathbf {R}^{m}\setminus{\partial\Omega}\) as was proved in [6].

In [6], the following Plemelj formulas hold for \(s\in\partial\Omega\):

$$ F^{+}(s)=\lim_{{x\in\Omega, x\rightarrow s}}f(x)= \frac{1}{2}f(s)+\frac{1}{A_{m}}\int_{\partial\Omega}E(y-s)\,d \sigma_{y}f(y) $$
(6)

and

$$ F^{-}(s)=\lim_{{x\in\mathbf{R}^{m}\setminus{\overline{\Omega}}, x\rightarrow s}}f(x)=- \frac{1}{2}f(s)+ \frac{1}{A_{m}}\int_{\partial \Omega}E(y-s)\,d \sigma_{y}f(y), $$
(7)

where \(\overline{\Omega}=\Omega\cup{\partial\Omega}\).

In order to obtain the main result in this section, we need the following lemma.

Lemma 3.1

Let \(x_{1}\) be a nonzero finite real number and \(\lambda\in{{C}}\). Then

$$ \operatorname{ker}\widetilde{D}_{\lambda} =e^{\lambda x_{1}} \operatorname{ker}\widetilde{D}, $$
(8)

where \(\operatorname{ker}\widetilde{D}_{\lambda}=\{f|f\in C^{1}(\Omega)\otimes R_{0,m-1},(\widetilde{D}- \lambda)f=0\}\), and \(\operatorname{ker}\widetilde{D}=\operatorname{ker}\widetilde{D}_{\lambda}\) for \(\lambda=0\).

Proof

Letting \(f\in \operatorname{ker}\widetilde{D}\), we have

$$\begin{aligned} \widetilde{D}_{\lambda} \bigl(e^{\lambda x_{1}}f \bigr) =&(\widetilde{D}- \lambda)e^{\lambda x_{1}}f \\ =&\lambda e^{\lambda x_{1}}f+e^{\lambda x_{1}}\widetilde{D}f-\lambda e^{\lambda x_{1}}f \\ =&0, \end{aligned}$$

which implies that \(e^{\lambda x_{1}}\operatorname{ker}\widetilde{D}\subset \operatorname{ker}\widetilde{D}_{\lambda}\).

On the contrary, for \(f\in\widetilde{D}_{\lambda}\), we can see that

$$\begin{aligned} \widetilde{D} \bigl(e^{-\lambda x_{1}}f \bigr) = e^{-\lambda x_{1}}(\widetilde{D}- \lambda)f =0, \end{aligned}$$

which means that \(\operatorname{ker}\widetilde{D}_{\lambda}\subset e^{\lambda x_{1}}\operatorname{ker}\widetilde{D} \). □

Therefore, we obtain the conclusion.

Theorem 3.2

Let f be a Hölder continuous function on Ω and let \(G\in Z(\mathbf{R}_{0,m-1})\) be invertible with inverse \(G^{-1}\). Then the Riemann type problem

$$ \left \{ \textstyle\begin{array}{@{}l} \widetilde{D}_{\lambda}\Phi(x)=0, \quad x\in\mathbf{R}^{m}\setminus {\partial\Omega},\\ \Phi^{+}(t)=G \Phi^{-}(t)+f(t), \quad t\in\partial\Omega,\\ \Phi^{-}(\infty)=0, \end{array}\displaystyle \right . $$
(I)

has a solution Φ given by

$$ \Phi(x)=\frac{X(x)}{A_{m}}\int_{\partial\Omega }e^{\lambda x_{1}}E(y-x) \,d\sigma_{y}G^{-1}f(y), $$
(9)

where

$$ X(x)=\left \{ \textstyle\begin{array}{@{}l@{\quad}l} G, & x\in\Omega,\\ 1, & x\in\mathbf{R}^{m}\setminus{\overline{\Omega}}. \end{array}\displaystyle \right . $$
(10)

Note that the center \(Z(\mathbf{R}_{0,m-1})\) of \(\mathbf{R}_{0,m-1}\) is the set of elements in \(\mathbf{R}_{0,m-1}\) which commute with all elements of \(\mathbf{R}_{0,m-1}\) (see e.g. [6])

Proof

First, it follows by Lemma 3.1 that the function \(\Phi(x)\) determined by (9) satisfies the equation

$$\widetilde{D}_{\lambda}\Phi(x)=0. $$

Secondly, let \(G\in Z(\mathbf{R}_{0,m-1})\) be invertible with inverse \(G^{-1}\). It follows by (6) and (7) that

$$\begin{aligned} & \Phi^{+}(s)-G\Phi^{-}(s) \\ &\quad= \biggl(\frac{X^{+}(s)}{2}G^{-1}f(s)+\frac{1}{A_{m}} \int _{\partial\Omega}e^{\lambda x_{1}}E(y-s)\,d\sigma_{y}G^{-1}f(y) \biggr) \\ &\qquad{} -G \biggl(-\frac{X^{-}(s)}{2}G^{-1}f(s)+\frac{1}{A_{m}} \int _{\partial\Omega}e^{\lambda x_{1}}E(y-s)\,d\sigma_{y}G^{-1}f(y) \biggr) \\ &\quad= f(s), \end{aligned}$$

where \(s\in\partial\Omega\).

Finally, it is obvious that the function \(\Phi(x)\) vanishes at infinity.

Thus, we obtain the conclusion. □

3.2 Riemann type boundary value problem (II)

In this section, using the Plemelj formulas, we consider the following Riemann type boundary value problem (II).

Suppose that f is a Hölder continuous function on Ω. Find a function \(\Psi\in C^{1}(\Omega)\otimes\mathbf{R}_{0,m-1}\) that satisfies

$$ \left \{ \textstyle\begin{array}{@{}l} \widetilde{D}_{\lambda}\Psi(x)=0, \quad x\in\mathbf{R}^{m}\setminus {\partial\Omega},\\ \Psi^{+}(t)a-\Psi^{-}(t)b=f(t), \quad t\in\partial\Omega,\\ \lim_{|x|\rightarrow+\infty}\frac{|\Psi (x)|}{|x|^{l}}< +\infty, \end{array}\displaystyle \right . $$
(II)

where

$$\Psi^{+}(t)=\lim_{x\rightarrow t}\Psi(x),\quad x\in\Omega; \qquad \Psi^{-}(t)=\lim_{x\rightarrow t}\Psi(x),\quad x\in{ \mathbf {R}^{m}\backslash\overline{\Omega}}. $$

a, b are given \(\mathbf{R}_{0,m-1}\) valued constants whose inverses are \(a^{-1}\), \(b^{-1}\).

Lemma 3.3

Let \(Z_{k}=x_{k}\varepsilon_{1}-x_{1}\varepsilon_{k}\), where \(2\leq k\leq m\). Then we have the polynomials of order p

$$V_{k_{1}\cdots{k_{p}}}=\frac{1}{p!}\sum_{\pi({k_{1}\cdots k_{p}})}Z_{k_{1}} \cdots Z_{k_{p}}\in \operatorname{ker}\widetilde{D}, $$

where the sum runs over all distinguishable permutations of all of \((k_{1},\ldots,k_{p})\).

The proof of Lemma 3.3 is similar to Proposition 11.2.3 in [3].

Theorem 3.4

The boundary value problem (II) has a solution.

Proof

We will prove the function

$$\Psi(x)=\left \{ \textstyle\begin{array}{@{}l@{\quad}l} (\int_{\partial\Omega}e^{\lambda x_{1}}E(y-x)\,d\sigma _{y}f(y) )a^{-1}\\ \quad{}+\sum_{p=1}^{l}\sum_{\pi({k_{1}\cdots k_{p}})}V_{k_{1}\cdots{k_{p}}}e^{\lambda x_{1}}a^{-1}, & x\in {\Omega}, l\in\mathbf{N},\\ (\int_{\partial\Omega}e^{\lambda x_{1}}E(y-x)\,d\sigma _{y}f(y) )b^{-1}\\ \quad{}+\sum_{p=1}^{l}\sum_{\pi({k_{1}\cdots k_{p}})}V_{k_{1}\cdots{k_{p}}}e^{\lambda x_{1}}b^{-1}, & x\in R^{m}\backslash\overline{\Omega}, l\in\mathbf{N}, \end{array}\displaystyle \right . $$

is a solution of the boundary value problem (II).

Denote

$$\Phi(x)=\left \{ \textstyle\begin{array}{@{}l@{\quad}l} \Psi(x)a, & x\in\Omega,\\ \Psi(x)b, & x\in{R^{m}\backslash\overline{\Omega}}. \end{array}\displaystyle \right . $$

Then

$$\left \{ \textstyle\begin{array}{@{}l@{\quad}l} \Phi^{+}(t)=\Psi^{+}(t)a, & t\in\partial\Omega,\\ \Phi^{-}(t)=\Psi^{-}(t)b, & t\in\partial\Omega, \end{array}\displaystyle \right . $$

where a, b have the inverses \(a^{-1}\), \(b^{-1}\), respectively. The boundary value problem (II) is equivalent to

$$\Phi^{+}(t)-\Phi^{-}(t)=f(t), \quad t\in\partial\Omega. $$

Note that

$$\bigl(T[f] \bigr) (x)=\int_{\partial\Omega}e^{\lambda x_{1}}E(y-x)\,d \sigma_{y}f(y) $$

for \(x\in\mathbf{R}^{m}\backslash\partial\Omega\) is meaningful and satisfies the boundary properties

$$\bigl(T[f] \bigr)^{+}(t)- \bigl(T[f] \bigr)^{-}(t)=f(t), \quad t\in \partial\Omega. $$

Thus

$$\Phi^{+}(t)- \bigl(T[f] \bigr)^{+}(t)=\Phi^{-}(t)- \bigl(T[f] \bigr)^{-}(t),\quad t\in\partial \Omega, $$

which means that \(\Phi(x)-(T[f])(x)=g(x)\in \operatorname{ker}\widetilde{D}_{\lambda}\) in \(\mathbf{R}^{m}\) by the Painlevé theorem. By Lemma 3.3, we put \(g(x)=\sum_{p=1}^{l}\sum_{\pi ({k_{1}\cdots k_{p}})}V_{k_{1}\cdots{k_{p}}}e^{\lambda x_{1}}\). Thus we have the conclusion. □

4 Almansi type expansion for the operator \(\widetilde{D}_{\lambda}^{k}\)

In 1899, the Almansi expansion for polyharmonic functions was established, which was equivalent to the Fischer decomposition for polynomials (see [8]). One can find important applications and generalizations of this result in the case of several complex variables in the monograph of Aronszajn et al. [9], e.g. concerning functions holomorphic in the neighborhood of the origin in \(C^{n}\). Also for the case of a Clifford analysis, one can refer to [10, 11]. But all these cases are limited to star-like domains. In this section, we consider the difficult case that Ω is some open domain in \(\mathbf{R}^{m}\) not limited to star-like domains.

Definition 4.1

We define the generalized Euler operator by

$$\mathbf{E}_{s}= s\mathbf{I} + \mathbf{E} =s\mathbf{I} + {\sum _{i=1}^{m}{{ x}_{i}} \partial_{{ x}_{i}}}, $$

where s is a complex constant, I is the identity operator, and E is the Euler operator.

Lemma 4.2

Let Ω be as stated before. For \(f( x)\in C^{2}(\Omega)\otimes\mathbf{R}_{0,m-1}\),

$$ \widetilde{D}\mathbf{E}_{s}f( x)=\mathbf{E}_{s+1} \widetilde{D}f( x), $$
(11)

where \(s\in{C}\).

Proof

For \(s=0\), from Definition 4.1 it follows that, for \(f( x)\in C^{2}(\Omega)\otimes\mathbf{R}_{0,m-1}\),

$$\begin{aligned} {\widetilde{D}\mathbf{E}f( x)} =& {\sum_{i=1}^{m}{ \varepsilon}_{i}{\partial _{{x}_{i}}} \Biggl(\sum _{j=1}^{m}{{ x}_{j}}\partial_{{x}_{j}}f(x) \Biggr)} = {\sum_{i=1}^{m}{\varepsilon}_{i}{ \partial _{{x}_{i}}}f(x)+\sum_{i,j=1}^{m}{ \varepsilon}_{i}{\partial_{{x}_{j}}} \frac{\partial^{2}f}{\partial x_{j}\,\partial x_{i}}} = {\widetilde{D}f(x)+\mathbf{E}\widetilde{D}f(x)}. \end{aligned}$$

This implies that \(\widetilde{D}\mathbf{E}=\mathbf{E_{1}}\widetilde{D}\). For \(s\neq0\),

$$\widetilde{D}\mathbf{E}_{s}=\widetilde{D}(s+\mathbf{E})=s\widetilde {D}+\mathbf{E_{1}}\widetilde{D}=\mathbf{E_{s+1}} \widetilde{D}. $$

This completes the lemma. □

Note that the proof of Lemma 4.2 is inspired by Ren in [11].

Lemma 4.3

If \(f\in \operatorname{ker}(\widetilde{D}_{\lambda})\), then

$$ C_{k}\widetilde{D}_{\lambda}^{k}\mathbf{E}_{\lambda}^{k}f=f, $$
(12)

where \(C_{k}=\frac{1}{k!\lambda^{k}}\) and \(k\in\mathbf{N}\).

Proof

Note that \(f\in \operatorname{ker}\widetilde{D}_{\lambda}\). For \(k=1\), Lemma 4.2 implies that

$$\begin{aligned} \widetilde{D}_{\lambda}\mathbf{E}_{\lambda}f =&\widetilde{D}_{\lambda} \mathbf{E}_{\lambda}f \\ =&\widetilde{D}\mathbf{E}_{\lambda}f-\lambda\mathbf{E}_{\lambda}f \\ =&\mathbf{E}_{\lambda+1}\widetilde{D}f-\lambda\mathbf{E}_{\lambda}f \\ =&\mathbf{E}_{\lambda+1}\widetilde{D}f-\lambda\mathbf{E}_{\lambda +1}f+ \lambda f \\ =&\mathbf{E}_{\lambda+1}(\widetilde{D}-\lambda)f+\lambda f \\ =&\lambda f. \end{aligned}$$

Suppose that, for \(k=l\),

$$C_{l}\widetilde{D}_{\lambda}^{l}\mathbf{E}_{\lambda}^{l}f=f, $$

where \(C_{l}=\frac{1}{l!\lambda^{l}}\). For \(k=l+1\),

$$\widetilde{D}_{\lambda}^{l+1}\mathbf{E}_{\lambda}^{l}f =\widetilde {D}_{\lambda}\widetilde{D}_{\lambda}^{l}\mathbf{\mathbf{E}}_{\lambda }^{l}f=\frac{1}{C_{l}}\widetilde{D}_{\lambda}f=0. $$

We calculate

$$\begin{aligned} \widetilde{D}_{\lambda}^{l+1}\mathbf{E}_{\lambda}^{l+1}f =&\widetilde {D}_{\lambda}^{l}\widetilde{D}_{\lambda} \mathbf{E}_{\lambda}\mathbf {E}_{\lambda}^{l}f \\ =&\widetilde{D}_{\lambda}^{l}(\mathbf{E}_{\lambda+1} \widetilde {D}_{\lambda}+\lambda)\mathbf{E}_{\lambda}^{l}f \\ =&\widetilde{D}_{\lambda}^{l}\mathbf{E}_{\lambda+1} \widetilde {D}_{\lambda}\mathbf{E}_{\lambda}^{l}f+\lambda D_{\lambda}^{l}\mathbf {E}_{\lambda}^{l}f \\ =&\widetilde{D}_{\lambda}^{l-1}\widetilde{D}_{\lambda} \mathbf {E}_{\lambda+1}\widetilde{D}_{\lambda}\mathbf{E}_{\lambda}^{l}f+ \frac{\lambda}{C_{l}}f \\ =&\widetilde{D}_{\lambda}^{l-1}\mathbf{E}_{\lambda+2} \widetilde {D}_{\lambda}^{2}\mathbf{E}_{\lambda}^{l}f+ \frac{2\lambda}{C_{l}}f \\ =&\cdots \\ =&\mathbf{E}_{\lambda+l+1}\widetilde{D}_{\lambda}^{l+1}\mathbf {E}_{\lambda}^{l}f+\frac{(l+1)\lambda}{C_{l}}f \\ =& \frac{1}{C_{l+1}} f, \end{aligned}$$

which implies the conclusion. □

Denote \(\operatorname{ker}\widetilde{D}_{\lambda}^{k}=\{f|f\in C^{k}(\Omega)\otimes R_{0,m-1}, (\widetilde{D}-\lambda)^{k}f=0, k\in\mathbf{N}\}\).

Theorem 4.4

If \(f(x)\in \operatorname{ker}\widetilde{D}_{\lambda}^{k}\), then there exist unique functions \(f_{0},\ldots, f_{k-1}\in \operatorname{ker}\widetilde{D}_{\lambda}\) such that

$$ f(x)=f_{0}(x)+\mathbf{E}_{\lambda}f_{1}(x)+ \mathbf{E}_{\lambda }^{2}f_{2}(x)+\cdots+ \mathbf{E}_{\lambda}^{k-1}f_{k-1}(x), $$
(13)

where \(f_{0},\ldots, f_{k-1}\) are given as follows:

$$ \left \{ \textstyle\begin{array}{@{}l} {f_{0}(x)= (\mathbf{I}-C_{1}\mathbf{E}_{\lambda}\widetilde {D}_{\lambda} ) (\mathbf{I} -C_{2}\mathbf{E}_{\lambda}^{2}\widetilde{D}_{\lambda}^{2} )\cdots (\mathbf{I}-C_{k-1}\mathbf{E}_{\lambda}^{k-1}\widetilde{D}_{\lambda }^{k-1} )f( x)},\\ {f_{1}(x)=C_{1}\widetilde{D}_{\lambda} (\mathbf{I} -C_{2}\mathbf{E}_{\lambda}^{2}\widetilde{D}_{\lambda}^{2} )\cdots (\mathbf{I}-C_{k-1}\mathbf{E}_{\lambda}^{k-1}\widetilde{D}_{\lambda }^{k-1} )f( x) },\\ \vdots\\ {f_{k-2}(x)=C_{k-2}\widetilde{D}_{\lambda}^{k-2}(\mathbf {I}-C_{k-1}\mathbf{E}_{\lambda}^{k-1}\widetilde{D}_{\lambda }^{k-1})f(x)},\\ { f_{k-1}( x)=C_{k-1}\widetilde{D}_{\lambda}^{k-1}f(x)}, \end{array}\displaystyle \right . $$
(14)

and \(C_{k}=\frac{1}{k!\lambda^{k}}\).

Conversely, if functions \(f_{0},\ldots, f_{k-1}\in \operatorname{ker}\widetilde {D}_{\lambda}\), then the function \(f(x)\) given by (13) satisfies the equation \(\widetilde{D}_{\lambda}^{k}f=0\).

Proof

If we let the operator \(\widetilde{D}_{\lambda}^{k-1}\) act on (13), then by Lemma 4.3, we have

$$\begin{aligned} & {\widetilde{D}_{\lambda}^{k-1}f(x)} \\ &\quad= {\widetilde{D}_{\lambda}^{k-1} \Biggl(f_{0}(x) +\sum_{i=1}^{k-1}\mathbf{E}_{\lambda}^{i}f_{i}(x) \Biggr)} \\ &\quad= {\widetilde{D}_{\lambda}^{k-1}\mathbf{E}_{\lambda }^{k-1}f_{k-1}(x)} \\ &\quad= {\frac{1}{C_{k-1}}f_{k-1}(x)}. \end{aligned}$$

Thus,

$$f_{k-1}( x)=C_{k-1}\widetilde{D}_{\lambda}^{k-1}f(x). $$

Similarly, if we let the operator \(\widetilde{D}_{\lambda}^{k-2}\) act on \(f(x)-\mathbf{E}_{\lambda}^{k-1}f_{k-1}(x)\), we have

$$\begin{aligned} & \widetilde{D}_{\lambda}^{k-2}\bigl[f(x)-\mathbf{E}_{\lambda }^{k-1}f_{k-1}(x) \bigr] \\ &\quad= \widetilde{D}_{\lambda}^{k-2} \Biggl(f_{0}(x)+ \sum_{i=1}^{k-2}\mathbf{E}_{\lambda}^{i}f_{i}(x) \Biggr) \\ &\quad= \widetilde{D}_{\lambda}^{k-2} \bigl(\mathbf{E}_{\lambda }^{k-2}f_{k-2}(x) \bigr) \\ &\quad= \frac{1}{C_{k-2}}f_{k-2}(x). \end{aligned}$$

Therefore, we have

$$f_{k-2}(x)=C_{k-2}\widetilde{D}_{\lambda}^{k-2} \bigl(\mathbf {I}-C_{k-1}\mathbf{E}_{\lambda}^{k-1} \widetilde{D}_{\lambda}^{k-1} \bigr)f(x). $$

By induction, we have

$$ \left \{ \textstyle\begin{array}{@{}l} f_{k-1}(x)=C_{k-1}\widetilde{D}_{\lambda}^{k-1}f(x),\\ {f_{k-2}(x)=C_{k-2}\widetilde{D}_{\lambda}^{k-2}(\mathbf {I}-C_{k-1}\mathbf{E}_{\lambda}^{k-1}\widetilde{D}_{\lambda }^{k-1})f(x)},\\ \vdots\\ {f_{1}(x)=C_{1}\widetilde{D}_{\lambda} (\mathbf{I} -C_{2}\mathbf{E}_{\lambda}^{2}\widetilde{D}_{\lambda}^{2} )\cdots (\mathbf{I}-C_{k-1}\mathbf{E}_{\lambda}^{k-1}\widetilde{D}_{\lambda }^{k-1} )f( x) },\\ {f_{0}(x)= (\mathbf{I}-C_{1}\mathbf{E}_{\lambda}\widetilde {D}_{\lambda} ) (\mathbf{I} -C_{2}\mathbf{E}_{\lambda}^{2}\widetilde{D}_{\lambda}^{2} )\cdots (\mathbf{I}-C_{k-1}\mathbf{E}_{\lambda}^{k-1}\widetilde{D}_{\lambda }^{k-1} )f( x)}. \end{array}\displaystyle \right . $$
(15)

Conversely, suppose that the functions \(f_{0},\ldots, f_{k-1}\in \operatorname{ker}\widetilde{D}_{\lambda}\). Applying Lemma 4.3, we obtain

$$\widetilde{D}_{\lambda}^{k}f(x)=\widetilde{D}_{\lambda}^{k} \Biggl[f_{0}(x)+\sum_{i=1}^{k-1} \mathbf{E}_{\lambda}^{i}f_{i}(x) \Biggr]=0, $$

which completes the proof. □

5 Boundary value problems for the operator \(\widetilde{D}_{\lambda}^{k}\)

5.1 Riemann type boundary value problem (III)

Now we consider the following Riemann type boundary value problem (III).

Suppose that \(g_{l}(t)\), \(l=0,\ldots,k-1\), are Hölder continuous functions on Ω. Find a function \(\Psi\in C^{k}(\Omega)\otimes\mathbf{R}_{0,m-1}\) that satisfies

$$ \left \{ \textstyle\begin{array}{@{}l} \widetilde{D}_{\lambda}^{k}\Psi(x)=0, \quad x\in\mathbf{R}^{m}\setminus{\partial\Omega},\\ \Psi^{+}(t)a-\Psi^{-}(t)b=g_{0}(t),\\ {[}\widetilde{D}_{\lambda}\Psi]^{+}(t)a-[\widetilde{D}_{\lambda}\Psi ]^{-}(t)b=g_{1}(t),\\ {\vdots}\\ {[}\widetilde{D}_{\lambda}^{k-1}\Psi]^{+}(t)a-[\widetilde{D}_{\lambda }^{k-1}\Psi]^{-}(t)b=g_{k-1}(t), \quad t\in\partial\Omega,\\ \lim_{|x|\rightarrow+\infty}\frac{|\Psi (x)|}{|x|^{l}}< +\infty, \end{array}\displaystyle \right . $$
(III)

where

$$\Psi^{+}(t)=\lim_{x\rightarrow t}x, \quad x\in{\Omega},\qquad \Psi^{-}(t)=\lim_{x\rightarrow t} x, \quad x\in \mathbf{R}^{m}\backslash\overline{\Omega}, $$

and a, b are given \(\mathbf{R}_{0,m}\) valued constants whose inverses are \(a^{-1}\), \(b^{-1}\).

Theorem 5.1

The boundary value problem (III) has a solution.

Proof

We will prove that the function

$$ \Phi( x)=F_{0}(x)+\mathbf{E}_{\lambda}F_{1}( x)+\cdots +\mathbf{E}_{\lambda}^{k-1}F_{k-1}(x), $$
(16)

where

$$ F_{i}(x)=\left \{ \textstyle\begin{array}{@{}l@{\quad}l} (\int_{\partial\Omega}e^{\lambda x_{1}}E(y-x)\,d\sigma _{y}g_{i}(y) )[a]^{-1}\\ \quad{}+\sum_{p=1}^{l}\sum_{\pi({k_{1}\cdots k_{p}})}V_{k_{1}\cdots{k_{p}}}e^{\lambda x_{1}}[a]^{-1}, \quad x\in {\Omega}, l\in\mathbf{N},\\ { (\int_{\partial\Omega}e^{\lambda x_{1}}E(y-x)\,d\sigma _{y}g_{i}(y) )[b]^{-1}}\\ \quad{}+\sum_{p=1}^{l}\sum_{\pi({k_{1}\cdots k_{p}})}V_{k_{1}\cdots{k_{p}}}e^{\lambda x_{1}}[b]^{-1}, \quad x\in {\mathbf{R}^{m}\setminus\overline{\Omega}}, l\in\mathbf{N}, \end{array}\displaystyle \right . $$
(17)

for \(0\leq i\leq k-1\), and

$$ F_{i}(t)= C_{i} \Biggl[g_{i}(t)- \widetilde{D}_{\lambda}^{i}{\sum_{l=i+1}^{k-i-1}}E_{\lambda}^{l}F_{l}(t) \Biggr],\quad t\in\partial \Omega, $$
(18)

is a solution of the boundary value problem (III).

From Theorem 3.4, we can see that \(F_{i}(x)\in\operatorname{ker}\widetilde {D}_{\lambda}\). It follows by Theorem 4.4 that \(\Phi(x)\in \operatorname{ker}\widetilde{D}_{\lambda}^{k}\).

Then, applying Lemma 4.3 and (18), we can see that

$$\begin{aligned} &{\bigl[\widetilde{D}_{\lambda}^{i}F\bigr]^{+}(t)a- \bigl[\widetilde {D}_{\lambda}^{i}F\bigr]^{-}(t)b} = {\frac{1}{C_{i}}F_{i}(t)+\widetilde{D}_{\lambda }^{i}{ \sum_{l=i+1}^{k-i-1}}E_{\lambda}^{l}F_{l}(t)} = {g_{i}(t)}, \end{aligned}$$

which completes the proof. □

5.2 Generalized Riquier problem for the operator \(\widetilde{D}_{\lambda}^{k}\)

In 1936, Nicolescu established Riquier problem for polyharmonic equations (see [12]). In 2003, applying the 0-normalized system of functions with respect to the Laplace operator, Karachik obtained a solution of the Riquier problem in harmonic analysis (see [13]). In this section, we will study the generalized Riquier problem for the operator \(\widetilde{D}_{\lambda}^{k}\) by the expansion (13), as follows:

Find a function Φ such that \(\widetilde{D}_{\lambda}^{i}\Phi\in {C(\overline{\Omega})\otimes\mathbf{R}_{0,m-1}}\), for \(i=0,\ldots ,k-1\), and

$$ \left \{ \textstyle\begin{array}{@{}l} \widetilde{D}_{\lambda}^{k}\Phi=0, \quad x\in\Omega,\\ \widetilde{D}_{\lambda}^{i}\Phi|_{\partial\Omega}=g_{i}(t),\quad t\in\partial\Omega. \end{array}\displaystyle \right . $$
(IV)

Theorem 5.2

Suppose that the functions \(f_{i}(x)\in{C(\overline{\Omega})\otimes \mathbf{R}_{0,m-1}}\), \(i=0,\ldots,k-1\). Then problem (IV) has a solution given by

$$ \Phi(x)=f_{0}(x)+{\sum_{i=1}^{k-1}} \mathbf{E}_{\lambda}^{i}f_{i}(x), $$
(19)

where the functions \(f_{i}(x)\) satisfy

$$ \left \{ \textstyle\begin{array}{@{}l} \widetilde{D}_{\lambda}f_{i}(x)=0, \quad x\in\Omega,\\ { f_{i}(x)|_{\partial\Omega}=C_{i} [g_{i}(t)-\widetilde {D}_{\lambda}^{i}{\sum_{j=i+1}^{k-i-1}} \mathbf{E}_{\lambda}^{j}f_{j}(x)|_{\partial\Omega} ]},\quad {t}\in\partial\Omega. \end{array}\displaystyle \right . $$
(20)

Proof

First, by Theorem 4.4, we can see that

$$\widetilde{D}_{\lambda}^{k}\Phi(x)=0. $$

Then, for \(0\leq i\leq k-1\), Lemma 4.3 implies that

$$\begin{aligned} \widetilde{D}_{\lambda}^{i}\Phi(x) = {\widetilde{D}_{\lambda}^{i} \Biggl(f_{0}(x) +{\sum_{j=1}^{k-1}} \mathbf{E}_{\lambda}^{j}f_{j}(x) \Biggr)} = {\frac{1}{C_{i}}f_{i}(x)+\widetilde{D}_{\lambda }^{i}{ \sum_{j=i+1}^{k-1-i}}\mathbf{E}_{\lambda}^{j}f_{j}(x)}. \end{aligned}$$

Letting \({x}\rightarrow t\), the formulas in (20) give \(\widetilde {D}_{\lambda}^{i}\Phi|_{\partial\Omega}=g_{i}(t)\), \(i=0,\ldots,k-1\), which implies the conclusion. □

References

  1. Obolashvili, E: Partially Differential Equations in Clifford Analysis. Pitman Monographs and Surveys in Pure and Applied Mathematics, vol. 96 (1999)

    Google Scholar 

  2. Obolashvili, E: Higher Order Partial Differential Equations in Clifford Analysis: Effective Solutions to Problems. Progress in Mathematical Physics, vol. 28. Birkhäuser, Boston (2003)

    Book  Google Scholar 

  3. Brackx, F, Delanghe, R, Sommen, F: Clifford Analysis. Res. Notes Math. Pitman, London (1982)

    MATH  Google Scholar 

  4. Huang, S, Qiao, YY, Wen, GC: Real and Complex Clifford Analysis. Springer, New York (2005)

    Google Scholar 

  5. Bernstein, S: On the index of Clifford algebra valued singular integral operators and the left linear Riemann problem. Complex Var. Theory Appl. 35, 33-64 (1998)

    Article  MATH  Google Scholar 

  6. Xu, ZY: Boundary value problems and function theory for Spin-invariant differential operators. Ph.D thesis, State University of Ghent (1989)

  7. Xu, ZY, Zhou, C: On boundary value problems of Riemann-Hilbert type for monogenic functions in a half space of \(R^{m}\), \(m\geq2\). Complex Var. Theory Appl. 22, 181-193 (1993)

    Article  MATH  MathSciNet  Google Scholar 

  8. Almansi, E: Sull’integrazione dell’equazione differenziale \(\Delta^{2m}u=0\). Ann. Mat. Pura Appl. 3(2), 1-51 (1899)

    Article  Google Scholar 

  9. Aronszajn, N, Creese, TM, Lipkin, LJ: Polyharmonic Functions. Oxford Mathematics Monographs. Clarendon, Oxford (1983)

    MATH  Google Scholar 

  10. Ryan, J: Iterated Dirac operators in \(C^{n}\). Z. Anal. Anwend. 9(5), 385-401 (1990)

    MATH  Google Scholar 

  11. Malonek, H, Ren, GB: Almansi-type theorems in Clifford analysis. Math. Methods Appl. Sci. 25, 1541-1552 (2002)

    Article  MATH  MathSciNet  Google Scholar 

  12. Nicolescu, M: Les fonctions polyharmoniques. Hermann, Paris (1936)

    Google Scholar 

  13. Karachik, VV: Normalized system of functions with respect to the Laplace operator and its applications. J. Math. Anal. Appl. 287, 577-592 (2003)

    Article  MATH  MathSciNet  Google Scholar 

Download references

Acknowledgements

This work was supported by the NNSF of China under Grant No. 11426082.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Hongfen Yuan.

Additional information

Competing interests

The author declares that they have no competing interests.

Author’s contributions

The author read and approved the final manuscript.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yuan, H. Boundary value problems for modified Dirac operators in Clifford analysis. Bound Value Probl 2015, 158 (2015). https://doi.org/10.1186/s13661-015-0426-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13661-015-0426-1

MSC

Keywords