Showing 1–50 of 106 results for author: Tian, X

Non-uniform Cocycles for Some Uniquely Ergodic Minimal Dynamical Systems on Connected Spaces

Authors: Wanshan Lin, Xueting Tian

Abstract: In this paper, we pay attention to a weaker version of Walters's question on the existence of non-uniform cocycles for uniquely ergodic minimal dynamical systems on non-degenerate connected spaces. We will classify such dynamical systems into three classes: not totally uniquely ergodic; totally uniquely ergodic but not topological weakly mixing; totally uniquely ergodic and topological weakly mixi… ▽ More In this paper, we pay attention to a weaker version of Walters's question on the existence of non-uniform cocycles for uniquely ergodic minimal dynamical systems on non-degenerate connected spaces. We will classify such dynamical systems into three classes: not totally uniquely ergodic; totally uniquely ergodic but not topological weakly mixing; totally uniquely ergodic and topological weakly mixing. We will give an affirmative answer to such question for the first two classes. Also, we will show the existence of such dynamical systems in the first class with arbitrary topological entropy. △ Less

Submitted 5 September, 2024; originally announced September 2024.

Numerical analysis of a class of penalty discontinuous Galerkin methods for nonlocal diffusion problems

Authors: Qiang Du, Lili Ju, Jianfang Lu, Xiaochuan Tian

Abstract: In this paper, we consider a class of discontinuous Galerkin (DG) methods for one-dimensional nonlocal diffusion (ND) problems. The nonlocal models, which are integral equations, are widely used in describing many physical phenomena with long-range interactions. The ND problem is the nonlocal analog of the classic diffusion problem, and as the interaction radius (horizon) vanishes, then the nonloc… ▽ More In this paper, we consider a class of discontinuous Galerkin (DG) methods for one-dimensional nonlocal diffusion (ND) problems. The nonlocal models, which are integral equations, are widely used in describing many physical phenomena with long-range interactions. The ND problem is the nonlocal analog of the classic diffusion problem, and as the interaction radius (horizon) vanishes, then the nonlocality disappears and the ND problem converges to the classic diffusion problem. Under certain conditions, the exact solution to the ND problem may exhibit discontinuities, setting it apart from the classic diffusion problem. Since the DG method shows its great advantages in resolving problems with discontinuities in computational fluid dynamics over the past several decades, it is natural to adopt the DG method to compute the ND problems. Based on [Du-Ju-Lu-Tian-CAMC2020], we develop the DG methods with different penalty terms, ensuring that the proposed DG methods have local counterparts as the horizon vanishes. This indicates the proposed methods will converge to the existing DG schemes as the horizon vanishes, which is crucial for achieving asymptotic compatibility. Rigorous proofs are provided to demonstrate the stability, error estimates, and asymptotic compatibility of the proposed DG schemes. To observe the effect of the nonlocal diffusion, we also consider the time-dependent convection-diffusion problems with nonlocal diffusion. We conduct several numerical experiments, including accuracy tests and Burgers' equation with nonlocal diffusion, and various horizons are taken to show the good performance of the proposed algorithm and validate the theoretical findings. △ Less

Submitted 13 August, 2024; originally announced August 2024.

MSC Class: 65M60; 65R20; 45A05

Uniqueness of positive radial solutions of Choquard type equations

Authors: Tao Wang, Xiaoyu Tian, Hui Guo

Abstract: In this paper, we consider the following Choquard type equation \begin{equation} \left\{\begin{aligned} &-Δu+λu=γ(Φ_N(|x|)\ast|u|^p)u \ \ \mbox{in $\mathbb{R}^N$}, \\ &\lim\limits_{|x|\to\infty}u(x)=0,\\ \end{aligned}\right. \end{equation} where $N\geq2,λ>0,γ>0, p\in[1,2]$ and $Φ_N(|x|)$ denotes the fundamental solution of the Laplacian $-Δ$ on $\mathbb{R}^N$. This equation does not have a varia… ▽ More In this paper, we consider the following Choquard type equation \begin{equation} \left\{\begin{aligned} &-Δu+λu=γ(Φ_N(|x|)\ast|u|^p)u \ \ \mbox{in $\mathbb{R}^N$}, \\ &\lim\limits_{|x|\to\infty}u(x)=0,\\ \end{aligned}\right. \end{equation} where $N\geq2,λ>0,γ>0, p\in[1,2]$ and $Φ_N(|x|)$ denotes the fundamental solution of the Laplacian $-Δ$ on $\mathbb{R}^N$. This equation does not have a variational frame when $p\neq 2.$ Instead of variational methods, we prove the existence and uniqueness of positive radial solutions of the above equation via the shooting method by establishing some new differential inequalities. The proofs are based on an analysis of the corresponding system of second-order differential equations, and our results extend the existing ones in the literature from $p=2$ to $p\in[1,2]$. △ Less

Submitted 1 August, 2024; originally announced August 2024.

Comments: 18pages

MSC Class: 35A02; 35B07; 35J08; 35J47

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

Authors: Qihao Ye, Xiaochuan Tian, Dong Wang

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

Submitted 21 July, 2024; originally announced July 2024.

Comments: 34 pages, 8 figures

MSC Class: 34K37; 44A35; 35Q84; 65D40; 33C10

Some aspects of semi-harmonious quasi-projection pairs

Authors: Xiaoyi Tian, Qingxiang Xu, Chunhong Fu

Abstract: A term called the quasi-projection pair $(P,Q)$ was introduced recently by the authors, where $P$ is a projection and $Q$ is an idempotent on a Hilbert $C^*$-module $H$ satisfying $Q^*=(2P-I)Q(2P-I)$, in which $Q^*$ is the adjoint operator of the idempotent $Q$ and $I$ is the identity operator on $H$. Some fundamental issues on quasi-projection pairs, such as the block matrix representations for q… ▽ More A term called the quasi-projection pair $(P,Q)$ was introduced recently by the authors, where $P$ is a projection and $Q$ is an idempotent on a Hilbert $C^*$-module $H$ satisfying $Q^*=(2P-I)Q(2P-I)$, in which $Q^*$ is the adjoint operator of the idempotent $Q$ and $I$ is the identity operator on $H$. Some fundamental issues on quasi-projection pairs, such as the block matrix representations for quasi-projection pairs and the $C^*$-morphisms associated with quasi-projection pairs, are worthwhile to be investigated. This paper aims to make some preparations. One object called the semi-harmonious quasi-projection pair is introduced in the general setting of the adjointable operators on Hilbert $C^*$-modules. Some related operator theories on the common similarity of operators and a norm equation associated with the Friedrichs angle are dealt with. △ Less

Submitted 16 August, 2024; v1 submitted 15 April, 2024; originally announced April 2024.

Comments: arXiv admin note: substantial text overlap with arXiv:2305.12984

MSC Class: 46L08; 47A05

Compactness results for a Dirichlet energy of nonlocal gradient with applications

Authors: Zhaolong Han, Tadele Mengesha, Xiaochuan Tian

Abstract: We prove two compactness results for function spaces with finite Dirichlet energy of half-space nonlocal gradients. In each of these results, we provide sufficient conditions on a sequence of kernel functions that guarantee the asymptotic compact embedding of the associated nonlocal function spaces into the class of square-integrable functions. Moreover, we will demonstrate that the sequence of no… ▽ More We prove two compactness results for function spaces with finite Dirichlet energy of half-space nonlocal gradients. In each of these results, we provide sufficient conditions on a sequence of kernel functions that guarantee the asymptotic compact embedding of the associated nonlocal function spaces into the class of square-integrable functions. Moreover, we will demonstrate that the sequence of nonlocal function spaces converges in an appropriate sense to a limiting function space. As an application, we prove uniform Poincaré-type inequalities for sequence of half-space gradient operators. We also apply the compactness result to demonstrate the convergence of appropriately parameterized nonlocal heterogeneous anisotropic diffusion problems. We will construct asymptotically compatible schemes for these type of problems. Another application concerns the convergence and robust discretization of a nonlocal optimal control problem. △ Less

Submitted 21 August, 2024; v1 submitted 13 April, 2024; originally announced April 2024.

The operator distances from projections to an idempotent

Authors: Xiaofeng Zhang, Xiaoyi Tian, Qingxiang Xu

Abstract: The main purpose of this paper is to give a full characterization of the operator distances from projections to an idempotent, which includes the minimum value, the maximum value and the intermediate values. Let $H$ be a Hilbert space and $\mathbb{B}(H)$ be the set of bounded linear operators on $H$. Given an arbitrary idempotent $Q\in \mathbb{B}(H)$, it is proved that… ▽ More The main purpose of this paper is to give a full characterization of the operator distances from projections to an idempotent, which includes the minimum value, the maximum value and the intermediate values. Let $H$ be a Hilbert space and $\mathbb{B}(H)$ be the set of bounded linear operators on $H$. Given an arbitrary idempotent $Q\in \mathbb{B}(H)$, it is proved that $$\|m(Q)-Q\|\le \|P-Q\|\le \|I-m(Q)-Q\|$$ for every projection $P$ on $H$, in which $I$ is the identity operator on $H$ and $m(Q)$ is a specific projection called the matched projection of $Q$. When $Q\in\mathbb{B}(H)$ is a non-projection idempotent, it is proved that for every number $α$ contained in the interval $\left[\|m(Q)-Q\|,\|I-m(Q)-Q\|\right]$, there exists a projection $P\in \mathbb{B}(H)$ such that $\|P-Q\|=α$. Two uniqueness problems concerning the projections that attain the minimum value or the maximum value are also dealt with. △ Less

Submitted 4 April, 2024; originally announced April 2024.

MSC Class: 47A05

Asymptotically compatible schemes for nonlinear variational models via Gamma-convergence and applications to nonlocal problems

Authors: Qiang Du, James M. Scott, Xiaochuan Tian

Abstract: We present a study on asymptotically compatible Galerkin discretizations for a class of parametrized nonlinear variational problems. The abstract analytical framework is based on variational convergence, or Gamma-convergence. We demonstrate the broad applicability of the theoretical framework by developing asymptotically compatible finite element discretizations of some representative nonlinear no… ▽ More We present a study on asymptotically compatible Galerkin discretizations for a class of parametrized nonlinear variational problems. The abstract analytical framework is based on variational convergence, or Gamma-convergence. We demonstrate the broad applicability of the theoretical framework by developing asymptotically compatible finite element discretizations of some representative nonlinear nonlocal variational problems on a bounded domain. These include nonlocal nonlinear problems with classically-defined, local boundary constraints through heterogeneous localization at the boundary, as well as nonlocal problems posed on parameter-dependent domains. △ Less

Submitted 12 February, 2024; originally announced February 2024.

Gradient stability of Caffarelli-Kohn-Nirenberg inequality involving weighted p-Laplace

Authors: Shengbing Deng, Xingliang Tian

Abstract: The best constant and extremal functions are well known of the following Caffarelli-Kohn-Nirenberg inequality \[ \int_{\mathbb{R}^N}|\nabla u|^p\frac{\mathrm{d}x}{|x|^μ}\geq \mathcal{S} \left(\int_{\mathbb{R}^N}|u|^r\frac{\mathrm{d}x}{|x|^s} \right)^{\frac{p}{r}}, \quad \mbox{for all}\quad u\in C^\infty_c(\mathbb{R}^N), \] where $1<p<p+μ<N$, $\fracμ{p}\leq \frac{s}{r}<\fracμ{p}+1$,… ▽ More The best constant and extremal functions are well known of the following Caffarelli-Kohn-Nirenberg inequality \[ \int_{\mathbb{R}^N}|\nabla u|^p\frac{\mathrm{d}x}{|x|^μ}\geq \mathcal{S} \left(\int_{\mathbb{R}^N}|u|^r\frac{\mathrm{d}x}{|x|^s} \right)^{\frac{p}{r}}, \quad \mbox{for all}\quad u\in C^\infty_c(\mathbb{R}^N), \] where $1<p<p+μ<N$, $\fracμ{p}\leq \frac{s}{r}<\fracμ{p}+1$, $r=\frac{p(N-s)}{N-p-μ}$. An important task is investigating the stability of extremals for this inequality. Firstly, we give the classification to the linearized problem related to the extremals which shows the extremals are non-degenerate. Then we investigate the gradient type remainder term of previous inequality by using spectral estimate combined with a compactness argument which partially extends the work of Wei and Wu [Math. Ann., 2022] to a general $p$-Laplace case, and also the work of Figalli and Zhang [Duke Math. J., 2022] to a weighted case. △ Less

Submitted 5 January, 2024; originally announced January 2024.

Comments: 38 pages. Any suggestions and comments are welcome! arXiv admin note: text overlap with arXiv:2308.04111

A note on the $L^p$-Sobolev inequality

Authors: Shengbing Deng, Xingliang Tian

Abstract: The usual Sobolev inequality in $\mathbb{R}^N$, asserts that $\|\nabla u\|_{L^p(\mathbb{R}^N)} \geq \mathcal{S}\|u\|_{L^{p^*}(\mathbb{R}^N)}$ for $1<p<N$ and $p^*=\frac{pN}{N-p}$, with $\mathcal{S}$ being the sharp constant. This note is concerned, instead, with function restricted to bounded domain $Ω\subset \mathbb{R}^N$. Based on the recent work of Figalli and Zhang [Duke Math. J., 2022], a rem… ▽ More The usual Sobolev inequality in $\mathbb{R}^N$, asserts that $\|\nabla u\|_{L^p(\mathbb{R}^N)} \geq \mathcal{S}\|u\|_{L^{p^*}(\mathbb{R}^N)}$ for $1<p<N$ and $p^*=\frac{pN}{N-p}$, with $\mathcal{S}$ being the sharp constant. This note is concerned, instead, with function restricted to bounded domain $Ω\subset \mathbb{R}^N$. Based on the recent work of Figalli and Zhang [Duke Math. J., 2022], a remainder term with weak norm is established \[ \frac{\|\nabla u\|_{L^p(Ω)}}{\|u\|_{L^{p^*}(Ω)}} -\mathcal{S} \geq \mathcal{C} \left(\frac{\|u\|_{L^{\bar{p}}_w(Ω)}} {\|u\|_{L^{p^*}(Ω)}}\right)^{\max\{2,p\}},\quad \forall u\in C^\infty_0(Ω)\setminus\{0\}, \] for some $\mathcal{C}=\mathcal{C}(N,p,Ω)>0$, where $\bar{p}=p^*(p-1)/p$ and $\|\cdot\|_{L^{\bar{p}}_w(Ω)}$ denotes the weak $L^{\bar{p}}$-norm. Furthermore, the weak norm can not be replaced by the strong norm. This result answers the long-standing open problem raised by Bianchi and Egnell [J. Funct. Anal., 1991]. △ Less

Submitted 31 December, 2023; originally announced January 2024.

Gaussian radial basis functions collocation for fractional PDEs: methodology and error analysis

Authors: Xiaochuan Tian, Yixuan Wu, Yanzhi Zhang

Abstract: The paper introduces a new meshfree pseudospectral method based on Gaussian radial basis functions (RBFs) collocation to solve fractional Poisson equations. Hypergeometric functions are used to represent the fractional Laplacian of Gaussian RBFs, enabling an efficient computation of stiffness matrix entries. Unlike existing RBF-based methods, our approach ensures a Toeplitz structure in the stiffn… ▽ More The paper introduces a new meshfree pseudospectral method based on Gaussian radial basis functions (RBFs) collocation to solve fractional Poisson equations. Hypergeometric functions are used to represent the fractional Laplacian of Gaussian RBFs, enabling an efficient computation of stiffness matrix entries. Unlike existing RBF-based methods, our approach ensures a Toeplitz structure in the stiffness matrix with equally spaced RBF centers, enabling efficient matrix-vector multiplications using fast Fourier transforms. We conduct a comprehensive study on the shape parameter selection, addressing challenges related to ill-conditioning and numerical stability. The main contribution of our work includes rigorous stability analysis and error estimates of the Gaussian RBF collocation method, representing a first attempt at the rigorous analysis of RBF-based methods for fractional PDEs to the best of our knowledge. We conduct numerical experiments to validate our analysis and provide practical insights for implementation. △ Less

Submitted 28 December, 2023; originally announced December 2023.

A neural network kernel decomposition for learning multiple steady states in parameterized dynamical systems

Authors: Yimeng Zhang, Alexander Cloninger, Bo Li, Xiaochuan Tian

Abstract: We develop a machine learning approach to identifying parameters with steady-state solutions, locating such solutions, and determining their linear stability for systems of ordinary differential equations and dynamical systems with parameters. Our approach begins with the construction of target functions that can be used to identify parameters with steady-state solution and the linear stability of… ▽ More We develop a machine learning approach to identifying parameters with steady-state solutions, locating such solutions, and determining their linear stability for systems of ordinary differential equations and dynamical systems with parameters. Our approach begins with the construction of target functions that can be used to identify parameters with steady-state solution and the linear stability of such solutions. We design a parameter-solution neural network (PSNN) that couples a parameter neural network and a solution neural network to approximate the target function, and develop efficient algorithms to train the PSNN and to locate steady-state solutions. We also present a theory of approximation of the target function by our PSNN based on the neural network kernel decomposition. Numerical results are reported to show that our approach is robust in identifying the phase boundaries separating different regions in the parameter space corresponding to no solution or different numbers of solutions and in classifying the stability of solutions. These numerical results also validate our analysis. Although the primary focus in this study centers on steady states of parameterized dynamical systems, our approach is applicable generally to finding solutions for parameterized nonlinear systems of algebraic equations. Some potential improvements and future work are discussed. △ Less

Submitted 15 December, 2023; originally announced December 2023.

The generalized polar decomposition, the weak complementarity and the parallel sum for adjointable operators on Hilbert $C^*$-modules

Authors: Xiaofeng Zhang, Xiaoyi Tian, Qingxiang Xu

Abstract: This paper deals mainly with some aspects of the adjointable operators on Hilbert $C^*$-modules. A new tool called the generalized polar decomposition for each adjointable operator is introduced and clarified. As an application, the general theory of the weakly complementable operators is set up in the framework of Hilbert $C^*$-modules. It is proved that there exists an operator equation which ha… ▽ More This paper deals mainly with some aspects of the adjointable operators on Hilbert $C^*$-modules. A new tool called the generalized polar decomposition for each adjointable operator is introduced and clarified. As an application, the general theory of the weakly complementable operators is set up in the framework of Hilbert $C^*$-modules. It is proved that there exists an operator equation which has a unique solution, whereas this unique solution fails to be the reduced solution. Some investigations are also carried out in the Hilbert space case. It is proved that there exist a closed subspace $M$ of certain Hilbert space $K$ and an operator $T\in \mathbb{B}(K)$ such that $T$ is $(M,M)$-weakly complementable, whereas $T$ fails to be $(M,M)$-complementable. The solvability of the equation $$A:B=X^*AX+(I-X)^*B(I-X) \quad (X\in\mathbb{B}(H))$$ is also dealt with in the Hilbert space case, where $A,B\in \mathbb{B}(H)$ are two general positive operators, and $A:B$ denotes their parallel sum. Among other things, it is shown that there exist certain positive operators $A$ and $B$ on the Hilbert space $\ell^2(\mathbb{N})\oplus \ell^2(\mathbb{N})$ such that the above equation has no solution. △ Less

Submitted 24 April, 2024; v1 submitted 12 December, 2023; originally announced December 2023.

Comments: This version is accepted for publication in Banach Journal of Mathematical Analysis

MSC Class: 46L08; 47A05

The Frobenious distances from projections to an idempotent matrix

Authors: Xiaoyi Tian, Qingxiang Xu, Chunhong Fu

Abstract: For each pair of matrices $A$ and $B$ with the same order, let $\|A-B\|_F$ denote their Frobenius distance. This paper deals mainly with the Frobenius distances from projections to an idempotent matrix. For every idempotent $Q\in \mathbb{C}^{n\times n}$, a projection $m(Q)$ called the matched projection can be induced. It is proved that $m(Q)$ is the unique projection whose Frobenius distance away… ▽ More For each pair of matrices $A$ and $B$ with the same order, let $\|A-B\|_F$ denote their Frobenius distance. This paper deals mainly with the Frobenius distances from projections to an idempotent matrix. For every idempotent $Q\in \mathbb{C}^{n\times n}$, a projection $m(Q)$ called the matched projection can be induced. It is proved that $m(Q)$ is the unique projection whose Frobenius distance away from $Q$ takes the minimum value among all the Frobenius distances from projections to $Q$, while $I_n-m(Q)$ is the unique projection whose Frobenius distance away from $Q$ takes the maximum value. Furthermore, it is proved that for every number $α$ between the minimum value and the maximum value, there exists a projection $P$ whose Frobenius distance away from $Q$ takes the value $α$. Based on the above characterization of the minimum distance, some Frobenius norm upper bounds and lower bounds of $\|P-Q\|_F$ are derived under the condition of $PQ=Q$ on a projection $P$ and an idempotent $Q$. △ Less

Submitted 17 December, 2023; v1 submitted 2 December, 2023; originally announced December 2023.

Journal ref: Linear Algebra Appl. 688 (2024), 21--43

Data-driven optimal control with neural network modeling of gradient flows

Authors: Xuping Tian, Baskar Ganapathysubramanian, Hailiang Liu

Abstract: Extracting physical laws from observation data is a central challenge in many diverse areas of science and engineering. We propose Optimal Control Neural Networks (OCN) to learn the laws of vector fields in dynamical systems, with no assumption on their analytical form, given data consisting of sampled trajectories. The OCN framework consists of a neural network representation and an optimal contr… ▽ More Extracting physical laws from observation data is a central challenge in many diverse areas of science and engineering. We propose Optimal Control Neural Networks (OCN) to learn the laws of vector fields in dynamical systems, with no assumption on their analytical form, given data consisting of sampled trajectories. The OCN framework consists of a neural network representation and an optimal control formulation. We provide error bounds for both the solution and the vector field. The bounds are shown to depend on both the training error and the time step between the observation data. We also demonstrate the effectiveness of OCN, as well as its generalization ability, by testing on several canonical systems, including the chaotic Lorenz system. △ Less

Submitted 2 December, 2023; originally announced December 2023.

Comments: 28 pages, 8 figures

MSC Class: 93C15; 49K15

Adaptive Preconditioned Gradient Descent with Energy

Authors: Hailiang Liu, Levon Nurbekyan, Xuping Tian, Yunan Yang

Abstract: We propose an adaptive step size with an energy approach for a suitable class of preconditioned gradient descent methods. We focus on settings where the preconditioning is applied to address the constraints in optimization problems, such as the Hessian-Riemannian and natural gradient descent methods. More specifically, we incorporate these preconditioned gradient descent algorithms in the recently… ▽ More We propose an adaptive step size with an energy approach for a suitable class of preconditioned gradient descent methods. We focus on settings where the preconditioning is applied to address the constraints in optimization problems, such as the Hessian-Riemannian and natural gradient descent methods. More specifically, we incorporate these preconditioned gradient descent algorithms in the recently introduced Adaptive Energy Gradient Descent (AEGD) framework. In particular, we discuss theoretical results on the unconditional energy-stability and convergence rates across three classes of objective functions. Furthermore, our numerical results demonstrate excellent performance of the proposed method on several test bed optimization problems. △ Less

Submitted 14 June, 2024; v1 submitted 10 October, 2023; originally announced October 2023.

Comments: 32 pages, 3 figures

MSC Class: 65K10; 90C26

Symmetry breaking of extremals for the high order Caffarelli-Kohn-Nirenberg type inequalities

Authors: Shengbing Deng, Xingliang Tian

Abstract: Based on the work of C.-S. Lin [Comm. Partial Differential Equations, 1986], an equivalent second-order Caffarelli-Kohn-Nirenberg type inequality will be established, namely \begin{equation*} \int_{\mathbb{R}^N}|x|^{-β}|\mathrm{div} (|x|^α\nabla u)|^2 \mathrm{d}x \geq \mathcal{S}\left(\int_{\mathbb{R}^N} |x|^β|u|^{p^*_{α,β}} \mathrm{d}x\right)^{\frac{2}{p^*_{α,β}}}, \quad \forall u\in C^\i… ▽ More Based on the work of C.-S. Lin [Comm. Partial Differential Equations, 1986], an equivalent second-order Caffarelli-Kohn-Nirenberg type inequality will be established, namely \begin{equation*} \int_{\mathbb{R}^N}|x|^{-β}|\mathrm{div} (|x|^α\nabla u)|^2 \mathrm{d}x \geq \mathcal{S}\left(\int_{\mathbb{R}^N} |x|^β|u|^{p^*_{α,β}} \mathrm{d}x\right)^{\frac{2}{p^*_{α,β}}}, \quad \forall u\in C^\infty_0(\mathbb{R}^N), \end{equation*} for some constant $\mathcal{S}=\mathcal{S}(N,α,β)>0$, where \begin{align*} N\geq 5,\quad α>2-N,\quad α-2<β\leq \frac{N}{N-2}α,\quad p^*_{α,β}=\frac{2(N+β)}{N-4+2α-β}. \end{align*} We obtain a symmetry breaking conclusion: when $α>0$ and $β_{\mathrm{FS}}(α)<β< \frac{N}{N-2}α$ where $β_{\mathrm{FS}}(α):= -N+\sqrt{N^2+α^2+2(N-2)α}$, then the extremal function for the best constant $\mathcal{S}$, if it exists, is nonradial. Furthermore, we give a symmetry result of extremal functions when $β=\frac{N}{N-2}α$ and $2-N<α<0$, i.e., \begin{equation*} \int_{\mathbb{R}^N}|x|^{-\frac{N}{N-2}α}|\mathrm{div} (|x|^α\nabla u)|^2 \mathrm{d}x \geq \left(1+\fracα{N-2}\right)^{4-\frac{4}{N}}\mathcal{S}_0 \left(\int_{\mathbb{R}^N} |x|^{\frac{N}{N-2}α}|u|^{\frac{2N}{N-4}} \mathrm{d}x\right)^{\frac{N-4}{N}}, \quad \forall u\in C^\infty_0(\mathbb{R}^N), \end{equation*} where $\mathcal{S}_0$ is the sharp constant of second-order Sobolev inequality, and equality holds only for some radial functions of the form $(1+|x|^{2+\frac{2α}{N-2}}) ^{-\frac{N-4}{2}}$ (up to scalings and multiplications). △ Less

Submitted 10 March, 2024; v1 submitted 15 August, 2023; originally announced August 2023.

Comments: arXiv admin note: text overlap with arXiv:2308.06014

MSC Class: 35P30; 35J30

Classification and non-degeneracy of positive radial solutions for a weighted fourth-order equation and its application

Authors: Shengbing Deng, Xingliang Tian

Abstract: This paper is devoted to radial solutions of the following weighted fourth-order equation \begin{equation*} \mathrm{div}(|x|^α\nabla(\mathrm{div}(|x|^α\nabla u)))=u^{2^{**}_α-1},\quad u>0\quad \mbox{in}\quad \mathbb{R}^N, \end{equation*} where $N\geq 2$, $\frac{4-N}{2}<α<2$ and $2^{**}_α=\frac{2N}{N-4+2α}$. It is obvious that the solutions of above equation are invariant under the scaling… ▽ More This paper is devoted to radial solutions of the following weighted fourth-order equation \begin{equation*} \mathrm{div}(|x|^α\nabla(\mathrm{div}(|x|^α\nabla u)))=u^{2^{**}_α-1},\quad u>0\quad \mbox{in}\quad \mathbb{R}^N, \end{equation*} where $N\geq 2$, $\frac{4-N}{2}<α<2$ and $2^{**}_α=\frac{2N}{N-4+2α}$. It is obvious that the solutions of above equation are invariant under the scaling $λ^{\frac{N-4+2α}{2}}u(λx)$ while they are not invariant under translation when $α\neq 0$. We characterize all the solutions to the related linearized problem about radial solutions, and obtain the conclusion of that if $α$ satisfies $(2-α)(2N-2+α)\neq4k(N-2+k)$ for all $k\in\mathbb{N}^+$ the radial solution is non-degenerate, otherwise there exist new solutions to the linearized problem that ``replace'' the ones due to the translations invariance. As applications, firstly we investigate the remainder terms of some inequalities related to above equation. Then when $N\geq 5$ and $0<α<2$, we establish a new type second-order Caffarelli-Kohn-Nirenberg inequality \begin{equation*} \int_{\mathbb{R}^N} |\mathrm{div}(|x|^α\nabla u)|^2 \mathrm{d}x \geq C \left(\int_{\mathbb{R}^N}|u|^{2^{**}_α} \mathrm{d}x\right)^{\frac{2}{2^{**}_α}},\quad \mbox{for all}\quad u\in C^\infty_0(\mathbb{R}^N), \end{equation*} and in this case we consider a prescribed perturbation problem by using Lyapunov-Schmidt reduction. △ Less

Submitted 20 December, 2023; v1 submitted 11 August, 2023; originally announced August 2023.

Comments: This version is the final one, corresponding to the paper now published in Nonlinear Analysis DOI: https://doi.org/10.1016/j.na.2023.113468

MSC Class: 35P30; 35J30

On the stability of Caffarelli-Kohn-Nirenberg inequality in R^2

Authors: Shengbing Deng, Xingliang Tian

Abstract: Dolbeault, Esteban and Loss [Invent. Math., 2016] obtained an optimal rigidity result, that is, when $a<0$ and $b_{\mathrm{FS}}(a)\leq b<a+1$ the extremal function for best constant $\mathcal{S}_{a,b}>0$ of the following Caffarelli-Kohn-Nirenberg inequality is symmetry, \[ \mathcal{S}_{a,b}\left(\int_{\mathbb{R}^2}|x|^{-qb}|u|^q \mathrm{d}x\right)^{\frac{2}{q}} \leq \int_{\mathbb{R}^2}|x|^{-2a}|… ▽ More Dolbeault, Esteban and Loss [Invent. Math., 2016] obtained an optimal rigidity result, that is, when $a<0$ and $b_{\mathrm{FS}}(a)\leq b<a+1$ the extremal function for best constant $\mathcal{S}_{a,b}>0$ of the following Caffarelli-Kohn-Nirenberg inequality is symmetry, \[ \mathcal{S}_{a,b}\left(\int_{\mathbb{R}^2}|x|^{-qb}|u|^q \mathrm{d}x\right)^{\frac{2}{q}} \leq \int_{\mathbb{R}^2}|x|^{-2a}|\nabla u|^2 \mathrm{d}x, \quad \mbox{for all}\quad u\in C^\infty_0(\mathbb{R}^2), \] where $b_{\mathrm{FS}}(a):=a-\frac{a}{\sqrt{a^2+1}}$, $q=\frac{2}{b-a}$. An important task is investigating the stability of extremal functions set $\mathcal{M}$ for this inequality. Firstly, we classify all solutions of the linearized problem related to the extremals which fills the work of Felli and Schneider [J. Diff. Equ., 2003]. When $b_{\mathrm{FS}}(a)< b<a+1$, we investigate the stability of previous inequality by using spectral estimate combined with a compactness argument that \begin{align*} \int_{\mathbb{R}^2}|x|^{-2a}|\nabla u|^2 \mathrm{d}x -\mathcal{S}_{a,b}\left(\int_{\mathbb{R}^2}|x|^{-qb}|u|^q \mathrm{d}x\right)^{\frac{2}{q}} \geq \mathcal{B} \mathrm{dist}(u,\mathcal{M})^2,\quad \mbox{for all}\quad u\in C^\infty_0(\mathbb{R}^2), \end{align*} for some $\mathcal{B}>0$, however it is false when $b=b_{\mathrm{FS}}(a)$, which extends the work of Wei and Wu [Math. Ann., 2022] to $\mathbb{R}^2$. Furthermore, we obtain the existence of minimizers for $\mathcal{B}$ which extends the recent work of König [J. Eur. Math. Soc., to appear]. △ Less

Submitted 22 May, 2024; v1 submitted 8 August, 2023; originally announced August 2023.

Comments: arXiv admin note: text overlap with arXiv:2211.14185 by other authors

MSC Class: 35P30; 26D10; 49J40

Boundness of Some Multilinear Fractional Integral Operators in Generalized Morrey Spaces on Stratified Lie Groups

Authors: J. Wu, X. Tian

Abstract: In this paper, the main aim is to consider the Spanne-type boundedness of the multiliinear fractional integral operator $\mathcal{I}_{α,m}$ and multiliinear fractional maximal operator $\mathcal{M}_{α,m}$ in the generalized Morrey spaces over some stratified Lie group $\mathbb{G}$. In this paper, the main aim is to consider the Spanne-type boundedness of the multiliinear fractional integral operator $\mathcal{I}_{α,m}$ and multiliinear fractional maximal operator $\mathcal{M}_{α,m}$ in the generalized Morrey spaces over some stratified Lie group $\mathbb{G}$. △ Less

Submitted 17 June, 2023; originally announced June 2023.

Comments: arXiv admin note: substantial text overlap with arXiv:2206.08340

MSC Class: 42B35; 43A80

Stability of Rellich-Sobolev type inequality involving Hardy term for bi-Laplacian

Authors: Shengbing Deng, Xingliang Tian

Abstract: For $N\geq 5$ and $0<μ<N-4$, we first show a non-degenerate result of the extremal for the following Rellich-Sobolev type inequality \begin{align*} & \int_{\mathbb{R}^N}|Δu|^2 \mathrm{d}x -C_{μ,1}\int_{\mathbb{R}^N}\frac{|\nabla u|^2}{|x|^2} \mathrm{d}x +C_{μ,2}\int_{\mathbb{R}^N}\frac{u^2}{|x|^4} \mathrm{d}x \geq \mathcal{S}_μ\left(\int_{\mathbb{R}^N}|u|^{\frac{2N}{N-4}} \mathrm{d}x\rig… ▽ More For $N\geq 5$ and $0<μ<N-4$, we first show a non-degenerate result of the extremal for the following Rellich-Sobolev type inequality \begin{align*} & \int_{\mathbb{R}^N}|Δu|^2 \mathrm{d}x -C_{μ,1}\int_{\mathbb{R}^N}\frac{|\nabla u|^2}{|x|^2} \mathrm{d}x +C_{μ,2}\int_{\mathbb{R}^N}\frac{u^2}{|x|^4} \mathrm{d}x \geq \mathcal{S}_μ\left(\int_{\mathbb{R}^N}|u|^{\frac{2N}{N-4}} \mathrm{d}x\right)^\frac{N-4}{N},\quad u\in C^\infty_0(\mathbb{R}^N), \end{align*} where $C_{μ,1}$, $C_{μ,1}$ and $\mathcal{S}_μ$ are constants depending on $μ$, furthermore equality only holds for some radial functions, which is a key ingredient in analyzing the blow-up phenomena of solutions to various elliptic equations on bounded or unbounded domain. Moreover, by using spectral estimate combined with a compactness argument, we give the remainder term of above inequality. △ Less

Submitted 2 January, 2024; v1 submitted 3 June, 2023; originally announced June 2023.

Remainder terms of a nonlocal Sobolev inequality1

Authors: Shengbing Deng, Xingliang Tian, Minbo Yang, Shunneng Zhao

Abstract: In this note we study a nonlocal version of the Sobolev inequality \begin{equation*} \int_{\mathbb{R}^N}|\nabla u|^2 dx \geq S_{HLS}\left(\int_{\mathbb{R}^N}\big(|x|^{-α} \ast u^{2_α^{\ast}}\big)u^{2_α^{\ast}} dx\right)^{\frac{1}{2_α^{\ast}}}, \quad \forall u\in \mathcal{D}^{1,2}(\mathbb{R}^N), \end{equation*} where $S_{HLS}$ is the best constant, $\ast$ denotes the standard convolution and… ▽ More In this note we study a nonlocal version of the Sobolev inequality \begin{equation*} \int_{\mathbb{R}^N}|\nabla u|^2 dx \geq S_{HLS}\left(\int_{\mathbb{R}^N}\big(|x|^{-α} \ast u^{2_α^{\ast}}\big)u^{2_α^{\ast}} dx\right)^{\frac{1}{2_α^{\ast}}}, \quad \forall u\in \mathcal{D}^{1,2}(\mathbb{R}^N), \end{equation*} where $S_{HLS}$ is the best constant, $\ast$ denotes the standard convolution and $\mathcal{D}^{1,2}(\mathbb{R}^N)$ denotes the classical Sobolev space with respect to the norm $\|u\|_{\mathcal{D}^{1,2}(\mathbb{R}^N)}=\|\nabla u\|_{L^2(\mathbb{R}^N)}$. By using the nondegeneracy property of the extremal functions, we prove that the existence of the gradient type remainder term and a reminder term in the weak $L^{\frac{N}{N-2}}$-norm of above inequality for all $0<α<N$. △ Less

Submitted 26 May, 2023; originally announced May 2023.

Comments: 15 pages

MSC Class: 35P30; 35J20

The matched projections of idempotents on Hilbert $C^*$-modules

Authors: Xiaoyi Tian, Qingxiang Xu, Chunhong Fu

Abstract: The aim of this paper is to give new characterizations of some fundamental issues about idempotents. In the general setting of adjointable operators on Hilbert $C^*$-modules, a new term of quasi-projection pair is introduced. For each idempotent $Q$, a projection $m(Q)$, called the matched projection of $Q$, is constructed. It is shown that $Q$ and $m(Q)$ as idempotents are homotopic, and… ▽ More The aim of this paper is to give new characterizations of some fundamental issues about idempotents. In the general setting of adjointable operators on Hilbert $C^*$-modules, a new term of quasi-projection pair is introduced. For each idempotent $Q$, a projection $m(Q)$, called the matched projection of $Q$, is constructed. It is shown that $Q$ and $m(Q)$ as idempotents are homotopic, and $\big(m(Q),Q\big)$ is a quasi-projection pair. Some formulas for $m(Q)$ are derived. Based on these formulas, representations and norm estimations associated with $m(Q)$ are dealt with. △ Less

Submitted 2 April, 2024; v1 submitted 22 May, 2023; originally announced May 2023.

Comments: The original submission entitled "Quasi-projection pairs on Hilbert $C^*$-modules" will be divided into three parts, and this is the first part. Some minor changes have been made in this revision, which is accepted for publication in Journal of Operator Theory

MSC Class: 46L08; 47A05

Stability of Hardy-Sobolev inequality involving p-Laplace

Authors: Shengbing Deng, Xingliang Tian

Abstract: This paper is devoted to considering the following Hardy-Sobolev inequality \[ \int_{\mathbb{R}^N}|\nabla u|^p \mathrm{d}x \geq \mathcal{S}_β\left(\int_{\mathbb{R}^N}\frac{|u|^{p^*_β}}{|x|^β} \mathrm{d}x\right)^\frac{p}{p^*_β},\quad \forall u\in C^\infty_0(\mathbb{R}^N), \] for some constant $\mathcal{S}_β>0$, where $1<p<N$, $0\leq β<p$, $p^*_β=\frac{p(N-β)}{N-p}$. Firstly, since this problem invo… ▽ More This paper is devoted to considering the following Hardy-Sobolev inequality \[ \int_{\mathbb{R}^N}|\nabla u|^p \mathrm{d}x \geq \mathcal{S}_β\left(\int_{\mathbb{R}^N}\frac{|u|^{p^*_β}}{|x|^β} \mathrm{d}x\right)^\frac{p}{p^*_β},\quad \forall u\in C^\infty_0(\mathbb{R}^N), \] for some constant $\mathcal{S}_β>0$, where $1<p<N$, $0\leq β<p$, $p^*_β=\frac{p(N-β)}{N-p}$. Firstly, since this problem involves quasilinear operator, we need to establish a compact embedding theorem for some suitable weighted spaces. Moreover, due to the Hardy term $|x|^{-β}$, some new estimates are established. Based on those works, we give the classification to the linearized problem related to the extremals which has its own interest such as in blow-up analysis. Then we investigate the gradient stability of above inequality by using spectral estimate combined with a compactness argument, which extends the work of Figalli and Zhang (Duke Math. J., 2022) to a weighted case. △ Less

Submitted 18 January, 2023; originally announced January 2023.

Comments: 49 pages. arXiv admin note: text overlap with arXiv:2212.05459; text overlap with arXiv:2003.04037 by other authors

Nonlocal half-ball vector operators on bounded domains: Poincaré inequality and its applications

Authors: Zhaolong Han, Xiaochuan Tian

Abstract: This work contributes to nonlocal vector calculus as an indispensable mathematical tool for the study of nonlocal models that arises in a variety of applications. We define the nonlocal half-ball gradient, divergence and curl operators with general kernel functions (integrable or fractional type with finite or infinite supports) and study the associated nonlocal vector identities. We study the non… ▽ More This work contributes to nonlocal vector calculus as an indispensable mathematical tool for the study of nonlocal models that arises in a variety of applications. We define the nonlocal half-ball gradient, divergence and curl operators with general kernel functions (integrable or fractional type with finite or infinite supports) and study the associated nonlocal vector identities. We study the nonlocal function space on bounded domains associated with zero Dirichlet boundary conditions and the half-ball gradient operator and show it is a separable Hilbert space with smooth functions dense in it. A major result is the nonlocal Poincaré inequality, based on which a few applications are discussed, and these include applications to nonlocal convection-diffusion, nonlocal correspondence model of linear elasticity, and nonlocal Helmholtz decomposition on bounded domains. △ Less

Submitted 16 March, 2024; v1 submitted 28 December, 2022; originally announced December 2022.

MSC Class: 45A05; 26B12; 47G10; 47G30; 35A23; 46E35; 74G65

Caffarelli-Kohn-Nirenberg-type inequalities related to weighted $p$-Laplace equations

Authors: Shengbing Deng, Xingliang Tian

Abstract: We use a suitable transform related to Sobolev inequality to investigate the sharp constants and optimizers for some Caffarelli-Kohn-Nirenberg-type inequalities which are related to the weighted $p$-Laplace equations. Moreover, we give the classification to the linearized problem related to the radial extremals. As an application, we investigate the gradient type remainder term of related inequali… ▽ More We use a suitable transform related to Sobolev inequality to investigate the sharp constants and optimizers for some Caffarelli-Kohn-Nirenberg-type inequalities which are related to the weighted $p$-Laplace equations. Moreover, we give the classification to the linearized problem related to the radial extremals. As an application, we investigate the gradient type remainder term of related inequality by using spectral estimate combined with a compactness argument which extends the work of Figalli and Zhang (Duke Math. J. 2022) at least for radial case. △ Less

Submitted 11 December, 2022; originally announced December 2022.

Comments: 39 pages, All comments are welcome

Monotone meshfree methods for linear elliptic equations in non-divergence form via nonlocal relaxation

Authors: Qihao Ye, Xiaochuan Tian

Abstract: We design a monotone meshfree finite difference method for linear elliptic equations in the non-divergence form on point clouds via a nonlocal relaxation method. The key idea is a novel combination of a nonlocal integral relaxation of the PDE problem with a robust meshfree discretization on point clouds. Minimal positive stencils are obtained through a local $l_1$-type optimization procedure that… ▽ More We design a monotone meshfree finite difference method for linear elliptic equations in the non-divergence form on point clouds via a nonlocal relaxation method. The key idea is a novel combination of a nonlocal integral relaxation of the PDE problem with a robust meshfree discretization on point clouds. Minimal positive stencils are obtained through a local $l_1$-type optimization procedure that automatically guarantees the stability and, therefore, the convergence of the meshfree discretization for linear elliptic equations. A major theoretical contribution is the existence of consistent and positive stencils for a given point cloud geometry. We provide sufficient conditions for the existence of positive stencils by finding neighbors within an ellipse (2d) or ellipsoid (3d) surrounding each interior point, generalizing the study for Poisson's equation by Seibold (Comput Methods Appl Mech Eng 198(3-4):592-601, 2008). It is well-known that wide stencils are in general needed for constructing consistent and monotone finite difference schemes for linear elliptic equations. Our result represents a significant improvement in the stencil width estimate for positive-type finite difference methods for linear elliptic equations in the near-degenerate regime (when the ellipticity constant becomes small), compared to previously known works in this area. Numerical algorithms and practical guidance are provided with an eye on the case of small ellipticity constant. At the end, we present numerical results for the performance of our method in both 2d and 3d, examining a range of ellipticity constants including the near-degenerate regime. △ Less

Submitted 5 August, 2023; v1 submitted 22 November, 2022; originally announced November 2022.

Comments: 36 pages, 16 figures

MSC Class: 65N06; 65N12; 65N35; 35B50; 35J15 (Primary) 65R20; 35J70; 45A05 (Secondary)

Anderson acceleration of gradient methods with energy for optimization problems

Authors: Hailiang Liu, Jia-Hao He, Xuping Tian

Abstract: Anderson acceleration (AA) as an efficient technique for speeding up the convergence of fixed-point iterations may be designed for accelerating an optimization method. We propose a novel optimization algorithm by adapting Anderson acceleration to the energy adaptive gradient method (AEGD) [arXiv:2010.05109]. The feasibility of our algorithm is examined in light of convergence results for AEGD, tho… ▽ More Anderson acceleration (AA) as an efficient technique for speeding up the convergence of fixed-point iterations may be designed for accelerating an optimization method. We propose a novel optimization algorithm by adapting Anderson acceleration to the energy adaptive gradient method (AEGD) [arXiv:2010.05109]. The feasibility of our algorithm is examined in light of convergence results for AEGD, though it is not a fixed-point iteration. We also quantify the accelerated convergence rate of AA for gradient descent by a factor of the gain at each implementation of the Anderson mixing. Our experimental results show that the proposed algorithm requires little tuning of hyperparameters and exhibits superior fast convergence. △ Less

Submitted 15 November, 2022; originally announced November 2022.

Comments: 18 pages, 4 figures

MSC Class: 65K10; 68Q25

Conditional intermediate entropy and Birkhoff average properties of hyperbolic flows

Authors: Xiaobo Hou, Xueting Tian

Abstract: Katok conjectured that every $C^{2}$ diffeomorphism $f$ on a Riemannian manifold has the intermediate entropy property, that is, for any constant $c \in[0, h_{top}(f))$, there exists an ergodic measure $μ$ of $f$ satisfying $h_μ(f)=c$. In this paper we consider a conditional intermediate metric entropy property and two conditional intermediate Birkhoff average properties for flows. For a basic set… ▽ More Katok conjectured that every $C^{2}$ diffeomorphism $f$ on a Riemannian manifold has the intermediate entropy property, that is, for any constant $c \in[0, h_{top}(f))$, there exists an ergodic measure $μ$ of $f$ satisfying $h_μ(f)=c$. In this paper we consider a conditional intermediate metric entropy property and two conditional intermediate Birkhoff average properties for flows. For a basic set $Λ$ of a flow $Φ$ and two continuous function $g,$ $h$ on $Λ,$ we obtain $$\mathrm{Int}\left\{h_μ(Φ):μ\in \mathcal{M}_{erg}(Φ,Λ)\text{ and }\int g dμ=α\right\}=\mathrm{Int}\left\{h_μ(Φ):μ\in \mathcal{M}(Φ,Λ) \text{ and }\int g dμ=α\right\},$$ $$\mathrm{Int}\left\{\int g dμ:μ\in \mathcal{M}_{erg}(Φ,Λ)\text{ and }h_μ(Φ)=c\right\}=\mathrm{Int}\left\{\int g dμ:μ\in \mathcal{M}(Φ,Λ) \text{ and }h_μ(Φ)=c\right\}$$ and $$\mathrm{Int}\left\{\int h dμ:μ\in \mathcal{M}_{erg}(Φ,Λ)\text{ and }\int g dμ=α\right\}=\mathrm{Int}\left\{\int h dμ:μ\in \mathcal{M}(Φ,Λ) \text{ and }\int g dμ=α\right\}$$ for any $α\in \left(\inf_{μ\in \in \mathcal{M}(Φ,Λ) }\int g dμ, \, \sup_{μ\in \in \mathcal{M}(Φ,Λ) }\int g dμ\right)$ and any $c\in (0,h_{top}(Λ)).$ In this process, we establish 'multi-horseshoe' entropy-dense property and use it to get the goal combined with conditional variational principles. We also obtain same result for singular hyperbolic attractors. △ Less

Submitted 2 November, 2022; v1 submitted 7 September, 2022; originally announced September 2022.

SGEM: stochastic gradient with energy and momentum

Authors: Hailiang Liu, Xuping Tian

Abstract: In this paper, we propose SGEM, Stochastic Gradient with Energy and Momentum, to solve a large class of general non-convex stochastic optimization problems, based on the AEGD method that originated in the work [AEGD: Adaptive Gradient Descent with Energy. arXiv: 2010.05109]. SGEM incorporates both energy and momentum at the same time so as to inherit their dual advantages. We show that SGEM featur… ▽ More In this paper, we propose SGEM, Stochastic Gradient with Energy and Momentum, to solve a large class of general non-convex stochastic optimization problems, based on the AEGD method that originated in the work [AEGD: Adaptive Gradient Descent with Energy. arXiv: 2010.05109]. SGEM incorporates both energy and momentum at the same time so as to inherit their dual advantages. We show that SGEM features an unconditional energy stability property, and derive energy-dependent convergence rates in the general nonconvex stochastic setting, as well as a regret bound in the online convex setting. A lower threshold for the energy variable is also provided. Our experimental results show that SGEM converges faster than AEGD and generalizes better or at least as well as SGDM in training some deep neural networks. △ Less

Submitted 3 August, 2022; originally announced August 2022.

Comments: 24 pages, 4 figures

MSC Class: 90C15; 65K10; 68Q25

Abundance of Smale's horseshoes and ergodic measures via multifractal analysis and various quantitative spectrums

Authors: Yiwei Dong, Xiaobo Hou, Xueting Tian

Abstract: In this article, we combine the perspectives of density, entropy, and multifractal analysis to investigate the structure of ergodic measures. We prove that for each transitive topologically Anosov system $(X,f)$, each continuous function $\varphi$ on $X$ and each $(a,h)\in \mathrm{Int}\{(\int \varphi dμ, h_μ(f)):μ\in M_f(X)\},$ the set $\{μ\in M_f^e(X): (\int \varphi dμ, h_μ(f))=(a,h)\}$ is non-em… ▽ More In this article, we combine the perspectives of density, entropy, and multifractal analysis to investigate the structure of ergodic measures. We prove that for each transitive topologically Anosov system $(X,f)$, each continuous function $\varphi$ on $X$ and each $(a,h)\in \mathrm{Int}\{(\int \varphi dμ, h_μ(f)):μ\in M_f(X)\},$ the set $\{μ\in M_f^e(X): (\int \varphi dμ, h_μ(f))=(a,h)\}$ is non-empty and contains a dense $G_δ$ subset of $\{μ\in M_f(X): (\int \varphi dμ, h_μ(f))=(a,h)\}.$ Meanwhile, combining the development of non-hyperbolic systems and cocycles we give a general framework and use it to obtain intermediate entropy property of ergodic measures with same Lyapunov exponent for non-hyperbolic step skew-products, elliptic $\operatorname{SL}(2, \mathbb{R})$ cocycles and robustly non-hyperbolic transitive diffeomorphisms. Moreover, we get generalized results on multiple functions and use them to obtain the intermediate Hausdorff dimension of ergodic measures for transitive average conformal or quasi-conformal Anosov diffeomorphisms, that is $\left\{\operatorname{dim}_H μ: μ\in M_f^e(M)\right\}= \left\{\operatorname{dim}_H μ: μ\in M_f(M)\right\}.$ In this process, we introduce and establish a 'multi-horseshoe' entropy-dense property and use it to get the goal combined with the well-known conditional variational principles. As applications, we also obtain many new observations on various other quantitative spectrums including Lyapunov exponents, first return rate, geometric pressure, unstable Hausdorff dimension, etc. △ Less

Submitted 20 February, 2024; v1 submitted 5 July, 2022; originally announced July 2022.

Strongly distributional chaos in the sets of twelve different types of non-recurrent points

Authors: An Chen, Xiaobo Hou, Wanshan Lin, Xueting Tian

Abstract: In present paper we mainly focus on non-recurrent dynamical orbits with empty syndetic center and show that twelve different statistical structures over mixing expanding maps or transitive Anosov diffeomorphisms all have dynamical complexity in the sense of strongly distributional chaos. In present paper we mainly focus on non-recurrent dynamical orbits with empty syndetic center and show that twelve different statistical structures over mixing expanding maps or transitive Anosov diffeomorphisms all have dynamical complexity in the sense of strongly distributional chaos. △ Less

Submitted 26 May, 2022; originally announced May 2022.

Some weighted fourth-order Hardy-Henon equations

Authors: Shengbing Deng, Xingliang Tian

Abstract: By using a suitable transform related to Sobolev inequality, we investigate the sharp constants and optimizers in radial space for the following weighted Caffarelli-Kohn-Nirenberg-type inequalities: \begin{equation*} \int_{\mathbb{R}^N}|x|^α|Δu|^2 dx \geq S^{rad}(N,α)\left(\int_{\mathbb{R}^N}|x|^{-α}|u|^{p^*_α} dx\right)^{\frac{2}{p^*_α}}, \quad u\in C^\infty_c(\mathbb{R}^N), \end{equation*} where… ▽ More By using a suitable transform related to Sobolev inequality, we investigate the sharp constants and optimizers in radial space for the following weighted Caffarelli-Kohn-Nirenberg-type inequalities: \begin{equation*} \int_{\mathbb{R}^N}|x|^α|Δu|^2 dx \geq S^{rad}(N,α)\left(\int_{\mathbb{R}^N}|x|^{-α}|u|^{p^*_α} dx\right)^{\frac{2}{p^*_α}}, \quad u\in C^\infty_c(\mathbb{R}^N), \end{equation*} where $N\geq 3$, $4-N<α<2$, $p^*_α=\frac{2(N-α)}{N-4+α}$. Then we obtain the explicit form of the unique (up to scaling) radial positive solution $U_{λ,α}$ to the weighted fourth-order Hardy (for $α>0$) or Hénon (for $α<0$) equation: \begin{equation*} Δ(|x|^αΔu)=|x|^{-α} u^{p^*_α-1},\quad u>0 \quad \mbox{in}\quad \mathbb{R}^N. \end{equation*} %Furthermore, we characterize all the solutions to the linearized problem related to above equation at $U_{1,α}$. For $α\neq 0$, it is known the solutions of above equation are invariant for dilations $λ^{\frac{N-4+α}{2}}u(λx)$ but not for translations. However we show that if $α$ is an even integer, there exist new solutions to the linearized problem, which related to above equation at $U_{1,α}$, that "replace" the ones due to the translations invariance. This interesting phenomenon was first shown by Gladiali, Grossi and Neves [Adv. Math. 249, 2013, 1-36] for the second-order Hénon problem. Finally, as applications, we investigate the reminder term of above inequality and also the existence of solutions to some related perturbed equations. △ Less

Submitted 2 November, 2022; v1 submitted 18 May, 2022; originally announced May 2022.

Journal ref: Journal of Functional Analysis,284(2023)109745

Dynamic behavior for a gradient algorithm with energy and momentum

Authors: Hailiang Liu, Xuping Tian

Abstract: This paper investigates a novel gradient algorithm, AGEM, using both energy and momentum, for addressing general non-convex optimization problems. The solution properties of the AGEM algorithm, including aspects such as uniformly boundedness and convergence to critical points, are examined. The dynamic behavior is studied through a comprehensive analysis of a high-resolution ODE system. This ODE s… ▽ More This paper investigates a novel gradient algorithm, AGEM, using both energy and momentum, for addressing general non-convex optimization problems. The solution properties of the AGEM algorithm, including aspects such as uniformly boundedness and convergence to critical points, are examined. The dynamic behavior is studied through a comprehensive analysis of a high-resolution ODE system. This ODE system, being nonlinear, is derived by taking the limit of the discrete scheme while preserving the momentum effect through a rescaling of the momentum parameter. The paper emphasizes the global well-posedness of the ODE system and the time-asymptotic convergence of solution trajectories. Furthermore, we establish a linear convergence rate for objective functions that adhere to the Polyak-Łojasiewicz condition. △ Less

Submitted 12 March, 2024; v1 submitted 23 March, 2022; originally announced March 2022.

Comments: 30 pages, 2 figures

MSC Class: 65K10 (Primary) 90C15 (Secondary)

An Adaptive Gradient Method with Energy and Momentum

Authors: Hailiang Liu, Xuping Tian

Abstract: We introduce a novel algorithm for gradient-based optimization of stochastic objective functions. The method may be seen as a variant of SGD with momentum equipped with an adaptive learning rate automatically adjusted by an 'energy' variable. The method is simple to implement, computationally efficient, and well suited for large-scale machine learning problems. The method exhibits unconditional en… ▽ More We introduce a novel algorithm for gradient-based optimization of stochastic objective functions. The method may be seen as a variant of SGD with momentum equipped with an adaptive learning rate automatically adjusted by an 'energy' variable. The method is simple to implement, computationally efficient, and well suited for large-scale machine learning problems. The method exhibits unconditional energy stability for any size of the base learning rate. We provide a regret bound on the convergence rate under the online convex optimization framework. We also establish the energy-dependent convergence rate of the algorithm to a stationary point in the stochastic non-convex setting. In addition, a sufficient condition is provided to guarantee a positive lower threshold for the energy variable. Our experiments demonstrate that the algorithm converges fast while generalizing better than or as well as SGD with momentum in training deep neural networks, and compares also favorably to Adam. △ Less

Submitted 23 March, 2022; originally announced March 2022.

Comments: 29 pages, 6 figures

MSC Class: 65K10 (Primary) 90C15; 68Q25 (Secondary)

Journal ref: Ann. Appl. Math. 38(2022)1-40

Nonlocal diffusion models with consistent local and fractional limits

Authors: Qiang Du, Xiaochuan Tian, Zhi Zhou

Abstract: For some spatially nonlocal diffusion models with a finite range of nonlocal interactions measured by a positive parameter $δ$, we review their formulation defined on a bounded domain subject to various conditions that correspond to some inhomogeneous data. We consider their consistency to similar inhomogeneous boundary value problems of classical partial differential equation (PDE) models as the… ▽ More For some spatially nonlocal diffusion models with a finite range of nonlocal interactions measured by a positive parameter $δ$, we review their formulation defined on a bounded domain subject to various conditions that correspond to some inhomogeneous data. We consider their consistency to similar inhomogeneous boundary value problems of classical partial differential equation (PDE) models as the nonlocal interaction kernel gets localized in the local $δ\to 0$ limit, and at the same time, for rescaled fractional type kernels, to corresponding inhomogeneous nonlocal boundary value problems of fractional equations in the global $δ\to \infty$ limit. Such discussions help to delineate issues related to nonlocal problems defined on a bounded domain with inhomogeneous data. △ Less

Submitted 24 December, 2022; v1 submitted 28 February, 2022; originally announced March 2022.

Comments: 31 pages, 5 figures

MSC Class: 47G10; 46E35; 35R11; 35J25

A Meshfree Peridynamic Model for Brittle Fracture in Randomly Heterogeneous Materials

Authors: Yiming Fan, Huaiqian You, Xiaochuan Tian, Xiu Yang, Xingjie Li, Naveen Prakash, Yue Yu

Abstract: In this work we aim to develop a unified mathematical framework and a reliable computational approach to model the brittle fracture in heterogeneous materials with variability in material microstructures, and to provide statistic metrics for quantities of interest, such as the fracture toughness. To depict the material responses and naturally describe the nucleation and growth of fractures, we con… ▽ More In this work we aim to develop a unified mathematical framework and a reliable computational approach to model the brittle fracture in heterogeneous materials with variability in material microstructures, and to provide statistic metrics for quantities of interest, such as the fracture toughness. To depict the material responses and naturally describe the nucleation and growth of fractures, we consider the peridynamics model. In particular, a stochastic state-based peridynamic model is developed, where the micromechanical parameters are modeled by a finite-dimensional random vector, or a combination of random variables truncating the Karhunen-Loève decomposition or the principle component analysis (PCA). To solve this stochastic peridynamic problem, probabilistic collocation method (PCM) is employed to sample the random field representing the micromechanical parameters. For each sample, the deterministic peridynamic problem is discretized with an optimization-based meshfree quadrature rule. We present rigorous analysis for the proposed scheme and demonstrate its convergence for a number of benchmark problems, showing that it sustains the asymptotic compatibility spatially and achieves an algebraic or sub-exponential convergence rate in the random space as the number of collocation points grows. Finally, to validate the applicability of this approach on real-world fracture problems, we consider the problem of crystallization toughening in glass-ceramic materials, in which the material at the microstructural scale contains both amorphous glass and crystalline phases. The proposed stochastic peridynamic solver is employed to capture the crack initiation and growth for glass-ceramics with different crystal volume fractions, and the averaged fracture toughness are calculated. The numerical estimates of fracture toughness show good consistency with experimental measurements. △ Less

Submitted 14 February, 2022; originally announced February 2022.

Efficient optimization-based quadrature for variational discretization of nonlocal problems

Authors: Marco Pasetto, Zhaoxiang Shen, Marta D'Elia, Xiaochuan Tian, Nathaniel Trask, David Kamensky

Abstract: Casting nonlocal problems in variational form and discretizing them with the finite element (FE) method facilitates the use of nonlocal vector calculus to prove well-posedeness, convergence, and stability of such schemes. Employing an FE method also facilitates meshing of complicated domain geometries and coupling with FE methods for local problems. However, nonlocal weak problems involve the comp… ▽ More Casting nonlocal problems in variational form and discretizing them with the finite element (FE) method facilitates the use of nonlocal vector calculus to prove well-posedeness, convergence, and stability of such schemes. Employing an FE method also facilitates meshing of complicated domain geometries and coupling with FE methods for local problems. However, nonlocal weak problems involve the computation of a double-integral, which is computationally expensive and presents several challenges. In particular, the inner integral of the variational form associated with the stiffness matrix is defined over the intersections of FE mesh elements with a ball of radius $δ$, where $δ$ is the range of nonlocal interaction. Identifying and parameterizing these intersections is a nontrivial computational geometry problem. In this work, we propose a quadrature technique where the inner integration is performed using quadrature points distributed over the full ball, without regard for how it intersects elements, and weights are computed based on the generalized moving least squares method. Thus, as opposed to all previously employed methods, our technique does not require element-by-element integration and fully circumvents the computation of element-ball intersections. This paper considers one- and two-dimensional implementations of piecewise linear continuous FE approximations, focusing on the case where the element size h and the nonlocal radius $δ$ are proportional, as is typical of practical computations. When boundary conditions are treated carefully and the outer integral of the variational form is computed accurately, the proposed method is asymptotically compatible in the limit of $h \sim δ\to 0$, featuring at least first-order convergence in L^2 for all dimensions, using both uniform and nonuniform grids. △ Less

Submitted 28 January, 2022; originally announced January 2022.

Comments: 59 pages, 21 figures

Report number: SAND2022-0834, Sandia National Laboratories

Journal ref: Computer Methods in Applied Mechanics and Engineering, 396 (2022), 115104

Boundary concentration of peak solutions for fractional Schrödinger-Poisson system

Authors: Shengbing Deng, Xingliang Tian

Abstract: The goal of this paper is to study the existence of peak solutions for the following fractional Schrödinger-Poisson system: \begin{eqnarray*} \left\{ \arraycolsep=1.5pt \begin{array}{ll} \varepsilon^{2s}(-Δ)^{s}u+u+φu=u^p,\ \ \ &\ \mbox{in}\ Ω,\\[2mm] (-Δ)^{s}φ=u^2,\ \ \ &\ \mbox{in}\ Ω,\\[2mm] u=φ=0,\ \ \ \ &\ \mbox{in}\ \mathbb{R}^N\setminus Ω, \end{array} \right. \end{eqnarray… ▽ More The goal of this paper is to study the existence of peak solutions for the following fractional Schrödinger-Poisson system: \begin{eqnarray*} \left\{ \arraycolsep=1.5pt \begin{array}{ll} \varepsilon^{2s}(-Δ)^{s}u+u+φu=u^p,\ \ \ &\ \mbox{in}\ Ω,\\[2mm] (-Δ)^{s}φ=u^2,\ \ \ &\ \mbox{in}\ Ω,\\[2mm] u=φ=0,\ \ \ \ &\ \mbox{in}\ \mathbb{R}^N\setminus Ω, \end{array} \right. \end{eqnarray*} where $s\in(0,1)$, $N>2s$, $p\in (1,\frac{N+2s}{N-2s})$, $Ω$ is a bounded domain in $\mathbb{R}^N$ with Lipschitz boundary, and $(-Δ)^{s}$ is the fractional Laplacian operator, $\varepsilon$ is a small positive parameter. By using the Lyapunov-Schmidt reduction method, we construct a single peak solution $(u_\varepsilon,φ_\varepsilon)$ such that the peak of $u_\varepsilon$ is in the domain but near the boundary. In order to characterize the boundary concentration of solutions, which concentrates at an approximate distance $\varepsilon^{2/3}$ away from the boundary $\partialΩ$ as $\varepsilon$ tends to 0, some new estimates and analytic technique are used. △ Less

Submitted 17 January, 2022; originally announced January 2022.

Existence and multiplicity of solutions to a Kirchhoff type elliptic system with Trudinger-Moser growth

Authors: Shengbing Deng, Xingliang Tian

Abstract: This paper deals with the existence and multiplicity of solutions for a class of Kirchhoff type elliptic system involving the Trudinger-Moser exponential growth nonlinearities. We first study the existence of solutions for the following system \begin{eqnarray*} \left\{ \arraycolsep=1.5pt \begin{array}{ll} -\big(a_1+b_1\|u\|^{2(θ_1-1)}\big)Δu= λH_u(x,u,v)\ \ \ &\ \mbox{in}\ \ \ Ω,\\[2mm] -\… ▽ More This paper deals with the existence and multiplicity of solutions for a class of Kirchhoff type elliptic system involving the Trudinger-Moser exponential growth nonlinearities. We first study the existence of solutions for the following system \begin{eqnarray*} \left\{ \arraycolsep=1.5pt \begin{array}{ll} -\big(a_1+b_1\|u\|^{2(θ_1-1)}\big)Δu= λH_u(x,u,v)\ \ \ &\ \mbox{in}\ \ \ Ω,\\[2mm] -\big(a_2+b_2\|v\|^{2(θ_2-1)}\big)Δv= λH_v(x,u,v)\ \ \ &\ \mbox{in}\ \ \ Ω,\\[2mm] u=0, v=0\ \ \ \ &\ \mbox{on}\ \ \ \partialΩ, \end{array} \right. \end{eqnarray*} where $Ω$ is a bounded domain in $\mathbb{R}^2$ with smooth boundary,\ $\|u\|=\big(\int_Ω|\nabla u|^2dx\big)^{1/2}$, $H_u$ and $H_v$ behave like $e^{β|s|^2}$ when $|s|\rightarrow \infty$ for some $β>0$, $a_1,\ a_2>0$, $b_1,\ b_2> 0$, $θ_1,\ θ_2> 1$ and $λ$ is a positive parameter. In the later part of the paper, we also discuss a new multiplicity result for the above system with a positive parameter induced by the nonlocal dependence. The Kirchhoff term and the lack of compactness of the associated energy functional due to the Trudinger-Moser embedding have to be overcome via some new techniques. △ Less

Submitted 5 October, 2022; v1 submitted 7 January, 2022; originally announced January 2022.

On a nonhomogeneous Kirchhoff type elliptic system with the singular Trudinger-Moser growth

Authors: Shengbing Deng, Xingliang Tian

Abstract: The aim of this paper is to study the multiplicity of solutions for the following Kirchhoff type elliptic systems \begin{eqnarray*} \left\{ \arraycolsep=1.5pt \begin{array}{ll} -m\left(\sum^k_{j=1}\|u_j\|^2\right)Δu_i=\frac{f_i(x,u_1,\ldots,u_k)}{|x|^β}+\varepsilon h_i(x),\ \ & \mbox{in}\ \ Ω, \ \ i=1,\ldots,k ,\\[2mm] u_1=u_2=\cdots=u_k=0,\ \ & \mbox{on}\ \ \partialΩ, \end{array} \r… ▽ More The aim of this paper is to study the multiplicity of solutions for the following Kirchhoff type elliptic systems \begin{eqnarray*} \left\{ \arraycolsep=1.5pt \begin{array}{ll} -m\left(\sum^k_{j=1}\|u_j\|^2\right)Δu_i=\frac{f_i(x,u_1,\ldots,u_k)}{|x|^β}+\varepsilon h_i(x),\ \ & \mbox{in}\ \ Ω, \ \ i=1,\ldots,k ,\\[2mm] u_1=u_2=\cdots=u_k=0,\ \ & \mbox{on}\ \ \partialΩ, \end{array} \right. \end{eqnarray*} where $Ω$ is a bounded domain in $\mathbb{R}^2$ containing the origin with smooth boundary, $β\in [0,2)$, $m$ is a Kirchhoff type function, $\|u_j\|^2=\int_Ω|\nabla u_j|^2dx$, $f_i$ behaves like $e^{βs^2}$ when $|s|\rightarrow \infty$ for some $β>0$, and there is $C^1$ function $F: Ω\times\mathbb{R}^k\to \mathbb{R}$ such that $\left(\frac{\partial F}{\partial u_1},\ldots,\frac{\partial F}{\partial u_k}\right)=\left(f_1,\ldots,f_k\right)$, $h_i\in \left(\big(H^1_0(Ω)\big)^*,\|\cdot\|_*\right)$. We establish sufficient conditions for the multiplicity of solutions of the above system by using variational methods with a suitable singular Trudinger-Moser inequality when $\varepsilon>0$ is small. △ Less

Submitted 7 January, 2022; originally announced January 2022.

Ergodic Optimization Restricted On Certain Subsets Of Invariant Measures

Authors: Wanshan Lin, Xueting Tian

Abstract: In this article, we pay attention to transitive dynamical systems having the shadowing property and the entropy functions are upper semicontinuous. As for these dynamical systems, when we consider ergodic optimization restricted on the subset of invariant measures whose metric entropy are equal or greater than a given constant, we prove that for generic real continuous functions the ergodic optimi… ▽ More In this article, we pay attention to transitive dynamical systems having the shadowing property and the entropy functions are upper semicontinuous. As for these dynamical systems, when we consider ergodic optimization restricted on the subset of invariant measures whose metric entropy are equal or greater than a given constant, we prove that for generic real continuous functions the ergodic optimization measure is unique, ergodic, full support and have metric entropy equal to the given constant. Similar results also hold for suspension flows over transitive subshift of finite type, Cr (r\geq 2)- generic geometric Lorenz attractors and C1-generic singular hyperbolic attractors. △ Less

Submitted 23 December, 2021; originally announced December 2021.

Uniqueness Of Ergodic Optimization Of Top Lyapunov Exponent For Typical Matrix Cocycles

Authors: Wanshan Lin, Xueting Tian

Abstract: In this article, we consider the ergodic optimization of the top Lyapunov exponent. We prove that there is a unique maximising measure of top Lyapunov expoent for typical matrix cocyles. By using the results we obtain, we prove that in any non-uniquely ergodic minimal dynamical system, the Lyapunov-irregular points are typical for typical matrix cocyles. In this article, we consider the ergodic optimization of the top Lyapunov exponent. We prove that there is a unique maximising measure of top Lyapunov expoent for typical matrix cocyles. By using the results we obtain, we prove that in any non-uniquely ergodic minimal dynamical system, the Lyapunov-irregular points are typical for typical matrix cocyles. △ Less

Submitted 23 December, 2021; originally announced December 2021.

Entropy of irregular points that are not uniformly hyperbolic

Authors: Xiaobo Hou, Xueting Tian

Abstract: In this article we prove that for a $C^{1+α}$ diffeomorphism on a compact Riemannian manifold, if there is a hyperbolic ergodic measure whose support is not uniformly hyperbolic, then the topological entropy of the set of irregular points that are not uniformly hyperbolic is larger than or equal to the metric entropy of the hyperbolic ergodic measure. In the process of proof, we give an abstract g… ▽ More In this article we prove that for a $C^{1+α}$ diffeomorphism on a compact Riemannian manifold, if there is a hyperbolic ergodic measure whose support is not uniformly hyperbolic, then the topological entropy of the set of irregular points that are not uniformly hyperbolic is larger than or equal to the metric entropy of the hyperbolic ergodic measure. In the process of proof, we give an abstract general mechanism to study topological entropy of irregular points provided that the system has a sequence of nondecreasing invariant compact subsets such that every subsystem has shadowing property and is transitive. △ Less

Submitted 16 November, 2021; originally announced November 2021.

Comments: 24 pages. arXiv admin note: substantial text overlap with arXiv:2111.06055

Strongly distributional chaos of irregular orbits that are not uniformly hyperbolic

Authors: Xiaobo Hou, Xueting Tian

Abstract: In this article we prove that for a diffeomorphism on a compact Riemannian manifold, if there is a nontrival homoclinic class that is not uniformly hyperbolic or the diffeomorphism is a $C^{1+α}$ and there is a hyperbolic ergodic measure whose support is not uniformly hyperbolic, then we find a type of strongly distributional chaos which is stronger than usual distributional chaos and Li-Yorke cha… ▽ More In this article we prove that for a diffeomorphism on a compact Riemannian manifold, if there is a nontrival homoclinic class that is not uniformly hyperbolic or the diffeomorphism is a $C^{1+α}$ and there is a hyperbolic ergodic measure whose support is not uniformly hyperbolic, then we find a type of strongly distributional chaos which is stronger than usual distributional chaos and Li-Yorke chaos in the set of irregular orbits that are not uniformly hyperbolic. Meanwhile, we prove that various fractal sets are strongly distributional chaotic, such as irregular sets, level sets, several recurrent level sets of points with different recurrent frequency, and some intersections of these fractal sets. In the process of proof, we give an abstract general mechanism to study strongly distributional chaos provided that the system has a sequence of nondecreasing invariant compact subsets such that every subsystem has exponential specification property, or has exponential shadowing property and transitivity. The advantage of this abstract framework is that it is not only applicable in systems with specification property including transitive Anosov diffeomorphisms, mixing expanding maps, mixing subshifts of finite type and mixing sofic subshifts but also applicable in systems without specification property including $β$-shifts, Katok map, generic systems in the space of robustly transitive diffeomorphisms and generic volume-preserving diffeomorphisms. △ Less

Submitted 11 November, 2021; originally announced November 2021.

Comments: 40 pages

On multifractal analysis and large deviations of singular-hyperbolic attractors

Authors: Yi Shi, Xueting Tian, Paulo Varandas, Xiaodong Wang

Abstract: In this paper we study the multifractal analysis and large derivations for singular hyperbolic attractors, including the geometric Lorenz attractors. For each singular hyperbolic homoclinic class whose periodic orbits are all homoclinically related and such that the space of ergodic probability measures is connected, we prove that: (i) level sets associated to continuous observables are dense in t… ▽ More In this paper we study the multifractal analysis and large derivations for singular hyperbolic attractors, including the geometric Lorenz attractors. For each singular hyperbolic homoclinic class whose periodic orbits are all homoclinically related and such that the space of ergodic probability measures is connected, we prove that: (i) level sets associated to continuous observables are dense in the homoclinic class and satisfy a variational principle; (ii) irregular sets are either empty or are Baire generic and carry full topological entropy. The assumptions are satisfied by $C^1$-generic singular hyperbolic attractors and $C^r$-generic geometric Lorenz attractors $(r\ge 2)$. Finally we prove level-2 large deviations bounds for weak Gibbs measures, which provide a large deviations principle in the special case of Gibbs measures. The main technique we apply is the horseshoe approximation property. △ Less

Submitted 6 July, 2023; v1 submitted 9 November, 2021; originally announced November 2021.

A Petrov-Galerkin method for nonlocal convection-dominated diffusion problems

Authors: Yu Leng, Xiaochuan Tian, Leszek Demkowicz, Hector Gomez, John T. Foster

Abstract: We present a Petrov-Gelerkin (PG) method for a class of nonlocal convection-dominated diffusion problems. There are two main ingredients in our approach. First, we define the norm on the test space as induced by the trial space norm, i.e., the optimal test norm, so that the inf-sup condition can be satisfied uniformly independent of the problem. We show the well-posedness of a class of nonlocal co… ▽ More We present a Petrov-Gelerkin (PG) method for a class of nonlocal convection-dominated diffusion problems. There are two main ingredients in our approach. First, we define the norm on the test space as induced by the trial space norm, i.e., the optimal test norm, so that the inf-sup condition can be satisfied uniformly independent of the problem. We show the well-posedness of a class of nonlocal convection-dominated diffusion problems under the optimal test norm with general assumptions on the nonlocal diffusion and convection kernels. Second, following the framework of Cohen et al.~(2012), we embed the original nonlocal convection-dominated diffusion problem into a larger mixed problem so as to choose an enriched test space as a stabilization of the numerical algorithm. In the numerical experiments, we use an approximate optimal test norm which can be efficiently implemented in 1d, and study its performance against the energy norm on the test space. We conduct convergence studies for the nonlocal problem using uniform $h$- and $p$-refinements, and adaptive $h$-refinements on both smooth manufactured solutions and solutions with sharp gradient in a transition layer. In addition, we confirm that the PG method is asymptotically compatible. △ Less

Submitted 14 July, 2021; originally announced July 2021.

An asymptotically compatible probabilistic collocation method for randomly heterogeneous nonlocal problems

Authors: Yiming Fan, Xiaochuan Tian, Xiu Yang, Xingjie Li, Clayton Webster, Yue Yu

Abstract: In this paper we present an asymptotically compatible meshfree method for solving nonlocal equations with random coefficients, describing diffusion in heterogeneous media. In particular, the random diffusivity coefficient is described by a finite-dimensional random variable or a truncated combination of random variables with the Karhunen-Loève decomposition, then a probabilistic collocation method… ▽ More In this paper we present an asymptotically compatible meshfree method for solving nonlocal equations with random coefficients, describing diffusion in heterogeneous media. In particular, the random diffusivity coefficient is described by a finite-dimensional random variable or a truncated combination of random variables with the Karhunen-Loève decomposition, then a probabilistic collocation method (PCM) with sparse grids is employed to sample the stochastic process. On each sample, the deterministic nonlocal diffusion problem is discretized with an optimization-based meshfree quadrature rule. We present rigorous analysis for the proposed scheme and demonstrate convergence for a number of benchmark problems, showing that it sustains the asymptotic compatibility spatially and achieves an algebraic or sub-exponential convergence rate in the random coefficients space as the number of collocation points grows. Finally, to validate the applicability of this approach we consider a randomly heterogeneous nonlocal problem with a given spatial correlation structure, demonstrating that the proposed PCM approach achieves substantial speed-up compared to conventional Monte Carlo simulations. △ Less

Submitted 15 July, 2021; v1 submitted 3 July, 2021; originally announced July 2021.

MSC Class: 45K05; 35R09; 35R60; 65N12; 35B65

Ergodic Average Of Typical Orbits And Typical Functions

Authors: Xiaobo Hou, Wanshan Lin, Xueting Tian

Abstract: In this article we mainly aim to know what kind of asymptotic behavior of typical orbits can display. For example, we show in any transitive system, the emprical measures of a typical orbit can cover all emprical measures of dense orbits and can intersect some physical-like measures. In particular, if the union set of emprical measures of all dense orbits is not singleton, then the typical orbit w… ▽ More In this article we mainly aim to know what kind of asymptotic behavior of typical orbits can display. For example, we show in any transitive system, the emprical measures of a typical orbit can cover all emprical measures of dense orbits and can intersect some physical-like measures. In particular, if the union set of emprical measures of all dense orbits is not singleton, then the typical orbit will display historic behavior simultaneously for typical continuous functions and the limit set of ergodic average along every continuous function equals to a closed interval composed by the union of limit sets of ergodic average on all dense orbits. Moreover, if the union set of emprical measures of all dense orbits contains all ergodic measures, the above interval equals to the rotation set. These results are not only suitable for systems with specification-like properties or minimal systems, but also suitable for many other systems including all general (not assumed uniformly hyperbolic) nontrivial homoclinic classes and Bowen eyes. Moreover, we introduce a new property called m-g-product property weaker than classical specification property and minimal property and nontrivial examples are constructed. △ Less

Submitted 12 November, 2021; v1 submitted 30 June, 2021; originally announced July 2021.

Comments: Reorganized version

Nonlocal Trace Spaces and Extension Results for Nonlocal Calculus

Authors: Qiang Du, Xiaochuan Tian, Cory Wright, Yue Yu

Abstract: For a given Lipschitz domain $Ω$, it is a classical result that the trace space of $W^{1,p}(Ω)$ is $W^{1-1/p,p}(\partialΩ)$, namely any $W^{1,p}(Ω)$ function has a well-defined $W^{1-1/p,p}(\partialΩ)$ trace on its codimension-1 boundary $\partialΩ$ and any $W^{1-1/p,p}(\partialΩ)$ function on $\partialΩ$ can be extended to a $W^{1,p}(Ω)$ function. Recently, Dyda and Kassmann (2019) characterize t… ▽ More For a given Lipschitz domain $Ω$, it is a classical result that the trace space of $W^{1,p}(Ω)$ is $W^{1-1/p,p}(\partialΩ)$, namely any $W^{1,p}(Ω)$ function has a well-defined $W^{1-1/p,p}(\partialΩ)$ trace on its codimension-1 boundary $\partialΩ$ and any $W^{1-1/p,p}(\partialΩ)$ function on $\partialΩ$ can be extended to a $W^{1,p}(Ω)$ function. Recently, Dyda and Kassmann (2019) characterize the trace space for nonlocal Dirichlet problems involving integrodifferential operators with infinite interaction ranges, where the boundary datum is provided on the whole complement of the given domain $\mathbb{R}^d\backslashΩ$. In this work, we study function spaces for nonlocal Dirichlet problems with a finite range of nonlocal interactions, which naturally serves a bridging role between the classical local PDE problem and the nonlocal problem with infinite interaction ranges. For these nonlocal Dirichlet problems, the boundary conditions are normally imposed on a region with finite thickness volume which lies outside of the domain. We introduce a function space on the volumetric boundary region that serves as a trace space for these nonlocal problems and study the related extension results. Moreover, we discuss the consistency of the new nonlocal trace space with the classical $W^{1-1/p,p}(\partialΩ)$ space as the size of nonlocal interaction tends to zero. In making this connection, we conduct an investigation on the relations between nonlocal interactions on a larger domain and the induced interactions on its subdomain. The various forms of trace, embedding and extension theorems may then be viewed as consequences in different scaling limits. △ Less

Submitted 15 July, 2021; v1 submitted 30 June, 2021; originally announced July 2021.

MSC Class: 46E35; 47G10; 35A23; 35R11