Showing 1–50 of 331 results for author: Liu, L

From bordisms of three-manifolds to domain walls between topological orders

Authors: Yu Leon Liu, Dalton A R Sakthivadivel

Abstract: We study a correspondence between spin three-manifolds and bosonic abelian topological orders. Let $N$ be a spin three-manifold. We can define a $(2+1)$-dimensional topological order $\mathrm{TO}_N$ as follows: its anyons are the torsion elements in $H_1(N)$, the braiding of anyons is given by the linking form, and their topological spins are given by the quadratic refinement of the linking form o… ▽ More We study a correspondence between spin three-manifolds and bosonic abelian topological orders. Let $N$ be a spin three-manifold. We can define a $(2+1)$-dimensional topological order $\mathrm{TO}_N$ as follows: its anyons are the torsion elements in $H_1(N)$, the braiding of anyons is given by the linking form, and their topological spins are given by the quadratic refinement of the linking form obtained from the spin structure. Under this correspondence, a surgery presentation of $N$ gives rise to a classical Chern--Simons description of the associated topological order $\mathrm{TO}_N$. We then extend the correspondence to spin bordisms between three-manifolds, and domain walls between topological orders. In particular, we construct a domain wall $\mathcal{D}_M$ between $\mathrm{TO}_N$ and $\mathrm{TO}_{N'}$, where $M$ is a spin bordism from $N$ to $N'$. This domain wall unfolds to a composition of a gapped boundary, obtained from anyon condensation, and a gapless Narain boundary CFT. △ Less

Submitted 15 July, 2024; originally announced July 2024.

Comments: 20+1 pages, five tikzpictures

MSC Class: 57K10; 57R56; 57R65; 81V27

Braiding on complex oriented Soergel bimodules

Authors: Yu Leon Liu

Abstract: In this note, we study U(n) Soergel bimodules in the context of stable homotopy theory. We define the $(\infty, 1)$-category $\mathrm{SBim}_E(n)$ of $E$-valued U(n) Soergel bimodules, where $E$ is a connective $\mathbb{E}_\infty$-ring spectrum, and assemble them into a monoidal locally additive $(\infty, 2)$-category $\mathrm{SBim}_E$. When $E$ has a complex orientation, we then construct a braidi… ▽ More In this note, we study U(n) Soergel bimodules in the context of stable homotopy theory. We define the $(\infty, 1)$-category $\mathrm{SBim}_E(n)$ of $E$-valued U(n) Soergel bimodules, where $E$ is a connective $\mathbb{E}_\infty$-ring spectrum, and assemble them into a monoidal locally additive $(\infty, 2)$-category $\mathrm{SBim}_E$. When $E$ has a complex orientation, we then construct a braiding, i.e. an $\mathbb{E}_2$-algebra structure, on the universal locally stable $(\infty, 2)$-category $\mathrm{K}^b_{\mathrm{loc}}(\mathrm{SBim}_E)$ associated to $\mathrm{SBim}_E$. Along the way, we also prove spectral analogs of standard splittings of Soergel bimodules. This is a topological generalization of the type $A$ Soergel bimodule theory developed in a previous paper. △ Less

Submitted 5 July, 2024; originally announced July 2024.

Comments: 31 pages, comments welcome

On a class of multi-fidelity methods for the semiclassical Schrödinger equation with uncertainties

Authors: Yiwen Lin, Liu Liu

Abstract: In this paper, we study the semiclassical Schrödinger equation with random parameters and develop several robust multi-fidelity methods. We employ the time-splitting Fourier pseudospectral (TSFP) method for the high-fidelity solver, and consider different low-fidelity solvers including the meshless method like frozen Gaussian approximation (FGA) and the level set (LS) method for the semiclassical… ▽ More In this paper, we study the semiclassical Schrödinger equation with random parameters and develop several robust multi-fidelity methods. We employ the time-splitting Fourier pseudospectral (TSFP) method for the high-fidelity solver, and consider different low-fidelity solvers including the meshless method like frozen Gaussian approximation (FGA) and the level set (LS) method for the semiclassical limit of the Schrödinger equation. With a careful choice of the low-fidelity model, we obtain an error estimate for the bi-fidelity method. We conduct numerous numerical experiments and validate the accuracy and efficiency of our proposed multi-fidelity methods, by comparing the performance of a class of bi-fidelity and tri-fidelity approximations. △ Less

Submitted 8 June, 2024; originally announced June 2024.

MSC Class: 35J10; 65M70

A specialisation theorem for Lang-Néron groups

Authors: Bruno Kahn, Long Liu

Abstract: We show that, for a polarised smooth projective variety $B \hookrightarrow \mathbb{P}^n_k$ of dimension $\geq 2$ over an infinite field $k$ and an abelian variety $A$ over the function field of $B$, there exists a dense Zariski open set of smooth geometrically connected hyperplane sections $h$ of $B$ such that $A$ has good reduction at $h$ and the specialisation homomorphism of Lang-Néron groups a… ▽ More We show that, for a polarised smooth projective variety $B \hookrightarrow \mathbb{P}^n_k$ of dimension $\geq 2$ over an infinite field $k$ and an abelian variety $A$ over the function field of $B$, there exists a dense Zariski open set of smooth geometrically connected hyperplane sections $h$ of $B$ such that $A$ has good reduction at $h$ and the specialisation homomorphism of Lang-Néron groups at $h$ is injective (up to a finite $p$-group in positive characteristic $p$). This gives a positive answer to a conjecture of the first author, which is used to deduce a negative definiteness result on his refined height pairing. This also sheds a new light on Néron's specialisation theorem. △ Less

Submitted 27 June, 2024; v1 submitted 9 May, 2024; originally announced May 2024.

Comments: Comments welcome

MSC Class: 11G99; 14K99

The Smith Fiber Sequence of Invertible Field Theories

Authors: Arun Debray, Sanath K. Devalapurkar, Cameron Krulewski, Yu Leon Liu, Natalia Pacheco-Tallaj, Ryan Thorngren

Abstract: Smith homomorphisms are maps between bordism groups that change both the dimension and the tangential structure. We give a completely general account of Smith homomorphisms, unifying the many examples in the literature. We provide three definitions of Smith homomorphisms, including as maps of Thom spectra, and show they are equivalent. Using this, we identify the cofiber of the spectrum-level Smit… ▽ More Smith homomorphisms are maps between bordism groups that change both the dimension and the tangential structure. We give a completely general account of Smith homomorphisms, unifying the many examples in the literature. We provide three definitions of Smith homomorphisms, including as maps of Thom spectra, and show they are equivalent. Using this, we identify the cofiber of the spectrum-level Smith map and extend the Smith homomorphism to a long exact sequence of bordism groups, which is a powerful computation tool. We discuss several examples of this long exact sequence, relating them to known constructions such as Wood's and Wall's sequences. Furthermore, taking Anderson duals yields a long exact sequence of invertible field theories, which has a rich physical interpretation. We developed the theory in this paper with applications in mind to symmetry breaking in quantum field theory, which we study in a companion paper. △ Less

Submitted 6 May, 2024; originally announced May 2024.

Comments: 73 pages, 3 figures. Originally posted as the mathematical sections of arXiv:2309.16749. This version features new content in sections 6 and 8

On the Sundman-Sperling estimates for the restricted one-center-two-body problem

Authors: Ku-Jung Hsu, Lei Liu

Abstract: In the past two decades, since the discovery of the figure-8 orbit by Chenciner and Montgomery, the variational method has became one of the most popular tools for constructing new solutions of the $N$-body problem and its extended problems. However, finding solutions to the restricted three-body problem, in particular, the two primaries form a collision Kepler system, remains a great difficulty.… ▽ More In the past two decades, since the discovery of the figure-8 orbit by Chenciner and Montgomery, the variational method has became one of the most popular tools for constructing new solutions of the $N$-body problem and its extended problems. However, finding solutions to the restricted three-body problem, in particular, the two primaries form a collision Kepler system, remains a great difficulty. One of the major reasons is the essential differences between two-body collisions and three-body collisions. In this paper, we consider a similar three-body system with less difficulty, i.e. the restricted one-center-two-body system, that is involving a massless particle and a collision Kepler system with one body fixed. It is an intermediate system between the restricted three-body problem and the two-center problem. By an in-depth analysis of the asymptotic behavior of the minimizer, and an argument of critical and infliction points, we prove the Sundman-Sperling estimates near the three-body collision for the minimizers. With these estimates, we provide a class of collision-free solutions with prescribed boundary angles. Finally, under the extended collision Kepler system from Gordon, we constructed a family of periodic and quasi-periodic solutions. △ Less

Submitted 21 April, 2024; originally announced April 2024.

Comments: 37 pages,5 figures

MSC Class: 70F07; 70F15; 70F16; 70G75; 70M20

Chow trace of 1-motives and the Lang-Néron groups

Authors: Long Liu

Abstract: We show that in the case of primary field extensions, the extension of scalars of Deligne $1$-motives admits a left adjoint, called Chow image, and a right adjoint, called Chow trace. This generalizes Chow's results on abelian varieties. Then we study the Chow trace in the framework of Voevodsky's triangulated categories of (étale) motives. With respect to the $1$-motivic $t$-structure on the cate… ▽ More We show that in the case of primary field extensions, the extension of scalars of Deligne $1$-motives admits a left adjoint, called Chow image, and a right adjoint, called Chow trace. This generalizes Chow's results on abelian varieties. Then we study the Chow trace in the framework of Voevodsky's triangulated categories of (étale) motives. With respect to the $1$-motivic $t$-structure on the category of Voevodsky's homological $1$-motives, the zero-th direct image of an abelian variety is given by the Chow trace, and the first direct image is the $0$-motive defined by the (geometric) Lang-Néron group. △ Less

Submitted 19 April, 2024; originally announced April 2024.

Comments: 54 pages, comments welcome

MSC Class: 14C15; 14F20; 14F42; 14K05

Urban Water Sprinkler Routing: A Multi-Depot Mixed Capacitated Arc Routing Problem Incorporating Real-Time Demands

Authors: Hongtai Yang, Luna Liu, Ke Han, Boyi Lei

Abstract: Fugitive road dust (FRD), as one of the major pollutants in the city, poses great harm to the environment and the physical health of citizens. A common countermeasure adopted by government agencies is employing on-road water trucks (sprinklers) to spray water (sprinkle) on urban streets to reduce the FRD. Currently, the traveling routes of sprinklers are usually planned based on drivers' experienc… ▽ More Fugitive road dust (FRD), as one of the major pollutants in the city, poses great harm to the environment and the physical health of citizens. A common countermeasure adopted by government agencies is employing on-road water trucks (sprinklers) to spray water (sprinkle) on urban streets to reduce the FRD. Currently, the traveling routes of sprinklers are usually planned based on drivers' experience, which may lead low operation efficiency and could not respond to the real-time sprinkling demands. To address these issues, this study formulates the routes planning of sprinklers as a multi-depot mixed capacitated arc routing problem with real-time demands with the aim of minimizing the sprinklers' travel distance. We develop an improved adaptive large neighborhood search (ALNS) algorithm that incorporates a tabu-list and a perturbation mechanism to solve this problem. Furthermore, a problem-specific acceleration mechanism is designed to reduce unnecessary search domains to improve the efficiency of the algorithm. Empirical experiments are conducted based on various scenarios and the results demonstrate that the proposed algorithm generates solutions that are superior or at least comparable to the solutions generated by the traditional ALNS algorithm but with significantly lower computation time. Sensitivity analysis is conducted to explore the effects of relevant parameters on the results. This study is the first to incorporate real-time FRD pollution information, gathered through multiple data sources via IoT technology, into urban sprinkling operations, extending the traditional CARP from a tactical planning to a real-time operational environment. A real-world implementation case is also presented. △ Less

Submitted 15 April, 2024; originally announced April 2024.

Comments: 29 pages, 5 figures

The paradigm of tax-reward and tax-punishment strategies in the advancement of public resource management dynamics

Authors: Lichen Wang, Yuyuan Liu, Ruqiang Guo, Liang Zhang, Linjie Liu, Shijia Hua

Abstract: In contemporary society, the effective utilization of public resources remains a subject of significant concern. A common issue arises from defectors seeking to obtain an excessive share of these resources for personal gain, potentially leading to resource depletion. To mitigate this tragedy and ensure sustainable development of resources, implementing mechanisms to either reward those who adhere… ▽ More In contemporary society, the effective utilization of public resources remains a subject of significant concern. A common issue arises from defectors seeking to obtain an excessive share of these resources for personal gain, potentially leading to resource depletion. To mitigate this tragedy and ensure sustainable development of resources, implementing mechanisms to either reward those who adhere to distribution rules or penalize those who do not, appears advantageous. We introduce two models: a tax-reward model and a tax-punishment model, to address this issue. Our analysis reveals that in the tax-reward model, the evolutionary trajectory of the system is influenced not only by the tax revenue collected but also by the natural growth rate of the resources. Conversely, the tax-punishment model exhibits distinct characteristics when compared to the tax-reward model, notably the potential for bistability. In such scenarios, the selection of initial conditions is critical, as it can determine the system's path. Furthermore, our study identifies instances where the system lacks stable points, exemplified by a limit cycle phenomenon, underscoring the complexity and dynamism inherent in managing public resources using these models. △ Less

Submitted 31 March, 2024; originally announced April 2024.

Comments: Accepted by Proceedings of the Royal Society B-Biological Sciences

Arc-transitive maps with coprime Euler characteristic and edge number

Authors: C. H. Li, Lu Yi Liu

Abstract: This is one of a series of papers which aim towards a classification of edge-transitive maps of which the Euler characteristic and the edge number are coprime. This one carries out the classification work for arc-transitive maps with nonsolvable automorphism groups, which illustrates how the edge number impacts on the Euler characteristic for maps. The classification is involved with the construct… ▽ More This is one of a series of papers which aim towards a classification of edge-transitive maps of which the Euler characteristic and the edge number are coprime. This one carries out the classification work for arc-transitive maps with nonsolvable automorphism groups, which illustrates how the edge number impacts on the Euler characteristic for maps. The classification is involved with the construction of some new and interesting arc-regular maps. △ Less

Submitted 27 March, 2024; originally announced March 2024.

MSC Class: 05C10; 20B25

DG singular equivalence and singular locus

Authors: Leilei Liu, Jieheng Zeng

Abstract: For a commutative Gorenstein Noetherian ring $R$, we construct an affine scheme $X$ solely from DG singularity category $S_{dg}(R)$ of $R$ such that there is a finite surjective morphism $X \rightarrow \mathrm{Spec}(R /I)$, where $\mathrm{Spec}(R /I)$ is the singular locus in $\mathrm{Spec}(R)$. As an application, for two such rings with equivalent DG singularity categories, we prove that the sing… ▽ More For a commutative Gorenstein Noetherian ring $R$, we construct an affine scheme $X$ solely from DG singularity category $S_{dg}(R)$ of $R$ such that there is a finite surjective morphism $X \rightarrow \mathrm{Spec}(R /I)$, where $\mathrm{Spec}(R /I)$ is the singular locus in $\mathrm{Spec}(R)$. As an application, for two such rings with equivalent DG singularity categories, we prove that the singular loci in their affine schemes have the same dimension. △ Less

Submitted 31 March, 2024; v1 submitted 20 March, 2024; originally announced March 2024.

Comments: 18 pages

Frequency-Reactive Power Optimization Strategy of Grid-forming Offshore Wind Farm Using DRU-HVDC Transmission

Authors: Zhekai Li, Kun Han, Xu Cai, Renxin Yang, Haotian Yu, Kepeng Xia, Lulu Liu

Abstract: The diode rectifier unit-based high voltage direct current (DRU-HVDC) transmission with grid-forming (GFM) wind turbine is becoming a promising scheme for offshore wind farm(OWF) integration due to its high reliability and low cost. In this scheme, the AC network of the OWF and the DRU has completely different synchronization mechanisms and power flow characteristics from the traditional power sys… ▽ More The diode rectifier unit-based high voltage direct current (DRU-HVDC) transmission with grid-forming (GFM) wind turbine is becoming a promising scheme for offshore wind farm(OWF) integration due to its high reliability and low cost. In this scheme, the AC network of the OWF and the DRU has completely different synchronization mechanisms and power flow characteristics from the traditional power system. To optimize the power flow and reduce the net loss, this paper carries out the power flow modeling and optimization analysis for the DRU-HVDC transmission system with grid-forming OWFs. The influence of the DRU and the GFM wind turbines on the power flow of the system is analyzed. On this basis, improved constraint conditions are proposed and an optimal power flow (OPF) method is established. This method can minimize the power loss by adjusting the reactive power output of each wind turbine and internal network frequency. Finally, based on MATLAB, this paper uses YALMIP toolkit and CPLEX mathematical solver to realize the programming solution of the OPF model proposed in this paper. The results show that the proposed optimization strategy can effectively reduce the power loss of the entire OWF and the transmission system with an optimization ratio of network losses exceeding 25.3%. △ Less

Submitted 16 March, 2024; originally announced March 2024.

Comments: 6 pages, 7 figures, to be published in the 7th IEEE Conference on Energy Internet and Energy System Integration (EI2 2023)

Second-Order Strong Optimality and Second-Order Duality for Nonsmooth Constrained Multiobjective Fractional Programming Problems

Authors: Jiawei Chen, Luyu Liu, Yibing Lv, Debdas Ghosh, Jen-Chih Yao

Abstract: This paper investigates constrained nonsmooth multiobjective fractional programming problem (NMFP) in real Banach spaces. It derives a quotient calculus rule for computing the first- and second-order Clarke derivatives of fractional functions involving locally Lipschitz functions. A novel second-order Abadie-type regularity condition is presented, defined with the help of the Clarke directional de… ▽ More This paper investigates constrained nonsmooth multiobjective fractional programming problem (NMFP) in real Banach spaces. It derives a quotient calculus rule for computing the first- and second-order Clarke derivatives of fractional functions involving locally Lipschitz functions. A novel second-order Abadie-type regularity condition is presented, defined with the help of the Clarke directional derivative and the Páles-Zeidan second-order directional derivative. We establish both first- and second-order strong necessary optimality conditions, which contain some new information on multipliers and imply the strong KKT necessary conditions, for a Borwein-type properly efficient solution of NMFP by utilizing generalized directional derivatives. Moreover, it derives second-order sufficient optimality conditions for NMFP under a second-order generalized convexity assumption. Additionally, we derive duality results between NMFP and its second-order dual problem under some appropriate conditions. △ Less

Submitted 15 March, 2024; originally announced March 2024.

Comments: 23 pages

MSC Class: 90C46 and 90C29 and 49J52

Topology and entanglement of molecular phase space

Authors: Victor V. Albert, Eric Kubischta, Mikhail Lemeshko, Lee R. Liu

Abstract: We formulate a quantum phase space for molecular rotational and nuclear-spin states. Taking in molecular geometry and nuclear-spin data, our framework yields admissible position and momentum states, inter-convertible via a generalized Fourier transform. We classify molecules into three types -- asymmetric, rotationally symmetric, and perrotationally symmetric -- with the last type having no macros… ▽ More We formulate a quantum phase space for molecular rotational and nuclear-spin states. Taking in molecular geometry and nuclear-spin data, our framework yields admissible position and momentum states, inter-convertible via a generalized Fourier transform. We classify molecules into three types -- asymmetric, rotationally symmetric, and perrotationally symmetric -- with the last type having no macroscopic analogue due to nuclear-spin statistics constraints. We identify two features in perrotationally symmetric state spaces that are Hamiltonian-independent and induced solely by symmetry and spin statistics. First, many molecular species are intrinsically rotation-spin entangled in a way that cannot be broken without transitioning to another species or breaking symmetry. Second, each molecular position state houses an internal pseudo-spin or "fiber" degree of freedom, and the fiber's Berry phase or matrix after adiabatic changes in position yields naturally robust operations, akin to braiding anyonic quasiparticles or realizing fault-tolerant quantum gates. We outline scenarios where these features can be experimentally probed. △ Less

Submitted 23 March, 2024; v1 submitted 7 March, 2024; originally announced March 2024.

Comments: 9 + 37 pages, 8 figures, 6 tables, 46 examples; v2 minor clarifications

Spectral convergence of a semi-discretized numerical system for the spatially homogeneous Boltzmann equation with uncertainties

Authors: Liu Liu, Kunlun Qi

Abstract: In this paper, we study the Boltzmann equation with uncertainties and prove that the spectral convergence of the semi-discretized numerical system holds in a combined velocity and random space, where the Fourier-spectral method is applied for approximation in the velocity space whereas the generalized polynomial chaos (gPC)-based stochastic Galerkin (SG) method is employed to discretize the random… ▽ More In this paper, we study the Boltzmann equation with uncertainties and prove that the spectral convergence of the semi-discretized numerical system holds in a combined velocity and random space, where the Fourier-spectral method is applied for approximation in the velocity space whereas the generalized polynomial chaos (gPC)-based stochastic Galerkin (SG) method is employed to discretize the random variable. Our proof is based on a delicate energy estimate for showing the well-posedness of the numerical solution as well as a rigorous control of its negative part in our well-designed functional space that involves high-order derivatives of both the velocity and random variables. This paper rigorously justifies the statement proposed in [Remark 4.4, J. Hu and S. Jin, J. Comput. Phys., 315 (2016), pp. 150-168]. △ Less

Submitted 6 May, 2024; v1 submitted 10 February, 2024; originally announced February 2024.

Comments: Revised version. To appear in SIAM/ASA Journal on Uncertainty Quantification

MSC Class: Primary 35Q20; 65M12; Secondary 65M70; 45G10

Data-induced multiscale losses and efficient multirate gradient descent schemes

Authors: Juncai He, Liangchen Liu, Yen-Hsi Richard Tsai

Abstract: This paper investigates the impact of multiscale data on machine learning algorithms, particularly in the context of deep learning. A dataset is multiscale if its distribution shows large variations in scale across different directions. This paper reveals multiscale structures in the loss landscape, including its gradients and Hessians inherited from the data. Correspondingly, it introduces a nove… ▽ More This paper investigates the impact of multiscale data on machine learning algorithms, particularly in the context of deep learning. A dataset is multiscale if its distribution shows large variations in scale across different directions. This paper reveals multiscale structures in the loss landscape, including its gradients and Hessians inherited from the data. Correspondingly, it introduces a novel gradient descent approach, drawing inspiration from multiscale algorithms used in scientific computing. This approach seeks to transcend empirical learning rate selection, offering a more systematic, data-informed strategy to enhance training efficiency, especially in the later stages. △ Less

Submitted 6 February, 2024; v1 submitted 5 February, 2024; originally announced February 2024.

Comments: 28 pages, 4 figures, submitted under review

MSC Class: 65F10; 65F45; 68T07 ACM Class: G.1.6; I.2.6

A characterization of extremal non-transmission-regular graphs by the distance (signless Laplacian) spectral radius

Authors: Jingfen Lan, Lele Liu

Abstract: Let $G$ be a simple connected graph of order $n$ and $\partial(G)$ is the spectral radius of the distance matrix $D(G)$ of $G$. The transmission $D_i$ of vertex $i$ is the $i$-th row sum of $D(G)$. Denote by $D_{\max}(G)$ the maximum of transmissions over all vertices of $G$, and $\partial^Q(G)$ is the spectral radius of the distance signless Laplacian matrix… ▽ More Let $G$ be a simple connected graph of order $n$ and $\partial(G)$ is the spectral radius of the distance matrix $D(G)$ of $G$. The transmission $D_i$ of vertex $i$ is the $i$-th row sum of $D(G)$. Denote by $D_{\max}(G)$ the maximum of transmissions over all vertices of $G$, and $\partial^Q(G)$ is the spectral radius of the distance signless Laplacian matrix $D(G)+\mbox{diag}(D_1,D_2,\ldots,D_n)$. In this paper, we present a sharp lower bound of $2D_{\max}(G)-\partial^Q(G)$ among all $n$-vertex connected graphs, and characterize the extremal graphs. Furthermore, we give the minimum values of respective $D_{\max}(G)-\partial(G)$ and $2D_{\max}(G)-\partial^Q(G)$ on trees and characterize the extremal trees. △ Less

Submitted 1 February, 2024; originally announced February 2024.

MSC Class: 05C50

On the supporting quasi-hyperplane and separation theorem of geodesic convex sets with applications on Riemannian manifolds

Authors: Li-wen Zhou, Ling-ling Liu, Chao Min, Yao-jia Zhang, Nan-Jing Huang

Abstract: In this paper, we first establish the separation theorem between a point and a locally geodesic convex set and then prove the existence of a supporting quasi-hyperplane at any point on the boundary of the closed locally geodesic convex set on a Riemannian manifold. As applications, some optimality conditions are obtained for optimization problems with constraints on Riemannian manifolds. In this paper, we first establish the separation theorem between a point and a locally geodesic convex set and then prove the existence of a supporting quasi-hyperplane at any point on the boundary of the closed locally geodesic convex set on a Riemannian manifold. As applications, some optimality conditions are obtained for optimization problems with constraints on Riemannian manifolds. △ Less

Submitted 24 January, 2024; originally announced January 2024.

A braided monoidal $(\infty,2)$-category of Soergel bimodules

Authors: Yu Leon Liu, Aaron Mazel-Gee, David Reutter, Catharina Stroppel, Paul Wedrich

Abstract: The Hecke algebras for all symmetric groups taken together form a braided monoidal category that controls all quantum link invariants of type A and, by extension, the standard canon of topological quantum field theories in dimension 3 and 4. Here we provide the first categorification of this Hecke braided monoidal category, which takes the form of an $\mathbb{E}_2$-monoidal $(\infty,2)$-category w… ▽ More The Hecke algebras for all symmetric groups taken together form a braided monoidal category that controls all quantum link invariants of type A and, by extension, the standard canon of topological quantum field theories in dimension 3 and 4. Here we provide the first categorification of this Hecke braided monoidal category, which takes the form of an $\mathbb{E}_2$-monoidal $(\infty,2)$-category whose hom-$(\infty,1)$-categories are $k$-linear, stable, idempotent-complete, and equipped with $\mathbb{Z}$-actions. This categorification is designed to control homotopy-coherent link homology theories and to-be-constructed topological quantum field theories in dimension 4 and 5. Our construction is based on chain complexes of Soergel bimodules, with monoidal structure given by parabolic induction and braiding implemented by Rouquier complexes, all modelled homotopy-coherently. This is part of a framework which allows to transfer the toolkit of the categorification literature into the realm of $\infty$-categories and higher algebra. Along the way, we develop families of factorization systems for $(\infty,n)$-categories, enriched $\infty$-categories, and $\infty$-operads, which may be of independent interest. As a service aimed at readers less familiar with homotopy-coherent mathematics, we include a brief introduction to the necessary $\infty$-categorical technology in the form of an appendix. △ Less

Submitted 6 February, 2024; v1 submitted 5 January, 2024; originally announced January 2024.

Comments: 142 pages, comments welcome, v2 with minor change in title

Variants of spectral Turán theorems and eigenvectors of graphs

Authors: Lele Liu, Bo Ning

Abstract: In 2002, Nikiforov proved that for an $n$-vertex graph $G$ with clique number $ω$ and edge number $m$, the spectral radius $λ(G)$ satisfies $λ(G) \leq \sqrt{2(1 - 1/ω) m}$, which confirmed a conjecture implicitly suggested by Edwards and Elphick. In this paper, we prove a local version of spectral Turán inequality, which states that $λ^2(G)\leq 2\sum_{e\in E(G)}\frac{c(e)-1}{c(e)}$, where $c(e)$ i… ▽ More In 2002, Nikiforov proved that for an $n$-vertex graph $G$ with clique number $ω$ and edge number $m$, the spectral radius $λ(G)$ satisfies $λ(G) \leq \sqrt{2(1 - 1/ω) m}$, which confirmed a conjecture implicitly suggested by Edwards and Elphick. In this paper, we prove a local version of spectral Turán inequality, which states that $λ^2(G)\leq 2\sum_{e\in E(G)}\frac{c(e)-1}{c(e)}$, where $c(e)$ is the order of the largest clique containing the edge $e$ in $G$. We also characterize the extremal graphs. We prove that our theorem implies Nikiforov's theorem and give an example to show that the difference of Nikiforov's bound and ours is $Ω(\sqrt{m})$ for some cases. Additionally, we establish a spectral counterpart to Ore's problem (1962) which asks for the maximum size of an $n$-vertex graph such that its complement is connected and does not contain $F$ as a subgraph. Our result leads to a new spectral Turán inequality applicable to graphs with connected complements. Finally, we disprove a conjecture of Gregory, asserting that for a connected $n$-vertex graph $G$ with chromatic number $k\geq 2$ and an independent set $S$, we have \[ \sum_{v\in S} x_v^2 \leq \frac{1}{2} - \frac{k-2}{2\sqrt{(k-2)^2 + 4(k-1)(n-k+1)}}, \] where $x_v$ is the component of the Perron vector of $G$ with respect to the vertex $v$. A modified version of Gregory's conjecture is proposed. △ Less

Submitted 19 March, 2024; v1 submitted 26 December, 2023; originally announced December 2023.

Comments: 20 pages. This is a new version of the previous paper titled "A local version of spectral Turán theorem". In this version, we add more results, including a variant of spectral Turán theorem and a disproof of a conjecture of Gregory

Nonlinear stability of planar shock wave to 3-D compressible Navier-Stokes equations in half space with Navier Boundary conditions

Authors: Lin Chang, Lingjun Liu, Lingda Xu

Abstract: In this paper, we consider the large time behavior of planar shock wave for 3-D compressible isentropic Navier-Stokes equations (CNS) in half space. Providing the strength of the shock wave and initial perturbations are small, we proved the planar shock wave for 3-D CNS is nonlinearly stable in half space with Navier boundary condition. The main difficulty comes from the compressibility of shock… ▽ More In this paper, we consider the large time behavior of planar shock wave for 3-D compressible isentropic Navier-Stokes equations (CNS) in half space. Providing the strength of the shock wave and initial perturbations are small, we proved the planar shock wave for 3-D CNS is nonlinearly stable in half space with Navier boundary condition. The main difficulty comes from the compressibility of shock wave, which leads to lower order terms with bad sign, see the third line in \cref{C17}. We apply a decomposition of the solution into zero and non-zero modes: we take the anti-derivative for the zero mode and obtain the space-time estimates for the energy of perturbation itself. Then combining the fact that the Poincaré inequality is available for the non-zero mode, we have successfully controlled the lower order terms with bad sign in \cref{C17}. To overcome the difficulty that comes from the boundary, we introduce the two crucial estimates on boundary \cref{CLem0} and fully utilize the property of Navier boundary conditions, which means that the normal velocity is zero on the boundary and the fluid tangential velocity is proportional to the tangential component of the viscous stress tensor on the boundary. Finally, the nonlinear stability is proved by the weighted energy method. △ Less

Submitted 9 December, 2023; originally announced December 2023.

MSC Class: 35Q30; 76L05; 76N06

Decay rate to the planar viscous shock wave for multi-dimensional scalar conservation laws

Authors: Lingjun Liu, Shu Wang, Lingda Xu

Abstract: In this paper, we study the time-decay rate toward the planar viscous shock wave for multi-dimensional (m-d) scalar viscous conservation law. We first decompose the perturbation into zero and non-zero mode, and then introduce the anti-derivative of the zero mode. Though an $L^p$ estimate and the area inequality introduced in \cite{DHS2020}, we obtained the decay rate for planar shock wave for n-d… ▽ More In this paper, we study the time-decay rate toward the planar viscous shock wave for multi-dimensional (m-d) scalar viscous conservation law. We first decompose the perturbation into zero and non-zero mode, and then introduce the anti-derivative of the zero mode. Though an $L^p$ estimate and the area inequality introduced in \cite{DHS2020}, we obtained the decay rate for planar shock wave for n-d scalar viscous conservation law for all $n\geq1$. The initial perturbations we studied are small, i.e., $\|Φ_0\|_{H^2}\bigcap\|Φ_0\|_{L^p}\le \varepsilon$, where $Φ_0$ is the anti-derivative of the zero mode of initial perturbation and $\varepsilon$ is a small constant, see \cref{antiderivative}. It is noted that there is no additional requirement on $Φ_0$, i.e., $Φ_0(x_1)$ only belongs to $H^2(\R)$. Thus, there are essential differences from previous results, in which the initial data is required to belong to some weighted Sobolev space, cf.\cite{Goo1989,KM1985}. Moreover, the exponential decay rate of the non-zero mode is also obtained. △ Less

Submitted 6 December, 2023; originally announced December 2023.

MSC Class: 35L65; 35L67; 35K15

Directional Differentiability of the Generalized Metric Projection in Hilbert spaces and Hilbertian Bochner spaces

Authors: Jinlu Li, Li Cheng, Lishan Liu, Linsen Xie

Abstract: Let $H$ be a real Hilbert space and $C$ a nonempty closed and convex subset of $H$. Let $P_C: H\rightarrow C$ denote the (standard) metric projection operator. In this paper, we study the Gâteaux directional differentiability of $P_C$ and investigate some of its properties. The Gâteaux directionally derivatives of $P_C$ are precisely given for the following cases of the considered subset $C$: 1. c… ▽ More Let $H$ be a real Hilbert space and $C$ a nonempty closed and convex subset of $H$. Let $P_C: H\rightarrow C$ denote the (standard) metric projection operator. In this paper, we study the Gâteaux directional differentiability of $P_C$ and investigate some of its properties. The Gâteaux directionally derivatives of $P_C$ are precisely given for the following cases of the considered subset $C$: 1. closed and convex subsets; 2. closed balls; 3. closed and convex cones (including proper closed subspaces). For special Hilbert spaces, we consider directional differentiability of $P_C$ for some Hilbert spaces with orthonormal bases and the real Hilbert space $L^2([-π,π])$ with the trigonometric orthonormal basis. △ Less

Submitted 24 October, 2023; originally announced October 2023.

Comments: This article has been accepted for publication

MSC Class: 49J50; 26A24; 47A58; 47J30; 49J40

$L^2$-Wasserstein contraction for Euler schemes of elliptic diffusions and interacting particle systems

Authors: Linshan Liu, Mateusz B. Majka, Pierre Monmarché

Abstract: We show the $L^2$-Wasserstein contraction for the transition kernel of a discretised diffusion process, under a contractivity at infinity condition on the drift and a sufficiently high diffusivity requirement. This extends recent results that, under similar assumptions on the drift but without the diffusivity restrictions, showed the $L^1$-Wasserstein contraction, or $L^p$-Wasserstein bounds for… ▽ More We show the $L^2$-Wasserstein contraction for the transition kernel of a discretised diffusion process, under a contractivity at infinity condition on the drift and a sufficiently high diffusivity requirement. This extends recent results that, under similar assumptions on the drift but without the diffusivity restrictions, showed the $L^1$-Wasserstein contraction, or $L^p$-Wasserstein bounds for $p > 1$ that were, however, not true contractions. We explain how showing the true $L^2$-Wasserstein contraction is crucial for obtaining the local Poincaré inequality for the transition kernel of the Euler scheme of a diffusion. Moreover, we discuss other consequences of our contraction results, such as concentration inequalities and convergence rates in KL-divergence and total variation. We also study the corresponding $L^2$-Wasserstein contraction for discretisations of interacting diffusions. As a particular application, this allows us to analyse the behaviour of particle systems that can be used to approximate a class of McKean-Vlasov SDEs that were recently studied in the mean-field optimization literature. △ Less

Submitted 24 October, 2023; originally announced October 2023.

Comments: 28 pages

Extreme Learning Machine-Assisted Solution of Biharmonic Equations via Its Coupled Schemes

Authors: Xi'an Li, Jinran Wu, Jiaxin Deng, Zhe Ding, You-Gan Wang, Xin Tai, Liang Liu

Abstract: Obtaining the solutions of partial differential equations based on various machine learning methods has drawn more and more attention in the fields of scientific computation and engineering applications. In this work, we first propose a coupled Extreme Learning Machine (called CELM) method incorporated with the physical laws to solve a class of fourth-order biharmonic equations by reformulating it… ▽ More Obtaining the solutions of partial differential equations based on various machine learning methods has drawn more and more attention in the fields of scientific computation and engineering applications. In this work, we first propose a coupled Extreme Learning Machine (called CELM) method incorporated with the physical laws to solve a class of fourth-order biharmonic equations by reformulating it into two well-posed Poisson problems. In addition, some activation functions including tangent, gauss, sine, and trigonometric (sin+cos) functions are introduced to assess our CELM method. Notably, the sine and trigonometric functions demonstrate a remarkable ability to effectively minimize the approximation error of the CELM model. In the end, several numerical experiments are performed to study the initializing approaches for both the weights and biases of the hidden units in our CELM model and explore the required number of hidden units. Numerical results show the proposed CELM algorithm is high-precision and efficient to address the biharmonic equation in both regular and irregular domains. △ Less

Submitted 26 October, 2023; v1 submitted 21 October, 2023; originally announced October 2023.

Optimal decay rates to the contact wave for 1-D compressible Navier-Stokes equations

Authors: Lingjun Liu, Shu Wang, Lingda Xu

Abstract: This paper investigates the decay rates of the contact wave in one-dimensional Navier-Stokes equations. We study two cases of perturbations, with and without zero mass condition, i.e., the integration of initial perturbations is zero and non-zero, respectively. For the case without zero mass condition, we obtain the optimal decay rate $(1+t)^{-\frac{1}{2}}$ for the perturbation in $L^\infty$ norm,… ▽ More This paper investigates the decay rates of the contact wave in one-dimensional Navier-Stokes equations. We study two cases of perturbations, with and without zero mass condition, i.e., the integration of initial perturbations is zero and non-zero, respectively. For the case without zero mass condition, we obtain the optimal decay rate $(1+t)^{-\frac{1}{2}}$ for the perturbation in $L^\infty$ norm, which provides a positive answer to the conjecture in \cite{HMX}. We applied the anti-derivative method, introducing the diffusion wave to carry the initial excess mass, diagonalizing the integrated system, and estimating the energy of perturbation in the diagonalized system. Precisely, due to the presence of diffusion waves, the decay rates for errors of perturbed system are too poor to get the optimal decay rate. We find the dissipation structural in the diagonalized system, see \cref{ds}. This observation makes us able to fully utilize the fact that the sign of the derivative of the contact wave is invariant and to control the terms with poor decay rates in energy estimates. For the case with zero mass condition, there are also terms with poor decay rates. In this case, note that there is a cancellation in the linearly degenerate field so that the terms with poor decay rates will not appear in the second equation of the diagonalized system. Thanks to this cancellation and a Poincaré type of estimate obtained by a critical inequality introduced by \cite{HLM}, we get the decay rate of $\ln^{\frac{1}{2}} (2+t)$ for $L^2$ norm of anti-derivatives of perturbation and $(1+t)^{-\frac{1}{2}}\ln^{\frac{1}{2}}(2+t)$ for the $L^2$ norm of perturbation itself, the decay rates are optimal, which is consistent with the results obtained by using pointwise estimate in \cite{XZ} for the system with artificial viscosity. △ Less

Submitted 19 October, 2023; originally announced October 2023.

Convolution quadratures based on block generalized Adams methods

Authors: Ling Liu, Junjie Ma

Abstract: This paper studies a family of convolution quadratures, a numerical technique for efficient evaluation of convolution integrals. We employ the block generalized Adams method to discretize the underlying initial value problem, departing from the well-established approaches that rely on linear multistep formulas or Runge-Kutta methods. The convergence order of the proposed convolution quadrature can… ▽ More This paper studies a family of convolution quadratures, a numerical technique for efficient evaluation of convolution integrals. We employ the block generalized Adams method to discretize the underlying initial value problem, departing from the well-established approaches that rely on linear multistep formulas or Runge-Kutta methods. The convergence order of the proposed convolution quadrature can be dynamically controlled without requiring grid point adjustments, enhancing exibility. Through strategic selection of the local interpolation polynomial and block size, the method achieves high-order convergence for calculation of convolution integrals with hyperbolic kernels. We provide a rigorous convergence analysis for the proposed convolution quadrature and numerically validate our theoretical findings for various convolution integrals. △ Less

Submitted 10 July, 2024; v1 submitted 15 October, 2023; originally announced October 2023.

Comments: 22 pages

MSC Class: 65D32; 65R10

A Long Exact Sequence in Symmetry Breaking: order parameter constraints, defect anomaly-matching, and higher Berry phases

Authors: Arun Debray, Sanath K. Devalapurkar, Cameron Krulewski, Yu Leon Liu, Natalia Pacheco-Tallaj, Ryan Thorngren

Abstract: We study defects in symmetry breaking phases, such as domain walls, vortices, and hedgehogs. In particular, we focus on the localized gapless excitations which sometimes occur at the cores of these objects. These are topologically protected by an 't Hooft anomaly. We classify different symmetry breaking phases in terms of the anomalies of these defects, and relate them to the anomaly of the broken… ▽ More We study defects in symmetry breaking phases, such as domain walls, vortices, and hedgehogs. In particular, we focus on the localized gapless excitations which sometimes occur at the cores of these objects. These are topologically protected by an 't Hooft anomaly. We classify different symmetry breaking phases in terms of the anomalies of these defects, and relate them to the anomaly of the broken symmetry by an anomaly-matching formula. We also derive the obstruction to the existence of a symmetry breaking phase with a local defect. We obtain these results using a long exact sequence of groups of invertible field theories, which we call the "symmetry breaking long exact sequence" (SBLES). The mathematical backbone of the SBLES is studied in a companion paper. Our work further develops the theory of higher Berry phase and its bulk-boundary correspondence, and serves as a new computational tool for classifying symmetry protected topological phases. △ Less

Submitted 8 May, 2024; v1 submitted 28 September, 2023; originally announced September 2023.

Comments: 31 pages, 4 figures. The mathematical sections have been removed and now appear in 2405.04649

On Deterministically Approximating Total Variation Distance

Authors: Weiming Feng, Liqiang Liu, Tianren Liu

Abstract: Total variation distance (TV distance) is an important measure for the difference between two distributions. Recently, there has been progress in approximating the TV distance between product distributions: a deterministic algorithm for a restricted class of product distributions (Bhattacharyya, Gayen, Meel, Myrisiotis, Pavan and Vinodchandran 2023) and a randomized algorithm for general product d… ▽ More Total variation distance (TV distance) is an important measure for the difference between two distributions. Recently, there has been progress in approximating the TV distance between product distributions: a deterministic algorithm for a restricted class of product distributions (Bhattacharyya, Gayen, Meel, Myrisiotis, Pavan and Vinodchandran 2023) and a randomized algorithm for general product distributions (Feng, Guo, Jerrum and Wang 2023). We give a deterministic fully polynomial-time approximation algorithm (FPTAS) for the TV distance between product distributions. Given two product distributions $\mathbb{P}$ and $\mathbb{Q}$ over $[q]^n$, our algorithm approximates their TV distance with relative error $\varepsilon$ in time $O\bigl( \frac{qn^2}{\varepsilon} \log q \log \frac{n}{\varepsilon Δ_{\text{TV}}(\mathbb{P},\mathbb{Q}) } \bigr)$. Our algorithm is built around two key concepts: 1) The likelihood ratio as a distribution, which captures sufficient information to compute the TV distance. 2) We introduce a metric between likelihood ratio distributions, called the minimum total variation distance. Our algorithm computes a sparsified likelihood ratio distribution that is close to the original one w.r.t. the new metric. The approximated TV distance can be computed from the sparsified likelihood ratio. Our technique also implies deterministic FPTAS for the TV distance between Markov chains. △ Less

Submitted 26 September, 2023; originally announced September 2023.

$L_{2,1}$-Norm Regularized Quaternion Matrix Completion Using Sparse Representation and Quaternion QR Decomposition

Authors: Juan Han, Kit Ian Kou, Jifei Miao, Lizhi Liu, Haojiang Li

Abstract: Color image completion is a challenging problem in computer vision, but recent research has shown that quaternion representations of color images perform well in many areas. These representations consider the entire color image and effectively utilize coupling information between the three color channels. Consequently, low-rank quaternion matrix completion (LRQMC) algorithms have gained significan… ▽ More Color image completion is a challenging problem in computer vision, but recent research has shown that quaternion representations of color images perform well in many areas. These representations consider the entire color image and effectively utilize coupling information between the three color channels. Consequently, low-rank quaternion matrix completion (LRQMC) algorithms have gained significant attention. We propose a method based on quaternion Qatar Riyal decomposition (QQR) and quaternion $L_{2,1}$-norm called QLNM-QQR. This new approach reduces computational complexity by avoiding the need to calculate the QSVD of large quaternion matrices. We also present two improvements to the QLNM-QQR method: an enhanced version called IRQLNM-QQR that uses iteratively reweighted quaternion $L_{2,1}$-norm minimization and a method called QLNM-QQR-SR that integrates sparse regularization. Our experiments on natural color images and color medical images show that IRQLNM-QQR outperforms QLNM-QQR and that the proposed QLNM-QQR-SR method is superior to several state-of-the-art methods. △ Less

Submitted 7 September, 2023; originally announced September 2023.

Global existence of spherically symmetry solutions for isothermal Euler-Poisson system outside a ball

Authors: Lingjun Liu

Abstract: In this paper, we consider an isothermal Euler-Poisson system with self-gravitational force, modeling a compact star such as strange quark star. We prove that there exists a global entropy solution with spherically symmetry outside a ball, through the fractional Lax-Friedrichs scheme and the theory of compensated compactness. In this paper, we consider an isothermal Euler-Poisson system with self-gravitational force, modeling a compact star such as strange quark star. We prove that there exists a global entropy solution with spherically symmetry outside a ball, through the fractional Lax-Friedrichs scheme and the theory of compensated compactness. △ Less

Submitted 24 August, 2023; originally announced August 2023.

Kropina metrics with isotropic scalar curvature

Authors: Liulin Liu, Xiaoling Zhang, Lili Zhao

Abstract: In this paper, we study Kropina metrics with isotropic scalar curvature. First, we obtain the expressions of Ricci curvature tensor and scalar curvature. Then, we characterize the Kropina metrics with isotropic scalar curvature on by tensor analysis. In this paper, we study Kropina metrics with isotropic scalar curvature. First, we obtain the expressions of Ricci curvature tensor and scalar curvature. Then, we characterize the Kropina metrics with isotropic scalar curvature on by tensor analysis. △ Less

Submitted 16 August, 2023; originally announced August 2023.

Vanishing viscosity limit to the planar rarefaction wave with vacuum for 3-D full compressible Navier-Stokes equations with temperature-dependent transport coefficients

Authors: Meichen Hou, Lingjun Liu, Shu Wang, Lingda Xu

Abstract: In this paper, we construct a family of global-in-time solutions of the 3-D full compressible Navier-Stokes (N-S) equations with temperature-dependent transport coefficients (including viscosity and heat-conductivity), and show that at arbitrary times {and arbitrary strength} this family of solutions converges to planar rarefaction waves connected to the vacuum as the viscosity vanishes in the sen… ▽ More In this paper, we construct a family of global-in-time solutions of the 3-D full compressible Navier-Stokes (N-S) equations with temperature-dependent transport coefficients (including viscosity and heat-conductivity), and show that at arbitrary times {and arbitrary strength} this family of solutions converges to planar rarefaction waves connected to the vacuum as the viscosity vanishes in the sense of $L^\infty(\R^3)$. We consider the Cauchy problem in $\R^3$ with perturbations of the infinite global norm, particularly, periodic perturbations. To deal with the infinite oscillation, we construct a suitable ansatz carrying this periodic oscillation such that the difference between the solution and the ansatz belongs to some Sobolev space and thus the energy method is feasible. The novelty of this paper is that the viscosity and heat-conductivity are temperature-dependent and degeneracies caused by vacuum. Thus the a priori assumptions and two Gagliardo-Nirenberg type inequalities are essentially used. Next, more careful energy estimates are carried out in this paper, by studying the zero and non-zero modes of the solutions, we obtain not only the convergence rate concerning the viscosity and heat conductivity coefficients but also the exponential time decay rate for the non-zero mode. △ Less

Submitted 23 February, 2024; v1 submitted 6 August, 2023; originally announced August 2023.

Comments: 40 pages. Accepted for publication in Mathematische Annalen

MSC Class: 35Q30; 35Q31; 35Q35; 76N06; 76N10

A symplectic dynamics approach to the spatial isosceles three-body problem

Authors: Xijun Hu, Lei Liu, Yuwei Ou, Pedro A. S. Salomão, Guowei Yu

Abstract: We study the spatial isosceles three-body problem from the perspective of Symplectic Dynamics. For certain choices of mass ratio, angular momentum, and energy, the dynamics on the energy surface is equivalent to a Reeb flow on the tight three-sphere. We find a Hopf link formed by the Euler orbit and a symmetric brake orbit, which spans an open book decomposition whose pages are annulus-like global… ▽ More We study the spatial isosceles three-body problem from the perspective of Symplectic Dynamics. For certain choices of mass ratio, angular momentum, and energy, the dynamics on the energy surface is equivalent to a Reeb flow on the tight three-sphere. We find a Hopf link formed by the Euler orbit and a symmetric brake orbit, which spans an open book decomposition whose pages are annulus-like global surfaces of section. In the case of large mass ratios, the Hopf link is non-resonant, forcing the existence of infinitely many periodic orbits. The rotation number of the Euler orbit plays a fundamental role in the existence of periodic orbits and their symmetries. We explore such symmetries in the Hill region and show that the Euler orbit is negative hyperbolic for an open set of parameters while it can never be positive hyperbolic. Finally, we address convexity and determine for each parameter whether the energy surface is strictly convex, convex, or non-convex. Dynamical consequences of this fact are then discussed. △ Less

Submitted 1 August, 2023; originally announced August 2023.

Comments: 66 pages, 15 figures

Spectral Turán-type problems on sparse spanning graphs

Authors: Lele Liu, Bo Ning

Abstract: Let $F$ be a graph and $\SPEX (n, F)$ be the class of $n$-vertex graphs which attain the maximum spectral radius and contain no $F$ as a subgraph. Let $\EX (n, F)$ be the family of $n$-vertex graphs which contain maximum number of edges and no $F$ as a subgraph. It is a fundamental problem in spectral extremal graph theory to characterize all graphs $F$ such that… ▽ More Let $F$ be a graph and $\SPEX (n, F)$ be the class of $n$-vertex graphs which attain the maximum spectral radius and contain no $F$ as a subgraph. Let $\EX (n, F)$ be the family of $n$-vertex graphs which contain maximum number of edges and no $F$ as a subgraph. It is a fundamental problem in spectral extremal graph theory to characterize all graphs $F$ such that $\SPEX (n, F)\subseteq \EX (n, F)$ when $n$ is sufficiently large. Establishing the conjecture of Cioabă, Desai and Tait [European J. Combin., 2022], Wang, Kang, and Xue [J. Combin. Theory Ser. B, 2023] prove that: for any graph $F$ such that the graphs in $\EX (n, F)$ are Turán graphs plus $O(1)$ edges, $\SPEX (n, F)\subseteq \EX (n, F)$ for sufficiently large $n$. In this paper, we prove that $\SPEX (n, F)\subseteq \EX (n, F)$ for sufficiently large $n$, where $F$ is an $n$-vertex graph with no isolated vertices and $Δ(F) \leq \sqrt{n}/40$. We also prove a signless Laplacian spectral radius version of the above theorem. These results give new contribution to the open problem mentioned above, and can be seen as spectral analogs of a theorem of Alon and Yuster [J. Combin. Theory Ser. B, 2013]. Furthermore, as immediate corollaries, we have tight spectral conditions for the existence of several classes of special graphs, including clique-factors, $k$-th power of Hamilton cycles and $k$-factors in graphs. The first special class of graphs gives a positive answer to a problem of Feng, and the second one extends a previous result of Yan et al. △ Less

Submitted 27 July, 2023; originally announced July 2023.

Quaternion tensor left ring decomposition and application for color image inpainting

Authors: Jifei Miao, Kit Ian Kou, Hongmin Cai, Lizhi Liu

Abstract: In recent years, tensor networks have emerged as powerful tools for solving large-scale optimization problems. One of the most promising tensor networks is the tensor ring (TR) decomposition, which achieves circular dimensional permutation invariance in the model through the utilization of the trace operation and equitable treatment of the latent cores. On the other hand, more recently, quaternion… ▽ More In recent years, tensor networks have emerged as powerful tools for solving large-scale optimization problems. One of the most promising tensor networks is the tensor ring (TR) decomposition, which achieves circular dimensional permutation invariance in the model through the utilization of the trace operation and equitable treatment of the latent cores. On the other hand, more recently, quaternions have gained significant attention and have been widely utilized in color image processing tasks due to their effectiveness in encoding color pixels by considering the three color channels as a unified entity. Therefore, in this paper, based on the left quaternion matrix multiplication, we propose the quaternion tensor left ring (QTLR) decomposition, which inherits the powerful and generalized representation abilities of the TR decomposition while leveraging the advantages of quaternions for color pixel representation. In addition to providing the definition of QTLR decomposition and an algorithm for learning the QTLR format, the paper further proposes a low-rank quaternion tensor completion (LRQTC) model and its algorithm for color image inpainting based on the defined QTLR decomposition. Finally, extensive experiments on color image inpainting demonstrate that the proposed LRQTC method is highly competitive. △ Less

Submitted 16 September, 2023; v1 submitted 20 July, 2023; originally announced July 2023.

Linear Regression on Manifold Structured Data: the Impact of Extrinsic Geometry on Solutions

Authors: Liangchen Liu, Juncai He, Richard Tsai

Abstract: In this paper, we study linear regression applied to data structured on a manifold. We assume that the data manifold is smooth and is embedded in a Euclidean space, and our objective is to reveal the impact of the data manifold's extrinsic geometry on the regression. Specifically, we analyze the impact of the manifold's curvatures (or higher order nonlinearity in the parameterization when the curv… ▽ More In this paper, we study linear regression applied to data structured on a manifold. We assume that the data manifold is smooth and is embedded in a Euclidean space, and our objective is to reveal the impact of the data manifold's extrinsic geometry on the regression. Specifically, we analyze the impact of the manifold's curvatures (or higher order nonlinearity in the parameterization when the curvatures are locally zero) on the uniqueness of the regression solution. Our findings suggest that the corresponding linear regression does not have a unique solution when the embedded submanifold is flat in some dimensions. Otherwise, the manifold's curvature (or higher order nonlinearity in the embedding) may contribute significantly, particularly in the solution associated with the normal directions of the manifold. Our findings thus reveal the role of data manifold geometry in ensuring the stability of regression models for out-of-distribution inferences. △ Less

Submitted 22 July, 2023; v1 submitted 5 July, 2023; originally announced July 2023.

Comments: 13 pages, 6 figures, accepted to TAGML23 workshop of ICML2023, to be published in PMLR

MSC Class: 53Z50 62J05 (Primary) 65D18 68T07 (Secondary) ACM Class: G.1.2; G.4

Large-time asymptotics for degenerate cross-diffusion population models with volume filling

Authors: Xiuqing Chen, Ansgar Jüngel, Xi Lin, Ling Liu

Abstract: The large-time asymptotics of the solutions to a class of degenerate parabolic cross-diffusion systems is analyzed. The equations model the interaction of an arbitrary number of population species in a bounded domain with no-flux boundary conditions. Compared to previous works, we allow for different diffusivities and degenerate nonlinearities. The proof is based on the relative entropy method, bu… ▽ More The large-time asymptotics of the solutions to a class of degenerate parabolic cross-diffusion systems is analyzed. The equations model the interaction of an arbitrary number of population species in a bounded domain with no-flux boundary conditions. Compared to previous works, we allow for different diffusivities and degenerate nonlinearities. The proof is based on the relative entropy method, but in contrast to usual arguments, the relative entropy and entropy production are not directly related by a logarithmic Sobolev inequality. The key idea is to apply convex Sobolev inequalities to modified entropy densities including "iterated degenerate" functions. △ Less

Submitted 30 June, 2023; originally announced June 2023.

MSC Class: 35K51; 35K59; 35K65; 35Q92; 92D25

QoS-based Beamforming and Compression Design for Cooperative Cellular Networks via Lagrangian Duality

Authors: Xilai Fan, Ya-Feng Liu, Liang Liu, Tsung-Hui Chang

Abstract: This paper considers the quality-of-service (QoS)-based joint beamforming and compression design problem in the downlink cooperative cellular network, where multiple relay-like base stations (BSs), connected to the central processor via rate-limited fronthaul links, cooperatively transmit messages to the users. The problem of interest is formulated as the minimization of the total transmit power o… ▽ More This paper considers the quality-of-service (QoS)-based joint beamforming and compression design problem in the downlink cooperative cellular network, where multiple relay-like base stations (BSs), connected to the central processor via rate-limited fronthaul links, cooperatively transmit messages to the users. The problem of interest is formulated as the minimization of the total transmit power of the BSs, subject to all users' signal-to-interference-plus-noise ratio (SINR) constraints and all BSs' fronthaul rate constraints. In this paper, we first show that there is no duality gap between the considered joint optimization problem and its Lagrangian dual by showing the tightness of its semidefinite relaxation (SDR). Then, we propose an efficient algorithm based on the above duality result for solving the considered problem. The proposed algorithm judiciously exploits the special structure of an enhanced Karush-Kuhn-Tucker (KKT) conditions of the considered problem and finds the solution that satisfies the enhanced KKT conditions via two fixed point iterations. Two key features of the proposed algorithm are: (1) it is able to detect whether the considered problem is feasible or not and find its globally optimal solution when it is feasible; (2) it is highly efficient because both of the fixed point iterations in the proposed algorithm are linearly convergent and evaluating the functions in the fixed point iterations are computationally cheap. Numerical results show the global optimality and efficiency of the proposed algorithm. △ Less

Submitted 24 June, 2023; originally announced June 2023.

Comments: 15 pages, 7 figures, submitted for possible publication

Assumption-lean falsification tests of rate double-robustness of double-machine-learning estimators

Authors: Lin Liu, Rajarshi Mukherjee, James M. Robins

Abstract: The class of doubly-robust (DR) functionals studied by Rotnitzky et al. (2021) is of central importance in economics and biostatistics. It strictly includes both (i) the class of mean-square continuous functionals that can be written as an expectation of an affine functional of a conditional expectation studied by Chernozhukov et al. (2022b) and (ii) the class of functionals studied by Robins et a… ▽ More The class of doubly-robust (DR) functionals studied by Rotnitzky et al. (2021) is of central importance in economics and biostatistics. It strictly includes both (i) the class of mean-square continuous functionals that can be written as an expectation of an affine functional of a conditional expectation studied by Chernozhukov et al. (2022b) and (ii) the class of functionals studied by Robins et al. (2008). The present state-of-the-art estimators for DR functionals $ψ$ are double-machine-learning (DML) estimators (Chernozhukov et al., 2018). A DML estimator $\widehatψ_{1}$ of $ψ$ depends on estimates $\widehat{p} (x)$ and $\widehat{b} (x)$ of a pair of nuisance functions $p(x)$ and $b(x)$, and is said to satisfy "rate double-robustness" if the Cauchy--Schwarz upper bound of its bias is $o (n^{- 1/2})$. Were it achievable, our scientific goal would have been to construct valid, assumption-lean (i.e. no complexity-reducing assumptions on $b$ or $p$) tests of the validity of a nominal $(1 - α)$ Wald confidence interval (CI) centered at $\widehatψ_{1}$. But this would require a test of the bias to be $o (n^{-1/2})$, which can be shown not to exist. We therefore adopt the less ambitious goal of falsifying, when possible, an analyst's justification for her claim that the reported $(1 - α)$ Wald CI is valid. In many instances, an analyst justifies her claim by imposing complexity-reducing assumptions on $b$ and $p$ to ensure "rate double-robustness". Here we exhibit valid, assumption-lean tests of $H_{0}$: "rate double-robustness holds", with non-trivial power against certain alternatives. If $H_{0}$ is rejected, we will have falsified her justification. However, no assumption-lean test of $H_{0}$, including ours, can be a consistent test. Thus, the failure of our test to reject is not meaningful evidence in favor of $H_{0}$. △ Less

Submitted 28 August, 2023; v1 submitted 18 June, 2023; originally announced June 2023.

Comments: corrected several extra typos and references

The dissipativity and controllability of input affine systems

Authors: Qin Xu, Liu Liu, Yufeng Lu

Abstract: The paper is considered with the dissipative theory and feedback control under the framework of dissipation with the supply rate is the inner product of input u and the derivation of output y for the linear/nonlinear time-invariant input affine system. A necessary and sufficient condition of dissipativity of the class of systems is presented. Furthermore, we obtain necessary and sufficient conditi… ▽ More The paper is considered with the dissipative theory and feedback control under the framework of dissipation with the supply rate is the inner product of input u and the derivation of output y for the linear/nonlinear time-invariant input affine system. A necessary and sufficient condition of dissipativity of the class of systems is presented. Furthermore, we obtain necessary and sufficient conditions of the negative imaginary property, input strictly negative imaginary property and output strictly negative imaginary property for the nonlinear input affine system, which can cover the case that the linear time-invariant system in the literature. Besides, under given conditions, we consider the equivalence between the negative imaginary property and input-output stability. In other words, the system with negative imaginary property may be input-output stable under some certain conditions; conversely, the input-output stable system may also have negative imaginary property. Last but not least, we consider the equivalence between the asymptotical stability and strict dissipativity with the supply rate and some dissipative matrices are given for a certain linear time-invariant system. △ Less

Submitted 14 June, 2023; originally announced June 2023.

The qualitative behavior at a vortex point for the Chern-Simon-Higgs equation

Authors: Jiayu Li, Lei Liu

Abstract: In this paper, we study the qualitative behavior at a vortex blow-up point for Chern-Simon-Higgs equation. Roughly speaking, we will establish an energy identity at a each such point, i.e. the local mass is the sum of the bubbles. Moreover, we prove that either there is only one bubble which is a singular bubble or there are more than two bubbles which contains no singular bubble. Meanwhile, we pr… ▽ More In this paper, we study the qualitative behavior at a vortex blow-up point for Chern-Simon-Higgs equation. Roughly speaking, we will establish an energy identity at a each such point, i.e. the local mass is the sum of the bubbles. Moreover, we prove that either there is only one bubble which is a singular bubble or there are more than two bubbles which contains no singular bubble. Meanwhile, we prove that the energies of these bubbles must satisfy a quadratic polynomial which can be used to prove the simple blow-up property when the multiplicity is small. As is well known, for many Liouville type system, Pohozaev type identity is a quadratic polynomial corresponding to energies which can be used directly to compute the local mass at a blow-up point. The difficulty here is that, besides the energy's integration, there is a additional term in the Pohozaev type identity of Chern-Simon-Higgs equation. We need some more detailed and delicated analysis to deal with it. △ Less

Submitted 6 June, 2023; originally announced June 2023.

A unified Bayesian inversion approach for a class of tumor growth models with different pressure laws

Authors: Yu Feng, Liu Liu, Zhennan Zhou

Abstract: In this paper, we use the Bayesian inversion approach to study the data assimilation problem for a family of tumor growth models described by porous-medium type equations. The models contain uncertain parameters and are indexed by a physical parameter $m$, which characterizes the constitutive relation between density and pressure. Based on these models, we employ the Bayesian inversion framework t… ▽ More In this paper, we use the Bayesian inversion approach to study the data assimilation problem for a family of tumor growth models described by porous-medium type equations. The models contain uncertain parameters and are indexed by a physical parameter $m$, which characterizes the constitutive relation between density and pressure. Based on these models, we employ the Bayesian inversion framework to infer parametric and nonparametric unknowns that affect tumor growth from noisy observations of tumor cell density. We establish the well-posedness and the stability theories for the Bayesian inversion problem and further prove the convergence of the posterior distribution in the so-called incompressible limit, $m \rightarrow \infty$. Since the posterior distribution across the index regime $m\in[2,\infty)$ can thus be treated in a unified manner, such theoretical results also guide the design of the numerical inference for the unknown. We propose a generic computational framework for such inverse problems, which consists of a typical sampling algorithm and an asymptotic preserving solver for the forward problem. With extensive numerical tests, we demonstrate that the proposed method achieves satisfactory accuracy in the Bayesian inference of the tumor growth models, which is uniform with respect to the constitutive relation. △ Less

Submitted 13 February, 2024; v1 submitted 3 June, 2023; originally announced June 2023.

Comments: 30 pages, 17 figures

MSC Class: 35R30 62F15 65M32 92-10

Coding information into all infinite subsets of a dense set

Authors: Matthew Harrison-Trainor, Lu Liu, Patrick Lutz

Abstract: Suppose you have an uncomputable set $X$ and you want to find a set $A$, all of whose infinite subsets compute $X$. There are several ways to do this, but all of them seem to produce a set $A$ which is fairly sparse. We show that this is necessary in the following technical sense: if $X$ is uncomputable and $A$ is a set of positive lower density then $A$ has an infinite subset which does not compu… ▽ More Suppose you have an uncomputable set $X$ and you want to find a set $A$, all of whose infinite subsets compute $X$. There are several ways to do this, but all of them seem to produce a set $A$ which is fairly sparse. We show that this is necessary in the following technical sense: if $X$ is uncomputable and $A$ is a set of positive lower density then $A$ has an infinite subset which does not compute $X$. We also prove an analogous result for PA degree: if $X$ is uncomputable and $A$ is a set of positive lower density then $A$ has an infinite subset which is not of PA degree. We will show that these theorems are sharp in certain senses and also prove a quantitative version formulated in terms of Kolmogorov complexity. Our results use a modified version of Mathias forcing and build on work by Seetapun, Liu, and others on the reverse math of Ramsey's theorem for pairs. △ Less

Submitted 11 August, 2023; v1 submitted 1 June, 2023; originally announced June 2023.

Comments: 37 pages. Paper updated to fix typos

MSC Class: 03D28; 03D32; 68Q30

On a neural network approach for solving potential control problem of the semiclassical Schrödinger equation

Authors: Yating Wang, Liu Liu

Abstract: Robust control design for quantum systems is a challenging and key task for practical technology. In this work, we apply neural networks to learn the control problem for the semiclassical Schrödinger equation, where the control variable is the potential given by an external field that may contain uncertainties. Inspired by a relevant work [29], we incorporate the sampling-based learning process in… ▽ More Robust control design for quantum systems is a challenging and key task for practical technology. In this work, we apply neural networks to learn the control problem for the semiclassical Schrödinger equation, where the control variable is the potential given by an external field that may contain uncertainties. Inspired by a relevant work [29], we incorporate the sampling-based learning process into the training of networks, while combining with the fast time-splitting spectral method for the Schrödinger equation in the semiclassical regime. The numerical results have shown the efficiency and accuracy of our proposed deep learning approach. △ Less

Submitted 30 May, 2023; originally announced May 2023.

Predict-then-Calibrate: A New Perspective of Robust Contextual LP

Authors: Chunlin Sun, Linyu Liu, Xiaocheng Li

Abstract: Contextual optimization, also known as predict-then-optimize or prescriptive analytics, considers an optimization problem with the presence of covariates (context or side information). The goal is to learn a prediction model (from the training data) that predicts the objective function from the covariates, and then in the test phase, solve the optimization problem with the covariates but without t… ▽ More Contextual optimization, also known as predict-then-optimize or prescriptive analytics, considers an optimization problem with the presence of covariates (context or side information). The goal is to learn a prediction model (from the training data) that predicts the objective function from the covariates, and then in the test phase, solve the optimization problem with the covariates but without the observation of the objective function. In this paper, we consider a risk-sensitive version of the problem and propose a generic algorithm design paradigm called predict-then-calibrate. The idea is to first develop a prediction model without concern for the downstream risk profile or robustness guarantee, and then utilize calibration (or recalibration) methods to quantify the uncertainty of the prediction. While the existing methods suffer from either a restricted choice of the prediction model or strong assumptions on the underlying data, we show the disentangling of the prediction model and the calibration/uncertainty quantification has several advantages. First, it imposes no restriction on the prediction model and thus fully unleashes the potential of off-the-shelf machine learning methods. Second, the derivation of the risk and robustness guarantee can be made independent of the choice of the prediction model through a data-splitting idea. Third, our paradigm of predict-then-calibrate applies to both (risk-sensitive) robust and (risk-neutral) distributionally robust optimization (DRO) formulations. Theoretically, it gives new generalization bounds for the contextual LP problem and sheds light on the existing results of DRO for contextual LP. Numerical experiments further reinforce the advantage of the predict-then-calibrate paradigm in that an improvement on either the prediction model or the calibration model will lead to a better final performance. △ Less

Submitted 10 May, 2024; v1 submitted 24 May, 2023; originally announced May 2023.

Comments: 30 pages, 8 figures

MSC Class: 90C15 (Primary); 90C05 (Secondary)

Unsolved Problems in Spectral Graph Theory

Authors: Lele Liu, Bo Ning

Abstract: Spectral graph theory is a captivating area of graph theory that employs the eigenvalues and eigenvectors of matrices associated with graphs to study them. In this paper, we present a collection of $20$ topics in spectral graph theory, covering a range of open problems and conjectures. Our focus is primarily on the adjacency matrix of graphs, and for each topic, we provide a brief historical overv… ▽ More Spectral graph theory is a captivating area of graph theory that employs the eigenvalues and eigenvectors of matrices associated with graphs to study them. In this paper, we present a collection of $20$ topics in spectral graph theory, covering a range of open problems and conjectures. Our focus is primarily on the adjacency matrix of graphs, and for each topic, we provide a brief historical overview. △ Less

Submitted 30 July, 2023; v1 submitted 17 May, 2023; originally announced May 2023.

Comments: v3, 30 pages, 1 figure, include comments from Clive Elphick, Xiaofeng Gu, William Linz, and Dragan Stevanović, respectively. Thanks! This paper will be published in Operations Research Transactions

Evans-Lewis inequalities in dimension three

Authors: Yi C. Huang, Li Liu

Abstract: In proving Rellich inequalities in the framework of equalities, N. Bez, S. Machihara, and T. Ozawa obtained some interesting norm inequalities in the spirit of Evans and Lewis that compare the standard Laplacian with its radial and spherical components. In this paper we give a simple unified proof and a strict improvement of these Evans-Lewis inequalities in the subtle dimension three case. Our ap… ▽ More In proving Rellich inequalities in the framework of equalities, N. Bez, S. Machihara, and T. Ozawa obtained some interesting norm inequalities in the spirit of Evans and Lewis that compare the standard Laplacian with its radial and spherical components. In this paper we give a simple unified proof and a strict improvement of these Evans-Lewis inequalities in the subtle dimension three case. Our approach is robust and explains clearly the occurrence of the sharp constant. △ Less

Submitted 23 May, 2023; v1 submitted 11 May, 2023; originally announced May 2023.

Comments: added a section which extends our machinery to the abstract framework; added more references

MSC Class: Primary 26D10; Secondary 46E35; 35A23

Root-n consistent semiparametric learning with high-dimensional nuisance functions under minimal sparsity

Authors: Lin Liu, Yuhao Wang

Abstract: Treatment effect estimation under unconfoundedness is a fundamental task in causal inference. In response to the challenge of analyzing high-dimensional datasets collected in substantive fields such as epidemiology, genetics, economics, and social sciences, many methods for treatment effect estimation with high-dimensional nuisance parameters (the outcome regression and the propensity score) have… ▽ More Treatment effect estimation under unconfoundedness is a fundamental task in causal inference. In response to the challenge of analyzing high-dimensional datasets collected in substantive fields such as epidemiology, genetics, economics, and social sciences, many methods for treatment effect estimation with high-dimensional nuisance parameters (the outcome regression and the propensity score) have been developed in recent years. However, it is still unclear what is the necessary and sufficient sparsity condition on the nuisance parameters for the treatment effect to be $\sqrt{n}$-estimable. In this paper, we propose a new Double-Calibration strategy that corrects the estimation bias of the nuisance parameter estimates computed by regularized high-dimensional techniques and demonstrate that the corresponding Doubly-Calibrated estimator achieves $1 / \sqrt{n}$-rate as long as one of the nuisance parameters is sparse with sparsity below $\sqrt{n} / \log p$, where $p$ denotes the ambient dimension of the covariates, whereas the other nuisance parameter can be arbitrarily complex and completely misspecified. The Double-Calibration strategy can also be applied to settings other than treatment effect estimation, e.g. regression coefficient estimation in the presence of diverging number of controls in a semiparametric partially linear model. △ Less

Submitted 12 April, 2024; v1 submitted 6 May, 2023; originally announced May 2023.

Comments: This is a preliminary draft intended as a reference for Yuhao Wang's presentation at the Online Causal Inference Seminar. The full article is currently writing in progress. In its first arXiv version, the proof of the lower bound was found to be incorrect. Consequently, all claims regarding the lower bound have been removed in this revised version. All other results remain unchanged

On the well-posedness of tracking Dirichlet data for Bernoulli free boundary problems

Authors: Wei Gong, Le Liu

Abstract: The aim of this paper is to study the shape optimization method for solving the Bernoulli free boundary problem, a well-known ill-posed problem that seeks the unknown free boundary through Cauchy data. Different formulations have been proposed in the literature that differ in the choice of the objective functional. Specifically, it was shown respectively in [14] and [16] that tracking Neumann data… ▽ More The aim of this paper is to study the shape optimization method for solving the Bernoulli free boundary problem, a well-known ill-posed problem that seeks the unknown free boundary through Cauchy data. Different formulations have been proposed in the literature that differ in the choice of the objective functional. Specifically, it was shown respectively in [14] and [16] that tracking Neumann data is well-posed but tracking Dirichlet data is not. In this paper we propose a new well-posed objective functional that tracks Dirichlet data at the free boundary. By calculating the Euler derivative and the shape Hessian of the objective functional we show that the new formulation is well-posed, i.e., the shape Hessian is coercive at the minimizers. The coercivity of the shape Hessian may ensure the existence of optimal solutions for the nonlinear Ritz-Galerkin approximation method and its convergence, thus is crucial for the formulation. As a summary, we conclude that tracking Dirichlet or Neumann data in its energy norm is not sufficient, but tracking it in a half an order higher norm will be well-posed. To support our theoretical results we carry out extensive numerical experiments. △ Less

Submitted 3 May, 2023; originally announced May 2023.

Comments: 27 pages, 6 figures, 4 tables, 27 references, will be resubmitted