Implementation and benchmark of a long-range corrected functional in the density functional based tight-binding method

Bridging the gap between first principles methods and empirical schemes, the density functional based tight-binding method (DFTB) has become a versatile tool in predictive atomistic simulations over the past years. One of the major restrictions of this method is the limitation to local or gradient corrected exchange-correlation functionals. This excludes the important class of hybrid or long-range corrected functionals, which are advantageous in thermochemistry, as well as in the computation of vibrational, photoelectron and optical spectra. The present work provides a detailed account of the implementation of DFTB for a long-range corrected functional in generalized Kohn-Sham theory. We apply the method to a set of organic molecules and compare ionization potentials and electron affinities with the original DFTB method and higher level theory. The new scheme cures the significant overpolarization in electric fields found for local DFTB, which parallels the functional dependence in first principles density functional theory (DFT). At the same time the computational savings with respect to full DFT calculations are not compromised as evidenced by numerical benchmark data

Using Density Functional Theory to Model Realistic TiO2 Nanoparticles, Their Photoactivation and Interaction with Water

Computational modeling of titanium dioxide nanoparticles of realistic size is extremely relevant for the direct comparison with experiments but it is also a rather demanding task. We have recently worked on a multistep/scale procedure to obtain global optimized minimum structures for chemically stable spherical titania nanoparticles of increasing size, with diameter from 1.5 nm (~300 atoms) to 4.4 nm (~4000 atoms). We use first self-consistent-charge density functional tight-binding (SCC-DFTB) methodology to perform thermal annealing simulations to obtain globally optimized structures and then hybrid density functional theory (DFT) to refine them and to achieve high accuracy in the description of structural and electronic properties. This allows also to assess SCC-DFTB performance in comparison with DFT(B3LYP) results. As a further step, we investigate photoexcitation and photoemission processes involving electron/hole pair formation, separation, trapping and recombination in the nanosphere of medium size by hybrid DFT. Finally, we show how a recently defined new set of parameters for SCC-DFTB allows for a proper description of titania/water multilayers interface, which paves the way for modeling large realistic nanoparticles in aqueous environment

Structure and electronic structure of Metal-Organic Frameworks within the Density-Functional based Tight-Binding method

Density-functional based tight-binding is a powerful method to describe large molecules and materials. Metal-Organic Frameworks (MOFs), materials with interesting catalytic properties and with very large surface areas have been developed and have become commercially available. Unit cells of MOFs typically include hundreds of atoms, which make the application of standard Density-Functional methods computationally very expensive, sometimes even unfeasible. The aim of this paper is to prepare and to validate the Self-Consistent Charge Density-Functional based Tight Binding (SCC-DFTB) method for MOFs containing Cu, Zn and Al metal centers. The method has been validated against full hybrid density-functional calculations for model clusters, against gradient corrected density-functional calculations for supercells, and against experiment. Moreover, the modular concept of MOF chemistry has been discussed on the basis of their electronic properties. We concentrate on MOFs comprising three common connector units: copper paddlewheels (HKUST-1), zinc oxide Zn4O tetrahedron (MOF-5, MOF-177, DUT-6 (MOF-205)) and aluminium oxide AlO4(OH)2 octahedron (MIL-53). We show that SCC-DFTB predicts structural parameters with a very good accuracy (with less than 5% deviation, even for adsorbed CO and H2O on HKUST-1), while adsorption energies differ by 12 kJ mol-1 or less for CO and water compared to DFT benchmark calculations.Comment: Submitted to Phys. Status Solidi

A comparative study of density functional and density functional tight binding calculations of defects in graphene

The density functional tight binding approach (DFTB) is well adapted for the study of point and line defects in graphene based systems. After briefly reviewing the use of DFTB in this area, we present a comparative study of defect structures, energies and dynamics between DFTB results obtained using the dftb+ code, and density functional results using the localised Gaussian orbital code, AIMPRO. DFTB accurately reproduces structures and energies for a range of point defect structures such as vacancies and Stone-Wales defects in graphene, as well as various unfunctionalised and hydroxylated graphene sheet edges. Migration barriers for the vacancy and Stone-Wales defect formation barriers are accurately reproduced using a nudged elastic band approach. Finally we explore the potential for dynamic defect simulations using DFTB, taking as an example electron irradiation damage in graphene

Mechanical Properties of Phosphorene Nanotubes: A Density Functional Tight-Binding Study

Using density functional tight-binding method, we studied the elastic properties, deformation and failure of armchair (AC) and zigzag (ZZ) phosphorene nano tubes (PNTs) under uniaxial tensile strain. We found that the deformation and failure of PNTs are very anisotropic. For ZZ PNTs, three deformation phases are recognized: The primary linear elastic phase, which is associated with the interactions between the neighboring puckers, succeeded by the bond rotation phase, where the puckered configuration of phosphorene is smoothed via bond rotation, and lastly the bond elongation phase, where the P-P bonds are directly stretched up to the maximally allowed limit and the failure is initiated by the rupture of the most stretched bonds

Quasiparticle energies for large molecules: a tight-binding GW approach

We present a tight-binding based GW approach for the calculation of quasiparticle energy levels in confined systems such as molecules. Key quantities in the GW formalism like the microscopic dielectric function or the screened Coulomb interaction are expressed in a minimal basis of spherically averaged atomic orbitals. All necessary integrals are either precalculated or approximated without resorting to empirical data. The method is validated against first principles results for benzene and anthracene, where good agreement is found for levels close to the frontier orbitals. Further, the size dependence of the quasiparticle gap is studied for conformers of the polyacenes ($C_{4n+2}H_{2n+4}$) up to n = 30.Comment: 10 pages, 5 eps figures submitted to Phys. Rev.

Comparison between Density Functional Theory and Density Functional Tight Binding approaches for finding the muon stopping site in organic molecular crystals

Finding the possible stopping sites for muons inside a crystalline sample is a key problem of muon spectroscopy. In a previous work, we suggested a computational approach to this problem, using Density Functional Theory software in combination with a random structure searching approach using a Poisson sphere distribution. In this work we test this methodology further by applying it to three organic molecular crystals model systems: durene, bithiophene, and tetracyanoquinodimethane (TCNQ). Using the same sets of random structures we compare the performance of Density Functional Theory software CASTEP and the much faster lower level approximation of Density Functional Tight Binding provided by DFTB+, combined with the use of the 3ob-3-1 parameter set. We show the benefits and limitations of such an approach and we propose the use of DFTB+ as a viable alternative to more cumbersome simulations for routine site-finding in organic materials. Finally, we introduce the Muon Spectroscopy Computational Project software suite, a library of Python tools meant to make these methods standardized and easy to use