-
Structure-Preserving Oscillation-Eliminating Discontinuous Galerkin Schemes for Ideal MHD Equations: Locally Divergence-Free and Positivity-Preserving
Authors:
Mengqing Liu,
Kailiang Wu
Abstract:
Numerically simulating magnetohydrodynamics (MHD) poses notable challenges, including the suppression of spurious oscillations near discontinuities (e.g., shocks) and preservation of essential physical structures (e.g., the divergence-free constraint of magnetic field and the positivity of density and pressure). This paper develops structure-preserving oscillation-eliminating discontinuous Galerki…
▽ More
Numerically simulating magnetohydrodynamics (MHD) poses notable challenges, including the suppression of spurious oscillations near discontinuities (e.g., shocks) and preservation of essential physical structures (e.g., the divergence-free constraint of magnetic field and the positivity of density and pressure). This paper develops structure-preserving oscillation-eliminating discontinuous Galerkin (OEDG) schemes for ideal MHD. The schemes leverage a locally divergence-free (LDF) oscillation-eliminating (OE) procedure to suppress spurious oscillations while retaining the LDF property of magnetic field and many desirable attributes of original DG schemes, such as conservation, local compactness, and optimal convergence rates. The OE procedure is based on the solution operator of a novel damping equation, a linear system of ordinary differential equations that are exactly solvable without any discretization. The OE procedure is performed after each Runge-Kutta stage and does not impact DG spatial discretization, facilitating its easy integration into existing DG codes as an independent module. Moreover, this paper presents a rigorous positivity-preserving (PP) analysis of the LDF OEDG schemes on Cartesian meshes, utilizing the optimal convex decomposition technique and the geometric quasi-linearization (GQL) approach. Efficient PP LDF OEDG schemes are derived by incorporating appropriate discretization of Godunov-Powell source terms into only the discrete equations of cell averages, under a condition achievable through a simple PP limiter. Several one- and two-dimensional MHD tests verify the accuracy, effectiveness, and robustness of the proposed structure-preserving OEDG schemes.
△ Less
Submitted 2 May, 2024; v1 submitted 25 April, 2024;
originally announced April 2024.
-
Improved impedance inversion by deep learning and iterated graph Laplacian
Authors:
Davide Bianchi,
Florian Bossmann,
Wenlong Wang,
Mingming Liu
Abstract:
Deep learning techniques have shown significant potential in many applications through recent years. The achieved results often outperform traditional techniques. However, the quality of a neural network highly depends on the used training data. Noisy, insufficient, or biased training data leads to suboptimal results.
We present a hybrid method that combines deep learning with iterated graph Lap…
▽ More
Deep learning techniques have shown significant potential in many applications through recent years. The achieved results often outperform traditional techniques. However, the quality of a neural network highly depends on the used training data. Noisy, insufficient, or biased training data leads to suboptimal results.
We present a hybrid method that combines deep learning with iterated graph Laplacian and show its application in acoustic impedance inversion which is a routine procedure in seismic explorations. A neural network is used to obtain a first approximation of the underlying acoustic impedance and construct a graph Laplacian matrix from this approximation. Afterwards, we use a Tikhonov-like variational method to solve the impedance inversion problem where the regularizer is based on the constructed graph Laplacian. The obtained solution can be shown to be more accurate and stable with respect to noise than the initial guess obtained by the neural network. This process can be iterated several times, each time constructing a new graph Laplacian matrix from the most recent reconstruction. The method converges after only a few iterations returning a much more accurate reconstruction.
We demonstrate the potential of our method on two different datasets and under various levels of noise. We use two different neural networks that have been introduced in previous works. The experiments show that our approach improves the reconstruction quality in the presence of noise.
△ Less
Submitted 25 April, 2024;
originally announced April 2024.
-
$L^p$-regularity of a geometrically nonlinear flat Cosserat micropolar model in supercritical dimensions
Authors:
Chang-Yu Guo,
Chang-Lin Xiang,
Ming-Lun Liu
Abstract:
In a recent work [Ann. Inst. H. Poincaré C Anal. Non Linéaire 2024], Gastel and Neff introduced an interesting system from a geometrically nonlinear flat cosserat micropolar model and established interior regularity in the critical dimension. Motived by this work, in this article, we establish both interior regularity and sharp $L^p$ regularity for their system in supercritical dimensions.
In a recent work [Ann. Inst. H. Poincaré C Anal. Non Linéaire 2024], Gastel and Neff introduced an interesting system from a geometrically nonlinear flat cosserat micropolar model and established interior regularity in the critical dimension. Motived by this work, in this article, we establish both interior regularity and sharp $L^p$ regularity for their system in supercritical dimensions.
△ Less
Submitted 24 April, 2024;
originally announced April 2024.
-
Word-length curve counting on the once-punctured torus
Authors:
David Fisac,
Mingkun Liu
Abstract:
We classify closed curves on a once-punctured torus with a single self-intersection from a combinatorial perspective. We determine the number of closed curves with given word-length and with zero, one, and arbitrary self-intersections.
We classify closed curves on a once-punctured torus with a single self-intersection from a combinatorial perspective. We determine the number of closed curves with given word-length and with zero, one, and arbitrary self-intersections.
△ Less
Submitted 14 April, 2024;
originally announced April 2024.
-
On the elliptic harmonic mappings and sense-preserving harmonic mappings
Authors:
Ming-Sheng Liu,
Hao XU
Abstract:
In this paper, we first establish two versions of Landau-Bloch type theorem for $(K,K')$-elliptic harmonic mappings with a bounded minimum distortion. Next, we provide several coefficient estimates and a conjecture for $(K,K')$-elliptic harmonic mappings. Then, we establish three new versions of Landau-Bloch type theorem for sense-preserving harmonic mappings. Finally, we establish two sharp versi…
▽ More
In this paper, we first establish two versions of Landau-Bloch type theorem for $(K,K')$-elliptic harmonic mappings with a bounded minimum distortion. Next, we provide several coefficient estimates and a conjecture for $(K,K')$-elliptic harmonic mappings. Then, we establish three new versions of Landau-Bloch type theorem for sense-preserving harmonic mappings. Finally, we establish two sharp versions of Landau-Bloch type theorem for certain harmonic mappings. These results are sharp in some given cases and improve the related results of different authors.
△ Less
Submitted 11 April, 2024;
originally announced April 2024.
-
A Structure-Guided Gauss-Newton Method for Shallow ReLU Neural Network
Authors:
Zhiqiang Cai,
Tong Ding,
Min Liu,
Xinyu Liu,
Jianlin Xia
Abstract:
In this paper, we propose a structure-guided Gauss-Newton (SgGN) method for solving least squares problems using a shallow ReLU neural network. The method effectively takes advantage of both the least squares structure and the neural network structure of the objective function. By categorizing the weights and biases of the hidden and output layers of the network as nonlinear and linear parameters,…
▽ More
In this paper, we propose a structure-guided Gauss-Newton (SgGN) method for solving least squares problems using a shallow ReLU neural network. The method effectively takes advantage of both the least squares structure and the neural network structure of the objective function. By categorizing the weights and biases of the hidden and output layers of the network as nonlinear and linear parameters, respectively, the method iterates back and forth between the nonlinear and linear parameters. The nonlinear parameters are updated by a damped Gauss-Newton method and the linear ones are updated by a linear solver. Moreover, at the Gauss-Newton step, a special form of the Gauss-Newton matrix is derived for the shallow ReLU neural network and is used for efficient iterations. It is shown that the corresponding mass and Gauss-Newton matrices in the respective linear and nonlinear steps are symmetric and positive definite under reasonable assumptions. Thus, the SgGN method naturally produces an effective search direction without the need of additional techniques like shifting in the Levenberg-Marquardt method to achieve invertibility of the Gauss-Newton matrix. The convergence and accuracy of the method are demonstrated numerically for several challenging function approximation problems, especially those with discontinuities or sharp transition layers that pose significant challenges for commonly used training algorithms in machine learning.
△ Less
Submitted 7 April, 2024;
originally announced April 2024.
-
Global strong solution to the inviscid liquid-gas two-phase flow model in $L^p$ framework
Authors:
Zhigang Wu,
Mengqian Liu,
Juanzi Cai
Abstract:
This paper is dedicated to the study of the inviscid liquid-gas two-phase flow model in $\mathbb{R}^d\ (d\geq1)$. We establish the global existence of strong solutions to this system with small initial data in hybrid Besov spaces based on general $L^p$-norms. Additionally, we obtain the decay estimates of solutions rely on the constructed Lyapunov functional.
This paper is dedicated to the study of the inviscid liquid-gas two-phase flow model in $\mathbb{R}^d\ (d\geq1)$. We establish the global existence of strong solutions to this system with small initial data in hybrid Besov spaces based on general $L^p$-norms. Additionally, we obtain the decay estimates of solutions rely on the constructed Lyapunov functional.
△ Less
Submitted 7 April, 2024;
originally announced April 2024.
-
An Efficient Quasi-Random Sampling for Copulas
Authors:
Sumin Wang,
Chenxian Huang,
Yongdao Zhou,
Min-Qian Liu
Abstract:
This paper examines an efficient method for quasi-random sampling of copulas in Monte Carlo computations. Traditional methods, like conditional distribution methods (CDM), have limitations when dealing with high-dimensional or implicit copulas, which refer to those that cannot be accurately represented by existing parametric copulas. Instead, this paper proposes the use of generative models, such…
▽ More
This paper examines an efficient method for quasi-random sampling of copulas in Monte Carlo computations. Traditional methods, like conditional distribution methods (CDM), have limitations when dealing with high-dimensional or implicit copulas, which refer to those that cannot be accurately represented by existing parametric copulas. Instead, this paper proposes the use of generative models, such as Generative Adversarial Networks (GANs), to generate quasi-random samples for any copula. GANs are a type of implicit generative models used to learn the distribution of complex data, thus facilitating easy sampling. In our study, GANs are employed to learn the mapping from a uniform distribution to copulas. Once this mapping is learned, obtaining quasi-random samples from the copula only requires inputting quasi-random samples from the uniform distribution. This approach offers a more flexible method for any copula. Additionally, we provide theoretical analysis of quasi-Monte Carlo estimators based on quasi-random samples of copulas. Through simulated and practical applications, particularly in the field of risk management, we validate the proposed method and demonstrate its superiority over various existing methods.
△ Less
Submitted 8 March, 2024;
originally announced March 2024.
-
Quandle Coloring Quivers of general Torus links by dihedral quandles
Authors:
Mohamed Elhamdadi,
Brooke Jones,
Minghui Liu
Abstract:
We completely characterize the coloring quivers of general torus links by dihedral quandles by first exhausting all possible numbers of colorings, followed by determining the interconnections between colorings in each case. The quiver is obtained as function of the number of colorings. The quiver always contains complete subgraphs, in particular a complete subgraph corresponding to the trivial col…
▽ More
We completely characterize the coloring quivers of general torus links by dihedral quandles by first exhausting all possible numbers of colorings, followed by determining the interconnections between colorings in each case. The quiver is obtained as function of the number of colorings. The quiver always contains complete subgraphs, in particular a complete subgraph corresponding to the trivial colorings, but the total number of subgraphs in the quiver and the weights of their edges varies depending on the number of colorings.
△ Less
Submitted 7 March, 2024;
originally announced March 2024.
-
Exact decay of the persistence probability in the Airy$_1$ process
Authors:
Patrik L. Ferrari,
Min Liu
Abstract:
We consider the Airy$_1$ process, which is the limit process in KPZ growth models with flat and non-random initial conditions. We study the persistence probability, namely the probability that the process stays below a given threshold $c$ for a time span of length $L$. This is expected to decay as $e^{-κ(c) L}$. We determine an analytic expression for $κ(c)$ for all $c\geq 3/2$ starting with the c…
▽ More
We consider the Airy$_1$ process, which is the limit process in KPZ growth models with flat and non-random initial conditions. We study the persistence probability, namely the probability that the process stays below a given threshold $c$ for a time span of length $L$. This is expected to decay as $e^{-κ(c) L}$. We determine an analytic expression for $κ(c)$ for all $c\geq 3/2$ starting with the continuum statistics formula for the persistence probability. As the formula is analytic only for $c>0$, we determine an analytic continuation of $κ(c)$ and numerically verify the validity for $c<0$ as well.
△ Less
Submitted 16 February, 2024;
originally announced February 2024.
-
Near optimal constructions of frameproof codes
Authors:
Miao Liu,
Zengjiao Ma,
Chong Shangguan
Abstract:
Frameproof codes are a class of secure codes that were originally introduced in the pioneering work of Boneh and Shaw in the context of digital fingerprinting. They can be used to enhance the security and credibility of digital content. Let $M_{c,l}(q)$ denote the largest cardinality of a $q$-ary $c$-frameproof code with length $l$. Based on an intriguing observation that relates $M_{c,l}(q)$ to t…
▽ More
Frameproof codes are a class of secure codes that were originally introduced in the pioneering work of Boneh and Shaw in the context of digital fingerprinting. They can be used to enhance the security and credibility of digital content. Let $M_{c,l}(q)$ denote the largest cardinality of a $q$-ary $c$-frameproof code with length $l$. Based on an intriguing observation that relates $M_{c,l}(q)$ to the renowned Erdős Matching Conjecture in extremal set theory, in 2003, Blackburn posed an open problem on the precise value of the limit $R_{c,l}=\lim_{q\rightarrow\infty}\frac{M_{c,l}(q)}{q^{\lceil l/c \rceil}}$. By combining several ideas from the probabilistic method, we present a lower bound for $M_{c,l}(q)$, which, together with an upper bound of Blackburn, completely determines $R_{c,l}$ for {\it all} fixed $c,l$, and resolves the above open problem in the full generality. We also present an improved upper bound for $M_{c,l}(q)$.
△ Less
Submitted 12 February, 2024;
originally announced February 2024.
-
Bilevel Optimization under Unbounded Smoothness: A New Algorithm and Convergence Analysis
Authors:
Jie Hao,
Xiaochuan Gong,
Mingrui Liu
Abstract:
Bilevel optimization is an important formulation for many machine learning problems. Current bilevel optimization algorithms assume that the gradient of the upper-level function is Lipschitz. However, recent studies reveal that certain neural networks such as recurrent neural networks (RNNs) and long-short-term memory networks (LSTMs) exhibit potential unbounded smoothness, rendering conventional…
▽ More
Bilevel optimization is an important formulation for many machine learning problems. Current bilevel optimization algorithms assume that the gradient of the upper-level function is Lipschitz. However, recent studies reveal that certain neural networks such as recurrent neural networks (RNNs) and long-short-term memory networks (LSTMs) exhibit potential unbounded smoothness, rendering conventional bilevel optimization algorithms unsuitable. In this paper, we design a new bilevel optimization algorithm, namely BO-REP, to address this challenge. This algorithm updates the upper-level variable using normalized momentum and incorporates two novel techniques for updating the lower-level variable: \textit{initialization refinement} and \textit{periodic updates}. Specifically, once the upper-level variable is initialized, a subroutine is invoked to obtain a refined estimate of the corresponding optimal lower-level variable, and the lower-level variable is updated only after every specific period instead of each iteration. When the upper-level problem is nonconvex and unbounded smooth, and the lower-level problem is strongly convex, we prove that our algorithm requires $\widetilde{\mathcal{O}}(1/ε^4)$ iterations to find an $ε$-stationary point in the stochastic setting, where each iteration involves calling a stochastic gradient or Hessian-vector product oracle. Notably, this result matches the state-of-the-art complexity results under the bounded smoothness setting and without mean-squared smoothness of the stochastic gradient, up to logarithmic factors. Our proof relies on novel technical lemmas for the periodically updated lower-level variable, which are of independent interest. Our experiments on hyper-representation learning, hyperparameter optimization, and data hyper-cleaning for text classification tasks demonstrate the effectiveness of our proposed algorithm.
△ Less
Submitted 17 January, 2024;
originally announced January 2024.
-
Viennot shadows and graded module structure in colored permutation groups
Authors:
Moxuan J. Liu
Abstract:
Let $\mathbf{x}_{n \times n}$ be a matrix of $n \times n$ variables, and let $\mathbb{C}[\mathbf{x}_{n \times n}]$ be the polynomial ring on these variables. Let $\mathfrak{S}_{n,r}$ be the group of colored permutations, consisting of $n \times n$ complex matrices with exactly one nonzero entry in each row and column, and each nonzero entry is a $r$-th root of unity. We associate an ideal…
▽ More
Let $\mathbf{x}_{n \times n}$ be a matrix of $n \times n$ variables, and let $\mathbb{C}[\mathbf{x}_{n \times n}]$ be the polynomial ring on these variables. Let $\mathfrak{S}_{n,r}$ be the group of colored permutations, consisting of $n \times n$ complex matrices with exactly one nonzero entry in each row and column, and each nonzero entry is a $r$-th root of unity. We associate an ideal $I_{\mathfrak{S}_{n,r}} \subseteq \mathbb{C}[\mathbf{x}_{n \times n}]$ with the group $\mathfrak{S}_{n,r}$, and use orbit harmonics to give an ideal-theoretic extension of the Viennot shadow line construction to $\mathfrak{S}_{n,r}$. This extension gives a standard monomial basis of $\mathbb{C}[\mathbf{x}_{n \times n}]/I_{\mathfrak{S}_{n,r}}$, and gives an analogous definition of "longest increasing subsequence" to the group $\mathfrak{S}_{n,r}$. We examine the extension of Chen's conjecture to this analogy. We also study the structure of $\mathbb{C}[\mathbf{x}_{n \times n}]/I_{\mathfrak{S}_{n,r}}$ as a graded $\mathfrak{S}_{n,r} \times \mathfrak{S}_{n,r}$ module, which then induces a graded $\mathfrak{S}_{n,r} \times \mathfrak{S}_{n,r}$ module structure on the $\mathbb{C}$-algebra $\mathbb{C}[\mathfrak{S}_{n,r}]$.
△ Less
Submitted 7 February, 2024; v1 submitted 15 January, 2024;
originally announced January 2024.
-
Decomposition of Hypergeometric SLE and Reversibility
Authors:
Mingchang Liu
Abstract:
In this paper, we consider hypergeometric SLE process for $κ\in (4,8)$ and $ν>\fracκ{2}-6$. Though the definition of hypergeometric SLE process is complicated, we show that given its hitting point on a specific boundary, its conditional law can be described by SLE$_κ(\underlineρ)$ process. Based on this observation, by constructing a pair of curves, we derive the reversibility of hypergeometric SL…
▽ More
In this paper, we consider hypergeometric SLE process for $κ\in (4,8)$ and $ν>\fracκ{2}-6$. Though the definition of hypergeometric SLE process is complicated, we show that given its hitting point on a specific boundary, its conditional law can be described by SLE$_κ(\underlineρ)$ process. Based on this observation, by constructing a pair of curves, we derive the reversibility of hypergeometric SLE for $κ\in(4,8)$ and $ν>-2$.
△ Less
Submitted 7 January, 2024;
originally announced January 2024.
-
Approximate generalized Steiner systems and near-optimal constant weight codes
Authors:
Miao Liu,
Chong Shangguan
Abstract:
Constant weight codes (CWCs) and constant composition codes (CCCs) are two important classes of codes that have been studied extensively in both combinatorics and coding theory for nearly sixty years. In this paper we show that for {\it all} fixed odd distances, there exist near-optimal CWCs and CCCs asymptotically achieving the classic Johnson-type upper bounds.
Let $A_q(n,w,d)$ denote the maxi…
▽ More
Constant weight codes (CWCs) and constant composition codes (CCCs) are two important classes of codes that have been studied extensively in both combinatorics and coding theory for nearly sixty years. In this paper we show that for {\it all} fixed odd distances, there exist near-optimal CWCs and CCCs asymptotically achieving the classic Johnson-type upper bounds.
Let $A_q(n,w,d)$ denote the maximum size of $q$-ary CWCs of length $n$ with constant weight $w$ and minimum distance $d$. One of our main results shows that for {\it all} fixed $q,w$ and odd $d$, one has $\lim_{n\rightarrow\infty}\frac{A_q(n,d,w)}{\binom{n}{t}}=\frac{(q-1)^t}{\binom{w}{t}}$, where $t=\frac{2w-d+1}{2}$. This implies the existence of near-optimal generalized Steiner systems originally introduced by Etzion, and can be viewed as a counterpart of a celebrated result of Rödl on the existence of near-optimal Steiner systems. Note that prior to our work, very little is known about $A_q(n,w,d)$ for $q\ge 3$. A similar result is proved for the maximum size of CCCs.
We provide different proofs for our two main results, based on two strengthenings of the well-known Frankl-Rödl-Pippenger theorem on the existence of near-optimal matchings in hypergraphs: the first proof follows by Kahn's linear programming variation of the above theorem, and the second follows by the recent independent work of Delcour-Postle, and Glock-Joos-Kim-Kühn-Lichev on the existence of near-optimal matchings avoiding certain forbidden configurations.
We also present several intriguing open questions for future research.
△ Less
Submitted 20 January, 2024; v1 submitted 1 January, 2024;
originally announced January 2024.
-
Eigenvalues of quantum Gelfand invariants
Authors:
Naihuan Jing,
Ming Liu,
Alexander Molev
Abstract:
We consider the quantum Gelfand invariants which first appeared in a landmark paper by Reshetikhin, Takhtadzhyan and Faddeev (1989). We calculate the eigenvalues of the invariants acting in irreducible highest weight representations of the quantized enveloping algebra for ${\mathfrak {gl}}_n$. The calculation is based on Liouville-type formulas relating two families of central elements in the quan…
▽ More
We consider the quantum Gelfand invariants which first appeared in a landmark paper by Reshetikhin, Takhtadzhyan and Faddeev (1989). We calculate the eigenvalues of the invariants acting in irreducible highest weight representations of the quantized enveloping algebra for ${\mathfrak {gl}}_n$. The calculation is based on Liouville-type formulas relating two families of central elements in the quantum affine algebras of type $A$.
△ Less
Submitted 15 March, 2024; v1 submitted 21 December, 2023;
originally announced December 2023.
-
Random surfaces with large systoles
Authors:
Mingkun Liu,
Bram Petri
Abstract:
We present two constructions, both inspired by ideas from graph theory, of sequences random surfaces of growing area, whose systoles grow logarithmically as a function of their area. This also allows us to prove a new lower bound on the maximal systole of a closed orientable hyperbolic surface of a given genus.
We present two constructions, both inspired by ideas from graph theory, of sequences random surfaces of growing area, whose systoles grow logarithmically as a function of their area. This also allows us to prove a new lower bound on the maximal systole of a closed orientable hyperbolic surface of a given genus.
△ Less
Submitted 1 March, 2024; v1 submitted 18 December, 2023;
originally announced December 2023.
-
Length spectrum of large genus random metric maps
Authors:
Simon Barazer,
Alessandro Giacchetto,
Mingkun Liu
Abstract:
We study the length of short cycles on uniformly random metric maps (also known as ribbon graphs) of large genus using a Teichmüller theory approach. We establish that, as the genus tends to infinity, the length spectrum converges to a Poisson point process with an explicit intensity. This result extends the work of Janson and Louf to the multi-faced case.
We study the length of short cycles on uniformly random metric maps (also known as ribbon graphs) of large genus using a Teichmüller theory approach. We establish that, as the genus tends to infinity, the length spectrum converges to a Poisson point process with an explicit intensity. This result extends the work of Janson and Louf to the multi-faced case.
△ Less
Submitted 16 December, 2023;
originally announced December 2023.
-
Compressible Navier-Stokes equations without heat conduction in Lp-framework
Authors:
Juanzi Cai,
Zhigang Wu,
Mengqian Liu
Abstract:
In this paper, we mainly consider global well-posedness and long time behavior of compressible Navier-Stokes equations without heat conduction in $L^p$-framework. This is a generalization of Peng and Zhai \cite{peng}(SIMA, 55(2023), no.2, 1439-1463), where they obtained the corresponding result in $L^2$-framework. Based on the key observation that we can release the regularity of non-dissipative e…
▽ More
In this paper, we mainly consider global well-posedness and long time behavior of compressible Navier-Stokes equations without heat conduction in $L^p$-framework. This is a generalization of Peng and Zhai \cite{peng}(SIMA, 55(2023), no.2, 1439-1463), where they obtained the corresponding result in $L^2$-framework. Based on the key observation that we can release the regularity of non-dissipative entropy $S$ in high frequency in \cite{peng}, we ultimately achieve the desired $L^p$ estimate in the high frequency via complicated calculations on the nonlinear terms. In addition, we get the $L^p$-decay rate of the solution.
△ Less
Submitted 14 December, 2023;
originally announced December 2023.
-
Ultra hypercyclicity
Authors:
Martin Liu,
David Walmsley
Abstract:
Recently, new topological properties that an operator acting on a topological vector space can have were introduced: strong hypercyclicity and hypermixing. We introduce a new property called ultra hypercyclicity, and we compare it to the other common topological properties, including the new ones. For the family of weighted backward shifts on $c_0$ and $\ell^p$, we characterize ultra hypercyclicit…
▽ More
Recently, new topological properties that an operator acting on a topological vector space can have were introduced: strong hypercyclicity and hypermixing. We introduce a new property called ultra hypercyclicity, and we compare it to the other common topological properties, including the new ones. For the family of weighted backward shifts on $c_0$ and $\ell^p$, we characterize ultra hypercyclicity in terms of the weight sequence, and frame the open question of whether or not there exists a strongly hypercyclic weighted backward shift which is not ultra hypercyclic in terms of properties that the weight sequence has.
△ Less
Submitted 19 December, 2023; v1 submitted 5 December, 2023;
originally announced December 2023.
-
On fiber and base decompositions in the Fukaya category of a symplectic Landau-Ginzburg model
Authors:
Haniya Azam,
Catherine Cannizzo,
Heather Lee,
Chiu-Chu Melissa Liu
Abstract:
We present some tools for computations in the Fukaya category of a symplectic Landau-Ginzburg model. Specifically, we prove that several computations for these fibrations split into base and fiber computations.
We present some tools for computations in the Fukaya category of a symplectic Landau-Ginzburg model. Specifically, we prove that several computations for these fibrations split into base and fiber computations.
△ Less
Submitted 19 December, 2023; v1 submitted 1 December, 2023;
originally announced December 2023.
-
Finite time blow up for systems of wave equations with multiple speeds posed on asymptotically Euclidean manifold
Authors:
Mengyun Liu
Abstract:
In this work, we study the finite time blow-up phenomenon of three types of semilinear wave systems with multiple speeds, posed on asymptotically Euclidean manifolds. We establish the upper bound estimates for the lifespan of solutions when the spatial dimension $n\geq 2$. In particular, for system related to the Glassey conjecture, we obtained the finite time blow up results under a new curve. Th…
▽ More
In this work, we study the finite time blow-up phenomenon of three types of semilinear wave systems with multiple speeds, posed on asymptotically Euclidean manifolds. We establish the upper bound estimates for the lifespan of solutions when the spatial dimension $n\geq 2$. In particular, for system related to the Glassey conjecture, we obtained the finite time blow up results under a new curve. This new curve is sharp at least for $n=3$.
△ Less
Submitted 29 November, 2023;
originally announced November 2023.
-
Low-lying geodesics on the modular surface and necklaces
Authors:
Ara Basmajian,
Mingkun Liu
Abstract:
The m-thick part of the modular surface X is the smallest compact subsurface of X with horocycle boundary containing all the closed geodesics which wind around the cusp at most m times. The m-thick parts form a compact exhaustion of X. We are interested in the geodesics that lie in the m-thick part (so called m low-lying geodesics). We produce a complete asymptotic expansion for the number of m lo…
▽ More
The m-thick part of the modular surface X is the smallest compact subsurface of X with horocycle boundary containing all the closed geodesics which wind around the cusp at most m times. The m-thick parts form a compact exhaustion of X. We are interested in the geodesics that lie in the m-thick part (so called m low-lying geodesics). We produce a complete asymptotic expansion for the number of m low-lying geodesics of length equal to 2n in the modular surface. In particular, we obtain the asymptotic growth rate of the m low-lying geodesics in terms of their word length using the natural generators of the modular group. After establishing a correspondence between this counting problem and the problem of counting necklaces with n beads, we perform a careful singularity analysis on the associated generating function of the sequence.
△ Less
Submitted 27 November, 2023;
originally announced November 2023.
-
The global strong solution to compressible system with fractional viscous term
Authors:
Mengqian Liu,
Zhigang Wu
Abstract:
This paper is dedicated to the study of compressible Navier-Stokes equations in $\mathbb{R}^d\ (d=2,3)$ with a fractional viscous term $(-Δ)^α$ and $α\in(\frac{1}{2},1)$. We establish the global existence of strong solutions to this system with small initial data restricted in Sobolev spaces $H^{s+1}(\mathbb{R}^d)(s>\frac{d}{2})$ by providing several commutators and using harmonic analysis tools (…
▽ More
This paper is dedicated to the study of compressible Navier-Stokes equations in $\mathbb{R}^d\ (d=2,3)$ with a fractional viscous term $(-Δ)^α$ and $α\in(\frac{1}{2},1)$. We establish the global existence of strong solutions to this system with small initial data restricted in Sobolev spaces $H^{s+1}(\mathbb{R}^d)(s>\frac{d}{2})$ by providing several commutators and using harmonic analysis tools (especially the Littlewood-Paley theory). This reduces the requirement of the regularity $H^{3α+1}(\mathbb{R}^3)$ for the global existence in Wang-Zhang [41]. Additionally, we derive $L^2$-decay rate for the highest derivative of the strong solution without decay loss by using a cancellation of a low-medium-frequency quantity.
△ Less
Submitted 20 November, 2023;
originally announced November 2023.
-
Fukaya category on a symplectic manifold with a B-field
Authors:
Haniya Azam,
Catherine Cannizzo,
Heather Lee,
Chiu-Chu Melissa Liu
Abstract:
We describe the formulation of Fukaya categories of symplectic manifolds with $B$-fields. In addition, we give a formula for how the $A_\infty$ structure maps change as we deform an object by a Lagrangian isotopy.
We describe the formulation of Fukaya categories of symplectic manifolds with $B$-fields. In addition, we give a formula for how the $A_\infty$ structure maps change as we deform an object by a Lagrangian isotopy.
△ Less
Submitted 21 November, 2023; v1 submitted 7 November, 2023;
originally announced November 2023.
-
Simultaneous Nonparametric Inference of M-regression under Complex Temporal Dynamics
Authors:
Miaoshiqi Liu,
Zhou Zhou
Abstract:
The paper considers simultaneous nonparametric inference for a wide class of M-regression models with time-varying coefficients. The covariates and errors of the regression model are tackled as a general class of nonstationary time series and are allowed to be cross-dependent. We construct $\sqrt{n}$-consistent inference for the cumulative regression function, whose limiting properties are disclos…
▽ More
The paper considers simultaneous nonparametric inference for a wide class of M-regression models with time-varying coefficients. The covariates and errors of the regression model are tackled as a general class of nonstationary time series and are allowed to be cross-dependent. We construct $\sqrt{n}$-consistent inference for the cumulative regression function, whose limiting properties are disclosed using Bahadur representation and Gaussian approximation theory. A simple and unified self-convolved bootstrap procedure is proposed. With only one tuning parameter, the bootstrap consistently simulates the desired limiting behavior of the M-estimators under complex temporal dynamics, even under the possible presence of breakpoints in time series. Our methodology leads to a unified framework to conduct general classes of Exact Function Tests, Lack-of-fit Tests, and Qualitative Tests for the time-varying coefficients under complex temporal dynamics. These tests enable one to, among many others, conduct variable selection procedures, check for constancy and linearity, as well as verify shape assumptions, including monotonicity and convexity. As applications, our method is utilized to study the time-varying properties of global climate data and Microsoft stock return, respectively.
△ Less
Submitted 26 February, 2024; v1 submitted 18 October, 2023;
originally announced October 2023.
-
Decomposition of global 2-SLE for $κ\in (4,8)$ and an application for critical FK-Ising model
Authors:
Yu Feng,
Mingchang Liu,
Hao Wu
Abstract:
We consider global 2-SLE$_κ$ $(η_1, η_2)$ in a topological rectangle with $κ\in (4,8)$. We derive the law of a random hitting point of the curves and show that, conditional on this random hitting point, the pair of two curves has the same law as Gaussian free field flow lines with proper boundary data. Using a similar idea, we derive the asymptotic of the probability for $η_1\capη_2=\emptyset$. As…
▽ More
We consider global 2-SLE$_κ$ $(η_1, η_2)$ in a topological rectangle with $κ\in (4,8)$. We derive the law of a random hitting point of the curves and show that, conditional on this random hitting point, the pair of two curves has the same law as Gaussian free field flow lines with proper boundary data. Using a similar idea, we derive the asymptotic of the probability for $η_1\capη_2=\emptyset$. As an application, we derive the asymptotic of the probability for the existence of two disjoint open paths in critical FK-Ising model.
△ Less
Submitted 16 October, 2023;
originally announced October 2023.
-
Scalable Statistical Inference in Non-parametric Least Squares
Authors:
Meimei Liu,
Zuofeng Shang,
Yun Yang
Abstract:
Stochastic approximation (SA) is a powerful and scalable computational method for iteratively estimating the solution of optimization problems in the presence of randomness, particularly well-suited for large-scale and streaming data settings. In this work, we propose a theoretical framework for stochastic approximation (SA) applied to non-parametric least squares in reproducing kernel Hilbert spa…
▽ More
Stochastic approximation (SA) is a powerful and scalable computational method for iteratively estimating the solution of optimization problems in the presence of randomness, particularly well-suited for large-scale and streaming data settings. In this work, we propose a theoretical framework for stochastic approximation (SA) applied to non-parametric least squares in reproducing kernel Hilbert spaces (RKHS), enabling online statistical inference in non-parametric regression models. We achieve this by constructing asymptotically valid pointwise (and simultaneous) confidence intervals (bands) for local (and global) inference of the nonlinear regression function, via employing an online multiplier bootstrap approach to functional stochastic gradient descent (SGD) algorithm in the RKHS. Our main theoretical contributions consist of a unified framework for characterizing the non-asymptotic behavior of the functional SGD estimator and demonstrating the consistency of the multiplier bootstrap method. The proof techniques involve the development of a higher-order expansion of the functional SGD estimator under the supremum norm metric and the Gaussian approximation of suprema of weighted and non-identically distributed empirical processes. Our theory specifically reveals an interesting relationship between the tuning of step sizes in SGD for estimation and the accuracy of uncertainty quantification.
△ Less
Submitted 1 October, 2023;
originally announced October 2023.
-
Open/closed Correspondence and Extended LG/CY Correspondence for Quintic Threefolds
Authors:
Konstantin Aleshkin,
Chiu-Chu Melissa Liu
Abstract:
We show that Walcher's disk potential for the quintic threefold can be represented as a central charge of a specific Gauged Linear Sigma Model which we call the extended quintic GLSM. This representation provides an open/closed correspondence for the quintic threefold since the central charge is a generating function of closed genus-zero GLSM invariants. We also explain how open Landau-Ginzburg/Ca…
▽ More
We show that Walcher's disk potential for the quintic threefold can be represented as a central charge of a specific Gauged Linear Sigma Model which we call the extended quintic GLSM. This representation provides an open/closed correspondence for the quintic threefold since the central charge is a generating function of closed genus-zero GLSM invariants. We also explain how open Landau-Ginzburg/Calabi-Yau correspondence and open mirror symmetry for the quintic are compatible with wall-crossing and mirror symmetry of the extended GLSM, respectively.
△ Less
Submitted 25 September, 2023;
originally announced September 2023.
-
Minimisation of 2-coverings of genus 2 Jacobians
Authors:
Tom Fisher,
Mengzhen Liu
Abstract:
An important problem in computational arithmetic geometry is to find changes of coordinates to simplify a system of polynomial equations with rational coefficients. This is tackled by a combination of two techniques, called minimisation and reduction. We give an algorithm for minimising certain pairs of quadratic forms, subject to the constraint that the first quadratic form is fixed. This has app…
▽ More
An important problem in computational arithmetic geometry is to find changes of coordinates to simplify a system of polynomial equations with rational coefficients. This is tackled by a combination of two techniques, called minimisation and reduction. We give an algorithm for minimising certain pairs of quadratic forms, subject to the constraint that the first quadratic form is fixed. This has applications to 2-descent on the Jacobian of a genus 2 curve.
△ Less
Submitted 12 September, 2023;
originally announced September 2023.
-
On Privacy Preservation of Electric Vehicle Charging Control via State Obfuscation
Authors:
Xiang Huo,
Mingxi Liu
Abstract:
The electric vehicle (EV) industry is rapidly evolving owing to advancements in smart grid technologies and charging control strategies. While EVs are promising in decarbonizing the transportation system and providing grid services, their widespread adoption has led to notable and erratic load injections that can disrupt the normal operation of power grid. Additionally, the unprotected collection…
▽ More
The electric vehicle (EV) industry is rapidly evolving owing to advancements in smart grid technologies and charging control strategies. While EVs are promising in decarbonizing the transportation system and providing grid services, their widespread adoption has led to notable and erratic load injections that can disrupt the normal operation of power grid. Additionally, the unprotected collection and utilization of personal information during the EV charging process cause prevalent privacy issues. To address the scalability and data confidentiality in large-scale EV charging control, we propose a novel decentralized privacy-preserving EV charging control algorithm via state obfuscation that 1) is scalable w.r.t. the number of EVs and ensures optimal EV charging solutions; 2) achieves privacy preservation in the presence of honest-but-curious adversaries and eavesdroppers; and 3) is applicable to eliminate privacy concerns for general multi-agent optimization problems in large-scale cyber-physical systems. The EV charging control is structured as a constrained optimization problem with coupled objectives and constraints, then solved in a decentralized fashion. Privacy analyses and simulations demonstrate the efficiency and efficacy of the proposed approach.
△ Less
Submitted 31 August, 2023;
originally announced September 2023.
-
Blow-up for semilinear wave equations with damping and potential in high dimensional Schwarzschild spacetime
Authors:
Mengliang Liu,
Mengyun Liu
Abstract:
In this work, we study the blow up results to power-type semilinear wave equation in the high dimensional Schwarzschild spacetime, with damping and potential terms. We can obtain the upper bound estimates of lifespan without the assumption that the support of the initial date should be far away from the black hole.
In this work, we study the blow up results to power-type semilinear wave equation in the high dimensional Schwarzschild spacetime, with damping and potential terms. We can obtain the upper bound estimates of lifespan without the assumption that the support of the initial date should be far away from the black hole.
△ Less
Submitted 3 August, 2023;
originally announced August 2023.
-
Least-Squares Neural Network (LSNN) Method For Linear Advection-Reaction Equation: Non-constant Jumps
Authors:
Zhiqiang Cai,
Junpyo Choi,
Min Liu
Abstract:
The least-squares ReLU neural network (LSNN) method was introduced and studied for solving linear advection-reaction equation with discontinuous solution in \cite{Cai2021linear,cai2023least}. The method is based on an equivalent least-squares formulation and employs ReLU neural network (NN) functions with $\lceil \log_2(d+1)\rceil+1$-layer representations for approximating solutions. In this paper…
▽ More
The least-squares ReLU neural network (LSNN) method was introduced and studied for solving linear advection-reaction equation with discontinuous solution in \cite{Cai2021linear,cai2023least}. The method is based on an equivalent least-squares formulation and employs ReLU neural network (NN) functions with $\lceil \log_2(d+1)\rceil+1$-layer representations for approximating solutions. In this paper, we show theoretically that the method is also capable of approximating non-constant jumps along discontinuous interfaces that are not necessarily straight lines. Numerical results for test problems with various non-constant jumps and interfaces show that the LSNN method with $\lceil \log_2(d+1)\rceil+1$ layers approximates solutions accurately with degrees of freedom less than that of mesh-based methods and without the common Gibbs phenomena along discontinuous interfaces.
△ Less
Submitted 5 February, 2024; v1 submitted 12 June, 2023;
originally announced June 2023.
-
Deep Ritz Method with Adaptive Quadrature for Linear Elasticity
Authors:
Min Liu,
Zhiqiang Cai,
Karthik Ramani
Abstract:
In this paper, we study the deep Ritz method for solving the linear elasticity equation from a numerical analysis perspective. A modified Ritz formulation using the $H^{1/2}(Γ_D)$ norm is introduced and analyzed for linear elasticity equation in order to deal with the (essential) Dirichlet boundary condition. We show that the resulting deep Ritz method provides the best approximation among the set…
▽ More
In this paper, we study the deep Ritz method for solving the linear elasticity equation from a numerical analysis perspective. A modified Ritz formulation using the $H^{1/2}(Γ_D)$ norm is introduced and analyzed for linear elasticity equation in order to deal with the (essential) Dirichlet boundary condition. We show that the resulting deep Ritz method provides the best approximation among the set of deep neural network (DNN) functions with respect to the ``energy'' norm. Furthermore, we demonstrate that the total error of the deep Ritz simulation is bounded by the sum of the network approximation error and the numerical integration error, disregarding the algebraic error. To effectively control the numerical integration error, we propose an adaptive quadrature-based numerical integration technique with a residual-based local error indicator. This approach enables efficient approximation of the modified energy functional. Through numerical experiments involving smooth and singular problems, as well as problems with stress concentration, we validate the effectiveness and efficiency of the proposed deep Ritz method with adaptive quadrature.
△ Less
Submitted 24 May, 2023;
originally announced May 2023.
-
Active Learning in the Predict-then-Optimize Framework: A Margin-Based Approach
Authors:
Mo Liu,
Paul Grigas,
Heyuan Liu,
Zuo-Jun Max Shen
Abstract:
We develop the first active learning method in the predict-then-optimize framework. Specifically, we develop a learning method that sequentially decides whether to request the "labels" of feature samples from an unlabeled data stream, where the labels correspond to the parameters of an optimization model for decision-making. Our active learning method is the first to be directly informed by the de…
▽ More
We develop the first active learning method in the predict-then-optimize framework. Specifically, we develop a learning method that sequentially decides whether to request the "labels" of feature samples from an unlabeled data stream, where the labels correspond to the parameters of an optimization model for decision-making. Our active learning method is the first to be directly informed by the decision error induced by the predicted parameters, which is referred to as the Smart Predict-then-Optimize (SPO) loss. Motivated by the structure of the SPO loss, our algorithm adopts a margin-based criterion utilizing the concept of distance to degeneracy and minimizes a tractable surrogate of the SPO loss on the collected data. In particular, we develop an efficient active learning algorithm with both hard and soft rejection variants, each with theoretical excess risk (i.e., generalization) guarantees. We further derive bounds on the label complexity, which refers to the number of samples whose labels are acquired to achieve a desired small level of SPO risk. Under some natural low-noise conditions, we show that these bounds can be better than the naive supervised learning approach that labels all samples. Furthermore, when using the SPO+ loss function, a specialized surrogate of the SPO loss, we derive a significantly smaller label complexity under separability conditions. We also present numerical evidence showing the practical value of our proposed algorithms in the settings of personalized pricing and the shortest path problem.
△ Less
Submitted 11 May, 2023;
originally announced May 2023.
-
Sequential good lattice point sets for computer experiments
Authors:
Xue-Ru Zhang,
Min-Qian Liu,
Dennis K. J. Lin,
Yong-Dao Zhou
Abstract:
Sequential Latin hypercube designs have recently received great attention for computer experiments. Much of the work has been restricted to invariant spaces. The related systematic construction methods are inflexible while algorithmic methods are ineffective for large designs. For such designs in space contraction, systematic construction methods have not been investigated yet. This paper proposes…
▽ More
Sequential Latin hypercube designs have recently received great attention for computer experiments. Much of the work has been restricted to invariant spaces. The related systematic construction methods are inflexible while algorithmic methods are ineffective for large designs. For such designs in space contraction, systematic construction methods have not been investigated yet. This paper proposes a new method for constructing sequential Latin hypercube designs via good lattice point sets in a variety of experimental spaces. These designs are called sequential good lattice point sets. Moreover, we provide fast and efficient approaches for identifying the (nearly) optimal sequential good lattice point sets under a given criterion. Combining with the linear level permutation technique, we obtain a class of asymptotically optimal sequential Latin hypercube designs in invariant spaces where the $L_1$-distance in each stage is either optimal or asymptotically optimal. Numerical results demonstrate that the sequential good lattice point set has a better space-filling property than the existing sequential Latin hypercube designs in the invariant space. It is also shown that the sequential good lattice point sets have less computational complexity and more adaptability.
△ Less
Submitted 16 May, 2023; v1 submitted 11 May, 2023;
originally announced May 2023.
-
Performative Federated Learning: A Solution to Model-Dependent and Heterogeneous Distribution Shifts
Authors:
Kun Jin,
Tongxin Yin,
Zhongzhu Chen,
Zeyu Sun,
Xueru Zhang,
Yang Liu,
Mingyan Liu
Abstract:
We consider a federated learning (FL) system consisting of multiple clients and a server, where the clients aim to collaboratively learn a common decision model from their distributed data. Unlike the conventional FL framework that assumes the client's data is static, we consider scenarios where the clients' data distributions may be reshaped by the deployed decision model. In this work, we levera…
▽ More
We consider a federated learning (FL) system consisting of multiple clients and a server, where the clients aim to collaboratively learn a common decision model from their distributed data. Unlike the conventional FL framework that assumes the client's data is static, we consider scenarios where the clients' data distributions may be reshaped by the deployed decision model. In this work, we leverage the idea of distribution shift mappings in performative prediction to formalize this model-dependent data distribution shift and propose a performative federated learning framework. We first introduce necessary and sufficient conditions for the existence of a unique performative stable solution and characterize its distance to the performative optimal solution. Then we propose the performative FedAvg algorithm and show that it converges to the performative stable solution at a rate of O(1/T) under both full and partial participation schemes. In particular, we use novel proof techniques and show how the clients' heterogeneity influences the convergence. Numerical results validate our analysis and provide valuable insights into real-world applications.
△ Less
Submitted 8 May, 2023;
originally announced May 2023.
-
Ideals and continuity for quantaloid-enriched categories
Authors:
Min Liu,
Shengwei Han,
Isar Stubbe
Abstract:
We study ideals in, and continuity of, quantaloid-enriched categories (Q-categories for short) as a 'many-valued and many-typed' generalization of domain theory. Abstractly, for any (saturated) class Phi of presheaves, we define and study the Phi-continuity of Q-categories. Concretely, we compute three examples of such saturated classes of presheaves - the class of flat ideals, the class of irredu…
▽ More
We study ideals in, and continuity of, quantaloid-enriched categories (Q-categories for short) as a 'many-valued and many-typed' generalization of domain theory. Abstractly, for any (saturated) class Phi of presheaves, we define and study the Phi-continuity of Q-categories. Concretely, we compute three examples of such saturated classes of presheaves - the class of flat ideals, the class of irreducible ideals and the class of conical ideals - which are proper generalizations of ideals in domain theory.
△ Less
Submitted 27 April, 2023;
originally announced April 2023.
-
On the relationship between shortlex order and $A_α$-spectral radii of graphs with starlike branch tree
Authors:
Haiying Shan,
Muhuo Liu
Abstract:
Let $P(n)$ denote the set of all partitions of $n$, whose elements are nondecreasing sequences of positive integers whose sum is $n$. For ${\bf a}=( n_{1}, n_{2},\ldots, n_{d}) \in P(n)$, let $G({\bf a},v)$ denote the graph obtained from connected graph $G$ appending $d$ paths with lengths $n_{1},n_{2},\ldots,n_{d}$ on vertex $v$ of $G$. We show that the ordering of graphs in…
▽ More
Let $P(n)$ denote the set of all partitions of $n$, whose elements are nondecreasing sequences of positive integers whose sum is $n$. For ${\bf a}=( n_{1}, n_{2},\ldots, n_{d}) \in P(n)$, let $G({\bf a},v)$ denote the graph obtained from connected graph $G$ appending $d$ paths with lengths $n_{1},n_{2},\ldots,n_{d}$ on vertex $v$ of $G$. We show that the ordering of graphs in $G_{n}(v)=\{ G({\bf a},v) \mid {\bf a} \in P(n) \}$ by $A_α$-spectral radii coincides with the shortlex ordering of $P(n)$.
△ Less
Submitted 19 April, 2023; v1 submitted 18 April, 2023;
originally announced April 2023.
-
Bohr-type inequalities for unimodular bounded analytic functions
Authors:
Kaixin Chen,
Ming-Sheng Liu,
Saminathan Ponnusamy
Abstract:
In this paper, we establish several new versions of Bohr-type inequalities for bounded analytic functions in the unit disk by allowing $\varphi=\{\varphi_n(r)\}^{\infty}_{n=0}$ in place of the $\{r^n\}^{\infty}_{n=0}$ in the power series representations of the functions involved with the Bohr sum and thereby introducing a single parameter, which generalize several related results of earlier author…
▽ More
In this paper, we establish several new versions of Bohr-type inequalities for bounded analytic functions in the unit disk by allowing $\varphi=\{\varphi_n(r)\}^{\infty}_{n=0}$ in place of the $\{r^n\}^{\infty}_{n=0}$ in the power series representations of the functions involved with the Bohr sum and thereby introducing a single parameter, which generalize several related results of earlier authors.
△ Less
Submitted 15 February, 2023;
originally announced February 2023.
-
Landau-type theorems for certain bounded bi-analytic functions and biharmonic mappings
Authors:
Ming-Sheng Liu,
Saminathan Ponnusamy
Abstract:
In this paper, we establish three new versions of Landau-type theorems for bounded bi-analytic functions of the form $F(z)=\bar{z}G(z)+H(z)$, where $G$ and $H$ are analytic in the unit disk $|z|<1$ with $G(0)=H(0)=0$ and $H'(0)=1$. In particular, two of them are sharp while the other one either generalizes or improves the corresponding result of Abdulhadi and Hajj. As consequences, several new sha…
▽ More
In this paper, we establish three new versions of Landau-type theorems for bounded bi-analytic functions of the form $F(z)=\bar{z}G(z)+H(z)$, where $G$ and $H$ are analytic in the unit disk $|z|<1$ with $G(0)=H(0)=0$ and $H'(0)=1$. In particular, two of them are sharp while the other one either generalizes or improves the corresponding result of Abdulhadi and Hajj. As consequences, several new sharp versions of Landau-type theorems for certain subclasses of bounded biharmonic mappings are proved.
△ Less
Submitted 15 February, 2023;
originally announced February 2023.
-
EPISODE: Episodic Gradient Clipping with Periodic Resampled Corrections for Federated Learning with Heterogeneous Data
Authors:
Michael Crawshaw,
Yajie Bao,
Mingrui Liu
Abstract:
Gradient clipping is an important technique for deep neural networks with exploding gradients, such as recurrent neural networks. Recent studies have shown that the loss functions of these networks do not satisfy the conventional smoothness condition, but instead satisfy a relaxed smoothness condition, i.e., the Lipschitz constant of the gradient scales linearly in terms of the gradient norm. Due…
▽ More
Gradient clipping is an important technique for deep neural networks with exploding gradients, such as recurrent neural networks. Recent studies have shown that the loss functions of these networks do not satisfy the conventional smoothness condition, but instead satisfy a relaxed smoothness condition, i.e., the Lipschitz constant of the gradient scales linearly in terms of the gradient norm. Due to this observation, several gradient clipping algorithms have been developed for nonconvex and relaxed-smooth functions. However, the existing algorithms only apply to the single-machine or multiple-machine setting with homogeneous data across machines. It remains unclear how to design provably efficient gradient clipping algorithms in the general Federated Learning (FL) setting with heterogeneous data and limited communication rounds. In this paper, we design EPISODE, the very first algorithm to solve FL problems with heterogeneous data in the nonconvex and relaxed smoothness setting. The key ingredients of the algorithm are two new techniques called \textit{episodic gradient clipping} and \textit{periodic resampled corrections}. At the beginning of each round, EPISODE resamples stochastic gradients from each client and obtains the global averaged gradient, which is used to (1) determine whether to apply gradient clipping for the entire round and (2) construct local gradient corrections for each client. Notably, our algorithm and analysis provide a unified framework for both homogeneous and heterogeneous data under any noise level of the stochastic gradient, and it achieves state-of-the-art complexity results. In particular, we prove that EPISODE can achieve linear speedup in the number of machines, and it requires significantly fewer communication rounds. Experiments on several heterogeneous datasets show the superior performance of EPISODE over several strong baselines in FL.
△ Less
Submitted 14 February, 2023;
originally announced February 2023.
-
Joint Edge-Model Sparse Learning is Provably Efficient for Graph Neural Networks
Authors:
Shuai Zhang,
Meng Wang,
Pin-Yu Chen,
Sijia Liu,
Songtao Lu,
Miao Liu
Abstract:
Due to the significant computational challenge of training large-scale graph neural networks (GNNs), various sparse learning techniques have been exploited to reduce memory and storage costs. Examples include \textit{graph sparsification} that samples a subgraph to reduce the amount of data aggregation and \textit{model sparsification} that prunes the neural network to reduce the number of trainab…
▽ More
Due to the significant computational challenge of training large-scale graph neural networks (GNNs), various sparse learning techniques have been exploited to reduce memory and storage costs. Examples include \textit{graph sparsification} that samples a subgraph to reduce the amount of data aggregation and \textit{model sparsification} that prunes the neural network to reduce the number of trainable weights. Despite the empirical successes in reducing the training cost while maintaining the test accuracy, the theoretical generalization analysis of sparse learning for GNNs remains elusive. To the best of our knowledge, this paper provides the first theoretical characterization of joint edge-model sparse learning from the perspective of sample complexity and convergence rate in achieving zero generalization error. It proves analytically that both sampling important nodes and pruning neurons with the lowest-magnitude can reduce the sample complexity and improve convergence without compromising the test accuracy. Although the analysis is centered on two-layer GNNs with structural constraints on data, the insights are applicable to more general setups and justified by both synthetic and practical citation datasets.
△ Less
Submitted 6 February, 2023;
originally announced February 2023.
-
Multiple solutions for a fractional Choquard problem with slightly subcritical exponents on bounded domains
Authors:
Marco G. Ghimenti,
Min Liu,
Zhongwei Tang
Abstract:
This paper is devoted to study a fractional Choquard problem with slightly subcritical exponents on bounded domains. When the exponent of the convolution type nonlinearity tends to the fractional critical one in the sense of Hardy-Littlewood-Sobolev inequality, we obtain the existence of multiple positive solutions via Lusternik-Schnirelmann category and nonlocal global compactness. Moreover, we p…
▽ More
This paper is devoted to study a fractional Choquard problem with slightly subcritical exponents on bounded domains. When the exponent of the convolution type nonlinearity tends to the fractional critical one in the sense of Hardy-Littlewood-Sobolev inequality, we obtain the existence of multiple positive solutions via Lusternik-Schnirelmann category and nonlocal global compactness. Moreover, we prove that the topology of the domain furnishes a lower bound for the number of positive solutions.
△ Less
Submitted 6 February, 2023;
originally announced February 2023.
-
Piecewise Temperleyan dimers and a multiple SLE$_8$
Authors:
Nathanaël Berestycki,
Mingchang Liu
Abstract:
We consider the dimer model on piecewise Temperleyan, simply connected domains, on families of graphs which include the square lattice as well as superposition graphs. We focus on the spanning tree $\mathcal{T}_δ$ associated to this model via Temperley's bijection, which turns out to be a Uniform Spanning Tree with singular alternating boundary conditions. Generalising the work of the second autho…
▽ More
We consider the dimer model on piecewise Temperleyan, simply connected domains, on families of graphs which include the square lattice as well as superposition graphs. We focus on the spanning tree $\mathcal{T}_δ$ associated to this model via Temperley's bijection, which turns out to be a Uniform Spanning Tree with singular alternating boundary conditions. Generalising the work of the second author with Peltola and Wu \cite{LiuPeltolaWuUST} we obtain a scaling limit result for $\mathcal{T}_δ$. For instance, in the simplest nontrivial case, the limit of $\mathcal{T}_δ$ is described by a pair of trees whose Peano curves are shown to converge jointly to a multiple SLE$_8$ pair. The interface between the trees is shown to be given by an SLE$_2(-1, \ldots, -1)$ curve. More generally we provide an equivalent description of the scaling limit in terms of imaginary geometry. This allows us to make use of the results developed by the first author and Laslier and Ray \cite{BLRdimers}. We deduce that, universally across these classes of graphs, the corresponding height function converges to a multiple of the Gaussian free field with boundary conditions that jump at each non-Temperleyan corner. After centering, this generalises a result of Russkikh \cite{RusskikhDimers} who proved it in the case of the square lattice. Along the way, we obtain results of independent interest on chordal hypergeometric SLE$_8$; for instance we show its law is equal to that of an SLE$_8 (\bar ρ)$ for a certain vector of force points, conditional on its hitting distribution on a specified boundary arc.
△ Less
Submitted 20 January, 2023;
originally announced January 2023.
-
Least-Squares Neural Network (LSNN) Method For Linear Advection-Reaction Equation: Discontinuity Interface
Authors:
Zhiqiang Cai,
Junpyo Choi,
Min Liu
Abstract:
We studied the least-squares ReLU neural network (LSNN) method for solving linear advection-reaction equation with discontinuous solution in [Cai, Zhiqiang, Jingshuang Chen, and Min Liu. ``Least-squares ReLU neural network (LSNN) method for linear advection-reaction equation.'' Journal of Computational Physics 443 (2021), 110514]. The method is based on a least-squares formulation and uses a new c…
▽ More
We studied the least-squares ReLU neural network (LSNN) method for solving linear advection-reaction equation with discontinuous solution in [Cai, Zhiqiang, Jingshuang Chen, and Min Liu. ``Least-squares ReLU neural network (LSNN) method for linear advection-reaction equation.'' Journal of Computational Physics 443 (2021), 110514]. The method is based on a least-squares formulation and uses a new class of approximating functions: ReLU neural network (NN) functions. A critical and additional component of the LSNN method, differing from other NN-based methods, is the introduction of a properly designed and physics preserved discrete differential operator.
In this paper, we study the LSNN method for problems with discontinuity interfaces. First, we show that ReLU NN functions with depth $\lceil \log_2(d+1)\rceil+1$ can approximate any $d$-dimensional step function on a discontinuity interface generated by a vector field as streamlines with any prescribed accuracy. By decomposing the solution into continuous and discontinuous parts, we prove theoretically that discretization error of the LSNN method using ReLU NN functions with depth $\lceil \log_2(d+1)\rceil+1$ is mainly determined by the continuous part of the solution provided that the solution jump is constant. Numerical results for both two- and three-dimensional test problems with various discontinuity interfaces show that the LSNN method with enough layers is accurate and does not exhibit the common Gibbs phenomena along discontinuity interfaces.
△ Less
Submitted 5 February, 2024; v1 submitted 15 January, 2023;
originally announced January 2023.
-
Privacy-Preserving Distributed Energy Resource Control with Decentralized Cloud Computing
Authors:
Xiang Huo,
Mingxi Liu
Abstract:
The rapidly growing penetration of renewable energy resources brings unprecedented challenges to power distribution networks - management of a large population of grid-tied controllable devices encounters control scalability crises and potential end-user privacy breaches. Despite the importance, research on privacy preservation of distributed energy resource (DER) control in a fully scalable manne…
▽ More
The rapidly growing penetration of renewable energy resources brings unprecedented challenges to power distribution networks - management of a large population of grid-tied controllable devices encounters control scalability crises and potential end-user privacy breaches. Despite the importance, research on privacy preservation of distributed energy resource (DER) control in a fully scalable manner is lacked. To fill the gap, this paper designs a novel decentralized privacy-preserving DER control framework that 1) achieves control scalability over DER population and heterogeneity; 2) eliminates peer-to-peer communications and secures the privacy of all participating DERs against various types of adversaries; and 3) enjoys higher computation efficiency and accuracy compared to state-of-the-art privacy-preserving methods. A strongly coupled optimization problem is formulated to control the power consumption and output of DERs, including solar photovoltaics and energy storage systems, then solved using the projected gradient method. Cloud computing and secret sharing are seamlessly integrated into the proposed decentralized computing to achieve privacy preservation. Simulation results prove the capabilities of the proposed approach in DER control applications.
△ Less
Submitted 5 January, 2023;
originally announced January 2023.
-
Higgs-Coulomb correspondence and Wall-Crossing in abelian GLSMs
Authors:
Konstantin Aleshkin,
Chiu-Chu Melissa Liu
Abstract:
We compute I-functions and central charges for abelian GLSMs using virtual matrix factorizations of Favero and Kim. In the Calabi-Yau case we provide analytic continuation for the central charges by explicit integral formulas. The integrals in question are called hemisphere partition functions and we call the integral representation Higgs-Coulomb correspondence. We then use it to prove GIT stabili…
▽ More
We compute I-functions and central charges for abelian GLSMs using virtual matrix factorizations of Favero and Kim. In the Calabi-Yau case we provide analytic continuation for the central charges by explicit integral formulas. The integrals in question are called hemisphere partition functions and we call the integral representation Higgs-Coulomb correspondence. We then use it to prove GIT stability wall-crossing for central charges.
△ Less
Submitted 3 January, 2023;
originally announced January 2023.
-
Boundary Regularity for an Even Order Elliptic System in the Critical Dimension
Authors:
Ming-Lun Liu,
Yao-Lan Tian
Abstract:
In this short note, we consider the Dirichlet problem associated to an even order elliptic system with antisymmetric first order potential. Given any continuous boundary data, we show that weak solutions are continuous up to boundary.
In this short note, we consider the Dirichlet problem associated to an even order elliptic system with antisymmetric first order potential. Given any continuous boundary data, we show that weak solutions are continuous up to boundary.
△ Less
Submitted 2 January, 2023;
originally announced January 2023.
-
Blow-up for semilinear wave equations on Kerr black hole backgrounds
Authors:
Mengyun Liu,
Chengbo Wang
Abstract:
We examine solutions to semilinear wave equations on black hole backgrounds and give a proof of an analog of the blow up part of the John theorem, with $F_p(u)=|u|^{p}$, on the Schwarzschild and Kerr black hole backgrounds. Concerning the case of Schwarzschild, we construct a class of small data, so that the solution blows up along the outgoing null cone, which applies for both $F_p(u)=|u|^{p}$ an…
▽ More
We examine solutions to semilinear wave equations on black hole backgrounds and give a proof of an analog of the blow up part of the John theorem, with $F_p(u)=|u|^{p}$, on the Schwarzschild and Kerr black hole backgrounds. Concerning the case of Schwarzschild, we construct a class of small data, so that the solution blows up along the outgoing null cone, which applies for both $F_p(u)=|u|^{p}$ and the focusing nonlinearity $F_p(u)=|u|^{p-1}u$. The proof suggests that the black hole does not have any essential influence on the formation of singularity, in the region away from the Cauchy horizon $r=r_-$ or the singularity $r=0$. Our approach is also robust enough to be adapted for general asymptotically flat space-time manifolds, possibly exterior to a compact domain, with spatial dimension $n\ge 2$. Typical examples include exterior domains, asymptotically Euclidean spaces, Reissner-Nördström space-times, and Kerr-Newman space-times.
△ Less
Submitted 29 December, 2022;
originally announced December 2022.