Search

In 1903, Epstein published his proof of meromorphic continuation and a functional equation for Dirichlet series associated with quadratic forms, now called Epstein zeta-functions. However, already in 1889 (or even earlier) Hurwitz was aware of these results as his mathematical diaries and some unpublished notes (in an almost final form) found in his estate at the ETH Zurich show. In this article we present and analyze Hurwitz's notes and compare his reasoning with Epstein's paper in detail.

The spectral gap conjecture of M. van den Berg [2, formula (65)] asserts that λ2 − λ1 ≥ 3π for all convex euclidean domains of diameter 1, where λ1 and λ2 denote the first two eigenvalues of the Dirichlet Laplacian. Notice that equality holds for the 1-dimensional unit interval, which can be regarded also as a degenerate n-dimensional rectangular box. The gap estimate is conjectured to hold more generally for Schrodinger operators with convex potentials, under Dirichlet boundary conditions; see the work of S.-T. Yau and collaborators [9, 11]. This Schrodinger gap conjecture was proved some time ago in 1 dimension by R. Lavine [8], and more recently in all dimensions by B. Andrews and J. Clutterbuck [1]. The proof in this journal by H. Donnelly [3] of the original gap conjecture in 2 dimensions (for the Dirichlet Laplacian with zero potential) is not correct. The Editors of Mathematische Zeitschrift have asked us to describe the flaws in the proof, in order to clarify the state of the literature. Donnelly’s approach to the problem is a natural one: first perform a shape optimization to rule out a non-degenerate minimizing domain, and then analyze the spectral gap for a sequence of domains degenerating to an interval, with the help of results by D. Jerison [5]. (For some history on this approach, and on the gap conjecture more generally, see the report on the AIM meeting “Low Eigenvalues of Laplace and Schrodinger Operators” [10], especially page 12 of the open problems list.) The error lies in the proof of the shape optimization step, as we now explain. Donnelly wishes to prove that no minimizing domain can exist for

In this paper, we analyze the complexity of an eco-epidemiological model for phytoplankton–zooplankton system in presence of toxicity and time delay. Holling type II function response is incorporated to address the predation rate as well as toxic substance distribution in zooplankton. It is also presumed that infected phytoplankton does recover from the viral infection. In the absence of time delay, stability and Hopf-bifurcation conditions are investigated to explore the system dynamics around all the possible equilibrium points. Further, in the presence of time delay, conditions for local stability are derived around the interior equilibria and the properties of the periodic solution are obtained by applying normal form theory and central manifold arguments. Computational simulation is performed to illustrate our theoretical findings. It is explored that system dynamics is very sensitive corresponding to carrying capacity and toxin liberation rate and able to generate chaos. Further, it is observed that time delay in the viral infection process can destabilize the phytoplankton density whereas zooplankton density remains in its old state. Incorporation of time delay also gives the scenario of double Hopf-bifurcation. Some control parameters are discussed to stabilize system dynamics. The effect of time delay on (i) growth rate of susceptible phytoplankton shows the extinction and double Hopf-bifurcation in the zooplankton population, (ii) a sufficiently large value of carrying capacity stabilizes the chaotic dynamics or makes the whole system chaotic with further increment.

Månsson and Shukur (Econ Model 28:1475–1481, 2011) proposed a Poisson ridge regression estimator (PRRE) to reduce the negative effects of multicollinearity. However, a weakness of the PRRE is its relatively large bias. Therefore, as a remedy, Türkan and Özel (J Appl Stat 43:1892–1905, 2016) examined the performance of almost unbiased ridge estimators for the Poisson regression model. These estimators will not only reduce the consequences of multicollinearity but also decrease the bias of PRRE and thus perform more efficiently. The aim of this paper is twofold. Firstly, to derive the mean square error properties of the Modified Almost Unbiased PRRE (MAUPRRE) and Almost Unbiased PRRE (AUPRRE) and then propose new ridge estimators for MAUPRRE and AUPRRE. Secondly, to compare the performance of the MAUPRRE with the AUPRRE, PRRE and maximum likelihood estimator. Using both simulation study and real-world dataset from the Swedish football league, it is evidenced that one of the proposed, MAUPRRE ($$ \hat{k}_{q4} $$ k ^ q 4 ) performed better than the rest in the presence of high to strong (0.80–0.99) multicollinearity situation.

© 2019, Pleiades Publishing, Ltd. It is well known that the classical and Sobolev wave fronts were extended to nonequivalent global versions by the use of the short-time Fourier transform. In this very short paper, we give complete characterizations of the former wave front sets in terms of the short-time Fourier transform.

In this short paper we combine the representability theorem introduced in [17, 18] with the theory of derived formal models introduced in [2] to prove the existence representability of the derived Hilbert space RHilb(X) for a separated k-analytic space X. Such representability results relies on a localization theorem stating that if X is a quasi-compact and quasi-separated formal scheme, then the \infty-category Coh^+(X^rig) of almost perfect complexes over the generic fiber can be realized as a Verdier quotient of the \infty-category Coh^+(X). Along the way, we prove several results concerning the the \infty-categories of formal models for almost perfect modules on derived k-analytic spaces., 28 pages

The "exact subgraph" approach was recently introduced as a hierarchical scheme to get increasingly tight semidefinite programming relaxations of several NP-hard graph optimization problems. Solving these relaxations is a computational challenge because of the potentially large number of violated subgraph constraints. We introduce a computational framework for these relaxations designed to cope with these difficulties. We suggest a partial Lagrangian dual, and exploit the fact that its evaluation decomposes into several independent subproblems. This opens the way to use the bundle method from non-smooth optimization to minimize the dual function. Finally computational experiments on the Max-Cut, stable set and coloring problem show the excellent quality of the bounds obtained with this approach., arXiv admin note: substantial text overlap with arXiv:1902.05345

We investigate the asymptotic properties of the trajectories generated by a second-order dynamical system with Hessian driven damping and a Tikhonov regularization term in connection with the minimization of a smooth convex function in Hilbert spaces. We obtain fast convergence results for the function values along the trajectories. The Tikhonov regularization term enables the derivation of strong convergence results of the trajectory to the minimizer of the objective function of minimum norm.

We propose a complete-search algorithm for solving a class of non-convex, possibly infinite-dimensional, optimization problems to global optimality. We assume that the optimization variables are in a bounded subset of a Hilbert space, and we determine worst-case run-time bounds for the algorithm under certain regularity conditions of the cost functional and the constraint set. Because these run-time bounds are independent of the number of optimization variables and, in particular, are valid for optimization problems with infinitely many optimization variables, we prove that the algorithm converges to an \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varepsilon $$\end{document}ε-suboptimal global solution within finite run-time for any given termination tolerance \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varepsilon > 0$$\end{document}ε>0. Finally, we illustrate these results for a problem of calculus of variations.

In this short paper we study the d-Hermite condition about stably free modules for skew $$\textit{PBW}$$ extensions. For this purpose, we estimate the stable rank of these non-commutative rings. In addition, and closely related with these questions, we will prove Kronecker’s theorem about the radical of finitely generated ideals for some particular types of skew $$\textit{PBW}$$ extensions.

Although planar quadratic differential systems and their applications have been studied in more than one thousand papers, we still have no complete understanding of these systems. In this paper we have two objectives. First we provide a brief bibliographical survey on the main results about quadratic systems. Here we do not consider the applications of these systems to many areas as in Physics, Chemist, Economics, Biology, … Second we characterize the new class of planar separable quadratic polynomial differential systems. For such class of systems we provide the normal forms which contain one parameter, and using the Poincare compactification and the blow up technique, we prove that there exist 10 non-equivalent topological phase portraits in the Poincare disc for the separable quadratic polynomial differential systems.

The singular value decomposition of a complex matrix is a fundamental concept in linear algebra and has proved extremely useful in many subjects. It is less clear what the situation is over a finite field. In this paper, we classify the orbits of GU m ( q ) × GU n ( q ) on M m × n ( q 2 ) (which is the analog of the singular value decomposition). The proof involves Kronecker’s theory of pencils and the Lang–Steinberg theorem for algebraic groups. Besides the motivation mentioned above, this problem came up in a recent paper of Guralnick et al. (2020) where a concept of character level for the complex irreducible characters of finite, general or special, linear and unitary groups was studied and bounds on the number of orbits was needed. A consequence of this work determines possible pairs of Jordan forms for nilpotent matrices of the form A A ∗ and A ∗ A over a finite field and A A ⊤ and A ⊤ A over arbitrary fields.

In this paper, we generalize the classical Freidlin-Wentzell’s theorem for random perturbations of Hamiltonian systems. In (Probability Theory and Related Fields 128 (2004) 441–466), M.Freidlin and M.Weber generalized the original result in the sense that the coefficient for the noise term is no longer the identity matrix but a state-dependent matrix and taking the drift term into consideration. In this paper, We generalize the result by adding a state-dependent matrix that converges uniformly to 0 on any compact sets as ϵ tends to 0 to a state-dependent noise and considering the drift term which contains two parts, the state-dependent mapping and a state-dependent mapping that converges uniformly to 0 on any compact sets as ϵ tends to 0. In the proof, we adapt a new way to prove the weak convergence inside the edge by constructing an auxiliary process and modify the proof in (Probability Theory and Related Fields 128 (2004) 441–466) when proving gluing condition.

In many applications it is desirable to consider not one random vector but a number of random vectors with the joint distribution. This paper is devoted to the integral and integral transformations connected with the joint vector Gaussian probability density function. Such integral and transformations arise in the statistical decision theory, particularly, in the dual control theory based on the statistical decision theory. One of the results represented in the paper is the integral of the joint Gaussian probability density function. The other results are the total probability formula and Bayes formula formulated in terms of the joint vector Gaussian probability density function. As an example the Bayesian estimations of the coefficients of the multiple regression function are obtained. The proposed integrals can be used as table integrals in various fields of research.

This paper considers a continuous-time quasi birth-death (qbd) process, which informally can be seen as a birth-death process of which the parameters are modulated by an external continuous-time Markov chain. The aim is to numerically approximate the time-dependent distribution of the resulting bivariate Markov process in an accurate and efficient way. An approach based on the Erlangization principle is proposed and formally justified. Its performance is investigated and compared with two existing approaches: one based on numerical evaluation of the matrix exponential underlying the qbd process, and one based on the uniformization technique. It is shown that in many settings the approach based on Erlangization is faster than the other approaches, while still being highly accurate. In the last part of the paper, we demonstrate the use of the developed technique in the context of the evaluation of the likelihood pertaining to a time series, which can then be optimized over its parameters to obtain the maximum likelihood estimator. More specifically, through a series of examples with simulated and real-life data, we show how it can be deployed in model selection problems that involve the choice between a qbd and its non-modulated counterpart.

This paper presents some new formulas and properties of the generalized fractional derivative with non-singular kernel that covers various types of fractional derivatives such as the Caputo–Fabrizio fractional derivative, the Atangana–Baleanu fractional derivative, and the weighted Atangana–Baleanu fractional derivative. These new properties extend many recent results existing in the literature. Furthermore, the paper proposes some interesting inequalities that estimate the generalized fractional derivatives of some specific functions. These inequalities can be used to construct Lyapunov functions with the aim to study the global asymptotic stability of several fractional-order systems arising from diverse fields of science and engineering.

Statistics are mathematical tools applying scientific investigations, such as engineering and medical and biological analyses. However, statistical methods are often improved. Nowadays, statisticians try to find an accurate way to solve a problem. One of these problems is estimation parameters, which can be expressed as an inverse problem when independent variables are highly correlated. This paper’s significant goal is to interpret the parameter estimates of double generalized Rayleigh distribution in a regression model using a wavelet basis. It is difficult to use the standard version of the regression methods in practical terms, which is obtained using the likelihood. Since a noise level usually makes the result of estimation unstable, multicollinearity leads to various estimates. This kind of problem estimates that features of the truth are complicated. So it is reasonable to use a mixed method that combines a fully Bayesian approach and a wavelet basis. The usual rule for wavelet approaches is to choose a wavelet basis, where it helps to compute the wavelet coefficients, and then, these coefficients are used to remove Gaussian noise. Recovering data is typically calculated by inverting the wavelet coefficients. Some wavelet bases have been considered, which provide a shift-invariant wavelet transform, simultaneously providing improvements in smoothness, in recovering, and in squared-error performance. The proposed method uses combining a penalized maximum likelihood approach, a penalty term, and wavelet tools. In this paper, real data are involved and modeled using double generalized Rayleigh distributions, as they are used to estimate the wavelet coefficients of the sample using numerical tools. In practical applications, wavelet approaches are recommended. They reduce noise levels. This process may be useful since the noise level is often corrupted in real data, as a significant cause of most numerical estimation problems. A simulation investigation is studied using the MCMC tool to estimate the underlying features as an essential task statistics.

Let $S \subset \Bbb R^n$ be a smooth compact hypersurface with a strictly positive second fundamental form, $E$ be the Fourier extension operator on $S$, and $X$ be a Lebesgue measurable subset of $\Bbb R^n$. If $X$ contains a ball of each radius, then the problem of determining the range of exponents $(p,q)$ for which the estimate $\| Ef \|_{L^q(X)} \leq C \| f \|_{L^p(S)}$ holds is equivalent to the restriction conjecture. In this paper, we study the estimate under the following assumption on the set $X$: there is a number $0 < \alpha \leq n$ such that $|X \cap B_R| \leq c \, R^\alpha$ for all balls $B_R$ in $\Bbb R^n$ of radius $R \geq 1$. On the left-hand side of this estimate, we are integrating the function $|Ef(x)|^q$ against the measure $\chi_X dx$. Our approach consists of replacing the characteristic function $\chi_X$ of $X$ by an appropriate weight function $H$, and studying the resulting estimate in three different regimes: small values of $\alpha$, intermediate values of $\alpha$, and large values of $\alpha$. In the first regime, we establish the estimate by using already available methods. In the second regime, we prove a weighted H\"{o}lder-type inequality that holds for general non-negative Lebesgue measurable functions on $\Bbb R^n$, and combine it with the result from the first regime. In the third regime, we borrow a recent fractal Fourier restriction theorem of Du and Zhang and combine it with the result from the second regime. In the opposite direction, the results of this paper improve on the Du-Zhang theorem in the range $0 < \alpha < n/2$., Comment: 31 pages. Minor revision

As more long-span bridges continue to be completed and opened to traffic, the safety of cars driving across the bridge has attracted more and more attention, especially when the car is suddenly affected by the crosswind, the car is likely to have direction deviation or even a rollover accident. In this paper, the large eddy simulation method is used to study the flow field characteristics and safety of the car on the bridge under the turbulent crosswind. The numerical simulation model is established by referring to the Donghai Bridge, and the correctness of the car model is validated by combining with the data of wind tunnel test. The influence of factors such as the porosity and height of the bridge guardrail and the Reynolds number of airflow on the flow field characteristics is analyzed. The study shows that, in order to ensure the safety of cars on the bridge, the bridge guardrail porosity should be small, 35.8% is more suitable, the guardrail height should be more suitable within the range of 1.5–1.625 meters, and the Reynolds number should not be 3.51e + 5. The research results of this paper will provide reference for the optimal design of bridge guardrail.

A classical theorem of Koksma states that for Lebesgue almost every $x>1$ the sequence $(x^n)_{n=1}^{\infty}$ is uniformly distributed modulo one. In the present paper we extend Koksma's theorem to the pair correlation setting. More precisely, we show that for Lebesgue almost every $x>1$ the pair correlations of the fractional parts of $(x^n)_{n=1}^{\infty}$ are asymptotically Poissonian. The proof is based on a martingale approximation method., Version 2: some minor changes. The paper will appear in the Israel Journal of Mathematics

PurposeThe purpose of this paper is to study the Bertotti–Kasner space-time and its geometric properties.Design/methodology/approachThis paper is based on the features of λ-tensor and the technique of six-dimensional formalism introduced by Pirani and followed by W. Borgiel, Z. Ahsan et al. and H.M. Manjunatha et al. This technique helps to describe both the geometric properties and the nature of the gravitational field of the space-times in the Segre characteristic.FindingsThe Gaussian curvature quantities specify the curvature of Bertotti–Kasner space-time. They are expressed in terms of invariants of the curvature tensor. The Petrov canonical form and the Weyl invariants have also been obtained.Originality/valueThe findings are revealed to be both physically and geometrically interesting for the description of the gravitational field of the cylindrical universe of Bertotti–Kasner type as far as the literature is concerned. Given the technique of six-dimensional formalism, the authors have defined the Weyl conformal λW-tensor and discussed the canonical form of the Weyl tensor and the Petrov scalars. To the best of the literature survey, this idea is found to be modern. The results deliver new insight into the geometry of the nonstatic cylindrical vacuum solution of Einstein's field equations.

In this paper we establish general weighted Hardy identities for several subelliptic settings including Hardy identities on the Heisenberg group, Carnot groups with respect to a homogeneous gauge and Carnot–Carathéodory metric, general nilpotent groups, and certain families of Hörmander vector fields. We also introduce new weighted uncertainty principles in these settings. This is done by continuing the program initiated by [N. Lam, G. Lu and L. Zhang, Factorizations and Hardy’s-type identities and inequalities on upper half spaces, Calc. Var. Partial Differential Equations 58 2019, 6, Paper No. 183; N. Lam, G. Lu and L. Zhang, Geometric Hardy’s inequalities with general distance functions, J. Funct. Anal. 279 2020, 8, Article ID 108673] of using the Bessel pairs introduced by [N. Ghoussoub and A. Moradifam, Functional Inequalities: New Perspectives and New Applications, Math. Surveys Monogr. 187, American Mathematical Society, Providence, 2013] to obtain Hardy identities. Using these identities, we are able to improve significantly existing Hardy inequalities in the literature in the aforementioned subelliptic settings. In particular, we establish the Hardy identities and inequalities in the spirit of [H. Brezis and J. L. Vázquez, Blow-up solutions of some nonlinear elliptic problems, Rev. Mat. Univ. Complut. Madrid 10 1997, 443–469] and [H. Brezis and M. Marcus, Hardy’s inequalities revisited. Dedicated to Ennio De Giorgi, Ann. Sc. Norm. Super. Pisa Cl. Sci. (4) 25 1997, 1–2, 217–237] in these settings.

The study in this paper is made on the nonlinear fractional differential equation whose nonlinearity involves the explicit fractional order D0+β u(t). The corresponding Green's function is derived first, and then the completely continuous operator is proved. Besides, based on the Schauder's fixed point theorem and the Krasnosel'skii's fixed point theorem, the sufficient conditions for at least one or two existence of positive solutions are established. Furthermore, several other sufficient conditions for at least three, n or 2n-1 positive solutions are also obtained by applying the generalized AveryHenderson fixed point theorem and the Avery-Peterson fixed point theorem. Finally, several simulation examples are provided to illustrate the main results of the paper. In particularly, a novel efficient iterative method is employed for simulating the examples mentioned above, that is, the interesting point of this paper is that the approximation graphics for the solutions are given by using the iterative method.

This paper concerns the existence of transonic shock solutions to the 2-D steady compressible Euler system in an almost flat finite nozzle ( in the sense that it is a generic small perturbation of a flat one ), under physical boundary conditions proposed by Courant-Friedrichs in \cite{CourantFriedrichs1948}, in which the receiver pressure is prescribed at the exit of the nozzle. In the resulting free boundary problem, the location of the shock-front is one of the most desirable information one would like to determine. However, the location of the normal shock-front in a flat nozzle can be anywhere in the nozzle so that it provides little information on the possible location of the shock-front when the nozzle's boundary is perturbed. So one of the key difficulties in looking for transonic shock solutions is to determine the shock-front. To this end, a free boundary problem for the linearized Euler system will be proposed, whose solution will be taken as an initial approximation for the transonic shock solution. In this paper, a sufficient condition in terms of the geometry of the nozzle and the given exit pressure is derived which yields the existence of the solutions to the proposed free boundary problem. Once an initial approximation is obtained, a further nonlinear iteration could be constructed and proved to lead to a transonic shock solution., 53 pages

We establish a theory of $Q$-valued functions minimizing a suitable generalization of the Dirichlet integral. In a second paper the theory will be used to approximate efficiently area minimizing currents $\mathrm{mod}(p)$ when $p=2Q$, and to establish a first general partial regularity theorem for every $p$ in any dimension and codimension., 37 pages. First part of a two-papers work aimed at establishing a first general partial regularity theory for area minimizing currents modulo p, for any p and in any dimension and codimension. v3 is the final version, to appear on Comm. Pure Appl. Math

This paper provides a rigorous derivation of a counterexample of Bourgain, related to a well-known question of pointwise a.e. convergence for the solution of the linear Schrödinger equation, for initial data in a Sobolev space. This counterexample combines ideas from analysis and number theory, and the present paper demonstrates how to build such counterexamples from first principles, and then optimize them.

In this paper, we consider a class of fractional-order systems described by the Caputo derivative. The behaviors of the dynamics of this particular class of fractional-order systems will be proposed and experienced by a numerical scheme to obtain the phase portraits. Before that, we will provide the conditions under which the considered fractional-order system’s solution exists and is unique. The fractional-order impact will be analyzed, and the advantages of the fractional-order derivatives in modeling chaotic systems will be discussed. How the parameters of the model influence the considered fractional-order system will be studied using the Lyapunov exponents. The topological changes of the systems and the detection of the chaotic and hyperchaotic behaviors at the assumed initial conditions and the considered fractional-order systems will also be investigated using the Lyapunov exponents. The investigations related to the Lyapunov exponents in the context of the fractional-order derivative will be the main novelty of this paper. The stability analysis of the model’s equilibrium points has been focused in terms of the Matignon criterion.

In this paper, an alternative approach is provided in the well-posedness analysis of elliptic variational–hemivariational inequalities in real Hilbert spaces. This includes the unique solvability and continuous dependence of the solution on the data. In most of the existing literature on elliptic variational–hemivariational inequalities, well-posedness results are obtained by using arguments of surjectivity for pseudomonotone multivalued operators, combined with additional compactness and pseudomonotonicity properties. In contrast, following (Han in Nonlinear Anal B Real World Appl 54:103114, 2020; Han in Numer Funct Anal Optim 42:371–395, 2021), the approach adopted in this paper is based on the fixed point structure of the problems, combined with minimization principles for elliptic variational–hemivariational inequalities. Consequently, only elementary results of functional analysis are needed in the approach, which makes the theory of elliptic variational–hemivariational inequalities more accessible to applied mathematicians and engineers. The theoretical results are illustrated on a representative example from contact mechanics.

In this paper, we provide the boundedness property of the Riesz transforms associated to the Schrodinger operator $${\cal L} = \Delta + {\bf{V}}$$ in a new weighted Morrey space which is the generalized version of many previous Morrey type spaces. The additional potential V considered in this paper is a non-negative function satisfying the suitable reverse Holder’s inequality. Our results are new and general in many cases of problems. As an application of the boundedness property of these singular integral operators, we obtain some regularity results of solutions to Schrodinger equations in the new Morrey space.

In this paper, we first use the distribution of the number of records to demonstrate that the right tail probabilities of counts of rare events are generally better approximated by the right tail probabilities of Poisson distribution than {those} of normal distribution. We then show the moderate deviations in Poisson approximation generally require an adjustment and, with suitable adjustment, we establish better error estimates of the moderate deviations in Poisson approximation than those in \cite{CFS}. Our estimates contain no unspecified constants and are easy to apply. We illustrate the use of the theorems in six applications: Poisson-binomial distribution, matching problem, occupancy problem, birthday problem, random graphs and 2-runs. The paper complements the works of \cite{CC92,BCC95,CFS}., 29 pages and 5 figures

It is well known that when the geometry and/or coefficients allow stable trapped rays, the solution operator of the Helmholtz equation (a.k.a. the resolvent of the Laplacian) grows exponentially through a sequence of real frequencies tending to infinity. In this paper we show that, even in the presence of the strongest-possible trapping, if a set of frequencies of arbitrarily small measure is excluded, the Helmholtz solution operator grows at most polynomially as the frequency tends to infinity. One significant application of this result is in the convergence analysis of several numerical methods for solving the Helmholtz equation at high frequency that are based on a polynomial-growth assumption on the solution operator (e.g. $hp$-finite elements, $hp$-boundary elements, certain multiscale methods). The result of this paper shows that this assumption holds, even in the presence of the strongest-possible trapping, for most frequencies.

In a series of papers published since 2017 the author introduced a simple alternative of the $n$-simplex type, to enhance the accuracy of approximations of second-order boundary value problems subject to Dirichlet boundary conditions, posed on smooth curved domains. This technique is based upon trial functions consisting of piecewise polynomials defined on straight-edged triangular or tetrahedral meshes, interpolating the Dirichlet boundary conditions at points of the true boundary. In contrast, the test functions are defined by the standard degrees of freedom associated with the underlying method for polytopic domains. While the mathematical analysis of the method for Lagrange and Hermite methods for two-dimensional second- and fourth-order problems was carried out in earlier paper by the author this paper is devoted to the study of the three-dimensional case. Well-posedness, uniform stability and optimal a priori error estimates in the energy norm are proved for a tetrahedron-based Lagrange family of finite elements. Novel error estimates in the $L^2$-norm, for the class of problems considered in this work, are also proved. A series of numerical examples illustrates the potential of the new technique. In particular, its superior accuracy at equivalent cost, as compared to the isoparametric technique, is highlighted.

An adaptive high-order difference solution about a 2D nonlinear degenerate singular reaction-diffusion equation with a convection term is initially proposed in the paper. After the first and the second central difference operator approximating the first-order and the second-order spatial derivative, respectively, the higher-order spatial derivatives are discretized by applying the Taylor series rule and the temporal derivative is discretized by using the Crank–Nicolson (CN) difference scheme. An alternating direction implicit (ADI) scheme with a nonuniform grid is built in this way. Meanwhile, accuracy analysis declares the second order in time and the fourth order in space under certain conditions. Sequentially, the high-order scheme is performed on an adaptive mesh to demonstrate quenching behaviors of the singular parabolic equation and analyse the influence of combustion chamber size on quenching. The paper displays rationally that the proposed scheme is practicable for solving the 2D quenching-type problem.

Recently, the authors have established $$L^p$$ -boundedness of vector-valued q-variational inequalities for averaging operators which take values in the Banach space satisfying the martingale cotype q property in Hong and Ma (Math Z 286(1–2):89–120, 2017). In this paper, we prove that the martingale cotype q property is also necessary for the vector-valued q-variational inequalities, which was a question left open in the previous paper. Moreover, we also prove that the UMD property and the martingale cotype q property can be characterized in terms of vector valued q-variational inequalities for the Hilbert transform.

In this paper, we show how to use the framework of mod-Gaussian convergence in order to study the fluctuations of certain models of random graphs, of random permutations and of random integer partitions. We prove that, in these three frameworks, a generic homogeneous observable of a generic random model is mod-Gaussian under an appropriate renormalisation. This implies a central limit theorem with an extended zone of normality, a moderate deviation principle, an estimate of the speed of convergence, a local limit theorem and a concentration inequality. The universal asymptotic behavior of the observables of these models gives rise to a notion of mod-Gaussian moduli space., Comment: New version: the paper has been slightly shortened, and a few references were added. 52 pages, 13 figures

In this paper we analyze a class of trace finite element methods for the discretization of vector-Laplace equations. A key issue in the finite element discretization of such problems is the treatment of the constraint that the unknown vector field must be tangential to the surface (‘tangent condition’). We study three different natural techniques for treating the tangent condition, namely a consistent penalty method, a simpler inconsistent penalty method and a Lagrange multiplier method. The main goal of the paper is to present an analysis that reveals important properties of these three different techniques for treating the tangent constraint. A detailed error analysis is presented that takes the approximation of both the geometry of the surface and the solution of the partial differential equation into account. Error bounds in the energy norm are derived that show how the discretization error depends on relevant parameters such as the degree of the polynomials used for the approximation of the solution, the degree of the polynomials used for the approximation of the level set function that characterizes the surface, the penalty parameter and the degree of the polynomials used for the approximation of the Lagrange multiplier.

International audience; In this paper, we consider a semi-linear fourth order parabolic equation in a bounded smooth domain Ω with homogeneous Dirichlet and Neumann boundary conditions. The main result of this paper is the null controllability and the exact controllability to the trajectories at any time T > 0 for the associated control system with a control function acting at the interior.; Dans ce papier, on considère uneéquation parabolique semi-linéaire de quatrième ordre dans un domaine borné régulier Ω avec des conditions aux limites de type Dirichlet et Neumann homogènes. Le résultat principal de ce papier concerne la contrôlabilitéà zéro et la contrôlabilité exacte pour tout T > 0 du système de contrôle associé avec un contrôle agissantà l'interieur.

This paper considers some multiscale radial basis function collocation methods for solving the two-dimensional Monge–Ampère equation with Dirichlet boundary. We discuss and study the performance of the three kinds of multiscale methods. The first method is the cascadic meshfree method, which was proposed by Liu and He (2013). The second method is the stationary multilevel method, which was proposed by Floater and Iske (1996), and is used to solve the fully nonlinear partial differential equation in the paper for the first time. The third is the hierarchical radial basis function method, which is constructed by employing successive refinement scattered data sets and scaled compactly supported radial basis functions with varying support radii. Compared with the first two methods, the hierarchical radial basis function method can not only solve the present problem on a single level with higher accuracy and lower computational cost but also produce highly sparse nonlinear discrete system. These observations are obtained by taking the direct approach of numerical experimentation.

In this paper, we establish the following Stein–Weiss inequality with the fractional Poisson kernel. Then we prove that there exist extremals for the Stein–Weiss inequality (⋆), and that the extremals must be radially decreasing about the origin. We also provide the regularity and asymptotic estimates of positive solutions to the integral systems which are the Euler–Lagrange equations of the extremals to the Stein–Weiss inequality (⋆) with the fractional Poisson kernel. Our result is inspired by the work of Hang, Wang and Yan, where the Hardy–Littlewood–Sobolev type inequality was first established when γ=2 and α=β=0. The proof of the Stein–Weiss inequality (⋆) with the fractional Poisson kernel in this paper uses recent work on the Hardy–Littlewood–Sobolev inequality with the fractional Poisson kernel by Chen, Lu and Tao, and the present paper is a further study in this direction.

Let (G, V) be a prehomogeneous vector space, let O be any G(F_q)-invariant subset of V(F_q), and let f be the characteristic function of O. In this paper we develop a method for explicitly and efficiently evaluating the Fourier transform of f, based on combinatorics and linear algebra. We then carry out these computations in full for each of five prehomogeneous vector spaces, including the 12-dimensional space of pairs of ternary quadratic forms. Our computations reveal that these Fourier transforms enjoy a great deal of structure, and sometimes exhibit more than square root cancellation on average. These Fourier transforms naturally arise in analytic number theory, where explicit formulas (or upper bounds) lead to sieve level of distribution results for related arithmetic sequences. We describe some examples, and in a companion paper we develop a new method to do so, designed to exploit the particular structure of these Fourier transforms., 28 pages; to appear in American Journal of Mathematics. Includes a table after the bibliography. Also includes relevant PARI/GP code, embedded as a comment in the LaTeX code

This paper is devoted to the study of dual singular integral equations with convolution kernels in the case of non-normal type. Via using the Fourier transforms, we transform such equations into Riemann boundary value problems. To solve the equation, we establish the regularity theory of solvability. The general solutions and the solvable conditions of the equation are obtained. Especially, we investigate the asymptotic property of solutions at nodes. This paper will have a significant meaning for the study of improving and developing complex analysis, integral equations and Riemann boundary value problems.

This paper contains an explanation of Ramanujan-type formulas with cubic radicals of cubic irrationalities in the situation when these radicals are contained in a pure cubic extension. We give a complete description of formulas of such type, answering the Zippel’s question. It turns out that Ramanujan-type formulas are in some sense unique in this situation. In particular, there must be no more than three summands in the right-hand side and the norm of the irrationality in question must be a cube. In this situation we associate cubic irrationalities with a cyclic cubic polynomial, which is reducible if and only if one can simplify the corresponding cubic radical. This correspondence is inverse to the so-called Ramanujan correspondence defined in the preceding papers, where one associates a pure cubic extension to some cyclic polynomial.

This paper aims to present a new pathwise approximation method, which gives approximate solutions of order $\begin{array}{} \displaystyle \frac{3}{2} \end{array}$ for stochastic differential equations (SDEs) driven by multidimensional Brownian motions. The new method, which assumes the diffusion matrix non-degeneracy, employs the Runge-Kutta method and uses the Itô-Taylor expansion, but the generating of the approximation of the expansion is carried out as a whole rather than individual terms. The new idea we applied in this paper is to replace the iterated stochastic integrals Iα by random variables, so implementing this scheme does not require the computation of the iterated stochastic integrals Iα. Then, using a coupling which can be found by a technique from optimal transport theory would give a good approximation in a mean square. The results of implementing this new scheme by MATLAB confirms the validity of the method.

[EN] In this paper, a novel technique is formed to obtain the solution of a fractional gas dynamics equation. Some reproducing kernel Hilbert spaces are defined. Reproducing kernel functions of these spaces have been found. Some numerical examples are shown to confirm the efficiency of the reproducing kernel Hilbert space method. The accurate pulchritude of the paper is arisen in its strong implementation of Caputo fractional order time derivative on the classical equations with the success of the highly accurate solutions by the series solutions. Reproducing kernel Hilbert space method is actually capable of reducing the size of the numerical work. Numerical results for different particular cases of the equations are given in the numerical section., This research was partially supported by Spanish Ministerio de Ciencia, Innovacion y Universidades PGC2018-095896-B-C22 and Generalitat Valenciana PROMETEO/2016/089.

The purpose of this paper is to show well-posedness results in L 2 -based Sobolev spaces for transmission, Dirichlet, Neumann, and mixed boundary value problems for the Brinkman system with L ∞ coefficients in Lipschitz domains on a compact Riemannian manifold of dimension m ≥ 2 . The Dirichlet, transmission, and mixed problems for the nonlinear Darcy-Forchheimer-Brinkman system with L ∞ coefficients are also analyzed. First, we focus on the well-posedness of linear transmission, Dirichlet and mixed boundary value problems for the Brinkman system with L ∞ coefficients in Lipschitz domains on compact Riemannian manifolds by using a variational approach that reduces such a boundary value problem to a mixed variational formulation defined in terms of two bilinear continuous forms, one of them satisfying a coercivity condition and another one the inf-sup condition. Further, we show the equivalence between each boundary value problem for the Brinkman system with L ∞ coefficients and its mixed variational counterpart, and then the well posedness in L 2 -based Sobolev spaces by using the Necas-Babuska-Brezzi technique. The second goal of this paper is the construction of the Newtonian and layer potential operators for the Brinkman system with L ∞ coefficients in Lipschitz domains on compact Riemannian manifolds by using the well-posedness results for the analyzed linear transmission problems. Various mapping properties of these operators are also obtained and used to describe the weak solutions of the Poisson problems with Dirichlet and Neumann conditions for the nonsmooth Brinkman system in terms of such potentials. Finally, we combine the well-posedness results of the Poisson problems of Dirichlet, transmission, and mixed type for the nonsmooth Brinkman system with a fixed point theorem in order to show the existence of a weak solution of the Poisson problem of Dirichlet, transmission, or mixed type for the (nonlinear) Darcy-Forchheimer-Brinkman system with L ∞ coefficients in L 2 -based Sobolev spaces in Lipschitz domains on compact Riemannian manifolds of dimension m ∈ { 2 , 3 } .

Though the boundary layer formation in the chemotactic process has been observed in experiment (cf. [63] ), the mathematical study on the boundary layer solutions of chemotaxis models is just in its infant stage. Apart from the sophisticated theoretical tools involved in the analysis, how to impose/derive physical boundary conditions is a state-of-the-art in studying the boundary layer problem of chemotaxis models. This paper will proceed with a previous work [24] in one dimension to establish the convergence of boundary layer solutions of the Keller–Segel model with singular sensitivity in a two-dimensional space (half-plane) with respect to the chemical diffusion rate denoted by e ≥ 0 . Compared to the one-dimensional boundary layer problem, there are many new issues arising from multi-dimensions such as possible Prandtl type degeneracy, curl-free preservation and well-posedness of large-data solutions. In this paper, we shall derive appropriate physical boundary conditions and gradually overcome these barriers and hence establish the convergence of boundary layer solutions of the singular Keller–Segel system in the half-plane as the chemical diffusion rate vanishes. Specially speaking, we justify that the boundary layer converges to the outer layer (solution with e = 0 ) plus the inner layer as e → 0 , where both outer and inner layer profiles are precisely derived and well understood. By doing this, the structure of boundary layer solutions is clearly characterized. We hope that our results and methods can shed lights on the understanding of underlying mechanisms of the boundary layer patterns observed in the experiment for chemotaxis such as the work by Tuval et al. [63] , and open a new window in the future theoretical study of chemotaxis models.

In some previous works, two of the authors introduced a technique to design high-order numerical methods for one-dimensional balance laws that preserve all their stationary solutions. The basis of these methods is a well-balanced reconstruction operator. Moreover, they introduced a procedure to modify any standard reconstruction operator, like MUSCL, ENO, CWENO, etc., in order to be well-balanced. This strategy involves a non-linear problem at every cell at every time step that consists in finding the stationary solution whose average is the given cell value. In a recent paper, a fully well-balanced method is presented where the non-linear problems to be solved in the reconstruction procedure are interpreted as control problems. The goal of this paper is to introduce a new technique to solve these local non-linear problems based on the application of the collocation RK methods. Special care is put to analyze the effects of computing the averages and the source terms using quadrature formulas. A general technique which allows us to deal with resonant problems is also introduced. To check the efficiency of the methods and their well-balance property, they have been applied to a number of tests, ranging from easy academic systems of balance laws consisting of Burgers equation with some non-linear source terms to the shallow water equations—without and with Manning friction—or Euler equations of gas dynamics with gravity effects.

This paper explains a method of constructing algebras, starting with the properties of discrimination in elementary discrete systems. We show how to use points of view about these systems to construct what we call iterant algebras and how these algebras naturally give rise to the complex numbers, Clifford algebras and matrix algebras. The paper discusses the structure of the Schrödinger equation, the Dirac equation and the Majorana Dirac equations, finding solutions via the nilpotent method initiated by Peter Rowlands.

This is a continuation of our previous works on entanglement entropy (EE) in interacting field theories. In arXiv:2103.05303, we have proposed the notion of $\mathbb{Z}_M$ gauge theory on Feynman diagrams to calculate EE in quantum field theories and shown that EE consists of two particular contributions from propagators and vertices. As shown in the next paper arXiv:2105.02598, the purely non-Gaussian contributions from interaction vertices can be interpreted as renormalized correlation functions of composite operators. In this paper, we will first provide a unified matrix form of EE containing both contributions from propagators and (classical) vertices, and then extract further non-Gaussian contributions based on the framework of the Wilsonian renormalization group. It is conjectured that the EE in the infrared is given by a sum of all the vertex contributions in the Wilsonian effective action., Comment: 29 pages, 10 figures; typos corrected, published version in Symmetry (v2)

Chinta and Gunnells introduced a rather intricate multi-parameter Weyl group action on rational functions on a torus, which, when the parameters are specialized to certain Gauss sums, describes the functional equations of Weyl group multiple Dirichlet series associated to metaplectic (n-fold) covers of algebraic groups. In subsequent joint work with Puskas, they extended this action to a "metaplectic" representation of the equal parameter affine Hecke algebra, which allowed them to obtain explicit formulas for the p-parts of these Dirichlet series. They have also verified by a computer check the remarkable fact that their formulas continue to define a group action for general (unspecialized) parameters. In the first part of paper we give a conceptual explanation of this fact, by giving a uniform and elementary construction of the "metaplectic" representation for generic Hecke algebras as a suitable quotient of a parabolically induced affine Hecke algebra module, from which the associated Chinta-Gunnells Weyl group action follows through localization. In the second part of the paper we extend the metaplectic representation to the double affine Hecke algebra, which provides a generalization of Cherednik's basic representation. This allows us to introduce a new family of "metaplectic" polynomials, which generalize nonsymmetric Macdonald polynomials. In this paper, we provide the details of the construction of metaplectic polynomials in type A; the general case will be handled in the sequel to this paper., 39 pages. Version 2 is a significant revision. Added second part introducing a new family of "metaplectic" polynomials, which generalize nonsymmetric Macdonald polynomials and metaplectic Iwahori-Whittaker functions. Title has been changed and the introduction has been expanded