Physics-informed line graph neural network for power flow calculation
Power flow calculation plays a significant role in the operation and planning of modern power systems. Traditional numerical calculation methods have good interpretability but high time complexity. They are unable to cope with increasing amounts of data in power systems; therefore, many machine learning based methods have been proposed for more efficient power flow calculation. Despite the good performance of these methods in terms of computation speed, they often overlook the importance of transmission lines and do not fully consider the physical mechanisms in the power systems, thereby weakening the prediction accuracy of power flow. Given the importance of the transmission lines as well as to comprehensively consider their mutual influence, we shift our focus from bus adjacency relationships to transmission line adjacency relationships and propose a physics-informed line graph neural network framework. This framework propagates information between buses and transmission lines by introducing the concepts of the incidence matrix and the line graph matrix. Based on the mechanics of the power flow equations, we further design a loss function by integrating physical information to ensure that the output results of the model satisfy the laws of physics and have better interpretability. Experimental results on different power grid datasets and different scenarios demonstrate the accuracy of our proposed model.
Decoding aging and cognitive functioning through spatiotemporal EEG patterns: Introducing spatiotemporal information-based similarity analysis
Exploring spatiotemporal patterns of high-dimensional electroencephalography (EEG) time series generated from complex brain system is crucial for deciphering aging and cognitive functioning. Analyzing high-dimensional EEG series poses challenges, particularly when employing distance-based methods for spatiotemporal dynamics. Therefore, we proposed an innovative methodology for multi-channel EEG data, termed as Spatiotemporal Information-based Similarity (STIBS) analysis. The core of this method is to first perform state space compression of multi-channel EEG time series using global field power, which can provide insight into the dynamic integration of spatiotemporal patterns between the steady states and non-steady states of brain. Subsequently, we quantify the pairwise differences and non-randomness of spatiotemporal patterns using an information-based similarity analysis. Results demonstrated that this method holds the potential to serve as a distinguishing marker between young and elderly on both pairwise differences and non-randomness indices. Young individuals and those with higher cognitive abilities exhibit more complex macrostructure and non-random spatiotemporal patterns, whereas both aging and cognitive decline lead to more randomized spatiotemporal patterns. We further extended the proposed analytics to brain regions adversarial STIBS (bra-STIBS), highlighting differences between young and elderly, as well as high and low cognitive groups. Furthermore, utilizing the STIBS-based XGBoost model yields superior recognition accuracy in aging (93.05%) and cognitive functioning (74.29%, 64.19%, and 80.28%, respectively, for attention, memory, and compatibility performance recognition). STIBS-based methodology not only contributes to the ongoing exploration of neurobiological changes in aging but also provides a powerful tool for characterizing the spatiotemporal nonlinear dynamics of the brain and their implications for cognitive functioning.
Orbits of a system of three point vortices and the associated chaotic mixing
We study the general periodic motion of a set of three point vortices in the plane, as well as the potentially chaotic motion of one or more tracer particles. While the motion of three vortices is simple in that it can only be periodic, the actual orbits can be surprisingly complex and varied. This rich behavior arises from the existence of both co-linear and equilateral relative equilibria (steady motion in a rotating frame of reference). Here, we start from a general (unsteady) co-linear array with arbitrary vortex circulations. The subsequent motion may take the vortices close to a distinct co-linear relative equilibrium or to an equilateral one. Both equilibrium states are necessarily unstable, as we demonstrate by a linear stability analysis. We go on to study mixing by examining Poincaré sections and finite-time Lyapunov exponents. Both indicate widespread chaotic motion in general, implying that the motion of three vortices efficiently mixes the nearby surrounding fluid outside of small regions surrounding each vortex.
An investigation of escape and scaling properties of a billiard system
We investigate some statistical properties of escaping particles in a billiard system whose boundary is described by two control parameters with a hole on its boundary. Initially, we analyze the survival probability for different hole positions and sizes. We notice that the survival probability follows an exponential decay with a characteristic power-law tail when the hole is positioned partially or entirely over large stability islands in phase space. We find that the survival probability exhibits scaling invariance with respect to the hole size. In contrast, the survival probability for holes placed in predominantly chaotic regions deviates from the exponential decay. We introduce two holes simultaneously and investigate the complexity of the escape basins for different hole sizes and control parameters by means of the basin entropy and the basin boundary entropy. We find a non-trivial relation between these entropies and the system's parameters and show that the basin entropy exhibits scaling invariance for a specific control parameter interval.
Mixed-mode oscillations and chaos in a complex chemical reaction network involving heterogeneous catalysis
The nonlinear dynamical behavior in a complex isothermal reaction network involving heterogeneous catalysis is studied. The method first determines the multiple steady states in the reaction network. This is followed by an analysis of bifurcation continuations to identify several kinds of bifurcations, including limit point, Bogdanov-Takens, generalized Hopf, period doubling, and generalized period doubling. Numerical simulations are performed around the period doubling and generalized period doubling bifurcations. Rich nonlinear behaviors are observed, including simple sustained oscillations, mixed-mode oscillations, non-mixed-mode chaotic oscillations, and mixed-mode chaotic oscillations. Concentration-time plots, 2D phase portraits, Poincaré maps, maximum Lyapunov exponents, frequency spectra, and cascade of bifurcations are reported. Period-doubling and period-adding routes leading to chaos are observed. Maximum Lyapunov exponents are positive for all the chaotic cases, but they are also positive for some non-chaotic orbits. This result diminishes the reliability of using maximum Lyapunov exponents as a tool for determining chaos in the network under study.
A novel discrete memristive hyperchaotic map with multi-layer differentiation, multi-amplitude modulation, and multi-offset boosting
In recent years, the introduction of memristors in discrete chaotic map has attracted much attention due to its enhancement of the complexity and controllability of chaotic maps, especially in the fields of secure communication and random number generation, which have shown promising applications. In this work, a three-dimensional discrete memristive hyperchaotic map (3D-DMCHM) based on cosine memristor is constructed. First, we analyze the fixed points of the map and their stability, showing that the map can either have a linear fixed point or none at all, and the stability depends on the parameters and initial state of the map. Then, phase diagrams, bifurcation diagrams, Lyapunov exponents, timing diagrams, and attractor basins are used to analyze the complex dynamical behaviors of the 3D-DMCHM, revealing that the 3D-DMCHM enters into a chaotic state through a period-doubling bifurcation path, and some special dynamical phenomena such as multi-layer differentiation, multi-amplitude control, and offset boosting behaviors are also observed. In particular, with the change of memristor initial conditions, there exists an offset that only homogeneous hidden chaotic attractors or a mixed state offset with coexistence of point attractors and chaotic attractors. Finally, we confirmed the high complexity of 3D-DMCHM through complexity tests and successfully implemented it using a digital signal processing circuit, demonstrating its hardware feasibility.
Chaotic dynamics in a class of generalized memristive maps
The memory effects of the memristors in nonlinear systems make the systems generate complicated dynamics, which inspires the development of the applications of memristors. In this article, the model of the discrete memristive systems with the generalized Ohm's law is introduced, where the classical Ohm's law is a linear relationship between voltage and current, and a generalized Ohm's law is a nonlinear relationship. To illustrate the rich dynamics of this model, the complicated dynamical behavior of three types of maps with three types of discrete memristances is investigated, where a cubic function representing a kind of generalized Ohm's law is used, and this cubic function is a simplified characteristic of the famous tunnel diode. The existence of attractors with one or two positive Lyapunov exponents (corresponding to chaotic or hyperchaotic dynamics) is obtained, and the coexistence of (infinitely) many attractors is observable. A hardware device is constructed to implement these maps and the analog voltage signals are experimentally acquired.
Templex-based dynamical units for a taxonomy of chaos
Discriminating different types of chaos is still a very challenging topic, even for dissipative three-dimensional systems for which the most advanced tool is the template. Nevertheless, getting a template is, by definition, limited to three-dimensional objects based on knot theory. To deal with higher-dimensional chaos, we recently introduced the templex combining a flow-oriented BraMAH cell complex and a directed graph (a digraph). There is no dimensional limitation in the concept of templex. Here, we show that a templex can be automatically reduced into a "minimal" form to provide a comprehensive and synthetic view of the main properties of chaotic attractors. This reduction allows for the development of a taxonomy of chaos in terms of two elementary units: the oscillating unit (O-unit) and the switching unit (S-unit). We apply this approach to various well-known attractors (Rössler, Lorenz, and Burke-Shaw) as well as a non-trivial four-dimensional attractor. A case of toroidal chaos (Deng) is also treated.
Evolution of cooperation in heterogeneous populations with asymmetric payoff distribution
The emergence and maintenance of cooperation is a complex and intriguing issue, especially in the context of widespread asymmetries in interactions that arise from individual differences in real-world scenarios. This study investigates how asymmetric payoff distribution affects cooperation in public goods games by considering a population composed of two types of individuals: strong and weak. The asymmetry is reflected in the fact that strong players receive a larger share of the public pool compared to weak players. Our results demonstrate that asymmetric payoff distribution can promote cooperation in well-mixed populations and trigger the co-evolution of cooperation between sub-populations of strong and weak players. In structured populations, however, the effect of asymmetric payoff distribution on cooperation is contingent on the proportion of strong players and the extent of their payoff share, which can either foster or inhibit cooperation. By adjusting the interaction probability between strong and weak players based on their spatial arrangement on lattice networks, we find that moderate interaction probabilities most effectively maintain cooperation. This study provides valuable insights into the dynamics of cooperation under asymmetric conditions, highlighting the complex role of asymmetrical interactions in the evolution of cooperation.
Dark gap solitons in bichromatic optical superlattices under cubic-quintic nonlinearities
We demonstrate the existence of two types of dark gap solitary waves-the dark gap solitons and the dark gap soliton clusters-in Bose-Einstein condensates trapped in a bichromatic optical superlattice with cubic-quintic nonlinearities. The background of these dark soliton families is different from the one in a common monochromatic linear lattice; namely, the background in our model is composed of two types of Gaussian-like pulses, whereas in the monochromatic linear lattice, it is composed of only one type of Gaussian-like pulses. Such a special background of dark soliton families is convenient for the manipulation of solitons by the parameters of bichromatic and chemical potentials. The dark soliton families in the first, second, and third bandgap in our model are studied. Their stability is assessed by the linear-stability analysis, and stable as well as unstable propagation of these gap solitons are displayed. The profiles, stability, and perturbed evolution of both types of dark soliton families are distinctly presented in this work.
Tunable disorder on the S-state majority-voter model
We investigate the nonequilibrium phase transition in the S-state majority-vote model for S=2,3, and 4. Each site, k, is characterized by a distinct noise threshold, qk, which indicates its resistance to adopting the majority state of its Nv nearest neighbors. Precisely, this noise threshold is governed by a hyperbolic distribution, P(k)∼1/k, bounded within the limits e-α/2
Model-free distributed state estimation with local measurements
In this paper, the state estimation problem of physical plants with unknown system dynamic is revisited from the perspective of limited output information measurement, which corresponds to those with characteristics of high-dimensional, wide-area coverage and scatter. Given this fact, a network of sensors are used to carry out the measurement with each one accessing only partial outputs of the targeted systems and a novel model-free state estimation approach, named distributed stochastic variational inference state estimation, is proposed. The key idea of this method is to compensate for the impacts of local output measurements by adding nearest-neighbor rule-based information interaction among estimators to complete the state estimation. It finds from the numerical experiments that the proposed method has clear advantages in both estimation accuracy and speed, and it also provides guidance on how to improve the efficiency of state estimation under local measurements.
Shrinking shrimp-shaped domains and multistability in the dissipative asymmetric kicked rotor map
An interesting feature in dissipative nonlinear systems is the emergence of characteristic domains in parameter space that exhibit periodic temporal evolution, known as shrimp-shaped domains. We investigate the parameter space of the dissipative asymmetric kicked rotor map and show that, in the regime of strong dissipation, the shrimp-shaped domains repeat themselves as the nonlinearity parameter increases while maintaining the same period. We analyze the dependence of the length of each periodic domain with the nonlinearity parameter, revealing that it follows a power law with the same exponent regardless of the dissipation parameter. Additionally, we find that the distance between adjacent shrimp-shaped domains is scaling invariant with respect to the dissipation parameter. Furthermore, we show that for weaker dissipation, a multistable scenario emerges within the periodic domains. We find that as the dissipation gets weaker, the ratio of multistable parameters for each periodic domain increases, and the area of the periodic basin decreases as the nonlinearity parameter increases.
Evolutionary dynamics of stochastic games in set-structured populations
In structured populations, the ecology of games may vary over neighborhoods. The effect of the ecological variations on population dynamics remains largely unknown. We here incorporate the ecological variations into the set-structured populations to explore the coevolutionary dynamics of the ecology and cooperation. Individuals of a population are distributed over sets. Interactions occur in the form of evolutionary games. When two individuals share more common sets, they play the weak prisoner's dilemma. Otherwise, they play the strong prisoner's dilemma. Both the set memberships and the strategy update in the evolutionary process. Changes in set memberships hold sway over the games to be played, which, in turn, influences the performance of strategies. Combining evolutionary set theory and random walks on graphs, we derived the conditions for cooperation to be selected under the weak selection limit. We find that a denser set-structured population increases the probability of individuals participating in a weak prisoner's dilemma, and thereby promoting the spread of cooperation. Properly modulating the population structure and the payoff feedback can further lower the critical benefit-cost ratio required for cooperation to be selected. Our results may help better understand the effects of ecological variations in enhancing cooperative behavior in set-structured populations.
Hidden multi-scroll and coexisting self-excited attractors in optical injection semiconductor laser system: Its electronic control
In this paper, we investigate hidden and coexisting self-excited multi-scroll attractors by using a modified rate equations model of semiconductor lasers (REM-SCLs) subjected to optical injection by exploring various quantifying analytical and numerical methods. The multi-leveled dynamics sticks out the existence of several sets of equilibria that asymptotically attract trajectories originating outside of them. Chaos topology based on the impact of equilibria allows the describing of the so-called stable or unstable multi-scroll chaotic attractors. Shaping of the new coexisting self-excited multi-scroll attractor, whose source is from coupling of equilibria, is analyzed, as well as its structural dynamics along with the dynamical emergence of the hidden multi-scroll attractor in the restricted interval, defined by an additional decisive parameter. Additionally, specific 3D plots with embedded contour plots obtained by harnessing two-parameter bifurcation analysis clarify structural dynamics of such a multi-scroll attractor and accurately circumscribe stretching of its fractal-like basin of attraction. Strange metamorphoses undergone by the fractal-like basin of attraction of the studied multi-scroll attractor are stepwisely parsed in the map of two-codimension bifurcation as its scroll number evolves. At last, an electronic circuit of equivalent REM-SCLs is designed and simulated in the PSpice environment alongside a tailored electronic controller. The achieved results align with the ones of numerical analysis; besides, temporal controlling of optical waves pertaining thereto is also fulfilled.
Chaotic dynamics and optimal therapeutic strategies for Caputo fractional tumor immune model in combination therapy
In this paper, a Caputo fractional tumor immune model of combination therapy is established. First, the stability and biological significance of each equilibrium point are analyzed, and it is demonstrated that chaos may arise under specific conditions. Combined with the mathematical definition of Caputo fractional differentiation (CFD), it is found that there is a high correlation between the chaotic phenomenon of the patient's condition and the sensitivity of the patient to the change in the state of the day. The bifurcation threshold of each parameter is determined through numerical simulation, and the Hopf bifurcation of direct competition coefficient and inhibition coefficient between tumor cells and host healthy cells is elaborated upon in detail. Subsequently, a novel method combining optimal control theory with the particle swarm optimization (PSO) algorithm is proposed for the optimal control of the tumor immune model in combination therapy. Finally, the Adams-Bashforth-Moulton (ABM) prediction correction method is utilized in numerical simulations which demonstrate that the introduction of the CFD alters the model dynamics. Furthermore, these results indicate that fractional calculus can effectively be applied to tumor immune models better to elucidate complex chaotic dynamics of tumor cell evolution. Concurrently, the PSO can be successfully integrated with optimal control theory to address optimization challenges in cancer treatment.
Co-evolutionary dynamics for two adaptively coupled Theta neurons
Natural and technological networks exhibit dynamics that can lead to complex cooperative behaviors, such as synchronization in coupled oscillators and rhythmic activity in neuronal networks. Understanding these collective dynamics is crucial for deciphering a range of phenomena from brain activity to power grid stability. Recent interest in co-evolutionary networks has highlighted the intricate interplay between dynamics on and of the network with mixed time scales. Here, we explore the collective behavior of excitable oscillators in a simple network of two Theta neurons with adaptive coupling without self-interaction. Through a combination of bifurcation analysis and numerical simulations, we seek to understand how the level of adaptivity in the coupling strength, a, influences the dynamics. We first investigate the dynamics possible in the non-adaptive limit; our bifurcation analysis reveals stability regions of quiescence and spiking behaviors, where the spiking frequencies mode-lock in a variety of configurations. Second, as we increase the adaptivity a, we observe a widening of the associated Arnol'd tongues, which may overlap and give room for multi-stable configurations. For larger adaptivity, the mode-locked regions may further undergo a period-doubling cascade into chaos. Our findings contribute to the mathematical theory of adaptive networks and offer insights into the potential mechanisms underlying neuronal communication and synchronization.
Two-parameter dynamics and multistability of a non-smooth railway wheelset system with dry friction damping
A deep understanding of non-smooth dynamics of vehicle systems, particularly with dry friction damping offer valuable insights into the design and optimization of railway vehicle systems, ultimately enhancing the safety and reliability of railway operations. In this paper, the two-parameter dynamics of a non-smooth railway wheelset system incorporating dry friction damping are investigated. The effect of the crucial parameters on the complexity of the evolution process is comprehensively exposed by identifying different dynamic responses in the two-parameter plane. In addition, the multistability and the various routes transition to chaos for the system are also discussed. It is found that dry friction induces highly complex dynamics in the system, encompassing a range of behaviors such as periodic, quasi-periodic, and chaotic motions. These intricate dynamics are a direct result of the interplay between multiple parameters, such as speed and damping coefficients, which are critical in determining the system's stability and performance. The presence of multistability further complicates the system, resulting in unpredictable transitions between different motion states.
Linear and nonlinear causality in financial markets
Identifying and quantifying co-dependence between financial instruments is a key challenge for researchers and practitioners in the financial industry. Linear measures such as the Pearson correlation are still widely used today, although their limited explanatory power is well known. In this paper, we present a much more general framework for assessing co-dependencies by identifying linear and nonlinear causalities in the complex system of financial markets. To do so, we use two different causal inference methods, transfer entropy and convergent cross-mapping, and employ Fourier transform surrogates to separate their linear and nonlinear contributions. We find that stock indices in Germany and the U.S. exhibit a significant degree of nonlinear causality and that correlation, while a very good proxy for linear causality, disregards nonlinear effects and hence underestimates causality itself. The presented framework enables the measurement of nonlinear causality, the correlation-causality fallacy, and motivates how causality can be used for inferring market signals, pair trading, and risk management of portfolios. Our results suggest that linear and nonlinear causality can be used as early warning indicators of abnormal market behavior, allowing for better trading strategies and risk management.
Inverse stochastic resonance in adaptive small-world neural networks
Inverse stochastic resonance (ISR) is a counterintuitive phenomenon where noise reduces the oscillation frequency of an oscillator to a minimum occurring at an intermediate noise intensity, and sometimes even to the complete absence of oscillations. In neuroscience, ISR was first experimentally verified with cerebellar Purkinje neurons [Buchin et al., PLOS Comput. Biol. 12, e1005000 (2016)]. These experiments showed that ISR enables a locally optimal information transfer between the input and output spike train of neurons. Subsequent studies have further demonstrated the efficiency of information processing and transfer in neural networks with small-world network topology. We have conducted a numerical investigation into the impact of adaptivity on ISR in a small-world network of noisy FitzHugh-Nagumo (FHN) neurons, operating in a bi-metastable regime consisting of a metastable fixed point and a metastable limit cycle. Our results show that the degree of ISR is highly dependent on the value of the FHN model's timescale separation parameter ε. The network structure undergoes dynamic adaptation via mechanisms of either spike-time-dependent plasticity (STDP) with potentiation-/depression-domination parameter P or homeostatic structural plasticity (HSP) with rewiring frequency F. We demonstrate that both STDP and HSP amplify the effect of ISR when ε lies within the bi-stability region of FHN neurons. Specifically, at larger values of ε within the bi-stability regime, higher rewiring frequencies F are observed to enhance ISR at intermediate (weak) synaptic noise intensities, while values of P consistent with depression-domination (potentiation-domination) consistently enhance (deteriorate) ISR. Moreover, although STDP and HSP control parameters may jointly enhance ISR, P has a greater impact on improving ISR compared to F. Our findings inform future ISR enhancement strategies in noisy artificial neural circuits, aiming to optimize local information transfer between input and output spike trains in neuromorphic systems and prompt venues for experiments in neural networks.
Impact of elevated carbon dioxide and temperature on crop yield: A study of autonomous and nonautonomous systems
Agricultural crops are vital for human survival, forming the backbone of global food supply. However, the rising atmospheric carbon dioxide (CO2) level and the increasing temperature relative to the pre-industrial level are poised to impact the yields of essential staple food crops significantly. In this research, we propose and analyze a nonlinear mathematical model to investigate the effects of elevated CO2 and temperature on crop yield. Our model assumes that rising CO2 levels elevate the global average temperature, and the surface temperature initially boosts the growth rate of crops until a threshold is reached, after which the growth rate declines. We also incorporate seasonal variations into the model and perform a comprehensive analytical and numerical analysis of both the autonomous and associated nonautonomous systems. Our findings reveal a critical threshold for anthropogenic CO2 emissions, beyond which the crop yield starts to decrease. Notably, crops with high-temperature tolerance demonstrate higher yields even under elevated CO2 conditions, suggesting a viable strategy for mitigating climate change impacts: developing or utilizing crop varieties with enhanced temperature tolerance. Moreover, our analysis of the nonautonomous system uncovers periodic solutions when the corresponding autonomous system is stable. The nonautonomous system also exhibits complex dynamics, including higher-period oscillations and chaos, when the autonomous system undergoes limit-cycle oscillations. This study provides valuable insights into the interplay between CO2 level, global average surface temperature, and crop yield, offering potential strategies for safeguarding agricultural productivity in the face of climate change.
