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.

Bistability and nonequilibrium condensation in a driven-dissipative Josephson array: a c-field model

Developing theoretical models for nonequilibrium quantum systems poses significant challenges. Here we develop and study a multimode model of a driven-dissipative Josephson junction chain of atomic Bose-Einstein condensates, as realised in the experiment of Labouvie et al.[Phys. Rev. Lett. 116, 235302 (2016)]. The model is based on c-field theory, a beyond-mean-field approach to Bose-Einstein condensates that incorporates fluctuations due to finite temperature and dissipation. We find the c-field model is capable of capturing all key features of the nonequilibrium phase diagram, including bistability and a critical slowing down in the lower branch of the bistable region. Our model is closely related to the so-called Lugiato-Lefever equation, and thus establishes new connections between nonequilibrium dynamics of ultracold atoms with nonlinear optics, exciton-polariton superfluids, and driven damped sine-Gordon systems.

Can Catastrophe Theory explain expansion and contagious of Covid-19?

Since SARS-Cov-2 started spreading in China and turned into a pandemic disease called Covid-19, many articles about prediction with mathematical model have appeared in the literature. In addition to models in specialized journals, a significant amount of software was made available, presenting with dashboards spreading of the pandemic for each new. These models are solved by computer simulation of traditional exponential models as a representation of the growth of cases and deaths. Some more accurate models are based on existing variations of SIR model (Susceptible, Infected and Recovered). A third class of study is developed in spatial or probabilistic models as a way of forecasting the effect of Covid-19 around the world. Data on the number of positive cases in all countries, show that SARS-Cov-2 shows great resistance even after strategies of lockdown or social distancing. The purpose of this article is to show how the bifurcation theory, known as Catastrophe Theory, can help to understand why Covid-19 expansion rates change so much and even with low values for a longtime trigger contagion quickly and abruptly. The Catastrophe Theory was conceived by the mathematician Rene Thom in the 60s with wide applications in works in the 70s. The outbreak of spruce budworm in Canada revealed a very interesting opportunity to test Catastrophe Theory whose explanation for the phenomenon was widely debated in the academic world. Inspired by the same mathematical approach to this phenomenon in Canada in the 1970s, we applied the Catastrophe Theory in the current Covid-19 pandemic. We observed that sudden outbreaks occur when the carrying capacity and the rate of expansion of the virus reach a region of bifurcation on the cusp surface. With actual Covid-19 data obtained from WHO, we fitted the dynamic model using the particle swarm technique and compared the results in the bifurcation plan with the Covid-19 outbreaks in different regions of the world. It is possible in many cases to observe the trajectory of the parameters between limit points in the bistable region and the consequent explosion of cases observed for each country assessed.

Biocrust as one of multiple stable states in global drylands

Biocrusts cover ~30% of global drylands with a prominent role in the biogeochemical cycles. Theoretically, biocrusts, vascular plants, and bare soil can represent multiple stable states in drylands. However, no empirical evidence for the existence of a biocrust stable state has been reported. Here, using a global drylands dataset, we found that biocrusts form an alternative stable state (biocrust cover, ~80%; vascular cover, ≤10%) besides bare soil (both biocrust and vascular cover, ≤10%) and vascular plants (vascular cover, >50%; biocrust cover, ~5%). The pattern of multiple stable states associated with biocrusts differs from the classic fold bifurcation, and values of the aridity index in the range of 0 to 0.6 define a bistable region where multiple stable states coexist. This study empirically demonstrates the existence and thresholds of multiple stable states associated with biocrusts along climatic gradients and thus may greatly contribute to conservation and restoration of global drylands.

Thermodynamic Analysis of Bistability in Rayleigh–Bénard Convection

Bistability is often encountered in association with dissipative systems far from equilibrium, such as biological, physical, and chemical phenomena. There have been various attempts to theoretically analyze the bistabilities of dissipative systems. However, there is no universal theoretical approach to determine the development of a bistable system far from equilibrium. This study shows that thermodynamic analysis based on entropy production can be used to predict the transition point in the bistable region during Rayleigh–Bénard convection using the experimental relationship between the thermodynamic flux and driving force. The bistable region is characterized by two distinct features: the flux of the second state is higher than that of the first state, and the entropy production of the second state is lower than that of the first state. This thermodynamic interpretation provides new insights that can be used to predict bistable behaviors in various dissipative systems.

Robustness and parameter geography in post-translational modification systems

Biological systems are acknowledged to be robust to perturbations but a rigorous understanding of this has been elusive. In a mathematical model, perturbations often exert their effect through parameters, so sizes and shapes of parametric regions offer an integrated global estimate of robustness. Here, we explore this “parameter geography” for bistability in post-translational modification (PTM) systems. We use the previously developed “linear framework” for timescale separation to describe the steady-states of a two-site PTM system as the solutions of two polynomial equations in two variables, with eight non-dimensional parameters. Importantly, this approach allows us to accommodate enzyme mechanisms of arbitrary complexity beyond the conventional Michaelis-Menten scheme, which unrealistically forbids product rebinding. We further use the numerical algebraic geometry tools Bertini, Paramotopy, and alphaCertified to statistically assess the solutions to these equations at ∼109 parameter points in total. Subject to sampling limitations, we find no bistability when substrate amount is below a threshold relative to enzyme amounts. As substrate increases, the bistable region acquires 8-dimensional volume which increases in an apparently monotonic and sigmoidal manner towards saturation. The region remains connected but not convex, albeit with a high visibility ratio. Surprisingly, the saturating bistable region occupies a much smaller proportion of the sampling domain under mechanistic assumptions more realistic than the Michaelis-Menten scheme. We find that bistability is compromised by product rebinding and that unrealistic assumptions on enzyme mechanisms have obscured its parametric rarity. The apparent monotonic increase in volume of the bistable region remains perplexing because the region itself does not grow monotonically: parameter points can move back and forth between monostability and bistability. We suggest mathematical conjectures and questions arising from these findings. Advances in theory and software now permit insights into parameter geography to be uncovered by high-dimensional, data-centric analysis.Author SummaryBiological organisms are often said to have robust properties but it is difficult to understand how such robustness arises from molecular interactions. Here, we use a mathematical model to study how the molecular mechanism of protein modification exhibits the property of multiple internal states, which has been suggested to underlie memory and decision making. The robustness of this property is revealed by the size and shape, or “geography,” of the parametric region in which the property holds. We use advances in reducing model complexity and in rapidly solving the underlying equations, to extensively sample parameter points in an 8-dimensional space. We find that under realistic molecular assumptions the size of the region is surprisingly small, suggesting that generating multiple internal states with such a mechanism is much harder than expected. While the shape of the region appears straightforward, we find surprising complexity in how the region grows with increasing amounts of the modified substrate. Our approach uses statistical analysis of data generated from a model, rather than from experiments, but leads to precise mathematical conjectures about parameter geography and biological robustness.

A chemically fueled non-enzymatic bistable network

One of the grand challenges in contemporary systems chemistry research is to mimic life-like functions using simple synthetic molecular networks. This is particularly true for systems that are out of chemical equilibrium and show complex dynamic behaviour, such as multi-stability, oscillations and chaos. We report here on thiodepsipeptide-based non-enzymatic networks propelled by reversible replication processes out of equilibrium, displaying bistability. Accordingly, we present quantitative analyses of the bistable behaviour, featuring a phase transition from the simple equilibration processes taking place in reversible dynamic chemistry into the bistable region. This behaviour is observed only when the system is continuously fueled by a reducing agent that keeps it far from equilibrium, and only when operating within a specifically defined parameter space. We propose that the development of biomimetic bistable systems will pave the way towards the study of more elaborate functions, such as information transfer and signalling.

Gradient theory of domain walls in thin, nematic liquid crystals films

In this paper, we describe domain walls appearing in a thin, nematic liquid crystal sample subject to an external field with intensity close to the Fréedericksz transition threshold. Using the gradient theory of the phase transition adapted to this situation, we show that depending on the parameters of the system, domain walls occur in the bistable region or at the border between the bistable and the monostable region.

Bistability and Triggering Analysis of Thermoacoustic Instability in a Horizontal Rijke Tube

Dynamical systems theory has been often employed to study nonlinear flow and flame dynamics in combustion systems. However, the corresponding studies using nonlinear dynamics to analyze the Rijke tube thermoacoustic system are still occasional. Little study has been performed to elucidate the characteristics of triggering phenomenon in the bistable region of the thermoacoustic system. In this regard, the main objectives of the present research are to contribute analysis for the lack of literature in these areas, especially to study the bistability and triggering properties of a thermoacoustic system. The thermoacoustic model of a horizontal Rijke tube is firstly established. The governing equations are expanded and solved by using Galerkin method. The analysis of the system is carried out by using nonlinear dynamics theory. Linear and nonlinear stability boundaries for the variation of non-dimensional heater power, damping coefficient and the relative heater location are obtained for different values of non-dimensional time lag in the system. Regions of global stability, global instability and bistability are characterized. The bistable region in the relative heater location is distributed symmetrically with xf=0.25. It is observed that the bistable region in the relative heater location firstly decreases with an increase in the non-dimensional time lag, reaching a minimum value at a certain critical value of τ, then increases. The situation for the bistable region in the damping coefficient presents a reverse variation, And the bistable region reach the maximum at τ=0.5. The triggering phenomenon and limit cycle of the system in the bistable region are studied. The critical triggering values are determined with the changes of the non-dimensional heater power, the damping coefficient and the relative heater location. The critical triggering value of velocity perturbation decreases with the increase of non-dimensional heater power, whereas an increasing trend is observed with the increase of damping coefficient. Interestingly, the critical triggering value firstly decreases and then increases with the increase of the relative heater location. The variation of critical triggering value for pressure perturbation is found to correspond with velocity perturbation. In the bistable region, the amplitude and frequency of the steady limit cycle oscillation of the system are independent of the initial perturbation values, but the perturbation value has strong effect on the duration needed to achieve the steady limit cycle, and the time required for the system to reach the limit cycle under the perturbation of U1=0.4 is about 3 times longer than that of U1=0.8.