2025-06-16 | | Total: 7
Self-organising systems demonstrate how simple local rules can generate complex stochastic patterns. Many natural systems rely on such dynamics, making self-organisation central to understanding natural complexity. A fundamental challenge in modelling such systems is solving the inverse problem: finding the unknown causal parameters from macroscopic observations. This task becomes particularly difficult when observations have a strong stochastic component, yielding diverse yet equivalent patterns. Traditional inverse methods fail in this setting, as pixel-wise metrics cannot capture feature similarities between variable outcomes. In this work, we introduce a novel inverse modelling method specifically designed to handle stochasticity in the observable space, leveraging the capacity of visual embeddings to produce robust representations that capture perceptual invariances. By mapping the pattern representations onto an invariant embedding space, we can effectively recover unknown causal parameters without the need for handcrafted objective functions or heuristics. We evaluate the method on two canonical models--a reaction-diffusion system and an agent-based model of social segregation--and show that it reliably recovers parameters despite stochasticity in the outcomes. We further apply the method to real biological patterns, highlighting its potential as a tool for both theorists and experimentalists to investigate the dynamics underlying complex stochastic pattern formation.
Advancements in Noisy Intermediate-Scale Quantum (NISQ) computing are steadily pushing these systems toward outperforming classical supercomputers on specific, well-defined computational tasks. In this work, we explore and control quantum chaos in NISQ systems using discrete-time quantum walks (DTQW) on cyclic graphs. To efficiently implement quantum walks on NISQ hardware, we employ the quantum Fourier transform (QFT) to diagonalize the conditional shift operator, optimizing circuit depth and fidelity. We experimentally realize the transition from quantum chaos to order via DTQW dynamics on both odd and even cyclic graphs, specifically 3- and 4-cycle graphs, using the counterintuitive Parrondo's paradox strategy across three different NISQ devices. While the 4-cycle graphs exhibit high-fidelity quantum evolution, the 3-cycle implementation shows significant fidelity improvement when augmented with dynamical decoupling pulses. Our results demonstrate a practical approach to probing and harnessing controlled chaotic dynamics on real quantum hardware, laying the groundwork for future quantum algorithms and cryptographic protocols based on quantum walks.
We present a coarse-grained stochastic model for axonal extension on periodic arrays of parallel micropatterns that integrates three key biophysical mechanisms: (i) the molecular clutch that couples actin retrograde flow to substrate adhesions, (ii) an active biopolymer-based mechanism generating traction-force fluctuations, and (iii) the mechanical interaction of the growth cone with the micropatterned substrate. Using the Shannon-Jaynes maximum entropy principle with constraints derived from experimental observations, we derive a unique probability distribution for the colored acceleration noise that enters the Langevin equation. The resulting stationary process exhibits power-law temporal correlations with negative exponent, which accounts for the observed superdiffusive dynamics of axons. For biologically relevant parameters the model predicts this exponent to be -1/2, in close quantitative agreement with measurements of cortical neurons cultured on patterned substrates.
We analyse a generalisation of Zöttl and Stark's model of active spherical particles [Phys. Rev. Lett. 108, 218104 (2012)] and prolate spheroidal particles [Eur. Phys. J. E 36(1), 4 (2013)] suspended in cylindrical Poiseuille flow, to particle dynamics in an arbitrary unidirectional steady laminar flow through a straight duct geometry. Our primary contribution is to describe a Hamiltonian formulation of these systems and provide explicit forms of the constants of motions in terms of the arbitrary fluid velocity field. The Hamiltonian formulation provides a convenient and robust approach to the computation of particle orbits whilst also providing new insights into the dynamics, specifically the way in which orbits are trapped within basins defined by a potential well. In addition to considering spherical and prolate spheroidal particles, we also illustrate that the model can be adapted to oblate spheroidal particles.
While synchronized states, and the dynamical pathways through which they emerge, are often regarded as the paradigm to understand the dynamics of information spreading on undirected networks of nonlinear dynamical systems, when we consider directed network architectures, dynamical stationary states can arise. To study this phenomenon we consider the simplest directed network, a single cycle, and excitable FitzHugh-Nagumo (FHN) neurons. We show numerically that a stationary dynamical state emerges in the form of a self-sustained traveling wave, through a saddle-point bifurcation of limit cycles that does not destabilize the global fixed point of the system. We then formulate an effective model for the dynamical steady state of the cycle in terms of a single-neuron Delay Differential Equation (DDE) featuring an explicitly delayed feedback, demonstrating numerically the possibility of mapping stationary solutions between the two models. The DDE based model is shown to reproduce the entire bifurcation, which also in this case does not destabilize the global fixed point, even though global properties differ in general between the systems. The discrete nature of the cycle graph is revealed as the origin of these coordinated states by the parametric analysis of solutions, and the DDE effective model is shown to preserve this feature accurately. Finally, the scaling of the inter-site propagation times hints to a solitonic nature of the wave state in the limit of large chain size.
Tipping points (TP) are abrupt transitions between metastable states in complex systems, most often described by a bifurcation or crisis of a multistable system induced by a slowly changing control parameter. An avenue for predicting TPs in real-world systems is critical slowing down (CSD), which is a decrease in the relaxation rate after perturbations prior to a TP that can be measured by statistical early warning signals (EWS) in the autocovariance of observational time series. In high-dimensional systems, we cannot expect a priori chosen scalar observables to show significant EWS, and some may even show an opposite signal. Thus, to avoid false negative or positive early warnings, it is desirable to monitor fluctuations only in observables that are designed to capture CSD. Here we propose that a natural observable for this purpose can be obtained by a data-driven approximation of the first non-trivial eigenfunction of the backward Fokker-Planck (or Kolmogorov) operator, using the diffusion map algorithm.
Periodic signals propagating along chains are common in biology, for example in locomotion and peristalsis, and are also of interest for continuum robots. In previous work we constructed such networks as 'feedforward lifts' of a central pattern generator (CPG). When the CPG undergoes periodic oscillations, created by Hopf bifurcation or other mechanisms, it can then transmit periodic signals along one or more feedforward chains in a synchronous or phase-synchronous manner. We proved necessary and sufficient conditions for the stability of these lifted periodic orbits, in several senses. Here we examine the implications of the resulting theory for chains of neurons, using several standard neuron models: FitzHugh-Nagumo, Morris-Lecar, Hindmarsh-Rose, and Hodgkin-Huxley. We compare different notions of transverse stability, and summarize some numerical simulations showing that for all these neuron models the propagating signal can be transversely Floquet stable. Finally we discuss implications for less idealized models.