New articles on Mathematics


[1] 2407.14511

Counting subgroups of the groups ${\Bbb Z}_{n_1} \times \cdots \times {\Bbb Z}_{n_k}$: a survey

We present a survey of exact and asymptotic formulas on the number of cyclic subgroups and total number of subgroups of the groups ${\Bbb Z}_{n_1} \times \cdots \times {\Bbb Z}_{n_k}$, where $k\ge 2$ and $n_1,\ldots,n_k$ are arbitrary positive integers.


[2] 2407.14512

Tetragonal intermediate modular curves

For every group $\{\pm1\}\subseteq \Delta\subseteq (\mathbb{Z}/N\mathbb{Z})^\times$, there exists an intermediate modular curve $X_\Delta(N)$. In this paper we determine all curves $X_\Delta(N)$ whose $\mathbb{Q}$-gonality is equal to $4$, all curves $X_\Delta(N)$ whose $\mathbb{C}$-gonality is equal to $4$, and all curves $X_\Delta(N)$ whose $\mathbb{Q}$-gonality is equal to $5$. We also determine the $\mathbb{Q}$-gonality of all curves $X_\Delta(N)$ for $N\leq 40$ and $\{\pm1\}\subsetneq \Delta \subsetneq (\mathbb{Z}/N\mathbb{Z})^\times$.


[3] 2407.14513

Modular Deutsch Entropic Uncertainty Principle

Khosravi, Drnov\v{s}ek and Moslehian [\textit{Filomat, 2012}] derived Buzano inequality for Hilbert C*-modules. Using this inequality we derive Deutsch entropic uncertainty principle for Hilbert C*-modules over commutative unital C*-algebras.


[4] 2407.14515

Newton-Okounkov bodies and toric degenerations

In recent times, a wide variety of combinatorics has been introduced in order to solve problems from algebraic geometry. Newton-Okounkov bodies and tropical geometry are two such combinatorial theories. As shown by Kaveh and Manon, there is a certain correspondence between these two. Building on this correspondence, and exploiting the link of both theories to toric degenerations, Harada and Escobar obtained their wall-crossing result for prime cones. This result states that moving between two adjacent prime maximal cones in a tropical variety corresponds to a mutation between the associated Newton-Okounkov bodies of these cones. In this thesis, we provide a method for applying the wall-crossing result to non-prime cones. Our approach uses a procedure developed by Bossinger, Lamboglia, Mincheva and Mohammadi in order to compute an embedding which changes a tropical variety in such a way that a non-prime cone becomes prime. Assuming that adjacent cones stay adjacent, the wall-crossing result can then be applied in this new embedding. Building on computations by Clarke, Mohammadi and Zaffalon, we show that for Gr(3,6), this approach works. We compute the new embedding and its tropicalization in this case, and study the relation between cones in the new embedding and the original embedding.


[5] 2407.14517

A group-theoretic condition equivalent to a condition on principal blocks

In this note, we give a group-theoretic condition which is equivalent to the fact that the trivial character is the only complex irreducible character of a finite group G which is contained in the principal p-block for each prime p in a specified set of prime divisors of the order of G. This character-theoretic condition has been studied previously by a number of authors.


[6] 2407.14526

A Random Matrix Model for a Family of Cusp Forms

The Katz-Sarnak philosophy states that statistics of zeros of $L$-function families near the central point as the conductors tend to infinity agree with those of eigenvalues of random matrix ensembles as the matrix size tends to infinity. While numerous results support this conjecture, S. J. Miller observed that for finite conductors, very different behavior can occur for zeros near the central point in elliptic curve $L$-function families. This led to the creation of the excised model of Due\~{n}ez, Huynh, Keating, Miller, and Snaith, whose predictions for quadratic twists of a given elliptic curve are well fit by the data. The key ingredients are relating the discretization of central values of the $L$-functions to excising matrices based on the value of the characteristic polynomials at 1 and using lower order terms (in statistics such as the one-level density and pair-correlation) to adjust the matrix size. We extended this model for a family of twists of an $L$-function associated to a given holomorphic cuspidal newform of odd prime level and arbitrary weight. We derive the corresponding "effective" matrix size for a given form by computing the one-level density and pair-correlation statistics for a chosen family of twists, and we show there is no repulsion for forms with weight greater than 2 and principal nebentype. We experimentally verify the accuracy of the model, and as expected, our model recovers the elliptic curve model.


[7] 2407.14545

Gauged linear sigma model in geometric phase. II. the virtual cycle

We provide the detailed construction of the virtual cycles needed for defining the cohomological field theory associated to a gauged linear sigma model in geometric phase.


[8] 2407.14547

Convexity and concavity of a class of functions related to the elliptic functions

We investigate the convexity property on $(0,1)$ of the function $$f_a(x)=\frac{{\cal K}{(\sqrt x)}}{a-(1/2)\log(1-x)}.$$ We show that $f_a$ is strictly convex on $(0,1)$ if and only if $a\geq a_c$ and $1/f_a$ is strictly convex on $(0,1)$ if and only if $a\leq\log 4$, where $a_c$ is some critical value. The second main result of the paper is to study the log-convexity and log-concavity of the function $$h_p(x)=(1-x)^p{\cal K}(\sqrt x).$$ We prove that $h_p$ is strictly log-concave on $(0,1)$ if and only if $p\geq 7/32$ and strictly log-convex if and only if $p\leq 0$. This solves some problems posed by Yang and Tian and complete their result and a result of Alzer and Richards that $f_a$ is strictly concave on $(0,1)$ if and only if $a=4/3$ and $1/f_a$ is strictly concave on $(0,1)$ if and only if $a\geq 8/5$. As applications of the convexity and concavity, we establish among other inequalities, that for $a\geq a_c$ and all $r\in(0,1)$ $$\frac{2\pi\sqrt\pi}{(2a+\log 2)\Gamma(3/4)^2}\leq \frac{{\cal K}(\sqrt r)}{a-\frac12\log (r)}+\frac{{\cal K}(\sqrt{1-r})}{a-\frac12\log (1-r)}<1+\frac\pi{2a},$$ and for $p\geq 3(2+\sqrt 2)/8$ and all $r\in(0,1)$ $$\sqrt{(r-r^2)^p{\cal K}(\sqrt{1-r}){\cal K}(\sqrt r)}< \frac{\pi\sqrt\pi}{2^{p+1}\Gamma(3/4)^2}<\frac{r^p{\cal K}(\sqrt{1-r})+(1-r)^p{\cal K}(\sqrt r)}{2}.$$


[9] 2407.14554

On the Distributions of Product and Quotient of two Independent $\hat{I}$-function variates

The study of probability distributions for random variables and their algebraic combinations has been a central focus driving the advancement of probability and statistics. Since the 1920s, the challenge of calculating the probability distributions of sums, differences, products, and quotients of independent random variables have drawn the attention of numerous statisticians and mathematicians who studied the algebraic properties and relationships of random variables. Statistical distributions are highly helpful in data science and machine learning, as they provide a range of possible values for the variables, aiding in the development of a deeper understanding of the underlying problem. In this paper, we have presented a new probability distribution based on the $\hat{I}$-function. Also, we have discussed the applications of the $\hat{I}$ function, particularly in deriving the distributions of product and the quotient involving two independent $\hat{I}$ function variates. Additionally, it has been shown that both the product and quotient of two independent $\hat{I}$-function variates also follow the $\hat{I}$-function distribution. Furthermore, the new distribution, known as the $\hat{I}$-function distribution, includes several well-known classical distributions such as the gamma, beta, exponential, normal H-function, and G-function distributions, among others, as special cases. Therefore, the $\hat{I}$-function distribution can be considered a characterization or generalization of the above-mentioned distributions.


[10] 2407.14566

Generalization Error Analysis of Deep Backward Dynamic Programming for Solving Nonlinear PDEs

We explore the application of the quasi-Monte Carlo (QMC) method in deep backward dynamic programming (DBDP) (Hure et al. 2020) for numerically solving high-dimensional nonlinear partial differential equations (PDEs). Our study focuses on examining the generalization error as a component of the total error in the DBDP framework, discovering that the rate of convergence for the generalization error is influenced by the choice of sampling methods. Specifically, for a given batch size $m$, the generalization error under QMC methods exhibits a convergence rate of $O(m^{-1+\varepsilon})$, where $\varepsilon$ can be made arbitrarily small. This rate is notably more favorable than that of the traditional Monte Carlo (MC) methods, which is $O(m^{-1/2+\varepsilon})$. Our theoretical analysis shows that the generalization error under QMC methods achieves a higher order of convergence than their MC counterparts. Numerical experiments demonstrate that QMC indeed surpasses MC in delivering solutions that are both more precise and stable.


[11] 2407.14569

On right $π$-inverse ordered semigroups

Here we introduce the notion of (left, right) $\pi$-$t$-simple, right $\pi$-inverse ordered semigroups and discuss characterizations and relationships concerning them. Semilattice decomposition of left $\pi$-$t$-simple ordered semigroups has been given here. Furthermore, we study an interrelation between the generalized Green's relations and the class of semigroups which are semilattices of right $\pi$-$t$-simple ordered semigroups.


[12] 2407.14610

On the Shapes of Rational Lemniscates

A rational lemniscate is a level set of $|r|$ where $r: \hat{\mathbb{C}} \rightarrow \hat{\mathbb{C}}$ is rational. We prove that any planar Euler graph can be approximated, in a strong sense, by a homeomorphic rational lemniscate. This generalizes Hilbert's lemniscate theorem; he proved that any Jordan curve can be approximated (in the same strong sense) by a polynomial lemniscate that is also a Jordan curve. As consequences, we obtain a sharp quantitative version of the classical Runge's theorem on rational approximation, and we give a new result on the approximation of planar continua by Julia sets of rational maps.


[13] 2407.14618

SOREL: A Stochastic Algorithm for Spectral Risks Minimization

The spectral risk has wide applications in machine learning, especially in real-world decision-making, where people are not only concerned with models' average performance. By assigning different weights to the losses of different sample points, rather than the same weights as in the empirical risk, it allows the model's performance to lie between the average performance and the worst-case performance. In this paper, we propose SOREL, the first stochastic gradient-based algorithm with convergence guarantees for the spectral risk minimization. Previous algorithms often consider adding a strongly concave function to smooth the spectral risk, thus lacking convergence guarantees for the original spectral risk. We theoretically prove that our algorithm achieves a near-optimal rate of $\widetilde{O}(1/\sqrt{\epsilon})$ in terms of $\epsilon$. Experiments on real datasets show that our algorithm outperforms existing algorithms in most cases, both in terms of runtime and sample complexity.


[14] 2407.14619

The curvature exponent of sub-Finsler Heisenberg groups

The curvature exponent $N_{\mathrm{curv}}$ of a metric measure space is the smallest number $N$ for which the measure contraction property $\mathsf{MCP}(0,N)$ holds. In this paper, we study the curvature exponent of sub-Finsler Heisenberg groups equipped with the Lebesgue measure. We prove that $N_{\mathrm{curv}} \geq 5$, and the equality holds if and only if the corresponding sub-Finsler Heisenberg group is actually sub-Riemannian. Furthermore, we show that for every $N\geq 5$, there is a sub-Finsler structure on the Heisenberg group such that $N_{\mathrm{curv}}=N$.


[15] 2407.14624

Waring Number of Ramified $p$-adic Rings

Let $p$ be prime and $e,k$ be positive integers, and let $R = {\mathbb Z}_p[\sqrt[e]{p}]$. We calculate the Waring numbers $g_R(k)$ for many values of $p, e$, and $k$, and investigate how the Waring numbers for $p=2$ change as $e$ and $k$ vary.


[16] 2407.14626

Independent GUE minor processes of perfect matchings on rail-yard graphs

We study perfect matchings on the rail-yard graphs in which the right boundary condition is given by the empty partition and the left boundary can be divided into finitely many alternating line segments where all the vertices along each line segment are either removed or remained. When the edge weights satisfy certain conditions, we show that the distributions of the locations of certain types of dimers near the right boundary converge to the spectra of independent GUE minor processes. The proof is based on new quantitative analysis of a formula to compute Schur functions at general points discovered in \cite{ZL18}.


[17] 2407.14637

An objective isogeometric mixed finite element formulation for nonlinear elastodynamic beams with incompatible warping strains

We present a stable mixed isogeometric finite element formulation for geometrically and materially nonlinear beams in transient elastodynamics, where a Cosserat beam formulation with extensible directors is used. The extensible directors yields a linear configuration space incorporating constant in-plane cross-sectional strains. Higher-order (incompatible) strains are introduced to correct stiffness, whose additional degrees-of-freedom are eliminated by an element-wise condensation. Further, the present discretization of the initial director field leads to the objectivity of approximated strain measures, regardless of the degree of basis functions. For physical stress resultants and strains, we employ a global patch-wise approxiation using B-spline basis functions, whose higher-order continuity enables to use much less degrees-of-freedom, compared to element-wise approximation. For time-stepping, we employ an implicit energy-momentum consistent scheme, which exhibits superior numerical stability in comparison to standard trapezoidal and mid-point rules. Several numerical examples are presented to verify the present method.


[18] 2407.14647

A combinatorial approach to categorical Möbius inversion and pseudoinversion

We use Cramer's formula for the inverse of a matrix and a combinatorial expression for the determinant in terms of paths of an associated digraph (which can be traced back to Coates) to give a combinatorial interpretation of M\"obius inversion whenever it exists. Every M\"obius coefficient is a quotient of two sums, each indexed by certain collections of paths in the digraph. Our result contains, as particular cases, previous theorems by Hall (for posets) and Leinster (for skeletal categories whose idempotents are identities). A byproduct is a novel expression for the magnitude of a metric space as sum over self-avoiding paths with finitely many terms. By means of Berg's formula, our main constructions can be extended to Moore-Penrose pseudoinverses, yielding an analogous combinatorial interpretation of M\"obius pseudoinversion and, consequently, of the magnitude of an arbitrary finite category.


[19] 2407.14652

The monomial expansion formula for Hall-Littlewood $P$-polynomials

We give a Hecke algebra derivation of Macdonald's expansion formula for Hall-Littlewood polynomials in terms of semistandard Young tableaux. This is accomplished by first obtaining a Hecke algebra lift of the expansion coefficients and then proving a generalization of Klostermann's recursions.


[20] 2407.14654

Uniform dispersion in growth models on homogeneous trees

We consider the dynamics of a population spatially structured in colonies that are vulnerable to catastrophic events occurring at random times, which randomly reduce their population size and compel survivors to disperse to neighboring areas. The dispersion behavior of survivors is critically significant for the survival of the entire species. In this paper, we consider an uniform dispersion scheme, where all possible survivor groupings are equally probable. The aim of the survivors is to establish new colonies, with individuals who settle in empty sites potentially initiating a new colony by themselves. However, all other individuals succumb to the catastrophe. We consider the number of dispersal options for surviving individuals in the aftermath of a catastrophe to be a fixed value $d$ within the neighborhood. In this context, we conceptualize the evolution of population dynamics occurring over a homogeneous tree. We investigate the conditions necessary for these populations to survive, presenting pertinent bounds for survival probability, the number of colonized vertices, the extent of dispersion within the population, and the mean time to extinction for the entire population.


[21] 2407.14656

On the Pair Correlation of Zeros of $L$-Functions for Non-CM Newforms in Shifted Ranges

We study the pair correlation between zeros of a shifted auxiliary $ L $-function attached to a non-CM newform, the scale of which is a fixed constant. We prove an unconditional asymptotic result for the pair correlation and introduce a simplicity hypothesis for the zeros of this function, which if true means that multiple zeros of the original $ L $-function cannot be separated by the same fixed distance. Our results provide macroscopic information in contrast to the pair correlation of the original $ L $-function which is of microscopic nature.


[22] 2407.14659

Equivariant cohomology and rings of functions

This submission is a PhD dissertation. It constitutes the summary of the author's work concerning the relations between cohomology rings of algebraic varieties and rings of functions on zero schemes and fixed point schemes. It includes the results from the co-authored article arXiv:2212.11836. They are complemented by: an introduction to the theory of group actions on algebraic varieties, with particular focus on vector fields; a historical overview of the field; a few newer results of the author. The fundamental theorem from arXiv:2212.11836 says that if the principal nilpotent has a unique zero, then the zero scheme over the Kostant section is isomorphic to the spectrum of the equivariant cohomology ring, remembering the grading in terms of a $\mathbb{C}$ action. In this thesis, we also tackle the case of a singular variety. As long as it is embedded in a smooth variety with regular action, we are able to study its cohomology as well by means of the zero scheme. In largest generality, this allows us to see geometrically a subring of the cohomology ring. We also show that the cohomology ring of spherical varieties appears as the ring of functions on the zero scheme. Lastly, the K-theory conjecture is studied, with some results attained for GKM spaces.


[23] 2407.14660

More on the sum-freedom of the multiplicative inverse function

In two papers entitled ``Two generalizations of almost perfect nonlinearity" and ``On the vector subspaces of $\mathbb F_{2^n}$ over which the multiplicative inverse function sums to zero", the first author has introduced and studied the notion of sum-freedom of vectorial functions, which expresses that a function sums to nonzero values over all affine subspaces of $\Bbb F_{2^n}$ of a given dimension $k\geq 2$, and he then focused on the $k$th order sum-freedom of the multiplicative inverse function $x\in \Bbb F_{2^n}\mapsto x^{2^n-2}$. Some general results were given for this function (in particular, the case of affine spaces that do not contain 0 was solved positively), and the cases of $k\in \{3,n-3\}$ and of $k$ not co-prime with $n$ were solved as well (negatively); but the cases of those linear subspaces of dimension $k\in [\![ 4;n-4]\!]$, co-prime with $n$, were left open. The present paper is a continuation of the previous work. After studying, from two different angles, the particular case of those linear subspaces that are stable under the Frobenius automorphism, we deduce from the second approach that, for $k$ small enough (approximately, $3\le k\leq n/10$), the multiplicative inverse function is not $k$th order sum-free. Finally, we extend a result previously obtained in the second paper mentioned above, and we deduce in particular that, for any even $n$ and every $2\leq k\leq n-2$, the multiplicative inverse function is not $k$th order sum-free.


[24] 2407.14674

Equivariant Smoothing Processes on Currents and Spaces with Bounded Curvature

We introduce actions of a compact Lie group in two regularization processes: in De Rham's approximation process of currents on a smooth manifold by smooth currents, and in a smoothing operator of Riemannian metrics of metric spaces with bounded curvature.


[25] 2407.14675

Dynamic Programming Techniques for Planar Orbital Transfer of Low Earth Orbit Satellites

In this paper, we present an application of the optimal control theory to orbital transfer of Low Earth Orbit satellites. The optimal control problem is treated with Dynamic Programming techniques which require solving the Hamilton--Jacobi--Bellman equation on a suitable state space, with the reconstruction of the optimal controls made in the form of a static feedback. In order to validate the numerical scheme without the complexity of the full model, this first study sets the problem in planar form, thus working in a four-dimensional state space. We will study various techniques to speed up the computation, and assess the accuracy of the numerical solution. This project is born from the attempt of evaluating and applying direct method of optimal control techniques based on Dynamic Programming as a complementary approach to the well known indirect methods, as Pontryagin or Lawden. In particular, the final aim is to treat the case of low thrust engines from real use cases, in the full 3-D problem.


[26] 2407.14678

Modular Character Sheaves on Reductive Lie Algebras

Mirkovi\'c introduced the notion of character sheaves on a Lie algebra. Due to their simple geometric characterization, character sheaves on Lie algebras can be thought of as a simplified model for Lusztig's theory of character sheaves on algebraic groups. We extend the theory to the case of modular coefficients. Along the way, we will reprove some of Mirkovi\'c's results and provide connections with the modular generalized Springer correspondence of Achar, Juteau, Riche, and Williamson.


[27] 2407.14692

Some three-dimensional non-nilpotent Leibniz algebras: automorphism groups

The article presents the structure of the automorphism groups of two types of non-nilpotent Leibniz algebras with a dimension of 3.


[28] 2407.14699

On the smoothing theory delooping of disc diffeomorphism and embedding spaces

The celebrated Morlet-Burghelea-Lashof-Kirby-Siebenmann smoothing theory theorem states that the group $\mathrm{Diff}_\partial(D^n)$ of diffeomorphisms of a disc $D^n$ relative to the boundary is equivalent to $\Omega^{n+1}\left(\mathrm{PL}_n/\mathrm{O}_n\right)$ for any $n\geq 1$ and to $\Omega^{n+1}\left(\mathrm{TOP}_n/\mathrm{O}_n\right)$ for $n\neq 4$. We revise smoothing theory results to show that the delooping generalizes to different versions of disc smooth embedding spaces relative to the boundary, namely the usual embeddings, those modulo immersions, and framed embeddings. The latter spaces deloop as $\mathrm{Emb}_\partial^{fr}(D^m,D^n)\simeq\Omega^{m+1}\left(\mathrm{O}_n\backslash\!\!\backslash\mathrm{PL}_n/\mathrm{PL}_{n,m}\right)\simeq \Omega^{m+1}\left(\mathrm{O}_n\backslash\!\!\backslash\mathrm{TOP}_n/\mathrm{TOP}_{n,m}\right)$ for any $n\geq m\geq 1$ ($n\neq 4$ for the second equivalence), where the left-hand side in the case $n-m=2$ or $(n,m)=(4,3)$ should be replaced by the union of the path-components of $\mathrm{PL}$-trivial knots (framing being disregarded). Moreover, we show that for $n\neq 4$, the delooping is compatible with the Budney $E_{m+1}$-action. We use this delooping to combine the Hatcher $\mathrm{O}_{m+1}$-action and the Budney $E_{m+1}$-action into a framed little discs operad $E_{m+1}^{\mathrm{O}_{m+1}}$-action on $\mathrm{Emb}_\partial^{fr}(D^m,D^n)$.


[29] 2407.14715

Analytic Structure of Stationary Flows of an Ideal Fluid with a Stagnation Point

The flow of an ideal fluid possesses a remarkable property: despite limited regularity of the velocity field, its particle trajectories are analytic curves. In our previous work, this fact was used to introduce the structure of an analytic Banach manifold in the set of 2D stationary flows having no stagnation points. The main feature of our description was to regard the stationary flow as a collection of its analytic flow lines, parameterized non-analytically by values of the stream function $\psi$. In this work, we adapt this description to the case of 2D stationary flows which have a single elliptic stagnation point. Namely, we consider flows in a domain bounded by the graph of analytic function $\rho = b(\varphi)$, where $(\rho,\varphi)$ are polar coordinates centred at the origin. The position $p$ of the stagnation point is an unknown and must be included in the solution. In polar coordinates $(r,\theta)$ centred at $p$, the flow lines are described by graphs of $r=a(\psi,\theta)$, where $a$ is a `partially-analytic' function (analytic in $\theta$, of finite regularity in $\psi$). The equation of stationary flow $\Delta \psi = F(\psi)$ is transformed to the quasilinear elliptic equation $\Xi(a) = F(\psi)$ for the flow lines. The analysis is complicated by the fact that the ellipticity of $\Xi$ degenerates at the stagnation point. We introduce function spaces for the partially-analytic family of flow lines, modelled on the weighted Kondratev spaces, appropriate for the degenerate setting. The equation of stationary flow is thus regarded as an analytic operator equation in complex Banach spaces, with local solution given by the implicit function theorem. In particular, we show that near the circular flow of constant vorticity, the equation has unique solution $p, a(\psi,\theta)$ depending analytically on parameters $b(\varphi)$ and $F(\psi)$.


[30] 2407.14718

A mimetic discretization of Westervelt's equation

A broad class of nonlinear acoustic wave models possess a Hamiltonian structure in their dissipation-free limit and a gradient flow structure for their dissipative dynamics. This structure may be exploited to design numerical methods which preserve the Hamiltonian structure in the dissipation-free limit, and which achieve the correct dissipation rate in the spatially-discrete dissipative dynamics. Moreover, by using spatial discretizations which preserve the de Rham cohomology, the non-evolving involution constraint for the vorticity may be exactly satisfied for all of time. Numerical examples are given using a mimetic finite difference spatial discretization.


[31] 2407.14724

Kernel-induced distance and its applications to Composition operators on Large Bergman spaces

In this paper, we obtain a complete characterization for the compact difference of two composition operators acting on Bergman spaces with a rapidly decreasing weight $\omega=e^{-\eta}$, $\Delta\eta>0$. In addition, we provide simple inducing maps which support our main result. We also study the topological path connected component of the space of all bounded composition operators on $A^2(\omega)$ endowed with the Hilbert-Schmidt norm topology.


[32] 2407.14729

Hochschild cohomology for free semigroup algebras

This paper focuses on the cohomology of operator algebras associated with the free semigroup generated by the set $\{z_{\alpha}\}_{\alpha\in\Lambda}$, with the left regular free semigroup algebra $\mathfrak{L}_{\Lambda}$ and the non-commutative disc algebra $\mathfrak{A}_{\Lambda}$ serving as two typical examples. We establish that all derivations of these algebras are automatically continuous. By introducing a novel computational approach, we demonstrate that the first Hochschild cohomology group of $\mathfrak{A}_{\Lambda}$ with coefficients in $\mathfrak{L}_{\Lambda}$ is zero. Utilizing the Ces\`aro operators and conditional expectations, we show that the first normal cohomology group of $\mathfrak{L}_{\Lambda}$ is trivial. Finally, we prove that the higher cohomology groups of the non-commutative disc algebras with coefficients in the complex field vanish when $|\Lambda|<\infty$. These methods extend to compute the cohomology groups of a specific class of operator algebras generated by the left regular representations of cancellative semigroups, which notably include Thompson's semigroup.


[33] 2407.14745

Randomized Radial Basis Function Neural Network for Solving Multiscale Elliptic Equations

To overcome these obstacles and improve computational accuracy and efficiency, this paper presents the Randomized Radial Basis Function Neural Network (RRNN), an innovative approach explicitly crafted for solving multiscale elliptic equations. The RRNN method commences by decomposing the computational domain into non-overlapping subdomains. Within each subdomain, the solution to the localized subproblem is approximated by a randomized radial basis function neural network with a Gaussian kernel. This network is distinguished by the random assignment of width and center coefficients for its activation functions, thereby rendering the training process focused solely on determining the weight coefficients of the output layer. For each subproblem, similar to the Petrov-Galerkin finite element method, a linear system will be formulated on the foundation of a weak formulation. Subsequently, a selection of collocation points is stochastically sampled at the boundaries of the subdomain, ensuring satisfying $C^0$ and $C^1$ continuity and boundary conditions to couple these localized solutions. The network is ultimately trained using the least squares method to ascertain the output layer weights. To validate the RRNN method's effectiveness, an extensive array of numerical experiments has been executed and the results demonstrate that the proposed method can improve the accuracy and efficiency well.


[34] 2407.14755

Remoteness in the category of bilocales

In locale theory, a sublocale is said to be remote in case it misses every nowhere dense sublocale. In this paper, we introduce and study a new class of sublocales in the category of bilocales, namely (i,j)-remote sublocales. These are bilocalic counterparts of remote sublocales and are the sublocales missing every (i,j)-nowhere dense sublocale, with (i,j)-nowhere dense sublocales being bilocalic counterparts of (\tau_{i},\tau_{j})-nowhere dense subsets in bitopological spaces. A comprehensive study of (i,j)-nowhere dense sublocales is given and we show that in the class of balanced bilocales, a sublocale is (i,j)-nowhere dense if and only if its bilocale closure is nowhere dense. We also consider weakly (i,j)-remote sublocales which are those sublocales missing every clopen (i,j)-nowhere dense sublocale. Furthermore, we extend (i,j)-remoteness to the categories of bitopological spaces as well as normed lattices. In the class of \sup-T_{D} bitopological spaces, a subset A of a bitopological space (X,\tau_{1},\tau_{2}) is (\tau_{i},\tau_{j})-remote if and only if the induced sublocale \widetilde{A} of \tau_{1}\vee\tau_{2} is (i,j)-remote.


[35] 2407.14756

Exponential bounds for the density of the law of the solution of a SDE with locally Lipschitz coefficients

Under the uniform H\"{o}rmander's hypothesis we study smoothness and exponential bounds of the density of the law of the solution of a stochastic differential equation (SDE) with locally Lipschitz drift that satisfy a monotonicity condition. To avoid non-integrability problems we use results about Malliavin differentiability based on the concepts of Ray Absolute Continuity and Stochastic Gate\^aux differentiability.


[36] 2407.14764

A Yannakakis-type theorem for lifts of affine semigroups

Yannakakis' theorem relating the extension complexity of a polytope to the size of a nonnegative factorization of its slack matrix is a seminal result in the study of lifts of convex sets. Inspired by this result and the importance of lifts in the setting of integer programming, we show that a similar result holds for the discrete analog of convex polyhedral cones-affine semigroups. We define the notions of the integer slack matrix and a lift of an affine semigroup. We show that many of the characterizations of the slack matrix in the convex cone setting have analogous results in the affine semigroup setting. We also show how slack matrices of affine semigroups can be used to obtain new results in the study of nonnegative integer rank of nonnegative integer matrices.


[37] 2407.14778

Minimax estimation of functionals in sparse vector model with correlated observations

We consider the observations of an unknown $s$-sparse vector ${\boldsymbol \theta}$ corrupted by Gaussian noise with zero mean and unknown covariance matrix ${\boldsymbol \Sigma}$. We propose minimax optimal methods of estimating the $\ell_2$ norm of ${\boldsymbol \theta}$ and testing the hypothesis $H_0: {\boldsymbol \theta}=0$ against sparse alternatives when only partial information about ${\boldsymbol \Sigma}$ is available, such as an upper bound on its Frobenius norm and the values of its diagonal entries to within an unknown scaling factor. We show that the minimax rates of the estimation and testing are leveraged not by the dimension of the problem but by the value of the Frobenius norm of ${\boldsymbol \Sigma}$.


[38] 2407.14780

Mating parabolic rational maps with Hecke groups

We prove that any degree $d$ rational map having a parabolic fixed point of multiplier $1$ with a fully invariant and simply connected immediate basin of attraction is mateable with the Hecke group $H_{d+1}$, with the mating realized by an algebraic correspondence. This solves the parabolic version of the Bullett-Freiberger Conjecture from 2003 on mateability between rational maps and Hecke groups. The proof is in two steps. The first is the construction of a pinched polynomial-like map which is a mating between a parabolic rational map and a parabolic circle map associated to the Hecke group. The second is lifting this pinched polynomial-like map to an algebraic correspondence via a suitable branched covering.


[39] 2407.14781

Bernstein-von Mises theorems for time evolution equations

We consider a class of infinite-dimensional dynamical systems driven by non-linear parabolic partial differential equations with initial condition $\theta$ modelled by a Gaussian process `prior' probability measure. Given discrete samples of the state of the system evolving in space-time, one obtains updated `posterior' measures on a function space containing all possible trajectories. We give a general set of conditions under which these non-Gaussian posterior distributions are approximated, in Wasserstein distance for the supremum-norm metric, by the law of a Gaussian random function. We demonstrate the applicability of our results to periodic non-linear reaction diffusion equations \begin{align*} \frac{\partial}{\partial t} u - \Delta u &= f(u) \\ u(0) &= \theta \end{align*} where $f$ is any smooth and compactly supported reaction function. In this case the limiting Gaussian measure can be characterised as the solution of a time-dependent Schr\"odinger equation with `rough' Gaussian initial conditions whose covariance operator we describe.


[40] 2407.14786

A Pseudo-Gradient Approach for Model-free Markov Chain Optimization

We develop a first-order (pseudo-)gradient approach for optimizing functions over the stationary distribution of discrete-time Markov chains (DTMC). We give insights into why solving this optimization problem is challenging and show how transformations can be used to circumvent the hard constraints inherent in the optimization problem. The optimization framework is model-free since no explicit model of the interdependence of the row elements of the Markov chain transition matrix is required. Upon the transformation we build an extension of Simultaneous Perturbation Stochastic Approximation (SPSA) algorithm, called stochastic matrix SPSA (SM-SPSA) to solve the optimization problem. The performance of the SM-SPSA gradient search is compared with a benchmark commercial solver. Numerical examples show that SM-SPSA scales better which makes it the preferred solution method for large problem instances. We also apply the algorithm to the maximization of web-page rankings in web-graphs based on a real-life data set. As we explain in the paper, when applying a first-order gradient search one typically encounters a phenomenon which we call ``infliction points," that is, jumps in the optimization trajectories between periods of almost stationary behavior that slow down the optimization. We propose a heuristic for avoiding such infliction points and present a metastudy on a wide range of networks showing the positive effect of our heuristic on the convergence properties of SM-SPSA gradient search.


[41] 2407.14803

Homologically smooth connected cochain DGAs

Let $\mathscr{A}$ be a connected cochain DG algebra such that $H(\mathscr{A})$ is a Noetherian graded algebra. We give some criteria for $\mathscr{A}$ to be homologically smooth in terms of the singularity category, the cone length of the canonical module $k$ and the global dimension of $\mathscr{A}$. For any cohomologically finite DG $\mathscr{A}$-module $M$, we show that it is compact when $\mathscr{A}$ is homologically smooth. If $\mathscr{A}$ is in addition Gorenstein, we get $$\mathrm{CMreg}M = \mathrm{depth}_{\mathscr{A}}\mathscr{A} + \mathrm{Ext.reg}\, M<\infty,$$ where $\mathrm{CMreg}M$ is the Castelnuovo-Mumford regularity of $M$, $\mathrm{depth}_{\mathscr{A}}\mathscr{A}$ is the depth of $\mathscr{A}$ and $ \mathrm{Ext.reg}\, M$ is the Ext-regularity of $M$.


[42] 2407.14805

Homological properties of homologically smooth connected cochain DGAs

Let $\mathscr{A}$ be a homologically smooth connected cochain DG algebra. We give some criteria for $\mathscr{A}$ to be Gorenstein and Calabi-Yau, respectively. It is proved that $\mathscr{A}$ is Gorenstein if and only if its $\mathrm{Ext}$-algebra $H(R{\rm Hom}_{\mathscr{A}}(\mathbb{k},\mathbb{k}))$ is a graded Frobenius algebra, and $\mathscr{A}$ is furthermore Calabi-Yau if and only if the Frobenius graded $\mathrm{Ext}$-algebra is symmetric. These results generalize respectively the corresponding results in [HW1] and [HM], where the connected cochain DG algebras studied are Koszul.


[43] 2407.14809

Central extensions, derivations, and automorphisms of semi-direct sums of the Witt algebra with its intermediate series modules

Lie algebras formed via semi-direct sums of the Witt algebra $\text{Der}(\mathbb{C}[t,t^{-1}])$ and its modules have become increasingly prominent in both physics and mathematics in recent years. In this paper, we complete the study of (Leibniz) central extensions, derivations and automorphisms of the Lie algebras formed from the semi-direct sum of the Witt algebra and its indecomposable intermediate series modules (that is, graded modules with one-dimensional graded components). Our techniques exploit the internal grading of the Witt algebra, which can be applied to a wider class of graded Lie algebras.


[44] 2407.14813

Ultraspherical Spectral Method for Block Copolymer Systems on Unit Disk

In this paper, we investigate ultraspherical spectral method for the Ohta-Kawasaki (OK) and Nakazawa-Ohta (NO) models in the disk domain, representing diblock and triblock copolymer systems, respectively. We employ ultraspherical spectral discretization for spatial variables in the disk domain and apply the second-order backward differentiation formula (BDF) method for temporal discretization. To our best knowledge, this is the first study to develop a numerical method for diblock and triblock copolymer systems with long-range interactions in disk domains. We show the energy stability of the numerical method in both semi-discrete and fully-discrete discretizations. In our numerical experiments, we verify the second-order temporal convergence rate and the energy stability of the proposed methods. Our numerical results show that the coarsening dynamics in diblock copolymers lead to bubble assemblies both inside and on the boundary of the disk. Additionally, in the triblock copolymer system, we observe several novel pattern formations, including single and double bubble assemblies in the unit disk. These findings are detailed through extensive numerical experiments.


[45] 2407.14815

Unified Far-Field and Near-Field in Holographic MIMO: A Wavenumber-Domain Perspective

This article conceives a unified representation for near-field and far-field holographic multiple-input multiple-output (HMIMO) channels, addressing a practical design dilemma: "Why does the angular-domain representation no longer function effectively?" To answer this question, we pivot from the angular domain to the wavenumber domain and present a succinct overview of its underlying philosophy. In re-examining the Fourier plane-wave series expansion that recasts spherical propagation waves into a series of plane waves represented by Fourier harmonics, we characterize the HMIMO channel employing these Fourier harmonics having different wavenumbers. This approach, referred to as the wavenumebr-domain representation, facilitates a unified view across the far-field and the near-field. Furthermore, the limitations of the DFT basis are demonstrated when identifying the sparsity inherent to the HMIMO channel, motivating the development of a wavenumber-domain basis as an alternative. We then present some preliminary applications of the proposed wavenumber-domain basis in signal processing across both the far-field and near-field, along with several prospects for future HMIMO system designs based on the wavenumber domain.


[46] 2407.14818

The Rescaled Dirac operator fDh and the noncommutative residue for 6-dimensional manifolds

In this paper, we compute the noncommutative residue for the rescaled Dirac operator fDh on 6-dimensional compact manifolds without boundary. And we give the proof of the Kastler-Kalau-Walze type theorem for the rescaled Dirac operator fDh on 6-dimensional compact manifolds with boundary. We also give some important special cases which can be solved by our calculation methods.


[47] 2407.14819

A Convex-Nonconvex Framework for Enhancing Minimization Induced Penalties

This paper presents a novel framework for nonconvex enhancement of minimization induced (MI) penalties while preserving the overall convexity of associated regularization models. MI penalties enable the adaptation to certain signal structures via minimization, but often underestimate significant components owing to convexity. To overcome this shortcoming, we design a generalized Moreau enhanced minimization induced (GME-MI) penalty by subtracting from the MI penalty its generalized Moreau envelope. While the proposed GME-MI penalty is nonconvex in general, we derive an overall convexity condition for the GME-MI regularized least-squares model. Moreover, we present a proximal splitting algorithm with guaranteed convergence to a globally optimal solution of the GME-MI model under the overall convexity condition. Numerical examples illustrate the effectiveness of the proposed framework.


[48] 2407.14832

On a fibre bundle version of the Caporaso-Harris formula

The Caporaso-Harris formula gives a recursive algorithm to enumerate delta nodal degree d curves in P^2. The recursion is obtained in terms of curves of lower degree that are tangent to a given divisor. This paper presents two generalizations of this method. The first result is on enumeration of one cuspidal curves on P^2, and the second result is an extension to the fiber bundle setting. We solve the question of counting the characteristic number planar nodal cubics in P^3 by extending the idea of Caporaso-Harris.


[49] 2407.14835

Sum of the exponential and a polynomial: Singular values and Baker wandering domains

This article studies the singular values of entire functions of the form $E^k (z)+P(z)$ where $E^k$ denotes the $k-$times composition of $e^z$ with itself and $P$ is any non-constant polynomial. It is proved that the full preimage of each neighborhood of $\infty$ is an infinitely connected domain without having any unbounded boundary component. Following the literature, the point at $\infty$ is called a Baker omitted value for the function in such a situation. More importantly, there are infinitely many critical values and no finite asymptotic value and in fact, the set of all critical values is found to be unbounded for these functions. We also investigate the iteration of three examples of entire functions with Baker omitted value and prove that these do not have any Baker wandering domain. There is a conjecture stating that the number of completely invariant domains of a transcendental entire function is at most one. How some of these maps are right candidates to work upon in view of this conjecture is demonstrated.


[50] 2407.14837

Assouad type dimensions of homogeneous Moran sets and Cantor-like sets

In this paper, we give the Assouad dimension formula and the upper bound of the lower dimension for homogeneous Moran sets under the condition $\sup_{k\ge 1}\{n_{k}\}<+\infty$. We also give the Assouad spectrum and the lower spectrum formulas for Cantor-like sets.


[51] 2407.14839

Optimizing over Multiple Distributions under Generalized Quasar-Convexity Condition

We study a typical optimization model where the optimization variable is composed of multiple probability distributions. Though the model appears frequently in practice, such as for policy problems, it lacks specific analysis in the general setting. For this optimization problem, we propose a new structural condition/landscape description named generalized quasar-convexity (GQC) beyond the realms of convexity. In contrast to original quasar-convexity \citep{hinder2020near}, GQC allows an individual quasar-convex parameter $\gamma_i$ for each variable block $i$ and the smaller of $\gamma_i$ implies less block-convexity. To minimize the objective function, we consider a generalized oracle termed as the internal function that includes the standard gradient oracle as a special case. We provide optimistic mirror descent (OMD) for multiple distributions and prove that the algorithm can achieve an adaptive $\tilde{\mathcal{O}}((\sum_{i=1}^d1/\gamma_i)\epsilon^{-1})$ iteration complexity to find an $epsilon$-suboptimal global solution without pre-known the exact values of $\gamma_i$ when the objective admits ``polynomial-like'' structural. Notably, it achieves iteration complexity that does not explicitly depend on the number of distributions and strictly faster $(\sum_{i=1}^d 1/\gamma_i \text{ v.s. } d\max_{i\in[1:d]} 1/\gamma_i)$ than mirror decent methods. We also extend GQC to the minimax optimization problem proposing the generalized quasar-convexity-concavity (GQCC) condition and a decentralized variant of OMD with regularization. Finally, we show the applications of our algorithmic framework on discounted Markov Decision Processes problem and Markov games, which bring new insights on the landscape analysis of reinforcement learning.


[52] 2407.14840

Exponential Fields: Lack of Generic Derivations

We investigate the existence of "generic derivations" in exponential fields. We show that exponential fields without additional compatibility conditions between derivation and exponentiation cannot support a generic derivation.


[53] 2407.14849

Conditional regularity for the compressible Navier-Stokes equations with potential temperature transport

We study conditional regularity for the compressible Navier-Stokes equations with potential temperature transport in a bounded domain $\Omega\subset\mathbb{R}^d$, $d\in\{2,3\}$, with no-slip boundary conditions. We first prove the existence and uniqueness of local-in-time strong solutions. Further, we prove a blow-up criterion for the strong solution in terms of $L^\infty$-norms for the density and the velocity.


[54] 2407.14866

Orientable sequences over non-binary alphabets

We describe new, simple, recursive methods of construction for orientable sequences over an arbitrary finite alphabet, i.e. periodic sequences in which any sub-sequence of n consecutive elements occurs at most once in a period in either direction. In particular we establish how two variants of a generalised Lempel homomorphism can be used to recursively construct such sequences, generalising previous work on the binary case. We also derive an upper bound on the period of an orientable sequence.


[55] 2407.14867

On a level analog of Selberg's result on $S(t)$

Let $S(t,f)=\pi^{-1}\arg L(1/2+it, f)$, where $f$ is a holomorphic Hecke cusp form of weight $2$ and prime level $q$. In this paper, we establish an asymptotic formula for the moments of $S(t,f)$ without assuming the GRH.


[56] 2407.14869

A total Solovay reducibility and totalizing of the notion of speedability

While the set of Martin-L\"of random left-c.e. reals is equal to the maximum degree of Solovay reducibility, Miyabe, Nies and Stephan(DOI:10.4115/jla.2018.10.3) have shown that the left-c.e. Schnorr random reals are not closed upwards under Solovay reducibility. Recall that for two left-c.e. reals~$\alpha$ and~$\beta$, the former is Solovay reducible to the latter in case there is a partially computable function $\varphi$ and constant~$c$ such that for all rational numbers $q < \beta$ we have \[\alpha - \varphi(q) < c(\beta - q).\] By requiring the translation function $\varphi$ to be total, we introduce a total version of Solovay reducibility that implies Schnorr reducibility. Accordingly, by Downey and Griffiths (DOI:10.2178/jsl/1082418542), the set of Schnorr random left-c.e. reals is closed upwards relative to total Solovay reducibility. Furthermore, we observe that the notion of speedability introduced by Merkle and Titov (DOI:10.1007/978-3-030-50026-9_22) can be equivalently characterized via partial computable translation functions in a way that resembles Solovay reducibility. By requiring the translation function to be total, we obtain the concept of total speedability. Like for speedability, this notion does not depend on the choice of the speeding constant.


[57] 2407.14870

On subspaces of Orlicz spaces spanned by independent copies of a mean zero function

We study subspaces of Orlicz spaces $L_M$ spanned by independent copies $f_k$, $k=1,2,\dots$, of a function $f\in L_M$, $\int_0^1 f(t)\,dt=0$. Any such a subspace $H$ is isomorphic to some Orlicz sequence space $\ell_\psi$. In terms of dilations of the function $f$, a description of strongly embedded subspaces of this type is obtained, and conditions, guaranteeing that the unit ball of such a subspace consists of functions with equicontinuous norms in $L_M$, are found. In particular, we prove that there is a wide class of Orlicz spaces $L_M$ (containing $L^p$-spaces, $1\le p< 2$), for which each of the above properties of $H$ holds if and only if the Matuszewska-Orlicz indices of the functions $M$ and $\psi$ satisfy the inequality: $\alpha_\psi^0>\beta_M^\infty$.


[58] 2407.14884

The Lions Derivative in Infinite Dimensions and Higher Order Expansion of Mean-Field SPDEs

In this paper we present a new interpretation of the Lions derivative as the Radon-Nikodym derivative of a vector measure, which provides a canonical extension of the Lions derivative for functions taking values in infinite dimensional Hilbert spaces. This is of particular relevance for the analysis of Hilbert space valued Mean-Field equations. As an illustration we derive a mild Ito-formula for Mean-Field SPDEs, which provides the basis for a higher order Taylor expansion and higher order numerical schemes.


[59] 2407.14889

Inverse spectral problem of Sturm-Liouville equation with many frozen arguments

This research was devoted to investigate the inverse spectral problem of Sturm-Liouville operator with many frozen arguments. Under some assumptions, the authors obtained uniqueness theorems. At the end, a numerical simulation for the inverse problem was presented.


[60] 2407.14890

Rate Splitting Multiple Access: Optimal Beamforming Structure and Efficient Optimization Algorithms

Joint optimization for common rate allocation and beamforming design have been widely studied in rate splitting multiple access (RSMA) empowered multiuser multi-antenna transmission networks. Due to the highly coupled optimization variables and non-convexity of the joint optimization problems, emerging algorithms such as weighted minimum mean square error (WMMSE) and successive convex approximation (SCA) have been applied to RSMA which typically approximate the original problem with a sequence of disciplined convex subproblems and solve each subproblem by an optimization toolbox. While these approaches are capable of finding a viable solution, they are unable to offer a comprehensive understanding of the solution structure and are burdened by high computational complexity. In this work, for the first time, we identify the optimal beamforming structure and common rate allocation for the weighted sum-rate (WSR) maximization problem of RSMA. We then propose a computationally efficient optimization algorithm that jointly optimizes the beamforming and common rate allocation without relying on any toolbox. Specifically, we first approximate the original WSR maximization problem with a sequence of convex subproblems based on fractional programming (FP). Numerical results show that the proposed algorithm achieves the same performance but takes only 0.5\% or less simulation time compared with the state-of-the-art WMMSE, SCA, and FP algorithms. Experimental results show that the proposed two methods have similar mean square error (MSE) performance compared with traditional method using CVX optimization tool, however, computational complexities are greatly reduced. The proposed algorithms pave the way for the practical and efficient optimization algorithm design for RSMA and its applications in 6G.


[61] 2407.14893

Critical dimensions for polyharmonic operators: The Pucci-Serrin conjecture for solutions of bounded energy

We prove a Pucci-Serrin conjecture on critical dimensions under a uniform bound on the energy. The method is based on the analysis of the Green's function of polyharmonic operators with "almost" Hardy potential.


[62] 2407.14901

Uniform K-stability of $G$-varieties of complexity 1

Let ${\rm k}$ be an algebraically closed field of characteristic 0 and $G$ a connect, reductive group over it. Let $X$ be a projective $G$-variety of complexity 1. We classify $G$-equivariant normal test configurations of $X$ with integral central fibre via the combinatorial data. We also give a formula of anti-canonical divisors on $X$. Based on this formula, when $X$ is $\mathbb Q$-Fano, we give an expression of the Futaki invariant, and derive a criterion of uniform K-stability in terms of the combinatorial data.


[63] 2407.14905

On the multicolor Turán conjecture for color-critical graphs

A {\it simple $k$-coloring} of a multigraph $G$ is a decomposition of the edge multiset as a disjoint sum of $k$ simple graphs which are referred as colors. A subgraph $H$ of a multigraph $G$ is called {\it multicolored} if its edges receive distinct colors in a given simple $k$-coloring of $G$. In 2004, Keevash-Saks-Sudakov-Verstra\"{e}te introduced the {\it $k$-color Tur\'{a}n number} $ex_k(n,H)$, which denotes the maximum number of edges in an $n$-vertex multigraph that has a simple $k$-coloring containing no multicolored copies of $H$. They made a conjecture for any $r\geq 3$ and $r$-color-critical graph $H$ that in the range of $k\geq \frac{r-1}{r-2}(e(H)-1)$, if $n$ is sufficiently large, then $ex_k(n, H)$ is achieved by the multigraph consisting of $k$ colors all of which are identical copies of the Tur\'{a}n graph $T_{r-1}(n)$. In this paper, we show that this holds in the range of $k\geq 2\frac{r-1}{r}(e(H)-1)$, significantly improving earlier results. Our proof combines the stability argument of Chakraborti-Kim-Lee-Liu-Seo with a novel graph packing technique for embedding multigraphs.


[64] 2407.14913

Symplectic Leibniz algebras as a non-commutative version of symplectic Lie algebras

We introduce symplectic left Leibniz algebras and symplectic right Leibniz algebras as generalizations of symplectic Lie algebras. These algebras possess a left symmetric product and are Lie-admissible. We describe completely symmetric Leibniz algebras that are symplectic as both left and right Leibniz algebras. Additionally, we show that symplectic left or right Leibniz algebras can be constructed from a symplectic Lie algebra and a vector space through a method that combines the double extension process and the $T^*$-extension. This approach allows us to generate a broad class of examples.


[65] 2407.14915

Black-box Optimization Algorithms for Regularized Least-squares Problems

We consider the problem of optimizing the sum of a smooth, nonconvex function for which derivatives are unavailable, and a convex, nonsmooth function with easy-to-evaluate proximal operator. Of particular focus is the case where the smooth part has a nonlinear least-squares structure. We adapt two existing approaches for derivative-free optimization of nonsmooth compositions of smooth functions to this setting. Our main contribution is adapting our algorithm to handle inexactly computed stationary measures, where the inexactness is adaptively adjusted as required by the algorithm (where previous approaches assumed access to exact stationary measures, which is not realistic in this setting). Numerically, we provide two extensions of the state-of-the-art DFO-LS solver for nonlinear least-squares problems and demonstrate their strong practical performance.


[66] 2407.14917

Distributed Model-Predictive Energy Management Strategy for Shipboard Power Systems Considering Battery Degradation

With the integration of loads such as pulse power loads, a new control challenge is presented in meeting their high ramp rate requirements. Existing onboard generators are ramp rate limited. The inability to meet the load power due to ramp rate limitation may lead to instability. The addition of energy storage elements in addition to the existing generators proves a viable solution in addressing the control challenges presented by high ramp rate loads. A distributed energy management strategy maximizing generator efficiency and minimizing energy storage degradation is developed that facilitates an optimal adaptive power split between generators and energy storage elements. The complex structure of the energy storage degradation model makes it tough for its direct integration into the optimization problem and is not practical for real-time implementation. A degradation heuristic to minimize absolute power extracted from the energy storage elements is proposed as a degradation heuristic measure. The designed strategy is tested through a numerical case study of a consolidated shipboard power system model consisting of a single generator, energy storage element, and load model. The results show the impact of the designed energy management strategy in effectively managing energy storage health.


[67] 2407.14921

AP-MIONet: Asymptotic-preserving multiple-input neural operators for capturing the high-field limits of collisional kinetic equations

In kinetic equations, external fields play a significant role, particularly when their strength is sufficient to balance collision effects, leading to the so-called high-field regime. Two typical examples are the Vlasov-Poisson-Fokker-Planck (VPFP) system in plasma physics and the Boltzmann equation in semiconductor physics. In this paper, we propose a generic asymptotic-preserving multiple-input DeepONet (AP-MIONet) method for solving these two kinetic equations with variable parameters in the high-field regime. Our method aims to tackle two major challenges in this regime: the additional variable parameters introduced by electric fields, and the absence of an explicit local equilibrium, which is a key component of asymptotic-preserving (AP) schemes. We leverage the multiple-input DeepONet (MIONet) architecture to accommodate additional parameters, and formulate the AP loss function by incorporating both the mass conservation law and the original kinetic system. This strategy can avoid reliance on the explicit local equilibrium, preserve the mass and adapt to non-equilibrium states. We demonstrate the effectiveness and efficiency of the proposed method through extensive numerical examples.


[68] 2407.14922

Geometric subfamily of functions convex in some direction and Blaschke products

Consider the family of locally univalent analytic functions $h$ in the unit disk $|z|<1$ with the normalization $h(0)=0$, $h'(0)=1$ and satisfying the condition $${\real} \left( \frac{z h''(z)}{\alpha h'(z)}\right) <\frac{1}{2} ~\mbox{ for $z\in \ID$,} $$ where $0<\alpha\leq1$. The aim of this article is to show that this family has several elegant properties such as involving Blaschke products, Schwarzian derivative and univalent harmonic mappings.


[69] 2407.14929

K-theory of rank one reductive p-adic groups and Bernstein blocks

We prove a colimit formula for the K-theory spectra of reductive p-adic groups of rank one with regular coefficients in terms of the K-theory of certain compact open subgroups. Furthermore, in the complex case, we show, using the construction of types provided by Roche, that this result can be improved to obtain a formula for the K-theory spectrum of every principal series Bernstein block if the group is split.


[70] 2407.14932

An invitation to singular foliations

These lecture notes attempt to invite the reader towards the theory of singular foliations, both smooth and holomorphic. In addition to a systematic review of the foundations, and an attempt to put in order examples and several elementary constructions, we detail several recent tools developed for non-commutative geometry, in particular the holonomy groupoid of Androulidakis and Skandalis and various methods for resolving singularities. We also introduce various homotopic notions, and end with a series of open questions.


[71] 2407.14941

Diffuse Interface Model for Two-Phase Flows on Evolving Surfaces with Different Densities: Local Well-Posedness

A Cahn-Hilliard-Navier-Stokes system for two-phase flow on an evolving surface with non-matched densities is derived using methods from rational thermodynamics. For a Cahn-Hilliard energy with a singular (logarithmic) potential short time well-posedness of strong solutions together with a separation property is shown, under the assumption of a priori prescribed surface evolution. The problem is reformulated with the help of a pullback to the initial surface. Then a suitable linearization and a contraction mapping argument for the pullback system are used. In order to deal with the linearized system, it is necessary to show maximal $L^2$-regularity for the surface Stokes operator in the case of variable viscosity and to obtain maximal $L^p$-regularity for the linearized Cahn-Hilliard system.


[72] 2407.14942

Concentration and limit of large random matrices with given margins

We study large random matrices with i.i.d. entries conditioned to have prescribed row and column sums (margin). This problem has rich connections to relative entropy minimization, Schr\"{o}dinger bridge, the enumeration of contingency tables, and random graphs with given degree sequences. We show that such margin-constrained random matrix is sharply concentrated around a certain deterministic matrix, which we call the \textit{typical table}. Typical tables have dual characterizations: (1) the expectation of the random matrix ensemble with minimum relative entropy from the base model constrained to have the expected target margin, and (2) the expectation of the maximum likelihood model obtained by rank-one exponential tilting of the base model. The structure of the typical table is dictated by two dual variables, which give the maximum likelihood estimates of the tilting parameters. Based on these results, for a sequence of "tame" margins that converges in \( L^{1} \) to a limiting continuum margin as the size of the matrix diverges, we show that the sequence of margin-constrained random matrices converges in cut norm to a limiting kernel, which is the $L^{2}$-limit of the corresponding rescaled typical tables. The rate of convergence is controlled by how fast the margins converge in $L^{1}$. We derive several new results for random contingency tables from our general framework.


[73] 2407.14947

A Distributionally Robust Optimization Framework for Stochastic Assessment of Power System Flexibility in Economic Dispatch

Given the complexity of power systems, particularly the high-dimensional variability of net loads, accurately depicting the entire operational range of net loads poses a challenge. To address this, recent methodologies have sought to gauge the maximum range of net load uncertainty across all buses. In this paper, we consider the stochastic nature of the net load and introduce a distributionally robust optimization framework that assesses system flexibility stochastically, accommodating a minimal extent of system violations. We verify the proposed method by solving the flexibility of the real-time economic dispatch problem on four IEEE standard test systems. Compared to traditional deterministic flexibility evaluations, our approach consistently yields less conservative flexibility outcomes.


[74] 2407.14948

One more proof about the spectrum of Transposition graph

A Transposition graph $T_n$ is defined as a Cayley graph over the symmetric group $Sym_n$ generated by all transpositions. This paper shows how the spectrum of $T_n$ can be obtained using the spectral properties of the Jucys-Murphy elements.


[75] 2407.14952

On the geometric side of the Jacquet-Rallis relative trace formula

We study some aspects of the geometric side of the Jacquet-Rallis relative trace formula. Globally, we compute each geometric term of the Jacquet-Rallis relative trace formula on the general linear group for regular supported test functions. We prove that it can be described by the regular orbital integral. Locally, we show that the regular orbital integral can be compared with the semisimple orbital integral on the unitary group.


[76] 2407.14956

Two-dimensional DtN-FEM scattering analysis of SH guided waves by an interface debonding in a double-layered plate

In this paper, a two-dimensional Dirichlet-to-Neumann (DtN) finite element method (FEM) is developed to analyze the scattering of SH guided waves due to an interface delamination in a bi-material plate. During the finite element analysis, it is necessary to determine the far-field DtN conditions at virtual boundaries where both displacements and tractions are unknown. In this study, firstly, the scattered waves at the virtual boundaries are represented by a superposition of guided waves with unknown scattered coefficients. Secondly, utilizing the mode orthogonality, the unknown tractions at virtual boundaries are expressed in terms of the unknown scattered displacements at virtual boundaries via scattered coefficients. Thirdly, this relationship at virtual boundaries can be finally assembled into the global DtN-FEM matrix to solve the problem. This method is simple and elegant, which has advantages on dimension reduction and needs no absorption medium or perfectly matched layer to suppress the reflected waves compared to traditional FEM. Furthermore, the reflection and transmission coefficients of each guided mode can be directly obtained without post-processing. This proposed DtN-FEM will be compared with boundary element method (BEM), and finally validated for several benchmark problems.


[77] 2407.14964

Projective geometries, $Q$-polynomial structures, and quantum groups

In 2023 we obtained a $Q$-polynomial structure for the projective geometry $L_N(q)$. In the present paper, we display a more general $Q$-polynomial structure for $L_N(q)$. Our new $Q$-polynomial structure is defined using a free parameter $\varphi$ that takes any positive real value. For $\varphi=1$ we recover the original $Q$-polynomial structure. We interpret the new $Q$-polynomial structure using the quantum group $U_{q^{1/2}}(\mathfrak{sl}_2)$ in the equitable presentation. We use the new $Q$-polynomial structure to obtain analogs of the four split decompositions that appear in the theory of $Q$-polynomial distance-regular graphs.


[78] 2407.14965

Commutation relations for functions of canonical conjugate operators

In this work, the commutator of any two reasonable functions of several pairs of canonical conjugate operators is obtained as a sum of terms of partial derivatives of those functions (equations 9, 10 or 11). When applied to quantum mechanics, first term in the sum is formally equivalent to Poisson bracket in classical mechanics, which is a well-known result. The novelty respect other papers is the type of functions of operators considered: equations mentioned are proved valid when each function is an infinite series of positive and negative powers of these operators, as long as every series converges. Therefore, this formalism can be applied to functions such as the Coulomb potential, where the inverse of radial distance cannot be expressed as a Taylor series.


[79] 2407.14978

Approximation of adelic divisors and equidistribution of small points

We study the asymptotic distribution of Galois orbits of algebraic points with small height in a projective variety over a number field. Our main result is a generalization of Yuan's equidistribution theorem that does not require Zhang's lower bound for the essential minimum to be attained. It extends to all projective varieties a theorem of Burgos Gil, Philippon, Rivera-Letelier and the second author for toric varieties, and implies the semiabelian equidistribution theorem of K\"uhne. We also generalize previous work of Chambert-Loir and Thuillier to prove new logarithmic equidistribution results. In the last section, we extend our main theorem to the quasi-projective setting recently introduced by Yuan and Zhang.


[80] 2407.14988

On the boundedness and Schatten class property of noncommutative martingale paraproducts and operator-valued commutators

We study the Schatten class membership of semicommutative martingale paraproducts and use the transference method to describe Schatten class membership of purely noncommutative martingale paraproducts, especially for CAR algebras and $\mathop{\otimes}\limits_{k=1}^{\infty}\mathbb{M}_{d}$ in terms of martingale Besov spaces. Using Hyt\"{o}nen's dyadic martingale technique, we also obtain sufficient conditions on the Schatten class membership and the boundedness of operator-valued commutators involving general singular integral operators. We establish the complex median method, which is applicable to complex-valued functions. We apply it to get the optimal necessary conditions on the Schatten class membership of operator-valued commutators associated with non-degenerate kernels in Hyt\"{o}nen's sense. This resolves the problem on the characterization of Schatten class membership of operator-valued commutators. Our results are new even in the scalar case. Our new approach is built on Hyt\"{o}nen's dyadic martingale technique and the complex median method. Compared with all the previous ones, this new one is more powerful in several aspects: $(a)$ it permits us to deal with more general singular integral operators with little smoothness; $(b)$ it allows us to deal with commutators with complex-valued kernels; $(c)$ it goes much further beyond the scalar case and can be applied to the semicommutative setting. By a weak-factorization type decomposition, we get some necessary but not optimal conditions on the boundedness of operator-valued commutators. In addition, we give a new proof of the boundedness of commutators still involving general singular integral operators concerning $BMO$ spaces in the commutative setting.


[81] 2407.14989

Nonparametric Estimation of Ordinary Differential Equations: Snake and Stubble

We study nonparametric estimation in dynamical systems described by ordinary differential equations (ODEs). Specifically, we focus on estimating the unknown function $f \colon \mathbb{R}^d \to \mathbb{R}^d$ that governs the system dynamics through the ODE $\dot{u}(t) = f(u(t))$, where observations $Y_{j,i} = u_j(t_{j,i}) + \varepsilon_{j,i}$ of solutions $u_j$ of the ODE are made at times $t_{j,i}$ with independent noise $\varepsilon_{j,i}$. We introduce two novel models -- the Stubble model and the Snake model -- to mitigate the issue of observation location dependence on $f$, an inherent difficulty in nonparametric estimation of ODE systems. In the Stubble model, we observe many short solutions with initial conditions that adequately cover the domain of interest. Here, we study an estimator based on multivariate local polynomial regression and univariate polynomial interpolation. In the Snake model we observe few long trajectories that traverse the domain on interest. Here, we study an estimator that combines univariate local polynomial estimation with multivariate polynomial interpolation. For both models, we establish error bounds of order $n^{-\frac{\beta}{2(\beta +1)+d}}$ for $\beta$-smooth functions $f$ in an infinite dimensional function class of H\"older-type and establish minimax optimality for the Stubble model in general and for the Snake model under some conditions via comparison to lower bounds from parallel work.


[82] 2407.14990

On the compactness of the Weyl operator in $\mathcal{S}_ω$

We characterize, using time-frequency analysis, the continuity and compactness of the Weyl operator in global classes of ultradifferentiable functions $\mathcal{S}_\omega$, for weight functions $\omega$ in the sense of Braun, Meise and Taylor. As a consequence, we give results about the compactness of the localization operator in $\mathcal{S}_\omega$, in relation with the spaces of $\omega$-multipliers and $\omega$-convolutors of $\mathcal{S}_\omega$. Moreover, we provide several examples that complement our investigation.


[83] 2407.14992

On the strength of Burer's lifted convex relaxation to quadratic programming with ball constraints

We study quadratic programs with $m$ ball constraints, and the strength of a lifted convex relaxation for it recently proposed by Burer (2024). Burer shows this relaxation is exact when $m=2$. For general $m$, Burer (2024) provides numerical evidence that this lifted relaxation is tighter than the Kronecker product based Reformulation Linearization Technique (RLT) inequalities introduced by Anstreicher (2017), and conjectures that this must be theoretically true as well. In this note, we provide an affirmative answer to this question and formally prove that this lifted relaxation indeed implies the Kronecker inequalities. Our proof is based on a decomposition of non-rank-one extreme rays of the lifted relaxation for each pair of ball constraints. Burer (2024) also numerically observes that for this lifted relaxation, an RLT-based inequality proposed by Zhen et al. (2021) is redundant, and conjectures this to be theoretically true as well. We also provide a formal proof that Zhen et al. (2021) inequalities are redundant for this lifted relaxation. In addition, we establish that Burer's lifted relaxation is a particular case of the moment-sum-of-squares hierarchy.


[84] 2407.14993

Lower Bounds for Nonparametric Estimation of Ordinary Differential Equations

We noisily observe solutions of an ordinary differential equation $\dot u = f(u)$ at given times, where $u$ lives in a $d$-dimensional state space. The model function $f$ is unknown and belongs to a H\"older-type smoothness class with parameter $\beta$. For the nonparametric problem of estimating $f$, we provide lower bounds on the error in two complementary model specifications: the snake model with few, long observed solutions and the stubble model with many short ones. The lower bounds are minimax optimal in some settings. They depend on various parameters, which in the optimal asymptotic regime leads to the same rate for the squared error in both models: it is characterized by the exponent $-2\beta/(2(\beta+1)+d)$ for the total number of observations $n$. To derive these results, we establish a master theorem for lower bounds in general nonparametric regression problems, which makes the proofs more comparable and seems to be a useful tool for future use.


[85] 2407.14998

Nilpotent symplectic alternating algebras III

In this paper we finish our classification of nilpotent symplectic alternating algebras of dimension 10 over any field F.


[86] 2407.14999

From sphere packing to Fourier interpolation

Viazovska's solution of the sphere packing problem in eight dimensions is based on a remarkable construction of certain special functions using modular forms. Great mathematics has consequences far beyond the problems that originally inspired it, and Viazovska's work is no exception. In this article, we'll examine how it has led to new interpolation theorems in Fourier analysis, specifically a theorem of Radchenko and Viazovska.


[87] 2407.15001

Integral and hypergeometric representations for multiple orthogonal polynomials

This paper addresses two primary objectives in the realm of classical multiple orthogonal polynomials with an arbitrary number of weights. Firstly, it establishes new and explicit hypergeometric expressions for type I Hahn multiple orthogonal polynomials. Secondly, applying the residue theorem and the Mellin transform, the paper derives contour integral representations for several families of orthogonal polynomials. Specifically, it presents contour integral formulas for both type I and type II multiple orthogonal polynomials in the Laguerre of the first kind, Jacobi-Pi\~neiro, and Hahn families. The evaluation of these integrals leads to explicit hypergeometric representations.


[88] 2407.15024

Algebraic relations of special $v$-adic arithmetic gamma values and their period interpretations

Let $v$ be a finite place of $\mathbb{F}_q(\theta)$. We show that algebraic relations of $v$-adic arithmetic gamma values over $\mathbb{F}_q(\theta)$ are explained by the standard functional equations together with Thakur's analogue of the Gross-Koblitz formula. A key step is working out a formula expressing $v$-adic crystalline-de Rham periods of Carlitz motives with Complex Multiplication as products of these special gamma values (a $v$-adic Chowla-Selberg formula). The algebraic independence of those periods in question results from determining the dimension of the motivic Galois groups, adapting and further developing existing tools.


[89] 2407.15027

Multifractal level sets and metric mean dimension with potential

Let $(X,f)$ be a dynamical system with the specification property and $\varphi$ be continuous functions. In this paper, we establish some conditional variational principles for the upper and lower Bowen/packing metric mean dimension with potential of multifractal level set $K_\alpha:=\{x\in X:\lim\limits_{n\to\infty}\dfrac{1}{n}\sum\limits_{i=0}^{n-1}\varphi(f^ix)=\alpha\}.$


[90] 2407.15030

Global well-posedness to the Cauchy problem of 2D nonhomogeneous magnetic Bénard system with large initial data and vacuum

This paper establishes the global well-posedness of strong solutions to the nonhomogeneous magnetic B\'enard system with positive density at infinity in the whole space $\mathbb{R}^2$. More precisely, we obtain the global existence and uniqueness of strong solutions for general large initial data. Our method relies on dedicate energy estimates and a logarithmic interpolation inequality.


[91] 2407.15033

Research on Real-time Diagnosis and Early Warning Technology of Automobile Faults Based on Fractional Order Differential Operator

Automobiles have become the main means of transportation for human beings, and their failures in the process of operation are directly related to the life and property safety of drivers. Therefore, automobile fault diagnosis and early warning technologies have become urgent problems in the current academic world. The premise of real-time accurate diagnosis and early warning of automobiles is to obtain high-quality information data in real time, but the automobile operating environment is complex and changeable, resulting in the measured information data under the influence of multiple factors, such as equipment performance and signal interference. There is an unpredictable measurement error, which greatly affects the reliability of fault diagnosis and early warning systems. In this paper, on the basis of studying the structure and operation characteristics of automobiles, we design a method that can be used for real-time diagnosis and early warning of automobile faults; through the study of fractional-order calculus theory, we establish a mathematical model of information data fusion based on fractional-order differential operators. By providing high-quality information data to automobile fault diagnosis and early warning systems, real-time and accurate diagnosis and early warning functions for automobile faults can be realized. The feasibility and effectiveness of the method were verified through an experiment applying the technology in automobile fault diagnosis and early warning. The research results are of great significance for promoting the development of the automobile industry and ensuring the safety of drivers' lives and property.


[92] 2407.15034

Failure of weak-type endpoint restriction estimates for quadratic manifolds

A paper of Beckner, Carbery, Semmes, and Soria proved that the Fourier extension operator associated to the sphere cannot be weak-type bounded at the restriction endpoint $q = 2d/(d-1)$. We generalize their approach to prove that the extension operator associated with any $n$-dimensional quadratic manifold in $\mathbb{R}^d$ cannot be weak-type bounded at $q = 2d/n$. The key step in generalizing the proof of Beckner, Carbery, Semmes, and Soria will be replacing Kakeya sets with what we will call $\mathcal{N}$-Kakeya sets, where $\mathcal{N}$ denotes a closed subset of the Grassmannian $\text{Gr}(d-n,d)$. We define $\mathcal{N}$-Kakeya sets to be subsets of $\mathbb{R}^d$ containing a translate of every $d-n$-plane segment in $\mathcal{N}$. We will prove that if $\mathcal{N}$ is closed and $n$-dimensional, then there exists compact, measure zero $\mathcal{N}$-Kakeya sets, generalizing the same result for standard Kakeya sets.


[93] 2407.15043

XI-DeepONet: An operator learning method for elliptic interface problems

Scientific computing has been an indispensable tool in applied sciences and engineering, where traditional numerical methods are often employed due to their superior accuracy guarantees. However, these methods often encounter challenges when dealing with problems involving complex geometries. Machine learning-based methods, on the other hand, are mesh-free, thus providing a promising alternative. In particular, operator learning methods have been proposed to learn the mapping from the input space to the solution space, enabling rapid inference of solutions to partial differential equations (PDEs) once trained. In this work, we address the parametric elliptic interface problem. Building upon the deep operator network (DeepONet), we propose an extended interface deep operator network (XI-DeepONet). XI-DeepONet exhibits three unique features: (1) The interface geometry is incorporated into the neural network as an additional input, enabling the network to infer solutions for new interface geometries once trained; (2) The level set function associated with the interface geometry is treated as the input, on which the solution mapping is continuous and can be effectively approximated by the deep operator network; (3) The network can be trained without any input-output data pairs, thus completely avoiding the need for meshes of any kind, directly or indirectly. We conduct a comprehensive series of numerical experiments to demonstrate the accuracy and robustness of the proposed method.


[94] 2407.15044

Singular perturbation in heavy ball dynamics

Given a $C^{1,1}_\mathrm{loc}$ lower bounded function $f:\mathbb{R}^n\rightarrow \mathbb{R}$ definable in an o-minimal structure on the real field, we show that the singular perturbation $\epsilon \searrow 0$ in the heavy ball system \begin{equation} \label{eq:P_eps} \tag{$P_\epsilon$} \epsilon\ddot{x}_\epsilon(t) + \gamma\dot{x}_\epsilon(t) + \nabla f(x_\epsilon(t)) = 0, ~~~ \forall t \geqslant 0, ~~~ x_\epsilon(0) = x_0, ~~~ \dot{x}_\epsilon(0) = \dot{x}_0, \end{equation} preserves boundedness of solutions, where $\gamma>0$ is the friction and $(x_0,\dot{x}_0) \in \mathbb{R}^n \times \mathbb{R}^n$ is the initial condition. This complements the work of Attouch, Goudou, and Redont which deals with finite time horizons. In other words, this work studies the asymptotic behavior of a ball rolling on a surface subject to gravitation and friction, without assuming convexity nor coercivity.


[95] 2407.15049

Accelerating Low-Rank Factorization-Based Semidefinite Programming Algorithms on GPU

In this paper, we address a long-standing challenge: how to achieve both efficiency and scalability in solving semidefinite programming problems. We propose breakthrough acceleration techniques for a wide range of low-rank factorization-based first-order methods using GPUs, making the computation much more efficient and scalable. To illustrate the idea and effectiveness of our approach, we use the low-rank factorization-based SDP solver, LoRADS, as an example, which involves both the classic Burer-Monterio method and a novel splitting scheme with a starting logarithmic rank. Our numerical results demonstrate that the accelerated GPU version of LoRADS, cuLoRADS, can solve huge-scale semidefinite programming problems with remarkable efficiency. By effectively leveraging GPU computational power, cuLoRADS exhibits outstanding performance. Specifically, it can solve a set of MaxCut problems with $10^7 \times 10^7$ matrix variables in 10 seconds to 1 minute each on an NVIDIA H100 GPU with 80GB memory, whereas previous solvers demonstrated the capability of handling problems of this scale, required at least dozens of hours per problem on CPUs. Additionally, cuLoRADS shows exceptional scalability by solving 1) a MaxCut problem with a $170 \text{ million} \times 170 \text{ million}$ matrix variable and 2) a Matrix Completion problem with a 20 million $\times$ 20 million matrix variable and approximately 200 million constraints, both in a matter of minutes.


[96] 2407.15052

The ring of differential operators on a quantized flag manifold

We establish some properties of the ring of differential operators on the quantized flag manifold. Especially, we give an explicit description of its localization on an affine open subset in terms of the quantum Weyl algebra ($q$-analogue of boson).


[97] 2407.15053

Stacked Intelligent Metasurfaces for Task-Oriented Semantic Communications

Semantic communication leveraging advanced deep learning (DL) technologies enhances the efficiency, reliability, and security of information transmission. Emerging stacked intelligent metasurface (SIM) having a diffractive neural network (DNN) architecture allows performing complex calculations at the speed of light. In this letter, we introduce an innovative SIM-aided semantic communication system for image recognition tasks. In the considered model, a SIM is positioned in front of the transmitting antenna. In contrast to conventional communication systems transmitting the modulated signals carrying the image information or compressed semantic information, the carrier electromagnetic (EM) wave is directly transmitted from the source in the proposed system. The input layer of the SIM is utilized for source encoding, while the remaining multi-layer architecture constitutes a DNN for semantic encoding. Specifically, the semantic encoder aims to transform the signals passing through the input layer of the SIM into a unique beam towards a receiving antenna corresponding to the image class. Remarkably, both the source and semantic encoding occur naturally as the EM waves propagate through the SIM. At the receiver, the image is recognized by probing the received signal magnitude across the receiving array. To this end, we develop an efficient algorithm to train the transmission coefficients of SIM's meta-atoms to learn the semantic representation of the image. Extensive numerical results verify the effectiveness of utilizing the SIM-based DNN for image recognition task-oriented semantic communications, achieving more than 90% recognition accuracy.


[98] 2407.15057

On a metric symplectization of a contact metric manifold

In this article, we investigate metric structures on the symplectization of a contact metric manifold and prove that there is a unique metric structure, which we call the metric symplectization, for which each slice of the symplectization has a natural induced contact metric structure. We then study the curvature properties of this metric structure and use it to establish equivalent formulations of the $(\kappa, \mu)$-nullity condition in terms of the metric symplectization. We also prove that isomorphisms of the metric symplectizations of $(\kappa, \mu)$-manifolds determine $(\kappa, \mu)$-manifolds up to D-homothetic transformations. These classification results show that the metric symplectization provides a unified framework to classify Sasakian manifolds, K-contact manifolds and $(\kappa, \mu)$-manifolds in terms of their symplectizations.


[99] 2407.15058

Exponential mixing for random nonlinear wave equations: weak dissipation and localized control

We establish a new criterion for exponential mixing of random dynamical systems. Our criterion is applicable to a wide range of systems, including in particular dispersive equations. Its verification is in nature related to several topics, i.e., asymptotic compactness in dynamical systems, global stability of evolution equations, and localized control problems. As an initial application, we exploit the exponential mixing of random nonlinear wave equations with degenerate damping, critical nonlinearity, and physically localized noise. The essential challenge lies in the fact that the weak dissipation and randomness interact in the evolution.


[100] 2407.15075

Mapping classes of Real Rational Surface Automorphisms

Let $\{F_n, n\ge 8\}$ be a family of diffeomorphisms on real rational surfaces that are birationally equivalent to birational maps on $\mathbf{P}^2(\mathbb{R})$. In this article, we investigate the mapping classes of the diffeomorphisms $F_n, n\ge 8$. These diffeomorphisms are reducible with unique invariant irreducible curves, and we determine the mapping classes of their restrictions, $\hat F_n, n \ge 8$, on the cut surfaces, showing that they are pseudo-Anosov and do not arise from Penner's construction. For $n=8$, Lehmer's number is realized as the stretch factor of $\hat F_8$, a pseudo-Anosov map on a once-punctured genus $5$ orientable surface. The diffeomorphism $\hat F_8$ is a new geometric realization of a Lehmer's number.


[101] 2407.15079

Speed of random walk on dynamical percolation in nonamenable transitive graphs

Let $G$ be a nonamenable transitive unimodular graph. In dynamical percolation, every edge in $G$ refreshes its status at rate $\mu>0$, and following the refresh, each edge is open independently with probability $p$. The random walk traverses $G$ only along open edges, moving at rate $1$. In the critical regime $p=p_c$, we prove that the speed of the random walk is at most $O(\sqrt{\mu \log(1/\mu)})$, provided that $\mu \le e^{-1}$. In the supercritical regime $p>p_c$, we prove that the speed on $G$ is of order 1 (uniformly in $\mu)$, while in the subcritical regime $p<p_c$, the speed is of order $\mu\wedge 1$.


[102] 2407.15088

Stationary Solitons in discrete NLS with non-nearest neighbour interactions

The aim of this paper is to provide a construction of stationary discrete solitons in an extended one-dimensional Discrete NLS model with non-nearest neighbour interactions. These models, models of the type with long-range interactions were studied in various other contexts. In particular, it was shown that, if the interaction strength decays sufficiently slowly as a function of distance, it gives rise to bistability of solitons, which may find applications in their controllable switching. Dynamical lattices with long-range interactions also serve as models for energy and charge transport in biological molecules. Using a dynamical systems method we are able to construct, with great accuracy, stationary discrete solitons for our model, for a large region of the parameter space.


[103] 2407.15091

Local models for smooth vector fields of the line

We present the local classification of singularities of smooth vector fields on the line, with respect to the equivalence relation of $C^1$--conjugacy. Along the way, we recall the analogous classification, up to $C^0$ and $C^{\infty}$ conjugacy. We also give the transversal unfoldings of the corresponding normal forms and treat the case where the changes of coordinates are tangent to the identity. Thus, a fairly complete description of the $1$--d case is achieved.


[104] 2407.15092

PFWNN: A deep learning method for solving forward and inverse problems of phase-field models

Phase-field models have been widely used to investigate the phase transformation phenomena. However, it is difficult to solve the problems numerically due to their strong nonlinearities and higher-order terms. This work is devoted to solving forward and inverse problems of the phase-field models by a novel deep learning framework named Phase-Field Weak-form Neural Networks (PFWNN), which is based on the weak forms of the phase-field equations. In this framework, the weak solutions are parameterized as deep neural networks with a periodic layer, while the test function space is constructed by functions compactly supported in small regions. The PFWNN can efficiently solve the phase-field equations characterizing the sharp transitions and identify the important parameters by employing the weak forms. It also allows local training in small regions, which significantly reduce the computational cost. Moreover, it can guarantee the residual descending along the time marching direction, enhancing the convergence of the method. Numerical examples are presented for several benchmark problems. The results validate the efficiency and accuracy of the PFWNN. This work also sheds light on solving the forward and inverse problems of general high-order time-dependent partial differential equations.


[105] 2407.15094

Determining a Time-Varying Potential in Time-Fractional Diffusion from Observation at a Single Point

We discuss the identification of a time-dependent potential in a time-fractional diffusion model from a boundary measurement taken at a single point. Theoretically, we establish a conditional Lipschitz stability for this inverse problem. Numerically, we develop an easily implementable iterative algorithm to recover the unknown coefficient, and also derive rigorous error bounds for the discrete reconstruction. These results are attained by using the (discrete) solution theory of direct problems, and applying error estimates that are optimal with respect to problem data regularity. Numerical simulations are provided to demonstrate the theoretical results.


[106] 2407.15095

On local solubility of Bao--Ratiu equations on surfaces related to the geometry of diffeomorphism group

We are concerned with the existence of asymptotic directions for the group of volume-preserving diffeomorphisms of a closed 2-dimensional surface $(\Sigma,g)$ within the full diffeomorphism group, described by the Bao--Ratiu equations, a system of second-order PDEs introduced in [On a non-linear equation related to the geometry of the diffeomorphism group, Pacific J. Math. 158 (1993); On a non-linear equation related to the geometry of the diffeomorphism group, Pacific J. Math. 158 (1993)]. It is known [The Bao--Ratiu equations on surfaces, Proc. R. Soc. Lond. A 449 (1995)] that asymptotic directions cannot exist globally on any $\Sigma$ with positive curvature. To complement this result, we prove that asymptotic directions always exist locally about a point $x_0 \in \Sigma$ in either of the following cases (where $K$ is the Gaussian curvature on $\Sigma$): (a), $K(x_0)>0$; (b) $K(x_0)<0$; or (c), $K$ changes sign cleanly at $x_0$, i.e., $K(x_0)=0$ and $\nabla K(x_0) \neq 0$. The key ingredient of the proof is the analysis following Han [On the isometric embedding of surfaces with Gauss curvature changing sign cleanly, Comm. Pure Appl. Math. 58 (2005)] of a degenerate Monge--Amp\`{e}re equation -- which is of the elliptic, hyperbolic, and mixed types in cases (a), (b), and (c), respectively -- locally equivalent to the Bao--Ratiu equations.


[107] 2407.15096

Difference of weighted composition operators on weighted Bergman spaces over the unit Ball

In this paper, we characterize the boundedness and compactness of differences of weighted composition operators from weighted Bergman spaces $A^p_\omega$ induced by a doubling weight $\omega$ to Lebesgue spaces $L^q_\mu$ on the unit ball for full $0<p,q<\infty$, which extend many results on the unit disk. As a byproduct, a new characterization of $q$-Carleson the measure for $A^p_\omega$ in terms of the Bergman metric ball is also presented.


[108] 2407.15103

Minimizing Schrödinger eigenvalues for confining potentials

We consider the problem of minimizing the lowest eigenvalue of the Schr\"odinger operator $-\Delta+V$ in $L^2(\mathbb R^d)$ when the integral $\int e^{-tV}\,dx$ is given for some $t>0$. We show that the eigenvalue is minimal for the harmonic oscillator and derive a quantitative version of the corresponding inequality.


[109] 2407.15104

Lifted Hamming and Simplex Codes and Their Designs

A generator matrix of a linear code $\C$ over $\gf(q)$ is also a matrix of the same rank $k$ over any extension field $\gf(q^\ell)$ and generates a linear code of the same length, same dimension and same minimum distance over $\gf(q^\ell)$, denoted by $\C(q|q^\ell)$ and called a lifted code of $\C$. Although $\C$ and their lifted codes $\C(q|q^\ell)$ have the same parameters, they have different weight distributions and different applications. Few results about lifted linear codes are known in the literature. This paper proves some fundamental theory for lifted linear codes, settles the weight distributions of lifted Hamming codes and lifted Simplex codes, and investigates the $2$-designs supported by the lifted Hamming and Simplex codes. Infinite families of $2$-designs are obtained. In addition, an infinite family of two-weight codes and an infinite family of three-weight codes are presented.


[110] 2407.15106

Semipositive line bundles on punctured Riemann surfaces: Bergman kernels and random zeros

We give an extensive study on the Bergman kernel expansions and the random zeros associated with the high tensor powers of a semipositive line bundle on a complete punctured Riemann surface. We prove several results for the zeros of Gaussian holomorphic sections in the semi-classical limit, including the equidistribution, large deviation estimates, central limit theorem, and number variances.


[111] 2407.15107

Momentum Space Feynman Integral for the Bound State Aharonov-Bohm Effect

We construct the Feynman integral for the Schr\"odinger propagator in the polar conjugate momentum space, which describes the bound state Aharonov-Bohm effect, as a well-defined white noise functional.


[112] 2407.15108

Global well-posedness of the energy-critical stochastic nonlinear Schrödinger equation on the three-dimensional torus

We study the Cauchy problem of the defocusing energy-critical stochastic nonlinear Schr\"odinger equation (SNLS) on the three dimensional torus, forced by an additive noise. We adapt the atomic spaces framework in the context of the energy-critical nonlinear Schr\"odinger equation, and employ probabilistic perturbation arguments in the context of stochastic PDEs, establishing the global well-posedness of the defocusing energy-critical quintic SNLS in $H^1(\mathbb{T}^3)$.


[113] 2407.15112

A dilation theoretic approach to Banach spaces

For a complex Banach space $\mathbb X$, we prove that $\mathbb X$ is a Hilbert space if and only if every strict contraction $T$ on $\mathbb X$ dilates to an isometry if and only if for every strict contraction $T$ on $\mathbb X$ the function $A_T: \mathbb X \rightarrow [0, \infty]$ defined by $A_T(x)=(\|x\|^2 -\|Tx\|^2)^{\frac{1}{2}}$ gives a norm on $\mathbb X$. We also find several other necessary and sufficient conditions in this thread such that a Banach sapce becomes a Hilbert space. We construct examples of strict contractions on non-Hilbert Banach spaces that do not dilate to isometries. Then we characterize all strict contractions on a non-Hilbert Banach space that dilate to isometries and find explicit isometric dilation for them. We prove several other results including characterizations of complemented subspaces in a Banach space, extension of a Wold isometry to a Banach space unitary and describing norm attainment sets of Banach space operators in terms of dilations.


[114] 2407.15114

Generalised spherical minors and their relations

Let H be a spherical subgroup of minimal rank of the semisimple simply connected complex algebraic group G. We define some functions on the homogeneous space G/H that we call generalised spherical minors. When G = H x H, we recover Fomin-Zelevinsky generalised minors. We prove that generalised spherical minors satisfy some integer coefficients polynomial relations, that generalise the identities of classical generalised minors due to Fomin and Zelevinsky, and that have the shape of exchange relations of LP-algebras.


[115] 2407.15115

The Aharonov-Bohm Hamiltonian: self-adjointness, spectral and scattering properties

This work provides an introduction and overview on some basic mathematical aspects of the single-flux Aharonov-Bohm Schr\"odinger operator. The whole family of admissible self-adjoint realizations is characterized by means of four different methods: von Neumann theory, boundary triplets, quadratic forms and Kre{\u\i}n's resolvent formalism. The relation between the different parametrizations thus obtained is explored, comparing the asymptotic behavior of functions in the corresponding operator domains close to the flux singularity. Special attention is devoted to those self-adjoint realizations which preserve the same rotational symmetry and homogeneity under dilations of the basic differential operator. The spectral and scattering properties of all the Hamiltonian operators are finally described.


[116] 2407.15118

A Cobham theorem for scalar multiplication

Let $\alpha,\beta \in \mathbb{R}_{>0}$ be such that $\alpha,\beta$ are quadratic and $\mathbb{Q}(\alpha)\neq \mathbb{Q}(\beta)$. Then every subset of $\mathbb{R}^n$ definable in both $(\mathbb{R},{<},+,\mathbb{Z},x\mapsto \alpha x)$ and $(\mathbb{R},{<},+,\mathbb{Z},x\mapsto \beta x)$ is already definable in $(\mathbb{R},{<},+,\mathbb{Z})$. As a consequence we generalize Cobham-Semenov theorems for sets of real numbers to $\beta$-numeration systems, where $\beta$ is a quadratic irrational.


[117] 2407.15120

Designs in compact symmetric spaces and applications of great antipodal sets

The theory of designs is an important branch of combinatorial mathematics. It is well-known in the theory of designs that a finite subset of a sphere is a tight spherical 1-design if and only if it is a pair of antipodal points. On the other hand, antipodal sets and 2-number for a Riemannian manifold are introduced by B.-Y. Chen and T. Nagano in 1982. An antipodal set is called a great antipodal set if its cardinality is equal to the 2-number. The main purpose of this paper is to provide a survey on important results in compact symmetric spaces with great antipodal sets as the designs. In the last two sections of this paper, we present some important applications of 2-number and great antipodal sets to topology and group theory.


[118] 2407.15121

Squared Distance Function on the Configuration Space of a planar Spider with Applications to Hooke Energy and Voronoi Distance

Spider mechanisms are the simplest examples of arachnoid mechanisms, they are one step more complicated than polygonal linkages. Their configuration spaces have been studied intensively, but are yet not completely understood. In the paper we study them using the Morse theory of the squared distance function from the "body" of the spider to some fixed point in the plane. Generically, it is a Morse-Bott function. We list its critical manifolds, describe them as products of polygon spaces, and derive a formula for their Morse-Bott indices. We apply the obtained results to Hooke energy and Voronoi distance.


[119] 2407.15128

A description of the integral depth-$r$ Bernstein center

In this paper we give a description of the depth-$r$ Bernstein center for non-negative integers $r$ of a reductive simply connected group $G$ over a non-archimedean local field as a limit of depth-$r$ standard parahoric Hecke algebras. Using the description, we construct maps from the algebra of stable functions on the $r$-th Moy-Prasad filtration quotient of hyperspecial parahorics to the depth-$r$ Bernstein center and use them to attach to each depth-$r$ irreducible representation $\pi$ an invariant $\theta(\pi)$, called the depth-$r$ Deligne-Lusztig parameter of $\pi$. We show that $\theta(\pi)$ is equal to the semi-simple part of minimal $K$-types of $\pi$.


[120] 2407.15135

The directional short-time fractional Fourier transform of distributions

We introduce the directional short-time fractional Fourier transform (DSTFRFT) and prove an extended Parseval's identity and a reconstruction formula for it. We also investigate the continuity of both the directional short-time fractional Fourier transform and its synthesis operator on the appropriate space of test functions. Using the obtained continuity results, we develop a distributional framework for the DSTFRFT on the space of tempered distributions $\mathcal{S}'(\mathbb{R}^n)$. We end the article with a desingularization formula.


[121] 2407.15146

On an Erdős-type conjecture on $\mathbb{F}_2[x]$

We generalize Erd\H{o}s's conjecture on integer covering systems to polynomial rings on finite field, for example, $\mathbb F_2[x]$, and proved the theorem: Every system of $n$ cosets of ideals in $\mathbb F_q[x]$ covering all polynomials whose degree is less than $n$ covers $\mathbb F_q[x]$. When $q=2$, our result is sharp.


[122] 2407.15156

Computational and analytical studies of a new nonlocal phase-field crystal model in two dimensions

A nonlocal phase-field crystal (NPFC) model is presented as a nonlocal counterpart of the local phase-field crystal (LPFC) model and a special case of the structural PFC (XPFC) derived from classical field theory for crystal growth and phase transition. The NPFC incorporates a finite range of spatial nonlocal interactions that can account for both repulsive and attractive effects. The specific form is data-driven and determined by a fitting to the materials structure factor, which can be much more accurate than the LPFC and previously proposed fractional variant. In particular, it is able to match the experimental data of the structure factor up to the second peak, an achievement not possible with other PFC variants studied in the literature. Both LPFC and fractional PFC (FPFC) are also shown to be distinct scaling limits of the NPFC, which reflects the generality. The advantage of NPFC in retaining material properties suggests that it may be more suitable for characterizing liquid-solid transition systems. Moreover, we study numerical discretizations using Fourier spectral methods, which are shown to be convergent and asymptotically compatible, making them robust numerical discretizations across different parameter ranges. Numerical experiments are given in the two-dimensional case to demonstrate the effectiveness of the NPFC in simulating crystal structures and grain boundaries.


[123] 2407.15159

A priori interior estimates for special Lagrangian curvature equations

We establish a priori interior curvature estimates for the special Lagrangian curvature equations in both the critical phase and convex case. Additionally, we prove a priori interior gradient estimates for any constant phases.


[124] 2407.15162

Random walk on dynamical percolation in Euclidean lattices: separating critical and supercritical regimes

We study the random walk on dynamical percolation of $\mathbb{Z}^d$ (resp., the two-dimensional triangular lattice $\mathcal{T}$), where each edge (resp., each site) can either be open or closed, refreshing its status at rate $\mu\in (0,1/e]$. The random walk moves along open edges in $\mathbb{Z}^d$ (resp., open sites in $\mathcal{T}$) at rate $1$. For the critical regime $p=p_c$, we prove the following two results: on $\mathcal{T}$, the mean squared displacement of the random walk from $0$ to $t$ is at most $O(t\mu^{5/132-\epsilon})$ for any $\epsilon>0$; on $\mathbb{Z}^d$ with $d\geq 11$, the corresponding upper bound for the mean squared displacement is $O(t \mu^{1/2}\log(1/\mu))$. For the supercritical regime $p>p_c$, we prove that the mean squared displacement on $\mathbb{Z}^d$ is at least $ct$ for some $c=c(d)>0$ that does not depend on $\mu$.


[125] 2407.15163

Limit cycles of piecewise smooth differential systems with nilpotent center and linear saddle

In this paper, we study the number of limit cycles of a piecewise smooth differential system separated by one or two parallel straight lines or rays formed by a nilpotent center or degenerate center and linear saddle. Piecewise linear differential systems separated by one or two parallel straight lines with one of the subsystems of type nilpotent center and other subsystems of type linear saddle can have at most two limit cycles and there are systems in these classes having one limit cycle. The limit cycle in particular consists of saddle separatrices of the subsystem.


[126] 2407.15180

Generalizing Trilateration: Approximate Maximum Likelihood Estimator for Initial Orbit Determination in Low-Earth Orbit

With the increase in the number of active satellites and space debris in orbit, the problem of initial orbit determination (IOD) becomes increasingly important, demanding a high accuracy. Over the years, different approaches have been presented such as filtering methods (for example, Extended Kalman Filter), differential algebra or solving Lambert's problem. In this work, we consider a setting of three monostatic radars, where all available measurements are taken approximately at the same instant. This follows a similar setting as trilateration, a state-of-the-art approach, where each radar is able to obtain a single measurement of range and range-rate. Differently, and due to advances in Multiple-Input Multiple-Output (MIMO) radars, we assume that each location is able to obtain a larger set of range, angle and Doppler shift measurements. Thus, our method can be understood as an extension of trilateration leveraging more recent technology and incorporating additional data. We formulate the problem as a Maximum Likelihood Estimator (MLE), which for some number of observations is asymptotically unbiased and asymptotically efficient. Through numerical experiments, we demonstrate that our method attains the same accuracy as the trilateration method for the same number of measurements and offers an alternative and generalization, returning a more accurate estimation of the satellite's state vector, as the number of available measurements increases.


[127] 2407.15183

Towards a solution of Archdeacon's conjecture on integer Heffter arrays

In this paper, we make significant progress on a conjecture proposed by Dan Archdeacon on the existence of integer Heffter arrays $H(m,n;s,k)$ whenever the necessary conditions hold, that is, $3\leqslant s \leqslant n$, $3\leqslant k\leqslant m$, $ms=nk$ and $nk\equiv 0,3 \pmod 4$. By constructing integer Heffter array sets, we prove the conjecture in the affirmative whenever $k\geqslant 7\cdot \gcd(s,k)$ is odd and $s\neq 3,5,6,10$.


[128] 2407.15194

Quasi-linear elliptic equations with superlinear convection

We discuss the existence and regularity of solutions to a quasi-linear elliptic equation involving a Leray-Lions operator and a convection term with superlinear growth. In particular, equations involving the p-Laplacian are covered. This paper generalizes some of the results in Boccardo et al. (2024).


[129] 2407.15195

Exact Convergence rate of the subgradient method by using Polyak step size

This paper studies the last iterate of subgradient method with Polyak step size when applied to the minimization of a nonsmooth convex function with bounded subgradients. We show that the subgradient method with Polyak step size achieves a convergence rate $\mathcal{O}\left(\tfrac{1}{\sqrt[4]{N}}\right)$ in terms of the final iterate. An example is provided to show that this rate is exact and cannot be improved. We introduce an adaptive Polyak step size for which the subgradient method enjoys a convergence rate $\mathcal{O}\left(\tfrac{1}{\sqrt{N}}\right)$ for the last iterate. Its convergence rate matches exactly the lower bound on the performance of any black-box method on the considered problem class. Additionally, we propose an adaptive Polyak method with a momentum term, where the step sizes are independent of the number of iterates. We establish that the algorithm also attains the optimal convergence rate. We investigate the alternating projection method. We derive a convergence rate $\left( \frac{2N }{ 2N+1 } \right)^N\tfrac{R}{\sqrt{2N+1}}$ for the last iterate, where $R$ is a bound on the distance between the initial iterate and a solution. An example is also provided to illustrate the exactness of the rate.


[130] 2407.15197

On fractional inequalities on metric measure spaces with polar decomposition

In this paper, we prove the fractional Hardy inequality on polarisable metric measure spaces. The integral Hardy inequality for $1<p\leq q<\infty$ is playing a key role in the proof. Moreover, we also prove the fractional Hardy-Sobolev type inequality on metric measure spaces. In addition, logarithmic Hardy-Sobolev and fractional Nash type inequalities on metric measure spaces are presented. In addition, we present applications on homogeneous groups and on the Heisenberg group.


[131] 2407.15198

Dots and Boxes on Certain Families of Graphs

We investigate the Dots and Boxes game, also known as ``Strings and Coins,'' for certain specific families of graphs. These include complete graphs, wheel graphs, and friendship graphs.


[132] 2407.15201

$q$-weighted analogue of the Trollope-Delange formula

Let $s(n)$ denote the number of "$1$"s in the dyadic representation of a positive integer $n$ and sequence $S(n) = s(1)+s(2)+\dots+s(n-1)$. The Trollope-Delange formula is a classic result that represents the sequence $S$ in terms of the Takagi function. In this work, we extend this result by introducing a weighted analog of $s(n)$, deriving a variant of the Trollope-Delange formula for this generalization.


[133] 2407.15206

The v-numbers and linear presentations of ideals of covers of graphs

Let $G$ be a graph and let $J=I_c(G)$ be its ideal of covers. The aims of this work are to study the {\rm v}-number ${\rm v}(J)$ of $J$ and to study when $J$ is linearly presented using combinatorics and commutative algebra. We classify when ${\rm v}(J)$ attains its minimum and maximum possible values in terms of the vertex covers of the graph that satisfy the exchange property. If the cover ideal of a graph has a linear presentation, we express its v-number in terms of the covering number of the graph. We show necessary and sufficient conditions for the graph $\mathcal{G}_J$ of $J$ to be connected. One of our main theorems shows that if $G$ has no induced 4-cycles, then $J$ is linearly presented. To prove this theorem, we show three results about unmixed K\"onig graphs related to some results already in the literature. For unmixed graphs without $3$- and $5$-cycles, we classify combinatorially when $J$ is linearly presented, and show that the columns of the linear syzygy matrix of $J$ are linearly independent if and only if $\mathcal{G}_J$ has no strong $3$-cycles.


[134] 2407.15210

Iteration of Exponentials with Sign Changes

In this paper we consider the iteration of infinitely many signed exponentials with the same base but the signs may vary. We show that for every base in an explicit interval this iteration converges for any sequence of signs and all the real numbers are possible limits. We give some more results for a base outside this interval.


[135] 2407.15215

Note on C*-algebras associated to boundary actions of hyperbolic 3-manifold groups

Using Kirchberg-Phillips' classification of purely infinite C*-algebras by K-theory, we prove that the isomorphism types of crossed product C*-algebras associated to certain hyperbolic 3-manifold groups acting on their Gromov boundary only depend on the manifold's homology. As a result, we obtain infinitely many pairwise non-isomorphic hyperbolic groups all of whose associated crossed products are isomorphic. These isomomorphisms are not of dynamical nature in the sense that they are not induced by isomorphisms of the underlying groupoids.


[136] 2407.15218

The finite products of shifted primes and Moreira's Theorem

Let $r\in\mathbb{N}$ and $\mathbb{N}=\bigcup_{i=1}^{r}C_{i}$. Do there exist $x,y\in\mathbb{N}$ and $r\in\left\{1,2,\ldots,r\right\}$ such that $\left\{x,y,xy,x+y\right\}\subseteq C_{i}$? This is still an unanswered question asked by N. Hindman. Joel Moreira in [J. Moreira, Monochromatic sums and products in N, Annals of Mathematics 185 (2017) 1069-1090] established a partial answer to this question and proved that for infinitely many $x,y\in\mathbb{N}$, $\left\{x,xy,x+y\right\}\subseteq C_{i}$ for some $i\in\left\{1,2,\ldots,r\right\}$, which is called Moreira's Theorem. Recently, H. Hindman and D. Strauss established a refinement of Moreira's Theorem and proved that for infinitely many $y$, $\left\{x\in\mathbb{N}:\left\{x,xy,x+y\right\}\subseteq C_{i}\right\}$ is a piecewise syndetic set. In this article, we will prove infinitely many $y\in FP\left(\mathbb{P}-1\right)$ such that $\left\{x\in\mathbb{N}:\left\{xy,x+f(y):f\in F\right\}\subseteq C_{i}\right\}$ is piecewise syndetic, where $F$ is a finite subset of $x\mathbb{Z}\left[x\right]$. $\mathbb{P}$ is the set of prime numbers in $\mathbb{N}$ and $FP\left(\mathbb{P}-1\right)$ is the set of all finite products of $\mathbb{P}-1$.


[137] 2407.15222

Fourier coefficients and cuspidality of modular forms: a new approach

We provide a simple and new induction based treatment of the problem of distinguishing cusp forms from the growth of the Fourier coefficients of modular forms. Our approach also gives the best possible ranges of the weights for this problem, and has wide adaptability.


[138] 2407.15230

On the boundary branching set of the one-phase problem

We consider minimizers of the one-phase Bernoulli free boundary problem in domains with analytic fixed boundary. In any dimension $d$, we prove that the branching set at the boundary has Hausdorff dimension at most $d-2$. As a consequence, we also obtain an analogous estimate on the branching set for solutions to the symmetric two-phase problem. The approach we use is based on the (almost-)monotonicity of a boundary Almgren-type frequency function, obtained via regularity estimates and a Calder\'on-Zygmund decomposition in the spirit of Almgren-De Lellis-Spadaro.


[139] 2407.15232

Transport equation driven by a stochastic measure

We consider the stochastic transport equation where the randomness is given by the symmetric integral with respect to stochastic measure. For stochastic measure, we assume only $\sigma$-additivity in probability and continuity of paths. The existence and uniqueness of the weak solution to the equation are proved.


[140] 2407.15244

On disjunction convex hulls by lifting

We study the natural extended-variable formulation for the disjunction of $n+1$ polytopes in $R^d$. We demonstrate that the convex hull $D$ in the natural extended-variable space $R^{d+n}$ is given by full optimal big-M lifting (i) when $d\leq 2$ (and that it is not generally true for $d\geq 3$), and also (ii) under some technical conditions, when the polytopes have a common facet-describing constraint matrix, for arbitrary $d\geq 1$ and $n\geq 1$. Additionally, we give further results on the polyhedral structure of $D$, and we demonstrate that all facets of $D$ can be enumerated in polynomial time.


[141] 2407.15245

Weyl Calculus and Exactly Solvable Schrödinger Bridges with Quadratic State Cost

Schr\"{o}dinger bridge--a stochastic dynamical generalization of optimal mass transport--exhibits a learning-control duality. Viewed as a stochastic control problem, the Schr\"{o}dinger bridge finds an optimal control policy that steers a given joint state statistics to another while minimizing the total control effort subject to controlled diffusion and deadline constraints. Viewed as a stochastic learning problem, the Schr\"{o}dinger bridge finds the most-likely distribution-valued trajectory connecting endpoint distributional observations, i.e., solves the two point boundary-constrained maximum likelihood problem over the manifold of probability distributions. Recent works have shown that solving the Schr\"{o}dinger bridge problem with state cost requires finding the Markov kernel associated with a reaction-diffusion PDE where the state cost appears as a state-dependent reaction rate. We explain how ideas from Weyl calculus in quantum mechanics, specifically the Weyl operator and the Weyl symbol, can help determine such Markov kernels. We illustrate these ideas by explicitly finding the Markov kernel for the case of quadratic state cost via Weyl calculus, recovering our earlier results but avoiding tedious computation with Hermite polynomials.


[142] 2407.15256

Weak-instrument-robust subvector inference in instrumental variables regression: A subvector Lagrange multiplier test and properties of subvector Anderson-Rubin confidence sets

We propose a weak-instrument-robust subvector Lagrange multiplier test for instrumental variables regression. We show that it is asymptotically size-correct under a technical condition. This is the first weak-instrument-robust subvector test for instrumental variables regression to recover the degrees of freedom of the commonly used Wald test, which is not robust to weak instruments. Additionally, we provide a closed-form solution for subvector confidence sets obtained by inverting the subvector Anderson-Rubin test. We show that they are centered around a k-class estimator. Also, we show that the subvector confidence sets for single coefficients of the causal parameter are jointly bounded if and only if Anderson's likelihood-ratio test rejects the hypothesis that the first-stage regression parameter is of reduced rank, that is, that the causal parameter is not identified. Finally, we show that if a confidence set obtained by inverting the Anderson-Rubin test is bounded and nonempty, it is equal to a Wald-based confidence set with a data-dependent confidence level. We explicitly compute this Wald-based confidence test.


[143] 2407.15262

A note about discrete Riesz potential on $\mathbb{Z}^n$

In this note we prove that the discrete Riesz potential $I_{\alpha}$ defined on $\mathbb{Z}^n$ is a bounded operator $H^p (\mathbb{Z}^n) \to \ell^q (\mathbb{Z}^n)$ for $0 < p \leq 1$ and $\frac{1}{q} = \frac{1}{p} - \frac{\alpha}{n}$, where $0 < \alpha < n$.


[144] 2407.15269

A short survey on biharmonic Riemannian submersions

The study of biharmonic submanifolds, initiated by B. Y. Chen and G. Y. Jiang independently, has received a great attention in the past 30 years with many important progress. This note attempts to give a short survey on the study of biharmonic Riemannian submersions which are a dual concept of biharmonic submanifolds (i.e., biharmonic isometric immersions).


[145] 2407.15271

Fixed Point Property in G-Complete Fuzzy Metric Space

Our purpose of this paper is to focus on fixed point property in fuzzy metric space. To achieve our objective, we will introduce a new contraction condition to examine the fixed point for multi-valued mapping, then we will be investigating the obtained result to ensure the existence and uniqueness of this property for single-valued mapping. To show the use of our main result, we will give the relative result in the ordinary metric space.


[146] 2407.15274

An iterable surgery formula on involutive knot lattice homotopy

In ``Knots in lattice homology", Ozsv\'ath, Stipsicz, and Szab\'o showed that knot lattice homology satisfies a surgery formula similar to the one relating knot Floer homology and Heegaard Floer homology, and in previous work, I showed that knot lattice homology is the persistent homology of a doubly filtered space. Here I provide an iterable version of the surgery formula that, provided the initial knot lattice space with flip map, produces a space isomorphic as a doubly-filtered space to the corresponding knot lattice space for the dual knot with the corresponding flip map. If we include the involutive data for the original knot and ambient three-manifold, we can also produce the corresponding involutive data on the new knot lattice space without assuming that the original three-manifold is an $L$-space. I construct $\infty$-categories where these operations are functorial. Finally, I use the surgery formula to compute some examples of knot lattice spaces including for the regular fiber of $\Sigma(2,3,7)$ and for a knot in a three-manifold that is not given by an almost rational graph.


[147] 2407.15280

Covering systems with the sum of the reciprocals of the moduli close to $1$

In 1952, H. Davenport posed the problem of determining a condition on the minimum modulus $m_{0}$ in a finite distinct covering system that would imply that the sum of the reciprocals of the moduli in the covering system is bounded away from $1$. In 1973, P. Erdos and J. Selfridge indicated that they believed that $m_{0} > 4$ would suffice. We provide a proof that this is the case.


[148] 2407.15287

Poisson bundles over unordered configurations

In this paper we construct a Poisson algebra bundle whose distributional sections are suitable to represent multilocal observables in classical field theory. To do this, we work with vector bundles over the unordered configuration space of a manifold $M$ and consider the structure of a $2$-monoidal category given by the usual (Hadamard) tensor product of bundles and a new (Cauchy) tensor product which provides a symmetrized version of the usual external tensor product of vector bundles on $M$. We use the symmetric algebras with respect to both products to obtain a Poisson 2-algebra bundle mimicking the construction of Peierls bracket from the causal propagator in field theory. The explicit description of observables from this Poisson algebra bundle will be carried out in a forthcoming paper.


[149] 2407.15292

Input-to-state stabilization of $1$-D parabolic equations with Dirichlet boundary disturbances under boundary fixed-time control

This paper addresses the problem of stabilization of $1$-D parabolic equations with destabilizing terms and Dirichlet boundary disturbances. By using the method of backstepping and the technique of splitting, a boundary feedback controller is designed to ensure the input-to-state stability (ISS) of the closed-loop system with Dirichlet boundary disturbances, while preserving fixed-time stability (FTS) of the corresponding disturbance-free system, for which the fixed time is either determined by the Riemann zeta function or freely prescribed. To overcome the difficulty brought by Dirichlet boundary disturbances, the ISS and FTS properties of the involved systems are assessed by applying the generalized Lyapunov method. Numerical simulations are conducted to illustrate the effectiveness of the proposed scheme of control design.


[150] 2407.15297

Distributional limits of graph cuts on discretized grids

Graph cuts are among the most prominent tools for clustering and classification analysis. While intensively studied from geometric and algorithmic perspectives, graph cut-based statistical inference still remains elusive to a certain extent. Distributional limits are fundamental in understanding and designing such statistical procedures on randomly sampled data. We provide explicit limiting distributions for balanced graph cuts in general on a fixed but arbitrary discretization. In particular, we show that Minimum Cut, Ratio Cut and Normalized Cut behave asymptotically as the minimum of Gaussians as sample size increases. Interestingly, our results reveal a dichotomy for Cheeger Cut: The limiting distribution of the optimal objective value is the minimum of Gaussians only when the optimal partition yields two sets of unequal volumes, while otherwise the limiting distribution is the minimum of a random mixture of Gaussians. Further, we show the bootstrap consistency for all types of graph cuts by utilizing the directional differentiability of cut functionals. We validate these theoretical findings by Monte Carlo experiments, and examine differences between the cuts and the dependency on the underlying distribution. Additionally, we expand our theoretical findings to the Xist algorithm, a computational surrogate of graph cuts recently proposed in Suchan, Li and Munk (arXiv, 2023), thus demonstrating the practical applicability of our findings e.g. in statistical tests.


[151] 2407.15306

Distribution of the zeros of polynomials near the unit circle

We estimate the number of zeros of a polynomial in $\mathbb{C}[z]$ within any small circular disc centered on the unit circle, which improves and comprehensively extends a result established by Borwein, Erd{\'e}lyi, and Littmann~\cite{BE1} in 2008. Furthermore, by combining this result with Euclidean geometry, we derive an upper bound on the number of zeros of such a polynomial within a region resembling a gear wheel. Additionally, we obtain a sharp upper bound on the annular discrepancy of such zeros near the unit circle. Our approach builds upon a modified version of the method described in \cite{BE1}, combined with the refined version of the best-known upper bound for angular discrepancy of zeros of polynomials.


[152] 2407.15307

The equidistant dimension of some graphs of convex polytopes

This paper is devoted to some rotationally symmetric classes of graphs denoted in literature as convex polytope graphs. Exact value of equidistant dimension is found for $T_n$. Next, for even $n$ exact values are found for $R''_n$ and $S''_n$, while for odd $n$ exact value is found for $S_n$. Finally, for odd $n$, lower bound are found for $R''_n$ and $S''_n$.


[153] 2407.15315

A Fast and Accurate Solver for the Fractional Fokker-Planck Equation with Dirac-Delta Initial Conditions

The classical Fokker-Planck equation (FPE) is a key tool in physics for describing systems influenced by drag forces and Gaussian noise, with applications spanning multiple fields. We consider the fractional Fokker-Planck equation (FFPE), which models the time evolution of probability densities for systems driven by L\'evy processes, relevant in scenarios where Gaussian assumptions fail. The paper presents an efficient and accurate numerical approach for the free-space FFPE with constant coefficients and Dirac-delta initial conditions. This method utilizes the integral representation of the solutions and enables the efficient handling of very high-dimensional problems using fast algorithms. Our work is the first to present a high-precision numerical solver for the free-space FFPE with Dirac-delta initial conditions. This opens the door for future research on more complex scenarios, including those with variable coefficients and other types of initial conditions.


[154] 2407.15332

Robust personalized pricing under uncertainty of purchase probabilities

This paper is concerned with personalized pricing models aimed at maximizing the expected revenues or profits for a single item. While it is essential for personalized pricing to predict the purchase probabilities for each consumer, these predicted values are inherently subject to unavoidable errors that can negatively impact the realized revenues and profits. To address this issue, we focus on robust optimization techniques that yield reliable solutions to optimization problems under uncertainty. Specifically, we propose a robust optimization model for personalized pricing that accounts for the uncertainty of predicted purchase probabilities. This model can be formulated as a mixed-integer linear optimization problem, which can be solved exactly using mathematical optimization solvers. We also develop a Lagrangian decomposition algorithm combined with line search to efficiently find high-quality solutions for large-scale optimization problems. Experimental results demonstrate the effectiveness of our robust optimization model and highlight the utility of our Lagrangian decomposition algorithm in terms of both computational efficiency and solution quality.


[155] 2407.15342

The finite basis problem for additively idempotent semirings of order four, I

We study the finite basis problem for 4-element additively idempotent semirings whose additive reducts are semilattices of height 1. Up to isomorphism, there are 58 such algebras. We show that 49 of them are finitely based and the remaining ones are nonfinitely based.


[156] 2407.15361

A vector-host epidemic model with spatial structure and seasonality

Recently, Li and Zhao [5] (Bull. Math. Biol., 83(5), 43, 25 pp (2021)) proposed and studied a periodic reaction-diffusion model of Zika virus with seasonality and spatial heterogeneous structure in host and vector populations. They found the basic reproduction ratio R0, which is a threshold parameter. In this short paper we shall use the upper and lower solutions method to study the model of [5] with Neumann boundary conditions replaced by general boundary conditions.


[157] 2407.15364

Hybrid PDE-Deep Neural Network Model for Calcium Dynamics in Neurons

Traditionally, calcium dynamics in neurons are modeled using partial differential equations (PDEs) and ordinary differential equations (ODEs). The PDE component focuses on reaction-diffusion processes, while the ODE component addresses transmission via ion channels on the cell's or organelle's membrane. However, analytically determining the underlying equations for ion channels is highly challenging due to the complexity and unknown factors inherent in biological processes. Therefore, we employ deep neural networks (DNNs) to model the open probability of ion channels, a task that can be intricate when approached with ODEs. This technique also reduces the number of unknowns required to model the open probability. When trained with valid data, the same neural network architecture can be used for different ion channels, such as sodium, potassium, and calcium. Furthermore, based on the given data, we can build more physiologically reasonable DNN models that can be customized. Subsequently, we integrated the DNN model into calcium dynamics in neurons with endoplasmic reticulum, resulting in a hybrid model that combines PDEs and DNNs. Numerical results are provided to demonstrate the flexibility and advantages of the PDE-DNN model.


[158] 2407.15365

Pseudo-Energy-Preserving Explicit Runge-Kutta Methods

Using a recent characterization of energy-preserving B-series, we derive the explicit conditions on the coefficients of a Runge-Kutta method that ensure energy preservation (for Hamiltonian systems) up to a given order in the step size, which we refer to as the pseudo-energy-preserving (PEP) order. We study explicit Runge-Kutta methods with PEP order higher than their classical order. We provide examples of such methods up to PEP order six, and test them on Hamiltonian ODE and PDE systems. We find that these methods behave similarly to exactly energy-conservative methods over moderate time intervals and exhibit significantly smaller errors, relative to other Runge-Kutta methods of the same order, for moderately long-time simulations.


[159] 2407.15367

Two edge contractions for the affine super Yangian

In the previous paper, we constructed two kinds of edge contractions for the affine super Yangian. In this article, we show that these two edge contractions are commutative with each other. As an application, we give a homomorphism from the affine super Yangian to some centralizer algebras of the universal enveloping algebra of $W$-algebras of type $A$.


[160] 2407.15368

Nearly Optimal $L_p$ Risk Minimization

Convex risk measures play a foundational role in the area of stochastic optimization. However, in contrast to risk neutral models, their applications are still limited due to the lack of efficient solution methods. In particular, the mean $L_p$ semi-deviation is a classic risk minimization model, but its solution is highly challenging due to the composition of concave-convex functions and the lack of uniform Lipschitz continuity. In this paper, we discuss some progresses on the design of efficient algorithms for $L_p$ risk minimization, including a novel lifting reformulation to handle the concave-convex composition, and a new stochastic approximation method to handle the non-Lipschitz continuity. We establish an upper bound on the sample complexity associated with this approach and show that this bound is not improvable for $L_p$ risk minimization in general.


[161] 2407.15372

TU-games with utilities: the prenucleolus and its characterization set

TU-games with utility functions are considered. Generalizations of the prenucleolus, essential coalitions and the core: the u-prenucleolus, u-essential coalitions and the u-core respectively are introduced. We show that u-essential coalitions form a characterisation set for the u-prenucleolus in case of games with nonempty u-core.


[162] 2407.15381

Prismatic Crystals for schemes in characteristic $p$

In this paper, we consider the (crystalline) prismatic crystals on a scheme $\mathfrak{X}$. We classify the crystals by $p$-connections on a certain ring and prove a cohomological comparison theorem. This equivalence is more general than the Ogus' equivalence. Moreover, we will prove the equivalence between the crystalline crystals and prismatic crystals by gluing local Ogus-Vologodsky correspondence.


[163] 2407.15384

Inversion Diameter and Treewidth

In an oriented graph $\overrightarrow{G}$, the inversion of a subset $X$ of vertices consists in reversing the orientation of all arcs with both end-vertices in $X$. The inversion graph of a graph $G$, denoted by $\mathcal{I}(G)$, is the graph whose vertices are orientations of $G$ in which two orientations $\overrightarrow{G_1}$ and $\overrightarrow{G_2}$ are adjacent if and only if there is an inversion $X$ transforming $\overrightarrow{G_1}$ into $\overrightarrow{G_2}$. The inversion diameter of a graph $G$ is the diameter of its inversion graph $\mathcal{I}(G)$ denoted by $diam(\mathcal{I}(G))$. Havet, H\"orsch, and Rambaud~(2024) first proved that for $G$ of treewidth $k$, $diam(\mathcal{I}(G)) \le 2k$, and there are graphs of treewidth $k$ with inversion diameter $k+2$. In this paper, we construct graphs of treewidth $k$ with inversion diameter $2k$, which implies that the previous upper bound $diam(\mathcal{I}(G)) \le 2k$ is tight. Moreover, for graphs with maximum degree $\Delta$, Havet, H\"orsch, and Rambaud~(2024) proved $diam(\mathcal{I}(G)) \le 2\Delta-1$ and conjectured that $diam(\mathcal{I}(G)) \le \Delta$. We prove the conjecture when $\Delta=3$ with the help of computer calculations.


[164] 2407.15393

On some recent quasi-copula problems and some new methods

The aim of this paper is to present three construction methods for quasi-copulas based on recent developments: a representation of multivariate quasi-copulas by means of infima and suprema of copulas, an extension of a classical result on shuffles of min to the setting of quasi-copulas, and a construction method for quasi-copulas obeying a given signed mass pattern on a patch.


[165] 2407.15400

A discrete crystal model in three dimensions: the line-tension limit for dislocations

We propose a discrete lattice model of the energy of dislocations in three-dimensional crystals which properly accounts for lattice symmetry and geometry, arbitrary harmonic interatomic interactions, elastic deformations and discrete crystallographic slip on the full complement of slip systems of the crystal class. Under the assumption of diluteness, we show that the discrete energy converges, in the sense of $\Gamma$-convergence, to a line-tension energy defined on Volterra line dislocations, regarded as integral vector-valued currents supported on rectifiable curves. Remarkably, the line-tension limit is of the same form as that derived from semi-discrete models of linear elastic dislocations based on a core cutoff regularization. In particular, the line-tension energy follows from a cell relaxation and differs from the classical ansatz, which is quadratic in the Burgers vector.


[166] 2407.15418

Mixing and CLT for Hénon-Sibony maps: plurisubharmonic observables

Let $f$ be a complex H\'enon map and $\mu$ its unique measure of maximal entropy. We prove that $\mu$ is exponentially mixing of all orders for all (not necessarily bounded) plurisubharmonic observables, and that all plurisubharmonic functions satisfy the central limit theorem with respect to $\mu$. Our results hold more generally for every H\'enon-Sibony map on $\mathbb{C}^k$.


[167] 2407.15430

Reduced theory of symmetric and antisymmetric exchange interactions in nanowires

We investigate the behavior of minimizers of perturbed Dirichlet energies supported on a wire generated by a regular simple curve $\gamma$ and defined in the space of $\mathbb{S}^2$-valued functions. The perturbation $K$ is represented by a matrix-valued function defined on $\mathbb{S}^2$ with values in $\mathbb{R}^{3 \times 3}$. Under natural regularity conditions on $K$, we show that the family of perturbed Dirichlet energies converges, in the sense of $\Gamma$-convergence, to a simplified energy functional on $\gamma$. The reduced energy unveils how part of the antisymmetric exchange interactions contribute to an anisotropic term whose specific shape depends on the curvature of $\gamma$. We also discuss the significant implications of our results for studies of ferromagnetic nanowires when Dzyaloshinskii-Moriya interaction (DMI) is present.


[168] 2407.15432

Elliptic curves and the residue-counts of $x^2+bx+c/x$ modulo $p$

For any prime $p>3$ and rational $p$-integers $b,c$ with $c(b^3-27c)\not\equiv 0\pmod p$ let $V_p(x^2+bx+\frac cx)$ be the residue-counts of $x^2+bx+\frac cx$ modulo $p$ as $x$ runs over $1,2,\ldots,p-1$. In this paper, we reveal the connection between $V_p(x^2+bx+\frac cx)$ and the number of points on certain elliptic curve over the field $\Bbb F_p$.


[169] 2407.15434

The Burgers equation driven by a stochastic measure

We study the class of one-dimensional equations driven by a stochastic measure $\mu$. For $\mu$ we assume only $\sigma$-additivity in probability. This class of equations include the Burgers equation and the heat equation. The existence and uniqueness of the solution are proved, and the averaging principle for the equation is studied.


[170] 2407.15437

The classification of links up to clasp-pass moves

We give a complete classification of links up to clasp-pass moves, which coincides with Habiro's $C_3$-equivalence. We also classify links up to band-pass and band-# moves, which are versions of the usual pass- and #-move, respectively, where each pair of parallel strands belong to the same component. This recovers and generalizes widely a number of partial results in the study of these local moves. The proofs make use of clasper theory.


[171] 2407.15444

Prime and maximal ideals in Hurwitz polynomial rings

In this paper we study prime and maximal ideals in a Hurwitz polynomial ring hR. It is well-known that to study many questions we may assume R is prime and consider just R-disjoint ideals. We give a characterization for an R-disjoint ideal to be prime. We study conditions under which there exists an R-disjoint ideal which is a maximal ideal and when this is the case how to determine all such maximal ideals. maximal ideal to be generated by polynomials of minimal degree.


[172] 2407.15449

Persistence-based Modes Inference

We consider the estimation of multiple modes of a (multivariate) density. We start by proposing an estimator of the $H_0$ persistence diagram. We then derive from it a procedure to estimate the number of modes, their locations and the associated local maxima. For large classes of piecewise-continuous functions, we show that these estimators achieve nearly minimax rates. These classes involve geometric control over the discontinuities set and differ from commonly considered function classes in mode(s) inference. Although the global regularity assumptions are stronger, we do not suppose regularity (or even continuity) in any neighborhood of the modes.


[173] 2407.15454

The Dowker theorem via discrete Morse theory

The Dowker theorem is a classical result in the topology of finite spaces, claiming that any binary relation between two finite spaces defines two homotopy-equivalent complexes (the Dowker complexes). Recently, Barmak strengthened this to a simple-homotopy-equivalence. We reprove Barmak's result using a combinatorial argument that constructs an explicit acyclic matching in the sense of discrete Morse theory.


[174] 2407.15456

Approximation of maps between real algebraic varieties

A nonsingular real algebraic variety Y is said to have the approximation property if for every real algebraic variety X the following holds: if f:X-->Y is a C^inf map that is homotopic to a regular map, then f can be approximated in the C^inf topology by regular maps. In this paper, we characterize the varieties Y with the approximation property. We also characterize the varieties Y with the approximation property combined with a suitable interpolation condition. Some of our results have variants concerning the regular approximation of continuous maps defined on (possibly singular) real algebraic varieties.


[175] 2407.15457

Cross-diffusion systems coupled via a moving interface

We propose and study a one-dimensional model which consists of two cross-diffusion systems coupled via a moving interface. The motivation stems from the modelling of complex diffusion processes in the context of the vapor deposition of thin films. In our model, cross-diffusion of the various chemical species can be respectively modelled by a size-exclusion system for the solid phase and the Stefan-Maxwell system for the gaseous phase. The coupling between the two phases is modelled by linear phase transition laws of Butler-Volmer type, resulting in an interface evolution. The continuous properties of the model are investigated, in particular its entropy variational structure and stationary states. We introduce a two-point flux approximation finite volume scheme. The moving interface is addressed with a moving-mesh approach, where the mesh is locally deformed around the interface. The resulting discrete nonlinear system is shown to admit a solution that preserves the main properties of the continuous system, namely: mass conservation, nonnegativity, volume-filling constraints, decay of the free energy and asymptotics. In particular, the moving-mesh approach is compatible with the entropy structure of the continuous model. Numerical results illustrate these properties and the dynamics of the model.


[176] 2407.15460

Invariance Times Transfer Properties

Invariance times are stopping times $\tau$ such that local martingales with respect to some reduced filtration and an equivalently changed probability measure, stopped before $\tau$ , are local martingales with respect to the original model filtration and probability measure. They arise naturally for modeling the default time of a dealer bank, in the mathematical finance context of counterparty credit risk. Assuming an invariance time endowed with an intensity and a positive Az{\'e}ma supermartingale, this work establishes a dictionary relating the semimartingale calculi in the original and reduced stochastic bases, regarding in particular conditional expectations, martingales, stochastic integrals, random measure stochastic integrals, martingale representation properties, semimartingale characteristics, Markov properties, transition semigroups and infinitesimal generators, and solutions of backward stochastic differential equations.


[177] 2407.15463

Integrated Access and Backhaul (IAB) in Low Altitude Platforms

In this paper, we explore the problem of utilizing Integrated Access and Backhaul (IAB) technology in Non-Terrestrial Networks (NTN), with a particular focus on aerial access networks. We consider an Uncrewed Aerial Vehicle (UAV)-based wireless network comprised of two layers of UAVs: (a) a lower layer consisting a number of flying users and a UAV Base Station (BS) that provides coverage for terrestrial users and, (b) an upper layer designated to provide both wireless access for flying users and backhaul connectivity for UAV BS. By adopting IAB technology, the backhaul and access links collaboratively share their resources, enabling aerial backhauling and the utilization of the same infrastructure and frequency resources for access links. A sum-rate maximization problem is formulated by considering aerial backhaul constraints to optimally allocate the frequency spectrum between aerial and terrestrial networks. We decompose the resulting non-convex optimization problem into two sub-problems of beamforming and spectrum allocation and then propose efficient solutions for each. Numerical results in different scenarios yield insightful findings about the effectiveness of using the IAB technique in aerial networks.


[178] 2407.15468

Efficient influence functions for Sobol' indices under two designs of experiments

In this note, we are interested in the asymptotic efficiency of Sobol' indices esti-mators. After recalling the basis of asymptotic efficiency, we compute the efficientinfluence functions for Sobol' indices in two different contexts: the Pick-Freeze andthe given-data settings.


[179] 2407.15471

Convergence of the Iterates for Momentum and RMSProp for Local Smooth Functions: Adaptation is the Key

Both accelerated and adaptive gradient methods are among state of the art algorithms to train neural networks. The tuning of hyperparameters is needed to make them work efficiently. For classical gradient descent, a general and efficient way to adapt hyperparameters is the Armijo backtracking. The goal of this work is to generalize the Armijo linesearch to Momentum and RMSProp, two popular optimizers of this family, by means of stability theory of dynamical systems. We establish convergence results, under the Lojasiewicz assumption, for these strategies. As a direct result, we obtain the first guarantee on the convergence of the iterates for RMSProp, in the non-convex setting without the classical bounded assumptions.


[180] 2407.15493

On rigidity of conformal submersions

Conformal submersions are generalizations of Riemannian submersions where the differential of the submersion map is a conformal isometry when restricted to the horizontal distribution. Riemannian submersions have been instrumental in the construction of new Riemannian metrics from existing ones. Therefore, being generalizations of Riemannian submersions, conformal submersions are potential tools for constructing new examples of (special) Riemannian manifolds. A conformal submersion is said to be rigid if it reduces to a Riemannian submersion up to homothety. In this paper, we study curvature conditions that force conformal submersions to be rigid. Moreover, under suitable circumstances, using these rigidity criteria for conformal submersions, we conclude the rigidity of certain quasi-Einstein metrics.


[181] 2407.15494

On the Particle Approximation of Lagged Feynman-Kac Formulae

In this paper we examine the numerical approximation of the limiting invariant measure associated with Feynman-Kac formulae. These are expressed in a discrete time formulation and are associated with a Markov chain and a potential function. The typical application considered here is the computation of eigenvalues associated with non-negative operators as found, for example, in physics or particle simulation of rare-events. We focus on a novel \emph{lagged} approximation of this invariant measure, based upon the introduction of a ratio of time-averaged Feynman-Kac marginals associated with a positive operator iterated $l \in\mathbb{N}$ times; a lagged Feynman-Kac formula. This estimator and its approximation using Diffusion Monte Carlo (DMC) have been extensively employed in the physics literature. In short, DMC is an iterative algorithm involving $N\in\mathbb{N}$ particles or walkers simulated in parallel, that undergo sampling and resampling operations. In this work, it is shown that for the DMC approximation of the lagged Feynman-Kac formula, one has an almost sure characterization of the $\mathbb{L}_1$-error as the time parameter (iteration) goes to infinity and this is at most of $\mathcal{O}(\exp\{-\kappa l\}/N)$, for $\kappa>0$. In addition a non-asymptotic in time, and time uniform $\mathbb{L}_1-$bound is proved which is $\mathcal{O}(l/\sqrt{N})$. We also prove a novel central limit theorem to give a characterization of the exact asymptotic in time variance. This analysis demonstrates that the strategy used in physics, namely, to run DMC with $N$ and $l$ small and, for long time enough, is mathematically justified. Our results also suggest how one should choose $N$ and $l$ in practice. We emphasize that these results are not restricted to physical applications; they have broad relevance to the general problem of particle simulation of the Feynman-Kac formula.


[182] 2407.15495

Ground states of a coupled pseudo-relativistic Hartree system: existence and concentration behavior

This paper is concerned with the ground states of a coupled pseudo-relativistic Hartree system in $\mathbb{R} ^{3} $ with trapping potentials, where the intraspecies and the interspecies interaction are both attractive. By investigating an associated constraint minimization problem, the existence and non-existence of ground states are classified completely. Under certain conditions on the trapping potentials, we present a precise analysis on the concentration behavior of the minimizers as the coupling coefficient goes to a critical value, where the minimizers blow up and the maximum point sequence concentrates at a global minima of the associated trapping potentials. We also identify an optimal blowing up rate under polynomial potentials by establishing some delicate estimates of energy functionals.


[183] 2407.15503

RGD-systems over $\mathbb{F}_2$

In this paper we prove that an RGD-system over $\mathbb{F}_2$ with prescribed commutation relations exists if and only if the commutation relations are Weyl-invariant and can be realized in the group $U_+$. This result gives us a machinery to produce new examples of RGD-systems with complicated commutation relations. We also discuss some applications of this result.


[184] 2407.15505

Global well-posedness of space-time fractional diffusion equation with Rockland operator on graded Lie group

In this article, we examine the general space-time fractional diffusion equation for left-invariant hypoelliptic homogeneous operators on graded Lie groups. Our study covers important examples such as the time-fractional diffusion equation, the space-time fractional diffusion equation when diffusion is under the influence of sub-Laplacian on the Heisenberg group, or general stratified Lie groups. We establish the global well-posedness of the Cauchy problem for the general space-time fractional diffusion equation for the Rockland operator on a graded Lie group in the associated Sobolev spaces. More precisely, we establish the existence and uniqueness results for both homogeneous and inhomogeneous fractional diffusion equations. In addition, we also develop some regularity estimates.


[185] 2407.15506

Construction of Commutator Blueprints

Commutator blueprints can be seen as blueprints for constructing RGD-systems over $\mathbb{F}_2$ with prescribed commutation relations. In this paper we construct several families of Weyl-invariant commutator blueprints, mostly of universal type. Together with the main result of \cite{BiRGD} we obtain new examples of exotic RGD-systems of universal type over $\mathbb{F}_2$.


[186] 2407.15515

Saddle-node bifurcations for concave in measure and d-concave in measure skewproduct flows with applications to population dynamics and circuits

Concave in measure and d-concave in measure nonautonomous scalar ordinary differential equations given by coercive and time-compactible maps have similar properties to equations satisfying considerably more restrictive hypotheses. This paper describes the generalized simple or double saddle-node bifurcation diagrams for one-parametric families of equations of these types, from which the dynamical possibilities for each of the equations follow. This new framework allows the analysis of ``almost stochastic" equations, whose coefficients vary in very large chaotic sets. The results also apply to the analysis of the occurrence of critical transitions for a range of models much larger than in previous approaches.


[187] 2407.15517

Well-posedness of the Stokes equations on a wedge with Navier-slip boundary conditions

We consider the incompressible and stationary Stokes equations on an infinite two-dimensional wedge with non-scaling invariant Navier-slip boundary conditions. We prove well-posedness and higher regularity of the Stokes problem in a certain class of weighted Sobolev spaces. The novelty of this work is the occurrence of two different scalings in the boundary condition, which is not treated so far for the Stokes system in unbounded wedge-type domains. These difficulties are overcome by first constructing a variational solution in a second order weighted Sobolev space and subsequently proving higher regularity up to the tip of the wedge by employing an iterative scheme. We believe that this method can be used for other problems with variational structure and multiple scales.


[188] 2407.15521

Phase space analysis of higher-order dispersive equations with point interactions

We investigate nonlinear, higher-order dispersive equations with measure (or even less regular) potentials and initial data with low regularity. Our approach is of distributional nature and relies on the phase space analysis (via Gabor wave packets) of the corresponding fundamental solution - in fact, locating the modulation/amalgam space regularity of such generalized Fresnel-type oscillatory functions is a problem of independent interest in harmonic analysis.


[189] 2407.15533

Self-repellent branching random walk

We consider a discrete-time binary branching random walk with independent standard normal increments subject to a penalty $\b$ for every pair of particles that get within distance $\e$ of each other at any time. We give a precise description of the most likely configurations of the particles under this law for $N$ large and $\b,\e$ fixed. Particles spread out over a distance $2^{2N/3}$, essentially in finite time, and subsequently arrange themselves so that at time $2N/3$ they cover a grid of width $\e$ with one particle per site. After time $2N/3$, the bulk of the particles and their descendants do not move anymore, while the particles in a boundary layer of width $2^{N/3}$ form a ``staircase" to the particles in the bulk. At time $N$, each site in the boundary layer is occupied by $2^{N/3}$ particles.


[190] 2407.15538

Representation theory of very non-standard quantum $so(2N-1)$

We classify the finite dimensional representations of the quantum symmetric pair coideal subalgebra $B_{\mathbf c}$ of type $DII$ corresponding to the symmetric pair $(so(2N),so(2N-1))$. For $B_{\mathbf c}$ defined over an arbitrary field $k$ and $q\in k$ not a root of unity we establish a one-to-one correspondence between finite dimensional, simple $B_{\mathbf c}$-modules and dominant integral weights for $so(2N-1)$. We use specialisation to show that the category of finite dimensional $B_{\mathbf c}$-modules is semisimple if $\mathrm{char}(k)=0$ and $q$ is transcendental over ${\mathbb Q}$. In this case the characters of simple $B_{\mathbf c}$-modules are given by Weyl's character formula. This means in particular that the quantum symmetric pair of type $DII$ can be used to obtain Gelfand-Tsetlin bases for irreducible representations of the Drinfeld-Jimbo quantum group $U_q(so(2N))$.


[191] 2407.15542

Fast second-order dynamics with slow vanishing damping approaching the zeros of a monotone and continuous operator

In this work, we approach the problem of finding the zeros of a continuous and monotone operator through a second-order dynamical system with a damping term of the form $1/t^{r}$, where $r\in [0, 1]$. The system features the time derivative of the operator evaluated along the trajectory, which is a Hessian-driven type damping term when the governing operator comes from a potential. Also entering the system is a time rescaling parameter $\beta(t)$ which satisfies a certain growth condition. We derive $o\left(\frac{1}{t^{2r}\beta(t)}\right)$ convergence rates for the norm of the operator evaluated along the generated trajectories as well as for a gap function which serves as a measure of optimality for the associated variational inequality. The parameter $r$ enters the growth condition for $\beta(t)$: when $r < 1$, the damping $1/t^{r}$ approaches zero at a slower speed than Nesterov's $1/t$ damping; in this case, we are allowed to choose $\beta(t)$ to be an exponential function, thus having linear convergence rates for the involved quantities. We also show weak convergence of the trajectories towards zeros of the governing operator. Through a particular choice for the operator, we establish a connection with the problem of minimizing a smooth and convex function with linear constraints. The convergence rates we derived in the operator case are inherited by the objective function evaluated at the trajectories and for the feasibility gap. We also prove weak convergence of the trajectories towards primal-dual solutions of the problem. A discretization of the dynamical system yields an implicit algorithm that exhibits analogous convergence properties to its continuous counterpart. We complement our theoretical findings with two numerical experiments.


[192] 2407.15547

On a space of functions with entire Laplace transforms and its connection with the optimality of the Ingham-Karamata theorem

We study approximation properties of the Fr\'{e}chet space of all continuously differentiable functions $\tau$ such that $\tau'(x)=o(1)$ and such that their Laplace transforms admit entire extensions to $\mathbb{C}$. As an application, these approximation results are combined with the open mapping theorem to show the optimality theorem for the Ingham-Karamata Tauberian theorem.


[193] 2407.15548

Correspondences on Riemann surfaces and non-uniform hyperbolicity

We consider certain analytic correspondences on a Riemann surface, and show that they admit a weak form of expansion. In terms of their algebraic encoding by bisets, this translates to contraction of group elements along sequences arising from iterated lifting. As an application, we show that for every non-exceptional rational map on $\mathbb{P}^1$ with $4$ post-critical points, there is a finite collection of isotopy classes of curves into which every curve eventually lands under iterated lifting.


[194] 2407.15557

Hierarchical Alternating Least Squares Methods for Quaternion Nonnegative Matrix Factorizations

In this report, we discuss a simple model for RGB color and polarization images under a unified framework of quaternion nonnegative matrix factorization (QNMF) and present a hierarchical nonnegative least squares method to solve the factor matrices. The convergence analysis of the algorithm is discussed as well. We test the proposed method in the polarization image and color facial image representation. Compared to the state-of-the-art methods, the experimental results demonstrate the effectiveness of the hierarchical nonnegative least squares method for the QNMF model.


[195] 2407.15559

Optimal synthesis control for evolution equations subject to nonlocal inputs

We consider the linear quadratic (LQ) optimal control problem for a class of evolution equations in infinite dimensions, in the presence of distributed and nonlocal inputs. Following the perspective taken in our previous research work on the LQ problem for integro-differential equations, where the memory term -- here involving the control actions -- is seen as a component of the state, we offer a full (closed-loop, Riccati-like) solution to the optimization problem.


[196] 2407.15561

Non-semisimple $\mathfrak{sl}_2$ quantum invariants of fibred links

The Akutsu-Deguchi-Ohtsuki (ADO) invariants are the most studied quantum link invariants coming from a non-semisimple tensor category. We show that, for fibered links in $S^3$, the degree of the ADO invariant is determined by the genus and the top coefficient is a root of unity. More precisely, we prove that the top coefficient is determined by the Hopf invariant of the plane field of $S^3$ associated to the fiber surface. Our proof is based on the genus bounds established in our previous work, together with a theorem of Giroux-Goodman stating that fiber surfaces in the three-sphere can be obtained from a disk by plumbing/deplumbing Hopf bands.


[197] 2407.15564

Non-parametric estimation of conditional quantiles for time series with heavy tails

We propose a modified weighted Nadaraya-Watson estimator for the conditional distribution of a time series with heavy tails. We establish the asymptotic normality of the proposed estimator. Simulation study is carried out to assess the performance of the estimator. We illustrate our method using a dataset.


[198] 2407.15570

Hybrid STAR-RIS Enabled Integrated Sensing and Communication

Integrated sensing and communication (ISAC) is recognized as one of the key enabling technologies for sixth-generation (6G) wireless communication networks, facilitating diverse emerging applications and services in an energy and cost-efficient manner. This paper proposes a multi-user multi-target ISAC system to enable full-space coverage for communication and sensing tasks. The proposed system employs a hybrid simultaneous transmission and reflection reconfigurable intelligent surface (STAR-RIS) comprising active transmissive and passive reflective elements. In the proposed scheme, the passive reflective elements support communication and sensing links for nearby communication users and sensing targets, while low-power active transmissive elements are deployed to improve sensing performance and overcome high path attenuation due to multi-hop transmission for remote targets. Moreover, to optimize the transmissive/reflective coefficients of the hybrid STAR-RIS, a semi-definite relaxation (SDR)-based algorithm is proposed. Furthermore, to evaluate sensing performance, signal-to-interference-noise ratio (SINR) and Cramer-Rao bound (CRB) metrics have been derived and investigated via conducting extensive computer simulations.


[199] 2407.15571

Numerical semigroups from rational matrices II: matricial dimension does not exceed multiplicity

We continue our study of exponent semigroups of rational matrices. Our main result is that the matricial dimension of a numerical semigroup is at most its multiplicity (the least generator), greatly improving upon the previous upper bound (the conductor). For many numerical semigroups, including all symmetric numerical semigroups, our upper bound is tight.


[200] 2407.15576

Curvature-dimension condition, rigidity theorems and entropy differential inequalities on Riemannian manifolds

In this paper, we use the information-theoretic approach to study curvature-dimension condition, rigidity theorems and entropy differential inequalities on Riemannian manifolds. We prove the equivalence of the ${\rm CD}(K, m)$-condition for $K\in \mathbb{R}$ and $m\in [n, \infty]$ and some entropy differential inequalities along the geodesics on the Wasserstein space over a Riemannian manifold. The rigidity models of the entropy differential inequalities are the $K$-Einstein manifolds and the $(K, m)$-Einstein manifolds with Hessian solitons. Moreover, we prove the monotonicity and rigidity theorem of the $W$-entropy along the geodesics on the Wasserstein space over Riemannian manifolds with CD$(0, m)$-condition. Finally, we compare our work with the synthetic geometry developed by Lott-Villani and Sturm.


[201] 2407.15578

The distance function to a finite set is a topological Morse function

In this short note, we show that the distance function to any finite set $X\subset \mathbb{R}^n$ is a topological Morse function, regardless of whether $X$ is in general position. We also precisely characterize its topological critical points and their indices, and relate them to the differential critical points of the function.


[202] 2407.15579

Sharp concentration phenomena in high-dimensional Orlicz balls

In this article, we present a precise deviation formula for the intersection of two Orlicz balls generated by Orlicz functions $V$ and $W$. Additionally, we establish a (quantitative) central limit theorem in the critical case and a strong law of large numbers for the "$W$-norm" of the uniform distribution on $\mathbb{B}^{(n,V)}$. Our techniques also enable us to derive a precise formula for the thin-shell concentration of uniformly distributed random vectors in high-dimensional Orlicz balls. In our approach we establish an Edgeworth-expansion using methods from harmonic analysis together with an exponential change of measure argument.


[203] 2407.15584

Large deviations for generalized backward stochastic differential equations

This work concerns generalized backward stochastic differential equations, which are coupled with a family of reflecting diffusion processes. First of all, we establish the large deviation principle for forward stochastic differential equations with reflecting boundaries under weak monotonicity conditions. Then based on the obtained result and the contraction principle, the large deviation principle for the generalized backward stochastic differential equations is proved. As a by-product, we obtain a limit result about parabolic partial differential equations with the nonlinear Neumann boundary conditions.


[204] 2407.15585

Competing DEA procedures: analysis, testing, and comparisons

Reducing the computational time to process large data sets in Data Envelopment Analysis (DEA) is the objective of many studies. Contributions include fundamentally innovative procedures, new or improved preprocessors, and hybridization between - and among - all these. Ultimately, new contributions are made when the number and size of the LPs solved is somehow reduced. This paper provides a comprehensive analysis and comparison of two competing procedures to process DEA data sets: BuildHull and Enhanced Hierarchical Decomposition (EHD). A common ground for comparison is made by examining their sequential implementations, applying to both the same preprocessors - when permitted - on a suite of data sets widely employed in the computational DEA literature. In addition to reporting on execution time, we discuss how the data characteristics affect performance and we introduce using the number and size of the LPs solved to better understand performances and explain differences. Our experiments show that the dominance of BuildHull can be substantial in large-scale and high-density datasets. Comparing and explaining performance based on the number and size of LPS lays the groundwork for a comparison of the parallel implementations of procedures BuildHull and EHD.


[205] 2407.15598

Shifted symplectic structure on Poisson Lie algebroid and generalized complex geometry

Generalized complex geometry was classically formulated by the language of differential geometry. In this paper, we reformulated a generalized complex manifold as a holomorphic symplectic differentiable formal stack in a homotopical sense. Meanwhile, by developing the machinery for shifted symplectic formal stack, we prove that the coisotropic intersection inherits shifted Poisson structure. Generalized complex branes are also studied.


[206] 2407.15601

Mokobodzki's intervals: an approach to Dynkin games when value process is not a semimartingale

We study Dynkin games governed by a nonlinear $\mathbb E^f$-expectation on a finite interval $[0,T]$, with payoff c\`adl\`ag processes $L,U$ of class (D) which are not imposed to satisfy (weak) Mokobodzki's condition - the existence of a c\`adl\`ag semimartingale between the barriers. For that purpose we introduce the notion of Mokobodzki's stochastic intervals $\mathscr M(\theta)$ (roughly speaking, maximal stochastic interval on which Mokobodzki's condition is satisfied when starting from the stopping time $\theta$) and the notion of reflected BSDEs without Mokobodzki's condition (this is a generalization and modification of the notion introduced by Hamad\'ene and Hassani (2005)). We prove an existence and uniqueness result for RBSDEs with driver $f$ that is non-increasing with respect to the value variable (no restrictions on the growth) and Lipschitz continuous with respect to the control variable, and with data in $L^1$ spaces. Next, by using RBSDEs, we show numerous results on Dynkin games: existence of the value process, saddle points, and convergence of the penalty scheme. We also show that the game is not played beyond $\mathscr M(\theta)$, when starting from $\theta$.


[207] 2407.15607

Constructions of Waldhausen categories via Grothendieck opfibrations

Given a Grothendieck opfibration $p: \mathcal{T} \to \mathcal{B}$, we describe a method to construct a Waldhausen category structure on the total category $\mathcal{T}$ via combining Waldhausen category structures on the fibers $\mathcal{T}_A$ for $A \in \mathrm{Ob}(\mathcal{B})$ and the basis category $\mathcal{B}$. As an application, we show that if $\mathsf{E}$ is a Waldhausen category with small coproducts such that the class of cofibrations is the left part of a weak factorization system in $\mathsf{E}$, then the representation category $\mathsf{Rep}(Q, \mathsf{coE})$ of a left rooted quiver $Q$ is a Waldhausen category, where $\mathsf{coE}$ is the subcategory of $\mathsf{E}$ whose morphisms are cofibrations.


[208] 2407.15619

Fractional Poisson Random Fields on $\mathbb{R}^2_+$

In this paper, we consider a fractional Poisson random field (FPRF) on positive plane. It is defined as a process whose one dimensional distribution is the solution of a system of fractional partial differential equations. A time-changed representation for the FPRF is given in terms of the composition of Poisson random field with a bivariate random process. Some integrals of the FPRF are introduced and studied. Using the Adomian decomposition method, a closed form expression for its probability mass function is obtained in terms of the generalized Wright function. Some results related to the order statistics of random numbers of random variables are presented. Also, we introduce a generalization of Poisson random field on $\mathbb{R}^d_+$, $d\ge1$ which reduces to the Poisson random field in a special case. For $d=1$, it further reduces to a generalized Poisson process (GPP). A time-changed representation for the GPP is established. Moreover, we construct a time-changed linear and planer random motions of a particle driven by the GPP. The conditional distribution of the random position of particle is derived. Later, we define the compound fractional Poisson random field via FPRF. Also, a generalized version of it on $\mathbb{R}^d_+$, $d\ge1$ is discussed.


[209] 2407.15628

A note on congruences for generalized cubic partitions modulo primes

Recently, Amdeberhan, Sellers and Singh introduced the notion of a generalized cubic partition function $a_c(n)$ and proved two isolated congruences via modular forms, namely $a_3(7n+4)\equiv 0\pmod{7}$ and $a_5(11n+10)\equiv 0\pmod{11}$. In this paper, we prove these congruences using $q$-series manipulations and provide additional congruences for $a_c(n)$ modulo certain primes via modular forms.


[210] 2407.15632

Partial Difference Sets with Denniston Parameters in Elementary Abelian $p$-Groups

Denniston \cite{D1969} constructed partial difference sets (PDS) with parameters $(2^{3m}, (2^{m+r}-2^m+2^r)(2^m-1), 2^m-2^r+(2^{m+r}-2^m+2^r)(2^r-2), (2^{m+r}-2^m+2^r)(2^r-1))$ in elementary abelian groups of order $2^{3m}$ for all $m\geq 2$ and $1 \leq r < m$. These PDS correspond to maximal arcs in the Desarguesian projective planes PG$(2, 2^m)$. Davis et al. \cite{DHJP2024} and also De Winter \cite{dewinter23} presented constructions of PDS with Denniston parameters $(p^{3m}, (p^{m+r}-p^m+p^r)(p^m-1), p^m-p^r+(p^{m+r}-p^m+p^r)(p^r-2), (p^{m+r}-p^m+p^r)(p^r-1))$ in elementary abelian groups of order $p^{3m}$ for all $m \geq 2$ and $r \in \{1, m-1\}$, where $p$ is an odd prime. The constructions in \cite{DHJP2024, dewinter23} are particularly intriguing, as it was shown by Ball, Blokhuis, and Mazzocca \cite{BBM1997} that no nontrivial maximal arcs in PG$(2, q^m)$ exist for any odd prime power $q$. In this paper, we show that PDS with Denniston parameters $(q^{3m}, (q^{m+r}-q^m+q^r)(q^m-1), q^m-q^r+(q^{m+r}-q^m+q^r)(q^r-2), (q^{m+r}-q^m+q^r)(q^r-1))$ exist in elementary abelian groups of order $q^{3m}$ for all $m \geq 2$ and $1 \leq r < m$, where $q$ is an arbitrary prime power.


[211] 2407.15637

The reciprocal complement of a polynomial ring in several variables over a field

The *reciprocal complement* $R(D)$ of an integral domain $D$ is the subring of its fraction field generated by the reciprocals of its nonzero elements. Many properties of $R(D)$ are determined when $D$ is a polynomial ring in $n\geq 2$ variables over a field. In particular, $R(D)$ is an $n$-dimensional, local, non-Noetherian, non-integrally closed, non-factorial, atomic G-domain, with infinitely many prime ideals at each height other than $0$ and $n$.


[212] 2407.15638

Orderings of the finite mixture with modified proportional hazard rate model

In this paper, we consider finite mixture models with modified proportional hazard rates. Sufficient conditions for the usual stochastic order and the hazard order are established under chain majorization. We study stochastic comparisons under different settings of T-transform for various values of chain majorization. We establish usual stochastic order and hazard rate order between two mixture random variables when a matrix of model parameters and mixing proportions changes to another matrix in some mathematical sense. Sufficient conditions for the star order and Lorenz order are established under weakly supermajorization. The results of this paper are illustrated with numerical examples.


[213] 2407.15644

Elliptic curves and spin

In the early 2000s, Ramakrishna asked the question: For the elliptic curve $$ E: y^2 = x^3 - x, $$ what is the density of primes $p$ for which the Fourier coefficient $a_p(E)$ is a cube modulo $p$? As a generalization of this question, Weston--Zaurova formulated conjectures concerning the distribution of power residues of degree $m$ of the Fourier coefficients of elliptic curves $E/\mathbb{Q}$ with complex multiplication. In this paper, we prove their conjecture for cubic residues using the analytic theory of spin. Our proof works for all elliptic curves $E$ with complex multiplication.


[214] 2407.15650

Sharp commutator estimates of all order for Coulomb and Riesz modulated energies

We prove functional inequalities in any dimension controlling the iterated derivatives along a transport of the Coulomb or super-Coulomb Riesz modulated energy in terms of the modulated energy itself. This modulated energy was introduced by the second author and collaborators in the study of mean-field limits and statistical mechanics of Coulomb/Riesz gases, where control of such derivatives by the energy itself is an essential ingredient. In this paper, we extend and improve such functional inequalities, proving estimates which are now sharp in their additive error term, in their density dependence, valid at arbitrary order of differentiation, and localizable to the support of the transport. Our method relies on the observation that these iterated derivatives are the quadratic form of a commutator. Taking advantage of the Riesz nature of the interaction, we identify these commutators as solutions to a degenerate elliptic equation with a right-hand side exhibiting a recursive structure in terms of lower-order commutators and develop a local regularity theory for the commutators, which may be of independent interest. These estimates have applications to obtaining sharp rates of convergence for mean-field limits, quasi-neutral limits, and in proving central limit theorems for the fluctuations of Coulomb/Riesz gases. In particular, we show here the expected $N^{\frac{\mathsf{s}}{\mathsf{d}}-1}$-rate in the modulated energy distance for the mean-field convergence of first-order Hamiltonian and gradient flows.


[215] 2407.15654

$K$-Positivity Preservers and their Generators

We study $K$-positivity preservers with given closed $K\subseteq\mathbb{R}^n$, i.e., linear maps $T:\mathbb{R}[x_1,\dots,x_n]\to\mathbb{R}[x_1,\dots,x_n]$ such that $T\mathrm{Pos}(K)\subseteq\mathrm{Pos}(K)$ holds, and their generators $A:\mathrm{R}[x_1,\dots,x_n]\to\mathrm{R}[x_1,\dots,x_n]$, i.e., $e^{tA}\mathrm{Pos}(K)\subseteq\mathrm{Pos}(K)$ holds for all $t\geq 0$. We characterize these maps $T$ for any closed $K\subseteq\mathbb{R}^n$ in Theorem 4.5. We characterize the maps $A$ in Theorem 5.8 for $K=\mathbb{R}^n$ and give partial results for general $K$. In Example 5.10 we give a map $A$ such that $e^{tA}$ is a positivity preserver for all $t\geq \tau$ for some $\tau>0$ but not for $t\in (0,\tau)$, i.e., we have an eventually positive semi-group.


[216] 2407.15664

Some new properties of the beta function and Ramanujan R-function

In this paper, the power series and hypergeometric series representations of the beta and Ramanujan functions \begin{equation*} \mathcal{B}\left( x\right) =\frac{\Gamma \left( x\right)^{2}}{\Gamma \left( 2x\right) }\text{ and }\mathcal{R}\left( x\right) =-2\psi \left( x\right) -2\gamma \end{equation*} are presented, which yield higher order monotonicity results related to $ \mathcal{B}(x)$ and $\mathcal{R}(x)$; the decreasing property of the functions $\mathcal{R}\left( x\right) /\mathcal{B}\left( x\right) $ and $[ \mathcal{B}(x) -\mathcal{R}(x)] /x^{2}$ on $\left( 0,\infty \right)$ are proved. Moreover, a conjecture put forward by Qiu et al. in [17] is proved to be true. As applications, several inequalities and identities are deduced. These results obtained in this paper may be helpful for the study of certain special functions. Finally, an interesting infinite series similar to Riemann zeta functions is observed initially.


[217] 2407.15667

On even $K$-groups over $p$-adic Lie extensions of global function fields

Let $p$ be a fixed prime number, and $F$ a global function field of characteristic not equal to $p$. In this paper, we shall study the growth of the Sylow $p$-subgroups of the even $K$-groups in a $p$-adic Lie extension of $F$, where the $p$-adic Lie extension is assumed to contain the cyclotomic $\mathbb{Z}_p$-extension of $F$. We also establish a duality between the direct limit and inverse limit of the even $K$-groups.


[218] 2407.15669

Delta-shock for the pressureless Euler-Poisson system

We study singularity formation for the pressureless Euler-Poisson system of cold ion dynamics. In contrast to the Euler-Poisson system with pressure, when its smooth solutions experience $C^1$ blow-up, the $L^\infty$ norm of the density becomes unbounded, which is often referred to as a delta-shock. We provide a constructive proof of singularity formation to obtain an exact blow-up profile and the detailed asymptotic behavior of the solutions near the blow-up point in both time and space. Our result indicates that at the blow-up time $t=T_\ast$, the density function is unbounded but is locally integrable with the profile of $\rho(x,T_\ast) \sim (x-x_*)^{-2/3}$ near the blow-up point $x=x_\ast$. This profile is not yet a Dirac measure. On the other hand, the velocity function has $C^{1/3}$ regularity at the blow-up point. Loosely following our analysis, we also obtain an exact blow-up profile for the pressureless Euler equations.


[219] 2407.15679

Lattice subsequences of fixed points of Toeplitz substitutions

We define the modulo-$m$ Toeplitz fixed point generated by Toeplitz substitution and study the lattice subsequence of such fixed point. Moreover, we provide a method to check whether one modulo-$m$ Toeplitz fixed point is a lattice subsequence of another.


[220] 2407.15681

Inverse random potential scattering for the polyharmonic wave equation using far-field patterns

This paper addresses the inverse scattering problem of a random potential associated with the polyharmonic wave equation in two and three dimensions. The random potential is represented as a centered complex-valued generalized microlocally isotropic Gaussian random field, where its covariance and relation operators are characterized as conventional pseudo-differential operators. Regarding the direct scattering problem, the well-posedness is established in the distribution sense for sufficiently large wavenumbers through analysis of the corresponding Lippmann--Schwinger integral equation. Furthermore, in the context of the inverse scattering problem, the uniqueness is attained in recovering the microlocal strengths of both the covariance and relation operators of the random potential. Notably, this is accomplished with only a single realization of the backscattering far-field patterns averaged over the high-frequency band.


[221] 2407.15684

A refinement of the Šidák-Khatri inequality and a strong Gaussian correlation conjecture

We prove two special cases of a strengthened Gaussian correlation conjecture, due to Tehranchi, and show that if the conjecture holds asymptotically, it holds for any dimension. Additionally, we use these special cases to prove a refined version of the \v{S}id\'ak-Khatri inequality.


[222] 2407.15693

Fisher-Rao Gradient Flow: Geodesic Convexity and Functional Inequalities

The dynamics of probability density functions has been extensively studied in science and engineering to understand physical phenomena and facilitate algorithmic design. Of particular interest are dynamics that can be formulated as gradient flows of energy functionals under the Wasserstein metric. The development of functional inequalities, such as the log-Sobolev inequality, plays a pivotal role in analyzing the convergence of these dynamics. The goal of this paper is to parallel the success of techniques using functional inequalities, for dynamics that are gradient flows under the Fisher-Rao metric, with various $f$-divergences as energy functionals. Such dynamics take the form of a nonlocal differential equation, for which existing analysis critically relies on using the explicit solution formula in special cases. We provide a comprehensive study on functional inequalities and the relevant geodesic convexity for Fisher-Rao gradient flows under minimal assumptions. A notable feature of the obtained functional inequalities is that they do not depend on the log-concavity or log-Sobolev constants of the target distribution. Consequently, the convergence rate of the dynamics (assuming well-posed) is uniform across general target distributions, making them potentially desirable dynamics for posterior sampling applications in Bayesian inference.


[223] 2407.15696

History of confluent Vandermonde matrices and inverting them algorithms

The author was encouraged to write this review by numerous enquiries from researchers all over the world, who needed a ready-to-use algorithm for the inversion of confluent Vandermonde matrices which works in quadratic time for any values of the parameters allowed by the definition, including the case of large root multiplicities of the characteristic polynomial. Article gives the history of the title special matrix since 1891 and surveys algorithms for solving linear systems with the title class matrix and inverting it. In particular, it presents, also by example, a numerical algorithm which does not use symbolic computations and is ready to be implemented in a general-purpose programming language or in a specific mathematical package.


[224] 2407.15699

Parahoric Hecke Ext-algebras in characteristic $p$

Let $\mathfrak{F}$ be a nonarchimedean local field of residual characteristic $p$, and let $G$ denote the group of $\mathfrak{F}$-points of a connected reductive group over $\mathfrak{F}$. For an open compact subgroup $\mathcal{U}$ of $G$ and a unital commutative ring $k$, we let $\mathbf{X}_{\mathcal{U}}$ denote the space of compactly supported $k$-valued functions on $G/\mathcal{U}$. Building on work of Ollivier--Schneider, we investigate the graded $\textrm{Ext}$-algebra $E_{\mathcal{U}}^* := \textrm{Ext}_G^*(\mathbf{X}_{\mathcal{U}},\mathbf{X}_{\mathcal{U}})^{\textrm{op}}$. In particular, we describe the Yoneda product, an involutive anti-automorphism, and (when $k$ is a field of characteristic $p$ and $\mathcal{U}$ has no $p$-torsion) a duality operation. We allow for the reductive group to be non-split, and for the open compact subgroup $\mathcal{U}$ to be non-pro-$p$. Specializing further to the case $G = \textrm{SL}_2(\mathbb{Q}_p)$ with $p \geq 5$ and a coefficient field of characteristic $p$, we obtain more precise results when $\mathcal{U}$ is equal to an Iwahori subgroup $J$ or a hyperspecial maximal compact subgroup $K$. In particular, we compute the structure of $E_J^*$ as an $E_J^0$-bimodule, obtain an explicit description of the center $\mathcal{Z}(E_J^*)$ of $E_J^*$, and construct a surjective morphism of algebras $\mathcal{Z}(E_J^*) \longrightarrow E_K^*$ (analogous to the compatibility between Bernstein and Satake isomorphisms in characteristic 0). From this we deduce the (somewhat surprising) fact that $E_K^*$ is not graded-commutative, contrary to what happens for almost all $\ell$-modular characteristics.


[225] 2407.15704

Distributions of consecutive level spacings of GUE and their ratio: an analytic treatment

In recent studies of many-body localization in nonintegrable quantum systems, the distribution of the ratio of two consecutive energy level spacings, $r_n=(E_{n+1}-E_n)/(E_{n}-E_{n-1})$ or $\tilde{r}_n=\min(r_n,r_n^{-1})$, has been used as a measure to quantify the chaoticity, alternative to the more conventional distribution of the level spacings, $s_n=\bar{\rho}(E_n)(E_{n+1}-E_n)$, as the former unnecessitates the unfolding required for the latter. Based on our previous work on the Tracy-Widom approach to the Janossy densities, we present analytic expressions for the joint probability distribution of two consecutive eigenvalue spacings and the distribution of their ratio for the Gaussian unitary ensemble (GUE) of random Hermitian $N\times N$ matrices at $N\to\infty$, in terms of a system of differential equations. As a showcase of efficacy of our results for characterizing an approach to quantum chaoticity, we contrast them to arguably the most ideal of all quantum-chaotic spectra: the zeroes of the Riemann $\zeta$ function on the critical line at increasing heights.


[226] 2407.15713

Inverse problems for coupled nonlocal nonlinear systems arising in mathematical biology

In this paper, we propose and study several inverse problems of determining unknown parameters in nonlocal nonlinear coupled PDE systems, including the potentials, nonlinear interaction functions and time-fractional orders. In these coupled systems, we enforce non-negativity of the solutions, aligning with realistic scenarios in biology and ecology. There are several salient features of our inverse problem study: the drastic reduction in measurement/observation data due to averaging effects, the nonlinear coupling between multiple equations, and the nonlocality arising from fractional-type derivatives. These factors present significant challenges to our inverse problem, and such inverse problems have never been explored in previous literature. To address these challenges, we develop new and effective schemes. Our approach involves properly controlling the injection of different source terms to obtain multiple sets of mean flux data. This allows us to achieve unique identifiability results and accurately determine the unknown parameters. Finally, we establish a connection between our study and practical applications in biology, further highlighting the relevance of our work in real-world contexts.


[227] 2407.15726

Classical discrete operators on variable $\ell^{p(\cdot)}(\mathbb{Z})$ spaces

We show, by applying discrete weighted norm inequalities and the Rubio de Francia algorithm, that the discrete Hilbert transform and discrete Riesz potential are bounded on variable $\ell^{p(\cdot)}(\mathbb{Z})$ spaces whenever the discrete Hardy-Littlewood maximal is bounded on $\ell^{p(\cdot)}(\mathbb{Z})$. We also obtain vector-valued inequalities for the discrete fractional maximal operator.


[228] 2407.15729

Self-Sustainable Metasurface-Assisted mmWave Indoor Communication System

In the design of a metasurface-assisted system for indoor environments, it is essential to take into account not only the performance gains and coverage extension provided by the metasurface but also the operating costs brought by its reconfigurability, such as powering and cabling. These costs can present challenges, particularly in indoor dense spaces (IDSs). A self-sustainable metasurface (SSM), which retains reconfigurability unlike a static metasurface (SMS), achieves a lower operating cost than a reconfigurable intelligent surface (RIS) by being self-sustainable through power harvesting. In this paper, in order to find a better trade-off between metasurface gain, coverage, and operating cost, the design and performance of an SSM-assisted indoor mmWave communication system are investigated. We first simplify the design of the SSM-assisted system by considering the use of SSMs in a preset-based manner and the formation of coverage groups by associating SSMs with the closest user equipments (UEs). We propose a two-stage iterative algorithm to maximize the minimum data rate in the system by jointly deciding the association between the UEs and the SSMs, the phase-shifts of the SSMs, and allocating time resources for each UE. The non-convexities that exist in the proposed optimization problem are tackled using the feasible point pursuit successive convex approximation method and the concave-convex procedure. To understand the best scenario for using SSM, the resulting performance is compared with that achieved with RIS and SMS. Our numerical results indicate that SSMs are best utilized in a small environment where self-sustainability is easier to achieve when the budget for operating costs is tight.


[229] 2407.15742

New solutions for the Lane-Emden problem on planar domains

We consider the Lane-Emden problem on planar domains. When the exponent is large, the existence and multiplicity of solutions strongly depend on the geometric properties of the domain, which also deeply affect their qualitative behavior. Remarkably, a wide variety of solutions, both positive and sign-changing, have been found when the exponent is sufficiently large. In this paper, we focus on this topic and fine new sign-changing solutions that exhibit an unexpected concentration phenomenon as the exponent approaches infinity.


[230] 2407.15743

Enhanced Achievable DoF Bounds for Cache-Aided MIMO Communication Systems

Integrating coded caching (CC) into multiple-input multiple-output (MIMO) communications may significantly enhance the achievable degrees of freedom (DoF) of the wireless networks. In this paper, we consider a cache-aided MIMO configuration with a CC gain $t$, where a server with $L$ Tx antennas communicates with $K$ users, each with $G$ Rx antennas. In the proposed content-aware MIMO strategy, we carefully adjust the number of users $\Omega$ and the number of parallel streams decoded by each user $\beta$ served in each transmission to maximize the DoF. As a result, we achieve a DoF of ${\max_{\beta, \Omega }}{\Omega \beta}$, where ${\beta \le \mathrm{min}\big(G,\frac{L \binom{\Omega-1}{t}}{1 + (\Omega - t-1)\binom{\Omega-1}{t}}\big)}$. To prove the achievability of the proposed DoF bound, we provide a novel transmission strategy based on the simultaneous unicasting of multiple data streams. In this strategy, the missing data packets are scheduled such that the number of parallel streams per transmission is maximized while the decodability of all useful terms by each target user is guaranteed. Numerical simulations validate the findings, confirming the enhanced DoF and improved performance of the proposed design.


[231] 2407.15744

A matrix algebra for graphical statistical models

Directed mixed graphs permit directed and bidirected edges between any two vertices. They were first considered in the path analysis developed by Sewall Wright and play an essential role in statistical modeling. We introduce a matrix algebra for walks on such graphs. Each element of the algebra is a matrix whose entries are sets of walks on the graph from the corresponding row to the corresponding column. The matrix algebra is then generated by applying addition (set union), multiplication (concatenation), and transpose to the two basic matrices consisting of directed and bidirected edges. We use it to formalize, in the context of Gaussian linear systems, the correspondence between important graphical concepts such as latent projection and graph separation with important probabilistic concepts such as marginalization and (conditional) independence. In two further examples regarding confounder adjustment and the augmentation criterion, we illustrate how the algebra allows us to visualize complex graphical proofs. A "dictionary" and LATEX macros for the matrix algebra are provided in the Appendix.


[232] 2407.15746

Complementation in Continuous Cohomology with Coefficients in Banach Modules

In this article, we introduce the concept of weakly uniquely stationary representations. This framework enables us to investigate the complementability of closed subspaces within the context of continuous cohomology with coeffcients in Banach modules. As an application, we extend and refine several cohomological results from the literature, particularly in the settings of nilpotent groups, products of groups, and lattices.


[233] 2407.15751

Rapidly convergent series expansions for a class of resolvents

Following advances in the theory of composites we develop rapidly converging series expansions about $z=\infty$ for the resolvent ${\bf R}(z)=[z{\bf I}-{\bf P}^\dagger{\bf Q}{\bf P}]^{-1}$ where ${\bf Q}$ is an orthogonal projection and ${\bf P}$ is such that ${\bf P}^\dagger{\bf P}={\bf I}$ while $\Gamma={\bf P}{\bf P}^\dagger$ is an orthogonal projection. It is assumed that the spectrum of ${\bf P}^\dagger{\bf Q}{\bf P}$ lies within the interval $[z^-,z^+]$ for some known $z^+\leq 1$ and $z^-\geq 0$. The series converges in the entire $z$-plane excluding the cut $[z^-,z^+]$. It is obtained using subspace substitution, where the desired resolvent is tied to a resolvent in a larger space and ${\bf Q}$ gets replaced by a projection $\underline{\bf Q}$ that is no longer orthogonal.


[234] 2407.15752

Broad and Spectral-Efficient Beamforming for the Uni-polarized Reconfigurable Intelligent Surfaces

A reconfigurable intelligent surface (RIS) is composed of low-cost elements that manipulate the propagation environment from a transmitter by intelligently applying phase shifts to incoming signals before they are reflected. This paper explores a uni-polarized RIS with linear shape aimed at transmitting a common signal to multiple user equipments (UEs) spread across a wide angular region. To achieve uniform coverage, the uni-polarized RIS is designed to emit a broad and spectral-efficient beam featuring a spatially flat-like array factor, diverging from the conventional narrow beam approach. To achieve this objective, we start by deriving probabilistic lower and upper bounds for the average spectral efficiency (SE) delivered to the UEs. Leveraging the insights from the lower bound, we focus on optimizing the minimum value of the power domain array factor (PDAF) across a range of azimuth angles from \(-\frac{\pi}{2}\) to \(\frac{\pi}{2}\). We employ the continuous genetic algorithm (CGA) for this optimization task, aiming to improve the SE delivered to the UEs while also creating a wide beam. Extensive simulation experiments are carried out to assess the performance of the proposed code, focusing on key metrics such as the minimum and average values of the PDAF and the SE delivered to the UEs. Our findings demonstrate that the proposed code enhances the minimum SE delivered to the UEs while maintaining the desired attribute of a broad beam. This performance is notably superior to that of established codes, including the Barker, Frank, and Chu codes.


[235] 2407.15760

A Hamilton-Jacobi approach to road-field reaction-diffusion models

We consider the road-field reaction-diffusion model introduced by Berestycki, Roquejoffre, and Rossi. By performing a "thin-front limit," we are able to deduce a Hamilton-Jacobi equation with a suitable effective Hamiltonian on the road that governs the front location of the road-field model. Our main motivation is to apply the theory of strong (flux-limited) viscosity solutions in order to determine a control formulation interpretation of the front location. In view of the ecological meaning of the road-field model, this is natural as it casts the invasion problem as one of finding optimal paths that balance the positive growth rate in the field with the fast diffusion on the road. Our main contribution is a nearly complete picture of the behavior on two-road conical domains. When the diffusivities on each road are the same, we show that the propagation speed in each direction in the cone can be computed via those associated with one-road half-space problem. When the diffusivities differ, we show that the speed along the faster road is unchanged, while the speed along the slower road can be enhanced. Along the way we provide a new proof of known results on the one-road half-space problem via our approach.


[236] 2407.15764

Huber means on Riemannian manifolds

This article introduces Huber means on Riemannian manifolds, providing a robust alternative to the Frechet mean by integrating elements of both square and absolute loss functions. The Huber means are designed to be highly resistant to outliers while maintaining efficiency, making it a valuable generalization of Huber's M-estimator for manifold-valued data. We comprehensively investigate the statistical and computational aspects of Huber means, demonstrating their utility in manifold-valued data analysis. Specifically, we establish minimal conditions for ensuring the existence and uniqueness of the Huber mean and discuss regularity conditions for unbiasedness. The Huber means are statistically consistent and enjoy the central limit theorem. Additionally, we propose a moment-based estimator for the limiting covariance matrix, which is used to construct a robust one-sample location test procedure and an approximate confidence region for location parameters. Huber means are shown to be highly robust and efficient in the presence of outliers or under heavy-tailed distribution. To be more specific, it achieves a breakdown point of at least 0.5, the highest among all isometric equivariant estimators, and is more efficient than the Frechet mean under heavy-tailed distribution. Numerical examples on spheres and the set of symmetric positive-definite matrices further illustrate the efficiency and reliability of the proposed Huber means on Riemannian manifolds.


[237] 2407.15765

Skolem, Gödel, and Hilbert fibrations

Grothendieck fibrations are fundamental in capturing the concept of dependency, notably in categorical semantics of type theory and programming languages. A relevant instance are Dialectica fibrations which generalise G\"odel's Dialectica proof interpretation and have been widely studied in recent years. We characterise when a given fibration is a generalised, dependent Dialectica fibration, namely an iterated completion of a fibration by dependent products and sums (along a given class of display maps). From a technical perspective, we complement the work of Hofstra on Dialectica fibrations by an internal viewpoint, categorifying the classical notion of quantifier-freeness. We also generalise both Hofstra's and Trotta et al.'s work on G\"odel fibrations to the dependent case, replacing the class of cartesian projections in the base category by arbitrary display maps. We discuss how this recovers a range of relevant examples in categorical logic and proof theory. Moreover, as another instance, we introduce Hilbert fibrations, providing a categorical understanding of Hilbert's $\epsilon$- and $\tau$-operators well-known from proof theory.


[238] 2407.15769

Hopf algebras and associative representations of two-dimensional evolution algebras

In this paper, we establish a connection between evolution algebras of dimension two and Hopf algebras, via the algebraic group of automorphisms of an evolution algebra. Initially, we describe the Hopf algebra associated with the automorphism group of a 2-dimensional evolution algebra. Subsequently, for a 2-dimensional evolution algebra $A$ over a field $K$, we detail the relation between the algebra associated with the (tight) universal associative and commutative representation of $A$, referred to as the (tight) $p$-algebra, and the corresponding Hopf algebra, $\mathcal{H}$, representing the affine group scheme $\text{Aut}(A)$. Our analysis involves the computation of the (tight) $p-$algebra associated with any 2-dimensional evolution algebra, whenever it exists. We find that $\text{Aut}(A)=1$ if and only if there is no faithful associative and commutative representation for $A$. Moreover, there is a faithful associative and commutative representation for $A$ if and only if $\mathcal{H}\not\cong K$ and $\text{char} (K)\neq 2$, or $\mathcal{H}\not\cong K(\epsilon)$ (the dual numbers algebra) and $\mathcal{H}\not\cong K$ in case of $\text{char} (K)= 2$. Furthermore, if $A$ is perfect and has a faithful tight $p$-algebra, then this $p$-algebra is isomorphic to $\mathcal{H}$ (as algebras). Finally, we derive implications for arbitrary finite-dimensional evolution algebras.


[239] 2407.15772

Siegel $\mathfrak{p}^2$ Vectors for Representations of $GSp(4)$

Let $F$ be a $p$-adic field and $(\pi, V)$ an irreducible complex representation of $G=GSp(4, F)$ with trivial central character. Let ${\rm Si}(\mathfrak{p}^2)\subset G$ denote the Siegel congruence subgroup of level $\mathfrak{p}^2$ and $u\in N_G({\rm Si}(\mathfrak{p}^2))$ the Atkin-Lehner element. We compute the dimension of the space of ${\rm Si}(\mathfrak{p}^2)$-fixed vectors in $V$ as well as the signatures of the involutions $\pi(u)$ acting on these spaces.


[240] 2407.15774

Metric mean dimension, Hölder regularity and Assouad spectrum

Metric mean dimension is a geometric invariant of dynamical systems with infinite topological entropy. We relate this concept with the fractal structure of the phase space and the H\"older regularity of the map. Afterwards we improve our general estimates in a family of interval maps by computing the metric mean dimension in a way similar to the Misiurewicz formula for the entropy, which in particular shows that our bounds are sharp. Of independent interest, we develop a dynamical analogue of Minkowski-Bouligand dimension for maps acting on Ahlfors regular spaces.


[241] 2407.15775

Improving Greedy Algorithms for Rational Approximation

When developing robust preconditioners for multiphysics problems, fractional functions of the Laplace operator often arise and need to be inverted. Rational approximation in the uniform norm can be used to convert inverting those fractional operators into inverting a series of shifted Laplace operators. Care must be taken in the approximation so that the shifted Laplace operators remain symmetric positive definite, making them better conditioned. In this work, we study two greedy algorithms for finding rational approximations to such fractional operators. The first algorithm improves the orthogonal greedy algorithm discussed in [Li et al., SISC, 2024] by adding one minimization step in the uniform norm to the procedure. The second approach employs the weak Chebyshev greedy algorithm in the uniform norm. Both methods yield non-increasing error. Numerical results confirm the effectiveness of our proposed algorithms, which are also flexible and applicable to other approximation problems. Moreover, with effective rational approximations to the fractional operator, the resulting algorithms show good performance in preconditioning a Darcy-Stokes coupled problem.


[242] 2407.15782

Reconfigurable Intelligent Surface Empowered Full Duplex Systems: Opportunities and Challenges

Reconfigurable intelligent surfaces (RISs) have emerged as a promising technology in wireless communications. Simultaneously transmitting and reflecting RIS (STAR-RISs) in particular have garnered significant attention due to their dual capabilities of simultaneous transmission and reflection, underscoring their potential applications in critical scenarios within the forthcoming sixth-generation (6G) technology landscape. Moreover, full-duplex (FD) systems have emerged as a breakthrough research direction in wireless transmission technology due to their high spectral efficiency. This paper explores the application potential of STAR-RIS in FD systems for future wireless communications, presenting an innovative technology that provides robust self-interference cancellation (SIC) capabilities for FD systems. We utilize the refraction functionality of STAR-RIS enhances the transmission capacity of FD systems, while its reflection functionality is used to eliminate self interference within the FD system. We delve into the applications of two different types of STAR-RIS in FD systems and compare their performance through simulations. Furthermore, we discuss the performance differences of STAR-RIS empowered FD systems under various configurations in a case study, and demonstrate the superiority of the proposed deep learning-based optimization algorithm. Finally, we discuss possible future research directions for STAR-RIS empowered FD systems.


[243] 2407.15785

Weak Freiman isomorphisms and sequencings of small sets

In this paper, we introduce a weakening of the Freiman isomorphisms between subsets of non necessarily abelian groups. Inspired by the breakthrough result of Kravitz, [14], on cyclic groups, as an application, we prove that any subset of size $k$ of the dihedral group $D_{2m}$ (and, more in general, of a class of semidirect products) is sequenceable, provided that the prime factors of $m$ are larger than $k!$. Also, a refined bound of $k!/2$ for the size of the prime factors of $m$ can be obtained for cyclic groups $\mathbb{Z}_m$, slightly improving the result of [14].


[244] 2407.15803

Unbounded operators and the uncertainty principle

We study a variant of the uncertainty principle in terms of the annihilation and creation operator on generalized Segal Bargmann spaces, which are used for the FBI-Bargmann transform. In addition, we compute the Berezin transform of these operators and indicate how to use spaces of entire functions in one variable to study the Szeg\H{o} kernel for hypersurfaces in $\mathbb C^2.$


[245] 2407.15812

On the stability of blowup solutions to the complex Ginzburg-Landau equation in R^d

Building upon the idea in \cite{HNWarXiv24}, we establish stability of the type-I blowup with log correction for the complex Ginzburg-Landau equation. In the amplitude-phase representation, a generalized dynamic rescaling formulation is introduced, with modulation parameters capturing the spatial translation and rotation symmetries of the equation and novel additional modulation parameters perturbing the scaling symmetry. This new formulation provides enough degrees of freedom to impose normalization conditions on the rescaled solution, completely eliminating the unstable and neutrally stable modes of the linearized operator around the blowup profile. It enables us to establish the full stability of the blowup by enforcing vanishing conditions via the choice of normalization and using weighted energy estimates, without relying on a topological argument or a spectrum analysis. The log correction for the blowup rate is captured by the energy estimates and refined estimates of the modulation parameters.


[246] 2407.15818

The connectivity of Vietoris-Rips complexes of spheres

We survey what is known and unknown about Vietoris-Rips complexes and thickenings of spheres. Afterwards, we show how to control the homotopy connectivity of Vietoris-Rips complexes of spheres in terms of coverings of spheres and projective spaces. Let $S^n$ be the $n$-sphere with the geodesic metric, and of diameter $\pi$, and let $\delta > 0$. Suppose that the first nontrivial homotopy group of the Vietoris-Rips complex $\mathrm{VR}(S^n;\pi-\delta)$ of the $n$-sphere at scale $\pi-\delta$ occurs in dimension $k$, i.e., suppose that the connectivity is $k-1$. Then $\mathrm{cov}_{S^n}(2k+2) \le \delta < 2\cdot \mathrm{cov}_{\mathbb{R}P^n}(k)$. In other words, there exist $2k+2$ balls of radius $\delta$ that cover $S^n$, and no set of $k$ balls of radius $\frac{\delta}{2}$ cover the projective space $\mathbb{R}P^n$. As a corollary, the homotopy type of $\mathrm{VR}(S^n;r)$ changes infinitely many times as the scale $r$ increases.


[247] 2407.15825

High-order numerical method for solving elliptic partial differential equations on unfitted node sets

In this paper, we present how high-order accurate solutions to elliptic partial differential equations can be achieved in arbitrary spatial domains using radial basis function-generated finite differences (RBF-FD) on unfitted node sets (i.e., not adjusted to the domain boundary). In this novel method, we only collocate on nodes interior to the domain boundary and enforce boundary conditions as constraints by means of Lagrange multipliers. This combination enables full geometric flexibility near boundaries without compromising the high-order accuracy of the RBF-FD method. The high-order accuracy and robustness of two formulations of this approach are illustrated by numerical experiments.


[248] 2407.15827

Kadec-type theorems for sampled group orbits

We extend the classical Kadec 1/4 theorem for systems of exponential functions on an interval to frames and atomic decompositions formed by sampling an orbit of a vector under an isometric group representation.


[249] 2407.15830

Stabilized measures of association: comparing the birational geometry of varieties having different dimensions

We extend some of the measures of association defined by Lazarsfeld and Martin, obtaining useful invariants to compare the birational geometry of two varieties having different dimensions. We explore such invariants providing examples and computing them in some cases of interest.


[250] 2407.11092

Configuration Spaces of Points: A User's Guide

This extensive survey is an invited contribution to the Encyclopedia of Mathematical Physics, 2nd edition. It covers both classical and more modern aspects of configuration spaces of points on a "ground space" $M$. Most results pertain to $M$ a manifold. Configuration spaces of points have become so omnipresent in so many areas of mathematics, physics, and even the applied sciences, that a survey can only cover a selection of topics. We review key ideas, constructions, and results.


[251] 2407.14510

A Review on Response Strategies in Infrastructure Network Restoration

This paper reviews the literature on response strategies for restoring infrastructure networks in the aftermath of a disaster. Our motivation for this review is twofold. First, the frequency and magnitude of natural and man-made disasters (e.g., wild fires, tornadoes, global pandemics, terrorist attacks) have been increasing. These events disrupt the operation of infrastructure networks, preventing the delivery of vital goods and services such as power and food. Therefore, it is critical to understand the state-of-the-art in responding to network disruptions in order to develop efficient strategies to mitigate their impacts. Second, it is critical to enable timely decisions in a rapidly changing and unpredictable environment while accounting for numerous interrelated factors. Because the vast majority of response decision problems are computationally challenging, quickly finding solutions that are compatible with real-time decision making is a difficult task. Hence, it is important to understand the nature of response activities and decisions, as well as the available solution methodologies and inherent trade-offs between computation time and solution quality. We review quantitative response methodologies developed for infrastructure network restoration, classifying relevant studies based on the properties of the underlying network. In particular, we focus on resource allocation, scheduling, routing and repair efforts within the domain of power, road, and water, oil and gas network restoration. We also discuss open research questions and future research directions.


[252] 2407.14518

Accurate Analysis of Sparse Random Projections

There has been recently a lot of research on sparse variants of random projections, faster adaptations of the state-of-the-art dimensionality reduction technique originally due to Johsnon and Lindenstrauss. Although the construction is very simple, its analyses are notoriously complicated. Meeting the demand for both simplicity and accuracy, this work establishes sharp sub-poissonian tail bounds for the distribution of sparse random projections. Compared to other works, this analysis provide superior numerical guarantees (exactly matching impossibility results) while being arguably less complicated (the technique resembles Bennet's Inequality and is of independent interest).


[253] 2407.14703

Generalizing and transporting causal inferences from randomized trials in the presence of trial engagement effects

Trial engagement effects are effects of trial participation on the outcome that are not mediated by treatment assignment. Most work on extending (generalizing or transporting) causal inferences from a randomized trial to a target population has, explicitly or implicitly, assumed that trial engagement effects are absent, allowing evidence about the effects of the treatments examined in trials to be applied to non-experimental settings. Here, we define novel causal estimands and present identification results for generalizability and transportability analyses in the presence of trial engagement effects. Our approach allows for trial engagement effects under assumptions of no causal interaction between trial participation and treatment assignment on the absolute or relative scales. We show that under these assumptions, even in the presence of trial engagement effects, the trial data can be combined with covariate data from the target population to identify average treatment effects in the context of usual care as implemented in the target population (i.e., outside the experimental setting). The identifying observed data functionals under these no-interaction assumptions are the same as those obtained under the stronger identifiability conditions that have been invoked in prior work. Therefore, our results suggest a new interpretation for previously proposed generalizability and transportability estimators; this interpretation may be useful in analyses under causal structures where background knowledge suggests that trial engagement effects are present but interactions between trial participation and treatment are negligible.


[254] 2407.14765

Data Augmentation in Graph Neural Networks: The Role of Generated Synthetic Graphs

Graphs are crucial for representing interrelated data and aiding predictive modeling by capturing complex relationships. Achieving high-quality graph representation is important for identifying linked patterns, leading to improvements in Graph Neural Networks (GNNs) to better capture data structures. However, challenges such as data scarcity, high collection costs, and ethical concerns limit progress. As a result, generative models and data augmentation have become more and more popular. This study explores using generated graphs for data augmentation, comparing the performance of combining generated graphs with real graphs, and examining the effect of different quantities of generated graphs on graph classification tasks. The experiments show that balancing scalability and quality requires different generators based on graph size. Our results introduce a new approach to graph data augmentation, ensuring consistent labels and enhancing classification performance.


[255] 2407.14858

An encryption algorithm using a generalization of the Markovski algorithm and a system of orthogonal operations based on T-quasigroups

Here is a more detailed description of the algorithm proposed in [1]. This algorithm simultaneously uses two cryptographic procedures: encryption using a generalization of the Markovski algorithm [2] and encryption using a system of orthogonal operations. In this paper, we present an implementation of this algorithm based on T-quasigroups, more precisely, based on medial quasigroups.


[256] 2407.14864

Computing Stokes flows in periodic channels via rational approximation

Rational approximation has proven to be a powerful method for solving two-dimensional (2D) fluid problems. At small Reynolds numbers, 2D Stokes flows can be represented by two analytic functions, known as Goursat functions. Xue, Waters and Trefethen [SIAM J. Sci. Comput., 46 (2024), pp. A1214-A1234] recently introduced the LARS algorithm (Lightning-AAA Rational Stokes) for computing 2D Stokes flows in general domains by approximating the Goursat functions using rational functions. In this paper, we introduce a new algorithm for computing 2D Stokes flows in periodic channels using trigonometric rational functions, with poles placed via the AAA-LS algorithm [Costa and Trefethen, European Congr. Math., 2023] in a conformal map of the domain boundary. We apply the algorithm to Poiseuille and Couette problems between various periodic channel geometries, where solutions are computed to at least 6-digit accuracy in less than 1 second. The applicability of the algorithm is highlighted in the computation of the dynamics of fluid particles in unsteady Couette flows.


[257] 2407.14882

Reduced Effectiveness of Kolmogorov-Arnold Networks on Functions with Noise

It has been observed that even a small amount of noise introduced into the dataset can significantly degrade the performance of KAN. In this brief note, we aim to quantitatively evaluate the performance when noise is added to the dataset. We propose an oversampling technique combined with denoising to alleviate the impact of noise. Specifically, we employ kernel filtering based on diffusion maps for pre-filtering the noisy data for training KAN network. Our experiments show that while adding i.i.d. noise with any fixed SNR, when we increase the amount of training data by a factor of $r$, the test-loss (RMSE) of KANs will exhibit a performance trend like $\text{test-loss} \sim \mathcal{O}(r^{-\frac{1}{2}})$ as $r\to +\infty$. We conclude that applying both oversampling and filtering strategies can reduce the detrimental effects of noise. Nevertheless, determining the optimal variance for the kernel filtering process is challenging, and enhancing the volume of training data substantially increases the associated costs, because the training dataset needs to be expanded multiple times in comparison to the initial clean data. As a result, the noise present in the data ultimately diminishes the effectiveness of Kolmogorov-Arnold networks.


[258] 2407.14986

Structured Input-Output Modeling and Robust Stability Analysis of Compressible Flows

The recently introduced structured input-output analysis is a powerful method for capturing nonlinear phenomena associated with incompressible flows, and this paper extends that method to the compressible regime. The proposed method relies upon a reformulation of the compressible Navier-Stokes equations, which allows for an exact quadratic formulation of the dynamics of perturbations about a steady base flow. To facilitate the structured input-output analysis, a pseudo-linear model for the quadratic nonlinearity is proposed and the structural information of the nonlinearity is embedded into a structured uncertainty comprising unknown `perturbations'. The structured singular value framework is employed to compute the input-output gain, which provides an estimate of the robust stability margin of the flow perturbations, as well as the forcing and response modes that are consistent with the nonlinearity structure. The analysis is then carried out on a plane, laminar compressible Couette flow over a range of Mach numbers. The structured input-output gains identify an instability mechanism, characterized by a spanwise elongated structure in the streamwise-spanwise wavenumber space at a subsonic Mach number, that evolves into an oblique structure at sonic and supersonic Mach numbers. In addition, the structured input-output forcing and response modes provide insight into the thermodynamic and momentum characteristics associated with a source of instability. Comparisons with a resolvent/unstructured analysis reveal discrepancies in the distribution of input-output gains over the wavenumber space as well as in the modal behavior of an instability, thus highlighting the strong correlation between the structural information of the nonlinearity and the underlying flow physics.


[259] 2407.15007

Is Behavior Cloning All You Need? Understanding Horizon in Imitation Learning

Imitation learning (IL) aims to mimic the behavior of an expert in a sequential decision making task by learning from demonstrations, and has been widely applied to robotics, autonomous driving, and autoregressive text generation. The simplest approach to IL, behavior cloning (BC), is thought to incur sample complexity with unfavorable quadratic dependence on the problem horizon, motivating a variety of different online algorithms that attain improved linear horizon dependence under stronger assumptions on the data and the learner's access to the expert. We revisit the apparent gap between offline and online IL from a learning-theoretic perspective, with a focus on general policy classes up to and including deep neural networks. Through a new analysis of behavior cloning with the logarithmic loss, we show that it is possible to achieve horizon-independent sample complexity in offline IL whenever (i) the range of the cumulative payoffs is controlled, and (ii) an appropriate notion of supervised learning complexity for the policy class is controlled. Specializing our results to deterministic, stationary policies, we show that the gap between offline and online IL is not fundamental: (i) it is possible to achieve linear dependence on horizon in offline IL under dense rewards (matching what was previously only known to be achievable in online IL); and (ii) without further assumptions on the policy class, online IL cannot improve over offline IL with the logarithmic loss, even in benign MDPs. We complement our theoretical results with experiments on standard RL tasks and autoregressive language generation to validate the practical relevance of our findings.


[260] 2407.15008

Control of Instability in a Vlasov-Poisson System Through an External Electric Field

Plasma instabilities are a major concern in plasma science, for applications ranging from particle accelerators to nuclear fusion reactors. In this work, we consider the possibility of controlling such instabilities by adding an external electric field to the Vlasov--Poisson equations. Our approach to determining the external electric field is based on conducting a linear analysis of the resulting equations. We show that it is possible to select external electric fields that completely suppress the plasma instabilities present in the system when the equilibrium distribution and the perturbation are known. In fact, the proposed strategy returns the plasma to its equilibrium with a rate that is faster than exponential in time. We further perform numerical simulations of the nonlinear two-stream and bump-on-tail instabilities to verify our theory and to compare the different strategies that we propose in this work.


[261] 2407.15069

Dipole Cosmology in $f(Q)$-gravity

Symmetric teleparallel $f(Q)$-gravity allows for the presence of a perfect fluid with a tilted velocity in the Kantowski-Sachs geometry. In this dipole model, we consider an ideal gas and we investigate the evolution of the physical parameters. The tilt parameter is constrained by the nonlinear function $f(Q)$ through the non-diagonal equations of the field equations. We find that the dynamics always reduce to the vacuum solutions of STEGR. This includes the Kasner universe, when no cosmological term is introduced by the $f(Q)$ function, and the isotropic de Sitter universe, where $f\left( Q\right) $ plays the role of the cosmological constant. In the extreme tilt limit, the universe is consistently anisotropic and accelerated. However, the final solution matches that of STEGR.


[262] 2407.15110

Practical multi-fidelity machine learning: fusion of deterministic and Bayesian models

Multi-fidelity machine learning methods address the accuracy-efficiency trade-off by integrating scarce, resource-intensive high-fidelity data with abundant but less accurate low-fidelity data. We propose a practical multi-fidelity strategy for problems spanning low- and high-dimensional domains, integrating a non-probabilistic regression model for the low-fidelity with a Bayesian model for the high-fidelity. The models are trained in a staggered scheme, where the low-fidelity model is transfer-learned to the high-fidelity data and a Bayesian model is trained for the residual. This three-model strategy -- deterministic low-fidelity, transfer learning, and Bayesian residual -- leads to a prediction that includes uncertainty quantification both for noisy and noiseless multi-fidelity data. The strategy is general and unifies the topic, highlighting the expressivity trade-off between the transfer-learning and Bayesian models (a complex transfer-learning model leads to a simpler Bayesian model, and vice versa). We propose modeling choices for two scenarios, and argue in favor of using a linear transfer-learning model that fuses 1) kernel ridge regression for low-fidelity with Gaussian processes for high-fidelity; or 2) deep neural network for low-fidelity with a Bayesian neural network for high-fidelity. We demonstrate the effectiveness and efficiency of the proposed strategies and contrast them with the state-of-the-art based on various numerical examples. The simplicity of these formulations makes them practical for a broad scope of future engineering applications.


[263] 2407.15196

Channel Shaping Using Beyond Diagonal Reconfigurable Intelligent Surface: Analysis, Optimization, and Enhanced Flexibility

This paper investigates the capability of a passive Reconfigurable Intelligent Surface (RIS) to redistribute the singular values of point-to-point Multiple-Input Multiple-Output (MIMO) channels for achieving power and rate gains. We depart from the conventional Diagonal (D)-RIS with diagonal phase shift matrix and adopt a Beyond Diagonal (BD) architecture that offers greater wave manipulation flexibility through element-wise connections. Specifically, we first provide shaping insights by characterizing the channel singular value regions attainable by D-RIS and BD-RIS via a novel geodesic optimization. Analytical singular value bounds are then derived to explore their shaping limits in typical deployment scenarios. As a side product, we tackle BD-RIS-aided MIMO rate maximization problem by a local-optimal Alternating Optimization (AO) and a shaping-inspired low-complexity approach. Results show that compared to D-RIS, BD-RIS significantly improves the dynamic range of all channel singular values, the trade-off in manipulating them, and thus the channel power and achievable rate. Those observations become more pronounced when the number of RIS elements and MIMO dimensions increase. Of particular interest, BD-RIS is shown to activate multi-stream transmission at lower transmit power than D-RIS, hence achieving the asymptotic Degrees of Freedom (DoF) at low Signal-to-Noise Ratio (SNR) thanks to its higher flexibility of shaping the distribution of channel singular values.


[264] 2407.15205

Backlund transformation of Kaup Kupershmidt equations with mutli soliton solutions un Darboux framewor

This article encloses the derivation of Darboux solutions for Kaup Kupershmidt equations with their generalization in determinantal form. One of the main focuses of this work is to construct the Backlund transformation for the different solutions of that equation through its associated Riccati equation and then that transformations further reduces to its algebraic analogue with the help of One-fold Darboux solution. Finally, its exact solutions upto three solitons are calculated with their graphical representations which reveal dynamical profiles of these solutions.


[265] 2407.15261

Pandora's Box Problem Over Time

The Pandora's Box problem models the search for the best alternative when evaluation is costly. In its simplest variant, a decision maker is presented with $n$ boxes, each associated with a cost of inspection and a distribution over the reward hidden within. The decision maker inspects a subset of these boxes one after the other, in a possibly adaptive ordering, and obtains as utility the difference between the largest reward uncovered and the sum of the inspection costs. While this classic version of the problem is well understood (Weitzman 1979), recent years have seen a flourishing of the literature on variants of the problem. In this paper, we introduce a general framework -- the Pandora's Box Over Time problem -- that captures a wide range of variants where time plays a role, e.g., as it might constrain the schedules of exploration and influence both costs and rewards. In the Pandora's Box Over Time problem, each box is characterized by time-dependent rewards and costs, and inspecting it might require a box-specific processing time. Moreover, once a box is inspected, its reward may deteriorate over time, possibly differently for each box. Our main result is an efficient $21.3$-approximation to the optimal strategy, which is NP-hard to compute in general. We further obtain improved results for the natural special cases where boxes have no processing time, or when costs and reward distributions do not depend on time (but rewards may deteriorate after inspecting).


[266] 2407.15276

Nonlinear Binscatter Methods

Binned scatter plots are a powerful statistical tool for empirical work in the social, behavioral, and biomedical sciences. Available methods rely on a quantile-based partitioning estimator of the conditional mean regression function to primarily construct flexible yet interpretable visualization methods, but they can also be used to estimate treatment effects, assess uncertainty, and test substantive domain-specific hypotheses. This paper introduces novel binscatter methods based on nonlinear, possibly nonsmooth M-estimation methods, covering generalized linear, robust, and quantile regression models. We provide a host of theoretical results and practical tools for local constant estimation along with piecewise polynomial and spline approximations, including (i) optimal tuning parameter (number of bins) selection, (ii) confidence bands, and (iii) formal statistical tests regarding functional form or shape restrictions. Our main results rely on novel strong approximations for general partitioning-based estimators covering random, data-driven partitions, which may be of independent interest. We demonstrate our methods with an empirical application studying the relation between the percentage of individuals without health insurance and per capita income at the zip-code level. We provide general-purpose software packages implementing our methods in Python, R, and Stata.


[267] 2407.15277

Conformal Predictions under Markovian Data

We study the split Conformal Prediction method when applied to Markovian data. We quantify the gap in terms of coverage induced by the correlations in the data (compared to exchangeable data). This gap strongly depends on the mixing properties of the underlying Markov chain, and we prove that it typically scales as $\sqrt{t_\mathrm{mix}\ln(n)/n}$ (where $t_\mathrm{mix}$ is the mixing time of the chain). We also derive upper bounds on the impact of the correlations on the size of the prediction set. Finally we present $K$-split CP, a method that consists in thinning the calibration dataset and that adapts to the mixing properties of the chain. Its coverage gap is reduced to $t_\mathrm{mix}/(n\ln(n))$ without really affecting the size of the prediction set. We finally test our algorithms on synthetic and real-world datasets.


[268] 2407.15289

Kerr--Newman Memory Effect

We bring the Kerr--Newman spacetime into the Bondi--Sachs gauge by means of zero angular momentum, null geodesics. We compute the memory effect produced at the black hole horizon by a transient gravitational shock wave, which from future null infinity is seen as a Bondi-Metzner-Sachs supertranslation. This results in a change of the supertransformation charges at infinity between the spacetime geometries defined by the black hole before, and after, the shockwave scattering. For an extremal Kerr--Newman black hole, we give the complementary description of this process in the near-horizon limit, as seen by an observer hovering over the horizon. In this limit, we compute the supertranformation charges and compare them to those calculated at null infinity. We analyze the effect of these transformations on the electromagnetic gauge field and explore the self-interaction between this and the angular momentum of the black hole.


[269] 2407.15301

U-learning for Prediction Inference via Combinatory Multi-Subsampling: With Applications to LASSO and Neural Networks

Epigenetic aging clocks play a pivotal role in estimating an individual's biological age through the examination of DNA methylation patterns at numerous CpG (Cytosine-phosphate-Guanine) sites within their genome. However, making valid inferences on predicted epigenetic ages, or more broadly, on predictions derived from high-dimensional inputs, presents challenges. We introduce a novel U-learning approach via combinatory multi-subsampling for making ensemble predictions and constructing confidence intervals for predictions of continuous outcomes when traditional asymptotic methods are not applicable. More specifically, our approach conceptualizes the ensemble estimators within the framework of generalized U-statistics and invokes the H\'ajek projection for deriving the variances of predictions and constructing confidence intervals with valid conditional coverage probabilities. We apply our approach to two commonly used predictive algorithms, Lasso and deep neural networks (DNNs), and illustrate the validity of inferences with extensive numerical studies. We have applied these methods to predict the DNA methylation age (DNAmAge) of patients with various health conditions, aiming to accurately characterize the aging process and potentially guide anti-aging interventions.


[270] 2407.15324

Cooperative Salvo Guidance over Leader-Follower Network with Free-Will Arbitrary Time Convergence

A cooperative salvo strategy is proposed in this paper which achieves consensus among the interceptors within a pre-defined arbitrary settling time. Considering non-linear engagement kinematics and a system lag to capture the effect of interceptor autopilot as present in realistic interception scenarios, the guidance schemes use the time-to-go estimates of the interceptors in order to achieve simultaneous interception of a stationary target at a pre-determined impact time. The guidance scheme ensures that consensus among the time-to-go estimates of the interceptors is achieved within a settling time whose upper bound can be pre-specified arbitrarily independent of the initial conditions or design parameters. The efficacy of the proposed guidance strategy is demonstrated using numerical simulations with varied conditions of initial position, velocities and heading angle errors of the interceptors as well as different desired impact times.


[271] 2407.15357

Lower bound for simulation cost of open quantum systems: Lipschitz continuity approach

Simulating quantum dynamics is one of the most promising applications of quantum computers. While the upper bound of the simulation cost has been extensively studied through various quantum algorithms, much less work has focused on establishing the lower bound, particularly for the simulation of open quantum system dynamics. In this work, we present a general framework to calculate the lower bound for simulating a broad class of quantum Markov semigroups. Given a fixed accessible unitary set, we introduce the concept of convexified circuit depth to quantify the quantum simulation cost and analyze the necessary circuit depth to construct a quantum simulation scheme that achieves a specific order. Our framework can be applied to both unital and non-unital quantum dynamics, and the tightness of our lower bound technique is illustrated by showing that the upper and lower bounds coincide in several examples.


[272] 2407.15417

Hemispheroidal parameterization and harmonic decomposition of simply connected open surfaces

Spectral analysis of open surfaces is gaining momentum for studying surface morphology in engineering, computer graphics, and medical domains. This analysis is enabled using proper parameterization approaches on the target analysis domain. In this paper, we propose the usage of customizable parameterization coordinates that allow mapping open surfaces into oblate or prolate hemispheroidal surfaces. For this, we proposed the usage of Tutte, conformal, area-preserving, and balanced mappings for parameterizing any given simply connected open surface onto an optimal hemispheroid. The hemispheroidal harmonic bases were introduced to spectrally expand these parametric surfaces by generalizing the known hemispherical ones. This approach uses the radius of the hemispheroid as a degree of freedom to control the size of the parameterization domain of the open surfaces while providing numerically stable basis functions. Several open surfaces have been tested using different mapping combinations. We also propose optimization-based mappings to serve various applications on the reconstruction problem. Altogether, our work provides an effective way to represent and analyze simply connected open surfaces.


[273] 2407.15448

Movable Antenna-Enhanced Wireless Communications: General Architectures and Implementation Methods

Movable antennas (MAs), traditionally explored in antenna design, have recently garnered significant attention in wireless communications due to their ability to dynamically adjust the antenna positions to changes in the propagation environment. However, previous research has primarily focused on characterizing the performance limits of various MA-assisted wireless communication systems, with less emphasis on their practical implementation. To address this gap, in this article, we propose several general MA architectures that extend existing designs by varying several key aspects to cater to different application scenarios and tradeoffs between cost and performance. Additionally, we draw from fields such as antenna design and mechanical control to provide an overview of candidate implementation methods for the proposed MA architectures, utilizing either direct mechanical or equivalent electronic control. Simulation results are finally presented to support our discussion.


[274] 2407.15455

Score matching for bridges without time-reversals

We propose a new algorithm for learning a bridged diffusion process using score-matching methods. Our method relies on reversing the dynamics of the forward process and using this to learn a score function, which, via Doob's $h$-transform, gives us a bridged diffusion process; that is, a process conditioned on an endpoint. In contrast to prior methods, ours learns the score term $\nabla_x \log p(t, x; T, y)$, for given $t, Y$ directly, completely avoiding the need for first learning a time reversal. We compare the performance of our algorithm with existing methods and see that it outperforms using the (learned) time-reversals to learn the score term. The code can be found at https://github.com/libbylbaker/forward_bridge.


[275] 2407.15466

Non-relativistic tachyons: a new representation of the Galilei group

An algebraic characterization of the contractions of the Poincar\'e group permits a proper construction of a non-relativistic limit of its tachyonic representation. We arrive at a consistent, nonstandard representation of the Galilei group which was disregarded long ago by supposedly unphysical properties. The corresponding quantum (and classical) theory shares with the relativistic one their fundamentals, and serves as a toy model to better comprehend the unusual behavior of the tachyonic representation. For instance, we see that evolution takes place in a spatial coordinate rather than time, as for relativistic tachyons, but the modulus of the three-momentum is the same for all Galilean observers, leading to a new dispersion relation for a Galilean system. Furthermore, the tachyonic objects described by the new representation cannot be regarded as localizable in the standard sense.


[276] 2407.15473

PyJama: Differentiable Jamming and Anti-Jamming with NVIDIA Sionna

Despite extensive research on jamming attacks on wireless communication systems, the potential of machine learning for amplifying the threat of such attacks, or our ability to mitigate them, remains largely untapped. A key obstacle to such research has been the absence of a suitable framework. To resolve this obstacle, we release PyJama, a fully-differentiable open-source library that adds jamming and anti-jamming functionality to NVIDIA Sionna. We demonstrate the utility of PyJama (i) for realistic MIMO simulations by showing examples that involve forward error correction, OFDM waveforms in time and frequency, realistic channel models, and mobility; and (ii) for learning to jam. Specifically, we use stochastic gradient descent to optimize jamming power allocation over an OFDM resource grid. The learned strategies are non-trivial, intelligible, and effective.


[277] 2407.15482

Quantum three-rotor problem in the identity representation

The quantum three-rotor problem concerns the dynamics of three equally massive particles moving on a circle subject to pairwise attractive cosine potentials and can model coupled Josephson junctions. Classically, it displays order-chaos-order behavior with increasing energy. The quantum system admits a dimensionless coupling with semiclassical behavior at strong coupling. We study stationary states with periodic `relative' wave functions. Perturbative and harmonic approximations capture the spectrum at weak coupling and that of low-lying states at strong coupling. More generally, the cumulative distribution of energy levels obtained by numerical diagonalization is well-described by a Weyl-like semiclassical estimate. However, the system has an $S_3 \times \mathbb{Z}_2$ symmetry that is obscured when working with relative angles. By exploiting a basis for invariant states, we obtain the purified spectrum in the identity representation. To uncover universal quantum hallmarks of chaos, we partition the spectrum into energy windows where the classical motion is regular, mixed or chaotic and unfold each separately. At strong coupling, we find striking signatures of transitions between regularity and chaos: spacing distributions morph from Poisson to Wigner-Dyson while the number variance shifts from linear to logarithmic behavior at small lengths. Some nonuniversal features are also examined. For instance, the number variance saturates and oscillates at large lengths for strong coupling while deviations from Poisson spacings at asymptotically low and high energies are well-explained by purified quantum harmonic and free-rotor spectra at strong and weak coupling. Interestingly, the degeneracy of free-rotor levels admits an elegant formula that we deduce using properties of Eisenstein primes.


[278] 2407.15492

Fast computation of 2-isogenies in dimension 4 and cryptographic applications

Dimension 4 isogenies have first been introduced in cryptography for the cryptanalysis of Supersingular Isogeny Diffie-Hellman (SIDH) and have been used constructively in several schemes, including SQIsignHD, a derivative of SQIsign isogeny based signature scheme. Unlike in dimensions 2 and 3, we can no longer rely on the Jacobian model and its derivatives to compute isogenies. In dimension 4 (and higher), we can only use theta-models. Previous works by Romain Cosset, David Lubicz and Damien Robert have focused on the computation of $\ell$-isogenies in theta-models of level $n$ coprime to $\ell$ (which requires to use $n^g$ coordinates in dimension $g$). For cryptographic applications, we need to compute chains of $2$-isogenies, requiring to use $\geq 3^g$ coordinates in dimension $g$ with state of the art algorithms. In this paper, we present algorithms to compute chains of $2$-isogenies between abelian varieties of dimension $g\geq 1$ with theta-coordinates of level $n=2$, generalizing a previous work by Pierrick Dartois, Luciano Maino, Giacomo Pope and Damien Robert in dimension $g=2$. We propose an implementation of these algorithms in dimension $g=4$ to compute endomorphisms of elliptic curve products derived from Kani's lemma with applications to SQIsignHD and SIDH cryptanalysis. We are now able to run a complete key recovery attack on SIDH when the endomorphism ring of the starting curve is unknown within a few seconds on a laptop for all NIST SIKE parameters.


[279] 2407.15504

Fundamental Limits of Prompt Compression: A Rate-Distortion Framework for Black-Box Language Models

We formalize the problem of prompt compression for large language models (LLMs) and present a framework to unify token-level prompt compression methods which create hard prompts for black-box models. We derive the distortion-rate function for this setup as a linear program, and provide an efficient algorithm to compute this fundamental limit via the dual of the linear program. Using the distortion-rate function as the baseline, we study the performance of existing compression schemes on a synthetic dataset consisting of prompts generated from a Markov chain, natural language queries, and their respective answers. Our empirical analysis demonstrates the criticality of query-aware prompt compression, where the compressor has knowledge of the downstream task/query for the black-box LLM. We show that there is a large gap between the performance of current prompt compression methods and the optimal strategy, and propose a query-aware, variable-rate adaptation of a prior work to close the gap. We extend our experiments to a small natural language dataset to further confirm our findings on our synthetic dataset.


[280] 2407.15518

An Efficient Regularity Lemma for Semi-Algebraic Hypergraphs

We use the polynomial method of Guth and Katz to establish stronger and {\it more efficient} regularity and density theorems for such $k$-uniform hypergraphs $H=(P,E)$, where $P$ is a finite point set in ${\mathbb R}^d$, and the edge set $E$ is determined by a semi-algebraic relation of bounded description complexity. In particular, for any $0<\epsilon\leq 1$ we show that one can construct in $O\left(n\log (1/\epsilon)\right)$ time, an equitable partition $P=U_1\uplus \ldots\uplus U_K$ into $K=O(1/\epsilon^{d+1+\delta})$ subsets, for any $0<\delta$, so that all but $\epsilon$-fraction of the $k$-tuples $U_{i_1},\ldots,U_{i_k}$ are {\it homogeneous}: we have that either $U_{i_1}\times\ldots\times U_{i_k}\subseteq E$ or $(U_{i_1}\times\ldots\times U_{i_k})\cap E=\emptyset$. If the points of $P$ can be perturbed in a general position, the bound improves to $O(1/\epsilon^{d+1})$, and the partition is attained via a {\it single partitioning polynomial} (albeit, at expense of a possible increase in worst-case running time). In contrast to the previous such regularity lemmas which were established by Fox, Gromov, Lafforgue, Naor, and Pach and, subsequently, Fox, Pach and Suk, our partition of $P$ does not depend on the edge set $E$ provided its semi-algebraic description complexity does not exceed a certain constant. As a by-product, we show that in any $k$-partite $k$-uniform hypergraph $(P_1\uplus\ldots\uplus P_k,E)$ of bounded semi-algebraic description complexity in ${\mathbb R}^d$ and with $|E|\geq \epsilon \prod_{i=1}^k|P_i|$ edges, one can find, in expected time $O\left(\sum_{i=1}^k|P_i|\log (1/\epsilon)+(1/\epsilon)\log(1/\epsilon)\right)$, subsets $Q_i\subseteq P_i$ of cardinality $|Q_i|\geq |P_i|/\epsilon^{d+1+\delta}$, so that $Q_1\times\ldots\times Q_k\subseteq E$.


[281] 2407.15567

A New Theoretical Perspective on Data Heterogeneity in Federated Optimization

In federated learning (FL), data heterogeneity is the main reason that existing theoretical analyses are pessimistic about the convergence rate. In particular, for many FL algorithms, the convergence rate grows dramatically when the number of local updates becomes large, especially when the product of the gradient divergence and local Lipschitz constant is large. However, empirical studies can show that more local updates can improve the convergence rate even when these two parameters are large, which is inconsistent with the theoretical findings. This paper aims to bridge this gap between theoretical understanding and practical performance by providing a theoretical analysis from a new perspective on data heterogeneity. In particular, we propose a new and weaker assumption compared to the local Lipschitz gradient assumption, named the heterogeneity-driven pseudo-Lipschitz assumption. We show that this and the gradient divergence assumptions can jointly characterize the effect of data heterogeneity. By deriving a convergence upper bound for FedAvg and its extensions, we show that, compared to the existing works, local Lipschitz constant is replaced by the much smaller heterogeneity-driven pseudo-Lipschitz constant and the corresponding convergence upper bound can be significantly reduced for the same number of local updates, although its order stays the same. In addition, when the local objective function is quadratic, more insights on the impact of data heterogeneity can be obtained using the heterogeneity-driven pseudo-Lipschitz constant. For example, we can identify a region where FedAvg can outperform mini-batch SGD even when the gradient divergence can be arbitrarily large. Our findings are validated using experiments.


[282] 2407.15580

Annealed Multiple Choice Learning: Overcoming limitations of Winner-takes-all with annealing

We introduce Annealed Multiple Choice Learning (aMCL) which combines simulated annealing with MCL. MCL is a learning framework handling ambiguous tasks by predicting a small set of plausible hypotheses. These hypotheses are trained using the Winner-takes-all (WTA) scheme, which promotes the diversity of the predictions. However, this scheme may converge toward an arbitrarily suboptimal local minimum, due to the greedy nature of WTA. We overcome this limitation using annealing, which enhances the exploration of the hypothesis space during training. We leverage insights from statistical physics and information theory to provide a detailed description of the model training trajectory. Additionally, we validate our algorithm by extensive experiments on synthetic datasets, on the standard UCI benchmark, and on speech separation.


[283] 2407.15678

Minimal-work protocols for inertial particles in non-harmonic traps

The progress of miniaturized technology allows controlling physical systems at nanoscale with remarkable precision in regimes where thermal fluctuations are non-negligible. Experimental advancements have sparked interest in control problems in stochastic thermodynamics, typically concerning a time-dependent potential applied to a nanoparticle to reach a target stationary state in a given time with minimal energy cost. We study this problem for a particle subject to thermal fluctuations in a regime that takes into account the effects of inertia, and, building on the results of [1], provide a numerical method to find optimal controls even for non-Gaussian initial and final conditions corresponding to non-harmonic confinements. We show that the momentum mean tends to a constant value along the trajectory except at the boundary and the evolution of the variance is non-trivial. Our results also support that the lower bound on the optimal entropy production computed from the overdamped case is tight in the adiabatic limit.


[284] 2407.15697

Voltage mapping in subcellular nanodomains using electro-diffusion modeling

Voltage distribution in sub-cellular micro-domains such as neuronal synapses, small protrusions or dendritic spines regulates the opening and closing of ionic channels, energy production and thus cellular homeostasis and excitability. Yet how voltage changes at such a small scale in vivo remains challenging due to the experimental diffraction limit, large signal fluctuations and the still limited resolution of fast voltage indicators. Here, we study the voltage distribution in nano-compartments using a computational approach based on the Poisson-Nernst-Planck equations for the electro-diffusion motion of ions, where inward and outward fluxes are generated between channels. We report a current-voltage (I-V) logarithmic relationship generalizing Nernst law that reveals how the local membrane curvature modulates the voltage. We further find that an influx current penetrating a cellular electrolyte can lead to perturbations from tens to hundreds of nanometers deep depending on the local channels organization. Finally, we show that the neck resistance of dendritic spines can be completely shunted by the transporters located on the head boundary, facilitating ionic flow. To conclude, we propose that voltage is regulated at a subcellular level by channels organization, membrane curvature and narrow passages.


[285] 2407.15730

Neural-based Video Compression on Solar Dynamics Observatory Images

NASA's Solar Dynamics Observatory (SDO) mission collects extensive data to monitor the Sun's daily activity. In the realm of space mission design, data compression plays a crucial role in addressing the challenges posed by limited telemetry rates. The primary objective of data compression is to facilitate efficient data management and transmission to work within the constrained bandwidth, thereby ensuring that essential information is captured while optimizing the utilization of available resources. This paper introduces a neural video compression technique that achieves a high compression ratio for the SDO's image data collection. The proposed approach focuses on leveraging both temporal and spatial redundancies in the data, leading to a more efficient compression. In this work, we introduce an architecture based on the Transformer model, which is specifically designed to capture both local and global information from input images in an effective and efficient manner. Additionally, our network is equipped with an entropy model that can accurately model the probability distribution of the latent representations and improves the speed of the entropy decoding step. The entropy model leverages a channel-dependent approach and utilizes checkerboard-shaped local and global spatial contexts. By combining the Transformer-based video compression network with our entropy model, the proposed compression algorithm demonstrates superior performance over traditional video codecs like H.264 and H.265, as confirmed by our experimental results.


[286] 2407.15740

The syzygy distinguisher

We present a new distinguisher for alternant and Goppa codes, whose complexity is subexponential in the code length. It does not suffer from the strong regime limitations of the previous distinguishers or structure recovery algorithms: in particular, it applies to the codes used in the Classic McEliece candidate for postquantum cryptography standardization. The invariants that allow us to distinguish are graded Betti numbers of the homogeneous coordinate ring of a shortening of the dual code. Since its introduction in 1978, this is the first time an analysis of the McEliece cryptosystem breaks the exponential barrier.


[287] 2407.15784

Diffusion Model Based Resource Allocation Strategy in Ultra-Reliable Wireless Networked Control Systems

Diffusion models are vastly used in generative AI, leveraging their capability to capture complex data distributions. However, their potential remains largely unexplored in the field of resource allocation in wireless networks. This paper introduces a novel diffusion model-based resource allocation strategy for Wireless Networked Control Systems (WNCSs) with the objective of minimizing total power consumption through the optimization of the sampling period in the control system, and blocklength and packet error probability in the finite blocklength regime of the communication system. The problem is first reduced to the optimization of blocklength only based on the derivation of the optimality conditions. Then, the optimization theory solution collects a dataset of channel gains and corresponding optimal blocklengths. Finally, the Denoising Diffusion Probabilistic Model (DDPM) uses this collected dataset to train the resource allocation algorithm that generates optimal blocklength values conditioned on the channel state information (CSI). Via extensive simulations, the proposed approach is shown to outperform previously proposed Deep Reinforcement Learning (DRL) based approaches with close to optimal performance regarding total power consumption. Moreover, an improvement of up to eighteen-fold in the reduction of critical constraint violations is observed, further underscoring the accuracy of the solution.