Efficient Calculation of Electronic Absorption Spectra by Means of Intensity-Selected Time-Dependent Density Functional Tight Binding

Palladium Metal Nanocomposites Based on PEI-Functionalized Nitrogen-Doped Graphene Quantum Dots: Synthesis, Characterization, Density Functional Theory Modeling, and Cell Cycle Arrest Effects on Human Ovarian Cancer Cells

In this study, the synthesis, characterization, density functional theory calculations (DFT), and effect of polyethylenimine (PEI)-functionalized nitrogen-doped graphene quantum dots (PEI N-GQDs) and their palladium metal nanoparticles nanocomposites (PdNPs/PEI N-GQDs) on cancer cells were extensively investigated. The focus also includes investigating their cytotoxic and apoptotic effects on ovarian cancer cells, which pose a serious risk to women's health and have high death rates from delayed diagnosis, inadequate response to treatment, and decreased survival. Graphene quantum dots and their palladium nanocomposites were differentially effective against ovarian cancer cell lines. In particular, the smaller particle size and morphology of PdNPs/PEI N-GQDs nanocomposites compared with PEI N-GQDs probably enhance their activity through highly improved uptake by cells. These findings emphasize the importance of particle size in composite drugs for efficient cancer treatment. DFT results revealed that the Pd-containing nanocomposite, with a smaller highest occupied molecular orbital-lowest unoccupied molecular orbital gap, exhibited higher reactivity and anticancer effects in human ovarian cancer cell line, OVCAR-3. Significantly, the application of nanocomposites to ovarian cancer cells initiated apoptosis, offering valuable insights into the intricate interplay between nanomaterials and cancer biology.

Theoretical characterization and biological activity investigation of indirubins, cyclin dependent kinases inhibitors

Up to now, significant research efforts have been directed towards investigating indirubin and its derivatives as potential candidates for developing new compounds with multiple biological activities. In the present work, natural indirubin and numerous of its chemical derivatives referred to as indirubins have been investigated computationally using DFT method with the B3LYP/6-311 + G(d,p) level of theory, in order to reveal structure- biological activity relationship. We started with a structural properties description. Results analysis indicated that extra interaction sites were provided through the set of substitutions in compounds (1): Indirubin-3'-monoxime, (2): Indirubin-5-sulfonic acid, (3): 5-Nitro-indirubinoxime, (4): 5'-OH-5-nitro-indirubinoxime (AGM130), (5): 7-Bromo-5'-carboxyindirubin-3'-oxime, and (6): 7 BIO and consequently, extra hydrogen bonds may be formed with the active sites of molecular targets, such as GSK-3, CDKs, and Aurora kinases, as well as the aryl hydrocarbon receptor. Subsequently, to get more information on the electronic properties of indirubin and its analogues, HOMO, LUMO, Egap, and further electronic parameters were carried out. The indirubin derivatives showed an easier interaction with its environment than indirubin, the parent compound. The UV-Visible spectra of indirubin and compounds 1-6 were also produced using TD-DFT with B3LYP functional and 6-311 + G(2d,p) basis set. The relationship between absorption and chemical structure is discussed. Two phototoxic brominated compounds showed important absorption spectra modifications. It was also found that the main absorption bands of all compounds derived from π→π*(HOMO→LUMO) transitions.Communicated by Ramaswamy H. Sarma.

Two excited-state datasets for quantum chemical UV-vis spectra of organic molecules

We present two open-source datasets that provide time-dependent density-functional tight-binding (TD-DFTB) electronic excitation spectra of organic molecules. These datasets represent predictions of UV-vis absorption spectra performed on optimized geometries of the molecules in their electronic ground state. The GDB-9-Ex dataset contains a subset of 96,766 organic molecules from the original open-source GDB-9 dataset. The ORNL_AISD-Ex dataset consists of 10,502,904 organic molecules that contain between 5 and 71 non-hydrogen atoms. The data reveals the close correlation between the magnitude of the gaps between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO), and the excitation energy of the lowest singlet excited state energies quantitatively. The chemical variability of the large number of molecules was examined with a topological fingerprint estimation based on extended-connectivity fingerprints (ECFPs) followed by uniform manifold approximation and projection (UMAP) for dimension reduction. Both datasets were generated using the DFTB+ software on the "Andes" cluster of the Oak Ridge Leadership Computing Facility (OLCF).

Efficient simulation of resonance Raman spectra with tight-binding approximations to Density Functional Theory

Resonance Raman spectroscopy has long been established as one of the most sensitive techniques for detection, structure characterization and probing the excited-state dynamics of biochemical systems. However, the analysis of resonance Raman spectra is much facilitated when measurements are accompanied by Density Functional Theory (DFT) calculations which are expensive for large biomolecules. In this work, resonance Raman spectra are therefore computed with the Density Functional Tight-Binding (DFTB) method in the time-dependent excited-state gradient approximation. To test the accuracy of the tight-binding approximations, this method is first applied to typical resonance Raman benchmark molecules like β-carotene and compared to results obtained with pure and range-separated exchange-correlation (xc) functionals. We then demonstrate the efficiency of the approach by considering a computationally challenging heme variation. Overall, we find that the vibrational frequencies and excited-state properties (energies and gradients) which are needed to simulate the spectra are reasonably accurate and suitable for interpretation of experiments. We can therefore recommend DFTB as a fast computational method to interpret resonance Raman spectra.

Absorption Properties of Large Complex Molecular Systems: The DFTB/Fluctuating Charge Approach

We report on the first formulation of a novel polarizable QM/MM approach, where the density functional tight binding (DFTB) is coupled to the fluctuating charge (FQ) force field. The resulting method (DFTB/FQ) is then extended to the linear response within the TD-DFTB framework and challenged to study absorption spectra of large condensed-phase systems.

Two-photon absorption and triplet excited state quenching of near-IR region aza-BODIPY photosensitizers via a triphenylamine moiety despite heavy bromine atoms

Aza-BODIPY compounds with methoxy groups at -3 and -5 positons and triphenylamine moieties at -1, -7 positions with and without heavy bromine atoms at -2, -6 positions have been designed and synthesized.

Excited states in hydrogenated single-layer MoS2

Our calculations of the excitation spectrum of single-layer MoS₂ at several hydrogen coverages, using a density-matrix based time-dependent density-functional theory (TDDFT) show that the fully hydrogenated system is metallic, while at lower coverages the spectrum consists of spin-polarized partially filled localized mid-gap states. The calculated absorption spectrum of the system reveals standard excitonic peaks corresponding to the bound valence-band hole and conduction-band electron, as well as excitonic peaks that involve the mid-gap states. Binding energies of the excitons of the hydrogenated system are found to be relatively large (few tens of meV), making their experimental detection facile and suggesting hydrogenation as a knob for tuning the optical properties of single-layer MoS₂. Importantly, we find hydrogenation to suppress visible light photoluminescence, in agreement with experimental observations. In contrast, both Li and Na atoms transform the system into an n-doped non-magnetic semiconductor that does not allow excitonic states.

TD-DFTB study of optical properties of silver nanoparticle homodimers and heterodimers

The absorption spectra for face-centered cubic nanoparticle dimers at various interparticle distances are investigated using time-dependent density functional tight binding. Both homodimers and heterodimers are investigated in this work. By studying nanoparticles at various interparticle distances and analyzing their vertical excitations, we found that as the interparticle distance decreases, a red shift arises from contributions of the transition dipole moment that are aligned along the z-axis with nondegenerate features; blue shifts occur for peaks that originate from transition dipole moment components in the x and y directions with double degeneracy. When the nanoparticles are similar in size, the features in the absorption spectra become more sensitive to the interparticle distances. The best-fit curves from vertical excitation energy in the form of AR^-b for ΔE_redshift/ΔE_blueshift vs R are determined. In this way, we determined trends for absorption peak shifts and how these depend on the interparticle distance.

Systematic Computational Design and Optimization of Light Absorbing Dyes

We present a workflow to aid the discovery of new dyes for the role of a photosensitive unit in the dye-sensitized photo-electrochemical cells (DS-PECs). New structures are generated in a fully automated way using the Compound Attachment Tool (CAT) introduced in this work. These structures are characterized with efficient approximate density functional theory (DFT) methods, and molecules with favorable optical properties are suggested for possible further use in DS-PECs. As around 2500 structures are generated in this work, and as we aim for still larger volumes of compounds to screen in subsequent applications, we have assessed the reliability of low-cost screening methods and show that simplified time-dependent density functional theory (sTDDFT) provides a satisfying accuracy/cost ratio. From the dyes considered, we propose a set that can be suitable for panchromatic sensitization of the photoelectrode in DS-PECs to further increase DS-PEC efficiency.

Density-functional tight-binding: basic concepts and applications to molecules and clusters

The scope of this article is to present an overview of the Density Functional based Tight Binding (DFTB) method and its applications. The paper introduces the basics of DFTB and its standard formulation up to second order. It also addresses methodological developments such as third order expansion, inclusion of non-covalent interactions, schemes to solve the self-interaction error, implementation of long-range short-range separation, treatment of excited states via the time-dependent DFTB scheme, inclusion of DFTB in hybrid high-level/low level schemes (DFT/DFTB or DFTB/MM), fragment decomposition of large systems, large scale potential energy landscape exploration with molecular dynamics in ground or excited states, non-adiabatic dynamics. A number of applications are reviewed, focusing on -(i)- the variety of systems that have been studied such as small molecules, large molecules and biomolecules, bare orfunctionalized clusters, supported or embedded systems, and -(ii)- properties and processes, such as vibrational spectroscopy, collisions, fragmentation, thermodynamics or non-adiabatic dynamics. Finally outlines and perspectives are given.

Time-Dependent Long-Range-Corrected Density-Functional Tight-Binding Method Combined with the Polarizable Continuum Model

In this study, excited-state free energies and geometries were efficiently evaluated using a linear-response time-dependent long-range-corrected density-functional tight-binding method integrated with the polarizable continuum model (TD-LC-DFTB2/PCM). Although the LC-DFTB method required the evaluation of the exchange-type term, which was moderately computationally expensive, a single evaluation of the excited-state gradient for a system consisting of more than 1000 atoms in a vacuum was completed within 30 min using one CPU core. Benchmark calculations were conducted for 3-hydroxyflavone, which exhibits dual emission: the absorption and enol-form emission wavelengths calculated by TD-LC-DFTB2/PCM agreed well with those predicted based on the density functional theory using a long-range corrected functional; however, there was a large error in the predicted keto-form emission wavelength. Further benchmark calculations for more than 20 molecules indicated that the conventional TD-DFTB method underestimated the absorption and 0-0 transition energies compared with those which were measured experimentally, whereas the TD-LC-DFTB2 method systematically overestimated these metrics. Nevertheless, the agreement of the results of the TD-LC-DFTB2 method with those obtained by the CAM-B3LYP method demonstrates the potential of the TD-LC-DFTB2/PCM method. Moreover, changing the range separation parameter to 0.15 minimized this deviation.

Frequency Range Selection Method for Vibrational Spectra

Theoretical calculations of vibrational properties are widely used to explain and predict experimental spectra. However, with standard quantum chemical methods all molecular motions are considered, which is rather time-consuming for large molecules. Because typically only a specific spectral region is of experimental interest, we propose here an efficient method that allows calculation of only a selected frequency interval. After a computationally cheap low-level estimate of the molecular motions, the computational time is proportional to the number of normal modes needed to describe this frequency range. Results for a medium-sized molecule show a reduction in computational time of up to 1 order of magnitude with negligible loss in accuracy. We also show that still larger computational savings are possible by using an additional intensity-selection procedure.

TD-DFT and TD-DFTB Investigation of the Optical Properties and Electronic Structure of Silver Nanorods and Nanorod Dimers

Here, we perform theoretical investigation using time-dependent density functional theory (TD-DFT) and time-dependent density functional tight binding (TD-DFTB) for the electronic structure and optical properties of silver nanorods. TD-DFTB generally performs well for the accurate description of optical properties with respect to the size and type of dimer assembly of silver nanorods compared to TD-DFT. However, the energies and intensities of the longitudinal and transverse peaks of the nanorods are somewhat underestimated with TD-DFTB compared to the values calculated at the TD-DFT level. By exploiting the computational efficiency of TD-DFTB, we also extend our investigation to longer nanorods and their dimers containing up to ∼2000 atoms. Our results show that the coupling between nanorods and the resulting optical properties of the dimer assemblies are quite dependent on the length of the monomers. In all cases, the energy shifts in dimers as a function of the gap distance deviate significantly from the dipole-dipole interaction model. Moreover, a comparison of the best-fit curves for the dependence of the fractional shifts (Δλ/λ₀) on nanorod length indicates that the parameters of the plasmon ruler equation depend on the length of the nanorods and the type of the assembly rather than approaching a universal value. These insights are enabled by the computational efficiency of TD-DFTB and its ability to treat quantum mechanical effects in large nanorod dimer systems.

Davydov-type excitonic effects on the absorption spectra of parallel-stacked and herringbone aggregates of pentacene: Time-dependent density-functional theory and time-dependent density-functional tight binding

Exciton formation leads to J-bands in solid pentacene. Describing these exciton bands represents a challenge for both time-dependent (TD) density-functional theory (DFT) and for its semi-empirical analog, namely, for TD density-functional tight binding (DFTB) for three reasons: (i) solid pentacene and pentacene aggregates are bound only by van der Waals forces which are notoriously difficult to describe with DFT and DFTB, (ii) the proper description of the long-range coupling between molecules, needed to describe Davydov splitting, is not easy to include in TD-DFT with traditional functionals and in TD-DFTB, and (iii) mixing may occur between local and charge transfer excitons, which may, in turn, require special functionals. We assess how far TD-DFTB has progressed toward a correct description of this type of exciton by including both a dispersion correction for the ground state and a range-separated hybrid functional for the excited state and comparing the results against corresponding TD-CAM-B3LYP/CAM-B3LYP+D3 results. Analytic results for parallel-stacked ethylene are derived which go beyond Kasha's exciton model [M. Kasha, H. R. Rawls, and A. El-Bayoumi, Pure Appl. Chem. 11, 371 (1965)] in that we are able to make a clear distinction between charge transfer and energy transfer excitons. This is further confirmed when it is shown that range-separated hybrids have a markedly greater effect on charge-transfer excitons than on energy-transfer excitons in the case of parallel-stacked pentacenes. TD-DFT calculations with the CAM-B3LYP functional and TD-lc-DFT calculations lead to negligible excitonic corrections for the herringbone crystal structure, possibly because of an overcorrection of charge-transfer effects (CAM refers to Coulomb attenuated method). In this case, TD-DFT calculations with the B3LYP functional or TD-DFTB calculations parameterized to B3LYP give the best results for excitonic corrections for the herringbone crystal structure as judged from comparison with experimental spectra and with Bethe-Salpeter equation calculations from the literature.

Tight-binding approximations to time-dependent density functional theory - A fast approach for the calculation of electronically excited states

We propose a new method of calculating electronically excited states that combines a density functional theory based ground state calculation with a linear response treatment that employs approximations used in the time-dependent density functional based tight binding (TD-DFTB) approach. The new method termed time-dependent density functional theory TD-DFT+TB does not rely on the DFTB parametrization and is therefore applicable to systems involving all combinations of elements. We show that the new method yields UV/Vis absorption spectra that are in excellent agreement with computationally much more expensive TD-DFT calculations. Errors in vertical excitation energies are reduced by a factor of two compared to TD-DFTB.

Chemical Control and Spectral Fingerprints of Electronic Coupling in Carbon Nanostructures

The optical and electronic properties of atomically thin materials such as single-walled carbon nanotubes and graphene are sensitively influenced by substrates, the degree of aggregation, and the chemical environment. However, it has been experimentally challenging to determine the origin and quantify these effects. Here we use time-dependent density-functional-theory calculations to simulate these properties for well-defined molecular systems. We investigate a series of core-shell structures containing C₆₀ enclosed in progressively larger carbon shells and their perhydrogenated or perfluorinated derivatives. Our calculations reveal strong electronic coupling effects that depend sensitively on the interparticle distance and on the surface chemistry. In many of these systems we predict considerable orbital mixing and charge transfer between the C₆₀ core and the enclosing shell. We predict that chemical functionalization of the shell can modulate the electronic coupling to the point where the core and shell are completely decoupled into two electronically independent chemical systems. Additionally, we predict that the C₆₀ core will oscillate within the confining shell, at a frequency directly related to the strength of the electronic coupling. This low-frequency motion should be experimentally detectable in the IR region.

Crystal structure of human interferon-γ receptor 2 reveals the structural basis for receptor specificity

Interferon-γ receptor 2 is a cell-surface receptor that is required for interferon-γ signalling and therefore plays a critical immunoregulatory role in innate and adaptive immunity against viral and also bacterial and protozoal infections. A crystal structure of the extracellular part of human interferon-γ receptor 2 (IFNγR2) was solved by molecular replacement at 1.8 Å resolution. Similar to other class 2 receptors, IFNγR2 has two fibronectin type III domains. The characteristic structural features of IFNγR2 are concentrated in its N-terminal domain: an extensive π-cation motif of stacked residues KWRWRH, a NAG-W-NAG sandwich (where NAG stands for N-acetyl-D-glucosamine) and finally a helix formed by residues 78-85, which is unique among class 2 receptors. Mass spectrometry and mutational analyses showed the importance of N-linked glycosylation to the stability of the protein and confirmed the presence of two disulfide bonds. Structure-based bioinformatic analysis revealed independent evolutionary behaviour of both receptor domains and, together with multiple sequence alignment, identified putative binding sites for interferon-γ and receptor 1, the ligands of IFNγR2.

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.