Papers

The simulation of nonlinear ordinary differential equations on quantum computers is inherently challenging, as quantum gates are linear operators on qubit states. In this paper, we put forth a Koopman-von Neumann (KvN) operator based algorithm for solving nonlinear ordinary differential equations on a quantum computer which overcomes the innate limitations of quantum operations. In this approach, a Liouville probability density is embedded into a wavefunction, the evolution of which is governed by an operator dual to the Koopman operator. The trajectories corresponding to the nonlinear differential equations are reconstructed from the average value of Koopman observables evaluated through quantum measurements. We specifically evaluate the computational limitations of solving nonlinear equations within this Liouville embedding framework. An Ehrenfest-type estimate is derived that relates the reconstruction error to the covariance of the transported density, highlighting the competition between linear stretching and Hessian-induced folding.This leads to a key stability criterion, which depends on the local Ehrenfest-Reynolds number. We discuss the effects of measurement-induced errors including those due to Hadamard-test sampling, amplitude estimation, and bias due to covariance in probabilistic Grover-type inference on the evaluation of trajectories. The limiting bounds on the accuracy of quantum computations are numerically verified for the Lotka-Volterra system and for the quartic oscillator. We observe a rapid growth in computational errors when the Ehrenfest-Reynolds number approaches the predicted threshold. The bounds resulting from our analysis provide key guidelines for selecting the width of the initial Gaussian probability distribution associated with a system of nonlinear ordinary differential equations.

We develop and analyze HDG methods for two central problems in finite element exterior calculus, the Hodge-Dirac problem and the Hodge-Laplace problem, in arbitrary dimension $n$. Our analysis allows for equal-order polynomial spaces with $\mathcal{O}(1)$ penalty parameters, by contrast with previous work on the Hodge-Laplace problem requiring an underlying conforming complex and $\mathcal{O}(h)$ or $\mathcal{O}(h^{-1})$ penalties. For the Hodge-Dirac problem with degree-$r$ polynomials, our error estimates give optimal order-$(r+1)$ convergence in all form degrees under suitable regularity hypotheses. For the Hodge-Laplace problem, we prove optimal order-$(r+1)$ convergence for the $k$-form solution and order-$(r+\frac{1}{2})$ convergence for the $(k \pm 1)$-form auxiliary variables. We obtain improved order-$(r+1)$ auxiliary-variable estimates when $k=1$ and/or $k=n-1$, which in particular sharpens some recent HDG error estimates for the two-dimensional vector Poisson equation. Our analysis also encompasses cases of lower regularity, including solutions with reentrant-corner singularities on non-convex domains, and requires only mild mesh regularity conditions. The results are illustrated by numerical experiments in dimensions two and three.

Let $R$ be a normalized root system in $\mathbb{R}^d$ with reflection group $G$, let $k$ be a $G$-invariant multiplicity function, and let $\mathcal{F}_k$ be the associated Dunkl transform. We write $k_\alpha=k(\alpha)>0$ and $\mathrm{d}\mu_k(x)=w(x)\,\mathrm{d}x$, where $w(x)=\prod_{\alpha\in R_+}|\langle x,\alpha\rangle|^{2k_\alpha}$. We study observability, Hölder-type interpolation, and spectral inequalities for the Dunkl heat equation on $\mathbb{R}^d$. We establish a Bernstein inequality for ordinary derivatives and a Logvinenko–Sereda theorem with a spectral constant of the form $e^{C(1+N)}$ for functions whose Dunkl transforms are supported in $\overline{B(0,N)}$. We characterize observable sets as the measurable sets that are thick with respect to $\mu_k$, and prove the equivalence of the observability, Hölder-type interpolation, and spectral inequalities.

In this paper, we study the asymptotic behavior of the extremal process of a standard one-dimensional branching Brownian motion (BBM) with drift $-\rho>-\sqrt2$ and absorbing barrier at level $-x$. We prove that the two-dimensional point process, with first component being the extremal process of the BBM and the second component being the running minimum of the BBM with drift, converges weakly to a decorated Poisson point process (DPPP) on $\mathbb{R} \times [0, \infty)$. This framework allows us to supply a common limit coupling of the extremal processes for different barrier levels and to explicitly derive the limit, as $t\to\infty$, of the extremal process of the killed BBM killed at level $-x$, demonstrating that the double limit, when $t\to\infty$ first and then $x\to\infty$, of the extremal process of the BBM killed at level $-x$ coincides with the limit of the extremal process of (un-killed) BBM.

This paper is the final installment in our series on the geometric theory of complex Monge-Ampere equations. We study singular Kahler metrics on compact normal Kahler spaces with klt singularities whose volume densities lie in $L^p$ for some $p>1$. We show that a lower bound for the Ricci current in the sense of pluripotential theory is equivalent to a synthetic Ricci lower bound in the sense of RCD theory. As a consequence, every Kahler class on a compact normal Kahler space with klt singularities admits an RCD structure, which makes differential geometric analysis available on singular complex spaces. As an application, we show that the fundamental group of a compact Kahler Calabi-Yau space with klt singularities is almost Abelian. We further establish compactness results for singular Kahler-Einstein spaces.

We study a multidimensional nonlinear filtering model whose coefficients depend on a given observation-adapted predictable process and whose observation drift may grow linearly in both the state and the random input. Due to the unboundedness of the observation drift, a global reference measure is not available. To overcome this hurdle, a localized entropy argument is adapted to prove the stopped likelihood to be a uniformly integrable martingale at each control-energy stopping level. The stopped Zakai equation, and hence, the stopped filtering equation is derived. The global filtering equation is then established by de-localization. The uniqueness of the solution to the stopped Zakai equation is obtained by a duality backward stochastic partial differential equation. This uniqueness then propagates to that of the global filtering equation through the stopped ones. Finally, a stability result is established in $W_1$-distance of measures.

This paper studies the antichain polynomials of $[k]\times P$, where $P$ is a connected minuscule poset. We give a formula for the number of antichains, counted by size, of an arbitrary poset. Using this formula, we present necessary and sufficient conditions for the palindromicity of antichain polynomials for two infinite families of connected minuscule posets. We show that, for every connected minuscule poset $P$, if the antichain polynomial of $[k]\times P$ is palindromic, then it has only real and strictly negative zeros. This result, in particular, gives an affirmative answer to Ding-Dong's conjecture about $\gamma$-positivity of the antichain polynomial of $[k]\times P$. By constructing a bijection between antichains and labeled Dyck paths, we also give a layer-refined enumeration of the antichains in $[2]\times[m]\times[n]$. We establish a relation between such antichains and Clar covers of the hexagonal flakes $O(2,m,n)$, thus prove a conjectured determinantal formula for a family of Zhang-Zhang polynomials. Real-rootedness and stability results are also obtained when the shortest chain has length at most two. Finally, we present infinitely many connected Peck posets whose antichain polynomials are not unimodal, disproving the log-concavity conjecture of Ding and Dong.

This paper is the continuation of our earlier work on Holder estimates for Kahler potentials on compact normal Kahler spaces. We establish a geometric counterpart for analytic Holder regularity of Kolodziej. If a singular metric $\omega_\phi$ has Holder continuous potentials with respect to a smooth background metric, then its distance function is Holder continuous with respect to a smooth background distance. The result holds on normal Kahler spaces and yields compactness of the metric completion as well as its identification with the underlying complex space if $\omega_\phi$ is a Kahler current. Under this positivity assumption, our estimates establish the Holder equivalence between the intrinsic canonical Kahler metrics and the extrinsic smooth metrics on Kahler-RCD spaces, particularly on smoothable Kahler-Einstein spaces.

Let $X$ be a compact Kahler manifold. The analytic stability theorem of Kolodziej for complex Monge-Ampere equation states that for any Kahler metrics $\omega$ and $\omega'$ in the same cohomology class, if their volume measures are bounded in $L^p(X)$ (for some $p>1$) and close in $L^1(X)$, then their Kahler potentials are close in $L^\infty(X)$. In this paper, we establish the geometric stability for complex Monge-Ampère equations that $L^1$-closeness of volume measures implies $L^\infty$-closeness for the induced distance functions by $\omega$ and $\omega'$. Consequently, we prove that any non-smooth Kahler current with volume measure bounded in $L^p$ (for some $p>1$) and Ricci current bounded below induces a unique metric space, which turns out to be a compact RCD space homeomorphic to $X$ itself.

We study a class of complex-valued KdV-type equations with cubic second-order derivative nonlinearities and arbitrary complex coefficients. For sufficiently small initial data in $L^2$, we prove local well-posedness in $H^s(\mathbb R)$ for $s\ge3/4$. The key ingredient is a family of dyadic resolution spaces $Z_k=X_k+Y_k$, designed to overcome the logarithmic divergence arising from high$*$low$*$low interactions in which both derivatives fall on the highest-frequency factor. For the completely integrable third-order Kaup–Newell flow, we establish global-in-time $H^s(\mathbb R)$ bounds for $0\le s<1$ and global well-posedness in $H^s(\mathbb R)$ for every $s\ge3/4$, with no smallness assumption on the initial data.

This paper is devoted to the non-diffusive computation of solutions to $N$-dimensional hyperbolic systems of conservation laws (HCL) and more specifically the accurate approximation of simple waves. We develop the LeafNet algorithm, which relies on: i) a reformulation of HCL as a coupled system involving smooth solutions, and ii) an accurate physics informed algorithms approximating smooth functions. An error estimate analysis and two-dimensional numerical experiments are proposed to illustrate the presented strategy.

Quasi-brittle fracture is prevalent in concrete, rock, ceramics, composites, and masonry, and its simulation faces a trade-off among accuracy, efficiency, and simplicity. The classical Lattice Spring Model (LSM) captures cracking via bond breakage without remeshing, but its elements are empirical and limited to a few tessellable shapes with fixed Poisson's ratios. We propose a tessellated Isotropic Elastic Lattice Spring Model (IELSM) that discretizes continua into polygonal elements with axial springs and a volumetric constraint, achieving isotropic elasticity on arbitrary polygonal tessellations. Macroscopic isotropy reduces to governing equations whose solvability gives a theoretical criterion for element admissibility, proving conventional tessellable elements and extending to arbitrary regular N-gons and concave elements. Exploiting boundary interpolation compatibility with finite elements, IELSM is assembled by direct node sharing, without interface elements or kinematic constraints. Coupled with an isotropic damage model, a pure bending test and four fracture benchmarks show that the coupling preserves displacement accuracy, yields crack paths and load-displacement curves agreeing with experiments and outperforming standard FEM, and is insensitive to mesh refinement. IELSM can also be restricted to damage-prone regions, with the remainder modeled by finite elements via node sharing. For the benchmarks, this reduces nodes by 41.2-80.9% and CPU time by 34.2-85.2% versus full-domain IELSM. The framework advances IELSM element construction from empirical trial and error to theoretical determination and simplifies coupling to node sharing, offering a balanced route for complex quasi-brittle fracture analysis.

Using refined absorption, we prove that for every integer $g\ge 1$ and real $\gamma > 0$, and for sufficiently large $n$, there exists an $n^{-1+\gamma}$-spread distribution on Latin squares of order $n$ and girth at least $g$ that have no proper subsquares. This implies a combinatorial proof of Hilton's conjecture from the 1970s (recently proved algebraically by Allsop and Wanless) that for all sufficiently large $n$, there exists a subsquare-free Latin square of order $n$; indeed, it implies there exist at least $n^{(1-o(1))n^2}$ subsquare-free squares. Simultaneously it also implies the existence of high girth Latin squares (recently proved by Kwan, Sah, Sawhney and Simkin) and even an $n^{-1+\gamma}$-spread distribution on high girth Latin squares.

A novel second-order consistent splitting scheme is proposed for the natural convection equations. The scheme is constructed using Taylor expansions about the time level $t^{n+k}$, where $k$ is a parameter to be determined. It is proved to be stable for all $k>4$, with the present study focusing primarily on the case $k=5$. Compared with the conventional splitting scheme based on Taylor expansions about $t^{n+1}$, the proposed scheme exhibits improved stability. By employing the Sobolev inequality and Gronwall's lemma, we rigorously establish stability of the proposed scheme and derive error estimates in both two and three spatial dimensions. Finally, several numerical experiments are presented to verify the stability and accuracy of the proposed scheme and to demonstrate its effectiveness.

Tambara-Yamagami categories are among the simplest fusion categories beyond the pointed case: one adds a single noninvertible simple object to a finite group of invertible ones. We ask what remains of their (diagrammatic) calculus when the group is replaced by a finite semigroup, and develop the resulting (diagrammatic) calculus.

A finite difference numerical scheme is proposed and analyzed for the Cahn-Hilliard-Navier-Stokes system, combined with a dynamical boundary condition. Such a physical system has potential applications in the moving contact line problem. The boundary profile is governed by a lower-dimensional energy potential, coupled with a non-homogeneous boundary condition for the phase variable. In the numerical design, a convex-splitting approach is applied to the chemical potential in both bulk and surface levels, which leads to a highly coupled nonlinear system. A semi-implicit discretization is taken in the nonlinear fluid convection, as well as the coupled terms between the fluid motion and phase variable evolution. A careful finite difference approximation and convexity analysis reveals that such a numerical system could be represented as a non-symmetric and monotone mapping associated with the fluid convection. In turn, the unique solvability is valid based on the monotonicity argument. The total energy stability analysis is obtained through a careful summation-by-part calculation. In particular, an optimal rate convergence analysis is theoretically established in this work. The discrete mass conservation of the exact solution is required to preserve the mean-zero property of the error function, so that the associated discrete H_h^{-1} norm is well-defined. A combination of the Fourier projection and an auxiliary function is applied to overcome this difficulty. Furthermore, an approach of rough and refined error estimates concludes the desired convergence result. Some numerical results are presented in this article, which demonstrate the robustness of the proposed numerical scheme. In our knowledge, this work provides a theoretical proof of convergence analysis and error estimate for a numerical scheme to the moving contact line problem, for the first time in the literature.

Anderson acceleration (AA) is an effective technique for accelerating fixed point iterations, but it generally lacks global convergence guarantees. We propose a nonmonotone globalized Anderson acceleration framework that retains the local acceleration of AA while ensuring global convergence. For contractive mappings, without requiring prior knowledge of the contraction factor, we prove that the proposed method reduces exactly to pure AA after finitely many iterations and enjoys global $r$-linear convergence for memory size $m\geq1$ and global residual $q$-linear convergence for $m=1$, with convergence factors no greater than the contraction factor of the underlying fixed point mapping. For nonexpansive mappings, we prove that the fixed point residuals converge globally to zero. A class-agnostic parameter selection rule is further developed to ensure these convergence properties in both settings. Numerical experiments demonstrate the robustness and efficiency of the proposed method.

We prove the converse of the Serre-Grothendieck finiteness theorem for the cohomology of coherent sheaves, originally due to Lipman. This result characterizes proper morphisms as those morphisms whose derived pushforward preserves pseudo-coherent complexes. Although the derived pushforward does not reflect pseudo-coherence, we prove that it does after twisting by all perfect complexes. This provides a new characterization of properness.

We determine the sharp weighted power-mean comparisons for the family $\{H_a(r)\}_{r\in (0,1)}$ with the logarithmic interpretation at $a=0$. In the parameter ranges considered by Barnard, Richards and Tiedeman, we give a complete characterization of all orders $\lambda,\mu\in\mathbb{R}$ for which $$A_\lambda(w;1,1-r)\leq H_a(r)\leq A_\mu(w;1,1-r)$$ holds for every $r\in(0,1)$. In particular, our results settle completely their two power-mean conjectures. The two sharp orders are determined by the second-order expansion at $r=0$ and the endpoint matching as $r\to1$. A weighted Wronskian identity and a sign analysis of its residual establish the global inequalities.

Let $F$ be the free group of rank two and, for a prime $p$, let $\mathbf H_p=\mathbf G_p*\mathbf{Ab}_{p-1}$ be the pseudovariety of finite groups having a normal $p$-subgroup with abelian quotient of exponent dividing $p-1$. For $H$ of rank two with $\ab_F(H)=F^{\ab}$ we determine the set $\mathfrak D(H)$ of primes $p$ where $H$ is $\mathbf H_p$-dense in $F$, and we show that this prime set is Frobenian in the sense of Serre: it is governed by a single Laurent polynomial $g$ attached to $H$. We prove that $p\in\mathfrak D(H)$ if and only if $g$ has no zero on the torus $(\F_p^\times)^2$, and hence that $\mathfrak D(H)$ possesses a computable natural density $d(H)$. Exactly one of three cases holds: $\mathfrak D(H)$ is the set of all primes, it is finite, or it is neither and $d(H)$ satisfies $\frac1{|G|}\le d(H)\le1-\frac1{|G|}$, where $G$ is the finite Galois group attached to $g$.

We prove that a hyperkähler manifold of OG6 type is determined, up to isomorphism, by its integral second cohomology endowed with its Hodge structure and Beauville–Bogomolov–Fujiki form. The proof combines the monodromy and chamber theorems of Mongardi–Rapagnetta with the bimeromorphic Torelli theorem. In particular, bimeromorphic OG6 manifolds are biholomorphic.

We say that the Poisson operator $P_r$ is Schwarz-contractive from a space $X$ to a space $Y$ if $\|P_r\|_{X\to Y}\le r$ holds for all $0<r<1$. The classical Schwarz lemma is the case $H^\infty_0\to H^\infty$, the subscript indicating zero mean. We are concerned with the non-holomorphic case $L^\infty_0 \to L^p$. For real functions Schwarz-contractivity holds up to the sharp exponent $p_{\mathbb R}=4.109\ldots$. For complex functions we prove it for $p\le 3$ and conjecture that the critical exponent is $4$.

We consider a stochastic parabolic equation on the $n$-dimensional flat torus $\mathbb{T}^n=(\mathbb{R}/(2\pi\mathbb Z))^n$, $n\in\mathbb N$, whose only control acts in the drift. The multiplicative-noise coefficient is adapted and may depend jointly on the sample point, time, and space. Its spatial profiles are assumed to be real analytic with a uniform positive radius of analyticity, while the corresponding analytic norm is only required to be square integrable in time, uniformly along sample paths. We prove a state-only observability inequality with an $L^1$-in-time mean-square observation norm and deduce null controllability from every measurable spatial set of positive measure by one adapted drift control. The same one-time interpolation estimate also yields approximate controllability to arbitrary square-integrable random terminal targets. The proof is based on a direct analytic energy estimate for the adjoint state, followed by propagation of smallness and a telescoping argument.

We prove that a compact immersed disk of constant mean curvature in $\R^3$ whose boundary is mapped diffeomorphically onto a round circle is a planar disk or a spherical cap.

Let $X$ be a Banach space and let $\mathcal{N}(X)$ denote the family of equivalent norms on $X$. We study Hilbert projective metrics on its projectivization $\mathcal{N}'(X)$ induced by ambient cones of nonnegative functions, with particular attention to the intrinsic Hilbert metric induced by the cone $\mathcal{N}(X)\cup\{0\}$. First, we show that the symmetric logarithmic metric on $\mathcal{N}'(X)$ is the Hilbert projective metric induced by the cone of nonnegative real-valued functions, and we characterize the ambient cones that induce the same metric. We then study the intrinsic order on $\mathcal{N}(X)$. For this purpose, we introduce the triangular defect $\Delta_p(x,y)=p(x)+p(y)-p(x+y)$ and prove that two norms $p,q\in\mathcal{N}(X)$ belong to the same intrinsic part if and only if their triangular defects are uniformly comparable, that is, $a\Delta_p\leq\Delta_q\leq b\Delta_p$ for some $a,b>0$. This yields an explicit formula for the intrinsic Hilbert metric in terms of the pointwise comparison of the norms and of their triangular defects. Finally, the map $\Phi(p)=(p,\Delta_p)$ realizes the intrinsic cone order of $\mathcal{N}(X)$ inside a canonical product cone of nonnegative functions and preserves the Hilbert metric on each intrinsic part.

The first examples of asymmetric L-space knots were found by Baker and Luecke. Among Baker-Luecke knots, the simplest one has braid index 12. Later, it turned out that there are just 9 asymmetric hyperbolic L-space knots in the SnapPy census, and their braid indices take the values 4, 5, 6 and 7. It is known that asymmetric L-space knots have braid index at least 4. Recently, Baker and the second author give infinitely many asymmetric hyperbolic L-space knots with braid index 4. In this paper, we construct an asymmetric hyperbolic L-space knot with arbitrary braid index bigger than 4. In fact, there exist infinitely many such knots for each braid index.

To solve an overdetermined polynomial system numerically, one first makes it square, usually by replacing the given equations with as many random linear combinations as there are unknowns. This is a provably safe step, but it can substantially enlarge the supports. The alternative is to keep that many of the given equations themselves. Selection preserves sparsity but risks geometry: a genuine solution can cease to be an isolated point of the subsystem's zero set. We show that deciding whether a safe choice exists is NP-complete, already for an explicit family of systems of degree three with radical ideal and exactly three simple rational solutions. For strong selection with the nondegenerate rational solutions supplied explicitly, three is the exact threshold when degrees are polynomially bounded: one or two solutions reduce to matroid intersection, three already give NP-completeness. Even without a degree bound, an arbitrarily long list never takes the decision problem beyond NP. On the hard family, five natural notions of a faithful subsystem coincide, and every failing choice fails visibly: its zero set contains an affine subspace through one of the three solutions. A degree-four variant shows that cost information does not help: every candidate that could possibly succeed has mixed volume exactly three, and the problem is NP-complete still. The construction realizes Karp's three-dimensional matching problem as the selection of a square subsystem from the given equations.

For a surface with an odd number of punctures, the moduli space of flat $SU(2)$ connections with traceless holonomy around each puncture is a symplectic manifold. When the moduli space is nonempty, there is a natural homomorphism from the mapping class group of the punctured surface to the symplectic mapping class group of this moduli space. It is shown that this homomorphism is injective if and only if the dimension of the moduli space is greater than $2$. This generalizes work of Seidel and Wehrheim–Woodward. Also given is a complete classification of Lagrangian spheres in the projective plane blown up at $5$ points with its monotone symplectic structure, which is the moduli space for the 5-punctured sphere. Furthermore, it is determined when two such Lagrangian spheres can be displaced by a symplectic isotopy. Results are also obtained regarding Lagrangian spheres in the intersection of two quadrics in $\mathbb{C}\mathbb{P}^5$. The proofs involve instanton Floer theory and results on Heegaard splittings. A main technical result establishes the approximation of any Hamiltonian isotopy of the $SU(2)$ moduli space by holonomy perturbations which are used in instanton homology.

This paper is concerned with the Weyl problem, i.e., the existence of isometric immersions or embeddings of two-spheres into the three-dimensional Euclidean space or general ambient three-manifolds. We establish two results on the degenerate Weyl problem, namely when the Gaussian curvature $K_g$ is only nonnegative rather than strictly positive. First, for a smooth metric $g$ on the two-sphere, if $K_g$ is strictly positive except at finitely many points where the Hessian of $K_g$ is positive definite, then $g$ admits a global $C^{2,1}$-isometric embedding into $\mathbb{R}^3$. It appears to be the first result on the degenerate Weyl problem with purely intrinsic conditions on $g$. Second, for a general simply-connected ambient three-manifold $(\mathcal{M},{\overline{g}})$, if $K_g \geq K_0 \geq {\rm sec}_{\overline{g}}$ for some constant $K_0$, $(K_g-K_0)^{-1/2} \in L^p$ with $p \geq 2$, and a certain uniform pinching condition holds for approximate nondegenerate isometric immersions, then there exists a $W^{3,p}$-isometric immersion. Alongside we also resolve the nondegenerate Weyl problem (i.e., when $K_g>0$) into general simply-connected ambient three-manifolds for $g$, ${\overline{g}} \in C^{2,1}$.

We establish an averaging principle for a class of multiscale mean-field forward-backward stochastic differential equations and identify several novel phenomena that are absent from classical fast-slow systems. In contrast with classical fast-slow systems, the effective dynamics cannot in general be obtained by simply freezing deterministic slow parameters and averaging against the invariant measure of the resulting fast equation. The appropriate averaging object is instead provided by a frozen fast dynamics in a random environment and its associated conditional invariant measures, which retain the coupling between the slow state and its distribution. The forward-backward structure creates a further obstruction: local averaging estimates need not remain stable when propagated over an arbitrary time horizon. We identify a uniform restart stability condition for the averaged system under which this obstruction can be overcome. Using a joint lifted semigroup for the state-law dynamics, together with a two-scale discretization and a Gordin-type decomposition, we prove strong averaging for both the forward and backward components with optimal convergence rate $O(\varepsilon^{1/2})$. As an application, we apply the general theory to a class of mean-field stochastic control problems and develop an efficient algorithm for solving such mean-field control problems.