The dynamics of a Leslie type predator–prey model with fear and Allee effect

In this article, we discuss the dynamics of a Leslie–Gower ratio-dependent predator–prey model incorporating fear in the prey population. Moreover, the Allee effect in the predator growth is added into account from both biological and mathematical points of view. We explore the influence of the Allee and fear effect on the existence of all positive equilibria. Furthermore, the local stability properties and possible bifurcation behaviors of the proposed system about positive equilibria are discussed with the help of trace and determinant values of the Jacobian matrix. With the help of Sotomayor’s theorem, the conditions for existence of saddle-node bifurcation are derived. Also, we show that the proposed system admits limit cycle dynamics, and its stability is discussed with the value of first Lyapunov coefficient. Moreover, the numerical simulations including phase portrait, one- and two-parameter bifurcation diagrams are performed to validate our important findings.


Introduction
One of the prominent topics for ecologists, biologists, and even mathematicians is the study of various mechanisms related to the interaction between prey and predators. The typical model of such interactions between two species was first proposed in the form of a simplest mathematical model [1,2]. The key factor in a predator-prey relationship is that the growth of predator depends on its functional response to the prey population, which describes the amount of prey each predator eats per unit of time. Many researchers have developed and implemented various forms of interaction terms to understand the more realistic situation. Thus the predation function depends on prey only such as Holling type classification [3,4], whereas ratio-dependent [5], Beddington [6], and Crowley Martin [7] are the functions of both prey and predator. In real-world situations, [8,9] clearly confirm that the ratio-dependent interaction has a greater description of predation rates than the prey-dependent interaction. The authors in [10] assumed that the predator growth function is different from the predator's predation function, and its reduction has a reciprocal relationship with per capita availability of its preferred food, the so-called Leslie type predator-prey model. Despite its rich and varying complexity, there has been greater focus paid by several researchers to modify and build the Leslie type model in several ways [7,11,12]. In some real situations, fractional-order systems are found to be more accurate than integer-order ones. For this reason, a fractional-order predator-prey system has been investigated by several researchers, and numerous outstanding results have been derived, for example, see [13,14].
The Allee effect is one of the most interesting phenomena in ecology, which is characterized by a positive correlation between the population size and the per capita growth rate for a population [15]. The Allee effect can be caused by various environmental factors including difficulties in finding mating partners at low density, genetic inbreeding, social felicitation of reproduction, low probability of successful mating, depletion in inbreeding rate, etc. Recently, there has been greater interest from many researchers to study the effect of incorporating the Allee effect in the predator-prey dynamics. In the ratiodependent predator-prey model, the authors in [16] compared the dynamical properties with and without the Allee effect in logistic prey growth term and showed that the inclusion of Allee effect can remove oscillation behavior of species densities. Further, the conditions for occurrence of various bifurcation behaviors were derived and verified numerically with the help of parametric diagrams and phase portraits. The same model with discrete time was investigated in [17], and the authors showed that varying integral step size can arise chaotic orbits via period-doubling and invariant circle. A Leslie type model with ratio-dependent type interaction term with state impulsive feedback control was considered, and the conditions for existence, uniqueness, and orbital asymptotical stability of periodic solution of order-1 were obtained by applying the geometry theory of differential equations [18]. The authors [19] considered a strong Allee effect in the prey growth function in the Holling type II model with density-dependent death rate for predators and examined various types of bifurcation behavior in the considered model. Since there have been fewer studies related to the double Allee effect in the predator-prey system, in [20] the authors attended the presence of double Allee effect in the ratio-dependent predatorprey system in the prey growth function, where they studied various types of bifurcation behavior and showed bi-stability behavior in the phase plane implying that the considered system dynamics is sensitive in choosing the initial conditions. Predators will actively affect the ecosystem and the prey reproduction process by consuming them and by affecting the behavior and physiology of the survivor indirectly [21]. Compared to direct predation, these forms of behavioral and physiological modifications related to fear of predation are more efficient and longer-lasting evolutionary effects [21,22]. For starters, birds respond with anti-predator defences to the sounds of predators and flee from their nests at the first sign of danger [22]. Fear of song sparrows during the breeding season by eliminating direct predation, recognizing that the understanding of predation risk alone reduces the amount of offspring produced by 40 percent per year [23]. Recently, many mathematical models have been studied with effect of fear in the prey population due to anti-predator defense, see for details [24][25][26].
The authors in [27] showed that the predator-prey model with Beddington-DeAngelis type interaction has both stabilizing and destabilizing effect on incorporating the fear term in the prey. In [28] the effect of fear of predator in the prey-predator model was investigated and the existence, stability of equilibria, and occurrence of limit cycle in the considered system were studied. The authors in [29] studied the dynamics of cost of fear and the Allee effect in the predator-prey model, also with disease in the prey population, while the dynamics of a prey-predator system with fear and group defense was considered in [30]. The predator populations in the natural world are more vulnerable to Allee effects than their prey, and because they are bigger, populations are more likely to undergo an Allee effect (predator populations are usually much smaller than prey populations). The model with Allee effect in the predator population was considered in [31,32]. Since there have been various dynamics and bifurcation behavior arising in the nonlinear population model, in [33] the authors attempted to study the predator-prey system with nonlinear harvesting and showed that the considered system has unstable periodic orbit by finding the value of first Lyapunov coefficient. They also showed the existence of saddle-node and transcritical bifurcation by Sotomayor's theorem, and Bogdanov-Takens bifurcation by finding universal unfolding near cusp. To understand the dynamics of fear and Allee effect, more works are required, and this fact has motivated our current study.
Encouraged by the above discussion, we aim to investigate the combined effect of incorporating both fear and Allee effect in a Leslie type predator-prey model with ratiodependent type interaction. To the best of our knowledge, this assumption has not been studied. We are interested in a way the system exhibits various complex dynamical behaviors. The main highlights of this article are summarized as follows: • The predator-prey model with fear effect in prey growth term and Allee effect in predator growth term is considered and analyzed. • The proposed models utilize the fact that predators have to go for alternative food when there is low density of their favorite food. • To confirm the long term coexistence of species of the proposed model, permanence analysis is performed. • The impact of introduced parameters on system dynamics is investigated through different types of bifurcations. • By taking fear and Allee effect constant as a bifurcation parameter, the proposed theoretical results are verified through numerical simulations. This article is organized as follows: In Sect. 2, we describe the formation of the system along with its assumptions, and some preliminaries of the system are given in Sect. 3. In Sect. 4, we address the existence of the system's equilibria and local stability with numerical examples. In Sect. 5, we explore the complex dynamical behavior of our system by varying the Allee parameter with a numerical example. We give numerical simulations in Sect. 6, and the conclusion is given in the last section.

The mathematical model
The mathematical model of the modified Leslie-Gower predator-prey model as in [11,12] is given by where x(t), y(t) are the prey and predator biomass at time t, f 1 (x, K 1 ) is a function representing the gross rate of increase in the prey population with carrying capacity K 1 in the lack of predators and satisfying f 1 (0, K 1 ) > 0, f 1 (K 1 , K 1 ) = 0, ∂f 1 ∂x ≤ 0, ∂f 1 ∂t | x=K 1 < 0, and ∂f 1 ∂K 1 > 0. And nx in (2.1b) is the prey-dependent carrying capacity of the predator, and it is assumed that c is a constant taken for the alternative food source, which in turns means that the predator proportion is small compared to the alternative food source. We consider the logistic growth of the form f 1 (x, K 1 ) = r 0ebx, with K 1 = r 0 -e b . Also, as in [5], the ratio-dependent (ratio between prey and predator) in Michaelis-Menten type interaction is implemented, which is of the form g( x y ) = px x+ay , where p, a are predators searching level and handling time spent for each prey. Then (2.1a) takes the form where r 0 , e, and b are positive constants, representing the birth rate, death rate, and death rate caused by intra-species competition. The dynamical analysis of a Leslie-Gower model with ratio-dependent interaction is performed in [34]. Since the prey production is indirectly influenced by the fear effect [35], we also attempted to include the fear effect in the prey growth term in (2.2). Hence, we multiply the function f 2 (k, y) = 1 1+ky with the birth rate r 0 , which is the cost of anti-predator defense due to fear [26], and for other fruitful results see [27][28][29]. Here, the quantity k measures the fear level on prey, and by the biological meanings of k, y, and f 2 (x, y), it is appropriate to take that Then (2.2) with fear term is of the following form: Since there have been fewer studies related to the presence of Allee effect on the predator population, in [12] the authors attempted a Leslie-Gower model with Holling type interaction model and studied the dynamics with effect to the Allee effect in both prey and predator, respectively. Hence, it is reasonable and more realistic to take the Allee effect into account for the predator population. We arrive at the following model with fear in the prey population and the Allee effect in the predator population: where y y+m is the term for the Allee effect and m is the Allee effect constant. The per capita growth rate of the predator is reduced from s to sy y+m . The bigger the m, the stronger the Allee effect, the slower the per capita growth rate of the predator population, and we always assume H1 r 0 > e, i.e., the reproduction rate of prey is always greater than its death rate.

Positivity
Biologically, the survival of the population is positivity. We follow Theorem A.4 in [36] to prove the positivity of system (2.4a)-(2.4b). We define the state space of (2.4a)-(2.4b) as X = {(x, y) ∈ R 2 + } due to the biological meaning. Now the functions g 1 (·, ·), g 2 (·, ·) are continuously differentiable and locally Lipschitz in X. Then, by Theorem A.4 in [36], we conclude that, for a sufficiently large number S, the solution of the initial value problem (2.4a)-(2.4b) exists and is unique in the interval [0, S).

Boundedness
Boundedness, whether there is a well-specified beginning or end of a case. As a consequence of limited resources, boundaries can be viewed as a natural constraint on growth.
Proposition 3.1 Assume that r 0 > e and r 0 -e b + n + 1 > qm, where 1 > 0 is an arbitrary small number. Then all the solutions (x(t), y(t)) of system (2.4a)-(2.4b) with the initial con- Since r 0 > e by H1, applying Lemma 3.1, we obtain Thus, there existsT 1 ∈ R,T 1 > 0 for some arbitrary constant 1 > 0 such that Also, (2.4b) can be written as Thus, for all t >T 1 , we have Thus,T 2 >T 1 ,T 2 ∈ R for some arbitrary constant 2 > 0 such that Therefore, the solutions of system (2.4a)-(2.4b) are bounded.

Permanence
The system's permanence is important as it ensures that the species present now will not be extinct in the future. System (2.4a)-(2.4b) is said to be permanent if there exist positive constants M 3 and M 4 such that each positive solution ( where M 2 is the upper bound of y population and w 0 = r 0 1+kM 2ep. Assume that w 0 > 0. Then, by applying Lemma 3.1, we obtain Also, from (2.4b), we have Assume that cm > 0, again by using Lemma 3.1, we obtain Then the condition of permanence of system (2.4a)-(2.4b) follows.

Equilibria and local stability analysis 4.1 Existence of equilibria
In order to obtain the equilibrium points of system (2.4a)-(2.4b), we consider the following equations: Simple computation yields the following equilibrium points of system (2.4a)-(2.4b): • Trivial equilibrium E 0 (0, 0) exists by [38].
, and y * is the positive root of the cubic equation where A 1 = bk(an + 2)(mc) + (an + 1)(b + ekn) + kn 2 p 3bk(an + 1) , By introducing a transformation z = y * + A 1 , (4.2) can be reduced to By [39], the existence of the positive roots of (4.3) can be stated in the following lemma. Moreover, algebraic computations show that (4.3) has two positive roots z 1 = It is noted that if (4.3) has one positive root, it must be the same as z 1 . Suppose p 2 > 0, p 1 < 0, and using the Descartes rule of sign change, we discuss the existence of possible interior equilibria of system (2.4a)-(2.4b) in the following cases.
(2) A unique positive interior equilibrium point of multiplicity two, i.e.,  (4.2) has a unique positive root, then the interior equilibrium is Note that the Allee effect term m plays an important role for the existence of interior equilibrium points. Moreover, if m < c, the two equilibria E 2 and E 3 collide to get E 1 and disappear if m crosses some particular threshold value, it states possible saddle-node bifurcation around E 1 . Furthermore, if m = c, then one of the equilibria E 2 , E 3 will disappear and one will collide with E 0 another with E 4 , then there may be possible existence of transcritical bifurcation around E 0 .
Hence, the expression for the critical threshold value m for the above scenario of the existence of positive equilibria is difficult to find, so we discuss it with the help of numerical example and nullcline analysis.  Fig. 1(a)) and, if m = 0.08586597105258, we get the unique positive interior equilibrium E 1 (0.07208, 0.18621) withÊ(0, 0.114134), and for m = 0.12, we get two positive interior equilibria E 2 (0.156899, 0.236899), E 3 (0.018512, 0.098512) withÊ(0, 0.08). If m = c = 0.2, we have a unique positive interior equilibrium E 4 (0.25, 0.25), bothÊ and E 3 collide with E 0 , and for m = 0.3 > c, E 3 disappears and we have E 5 (0.341366, 0.241366). Finally, in all cases we have E 0 (0, 0) andĒ(0.95, 0). This is shown in Fig. 1(b). Since the value of a fear function decreases as the level of fear increases as in [40], the parameter k also has the same capability for the possible existence of saddle-node bifurcation, see Fig. 1(c). The saddle-node bifurcation plot for varying m is given in Fig. 2(a) and for k in Fig. 2(b). The region for existence of interior equilibria in km plane is given in Fig. 2(c).
Thus, the trace(Tr) and determinant (Det) values of the Jacobian matrix play a prominent role in studying the local stability of system (2.4a)-(2.4b), we discuss the local stability in the following subsection.

Local stability analysis
The Jacobian matrix of the system at any equilibrium point (x, y) is given by Note that the Jacobian matrix cannot be evaluated at E 0 (0, 0). To study the stability of system (2.4a)-(2.4b) at E 0 , we transform the ratio-dependent predator-prey system to an equivalent system by making use of transformation x = x, y = y x , then the resultant system takes the form: Thus, the above system has the same equilibrium E 0 = (0, 0) and axial equilibriumĒ(x, 0) = ( r 0 -e b , 0) after re-transfer to its original variables. The Jacobian matrix of system (2.4a)-(2.4b) at E 0 from (4.5) is given as The Jacobian matrix of system (2.4a)-(2.4b) at the axial equilibria pointsĒ andÊ is given by Since the sign of the eigenvalues of J E 0 , JĒ, JÊ determines the local stability of system (2.4a)-(2.4b) around the equilibrium points E 0 ,Ē,Ê, thus the conditions for local stability around the E 0 ,Ē,Ê are summarized in the following proposition. Let E * (x * , y * ) be an arbitrary interior equilibrium whose Jacobian matrix is given as follows: (4.7) Thus, the characteristic equation of the above matrix is where Thus, E * is locally asymptotically stable if it satisfies Routh-Hurwitz criteria that are Tr(J E * ) < 0 and Det(J E * ) > 0. Thus, we have the following proposition.

Bifurcation analysis
In this section, we focus on the possible existence of various bifurcation behaviors of system (2.4a)-(2.4b). In the nonlinear population interaction model, the most common phenomenon is the existence of periodic solutions around the equilibrium point, that is, the local birth or death of a periodic solution, the so-called Hopf bifurcation. Moreover, the two branches of stable and unstable equilibria created collide and disappear and are defined as saddle-node bifurcation. Finally, the exchange of stability properties of two equilibria is known as transcritical bifurcation. Hence, the bifurcation behavior of system (2.4a)-(2.4b) is discussed in what follows. Proof By taking m as the bifurcation parameter, we use Sotomayor's theorem as in [41] to prove that system (2.4a)-(2.4b) admits a saddle-node bifurcation. According to [41], since Det(J E 1 ) = λ 1 λ 2 = 0, either of the eigenvalues must be zero and another is negative, also Tr(J E 1 ) < 0. Let g = (g 1 , g 2 ) T , where g 1 , g 2 are given in Sect. 2. Then the matrix J E * at E 1 (x 1 , y 1 ) is written by
Example 5.1 Let us take the parameter values as in Example 4.1, when k = 1 the two equilibria E 2 , E 3 collide at some critical magnitude of m, i.e., m = m s = 0.08586597105258 and a unique equilibrium E 1 emerges and disappears for m < 0.08586597105258 = m s , which ensures the possible existence of saddle-node bifurcation near E 1 (see Fig. 2(a)). In addition, by [41] the values 1 = -0.0302379 = 0, 2 = 7.70607 = 0 confirm that system (2.4a)-(2.4b) undergoes saddle-node bifurcation near E 1 .

Hopf bifurcation
In this subsection, the possible occurrence of limit cycle dynamics (Hopf bifurcation) near the arbitrary interior equilibrium E * (x * , y * ) is to be analyzed by varying m and fixing all other parameters as constant.
Proof To prove that system (2.4a)-(2.4b) exhibits limit cycle dynamics, we need to consider that at some threshold value m = m h the trace value of J E * should be zero. The threshold value m h is calculated by Thus, the characteristic equation (4.8) in Sect. 4 becomes

Direction and stability
To study the direction of the possible occurrence of limit cycle oscillation and its stability nature originating around the arbitrary equilibrium E * , we derive the first Lyapunov coefficient (l 1 ) and use the results in [42].
Let us use the transformation Thus, the above system takes the form where S 1 and S 2 are nonlinear in u and v, given by By using the normal form theory in [42], we need the following quantity l 1 to study the properties of the limit cycle, which is calculated as follows: where the partial derivative values are calculated at the critical point, i.e., (v 1 , v 2 ; m) = (0, 0; m h ). By applying the result given in [42], the existing limit cycle (Hopf bifurcation) is supercritical and subcritical if l 1 < 0 and l 1 > 0 correspondingly. Tr(x 4 (m), y 4 (m))] m=m h = 0.268271 = 0 ensures that system (2.4a)-(2.4b) has limit cycle oscillation near E 4 at m h . Also, system (2.4a)-(2.4b) has θ 0 = 0.185849 and the first Lyapunov coefficient l 1 = -1.53282 < 0. Therefore, the existing Hopf bifurcation around E 4 is supercritical (unstable and surrounded by a stable limit cycle), see Figs. 3, 4, 5, and 6.

Bogdanov-Takens bifurcation
In the previous subsection, we discussed the existence of one-parameter bifurcation for system (2.4a)-(2.4b). In this subsection the possible existence of a codimension-2 bifurcation, namely Bogdanov-Takens (BT) bifurcation for (2.4a)-(2.4b), is analyzed. For this, we have to take two parameter values m and k as the varying parameters and fix all other parameter constants. Here, we only discuss it numerically. Remark 5.3 In most of the existing literature, authors have incorporated fear and Allee effect separately on predator-prey models, see [27,28,32,40,43]. For example, the issue on the impact of fear effect in the predator-prey model with Holling type II predator-prey interaction was reported in [40], where they stated that fear constant is able to stabilize the proposed system at an interior steady state, that is, both the species can exist together, or it can create the oscillatory coexistence of species populations. Weak Allee effect on the stability of a discrete-time predator-prey model was given in [43], and they observed that the system becomes stable from chaotic dynamics as the Allee parameter increases. The red line represents prey nullcline, the black line represents predator nullcline, the red and green arrow lines represent stable trajectories along x and y axes respectively, which approachÊ andĒ. The black dashed arrow line represents the separatrix curve separating the trajectories which approachĒ and interior equilibria E 2,4,5 in the x-y plane Authors in [44] discussed a food chain model that experienced strange dynamics by varying Allee effect parameters, that is, the system state changes from stable to chaotic via period doubling as Allee effect increases. The effects of both Allee and fear parameters in predator-prey interaction model with prey-dependent functional response were considered in [29]. It is worth mentioning that our derived ratio-dependent modified Leslie type model (2.4a)-(2.4b) incorporates fear effect in the prey population and Allee effect in the predators growth. Moreover, system (2.4a)-(2.4b) exhibits complex dynamics with respect to both fear and Allee parameters, which is explained theoretically through bifurcation analysis and is verified in terms of numerical simulations.

Numerical simulations
In this section, in order to obtain a better visualization of how different parameter selection values influence the system (2.4a)-(2.4b) dynamics, numerical simulations are shown. Since we aim to investigate the effect of fear(k) and the Allee effect(m), all other parame-  The scenario of no equilibrium, unique interior equilibrium, and two interior equilibria is given in Example 4.1 and showed in Fig. 1 with respect to the parameters m and k. The saddle-node bifurcation behavior for system (2.4a)-(2.4b) for parameter m is shown in Fig. 2(a) and for k in Fig. 2(b), where the blue dashed line represents the stable equilibrium branch and the red line represents the unstable equilibrium branch. The regions for existence of interior equilibria are shown in Fig. 2(c), where Region I is the space below the black darker line, here interior equilibrium point is not possible. Region II is the space between the black dashed line and the black thick line, here two interior equilibrium points arise. On and above the black dashed line representing Region III, the unique interior equilibrium point appears. In Fig. 3, the two-parameter bifurcation diagram is plotted, since the trace value is zero on the blue line, which separates the stable and unstable region, and the red dot denotes the Bogdanov-Takens bifurcation point, i.e., the intersection point of trace and determinant equals zero.
In order to show the birth of periodic solution of system (2.4a)-(2.4b), fixing k = 2.5 and varying m in (0, 0.2), the one-parameter bifurcation diagram is plotted in Fig. 4. When 0.16 < m < 0.174, the only equilibriumĒ exists and is stable. For 0.174 < m < 0.185, the periodic solution around the interior equilibrium E 2 exists and it is stable in 0.185 < m < 0.2. Hence, it states that the Allee parameter has a stabilizing effect on system (2.4a)- (2.4b). Also, in Fig. 5, the one-parameter bifurcation diagram is plotted with respect to the parameter k considering the cases in Lemma 4.1: For case (i) m < c = 0.18 in red, for case (ii) m = c = 0.2 in blue, and for case (iii) m > c = 0.22 in green, which shows the birth of periodic solution on increasing k. Therefore, the fear parameter k has a destabilizing role for system (2.4a)-(2.4b). The existence of stable, periodic, and unstable phase portraits around the equilibrium E 2 for different choice of m and k is shown in Fig. 6.
The time trajectories may tend to eitherĒ andÊ or E 2 , it depends on the choice of the initial condition. As in Fig. 6, there is a separatrix curve separating the trajectories that approachĒ = (0, 0.3) and E 2 = (0.137202, 0.167202), thus system (2.4a)-(2.4b) has three stable equilibria at m = 0.17, k = 2.14 and the remaining parameters are the same as in Example 4.1. Since the system is sensitive to the initial condition, we attempt to evaluate the 133 × 133 initial conditions in prey vs predator plane and mark in different colors in 7(c), red color approachesÊ, green color approachesĒ, and blue color approaches E 2 . In Figs. 7(a) and 7(b), it is shown that the different initial conditions approach different equilibria in time trajectories for both prey and predator, respectively.
The local amplification of two-parameter spaces k and m in Fig. 3 around the BT point (Bogdanov-Takens bifurcation point) is shown in Fig. 8(a) for clear view, which is (k, m) = Figure 10 Two-parameter bifurcation diagram for system (2.4a)-(2.4b) with different choice of parameter values, the blue curve represents Hopf bifurcation curve, the black curve represents determinant zero (below in (b), (c) and above in (d), system (2.4a)-(2.4b) has no interior equilibria), and the black dashed curve represents m = c (1.6083, 0.1175466), and the corresponding phase portrait around E 1 = (0.0585234, 0.140977) is shown Fig. 8(b), which is the cusp of co-dimension 2.
The values of trace and determinant are shown in the colored density graph for E 2 in Figs. 9(a), 9(b) and for E 3 in Figs. 9(c), 9(d). In Fig. 9(a) the black dashed line (Tr(E 2 ) = 0) states that the trace value changes from negative to positive, also in Fig. 9(b) the determinant value is always positive. Hence system (2.4a)-(2.4b) shows limit cycle oscillations. However, in Fig. 9(c), even the trace values change from negative to positive, the determinant values are always negative in Fig. 9(d). Therefore, system (2.4a)-(2.4b) is stable near E 2 and unstable near E 3 . Finally, in Fig. 10, the two-parameter bifurcation diagram is plotted for the equilibria E 2 , E 4 , and E 5 , which includes Hopf bifurcation curve, stable and unstable regions with different choices of system parameters. Note that if E 3 exists, it is always unstable.

Conclusion
In the present study, the fear effect in the prey population and the Allee effect in the predator population were considered in the modified Leslie type predator-prey model with ratio-dependent interaction term. First, we discuss the positivity, boundedness, and permanence analysis for system (2.4a)-(2.4b). Then the condition for existence of possible equilibria of system (2.4a)-(2.4b) is derived which showed that two equilibria arise, annihilate, and disappear with influence of both fear and Allee effects. Furthermore, we studied the local stability of all possible equilibria. For certain threshold values of Allee effect term, system (2.4a)-(2.4b) has two interior equilibria. Hence, we discussed the an-alytical conditions for local stability of arbitrary interior equilibria and plotted the stable and unstable regions using two-parameter bifurcation diagram. The colored map of trace and determinant values supports the study of stability of two possible interior equilibria. In addition, the considered system experiences the various bifurcation behaviors by varying a suitable parameter, and we observed that the increase of fear effect can change the system dynamics from stable to unstable via Hopf bifurcation. Also, we found that the existing Hopf bifurcation is supercritical by calculating the first Lyapunov coefficient. We showed that the Allee effect in predator helps to stabilize the unstable behavior caused by fear effect. Additionally, the existence of saddle-node bifurcation is discussed with the help of Sotomayor's theorem. It is important to analyze the effect of fear and Allee effects on other system parameters, for this we have plotted various two-parameter bifurcation diagrams with the choice of different parameter values. Hence, the fear and the Allee effects play a vital role in the system dynamical behavior, it is important to analyze the long term survival for both species. Thus, it could be interesting and meaningful to study the dynamics of a predator-prey model with the Allee effects in both prey and predator. However, these terms will increase the complexity of the system, and we will leave this as future research.