Efficient and accurate solver of the three-dimensional screened and unscreened Poisson's equation with generic boundary conditions

Abstract

We present an explicit solver of the three-dimensional screened and unscreened Poisson's equation which combines accuracy, computational efficiency and versatility. The solver, based on a mixed plane-wave / interpolating scaling function representation, can deal with any kind of periodicity (along one, two, or three spatial axes) as well as with fully isolated boundary conditions. It can seamlessly accommodate a finite screening length, non-orthorhombic lattices and charged systems. This approach is particularly advantageous because convergence is attained by simply refining the real space grid, namely without any adjustable parameter. At the same time, the numerical method features O(N log N) scaling of the computational cost (N being the number of grid points) very much like plane-wave methods. The methodology, validated on model systems, is tailored for leading-edge computer simulations of materials (including ab initio electronic structure computations), but it might as well be beneficial for other research domains.

10 Pages, 14 figures. Submitted to The Journal of Chemical Physics

References in corpus (8)

• An exact Coulomb cutoff technique for supercell calculations
• Daubechies wavelets as a basis set for density functional pseudopotential calculations
• Efficient solution of Poisson's equation with free boundary conditions
• Electrostatics in Periodic Boundary Conditions and Real-space Corrections
• Efficient and accurate three dimensional Poisson solver for surface problems
• Supercell technique for total-energy calculations of finite charged and polar systems
• Electrostatic Interactions in Finite Systems treated with Periodic Boundary Conditions: Application to Linear-Scaling Density Functional Theory
• Wavelet-Based Linear-Response Time-Dependent Density-Functional Theory

Cited by in corpus (13)

• Challenges in Large Scale Quantum Mechanical Calculations
• A generalized Poisson and Poisson-Boltzmann solver for electrostatic environments
• Real-space density functional theory on graphical processing units: computational approach and comparison to Gaussian basis set methods
• Common workflows for computing material properties using different quantum engines
• A wavelet-based Projector Augmented-Wave (PAW) method: reaching frozen-core all-electron precision with a systematic, adaptive and localized wavelet basis set
• The CECAM Electronic Structure Library and the modular software development paradigm
• Accurate and efficient computation of nonlocal potentials based on Gaussian-sum approximation
• Fragment Approach to Constrained Density Functional Theory Calculations using Daubechies Wavelets
• A Customized 3D GPU Poisson Solver for Free Boundary Conditions
• Boundary conditions for the solution of the 3-dimensional Poisson equation in open metallic enclosures
• A kernel-independent sum-of-Gaussians method by de la Vallée-Poussin sums
• Fast calculation of retarded potentials in multi-domain TDDFT
• An integral equation method for the advection-diffusion equation on time-dependent domains in the plane