Computational enzymology for degradation of chemical warfare agents: promising technologies for remediation processes

1.   Introduction

Mathematical models of physics, chemistry, ecology, physiology, psychology, engineering and social sciences have been governed by differential equation and difference equation. With the development of computers, compared with continuous-time model, discrete-time models described by difference equations are better formulated and analyzed in the past decades. Recently laser model has vast application in medical sciences, industries, highly security areas in army ^{[1,2,3,4,5,6,7]}. Laser, whose basic principal lies on the Einstein theory of light proposed in 1916, is a device that produces intense beam of monochromatic and coherent light. Since then it is developed by Gordon Gould in 1957. In 1960, the first working ruby laser was invented by Theodore Maiman. Laser light is coherent, highly directional and monochromatic which makes it different from ordinary light. The working principle of laser is based on the spontaneous absorption, spontaneous emission, stimulated emission, and population inversion are essential for the laser formation. The readers can refer to ^[8,9,10]. For instance, Hakin ^[11] proposed a simple continuous-time laser model in 1983. Khan and Sharif ^[12] proposed a discrete-time laser model and studied extensively dynamical properties about fixed points, the existence of prime period and periodic points, and transcritical bifurcation of a one-dimensional discrete-time laser model in $R^+$.

In fact, the identification of the parameters of the model is usually based on statistical method, starting from data experimentally obtained and on the choice of some method adapted to the identification. These models, even the classic deterministic approach, are subjected to inaccuracies (fuzzy uncertainty) that can be caused by the nature of the state variables, by parameters as coefficients of the model and by initial conditions. In fact, fuzzy difference equation is generation of difference equation whose parameters or the initial values are fuzzy numbers, and its solutions are sequences of fuzzy numbers. It has been used to model a dynamical systems under possibility uncertainty. Due to the applicability of fuzzy difference equation for the analysis of phenomena where imprecision is inherent, this class of difference equation is a very important topic from theoretical point of view and also its applications. Recently there has been an increasing interest in the study of fuzzy difference equations (see ^{[13,14,15,16,17,18,19,20,21,22,23,24,25]}).

In this paper, by virtue of the theory of fuzzy sets, we consider the following discrete-time laser model with fuzzy uncertainty parameters and initial conditions.

where $x_n$ is the number of laser photon at the $n$th time, $A, B, C, H$ and the initial value $x_0$ are positive fuzzy numbers.

The main aim of this work is to study the existence of positive solutions of discrete-time laser model (1.1). Furthermore, according to a generation of division (g-division) of fuzzy numbers, we derive some conditions so that every positive solution of discrete-time laser model (1.1) is bounded. Finally, under some conditions we prove that discrete-time laser model (1.1) has a fixed point $0$ which is asymptotically stable, and a unique positive fixed point $x$.

2.   Preliminaries and definitions

Firstly, we give the following definitions and lemma needed in the sequel.

Definition 2.1.^[26] $u: R\rightarrow [0, 1]$ is said to be a fuzzy number if it satisfies conditions (i)–(iv) written below:

(i) $u$ is normal, i. e., there exists an $x\in R$ such that $u(x) = 1$;

(ii) $u$ is fuzzy convex, i. e., for all $t\in[0, 1]$ and $x_1, x_2\in R$ such that

(iii) $u$ is upper semicontinuous;

(iv) The support of $u$, $\mbox{supp}u = \overline{\bigcup_{\alpha\in(0, 1]}[u]^{\alpha}} = \overline{\{x:u(x) > 0\}}$ is compact.

For $\alpha\in(0, 1]$, the $\alpha-$cuts of fuzzy number $u$ is $[u]^\alpha = \{x\in R: u(x)\geq\alpha\}$ and for $\alpha = 0$, the support of $u$ is defined as $\mbox{supp}u = [u]^0 = \overline{\{x\in R|u(x) > 0\}}$.

A fuzzy number can also be described by a parametric form.

Definition 2.2. ^[26] A fuzzy number $u$ in a parametric form is a pair $(u_l, u_r)$ of functions $u_l, u_r, 0\le \alpha\le 1$, which satisfies the following requirements:

(i) $u_l(\alpha)$ is a bounded monotonic increasing left continuous function,

(ii) $u_r(\alpha)$ is a bounded monotonic decreasing left continuous function,

(iii) $u_l(\alpha)\le u_r(\alpha), 0\le \alpha\le 1.$

A crisp (real) number $x$ is simply represented by $(u_l(\alpha), u_r(\alpha)) = (x, x), 0\le \alpha\le 1.$ The fuzzy number space $\{(u_l(\alpha), u_r(\alpha))\}$ becomes a convex cone $E^1$ which could be embedded isomorphically and isometrically into a Banach space (see ^[26]).

Definition 2.3.^[26] The distance between two arbitrary fuzzy numbers $u$ and $v$ is defined as follows:

It is clear that $(E^1, D)$ is a complete metric space.

Definition 2.4.^[26] Let $u = (u_l(\alpha), u_r(\alpha)), v = (v_l(\alpha), v_r(\alpha))\in E^1, 0\le \alpha\le 1,$ and $k\in R.$ Then

(i) $u = v$ iff $u_l(\alpha) = v_l(\alpha), u_r(\alpha) = v_r(\alpha)$,

(ii) $u+v = (u_l(\alpha)+v_l(\alpha), u_r(\alpha)+v_r(\alpha))$,

(iii) $u-v = (u_l(\alpha)-v_r(\alpha), u_r(\alpha)-v_l(\alpha))$,

(iv) $ku = \left\{ \begin{array}{ll} (ku_l(\alpha), ku_r(\alpha)), & k\ge 0;\\ (ku_r(\alpha), ku_l(\alpha)), & k < 0, \end{array} \right.$

(v) $uv = (\min\{u_l(\alpha)v_l(\alpha), u_l(\alpha)v_r(\alpha), u_r(\alpha)v_l(\alpha), u_r(\alpha)v_r(\alpha)\}, \max\{u_l(\alpha)v_l(\alpha), u_l(\alpha)v_r(\alpha), u_r(\alpha)v_l(\alpha), u_r(\alpha)v_r(\alpha)\})$.

Definition 2.5. ^[27] Suppose that $u, v\in E^1$ have $\alpha$-cuts $[u]^\alpha = [u_{l, \alpha}, u_{r, \alpha}], [v]^\alpha = [v_{l, \alpha}, v_{r, \alpha}]$, with $0\notin [v]^\alpha, \forall\alpha\in[0, 1]$. The generation of division (g-division) $\div_g$ is the operation that calculates the fuzzy number $s = u\div_g v$ having level cuts $[s]^\alpha = [s_{l, \alpha}, s_{r, \alpha}]$(here ${[u]^\alpha}^{-1} = [1/u_{r, \alpha}, 1/u_{l, \alpha}]$) defined by

provided that $s$ is a proper fuzzy number $s_{l, \alpha}$ is nondecreasing, $s_{r, \alpha}$ is nonincreasing, $s_{l, 1}\le s_{r, 1}$.

Remark 2.1. According to ^[27], in this paper the fuzzy number is positive, if $u\div_g v = s\in E^1$ exists, then the following two cases are possible

Case I. if $u_{l, \alpha}v_{r, \alpha}\le u_{r, \alpha}v_{l, \alpha}, \forall\alpha\in[0, 1],$ then $s_{l, \alpha} = \frac{u_{l, \alpha}}{v_{l, \alpha}}, s_{r, \alpha} = \frac{u_{r, \alpha}}{v_{r, \alpha}},$

Case II. if $u_{l, \alpha}v_{r, \alpha}\ge u_{r, \alpha}v_{l, \alpha}, \forall\alpha\in[0, 1],$ then $s_{l, \alpha} = \frac{u_{r, \alpha}}{v_{r, \alpha}}, s_{r, \alpha} = \frac{u_{l, \alpha}}{v_{l, \alpha}}.$

Definition 2.6. ^[26] A triangular fuzzy number (TFN) denoted by $A$ is defined as $(a, b, c)$ where the membership function

The $\alpha-$cuts of $A = (a, b, c)$ are described by $[A]^\alpha = \{x\in R: A(x)\ge \alpha\} = [a+\alpha(b-a), c-\alpha(c-b)] = [A_{l, \alpha}, A_{r, \alpha}]$, $\alpha\in[0, 1]$, it is clear that $[A]^\alpha$ are closed interval. A fuzzy number $A$ is positive if $\mbox{supp} A\subset(0, \infty).$

The following proposition is fundamental since it characterizes a fuzzy set through the $\alpha$-levels.

Proposition 2.1.^[26] If $\{A^\alpha: \alpha\in[0, 1]\}$ is a compact, convex and not empty subset family of $R^n$ such that

(i) $\overline{\bigcup A^\alpha}\subset A^0.$

(ii) $A^{\alpha_2}\subset A^{\alpha_1}$ if $\alpha_1\le\alpha_2.$

(iii) $A^\alpha = \bigcap_{k\ge 1}A^{\alpha_k}$ if $\alpha_k\uparrow\alpha > 0.$

Then there is $u\in E^n (E^n$ denotes $n$ dimensional fuzzy number space) such that $[u]^\alpha = A^\alpha$ for all $\alpha\in(0, 1]$ and $[u]^0 = \overline{\bigcup_{0 < \alpha\le 1}A^\alpha}\subset A^0.$

The fuzzy analog of the boundedness and persistence (see ^[15,16]) is as follows:

Definition 2.7. A sequence of positive fuzzy numbers $(x_n)$ is persistence (resp. bounded) if there exists a positive real number $M$ (resp. $N$) such that

A sequence of positive fuzzy numbers $(x_n)$ is bounded and persistence if there exist positive real numbers $M, N > 0$ such that

A sequence of positive fuzzy numbers $(x_n), n = 1, 2, \cdots$, is an unbounded if the norm $\|x_n\|, n = 1, 2, \cdots,$ is an unbounded sequence.

Definition 2.8. $x_n$ is a positive solution of (1.1) if $(x_n)$ is a sequence of positive fuzzy numbers which satisfies (1.1). A positive fuzzy number $x$ is called a positive equilibrium of (1.1) if

Let $(x_n)$ be a sequence of positive fuzzy numbers and $x$ is a positive fuzzy number, $x_n \rightarrow x$ as $n\rightarrow\infty$ if $\lim_{n\rightarrow\infty}D(x_n, x) = 0$.

Lemma 2.1. ^[26] Let $f: R^+\times R^+\times R^+\times R^+\rightarrow R^+$ be continuous, $A, B, C, D$ are fuzzy numbers. Then

3.   Main results

3.1.   Existence of positive solution of discrete-time laser model

Firstly we give the existence of positive solutions of discrete-time laser model with fuzzy environment.

Theorem 3.1. Let parameters $A, B, C, H$ and initial value $x_0$ of discrete-time laser model be fuzzy numbers. Then, for any positive fuzzy number $x_0$, there exists a unique positive solution $x_n$ of discrete-time laser model with initial conditions $x_0$.

Proof. The proof is similar to those of Proposition 2.1 ^[25]. So we omit the proof of Theorem 3.1.

Noting Remark 2.1, taking $\alpha$-cuts, one of the following two cases holds

Case I

Case II

3.2.   Dynamics of discrete-time laser model (1.1)

To study the dynamical behavior of the positive solutions of discrete-time laser model (1.1), according to Definition 2.5, we consider two cases.

First, if Case I holds true, we need the following lemma.

Lemma 3.1. Consider the following difference equation

where $a\in(0, 1), b, c, h\in(0, +\infty), y_0\in(0, +\infty)$, then the following statements are true:

(i) Every positive solution of (3.3) satisfies

(ii) The equation has a fixed point $y^* = 0$ if $b\le (1-a)h.$

(iii) The equation has two fixed points $y^* = 0, y^* = \frac{b}{c(1-a)}-\frac{h}{c}$ if $b > (1-a)h$.

Proof. (i) Let $y_n$ be a positive solution of (3.3). It follows from (3.3) that, for $n\ge 0,$

From which we have

This completes the proof of (i).

If $y^*$ is fixed point of (3.3), i.e., $y_n = y^*.$ So from (3.3), we have

After some manipulation, from (3.5), we can get

From (3.3), we can summarized the existence of fixed points as follows

(ii) If $b < h(1-a)$, then $y^* = \frac{b}{c(1-a)}-\frac{h}{c}$ is not a positive number. And hence if $b\le h(1-a)$ then (3.3) has a fixed point $y^* = 0$.

(iii) If $b > h(1-a)$, then $y^* = \frac{b}{c(1-a)}-\frac{h}{c}$ is a positive number. And hence if $b > h(1-a)$ then (3.3) has two fixed points $y^* = 0, y^* = \frac{b}{c(1-a)}-\frac{h}{c}$.

Proposition 3.1. The following statements hold true

(i) The fixed point $y^* = 0$ of (3.3) is stable if $b < (1-a)h$.

(ii) The fixed point $y^* = 0$ of (3.3) is unstable if $b > (1-a)h$.

(iii) The fixed point $y^* = 0$ of (3.3) is non-hyperbolic if $b = (1-a)h$.

Proof. From (3.3), let

From (3.7), it follows that

From (3.8), it can get

Therefore from (3.9), it can conclude that $y^* = 0$ is stable if $b < (1-a)h$, unstable if $b > (1-a)h$, non-hyperbolic if $b = (1-a)h$.

Proposition 3.2. The fixed point $y^* = \frac{b}{c(1-a)}-\frac{h}{c}$ of (3.3) is globally asymptotically stable if $b > (1-a)h$.

Proof. From (3.8), it can get

Therefore from (3.10), it can conclude that, if $a+\frac{h(1-a)^2}{b} < 1$, i.e., $b > (1-a)h$ then the fixed point $y^* = \frac{b}{c(1-a)}-\frac{h}{c}$ is stable.

On the other hand, it follows from (3.4) that $(y_n)$ is bounded. And from (3.8), we have $f'(y) > 0$. Namely $(y_n)$ is monotone increasing. So we have

Therefore the fixed point $y^* = \frac{b}{c(1-a)}-\frac{h}{c}$ is globally asymptotically stable.

Proposition 3.3. The fixed point $y^* = 0$ of (3.3) is globally asymptotically stable if $b < (1-a)h$.

Proof. From (3.3), we can get that

From (3.12), it follows that

Since $b < (1-a)h$, then $\lim_{n\rightarrow\infty}y_n = 0$. Therefore the fixed point $y^* = 0$ of (3.3) is globally asymptotically stable.

Theorem 3.2. Consider discrete-time laser model (1.1), where $A, B, C, H$ and initial value $x_0$ are positive fuzzy numbers. There exists positive number $N_A$, $\forall \alpha\in(0, 1]$, $A_{r, \alpha} < N_A < 1$. If

and

Then (1.1) has a fixed point $x^* = 0$ which is globally asymptotically stable.

Proof. Since $A, B, C, H$ are positive fuzzy numbers and (3.14) holds true, taking $\alpha$-cuts of model (1.1) on both sides, we can have the following difference equation system with parameters

Since (3.15) holds true, using Proposition 3.2, we can get

On the other hand, let $x_n = x^*$, where $[x_n]^\alpha = [L_{n, \alpha}, R_{n, \alpha}] = [L_\alpha, R_\alpha] = [x^*]^\alpha, \alpha\in(0, 1]$, be the fixed point of (1.1). From (3.16), one can get

Since (3.15) is satisfied, from (3.16), it follows that

Therefore, by virtue of Proposition 3.3, the fixed point $x^* = 0$ is globally asymptotically stable.

Theorem 3.3. Consider discrete-time laser model (1.1), where $A, B, C, H$ and initial value $x_0$ are positive fuzzy numbers, there exists positive number $N_A$, $\forall \alpha\in(0, 1]$, $A_{r, \alpha} < N_A < 1$. If (3.14) holds true, and

then the following statements are true.

(i) Every positive solution of (1.1) is bounded.

(ii) Equation (1.1) has a unique positive fixed point $x^*$ which is asymptotically stable.

Proof. (i) Since $A, B, C, H$ and $x_0$ are positive fuzzy numbers, there exist positive constants $M_A, N_A, M_B,$ $N_B, M_C, N_C, M_H, N_H, M_0, N_0$ such that

Using (i) of Lemma 3.1, we can get that

From (3.21) and (3.22), we have that for $\alpha\in(0, 1]$

where $N = \frac{N_B}{M_C(1-N_A)}+N_0.$ From (3.22), we have for $n\ge 1, \bigcup_{\alpha\in(0, 1]}[L_{n, \alpha}, R_{n, \alpha}]\subset (0, N]$, so $\overline{\bigcup_{\alpha\in(0, 1]}[L_{n, \alpha}, R_{n, \alpha}]}\subset (0, N].$ Thus the proof of (i) is completed.

(ii) We consider system (3.18), then the positive solution of (3.18) is given by

Let $x_n$ be a positive solution of (1.1) such that $[x_n]^\alpha = [L_{n, \alpha}, R_{n, \alpha}], \alpha\in(0, 1], n = 0, 1, 2, \cdots.$ Then applying Proposition 3.2 to system (3.16), we have

From (3.23) and (3.25), we have, for $0 < \alpha_1 < \alpha_2 < 1,$

Since $A_{l, \alpha}, A_{r, \alpha}, B_{l, \alpha}, B_{r, \alpha}, C_{l, \alpha}, C_{r, \alpha}, H_{l, \alpha}, H_{r, \alpha}$ are left continuous. It follows from (3.24) that $L_\alpha, R_\alpha$ are also left continuous.

From (3.21) and (3.24), it follows that

Therefore (3.27) implies that $[L_\alpha, R_\alpha]\subset[c, d],$ and so $\bigcup_{\alpha\in(0, 1]}[L_\alpha, R_\alpha]\subset[c, d].$ It is clear that

So from Definition 2.2, (3.24), (3.26), (3.28) and since $L_\alpha, R_\alpha, \alpha\in(0, 1]$ determine a fuzzy number $x^*$ such that

Suppose that there exists another positive fixed point $\bar{x}$ of (1.1), then there exist functions $\overline{L}_\alpha, \overline{R}_\alpha: (0, 1)\rightarrow(0, \infty)$ such that

From which we have

So $\overline{L}_\alpha = L_\alpha, \overline{R}_\alpha = R_\alpha, \alpha\in(0, 1]$. Hence $\bar{x} = x^*$, namely $x^*$ is a unique positive fixed point of (1.1).

From (3.25), we have

Namely, every positive solution $x_n$ of (1.1) converges the unique fixed point $x^*$ with respect to $D$ as $n\rightarrow\infty.$ Applying Proposition 3.2, it can obtain that the positive fixed point $x^*$ is globally asymptotically stable.

Secondly, if Case II holds true, it follows that for $n\in\{0, 1, 2, \cdots\}, \alpha\in(0, 1]$

We need the following lemma.

Lemma 3.2. Consider the system of difference equations

where $a_i\in(0, 1), b_i, c_i, h_i\in(0, +\infty) (i = 1, 2), y_0, z_0\in(0, +\infty)$. If

and

Then the following statements are true.

(i) Every positive solution $(y_n, z_n)$ of (3.32) satisfy

(ii) System (3.32) has fixed point $(y, z) = (0, 0)$ which is asymptotically stable.

(iii) System (3.32) has a unique fixed point

where $K = \frac{b_1b_2-h_1h_2(1-a_1)(1-a_2)}{(1-a_1)(1-a_2)}.$

Proof. (i) Let $(y_n, z_n)$ be a positive solution of (3.32). It follows from (3.32) that, for $n\ge 0$,

From which, we have

This completes the proof of (i).

(ii) It is clear that $(0, 0)$ is a fixed point of (3.32). We can obtain that the linearized system of (3.32) about the fixed point $(0, 0)$ is

where $X_n = (x_n, y_n)^T$ and

Thus the characteristic equation of (3.37) is

Since (3.33) and (3.34) hold true, we have

By virtue of Theorem 1.3.7 ^[28], we obtain that the fixed point $(0, 0)$ is asymptotically stable.

(iii) Let $(y_n, z_n) = (y, z)$ be fixed point of (3.32). We consider the following system

It is clear that the positive fixed point $(y, z)$ can be written by (3.36).

Theorem 3.4. Consider the difference Eq (1.1), where $A, B, C, H$ are positive fuzzy numbers. There exists positive number $N_A, \forall \alpha\in(0, 1], A_{r, \alpha} < N_A < 1.$ If

and

Then the following statements are true

(i) Every positive solution of (1.1) is bounded.

(ii) The Eq (1.1) has a fixed point $0$ which is globally asymptotically stable.

(iii) The Eq (1.1) has a unique positive fixed point $x$ such that

where

Proof. (i) Let $x_n$ be a positive solution of (1.1). Applying (i) of Lemma 3.2, we have

Next, the proof is similar to (i) of Theorem 3.3. So we omit it.

(ii) The proof is similar to those of Theorem 3.2. We omit it.

(iii) Let $x_n = x$ be a fixed point of (1.1), then

Taking $\alpha$-cuts on both sides of (3.46), since (3.41) holds true, one gets the following system

From which we obtain that

Next, we can show that $L_\alpha, R_\alpha$ constitute a positive fuzzy number $x$ such that $[x]^\alpha = [L_\alpha, R_\alpha], \alpha\in(0, 1].$ The proof is similar to (ii) of Theorem 3.3. We omit it.

Remark 3.1. In dynamical system model, the parameters of model derived from statistic data with vagueness or uncertainty. It corresponds to reality to use fuzzy parameters in dynamical system model. Compared with classic discrete time laser model, the solution of discrete time fuzzy laser model is within a range of value (approximate value), which are taken into account fuzzy uncertainties. Furthermore the global asymptotic behaviour of discrete time laser model are obtained in fuzzy context. The results obtained is generation of discrete time Beverton-Holt population model with fuzzy environment ^[25].

4.   Numerical simulations

In this section, we give two numerical examples to verify the effectiveness of theoretic results obtained.

Example 4.1. Consider the following fuzzy discrete time laser model

we take $A, B, C, H$ and the initial values $x_0$ such that

From (4.2), we get

From (4.3) and (4.4), we get

Therefore, it follows that

From (4.1), it results in a coupled system of difference equations with parameter $\alpha,$

Therefore, it is clear that $A_{r, \alpha} < 1, \forall\alpha\in(0, 1],$ (3.14) and (3.15) hold true. so from Theorem 3.2, we have that every positive solution $x_n$ of Eq (4.1) is bounded In addition, from Theorem 3.2, Eq (4.1) has a fixed point $0$. Moreover every positive solution $x_n$ of Eq (4.1) converges the fixed point $0$ with respect to $D$ as $n\rightarrow \infty.$ (see Figures 1–3).

Example 4.2. Consider the following fuzzy discrete time laser model (4.1). where $A, C, H$ and the initial values $x_0$ are same as Example 4.1.

From (4.9), we get

Therefore, it follows that

It is clear that (3.20) is satisfied, so from Theorem 3.3, Eq (4.1) has a unique positive equilibrium $\overline{x} = (0.341, 1.667, 3)$. Moreover every positive solution $x_n$ of Eq (4.1) converges the unique equilibrium $\overline{x}$ with respect to $D$ as $n\rightarrow \infty.$ (see Figures 4–6)

5.   Conclusions

In this work, according to a generalization of division (g-division) of fuzzy number, we study the fuzzy discrete time laser model $x_{n+1} = Ax_n+\frac{Bx_{n}}{Cx_{n}+H}$. The existence of positive solution and qualitative behavior to (1.1) are investigated. The main results are as follows

(i) Under Case I, the positive solution is bounded if $B_{l, \alpha} < H_{l, \alpha}(1-A_{l, \alpha}), B_{r, \alpha} < H_{r, \alpha}(1-A_{r, \alpha}), \alpha\in(0, 1]$. Moreover system (1.1) has a fixed point $0$ which is globally asymptotically stable. Otherwise, if $B_{l, \alpha}\ge H_{l, \alpha}(1-A_{l, \alpha}), B_{r, \alpha}\ge H_{r, \alpha}(1-A_{r, \alpha}), \alpha\in(0, 1]$. Then system (1.1) has a unique positive fixed point $x^*$ which is asymptotically stable.

(ii) Under Case II, if $A_{l, \alpha}+A_{r, \alpha} < 1$ and $B_{l, \alpha}B_{r, \alpha} > H_{l, \alpha}H_{r, \alpha}(1-A_{l, \alpha})(1-A_{r, \alpha}), \alpha\in(0, 1]$, then the positive solution is bounded. Moreover system (1.1) has a unique positive fixed point $x$ and fixed point $0$ which is global asymptotically stable.

Acknowledgments

The authors would like to thank the Editor and the anonymous Reviewers for their helpful comments and valuable suggestions to improve the paper. The work is partially supported by National Natural Science Foundation of China (11761018), Scientific Research Foundation of Guizhou Provincial Department of Science and Technology([2020]1Y008, [2019]1051), Priority Projects of Science Foundation at Guizhou University of Finance and Economics (2018XZD02), and Scientific Climbing Programme of Xiamen University of Technology(XPDKQ20021).

Conflict of interest

The authors declare that they have no competing interests.