Ask Hjorth Larsen

Electronic structure calculations with GPAW: A real-space implementation of the projector augmented-wave method

Electronic structure calculations have become an indispensable tool in many areas of materials science and quantum chemistry. Even though the Kohn-Sham formulation of the density-functional theory (DFT) simplifies the many-body problem significantly, one is still confronted with several numerical challenges. In this article we present the projector augmented-wave (PAW) method as implemented in the GPAW program package (https://wiki.fysik.dtu.dk/gpaw) using a uniform real-space grid representation of the electronic wavefunctions. Compared to more traditional plane wave or localized basis set approaches, real-space grids offer several advantages, most notably good computational scalability and systematic convergence properties. However, as a unique feature GPAW also facilitates a localized atomic-orbital basis set in addition to the grid. The efficient atomic basis set is complementary to the more accurate grid, and the possibility to seamlessly switch between the two representations provides great flexibility. While DFT allows one to study ground state properties, time-dependent density-functional theory (TDDFT) provides access to the excited states. We have implemented the two common formulations of TDDFT, namely the linear-response and the time propagation schemes. Electron transport calculations under finite-bias conditions can be performed with GPAW using non-equilibrium Green functions and the localized basis set. In addition to the basic features of the real-space PAW method, we also describe the implementation of selected exchange-correlation functionals, parallelization schemes, ΔSCF-method, x-ray absorption spectra, and maximally localized Wannier orbitals.

Investigation of Catalytic Finite-Size-Effects of Platinum Metal Clusters

In this paper, we use density functional theory (DFT) calculations on highly parallel computing resources to study size-dependent changes in the chemical and electronic properties of platinum (Pt) for a number of fixed freestanding clusters ranging from 13 to 1415 atoms, or 0.7-3.5 nm in diameter. We find that the surface catalytic properties of the clusters converge to the single crystal limit for clusters with as few as 147 atoms (1.6 nm). Recently published results for gold (Au) clusters showed analogous convergence with size. However, this convergence happened at larger sizes, because the Au d-states do not contribute to the density of states around the Fermi-level, and the observed level fluctuations were not significantly damped until the cluster reached ca. 560 atoms (2.7 nm) in size.

Localized atomic basis set in the projector augmented wave method

We present an implementation of localized atomic-orbital basis sets in the projector augmented wave PAW formalism within the density-functional theory. The implementation in the real-space GPAW code provides a complementary basis set to the accurate but computationally more demanding grid representation. The possibility to switch seamlessly between the two representations implies that simulations employing the local basis can be fine tuned at the end of the calculation by switching to the grid, thereby combining the strength of the two representations for optimal performance. The implementation is tested by calculating atomization energies and equilibrium bulk properties of a variety of molecules and solids, comparing to the grid results. Finally, it is demonstrated how a grid-quality structure optimization can be performed with significantly reduced computational effort by switching between the grid and basis representations.

Real-space grids and the Octopus code as tools for the development of new simulation approaches for electronic systems

Real-space grids are a powerful alternative for the simulation of electronic systems. One of the main advantages of the approach is the flexibility and simplicity of working directly in real space where the different fields are discretized on a grid, combined with competitive numerical performance and great potential for parallelization. These properties constitute a great advantage at the time of implementing and testing new physical models. Based on our experience with the Octopus code, in this article we discuss how the real-space approach has allowed for the recent development of new ideas for the simulation of electronic systems. Among these applications are approaches to calculate response properties, modeling of photoemission, optimal control of quantum systems, simulation of plasmonic systems, and the exact solution of the Schrödinger equation for low-dimensionality systems.

Density functional theory based screening of ternary alkali-transition metal borohydrides: A computational material design project

We present a computational screening study of ternary metal borohydrides for reversible hydrogen storage based on density functional theory. We investigate the stability and decomposition of alloys containing 1 alkali metal atom, Li, Na, or K (M(1)); and 1 alkali, alkaline earth or 3d/4d transition metal atom (M(2)) plus two to five (BH(4))(-) groups, i.e., M(1)M(2)(BH(4))(2-5), using a number of model structures with trigonal, tetrahedral, octahedral, and free coordination of the metal borohydride complexes. Of the over 700 investigated structures, about 20 were predicted to form potentially stable alloys with promising decomposition energies. The M(1)(Al/Mn/Fe)(BH(4))(4), (Li/Na)Zn(BH(4))(3), and (Na/K)(Ni/Co)(BH(4))(3) alloys are found to be the most promising, followed by selected M(1)(Nb/Rh)(BH(4))(4) alloys.

Localized surface plasmon resonance in silver nanoparticles: Atomistic first-principles time-dependent density-functional theory calculations

We observe using ab initio methods that localized surface plasmon resonances in icosahedral silver nanoparticles enter the asymptotic region already between diameters of 1 and 2 nm, converging close to the classical quasistatic limit around 3.4 eV. We base the observation on time-dependent density-functional theory simulations of the icosahedral silver clusters Ag-55 (1.06 nm), Ag-147 (1.60 nm), Ag-309 (2.14 nm), and Ag-561 (2.68 nm). The simulation method combines the adiabatic GLLB-SC exchange-correlation functional with real time propagation in an atomic orbital basis set using the projector-augmented wave method. The method has been implemented for the electron structure code GPAW within the scope of this work. We obtain good agreement with experimental data and modeled results, including photoemission and plasmon resonance. Moreover, we can extrapolate the ab initio results to the classical quasistatically modeled icosahedral clusters.

Electronic shell structure and chemisorption on gold nanoparticles

We use density functional theory (DFT) to investigate the electronic structure and chemical properties of gold nanoparticles. Different structural families of clusters are compared. For up to 60 atoms we optimize structures using DFT-based simulated annealing. Cluster geometries are found to distort considerably, creating large band gaps at the Fermi level. For up to 200 atoms we consider structures generated with a simple EMT potential and clusters based on cuboctahedra and icosahedra. All types of cluster geometry exhibit jellium-like electronic shell structure. We calculate adsorption energies of several atoms on the cuboctahedral clusters. Adsorption energies are found to vary abruptly at magic numbers. Using a Newns-Anderson model we find that the effect of magic numbers on adsorption energy can be understood from the location of adsorbate-induced states with respect to the cluster Fermi level.

Modeling electron dynamics coupled to continuum states in finite volumes with absorbing boundaries

Absorbing boundaries are frequently employed in real-time propagation of the Schrödinger equation to remove spurious reflections and efficiently emulate outgoing boundary conditions. These conditions are a fundamental ingredient for the calculation of observables involving infinitely extended continuum states in finite volumes. In the literature, several boundary absorbers have been proposed. They mostly fall into three main families: mask function absorbers, complex absorbing potentials, and exterior complex-scaled potentials. To date none of the proposed absorbers is perfect, and all present a certain degree of reflections. Characterization of such reflections is thus a critical task with strong implications for time-dependent simulations of atoms and molecules. We introduce a method to evaluate the reflection properties of a given absorber and present a comparison of selected samples for each family of absorbers. Further, we discuss the connections between members of each family and show how the same reflection curves can be obtained with very different absorption schemes.

Stark Ionization of Atoms and Molecules within Density Functional Resonance Theory

We show that the energetics and lifetimes of resonances of finite systems under an external electric field can be captured by Kohn–Sham density functional theory (DFT) within the formalism of uniform complex scaling. Properties of resonances are calculated self-consistently in terms of complex densities, potentials, and wave functions using adapted versions of the known algorithms from DFT. We illustrate this new formalism by calculating ionization rates using the complex-scaled local density approximation and exact exchange. We consider a variety of atoms (H, He, Li, and Be) as well as the H2 molecule. Extensions are briefly discussed.

libvdwxc: a library for exchange-correlation functionals in the vdW-DF family

We present libvdwxc, a general library for evaluating the energy and potential for the family of vdW-DF exchange-correlation functionals. libvdwxc is written in C and provides an efficient implementation of the vdW-DF method and can be interfaced with various general-purpose DFT codes. Currently, the GPAW and Octopus codes implement interfaces to libvdwxc. The present implementation emphasizes scalability and parallel performance, and thereby enables ab initio calculations of nanometer-scale complexes. The numerical accuracy is bench-marked on the S22 test set whereas parallel performance is benchmarked on ligand-protected gold nanoparticles (Au-144(SC11NH25)(60)) up to 9696 atoms.