Time-dependent Hermite-Galerkin spectral method and its applications

A time-dependent Hermite-Galerkin spectral method (THGSM) is investigated in this paper for the nonlinear convection-diffusion equations in the unbounded domains. The time-dependent scaling factor and translating factor are introduced in the definition of the generalized Hermite functions (GHF). As a consequence, the THGSM based on these GHF has many advantages, not only in theorethical proofs, but also in numerical implementations. The stability and spectral convergence of our proposed method have been established in this paper. The Korteweg-de Vries-Burgers (KdVB) equation and its special cases, including the heat equation and the Burgers' equation, as the examples, have been numerically solved by our method. The numerical results are presented, and it surpasses the existing methods in accuracy. Our theoretical proof of the spectral convergence has been supported by the numerical results.Comment: 16 pages, 7 table

Whole genome single nucleotide polymorphism genotyping of Staphylococcus aureus

Next-generation sequencing technology enables routine detection of bacterial pathogens for clinical diagnostics and genetic research. Whole genome sequencing has been of importance in the epidemiologic analysis of bacterial pathogens. However, few whole genome sequencing-based genotyping pipelines are available for practical applications. Here, we present the whole genome sequencing-based single nucleotide polymorphism (SNP) genotyping method and apply to the evolutionary analysis of methicillin-resistant Staphylococcus aureus. The SNP genotyping method calls genome variants using next-generation sequencing reads of whole genomes and calculates the pair-wise Jaccard distances of the genome variants. The method may reveal the high-resolution whole genome SNP profiles and the structural variants of different isolates of methicillin-resistant S. aureus (MRSA) and methicillin-susceptible S. aureus (MSSA) strains. The phylogenetic analysis of whole genomes and particular regions may monitor and track the evolution and the transmission dynamic of bacterial pathogens. The computer programs of the whole genome sequencing-based SNP genotyping method are available to the public at https://github.com/cyinbox/NGS.Comment: 6 figures, 2 table

Hermite spectral method to 1D forward Kolmogorov equation and its application to nonlinear filtering problems

In this paper, we investigate the Hermite spectral method (HSM) to numerically solve the forward Kolmogorov equation (FKE). A useful guideline of choosing the scaling factor of the generalized Hermite functions is given in this paper. It greatly improves the resolution of HSM. The convergence rate of HSM to FKE is analyzed in the suitable function space and has been verified by the numerical simulation. As an important application and our primary motivation to study the HSM to FKE, we work on the implementation of the nonlinear filtering (NLF) problem with a real-time algorithm developed in [17]. The HSM to FKE is served as the off-line computation in this algorithm. The translating factor of the generalized Hermite functions and the moving-window technique are introduced to deal with the drifting of the posterior conditional density function of the states in the on-line experiments. Two numerical experiments of NLF problems are carried out to illustrate the feasibility of our algorithm. Moreover, our algorithm surpasses the particle filter as a real-time solver to NLF.Comment: 14 pages, 6 figures, 1 tabl

On the quenching behavior of the MEMS with fringing field

The singular parabolic problem $u_t-\triangle u=\lambda{\frac{1+\delta|\nabla u|^2}{(1-u)^2}}$ on a bounded domain $\Omega$ of $\mathbb{R}^n$ with Dirichlet boundary condition, models the Microelectromechanical systems (MEMS) device with fringing field. In this paper, we focus on the quenching behavior of the solution to this equation. We first show that there exists a critical value $\lambda_\delta^*>0$ such that if $0<\lambda<\lambda_\delta^*$, all solutions exist globally; while for $\lambda>\lambda_\delta^*$, all the solution will quench in finite time. The estimate of the quenching time in terms of large voltage $\lambda$ is investigated. Furthermore, the quenching set is a compact subset of $\Omega$, provided $\Omega$ is a convex bounded domain in $\mathbb{R}^n$. In particular, if the domain $\Omega$ is radially symmetric, then the origin is the only quenching point. We not only derive the one-side estimate of the quenching rate, but also further study the refined asymptotic behavior of the finite quenching solution.Comment: 27 pages, 3 figures, 3 tables; accepted by Quart. Appl. Math. arXiv admin note: text overlap with arXiv:0712.3071 by other author

Classification of 3-dimensional isolated rational hypersurface singularities with C*-action

In the paper "Algebraic classification of rational CR structures on topological 5-sphere with transversal holomorphic S^1-action in C^4" (Yau and Yu, Math. Nachrichten 246-247(2002), 207-233), we give algebraic classification of rational CR structures on the topological 5-sphere with transversal holomorphic S^1-action in C^4. Here, algebraic classification of compact strongly pseudoconvex CR manifolds X means classification up to algebraic equivalence, i.e. roughly up to isomorphism of the normalization of the complex analytic variety V which has X as boundary. The problem is intimately related to the study of 3-dimensional isolated rational weighted homogeneous hypersurface singularities with link homeomorphic to S^5. For this, we need the classification of 3-dimensional isolated rational hypersurface singularities with a C*-action. This list is only available at the homepage of one of us. Since there is a desire for a complete list of this classification (cf. Theorem 3.3), we decide to publish it for the convenience of readers

4d N=2 SCFT from Complete Intersection Singularity

Detailed studies of four dimensional N=2 superconformal field theories (SCFT) defined by isolated complete intersection singularities are performed: we compute the Coulomb branch spectrum, Seiberg-Witten solutions and central charges. Most of our theories have exactly marginal deformations and we identify the weakly coupled gauge theory descriptions for many of them, which involve (affine) D and E shaped quiver gauge theories and theories formed from Argyres-Douglas matters. These investigations provide strong evidence for the singularity approach in classifying 4d N=2 SCFTs.Comment: 46 pages, 85 figure

Counterexample to boundary regularity of a strongly pseudoconvex CR submanifold: An addendum to the paper of Harvey-Lawson

The purpose of this paper is to give a counterexample of Theorem 10.4 in [Ann. of Math. 102 (1975), 223-290]. In the Harvey-Lawson paper, a global result is claimed, but only a local result is proven. This theorem has had a big impact on CR geometry for almost a quarter of a century because one can use the theory of isolated singularities to study the theory of CR manifolds and vice versa.Comment: 2 pages, published versio

4d N=2 SCFT and singularity theory Part II: Complete intersection

We classify three dimensional isolated weighted homogeneous rational complete intersection singularities, which define many new four dimensional N=2 superconformal field theories. We also determine the mini-versal deformation of these singularities, and therefore solve the Coulomb branch spectrum and Seiberg-Witten solution.Comment: 64 pages, 1 figur

4d N=2 SCFT and singularity theory Part III: Rigid singularity

We classify three fold isolated quotient Gorenstein singularity $C^3/G$. These singularities are rigid, i.e. there is no non-trivial deformation, and we conjecture that they define 4d $\mathcal{N}=2$ SCFTs which do not have a Coulomb branch.Comment: 17 pages, 1 figur

Complete Weight Distribution and MacWilliams Identities for Asymmetric Quantum Codes

In 1997, Shor and Laflamme defined the weight enumerators for quantum error-correcting codes and derived a MacWilliams identity. We extend their work by introducing our double weight enumerators and complete weight enumerators. The MacWilliams identities for these enumerators can be obtained similarly. With the help of MacWilliams identities, we obtain various bounds for asymmetric quantum codes.Comment: 15 page