2026-10-20 | | Total: 779
Let $0<γ_1\leq γ_2\leq \cdots$ be the ordinates of the non-trivial zeros of the Riemann zeta function in the upper half-plane, counted with multiplicity, and fix $c>0$. We say that a non-trivial zero $ρ_n=β_n+iγ_n$ is $c$-isolated if $\min(γ_{n+1}-γ_n, γ_n-γ_{n-1})>2πc/\log γ_n$. Note that any isolated zero is necessarily simple and on the critical line. We prove unconditionally that more than $50.87\%$ half of the non-trivial zeros are $1/6$-isolated. Previously, it was not known unconditionally that there are infinitely many isolated zeros. By allowing $c$ to be extremely small, we can push the above proportion to $13/20$. We also show that there are positive proportions of $k$ consecutive ``moderate'' gaps (gaps larger than a constant times the average local gap) between the ordinates of the zeta zeros, for $k=3, 4, 5$. As an application we combine our work on isolated zeros with an unconditional refinement of the Bombieri--Hejhal method to show unconditionally that at least $31\%$ of the zeros of a linear combination of two distinct primitive Dirichlet $L$-functions with the same conductor and parity lie on the critical line. The previous unconditional record for this proportion was $7.38\times 10^{-9}$, which was obtained by Dousselin who refined and made explicit a result of Selberg.
Bilevel optimization provides a framework for machine learning tasks such as hyperparameter optimization, data reweighting, and meta-learning. Existing Newton-based bilevel methods typically rely on the classical hypergradient formula, which is generally unavailable when lower-level solutions are nonunique and the lower-level Hessian is singular. We develop two Newton-based methods beyond lower-level strong convexity by exploiting curvature shared by the lower-level and multiplier updates. Newton-based GALET handles potentially singular lower-level curvature under a positive spectral-gap condition, exploiting its range structure to avoid the squared-Hessian multiplier system used in GALET. Newton-based BAMM uses a strongly convex aggregation model when the lower-level objective is convex and the upper-level objective is strongly convex in the lower-level variable. We establish convergence guarantees and improved complexity bounds relative to GALET and sl-BAMM under the corresponding assumptions. We further introduce GNBO, a generalized Newton-based framework that captures the shared-curvature structure of the proposed methods and several existing Newton and quasi-Newton bilevel methods. Experiments on synthetic problems, data hyper-cleaning, and GPT-2 data reweighting evaluate the proposed methods across different bilevel learning tasks.
In [Liu26a], varying the length and moment parameters in the cubic-theta argument of [Ope26a] gave the zero-free half-plane $Re(s)>B_{\mathrm{new}}=0.874957019420098946\ldots$ for finite-order Hecke $L$-functions over $\mathbb Q(\sqrt{-3})$ and for Dirichlet $L$-functions, and its optimality in the specified exponent comparison system. In this paper, we extend the argument to finite-order Hecke $L$-functions over cyclotomic fields containing $μ_{12}$ and, by class-field descent, to every pair $(K,η)$ whose cyclic class field has metabelian normal closure over $\mathbb Q$, including every finite-order character over every abelian number field. Then, adapting the averaged numerator estimate developed by Xu and Sang [XS26] for Hecke $L$-functions over $\mathbb Q(\sqrt{-3})$ and combining it with a limiting contradiction argument, we obtain the boundary $B_\flat=0.874886647039964\ldots$, the root near $0.8749$ of $133488B^4-376596B^3+384534B^2-167541B+26231$; it is the exact infimum of the refined comparison system, unattained by finite chains starting above it. Since the analytic argument uses only $μ_{12}\subset F$, we run it over $F=K(ζ_{12})$ for an arbitrary number field $K$ and descend: every finite-order Hecke $L$-function over every number field is zero-free for $Re(s)>B_\flat$. The analytic argument uses the exceptional-theta constructions and local calculations of the preprint [Ope26b].
Sontag conjectured in 1988 that global controllability of polynomial control systems with rational coefficients is undecidable, but that the controllable systems form a recursively enumerable set. We prove the first part of the conjecture and disprove the second. For every homogeneous polynomial vector field $F$ of odd degree with rational coefficients, we construct an explicit polynomial control system $Σ_F$, with one more state than $F$, that is globally controllable if and only if the origin is \emph{not} globally asymptotically stable for $\dot x=F(x)$. By a recent result, the homogeneous polynomial vector fields with a globally asymptotically stable origin form a set that is recursively enumerable but not decidable. Consequently, global controllability of polynomial systems is undecidable, and the controllable polynomial systems are not recursively enumerable. The same conclusions hold for systems of the form $\dotξ=f(ξ)+Bu$ with a polynomial drift and a constant input matrix.
We study bending-minimizing sequences of closed curves in a fixed knot class and the self-contact of their limits. A positive matrix-valued measure records a possible loss of curvature energy. Ambient variations, localized on arbitrary coincidence sets, force this measure to vanish and give strong convergence in $H^2$. We then derive the elastic equation on isolated strands and overlap bundles, describe the distribution that remains at contact, and prove balance and rigidity statements under their precise additional hypotheses. Equality in the bridge bound gives constant curvature. Combining strong compactness with quantitative competitors gives an alternative proof that elastic limits in classes with equal bridge and braid indices are covered circles. An independent recovery theorem identifies the strong $H^2$ closure with the regular $H^2$ part of the $C^1$ closure and also implies strong compactness. Finally, a whole-line contact model illustrates which estimates are needed for finer shape selection.
On the real hyperbolic space $\mathbb{H}^n$, the classical Muckenhoupt condition does not describe the weights $w$ for which the Hardy-Littlewood maximal operator $M$ is bounded on $L^p(w)$. We show that, for $1<p<\infty$, $M$ is bounded on $L^p(w)$ if and only if $w$ is a local $A_p$ weight and satisfies a two-layer condition, which compares $w$ and its dual weight on two layers of the shadow of a point. We also treat the weak type case, radial weights and the range $p\le1$, and we prove a weighted Herz principle for radial convolution operators.
We determine the dimension-independent first-order query complexity of minimizing a Lipschitz geodesically convex (g-convex) function on hyperbolic space. Given a ball of radius $r$ containing a global minimizer, we prove that $Ω(ζ_r/ε^2)$ queries are necessary to attain accuracy $εMr$, where $M$ is the Lipschitz constant and $ζ_r=r/\tanh r$ captures the effect of curvature. This matches the known $O(ζ_r/ε^2)$ upper bound for the Riemannian subgradient method (Zhang and Sra, 2016). This is the first query-complexity result for g-convex optimization that is tight simultaneously in curvature and target accuracy. The standard Euclidean recipe of constructing hard instances as maxima of affine functions does not transfer to hyperbolic space, so the lower bound requires a new mechanism. Our starting point is the worst-function construction of Criscitiello and Boumal (2023, Appendix A). It gives the same lower bound for a subclass of algorithms, including the subgradient method, whose queries remain in prescribed g-convex regions. We retain this geometric setup and carefully tune and augment their hard function to accommodate every deterministic first-order algorithm. The key new mechanism is an advance-or-stall scheme: when a query lies in these regions, the oracle reveals at most one new information direction (advance); when a query leaves them, we show that the oracle can answer without revealing any new hidden direction (stall). Unlike earlier curvature-dependent lower bounds, the proof uses no ball-packing argument.
In this paper, we generalize the code basis reduction to non-Hamming metrics. We propose a block basis reduction (BKZ) algorithm for codes that works for all coordinate-additive metrics. Our algorithm is a generalization of the algorithm of Ghentiyala and Stephens-Davidowitz, which in turn is an adaptation of the famous BKZ algorithm for lattices to linear codes and the Hamming metric. For block size~$2$ (LLL), we propose a polynomial-time algorithm with complexity~$\mathcal{O}(n^4 q^3)$ operations in~$\mathbb{F}_q$, where~$n$ is the length of the code and~$q$ is the alphabet size. For the Lee metric we provide an implementation of our algorithms and present experimental results on finding low Lee weight codewords.
We construct a new type of congruence from irregular weight mod $p$ Picard modular forms (for $p$ inert) to ones of regular weight, analogous to the one obtained by Kansal--Levin--Savitt in the Hilbert setting. The congruence is obtained by using a geometric Jacquet--Langlands phenomenon, as first established by Vollaard, to transfer irregular weight mod $p$ Picard modular forms from the Picard modular surface to a Shimura set. We call this an exotic congruence, since it relates mod $p$ automorphic forms associated to different groups. This congruence will be used in upcoming work of the author proving mod $p$ local-global compatibility for irregular weight Picard modular forms (in the $p$ inert case).
We study directional expansiveness of the chair tiling. We show that for every direction $V$ and every $r>0$, there exist distinct chair tilings that agree on $\mathbb{Z}^2\cap V^r$. Hence the chair tiling has no expansive direction.
We prove uniform convergence of the canonical discrete uniformization maps for every sequence of equilateral polygonal disks whose mesh tends to zero and whose marked Jordan boundaries converge uniformly. The limit is the normalized Riemann map, with convergence on the entire closed reference triangle. This proves Conjecture 7.3 of Luo, Sun and Wu for the original construction with Delaunay flips. The proof combines closed-domain compactness, marked lattice limits, a packing and resistance estimate for logarithmic cusp scales, and an envelope comparison principle, followed by a local affine energy argument.
We introduce a metric geometry on the space of positive, linearly growing quasimorphisms on a mixing subshift of finite type. Two such quasimorphisms are identified when their homogenizations agree up to multiplicative scaling. The resulting moduli space carries a natural symmetrised Thurston-type metric, defined from the extremal ratios of their periodic orbit lengths. We prove that this is a genuine metric and determine several global features of the resulting space: it is unbounded, separable, incomplete, and has geodesics that are canonically generated by thermodynamic pressure functions. We then describe the endpoints of these geodesics and develop a boundary theory for our moduli space. We apply this theory to metric structures on right-angled Artin and Coxeter groups. In particular we construct isometric embeddings of geometric moduli spaces into our quasimorphism moduli spaces. As a consequence, Manhattan curves for metrics on these groups are $C^1$, with equality cases characterised by rough similarity of the underlying metrics. This provides insights into Manhattan curves in settings where neither classical thermodynamic formalism for Hölder potentials nor Patterson-Sullivan theory apply. The same construction recovers the corresponding results for hyperbolic groups and places these examples within a common dynamical-geometric framework.
The Knizhnik--Zamolodchikov (KZ) connection is defined on the constant bundle determined by $n$ finite-dimensional irreducible representations of a simple Lie algebra $\mathfrak{g}$ and a nonzero complex parameter $κ$, over the configuration space of $n$ distinct points on the affine line. Work of Schechtman--Varchenko and Looijenga shows that, when $κ$ is rational, the resulting local systems admit natural motivic realizations, called KZ motives. The Kazhdan--Lusztig equivalence identifies the local systems underlying these motivic variations with certain representations of the pure braid group arising from the corresponding Lusztig quantum group at a root of unity. We formulate conjectures relating the weight filtrations of KZ motives to structures in quantum-group representation theory, particularly tilting modules and the Jantzen and Andersen filtrations, and provide evidence for these conjectures. Our methods use a new deformation of the representation ring that records the first moments of the weight filtration, together with Schmid's theorem on limiting mixed Hodge structures and results on the Hodge filtration obtained in earlier joint work with S. Mukhopadhyay.
We derive a canonical functional model for a pair of isometric operators $(V_1, V_2)$ acting on a Hilbert space $\mathcal{H}$ that satisfy the automorphic covariance relation $V_1 V_2 = V_2 \varphi(V_1)$, where $\varphi$ is an automorphism of the unit disk $\mathbb{D}$. This extends the classical Bercovici--Douglas--Foias functional model for commuting isometric pairs to a dynamically twisted setting. When both $V_1$ and $V_2$ are pure isometries (shifts), the model simplifies to $(M_z, M_Θ W_\varphi^*)$ acting on $H^2(\mathcal{F})$, where $W_\varphi$ is a unitary weighted composition operator on the Hardy space $H^2$ and $Θ\in H^\infty(\mathbb{D}, \mathcal{B}(\mathcal{F}))$ is an operator-valued inner function such that $M_Θ$ is a shift. As an application, we establish a complete unitary classification theorem for automorphically covariant isometric pairs by identifying a complete set of joint unitary invariants. Notably, our results are novel even in the special case of rotation automorphisms ($V_1 V_2 = q V_2 V_1$ with $|q| = 1$).
We exhibit the $p$-complete $\mathbb{C}$-motivic Eilenberg-Mac Lane spectra $\mathrm{M}\mathbb{F}_p$ and $(\mathrm{M}\mathbb{Z})_p^{\wedge}$ as $\mathbb{E_3}$-$\mathrm{MGL}_p^{\wedge}$-motivic Thom spectra over the motivic $\mathbb{G}_m$-loop spaces $Ω_{\mathbb{G}_m}\mathrm{BGL}$ and $Ω_{\mathbb{G}_m}\mathrm{BSL}$, respectively. Furthermore, for $p=2$, we compute the $\mathrm{MGL}_2^{\wedge}$-relative $\mathbb{G}_m$-Hochschild Homology $\mathbb{G}_m \odot^{\mathrm{MGL}_2^{\wedge}} {-}$, recently introduced by Bachmann-Burklund-Dundas-Østvær, of $\mathrm{M}\mathbb{F}_2$ and $(\mathrm{M}\mathbb{Z})_2^{\wedge}$ to be $\mathrm{M}\mathbb{F}_2[\mathrm{BGL}]$ and $(\mathrm{M}\mathbb{Z})_2^{\wedge}[\mathrm{BSL}]$ as $\mathbb{E}_2$-$\mathrm{M}\mathbb{F}_2$-algebras and $\mathbb{E}_2$-$(\mathrm{M}\mathbb{Z})_2^{\wedge}$-algebras, respectively.
We study differential equations -- not necessarily monic -- for Jacobi modular forms. As an application, we derive a second-order differential equation for the elliptic genus of a $K3$ surface involving the logarithmic derivative of the modular invariant $j(τ)$, and we explicitly find a third-order equation with a leading coefficient of weight 12 for the elliptic genus of any four-dimensional Calabi--Yau variety. Furthermore, we propose a simple recursive algorithm for investigating second-order differential equations with leading coefficient $E_4(τ)$ or $E_6(τ)$.
Let $F$ be a non-archimedean local field. For any irreducible component $C$ of a nilpotent variety of type $A$ let $π(C)$ be the irreducible representation of a general linear group over $F$ attached to it. We show that if the multisegment of $C$ is multiplicity free, $C$ is rigid if and only if $π(C)$ is $\square$-irreducible. In this case, \[ π(C*D)=soc\bigl(π(C)\timesπ(D)\bigr),\qquad π(D*C)=soc\bigl(π(D)\timesπ(C)\bigr) \] for every irreducible component $D$, where $*$ is the binary operation defined by Aizenbud and Lapid, in accordance with a conjecture of Lapid and Mínguez. We also give a computer-assisted example in $A_7$ of a non-simple rigid module $z\ne0$ which does not admit a proper rigid submodule $y\ne0$ such that $Ext_Π^1(z,z/y)=0$, answering negatively a question of Aizenbud and Lapid.
In 1991, Erdős, Faudree, Rousseau and Schelp asked to estimate the smallest $u=u(n,d,s)$ such that every $n$-vertex graph with average degree at least $d$ has a subgraph on at most $u$ vertices with minimum degree at least $s$. We resolve this problem up to polylogarithmic factors for all even values of $s$. More precisely, we prove that every $n$-vertex graph with minimum degree $d\geq s$ has a subgraph on $O_s\left(\max\left\{\frac{n}{d^{s/(s-2)}}(\log n)^{180},\log n\right\}\right)$ vertices with minimum degree at least $s$, whenever $s\ge 4$ is even. This improves a result of Janzer, Sudakov and Tomon, who obtained a similar conclusion for subgraphs of average degree at least $s$. We also prove that for every positive integer $s$ and every $\varepsilon>0$, there exists a finite family $\mathcal{H}$ of graphs with minimum degree $s$ such that $\operatorname{ex}(n,\mathcal{H})=O(n^{2-2/s+\varepsilon})$. This is tight and is the first result of similar kind that applies to odd values of $s$ greater than 3.
We study an inverse problem for the variable coefficient Moore Gibson Thompson (MGT) equation in a bounded domain, where the objective is to recover an unknown initial acceleration from boundary measurements of the acoustic field. We establish both injectivity and stability of the associated measurement operator. Injectivity is established through a unique continuation property by applying a theorem of Tataru for analytically principally normal operators. The argument relies on a factorization of the third order operator into wave type components, combined with Carleman estimates and propagation techniques in the spirit of Tataru. Stability is derived through a microlocal analysis of the forward operator, based on the construction of geometric optics approximate solutions and the identification of an elliptic pseudodifferential operator arising from microlocal back propagation of the data. These results connect the MGT inverse problem with the classical theory of hyperbolic inverse problems, such as photoacoustic tomography, while highlighting new features induced by the third order dynamics and damping.
We prove a square root bound for the exceptional sets in the binary Goldbach and prime difference problems under quasi Riemann Hypothesis. If every nontrivial zero of every primitive Dirichlet $L$-function has real part at most a fixed $1/2\leqΘ<11/12$, then $E_+(N)\ll_ΘN^{1/2}(\log N)^{100}$ and $E_-(N)\ll_ΘN^{1/2}(\log N)^{100}$. Here $E_+(N)$ counts the even integers greater than $2$ and at most $N$ that are not sums of two primes, and $E_-(N)$ counts the positive even integers at most $N$ that are not differences of two primes. The common analytic estimate is a fourth moment bound of the minor arcs for a smoothed prime exponential sum.
The water removal from refractory castables is a complex process involving heat and mass transport and thermohygromechanical stresses. In order to reduce the risk of explosive spalling by improving the heat-up curves, numerical models are critical to predict the castable behavior during drying. However, most of the existing models either neglect the thermomechanical stresses or consider them as an \textit{a posteriori} assessment of the likelihood of failure. The present work proposes an open-source thermohygromechanical phase-field model to simulate the drying of refractory castables and predict the occurrence of explosive spalling. The model is based on an adaptive two-scale time-stepping strategy that allows for a significant reduction of the computational cost associated with the simulation. The model was validated against literature results and qualitative comparisons with experimental data were carried out to assess its ability to predict the failure behavior of partially saturated porous media during heating. The numerical simulations were also applied to investigate the impact of different heating methods on the risk of failure during drying. It was found that even at the significantly lower heating rates typical of refractory castable drying - compared to the fire accidents that motivate most thermomechanical studies on concrete structures - thermomechanical stresses still play a significant role in the explosive spalling, a finding well established for concrete under fire but novel for refractory castables. As a consequence, it was also concluded that the use of heterogeneous heating methods can increase the likelihood of explosive spalling. Ultimately, feasible solutions to increase the drying efficiency and safety that are compatible with the industrial setting were proposed and a free, open-source implementation of the model is provided.
We show that gradient descent (GD) with step size $η$ can achieve $\varepsilon$ optimality gap on the multiclass logistic regression loss in a number of iterations $T$ that is independent of $\varepsilon$. To do so, the step size $η$ must be chosen large enough as a function of $\varepsilon$ (and properties of the linearly separable dataset), while $T$ is a constant that only depends on the data. The relevant properties of the data include the number of points and their margin (for $T$), and one additional quantity (for $η$) that turns out to be necessary in dimension $d \geq 3$. The proofs build on the known fact that, as $η\to \infty$, suitably rescaled GD becomes the perceptron algorithm, and the latter finds a perfect classifier in $T$ iterations. A similar claim was shown for two classes in dimension $d = 2$ using other proof techniques. This result adds to previous works that showed GD on the logistic regression loss can be accelerated by taking large step sizes. Those provoke a wild, oscillatory initial phase, followed by a stable phase. We provide results of this (qualitatively different) nature, too.
We present a boundary integral method for evaluating electron repulsion integrals (ERIs) in dielectric environments of piecewise constant permittivity. This generalizes conventional bare-Coulomb ERIs by incorporating screening from structured dielectric surroundings, such as the substrates and/or encapsulating layers present in modern two-dimensional-material devices. We formulate the problem as a second-kind boundary integral equation for the induced surface charge on the dielectric interfaces, and discretize it with a Nyström method on adaptively refined panels. The incident potential generated by each smooth orbital source density is evaluated by a Fourier truncated-kernel method in or near to its source box, and by fast multipole method (FMM) accelerated quadrature farther away. The linear system is solved by a Krylov iteration also accelerated by the FMM. We focus on cuboidal device geometries, and retain efficiency at high aspect ratios via close-to-touching panel quadratures. We further develop a distributed pipeline that computes the full ERI tensor, comprising thousands of orbital pairs, embedded in a single dielectric geometry. We demonstrate the method on an effective three-layer graphene slab over a Si/SiO$_2$ heterojunction, focusing on the central graphene layer represented by $198$ $p_z$-like orbitals, retaining the $2613^2$ non-negligible entries of the $198^4$-sized ERI tensor. The model is resolved on $1.3$ million surface unknowns, achieving a relative error of $0.004\%$, at a cost of $22$ node-hours.
Let $A\subset\mathbb{R}^n$ be a nonempty compact set and let $K\subset\mathbb{R}^n$ be a nonempty compact convex set. We study the problem of when there exists $λ>0$ such that the Minkowski sum $A+λK$ is convex. First, we show that this finite convexification is equivalent to a uniform covering condition on the exposed sets of $A$ and $K$. We then show that another characterization is possible in terms of the covering radii of the iterated exposed sets after translating each exposed face of $K$ so that its relative centroid is $0$. In the case of convex polytopes, we show that it suffices to check the exposed sets in the first level of the iteration. We also study several special cases, such as strictly convex convexifiers, lower-dimensional convexifiers, strongly isomorphic polytopes, and the planar case.
We present FOMTools, an open-source Python package for developing, analyzing, and deploying first-order algorithms for parametric linear and convex quadratic optimization. The user specifies an iteration in Python, with syntax close to its mathematical description. The package represents the iteration as a computation graph of affine and piecewise affine operations. To analyze worst-case performance, FOMTools formulates a mixed-integer verification problem that maximizes a convergence residual over a specified set of problem instances and initial iterates. For probabilistic performance analysis, a batched GPU simulator collects algorithm trajectories and computes statistical bounds on the probability of exceeding a specified convergence tolerance. From the same algorithm description, FOMTools generates specialized solver code for embedded hardware, multicore CPUs, and GPUs. The package first generates the numerical data, factorizations, and matrix-vector routines without any agent. A coding agent then translates the iteration into C or CUDA code and iteratively optimizes the implementation. After each phase, the package runs the generated code on sampled inputs and compares it with a reference execution. FOMTools also serves as a CVXPYgen backend, generating and optimizing solvers for parametric CVXPY problems. We demonstrate FOMTools on sparse coding, network optimization, optimal control, and portfolio rebalancing. The verification bounds reveal worst-case behavior that sampling misses. The generated solvers run up to two orders of magnitude faster than NumPy implementations, and they outperform hand-written GPU code and the OSQP solvers generated by CVXPYgen.