Time-Dependent Density-Functional Theory

How do spin-scaled double hybrids designed for excitation energies perform for noncovalent excited-state interactions? An investigation on aromatic excimer models

Time-dependent double hybrids with spin-component or spin-opposite scaling to their second-order perturbative correlation correction have demonstrated competitive robustness in the computation of electronic excitation energies. Some of the most robust are those recently published by our group (M. Casanova-Páez, L. Goerigk, J. Chem. Theory Comput. 2021, 20, 5165). So far, the implementation of these functionals has not allowed correctly calculating their ground-state total energies. Herein, we define their correct spin-scaled ground-state energy expressions which enables us to test our methods on the noncovalent excited-state interaction energies of four aromatic excimers. A range of 22 double hybrids with and without spin scaling are compared to the reasonably accurate wavefunction reference from our previous work (A. C. Hancock, L. Goerigk, RSC Adv. 2023, 13, 35964). The impact of spin scaling is highly dependent on the underlying functional expression, however, the smallest overall errors belong to spin-scaled functionals with range separation: SCS- and SOS- ω PBEPP86, and SCS-RSX-QIDH. We additionally determine parameters for DFT-D3(BJ)/D4 ground-state dispersion corrections of these functionals, which reduce errors in most cases. We highlight the necessity of dispersion corrections for even the most robust TD-DFT methods but also point out that ground-state based corrections are insufficient to completely capture dispersion effects for excited-state interaction energies.

Ultrafast fragmentation of highly-excited doubly-ionized deoxyribose: role of the liquid water environment

Ab initio molecular dynamics simulations are used to investigate the fragmentation dynamics following the double ionization of 2-deoxy-D-ribose (DR), a major component in the DNA chain. Different ionization scenarios are considered to provide a complete picture. First focusing on isolated DR2+, fragmentation patterns are determined for the ground electronic state, adding randomly distributed excitation energy to the nuclei. These patterns differ for the two isomers studied. To compare thermal and electronic excitation effects, Ehrenfest dynamics are also performed, allowing to remove the two electrons from selected molecular orbitals. Two intermediate-energy orbitals, localized on the carbon chain, were selected. The dissociation pattern corresponds to the most frequent pattern obtained when adding thermal excitation. On the contrary, targeting the four deepest orbitals, localized on the oxygen atoms, leads to selective ultrafast C-O and/or O-H bond dissociation. To probe the role of environment, a system consisting of a DR molecule embedded in liquid water is then studied. The two electrons are removed from either the DR or the water molecules directly linked to the sugar through hydrogen bonds. Although the dynamics onset is similar to that of isolated DR when removing the same deep orbitals localized on the sugar oxygen atoms, the subsequent fragmentation patterns differ. Sugar damage also occurs following the Coulomb explosion of neighboring H2O2+ molecules due to interaction with the emitted O or H atoms.

Acenaphthylene-Based Chromophores for Dye-Sensitized Solar Cells: Synthesis, Spectroscopic Properties, and Theoretical Calculations

A set of acenaphthylene dyes with arylethynyl π-bridges was tested for dye-sensitized solar cells (DSSCs). Crucial steps for the extension of the conjugated system from the acenaphylene core involved Sonogashira coupling reactions. Phenyl, thiophene, benzotriazole, and thieno-[3,2-b]thiophene moieties were employed to extend the conjugation of the π-bridges. The systems were characterized by cyclic voltammetry and by UV-vis absorption and emission. The spectroscopic characterization showed that the last three bridges resulted in red-shifted absorption and emission spectra relative to the parent phenyl-bridged compound, in accordance with TD-DFT calculations. The phenylethynyl derivative 6a achieved a conversion efficiency of 2.51% with Voc, Jsc, and FF values of 0.365 V, 13.32 mA/cm2, and 0.52, respectively. The efficiency of this compound improved to 3.15% with the addition of CDCA (10 mM), representing the best efficiency result in this study. The overall conversion efficiency of the other aryl derivatives 6b-d proved to be significantly inferior (14-40%) to that of 6a due to a significant decrease of Jsc.

Heavy element incorporation in nitroimidazole radiosensitizers: molecular-level insights into fragmentation dynamics

The present study investigates the photofragmentation behavior of iodine-enhanced nitroimidazole-based radiosensitizer model compounds in their protonated form using near-edge X-ray absorption mass spectrometry and quantum mechanical calculations. These molecules possess dual functionality: improved photoabsorption capabilities and the ability to generate species that are relevant to cancer sensitization upon photofragmentation. Four samples were investigated by scanning the generated fragments in the energy regions around C 1s, N 1s, O 1s, and I 3d-edges with a particular focus on NO2+ production. The experimental summed ion yield spectra are explained using the theoretical near-edge X-ray absorption fine structure spectrum based on density functional theory. Born-Oppenheimer-based molecular dynamics simulations were performed to investigate the fragmentation processes.

Mechanistic elucidation of the degradation and transformation of hydroxy-α-sanshool and its conformers as the pungent dietary components in Sichuan pepper: A DFT study

To better understand the structural changes of sanshool pungent dietary components during the process and preservation of Sichuan pepper and pungent foods, the mechanistic insights into the intrinsic degradation and transformation of 16 hydroxy-α-sanshool conformers have been explored computationally. Our results have revealed that increasing the cis-CC bond numbers in the most stable all-trans hydroxy-β-sanshool structure causes the maximum 34.21 kJ/mol conformational energetic difference, and the existent probability of C2nnn would be lower than that of C1nnn (n = 1,2). The isomerization between the conformers could be much easier when they are excited by light radiation, as the strength of the CC bonds and their connected CC bonds becomes significantly close, and the relative energies among conformers are largely reduced. Besides, the different combination of cis/trans-CC bonds changes the delocalization degree of molecular Frontier orbitals, which consequently causes the different photochemical stability. Finally, the possible molecular oxidation degradation mechanism is proposed.

Theoretical design of durable and strong polycarbonates against photodegradation

The photodegradation mechanism of polycarbonate (PC) was investigated by quantum chemistry, and a novel antidegradation molecular design using substituents was proposed. It was demonstrated that electron-withdrawing substituents in the phenyl moiety controlled bond alternation, leading to inhibition of the O-C bond cleavage in the carbonate moiety. These results provide a promising alternative for durable PC synthesis.

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.

Effects of Oxygen Adsorption on the Optical Properties of Ag Nanoparticles

Plasmonic metal nanoparticles are efficient light harvesters with a myriad of sensing- and energy-related applications. For such applications, the optical properties of nanoparticles of metals such as Cu, Ag, and Au can be tuned by controlling the composition, particle size, and shape, but less is known about the effects of oxidation on the plasmon resonances. In this work, we elucidate the effects of O adsorption on the optical properties of Ag particles by evaluating the thermodynamic properties of O-decorated Ag particles with calculations based on the density functional theory and subsequently computing the photoabsorption spectra with a computationally efficient time-dependent density functional theory approach. We identify stable Ag nanoparticle structures with oxidized edges and a quenching of the plasmonic character of the metal particles upon oxidation and trace back this effect to the sp orbitals (or bands) of Ag particles being involved both in the plasmonic excitation and in the hybridization to form bonds with the adsorbed O atoms. Our work has important implications for the understanding and application of plasmonic metal nanoparticles and plasmon-mediated processes under oxidizing environments.

Noncovalently bound excited-state dimers: a perspective on current time-dependent density functional theory approaches applied to aromatic excimer models

Excimers are supramolecular systems whose binding strength is influenced by many factors that are ongoing challenges for computational methods, such as charge transfer, exciton coupling, and London dispersion interactions. Treating the various intricacies of excimer binding at an adequate level is expected to be particularly challenging for time-dependent Density Functional Theory (TD-DFT) methods. In addition to well-known limitations for some TD-DFT methods in the description of charge transfer or exciton coupling, the inherent London dispersion problem from ground-state DFT translates to TD-DFT. While techniques to appropriately treat dispersion in DFT are well-developed for electronic ground states, these dispersion corrections remain largely untested for excited states. Herein, we aim to shed light on current TD-DFT methods, including some of the newest developments. The binding of four model excimers is studied across nine density functionals with and without the application of additive dispersion corrections against a wave function reference of SCS-CC2/CBS(3,4) quality, which approximates select CCSDR(3)/CBS data adequately. To our knowledge, this is the first study that presents single-reference wave function dissociation curves at the complete basis set level for the assessed model systems. It is also the first time range-separated double-hybrid density functionals are applied to excimers. In fact, those functionals turn out to be the most promising for the description of excimer binding followed by global double hybrids. Range-separated and global hybrids-particularly with large fractions of Fock exchange-are outperformed by double hybrids and yield worse dissociation energies and inter-molecular equilibrium distances. The deviation between each assessed functional and reference increases with system size, most likely due to missing dispersion interactions. Additive dispersion corrections of the DFT-D3(BJ) and DFT-D4 types reduce the average errors for TD-DFT methods but do so inconsistently and therefore do not offer a black-box solution in their ground-state parametrised form. The lack of appropriate description of dispersion effects for TD-DFT methods is likely hindering the practical application of the herein identified more efficient methods. Dispersion corrections parametrised for excited states appear to be an important next step to improve the applicability of TD-DFT methods and we hope that our work assists with the future development of such corrections.

Excited State Properties of Point Defects in Semiconductors and Insulators Investigated with Time-Dependent Density Functional Theory

We present a formulation of spin-conserving and spin-flip hybrid time-dependent density functional theory (TDDFT), including the calculation of analytical forces, which allows for efficient calculations of excited state properties of solid-state systems with hundreds to thousands of atoms. We discuss an implementation on both GPU- and CPU-based architectures along with several acceleration techniques. We then apply our formulation to the study of several point defects in semiconductors and insulators, specifically the negatively charged nitrogen-vacancy and neutral silicon-vacancy centers in diamond, the neutral divacancy center in 4H silicon carbide, and the neutral oxygen-vacancy center in magnesium oxide. Our results highlight the importance of taking into account structural relaxations in excited states in order to interpret and predict optical absorption and emission mechanisms in spin defects.

Time-Dependent Density Functional Theory for X-ray Absorption Spectra: Comparing the Real-Time Approach to Linear Response

We simulate X-ray absorption spectra at elemental K-edges using time-dependent density functional theory (TDDFT) in both its conventional linear-response implementation and its explicitly time-dependent or "real-time" formulation. Real-time TDDFT simulations enable broadband spectra calculations without the need to invoke frozen occupied orbitals ("core/valence separation"), but we find that these spectra are often contaminated by transitions to the continuum that originate from lower-energy core and semicore orbitals. This problem becomes acute in triple-ζ basis sets, although it is sometimes sidestepped in double-ζ basis sets. Transitions to the continuum acquire surprisingly large dipole oscillator strengths, leading to spectra that are difficult to interpret. Meaningful spectra can be recovered by means of a filtering technique that decomposes the spectrum into contributions from individual occupied orbitals, and the same procedure can be used to separate L- and K-edge spectra arising from different elements within a given molecule. In contrast, conventional linear-response TDDFT requires core/valence separation but is free of these artifacts. It is also significantly more efficient than the real-time approach, even when hundreds of individual states are needed to reproduce near-edge absorption features and even when Padé approximants are used to reduce the real-time simulations to just 2-4 fs of time propagation. Despite the cost, the real-time approach may be useful to examine the validity of the core/valence separation approximation.

Tautomeric equilibrium and spectroscopic properties of 8-azaguanine revealed by quantum chemistry methods

8-azaguanine is a triazolopyrimidine nucleobase analog possessing potent antibacterial and antitumor activities, and it has been implicated as a lead molecule in cancer and malaria therapy. Its intrinsic fluorescence properties can be utilized for monitoring its interactions with biological polymers like proteins or nucleic acids. In order to better understand these interactions, it is important to know the tautomeric equilibrium of this compound. In this work, the tautomeric equilibrium of all natural neutral and anionic compound forms (except highly improbable imino-enol tautomers) as well as their methyl derivatives and ribosides was revealed by quantum chemistry methods. It was shown that, as expected, tautomers protonated at positions 1 and 9 dominate neutral forms both in gas phase and in aqueous solution. 8-azaguanines methylated at any position of the triazole ring are protonated at position 1. The computed vertical absorption and emission energies are in very good agreement with the experimental data. They confirm the validity of the assumption that replacing the proton with the methyl group does not significantly change the positions of absorption and fluorescence peaks.

Fused Indacene Dimers

Polycyclic hydrocarbons consisting of two or more directly fused antiaromatic subunits are rare due to their high reactivity. However, it is important to understand how the interactions between the antiaromatic subunits influence the electronic properties of the fused structure. Herein, we present the synthesis of two fused indacene dimer isomers: s-indaceno[2,1-a]-s-indacene (s-ID) and as-indaceno[3,2-b]-as-indacene (as-ID), containing two fused antiaromatic s-indacene or as-indacene units, respectively. Their structures were confirmed by X-ray crystallographic analysis. 1 H NMR/ESR measurements and DFT calculations revealed that both s-ID and as-ID have an open-shell singlet ground state. However, while localized antiaromaticity was observed in s-ID, as-ID showed weak global aromaticity. Moreover, as-ID exhibited a larger diradical character and a smaller singlet-triplet gap than s-ID. All the differences can be attributed to their distinct quinoidal substructures.

Mechanochemical preparation of strongly emissive monosubstituted triarylphosphane gold(i) compounds activated by hydrogen bonding driven aggregations

Gold(i) triarylphosphane compounds are a well-known class of coordination compounds displaying from mild to strong emissive properties. Mechanochemical approaches to the preparation, spectroscopic characterization, X-ray diffraction structural determination, and photophysical studies of green emissive neutral linear monophosphane or neutral pseudo-T-shaped or cationic bis-phosphane gold(i) compounds, are herein discussed. The mechanochemical approach to the preparation of gold(i) derivatives was particularly successful for ligands bearing the carboxylic group, while the preparation with esterified ligands yields better results with solvent-mediated methods. The introduction of carboxyl or ester substituents in one aryl group favors the ligand-centered emissions. The analysis of the origin of the emissions was elucidated on the basis of DFT calculations, addressing the emissive behavior to ligand-centered excited states, strongly affected by supramolecular reversible hydrogen bonding aggregation. The study indicates that the ligand with the carboxylic group is particularly suitable for the mechanochemical preparation of emissive gold(i) complexes for material science applications.

Activating Hydrogen Evolution Reaction on Carbon Nanotube via Aryl Functionalisation: The Role of Hybrid sp2-sp3 Interface and Curvature

The hydrogen evolution reaction (HER) is a remarkable mechanism which yields the production of hydrogen through a process of water electrolysis. However, the evolution of hydrogen requires highly conductive and stable catalysts, such as the noble metal platinum (Pt). However, the problem lies in the limitations that this catalyst and others of its kind present. Due to limited availability, as well as the costs involved in acquiring such catalysts, researchers are challenged to manufacture catalysts that do not present these limitations. Carbon nanotubes (CNTs), which are nanomaterials, are known to have a wide range of applications. However, specifically, the pristine carbon nanotube is not suitable for the HER due to the binding free energy of its positive H-atoms. Hence, for the first time, we demonstrated the use of the proposed aryl-functionalised catalysts, i.e., Aryl-L@SWCNT (L = Br, CCH, Cl, CO₂CH₃, F, I, NO₂, or t-butyl), along with the effect of the sp2-sp3 hybridised interface through the density functional theory (DFT). We performed calculations of single-walled carbon nanotubes with multiple aryl functional groups. By employing the DFT calculations, we proved that the curvature of the nanotubes along with the proposed aryl-functionalised catalysts had a noteworthy effect on the performance of the HER. Our study opens the door to investigating a promising group of catalysts for sustainable hydrogen production.

An insight into the mixed quantum mechanical-molecular dynamic simulation of a ZnII-Curcumin complex with a chosen DNA sequence that supports experimental DNA binding investigations

An important aspect of research pertaining to Curcumin (HCur) is the need to arrest its degradation in aqueous solution and in biological milieu. This may be achieved through complex formation with metal ions. For this reason, a complex of HCur was prepared with Zn^II, that is not likely to be active in redox pathways, minimizing further complications. The complex is monomeric, tetrahedral, with one HCur, an acetate and a molecule of water bound to Zn^II. It arrests degradation of HCur to a considerable extent that was realized by taking it in phosphate buffer and in biological milieu. The structure was obtained by DFT calculations. Stable adduct formation was identified between optimized structures of HCur and [Zn(Cur)] with DNA (PDB ID: 1BNA) through experiments validated with multiscale modeling approach. Molecular docking studies provide 2D and 3D representations of binding of HCur and [Zn(Cur)] through different non-covalent interactions with the nucleotides of the chosen DNA. Through molecular dynamics simulation, a detailed understanding of binding pattern and key structural characteristics of the generated DNA-complex was obtained following analysis by RMSD, RMSF, radius of gyration, SASA and aspects like formation of hydrogen bonds. Experimental studies provide binding constants for [Zn(Cur)] with calf thymus DNA at 25 °C that effectively helps one to realize its high affinity towards DNA. In the absence of an experimental binding study of HCur with DNA, owing to its tendency to degrade in solution, a theoretical analysis of the binding of HCur to DNA is extremely helpful. Besides, both experimental and simulated binding of [Zn(Cur)] to DNA may be considered as a case of pseudo-binding of HCur to DNA. In a way, such studies on interaction with DNA helps one to identify HCur's affinity for cellular target DNA, not realized through experiments. The entire investigation is an understanding of experimental and theoretical approaches that has been compared continuously, being particularly useful when a molecule's interaction with a biological target cannot realized experimentally.

Molecular Design of Luminescent Complexes of Eu(III): What Can We Learn from the Ligands

The luminescent metal-organic complexes of rare earth metals are advanced materials with wide application potential in chemistry, biology, and medicine. The luminescence of these materials is due to a rare photophysical phenomenon called antenna effect, in which the excited ligand transmits its energy to the emitting levels of the metal. However, despite the attractive photophysical properties and the intriguing from a fundamental point of view antenna effect, the theoretical molecular design of new luminescent metal-organic complexes of rare earth metals is relatively limited. Our computational study aims to contribute in this direction, and we model the excited state properties of four new phenanthroline-based complexes of Eu(III) using the TD-DFT/TDA approach. The general formula of the complexes is EuL₂A₃, where L is a phenanthroline with -2-CH₃O-C₆H₄, -2-HO-C₆H₄, -C₆H₅ or -O-C₆H₅ substituent at position 2 and A is Cl^- or NO₃^-. The antenna effect in all newly proposed complexes is estimated as viable and is expected to possess luminescent properties. The relationship between the electronic properties of the isolated ligands and the luminescent properties of the complexes is explored in detail. Qualitative and quantitative models are derived to interpret the ligand-to-complex relation, and the results are benchmarked with respect to available experimental data. Based on the derived model and common molecular design criteria for efficient antenna ligands, we choose phenanthroline with -O-C₆H₅ substituent to perform complexation with Eu(III) in the presence of NO₃¯. Experimental results for the newly synthesized Eu(III) complex are reported with a luminescent quantum yield of about 24% in acetonitrile. The study demonstrates the potential of low-cost computational models for discovering metal-organic luminescent materials.

Exploring Proton-Coupled Electron Transfer at Multiple Scales

The coupling of electron and proton transfer is critical for chemical and biological processes spanning a wide range of length and time scales and often occurring in complex environments. Thus, diverse modeling strategies, including analytical theories, quantum chemistry, molecular dynamics, and kinetic modeling, are essential for a comprehensive understanding of such proton-coupled electron transfer reactions. Each of these computational methods provides one piece of the puzzle, and all these pieces must be viewed together to produce the full picture.

Photo-induced phase-transitions in complex solids

Photo-induced phase-transitions (PIPTs) driven by highly cooperative interactions are of fundamental interest as they offer a way to tune and control material properties on ultrafast timescales. Due to strong correlations and interactions, complex quantum materials host several fascinating PIPTs such as light-induced charge density waves and ferroelectricity and have become a desirable setting for studying these PIPTs. A central issue in this field is the proper understanding of the underlying mechanisms driving the PIPTs. As these PIPTs are highly nonlinear processes and often involve multiple time and length scales, different theoretical approaches are often needed to understand the underlying mechanisms. In this review, we present a brief overview of PIPTs realized in complex materials, followed by a discussion of the available theoretical methods with selected examples of recent progress in understanding of the nonequilibrium pathways of PIPTs.

Predicting DNA-Reactivity of N-Nitrosamines: A Quantum Chemical Approach

N-Nitrosamines (NAs) are a class of reactive organic chemicals that humans may be exposed to from environmental sources, food but also impurities in pharmaceutical preparations. Some NAs were identified as DNA-reactive mutagens and many of those have been classified as probable human carcinogens. Beyond high-potency mutagenic carcinogens that need to be strictly controlled, NAs of low potency need to be considered for risk assessment as well. NA impurities and nitrosylated products of active pharmaceutical ingredients (APIs) often arise from production processes or degradation. Most NAs require metabolic activation to ultimately become carcinogens, and their activation can be appropriately described by first-principles computational chemistry approaches. To this end, we treat NA-induced DNA alkylation as a series of subsequent association and dissociation reaction steps that can be calculated stringently by density functional theory (DFT), including α-hydroxylation, proton transfer, hydroxyl elimination, direct S_N2/S_NAr DNA alkylation, competing hydrolysis and S_N1 reactions. Both toxification and detoxification reactions are considered. The activation reactions are modeled by DFT at a high level of theory with an appropriate solvent model to compute Gibbs free energies of the reactions (thermodynamical effects) and activation barriers (kinetic effects). We study congeneric series of aliphatic and cyclic NAs to identify trends. Overall, this work reveals detailed insight into mechanisms of activation for NAs, suggesting that individual steric and electronic factors have directing and rate-determining influence on the formation of carbenium ions as the ultimate pro-mutagens and thus carcinogens. Therefore, an individual risk assessment of NAs is suggested, as exemplified for the complex API-like 4-(N-nitroso-N-methyl)aminoantipyrine which is considered as low-potency NA by in silico prediction.

Machine learning the Hohenberg-Kohn map for molecular excited states

The Hohenberg-Kohn theorem of density-functional theory establishes the existence of a bijection between the ground-state electron density and the external potential of a many-body system. This guarantees a one-to-one map from the electron density to all observables of interest including electronic excited-state energies. Time-Dependent Density-Functional Theory (TDDFT) provides one framework to resolve this map; however, the approximations inherent in practical TDDFT calculations, together with their computational expense, motivate finding a cheaper, more direct map for electronic excitations. Here, we show that determining density and energy functionals via machine learning allows the equations of TDDFT to be bypassed. The framework we introduce is used to perform the first excited-state molecular dynamics simulations with a machine-learned functional on malonaldehyde and correctly capture the kinetics of its excited-state intramolecular proton transfer, allowing insight into how mechanical constraints can be used to control the proton transfer reaction in this molecule. This development opens the door to using machine-learned functionals for highly efficient excited-state dynamics simulations.

Nonlinear optical properties and optimization strategies of D-π-A type phenylamine derivatives in the near-infrared region

Modifying simple molecular structures to significantly improve nonlinear optical (NLO) performance is a primary prerequisite for scientific research. Based on the four phenylamine derivatives reported in previous studies, we designed four organic nonlinear molecules by changing the acceptor group and π-linker. (Time-dependent) density functional theory (DFT/TD-DFT) was performed on molecular geometry optimization, the contribution of π electrons to the bond order, linear and two-photon absorption (TPA) spectra, the intra-molecular charge transfer matrix (CTM), and NLO coefficients. These aspects were considered to analyze in detail how the structural modification of acceptors and π-linkers affects NLO characteristics. The three modification methods were: adding a carbonyl group at the junction of the π-linker and the acceptor group, adding a carbonyl group and a nitrogen atom to the acceptor group, and replacing the quinolinone with a pyrenyl group as the π-linker. The latter two methods can significantly reduce the excitation energy and enhance the intensity of intra-molecular charge transfer during the two-photon transition. The maximum TPA cross-sections and wavelengths of the designed molecules are DPPM (84722.6 GM, 815.7 nm) and DDPM (21600.6 GM, 781.3 nm). These two molecules have large TPA cross-sections in the near-infrared region, which renders them as possible NLO materials with broad application prospects.

Optimization of density fitting auxiliary Slater-type basis functions for time-dependent density functional theory

A new set of auxiliary basis function suitable to fit the induced electron density is presented. Such set has been optimized in order to furnish accurate absorption spectra using the complex polarizability algorithm of time-dependent density functional theory (TDDFT). An automatic procedure has been set up, able, thanks to the definition of suitable descriptors, to evaluate the resemblance of the auxiliary basis-dependent calculated spectra with respect to a reference. In this way, it has been possible to reduce the size of the basis set maximizing the basis set accuracy. Thanks to the choice to employ a collection of molecules for each element, such basis has proven transferable to molecules outside the collection. The final sets are therefore much more accurate and smaller than the previously optimized ones and have been already included in the database of the last release of the AMS suite of programs. The availability of the present new set will allow to improve drastically the applicability range of the polTDDFT method with higher accuracy and less computational effort.

Time-dependent density functional theory calculations of electronic friction in non-homogeneous media

The excitation of low-energy electron-hole pairs is one of the most relevant processes in the gas-surface interaction. An efficient tool to account for these excitations in simulations of atomic and molecular dynamics at surfaces is the so-called local density friction approximation (LDFA). The LDFA is based on a strong approximation that simplifies the dynamics of the electronic system: a local friction coefficient is defined using the value of the electronic density for the unperturbed system at each point of the dynamics. In this work, we apply real-time time-dependent density functional theory to the problem of the electronic friction of a negative point charge colliding with spherical jellium metal clusters. Our non-adiabatic, parameter-free results provide a benchmark for the widely used LDFA approximation and allow the discussion of various processes relevant to the electronic response of the system in the presence of the projectile.

Coupling between Plasmonic and Molecular Excitations: TDDFT Investigation of an Ag-Nanorod/BODIPY-Dye Interaction

A time-dependent density functional theory (TDDFT) computational approach is employed to study the optical coupling between a plasmonic system (a Ag₅₀ nanorod) and a fluorescent dye (BODIPY). It is found that the BODIPY dye can interact with a plasmonic system in a rather different and selective way according to the mutual orientation of the fragments. Indeed, (i) the plasmon excitation turns out to be sensitive to the presence of the BODIPY transition and (ii) this can lead to amplify or suppress the resonance accordingly to the relative orientation of the corresponding transition dipoles. To understand the coupling mechanism, we analyze the shape of the induced density in real space and the Individual Component Map of the Oscillator Strength (ICM-OS) plots and achieve a simple rationalization and insight on the origin and features of the coupling. The resulting possibility of understanding plasmon/fluorophore interactions by simple qualitative arguments opens the way to a rational design of hybrid (plasmon + dye) systems with the desired optical behavior.

Going beyond the Electric-Dipole Approximation in the Calculation of Absorption and (Magnetic) Circular Dichroism Spectra including Scalar Relativistic and Spin-Orbit Coupling Effects

In this work, a time-dependent density functional theory (TD-DFT) scheme for computing optical spectroscopic properties in the framework of linearly and circularly polarized light is presented. The scheme is based on a previously formulated theory for predicting optical absorption and magnetic circular dichroism (MCD) spectra. The scheme operates in the framework of the full semi-classical field-matter interaction operator, thus generating a powerful and general computational scheme capable of computing the absorption (ABS), circular dichroism (CD), and MCD spectra. In addition, our implementation includes the treatment of relativistic effects in the framework of quasidegenerate perturbation theory, which accounts for scalar relativistic effects (in the self-consistent field step) and spin-orbit coupling (in the TD-DFT step), as well as external magnetic field perturbations. Hence, this formalism is also able to probe spin-forbidden transitions. The random orientations of molecules are taken into account by a semi-numerical approach involving a Lebedev numerical quadrature alongside analytical integration. It is demonstrated the numerical quadrature requires as few as 14 points for satisfactory converged results thus leading to a highly efficient scheme, while the calculation of the exact transition moments creates no computational bottlenecks. It is demonstrated that at zero magnetic field, the CD spectrum is recovered while the sum of left and right circularly polarized light contributions provides the linear absorption spectrum. The virtues of this efficient and general protocol are demonstrated on a selected set of organic molecules where the various contributions to the spectral intensities have been analyzed in detail.

Azo-oximate metal-carbonyl to metallocarboxylic acid via the intermediate Ir(III) radical congener: quest for co-ligand driven stability of open- and closed-shell complexes

The redox non-innocent behavior of the diaryl-azo-oxime ligand L^NOH1 has been accentuated via the synthesis of metastable anion radical complexes of type trans-[Ir(L^NO˙-)Cl(CO)(PPh₃)₂] 2 (CO is trans to azo group of the ligand) by the oxidative coordination reaction of 1 with Vaska's complex. The stereochemical role of co-ligands vis-à-vis the interplay of π-bonding has been found to be decisive in controlling the aptitude of the coordinated redox non-innocent ligand to accept or reject an electron. This has been clarified via the isolation of quite a few complexes as well as the failure to synthesize some others. The oxidized analogues of type trans-[Ir(L^NO-)Cl(CO)(PPh₃)₂]⁺2+ (CO and azo group of the ligand are trans) as well as its cis isomer cis-[Ir(L^NO-)Cl(CO)(PPh₃)₂]⁺3+ (CO and azo group of the ligand are cis) have been structurally characterized but the radical anion congener of the latter could not be synthesized. Furthermore, the closed shell complexes [Ir(L^NO-)Cl₂(PPh₃)₂] 4 and [Ir(L^NO-)₂Cl(PPh₃)] 5 have been well characterized by diffraction as well as spectral techniques but their corresponding azo anion radical complexes could not be isolated and this is attributed to the trans influence of ancillary ligands. The anion radical complexes trans-[Ir(L^NO˙-)Cl(CO)(PPh₃)₂] 2 may be rapidly transformed to the metallocarboxylic acids trans-[Ir(L^NO-)Cl(CO₂H)(PPh₃)₂] 6via a proton-coupled electron transfer (PCET) process, thereby demonstrating the role of odd electron over the coordinated ligand framework to trigger metal-mediated carbonyl to carboxylic acid functionalization. Complexes 6 are further stabilized via intramolecular -CO₂H⋯ON- (carboxylic acid⋯oximato) H-bonding. The optoelectronic properties as well as the origin of transitions in the complexes were analyzed by TD-DFT and theoretical analysis, which further disclose that the odd electron in trans-[Ir(L^NO˙-)Cl(CO)(PPh₃)₂] 2 is primarily azo-oxime centric with very low contribution from the iridium center.

Recommendation of Orbitals for G0W0 Calculations on Molecules and Crystals

The many-body GW approximation, especially the G₀W₀ method, has been widely used for condensed matter and molecules to calculate quasiparticle energies for ionization, electron attachment, and band gaps. Because G₀W₀ calculations are well-known to have a strong dependence on the orbitals, the goal of the present work is to provide guidance on the choice of density functional used to generate orbitals and to recommend a choice that gives the most broadly accurate results. We have systematically investigated the dependence of G₀W₀ calculations on the orbitals for 100 molecules and 8 crystals by considering orbitals obtained with a diverse set of Kohn-Sham (KS) and generalized KS (GKS) functionals (63 functionals plus Hartree-Fock). The percentage of Hartree-Fock exchange employed in density functionals has been found to have strong influence on the predicted molecular ionization energy and crystal fundamental band gaps (with optimum values between 40 and 56%), but to have less effect on predicting molecular electron affinities. The low cost of the Gaussian implementation, even with hybrid functionals in periodic calculations, the better performance of global hybrids as compared to range-separated hybrids of either than screened exchange or long-range-corrected type, and the relatively low cost of global-hybrid-functional periodic calculations using Gaussians means that one can employ global-hybrid functionals at a very reasonable cost and obtain more accurate band gaps of semiconductors than are obtained by the methods currently widely employed, namely local gradient approximations. We single out three global-hybrid functionals that give especially good results for both molecules (100 in the test set) and crystals (8 in the test set, for all of which our benchmark data are the proper band gap rather than an optical band gap uncorrected for exciton effects).

Noncovalently bound excited-state dimers: a perspective on current time-dependent density functional theory approaches applied to aromatic excimer models

Excimers are supramolecular systems whose binding strength is influenced by many factors that are ongoing challenges for computational methods, such as charge transfer, exciton coupling, and London dispersion interactions. Treating the various intricacies of excimer binding at an adequate level is expected to be particularly challenging for Time-Dependent Density Functional Theory (TD-DFT) methods. In addition to well-known limitations for some TD-DFT methods in the description of charge transfer or exciton coupling, the inherent London dispersion problem from ground-state DFT translates to TD-DFT. While techniques to appropriately treat dispersion in DFT are well-developed for electronic ground states, these dispersion corrections remain largely untested for excited states. Herein, we aim to shed light on current TD-DFT methods, including some of the newest developments. The binding of four model excimers is studied across nine density functionals with and without the application of additive dispersion corrections against a wave function reference of SCS-CC2/CBS(3,4) quality, which approximates select CCSDR(3)/CBS data adequately. To our knowledge, this is the first study that presents single-reference wave function dissociation curves at the complete basis set level for the assessed model systems. It is also the first time range-separated double-hybrid density functionals are applied to excimers. In fact, those functionals turn out to be the most promising for the description of excimer binding followed by global double hybrids. Range-separated and global hybrids-particularly with large fractions of Fock exchange-are outperformed by double hybrids and yield worse dissociation energies and inter-molecular equilibrium distances. The deviation between each assessed functional and reference increases with system size, most likely due to missing dispersion interactions. Additive dispersion corrections of the DFT-D3(BJ) and DFT-D4 types reduce the average errors for TD-DFT methods but do so inconsistently and therefore do not offer a black-box solution in their ground-state parametrised form. The lack of appropriate description of dispersion effects for TD-DFT methods is likely hindering the practical application of the herein identified more efficient methods. Dispersion corrections parametrised for excited states appear to be an important next step to improve the applicability of TD-DFT methods and we hope that our work assists with the future development of such corrections.

Cation-anion interaction effect on the nonlinear optical behavior of pyridinium-based ionic liquids

Influence of the cation-anion interaction (expressed as internal hydrogen bonds) on the nonlinear optical properties (dipole moment, polarizability, anisotropy of polarizability and first-order hyperpolarizability) for the ionic liquids pyridinium hydrogen sulfate (PHS) and pyridinium bis(dihydrogen phosphate) (P2HP) was investigated in this paper. It was achieved analyzing vibrational modes (FT-IR) and electronic transitions (UV-vis) in PHS and P2HP. Molecular geometry of the title compounds was optimized by means of density functional theory calculations at B3LYP/6-311++G(2d,2p) level. Structural analysis exposed stronger internal hydrogen bonds (H⁷···O³ and H¹⁸···O⁴) between cationic and anionic fragments in the pyridinium hydrogen sulfate than these (H²⁶···O⁴ and H²⁵···O¹²) in pyridinium bis(dihydrogen phosphate). These observations are confirmed by FT-IR and UV-vis methods, where: (i) the infrared bands due to H···O and N-C interactions in P2HP appear at lower frequencies than the same ones in the spectrum of PHS and (ii) the n → π* electron transitions in PHS are more pronounced in comparison with these in P2HP. Based on the weaker cation-anion interaction in P2HP with respect to PHS, the former one is characterized with a remarkably higher first-order hyperpolarizability value. Hence, the ionic liquid P2HP was referred as a better nonlinear optical material than PHS.

Functionalized Crystalline N-Trimethyltriindoles: Counterintuitive Influence of Peripheral Substituents on Their Semiconducting Properties

Three crystalline N-trimethyltriindoles endowed with different functionalities at 3, 8 and 13 positions (either unsubstituted or with three methoxy or three acetyl groups attached) are investigated, and clear correlations between the electronic nature of the substituents and their solid-state organization, electronic properties and semiconductor behavior are established. The three compounds give rise to similar columnar hexagonal crystalline structures; however, the insertion of electron-donor methoxy groups results in slightly shorter stacking distances when compared with the unsubstituted derivative, whereas the insertion of electron-withdrawing acetyl groups lowers the crystallinity of the system. Functionalization significantly affects hole mobilities with the triacetyl derivative showing the lowest mobility within the series in agreement with the lower degree of order. However, attaching three methoxy groups also results in lower hole mobility values in the OFETs (0.022 vs. 0.0014 cm² V^-1 s^-1) in spite of the shorter stacking distances. This counterintuitive behavior has been explained with the help of DFT calculations performed to rationalize the interplay between the intramolecular and intermolecular properties, which point to lower transfer integrals in the trimethoxy derivative due to the HOMO wave function extension over the peripheral methoxy groups. The results of this study provide useful insights into how peripheral substituents influence the fundamental charge transport parameters of chemically modified triindole platforms of fundamental importance to design new derivatives with improved semiconducting performance.

Theoretical Study of Light-Induced Crosslinking Reaction Between Pyrimidine DNA Bases and Aromatic Amino Acids

Low-lying electronic excited states and their relaxation pathways as well as energetics of the crosslinking reaction between uracil as a model system for pyrimidine-type building blocks of DNA and RNA and benzene as a model system for aromatic groups of tyrosine (Tyr) and phenylalanine (Phe) amino acids have been studied in the framework of density functional theory. The equilibrium geometries of the ground and electronic excited states as well as the crossing points between the potential energy surfaces of the uracil–benzene complex were computed. Based on these results, different relaxation pathways of the electronic excited states that lead to either back to the initial geometry configuration or the dimerization between the six-membered rings of the uracil–benzene complex have been identified, and the energetic conditions for their occurrence are discussed. It can be concluded that the DNA–protein crosslinking reaction can be induced by the external electromagnetic field via the dimerization reaction between the six-membered rings of the uracil–benzene pair at the electronic excited-state level of the complex. In the case of the uracil–phenol complex, the configuration of the cyclic adduct (dimerized) conformation is less likely to be formed.

Fragment-Based Quantum Mechanical Calculation of Excited-State Properties of Fluorescent RNAs

Fluorescent RNA aptamers have been successfully applied to track and tag RNA in a biological system. However, it is still challenging to predict the excited-state properties of the RNA aptamer–fluorophore complex with the traditional electronic structure methods due to expensive computational costs. In this study, an accurate and efficient fragmentation quantum mechanical (QM) approach of the electrostatically embedded generalized molecular fractionation with conjugate caps (EE-GMFCC) scheme was applied for calculations of excited-state properties of the RNA aptamer–fluorophore complex. In this method, the excited-state properties were first calculated with one-body fragment quantum mechanics/molecular mechanics (QM/MM) calculation (the excited-state properties of the fluorophore) and then corrected with a series of two-body fragment QM calculations for accounting for the QM effects from the RNA on the excited-state properties of the fluorophore. The performance of the EE-GMFCC on prediction of the absolute excitation energies, the corresponding transition electric dipole moment (TEDM), and atomic forces at both the TD-HF and TD-DFT levels was tested using the Mango-II RNA aptamer system as a model system. The results demonstrate that the calculated excited-state properties by EE-GMFCC are in excellent agreement with the traditional full-system time-dependent ab initio calculations. Moreover, the EE-GMFCC method is capable of providing an accurate prediction of the relative conformational excited-state energies for different configurations of the Mango-II RNA aptamer system extracted from the molecular dynamics (MD) simulations. The fragmentation method further provides a straightforward approach to decompose the excitation energy contribution per ribonucleotide around the fluorophore and then reveals the influence of the local chemical environment on the fluorophore. The applications of EE-GMFCC in calculations of excitation energies for other RNA aptamer–fluorophore complexes demonstrate that the EE-GMFCC method is a general approach for accurate and efficient calculations of excited-state properties of fluorescent RNAs.

Zn(ii) complexes with thiazolyl–hydrazones: structure, intermolecular interactions, photophysical properties, computational study and anticancer activity

Title ligands and their symmetrical octahedral complexes are not photoluminescent, contrary to other synthesized asymmetrical complexes. In comparison to the ligands, the complexes showed improved antiproliferative activity and lower toxicity.

A theoretical investigation of uranyl covalency via symmetry-preserving excited state structures

The structures of electronically excited states of uranyl are probed via density-based analysis to deepen understanding of uranium bonding.

Exploration of Photophysical and Nonlinear Properties of Salicylaldehyde-Based Functionalized Materials: A Facile Synthetic and DFT Approach

The current research presents the synthesis of novel salicylaldehyde thiosemicarbazones (1-6) and their spectroscopic characterization employing UV-visible, Fourier transform infrared spectroscopy, and NMR techniques. Experimental results are compared and validated with the results obtained theoretically by employing density functional theory at the M06 level with the 6-311G (d,p) basis set. Further, various parameters [natural bond orbital (NBO)], linear and nonlinear optical (NLO) properties, and global reactivity parameters (GRPs) are computationally calculated. The NBO approach has confirmed the stability of compounds on account of charge delocalization and hyper conjugative interaction network. Frontier molecular orbital analysis has explained the charge transfer and chemical reactivity capability, while GRPs have led to the analysis of kinetic stability of the studied molecules. Further, the probability of being NLO-active has been theoretically proved by the HOMO/LUMO energy difference (4.133-4.186 eV) and β values (192.778-501.709 a.u). These findings suggest that the studied compounds possess potential NLO applications as they have shown larger NLO values in comparison with that of the urea molecule, and such distinct properties prove their technological importance.

Simulation of Low-Lying Singlet and Triplet Excited States of Multiple-Resonance-Type Thermally Activated Delayed Fluorescence Emitters by Delta Self-Consistent Field (ΔSCF) Method

The delta self-consistent field (ΔSCF) method was applied to the simulation of low-lying singlet and triplet excited states of multiple-resonance (MR)-type thermally activated delayed fluorescence (TADF) molecules, which form a promising group for organic light-emitting diode (OLED) emitters. A comparison with the experimental values of 13 emitters from the literature showed that ΔSCF gave fairly accurate S₁ and T₁ excitation energies (mean absolute errors (MAEs) of 0.092 and 0.055 eV, respectively) as well as quite accurate ΔE_ST (S₁-T₁ gap, MAE of 0.041 eV), which could not be calculated with sufficient accuracy by the conventional time-dependent density functional theory (TDDFT). ΔSCF also demonstrated its utility for the analysis of photophysical properties through a simulation of the reverse intersystem crossing (RISC) process of an MR-type emitter.

Synthesis and Spectral Characterisation of (E)-3-(3-(4 (Dimethylamino)Phenyl)Acrylo-yl)-4-Hydroxy-2H-Chromen-2-One and Their Antibacterial Activity and Acetylcholinesterase Inhibition

A new coumarin derivative, (E)-3-(3-(4-(dimethylamino) phenyl) acrylo-yl)-4-hydroxy-2H-chromen-2-one (3), was synthesized by the condensation of 3-acetyl-4-hydroxycoumarin (1) with 4-N,N-dimethylaminobenzaldehyde (2) in the presence of piperidine in ethanol. The structure of the synthesized compound was characterized using spectroscopic data (IR and 1H NMR) and elemental analysis. The antimicrobial properties and acetylcholinesterase inhibition activity (AChEI) of coumarin 3 were investigated, with the highest observed AChEI activity providing 48.25% inhibition. The electronic absorption and emission spectra revealed that 3 exists as two, main keto-enol tautomers. The ratios of these tautomers in both protic and aprotic solvents with different polarities and dielectric constants were calculated. The fluorescence of coumarin 3 was enhanced upon increasing the medium viscosity, which was due to the resultant molecular rigidity. This criterion was further investigated using DNA, whereby 3 showed enhanced fluorescence upon its uptake in DNA grooves and was therefore tested as a novel DNA fluorescent stain.

Application of TD-DFT Theory to Studying Porphyrinoid-Based Photosensitizers for Photodynamic Therapy: A Review

An important focus for innovation in photodynamic therapy (PDT) is theoretical investigations. They employ mostly methods based on Time-Dependent Density Functional Theory (TD-DFT) to study the photochemical properties of photosensitizers. In the current article we review the existing state-of-the-art TD-DFT methods (and beyond) which are employed to study the properties of porphyrinoid-based systems. The review is organized in such a way that each paragraph is devoted to a separate aspect of the PDT mechanism, e.g., correct prediction of the absorption spectra, determination of the singlet-triplet intersystem crossing, and interaction with molecular oxygen. Aspects of the calculation schemes are discussed, such as the choice of the most suitable functional and inclusion of a solvent. Finally, quantitative structure-activity relationship (QSAR) methods used to explore the photochemistry of porphyrinoid-based systems are discussed.

cis-[(η5-C5H5)Fe(η1-CO)(μ-CO)]2, the poor relative between cis and trans tautomers. A theoretical study of the gas-phase Fe L3-edge and C and O K-edge XAS of trans-/cis-[(η5-C5H5)Fe(η1-CO)(μ-CO)]2

The relative stability of trans-[(η⁵-C₅H₅)Fe(η¹-CO)(μ-CO)]₂ (trans-I) and cis-I tautomers in a vacuum and in solvents with different dielectric constants (ε) has been investigated by exploiting density functional theory (DFT). Theoretical results indicate that, in agreement with experimental evidence, trans-I is more stable than cis-I in a vacuum (∼1.5 kcal mol^-1; ε = 1), while the opposite is true in media with ε > 7. Differently from solution, DFT outcomes pertaining to the vapor-phase cis-I ⇆ trans-I equilibrium at T = 368 K, the temperature at which the Fe L_2,3-edges and the C and O K-edge X-ray absorption spectroscopy (XAS) data of I have been recorded, ultimately indicate the trans-I predominance (∼93%). Compositions, oscillator strengths (f) and excitation energy (EE) values of cis-I transitions substantially mirror those of trans-I; nevertheless, the weighted cis-If(EE) distributions negligibly contribute to the diverse simulated XA spectra of I.

Electronic Exciton-Plasmon Coupling in a Nanocavity Beyond the Electromagnetic Interaction Picture

The optical response of a system formed by a quantum emitter and a plasmonic gap nanoantenna is theoretically addressed within the frameworks of classical electrodynamics and the time-dependent density functional theory (TDDFT). A fully quantum many-body description of the electron dynamics within TDDFT allows for analyzing the effect of electronic coupling between the emitter and the nanoantenna, usually ignored in classical descriptions of the optical response. We show that the hybridization between the electronic states of the quantum emitter and those of the metallic nanoparticles strongly modifies the energy, the width, and the very existence of the optical resonances of the coupled system. We thus conclude that the application of a quantum many-body treatment that correctly addresses charge-transfer processes between the emitter and the nanoantenna is crucial to address complex electronic processes involving plasmon-exciton interactions directly impacting optoelectronic applications.