An accurate and reliable real space method for the ab initio calculation of electronic-structures of materials has been desired. Historically, the most popular method in this field has been the Plane Wave method. However, because the basis functions of the Plane Wave method are not local in real space, it is inefficient to represent the highly localized inner-shell electron state and it generally give rise to a large dense potential matrix which is difficult to deal with. Moreover, it is not suitable for parallel computers, because it requires Fourier transformations. These limitations of the Plane Wave method have led to the development of various real space methods including finite difference method and finite element method, and studies are still in progress. Recently, we have proposed a new numerical method, the CIP-Basis Set (CIP-BS) method [1], by generalizing the concept of the Constrained Interpolation Profile (CIP) method from the viewpoint of the basis set. This method uses a simple polynomial basis set that is easily extendable to any desired higher-order accuracy. The interpolating profile is chosen so that the sub-grid scale solution approaches the local real solution by the constraints from the spatial derivative of the original equation. Thus the solution even on the sub-grid scale becomes consistent with the master equation. By increasing the order of the polynomial, this solution quickly converges. The governing equations are unambiguously discretized into matrix form equations requiring the residuals to be orthogonal to the basis functions via the same procedure as the Galerkin method. We have already demonstrated that the method can be applied to calculations of the band structures for crystals with pseudopotentials. It has been certified that the method gives accurate solutions in the very coarse meshes and the errors converge rapidly when meshes are refined. Although, we have dealt with problems in which potentials are represented analytically, in Kohn-Sham equation the potential is obtained by solving Poisson equation, where the charge density is determined by using wave functions. In this paper, we present the CIP-BS method gives accurate solutions for Poisson equation. Therefore, we believe that the method would be a promising method for solving self-consistent eigenvalue problems in real space.
Skip Nav Destination
ASME 2009 InterPACK Conference collocated with the ASME 2009 Summer Heat Transfer Conference and the ASME 2009 3rd International Conference on Energy Sustainability
July 19–23, 2009
San Francisco, California, USA
Conference Sponsors:
- Electronic and Photonic Packaging Division
ISBN:
978-0-7918-4360-4
PROCEEDINGS PAPER
Numerical Solutions of Poisson Equation by the CIP-Basis Set Method Available to Purchase
Daisuke Shibata,
Daisuke Shibata
Tokyo University of Science Yamaguchi, Sanyo-Onoda, Japan
Search for other works by this author on:
Takayuki Utsumi
Takayuki Utsumi
Tokyo University of Science Yamaguchi, Sanyo-Onoda, Japan
Search for other works by this author on:
Daisuke Shibata
Tokyo University of Science Yamaguchi, Sanyo-Onoda, Japan
Takayuki Utsumi
Tokyo University of Science Yamaguchi, Sanyo-Onoda, Japan
Paper No:
InterPACK2009-89150, pp. 91-98; 8 pages
Published Online:
December 24, 2010
Citation
Shibata, D, & Utsumi, T. "Numerical Solutions of Poisson Equation by the CIP-Basis Set Method." Proceedings of the ASME 2009 InterPACK Conference collocated with the ASME 2009 Summer Heat Transfer Conference and the ASME 2009 3rd International Conference on Energy Sustainability. ASME 2009 InterPACK Conference, Volume 2. San Francisco, California, USA. July 19–23, 2009. pp. 91-98. ASME. https://doi.org/10.1115/InterPACK2009-89150
Download citation file:
4
Views
Related Proceedings Papers
Related Articles
An Improved Three-Dimensional Level Set Method for Gas-Liquid Two-Phase Flows
J. Fluids Eng (July,2004)
Ritz Legendre Multiwavelet Method for the Damped Generalized Regularized Long-Wave Equation
J. Comput. Nonlinear Dynam (January,2012)
Related Chapters
Conclusions
Introduction to Finite Element, Boundary Element, and Meshless Methods: With Applications to Heat Transfer and Fluid Flow
Approximate Analysis of Plates
Design of Plate and Shell Structures
On Interfacial Ribbons of ζ- and γ-Hydride Phases Surrounding δ Precipitates in Zircaloy-2
Zirconium in the Nuclear Industry: 20th International Symposium