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.

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.

Assessment of the nonempirical r2SCAN-QIDH double-hybrid density functional against large and diverse datasets

We update the Quadratic Integrand Double-Hybrid (QIDH) model [J. Chem. Phys. 141, 031101 (2014)] by incorporating the nonempirical restored-regularized Strongly Constrained and Appropriately Normed (r2SCAN) meta-generalized gradient approximation exchange-correlation functional, thus devising a robust density functional approximation free of any empirical parameter and incorporating all the constraints so far known for the exchange-correlation kernel. We assessed the new r2SCAN-QIDH expression on the GMTKN55 database and further extend its application to various types of non-covalent interactions (e.g., S66 × 8, O24 × 5). The assessment done shows that the model becomes very competitive in accuracy with respect to parent exchange-correlation functionals of any type, but without relying on any fitted parameter or numerical training.

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.

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 [H₂L]⁴⁺/[H₃L]⁵⁺ being borderline, whereas complexation appears to weaken the desired photoactivity.

Excitation energies of polycylic aromatic hydrocarbons by double-hybrid functionals: Assessing the PBE0-DH and PBE-QIDH models and their range-separated versions

A family of non-empirical double-hybrid (DH) density functionals, such as Perdew-Burke-Ernzerhof (PBE)0-DH, PBE-QIDH, and their range-separated exchange (RSX) versions RSX-0DH and RSX-QIDH, all using Perdew-Burke-Ernzerhof(PBE) exchange and correlationfunctionals, is applied here to calculate the excitation energies for increasingly longer linear and cyclic acenes as part of their intense benchmarking for excited states of all types. The energies for the two lowest-lying singlet ¹L_a and ¹L_b states of linear oligoacenes as well as the triplet ³L_a and ³L_b states, are calculated and compared with experimental results. These functionals clearly outperform the results obtained from hybrid functionals and favorably compare with other double-hybrid expressions also tested here, such as B2-PLYP, B2GP-PLYP, ωB2-PLYP, and ωB2GP-PLYP. The study is complemented by the computation of adiabatic S₀-T₁ singlet-triplet energy difference for linear acenes as well as the extension of the study to strained cyclic oligomers, showing how the family of non-empirical expressions robustly leads to competitive results.

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.

Stability of the polyynic form of C18, C22, C26, and C30 nanorings: a challenge tackled by range-separated double-hybrid density functionals

We calculate the relative energy between the cumulene and polyyne structures of a set of C_4k+2 (k = 4-7) rings (C₁₈, C₂₂, C₂₆, and C₃₀ prompted by the recent synthesis of the cyclo[18]carbon (or simply C₁₈) compounds. Reference results were obtained by a costly Quantum Monte-Carlo (QMC) approach, providing thus very accurate values allowing to systematically compare the performance of a variety of wavefunction methods [(i.e., MP2, SCS-MP2, SOS-MP2, DLPNO-CCSD, and DLPNO-CCSD(T)] as well as DFT approaches, applying for the latter a diversity of density functionals covering global and range-separated hybrid and double-hybrid models. The influence of the use of a range-separation scheme for density functionals, for both hybrid and double-hybrid expressions, is discussed according to its key role. Overall, range-separated double-hybrid functionals (e.g., RSX-QIDH) behave very accurately and provide competitive results compared with DLPNO-CCSD(T), at a more reasonable computational cost.

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.

Pairing double hybrid functionals with a tailored basis set for an accurate thermochemistry of hydrocarbons

A collection of five challenging datasets, including noncovalent interactions, reaction barriers and electronic rearrangements of medium-sized hydrocarbons, has been selected to verify the robustness of double-hybrid functionals used in conjunction with the small DH-SVPD basis set, especially developed for noncovalent interactions. The analysis is completed by other, more standard functionals, for a total of 17 models, including also empirical corrections for dispersion. The obtained results show that the chemical accuracy threshold, that is an error lower than 1.0 kcal mol⁻¹, can be obtained by pairing the nonempirical PBE-QIDH functional with the DH-SVPD basis set, as well as by other semi-empirical functionals, such as DSD-PBEP86, using larger basis sets and empirical corrections. More in general, a significant improvement can be obtained using the DH-SVPD basis set with DHs, without resorting to any empirical corrections. This choice leads to a fast computational protocol that, avoiding any empirical potential, remains on a fully quantum ground.

The pairing of the PBE-QIDH double-hybrid functional with a tailored split-valence basis set leads to a fast computational protocol for the accurate evaluation of hydrocarbon thermochemistry, without resorting to any empirical correction.

Chasing unphysical TD-DFT excited states in transition metal complexes with a simple diagnostic tool

Transition Metal Complexes (TMCs) are known for the rich variety of their excited states showing different nature and degrees of locality. Describing the energies of these excited states with the same degree of accuracy is still problematic when using time-dependent density functional theory in conjunction with the most current density functional approximations. In particular, the presence of unphysically low lying excited states possessing a relevant Charge Transfer (CT) character may significantly affect the spectra computed at such a level of theory and, more relevantly, the interpretation of their photophysical behavior. In this work, we propose an improved version of the M_AC index, recently proposed by the authors and collaborators, as a simple and computationally inexpensive diagnostic tool that can be used for the detection and correction of the unphysically predicted low lying excited states. The analysis, performed on five prototype TMCs, shows that spurious and ghost states can appear in a wide spectral range and that it is difficult to detect them only on the basis of their CT extent. Indeed, both delocalization of the excited state and CT extent are criteria that must be combined, as in the M_AC index, to detect unphysical states.

Aggregation-Induced Emission: A Challenge for Computational Chemistry Taking TPA-BMO as an Example*

A multi-environment computational approach is proposed to study the modulation of the emission behavior of the triphenylamine (Z)-4-benzylidene-2-methyloxazol-5(4H)-one (TPA-BMO) molecule [Tang et al., J. Phys. Chem. C 119, 21875 (2015)]. We aim at (1) proposing a realistic description of the molecule in several environments (solution, aggregate, polymer matrix), (2) modelling its absorption and emission properties, and (3) providing a qualitative understanding of the experimental observations by highlighting the photophysical phenomena leading to the emission modulation. To this purpose, we rely on (TD-)DFT calculations and classical Molecular Dynamics simulations, but also on the hybrid ONIOM QM/QM' approach and the in situ chemical polymerization methodology. In low-polar solvents, the investigation of the potential energy surfaces and the modulation of the emission quantum yield can be attributed to possible photophysical energy dissipation caused by low-frequency vibrational modes. In the aggregate and in the polymer matrix, the emission modulation can be qualitatively interpreted in terms of the possible restriction of the intramolecular vibrations. For these two systems, our study highlights that a careful modelling of the environment is far from trivial but is fundamental to model the optical properties of the fluorophore.

Beyond Chemical Accuracy for Alkane Thermochemistry: The DHthermo Approach

The so-called protobranching phenomenon, that is the greater stability of branched alkanes with respect to their linear isomers, represents an interesting challenge for approaches based on density functional theory (DFT), since it requires a balanced description of several electronic effects, including (intramolecular) dispersion forces. Here, we investigate this problem using a protocol recently developed based on double-hybrid functionals and a small basis set, DH-SVPD, suited for noncovalent interactions. The energies of bond separation reactions (BSR), defined on the basis of an isodesmic principle, are taken as reference properties for the evaluation of 15 DFT approaches. The obtained results show that error lower than the so-called "chemical accuracy" (<1.0 kcal/mol) can be obtained by the proposed protocol on both relative reaction energies and enthalpies. These results are then verified on the standard BSR36 data set and support the proposition of our computational protocol, named DHthermo, where any DH functional, such as PBE-QIDH or B2PLYP, provides accurate results when coupled to an empirical dispersion correction and the DH-SVPD basis set. This protocol not only gives subchemical accuracy on the thermochemistry of alkanes but it is extremely easy to use with common quantum-chemistry codes.

A combined Monte Carlo/DFT approach to simulate UV-vis spectra of molecules and aggregates: Merocyanine dyes as a case study

The combination of a Monte Carlo (MC) sampling of the configurational space with time dependent-density functional theory (TD-DFT) to estimate vertical excitations energies has been applied to compute the absorption spectra of a family of merocyanine dyes in both their monomeric and dimeric forms. These results have been compared to those obtained using a static DFT/TD-DFT approach as well as to the available experimental spectra. Though suffering of the limitations related to the use of DFT and TD-DFT for this type of systems, our data clearly show that the classical MC sampling provides a suitable alternative to classical molecular dynamics to explore the structural flexibility of these donor-acceptor (D-π-A) chromophores enabling a realistic description of the potential energy surface of both their monomers and aggregates (here dimers) and thus of their spectra. Overall, the combination of MC sampling with quantum mechanics (TD-DFT) calculations, carried out in implicit dioxane solvent on random snapshots, provides a workable compromise to solve the combined challenge of accuracy and time-consuming problem not only for merocyanines momers, but also for their dimers, up to now less investigated. Indeed, the simulated absorption spectra fairly agree with the experimental ones, suggesting the general reliability of the method.

Assessing challenging intra- and inter-molecular charge-transfer excitations energies with double-hybrid density functionals

We investigate the performance of a set of recently introduced range-separated double-hybrid functionals, namely ωB2-PLYP, ωB2GP-PLYP, RSX-0DH, and RSX-QIDH models for hard-to-calculate excitation energies. We compare with the parent (B2-PLYP, B2GP-PLYP, PBE0-DH, and PBE-QIDH) and other (DSD-PBEP86) double-hybrid models as well as with some of the most widely employed hybrid functionals (B3LYP, PBE0, M06-2X, and ωB97X). For this purpose, we select a number of medium-sized intra- and inter-molecular charge-transfer excitations, which are known to be challenging to calculate using time-dependent density-functional theory (TD-DFT) and for which accurate reference values are available. We assess whether the high accuracy shown by the newest double-hybrid models is also confirmed for those cases too. We find that asymptotically corrected double-hybrid models yield a superior performance, especially for the inter-molecular charge-transfer excitation energies, as compared to standard double-hybrid models. Overall, the PBE-QIDH and its corresponding range-separated RSX-QIDH functional are recommended for general-purpose TD-DFT applications, depending on whether long-range effects are expected to play a significant role.

Theoretical insights on acceptor-donor dyads for organic photovoltaics

The field of organic photovoltaics has witnessed a steady growth in the last few decades and a recent renewal with the blossoming of single-material organic solar cells (SMOSCs). However, due to the intrinsic complexity of these devices (both in terms of their size and of the condensed phases involved), computational approaches to accurately predict their geometrical and electronic structure and to link their microscopic properties to the observed macroscopic behaviour are still lacking. In this work, we have focused on the rationalization of transport dynamics and we have set up a computational approach that makes a combined use of classical simulations and Density Functional Theory with the aim of disclosing the most relevant electronic and structural features of dyads used for SMOSC applications. As a prototype dyad, we have considered a molecule that consists in a dithiafulvalene-functionalized diketopyrrolopyrrole (DPP), acting as an electron donor, covalently linked to a fulleropyrrolidine (Ful), the electron acceptor. Our results, beside a quantitative agreement with experiments, show that the overall observed mobilities result from the competing packing mechanisms of the constituting units within the dyad both in the case of crystalline and amorphous phases. As a consequence, not all stable polymorphs have the same efficiency in transporting holes or electrons which often results in a highly directional carrier transport that is not, in general, a desirable feature for polycrystalline thin-films. The present work, linking microscopic packing to observed transport, thus opens the route for the in silico design of new dyads with enhanced and controlled structural and electronic features.

Subterahertz Momentum Drag and Violation of Matthiessen's Rule in an Ultraclean Ferromagnetic SrRuO_{3} Metallic Thin Film

SrRuO_{3}, a ferromagnet with an approximately 160 K Curie temperature, exhibits a T^{2}-dependent dc resistivity below ≈30  K. Nevertheless, previous optical studies in the infrared and terahertz range show non-Drude dynamics at low temperatures, which seem to contradict Fermi-liquid predictions. In this work, we measure the low-frequency THz range response of thin films with residual resistivity ratios, ρ_{300K}/ρ_{4K}≈74. At temperatures below 30 K, we find both a sharp zero frequency mode which has a width narrower than k_{B}T/ℏ as well as a broader zero frequency Lorentzian that has at least an order of magnitude larger scattering. Both features have temperature dependences consistent with a Fermi liquid with the wider feature explicitly showing a T^{2} scaling. Above 30 K, there is a crossover to a regime described by a single Drude peak that we believe arises from strong interband electron-electron scattering. Such two channel Drude transport sheds light on reports of the violation of Matthiessen's rule and extreme sensitivity to disorder in metallic ruthenates.

Modeling the Electron Transfer Chain in an Artificial Photosynthetic Machine

The development of efficient artificial leaves relies on the subtle combination of molecular assemblies able to absorb sunlight, converting light energy into electrochemical potential energy and finally transducing it into accessible chemical energy. The electronic design of these charge transfer molecular machines is crucial to build a complex supramolecular architecture for the light energy conversion. Here, we present an ab initio simulation of the whole decay pathways of a recently proposed artificial molecular reaction center. A complete structural and energetic characterization has been carried out with methods based on density functional theory, its time-dependent version, and a broken-symmetry approach. On the basis of our findings we provide a revision of the pathway only indirectly postulated from an experimental point of view, along with unprecedented and significant insights on the electronic and nuclear structure of intramolecular charge-separated states, which are fundamental for the application of this molecular assembly in photoelectrochemical cells. Importantly, we unravel the molecular driving forces of the various charge transfer steps, in particular those leading to the proton-coupled electron transfer final product, highlighting key elements for the future design strategies of such molecular assays.

Mechanochromic LLDPE Films Doped with NIR Reflective Paliogen Black

The perylene bisimide derivative Paliogen Black (P-black) is proposed as a new chromogenic probe that shows visible (vis) and near-infrared (NIR) responses after mechanical solicitations of host linear low-density polyethylene (LLDPE) films. P-black is reported to display strong absorption in the vis spectrum and unusual reflective and cooling features in the NIR region. Uniaxial deformation of the 2.5, 5, and 10 wt% P-black/LLDPE films yields a dichroic absorption under polarized light with color variations attributed by the computational analysis to the distinct anisotropic behavior of the transition dipole moments of P-black chromophores. When LLDPE films are deformed, P-black aggregates reduce their size from ≈30-40 µm to ≈5-10 µm that, in turn, causes reflectivity losses of about 30-40% at the maximum elongation. This gives rise to warming of 5-6 °C of the locally oriented film placed in contact with a black substrate under the illumination with an IR lamp for 5 s. These features combined with the high sensitivity of the vis-NIR response toward mechanical solicitations render P-black as a new solution to detect uniaxial deformations of plastic films through both optical and thermal outputs.

Range-separated hybrid and double-hybrid density functionals: A quest for the determination of the range-separation parameter

We recently derived a new and simple route to the determination of the range-separation parameter in range-separated exchange hybrid and double-hybrid density functionals by imposing an additional constraint to the exchange-correlation energy to recover the total energy of the hydrogen atom [Brémond et al., J. Chem. Phys. 15, 201102 (2019)]. Here, we thoroughly assess this choice by statistically comparing the derived values of the range-separation parameters to the ones obtained using the optimal tuning (OT) approach. We show that both approaches closely agree, thus, confirming the reliability of ours. We demonstrate that it provides very close performances in the computation of properties particularly prone to the one- and many-electron self-interaction errors (i.e., ionization potentials). Our approach arises as an alternative to the OT procedure, conserving the accuracy and efficiency of a standard Kohn-Sham approach to density-functional theory computation.

Understanding the properties of dithienylethenes functionalized for supramolecular self-assembly: a molecular modeling study

A dithienylethene (DTE) photochromic compound functionalized by ureidopyrimidinone (UPy) quadruple hydrogen bonding blocks was synthesized by Takeshita and coworkers [Takeshita et al., Chem. Commun., 2005, 761] in order to form a light-responsive supramolecular self-assembling system. In solution, the formation of supramolecular assemblies was only observed for one DTE isomer, namely the closed-form isomer. To rationalize this experimental finding, with the help of Molecular Dynamics (MD) and (time-dependent) DFT calculations, the behaviour of open-form and closed-form monomers, dimers, hexamers and π-stacked dimers in solution is investigated. Our simulations show that, for the open-form oligomers, the progression of the supramolecular assembly is hindered due to (i) the possible formation of a very stable cyclic dimer for the open-form parallel isomer, (ii) the relative flexibility of the open-form oligomers compared to their closed-form counterparts, and (iii) the possible existence of π-stacked dimers that constitute bottlenecks blocking the progression of the supramolecular self-assembly.

Toward the Mechanistic Understanding of the Additives' Role on Ammonium Nitrate Decomposition: Calcium Carbonate and Calcium Sulfate as Case Studies

The reaction mechanism involved in the decomposition of ammonium nitrate (AN) in the presence of CaCO3 and CaSO4, commonly used for stabilization and the reduction of explosivity properties of AN, was theoretically investigated using a computational approach based on density functional theory. The presented computational results suggest that both carbonate and sulfate anions can intercept an acid proton from nitric acid issued from the first step of decomposition of AN, thus inhibiting its runaway decomposition and the generation of reactive species (radicals). The reaction then leads to the production of stable products, as experimentally observed. Our modeling outcomes allow for tracing a relationship between the capability of proton acceptance of both carbonate and sulfate anions and the macroscopic behavior of these two additives as inhibitor or inert in the AN mixture.

sp-hybridized carbon allotrope molecular structures: An ongoing challenge for density-functional approximations

The recent synthesis of a C₁₈ monocyclic ring constitutes a major breakthrough as a new all-carbon disclosed form. However, modern density functional theory approaches do not lead to the correct experimental polyynic structure and favor the cumulenic one instead. We demonstrate here that this serious drawback can be solved by recently developed range-separated nonempirical schemes, independently of which kind of functional is being applied (i.e., semilocal, hybrid, or double-hybrid).

Aggregation Effects on Pigment Coatings: Pigment Red 179 as a Case Study

Here, we have studied, with a combined experimental and computational approach, the effect of the crystal environment and aggregation on the electronic properties of Pigment Red 179, which affect both its color and optical energy gap. Spectra acquired in the near-infrared and visible range of energies suggest that this molecule is indeed a "cool" dye, which can be employed as a red pigment that provides effective color coverage to different substrates without contributing to their heating during light irradiation. Spectra acquired on different polymer mixtures at different pigment concentrations (i.e., 2.5-10 wt %) suggest that absorption features depend on chromophoric arrangements promoted by the strong intermolecular π-π interactions. Calculations, performed at the time-dependent density functional theory level, allowed to both attribute the nature of the electronic transitions causing the observed spectra involved and understand the effect of the environment. Indeed, the visible spectra of the pigment is dominated by two localized transitions, with negligible charge transfer for both a dye monomer and dimer either in vacuum or acetonitrile solution. Instead, models including the crystal environment of the pigment show the presence of a high-wavelength S1 ← S0 charge transfer transition between two adjacent molecules, in quantitative agreement with the experimental absorption energy of the crystal pigment.

Theoretical insights into inorganic-organic intercalation products of the layered perovskite HLaNb2O7: perspectives for hybrid proton conductors

The modification of metal oxide surfaces with organic moieties has been widely studied as a method of preparing organic-inorganic hybrid materials for various applications. Among the inorganic oxides, ion-exchangeable layered perovskites are particularly interesting, because of their appealing electronic and reactive properties. In particular, their protonated interlayer surface can be easily functionalized with organic groups allowing the production of stable hybrid materials. As a further step in the design of new inorganic-organic hybrid proton conductors, a combined experimental and theoretical study of two intercalated compounds (propanol and imidazole) in HLaNb2O7 is presented here. A generally very good agreement with the available experimental data is found in reproducing both structural features and 13C-NMR chemical shifts, and marked differences between the two considered intercalated compounds are evidenced, with possible important outcomes for proton conduction. Notably, the free imidazole molecules are easily protonated by the acidic protons present in the interlayer spacing, thus inhibiting an efficient charge transport mechanism. In order to overcome this problem, a model system has been considered, where the imidazoles are bound to the end of a butyl chain, the whole being intercalated between two perovskite layers. The obtained theoretical data suggest that, in such a system, proton transfer between two adjacent imidazoles is a barrierless process. These results could then open new perspectives for such hybrid proton conductors.

Range-separated hybrid density functionals made simple

In this communication, we present a new and simple route to derive range-separated exchange (RSX) hybrid and double hybrid density functionals in a nonempirical fashion. In line with our previous developments [Brémond et al., J. Chem. Theory Comput. 14, 4052 (2018)], we show that by imposing an additional physical constraint to the exchange-correlation energy, i.e., by enforcing to reproduce the total energy of the hydrogen atom, we are able to generalize the nonempirical determination of the range-separation parameter to a family of RSX hybrid density functionals. The success of the resulting models is illustrated by an accurate modeling of several molecular systems and properties, like ionization potentials, particularly prone to the one- and many-electron self-interaction errors.

Strain-engineering Mott-insulating La2CuO4

The transition temperature T_c of unconventional superconductivity is often tunable. For a monolayer of FeSe, for example, the sweet spot is uniquely bound to titanium-oxide substrates. By contrast for La_2-xSr_xCuO₄ thin films, such substrates are sub-optimal and the highest T_c is instead obtained using LaSrAlO₄. An outstanding challenge is thus to understand the optimal conditions for superconductivity in thin films: which microscopic parameters drive the change in T_c and how can we tune them? Here we demonstrate, by a combination of x-ray absorption and resonant inelastic x-ray scattering spectroscopy, how the Coulomb and magnetic-exchange interaction of La₂CuO₄ thin films can be enhanced by compressive strain. Our experiments and theoretical calculations establish that the substrate producing the largest T_c under doping also generates the largest nearest neighbour hopping integral, Coulomb and magnetic-exchange interaction. We hence suggest optimising the parent Mott state as a strategy for enhancing the superconducting transition temperature in cuprates.

Aggregation-caused quenching versus crystallization induced emission in thiazolo[5,4-b]thieno[3,2-e]pyridine (TTP) derivatives: theoretical insights

We report a QM (TD-DFT) and QM/QM' (ONIOM) study of the modulation of emission in a series of thiazolo[5,4,b]thieno[3,2-e]pyridine (TTP) derivatives [Huang et al., J. Mater. Chem. C, 2017, 14, 3456]. By computing the excitation energy transfer couplings and the Huang-Rhys (HR) factors, we rationalize the aggregation-caused quenching (ACQ) observed for the parent molecule and the crystallization-induced emission (CIE) observed for the derivatives presenting intra-molecular H-bonding. We also show that the CIE strategy relying on the rigidification of the arch-bridge-like stator should be considered with caution since it can promote the energy dissipation through vibrational motions.

Photophysical properties of fluorescent imaging biological probes of nucleic acids: SAC-CI and TD-DFT Study

Recently, exciton-controlled hybridization-sensitive fluorescent oligonucleotide (ECHO) probe, which shows strong emission in the near-infrared region via hybridization to the target DNA and/or RNA strand, has been developed. In this work, photophysical properties of the chromophores of these probes and the fluorescent mechanism have been investigated by the SAC-CI and TD-DFT calculations. Three fluorescent cyanine chromophores whose excitation is challenging for TD-DFT methods, have been examined regarding the photo-absorption and emission spectra. The SAC-CI method well reproduces the experimental values with respect to transition energies, while the quantitative prediction by TD-DFT calculations is difficult for these chromophores. Some stable structures of H-aggregate system were computationally located and two of the configurations were examined for the photo-absorption. The present results support for the assumption based on experimental measurement in which strong fluorescence is due to the monomer unit in nearly planar structure and its suppression of probes is to the H-aggregates of two exciton units. Stokes shifts of these three chromophores were qualitatively reproduced by the theoretical calculations, while the energy splitting due to H-aggregate in the hybridized probe was slightly overestimated. © 2018 Wiley Periodicals, Inc.

Communication: Accurate description of interaction energies and three-body effects in weakly bound molecular complexes by PBE-QIDH models

We apply a recently developed parameter-free double-hybrid density functional belonging to the quadratic-integrand double-hybrid model to calculate association energies (ΔE) and three-body effects (Δ³E) arising from intermolecular interactions in weakly bound supramolecular complexes (i.e., the dataset 3B-69). The model behaves very accurately for trimer association energies and is found to outperform widely used density functional approximations while approaching the accuracy of more costly ab initio methods for three-body effects. The results are further improved when we add some specific corrections for the remaining dispersion interactions, D3(BJ) or VV10 for two-body effects and Axilrod-Teller-Muto for three-body effects, leading to marginal deviations (less than 1 kcal/mol for ΔE and around 0.03-0.04 kcal/mol for Δ³E) with respect to benchmark results.

Range-Separated Double-Hybrid Functional from Nonempirical Constraints

On the basis of our previous developments in the field of nonempirical double hybrids, we present here a new exchange-correlation functional based on a range-separated model for the exchange part and integrating a nonlocal perturbative correction to the electron correlation contribution. Named RSX-QIDH, the functional is free from any kind of empirical parametrization. Its range-separation parameter is set to recover the total energy of the hydrogen atom, thus eliminating the self-interaction error for this one-electron system. Subsequent tests on some relevant benchmark data sets confirm that the self-interaction error is particularly low for RSX-QIDH. This new functional provides also correct dissociation profiles for charged rare-gas dimers and very accurate ionization potentials directly from Kohn-Sham orbital energies. Above all, these good results are not obtained at the expense of other properties. Indeed, further tests on standard benchmarks show that RSX-QIDH is competitive with the more empirical ωB97X-2 double hybrid and outperforms the parent LC-PBE long-range corrected hybrid, thus underlining the important role of the nonlocal perturbative correlation.

Interaction of osmium(ii) redox probes with DNA: insights from theory

In the course of developing ultrasensitive and quantitative electrochemical point-of-care analytical tools for genetic detection of infectious diseases, osmium(ii) metallointercalators were revealed to be suitable and efficient redox probes to monitor the in vitro DNA amplification [Defever etal, Anal. Chem., 2011, 83, 1815-1821]. In this work, we thus propose a complete computational protocol in order to evaluate the affinity between Os(ii) complexes with double-stranded DNA. This protocol is based on molecular dynamics, with the parametrization of the GAFF force field for the Os(ii) complexes presenting an octahedral environment with polypyridine ligands, and QM/QM' calculations to evaluate the binding energy. For three Os(ii) probes and different binding sites, molecular dynamics simulations and interaction energies calculated at the QM/QM' level are successively discussed and compared to experimental data in order to identify the most stable binding sites. The computational protocol we propose should then be used to design more efficient Os(ii) metallointercalators.