-
Algebraic solitons in the massive Thirring model
Authors:
Jiaqi Han,
Cheng He,
Dmitry E. Pelinovsky
Abstract:
We present exact solutions describing dynamics of two algebraic solitons in the massive Thirring model. Each algebraic soliton corresponds to a simple embedded eigenvalue in the Kaup--Newell spectral problem and attains the maximal mass among the family of solitary waves traveling with the same speed. By coalescence of speeds of the two algebraic solitons, we find a new solution for an algebraic d…
▽ More
We present exact solutions describing dynamics of two algebraic solitons in the massive Thirring model. Each algebraic soliton corresponds to a simple embedded eigenvalue in the Kaup--Newell spectral problem and attains the maximal mass among the family of solitary waves traveling with the same speed. By coalescence of speeds of the two algebraic solitons, we find a new solution for an algebraic double-soliton which corresponds to a double embedded eigenvalue. We show that the double-soliton attains the double mass of a single soliton and describes a slow interaction of two identical algebraic solitons.
△ Less
Submitted 10 June, 2024;
originally announced June 2024.
-
Riemannian Accelerated Zeroth-order Algorithm: Improved Robustness and Lower Query Complexity
Authors:
Chang He,
Zhaoye Pan,
Xiao Wang,
Bo Jiang
Abstract:
Optimization problems with access to only zeroth-order information of the objective function on Riemannian manifolds arise in various applications, spanning from statistical learning to robot learning. While various zeroth-order algorithms have been proposed in Euclidean space, they are not inherently designed to handle the challenging constraints imposed by Riemannian manifolds. The proper adapta…
▽ More
Optimization problems with access to only zeroth-order information of the objective function on Riemannian manifolds arise in various applications, spanning from statistical learning to robot learning. While various zeroth-order algorithms have been proposed in Euclidean space, they are not inherently designed to handle the challenging constraints imposed by Riemannian manifolds. The proper adaptation of zeroth-order techniques to Riemannian manifolds remained unknown until the pioneering work of \cite{li2023stochastic}. However, zeroth-order algorithms are widely observed to converge slowly and be unstable in practice. To alleviate these issues, we propose a Riemannian accelerated zeroth-order algorithm with improved robustness. Regarding efficiency, our accelerated algorithm has the function query complexity of $\mathcal{O}(ε^{-7/4}d)$ for finding an $ε$-approximate first-order stationary point. By introducing a small perturbation, it exhibits a function query complexity of $\tilde{\mathcal{O}}(ε^{-7/4}d)$ for seeking a second-order stationary point with a high probability, matching state-of-the-art result in Euclidean space. Moreover, we further establish the almost sure convergence in the asymptotic sense through the Stable Manifold Theorem. Regarding robustness, our algorithm requires larger smoothing parameters in the order of $\tilde{\mathcal{O}}(ε^{7/8}d^{-1/2})$, improving the existing result by a factor of $\tilde{\mathcal{O}}(ε^{3/4})$.
△ Less
Submitted 9 May, 2024;
originally announced May 2024.
-
Newton-CG methods for nonconvex unconstrained optimization with Hölder continuous Hessian
Authors:
Chuan He,
Zhaosong Lu
Abstract:
In this paper we consider a nonconvex unconstrained optimization problem minimizing a twice differentiable objective function with Hölder continuous Hessian. Specifically, we first propose a Newton-conjugate gradient (Newton-CG) method for finding an approximate first-order stationary point (FOSP) of this problem, assuming the associated the Hölder parameters are explicitly known. Then we develop…
▽ More
In this paper we consider a nonconvex unconstrained optimization problem minimizing a twice differentiable objective function with Hölder continuous Hessian. Specifically, we first propose a Newton-conjugate gradient (Newton-CG) method for finding an approximate first-order stationary point (FOSP) of this problem, assuming the associated the Hölder parameters are explicitly known. Then we develop a parameter-free Newton-CG method without requiring any prior knowledge of these parameters. To the best of our knowledge, this method is the first parameter-free second-order method achieving the best-known iteration and operation complexity for finding an approximate FOSP of this problem. Furthermore, we propose a Newton-CG method for finding an approximate second-order stationary point (SOSP) of the considered problem with high probability and establish its iteration and operation complexity. Finally, we present preliminary numerical results to demonstrate the superior practical performance of our parameter-free Newton-CG method over a well-known regularized Newton method.
△ Less
Submitted 21 November, 2023;
originally announced November 2023.
-
A Universal Trust-Region Method for Convex and Nonconvex Optimization
Authors:
Yuntian Jiang,
Chang He,
Chuwen Zhang,
Dongdong Ge,
Bo Jiang,
Yinyu Ye
Abstract:
This paper presents a universal trust-region method simultaneously incorporating quadratic regularization and the ball constraint. We introduce a novel mechanism to set the parameters in the proposed method that unifies the analysis for convex and nonconvex optimization. Our method exhibits an iteration complexity of $\tilde O(ε^{-3/2})$ to find an approximate second-order stationary point for non…
▽ More
This paper presents a universal trust-region method simultaneously incorporating quadratic regularization and the ball constraint. We introduce a novel mechanism to set the parameters in the proposed method that unifies the analysis for convex and nonconvex optimization. Our method exhibits an iteration complexity of $\tilde O(ε^{-3/2})$ to find an approximate second-order stationary point for nonconvex optimization. Meanwhile, the analysis reveals that the universal method attains an $O(ε^{-1/2})$ complexity bound for convex optimization and can be accelerated. These results are complementary to the existing literature as the trust-region method was historically conceived for nonconvex optimization. Finally, we develop an adaptive universal method to address practical implementations. The numerical results show the effectiveness of our method in both nonconvex and convex problems.
△ Less
Submitted 12 March, 2024; v1 submitted 19 November, 2023;
originally announced November 2023.
-
Federated Learning with Convex Global and Local Constraints
Authors:
Chuan He,
Le Peng,
Ju Sun
Abstract:
In practice, many machine learning (ML) problems come with constraints, and their applied domains involve distributed sensitive data that cannot be shared with others, e.g., in healthcare. Collaborative learning in such practical scenarios entails federated learning (FL) for ML problems with constraints, or FL with constraints for short. Despite the extensive developments of FL techniques in recen…
▽ More
In practice, many machine learning (ML) problems come with constraints, and their applied domains involve distributed sensitive data that cannot be shared with others, e.g., in healthcare. Collaborative learning in such practical scenarios entails federated learning (FL) for ML problems with constraints, or FL with constraints for short. Despite the extensive developments of FL techniques in recent years, these techniques only deal with unconstrained FL problems or FL problems with simple constraints that are amenable to easy projections. There is little work dealing with FL problems with general constraints. To fill this gap, we take the first step toward building an algorithmic framework for solving FL problems with general constraints. In particular, we propose a new FL algorithm for constrained ML problems based on the proximal augmented Lagrangian (AL) method. Assuming convex objective and convex constraints plus other mild conditions, we establish the worst-case complexity of the proposed algorithm. Our numerical experiments show the effectiveness of our algorithm in performing Neyman-Pearson classification and fairness-aware learning with nonconvex constraints, in an FL setting.
△ Less
Submitted 1 May, 2024; v1 submitted 16 October, 2023;
originally announced October 2023.
-
On the sum of the first two largest signless Laplacian eigenvalues of a graph
Authors:
Zi-Ming Zhou,
Chang-Xiang He,
Hai-Ying Shan
Abstract:
For a graph $G$, let $S_2(G)$ be the sum of the first two largest signless Laplacian eigenvalues of $G$, and $f(G)=e(G)+3-S_2(G)$. Oliveira, Lima, Rama and Carvalho conjectured that $K^+_{1,n-1}$ (the star graph with an additional edge) is the unique graph with minimum value of $f(G)$ on $n$ vertices. In this paper, we prove this conjecture, which also confirm a conjecture for the upper bound of…
▽ More
For a graph $G$, let $S_2(G)$ be the sum of the first two largest signless Laplacian eigenvalues of $G$, and $f(G)=e(G)+3-S_2(G)$. Oliveira, Lima, Rama and Carvalho conjectured that $K^+_{1,n-1}$ (the star graph with an additional edge) is the unique graph with minimum value of $f(G)$ on $n$ vertices. In this paper, we prove this conjecture, which also confirm a conjecture for the upper bound of $S_2(G)$ proposed by Ashraf et al.
△ Less
Submitted 13 June, 2024; v1 submitted 13 October, 2023;
originally announced October 2023.
-
Massive Thirring Model: Inverse Scattering and Soliton Resolution
Authors:
Cheng He,
Jiaqi Liu,
Changzheng Qu
Abstract:
In this paper the long-time dynamics of the massive Thirring model is investigated. Firstly the nonlinear steepest descent method for Riemann-Hilbert problem is explored to obtain the soliton resolution of the solutions to the massive Thirring model whose initial data belong to some weighted-Sobolev spaces. Secondly, the asymptotic stability of multi-solitons follow as a corollary. The main diffic…
▽ More
In this paper the long-time dynamics of the massive Thirring model is investigated. Firstly the nonlinear steepest descent method for Riemann-Hilbert problem is explored to obtain the soliton resolution of the solutions to the massive Thirring model whose initial data belong to some weighted-Sobolev spaces. Secondly, the asymptotic stability of multi-solitons follow as a corollary. The main difficulty in studying the massive Thirring model through inverse scattering is that the corresponding Lax pair has singularities at the origin and infinity. We overcome this difficulty by making use of two transforms that separate the singularities.
△ Less
Submitted 28 July, 2023;
originally announced July 2023.
-
Best approximation results and essential boundary conditions for novel types of weak adversarial network discretizations for PDEs
Authors:
Silvia Bertoluzza,
Erik Burman,
Cuiyu He
Abstract:
In this paper, we provide a theoretical analysis of the recently introduced weakly adversarial networks (WAN) method, used to approximate partial differential equations in high dimensions. We address the existence and stability of the solution, as well as approximation bounds. More precisely, we prove the existence of discrete solutions, intended in a suitable weak sense, for which we prove a quas…
▽ More
In this paper, we provide a theoretical analysis of the recently introduced weakly adversarial networks (WAN) method, used to approximate partial differential equations in high dimensions. We address the existence and stability of the solution, as well as approximation bounds. More precisely, we prove the existence of discrete solutions, intended in a suitable weak sense, for which we prove a quasi-best approximation estimate similar to Cea's lemma, a result commonly found in finite element methods. We also propose two new stabilized WAN-based formulas that avoid the need for direct normalization. Furthermore, we analyze the method's effectiveness for the Dirichlet boundary problem that employs the implicit representation of the geometry. The key requirement for achieving the best approximation outcome is to ensure that the space for the test network satisfies a specific condition, known as the inf-sup condition, essentially requiring that the test network set is sufficiently large when compared to the trial space. The method's accuracy, however, is only determined by the space of the trial network. We also devise a pseudo-time XNODE neural network class for static PDE problems, yielding significantly faster convergence results than the classical DNN network.
△ Less
Submitted 29 January, 2024; v1 submitted 11 July, 2023;
originally announced July 2023.
-
Homogeneous Second-Order Descent Framework: A Fast Alternative to Newton-Type Methods
Authors:
Chang He,
Yuntian Jiang,
Chuwen Zhang,
Dongdong Ge,
Bo Jiang,
Yinyu Ye
Abstract:
This paper proposes a homogeneous second-order descent framework (HSODF) for nonconvex and convex optimization based on the generalized homogeneous model (GHM). In comparison to the Newton steps, the GHM can be solved by extremal symmetric eigenvalue procedures and thus grant an advantage in ill-conditioned problems. Moreover, GHM extends the ordinary homogeneous model (OHM) (Zhang et al. 2022) to…
▽ More
This paper proposes a homogeneous second-order descent framework (HSODF) for nonconvex and convex optimization based on the generalized homogeneous model (GHM). In comparison to the Newton steps, the GHM can be solved by extremal symmetric eigenvalue procedures and thus grant an advantage in ill-conditioned problems. Moreover, GHM extends the ordinary homogeneous model (OHM) (Zhang et al. 2022) to allow adaptiveness in the construction of the aggregated matrix. Consequently, HSODF is able to recover some well-known second-order methods, such as trust-region methods and gradient regularized methods, while maintaining comparable iteration complexity bounds. We also study two specific realizations of HSODF. One is adaptive HSODM, which has a parameter-free $O(ε^{-3/2})$ global complexity bound for nonconvex second-order Lipschitz continuous objective functions. The other one is homotopy HSODM, which is proven to have a global linear rate of convergence without strong convexity. The efficiency of our approach to ill-conditioned and high-dimensional problems is justified by some preliminary numerical results.
△ Less
Submitted 6 May, 2024; v1 submitted 30 June, 2023;
originally announced June 2023.
-
A Parameter-Free Conditional Gradient Method for Composite Minimization under Hölder Condition
Authors:
Masaru Ito,
Zhaosong Lu,
Chuan He
Abstract:
In this paper we consider a composite optimization problem that minimizes the sum of a weakly smooth function and a convex function with either a bounded domain or a uniformly convex structure. In particular, we first present a parameter-dependent conditional gradient method for this problem, whose step sizes require prior knowledge of the parameters associated with the Hölder continuity of the gr…
▽ More
In this paper we consider a composite optimization problem that minimizes the sum of a weakly smooth function and a convex function with either a bounded domain or a uniformly convex structure. In particular, we first present a parameter-dependent conditional gradient method for this problem, whose step sizes require prior knowledge of the parameters associated with the Hölder continuity of the gradient of the weakly smooth function, and establish its rate of convergence. Given that these parameters could be unknown or known but possibly conservative, such a method may suffer from implementation issue or slow convergence. We therefore propose a parameter-free conditional gradient method whose step size is determined by using a constructive local quadratic upper approximation and an adaptive line search scheme, without using any problem parameter. We show that this method achieves the same rate of convergence as the parameter-dependent conditional gradient method. Preliminary experiments are also conducted and illustrate the superior performance of the parameter-free conditional gradient method over the methods with some other step size rules.
△ Less
Submitted 29 May, 2023;
originally announced May 2023.
-
Scarf's algorithm and stable marriages
Authors:
Yuri Faenza,
Chengyue He,
Jay Sethuraman
Abstract:
Scarf's algorithm gives a pivoting procedure to find a special vertex -- a dominating vertex -- in down-monotone polytopes. This paper studies the behavior of Scarf's algorithm when employed to find stable matchings in bipartite graphs. First, it proves that Scarf's algorithm can be implemented to run in polynomial time, showing the first positive result on its runtime in significant settings. Sec…
▽ More
Scarf's algorithm gives a pivoting procedure to find a special vertex -- a dominating vertex -- in down-monotone polytopes. This paper studies the behavior of Scarf's algorithm when employed to find stable matchings in bipartite graphs. First, it proves that Scarf's algorithm can be implemented to run in polynomial time, showing the first positive result on its runtime in significant settings. Second, it shows an infinite family of instances where, no matter the pivoting rule and runtime, Scarf's algorithm outputs a matching from an exponentially small subset of all stable matchings, thus showing a structural weakness of the approach.
△ Less
Submitted 1 March, 2023;
originally announced March 2023.
-
A Newton-CG based barrier-augmented Lagrangian method for general nonconvex conic optimization
Authors:
Chuan He,
Heng Huang,
Zhaosong Lu
Abstract:
In this paper we consider finding an approximate second-order stationary point (SOSP) of general nonconvex conic optimization that minimizes a twice differentiable function subject to nonlinear equality constraints and also a convex conic constraint. In particular, we propose a Newton-conjugate gradient (Newton-CG) based barrier-augmented Lagrangian method for finding an approximate SOSP of this p…
▽ More
In this paper we consider finding an approximate second-order stationary point (SOSP) of general nonconvex conic optimization that minimizes a twice differentiable function subject to nonlinear equality constraints and also a convex conic constraint. In particular, we propose a Newton-conjugate gradient (Newton-CG) based barrier-augmented Lagrangian method for finding an approximate SOSP of this problem. Under some mild assumptions, we show that our method enjoys a total inner iteration complexity of $\widetilde{\cal O}(ε^{-11/2})$ and an operation complexity of $\widetilde{\cal O}(ε^{-11/2}\min\{n,ε^{-5/4}\})$ for finding an $(ε,\sqrtε)$-SOSP of general nonconvex conic optimization with high probability. Moreover, under a constraint qualification, these complexity bounds are improved to $\widetilde{\cal O}(ε^{-7/2})$ and $\widetilde{\cal O}(ε^{-7/2}\min\{n,ε^{-3/4}\})$, respectively. To the best of our knowledge, this is the first study on the complexity of finding an approximate SOSP of general nonconvex conic optimization. Preliminary numerical results are presented to demonstrate superiority of the proposed method over first-order methods in terms of solution quality.
△ Less
Submitted 10 January, 2023;
originally announced January 2023.
-
A Newton-CG based augmented Lagrangian method for finding a second-order stationary point of nonconvex equality constrained optimization with complexity guarantees
Authors:
Chuan He,
Zhaosong Lu,
Ting Kei Pong
Abstract:
In this paper we consider finding a second-order stationary point (SOSP) of nonconvex equality constrained optimization when a nearly feasible point is known. In particular, we first propose a new Newton-CG method for finding an approximate SOSP of unconstrained optimization and show that it enjoys a substantially better complexity than the Newton-CG method [56]. We then propose a Newton-CG based…
▽ More
In this paper we consider finding a second-order stationary point (SOSP) of nonconvex equality constrained optimization when a nearly feasible point is known. In particular, we first propose a new Newton-CG method for finding an approximate SOSP of unconstrained optimization and show that it enjoys a substantially better complexity than the Newton-CG method [56]. We then propose a Newton-CG based augmented Lagrangian (AL) method for finding an approximate SOSP of nonconvex equality constrained optimization, in which the proposed Newton-CG method is used as a subproblem solver. We show that under a generalized linear independence constraint qualification (GLICQ), our AL method enjoys a total inner iteration complexity of $\widetilde{\cal O}(ε^{-7/2})$ and an operation complexity of $\widetilde{\cal O}(ε^{-7/2}\min\{n,ε^{-3/4}\})$ for finding an $(ε,\sqrtε)$-SOSP of nonconvex equality constrained optimization with high probability, which are significantly better than the ones achieved by the proximal AL method [60]. Besides, we show that it has a total inner iteration complexity of $\widetilde{\cal O}(ε^{-11/2})$ and an operation complexity of $\widetilde{\cal O}(ε^{-11/2}\min\{n,ε^{-5/4}\})$ when the GLICQ does not hold. To the best of our knowledge, all the complexity results obtained in this paper are new for finding an approximate SOSP of nonconvex equality constrained optimization with high probability. Preliminary numerical results also demonstrate the superiority of our proposed methods over the ones in [56,60].
△ Less
Submitted 8 January, 2023;
originally announced January 2023.
-
Perturbation results for distance-edge-monitoring numbers
Authors:
Chenxu Yang,
Ralf Klasing,
Changxiang He,
Yaping Mao
Abstract:
Foucaud et al. recently introduced and initiated the study of a new graph-theoretic concept in the area of network monitoring. Given a graph $G=(V(G), E(G))$, a set $M \subseteq V(G)$ is a distance-edge-monitoring set if for every edge $e \in E(G)$, there is a vertex $x \in M$ and a vertex $y \in V(G)$ such that the edge $e$ belongs to all shortest paths between $x$ and $y$. The smallest size of s…
▽ More
Foucaud et al. recently introduced and initiated the study of a new graph-theoretic concept in the area of network monitoring. Given a graph $G=(V(G), E(G))$, a set $M \subseteq V(G)$ is a distance-edge-monitoring set if for every edge $e \in E(G)$, there is a vertex $x \in M$ and a vertex $y \in V(G)$ such that the edge $e$ belongs to all shortest paths between $x$ and $y$. The smallest size of such a set in $G$ is denoted by $\operatorname{dem}(G)$. Denoted by $G-e$ (resp. $G \backslash u$) the subgraph of $G$ obtained by removing the edge $e$ from $G$ (resp. a vertex $u$ together with all its incident edges from $G$). In this paper, we first show that $\operatorname{dem}(G-e)- \operatorname{dem}(G)\leq 2$ for any graph $G$ and edge $e \in E(G)$. Moreover, the bound is sharp. Next, we construct two graphs $G$ and $H$ to show that $\operatorname{dem}(G)-\operatorname{dem}(G\setminus u)$ and $\operatorname{dem}(H\setminus v)-\operatorname{dem}(H)$ can be arbitrarily large, where $u \in V(G)$ and $v \in V(H)$. We also study the relation between $\operatorname{dem}(H)$ and $\operatorname{dem}(G)$, where $H$ is a subgraph of $G$. In the end, we give an algorithm to judge whether the distance-edge monitoring set still remain in the resulting graph when any edge of the graph $G$ is deleted.
△ Less
Submitted 10 May, 2024; v1 submitted 6 January, 2023;
originally announced January 2023.
-
A Homogeneous Second-Order Descent Method for Nonconvex Optimization
Authors:
Chuwen Zhang,
Dongdong Ge,
Chang He,
Bo Jiang,
Yuntian Jiang,
Chenyu Xue,
Yinyu Ye
Abstract:
In this paper, we introduce a Homogeneous Second-Order Descent Method (HSODM) using the homogenized quadratic approximation to the original function. The merit of homogenization is that only the leftmost eigenvector of a gradient-Hessian integrated matrix is computed at each iteration. Therefore, the algorithm is a single-loop method that does not need to switch to other sophisticated algorithms a…
▽ More
In this paper, we introduce a Homogeneous Second-Order Descent Method (HSODM) using the homogenized quadratic approximation to the original function. The merit of homogenization is that only the leftmost eigenvector of a gradient-Hessian integrated matrix is computed at each iteration. Therefore, the algorithm is a single-loop method that does not need to switch to other sophisticated algorithms and is easy to implement. We show that HSODM has a global convergence rate of $O(ε^{-3/2})$ to find an $ε$-approximate second-order stationary point, and has a local quadratic convergence rate under the standard assumptions. The numerical results demonstrate the advantage of the proposed method over other second-order methods.
△ Less
Submitted 6 May, 2024; v1 submitted 15 November, 2022;
originally announced November 2022.
-
DRSOM: A Dimension Reduced Second-Order Method
Authors:
Chuwen Zhang,
Dongdong Ge,
Chang He,
Bo Jiang,
Yuntian Jiang,
Yinyu Ye
Abstract:
In this paper, we propose a Dimension-Reduced Second-Order Method (DRSOM) for convex and nonconvex (unconstrained) optimization. Under a trust-region-like framework, our method preserves the convergence of the second-order method while using only curvature information in a few directions. Consequently, the computational overhead of our method remains comparable to the first-order such as the gradi…
▽ More
In this paper, we propose a Dimension-Reduced Second-Order Method (DRSOM) for convex and nonconvex (unconstrained) optimization. Under a trust-region-like framework, our method preserves the convergence of the second-order method while using only curvature information in a few directions. Consequently, the computational overhead of our method remains comparable to the first-order such as the gradient descent method. Theoretically, we show that the method has a local quadratic convergence and a global convergence rate of $O(ε^{-3/2})$ to satisfy the first-order and second-order conditions if the subspace satisfies a commonly adopted approximated Hessian assumption. We further show that this assumption can be removed if we perform a corrector step using a Krylov-like method periodically at the end stage of the algorithm. The applicability and performance of DRSOM are exhibited by various computational experiments, including $L_2 - L_p$ minimization, CUTEst problems, and sensor network localization.
△ Less
Submitted 2 July, 2023; v1 submitted 30 July, 2022;
originally announced August 2022.
-
A Newton-CG based barrier method for finding a second-order stationary point of nonconvex conic optimization with complexity guarantees
Authors:
Chuan He,
Zhaosong Lu
Abstract:
In this paper we consider finding an approximate second-order stationary point (SOSP) of nonconvex conic optimization that minimizes a twice differentiable function over the intersection of an affine subspace and a convex cone. In particular, we propose a Newton-conjugate gradient (Newton-CG) based barrier method for finding an $(ε,\sqrtε)$-SOSP of this problem. Our method is not only implementabl…
▽ More
In this paper we consider finding an approximate second-order stationary point (SOSP) of nonconvex conic optimization that minimizes a twice differentiable function over the intersection of an affine subspace and a convex cone. In particular, we propose a Newton-conjugate gradient (Newton-CG) based barrier method for finding an $(ε,\sqrtε)$-SOSP of this problem. Our method is not only implementable, but also achieves an iteration complexity of ${\cal O}(ε^{-3/2})$, which matches the best known iteration complexity of second-order methods for finding an $(ε,\sqrtε)$-SOSP of unconstrained nonconvex optimization. The operation complexity, consisting of ${\cal O}(ε^{-3/2})$ Cholesky factorizations and $\widetilde{\cal O}(ε^{-3/2}\min\{n,ε^{-1/4}\})$ other fundamental operations, is also established for our method.
△ Less
Submitted 11 October, 2022; v1 submitted 12 July, 2022;
originally announced July 2022.
-
Normalized solutions for Schrödinger-Bopp-Podolsky system
Authors:
Chuan-Min He,
Lin Li,
Shang-Jie Chen
Abstract:
In this paper, we study the following energy functional originates from the Schrödinger-Bopp-Podolsky system $$I(u)=\frac{1}{2}\int_{\mathbb{R}^{3}}|\nabla u|^{2}dx+\frac{1}{4}\int_{\mathbb{R}^{3}} φ_{u}u^{2}dx-\frac{1}{p}\int_{\mathbb{R}^{3}}|u|^{p}dx$$ constrained on $B_ρ=\left\{u\in H^{1}(\mathbb{R}^{3},C):\ \left\|u\right\|_{2}=ρ\right\},$ where $ρ>0.$ As such constrained problem $I(u)$ is bou…
▽ More
In this paper, we study the following energy functional originates from the Schrödinger-Bopp-Podolsky system $$I(u)=\frac{1}{2}\int_{\mathbb{R}^{3}}|\nabla u|^{2}dx+\frac{1}{4}\int_{\mathbb{R}^{3}} φ_{u}u^{2}dx-\frac{1}{p}\int_{\mathbb{R}^{3}}|u|^{p}dx$$ constrained on $B_ρ=\left\{u\in H^{1}(\mathbb{R}^{3},C):\ \left\|u\right\|_{2}=ρ\right\},$ where $ρ>0.$ As such constrained problem $I(u)$ is bounded from below on $B_ρ$ when $p\in(2,\frac{10}{3}).$ We use minimizing method to get a normalized solution.
△ Less
Submitted 8 June, 2022;
originally announced June 2022.
-
Decision Rule Approaches for Pessimistic Bilevel Linear Programs under Moment Ambiguity with Facility Location Applications
Authors:
Akshit Goyal,
Yiling Zhang,
Chuan He
Abstract:
We study a pessimistic stochastic bilevel program in the context of sequential two-player games, where the leader makes a binary here-and-now decision, and the follower responds a continuous wait-and-see decision after observing the leader's action and revelation of uncertainty. Only the information of the mean, covariance, and support is known. We formulate the problem as a distributionally robus…
▽ More
We study a pessimistic stochastic bilevel program in the context of sequential two-player games, where the leader makes a binary here-and-now decision, and the follower responds a continuous wait-and-see decision after observing the leader's action and revelation of uncertainty. Only the information of the mean, covariance, and support is known. We formulate the problem as a distributionally robust (DR) two-stage problem. The pessimistic DR bilevel program is shown to be equivalent to a generic two-stage distributionally robust stochastic (nonlinear) program with both a random objective and random constraints under proper conditions of ambiguity sets. Under continuous distributions, using linear decision rule approaches, we construct upper bounds on the pessimistic DR bilevel program based on (1) 0-1 semidefinite programming (SDP) approximation and (2) an exact 0-1 copositive programming reformulations. When the ambiguity set is restricted to discrete distributions, an exact 0-1 SDP reformulation is developed, and explicit construction of the worst-case distribution is derived. To further improve the computation of the proposed 0-1 SDPs, a cutting-plane framework is developed. Moreover, based on a mixed-integer linear programming approximation, another cutting-plane algorithm is proposed. Extensive numerical studies are conducted to demonstrate the effectiveness of the proposed approaches on a facility location problem.
△ Less
Submitted 7 June, 2022;
originally announced June 2022.
-
Maximum principal ratio of the signless Laplacian of graphs
Authors:
Lele Liu,
Shengming Hu,
Changxiang He
Abstract:
Let $G$ be a connected graph and $Q(G)$ be the signless Laplacian of $G$. The principal ratio $γ(G)$ of $Q(G)$ is the ratio of the maximum and minimum entries of the Perron vector of $Q(G)$. In this paper, we consider the maximum principal ratio $γ(G)$ among all connected graphs of order $n$, and show that for sufficiently large $n$ the extremal graph is a kite graph obtained by identifying an end…
▽ More
Let $G$ be a connected graph and $Q(G)$ be the signless Laplacian of $G$. The principal ratio $γ(G)$ of $Q(G)$ is the ratio of the maximum and minimum entries of the Perron vector of $Q(G)$. In this paper, we consider the maximum principal ratio $γ(G)$ among all connected graphs of order $n$, and show that for sufficiently large $n$ the extremal graph is a kite graph obtained by identifying an end vertex of a path to any vertex of a complete graph.
△ Less
Submitted 10 April, 2022;
originally announced April 2022.
-
Equivariant formality of corank-one isotropy actions and products of rational spheres
Authors:
Jeffrey Carlson,
Chen He
Abstract:
We completely characterize the pairs of connected Lie groups $G > K$ such that $\mathrm{rank}(G) - \mathrm{rank}(K) = 1$ and the left action of $K$ on $G/K$ is equivariantly formal. The analysis requires us to correct and extend an existing partial classification of homogeneous quotients $G/K$ with the rational homotopy type of a product of an odd- and an even-dimensional sphere.
We completely characterize the pairs of connected Lie groups $G > K$ such that $\mathrm{rank}(G) - \mathrm{rank}(K) = 1$ and the left action of $K$ on $G/K$ is equivariantly formal. The analysis requires us to correct and extend an existing partial classification of homogeneous quotients $G/K$ with the rational homotopy type of a product of an odd- and an even-dimensional sphere.
△ Less
Submitted 27 November, 2023; v1 submitted 31 March, 2022;
originally announced April 2022.
-
Some multivariable Rado numbers
Authors:
Gang Yang,
Yaping Mao,
Changxiang He,
Zhao Wang
Abstract:
The Rado number of an equation is a Ramsey-theoretic quantity associated to the equation. Let $\mathcal{E}$ be a linear equation. Denote by $\operatorname{R}_r(\mathcal{E})$ the minimal integer, if it exists, such that any $r$-coloring of $[1,\operatorname{R}_r(\mathcal{E})]$ must admit a monochromatic solution to $\mathcal{E}$. In this paper, we give upper and lower bounds for the Rado number of…
▽ More
The Rado number of an equation is a Ramsey-theoretic quantity associated to the equation. Let $\mathcal{E}$ be a linear equation. Denote by $\operatorname{R}_r(\mathcal{E})$ the minimal integer, if it exists, such that any $r$-coloring of $[1,\operatorname{R}_r(\mathcal{E})]$ must admit a monochromatic solution to $\mathcal{E}$. In this paper, we give upper and lower bounds for the Rado number of $\sum_{i=1}^{m-2}x_i+kx_{m-1}=\ell x_{m}$, and some exact values are also given. Furthermore, we derive some results for the cases that $\ell=m=4$ and $m=5, \ell=k+i \ (1\leq i\leq 5)$. As a generalization, the \emph{$r$-color Rado numbers} for linear equations $\mathcal{E}_1,\mathcal{E}_2,...,\mathcal{E}_r$ is defined as the minimal integer, if it exists, such that any $r$-coloring of $[1,\operatorname{R}_r(\mathcal{E}_1,\mathcal{E}_2,...,\mathcal{E}_r)]$ must admit a monochromatic solution to some $\mathcal{E}_i$, where $1\leq i\leq r$. A lower bound for $\operatorname{R}_r(\mathcal{E}_1,\mathcal{E}_2,...,\mathcal{E}_r)$ and the exact values of $\operatorname{R}_2(x+y=z,\ell x=y)=5k$ and $\operatorname{R}_2(x+y=z, x+a=y)$ was given by Lovász Local Lemma.
△ Less
Submitted 10 March, 2022; v1 submitted 5 March, 2022;
originally announced March 2022.
-
Existence of Large-Data Global Weak Solutions to Kinetic Models of Nonhomogeneous Dilute Polymeric Fluids
Authors:
Chuhui He,
Endre Süli
Abstract:
We prove the existence of large-data global-in-time weak solutions to a general class of coupled bead-spring chain models with finitely extensible nonlinear elastic (FENE) type spring potentials for nonhomogeneous incompressible dilute polymeric fluids in a bounded domain in $\mathbb{R}^d$, $d=2$ or $3$. The class of models under consideration involves the Navier--Stokes system with variable densi…
▽ More
We prove the existence of large-data global-in-time weak solutions to a general class of coupled bead-spring chain models with finitely extensible nonlinear elastic (FENE) type spring potentials for nonhomogeneous incompressible dilute polymeric fluids in a bounded domain in $\mathbb{R}^d$, $d=2$ or $3$. The class of models under consideration involves the Navier--Stokes system with variable density, where the viscosity coefficient depends on both the density and the polymer number density, coupled to a Fokker--Planck equation with a density-dependent drag coefficient. The proof is based on combining a truncation of the probability density function with a two-stage Galerkin approximation and weak compactness and compensated compactness techniques to pass to the limits in the sequence of Galerkin approximations and in the truncation level.
△ Less
Submitted 13 February, 2022;
originally announced February 2022.
-
Internal symmetry of the $L_{\leqslant 3}$ algebra arising from a Lie pair
Authors:
Dadi Ni,
Jiahao Cheng,
Zhuo Chen,
Chen He
Abstract:
A Lie pair is an inclusion $A$ to $L$ of Lie algebroids over the same base manifold. In an earlier work, the third author with Bandiera, Stiénon, and Xu introduced a canonical $L_{\leqslant 3}$ algebra $Γ(\wedge^\bullet A^\vee \otimes L/A)$ whose unary bracket is the Chevalley-Eilenberg differential arising from every Lie pair $(L,A)$. In this note, we prove that to such a Lie pair there is an ass…
▽ More
A Lie pair is an inclusion $A$ to $L$ of Lie algebroids over the same base manifold. In an earlier work, the third author with Bandiera, Stiénon, and Xu introduced a canonical $L_{\leqslant 3}$ algebra $Γ(\wedge^\bullet A^\vee \otimes L/A)$ whose unary bracket is the Chevalley-Eilenberg differential arising from every Lie pair $(L,A)$. In this note, we prove that to such a Lie pair there is an associated Lie algebra action by $\mathrm{Der}(L)$ on the $L_{\leqslant 3}$ algebra $Γ(\wedge^\bullet A^\vee \otimes L/A)$. Here $\mathrm{Der}(L)$ is the space of derivations on the Lie algebroid $L$, or infinitesimal automorphisms of $L$. The said action gives rise to a larger scope of gauge equivalences of Maurer-Cartan elements in $Γ(\wedge^\bullet A^\vee \otimes L/A)$, and for this reason we elect to call the $\mathrm{Der}(L)$-action internal symmetry of $Γ(\wedge^\bullet A^\vee \otimes L/A)$.
△ Less
Submitted 4 July, 2022; v1 submitted 31 January, 2022;
originally announced January 2022.
-
Nontrivial solution for Klein-Gordon equation coupled with Born-Infeld theory with critical growth
Authors:
Chuan-Min He,
Lin Li,
Shang-Jie Chen
Abstract:
In this paper, we study the following system \begin{eqnarray*}
\left\{ \begin{array}{ll}
-Δu + V(x)u-(2ω+φ)φu=λf(u)+|u|^{4}u, \ & \text{in} \ \mathbb{R}^{3},
Δφ+ βΔ_4φ= 4π(ω+φ) u^{2}, \ & \text{in}\ \mathbb{R}^{3},\\ \end{array} \right.
\end{eqnarray*} where $f(u)$ without any growth and Ambrosetti-Rabinowitz conditions. We use cut-off function and Moser iteration to obtain the existence of…
▽ More
In this paper, we study the following system \begin{eqnarray*}
\left\{ \begin{array}{ll}
-Δu + V(x)u-(2ω+φ)φu=λf(u)+|u|^{4}u, \ & \text{in} \ \mathbb{R}^{3},
Δφ+ βΔ_4φ= 4π(ω+φ) u^{2}, \ & \text{in}\ \mathbb{R}^{3},\\ \end{array} \right.
\end{eqnarray*} where $f(u)$ without any growth and Ambrosetti-Rabinowitz conditions. We use cut-off function and Moser iteration to obtain the existence of nontrivial solution. Finally, as a by-product of our approaches, we get the same result for Klein-Gordon-Maxwell system.
△ Less
Submitted 9 December, 2021;
originally announced December 2021.
-
Orbital stability of two-component peakons
Authors:
Cheng He,
Xiaochuan Liu,
Changzheng Qu
Abstract:
We prove that the two-component peakon solutions are orbitally stable in the energy space. The system concerned here is a two-component Novikov system, which is an integrable multicomponent extension of the integrable Novikov equation. We improve the method for the scalar peakons to the two-component case with genuine nonlinear interactions by establishing optimal inequalities for the conserved qu…
▽ More
We prove that the two-component peakon solutions are orbitally stable in the energy space. The system concerned here is a two-component Novikov system, which is an integrable multicomponent extension of the integrable Novikov equation. We improve the method for the scalar peakons to the two-component case with genuine nonlinear interactions by establishing optimal inequalities for the conserved quantities involving the coupled structures. Moreover, we also establish the orbital stability for the train-profiles of these two-component peakons by using the refined analysis based on monotonicity of the local energy and an induction method.
△ Less
Submitted 5 January, 2023; v1 submitted 30 November, 2021;
originally announced November 2021.
-
Towards fast weak adversarial training to solve high dimensional parabolic partial differential equations using XNODE-WAN
Authors:
Paul Valsecchi Oliva,
Yue Wu,
Cuiyu He,
Hao Ni
Abstract:
Due to the curse of dimensionality, solving high dimensional parabolic partial differential equations (PDEs) has been a challenging problem for decades. Recently, a weak adversarial network (WAN) proposed in (Y.Zang et al., 2020) offered a flexible and computationally efficient approach to tackle this problem defined on arbitrary domains by leveraging the weak solution. WAN reformulates the PDE pr…
▽ More
Due to the curse of dimensionality, solving high dimensional parabolic partial differential equations (PDEs) has been a challenging problem for decades. Recently, a weak adversarial network (WAN) proposed in (Y.Zang et al., 2020) offered a flexible and computationally efficient approach to tackle this problem defined on arbitrary domains by leveraging the weak solution. WAN reformulates the PDE problem as a generative adversarial network, where the weak solution (primal network) and the test function (adversarial network) are parameterized by the multi-layer deep neural networks (DNNs). However, it is not yet clear whether DNNs are the most effective model for the parabolic PDE solutions as they do not take into account the fundamentally different roles played by time and spatial variables in the solution. To reinforce the difference, we design a novel so-called XNODE model for the primal network, which is built on the neural ODE (NODE) model with additional spatial dependency to incorporate the a priori information of the PDEs and serve as a universal and effective approximation to the solution. The proposed hybrid method (XNODE-WAN), by integrating the XNODE model within the WAN framework, leads to significant improvement in the performance and efficiency of training. Numerical results show that our method can reduce the training time to a fraction of that of the WAN model.
△ Less
Submitted 14 October, 2021;
originally announced October 2021.
-
Quaternion matrix decomposition and its theoretical implications
Authors:
Chang He,
Bo Jiang,
Xihua Zhu
Abstract:
This paper proposes a novel matrix rank-one decomposition for quaternion Hermitian matrices, which admits a stronger property than the previous results in (sturm2003cones,huang2007complex,ai2011new). The enhanced property can be used to drive some improved results in joint numerical range, $\mathcal{S}$-Procedure and quadratically constrained quadratic programming (QCQP) in the quaternion domain,…
▽ More
This paper proposes a novel matrix rank-one decomposition for quaternion Hermitian matrices, which admits a stronger property than the previous results in (sturm2003cones,huang2007complex,ai2011new). The enhanced property can be used to drive some improved results in joint numerical range, $\mathcal{S}$-Procedure and quadratically constrained quadratic programming (QCQP) in the quaternion domain, demonstrating the capability of our new decomposition technique.
△ Less
Submitted 11 September, 2021;
originally announced September 2021.
-
Revisiting the Maximum Principal Ratio of Graphs
Authors:
Lele Liu,
Changxiang He
Abstract:
Let $G$ be a connected graph, the principal ratio of $G$ is the ratio of the maximum and minimum entries of its Perron eigenvector. In 2007, Cioab\v a and Gregory conjectured that among all connected graphs on $n$ vertices, the kite graph attains the maximum principal ratio. In 2018, Tait and Tobin confirmed the conjecture for sufficientlty large $n$. In this article, we show the conjecture is tru…
▽ More
Let $G$ be a connected graph, the principal ratio of $G$ is the ratio of the maximum and minimum entries of its Perron eigenvector. In 2007, Cioab\v a and Gregory conjectured that among all connected graphs on $n$ vertices, the kite graph attains the maximum principal ratio. In 2018, Tait and Tobin confirmed the conjecture for sufficientlty large $n$. In this article, we show the conjecture is true for all $n\geq 5000$.
△ Less
Submitted 21 June, 2021;
originally announced June 2021.
-
Flux recovery for Cut finite element method and its application in a posteriori error estimation
Authors:
Daniela Capatina,
Cuiyu He
Abstract:
In this article, we aim to recover locally conservative and $H(div)$ conforming fluxes for the linear Cut Finite Element Solution with Nitsche's method for Poisson problems with Dirichlet boundary condition. The computation of the conservative flux in the Raviart-Thomas space is completely local and does not require to solve any mixed problem. The $L^2$-norm of the difference between the numerical…
▽ More
In this article, we aim to recover locally conservative and $H(div)$ conforming fluxes for the linear Cut Finite Element Solution with Nitsche's method for Poisson problems with Dirichlet boundary condition. The computation of the conservative flux in the Raviart-Thomas space is completely local and does not require to solve any mixed problem. The $L^2$-norm of the difference between the numerical flux and the recovered flux can then be used as a posteriori error estimator in the adaptive mesh refinement procedure. Theoretically we are able to prove the global reliability and local efficiency. The theoretical results are verified in the numerical results. Moreover, in the numerical results we also observe optimal convergence rate for the flux error.
△ Less
Submitted 4 June, 2021;
originally announced June 2021.
-
Some $α$-spectral extremal results for some digraphs
Authors:
Haiying Shan,
Feifei Wang,
Changxiang He
Abstract:
In this paper, we characterize the extremal digraphs with the maximal or minimal $α$-spectral radius among some digraph classes such as rose digraphs, generalized theta digraphs and tri-ring digraphs with given size $m$. These digraph classes are denoted by $\mathcal{R}_{m}^k$, $\widetilde{\boldsymbolΘ}_k(m)$ and $\INF(m)$ respectively. The main results about spectral extremal digraph by Guo and L…
▽ More
In this paper, we characterize the extremal digraphs with the maximal or minimal $α$-spectral radius among some digraph classes such as rose digraphs, generalized theta digraphs and tri-ring digraphs with given size $m$. These digraph classes are denoted by $\mathcal{R}_{m}^k$, $\widetilde{\boldsymbolΘ}_k(m)$ and $\INF(m)$ respectively. The main results about spectral extremal digraph by Guo and Liu in \cite{MR2954483} and Li and Wang in \cite{MR3777498} are generalized to $α$-spectral graph theory. As a by-product of our main results, an open problem in \cite{MR3777498} is answered.
Furthermore, we determine the digraphs with the first three minimal $α$-spectral radius among all strongly connected digraphs. Meanwhile, we determine the unique digraph with the fourth minimal $α$-spectral radius among all strongly connected digraphs for $0\le α\le \frac{1}{2}$.
△ Less
Submitted 7 May, 2021;
originally announced May 2021.
-
Homothetic covering of convex hulls of compact convex sets
Authors:
Senlin Wu,
Keke Zhang,
Chan He
Abstract:
Let $K$ be a compact convex set and $m$ be a positive integer. The covering functional of $K$ with respect to $m$ is the smallest $λ\in[0,1]$ such that $K$ can be covered by $m$ translates of $λK$. Estimations of the covering functionals of convex hulls of two or more compact convex sets are presented. It is proved that, if a three-dimensional convex body $K$ is the convex hull of two compact conv…
▽ More
Let $K$ be a compact convex set and $m$ be a positive integer. The covering functional of $K$ with respect to $m$ is the smallest $λ\in[0,1]$ such that $K$ can be covered by $m$ translates of $λK$. Estimations of the covering functionals of convex hulls of two or more compact convex sets are presented. It is proved that, if a three-dimensional convex body $K$ is the convex hull of two compact convex sets having no interior points, then the least number $c(K)$ of smaller homothetic copies of $K$ needed to cover $K$ is not greater than $8$ and $c(K)=8$ if and only if $K$ is a parallelepiped.
△ Less
Submitted 14 November, 2021; v1 submitted 18 April, 2021;
originally announced April 2021.
-
Generalized Opinion Dynamics Model for Social Trust Networks in Opinion Maximization
Authors:
Changxiang He,
Jiayuan Zeng,
Shuting Liu,
Guang Zhang,
Xiaofei Qin,
Xuedian Zhang,
Lele Liu
Abstract:
In this paper, we propose a generalized opinion dynamics model (GODM), which can dynamically compute each person's expressed opinion, to solve the internal opinion maximization problem for social trust networks. In the model, we propose a new, reasonable and interpretable confidence index, which is determined by both person's social status and the evaluation around him. By using the theory of diag…
▽ More
In this paper, we propose a generalized opinion dynamics model (GODM), which can dynamically compute each person's expressed opinion, to solve the internal opinion maximization problem for social trust networks. In the model, we propose a new, reasonable and interpretable confidence index, which is determined by both person's social status and the evaluation around him. By using the theory of diagonally dominant, we obtain the optimal analytic solution of the Nash equilibrium with maximum overall opinion. We design a novel algorithm to maximize the overall with given budget by modifying the internal opinions of people in the social trust network, and prove its optimality both from the algorithm itself and the traditional optimization algorithm-ADMM algorithms with $l_1$-regulations. A series of experiments are conducted, and the experimental results show that our method is superior to the state-of-the-art in four datasets. The average benefit has promoted $67.5\%$, $83.2\%$, $31.5\%$, and $33.7\%$ on four datasets, respectively.
△ Less
Submitted 9 March, 2021;
originally announced March 2021.
-
Optimal Control of Convection-Cooling and Numerical Implementation
Authors:
Cuiyu He,
Weiwei Hu,
Lin Mu
Abstract:
This paper is concerned with the problem of enhancing convection-cooling via active control of the incompressible velocity field, described by a stationary diffusion-convection model. This essentially leads to a bilinear optimal control problem. A rigorous proof of the existence of an optimal control is presented and the first order optimality conditions are derived for solving the control using a…
▽ More
This paper is concerned with the problem of enhancing convection-cooling via active control of the incompressible velocity field, described by a stationary diffusion-convection model. This essentially leads to a bilinear optimal control problem. A rigorous proof of the existence of an optimal control is presented and the first order optimality conditions are derived for solving the control using a variational inequality. Moreover, the second order sufficient conditions are established to characterize the local minimizer. Finally, numerical experiments are conducted utilizing finite elements methods together with nonlinear iterative schemes, to demonstrate and validate the effectiveness of our control design.
△ Less
Submitted 23 March, 2021; v1 submitted 3 September, 2020;
originally announced September 2020.
-
A proof of a conjecture on the distance spectral radius and maximum transmission of graphs
Authors:
Lele Liu,
Haiying Shan,
Changxiang He
Abstract:
Let $G$ be a simple connected graph, and $D(G)$ be the distance matrix of $G$. Suppose that $D_{\max}(G)$ and $λ_1(G)$ are the maximum row sum and the spectral radius of $D(G)$, respectively. In this paper, we give a lower bound for $D_{\max}(G)-λ_1(G)$, and characterize the extremal graphs attaining the bound. As a corollary, we solve a conjecture posed by Liu, Shu and Xue.
Let $G$ be a simple connected graph, and $D(G)$ be the distance matrix of $G$. Suppose that $D_{\max}(G)$ and $λ_1(G)$ are the maximum row sum and the spectral radius of $D(G)$, respectively. In this paper, we give a lower bound for $D_{\max}(G)-λ_1(G)$, and characterize the extremal graphs attaining the bound. As a corollary, we solve a conjecture posed by Liu, Shu and Xue.
△ Less
Submitted 3 September, 2020; v1 submitted 29 August, 2020;
originally announced August 2020.
-
Comparison of Shape Derivatives using CutFEM for Ill-posed Bernoulli Free Boundary Problem
Authors:
Erik Burman,
Cuiyu He,
Mats G. Larson
Abstract:
In this paper we discuss a level set approach for the identification of an unknown boundary in a computational domain. The problem takes the form of a Bernoulli problem where only the Dirichlet datum is known on the boundary that is to be identified, but additional information on the Neumann condition is available on the known part of the boundary. The approach uses a classical constrained optimiz…
▽ More
In this paper we discuss a level set approach for the identification of an unknown boundary in a computational domain. The problem takes the form of a Bernoulli problem where only the Dirichlet datum is known on the boundary that is to be identified, but additional information on the Neumann condition is available on the known part of the boundary. The approach uses a classical constrained optimization problem, where a cost functional is minimized with respect to the unknown boundary, the position of which is defined implicitly by a level set function. To solve the optimization problem a steepest descent algorithm using shape derivatives is applied. In each iteration the cut finite element method is used to obtain high accuracy approximations of the pde-model constraint for a given level set configuration without re-meshing. We consider three different shape derivatives. First the classical one, derived using the continuous optimization problem (optimize then discretize). Then the functional is first discretized using the CutFEM method and the shape derivative is evaluated on the finite element functional (discretize then optimize). Finally we consider a third approach, also using a discretized functional. In this case we do not perturb the domain, but consider a so-called boundary value correction method, where a small correction to the boundary position may be included in the weak boundary condition. Using this correction the shape derivative may be obtained by perturbing a distance parameter in the discrete variational formulation. The theoretical discussion is illustrated with a series of numerical examples showing that all three approaches produce similar result on the proposed Bernoulli problem.
△ Less
Submitted 21 August, 2020;
originally announced August 2020.
-
An a posteriori error estimate of the outer normal derivative using dual weights
Authors:
Silvia Bertoluzza,
Erik Burman,
Cuiyu He
Abstract:
We derive a residual based a-posteriori error estimate for the outer normal flux of approximations to {the diffusion problem with variable coefficient}. By analyzing the solution of the adjoint problem, we show that error indicators in the bulk may be defined to be of higher order than those close to the boundary, which lead to more economic meshes. The theory is illustrated with some numerical ex…
▽ More
We derive a residual based a-posteriori error estimate for the outer normal flux of approximations to {the diffusion problem with variable coefficient}. By analyzing the solution of the adjoint problem, we show that error indicators in the bulk may be defined to be of higher order than those close to the boundary, which lead to more economic meshes. The theory is illustrated with some numerical examples.
△ Less
Submitted 25 October, 2021; v1 submitted 17 August, 2020;
originally announced August 2020.
-
On the index of unbalanced signed bicyclic graphs
Authors:
Changxiang He,
Yuying Li,
Haiying Shan,
Wenyan Wang
Abstract:
In this paper, we focus on the index ( largest eigenvalue) of the adjacency matrix of connected signed graphs. We give some general results on the index when the corresponding signed graph is perturbed. As applications, we determine the first five largest index among all unbalanced bicyclic graphs on n >= 36 vertices together with the corresponding extremal signed graphs whose index attain these v…
▽ More
In this paper, we focus on the index ( largest eigenvalue) of the adjacency matrix of connected signed graphs. We give some general results on the index when the corresponding signed graph is perturbed. As applications, we determine the first five largest index among all unbalanced bicyclic graphs on n >= 36 vertices together with the corresponding extremal signed graphs whose index attain these values.
△ Less
Submitted 28 May, 2020;
originally announced May 2020.
-
A Mesh-free Method Using Piecewise Deep Neural Network for Elliptic Interface Problems
Authors:
Cuiyu He,
Xiaozhe Hu,
Lin Mu
Abstract:
In this paper, we propose a novel mesh-free numerical method for solving the elliptic interface problems based on deep learning. We approximate the solution by the neural networks and, since the solution may change dramatically across the interface, we employ different neural networks in different sub-domains. By reformulating the interface problem as a least-squares problem, we discretize the obj…
▽ More
In this paper, we propose a novel mesh-free numerical method for solving the elliptic interface problems based on deep learning. We approximate the solution by the neural networks and, since the solution may change dramatically across the interface, we employ different neural networks in different sub-domains. By reformulating the interface problem as a least-squares problem, we discretize the objective function using mean squared error via sampling and solve the proposed deep least-squares method by standard training algorithms such as stochastic gradient descent. The discretized objective function utilizes only the point-wise information on the sampling points and thus no underlying mesh is required. Doing this circumvents the challenging meshing procedure as well as the numerical integration on the complex interface. To improve the computational efficiency for more challenging problems, we further design an adaptive sampling strategy based on the residual of the least-squares function and propose an adaptive algorithm. Finally, we present several numerical experiments in both 2D and 3D to show the flexibility, effectiveness, and accuracy of the proposed deep least-square method for solving interface problems.
△ Less
Submitted 10 May, 2020;
originally announced May 2020.
-
Generalized Prager-Synge Inequality and Equilibrated Error Estimators for Discontinuous Elements
Authors:
Cuiyu He,
Zhiqiang Cai,
Shun Zhang
Abstract:
The well-known Prager-Synge identity is valid in $H^1(Ω)$ and serves as a foundation for developing equilibrated a posteriori error estimators for continuous elements. In this paper, we introduce a new inequality, that may be regarded as a generalization of the Prager-Synge identity, to be valid for piecewise $H^1(Ω)$ functions for diffusion problems. The inequality is proved to be identity in two…
▽ More
The well-known Prager-Synge identity is valid in $H^1(Ω)$ and serves as a foundation for developing equilibrated a posteriori error estimators for continuous elements. In this paper, we introduce a new inequality, that may be regarded as a generalization of the Prager-Synge identity, to be valid for piecewise $H^1(Ω)$ functions for diffusion problems. The inequality is proved to be identity in two dimensions.
For nonconforming finite element approximation of arbitrary odd order, we propose a fully explicit approach that recovers an equilibrated flux in $H(div; Ω)$ through a local element-wise scheme and that recovers a gradient in $H(curl;Ω)$ through a simple averaging technique over edges. The resulting error estimator is then proved to be globally reliable and locally efficient. Moreover, the reliability and efficiency constants are independent of the jump of the diffusion coefficient regardless of its distribution.
△ Less
Submitted 24 January, 2020;
originally announced January 2020.
-
Robust concurrent topology optimization of structure and its composite material considering uncertainty with imprecise probability
Authors:
Y. Wu,
Eric Li,
Z. C. He,
X. Y. Lin,
H. X. Jiang
Abstract:
This paper studied a robust concurrent topology optimization (RCTO) approach to design the structure and its composite materials simultaneously. For the first time, the material uncertainty with imprecise probability is integrated into the multi-scale concurrent topology optimization (CTO) framework. To describe the imprecise probabilistic uncertainty efficiently, the type I hybrid interval random…
▽ More
This paper studied a robust concurrent topology optimization (RCTO) approach to design the structure and its composite materials simultaneously. For the first time, the material uncertainty with imprecise probability is integrated into the multi-scale concurrent topology optimization (CTO) framework. To describe the imprecise probabilistic uncertainty efficiently, the type I hybrid interval random model is adopted. An improved hybrid perturbation analysis (IHPA) method is formulated to estimate the expectation and stand variance of the objective function in the worst case. Combined with the bi-directional evolutionary structural optimization (BESO) framework, the robust designs of the structure and its composite material are carried out. Several 2D and 3D numerical examples are presented to illustrate the effectiveness of the proposed method. The results show that the proposed method has high efficiency and low precision loss. In addition, the proposed RCTO approach remains efficient in both of linear static and dynamic structures, which shows its extensive adaptability.
△ Less
Submitted 21 October, 2019;
originally announced October 2019.
-
Residual-based a posteriori error estimation for immersed finite element methods
Authors:
Cuiyu He,
Xu Zhang
Abstract:
In this paper we introduce and analyze the residual-based a posteriori error estimation of the partially penalized immersed finite element method for solving elliptic interface problems. The immersed finite element method can be naturally utilized on interface-unfitted meshes. Our a posteriori error estimate is proved to be both reliable and efficient with reliability constant independent of the l…
▽ More
In this paper we introduce and analyze the residual-based a posteriori error estimation of the partially penalized immersed finite element method for solving elliptic interface problems. The immersed finite element method can be naturally utilized on interface-unfitted meshes. Our a posteriori error estimate is proved to be both reliable and efficient with reliability constant independent of the location of the interface. Numerical results indicate that the efficiency constant is independent of the interface location and that the error estimation is robust with respect to the coefficient contrast.
△ Less
Submitted 17 October, 2019;
originally announced October 2019.
-
A Posteriori Error Estimates with Boundary Correction for a Cut Finite Element Method
Authors:
Erik Burman,
Cuiyu He,
Mats G. Larson
Abstract:
In this work we study a residual based a posteriori error estimation for the CutFEM method applied to an elliptic model problem. We consider the problem with non-polygonal boundary and the analysis takes into account the geometry and data approximation on the boundary. The reliability and efficiency are theoretically proved. Moreover, constants are robust with respect to how the domain boundary cu…
▽ More
In this work we study a residual based a posteriori error estimation for the CutFEM method applied to an elliptic model problem. We consider the problem with non-polygonal boundary and the analysis takes into account the geometry and data approximation on the boundary. The reliability and efficiency are theoretically proved. Moreover, constants are robust with respect to how the domain boundary cuts the mesh.
△ Less
Submitted 3 June, 2019;
originally announced June 2019.
-
A new understanding of $ζ(k)$
Authors:
Chenfeng He
Abstract:
In this paper, by introducing a new operation in the vector space of Laurent series, the author derived explicit series for the values of $ζ$-funtion at positive integers, where $ζ$ denotes the Riemann zeta function. The values of $ζ(k),\ k>1$ are largely connected with Bernoulli numbers and binomial numbers. The method in this paper seems new, and the resluts are about divergent series. Using Bor…
▽ More
In this paper, by introducing a new operation in the vector space of Laurent series, the author derived explicit series for the values of $ζ$-funtion at positive integers, where $ζ$ denotes the Riemann zeta function. The values of $ζ(k),\ k>1$ are largely connected with Bernoulli numbers and binomial numbers. The method in this paper seems new, and the resluts are about divergent series. Using Borel summation for these divergent series one can connect $ζ$ function, Bernoulli numbers, and most series representations of Riemann zeta function.
△ Less
Submitted 12 March, 2019; v1 submitted 16 December, 2018;
originally announced December 2018.
-
A new formula for $ζ(s)$
Authors:
Chenfeng He
Abstract:
In this paper, by introducing a new operation in the vector space of analytic functions, the author presents a method for derivating the well-known formulas: $ζ(1-k)=-\frac{B_k}{k}$ and $ζ(1-n,a)=-\frac{B_n(a)}{n}$ , where $ζ$, $ζ(1-n,a)$ denote the Riemann zeta function and the Hurwitz zeta function respectively. $B_k$ is the $k$-th Bernoulli number. Also the author steps further to deduce some i…
▽ More
In this paper, by introducing a new operation in the vector space of analytic functions, the author presents a method for derivating the well-known formulas: $ζ(1-k)=-\frac{B_k}{k}$ and $ζ(1-n,a)=-\frac{B_n(a)}{n}$ , where $ζ$, $ζ(1-n,a)$ denote the Riemann zeta function and the Hurwitz zeta function respectively. $B_k$ is the $k$-th Bernoulli number. Also the author steps further to deduce some identities related to Bernoulli number and Bernoulli polynomial. Moreover, when combining the operation with forward difference, we can show a new formula for Riemann zeta function, i.e. \[ζ(s)=e\sum_{n=0}^{\infty}\sum_{i=0}^{n}(-1)^{n-i}\frac{1}{(n-i)!(1+i)^{s}}.\]
△ Less
Submitted 12 March, 2019; v1 submitted 22 November, 2018;
originally announced November 2018.
-
Primal dual mixed finite element methods for indefinite advection--diffusion equations
Authors:
Erik Burman,
Cuiyu He
Abstract:
We consider primal-dual mixed finite element methods for the advection--diffusion equation. For the primal variable we use standard continuous finite element space and for the flux we use the Raviart-Thomas space. We prove optimal a priori error estimates in the energy- and the $L^2$-norms for the primal variable in the low Peclet regime. In the high Peclet regime we also prove optimal error estim…
▽ More
We consider primal-dual mixed finite element methods for the advection--diffusion equation. For the primal variable we use standard continuous finite element space and for the flux we use the Raviart-Thomas space. We prove optimal a priori error estimates in the energy- and the $L^2$-norms for the primal variable in the low Peclet regime. In the high Peclet regime we also prove optimal error estimates for the primal variable in the $H(div)$ norm for smooth solutions. Numerically we observe that the method eliminates the spurious oscillations close to interior layers that pollute the solution of the standard Galerkin method when the local Peclet number is high. This method, however, does produce spurious solutions when outflow boundary layer presents. In the last section we propose two simple strategies to remove such numerical artefacts caused by the outflow boundary layer and validate them numerically.
△ Less
Submitted 2 November, 2018;
originally announced November 2018.
-
Limiting Measure of Lee--Yang Zeros for the Cayley Tree
Authors:
Ivan Chio,
Caleb He,
Anthony L. Ji,
Roland K. W. Roeder
Abstract:
This paper is devoted to an in-depth study of the limiting measure of Lee--Yang zeroes for the Ising Model on the Cayley Tree. We build on previous works of Müller-Hartmann-Zittartz (1974 and 1977), Barata--Marchetti (1997), and Barata--Goldbaum (2001), to determine the support of the limiting measure, prove that the limiting measure is not absolutely continuous with respect to Lebesgue measure, a…
▽ More
This paper is devoted to an in-depth study of the limiting measure of Lee--Yang zeroes for the Ising Model on the Cayley Tree. We build on previous works of Müller-Hartmann-Zittartz (1974 and 1977), Barata--Marchetti (1997), and Barata--Goldbaum (2001), to determine the support of the limiting measure, prove that the limiting measure is not absolutely continuous with respect to Lebesgue measure, and determine the pointwise dimension of the measure at Lebesgue a.e. point on the unit circle and every temperature. The latter is related to the critical exponents for the phase transitions in the model as one crosses the unit circle at Lebesgue a.e. point, providing a global version of the "phase transition of continuous order" discovered by Müller-Hartmann-Zittartz. The key techniques are from dynamical systems because there is an explicit formula for the Lee-Yang zeros of the finite Cayley Tree of level $n$ in terms of the $n$-th iterate of an expanding Blaschke Product. A subtlety arises because the conjugacies between Blaschke Products at different parameter values are not absolutely continuous.
△ Less
Submitted 7 January, 2019; v1 submitted 1 June, 2018;
originally announced June 2018.
-
Uniqueness of completions and related topics
Authors:
Chan He,
Horst Martini,
Senlin Wu
Abstract:
A bounded subset of a normed linear space is said to be (diametrically) complete if it cannot be enlarged without increasing the diameter. A complete super set of a bounded set $K$ having the same diameter as $K$ is called a completion of $K$. In general, a bounded set may have different completions. We study normed linear spaces having the property that there exists a nontrivial segment with a un…
▽ More
A bounded subset of a normed linear space is said to be (diametrically) complete if it cannot be enlarged without increasing the diameter. A complete super set of a bounded set $K$ having the same diameter as $K$ is called a completion of $K$. In general, a bounded set may have different completions. We study normed linear spaces having the property that there exists a nontrivial segment with a unique completion. It turns out that this property is strictly weaker than the property that each complete set is a ball, and it is strictly stronger than the property that each set of constant width is a ball. Extensions of this property are also discussed.
△ Less
Submitted 25 February, 2018;
originally announced February 2018.
-
The equivariant cohomology ring of a cohomogeneity-one action
Authors:
Jeffrey D. Carlson,
Oliver Goertsches,
Chen He,
Augustin-Liviu Mare
Abstract:
We compute the rational Borel equivariant cohomology ring of a cohomogeneity-one action of a compact Lie group.
We compute the rational Borel equivariant cohomology ring of a cohomogeneity-one action of a compact Lie group.
△ Less
Submitted 26 February, 2019; v1 submitted 6 February, 2018;
originally announced February 2018.
-
The energy change of the complete multipartite graph
Authors:
Hai-Ying Shan,
Chang-Xiang He,
Zhen-Sheng Yu
Abstract:
The energy of a graph is defined as the sum of the absolute values of all eigenvalues of the graph. Akbari et al. \cite{S. Akbari} proved that for a complete multipartite graph $K_{t_1 ,\ldots,t_k}$, if $t_i\geq 2 \ (i=1,\ldots,k)$, then deleting any edge will increase the energy. A natural question is how the energy changes when $\min\{t_1 ,\ldots,t_k\}=1$. In this paper, we will answer this ques…
▽ More
The energy of a graph is defined as the sum of the absolute values of all eigenvalues of the graph. Akbari et al. \cite{S. Akbari} proved that for a complete multipartite graph $K_{t_1 ,\ldots,t_k}$, if $t_i\geq 2 \ (i=1,\ldots,k)$, then deleting any edge will increase the energy. A natural question is how the energy changes when $\min\{t_1 ,\ldots,t_k\}=1$. In this paper, we will answer this question and completely determine how the energy of a complete multipartite graph changes when one edge is removed.
△ Less
Submitted 11 November, 2017;
originally announced November 2017.