This article focuses on numerically studying the eigenstructure behavior of generalized eigenvalue problems (GEPs) arising in three dimensional (3D) source-free Maxwell's equations with magnetoelectric coupling effects which model 3D reciprocal chiral media. It is challenging to solve such a large-scale GEP efficiently. We combine the null-space free method with the inexact shift-invert residual Arnoldi method and MINRES linear solver to solve the GEP with a matrix dimension as large as 5,308,416. The eigenstructure is heavily determined by the chirality parameter $\gamma$. We show that all the eigenvalues are real and finite for a small chirality $\gamma$. For a critical value $\gamma = \gamma^*$, the GEP has $2 \times 2$ Jordan blocks at infinity eigenvalues. Numerical results demonstrate that when $\gamma$ increases from $\gamma^*$, the $2 \times 2$ Jordan block will first split into a complex conjugate eigenpair, then rapidly collide with the real axis and bifurcate into positive (resonance) and negative eigenvalues with modulus smaller than the other existing positive eigenvalues. The resonance band also exhibits an anticrossing interaction. Moreover, the electric and magnetic fields of the resonance modes are localized inside the structure, with only a slight amount of field leaking into the background (dielectric) material.

Szemerédi's regularity lemma and its variants are some of the most powerful tools in combinatorics. In this paper, we establish several results around the regularity lemma. First, we prove that whether or not we include the condition that the desired vertex partition in the regularity lemma is equitable has a minimal effect on the number of parts of the partition. Second, we use an algorithmic version of the (weak) Frieze–Kannan regularity lemma to give a substantially faster deterministic approximation algorithm for counting subgraphs in a graph. Previously, only an exponential dependence for the running time on the error parameter was known, and we improve it to a polynomial dependence. Third, we revisit the problem of finding an algorithmic regularity lemma, giving approximation algorithms for several co-NP-complete problems. We show how to use the weak Frieze–Kannan regularity lemma to approximate the regularity of a pair of vertex subsets. We also show how to quickly find, for each ε′>ε, an ε′-regular partition with k parts if there exists an ε-regular partition with k parts. Finally, we give a simple proof of the permutation regularity lemma which improves the tower-type bound on the number of parts in the previous proofs to a single exponential bound.

In this note, we prove the sharp Davies-Gaffney-Grigor’yan lemma for minimal heat kernels on graphs.

We study the Newton polygon jumping locus of a Mumford family in char p. Our main result says that, under a mild assumption on p, the jumping locus consists of only supersingular points and its cardinality is equal to (p^r − 1)(g − 1), where r is the degree of the defining field of the base curve of a Mumford family in char p and g is the genus of the curve. The underlying technique is the p-adic Hodge theory.

For each natural odd number n ≥ 3, we exhibit a maximal family of n-dimensional Calabi–Yau manifolds whose Yukawa coupling length is 1. As a consequence, Shafarevich’s conjecture holds true for these families. Moreover, it follows from Deligne and Mostow (Publ. Math. IHÉS, 63:5–89, 1986) and Mostow (Publ. Math. IHÉS, 63:91–106, 1986; J. Am. Math. Soc., 1(3):555–586, 1988) that, for n = 3, it can be partially compactified to a Shimura family of ball type, and for n = 5,9, there is a sub Q-PVHS of the family uniformizing a Zariski open subset of an arithmetic ball quotient

We show that $\lim_{t \to 0} e^{it\Delta}f(x) = f(x)$ almost everywhere for all $f \in H^s (\mathbb{R}^2)$ provided that $s>1/3$. This result is sharp up to the endpoint. The proof uses polynomial partitioning and decoupling.

Given a compact Riemannian manifold $M$, we consider a warped product manifold $\bar M = I \times_h M$, where $I$ is an open interval in $\Bbb R$. For a positive function $\psi$ defined on $\bar M$, we generalize the arguments in \cite{GRW2015} and \cite{RW16}, to obtain the curvature estimates for Hessian equations $\sigma_k(\kappa)=\psi(V,\nu(V))$. We also obtain some existence results for the starshaped compact hypersurface $\Sigma$ satisfying the above equation with various assumptions.

We consider the quermassintegral preserving flow of closed \emph{h-convex} hypersurfaces in hyperbolic space with the speed given by any positive power of a smooth symmetric, strictly increasing, and homogeneous of degree one function $f$ of the principal curvatures which is inverse concave and has dual $f_*$ approaching zero on the boundary of the positive cone. We prove that if the initial hypersurface is \emph{h-convex}, then the solution of the flow becomes strictly \emph{h-convex} for $t>0$, the flow exists for all time and converges to a geodesic sphere exponentially in the smooth topology.

No abstract uploaded!

No abstract uploaded!

In this paper, we establish various L 2 -estimates for the exterior differential operator on p-convex Riemannian manifolds in the sense of Harvey and Lawson. As applications, we establish a Carleman type estimate which is uniform with respect to both weight functions and domains, and we also obtain topological restrictions for a Riemannian manifold to be p-convex.

We supply a detailed proof of the result by P.S. Green and T. H$\ddot{\text{u}}$bsch that all complete intersection Calabi--Yau 3-folds in product of projective spaces are connected through projective conifold transitions (known as the standard web). We also introduce a subclass of small transitions which we call \emph{primitive} small transitions and study such subclass. More precisely, given a small projective resolution $\pi : \widehat{X} \rightarrow X$ of a Calabi--Yau 3-fold $X$, we show that if the natural closed immersion $\Def(\widehat{X}) \hookrightarrow \Def(X)$ is an isomorphism then $X$ has only ODPs as singularities.

We apply quantum gravitational results to spatially unbounded Friedmann universes and try to answer some questions related to dark energy, dark matter, inflation and the missing antimatter.

We discuss Donaldson-Thomas (DT) invariants of torsion sheaves with 2 dimensional support on a smooth projective surface in an ambient non-compact Calabi Yau fourfold given by the total space of a rank 2 bundle on the surface. We prove that in certain cases, when the rank 2 bundle is chosen appropriately, the universal truncated Atiyah class of these codimension 2 sheaves reduces to one, defined over the moduli space of such sheaves realized as torsion codimension 1 sheaves in a noncompact divisor (threefold) embedded in the ambient fourfold. Such reduction property of universal Atiyah class enables us to relate our fourfold DT theory to a reduced DT theory of a threefold and subsequently then to the moduli spaces of sheaves on the base surface using results in arXiv:1701.08899 and arXiv:1701.08902 and . We finally make predictions about modularity of such fourfold invariants when the base surface is an elliptic K3.

We propose a novel framework for reconstructing lightweight polygonal surfaces from point clouds. Unlike traditional methods that focus on either extracting good geometric primitives or obtaining proper arrangements of primitives, the emphasis of this work lies in intersecting the primitives (planes only) and seeking for an appropriate combination of them to obtain a manifold polygonal surface model without boundary. We show that reconstruction from point clouds can be cast as a binary labeling problem. Our method is based on a hypothesizing and selection strategy. We first generate a reasonably large set of face candidates by intersecting the extracted planar primitives. Then an optimal subset of the candidate faces is selected through optimization. Our optimization is based on a binary linear programming formulation under hard constraints that enforce the final polygonal surface model to be manifold and watertight. Experiments on point clouds from various sources demonstrate that our method can generate lightweight polygonal surface models of arbitrary piecewise planar objects. Besides, our method is capable of recovering sharp features and is robust to noise, outliers, and missing data.

We prove trace class estimates for self-adjoint elliptic operators defined in $\R[n]$ or $\R[]_+$. These estimates are also applicable when a physical system is governed by a wave equation by employing separation of variables to obtain corresponding temporal and spatial Hamiltonians. It is shown, in one important example, that the resulting Hamiltonians are of trace class such that quantum statistics can be applied to the system.

Full spectrum of a large sparse $\top$-palindromic quadratic eigenvalue problem ($\top$-PQEP) is considered arguably for the first time in this article. Such a problem is posed by calculation of surface Green's functions (SGFs) of mesoscopic transistors with a tremendous non-periodic cross-section. For this problem, general purpose eigensolvers are not efficient, nor is advisable to resort to the decimation method \etc to obtain the Wiener-Hopf factorization. After reviewing some rigorous understanding of SGF calculation from the perspective of $\top$-PQEP and nonlinear matrix equation, we present our new approach to this problem. In a nutshell, the unit disk where the spectrum of interest lies is broken down adaptively into pieces small enough that they each can be locally tackled by the generalized $\top$-skew-Hamiltonian implicitly restarted shift-and-invert Arnoldi (G$\top$SHIRA) algorithm with suitable shifts and other parameters, and the eigenvalues missed by this divide-and-conquer strategy can be recovered thanks to the accurate estimation provided by our newly developed scheme. Notably the novel non-equivalence deflation is proposed to avoid as much as possible duplication of nearby known eigenvalues when a new shift of G$\top$SHIRA is determined. We demonstrate our new approach by calculating the SGF of a realistic nanowire whose unit cell is described by a matrix of size $4000\times 4000$ at the density functional tight binding level, corresponding to a $8\times 8\; \mbox{nm}^{2}$ cross-section. We believe that quantum transport simulation of realistic nano-devices in the mesoscopic regime will greatly benefit from this work.

Calculation of band structures of three dimensional photonic crystals amounts to solving large-scale Maxwell eigenvalue problems, which are notoriously challenging due to high multiplicity of zero eigenvalues. In this paper, we try to address this problem in such a broad context that band structures of three dimensional isotropic photonic crystals in all $14$ Bravais lattices can be efficiently computed in a unified framework. In this work, we uncover the delicate machinery behind several key results of our framework and on the basis of this new understanding we drastically simplify the derivations, proofs and arguments. Particular effort is made on reformulating the Bloch condition for all $14$ Bravais lattices in the redefined orthogonal coordinate system, and establishing eigen-decomposition of discrete partial derivative operators by identifying the hierarchical structure of the underlying normal (block) companion matrix, and reducing the eigen-decomposition of the double-curl operator to a simple factorization of a $3$-by-$3$ complex skew-symmetric matrix. With the validity of the novel nullspace free method in the broad context, we perform some calculations on one benchmark system to demonstrate the accuracy and efficiency of our algorithm to solve Maxwell eigenvalue problems.

We study the heat kernel expansion of the Laplacian on n-forms defined on a subgraph of a directed complete graph. We derive two expressions for the subgraph heat kernel on 0-forms and compute the coefficients of the expansion. We also obtain the subgraph heat kernel of the Laplacian on 1-forms.

We set up a framework to study one-dimensional heat equations defined by fractal Laplacians associated with self-similar measures with overlaps. We show that for a class of such self-similar measures, a heat equation can be discretized and the finite element method can be applied to yield a system of linear differential equations. We show that the numerical solutions converge to the actual solution and obtain the rate of convergence. We also study some properties of the solutions of the heat equation.

We study infinite iterated functions systems (IIFSs) consisting of bi-Lipschitz mappings instead of conformal contractions, focusing on IFSs that do not satisfy the open set condition. By assuming the logarithmic distortion property and some cardinality growth condition, we obtain a formula for the Hausdorff, box, and packing dimensions of the limit set in terms of certain topological pressure. By assuming, in addition, the weak separation condition, we show that these dimensions are equal to the growth dimension of the limit set.

No abstract uploaded!

Motivated by Kapranov's discovery of an $L_\infty$ algebra structure on the tangent complex of a K\"{a}hler manifold and Chen-Sti\'{e}non-Xu's construction of a Leibniz$_\infty[1]$ algebra associated with a Lie pair, we find a general method to construct Leibniz$_\infty[1]$ algebras ---from a DG derivation $\mathscr{A} \xrightarrow{\delta} \Omega$ of a commutative differential graded algebra $\mathscr{A}$ valued in a DG $\mathscr{A}$-module $\Omega$. We prove that for any $\delta$-connection $\nabla$ on $\mathcal{B}$, the $\A$-dual of $\Omega$, there associates a Leibniz$_\infty[1]$ $\mathscr{A}$-algebra $(\mathcal{B},\{\mathcal{R}^\nabla_k\}_{k\geq 1})$. Moreover, this construction is canonical, i.e., the isomorphism class of $(\mathcal{B},\{\mathcal{R}^\nabla_k\}_{k\geq 1})$ only depends on the homotopy class of $\delta$.

In this paper, we study the Atiyah class and Todd class of the DG manifold $(F[1],d_F)$ corresponding to an integrable distribution $F \subset T_\mathbb{K} M = TM \otimes_{\mathbb{R}} \mathbb{K}$, where $\mathbb{K} = \mathbb{R}$ or $\\mathbb{C}$. We show that these two classes are canonically identical to those of the Lie pair $(T_{\mathbb{K}} M, F)$. As a consequence, the Atiyah class of a complex manifold $X$ is isomorphic to the Atiyah class of the corresponding DG manifold $(T^{0,1}_X[1],\bar{\partial})$. Moreover, if $X$ is a compact K\"{a}hler manifold, then the Todd class of $X$ is also isomorphic to the Todd class of the corresponding DG manifold $(T^{0,1}_X[1],\bar{\partial})$.

The doubling algorithms are very efficient iterative methods for computing the unique minimal nonnegative solution to an M-matrix algebraic Riccati equation (MARE). They are globally and quadratically convergent, except for MARE in the critical case where convergence is linear with the linear rate 1/2. However, the initialization phase and the doubling iteration kernel of any doubling algorithm involve inverting nonsingular M-matrices. In particular, for MARE in the critical case, the M-matrices in the doubling iteration kernel, although nonsingular, move towards singular M-matrices at convergence. These inversions are causes of concerns on entrywise relative accuracy of the eventually computed minimal nonnegative solution. Fortunately, a nonsingular M-matrix can be inverted by the GTH-like algorithm of Alfa, Xue, and Ye [Math. Comp., 71:217--236, 2002] to almost full entrywise relative accuracy, provided a triplet representation of the matrix is known. Recently, Nguyen and Poloni [Numer. Math., 130(4):763--792, 2015] discovered a way to construct triplet representations in a cancellation-free manner for all involved M-matrices in the doubling iteration kernel, for a special class of MAREs arising from Markov-modulated fluid queues. In this paper, we extend Nguyen's and Poloni's work to all MAREs by also devising a way to construct the triplet representations cancellation-free. Our construction, however, is not a straightforward extension of theirs. It is made possible by an introduction of novel recursively computable auxiliary nonnegative vectors. As the second contribution, we propose an entrywise relative residual for an approximate solution. The residual has an appealing feature of being able to reveal the entrywise relative accuracies of all entries, large and small, of the approximation. This is in marked contrast to the usual legacy normalized residual which reflects relative accuracies of large entries well but not so much those of very tiny entries. Numerical examples are presented to demonstrate and confirm our claims.