New articles on Mathematics


[1] 2106.09060

On the stability of the $L^{2}$ projection and the quasiinterpolant in the space of smooth periodic splines

In this paper we derive stability estimates in $L^{2}$- and $L^{\infty}$- based Sobolev spaces for the $L^{2}$ projection and a family of quasiinterolants in the space of smooth, 1-periodic, polynomial splines defined on a uniform mesh in $[0,1]$. As a result of the assumed periodicity and the uniform mesh, cyclic matrix techniques and suitable decay estimates of the elements of the inverse of a Gram matrix associated with the standard basis of the space of splines, are used to establish the stability results.


[2] 2106.09062

Local well-posedness for the Boltzmann equation with very soft potential and polynomially decaying initial data

In this paper, we address the local well-posedness of the spatially inhomogeneous non-cutoff Boltzmann equation when the initial data decays polynomially in the velocity variable. We consider the case of very soft potentials $\gamma + 2s < 0$. Our main result completes the picture for local well-posedness in this decay class by removing the restriction $\gamma + 2s > -3/2$ of previous works. % As an application, we can combine our existence result with the recent conditional regularity estimates of Imbert-Silvestre (arXiv:1909.12729 [math.AP]) to conclude solutions can be continued for as long as the mass, energy, and entropy densities remain under control. Our approach is entirely based on the Carleman decomposition of the collision operator into a lower order term and an integro-differential operator similar to the fractional Laplacian. Interestingly, this yields a very short proof of local well-posedness when $\gamma \in (-3,0]$ and $s \in (0,1/2)$ in a weighted $C^1$ space that we include as well.


[3] 2106.09066

Asymptotic shape of the concave majorant of a Lévy process

We establish distributional limit theorems for the shape statistics of a concave majorant (i.e. the fluctuations of its length, its supremum, the time it is attained and its value at $T$) of any L\'evy process on $[0,T]$ as $T\to\infty$. The scale of the fluctuations of the length and other statistics, as well as their asymptotic dependence, vary significantly with the tail behaviour of the L\'evy measure. The key tool in the proofs is the recent representation of the concave majorant for all L\'evy processes using a stick-breaking representation.


[4] 2106.09067

Very Well-Covered Graphs with the Erdös-Ko-Rado Property

A family of independent $r$-sets of a graph $G$ is an $r$-star if every set in the family contains some fixed vertex $v$. A graph is $r$-EKR if the maximum size of an intersecting family of independent $r$-sets is an $r$-star. Holroyd and Talbot conjecture that a graph is $r$-EKR as long as $1\leq r\leq\frac{\mu(G)}{2}$, where $\mu(G)$ is the minimum size of a maximal independent set. It is suspected that the smallest counterexample to this conjecture is a well-covered graph. Here we consider the class of very well-covered graphs $G^*$ obtained by appending a single pendant edge to each vertex of $G$. We prove that the pendant complete graph $K_n^*$ is $r$-EKR when $n \geq 2r$ and strictly so when $n>2r$. Pendant path graphs $P_n^*$ are also explored and the vertex whose $r$-star is of maximum size is determined.


[5] 2106.09074

Robust a posteriori error analysis for rotation-based formulations of the elasticity/poroelasticity coupling

We develop the \textit{a posteriori} error analysis of three mixed finite element formulations for rotation-based equations in elasticity, poroelasticity, and interfacial elasticity-poroelasticity. The discretisations use $H^1$-conforming finite elements of degree $k+1$ for displacement and fluid pressure, and discontinuous piecewise polynomials of degree $k$ for rotation vector, total pressure, and elastic pressure. Residual-based estimators are constructed, and upper and lower bounds (up to data oscillations) for all global estimators are rigorously derived. The methods are all robust with respect to the model parameters (in particular, the Lam\'e constants), they are valid in 2D and 3D, and also for arbitrary polynomial degree $k\geq 0$. The error behaviour predicted by the theoretical analysis is then demonstrated numerically on a set of computational examples including different geometries on which we perform adaptive mesh refinement guided by the \textit{a posteriori} error estimators.


[6] 2106.09077

Sharp upper and lower bounds of the attractor dimension for 3D damped Euler-Bardina equations

The dependence of the fractal dimension of global attractors for the damped 3D Euler--Bardina equations on the regularization parameter $\alpha>0$ and Ekman damping coefficient $\gamma>0$ is studied. We present explicit upper bounds for this dimension for the case of the whole space, periodic boundary conditions, and the case of bounded domain with Dirichlet boundary conditions. The sharpness of these estimates when $\alpha\to0$ and $\gamma\to0$ (which corresponds in the limit to the classical Euler equations) is demonstrated on the 3D Kolmogorov flows on a torus.


[7] 2106.09082

Zeroth-Order Methods for Convex-Concave Minmax Problems: Applications to Decision-Dependent Risk Minimization

Min-max optimization is emerging as a key framework for analyzing problems of robustness to strategically and adversarially generated data. We propose a random reshuffling-based gradient free Optimistic Gradient Descent-Ascent algorithm for solving convex-concave min-max problems with finite sum structure. We prove that the algorithm enjoys the same convergence rate as that of zeroth-order algorithms for convex minimization problems. We further specialize the algorithm to solve distributionally robust, decision-dependent learning problems, where gradient information is not readily available. Through illustrative simulations, we observe that our proposed approach learns models that are simultaneously robust against adversarial distribution shifts and strategic decisions from the data sources, and outperforms existing methods from the strategic classification literature.


[8] 2106.09083

Anisotropic percolation in high dimensions: the non-oriented case

We consider inhomogeneous non-oriented Bernoulli bond percolation on $\Z^d$, where each edge has a parameter, depending on its direction. We prove that, under certain conditions, if the sum of the parameters is strictly greater than 1/2, we have percolation for sufficiently high dimensions. The main tool is a dynamical coupling between the models in different dimensions with different set of parameters.


[9] 2106.09092

Norm inequalities for the spectral spread of Hermitian operators

In this work we introduce a new measure of dispersion of the eigenvalues of a Hermitian (self-adjoint) compact operator, that we call spectral spread. Then, we obtain several inequalities for unitarily invariant norms involving the spectral spread of self-adjoint compact operators, that are related with Bhatia-Davis's and Corach-Porta-Recht's work on Arithmetic-Geometric mean inequalities, Zhan's inequality for the difference of positive compact operators, Tao's inequalities for anti-diagonal blocks of positive compact operators and Kittaneh's commutator inequalities for positive operators.


[10] 2106.09099

Expanding measure has nonuniform specification property on random dynamical system

In the present paper, we study the distribution of the return points inthe fibers for a RDS (random dynamical systems) nonuniformly expanding preserving an ergodic probability, we also show the abundance of nonla-cunarity of hyperbolic times that are obtained along the orbits throughthe fibers. We conclude that any ergodic measure with positive Lyapunov exponents satisfies the nonuniform specification property between fibers.As consequences, we prove that any expanding measure is the limit of invariant measure whose measures of disintegration on the fibers are sup-ported by a finite number of return points and we prove that the average of the measures on the fibers corresponding to a disintegration, along an orbit in the base dynamics is the limit of Dirac measures supported inreturn orbits on the fibers.


[11] 2106.09101

Convex geometry of finite exchangeable laws and de Finetti style representation with universal correlated corrections

We present a novel analogue for finite exchangeable sequences of the de Finetti, Hewitt and Savage theorem and investigate its implications for multi-marginal optimal transport (MMOT) and Bayesian statistics. If $(Z_1,...,Z_N)$ is a finitely exchangeable sequence of $N$ random variables taking values in some Polish space $X$, we show that the law $\mu_k$ of the first $k$ components has a representation of the form $\mu_k=\int_{{\mathcal P}_{\frac{1}{N}}(X)} F_{N,k}(\lambda) \, \mbox{d} \alpha(\lambda)$ for some probability measure $\alpha$ on the set of $1/N$-quantized probability measures on $X$ and certain universal polynomials $F_{N,k}$. The latter consist of a leading term $N^{k-1}\! /{\small \prod_{j=1}^{k-1}(N\! -\! j)\, \lambda^{\otimes k}}$ and a finite, exponentially decaying series of correlated corrections of order $N^{-j}$ ($j=1,...,k$). The $F_{N,k}(\lambda)$ are precisely the extremal such laws, expressed via an explicit polynomial formula in terms of their one-point marginals $\lambda$. Applications include novel approximations of MMOT via polynomial convexification and the identification of the remainder which is estimated in the celebrated error bound of Diaconis-Freedman between finite and infinite exchangeable laws.


[12] 2106.09103

Approximately invertible elements in non-unital normed algebras

We introduce a concept of approximately invertible elements in non-unital normed algebras which is, on one side, a natural generalization of invertibility when having approximate identities at hand, and, on the other side, it is a direct extension of topological invertibility to non-unital algebras. Basic observations relate approximate invertibility with concepts of topological divisors of zero and density of (modular) ideals. We exemplify approximate invertibility in the group algebra, Wiener algebras, and operator ideals. For Wiener algebras with approximate identities (in particular, for the Fourier image of the convolution algebra), the approximate invertibility of an algebra element is equivalent to the property that it does not vanish. We also study approximate invertibility and its deeper connection with the Gelfand and representation theory in non-unital abelian Banach algebras as well as abelian and non-abelian C*-algebras.


[13] 2106.09107

Right triangulated categories: As extriangulated categories, aisles and co-aisles

Right triangulated categories can be thought of as triangulated categories whose shift functor is not an equivalence. We give intrinsic characterisations of when such categories have a natural extriangulated structure and are appearing as the (co-)aisle of a (co-)t-structure in an associated triangulated category.


[14] 2106.09116

Periodic Points of Ward-Veech Surfaces

For a non-arithmetic Veech surface, it is known that the set points having finite orbit under the Veech group, called the set of periodic points, is finite. However, few examples of these periodic point sets have been computed. In what follows, we compute new examples of periodic point sets for a family of non-hyperelliptic Veech surfaces with unbounded complexity and consider applications to classifying the holomorphic sections of certain surface bundles as well as applications to the study of infinitely generated Veech groups.


[15] 2106.09118

Locally compact sofic groups

We introduce the notion of soficity for locally compact groups and list a number of open problems.


[16] 2106.09120

A twisted Yu construction, Harish-Chandra characters, and endoscopy

We give a modification of Yu's construction of supercuspidal representations of a connected reductive group over a non-archimedean local field. This modification restores the validity of certain key intertwining property claims made by Yu, which were recently proven to be false for the original construction. This modification is also an essential ingredient in the explicit construction of supercuspidal L-packets. As further applications, we prove the stability and many instances of endoscopic character identities of these supercuspidal L-packets, subject to some conditions on the base field. In particular, for regular supercuspidal parameters we prove all instances of standard endoscopy. In addition, we prove that these supercuspidal L-packets satisfy a certain property, which, together with standard endoscopy, uniquely characterizes the local Langlands correspondence for supercuspidal L-packets (again subject to the above mentioned conditions on the base field). These results are based on a statement of the Harish-Chandra character formula for the supercuspidal representations arising from the twisted Yu construction.


[17] 2106.09123

Second-Order Conic and Polyhedral Approximations of the Exponential Cone: Application to Mixed-Integer Exponential Conic Programs

Exponents and logarithms exist in many important applications such as logistic regression, maximum likelihood, relative entropy and so on. Since the exponential cone can be viewed as the epigraph of perspective of the natural exponential function or the hypograph of perspective of the natural logarithm function, many mixed-integer nonlinear convex programs involving exponential or logarithm functions can be recast as mixed-integer exponential conic programs (MIECPs). Recently, solver MOSEK is able to solve large-scale continuous exponential conic programs (ECPs). However, unlike mixed-integer linear programs (MILPs) and mixed-integer second-order conic programs (MISOCPs), MIECPs are far beyond development. To harvest the past efforts on MILPs and MISOCPs, this paper presents second-order conic (SOC) and poly-hedral approximation schemes for the exponential cone with application to MIECPs. To do so, we first extend and generalize existing SOC approximation approaches in the extended space, propose new scaling and shifting methods, prove approximation accuracies, and derive lower bounds of approximations. We then study the polyhedral outer approximation of the exponential cones in the original space using gradient inequalities, show its approximation accuracy, and derive a lower bound of the approximation. When implementing SOC approximations, we suggest testing smaller cases by learning the approximation pattern and then applying to the large-scale cases; and for the polyhedral approximation, we suggest using the cutting plane method when solving the continuous ECP and branch and cut method for MIECPs. Our numerical study shows that the proposed scaling, shifting, and polyhedral outer approximation methods outperform solver MOSEK for both continuous ECPs and MIECPs and can achieve up to 20 times speed-ups compared to solver MOSEK when solving MIECPs.


[18] 2106.09125

Convex Optimization for Trajectory Generation

Reliable and efficient trajectory generation methods are a fundamental need for autonomous dynamical systems of tomorrow. The goal of this article is to provide a comprehensive tutorial of three major convex optimization-based trajectory generation methods: lossless convexification (LCvx), and two sequential convex programming algorithms known as SCvx and GuSTO. In this article, trajectory generation is the computation of a dynamically feasible state and control signal that satisfies a set of constraints while optimizing key mission objectives. The trajectory generation problem is almost always nonconvex, which typically means that it is not readily amenable to efficient and reliable solution onboard an autonomous vehicle. The three algorithms that we discuss use problem reformulation and a systematic algorithmic strategy to nonetheless solve nonconvex trajectory generation tasks through the use of a convex optimizer. The theoretical guarantees and computational speed offered by convex optimization have made the algorithms popular in both research and industry circles. To date, the list of applications includes rocket landing, spacecraft hypersonic reentry, spacecraft rendezvous and docking, aerial motion planning for fixed-wing and quadrotor vehicles, robot motion planning, and more. Among these applications are high-profile rocket flights conducted by organizations like NASA, Masten Space Systems, SpaceX, and Blue Origin. This article aims to give the reader the tools and understanding necessary to work with each algorithm, and to know what each method can and cannot do. A publicly available source code repository supports the provided numerical examples. By the end of the article, the reader should be ready to use the methods, to extend them, and to contribute to their many exciting modern applications.


[19] 2106.09126

The minimal ramification problem for rational function fields over finite fields

We study the minimal number of ramified primes in Galois extensions of rational function fields over finite fields with prescribed finite Galois group. In particular, we obtain a general conjecture in analogy with the well studied case of number fields, which we establish for abelian, symmetric and alternating groups in many cases.


[20] 2106.09131

Compactness for $Ω$-Yang-Mills connections

On a Riemannian manifold of dimension $n$ we extend the known analytic results on Yang-Mills connections to the class of connections called $\Omega$-Yang-Mills connections, where $\Omega$ is a smooth, not necessarily closed, $(n-4)$-form. Special cases include $\Omega$-anti-self-dual connections and Hermitian-Yang-Mills connections over general complex manifolds. By a key observation, a weak compactness result is obtained for moduli space of smooth $\Omega$-Yang-Mills connections with uniformly $L^2$ bounded curvature, and it can be improved in the case of Hermitian-Yang-Mills connections over general complex manifolds. A removable singularity theorem for singular $\Omega$-Yang-Mills connections on a trivial bundle with small energy concentration is also proven. As an application, it is shown how to compactify the moduli space of smooth Hermitian-Yang-Mills connections on unitary bundles over a class of balanced manifolds of Hodge-Riemann type. This class includes the metrics coming from multipolarizations, and in particular, the Kaehler metrics. In the case of multipolarizations on a projective algebraic manifold, the compactification of smooth irreducible Hermitian-Yang-Mills connections with fixed determinant modulo gauge transformations inherits a complex structure from algebro-geometric considerations.


[21] 2106.09133

The nonabelian Hodge correspondence for balanced hermitian metrics of Hodge-Riemann type

This paper extends the nonabelian Hodge correspondence for Kaehler manifolds to a larger class of hermitian metrics on complex manifolds called balanced of Hodge-Riemann type. Essentially, it grows out of a few key observations so that the known results, especially the Donaldson-Uhlenbeck-Yau theorem and Corlette's theorem, can be applied in our setting. Though not necessarily Kaehler, we show that the Sampson-Siu Theorem proving that harmonic maps are pluriharmonic remains valid for a slightly smaller class by using the known argument. Special important examples include those balanced metrics arising from multipolarizations.


[22] 2106.09136

Binary classification with corrupted labels

In a binary classification problem where the goal is to fit an accurate predictor, the presence of corrupted labels in the training data set may create an additional challenge. However, in settings where likelihood maximization is poorly behaved-for example, if positive and negative labels are perfectly separable-then a small fraction of corrupted labels can improve performance by ensuring robustness. In this work, we establish that in such settings, corruption acts as a form of regularization, and we compute precise upper bounds on estimation error in the presence of corruptions. Our results suggest that the presence of corrupted data points is beneficial only up to a small fraction of the total sample, scaling with the square root of the sample size.


[23] 2106.09139

Blow-up and scattering for the 1D NLS with point nonlinearity above the mass-energy threshold

In this paper, we study the nonlinear Schr\"odinger equation with focusing point nonlinearity in dimension one. First, we establish a scattering criterion for the equation based on Kenig-Merle's compactness-rigidity argument. Then we prove the energy scattering below and above the mass-energy threshold. We also describe the dynamics of solutions with data at the ground state threshold. Finally, we prove a blow-up criteria for the equation with initial data with arbitrarily large energy.


[24] 2106.09142

A Real-World Markov Chain arising in Recreational Volleyball

Card shuffling models have provided simple motivating examples for the mathematical theory of mixing times for Markov chains. As a complement, we introduce a more intricate realistic model of a certain observable real-world scheme for mixing human players onto teams. We quantify numerically the effectiveness of this mixing scheme over the 7 or 8 steps performed in practice. We give a combinatorial proof of the non-trivial fact that the chain is indeed irreducible.


[25] 2106.09143

Staircase symmetries in Hirzebruch surfaces

This paper continues the investigation of staircases in the family of Hirzebruch surfaces formed by blowing up the projective plane with weight b, that was started in Bertozzi, Holm et al. in arXiv:2010.08567. We explain the symmetries underlying the structure of the set of b that admit staircases, and show how the properties of these symmetries arise from a governing Diophantine equation. We also greatly simplify the techniques needed to show that a family of steps does form a staircase by using arithmetic properties of the accumulation function. There should be analogous results about both staircases and mutations for the other rational toric domains considered, for example, by Cristofaro-Gardiner et al. in arXiv:2004.07829 and by Casals--Vianna in arXiv:2004.13232.


[26] 2106.09145

Topological full groups of minimal subshifts and quantifying local embeddings into finite groups

We investigate quantitative aspects of the LEF property for subgroups of the topological full group $[[ \sigma ]]$ of a two-sided minimal subshift over a finite alphabet, measured via the LEF growth function. We show that the LEF growth of $[[ \sigma ]]^{\prime}$ may be bounded from above and below in terms of the recurrence function and the complexity function of the subshift, respectively. As an application, we construct groups of previously unseen LEF growth types, and exhibit a continuum of finitely generated LEF groups which may be distinguished from one another by their LEF growth.


[27] 2106.09147

Efficient recurrence for the enumeration of permutations with fixed pinnacle set

Initiated by Davis, Nelson, Petersen and Tenner (2018), the enumerative study of pinnacle sets of permutations has attracted a fair amount of attention recently. In this article, we provide a recurrence that can be used to compute efficiently the number $|\mathfrak{S}_n(P)|$ of permutations of size $n$ with a given pinnacle set $P$, with arithmetic complexity $O(k^4 + k\log n)$ for $P$ of size $k$. A symbolic expression can also be computed in this way for pinnacle sets of fixed size. A weighted sum $q_n(P)$ of $|\mathfrak{S}_n(P)|$ proposed in Davis, Nelson, Petersen and Tenner (2018) seems to have a simple form, and a conjectural form is given recently by Flaque, Novelli and Thibon (2021+). We settle the problem by providing and proving an alternative form of $q_n(P)$, which has a strong combinatorial flavor.


[28] 2106.09149

Fréchet derivatives of expected functionals of solutions to stochastic differential equations

In the analysis of stochastic dynamical systems described by stochastic differential equations (SDEs), it is often of interest to analyse the sensitivity of the expected value of a functional of the solution of the SDE with respect to perturbations in the SDE parameters. In this paper, we consider path functionals that depend on the solution of the SDE up to a stopping time. We derive formulas for Fr\'{e}chet derivatives of the expected values of these functionals with respect to bounded perturbations of the drift, using the Cameron-Martin-Girsanov theorem for the change of measure. Using these derivatives, we construct an example to show that the map that sends the change of drift to the corresponding relative entropy is not in general convex. We then analyse the existence and uniqueness of solutions to stochastic optimal control problems defined on possibly random time intervals, as well as gradient-based numerical methods for solving such problems.


[29] 2106.09150

$k$-positivity of dual canonical basis elements from 1324- and 2143-avoiding Kazhdan-Lusztig immanants

In this note, we show that certain dual canonical basis elements of $\mathbb{C}[SL_m]$ are positive when evaluated on $k$-positive matrices, matrices whose minors of size $k \times k$ and smaller are positive. Skandera showed that all dual canonical basis elements of $\mathbb{C}[SL_m]$ can be written in terms of Kazhdan-Lusztig immanants, which were introduced by Rhoades and Skandera. We focus on the basis elements which are expressed in terms of Kazhdan-Lusztig immanants indexed by 1324- and 2143-avoiding permutations. This extends previous work of the authors on Kazhdan-Lusztig immanants and uses similar tools, namely Lewis Carroll's identity (also known as the Desnanot-Jacobi identity).


[30] 2106.09151

Recovery Guarantees for Time-varying Pairwise Comparison Matrices with Non-transitivity

Pairwise comparison matrices have received substantial attention in a variety of applications, especially in rank aggregation, the task of flattening items into a one-dimensional (and thus transitive) ranking. However, non-transitive preference cycles can arise in practice due to the fact that making a decision often requires a complex evaluation of multiple factors. In some applications, it may be important to identify and preserve information about the inherent non-transitivity, either in the pairwise comparison data itself or in the latent feature space. In this work, we develop structured models for non-transitive pairwise comparison matrices that can be exploited to recover such matrices from incomplete noisy data and thus allow the detection of non-transitivity. Considering that individuals' tastes and items' latent features may change over time, we formulate time-varying pairwise comparison matrix recovery as a dynamic skew-symmetric matrix recovery problem by modeling changes in the low-rank factors of the pairwise comparison matrix. We provide theoretical guarantees for the recovery and numerically test the proposed theory with both synthetic and real-world data.


[31] 2106.09156

Separated and complete adelic models for one-dimensional Noetherian tensor-triangulated categories

We prove the existence of various adelic-style models for rigidly small-generated tensor-triangulated categories whose Balmer spectrum is a one-dimensional Noetherian topological space. This special case of our general programme of giving adelic models is particularly concrete and accessible, and we illustrate it with examples from algebra, geometry, topology and representation theory.


[32] 2106.09158

Mixed-Integer Nonlinear Programming for State-based Non-Intrusive Load Monitoring

Energy disaggregation, known in the literature as Non-Intrusive Load Monitoring (NILM), is the task of inferring the energy consumption of each appliance given the aggregate signal recorded by a single smart meter. In this paper, we propose a novel two-stage optimization-based approach for energy disaggregation. In the first phase, a small training set consisting of disaggregated power profiles is used to estimate the parameters and the power states by solving a mixed integer programming problem. Once the model parameters are estimated, the energy disaggregation problem is formulated as a constrained binary quadratic optimization problem. We incorporate penalty terms that exploit prior knowledge on how the disaggregated traces are generated, and appliance-specific constraints characterizing the signature of different types of appliances operating simultaneously. Our approach is compared with existing optimization-based algorithms both on a synthetic dataset and on three real-world datasets. The proposed formulation is computationally efficient, able to disambiguate loads with similar consumption patterns, and successfully reconstruct the signatures of known appliances despite the presence of unmetered devices, thus overcoming the main drawbacks of the optimization-based methods available in the literature.


[33] 2106.09160

Large-scale regularity for the stationary Navier-Stokes equations over non-Lipschitz boundaries

In this paper we address the large-scale regularity theory for the stationary Navier-Stokes equations in highly oscillating bumpy John domains. These domains are very rough, possibly with fractals or cusps, at the microscopic scale, but are amenable to the mathematical analysis of the Navier-Stokes equations. We prove: (i) a large-scale Calder\'on-Zygmund estimate, (ii) a large-scale Lipschitz estimate, (iii) large-scale higher-order regularity estimates, namely, $C^{1,\gamma}$ and $C^{2,\gamma}$ estimates. These nice regularity results are inherited only at mesoscopic scales, and clearly fail in general at the microscopic scales. We emphasize that the large-scale $C^{1,\gamma}$ regularity is obtained by using first-order boundary layers constructed via a new argument. The large-scale $C^{2,\gamma}$ regularity relies on the construction of second-order boundary layers, which allows for certain boundary data with linear growth at spatial infinity. To the best of our knowledge, our work is the first to carry out such an analysis. In the wake of many works in quantitative homogenization, our results strongly advocate in favor of considering the boundary regularity of the solutions to fluid equations as a multiscale problem, with improved regularity at or above a certain scale.


[34] 2106.09165

Stability properties of complete self-shrinking surfaces in $\mathbb{R}^3$

This paper studies rigidity for immersed self-shrinkers of the mean curvature flow of surfaces in the three-dimensional Euclidean space $\mathbb{R}^3.$ We prove that an immersed self-shrinker with finite $L$-index must be proper and of finite topology. As one of consequences, there is no stable two-dimensional self-shrinker in $\mathbb{R}^3$ without assuming properness. We conclude the paper by giving an affirmative answer to a question of Mantegazza.


[35] 2106.09172

On integral conditions for the existence of first integrals analytic saddle singularities

We study one-parameter analytic integrable deformations of the germ of $2\times(n-2)$-type complex saddle singularity given by $d(xy)=0$ at the origin $0 \in \mathbb C^2\times \mathbb C^{n-2}$. Such a deformation writes ${\omega}^t=d(xy) + \sum\limits_{j=1}^\infty t^j \omega_j$ where $t\in \mathbb C,0$ is the parameter of the deformation and the coefficients $\omega_j$ are holomorphic one-forms in some neighborhood of the origin $0\in \mathbb C^n$. We prove that, under a nondegeneracy condition of the singular set of the deformation, with respect to the fibration $d(xy)=0$, the existence of a holomorphic first integral for each element ${\omega}^t$ of the deformation is equivalent to the vanishing of certain line integrals $\oint_{\gamma_c}{\omega}^t=0, \forall \gamma_c, \forall t$ calculated on cycles $\gamma_c$ contained in the fibers $xy=c, \,0 \ne c \in \mathbb C,0$. This result is quite sharp regarding the conditions of the singular set and on the vanishing of the integrals in cycles. It is also not valid for ramified saddles, i.e., for deformations of saddles of the form $x^ny^m=c$ where $n+m>2$. As an application of our techniques we obtain a criteria for the existence of first integrals for integrable codimension one deformations of quadratic analytic center-cylinder type singularities in terms of the vanishing of some easy to compute line integrals.


[36] 2106.09175

Efficient and accurate KAM tori construction for the dissipative spin-orbit problem using a map reduction

We consider the dissipative spin-orbit problem in Celestial Mechanics, which describes the rotational motion of a triaxial satellite moving on a Keplerian orbit subject to tidal forcing and "drift". Our goal is to construct quasi-periodic solutions with fixed frequency, satisfying appropriate conditions. With the goal of applying rigorous KAM theory, we compute such quasi-periodic solution with very high precision. To this end, we have developed a very efficient algorithm. The first step is to compute very accurately the return map to a surface of section (using a high order Taylor's method with extended precision). Then, we find an invariant curve for the return map using recent algorithms that take advantage of the geometric features of the problem. This method is based on a rapidly convergent Newton's method which is guaranteed to converge if the initial error is small enough. So, it is very suitable for a continuation algorithm. The resulting algorithm is quite efficient. We only need to deal with a one dimensional function. If this function is discretized in $N$ points, the algorithm requires $O(N \log N) $ operations and $O(N) $ storage. The most costly step (the numerical integration of the equation along a turn) is trivial to parallelize. The main goal of the paper is to present the algorithms, implementation details and several sample results of runs. We also present both a rigorous and a numerical comparison of the results of averaged and not averaged models.


[37] 2106.09181

A cyclotomic family of thin hypergeometric monodromy groups in ${Sp}_4(\mathbb{R})$

We exhibit an infinite family of discrete subgroups of ${Sp}_4(\mathbb R)$ which have a number of remarkable properties. Our results are established by showing that each group plays ping-pong on an appropriate set of cones. The groups arise as the monodromy of hypergeometric differential equations with parameters $\left(\tfrac{N-3}{2N},\tfrac{N-1}{2N}, \tfrac{N+1}{2N}, \tfrac{N+3}{2N}\right)$ at infinity and maximal unipotent monodromy at zero, for any integer $N\geq 4$. Additionally, we relate the cones used for ping-pong in $\mathbb R^4$ with crooked surfaces, which we then use to exhibit domains of discontinuity for the monodromy groups in the Lagrangian Grassmannian.


[38] 2106.09182

The cohomology of left-invariant CR structures of hypersurface type on compact Lie groups

Pittie [Pit88] proved that the Dolbeault cohomology of all left-invariant complex structures on compact Lie groups can be computed by looking at the Dolbeault cohomology induced on a conveniently chosen maximal torus. We use the algebraic classification of left-invariant CR structures of maximal rank on compact Lie groups [CK04] to generalize Pittie's result to left-invariant Levi-flat CR structures of maximal rank on compact Lie groups.


[39] 2106.09183

Global Dynamics of a Predator-Prey Model with State-Dependent Maturation-Delay

In this paper, a stage structured predator-prey model with general nonlinear type of functional response is established and analyzed. The state-dependent time delay (hereafter SDTD) is the time taken from predator's birth to its maturity, formatted as a monotonical (ly) increasing, continuous(ly) differentiable and bounded function on the number of mature predator. The model is quite different from many previous models with SDTD, in the sense that the derivative of delay on the time is involved in the model. First, we have shown that for a large class of commonly used types of functional responses, including Holling types I, II and III, Beddington-DeAngelis-type (hereafter BD-type), etc, the predator coexists with the prey permanently if and only if the predator's net reproduction number is larger than one unit; Secondly, we have discussed the local stability of the equilibria of the model; Finally, for the special case of BD-type functional response, we claim that if the system is permanent, that is, the derivative of SDTD on the state is small enough and the predator interference is large enough, then the coexistence equilibrium is globally asymptotically stable.


[40] 2106.09184

A fourth-order compact time-splitting method for the Dirac equation with time-dependent potentials

In this paper, we present an approach to deal with the dynamics of the Dirac equation with time-dependent electromagnetic potentials using the fourth-order compact time-splitting method ($S_\text{4c}$). To this purpose, the time-ordering technique for time-dependent Hamiltonians is introduced, so that the influence of the time-dependence could be limited to certain steps which are easy to treat. Actually, in the case of the Dirac equation, it turns out that only those steps involving potentials need to be amended, and the scheme remains efficient, accurate, as well as easy to implement. Numerical examples in 1D and 2D are given to validate the scheme.


[41] 2106.09191

The Biot-Stokes coupling using total pressure: formulation, analysis and application to interfacial flow in the eye

We consider a multiphysics model for the flow of Newtonian fluid coupled with Biot consolidation equations through an interface, and incorporating total pressure as an unknown in the poroelastic region. A new mixed-primal finite element scheme is proposed solving for the pairs fluid velocity - pressure and displacement - total poroelastic pressure using Stokes-stable elements, and where the formulation does not require Lagrange multipliers to set up the usual transmission conditions on the interface. The stability and well-posedness of the continuous and semi-discrete problems are analysed in detail. Our numerical study {is framed in} the context of different interfacial flow regimes in Cartesian and axisymmetric coordinates that could eventually help describe early morphologic changes associated with glaucoma development in canine species.


[42] 2106.09195

On the Classifying Space for Commutativity in $U(3)$

For a Lie group $G$, let $B_{com} G$ be the classifying space for commutativity. Let $E_{com} G$ be the total space of the principal $G$-bundle associated to $B_{com} G$. In this article, we present a computation of the cohomology of $E_{com} U(3)$ using the spectral sequence associated to a homotopy colimit. As a part of our computation we will also compute the integral cohomology of $U(3)/N(T)$ and $(U(3)/T) \times_{\Sigma_3} (U(3)/T)$ where $T$ is a maximal torus of $U(3)$ with normalizer $N(T)$.


[43] 2106.09200

Automorphism groups over a hyperimaginary

In this paper we study the Lascar group over a hyperimaginary e. We verify that various results about the group over a real set still hold when the set is replaced by e. First of all, there is no written proof in the available literature that the group over e is a topological group. We present an expository style proof of the fact, which even simplifies existing proofs for the real case. We further extend a result that the orbit equivalence relation under a closed subgroup of the Lascar group is type-definable. On the one hand, we correct errors appeared in the book, "Simplicity Theory" [6, 5.1.14-15] and produce a counterexample. On the other, we extend Newelski's Theorem in "The diameter of a Lascar strong type" [12] that `a G-compact theory over a set has a uniform bound for the Lascar distances' to the hyperimaginary context. Lastly, we supply a partial positive answer to a question raised in "The relativized Lascar groups, type-amalgamations, and algebraicity" [4, 2.11], which is even a new result in the real context.


[44] 2106.09205

Improved bounds in Weaver's ${\rm KS}_r$ conjecture for high rank positive semidefinite matrices

Recently Marcus, Spielman and Srivastava proved Weaver's ${\rm{KS}}_r$ conjecture, which gives a positive solution to the Kadison-Singer problem. Cohen and Br\"and\'en independently extended this result to obtain the arbitrary-rank version of Weaver's ${\rm{KS}}_r$ conjecture. In this paper, we present a new bound in Weaver's ${\rm{KS}}_r$ conjecture for the arbitrary-rank case. To do that, we introduce the definition of $(k,m)$-characteristic polynomials and employ it to improve the previous estimate on the largest root of the mixed characteristic polynomials. For the rank-one case, our bound agrees with the Bownik-Casazza-Marcus-Speegle's bound when $r=2$ and with the Ravichandran-Leake's bound when $r>2$. For the higher-rank case, we sharpen the previous bounds from Cohen and from Br\"and\'en .


[45] 2106.09209

Some tight bounds on the minimum and maximum forcing numbers of graphs

Let $G$ be a simple graph with $2n$ vertices and a perfect matching. The forcing number of a perfect matching $M$ of $G$ is the smallest cardinality of a subset of $M$ that is contained in no other perfect matching of $G$. Let $f(G)$ and $F(G)$ denote the minimum and maximum forcing number of $G$ among all perfect matchings, respectively. Hetyei obtained that the maximum number of edges of graphs $G$ with a unique perfect matching is $n^2$ (see Lov\'{a}sz [20]). We know that $G$ has a unique perfect matching if and only if $f(G)=0$. Along this line, we generalize the classical result to all graphs $G$ with $f(G)=k$ for $0\leq k\leq n-1$, and obtain that the number of edges is at most $n^2+2nk-k^2-k$ and characterize the extremal graphs as well. Conversely, we get a non-trivial lower bound of $f(G)$ in terms of the order and size. For bipartite graphs, we gain corresponding stronger results. Further, we obtain a new upper bound of $F(G)$. Finally some open problems and conjectures are proposed.


[46] 2106.09213

The Affine Shape of a Figure-Eight under the Curve Shortening Flow

We give a precise account of the global limiting shape of a natural class of figure-eight curves under the curve-shortening flow. In addition to proving the rescaling convergence (on a certain subarc) to the Grim Reaper soliton, we take the novel approach of applying a family of non-conformal affine transformations to the solution so that the resulting images always have the unit square for a bounding box. We then prove that the limit converges in the Hausdorff metric to a quadrilateral called a bowtie. Our results enhance those of S. Angenent and M. Grayson on this topic.


[47] 2106.09220

Infinite-time blowing-up solutions to small perturbations of the Yamabe flow

Under the validity of the positive mass theorem, the Yamabe flow on a smooth compact Riemannian manifold of dimension $N \ge 3$ is known to exist for all time $t$ and converges to a solution to the Yamabe problem as $t \to \infty$. We prove that if a suitable perturbation, which may be smooth and arbitrarily small, is imposed on the linear term of the Yamabe flow on any given Riemannian manifold $M$ of dimension $N \ge 5$, the resulting flow may blow-up in the infinite time, forming singularities each of which looks like a solution of the Yamabe problem on the unit sphere $\mathbb{S}^N$. This shows that the Yamabe flow is an equation at the borderline guaranteeing the global existence and uniformly boundness of solutions. In addition, we examine the stability of the blow-up phenomena, imposing some negativity condition on the Ricci curvature at blow-up points.


[48] 2106.09228

Regularity structure of conservative solutions to the Hunter-Saxton equation

n this paper we characterize the regularity structure, as well as show the global-in-time existence and uniqueness, of (energy) conservative solutions to the Hunter-Saxton equation by using the method of characteristics. The major difference between the current work and previous results is that we are able to characterize the singularities of energy measure and their nature in a very precise manner. In particular, we show that singularities, whose temporal and spatial locations are also explicitly given in this work, may only appear at at most countably many times, and are completely determined by the absolutely continuous part of initial energy measure. Our mathematical analysis is based on using the method of characteristics in a generalized framework that consists of the evolutions of solution to the Hunter-Saxton equation and the energy measure. This method also provides a clear description of the semi-group property for the solution and energy measure for all times.


[49] 2106.09238

On the spectral radius of unicyclic and bicyclic graphs with a fixed diameter

The $\alpha$-spectral radius of a connected graph $G$ is the spectral radius of $A_\alpha$-matrix of $G$. In this paper, we discuss the methods for comparing $\alpha$-spectral radius of graphs. As applications, we characterize the graphs with the maximal $\alpha$-spectral radius among all unicyclic and bicyclic graphs of order $n$ with diameter $d$, respectively. Finally, we determine the unique graph with maximal signless Laplacian spectral radius among bicyclic graphs of order $n$ with diameter $d$.


[50] 2106.09247

Anders Wiman's "On the Application of Tschirnhaus Transformations to the Reduction of Algebraic Equations''

This is an English translation of Anders Wiman's classical paper "\"{U}ber die Anwendung der Tschirnhausen Transformation auf die Reduktion algebraischer Gleichungen" from 1927. The original work first appeared in the 1927 extraordinary edition of Nova Acta Regiae Societatis Scientiarum Upsaliensis. In this paper, Wiman gives an argument that the general polynomial of degree nine can be solved using one algebraic function of four variables and accessory irrationalities of degree at most five. However, his argument assumes certain intersections in affine space are generic without proof.


[51] 2106.09253

Stability of Caffarelli-Kohn-Nirenberg inequality

In this paper, we consider the Caffarelli-Kohn-Nirenberg (CKN) inequality: \begin{eqnarray*} \bigg(\int_{{\mathbb R}^N}|x|^{-b(p+1)}|u|^{p+1}dx\bigg)^{\frac{2}{p+1}}\leq C_{a,b,N}\int_{{\mathbb R}^N}|x|^{-2a}|\nabla u|^2dx \end{eqnarray*} where $N\geq3$, $a<\frac{N-2}{2}$, $a\leq b\leq a+1$ and $p=\frac{N+2(1+a-b)}{N-2(1+a-b)}$. It is well-known that up to dilations $\tau^{\frac{N-2}{2}-a}u(\tau x)$ and scalar multiplications $Cu(x)$, the CKN inequality has a unique extremal function $W(x)$ which is positive and radially symmetric in the parameter region $b_{FS}(a)\leq b<a+1$ with $a<0$ and $a\leq b<a+1$ with $a\geq0$ and $a+b>0$, where $b_{FS}(a)$ is the Felli-Schneider curve. We prove that in the above parameter region the following stabilities hold: \begin{enumerate} \item[$(1)$] \quad stability of CKN inequality in the functional inequality setting $$dist_{D^{1,2}_{a}}^2(u, \mathcal{Z})\lesssim\|u\|^2_{D^{1,2}_a({\mathbb R}^N)}-C_{a,b,N}^{-1}\|u\|^2_{L^{p+1}(|x|^{-b(p+1)},{\mathbb R}^N)}$$ where $\mathcal{Z}= \{ c W_\tau\mid c\in\bbr\backslash\{0\}, \tau>0\}$; \item[$(2)$]\quad stability of CKN inequality in the critical point setting (in the class of nonnegative functions) \begin{eqnarray*} dist_{D_a^{1,2}}(u, \mathcal{Z}_0^\nu)\lesssim\left\{\aligned &\Gamma(u),\quad p>2\text{ or }\nu=1,\\ &\Gamma(u)|\log\Gamma(u)|^{\frac12},\quad p=2\text{ and }\nu\geq2,\\ &\Gamma(u)^{\frac{p}{2}},\quad 1<p<2\text{ and }\nu\geq2, \endaligned\right. \end{eqnarray*} where $\Gamma (u)=\|div(|x|^{-a}\nabla u)+|x|^{-b(p+1)}|u|^{p-1}u\|_{(D^{1,2}_a)^{'}}$ and $$\mathcal{Z}_0^\nu=\{(W_{\tau_1},W_{\tau_2},\cdots,W_{\tau_\nu})\mid \tau_i>0\}.$$


[52] 2106.09254

On Hook Formulas for Cylindric Skew Diagrams

We present a conjectual hook formula concerning the number of the standard tableaux on "cylindric" skew diagrams. Our formula can be seen as an extension of Naruse's hook formula for skew diagrams. Moreover, we prove our conjecture in some special cases.


[53] 2106.09262

Componentwise linear ideals in Veronese rings

In this article, we study the componentwise linear ideals in the Veronese subrings of $R=K[x_1,\ldots,x_n]$. If char$(K)=0$, then we give a characterization for graded ideals in the $c^{th}$ Veronese ring $R^{(c)}$ to be componentwise linear. This characterization is an analogue of that over $R$ due to Aramova, Herzog and Hibi in \cite{ah00} and Conca, Herzog and Hibi in \cite{chh04}.


[54] 2106.09264

Fano 4-folds with a small contraction

Let X be a smooth complex Fano 4-fold. We show that if X has a small elementary contraction, then the Picard number rho(X) of X is at most 12. This result is based on a careful study of the geometry of X, on which we give a lot of information. We also show that in the boundary case rho(X)=12 an open subset of X has a smooth fibration with fiber the projective line. Together with previous results, this implies if X is a Fano 4-fold with rho(X)>12, then every elementary contraction of X is divisorial and sends a divisor to a surface. The proof is based on birational geometry and the study of families of rational curves. More precisely the main tools are: the study of families of lines in Fano 4-folds and the construction of divisors covered by lines, a detailed study of fixed prime divisors, the properties of the faces of the effective cone, and a detailed study of rational contractions of fiber type.


[55] 2106.09268

Heat kernel asymptotics for Kohn Laplacians on CR manifolds

Let $X$ be an abstract orientable not necessarily compact CR manifold of dimension $2n+1$, $n\geq1$, and let $L^k$ be the $k$-th tensor power of a CR complex line bundle $L$ over $X$. Suppose that condition $Y(q)$ holds at each point of $X$, we establish asymptotics of the heat kernel of Kohn Laplacian with values in $L^k$. As an application, we give a heat kernel proof of Morse inequalities on compact CR manifolds. When $X$ admits a transversal CR $\mathbb R$-action, we also establish asymptotics of the $\mathbb R$-equivariant heat kernel of Kohn Laplacian with values in $L^k$. As an application, we get $\mathbb R$-equivariant Morse inequalities on compact CR manifolds with transversal CR $\mathbb R$-action.


[56] 2106.09270

Differential models for the Anderson dual to bordism theories and invertible QFT's

In this paper, we construct new models for the Anderson duals $(I\Omega^G)^*$ to the stable tangential $G$-bordism theories and their differential extensions. The cohomology theory $(I\Omega^G)^*$ is conjectured by Freed and Hopkins \cite{Freed:2016rqq} to classify deformation classes of possibly non-topological invertible quantum field theories (QFT's). Our model is made by abstractizing certain properties of invertible QFT's, thus supporting their conjecture.


[57] 2106.09284

Reconstructing simplicial polytopes from their graphs and affine $2$-stresses

A conjecture of Kalai from 1994 posits that for an arbitrary $2\leq k\leq \lfloor d/2 \rfloor$, the combinatorial type of a simplicial $d$-polytope $P$ is uniquely determined by the $(k-1)$-skeleton of $P$ (given as an abstract simplicial complex) together with the space of affine $k$-stresses on $P$. We establish the first non-trivial case of this conjecture, namely, the case of $k=2$. We also prove that for a general $k$, Kalai's conjecture holds for the class of $k$-neighborly polytopes.


[58] 2106.09286

Sub-linear convergence of a tamed stochastic gradient descent method in Hilbert space

In this paper, we introduce the tamed stochastic gradient descent method (TSGD) for optimization problems. Inspired by the tamed Euler scheme, which is a commonly used method within the context of stochastic differential equations, TSGD is an explicit scheme that exhibits stability properties similar to those of implicit schemes. As its computational cost is essentially equivalent to that of the well-known stochastic gradient descent method (SGD), it constitutes a very competitive alternative to such methods. We rigorously prove (optimal) sub-linear convergence of the scheme for strongly convex objective functions on an abstract Hilbert space. The analysis only requires very mild step size restrictions, which illustrates the good stability properties. The analysis is based on a priori estimates more frequently encountered in a time integration context than in optimization, and this alternative approach provides a different perspective also on the convergence of SGD. Finally, we demonstrate the usability of the scheme on a problem arising in a context of supervised learning.


[59] 2106.09288

The Stark problem as a concave toric domain

The Stark problem is a completely integrable system which describes the motion of an electron in a constant electric field and subject to the attraction of a proton. In this paper we show that in the planar case after Levi-Civita regularization the bounded component of the energy hypersurfaces of the Stark problem for energies below the critical value can be interpreted as boundaries of concave toric domains.


[60] 2106.09294

Prescribing Morse scalar curvatures: incompatibility of non existence

Given a closed manifold of positive Yamabe invariant and for instance positive Morse functions upon it, the conformally prescribed scalar curvature problem raises the question, whether or not such functions can by conformally changing the metric be realised as the scalar curvature of this manifold. As we shall quantify, the answer is depending on the complexity of these functions most likely yes and, if one such Morse function happens not to be conformally prescribable, most others can and this in a controllable way.


[61] 2106.09299

Relaxed Lagrangian duality in convex infinite optimization: reverse strong duality and optimality

We associate with each convex optimization problem posed on some locally convex space with an infinite index set T, and a given non-empty family H formed by finite subsets of T, a suitable Lagrangian-Haar dual problem. We provide reverse H-strong duality theorems, H-Farkas type lemmas and optimality theorems. Special attention is addressed to infinite and semi-infinite linear optimization problems.


[62] 2106.09313

Counting Discrete, Level-$1$, Quaternionic Automorphic Representations on $G_2$

Quaternionic automorphic representations are one attempt to generalize to other groups the special place holomorphic modular forms have among automorphic representations of $\mathrm{GL}_2$. Here, we study quaternionic automorphic representations on the exceptional group $G_2$. Using "hyperendoscopy" techniques from arXiv:1910.10800, we develop for quaternionic, $G_2$-representations an analog of the Eichler-Selberg trace formula for classical modular forms. We then use this together with some techniques of Chenevier, Renard, and Ta\"ibi to compute dimensions of spaces of level-$1$ quaternionic representations. On the way, we prove a Jacquet-Langlands-style result describing them in terms of classical modular forms and automorphic representations on the compact-at-infinity form $G_2^c$. The main technical difficulty is that the quaternionic discrete series that quaternionic autmorphic representations are defined in terms of do not satisfy a condition of being "regular". Using some computations from arXiv:2010.02712, we show that this miraculously does not matter for specifically the case of quaternionic discrete series on $G_2$. We hope that this project serves as a guiding example for anyone interested in computing exact information about discrete-at-infinity automorphic representations on arbitrary reductive groups instead of just classical ones.


[63] 2106.09316

Optimized Power Control Design for Over-the-Air Federated Edge Learning

This paper investigates the transmission power control in over-the-air federated edge learning (Air-FEEL) system. Different from conventional power control designs (e.g., to minimize the individual mean squared error (MSE) of the over-the-air aggregation at each round), we consider a new power control design aiming at directly maximizing the convergence speed. Towards this end, we first analyze the convergence behavior of Air-FEEL (in terms of the optimality gap) subject to aggregation errors at different communication rounds. It is revealed that if the aggregation estimates are unbiased, then the training algorithm would converge exactly to the optimal point with mild conditions; while if they are biased, then the algorithm would converge with an error floor determined by the accumulated estimate bias over communication rounds. Next, building upon the convergence results, we optimize the power control to directly minimize the derived optimality gaps under both biased and unbiased aggregations, subject to a set of average and maximum power constraints at individual edge devices. We transform both problems into convex forms, and obtain their structured optimal solutions, both appearing in a form of regularized channel inversion, by using the Lagrangian duality method. Finally, numerical results show that the proposed power control policies achieve significantly faster convergence for Air-FEEL, as compared with benchmark policies with fixed power transmission or conventional MSE minimization.


[64] 2106.09328

Polaron models with regular interactions at strong coupling

We study a class of polaron-type Hamiltonians with sufficiently regular form factor in the interaction term. We investigate the strong-coupling limit of the model, and prove suitable bounds on the ground state energy as a function of the total momentum of the system. These bounds agree with the semiclassical approximation to leading order. The latter corresponds here to the situation when the particle undergoes harmonic motion in a potential well whose frequency is determined by the corresponding Pekar functional. We show that for all such models the effective mass diverges in the strong coupling limit, in all spatial dimensions. Moreover, for the case when the phonon dispersion relation grows at least linearly with momentum, the bounds result in an asymptotic formula for the effective mass quotient, a quantity generalizing the usual notion of the effective mass. This asymptotic form agrees with the semiclassical Landau--Pekar formula and can be regarded as the first rigorous confirmation, in a slightly weaker sense than usually considered, of the validity of the semiclassical formula for the effective mass.


[65] 2106.09332

On first and second order linear Stieltjes differential equations

This work deals with the obtaining of solutions of first and second order Stieltjes differential equations. We define the notions of Stieltjes derivative on the whole domain of the functions involved, provide a notion of n-times continuously Stieltjes-differentiable functions and prove existence and uniqueness results of Stieltjes differential equations in those spaces. We also present the Green's functions associated to the different problems and an application to the Stieltjes harmonic oscillator.


[66] 2106.09339

Optimal explicit stabilized postprocessed $τ$-leap method for the simulation of chemical kinetics

The simulation of chemical kinetics involving multiple scales constitutes a modeling challenge (from ordinary differential equations to Markov chain) and a computational challenge (multiple scales, large dynamical systems, time step restrictions). In this paper we propose a new discrete stochastic simulation algorithm: the postprocessed second kind stabilized orthogonal $\tau$-leap Runge-Kutta method (PSK-$\tau$-ROCK). In the context of chemical kinetics this method can be seen as a stabilization of Gillespie's explicit $\tau$-leap combined with a postprocessor. The stabilized procedure allows to simulate problems with multiple scales (stiff), while the postprocessing procedure allows to approximate the invariant measure (e.g. mean and variance) of ergodic stochastic dynamical systems. We prove stability and accuracy of the PSK-$\tau$-ROCK. Numerical experiments illustrate the high reliability and efficiency of the scheme when compared to other $\tau$-leap methods.


[67] 2106.09340

A trust region-type normal map-based semismooth Newton method for nonsmooth nonconvex composite optimization

We propose a novel trust region method for solving a class of nonsmooth and nonconvex composite-type optimization problems. The approach embeds inexact semismooth Newton steps for finding zeros of a normal map-based stationarity measure for the problem in a trust region framework. Based on a new merit function and acceptance mechanism, global convergence and transition to fast local q-superlinear convergence are established under standard conditions. In addition, we verify that the proposed trust region globalization is compatible with the Kurdyka-{\L}ojasiewicz (KL) inequality yielding finer convergence results. We further derive new normal map-based representations of the associated second-order optimality conditions that have direct connections to the local assumptions required for fast convergence. Finally, we study the behavior of our algorithm when the Hessian matrix of the smooth part of the objective function is approximated by BFGS updates. We successfully link the KL theory, properties of the BFGS approximations, and a Dennis-Mor{\'e}-type condition to show superlinear convergence of the quasi-Newton version of our method. Numerical experiments on sparse logistic regression and image compression illustrate the efficiency of the proposed algorithm.


[68] 2106.09341

Positivity for the clamped plate equation under high tension

In this article we consider positivity issues for the clamped plate equation with high tension $\gamma>0$. This equation is given by $\Delta^2u - \gamma\Delta u=f$ under clamped boundary conditions. Here we show, that given a positive $f$, i.e. upwards pushing, we find a $\gamma_0>0$ such that for all $\gamma\geq \gamma_0$ the bending $u$ is indeed positive. In contrast to a recent result by Cassani\&Tarsia, our approach is valid in all dimensions.


[69] 2106.09342

On the Transcendence of Period Images

Let $f : X \to S$ be a family of smooth projective algebraic varieties over a smooth connected base $S$, with everything defined over $\overline{\mathbb{Q}}$. Denote by $\mathbb{V} = R^{2i} f_{*} \mathbb{Z}(i)$ the associated integral variation of Hodge structure on the degree $2i$ cohomology. We consider the following question: when can a fibre $\mathbb{V}_{s}$ above an algebraic point $s \in S(\overline{\mathbb{Q}})$ be isomorphic to a transcendental fibre $\mathbb{V}_{s'}$ with $s' \in S(\mathbb{C}) \setminus S(\overline{\mathbb{Q}})$? When $\mathbb{V}$ induces a quasi-finite period map $\varphi : S \to \Gamma \backslash D$, conjectures in Hodge theory predict that such isomorphisms cannot exist. We introduce new differential-algebraic techniques to show this is true for all points $s \in S(\overline{\mathbb{Q}})$ outside of an explicit proper closed algebraic subset of $S$. As a corollary we establish the existence of a canonical $\overline{\mathbb{Q}}$-algebraic model for normalizations of period images.


[70] 2106.09345

Scaling limits of planar symplectic ensembles

We consider various asymptotic scaling limits $N\to\infty$ for the $2N$ complex eigenvalues of non-Hermitian random matrices in the symmetry class of the symplectic Ginibre ensemble. These are known to be integrable, forming Pfaffian point processes, and we obtain limiting expressions for the corresponding kernel for different potentials. The first part is devoted to the symplectic Ginibre ensemble with a Gaussian potential. We obtain the asymptotic at the edge of the spectrum in the vicinity of the real line. The unifying form of the kernel allows us to make contact with the bulk scaling along the real line and with the edge scaling away from the real line, where we recover the known determinantal process of the complex Ginibre ensemble. Part two covers ensembles of Mittag-Leffler type with a singularity at the origin. For potentials $Q(\zeta)=|\zeta|^{2\lambda}-(2c/N)\log|\zeta|$, with $\lambda>0$ and $c>-1$, the limiting kernel obeys a linear differential equation of fractional order $1/\lambda$ at the origin. For integer $m=1/\lambda$ it can be solved in terms of Mittag-Leffler functions. In the last part, we derive the Ward's equation for a general class of potentials as a tool to investigate universality. This allows us to determine the functional form of kernels that are translation invariant up to its integration domain.


[71] 2106.09348

Hybrid high-order methods. A primer with application to solid mechanics

This book is organized into eight chapters. The first three gently introduce the basic principles of hybrid high-order methods on a linear diffusion problem, the key ideas underlying the mathematical analysis, and some useful variants of the method as well as links to other methods from the literature. The following four present various challenging applications to solid mechanics, including linear elasticity and hyperelasticity, elastodynamics, contact/friction, and plasticity. The last chapter reviews implementation aspects. This book is primarily intended for graduate students, researchers (in applied mathematics, numerical analysis, and computational mechanics), and engineers working in related fields of application. Basic knowledge of the devising and analysis of finite element methods is assumed. Special effort was made to streamline the presentation so as to pinpoint the essential ideas, address key mathematical aspects, present examples, and provide bibliographic pointers.


[72] 2106.09360

A recursive Lovász theta number for simplex-avoiding sets

We recursively extend the Lov\'asz theta number to geometric hypergraphs on the unit sphere and on Euclidean space, obtaining an upper bound for the independence ratio of these hypergraphs. As an application we reprove a result in Euclidean Ramsey theory in the measurable setting, namely that every $k$-simplex is exponentially Ramsey, and we improve existing bounds for the base of the exponential.


[73] 2106.09364

Wigner transform and quasicrystals

Quasicrystals are tempered distributions $\mu$ which satisfy symmetric conditions on $\mu$ and $\widehat \mu$. This suggests that techniques from time-frequency analysis could possibly be useful tools in the study of such structures. In this paper we explore this direction considering quasicrystals type conditions on time-frequency representations instead of separately on the distribution and its Fourier transform. More precisely we prove that a tempered distribution $\mu$ on ${\mathbb R}^d$ whose Wigner transform, $W(\mu)$, is supported on a product of two uniformly discrete sets in ${\mathbb R}^d$ is a quasicrystal. This result is partially extended to a generalization of the Wigner transform, called matrix-Wigner transform which is defined in terms of the Wigner transform and a linear map $T$ on ${\mathbb R}^{2d}$.


[74] 2106.09366

Boundary path groupoids of generalized Boolean dynamical systems and their C*-algebras

In this paper, we provide two types of boundary path groupoids from a generalized Boolean dynamical system $(\mathcal{B},\mathcal{L}, \theta, \mathcal{I}_{\alpha})$. For the first groupoid, we associate an inverse semigroup to a generalized Boolean dynamical system and use the tight spectrum $\mathsf{T}$ as the unit space of a groupoid $\Gamma(\mathcal{B},\mathcal{L}, \theta, \mathcal{I}_{\alpha})$ that is isomorphic to the tight groupoid $\mathcal{G}_{tight}$. The other one is defined as the Renault-Deaconu groupoid $\Gamma(\partial E, \sigma_E)$ arising from a topological correspondence $E$ associated with a generalized Boolean dynamical system. We then prove that the tight spectrum $\mathsf{T} $ is homeomorphic to the boundary path space $\partial E$ obtained from the topological correspondence. Using this result, we prove that the groupoid $\Gamma(\mathcal{B},\mathcal{L}, \theta, \mathcal{I}_{\alpha})$ equipped with the topology induced from the topology on $\mathcal{G}_{tight}$ is isomorphic to $\Gamma(\partial E, \sigma_E)$ as a topological groupoid. Finally, we show that their $C^*$-algebras are isomorphic to the $C^*$-algebra of the generalized Boolean dynamical system.


[75] 2106.09387

Taming Nonconvexity in Kernel Feature Selection---Favorable Properties of the Laplace Kernel

Kernel-based feature selection is an important tool in nonparametric statistics. Despite many practical applications of kernel-based feature selection, there is little statistical theory available to support the method. A core challenge is the objective function of the optimization problems used to define kernel-based feature selection are nonconvex. The literature has only studied the statistical properties of the \emph{global optima}, which is a mismatch, given that the gradient-based algorithms available for nonconvex optimization are only able to guarantee convergence to local minima. Studying the full landscape associated with kernel-based methods, we show that feature selection objectives using the Laplace kernel (and other $\ell_1$ kernels) come with statistical guarantees that other kernels, including the ubiquitous Gaussian kernel (or other $\ell_2$ kernels) do not possess. Based on a sharp characterization of the gradient of the objective function, we show that $\ell_1$ kernels eliminate unfavorable stationary points that appear when using an $\ell_2$ kernel. Armed with this insight, we establish statistical guarantees for $\ell_1$ kernel-based feature selection which do not require reaching the global minima. In particular, we establish model-selection consistency of $\ell_1$-kernel-based feature selection in recovering main effects and hierarchical interactions in the nonparametric setting with $n \sim \log p$ samples.


[76] 2106.09391

Convergence of Dynamic Programming on the Semidefinite Cone

The goal of this paper is to investigate new and simple convergence analysis of dynamic programming for linear quadratic regulator problem of discrete-time linear time-invariant systems. In particular, bounds on errors are given in terms of both matrix inequalities and matrix norm. Under a mild assumption on the initial parameter, we prove that the Q-value iteration exponentially converges to the optimal solution. Moreover, a global asymptotic convergence is also presented. These results are then extended to the policy iteration. We prove that in contrast to the Q-value iteration, the policy iteration always converges exponentially fast. An example is given to illustrate the results.


[77] 2106.09394

On temporal homogenization in the numerical simulation of atherosclerotic plaque growth

A temporal homogenization approach for the numerical simulation of atherosclerotic plaque growth is extended to fully coupled fluid-structure interaction (FSI) simulations. The numerical results indicate that the two-scale approach yields significantly different results compared to a simple heuristic averaging, where only stationary long-scale FSI problems are solved, confirming the importance of incorporating stress variations on small time-scales. In the homogenization approach, a periodic fine-scale problem, which is periodic with respect to the heart beat, has to be solved for each long-scale time step. Even if no exact initial conditions are available, periodicity can be achieved within only 2-3 heart beats by simple time-stepping.


[78] 2106.09397

Quantized Federated Learning under Transmission Delay and Outage Constraints

Federated learning (FL) has been recognized as a viable distributed learning paradigm which trains a machine learning model collaboratively with massive mobile devices in the wireless edge while protecting user privacy. Although various communication schemes have been proposed to expedite the FL process, most of them have assumed ideal wireless channels which provide reliable and lossless communication links between the server and mobile clients. Unfortunately, in practical systems with limited radio resources such as constraint on the training latency and constraints on the transmission power and bandwidth, transmission of a large number of model parameters inevitably suffers from quantization errors (QE) and transmission outage (TO). In this paper, we consider such non-ideal wireless channels, and carry out the first analysis showing that the FL convergence can be severely jeopardized by TO and QE, but intriguingly can be alleviated if the clients have uniform outage probabilities. These insightful results motivate us to propose a robust FL scheme, named FedTOE, which performs joint allocation of wireless resources and quantization bits across the clients to minimize the QE while making the clients have the same TO probability. Extensive experimental results are presented to show the superior performance of FedTOE for a deep learning-based classification task with transmission latency constraints.


[79] 2106.09401

Asymptotic normality for $m$-dependent and constrained $U$-statistics, with applications to pattern matching in random strings and permutations

We study (asymmetric) $U$-statistics based on a stationary sequence of $m$-dependent variables; moreover, we consider constrained $U$-statistics, where the defining multiple sum only includes terms satisfying some restrictions on the gaps between indices. Results include a law of large numbers and a central limit theorem. Special attention is paid to degenerate cases where, after the standard normalization, the asymptotic variance vanishes; in these cases non-normal limits occur after a different normalization. The results are motivated by applications to pattern matching in random strings and permutations. We obtain both new results and new proofs of old results.


[80] 2106.09403

Density of Free Modules over Finite Chain Rings

In this paper we focus on modules over a finite chain ring $\mathcal{R}$ of size $q^s$. We compute the density of free modules of $\mathcal{R}^n$, where we separately treat the asymptotics in $n,q$ and $s$. In particular, we focus on two cases: one where we fix the type of the module and one where we fix the rank of the module. In both cases, the density results can be bounded by the Andrews-Gordon identities. We also study the asymptotic behaviour of modules generated by random matrices over $\mathcal{R}$. Since linear codes over $\mathcal{R}$ are submodules of $\mathcal{R}^n$ we get direct implications for coding theory. For example, we show that random codes achieve the Gilbert-Varshamov bound with high probability.


[81] 2106.09404

The tiny trace ideals of the canonical modules in Cohen-Macaulay rings of dimension one

We study one-dimensional Cohen-Macaulay rings whose trace ideal of the canonical module is as small as possible. In this paper we call such rings far-flung Gorenstein rings. We investigate far-flung Gorenstein rings in relation with the endomorphism algebras of the maximal ideals and numerical semigroup rings. We show that the solution of the Rohrbach problem in additive number theory provides an upper bound for the multiplicity of far-flung Gorenstein numerical semigroup rings. Reflexive modules over far-flung Gorenstein rings are also studied.


[82] 2106.09405

Mertens conjectures in absorbing games with incomplete information

In a zero-sum stochastic game with signals, at each stage, two adversary players take decisions and receive a stage payoff determined by these decisions and a variable called state. The state follows a Markov chain, that is controlled by both players. Actions and states are imperfectly observed by players, who receive a private signal at each stage. Mertens (ICM 1986) conjectured two properties regarding games with long duration: first, that limit value always exists, second, that when Player 1 is more informed than Player 2, she can guarantee uniformly the limit value. These conjectures were disproved recently by the author, but remain widely open in many subclasses. A well-known particular subclass is the one of absorbing games with incomplete information on both sides, in which the state can move at most once during the game, and players get a private signal about it at the outset of the game. This paper proves Mertens conjectures in this particular model, by introducing a new approximation technique of belief dynamics, that is likely to generalize to many other frameworks. In particular, this makes a significant step towards the understanding of the following broad question: in which games do Mertens conjectures hold?


[83] 2106.09414

Internal structures in extended affine Lie algebras

We construct certain integral structures for the cores of reduced tame extended affine Lie algebras of rank at least 2. One of the main tools to achieve this is a generalization of Chevalley automorphisms in the context of extended affine Lie algebras. As an application, groups of extended affine Lie type associated to the adjoint representation are defined over arbitrary fields.


[84] 2106.09421

State Estimation with Model Reduction and Shape Variability. Application to biomedical problems

We develop a mathematical and numerical framework to solve state estimation problems for applications that present variations in the shape of the spatial domain. This situation arises typically in a biomedical context where inverse problems are posed on certain organs or portions of the body which inevitably involve morphological variations. If one wants to provide fast reconstruction methods, the algorithms must take into account the geometric variability. We develop and analyze a method which allows to take this variability into account without needing any a priori knowledge on a parametrization of the geometrical variations. For this, we rely on morphometric techniques involving Multidimensional Scaling, and couple them with reconstruction algorithms that make use of reduced model spaces pre-computed on a database of geometries. We prove the potential of the method on a synthetic test problem inspired from the reconstruction of blood flows and quantities of medical interest with Doppler ultrasound imaging.


[85] 2106.09425

Partially multiplicative quandles and simplicial Hurwitz spaces

This is the first article in a series about Hurwitz spaces. We introduce the notion of partially multiplicative quandle (PMQ), which is a generalisation of the classical notions of partial abelian monoid and of quandle. We give a first, simplicial definition of Hurwitz spaces of configurations of points in the plane with monodromies in a PMQ.


[86] 2106.09428

Equivalence of cubical and simplicial approaches to $(\infty,n)$-categories

We prove that the marked triangulation functor from the category of marked cubical sets equipped with a model structure for ($n$-trivial, saturated) comical sets to the category of marked simplicial set equipped with a model structure for ($n$-trivial, saturated) complicial sets is a Quillen equivalence. Our proof is based on the theory of cones, previously developed by the first two authors together with Lindsey and Sattler.


[87] 2106.09441

Traveling waves for 2D parabolic Allen-Cahn systems

This paper is devoted to the study of the existence of traveling waves $w: [0,+\infty) \times \mathbb{R}^2 \to \mathbb{R}^k$ ($k \geq 2$) for the parabolic Allen-Cahn system \begin{equation} \partial_t w - \Delta w = -\nabla_u V(w) \mbox{ in } [0,+\infty) \times \mathbb{R}^2, \end{equation} satisfying some conditions at infinity. The potential $V$ is a non-negative and smooth multi-well potential, meaning its null set is finite and contains at least two elements. We ask $w$ to propagate in one direction with positive speed $c$ and with profile $\mathfrak{U}$. The pair $(c,\mathfrak{U})$ is then expected to solve the resulting entire elliptic system, which is supplemented with some conditions at infinity for $\mathfrak{U}$. In particular, we ask it to join at infinity (in a suitable sense) two locally minimizing 1D heteroclinics which have different energy. The $c=0$ case (stationary problem) yields the well-known solution of Alama, Bronsard and Gui (1997) in the situation of two distinct globally minimizing 1D heteroclinics. However, to our knowledge, the non-stationary problem has not been addressed before. In this paper, we show the existence of such a pair $(c,\mathfrak{U})$ under the suitable assumptions. In particular, our proof is limited to potentials $V$ which are a suitable small perturbation of those considered for the stationary problem. Our strategy consists on adapting a result of Alikakos and Katzourakis (2011) for curves which take values in a possibly infinite-dimensional Hilbert space and showing that it applies to our initial setting.


[88] 2106.09442

Dynamic Metasurface Antennas for Energy Efficient Massive MIMO Uplink Communications

Future wireless communications are largely inclined to deploy a massive number of antennas at the base stations (BS) by exploiting energy-efficient and environmentally friendly technologies. An emerging technology called dynamic metasurface antennas (DMAs) is promising to realize such massive antenna arrays with reduced physical size, hardware cost, and power consumption. This paper aims to optimize the energy efficiency (EE) performance of DMAs-assisted massive MIMO uplink communications. We propose an algorithmic framework for designing the transmit precoding of each multi-antenna user and the DMAs tuning strategy at the BS to maximize the EE performance, considering the availability of the instantaneous and statistical channel state information (CSI), respectively. Specifically, the proposed framework includes Dinkelbach's transform, alternating optimization, and deterministic equivalent methods. In addition, we obtain a closed-form solution to the optimal transmit signal directions for the statistical CSI case, which simplifies the corresponding transmission design. The numerical results show good convergence performance of our proposed algorithms as well as considerable EE performance gains of the DMAs-assisted massive MIMO uplink communications over the baseline schemes.


[89] 2106.09445

A structure-preserving surrogate model for the closure of the moment system of the Boltzmann equation using convex deep neural networks

Direct simulation of physical processes on a kinetic level is prohibitively expensive in aerospace applications due to the extremely high dimension of the solution spaces. In this paper, we consider the moment system of the Boltzmann equation, which projects the kinetic physics onto the hydrodynamic scale. The unclosed moment system can be solved in conjunction with the entropy closure strategy. Using an entropy closure provides structural benefits to the physical system of partial differential equations. Usually computing such closure of the system spends the majority of the total computational cost, since one needs to solve an ill-conditioned constrained optimization problem. Therefore, we build a neural network surrogate model to close the moment system, which preserves the structural properties of the system by design, but reduces the computational cost significantly. Numerical experiments are conducted to illustrate the performance of the current method in comparison to the traditional closure.


[90] 2106.09448

Minimizing under relaxed symmetry constraints: Triple and $N$-junctions

We consider a nonnegative potential $W:\mathbb{R}^2\rightarrow\mathbb{R}$ invariant under the action of the rotation group $C_N$ of the regular polygon with $N$ sides, $N\geq 3$. We assume that $W$ has $N$ nondegenerate zeros and prove the existence of a $N$-junction solution to the vector Allen-Cahn equation. The proof is variational and is based on sharp lower and upper bounds for the energy and on a new pointwise estimate for vector minimizers.


[91] 2106.09450

Simultaneous Transmission and Reflection Reconfigurable Intelligent Surface Assisted MIMO Systems

In this work, we investigate a novel simultaneous transmission and reflection reconfigurable intelligent surface (RIS)-assisted multiple-input multiple-output downlink system, where three practical transmission protocols, namely, energy splitting (ES), mode selection (MS), and time splitting (TS), are studied. For the system under consideration, we maximize the weighted sum rate with multiple coupled variables. To solve this optimization problem, a block coordinate descent algorithm is proposed to reformulate this problem and design the precoding matrices and the transmitting and reflecting coefficients (TARCs) in an alternate manner. Specifically, for the ES scheme, the precoding matrices are solved using the Lagrange dual method, while the TARCs are obtained using the penalty concave-convex method. Additionally, the proposed method is extended to the MS scheme by solving a mixed-integer problem. Moreover, we solve the formulated problem for the TS scheme using a one-dimensional search and the Majorization-Minimization technique. Our simulation results reveal that: 1) Simultaneous transmission and reflection RIS (STAR-RIS) can achieve better performance than reflecting-only RIS; 2) In unicast communication, TS scheme outperforms the ES and MS schemes, while in broadcast communication, ES scheme outperforms the TS and MS schemes.


[92] 2106.09452

Spectral convergence of high-dimensional spheres to Gaussian spaces

We prove that the spectral structure on the $N$-dimensional standard sphere of radius $(N-1)^{1/2}$ compatible with a projection onto the first $n$-coordinates converges to the spectral structure on the $n$-dimensional Gaussian space with variance $1$ as $N\to \infty$. We also show the analogue for the first Dirichlet eigenvalue and its eigenfunction on a ball in the sphere and on a half-space in the Gaussian space.


[93] 2106.09466

A Rigorous Derivation of the Functional Renormalisation Group Equation

The functional renormalisation group equation is derived in a mathematically rigorous fashion in a framework suitable for the Osterwalder-Schrader formulation of quantum field theory. To this end, we devise a very general regularisation scheme and give precise conditions for the involved regulators guaranteeing physical boundary conditions. Furthermore, it is shown how the classical limit is altered by the regularisation process leading to an inevitable breaking of translation invariance. We also give precise conditions for the convergence of the obtained theories upon removal of the regularisation.


[94] 2106.09468

Vertex-regular $1$-factorizations in infinite graphs

The existence of $1$-factorizations of an infinite complete equipartite graph $K_m[n]$ (with $m$ parts of size $n$) admitting a vertex-regular automorphism group $G$ is known only when $n=1$ and $m$ is countable (that is, for countable complete graphs) and, in addition, $G$ is a finitely generated abelian group $G$ of order $m$. In this paper, we show that a vertex-regular $1$-factorization of $K_m[n]$ under the group $G$ exists if and only if $G$ has a subgroup $H$ of order $n$ whose index in $G$ is $m$. Furthermore, we provide a sufficient condition for an infinite Cayley graph to have a regular $1$-factorization. Finally, we construct 1-factorizations that contain a given subfactorization, both having a vertex-regular automorphism group.


[95] 2106.09469

A posteriori estimator for the adaptive solution of a quasi-static fracture phase-field model with irreversibility constraints

Within this article, we develop a residual type a posteriori error estimator for a time discrete quasi-static phase-field fracture model. Particular emphasize is given to the robustness of the error estimator for the variational inequality governing the phase-field evolution with respect to the phase-field regularization parameter $\epsilon$. The article concludes with numerical examples highlighting the performance of the proposed a posteriori error estimators on three standard test cases; the single edge notched tension and shear test as well as the L-shaped panel test.


[96] 2106.09471

A skeleton model to enumerate standard puzzle sequences

Guo-Niu Han [arXiv:2006.14070 [math.CO]] has introduced a new combinatorial object named standard puzzle. We use digraphs to show the relations between numbers in standard puzzles and propose a skeleton model. By this model, we solve the enumeration problem of over fifty thousand standard puzzle sequences. Most of them can be represented by classical numbers, such as Catalan numbers, double factorials, Secant numbers and so on. Also, we prove several identities in standard puzzle sequences.


[97] 2106.09481

Stochastic Bias-Reduced Gradient Methods

We develop a new primitive for stochastic optimization: a low-bias, low-cost estimator of the minimizer $x_\star$ of any Lipschitz strongly-convex function. In particular, we use a multilevel Monte-Carlo approach due to Blanchet and Glynn to turn any optimal stochastic gradient method into an estimator of $x_\star$ with bias $\delta$, variance $O(\log(1/\delta))$, and an expected sampling cost of $O(\log(1/\delta))$ stochastic gradient evaluations. As an immediate consequence, we obtain cheap and nearly unbiased gradient estimators for the Moreau-Yoshida envelope of any Lipschitz convex function, allowing us to perform dimension-free randomized smoothing. We demonstrate the potential of our estimator through four applications. First, we develop a method for minimizing the maximum of $N$ functions, improving on recent results and matching a lower bound up logarithmic factors. Second and third, we recover state-of-the-art rates for projection-efficient and gradient-efficient optimization using simple algorithms with a transparent analysis. Finally, we show that an improved version of our estimator would yield a nearly linear-time, optimal-utility, differentially-private non-smooth stochastic optimization method.


[98] 2106.09485

Secure Multi-Function Computation with Private Remote Sources

We consider a distributed function computation problem in which parties observing noisy versions of a remote source facilitate the computation of a function of their observations at a fusion center through public communication. The distributed function computation is subject to constraints, including not only reliability and storage but also privacy and secrecy. Specifically, 1) the remote source should remain private from an eavesdropper and the fusion center, measured in terms of the information leaked about the remote source; 2) the function computed should remain secret from the eavesdropper, measured in terms of the information leaked about the arguments of the function, to ensure secrecy regardless of the exact function used. We derive the exact rate regions for lossless and lossy single-function computation and illustrate the lossy single-function computation rate region for an information bottleneck example, in which the optimal auxiliary random variables are characterized for binary-input symmetric-output channels. We extend the approach to lossless and lossy asynchronous multiple-function computations with joint secrecy and privacy constraints, in which case inner and outer bounds for the rate regions differing only in the Markov chain conditions imposed are characterized.


[99] 2106.09497

On ergodic control problem for viscous Hamilton--Jacobi equations for weakly coupled elliptic systems

In this article we study ergodic problems in the whole space $\mathbb{R}^N$ for a weakly coupled systems of viscous Hamilton-Jacobi equations with coercive right-hand sides. The Hamiltonians are assumed to have a fairly general structure and the switching rates need not be constant. We prove the existence of a critical value $\lambda^*$ such that the ergodic eigenvalue problem has a solution for every $\lambda\leq\lambda^*$ and no solution for $\lambda>\lambda^*$. Moreover, the existence and uniqueness of non-negative solutions corresponding to the value $\lambda^*$ are also established. We also exhibit the implication of these results to the ergodic optimal control problems of controlled switching diffusions.


[100] 2106.09507

On endomorphism algebras of Gelfand-Graev representations

For a connected reductive group $G$ defined over $\mathbb{F}_q$ and equipped with the induced Frobenius endomorphism $F$, we study the relation among the following three $\mathbb{Z}$-algebras: (i) the $\mathbb{Z}$-model $\mathsf{E}_G$ of endomorphism algebras of Gelfand-Graev representations of $G^F$; (ii) the Grothendieck group $\mathsf{K}_{G^\ast}$ of the category of representations of $G^{\ast F^\ast}$ over $\overline{\mathbb{F}_q}$ (Deligne-Lusztig dual side); (iii) the ring $\mathsf{B}_{G^\vee}$ of the scheme $(T^\vee/\!\!/ W)^{F^\vee}$ over $\mathbb{Z}$ (Langlands dual side). The comparison between (i) and (iii) is motivated by recent advances in the local Langlands program.


[101] 2106.09510

Extremal mild solutions for Hilfer fractional evolution equation with mixed monotone Impulsive conditions

The well established mixed monotone iterative technique that is used to study the existence and uniqueness of fractional order system is studied explicitly for impulsive system with Hilfer fractional order in this paper. The procedure of finding mild $L$-quasi solution of such impulsive evolution equation with noncomapct semigroups involves measure of non-compactness and Sadovskii's fixed point theorem as well. An example is provided to illustrate the main results.


[102] 2106.09511

Gevrey well posedness for $3$-evolution equations with variable coefficients

We study the Cauchy problem for a class of third order linear anisotropic evolution equations with complex valued lower order terms depending both on time and space variables. Under suitable decay assumptions for $|x| \to \infty$ on these coefficients, we prove a well posedness result in Gevrey-type spaces.


[103] 2106.09519

The Quasi-Zariski topology on the graded quasi-primary spectrum of a graded module over a graded commutative ring

Let $G$ be a group with identity $e$. Let $R$ be a $G$-graded commutative ring and $M$ a graded $R$-module. A proper graded submodule $Q$ of $M$ is called a graded quasi-primary submodule if whenever $r\in h(R)$ and $m\in h(M)$ with $rm\in Q$, then either $r\in Gr((Q:_{R}M))$ or $m\in Gr_{M}(Q)$. The graded quasi primary spectrum $qp.Spec_{g}(M)$ is defined to be the set of all graded quasi primary submodules of $M$. In this paper, we introduce and study a topology on $qp.Spec_{g}(M)$, called the Quasi-Zariski Topology, and investigate properties of this topology and some conditions under which (% $qp.Spec_{g}(M)$, $q.\tau ^{g}$) is a Noetherian, spctral space.


[104] 2106.09528

Optimal control strategies to tailor antivirals for acute infectious diseases in the host

Several mathematical models in SARS-CoV-2 have shown how target-cell model can help to understand the spread of the virus in the host and how potential candidates of antiviral treatments can help to control the virus. Concepts as equilibrium and stability show to be crucial to qualitatively determine the best alternatives to schedule drugs, according to effectivity in inhibiting the virus infection and replication rates. Important biological events such as rebounds of the infections (when antivirals are incorrectly interrupted) can also be explained by means of a dynamic study of the target-cell model. In this work, a full characterization of the dynamical behavior of the target-cell models under control actions is made and, based on this characterization, the optimal fixed-dose antiviral schedule that produces the smallest amount of dead cells (without viral load rebounds) is computed. Several simulation results - performed by considering real patient data - show the potential benefits of both, the model characterization and the control strategy.


[105] 2106.09549

Optimal thresholds for preserving embeddedness of elastic flows

We consider elastic flows of closed curves in Euclidean space. We obtain optimal energy thresholds below which elastic flows preserve embeddedness of initial curves for all time. The obtained thresholds take different values between codimension one and higher. The main novelty lies in the case of codimension one, where we obtain the variational characterization that the thresholding shape is a minimizer of the bending energy (normalized by length) among all nonembedded planar closed curves of unit rotation number. It turns out that a minimizer is uniquely given by a nonclassical shape, which we call ``elastic two-teardrop''.


[106] 2106.09552

Mixing of the Averaging process and its discrete dual on finite-dimensional geometries

We analyze the $L^1$-mixing of a generalization of the Averaging process introduced by Aldous. The process takes place on a growing sequence of graphs which we assume to be finite-dimensional, in the sense that the random walk on those geometries satisfies a family of Nash inequalities. As a byproduct of our analysis, we provide a complete picture of the total variation mixing of a discrete dual of the Averaging process, which we call Binomial Splitting process. A single particle of this process is essentially the random walk on the underlying graph. When several particles evolve together, they interact by synchronizing their jumps when placed on neighboring sites. We show that, given $k$ the number of particles and $n$ the (growing) size of the underlying graph, the system exhibits cutoff in total variation if $k\to\infty$ and $k=O(n^2)$. Finally, we exploit the duality between the two processes to show that the Binomial Splitting satisfies a version of Aldous' spectral gap identity, namely, the relaxation time of the process is independent of the number of particles.


[107] 2106.09554

A Unique Perfect Power Decagonal Number

Let $\mathcal{P}_s(n)$ denotes the $n$th $s$-gonal number. We consider the equation \[\mathcal{P}_s(n) = y^m, \] for integers $n,s,y,$ and $m$. All solutions to this equation are known for $m>2$ and $s \in \{3,5,6,8,20 \}$. We consider the case $s=10$, that of decagonal numbers. Using a descent argument and the modular method, we prove that the only decagonal number $>1$ expressible as a perfect $m$th power with $m>1$ is $\mathcal{P}_{10}(3) = 3^3$.


[108] 2106.09561

Cayley (Di)Hypergraphs

A Cayley (di)hypergraph is a hypergraph that its automorphism group contains a subgroup acting regularly on (hyper)vertices. In this paper, we study Cayley (di)hypergraph and its automorphism group.


[109] 2106.09569

Basic PU(1,1)-representations of the hyperelliptic group are discrete

We show that a PU(1,1)-representation of the hyperelliptic group $H_{n}$ is basic if and only if it is discrete and faithful, thus partially proving a conjecture by S. Anan'in and E. Bento Gon\c{c}alves in the case of the Poincar\'{e} disc.


[110] 2106.09570

Noise sensitivity for the top eigenvector of a sparse random matrix

We investigate the noise sensitivity of the top eigenvector of a sparse random symmetric matrix. Let $v$ be the top eigenvector of an $N\times N$ sparse random symmetric matrix with an average of $d$ non-zero centered entries per row. We resample $k$ randomly chosen entries of the matrix and obtain another realization of the random matrix with top eigenvector $v^{[k]}$. Building on recent results on sparse random matrices and a noise sensitivity analysis previously developed for Wigner matrices, we prove that, if $d\geq N^{2/9}$, with high probability, when $k \ll N^{5/3}$, the vectors $v$ and $v^{[k]}$ are almost collinear and, on the contrary, when $k\gg N^{5/3}$, the vectors $v$ and $v^{[k]}$ are almost orthogonal. A similar result holds for the eigenvector associated to the second largest eigenvalue of the adjacency matrix of an Erd\H{o}s-R\'enyi random graph with average degree $d \geq N^{2/9}$.


[111] 2106.09582

Improvement of generalization of Larman-Rogers-Seidel's theorem

A finite set $X$ in the $d$-dimensional Euclidean space is called an $s$-distance set if the set of distances between any two distinct points of $X$ has size $s$. In 1977, Larman-Rogers-Seidel proved that if the cardinality of an two-distance set is large enough, then there exists an integer $k$ such that the two distances $\alpha$, $\beta$ $(\alpha < \beta)$ having the integer condition, namely, $\frac{\alpha^2}{\beta^2}=\frac{k-1}{k}$. In 2011, Nozaki generalized Larman-Rogers-Seidel's theorem to the case of $s$-distance sets, i.e. if the cardinality of an $s$-distance set $|X|\geqslant 2N$ with distances $\alpha_1,\alpha_2,\cdots,\alpha_s$, where $N=\binom{d+s-1}{s-1}+\binom{d+s-2}{s-2}$, then the numbers $k_i=\prod_{j=1,2,\cdots,s,\text{ }j\neq i}\frac{\alpha_{j}^{2}}{\alpha_{j}^{2}-\alpha_{i}^{2}}$ are integers. In this note, we reduce the lower bound of the requirement of integer condition of $s$-distance sets in $\mathbb{R}^d$. Furthermore, we can show that there are only finitely many $s$-distance sets $X$ in $\mathbb{R}^d$ with $|X|\geqslant 2\binom{d+s-1}{s-1}.$


[112] 2106.09585

Another Riemann Hypothesis Equivalent

Identities involving Mobius function values (u(j),u(k)) are used to generate a Riemann Hypothesis equivalent.


[113] 2106.09591

Notes on regularity of Anosov splitting

In these expository notes, we give a proof of regularity of Anosov splitting for Anosov diffeomorphisms in a torus. We also generalize the idea to higher dimensions and to Anosov flows.


[114] 2106.09605

Asymptotic stability of the sine-Gordon kink under odd perturbations

We establish the asymptotic stability of the sine-Gordon kink under odd perturbations that are sufficiently small in a weighted Sobolev norm. Our approach is perturbative and does not rely on the complete integrability of the sine-Gordon model. Key elements of our proof are a specific factorization property of the linearized operator around the sine-Gordon kink, a remarkable non-resonance property exhibited by the quadratic nonlinearity in the Klein-Gordon equation for the perturbation, and a variable coefficient quadratic normal form introduced in [52]. We emphasize that the restriction to odd perturbations does not bypass the effects of the odd threshold resonance of the linearized operator. Our techniques have applications to soliton stability questions for several well-known non-integrable models, for instance, to the asymptotic stability problem for the kink of the $\phi^4$ model as well as to the conditional asymptotic stability problem for the solitons of the focusing quadratic and cubic Klein-Gordon equations in one space dimension.


[115] 2106.09606

Cardinality Minimization, Constraints, and Regularization: A Survey

We survey optimization problems that involve the cardinality of variable vectors in constraints or the objective function. We provide a unified viewpoint on the general problem classes and models, and give concrete examples from diverse application fields such as signal and image processing, portfolio selection, or machine learning. The paper discusses general-purpose modeling techniques and broadly applicable as well as problem-specific exact and heuristic solution approaches. While our perspective is that of mathematical optimization, a main goal of this work is to reach out to and build bridges between the different communities in which cardinality optimization problems are frequently encountered. In particular, we highlight that modern mixed-integer programming, which is often regarded as impractical due to commonly unsatisfactory behavior of black-box solvers applied to generic problem formulations, can in fact produce provably high-quality or even optimal solutions for cardinality optimization problems, even in large-scale real-world settings. Achieving such performance typically draws on the merits of problem-specific knowledge that may stem from different fields of application and, e.g., shed light on structural properties of a model or its solutions, or lead to the development of efficient heuristics; we also provide some illustrative examples.


[116] 2106.09615

On multiplicative Chung--Diaconis--Graham process

We study the lazy Markov chain on $\mathbf{F}_p$ defined as $X_{n+1}=X_n$ with probability $1/2$ and $X_{n+1}=f(X_n) \cdot \varepsilon_{n+1}$, where $\varepsilon_n$ are random variables distributed uniformly on $\{ \gamma^{}, \gamma^{-1}\}$, $\gamma$ is a primitive root and $f(x) = \frac{x}{x-1}$ or $f(x)=\mathrm{ind} (x)$. Then we show that the mixing time of $X_n$ is $\exp(O(\log p / \log \log p))$. Also, we obtain an application to an additive--combinatorial question concerning a certain Sidon--type family of sets.


[117] 2106.09617

On Tutte cycles containing three prescribed edges

A cycle $C$ in a graph $G$ is called a Tutte cycle if, after deleting $C$ from $G$, each component has at most three neighbors on $C$. Tutte cycles play an important role in the study of Hamiltonicity of planar graphs. Thomas and Yu and independently Sanders proved the existence of Tutte cycles containining three specified edges of a facial cycle in a 2-connected plane graph. We prove a quantitative version of this result, bounding the number of components of the graph obtained by deleting a Tutte cycle. As a corollary, we can find long cycles in essentially 4-connected plane graphs that also contain three prescribed edges of a facial cycle.


[118] 2106.09619

Asymptotic bounds for cycle integrals of the j-function on Markov geodesics

We give asymptotic upper and lower bounds for the real and imaginary parts of cycle integrals of the classical modular j-function along geodesics that correspond to Markov irrationalities.


[119] 2106.09625

Poisson-Dirichlet asymptotics in condensing particle systems

We study measures on random partitions, arising from condensing stochastic particle systems with stationary product distributions. We provide fairly general conditions on the stationary weights, which lead to Poisson-Dirichlet statistics of the condensed phase in the thermodynamic limit. The Poisson-Dirichlet distribution is known to be the unique reversible measure of split-merge dynamics for random partitions, which we use to characterize the limit law. We also establish concentration results for the macroscopic phase, using size-biased sampling techniques and the equivalence of ensembles to characterize the bulk distribution of the system.


[120] 2106.09627

Towards Prevention of Sportsmen Burnout: Formal Analysis of Sub-Optimal Tournament Scheduling

Scheduling a sports tournament is a complex optimization problem, which requires a large number of hard constraints to satisfy. Despite the availability of several such constraints in the literature, there remains a gap since most of the new sports events pose their own unique set of requirements, and demand novel constraints. Specifically talking of the strictly time bound events, ensuring fairness between the different teams in terms of their rest days, traveling, and the number of successive games they play, becomes a difficult task to resolve, and demands attention. In this work, we present a similar situation with a recently played sports event, where a suboptimal schedule favored some of the sides more than the others. We introduce various competitive parameters to draw a fairness comparison between the sides and propose a weighting criterion to point out the sides that enjoyed this schedule more than the others. Furthermore, we use root mean squared error between an ideal schedule and the actual ones for each side to determine unfairness in the distribution of rest days across their entire schedules. The latter is crucial, since successively playing a large number of games may lead to sportsmen burnout, which must be prevented.


[121] 2106.09641

Local non-periodic order and diam-mean equicontinuity on cellular automata

Diam-mean equicontinuity is a dynamical property that has been of use in the study of non-periodic order. Using some type of "local" skew product between a shift and an odometer looking cellular automaton (CA) we will show there exists an almost diam-mean equicontinuous CA that is not almost equicontinuous, and hence not almost locally periodic.


[122] 2106.09661

A Fubini-type theorem for Hausdorff dimension

It is well known that a classical Fubini theorem for Hausdorff dimension cannot hold; that is, the dimension of the intersections of a fixed set with a parallel family of planes do not determine the dimension of the set. Here we prove that a Fubini theorem for Hausdorff dimension does hold modulo sets that are small on all Lipschitz curves/surfaces. We say that $G\subset \mathbb{R}^k\times \mathbb{R}^n$ is $\Gamma_k$-null if for every Lipschitz function $f:\mathbb{R}^k\to \mathbb{R}^n$ the set $\{t\in \mathbb{R}^k\,:\,(t,f(t))\in G\}$ has measure zero. We show that for every compact set $E\subset \mathbb{R}^k\times \mathbb{R}^n$ there is a $\Gamma_k$-null subset $G\subset E$ such that $$\dim (E\setminus G) = k+\text{ess-}\sup(\dim E_t)$$ where $\text{ess-}\sup(\dim E_t)$ is the essential supremum of the Hausdorff dimension of the vertical sections $\{E_t\}_{t\in \mathbb{R}^k}$ of $E$, assuming that $proj_{\mathbb{R}^k} E$ has positive measure. We also obtain more general results by replacing $\mathbb{R}^k$ by an Ahlfors regular set. Applications of our results include Fubini-type results for unions of affine subspaces and related projection theorems.


[123] 2106.09663

A Short Note of PAGE: Optimal Convergence Rates for Nonconvex Optimization

In this note, we first recall the nonconvex problem setting and introduce the optimal PAGE algorithm (Li et al., ICML'21). Then we provide a simple and clean convergence analysis of PAGE for achieving optimal convergence rates. Moreover, PAGE and its analysis can be easily adopted and generalized to other works. We hope that this note provides the insights and is helpful for future works.


[124] 2106.09673

Equivariant maps to subshifts whose points have small stabilizers

Let $\Gamma$ be a countably infinite group. Given $k \in \mathbb{N}$, we use $\mathrm{Free}(k^\Gamma)$ to denote the free part of the Bernoulli shift action of $\Gamma$ on $k^\Gamma$. Seward and Tucker-Drob showed that there exists a free subshift $\mathcal{S} \subseteq \mathrm{Free}(2^\Gamma)$ such that every free Borel action of $\Gamma$ on a Polish space admits a Borel $\Gamma$-equivariant map to $\mathcal{S}$. Here we generalize this result as follows. Let $\mathcal{S}$ be a subshift of finite type (for example, $\mathcal{S}$ could be the set of all proper colorings of the Cayley graph of $\Gamma$ with some finite number of colors). Suppose that $\pi \colon \mathrm{Free}(k^\Gamma) \to \mathcal{S}$ is a continuous $\Gamma$-equivariant map and let $\mathrm{Stab}(\pi)$ be the set of all group elements that fix every point in the image of $\pi$. Unless $\pi$ is constant, $\mathrm{Stab}(\pi)$ is a finite normal subgroup of $\Gamma$. We prove that there exists a subshift $\mathcal{S}' \subseteq \mathcal{S}$ such that the stabilizer of every point in $\mathcal{S}'$ is $\mathrm{Stab}(\pi)$ and every free Borel action of $\Gamma$ on a Polish space admits a Borel $\Gamma$-equivariant map to $\mathcal{S}'$.


[125] 2106.09676

On ECP-Groups

According to T. Foguel a subgroup $H$ of a group $G$ is called conjugate-permutable if $ HH^x=H^xH$ for every $x\in G$. Mingyao Xu and Qinhai Zhang studied finite groups with every subgroup conjugate-permutable (ECP-groups) and asked three questions about them. We gave the answers on these questions. In particular, every group of exponent 3 is ECP-group, there exist non-regular ECP-$3$-groups and the class of all finite ECP-groups is neither formation nor variety.


[126] 2106.09682

On the chaoticity of derivatives

We show that the $n$th derivative with maximal domain is a chaotic operator in $L_p(a,b)$ ($1\le p<\infty$, $-\infty<a<b<\infty$) for each $n\in \mathbb{N}$.


[127] 2106.09687

Super-Acceleration with Cyclical Step-sizes

Cyclical step-sizes are becoming increasingly popular in the optimization of deep learning problems. Motivated by recent observations on the spectral gaps of Hessians in machine learning, we show that these step-size schedules offer a simple way to exploit them. More precisely, we develop a convergence rate analysis for quadratic objectives that provides optimal parameters and shows that cyclical learning rates can improve upon traditional lower complexity bounds. We further propose a systematic approach to design optimal first order methods for quadratic minimization with a given spectral structure. Finally, we provide a local convergence rate analysis beyond quadratic minimization for the proposed methods and illustrate our findings through benchmarks on least squares and logistic regression problems.


[128] 2106.09688

A Ramsey-Turán theory for tilings in graphs

For a $k$-vertex graph $F$ and an $n$-vertex graph $G$, an $F$-tiling in $G$ is a collection of vertex-disjoint copies of $F$ in $G$. For $r\in \mathbb{N}$, the $r$-independence number of $G$, denoted $\alpha_r(G)$ is the largest size of a $K_r$-free set of vertices in $G$. In this paper, we discuss Ramsey--Tur\'an-type theorems for tilings where one is interested in minimum degree and independence number conditions (and the interaction between the two) that guarantee the existence of optimal $F$-tilings. For cliques, we show that for any $k\geq 3$ and $\eta>0$, any graph $G$ on $n$ vertices with $\delta(G)\geq \eta n$ and $\alpha_k(G)=o(n)$ has a $K_k$-tiling covering all but $\lfloor\tfrac{1}{\eta}\rfloor(k-1)$ vertices. All conditions in this result are tight; the number of vertices left uncovered can not be improved and for $\eta<\tfrac{1}{k}$, a condition of $\alpha_{k-1}(G)=o(n)$ would not suffice. When $\eta>\tfrac{1}{k}$, we then show that $\alpha_{k-1}(G)=o(n)$ does suffice, but not $\alpha_{k-2}(G)=o(n)$. These results unify and generalise previous results of Balogh-Molla-Sharifzadeh, Nenadov-Pehova and Balogh-McDowell-Molla-Mycroft on the subject. We further explore the picture when $F$ is a tree or a cycle and discuss the effect of replacing the independence number condition with $\alpha^*(G)=o(n)$ (meaning that any pair of disjoint linear sized sets induce an edge between them) where one can force perfect $F$-tilings covering all the vertices. Finally we discuss the consequences of these results in the randomly perturbed setting.


[129] 2106.09691

Interactive Change Point Detection using optimisation approach and Bayesian statistics applied to real world applications

Change point detection becomes more and more important as datasets increase in size, where unsupervised detection algorithms can help users process data. To detect change points, a number of unsupervised algorithms have been developed which are based on different principles. One approach is to define an optimisation problem and minimise a cost function along with a penalty function. In the optimisation approach, the choice of the cost function affects the predictions made by the algorithm. In extension to the existing studies, a new type of cost function using Tikhonov regularisation is introduced. Another approach uses Bayesian statistics to calculate the posterior probability distribution of a specific point being a change point. It uses a priori knowledge on the distance between consecutive change points and a likelihood function with information about the segments. The optimisation and Bayesian approaches for offline change point detection are studied and applied to simulated datasets as well as a real world multi-phase dataset. The approaches have previously been studied separately and a novelty lies in comparing the predictions made by the two approaches in a specific setting, consisting of simulated datasets and a real world example. The study has found that the performance of the change point detection algorithms are affected by the features in the data.


[130] 2106.09697

The non-resonant bilinear Hilbert--Carleson operator

In this paper we introduce the class of bilinear Hilbert--Carleson operators $\{BC^a\}_{a>0}$ defined by $$ BC^{a}(f,g)(x):= \sup_{\lambda\in {\mathbb R}} \Big|\int f(x-t)\, g(x+t)\, e^{i\lambda t^a} \, \frac{dt}{t} \Big| $$ and show that in the non-resonant case $a\in (0,\infty)\setminus\{1,2\}$ the operator $BC^a$ extends continuously from $L^p({\mathbb R})\times L^q({\mathbb R})$ into $L^r({\mathbb R})$ whenever $\frac{1}{p}+\frac{1}{q}=\frac{1}{r}$ with $1<p,\,q\leq\infty$ and $\frac{2}{3}<r<\infty$. A key novel feature of these operators is that -- in the non-resonant case -- $BC^{a}$ has a \emph{hybrid} nature enjoying both (1) ``zero curvature'' features inherited from the modulation invariance property of the classical bilinear Hilbert transform (BHT), and (2) ``non-zero curvature'' features arising from the Carleson-type operator with nonlinear phase $\lambda t^a$.


[131] 2106.09709

Independent sets of a given size and structure in the hypercube

We determine the asymptotics of the number of independent sets of size $\lfloor \beta 2^{d-1} \rfloor$ in the discrete hypercube $Q_d = \{0,1\}^d$ for any fixed $\beta \in [0,1]$ as $d \to \infty$, extending a result of Galvin for $\beta \in [1-1/\sqrt{2},1]$. Moreover, we prove a multivariate local central limit theorem for structural features of independent sets in $Q_d$ drawn according to the hard core model at any fixed fugacity $\lambda>0$. In proving these results we develop several general tools for performing combinatorial enumeration using polymer models and the cluster expansion from statistical physics along with local central limit theorems.


[132] 2106.09027

On Causal State Updates in Quantum Field Theory

In relativistic Quantum Field Theory (QFT) ideal measurements of certain observables are physically impossible without violating causality. This prompts two questions: i) can a given observable be ideally measured in QFT, and ii) if not, in what sense can it be measured? Here we formulate a necessary and sufficient condition that any measurement, and more generally any state update (quantum operation), must satisfy to respect causality. Our focus is scalar QFT, although our results should be applicable to observables in fermionic QFT. We argue that for unitary `kicks' and operations involving 1-parameter families of Kraus operators, e.g. Gaussian measurements, the only causal observables are smeared fields and the identity - the basic observables in QFT. We provide examples with more complicated operators such as products of smeared fields, and show that the associated state updates are acausal, and hence impossible. Despite this, one can still recover expectation values of such operators, and we show how to do this using only causal measurements of smeared fields.


[133] 2106.09037

Conformal Rigidity from Focusing

The null curvature condition (NCC) is the requirement that the Ricci curvature of a Lorentzian manifold be nonnegative along null directions, which ensures the focusing of null geodesic congruences. In this note, we show that the NCC together with the causal structure significantly constrain the metric. In particular, we prove that any conformal rescaling of a vacuum spacetime introduces either geodesic incompleteness or negative null curvature, provided the conformal factor is non-constant on at least one complete null geodesic. In the context of bulk reconstruction in AdS/CFT, our results combined with the technique of light-cone cuts can be used in vacuum spacetimes to reconstruct the full metric in regions probed by complete null geodesics reaching the boundary. For non-vacuum spacetimes, our results constrain the conformal factor, giving an approximate reconstruction of the metric.


[134] 2106.09049

Hamiltonian analysis of fermions coupled to the Holst action

We report three manifestly Lorentz-invariant Hamiltonian formulations of minimally and nonminimally coupled fermion fields to the Holst action. These formulations are achieved by making a suitable parametrization of both the tetrad and the Lorentz connection, which allows us to integrate out some auxiliary fields without spoiling the local Lorentz symmetry. They have the peculiarity that their noncanonical symplectic structures as well as the phase-space variables for the gravitational sector are real. Moreover, two of these Hamiltonian formulations involve half-densitized fermion fields. We also impose the time gauge on these formulations, which leads to real connections for the gravitational configuration variables. Finally, we perform a symplectomorphism in one of the manifestly Lorentz-invariant Hamiltonian formulations and analyze the resulting formulation, which becomes the Hamiltonian formulation of fermion fields minimally coupled to the Palatini action for particular values of the coupling parameters.


[135] 2106.09071

Pre-processing with Orthogonal Decompositions for High-dimensional Explanatory Variables

Strong correlations between explanatory variables are problematic for high-dimensional regularized regression methods. Due to the violation of the Irrepresentable Condition, the popular LASSO method may suffer from false inclusions of inactive variables. In this paper, we propose pre-processing with orthogonal decompositions (PROD) for the explanatory variables in high-dimensional regressions. The PROD procedure is constructed based upon a generic orthogonal decomposition of the design matrix. We demonstrate by two concrete cases that the PROD approach can be effectively constructed for improving the performance of high-dimensional penalized regression. Our theoretical analysis reveals their properties and benefits for high-dimensional penalized linear regression with LASSO. Extensive numerical studies with simulations and data analysis show the promising performance of the PROD.


[136] 2106.09115

Clustering inference in multiple groups

Inference in clustering is paramount to uncovering inherent group structure in data. Clustering methods which assess statistical significance have recently drawn attention owing to their importance for the identification of patterns in high dimensional data with applications in many scientific fields. We present here a U-statistics based approach, specially tailored for high-dimensional data, that clusters the data into three groups while assessing the significance of such partitions. Because our approach stands on the U-statistics based clustering framework of the methods in R package uclust, it inherits its characteristics being a non-parametric method relying on very few assumptions about the data, and thus can be applied to a wide range of dataset. Furthermore our method aims to be a more powerful tool to find the best partitions of the data into three groups when that particular structure is present. In order to do so, we first propose an extension of the test U-statistic and develop its asymptotic theory. Additionally we propose a ternary non-nested significance clustering method. Our approach is tested through multiple simulations and found to have more statistical power than competing alternatives in all scenarios considered. Applications to peripheral blood mononuclear cells and to image recognition shows the versatility of our proposal, presenting a superior performance when compared with other approaches.


[137] 2106.09121

Scaling-up Diverse Orthogonal Convolutional Networks with a Paraunitary Framework

Enforcing orthogonality in neural networks is an antidote for gradient vanishing/exploding problems, sensitivity by adversarial perturbation, and bounding generalization errors. However, many previous approaches are heuristic, and the orthogonality of convolutional layers is not systematically studied: some of these designs are not exactly orthogonal, while others only consider standard convolutional layers and propose specific classes of their realizations. To address this problem, we propose a theoretical framework for orthogonal convolutional layers, which establishes the equivalence between various orthogonal convolutional layers in the spatial domain and the paraunitary systems in the spectral domain. Since there exists a complete spectral factorization of paraunitary systems, any orthogonal convolution layer can be parameterized as convolutions of spatial filters. Our framework endows high expressive power to various convolutional layers while maintaining their exact orthogonality. Furthermore, our layers are memory and computationally efficient for deep networks compared to previous designs. Our versatile framework, for the first time, enables the study of architecture designs for deep orthogonal networks, such as choices of skip connection, initialization, stride, and dilation. Consequently, we scale up orthogonal networks to deep architectures, including ResNet, WideResNet, and ShuffleNet, substantially increasing the performance over the traditional shallow orthogonal networks.


[138] 2106.09130

Tetrahedron maps, Yang-Baxter maps, and partial linearisations

We study tetrahedron maps, which are set-theoretical solutions to Zamolodchikov's functional tetrahedron equation, and their relations with Yang-Baxter maps, which are set-theoretical solutions to the quantum Yang-Baxter equation. In particular, we clarify the structure of the nonlinear algebraic relations which define linear (parametric) tetrahedron maps (with nonlinear dependence on parameters), and we present several transformations which allow one to obtain new such maps from known ones. Also, we prove that the differential of a (nonlinear) tetrahedron map on a manifold is a tetrahedron map as well. Using the obtained general results, we construct new examples of (parametric) Yang-Baxter and tetrahedron maps. Considered examples include maps associated with integrable systems and matrix groups. In particular, we obtain a parametric family of new linear tetrahedron maps, which are linear approximations for the nonlinear tetrahedron map constructed by Dimakis and M\"uller-Hoissen [arXiv:1708.05694] in a study of soliton solutions of matrix Kadomtsev-Petviashvili (KP) equations. Also, we present invariants for this nonlinear tetrahedron map.


[139] 2106.09148

Efficient Numerical Optimal Control for Ground-State Reset of Open Quantum Systems

This paper presents a numerical optimization framework for solving the ground state reset problem where a number of qubits are dispersively coupled to a readout cavity. We model the open system quantum dynamics using Lindblad's master equation, driven by external control pulses. We develop a basis of density matrices (a parameterization) where each basis element is a density matrix itself. Together with a new objective function, we show how a superposition of the basis elements can be used in such a way that the objective function can be evaluated by solving the master equation for one initial condition only - independent of the system dimension. This enables efficient ground state reset allowing for an increasing number of qubits. Numerical results demonstrate efficient optimal control for ground-state reset of one and two qubits coupled to a readout cavity.


[140] 2106.09167

Generalization of waving-plate theory to multiple interacting swimmers

Early research in aerodynamics and biological propulsion was dramatically advanced by the analytical solutions of Theodorsen, von K\'{a}rm\'{a}n, Wu and others. While these classical solutions apply only to isolated swimmers, the flow interactions between multiple swimmers are relevant to many practical applications, including the schooling and flocking of animal collectives. In this work, we derive a class of solutions that describe the hydrodynamic interactions between an arbitrary number of swimmers in a two-dimensional inviscid fluid. Our approach is rooted in multiply-connected complex analysis and exploits several recent results. Specifically, the transcendental (Schottky-Klein) prime function serves as the basic building block to construct the appropriate conformal maps and leading-edge-suction functions, which allows us to solve the modified Schwarz problem that arises. As such, our solutions generalize classical thin aerofoil theory, specifically Wu's waving-plate analysis, to the case of multiple swimmers. For the case of a pair of interacting swimmers, we develop an efficient numerical implementation that allows rapid computations of the forces on each swimmer. We investigate flow-mediated equilibria and find excellent agreement between our new solutions and previously reported experimental results. Our solutions recover and unify disparate results in the literature, thereby opening the door for future studies into the interactions between multiple swimmers.


[141] 2106.09194

RG Flows and Dynamical Systems

In the context of Wilsonian Renormalization, renormalization group (RG) flows are a set of differential equations that defines how the coupling constants of a theory depend on an energy scale. These equations closely resemble thermodynamical equations and how thermodynamical systems are related to temperature. In this sense, it is natural to look for structures in the flows that show a thermodynamics-like behaviour. The mathematical theory to study these equations is called Dynamical Systems, and applications of that have been used to study RG flows. For example, the classical Zamolodchikov's C-Theorem and its higher-dimensional counterparts, that show that there is a monotonically decreasing function along the flow and it is a property that resembles the second-law of thermodynamics, is related to the Lyapunov function in the context of Dynamical Systems. It can be used to rule out exotic asymptotic behaviours like periodic flows (also known as limit cycles). We also study bifurcation theory and index theories, which have been proposed to be useful in the study of RG flows, the former can be used to explain couplings crossing through marginality and the latter to extract global information about the space the flows lives in. In this dissertation, we also look for applications in holographic RG flows and we try to see if the structural behaviours in holographic theories are the same as the ones in the dual field theory side.


[142] 2106.09207

On the Power of Preconditioning in Sparse Linear Regression

Sparse linear regression is a fundamental problem in high-dimensional statistics, but strikingly little is known about how to efficiently solve it without restrictive conditions on the design matrix. We consider the (correlated) random design setting, where the covariates are independently drawn from a multivariate Gaussian $N(0,\Sigma)$ with $\Sigma : n \times n$, and seek estimators $\hat{w}$ minimizing $(\hat{w}-w^*)^T\Sigma(\hat{w}-w^*)$, where $w^*$ is the $k$-sparse ground truth. Information theoretically, one can achieve strong error bounds with $O(k \log n)$ samples for arbitrary $\Sigma$ and $w^*$; however, no efficient algorithms are known to match these guarantees even with $o(n)$ samples, without further assumptions on $\Sigma$ or $w^*$. As far as hardness, computational lower bounds are only known with worst-case design matrices. Random-design instances are known which are hard for the Lasso, but these instances can generally be solved by Lasso after a simple change-of-basis (i.e. preconditioning). In this work, we give upper and lower bounds clarifying the power of preconditioning in sparse linear regression. First, we show that the preconditioned Lasso can solve a large class of sparse linear regression problems nearly optimally: it succeeds whenever the dependency structure of the covariates, in the sense of the Markov property, has low treewidth -- even if $\Sigma$ is highly ill-conditioned. Second, we construct (for the first time) random-design instances which are provably hard for an optimally preconditioned Lasso. In fact, we complete our treewidth classification by proving that for any treewidth-$t$ graph, there exists a Gaussian Markov Random Field on this graph such that the preconditioned Lasso, with any choice of preconditioner, requires $\Omega(t^{1/20})$ samples to recover $O(\log n)$-sparse signals when covariates are drawn from this model.


[143] 2106.09211

Square Root Principal Component Pursuit: Tuning-Free Noisy Robust Matrix Recovery

We propose a new framework -- Square Root Principal Component Pursuit -- for low-rank matrix recovery from observations corrupted with noise and outliers. Inspired by the square root Lasso, this new formulation does not require prior knowledge of the noise level. We show that a single, universal choice of the regularization parameter suffices to achieve reconstruction error proportional to the (a priori unknown) noise level. In comparison, previous formulations such as stable PCP rely on noise-dependent parameters to achieve similar performance, and are therefore challenging to deploy in applications where the noise level is unknown. We validate the effectiveness of our new method through experiments on simulated and real datasets. Our simulations corroborate the claim that a universal choice of the regularization parameter yields near optimal performance across a range of noise levels, indicating that the proposed method outperforms the (somewhat loose) bound proved here.


[144] 2106.09235

An introduction to the relativistic kinetic theory on curved spacetimes

This article provides a self-contained pedagogical introduction to the relativistic kinetic theory of a dilute gas propagating on a curved spacetime manifold (M,g) of arbitrary dimension. Special emphasis is made on geometric aspects of the theory in order to achieve a formulation which is manifestly covariant on the relativistic phase space. Whereas most previous work has focussed on the tangent bundle formulation, here we work on the cotangent bundle associated with (M,g) which is more naturally adapted to the Hamiltonian framework of the theory. In the first part of this work we discuss the relevant geometric structures of the cotangent bundle T*M, starting with the natural symplectic form on T*M, the one-particle Hamiltonian and the Liouville vector field, defined as the corresponding Hamiltonian vector field. Next, we discuss the Sasaki metric on T*M and its most important properties, including the role it plays for the physical interpretation of the one-particle distribution function. In the second part of this work we describe the general relativistic theory of a collisionless gas, starting with the derivation of the collisionless Boltzmann equation for a neutral simple gas. Subsequently, the description is generalized to a charged gas consisting of several species of particles and the general relativistic Vlasov-Maxwell equations are derived for this system. The last part of this work is devoted to a transparent derivation of the collision term, leading to the general relativistic Boltzmann equation on (M,g). The meaning of global and local equilibrium and the strident restrictions for the existence of the former on a curved spacetime are discussed. We close this article with an application of our formalism to the expansion of a homogeneous and isotropic universe filled with a collisional simple gas and its behavior in the early and late epochs. [abbreviated version]


[145] 2106.09252

Temporal Logic Planning for Minimum-Time Positioning of Multiple Threat-Seduction Decoys

Reusable decoys offer a cost-effective alternative to the single-use hardware commonly applied to protect surface assets from threats. Such decoys portray fake assets to lure threats away from the true asset. To deceive a threat, a decoy first has to position itself such that it can break the radar lock. Considering multiple simultaneous threats, this paper introduces an approach for controlling multiple decoys to minimise the time required to break the locks of all the threats. The method includes the optimal allocation of one decoy to every threat with an assignment procedure that provides local position constraints to guarantee collision avoidance and thereby decouples the control of the decoys. A crude model of a decoy with uncertainty is considered for motion planning. The task of a decoy reaching a state in which the lock of the assigned threat can be broken is formulated as a temporal logic specification. To this end, the requirements to complete the task are modelled as time-varying set-membership constraints. The temporal and logical combination of the constraints is encoded in a mixed-integer optimisation problem. To demonstrate the results a simulated case study is provided.


[146] 2106.09267

Trading with the Crowd

We formulate and solve a multi-player stochastic differential game between financial agents who seek to cost-efficiently liquidate their position in a risky asset in the presence of jointly aggregated transient price impact, along with taking into account a common general price predicting signal. The unique Nash-equilibrium strategies reveal how each agent's liquidation policy adjusts the predictive trading signal to the aggregated transient price impact induced by all other agents. This unfolds a quantitative relation between trading signals and the order flow in crowded markets. We also formulate and solve the corresponding mean field game in the limit of infinitely many agents. We prove that the equilibrium trading speed and the value function of an agent in the finite $N$-player game converges to the corresponding trading speed and value function in the mean field game at rate $O(N^{-2})$. In addition, we prove that the mean field optimal strategy provides an approximate Nash-equilibrium for the finite-player game.


[147] 2106.09276

Uniform Convergence of Interpolators: Gaussian Width, Norm Bounds, and Benign Overfitting

We consider interpolation learning in high-dimensional linear regression with Gaussian data, and prove a generic uniform convergence guarantee on the generalization error of interpolators in an arbitrary hypothesis class in terms of the class's Gaussian width. Applying the generic bound to Euclidean norm balls recovers the consistency result of Bartlett et al. (2020) for minimum-norm interpolators, and confirms a prediction of Zhou et al. (2020) for near-minimal-norm interpolators in the special case of Gaussian data. We demonstrate the generality of the bound by applying it to the simplex, obtaining a novel consistency result for minimum l1-norm interpolators (basis pursuit). Our results show how norm-based generalization bounds can explain and be used to analyze benign overfitting, at least in some settings.


[148] 2106.09315

Fast evaluation of some p-adic transcendental functions

We design algorithms for computing values of many p-adic elementary and special functions, including logarithms, exponentials, polylogarithms, and hypergeometric functions. All our algorithms feature a quasi-linear complexity with respect to the target precision and most of them are based on an adaptation to the-adic setting of the binary splitting and bit-burst strategies.


[149] 2106.09324

Solving the Bose-Hubbard model in new ways

We introduce a new method for analysing the Bose-Hubbard model for an array of bosons with nearest neighbor interactions. It is based on a number-theoretic implementation of the creation and annihilation operators that constitute the model. One of the advantages of this approach is that it facilitates computation with arbitrary accuracy, enabling nearly perfect numerical experimentation. In particular, we provide a rigorous computer assisted proof of quantum phase transitions in finite systems of this type. Furthermore, we investigate properties of the infinite array via harmonic analysis on the multiplicative group of positive rationals. This furnishes an isomorphism that recasts the underlying Fock space as an infinite tensor product of Hecke spaces, i.e., spaces of square-integrable periodic functions that are a superposition of non-negative frequency harmonics. Under this isomorphism, the number-theoretic creation and annihilation operators are mapped into the Kastrup model of the harmonic oscillator on the circle. It also enables us to highlight a kinship of the model at hand with an array of spin moments with a local anisotropy field. This identifies an interesting physical system that can be mapped into the model at hand.


[150] 2106.09327

Minimax Estimation of Partially-Observed Vector AutoRegressions

To understand the behavior of large dynamical systems like transportation networks, one must often rely on measurements transmitted by a set of sensors, for instance individual vehicles. Such measurements are likely to be incomplete and imprecise, which makes it hard to recover the underlying signal of interest.Hoping to quantify this phenomenon, we study the properties of a partially-observed state-space model. In our setting, the latent state $X$ follows a high-dimensional Vector AutoRegressive process $X_t = \theta X_{t-1} + \varepsilon_t$. Meanwhile, the observations $Y$ are given by a noise-corrupted random sample from the state $Y_t = \Pi_t X_t + \eta_t$. Several random sampling mechanisms are studied, allowing us to investigate the effect of spatial and temporal correlations in the distribution of the sampling matrices $\Pi_t$.We first prove a lower bound on the minimax estimation error for the transition matrix $\theta$. We then describe a sparse estimator based on the Dantzig selector and upper bound its non-asymptotic error, showing that it achieves the optimal convergence rate for most of our sampling mechanisms. Numerical experiments on simulated time series validate our theoretical findings, while an application to open railway data highlights the relevance of this model for public transport traffic analysis.


[151] 2106.09353

Supergroup Structure of Jackiw-Teitelboim Supergravity

We develop the gauge theory formulation of $\mathcal{N}=1$ Jackiw-Teitelboim supergravity in terms of the underlying $\text{OSp}(1|2, \mathbb{R})$ supergroup, focusing on boundary dynamics and the exact structure of gravitational amplitudes. We prove that the BF description reduces to a super-Schwarzian quantum mechanics on the holographic boundary, where boundary-anchored Wilson lines map to bilocal operators in the super-Schwarzian theory. A classification of defects in terms of monodromies of $\text{OSp}(1|2, \mathbb{R})$ is carried out and interpreted in terms of character insertions in the bulk. From a mathematical perspective, we construct the principal series representations of $\text{OSp}(1|2, \mathbb{R})$ and show that whereas the corresponding Plancherel measure does not match the density of states of $\mathcal{N}=1$ JT supergravity, a restriction to the positive subsemigroup $\text{OSp}^+(1|2, \mathbb{R})$ yields the correct density of states, mirroring the analogous results for bosonic JT gravity. We illustrate these results with several gravitational applications, in particular computing the late-time complexity growth in JT supergravity.


[152] 2106.09367

Towards sampling complex actions

Path integrals with complex actions are encountered for many physical systems ranging from spin- or mass-imbalanced atomic gases and graphene to quantum chromo-dynamics at finite density to the non-equilibrium evolution of quantum systems. Many computational approaches have been developed for tackling the sign problem emerging for complex actions. Among these, complex Langevin dynamics has the appeal of general applicability. One of its key challenges is the potential convergence of the dynamics to unphysical fixed points. The statistical sampling process at such a fixed point is not based on the physical action and hence leads to wrong predictions. Moreover, its unphysical nature is hard to detect due to the implicit nature of the process. In the present work we set up a general approach based on a Markov chain Monte Carlo scheme in an extended state space. In this approach we derive an explicit real sampling process for generalized complex Langevin dynamics. Subject to a set of constraints, this sampling process is the physical one. These constraints originate from the detailed-balance equations satisfied by the Monte Carlo scheme. This allows us to re-derive complex Langevin dynamics from a new perspective and establishes a framework for the explicit construction of new sampling schemes for complex actions.


[153] 2106.09372

Amplituhedron-like geometries

We consider amplituhedron-like geometries which are defined in a similar way to the intrinsic definition of the amplituhedron but with non-maximal winding number. We propose that for the cases with minimal number of points the canonical form of these geometries corresponds to the product of parity conjugate amplitudes at tree as well as loop level. The product of amplitudes in superspace lifts to a star product in bosonised superspace which we give a precise definition of. We give an alternative definition of amplituhedron-like geometries, analogous to the original amplituhedron definition, and also a characterisation as a sum over pairs of on-shell diagrams that we use to prove the conjecture at tree level. The union of all amplituhedron-like geometries has a very simple definition given by only physical inequalities. Although such a union does not give a positive geometry, a natural extension of the standard definition of canonical form, the globally oriented canonical form, acts on this union and gives the square of the amplitude.


[154] 2106.09375

Recovery under Side Constraints

This paper addresses sparse signal reconstruction under various types of structural side constraints with applications in multi-antenna systems. Side constraints may result from prior information on the measurement system and the sparse signal structure. They may involve the structure of the sensing matrix, the structure of the non-zero support values, the temporal structure of the sparse representationvector, and the nonlinear measurement structure. First, we demonstrate how a priori information in form of structural side constraints influence recovery guarantees (null space properties) using L1-minimization. Furthermore, for constant modulus signals, signals with row-, block- and rank-sparsity, as well as non-circular signals, we illustrate how structural prior information can be used to devise efficient algorithms with improved recovery performance and reduced computational complexity. Finally, we address the measurement system design for linear and nonlinear measurements of sparse signals. Moreover, we discuss the linear mixing matrix design based on coherence minimization. Then we extend our focus to nonlinear measurement systems where we design parallel optimization algorithms to efficiently compute stationary points in the sparse phase retrieval problem with and without dictionary learning.


[155] 2106.09377

A New Dissipativity Condition for Asymptotic Stability of Discounted Economic MPC

Economic Model Predictive Control has recently gained popularity due to its ability to directly optimize a given performance criterion, while enforcing constraint satisfaction for nonlinear systems. Recent research has developed both numerical algorithms and stability analysis for the undiscounted case. The introduction of a discount factor in the cost, however, can be desirable in some cases of interest, e.g., economics, stochastically terminating processes, Markov decision processes, etc. Unfortunately, the stability theory in this case is still not fully developed. In this paper we propose a new dissipativity condition to prove asymptotic stability in the infinite horizon case and we connect our results with existing ones in the literature on discounted economic optimal control. Numerical examples are provided to illustrate the theoretical results.


[156] 2106.09429

Convergence of time-averaged observables towards their steady values for Markov trajectories with application to Skew-Detailed-Balance Lifted-Markov processes

Among the Markov chains breaking detailed-balance that have been proposed in the field of Monte-Carlo sampling in order to accelerate the convergence towards the steady state with respect to the detailed-balance dynamics, the idea of 'Lifting' consists in duplicating the configuration space into two copies $\sigma=\pm$ and in imposing directed flows in each copy in order to explore the configuration space more efficiently. The skew-detailed-balance Lifted-Markov-chain introduced by K. S. Turitsyn, M. Chertkov and M. Vucelja [Physica D Nonlinear Phenomena 240 , 410 (2011)] is revisited for the Curie-Weiss mean-field ferromagnetic model, where the dynamics for the magnetization is closed. The typical convergence properties and the large deviations at various levels for empirical time-averaged observables are analyzed and compared with their detailed-balance counterparts, both for the discrete extensive magnetization $M$ and for the continuous intensive magnetization $m=\frac{M}{N}$ for large system-size $N$.


[157] 2106.09518

Multi-Layered Blockchain Governance Game

This paper deals with design of an integrated secure Blockchain network framework to prevent damages from attackers. The multi-layer concept which could handle multiple number of networks is adapted on the top of Blockchain Governance Game frameworks. This new integrated theoretical model is designed to find the best strategies toward preparation for preventing whole network systems malfunction from attackers and it is developed based on the combination of the Blockchain Governance Game and the Strategic Alliance for Blockchain Governance Game. Analytically tractable results for executing a safety mode are fully obtained and simulated results are demonstrated to obtain the optimal values of hyper parameters of a Blockchain based security network. This research helps those whom are constructing a multiple layer network by enhancing security features through multi-layer framework in a decentralized network.


[158] 2106.09586

Prevalence and Propagation of Fake News

In recent years, scholars have raised concerns on the effects that unreliable news, or "fake news," has on our political sphere, and our democracy as a whole. For example, the propagation of fake news on social media is widely believed to have influenced the outcome of national elections, including the 2016 U.S. Presidential Election, and the 2020 COVID-19 pandemic. What drives the propagation of fake news on an individual level, and which interventions could effectively reduce the propagation rate? Our model disentangles bias from truthfulness of an article and examines the relationship between these two parameters and a reader's own beliefs. Using the model, we create policy recommendations for both social media platforms and individual social media users to reduce the spread of untruthful or highly biased news. We recommend that platforms sponsor unbiased truthful news, focus fact-checking efforts on mild to moderately biased news, recommend friend suggestions across the political spectrum, and provide users with reports about the political alignment of their feed. We recommend that individual social media users fact check news that strongly aligns with their political bias and read articles of opposing political bias.


[159] 2106.09597

Hierarchical surrogate-based Approximate Bayesian Computation for an electric motor test bench

Inferring parameter distributions of complex industrial systems from noisy time series data requires methods to deal with the uncertainty of the underlying data and the used simulation model. Bayesian inference is well suited for these uncertain inverse problems. Standard methods used to identify uncertain parameters are Markov Chain Monte Carlo (MCMC) methods with explicit evaluation of a likelihood function. However, if the likelihood is very complex, such that its evaluation is computationally expensive, or even unknown in its explicit form, Approximate Bayesian Computation (ABC) methods provide a promising alternative. In this work both methods are first applied to artificially generated data and second on a real world problem, by using data of an electric motor test bench. We show that both methods are able to infer the distribution of varying parameters with a Bayesian hierarchical approach. But the proposed ABC method is computationally much more efficient in order to achieve results with similar accuracy. We suggest to use summary statistics in order to reduce the dimension of the data which significantly increases the efficiency of the algorithm. Further the simulation model is replaced by a Polynomial Chaos Expansion (PCE) surrogate to speed up model evaluations. We proof consistency for the proposed surrogate-based ABC method with summary statistics under mild conditions on the (approximated) forward model.


[160] 2106.09601

A gravitational action with stringy $Q$ and $R$ fluxes via deformed differential graded Poisson algebras

We study a deformation of a $2$-graded Poisson algebra where the functions of the phase space variables are complemented by linear functions of parity odd velocities. The deformation is carried by a $2$-form $B$-field and a bivector $\Pi$, that we consider as gauge fields of the geometric and non-geometric fluxes $H$, $f$, $Q$ and $R$ arising in the context of string theory compactification. The technique used to deform the Poisson brackets is widely known for the point particle interacting with a $U(1)$ gauge field, but not in the case of non-abelian or higher spin fields. The construction is closely related to Generalized Geometry: With an element of the algebra that squares to zero, the graded symplectic picture is equivalent to an exact Courant algebroid over the generalized tangent bundle $E \cong TM \oplus T^{*}M$, and to its higher gauge theory. A particular idempotent graded canonical transformation is equivalent to the generalized metric. Focusing on the generalized differential geometry side we construct an action functional with the Ricci tensor of a connection on covectors, encoding the dynamics of a gravitational theory for a contravariant metric tensor and $Q$ and $R$ fluxes. We also extract a connection on vector fields and determine a non-symmetric metric gravity theory involving a metric and $H$-flux.


[161] 2106.09602

Integrable systems of the intermediate long wave type in 2+1 dimensions

We classify 2+1 dimensional integrable systems with nonlocality of the intermediate long wave type. Links to the 2+1 dimensional waterbag system are established. Dimensional reductions of integrable systems constructed in this paper provide dispersive regularisations of hydrodynamic equations governing propagation of long nonlinear waves in a shear flow with piecewise linear velocity profile (for special values of vorticities).


[162] 2106.09646

Impurity-induced increase in the thermal quantum correlations and teleportation in an Ising-XXZ diamond chain

In this work we analyze the quantum correlations in a spin-1/2 Ising-XXZ diamond chain with one plaquette distorted impurity. We have shown that the introduction of impurity into the chain can significantly increase entanglement as well as quantum correlations compared to the original model, without impurity. Due to the great flexibility in the choice of impurity parameters, the model presented is very general and this fact can be very useful for future experimental measurements. In addition to entanglement and quantum coherence, we studied quantum teleportation through a quantum channel composed by a coupled of Heisenberg dimers with distorted impurity in an Ising-XXZ diamond chain, as well as fidelity in teleportation. Our analysis shows that the appropriate choice of parameters can greatly increase all the measures analyzed. For comparison purposes, we present all our results together with the results of the measurements made for the original model, without impurity, studied in previous works.


[163] 2106.09658

Non-intrusive Nonlinear Model Reduction via Machine Learning Approximations to Low-dimensional Operators

Although projection-based reduced-order models (ROMs) for parameterized nonlinear dynamical systems have demonstrated exciting results across a range of applications, their broad adoption has been limited by their intrusivity: implementing such a reduced-order model typically requires significant modifications to the underlying simulation code. To address this, we propose a method that enables traditionally intrusive reduced-order models to be accurately approximated in a non-intrusive manner. Specifically, the approach approximates the low-dimensional operators associated with projection-based reduced-order models (ROMs) using modern machine-learning regression techniques. The only requirement of the simulation code is the ability to export the velocity given the state and parameters as this functionality is used to train the approximated low-dimensional operators. In addition to enabling nonintrusivity, we demonstrate that the approach also leads to very low computational complexity, achieving up to $1000\times$ reduction in run time. We demonstrate the effectiveness of the proposed technique on two types of PDEs.


[164] 2106.09659

Robustness and Consistency in Linear Quadratic Control with Predictions

We study the problem of learning-augmented predictive linear quadratic control. Our goal is to design a controller that balances consistency, which measures the competitive ratio when predictions are accurate, and robustness, which bounds the competitive ratio when predictions are inaccurate. We propose a novel $\lambda$-confident controller and prove that it maintains a competitive ratio upper bound of $1+\min\{O(\lambda^2\varepsilon)+ O(1-\lambda)^2,O(1)+O(\lambda^2)\}$ where $\lambda\in [0,1]$ is a trust parameter set based on the confidence in the predictions, and $\varepsilon$ is the prediction error. Further, we design a self-tuning policy that adaptively learns the trust parameter $\lambda$ with a regret that depends on $\varepsilon$ and the variation of perturbations and predictions.


[165] 2106.09683

PAC-Bayes, MAC-Bayes and Conditional Mutual Information: Fast rate bounds that handle general VC classes

We give a novel, unified derivation of conditional PAC-Bayesian and mutual information (MI) generalization bounds. We derive conditional MI bounds as an instance, with special choice of prior, of conditional MAC-Bayesian (Mean Approximately Correct) bounds, itself derived from conditional PAC-Bayesian bounds, where `conditional' means that one can use priors conditioned on a joint training and ghost sample. This allows us to get nontrivial PAC-Bayes and MI-style bounds for general VC classes, something recently shown to be impossible with standard PAC-Bayesian/MI bounds. Second, it allows us to get faster rates of order $O \left(({\text{KL}}/n)^{\gamma}\right)$ for $\gamma > 1/2$ if a Bernstein condition holds and for exp-concave losses (with $\gamma=1$), which is impossible with both standard PAC-Bayes generalization and MI bounds. Our work extends the recent work by Steinke and Zakynthinou [2020] who handle MI with VC but neither PAC-Bayes nor fast rates, the recent work of Hellstr\"om and Durisi [2020] who extend the latter to the PAC-Bayes setting via a unifying exponential inequality, and Mhammedi et al. [2019] who initiated fast rate PAC-Bayes generalization error bounds but handle neither MI nor general VC classes.


[166] 2106.09689

Statistical Query Lower Bounds for List-Decodable Linear Regression

We study the problem of list-decodable linear regression, where an adversary can corrupt a majority of the examples. Specifically, we are given a set $T$ of labeled examples $(x, y) \in \mathbb{R}^d \times \mathbb{R}$ and a parameter $0< \alpha <1/2$ such that an $\alpha$-fraction of the points in $T$ are i.i.d. samples from a linear regression model with Gaussian covariates, and the remaining $(1-\alpha)$-fraction of the points are drawn from an arbitrary noise distribution. The goal is to output a small list of hypothesis vectors such that at least one of them is close to the target regression vector. Our main result is a Statistical Query (SQ) lower bound of $d^{\mathrm{poly}(1/\alpha)}$ for this problem. Our SQ lower bound qualitatively matches the performance of previously developed algorithms, providing evidence that current upper bounds for this task are nearly best possible.