No Arabic abstract
Estimation of parameters in differential equation models can be achieved by applying learning algorithms to quantitative time-series data. However, sometimes it is only possible to measure qualitative changes of a system in response to a controlled condition. In dynamical systems theory, such change points are known as bifurcations and lie on a function of the controlled condition called the bifurcation diagram. In this work, we propose a gradient-based semi-supervised approach for inferring the parameters of differential equations that produce a user-specified bifurcation diagram. The cost function contains a supervised error term that is minimal when the model bifurcations match the specified targets and an unsupervised bifurcation measure which has gradients that push optimisers towards bifurcating parameter regimes. The gradients can be computed without the need to differentiate through the operations of the solver that was used to compute the diagram. We demonstrate parameter inference with minimal models which explore the space of saddle-node and pitchfork diagrams and the genetic toggle switch from synthetic biology. Furthermore, the cost landscape allows us to organise models in terms of topological and geometric equivalence.
Asymptotic state of an open quantum system can undergo qualitative changes upon small variation of system parameters. We demonstrate it that such quantum bifurcations can be appropriately defined and made visible as changes in the structure of the asymptotic density matrix. By using an $N$-boson open quantum dimer, we present quantum diagrams for the pitchfork and saddle-node bifurcations in the stationary case and visualize a period-doubling transition to chaos for the periodically modulated dimer. In the latter case, we also identify a specific bifurcation of purely quantum nature.
Empirical Dynamic Modeling (EDM) is a nonlinear time series causal inference framework. The latest implementation of EDM, cppEDM, has only been used for small datasets due to computational cost. With the growth of data collection capabilities, there is a great need to identify causal relationships in large datasets. We present mpEDM, a parallel distributed implementation of EDM optimized for modern GPU-centric supercomputers. We improve the original algorithm to reduce redundant computation and optimize the implementation to fully utilize hardware resources such as GPUs and SIMD units. As a use case, we run mpEDM on AI Bridging Cloud Infrastructure (ABCI) using datasets of an entire animal brain sampled at single neuron resolution to identify dynamical causation patterns across the brain. mpEDM is 1,530 X faster than cppEDM and a dataset containing 101,729 neuron was analyzed in 199 seconds on 512 nodes. This is the largest EDM causal inference achieved to date.
Motivated by questions originating from the study of a class of shallow student-teacher neural networks, methods are developed for the analysis of spurious minima in classes of gradient equivariant dynamics related to neural nets. In the symmetric case, methods depend on the generic equivariant bifurcation theory of irreducible representations of the symmetric group on $n$ symbols, $S_n$; in particular, the standard representation of $S_n$. It is shown that spurious minima do not arise from spontaneous symmetry breaking but rather through a complex deformation of the landscape geometry that can be encoded by a generic $S_n$-equivariant bifurcation. We describe minimal models for forced symmetry breaking that give a lower bound on the dynamic complexity involved in the creation of spurious minima when there is no symmetry. Results on generic bifurcation when there are quadratic equivariants are also proved; this work extends and clarifies results of Ihrig & Golubitsky and Chossat, Lauterback & Melbourne on the instability of solutions when there are quadratic equivariants.
An epidemic model with distributed time delay is derived to describe the dynamics of infectious diseases with varying immunity. It is shown that solutions are always positive, and the model has at most two steady states: disease-free and endemic. It is proved that the disease-free equilibrium is locally and globally asymptotically stable. When an endemic equilibrium exists, it is possible to analytically prove its local and global stability using Lyapunov functionals. Bifurcation analysis is performed using DDE-BIFTOOL and traceDDE to investigate different dynamical regimes in the model using numerical continuation for different values of system parameters and different integral kernels.
The stability of functional differential equations under delayed feedback is investigated near a Hopf bifurcation. Necessary and sufficient conditions are derived for the stability of the equilibrium solution using averaging theory. The results are used to compare delayed versus undelayed feedback, as well as discrete versus distributed delays. Conditions are obtained for which delayed feedback with partial state information can yield stability where undelayed feedback is ineffective. Furthermore, it is shown that if the feedback is stabilizing (respectively, destabilizing), then a discrete delay is locally the most stabilizing (resp., destabilizing) one among delay distributions having the same mean. The result also holds globally if one considers delays that are symmetrically distributed about their mean.