Ultrafast internal conversion in 4-aminobenzonitrile occurs sequentially along the seam

A mountaineering strategy to excited states: Accurate vertical transition energies and benchmarks for substituted benzenes

In an effort to expand the existing QUEST database of accurate vertical transition energies [Véril et al. WIREs Comput. Mol. Sci. 2021, 11, e1517], we have modeled more than 100 electronic excited states of different natures (local, charge-transfer, Rydberg, singlet, and triplet) in a dozen of mono- and di-substituted benzenes, including aniline, benzonitrile, chlorobenzene, fluorobenzene, nitrobenzene, among others. To establish theoretical best estimates for these vertical excitation energies, we have employed advanced coupled-cluster methods including iterative triples (CC3 and CCSDT) and, when technically possible, iterative quadruples (CC4). These high-level computational approaches provide a robust foundation for benchmarking a series of popular wave function methods. The evaluated methods all include contributions from double excitations (ADC(2), CC2, CCSD, CIS(D), EOM-MP2, STEOM-CCSD), along with schemes that also incorporate perturbative or iterative triples (ADC(3), CCSDR(3), CCSD(T)(a)⋆ $$ {}^{\star } $$ , and CCSDT-3). This systematic exploration not only broadens the scope of the QUEST database but also facilitates a rigorous assessment of different theoretical approaches in the framework of a homologous chemical series, offering valuable insights into the accuracy and reliability of these methods in such cases. We found that both ADC(2.5) and CCSDT-3 can provide very consistent estimates, whereas among less expensive methods SCS-CC2 is likely the most effective approach. Importantly, we show that some lower order methods may offer reasonable trends in the homologous series while providing quite large average errors, and vice versa. Consequently, benchmarking the accuracy of a model based solely on absolute transition energies may not be meaningful for applications involving a series of similar compounds.

Comparison of Spin-Flip TDDFT-Based Conical Intersection Approaches with XMS-CASPT2

Determining conical intersection geometries is of key importance to understanding the photochemical reactivity of molecules. While many small- to medium-sized molecules can be treated accurately using multireference approaches, larger molecules require a less computationally demanding approach. In this work, minimum energy crossing point conical intersection geometries for a series of molecules have been studied using spin-flip TDDFT (SF-TDDFT), within the Tamm-Dancoff Approximation, both with and without explicit calculation of nonadiabatic coupling terms, and compared with both XMS-CASPT2 and CASSCF calculated geometries. The less computationally demanding algorithms, which do not require explicit calculation of the nonadiabatic coupling terms, generally fare well with the XMS-CASPT2 reference structures, while the relative energetics are only reasonably replicated with the MECP structure as calculated with the BHHLYP functional and full nonadiabatic coupling terms. We also demonstrate that, occasionally, CASSCF structures deviate quantitatively from the XMS-CASPT2 structures, showing the importance of including dynamical correlation.

Simulating the Nonadiabatic Relaxation Dynamics of 4-(N,N-Dimethylamino)benzonitrile (DMABN) in Polar Solution

The compound 4-(N,N-dimethylamino)benzonitrile (DMABN) represents the archetypal system for dual fluorescence, a rare photophysical phenomenon in which a given fluorophore shows two distinct emission bands. Despite extensive studies, the underlying mechanism remains the subject of debate. In the present contribution, we address this issue by simulating the excited-state relaxation process of DMABN as it occurs in polar solution. The potential energy surfaces for the system are constructed with the use of the additive quantum mechanics/molecular mechanics (QM/MM) method, and the coupled dynamics of the electronic wave function and the nuclei is propagated with the semiclassical fewest switches surface hopping method. The DMABN molecule, which comprises the QM subsystem, is treated with the use of the second-order algebraic diagrammatic construction (ADC(2)) method with the imposition of spin-opposite scaling (SOS). It is verified that this level of theory achieves a realistic description of the excited-state potential energy surfaces of DMABN. The simulation results qualitatively reproduce the main features of the experimentally observed fluorescence spectrum, thus allowing the unambiguous assignment of the two fluorescence bands: the normal band is due to the near-planar locally excited (LE) structure of DMABN, while the so-called "anomalous" second band arises from the twisted intramolecular charge transfer (TICT) structure. The transformation of the LE structure into the TICT structure takes place directly via intramolecular rotation, and is not mediated by another excited-state structure. In particular, the oft-discussed rehybridized intramolecular charge transfer (RICT) structure, which is characterized by a bent nitrile group, does not play a role in the relaxation process.

A molecular perspective on Tully models for nonadiabatic dynamics

We present a series of standardized molecular tests for nonadiabatic dynamics, reminiscent of the one-dimensional Tully models proposed in 1990.

Using density based indexes to characterize excited states evolution

With the aim of offering new computational tools helping in the description of photochemical reactions and phenomena occurring at the excited state, we present in this work the capability of a density based index (Π) in locating decay channels from higher to lower excited states. The Π index, previously applied to disclose non-radiative decay channels from the first excited state to the ground state, is very simple in its formulation and can be evaluated, practically with no extra computational cost, and coupled to any quantum method able to provide excited states densities. Indeed, this index relies only on the knowledge of energetics and electron densities of the different electronic states involved in the decay. In the present work, we show the proficiency of the Π index in the general case of decay between excited states by applying it to two model systems well characterized both theoretically and experimentally. In both cases, this descriptor was successful in spotting the regions where excited states are more likely to decay, thus suggesting its potential interest for further application in the design of new compounds. © 2018 Wiley Periodicals, Inc.