Electrostatic Interactions in Finite Systems treated with Periodic Boundary Conditions: Application to Linear-Scaling Density Functional Theory

Abstract

We present a comparison of methods for treating the electrostatic interactions of finite, isolated systems within periodic boundary conditions (PBCs), within Density Functional Theory (DFT), with particular emphasis on linear-scaling (LS) DFT. Often, PBCs are not physically realistic but are an unavoidable consequence of the choice of basis set and the efficacy of using Fourier transforms to compute the Hartree potential. In such cases the effects of PBCs on the calculations need to be avoided, so that the results obtained represent the open rather than the periodic boundary. The very large systems encountered in LS-DFT make the demands of the supercell approximation for isolated systems more difficult to manage, and we show cases where the open boundary (infinite cell) result cannot be obtained from extrapolation of calculations from periodic cells of increasing size. We discuss, implement and test three very different approaches for overcoming or circumventing the effects of PBCs: truncation of the Coulomb interaction combined with padding of the simulation cell, approaches based on the minimum image convention, and the explicit use of Open Boundary Conditions (OBCs). We have implemented these approaches in the ONETEP LS-DFT program and applied them to a range of systems, including a polar nanorod and a protein. We compare their accuracy, complexity, and rate of convergence with simulation cell size. We demonstrate that corrective approaches within PBCs can achieve the OBC result more efficiently and accurately than pure OBC approaches.

39 pages, 10 figures

References in corpus (7)

• Truncation of Periodic Image Interactions for Confined Systems
• 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
• Accurate ionic forces and geometry optimisation in linear scaling density-functional theory with local orbitals
• Minimal parameter implicit solvent model for ab initio electronic structure calculations
• Linear Algebraic Calculation of Green's function for Large-Scale Electronic Structure Theory

Cited by in corpus (19)

• Electrostatics of solvated systems in periodic boundary conditions
• The role of spin in the calculation of Hubbard $U$ and Hund's $J$ parameters from first principles
• Efficient and accurate solver of the three-dimensional screened and unscreened Poisson's equation with generic boundary conditions
• Linear-scaling time-dependent density-functional theory (TDDFT) beyond the Tamm-Dancoff approximation: obtaining efficiency and accuracy with in situ optimised local orbitals
• A wavelet-based Projector Augmented-Wave (PAW) method: reaching frozen-core all-electron precision with a systematic, adaptive and localized wavelet basis set
• Ligand Discrimination in Myoglobin from Linear-Scaling DFT+U
• Supercell convergence of charge-transfer energies in pentacene molecular crystals from constrained DFT
• Fermi level pinning can determine polarity in semiconductor nanorods
• Charged surfaces and slabs in periodic boundary conditions
• Band alignment and interlayer hybridisation in transition metal dichalcogenide/hexagonal boron nitride heterostructures
• Dopant levels in large nanocrystals using stochastic optimally tuned range-separated hybrid density functional theory
• Tracing potential energy surfaces of electronic excitations via their transition origins: application to Oxirane
• Simulations of Nanocrystals Under Pressure: Combining Electronic Enthalpy and Linear-Scaling Density-Functional Theory
• Bromophenyl functionalization of carbon nanotubes : an ab initio study
• Wannier-function-based constrained DFT with nonorthogonality-correcting Pulay forces in application to the reorganization effects in graphene-adsorbed pentacene
• Strain-dependent one-dimensional confinement channels in twisted bilayer 1T${}^{\prime }$-WTe${}_2$
• Electrostatic treatment of charged interfaces in classical atomistic simulations
• Density functional calculations of atomic structure, charging effect, and static dielectric constant of two-dimensional systems based on B-splines
• Molecules with dipoles in periodic boundary conditions in a tetragonal cell