Patterns and Stability of Coupled Multi-Stable Nonlinear Oscillators

Nonlinear isolated and coupled oscillators are extensively studied as prototypical nonlinear dynamics models. Much attention has been devoted to oscillator synchronization or the lack thereof. Here, we study the synchronization and stability of coupled driven-damped Helmholtz-Duffing oscillators in bi-stability regimes. We find that despite the fact that the system parameters and the driving force are identical, the stability of the two states to spatially non-uniform perturbations is very different. Moreover, the final stable states, resulting from these spatial perturbations, are not solely dictated by the wavelength of the perturbing mode and take different spatial configurations in terms of the coupled oscillator phases.

A coupled algebraic-delay differential system modeling tick-host interactive behavioural dynamics and multi-stability

We propose a coupled system of delay-algebraic equations to describe tick attaching and host grooming behaviors in the tick-host interface, and use the model to understand how this tick-host interaction impacts the tick population dynamics. We consider two critical state variables, the loads of feeding ticks on host and the engorged ticks on the ground for ticks in a particular development stage (nymphal stage) and show that the model as a coupled system of delay differential equation and an algebraic (integral) equation may have rich structures of equilibrium states, leading to multi-stability. We perform asymptotic analyses and use the implicit function theorem to characterize the stability of these equilibrium states, and show that bi-stability and quadri-stability occur naturally in several combinations of tick attaching and host grooming behaviours. In particular, we show that in the case when host grooming is triggered by the tick biting, the system will have three stable equilibrium states including the extinction state, and two unstable equilibrium states. In addition, the two nontrivial stable equilibrium states correspond to a low attachment rate and a large number of feeding ticks, and a high attachment rate and a small number of feeding ticks, respectively.

Bifurcations and hydra effects in Bazykin's predator-prey model

We consider Bazykin's model to address harvesting induced stability exchanges through bifurcation analysis. We examine the existence of hydra effects and analyze the stock pattern under predator harvesting. Prey harvesting cannot produce hydra effects in our model, whereas predator harvesting may cause multiple hydra effects. Our study reveals that type II response function and mutual interference among predators jointly induce multiple hydra effects and bistability. Bifurcations such as single Hopf-bifurcation, multiple Hopf-bifurcations and multiple saddle-node bifurcations appear for increasing harvesting rate on the predators. However, over-exploitation of the predators cannot generate any such bifurcation in our study. In simulations, the maximum sustainable yield (MSY) exists at a globally stable state. When predator is culled under increasing effort, basin of attraction of the equilibrium corresponding to the higher predator stock gets expanded, which alternatively is in favor of stock benefit for predators. The ecological theory developed in this study might be useful to understand conservation policy and fishery management.

Mathematical Modeling of Plasticity and Heterogeneity in EMT

The epithelial-mesenchymal transition (EMT) and the corresponding reverse process, mesenchymal-epithelial transition (MET), are dynamic and reversible cellular programs orchestrated by many changes at both biochemical and morphological levels. A recent surge in identifying the molecular mechanisms underlying EMT/MET has led to the development of various mathematical models that have contributed to our improved understanding of dynamics at single-cell and population levels: (a) multi-stability-how many phenotypes can cells attain during an EMT/MET?, (b) reversibility/irreversibility-what time and/or concentration of an EMT inducer marks the "tipping point" when cells induced to undergo EMT cannot revert?, (c) symmetry in EMT/MET-do cells take the same path when reverting as they took during the induction of EMT?, and (d) non-cell autonomous mechanisms-how does a cell undergoing EMT alter the tendency of its neighbors to undergo EMT? These dynamical traits may facilitate a heterogenous response within a cell population undergoing EMT/MET. Here, we present a few examples of designing different mathematical models that can contribute to decoding EMT/MET dynamics.

Breaking Away: The Role of Homeostatic Drive in Perpetuating Depression

We propose that the complexity of regulatory interactions modulating brain neurochemistry and behavior is such that multiple stable responses may be supported, and that some of these alternate regulatory programs may play a role in perpetuating persistent psychological dysfunction. To explore this, we constructed a model network representing major neurotransmission and behavioral mechanisms reported in literature as discrete logic circuits. Connectivity and information flow through this biobehavioral circuitry supported two distinct and stable regulatory programs. One such program perpetuated a depressive state with a characteristic neurochemical signature including low serotonin. Further analysis suggested that small irregularities in glutamate levels may render this pathology more directly accessible. Computer simulations mimicking selective serotonin reuptake inhibitor (SSRI) therapy in the presence of everyday stressors predicted recidivism rates similar to those reported clinically and highlighted the potentially significant benefit of concurrent behavioral stress management therapy.

Regime shifts under forcing of non-stationary attractors: Conceptual model and case studies in hydrologic systems

We present here a conceptual model and analysis of complex systems using hypothetical cases of regime shifts resulting from temporal non-stationarity in attractor strengths, and then present selected published cases to illustrate such regime shifts in hydrologic systems (shallow aquatic ecosystems; water table shifts; soil salinization). Complex systems are dynamic and can exist in two or more stable states (or regimes). Temporal variations in state variables occur in response to fluctuations in external forcing, which are modulated by interactions among internal processes. Combined effects of external forcing and non-stationary strengths of alternative attractors can lead to shifts from original to alternate regimes. In systems with bi-stable states, when the strengths of two competing attractors are constant in time, or are non-stationary but change in a linear fashion, regime shifts are found to be temporally stationary and only controlled by the characteristics of the external forcing. However, when attractor strengths change in time non-linearly or vary stochastically, regime shifts in complex systems are characterized by non-stationary probability density functions (pdfs). We briefly discuss implications and challenges to prediction and management of hydrologic complex systems.