Electronic spectroscopy of 1,3‐cyclopentadiene, 1,3‐cyclohexadiene and 1,3‐cycloheptadiene by electron impact

Mapping the photochemistry of cyclopentadiene: from theory to ultrafast X-ray scattering

The photoinduced ring-conversion reaction when cyclopentadiene (CP) is excited at 5.10 eV is simulated using surface-hopping semiclassical trajectories with XMS(3)-CASPT2(4,4)/cc-pVDZ electronic structure theory. In addition, PBE0/def2-SV(P) is employed for ground state propagation of the trajectories. The dynamics is propagated for 10 ps, mapping both the nonadiabatic short-time dynamics (<300 fs) and the increasingly statistical dynamics on the electronic ground state. The short-time dynamics yields a mixture of hot CP and bicyclo[2.1.0]pentene (BP), with the two products reached via different regions of the same conical intersection seam. On the ground state, we observe slow conversion from BP to CP which is modelled by RRKM theory with a transition state determined using PBE0/def2-TZVP. The CP products are furthermore associated with ground state hydrogen shifts and some H-atom dissociation. Finally, the prospects for detailed experimental mapping using novel ultrafast X-ray scattering experiments are discussed and observables for such experiments are predicted. In particular, we assess the possibility of retrieving electronic states and their populations alongside the structural dynamics.

A hybrid approach to excited-state-specific variational Monte Carlo and doubly excited states

We extend our hybrid linear-method/accelerated-descent variational Monte Carlo optimization approach to excited states and investigate its efficacy in double excitations. In addition to showing a superior statistical efficiency when compared to the linear method, our tests on small molecules show good energetic agreement with benchmark methods. We also demonstrate the ability to treat double excitations in systems that are too large for a full treatment by using selected configuration interaction methods via an application to 4-aminobenzonitrile. Finally, we investigate the stability of state-specific variance optimization against collapse to other states' variance minima and find that symmetry, Ansatz quality, and sample size all have roles to play in achieving stability.

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.

Multiconfiguration Pair-Density Functional Theory

We present a new theoretical framework, called Multiconfiguration Pair-Density Functional Theory (MC-PDFT), which combines multiconfigurational wave functions with a generalization of density functional theory (DFT). A multiconfigurational self-consistent-field (MCSCF) wave function with correct spin and space symmetry is used to compute the total electronic density, its gradient, the on-top pair density, and the kinetic and Coulomb contributions to the total electronic energy. We then use a functional of the total density, its gradient, and the on-top pair density to calculate the remaining part of the energy, which we call the on-top-density-functional energy in contrast to the exchange-correlation energy of Kohn-Sham DFT. Because the on-top pair density is an element of the two-particle density matrix, this goes beyond the Hohenberg-Kohn theorem that refers only to the one-particle density. To illustrate the theory, we obtain first approximations to the required new type of density functionals by translating conventional density functionals of the spin densities using a simple prescription, and we perform post-SCF density functional calculations using the total density, density gradient, and on-top pair density from the MCSCF calculations. Double counting of dynamic correlation or exchange does not occur because the MCSCF energy is not used. The theory is illustrated by applications to the bond energies and potential energy curves of H2, N2, F2, CaO, Cr2, and NiCl and the electronic excitation energies of Be, C, N, N(+), O, O(+), Sc(+), Mn, Co, Mo, Ru, N2, HCHO, C4H6, c-C5H6, and pyrazine. The method presented has a computational cost and scaling similar to MCSCF, but a quantitative accuracy, even with the present first approximations to the new types of density functionals, that is comparable to much more expensive multireference perturbation theory methods.

LONG-RANGE-CORRECTED TIME-DEPENDENT DENSITY FUNCTIONAL STUDY ON ELECTRONIC SPECTRA OF FIVE-MEMBERED RING COMPOUNDS AND FREE-BASE PORPHYRIN

Long-range-corrected time-dependent density functional theory (LC-TDDFT) was applied to five-membered ring compounds (cyclopentadiene, pyrrole, and furan molecules) and free-base porphyrin. The vertical π-π* and Rydberg excitation energies and corresponding oscillator strengths were calculated by LC-TDDFT. The LC scheme obviously improved the Rydberg excitation energies and oscillator strengths of these systems, which have been fairly underestimated by TDDFT with conventional pure and hybrid B3LYP functionals. On the whole, LC-TDDFT results were very close to the results of the ab initio symmetry-adapted cluster configuration interaction (SAC-CI) method for most excitations. However, LC-TDDFT is poor in describing doubly excited states such as the 1 1A1 state of five-membered ring compounds.