Showing 1–50 of 83 results for author: Cheng, M

Averaging principle for SDEs with singular drifts driven by $α$-stable processes

Authors: Mengyu Cheng, Zimo Hao, Xicheng Zhang

Abstract: In this paper, we investigate the convergence rate of the averaging principle for stochastic differential equations (SDEs) with $β$-Hölder drift driven by $α$-stable processes. More specifically, we first derive the Schauder estimate for nonlocal partial differential equations (PDEs) associated with the aforementioned SDEs, within the framework of Besov-Hölder spaces. Then we consider the case whe… ▽ More In this paper, we investigate the convergence rate of the averaging principle for stochastic differential equations (SDEs) with $β$-Hölder drift driven by $α$-stable processes. More specifically, we first derive the Schauder estimate for nonlocal partial differential equations (PDEs) associated with the aforementioned SDEs, within the framework of Besov-Hölder spaces. Then we consider the case where $(α,β)\in(0,2)\times(1-\tfracα{2},1)$. Using the Schauder estimate, we establish the strong convergence rate for the averaging principle. In particular, under suitable conditions we obtain the optimal rate of strong convergence when $(α,β)\in(\tfrac{2}{3},1]\times(2-\tfrac{3α}{2},1)\cup(1,2)\times(\tfracα{2},1)$. Furthermore, when $(α,β)\in(0,1]\times(1-α,1-\tfracα{2}]\cup(1,2)\times(\tfrac{1-α}{2},1-\tfracα{2}]$, we show the convergence of the martingale solutions of original systems to that of the averaged equation. When $α\in(1,2)$, the drift can be a distribution. △ Less

Submitted 19 September, 2024; originally announced September 2024.

Comments: 30 pages

MSC Class: 60H10; 34C29

Slow entropy and variational dynamical systems

Authors: Minhua Cheng, Carlos Ospina, Kurt Vinhage, Yibo Zhai

Abstract: We define variational properties for dynamical systems with subexponential complexity, and study these properties in certain specific examples. By computing the value of slow entropy directly, we show that some subshifts are not variational, while a class of interval exchange transformations are variational. We define variational properties for dynamical systems with subexponential complexity, and study these properties in certain specific examples. By computing the value of slow entropy directly, we show that some subshifts are not variational, while a class of interval exchange transformations are variational. △ Less

Submitted 19 October, 2024; v1 submitted 17 May, 2024; originally announced May 2024.

Recovery of transversely-isotropic elastic material parameters in induction motor rotors

Authors: Hanz Martin Cheng, Tapio Helin, Ville-Petteri Manninen, Timo Holopainen, Juha Jokinen, Samu Sorvari, Andreas Rupp

Abstract: We propose numerical algorithms for recovering parameters in eigenvalue problems for linear elasticity of transversely isotropic materials. Specifically, the algorithms are used to recover the elastic constants of a rotor core. Numerical tests show that in the noiseless setup, two pairs of bending modes are sufficient for recovering one to four parameters accurately. To recover all five parameters… ▽ More We propose numerical algorithms for recovering parameters in eigenvalue problems for linear elasticity of transversely isotropic materials. Specifically, the algorithms are used to recover the elastic constants of a rotor core. Numerical tests show that in the noiseless setup, two pairs of bending modes are sufficient for recovering one to four parameters accurately. To recover all five parameters that govern the elastic properties of electric engines accurately, we require three pairs of bending modes and one torsional mode. Moreover, we study the stability of the inversion method against multiplicative noise; for tests in which the data contained multiplicative noise of at most $1\%$, we find that all parameters can be recovered with an error less than $10\%$. △ Less

Submitted 10 May, 2024; originally announced May 2024.

MSC Class: 65Z05; 65C20

Fast OMP for Exact Recovery and Sparse Approximation

Authors: Huiyuan Yu, Jia He, Maggie Cheng

Abstract: Orthogonal Matching Pursuit (OMP) has been a powerful method in sparse signal recovery and approximation. However OMP suffers computational issue when the signal has large number of non-zeros. This paper advances OMP in two fronts: it offers a fast algorithm for the orthogonal projection of the input signal at each iteration, and a new selection criterion for making the greedy choice, which reduce… ▽ More Orthogonal Matching Pursuit (OMP) has been a powerful method in sparse signal recovery and approximation. However OMP suffers computational issue when the signal has large number of non-zeros. This paper advances OMP in two fronts: it offers a fast algorithm for the orthogonal projection of the input signal at each iteration, and a new selection criterion for making the greedy choice, which reduces the number of iterations it takes to recover the signal. The proposed modifications to OMP directly reduce the computational complexity. Experiment results show significant improvement over the classical OMP in computation time. The paper also provided a sufficient condition for exact recovery under the new greedy choice criterion. For general signals that may not have sparse representations, the paper provides a bound for the approximation error. The approximation error is at the same order as OMP but is obtained within fewer iterations and less time. △ Less

Submitted 29 March, 2024; originally announced April 2024.

3d Modularity Revisited

Authors: Miranda C. N. Cheng, Ioana Coman, Piotr Kucharski, Davide Passaro, Gabriele Sgroi

Abstract: The three-manifold topological invariants $\hat Z$ capture the half-index of the three-dimensional theory with ${\cal N}=2$ supersymmetry obtained by compactifying the M5 brane theory on the closed three-manifold. In 2019, surprising general relations between the $\hat Z$-invariants, quantum modular forms, and vertex algebras, have been proposed. In the meanwhile, an extensive array of examples ha… ▽ More The three-manifold topological invariants $\hat Z$ capture the half-index of the three-dimensional theory with ${\cal N}=2$ supersymmetry obtained by compactifying the M5 brane theory on the closed three-manifold. In 2019, surprising general relations between the $\hat Z$-invariants, quantum modular forms, and vertex algebras, have been proposed. In the meanwhile, an extensive array of examples have been studied, but several general important structural questions remain. First, for many three-manifolds it was observed that the different $\hat Z$-invariants for the same three-manifolds are quantum modular forms that span a subspace of a Weil representation for the modular group $SL_2(Z)$, corresponding to the structure of vector-valued quantum modular forms. We elucidate the meaning of this vector-valued quantum modular form structure by first proposing the analogue $\hat Z$-invariants with supersymmetric defects, and subsequently showing that the full vector-valued quantum modular form is precisely the object capturing all the $\hat Z$-invariants, with and without defects, of a given three-manifold. Second, it was expected that matching radial limits is a key feature of $\hat Z$-invariants when changing the orientation of the plumbed three-manifold, suggesting the relevance of mock modularity. We substantiate the conjecture by providing explicit proposals for such $\hat Z$-invariants for an infinite family of three-manifolds and verify their mock modularity and limits. Third, we initiate the study of the vertex algebra structure of the mock type invariants by showcasing a systematic way to construct cone vertex operator algebras associated to these invariants, which can be viewed as the partner of logarithmic vertex operator algebras in this context. △ Less

Submitted 25 March, 2024; v1 submitted 21 March, 2024; originally announced March 2024.

Comments: 59 pages, typos corrected

Learning Domain-Independent Green's Function For Elliptic Partial Differential Equations

Authors: Pawan Negi, Maggie Cheng, Mahesh Krishnamurthy, Wenjun Ying, Shuwang Li

Abstract: Green's function characterizes a partial differential equation (PDE) and maps its solution in the entire domain as integrals. Finding the analytical form of Green's function is a non-trivial exercise, especially for a PDE defined on a complex domain or a PDE with variable coefficients. In this paper, we propose a novel boundary integral network to learn the domain-independent Green's function, ref… ▽ More Green's function characterizes a partial differential equation (PDE) and maps its solution in the entire domain as integrals. Finding the analytical form of Green's function is a non-trivial exercise, especially for a PDE defined on a complex domain or a PDE with variable coefficients. In this paper, we propose a novel boundary integral network to learn the domain-independent Green's function, referred to as BIN-G. We evaluate the Green's function in the BIN-G using a radial basis function (RBF) kernel-based neural network. We train the BIN-G by minimizing the residual of the PDE and the mean squared errors of the solutions to the boundary integral equations for prescribed test functions. By leveraging the symmetry of the Green's function and controlling refinements of the RBF kernel near the singularity of the Green function, we demonstrate that our numerical scheme enables fast training and accurate evaluation of the Green's function for PDEs with variable coefficients. The learned Green's function is independent of the domain geometries, forcing terms, and boundary conditions in the boundary integral formulation. Numerical experiments verify the desired properties of the method and the expected accuracy for the two-dimensional Poisson and Helmholtz equations with variable coefficients. △ Less

Submitted 30 January, 2024; originally announced January 2024.

The Braids on your Blanket

Authors: Michelle Cheng, Robert Laugwitz

Abstract: In this expositional essay, we introduce some elements of the study of groups by analysing the braid pattern on a knitted blanket. We determine that the blanket features pure braids with a minimal number of crossings. Moreover, we determine polynomial invariants associated to the links obtained by closing the braid patterns of the blanket. In this expositional essay, we introduce some elements of the study of groups by analysing the braid pattern on a knitted blanket. We determine that the blanket features pure braids with a minimal number of crossings. Moreover, we determine polynomial invariants associated to the links obtained by closing the braid patterns of the blanket. △ Less

Submitted 18 October, 2023; originally announced October 2023.

Comments: Expositional article for a general readership. 32 pages, several figures

MSC Class: 00A66 (Primary) 00-01; 20F36; 57K10 (Secondary)

Journal ref: Journal of Humanistic Mathematics, Volume 14 Issue 2 (July 2024), pages 286-337. Available at: https://scholarship.claremont.edu/jhm/vol14/iss2/10

Averaging principle and normal deviation for multi-scale SDEs with polynomial nonlinearity

Authors: Mengyu Cheng, Zhenxin Liu, Michael Röckner

Abstract: We investigate three types of averaging principles and the normal deviation for multi-scale stochastic differential equations (in short, SDEs) with polynomial nonlinearity. More specifically, we first demonstrate the strong convergence of the solution of SDEs, which involves highly oscillating components and fast processes, to that of the averaged equation. Then we investigate the small fluctuatio… ▽ More We investigate three types of averaging principles and the normal deviation for multi-scale stochastic differential equations (in short, SDEs) with polynomial nonlinearity. More specifically, we first demonstrate the strong convergence of the solution of SDEs, which involves highly oscillating components and fast processes, to that of the averaged equation. Then we investigate the small fluctuations of the system around its average, and show that the normalized difference weakly converges to an Ornstein-Uhlenbeck type process, which can be viewed as a functional central limit theorem. Additionally, we show that the attractor of the original system tends to that of the averaged equation in probability measure space as the time scale $\varepsilon$ goes to zero. Finally, we establish the second Bogolyubov theorem; that is to say, we prove that there exists a quasi-periodic solution in a neighborhood of the stationary solution of the averaged equation when the $\varepsilon$ is small. △ Less

Submitted 21 August, 2023; originally announced August 2023.

Comments: 39 pages

Natural Actor-Critic for Robust Reinforcement Learning with Function Approximation

Authors: Ruida Zhou, Tao Liu, Min Cheng, Dileep Kalathil, P. R. Kumar, Chao Tian

Abstract: We study robust reinforcement learning (RL) with the goal of determining a well-performing policy that is robust against model mismatch between the training simulator and the testing environment. Previous policy-based robust RL algorithms mainly focus on the tabular setting under uncertainty sets that facilitate robust policy evaluation, but are no longer tractable when the number of states scales… ▽ More We study robust reinforcement learning (RL) with the goal of determining a well-performing policy that is robust against model mismatch between the training simulator and the testing environment. Previous policy-based robust RL algorithms mainly focus on the tabular setting under uncertainty sets that facilitate robust policy evaluation, but are no longer tractable when the number of states scales up. To this end, we propose two novel uncertainty set formulations, one based on double sampling and the other on an integral probability metric. Both make large-scale robust RL tractable even when one only has access to a simulator. We propose a robust natural actor-critic (RNAC) approach that incorporates the new uncertainty sets and employs function approximation. We provide finite-time convergence guarantees for the proposed RNAC algorithm to the optimal robust policy within the function approximation error. Finally, we demonstrate the robust performance of the policy learned by our proposed RNAC approach in multiple MuJoCo environments and a real-world TurtleBot navigation task. △ Less

Submitted 10 December, 2023; v1 submitted 17 July, 2023; originally announced July 2023.

Comments: Neurips 2023

Matroidal Entropy Functions: Constructions, Characterizations and Representations

Authors: Qi Chen, Minquan Cheng, Baoming Bai

Abstract: Matroidal entropy functions are entropy functions in the form $\mathbf{h} = \log v \cdot \mathbf{r}_M$ , where $v \ge 2$ is an integer and $\mathbf{r}_M$ is the rank function of a matroid $M$. They can be applied into capacity characterization and code construction of information theory problems such as network coding, secret sharing, index coding and locally repairable code. In this paper, by con… ▽ More Matroidal entropy functions are entropy functions in the form $\mathbf{h} = \log v \cdot \mathbf{r}_M$ , where $v \ge 2$ is an integer and $\mathbf{r}_M$ is the rank function of a matroid $M$. They can be applied into capacity characterization and code construction of information theory problems such as network coding, secret sharing, index coding and locally repairable code. In this paper, by constructing the variable strength arrays of some matroid operations, we characterized matroidal entropy functions induced by regular matroids and some matroids with the same p-characteristic set as uniform matroid $U_{2,4}$. △ Less

Submitted 30 January, 2024; v1 submitted 29 June, 2023; originally announced June 2023.

Comments: 13 pages, 3 figures, accepted by IEEE Transactions on Information Theory

Class Numbers, Congruent Numbers and Umbral Moonshine

Authors: Miranda C. N. Cheng, John F. R. Duncan, Michael H. Mertens

Abstract: In earlier work we initiated a program to study relationships between finite groups and arithmetic geometric invariants of modular curves in a systematic way. In the present work we continue this program, with a focus on the two smallest sporadic simple Mathieu groups. To do this we first elucidate a connection between a special case of umbral moonshine and the imaginary quadratic class numbers. T… ▽ More In earlier work we initiated a program to study relationships between finite groups and arithmetic geometric invariants of modular curves in a systematic way. In the present work we continue this program, with a focus on the two smallest sporadic simple Mathieu groups. To do this we first elucidate a connection between a special case of umbral moonshine and the imaginary quadratic class numbers. Then we use this connection to classify a distinguished set of modules for the smallest sporadic Mathieu group. Finally we establish consequences of the classification for the congruent number problem of antiquity. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 36 pages, 5 tables

MSC Class: 11F22; 11F37; 11G05; 11G40; 20C34

The Average Number of Goldbach Representations and Zero-Free Regions of the Riemann Zeta-Function

Authors: Keith Billington, Maddie Cheng, Jordan Schettler, Ade Irma Suriajaya

Abstract: In this paper, we prove an unconditional form of Fujii's formula for the average number of Goldbach representations and show that the error in this formula is determined by a general zero-free region of the Riemann zeta-function, and vice versa. In particular, we describe the error in the unconditional formula in terms of the remainder in the Prime Number Theorem which connects the error to zero-f… ▽ More In this paper, we prove an unconditional form of Fujii's formula for the average number of Goldbach representations and show that the error in this formula is determined by a general zero-free region of the Riemann zeta-function, and vice versa. In particular, we describe the error in the unconditional formula in terms of the remainder in the Prime Number Theorem which connects the error to zero-free regions of the Riemann zeta-function. △ Less

Submitted 11 September, 2023; v1 submitted 15 June, 2023; originally announced June 2023.

Comments: 22 pages (content in 20 pages), a student project conducted at SJSU under Kyushu University SENTAN-Q

Report number: RIKEN-iTHEMS-Report-23 MSC Class: 11P32; 11N05; 11N37; 11M26

Expanding solutions near unstable Lane-Emden stars

Authors: Ming Cheng, Xing Cheng, Zhiwu Lin

Abstract: We consider the compressible Euler-Poisson equations for polytropes $P(ρ)=Kρ^γ$ with $γ\in \left(\frac{6}{5},\frac{4}{3} \right]$ and the white dwarf stars. For $γ=\frac{4}{3},$ we show that there exists a global weak solution for the spherically symmetric initial data with mass less than the mass of the Lane-Emden stars (i.e. non-rotating polytropes). For… ▽ More We consider the compressible Euler-Poisson equations for polytropes $P(ρ)=Kρ^γ$ with $γ\in \left(\frac{6}{5},\frac{4}{3} \right]$ and the white dwarf stars. For $γ=\frac{4}{3},$ we show that there exists a global weak solution for the spherically symmetric initial data with mass less than the mass of the Lane-Emden stars (i.e. non-rotating polytropes). For $γ\in \left(\frac{6}{5},\frac{4}{3} \right)$, we show the existence of global weak solution for spherical symmetric initial data in an invariant set containing a neighborhood of Lane-Emden stars. Moreover, the support of these solution expands to infinity. As a corollary, this proves the strong instability of the Lane-Emden stars for $γ\in \left( \frac{6}{5},\frac{4}{3}\right] $. For $γ\in \left(\frac{6}{5},\frac{4}{3} \right),$ our results provide the first example of expanding solutions. For white dwarf stars, we prove that the solution cannot collapse if the mass of initial data is less than the Chandrasekhar limit mass, which is the supremum of the mass of the non-rotating white dwarf stars. Our proof strongly uses the variational characterization of the Lane-Emden stars. First, we relate the best constant of a Hardy-Littlewood type inequality with the mass of the Lane-Emden stars with $γ=\frac{4}{3}$, which is further shown to equal the Chandrasekhar limit mass. For $γ\in\left( \frac{6}{5},\frac{4}{3}\right) $, we show that the Lane-Emden stars are minimizers of an energy-mass functional subject to a Pohozaev type constraint. This is crucial in the construction of the invariant set of expanding solutions. △ Less

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

MSC Class: 35Q85; 35R35; 35A15; 85A05

Quantum Modular $\widehat Z{}^G$-Invariants

Authors: Miranda C. N. Cheng, Ioana Coman, Davide Passaro, Gabriele Sgroi

Abstract: We study the quantum modular properties of $\widehat Z{}^G$-invariants of closed three-manifolds. Higher depth quantum modular forms are expected to play a central role for general three-manifolds and gauge groups $G$. In particular, we conjecture that for plumbed three-manifolds whose plumbing graphs have $n$ junction nodes with definite signature and for rank $r$ gauge group $G$, that… ▽ More We study the quantum modular properties of $\widehat Z{}^G$-invariants of closed three-manifolds. Higher depth quantum modular forms are expected to play a central role for general three-manifolds and gauge groups $G$. In particular, we conjecture that for plumbed three-manifolds whose plumbing graphs have $n$ junction nodes with definite signature and for rank $r$ gauge group $G$, that $\widehat Z{}^G$ is related to a quantum modular form of depth $nr$. We prove this for $G={\rm SU}(3)$ and for an infinite class of three-manifolds (weakly negative Seifert with three exceptional fibers). We also investigate the relation between the quantum modularity of $\widehat Z{}^G$-invariants of the same three-manifold with different gauge group $G$. We conjecture a recursive relation among the iterated Eichler integrals relevant for $\widehat Z{}^G$ with $G={\rm SU}(2)$ and ${\rm SU}(3)$, for negative Seifert manifolds with three exceptional fibers. This is reminiscent of the recursive structure among mock modular forms playing the role of Vafa-Witten invariants for ${\rm SU}(N)$. We prove the conjecture when the three-manifold is moreover an integral homological sphere. △ Less

Submitted 9 March, 2024; v1 submitted 8 April, 2023; originally announced April 2023.

Journal ref: SIGMA 20 (2024), 018, 52 pages

Turning point principle for stability of viscous gaseous stars

Authors: Ming Cheng, Zhiwu Lin, Yucong Wang

Abstract: We consider stability of non-rotating viscous gaseous stars modeled by the Navier-Stokes-Poisson system. Under general assumptions on the equations of states, we proved that the number of unstable modes for the linearized Navier-Stokes-Poisson system equals that of the linearized Euler-Poisson system modeling inviscid gaseous stars. In particular, the turning point principle holds true for non-rot… ▽ More We consider stability of non-rotating viscous gaseous stars modeled by the Navier-Stokes-Poisson system. Under general assumptions on the equations of states, we proved that the number of unstable modes for the linearized Navier-Stokes-Poisson system equals that of the linearized Euler-Poisson system modeling inviscid gaseous stars. In particular, the turning point principle holds true for non-rotating stars with or without viscosity. That is, the transition of stability only occurs at the extrema of the total mass and the number of unstable modes is determined by the mass-radius curve. For the proof, we establish an infinite dimensional Kelvin-Tait-Chetaev theorem for a class of linear second order PDEs with dissipation. Moreover, we prove that linear stability implies nonlinear asymptotic stability and linear instability implies nonlinear instability for Navier-Stokes-Poisson system. △ Less

Submitted 18 January, 2023; originally announced January 2023.

A Light-speed Linear Program Solver for Personalized Recommendation with Diversity Constraints

Authors: Haoyue Wang, Miao Cheng, Kinjal Basu, Aman Gupta, Keerthi Selvaraj, Rahul Mazumder

Abstract: We study a structured linear program (LP) that emerges in the need of ranking candidates or items in personalized recommender systems. Since the candidate set is only known in real time, the LP also needs to be formed and solved in real time. Latency and user experience are major considerations, requiring the LP to be solved within just a few milliseconds. Although typical instances of the problem… ▽ More We study a structured linear program (LP) that emerges in the need of ranking candidates or items in personalized recommender systems. Since the candidate set is only known in real time, the LP also needs to be formed and solved in real time. Latency and user experience are major considerations, requiring the LP to be solved within just a few milliseconds. Although typical instances of the problem are not very large in size, this stringent time limit appears to be beyond the capability of most existing (commercial) LP solvers, which can take $20$ milliseconds or more to find a solution. Thus, reliable methods that address the real-world complication of latency become necessary. In this paper, we propose a fast specialized LP solver for a structured problem with diversity constraints. Our method solves the dual problem, making use of the piece-wise affine structure of the dual objective function, with an additional screening technique that helps reduce the dimensionality of the problem as the algorithm progresses. Experiments reveal that our method can solve the problem within roughly 1 millisecond, yielding a 20x improvement in speed over efficient off-the-shelf LP solvers. This speed-up can help improve the quality of recommendations without affecting user experience, highlighting how optimization can provide solid orthogonal value to machine-learned recommender systems. △ Less

Submitted 22 November, 2022; originally announced November 2022.

Multi-access Coded Caching with Optimal Rate and Linear Subpacketization under PDA and Consecutive Cyclic Placement

Authors: Jinyu Wang, Minquan Cheng, Youlong Wu

Abstract: This work considers the multi-access caching system proposed by Hachem et al., where each user has access to L neighboring caches in a cyclic wrap-around fashion. We first propose a placement strategy called the consecutive cyclic placement, which achieves the maximal local caching gain. Then under the consecutive cyclic placement, we derive the optimal coded caching gain from the perspective of P… ▽ More This work considers the multi-access caching system proposed by Hachem et al., where each user has access to L neighboring caches in a cyclic wrap-around fashion. We first propose a placement strategy called the consecutive cyclic placement, which achieves the maximal local caching gain. Then under the consecutive cyclic placement, we derive the optimal coded caching gain from the perspective of Placement Delivery Array (PDA), thus obtaining a lower bound on the rate of PDA. Finally, under the consecutive cyclic placement, we construct a class of PDA, leading to a multi-access coded caching scheme with linear subpacketization, which achieves our derived lower bound for some parameters; while for other parameters, the achieved coded caching gain is only 1 less than the optimal one. Analytical and numerical comparisons of the proposed scheme with existing schemes are provided to validate the performance. △ Less

Submitted 7 September, 2022; originally announced September 2022.

Comments: 30 pages, 7 figures

Strong and weak convergence for averaging principle of DDSDE with singular drift

Authors: Mengyu Cheng, Zimo Hao, Michael Röckner

Abstract: In this paper, we study the averaging principle for distribution dependent stochastic differential equations with drift in localized $L^p$ spaces. Using Zvonkin's transformation and estimates for solutions to Kolmogorov equations, we prove that the solutions of the original system strongly and weakly converge to the solution of the averaged system as the time scale $\eps$ goes to zero. Moreover, w… ▽ More In this paper, we study the averaging principle for distribution dependent stochastic differential equations with drift in localized $L^p$ spaces. Using Zvonkin's transformation and estimates for solutions to Kolmogorov equations, we prove that the solutions of the original system strongly and weakly converge to the solution of the averaged system as the time scale $\eps$ goes to zero. Moreover, we obtain rates of the strong and weak convergence that depend on $p$ respectively. △ Less

Submitted 26 October, 2022; v1 submitted 25 July, 2022; originally announced July 2022.

Comments: 41 pages

Gauging Lie group symmetry in (2+1)d topological phases

Authors: Meng Cheng, Po-Shen Hsin, Chao-Ming Jian

Abstract: We present a general algebraic framework for gauging a 0-form compact, connected Lie group symmetry in (2+1)d topological phases. Starting from a symmetry fractionalization pattern of the Lie group $G$, we first extend $G$ to a larger symmetry group $\tilde{G}$, such that there is no fractionalization with respect to $\tilde{G}$ in the topological phase, and the effect of gauging $\tilde{G}$ is to… ▽ More We present a general algebraic framework for gauging a 0-form compact, connected Lie group symmetry in (2+1)d topological phases. Starting from a symmetry fractionalization pattern of the Lie group $G$, we first extend $G$ to a larger symmetry group $\tilde{G}$, such that there is no fractionalization with respect to $\tilde{G}$ in the topological phase, and the effect of gauging $\tilde{G}$ is to tensor the original theory with a $\tilde{G}$ Chern-Simons theory. To restore the desired gauge symmetry, one then has to gauge an appropriate one-form symmetry (or, condensing certain Abelian anyons) to obtain the final result. Studying the consistency of the gauging procedure leads to compatibility conditions between the symmetry fractionalization pattern and the Hall conductance. When the gauging can not be consistently done (i.e. the compatibility conditions can not be satisfied), the symmetry $G$ with the fractionalization pattern has an 't Hooft anomaly and we present a general method to determine the (3+1)d topological term for the anomaly. We provide many examples, including projective simple Lie groups and unitary groups to illustrate our approach. △ Less

Submitted 30 November, 2022; v1 submitted 30 May, 2022; originally announced May 2022.

Comments: 31+5 pages

Journal ref: SciPost Phys. 14, 100 (2023)

The generalized reciprocal distance matrix of graphs

Authors: Gui-Xian Tian, Mei-Jiao Cheng, Shu-Yu Cui

Abstract: Let $G$ be a simple undirected connected graph with the Harary matrix $RD(G)$, which is also called the reciprocal distance matrix of $G$. The reciprocal distance signless Laplacian matrix of $G$ is $RQ(G)=RT(G)+RD(G)$, where $RT(G)$ denotes the diagonal matrix of the vertex reciprocal transmissions of graph $G$. This paper intends to introduce a new matrix $RD_α(G)=αRT(G)+(1-α)RD(G)$,… ▽ More Let $G$ be a simple undirected connected graph with the Harary matrix $RD(G)$, which is also called the reciprocal distance matrix of $G$. The reciprocal distance signless Laplacian matrix of $G$ is $RQ(G)=RT(G)+RD(G)$, where $RT(G)$ denotes the diagonal matrix of the vertex reciprocal transmissions of graph $G$. This paper intends to introduce a new matrix $RD_α(G)=αRT(G)+(1-α)RD(G)$, $α\in [0,1]$, to track the gradual change from $RD(G)$ to $RQ(G)$. First, we describe completely the eigenvalues of $RD_α(G)$ of some special graphs. Then we obtain serval basic properties of $RD_α(G)$ including inequalities that involve the spectral radii of the reciprocal distance matrix, reciprocal distance signless Laplacian matrix and $RD_α$-matrix of $G$. We also provide some lower and upper bounds of the spectral radius of $RD_α$-matrix. Finally, we depict the extremal graphs with maximal spectral radius of the $RD_α$-matrix among all connected graphs of fixed order and precise vertex connectivity, edge connectivity, chromatic number and independence number, respectively. △ Less

Submitted 7 April, 2022; originally announced April 2022.

Comments: 21 pages

MSC Class: 05C50; 15A18

Complete flux scheme for variable velocity fields: coupling between the advection-diffusion equation and the Poisson equation for the velocity field

Authors: Hanz Martin Cheng, Jan ten Thije Boonkkamp

Abstract: In this work, we consider an advection-diffusion equation, coupled to a Poisson equation for the velocity field. This type of coupling is typically encountered in models arising from plasma physics or porous media flow. The aim of this work is to build upon the complete flux scheme (an improvement over the Scharfetter-Gummel scheme by considering the contribution of the source term), so that its s… ▽ More In this work, we consider an advection-diffusion equation, coupled to a Poisson equation for the velocity field. This type of coupling is typically encountered in models arising from plasma physics or porous media flow. The aim of this work is to build upon the complete flux scheme (an improvement over the Scharfetter-Gummel scheme by considering the contribution of the source term), so that its second-order convergence, which is uniform in Péclet numbers, carries over to these models. This is done by considering a piecewise linear approximation of the velocity field, which is then used for defining upwind-adjusted Péclet numbers. △ Less

Submitted 13 February, 2023; v1 submitted 7 April, 2022; originally announced April 2022.

MSC Class: 65N08

Class Numbers, Cyclic Simple Groups and Arithmetic

Authors: Miranda C. N. Cheng, John F. R. Duncan, Michael H. Mertens

Abstract: Here we initiate a program to study relationships between finite groups and arithmetic-geometric invariants in a systematic way. To do this we first introduce a notion of optimal module for a finite group in the setting of holomorphic mock Jacobi forms. Then we classify optimal modules for the cyclic groups of prime order, in the special case of weight two and index one, where class numbers of ima… ▽ More Here we initiate a program to study relationships between finite groups and arithmetic-geometric invariants in a systematic way. To do this we first introduce a notion of optimal module for a finite group in the setting of holomorphic mock Jacobi forms. Then we classify optimal modules for the cyclic groups of prime order, in the special case of weight two and index one, where class numbers of imaginary quadratic fields play an important role. Finally we exhibit a connection between the classification we establish and the arithmetic geometry of imaginary quadratic twists of modular curves of prime level. △ Less

Submitted 13 March, 2023; v1 submitted 2 April, 2022; originally announced April 2022.

Comments: 36 pages, v2 contains minor changes according to referee's comments

MSC Class: 11F22; 11F37; 11F50; 11G05; 11G40

Cone Vertex Algebras, Mock Theta Functions, and Umbral Moonshine Modules

Authors: Miranda C. N. Cheng, Gabriele Sgroi

Abstract: We describe a family of indefinite theta functions of signature $(1,1)$ that can be expressed in terms of trace functions of vertex algebras built from cones in lattices. The family of indefinite theta functions considered has interesting connections with mock theta functions and Appell-Lerch sums. We use these relations to write the McKay-Thompson series of umbral moonshine at lambency… ▽ More We describe a family of indefinite theta functions of signature $(1,1)$ that can be expressed in terms of trace functions of vertex algebras built from cones in lattices. The family of indefinite theta functions considered has interesting connections with mock theta functions and Appell-Lerch sums. We use these relations to write the McKay-Thompson series of umbral moonshine at lambency $\ell=8,12,16$ in terms of trace functions of vertex algebras modules, and thereby provide the modules for these instances of umbral moonshine. △ Less

Submitted 6 March, 2022; originally announced March 2022.

Averaging principle for stochastic complex Ginzburg-Landau equations

Authors: Mengyu Cheng, Zhenxin Liu, Michael Röckner

Abstract: Averaging principle is an effective method for investigating dynamical systems with highly oscillating components. In this paper, we study three types of averaging principle for stochastic complex Ginzburg-Landau equations. Firstly, we prove that the solution of the original equation converges to that of the averaged equation on finite intervals as the time scale $\varepsilon$ goes to zero when th… ▽ More Averaging principle is an effective method for investigating dynamical systems with highly oscillating components. In this paper, we study three types of averaging principle for stochastic complex Ginzburg-Landau equations. Firstly, we prove that the solution of the original equation converges to that of the averaged equation on finite intervals as the time scale $\varepsilon$ goes to zero when the initial data are the same. Secondly, we show that there exists a unique recurrent solution (in particular, periodic, almost periodic, almost automorphic, etc.) to the original equation in a neighborhood of the stationary solution of the averaged equation when the time scale is small. Finally, we establish the global averaging principle in weak sense, i.e. we show that the attractor of original system tends to that of the averaged equation in probability measure space as $\varepsilon$ goes to zero. △ Less

Submitted 21 November, 2022; v1 submitted 4 March, 2022; originally announced March 2022.

Comments: 35 pages

MSC Class: 35Q56; 60H15; 37B20

Coded Caching for Two-Dimensional Multi-Access Networks

Authors: Mingming Zhang, Kai Wan, Minquan Cheng, Giuseppe Caire

Abstract: This paper studies a novel multi-access coded caching (MACC) model in the two-dimensional (2D) topology, which is a generalization of the one-dimensional (1D) MACC model proposed by Hachem et al. The 2D MACC model is formed by a server containing $N$ files, $K_1\times K_2$ cache-nodes with $M$ files located at a grid with $K_1$ rows and $K_2$ columns, and $K_1\times K_2$ cache-less users where eac… ▽ More This paper studies a novel multi-access coded caching (MACC) model in the two-dimensional (2D) topology, which is a generalization of the one-dimensional (1D) MACC model proposed by Hachem et al. The 2D MACC model is formed by a server containing $N$ files, $K_1\times K_2$ cache-nodes with $M$ files located at a grid with $K_1$ rows and $K_2$ columns, and $K_1\times K_2$ cache-less users where each user is connected to $L^2$ nearby cache-nodes. The server is connected to the users through an error-free shared link, while the users can retrieve the cached content of the connected cache-nodes without cost. Our objective is to minimize the worst-case transmission load over all possible users' demands. In this paper, we first propose a grouping scheme for the case where $K_1$ and $K_2$ are divisible by $L$. By partitioning the cache-nodes and users into $L^2$ groups such that no two users in the same group share any cache-node, we use the shared-link coded caching scheme proposed by Maddah-Ali and Niesen for each group. Then for any model parameters satisfying $\min\{K_1,K_2\}>L$, we propose a transformation approach which constructs a 2D MACC scheme from two classes of 1D MACC schemes in vertical and horizontal projections, respectively. As a result, we can construct 2D MACC schemes that achieve maximum local caching gain and improved coded caching gain, compared to the baseline scheme by a direct extension from 1D MACC schemes. △ Less

Submitted 23 June, 2022; v1 submitted 27 January, 2022; originally announced January 2022.

Comments: 40 pages, 11 figures

Gauging U(1) symmetry in (2+1)d topological phases

Authors: Meng Cheng, Chao-Ming Jian

Abstract: We study the gauging of a global U(1) symmetry in a gapped system in (2+1)d. The gauging procedure has been well-understood for a finite global symmetry group, which leads to a new gapped phase with emergent gauge structure and can be described algebraically using the mathematical framework of modular tensor category (MTC). We develop a categorical description of U(1) gauging in an MTC, taking int… ▽ More We study the gauging of a global U(1) symmetry in a gapped system in (2+1)d. The gauging procedure has been well-understood for a finite global symmetry group, which leads to a new gapped phase with emergent gauge structure and can be described algebraically using the mathematical framework of modular tensor category (MTC). We develop a categorical description of U(1) gauging in an MTC, taking into account the dynamics of U(1) gauge field absent in the finite group case. When the ungauged system has a non-zero Hall conductance, the gauged theory remains gapped and we determine the complete set of anyon data for the gauged theory. On the other hand, when the Hall conductance vanishes, we argue that gauging has the same effect of condensing a special Abelian anyon nucleated by inserting $2π$ U(1) flux. We apply our procedure to the SU(2)$_k$ MTCs and derive the full MTC data for the $\mathbb{Z}_k$ parafermion MTCs. We also discuss a dual U(1) symmetry that emerges after the original U(1) symmetry of an MTC is gauged. △ Less

Submitted 31 May, 2022; v1 submitted 18 January, 2022; originally announced January 2022.

Comments: 11+10pages; 0+3 figures

Journal ref: SciPost Phys. 12, 202 (2022)

3-Manifolds and VOA Characters

Authors: Miranda C. N. Cheng, Sungbong Chun, Boris Feigin, Francesca Ferrari, Sergei Gukov, Sarah M. Harrison, Davide Passaro

Abstract: By studying the properties of $q$-series $\widehat Z$-invariants, we develop a dictionary between 3-manifolds and vertex algebras. In particular, we generalize previously known entries in this dictionary to Lie groups of higher rank, to 3-manifolds with toral boundaries, and to BPS partition functions with line operators. This provides a new physical realization of logarithmic vertex algebras in t… ▽ More By studying the properties of $q$-series $\widehat Z$-invariants, we develop a dictionary between 3-manifolds and vertex algebras. In particular, we generalize previously known entries in this dictionary to Lie groups of higher rank, to 3-manifolds with toral boundaries, and to BPS partition functions with line operators. This provides a new physical realization of logarithmic vertex algebras in the framework of the 3d-3d correspondence and opens new avenues for their future study. For example, we illustrate how invoking a knot-quiver correspondence for $\widehat{Z}$-invariants leads to many infinite families of new fermionic formulae for VOA characters. △ Less

Submitted 4 November, 2022; v1 submitted 12 January, 2022; originally announced January 2022.

Comments: 85 pages, 3 figures, 6 tables

Riemann-Hurwitz theorem and second main theorem for harmonic morphisms on graphs or metrized complexes

Authors: Tingbin Cao, Mengnan Cheng

Abstract: In this article, we mainly obtain the Riemann-Hurwitz theorems for harmonic morphisms on (vertex-weighted) metric graphs or metrized complexes of algebraic curves, inspired of the recent work on harmonic morphisms of graphs or metrized complexes due to many researchers. By making use of these Riemann-Hurwitz theorems, we then systematically establish the second main theorems for harmonic morphisms… ▽ More In this article, we mainly obtain the Riemann-Hurwitz theorems for harmonic morphisms on (vertex-weighted) metric graphs or metrized complexes of algebraic curves, inspired of the recent work on harmonic morphisms of graphs or metrized complexes due to many researchers. By making use of these Riemann-Hurwitz theorems, we then systematically establish the second main theorems for harmonic morphisms on finite graphs, vertex-weighted graphs, (vertex-weighted) metric graphs or metrized complexes of algebraic curves, from the viewpoint of Nevanlinna theory. △ Less

Submitted 12 January, 2022; originally announced January 2022.

Comments: 36 pages

MSC Class: 14N10; 05E14; 05C10; 30D35

Global existence for systems of nonlinear wave and Klein-Gordon equations in two space dimensions under a kind of the weak null condition

Authors: Minggang Cheng

Abstract: We consider the coupled systems of nonlinear wave and Klein-Gordon equations in two space dimensions with cubic nonlinearity. For this kind of systems, the small data global existence is already known if the cubic nonlinearity satisfies a certain condition related to the null condition. In this article, our aim is to investigate the small data global existence under a condition related to the weak… ▽ More We consider the coupled systems of nonlinear wave and Klein-Gordon equations in two space dimensions with cubic nonlinearity. For this kind of systems, the small data global existence is already known if the cubic nonlinearity satisfies a certain condition related to the null condition. In this article, our aim is to investigate the small data global existence under a condition related to the weak null condition. We also make a remark on the asymptotic behavior of global solutions. △ Less

Submitted 26 November, 2021; originally announced November 2021.

MSC Class: 35L70 (Primary); 35L52 (Secondary)

Journal ref: Journal of Evolution Equations volume 22, Article number: 49 (2022)

Moderate and $L^p$ maximal inequalities for diffusion processes and conformal martingales

Authors: Xian Chen, Yong Chen, Mumien Cheng, Chen Jia

Abstract: The $L^p$ maximal inequalities for martingales are one of the classical results in the theory of stochastic processes. Here we establish the sharp moderate maximal inequalities for one-dimensional diffusion processes, which include the $L^p$ maximal inequalities as special cases. Moreover, we apply our theory to many specific examples, including the Ornstein-Uhlenbeck (OU) process, Brownian motion… ▽ More The $L^p$ maximal inequalities for martingales are one of the classical results in the theory of stochastic processes. Here we establish the sharp moderate maximal inequalities for one-dimensional diffusion processes, which include the $L^p$ maximal inequalities as special cases. Moreover, we apply our theory to many specific examples, including the Ornstein-Uhlenbeck (OU) process, Brownian motion with drift, reflected Brownian motion with drift, Cox-Ingersoll-Ross process, radial OU process, and Bessel process. The results are further applied to establish the moderate maximal inequalities for some high-dimensional processes, including the complex OU process and general conformal local martingales. △ Less

Submitted 4 November, 2021; originally announced November 2021.

The second Bogolyubov theorem and global averaging principle for SPDEs with monotone coefficients

Authors: Mengyu Cheng, Zhenxin Liu

Abstract: In this paper, we establish the second Bogolyubov theorem and global averaging principle for stochastic partial differential equations (in short, SPDEs) with monotone coefficients. Firstly, we prove that there exists a unique $L^{2}$-bounded solution to SPDEs with monotone coefficients and this bounded solution is globally asymptotically stable in square-mean sense. Then we show that the $L^{2}$-b… ▽ More In this paper, we establish the second Bogolyubov theorem and global averaging principle for stochastic partial differential equations (in short, SPDEs) with monotone coefficients. Firstly, we prove that there exists a unique $L^{2}$-bounded solution to SPDEs with monotone coefficients and this bounded solution is globally asymptotically stable in square-mean sense. Then we show that the $L^{2}$-bounded solution possesses the same recurrent properties (e.g. periodic, quasi-periodic, almost periodic, almost automorphic, Birkhoff recurrent, Levitan almost periodic, etc.) in distribution sense as the coefficients. Thirdly, we prove that the recurrent solution of the original equation converges to the stationary solution of averaged equation under the compact-open topology as the time scale goes to zero--in other words, there exists a unique recurrent solution to the original equation in a neighborhood of the stationary solution of averaged equation when the time scale is small. Finally, we establish the global averaging principle in weak sense, i.e. we show that the attractor of original system tends to that of the averaged equation in probability measure space as the time scale goes to zero. For illustration of our results, we give two applications, including stochastic reaction diffusion equations and stochastic generalized porous media equations. △ Less

Submitted 9 August, 2022; v1 submitted 1 September, 2021; originally announced September 2021.

Comments: 41 pages

MSC Class: 70K65; 60H15; 37B20; 37L15

Combining the hybrid mimetic mixed method with the Scharfetter-Gummel scheme for magnetised transport in plasmas

Authors: Hanz Martin Cheng, Jan ten thije Boonkkamp, Jesper Janssen, Diana Mihailova, Jan van Dijk

Abstract: In this paper, we propose a numerical scheme for fluid models of magnetised plasmas. One important feature of the numerical scheme is that it should be able to handle the anisotropy induced by the magnetic field. In order to do so, we propose the use of the hybrid mimetic mixed (HMM) scheme for diffusion. This is combined with a hybridised variant of the Scharfetter-Gummel (SG) scheme for advectio… ▽ More In this paper, we propose a numerical scheme for fluid models of magnetised plasmas. One important feature of the numerical scheme is that it should be able to handle the anisotropy induced by the magnetic field. In order to do so, we propose the use of the hybrid mimetic mixed (HMM) scheme for diffusion. This is combined with a hybridised variant of the Scharfetter-Gummel (SG) scheme for advection. The proposed hybrid scheme can be implemented very efficiently via static condensation. Numerical tests are then performed to show the applicability of the combined HMM-SG scheme, even for highly anisotropic magnetic fields. △ Less

Submitted 8 September, 2022; v1 submitted 23 July, 2021; originally announced July 2021.

Singular metrics with negative scalar curvature

Authors: Man-Chuen Cheng, Man-Chun Lee, Luen-Fai Tam

Abstract: Motivated by the work of Li and Mantoulidis, we study singular metrics which are uniformly Euclidean $(L^\infty)$ on a compact manifold $M^n$ ($n\ge 3$) with negative Yamabe invariant $σ(M)$. It is well-known that if $g$ is a smooth metric on $M$ with unit volume and with scalar curvature $R(g)\ge σ(M)$, then $g$ is Einstein. We show, in all dimensions, the same is true for metrics with edge singu… ▽ More Motivated by the work of Li and Mantoulidis, we study singular metrics which are uniformly Euclidean $(L^\infty)$ on a compact manifold $M^n$ ($n\ge 3$) with negative Yamabe invariant $σ(M)$. It is well-known that if $g$ is a smooth metric on $M$ with unit volume and with scalar curvature $R(g)\ge σ(M)$, then $g$ is Einstein. We show, in all dimensions, the same is true for metrics with edge singularities with cone angles $\leq 2π$ along codimension-2 submanifolds. We also show in three dimension, if the Yamabe invariant of connected sum of two copies of $M$ attains its minimum, then the same is true for $L^\infty$ metrics with isolated point singularities. △ Less

Submitted 18 July, 2021; originally announced July 2021.

Comments: 27 pages

MSC Class: Primary 53C20

A fully local hybridised second-order accurate scheme for advection-diffusion equations

Authors: Hanz Martin Cheng

Abstract: In this paper, we present a fully local second-order upwind scheme, applicable on generic meshes. This is done by hybridisation, which is achieved by introducing unknowns on each edge of the mesh. By doing so, fluxes only depend on values associated to a single cell, and thus, this scheme can easily be applied even on cells near the boundary of the domain. Another advantage of hybridised schemes i… ▽ More In this paper, we present a fully local second-order upwind scheme, applicable on generic meshes. This is done by hybridisation, which is achieved by introducing unknowns on each edge of the mesh. By doing so, fluxes only depend on values associated to a single cell, and thus, this scheme can easily be applied even on cells near the boundary of the domain. Another advantage of hybridised schemes is that static condensation can be employed, leading to a very efficient implementation. A convergence analysis, which also covers a flux-limited TVD variant of the scheme, is then presented. Numerical results are also given in order to compare this with a hybridised first-order upwind scheme and a classical cell-centered second-order upwind type scheme. △ Less

Submitted 15 March, 2021; originally announced March 2021.

MSC Class: 65N08

A generalised complete flux scheme for anisotropic advection-diffusion equations

Authors: Hanz Martin Cheng, Jan ten Thije Boonkkamp

Abstract: In this paper, we consider separating the discretisation of the diffusive and advective fluxes in the complete flux scheme. This allows the combination of several discretisation methods for the homogeneous flux with the complete flux (CF) method. In particular, we explore the combination of the hybrid mimetic mixed (HMM) method and the CF method, in order to utilize the advantages of each of these… ▽ More In this paper, we consider separating the discretisation of the diffusive and advective fluxes in the complete flux scheme. This allows the combination of several discretisation methods for the homogeneous flux with the complete flux (CF) method. In particular, we explore the combination of the hybrid mimetic mixed (HMM) method and the CF method, in order to utilize the advantages of each of these methods. The usage of HMM allows us to handle anisotropic diffusion tensors on generic polygonal (polytopal) grids; whereas the CF method provides a framework for the construction of a uniformly second order method, even when the problem is advection dominated. △ Less

Submitted 22 December, 2020; v1 submitted 23 March, 2020; originally announced March 2020.

MSC Class: 65N08

Journal ref: Advances in Computational Mathematics 2021

Modular Exercises for Four-Point Blocks -- I

Authors: Miranda C. N. Cheng, Terry Gannon, Guglielmo Lockhart

Abstract: The well-known modular property of the torus characters and torus partition functions of (rational) vertex operator algebras (VOAs) and 2d conformal field theories (CFTs) has been an invaluable tool for studying this class of theories. In this work we prove that sphere four-point chiral blocks of rational VOAs are vector-valued modular forms for the groups $Γ(2)$, $Γ_0(2)$, or… ▽ More The well-known modular property of the torus characters and torus partition functions of (rational) vertex operator algebras (VOAs) and 2d conformal field theories (CFTs) has been an invaluable tool for studying this class of theories. In this work we prove that sphere four-point chiral blocks of rational VOAs are vector-valued modular forms for the groups $Γ(2)$, $Γ_0(2)$, or $\text{SL}_2(\mathbb{Z})$. Moreover, we prove that the four-point correlators, combining the holomorphic and anti-holomorphic chiral blocks, are modular invariant. In particular, in this language the crossing symmetries are simply modular symmetries. This gives the possibility of exploiting the available techniques and knowledge about modular forms to determine or constrain the physically interesting quantities such as chiral blocks and fusion coefficients, which we illustrate with a few examples. We also highlight the existence of a sphere-torus correspondence equating the sphere quantities of certain theories $\mathcal{T}_s$ with the torus quantities of another family of theories $\mathcal{T}_t$. A companion paper will delve into more examples and explore more systematically this sphere-torus duality. △ Less

Submitted 25 February, 2020; originally announced February 2020.

Comments: 72 pages, 4 figures

Three-Manifold Quantum Invariants and Mock Theta Functions

Authors: Miranda C. N. Cheng, Francesca Ferrari, Gabriele Sgroi

Abstract: Mock modular forms have found applications in numerous branches of mathematical sciences since they were first introduced by Ramanujan nearly a century ago. In this proceeding we highlight a new area where mock modular forms start to play an important role, namely the study of three-manifold invariants. For a certain class of Seifert three-manifolds, we describe a conjecture on the mock modular pr… ▽ More Mock modular forms have found applications in numerous branches of mathematical sciences since they were first introduced by Ramanujan nearly a century ago. In this proceeding we highlight a new area where mock modular forms start to play an important role, namely the study of three-manifold invariants. For a certain class of Seifert three-manifolds, we describe a conjecture on the mock modular properties of a recently proposed quantum invariant. As an illustration, we include concrete computations for a specific three-manifold, the Brieskorn sphere $Σ(2,3,7)$. This note is partially based on the talk by the first author in the conference "Srinivasa Ramanujan: in celebration of the centenary of his election as FRS" held at the Royal Society in 2018. △ Less

Submitted 16 March, 2020; v1 submitted 17 December, 2019; originally announced December 2019.

Comments: 19 pages

Journal ref: Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences. Volume 378, Issue 2163 (2019)

Periodic, almost periodic and almost automorphic solutions for SPDEs with monotone coefficients

Authors: Mengyu Cheng, Zhenxin Liu

Abstract: In this paper, we use the variational approach to investigate recurrent properties of solutions for stochastic partial differential equations, which is in contrast to the previous semigroup framework. Consider stochastic differential equations with monotone coefficients. Firstly, we establish the continuous dependence on initial values and coefficients for solutions. Secondly, we prove the existen… ▽ More In this paper, we use the variational approach to investigate recurrent properties of solutions for stochastic partial differential equations, which is in contrast to the previous semigroup framework. Consider stochastic differential equations with monotone coefficients. Firstly, we establish the continuous dependence on initial values and coefficients for solutions. Secondly, we prove the existence of recurrent solutions, which include periodic, almost periodic and almost automorphic solutions. Then we show that these recurrent solutions are globally asymptotically stable in square-mean sense. Finally, for illustration of our results we give two applications, i.e. stochastic reaction diffusion equations and stochastic porous media equations. △ Less

Submitted 5 November, 2019; originally announced November 2019.

A framework of constructing placement delivery arrays for centralized coded caching

Authors: Minquan Cheng, Jinyu Wang, Xi Zhong, Qiang Wang

Abstract: In caching system, it is desirable to design a coded caching scheme with the transmission load $R$ and subpacketization $F$ as small as possible, in order to improve efficiency of transmission in the peak traffic times and to decrease implementation complexity. Yan et al. reformulated the centralized coded caching scheme as designing a corresponding $F\times K$ array called placement delivery arra… ▽ More In caching system, it is desirable to design a coded caching scheme with the transmission load $R$ and subpacketization $F$ as small as possible, in order to improve efficiency of transmission in the peak traffic times and to decrease implementation complexity. Yan et al. reformulated the centralized coded caching scheme as designing a corresponding $F\times K$ array called placement delivery array (PDA), where $F$ is the subpacketization and $K$ is the number of users. Motivated by several constructions of PDAs, we introduce a framework for constructing PDAs, where each row is indexed by a row vector of some matrix called row index matrix and each column's index is labelled by an element of a direct product set. Using this framework, a new scheme is obtained, which can be regarded as a generalization of some previously known schemes. When $K$ is equal to ${m\choose t}q^t$ for positive integers $m$, $t$ with $t<m$ and $q\geq 2$, we show that the row index matrix must be an orthogonal array if all the users have the same memory size. Furthermore, the row index matrix must be a covering array if the coded gain is ${m\choose t}$, which is the maximal coded gain under our framework. Consequently the lower bounds on the transmission load and subpacketization of the schemes are derived under our framework. Finally, using orthogonal arrays as the row index matrix, we obtain two more explicit classes of schemes which have significantly advantages on the subpacketization while the transmission load is equal or close to that of the schemes constructed by Shangguan et al. (IEEE Trans. Inf. Theory, 64, 5755-5766, 2018) for the same number of users and memory size. △ Less

Submitted 30 May, 2021; v1 submitted 16 August, 2019; originally announced August 2019.

Comments: 13 pages

An efficient implementation of mass conserving characteristic-based schemes in 2D and 3D

Authors: Hanz Martin Cheng, Jerome Droniou

Abstract: In this paper, we develop the ball-approximated characteristics (B-char) method, which is an algorithm for efficiently implementing characteristic-based schemes in 2D and 3D. Core to the implementation of numerical schemes is the evaluation of integrals, which in the context of characteristic-based schemes with piecewise constant approximations boils down to computing the intersections between two… ▽ More In this paper, we develop the ball-approximated characteristics (B-char) method, which is an algorithm for efficiently implementing characteristic-based schemes in 2D and 3D. Core to the implementation of numerical schemes is the evaluation of integrals, which in the context of characteristic-based schemes with piecewise constant approximations boils down to computing the intersections between two regions. In the literature, these regions are approximated by polytopes (polygons in 2D and polyhedra in 3D) and, due to this, the implementation in 3D is nontrivial. The main novelty in this paper is the approximation of the regions by balls, whose intersections are much cheaper to compute than those of polytopes. Of course, balls cannot fully tessellate a region, and hence some mass may be lost. We perform some adjustments, and also solve an optimisation problem, in order to yield a scheme that is both locally and globally mass conserving. This algorithm can achieve results that are similar to those obtained from an implementation which uses polytopal intersections, with a much cheaper computational cost. △ Less

Submitted 22 January, 2020; v1 submitted 15 August, 2019; originally announced August 2019.

MSC Class: 65M08; 65M25

Probabilistic Existence Results for Parent-Identifying Schemes

Authors: Yujie Gu, Minquan Cheng, Grigory Kabatiansky, Ying Miao

Abstract: Parent-identifying schemes provide a way to identify causes from effects for some information systems such as digital fingerprinting and group testing. In this paper, we consider combinatorial structures for parent-identifying schemes. First, we establish an equivalent relationship between parent-identifying schemes and forbidden configurations. Based on this relationship, we derive probabilistic… ▽ More Parent-identifying schemes provide a way to identify causes from effects for some information systems such as digital fingerprinting and group testing. In this paper, we consider combinatorial structures for parent-identifying schemes. First, we establish an equivalent relationship between parent-identifying schemes and forbidden configurations. Based on this relationship, we derive probabilistic existence lower bounds for two related combinatorial structures, that is, $t$-parent-identifying set systems ($t$-IPPS) and $t$-multimedia parent-identifying codes ($t$-MIPPC), which are used in broadcast encryption and multimedia fingerprinting respectively. The probabilistic lower bound for the maximum size of a $t$-IPPS has the asymptotically optimal order of magnitude in many cases, and that for $t$-MIPPC provides the asymptotically optimal code rate when $t=2$ and the best known asymptotic code rate when $t\geq 3$. Furthermore, we analyze the structure of $2$-IPPS and prove some bounds for certain cases. △ Less

Submitted 3 June, 2019; originally announced June 2019.

Comments: 14 pages

On the moduli spaces of commuting elements in the projective unitary groups

Authors: Alejandro Adem, Man Chuen Cheng

Abstract: We provide descriptions for the moduli spaces $\text{Rep}(Γ, PU(m))$, where $Γ$ is any finitely generated abelian group and $PU(m)$ is the group of $m\times m$ projective unitary matrices. As an application we show that for any connected CW-complex $X$ with $π_1(X)\cong \mathbf{Z}^n$, the natural map $π_0(\text{Rep}(π_1(X), PU(m)))\to [X, BPU(m)]$ is injective, hence providing a complete enumerati… ▽ More We provide descriptions for the moduli spaces $\text{Rep}(Γ, PU(m))$, where $Γ$ is any finitely generated abelian group and $PU(m)$ is the group of $m\times m$ projective unitary matrices. As an application we show that for any connected CW-complex $X$ with $π_1(X)\cong \mathbf{Z}^n$, the natural map $π_0(\text{Rep}(π_1(X), PU(m)))\to [X, BPU(m)]$ is injective, hence providing a complete enumeration of the isomorphism classes of flat principal $PU(m)$-bundles over $X$. △ Less

Submitted 2 July, 2019; v1 submitted 9 March, 2019; originally announced March 2019.

Comments: 10 pages. Minor typos fixed. To appear in the Journal of Mathematical Physics

MSC Class: 20C25; 55N91; 55P10

An Efficient Solver for Cumulative Density Function-based Solutions of Uncertain Kinematic Wave Models

Authors: Ming Cheng, Yi Qin, Akil Narayan, Xinghui Zhong, Xueyu Zhu, Peng Wang

Abstract: We develop a numerical framework to implement the cumulative density function (CDF) method for obtaining the probability distribution of the system state described by a kinematic wave model. The approach relies on Monte Carlo Simulations (MCS) of the fine-grained CDF equation of system state, as derived by the CDF method. This fine-grained CDF equation is solved via the method of characteristics.… ▽ More We develop a numerical framework to implement the cumulative density function (CDF) method for obtaining the probability distribution of the system state described by a kinematic wave model. The approach relies on Monte Carlo Simulations (MCS) of the fine-grained CDF equation of system state, as derived by the CDF method. This fine-grained CDF equation is solved via the method of characteristics. Each method of characteristics solution is far more computationally efficient than the direct solution of the kinematic wave model, and the MCS estimator of the CDF converges relatively quickly. We verify the accuracy and robustness of our procedure via comparison with direct MCS of a particular kinematic wave system, the Saint-Venant equation. △ Less

Submitted 24 January, 2019; originally announced January 2019.

Comments: 19 pages, 6 figures, 3 tables

MSC Class: 65M22; 65M25

Journal ref: Journal of Computational Physics, 2018

3d Modularity

Authors: Miranda C. N. Cheng, Sungbong Chun, Francesca Ferrari, Sergei Gukov, Sarah M. Harrison

Abstract: We find and propose an explanation for a large variety of modularity-related symmetries in problems of 3-manifold topology and physics of 3d $\mathcal{N}=2$ theories where such structures a priori are not manifest. These modular structures include: mock modular forms, $SL(2,\mathbb{Z})$ Weil representations, quantum modular forms, non-semisimple modular tensor categories, and chiral algebras of lo… ▽ More We find and propose an explanation for a large variety of modularity-related symmetries in problems of 3-manifold topology and physics of 3d $\mathcal{N}=2$ theories where such structures a priori are not manifest. These modular structures include: mock modular forms, $SL(2,\mathbb{Z})$ Weil representations, quantum modular forms, non-semisimple modular tensor categories, and chiral algebras of logarithmic CFTs. △ Less

Submitted 27 May, 2020; v1 submitted 26 September, 2018; originally announced September 2018.

Comments: 119 pages, 10 figures and 20 tables

Report number: CALT-TH-2018-037

TASI Lectures on Moonshine

Authors: Vassilis Anagiannis, Miranda C. N. Cheng

Abstract: The word moonshine refers to unexpected relations between the two distinct mathematical structures: finite group representations and modular objects. It is believed that the key to understanding moonshine is through physical theories with special symmetries. Recent years have seen a varieties of new ways in which finite group representations and modular objects can be connected to each other, and… ▽ More The word moonshine refers to unexpected relations between the two distinct mathematical structures: finite group representations and modular objects. It is believed that the key to understanding moonshine is through physical theories with special symmetries. Recent years have seen a varieties of new ways in which finite group representations and modular objects can be connected to each other, and these developments have brought promises and also puzzles into the string theory community. These lecture notes aim to bring graduate students in theoretical physics and mathematical physics to the forefront of this active research area. In Part II of this note, we review the various cases of moonshine connections, ranging from the classical monstrous moonshine established in the last century to the most recent findings. In Part III, we discuss the relation between the moonshine connections and physics, especially string theory. After briefly reviewing a recent physical realisation of monstrous moonshine, we will describe in some details the mystery of the physical aspects of umbral moonshine, and also mention some other setups where string theory black holes can be connected to moonshine. To make the exposition self-contained, we also provide the relevant background knowledge in Part I, including sections on finite groups, modular objects, and two-dimensional conformal field theories. This part occupies half of the pages of this set of notes and can be skipped by readers who are already familiar with the relevant concepts and techniques. △ Less

Submitted 2 July, 2018; originally announced July 2018.

Comments: summer school lecture notes, 78 pages

Stochastic Zeroth-order Optimization via Variance Reduction method

Authors: Liu Liu, Minhao Cheng, Cho-Jui Hsieh, Dacheng Tao

Abstract: Derivative-free optimization has become an important technique used in machine learning for optimizing black-box models. To conduct updates without explicitly computing gradient, most current approaches iteratively sample a random search direction from Gaussian distribution and compute the estimated gradient along that direction. However, due to the variance in the search direction, the convergenc… ▽ More Derivative-free optimization has become an important technique used in machine learning for optimizing black-box models. To conduct updates without explicitly computing gradient, most current approaches iteratively sample a random search direction from Gaussian distribution and compute the estimated gradient along that direction. However, due to the variance in the search direction, the convergence rates and query complexities of existing methods suffer from a factor of $d$, where $d$ is the problem dimension. In this paper, we introduce a novel Stochastic Zeroth-order method with Variance Reduction under Gaussian smoothing (SZVR-G) and establish the complexity for optimizing non-convex problems. With variance reduction on both sample space and search space, the complexity of our algorithm is sublinear to $d$ and is strictly better than current approaches, in both smooth and non-smooth cases. Moreover, we extend the proposed method to the mini-batch version. Our experimental results demonstrate the superior performance of the proposed method over existing derivative-free optimization techniques. Furthermore, we successfully apply our method to conduct a universal black-box attack to deep neural networks and present some interesting results. △ Less

Submitted 2 August, 2018; v1 submitted 30 May, 2018; originally announced May 2018.

A combined GDM--ELLAM--MMOC scheme for advection dominated PDEs

Authors: Hanz Martin Cheng, Jérôme Droniou, Kim-Ngan Le

Abstract: We propose a combination of the Eulerian Lagrangian Localised Adjoint Method (ELLAM) and the Modified Method of Characteristics (MMOC) for time-dependent advection-domina\-ted PDEs. The combined scheme, so-called GEM scheme, takes advantages of both ELLAM scheme (mass conservation) and MMOC scheme (easier computations), while at the same time avoids their disadvantages (respectively, harder tracki… ▽ More We propose a combination of the Eulerian Lagrangian Localised Adjoint Method (ELLAM) and the Modified Method of Characteristics (MMOC) for time-dependent advection-domina\-ted PDEs. The combined scheme, so-called GEM scheme, takes advantages of both ELLAM scheme (mass conservation) and MMOC scheme (easier computations), while at the same time avoids their disadvantages (respectively, harder tracking around the injection regions, and loss of mass). We present a precise analysis of mass conservation properties for these three schemes, and after achieving global mass balance, an adjustment yielding local volume conservation is then proposed. Numerical results for all three schemes are then compared, illustrating the advantages of the GEM scheme. A convergence result of the MMOC scheme, motivated by our previous work on the convergence of ELLAM schemes, is provided, which can be extended to obtain the convergence of GEM scheme. △ Less

Submitted 13 October, 2021; v1 submitted 15 May, 2018; originally announced May 2018.

Convergence analysis of a family of ELLAM schemes for a fully coupled model of miscible displacement in porous media

Authors: Hanz Martin Cheng, Jerome Droniou, Kim-Ngan Le

Abstract: We analyse the convergence of numerical schemes in the GDM-ELLAM (Gradient Discretisation Method-Eulerian Lagrangian Localised Adjoint Method) framework for a strongly coupled elliptic-parabolic PDE which models miscible displacement in porous media. These schemes include, but are not limited to Mixed Finite Element-ELLAM and Hybrid Mimetic Mixed-ELLAM schemes. A complete convergence analysis is p… ▽ More We analyse the convergence of numerical schemes in the GDM-ELLAM (Gradient Discretisation Method-Eulerian Lagrangian Localised Adjoint Method) framework for a strongly coupled elliptic-parabolic PDE which models miscible displacement in porous media. These schemes include, but are not limited to Mixed Finite Element-ELLAM and Hybrid Mimetic Mixed-ELLAM schemes. A complete convergence analysis is presented on the coupled model, using only weak regularity assumptions on the solution (which are satisfied in practical applications), and not relying on $L^\infty$ bounds (which are impossible to ensure at the discrete level given the anisotropic diffusion tensors and the general grids used in applications). △ Less

Submitted 14 November, 2017; v1 submitted 5 October, 2017; originally announced October 2017.

K3 Elliptic Genus and an Umbral Moonshine Module

Authors: Vassilis Anagiannis, Miranda C. N. Cheng, Sarah M. Harrison

Abstract: Umbral moonshine connects the symmetry groups of the 23 Niemeier lattices with 23 sets of distinguished mock modular forms. The 23 cases of umbral moonshine have a uniform relation to symmetries of $K3$ string theories. Moreover, a supersymmetric vertex operator algebra with Conway sporadic symmetry also enjoys a close relation to the $K3$ elliptic genus. Inspired by the above two relations betwee… ▽ More Umbral moonshine connects the symmetry groups of the 23 Niemeier lattices with 23 sets of distinguished mock modular forms. The 23 cases of umbral moonshine have a uniform relation to symmetries of $K3$ string theories. Moreover, a supersymmetric vertex operator algebra with Conway sporadic symmetry also enjoys a close relation to the $K3$ elliptic genus. Inspired by the above two relations between moonshine and $K3$ string theory, we construct a chiral CFT by orbifolding the free theory of 24 chiral fermions and two pairs of fermionic and bosonic ghosts. In this paper we mainly focus on the case of umbral moonshine corresponding to the Niemeier lattice with root system given by 6 copies of $D_4$ root system. This CFT then leads to the construction of an infinite-dimensional graded module for the umbral group $G^{D_4^{\oplus 6}}$ whose graded characters coincide with the umbral moonshine functions. We also comment on how one can recover all umbral moonshine functions corresponding to the Niemeier root systems $A_5^{\oplus 4}D_4$, $A_7^{\oplus 2}D_5^{\oplus 2}$ , $A_{11}D_7 E_6$, $A_{17}E_7$, and $D_{10}E_7^{\oplus 2}$. △ Less

Submitted 6 September, 2017; originally announced September 2017.

Comments: 41 pages, 3 appendices

Attractive Strings and Five-Branes, Skew-Holomorphic Jacobi Forms and Moonshine

Authors: Miranda C. N. Cheng, John F. R. Duncan, Sarah M. Harrison, Jeffrey A. Harvey, Shamit Kachru, Brandon C. Rayhaun

Abstract: We show that certain BPS counting functions for both fundamental strings and strings arising from fivebranes wrapping divisors in Calabi--Yau threefolds naturally give rise to skew-holomorphic Jacobi forms at rational and attractor points in the moduli space of string compactifications. For M5-branes wrapping divisors these are forms of weight negative one, and in the case of multiple M5-branes sk… ▽ More We show that certain BPS counting functions for both fundamental strings and strings arising from fivebranes wrapping divisors in Calabi--Yau threefolds naturally give rise to skew-holomorphic Jacobi forms at rational and attractor points in the moduli space of string compactifications. For M5-branes wrapping divisors these are forms of weight negative one, and in the case of multiple M5-branes skew-holomorphic mock Jacobi forms arise. We further find that in simple examples these forms are related to skew-holomorphic (mock) Jacobi forms of weight two that play starring roles in moonshine. We discuss examples involving M5-branes on the complex projective plane, del Pezzo surfaces of degree one, and half-K3 surfaces. For del Pezzo surfaces of degree one and certain half-K3 surfaces we find a corresponding graded (virtual) module for the degree twelve Mathieu group. This suggests a more extensive relationship between Mathieu groups and complex surfaces, and a broader role for M5-branes in the theory of Jacobi forms and moonshine. △ Less

Submitted 23 July, 2018; v1 submitted 24 August, 2017; originally announced August 2017.

Comments: 36 pages, LaTeX; minor typos corrected, footnote added at bottom of page 9 to accommodate JHEP editor's suggestion