2026-04-15 | | Total: 224
We study a nonlocal diffusion equation of porous medium type featuring a generalised fractional pressure with spatial anisotropy. We construct a finite element method for the numerical solution of the equation on a bounded open Lipschitz polytopal domain $Ω\subset \mathbb{R}^{d}$, where $d = 2$ or $3$. The pressure in the model is defined as the solution of fractional elliptic problem involving the fractional power of a second order differential operator, in terms of its spectral definition. Under suitable assumptions on the fractional order and the coefficients of the operator, we rigorously prove convergence of the numerical scheme. The analysis is carried out in two stages: first passing to the limit in the spatial discretization, and then in the time step, ultimately showing that a subsequence of the sequence of finite element approximations defined by the proposed numerical method converges to a bounded and nonnegative weak solution of the initial-boundary-value problem under consideration. Finally, we present numerical experiments in two dimensions illustrating the computational aspects of the method and highlighting the interplay between nonlocal effects and spatial anisotropy under different configurations. We also show numerically the failure of the comparison principle and exponential decay of the numerical solution to a steady state.
The Bell colouring graph $\mathcal{B}(G)$ of a graph $G$ is the graph whose vertices are the partitions of the vertex set of $G$ into independent sets, with an edge between two partitions if and only if one can be obtained from the other by changing the part of a single vertex of $G$. Given a natural number $k$, the Bell $k$-colouring graph $\mathcal{B}_k(G)$ and the upper-Bell $k$-colouring graph $\mathcal{B}_{\geq k}(G)$ are the induced subgraphs of $\mathcal{B}(G)$ consisting of all partitions with at most $k$ parts and at least $k$ parts, respectively. We determine precisely when two finite graphs have isomorphic Bell colouring graphs. In particular, we show that every $n$-vertex graph $G$ with no vertices of degree $n-1$ is uniquely determined by its Bell colouring graph $\mathcal{B}(G)$, and by its upper-Bell colouring graph $\mathcal{B}_{\geq k}(G)$ if $k\leq n-2$. We also show that every $n$-vertex graph with maximum degree $Δ(G)< \frac{1}{9}n-\frac{1}{3}$ is uniquely determined by its Bell $k$-colouring graph $\mathcal{B}_k(G)$ if $k>χ(G)$. By taking graph complements, each of these results can be restated in terms of partitions into cliques.
We establish sufficient conditions on discrete subsets of $\mathbb{R}^d$ for them to form a uniqueness or a non-uniqueness pair for the fractional Laplacian. Specifically, assuming that $f=0$ on $Λ$ and that $(-Δ)^sf=0$ on $M$, where $Λ, M \subset \mathbb{R}^d$ are discrete, we find sufficient conditions on these sets that force $f$ to vanish identically, and we provide examples in which non-uniqueness occurs. Some of the ideas used in the proofs also extend to a broader class of multiplier operators.
This work establishes that an optimal transport~(OT) problem regularized by a given $f$-divergence admits the same solution as another OT problem regularized by a different $g$-divergence, under an appropriate transformation of the cost function. This structural equivalence between OT problems regularized by distinct divergences, in the sense of sharing the same unique minimizer, is demonstrated within the framework of Polish spaces with bounded cost functions.
We develop a variationally consistent mesoscopic extension of Cosserat elasticity motivated by the breakdown of compatibility in classical formulations. By admitting compatibility-breaking perturbations, the classical theory ceases to remain closed under admissible variations, necessitating an enlargement of the constitutive framework. This leads naturally to a formulation in which torsion and curvature are treated as independent distributed measures of defects. The theory is constructed using a Palatini-type variational approach, with the coframe and connection as independent fields. The resulting Euler--Lagrange equations yield both the standard balance laws and defect-related excitation fields. Material invariance gives rise to configurational forces and moments, which emerge as Noether currents and are directly linked to defect transport governed by the Bianchi identities. The framework provides a unified description of defect kinematics, configurational mechanics, and microstructural evolution. Illustrative examples and numerical evaluations demonstrate how defect transport generates configurational forces and highlight the underlying Maxwell-type structure of the theory. The proposed formulation offers a consistent geometric foundation for the analysis of structured solids with evolving internal geometry and provides a basis for future developments in defect dynamics and dissipative processes.
Recently, the authors of this paper introduced logarithmic Hochschild (co)homology of logarithmic spaces in a geometric way using formality of derived intersections. In this paper, the authors extend the decomposition theorem for the logarithmic Hochschild (co)homology of firm orbifolds to general logarithmic orbifolds and consider two applications of the decomposition theorem. First, we consider two versions of a symmetric product and compute the logarithmic Hochschild homology of them. Second, we show that logarithmic Hochschild homology is invariant under root stack operations.
Drawing inspiration from [BBF15], we construct a family of unbounded quasi-trees for a connected closed oriented surface \(S_g\) of genus \(g\geq 2\), upon which the group \(\Homeo_0(S_g)\) acts coboundedly by isometries. As an application, we show that some surface homeomorphisms preserving a non-sporadic essential subsurface or an essential subsurface homeomorphic to a once-bordered torus can have positive stable commutator length in \(\Homeo_0(S_g)\). Moreover, we provide a version of projection complex that does not require the finiteness conditions.
For stochastic Hamilton-Jacobi (SHJ) equations, instability points are the space-time locations where two eternal solutions with the same asymptotic velocity differ. Another fundamental structure in such equations is shocks, which are the space-time locations where the velocity field is discontinuous. In this work, we study the KPZ fixed point, the central object of the KPZ universality class, which can be viewed as a prototype--albeit degenerate--of an inviscid SHJ equation in one spatial dimension. We describe the geometric structure of the instability region and give a detailed and precise analysis of its interplay with the shock structures of the two eternal solutions. We show that these shock structures allow one to reconstruct the instability region. Along the way, we obtain a complete classification of all possible configurations of semi-infinite geodesics emanating from arbitrary space-time points, in the directed landscape--the random environment in which the KPZ fixed point evolves.
We consider steady states of the incompressible Euler equation on two-dimensional domains. For non-radial analytic steady states on bounded simply connected domains, it was shown previously that there must be a global functional relationship between the stream function and the vorticity. We show that this does not extend to smooth functions even when further structural assumptions are imposed, such as the Morse condition or Arnold's stability criterion. In fact, we show that generic steady states with multiple critical points can be perturbed to a smooth steady state with no (single-valued) functional relation between the stream function and the vorticity. As a consequence, there are "branches" of the set of smooth steady states that are isolated from analytic steady states. In some cases, these branches can even consist entirely of (linearly) stable steady states.
Recently, the statistical properties of empirical Entropic Optimal Transport (EOT) have attracted great interest, as this quantity has been shown to be useful for complex data analysis, among other reasons due to its computational efficiency. In several applications, it has been observed that the EOT plan provides valuable information beyond just the optimal value. For example, in cell biology, colocalization analysis based on the EOT plan has been introduced as a measure for quantification of spatial proximity of different protein assemblies. Despite recent progress in the analysis of its risk properties, a precise understanding of its statistical fluctuations to make it accessible for inference remains elusive to a large extent. In this paper, we derive asymptotic weak convergence result for a large class of functionals of the EOT plan, in which the colocalization process is included. The proof is based on Hadamard differentiability and the extended delta method. As an application, we obtain uniform confidence bands for colocalization curves and bootstrap consistency. Our theory is supported by simulation studies and is illustrated by real world data analysis from mitochondrial protein colocalization.
Let $V$ be a grading-restricted vertex algebra and let $A^\infty(V)=U^\infty(V)/Q^\infty(V)$ be the associative algebra constructed by Huang, where $U^\infty(V)$ is the space of column-finite infinite matrices with entries in V and $Q^\infty(V)$ is an ideal of a (nonassociative) algebra structure on $U^\infty(V)$ defined by Huang. Huang introduced families of elements in $Q^\infty(V)$ and conjectured that these elements generate $Q^\infty(V)$. We discover and prove that Huang's elements all satisfy what we call ``the diagonal shift property". On the other hand, in the case that $V$ is the rank one Heisenberg vertex operator algebra, we construct infinitely many linearly independent elements in $Q^\infty(V)$ that do not satisfy the diagonal shift property. As a corollary, we disprove Huang's conjecture.
The LogQ algorithm encodes Quadratic Unconstrained Binary Optimization (QUBO) problems, which are often encountered in the industry (portfolio optimization, fleet optimization, charging stations, etc.). It was developed within the framework of quantum computing, designed as a pragmatic approach to quantum combinatorial optimization that drastically reduces the number of required qubits and quantum circuit depth. While LogQ has recently been made compliant with gradient-inspired methods, greatly improving parameter optimization efficiency, it still faced hurdles regarding Pauli decomposition and measurement overhead. We here demonstrate that LogQ can be fully reformulated within a classical framework, which effectively eliminates the need for Pauli decomposition and bypasses the measurement challenges altogether. This finally leads to a classical heuristic based on a non-linear continuous relaxation of variables and is, to the best of our knowledge, novel. The LogQ story illustrates how quantum computing can inspire classical algorithms, leading to so-called "quantum-inspired" methods.
We prove that a Le Page-type inequality is also valid for metrically characterizing those JB$^*$-triples that are commutative. More precisely, we establish that the following statements are equivalent for any JB$^*$-triple $E$: $(a)$ $E$ is commutative. $(b)$ There exists $γ>0$ satisfying $$\big\|\{a,b,\{x,y,z\}\}\big\|\leq γ\ \! \big\|\{x,y,\{a,b,z\}\}\big\|, \hbox{ for all } a,b,x,y,z\in E.$$
We consider nonlinear GMRES (NGMRES) as an acceleration technique for the Navier-Stokes Picard iteration, a direction that has not previously been explored. We identify the optimal norm for the least squares optimization problem arising in the NGMRES algorithm, and establish a convergence analysis for NGMRES with general depth that proves NGMRES scales the Picard Lipschitz constant by the gain of the optimization problem. To our knowledge, this is the first convergence proof for NGMRES that identifies the mechanism responsible for convergence acceleration. Numerical experiments demonstrate that the convergence estimates are remarkably sharp. In addition, NGMRES greatly improves the performance of the Picard iteration, even in cases where the unaccelerated iteration diverges.
In 2005, Parreau proved that if a measure preserving system is not strongly mixing then it contains a non-trivial factor that is disjoint from every strongly mixing system. Taking this construction as the starting point, we develop the complementary notions of $\mathcal U$-generated and $\mathcal U$-mixing systems, for a set $\mathcal U$ of ultrafilters, and use them to recover several classical results in ergodic theory as special cases of a unified framework. We prove that a system is $\mathcal U$-mixing if and only if it is disjoint from all $\mathcal U$-generated systems. In fact, we show that if $\mathcal Y$ is a $\mathcal U$-generated system and $\mathcal Z$ is disjoint from every $\mathcal U$-mixing system, then any joining of $\mathcal Y$ and $\mathcal Z$ remains disjoint from all $\mathcal U$-mixing systems. We also show that every partially rigid system is a finite extension of some $\mathcal{U}$-generated system.
The doubling conjecture predicts that a manifold admits positive scalar curvature with mean convex boundary if and only if its double admits positive scalar curvature. We show that it holds true for manifolds where the inclusion of the boundary satisfies a certain split-condition on fundamental groups. Our proof is based on surgery-techniques for positive scalar and mean curvature. If the boundary is non-connected, we use existence of area-minimizing hypersurfaces and the monotonicity-formula. Furthermore, we investigate if a psc-metric on a closed manifold can be adjusted so that a given embedded hypersurface is minimal, stable minimal or totally geodesic. While not true in general, such an adjustment is possible in many cases.
The attitude dynamics of aircraft and spacecraft exhibit significantly nonlinear behaviour. In spacecraft, torque is generated through reaction wheels and control moment gyros. In aircraft, torque is generated using lift on control surfaces. In both cases, complex geometries, unique configurations, and internal/environmental changes imply that models must be identified, verified, and updated using in-flight experimental data. However, this data is often noisy, sparsely sampled, and partial in that modeled states may not be directly measurable. In this paper, we propose a method for estimating key parameters in realistic Ordinary Differential Equation (ODE) models of both spacecraft and aircraft dynamics. This method avoids the need to directly measure state derivatives by coupling sampled outputs using the Picard mapping -- an integral constraint on the solution of the parameterized ODE. This constraint is then enforced, and optimal parameter estimates are found using a gradient contraction algorithm. This algorithm is applied to well-studied models of spacecraft and aircraft motion. First, the algorithm is used to estimate the inertia tensor in a 4 control-moment gyro (CMG) model of spacecraft motion. Second, we estimate the 28 higher-order control surface coefficients in a model of the F/A-18 aircraft.
Let \(H(n)\) denote the Hilbert number, i.e.\ the maximal number of limit cycles of planar polynomial vector fields of degree \(\le n\). A classical lower-bound mechanism for \(H(n)\) is \emph{replication}: one pulls back a vector field by a polynomial map and lifts each existing limit cycle to several disjoint copies while controlling the resulting degree. In this paper we give a fully self-contained replication theorem based on the separable Chebyshev covering \[ Φ(u,v)=(T_m(u),T_m(v)). \] Using the \(m\) monotone full branches of \(T_m\) on \((-1,1)\), we prove that every degree-\(\le n\) polynomial vector field with \(k\) limit cycles gives rise to a degree-\(\le nm+m-1\) polynomial vector field with at least \(m^2k\) limit cycles. Consequently, \[ H(nm+m-1)\ge m^2H(n)\qquad (m\ge 2). \] We then extend the construction to general separable pullbacks \((u,v)\mapsto (p(u),p(v))\), show that Chebyshev attains the maximal possible branch count among degree-\(m\) separable pullbacks, and prove a quadratic ceiling for replication-only schemes: if one iterates separable pullbacks and no additional limit cycles are created beyond those forced by lifting, then the number of resulting limit cycles is at most quadratic in the final degree. This shows that superquadratic lower bounds, such as the known \(n^2\log n\)-type bounds, necessarily require mechanisms beyond pure separable replication. Finally, combining our replication theorem with the strongest currently published seed bounds, we obtain new explicit lower estimates in several degrees, including \begin{gather*} H(14)\ge 252,\qquad H(29)\ge 1080,\\ H(31)\ge 1380,\qquad H(39)\ge 2012. \end{gather*}
Hurwitz theory provides a large variety of enumerative problems related to algebraic geometry, mathematical physics, and combinatorics. We give a general framework to approach the large genus asymptotics of Hurwitz theory using only elementary methods and apply it to several types of Hurwitz numbers: single, double, or with an arbitrary numbers of fixed ramifications; simple and / or including completed cycles type of ramification and / or finitely many blocks of weakly monotone and / or strictly monotone types of ramifications. These, to the best of our knowledge, cover most of the Hurwitz numbers studied, and include for instance correlators of the HCIZ matrix model, Grothendieck dessins d'enfant, weighted Hurwitz numbers, and Gromov-Witten invariants of the Riemann sphere. We also apply our method to b-content Hurwitz numbers. As a specialisation, we recover some previously known about the large genus asymptotics of Hurwitz theory, namely classical results by Hurwitz and recent results of Do-He-Robertson, C. Yang, and results connected to recent work of X. Li.
We prove a general Fueter Theorem over real alternative *-algebras. We show that a suitable power of the Laplacian maps Dunkl-regular functions to Dunkl monogenic functions with axial symmetries. Using the embedding of hypercomplex function theories in the class of Dunkl monogenic functions, we subsume several Fueter-type results known in the literature and obtain the most general form for the action of the Laplacian on function spaces over hypercomplex subspaces. We show that Fueter Theorems are in a one-to-one correspondence with a class of graphs, the Fueter trees, that describe the interactions between Dunkl-regular function spaces and the relation with the iterated Laplacian. We obtain that the number of distinct Fueter trees on a hypercomplex space of dimension $n+1$ is equal to the number of partitions in odd parts of the integer $n$.
We consider the problem of reconstructing missing data on a smooth manifold from incomplete and nonuniform samples. While classical methods for manifold approximation typically assume quasi-uniform data, their performance deteriorates significantly in the presence of large gaps or holes. We propose a unified framework for manifold data imputation that reduces the problem to function reconstruction on locally defined tangent spaces. The approach combines two complementary strategies. The first is a Fourier-based method that determines missing values by prescribing a decay rate of the discrete Fourier coefficients, thereby enforcing high-order smoothness through a global spectral criterion. The second is a local variational method based on minimizing high-order central differences, leading to sparse least-squares systems with favorable stability and conditioning properties. We establish a discrete inverse estimate linking decay of Fourier coefficients to uniform bounds on high-order divided differences, providing a theoretical foundation for the spectral approach. For the variational method, we analyze existence, uniqueness, and scaling behavior, showing that conditioning depends primarily on the geometry of the missing region. These functional reconstruction techniques are integrated with a moving least-squares projection framework to yield a practical algorithm for manifold completion. Numerical experiments, including reconstruction on surfaces with significant missing regions, demonstrate accurate and stable recovery without requiring a global parameterization. The proposed framework provides a flexible and effective approach to manifold data imputation in challenging settings with incomplete data.
Let $d(\cdot)$ denote the natural density on the positive integers. We characterize all sets $A,B$ with positive density satisfying $d(A+B)=d(A)+d(B)$, under the assumption that the two sets are not both contained in a proper finite union of residue classes. This gives a new inverse theorem for Kneser's sumset inequality in the integers, and provides a partial answer to a long-standing open question of Erdős and Graham.
We develop the interpolatory $\mathcal{H}_2$ optimal model reduction framework for linear control systems posed on infinite dimensional state, input and output spaces. Specifically, we consider linear systems formulated as controlled abstract Cauchy problems on a Banach space and approximate them via Petrov-Galerkin projection onto finite dimensional trial and test subspaces. We show that the resulting reduced order transfer function interpolates the original at prescribed points, and we characterize precisely how the projection subspaces must be constructed to enforce this interpolation. Building on this, we develop a data-driven realization framework -- an infinite dimensional analogue of the Loewner approach -- that recovers the system behavior directly from input-output data without requiring access to the underlying operators. Finally, we derive $\mathcal{H}_2$ optimality conditions for the reduced model and show that the classical interpolatory characterization persists in this infinite dimensional setting: first-order optimality requires Hermite interpolation of the transfer function at the mirror images of the reduced model's poles. Taken together, these results establish that the interpolatory $\mathcal{H}_2$ optimal model reduction theory extends naturally and completely to infinite dimensional linear control systems with infinite dimensional input and output spaces.
We consider the Schrödinger equation with a multipoint potential of Bethe-Peierls-Thomas-Fermi type. For this singular potential, we develop scattering and inverse scattering at high energies. In particular, in this framework, our results include analogs of the "regular" Born-Faddeev formula for the scattering amplitude and analogs of related "regular" inverse scattering reconstructions at high energies. Related results for scattering solutions at high energies are also presented.
Consider time-homogeneous discrete-time Markov chains $X$, $Y$, and $Z$ on countable state spaces, considered as stochastic processes with specified initial distributions. Suppose for maps $f$ and $g$ that $(f(X_t))_{t \ge 0}$ and $(g(Y_t))_{t \ge 0}$ are both equal in law to $Z$. We prove that $X$ and $Y$ can be coupled so that $(X_t, Y_t)_{t \ge 0}$ is a homogeneous Markov chain with $f(X_t) = g(Y_t)$ for all $t \ge 0$. Without the assumption that $Z$ is Markov, no such Markov coupling exists in general, even an inhomogeneous one. Moreover, we give an explicit construction of such a coupling, with the additional property that $X$ and $Y$ are conditionally independent given the entire trajectory $(f(X_t))_{t \ge 0}$. Under the further assumption that $X$ and $Y$ are stationary, we construct a coupling having the above properties that is also stationary. In this case, conditional independence holds for the corresponding two-sided chains indexed by $\mathbb{Z}$ (but not necessarily for the one-sided versions). We prove further properties of our couplings in special cases where $f$ or $g$ satisfies the strong lumping condition (also known as Dynkin's condition) or the exact lumping condition (also known as the Pitman-Rogers condition). When $f$ is a strong lumping and $g$ is an exact lumping, we show that our coupling coincides with an intertwining of Markov chains as constructed by Diaconis and Fill.