Existence and analysis of Hilfer-Hadamard fractional differential equations in RLC circuit models

1.   Introduction

It is widely accepted that the delayed dynamical model is an important tool for answering many issues in biology. Numerous academics focus extensively on the development of diverse predator-prey models to elucidate the dynamics and internal structure within biological populations. We can successfully regulate the concentrations of predators and preys in the natural world by investigating the various dynamical characteristics of predator-prey models. Numerous studies on predator-prey models have been submitted and released in recent years, and several outstanding studies have been discussed. For example, Din et al. ^[1] studied the discrete predator-prey interactions with chaos control and bifurcation analysis. Liu and Guo ^[2] examined the dynamics of a predator-prey system with nonlinear prey-taxis, using the Lyapunov-Schmidt reduction approach, the fold-Hopf singularity, Hopf bifurcation, and steady state bifurcation are examined in depth. Al-Kaff et al. ^[3] investigated bifurcation and chaos in a discrete predator-prey system using a linked logistic map, a wide range of the system's behavior is thoroughly examined in the study. Pita et al. ^[4] investigated a review of some recent advances in predator-prey models. For more extensive investigations, see ^[5,6,7,8,9].

In particular, Lotka-Volterra models are vital predator-prey models in biology. They play an important role in describing the correlation between predators and prey. During the past decades, many works on this topic have been reported. For example, Hsu et al. ^[10] studied three different species of omnivorous Lotka-Volterra food web models. Bunin ^[11] studied the Lotka-Volterra model for ecological communities exhibiting symmetry and connected it to other well-known models. Wu et al. ^[12] investigated the Grey Lotka-Volterra model and its application, effectively analyzing the link between the two variables and forecasting their values using the gray Lotka-Volterra model. Kloppers et al. ^[13] estimated the parameters of the Lotka-Volterra model from empirical data. Marasco et al. ^[14] found the development of market share dynamics through the application of Lotka-Volterra models, and Wang et al. ^[15] examined the problems with free boundaries in a Lotka-Volterra competition system. Zhou ^[16] investigated a two-species Lotka-Volterra competition system in one-dimensional advective environments. Cherniha et al. ^[17] studied a review and new findings regarding the construction and application of exact solutions to the diffusive Lotka-Volterra system.

In 2022, Prabir et al. ^[18] suggested the subsequent Lotka-Volterra mutualistic symbiotic relationship:

where $\mho_1$ and $\mho_2$ indicate the density of two logistically increasing competing species, $\mho_3$ represents the density of a predatory species of $\mho_2$, ${r_1}, {r_2}$ are the natural rates of growth for species $\mho_1$ and $\mho_2$, ${K_1}, {K_2}$ are the carrying capacity of species $\mho_1$ and $\mho_2$ respectively, ${\alpha _{12}}, {\alpha _{21}}$ represent the fight between two species $\mho_1$ and $\mho_2$, $\delta _{13}$ is the commensal coefficient of $\mho_1$ over $\mho_3$ species, $a$ is the frequency at which $\mho_3$ species attack $\mho_2$ species, $p$ is the refuge rate of $\mho_2$ species, $c$ refers to the rate of conservation of $\mho_2$ species, $b$ is the steady value for half saturation for the Holling type Ⅱ function, and $e$ is the rate of extinction for species $\mho_3$. One might refer to ^[18] for a more precise explanation of the meaning of system (1.1).

In order to make the suggested mathematical model nondimensional, we have now added a few variables and constants as $t = {r_1}T, {u_1} = \frac{{\mho_1}}{{{K_1}}}, {u_2} = \frac{{\mho_2}}{{{K_2}}}, {u_3} = \frac{{a{\mho_3}}}{{{r_1}{K_2}}}, {\gamma _{12}} = \frac{{{\alpha _{12}}K_2^2}}{{{K_1}}}$, $\gamma = \frac{{{\delta _{13}}{K_2}}}{a}, r = \frac{{{r_2}}}{{{r_1}}}, {\gamma _{21}} = \frac{{{\alpha _{21}}K_1^2}}{{{K_2}}}, {v_1} = \frac{b}{{{K_2}}}, {v_2} = \frac{e}{{{r_1}}}, {v_3} = \frac{{ac}}{{{r_1}}}$.

Then model (1.1) can be simplified to the following form

with initial conditions ${u_1}(0) \ge 0, {u_2}(0) \ge 0, {u_3}(0) \ge 0$.

In this instance we want to call your attention to the truth that the evolution of species frequently depends on both historical and present times; hence, delays must be incorporated into biological models. Notice that the development of the density of two logistically increasing competing species $\mho_1$ and $\mho_2$ rely on not only the current time but also the past time. In other words, the development of the density of two logistically increasing competing species $\mho_1$ and $\mho_2$ shall have feedback delays. Based on this idea, we assume that there are two self-feedback delays in model (1.2): One is the self-feedback delay from the first species $u_1$ to the first species $u_1$, and the other the self-feedback delay from the first species $u_2$ to the first species $u_2$. The points that follow delayed Lotka-Volterra commercially available symbiosis system may then be loosely formulated:

where $\theta > 0$ is a time delay and all the parameters ${\gamma _{12}}, \gamma, r, {\gamma _{21}}, {v_1}, {v_2}, {v_3}$ are positive constants. The initial conditions of system (1.3) are given by $u_i(s) = \phi_i(s)\geq0, i = 1, 2, 3, s\in[-\theta, 0]$ and $\phi_i(s)\in C([-\theta, 0], R^{+}), i = 1, 2, 3.$

Mathematically speaking, delays play a key role in determining how different differential systems behave dynamically ^[19,20,21]. Delays can lead to a variety of consequences, including changes in stability, the formation of bifurcations, and the start of chaotic behavior ^[22,23,24]. During the past decades, many authors deal with this topic, see ^[25,26,27]. Specifically, the dynamic phenomena of delays-induced Hopf bifurcation is quite significant ^[28,29,30]. From a biological perspective, delay-induced Hopf bifurcation nicely captures the equilibrium between densities of different living populations ^[31]. We contend that in order to understand the dynamics of interactions across various biological populations, it is critical to investigate delay-induced Hopf bifurcation in a variety of biological models. Motivated by this idea, we will concentrate on the Hopf bifurcation and its bifurcation control mechanisms. Specifically, we will address three critical questions: (1) Assess the distinctive features of the solution to system (1.3), including its nonnegativity, existence and uniqueness, and boundedness. (2) Examine the system's stability problem and the formation of the Hopf bifurcation phenomena (1.3). (3) Build both controllers for adjusting the range of stability and the generation time of the Hopf bifurcation in system (1.3).

The following is an introduction of the study's primary luminous spot: (ⅰ) A completely new delay-independent bifurcation and stability criterion for system (1.3) is established based on the earlier studies. (ⅱ) Using separate controllers, the realm of stability and the time of Hopf bifurcation in system (1.3) may be regulated well. (ⅲ) This study examines how delay affects Hopf bifurcation and the stability of first and second species densities in system (1.3).

The stated arrangement of the article is as follows. The "Feature of solution" section discusses the unique characteristics of system (1.3), such as boundedness, nonnegativeness, existence, and uniqueness. The "Bifurcation analysis" section discusses the bifurcation phenomena and the system's stability (1.3). The "Bifurcation domination using extended hybrid controller Ⅰ" section centers around the control issue associated with bifurcation phenomena for system (1.3) by formulating a plausible hybrid delayed feedback controller that involves parameter perturbation is accompanied by latency and state feedback. "Bifurcation domination using extended hybrid controller Ⅱ" section addresses the control challenge of the bifurcation phenomena in system (1.3). Creating a suitable hybrid delayed feedback controller with parameter perturbation for delay and state feedback. The "Software experiments" section provides Matlab software simulation results to assess the veracity of the obtained main results. "Conclusions" concludes this study concisely.

2.   Feature of solution

This section will examine the characteristics of the solution to system (1.3), such as its existence and uniqueness, nonnegativeness, and boundedness. To do this, we will employ fixed point theory, inequality methods, and the design of an appropriate function.

Theorem 2.1. Let $\Delta = \{(u_1, u_2, u_3)\in R^3: \max\{|u_1|, |u_2|, |u_3|\}\leq M\}$, where $M > 0$ denotes a constant. For every $(u_{10}, u_{20}, u_{30})\in \Delta$, system (1.3) under the initial value $(u_{10}, u_{20}, u_{30})$ owns a unique solution $U = (u_1, u_2, u_3)\in \Delta.$

Proof. Set

where

For arbitrary $U, \tilde{U}\in\Delta$, one gains

where

Let

Then it follows from (2.3) that

Consequently, $f (U)$ adheres to the Lipschitz condition for $U$. By applying the fixed point theorem, it's possible to infer with ease that Theorem 2.1 is right.

Theorem 2.2. Each solution to system (1.3) beginning with $R_{+}^3$ maintains a nonnegative value.

Proof. Considering the initial equation of system (1.3), it's possible to achieve

then

which leads to

then one gets

thus

In the same way, we know

Therefore, the validity of Theorem 2.2 is affirmed.

Theorem 2.3. Each solution to system (1.3), starting with $R_{+}^3$, is invariably restricted.

Proof. From the second equation of system (1.3), we have

Then

Integrating from $(t - \theta)$ to $t$ on both sides of the Eq (2.15) leads to

We get

then

from (2.14), we can get

So

Hence

Now, we define the following function

Then

Thus

Next we define another function $W_2$ as follows

Then, we get

where

Stemming from Eq (2.26), it is deduced that

By Eq (2.28), we obtain

Consequently, every solution to the system (1.3) is consistently limited.

3.   Bifurcation analysis

This part delves into examining the bifurcation and stability aspects of model (1.3). Initially, our assumption is that $E (u_{1\star}, u_{2\star}, u_{3\star})$ represents the balance point in model (1.3); subsequently, $u_{1\star}, u_{2\star}, u_{3\star}$ adhere to this stipulation:

Let

Substituting system (3.2) into system (1.3), we achieve the linear configuration of model (1.3) at $E (u_{1\star}, u_{2\star}, u_{3\star})$ (denote ${\bar u}_1$ as $u_1$, ${\bar u}_2$ as $u_2$, ${\bar u}_3$ as $u_3$)

where

The characteristic Eq (3.3) of the system owns the following expressions:

which leads to

where

If $\theta = 0,$ then Eq (3.6) becomes

If

is fulfilled, then the three roots $\lambda_1, \lambda_2, \lambda_3$ of Eq (3.8) have negative real parts. Thus, the equilibrium point $E(u_{1\star}, u_{2\star}, u_{3\star})$ of system (1.3) with $\theta = 0$ is locally asymptotically stable.

Assume that $\lambda = i\varepsilon$ is the root of Eq (3.6), then Eq (3.6) becomes

It follows from (3.9) that

where

and

It follows from (3.10) that

Because of ${\sin ^2}(\varepsilon \theta) + {\cos ^2}(\varepsilon \theta) = 1$, we can get

So,

Using (3.11) and (3.15), we know

Therefore, the results can be obtained as follows:

where

Let

Assume that

By virtue of $({\mathcal{A}}_2)$, we know ${\Re _1}(0) = - (c_6^2g_3^2 - g_2^2g_2^2) < 0$, and since $\mathop {\lim }\limits_{\varepsilon \to \infty } {\Re _1}(\varepsilon) > 0$, then we will know Eq (3.17) has at least one positive real root. Therefore Eq (3.6) has at least one pair of pure roots. Without loss of generality, we can assume that Eq (3.17) has twelve positive real roots (say, $\varepsilon _j, j = 1, 2, 3, \cdots, 12$). Relying on (3.12), we know

where $j = 1, 2, \cdots, 12; n = 0, 1, 2, \cdots$;

Assume ${\theta _0} = {\min _{\{ j = 1, 2, \cdots, 12;n = 0, 1, 2, \cdots \} }}\{ \theta _j^{(n)}\}$ and suppose that when $\theta = {\theta _0}$, Eq (3.6) has a pair of imaginary roots $\pm i{\varepsilon _0}$. Next we present the following assumption:

where

Lemma 3.1. Suppose that $\lambda(\theta) = \epsilon_1(\theta)+i\epsilon_2(\theta)$ is the root of Eq (3.6) at $\theta = \theta_0$ such that $\epsilon_1(\theta_0) = 0$, $\epsilon_2(\theta_0) = \varepsilon _0$, then $Re\left(\frac{d\lambda}{d\theta}\right)\Big|_{\theta = \theta_0, \varepsilon = \varepsilon_0} > 0.$

Proof. By Eq (3.6), we can get

It means that

where

Hence

By the assumption $({\mathcal{A}}_3)$, we get

This concludes the proof. According to the preceding reasoning, the following result is simply deduced.

Theorem 3.1. Assume that $({\mathcal{A}}_1)$–$({\mathcal{A}}_3)$ hold, then the equilibrium point $E(u_{1\star}, u_{2\star}, u_{3\star})$ of model (1.3) holds a locally asymptotically stable state if $\theta \in [0, \theta_{0})$ and model (1.3) generates Hopf bifurcations around the equilibrium point $E(u_{1\star}, u_{2\star}, u_{3\star})$ when $\theta = \theta_{0}.$

4.   Bifurcation domination using extended hybrid controller Ⅰ

Motivated by the works of ^{[30,32,33,34]}, we design the hybrid controller in Sections 4 and 5. In fact, we add a perturbation by adjusting the rate of change of population and then we can check the controllability via theoretical analysis and computer simulations. In this part, we will look at the Hopf bifurcation control issue in system (1.3) using a suitable hybrid controller that combines state feedback and parameter perturbation with delay. Using the ideas from ^[30,32], we achieve the following controlled 3D Lotka-Volterra commensal symbiosis system:

where $\delta _1, \delta _2$ stands for feedback gain parameters. Systems (4.1) and (1.3) own the same equilibrium points $E(u_{1\star}, u_{2\star}, u_{3\star})$. Let

denote ${\bar u}_1$ as $u_1$, ${\bar u}_2$ as $u_2$, ${\bar u}_3$ as $u_3$, and the linear system of system (4.1) around $E(u_{1\star}, u_{2\star}, u_{3\star})$ takes the following expression:

where

The following expression is owned by the system $(4.4)^{, }$s characteristic equation:

which leads to

where

If $\theta = 0,$ then Eq (4.6) reads as:

If

holds, then the three roots $\lambda_1, \lambda_2, \lambda_3$ of Eq (4.8) have negative real portions. Thus, the equilibrium point $E(u_{1\star}, u_{2\star}, u_{3\star})$ of model (4.1) remains locally asymptotically stable at $\theta = 0$. From (4.6), we can get

Suppose that $\lambda = i\zeta$ is the root of Eq (4.9). Then, Eq (4.9) takes

which results in

It follows from (4.11) that

where

and

It follows from Cramers rule that

In view of $\cos^2(\zeta\theta)+\sin^2(\zeta\theta) = 1$, we get

Then,

Using (4.13) and (4.17), we know

where

Let

Suppose that

holds, noticing that $\lim_{\zeta\rightarrow +\infty}{\Re _2}(\zeta) = +\infty > 0$, then we find that Eq (4.18) owns at least one positive real root. Thus, Eq (4.6) owns at least one pair of pure roots. Without loss of generality, here we assume that Eq (4.18) admits twelve positive real roots (say, $\zeta_\iota, \iota = 1, 2, 3, \cdots, 12)$. According to (4.15), one gets

where $\iota = 1, 2, 3, \cdots\, 12; k = 0, 1, 2\cdots;$

Denote $\theta_{\star} = \min_{\{\iota = 1, 2, 3, \cdots, 12; k = 0, 1, 2, \cdots\}}\{\theta_{\iota}^{(k)}\}$ and suppose that when $\theta = \theta_\star$, (4.6) owns a pair of imaginary roots $\pm i\zeta_{0}$.

Now, the following condition is presented:

where

Lemma 4.1. Let $\lambda (\theta) = {\eta _1}(\theta) + i{\eta _2}(\theta)$ be the root of Eq. (4.9) at $\theta = \theta_\star$ obeying $\eta_1(\theta_\star) = 0, \eta_2(\theta_\star) = \zeta_0$, then $Re\left(\frac{d\lambda}{d\theta}\right)\Big|_{\theta = \theta_\star, \zeta = \zeta_0} > 0.$

which leads to

where

Hence,

By $({\mathcal{A}}_6)$, one gets

which completes the proof.

Based on the study above, the following conclusion is lightly acquired.

Theorem 4.1. Suppose that $({\mathcal{A}}_4)$–$({\mathcal{A}}_6)$ hold, then the equilibrium point $E(u_{1\star}, u_{2\star}, u_{3\star})$ of model (4.1) holds locally asymptotically stable if $\theta\in [0, \theta_\star)$ and model (4.1) produces Hopf bifurcations at the equilibrium point $E(u_{1\star}, u_{2\star}, u_{3\star})$ when $\theta = \theta_\star.$

5.   Bifurcation domination using extended hybrid controller Ⅱ

In this part, we will investigate the Hopf bifurcation problem of system (1.3) using a suitable extended delayed feedback controller consisting of parameter perturbation with delay. In accordance with ^[33,34], we propose the following controlled predator-prey model:

where ${\rho _1}-{\rho _4}$ stands for feedback gain parameters. System (5.1) owns the same equilibrium point $E(u_{1\star}, u_{2\star}, u_{3\star})$ as that of system (1.3). Let

denote ${\bar u}_1$ as $u_1$, ${\bar u}_2$ as $u_2$, ${\bar u}_3$ as $u_3$, and the linear system of system (5.1) around $E(u_{1\star}, u_{2\star}, u_{3\star})$ takes the following expression:

where

The characteristic equation of system (5.3) owns the following expression:

which leads to

that is,

where

If $\delta = 0,$ then Eq (5.6) reads as:

If

is fulfilled, the three roots of Eq (5.6), \(\lambda_1, \lambda_2, \lambda_3\), have negative real components. Therefore, the equilibrium point $E(u_{1\star}, u_{2\star}, u_{3\star})$ of system (5.1) with $\theta = 0$ is locally asymptotically stable.

Suppose that $\lambda = i\omega$ is the root of Eq (5.7), then Eq (5.7) becomes:

By (5.10), we have

and from (5.11), we can get

where

So, we can get

where

From (5.14), we can suppose that $\cos (\omega \theta) = y$, and we have

According to the computer software, we can get

It follows from (5.17) that

where $t = 1, 2, 3, 4$; $n = 0, 1, 2, 3, \cdots$.

Let $\theta^\star = \min_{\{t = 1, 2, 3, 4; n = 1, 2, 3, \cdots\}}\{\theta_{t}^{(n)}\}$, and assume that when $\theta = \theta^\star$, Eq (5.6) has at least one pair of pure real roots $\pm i\omega_0$. Next, the following assumption is needed:

where

Lemma 5.1. Suppose that $\lambda(\theta) = \psi_1(\theta)+i\psi_2(\theta)$ is the root of Eq (5.7) at $\theta = \theta ^\star$ such that $\psi_1(\theta ^\star) = 0$, $\psi_2(\theta^\star) = \omega_0$, then $Re\left(\frac{d\lambda}{d\theta}\right)\Big|_{\theta = \theta ^\star, \omega = \omega_0} > 0.$

Proof. By Eq (5.7), one gets

which implies

where

Hence

By $({\mathcal{A}}_{8})$, we have

which concludes the proof.

Based on the research presented above, the following conclusion is loosely drawn.

Theorem 5.1. Suppose that $({\mathcal{A}}_{7})-({\mathcal{A}}_{8})$ hold, then the equilibrium point $E(u_{1\star}, u_{2\star}, u_{3\star})$ of model (5.1) is locally asymptotically stable if $\theta\in [0, \theta^\star)$ and model (5.1) generates Hopf bifurcations at the equilibrium point $E(u_{1\star}, u_{2\star}, u_{3\star})$ when $\theta = \theta^\star.$

Remark 5.1. In Section 4, the controller is called the hybrid controller that includes state feedback and parameter perturbation with delay. This controller is only added to the first equation of system (1.3). In Section 5, the controller is called the extended hybrid controller that includes state feedback and parameter perturbation with delay. This controller is added to the first equation and the third equation of system (1.3). Hybrid controller Ⅱ owns more control parameters than those in hybrid controller Ⅰ and has greater adjustment flexibility in controlling the stability domain and the onset of Hopf bifurcation of system (1.3).

6.   Software experiments

Example 6.1. Think about the following Lotka-Volterra commensal symbiosis system:

It is straightforward to see that system (6.1) has a single positive equilibrium point $E(0.7092, 0.5185, 0.1224)$. One can easily verify that the conditions $({\mathcal{A}}_1)$–$({\mathcal{A}}_3)$ of Theorem 3.1 hold true. Using Matlab software, we can obtain $\theta_{0}\approx1.3$. To validate the accuracy of Theorem 3.1, we use two distinct delay values: $\theta = 0.8$ and $\theta = 1.9$. For $\theta = 0.8 < \theta_{0}\approx1.3$, simulation graphs are provided in Figure 1. Figure 1 shows that $u_1\rightarrow0.7092, u_2\rightarrow0.5185, u_3\rightarrow0.1224$ as $t\rightarrow +\infty$. In this case, the equilibrium point $E(0.7092, 0.5185, 0.1224)$ of model (6.1) has a locally asymptotically stable state. For $\theta = 1.9 > \theta_{0}\approx1.3$, we receive simulation graphs, as shown in Figure 2. Figure 2 shows that $u_1$ maintains a periodic vibrating level around $0.7092$, whereas $u_2$ maintains a level around $0.5185$ and $u_3$ maintains a periodic vibrating level around $0.1224$. That is, a set of periodic solutions (known as Hopf bifurcations) arise at the equilibrium point $E(0.7092, 0.5185, 0.1224)$.

Example 6.2. Think about the following controlled Lotka-Volterra commensal symbiosis system:

It is straightforward to see that system (6.2) has a single positive equilibrium point $E(0.7092, 0.5185, 0.1224)$. Let ${\delta _1} = 0.5, {\delta _2} = -0.1$. One can easily verify that the conditions $({\mathcal{A}}_4)$–$({\mathcal{A}}_6)$ of Theorem 4.1 hold true. Using Matlab software, we can obtain $\theta_{\star}\approx1.8$. To validate the correctness of the acquired assertions of Theorem 4.1, we use two distinct delay values: $\theta = 1.6$ and $\theta = 2.25.$ For $\theta = 1.6 < \theta_{\star}\approx1.8$, we get simulation diagrams which are presented in Figure 3. Based on Figure 3, we find that $u_1\rightarrow0.7092, u_2\rightarrow0.5185, u_3\rightarrow0.1224$ when $t\rightarrow +\infty.$ In other words, the equilibrium point $E(0.7092, 0.5185, 0.1224)$ of model (6.2) holds a locally asymptotically stable state. For $\theta = 2.25 > \theta_{\star}\approx1.8$, we get simulation diagrams which are presented in Figure 4. Based on Figure 4, we find that $u_1$ maintains a periodic vibrating level around $0.7092$, whereas $u_2$ maintains a level around $0.5185$ and $u_3$ maintains a periodic vibrating level around $0.1224$. That is to say, a family of periodic solutions (namely, Hopf bifurcations) appear near the equilibrium point $E(0.7092, 0.5185, 0.1224)$.

Example 6.3. Think about the following controlled Lotka-Volterra commensal symbiosis system:

It is straightforward to see that system (6.3) has a single positive equilibrium point $E(0.7092, 0.5185, 0.1224)$. Let ${\rho _1} = 0.5, {\rho _2} = -0.1, {\rho _3} = 0.6, {\rho _4} = -0.1$. One can easily verify that the conditions $({\mathcal{A}}_7)$ and $({\mathcal{A}}_8)$ of Theorem 5.1 hold true. By applying Matlab software, one can get $\theta^\star\approx2.20$. To validate the correctness of the acquired assertions of Theorem 5.1, we choose both different delay values: $\theta = 2.00$ and $\theta = 2.75.$ For $\theta = 2.00 < \theta^\star\approx2.20$, we get simulation diagrams which are presented in Figure 5. Based on Figure 5, we find that $u_1\rightarrow0.7092, u_2\rightarrow0.5185, u_3\rightarrow0.1224$ when $t\rightarrow +\infty.$ In other words, the equilibrium point $E(0.7092, 0.5185, 0.1224)$ of model (6.3) holds a locally asymptotically stable state. For $\theta = 2.75 > \theta^\star\approx2.20$, we get simulation diagrams which are presented in Figure 6. Based on Figure 6, we find that $u_1$ maintains a periodic vibrating level around $0.7092$, whereas $u_2$ maintains a level around $0.5185$ and $u_3$ maintains a periodic vibrating level around $0.1224$. That is to say, a family of periodic solutions (namely, Hopf bifurcations) appear near the equilibrium point $E(0.7092, 0.5185, 0.1224)$.

Remark 6.1. It follows from the Matlab simulation results of Examples 6.1–6.3 that one can know that the bifurcation value of system (6.1) is equal to $1.3$, the bifurcation value of system (6.2) is equal to $1.8$ and the bifurcation value of system (6.3) is equal to $2.20$, which indicates that we can expand the domain of stability of system (6.1) and postpone the time of emergence of Hopf bifurcation of system (6.1) via the formulated two hybrid delayed feedback controllers.

7.   Conclusions

It is generally recognized that the delayed dynamical model is an important tool for understanding the interactions of many biological populations in the natural environment ^[35,36,37]. Many studies on predator-prey models have conducted and yielded numerous results over the last few decades ^[38,39,40]. In this study, we provide a novel delayed Lotka-Volterra commensal symbiosis model. This paper discusses the uniqueness, nonnegativeness, and boundedness of the delayed Lotka-Volterra commensal symbiosis solution. The Hopf bifurcation issue is addressed. Then, the critical delay value $\theta_0$ is retrieved. In order to modify the domain of stability and the time of the bifurcation phenomenon in this model, we have successfully developed two distinct hybrid delayed feedback controllers. Two critical delay values, $\theta_{\star}, \theta^\star,$ are acquired. In these two controllers, the role of delay is displayed. Theoretically, the exploration fruits are very useful for managing and balancing the populations of two species. Furthermore, the exploratory concepts may be used for other fractional-order and integer-order dynamical systems in a wide range of disciplines to dominate the bifurcation phenomena, stability, and chaos ^[41,42,43]. During the past decades, many works on this topis is explored, see ^[44,45,46]. In this paper, we only deal with the Hopf bifurcation onset and Hopf bifurcation control in this paper. We leave the stability and direction of Hopf bifurcation periodic solutions for future work and we will refer to the works in ^[47,48,49]. In addition, we will explore the Hopf bifurcation of fractional-order dynamical models ^[50,51,52].

Author contributions

Yingyan Zhao: Conceptualization, formal analysis, investigation, methodology, software, writing-original draft, writing-review & editing; Changjin Xu: Conceptualization, formal analysis, investigation, methodology, software, writing-original draft, writing-review & editing; Yiya Xu: Formal analysis, investigation, software, writing-review & editing; Jinting Lin: Conceptualization, formal analysis, investigation, methodology, writing-original draft, writing-review & editing; Yicheng Pang: Conceptualization, investigation, methodology, writing-review & editing; Zixin Liu: Conceptualization, investigation, methodology, writing-review & editing; Jianwei Shen: Conceptualization, investigation, methodology, writing-review & editing. All authors have read and approved the final version of the manuscript for publication.

Acknowledgments

This work is supported by National Natural Science Foundation of China (No. 12261015, No. 62062018), Project of High-level Innovative Talents of Guizhou Province ([2016]5651), Guizhou Provincial Science and Technology Projects (ZK[2022]025) and Science Project of Guizhou University of Finance and Economics (2022XSXMB08).

The authors would like to thank the referees and the editor for helpful suggestions incorporated into this paper.

Conflict of interest

The authors declare that they have no conflict of interest.