2026-07-17 | | Total: 18
Here we study the core entropy function $h: PM(d) \rightarrow \mathbb{R}$ defined on the space of degree $d$ primitive majors, a combinatorial model for the space of degree $d$ polynomials. In particular, we introduce techniques to classify the global maximum for each $d$ and the maxima along certain strata. We show that there are $d-1$ global maxima over all of $PM(d)$ and that the maxima on the unicritical stratum is a Cantor set.
Given a scalar observable of an ergodic dynamical system with a low-dimensional attractor, two families of methods reconstruct and predict the underlying state: recurrence-based methods (the method of analogues and its descendants), which wait for the trajectory to return to an $\varepsilon$-neighborhood of a previously observed state, and observer-based methods, which fit a converging state estimator on the delay reconstruction. We formalize and empirically verify an exponential separation between the two: the expected cost of recurrence scales as $\varepsilon^{-d}$, where $d$ is the pointwise dimension of the invariant measure (a consequence of the Kac lemma and quantitative Poincare recurrence), whereas a detectable linear observer converges in $Θ(\log(1/\varepsilon)/(1-ρ(A_{cl})^2))$ steps, where $ρ(A_{cl})$ is the closed-loop spectral radius of the Riccati fixed point. Both laws are verified numerically (return-time exponent $-1.8$ on the Lorenz attractor against the theoretical $-2.05$; observer cost linear in $\log(1/\varepsilon)$ with $R^2=1.000$ and in $(1-ρ^2)^{-1}$ with $R^2=0.985$), yielding a measured cost gap of $\sim 10^{9}$ at $\varepsilon=10^{-6}$ for $d\approx 2$. We complement the theorem with an admission protocol (the Kac-Riccati gate) deciding whether a signal lies inside the theorem's class, via surrogate-data prediction gating; it also explains the folklore of "universal" fractal dimensions as a dataset-size artifact bounded by $2\log_{10}N$. On real data the gate admits the Santa Fe laser benchmark ($\hat D_2=2.0$) and refuses the monthly sunspot series, reproducing the settled resolution of historical low-dimensionality claims. All results reproduce from a single verification script (17/17 checks).
We consider a planar ODE system from an important class of models for gene regulatory dynamics. The system depends singularly on the steepness parameters $0<\varepsilon_1,\varepsilon_2 \ll 1$ and converges to a piecewise-smooth system as these parameters tend to zero. Unlike previous studies, we do not assume that $\varepsilon_1 = \varepsilon_2$. As a consequence, the dynamics when $(\varepsilon_1, \varepsilon_2) \to (0,0)$ depends upon how the limit is taken. Using a preliminary blow-up in parameter space, we identify three distinct singular limits. We perform a two-parameter bifurcation analysis in each case, and apply multiple geometric blow-ups in variable and parameter space to determine the bifurcation structure and the associated global dynamics. Bogdanov-Takens bifurcations are revealed in two of three cases, and in one case in particular, a regularised visible-invisible two-fold singularity is shown to organise the unfolding of singular bifurcations in the vicinity of canards. Our results show that the qualitative dynamics and overall sensitivity of the system to parameter variation depends on the relative size of the steepness parameters. More generally, the analytical framework developed herein provides a systematic approach to singular perturbation problems with multiple independent small parameters that should apply well beyond gene regulatory network models.
When $X$ is a locally compact Euclidean building, the isometry group of $X$ acts by homeomorphisms on the space of $k$-simplices in the visual boundary of $X$. We consider perturbations of these actions for discrete groups of isometries acting with compact quotient on $X$, showing that all small enough perturbations are semi-conjugate to the original action. This proves in particular that, when $Q$ is any parabolic subgroup in a semisimple $p$-adic Lie group $G$, the induced action of a cocompact lattice in $G$ has a topologically stable action on $G/Q$.
Contemporary machine learning is founded almost entirely on Archimedean mathematics: data are embedded in real vector spaces, similarity is measured by Euclidean-type metrics, learning is formulated through real-valued loss functions, and optimization is driven by differential or gradient-based methods. Recent work has shown that p-adic and ultrametric methods can enrich neural architectures, hierarchical representation, classification, and information processing. Rather than constructing a p-adic analogue of existing neural networks, we advocate a more fundamental viewpoint. We propose that a genuinely non-Archimedean theory of learning should be derived intrinsically from the topology, geometry, algebra, and measure structure of ultrametric spaces. To this end, we formulate the Principle of Intrinsic Non-Archimedean Learning and introduce an axiomatic framework in which information states are organized by nested ultrametric balls, learning is represented by hierarchical refinement and redistribution of information, and scale is determined by valuation depth. In this framework, neurons, layers, activation functions, gradients, and backpropagation are not assumed to be primitive concepts but possible computational realizations of a deeper mathematical theory. We prove a basic structural result showing that, on \(\Zp\) refinement depth coincides with common-prefix depth and that the finite quotients \(\mathbb Z/p^N \mathbb Z \) preserve the rooted ball hierarchy through level N. Consequently, hierarchy is not an auxiliary structure to be learned but is already encoded in the intrinsic geometry of the underlying information space
As a continuation of our previous work, where bifurcations of periodic and antiperiodic orbits near an equilibrium were studied for autonomous differential delay systems with one or two delays, this paper investigates bifurcations of periodic and antiperiodic solutions for five classes of non-autonomous potential delay differential systems with multiple delays. By reformulating these systems as generalized Hamiltonian systems on Euclidean spaces, we apply the Hamiltonian bifurcation theory recently developed by the author to establish alternative bifurcation results of both Fadell--Rabinowitz and Rabinowitz type.
We investigate bifurcations of periodic and antiperiodic solutions near equilibria for four classes of potential differential-delay equations involving one or two delays. By reformulating each system as a (generalized) Hamiltonian system, and applying the Hamiltonian bifurcation theory recently developed by the author we obtain bifurcation results of both Fadell--Rabinowitz and Rabinowitz type. The analysis relies on the computation of Maslov--type indices for fundamental solutions of the associated linear Hamiltonian systems.
We prove the analogue of Belinskaya's theorem for measure-preserving flows: two free ergodic measure-preserving flows whose $\mathrm{L}^1$ full groups are isomorphic as abstract groups are conjugate up to a scalar time change. This answers a question posed by François Le Maître and the author. We show that whenever two free ergodic flows generate the same orbit equivalence relation and one is contained in the other's $\mathrm{L}^1$ full group, their positive half-orbits are commensurate after possibly reversing time. Katznelson's criterion then yields conjugacy after a scalar time change. The key new ingredient is a commensuration criterion asserting that a measurable subset of the real line whose symmetric differences with its translates have finite average measure over the unit interval is commensurate with exactly one of the empty set, the whole line, and the two half-lines. This criterion and its application were discovered autonomously by a two-agent AI system. The author independently verified the proofs and prepared the final text.
Given a pair of rational maps $(f, g)$, of degrees $p$ and $q$, each with a parabolic fixed point having a fully invariant simply-connected basin of attraction, we construct an algebraic correspondence $F$ on the Riemann sphere, of bidegree $(pq, pq)$, realizing a mating between the two compositions $g\circ f$ and $f\circ g$ of the maps, and the parabolic faithful discrete representation of the free product of cyclic groups of orders $p + 1$ and $q + 1$. We also show that $F$ is the composition of a pair of deleted covering correspondences of rational maps which are conjugated to polynomials of degrees $p + 1$ and $q + 1$. We generalize our method to construct matings between compositions of pairs of polynomials and (non-parabolic) faithful Kleinian representations of the same group, now with connected regular set. As far as we are aware, these matings between pairs of maps and groups are the first examples that are not time-reversible (that is, they are not conjugate to their own inverses).
The $p$-adic method for Dynamical Mordell--Lang often reduces a residue class of an orbit to the zero set of a locally analytic interpolating function. This paper assumes the standard interpolation, Strassmann, Mahler, and Weierstrass tools, and studies the effective local zero-bound problem that remains after interpolation: certifying the Strassmann index of the resulting one-variable analytic function. We give finite certificates for this index, including finite-data, finite-precision, refined-tail, adaptive residue-class, one-shot, and first-order escape criteria. Since the Strassmann index is a rigorous upper bound for zeros in $\Zp$ and, through Weierstrass preparation, a root count on the closed disc over $\Cp$, these certificates give checkable stopping criteria for local orbit-intersection computations. A residue-class zooming principle replaces a congruence class of times by the iterate $f^{p^h}$, gaining $h$ additional powers of $p$ in the certificate tails. We also introduce an arc-ideal viewpoint for target varieties defined by several equations, replacing a chosen hypersurface bound by the one-variable gcd of all defining equations along the interpolated orbit. In dimension one, the method identifies the certified Strassmann index with the corresponding local Weierstrass root count in the orbit ball. Applications include certified bounds for intersections of non-fixed power-map orbits with finite target sets, and root-of-unity avoidance for maps tangent to the identity at torsion units.
Finite-time blow-up is one of the ways in which a dynamical model can become singular, often signaling the breakdown of either the modeled physical system or the model itself. Determining whether blow-up occurs, and for which parameter values and initial conditions, is therefore a fundamental problem in the analysis of nonlinear dynamical systems. We develop a numerical framework for identifying regions of parameter space in which a dynamical system governed by a system of first-order ordinary differential equations with polynomial right-hand sides exhibits finite-time blow-up for at least one initial condition. The approach combines compactification of the phase space with computational algebraic techniques, producing partitioned parameter landscapes that reveal blow-up and non-blow-up regimes. Through several examples, we show that the method replaces problem-specific hand calculations with an automated computational tool for analyzing blow-up regions in parameter-dependent dynamical systems.
The notion of circuit refers to a cyclic oriented influence between the elements of a dynamical system. There are two classes of circuit: positive and negative. R. Thomas conjectured that a necessary condition of multi stationarity is the existence of positive circuits. In this paper we use dynamical system tools and planar analysis to find conditions for which the conjecture holds for planar systems.
Recent foundation models (FMs) for zero-shot reconstruction of dynamical systems (DS) achieve strong out-of-domain generalization but provide little insight into the mechanisms that underlie their forecasts. Such an understanding could help to strip down overladen FM architectures to their bare essence and expose the minimal requirements for in-context learning in the DS domain. Toward this goal, here we iteratively reduce a recent powerful SOTA model for DS reconstruction, DynaMix (Hemmer & Durstewitz, 2025), to a minimal interpretable two-parameter form, which we call DynaBase. DynaBase produces forecasts through a linear blend of the current latent state and the nearest in-context neighbor and its temporal successor. Surprisingly, despite its extreme simplicity, DynaBase produces highly competitive zero-shot DS reconstructions across chaotic and cyclic systems, with a negligible parameter load, many orders of magnitude below that of other FMs. Even more, this extreme simplicity permits direct model optimization on DS reconstruction measures, as well as closed-form one-step analytical solutions on prediction MSE. Theoretical and empirical analysis of DynaBase further leads to a 1-parameter family of maps, with the context-parroting algorithm of (Zhang & Gilpin, 2026) recovered at one end, and chaotic (divergent but bounded) behavior at the other. We further show how different training strategies lead to models either optimal for short-term prediction or for DS reconstruction. Thus, DynaBase not only exposes the minimal mechanisms required for producing zero-shot DS reconstruction, but also reconciles within an accessible mathematical frame divergent observations in the literature.
We prove an all-directions Marstrand-Mattila projection theorem for self-affine measures and sets in $\mathbb{R}^d$. Under exponential separation, together with proximality and strong irreducibility assumptions on the linear parts, the projection of a self-affine measure onto every line has the expected Hausdorff dimension. If the proximality assumption is strengthened to strong pinching, then the same conclusion holds for the self-affine set $X$ itself, without any separation assumption. In the plane, strong irreducibility of the linear parts alone suffices, and this is sharp. As a corollary, if $X$ additionally has upper Minkowski dimension at most one, then its Minkowski dimension exists and equals its Hausdorff dimension, giving a partial affirmative answer to the folklore question of whether the Minkowski dimension exists for every self-affine set.
We give an explicit algorithm to construct aperiodic tile sets based on Sturmian words of quadratic slopes. The method works for any quadratic irrational slope, and we can produce an aperiodic tile set whose underlying scaling constant is a unit of any real quadratic field. There are two key ingredients in our construction. The first one is the ``Sturmian lattices'', an interesting grid structure generated by Sturmian words that emerged in an aperiodic monotile called Smith Turtle. The second is the bounded displacement equivalence of Delone sets, which plays a central role in this construction. A classification of Sturmian lattices and complete proofs are given in the full version.
Let $S_g$ be a closed orientable surface and let $Ψ\colon \mathrm{Mod}(S_g)\to \mathrm{Sp}(2g,\mathbb{Z})$ be the representation induced by the action on first homology. We investigate the characteristic polynomials of integral symplectic matrices arising from mapping classes of algebraically finite type and give a complete characterization in the cyclotomic case: for $n\geq 3$, the polynomial $\varphi_n(x)$ is realized by a mapping class of algebraically finite type if and only if $n$ has at most two distinct prime divisors. Consequently, if $n$ is square-free and has at least three distinct prime divisors, then every mapping class with characteristic polynomial $\varphi_n(x)$ is pseudo-Anosov. This gives a cyclotomic complement to the Casson--Bleiler homological criterion and yields a complete criterion for a symplectic polynomial to be realized only by pseudo-Anosov mapping classes.
We study a mixed quantum Ising-$XY$ model on the semi-infinite rooted Cayley tree of order two. For every vertex $u$, the edge $\langle u,(u,1)\rangle$ carries an $XY$ interaction and the edge $\langle u,(u,2)\rangle$ carries an Ising interaction. Using the compatibility criterion for tree-indexed quantum Markov chains and consistently working with the normalized trace, we derive the translation-invariant boundary equation and compute explicitly the associated local transfer operator, namely the one-step partial-trace map which propagates successor boundary data to the parent vertex. We prove that the boundary equation has a unique positive translation-invariant solution for all $J_I,J_{XY}\in\mathbb R$ and $β>0$. Hence the model admits a unique translation-invariant quantum Markov chain generated by a positive translation-invariant boundary condition. We also show that the reduced boundary-law dynamics, i.e. the induced finite-dimensional recursion for the boundary-law parameters, has no admissible periodic points of period greater than one and compute the local two-site entanglement on the natural three-site cluster of the tree.
Flow matching has emerged as an effective framework for learning complex data distributions, but adapting pretrained flow models to new tasks often requires computationally expensive retraining. Post-training guidance provides a more efficient alternative, but existing methods are largely heuristic and offer no explicit stability guarantees. We address this limitation by proposing LyaGuide, a unified Lyapunov-guided framework that formulates flow guidance as a Lyapunov control problem. Our main theoretical result establishes an equivalence between guided flow matching and Lyapunov control, thereby unifying common guidance strategies, such as classifier guidance, reward guidance, and energy-based guidance, within a single control-theoretic framework. To enforce the Lyapunov condition, we introduce a pseudo-projection operator with a closed-form expression that endows learned or heuristic guidance terms with explicit stability guarantees. LyaGuide supports two practical settings: a model-driven setting, where the target guidance distribution is specified through a known Lyapunov function, and a data-driven setting, where the guidance is adapted from task-specific downstream data. LyaGuide is compatible with existing guidance methods, introduces minimal additional computational overhead, and is straightforward to integrate in practice. Extensive experiments on synthetic benchmarks, image inverse problems, reinforcement learning planning, and energy-based modeling demonstrate consistent improvements in sample quality, guidance fidelity, and robustness, while maintaining computational efficiency.