Numerical simulation of the space fractional $(3+1)$-dimensional Gray-Scott models with the Riesz fractional derivative

1.   Introduction

Let $\mathbb{D}$ be the unit disk in the complex plane $\mathbb{C}$ centered at the origin and denote $\mathbb{T} = \partial\mathbb{D}.$ The symbol $\mathcal{C}^{n}(\Omega)$ stands for the class of all complex-valued $n$-times continuously differentiable functions from the domain $\Omega$ into $\mathbb{C}$. In particular, we let $\mathcal{C}(\Omega) = \mathcal{C}^{0}(\Omega).$ A real-valued function $u\in\mathcal{C}^{2}(\Omega)$ is called real harmonic if it satisfies the following Laplace's equation:

A complex-valued function $\omega = u + iv$ is harmonic if both $u$ and $v$ are real harmonic. We refer the readers to ^[12] for an excellent discussion on harmonic mappings in the plane.

Suppose that $g\in\mathcal{C}(\mathbb{D}).$ Then, it is well known that the solution of the Poisson's equation $\Delta\omega = g$ in $\mathbb{D}$ satisfying the boundary condition $\omega|_{\mathbb{T}} = f\in L^{1}(\mathbb{T})$ is given by

where

and $m$ denotes the Lebesgue measure in the plane.

The hyperbolic metric on $\mathbb{D},$ the Gaussian curvature of which is equal to -4, is given by $|dz|/(1-|z|^{2}).$ Then, the hyperbolic distance between two points $z_{1}, z_{2}\in \mathbb{D}$ is defined by

where the infimum is taken over all the rectifiable curves $\gamma$ which connect $z_{1}$ and $z_{2}$ with $\gamma\in\mathbb{D}.$ It is well known that

Also, we can consider the quasihyperbolic metric $|dz|/(1-|z|)$ on $\mathbb{D}.$ Then, it gives rise to the quasihyperbolic distance between two points $z_{1}, z_{2}\in \mathbb{D},$ which is defined by

where the infimum is taken over all rectifiable curves $\gamma$ which connect $z_{1}$ and $z_{2}$ with $\gamma\in\mathbb{D}.$ Actually, we can consider the general $q$-metric $|dz|/(1-|z|^{q})$ on $\mathbb{D}, q\in[1,2],$ which leads to the following $q$-pseudo distance

where the infimum is taken over all rectifiable curves $\gamma$ which connect $z_{1}$ and $z_{2}$ with $\gamma\in\mathbb{D}.$ For completeness, we will prove that this $q$-pseudo distance is a distance in Lemma 7.

A self-mapping $f$ of $\mathbb{D}$ is said to be Lipschitz continuous with respect to the $q$-distance if there exists a constant $L > 0,$ such that the inequality

holds for any $z_{1}, z_{2}\in \mathbb{D}.$

The classic Schwarz-Pick lemma plays an important role in complex analysis. It is stated as follows: If $f$ is a holomorphic self-mapping of $\mathbb{D},$ then we have

or

The Schwarz-Pick lemmas (1.6) or (1.7) have many generalizations. See the references ^{[1,6,13,30,31,32,34,35,39,46]} on this topic. The first work of this paper is to extend (1.7) to the case for the $q$-distance. It is presented as follows:

Theorem 1. Suppose that $q\in [1,2].$ If $f$ is a holomorphic self-mapping of $\mathbb{D}$ satisfying $f(0) = 0,$ then the inequality

holds for any $z_{1}, z_{2}\in\mathbb{D}.$

Suppose that $f$ satisfies the following the Poisson differential inequality:

where $|Df| = |f_{z}|+|f_{\overline{z}}|.$ This kind of mapping contains many classical mappings. For example, it includes a holomorphic mapping, a harmonic mapping, the mapping satisfying Poisson's equation $\Delta f = g$ (where $g$ is continuous on $\mathbb{\overline{D}}$), the harmonic mapping between Riemannian surfaces satisfying some inequalities and so on. In addition, those mappings that satisfy (1.9) are also closely related to the following partial differential equations

where $Q = (Q_{1}, Q_{2}, \ldots, Q_{m}),$ $Q_{j}$ are quadratic polynomial on $\partial u_{i}/\partial x_{k},$ $i = 1, \ldots, m, k = 1, 2.$ This equation has a deep connection with the average curvature problem of surfaces in differential geometry, the conjugate isothermal coordinates of surfaces, and the Monge-Ampère equation and so on. Based on this backgrounds, the study of mappings that satisfy the Poisson differential inequality (1.9) has attracted much attention from researchers. See the articles ^{[4,18,20,42,45]} and the references therein. Recently, Kalaj ^[21] extended the results of Bernstein ^[4] and Heinz ^[18] to the case for spaces and applied the relevant results to the theory of harmonic quasiconformal mappings. The study of mappings satisfying the Poisson differential inequality that are also quasiconformal or quasiregular have also attracted many authors' interest, see ^[3,10,22].

Let $f$ be a quasiconformal self-mapping of the unit disk satisfying the following Poisson differential equation

In ^[23], Kalaj obtained the following result.

Lemma 1. ^[23] Suppose that $f$ is a $K$-quasiconformal self-mapping of $\mathbb{D}$ satisfying the Poisson differential equation (1.11), and that $f(0) = 0$. Then there exists a constant $C(B, K),$ which depends only on $B$ and $K,$ such that

The inequality (1.12) can be seen as a kind of Schwarz-Pick type inequality for the $K$-quasiconformal self-mapping of $\mathbb{D}$ satisfying the Poisson differential inequality (1.11). Generally, the condition for quasiconformality in Lemma 1 can not be removed. Recently, Zhong et al. ^[47] showed the following result.

Theorem 2. [47, Theorem 1.4] For a given $q\in\{1\}\cup[2, +\infty)$. Suppose that $f: \overline{\mathbb{D}}\rightarrow \overline{\mathbb{D}}$ is continuous on $\overline{\mathbb{D}},$ $f|_{\mathbb{D}}\in \mathcal{C}^{2}, f|_{\mathbb{T}}\in\mathcal{C}^{2}$ and $f(0) = 0$, and that it satisfies the following Poisson differential inequality

as well as satisfying that $\left|\frac{\partial^{2} f(e^{i\varphi})}{\partial \varphi^{2}}\right|\leq K,$ where $0 < a < 1/2, 0 < b, K < \infty.$ If

then, we have

In the next result of this paper, we will improve Theorem 2 by giving the explicit constant $L(a, b, K)$ and showing that the result is still valid when we consider all of the cases for $q > 0.$ Our result is presented as follows.

Theorem 3. For given $q > 0,$ suppose that $f: \overline{\mathbb{D}}\rightarrow \overline{\mathbb{D}}$ is continuous in $\overline{\mathbb{D}},$ $f|_{\mathbb{D}}\in \mathcal{C}^{2}, f|_{\mathbb{T}}\in\mathcal{C}^{2}$ and $f(0) = 0$, and that it satisfies the Poisson differential inequality $|\Delta f(z)|\leq a\cdot|D f(z)|^{2}+b$ and $\left|\frac{\partial^{2} f(e^{i\varphi})}{\partial \varphi^{2}}\right|\leq K,$ where $0 < a < 1/2, 0 < b, K < \infty.$ If, in addition, $48ab+640a < 9$ and $\max\{a, b\}\cdot\left(L^{2}(a, b, K)+1\right) < 2/\pi,$ then we have

where

As an application of Theorem 3, we provide the following more general form of [47, Corollary 1.5].

Corollary 1. For a given $q > 0,$ suppose that $f: \overline{\mathbb{D}}\rightarrow \overline{\mathbb{D}}$ is continuous in $\overline{\mathbb{D}},$ $f|_{\mathbb{D}}\in \mathcal{C}^{2}, f|_{\mathbb{T}}\in\mathcal{C}^{2}$ and $f(0) = 0$, and that it satisfies the Poisson differential inequality $|\Delta f(z)|\leq a\cdot|D f(z)|^{2}+b$ and $\left|\frac{\partial^{2} f(e^{i\varphi})}{\partial \varphi^{2}}\right|\leq K,$ where $0 < a < 1/2, 0 < b, K < \infty.$ If, in addition, $48ab+640a < 9$ and $\max\{a, b\}\cdot\left(L^{2}(a, b, K)+1\right) < 2/\pi,$ then we have

Proof. According to Lemma 9, we see that there exists a constant $L(a, b, K)$ such that the inequality

holds for any $z\in\mathbb{D}.$ Combining this with Theorem 3, we get

By the definition of $q$-distance, we derive that inequality (1.15) is valid. □

The classic boundary Schwarz lemma for the holomorphic mapping states the following:

Lemma 2. (^[14]) Let $f$ be a holomorphic self-mapping of $\mathbb{D}$ satisfying that $f(0) = 0.$ If $f$ is differentiable at the point $z = 1$ and satisfies that $f(1) = 1,$ then we have the follwoing:

$({\rm{i}})$ $f'(1)\geq1;$

$({\rm{ii}})$ $f'(1) = 1$ if and only if $f(z) = z.$

The boundary Schwarz lemma is a basic and important result in complex analysis. Its generalizations and applications have received widespread attention; see ^{[5,7,9,15,16,17,25,26,27,28,29,40,43,44]}. In ^[25], Lemma 2 was generalized to the following form:

Lemma 3. (^[25]) Let $f$ be a holomorphic self-mapping of $\mathbb{D}$ satisfying that $f(0) = 0.$ If $f$ is differentiable at the point $z = \alpha\in\mathbb{T}$ and satisfies that $f(\alpha) = \beta\in\mathbb{T}.$ Then we have the following

$({\rm{i}} ')$ $\overline{\beta} f'(\alpha)\alpha\geq1;$

$({\rm{ii}} ')$ $\overline{\beta} f'(\alpha)\alpha = 1$ if and only if $f(z) = e^{i\theta}z,$ where $e^{i\theta} = \beta\alpha^{-1}$ and $\theta\in\mathbb{R}.$

Recently, the extension of the boundary Schwarz lemma to harmonic mappings has also attracted scholars' attention. In ^[50], Zhu and Wang established the following harmonic version of the boundary Schwarz lemma.

Lemma 4. ([50, Theorem 1.2]) Suppose that $g\in \mathcal{C}(\overline{\mathbb{D}}).$ If $f\in \mathcal{C}^{2}(\mathbb{D})\bigcap\mathcal{C}(\mathbb{T})$ is a self-mapping of $\mathbb{D}$ satisfying the Poisson equation $\Delta f = g$ and $f(0) = 0,$ and if $f$ is differentiable at the point $z = 1$ and satisfies that $f(1) = 1,$ then we have

where $||g||_{\infty}: = \sup _{z\in\mathbb{D}}|g(z)|.$

This result was extended by Mohapatra et al. ^[37] and Mohapatra ^[38]. See the high dimensional version of the boundary Schwarz lemma for harmonic mappings in ^{[24,33,36,49]}. As a generalization of planar harmonic maps, those self-mappings that satisfy the Poisson differential inequality (1.11) can be naturally considered as analogous to the boundary Schwarz lemma. Unfortunately, those types of mappings generally are not associated with the boundary Schwarz lemma of types (1.16). This is because we have established the following results:

Theorem 4. There exists a self-mapping $f$ of $\mathbb{D},$ which is differentiable at $z = 1$ and satisfies that $f(0) = 0$ and $f(1) = 1$, in addition to the Poisson differential inequality $|\Delta f(z)|\leq a\cdot|D f(z)|^{2}.$ But it satisfies the equation ${\rm Re}(f_{x}(1)) = 0.$

According to Lemma 4, we see that there exists a self-mapping of $\mathbb{D}$ that satisfies the Poisson differential inequality (1.9), which is associated with a boundary Schwarz lemma of type (1.16). Therefore, we naturally ask the following question: which subset of the family

has boundary Schwarz lemma of type (1.16)? Next, we establish a result in this regard.

Theorem 5. Suppose that $f: \overline{\mathbb{D}}\rightarrow \overline{\mathbb{D}}$ is continuous in $\overline{\mathbb{D}},$ $f|_{\mathbb{D}}\in \mathcal{C}^{2}, f|_{\mathbb{T}}\in\mathcal{C}^{2}$ and $f(0) = 0$, and that it satisfies the Poisson differential inequality $|\Delta f(z)|\leq a\cdot|D f(z)|^{2}+b$ and $\left|\frac{\partial^{2} f(e^{i\varphi})}{\partial \varphi^{2}}\right|\leq K,$ where $0 < a < 1/2, 0 < b, K < \infty.$ If, in addition, $48ab+640a < 9$ and $\max\{a, b\}\cdot\left(L^{2}(a, b, K)+1\right) < 8/(3\pi),$ then we have

where

2.   Some preparations

L' Hospital's rule for monotonicity is important to derive some inequalities and it will be used to prove some results.

Lemma 5. ([2, Theorem 2]) Let $-\infty < a < b < +\infty$ and let $f, g:[a, b]\rightarrow \mathbb{R}$ be continuous functions that are differentiable on $(a, b)$ with $f(a) = g(a) = 0$ or $f(b) = g(b) = 0.$ Suppose that $g'(x)\neq0$ for each $x\in(a, b).$ If $f'/g'$ is increasing $($decreasing$)$ on $(a, b),$ then so is $f/g.$

The following lemma can be obtained in ^[48] by letting $x = t^{2}.$

Lemma 6. For any $q > 0,$ the inequalities

hold for any $x\in[0, 1).$

Remark 1. In the rest of this paper, we will always use the constant $c_{q},$ which is defined by

In Section 1, we defined the following general $q$-pseudo distance in $\mathbb{D}$:

where the infimum is taken over all rectifiable curves $\gamma$ in $\mathbb{D}$ connected $z_{1}$ and $z_{2}$. Next, we will prove that it is a distance.

Lemma 7. The function $d_{q}(\cdot, \cdot)$ defined by (2.3) is a distance.

Proof. Obviously, for any $z_{1}, z_{2}\in\mathbb{D},$ we have that $d_{q}(z_{1}, z_{2})\geq0$ and $d_{q}(z_{1}, z_{2}) = d_{q}(z_{2}, z_{1}).$ Hence, we only prove that if $d_{q}(z_{1}, z_{2}) = 0,$ then $z_{1} = z_{2}.$ This is because, by taking $x = |z|^{2}$ in Lemma 6, we get

Hence, by the definitions of hyperbolic distance and $q$-pseudo distance, we have

Therefore, it must be that $d_{h_{\mathbb{D}}}(z_{1}, z_{2}) = 0$ if $d_{q}(z_{1}, z_{2}) = 0.$ Since $d_{h_{\mathbb{D}}}(\cdot, \cdot)$ is a distance, we get that $z_{1} = z_{2}.$ This shows that the $q$-pseudo distance $d_{q}(\cdot, \cdot)$ defined in $\mathbb{D}$ is a distance. □

We also need the following general Schwarz lemma for harmonic self-mappings of the unit disk, which was obtained by Hethcote ^[19].

Lemma 8. ([19, Theorem 1]) Suppose that $f$ is a harmonic self-mapping of $\mathbb{D}.$ Then we have

In ^[18], Heinz obtained a gradient estimate for self-mappings of $\mathbb{D}$ satisfying the Poisson differential inequality. The upper bound he obtained was abstract. In what follows, we give an explicitly upper bound based on the proof in ^[18]. The disadvantage is that the range of $a, b$ becomes smaller.

Lemma 9. Suppose that $f: \overline{\mathbb{D}}\rightarrow \overline{\mathbb{D}}$ is continuous in $\overline{\mathbb{D}},$ $f|_{\mathbb{D}}\in \mathcal{C}^{2}, f|_{\mathbb{T}}\in\mathcal{C}^{2}$ and $f(0) = 0$, and that it satisfies the Poisson differential inequality $|\Delta f(z)|\leq a\cdot|D f(z)|^{2}+b$ and $\left|\frac{\partial^{2} f(e^{i\varphi})}{\partial \varphi^{2}}\right|\leq K,$ where $0 < a < 1/2, 0 < b, K < \infty.$ If $48ab+640a < 9,$ then the inequality

holds for any $z\in\mathbb{D},$ where

Proof. First, we calculate $c'_{10}(a, b, K)$ in [18, P. 239]. According to [18, Lemma 11] and given

in [18, P. 239] and [18, (2.2.45), P. 238], we get

Hence, by [18, (2.2.47), P. 239], we can take

Next, we will estimate the constant $c_{5}(a, b)$ in ^[18]. According to the proof of [18, Theorem 2], we can take $\alpha = 2, \; \beta = 2(1-a)$ and $\; \gamma = 2b$ in [18, (2.1.27)]. Therefore, the formula given by [18, (2.1.27)] can be rewritten as

where $\lambda$ and $\theta$ are constants satisfying that $0 < \lambda < \theta < 1.$ Now, we take $\theta = 2\lambda$ and assume that $0 < \lambda < 1/4.$ Since $\ln\frac{1}{\theta}\geq1-\theta$ and $\frac{a^{2}}{1-a}\leq a$ (as $a < 1/2$), we have

As the functions $\varphi_{1}(\lambda) = 8a+4ab\lambda+\frac{a(b+16)}{1-2\lambda}, \varphi_{2}(\lambda) = -(\lambda-1)^{2}$ are monotonically increasing on $(0, 1/4),$ we see that the inequality

holds for any $0 < \lambda < 1/4$ when $48ab+640a < 9.$ Especially, inequality (2.10) holds for any $0 < \lambda\leq1/8.$ Now, we take $\lambda_{0} = 1/8$ and $\theta_{0} = 1/4.$ Then by [18, (2.1.31), P. 227], we have

Hence, by the inequality given by[18, (2.1.31)], we get

Combining [18, P. 239] with (2.7) and (2.11), we have

This completes the proof. □

The proof of following estimate can be found in ^[8].

Lemma 10. (^[8]) Suppose that $g\in\mathcal{C}(\mathbb{D}).$ Then, $|G[g](z)|\leq\frac{(1-|z|^{2})\cdot||g||_{\infty}}{4},$ where $||g||_{\infty}: = \sup_{z\in\mathbb{D}}|g(z)|.$

Lemma 11. For any $t\in[0, 1),$ we have

Proof. Let $\varphi_{1}(t) = 1-\frac{4}{\pi}\cdot\arctan t$ and $\varphi_{2}(t) = 1-t, t\in[0, 1).$ Then, $\varphi_{1}(1) = \varphi_{2}(1) = 0$ and the function

is monotonically decreasing on $t\in[0, 1).$ By L' Hospital's rule for monotonicity ^[2], we see that the function $\varphi_{1}/\varphi_{2}$ is monotonically decreasing on $t\in[0, 1).$ Hence, by L'Hospital's rule, we get

This finishes the proof. □

The following result was proved by Ruscheweyh ^[41].

Lemma 12. (^[41]) Suppose that $f$ is a holomorphic self-mapping of $\mathbb{D}.$ Then,

In what follows, we generalize Lemma 12 to the following form.

Proposition 1. For any $q\in [1,2]$, suppose that $f$ is a holomorphic self-mapping from the unit disk $\mathbb{D}$ into itself; then, we have

where $n\in\{1, 2, \cdot\cdot\cdot\}$.

Proof. Let

Since $\frac{(1-t^{2})'}{(1-t^{q})'} = \frac{2}{q}\cdot t^{2-q}$ is monotonically increasing on $[0, 1),$ the function $\varphi$ is also monotonically increasing on $[0, 1).$ By Lindelöf's inequality, we have that $|f(z)|\leq\frac{|z|+|f(0)|}{1+|z|\cdot |f(0)|}, z\in\mathbb{D}.$ Therefore,

Combining this with Lemma 12, we get

This has finished the proof. □

Remark 2. If $q = 2,$ then Proposition 1 becomes Lemma 12.

3.   Proof of Theorem 1

Proof. By using Proposition 1 for $n = 1$ and $f(0) = 0,$ we get

Namely,

According to the definition of $q$-distance, we have that $d_{q}(f(z_{1}), f(z_{2}))\leq d_{q}(z_{1}, z_{2}).$ This finishes the proof. □

Remark 3. Considering the case for $n = 1$ and $f(0) = 0$ in Proposition 1, we get the following Schwarz-Pick lemma: Suppose that $f$ is a holomorphic self-mapping from unit disk $\mathbb{D}$ into itself with $f(0) = 0$; then, we have

From the proof of Proposition 1, we see that the inequality (3.2) is sharp. The equation holds if and only if $|f(z)| = |z|,$ which is subject to the condition if and only if $f(z) = e^{i\theta}\cdot z, \theta\in\mathbb{R}$ is true.

4.   Proof of Theorem 3

Proof. According to Lemma 9, we see that there exists a constant $L(a, b, K),$ such that the inequality

holds for any $z\in\mathbb{D},$ where

Now, let

where $h(z) = l(z)(|Df(z)|^{2}+1)$ and $||l||_{\infty}\leq\max\{a, b\}.$ For example, we can define $l$ as follows:

Since $f\in \mathcal{C}^{2}(\mathbb{D})$, we see that $h\in\mathcal{C}(\mathbb{D}).$ Hence, by Lemma 9, we get

By virtue of the solution of Poisson's equation, we have

where $f|_{\mathbb{T}}: = k.$ Since $P[k]$ is a harmonic self-mapping of $\mathbb{D}$, we get the following by Lemma 8

Combining Lemmas 6, 8, 10 and 11, we get

This completes the proof. □

5.   Proof of Theorem 4

Proof. Consider the function

then, $|f(z)| < 1, f(e^{i\theta}) = e^{i\theta}, f_{z}(z) = \frac{3-2|z|^{2}}{2}, f_{\overline{z}}(z) = -\frac{z^{2}}{2}, \Delta f(z) = 4f_{z\overline{z}}(z) = -4z$ and $|f_{z}(z)|+|f_{\overline{z}}(z)| = \frac{3-|z|^{2}}{2}.$ Therefore, $f$ is a self-mapping of $\mathbb{D}$ satisfying that $|\Delta f(z)|\leq4(Df(z))^{2}$ and $f(0) = 0.$ But ${\rm Re}(f_{x}(1)) = {\rm Re}(f_{z}(1)+f_{\overline{z}}(1)) = 0.$ □

6.   Proof of Theorem 5

Proof. Set $\Delta f(z) = h(z), z\in\mathbb{D},$ where $h(z) = l(z)(|Df(z)|^{2}+1)$ and $||l||_{\infty}\leq\max\{a, b\}.$ Since $h\in\mathcal{C}(\mathbb{D}),$ we get, by Lemma 9 that $||h||_{\infty}\leq\max\{a, b\}\cdot(L^{2}(a, b, K)+1).$ Then,

where $f|_{\mathbb{T}}: = k.$ Next, the method is similar to the one given in ^[50]. But, here, we do not demand that $h\in\mathcal{C}(\mathbb{\overline{D}}),$ as we only need $h\in\mathcal{C}(\mathbb{D}).$ By virtue of (6.1) and Lemmas 8 and 10, we have

Since $f$ is differential at $z = 1$ and satisfies that $f(1) = 1,$ we get

Since $|f(z)|^{2}\leq M^{2}(z),$

By letting $z = r\in(0, 1)$ and $r\rightarrow1^{-},$ we obtain

This finishes the proof. □

7.   Conclusions

By establishing the general Schwarz-Pick lemma, which is a generalization of the result of Ruscheweyh ^[41] and sharp when $n = 1$ and $f(0) = 0,$ we have obtained a Schwarz-Pick-type lemma for the holomorphic self-mapping of the unit disk with respect to the $q$-distance. We also have established the general Schwarz-Pick lemma for self-mappings of the unit disk satisfying the Poisson differential inequality. This obtained result improved the one in ^[47] by giving the explicit constant $L(a, b, K)$, and it showed that the result is still valid when we consider all of the cases for $q > 0.$ As an application, it has been proven that this mapping is Lipschitz continuous with respect to the $q$-distance under certain conditions; Moreover, the corresponding explicit Lipschitz constant has been given. Next, we investigated the boundary Schwarz lemma for self-mappings of the unit disk satisfying the Poisson differential inequality. We found that those types of mappings are generally not associated with the boundary Schwarz lemma. With some additional conditions, we have established a boundary Schwarz lemma for self-mappings of the unit disk satisfying the Poisson differential inequality.

Use of AI tools declaration

The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

Acknowledgments

The authors would like to express their sincere thanks to the referees for their great efforts to improve this paper. This work was supported by the Guangdong Basic and Applied Basic Research Foundation (No. 2022A1515110967; No. 2023A1515011809) and Research Foundation of Shenzhen Polytechnic University (No. 6023312032K).

Conflict of interest

The authors declare that they have no competing interests.