Dynamical Systems
Papers
Fresh from arXiv, typeset properly. Highlighted words trade on the market.
Thursday 8 October
We study the decision problem of global asymptotic stability (GAS) of the origin for homogeneous polynomial vector fields with rational coefficients. In the first part of the paper we give, for every fixed dimension $n$ and odd degree $2k+1$, a computable map which sends a homogeneous vector field of degree $2k+1$ on $\mathbb{R}^n$ to a homogeneous cubic vector field on $\mathbb{R}^m$, $m=\binom{n+k-1}{k}$, and preserves both GAS and Lyapunov stability. The map passes to new variables $y_\alpha=x^\alpha$ (one for each monomial $x^\alpha$ of degree $k$), in which the vector field becomes cubic, and adds a cubic penalty term transverse to the set of points of the form $(x^\alpha)_{|\alpha|=k}$; the main point is that this controls every initial state of the new system. As a consequence, the GAS problems for all degrees, for each fixed odd degree $d\ge3$, and for degree three are pairwise many-one equivalent. In the second part we prove that there are a fixed dimension $N_*$ and a fixed odd degree $d_*$ such that the set of homogeneous vector fields on $\mathbb{R}^{N_*}$ of degree $d_*$ whose origin is GAS is many-one complete. The reduction encodes an undecidable family of polynomial equations in integer unknowns into the constant term of a product of Jacobi theta series, realizes these series as outputs of a polynomial differential system with rational coefficients, and turns the existence of an integer solution into the sign of a maximal invariant average, which homogeneous GAS detects. Combining the two parts shows that GAS of rational homogeneous vector fields on $\mathbb{R}^N$ of degree $d$ is many-one complete for every odd $d\ge3$ and every sufficiently large $N$; in particular it is undecidable for cubic vector fields in a fixed dimension. Building on recent work, some analogous results are proven for Lyapunov stability of homogeneous vector fields.
Combinatorics
The Hilton-Milner theorem in the $\ell_2$-norm
The Hilton-Milner theorem is a classical theorem in extremal set theory which determines the maximum size of a nontrivial $k$-uniform intersecting family. In this paper, we prove an exact version of the Hilton-Milner theorem in the $\ell_2$-norm. Our method also gives an asymptotic version of the Hilton-Milner theorem in the $(j, p)$-norm.
Read the full abstract
Analysis of PDEs
Estimates on Fully Nonlinear Elliptic Equations with Unrectifiable Exceptional Gradient Sets
We study fully nonlinear uniformly elliptic equations that are imposed only when the gradients of solutions lie outside a closed ball and a closed purely one-unrectifiable set that can be unbounded. We prove Harnack inequalities and Hölder regularity estimates for viscosity solutions, with constants independent of the exceptional set. The proofs use the geometry of the contact sets in the cusp and barrier arguments to exclude exceptional gradients. We also provide examples and counterexamples showing the role of pure one-unrectifiability. Materials and proofs are presented in detail to ensure the paper is fully self-contained and accessible to a broad audience.
Read the full abstract
Probability
On Large Drift Enhancement for the Relaxation Time of Averages
In this note we prove in a broad analytic set-up a general result characterizing the asymptotic vanishing of the $L^2$-relaxation time of averages in the presence of a large incompressible drift-like term. The broad analytic set-up is analogous to that of Constantin-Kiselev-Ryzhik-Zlatos '08 or of Coti Zelati-Delgadino-Elgindi '20. A crucial role is played by the singular value gap method of Chatterjee '25, further developed by Eberle-Lörler '26.
Read the full abstract
Information Theory
Commitment over Gaussian channels with deterministic identification codes
We establish a previously unnoticed connection between commitment (BC), a fundamental cryptographic primitive, and deterministic identification (DI), a post-Shannon communication setting. Specifically, for additive white Gaussian noise channels $\mathcal{G}$ we show that a reliable commitment protocol can be obtained from deterministic identification codes in both the honest but curious and fully dishonest security settings. This viewpoint allows BC schemes to inherit the asymptotic performance of DI codes. Indeed, we show that commitment is naturally achievable in the linearithmic regime, i.e., with message sets of size $N_n=\exp [\Theta(n\log n)]$, and we establish a lower bound on the linearithmic-scale BC capacity of $\dot C_{\text{hBC}}(\mathcal G)\geq\frac12$ in the honest but curious case and $\dot C_{\text{BC}}(\mathcal G)\geq\frac14$ in the fully dishonest picture. This refines the best previously known characterisation of the commitment capacity over $\mathcal{G}$, which established only an infinite linear-scale capacity.
Read the full abstract
Numerical Analysis
A neural characteristic mapping method: Lagrangian PINNs based on flow maps for transport-dominated problems
We propose a novel numerical method for transport-dominated problems using Physics-Informed Neural Networks (PINNs). In such problems, the solution can develop large gradients and fine structures over time. This is particularly challenging for classical PINNs, which represent the solution itself. In our approach, rather than the solution, we choose to approximate the flow map of the equation. The solution is then recovered by pulling back the exact initial condition, so that its fine structures are produced by the map rather than represented by the network. Our method is constructed by approximating the flow map of the equation with multiple PINNs, one per time subinterval, and composing them using the semigroup property of the flow, so that each network only has to represent a map close to the identity. To highlight the advantages of our method compared to classical PINNs or classical numerical schemes, we present numerical results on the linear advection equation, the incompressible Euler equation in vorticity formulation, the Vlasov-Poisson equation, and finally a drift-kinetic equation.
Read the full abstract
Combinatorics · Cross-listed from cs.CG
Tiling 3D by Translates of a Single Polycube is Undecidable
We prove co-RE-completeness, and thus undecidability, of the following problem: given a single (connected) polycube, decide whether it tiles 3D Euclidean space by translations. We reduce from Wang tiling using the decorated two-prime Sudoku construction of Greenfeld and Tao and a cyclic encoding adapted from OpenAI's 3D aperiodic tile, and apply a reduction of Kim to make the prototile connected (via faces). Dimension three is optimal: translational monotiling is known to be decidable in $\mathbb{Z}^2$ and for a single (possibly disconnected) polyomino in $\mathbb{R}^2$.
Read the full abstract
Analysis of PDEs
Endpoint Compactness and Finite-Energy Weak Solutions for the 2D Isothermal Compressible Navier-Stokes
We prove global existence of finite-energy renormalized weak solutions to the 2D isothermal compressible Navier–Stokes equations in bounded $C^2$ domains with no-slip boundary conditions. The initial data are arbitrary in the natural energy class, allowing vacuum and unbounded density. This resolves the isothermal endpoint $\gamma=1$ of the classical 2D finite-energy renormalized existence theory for power-law pressures. We develop a new endpoint compactness argument based on a truncation-independent integrable bound for centered density oscillations. This argument yields strong $L^1$ compactness of the approximating densities at the natural energy level.
Read the full abstract
Mathematical Physics · Cross-listed from hep-th
Platonic self-dual solitons in a (1+1)-dimensional $\mathbb{CP}^1$ model
We construct a family of self-dual $(1+1)$-dimensional $\mathbb{CP}^1$ sigma models whose potentials $V$ are determined from pre-potentials $U$ by Platonic geometry. Using the round metric on $\mathbb{S}^2\simeq\mathbb{CP}^1$ and the generalized self-duality formalism, we build symmetry-invariant pre-potentials $U$ from the unit vertices of each Platonic solid intersecting the unit circumsphere. The regular-lattice sites associated with vertex, face-center, and edge-midpoint directions are critical points of $U$ and degenerate vacua of the scalar potential $V$. Dual solids share the same regular lattice but generate different pre-potentials, charges, potentials, and metric-gradient flows. Numerical integration reveals distinct self-dual trajectories and energy-density profiles, including multi-peak structures for paths passing near saddle points, while solutions with the same endpoint difference $\Delta U$ have the same total energy. The positive maxima of $V$ are also organized into symmetry orbits. This construction provides a geometric prescription for controlling the vacuum structure and self-dual sectors of $\mathbb{CP}^1$ models with a potential.
Read the full abstract
Algebraic Geometry
Local Constancy of the Betti Automorphic Category over the Moduli of Curves
Nadler and Yun conjectured that the automorphic side of the Betti geomet- ric Langlands correspondence forms a local system of categories over the moduli of curves. We prove this local constancy conjecture using purely microlocal methods, independent of the recently estab- lished Betti geometric Langlands conjecture. In particular, it follows that the Betti automorphic category depends only on the genus of the curve, verifying a conjecture of Ben-Zvi and Nadler.
Read the full abstract
Mathematical Physics · Cross-listed from hep-th
Differential $T_2$-Duality and Spans of Principal 3-Bundles with Connections
We present a generalization of T-duality, called $T_2$-duality, that captures certain aspects of U-duality in M-theory. In particular, it features backgrounds that contain a 2-gerbe over a principal torus bundle, structures familiar from eleven-dimensional supergravity backgrounds. In special cases, it also incorporates S-duality in a precise sense. Our main result is a theorem characterizing $T_2$-dual pairs through spans of categorified principal bundles with connections, giving explicit rules for constructing a $T_2$-dual for a given background. A number of explicit examples illustrate how the various characteristic classes are mapped between the left and right backgrounds. For two-dimensional torus fibers, we also relate our construction to the combination of dimensional reduction and T-duality that links eleven-dimensional supergravity to type IIB supergravity. Finally, we discuss how $T_2$-duality arises from a Buscher-like span of 3d sigma models.
Read the full abstract
Information Theory · Cross-listed from cs.GR
FeatureZ: A General Framework for Feature-Preserving Compression via Pointwise Bounds and Star Classification
Geometric and topological features, such as isosurfaces, quantiles, merge trees, and Morse-Smale complexes, are central to the analysis and visualization of scientific data across diverse domains, including medical imaging, climate science, materials science, and astronomy. However, most lossy compressors for scientific data provide only pointwise error guarantees and do not preserve derived features. Existing feature-preserving compressors are often difficult to develop and typically tailored to a single feature. In this paper, we introduce FeatureZ, a lossy compression framework for structured volumetric scalar fields that can preserve a wide class of geometric and topological features. In particular, FeatureZ frames feature-preserving compression as preserving pointwise upper and lower bounds together with a consistent classification of each point's star (i.e., its incident cells) in the underlying structured mesh. Although not all notions of feature preservation fit this framework, many features and topological descriptors targeted by existing compressors do. FeatureZ provides an efficient implementation of feature-preserving compression under this formulation. In particular, FeatureZ operates as an augmentation layer that refines the output of an existing lossy compressor. It first applies quantization to enforce pointwise bounds, and then employs an iterative procedure to ensure consistent star classification. We demonstrate that FeatureZ preserves diverse features during compression with minimal overhead, achieving compression ratios comparable to or better than methods specialized for individual features.
Read the full abstract
Dynamical Systems · Cross-listed from cs.LG
Bilevel optimization for data-driven learning of Koopman embeddings using kernel-based autoencoders
Koopman operator theory provides a linear framework for analyzing nonlinear dynamical systems and has become a major tool for data-driven modeling. A central challenge, however, is that finite-dimensional approximations computed by methods such as extended dynamic mode decomposition (EDMD) require the dictionary to be specified a priori. Recent machine-learning approaches address this limitation by learning the dictionary from data, predominantly using artificial neural network (ANN) autoencoder architectures. Although kernel methods offer an alternative with greater interpretability and tractability for theoretical analysis, they have received little attention in this setting. We introduce extended dynamic mode decomposition with kernel-based dictionary learning (EDMD-kDL), a kernel-based method for learning finite-dimensional Koopman embeddings directly from data. The method combines ideas from collocation methods and bilevel optimization to simultaneously learn a kernel dictionary and the corresponding Koopman approximation. We evaluate EDMD-kDL against state-of-the-art ANN-based approaches on a range of numerical experiments, including global sea-surface-temperature forecasting and learning directly from video data. Across all tested settings, EDMD-kDL achieves performance comparable to or better than the ANN-based methods. Moreover, in contrast to standard kernel methods, the proposed approach is scalable to large datasets by design since the size of the required kernel matrices depends on the number of collocation points rather than the size of the training dataset.
Read the full abstract
Probability
Directed fractal percolation and non-Lipschitz variants
While models of undirected and directed geometry in random i.i.d. disorder with rapidly decaying tails are expected to behave similarly and exhibit features of the Kardar–Parisi–Zhang universality class, the scenario when the background noise is fractal is expected to be significantly different. A canonical example of the latter is Liouville Quantum Gravity, where the planar Gaussian free field forms the fractal background. A toy model capturing some of the essential features is given by Mandelbrot's fractal percolation, a random Cantor set. In the simplest setting, dyadic cubes of different scales are retained independently with probability $p$ (a parameter of the model), and the intersection of all retained cubes forms the set of open sites. Connectivity properties of such sets have been intensely studied. Across [CCD88, Cha95], a surprising result was proven. Namely, while there is a phase transition for connectivity of the fractal percolation set (like in usual bond percolation), directed percolation never occurs. Results of a similar spirit have been obtained in geometric measure theory as well. As a step towards understanding last passage percolation (LPP) driven by the Gaussian free field (initiated in [GGN24]), we study the problem of quantifying LPP in fractal percolation. Developing a multi-scale framework, we show an almost polynomial gap from linear growth of the passage time. The obstruction to high LPP values is the Lipschitz nature of directed paths. Along these lines, [Cha96] showed that any connected subset of fractal percolation must have Hausdorff dimension strictly larger than one. We show that one can construct directed paths passing through open sites in fractal percolation provided they move exponentially fast in the spatial direction—a model we introduce and term as unrectifiable directed percolation.
Read the full abstract
Algebraic Geometry
The moduli space of curves of genus 16 is uniruled
Using in a novel way non-abelian Brill-Noether theory, we show that the moduli space of curves of genus 16 is uniruled. This is the highest genus for which such a result is known.
Probability
Subspace Uncertainty and Sharp Sampling Thresholds on the Boolean Cube
We study Gaussian regression under squared population $L_2$ loss in a known $m$-dimensional subspace of degree-at-most-$k$ functions on the $d$-dimensional Boolean cube. Random inputs can undersample regions essential for prediction, delaying the parametric rate even when the model is known. For fixed $q_0<1/2$, $1\le k\le q_0d$, and sufficiently large fixed $A$, the worst-subspace sample threshold for minimax error $A\sigma^2(m+t)/n$ with confidence $1-e^{-t}$, $t\ge\log4$, is \[ N=(m+t)\exp\{E_{d,k}+O(k^{1/3})\}, \quad E_{d,k}=d\Psi(k/d), \] where $\Psi(q)=\log2-\mathsf H(\tfrac12-\sqrt{q(1-q)})$ and $\mathsf H$ is binary entropy with natural logarithms. The upper bound holds for every feasible $m$; the matching lower bound holds when $m\le\binom d{\lfloor k^{1/3}\rfloor}$ or $t\ge m$. We sharpen the Polyanskiy–Samorodnitsky uncertainty principle in two respects. First, for fixed leakage $\rho\in(0,1)$, the smallest set carrying a fraction $1-\rho$ of a nonzero degree-at-most-$k$ polynomial's energy has probability $\exp\{-E_{d,k}+O_{\rho,q_0}(k^{1/3})\}$. An Airy-kernel construction proves that the remainder cannot be $o(k^{1/3})$ in general. Second, we construct a subspace of dimension $\binom d{\lfloor k^{1/3}\rfloor}$ such that every function in the subspace has at least a fraction $1-\rho$ of its energy on the same set, whose probability is at most $\exp\{-E_{d,k}+C_{\rho,q_0}k^{1/3}\}$. For sufficiently large $k$, this set is a Hamming ball. A striking consequence is an exponential cost of noise: the parametric rate can require $(m+t)4^k\exp\{-O(k^{1/3})\}$ samples, whereas $O((m+t)2^k)$ suffice for noiseless identification. As $k\to\infty$ with $k/d\to0$, the noisy threshold is $(m+t)\exp\{2k+o(k)\}$.
Read the full abstract
Information Theory · Cross-listed from quant-ph
Improved Local Leakage Resilience of Shamir Secret Sharing and Worst-Case Optimal Polynomial Intersection
We study two problems: Local Leakage Resilience (LLR) for Shamir secret sharing, and worst-case Optimal Polynomial Intersection (OPI). Both problems concern polynomials $Q(X)$ of degree less than $k$, over a prime-order finite field $\mathbb{F}_p$. In LLR for Shamir secret sharing, one asks how much one can learn about $Q(0)$ given a few bits leaked from each of $Q(\alpha_1), \ldots, Q(\alpha_n)$, for distinct non-zero evaluation points $\alpha_i \in \mathbb{F}_p$. In OPI, one is given input list $S_1, \ldots, S_n \subset \mathbb{F}_p$, and wants to find a polynomial $Q(X)$ of degree less than $k$ so that $Q(\alpha_i) \in S_i$ for as many $i$ as possible. Leveraging recent connection between these two problems due to (Sun, Wootters 2026), we improve the state-of-the-art for both problems. For LLR, we show that there is some constant $\delta > 0$ so that, as long as $R := k/n \geq 1/2 - \delta$, Shamir secret-sharing is one-bit locally leakage resilient (meaning that one can learn only a negligible amount about $Q(0)$). This is the first result to break the so-called "one-half barrier" for LLR, and improves over the previous best known result, requiring $R \geq 0.668$ (Kasser, 2025). For OPI, we give a quantum algorithm that finds a polynomial $Q(X)$ that agrees with at least a $\mathsf{SCL}_\rho(R)-\varepsilon$ fraction of the lists in expectation, for every fixed $\varepsilon>0$, where $\mathsf{SCL}_\rho$ is the semicircle law of (Jordan et al., 2025). This improves previous algorithmic (and existential) results of (Jo, 2026) and (Horinaga, Yamakawa, 2026). We also give further improved existential results. We also adapt the hardness result of (Yamakawa, Zhandry, 2024) to apply to OPI (rather than a folded version); over large fields, this gives an unconditional separation between the quantum and classical hardness of OPI relative to a membership oracle.
Read the full abstract
Probability
Brownian motion with drift conditioned to have restricted $L^2$-norm
We study the long-time behaviour of Brownian motion with drift under a rare-event conditioning of the quadratic additive functional $Z_T=\int_0^T W_s^2ds$, where $(W_t)_{t\geq 0}$ is a Brownian motion with drift $\mu$. More precisely, we consider the conditional law of $W$ given the event $Z_T\le \theta T$, as $T\to\infty$, allowing the drift $\mu$ to depend on $T$. We derive sharp, uniform small-deviation asymptotics for $Z_T$, including the exact prefactor, and use them to show that the conditioned process converges weakly to an Ornstein–Uhlenbeck process. Remarkably, the limiting dynamics are independent of the original drift.
Read the full abstract
Optimization and Control
Asymptotically Vanishing Excitation Without Loss of Parameter Convergence
Parameter convergence in recursive identification requires persistent excitation, which is typically enforced with a probing signal that degrades performance if retained indefinitely. Regulating the probing amplitude with the one-step prediction error is insufficient, since a small prediction error does not imply parameter convergence. This paper proposes a self-regulating excitation law that scales the probing amplitude with the associated covariance matrix, so that the excitation vanishes only as fast as the estimator's own measure of remaining parameter uncertainty allows. The key idea is that the excitation amplitude must decay more slowly than the covariance measure, the largest eigenvalue of the covariance matrix in this work, for parameter convergence to be preserved. Under the proposed excitation law, the excitation amplitude and the covariance measure are shown to jointly converge to zero while the parameter estimate converges to the true parameter, with the excitation vanishing more slowly than the covariance measure throughout. A scalar exponent in the law trades the rate of covariance decay against the amount of excitation retained. Numerical results confirm both properties.
Read the full abstract
Geometric Topology
The Colored Jones Function does not determine any colored Links–Gould polynomial
We prove that there is no uniform description of $\mathrm{LG}^{(n)}$ in terms of colored Jones polynomials for all knots, even when allowing infinitely many colors and arbitrary operations. We prove the $n=1$ case, using an existing cabling formula and the separation of the Conway and Kinoshita–Terasaka knots by $\mathrm{LG}^{(2)}$. We develop an odd cabling formula and use a Vandermonde argument to show that for each $n\geq 2$ some satellite of each of these knots is distinguished by $\mathrm{LG}^{(n)}$.
Read the full abstract
Differential Geometry
Prevalence of sparsity for negatively curved metrics
Let $M$ be a closed manifold of arbitrary dimension. The length spectrum of a negatively curved metric $g$ on $M$ is the set of lengths of all possible closed geodesics of $g$. We prove that metrics whose length spectrum is exponentially sparse, i.e., the gaps between distinct lengths are bounded below by a quantity decaying exponentially in the length, are prevalent, and hence dense, among $C^k$-smooth negatively curved metrics on $M$ for sufficiently large $k$. Moreover, unlike in all known constructions, the exponent in our bound grows only sublinearly in the regularity $k$. This work is motivated by our recent result on prevalent unmarked length spectral rigidity for expanding circle maps and a very recent result by DeWitt, Durham, Reber, and O'Hare that $C^k$-smooth negatively curved metrics with exponentially sparse spectra are locally rigid.
Read the full abstract
Geometric Topology
Coloring graphs of homologous spheres
The curve graph of a surface is the graph whose vertices are isotopy classes of essential, closed curves on the surface. Edges join vertices with disjoint representatives. The curve graph plays a central role in the theory of mapping class groups of surfaces. Gaster, Greene, and Vlamis have shown that the subgraph of the curve graph induced by curves in a single primitive homology class is uniquely and finitely colorable. In the analogy between the theory of mapping class groups of surfaces and outer automorphisms of free groups, the sphere graph of $M_r$, a connect sum of $r$ copies of $S_1 \times S_2$ plays an analogous role in the theory. Motivated by Gaster, Green, and Vlamis' result, we investigate $\mathcal{S}_v(M_r)$, the subgraph of the sphere graph of $M_r$ induced by spheres representing a single primitive homology class $[v] \in H_2(M_r;\mathbb{Z})$. We show that this graph is uniquely $r-1$ colorable. As a corollary, with $r\ge 3$ we obtain a nontrivial homomorphism from the Torelli subgroup $IO_r \le \operatorname{Out}(F_r)$ to the symmetric group on $r-1$ symbols; this parallels Gaster, Greene, and Vlamis' connection between their coloring and the Chillingworth homomorphism. In the course of our proof we find new smaller generating sets for mapping class groups of $M_r$ with spheres deleted.
Read the full abstract
Numerical Analysis
Volume Sampling and Spectral Equalization in Randomized Alternating Projections
We study a family of randomized alternating projection methods that alternate among subspaces spanned by subsets of \(n\) vectors from a prescribed set. We derive an explicit characterization of performance bounds when these subspaces are sampled with probabilities proportional to the volumes subtended by the vectors in their corresponding subsets. For each \(n\), the bound is obtained through an explicit nonlinear transformation of the spectrum of an associated matrix. Our analysis reveals an explicit spectral equalization mechanism that drives the spectrum toward improved conditioning as $n$ increases, establishing an unexpected connection between volume sampling and the Cayley-Hamilton theorem. Furthermore, we introduce a computationally efficient uniform sampling scheme that achieves comparable theoretical guarantees through a specific relaxation strategy. Empirical results show that this relaxation method achieves convergence rates comparable to, and often better than, volume sampling, while avoiding the infeasibility of volume sampling at large scale. Besides the randomized Kaczmarz, these results also directly apply to randomized Gauss-Seidel and coordinate descent methods.
Read the full abstract
Number Theory
Geometrically Abelian Sections on Higher-Dimensional Varieties
We investigate how much arithmetic information about a variety can be detected by sections of the geometrically abelian fundamental groups of sufficiently small open subvarieties. Over a \(p\)-adic field, such a section determines a rational point on the Albanese torsor, and, under a conjectural assumption, this point is represented by a \(0\)-cycle of degree \(1\) on the original variety if the variety has index \(1\). For curves, the latter statement holds unconditionally. Over finitely generated fields, such sections still force the relative Brauer group to vanish away from the characteristic. The key input is an identification of the Chern class obstruction attached to a section with the obstruction arising from local Tate duality. As a byproduct, we give a new proof of an open-subvariety strengthening of a theorem of Esnault and Wittenberg.
Read the full abstract
Optimization and Control
Constraint Enforcement via Reference Governor in Dynamic Mode Adaptive Control
This paper integrates reference governance with Dynamic Mode Adaptive Control (DMAC) to enforce constraints without a prior plant model. DMAC identifies and updates a discrete-time model online from measured state and input data, and the reference governor (RG) uses this identified model, rather than the true plant dynamics, to predict constrained closed-loop behavior. Because the prediction model changes as DMAC adapts, a governed reference that was previously feasible may become infeasible, and the classical interpolation-based RG may be unable to move the reference in the direction required to recover predicted constraint satisfaction. The proposed DMAC-RG framework instead searches directly over the governed reference rather than over an interpolation factor. When feasible references exist, the governed reference is selected to track the desired command, with an optional penalty on deviation from the previously applied reference. When the finite-horizon feasible set is empty, the reference is selected to minimize predicted constraint violation while retaining the same tracking preference. This work is a numerical investigation and does not provide theoretical guarantees of recursive feasibility, constraint satisfaction, or closed- loop stability. Numerical examples on a linear system, a nonlinear Van der Pol oscillator, and a high-dimensional Burgers equation under partial-state measurement demonstrate that the proposed RG enforces constraints using only the model identified online by DMAC, without access to the true plant dynamics or the complete system state.
Read the full abstract
Numerical Analysis · Cross-listed from cs.MS
asdex: Automatic Sparse Differentiation in JAX
Many tasks in scientific computing and machine learning require the Jacobian or Hessian matrix of a function. Automatic differentiation (AD) computes these derivatives to machine precision, but materializing a dense $m \times n$ Jacobian requires $n$ forward-mode or $m$ reverse-mode AD passes, one per column or row. For a large class of functions, each output depends on only a few inputs, making the derivative matrix sparse. Automatic sparse differentiation (ASD) exploits this structure in four steps: detection of the input-agnostic sparsity pattern, coloring of a graph to group columns or rows that can share an AD pass, compressed differentiation to compute a compressed derivative matrix with one AD pass per color, and finally decompression into the original sparsity pattern. The number of colors, and hence of AD passes, is often independent of the problem dimension: a banded Jacobian with $b$ contiguous bands, for instance, only ever requires $b$ colors, regardless of its size. asdex offers the first standalone ASD toolkit in the popular JAX ecosystem. With $\texttt{asdex.jacobian}$ and $\texttt{asdex.hessian}$, it provides sparse drop-in replacements for $\texttt{jax.jacobian}$ and $\texttt{jax.hessian}$.
Read the full abstract
Probability
Anomalous scaling limit of a Brownian particle in a log-correlated potential
We study a Brownian particle in $\mathbb{R}^d$, $d\geq2$, with drift given by the gradient of a log-correlated Gaussian potential. At weak disorder, we prove convergence to a scaling limit whose law is singular with respect to Brownian motion and whose invariant measure is given by Gaussian multiplicative chaos. The proof uses a renormalization group induction: at each scale, the elliptic operator is well approximated by a Laplacian with effective diffusivity decaying as a power of scale. We compute this exponent exactly in dimension two.
Read the full abstract
Mathematical Physics
Wave packets asymptotics in $\kappa$-minkowski space
In this paper, we study the nature of wave packets and their decay behavior in the $\kappa$-Minkowski space. We employ techniques from singularity theory throughout our analysis. Specifically, we use methods from the theory of Lagrangian singularities and caustics to study the topology of the energy-level manifold. The Newton polyhedron method from the theory of oscillatory integrals is applied to analyze the decay behavior of wave packets. We find that for a generic Hamiltonian in the $\kappa$-Minkowski space, the associated wave packets exhibit different decay rates depending on the control parameter of the system. When the control parameter takes values on the cusp, the system exhibits decay of order $t^{-1/4}$; otherwise, it exhibits decay of order $t^{-1/2}$. This distinctive behavior of wave packets in the $\kappa$-Minkowski makes it fundamentally different from commutative spacetime.
Read the full abstract
Algebraic Geometry
Characteristic cycles and the Gevrey filtration of irregularity in hypergeometric families
We show that the characteristic cycle multiplicities of $A$-hypergeometric systems are upper semicontinuous in the parameter for every rational projective weight, thereby proving a conjecture of Schulze and Walther [Duke Math. J. 142 (2008), 465–509]. The proof uses a specialization formula for filtered families satisfying a fixed product support condition, expressing the difference between the special and generic characteristic cycles as the characteristic cycle of the first Tor module. Combined with the Gevrey index theorem, this yields upper semicontinuity of generic Gevrey irregularity and of each of its graded dimensions. We also construct non-Cohen–Macaulay semigroup rings whose Gevrey dimensions are constant in the parameter despite the presence of a genuine irregular slope, showing that the Gevrey jump locus need not coincide with the full exceptional arrangement defined by Ext.
Read the full abstract
Group Theory
Ordinary $G$-permutable subgroups of finite simple groups
We completely solve part (a) of Problem 17.112 in the Kourovka Notebook. We prove that a non-abelian finite simple group $G$ has a proper non-trivial ordinary $G$-permutable subgroup if and only if $G$ is isomorphic to ${\rm PSL_2(q)}$ where $q\geq4$ is a prime power and $q\not\equiv3\pmod4$, ${}^2B_2(2^{2a+1})$ where $a\geq1$, $J_1$, ${\rm Sp}_4(4)$, or ${\rm PSp}_4(p)$ where $p$ is prime and $p\equiv1,13,17\pmod{24}$. Moreover, we classify all proper non-trivial ordinary $G$-permutable subgroups of these groups. Every other non-abelian finite simple group $G$ has no proper non-trivial ordinary $G$-permutable subgroup.