Perspective: Methods for large-scale density functional calculations on metallic systems

Abstract

Perspective: Explicitly correlated electronic structure theory for complex systems The Journal of Chemical Physics 146, 080901 (2017) Current research challenges in areas such as energy and bioscience have created a strong need for Density Functional Theory (DFT) calculations on metallic nanostructures of hundreds to thousands of atoms to provide understanding at the atomic level in technologically important processes such as catalysis and magnetic materials. Linear-scaling DFT methods for calculations with thousands of atoms on insulators are now reaching a level of maturity. However such methods are not applicable to metals, where the continuum of states through the chemical potential and their partial occupancies provide significant hurdles which have yet to be fully overcome. Within this perspective we outline the theory of DFT calculations on metallic systems with a focus on methods for large-scale calculations, as required for the study of metallic nanoparticles. We present early approaches for electronic energy minimization in metallic systems as well as approaches which can impose partial state occupancies from a thermal distribution without access to the electronic Hamiltonian eigenvalues, such as the classes of Fermi operator expansions and integral expansions. We then focus on the significant progress which has been made in the last decade with developments which promise to better tackle the lengthscale problem in metals. We discuss the challenges presented by each method, the likely future directions that could be followed and whether an accurate linear-scaling DFT method for metals is in sight. Published by AIP Publishing. [http://dx

On nearsightedness in metallic systems for O(N) Density Functional Theory calculations: A case study on aluminum

“…The truncation that can be employed in practice, while maintaining the desired accuracy, determines the prefactor of O(N) methods as well as the system sizes at which linear scaling can be achieved in practical calculations. While this has been studied previously for insulators [22,23], the dependence of truncation region size on smearing for a given accuracy in metallic systems has not been carefully studied heretofore [24].…”

Predicting the Oxygen-Binding Properties of Platinum Nanoparticle Ensembles by Combining High-Precision Electron Microscopy and Density Functional Theory

“…The major challenge in developing linear-scaling DFT approaches is being able to retain the same high level of accuracy as conventional cubic-scaling DFT. DFT calculations for metallic systems are more challenging than those for insulators and the development of reduced-scaling, or even linearscaling methods, for metallic systems is an open area of research 25 . While a linear-scaling DFT method for metals is not yet available, there has been significant progress in the development of reduced-scaling DFT methods for large-scale calculations on metallic systems.…”

SQDFT: Spectral Quadrature method for large-scale parallel O(N) Kohn–Sham calculations at high temperature

“…In order to update the positions of the ions during the Born-Oppenheimer quantum molecular dynamics (QMD) simulations, the Hellmann-Feynman force on the I th nucleus can be written as [38] 4 The electronic temperature/smearing is typically set to be equal to the ionic temperature in QMD simulations, particularly for those performed at high temperature. 5 The repulsive energy correction for overlapping pseudocharges [45,46] has been explicitly included in the free energy expression. This term plays a particularly important role in high-temperature simulations since the ions get significantly closer compared to ambient temperature.…”