Numerical Analysis
Papers
Fresh from arXiv, typeset properly. Highlighted words trade on the market.
Tuesday 21 July
Refinement operators generate many functions used in wavelet constructions, subdivision schemes, and geometric modeling. Their finite iterates can develop rapidly increasing numbers of linear pieces, making them a natural test case for the expressive power of deep neural networks. Earlier work showed that, for scalar binary refinement with a finitely supported mask, every compactly supported continuous piecewise linear seed has finite refinement iterates that admit exact ReLU realizations of fixed width and depth growing linearly with the number of refinement steps. The present paper gives a new construction of this known theorem. The difficulty is that the refinement cascade is driven by discontinuous binary digit choices, whereas ReLU networks produce continuous piecewise linear maps. We represent the residual dynamics on a polygonal model of the circle and describe each residual position in two overlapping coordinate systems, one ordinary and one shifted by one half. Their discontinuities occur at different points. The network switches between the two descriptions only where both are valid and the corresponding fixed linear cascade updates agree, so the switch is exact and requires no multiplication by a variable selector. The construction also gives exact readout of every continuous piecewise linear circle function satisfying the natural endpoint compatibility condition. Localized seeds are handled by a two-pass network, and translation covariance, finite decomposition, and gluing extend the result to arbitrary compactly supported continuous piecewise linear seeds in a preserved support window.
Algebraic Topology
Two-Sided Dimension Bounds for the Peterson Hit Problem via Projections and Matrix Minors
Let $\mathcal P_k=\mathbb F_2[x_1,\ldots,x_k]$ be the polynomial algebra over the prime field $\mathbb F_2$, viewed as an unstable module over the mod-$2$ Steenrod algebra $\mathcal A$. The well-known Peterson hit problem asks for a minimal set of generators for the $\mathcal A$-module $\mathcal P_k$. This is equivalent to determining the dimension of the cohit space $(Q\mathcal P_k)_d=(\mathcal P_k/\mathcal A^{+}\mathcal P_k)_d$, where $\mathcal A^{+}$ denotes the augmentation ideal of $\mathcal A$, for every $k\geq1$ and positive degree $d$. Although solved in every degree for at most four variables, it remains a difficult open problem in general. Furthermore, given the limitations of current tools, explicitly determining the dimension of $(Q\mathcal P_k)_d$ in the general case appears out of reach. Motivated by these limitations, we establish explicit upper and lower bounds for this dimension for arbitrary positive integers $k$ and $d.$ Our method combines binary combinatorics, linear algebra, and graph and simplicial structures associated with the generating Steenrod squares. We characterize zero rows, count zero columns, and refine rank estimates using Adem relations. Minors and zero rows of the resulting smaller matrix yield further two-sided cohit bounds without determining a complete basis or computing the full hit rank.
Read the full abstract
Mathematical Physics
Entropy Geometry and Normalized Means on Infinite-Dimensional Hamiltonian Manifolds
We propose a geometric–analytic framework for equilibrium statistical mechanics on infinite-dimensional Hamiltonian systems. In situations where no suitable $\sigma$-additive invariant measure is available, we use normalized means, which generalize probability measures and normalized integrals. This construction yields entropy and free-energy functionals on weak symplectic Fréchet manifolds and gives existence and uniqueness of exponential-family equilibrium states under explicit admissibility and separation assumptions. These states are stationary under Hamiltonian flows preserving both the reference mean and the equilibrium weight, and satisfy a classical Poisson–KMS identity when the reference mean is Poisson invariant. Under a local exponential regularity assumption, the logarithmic partition functional is smooth and convex, with Hessian given by the covariance form. It is strictly convex modulo thermodynamically null directions and, through Legendre–Fenchel duality, induces a concave entropy on the domain of extensive variables. We illustrate the framework with $H^s$-geodesic equations on current groups $\operatorname{Map}(M,G)$ and diffeomorphism groups $\operatorname{Diff}(M)$, including hydrodynamic and field-theoretic examples.
Read the full abstract
Optimization and Control · Cross-listed from cs.DS
Stronger Memory-Query Tradeoffs for Convex Optimization: The Limitations of Subquadratic Memory
We prove two lower bounds for the first order oracle complexity of minimizing a $d$-dimensional $1$-Lipschitz convex function over the unit ball with $m$ bits of memory. We first show that any such (possibly randomized) algorithm must make $\tilde{\Omega}(\frac{d^2}{\sqrt{m}})$ oracle queries. For deterministic optimization algorithms, we show that $\tilde{\Omega}(\min\{d^{1.6},\frac{d^{8/3}}{m^{2/3}}\})$ queries are required. For all memory regimes of interest, these improves upon the previous best known lower bounds of $\tilde{\Omega}(\max\{\frac{d^{8/3}}{m^{4/3}},\frac{d^{4/3}}{m^{1/6}}\})$ and $\tilde{\Omega}(\frac{d^{5/3}}{m^{1/3}})$ for randomized and deterministic algorithms respectively. Notably, due to existing upper bounds, our lower bound for deterministic algorithms is the first to show a sharp oracle complexity phase transition around $m\approx d^2$, where a polylogarithmic change in memory leads to a $\mathsf{poly}(d)$ change in the number of required oracle calls. Further, when the suboptimality is polynomially small in $d$, our lower bound randomized algorithms is the first to show that $\tilde{\Omega}(d^2)$ memory is necessary to nearly match the optimal query complexity among algorithms without memory constraints. Previously, such a result was only known for the regime where the suboptimality is quasipolynomially small in $d$.
Read the full abstract
Combinatorics
Compactness of abundance in asymmetric hypergraph removal lemmas
Fix an integer $r\ge 2$ and a finite simple $r$-uniform hypergraph $F$ with at least one edge and no isolated vertices. An $n$-vertex $r$-graph is $\epsilon$-far from being $F$-free if at least $\epsilon n^r$ edges must be deleted to destroy every copy of $F$. A finite $r$-graph $H$ is $F$-abundant if there are constants $c,C>0$ such that every sufficiently large $\epsilon$-far host contains at least $c\epsilon^C n^{v(H)}$ labelled copies of $H$. A family is $F$-abundant when one member has this lower bound in each host, although the member may depend on the host and on $\epsilon$, while $c$ and $C$ are common to the family. We prove that every abundant family contains an abundant member. We prove the analogous coloured theorem for $F$-partite hosts containing edge-disjoint part-respecting copies of $F$ such that every vertex lies in at least $\epsilon n^{r-1}$ of them. The case $r=2$ yields the coloured and uncoloured graph compactness theorems, answers Question 5.2 of Girão, Hurley, Illingworth and Michel, and proves their Conjecture 5.1 [J. Lond. Math. Soc., 2024]. We also obtain an explicit bound $\lfloor 2r(C+1)\rfloor$ for the order of the non-isolated core of a selected witness. Moreover, we give several applications. For example, we construct translation-invariant linear systems from abundant coloured hypergraphs, obtain a square-root bound for an equation associated with a cycle of bounded length, give a one-sided tester based on one fixed graph when distance from the property gives a polynomial lower bound on distance from being $F$-free, and prove that no algorithm decides whether a family of finite simple graphs enumerated by a Turing machine is $K_3$-abundant.
Read the full abstract
Mathematical Physics
Probability, Curvature and Spectrum on Graphs
We show that, for a metric graph equipped with the Laplacian operator $\Delta=-\frac{d^2}{dx^2}$, the graph trace formula admits a new interpretation in terms of quantum probability and curvature. Our approach is based on a notion of graph curvature inspired by Wilson-lines and holonomy, together with von-Neumann's ergodic theorem.
Read the full abstract
Optimization and Control · Cross-listed from cs.LG
Quasi-Monte Carlo Initialization for Meta-Reinforcement Learning
This paper explores the efficacy of quasi-Monte Carlo (QMC) weight initialization for meta-reinforcement learning within modern benchmark environments. Various sampling methods are used to bound a population-based search and aggregate an optimal prior from a baseline set of tasks. The QMC meta-priors show improvements in training convergence compared to modern orthogonal (SB3) defaults when extrapolated to similar unseen continuous control environments. In dissimilar tasks, the orthogonal orientation was globally superior for an unbiased search.
Read the full abstract
Number Theory
A problem on sumset sizes of sets of lattice points
A central problem in additive number theory is to understand the set of sizes of $h$-fold sums of finite subsets of an additive abelian semigroup. It is proved that the "range of sumset sizes” is the same for finite sets of integers and for finite sets of $n$-dimensional lattice points. The associated geometrical problem is to determine if these sumset size sets can be computed more efficiently with lattice points than with integers.
Read the full abstract
Group Theory
Some answers regarding factorizations of finite groups
In this note we answer some questions of G. M. Bergman concerning factorizations of finite groups. Such questions originated when M. H. Hooshmand asked whether, given a finite group $G$ and a factorization $|G| = n_1 \cdots n_k$, one can always find subsets $A_1, \dots , A_k$ of $G$ with $|A_i| = n_i$ such that $G = A_1 \cdots A_k$. This was Question 19.35 of the Kourovka Notebook. G. M. Bergman provided a counterexample for $k=3$, and Kabenyuk provided counterexamples for all $k \geq 3$. We present these counterexamples with direct proofs that mirror Bergman's original $k=3$ argument, and we answer another recent question of Bergman.
Read the full abstract
Algebraic Topology
The rank-five Peterson hit problem, the fifth Singer transfer, and a geometric generator in unoriented cobordism
The Peterson hit problem seeks a minimal set of generators for the polynomial algebra $P_s = \mathbb{F}_2[x_1,\dots,x_s]$ as an unstable module over the mod-2 Steenrod algebra $\mathcal{A}$. For rank five, general admissible bases fail, and the interplay between Kameko periodicity and modular invariants becomes computationally complex. In this paper, we study the rank-five cohit module in the generic family $N_d = 27\cdot 2^d - 5$. Exact sparse elimination in degree $49$ processes $292825$ monomials, yielding a hit rank of $289969$ and a cohit dimension of $2856$. We determine the exact weight summands and prove that the weight-$(3,3,2,2,1)$ summand is exactly the kernel of Kameko's operation, with dimension $1891$. These exact values systematically correct the corresponding rank-five kernel and dimension assertions in Nguyen Khac Tin's previous paper. An exact invariant calculation shows that the general linear group invariants in degree $49$ form a one-dimensional space generated by a $283$-term polynomial, and we prove that the fifth Singer cohomological transfer is an isomorphism in this family. Geometrically, the Hilbert-Poincare series of the unoriented cobordism ring gives the dimension of the degree-$49$ cobordism group as $5692$. We prove that the Milnor hypersurface $H_{2,48} \subset \mathbb{R}P^2 \times \mathbb{R}P^{48}$ represents the unique nonzero indecomposable class by computing a tangential Stiefel-Whitney number, providing an explicit geometric generator. However, the evident map from $H_{2,48}$ to the classifying space $B(\mathbb{Z}/2)^5$ sends its fundamental class to a homology class with nonzero $Sq^2_*$. Consequently, this geometric generator cannot be identified with the functional dual of the algebraic invariant, establishing a precise boundary between the Steenrod-theoretic invariant line and the geometric cobordism generator.
Read the full abstract
Combinatorics
Counting spanning quasi-trees of ribbon graphs: determinants and #P-completeness
A quasi-tree of a connected ribbon graph is a spanning ribbon subgraph with exactly one boundary component; quasi-trees play the role of spanning trees in the topological graph theory of embedded graphs. We prove that counting them is #P-complete under polynomial-time Turing reductions, already for bouquets. The proof identifies every nonempty framed chord diagram, up to natural identifications, with a 4-regular map equipped with a distinguished A-trail, in such a way that quasi-trees correspond to A-trails, whose counting is #P-complete by a theorem of Ge and Štefankovič. Through the framed Cohn-Lempel equality the count is also an interlace-polynomial evaluation - $q(H;2,1)$, the number of full-rank induced subgraphs of the looped circle graph $H$ of the diagram - placing it on the line $y=1$ left open in the complexity classification of Bläser and Hoffmann; a cloning argument then makes every fixed rational point of that line, other than the trivial $(1,1)$, #P-hard on looped circle graphs, even when a framed chord representation is supplied. On the tractable side, the same GF(2) model yields short proofs of the known determinantal cases: for orientable ribbon graphs the count is a determinant, essentially the Matrix-Quasi-tree Theorem of Merino, Moffatt and Noble, proved here via Bouchet's principal unimodularity, and for bouquets with exactly one non-orientable loop it is a sum of two orientable determinants, equivalent by a rank-one determinant identity to the determinant formula of Deng, Jin and Yan.
Read the full abstract
Combinatorics
Negative-index $u$-deformed homogeneous functions and Ramanujan partial theta functions
We introduce and study a class of \(u\)-deformed homogeneous functions of complex order \(\rr_{\alpha}(x,y;u\,|\,q)\), obtained by extending \(u\)-deformed homogeneous polynomials beyond the nonnegative integer regime. We establish their basic analytic and algebraic properties, including convergence criteria, recurrence relations, \(q\)-derivative formulas, and an operator-theoretic realization through a deformed \(q\)-exponential operator. A central feature of the theory appears in the negative-index sector: for every nonnegative integer \(n\), the functions \(\rr_{\minus n\minus1}(x,y;u\,|\,q)\) admit finite triangular decompositions in terms of shifted Ramanujan partial theta functions. We also prove the inverse triangular transformation, showing that the corresponding partial theta functions can be recovered from the negative-index \(u\)-deformed homogeneous functions. We also derive pantograph-type \(q\)-difference equations for these functions. Finally, we derive specialized forms of the negative-index theory in Cauchy and Stieltjes–Wigert-type regimes, including a boundary geometric case and genuine partial-theta regimes.
Read the full abstract
Analysis of PDEs
Inverse Modified Scattering for Cubic NLS with a Repulsive Delta Potential
We study the one-dimensional cubic nonlinear Schrödinger equation with a repulsive delta potential and a localized inhomogeneous coefficient. We prove small-data modified scattering and construct the associated vector-valued modified scattering map, whose two components encode the coupling of the frequencies $\xi$ and $-\xi$ induced by the point interaction. We show that this map determines both the strength of the delta potential and the inhomogeneous coefficient. Quantitatively, the delta strength is recovered with Lipschitz stability, while the inhomogeneous coefficients satisfy a Hölder stability estimate. To our knowledge, these are the first recovery and stability results for a modified scattering map in the presence of an external potential.
Read the full abstract
Information Theory
CRB-Driven Beamforming and Trajectory Optimization for UAV-assisted ISAC System
In this paper, we study an unmanned aerial vehicle (UAV)-assisted integrated sensing and communication (ISAC) system, where a UAV enhances the sensing capability of a base station (BS) towards a target while ensuring reliable communication towards a downlink user. This architecture is practically attractive for future wireless networks due to the UAV's controllable mobility and adaptive sensing coverage in wireless environments. The sensing performance is characterized by the average Cramér-Rao bound (CRB), which quantifies the minimum variance of the unbiased angle-of-arrival estimation. To enhance the sensing performance, the UAV trajectory and beamforming parameters are jointly optimized under power and mobility constraints, while satisfying communication requirements to the downlink user. To address the resulting non-convex problem, we employ null-space projection for beamforming design and adopt deep reinforcement learning for the trajectory optimization over a discrete-time scale. In each time slot, beamforming is optimized based on the channel state information to improve CRB performance while mitigating interference between the BS and the communication user. Simulation results demonstrate that the proposed method significantly reduces the time-averaged CRB by over 10%, compared with the ISAC system without UAV assistance, and also achieves a higher sensing accuracy than both the fixed-UAV-trajectory and the maximum-ratio-transmission-based beamforming benchmarks.
Read the full abstract
Combinatorics
Rational exponents near 3/2
Given a graph $H$, the extremal number $ex(n,H)$ is the maximum number of edges in an $n$-vertex graph not containing $H$ as a subgraph. The well-known rational exponents conjecture of Erdős and Simonovits states that for any rational $\gamma\in (1,2)$ there exists a single bipartite graph $H$ satisfying $ex(n,H)=\Theta(n^\gamma)$. Among other results, the conjecture has been verified for all $\gamma=1+a/b$, where $b>a^2$, by Jiang and Qiu and for all $\gamma=2-a/b$, where $b>\max\{a, (a-1)^2\}$, by Conlon and Janzer. In this paper, we establish the rational exponents conjecture for many $\gamma$ near the center of the interval, namely, for all $\gamma=1+\frac{rt-1}{2rt+2r}$, where $r,t$ are natural numbers satisfying $t\geq 2$, $r\geq 2t+3$.
Read the full abstract
Numerical Analysis · Cross-listed from cs.GR
A Splitting Architecture for Exact Reduced Coulomb Friction
Existing approaches to frictional contact dynamics typically either modify the Coulomb law to improve numerical robustness or solve the exact law in a fully coupled monolithic form. However, in its reduced form, exact Coulomb friction can be written as a cone complementarity problem with an augmented velocity, which reveals a natural split between a cone-constrained linear response and a scalar non-associated coupling induced by tangential velocity. We exploit this structure in the solver design. Our method uses an outer iteration to update the non-associated coupling explicitly, and an inner solve for a strongly convex cone-constrained quadratic program. This separation also makes the inner solver modular, so different numerical schemes can be used without changing the outer iteration. We evaluate the method on rigid-body benchmarks with stick-slip transitions and frictional stacking, and show that it reproduces exact Coulomb complementarity without smoothing or relaxing the friction law.
Read the full abstract
Algebraic Geometry
Maximum Likelihood Estimation on the Grassmannian of Lines
We study the positive Grassmannian through the lens of algebraic statistics. A closed formula is presented for the maximum likelihood degree of the Grassmannian of lines. We study the real and positive critical points for the Grassmannian of lines in 3-space.
Information Theory
Learning to Transmit: Volatility-Aware Predictive Communication for Energy-Efficient IoT Networks
Communication is the dominant source of energy consumption in Internet-of-Things (IoT) networks, yet many sensed measurements exhibit strong temporal correlations and provide little new information to the receiver. This paper introduces ADAPTIVEML, a volatility-aware predictive communication framework that enables IoT devices to intelligently decide when communication is necessary. Each sensor maintains a lightweight machine learning predictor and transmits only when the prediction residual exceeds an adaptive threshold proportional to the local signal volatility. By normalizing prediction errors using a rolling estimate of signal variability, the proposed transmission policy automatically adapts to changing environmental conditions, seasonal variations, and deployment-specific dynamics without manual threshold tuning. To address long-term non-stationarity, we further propose ADAPTIVEML-RLS, an online learning extension based on Recursive Least Squares (RLS) with exponential forgetting, allowing continuous adaptation to sensor drift and evolving signal characteristics. Extensive experiments are conducted on three heterogeneous real-world datasets comprising more than 2.4 million sensor observations from outdoor environmental monitoring, indoor wireless sensor networks, and urban air-quality sensing. Compared with six representative baselines, including periodic transmission, static-threshold suppression, ARIMA, Kalman filtering, EMA, and LMS filtering, ADAPTIVEML achieves up to 94.7% transmission reduction while maintaining a reconstruction error of 0.352$^\circ$C. ADAPTIVEML-RLS further reduces reconstruction error by 12–18% under drift conditions while preserving transmission reduction above 93%. These results demonstrate the effectiveness of volatility-aware predictive communication for energy-efficient and adaptive IoT networks.
Read the full abstract
Group Theory
On a function related to the Chermak-Delgado measure of a finite group
In this note, we give some conditions for commutativity and (super)solvability of a finite group $G$ based on a function related to the Chermak-Delgado measure $m(G)$. The possible values of this function are also investigated.
Category Theory
A Syntactic Approach to Ulmer's Bialgebras
Ulmer introduced a semantic notion of bialgebras that unifies a broad class of algebraic and coalgebraic structures. We develop a syntactic counterpart by introducing signature pairs $(\Sigma,\sigma)$ and bialgebraic theories $T$, providing a uniform language for constructing internal bialgebras in a $2$-categorical setting. For every bialgebraic theory $T$ and $\Sigma$-model $M$ within a $2$-category with PIE limits, we construct the object $M^T$ of internal $T$-bialgebras. Our approach to bialgebras admits a general Induced Functor of Algebras Theorem extending the classical lifting of lax monoidal functors to the categories of internal monoids. Since the construction of $M^T$ is expressed entirely in terms of PIE limits, accessibility, local presentability, orthogonal factorization systems, regularity, and exactness lift along the construction $M \mapsto M^T$ under suitable assumptions.
Read the full abstract
Probability
Density-Dependent McKean–Vlasov Diffusions: Subgaussian Occupancy Bounds and Polynomial Propagation of Chaos
We study the local density-dependent diffusion $dY_t=-\Xi(p_t(Y_t))\nabla\Phi(Y_t)\,dt+\sqrt2\,dW_t$ and a clipped, randomly shifted histogram particle approximation on $\mathbb{R}^d$. The central difficulty is that the empirical density is evaluated at the particles' locations and re-enters their drift, while the confining force $\nabla\Phi$ may be unbounded. We provide a path-space entropy proof under two verifiable analytic conditions: a uniform pointwise Gaussian envelope for the true density $p_t$, and a Gaussian–polynomial bound for its spatial gradient $\nabla p_t$. The potential is allowed to have a gradient of at most linear growth. The probabilistic input is a weighted exponential occupancy estimate under the independent product law. It is proved by Poissonizing the system at total intensity $N-1$, performing a one-cell leave-one-out estimate bounded via Poisson information, using Gaussian cell summability, and de-Poissonizing. For every fixed time horizon $T$, we obtain $\operatorname{Ent}(P_t^{N,k}|p_t^{\otimes k})\leq C_T k(h^2(1+|\log h|)+(h^{-d}+\log N)/N)$. Consequently, selecting the optimally balanced bandwidth $h\asymp (N\log N)^{-1/(d+2)}$ yields a total variation error of $\Vert P_t^{N,k}-p_t^{\otimes k}\Vert_{\operatorname{TV}}\leq C_T\sqrt{k}\,N^{-1/(d+2)}(\log N)^{d/[2(d+2)]}$ for fixed $k$. This includes the usual Ornstein–Uhlenbeck density and the density-dependent OU model whenever the PDE estimates hold on the considered interval. Furthermore, the histogram estimator offers a scalable approach for particle approximations. Using occupied-cell hashing, one algorithm step evaluates in expected $O(dLN)$ operations under standard constant-time hashing assumptions. For a fixed dimension and number of shifts, this requires expected $O(N)$ time, avoiding the $O(N^2)$ evaluation cost typical of standard kernel density estimators.
Read the full abstract
Optimization and Control
Data informativity for feedback rectification of switched linear systems
For switched linear systems, we leverage the concept of data informativity for data-driven feedback rectification to design a stabilizing controller. The established rectified eigenstructure means that all closed-loop subsystems share a common set of eigenvectors and enables us to compute a common quadratic Lyapunov function to ensure global uniform exponential stability of the closed-loop system. The proposed methodology is illustrated for two examples including a model of an unmanned aerial vehicle.
Read the full abstract
Information Theory
Bounds and Limitations on Codes Achieving List Recovery Capacity
In coding theory, list recoverability is a fundamental concept which robustly captures how “spread-out” codewords are in a code. More formally, given a code $C \subseteq \Sigma^n$ and input lists $S_1, \dots, S_n \subseteq \Sigma$ of size at most $\ell$, list recoverability requires that there are at most $L$ codewords $c \in C$ such that $c_i \in S_i$ for at least $(1-\rho)n$ choices of $i \in [n]$. List recovery is an important question which has found applications in many areas, including complexity theory, property testing, compressed sensing, streaming algorithms, and cryptography. As our first main result, we establish a tight “generalized singleton bound”. Formally, we show that for constant $\ell, L,\rho$ and sufficiently large alphabets $\Sigma$, if we define $R^*=\frac{L+1-\ell}{L}-\frac{L+1}{L}\rho$, it is possible for a $(\rho,\ell,L)$ list-recoverable code to have rate $R^*-\epsilon$ but impossible to have rate $R^*+\epsilon$. One direction of our result already directly generalizes and improves a weaker impossibility result due to Goldberg, Shangguan, and Tamo. For our second main result, we prove that there is a fundamental shortcoming in existing methods that aim to construct explicit, optimal list-recoverable codes. Indeed, recent work has constructed explicit codes achieving list-decoding capacity (along with other related properties) using a framework introduced in the work of Alon–Edmonds–Luby (AEL). We give a meta-analysis of such constructions by presenting an “AEL framework” which captures all such recent constructions in the literature. Within this framework, we show that no AEL-based code can break a recently-identified list-recovery barrier for additive and linear codes.
Read the full abstract
Metric Geometry
No Prime Tiling of an Isosceles Triangle
No isosceles triangle can be cut into a prime number (greater than three) of congruent triangles, and only an equilateral triangle can be cut into three congruent triangles.
Optimization and Control
Optimization models and algorithms for the Unit Commitment problem
The unit commitment problem determines the optimal strategy to meet the electricity demand at minimum cost by committing power generation units at each point of time. Solving the unit commitment problem gives rise to a challenging optimization problem due to its combinatorial complexity and potentially long solution time requirements. Our proposed solution approach utilizes a decomposition method in conjunction with alternative models from the EGRET library. Results of this decomposition approach tested against four benchmarking systems show that significant computational speed ups are achieved.
Read the full abstract
Mathematical Physics · Cross-listed from physics.flu-dyn
Cavitation Acoustic Perturbation Equations: A Computational Framework for Source-Resolved Multiphase Hydroacoustics
This work develops a cavitation-consistent acoustic perturbation framework for predicting sound generation and propagation in cavitating flows. Unlike conventional acoustic perturbation equations for single-phase or weakly compressible flows, the proposed formulation embeds cavitation physics directly into the acoustic equations. The cavitation acoustic perturbation equations (CAPE) incorporate vapor mass transfer, mixture compressibility, and pressure-rate effects within a unified formulation, allowing cavitation-induced noise sources to be resolved in the computational domain. The numerical framework is verified using one-dimensional wave-propagation problems. The solutions become insensitive to further mesh and time-step refinement, the perfectly matched layer suppresses boundary reflections, and the predicted attenuation over a range of source frequencies follows Stokes' sound attenuation law. The framework is then applied to cavitating flow past a circular cylinder and a NACA hydrofoil. The non-cavitating benchmark shows dipole-like radiation associated with unsteady loading, whereas cavitating cases exhibit monopole-like or geometry-modulated radiation caused by volumetric phase change. Source-term analyses identify tonal frequencies associated with vortex shedding, cavity shedding, and collapse-induced excitation. The phase-change terms provide a direct volumetric contribution to the monopole-like source, while localized collapse events appear through amplification of the pressure-rate source. The proposed framework extends acoustic perturbation methods to cavitating multiphase flows and provides an efficient tool for hydroacoustic prediction, source localization, and mechanism analysis in marine and hydraulic applications.
Read the full abstract
Information Theory
Extinction Depth and q-ary Error-Correcting Codes for the Limited Permutation Channel
In the radius-one limited permutation channel, errors consist of disjoint adjacent transpositions. A correcting code must separate distinct codewords: their error balls may not contain a common received word. Hamming distance does not ensure this, because disjoint swaps can make words differing in many positions confusable. For block-concatenation codes, earlier work tested each possible collision only through the longer initial block. This sufficient condition is not necessary: we exhibit a valid ternary block set that it does not certify. We introduce extinction depth, which tracks unresolved first-block pairs through later block extensions, and prove that their extinction at one common finite horizon certifies correction at every length. The criterion gives explicit $q=3$ and $q=5$ block sets with rates above $0.6777475$ and $0.6694926$. No block-set-independent depth bound exists: we give an exact linear family, verify a quadratic formula for every $3\leq k\leq 60$, and derive a polynomial-time finite-graph test. For growing alphabets, the normalized correction loss lies between $\ln\varphi$ and $\ln(1.82560995)$, while explicit finite-length covers improve finite-alphabet upper bounds. We develop a parallel directed-extinction theory for detection, including an all-length block criterion, a set unresolved by the earlier test, and exact corridor depth $4k$. Stable type lifting yields optimal first-order loss $\sqrt{2}$, and weak-zigzag codes improve the $q=3,4$ lower bounds. Finally, window restriction, cancellation, and a bounded pending-input frontier extend the criterion and its finite-state verification to every fixed displacement radius $r$.
Read the full abstract
Logic
Model Theory of Generic Vector Space Endomorphisms IV: Preservation of NATP
This paper further studies the model companion of an endomorphism acting on a vector space, possibly with extra structure. Let $T$ be a model-complete theory that $\varnothing$-defines an infinite $K$-vector space $\mathbb{V}$. In previous work, we introduced a family $\{T^C_\theta : C \in \mathcal{C}\}$ of extensions of the theory $T_\theta := T \cup \{\text{``$\theta$ is an endomorphism of $\mathbb{V}$''}\}$ that parameterizes all consistent extensions of the form $$ T_\theta \cup \left\{\sum\nolimits_{k}\bigcap\nolimits_{l}\operatorname{Ker}(\rho_{j, k, l}[\theta]) = \sum\nolimits_{k}\bigcap\nolimits_{l} \operatorname{Ker}(\eta_{j, k, l}[\theta]) : j \in \mathcal{J}\right\}, $$ where all sums and intersections are finite, all the $\rho[\theta]$'s and $\eta[\theta]$'s are polynomials over $K$ with $\theta$ plugged in, and $\mathcal{J}$ is some possibly infinite index set. We also presented a sufficient condition that implies that every $T^C_\theta$ has a model companion $T\theta^C$. In this paper, we show that, under this sufficient condition, the model companion $T\theta^C$ has $\operatorname{NATP}$, a neostability property recently introduced by Ahn and Kim, whenever $T$ does.
Read the full abstract
Category Theory
Twisted Deligne products of semisimple tensor categories
We discuss the classification of twisted Deligne products of two semisimple tensor categories $\mathcal C,\mathcal D$, i.e., categorifications of the tensor product of their Grothendieck rings in which the factors are categorified by $\mathcal C$ and $\mathcal D$. In particular, we show that if both factors have no non-trivial gradings, or if one factor has neither non-trivial gradings nor tensor structures on the identity functor, then the only twisted Deligne product is the ordinary one. Using the work arXiv:2405.10207 by Müller, Peña Pollastri and Plavnik, this gives, in principle, a group-theoretical classification of twisted Deligne products and, more generally, exact factorizations of arbitrary fusion categories. In the Appendix we introduce the notion of categorical $n$-cocycles for $n=2,3,4$ and show that they are all pullbacks of group $n$-cocycles from the universal grading group of the underlying based ring. In the case of $4$-cocycles, this answers a question of Johnson-Freyd, Ostrik and Yu from arXiv:2601.09060.
Read the full abstract
Information Theory
Combinatorial Capacity Bounds for the $q$-ary Deletion Channel
We study the \(q\)-ary deletion channel via the pattern-count scalar \(N_n(x,y)\), the number of deletion subsets mapping \(x\in\Sigma_q^n\) to \(y\in\Sigma_q^k\), which factorizes the transition probability. Two sum identities on \(N_n\) certify stochastic normalization and, under uniform input, yield an exact closed-form output entropy. These give the finite-block capacity sandwich \( (1-d)\log_2 q-h_2(d)\;\le\; C_{q,n}\;\le\;(1-d)\log_2 q. \) The exact uniform-input rate is \( \frac{1}{n}I_U(X;Y) =(1-d)\log_2 q+\frac{1}{n}H_{\mathrm{Bin}}(n,1-d)-h_2(d)+\frac{\Delta_n(d)}{n}, \) from which the simpler certified bound \( C_{q,n}\ge (1-d)\log_2 q-h_2(d)+\frac{\Delta_n(d)}{n} \) follows. The small-\(d\) bound \(C_q(d)\ge\log_2 q+d\log_2 d+O(d)\) follows for all \(q\ge 2\). Numerical experiments at \(n=3,5,10\) and \(q=2,3\) confirm all bounds.