Analysis and visualization of energy densities. II. Insights from linear-response time-dependent density functional theory calculations.
Ontology highlight
ABSTRACT:
SUBMITTER: Pei Z
PROVIDER: S-EPMC8258743 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:In this article, we report a scheme to analyze and visualize the energy density fluctuations during the real-time time-dependent density functional theory (RT-TDDFT) simulations. Using Ag4-N2 complexes as examples, it is shown that the grid-based Kohn-Sham energy density can be computed at each time step using a procedure from Nakai and coworkers. Then the instantaneous energy of each molecular fragment (such as Ag4 and N2) can be obtained by partitioning the Kohn-Sham energy densities using Becke or fragment-based Hirshfeld (FBH) scheme. A strong orientation-dependence is observed for the energy flow between the Ag4 cluster and a nearby N2 molecule in the RT-TDDFT simulations. Future applications of such an energy density analysis in electron dynamics simulations are discussed.
| S-EPMC7722154 | biostudies-literature
Project description:Linear scaling density functional theory (DFT) approaches to the electronic structure of materials are often based on the tendency of electrons to localize in large atomic and molecular systems. However, in many cases of actual interest, such as semiconductor nanocrystals, system sizes can reach a substantial extension before significant electron localization sets in, causing a considerable deviation from linear scaling. Herein, we address this class of systems by developing a massively parallel DFT approach which does not rely on electron localization and is formally quadratic scaling yet enables highly efficient linear wall-time complexity in the weak scalability regime. The method extends from the stochastic DFT approach described in Fabian et al. ( WIRES: Comp. Mol. Sci. 2019, e1412) but is entirely deterministic. It uses standard quantum chemical atom-centered Gaussian basis sets to represent the electronic wave functions combined with Cartesian real-space grids for some operators and enables a fast solver for the Poisson equation. Our main conclusion is that when a processor-abundant high-performance computing (HPC) infrastructure is available, this type of approach has the potential to allow the study of large systems in regimes where quantum confinement or electron delocalization prevents linear scaling.
| S-EPMC9009081 | biostudies-literature
Project description:We develop a second-order correction to commonly used density functional approximations (DFAs) to eliminate the systematic delocalization error. The method, based on the previously developed global scaling correction (GSC), is an exact quadratic correction to the DFA for the fractional charge behavior and uses the analytical second derivatives of the total energy with respect to fractional occupation numbers of the canonical molecular orbitals. For small and medium-size molecules, this correction leads to ground-state orbital energies that are a highly accurate approximation to the corresponding quasiparticle energies. It provides excellent predictions of ionization potentials, electron affinities, photoemission spectrum, and photoexcitation energies beyond previous approximate second-order approaches, thus showing potential for broad applications in computational spectroscopy.
| S-EPMC9367128 | biostudies-literature
Project description:Recently, exchange-correlation potentials in density functional theory were developed with the goal of providing improved band gaps in solids. Among them, the semilocal potentials are particularly interesting for large systems since they lead to calculations that are much faster than with hybrid functionals or methods like GW. We present an exhaustive comparison of semilocal exchange-correlation potentials for band gap calculations on a large test set of solids, and particular attention is paid to the potential HLE16 proposed by Verma and Truhlar. It is shown that the most accurate potential is the modified Becke-Johnson potential, which, most noticeably, is much more accurate than all other semilocal potentials for strongly correlated systems. This can be attributed to its additional dependence on the kinetic energy density. It is also shown that the modified Becke-Johnson potential is at least as accurate as the hybrid functionals and more reliable for solids with large band gaps.
| S-EPMC5423078 | biostudies-literature
Project description:Linear-scaling density functional theory (DFT) is an efficient method to describe the electronic structures of molecules, semiconductors, and insulators to avoid the high cubic-scaling cost in conventional DFT calculations. Here, we present a parallel implementation of linear-scaling density matrix trace correcting (TC) purification algorithm to solve the Kohn-Sham (KS) equations with the numerical atomic orbitals in the HONPAS package. Such a linear-scaling density matrix purification algorithm is based on the Kohn's nearsightedness principle, resulting in a sparse Hamiltonian matrix with localized basis sets in the DFT calculations. Therefore, sparse matrix multiplication is the most time-consuming step in the density matrix purification algorithm for linear-scaling DFT calculations. We propose to use the MPI_Allgather function for parallel programming to deal with the sparse matrix multiplication within the compressed sparse row (CSR) format, which can scale up to hundreds of processing cores on modern heterogeneous supercomputers. We demonstrate the computational accuracy and efficiency of this parallel density matrix purification algorithm by performing large-scale DFT calculations on boron nitrogen nanotubes containing tens of thousands of atoms.
| S-EPMC7726133 | biostudies-literature
Project description:The absorption and emission spectra of three azo sulfonamide compounds in different solvents were investigated theoretically by using response functions combined with density functional theory (DFT), while the solvent effect on the structure and the electronic transitions was determined using the integral equation formalism for the polarizable continuum model (IEF-PCM). The results show that the applied different exchange-correlation functionals can reproduce the experimental values well. DFT calculations of the title compounds showed that the H-bond formed between the solute and solvent molecules is one of the major causes of the reversible solvatochromism observed in measured spectra. This is due to a better stabilization of the neutral form than the zwitterionic form in the polar protic solvents, which is characteristic of the hypsochromic shift. On the other hand, the molecules considered exhibit a monotonic behavior regarding the polarity of the low-lying excited state (Δμg-CT) as a function of the solvent polarity. This dependence occurs in the case of the positive solvatochromism and confirms the thesis regarding the H-bond solute-solvent interactions. Theoretically determined values of the two-photon cross section revealed that the (σOF(2)) shows similar trends with changes in λabs, in contrast to <δ(OF)> values. In conclusion, the results demonstrate that the investigated molecules can be used successfully as fluorochromes in bioimaging.
| S-EPMC4427650 | biostudies-literature
Project description:Kynurenines, the main products of tryptophan catabolism, possess both prooxidant and anioxidant effects. Having multiple neuroactive properties, kynurenines are implicated in the development of neurological and cognitive disorders, such as Alzheimer's, Parkinson's, and Huntington's diseases. Autoxidation of 3-hydroxykynurenine (3HOK) and its derivatives, 3-hydroxyanthranilic acid (3HAA) and xanthommatin (XAN), leads to the hyperproduction of reactive oxygen species (ROS) which damage cell structures. At the same time, 3HOK and 3HAA have been shown to be powerful ROS scavengers. Their ability to quench free radicals is believed to result from the presence of the aromatic hydroxyl group which is able to easily abstract an electron and H-atom. In this study, the redox properties for kynurenines and several natural and synthetic antioxidants have been calculated at different levels of density functional theory in the gas phase and water solution. Hydroxyl bond dissociation enthalpy (BDE) and ionization potential (IP) for 3HOK and 3HAA appear to be lower than for xanthurenic acid (XAA), several phenolic antioxidants, and ascorbic acid. BDE and IP for the compounds with aromatic hydroxyl group are lower than for their precursors without hydroxyl group. The reaction rate for H donation to *O-atom of phenoxyl radical (Ph-O*) and methyl peroxy radical (Met-OO*) decreases in the following rankings: 3HOK ~ 3HAA > XAAOXO > XAAENOL. The enthalpy absolute value for Met-OO* addition to the aromatic ring of the antioxidant radical increases in the following rankings: 3HAA* < 3HOK* < XAAOXO* < XAAENOL*. Thus, the high free radical scavenging activity of 3HAA and 3HOK can be explained by the easiness of H-atom abstraction and transfer to O-atom of the free radical, rather than by Met-OO* addition to the kynurenine radical.
| S-EPMC5115656 | biostudies-literature
Project description:Herein, we report a pair of regioselective N 1- and N 2 -alkylations of a versatile indazole, methyl 5-bromo-1H-indazole-3-carboxylate (6) and the use of density functional theory (DFT) to evaluate their mechanisms. Over thirty N 1- and N 2-alkylated products were isolated in over 90% yield regardless of the conditions. DFT calculations suggest a chelation mechanism produces the N 1-substituted products when cesium is present and other non-covalent interactions (NCIs) drive the N 2-product formation. Methyl 1H-indazole-7-carboxylate (18) and 1H-indazole-3-carbonitrile (21) were also subjected to the reaction conditions and their mechanisms were evaluated. The N 1- and N 2-partial charges and Fukui indices were calculated for compounds 6, 18, and 21 via natural bond orbital (NBO) analyses which further support the suggested reaction pathways.
| S-EPMC11318628 | biostudies-literature
Project description:The double-hybrid (DH) time-dependent density functional theory is extended to core excitations. Two different DH formalisms are presented utilizing the core-valence separation (CVS) approximation. First, a CVS-DH variant is introduced relying on the genuine perturbative second-order correction, while an iterative analogue is also presented using our second-order algebraic-diagrammatic construction [ADC(2)]-based DH ansatz. The performance of the new approaches is tested for the most popular DH functionals using the recently proposed XABOOM [J. Chem. Theory Comput.2021, 17, 1618] benchmark set. In order to make a careful comparison, the accuracy and precision of the methods are also inspected. Our results show that the genuine approaches are highly competitive with the more advanced CVS-ADC(2)-based methods if only excitation energies are required. In contrast, as expected, significant differences are observed in oscillator strengths; however, the precision is acceptable for the genuine functionals as well. Concerning the performance of the CVS-DH approaches, the PBE0-2/CVS-ADC(2) functional is superior, while its spin-opposite-scaled variant is also recommended as a cost-effective alternative. For these approaches, significant improvements are realized in the error measures compared with the popular CVS-ADC(2) method.
| S-EPMC9979613 | biostudies-literature
Project description:In this paper we review recent advances in computational chemistry and specifically focus on the chemical description of heme proteins and synthetic porphyrins that act as both mimics of natural processes and technological uses. These are challenging biochemical systems involved in electron transfer as well as biocatalysis processes. In recent years computational tools have improved considerably and now can reproduce experimental spectroscopic and reactivity studies within a reasonable error margin (several kcal·mol(-1)). This paper gives recent examples from our groups, where we investigated heme and synthetic metal-porphyrin systems. The four case studies highlight how computational modelling can correctly reproduce experimental product distributions, predicted reactivity trends and guide interpretation of electronic structures of complex systems. The case studies focus on the calculations of a variety of spectroscopic features of porphyrins and show how computational modelling gives important insight that explains the experimental spectra and can lead to the design of porphyrins with tuned properties.
| S-EPMC4848975 | biostudies-literature