-
Global well-posedness of the free boundary problem for incompressible viscous resistive MHD in critical Besov spaces
Authors:
Wei Zhang,
Jie Fu,
Chengchun Hao,
Siqi Yang
Abstract:
This paper aims to establish the global well-posedness of the free boundary problem for the incompressible viscous resistive magnetohydrodynamic (MHD) equations. Under the framework of Lagrangian coordinates, a unique global solution exists in the half-space provided that the norm of the initial data in the critical homogeneous Besov space $\dot{B}_{p, 1}^{-1+N/p}(\mathbb{R}_{+}^N)$ is sufficientl…
▽ More
This paper aims to establish the global well-posedness of the free boundary problem for the incompressible viscous resistive magnetohydrodynamic (MHD) equations. Under the framework of Lagrangian coordinates, a unique global solution exists in the half-space provided that the norm of the initial data in the critical homogeneous Besov space $\dot{B}_{p, 1}^{-1+N/p}(\mathbb{R}_{+}^N)$ is sufficiently small, where $p \in [N, 2N-1)$. Building upon prior work such as (Danchin and Mucha, J. Funct. Anal. 256 (2009) 881--927) and (Ogawa and Shimizu, J. Differ. Equations 274 (2021) 613--651) in the half-space setting, we establish maximal $L^{1}$-regularity for both the Stokes equations without surface stress and the linearized equations of the magnetic field with zero boundary condition. The existence and uniqueness of solutions to the nonlinear problems are proven using the Banach contraction mapping principle.
△ Less
Submitted 19 August, 2024;
originally announced August 2024.
-
A Derivative-Free Martingale Neural Network Soc-Martnet For The Hamilton-Jacobi-Bellman Equations In Stochastic Optimal Controls
Authors:
Wei Cai,
Shuixin Fang,
Wenzhong Zhang,
Tao Zhou
Abstract:
In this paper, we propose an efficient derivative-free version of a martingale neural network SOC-MartNet proposed in Cai et al. [2] for solving high-dimensional Hamilton-Jacobi-Bellman (HJB) equations and stochastic optimal control problems (SOCPs) with controls on both drift and volatility. The solution of the HJB equation consists of two steps: (1) finding the optimal control from the value fun…
▽ More
In this paper, we propose an efficient derivative-free version of a martingale neural network SOC-MartNet proposed in Cai et al. [2] for solving high-dimensional Hamilton-Jacobi-Bellman (HJB) equations and stochastic optimal control problems (SOCPs) with controls on both drift and volatility. The solution of the HJB equation consists of two steps: (1) finding the optimal control from the value function, and (2) deriving the value function from a linear PDE characterized by the optimal control. The linear PDE is reformulated into a weak form of a new martingale formulation from the original SOC-MartNet where all temporal and spatial derivatives are replaced by an univariate, first-order random finite difference operator approximation, giving the derivative free version of the SOC-MartNet. Then, the optimal feedback control is identified by minimizing the mean of the value function, thereby avoiding the need for pointwise minimization on the Hamiltonian. Finally, the optimal control and value function are approximated by neural networks trained via adversarial learning using the derivative-free formulation. This method eliminates the reliance on automatic differentiation for computing temporal and spatial derivatives, offering significant efficiency in solving high-dimensional HJB equations and SOCPs.
△ Less
Submitted 26 August, 2024; v1 submitted 26 August, 2024;
originally announced August 2024.
-
A Drinfeld type presentation of twisted Yangians of quasi-split type
Authors:
Kang Lu,
Weinan Zhang
Abstract:
We formulate a family of algebras, twisted Yangians (of simply-laced quasi-split type) in Drinfeld type current generators and defining relations. These new algebras admit PBW type bases and are shown to be a deformation of twisted current algebras. For all quasi-split type excluding the even rank case in type AIII, we show that the twisted Yangians can be realized via a degeneration on the Drinfe…
▽ More
We formulate a family of algebras, twisted Yangians (of simply-laced quasi-split type) in Drinfeld type current generators and defining relations. These new algebras admit PBW type bases and are shown to be a deformation of twisted current algebras. For all quasi-split type excluding the even rank case in type AIII, we show that the twisted Yangians can be realized via a degeneration on the Drinfeld type presentation of affine $\imath$quantum groups. For both even and odd rank cases in type AIII, we use the Gauss decomposition method to show that these new algebras are isomorphic to Molev-Ragoucy's reflection algebras defined in the R-matrix presentation.
△ Less
Submitted 13 August, 2024;
originally announced August 2024.
-
Higher $q$-Continued Fractions
Authors:
Amanda Burcroff,
Nicholas Ovenhouse,
Ralf Schiffler,
Sylvester W. Zhang
Abstract:
We introduce a $q$-analog of the higher continued fractions introduced by the last three authors in a previous work (together with Gregg Musiker), which are simultaneously a generalization of the $q$-rational numbers of Morier-Genoud and Ovsienko. They are defined as ratios of generating functions for $P$-partitions on certain posets. We give matrix formulas for computing them, which generalize pr…
▽ More
We introduce a $q$-analog of the higher continued fractions introduced by the last three authors in a previous work (together with Gregg Musiker), which are simultaneously a generalization of the $q$-rational numbers of Morier-Genoud and Ovsienko. They are defined as ratios of generating functions for $P$-partitions on certain posets. We give matrix formulas for computing them, which generalize previous results in the $q=1$ case. We also show that certain properties enjoyed by the $q$-rationals are also satisfied by our higher versions.
△ Less
Submitted 13 August, 2024;
originally announced August 2024.
-
UpLIF: An Updatable Self-Tuning Learned Index Framework
Authors:
Alireza Heidari,
Amirhossein Ahmadi,
Wei Zhang
Abstract:
The emergence of learned indexes has caused a paradigm shift in our perception of indexing by considering indexes as predictive models that estimate keys' positions within a data set, resulting in notable improvements in key search efficiency and index size reduction; however, a significant challenge inherent in learned index modeling is its constrained support for update operations, necessitated…
▽ More
The emergence of learned indexes has caused a paradigm shift in our perception of indexing by considering indexes as predictive models that estimate keys' positions within a data set, resulting in notable improvements in key search efficiency and index size reduction; however, a significant challenge inherent in learned index modeling is its constrained support for update operations, necessitated by the requirement for a fixed distribution of records. Previous studies have proposed various approaches to address this issue with the drawback of high overhead due to multiple model retraining. In this paper, we present UpLIF, an adaptive self-tuning learned index that adjusts the model to accommodate incoming updates, predicts the distribution of updates for performance improvement, and optimizes its index structure using reinforcement learning. We also introduce the concept of balanced model adjustment, which determines the model's inherent properties (i.e. bias and variance), enabling the integration of these factors into the existing index model without the need for retraining with new data. Our comprehensive experiments show that the system surpasses state-of-the-art indexing solutions (both traditional and ML-based), achieving an increase in throughput of up to 3.12 times with 1000 times less memory usage.
△ Less
Submitted 7 August, 2024;
originally announced August 2024.
-
More results on the spectral radius of graphs with no odd wheels
Authors:
Wenqian Zhang
Abstract:
For a graph $G$, the spectral radius $λ_{1}(G)$ of $G$ is the largest eigenvalue of its adjacency matrix. An odd wheel $W_{2k+1}$ with $k\geq2$ is a graph obtained from a cycle of order $2k$ by adding a new vertex connecting to all the vertices of the cycle. Let ${\rm SPEX}(n,W_{2k+1})$ be the set of $W_{2k+1}$-free graphs of order $n$ with the maximum spectral radius. Very recently, Cioabă, Desai…
▽ More
For a graph $G$, the spectral radius $λ_{1}(G)$ of $G$ is the largest eigenvalue of its adjacency matrix. An odd wheel $W_{2k+1}$ with $k\geq2$ is a graph obtained from a cycle of order $2k$ by adding a new vertex connecting to all the vertices of the cycle. Let ${\rm SPEX}(n,W_{2k+1})$ be the set of $W_{2k+1}$-free graphs of order $n$ with the maximum spectral radius. Very recently, Cioabă, Desai and Tait \cite{CDT2} characterized the graphs in ${\rm SPEX}(n,W_{2k+1})$ for sufficiently large $n$, where $k\geq2$ and $k\neq4,5$. And they left the case $k=4,5$ as a problem. In this paper, we settle this problem. Moreover, we completely characterize the graphs in ${\rm SPEX}(n,W_{2k+1})$ when $k\geq4$ is even and $n\equiv2~(\mod4)$ is sufficiently large. Consequently, the graphs in ${\rm SPEX}(n,W_{2k+1})$ are characterized completely for any $k\geq2$ and sufficiently large $n$.
△ Less
Submitted 7 August, 2024;
originally announced August 2024.
-
Gaussian Approximations for the $k$th coordinate of sums of random vectors
Authors:
Yixi Ding,
Qizhai Li,
Yuke Shi,
Wei Zhang
Abstract:
We consider the problem of Gaussian approximation for the $κ$th coordinate of a sum of high-dimensional random vectors. Such a problem has been studied previously for $κ=1$ (i.e., maxima). However, in many applications, a general $κ\geq1$ is of great interest, which is addressed in this paper. We make four contributions: 1) we first show that the distribution of the $κ$th coordinate of a sum of ra…
▽ More
We consider the problem of Gaussian approximation for the $κ$th coordinate of a sum of high-dimensional random vectors. Such a problem has been studied previously for $κ=1$ (i.e., maxima). However, in many applications, a general $κ\geq1$ is of great interest, which is addressed in this paper. We make four contributions: 1) we first show that the distribution of the $κ$th coordinate of a sum of random vectors, $\boldsymbol{X}= (X_{1},\cdots,X_{p})^{\sf T}= n^{-1/2}\sum_{i=1}^n \boldsymbol{x}_{i}$, can be approximated by that of Gaussian random vectors and derive their Kolmogorov's distributional difference bound; 2) we provide the theoretical justification for estimating the distribution of the $κ$th coordinate of a sum of random vectors using a Gaussian multiplier procedure, which multiplies the original vectors with i.i.d. standard Gaussian random variables; 3) we extend the Gaussian approximation result and Gaussian multiplier bootstrap procedure to a more general case where $κ$ diverges; 4) we further consider the Gaussian approximation for a square sum of the first $d$ largest coordinates of $\boldsymbol{X}$. All these results allow the dimension $p$ of random vectors to be as large as or much larger than the sample size $n$.
△ Less
Submitted 6 August, 2024;
originally announced August 2024.
-
Implementable Semismooth* Newton Methods for Generalized Equations are G-Semismooth Newton Methods
Authors:
Liang Chen,
Defeng Sun,
Wangyongquan Zhang
Abstract:
Semismooth* Newton methods have been proposed in recent years targeting multi-valued inclusion problems and have been successfully implemented to deal with several concrete generalized equations. In this paper, we show that these executable implementations are exactly the applications of G-semismooth Newton methods for solving nonsmooth equations localized from these generalized equations. This ne…
▽ More
Semismooth* Newton methods have been proposed in recent years targeting multi-valued inclusion problems and have been successfully implemented to deal with several concrete generalized equations. In this paper, we show that these executable implementations are exactly the applications of G-semismooth Newton methods for solving nonsmooth equations localized from these generalized equations. This new understanding expands the breadth of G-semismooth Newton methods in theory, and more importantly, facilitates the design and implementation of practical Newton-type algorithms for solving generalized equations.
△ Less
Submitted 19 July, 2024;
originally announced July 2024.
-
PBW bases for $\imath$quantum groups
Authors:
Ming Lu,
Ruiqi Yang,
Weinan Zhang
Abstract:
We establish PBW type bases for $\imath$quantum groups of arbitrary finite type, using the relative braid group symmetries. Explicit formulas for root vectors are provided for $\imath$quantum groups of each rank 1 type. We show that our PBW type bases give rise to integral bases for the modified $\imath$quantum groups. The leading terms of our bases can be identified with the usual PBW bases in th…
▽ More
We establish PBW type bases for $\imath$quantum groups of arbitrary finite type, using the relative braid group symmetries. Explicit formulas for root vectors are provided for $\imath$quantum groups of each rank 1 type. We show that our PBW type bases give rise to integral bases for the modified $\imath$quantum groups. The leading terms of our bases can be identified with the usual PBW bases in the theory of quantum groups.
△ Less
Submitted 17 July, 2024;
originally announced July 2024.
-
Normalized solution for $p$-Laplacian equation in exterior domain
Authors:
Weiqiang Zhang,
Yanyun Wen
Abstract:
We are devoted to the study of the following nonlinear $p$-Laplacian Schrödinger equation with $L^{p}$-norm constraint \begin{align*} \begin{cases} &-Δ_{p} u=λ|u|^{p-2}u +|u|^{r-2}u\quad\mbox{in}\quadΩ,\\ &u=0\quad\mbox{on}\quad \partialΩ,\\ &\int_Ω|u|^{p}dx=a, \end{cases} \end{align*} where $Δ_{p}u=\text{div} (|\nabla u|^{p-2}\nabla u)$, $Ω\subset\mathbb{R}^{N}$ is an exterior domain with smooth…
▽ More
We are devoted to the study of the following nonlinear $p$-Laplacian Schrödinger equation with $L^{p}$-norm constraint \begin{align*} \begin{cases} &-Δ_{p} u=λ|u|^{p-2}u +|u|^{r-2}u\quad\mbox{in}\quadΩ,\\ &u=0\quad\mbox{on}\quad \partialΩ,\\ &\int_Ω|u|^{p}dx=a, \end{cases} \end{align*} where $Δ_{p}u=\text{div} (|\nabla u|^{p-2}\nabla u)$, $Ω\subset\mathbb{R}^{N}$ is an exterior domain with smooth boundary $\partialΩ\neq\emptyset$ satisfying that $\R^{N}\setminusΩ$ is bounded, $N\geq3$, $2\leq p<N$, $p<r<p+\frac{p^{2}}{N}$, $a>0$ and $λ\in\R$ is an unknown Lagrange multiplier. First, by using the splitting techniques and the Gagliardo-Nirenberg inequality, the compactness of Palais-Smale sequence of the above problem at higher energy level is established. Then, exploiting barycentric function methods, Brouwer degree and minimax principle, we obtain a solution $(u,\la)$ with $u>0$ in $\R^{N}$ and $\la<0$ when $\R^{N}\setminusΩ$ is contained in a small ball. Moreover, we give a similar result if we remove the restriction on $Ω$ and assume $a>0$ small enough. Last, with the symmetric assumption on $Ω$, we use genus theory to consider infinite many solutions.
△ Less
Submitted 16 July, 2024;
originally announced July 2024.
-
Noncommutative Chern form on étale groupoid is closed
Authors:
Wen Zhang
Abstract:
We use bisection to provide an algebraic proof that the Chern form on the convolution algebra of an étale groupoid is closed.
We use bisection to provide an algebraic proof that the Chern form on the convolution algebra of an étale groupoid is closed.
△ Less
Submitted 21 July, 2024; v1 submitted 12 July, 2024;
originally announced July 2024.
-
The multiple birth properties of multi-type Markov branching processes
Authors:
Junping Li,
Wanting Zhang
Abstract:
The main purpose of this paper is to consider the multiple birth properties for multi-type Markov branching processes. We first construct a new multi-dimensional Markov process based on the multi-type Markov branching process, which can reveal the multiple birth characteristics. Then the joint probability distribution of multiple birth of multi-type Markov branching process until any time $t$ is o…
▽ More
The main purpose of this paper is to consider the multiple birth properties for multi-type Markov branching processes. We first construct a new multi-dimensional Markov process based on the multi-type Markov branching process, which can reveal the multiple birth characteristics. Then the joint probability distribution of multiple birth of multi-type Markov branching process until any time $t$ is obtained by using the new process. Furthermore, the probability distribution of multiple birth until the extinction of the process is also given.
△ Less
Submitted 6 July, 2024;
originally announced July 2024.
-
Extended alternating structure-adapted proximal gradient algorithm for nonconvex nonsmooth problems
Authors:
Ying Gao,
Chunfeng Cui,
Wenxing Zhang,
Deren Han
Abstract:
Alternating structure-adapted proximal (ASAP) gradient algorithm (M. Nikolova and P. Tan, SIAM J Optim, 29:2053-2078, 2019) has drawn much attention due to its efficiency in solving nonconvex nonsmooth optimization problems. However, the multiblock nonseparable structure confines the performance of ASAP to far-reaching practical problems, e.g., coupled tensor decomposition. In this paper, we propo…
▽ More
Alternating structure-adapted proximal (ASAP) gradient algorithm (M. Nikolova and P. Tan, SIAM J Optim, 29:2053-2078, 2019) has drawn much attention due to its efficiency in solving nonconvex nonsmooth optimization problems. However, the multiblock nonseparable structure confines the performance of ASAP to far-reaching practical problems, e.g., coupled tensor decomposition. In this paper, we propose an extended ASAP (eASAP) algorithm for nonconvex nonsmooth optimization whose objective is the sum of two nonseperable functions and a coupling one. By exploiting the blockwise restricted prox-regularity, eASAP is capable of minimizing the objective whose coupling function is multiblock nonseparable. Moreover, we analyze the global convergence of eASAP by virtue of the Aubin property on partial subdifferential mapping and the Kurdyka-Łojasiewicz property on the objective. Furthermore, the sublinear convergence rate of eASAP is built upon the proximal point algorithmic framework under some mild conditions. Numerical simulations on multimodal data fusion demonstrate the compelling performance of the proposed method.
△ Less
Submitted 24 June, 2024;
originally announced June 2024.
-
Toughness and spectral radius in graphs
Authors:
Sufang Wang,
Wei Zhang
Abstract:
Let $t$ be a positive integer, and let $G$ be a connected graph of order $n$ with $n\geq t+2$. A graph $G$ is said to be $\frac{1}{t}$-tough if $|S|\geq\frac{1}{t}c(G-S)$ for every subset $S$ of $V(G)$ with $c(G-S)\geq2$, where $c(G-S)$ is the number of connected components in $G-S$. The adjacency matrix of $G$ is denoted by $A(G)$. Let $λ_1(G)\geqλ_2(G)\geq\dots\geqλ_n(G)$ be the eigenvalues of…
▽ More
Let $t$ be a positive integer, and let $G$ be a connected graph of order $n$ with $n\geq t+2$. A graph $G$ is said to be $\frac{1}{t}$-tough if $|S|\geq\frac{1}{t}c(G-S)$ for every subset $S$ of $V(G)$ with $c(G-S)\geq2$, where $c(G-S)$ is the number of connected components in $G-S$. The adjacency matrix of $G$ is denoted by $A(G)$. Let $λ_1(G)\geqλ_2(G)\geq\dots\geqλ_n(G)$ be the eigenvalues of $A(G)$. In particular, the eigenvalue $λ_1(G)$ is called the spectral radius of $G$. In this paper, we prove that $G$ is a $\frac{1}{t}$-tough graph unless $G=K_1\vee(K_{n-t-1}\cup tK_1)$ if $λ_1(G)\geqη(t,n)$, where $η(t,n)$ is the largest root of $x^{3}-(n-t-2)x^{2}-(n-1)x+t(n-t-2)=0$.
△ Less
Submitted 12 June, 2024;
originally announced June 2024.
-
Walks, infinite series and spectral radius of graphs
Authors:
Wenqian Zhang
Abstract:
For a graph G, the spectral radius \r{ho}(G) of G is the largest eigenvalue of its adjacency matrix. In this paper, we seek the relationship between \r{ho}(G) and the walks of the subgraphs of G. Especially, if G contains a complete multi-partite graph as a spanning subgraph, we give a formula for \r{ho}(G) by using an infinite series on walks of the subgraphs of G. These results are useful for th…
▽ More
For a graph G, the spectral radius \r{ho}(G) of G is the largest eigenvalue of its adjacency matrix. In this paper, we seek the relationship between \r{ho}(G) and the walks of the subgraphs of G. Especially, if G contains a complete multi-partite graph as a spanning subgraph, we give a formula for \r{ho}(G) by using an infinite series on walks of the subgraphs of G. These results are useful for the current popular spectral extremal problem.
△ Less
Submitted 11 June, 2024;
originally announced June 2024.
-
Spanning trees and signless Laplacian spectral radius in graphs
Authors:
Sufang Wang,
Wei Zhang
Abstract:
Let $G$ be a connected graph and let $k$ be a positive integer. Let $T$ be a spanning tree of $G$. The leaf degree of a vertex $v\in V(T)$ is defined as the number of leaves adjacent to $v$ in $T$. The leaf degree of $T$ is the maximum leaf degree among all the vertices of $T$. Let $A(G)$ be the adjacency matrix of $G$ and $D(G)$ be the diagonal degree matrix of $G$. Let $Q(G)=D(G)+A(G)$ be the si…
▽ More
Let $G$ be a connected graph and let $k$ be a positive integer. Let $T$ be a spanning tree of $G$. The leaf degree of a vertex $v\in V(T)$ is defined as the number of leaves adjacent to $v$ in $T$. The leaf degree of $T$ is the maximum leaf degree among all the vertices of $T$. Let $A(G)$ be the adjacency matrix of $G$ and $D(G)$ be the diagonal degree matrix of $G$. Let $Q(G)=D(G)+A(G)$ be the signless Laplacian matrix of $G$. The largest eigenvalue of $Q(G)$, denoted by $q(G)$, is called the signless Laplacian spectral radius of $G$. In this paper, we investigate the connection between the spanning tree and the signless Laplacian spectral radius of $G$, and put forward a sufficient condition based upon the signless Laplacian spectral radius to guarantee that a graph $G$ contains a spanning tree with leaf degree at most $k$. Finally, we construct some extremal graphs to claim all the bounds obtained in this paper are sharp.
△ Less
Submitted 11 June, 2024;
originally announced June 2024.
-
An efficient branch-and-cut approach for large-scale competitive facility location problems with limited choice rule
Authors:
Wei-Kun Chen,
Wei-Yang Zhang,
Yan-Ru Wang,
Shahin Gelareh,
Yu-Hong Dai
Abstract:
In the paper, we consider the competitive facility location problem with limited choice rule (CFLPLCR), which attempts to open a subset of facilities to maximize the net profit of a newcomer company, requiring customers to patronize only a limited number of opening facilities and an outside option. We propose an efficient branch-and-cut (B&C) approach for the CFLPLCR based on newly proposed mixed…
▽ More
In the paper, we consider the competitive facility location problem with limited choice rule (CFLPLCR), which attempts to open a subset of facilities to maximize the net profit of a newcomer company, requiring customers to patronize only a limited number of opening facilities and an outside option. We propose an efficient branch-and-cut (B&C) approach for the CFLPLCR based on newly proposed mixed integer linear programming (MILP) formulations. Specifically, by establishing the submodularity of the probability function, we develop an MILP formulation for the CFLPLCR using the submodular inequalities. For the special case where each customer patronizes at most one open facility and the outside option, we show that the submodular inequalities can characterize the convex hull of the considered set and provide a compact MILP formulation. Moreover, for the general case, we strengthen the submodular inequalities by sequential lifting, resulting in a class of facet-defining inequalities. The proposed lifted submodular inequalities are shown to be stronger than the classic submodular inequalities, enabling to obtain another MILP formulation with a tighter linear programming (LP) relaxation. By extensive numerical experiments, we show that the proposed B&C approach outperforms the state-of-the-art generalized Benders decomposition approach by at least one order of magnitude. Furthermore, it enables to solve CFLPLCR instances with 10000 customers and 2000 facilities.
△ Less
Submitted 9 June, 2024;
originally announced June 2024.
-
Affine $\imath$quantum groups and twisted Yangians in Drinfeld presentations
Authors:
Kang Lu,
Weiqiang Wang,
Weinan Zhang
Abstract:
We formulate a family of algebras, twisted Yangians (of split type) in current generators and relations, via a degeneration of the Drinfeld presentation of affine $\imath$quantum groups (associated with split Satake diagrams). These new algebras admit PBW type bases and are shown to be a deformation of twisted current algebras; presentations for twisted current algebras are also provided. For type…
▽ More
We formulate a family of algebras, twisted Yangians (of split type) in current generators and relations, via a degeneration of the Drinfeld presentation of affine $\imath$quantum groups (associated with split Satake diagrams). These new algebras admit PBW type bases and are shown to be a deformation of twisted current algebras; presentations for twisted current algebras are also provided. For type AI, it matches with the Drinfeld presentation of twisted Yangian obtained via Gauss decomposition. We conjecture that our split twisted Yangians are isomorphic to the corresponding ones in RTT presentation.
△ Less
Submitted 7 June, 2024;
originally announced June 2024.
-
A novel model reduction method to solve inverse problems of parabolic type
Authors:
Wenlong Zhang,
Zhiwen Zhang
Abstract:
In this paper, we propose novel proper orthogonal decomposition (POD)--based model reduction methods that effectively address the issue of inverse crime in solving parabolic inverse problems. Both the inverse initial value problems and inverse source problems are studied. By leveraging the inherent low-dimensional structures present in the data, our approach enables a reduction in the forward mode…
▽ More
In this paper, we propose novel proper orthogonal decomposition (POD)--based model reduction methods that effectively address the issue of inverse crime in solving parabolic inverse problems. Both the inverse initial value problems and inverse source problems are studied. By leveraging the inherent low-dimensional structures present in the data, our approach enables a reduction in the forward model complexity without compromising the accuracy of the inverse problem solution. Besides, we prove the convergence analysis of the proposed methods for solving parabolic inverse problems. Through extensive experimentation and comparative analysis, we demonstrate the effectiveness of our method in overcoming inverse crime and achieving improved inverse problem solutions. The proposed POD model reduction method offers a promising direction for improving the reliability and applicability of inverse problem-solving techniques in various domains.
△ Less
Submitted 4 June, 2024;
originally announced June 2024.
-
On a multi-dimensional transport equation with nonlocal velocity and fractional dissipation
Authors:
Wanwan Zhang
Abstract:
This paper aims to investigate a multi-dimensional transport equation with nonlocal velocity and fractional dissipation.The balance between the nonlinearity and dissipation gives rise to three different cases, namely the subcritical, critical and supercritical ranges. We study those three cases and obtain a set of results containing local well-posedness, global smoothness, eventual regularity and…
▽ More
This paper aims to investigate a multi-dimensional transport equation with nonlocal velocity and fractional dissipation.The balance between the nonlinearity and dissipation gives rise to three different cases, namely the subcritical, critical and supercritical ranges. We study those three cases and obtain a set of results containing local well-posedness, global smoothness, eventual regularity and finite-time blowup of smooth solutions.
△ Less
Submitted 25 May, 2024;
originally announced May 2024.
-
Restarted Primal-Dual Hybrid Conjugate Gradient Method for Large-Scale Quadratic Programming
Authors:
Yicheng Huang,
Wanyu Zhang,
Hongpei Li,
Weihan Xue,
Dongdong Ge,
Huikang Liu,
Yinyu Ye
Abstract:
Convex quadratic programming (QP) is an essential class of optimization problems with broad applications across various fields. Traditional QP solvers, typically based on simplex or barrier methods, face significant scalability challenges. In response to these limitations, recent research has shifted towards matrix-free first-order methods to enhance scalability in QP. Among these, the restarted a…
▽ More
Convex quadratic programming (QP) is an essential class of optimization problems with broad applications across various fields. Traditional QP solvers, typically based on simplex or barrier methods, face significant scalability challenges. In response to these limitations, recent research has shifted towards matrix-free first-order methods to enhance scalability in QP. Among these, the restarted accelerated primal-dual hybrid gradient (rAPDHG) method, proposed by H.Lu(2023), has gained notable attention due to its linear convergence rate to an optimal solution and its straightforward implementation on Graphics Processing Units (GPUs). Building on this framework, this paper introduces a restarted primal-dual hybrid conjugate gradient (PDHCG) method, which incorporates conjugate gradient (CG) techniques to address the primal subproblems inexactly. We demonstrate that PDHCG maintains a linear convergence rate with an improved convergence constant and is also straightforward to implement on GPUs. Extensive numerical experiments affirm that, compared to rAPDHG, our method could significantly reduce the number of iterations required to achieve the desired accuracy and offer a substantial performance improvement in large-scale problems. These findings highlight the significant potential of our proposed PDHCG method to boost both the efficiency and scalability of solving complex QP challenges.
△ Less
Submitted 25 May, 2024;
originally announced May 2024.
-
On finding optimal collective variables for complex systems by minimizing the deviation between effective and full dynamics
Authors:
Wei Zhang,
Christof Schütte
Abstract:
This paper is concerned with collective variables, or reaction coordinates, that map a discrete-in-time Markov process $X_n$ in $\mathbb{R}^d$ to a (much) smaller dimension $k \ll d$. We define the effective dynamics under a given collective variable map $ξ$ as the best Markovian representation of $X_n$ under $ξ$. The novelty of the paper is that it gives strict criteria for selecting optimal coll…
▽ More
This paper is concerned with collective variables, or reaction coordinates, that map a discrete-in-time Markov process $X_n$ in $\mathbb{R}^d$ to a (much) smaller dimension $k \ll d$. We define the effective dynamics under a given collective variable map $ξ$ as the best Markovian representation of $X_n$ under $ξ$. The novelty of the paper is that it gives strict criteria for selecting optimal collective variables via the properties of the effective dynamics. In particular, we show that the transition density of the effective dynamics of the optimal collective variable solves a relative entropy minimization problem from certain family of densities to the transition density of $X_n$. We also show that many transfer operator-based data-driven numerical approaches essentially learn quantities of the effective dynamics. Furthermore, we obtain various error estimates for the effective dynamics in approximating dominant timescales / eigenvalues and transition rates of the original process $X_n$ and how optimal collective variables minimize these errors. Our results contribute to the development of theoretical tools for the understanding of complex dynamical systems, e.g. molecular kinetics, on large timescales. These results shed light on the relations among existing data-driven numerical approaches for identifying good collective variables, and they also motivate the development of new methods.
△ Less
Submitted 3 May, 2024;
originally announced May 2024.
-
A Double Maximization Approach for Optimizing the LM Rate of Mismatched Decoding
Authors:
Lingyi Chen,
Shitong Wu,
Xinwei Li,
Huihui Wu,
Hao Wu,
Wenyi Zhang
Abstract:
An approach is established for maximizing the Lower bound on the Mismatch capacity (hereafter abbreviated as LM rate), a key performance bound in mismatched decoding, by optimizing the channel input probability distribution. Under a fixed channel input probability distribution, the computation of the corresponding LM rate is a convex optimization problem. When optimizing the channel input probabil…
▽ More
An approach is established for maximizing the Lower bound on the Mismatch capacity (hereafter abbreviated as LM rate), a key performance bound in mismatched decoding, by optimizing the channel input probability distribution. Under a fixed channel input probability distribution, the computation of the corresponding LM rate is a convex optimization problem. When optimizing the channel input probability distribution, however, the corresponding optimization problem adopts a max-min formulation, which is generally non-convex and is intractable with standard approaches. To solve this problem, a novel dual form of the LM rate is proposed, thereby transforming the max-min formulation into an equivalent double maximization formulation. This new formulation leads to a maximization problem setup wherein each individual optimization direction is convex. Consequently, an alternating maximization algorithm is established to solve the resultant maximization problem setup. Each step of the algorithm only involves a closed-form iteration, which is efficiently implemented with standard optimization procedures. Numerical experiments show the proposed approach for optimizing the LM rate leads to noticeable rate gains.
△ Less
Submitted 1 May, 2024;
originally announced May 2024.
-
Llarull's theorem on odd dimensional manifolds: the noncompact case
Authors:
Yihan Li,
Guangxiang Su,
Xiangsheng Wang,
Weiping Zhang
Abstract:
Let $(M,g^{TM})$ be an odd dimensional ($\dim M\geq 3$) connected oriented noncompact complete spin Riemannian manifold. Let $k^{TM}$ be the associated scalar curvature. Let $f:M\to S^{\dim M}(1)$ be a smooth area decreasing map which is locally constant near infinity and of nonzero degree. Suppose $k^{TM}\geq ({\dim M})({\dim M}-1)$ on the support of ${\rm d}f$, we show that $\inf(k^{TM})<0$. Thi…
▽ More
Let $(M,g^{TM})$ be an odd dimensional ($\dim M\geq 3$) connected oriented noncompact complete spin Riemannian manifold. Let $k^{TM}$ be the associated scalar curvature. Let $f:M\to S^{\dim M}(1)$ be a smooth area decreasing map which is locally constant near infinity and of nonzero degree. Suppose $k^{TM}\geq ({\dim M})({\dim M}-1)$ on the support of ${\rm d}f$, we show that $\inf(k^{TM})<0$. This answers a question of Gromov.
△ Less
Submitted 28 April, 2024;
originally announced April 2024.
-
Stability and Bounded Real Lemmas of Discrete-Time MJLSs with the Markov Chain on a Borel Space
Authors:
Chunjie Xiao,
Ting Hou,
Weihai Zhang
Abstract:
In this paper, exponential stability of discrete-time Markov jump linear systems (MJLSs) with the Markov chain on a Borel space $(Θ, \mathcal{B}(Θ))$ is studied, and bounded real lemmas (BRLs) are given. The work generalizes the results from the previous literature that considered only the Markov chain taking values in a countable set to the scenario of an uncountable set and provides unified appr…
▽ More
In this paper, exponential stability of discrete-time Markov jump linear systems (MJLSs) with the Markov chain on a Borel space $(Θ, \mathcal{B}(Θ))$ is studied, and bounded real lemmas (BRLs) are given. The work generalizes the results from the previous literature that considered only the Markov chain taking values in a countable set to the scenario of an uncountable set and provides unified approaches for describing exponential stability and $H_{\infty}$ performance of MJLSs. This paper covers two kinds of exponential stabilities: one is exponential mean-square stability with conditioning (EMSSy-C), and the other is exponential mean-square stability (EMSSy). First, based on the infinite-dimensional operator theory, the equivalent conditions for determining these two kinds of stabilities are shown respectively by the exponentially stable evolutions generated by the corresponding bounded linear operators on different Banach spaces, which turn out to present the spectral criteria of EMSSy-C and EMSSy. Furthermore, the relationship between these two kinds of stabilities is discussed. Moreover, some easier-to-check criteria are established for EMSSy-C of MJLSs in terms of the existence of uniformly positive definite solutions of Lyapunov-type equations or inequalities. In addition, BRLs are given separately in terms of the existence of solutions of the $Θ$-coupled difference Riccati equation for the finite horizon case and algebraic Riccati equation for the infinite horizon case, which facilitates the $H_{\infty}$ analysis of MJLSs with the Markov chain on a Borel space.
△ Less
Submitted 24 April, 2024;
originally announced April 2024.
-
Local cohomology and Segre products
Authors:
Jiamin Li,
Wenliang Zhang
Abstract:
We prove a Künneth formula for local cohomology of a Segré product of graded modules supported in a Segré product of ideals. In order to apply our formula to the study of cohomological dimension, we also investigate asymptotic behaviors of Eulerian graded $\scr{D}$-modules.
We prove a Künneth formula for local cohomology of a Segré product of graded modules supported in a Segré product of ideals. In order to apply our formula to the study of cohomological dimension, we also investigate asymptotic behaviors of Eulerian graded $\scr{D}$-modules.
△ Less
Submitted 22 May, 2024; v1 submitted 4 April, 2024;
originally announced April 2024.
-
Symmetry of Convex Solutions to Fully Nonlinear Elliptic Systems: Bounded Domains
Authors:
Weijun Zhang,
Zhitao Zhang
Abstract:
In this paper, we are concerned with the monotonic and symmetric properties of convex solutions to fully nonlinear elliptic systems. We mainly discuss Monge-Ampère type systems for instance, considering
\begin{equation*}
\det(D^2u^i)=f^i(x,{\bf u},\nabla u^i), \ 1\leq i\leq m,
\end{equation*} over bounded domains of various cases, including the bounded smooth simply connected domains and bou…
▽ More
In this paper, we are concerned with the monotonic and symmetric properties of convex solutions to fully nonlinear elliptic systems. We mainly discuss Monge-Ampère type systems for instance, considering
\begin{equation*}
\det(D^2u^i)=f^i(x,{\bf u},\nabla u^i), \ 1\leq i\leq m,
\end{equation*} over bounded domains of various cases, including the bounded smooth simply connected domains and bounded tube shape domains in $\mathbb{R}^n$. We obtain monotonic and symmetric properties of the solutions to the problem with respect to the geometry of domains and the monotonic and symmetric properties of right-hand side terms. The proof is based on carefully using the moving plane method together with various maximum principles and Hopf's lemmas. The existence and uniqueness to an interesting example of such system is also discussed as an application of our results.
△ Less
Submitted 1 April, 2024;
originally announced April 2024.
-
Quasi-canonical AFL and Arithmetic Transfer conjectures at parahoric levels
Authors:
Chao Li,
Michael Rapoport,
Wei Zhang
Abstract:
In the first part of the paper, we formulate several arithmetic transfer conjectures, which are variants of the arithmetic fundamental lemma conjecture in the presence of ramification. The ramification comes from the choice of non-hyperspecial parahoric level structure. We prove a graph version of these arithmetic transfer conjectures, by relating it to the quasi-canonical arithmetic fundamental l…
▽ More
In the first part of the paper, we formulate several arithmetic transfer conjectures, which are variants of the arithmetic fundamental lemma conjecture in the presence of ramification. The ramification comes from the choice of non-hyperspecial parahoric level structure. We prove a graph version of these arithmetic transfer conjectures, by relating it to the quasi-canonical arithmetic fundamental lemma, which we also establish. We relate some of the arithmetic transfer conjectures to the arithmetic fundamental lemma conjecture for the whole Hecke algebra in our recent paper arXiv:2305.14465. As a consequence, we prove these conjectures in some simple cases. In the second part of the paper, we elucidate the structure of an integral model of a certain member of the almost selfdual Rapoport-Zink tower, thereby proving conjectures of Kudla and the second author. This result allows us verify the hypotheses of the graph version of the arithmetic transfer conjectures in a particular case.
△ Less
Submitted 2 April, 2024;
originally announced April 2024.
-
Bifurcation on Fully Nonlinear Elliptic Equations and Systems
Authors:
Jing Gao,
Weijun Zhang,
Zhitao Zhang
Abstract:
In this paper, we study the following fully nonlinear elliptic equations \begin{equation*} \left\{\begin{array}{rl} \left(S_{k}(D^{2}u)\right)^{\frac1k}=λf(-u) & in\quadΩ\\ u=0 & on\quad \partialΩ\\ \end{array} \right. \end{equation*} and coupled systems \begin{equation*} \left\{\begin{array}{rl} (S_{k}(D^{2}u))^\frac1k=λg(-u,-v) & in\quadΩ\\ (S_{k}(D^{2}v))^\frac1k=λh(-u,-v) & in\quadΩ\\ u=v=0 &…
▽ More
In this paper, we study the following fully nonlinear elliptic equations \begin{equation*} \left\{\begin{array}{rl} \left(S_{k}(D^{2}u)\right)^{\frac1k}=λf(-u) & in\quadΩ\\ u=0 & on\quad \partialΩ\\ \end{array} \right. \end{equation*} and coupled systems \begin{equation*} \left\{\begin{array}{rl} (S_{k}(D^{2}u))^\frac1k=λg(-u,-v) & in\quadΩ\\ (S_{k}(D^{2}v))^\frac1k=λh(-u,-v) & in\quadΩ\\ u=v=0 & on\quad \partialΩ\\ \end{array} \right. \end{equation*} dominated by $k$-Hessian operators, where $Ω$ is a $(k$-$1)$-convex bounded domain in $\mathbb{R}^{N}$, $λ$ is a non-negative parameter, $f:\left[0,+\infty\right)\rightarrow\left[0,+\infty\right)$ is a continuous function with zeros only at $0$ and $g,h:\left[0,+\infty\right)\times \left[0,+\infty\right)\rightarrow \left[0,+\infty\right)$ are continuous functions with zeros only at $(\cdot,0)$ and $(0,\cdot)$. We determine the interval of $λ$ about the existence, non-existence, uniqueness and multiplicity of $k$-convex solutions to the above problems according to various cases of $f,g,h$, which is a complete supplement to the known results in previous literature. In particular, the above results are also new for Laplacian and Monge-Ampère operators. We mainly use bifurcation theory, a-priori estimates, various maximum principles and technical strategies in the proof.
△ Less
Submitted 1 April, 2024;
originally announced April 2024.
-
Symmetry of Convex Solutions to Fully Nonlinear Elliptic Systems: Unbounded Domains
Authors:
Weijun Zhang,
Zhitao Zhang
Abstract:
In this paper, we are concerned with the monotonic and symmetric properties of convex solutions Monge-Ampère systems for instance, considering
\begin{equation*}
\det(D^2u^i)=f^i(x,{\bf u},\nabla u^i), \ 1\leq i\leq m,
\end{equation*} over unbounded domains of various cases, including the whole spaces $\mathbb{R}^n$, the half spaces $\mathbb{R}^n_+$ and the unbounded tube shape domains in…
▽ More
In this paper, we are concerned with the monotonic and symmetric properties of convex solutions Monge-Ampère systems for instance, considering
\begin{equation*}
\det(D^2u^i)=f^i(x,{\bf u},\nabla u^i), \ 1\leq i\leq m,
\end{equation*} over unbounded domains of various cases, including the whole spaces $\mathbb{R}^n$, the half spaces $\mathbb{R}^n_+$ and the unbounded tube shape domains in $\mathbb{R}^n$. We obtain monotonic and symmetric properties of the solutions to the problem with respect to the geometry of domains and the monotonic and symmetric properties of right-hand side terms. The proof is based on carefully using the moving plane method together with various maximum principles and Hopf's lemmas.
△ Less
Submitted 1 April, 2024;
originally announced April 2024.
-
Variational method for fractional Hamiltonian system in bounded domain
Authors:
Weimin Zhang
Abstract:
Here we consider the following fractional Hamiltonian system \begin{equation*} \begin{cases} \begin{aligned} (-Δ)^{s} u&=H_v(u,v) \;\;&&\text{in}~Ω,\\ (-Δ)^{s} v&=H_u(u,v) &&\text{in}~Ω,\\ u &= v = 0 &&\text{in} ~ \mathbb{R}^N\setminusΩ, \end{aligned} \end{cases} \end{equation*} where $s\in (0,1)$, $N>2s$, $H \in C^1(\mathbb{R}^2, \mathbb{R})$ and $Ω\subset \mathbb{R}^N$ is a smooth bounded domain…
▽ More
Here we consider the following fractional Hamiltonian system \begin{equation*} \begin{cases} \begin{aligned} (-Δ)^{s} u&=H_v(u,v) \;\;&&\text{in}~Ω,\\ (-Δ)^{s} v&=H_u(u,v) &&\text{in}~Ω,\\ u &= v = 0 &&\text{in} ~ \mathbb{R}^N\setminusΩ, \end{aligned} \end{cases} \end{equation*} where $s\in (0,1)$, $N>2s$, $H \in C^1(\mathbb{R}^2, \mathbb{R})$ and $Ω\subset \mathbb{R}^N$ is a smooth bounded domain. %As the problem remains unchanged if $H(u, v)$ is replaced by $H(u, v)-H(0, 0)$, hence we always assume $H(0,0)=0$.
To apply the variational method for this problem, the key question is to find a suitable functional setting. Instead of usual fractional Sobolev spaces, we use the solutions space of $(-Δ)^{s}u=f\in L^r(Ω)$ for $r\ge 1$, for which we show the (compact) embedding properties. When $H$ has subcritical and superlinear growth, we construct two frameworks, respectively with interpolation space method and dual method, to show the existence of nontrivial solution. As byproduct, we revisit the fractional Lane-Emden system, i.e. $H(u, v)=\frac{1}{p+1}|u|^{p+1}+\frac{1}{q+1}|v|^{q+1}$, and consider the existence, uniqueness of (radial) positive solutions under subcritical assumption.
△ Less
Submitted 17 August, 2024; v1 submitted 31 March, 2024;
originally announced April 2024.
-
Stable solution and extremal solution for the fractional $p$-Laplacian equation
Authors:
Weimin Zhang
Abstract:
To our knowledge, this paper is the first attempt to consider the stable solution and extremal solution for the fractional $p$-Laplacian equation: $(-Δ)_p^s u= λf(u),\; u> 0 ~\text{in}~Ω;\; u=0\;\text{in}~ \mathbb{R}^N\setminusΩ$, where $p>1$, $s\in (0,1)$, $N>sp$, $λ>0$ and $Ω$ is a bounded domain with continuous boundary. We first construct the notion of stable solution, and then we prove that w…
▽ More
To our knowledge, this paper is the first attempt to consider the stable solution and extremal solution for the fractional $p$-Laplacian equation: $(-Δ)_p^s u= λf(u),\; u> 0 ~\text{in}~Ω;\; u=0\;\text{in}~ \mathbb{R}^N\setminusΩ$, where $p>1$, $s\in (0,1)$, $N>sp$, $λ>0$ and $Ω$ is a bounded domain with continuous boundary. We first construct the notion of stable solution, and then we prove that when $f$ is of class $C^1$, nondecreasing and such that $f(0)>0$ and $\underset{t\to \infty}{\lim}\frac{f(t)}{t^{p-1}}=\infty$, there exists an extremal parameter $λ^*\in (0, \infty)$ such that a bounded minimal solution $u_λ$ exists if $λ\in (0, λ^*)$, and no bounded solution exists if $λ>λ^*$, no $W_0^{s,p}(Ω)$ solution exists if in addition $f(t)^{\frac{1}{p-1}}$ is convex. Moreover, this family of minimal solutions are stable, and nondecreasing in $λ$, therefore the extremal function $u^*:=\underset{λ\toλ^*}{\lim}u_λ$ exists.
For the regularity of the extremal function, we first show the $L^r$-estimates for the equation $(-Δ)_p^su=g$ with $g\in W_0^{s, p}(Ω)^*\cap L^q(Ω)$, $q\ge 1$. When $f$ is a power-like nonlinearity, we derive the $W_0^{s,p}(Ω)$ regularity of $u^*$ in all dimension and $L^{\infty}(Ω)$ regularity of $u^*$ in some low dimensions. For more general nonlinearities, when $f$ is class of $C^2$ and such that some convexity assumptions, then $u^*\in W_0^{s,p}(Ω)$ if $N<sp(1+\frac{p}{p-1})$ and $u^*\in L^{\infty}(Ω)$ if $N<\frac{sp^2}{p-1}$. Furthermore, when the limit $τ=\underset{t\to\infty}{\lim}\frac{f(t)f''(t)}{f'(t)^2}$ exists and $τ>\frac{p-2}{p-1}$, the results above can be improved as: $u^*\in W_0^{s,p}(Ω)$ for all dimensions, and $u^*\in L^{\infty}(Ω)$ if $N<sp+\frac{4sp}{p-1}$.
△ Less
Submitted 25 March, 2024;
originally announced March 2024.
-
Finite element method coupled with multiscale finite element method for the non-stationary Stokes-Darcy model
Authors:
Yachen Hong,
Wenhan Zhang,
Lina Zhao,
Haibiao Zheng
Abstract:
In this paper, we combine the multiscale flnite element method to propose an algorithm for solving the non-stationary Stokes-Darcy model, where the permeability coefflcient in the Darcy region exhibits multiscale characteristics. Our algorithm involves two steps: first, conducting the parallel computation of multiscale basis functions in the Darcy region. Second, based on these multiscale basis fu…
▽ More
In this paper, we combine the multiscale flnite element method to propose an algorithm for solving the non-stationary Stokes-Darcy model, where the permeability coefflcient in the Darcy region exhibits multiscale characteristics. Our algorithm involves two steps: first, conducting the parallel computation of multiscale basis functions in the Darcy region. Second, based on these multiscale basis functions, we employ an implicitexplicit scheme to solve the Stokes-Darcy equations. One signiflcant feature of the algorithm is that it solves problems on relatively coarse grids, thus signiflcantly reducing computational costs. Moreover, under the same coarse grid size, it exhibits higher accuracy compared to standard flnite element method. Under the assumption that the permeability coefflcient is periodic and independent of time, this paper demonstrates the stability and convergence of the algorithm. Finally, the rationality and effectiveness of the algorithm are verifled through three numerical experiments, with experimental results consistent with theoretical analysis.
△ Less
Submitted 18 March, 2024;
originally announced March 2024.
-
Convergence Rates For Tikhonov Regularization of Coefficient Identification Problems in Robin-Boundary Equation
Authors:
Huimin Huang,
Wensheng Zhang
Abstract:
This paper investigates the convergence rate for Tikhonov regularization of the problem of identifying the coefficient $a \in L^{\infty}(Ω)$ in the Robin-boundary equation $-\mathrm{div}(a\nabla u)-bu=f,~ x \in Ω\subset \mathbb R^M,~ M \geq 1$ and $u=0,~ x ~on~ \partialΩ$, where $f(x)\in L^{\infty}(Ω)$. Assume we only know the imprecise values of $u$ in the subset $Ω_1 \subset Ω$ given by…
▽ More
This paper investigates the convergence rate for Tikhonov regularization of the problem of identifying the coefficient $a \in L^{\infty}(Ω)$ in the Robin-boundary equation $-\mathrm{div}(a\nabla u)-bu=f,~ x \in Ω\subset \mathbb R^M,~ M \geq 1$ and $u=0,~ x ~on~ \partialΩ$, where $f(x)\in L^{\infty}(Ω)$. Assume we only know the imprecise values of $u$ in the subset $Ω_1 \subset Ω$ given by $z^δ \in {H}^1(Ω_1)$, satisfies $\|u-z^δ\|_{H^1(Ω_1)}\leq δ$. We assume $u$ satisfy the following boundary conditions on $\partialΩ_1$: \begin{align*} \nabla u \cdot \vec{n}+γu =0~on~\partialΩ_1, \end{align*} where $\vec{n}$ is the normal vector of $\partialΩ_1$ and $γ>0$ is a constant. We regularize this problem by correspondingly minimizing the strictly convex functional:
\begin{align*}
\min \limits_{a \in \mathbb A} &\frac12 \int_{Ω_1} a | {\nabla(U(a)-z^δ)}|^2 +\frac12\int_{\partialΩ_1} aγ[U(a)-z^δ]^2-\frac12 \int_{Ω_1} b [U(a)-z^δ]^2\\ &+ ρ\| a-a^* \|^2_{L^2(Ω)},
\end{align*}
where $U(a)$ is a map for $a$ to the solution of the Robin-boundary problem, $ρ> 0$ is the regularization parameter and $a^*$ is a priori estimate of $a$. We prove that the functional attain a unique global minimizer on the admissible set. Further, we give very simple source condition without the smallness requirement on the source function which provide the convergence rate $O(\sqrtδ)$ for the regularized solution.
△ Less
Submitted 15 March, 2024;
originally announced March 2024.
-
On the blowup of solutions for a nonlocal multi-dimensional transport equation
Authors:
Wanwan Zhang
Abstract:
In this paper, we revisit the problem of finite-time blowup for a multi-dimensional nonlocal transport equation studied in [Dong, Adv. Math. 264 (2014) 747-761]. Inspired by a one-dimensional analogous model considered in [Li-Rodrigo, Adv. Math. 374 (2020) 1-26], we establish a new weighted nonlinear inequality implying the blow-up by a completely real variable based technique.
In this paper, we revisit the problem of finite-time blowup for a multi-dimensional nonlocal transport equation studied in [Dong, Adv. Math. 264 (2014) 747-761]. Inspired by a one-dimensional analogous model considered in [Li-Rodrigo, Adv. Math. 374 (2020) 1-26], we establish a new weighted nonlinear inequality implying the blow-up by a completely real variable based technique.
△ Less
Submitted 11 March, 2024;
originally announced March 2024.
-
Model-free $H_{\infty}$ control of Itô stochastic system via off-policy reinforcement learning
Authors:
Jing Guo Jing Guo,
Xiushan Jiang,
Weihai Zhang
Abstract:
The stochastic $H_{\infty}$ control is studied for a linear stochastic Itô system with an unknown system model. The linear stochastic $H_{\infty}$ control issue is known to be transformable into the problem of solving a so-called generalized algebraic Riccati equation (GARE), which is a nonlinear equation that is typically difficult to solve analytically. Worse, model-based techniques cannot be ut…
▽ More
The stochastic $H_{\infty}$ control is studied for a linear stochastic Itô system with an unknown system model. The linear stochastic $H_{\infty}$ control issue is known to be transformable into the problem of solving a so-called generalized algebraic Riccati equation (GARE), which is a nonlinear equation that is typically difficult to solve analytically. Worse, model-based techniques cannot be utilized to approximately solve a GARE when an accurate system model is unavailable or prohibitively expensive to construct in reality. To address these issues, an off-policy reinforcement learning (RL) approach is presented to learn the solution of a GARE from real system data rather than a system model; its convergence is demonstrated, and the robustness of RL to errors in the learning process is investigated. In the off-policy RL approach, the system data may be created with behavior policies rather than the target policies, which is highly significant and promising for use in actual systems. Finally, the proposed off-policy RL approach is validated on a stochastic linear F-16 aircraft system.
△ Less
Submitted 7 March, 2024;
originally announced March 2024.
-
Sublinear expectation structure under finite states space
Authors:
Shuzhen Yang,
Wenqing Zhang
Abstract:
In this study, we propose the sublinear expectation structure under finite states space. To describe an interesting "nonlinear randomized" trial, based on a convex closed domain, we introduce a family of probability measures under finite states space. Corresponding the sublinear expectation operator introduced by S. Peng, we consider the related notation under finite states space. Within the finit…
▽ More
In this study, we propose the sublinear expectation structure under finite states space. To describe an interesting "nonlinear randomized" trial, based on a convex closed domain, we introduce a family of probability measures under finite states space. Corresponding the sublinear expectation operator introduced by S. Peng, we consider the related notation under finite states space. Within the finite states framework, the sublinear expectation can be explicitly calculated by a novel repeated summation formula, and some interesting examples are given. Furthermore, we establish Monotone convergence theorem, Fatou's lemma and Dominated convergence theorem under finite states space. Afterwards, we consider the independence under each probability measure, upon which we establish the nonlinear law of large numbers and obtain the maximal distribution under sublinear expectation.
△ Less
Submitted 7 March, 2024;
originally announced March 2024.
-
How Multimodal Integration Boost the Performance of LLM for Optimization: Case Study on Capacitated Vehicle Routing Problems
Authors:
Yuxiao Huang,
Wenjie Zhang,
Liang Feng,
Xingyu Wu,
Kay Chen Tan
Abstract:
Recently, large language models (LLMs) have notably positioned them as capable tools for addressing complex optimization challenges. Despite this recognition, a predominant limitation of existing LLM-based optimization methods is their struggle to capture the relationships among decision variables when relying exclusively on numerical text prompts, especially in high-dimensional problems. Keeping…
▽ More
Recently, large language models (LLMs) have notably positioned them as capable tools for addressing complex optimization challenges. Despite this recognition, a predominant limitation of existing LLM-based optimization methods is their struggle to capture the relationships among decision variables when relying exclusively on numerical text prompts, especially in high-dimensional problems. Keeping this in mind, we first propose to enhance the optimization performance using multimodal LLM capable of processing both textual and visual prompts for deeper insights of the processed optimization problem. This integration allows for a more comprehensive understanding of optimization problems, akin to human cognitive processes. We have developed a multimodal LLM-based optimization framework that simulates human problem-solving workflows, thereby offering a more nuanced and effective analysis. The efficacy of this method is evaluated through extensive empirical studies focused on a well-known combinatorial optimization problem, i.e., capacitated vehicle routing problem. The results are compared against those obtained from the LLM-based optimization algorithms that rely solely on textual prompts, demonstrating the significant advantages of our multimodal approach.
△ Less
Submitted 4 March, 2024;
originally announced March 2024.
-
High dimensional Gross--Zagier formula: a survey
Authors:
Wei Zhang
Abstract:
We survey recent developments on generalizing the Gross--Zagier formula to high dimensional Shimura varieties, with an emphasis on the
Arithmetic Gan--Gross--Prasad conjecture and the relative trace formula approach.
We survey recent developments on generalizing the Gross--Zagier formula to high dimensional Shimura varieties, with an emphasis on the
Arithmetic Gan--Gross--Prasad conjecture and the relative trace formula approach.
△ Less
Submitted 27 February, 2024;
originally announced February 2024.
-
An equation in nonlinear combination of iterates
Authors:
Chaitanya Gopalakrishna,
Weinian Zhang
Abstract:
In this paper we deal with an equation in nonlinear combination of iterates. Although it can be reduced by the logarithm conjugacy to a form for application of Schauder's or Banach's fixed point theorems, a difficulty called Zero Problem is encountered for continuous solutions because the domain does not contain $0$. So we consider solutions with weaker regularity, using the Knaster-Tarski fixed p…
▽ More
In this paper we deal with an equation in nonlinear combination of iterates. Although it can be reduced by the logarithm conjugacy to a form for application of Schauder's or Banach's fixed point theorems, a difficulty called Zero Problem is encountered for continuous solutions because the domain does not contain $0$. So we consider solutions with weaker regularity, using the Knaster-Tarski fixed point theorem for complete lattices to give order-preserving solutions. Then we give semi-continuous solutions and integrable solutions.
△ Less
Submitted 12 January, 2024;
originally announced January 2024.
-
Stability of Rarefaction Waves Under Periodic Perturbation for A Rate-Type Viscoelastic System
Authors:
lin Chang,
Duo Liu,
Weiqiang Zhang
Abstract:
In this paper, a rarefaction wave under space-periodic perturbation for the 3 times 3 rate-type viscoelastic system is considered. It is shown that if the initial perturbation around the rarefaction wave is suitably small, then the solution of the rate-type viscoelastic system tends to the rarefaction wave. The stability of solutions under periodic perturbations is an interesting and important pro…
▽ More
In this paper, a rarefaction wave under space-periodic perturbation for the 3 times 3 rate-type viscoelastic system is considered. It is shown that if the initial perturbation around the rarefaction wave is suitably small, then the solution of the rate-type viscoelastic system tends to the rarefaction wave. The stability of solutions under periodic perturbations is an interesting and important problem since the perturbation keeps oscillating at the far fields. That is, the perturbation is not integral in space. The key of proof is to construct a suitable ansatz carrying the same oscillation as the solution. Then we can find cancellations between solutions and ansatz such that the perturbation belongs to some Sobolev space. The nonlinear stability can be obtained by the weighted energy method.
△ Less
Submitted 8 January, 2024;
originally announced January 2024.
-
Multiscale finite element method for Stokes-Darcy model
Authors:
Yachen Hong,
Wenhan Zhang,
Lina Zhao,
Haibiao Zheng
Abstract:
This paper explores the application of the multiscale finite element method (MsFEM) to address steady-state Stokes-Darcy problems with BJS interface conditions in highly heterogeneous porous media. We assume the existence of multiscale features in the Darcy region and propose an algorithm for the multiscale Stokes-Darcy model. During the offline phase, we employ MsFEM to construct permeability-dep…
▽ More
This paper explores the application of the multiscale finite element method (MsFEM) to address steady-state Stokes-Darcy problems with BJS interface conditions in highly heterogeneous porous media. We assume the existence of multiscale features in the Darcy region and propose an algorithm for the multiscale Stokes-Darcy model. During the offline phase, we employ MsFEM to construct permeability-dependent offline bases for efficient coarse-grid simulation, with this process conducted in parallel to enhance its efficiency. In the online phase, we use the Robin-Robin algorithm to derive the model's solution. Subsequently, we conduct error analysis based on $L^2$ and $H^1$ norms, assuming certain periodic coefficients in the Darcy region. To validate our approach, we present extensive numerical tests on highly heterogeneous media, illustrating the results of the error analysis.
△ Less
Submitted 8 January, 2024;
originally announced January 2024.
-
Legendre-Moment Transform for Linear Ensemble Control and Computation
Authors:
Xin Ning,
Gong Cheng,
Wei Zhang,
Jr-Shin Li
Abstract:
Ensemble systems, pervasive in diverse scientific and engineering domains, pose challenges to existing control methods due to their massive scale and underactuated nature. This paper presents a dynamic moment approach to addressing theoretical and computational challenges in systems-theoretic analysis and control design for linear ensemble systems. We introduce the Legendre-moments and Legendre-mo…
▽ More
Ensemble systems, pervasive in diverse scientific and engineering domains, pose challenges to existing control methods due to their massive scale and underactuated nature. This paper presents a dynamic moment approach to addressing theoretical and computational challenges in systems-theoretic analysis and control design for linear ensemble systems. We introduce the Legendre-moments and Legendre-moment transform, which maps an ensemble system defined on the $L^2$-space to a Legendre-moment system defined on the $\ell^2$-space. We show that this pair of systems is of one-to-one correspondence and shares the same controllability property. This equivalence admits the control of an ensemble system through the control of the corresponding Legendre-moment system and inspires a unified control design scheme for linear ensemble systems using structured truncated moment systems. In particular, we develop a sampling-free ensemble control design algorithm, then conduct error analysis for control design using truncated moment systems and derive error bounds with respect to the truncation orders, which are illustrated with numerical examples.
△ Less
Submitted 3 January, 2024;
originally announced January 2024.
-
Normalized solutions for nonlinear Schrödinger equation involving potential and Sobolev critical exponent
Authors:
Zhen-Feng Jin,
Weimin Zhang
Abstract:
In this paper, we consider the existence of positive solutions with prescribed $L^2$-norm for the following nonlinear Schrödinger equation involving potential and Sobolev critical exponent \begin{equation*} \begin{cases} -Δu+V(x)u=λu+μ|u|^{p-2}u+|u|^{\frac{4}{N-2}}u \;\;\text { in } \mathbb{R}^N, \\ \|u\|_2=a>0,\\ \end{cases} \end{equation*} where $N\ge 3$, $μ>0$,…
▽ More
In this paper, we consider the existence of positive solutions with prescribed $L^2$-norm for the following nonlinear Schrödinger equation involving potential and Sobolev critical exponent \begin{equation*} \begin{cases} -Δu+V(x)u=λu+μ|u|^{p-2}u+|u|^{\frac{4}{N-2}}u \;\;\text { in } \mathbb{R}^N, \\ \|u\|_2=a>0,\\ \end{cases} \end{equation*} where $N\ge 3$, $μ>0$, $p\in [2+\frac{4}{N}, \frac{2N}{N-2})$ and $V\in C^1(\mathbb{R}^N)$. Under different assumptions on $V$, we derive two different Pohozaev identities. Based on these two cases, we respectively obtain the existence of positive solution. As far as we are aware, we did not find any works on normalized solutions with Sobolev critical growth and potential $V \not\equiv 0$. Our results extend some results of Wei and Wu [J. Funct. Anal. 283(2022)] to the potential case.
△ Less
Submitted 26 December, 2023;
originally announced December 2023.
-
The Rate-Distortion-Perception-Classification Tradeoff: Joint Source Coding and Modulation via Inverse-Domain GANs
Authors:
Junli Fang,
João F. C. Mota,
Baoshan Lu,
Weicheng Zhang,
Xuemin Hong
Abstract:
The joint source-channel coding (JSCC) framework leverages deep learning to learn from data the best codes for source and channel coding. When the output signal, rather than being binary, is directly mapped onto the IQ domain (complex-valued), we call the resulting framework joint source coding and modulation (JSCM). We consider a JSCM scenario and show the existence of a strict tradeoff between c…
▽ More
The joint source-channel coding (JSCC) framework leverages deep learning to learn from data the best codes for source and channel coding. When the output signal, rather than being binary, is directly mapped onto the IQ domain (complex-valued), we call the resulting framework joint source coding and modulation (JSCM). We consider a JSCM scenario and show the existence of a strict tradeoff between channel rate, distortion, perception, and classification accuracy, a tradeoff that we name RDPC. We then propose two image compression methods to navigate that tradeoff: the RDPCO algorithm which, under simple assumptions, directly solves the optimization problem characterizing the tradeoff, and an algorithm based on an inverse-domain generative adversarial network (ID-GAN), which is more general and achieves extreme compression. Simulation results corroborate the theoretical findings, showing that both algorithms exhibit the RDPC tradeoff. They also demonstrate that the proposed ID-GAN algorithm effectively balances image distortion, perception, and classification accuracy, and significantly outperforms traditional separation-based methods and recent deep JSCM architectures in terms of one or more of these metrics.
△ Less
Submitted 6 June, 2024; v1 submitted 22 December, 2023;
originally announced December 2023.
-
Snake Graphs for Graph LP Algebras
Authors:
Esther Banaian,
Sunita Chepuri,
Elizabeth Kelley,
Sylvester W. Zhang
Abstract:
Graph LP algebras are a generalization of cluster algebras introduced by Lam and Pylyavskyy. We provide a combinatorial proof of positivity for certain cluster variables in these algebras. This proof uses a hypergraph generalization of snake graphs, a class of planar graphs which were used by Musiker, Schiffler, and Williams to prove positivity for cluster algebras from surfaces. These results ext…
▽ More
Graph LP algebras are a generalization of cluster algebras introduced by Lam and Pylyavskyy. We provide a combinatorial proof of positivity for certain cluster variables in these algebras. This proof uses a hypergraph generalization of snake graphs, a class of planar graphs which were used by Musiker, Schiffler, and Williams to prove positivity for cluster algebras from surfaces. These results extend those given in our previous paper, where we used a related combinatorial object known as a $T$-path.
△ Less
Submitted 19 December, 2023;
originally announced December 2023.
-
Braid group action and quasi-split affine $\imath$quantum groups II: higher rank
Authors:
Ming Lu,
Weiqiang Wang,
Weinan Zhang
Abstract:
This paper studies quantum symmetric pairs $(\widetilde{\mathbf U}, \widetilde{\mathbf U}^\imath )$ associated with quasi-split Satake diagrams of affine type $A_{2r-1}, D_r, E_{6}$ with a nontrivial diagram involution fixing the affine simple node. Various real and imaginary root vectors for the universal $\imath$quantum groups $\widetilde{\mathbf U}^\imath$ are constructed with the help of the r…
▽ More
This paper studies quantum symmetric pairs $(\widetilde{\mathbf U}, \widetilde{\mathbf U}^\imath )$ associated with quasi-split Satake diagrams of affine type $A_{2r-1}, D_r, E_{6}$ with a nontrivial diagram involution fixing the affine simple node. Various real and imaginary root vectors for the universal $\imath$quantum groups $\widetilde{\mathbf U}^\imath$ are constructed with the help of the relative braid group action, and they are used to construct affine rank one subalgebras of $\widetilde{\mathbf U}^\imath$. We then establish relations among real and imaginary root vectors in different affine rank one subalgebras and use them to give a Drinfeld type presentation of $\widetilde{\mathbf U}^\imath$.
△ Less
Submitted 29 March, 2024; v1 submitted 16 November, 2023;
originally announced November 2023.
-
Stochastic convergence of regularized solutions for backward heat conduction problems
Authors:
Zhongjian Wang,
Wenlong Zhang,
Zhiwen Zhang
Abstract:
In this paper, we study the stochastic convergence of regularized solutions for backward heat conduction problems. These problems are recognized as ill-posed due to the exponential decay of eigenvalues associated with the forward problems. We derive an error estimate for the least-squares regularized minimization problem within the framework of stochastic convergence. Our analysis reveals that the…
▽ More
In this paper, we study the stochastic convergence of regularized solutions for backward heat conduction problems. These problems are recognized as ill-posed due to the exponential decay of eigenvalues associated with the forward problems. We derive an error estimate for the least-squares regularized minimization problem within the framework of stochastic convergence. Our analysis reveals that the optimal error of the Tikhonov-type least-squares optimization problem depends on the noise level, the number of sensors, and the underlying ground truth. Moreover, we propose a self-adaptive algorithm to identify the optimal regularization parameter for the optimization problem without requiring knowledge of the noise level or any other prior information, which will be very practical in applications. We present numerical examples to demonstrate the accuracy and efficiency of our proposed method. These numerical results show that our method is efficient in solving backward heat conduction problems.
△ Less
Submitted 6 November, 2023;
originally announced November 2023.
-
Commutators for certain fractional type operators on weighted spaces and Orlicz-Morrey spaces
Authors:
Huoxiong Wu Tong Zhang
Abstract:
In this paper, we focus on a class of fractional type integral operators that can be served as extensions of Riesz potential with kernels $$K(x,y)=\frac{Ω_1(x-A_1 y)}{|x-A_1 y |^{\frac{n}{q_1}}} \cdots \frac{Ω_m(x-A_m y)}{|x-A_m y |^{\frac{n}{q_m}}},$$ where $α\in [0,n), m\geqslant1, \sum_{i=1}^m\frac{n}{q_i}=n-α$, $\{A_i\}^m_{i=1}$ are invertible matrixes, $Ω_i$ is homogeneous of degree 0 on…
▽ More
In this paper, we focus on a class of fractional type integral operators that can be served as extensions of Riesz potential with kernels $$K(x,y)=\frac{Ω_1(x-A_1 y)}{|x-A_1 y |^{\frac{n}{q_1}}} \cdots \frac{Ω_m(x-A_m y)}{|x-A_m y |^{\frac{n}{q_m}}},$$ where $α\in [0,n), m\geqslant1, \sum_{i=1}^m\frac{n}{q_i}=n-α$, $\{A_i\}^m_{i=1}$ are invertible matrixes, $Ω_i$ is homogeneous of degree 0 on $\R^n$ and $Ω_i\in L^{p_i}(S^{n-1})$ for some $p_i\in [1,\infty)$. Under appropriate assumptions, we obtain the weighted $L^p$ estimates as well as weighted Hardy estimates of the commutator for such operators with $BMO$-type function. In addition, we acquire the boundedness of these operators and their commutators with a function in Campanato space on Orcliz-Morrey spaces as well as the compactness for such commutators in a special case: $m=1$ and $A=I$.
△ Less
Submitted 5 November, 2023;
originally announced November 2023.
-
Vanishing results for the modified diagonal cycles II: Shimura curves
Authors:
Congling Qiu,
Wei Zhang
Abstract:
We prove vanishing results for the modified diagonal cycles in the Chow groups of the triple products of Shimura curves and their motivic direct summands. In particular we find examples of curves with trivial automorphism groups and vanishing modified diagonal cycles.
We prove vanishing results for the modified diagonal cycles in the Chow groups of the triple products of Shimura curves and their motivic direct summands. In particular we find examples of curves with trivial automorphism groups and vanishing modified diagonal cycles.
△ Less
Submitted 30 October, 2023;
originally announced October 2023.