In this paper we consider operators of the form$H$=λ(-$i$∇), with λ analytic in a strip and with some specific growth conditions at infinity, and prove Hardy type estimates in$L$^{2}(ℝ^{$n$}) with exponential weights. In fact we extend our previous results [19] from bounded analytic functions on a strip to analytic functions with polynomial growth in that strip.

This paper deals with a degenerate diffusion Patlak–Keller–Segel system in $n\geq 3$ dimension. The main difference between the current work and many other recent studies on the same model is that we study the diffusion exponent $m=2n/(n+2)$, which is smaller than the usual exponent $m^{*}=2-2/n$ used in other studies. With the exponent $m=2n/(n+2)$, the associated free energy is conformal invariant, and there is a family of stationary solutions $U_{\lambda,x_0}(x)=C(\frac{\lambda} {\lambda^2+|x-x_0|^2})^{\frac{n+2}{2}}$ $\forall \lambda>0$, $x_0\in {\mathbb R}^n$. For radially symmetric solutions, we prove that if the initial data are strictly below $U_{\lambda,0}(x)$ for some $\lambda$, then the solution vanishes in $L^1_{loc}$ as $t\to\infty$; if the initial data are strictly above $U_{\lambda,0}(x)$ for some $\lambda$, then the solution either blows up at a finite time or has a mass concentration at $r=0$ as time goes to infinity. For general initial data, we prove that there is a global weak solution provided that the $L^m$ norm of initial density is less than a universal constant, and the weak solution vanishes as time goes to infinity. We also prove a finite time blow-up of the solution if the $L^m$ norm for initial data is larger than the $L^m$ norm of $U_{\lambda,x_0}(x)$, which is constant independent of $\lambda$ and $x_0$, and the free energy of initial data is smaller than that of $U_{\lambda,x_0}(x)$.

In this paper, we develop high-order bound-preserving (BP) discontinuous Galerkin (DG) methods for the coupled system of compressible miscible displacements on triangular meshes. We consider the problem with multi-component fluid mixture and the (volumetric) concentration of the $j$th component, $c_j$, should be between 0 and 1. There are three main difficulties. Firstly, $c_j$ does not satisfy a maximum-principle. Therefore, the numerical techniques introduced in (X. Zhang and C.-W. Shu, Journal of Computational Physics, 229 (2010), 3091-3120) cannot be applied directly. The main idea is to apply the positivity-preserving techniques to all $c_j's$ and enforce $\sum_jc_j=1$ simultaneously to obtain physically relevant approximations. By doing so, we have to treat the time derivative of the pressure $dp/dt$ as a source in the concentration equation and choose suitable fluxes in the pressure and concentration equations. Secondly, it is not easy to construct first-order numerical fluxes for interior penalty DG methods on triangular meshes. One of the key points in the high-order BP technique applied in this paper is the combination of high-order and lower-order numerical fluxes. We will construct second-order BP schemes and use the second-order numerical fluxes as the lower-order one. Finally, the classical slope limiter cannot be applied to $c_j$. To construct the BP technique, we will not approximate $c_j$ directly. Therefore, a new limiter will be introduced. Numerical experiments will be given to demonstrate the high-order accuracy and good performance of the numerical technique.

Wave dynamics in topological materials has been widely studied recently. A striking feature is the existence of robust and chiral wave propagations that have potential applications in many fields. A common way to realize such wave patterns is to utilize Dirac points which carry topological indices and is supported by the symmetries of the media. In this work, we investigate these phenomena in photonic media. Starting with Maxwell's equations with a honeycomb material weight as well as the nonlinear Kerr effect, we first prove the existence of Dirac points in the dispersion surfaces of transverse electric and magnetic Maxwell operators under very general assumptions of the material weight. Our assumptions on the material weight are almost the minimal requirements to ensure the existence of Dirac points in a general hexagonal photonic crystal. We then derive the associated wave packet dynamics in the scenario where the honeycomb structure is weakly modulated. It turns out the reduced envelope equation is generally a two-dimensional nonlinear Dirac equation with a spatially varying mass. By studying the reduced envelope equation with a domain-wall-like mass term, we realize the subtle wave motions which are chiral and immune to local defects. The underlying mechanism is the existence of topologically protected linear line modes, also referred to as edge states. However, we show that these robust linear modes do not survive with nonlinearity. We demonstrate the existence of nonlinear line modes, which can propagate in the nonlinear media based on high-accuracy numerical computations. Moreover, we also report a new type of nonlinear modes which are localized in both directions.

We consider the ensemble of adjacency matrices of Erd{\H o}s-R\'enyi random graphs, i.e.\ graphs on $N$ vertices where every edge is chosen independently and with probability $p \equiv p(N)$. We rescale the matrix so that its bulk eigenvalues are of order one. Under the assumption $p N \gg N^{2/3}$, we prove the universality of eigenvalue distributions both in the bulk and at the edge of the spectrum. More precisely, we prove (1) that the eigenvalue spacing of the Erd{\H o}s-R\'enyi graph in the bulk of the spectrum has the same distribution as that of the Gaussian orthogonal ensemble; and (2) that the second largest eigenvalue of the Erd{\H o}s-R\'enyi graph has the same distribution as the largest eigenvalue of the Gaussian orthogonal ensemble. As an application of our method, we prove the bulk universality of generalized Wigner matrices under the assumption that the matrix entries have at least $4 + \epsilon$ moments.