We study a robust and efficient eigensolver for computing the positive dense spectrum of the two-dimensional transmission eigenvalue problem (TEP) which is derived from the Maxwell’s equation with complex media in pseudo-chiral model and the transverse magnetic mode. The discretized governing equations by the N ́ed ́elec edge element result in a large-scale quadratic eigenvalue problem (QEP). We estimate half of the positive eigenvalues of the QEP are on some interval which forms a dense spectrum of the QEP. The quadratic Jacobi-Davidson method with a so-called non-equivalence deflation technique is proposed to compute the dense spectrum of the QEP. Intensive numerical experiments show that our proposed method makes the convergence efficiently and robustly even it needs to compute more than 5000 desired eigenpairs. Numerical results also illustrate that the computed eigenvalue curves can be approximated by the non- linear functions which can be applied to estimate the density of the eigenvalues for the TEP.

We consider Hermitian and symmetric random band matrices $H = (h_{xy})$ in $d \geq 1$ dimensions. The matrix entries $h_{xy}$, indexed by $x,y \in (\bZ/L\bZ)^d$, are independent, centred random variables with variances $s_{xy} = \E |h_{xy}|^2$. We assume that $s_{xy}$ is negligible if $|x-y|$ exceeds the band width $W$. In one dimension we prove that the eigenvectors of $H$ are delocalized if $W\gg L^{4/5}$. We also show that the magnitude of the matrix entries $\abs{G_{xy}}^2$ of the resolvent $G=G(z)=(H-z)^{-1}$ is self-averaging and we compute $\E \abs{G_{xy}}^2$. We show that, as $L\to\infty$ and $W\gg L^{4/5}$, the behaviour of $\E |G_{xy}|^2$ is governed by a diffusion operator whose diffusion constant we compute. Similar results are obtained in higher dimensions.

No abstract uploaded!

In this paper, we apply the local discontinuous Galerkin (LDG) method to 2D Keller--Segel (KS) chemotaxis model. We improve the results upon (Y. Epshteyn and A. Kurganov, SIAM Journal on Numerical Analysis, 47 (2008), 368-408) and give optimal rate of convergence under special finite element spaces before the blow-up occurs (the exact solutions are smooth). Moreover, to construct physically relevant numerical approximations, we consider $P^1$ LDG scheme and develop a positivity-preserving limiter to the scheme, extending the idea in (Y. Zhang, X. Zhang and C.-W. Shu, Journal of Computational Physics, 229 (2010), 8918-8934). With this limiter, we can prove the $L^1$-stability of the numerical scheme. Numerical experiments are performed to demonstrate the good performance of the positivity-preserving LDG scheme. Moreover, it is known that the chemotaxis model will yield blow-up solutions under certain initial conditions. We numerically demonstrate how to find the approximate blow-up time by using the $L^2$-norm of the $L^1$-stable numerical solution.

The aim of this paper is twofold. First we give an explicit construction of the in¯nitesimal deformations of the category Coh(X) of coherent sheaves on a smooth projective variety X. Secondly,we show that any Fourier-Mukai transform ©: Db(X) ! Db(Y ) extends to an equivalence between the derived categories of the deformed Abelian categories.