Long-range correction for density functional theory

Condensed Matter Systems Exposed to Radiation: Multiscale Theory, Simulations, and Experiment

This roadmap reviews the new, highly interdisciplinary research field studying the behavior of condensed matter systems exposed to radiation. The Review highlights several recent advances in the field and provides a roadmap for the development of the field over the next decade. Condensed matter systems exposed to radiation can be inorganic, organic, or biological, finite or infinite, composed of different molecular species or materials, exist in different phases, and operate under different thermodynamic conditions. Many of the key phenomena related to the behavior of irradiated systems are very similar and can be understood based on the same fundamental theoretical principles and computational approaches. The multiscale nature of such phenomena requires the quantitative description of the radiation-induced effects occurring at different spatial and temporal scales, ranging from the atomic to the macroscopic, and the interlinks between such descriptions. The multiscale nature of the effects and the similarity of their manifestation in systems of different origins necessarily bring together different disciplines, such as physics, chemistry, biology, materials science, nanoscience, and biomedical research, demonstrating the numerous interlinks and commonalities between them. This research field is highly relevant to many novel and emerging technologies and medical applications.

Theoretical description of photoinduced electron transfer in donor-acceptor supramolecular complexes based on carbon buckybowls

Herein, we explore, from a theoretical perspective, the nonradiative photoinduced processes (charge separation and energy transfer) within a family of donor-acceptor supramolecular complexes based on the electron-donor truxene-tetrathiafulvalene (truxTTF) derivative and a series of curved fullerene fragments (buckybowls) of different shapes and sizes (C30H12, C32H12, and C38H14) as electron acceptors that successfully combine with truxTTF via non-covalent interactions. The resulting supramolecular complexes (truxTTF·C30H12, truxTTF·C32H12, and truxTTF·C38H14) undergo charge-separation processes upon photoexcitation through charge-transfer states involving the donor and acceptor units. Despite the not so different size of the buckybowls, they present noticeable differences in the charge-separation efficiency owing to a complex decay post-photoexcitation mechanism involving several low-lying excited states of different natures (local and charge-transfer excitations), all closely spaced in energy. In this intricate scenario, we have adopted a theoretical approach combining electronic structure calculations at (time-dependent) density functional theory, a multistate multifragment diabatization method, the Marcus-Levitch-Jortner semiclassical rate expression, and a kinetic model to estimate the charge separation rate constants of the supramolecular heterodimers. Our outcomes highlight that the efficiency of the photoinduced charge-separation process increases with the extension of the buckybowl backbone. The supramolecular heterodimer with the largest buckybowl (truxTTF·C38H14) displays multiple and efficient electron-transfer pathways, providing a global photoinduced charge separation in the ultrafast time scale in line with the experimental findings. The study reported indicates that modifications in the shape and size of buckybowl systems can give rise to attractive novel acceptors for potential photovoltaic applications.

Singlet Fission as the Gateway to Triplet Generation in Heavy Atom-Free Organic Molecules

Triplet generations in heavy atom-free organic molecules are primarily revealed to proceed through singlet fissions (SFs) by investigating the contributions of SFs and intersystem crossings to the generation rates. The spin-flip long-range corrected time-dependent density functional theory calculations on 11 organic molecules known for triplet generation under photoirradiation are performed. The correlation between the descriptors for SF and the experimental singlet-to-triplet conversion rates strongly supports the predominance of SF progressions in all these molecules, corroborated by experimental observations of their triplet-triplet annihilations. Based on these findings, we propose updated conditions for SF progression: There is a high-absorption singlet state just above the triplet-triplet excitation of the chromophore dimer, or the singlet (triplet-triplet) excitation itself is responsible for photoabsorption. To the best of our knowledge, all organic molecules known for rapid triplet state generation fulfill these conditions.

Optoelectronic Properties of Nitrogen-Doped Hexagonal Graphene Quantum Dots: A First-Principles Study

Graphene quantum dots have been widely studied owing to their unique optical, electrical, and optoelectrical properties for various applications in solar devices. Here, we investigate the optoelectronic properties of hexagonal and nitrogen-doped graphene quantum dots using the first-principles method. We find that doping nitrogen atoms to hexagonal graphene quantum dots results in a significant red shift toward the visible light range as compared to that of the pristine graphene quantum dots, and the doped nitrogen atoms also induce a clear signature of anisotropy of the frontier orbitals induced by the electron correlation between the doped nitrogen atoms and their adjacent carbon atoms. Moreover, time-dependent density functional theory calculations with the M06-2X functional and 6-311++G(d,p) basis set reproduce well the experimental absorption spectra reported recently. These results provide us with a novel approach for more systematic investigations on next-generation solar devices with assembled quantum dots to improve their light selectivity as well as efficiency.

Designing Thiadiazoloquinoxaline-Based Conjugated Polymers for Efficient Organic Photovoltaics: A DFT/TDDFT Study

Clean and renewable energy development is becoming frontier research for future energy resources, as renewable energy offers sustainable and environmentally friendly alternatives to non-renewable sources such as fossil fuels. Among various renewable energy sources, tremendous progress has been made in converting solar energy to electric energy by developing efficient organic photovoltaics. Organic photovoltaic materials comprising conjugated polymers (CP) with narrow optical energy gaps are promising candidates for developing sustainable sources due to their potentially lower manufacturing costs. Organic semiconductor materials with a high electron affinity are required for many optoelectronic applications. We have designed a series of organic semiconductors comprised of cyclopentadithiophene as a donor and thiadiazoloquinoxaline (TQ) as an acceptor, varying the π-conjugation and TQ-derivatives. We have employed density functional theory (DFT) and time-dependent DFT (TDDFT) to evaluate the designed CP's optoelectronic properties, such as optical energy gap, dipole moment, and absorption spectra. Our DFT/TDDFT result shows that the energy gap of CPs is lowered and redshifted in the absorption spectra if there is no insertion of conjugation units such as thiophene and selenophene between donor and acceptor. In addition, selenophene shows relatively better redshift behavior compared to thiophene. Our work also provides rational insight into designing donor/acceptor-based CPs for organic solar cells.

Singlet fission initiating organic photosensitizations

The feasibility of singlet fission (SF) in organic photosensitizers is investigated through spin-flip long-range corrected time-dependent density functional theory. This study focuses on four major organic photosensitizer molecules: benzophenone, boron-dipyrromethene, methylene blue, and rose bengal. Calculations demonstrate that all these molecules possess moderate [Formula: see text]-stacking energies and closely-lying singlet (S) and quintet (triplet-triplet, TT) excitations, satisfying the essential conditions for SF: (1) Near-degenerate low-lying S and (TT) excitations with a significant S-T energy gap, and (2) Moderate [Formula: see text]-stacking energy of chromophores, slightly higher than solvation energy, enabling dissociation for triplet-state chromophore generation. Moreover, based on the El-Sayed rule, intersystem crossing is found to simultaneously proceed at very slow rates in all these photosensitizers. This is attributed to the fact that the lowest singlet excitation of the monomers partly involves [Formula: see text] transitions alongside the main [Formula: see text] transitions. The proposed mechanisms are strongly substantiated by comparisons with experimental studies.

Roles of Singlet Fission in the Photosensitization of Silicon Phthalocyanine

The roles of singlet fission in the triplet generation of silicon phthalocyanine (SiPc), a compound analogous to the IRDye700DX photosensitizer used in near-infrared photoimmunotherapy, are investigated by considering the energetical relation between the excitations of this compound. These excitations are obtained through spin-flip long-range corrected time-dependent density functional theory calculations. To initiate singlet fission, chromophores must meet two conditions: (1) near-degenerate low-lying singlet and quintet (triplet-triplet) excitations with a considerable energy gap of the lowest singlet and triplet excited states and (2) moderate π-stacking energy of chromophores, which is higher than but not far from the solvation energy, to facilitate the dissociation and generation of triplet-state chromophores. The present calculations demonstrate that SiPc satisfies both of these conditions after the formation of π-stacking irrespective of the presence of an axial ligand(s), suggesting that singlet fission plays a crucial role in the triplet generation process, although intersystem crossing occurs simultaneously at a very slow rate.

Total and orbital density-based analyses of molecules revealing long-range interaction regions

Total and orbital electron densities of molecules are explored for the effect of the long-range correction (LC) for density functional theory (DFT) exchange functionals by comparing to the effect of the ab initio coupled cluster singles and doubles (CCSD) method. Calculating the LC effect on the total electron densities shows that the LC stabilizes the electrons around the long-range interaction regions of kinetic energy density, which are assumed to be electrons other than free electrons and self-interacting electrons, while the CCSD method stabilizes the electrons in the long-range interaction regions in the vertical molecular planes. As a more precise test, the LC effect on orbital densities are compared to the CCSD effect on Dyson orbital densities. Surprisingly, these effects are similar for the unoccupied orbitals, indicating that the LC covers the effects required to reproduce the CCSD Dyson unoccupied orbitals. For exploring the discrepancies between these effects on the occupied orbitals, the photoionization cross sections are calculated as a direct test for the shapes of the HOMOs to investigate the differences between these effects on the occupied orbitals. Consequently, the LC clearly produces the canonical HOMOs close to the CCSD Dyson and experimental ones, except for the HOMO of benzene molecule that mixes with the HOMO - 1 for the CCSD Dyson orbitals. This indicates that the orbital analyses using the photoionization cross sections are available as a direct test for the quality of DFT functionals.

Understanding Photophysical Properties of Molecules Relevant in Organic Semiconductor Laser Diodes from Electron Localization Function-Tuned and Solvent-Tuned Range-Separated Functionals

Organic semiconductor laser diodes (OSLDs) are prevalent in optoelectronics because of their sustainable energy applications. Organic molecules used in such diodes are usually large; hence, their studies are computationally challenging with high-end benchmark methods. Computational methods with reliable accuracy and efficiency are always indispensable. In the present work, we have applied our computationally inexpensive, nonempirically tuned [electron localization function (ELF*) and solvent (Sol*)] range-separated (RS) functionals to study five molecules used in OSLDs. The emission energies in three different environments [toluene, CBP (4,4'-bis(n-carbazolyl)-1,1'-biphenyl) film, and gas] have been computed with the tuned functionals and compared with the experimental emission energies. ELF* and Sol* functionals can accurately reproduce emission energies in toluene and CBP film environments. On the other hand, both ELF* and IP-tuned functionals with excited-state geometry (IP*) perform better in the gas phase. In addition, a comparative study is performed between time-dependent density functional theory and the Tamm-Dancoff approximation. Along with the emission energy, oscillator strength values have also been reported. Different IP-tuned RS parameters were obtained with the ground- and excited-state geometries. Interestingly, it has been observed that the optimally tuned RS parameter with excited-state geometry (IP*) performs better compared to that with ground-state geometries (IP). Fractional occupation calculations show that the tuned functionals exhibit less localization and delocalization error. The study envisages that ELF* and Sol* functionals can be used to design future candidates for OSLDs.

Wetting vs Droplet Aggregation: A Broadband Rotational Spectroscopic Study of 3-Methylcatechol⋅⋅⋅Water Clusters

Two competing solvation pathways of 3-methylcatechol (MC), an atmospherically relevant aromatic molecule, with up to five water molecules were explored in detail by using a combination of broadband rotational spectroscopy and computational chemistry. Theoretically, two different pathways of solvation emerge: the commonly observed droplet pathway which involves preferential binding among the water molecules while the solute serves as an anchor point for the formation of a water cluster, and an unexpected wetting pathway which involves interactions between the water molecules and the aromatic face of MC, i.e., a wetting of the π-surface. Conclusive identification of the MC hydrate structures, and therefore the wetting pathway, was facilitated by rotational spectra of the parent MC hydrates and several H2 18 O and 13 C isotopologues which exhibit splittings associated with methyl internal rotation and/or water tunneling motions. Theoretical modelling and analyses offer insights into the tunneling and conversion barriers associated with the observed hydrate conformers and the nature of the non-covalent interactions involved in choosing the unusual wetting pathway.

Comparison of different density functional theory methods for the calculation of vibrational circular dichroism spectra

The determination of the absolute configuration (AC) of an organic molecule is still a challenging task for which the combination of spectroscopic with quantum-mechanical methods has become a promising approach. In this study, we investigated the accuracy of DFT methods (480 overall combinations of 15 functionals, 16 basis sets, and 2 solvation models) to calculate the VCD spectra of six chiral organic molecules in order to benchmark their capability to facilitate the determination of the AC.

Talarolide A and Talaropeptides A-D: Potential Marine-Derived Therapeutic Peptides with Interesting Chemistry and Biological Activity Studied through Density Functional Theory (DFT) and Conceptual DFT

Molecules sourced from marine environments hold immense promise for the development of novel therapeutic drugs, owing to their distinctive chemical compositions and valuable medicinal attributes. Notably, Talarolide A and Talaropeptides A-D have gained recent attention as potential candidates for pharmaceutical applications. This study aims to explore the chemical reactivity of Talarolide A and Talaropeptides A-D through the application of molecular modeling and computational chemistry techniques, specifically employing Conceptual Density Functional Theory (CDFT). By investigating their chemical behaviors, the study seeks to contribute to the understanding of the potential pharmacological uses of these marine-derived compounds. The molecular geometry optimizations and frequency calculations were conducted using the Density Functional Tight Binding (DFTBA) method. This was followed by a subsequent round of geometry optimization, frequency analysis, and computation of electronic properties and chemical reactivity descriptors. We employed the MN12SX/Def2TZVP/H2O model chemistry, utilizing the Gaussian 16 program and the SMD solvation model. The analysis of the global reactivity descriptors arising from CDFT was achieved as well as the graphical comparison of the dual descriptor DD revealing the areas of the molecules with more propensity to suffer a nucleophilic or electrophilic attack. Additionally, Molinspiration and SwissTargetPrediction were considered for the calculation of molecular characteristics and predicted biological targets. These include enzymes, nuclear receptors, kinase inhibitors, GPCR ligands, and ion channel modulators. The graphical results show that Talarolide A and the Talaropeptides A-D are likely to behave as protease inhibitors.

Irradiation of Plutonium Tributyl Phosphate Complexes by Ionizing Alpha Particles: A Computational Study

The PUREX solvent extraction process, widely used for recovering uranium and plutonium from spent nuclear fuel, utilizes an organic solvent composed of tributyl phosphate (TBP). The emission of ionizing particles such as alpha particles, resulting from the decay of plutonium, makes the organic solvent vulnerable to degradation. Here, we study the ultrashort time alpha irradiation of tributylphosphate (TBP) and Pu(NO3)4(TBP)2 complex formed in the PUREX process. Electron dynamics is propagated by Real-Time-Dependent Auxiliary Density Functional Theory (RT-TD-ADFT). We investigate the use of previously proposed absorption boundary conditions (ABC) in the molecular orbital space to treat secondary electron emission. Basis set and exchange correlation functional effects with ABC are reported as well as a detailed analysis of the ABC parametrization. Preliminary results on the water molecule and then on TBP show that the phenomenological nature of the ABC parameters necessitates selecting appropriate values for each system under study. Irradiation of free and complexed TBP shows an influence of the ligands on the variation of atomic charges on the femtosecond time scale. An accumulation of atomic charges in the alkyl chains of TBP is observed in the case where the nitrate groups are predominantly irradiated. In addition, we find that the Pu atom regains its electric charge very rapidly after being hit by the projectile, with the coordination sphere serving as an electron reservoir to preserve its formal redox state. This study paves the road toward a full understanding of the degradation of organic extracants employed in the nuclear industry.

Building on the strengths of a double-hybrid density functional for excitation energies and inverted singlet-triplet energy gaps

It is demonstrated that a double hybrid density functional approximation, ωB88PTPSS, that incorporates equipartition of density functional theory and the non-local correlation, however with a meta-generalized gradient approximation correlation functional, as well as with the range-separated exchange of ωB2PLYP, provides accurate excitation energies for conventional systems, as well as correct prescription of negative singlet-triplet gaps for non-conventional systems with inverted gaps, without any necessity for parametric scaling of the same-spin and opposite-spin non-local correlation energies. Examined over "safe" excitations of the QUESTDB set, ωB88PTPSS performs quite well for open-shell systems, correctly and fairly accurately [relative to equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) reference] predicts negative gaps for 50 systems with inverted singlet-triplet gaps, and is one of the leading performers for intramolecular charge-transfer excitations and achieves near-second-order approximate coupled cluster (CC2) and second-order algebraic diagrammatic construction quality for the Q1 and Q2 subsets. Subsequently, we tested ωB88PTPSS on two sets of real-life examples from recent computational chemistry literature-the low energy bands of chlorophyll a (Chl a) and a set of thermally activated delayed fluorescence (TADF) systems. For Chl a, ωB88PTPSS qualitatively and quantitatively achieves DLPNO-STEOM-CCSD-level performance and provides excellent agreement with experiment. For TADF systems, ωB88PTPSS agrees quite well with spin-component-scaled CC2 (SCS-CC2) excitation energies, as well as experimental values, for the gaps between the S₁ and T₁ excited states.

Singlet fission initiating triplet generations of BODIPY derivatives through [Formula: see text]-stacking: a theoretical study

The role of singlet fission (SF) in the triplet-state generation mechanism of 1,3,5,7-tetramethyl-boron-dipyrromethene derivatives is revealed by exploring the cause for the solvent dependence of the generation rate. Comparing the adsorption energy calculations of solvent molecules, i.e., cyclohexane, chloroform and acetonitrile molecules, to the derivatives with the [Formula: see text]-stacking energies of these derivatives surprisingly show that the hierarchy of the solvation energies and [Formula: see text]-stacking energies strongly correlates with the experimentally-suggested solvent dependence of the triplet-state generation of these derivatives for five and more adsorbing solvent molecules. Following this finding, the excitation spectra of these derivatives in acetonitrile solvent are explored using the proprietary spin-flip long-range corrected time-dependent density functional theory. It is, consequently, confirmed that the [Formula: see text]-stacking activates the second lowest singlet excitation to trigger the spin-allowed transition to the singlet doubly-excited tetraradical (TT)[Formula: see text] state, which generates the long-lived quintet (TT)[Formula: see text] state causing the SF. However, it is also found that the [Formula: see text]-stacking also get a slow intersystem crossing active around the intersections of the lowest singlet excitations with the lowest triplet T[Formula: see text] excitations in parallel with the SF due to the charge transfer characters of the lowest singlet excitations. These results suggest that SF initiates the triplet-state generations through near-degenerate low-lying singlet and (TT) excitations with a considerable singlet-triplet energy gap after the [Formula: see text]-stacking of chromophores stronger than but not far from the solvation. Since these derivatives are organic photosensitizers, this study proposes that SF should be taken into consideration in developing novel heavy atom-free organic photosensitizers, which will contribute to a variety of research fields such as medical care, photobiology, energy science, and synthetic chemistry.

Elbehairi S, Achar RR, Silina E, Stupin V, Manturova N, Glossman-Mitnik D, Shivamallu C, Kollur SP. Anticholinesterase activity of Areca Catechu: In Vitro and in silico green synthesis approach in search for therapeutic agents against Alzheimer's disease

For many years, the primary focus has been on finding effective treatments for Alzheimer's disease (AD), which has led to the identification of promising therapeutic targets. The necessity for AD stage-dependent optimal settings necessitated a herbal therapy strategy. The plant species Areca Catechu L. (AC) was selected based on the traditional uses against CNS-related diseases. AC leaf extract were prepared using a Soxhlet extraction method and hydroxyapatite nanoparticles (HAp-NPs) were synthesized from the same (AC-HAp-NPs). Powder X-ray diffractometer (XRD), scanning electron microscopy (SEM), transmission electron microscopy (TEM), selected area electron diffraction (SAED) and fourier transform infrared spectroscopy (FTIR) were used to confirm the structure and morphology of the as-prepared AC-HAp-NPs. The crystalline character of the AC-HAp-NPs was visible in the XRD pattern. The synthesized material was found to be nanoflake, with an average diameter of 15-20 nm, according to SEM analysis. The TEM and SAED pictures also revealed the form and size of AC-HAp-NPs. In vitro anti-acetylcholinesterase and butyrylcholinesterase (AChE and BChE) activities of hydroxyapatite nanoparticles produced from an AC leaf extract was tested in this study. When compared to control, AC-HAp-NPs had higher anti-AChE and BChE activity. The anti-acetylcholinesterase action of phytoconstituents generated from AC leaf extract was mediated by 4AQD and 4EY7, according to a mechanistic study conducted utilizing in silico research. The global and local descriptors, which are the underpinnings of Conceptual Density Functional Theory (CDFT), have been predicted through the MN12SX/Def2TZVP/H2O model chemistry to help in the comprehension of the chemical reactivity properties of the five ligands considered in this study. The CDFT experiments are supplemented by the calculation of several useful calculated pharmacokinetics indices, their expected biological targets connected to the bioavailability of the five ligands in order to further the goal of studying their bioactivity.

Reactive orbital energy theory serving a theoretical foundation for the electronic theory of organic chemistry

It is established that the reactive orbital energy theory (ROET) theoretically reproduces the rule-based electronic theory diagrams of organic chemistry by a comparative study on the charge transfer natures of typical organic carbon-carbon and carbon-heteroatom bond formation reactions: aldol, Mannich, α-aminooxylation, and isogyric reactions. The ROET, which is an expansion of the reaction electronic theories (e.g., the frontier orbital theory) in terms of orbital energies, elucidates the reactive orbitals driving reactions and the charge transferability indices of the reactions. Performing the ROET analyses of these reactions shows that the charge transfer directions given in the rule-based diagrams of the electronic theory are reproduced even for the functional groups of charge transfer destinations in all but only two processes for 38 reaction processes. The ROET analyses also make clear the detailed orbital-based pictures of these bond formation reactions: that is, the use of the out-of-plane antibonding π orbitals in acidic conditions (enol-mode) and in-plane antibonding π orbitals in basic conditions (enolate-mode), which explain the experimentally assumed mechanisms such as the π-bond formations in acidic conditions and σ-bond formations at α-carbons in basic conditions. Furthermore, the ROET analyses explicate that the methyl group initially accepts electrons and then donates them to the bond formations in the target reactions. It is, consequently, suggested that the ROET serves a theoretical foundation for the electronic theory of organic chemistry.

In Vitro Anticancer Screening, Molecular Docking and Antimicrobial Studies of Triazole-Based Nickel(II) Metal Complexes

Herein we describe the synthesis of a series of nickel(II) complexes (C1–C3) with Schiff bases (HL1–HL3) derived from 4-amino-5-mercapto-3-methyl-1,2,4-triazole and ortho/meta/para-nitrobenzaldehyde having composition [Ni(L)₂(H₂O)₂]. The obtained ligands and their complexes were characterized using physico-chemical techniques viz., elemental analysis, magnetic moment study, spectral (electronic, FT-IR, 1H-NMR) and thermal analysis. The elemental analysis and spectral analysis revealed that Schiff bases behave as monoanionic bidentate ligands towards the Ni(II) ion. Whereas, the magnetic moment study suggested the octahedral geometry of all the Ni(II) complexes. The thermal behavior of the complexes has been studied by thermogravimetric analysis and agrees well with the composition of complexes. Further, the biological activities such as antimicrobial and antifungal studies of the Schiff bases and Ni(II) complexes have been screened against bacterial species (Staphylococcus aureus and Pseudomonas aeruginosa) and fungal species (Aspergillus niger and Candida albicans) activity by MIC method, the results of which revealed that metal complexes exhibited significant antimicrobial activities than their respective ligands against the tested microbial species. Furthermore, the molecular docking technique was employed to investigate the active sites of the selected protein, which indeed helped us to screen the potential anticancer agents among the synthesized ligand and complexes. Further, these compounds have been screened for their in vitro anticancer activity using OVCAR-3 cell line. The results revealed that the complexes are more active than the ligands.

An insight into reactivity and bioactivity properties of quorum sensing peptides against PDE10A: a computational peptidology approach

Peptides are currently the most promising lead molecules. Quorum sensing peptides have a variety of structural features and are regularly exposed to post-translational modifications. Antiparkinsonian drugs lose their efficacy after a long period of use, and patients develop motor problems such as drug-induced dyskinesia (DIDs). The interaction between PDE10A and cAMP is necessary for dopamine neurotransmission and may play a role in Parkinson's disease pathogenesis. cAMP and cGMP are cyclic nucleotides that act as secondary messengers in the signal transduction pathway, influencing a range of CNS activities. PDE enzymes hydrolyze phosphodiester bonds to break down cAMP and cGMP, allowing them to control intracellular levels of these second messengers effectively. PDE expression, and hence cyclic nucleotide levels and their downstream targets, may change with age and in numerous age-related illnesses, including Parkinson's disease, according to mounting evidence. At the peak of dyskinesias, cyclic nucleotide levels were lower, and using phosphodiesterase inhibitors before antiparkinsonian medicines reduced the severity of dyskinesias. In a recent study, PapRIV was found to have the ability to activate BV-2 microglia cells, indicating that this quorum sensing peptide may play a role in gut-brain contact. As a result of the current in silico work, mainly focused on QSPs as a lead molecule for inhibiting PDE10A, the SRNAT QSP sequence has been a potent molecule in molecular docking and molecular dynamics simulations. Furthermore, we can test the efficiency of therapeutic components in vitro and in vivo utilizing this computational approach against PDE10A.

Excited-State Properties of Some Thermally Activated Delayed Fluorescence Emitters: Quest for an Accurate and Reliable Computational Method

Thermally activated delayed fluorescence (TADF) finds application in organic light-emitting diodes. The molecules exhibiting TADF are characterized by small singlet-triplet energy gaps that help reverse intersystem crossing. Recently, ionization potential (IP)-tuned range-separated (RS) density functionals have been well accepted for studying excited-state properties. In the present work, two efficient descriptor-based tuning schemes [electron localization function (ELF) and Sol] of RS density functionals have been used to accurately reproduce the excited-state properties of TADF emitters by performing a single self-consistent field calculation. The lowest singlet vertical excitation energies (E_VA(S₁)) and the vertical singlet-triplet energy gaps (ΔE_VST) are computed with ELF-, Sol-, and IP-tuned RS functionals (LC-BLYP, ωB97, ωB97X, and ωB97XD). Encouraging mean absolute deviations from the experimental values with ELF*-, Sol*-, and IP-tuned functionals are observed. Consistent performance of the non-empirical tuned functionals is noted in different solvent dielectrics. In addition to these, fractional occupation calculations have shown that our tuned functionals almost satisfy the energy linearity curve. Thus, ELF*- and Sol*-tuned functionals are promising and reliable alternatives in computing the excited-state properties. Considering the small experimental singlet-triplet gap, we recommend ELF* to calculate E_VA(S₁) and Sol* to calculate ΔE_VST.

Virtual Prospection of Marine Cyclopeptides as Therapeutics by Means of Conceptual DFT and Computational ADMET

Bioactive peptides are chemical compounds created through the covalent bonding of amino acids, known as amide or peptide bonds. Due to their unusual chemistry and various biological effects, marine bioactive peptides have garnered considerable research. The effectiveness of a bioactive marine peptide is attributed to its structural features, such as amino acid content and sequence, which vary depending on the degree of action. Cyclic peptides combine several favorable properties such as good binding affinity, target selectivity and low toxicity that render them an attractive modality for the development of therapeutics. The apratoxins are a class of molecules formed by a series of cyclic depsipeptides with potent cytotoxic activities. The objective of this research is to pursue a computational prospection of the molecular structures and properties of several cylopeptides of marine origin with potential therapeutic applications. The methodology will be based on the determination of the chemical reactivity descriptors of the studied molecules through the consideration of the Conceptual DFT model and validation of a particular model chemistry, MN12SX/Def2TZVP/H₂O. These studies will be complemented by a determination of the pharmacokinetics and ADMET parameters by resorting to certain cheminformatics tools.

Computational Chemistry Study of Natural Apocarotenoids and Their Synthetic Glycopeptide Conjugates as Therapeutic Drugs

The objective of the research to be presented in the chapter is the determination of the chemical reactivity properties of some natural apocarotenoids and their synthetic glycopeptide conjugates that could have the ability to inhibit SARS-CoV-2 replication. The study will be based on the consideration of the Conceptual DFT branch of Density Functional Theory (DFT) through the consideration of particular successful model chemistry which has been demonstrated as satisfying the Janak and Ionization Energy theorems within Generalized Gradient Approximation (GGA) theory. The research will be complemented by a report of the ADMET and pharmacokinetic properties hoping that this information could be of help in the development of new pharmaceutical drugs for fighting COVID-19.

Water-oxidation mechanism of cobalt phosphate co-catalyst in artificial photosynthesis: a theoretical study

The initial water-oxidation reaction mechanism of the hydrated cobalt phosphate (CoPi) co-catalyst, which is consistent with conventional experimental findings, is investigated for O-O bond and OOH formation in this study. Theoretical calculations of hydrated CoPi cluster models, which are validated by vibrational spectrum calculations, elucidate the roles of phosphate as a source of oxygen and deliverer of protons, both of which result in the spontaneous formation of an O-O bond after the release of two electrons and two protons. The calculations also show that OOH formation proceeds subsequently depending on the spin electronic states of the hydrated CoPi surface, and O₂ formation then spontaneously progresses after the release of two electrons and two protons. By theoretically tracing these processes, the initial water-oxidation reaction mechanism of the hydrated CoPi co-catalyst is proposed.

Conceptual DFT-Based Computational Peptidology, Pharmacokinetics Study and ADMET Report of the Veraguamides A–G Family of Marine Natural Drugs

As a continuation of our research on the chemical reactivity, pharmacokinetics and ADMET properties of cyclopeptides of marine origin with potential therapeutic abilities, in this work our already presented integrated molecular modeling protocol has been used for the study of the chemical reactivity and bioactivity properties of the Veraguamides A–G family of marine natural drugs. This protocol results from the estimation of the conceptual density functional theory (CDFT) chemical reactivity descriptors together with several chemoinformatics tools commonly considered within the process of development of new therapeutic drugs. CP-CDFT is a branch of computational chemistry and molecular modeling dedicated to the study of peptides, and it is a protocol that allows the estimation with great accuracy of the CDFT-based reactivity descriptors and the associated physical and chemical properties, which can aid in determining the ability of the studied peptides to behave as potential useful drugs. Moreover, the superiority of the MN12SX density functional over other long-range corrected density functionals for the prediction of chemical and physical properties in the presence of water as the solvent is clearly demonstrated. The research was supplemented with an investigation of the bioactivity of the molecular systems and their ADMET (absorption, distribution, metabolism, excretion, and toxicity) parameters, as is customary in medicinal chemistry. Some instances of the CDFT-based chemical reactivity descriptors’ capacity to predict the pKas of peptides as well as their potential as AGE inhibitors are also shown.

Computational peptidology approach to the study of the chemical reactivity and bioactivity properties of Aspergillipeptide D, a cyclopentapeptide of marine origin

Aspergillipeptide D is a cyclic pentapeptide isolated from the marine gorgonian Melitodes squamata-derived fungus Aspergillus sp. SCSIO 41501 that it has been shown to present moderate activity against herpes virus simplex type 1 (HSV-1). Thus, this paper presents the results of a computational study of this cyclopentapeptide's chemical reactivity and bioactivity properties using a CDFT-based computational peptidology (CDFT-CP) methodology, which is derived from combining chemical reactivity descriptors derived from Conceptual Density Functional Theory (CDFT) and some Cheminformatics tools which may be used. This results in an improvement of the virtual screening procedure by a similarity search allowing the identification and validation of the known ability of the peptide to act as a possible useful drug. This was followed by an examination of the drug's bioactivity and pharmacokinetics indices in relation to the ADMET (Absorption, Distribution, Metabolism, Excretion, and Toxicity) characteristics. The findings provide further evidence of the MN12SX density functional's superiority in proving the Janak and Ionization Energy theorems using the proposed KID approach. This has proven to be beneficial in accurately predicting CDFT reactivity characteristics, which aid in the understanding of chemical reactivity. The Computational Pharmacokinetics study revealed the potential ability of Aspergillipeptide D as a therapeutic drug through the interaction with different target receptors. The ADMET indices confirm this assertion through the absence of toxicity and good absorption and distribution properties.

Taking Advantage of a Systematic Energy Non-linearity Error in Density Functional Theory for the Calculation of Electronic Energy Levels

We present an approximate approach for the calculation of ionization potential (IP) and electron affinity (EA) by exploiting the complementary energy non-linearity errors for a species M and its one-electron-ionized counterpart (M⁺). Reasonable IPs and EAs are thus obtained by averaging the orbital energies of M and M⁺, even with a low-level method such as BLYP/6-31G(d). By combining the corrected IPs and EAs, we can further obtain reasonable excitation energies. The errors in uncorrected valence IPs and uncorrected virtual-orbital energies show systematic trends. These characteristics provide a convenient and computationally efficient avenue for qualitative estimation of these properties with single corrections for multiple IPs and excitation energies.

Charge-Separation and Charge-Recombination Rate Constants in a Donor-Acceptor Buckybowl-Based Supramolecular Complex: Multistate and Solvent Effects

The kinetics of the nonradiative photoinduced processes (charge-separation and charge-recombination) experimented in solution by a supramolecular complex formed by an electron-donating bowl-shaped truxene-tetrathiafulvalene (truxTTF) derivative and an electron-accepting fullerene fragment (hemifullerene, C30H12) has been theoretically investigated. The truxTTF·C30H12 heterodimer shows a complex decay mechanism after photoexcitation with the participation of several low-lying excited states of different nature (local and charge-transfer excitations) all close in energy. In this scenario, the absolute rate constants for all of the plausible charge-separation (CS) and charge-recombination (CR) channels have been successfully estimated using the Marcus-Levich-Jortner (MLJ) rate expression, electronic structure calculations, and a multistate diabatization method. The outcomes suggest that for a reasonable estimate of the CS and CR rate constants, it is necessary to include the following: (i) optimally tuned long-range (LC) corrected density functionals, to predict a correct energy ordering of the low-lying excited states; (ii) multistate effects, to account for the electronic couplings; and (iii) environmental solvent effects, to provide a proper stabilization of the charge-transfer excited states and accurate external reorganization energies. The predicted rate constants have been incorporated in a simple but insightful kinetic model that allows estimating global CS and CR rate constants in line with the most generalized three-state model used for the CS and CR processes. The values computed for the global CS and CR rates of the donor-acceptor truxTTF·C30H12 supramolecular complex are found to be in good agreement with the experimental values.

Enhancement of absorption capacity, optical and non-linear optical properties of graphene oxide nanosheet

We studied the absorption capacity, optoelectronic and non-linear optical (NLO) properties of graphene oxide nanosheet (GON) and its doped derivatives with aluminum (-Al) atoms. The investigations have been performed using three functionals (B3LYP, B3LYP-D3 and ωB97XD) of the density functional theory (DFT) associated to the basis set 6-31+G(d,p). Aluminum atoms were incorporated into GON at different sites in order to search for suitable candidates that could lead to the enhancement of NLO properties and decrease the band gap value of pristine graphene oxide. As per our molecular investigations, several doped molecular design schemes based on push-pull models of GON were proposed. The best electronic and NLO configurations responses highlight the doped derivatives which were obtained by replacement of carbon atoms which support the functional groups present on the honeycomb lattice of GON with -Al atoms.

One-to-One Correspondence between Reaction Pathways and Reactive Orbitals

The one-to-one correspondence between reaction pathways in the potential energy theory and reactive orbitals in the electronic theory for reactions is presented. In this study, the reactive orbital energy method is applied to the intrinsic reaction coordinates of the global reaction route map generated by an automated reaction path search method. The reactive orbital energy method specifies the pairs of occupied and unoccupied reactive orbitals driving chemical reactions and determines whether the reactions are electron transfer-driven or dynamics-driven. Surprisingly, it is found that the reactive orbital pairs are determined one by one for the electron transfer-driven reaction pathways from an identical molecule. The reactive orbital energy method is also found to provide the sophisticated interpretations of reactions for the electronic motions. This one-to-one correspondence is expected to trigger the unification of the potential energy theory and the electronic theory for reactions that have been independently developed.

General Design Rules for Bimetallic Platinum(II) Complexes

A series of platinum(II) bimetallic complexes were studied to investigate the effects of ligands on both the geometric and electronic structure. Modulating the Pt-Pt distance through the bridging ligand architecture was found to dictate the nature of the lowest energy electronic transitions, localized in one-half of the molecule or delocalized across the entire molecule. By reducing the separation between the platinum atoms, the lowest energy electronic transitions will be dominated by the metal-metal-to-ligand charge transfer transition. Conversely, by increasing the distance between the platinum atoms, the lowest electronic transition will be largely localized metal-to-ligand charge transfer or ligand centered in nature. Additionally, the cyclometalating ligands were observed to have a noticeable stabilizing effect on the triplet excited states as the conjugation increased, arising from geometric reorientation and increased electron delocalization of the ligands. Such stabilization of the triplet state energy has been shown to alter the excited state potential energy landscape as well as the excited state trajectory.

Analytical Gradients for Nuclear-Electronic Orbital Time-Dependent Density Functional Theory: Excited-State Geometry Optimizations and Adiabatic Excitation Energies

The computational investigation of photochemical processes often entails the calculation of excited-state geometries, energies, and energy gradients. The nuclear-electronic orbital (NEO) approach treats specified nuclei, typically protons, quantum mechanically on the same level as the electrons, thereby including the associated nuclear quantum effects and non-Born-Oppenheimer behavior into quantum chemistry calculations. The multicomponent density functional theory (NEO-DFT) and time-dependent DFT (NEO-TDDFT) methods allow efficient calculations of ground and excited states, respectively. Herein, the analytical gradients are derived and implemented for the NEO-TDDFT method and the associated Tamm-Dancoff approximation (NEO-TDA). The programmable equations for these analytical gradients as well as the NEO-DFT analytical Hessian are provided. The NEO approach includes the anharmonic zero-point energy (ZPE) and density delocalization associated with the quantum protons as well as vibronic mixing in geometry optimizations and energy calculations of ground and excited states. The harmonic ZPE associated with the other nuclei can be computed via the NEO Hessian. This approach is used to compute the 0-0 adiabatic excitation energies for a set of nine small molecules with all protons quantized, exhibiting slight improvement over the conventional electronic approach. Geometry optimizations of two excited-state intramolecular proton-transfer systems, [2,2'-bipyridyl]-3-ol and [2,2'-bipyridyl]-3,3'-diol, are performed with one and two quantized protons, respectively. The NEO calculations for these systems produce electronically excited-state geometries with stronger intramolecular hydrogen bonds and similar relative stabilities compared to conventional electronic methods. This work provides the foundation for nonadiabatic dynamics simulations of fundamental processes such as photoinduced proton transfer and proton-coupled electron transfer.

Spin-Scaled Range-Separated Double-Hybrid Density Functional Theory for Excited States

Our recently presented range-separated (RS) double-hybrid (DH) time-dependent density functional approach [J. Chem. Theory Comput. 17, 927 (2021)] is combined with spin-scaling techniques. The proposed spin-component-scaled (SCS) and scaled-opposite-spin (SOS) variants are thoroughly tested for almost 500 excitations including the most challenging types. This comprehensive study provides useful information not only about the new approaches but also about the most prominent methods in the DH class. The benchmark calculations confirm the robustness of the RS-DH ansatz, while several tendencies and deficiencies are pointed out for the existing functionals. Our results show that the SCS variant consistently improves the results, while the SOS variant preserves the benefits of the original RS-DH method reducing its computational expenses. It is also demonstrated that, besides our approaches, only the nonempirical functionals provide balanced performance for general applications, while particular methods are only suggested for certain types of excitations.

Assessing locally range-separated hybrid functionals from a gradient expansion of the exchange energy density

Locally range-separated hybrid (LRSH) functionals feature a real-space-dependent range separation function (RSF) instead of a system-independent range-separation parameter, which thus enables a more flexible admixture of exact exchange than conventional range-separated hybrid functionals. In particular, the development of suitable RSF models and exploring the capabilities of the LRSH approach, in general, are tasks that require further investigations and will be addressed in this work. We propose a non-empirical scheme based on a detailed scaling analysis with respect to a uniform coordinate scaling and on a short-range expansion of the range-separated exchange energy density to derive new RSF models from a gradient expansion of the exchange energy density. After optimizing a small set of empirical parameters introduced to enhance their flexibility, the resulting second- and fourth-order RSFs are evaluated with respect to atomic exchange energies, atomization energies, and transition barrier heights.

A Simple Range-Separated Double-Hybrid Density Functional Theory for Excited States

A simple and robust range-separated (RS) double-hybrid (DH) time-dependent density functional approach is presented for the accurate calculation of excitation energies of molecules within the Tamm-Dancoff approximation. The scheme can be considered as an excited-state extension of the ansatz proposed by Toulouse and co-workers [J. Chem. Phys. 2018, 148, 164105], which is based on the two-parameter decomposition of the Coulomb potential, for which both the exchange and correlation contributions are range-separated. A flexible and efficient implementation of the new scheme is also presented, which facilitates its extension to any combination of exchange and correlation functionals. The performance of the new approximation is tested for singlet excitations on several benchmark compilations and thoroughly compared to that of representative DH, RS hybrid, and RS DH functionals. The one-electron basis set dependence and computation times are also assessed. Our results show that the new approach improves on standard DHs in most cases, and it can provide a more robust and accurate alternative. In addition, on average, it noticeably surpasses the existing RS hybrid and RS DH functionals.

Fenton reaction mechanism generating no OH radicals in Nafion membrane decomposition

Mechanism of Fenton reaction, which is a most widely-used degradation test for organic materials using hydrogen peroxide (H\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_2$$\end{document}2O\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_2$$\end{document}2) and iron (Fe) cations, is revealed for the decomposition of hydrated Nafion membrane. This reaction mechanism has been assumed to generate OH radicals. For a doubly-hydrated Nafion membrane model, Fenton reaction with divalent and monovalent Fe (Fe\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{2+}$$\end{document}2+ and Fe\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^+$$\end{document}+) cation hydration complexes is explored for experimentally-supported hydration numbers using long-range correction for density functional theory. As a result, it is found that H\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_2$$\end{document}2O\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_2$$\end{document}2 coordinating to the Fe\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{2+}$$\end{document}2+ hydration complexes first approaches Nafion side chains in high humidity, then leads to the C–S bond dissociation of the side chain to produce carbonic acid group and sulfonic acid ion. On the other hand, once electron transfer proceeds between iron ions, the O–O bond of the coordinating H\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_2$$\end{document}2O\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_2$$\end{document}2 is extended, then the C–S bond is dissociated to produce trihydroxymethyl group and sulfur trioxide, which are rapidly transformed to carboxyl group and sulfonic acid ion in aquo. This mechanism is confirmed by the vibrational spectrum analysis of the decomposed product. Collective Nafion decomposition mechanisms also suggest that the decomposition reaction uses the recycle of generated Fe cation hydration complexes under acidic condition near membrane surface.

Rapid Prediction of Ultraviolet-Visible Spectra from Conventional (Non-Time-Dependent) Density Functional Theory Calculations

We present an approximate approach for the simulation of UV/vis spectra using conventional [non-time-dependent (non-TD)] DFT computations. It uses Kohn-Sham orbitals and orbital energies to estimate both the excitation energies and the associated oscillator strengths. For a wide range of systems from small molecules to large molecular dyes used in electrochromic and solar-cell applications, reasonable UV/vis spectra are generated, each with just two conventional DFT computations. The accuracy is generally comparable to what one would expect from TD-DFT calculations. In comparison to TD-DFT, the protocol of the present study provides an intuitive and notably more rapid means for simulating electronic absorption properties. It enables efficient screening of materials for a wide range of relevant applications.

Elucidating the binding mechanism of thione-containing mercaptopurine and thioguanine drugs to small gold clusters

Density functional theory methods were employed to clarify the adsorption/desorption behaviors of the thione-containing mercaptopurine and thioguanine drugs on the gold surface using both small Au₆ and Au₈ clusters as model reactants. Structural features, thermodynamic parameters, bonding characteristics, and electronic properties of the resulting complexes were investigated using the Perdew-Burke-Ernzerhof (PBE) and LC-BLYP functionals along with correlation-consistent basis sets, namely cc-pVDZ-PP for gold and cc-pVTZ for non-metals. Computed results show that the drug molecules tend to anchor on the gold cluster at the S atom with binding energies around -34 to -40 kcal/mol (in vacuum) and - 28 to -32 kcal/mol (in aqueous solution). As compared to Au₈ , Au₆ undergoes a shorter recovery time and a larger change of energy gap that could be converted to an electrical signal for selective detection of the drugs. Furthermore, interactions between the drugs and gold clusters are reversible processes and a drug release mechanism was also proposed. Accordingly, the drugs are able to separate from the gold surface due to either a slight change of pH in tumor cells or the presence of cysteine residues in protein matrices.

Unveiling the role of short-range exact-like exchange in the optimally tuned range-separated hybrids for fluorescence lifetime modeling

We propose and validate several variants of the optimally tuned range-separated hybrid functionals (OT-RSHs) including different density functional approximations for predicting the fluorescence lifetimes of different categories of fluorophores within the time-dependent density functional theory (TD-DFT) framework using both the polarizable continuum and state-specific solvation models. Our main idea originates from performing the optimal tuning in the presence of a contribution of the exact-like exchange at the short-range part, which, in turn, leads to the small values of the range-separation parameter, and computing the fluorescence lifetimes using the models including no or small portions of the short-range exact-like exchange. Particular attention is also paid to the influence of the geometries of emitters on fluorescence lifetime computations. It is shown that our developed OT-RSHs along with the polarizable continuum model can be considered as the promising candidates within the TD-DFT framework for the prediction of fluorescence lifetimes for various fluorophores. We find that the proposed models not only outperform their standard counterparts but also provide reliable data better than or comparable to the conventional hybrid functionals with both the fixed and interelectronic distance-dependent exact-like exchanges. Furthermore, it is also revealed that when the excited state geometries come into play, more accurate descriptions of the fluorescence lifetimes can be achieved. Hopefully, our findings can give impetus for future developments of OT-RSHs for computational modeling of other characteristics in fluorescence spectroscopy as well as for verification of the related experimental observations.

Synthesis and Synthetic Application of Chloro- and Bromofuroxans

Furoxans are potentially useful heteroaromatic units in pharmaceuticals and agrichemicals. However, the applications for furoxan-based compounds have been hampered due to the underdevelopment of their synthetic methods. Herein, we report a new synthetic approach for the synthesis of chloro- and bromofuroxans. The starting materials were dichloro- and dibromofuroxans, and the substituents were directly introduced to the furoxan ring in a modular fashion. The synthesized monohalofuroxans served as substrates for the installation of a second substituent to prepare further functionalized furoxans.

Vibrationally Mediated Stabilization of Electrons in Nonpolar Matter

We explore solvation of electrons in nonpolar matter, here represented by butadiene clusters. Isolated butadiene supports only the existence of transient anions (resonances). Two-dimensional electron energy loss spectroscopy shows that the resonances lead to an efficient vibrational excitation of butadiene, which can result into the almost complete loss of energy of the interacting electron. Cluster-beam experiments show that molecular clusters of butadiene form stable anions, however only at sizes of more than 9 molecular units. We have calculated the distribution of electron affinities of clusters using classical and path integral molecular dynamics simulations. There is almost a continuous transition from the resonant to the bound anions with an increase in cluster size. The comparison of the classical and quantum dynamics reveals that the electron binding is strongly supported by molecular vibrations, brought about by nuclear zero-point motion and thermal agitation. We also inspected the structure of the solvated electron, finding it well localized.

Performance of TDDFT Vertical Excitation Energies of Core-Substituted Naphthalene Diimides

We have evaluated the performance of various density functionals, covering generalized gradient approximation (GGA), global hybrid (GH) and range‐separated hybrid (RSH), using time dependent density functional theory (TDDFT) for computing vertical excitation energies against experimental absorption maximum (λ_max) for a set of 10 different core‐substituted naphthalene diimides (cNDI) recorded in dichloromethane. The computed excitation in case of GH PBE0 is most accurate while the trend is most systematic with RSH LCY‐BLYP compared to λ_max. We highlight the importance of including solvent effects for optimal agreement with the λ_max. Increasing the basis set size from TZ2P to QZ4P has a negligible influence on the computed excitation energies. Notably, RSH CAMY‐B3LYP gave the least error for charge‐transfer excitation. The poorest agreement with λ_max is obtained with semi‐local GGA functionals. Use of the optimally‐tuned RSH LCY‐BLYP* is not recommended because of the high computational cost and marginal improvement in results.

Structural Evolution and Stability Trend of Small-Sized Gold Clusters Aun (n = 20-30)

Structural evolution and stability pattern of pure neutral gold clusters Au_n in the small size range of n = 20-30 are examined using density functional theory (DFT) calculations. The equilibrium geometries are either confirmed or determined, and some new ground state structures are identified. The most stable configurations of Au₂₁-Au₂₃ sizes are formed by adding extra gold atoms to the highly stable pyramidal structure of Au₂₀, while flat-cage shapes are the best candidates for the global minima of both Au₂₄ and Au₂₅. For larger sizes of n = 26-30, pyramidal motifs tend to dominate the lower-lying population rather than tubular conformations as previously reported. The energy gaps, excitation energies, and exciton binding energies are also computed to test out the performance of the computational methods employed. Accordingly, a density functional with long-range exchange effects is highly recommended to quantitatively investigate both the ground and excited states of pure gold clusters.