Papers

In this survey article, we give an introduction to the relatively young subject of $2$-representation theory, which studies categorifications of important players in classical representation theory. In particular, we provide a streamlined exposition of the main results leading to the classification of simple $2$-representations for Soergel bimodules associated to finite Weyl groups in characteristic $0$, obtained in [M. Mackaay, V. Mazorchuk, V. Miemietz, D. Tubbenhauer and X. Zhang, Simple transitive $2$-representations of Soergel bimodules for finite Coxeter types. Proc. London Math. Soc. 126 (2023)].

We explain the relation between the relative topological Hochschild homology $\THH(R/\Sphere_{W(k)}[x_0,\ldots,x_n])$ and the Nygaard-completed Frobenius-twisted relative prismatic cohomology $\widehat{\Prism}^{(1)}_{R/W(k)[x_0,\ldots,x_n]^\wedge}$, where $W(k)[x_0,x_1,\ldots,x_n]\rightarrow R$ is relatively quasiregular semiperfectoid. As an application, for $R=\Z_p[x]/(px)$, we compute $\pi_*\THH(R)^\wedge_p$ by descent along $\THH(R)^\wedge_p\rightarrow \THH(R/\Sphere_p[z,x])$, where $R$ is regarded as an $\Einfty$-$\Sphere_p[z,x]$-algebra via $\Sphere_p[z,x]\xrightarrow{z\mapsto p,x\mapsto x}\Z_p[x]/(px)$.

Computing eigenvalues and eigenvectors is at the heart of the solution of many non-linear problems. For instance, finding the roots of polynomial systems reduces to computing joint eigenvectors of operators of multiplication. Similarly, tensor decomposition can be performed via the joint diagonalization of submatrices of the Catalecticant of the tensor. We describe and illustrate symbolic-numeric methods for computing the solutions of these algebraic problems from the computation of joint eigenvectors of commuting operators, and for analysing their multiplicity structure, as well as their implementation in the package AlgebraicSolvers.jl.

Reliability set estimation identifies input regions where a response probability exceeds a target level, bridging estimation and safety-critical decisions. Practitioners typically start with a working model, an imperfect approximation of the true response surface. Relying on this imperfect model may incur decision risk, potentially certifying unsafe regions as safe. We develop a unified framework that turns such a working model into a trustworthy decision rule. First, a modeling-then-calibration procedure decouples estimation from decision. Since the true set is unobservable, we introduce an asymmetric, observable surrogate loss and use a separate calibration set to select a bias-correcting threshold, reducing decision risk and achieving $O_P(1/n)$ volume convergence. Second, we leverage conformal risk control with the surrogate loss to control false inclusion risk, which is the most safety-critical error, at a pre-specified level regardless of working model quality. Together, these calibration procedures show that a separate calibration set is necessary for risk control. Third, an adaptive design concentrates observations on the reliability set and its boundary, improving model quality where errors most affect decisions while controlling budget elsewhere. Numerical studies show not only more accurate set estimates but also calibrated finite-sample risk control that classical plug-in methods lack.

Freidlin–Wentzell theory converts weak-noise large deviations into a Hamiltonian variational problem. We study how a $k$-dimensional normally hyperbolic invariant manifold (NHIM) $N$ of the deterministic dynamics appears in this Hamiltonian system. Its zero-momentum copy $N_0=N\times\{0\}$ is invariant, but the Hamiltonian dynamics has $2k$ center directions near $N_0$: $k$ tangent to $N$ and $k$ conjugate covector directions. We construct the resulting local symplectic geometry and show that fluctuation extremals approaching $N_0$ backward in time at the strong normal rate form an $n$-dimensional exact Lagrangian invariant manifold carrying a single-valued action. A counterexample shows that a corresponding zero-energy section need not be normally hyperbolic within $H_{\mathrm{FW}}^{-1}(0)$. We then consider reaction dynamics. If a parameter moves a deterministic trajectory toward a codimension-one reactivity boundary, the minimum Freidlin–Wentzell action required to reach the boundary is quadratic in the distance from the threshold parameter. Its coefficient is determined by the relative motion of trajectory and boundary and by how effectively the available noise acts transversely. This deterministic boundary is distinct from a noise-dependent stochastic transition state or a committor surface. In a solvent–solute model, varying solvent mass moves the phase-space reactivity boundary while leaving the potential-energy surface fixed, changing the rare-event cost without changing the potential-energy barrier.

We consider the natural family of Fourier inequalities for the Paley–Wiener space $\mathrm{PW}^q(\Omega)$, consisting of $L^q$-functions with Fourier support in a convex set $\Omega \subset \mathbb{R}^n$, $n \geq 2$, free of affine lines. Namely, \[ \int_{\Omega}\dfrac{|\hat{f}(x)|^p}{\omega_{\Omega}^d(x)}dx\leq C\|f\|_{L^q}^p,\quad f\in \mathrm{PW}^q(\Omega). \] Here $\hat{f}$ is the Fourier transform of $f$, $1 \leq p, q < \infty$, $d \in \mathbb{R}$, and $\omega_\Omega$ is the frequency multiplier weight associated with the Paley–Wiener space of $\Omega$, \[ \omega_{\Omega}(x)=m(\Omega\cap (2x-\Omega)), \qquad x \in \Omega. \] For an arbitrary polyhedron $P$, we completely characterize the triples $(p,q,d)$ which yield valid Fourier inequalities. For a ball $B$, we characterize the valid triples when $p \geq 2$. When $p < 2$, the situation is different for the ball, and natural critical inequalities fail. However, we show that the spherical restriction conjecture implies a family of subcritical Fourier inequalities for the ball, which in turn imply the Kakeya conjecture (in its Minkowski-form). Finally, we link our family of Fourier inequalities to the theory of truncated Hankel operators acting on the Paley–Wiener space of $\Omega$.

Cells for the canonical and $p$-canonical bases organise the representation theory and geometry of Hecke categories. We determine the relevant $p$-canonical basis elements and the resulting $p$-cell structure for the subregular and submaximal cells in all classical types.

Shallow networks with prescribed or randomly sampled hidden parameters are widely used as numerical trial spaces, yet their optimal Sobolev approximation power with standard smooth sigmoidal activations in general dimension remains unresolved. We establish the corresponding optimal rates for a class of smooth sigmoidal activations with Schwartz-class derivative decay, including $\tanh$, the logistic sigmoid, and the error function $erf$. We first construct deterministic direction–offset dictionaries with $M$ features such that every $u\in H^k(\Omega)$ can be approximated with error of order $M^{-(k-m)/d}$ in $H^m(\Omega)$ for all $0\le m\le k$. This rate is optimal in the sense of Kolmogorov widths for Sobolev balls. We further prove that dictionaries obtained by independent parameter sampling from any prescribed density bounded away from zero attain the same approximation exponent with high probability, up to logarithmic oversampling. The analysis develops a sigmoidal ridge representation and combines it with deterministic or probabilistic quadrature in direction–offset space while retaining polynomial control of the output coefficients. Numerical experiments across a broad range of dimensions, target regularities, and Sobolev error norms recover the predicted algebraic rates for both deterministic and random feature dictionaries.

Open-porous materials exhibit pronounced compressibility, nonlinear densification, and power-law scaling of stiffness with density. In this work, we propose a thermodynamically consistent compressible neo-Hookean constitutive model for open-porous solids in which porosity serves as the primary governing variable. The strain-energy density is formulated to couple distortional elasticity of the solid skeleton with a volumetric response governed by deformation-induced porosity evolution, including a bounded representation of pore collapse. The formulation introduces a minimal set of parameters, namely the initial porosity, intrinsic skeleton moduli, and a scalar parameter controlling the onset of densification. A key feature of the model is a modified volumetric term in which the response is normalised by the current porosity, ensuring a physically consistent transition from a porous to a densified state without artificial stiffening. In the small-strain limit, the model recovers classical linear elasticity with effective moduli that may be chosen either from homogenisation bounds, such as the Hashin-Shtrikman estimates, or from Gibson-Ashby-type power-law scaling to capture topology-dependent behaviour. At finite strains, the formulation captures the characteristic nonlinear stiffening and convex stress-stretch response associated with progressive pore collapse. The proposed framework thus provides a compact, flexible, and extensible constitutive description that unifies effective-medium consistency with experimentally observed scaling behaviour, and is well suited for finite element implementation and multiscale modelling of highly compressible open-porous materials. The model is finally validated against available experimental data.

Suppose that $p$ is a prime and $X$ is a finite group with a strongly $p$-embedded subgroup, for example a rank $1$ group of Lie type in characteristic $p$. Let $m$ denote the $p$-rank of $X$ and assume that $m \ge 2$. We say that a faithful $\FF_p X$-module is $k$-active if some element of order $p$ in $X$ acts with exactly $k$ non-trivial Jordan blocks. In this paper, we determine the non-trivial composition factors of the faithful $\FF_pX$-modules which are $k$-active for some $k\leq m$.

Discovering dynamical systems from trajectory data is a central problem in applied mathematics and engineering. Whilst recent advances in machine learning have led to strong progress in data-driven system identification, much less attention has been given to systems with discontinuous dynamics. These systems are nevertheless highly relevant in applications, including climate dynamics and mechanical systems with friction. In this work, we consider the problem of identifying piecewise-smooth dynamical systems directly from trajectory data. Compared with the smooth setting, this requires recovering the governing equations and detecting the switching hyperplanes that separate different dynamical regimes and characterising their behaviour, such as sliding motion. We present a modular framework for discovering such systems by first estimating switching hyperplanes from data and then learning smooth dynamics within each region using geometry-constrained neural networks. The geometry-learning phase is studied from a statistical perspective, analysing the identifiability of the discontinuities and the robustness of the procedure. We also introduce a novel neural network architecture with a prescribed discontinuity set, and provide a theoretical analysis of its approximation properties. The approach is tested on low-dimensional benchmark problems, including dry-friction oscillators and the PP04 climate model for the ice ages.

This paper establishes a new upper bound for the irrationality measure of arctan(1/2), namely 8.585166... Following the strategy and proof techniques of Zudilin and Zeilberger (who improved the bound for pi), we use a family of complex contour integrals with symmetric integrands. The integrality and divisibility of the partial-fraction coefficients are derived via p-adic valuations and a specially defined prime set P_n. Combined with saddle-point asymptotics, the growth rates of the integral sequence and its rational/logarithmic components yield the desired irrationality measure.

Let $G$ be a connected bridgeless graph. In 2011, Kung and Royle showed that all roots of the flow polynomial $F(G,\lambda)$ of $G$ are integers if and only if $G$ is the dual of a chordal plane graph. In this article, we further prove that if $F(G,\lambda)$ has real roots only, then $G$ is the dual of a chordal plane graph and each root of $F(G,\lambda)$ is an integer in the set $\{1,2,3\}$.

We construct 2d $\mathcal{N}=(2,2)$ gauged linear sigma models with $\mathrm{U}(1|1)^N$ supergauge group possibly with superpotential. Despite being nonunitary, one can still study their space of supersymmetric states and explore their applications to mathematics. In particular, we find a relation between a nonlinear sigma model on a Calabi-Yau complete intersection of hypersurfaces in a super-Grassmannian and a supergauged Landau-Ginzburg orbifold, which can reduce to a regular Calabi-Yau/Landau-Ginzburg correspondence for complete intersections. This defines a super-Grassmannian/supergroup generalization of the correspondence proved by Clader [1] and Zhao [2]. Similarly, we find a relation between a nonlinear sigma model on a Calabi-Yau hypersurface in a product of super-Grassmannians and a hybrid NLSM/supergauged Landau-Ginzburg orbifold, which can reduce to a regular hybrid Calabi-Yau/Landau-Ginzburg correspondence for hypersurfaces in product space. This defines a super-Grassmannian/supergroup generalization of the correspondence proved by Fan-Jarvis-Ruan [3]. We also find that Calabi-Yau supervector bundles over a super-Grassmannian can undergo a physically related mild topology change which is reducible to a regular Atiyah-type flop transition. This defines a super-Grassmannian generalization of a birational equivalence of Calabi-Yau vector bundles in mathematics. Similarly, we find that a Calabi-Yau complete intersection of quadrics in a super-Grassmannian can also undergo a physically related topology change which is reducible to a regular conifold transition. This defines a super-Grassmannian generalization of a homological projective duality for Calabi-Yau quadrics by Kuznetsov-Perry [4] in mathematics.

For any graph $G$, let $P(G,k)$ and $P_{\ell}(G,k)$ denote the chromatic polynomial and the list-color function of $G$, respectively. It remains an open problem whether, for every graph $G$ and integer $k$, the equality $P(G,k)=P_{\ell}(G,k)>0$ implies that $P(G,k+1)=P_{\ell}(G,k+1)$ also holds. In this paper, we answer this question in the negative. For every integer $k\ge 3$, we construct an infinite family of graphs $G$ such that $P(G,k)=P_{\ell}(G,k)>0$ while $P(G,k+1)>P_{\ell}(G,k+1)$. Moreover, using this infinite family of graphs as attachment gadgets, we further show that any graph $H$ with $P(H,k)=P_{\ell}(H,k)>0$ can be developed into an infinite family of graphs $H'$ with $P(H',k)=P_{\ell}(H',k)>0$ and $P(H',k+1)>P_{\ell}(H',k+1)$.

A minibatch can influence training beyond the update at which it is observed because AdamW stores past gradient information in its optimizer states. We study this delayed effect through paired trajectories that differ only in one gradient update and share the same subsequent training sequence. We formulate AdamW as a finite-horizon input–state–output (ISO) system whose state contains the model parameters and first- and second-moment estimates. Linearizing the joint dynamics yields a signed response operator that maps a localized gradient perturbation to its future loss effects, revealing how optimizer memory shapes their magnitude, timing, and sign. We further derive an exact multistep error decomposition and establish first-order finite-horizon accuracy under local smoothness and controlled activation switching. Experiments validate the response mechanism and optimizer-state effects, while repeated-future analyses reveal substantial prospective structure in delayed influence that can be partially recovered from ISO approximations. Code is available at https://github.com/Kanyooo/Loss_ISO.

Exceptional points (EPs) can produce striking responses, but locating one does not reveal how much of its dynamics an excitation accesses, which responses a detector distinguishes, or whether a missing signal is absent or undetected. We introduce Krylov tomography, a preparation- and measurement-aware framework that uses time-resolved data and models with stated uncertainty limits to answer these questions, while also bounding offset-induced departures from exact-EP behavior. As an explicit illustration of the general framework, we present finite-precision simulations of a red-sideband optomechanical model, whose second-moment coherence sector contains a third-order EP, certifying preparation-sensitive access to two and three directions. Krylov tomography thus bridges EP structure and finite-precision measurements, providing a general framework for non-Hermitian systems.

This paper proposes a controller design method for shaping the controllability Gramian into a desired form to design the effect from exogenous inputs to the system state. Using the Bures–Wasserstein distance, we formulate the shaping problem as the minimization of the distance between the system Gramian and a desired Gramian, and the objective function is shown to be strictly convex on the set of symmetric positive definite matrices. In addition, by deriving a semidefinite programming formulation via a linear matrix inequality (LMI), computational efficiency is improved and additional LMI constraints can be incorporated. When the exogenous input is modeled as Gaussian white noise, the proposed framework is closely related to $H_2$ control, which can be interpreted as a special case of optimal transport. Numerical examples demonstrate anisotropic controllability design for a guidance robot and verify the ability to impose additional directional constraints through LMIs. The numerical examples also confirm that the proposed method approaches $H_2$ control as the desired Gramian tends to zero.

We study small-noise large deviations for the stochastic heat equation (SHE) on the torus with unbounded, multiplicative space-time white noise. We establish uniform large deviation principles (ULDPs) over unbounded sets of initial data, alongside novel well-posedness and regularity results. The ULDPs are obtained for both continuous-in-space and $L^p$ initial data. For these two classes, the ULDP holds uniformly over $L^q$-bounded subsets, with $q<\infty$ and $q<p$ allowed respectively, so these sets are highly unbounded in the state space. The admissible range of $q$ is dictated by the growth of the noise coefficient in the SHE. Crucially, our methods allow us to reach down to $q=1$. This yields ULDPs over $L^1$-bounded subsets, enabling the study of exit times in physical systems with mass conservation.

Meissner polyhedra are constant-width bodies obtained from extremal finite sets of unit diameter. Such a generating set may contain dangling points, namely points having exactly two diametric neighbors. This article studies whether these points can play an essential role in surface-area minimization. Given an extremal set we show that the smallest surface area among the Meissner polyhedra based on it cannot increase by deleting a dangling point. Adding a dangling point cannot decrease the smallest achievable surface area. This reduces the search for area-minimizing Meissner polyhedra to generating sets without dangling points.

We study Brownian motion in R^d conditioned so that the time average of a continuous confining potential remains below a fixed level. On every fixed initial time interval, we prove that the conditioned process converges in total variation to the ground-state diffusion associated with a suitable Schrödinger operator. We also obtain sharp asymptotics for the probability of the conditioning event, including bounded perturbations of the constraint. The proof is based on a local limit theorem for the corresponding Feynman-Kac measures. Our results extend the previously known one-dimensional quadratic case to arbitrary finite dimension and a broad class of confining potentials, therefore resolving a conjecture of Aurzada, Lifshits and Schickentanz. The presented approach also works when Brownian motion is replaced by suitable reversible Markov processes, including multidimensional Ornstein-Uhlenbeck processes, CIR processes and continuous-time Markov chains.

Silverman conjectured that the geometric divisibility sequence attached to a Zariski-dense point on an irreducible commutative algebraic group of dimension at least two, with no unipotent part, returns to its initial value infinitely often. We construct, for the first time, unconditional examples of this phenomenon on geometrically nonsplit semiabelian varieties over number fields using torsion translates of normalized Ribet points. More precisely, let $A/K$ be a positive-dimensional abelian variety over number field $K$, let \[ 1\longrightarrow\mathbf G_m\xrightarrow{\iota}G_q\xrightarrow{\pi}A \longrightarrow0 \] be the extension represented by $q\in A^\vee(K)$, and let $R_\beta(q)\in G_q(K)$ be the normalized Ribet point associated with a homomorphism $\beta:A^\vee\to A$. We set $\delta:=\beta-\widehat\beta$ and assume that $\delta$ is an isogeny and that $\mathbf Z(\delta q)$ is Zariski dense in $A$. For a torsion point $t\in\mathbf G_m(K)_{\mathrm{tors}}$, identify $t$ with $\iota(t)$ and set $P=R_\beta(q)+t$. Then $P$ has Zariski-dense cyclic orbit in the geometrically nonsplit extension $G_q$. There exists an explicit integer $N_{\delta,t}$ such that if $N_{\delta,t}>1$, then $N_{\delta,t}\mid \mathrm{ord}(\overline P_v)$ at all but finitely many places $v$. Consequently, there is a squarefree integer $Q_P>1$ such that \[ (n,Q_P)=1 \quad\Longrightarrow\quad \mathfrak d_{\mathcal N}(nP)=\mathfrak d_{\mathcal N}(P), \] where $\mathfrak d_{\mathcal N}$ denotes the full denominator ideal on the Néron lft-model $\mathcal N$. In particular, we construct explicitly a geometrically nonsplit semiabelian surface $G/\mathbf Q$ and a semiabelian threefold over $\mathbf Q$ satisfying the Silverman conjecture. It follows that we can construct instances satisfying the Silverman conjecture for every dimension at least two.

Transmissive reconfigurable intelligent surfaces (T-RISs) integrated into the transmitter provide a viable realization of tri-hybrid multiple-input multiple-output (MIMO), where spatial processing is distributed across the digital, analog radio-frequency (RF), and electromagnetic (EM) domains. However, unlike conventional RIS-assisted links, a transmitter-native T-RIS directly participates in radiation formation, making T-RIS front-end modeling and weighted sum-rate (WSR)-oriented joint precoding challenging under practical hardware constraints. This paper develops a unified modeling and precoding framework for transmitter-native T-RIS tri-hybrid multi-user (MU) downlink transmission. Specifically, starting from a continuous-field description, the received field is characterized by the interaction among the feed array, the programmable T-RIS aperture, and the user-side propagation, leading to a cascaded baseband input-output model for MU precoding. Furthermore, the same representation is instantiated in the Fresnel, Fraunhofer, and mixed-field regimes, so that near-field focusing and far-field angular steering can be handled within one front-end model. Additionally, a WSR maximization problem is formulated over the digital precoder, analog network, and T-RIS coefficients under power and quantization constraints. A two-level solver is then developed by coupling outer weighted minimum mean-square error (WMMSE) updates with WMMSE-induced aperture-field shaping and hardware projection. Simulations validate the modeling accuracy and convergence, and show that the proposed full tri-hybrid design improves WSR over baselines while suppressing mixed-field cross-regime leakage.

Following \cite {KSW} we continue the study the cellular complexes formed by polynomials of a given degree having a given sequence of multiplicities of real roots. A computer-assisted calculation disproves the earlier conjecture that homology of one point compactification of the closure of such cell is concentrated in at most one degree. Namely, for $\omega=(3,1,1,3)$ in degree $d=18$, the reduced homology is $\ZZ^2$ in degree $7$. The obstruction is already visible in a signed cell count, whose value is $-2$. We relate that count to the rational signed weight enumerator $F_\omega(t)=\sum_{\eta\preceq\omega}(-1)^{\elln(\eta)}t^{|\eta|}$. For the counterexample, $F_\omega(t)=-t^{14}/(1+t^2)^2$, which gives an exact linear formula for the Euler characteristic and forces the total rational Betti number to be unbounded. We prove a number of results and formulate a complete conjecture describing the above homology.

We establish a quantitative stability estimate for the critical Hardy inequality on bounded domains containing the origin. Our result improves the existing quantitative stability estimate by reducing the exponent in the distance function from $N^{2}$ to $N$, replacing the Lorentz-Zygmund framework with the Luxemburg norm of the critical exponential Orlicz space $\operatorname{Exp}L^{\frac{N}{N-1}}(\Omega)$, and avoiding any cut-off modification of the virtual extremizers. As a consequence, we also obtain an improved quantitative stability estimate for the critical Hardy inequality with the logarithmic weight considered by Cianchi and Ferone, where the exponent is likewise reduced from $N^{2}$ to $N$ and the distance is measured directly from the virtual extremizer without truncation. The proof is completely rearrangement-free and relies on a critical Hardy inequality with a remainder term, scale-invariant Sobolev inequalities, and a refined dyadic summation argument. These ingredients yield stronger quantitative stability estimates for both forms of the critical Hardy inequality.

This study examines Cauchy problems governed by highly singular, time-measurable pseudo-differential operators (singular measurable families of Fourier multipliers). We show that the symbols of these operators can exhibit arbitrary blow-up behavior. In particular, we prove the existence and uniqueness of solutions even when the symbols grow super-exponentially in time and frequency. As a concrete application, we solve evolutionary equations driven by fractional Laplacians of any negative order. Additionally, we establish unique strong solutions under the sole condition that the symbol is locally integrable in frequency, even in the presence of severe blow-up at the initial time.

We conduct a retrospective analysis of COVID-19 transmission dynamics across U.S. states using a modified population-based Susceptible-Exposed-Infectious-Asymptomatic-Recovered-Deceased (SEIARD) compartmental model. The proposed framework introduces time-varying transmission, reporting, and mortality rates to capture temporal variations in public behavior and policy interventions during the pandemic. In particular, the transmission rate is modeled as a function of population mobility (derived from Google Mobility Reports), with a residual time-decay term capturing the net effect of unobserved factors such as behavioral adaptation and control measures, while reporting is linked to nationwide testing strategies. We employ a Bayesian approach to integrate multiple data sources and quantify uncertainties in model parameters. The model explicitly distinguishes between symptomatic and asymptomatic infectious individuals and links the latent epidemic states to observable quantities, including reported cases and deaths, through a dynamic reporting function. This retrospective modeling framework provides insights into state-level epidemic trajectories and supports data-driven decision-making for optimal allocation of healthcare resources and evaluation of public health interventions during future pandemics. We further apply a clustering analysis to the posterior parameter estimates to identify groups of U.S. states exhibiting similar epidemiological characteristics, revealing substantial regional heterogeneity in transmission intensity, reproduction dynamics, and mortality burden.

High-multiplicity flexible job shops involve many copies of a small set of job types that must be scheduled on alternative machines. Fluid relaxations provide scalable workload lower bounds, but their fractional machine allocations do not define feasible schedules for individual jobs. We show that, after a suitable finite scaling, an optimal fluid allocation can be realized exactly by a feasible repeating discrete schedule. The construction scales the fluid allocation to integer operation counts, places the resulting operations in nonoverlapping machine intervals, repeats this arrangement, and links operations across repetitions into individual jobs without moving any interval, thereby enforcing job precedence while preserving machine feasibility. For growing finite instances with the same job-type composition, even with a fixed number of extra jobs, the gap between the optimal makespan and the fluid lower bound remains bounded by a constant; hence the relative gap vanishes as the instance grows. Guided by this repeated structure, we develop type-based cyclic template replay (TCTR), which searches job-type templates using an optimization model whose size does not grow with the number of job copies and replays the selected template on the full instance. On 676 multiplicity-expanded public flexible-job-shop instances, TCTR is feasible in every case and achieves a 1.54% mean gap to the fluid lower bound, compared with 5.32% for a job-indexed adaptive large-neighborhood search and 6.18% for a hybrid genetic algorithm.

The four orders of the quaternion algebra B6 carry four reflective wall data on lattices of signature (3,2); three integrate to Borcherds denominators and one is obstructed. The failed denominator survives as a weakly harmonic Maass form, and we prove its shadow is a Hecke eigenform on the line of the newform 6.4.a.a, with zero twist component. The mechanism is invariance selection: the obstruction functional is invariant under the discriminant isometry group, whose invariants in S_{5/2} are one-dimensional; the same mechanism, verified at quaternion discriminants 10 and 22, places the shadows there on 10.4.a.a and 22.4.a.c. On the weight-1/2 layer we prove a determination theorem: the canonical form exists and is unique precisely when the obstruction space vanishes, and among the 71 discriminants below 230 this happens exactly for D in {6, 10, 22}, the genus-zero compact Shimura curves, whose maximal-order ternaries are reflective with integral Weyl chambers of ranks 3, 4, 4; completeness beyond that range is reduced to an estimate on a quadratic Dedekind-type sum, given a bound on the Gauss-sum term. Parity confines this layer to odd channels; on the obstructed orientation the section layer is obstructed outside an explicit 40-element locus of orientations (a double shadow: CM, 36.2.a.a, at weight 3/2 and newform at weight 5/2), while the deck-symmetric directions instead carry a unique canonical weight-1/2 form. The defect invariant satisfies ||Xi||^2 <v_+,v_+> = 144 exactly and equals L(f,2)/48 pi^2 <f,f> to 31 digits. On the section layer the Petersson geometry is rigid: the Gram matrix of S_{3/2}(rho_4) is a single transcendental multiple of an exact rational form, and that transcendental is identified, to 40 digits, as 3 Gamma(1/3)^3 / (2^{7/3} pi^2): the weight-3/2 shadow norms lie in the Chowla-Selberg ring; the weight-5/2 norm is numerically excluded from it.

We study a class of $p$-Laplace equations $$\Delta_p u-\lambda u^{p-1}+ u^{q-1}=0$$ on a closed $n$-dimensional Riemannian manifold $(M,g)$ with $\operatorname{Ric}\geqslant(n-1)g$. For $1<p<2$, $p<q<p^*$, and $0<\lambda<S_{p,q}^{-1}$, where $$S_{p,q}=\frac{q-p}{2}(\frac p n)^{\frac p 2}\Big(\frac{(p^*-1)^2(2-p)}{(p_*-1)(q-1)(p^*-q)}\Big)^{\frac{2-p}{2}},$$ with $p_*=\frac{(n-1)p}{n-p}$ and $p^*=\frac{np}{n-p}$, we prove that the constant $\lambda^{\frac{1}{q-p}}$ is the unique positive solution of the equation. In contrast, for $p>2$ and $p<q<p^*$, the uniqueness fails for every $\lambda>0$; aside from the constant solution, the equation admits a positive nonconstant solution. This answers Véron's problem raised in \cite{Ver92}.