-
Accelerating Eigenvalue Computation for Nuclear Structure Calculations via Perturbative Corrections
Authors:
Dong Min Roh,
Esmond Ng,
Chao Yang,
Dean Lee,
Pieter Maris,
James P. Vary
Abstract:
We present a new method for computing the lowest few eigenvalues and the corresponding eigenvectors of a nuclear many-body Hamiltonian represented in a truncated configuration interaction subspace, i.e., the no-core shell model (NCSM). The method uses the hierarchical structure of the NCSM Hamiltonian to partition the Hamiltonian as the sum of two matrices. The first matrix corresponds to the Hami…
▽ More
We present a new method for computing the lowest few eigenvalues and the corresponding eigenvectors of a nuclear many-body Hamiltonian represented in a truncated configuration interaction subspace, i.e., the no-core shell model (NCSM). The method uses the hierarchical structure of the NCSM Hamiltonian to partition the Hamiltonian as the sum of two matrices. The first matrix corresponds to the Hamiltonian represented in a small configuration space, whereas the second is viewed as the perturbation to the first matrix. Eigenvalues and eigenvectors of the first matrix can be computed efficiently. Perturbative corrections to the eigenvectors of the first matrix can be obtained from the solutions of a sequence of linear systems of equations defined in the small configuration space. These correction vectors can be combined with the approximate eigenvectors of the first matrix to construct a subspace from which more accurate approximations of the desired eigenpairs can be obtained. We call this method a Subspace Projection with Perturbative Corrections (SPPC) method. We show by numerical examples that the SPPC method can be more efficient than conventional iterative methods for solving large-scale eigenvalue problems such as the Lanczos, block Lanczos and the locally optimal block preconditioned conjugate gradient (LOBPCG) method. The method can also be combined with other methods to avoid convergence stagnation.
△ Less
Submitted 11 July, 2024;
originally announced July 2024.
-
Solving General Natural-Language-Description Optimization Problems with Large Language Models
Authors:
Jihai Zhang,
Wei Wang,
Siyan Guo,
Li Wang,
Fangquan Lin,
Cheng Yang,
Wotao Yin
Abstract:
Optimization problems seek to find the best solution to an objective under a set of constraints, and have been widely investigated in real-world applications. Modeling and solving optimization problems in a specific domain typically require a combination of domain knowledge, mathematical skills, and programming ability, making it difficult for general users and even domain professionals. In this p…
▽ More
Optimization problems seek to find the best solution to an objective under a set of constraints, and have been widely investigated in real-world applications. Modeling and solving optimization problems in a specific domain typically require a combination of domain knowledge, mathematical skills, and programming ability, making it difficult for general users and even domain professionals. In this paper, we propose a novel framework called OptLLM that augments LLMs with external solvers. Specifically, OptLLM accepts user queries in natural language, convert them into mathematical formulations and programming codes, and calls the solvers to calculate the results for decision-making. In addition, OptLLM supports multi-round dialogues to gradually refine the modeling and solving of optimization problems. To illustrate the effectiveness of OptLLM, we provide tutorials on three typical optimization applications and conduct experiments on both prompt-based GPT models and a fine-tuned Qwen model using a large-scale selfdeveloped optimization dataset. Experimental results show that OptLLM works with various LLMs, and the fine-tuned model achieves an accuracy boost compared to the promptbased models. Some features of OptLLM framework have been available for trial since June 2023 (https://opt.alibabacloud.com/chat or https://opt.aliyun.com/chat).
△ Less
Submitted 9 July, 2024;
originally announced July 2024.
-
On the boundedness of degenerate hypergraphs
Authors:
Jianfeng Hou,
Caiyun Hu,
Heng Li,
Xizhi Liu,
Caihong Yang,
Yixiao Zhang
Abstract:
We investigate the impact of a high-degree vertex in Turán problems for degenerate hypergraphs (including graphs). We say an $r$-graph $F$ is bounded if there exist constants $α, β>0$ such that for large $n$, every $n$-vertex $F$-free $r$-graph with a vertex of degree at least $α\binom{n-1}{r-1}$ has fewer than $(1-β) \cdot \mathrm{ex}(n,F)$ edges. The boundedness property is crucial for recent wo…
▽ More
We investigate the impact of a high-degree vertex in Turán problems for degenerate hypergraphs (including graphs). We say an $r$-graph $F$ is bounded if there exist constants $α, β>0$ such that for large $n$, every $n$-vertex $F$-free $r$-graph with a vertex of degree at least $α\binom{n-1}{r-1}$ has fewer than $(1-β) \cdot \mathrm{ex}(n,F)$ edges. The boundedness property is crucial for recent works~\cite{HHLLYZ23a,DHLY24} that aim to extend the classical Hajnal--Szemerédi Theorem and the anti-Ramsey theorems of Erdős--Simonovits--Sós.
We show that many well-studied degenerate hypergraphs, such as all even cycles, most complete bipartite graphs, and the expansion of most complete bipartite graphs, are bounded. In addition, to prove the boundedness of the expansion of complete bipartite graphs, we introduce and solve a Zarankiewicz-type problem for $3$-graphs, strengthening a theorem by Kostochka--Mubayi--Verstraëte~\cite{KMV15}.
△ Less
Submitted 29 June, 2024;
originally announced July 2024.
-
Tight bounds for rainbow partial $F$-tiling in edge-colored complete hypergraphs
Authors:
Jinghua Deng,
Jianfeng Hou,
Xizhi Liu,
Caihong Yang
Abstract:
For an $r$-graph $F$ and integers $n,t$ satisfying $t \le n/v(F)$, let $\mathrm{ar}(n,tF)$ denote the minimum integer $N$ such that every edge-coloring of $K_{n}^{r}$ using $N$ colors contains a rainbow copy of $tF$, where $tF$ is the $r$-graphs consisting of $t$ vertex-disjoint copies of $F$. The case $t=1$ is the classical anti-Ramsey problem proposed by Erdős--Simonovits--Sós~\cite{ESS75}. When…
▽ More
For an $r$-graph $F$ and integers $n,t$ satisfying $t \le n/v(F)$, let $\mathrm{ar}(n,tF)$ denote the minimum integer $N$ such that every edge-coloring of $K_{n}^{r}$ using $N$ colors contains a rainbow copy of $tF$, where $tF$ is the $r$-graphs consisting of $t$ vertex-disjoint copies of $F$. The case $t=1$ is the classical anti-Ramsey problem proposed by Erdős--Simonovits--Sós~\cite{ESS75}. When $F$ is a single edge, this becomes the rainbow matching problem introduced by Schiermeyer~\cite{Sch04} and Özkahya--Young~\cite{OY13}. We conduct a systematic study of $\mathrm{ar}(n,tF)$ for the case where $t$ is much smaller than $\mathrm{ex}(n,F)/n^{r-1}$. Our first main result provides a reduction of $\mathrm{ar}(n,tF)$ to $\mathrm{ar}(n,2F)$ when $F$ is bounded and smooth, two properties satisfied by most previously studied hypergraphs. Complementing the first result, the second main result, which utilizes gaps between Turán numbers, determines $\mathrm{ar}(n,tF)$ for relatively smaller $t$. Together, these two results determine $\mathrm{ar}(n,tF)$ for a large class of hypergraphs. Additionally, the latter result has the advantage of being applicable to hypergraphs with unknown Turán densities, such as the famous tetrahedron $K_{4}^{3}$.
△ Less
Submitted 21 June, 2024; v1 submitted 20 June, 2024;
originally announced June 2024.
-
CPAFT: A Consistent Parallel Advancing Front Technique for Unstructured Triangular/Tetrahedral Mesh Generation
Authors:
Chengdi Ma,
Jizu Huang,
Hao Luo,
Chao Yang
Abstract:
Compared with the remarkable progress made in parallel numerical solvers of partial differential equations,the development of algorithms for generating unstructured triangular/tetrahedral meshes has been relatively sluggish. In this paper, we propose a novel, consistent parallel advancing front technique (CPAFT) by combining the advancing front technique, the domain decomposition method based on s…
▽ More
Compared with the remarkable progress made in parallel numerical solvers of partial differential equations,the development of algorithms for generating unstructured triangular/tetrahedral meshes has been relatively sluggish. In this paper, we propose a novel, consistent parallel advancing front technique (CPAFT) by combining the advancing front technique, the domain decomposition method based on space-filling curves, the distributed forest-of-overlapping-trees approach, and the consistent parallel maximal independent set algorithm. The newly proposed CPAFT algorithm can mathematically ensure that the generated unstructured triangular/tetrahedral meshes are independent of the number of processors and the implementation of domain decomposition. Several numerical tests are conducted to validate the parallel consistency and outstanding parallel efficiency of the proposed algorithm, which scales effectively up to two thousand processors. This is, as far as we know, the first parallel unstructured triangular/tetrahedral mesh generator with scalability to O(1,000) CPU processors.
△ Less
Submitted 31 May, 2024;
originally announced May 2024.
-
A simple inverse power method for balanced graph cut
Authors:
Sihong Shao,
Chuan Yang
Abstract:
The existing inverse power ($\mathbf{IP}$) method for solving the balanced graph cut lacks local convergence and its inner subproblem requires a nonsmooth convex solver. To address these issues, we develop a simple inverse power ($\mathbf{SIP}$) method using a novel equivalent continuous formulation of the balanced graph cut, and its inner subproblem allows an explicit analytic solution, which is…
▽ More
The existing inverse power ($\mathbf{IP}$) method for solving the balanced graph cut lacks local convergence and its inner subproblem requires a nonsmooth convex solver. To address these issues, we develop a simple inverse power ($\mathbf{SIP}$) method using a novel equivalent continuous formulation of the balanced graph cut, and its inner subproblem allows an explicit analytic solution, which is the biggest advantage over $\mathbf{IP}$ and constitutes the main reason why we call it $\mathit{simple}$. By fully exploiting the closed-form of the inner subproblem solution, we design a boundary-detected subgradient selection with which $\mathbf{SIP}$ is proved to be locally converged. We show that $\mathbf{SIP}$ is also applicable to a new ternary valued $θ$-balanced cut which reduces to the balanced cut when $θ=1$. When $\mathbf{SIP}$ reaches its local optimum, we seamlessly transfer to solve the $θ$-balanced cut within exactly the same iteration algorithm framework and thus obtain $\mathbf{SIP}$-$\mathbf{perturb}$ -- an efficient local breakout improvement of $\mathbf{SIP}$, which transforms some ``partitioned" vertices back to the ``un-partitioned" ones through the adjustable $θ$. Numerical experiments on G-set for Cheeger cut and Sparsest cut demonstrate that $\mathbf{SIP}$ is significantly faster than $\mathbf{IP}$ while maintaining approximate solutions of comparable quality, and $\mathbf{SIP}$-$\mathbf{perturb}$ outperforms $\mathtt{Gurobi}$ in terms of both computational cost and solution quality.
△ Less
Submitted 28 May, 2024;
originally announced May 2024.
-
Multi-qubit Lattice Surgery Scheduling
Authors:
Allyson Silva,
Xiangyi Zhang,
Zak Webb,
Mia Kramer,
Chan Woo Yang,
Xiao Liu,
Jessica Lemieux,
Ka-Wai Chen,
Artur Scherer,
Pooya Ronagh
Abstract:
Fault-tolerant quantum computation using two-dimensional topological quantum error correcting codes can benefit from multi-qubit long-range operations. By using simple commutation rules, a quantum circuit can be transpiled into a sequence of solely non-Clifford multi-qubit gates. Prior work on fault-tolerant compilation avoids optimal scheduling of such gates since they reduce the parallelizabilit…
▽ More
Fault-tolerant quantum computation using two-dimensional topological quantum error correcting codes can benefit from multi-qubit long-range operations. By using simple commutation rules, a quantum circuit can be transpiled into a sequence of solely non-Clifford multi-qubit gates. Prior work on fault-tolerant compilation avoids optimal scheduling of such gates since they reduce the parallelizability of the circuit. We observe that the reduced parallelization potential is outweighed by the significant reduction in the number of gates. We therefore devise a method for scheduling multi-qubit lattice surgery using an earliest-available-first policy, solving the associated forest packing problem using a representation of the multi-qubit gates as Steiner trees. Our extensive testing on random and application-inspired circuits demonstrates the method's scalability and performance. We show that the transpilation significantly reduces the circuit length on the set of circuits tested, and that the resulting circuit of multi-qubit gates has a further reduction in the expected circuit execution time compared to serial execution.
△ Less
Submitted 10 June, 2024; v1 submitted 27 May, 2024;
originally announced May 2024.
-
APTT: An accuracy-preserved tensor-train method for the Boltzmann-BGK equation
Authors:
Zhitao Zhu,
Chuanfu Xiao,
Kejun Tang,
Jizu Huang,
Chao Yang
Abstract:
Solving the Boltzmann-BGK equation with traditional numerical methods suffers from high computational and memory costs due to the curse of dimensionality. In this paper, we propose a novel accuracy-preserved tensor-train (APTT) method to efficiently solve the Boltzmann-BGK equation. A second-order finite difference scheme is applied to discretize the Boltzmann-BGK equation, resulting in a tensor a…
▽ More
Solving the Boltzmann-BGK equation with traditional numerical methods suffers from high computational and memory costs due to the curse of dimensionality. In this paper, we propose a novel accuracy-preserved tensor-train (APTT) method to efficiently solve the Boltzmann-BGK equation. A second-order finite difference scheme is applied to discretize the Boltzmann-BGK equation, resulting in a tensor algebraic system at each time step. Based on the low-rank TT representation, the tensor algebraic system is then approximated as a TT-based low-rank system, which is efficiently solved using the TT-modified alternating least-squares (TT-MALS) solver. Thanks to the low-rank TT representation, the APTT method can significantly reduce the computational and memory costs compared to traditional numerical methods. Theoretical analysis demonstrates that the APTT method maintains the same convergence rate as that of the finite difference scheme. The convergence rate and efficiency of the APTT method are validated by several benchmark test cases.
△ Less
Submitted 21 May, 2024;
originally announced May 2024.
-
Extremal oriented graphs avoiding 1-subdivision of an in-star
Authors:
Zejun Huang,
Chenxi Yang
Abstract:
An oriented graph is a digraph obtained from an undirected graph by choosing an orientation for each edge. Given a positive integer $n$ and an oriented graph $F$, the oriented Tur$\acute{\rm a}$n number $ex_{ori}(n,F)$ is the maximum number of arcs in an $F$-free oriented graph of order $n$. In this paper, we investigate the oriented Tur$\acute{\rm a}$n number…
▽ More
An oriented graph is a digraph obtained from an undirected graph by choosing an orientation for each edge. Given a positive integer $n$ and an oriented graph $F$, the oriented Tur$\acute{\rm a}$n number $ex_{ori}(n,F)$ is the maximum number of arcs in an $F$-free oriented graph of order $n$. In this paper, we investigate the oriented Tur$\acute{\rm a}$n number $ex_{ori}(n, \overrightarrow{S_{k,1}} )$, where $\overrightarrow{S_{k,1}}$ is the $1$-subdivision of the in-star of order $k+1$. We determine $ex_{ori}(n,\overrightarrow{S_{k,1}}) $ for $k=2,3$ as well as the extremal oriented graphs. For $k\ge 4$, we establish a lower bound and an upper bound on $ex_{ori}(n,\overrightarrow{S_{k,1}})$.
△ Less
Submitted 20 May, 2024;
originally announced May 2024.
-
NP-completeness of Tiling Finite Simply Connected Regions with a Fixed Set of Wang Tiles
Authors:
Chao Yang,
Zhujun Zhang
Abstract:
The computational complexity of tiling finite simply connected regions with a fixed set of tiles is studied in this paper. We show that the problem of tiling simply connected regions with a fixed set of $23$ Wang tiles is NP-complete. As a consequence, the problem of tiling simply connected regions with a fixed set of $111$ rectangles is NP-complete. Our results improve that of Igor Pak and Jed Ya…
▽ More
The computational complexity of tiling finite simply connected regions with a fixed set of tiles is studied in this paper. We show that the problem of tiling simply connected regions with a fixed set of $23$ Wang tiles is NP-complete. As a consequence, the problem of tiling simply connected regions with a fixed set of $111$ rectangles is NP-complete. Our results improve that of Igor Pak and Jed Yang by using fewer numbers of tiles. Notably in the case of Wang tiles, the number has decreased by more than one third from $35$ to $23$.
△ Less
Submitted 2 May, 2024;
originally announced May 2024.
-
Critical grid method: An extensible Smoothed Particle Hydrodynamics fluid general interpolation method for Fluid-Structure Interaction surface coupling based on preCICE
Authors:
Sifan Long,
Xiaowei Guo,
Xiaokang Fan,
Canqun Yang
Abstract:
Solving Fluid-Structure Interaction (FSI) problems using traditional methods is a big challenge in the field of numerical simulation. As a powerful multi-physical field coupled library, preCICE has a bright application prospect for solving FSI, which supports many open/closed source software and commercial CFD solvers to solve FSI problems in the form of a black box. However, this library currentl…
▽ More
Solving Fluid-Structure Interaction (FSI) problems using traditional methods is a big challenge in the field of numerical simulation. As a powerful multi-physical field coupled library, preCICE has a bright application prospect for solving FSI, which supports many open/closed source software and commercial CFD solvers to solve FSI problems in the form of a black box. However, this library currently only supports mesh-based coupling schemes. This paper proposes a critical grid (mesh) as an intermediate medium for the particle method to connect a bidirectional coupling tool named preCICE. The particle and critical mesh are used to interpolate the displacement and force so that the pure Lagrangian Smoothed Particle Hydrodynamic (SPH) method can also solve the FSI problem. This method is called the particle mesh coupling (PMC) method, which theoretically solves the mesh mismatch problem based on the particle method to connect preCICE. In addition, we conduct experiments to verify the performance of the PMC method, in which the fluid and the structure is discretized by SPH and the Finite Element Method (FEM), respectively. The results show that the PMC method given in this paper is effective for solving FSI problems. Finally, our source code for the SPH fluid adapter is open-source and available on GitHub for further developing preCICE compatibility with more meshless methods.
△ Less
Submitted 28 April, 2024;
originally announced April 2024.
-
On the local solvability and stability of the partial inverse problems for the non-self-adjoint Sturm-Liouville operators with a discontinuity
Authors:
Xiao-Chuan Xu,
Chuan-Fu Yang,
Natalia Pavlovna Bondarenko
Abstract:
In this work, we study the inverse spectral problems for the Sturm-Liouville operators on [0,1] with complex coefficients and a discontinuity at $x=a\in(0,1)$. Assume that the potential on (a,1) and some parameters in the discontinuity and boundary conditions are given. We recover the potential on (0,a) and the other parameters from the eigenvalues. This is the so-called partial inverse problem. T…
▽ More
In this work, we study the inverse spectral problems for the Sturm-Liouville operators on [0,1] with complex coefficients and a discontinuity at $x=a\in(0,1)$. Assume that the potential on (a,1) and some parameters in the discontinuity and boundary conditions are given. We recover the potential on (0,a) and the other parameters from the eigenvalues. This is the so-called partial inverse problem. The local solvability and stability of the partial inverse problems are obtained for $a\in(0,1)$, in which the error caused by the given partial potential is considered. As a by-product, we also obtain two new uniqueness theorems for the partial inverse problem.
△ Less
Submitted 12 July, 2024; v1 submitted 20 April, 2024;
originally announced April 2024.
-
FCNCP: A Coupled Nonnegative CANDECOMP/PARAFAC Decomposition Based on Federated Learning
Authors:
Yukai Cai,
Hang Liu,
Xiulin Wang,
Hongjin Li,
Ziyi Wang,
Chuanshuai Yang,
Fengyu Cong
Abstract:
In the field of brain science, data sharing across servers is becoming increasingly challenging due to issues such as industry competition, privacy security, and administrative procedure policies and regulations. Therefore, there is an urgent need to develop new methods for data analysis and processing that enable scientific collaboration without data sharing. In view of this, this study proposes…
▽ More
In the field of brain science, data sharing across servers is becoming increasingly challenging due to issues such as industry competition, privacy security, and administrative procedure policies and regulations. Therefore, there is an urgent need to develop new methods for data analysis and processing that enable scientific collaboration without data sharing. In view of this, this study proposes to study and develop a series of efficient non-negative coupled tensor decomposition algorithm frameworks based on federated learning called FCNCP for the EEG data arranged on different servers. It combining the good discriminative performance of tensor decomposition in high-dimensional data representation and decomposition, the advantages of coupled tensor decomposition in cross-sample tensor data analysis, and the features of federated learning for joint modelling in distributed servers. The algorithm utilises federation learning to establish coupling constraints for data distributed across different servers. In the experiments, firstly, simulation experiments are carried out using simulated data, and stable and consistent decomposition results are obtained, which verify the effectiveness of the proposed algorithms in this study. Then the FCNCP algorithm was utilised to decompose the fifth-order event-related potential (ERP) tensor data collected by applying proprioceptive stimuli on the left and right hands. It was found that contralateral stimulation induced more symmetrical components in the activation areas of the left and right hemispheres. The conclusions drawn are consistent with the interpretations of related studies in cognitive neuroscience, demonstrating that the method can efficiently process higher-order EEG data and that some key hidden information can be preserved.
△ Less
Submitted 18 April, 2024;
originally announced April 2024.
-
Undecidability of tiling the plane with a fixed number of Wang bars
Authors:
Chao Yang,
Zhujun Zhang
Abstract:
To study the fixed parameter undecidability of tiling problem for a set of Wang tiles, Jeandel and Rolin show that the tiling problem for a set of 44 Wang bars is undecidable. In this paper, we improve their result by proving that whether a set of 29 Wang bars can tile the plane is undecidable. As a consequence, the tiling problem for a set of Wang tiles with color deficiency of 25 is also undecid…
▽ More
To study the fixed parameter undecidability of tiling problem for a set of Wang tiles, Jeandel and Rolin show that the tiling problem for a set of 44 Wang bars is undecidable. In this paper, we improve their result by proving that whether a set of 29 Wang bars can tile the plane is undecidable. As a consequence, the tiling problem for a set of Wang tiles with color deficiency of 25 is also undecidable.
△ Less
Submitted 6 April, 2024;
originally announced April 2024.
-
Hook-Lengths, Symplectic/Orthogonal Contents and Amdeberhan's Conjectures
Authors:
Chenglang Yang
Abstract:
The symplectic/orthogonal contents of partitions are related to the dimensions of irreducible representations of symplectic/orthogonal groups. In 2012, motivated by Nekrasov-Okounkov's hook-length formula and Stanley's hook-content formula, Amdeberhan proposed several conjectures about infinite product formulas for certain generating functions of hook-lengths and symplectic/orthogonal contents. So…
▽ More
The symplectic/orthogonal contents of partitions are related to the dimensions of irreducible representations of symplectic/orthogonal groups. In 2012, motivated by Nekrasov-Okounkov's hook-length formula and Stanley's hook-content formula, Amdeberhan proposed several conjectures about infinite product formulas for certain generating functions of hook-lengths and symplectic/orthogonal contents. Some special cases of his conjectures were recently proved by Amdeberhan, Andrews and Ballantine. In this paper, we prove the general cases of Amdeberhan's conjectures.
△ Less
Submitted 24 April, 2024; v1 submitted 2 April, 2024;
originally announced April 2024.
-
On the Korteweg-de Vries limit for the Boussinesq equation
Authors:
Younghun Hong,
Changhun Yang
Abstract:
The Korteweg-de Vries (KdV) equation is known as a universal equation describing various long waves in dispersive systems. In this article, we prove that in a certain scaling regime, a large class of rough solutions to the Boussinesq equation are approximated by the sums of two counter-propagating waves solving the KdV equations. It extends the earlier result by \cite{Schneider1998} to slightly mo…
▽ More
The Korteweg-de Vries (KdV) equation is known as a universal equation describing various long waves in dispersive systems. In this article, we prove that in a certain scaling regime, a large class of rough solutions to the Boussinesq equation are approximated by the sums of two counter-propagating waves solving the KdV equations. It extends the earlier result by \cite{Schneider1998} to slightly more regular than $L^2$-solutions. Our proof is based on robust Fourier analysis methods developed for the low regularity theory of nonlinear dispersive equations.
△ Less
Submitted 11 April, 2024; v1 submitted 25 March, 2024;
originally announced March 2024.
-
Michell Truss and From 1-beam to k-beam
Authors:
Chengcheng Yang
Abstract:
This paper generalizes the Michell Truss problem and Gangbo's paper from 1-dimension to higher dimensions using geometric measure theory.
Given an elastic surface $S$ made of $(k-1)$-beams under an equilibriated system $F$ of external forces, then we ask the following two questions:
1. What are the necessary and sufficient conditions for the existence of an elastic body made of $k$-beams whose…
▽ More
This paper generalizes the Michell Truss problem and Gangbo's paper from 1-dimension to higher dimensions using geometric measure theory.
Given an elastic surface $S$ made of $(k-1)$-beams under an equilibriated system $F$ of external forces, then we ask the following two questions:
1. What are the necessary and sufficient conditions for the existence of an elastic body made of $k$-beams whose forces on the surface balance $F$ and whose surfaces consist of $S$.
2. What is an optimal design so that the total cost is a minimum?
We've solved the existence question completely; and research is still in progress for the minimal question. In particular when $k=1$, it involves a system of beams joining a given finite collection of pointed forces. It was first introduced by A. Michell in 1904, then used in mechanical engineering, and recently popularized in many pure mathematics works by W. Gangbo, Prager, and others. Here we are going to generalize them to higher dimensional cases. We have already found the minimal solutions in terms of the flat chain complex and vector-valued currents. Right now we are studying the Calibration theory for future directions. I appreciate the discussion with Prof. Robert Hardt!
△ Less
Submitted 23 March, 2024;
originally announced March 2024.
-
A proof of Ollinger's conjecture: undecidability of tiling the plane with a set of $8$ polyominoes
Authors:
Chao Yang,
Zhujun Zhang
Abstract:
We give a proof of Ollinger's conjecture that the problem of tiling the plane with translated copies of a set of $8$ polyominoes is undecidable. The techniques employed in our proof include a different orientation for simulating the Wang tiles in polyomino and a new method for encoding the colors of Wang tiles.
We give a proof of Ollinger's conjecture that the problem of tiling the plane with translated copies of a set of $8$ polyominoes is undecidable. The techniques employed in our proof include a different orientation for simulating the Wang tiles in polyomino and a new method for encoding the colors of Wang tiles.
△ Less
Submitted 20 March, 2024;
originally announced March 2024.
-
λ-Biharmonic hypersurfaces in the product space L^{m}\times \mathbb{R}
Authors:
Chao Yang,
Zhen Zhao
Abstract:
In this paper, we study λ-biharmonic hypersurfaces in the product space L^{m}\times\mathbb{R}, where L^{m} is an Einstein space and \mathbb{R} is a real line. We prove that λ-biharmonic hypersurfaces with constant mean curvature in L^{m}\times\mathbb{R} are either minimal or vertical cylinders, and obtain some classification results for λ$-biharmonic hypersurfaces under various constraints. Furthe…
▽ More
In this paper, we study λ-biharmonic hypersurfaces in the product space L^{m}\times\mathbb{R}, where L^{m} is an Einstein space and \mathbb{R} is a real line. We prove that λ-biharmonic hypersurfaces with constant mean curvature in L^{m}\times\mathbb{R} are either minimal or vertical cylinders, and obtain some classification results for λ$-biharmonic hypersurfaces under various constraints. Furthermore, we investigate λ-biharmonic hypersurfaces in the product space L^{m}(c)\times\mathbb{R}, where L^{m}(c) is a space form with constant sectional curvature c, and categorize hypersurfaces that are either totally umbilical or semi-parallel.
△ Less
Submitted 16 March, 2024;
originally announced March 2024.
-
PMCV hypersurfaces in non-flat pseudo-Riemannian space forms
Authors:
Chao Yang,
Jiancheng Liu,
Li Du
Abstract:
In this paper, we prove that PMCV (i.e. Δ\vec{H} is proportional to \vec{H}) hypersurface M^n_r of a non-flat pseudo-Riemannian space form N^{n+1}_s(c) with at most two distinct principal curvatures is minimal or locally isoparametric, and compute the mean curvature for the isoparametric ones. As an application, we give full classification results of such non-minimal Lorentzian hypersurfaces of no…
▽ More
In this paper, we prove that PMCV (i.e. Δ\vec{H} is proportional to \vec{H}) hypersurface M^n_r of a non-flat pseudo-Riemannian space form N^{n+1}_s(c) with at most two distinct principal curvatures is minimal or locally isoparametric, and compute the mean curvature for the isoparametric ones. As an application, we give full classification results of such non-minimal Lorentzian hypersurfaces of non-flat Lorentz space forms.
△ Less
Submitted 12 March, 2024;
originally announced March 2024.
-
HOSCF: Efficient decoupling algorithms for finding the best rank-one approximation of higher-order tensors
Authors:
Chuanfu Xiao,
Zeyu Li,
Chao Yang
Abstract:
Best rank-one approximation is one of the most fundamental tasks in tensor computation. In order to fully exploit modern multi-core parallel computers, it is necessary to develop decoupling algorithms for computing the best rank-one approximation of higher-order tensors at large scales. In this paper, we first build a bridge between the rank-one approximation of tensors and the eigenvector-depende…
▽ More
Best rank-one approximation is one of the most fundamental tasks in tensor computation. In order to fully exploit modern multi-core parallel computers, it is necessary to develop decoupling algorithms for computing the best rank-one approximation of higher-order tensors at large scales. In this paper, we first build a bridge between the rank-one approximation of tensors and the eigenvector-dependent nonlinear eigenvalue problem (NEPv), and then develop an efficient decoupling algorithm, namely the higher-order self-consistent field (HOSCF) algorithm, inspired by the famous self-consistent field (SCF) iteration frequently used in computational chemistry. The convergence theory of the HOSCF algorithm and an estimation of the convergence speed are further presented. In addition, we propose an improved HOSCF (iHOSCF) algorithm that incorporates the Rayleigh quotient iteration, which can significantly accelerate the convergence of HOSCF. Numerical experiments show that the proposed algorithms can efficiently converge to the best rank-one approximation of both synthetic and real-world tensors and can scale with high parallel scalability on a modern parallel computer.
△ Less
Submitted 4 March, 2024;
originally announced March 2024.
-
Atropos-k is PSPACE-complete
Authors:
Chao Yang,
Zhujun Zhang
Abstract:
Burke and Teng introduced a two-player combinatorial game Atropos based on Sperner's lemma, and showed that deciding whether one has a winning strategy for Atropos is PSPACE-complete. In the original Atropos game, the players must color a node adjacent to the last colored node. Burke and Teng also mentioned a variant Atropos-k in which each move is at most of distance k of the previous move, and a…
▽ More
Burke and Teng introduced a two-player combinatorial game Atropos based on Sperner's lemma, and showed that deciding whether one has a winning strategy for Atropos is PSPACE-complete. In the original Atropos game, the players must color a node adjacent to the last colored node. Burke and Teng also mentioned a variant Atropos-k in which each move is at most of distance k of the previous move, and asked a question on determining the computational complexity of this variant. In this paper, we answer this question by showing that for any fixed integer k (k>=2), Atropos-k is PSPACE-complete by reduction from True Quantified Boolean Formula (TQBF).
△ Less
Submitted 3 March, 2024;
originally announced March 2024.
-
A class of higher order inverse spectral problems
Authors:
Ai-Wei Guan,
Chuan-Fu Yang,
Natalia P. Bondarenko
Abstract:
In this paper, we consider the recovery of third-order differential operators from two spectra, as well as fourth-order or fifth-order differential operators from three spectra, where these differential operators are endowed with complex-valued distributional coefficients. For the case of multiple spectra, we first establish the relationship between spectra and the Weyl-Yurko matrix. Secondly, we…
▽ More
In this paper, we consider the recovery of third-order differential operators from two spectra, as well as fourth-order or fifth-order differential operators from three spectra, where these differential operators are endowed with complex-valued distributional coefficients. For the case of multiple spectra, we first establish the relationship between spectra and the Weyl-Yurko matrix. Secondly, we prove the uniqueness theorem for the solution of the inverse problems. Our approach allows us to obtain results for the general case of complex-valued distributional coefficients.
△ Less
Submitted 28 February, 2024;
originally announced February 2024.
-
The sign of linear periods
Authors:
U. K. Anandavardhanan,
Hengfei Lu,
Nadir Matringe,
Vincent Sécherre,
Chang Yang
Abstract:
Let $G$ be a group with subgroup $H$, and let $(π,V)$ be a complex representation of $G$. The natural action of the normalizer $N$ of $H$ in $G$ on the space $\mathrm{Hom}_H(π,\mathbb{C})$ of $H$-invariant linear forms on $V$, provides a representation $χ_π$ of $N$ trivial on $H$, which is a character when $\mathrm{Hom}_H(π,\mathbb{C})$ is one dimensional. If moreover $G$ is a reductive group over…
▽ More
Let $G$ be a group with subgroup $H$, and let $(π,V)$ be a complex representation of $G$. The natural action of the normalizer $N$ of $H$ in $G$ on the space $\mathrm{Hom}_H(π,\mathbb{C})$ of $H$-invariant linear forms on $V$, provides a representation $χ_π$ of $N$ trivial on $H$, which is a character when $\mathrm{Hom}_H(π,\mathbb{C})$ is one dimensional. If moreover $G$ is a reductive group over a local field, and $π$ is smooth irreducible, it is an interesting problem to express $χ_π$ in terms of the possibly conjectural Langlands parameter $φ_π$ of $π$. In this paper we consider the following situation: $G=\mathrm{GL}_m(D)$ for $D$ a central division algebra of dimension $d^2$ over a local field $F$ of characteristic zero, $H$ is the centralizer of a non central element $δ\in G$ such that $δ^2$ is in the center of $G$, and $π$ has generic Jacquet-Langlands transfer to $\mathrm{GL}_{md}(F)$. In this setting the space $\mathrm{Hom}_H(π,\mathbb{C})$ is at most one dimensional. When $\mathrm{Hom}_H(π,\mathbb{C})\simeq \mathbb{C}$ and $H\neq N$, we prove that the value of the $χ_π$ on the non trivial class of $\frac{N}{H}$ is $(-1)^mε(φ_π)$ where $ε(φ_π)$ is the root number of $φ_π$. Along the way we extend many useful multiplicity one results for linear and Shalika models to the case of non split $G$. When $F$ is $p$-adic we also classify standard modules with linear periods and Shalika models, which are new results even when $D=F$.
△ Less
Submitted 16 July, 2024; v1 submitted 19 February, 2024;
originally announced February 2024.
-
Deep adaptive sampling for surrogate modeling without labeled data
Authors:
Xili Wang,
Kejun Tang,
Jiayu Zhai,
Xiaoliang Wan,
Chao Yang
Abstract:
Surrogate modeling is of great practical significance for parametric differential equation systems. In contrast to classical numerical methods, using physics-informed deep learning methods to construct simulators for such systems is a promising direction due to its potential to handle high dimensionality, which requires minimizing a loss over a training set of random samples. However, the random s…
▽ More
Surrogate modeling is of great practical significance for parametric differential equation systems. In contrast to classical numerical methods, using physics-informed deep learning methods to construct simulators for such systems is a promising direction due to its potential to handle high dimensionality, which requires minimizing a loss over a training set of random samples. However, the random samples introduce statistical errors, which may become the dominant errors for the approximation of low-regularity and high-dimensional problems. In this work, we present a deep adaptive sampling method for surrogate modeling ($\text{DAS}^2$), where we generalize the deep adaptive sampling (DAS) method [62] [Tang, Wan and Yang, 2023] to build surrogate models for low-regularity parametric differential equations. In the parametric setting, the residual loss function can be regarded as an unnormalized probability density function (PDF) of the spatial and parametric variables. This PDF is approximated by a deep generative model, from which new samples are generated and added to the training set. Since the new samples match the residual-induced distribution, the refined training set can further reduce the statistical error in the current approximate solution. We demonstrate the effectiveness of $\text{DAS}^2$ with a series of numerical experiments, including the parametric lid-driven 2D cavity flow problem with a continuous range of Reynolds numbers from 100 to 1000.
△ Less
Submitted 17 February, 2024;
originally announced February 2024.
-
An Efficient Quantum Circuit for Block Encoding a Pairing Hamiltonian
Authors:
Diyi Liu,
Weijie Du,
Lin Lin,
James P. Vary,
Chao Yang
Abstract:
We present an efficient quantum circuit for block encoding pairing Hamiltonian often studied in nuclear physics. Our block encoding scheme does not require mapping the creation and annihilation operators to the Pauli operators and representing the Hamiltonian as a linear combination of unitaries. Instead, we show how to encode the Hamiltonian directly using controlled swap operations. We analyze t…
▽ More
We present an efficient quantum circuit for block encoding pairing Hamiltonian often studied in nuclear physics. Our block encoding scheme does not require mapping the creation and annihilation operators to the Pauli operators and representing the Hamiltonian as a linear combination of unitaries. Instead, we show how to encode the Hamiltonian directly using controlled swap operations. We analyze the gate complexity of the block encoding circuit and show that it scales polynomially with respect to the number of qubits required to represent a quantum state associated with the pairing Hamiltonian. We also show how the block encoding circuit can be combined with the quantum singular value transformation to construct an efficient quantum circuit for approximating the density of states of a pairing Hamiltonian. The techniques presented can be extended to encode more general second-quantized Hamiltonians.
△ Less
Submitted 21 February, 2024; v1 submitted 17 February, 2024;
originally announced February 2024.
-
Friends-and-strangers is PSPACE-complete
Authors:
Chao Yang,
Zhujun Zhang
Abstract:
In this paper, we show that the friends-and-strangers problem is PSPACE-complete by reduction from the Ncl (non-deterministic constraint logic) problem.
In this paper, we show that the friends-and-strangers problem is PSPACE-complete by reduction from the Ncl (non-deterministic constraint logic) problem.
△ Less
Submitted 5 February, 2024;
originally announced February 2024.
-
Enhancing Stochastic Gradient Descent: A Unified Framework and Novel Acceleration Methods for Faster Convergence
Authors:
Yichuan Deng,
Zhao Song,
Chiwun Yang
Abstract:
Based on SGD, previous works have proposed many algorithms that have improved convergence speed and generalization in stochastic optimization, such as SGDm, AdaGrad, Adam, etc. However, their convergence analysis under non-convex conditions is challenging. In this work, we propose a unified framework to address this issue. For any first-order methods, we interpret the updated direction $g_t$ as th…
▽ More
Based on SGD, previous works have proposed many algorithms that have improved convergence speed and generalization in stochastic optimization, such as SGDm, AdaGrad, Adam, etc. However, their convergence analysis under non-convex conditions is challenging. In this work, we propose a unified framework to address this issue. For any first-order methods, we interpret the updated direction $g_t$ as the sum of the stochastic subgradient $\nabla f_t(x_t)$ and an additional acceleration term $\frac{2|\langle v_t, \nabla f_t(x_t) \rangle|}{\|v_t\|_2^2} v_t$, thus we can discuss the convergence by analyzing $\langle v_t, \nabla f_t(x_t) \rangle$. Through our framework, we have discovered two plug-and-play acceleration methods: \textbf{Reject Accelerating} and \textbf{Random Vector Accelerating}, we theoretically demonstrate that these two methods can directly lead to an improvement in convergence rate.
△ Less
Submitted 2 February, 2024;
originally announced February 2024.
-
A parameter-free approach for solving SOS-convex semi-algebraic fractional programs
Authors:
Chengmiao Yang,
Liguo Jiao,
Jae Hyoung Lee
Abstract:
In this paper, we study a class of nonsmooth fractional programs {\rm (FP, for short)} with SOS-convex semi-algebraic functions. Under suitable assumptions, we derive a strong duality result between the problem (FP) and its semidefinite programming (SDP) relaxations. Remarkably, we extract an optimal solution of the problem (FP) by solving one and only one associated SDP problem. Numerical example…
▽ More
In this paper, we study a class of nonsmooth fractional programs {\rm (FP, for short)} with SOS-convex semi-algebraic functions. Under suitable assumptions, we derive a strong duality result between the problem (FP) and its semidefinite programming (SDP) relaxations. Remarkably, we extract an optimal solution of the problem (FP) by solving one and only one associated SDP problem. Numerical examples are also given.
△ Less
Submitted 29 January, 2024;
originally announced January 2024.
-
On A Proof of the ADKMV Conjecture
Authors:
Zhiyuan Wang,
Chenglang Yang,
Jian Zhou
Abstract:
We present a mathematical proof of a conjectural formula due to Aganagic, Dijkgraaf, Klemm, Mariño and Vafa, expressing the topological vertex as a Bogoliubov transform of the fermionic vacuum. In our proof we introduce a boson-fermionic field assignment which generalizes the well-known boson-fermion correspondence. The proof also works for the generalization to the framed topological vertex made…
▽ More
We present a mathematical proof of a conjectural formula due to Aganagic, Dijkgraaf, Klemm, Mariño and Vafa, expressing the topological vertex as a Bogoliubov transform of the fermionic vacuum. In our proof we introduce a boson-fermionic field assignment which generalizes the well-known boson-fermion correspondence. The proof also works for the generalization to the framed topological vertex made by Deng and Zhou. As a consequence, partition functions of toric Calabi-Yau threefolds are related to tau-functions of multi-component KP hierarchy.
△ Less
Submitted 23 January, 2024;
originally announced January 2024.
-
Scattering for 2d semi-relativistic Hartree equations with short range potential
Authors:
Changhun Yang
Abstract:
We study the long time behavior of small solutions to the semi-relativistic Hartree equations in two dimension. The nonlinear term is convolved with the singular potential $|x|^{-γ}$ for $1<γ<2$, which is referred to as short-range interaction potential in the sense of scattering phenomenon. We establish the scattering results for small solutions in a weighted space, in other words, we prove that…
▽ More
We study the long time behavior of small solutions to the semi-relativistic Hartree equations in two dimension. The nonlinear term is convolved with the singular potential $|x|^{-γ}$ for $1<γ<2$, which is referred to as short-range interaction potential in the sense of scattering phenomenon. We establish the scattering results for small solutions in a weighted space, in other words, we prove that the nonlinear solutions exist globally and behave asymptotically like a linear solution whenever the initial data is sufficiently small. To achieve this, we should obtain time decay estimates for the nonlinear term which is integrable. A key observation is that the loss in time in the course of weighted energy estimates can be recovered by the space resonance method with the help of null structure.
△ Less
Submitted 21 December, 2023;
originally announced December 2023.
-
Ext-distinction for $p$-adic symmetric spaces
Authors:
Chang Yang
Abstract:
Let $G/H$ be a $p$-adic symmetric space. We compute explicitly the higher relative extension groups for all discrete series representations of $G$ in two examples: the symplectic case and the linear case. The results have immediate applications to the computation of the Euler-Poincaré pairing, the alternating sum of the dimensions of the Ext-groups. In the linear case we confirm a conjecture of Wa…
▽ More
Let $G/H$ be a $p$-adic symmetric space. We compute explicitly the higher relative extension groups for all discrete series representations of $G$ in two examples: the symplectic case and the linear case. The results have immediate applications to the computation of the Euler-Poincaré pairing, the alternating sum of the dimensions of the Ext-groups. In the linear case we confirm a conjecture of Wan which asserts the equality of the Euler-Poincaré pairing with the geometric multiplicity for any irreducible representation of $G$. In the symplectic case we reduce the verification of Wan's conjecture to the case of discrete series representations. We also determine all relatively supercuspidal representations in both cases.
△ Less
Submitted 18 December, 2023;
originally announced December 2023.
-
On two families of Nekrasov-Okounkov type formulas
Authors:
Chenglang Yang
Abstract:
In this paper, we use the vacuum expectation value formula of the topological vertex and its rotation symmetry to derive two families of Nekrasov-Okounkov type formulas. Each family of formulas depends on $2N+1$ parameters for a positive integer $N$.
In this paper, we use the vacuum expectation value formula of the topological vertex and its rotation symmetry to derive two families of Nekrasov-Okounkov type formulas. Each family of formulas depends on $2N+1$ parameters for a positive integer $N$.
△ Less
Submitted 13 December, 2023;
originally announced December 2023.
-
A remark on certain restricted plane partitions and crystal melting model
Authors:
Chenglang Yang
Abstract:
In this paper, we provide formulas calculating the partition functions of two types of plane partitions using the crystal melting model method introduced by Okounkov, Reshetikhin and Vafa. As applications, we obtain a product formula for the partition function of the plane partitions with a limit shape boundary. A corollary of this formula is the demonstration of the equivalence between this parti…
▽ More
In this paper, we provide formulas calculating the partition functions of two types of plane partitions using the crystal melting model method introduced by Okounkov, Reshetikhin and Vafa. As applications, we obtain a product formula for the partition function of the plane partitions with a limit shape boundary. A corollary of this formula is the demonstration of the equivalence between this partition function and the open-closed string amplitude of the double$-\mathbb{P}^1$ model. We also derive a product formula for the partition function of symmetric plane partitions with a limit shape boundary.
△ Less
Submitted 5 December, 2023;
originally announced December 2023.
-
Many vertex-disjoint even cycles of fixed length in a graph
Authors:
Jianfeng Hou,
Caiyun Hu,
Heng Li,
Xizhi Liu,
Caihong Yang,
Yixiao Zhang
Abstract:
For every integer $k \ge 3$, we determine the extremal structure of an $n$-vertex graph with at most $t$ vertex-disjoint copies of $C_{2k}$ when $n$ is sufficiently large and $t$ lies in the interval $\left[\frac{\mathrm{ex}(n,C_{2k})}{\varepsilon n}, \varepsilon n\right]$, where $\varepsilon>0$ is a constant depending only on $k$. The question for $k = 2$ and…
▽ More
For every integer $k \ge 3$, we determine the extremal structure of an $n$-vertex graph with at most $t$ vertex-disjoint copies of $C_{2k}$ when $n$ is sufficiently large and $t$ lies in the interval $\left[\frac{\mathrm{ex}(n,C_{2k})}{\varepsilon n}, \varepsilon n\right]$, where $\varepsilon>0$ is a constant depending only on $k$. The question for $k = 2$ and $t = o\left(\frac{\mathrm{ex}(n,C_{2k})}{n}\right)$ was explored in prior work~\cite{HHLLYZ23a}, revealing different extremal structures in these cases. Our result can be viewed as an extension of the theorems by Egawa~\cite{Ega96} and Verstraëte~\cite{Ver03}, where the focus was on the existence of many vertex-disjoint cycles of the same length without any length constraints.
△ Less
Submitted 25 November, 2023;
originally announced November 2023.
-
Toward a density Corrádi--Hajnal theorem for degenerate hypergraphs
Authors:
Jianfeng Hou,
Caiyun Hu,
Heng Li,
Xizhi Liu,
Caihong Yang,
Yixiao Zhang
Abstract:
Given an $r$-graph $F$ with $r \ge 2$, let $\mathrm{ex}(n, (t+1) F)$ denote the maximum number of edges in an $n$-vertex $r$-graph with at most $t$ pairwise vertex-disjoint copies of $F$. Extending several old results and complementing prior work [J. Hou, H. Li, X. Liu, L.-T. Yuan, and Y. Zhang. A step towards a general density Corrádi--Hajnal theorem. arXiv:2302.09849, 2023.] on nondegenerate hyp…
▽ More
Given an $r$-graph $F$ with $r \ge 2$, let $\mathrm{ex}(n, (t+1) F)$ denote the maximum number of edges in an $n$-vertex $r$-graph with at most $t$ pairwise vertex-disjoint copies of $F$. Extending several old results and complementing prior work [J. Hou, H. Li, X. Liu, L.-T. Yuan, and Y. Zhang. A step towards a general density Corrádi--Hajnal theorem. arXiv:2302.09849, 2023.] on nondegenerate hypergraphs, we initiate a systematic study on $\mathrm{ex}(n, (t+1) F)$ for degenerate hypergraphs $F$. For a broad class of degenerate hypergraphs $F$, we present near-optimal upper bounds for $\mathrm{ex}(n, (t+1) F)$ when $n$ is sufficiently large and $t$ lies in intervals $\left[0, \frac{\varepsilon \cdot \mathrm{ex}(n,F)}{n^{r-1}}\right]$, $\left[\frac{\mathrm{ex}(n,F)}{\varepsilon n^{r-1}}, \varepsilon n \right]$, and $\left[ (1-\varepsilon)\frac{n}{v(F)}, \frac{n}{v(F)} \right]$, where $\varepsilon > 0$ is a constant depending only on $F$. Our results reveal very different structures for extremal constructions across the three intervals, and we provide characterizations of extremal constructions within the first interval. Additionally, for graphs, we offer a characterization of extremal constructions within the second interval. Our proof for the first interval also applies to a special class of nondegenerate hypergraphs, including those with undetermined Turán densities, partially improving a result in [J. Hou, H. Li, X. Liu, L.-T. Yuan, and Y. Zhang. A step towards a general density Corrádi--Hajnal theorem. arXiv:2302.09849, 2023.]
△ Less
Submitted 25 November, 2023;
originally announced November 2023.
-
Theory of Infectious Diseases with Testing and Testing-less Covid-19 Endemic
Authors:
Bo Deng,
Chayu Yang
Abstract:
What is the long term dynamics of the Covid-19 pandemic? How will it end? Here we constructed an infectious disease model with testing and analyzed the existence and stability of its endemic states. For a large parameter set, including those relevant to the SARS-CoV-2 virus, we demonstrated the existence of one endemic equilibrium without testing and one endemic equilibrium with testing and proved…
▽ More
What is the long term dynamics of the Covid-19 pandemic? How will it end? Here we constructed an infectious disease model with testing and analyzed the existence and stability of its endemic states. For a large parameter set, including those relevant to the SARS-CoV-2 virus, we demonstrated the existence of one endemic equilibrium without testing and one endemic equilibrium with testing and proved their local and global stabilities for some cases. Our results suggest that the pandemic is to end with a testing-less endemic state through a novel and surprising mechanism called stochastic trapping.
△ Less
Submitted 14 November, 2023;
originally announced November 2023.
-
Nash's Existence Theorem for Non-compact Strategy Sets
Authors:
Xinyu Zhang,
Cheng Chen,
Chunyan Yang
Abstract:
This paper generalizes the Fan-Knaster-Kuratowski-Mazurkiewicz (FKKM) lemma to the case of weak topology, and obtains the Ky Fan minimax inequality defined on non-empty non-compact convex subsets in reflexive Banach spaces, then we apply it to game theory and obtain Nash's existence theorem for non-compact strategy sets, together with John von Neumann's existence theorem in two-player zero-sum gam…
▽ More
This paper generalizes the Fan-Knaster-Kuratowski-Mazurkiewicz (FKKM) lemma to the case of weak topology, and obtains the Ky Fan minimax inequality defined on non-empty non-compact convex subsets in reflexive Banach spaces, then we apply it to game theory and obtain Nash's existence theorem for non-compact strategy sets, together with John von Neumann's existence theorem in two-player zero-sum games.
△ Less
Submitted 12 November, 2023;
originally announced November 2023.
-
Asymptotic analysis for Bloch electrons with Weyl nodes
Authors:
Jianfeng Lu,
Changhe Yang,
Zhennan Zhou
Abstract:
In this paper, we study the semiclassical behavior of Bloch electrons in the presence of Weyl nodes, which are singular points in the band structure of certain materials. We carry out asymptotic analysis and present a rigorous derivation of the semiclassical asymptotic expansion of the current of Bloch electrons with the presence of Weyl nodes. The analysis shows that the current contains two part…
▽ More
In this paper, we study the semiclassical behavior of Bloch electrons in the presence of Weyl nodes, which are singular points in the band structure of certain materials. We carry out asymptotic analysis and present a rigorous derivation of the semiclassical asymptotic expansion of the current of Bloch electrons with the presence of Weyl nodes. The analysis shows that the current contains two parts, one independent of the Weyl nodes and the other a contribution from the singular points. This work provides a theoretical foundation towards a rigorous justification of recent scientific discoveries in Weyl semimetals. The main innovation of this paper is a new strategy to deal with the singular points with quantitative estimates, which may have broader applications in multiscale models with singularities.
△ Less
Submitted 4 November, 2023;
originally announced November 2023.
-
The free boundary for a semilinear non-homogeneous Bernoulli problem
Authors:
Lili Du,
Chunlei Yang
Abstract:
In the classical homogeneous one-phase Bernoulli-type problem, the free boundary consists of a "regular" part and a "singular" part, as Alt and Caffarelli have shown in their pioneer work (J. Reine Angew. Math., 325, 105-144, 1981) that regular points are $C^{1,γ}$ in two-dimensions. Later, Weiss (J. Geom. Anal., 9, 317-326, 1999) first realized that in higher dimensions a critical dimension…
▽ More
In the classical homogeneous one-phase Bernoulli-type problem, the free boundary consists of a "regular" part and a "singular" part, as Alt and Caffarelli have shown in their pioneer work (J. Reine Angew. Math., 325, 105-144, 1981) that regular points are $C^{1,γ}$ in two-dimensions. Later, Weiss (J. Geom. Anal., 9, 317-326, 1999) first realized that in higher dimensions a critical dimension $d^{*}$ exists so that the singularities of the free boundary can only occur when $d\geqslant d^{*}$.
In this paper, we consider a non-homogeneous semilinear one-phase Bernoulli-type problem, and we show that the free boundary is a disjoint union of a regular and a singular set. Moreover, the regular set is locally the graph of a $C^{1,γ}$ function for some $γ\in(0,1)$. In addition, there exists a critical dimension $d^{*}$ so that the singular set is empty if $d<d^{*}$, discrete if $d=d^{*}$ and of locally finite $\mathcal{H}^{d-d^{*}}$ Hausdorff measure if $d>d^{*}$. As a byproduct, we relate the existence of viscosity solutions of a non-homogeneous problem to the Weiss-boundary adjusted energy, which provides an alternative proof to existence of viscosity solutions for non-homogeneous problems.
△ Less
Submitted 8 May, 2024; v1 submitted 31 October, 2023;
originally announced November 2023.
-
Counting triangles in graphs without vertex disjoint odd cycles
Authors:
Jianfeng Hou,
Caihong Yang,
Qinghou Zeng
Abstract:
Given two graphs $H$ and $F$, the maximum possible number of copies of $H$ in an $F$-free graph on $n$ vertices is denoted by $\mathrm{ex}(n, H, F)$. Let $(\ell+1) \cdot F$ denote $\ell+1$ vertex disjoint copies of $F$. In this paper, we determine the exact value of $\mathrm{ex}(n, C_3, (\ell+1)\cdot C_{2k+1})$ and its extremal graph, which generalizes some known results.
Given two graphs $H$ and $F$, the maximum possible number of copies of $H$ in an $F$-free graph on $n$ vertices is denoted by $\mathrm{ex}(n, H, F)$. Let $(\ell+1) \cdot F$ denote $\ell+1$ vertex disjoint copies of $F$. In this paper, we determine the exact value of $\mathrm{ex}(n, C_3, (\ell+1)\cdot C_{2k+1})$ and its extremal graph, which generalizes some known results.
△ Less
Submitted 29 October, 2023;
originally announced October 2023.
-
Constructing disjoint Steiner trees in Sierpiński graphs
Authors:
Chenxu Yang,
Ping Li,
Yaping Mao,
Eddie Cheng,
Ralf Klasing
Abstract:
Let $G$ be a graph and $S\subseteq V(G)$ with $|S|\geq 2$. Then the trees $T_1, T_2, \cdots, T_\ell$ in $G$ are \emph{internally disjoint Steiner trees} connecting $S$ (or $S$-Steiner trees) if $E(T_i) \cap E(T_j )=\emptyset$ and $V(T_i)\cap V(T_j)=S$ for every pair of distinct integers $i,j$, $1 \leq i, j \leq \ell$. Similarly, if we only have the condition $E(T_i) \cap E(T_j )=\emptyset$ but wit…
▽ More
Let $G$ be a graph and $S\subseteq V(G)$ with $|S|\geq 2$. Then the trees $T_1, T_2, \cdots, T_\ell$ in $G$ are \emph{internally disjoint Steiner trees} connecting $S$ (or $S$-Steiner trees) if $E(T_i) \cap E(T_j )=\emptyset$ and $V(T_i)\cap V(T_j)=S$ for every pair of distinct integers $i,j$, $1 \leq i, j \leq \ell$. Similarly, if we only have the condition $E(T_i) \cap E(T_j )=\emptyset$ but without the condition $V(T_i)\cap V(T_j)=S$, then they are \emph{edge-disjoint Steiner trees}. The \emph{generalized $k$-connectivity}, denoted by $κ_k(G)$, of a graph $G$, is defined as $κ_k(G)=\min\{κ_G(S)|S \subseteq V(G) \ \textrm{and} \ |S|=k \}$, where $κ_G(S)$ is the maximum number of internally disjoint $S$-Steiner trees. The \emph{generalized local edge-connectivity} $λ_{G}(S)$ is the maximum number of edge-disjoint Steiner trees connecting $S$ in $G$. The {\it generalized $k$-edge-connectivity} $λ_k(G)$ of $G$ is defined as $λ_k(G)=\min\{λ_{G}(S)\,|\,S\subseteq V(G) \ and \ |S|=k\}$. These measures are generalizations of the concepts of connectivity and edge-connectivity, and they and can be used as measures of vulnerability of networks. It is, in general, difficult to compute these generalized connectivities. However, there are precise results for some special classes of graphs. In this paper, we obtain the exact value of $λ_{k}(S(n,\ell))$ for $3\leq k\leq \ell^n$, and the exact value of $κ_{k}(S(n,\ell))$ for $3\leq k\leq \ell$, where $S(n, \ell)$ is the Sierpiński graphs with order $\ell^n$. As a direct consequence, these graphs provide additional interesting examples when $λ_{k}(S(n,\ell))=κ_{k}(S(n,\ell))$. We also study the some network properties of Sierpiński graphs.
△ Less
Submitted 25 October, 2023;
originally announced October 2023.
-
The free boundary of steady axisymmetric inviscid flow with vorticity II: near the non-degenerate point
Authors:
Lili Du,
Chunlei Yang
Abstract:
This is the sequel of the recent work (Du, Huang, Pu, Commun. Math. Phys, 2023, doi: 10.1007/s00220-023-04651-7) on axially symmetric gravity water waves with general vorticities, which has investigated the singular wave profile of the free boundary near the degenerate points. In this companion paper, we are interested in the regularity of the free surface of the water wave near the non-degenerate…
▽ More
This is the sequel of the recent work (Du, Huang, Pu, Commun. Math. Phys, 2023, doi: 10.1007/s00220-023-04651-7) on axially symmetric gravity water waves with general vorticities, which has investigated the singular wave profile of the free boundary near the degenerate points. In this companion paper, we are interested in the regularity of the free surface of the water wave near the non-degenerate point. Precisely, we showed that the free boundary is $C^{1,γ}$ smooth for some $γ\in(0,1)$ near all non-degenerate points. The problem is intrinsically intertwined with the regularity theory of the semilinear Bernoulli-type free boundary problem. Our approach is closely related to the monotonicity formula developed by Weiss in his celebrated work (Weiss, J. Geom. Anal. 9: 317-326, 1999), and to a partial boundary Harnack inequality for the one-phase free boundary problem, which is dedicated to De Silva (De Silva, Interfaces Free Bound. 13, 223-238, 2011). Mathematically, we associate the existence of viscosity solutions with the Weiss boundary-adjusted energy. Compared to the classical approach of Caffarelli (Caffarelli, Ann. Scuola Norm. Sup. Pisa, 15, 583-602, 1988), we provide an alternative proof of the existence of viscosity solutions for a large class of semilinear free boundary problems.
△ Less
Submitted 13 October, 2023;
originally announced October 2023.
-
Learning nonlinear integral operators via Recurrent Neural Networks and its application in solving Integro-Differential Equations
Authors:
Hardeep Bassi,
Yuanran Zhu,
Senwei Liang,
Jia Yin,
Cian C. Reeves,
Vojtech Vlcek,
Chao Yang
Abstract:
In this paper, we propose using LSTM-RNNs (Long Short-Term Memory-Recurrent Neural Networks) to learn and represent nonlinear integral operators that appear in nonlinear integro-differential equations (IDEs). The LSTM-RNN representation of the nonlinear integral operator allows us to turn a system of nonlinear integro-differential equations into a system of ordinary differential equations for whic…
▽ More
In this paper, we propose using LSTM-RNNs (Long Short-Term Memory-Recurrent Neural Networks) to learn and represent nonlinear integral operators that appear in nonlinear integro-differential equations (IDEs). The LSTM-RNN representation of the nonlinear integral operator allows us to turn a system of nonlinear integro-differential equations into a system of ordinary differential equations for which many efficient solvers are available. Furthermore, because the use of LSTM-RNN representation of the nonlinear integral operator in an IDE eliminates the need to perform a numerical integration in each numerical time evolution step, the overall temporal cost of the LSTM-RNN-based IDE solver can be reduced to $O(n_T)$ from $O(n_T^2)$ if a $n_T$-step trajectory is to be computed. We illustrate the efficiency and robustness of this LSTM-RNN-based numerical IDE solver with a model problem. Additionally, we highlight the generalizability of the learned integral operator by applying it to IDEs driven by different external forces. As a practical application, we show how this methodology can effectively solve the Dyson's equation for quantum many-body systems.
△ Less
Submitted 13 October, 2023;
originally announced October 2023.
-
Hopf actions on vertex algebras
Authors:
Chongying Dong,
Li Ren,
Chao Yang
Abstract:
In this article, we investigate Hopf actions on vertex algebras. Our first main result is that every finite-dimensional Hopf algebra that inner faithfully acts on a given π_2-injective vertex algebra must be a group algebra. Secondly, under suitable assumptions, we establish a Schur-Weyl type duality for semisimple Hopf actions on Hopf modules of vertex algebras.
In this article, we investigate Hopf actions on vertex algebras. Our first main result is that every finite-dimensional Hopf algebra that inner faithfully acts on a given π_2-injective vertex algebra must be a group algebra. Secondly, under suitable assumptions, we establish a Schur-Weyl type duality for semisimple Hopf actions on Hopf modules of vertex algebras.
△ Less
Submitted 12 October, 2023;
originally announced October 2023.
-
Rationality of Néron-Tate height over function fields
Authors:
Chengyuan Yang
Abstract:
We prove that the Néron-Tate height of subvarieties are always rational numbers. We use the induction formula, and characterize the canonical metric by theta functions.
We prove that the Néron-Tate height of subvarieties are always rational numbers. We use the induction formula, and characterize the canonical metric by theta functions.
△ Less
Submitted 22 September, 2023;
originally announced September 2023.
-
The non-existence of horizontally flat singularity for steady axisymmetric free surface flows near stagnation points
Authors:
Lili Du,
Chunlei Yang
Abstract:
In a recent research on degenerate points of steady axisymmetric gravity flows with general vorticity, it has been shown that the possible asymptotics near any stagnation point must be the "Stokes corner", the "horizontal cusp", or the "horizontal flatness" (Theorem 1.1, Du, Huang, Pu, Commun. Math. Phys., 400, 2137-2179, 2023). In this paper, we focus on the horizontally flat singularity and show…
▽ More
In a recent research on degenerate points of steady axisymmetric gravity flows with general vorticity, it has been shown that the possible asymptotics near any stagnation point must be the "Stokes corner", the "horizontal cusp", or the "horizontal flatness" (Theorem 1.1, Du, Huang, Pu, Commun. Math. Phys., 400, 2137-2179, 2023). In this paper, we focus on the horizontally flat singularity and show that it is not possible, and therefore the "Stokes corner" and the "cusp" are the only possible asymptotics at the stagnation points. The basic idea of our proof relies on a perturbation of the frequency formula for the two-dimensional problem (Varvaruca, Weiss, Acta Math., 206, 363-403, 2011). Our analysis also suggests that, for steady axisymmetric rotational gravity flows, the singular asymptotic profiles at stagnation points are similar to the scenario observed in two-dimensional waves with vorticity (Varvaruca, Weiss, Ann. I. H. Poincare-AN, 29, 861-885, 2012).
△ Less
Submitted 6 December, 2023; v1 submitted 17 September, 2023;
originally announced September 2023.
-
AONN-2: An adjoint-oriented neural network method for PDE-constrained shape optimization
Authors:
Xili Wang,
Pengfei Yin,
Bo Zhang,
Chao Yang
Abstract:
Shape optimization has been playing an important role in a large variety of engineering applications. Existing shape optimization methods are generally mesh-dependent and therefore encounter challenges due to mesh deformation. To overcome this limitation, we present a new adjoint-oriented neural network method, AONN-2, for PDE-constrained shape optimization problems. This method extends the capabi…
▽ More
Shape optimization has been playing an important role in a large variety of engineering applications. Existing shape optimization methods are generally mesh-dependent and therefore encounter challenges due to mesh deformation. To overcome this limitation, we present a new adjoint-oriented neural network method, AONN-2, for PDE-constrained shape optimization problems. This method extends the capabilities of the original AONN method [1], which is developed for efficiently solving parametric optimal control problems. AONN-2 inherits the direct-adjoint looping (DAL) framework for computing the extremum of an objective functional and the neural network methods for solving complicated PDEs from AONN. Furthermore, AONN-2 expands the application scope to shape optimization by taking advantage of the shape derivatives to optimize the shape represented by discrete boundary points. AONN-2 is a fully mesh-free shape optimization approach, naturally sidestepping issues related to mesh deformation, with no need for maintaining mesh quality and additional mesh corrections. A series of experimental results are presented, highlighting the flexibility, robustness, and accuracy of AONN-2.
△ Less
Submitted 15 September, 2023;
originally announced September 2023.
-
The modified scattering of 2 dimensional semi-relativistic Hartree equations
Authors:
Soonsik Kwon,
Kiyeon Lee,
Changhun Yang
Abstract:
In this paper, we consider the asymptotic behaviors of small solutions to the semi-relativistic Hartree equations in two dimension. The nonlinear term is convolved with the Coulomb potential 1/|x|, and it produces the long-range interaction in the sense of scattering phenomenon. From this observation, one anticipates that small solutions converge to a modified scattering states, although they deca…
▽ More
In this paper, we consider the asymptotic behaviors of small solutions to the semi-relativistic Hartree equations in two dimension. The nonlinear term is convolved with the Coulomb potential 1/|x|, and it produces the long-range interaction in the sense of scattering phenomenon. From this observation, one anticipates that small solutions converge to a modified scattering states, although they decay as linear solutions. We show the global well-posedness and the modified scattering for small solutions in weighted Sobolev spaces. Our proof follows a road map of exploiting the space-time resonance developed by Germain, Masmoudi, and Shatah. Compared to the result in three dimensional case by Pusateri, weaker time decay in two dimension is one of the main obstacles.
△ Less
Submitted 29 July, 2023;
originally announced July 2023.
-
Optimal Surrogate Boundary Selection and Scalability Studies for the Shifted Boundary Method on Octree Meshes
Authors:
Cheng-Hau Yang,
Kumar Saurabh,
Guglielmo Scovazzi,
Claudio Canuto,
Adarsh Krishnamurthy,
Baskar Ganapathysubramanian
Abstract:
The accurate and efficient simulation of Partial Differential Equations (PDEs) in and around arbitrarily defined geometries is critical for many application domains. Immersed boundary methods (IBMs) alleviate the usually laborious and time-consuming process of creating body-fitted meshes around complex geometry models (described by CAD or other representations, e.g., STL, point clouds), especially…
▽ More
The accurate and efficient simulation of Partial Differential Equations (PDEs) in and around arbitrarily defined geometries is critical for many application domains. Immersed boundary methods (IBMs) alleviate the usually laborious and time-consuming process of creating body-fitted meshes around complex geometry models (described by CAD or other representations, e.g., STL, point clouds), especially when high levels of mesh adaptivity are required. In this work, we advance the field of IBM in the context of the recently developed Shifted Boundary Method (SBM). In the SBM, the location where boundary conditions are enforced is shifted from the actual boundary of the immersed object to a nearby surrogate boundary, and boundary conditions are corrected utilizing Taylor expansions. This approach allows choosing surrogate boundaries that conform to a Cartesian mesh without losing accuracy or stability. Our contributions in this work are as follows: (a) we show that the SBM numerical error can be greatly reduced by an optimal choice of the surrogate boundary, (b) we mathematically prove the optimal convergence of the SBM for this optimal choice of the surrogate boundary, (c) we deploy the SBM on massively parallel octree meshes, including algorithmic advances to handle incomplete octrees, and (d) we showcase the applicability of these approaches with a wide variety of simulations involving complex shapes, sharp corners, and different topologies. Specific emphasis is given to Poisson's equation and the linear elasticity equations.
△ Less
Submitted 4 July, 2023;
originally announced July 2023.