Tight-binding Density Functional Theory Research Articles

A tight-binding study of the electron transport through single-walled carbon nanotube-graphene hybrid nanostructures.

Hybrid carbon nanostructures based on the sp2 hybridized allotropes of carbon, such as graphene and single-walled carbon nanotubes (SWCNTs), hold vast potential for applications in electronics of various forms. Electronic properties of such hybrid structures are modified due to the interaction between atoms of the components, which can be utilized to tailor the properties of the hybrid structures to suite the application. In this study, we have explored charge (electron) transport through the hybrid structures of single-layer graphene (SLG) and SWCNTs (both metallic and semiconducting) using the nonequilibrium Green's function formalism within the framework of tight-binding density functional theory. Our calculations show that the electronic transport in hybrid nanostructures is affected by the interactions between SWCNT and SLG in comparison to the individual components. The changes in the electronic structure and the transport properties with increasing interaction in hybrids (captured by decreasing the separation between SWCNT and SLG) are discussed, and it is demonstrated from this analysis that the hybrids with semiconducting SWCNTs and metallic SWCNTs show different behavior in the low bias regime while they show similar behavior at higher biases. The difference in the transport properties of hybrids with semiconducting and metallic SWCNTs is explained in terms of changes in the electronic structure, the local density of states, and the energy dispersion for electrons due to the interaction between atoms of the two components.

Overlap energy role investigation on location of π and π* electronic bands in α-graphyne

Graphyne with outstanding properties is a promising material in technological applications such as sensors and optoelectronic devices. In this study, relaxation effect in the form of overlap energy variation on π and π* bands location of α-graphyne is investigated. Overlap energy can be controlled by external parameters such as stress and applied electric field. The Eigen values of band energy for eight possible cases based on the tight-binding model is calculated and the intertwined effect of single and double bonds energies as an overlap energy parameter is discussed. Moreover, a comparative study is done through tight-binding model and density functional theory method, which supports the accuracy of the presented model. The remarkable properties of α-graphyne mentioned in this research and obtained results bring new hopes for using of this material in nanoelectronics and optoelectronics applications.

Combined experimental and theoretical studies on enlarged bandgap and improved photoelectrochemical properties of reduced graphene oxide film by hydrogen annealing

• Relationship between thermal annealing, sp 2 domain, bandgap and properties in RGO. • Correlation between experimental data with calculations (GFN2-xTB, sTD-DFT). • The significant increase of 45.8% in the PEC hydrogen evolution for RGO1100 films. The effect of thermal annealing on the recovery of sp 2 carbon networks in the reduced graphene oxide (RGO) and its consequent bandgap enlargement and electrical conductivity improvement were thoroughly studied, experimentally and theoretically correlated. The annealed RGO film at 1100 °C (designated as RGO1100) exhibited a low sheet resistance (1250 Ω.cm −2 ) with a wide bandgap of 2.82 eV. The photoelectrochemical (PEC) experiments in 1 M H 2 SO 4 solution under solar simulation light (1.5 AM) showed photocurrent densities of 48 and 70 μA.cm −2 for non-annealed RGO (72.5% sp 2 fraction) and RGO1100 (87% sp 2 fraction) respectively, equivalent to an increase of 45.8% for RGO1100. The bandgap widening as well as electrical conductivity increase of RGO1100 were closely interpreted on the basis of sp 2 domain restoration, electronic and molecular orbital structures obtained from extensive experimental data as well as elaborated Geometry, Frequency, Noncovalent Extended Tight-binding (GFN2-xTB) and simplified Tamm-Dancoff density functional theory (sTD-DFT) calculations. The correlation between the electronic structure and the PEC performance of the RGO was also evaluated when the impact of the spatial distribution of oxidation and reduction regions on the recombination rate between photogenerated electrons and holes had been taken into account. The findings shed light on the structure-PEC relationship of RGO1100, which will definitely assist in designing and tailoring RGO-based, highly active materials for PEC applications.

Atom Equivalent Energies for the Rapid Estimation of the Heat of Formation of Explosive Molecules from Density Functional Tight Binding Theory.

Atom equivalent energies have been derived from which the gas-phase heat of formation of explosive molecules can be estimated from fast, semiempirical density functional tight binding total energy calculations. The root-mean-square deviation and maximum deviation of the heats of formation from the experimental values for the set of 45 energetic molecules compiled by Byrd and Rice [ J. Phys. Chem. A, 2006, 110, 1005-1013] are 10.4 and 25.5 kcal/mol, respectively, using 4 atom equivalent energies and 7.4 and 15.0 kcal/mol, respectively, using 7 atom equivalent energies. These errors are around a factor of 2-3 larger than those obtained from density functional theory calculations but are smaller than those obtained from other semiempirical electronic structure methods. Heats of formation calculated with density functional tight binding theory using the 4 and 7 atom equivalent energies, the Byrd and Rice scheme, and the atom pair contribution method for a new set of 531 energetic molecules that contain only carbon, hydrogen, nitrogen, and oxygen are provided.

Theoretical study of complex susceptibility of Pd and Pt

Spin transfer torque in a multilayer magnetic system can be understood in terms of interaction between a normal metal and a ferromagnetic layer on the multilayer, in which the interaction is characterized by the magnetic susceptibility. Most of experimental studies on this subject utilizes heavy metals, such as Pd and Pt, while a recent theory assumes a light metal as the normal metal. We aim to reconcile this discrepancy and achieve consistency with the experiments by studying magnetic susceptibility of Pd and Pt. We theoretically investigate the complex magnetic susceptibility of Pd and Pt by combining Tight-binding Approximation, Hubbard Model, and Density Functional Theory. We calculate a single band susceptibility for Pd and Pt. Our results show that the band structure is dominated by the d-orbitals and the susceptibility is enhanced by its large repulsive interaction.

UV-to-IR Absorption of Molecularly p-Doped Polythiophenes with Alkyl and Oligoether Side Chains: Experiment and Interpretation Based on Density Functional Theory.

The UV-to-IR transitions in p-doped poly(3-hexylthiophene) (P3HT) with alkyl side chains and polar polythiophene with tetraethylene glycol side chains are studied experimentally by means of the absorption spectroscopy and computationally using density functional theory (DFT) and tight-binding DFT. The evolution of electronic structure is calculated as the doping level is varied, while the roles of dopant ions, chain twisting, and π–π stacking are also considered, each of these having the effect of broadening the absorption peaks while not significantly changing their positions. The calculated spectra are found to be in good agreement with experimental spectra obtained for the polymers doped with a molybdenum dithiolene complex. As in other DFT studies of doped conjugated polymers, the electronic structure and assignment of optical transitions that emerge are qualitatively different from those obtained through earlier “traditional” approaches. In particular, the two prominent bands seen for the p-doped materials are present for both polarons and bipolarons/polaron pairs. The lowest energy of these transitions is due to excitation from the valence band to a spin-resolved orbitals located in the gap between the bands. The higher-energy band is a superposition of excitation from the valence band to a spin-resolved orbitals in the gap and an excitation between bands.

Thermochemical stabilities of giant fullerenes using density functional tight binding theory and isodesmic-type reactions.

We present a systematic assessment of the density functional tight binding (DFTB) method for calculating heats of formation of fullerenes with isodesmic-type reaction schemes. We show that DFTB3-D/3ob can accurately predict Δf H values of the 1812 structural isomers of C60 , reproduce subtle trends in Δf H values for 24 isolated pentagon rule (IPR) isomers of C84 , and predict Δf H values of giant fullerenes that are in effectively exact agreement with benchmark DSD-PBEP86/def2-QZVPP calculations. For fullerenes up to C320 , DFTB Δf H values are within 1.0 kJ mol-1 of DSD-PBEP86/def2-QZVPP values per carbon atom, and on a per carbon atom basis DFTB3-D/3ob yields exactly the same numerical trend of (Δf H [per carbon] = 722n-0.72 + 5.2 kJ mol-1 ). DFTB3-D/3ob is therefore an accurate replacement for high-level DHDFT and composite thermochemical methods in predicting of thermochemical stabilities of giant fullerenes and analogous nanocarbon architectures.

Comprehensive Assessment of GFN Tight-Binding and Composite Density Functional Theory Methods for Calculating Gas-Phase Infrared Spectra.

Vibrational spectroscopy is a valuable and widely used analytical tool for the characterization of chemical substances. We investigate the performance of semiempirical quantum mechanical GFN tight-binding and force-field methods for the calculation of gas-phase infrared spectra in comparison to experiment and low-cost (B3LYP-3c) density functional theory. A data set of 7247 experimental references was used to evaluate method performance based on automatic spectra comparison. Various quantitative spectral similarity measures were employed for the comparison between theory and experiment and for determining empirical scaling factors. It is shown that the scaling of atomic masses provides an accurate yet simple alternative to standard global frequency scaling in density functional theory (DFT) and semiempirical calculations. Furthermore, the method performance for 58 exemplary transition metal complexes was investigated. The efficient DFT composite method B3LYP-3c, that was introduced in the course of this work, was found to be excellently suited for general IR spectra calculations. The GFN1- and GFN2-xTB tight-binding methods clearly outperformed the PMx competitors. Conformational changes were investigated for a subset of the data and are found to have a mediocre strong influence on the simulated spectra suggesting that the corresponding elaborate sampling steps may be neglected in automated compound identification workflows.

Density Functional Tight Binding Theory Approach for the CO2 Reduction Reaction Paths on Anatase TiO2 Surfaces.

Herein, we have investigated the CO2 reduction paths on the (101) anatase TiO2 surface using an approach based on the density functional tight binding (DFTB) theory. We analyzed the reaction paths for the conversion of carbon dioxide to methane by performing a large number of calculations with intermediates placed in various orientations and locations at the surface. Our results show that the least stable intermediate is CO2H and therefore a key bottleneck is the reduction of CO2 to formic acid. Hydrogen adsorption is also weak and would also be a limiting factor, unless very high pressures of hydrogen are used. The results from our DFTB approach are in good agreement with the hybrid functional based density functional theory calculations presented in the literature.

Enantiomeric Transitions in the Chiral Cluster Au15 Studied by a Reaction Coordinate Deduced from Molecular Dynamics Simulations

A recently developed modified basin hopping (MBH) optimization algorithm, combined with an energy function calculated by the semiempirical density functional tight-binding (DFTB) theory, was applied to determine the lowest-energy structures of Aun clusters with size n = 3-20. It was predicted from the DFTB/MBH optimization algorithm calculations that clusters Au10, Au15, and Au18 exhibit chiral properties; i.e., each of these three clusters possesses the same energy value and associated with it are two nonsuperposable mirror-image clusters. In the potential energy landscape, there thus exist multidimensional barriers separating the two enantiomers, and this lowest-energy double-well morphology is surrounded by potential-energy minima of higher energies. In this paper, we have chosen to study the chiral cluster Au15 by employing an isothermal Brownian-type molecular dynamics simulation to discern in greater detail its conformational transition from one enantiomer, say left, to its right counterpart. To facilitate our analysis of the simulation data, we transpose the multidimensional configurational space description to a lower dimensional collective variable (CV) space spanned by two geometry-relevant CVs. The thermally driven progression and mechanism of enantiomeric transitions between the left and right enantiomers will be our main focus, and the strategy is to dissect the time development of the CVs collected from different sets of independent simulation runs. From simulation data, we found that an understanding of the dynamics of enantiomeric transitions needs first to seek out the left and right enantiomers through a molecular modeling and visualizing program, then to ferret out and identify between the left and right enantiomers a symmetrical structure, and finally to define from the latter a reaction coordinate. We showed in this work that this single reaction coordinate is predictive in unraveling the left ⇌ right enantiomeric transition events, providing a specific inkling of the transition time span and its associated distribution which can be checked further for its reasonableness by the autocorrelation function and a vibrational analysis, all of which shed light on the mechanisms of transition.

A theoretical and experimental study of castor oil-based inhibitor for corrosion inhibition of mild steel in acidic medium at elevated temperatures

The castor oil as a cheap and environmentally friendly source was used to prepare a novel green corrosion inhibitor for mild steel in acidic medium. The poor performance at high temperatures and low biodegradability are two important challenges of corrosion inhibitors that castor oil-based corrosion inhibitor (COCI) is designed to overcome them. The inhibition performance of COCI was evaluated by weight loss, open circuit potential (OCP), electrochemical impedance spectroscopy (EIS) and polarization techniques. The morphology of MS was examined by field emission scanning electron microscope (FESEM), energy dispersive X-ray (EDX) and atomic force microscopy (AFM). Furthermore, a comprehensive computational study was done to clarify the anticorrosive mechanism of COCI by molecular dynamics (MD) simulation, density functional based tight-binding (DFTB) approach and density functional theory (DFT). Maximum inhibition efficiency of 85 % and 91 % were achieved at 80 °C using 140 μM of COCI in EIS and polarization tests, respectively. COCI acted as a mixed type inhibitor and the adsorption of the inhibitor on the mild steel surface was chemisorption. The obtained results demonstrate that the COCI has a considerable ability to suppress acidic corrosion at all concentrations were studied here, particularly at high temperatures. Given the fact that several articles were reported a significant reduction in performance of corrosion inhibitors by increasing temperature, in this study was found out not only was there no decrease in inhibition efficiency of COCI at higher temperatures, but it depicted good performance at 80 °C. In addition, the high adsorption energy of COCI confirmed a strong interaction between Fe surface and inhibitor molecules and showed a significant effect of the urethane bonds and the triglyceride groups on molecular activity. The Fermi energy state indicated that the chemical adsorption occurred in the interaction of COCI-Fe surface by electron transition from inhibitor to Fe 3d band state. We expect the results of this work provide new opportunities for the design and synthesis of efficient corrosion inhibitors based on vegetable oils.

ELSI — An open infrastructure for electronic structure solvers

Routine applications of electronic structure theory to molecules and periodic systems need to compute the electron density from given Hamiltonian and, in case of non-orthogonal basis sets, overlap matrices. System sizes can range from few to thousands or, in some examples, millions of atoms. Different discretization schemes (basis sets) and different system geometries (finite non-periodic vs. infinite periodic boundary conditions) yield matrices with different structures. The ELectronic Structure Infrastructure (ELSI) project provides an open-source software interface to facilitate the implementation and optimal use of high-performance solver libraries covering cubic scaling eigensolvers, linear scaling density-matrix-based algorithms, and other reduced scaling methods in between. In this paper, we present recent improvements and developments inside ELSI, mainly covering (1) new solvers connected to the interface, (2) matrix layout and communication adapted for parallel calculations of periodic and/or spin-polarized systems, (3) routines for density matrix extrapolation in geometry optimization and molecular dynamics calculations, and (4) general utilities such as parallel matrix I/O and JSON output. The ELSI interface has been integrated into four electronic structure code projects (DFTB+, DGDFT, FHI-aims, SIESTA), allowing us to rigorously benchmark the performance of the solvers on an equal footing. Based on results of a systematic set of large-scale benchmarks performed with Kohn–Sham density-functional theory and density-functional tight-binding theory, we identify factors that strongly affect the efficiency of the solvers, and propose a decision layer that assists with the solver selection process. Finally, we describe a reverse communication interface encoding matrix-free iterative solver strategies that are amenable, e.g., for use with planewave basis sets. Program summaryProgram title: ELSI InterfaceCPC Library link to program files:http://dx.doi.org/10.17632/473mbbznrs.1Licensing provisions: BSD 3-clauseProgramming language: Fortran 2003, with interface to C/C++External routines/libraries: BLACS, BLAS, BSEPACK (optional), EigenExa (optional), ELPA, FortJSON, LAPACK, libOMM, MPI, MAGMA (optional), MUMPS (optional), NTPoly, ParMETIS (optional), PETSc (optional), PEXSI, PT-SCOTCH (optional), ScaLAPACK, SLEPc (optional), SuperLU_DISTNature of problem: Solving the electronic structure from given Hamiltonian and overlap matrices in electronic structure calculations.Solution method: ELSI provides a unified software interface to facilitate the use of various electronic structure solvers including cubic scaling dense eigensolvers, linear scaling density matrix methods, and other approaches.

DFTB Investigations on Transition Metals Doped TiO2 Quantum Dots

The structural, optical and electronic properties of different-sized pure as well as transition metals doped TiO2 quantum dots were investigated using density functional tight binding (DFTB) methods. The self-consistent charge density functional tight binding theory (SCC-DFTB) was used to model the TiO2 quantum dots (QDs) of increasing size up to 3.03 nm. The size dependence of density of states, position of molecular orbitals, and gap between highest orbital molecular orbitals (HOMO) and lowest unoccupied molecular orbitals (LUMO) of the QDs were studied. It appeared that, when the size of QDs approaches 3.03 nm, the energy gap narrowed down due to quantum confinement effects. The QD Ti27O54 was doped with 3d transition metals to further explore the possibility of tailoring the properties. The dopants introduced 3d impurity gap states which opens the opportunity to modify the band gap or positions of principal bands. The position of a shifted Fermi level was monitored to explore the changes in conductivity and n- or p-type behaviors of the doped QDs. It was observed that TiO2 QDs doped with V and Cr showed the n-type behavior while those doped with Sc, Mn and Fe showed p-type behavior. The findings are helpful to enhance the photocurrent efficiency of dye-sensitized solar cells and use the materials in electronic and optoelectronic properties.

Exploring the Polymorphism of Drostanolone Propionate.

2α-Methyl-4,5α-dihydrotestosterone 17β-propionate, known as drostanolone propionate or masteron, is a synthetic anabolic-androgenic steroid derived from dihydrotestosterone. The crystal structures of two polymorphs of drostanolone propionate have been determined by single crystal X-ray diffraction and both crystallizes in the monoclinic crystal system. One is belonging to the P21 space group, Z = 2, and has one molecule in the asymmetric unit while the second belongs to the I2 space group, Z = 4, and contains two molecules in the asymmetric unit. Another polymorph has been investigated by an X-ray powder diffraction method and solved by Parallel tempering/Monte Carlo technique and refined with the Rietveld method. This polymorph crystallizes in the orthorhombic P212121 space group, Z = 4 having one molecule in the asymmetric unit. The structural configuration analysis shows that the A, B, and C steroid rings exist as chair geometry, while ring D adopts a C13 distorted envelope configuration in all structures. For all polymorphs, the lattice energy has been computed by CLP (Coulomb-London-Pauli), and tight-binding density functional theory methods. Local electron correlation methods were used to estimate the role of electron correlation in the magnitude of the dimer energies. The nature of the intermolecular interactions has been analyzed by the SAPT0 energy decomposition methods as well as by Hirshfeld surfaces.

Structural rigidity accelerates quantum decoherence and extends carrier lifetime in porphyrin nanoballs: a time domain atomistic simulation.

Nonradiative electron–hole (e–h) recombination is the primary source of energy loss in photovoltaic cells and inevitably, it competes with the charge transfer process, leading to poor device performance. Therefore, much attention has to be paid for delaying such processes; increasing the excitonic lifetime may be a solution for this. Using the real-time, density functional tight-binding theory (DFTB) combined with nonadiabatic molecular dynamics (NAMD) simulations, we demonstrate the exciton relaxation phenomena of different metal-centered porphyrin nanoballs, which are supposed to be very important for the light-harvesting process. It has been revealed that the carrier recombination rate gradually decreases with the increase in the molecular stiffness by introducing metal-coordinating templating agents into the nanoball. Our simulation demonstrates that the lower atomic fluctuations lead to poorer electron–phonon nonadiabatic coupling in association with weak phonon modes and these as a whole are responsible for shorter quantum coherence and hence delayed recombination events. Our analysis is in good agreement with the recent experimental observation. By replacing the Zn metal center with a heavier Cd atom, a similar trend is observed; however, the rate slows down abruptly. The present simulation study provides the fundamental mechanism in detail behind the undesired energy loss during exciton recombination and suggests a rational design of impressive nanosystems for future device fabrication.

Coherent Real-Space Charge Transport Across a Donor-Acceptor Interface Mediated by Vibronic Couplings.

There is growing experimental and theoretical evidence that vibronic couplings, couplings between electronic and nuclear degrees of freedom, play a fundamental role in ultrafast excited-state dynamics in organic donor-acceptor hybrids. Whereas vibronic coupling has been shown to support charge separation at donor-acceptor interfaces, so far, little is known about its role in the real-space transport of charges in such systems. Here we theoretically study charge transport in thiophene:fullerene stacks using time-dependent density functional tight-binding theory combined with Ehrenfest molecular dynamics for open systems. Our results reveal coherent oscillations of the charge density between neighboring donor sites, persisting for ∼200 fs and promoting charge transport within the polymer stacks. At the donor-acceptor interface, vibronic wave packets are launched, propagating coherently over distances of more than 3 nm into the acceptor region. This supports previous experimental observations of long-range ballistic charge-carrier motion in organic photovoltaic systems and highlights the importance of vibronic coupling engineering as a concept for tailoring the functionality of hybrid organic devices.

N-Carbophenes: two-dimensional covalent organic frameworks derived from linear N-phenylenes

N-Carbophene (carbophene) is a novel class of two-dimensional covalent organic frameworks (2DCOF), based on linear N-phenylenes, that have moderate band gaps and low-mobility bands surrounding the Fermi energy; the simplest of which may have been recently synthesized. Using tight-binding density functional theory, the ground state configurations of carbophene single layers, bilayers, and bulk systems were determined. This work finds that carbophenes have formation energies per carbon atom similar to that of graphenylene. The similarity of formation energies between graphenylene and carbophene suggests that when trying to synthesize one, the other may also be synthesized. The formation energies could explain why the first reported synthesis of graphenylene also indicated that they may have synthesized 3-carbophene. Results contained in this work suggests that a carbophene was synthesized instead of graphenylene. The projected density of states (PDOS) demonstrates that the anti-aromatic nature of the cyclobutene units play a direct role in the creation of bands around the Fermi level, making this an exciting material in the theoretical understanding of the nature of aromatic bonds.

Geometry Optimization Skills of DFTB: Simple to Complex Molecular Systems

As the Density-functional theory (DFT) is an exact theory in principle for computing ground state electronic structures of the multi-electron and many-body systems, it's approximate variants currently being used are far from fail-safe. One of the very fundamental problem that has become apparent is its inability to account the substantial effect of van der Waals (vdW) type interactions exist in large molecular assemblies. To overcome such problem, a density functional tight binding (DFTB) theory whose fundamental formulation is based on the DFT but implements Slater−Kirkwood model and Slater−Koster files with a focus on solid state systems having vdW interactions as a binding force has been widely used recently. Its self-consistent charge ( SCC) approach is more promising theoretical model due to introducing self-consistent calculation of Mulliken charges. Present work is aimed at evaluating the geometry optimization skills of such DFTB method while applying to very simple to quite complex molecular systems of the order: water, benzene, crystalline 1,4-bis (tri-methylsilyl) benzene, and crystalline siloxaalkane. We fully optimized the isolated molecules of each of them as well as the unit cell geometries of the last two specimens and measured the dimensions of the particular sets of bond lengths, bond angles, and torsional angles in each optimized geometry. These values are found to be in an excellent agreement with the concerned experimental values. It makes the DFTB method very versatile and superb quantum mechanical model for computing ground state electronic structures.

Effects of surface hydration on the electron injection rate from graphene to anatase and rutile TiO2 surfaces

The effect of hydration of anatase (0 0 1), and rutile (0 0 1), (1 0 1) and (1 1 0) surfaces of TiO2 on the rate of electron capture from an adsorbed graphene sheet is elucidated by means of computation at a level of tight-binding density functional theory. Initially, the clean surfaces are allowed to interact with bulk water. Following an equilibration period, the first solvation shell is exposed to vacuum at 300 K for a sufficiently long period of time to allow any weakly bound water molecules to evaporate or be displaced from the surface. The resultant hydroxylated and originally clean surfaces are used for studying the interaction with a graphene sheet. Comparison of the calculated surface-graphene interaction energies and electron transfer probabilities of a cold electron to clean and hydroxylated surfaces shows that the water-sourced OH layer of these surfaces acts as an insulator in interfacial charge transport of the cold electron.

Atomistic simulations of charge transport in photoswitchable organic-graphene hybrids

Photoswitchable self-assembled monolayers (SAMs) in contact with a conductive or semiconductive layer can be used to remotely trigger changes in electrical current using light. In this study, we apply full-atomistic simulations to assess the changes in electronic structure and charge-transport properties of a graphene sheet in contact with an amorphous silica dielectric decorated by an azobenzene SAM. The simulations explicitly account for the structural and electrostatic disorder sourced by the dielectric, which turns out to be weakly affected by photoisomerization and spatially correlated over a length scale of 4–5 nm. Most interestingly, by combining large-scale (tight binding) density functional theory with Kubo–Greenwood quantum transport calculations, we predict that the trans-cis isomerization should induce a shift in surface electrostatic potential by a few tenths of a volt, accompanied by a variation in conductivity by a factor of about 3.