Bálint Aradi

An Improved Self-Consistent-Charge Density-Functional Tight-Binding (SCC-DFTB) Set of Parameters for Simulation of Bulk and Molecular Systems Involving Titanium

A new self-consistent-charge density-functional tight-binding (SCC-DFTB) set of parameters for Ti-X pairs of elements (X = Ti, H, C, N, O, S) has been developed. The performance of this set has been tested with respect to TiO2 bulk phases and small molecular systems. It has been found that the band structures, geometric parameters, and cohesive energies of rutile and anatase polymorphs are in good agreement with the reference DFT data and with experiment. Low-index rutile and anatase surfaces were also tested. For molecular systems, binding and atomization energies close to their DFT analogues have been achieved. Large errors, however, have been found for systems in high-spin states and/or having multireference character of their wave functions. The correct performance of SCC-DFTB for surface reactions has been demonstrated via the water splitting on anatase (001) surface. The current SCC-DFTB set is a suitable tool for future in-depth investigation of chemical processes occurring on the surfaces of TiO2 polymorphs as well as for other processes of physicochemical interest.

Parametrization of the SCC-DFTB Method for Halogens

Parametrization of the approximative DFT method SCC-DFTB for halogen elements is presented. The new parameter set is intended to describe halogenated organic as well as inorganic molecules, and it is compatible with the established parametrization of SCC-DFTB for carbon, hydrogen, oxygen, and nitrogen. The performance of the parameter set is tested on a representative set of molecules and discussed.

Extensions of the Time-Dependent Density Functional Based Tight-Binding Approach

The time-dependent density functional based tight-binding (TD-DFTB) approach is generalized to account for fractional occupations. In addition, an on-site correction leads to marked qualitative and quantitative improvements over the original method. Especially, the known failure of TD-DFTB for the description of σ → π* and n → π* excitations is overcome. Benchmark calculations on a large set of organic molecules also indicate a better description of triplet states. The accuracy of the revised TD-DFTB method is found to be similar to first principles TD-DFT calculations at a highly reduced computational cost. As a side issue, we also discuss the generalization of the TD-DFTB method to spin-polarized systems. In contrast to an earlier study, we obtain a formalism that is fully consistent with the use of local exchange-correlation functionals in the ground state DFTB method.

Nonadiabatic Molecular Dynamics for Thousand Atom Systems: A Tight-Binding Approach toward PYXAID

Excited state dynamics at the nanoscale requires treatment of systems involving hundreds and thousands of atoms. In the majority of cases, depending on the process under investigation, the electronic structure component of the calculation constitutes the computation bottleneck. We developed an efficient approach for simulating nonadiabatic molecular dynamics (NA-MD) of large systems in the framework of the self-consistent charge density functional tight binding (SCC-DFTB) method. SCC-DFTB is combined with the fewest switches surface hopping (FSSH) and decoherence induced surface hopping (DISH) techniques for NA-MD. The approach is implemented within the Python extension for the ab initio dynamics (PYXAID) simulation package, which is an open source NA-MD program designed to handle nanoscale materials. The accuracy of the developed approach is tested with ab initio DFT and experimental data, by considering intraband electron and hole relaxation, and nonradiative electron-hole recombination in a CdSe quantum dot and the (10,5) semiconducting carbon nanotube. The technique is capable of treating accurately and efficiently excitation dynamics in large, realistic nanoscale materials, employing modest computational resources.

Proper Surface Termination for Luminescent Near-Surface NV Centers in Diamond

By accurate quantum mechanical simulations, we show that typical diamond surfaces possess image states with sub-bandgap energies, and compromise the photostability of NV centers placed within a few nm of the surface. This occurs due to the mixture of the NV-related gap states and the surface image states, which is a novel and distinct process from the well-established band bending effect. We also find that certain types of coverages on the diamond surface may lead to blinking or bleaching due to the presence of acceptor surface states. We identify a combination of surface terminators that is perfect for NV-center based nanoscale sensing.

Comparison of Nb- and Ta-doping of anatase TiO2 for transparent conductor applications

Nb- or Ta-doped anatase TiO2 was shown to be a viable candidate for replacing indium-tin-oxide as a transparent conductive oxide. Calculating the electronic structures we find that Ta has the considerably higher solubility and lower optical effective mass of the two dopants. Our calculations also show that a reducing atmosphere is necessary for efficient dopant incorporation, and oxygen vacancies do not necessarily play a role in their activation.

Prediction of energetically optimal single-walled carbon nanotubes for hydrogen physisorption

Hydrogen storage by carbon nanotubes (CNTs) is a challenging issue still in debate. Using an approximate density functional method augmented with a van der Waals dispersion term, we have shown that there are binding maxima for H2/single-walled carbon nanotube (SWCNT) complexes at (5, 5) and (8, 0) tubes for armchair and zigzag CNTs, respectively, with binding energies around three times as large as that of H2 on graphene surface. We predict that SWCNTs with diameters of 6–7 A are energetically optimal candidates for physisorption of molecular hydrogen.

SCC-DFTB parameters for simulating hybrid gold-thiolates compounds

We present a parametrization of a self‐consistent charge density functional‐based tight‐binding scheme (SCC‐DFTB) to describe gold‐organic hybrid systems by adding new Au‐X (X = Au, H, C, S, N, O) parameters to a previous set designed for organic molecules. With the aim of describing gold‐thiolates systems within the DFTB framework, the resulting parameters are successively compared with density functional theory (DFT) data for the description of Au bulk, Aun gold clusters (n = 2, 4, 8, 20), and AunSCH3 (n = 3 and 25) molecular‐sized models. The geometrical, energetic, and electronic parameters obtained at the SCC‐DFTB level for the small Au3SCH3 gold–thiolate compound compare very well with DFT results, and prove that the different binding situations of the sulfur atom on gold are correctly described with the current parameters. For a larger gold–thiolate model, Au25SCH3, the electronic density of states and the potential energy surfaces resulting from the chemisorption of the molecule on the gold aggregate obtained with the new SCC‐DFTB parameters are also in good agreement with DFT results.

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.

Graphene nucleation on a surface-molten copper catalyst: quantum chemical molecular dynamics simulations

Chemical vapor deposition (CVD) growth of graphene on Cu(111) has been modeled with quantum chemical molecular dynamics (QM/MD) simulations. These simulations demonstrate at the atomic level how graphene forms on copper surfaces. In contrast to other popular catalysts, such as nickel and iron, copper is in a surface molten state throughout graphene growth at CVD-relevant temperatures, and graphene growth takes place without subsurface diffusion of carbon. Surface Cu atoms have remarkably high mobilities on the Cu(111) surface, both before and after graphene nucleation. This surface mobility drives “defect healing” processes in the nucleating graphene structure that convert defects such as pentagons and heptagons into carbon hexagons. Consequently, the graphene defects that become “kinetically trapped” using other catalysts, such as Ni and Fe, are less commonly observed in the case of Cu. We propose this mechanism to be the basis of coppers ability to form high-quality, large-domain graphene flakes.