Bifurcation Scenario Research Articles

Rich dynamics of a hepatitis C virus infection model with logistic proliferation and time delays

AbstractIn this paper, we study a dynamic model of hepatitis C virus (HCV) infection with density‐dependent proliferation of uninfected and infected hepatocytes and two time delays, which is derived from a three‐dimensional model by the quasi‐steady‐state approximation. The model can exhibit forward bifurcation or backward bifurcation, and an explicit control threshold parameter is obtained for the case of backward bifurcation. It is shown that if the proliferation rate of infected hepatocytes is greater than the proliferation rate of uninfected hepatocytes by a certain amount, it becomes more difficult for the virus to be removed. The model has rich dynamical properties: (i) In some parameter regions, bistability can occur; (ii) both time delays (virus‐to‐cell delay) and (cell‐to‐cell delay) can lead to Hopf bifurcations; (iii) same length of time delays and can lead to at most one stability switch, but different time delays can lead to multiple stability switches. Several sufficient conditions for the global stability of the disease‐free equilibrium and the endemic equilibrium are obtained for both forward and backward bifurcation scenarios. In particular, several sharp results on global stability are obtained. Theoretical and numerical results portray the complexity of viral evolutionary dynamics in chronic HCV‐infected patients.

Seasonal refuge patterns of phytoplankton trigger irregular bloom events in a contaminated environment

Several experimental evidences and field data documented that zooplankton may alter its behavioral response in the presence of toxic phytoplankton, reducing its consumption to the point of starvation. This paper is devoted to the mathematical study of such interactions of toxic phytoplankton with grazer zooplankton. The non-toxic phytoplankton is assumed to adopt a density-dependent refuge strategy to avoid over-predation by zooplankton. Both groups of phytoplankton are assumed to suffer direct harm from anthropogenic toxicants, while zooplankton is affected indirectly by ingesting contaminated phytoplankton. We calibrate the proposed model with the field data from Talsari and Digha Mohana, India, and estimate some crucial model parameters consistent with the behavior of the observed data. Our results demonstrate that zooplankton grazing on toxic phytoplankton plays a key role in the emergence or mitigation of plankton blooms. We also highlight the system’s potential to exhibit multiple stable configurations under the same ecological conditions. The plankton system experiences significant regime shifts, which are explored through various bifurcation scenarios, such as transcritical and saddle-node bifurcations. These shifts are influenced by changes in refuge capacity, species growth rates, and environmental carrying capacity. Furthermore, we incorporate environmental variations due to seasonal periodic or almost periodic changes, allowing the refuge parameter to be time-dependent. We observe that the forced system exhibits double periodic solutions. Moreover, stronger seasonal variations in the refuge pattern lead to irregular chaotic blooms. In conclusion, the results offer valuable insights into the sustainability of biodiversity, potentially shedding light on the origin of diverse plankton bloom phenomena.

Stability and Bifurcation Analysis of a Symmetric Fractional-Order Epidemic Mathematical Model with Time Delay and Non-Monotonic Incidence Rates for Two Viral Strains

This study investigates a symmetric fractional-order epidemic model with time delays and non-monotonic incidence rates, considering two viral strains. By confirming the existence, uniqueness, and boundedness of the system’s solutions, the research ensures the model’s well-posedness, guaranteeing its mathematical soundness and practical relevance. The study calculates and evaluates the equilibrium points and the basic reproduction numbers R01 and R02 to understand the dynamic behavior of the model under different parameter settings. Through the application of the Lyapunov method, the research examines the asymptotic global stability of the system, determining whether it will converge to a particular equilibrium state over time. Furthermore, Hopf bifurcation theory is employed to investigate potential periodic solutions and bifurcation scenarios, highlighting how the system might shift from stability to periodic oscillations under certain conditions. By utilizing the Adams-Bashforth-Moulton numerical simulation method, the theoretical results are validated, reinforcing the conclusions and demonstrating the model’s applicability in real-world contexts. It emphasizes the importance of fractional-order models in addressing epidemiological issues related to time delays (τ), individual heterogeneity (m, k), and memory effects (θ), offering greater accuracy compared with traditional integer-order models. In summary, this research provides a theoretical foundation and practical insights, enhancing the understanding and management of epidemic dynamics through fractional-order epidemic models.

Influences of short-term and long-term plasticity of memristive synapse on firing activity of neuronal network

Abstract Synaptic plasticity can greatly affect the firing behavior of neural networks, and it specifically refers to changes in the strength, morphology, and function of synaptic connections. In this paper, a novel memristor model, which can be configured as a volatile and nonvolatile memristor by adjusting its internal parameter, is proposed to mimic the short-term and long-term synaptic plasticity. Then, a bi-neuron network model, with the proposed memristor serving as a coupling synapse and the external electromagnetic radiation being emulated by the flux-controlled memristors, is established to elucidate the effects of short-term and long-term synaptic plasticity on firing activity of the neuron network. The resultant 7D neuron network has no equilibrium point and its hidden dynamical behavior is revealed by phase diagram, time series, bifurcation diagram, Lyapunov exponent spectrum and two-dimensional dynamic map. Our results show the short-term and long-term plasticity can induce different bifurcation scenarios when the coupling strength increases. In addition, memristor synaptic plasticity has a great influence on the distribution of firing patterns in the parameter space. More interestingly, when exploring the synchronous firing behavior of two neurons, the two neurons can gradually achieve phase synchronization as the coupling strength increases along the opposite directions under two different memory attributes. Finally, a microcontroller-based hardware system is implemented to verify the numerical simulation results.

Cascades of heterodimensional cycles via period doubling

A heterodimensional cycle is formed by the intersection of stable and unstable manifolds of two saddle periodic orbits that have unstable manifolds of different dimensions: connecting orbits exist from one periodic orbit to the other, and vice versa. The difference in dimensions of the invariant manifolds can only be achieved in vector fields of dimension at least four. At least one of the connecting orbits of the heterodimensional cycle will necessarily be structurally unstable, meaning that is does not persist under small perturbations. Nevertheless, the theory states that the existence of a heterodimensional cycle is generally a robust phenomenon: any sufficiently close vector field (in the C1-topology) also has a heterodimensional cycle.We investigate a particular four-dimensional vector field that is known to have a heterodimensional cycle. We continue this cycle as a codimension-one invariant set in a two-parameter plane. Our investigations make extensive use of advanced numerical methods that prove to be an important tool for uncovering the dynamics and providing insight into the underlying geometric structure. We study changes in the family of connecting orbits as two parameters vary and Floquet multipliers of the periodic orbits in the heterodimensional cycle change. In particular the Floquet multipliers of one of the periodic orbits change from real positive to real negative prior to a period-doubling bifurcation. We then focus on the transitions that occur near this period-doubling bifurcation and find that it generates new families of heterodimensional cycles with different geometric properties. Our careful numerical study suggests that further two-parameter continuation of the ‘period-doubled heterodimensional cycles’ gives rise to an abundance of heterodimensional cycles of different types in the limit of a period-doubling cascade.Our results for this particular example vector field make a contribution to the emerging bifurcation theory of heterodimensional cycles. In particular, the bifurcation scenario we present can be viewed as a specific mechanism behind so-called stabilisation of a heterodimensional cycle via the embedding of one of its constituent periodic orbits into a more complex invariant set.

Effect of Hopf-Hopf bifurcation on the post-flutter behavior of a three-degree-of-freedom airfoil

The post-flutter dynamics of a three-degree-of-freedom nonlinear airfoil with unsteady aerodynamics are investigated based on the Hopf-Hopf bifurcation theory. Many prior works have relied on Hopf bifurcation theory to predict flutter behavior in airfoil systems. Although this approach facilitates flutter prediction and characterization of post-flutter behavior in numerous scenarios, it may be invalid in specific instances, such as when the Hopf bifurcation of the system degenerates. Therefore, this study focuses on a classical degenerate case of Hopf bifurcation in airfoil systems, specifically the Hopf-Hopf bifurcation. We show that the system undergoes various Hopf-Hopf bifurcations under specific parameter conditions as the center of gravity shifts. The local dynamics near the Hopf-Hopf bifurcation points are represented, including quasiperiodic oscillations on a three-dimensional torus. The airfoil begins to oscillate quasiperiodically after the airflow speed crosses specific Hopf-Hopf bifurcation points. The study also uncovers complex quasiperiodic crises and quasiperiodic hysteresis loops, which have not been reported in previous studies of aeroelastic systems. Then, many singularities and bifurcation curves are obtained near the Hopf-Hopf bifurcation point by semiglobal unfolding. Furthermore, the influence of stall effects on the bifurcation structure of the system is represented. It is shown that the types of Hopf-Hopf bifurcations may vary with the changes of stall effects, influencing the system's semiglobal bifurcation structures consequently. For all Hopf-Hopf bifurcation scenarios, stall effects affect one of the Neimark-Sacker bifurcation curve structures unfolded from the Hopf-Hopf bifurcation point significantly, while the other Neimark-Sacker bifurcation curve experiences minimal impact from stall effects. Moreover, a large nonlinear stall coefficient will postpone the onset of quasiperiodic/chaotic oscillations.

Analyzing bifurcation, stability, and wave solutions in nonlinear telecommunications models using transmission lines, Hamiltonian and Jacobian techniques

This study presents a comprehensive analysis of a nonlinear telecommunications model, exploring bifurcation, stability, and wave solutions using Hamiltonian and Jacobian techniques. The investigation begins with a thorough examination of bifurcation behavior, identifying critical points and their stability characteristics, leading to the discovery of diverse bifurcation scenarios. The stability of critical points is further assessed through graphical and numerical methods, highlighting the sensitivity to parameter variations. The study delves into the derivation of both numerical and analytical wave solutions, aligning them with energy orbits depicted in phase portraits, revealing a spectrum of wave behaviors. Additionally, the analysis extends to traveling wave solutions, providing insights into wave propagation dynamics. Notably, the study underscores the efficacy of the planar dynamical approach in capturing system behavior in harmony with phase portrait orbits. The findings have significant implications for telecommunications engineers and researchers, offering insights into system behavior, stability, and signal propagation, ultimately advancing our understanding of complex nonlinear dynamics in telecommunications networks.

Prey group defense and hunting cooperation among generalist-predators induce complex dynamics: a mathematical study.

Group defense in prey and hunting cooperation in predators are two important ecological phenomena and can occur concurrently. In this article, we consider cooperative hunting in generalist predators and group defense in prey under a mathematical framework to comprehend the enormous diversity the model could capture. To do so, we consider a modified Holling-Tanner model where we implement Holling type IV functional response to characterize grazing pattern of predators where prey species exhibit group defense. Additionally, we allow a modification in the attack rate of predators to quantify the hunting cooperation among them. The model admits three boundary equilibria and up to three coexistence equilibrium points. The geometry of the nontrivial prey and predator nullclines and thus the number of coexistence equilibria primarily depends on a specific threshold of the availability of alternative food for predators. We use linear stability analysis to determine the types of hyperbolic equilibrium points and characterize the non-hyperbolic equilibrium points through normal form and center manifold theory. Change in the model parameters leading to the occurrences of a series of local bifurcations from non-hyperbolic equilibrium points, namely, transcritical, saddle-node, Hopf, cusp and Bogdanov-Takens bifurcation; there are also occurrences of global bifurcations such as homoclinic bifurcation and saddle-node bifurcation of limit cycles. We observe two interesting closed 'bubble' form induced by global bifurcations due to change in the strength of hunting cooperation and the availability of alternative food for predators. A three dimensional bifurcation diagram, concerning the original system parameters, captures how the alternation in model formulation induces gradual changes in the bifurcation scenarios. Our model highlights the stabilizing effects of group or gregarious behaviour in both prey and predator, hence supporting the predator-herbivore regulation hypothesis. Additionally, our model highlights the occurrence of "saltatory equilibria" in ecological systems and capture the dynamics observed for lion-herbivore interactions.

Unusual bifurcation scenario in a stably stratified, valley-shaped enclosure heated from below

Unusual bifurcation scenario in a stably stratified, valley-shaped enclosure heated from below

Spiral attractors in a reduced mean-field model of neuron-glial interaction.

This paper investigates various bifurcation scenarios of the appearance of bursting activity in the phenomenological mean-field model of neuron-glial interactions. In particular, we show that the homoclinic spiral attractors in this system can be the source of several types of bursting activity with different properties.

Bifurcation analysis in the regularized four-sided cavity flow: Equilibrium states in a D2-symmetric fluid system

In this study, we investigate a handful of equilibrium states and their hydrodynamic stability for the incompressible flow in a square cavity driven by the tangential simultaneous motion of all of its lids at the same speed. This problem has been investigated in the recent past, albeit only partially and always disregarding the singular boundary conditions at the corners. The implications of the system symmetries have not been rigorously addressed, either. In our study, we introduce a regularized version of the boundary conditions to avoid the corner singularities, in a way that preserves the main features of the original setup. In addition, the Navier–Stokes equations are discretized by means of highly accurate Chebyshev spectral methods that provide exponential convergence of all computed flows and consistently eliminate any potential source of structural instability of the bifurcation scenario. We employ Newton–Krylov solvers, implemented within continuation algorithms, to accurately compute equilibrium solutions. Linear stability analysis of both the primary symmetric base flow and the secondary asymmetric states uncovers new branches of fully asymmetric steady states. The analysis has allowed identification of six previously undetected bifurcations, all of which are associated with the disruption of either the rotational invariance or the reflection symmetry. Some of these bifurcations have been found to be quite clustered in some regions of the parameter space, which points at the underlying action of higher codimension mechanisms. Notably, all the bifurcations reported occur within the range of low to moderate Reynolds numbers, making the regularized four-sided lid-driven cavity flow a reliable benchmark for assessing different numerical schemes in the context of Navier–Stokes equivariant bifurcation theory.

Bifurcation of solutions of Ornstein–Zernike equation in the gas–liquid coexistence region of simple fluids

One objective of this paper is to investigate bifurcation of solutions of Ornstein–Zernike equation in liquid state theory for the gas–liquid coexistence region of simple fluids. The analysis uses a discrete form of the Ornstein–Zernike equation in real space, together with a closure relation which provides nearly self consistent thermodynamic properties of fluids governed by interaction potentials like the Lennard-Jones potential. Accurate asymptotic forms of correlation functions are incorporated in the real space algorithm so as to mitigate the influence of cutoff length used in the discrete approximations. The global error reduction of this algorithm is similar to that of Simpson’s rule. It is found that there is no spinodal on the vapor side on a sub-critical isotherm. However, there is a density at which compressibility vanishes, but this lies on a nonphysical branch. There are fold bifurcation points on the vapor and liquid sides together with a ‘no-real-solution-region’ in between. This is followed by a region of negative compressibility extending up to the spinodal on the liquid side. Importantly, it is found that in regions where compressibility is non-positive, there are infinite number of solutions to OZE which violate an integral constraint that is needed for applications to thermodynamics. The emerging picture is somewhat similar to that in the case of the hypernetted-chain-closure. It is also found that complex solutions connect the physical solutions that exist outside the ‘no-real-solution-region’ as well as in region of negative compressibility. The second objective of the paper, in the light of bifurcation scenario, is to evolve an interpretation of the coupling-parameter expansion of correlation functions pertaining to the region of negative compressibility. It is shown that the expansion provides an accurate description of results of the restricted ensembles, implied in mean field theories and simulation methods. Bifurcation of singlet density to an in-homogeneous distribution at a spinodal is also established using a simple equation, which uses the direct correlation function. Thus the paper re-establishes the fact that restriction to a homogeneous singlet density is the root cause of all inconsistencies found in the gas–liquid coexistence region.

A comprehensive study of spatial spread and multiple time delay in an eco-epidemiological model with infected prey

This paper studies the dynamics of interacting Tilapia fish and Pelican bird population in the Salton Sea. We assume that the diseases spread in Tilapia fish follows the Holling type II response function, and the interaction between Tilapia and Pelican follows the Beddington–DeAngelis response function. The dynamics of diffusive and delayed system are discussed separately. Analytically, all the feasible equilibria and their stability are discussed. The criterion for Turing instability is derived. Based on the normal form theory and center manifold arguments, the existence of stability criterion and the direction of Hopf bifurcation are obtained. Numerical simulation shows the occurrence Hopf bifurcation, double Hopf bifurcation and transcritical bifurcation scenarios. The snap shot shows the spot, spot-strip mix patterns in the whole domain. Further, the stability switching phenomena is observed in the delayed system. Our comprehensive study highlights the effect of different parameters, multiple time delay and extinction in Pelican populations.

Breathing and switching cyclops states in Kuramoto networks with higher-mode coupling.

Cyclops states are intriguing cluster patterns observed in oscillator networks, including neuronal ensembles. The concept of cyclops states formed by two distinct, coherent clusters and a solitary oscillator was introduced by Munyaev etal. [Phys. Rev. Lett. 130, 107201 (2023)0031-900710.1103/PhysRevLett.130.107201], where we explored the surprising prevalence of such states in repulsive Kuramoto networks of rotators with higher-mode harmonics in the coupling. This paper extends our analysis to understand the mechanisms responsible for destroying the cyclops' states and inducing dynamical patterns called breathing and switching cyclops states. We first analytically study the existence and stability of cyclops states in the Kuramoto-Sakaguchi networks of two-dimensional oscillators with inertia as a function of the second coupling harmonic. We then describe two bifurcation scenarios that give birth to breathing and switching cyclops states. We demonstrate that these states and their hybrids are prevalent across a wide coupling range and are robust against a relatively large intrinsic frequency detuning. Beyond the Kuramoto networks, breathing and switching cyclops states promise to strongly manifest in other physical and biological networks, including coupled theta neurons.

Planar quartic–quadratic fold–fold singularity of Filippov systems and its bifurcation

In this paper it is exhibited the bifurcation scenario concerning a typical singularity of planar piecewise smooth vector fields in two zones. This singularity is characterized by a quadratic contact of one vector field and a quartic contact of another vector field at the same point of the switching manifold. By means of a three parameter perturbation, we observe the presence of bifurcations like: Hopf Bifurcation, Sliding Hopf Bifurcation, the birth and annihilation of chaotic scenarios, among others. The two and three dimensional bifurcation diagrams are included.

Scenarios for the appearance of strange attractors in a model of three interacting microbubble contrast agents

We study nonlinear dynamics in a model of three interacting encapsulated gas bubbles in a liquid. The model is a system of three coupled nonlinear oscillators with an external periodic force. Such bubbles have numerous applications, for instance, they are used as contrast agents in ultrasound visualization. Certain types of bubbles oscillations may be beneficial or undesirable depending on a given application and, hence, the dependence of the regimes of bubbles oscillations on the control parameters is worth studying. We demonstrate that there is a wide variety of types of dynamics in the model by constructing a chart of dynamical regimes in the control parameters space. Here we focus on hyperchaotic attractors characterized by three positive Lyapunov exponents and strange attractors with one or two positive Lyapunov exponents possessing an additional zero Lyapunov exponent, which have not been observed previously in the context of bubbles oscillations. We also believe that we provide a first example of a hyperchaotic attractor with additional zero Lyapunov exponent. Furthermore, the mechanisms of the onset of these types of attractors are still not well studied. We identify two-parametric regions in the control parameter space where these hyperchaotic and chaotic attractors appear and study one-parametric routes leading to them. We associate the appearance of hyperchaotic attractors with three positive Lyapunov exponents with the inclusion of a periodic orbit with a three-dimensional unstable manifold, while the onset of chaotic oscillations with an additional zero Lyapunov exponent is connected to the partial synchronization of bubbles oscillations. We propose several underlying bifurcation mechanisms that explain the emergence of these regimes. We believe that these bifurcation scenarios are universal and can be observed in other systems of coupled oscillators.

Regularization of the Boundary Equilibrium Bifurcation in Filippov System with Rich Discontinuity Boundaries

This paper studies a particular type of planar Filippov system that consists of two discontinuity boundaries separating the phase plane into three disjoint regions with different dynamics. This type of system has wide applications in various subjects. As an illustration, a plant disease model and an avian-only model are presented, and their bifurcation scenarios are investigated. By means of the regularization approach, the blowing up method, and the singular perturbation theory, we provide a different way to analyze the dynamics of this type of Filippov system. In particular, the boundary equilibrium bifurcations of such systems are studied. As a consequence, the nonsmooth fold bifurcation becomes a saddle-node bifurcation, while the persistence bifurcation disappears after regularization.

Catastrophic and noncatastrophic population crashes in a bitrophic system with dynamic additional food provision to cooperative predators.

In this article we contemplate the dynamics of an additional food-provided prey-predator system. We assume that the behavior of cooperative predators induces fear in prey, which radically affects the prey's birth and death rates. We observe that the structural instability imposed by strong cooperative hunting among predators goes away with higher intensities of fear levels affecting the prey's reproductive output and mortality. High levels of prey refuge are not conducive to the survival of predators. In such a situation, adequate supply of high-quality additional food is favorable regarding the persistence and stability of the system. Interestingly, the system potentially exhibits two stable configurations under identical ecological conditions by allowing different bifurcation scenarios, including saddle-node and backward bifurcations, and associated hysteresis effects with prey refuge along with additional food quantity and quality. In the stochastic environment, the system experiences critical transitions through bifurcation-induced tipping events with time-varying additional food for predators. Enhanced disturbance events promote noise-induced switching and tipping events. Finally, our investigation explores whether impending population crashes resulting from the variability of additional food quantity and quality can reliably be predicted using early warning signals in the context of redshifted noise. Overall, our results may provide insights for finding control strategies in the context of community ecology.

Bifurcation analysis of autonomous and nonautonomous modified Leslie-Gower models.

In ecological systems, the predator-induced fear dampens the prey's birth rate; yet, it fails to extinguish their population, as they endure and survive even under significant fear-induced costs. In this study, we unveil a modified Leslie-Gower predator-prey model by incorporating the fear of predators, cooperative hunting, and predator-taxis sensitivity. We embark upon an exploration of the positivity and boundedness of solutions, unearthing ecologically viable equilibrium points and their stability conditions governed by the model parameters. Delving deeper, we unravel the scenario of transcritical, saddle-node, Hopf, Bogdanov-Takens, and generalized-Hopf bifurcations within the system's intricate dynamics. Additionally, we observe the bistable nature of the system under some parametric conditions. Further, the nonautonomous extension of our model introduces the intriguing interplay of seasonality in some crucial parameters. We establish a set of sufficient conditions that guarantee the permanence of the seasonally driven system. By conducting a numerical study on the seasonally forced model, we observe a myriad of phenomena manifesting the predator-prey dynamics. Notably, periodic solutions, higher periodic solutions, and bursting patterns emerge, alongside intriguing chaotic dynamics. Specifically, seasonal variations of the predator-taxis sensitivity and hunting cooperation can lead to the extinction of prey species and even the control of chaotic (higher periodic) solutions through the generation of a simple periodic solution. Remarkably, the seasonal forcing has the capacity to govern the chaotic behavior, leading to an exceptionally quasi-periodic arrangement in both prey and predator populations.

On stable and quasi-chaotic regimes in a one-dimensional unimodal mapping obtained by modeling the dynamics of a biological population

The paper considers the properties of the difference equation describing the dynamics of the animal population, obtained earlier in the framework of studies of tundra communities. We have considered a special case in which the model is represented by a one-parameter difference equation that defines a one-dimensional unimodal mapping of a segment into itself, similar to the well-known triangular (tent) mapping, supplemented by a region with a constant value. A change in the mapping parameter generates a bifurcation scenario, in which stability zones arise, characterized by orbits of a constant period, interspersed with zones with more complicated, “quasi-chaotic” regimes. Based on the properties of the n-iterated triangular mapping, a necessary and sufficient condition for the localization of cyclic orbits in the considered type of unimodal mappings is formulated, which makes it possible to identify stability regions for any given period n. On the basis of this condition an algorithm for detecting stability zones is proposed. The main subject of the study is the fractal properties of the set, which is the complement of the obtained set of stability regions to the entire domain of mapping definition. The dynamics of the DH (n) value is obtained, the limit of which at n → ∞ is equal to the fractal dimension dH . It is shown that in the studied range of n (2 ≤ n ≤ 22), the DH (n) < 0.9, which suggest that dH < 1. If so, then according to the definition of a fractal set, its topological dimension is dT = 0, which means that the complement of the set of stability regions consists of isolated points.