Enabling Ab Initio Molecular Dynamics under Bias: The CP2K+SMEAGOL Interface for Integrating Density Functional Theory and Non-Equilibrium Green Functions.

Density functional theory (DFT) combined with non-equilibrium Green's functions (NEGF) is a powerful approach to model quantum transport under external bias potentials at reasonable computational cost. In this work, we present a new interface between the popular mixed Gaussian/plane waves electronic structure package, CP2K, and the NEGF, code SMEAGOL, the most feature-rich implementation of DFT-NEGF available for CP2K to date. The CP2K+SMEAGOL interface includes the implementation of current induced forces. We verify this implementation for a variety of systems: an infinite 1D Au wire, a parallel-plate capacitor, and a Au-H2-Au junction. We find good agreement with SMEAGOL calculations performed with SIESTA for the same systems and with the example of a solvated Au wire demonstrating for the first time that DFT-NEGF can be used to perform molecular dynamics simulations under bias of large-scale condensed phase systems under realistic operating conditions.

Excited-State Properties for Extended Systems: Efficient Hybrid Density Functional Methods.

Time-dependent density functional theory has become state-of-the-art for describing photophysical and photochemical processes in extended materials because of its affordable cost. The inclusion of exact exchange was shown to be essential for the correct description of the long-range asymptotics of electronic interactions and thus a well-balanced description of valence, Rydberg, and charge-transfer excitations. Several approaches for an efficient treatment of exact exchange have been established for the ground state, while implementations for excited-state properties are rare. Furthermore, the high computational costs required for excited-state properties in comparison to ground-state computations often hinder large-scale applications on periodic systems with hybrid functional accuracy. We therefore propose two approximate schemes for improving computational efficiency for the treatment of exact exchange. Within the auxiliary density matrix method (ADMM), exact exchange is estimated using a relatively small auxiliary basis and the introduced basis set incompleteness error is compensated by an exchange density functional correction term. Benchmark results for a test set of 35 molecules demonstrate that the mean absolute error introduced by ADMM is smaller than 0.3 pm for excited-state bond lengths and in the range of 0.02-0.04 eV for vertical excitation, adiabatic excitation, and fluorescence energies. Computational timings for a series of covalent-organic frameworks demonstrate that a speed-up of at least 1 order of magnitude can be achieved for excited-state geometry optimizations in comparison to conventional hybrid functionals. The second method is to use a semiempirical tight binding approximation for both Coulomb and exchange contributions to the excited-state kernel. This simplified Tamm-Dancoff approximation (sTDA) achieves an accuracy comparable to approximated hybrid density functional theory when referring to highly accurate coupled-cluster reference data. We find that excited-state bond lengths deviate by 1.1 pm on average and mean absolute errors in vertical excitation, adiabatic excitation, and fluorescence energies are in the range of 0.2-0.5 eV. In comparison to ADMM-approximated hybrid functional theory, sTDA accelerates the computation of broad-band excitation spectra by 1 order of magnitude, suggesting its potential use for large-scale screening purposes.

The Role of Cation-Vacancies for the Electronic and Optical Properties of Aluminosilicate Imogolite Nanotubes: A Non-local, Linear-Response TDDFT Study.

We report a combined non-local (PBE-TC-LRC) Density Functional Theory (DFT) and linear-response time-dependent DFT (LR-TDDFT) study of the structural, electronic, and optical properties of the cation-vacancy based defects in aluminosilicate (AlSi) imogolite nanotubes (Imo-NTs) that have been recently proposed on the basis of Nuclear Magnetic Resonance (NMR) experiments. Following numerical determination of the smallest AlSi Imo-NT model capable of accommodating the defect-induced relaxation with negligible finite-size errors, we analyse the defect-induced structural deformations in the NTs and ensuing changes in the NTs' electronic structure. The NMR-derived defects are found to introduce both shallow and deep occupied states in the pristine NTs' band gap (BG). These BG states are found to be highly localized at the defect site. No empty defect-state is modeled for any of the considered systems. LR-TDDFT simulation of the defects reveal increased low-energy optical absorbance for all but one defects, with the appearance of optically active excitations at energies lower than for the defect-free NT. These results enable interpretation of the low-energy tail in the experimental UV-vis spectra for AlSi NTs as being due to the defects. Finally, the PBE-TC-LRC-approximated exciton binding energy for the defects' optical transitions is found to be substantially lower (up to 0.8 eV) than for the pristine defect-free NT's excitations (1.1 eV).

First principles calculations of optical properties for oxygen vacancies in binary metal oxides.

Using an advanced computational methodology implemented in CP2K, a non-local PBE0-TC-LRC density functional and the recently implemented linear response formulation of the Time-dependent Density Functional Theory equations, we test the interpretation of the optical absorption and photoluminescence signatures attributed by previous experimental and theoretical studies to O-vacancies in two widely used oxides-cubic MgO and monoclinic (m)-HfO2. The results obtained in large periodic cells including up to 1000 atoms emphasize the importance of accurate predictions of defect-induced lattice distortions. They confirm that optical transitions of O-vacancies in 0, +1, and +2 charge states in MgO all have energies close to 5 eV. We test the models of photoluminescence of O-vacancies proposed in the literature. The photoluminescence of VO +2 centers in m-HfO2 is predicted to peak at 3.7 eV and originate from radiative tunneling transition between a VO +1 center and a self-trapped hole created by the 5.2 eV excitation.

A fragment method for systematic improvement of anharmonic adsorbate vibrational frequencies: acetylene on Cu(001).

We suggest a novel method for systematic improvement of anharmonic adsorbate frequencies based on a fragment approach. The calculations are carried out by considering the adsorbed molecule separately and computing an energy correction using high-level ab initio method in addition to a standard calculation of the whole adsorbed system using quantum mechanical techniques with periodic boundary conditions. We demonstrate its reliability for a C2H2 molecule chemisorbed on a Cu(001) surface. We also show that the accuracy of the presented approach with a suitable description of the periodic surface depends mainly on the accuracy of the high-level ab initio method used to describe the adsorbate molecule. Moreover, our technique potentially allows to predict adsorbate vibrational spectra with spectroscopic accuracy.

Reconstruction of Density Functionals from Kohn-Sham Potentials by Integration along Density Scaling Paths.

We demonstrate by specific examples that if a Kohn-Sham exchange-correlation potential is given explicitly in terms of the electron density and its derivatives, then one can easily reconstruct the parent density functional by evaluating analytically (or numerically with one-dimensional quadratures) the van Leeuwen-Baerends line integral (Phys. Rev. A 1995, 51, 170-178) along a path of (coordinate)-scaled densities. The choice of a density scaling path amounts to defining the gauge of the resultant exchange-correlation energy density. The well-known Levy-Perdew virial relation for exchange potentials can be viewed as an analytical line integral along the electron-number-conserving uniform density scaling path. Energies obtained from model exchange-correlation potentials should be interpreted with caution because the reconstructed density functional is unique (up to a gauge transformation) only if the model Kohn-Sham potential is a functional derivative.