Showing 1–50 of 532 results for author: Zhao, Y

Global well-posedness and large time behavior of solutions to the compressible Oldroyd-B model without stress diffusion

Authors: Yajuan Zhao, Yongsheng Li, Tao Liang, Xiaoping Zhai

Abstract: We consider the Cauchy problem ($\mathbb{R}^d, d=2,3$) and the initial boundary values problem ($\mathbb{T}^d, d=2,3$)associated to the compressible Oldroyd-B model which is first derived by Barrett, Lu and Süli [Existence of large-data finite-energy global weak solutions to a compressible Oldroyd-B model, Commun. Math. Sci., 15 (2017), 1265--1323] through micro-macro-analysis of the compressible… ▽ More We consider the Cauchy problem ($\mathbb{R}^d, d=2,3$) and the initial boundary values problem ($\mathbb{T}^d, d=2,3$)associated to the compressible Oldroyd-B model which is first derived by Barrett, Lu and Süli [Existence of large-data finite-energy global weak solutions to a compressible Oldroyd-B model, Commun. Math. Sci., 15 (2017), 1265--1323] through micro-macro-analysis of the compressible Navier-Stokes-Fokker-Planck system.Due to lack of stress diffusion, the problems considered here are very difficult. Exploiting tools from harmonic analysis,notably the Littlewood Paley theory,we first establish the global well-posedness and time-decay rates for solutions of the model with small initial data in Besov spaces with critical regularity.Then, through deeply exploring and fully utilizing the structure of the perturbation system,we obtain the global well-posedness and exponential decay rates for solutions of the model with small initial data in the Soboles spaces $H^3(\mathbb{T}^d)$.Our obtained results improve considerably the recent results by Lu, Pokorný [Anal. Theory Appl., 36 (2020), 348--372],Wang, Wen [Math. Models Methods Appl. Sci., 30 (2020), 139--179],and Liu, Lu, Wen [SIAM J. Math. Anal., 53 (2021), 6216--6242]. △ Less

Submitted 28 August, 2024; originally announced August 2024.

Comments: 51pages

Analytic properties of Speyer's $g$-polynomial of uniform matroids

Authors: Rong Zhang, James Jing Yu Zhao

Abstract: Let $U_{n,d}$ denote the uniform matroid of rank $d$ on $n$ elements. We obtain some recurrence relations satisfied by Speyer's $g$-polynomials $g_{U_{n,d}}(t)$ of $U_{n,d}$. Based on these recurrence relations, we prove that the polynomial $g_{U_{n,d}}(t)$ has only real zeros for any $n-1\geq d\geq 1$. Furthermore, we show that the coefficient of $g_{U_{n,[n/2]}}(t)$ is asymptotically normal by l… ▽ More Let $U_{n,d}$ denote the uniform matroid of rank $d$ on $n$ elements. We obtain some recurrence relations satisfied by Speyer's $g$-polynomials $g_{U_{n,d}}(t)$ of $U_{n,d}$. Based on these recurrence relations, we prove that the polynomial $g_{U_{n,d}}(t)$ has only real zeros for any $n-1\geq d\geq 1$. Furthermore, we show that the coefficient of $g_{U_{n,[n/2]}}(t)$ is asymptotically normal by local and central limit theorems. △ Less

Submitted 27 August, 2024; originally announced August 2024.

Comments: 14 pages

MSC Class: Primary 05B35; 26C10; Secondary 60F05; 62E20

The Degeneracy Loci for Smooth Moduli of Sheaves

Authors: Yu Zhao

Abstract: Let S be a smooth projective surface over $\mathbb{C}$. We prove that, under certain technical assumptions, the degeneracy locus of the universal sheaf over the moduli space of stable sheaves is either empty or an irreducible Cohen-Macaulay variety of the expected dimension. We also provide a criterion for when the degeneracy locus is non-empty. This result generalizes the work of Bayer, Chen, and… ▽ More Let S be a smooth projective surface over $\mathbb{C}$. We prove that, under certain technical assumptions, the degeneracy locus of the universal sheaf over the moduli space of stable sheaves is either empty or an irreducible Cohen-Macaulay variety of the expected dimension. We also provide a criterion for when the degeneracy locus is non-empty. This result generalizes the work of Bayer, Chen, and Jiang for the Hilbert scheme of points on surfaces. The above result is a special case of a general phenomenon: for a perfect complex of Tor-amplitude [0,1], the geometry of the degeneracy locus is closely related to the geometry of the derived Grassmannian. We analyze their birational geometry and relate it to the incidence varieties of derived Grassmannians. As a corollary, we prove a statement previously claimed by the author in arXiv:2408.06860. △ Less

Submitted 26 August, 2024; originally announced August 2024.

Conditional nonparametric variable screening by neural factor regression

Authors: Jianqing Fan, Weining Wang, Yue Zhao

Abstract: High-dimensional covariates often admit linear factor structure. To effectively screen correlated covariates in high-dimension, we propose a conditional variable screening test based on non-parametric regression using neural networks due to their representation power. We ask the question whether individual covariates have additional contributions given the latent factors or more generally a set of… ▽ More High-dimensional covariates often admit linear factor structure. To effectively screen correlated covariates in high-dimension, we propose a conditional variable screening test based on non-parametric regression using neural networks due to their representation power. We ask the question whether individual covariates have additional contributions given the latent factors or more generally a set of variables. Our test statistics are based on the estimated partial derivative of the regression function of the candidate variable for screening and a observable proxy for the latent factors. Hence, our test reveals how much predictors contribute additionally to the non-parametric regression after accounting for the latent factors. Our derivative estimator is the convolution of a deep neural network regression estimator and a smoothing kernel. We demonstrate that when the neural network size diverges with the sample size, unlike estimating the regression function itself, it is necessary to smooth the partial derivative of the neural network estimator to recover the desired convergence rate for the derivative. Moreover, our screening test achieves asymptotic normality under the null after finely centering our test statistics that makes the biases negligible, as well as consistency for local alternatives under mild conditions. We demonstrate the performance of our test in a simulation study and two real world applications. △ Less

Submitted 20 August, 2024; originally announced August 2024.

Comments: 80 pages, 14 figures

MSC Class: 62P20 (Primary); 62G20 (Secondary)

Stability for the $2\frac12$-D compressible viscous non-resistive and heat-conducting magnetohydrodynamic flow

Authors: Xiaoping Zhai, Yongsheng Li, Yajuan Zhao

Abstract: In this paper, we are concerned with the initial boundary values problem associated to the compressible viscous non-resistive and heat-conducting magnetohydrodynamic flow, where the magnetic field is vertical. More precisely, by exploiting the intrinsic structure of the system and introducing several new unknown quantities, we overcome the difficulty stemming from the lack of dissipation for densi… ▽ More In this paper, we are concerned with the initial boundary values problem associated to the compressible viscous non-resistive and heat-conducting magnetohydrodynamic flow, where the magnetic field is vertical. More precisely, by exploiting the intrinsic structure of the system and introducing several new unknown quantities, we overcome the difficulty stemming from the lack of dissipation for density and magnetic field, and prove the global well-posedness of strong solutions in the framework of Soboles spaces $H^3$. In addition, we also get the exponential decay for this non-resistive system. Different from the known results [23], [24], [42], we donot need the assumption that the background magnetic field is positive here. △ Less

Submitted 14 August, 2024; originally announced August 2024.

Comments: 25pages

Hilbert schemes on blowing ups and the free Boson

Authors: Yu Zhao

Abstract: For different cohomology theories (including the Hochschild homology, Hodge cohomology, and Chow groups), we identify the cohomology of moduli space of rank 1 perverse coherent sheaves on the blow-up of a surface by Nakajima-Yoshioka as the tensor product of cohomology of Hilbert schemes with a Fermionic Fock space. As the stable limit, we identify the cohomology of Hilbert schemes of blow-up of a… ▽ More For different cohomology theories (including the Hochschild homology, Hodge cohomology, and Chow groups), we identify the cohomology of moduli space of rank 1 perverse coherent sheaves on the blow-up of a surface by Nakajima-Yoshioka as the tensor product of cohomology of Hilbert schemes with a Fermionic Fock space. As the stable limit, we identify the cohomology of Hilbert schemes of blow-up of a surface as the tensor product of cohomology of Hilbert schemes on the surface with a Bosonic Fock space. The actions of the infinite-dimensional Clifford algebra and Heisenberg algebra are all given by geometric correspondences. △ Less

Submitted 13 August, 2024; originally announced August 2024.

Ultraspherical Spectral Method for Block Copolymer Systems on Unit Disk

Authors: Wangbo Luo, Yanxiang Zhao

Abstract: In this paper, we investigate ultraspherical spectral method for the Ohta-Kawasaki (OK) and Nakazawa-Ohta (NO) models in the disk domain, representing diblock and triblock copolymer systems, respectively. We employ ultraspherical spectral discretization for spatial variables in the disk domain and apply the second-order backward differentiation formula (BDF) method for temporal discretization. To… ▽ More In this paper, we investigate ultraspherical spectral method for the Ohta-Kawasaki (OK) and Nakazawa-Ohta (NO) models in the disk domain, representing diblock and triblock copolymer systems, respectively. We employ ultraspherical spectral discretization for spatial variables in the disk domain and apply the second-order backward differentiation formula (BDF) method for temporal discretization. To our best knowledge, this is the first study to develop a numerical method for diblock and triblock copolymer systems with long-range interactions in disk domains. We show the energy stability of the numerical method in both semi-discrete and fully-discrete discretizations. In our numerical experiments, we verify the second-order temporal convergence rate and the energy stability of the proposed methods. Our numerical results show that the coarsening dynamics in diblock copolymers lead to bubble assemblies both inside and on the boundary of the disk. Additionally, in the triblock copolymer system, we observe several novel pattern formations, including single and double bubble assemblies in the unit disk. These findings are detailed through extensive numerical experiments. △ Less

Submitted 20 July, 2024; originally announced July 2024.

SFC: Achieve Accurate Fast Convolution under Low-precision Arithmetic

Authors: Liulu He, Yufei Zhao, Rui Gao, Yuan Du, Li Du

Abstract: Fast convolution algorithms, including Winograd and FFT, can efficiently accelerate convolution operations in deep models. However, these algorithms depend on high-precision arithmetic to maintain inference accuracy, which conflicts with the model quantization. To resolve this conflict and further improve the efficiency of quantized convolution, we proposes SFC, a new algebra transform for fast co… ▽ More Fast convolution algorithms, including Winograd and FFT, can efficiently accelerate convolution operations in deep models. However, these algorithms depend on high-precision arithmetic to maintain inference accuracy, which conflicts with the model quantization. To resolve this conflict and further improve the efficiency of quantized convolution, we proposes SFC, a new algebra transform for fast convolution by extending the Discrete Fourier Transform (DFT) with symbolic computing, in which only additions are required to perform the transformation at specific transform points, avoiding the calculation of irrational number and reducing the requirement for precision. Additionally, we enhance convolution efficiency by introducing correction terms to convert invalid circular convolution outputs of the Fourier method into effective ones. The numerical error analysis is presented for the first time in this type of work and proves that our algorithms can provide a 3.68x multiplication reduction for 3x3 convolution, while the Winograd algorithm only achieves a 2.25x reduction with similarly low numerical errors. Experiments carried out on benchmarks and FPGA show that our new algorithms can further improve the computation efficiency of quantized models while maintaining accuracy, surpassing both the quantization-alone method and existing works on fast convolution quantization. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Comments: ICML 2024

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

Authors: Yuying Zhao, Xiaojie Wang, Zhongqiang Zhang

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

Submitted 20 June, 2024; originally announced June 2024.

The extended reverse ultra log-concavity of transposed Boros-Moll sequences

Authors: James J. Y. Zhao

Abstract: The Boros-Moll sequences $\{d_\ell(m)\}_{\ell=0}^m$ arise in the study of evaluation of a quartic integral. After the infinite log-concavity conjecture of the sequence $\{d_\ell(m)\}_{\ell=0}^m$ was proposed by Boros and Moll, a lot of interesting inequalities on $d_\ell(m)$ were obtained, although the conjecture is still open. Since $d_\ell(m)$ has two parameters, it is natural to consider the pr… ▽ More The Boros-Moll sequences $\{d_\ell(m)\}_{\ell=0}^m$ arise in the study of evaluation of a quartic integral. After the infinite log-concavity conjecture of the sequence $\{d_\ell(m)\}_{\ell=0}^m$ was proposed by Boros and Moll, a lot of interesting inequalities on $d_\ell(m)$ were obtained, although the conjecture is still open. Since $d_\ell(m)$ has two parameters, it is natural to consider the properties for the sequences $\{d_\ell(m)\}_{m\ge \ell}$, which are called the \emph{transposed Boros-Moll sequences} here. In this paper, we mainly prove the extended reverse ultra log-concavity of the transposed Boros-Moll sequences $\{d_\ell(m)\}_{m\ge \ell}$, and hence give an upper bound for the ratio ${d_\ell^2(m)}/{(d_\ell(m-1)d_\ell(m+1))}$. A lower bound for this ratio is also established which implies a result stronger than the log-concavity of the sequences $\{d_\ell(m)\}_{m\ge \ell}$. As a consequence, we also show that the transposed Boros-Moll sequences possess a stronger log-concave property than the Boros-Moll sequences do. At last, we propose some conjectures on the Boros-Moll sequences and their transposes. △ Less

Submitted 19 June, 2024; originally announced June 2024.

Comments: 15 pages

MSC Class: 05A20; 11B83

Kato's chaos of multiple mappings and its continuous self-maps

Authors: Yingcui Zhao

Abstract: In 2016, Hou and Wang introduced the concept of multiple mappings based on iterated function system, which is an important branch of fractal theory. In this paper, we introduce the definitions of sensitivity, accessibility, and Kato's chaos of multiple mappings from a set-valued perspective. We show that multiple mappings and its continuous self-maps do not imply each other in terms of sensitivity… ▽ More In 2016, Hou and Wang introduced the concept of multiple mappings based on iterated function system, which is an important branch of fractal theory. In this paper, we introduce the definitions of sensitivity, accessibility, and Kato's chaos of multiple mappings from a set-valued perspective. We show that multiple mappings and its continuous self-maps do not imply each other in terms of sensitivity and accessibility. While a sufficient condition for multiple mappings to be sensitive, accessible and Kato's chaotic is provided, respectively. And the sensitivity, accessibility, and Kato's chaos of multiple mappings are preserved under topological conjugation. △ Less

Submitted 14 March, 2024; originally announced June 2024.

Comments: arXiv admin note: substantial text overlap with arXiv:2311.03715, arXiv:2311.02932

Multilinear Fractional Integral Operators with Generalized Kernels

Authors: Yan Lin, Yuhang Zhao, Shuhui Yang

Abstract: In this article, we introduce a class of multilinear fractional integral operators with generalized kernels that are weaker than the Dini kernel condition. We establish the boundedness of multilinear fractional integral operators with generalized kernels on weighted Lebesgue spaces and variable exponent Lebesgue spaces, as well as the boundedness of multilinear commutators generated by multilinear… ▽ More In this article, we introduce a class of multilinear fractional integral operators with generalized kernels that are weaker than the Dini kernel condition. We establish the boundedness of multilinear fractional integral operators with generalized kernels on weighted Lebesgue spaces and variable exponent Lebesgue spaces, as well as the boundedness of multilinear commutators generated by multilinear fractional integral operators with generalized kernels and $BMO$ functions. Even when the generalized kernels condition goes back to the Dini kernel condition, the conclusions on the commutators remain new. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 20 pages

MSC Class: 42B25; 26A33; 42B35

Synchronized Optimal Transport for Joint Modeling of Dynamics Across Multiple Spaces

Authors: Zixuan Cang, Yanxiang Zhao

Abstract: Optimal transport has been an essential tool for reconstructing dynamics from complex data. With the increasingly available multifaceted data, a system can often be characterized across multiple spaces. Therefore, it is crucial to maintain coherence in the dynamics across these diverse spaces. To address this challenge, we introduce Synchronized Optimal Transport (SyncOT), a novel approach to join… ▽ More Optimal transport has been an essential tool for reconstructing dynamics from complex data. With the increasingly available multifaceted data, a system can often be characterized across multiple spaces. Therefore, it is crucial to maintain coherence in the dynamics across these diverse spaces. To address this challenge, we introduce Synchronized Optimal Transport (SyncOT), a novel approach to jointly model dynamics that represent the same system through multiple spaces. With given correspondence between the spaces, SyncOT minimizes the aggregated cost of the dynamics induced across all considered spaces. The problem is discretized into a finite-dimensional convex problem using a staggered grid. Primal-dual algorithm-based approaches are then developed to solve the discretized problem. Various numerical experiments demonstrate the capabilities and properties of SyncOT and validate the effectiveness of the proposed algorithms. △ Less

Submitted 5 June, 2024; originally announced June 2024.

Muckenhoupt Weights Meet Brezis--Seeger--Van Schaftingen--Yung Formulae in Ball Banach Function Spaces

Authors: Yinqin Li, Dachun Yang, Wen Yuan, Yangyang Zhang, Yirui Zhao

Abstract: In this article, via first establishing a weighted variant of the profound and far-reaching inequality obtained by A. Cohen, W. Dahmen, I. Daubechies, and R. DeVore in 2003, the authors give two new characterizations of Muckenhoupt weights. As an application, the authors further establish a representation formula of gradients with sharp parameters in ball Banach function spaces, which extends the… ▽ More In this article, via first establishing a weighted variant of the profound and far-reaching inequality obtained by A. Cohen, W. Dahmen, I. Daubechies, and R. DeVore in 2003, the authors give two new characterizations of Muckenhoupt weights. As an application, the authors further establish a representation formula of gradients with sharp parameters in ball Banach function spaces, which extends the famous formula obtained by H. Brezis, A. Seeger, J. Van Schaftingen, and P.-L. Yung in 2021 from classical Sobolev spaces to various different Sobolev-type spaces and gives an affirmative answer to the question in page 29 of [Calc. Var. Partial Differential Equations 62 (2023), Paper No. 234]. The most novelty of this article exists in subtly revealing the mutual equivalences among the Muckenhoupt weight, the weighted variant of the inequality of Cohen et al., and the weighted upper estimate of the formula of Brezis et al. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Comments: 48 pages; Submitted

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

Turán number of complete multipartite graphs in multipartite graphs

Authors: Jie Han, Yi Zhao

Abstract: In this paper we study a multi-partite version of the Erdős--Stone theorem. Given integers $r<k$ and $t\ge 1$, let $ex_k(n, K_{r+1}(t))$ be the maximum number of edges of $K_{r+1}(t)$-free $k$-partite graphs with $n$ vertices in each part, where $K_{r+1}(t)$ is the $t$-blowup of $K_{r+1}$. An easy consequence of the supersaturaion result gives that $ex_k(n, K_{r+1}(t)) = ex_k(n, K_{r+1})+o(n^2)$.… ▽ More In this paper we study a multi-partite version of the Erdős--Stone theorem. Given integers $r<k$ and $t\ge 1$, let $ex_k(n, K_{r+1}(t))$ be the maximum number of edges of $K_{r+1}(t)$-free $k$-partite graphs with $n$ vertices in each part, where $K_{r+1}(t)$ is the $t$-blowup of $K_{r+1}$. An easy consequence of the supersaturaion result gives that $ex_k(n, K_{r+1}(t)) = ex_k(n, K_{r+1})+o(n^2)$. Similar to a result of Erd\H os and Simonovits for the non-partite case, we find that the error term is closely related to the (multi-partite) Zarankiewicz problem. Using such Zarankiewicz numbers, for $t=2,3$, we determine the error term up to an additive linear term; using some natural assumptions on such Zarankiewicz numbers, we determine the error term up to an additive constant depending on $k$, $r$ and $t$. We actually obtain exact results in many cases, for example, when $k\equiv 0, 1 \pmod{r}$. Our proof uses the stability method and starts by proving a stability result for $K_{r+1}$-free multi-partite graphs. △ Less

Submitted 26 May, 2024; originally announced May 2024.

Comments: 33 pages

The number of cliques in hypergraphs with forbidden subgraphs

Authors: Ayush Basu, Vojtech Rodl, Yi Zhao

Abstract: We study the maximum number of $r$-vertex cliques in $(r-1)$-uniform hypergraphs not containing complete $r$-partite hypergraphs $K_r^{(r-1)}(a_1, \dots, a_r)$. By using the hypergraph removal lemma, we show that this maximum is $o( n^{r - 1/(a_1 \cdots a_{r-1})} )$. This immediately implies the corresponding results of Mubayi and Mukherjee and of Balogh, Jiang, and Luo for graphs. We also provide… ▽ More We study the maximum number of $r$-vertex cliques in $(r-1)$-uniform hypergraphs not containing complete $r$-partite hypergraphs $K_r^{(r-1)}(a_1, \dots, a_r)$. By using the hypergraph removal lemma, we show that this maximum is $o( n^{r - 1/(a_1 \cdots a_{r-1})} )$. This immediately implies the corresponding results of Mubayi and Mukherjee and of Balogh, Jiang, and Luo for graphs. We also provide a lower bound by using hypergraph Turán numbers. △ Less

Submitted 13 May, 2024; originally announced May 2024.

Exponentially Weighted Algorithm for Online Network Resource Allocation with Long-Term Constraints

Authors: Ahmed Sid-Ali, Ioannis Lambadaris, Yiqiang Q. Zhao, Gennady Shaikhet, Amirhossein Asgharnia

Abstract: This paper studies an online optimal resource reservation problem in communication networks with job transfers where the goal is to minimize the reservation cost while maintaining the blocking cost under a certain budget limit. To tackle this problem, we propose a novel algorithm based on a randomized exponentially weighted method that encompasses long-term constraints. We then analyze the perform… ▽ More This paper studies an online optimal resource reservation problem in communication networks with job transfers where the goal is to minimize the reservation cost while maintaining the blocking cost under a certain budget limit. To tackle this problem, we propose a novel algorithm based on a randomized exponentially weighted method that encompasses long-term constraints. We then analyze the performance of our algorithm by establishing an upper bound for the associated regret and the cumulative constraint violations. Finally, we present numerical experiments where we compare the performance of our algorithm with those of reinforcement learning where we show that our algorithm surpasses it. △ Less

Submitted 3 May, 2024; originally announced May 2024.

Comments: arXiv admin note: text overlap with arXiv:2305.15558

Analysing 2-$(v,k,2)$ designs admitting a flag-transitive almost simple automorphism group with socle $PSL(2,q)$ by means of conics and hyperovals of $PG(2,q)$

Authors: Alessandro Montinaro, Yanwei Zhao, Zhilin Zhang, Shenglin Zhou

Abstract: The classification of the $2$-designs with $λ=2$ admitting a flag-transitive automorphism groups with socle $PSL(2,q)$ is completed by settling the two open cases in \cite{ABDT}. The result is achieved by using conics and hyperovals of $PG(2,q)$. The classification of the $2$-designs with $λ=2$ admitting a flag-transitive automorphism groups with socle $PSL(2,q)$ is completed by settling the two open cases in \cite{ABDT}. The result is achieved by using conics and hyperovals of $PG(2,q)$. △ Less

Submitted 29 April, 2024; originally announced April 2024.

MSC Class: 05B05; 05B25; 20B25; 51E15; 51E21

Efficient bound preserving and asymptotic preserving semi-implicit schemes for the fast reaction-diffusion system

Authors: Yu Zhao, Zhennan Zhou

Abstract: We consider a special type of fast reaction-diffusion systems in which the coefficients of the reaction terms of the two substances are much larger than those of the diffusion terms while the diffusive motion to the substrate is negligible. Specifically speaking, the rate constants of the reaction terms are $O(1/ε)$ while the diffusion coefficients are $O(1)$ where the parameter $ε$ is small. When… ▽ More We consider a special type of fast reaction-diffusion systems in which the coefficients of the reaction terms of the two substances are much larger than those of the diffusion terms while the diffusive motion to the substrate is negligible. Specifically speaking, the rate constants of the reaction terms are $O(1/ε)$ while the diffusion coefficients are $O(1)$ where the parameter $ε$ is small. When the rate constants of the reaction terms become highly large, i.e. $ε$ tends to 0, the singular limit behavior of such a fast reaction-diffusion system is inscribed by the Stefan problem with latent heat, which brings great challenges in numerical simulations. In this paper, we adopt a semi-implicit scheme, which is first-order accurate in time and can accurately approximate the interface propagation even when the reaction becomes extremely fast, that is to say, the parameter $ε$ is sufficiently small. The scheme satisfies the positivity, bound preserving properties and has $L^2$ stability and the linearized stability results of the system. For better performance on numerical simulations, we then construct a semi-implicit Runge-Kutta scheme which is second-order accurate in time. Numerous numerical tests are carried out to demonstrate the properties, such as the order of accuracy, positivity and bound preserving, the capturing of the sharp interface with various $ε$ and to simulate the dynamics of the substances and the substrate, and to explore the heat transfer process, such as solid melting or liquid solidification in two dimensions. △ Less

Submitted 29 April, 2024; originally announced April 2024.

MSC Class: 35Q92; 65M06; 65M12; 92E20

Optimal gradient estimates for the insulated conductivity problem with general convex inclusions case

Authors: Haigang Li, Yan Zhao

Abstract: We study the insulated conductivity problem which involves two adjacent convex insulators embedded in a bounded domain. It is known that the gradient of solutions may blow up as the distance between the two inclusions tends to zero. However, the sharpness of the blow up rate for general convex insulator case in dimension $n\geq3$ has remained open. The novelty of this paper is that we answer this… ▽ More We study the insulated conductivity problem which involves two adjacent convex insulators embedded in a bounded domain. It is known that the gradient of solutions may blow up as the distance between the two inclusions tends to zero. However, the sharpness of the blow up rate for general convex insulator case in dimension $n\geq3$ has remained open. The novelty of this paper is that we answer this problem affirmatively by establishing a pointwise upper bound of the gradient for general convex insulators, along with a corresponding lower bound that achieves optimal blow up rates. These rates are associated with the first nonzero eigenvalue of an elliptic operator determined by the geometry of insulators. Our results improve and make complete the previous result for ball insulators case studied in \cite{DLY}. △ Less

Submitted 28 August, 2024; v1 submitted 26 April, 2024; originally announced April 2024.

Comments: arXiv admin note: text overlap with arXiv:2203.10081 by other authors

Uncommon linear systems of two equations

Authors: Dingding Dong, Anqi Li, Yufei Zhao

Abstract: A system of linear equations $L$ is common over $\mathbb{F}_p$ if, as $n\to\infty$, any 2-coloring of $\mathbb{F}_p^n$ gives asymptotically at least as many monochromatic solutions to $L$ as a random 2-coloring. The notion of common linear systems is analogous to that of common graphs, i.e., graphs whose monochromatic density in 2-edge-coloring of cliques is asymptotically minimized by the random… ▽ More A system of linear equations $L$ is common over $\mathbb{F}_p$ if, as $n\to\infty$, any 2-coloring of $\mathbb{F}_p^n$ gives asymptotically at least as many monochromatic solutions to $L$ as a random 2-coloring. The notion of common linear systems is analogous to that of common graphs, i.e., graphs whose monochromatic density in 2-edge-coloring of cliques is asymptotically minimized by the random coloring. Saad and Wolf initiated a systematic study on identifying common linear systems, built upon the earlier work of Cameron-Cilleruelo-Serra. When $L$ is a single equation, Fox-Pham-Zhao gave a complete characterization of common linear equations. When $L$ consists of two equations, Kamčev-Liebenau-Morrison showed that irredundant $2\times 4$ linear systems are always uncommon. In this work, (1) we determine commonness of all $2\times 5$ linear systems up to a small number of cases, and (2) we show that all $2\times k$ linear systems with $k$ even and girth (minimum number of nonzero coefficients of a nonzero equation spanned by the system) $k-1$ are uncommon, answering a question of Kamčev-Liebenau-Morrison. △ Less

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

Comments: 59 pages, 1 figure

On length-preserving and area-preserving inverse curvature flow of planar curves with singularities

Authors: Yunlong Yang, Yanwen Zhao, Jianbo Fang, Yanlong Zhang

Abstract: This paper aims to investigate the evolution problem for planar curves with singularities. Motivated by the inverse curvature flow introduced by Li and Wang (Calc. Var. Partial Differ. Equ. 62 (2023), No. 135), we intend to consider the area-preserving and length-preserving inverse curvature flow with nonlocal term for $\ell$-convex Legendre curves. For the area-preserving flow, an $\ell$-convex L… ▽ More This paper aims to investigate the evolution problem for planar curves with singularities. Motivated by the inverse curvature flow introduced by Li and Wang (Calc. Var. Partial Differ. Equ. 62 (2023), No. 135), we intend to consider the area-preserving and length-preserving inverse curvature flow with nonlocal term for $\ell$-convex Legendre curves. For the area-preserving flow, an $\ell$-convex Legendre curve %of with initial algebraic area $A_0>0$ evolves to a circle of radius $\sqrt{\frac{A_0}π}$. For the length-preserving flow, an $\ell$-convex Legendre curve %of with initial algebraic length $L_0$ evolves to a circle of radius $\frac{L_0}{2π}$. As the by-product, we obtain some geometric inequalities for $\ell$-convex Legendre curves through the length-preserving flow. △ Less

Submitted 17 April, 2024; originally announced April 2024.

MSC Class: 53E99; 53A04

Covariance Regression with High-Dimensional Predictors

Authors: Yuheng He, Changliang Zou, Yi Zhao

Abstract: In the high-dimensional landscape, addressing the challenges of covariance regression with high-dimensional covariates has posed difficulties for conventional methodologies. This paper addresses these hurdles by presenting a novel approach for high-dimensional inference with covariance matrix outcomes. The proposed methodology is illustrated through its application in elucidating brain coactivatio… ▽ More In the high-dimensional landscape, addressing the challenges of covariance regression with high-dimensional covariates has posed difficulties for conventional methodologies. This paper addresses these hurdles by presenting a novel approach for high-dimensional inference with covariance matrix outcomes. The proposed methodology is illustrated through its application in elucidating brain coactivation patterns observed in functional magnetic resonance imaging (fMRI) experiments and unraveling complex associations within anatomical connections between brain regions identified through diffusion tensor imaging (DTI). In the pursuit of dependable statistical inference, we introduce an integrative approach based on penalized estimation. This approach combines data splitting, variable selection, aggregation of low-dimensional estimators, and robust variance estimation. It enables the construction of reliable confidence intervals for covariate coefficients, supported by theoretical confidence levels under specified conditions, where asymptotic distributions are provided. Through various types of simulation studies, the proposed approach performs well for covariance regression in the presence of high-dimensional covariates. This innovative approach is applied to the Lifespan Human Connectome Project (HCP) Aging Study, which aims to uncover a typical aging trajectory and variations in the brain connectome among mature and older adults. The proposed approach effectively identifies brain networks and associated predictors of white matter integrity, aligning with established knowledge of the human brain. △ Less

Submitted 9 April, 2024; originally announced April 2024.

Identification of High-dimensional ARMA Models with Binary-Valued Observations

Authors: Xin Li, Ting Wang, Jin Guo, Yanlong Zhao

Abstract: This paper studies system identification of high-dimensional ARMA models with binary-valued observations. Compared with existing quantized identification of ARMA models, this problem is more challenging since the accessible information is much less. Different from the identification of FIR models with binary-valued observations, the prediction of original system output and the parameter both need… ▽ More This paper studies system identification of high-dimensional ARMA models with binary-valued observations. Compared with existing quantized identification of ARMA models, this problem is more challenging since the accessible information is much less. Different from the identification of FIR models with binary-valued observations, the prediction of original system output and the parameter both need to be estimated in ARMA models. We propose an online identification algorithm consisting of parameter estimation and prediction of original system output. The parameter estimation and the prediction of original output are strongly coupled but mutually reinforcing. By analyzing the two estimates at the same time instead of analyzing separately, we finally prove that the parameter estimate can converge to the true parameter with convergence rate $O(1/k)$ under certain conditions. Simulations are given to demonstrate the theoretical results. △ Less

Submitted 28 May, 2024; v1 submitted 1 April, 2024; originally announced April 2024.

Stability for a multi-frequency inverse random source problem

Authors: Tianjiao Wang, Xiang Xu, Yue Zhao

Abstract: We present stability estimates for the inverse source problem of the stochastic Helmholtz equation in two and three dimensions by either near-field or far-field data. The random source is assumed to be a microlocally isotropic generalized Gaussian random function. For the direct problem, by exploring the regularity of the Green function, we demonstrate that the direct problem admits a unique bound… ▽ More We present stability estimates for the inverse source problem of the stochastic Helmholtz equation in two and three dimensions by either near-field or far-field data. The random source is assumed to be a microlocally isotropic generalized Gaussian random function. For the direct problem, by exploring the regularity of the Green function, we demonstrate that the direct problem admits a unique bounded solution with an explicit integral representation, which enhances the existing regularity result. For the case using near-field data, the analysis of the inverse problem employs microlocal analysis to achieve an estimate for the Fourier transform of the micro-correlation strength by the near-field correlation data and a high-frequency tail. The stability follows by showing the analyticity of the data and applying a novel analytic continuation principle. The stability estimate by far-field data is derived by investigating the correlation of the far-field data. The stability estimate consists of the Lipschitz type data discrepancy and the logarithmic stability. The latter decreases as the upper bound of the frequency increases, which exhibits the phenomenon of increasing stability. △ Less

Submitted 19 March, 2024; originally announced March 2024.

Comments: 19pages

MSC Class: 35Q74; 35R30; 78A46

Global Solutions and Asymptotic Behavior for the Three-dimensional Viscous Non-resistive MHD System with Some Large Perturbations

Authors: Youyi Zhao

Abstract: We revisit the global existence of solutions with some large perturbations to the incompressible, viscous, and non-resistive MHD system in a three-dimensional periodic domain, where the impressed magnetic field satisfies the Diophantine condition, and the intensity of the impressed magnetic field, denoted by $m$, is large compared to the perturbations. It was proved by Jiang--Jiang that the highes… ▽ More We revisit the global existence of solutions with some large perturbations to the incompressible, viscous, and non-resistive MHD system in a three-dimensional periodic domain, where the impressed magnetic field satisfies the Diophantine condition, and the intensity of the impressed magnetic field, denoted by $m$, is large compared to the perturbations. It was proved by Jiang--Jiang that the highest-order derivatives of the velocity increase with $m$, and the convergence rate of the nonlinear system towards a linearized problem is of $m^{-1/2}$ in [F. Jiang and S. Jiang, Arch. Ration. Mech. Anal., 247 (2023), 96]. In this paper, we adopt a different approach by leveraging vorticity estimates to establish the highest-order energy estimate. This strategy prevents the appearance of terms that grow with $m$, and thus the increasing behavior of the highest-order derivatives of the velocity with respect to $m$ does not appear. Additionally, we use the vorticity estimate to demonstrate the convergence rate of the nonlinear system towards a linearized problem as time or $m$ approaches infinity. Notably, our analysis reveals that the convergence rate in $m$ is faster compared to the finding of Jiang--Jiang. Finally, a key contribution of our work is the identification of an integrable time-decay of the lower dissipation, which can replace the time-decay of lower energy in closing the highest-order energy estimate. This finding significantly relaxes the regularity requirements for the initial perturbations. △ Less

Submitted 8 March, 2024; originally announced March 2024.

Analysis of a Leslie-Gower model with Alle effects, cooperative hunting, and constant placement rates

Authors: Yonghui Zhao

Abstract: This paper investigates the dynamical properties of the Leslie-Gower model with Alle effects, cooperative hunting, and constant placement rates. The conditions for the existence of the triple equilibrium point of the model are first analyzed. Subsequently, the canonical type theory and the qualitative theory of planar systems are applied to obtain that the triple equilibrium point can be a node wi… ▽ More This paper investigates the dynamical properties of the Leslie-Gower model with Alle effects, cooperative hunting, and constant placement rates. The conditions for the existence of the triple equilibrium point of the model are first analyzed. Subsequently, the canonical type theory and the qualitative theory of planar systems are applied to obtain that the triple equilibrium point can be a node with a residual dimension of 2 and an equilibrium point with a residual dimension of 3 under different parameter conditions. Finally, it is proved that the system bifurcates with a residual dimension of 2 in the vicinity of the node with cooperative hunting and placement rate as branching parameters. △ Less

Submitted 7 March, 2024; originally announced March 2024.

On the convergence of conditional gradient method for unbounded multiobjective optimization problems

Authors: Wang Chen, Yong Zhao, Liping Tang, Xinmin Yang

Abstract: This paper focuses on developing a conditional gradient algorithm for multiobjective optimization problems with an unbounded feasible region. We employ the concept of recession cone to establish the well-defined nature of the algorithm. The asymptotic convergence property and the iteration-complexity bound are established under mild assumptions. Numerical examples are provided to verify the algori… ▽ More This paper focuses on developing a conditional gradient algorithm for multiobjective optimization problems with an unbounded feasible region. We employ the concept of recession cone to establish the well-defined nature of the algorithm. The asymptotic convergence property and the iteration-complexity bound are established under mild assumptions. Numerical examples are provided to verify the algorithmic performance. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Direct and Inverse scattering in a three-dimensional planar waveguide

Authors: Yan Chang, Yukun Guo, Yue Zhao

Abstract: In this paper, we study the direct and inverse scattering of the Schrödinger equation in a three-dimensional planar waveguide. For the direct problem, we derive a resonance-free region and resolvent estimates for the resolvent of the Schrödinger operator in such a geometry. Based on the analysis of the resolvent, several inverse problems are investigated. First, given the potential function, we pr… ▽ More In this paper, we study the direct and inverse scattering of the Schrödinger equation in a three-dimensional planar waveguide. For the direct problem, we derive a resonance-free region and resolvent estimates for the resolvent of the Schrödinger operator in such a geometry. Based on the analysis of the resolvent, several inverse problems are investigated. First, given the potential function, we prove the uniqueness of the inverse source problem with multi-frequency data. We also develop a Fourier-based method to reconstruct the source function. The capability of this method is numerically illustrated by examples. Second, the uniqueness and increased stability of an inverse potential problem from data generated by incident waves are achieved in the absence of the source function. To derive the stability estimate, we use an argument of quantitative analytic continuation in complex theory. Third, we prove the uniqueness of simultaneously determining the source and potential by active boundary data generated by incident waves. In these inverse problems, we only use the limited lateral Dirichlet boundary data at multiple wavenumbers within a finite interval. △ Less

Submitted 24 February, 2024; originally announced February 2024.

Comments: 22 pages, 7 figures

Small positive values and limit theorems for supercritical branching processes with immigration in random environment

Authors: Yinxuan Zhao, Mei Zhang

Abstract: Let $(Z_n)$ be a supercritical branching process with immigration in a random environment. The small positive values and some lower deviation inequalities for $Z$ are investigated. Based on these results, the central limit theorem of $\log Z_n$ and the Edgeworth expansion are obtained. The study is taken under the assumption that each individual produces $0$ offspring with a positive probability. Let $(Z_n)$ be a supercritical branching process with immigration in a random environment. The small positive values and some lower deviation inequalities for $Z$ are investigated. Based on these results, the central limit theorem of $\log Z_n$ and the Edgeworth expansion are obtained. The study is taken under the assumption that each individual produces $0$ offspring with a positive probability. △ Less

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

The Briggs inequality of Boros-Moll sequences

Authors: Zhong-Xue Zhang, James Jing Yu Zhao

Abstract: Briggs conjectured that if a polynomial $a_0+a_1x+\cdots+a_nx^n$ with real coefficients has only negative zeros, then $$a^2_k(a^2_k - a_{k-1}a_{k+1}) > a^2_{k-1}(a^2_{k+1} - a_ka_{k+2})$$ for any $1\leq k\leq n-1$. The Boros-Moll sequence $\{d_i(m)\}_{i=0}^m$ arises in the study of evaluation of certain quartic integral, and a lot of interesting inequalities for this sequence have been obtained. I… ▽ More Briggs conjectured that if a polynomial $a_0+a_1x+\cdots+a_nx^n$ with real coefficients has only negative zeros, then $$a^2_k(a^2_k - a_{k-1}a_{k+1}) > a^2_{k-1}(a^2_{k+1} - a_ka_{k+2})$$ for any $1\leq k\leq n-1$. The Boros-Moll sequence $\{d_i(m)\}_{i=0}^m$ arises in the study of evaluation of certain quartic integral, and a lot of interesting inequalities for this sequence have been obtained. In this paper we show that the Boros-Moll sequence $\{d_i(m)\}_{i=0}^m$, its normalization $\{d_i(m)/i!\}_{i=0}^m$, and its transpose $\{d_i(m)\}_{m\ge i}$ satisfy the Briggs inequality. For the first two sequences, we prove the Briggs inequality by using a lower bound for $(d_{i-1}(m)d_{i+1}(m))/d_i^2(m)$ due to Chen and Gu and an upper bound due to Zhao. For the transposed sequence, we derive the Briggs inequality by establishing its strict ratio-log-convexity. As a consequence, we also obtain the strict log-convexity of the sequence $\{\sqrt[n]{d_i(i+n)}\}_{n\ge 1}$ for $i\ge 1$. △ Less

Submitted 18 February, 2024; originally announced February 2024.

Comments: 18 pages

MSC Class: 05A20; 11B83

Conformal Predictive Programming for Chance Constrained Optimization

Authors: Yiqi Zhao, Xinyi Yu, Jyotirmoy V. Deshmukh, Lars Lindemann

Abstract: Motivated by the advances in conformal prediction (CP), we propose conformal predictive programming (CPP), an approach to solve chance constrained optimization (CCO) problems, i.e., optimization problems with nonlinear constraint functions affected by arbitrary random parameters. CPP utilizes samples from these random parameters along with the quantile lemma -- which is central to CP -- to transfo… ▽ More Motivated by the advances in conformal prediction (CP), we propose conformal predictive programming (CPP), an approach to solve chance constrained optimization (CCO) problems, i.e., optimization problems with nonlinear constraint functions affected by arbitrary random parameters. CPP utilizes samples from these random parameters along with the quantile lemma -- which is central to CP -- to transform the CCO problem into a deterministic optimization problem. We then present two tractable reformulations of CPP by: (1) writing the quantile as a linear program along with its KKT conditions (CPP-KKT), and (2) using mixed integer programming (CPP-MIP). CPP comes with marginal probabilistic feasibility guarantees for the CCO problem that are conceptually different from existing approaches, e.g., the sample approximation and the scenario approach. While we explore algorithmic similarities with the sample approximation approach, we emphasize that the strength of CPP is that it can easily be extended to incorporate different variants of CP. To illustrate this, we present robust conformal predictive programming to deal with distribution shifts in the uncertain parameters of the CCO problem. △ Less

Submitted 11 February, 2024; originally announced February 2024.

Improved upper bounds for wide-sense frameproof codes

Authors: Yuhao Zhao, Xiande Zhang

Abstract: Frameproof codes have been extensively studied for many years due to their application in copyright protection and their connection to extremal set theory. In this paper, we investigate upper bounds on the cardinality of wide-sense $t$-frameproof codes. For $t=2$, we apply results from Sperner theory to give a better upper bound, which significantly improves a recent bound by Zhou and Zhou. For… ▽ More Frameproof codes have been extensively studied for many years due to their application in copyright protection and their connection to extremal set theory. In this paper, we investigate upper bounds on the cardinality of wide-sense $t$-frameproof codes. For $t=2$, we apply results from Sperner theory to give a better upper bound, which significantly improves a recent bound by Zhou and Zhou. For $t\geq 3$, we provide a general upper bound by establishing a relation between wide-sense frameproof codes and cover-free families. Finally, when the code length $n$ is at most $\frac{15+\sqrt{33}}{24}(t-1)^2$, we show that a wide-sense $t$-frameproof code has at most $n$ codewords, and the unique optimal code consists of all weight-one codewords. As byproducts, our results improve several best known results on binary $t$-frameproof codes. △ Less

Submitted 8 February, 2024; originally announced February 2024.

Nonconvex weighted variational metal artifacts removal via convergent primal-dual algorithms

Authors: Lianfang Wang, Zhangling Chen, Zhifang Liu, Yutong Li, Yunsong Zhao, Hongwei Li, Huibin Chang

Abstract: Direct reconstruction through filtered back projection engenders metal artifacts in polychromatic computed tomography images, attributed to highly attenuating implants, which further poses great challenges for subsequent image analysis. Inpainting the metal trace directly in the Radon domain for the extant variational method leads to strong edge diffusion and potential inherent artifacts. With nor… ▽ More Direct reconstruction through filtered back projection engenders metal artifacts in polychromatic computed tomography images, attributed to highly attenuating implants, which further poses great challenges for subsequent image analysis. Inpainting the metal trace directly in the Radon domain for the extant variational method leads to strong edge diffusion and potential inherent artifacts. With normalization based on pre-segmentation, the inpainted outcome can be notably ameliorated. However, its reconstructive fidelity is heavily contingent on the precision of the presegmentation, and highly accurate segmentation of images with metal artifacts is non-trivial in actuality. In this paper, we propose a nonconvex weighted variational approach for metal artifact reduction. Specifically, in lieu of employing a binary function with zeros in the metal trace, an adaptive weight function is designed in the Radon domain, with zeros in the overlapping regions of multiple disjoint metals as well as areas of highly attenuated projections, and the inverse square root of the measured projection in other regions. A nonconvex L1-alpha L2 regularization term is incorporated to further enhance edge contrast, alongside a box-constraint in the image domain. Efficient first-order primal-dual algorithms, proven to be globally convergent and of low computational cost owing to the closed-form solution of all subproblems, are devised to resolve such a constrained nonconvex model. Both simulated and real experiments are conducted with comparisons to other variational algorithms, validating the superiority of the presented method. Especially in comparison to Reweighted JSR, our proposed algorithm can curtail the total computational cost to at most one-third, and for the case of inaccurate pre-segmentation, the recovery outcomes by the proposed algorithms are notably enhanced. △ Less

Submitted 2 February, 2024; originally announced February 2024.

A parallel preconditioner for the all-at-once linear system from evolutionary PDEs with Crank-Nicolson discretization

Authors: Yong-Liang Zhao, Xian-Ming Gu, Cornelis W. Oosterlee

Abstract: The Crank-Nicolson (CN) method is a well-known time integrator for evolutionary partial differential equations (PDEs) arising in many real-world applications. Since the solution at any time depends on the solution at previous time steps, the CN method is inherently difficult to parallelize. In this paper, we consider a parallel method for the solution of evolutionary PDEs with the CN scheme. Using… ▽ More The Crank-Nicolson (CN) method is a well-known time integrator for evolutionary partial differential equations (PDEs) arising in many real-world applications. Since the solution at any time depends on the solution at previous time steps, the CN method is inherently difficult to parallelize. In this paper, we consider a parallel method for the solution of evolutionary PDEs with the CN scheme. Using an all-at-once approach, we can solve for all time steps simultaneously using a parallelizable over time preconditioner within a standard iterative method. Due to the diagonalization of the proposed preconditioner, we can prove that most eigenvalues of preconditioned matrices are equal to 1 and the others lie in the set: $\left\{z\in\mathbb{C}: 1/(1 + α) < |z| < 1/(1 - α)~{\rm and}~\Re{\rm e}(z) > 0\right\}$, where $0 < α< 1$ is a free parameter. Besides, the efficient implementation of the proposed preconditioner is described. Given certain conditions, we prove that the preconditioned GMRES method exhibits a mesh-independent convergence rate. Finally, we will verify both theoretical findings and the efficacy of the proposed preconditioner via numerical experiments on financial option pricing PDEs. △ Less

Submitted 11 February, 2024; v1 submitted 29 January, 2024; originally announced January 2024.

Comments: 18 pages, 5 figures and 4 tables (update some contexts)

MSC Class: 65L05; 65N22; 65F10

Time-harmonic elastic scattering by unbounded deterministic and random rough surfaces in three dimensions

Authors: Guanghui Hu, Tianjiao Wang, Xiang Xu, Yue Zhao

Abstract: In this paper, we investigate well-posedness of time-harmonic scattering of elastic waves by unbounded rigid rough surfaces in three dimensions. The elastic scattering is caused by an $L^2$ function with a compact support in the $x_3$-direction, and both deterministic and random surfaces are investigated via the variational approach. The rough surface in a deterministic setting is assumed to be Li… ▽ More In this paper, we investigate well-posedness of time-harmonic scattering of elastic waves by unbounded rigid rough surfaces in three dimensions. The elastic scattering is caused by an $L^2$ function with a compact support in the $x_3$-direction, and both deterministic and random surfaces are investigated via the variational approach. The rough surface in a deterministic setting is assumed to be Lipschitz and lie within a finite distance of a flat plane, and the scattering is caused by an inhomogeneous term in the elastic wave equation whose support lies within some finite distance of the boundary. For the deterministic case, a stability estimate of elastic scattering by rough surface is shown at an arbitrary frequency. It is noticed that all constants in {\it a priori} bounds are bounded by explicit functions of the frequency and geometry of rough surfaces. Furthermore, based on this explicit dependence on the frequency together with the measurability and $\mathbb{P}$-essentially separability of the randomness, we obtain a similar bound for the solution of the scattering by random surfaces. △ Less

Submitted 28 January, 2024; originally announced January 2024.

MSC Class: 35A15; 35P25; 74J20

Spectral non-concentration near the top for unimodular random graphs

Authors: Mikolaj Fraczyk, Ben Hayes, Madhu Sudan, Yufei Zhao

Abstract: In recent work on equiangular lines, Jiang, Tidor, Yuan, Zhang, and Zhao showed that a connected bounded degree graph has sublinear second eigenvalue multiplicity. More generally they show that there cannot be too many eigenvalues near the top of the spectrum. We extend this result to infinite unimodular random graphs. As a corollary, the spectral distribution of the adjacency operator cannot have… ▽ More In recent work on equiangular lines, Jiang, Tidor, Yuan, Zhang, and Zhao showed that a connected bounded degree graph has sublinear second eigenvalue multiplicity. More generally they show that there cannot be too many eigenvalues near the top of the spectrum. We extend this result to infinite unimodular random graphs. As a corollary, the spectral distribution of the adjacency operator cannot have an atom at the top. For an infinite regular expander, we deduce that the singularity of the spectral measure at the top satisfies $μ_G[(1-θ)ρ,ρ] \lesssim θ^c$ for some constant $c>0$, where $ρ$ is the spectral radius of the adjacency operator of the graph. This implies new general estimates on the return probabilities of random walks. △ Less

Submitted 13 January, 2024; originally announced January 2024.

Comments: 19 pages

Supervised Gromov-Wasserstein Optimal Transport

Authors: Zixuan Cang, Yaqi Wu, Yanxiang Zhao

Abstract: We introduce the supervised Gromov-Wasserstein (sGW) optimal transport, an extension of Gromov-Wasserstein by incorporating potential infinity patterns in the cost tensor. sGW enables the enforcement of application-induced constraints such as the preservation of pairwise distances by implementing the constraints as an infinity pattern. A numerical solver is proposed for the sGW problem and the eff… ▽ More We introduce the supervised Gromov-Wasserstein (sGW) optimal transport, an extension of Gromov-Wasserstein by incorporating potential infinity patterns in the cost tensor. sGW enables the enforcement of application-induced constraints such as the preservation of pairwise distances by implementing the constraints as an infinity pattern. A numerical solver is proposed for the sGW problem and the effectiveness is demonstrated in various numerical experiments. The high-order constraints in sGW are transferred to constraints on the coupling matrix by solving a minimal vertex cover problem. The transformed problem is solved by the Mirror-C descent iteration coupled with the supervised optimal transport solver. In the numerical experiments, we first validate the proposed framework by applying it to matching synthetic datasets and investigating the impact of the model parameters. Additionally, we successfully apply sGW to real single-cell RNA sequencing data. Through comparisons with other Gromov-Wasserstein variants on real data, we demonstrate that sGW offers the novel utility of controlling distance preservation, leading to the automatic estimation of overlapping portions of datasets, which brings improved stability and flexibility in data-driven applications. △ Less

Submitted 11 January, 2024; originally announced January 2024.

MSC Class: 49Q22

Inverse source problem of the biharmonic equation from multi-frequency phaseless data

Authors: Yan Chang, Yukun Guo, Yue Zhao

Abstract: This work deals with an inverse source problem for the biharmonic wave equation. A two-stage numerical method is proposed to identify the unknown source from the multi-frequency phaseless data. In the first stage, we introduce some artificially auxiliary point sources to the inverse source system and establish a phase retrieval formula. Theoretically, we point out that the phase can be uniquely de… ▽ More This work deals with an inverse source problem for the biharmonic wave equation. A two-stage numerical method is proposed to identify the unknown source from the multi-frequency phaseless data. In the first stage, we introduce some artificially auxiliary point sources to the inverse source system and establish a phase retrieval formula. Theoretically, we point out that the phase can be uniquely determined and estimate the stability of this phase retrieval approach. Once the phase information is retrieved, the Fourier method is adopted to reconstruct the source function from the phased multi-frequency data. The proposed method is easy-to-implement and there is no forward solver involved in the reconstruction. Numerical experiments are conducted to verify the performance of the proposed method. △ Less

Submitted 5 January, 2024; originally announced January 2024.

Comments: 22 pages, 3 figures, 5 tables

MSC Class: 65N21; 35R30; 74J25

Closest targets in Russell graph measure of strongly monotonic efficiency for an extended facet production possibility set

Authors: Kazuyuki Sekitani, Yu Zhao

Abstract: This study modifies the Russell graph measure of efficiency to find the closest target from the assessed decision-making unit to the efficient frontier of an empirical production possibility set using an extended facet approach. The closest target requires a single improvement in either an output term or an input term, and needs some positive levels of all inputs to achieve the output target. The… ▽ More This study modifies the Russell graph measure of efficiency to find the closest target from the assessed decision-making unit to the efficient frontier of an empirical production possibility set using an extended facet approach. The closest target requires a single improvement in either an output term or an input term, and needs some positive levels of all inputs to achieve the output target. The practical advantage of the proposed efficiency measure is demonstrated numerically in comparison to other existing non-radial efficiency measures. △ Less

Submitted 24 December, 2023; originally announced December 2023.

Comments: 27 pages, 2 figures

MSC Class: 90-10

Half-integral levels

Authors: Yifei Zhao

Abstract: We construct equivalences among four notions associated to a reductive group scheme $G$: factorization super central extensions of the loop group of $G$ by $\mathbb G_m$ subject to a condition on the commutator, factorization super line bundles on the affine Grassmannian of $G$, rigidified sections of a quotient of $2$-truncated K-theory over the Zariski classifying stack of $G$, and combinatorial… ▽ More We construct equivalences among four notions associated to a reductive group scheme $G$: factorization super central extensions of the loop group of $G$ by $\mathbb G_m$ subject to a condition on the commutator, factorization super line bundles on the affine Grassmannian of $G$, rigidified sections of a quotient of $2$-truncated K-theory over the Zariski classifying stack of $G$, and combinatorial data defined by Brylinski and Deligne in a conjectural extension of their classification theorem. △ Less

Submitted 18 December, 2023; originally announced December 2023.

Comments: 49 pages

MSC Class: 14L15; 19C09; 81T40

A category equivalence on the Lie algebra of polynomial vector fields

Authors: Genqiang Liu, Yufang Zhao

Abstract: For any positive integer $n$, let $W_n=\text{Der}(\mathbb{C}[t_1,\dots,t_n])$. The subspaces $\mathfrak{h}_n=\text{Span}\{t_1\frac{\partial}{\partial{t_1}},\dots,t_n\frac{\partial}{\partial{t_n}}\}$ and $Δ_n=\text{Span}\{\frac{\partial}{\partial{t_1}},\dots,\frac{\partial}{\partial{t_n}}\}$ are two abelian subalgebras of $W_n$. We show that a full subcategory $Ω_{\mathbf{1}}$ of the category of… ▽ More For any positive integer $n$, let $W_n=\text{Der}(\mathbb{C}[t_1,\dots,t_n])$. The subspaces $\mathfrak{h}_n=\text{Span}\{t_1\frac{\partial}{\partial{t_1}},\dots,t_n\frac{\partial}{\partial{t_n}}\}$ and $Δ_n=\text{Span}\{\frac{\partial}{\partial{t_1}},\dots,\frac{\partial}{\partial{t_n}}\}$ are two abelian subalgebras of $W_n$. We show that a full subcategory $Ω_{\mathbf{1}}$ of the category of $W_n$-modules $M$ which are locally finite over $Δ_n$ is equivalent to some full subcategory of weight $W_n$-modules $M$ which are cuspidal modules when restricted to the subalgebra $\mathfrak{sl}_{n+1}$ of $W_n$. △ Less

Submitted 17 December, 2023; originally announced December 2023.

Positivity is undecidable in tensor products of free algebras

Authors: Arthur Mehta, William Slofstra, Yuming Zhao

Abstract: It is well known that an element of the algebra of noncommutative *-polynomials is positive in all *-representations if and only if it is a sum of squares. This provides an effective way to determine if a given *-polynomial is positive, by searching through sums of squares decompositions. We show that no such procedure exists for the tensor product of two noncommutative *-polynomial algebras: dete… ▽ More It is well known that an element of the algebra of noncommutative *-polynomials is positive in all *-representations if and only if it is a sum of squares. This provides an effective way to determine if a given *-polynomial is positive, by searching through sums of squares decompositions. We show that no such procedure exists for the tensor product of two noncommutative *-polynomial algebras: determining whether a *-polynomial of such an algebra is positive is coRE-hard. We also show that it is coRE-hard to determine whether a noncommutative *-polynomial is trace-positive. Our results hold if noncommutative *-polynomial algebras are replaced by other sufficiently free algebras such as group algebras of free groups or free products of cyclic groups. △ Less

Submitted 9 December, 2023; originally announced December 2023.

Comments: 45 pages

Conjugate gradient methods without line search for multiobjective optimization

Authors: Wang Chen, Yong Zhao, Liping Tang, Xinmin Yang

Abstract: This paper addresses an unconstrained multiobjective optimization problem where two or more continuously differentiable functions have to be minimized. We delve into the conjugate gradient methods proposed by Lucambio Pérez and Prudente (SIAM J Optim, 28(3): 2690--2720, 2018) for this type of problem. Instead of the Wolfe-type line search procedure used in their work, we employ a fixed stepsize fo… ▽ More This paper addresses an unconstrained multiobjective optimization problem where two or more continuously differentiable functions have to be minimized. We delve into the conjugate gradient methods proposed by Lucambio Pérez and Prudente (SIAM J Optim, 28(3): 2690--2720, 2018) for this type of problem. Instead of the Wolfe-type line search procedure used in their work, we employ a fixed stepsize formula (or no-line-search scheme), which can mitigate the selection pressure caused by multiple inequalities and avoid the computational cost associated with objective function evaluations in specific applications. The no-line-search scheme is utilized to derive the condition of Zoutendijk's type. Global convergence encompasses the vector extensions of Fletcher-Reeves, conjugate descent, Dai-Yuan, Polak-Ribière-Polyak and Hestenes-Stiefel parameters, subject to certain mild assumptions. Additionally, numerical experiments are conducted to demonstrate the practical performance of the proposed stepsize rule, and comparative analyses are made with the multiobjective steepest descent methods using the Armijo line search and the multiobjective conjugate gradient methods using the Wolfe-type line search. △ Less

Submitted 20 June, 2024; v1 submitted 4 December, 2023; originally announced December 2023.

Learning Free Terminal Time Optimal Closed-loop Control of Manipulators

Authors: Wei Hu, Yue Zhao, Weinan E, Jiequn Han, Jihao Long

Abstract: This paper presents a novel approach to learning free terminal time closed-loop control for robotic manipulation tasks, enabling dynamic adjustment of task duration and control inputs to enhance performance. We extend the supervised learning approach, namely solving selected optimal open-loop problems and utilizing them as training data for a policy network, to the free terminal time scenario. Thr… ▽ More This paper presents a novel approach to learning free terminal time closed-loop control for robotic manipulation tasks, enabling dynamic adjustment of task duration and control inputs to enhance performance. We extend the supervised learning approach, namely solving selected optimal open-loop problems and utilizing them as training data for a policy network, to the free terminal time scenario. Three main challenges are addressed in this extension. First, we introduce a marching scheme that enhances the solution quality and increases the success rate of the open-loop solver by gradually refining time discretization. Second, we extend the QRnet in Nakamura-Zimmerer et al. (2021b) to the free terminal time setting to address discontinuity and improve stability at the terminal state. Third, we present a more automated version of the initial value problem (IVP) enhanced sampling method from previous work (Zhang et al., 2022) to adaptively update the training dataset, significantly improving its quality. By integrating these techniques, we develop a closed-loop policy that operates effectively over a broad domain with varying optimal time durations, achieving near globally optimal total costs. △ Less

Submitted 29 November, 2023; originally announced November 2023.

Uniacute Spherical Codes

Authors: Saba Lepsveridze, Aleksandre Saatashvili, Yufei Zhao

Abstract: A spherical $L$-code, where $L \subseteq [-1,\infty)$, consists of unit vectors in $\mathbb{R}^d$ whose pairwise inner products are contained in $L$. Determining the maximum cardinality $N_L(d)$ of an $L$-code in $\mathbb{R}^d$ is a fundamental question in discrete geometry and has been extensively investigated for various choices of $L$. Our understanding in high dimensions is generally quite poo… ▽ More A spherical $L$-code, where $L \subseteq [-1,\infty)$, consists of unit vectors in $\mathbb{R}^d$ whose pairwise inner products are contained in $L$. Determining the maximum cardinality $N_L(d)$ of an $L$-code in $\mathbb{R}^d$ is a fundamental question in discrete geometry and has been extensively investigated for various choices of $L$. Our understanding in high dimensions is generally quite poor. Equiangular lines, corresponding to $L = \{-α, α\}$, is a rare and notable solved case. Bukh studied an extension of equiangular lines and showed that $N_L(d) = O_L(d)$ for $L = [-1, -β] \cup \{α\}$ with $α,β> 0$ (we call such $L$-codes "uniacute"), leaving open the question of determining the leading constant factor. Balla, Dräxler, Keevash, and Sudakov proved a "uniform bound" showing $\limsup_{d\to\infty} N_L(d)/d \le 2p$ for $L = [-1, -β] \cup \{α\}$ and $p = \lfloor α/β\rfloor + 1$. For which $(α,β)$ is this uniform bound tight? We completely answer this question. We develop a framework for studying uniacute codes, including a global structure theorem showing that the Gram matrix has an approximate $p$-block structure. We also formulate a notion of "modular codes," which we conjecture to be optimal in high dimensions. △ Less

Submitted 29 November, 2023; v1 submitted 29 November, 2023; originally announced November 2023.

Comments: 26 pages, 8 figures

Asymptotically Compatible Schemes for Nonlocal Ohta Kawasaki Model

Authors: Wangbo Luo, Yanxiang Zhao

Abstract: We study the asymptotical compatibility of the Fourier spectral method in multidimensional space for the Nonlocal Ohta-Kawasaka (NOK) model, which is proposed in our previous work. By introducing the Fourier collocation discretization for the spatial variable, we show that the asymptotical compatibility holds in 2D and 3D over a periodic domain. For the temporal discretization, we adopt the second… ▽ More We study the asymptotical compatibility of the Fourier spectral method in multidimensional space for the Nonlocal Ohta-Kawasaka (NOK) model, which is proposed in our previous work. By introducing the Fourier collocation discretization for the spatial variable, we show that the asymptotical compatibility holds in 2D and 3D over a periodic domain. For the temporal discretization, we adopt the second-order backward differentiation formula (BDF) method. We prove that for certain nonlocal kernels, the proposed time discretization schemes inherit the energy dissipation law. In the numerical experiments, we verify the asymptotical compatibility, the second-order temporal convergence rate, and the energy stability of the proposed schemes. More importantly, we discover a novel square lattice pattern when certain nonlocal kernel are applied in the model. In addition, our numerical experiments confirm the existence of an upper bound for the optimal number of bubbles in 2D for some specific nonlocal kernels. Finally, we numerically explore the promotion/demotion effect induced by the nonlocal horizon, which is consistent with the theoretical studies presented in our earlier work. △ Less

Submitted 25 November, 2023; originally announced November 2023.

The point insertion technique and open $r$-spin theories II: intersection theories in genus-zero

Authors: Ran J. Tessler, Yizhen Zhao

Abstract: The papers [5, 3, 6, 19, 20] initiated the study of open $r$-spin and open FJRW intersection theories, and related them to integrable hierarchies and mirror symmetry. This paper uses a new technique, the point insertion technique, developed in the prequel [36], to define new open r-spin and open FJRW intersection theories. These new constructions provide potential candidates for theories whose exi… ▽ More The papers [5, 3, 6, 19, 20] initiated the study of open $r$-spin and open FJRW intersection theories, and related them to integrable hierarchies and mirror symmetry. This paper uses a new technique, the point insertion technique, developed in the prequel [36], to define new open r-spin and open FJRW intersection theories. These new constructions provide potential candidates for theories whose existence was conjectured before: $\bullet$ K. Hori [23] predicted the existence of open $r$-spin theory with $\lfloor\frac{r}{2}\rfloor$ types of boundary states. The one constructed in [5, 3] has only one type of boundary state. In this work we describe $\lfloor\frac{r}{2}\rfloor$ open $r$-spin theories, labelled by $\mathfrak{h}\in\{0,\ldots,\lfloor\frac{r}{2}\rfloor-1\},$ where the $\mathfrak{h}$-th one has $\mathfrak{h}+1$ boundary states. We prove that the $\mathfrak{h}=0$ theory is equivalent to the [5, 3] construction, and calculate all intersection numbers for all these theories. $\bullet$ In [1] K. Aleshkin and C.C.M. Liu conjectured the existence of a quintic Fermat FJRW theory. We construct such an FJRW theory, and provide evidence that this is the conjectured theory. We also explain how the point insertion technique can be used for constructing other open enumerative theories, satisfying the same universal recursions. △ Less

Submitted 20 November, 2023; originally announced November 2023.

Devaney chaos of multiple mappings from a set-valued view

Authors: Yingcui Zhao

Abstract: In this paper, we introduce the definitions of periodic point, transitivity, sensitivity and Devaney chaos of multiple mappings from a set-valued perspective. We study the relation between multiple mappings and its continuous self-maps and show that multiple mappings and its continuous self-maps do not imply each other in terms of periodic point or transitivity. While the sufficient condition for… ▽ More In this paper, we introduce the definitions of periodic point, transitivity, sensitivity and Devaney chaos of multiple mappings from a set-valued perspective. We study the relation between multiple mappings and its continuous self-maps and show that multiple mappings and its continuous self-maps do not imply each other in terms of periodic point or transitivity. While the sufficient condition for multiple mappings to have periodic points, to be transitive and sensitive is provided separately. And we show that transitivity plus dense periodic point set implies sensitivity. Also, a sufficient condition for multiple mappings to be Devaney chaotic is given. △ Less

Submitted 6 November, 2023; originally announced November 2023.

Comments: arXiv admin note: substantial text overlap with arXiv:2311.02932

Shadowing, average shadowing and transitive properties of multiple mappings

Authors: Yingcui Zhao

Abstract: In 2016, Hou and Wang introduced the concept of multiple mappings based on iterated function system, which is an important branch of fractal theory. In this paper, we introduce the definitions of shadowing, average shadowing, transitive, weakly mixing, mixing, chain transitive and chain mixing properties of multiple mappings from a set-valued perspective. We show that both two continuous self-maps… ▽ More In 2016, Hou and Wang introduced the concept of multiple mappings based on iterated function system, which is an important branch of fractal theory. In this paper, we introduce the definitions of shadowing, average shadowing, transitive, weakly mixing, mixing, chain transitive and chain mixing properties of multiple mappings from a set-valued perspective. We show that both two continuous self-maps have shadowing (respectively, average shadowing, transitive, weakly mixing, mixing) property may not imply the multiple mappings they form has the corresponding property. While a sufficient condition for multiple mappings to be transitive (respective, weakly mixing, mixing, chain transitive, chain mixing) is given. Both shadowing and average shadowing properties are invariant under the iterative action of multiple mappings. Also, we study that for multiple mappings shadowing property plus chain mixing implies mixing and average shadowing properties implies chain transitivity. △ Less

Submitted 6 November, 2023; originally announced November 2023.

MSC Class: 54F17; 54E45