Fast solution of elliptic partial differential equations using linear combinations of plane waves

Abstract

Given an arbitrary elliptic partial differential equation (PDE), a procedure for obtaining its solution is proposed based on the method of Ritz: the solution is written as a linear combination of plane waves and the coefficients are obtained by variational minimization. The PDE to be solved is cast as a system of linear equations Ax = b, where the matrix A is not sparse, which prevents the straightforward application of standard iterative methods in order to solve it. This sparseness problem can be circumvented by means of a recursive bisection approach based on the fast Fourier transform, which makes it possible to implement fast versions of some stationary iterative methods (such as Gauss-Seidel) consuming O(N log N ) memory and executing an iteration in O(N log 2 N ) time, N being the number of plane waves used. In a similar way, fast versions of Krylov subspace methods and multigrid methods can also be implemented. These procedures are tested on Poisson's equation expressed in adaptive coordinates. It is found that the best results are obtained with the GMRES method using a multigrid preconditioner with Gauss-Seidel relaxation steps.

Spin-orbit coupling effects and g-factors in zinc-blende InSb and wurtzite InAs nanowires using realistic multiband k · p method

“…These are the same axes as used in the confinement geometry a).Source: By the author. expansion,126, 172,[177][178][179][180] i. e., a Fourier transformation of the space dependent parameters.To simulate the vacuum region, we replaced the bulk parameters by a set of degenerate free electron bands with band mismatch of about 5.5 eV for conduction band and −2.5 eV for valence band, seeFig. 14.…”

“…capable of solving this type of equation such as finite differences,174,175 finite elements176 and spectral methods126, 172,[177][178][179][180] to cite a few. The solution of equations like Eq.…”