Showing 1–50 of 1,227 results for author: Wang, X

A Distance for Geometric Graphs via the Labeled Merge Tree Interleaving Distance

Authors: Erin Wolf Chambers, Elizabeth Munch, Sarah Percival, Xinyi Wang

Abstract: Geometric graphs appear in many real-world data sets, such as road networks, sensor networks, and molecules. We investigate the notion of distance between embedded graphs and present a metric to measure the distance between two geometric graphs via merge trees. In order to preserve as much useful information as possible from the original data, we introduce a way of rotating the sublevel set to obt… ▽ More Geometric graphs appear in many real-world data sets, such as road networks, sensor networks, and molecules. We investigate the notion of distance between embedded graphs and present a metric to measure the distance between two geometric graphs via merge trees. In order to preserve as much useful information as possible from the original data, we introduce a way of rotating the sublevel set to obtain the merge trees via the idea of the directional transform. We represent the merge trees using a surjective multi-labeling scheme and then compute the distance between two representative matrices. We show some theoretically desirable qualities and present two methods of computation: approximation via sampling and exact distance using a kinetic data structure, both in polynomial time. We illustrate its utility by implementing it on two data sets. △ Less

Submitted 12 July, 2024; originally announced July 2024.

Regularity of the $p$-Gauss curvature flow with flat side

Authors: Genggeng Huang, Xu-Jia Wang, Yang Zhou

Abstract: We study the regularity of the $p$-Gauss curvature flow with flat side. In our previous paper(arxiv:2403.12292), we obtained the regularity of the interface, namely the boundary of the flat part. In this paper, we study the regularity of the convex hypersurface near the interface. We study the regularity of the $p$-Gauss curvature flow with flat side. In our previous paper(arxiv:2403.12292), we obtained the regularity of the interface, namely the boundary of the flat part. In this paper, we study the regularity of the convex hypersurface near the interface. △ Less

Submitted 8 July, 2024; originally announced July 2024.

MSC Class: 53E40; 35K96; 35R35; 35B65

Hidden Convexity-Based Distributed Operation of Integrated Electricity-Gas Systems

Authors: Rong-Peng Liu, Yue Song, Junhong Liu, Xiaozhe Wang, Jinpeng Guo, Yunhe Hou

Abstract: We propose a hidden convexity-based method to address distributed optimal energy flow (OEF) problems for transmission-level integrated electricity-gas systems. First, we develop a node-wise decoupling method to de-compose an OEF problem into multiple OEF subproblems. Then, we propose a hidden convexity-based method to equivalently reformulate nonconvex OEF subproblems as semi-definite programs. Th… ▽ More We propose a hidden convexity-based method to address distributed optimal energy flow (OEF) problems for transmission-level integrated electricity-gas systems. First, we develop a node-wise decoupling method to de-compose an OEF problem into multiple OEF subproblems. Then, we propose a hidden convexity-based method to equivalently reformulate nonconvex OEF subproblems as semi-definite programs. This method differs from any ap-proximation and convexification methods that may incur infeasible solutions. Since all OEF subproblems are origi-nally convex or equivalently convexified, we adopt an ADMM to solve the hidden convexity-based distributed OEF problem with convergence analysis. Test results validate the effectiveness of the proposed method, especially in handling a large number of agents. △ Less

Submitted 7 July, 2024; originally announced July 2024.

Comments: 7 pages

ScaleBiO: Scalable Bilevel Optimization for LLM Data Reweighting

Authors: Rui Pan, Jipeng Zhang, Xingyuan Pan, Renjie Pi, Xiaoyu Wang, Tong Zhang

Abstract: Bilevel optimization has shown its utility across various machine learning settings, yet most algorithms in practice require second-order information, making it challenging to scale them up. Only recently, a paradigm of first-order algorithms emerged, capable of effectively addressing bilevel optimization problems. Nevertheless, the practical efficiency of this paradigm remains unverified, particu… ▽ More Bilevel optimization has shown its utility across various machine learning settings, yet most algorithms in practice require second-order information, making it challenging to scale them up. Only recently, a paradigm of first-order algorithms emerged, capable of effectively addressing bilevel optimization problems. Nevertheless, the practical efficiency of this paradigm remains unverified, particularly in the context of large language models (LLMs). This paper introduces the first scalable instantiation of this paradigm called ScaleBiO, focusing on bilevel optimization for large-scale LLM data reweighting. By combining with a recently proposed memory-efficient training technique called LISA, our novel algorithm allows the paradigm to scale to 34-billion-parameter LLMs on eight A40 GPUs, marking the first successful application of bilevel optimization under practical scenarios for large-sized LLMs. Empirically, extensive experiments on data reweighting verify the effectiveness of ScaleBiO for different-scaled models, including GPT-2, LLaMA-3-8B, GPT-NeoX-20B, and Yi-34B, where bilevel optimization succeeds in filtering irrelevant data samples and selecting informative samples. Theoretically, ScaleBiO ensures the optimality of the learned data weights, along with a convergence guarantee matching the conventional first-order bilevel optimization paradigm on smooth and strongly convex objectives. △ Less

Submitted 28 June, 2024; originally announced June 2024.

A note on the stability of two families of two-step schemes

Authors: Xiaoming Wang, Yinqian Yu

Abstract: We investigate the stability of two families of three-level two-step schemes that extend the classical second order BDF (BDF2) and second order Adams-Moulton (AM2) schemes. For a free parameter restricted to an appropriate range that covers the classical case, we show that both the generalized BDF2 and the generalized AM2 schemes are A-stable. We also introduce the concept of uniform-in-time stabi… ▽ More We investigate the stability of two families of three-level two-step schemes that extend the classical second order BDF (BDF2) and second order Adams-Moulton (AM2) schemes. For a free parameter restricted to an appropriate range that covers the classical case, we show that both the generalized BDF2 and the generalized AM2 schemes are A-stable. We also introduce the concept of uniform-in-time stability which characterizes a scheme's ability to inherit the uniform boundedness over all time of the solution of damped and forced equation with the force uniformly bounded in time. We then demonstrate that A-stability and uniform-in-time stability are equivalent for three-level two-step schemes. Next, these two families of schemes are utilized to construct efficient and unconditionally stable IMEX schemes for systems that involve a damping term, a skew symmetric term, and a forcing term. These novel IMEX schemes are shown to be uniform-in-time energy stable in the sense that the norm of any numerical solution is bounded uniformly over all time, provided that the forcing term is uniformly bounded time, the skew symmetric term is dominated by the dissipative term, together with a mild time-step restriction. Numerical experiments verify our theoretical results. They also indicate that the generalized schemes could be more accurate and/or more stable than the classical ones for suitable choice of the parameter. △ Less

Submitted 26 June, 2024; originally announced June 2024.

MSC Class: 65L06; 65L05; 65L20; 65M12

Dependence on parameters of solutions for a generalized poly-Laplacian system on weighted graphs

Authors: Xiaoyu Wang, Junping Xie, Xingyong Zhang, Xin Ou

Abstract: We mainly investigate the continuous dependence on parameters of nontrivial solutions for a generalized poly-Laplacian system on the weighted finite graph $G=(V, E)$. We firstly present an existence result of mountain pass type nontrivial solutions when the nonlinear term $F$ satisfies the super-$(p, q)$ linear growth condition which is a simple generalization of those results in [28]. Then we mai… ▽ More We mainly investigate the continuous dependence on parameters of nontrivial solutions for a generalized poly-Laplacian system on the weighted finite graph $G=(V, E)$. We firstly present an existence result of mountain pass type nontrivial solutions when the nonlinear term $F$ satisfies the super-$(p, q)$ linear growth condition which is a simple generalization of those results in [28]. Then we mainly show that the mountain pass type nontrivial solutions of the poly-Laplacian system are uniformly bounded for parameters and the concrete upper and lower bounds are given, and are continuously dependent on parameters. Similarly, we also present the existence result, the concrete upper and lower bounds, uniqueness, and dependence on parameters for the locally minimum type nontrivial solutions. Subsequently, we present an example on optimal control as an application of our results. Finally, we give a nonexistence result and some results for the corresponding scalar equation. △ Less

Submitted 26 June, 2024; originally announced June 2024.

A projected Euler Method for Random Periodic Solutions of Semi-linear SDEs with non-globally Lipschitz coefficients

Authors: Yujia Guo, Xiaojie Wang, Yue Wu

Abstract: The present work introduces and investigates an explicit time discretization scheme, called the projected Euler method, to numerically approximate random periodic solutions of semi-linear SDEs under non-globally Lipschitz conditions. The existence of the random periodic solution is demonstrated as the limit of the pull-back of the discretized SDE. Without relying on a priori high-order moment boun… ▽ More The present work introduces and investigates an explicit time discretization scheme, called the projected Euler method, to numerically approximate random periodic solutions of semi-linear SDEs under non-globally Lipschitz conditions. The existence of the random periodic solution is demonstrated as the limit of the pull-back of the discretized SDE. Without relying on a priori high-order moment bounds of the numerical approximations, the mean square convergence rate is proved to be order 0.5 for SDEs with multiplicative noise and order 1 for SDEs with additive noise. Numerical examples are also provided to validate our theoretical findings. △ Less

Submitted 27 June, 2024; v1 submitted 23 June, 2024; originally announced June 2024.

Comments: 25 pages,5 figures

MSC Class: 37H99; 60H10; 60H35; 65C30

Towards Dynamic Resource Allocation and Client Scheduling in Hierarchical Federated Learning: A Two-Phase Deep Reinforcement Learning Approach

Authors: Xiaojing Chen, Zhenyuan Li, Wei Ni, Xin Wang, Shunqing Zhang, Yanzan Sun, Shugong Xu, Qingqi Pei

Abstract: Federated learning (FL) is a viable technique to train a shared machine learning model without sharing data. Hierarchical FL (HFL) system has yet to be studied regrading its multiple levels of energy, computation, communication, and client scheduling, especially when it comes to clients relying on energy harvesting to power their operations. This paper presents a new two-phase deep deterministic p… ▽ More Federated learning (FL) is a viable technique to train a shared machine learning model without sharing data. Hierarchical FL (HFL) system has yet to be studied regrading its multiple levels of energy, computation, communication, and client scheduling, especially when it comes to clients relying on energy harvesting to power their operations. This paper presents a new two-phase deep deterministic policy gradient (DDPG) framework, referred to as ``TP-DDPG'', to balance online the learning delay and model accuracy of an FL process in an energy harvesting-powered HFL system. The key idea is that we divide optimization decisions into two groups, and employ DDPG to learn one group in the first phase, while interpreting the other group as part of the environment to provide rewards for training the DDPG in the second phase. Specifically, the DDPG learns the selection of participating clients, and their CPU configurations and the transmission powers. A new straggler-aware client association and bandwidth allocation (SCABA) algorithm efficiently optimizes the other decisions and evaluates the reward for the DDPG. Experiments demonstrate that with substantially reduced number of learnable parameters, the TP-DDPG can quickly converge to effective polices that can shorten the training time of HFL by 39.4% compared to its benchmarks, when the required test accuracy of HFL is 0.9. △ Less

Submitted 21 June, 2024; originally announced June 2024.

On one-step numerical schemes of weak convergence for SDEs with super-linear coefficients

Authors: Yuying Zhao, Xiaojie Wang, Zhongqiang Zhang

Abstract: We consider weak convergence of one-step schemes for solving stochastic differential equations (SDEs) with one-sided Lipschitz conditions. It is known that the super-linear coefficients may lead to a blowup of moments of solutions and their numerical solutions. When solutions to SDEs have all finite moments, weak convergence of numerical schemes has been investigated in [Wang et al (2023), Weak er… ▽ More We consider weak convergence of one-step schemes for solving stochastic differential equations (SDEs) with one-sided Lipschitz conditions. It is known that the super-linear coefficients may lead to a blowup of moments of solutions and their numerical solutions. When solutions to SDEs have all finite moments, weak convergence of numerical schemes has been investigated in [Wang et al (2023), Weak error analysis for strong approximation schemes of SDEs with super-linear coefficients, IMA Journal numerical analysis]. Some modified Euler schemes have been analyzed for weak convergence. In this work, we present a family of explicit schemes of first and second-order weak convergence based on classical schemes for SDEs. We explore the effects of limited moments on these schemes. We provide a systematic but simple way to establish weak convergence orders for schemes based on approximations/modifications of drift and diffusion coefficients. We present several numerical examples of these schemes and show their weak convergence orders. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Koebe uniformization for infinitely connected attracting Fatou domains

Authors: Xiaoguang Wang, Yi Zhong

Abstract: This paper works on the structure of infinitely connected Fatou damains of rational maps in terms of Koebe uniformization. Due to the complicated boundary behavior, the existing uniformization results are failed to apply in general. We proved that if the rational map is geometrically finite, then its infinitely connected attracting Fatou damain is conformally homeomorphic to a circle domain. This paper works on the structure of infinitely connected Fatou damains of rational maps in terms of Koebe uniformization. Due to the complicated boundary behavior, the existing uniformization results are failed to apply in general. We proved that if the rational map is geometrically finite, then its infinitely connected attracting Fatou damain is conformally homeomorphic to a circle domain. △ Less

Submitted 19 June, 2024; originally announced June 2024.

Comments: 13 pages

MSC Class: 30C20(Primary); 30C35(Secondary)

Regular algebras of dimension four associated to coordinate rings of rank-two quadrics

Authors: R. G. Chandler, H. Tran, P. Veerapen, X. Wang

Abstract: In this paper, we classify connected graded quadratic Artin-Schelter regular (AS-regular, henceforth) algebras of global dimension four that have a Hilbert series the same as that of the polynomial ring on four generators and that map onto a twisted homogeneous coordinate ring of a rank-two quadric. A twisted homogeneous coordinate ring is a construction that was defined by Artin, Tate, and Van de… ▽ More In this paper, we classify connected graded quadratic Artin-Schelter regular (AS-regular, henceforth) algebras of global dimension four that have a Hilbert series the same as that of the polynomial ring on four generators and that map onto a twisted homogeneous coordinate ring of a rank-two quadric. A twisted homogeneous coordinate ring is a construction that was defined by Artin, Tate, and Van den Bergh in \cite{ATV1, ATV2, AVdB1990} in the context of the classification of AS-regular algebras of global dimension three. In \cite{SV99,VVr}, the authors classified AS-regular algebras of global dimension four that map onto the twisted homogeneous coordinate ring of a rank-three and a rank-four quadric, respectively. We expand on their work to include the case of coordinate rings of a rank-two quadric. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Comments: Comments are welcome!

MSC Class: 14A22; 16S37; 16S38

Strong convergence rates for long-time approximations of SDEs with non-globally Lipschitz continuous coefficients

Authors: Xiaoming Wu, Xiaojie Wang

Abstract: This paper is concerned with long-time strong approximations of SDEs with non-globally Lipschitz coefficients.Under certain non-globally Lipschitz conditions, a long-time version of fundamental strong convergence theorem is established for general one-step time discretization schemes. With the aid of the fundamental strong convergence theorem, we prove the expected strong convergence rate over inf… ▽ More This paper is concerned with long-time strong approximations of SDEs with non-globally Lipschitz coefficients.Under certain non-globally Lipschitz conditions, a long-time version of fundamental strong convergence theorem is established for general one-step time discretization schemes. With the aid of the fundamental strong convergence theorem, we prove the expected strong convergence rate over infinite time for two types of schemes such as the backward Euler method and the projected Euler method in non-globally Lipschitz settings. Numerical examples are finally reported to confirm our findings. △ Less

Submitted 15 June, 2024; originally announced June 2024.

Comments: 30 pages,4 figures

MSC Class: 60H35; 65C30

Event-Triggered Optimal Tracking Control for Strict-Feedback Nonlinear Systems With Non-Affine Nonlinear Faults

Authors: Ling Wang, Xin Wang, Ziming Wang

Abstract: This article studies the control ideas of the optimal backstepping technique, proposing an event-triggered optimal tracking control scheme for a class of strict-feedback nonlinear systems with non-affine and nonlinear faults. A simplified identifier-critic-actor framework is employed in the reinforcement learning algorithm to achieve optimal control. The identifier estimates the unknown dynamic fu… ▽ More This article studies the control ideas of the optimal backstepping technique, proposing an event-triggered optimal tracking control scheme for a class of strict-feedback nonlinear systems with non-affine and nonlinear faults. A simplified identifier-critic-actor framework is employed in the reinforcement learning algorithm to achieve optimal control. The identifier estimates the unknown dynamic functions, the critic evaluates the system performance, and the actor implements control actions, enabling modeling and control of anonymous systems for achieving optimal control performance. In this paper, a simplified reinforcement learning algorithm is designed by deriving update rules from the negative gradient of a simple positive function related to the Hamilton-Jacobi-Bellman equation, and it also releases the stringent persistent excitation condition. Then, a fault-tolerant control method is developed by applying filtered signals for controller design. Additionally, to address communication resource reduction, an event-triggered mechanism is employed for designing the actual controller. Finally, the proposed scheme's feasibility is validated through theoretical analysis and simulation. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Efficient algorithm for the oscillatory matrix functions

Authors: Dongping Li, Xue Wang, Xiuying Zhang

Abstract: This paper introduces an efficient algorithm for computing the general oscillatory matrix functions. These computations are crucial for solving second-order semi-linear initial value problems. The method is exploited using the scaling and restoring technique based on a quadruple angle formula in conjunction with a truncated Taylor series. The choice of the scaling parameter and the degree of the T… ▽ More This paper introduces an efficient algorithm for computing the general oscillatory matrix functions. These computations are crucial for solving second-order semi-linear initial value problems. The method is exploited using the scaling and restoring technique based on a quadruple angle formula in conjunction with a truncated Taylor series. The choice of the scaling parameter and the degree of the Taylor polynomial relies on a forward error analysis. Numerical experiments show that the new algorithm behaves in a stable fashion and performs well in both accuracy and efficiency. △ Less

Submitted 10 June, 2024; originally announced June 2024.

Comments: 12 pages

MSC Class: 65F30; 65F60 ACM Class: G.1.3

Expressive Power of Graph Neural Networks for (Mixed-Integer) Quadratic Programs

Authors: Ziang Chen, Xiaohan Chen, Jialin Liu, Xinshang Wang, Wotao Yin

Abstract: Quadratic programming (QP) is the most widely applied category of problems in nonlinear programming. Many applications require real-time/fast solutions, though not necessarily with high precision. Existing methods either involve matrix decomposition or use the preconditioned conjugate gradient method. For relatively large instances, these methods cannot achieve the real-time requirement unless the… ▽ More Quadratic programming (QP) is the most widely applied category of problems in nonlinear programming. Many applications require real-time/fast solutions, though not necessarily with high precision. Existing methods either involve matrix decomposition or use the preconditioned conjugate gradient method. For relatively large instances, these methods cannot achieve the real-time requirement unless there is an effective precondition. Recently, graph neural networks (GNNs) opened new possibilities for QP. Some promising empirical studies of applying GNNs for QP tasks show that GNNs can capture key characteristics of an optimization instance and provide adaptive guidance accordingly to crucial configurations during the solving process, or directly provide an approximate solution. Despite notable empirical observations, theoretical foundations are still lacking. In this work, we investigate the expressive or representative power of GNNs, a crucial aspect of neural network theory, specifically in the context of QP tasks, with both continuous and mixed-integer settings. We prove the existence of message-passing GNNs that can reliably represent key properties of quadratic programs, including feasibility, optimal objective value, and optimal solution. Our theory is validated by numerical results. △ Less

Submitted 9 June, 2024; originally announced June 2024.

Spectral properties of the Kramers-Fokker-Planck operator with a long-range potential

Authors: Xue Ping Wang

Abstract: We study real resonances and embedded eigenvalues of the Kramers--Fokker--Planck operator with a long-range potential. We prove that thresholds are only possible accumulation points of eigenvalues and that the limiting absorption principle holds true for energies outside an exceptional set. We also prove that the eigenfunctions associated with discrete eigenvalues decay exponentially and those ass… ▽ More We study real resonances and embedded eigenvalues of the Kramers--Fokker--Planck operator with a long-range potential. We prove that thresholds are only possible accumulation points of eigenvalues and that the limiting absorption principle holds true for energies outside an exceptional set. We also prove that the eigenfunctions associated with discrete eigenvalues decay exponentially and those associated with embedded non-threshold ones decay polynomially. △ Less

Submitted 7 June, 2024; originally announced June 2024.

MSC Class: 35J10; 35P15; 47A55

Model order reduction for discrete time-delay systems with inhomogeneous initial conditions

Authors: Xiaolong Wang, Kejia Xu

Abstract: We propose two kinds of model order reduction methods for discrete time-delay systems with inhomogeneous initial conditions. The peculiar properties of discrete Walsh functions are directly utilized to compute the Walsh coefficients of the systems, and the projection matrix is defined properly to generate reduced models by taking into account the non-zero initial conditions. It is shown that reduc… ▽ More We propose two kinds of model order reduction methods for discrete time-delay systems with inhomogeneous initial conditions. The peculiar properties of discrete Walsh functions are directly utilized to compute the Walsh coefficients of the systems, and the projection matrix is defined properly to generate reduced models by taking into account the non-zero initial conditions. It is shown that reduced models can preserve some Walsh coefficients of the expansion of the original systems. Further, the superposition principle is exploited to achieve a decomposition of the original systems, and a new definition of Gramians is proposed by combining the individual Gramians of each subsystem. As a result, the balanced truncation method is applied to systems with inhomogeneous initial conditions. We also provide a low-rank approximation to Gramians based on the discrete Laguerre polynomials, which enables an efficient execution of our approach. Numerical examples confirm the feasibility and effectiveness of the proposed methods. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Five-dimensional compatible systems and the Tate conjecture for elliptic surfaces

Authors: Lian Duan, Xiyuan Wang, Ariel Weiss

Abstract: Let $(ρ_λ\colon G_{\mathbb Q}\to \operatorname{GL}_5(\overline{E}_λ))_λ$ be a strictly compatible system of Galois representations such that no Hodge--Tate weight has multiplicity $5$. We show that if $ρ_{λ_0}$ is irreducible for some $λ_0$, then $ρ_λ$ is irreducible for all but finitely many $λ$. More generally, if $(ρ_λ)_λ$ is essentially self-dual, we show that either $ρ_λ$ is irreducible for a… ▽ More Let $(ρ_λ\colon G_{\mathbb Q}\to \operatorname{GL}_5(\overline{E}_λ))_λ$ be a strictly compatible system of Galois representations such that no Hodge--Tate weight has multiplicity $5$. We show that if $ρ_{λ_0}$ is irreducible for some $λ_0$, then $ρ_λ$ is irreducible for all but finitely many $λ$. More generally, if $(ρ_λ)_λ$ is essentially self-dual, we show that either $ρ_λ$ is irreducible for all but finitely many $λ$, or the compatible system $(ρ_λ)_λ$ decomposes as a direct sum of lower-dimensional compatible systems. We apply our results to study the Tate conjecture for elliptic surfaces. For example, if $X_0\colon y^2 + (t+3)xy + y= x^3$, we prove the codimension one $\ell$-adic Tate conjecture for all but finitely many $\ell$, for all but finitely many general, degree $3$, genus $2$ branched multiplicative covers of $X_0$. To prove this result, we classify the elliptic surfaces into six families, and prove, using perverse sheaf theory and a result of Cadoret--Tamagawa, that if one surface in a family satisfies the Tate conjecture, then all but finitely many do. We then verify the Tate conjecture for one representative of each family by making our irreducibility result explicit: for the compatible system arising from the transcendental part of $H^2_{\mathrm{et}}(X_{\overline{\mathbb Q}}, \mathbb{Q}_\ell(1))$ for a representative $X$, we formulate an algorithm that takes as input the characteristic polynomials of Frobenius, and terminates if and only if the compatible system is irreducible. △ Less

Submitted 5 June, 2024; originally announced June 2024.

Comments: 52 pages. Comments welcome!

MSC Class: 11F80; 11F70; 14C25; 14D05; 14J27

Skew Knörrer's periodicity Theorem

Authors: Yang Liu, Yuan Shen, Xin Wang

Abstract: In this paper, we introduce a class of twisted matrix algebras of $M_2(E)$ and twisted direct products of $E\times E$ for an algebra $E$. Let $A$ be a noetherian Koszul Artin-Schelter regular algebra, $z\in A_2$ be a regular central element of $A$ and $B=A_P[y_1,y_2;σ]$ be a graded double Ore extension of $A$. We use the Clifford deformation $C_{A^!}(z)$ of Koszul dual $A^!$ to study the noncommut… ▽ More In this paper, we introduce a class of twisted matrix algebras of $M_2(E)$ and twisted direct products of $E\times E$ for an algebra $E$. Let $A$ be a noetherian Koszul Artin-Schelter regular algebra, $z\in A_2$ be a regular central element of $A$ and $B=A_P[y_1,y_2;σ]$ be a graded double Ore extension of $A$. We use the Clifford deformation $C_{A^!}(z)$ of Koszul dual $A^!$ to study the noncommutative quadric hypersurface $B/(z+y_1^2+y_2^2)$. We prove that the stable category of graded maximal Cohen-Macaulay modules over $B/(z+y_1^2+y_2^2)$ is equivalent to certain bounded derived categories, which involve a twisted matrix algebra of $M_2(C_{A^!}(z))$ or a twisted direct product of $C_{A^!}(z)\times C_{A^!}(z)$ depending on the values of $P$. These results are presented as skew versions of Knörrer's periodicity theorem. Moreover, we show $B/(z+y_1^2+y_2^2)$ may not be a noncommutative graded isolated singularity even if $A/(z)$ is. △ Less

Submitted 5 June, 2024; originally announced June 2024.

Large Time Behavior and Sharp Interface Limit of Compressible Navier-Stokes/Allen-Cahn System for Interacting Shock Waves

Authors: Yazhou Chen, Qiaolin He, Xiaoding Shi, Xiaoping Wang

Abstract: In this paper, we study the large time behavior and sharp interface limit of the Cauchy problem for compressible Navier-Stokes/Allen-Cahn system with interaction shock waves in the same family. This system is an important mathematical model for describing the motion of immiscible two-phase flow. The results show that, if the initial density and velocity are near the superposition of two shock wave… ▽ More In this paper, we study the large time behavior and sharp interface limit of the Cauchy problem for compressible Navier-Stokes/Allen-Cahn system with interaction shock waves in the same family. This system is an important mathematical model for describing the motion of immiscible two-phase flow. The results show that, if the initial density and velocity are near the superposition of two shock waves in the same family, then there exists a unique global solution to the compressible Navier-Stokes/Allen-Cahn system, and this solution asymptotically converges to the superposition of the viscous shock wave and rarefaction wave which moving in opposite directions. Moreover, this global-in-time solution converges to the entropy solution of $p$-system in $L^\infty$-norm as the thickness of the diffusion interface tends to zero. △ Less

Submitted 3 June, 2024; originally announced June 2024.

Comments: 41pages, 2 figures

MSC Class: 35Q35; 35B65; 76N10; 35M10; 35B40; 35C20; 76T30

Level proximal subdifferential, variational convexity, and pointwise Lipschitz smoothness

Authors: Honglin Luo, Xianfu Wang, Ziyuan Wang, Xinmin Yang

Abstract: Level proximal subdifferential was introduced by Rockafellar recently as a tool for studying proximal mappings of possibly nonconvex functions. In this paper we give a systematic study of level proximal subdifferntial, characterize variational convexity of the function by locally firm nonexpansiveness of proximal mappings or locally relative monotonicity of level proximal subdifferential, and inve… ▽ More Level proximal subdifferential was introduced by Rockafellar recently as a tool for studying proximal mappings of possibly nonconvex functions. In this paper we give a systematic study of level proximal subdifferntial, characterize variational convexity of the function by locally firm nonexpansiveness of proximal mappings or locally relative monotonicity of level proximal subdifferential, and investigate pointwise Lipschitz smoothness of the function. Integration and single-valuedness of level proximal subdifferential are also examined. △ Less

Submitted 2 June, 2024; originally announced June 2024.

Comments: 25 pages, comments welcomed

MSC Class: Primary 49J53; 49J52; 47H05; Secondary 90C26; 47H09; 26B25

Removable edges in near-bipartite bricks

Authors: Yipei Zhang, Fuliang Lu, Xiumei Wang, Jinjiang Yuan

Abstract: An edge $e$ of a matching covered graph $G$ is removable if $G-e$ is also matching covered. The notion of removable edge arises in connection with ear decompositions of matching covered graphs introduced by Lovász and Plummer. A nonbipartite matching covered graph $G$ is a brick if it is free of nontrivial tight cuts. Carvalho, Lucchesi, and Murty proved that every brick other than $K_4$ and… ▽ More An edge $e$ of a matching covered graph $G$ is removable if $G-e$ is also matching covered. The notion of removable edge arises in connection with ear decompositions of matching covered graphs introduced by Lovász and Plummer. A nonbipartite matching covered graph $G$ is a brick if it is free of nontrivial tight cuts. Carvalho, Lucchesi, and Murty proved that every brick other than $K_4$ and $\overline{C_6}$ has at least $Δ-2$ removable edges. A brick $G$ is near-bipartite if it has a pair of edges $\{e_1,e_2\}$ such that $G-\{e_1,e_2\}$ is a bipartite matching covered graph. In this paper, we show that in a near-bipartite brick $G$ with at least six vertices, every vertex of $G$, except at most six vertices of degree three contained in two disjoint triangles, is incident with at most two nonremovable edges; consequently, $G$ has at least $\frac{|V(G)|-6}{2}$ removable edges. Moreover, all graphs attaining this lower bound are characterized. △ Less

Submitted 31 May, 2024; originally announced June 2024.

Comments: 23 pages, 1 figure

On a problem of Pavlović involving harmonic quasiconformal mappings

Authors: Zhi-Gang Wang, Xiao-Yuan Wang, Antti Rasila, Jia-Le Qiu

Abstract: We obtain a sharp result on order of certain affine and linear invariant families of harmonic quasiconformal mappings with bounded Schwarzian norm. This problem is motivated by the work of Chuaqui, Hernández and Martín [Math. Ann. 367: 1099--1122, 2017]. Firstly, for $K\ge1$, we construct a harmonic $K$-quasiconformal counterpart of the classical Koebe function and use it to formulate the correspo… ▽ More We obtain a sharp result on order of certain affine and linear invariant families of harmonic quasiconformal mappings with bounded Schwarzian norm. This problem is motivated by the work of Chuaqui, Hernández and Martín [Math. Ann. 367: 1099--1122, 2017]. Firstly, for $K\ge1$, we construct a harmonic $K$-quasiconformal counterpart of the classical Koebe function and use it to formulate the corresponding conjectures. Then we consider Hardy spaces $H^p$ of harmonic quasiconformal mappings by applying results for quasiconformal mappings obtained by Astala and Koskela [Pure Appl. Math. Q. 7: 19--50, 2011]. In particular, we determine the optimal order of the family of harmonic quasiconformal mappings with bounded Schwarzian norm to belong to a harmonic Hardy space. This partially solves an open problem posed by Pavlović in 2014. Finally, we derive pre-Schwarzian and Schwarzian norm estimates of certain harmonic mappings. △ Less

Submitted 14 July, 2024; v1 submitted 30 May, 2024; originally announced May 2024.

Comments: 20 pages, 6 figures. Comments are welcome

MSC Class: 30C55; 30C62; 30H10; 31A05

Decentralized Directed Collaboration for Personalized Federated Learning

Authors: Yingqi Liu, Yifan Shi, Qinglun Li, Baoyuan Wu, Xueqian Wang, Li Shen

Abstract: Personalized Federated Learning (PFL) is proposed to find the greatest personalized models for each client. To avoid the central failure and communication bottleneck in the server-based FL, we concentrate on the Decentralized Personalized Federated Learning (DPFL) that performs distributed model training in a Peer-to-Peer (P2P) manner. Most personalized works in DPFL are based on undirected and sy… ▽ More Personalized Federated Learning (PFL) is proposed to find the greatest personalized models for each client. To avoid the central failure and communication bottleneck in the server-based FL, we concentrate on the Decentralized Personalized Federated Learning (DPFL) that performs distributed model training in a Peer-to-Peer (P2P) manner. Most personalized works in DPFL are based on undirected and symmetric topologies, however, the data, computation and communication resources heterogeneity result in large variances in the personalized models, which lead the undirected aggregation to suboptimal personalized performance and unguaranteed convergence. To address these issues, we propose a directed collaboration DPFL framework by incorporating stochastic gradient push and partial model personalized, called \textbf{D}ecentralized \textbf{Fed}erated \textbf{P}artial \textbf{G}radient \textbf{P}ush (\textbf{DFedPGP}). It personalizes the linear classifier in the modern deep model to customize the local solution and learns a consensus representation in a fully decentralized manner. Clients only share gradients with a subset of neighbors based on the directed and asymmetric topologies, which guarantees flexible choices for resource efficiency and better convergence. Theoretically, we show that the proposed DFedPGP achieves a superior convergence rate of $\mathcal{O}(\frac{1}{\sqrt{T}})$ in the general non-convex setting, and prove the tighter connectivity among clients will speed up the convergence. The proposed method achieves state-of-the-art (SOTA) accuracy in both data and computation heterogeneity scenarios, demonstrating the efficiency of the directed collaboration and partial gradient push. △ Less

Submitted 28 May, 2024; originally announced May 2024.

Comments: CVPR 2024. arXiv admin note: text overlap with arXiv:2305.15157

Modified Jost solutions of Schrödinger operators with locally $H^{-1}$ potentials

Authors: Milivoje Lukić, Xingya Wang

Abstract: We study Jost solutions of Schrödinger operators with potentials which decay with respect to a local $H^{-1}$ Sobolev norm; in particular, we generalize to this setting the results of Christ--Kiselev for potentials between the integrable and square-integrable rates of decay, proving existence of solutions with WKB asymptotic behavior on a large set of positive energies. This applies to new classes… ▽ More We study Jost solutions of Schrödinger operators with potentials which decay with respect to a local $H^{-1}$ Sobolev norm; in particular, we generalize to this setting the results of Christ--Kiselev for potentials between the integrable and square-integrable rates of decay, proving existence of solutions with WKB asymptotic behavior on a large set of positive energies. This applies to new classes of potentials which are not locally integrable, or have better decay properties with respect to the $H^{-1}$ norm due to rapid oscillations. △ Less

Submitted 27 May, 2024; originally announced May 2024.

MSC Class: 34L40 (Primary) 35J10 (Secondary)

Claw-free minimal matching covered graphs

Authors: Yipei Zhang, Xiumei Wang, Jinjiang Yuan, C. T. Ng, T. C. E. Cheng

Abstract: A matching covered graph $G$ is minimal if for each edge $e$ of $G$, $G-e$ is not matching covered. An edge $e$ of a matching covered graph $G$ is removable if $G-e$ is also matching covered. Thus a matching covered graph is minimal if and only if it is free of removable edges. For bipartite graphs, Lovász and Plummer gave a characterization of bipartite minimal matching covered graphs. For bricks… ▽ More A matching covered graph $G$ is minimal if for each edge $e$ of $G$, $G-e$ is not matching covered. An edge $e$ of a matching covered graph $G$ is removable if $G-e$ is also matching covered. Thus a matching covered graph is minimal if and only if it is free of removable edges. For bipartite graphs, Lovász and Plummer gave a characterization of bipartite minimal matching covered graphs. For bricks, Lovász showed that the only bricks that are minimal matching covered are $K_4$ and $\overline{C_6}$. In this paper, we present a complete characterization of minimal matching covered graphs that are claw-free. Moreover, for cubic claw-free matching covered graphs that are not minimal matching covered, we obtain the number of their removable edges (with respect to their bricks), and then prove that they have at least 12 removable edges (the bound is sharp). △ Less

Submitted 27 May, 2024; originally announced May 2024.

MSC Class: 05C70; 05C75

Dual-Delayed Asynchronous SGD for Arbitrarily Heterogeneous Data

Authors: Xiaolu Wang, Yuchang Sun, Hoi-To Wai, Jun Zhang

Abstract: We consider the distributed learning problem with data dispersed across multiple workers under the orchestration of a central server. Asynchronous Stochastic Gradient Descent (SGD) has been widely explored in such a setting to reduce the synchronization overhead associated with parallelization. However, the performance of asynchronous SGD algorithms often depends on a bounded dissimilarity conditi… ▽ More We consider the distributed learning problem with data dispersed across multiple workers under the orchestration of a central server. Asynchronous Stochastic Gradient Descent (SGD) has been widely explored in such a setting to reduce the synchronization overhead associated with parallelization. However, the performance of asynchronous SGD algorithms often depends on a bounded dissimilarity condition among the workers' local data, a condition that can drastically affect their efficiency when the workers' data are highly heterogeneous. To overcome this limitation, we introduce the \textit{dual-delayed asynchronous SGD (DuDe-ASGD)} algorithm designed to neutralize the adverse effects of data heterogeneity. DuDe-ASGD makes full use of stale stochastic gradients from all workers during asynchronous training, leading to two distinct time lags in the model parameters and data samples utilized in the server's iterations. Furthermore, by adopting an incremental aggregation strategy, DuDe-ASGD maintains a per-iteration computational cost that is on par with traditional asynchronous SGD algorithms. Our analysis demonstrates that DuDe-ASGD achieves a near-minimax-optimal convergence rate for smooth nonconvex problems, even when the data across workers are extremely heterogeneous. Numerical experiments indicate that DuDe-ASGD compares favorably with existing asynchronous and synchronous SGD-based algorithms. △ Less

Submitted 27 May, 2024; originally announced May 2024.

Distributed Riemannian Stochastic Gradient Tracking Algorithm on the Stiefel Manifold

Authors: Jishu Zhao, Xi Wang, Jinlong Lei

Abstract: This paper focus on investigating the distributed Riemannian stochastic optimization problem on the Stiefel manifold for multi-agent systems, where all the agents work collaboratively to optimize a function modeled by the average of their expectation-valued local costs. Each agent only processes its own local cost function and communicate with neighboring agents to achieve optimal results while en… ▽ More This paper focus on investigating the distributed Riemannian stochastic optimization problem on the Stiefel manifold for multi-agent systems, where all the agents work collaboratively to optimize a function modeled by the average of their expectation-valued local costs. Each agent only processes its own local cost function and communicate with neighboring agents to achieve optimal results while ensuring consensus. Since the local Riemannian gradient in stochastic regimes cannot be directly calculated, we will estimate the gradient by the average of a variable number of sampled gradient, which however brings about noise to the system. We then propose a distributed Riemannian stochastic optimization algorithm on the Stiefel manifold by combining the variable sample size gradient approximation method with the gradient tracking dynamic. It is worth noticing that the suitably chosen increasing sample size plays an important role in improving the algorithm efficiency, as it reduces the noise variance. In an expectation-valued sense, the iterates of all agents are proved to converge to a stationary point (or neighborhood) with fixed step sizes. We further establish the convergence rate of the iterates for the cases when the sample size is exponentially increasing, polynomial increasing, or a constant, respectively. Finally, numerical experiments are implemented to demonstrate the theoretical results. △ Less

Submitted 27 May, 2024; originally announced May 2024.

Quantum-symmetric equivalence is a graded Morita invariant

Authors: Hongdi Huang, Van C. Nguyen, Padmini Veerapen, Kent B. Vashaw, Xingting Wang

Abstract: We show that if two $m$-homogeneous algebras have Morita equivalent graded module categories, then they are quantum-symmetrically equivalent, that is, there is a monoidal equivalence between the categories of comodules for their associated universal quantum groups (in the sense of Manin) which sends one algebra to the other. As a consequence, any Zhang twist of an $m$-homogeneous algebra is a 2-co… ▽ More We show that if two $m$-homogeneous algebras have Morita equivalent graded module categories, then they are quantum-symmetrically equivalent, that is, there is a monoidal equivalence between the categories of comodules for their associated universal quantum groups (in the sense of Manin) which sends one algebra to the other. As a consequence, any Zhang twist of an $m$-homogeneous algebra is a 2-cocycle twist by some 2-cocycle from its Manin's universal quantum group. △ Less

Submitted 20 May, 2024; originally announced May 2024.

MSC Class: 16T05; 16W50; 17B37

Optimal balanced-norm error estimate of the LDG method for reaction-diffusion problems II: the two-dimensional case with layer-upwind flux

Authors: Yao Cheng, Xuesong Wang, Martin Stynes

Abstract: A singularly perturbed reaction-diffusion problem posed on the unit square in $\mathbb{R}^2$ is solved numerically by a local discontinuous Galerkin (LDG) finite element method. Typical solutions of this class of problem exhibit boundary layers along the sides of the domain; these layers generally cause difficulties for numerical methods. Our LDG method handles the boundary layers by using a Shish… ▽ More A singularly perturbed reaction-diffusion problem posed on the unit square in $\mathbb{R}^2$ is solved numerically by a local discontinuous Galerkin (LDG) finite element method. Typical solutions of this class of problem exhibit boundary layers along the sides of the domain; these layers generally cause difficulties for numerical methods. Our LDG method handles the boundary layers by using a Shishkin mesh and also introducing the new concept of a ``layer-upwind flux" -- a discrete flux whose values are chosen on the fine mesh (which lies inside the boundary layers) in the direction where the layer weakens. On the coarse mesh, one can use a standard central flux. No penalty terms are needed with these fluxes, unlike many other variants of the LDG method. Our choice of discrete flux makes it feasible to derive an optimal-order error analysis in a balanced norm; this norm is stronger than the usual energy norm and is a more appropriate measure for errors in computed solutions for singularly perturbed reaction-diffusion problems. It will be proved that the LDG method is usually convergent of order $O((N^{-1}\ln N)^{k+1})$ in the balanced norm, where $N$ is the number of mesh intervals in each coordinate direction and tensor-product piecewise polynomials of degree~$k$ in each coordinate variable are used in the LDG method. This result is the first of its kind for the LDG method applied to this class of problem and is optimal for convergence on a Shishkin mesh. Its sharpness is confirmed by numerical experiments. △ Less

Submitted 20 May, 2024; originally announced May 2024.

Comments: 46 pages, 14 figures

MSC Class: 65N15; 65N30

A class of new linear, efficient and high-order implicit-explicit methods for the coupled free flow-porous media system based on nonlinear Lions interface condition

Authors: Xinhui Wang, Xu Guo, Xiaoli Li

Abstract: In this paper, we construct and analyze new first- and second-order implicit-explicit (IMEX) schemes for the unsteady Navier-Stokes-Darcy model to describe the coupled free flow-porous media system, which is based on the scalar auxiliary variable (SAV) approach in time and finite element method in space. The constructed schemes are linear, only require solving a sequence of linear differential equ… ▽ More In this paper, we construct and analyze new first- and second-order implicit-explicit (IMEX) schemes for the unsteady Navier-Stokes-Darcy model to describe the coupled free flow-porous media system, which is based on the scalar auxiliary variable (SAV) approach in time and finite element method in space. The constructed schemes are linear, only require solving a sequence of linear differential equations with constant coefficients at each time step, and can decouple the Navier-Stokes and Darcy systems. The unconditional stability of both the first- and second-order IMEX schemes can be derived for the coupled system equipped with the Lions interface condition, where the key point is that we should construct a new trilinear form to balance the fully explicit discretizations of the nonlinear terms in the complex system. We can also establish rigorous error estimates for the velocity and hydraulic head of the first-order scheme without any time step restriction. Numerical examples are presented to validate the proposed schemes. △ Less

Submitted 18 May, 2024; originally announced May 2024.

On the transcendentality condition for Gaussian Gabor frames

Authors: Franz Luef, Johannes Testorf, Xu Wang

Abstract: We give a criterion for higher-dimensional Gaussian Gabor frames, which is a reformulation of one of the main results in a previous article by the first and last authors in more explicit terms. We also show that this density criterion for Gaussian Gabor frames is generic in a certain sense. We give a criterion for higher-dimensional Gaussian Gabor frames, which is a reformulation of one of the main results in a previous article by the first and last authors in more explicit terms. We also show that this density criterion for Gaussian Gabor frames is generic in a certain sense. △ Less

Submitted 29 May, 2024; v1 submitted 14 May, 2024; originally announced May 2024.

Comments: 7 pages. Comments welcome

MSC Class: 42C15(Primary) 32Q99 (Secondary)

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.

Comments: Accepted by ICML 2024

Nonlinear Landau damping and wave operators in sharp Gevrey spaces

Authors: A. D. Ionescu, B. Pausader, X. Wang, K. Widmayer

Abstract: We prove nonlinear Landau damping in optimal weighted Gevrey-3 spaces for solutions of the confined Vlasov-Poisson system on $\T^d\times\R^d$ which are small perturbations of homogeneous Penrose-stable equilibria. We also prove the existence of nonlinear scattering operators associated to the confined Vlasov-Poisson evolution, as well as suitable injectivity properties and Lipschitz estimates (a… ▽ More We prove nonlinear Landau damping in optimal weighted Gevrey-3 spaces for solutions of the confined Vlasov-Poisson system on $\T^d\times\R^d$ which are small perturbations of homogeneous Penrose-stable equilibria. We also prove the existence of nonlinear scattering operators associated to the confined Vlasov-Poisson evolution, as well as suitable injectivity properties and Lipschitz estimates (also in weighted Gevrey-3 spaces) on these operators. Our results give definitive answers to two well-known open problems in the field, both of them stated in the recent review of Bedrossian [4, Section 6]. △ Less

Submitted 7 May, 2024; originally announced May 2024.

Comments: 38 pages

On the Kauffman bracket skein module of $(S^1 \times S^2) \ \# \ (S^1 \times S^2)$

Authors: Rhea Palak Bakshi, Seongjeong Kim, Xiao Wang

Abstract: Determining the structure of the Kauffman bracket skein module of all $3$-manifolds over the ring of Laurent polynomials $\mathbb Z[A^{\pm 1}]$ is a big open problem in skein theory. Very little is known about the skein module of non-prime manifolds over this ring. In this paper, we compute the Kauffman bracket skein module of the $3$-manifold $(S^1 \times S^2) \ \# \ (S^1 \times S^2)$ over the ri… ▽ More Determining the structure of the Kauffman bracket skein module of all $3$-manifolds over the ring of Laurent polynomials $\mathbb Z[A^{\pm 1}]$ is a big open problem in skein theory. Very little is known about the skein module of non-prime manifolds over this ring. In this paper, we compute the Kauffman bracket skein module of the $3$-manifold $(S^1 \times S^2) \ \# \ (S^1 \times S^2)$ over the ring $\mathbb Z[A^{\pm 1}]$. We do this by analysing the submodule of handle sliding relations, for which we provide a suitable basis. Along the way we also compute the Kauffman bracket skein module of $(S^1 \times S^2) \ \# \ (S^1 \times D^2)$. △ Less

Submitted 13 May, 2024; v1 submitted 7 May, 2024; originally announced May 2024.

Comments: 30 pages, 20 figures

MSC Class: 57K31; 57K10

A Symplectic Analysis of Alternating Mirror Descent

Authors: Jonas Katona, Xiuyuan Wang, Andre Wibisono

Abstract: Motivated by understanding the behavior of the Alternating Mirror Descent (AMD) algorithm for bilinear zero-sum games, we study the discretization of continuous-time Hamiltonian flow via the symplectic Euler method. We provide a framework for analysis using results from Hamiltonian dynamics, Lie algebra, and symplectic numerical integrators, with an emphasis on the existence and properties of a co… ▽ More Motivated by understanding the behavior of the Alternating Mirror Descent (AMD) algorithm for bilinear zero-sum games, we study the discretization of continuous-time Hamiltonian flow via the symplectic Euler method. We provide a framework for analysis using results from Hamiltonian dynamics, Lie algebra, and symplectic numerical integrators, with an emphasis on the existence and properties of a conserved quantity, the modified Hamiltonian (MH), for the symplectic Euler method. We compute the MH in closed-form when the original Hamiltonian is a quadratic function, and show that it generally differs from the other conserved quantity known previously in that case. We derive new error bounds on the MH when truncated at orders in the stepsize in terms of the number of iterations, $K$, and use these bounds to show an improved $\mathcal{O}(K^{1/5})$ total regret bound and an $\mathcal{O}(K^{-4/5})$ duality gap of the average iterates for AMD. Finally, we propose a conjecture which, if true, would imply that the total regret for AMD scales as $\mathcal{O}\left(K^{\varepsilon}\right)$ and the duality gap of the average iterates as $\mathcal{O}\left(K^{-1+\varepsilon}\right)$ for any $\varepsilon>0$, and we can take $\varepsilon=0$ upon certain convergence conditions for the MH. △ Less

Submitted 28 May, 2024; v1 submitted 6 May, 2024; originally announced May 2024.

Comments: 94 pages, 3 figures

Convergence analysis of a second order numerical scheme for the Flory-Huggins-Cahn-Hilliard-Navier-Stokes system

Authors: Wenbin Chen, Jianyu Jing, Qianqian Liu, Cheng Wang, Xiaoming Wang

Abstract: We present an optimal rate convergence analysis for a second order accurate in time, fully discrete finite difference scheme for the Cahn-Hilliard-Navier-Stokes (CHNS) system, combined with logarithmic Flory-Huggins energy potential. The numerical scheme has been recently proposed, and the positivity-preserving property of the logarithmic arguments, as well as the total energy stability, have been… ▽ More We present an optimal rate convergence analysis for a second order accurate in time, fully discrete finite difference scheme for the Cahn-Hilliard-Navier-Stokes (CHNS) system, combined with logarithmic Flory-Huggins energy potential. The numerical scheme has been recently proposed, and the positivity-preserving property of the logarithmic arguments, as well as the total energy stability, have been theoretically justified. In this paper, we rigorously prove second order convergence of the proposed numerical scheme, in both time and space. Since the CHNS is a coupled system, the standard $\ell^\infty (0, T; \ell^2) \cap \ell^2 (0, T; H_h^2)$ error estimate could not be easily derived, due to the lack of regularity to control the numerical error associated with the coupled terms. Instead, the $\ell^\infty (0, T; H_h^1) \cap \ell^2 (0, T; H_h^3)$ error analysis for the phase variable and the $\ell^\infty (0, T; \ell^2)$ analysis for the velocity vector, which shares the same regularity as the energy estimate, is more suitable to pass through the nonlinear analysis for the error terms associated with the coupled physical process. Furthermore, the highly nonlinear and singular nature of the logarithmic error terms makes the convergence analysis even more challenging, since a uniform distance between the numerical solution and the singular limit values of is needed for the associated error estimate. Many highly non-standard estimates, such as a higher order asymptotic expansion of the numerical solution (up to the third order accuracy in time and fourth order in space), combined with a rough error estimate (to establish the maximum norm bound for the phase variable), as well as a refined error estimate, have to be carried out to conclude the desired convergence result. △ Less

Submitted 4 May, 2024; originally announced May 2024.

MSC Class: 35K35; 35K55; 49J40; 65M06; 65M12

Inexact Adaptive Cubic Regularization Algorithms on Riemannian Manifolds and Application

Authors: Z. Y. Li, X. M. Wang

Abstract: The adaptive cubic regularization algorithm employing the inexact gradient and Hessian is proposed on general Riemannian manifolds, together with the iteration complexity to get an approximate second-order optimality under certain assumptions on accuracies about the inexact gradient and Hessian. The algorithm extends the inexact adaptive cubic regularization algorithm under true gradient in [Math.… ▽ More The adaptive cubic regularization algorithm employing the inexact gradient and Hessian is proposed on general Riemannian manifolds, together with the iteration complexity to get an approximate second-order optimality under certain assumptions on accuracies about the inexact gradient and Hessian. The algorithm extends the inexact adaptive cubic regularization algorithm under true gradient in [Math. Program., 184(1-2): 35-70, 2020] to more general cases even in Euclidean settings. As an application, the algorithm is applied to solve the joint diagonalization problem on the Stiefel manifold. Numerical experiments illustrate that the algorithm performs better than the inexact trust-region algorithm in [Advances of the neural information processing systems, 31, 2018]. △ Less

Submitted 4 May, 2024; originally announced May 2024.

Comments: 15 pages, 1 table

MSC Class: 53C20(Primary); 53C22(Secondary)

Maximum bound principle and original energy dissipation of arbitrarily high-order rescaled exponential time differencing Runge-Kutta schemes for Allen--Cahn equations

Authors: Chaoyu Quan, Xiaoming Wang, Pinzhong Zheng, Zhi Zhou

Abstract: The energy dissipation law and the maximum bound principle are two critical physical properties of the Allen--Cahn equations. While many existing time-stepping methods are known to preserve the energy dissipation law, most apply to a modified form of energy. In this work, we demonstrate that, when the nonlinear term of the Allen--Cahn equation is Lipschitz continuous, a class of arbitrarily high-o… ▽ More The energy dissipation law and the maximum bound principle are two critical physical properties of the Allen--Cahn equations. While many existing time-stepping methods are known to preserve the energy dissipation law, most apply to a modified form of energy. In this work, we demonstrate that, when the nonlinear term of the Allen--Cahn equation is Lipschitz continuous, a class of arbitrarily high-order exponential time differencing Runge--Kutta (ETDRK) schemes preserve the original energy dissipation property, under a mild step-size constraint. Additionally, we guarantee the Lipschitz condition on the nonlinear term by applying a rescaling post-processing technique, which ensures that the numerical solution unconditionally satisfies the maximum bound principle. Consequently, our proposed schemes maintain both the original energy dissipation law and the maximum bound principle and can achieve arbitrarily high-order accuracy. We also establish an optimal error estimate for the proposed schemes. Some numerical experiments are carried out to verify our theoretical results. △ Less

Submitted 29 April, 2024; originally announced April 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.

Multifractal analysis of the power-2-decaying Gauss-like expansion

Authors: Xue-Jiao Wang

Abstract: Each real number $x\in[0,1]$ admits a unique power-2-decaying Gauss-like expansion (P2GLE for short) as $x=\sum_{i\in\mathbb{N}} 2^{-(d_1(x)+d_2(x)+\cdots+d_i(x))}$, where $d_i(x)\in\mathbb{N}$. For any $x\in(0,1]$, the Khintchine exponent $γ(x)$ is defined by $γ(x):=\lim_{n\to\infty}\frac{1}{n}\sum_{j=1}^nd_j(x)$ if the limit exists. We investigate the sizes of the level sets… ▽ More Each real number $x\in[0,1]$ admits a unique power-2-decaying Gauss-like expansion (P2GLE for short) as $x=\sum_{i\in\mathbb{N}} 2^{-(d_1(x)+d_2(x)+\cdots+d_i(x))}$, where $d_i(x)\in\mathbb{N}$. For any $x\in(0,1]$, the Khintchine exponent $γ(x)$ is defined by $γ(x):=\lim_{n\to\infty}\frac{1}{n}\sum_{j=1}^nd_j(x)$ if the limit exists. We investigate the sizes of the level sets $E(ξ):=\{x\in(0,1]:γ(x)=ξ\}$ for $ξ\geq 1$. Utilizing the Ruelle operator theory, we obtain the Khintchine spectrum $ξ\mapsto\dim_H E(ξ)$, where $\dim_H$ denotes the Hausdorff dimension. We establish the remarkable fact that the Khintchine spectrum has exactly one inflection point, which was never proved for the corresponding spectrum in continued fractions. As a direct consequence, we also obtain the Lyapunov spectrum. Furthermore, we find the Hausdorff dimensions of the level sets $\{x\in(0,1]:\lim_{n\to\infty}\frac{1}{n}\sum_{j=1}^{n}\log(d_j(x))=ξ\}$ and $\{x\in(0,1]:\lim_{n\to\infty}\frac{1}{n}\sum_{j=1}^{n}2^{d_j(x)}=ξ\}$. △ Less

Submitted 26 April, 2024; originally announced April 2024.

High-accurate and efficient numerical algorithms for the self-consistent field theory of liquid-crystalline polymers

Authors: Liwei Tan, Zhijuan He, Xin Wang, Kai Jiang

Abstract: In this paper, we develop and investigate numerical methods for the self-consistent field theory (SCFT) of liquid crystalline polymers. Both the Flory-Huggins interaction potential and the Maier-Saupe orientational interaction are considered, enabling simultaneous exploration of microphase separation and liquid crystalline order in these systems. The main challenge in numerically solving this comp… ▽ More In this paper, we develop and investigate numerical methods for the self-consistent field theory (SCFT) of liquid crystalline polymers. Both the Flory-Huggins interaction potential and the Maier-Saupe orientational interaction are considered, enabling simultaneous exploration of microphase separation and liquid crystalline order in these systems. The main challenge in numerically solving this complex system lies in solving 6-dimensional (3D spatial + 2D orientation + 1D time) partial differential equations and performing nonlinear iterations to optimize the fields. We present ten numerical algorithms for solving high-dimensional PDEs and give an evaluation criterion of choosing the most appropriate PDE solver both from theoretical and numerical performance. Results demonstrate that coupling the fourth-order Runge-Kutta scheme with Fourier pseudo-spectral methods and spherical harmonic transformations is superior compared to other algorithms. We develop two nonlinear iteration schemes, alternating direction iteration method and Anderson mixing method, and design an adaptive technique to enhance the stability of the Anderson mixing method. Moreover, the cascadic multi-level (CML) technique further accelerates the SCFT computation when applied to these iteration methods. Meanwhile, an approach to optimize the computational domain is developed. To demonstrate the power of our approaches, we investigate the self-assembled behaviors of flexible-semiflexible diblock copolymers in 4,5,6-dimensional simulations. Numerical results illustrate the efficiency of our proposed method. △ Less

Submitted 18 April, 2024; originally announced April 2024.

Comments: 25 pages, 13figures

Average energy dissipation rates of explicit exponential Runge-Kutta methods for gradient flow problems

Authors: Hong-lin Liao, Xuping Wang

Abstract: We propose a unified theoretical framework to examine the energy dissipation properties at all stages of explicit exponential Runge-Kutta (EERK) methods for gradient flow problems. The main part of the novel framework is to construct the differential form of EERK method by using the difference coefficients of method and the so-called discrete orthogonal convolution kernels. As the main result, we… ▽ More We propose a unified theoretical framework to examine the energy dissipation properties at all stages of explicit exponential Runge-Kutta (EERK) methods for gradient flow problems. The main part of the novel framework is to construct the differential form of EERK method by using the difference coefficients of method and the so-called discrete orthogonal convolution kernels. As the main result, we prove that an EERK method can preserve the original energy dissipation law unconditionally if the associated differentiation matrix is positive semi-definite. A simple indicator, namely average dissipation rate, is also introduced for these multi-stage methods to evaluate the overall energy dissipation rate of an EERK method such that one can choose proper parameters in some parameterized EERK methods or compare different kinds of EERK methods. Some existing EERK methods in the literature are evaluated from the perspective of preserving the original energy dissipation law and the energy dissipation rate. Some numerical examples are also included to support our theory. △ Less

Submitted 23 April, 2024; originally announced April 2024.

Comments: 35 pages, 44 figures

MSC Class: 35K58; 65L20; 65M06; 65M12

Intuitionistic Quantum Logic Perspective: Static and Dynamic Revision Operators

Authors: Heng Zhou, Yongjun Wang, Baoshan Wang, Jian Yan, Xiaoyang Wang

Abstract: The classical belief revision framework, as proposed by Alchourron, Gardenfors, and Makinson, involves the revision of a theory based on eight postulates. In this paper, we focus on the exploration of a revision theory grounded in quantum mechanics, referred to as the natural revision theory. There are two reasoning modes in quantum systems: static intuitionistic reasoning, which incorporates co… ▽ More The classical belief revision framework, as proposed by Alchourron, Gardenfors, and Makinson, involves the revision of a theory based on eight postulates. In this paper, we focus on the exploration of a revision theory grounded in quantum mechanics, referred to as the natural revision theory. There are two reasoning modes in quantum systems: static intuitionistic reasoning, which incorporates contextuality, and dynamic reasoning, which is achieved through projection measurement. We combine the advantages of two intuitionistic quantum logic frameworks, as proposed by D{ö}ring and Coecke, respectively. Our goal is to establish a truth-value assignment for intuitionistic quantum logic that not only aligns with the inherent characteristics of quantum mechanics but also supports truth-value reasoning. The natural revision theory is then investigated based on this approach. We introduce two types of revision operators that correspond to the two reasoning modes in quantum systems: static and dynamic revision. Furthermore, we highlight the distinctions between these two operators. Shifting away from classical revision paradigms, we consider the revision of consequence relations in intuitionistic quantum logic. We demonstrate how, within the natural revision theory framework, both revision operators collectively influence the consequence relations. Notably, the outcomes of revision process are impacted by the sequence in which these interweaved operators are deployed. △ Less

Submitted 30 May, 2024; v1 submitted 21 April, 2024; originally announced April 2024.

FCNCP: A Coupled Nonnegative CANDECOMP/PARAFAC Decomposition Based on Federated Learning

Authors: Yukai Cai, Hang Liu, Xiulin Wang, Hongjin Li, Ziyi Wang, Chuanshuai Yang, Fengyu Cong

Abstract: In the field of brain science, data sharing across servers is becoming increasingly challenging due to issues such as industry competition, privacy security, and administrative procedure policies and regulations. Therefore, there is an urgent need to develop new methods for data analysis and processing that enable scientific collaboration without data sharing. In view of this, this study proposes… ▽ More In the field of brain science, data sharing across servers is becoming increasingly challenging due to issues such as industry competition, privacy security, and administrative procedure policies and regulations. Therefore, there is an urgent need to develop new methods for data analysis and processing that enable scientific collaboration without data sharing. In view of this, this study proposes to study and develop a series of efficient non-negative coupled tensor decomposition algorithm frameworks based on federated learning called FCNCP for the EEG data arranged on different servers. It combining the good discriminative performance of tensor decomposition in high-dimensional data representation and decomposition, the advantages of coupled tensor decomposition in cross-sample tensor data analysis, and the features of federated learning for joint modelling in distributed servers. The algorithm utilises federation learning to establish coupling constraints for data distributed across different servers. In the experiments, firstly, simulation experiments are carried out using simulated data, and stable and consistent decomposition results are obtained, which verify the effectiveness of the proposed algorithms in this study. Then the FCNCP algorithm was utilised to decompose the fifth-order event-related potential (ERP) tensor data collected by applying proprioceptive stimuli on the left and right hands. It was found that contralateral stimulation induced more symmetrical components in the activation areas of the left and right hemispheres. The conclusions drawn are consistent with the interpretations of related studies in cognitive neuroscience, demonstrating that the method can efficiently process higher-order EEG data and that some key hidden information can be preserved. △ Less

Submitted 18 April, 2024; originally announced April 2024.

Statistical Independence and the Brockwell Transform -- From an Integral Equation Perspective

Authors: Xingzhi Wang

Abstract: Statistical independence is a notion ubiquitous in various fields such as in statistics, probability, number theory and physics. We establish the stability of independence for any pair of random variables by their corresponding Brockwell transforms (Brockwell, 2007) beyond the non-atomic condition that is naturally imposed on their distributions, thereby generalizing the proposition originated by… ▽ More Statistical independence is a notion ubiquitous in various fields such as in statistics, probability, number theory and physics. We establish the stability of independence for any pair of random variables by their corresponding Brockwell transforms (Brockwell, 2007) beyond the non-atomic condition that is naturally imposed on their distributions, thereby generalizing the proposition originated by Cai et al. (2022). A central novelty in our work is to formulate the problem as a possibly new type of mathematical inverse problem, which aims to claim that an integral equation has a unique solution in terms of its stochastic kernel -- not under-determined contrary to the usual cases, and also to design a recursive constructive scheme, combined with the property of the quantile function, that solves the aforementioned integral equation in an iterative manner. △ Less

Submitted 11 April, 2024; originally announced April 2024.

A Calabi-Yau-to-Curve Correspondence for Feynman Integrals

Authors: Hans Jockers, Sören Kotlewski, Pyry Kuusela, Andrew J. McLeod, Sebastian Pögel, Maik Sarve, Xing Wang, Stefan Weinzierl

Abstract: It has long been known that the maximal cut of the equal-mass four-loop banana integral is a period of a family of Calabi-Yau threefolds that depends on the kinematic variable $z=m^2/p^2$. We show that it can also be interpreted as a period of a family of genus-two curves. We do this by introducing a general Calabi-Yau-to-curve correspondence, which in this case locally relates the original period… ▽ More It has long been known that the maximal cut of the equal-mass four-loop banana integral is a period of a family of Calabi-Yau threefolds that depends on the kinematic variable $z=m^2/p^2$. We show that it can also be interpreted as a period of a family of genus-two curves. We do this by introducing a general Calabi-Yau-to-curve correspondence, which in this case locally relates the original period of the family of Calabi-Yau threefolds to a period of a family of genus-two curves that varies holomorphically with the kinematic variable $z$. In addition to working out the concrete details of this correspondence for the equal-mass four-loop banana integral, we outline when we expect a correspondence of this type to hold. △ Less

Submitted 8 April, 2024; originally announced April 2024.

Report number: MITP/24-038, TUM-HEP-1503/24

Global solution of 2D hyperbolic liquid crystal system for small initial data

Authors: Xuecheng Wang

Abstract: We prove the global stability of small perturbation near the the constant equilibrium for the two dimensional simplified Ericksen-Leslie's hyperbolic system for incompressible liquid crystal model, where the direction function of liquid crystal molecules satisfies a wave map equation with an acoustical metric. This improves the almost global existence result by Huang-Jiang-Zhao. As byproducts, we… ▽ More We prove the global stability of small perturbation near the the constant equilibrium for the two dimensional simplified Ericksen-Leslie's hyperbolic system for incompressible liquid crystal model, where the direction function of liquid crystal molecules satisfies a wave map equation with an acoustical metric. This improves the almost global existence result by Huang-Jiang-Zhao. As byproducts, we obtain the sharp (same as the linear solution) decay estimates for the nonlinear velocity and the nonlinear wave part. Moreover the nonlinear wave part scatters to a linear solution as time goes to infinity. The main novelty of this paper is that we uncover a null structure inside the velocity equation on the Fourier side for the nonlinear interaction between nonlinear heat equation and nonlinear wave equation. △ Less

Submitted 27 March, 2024; originally announced March 2024.

Comments: 24pages, comments are welcome!

Unconditionally positivity-preserving approximations of the Ait-Sahalia type model: Explicit Milstein-type schemes

Authors: Yingsong Jiang, Ruishu Liu, Xiaojie Wang, Jinghua Zhuo

Abstract: The present article aims to design and analyze efficient first-order strong schemes for a generalized Aït-Sahalia type model arising in mathematical finance and evolving in a positive domain $(0, \infty)$, which possesses a diffusion term with superlinear growth and a highly nonlinear drift that blows up at the origin. Such a complicated structure of the model unavoidably causes essential difficul… ▽ More The present article aims to design and analyze efficient first-order strong schemes for a generalized Aït-Sahalia type model arising in mathematical finance and evolving in a positive domain $(0, \infty)$, which possesses a diffusion term with superlinear growth and a highly nonlinear drift that blows up at the origin. Such a complicated structure of the model unavoidably causes essential difficulties in the construction and convergence analysis of time discretizations. By incorporating implicitness in the term $α_{-1} x^{-1}$ and a corrective mapping $Φ_h$ in the recursion, we develop a novel class of explicit and unconditionally positivity-preserving (i.e., for any step-size $h>0$) Milstein-type schemes for the underlying model. In both non-critical and general critical cases, we introduce a novel approach to analyze mean-square error bounds of the novel schemes, without relying on a priori high-order moment bounds of the numerical approximations. The expected order-one mean-square convergence is attained for the proposed scheme. The above theoretical guarantee can be used to justify the optimal complexity of the Multilevel Monte Carlo method. Numerical experiments are finally provided to verify the theoretical findings. △ Less

Submitted 12 July, 2024; v1 submitted 25 March, 2024; originally announced March 2024.

Comments: 19 pages, 3 figures

Long time regularity of the $p$-Gauss curvature flow with flat side

Authors: G. Huang, X. -J. Wang, Y. Zhou

Abstract: In this paper, we prove the long time regularity of the interface in the $p$-Gauss curvature flow with flat side in all dimensions for $p>\frac1n$. Here the interface is the boundary of the flat part in the flow. In dimension $2$, this problem was solved in \cite{DL2004} for $p=1$ and in \cite{KimLeeRhee2013} for $p\in(1/2,1)$. We utilize the duality method to transform the Gauss curvature flow to… ▽ More In this paper, we prove the long time regularity of the interface in the $p$-Gauss curvature flow with flat side in all dimensions for $p>\frac1n$. Here the interface is the boundary of the flat part in the flow. In dimension $2$, this problem was solved in \cite{DL2004} for $p=1$ and in \cite{KimLeeRhee2013} for $p\in(1/2,1)$. We utilize the duality method to transform the Gauss curvature flow to a singular parabolic Monge-Ampère equation, and prove the regularity of the interface by studying the asymptotic cone of the parabolic Monge-Ampère equation in the polar coordinates. △ Less

Submitted 18 March, 2024; originally announced March 2024.