Behind the glow: unveiling the nature of NanoLuc reactants and products.

Due to the largely recognized utility of bioluminescence in many fields, a wide variety of luciferase-luciferin systems have been investigated in order to find the best-suited for a number of different applications. The collected knowledge has allowed the identification of a few necessary, or at least desirable, properties, such as bright luminescence, low background signal and small dimension of the enzyme that must exhibit structural stability at operating conditions. The NanoLuc-furimazine pair seems to meet all these requirements, but the mechanism of the reaction and the characteristics of the species responsible for the emission remain unknown. The aim of this study is to identify the luminescent product among the possible forms of oxidized furimazine and to understand how the chemical form and structure of the system, before and after the oxidation, are involved into the reaction mechanism and determine emission. To do this, we consider two possible forms of furimazine, the keto and the enol one, and test which of them is the most plausible candidate in the bioluminescence process on the basis of enzyme-substrate interactions from docking calculations. A similar procedure is repeated for three possible forms of the furimamide luminescent product, and their properties in the protein environment are then evaluated via QM/MM calculations. In contrast with previous indications, our simulations well support the involvement of the enol form of furimazine as reagent and point to the zwitterionic forms of furimamide as emissive species.

Optimal-tuning of range-separated density functionals to describe the optical and photophysical properties of rhodamine B dimers.

Setting up computational approaches enabling the correct prediction of the photophysical properties of rhodamine B (RB) derivatives and their aggregates (referred to as dimers here) is of fundamental importance to rationally drive the design of novel systems of applicative relevance, such as artificial light-harvesting nanosystems. Currently, approaches rooted in time-dependent density functional theory (TD-DFT), which are appealing for their relatively low computational cost, nonetheless have limitations in terms of accuracy, especially while considering RB dimeric species. In this work, we investigated the performances of optimally tuned range-separated hybrid functionals for describing the excited states of RB and its H dimer, focusing on dimeric charge-transfer (CT) states. We compared different optimal-tuning (OT) procedures including or not solvent screening. The results show that the properties of CT states, such as their CT extent, brightness and relative energy ordering, are crucially affected by the Hartree-Fock exchange amount dictated by OT, which is in turn driven by the chosen asymptotic behavior. Finally, to understand the dissimilar pictures provided for the CT states by different tuning approaches, we performed an extensive analysis aimed at elucidating how CT states are affected by range-separation parameters. As a result, a simple procedure is finally provided for easily achieving a unique functional to be applied on both monomer and dimers, its current limitations are highlighted and possible perspectives for future development are envisaged.

Study of sterically crowded alkanes: assessment of non-empirical density functionals including double-hybrid cost-effective methods.

We theoretically study the homolytic dissociation reactions of sterically crowded alkanes of increasing size, carrying three different (bulky) substituents such as tert-butyl, adamantane, and [1.1.1]propellanyl, employing a family of parameter-free functionals ranging from semi-local, to hybrid and double-hybrid models. The study is complemented with the interaction between a pair of HC(CH3)3 molecules at repulsive and attractive regions, as an example of a system composed by a pair of weakly bound sterically crowded alkanes. We also assessed the effect of incorporating reliable dispersion corrections (i.e., D4 or NL) for all the functionals assessed, as well as the use of a tailored basis set (DH-SVPD) for non-covalent interactions, which provides the best trade-off between accuracy and computational cost for a seemingly extended applications to branched or crowded systems. Overall, the PBE-QIDH/DH-SVPD and r2SCAN-QIDH/DH-SVPD methods represent an excellent compromise providing relatively low, and thus very competitive, errors at a fraction of the cost of other quantum-chemical methods in use.

Evaluating Noncovalent Interactions in Halogenated Molecules with Double-Hybrid Functionals and a Dedicated Small Basis Set.

We present here an extension of our recently developed PBE-QIDH/DH-SVPD basis set to halogen atoms, with the aim of obtaining, for weakly interacting halogenated molecules, interaction energies close to those provided by a large basis set (def2-TZVPP) coupled to empirical dispersion potential. The core of our approach is the split-valence basis set, DH-SVPD, that has been developed for F, Cl, Br, and I atoms using a self-consistent formula, containing only energy terms computed for dimers and the corresponding monomers at the same level of theory. The basis set developed considering four systems, one for each halogen atoms, has been then tested on the X40, X4 × 10 benchmarks as well as on other two, less standard, data sets. Finally, a large system (380 atoms) has been also considered as a "crash" test. Our results show that the simple and nonempirical PBE-QIDH/DH-SVPD approach is able to provide accurate results for interaction energies of all the considered systems and can thus be considered as a cheaper alternative to DH functionals paired with empirical dispersion corrections and a large basis set of triple-ζ quality.

Monitoring Density Redistribution at the Excited State in a Dual Emitting Molecule: An Analysis Based on Real-Time Density Functional Theory and Density Descriptors.

In this work, we computed and analyzed, by means of density-based descriptors, the real-time evolution of both the locally excited (LE) and charge-transfer (CT) excited states for the planar and twisted conformations of the DMABN (4-(N,N-dimethylamino)benzonitrile) molecule using real-time time-dependent density functional theory (DFT) and three different exchange-correlation energy functionals (EXC) belonging to the same family (the PBE one). Our results based on the analysis of density-based descriptors show that the underlying EXC modifies the evolution in time of the density. In particular, comparing the frequency of density reorganization computed with the three functionals (PBE, PBE0, and LC-PBE), we found that the frequency of electronic interconversion of the individual determinants involved during the dynamics increases from PBE to PBE0 and to LC-PBE. This allows us to show that there is a correlation between the delocalization of the electronic density and the frequency of reorganization. In particular, the greater the mean hole-electron distance during the dynamics, the lower is the frequency of density reorganization.

Modeling Multi-Step Organic Reactions: Can Density Functional Theory Deliver Misleading Chemistry?

Many organic reactions are characterized by a complex mechanism with a variety of transition states and intermediates of different chemical natures. Their correct and accurate theoretical characterization critically depends on the accuracy of the computational method used. In this work, we study a complex ambimodal cycloaddition with five transition states, two intermediates, and three products, and we ask whether density functional theory (DFT) can provide a correct description of this type of complex and multifaceted reaction. Our work fills a gap in that most systematic benchmarks of DFT for chemical reactions have considered much simpler reactions. Our results show that many density functionals not only lead to seriously large errors but also differ from one another in predicting whether the reaction is ambimodal. Only a few of the available functionals provide a balanced description of the complex and multifaceted reactions. The parameters varied in the tested functionals are the ingredients, the treatment of medium-range and nonlocal correlation energy, and the inclusion of Hartree-Fock exchange. These results show a clear need for more benchmarks on the mechanisms of large molecules in complex reactions.

Axial-equatorial equilibrium in substituted cyclohexanes: a DFT perspective on a small but complex problem.

In Chemistry, complexity is not necessarily associated to large systems, as illustrated by the textbook example of axial-equatorial equilibrium in mono-substituted cyclohexanes. The difficulty in modelling such a simple isomerization is related to the need for reproducing the delicate balance between two forces, with opposite effects, namely the attractive London dispersion and the repulsive steric interactions. Such balance is a stimulating challenge for density-functional approximations and it is systematically explored here by considering 20 mono-substituted cyclohexanes. In comparison to highly accurate CCSD(T) reference calculations, their axial-equatorial equilibrium is studied with a large set of 48 exchange-correlation approximations, spanning from semilocal to hybrid to more recent double hybrid functionals. This dataset, called SAV20 (as Steric A-values for 20 molecules), allows to highlight the difficulties encountered by common and more original DFT approaches, including those corrected for dispersion with empirical potentials, the 6-31G*-ACP model, and our cost-effective PBE-QIDH/DH-SVPD protocol, in modeling these challenging interactions. Interestingly, the performance of the approaches considered in this contribution on the SAV20 dataset does not correlate with that obtained with other more standard datasets, such as S66, IDISP or NC15, thus indicating that SAV20 covers physicochemical features not already considered in previous noncovalent interaction benchmarks.

SOS1-RSX-QIDH: A spin-opposite-scaled range-separated-exchange quadratic-integrand double-hybrid density functional.

We develop and validate the SOS1-RSX-QIDH density functional, a one-parameter spin-opposite-scaled variant of the range-separated-exchange quadratic-integrand double-hybrid (RSX-QIDH) model. By entering into the family of spin-biased double hybrids, this new density functional benefits from an improved computational scaling that rivals with the one of hybrids, still conserving the accuracy of its RSX-QIDH version. As part of the latter family, this density functional is well-adapted to treat molecular systems that are particularly prone to self-interaction errors in their ground and excited states. In particular, we show that the SOS1-RSX-QIDH model is a good compromise to treat ground-state problems dealing with kinetics and has a real added value when applied to the evaluation of the excited-state properties of equilibrium and out-of-equilibrium molecular complexes. Even if spin-biased double hybrids are recognized to strongly underestimate noncovalent interactions, we notice and recommend coupling SOS1-RSX-QIDH with a nonlocal van der Waals potential, a combination that is here proved to compete with the best density-functional approximations currently in use.

Effect of Polymer Composition on the Optical Properties of a New Aggregation-Induced Emission Fluorophore: A Combined Experimental and Computational Approach.

Nowadays, fluorophores with a tetraphenylethylene (TPE) core are considered interesting due to the aggregation-induced emission (AIE) behavior that enables their effective use in polymer films. We propose a novel TPE fluorophore (TPE-BPAN) bearing two dimethylamino push and a 4-biphenylacetonitrile pull moieties with the typical AIE characteristics in solution and in the solid state, as rationalized by DFT calculations. Five different host polymer matrices with different polarity have been selected: two homopolymers of poly(methylmethacrylate) (PMMA) and poly(cyclohexyl methacrylate) (PCHMA) and three copolymers at different compositions (P(MMA-co-CHMA) 75:25, 50:50, and 25:75 mol%). The less polar comonomer of CHMA appeared to enhance TPE-BPAN emission with the highest quantum yield (QY) of about 40% measured in P(MMA-co-CHMA) 75:25. Further reduction in polymer polarity lowered QY and decreased the film stability and adhesion to the glass surface. LSC performances were not significantly affected by the matrix's polarity and resulted in around one-third of the state-of-the-art due to the reduced QY of TPE-BPAN. The theoretical investigation based on density functional theory (DFT) calculations clarified the origin of the observed AIE and the role played by the environment in modulating the photophysical behavior.

Concerted versus stepwise proton transfer reactions in the [2, 2'-bipyridyl]-3-3'-diol molecule: A static and dynamic ab-initio investigation.

The double proton transfer (PT) reaction has been investigated in the [2,2'-bipyridyl]-3-3'-diol, a complex molecule where the proton movements is coupled to significant rearrangement of the electronic structure. Moreover, the reaction could be concerted, that is the two protons are exchanged simultaneously, or stepwise, where the two protons are transferred sequentially. To this end, a static exploration of the potential energy surface (PES) was carried together with the analysis of the free-energy surface (FES), both surfaces being evaluated at density functional theory level and different exchange-correlation functionals. While the concerted mechanism has been clearly discharged, the characteristics of the stepwise PT significantly depends on the chosen functionals, some suggesting a clear stepwise mechanism characterized by a stable reaction intermediates and two transitions states, whereas other approaches propend for a asynchronous PT, with a single TS. These features appear on both PES and FES, albeit some differences appears due to their different nature.

Barrier Heights for Diels-Alder Transition States Leading to Pentacyclic Adducts: A Benchmark Study of Crowded, Strained Transition States of Large Molecules.

Theoretical characterization of reactions of complex molecules depends on providing consistent accuracy for the relative energies of intermediates and transition states. Here we employ the DLPNO-CCSD(T) method with core-valence correlation, large basis sets, and extrapolation to the CBS limit to provide benchmark values for Diels-Alder transition states leading to competitive strained pentacyclic adducts. We then used those benchmarks to test a diverse set of wave function and density functional methods for the absolute and relative barrier heights of these transition states. Our results show that only a few of the tested density functionals can predict the absolute barrier heights satisfactorily, although relative barrier heights are more accurate. The most accurate functionals tested are ωB97M-V, M11plus, ωB97X-V, PBE-D3(0), M11, and MN15 with MUDs from best estimates less than 3.0 kcal. These findings can guide selection of density functionals for future studies of crowded, strained transition states of large molecules.

Organic compounds for solid state luminescence enhancement/aggregation induced emission: a theoretical perspective.

Organic luminophores displaying one or more forms of luminescence enhancement in solid state are extremely promising for the development and performance optimization of functional materials essential to many modern key technologies. Yet, the effort to harness their huge potential is riddled with hurdles that ultimately come down to a limited understanding of the interactions that result in the diverse molecular environments responsible for the macroscopic response. In this context, the benefits of a theoretical framework able to provide mechanistic explanations to observations, supported by quantitative predictions of the phenomenon, are rather apparent. In this perspective, we review some of the established facts and recent developments about the current theoretical understanding of solid-state luminescence enhancement (SLE) with an accent on aggregation-induced emission (AIE). A description of the macroscopic phenomenon and the questions it raises is accompanied by a discussion of the approaches and quantum chemistry methods that are more apt to model these molecular systems with the inclusion of an accurate yet efficient simulation of the local environment. A sketch of a general framework, building from the current available knowledge, is then attempted via the analysis of a few varied SLE/AIE molecular systems from literature. A number of fundamental elements are identified offering the basis for outlining design rules for molecular architectures exhibiting SLE that involve specific structural features with the double role of modulating the optical response of the luminophores and defining the environment they experience in solid state.

Tuning the Photophysical Properties of Ru(II) Photosensitizers for PDT by Protonation and Metallation: A DFT Study.

Due to their activity, photosensitizers with the Ru(II)-polypyridyl complex structure represent an intriguing class of photodynamic therapy agents used to treat neoplasms. However, their solubility is poor, intensifying experimental research into improving this property. One recently proposed solution is to attach a polyamine macrocycle ring. In this paper, the density functional theory (DFT) and time-dependent DFT (TD-DFT) studies on such derivative were performed to assess the impact of the protonation-capable macrocycle and its ability to chelate transition state metals, as illustrated by the Cu(II) ion, on the expected photophysical activity. These properties were determined by examining ultraviolet-visible (UV-vis) spectra, intersystem conversion, and type I and II photoreactions of all species possibly present in a tumor cell. For comparison purposes the structure devoid of the macrocycle was also examined. The results show that the subsequent protonation of amine groups improves the reactivity, with [H2L]4+/[H3L]5+ being borderline, whereas complexation appears to weaken the desired photoactivity.

Absorption Spectra of Flexible Fluorescent Probes by a Combined Computational Approach: Molecular Dynamics Simulations and Time-Dependent Density Functional Theory.

A detailed understanding and interpretation of absorption spectra of molecular systems, especially in condensed phases, requires computational models that allow their structural and electronic features to be connected to the observed macroscopic spectra. This work is focused on modeling the electronic absorption spectrum of a fluorescent probe, namely, the 9-(4-((bis(2-((2-(ethylthio)ethyl)thio)ethyl)amino)methyl)phenyl)-6-(pyrrolidin-1-yl)-3H-xanthen-3-one molecule, depicted by a combined classical-quantum chemical approach. Particularly, first classical molecular dynamics (MD) has been used to explore the configurational space, and next, the absorption spectrum has been reconstructed by averaging the results of time-dependent density functional theory (TD-DFT) calculations performed on equispaced molecular conformations extracted from MD to properly sample the configurational space explored at finite temperature. To verify the effect of molecular conformation on the spectral profile, the generated electronic absorption spectra were compared with those obtained considering a single structure corresponding to the optimized one, an approach also referred to as static. This comparison allows one to highlight a sizable though small shift between the maxima of the corresponding reconstructed absorption spectra, highlighting the importance of conformational sampling in the case of this rather flexible molecule. Four different exchange and correlation functionals (PBE, BLYP, PBE0, B3LYP) were considered to compute vertical transition via TD-DFT calculations. From the results obtained in gas and in condensed, here solution, phases, it appears that the magnitude of the shift is actually more affected by the phase in which the system is found than by the functional used. This fact underlines the central importance of conformational mobility, that is flexibility, of this molecule. From a more quantitative point of view, a comparison with available experimental data shows that hybrid functionals, such as PBE0 and B3LYP, enable one to faithfully reproduce the observed absorption maxima.

DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science.

In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.

Tackling an accurate description of molecular reactivity with double-hybrid density functionals.

In this Communication, we assess a panel of 18 double-hybrid density functionals for the modeling of the thermochemical and kinetic properties of an extended dataset of 449 organic chemistry reactions belonging to the BH9 database. We show that most of DHs provide a statistically robust performance to model barrier height and reaction energies in reaching the "chemical accuracy." In particular, we show that nonempirical DHs, such as PBE0-DH and PBE-QIDH, or minimally parameterized alternatives, such as ωB2PLYP and B2K-PLYP, succeed to accurately model both properties in a balanced fashion. We demonstrate, however, that parameterized approaches, such as ωB97X-2 or DSD-like DHs, are more biased to only one of both properties.

Double Hybrids and Noncovalent Interactions: How Far Can We Go?

The accurate evaluation of weak noncovalent interactions in large, that is those containing up to thousand atoms, molecular systems represents a difficult challenge for any quantum chemical method. Indeed, some approximations are often introduced to render affordable these calculations. Here, we consider the PBE-QIDH/DH-SVPD protocol, combining a nonempirical double hybrid functional (PBE-QIDH) with a small basis set (DH-SVPD) tailored for noncovalent interactions with a double aim: (i) explore the robustness and accuracy of this protocol with respect to other Density Functional Approximations; (ii) illustrate how its performances are affected by the computational parameters underlying the calculation of the exact exchange and the Coulomb contribution, as well as the perturbative term. To this end, we consider three data sets, namely S66, L7, and CiM13, incorporating molecules of increasing size. On the bright side, our results suggest that the PBE-QIDH/DH-SVPD protocol is particularly accurate for large systems such as those contained in the CiM13 set (up to more than 1000 atoms and 14 000 basis functions), for which the DLPNO approximation leads to a significant speed-up for the evaluation of the perturbative correlation term. However, our analysis also points out the limit of this statistical exercise, when the quality of the reference data cannot be easily assessed, due to the size of the molecular complexes involved, and when the number of molecules is limited.

Free Energy Profiles of Proton Transfer Reactions: Density Functional Benchmark from Biased Ab Initio Dynamics.

By coupling an enhanced sampling algorithm with an orbital-localized variant of Car-Parrinello molecular dynamics, the so-called atomic centered density matrix propagation model, we reconstruct the free energy profiles along reaction pathways using different density functional approximations (DFAs) ranging from locals to hybrids. In particular, we compare the computed free energy barrier height of proton transfer (PT) reactions to those obtained by a more traditional static approach, based on the intrinsic reaction coordinate (IRC), for two case systems, namely malonaldehyde and formic acid dimer. The obtained results show that both the IRC profiles and the potentials of mean force, derived from biased dynamic trajectories, are very sensitive to the density functional approximation applied. More precisely, we observe that, with the notable exception of M06-L, local density functionals always strongly underestimate the reaction barrier heights. More generally, we find that also the shape of the free energy profile is very sensitive to the density functional choice, thus highlighting the effect, often neglected, that the choice of DFA has also in the case of dynamics simulations.

Electronic Energy and Local Property Errors at QTAIM Critical Points while Climbing Perdew's Ladder of Density-Functional Approximations.

We investigate the relationships between electron-density and electronic-energy errors produced by modern exchange-correlation density-functional approximations belonging to all of the rungs of Perdew's ladder. To this aim, a panel of relevant (semi)local properties evaluated at critical points of the electron-density field (as defined within the framework of Bader's atoms-in-molecules theory) are computed on a large selection of molecular systems involved in thermodynamic, kinetic, and noncovalent interaction chemical databases using density functionals developed in a nonempirical and minimally and highly parametrized fashion. The comparison of their density- and energy-based performance, also discussed in terms of density-driven errors, casts light on the strengths and weaknesses of the most recent and efficient density-functional approximations.