Research

Convergence analysis of planewave expansion methods for 2D Schrodinger operators with discontinuous periodic potentials


Reference:

Norton, R. and Scheichl, R., 2010. Convergence analysis of planewave expansion methods for 2D Schrodinger operators with discontinuous periodic potentials. SIAM Journal on Numerical Analysis (SINUM), 47 (6), pp. 4356-4380.

Related documents:

[img]
Preview
PDF (SIAMJNumAnalysis47(6)p4356.pdf) - Requires a PDF viewer such as GSview, Xpdf or Adobe Acrobat Reader
Download (446kB) | Preview

    Official URL:

    http://dx.doi.org/10.1137/090756545

    Abstract

    In this paper we consider the problem of computing the spectrum of a two-dimensional Schrodinger operator with discontinuous, periodic potential in two dimensions using Fourier (or planewave expansion) methods. Problems of this kind are currently of great interest in the design of new optical devices to determine band gaps and to compute localized modes in photonic crystal materials. Although Fourier methods may not be every applied mathematician's first choice for this problem because of the discontinuities in the potential, we will show here that, even though (as expected) the convergence is not exponential, the method has several desirable features that make it competitive with other discretization techniques, such as finite element methods, both with respect to implementation and convergence properties. In particular, we will prove that simple preconditioners for the system matrix are optimal leading to a computational complexity of O(N logN) in the number of planewaves N (using the fast Fourier transform). Moreover, we derive sharp error estimates that show that the method is essentially third order in the eigenvalues and of order 3 2 in the eigenfunctions in the H-1-norm and 5/2 in the L-2-norm. To improve the planewave expansion method in the case of discontinuous potentials, it has been proposed in the physics literature to replace the discontinuous potential with an effective potential that is smooth, despite the additional error this incurs. We will here answer the question whether this smoothing is worth it. In fact, our convergence analysis of the modified method provides an optimal choice for the smoothing parameter, but it also shows that the overall rate of convergence is no faster than before and so smoothing does not seem to be worth it. All the theoretical results are confirmed in our numerical experiments.

    Details

    Item Type Articles
    CreatorsNorton, R.and Scheichl, R.
    DOI10.1137/090756545
    Uncontrolled Keywordsdiscontinuous coefficients, smoothing, pde eigenproblems, error analysis, spectral approximation, fourier methods, optimal preconditioners
    DepartmentsFaculty of Science > Mathematical Sciences
    Publisher StatementSIAMJNumAnalysis47(6)p4356.pdf: © SIAM
    RefereedYes
    StatusPublished
    ID Code19094

    Export

    Actions (login required)

    View Item

    Document Downloads

    More statistics for this item...