Bistable Region Research Articles

Electrical spectroscopy of the spin-wave dispersion and bistability in gallium-doped yttrium iron garnet

Yttrium iron garnet (YIG) is a magnetic insulator with record-low damping, allowing spin-wave transport over macroscopic distances. Doping YIG with gallium ions greatly reduces the demagnetizing field and introduces a perpendicular magnetic anisotropy, which leads to an isotropic spin-wave dispersion that facilitates spin-wave optics and spin-wave steering. Here, we characterize the dispersion of a gallium-doped YIG (Ga:YIG) thin film using electrical spectroscopy. We determine the magnetic anisotropy parameters and Gilbert damping from the frequency and linewidth of the ferromagnetic resonance (FMR). Next, we use propagating spin wave spectroscopy in the Damon–Eshbach configuration to detect the small spin-wave magnetic fields of this ultrathin weak magnet over a wide range of wavevectors, enabling the extraction of the exchange constant αex=1.32×10−12 J/m. We observe foldover of the FMR with increasing drive power, leading to frequency shifts of the spin-wave modes and a bistable region in the spin-wave spectra. Our results shed light on isotropic spin-wave transport in Ga:YIG and highlight the potential of electrical spectroscopy to map out the dispersion and bistability of propagating spin waves in magnets with a low saturation magnetization.

Modeling the Dynamics of a Gyroscopic Rigid Rotor with Linear and Nonlinear Damping and Nonlinear Stiffness of the Elastic Support

This study analytically and numerically modeled the dynamics of a gyroscopic rigid rotor with linear and nonlinear cubic damping and nonlinear cubic stiffness of an elastic support. It has been shown that (i) joint linear and nonlinear cubic damping significantly suppresses the vibration amplitude (including the maximum) in the resonant velocity region and beyond it, and (ii) joint linear and nonlinear cubic damping more effectively affects the boundaries of the bistability region by its narrowing than linear damping. A methodology is proposed for determining and identifying the coefficients of nonlinear stiffness, linear damping, and nonlinear cubic damping of the support material, where jump-like effects are eliminated. Damping also affects the stability of motion; if linear damping shifts the left boundary of the instability region towards large amplitudes and speeds of rotation of the shaft, then nonlinear cubic damping can completely eliminate it. The varying amplitude (VAM) method is used to determine the nature of the system response, supplemented with the concept of “slow” time, which allows us to investigate and analyze the effect of nonlinear cubic damping and nonlinear rigidity of cubic order on the frequency response at a nonstationary resonant transition.

Bistabilities in two parallel Kármán wakes

Bistabilities of two equilibrium states discovered in the coupled side-by-side Kármán wakes are investigated through Floquet analysis and direct numerical simulation (DNS) with different initial conditions over a range of gap-to-diameter ratio ( $g^*= 0.2\text {--}3.5$ ) and Reynolds number ( $Re = 47\text {--}100$ ). Two bistabilities are found in the transitional $g^*-Re$ regions from in-phase (IP) to anti-phase (AP) vortex shedding states. By initialising the flow in DNS with zero initial conditions, the flow in the first bistable region (i.e. bistable IP/FF $_C$ at $g^*= 1.4 \text {--} 2.0$ , where FF $_C$ denotes the conditional flip-flop flow) attains flip-flop (FF) flow, it settles into the IP state by initialising the flow with an IP flow. The second bistability is observed between cylinder-scale IP and AP states at large $g^*$ ( $=$ 2.0–3.5). The transition from the FF $_C$ to IP is dependent on initial conditions and irreversible over the parameter space, meaning that the flow cannot revert back to the FF $_C$ state once it jumps to the IP state irrespective of the direction of $Re$ variations. Its counterpart for the bistable IP/AP state is reversible. We also found that the FF $_C$ flow in the first bistable region is primarily bifurcated from synchronised AP with cluster-scale features, possibly because the cluster-scale AP flow is inherently unstable to FF flow instabilities. It is demonstrated that the irreversible bistability exists in other interacting wakes around multiple cylinders. A good understanding of flow bistabilities is pivotal to flow control applications and the interpretation of desynchronised flow features observed at high $Re$ values.

Modeling and Online Adaptation of ALOHA for Low-Power Wide-Area Networks (LPWANs)

Unslotted ALOHA protocol has been adopted as a channel access mechanism in commercial low-power wide-area networks (LPWANs), such as Sigfox and long-range (LoRa) alliance. This work examines the throughput and random access (RA) delay distribution of unslotted ALOHA systems by considering exponential random backoff (ERB) or uniform random backoff (URB) algorithm. We further characterize the operating region of the systems as unsaturated stable, bistable, and saturated regions in terms of the new packet arrival and retransmission rates. To run the system stably with the maximum throughput, we propose a Bayesian online backoff algorithm that estimates the number of backlogged devices. Its performance is compared with other algorithms, such as particle filter (PF)-based algorithm, binary exponential backoff (BEB) algorithm, and the algorithm of exploiting exact backlog size information. Through extensive simulations, it is demonstrated that the performance of the proposed algorithm is very close to the upper bound and robust to time-varying traffic condition.

The optical response of a complex nanodimer composed of a metal nanoparticle and a quantum dot: The effect of nanoparticle geometry

The optical response of a hybrid structure consisting of a metal nanoparticle in the vicinity of a quantum dot under an external field is studied. The effect of the nanoparticle shape as it evolves from a sphere to a rod is investigated for various structural parameters. Significant changes are observed in the phase diagram of the optical response at the borders between Fano, exciton-induced transparency, suppression, and bistable regions. In the appropriate limits, the configuration exhibits a response where notable variations are particularly observed in suppression and bistable regions. In the presence of the nanorod, the suppression region is found as the beginning of bistability. Our calculations show that the nanorod can generate significant nonlinearities at the site of quantum dot in the bistable region, making the hybrid structure suitable for devices involving optical switching or information transfer. By considering the effect of other structural parameters, for example the configuration state, background dielectric constant, and interparticle distance, one can obtain additional control of the optical response of the system, engineering it to achieve particular objectives.

Dynamical analysis of early afterdepolarization patterns in a biophysically detailed cardiac model.

Arrhythmogenic early afterdepolarizations (EADs) are investigated in a biophysically detailed mathematical model of a rabbit ventricular myocyte, providing their location in the parameter phase space and describing their dynamical mechanisms. Simulations using the Sato model, defined by 27 state variables and 177 parameters, are conducted to generate electrical action potentials (APs) for different values of the pacing cycle length and other parameters related to sodium and calcium concentrations. A detailed study of the different AP patterns with or without EADs is carried out, showing the presence of a high variety of temporal AP configurations with chaotic and quasiperiodic behaviors. Regions of bistability are identified and, importantly, linked to transitions between different behaviors. Using sweeping techniques, one-, two-, and three-parameter phase spaces are provided, allowing ascertainment of the role of the selected parameters as well as location of the transition regions. A Devil's staircase, with symbolic sequence analysis, is proposed to describe transitions in the ratio between the number of voltage (EAD and AP) peaks and the number of APs. To conclude, the obtained results are linked to recent studies for low-dimensional models and a conjecture is made for the internal dynamical structure of the transition region from non-EAD to EAD behavior using fold and cusp bifurcations and maximal canards.

Multilayer social reinforcement induces bistability on multiplex networks

The social reinforcement mechanism, which characterizes the promoting effects when exposed to multiple sources in the social contagion process, is ubiquitous in information technology ecosystems and has aroused great attention in recent years. While the impacts of social reinforcement on single-layer networks are well documented, extension to multilayer networks is needed to study how reinforcement from different social circles influences the spreading dynamics. To this end, we incorporate multilayer social reinforcement into an ignorant–spreader–ignorant model on multiplex networks. Our theoretical analysis combines the pairwise method and mean-field theory and agrees well with large-scale simulations. Surprisingly, we find this complex social contagion mechanism triggers the emergence of bistability phenomena, where extinction and outbreak states coexist. In particular, the hysteresis loop of stationary prevalence occurs in this bistable region, explaining why the fight against the spread of rumors is protracted and difficult in modern society. Further, we show that the final state of bistable regions depends on the initial density of adopters, the critical value of which decreases as the contagion transmissibility or the multilayer reinforcement increases. In particular, we highlight two possible conditions for the outbreak of social contagion: to possess large contagion transmissibility, or to possess a large initial density of adopters with strong multilayer reinforcement. Our results unveil the non-negligible power of social reinforcement on multiplex networks, which sheds lights on designing efficient strategies in spreading behaviors such as marketing and promoting innovations.

Network clique cover approximation to analyze complex contagions through group interactions

Contagion processes have been proven to fundamentally depend on the structural properties of the interaction networks conveying them. Many real networked systems are characterized by clustered substructures representing either collections of all-to-all pair-wise interactions (cliques) and/or group interactions, involving many of their members at once. In this work, focusing on interaction structures represented as simplicial complexes, we present a discrete-time microscopic model of complex contagion for a susceptible-infected-susceptible dynamics. Introducing a particular edge clique cover and a heuristic to find it, the model accounts for the higher-order dynamical correlations among the members of the substructures (cliques/simplices). The analytical computation of the critical point reveals that higher-order correlations are responsible for its dependence on the higher-order couplings. While such dependence eludes any mean-field model, the possibility of a bi-stable region is extended to structured populations.

ROLE OF REPEAT INFECTION IN THE DYNAMICS OF A SIMPLE MODEL OF WANING AND BOOSTING IMMUNITY

Some infectious diseases produce lifelong immunity while others only produce temporary immunity. In the case of short-lived immunity, the level of protection wanes over time and may be boosted upon re-exposure, via infection or vaccination. Previous work developed a simple model capturing waning and boosting immunity, known as the Susceptible-Infectious-Recovered-Waned-Susceptible (SIRWS) model, which exhibits rich dynamical behavior including supercritical and subcritical Hopf bifurcations among other structures. Here, we extend the bifurcation analyses of the SIRWS model to examine the influence of all parameters on these bifurcation structures. We show that the bistable region, involving both a stable fixed point and a stable limit cycle, exists only for a small region of biologically realistic parameter space. Furthermore, we contrast the SIRWS model with a modified version, where immune boosting may involve the occurrence of a secondary infection. Analysis of this extended model shows that oscillations and bistability, as found in the SIRWS model, depend on strong assumptions about infectivity and recovery rate from secondary infection. Understanding the dynamics of models of waning and boosting immunity is important for accurately assessing epidemiological data.

Effect of time-delayed feedback in a bistable system inferred by logic operation

• The effect of time delay in bistable system can be inferred by logic operation. • Time-delayed negative feedback can reduce the width of bistable region. • There is an unstable limit cycle on the directly connected boundary between the oscillatory and suprathreshold region. • In the alternate region, the system is sensitive to parameters and can enter two potential wells alternately with the slight variation in sensitive parameters. The effect of time delay is inferred by studying the influence of time-delayed feedback on the signal threshold of logic input signal, which reflects the width of the bistable region. Time-delayed positive feedback does not alter the width of the bistable region but increases the switching time. Time-delayed negative feedback leads to decreases in the depth of two potential wells and the width of the bistable region of the system, resulting in instability and oscillation of the system. Then, the parameter plane of the system can be divided into suprathreshold, subthreshold, oscillatory and alternate regions. With increasing time delay and negative feedback intensity, the suprathreshold and oscillatory regions continuously squeeze the middle subthreshold region until they are connected. One part of the oscillatory region is directly connected with the suprathreshold region, while the other part is transiently connected with the alternate region. There is an unstable limit cycle on the directly connected boundary. In the alternate region, the system is sensitive to initial values and system parameters and can enter two potential wells alternately with slight variations in sensitive parameters.

A model for the regulation of apoptosis intrinsic pathway: The potential role of the transcriptional regulator E2F in the point of no return

Apoptosis has been extensively characterized by both experimental approaches and model simulations. However, it is still not fully understood how the regulation occurs, especially in the intrinsic pathway, which can be activated by a great variety of signals. In addition, the conditions in which a point of no return could be reached remain elusive. In this work, we use differential equations models to approach these issues. Our starting point was the model for caspase activation of Legewie et al. (Legewie S, et al., PLoS Computational Biology 2006, 2(9): e120), which exhibits irreversible bistability. We added an activation module to this model, with the main events related to mitochondrial outer membrane permeabilization, which includes cytochrome C release by the mitochondria and its effects on caspase activation and respiratory chain disruption. This “Extended Legewie Model” (ELM) uses BAK as the apoptotic stimulus and active caspase 3 as a measure of apoptosis activation. Unexpectedly, in the extended model, BAK cannot trigger apoptosis activation using physiologically sound initial values of the variables, due to limitations in apoptosome concentration increase. Therefore, the next step was to find a regulatory mechanism, allowing apoptosis activation in the ELM, starting from physiological initial concentrations. For this aim, we performed a sensitivity analysis on the 61 parameters of the system, finding that those producing the most relevant changes in the qualitative behaviour were the rates of synthesis of caspase 3, caspase 9 and XIAP. Based on these results, the transcription factor E2F was included in the ELM because it directly regulates the rate of synthesis of caspase 3 and 9. Depending on the concentration of E2F, the ELM shows different qualitative behaviours. On one hand, for low E2F apoptosis is impossible and for high E2F apoptosis is inevitable. Therefore, if E2F is sufficiently increased, the point of no return is crossed. On the other hand, for intermediate values of E2F there is a bistable region where the fate of the system also depends on the concentration of BAK and other signalling species.

The study of optical bistabiltiy in a hydrodynamic metallic nanoparticle

In this paper, the optical bistability in a plasmonic nanoparticle by using the Hydrodynamic model for conduction electrons is investigated. It is shown that the optical bistable threshold, the bistable region and the shape of hysteresis loop are dependent on the permittivity of the surrounding medium, and the composition of plasmonic nanoparticle.

Simplicial SIRS epidemic models with nonlinear incidence rates.

Mathematical epidemiology that describes the complex dynamics on social networks has become increasingly popular. However, a few methods have tackled the problem of coupling network topology with complex incidence mechanisms. Here, we propose a simplicial susceptible-infected-recovered-susceptible (SIRS) model to investigate the epidemic spreading via combining the network higher-order structure with a nonlinear incidence rate. A network-based social system is reshaped to a simplicial complex, in which the spreading or infection occurs with nonlinear reinforcement characterized by the simplex dimensions. Compared with the previous simplicial susceptible-infected-susceptible (SIS) models, the proposed SIRS model can not only capture the discontinuous transition and the bistability of a complex system but also capture the periodic phenomenon of epidemic outbreaks. More significantly, the two thresholds associated with the bistable region and the critical value of the reinforcement factor are derived. We further analyze the stability of equilibrium points of the proposed model and obtain the condition of existence of the bistable states and limit cycles. This work expands the simplicial SIS models to SIRS models and sheds light on a novel perspective of combining the higher-order structure of complex systems with nonlinear incidence rates.

Effect of time delay in a bistable synthetic gene network.

The essence of logical stochastic resonance is the dynamic manipulation of potential wells. The effect of time delay on the depth of potential wells and the width of a bistable region can be inferred by logic operations in the bistable system with time delay. In a time-delayed synthetic gene network, time delay in the synthesis process can increase the depth of the potential wells, while that in the degradation process, it can reduce the depth of the potential wells, which will result in a decrease in the width of the bistable region (the reason for time delay to induce logic operations without external driving force) and the instability of the system (oscillation). These two opposite effects imply stretching and folding, leading to complex dynamical behaviors of the system, including period, chaos, bubble, chaotic bubble, forward and reverse period doubling bifurcation, intermittency, and coexisting attractors.

Tristability in Viscoelastic Flow Past Side-by-Side Microcylinders.

Viscoelastic flows through microscale porous arrays exhibit complex path selection and switching phenomena. However, understanding this process is limited by a lack of studies linking between a single object and large arrays. Here, we report experiments on viscoelastic flow past side-by-side microcylinders with variable intercylinder gap. With increasing flow rate, a sequence of two imperfect symmetry-breaking bifurcations forces selection of either one or two of the three possible flow paths around the cylinders. Tuning the gap length through the value where the first bifurcation becomes perfect reveals regions of bistability and tristability in a dimensionless flow rate-gap length phase diagram.

Dynamical behaviors and optimal harvesting of an intraguild prey-predator fishery model with Michaelis-Menten type predator harvesting

The present paper discusses the dynamics and optimal harvesting of an intraguild prey-predator fishery model by incorporating the nonlinear Michaelis-Menten type of harvesting in predator. To our knowledge, there is limited literature working on Michaelis-Menten type of harvesting in a three species intraguild model with variable carrying capacity. The carrying capacity is proportional to the density of biotic resource. The existence of the possible equilibria has been studied along with the stability criteria. We consider the impact of predator fish harvesting as the bifurcation parameter to analyze the long time behavior of the proposed system. From the economic perspective, bionomic equilibrium of the system is studied and optimal harvesting policy is derived with the assistance of Pontryagin Maximum Principle. Finally, numerical results are presented to verify our analytical results. It is shown that in the bifurcation diagrams, the system can exhibit transcritical and Hopf bifurcations in the neighborhood of coexistence equilibrium at low and relatively high level of predator harvesting respectively. Interestingly, the system enters to a bistable region where both the coexistence and predator-free equilibria can be stable depending on the initial values. This bistable behavior might be novel to the existing literature that studied intraguild models with variable carrying capacity. The objective of this study is to derive the optimal threshold for the predator harvesting that gives maximum financial profit while sustaining the fishery resources.

Hexagon Invasion Fronts Outside the Homoclinic Snaking Region in the Planar Swift--Hohenberg Equation

Stationary fronts connecting the trivial state and a cellular (distorted) hexagonal pattern in the Swift-Hohenberg equation with a quadratic-cubic nonlinearity are known to undergo a process of infinitely many folds as a parameter is varied, known as homoclinic snaking, where new hexagon cells are added to the core, leading to a region of infinitely-many, co-existing localized states. Outside the homoclinic snaking region, the hexagon fronts can invade the trivial state in a bursting fashion. In this paper, we use a far-field core decomposition to set up a numerical path-following routine to trace out the bifurcation diagrams of hexagon fronts for the two main orientations of cellular hexagon pattern with respect to the interface in the bistable region. We find for one orientation that the hexagon fronts can destabilize as the distorted hexagons are stretched in the transverse direction leading to defects occurring in the deposited cellular pattern. We then plot diagrams showing when the selected fronts for the two main orientations, aligned perpendicular to each other, are compatible leading to a hexagon wavenumber selection prediction for hexagon patches on the plane. Finally, we verify the compatibility criterion for hexagon patches in the Swift-Hohenberg equation with a quadratic-cubic nonlinearity and a large non-variational perturbation. The numerical algorithms presented in this paper can be adapted to general reaction-diffusion systems.

Analysis of Two-Step Random Access Procedure for Cellular Ultra-Reliable Low Latency Communications

Due to the emergence of Internet of Things (IoTs), it can be expected that the bandwidth provided by cellular systems might be consumed up soon. Some applications of them are delay-sensitive such that it would be critical to guarantee random access (RA) delay less than a threshold. Since the existing Long-Term Evolution-Advanced (LTE-A) RA procedure is a four-step signaling procedure, it may not be suitable for such delay-sensitive applications due to its time-consuming procedure. This work investigates a two-step RA procedure for 5G New Radio systems, where RA preamble and bandwidth request message are transmitted at the same time. First we show that the operating region of two-step RA procedure can be divided into three regions such as unsaturated stable, bistable, and saturated regions in terms of a packet generation probability, retransmission probability, the number of devices, and the number of RA preambles. To see whether RA delay requirement of delay-sensitive applications can be guaranteed, this work shows that the system should run under unsaturated region and derives RA delay distribution when IoT devices employ geometric probability backoff (GPB) or uniform window backoff (UWB) algorithm. We then examine the probability that the RA delay would be larger than some threshold depending on the operation regions.

Asynchronous Chaos and Bifurcations in a Model of Two Coupled Identical Hindmarsh – Rose Neurons

We study a minimal network of two coupled neurons described by the Hindmarsh – Rose model with a linear coupling. We suppose that individual neurons are identical and study whether the dynamical regimes of a single neuron would be stable synchronous regimes in the model of two coupled neurons. We find that among synchronous regimes only regular periodic spiking and quiescence are stable in a certain range of parameters, while no bursting synchronous regimes are stable. Moreover, we show that there are no stable synchronous chaotic regimes in the parameter range considered. On the other hand, we find a wide range of parameters in which a stable asynchronous chaotic regime exists. Furthermore, we identify narrow regions of bistability, when synchronous and asynchronous regimes coexist. However, the asynchronous attractor is monostable in a wide range of parameters. We demonstrate that the onset of the asynchronous chaotic attractor occurs according to the Afraimovich – Shilnikov scenario. We have observed various asynchronous firing patterns: irregular quasi-periodic and chaotic spiking, both regular and chaotic bursting regimes, in which the number of spikes per burst varied greatly depending on the control parameter.

Two-parameter bifurcation analysis of an aircraft nose landing gear model

The Hopf–Hopf bifurcations of an aircraft nose landing gear model are investigated in this paper based on two pairs of continuation parameters. It is shown that due to the effect of higher-order terms, there exists a deviation between the bifurcation curve of the original system and the bifurcation curve of its truncated amplitude system, which results in a contraction in the bistable region compared with the previous work. Two codimension 3 strong resonant points can be detected on the Neimark–Sacker bifurcation curve. Period-3 and period-4 limit cycles with torsional-lateral shimmy can emerge when the vertical force and the taxiing velocity are chosen near the resonant points. The occurrence of 1:3 and 1:4 resonant Hopf–Hopf bifurcations can induce tire-ground vibrations in torsional and lateral direction, which can greatly degrade the performance and safety and possibly destroy the airstrip. So it is of significance to investigate the resonate behaviors of the aircraft nose landing gear to ensure the safe landing.