A theoretical study of the electronic spectra of pyridine and phosphabenzene

Dispersion Interactions in Exciton-Localized States. Theory and Applications to π-π* and n-π* Excited States

We address the problem of intermolecular interaction energy calculations in molecular complexes with localized excitons. Our focus is on the correct representation of the dispersion energy. We derive an extended Casimir-Polder formula for direct computation of this contribution through second order in the intermolecular interaction operator V̂. An alternative formula, accurate to infinite order in V̂, is derived within the framework of the adiabatic connection (AC) theory. We also propose a new parametrization of the VV10 nonlocal correlation density functional, so that it corrects the CASSCF energy for the dispersion contribution and can be applied to excited-state complexes. A numerical investigation is carried out for benzene, pyridine, and peptide complexes with the local exciton corresponding to the lowest π–π* or n– π* states. The extended Casimir-Polder formula is implemented in the framework of multiconfigurational symmetry-adapted perturbation theory, SAPT(MC). A SAPT(MC) analysis shows that the creation of a localized exciton affects mostly the electrostatic component of the interaction energy of investigated complexes. Nevertheless, the changes in Pauli repulsion and dispersion energies cannot be neglected. We verify the performance of several perturbation- and AC-based methods. Best results are obtained with a range-separated variant of an approximate AC approach employing extended random phase approximation and CASSCF wave functions.

Assessing the Performances of CASPT2 and NEVPT2 for Vertical Excitation Energies

Methods able to simultaneously account for both static and dynamic electron correlations have often been employed, not only to model photochemical events but also to provide reference values for vertical transition energies, hence allowing benchmarking of lower-order models. In this category, both the complete-active-space second-order perturbation theory (CASPT2) and the N-electron valence state second-order perturbation theory (NEVPT2) are certainly popular, the latter presenting the advantage of not requiring the application of the empirical ionization-potential-electron-affinity (IPEA) and level shifts. However, the actual accuracy of these multiconfigurational approaches is not settled yet. In this context, to assess the performances of these approaches, the present work relies on highly accurate (±0.03 eV) aug-cc-pVTZ vertical transition energies for 284 excited states of diverse character (174 singlet, 110 triplet, 206 valence, 78 Rydberg, 78 n → π*, 119 π → π*, and 9 double excitations) determined in 35 small- to medium-sized organic molecules containing from three to six non-hydrogen atoms. The CASPT2 calculations are performed with and without IPEA shift and compared to the partially contracted (PC) and strongly contracted (SC) variants of NEVPT2. We find that both CASPT2 with IPEA shift and PC-NEVPT2 provide fairly reliable vertical transition energy estimates, with slight overestimations and mean absolute errors of 0.11 and 0.13 eV, respectively. These values are found to be rather uniform for the various subgroups of transitions. The present work completes our previous benchmarks focused on single-reference wave function methods ( J. Chem. Theory Comput. 2018, 14, 4360; J. Chem. Theory Comput. 2020, 16, 1711), hence allowing for a fair comparison between various families of electronic structure methods. In particular, we show that ADC(2), CCSD, and CASPT2 deliver similar accuracies for excited states with a dominant single-excitation character.

An Ab Initio Multiple Cloning Method for Non-Adiabatic Excited-State Molecular Dynamics in NWChem

The recently developed ab initio multiple cloning (AIMC) approach based on the multiconfigurational Ehrenfest (MCE) method provides a powerful and accurate way of describing the excited-state dynamics of molecular systems. The AIMC method is a controlled approximation to nonadiabatic dynamics with a particular strength in the proper description of decoherence effects because of the branching of vibrational wavepackets at a level crossing. Here, we report a new implementation of the AIMC algorithm in the open source NWChem computational chemistry program. The framework combines linear-response time-dependent density functional theory with Ehrenfest mean-field theory to determine the equations of motion for classical trajectories. The multidimensional wave function is decomposed into a superposition of Gaussian coherent states guided by Ehrenfest trajectories (i.e., MCE approach), which can clone with fully quantum mechanical amplitudes and phases. By using an efficient time-derivative based nonadiabatic coupling approach within the AIMC method, all observables are calculated on-the-fly in the nonadiabatic molecular dynamics process. As a representative example, we apply our implementation to study the ultrafast photoinduced electronic and vibrational energy transfer in a pyridine molecule. The effects of the cloning procedure on electronic and vibrational coherence, relaxation and unidirectional energy transfer are discussed. This new AIMC implementation provides a high-level nonadiabatic molecular dynamics framework for simulating photoexcited dynamics in complex molecular systems and experimentally relevant ultrafast spectroscopic probes, such as nonlinear coherent optical and X-ray signals.

Koopmans'-Type Theorem in Kohn-Sham Theory with Optimally Tuned Long-Range-Corrected (LC) Functionals

In the present study, we have investigated the applicability of long-range-corrected (LC) functionals to a Kohn-Sham (KS) Koopmans'-type theorem. Specifically, we have examined the performance of optimally tuned LCgau-core functionals (in combination with BOP and PW86-PW91 exchange-correlation functionals) by calculating the ionization potential (IP) within the context of Koopmans' prediction. In the LC scheme, the electron repulsion operator, 1/r₁₂, is divided into short-range and long-range components using a standard error function, with a range separation parameter μ determining the weight of the two ranges. For each system that we have examined (H₂O, CO, benzene, N₂, HF, H₂CO, C₂H₄, and five-membered ring compounds cyclo-C₄H₄X, with X = CH₂, NH, O, and S, and pyridine), the value of μ is optimized to minimize the deviation of the negative HOMO energy from the experimental IP. Our results demonstrate the utility of optimally tuned LC functionals in predicting the IP of outer valence levels. The accuracy is comparable to that of highly accurate ab initio theory. However, our Koopmans' method is less accurate for the inner valence and core levels. Overall, our results support the notion that orbitals in KS-DFT, when obtained with the LC functional, provide an accurate one-electron energy spectrum. This method represents a one-electron orbital theory that is attractive in its simple formulation and effective in its practical application.

Simultaneous observation of nuclear and electronic dynamics by ultrafast electron diffraction

Simultaneous observation of nuclear and electronic motion is crucial for a complete understanding of molecular dynamics in excited electronic states. It is challenging for a single experiment to independently follow both electronic and nuclear dynamics at the same time. Here we show that ultrafast electron diffraction can be used to simultaneously record both electronic and nuclear dynamics in isolated pyridine molecules, naturally disentangling the two components. Electronic state changes (S₁→S₀ internal conversion) were reflected by a strong transient signal in small-angle inelastic scattering, and nuclear structural changes (ring puckering) were monitored by large-angle elastic diffraction. Supported by ab initio nonadiabatic molecular dynamics and diffraction simulations, our experiment provides a clear view of the interplay between electronic and nuclear dynamics of the photoexcited pyridine molecule.

A Mountaineering Strategy to Excited States: Highly Accurate Energies and Benchmarks for Medium Sized Molecules

Following our previous work focusing on compounds containing up to 3 non-hydrogen atoms [J. Chem. Theory Comput. 2018, 14, 4360-4379], we present here highly accurate vertical transition energies obtained for 27 molecules encompassing 4, 5, and 6 non-hydrogen atoms: acetone, acrolein, benzene, butadiene, cyanoacetylene, cyanoformaldehyde, cyanogen, cyclopentadiene, cyclopropenone, cyclopropenethione, diacetylene, furan, glyoxal, imidazole, isobutene, methylenecyclopropene, propynal, pyrazine, pyridazine, pyridine, pyrimidine, pyrrole, tetrazine, thioacetone, thiophene, thiopropynal, and triazine. To obtain these energies, we use equation-of-motion/linear-response coupled cluster theory up to the highest technically possible excitation order for these systems (CC3, EOM-CCSDT, and EOM-CCSDTQ) and selected configuration interaction (SCI) calculations (with tens of millions of determinants in the reference space), as well as the multiconfigurational n-electron valence state perturbation theory (NEVPT2) method. All these approaches are applied in combination with diffuse-containing atomic basis sets. For all transitions, we report at least CC3/aug-cc-pVQZ vertical excitation energies as well as CC3/aug-cc-pVTZ oscillator strengths for each dipole-allowed transition. We show that CC3 almost systematically delivers transition energies in agreement with higher-level methods with a typical deviation of ±0.04 eV, except for transitions with a dominant double excitation character where the error is much larger. The present contribution gathers a large, diverse, and accurate set of more than 200 highly accurate transition energies for states of various natures (valence, Rydberg, singlet, triplet, n → π*, π → π*, ...). We use this series of theoretical best estimates to benchmark a series of popular methods for excited state calculations: CIS(D), ADC(2), CC2, STEOM-CCSD, EOM-CCSD, CCSDR(3), CCSDT-3, CC3, and NEVPT2. The results of these benchmarks are compared to the available literature data.

A benchmark study of electronic excitation energies, transition moments, and excited-state energy gradients on the nicotine molecule

In this work, we report a comparative study of computed excitation energies, oscillator strengths, and excited-state energy gradients of (S)-nicotine, chosen as a test case, using multireference methods, coupled cluster singles and doubles, and methods based on time-dependent density functional theory. This system was chosen because its apparent simplicity hides a complex electronic structure, as several different types of valence excitations are possible, including n-π(*), π-π(*), and charge-transfer states, and in order to simulate its spectrum it is necessary to describe all of them consistently well by the chosen method.

Multiconfigurational Second-Order Perturbation Theory Restricted Active Space (RASPT2) Method for Electronic Excited States: A Benchmark Study

The recently developed second-order perturbation theory restricted active space (RASPT2) method has been benchmarked versus the well-established complete active space (CASPT2) approach. Vertical excitation energies for valence and Rydberg excited states of different groups of organic (polyenes, acenes, heterocycles, azabenzenes, nucleobases, and free base porphin) and inorganic (nickel atom and copper tetrachloride dianion) molecules have been computed at the RASPT2 and multistate (MS) RASPT2 levels using different reference spaces and compared with CASPT2, CCSD, and experimental data in order to set the accuracy of the approach, which extends the applicability of multiconfigurational perturbation theory to much larger and complex systems than previously. Relevant aspects in multiconfigurational excited state quantum chemistry such as the valence-Rydberg mixing problem in organic molecules or the double d-shell effect for first-row transition metals have also been addressed.

Prediction of Excitation Energies for Conjugated Oligomers and Polymers from Time-Dependent Density Functional Theory

With technological advances, light-emitting conjugated oligomers and polymers have become competitive candidates in the commercial market of light-emitting diodes for display and other technologies, due to the ultralow cost, light weight, and flexibility. Prediction of excitation energies of these systems plays a crucial role in the understanding of their optical properties and device design. In this review article, we discuss the calculation of excitation energies with time-dependent density functional theory, which is one of the most successful methods in the investigation of the dynamical response of molecular systems to external perturbation, owing to its high computational efficiency.

Theoretical study of the vertical excited states of benzene, pyrimidine, and pyrazine by the symmetry adapted cluster—Configuration interaction method

The ground state and the excited states of benzene, pyrimidine, and pyrazine have been examined by using the symmetry adapted cluster-configuration interaction (SAC-CI) method. Detailed characterizations and the structures of the absorption peaks in the vacuum ultraviolet (VUV), low energy electron impact (LEEI), and electron energy loss (EEL) spectra were theoretically clarified by calculating the excitation energy and the oscillator strength for each excited state. We show that SAC-CI has the power to well reproduce the electronic excitation spectra (VUV, LEEI, and EEL) simultaneously to an accuracy for both the singlet and the triplet excited states originated from the low-lying pi --> pi*, n --> pi*, pi --> sigma* and n --> sigma* excited states of the titled compounds. The present results are compared with those of the previous theoretical studies by methods, such as EOM-CCSD(T), STEOM-CCSD, CASPT2 and TD-B3LYP, etc.

Ab initio (CASPT2) excited state calculations, including circular dichroism, of helically twisted cyanine dyes

Ab initio calculations at the CASSCF/CASPT2 level were performed on helically twisted mono-, tri-, and pentamethine cyanine dyes in the all-Z-configurations. Excitation energies and oscillator and rotatory strengths were calculated for the five lowest energy singlet states. Both the long wavelength methine band and the cis-band could be identified unambiguously from their configurational parentage. The calculated state energies are within 0.09 eV of the experimental value for the methine band and within 0.16 eV for the cis-band. The calculated rotatory strengths of the methine band shows sign inversion as the length of the chromophore increases: negative for the short monomethine, strongly positive for the pentamethine. The trimethine presents a borderline case: the measured rotatory strength is almost nil, the calculated one depends on the geometry. There is good agreement between rotatory strengths calculated in the velocity and in the length formalism.

Theoretical characterization of photoisomerization channels of dimethylpyridines on the singlet and triplet potential energy surfaces

Photoexcitations and photoisomerizations due to low-lying n pi* and pi pi* excited states of dimethylpyridines are investigated by density functional theory, CASSCF, CASPT2 and MRCI methodologies. Mechanistic details for the formation of Dewar dimethylpyridines and the interconversions of dimethylpyridines are rationalized through the characterization of minima and transition states on the singlet and triplet potential energy surfaces of relevant intermediates. Our present theoretical schemes suggest that Mobius dimethylpyridine intermediate 14 and azabenzvalene intermediate 10 can serve as possible precursors to Dewar dimethylpyridines and singlet phototransposition products, respectively. The calculations suggest that an S1(pi pi*)/S0 conical intersection in dimethylpyridines 2 is involved in the formation of 14. An azabenzvalene 10 might be formed through S2(pi pi*)/S1(n pi*) interaction followed by an S1/S0 decay in dimethylpyridine 6. Calculated barriers of isomerizations from 14 to Dewar dimethylpyridine 7 and from 10 to 4 are 8.4 and 28.5 kcal mol(-1) at the B3LYP/6-311 G** level, respectively. In the suggested triplet multistage transposition mechanism, an out-of-plane distorted geometry 19 due to vibrational relaxation of the T1(3B1) excited state of 3,5-dimethylpyridine 6 is a precursor of the interconversion of 6 to 2.4-dimethylpyridine 4. The formation of a triplet azaprefulvene 21 with a barrier of 20.7 kcal mol(-1) is a key step during the triplet migration process leading to another out-of-plane distorted structure 27. Subsequent rearomatization of 27 completes the interconversion of 6 with 4. Present calculations provide some insight into the photochemistry of dimethylpyridines at 254 nm.