Calculation of quasi-diabatic states within the DFT/MRCI(2) framework: The QD-DFT/MRCI(2) method.

We describe a procedure for the calculation of quasi-diabatic states within the recently introduced DFT/MRCI(2) framework [S. P. Neville and M. S. Schuurman, J. Chem. Phys. 157, 164103 (2022)]. Based on an effective Hamiltonian formalism, the proposed procedure, which we term QD-DFT/MRCI(2), has the advantageous characteristics of being simultaneously highly efficient and effectively black box in nature while directly yielding both quasi-diabatic potentials and wave functions of high quality. The accuracy and efficiency of the QD-DFT/MRCI(2) formalism are demonstrated via the simulation of the vibronic absorption spectra of furan and chlorophyll a.

A DFT/MRCI Hamiltonian parameterized using only ab initio data: I. valence excited states.

A new combined density functional theory and multi-reference configuration interaction (DFT/MRCI) Hamiltonian parameterized solely using the benchmark ab initio vertical excitation energies obtained from the QUEST databases is presented. This new formulation differs from all previous versions of the method in that the choice of the underlying exchange-correlation (XC) functional employed to construct the one-particle (orbital) basis is considered, and a new XC functional, QTP17, is chosen for its ability to generate a balanced description of core and valence vertical excitation energies. The ability of the new DFT/MRCI Hamiltonian, termed QE8, to furnish accurate excitation energies is confirmed using benchmark quantum chemistry computations, and a mean absolute error of 0.16 eV is determined for the wide range of electronic excitations included in the validation dataset. In particular, the QE8 Hamiltonian dramatically improves the performance of DFT/MRCI for doubly excited states. The performance of fast approximate DFT/MRCI methods, p-DFT/MRCI and DFT/MRCI(2), is also evaluated using the QE8 Hamiltonian, and they are found to yield excitation energies in quantitative agreement with the parent DFT/MRCI method, with the two methods exhibiting a mean difference of 0.01 eV with respect to DFT/MRCI over the entire benchmark set.

Ultrafast Molecular Frame Quantum Tomography.

We develop and experimentally demonstrate a methodology for a full molecular frame quantum tomography (MFQT) of dynamical polyatomic systems. We exemplify this approach through the complete characterization of an electronically nonadiabatic wave packet in ammonia (NH_{3}). The method exploits both energy and time-domain spectroscopic data, and yields the lab frame density matrix (LFDM) for the system, the elements of which are populations and coherences. The LFDM fully characterizes electronic and nuclear dynamics in the molecular frame, yielding the time- and orientation-angle dependent expectation values of any relevant operator. For example, the time-dependent molecular frame electronic probability density may be constructed, yielding information on electronic dynamics in the molecular frame. In NH_{3}, we observe that electronic coherences are induced by nuclear dynamics which nonadiabatically drive electronic motions (charge migration) in the molecular frame. Here, the nuclear dynamics are rotational and it is nonadiabatic Coriolis coupling which drives the coherences. Interestingly, the nuclear-driven electronic coherence is preserved over longer timescales. In general, MFQT can help quantify entanglement between electronic and nuclear degrees of freedom, and provide new routes to the study of ultrafast molecular dynamics, charge migration, quantum information processing, and optimal control schemes.

Time-resolved photoelectron spectroscopy: the continuing evolution of a mature technique.

Time-resolved photoelectron spectroscopy (TRPES) has become one of the most widespread techniques for probing nonadiabatic dynamics in the excited electronic states of molecules. Furthermore, the complementary development of ab initio approaches for the simulation of TRPES signals has enabled the interpretation of these transient spectra in terms of underlying coupled electronic-nuclear dynamics. In this perspective, we discuss the current state-of-the-art approaches, including efforts to push femtosecond pulses into vacuum ultraviolet and soft X-ray regimes as well as the utilization of novel polarizations to use time-resolved optical activity as a probe of nonadiabatic dynamics. We close this perspective with a forward-looking prospectus on the new areas of application for this technique.

Resolving competing conical intersection pathways: time-resolved X-ray absorption spectroscopy of trans-1,3-butadiene.

Time-resolved X-ray absorption spectroscopy is emerging as a uniquely powerful tool to probe coupled electronic-nuclear dynamics in photo-excited molecules. Theoretical studies to date have established that time-resolved X-ray absorption spectroscopy is an atom-specific probe of excited-state wave packet passage through a seam of conical intersections (CIs). However, in many molecular systems, there are competing dynamical pathways involving CIs of different electronic and nuclear character. Discerning these pathways remains an important challenge. Here, we demonstrate that time-resolved X-ray absorption spectroscopy (TRXAS) has the potential to resolve competing channels in excited-state non-adiabatic dynamics. Using the example of 1,3-butadiene, we show how TRXAS discerns the different electronic structures associated with passage through multiple conical intersections. trans-1,3-Butadiene exhibits a branching between polarized and radicaloid pathways associated with ethylenic "twisted-pyramidalized" and excited-state cis-trans isomerization dynamics, respectively. The differing electronic structures along these pathways give rise to different XAS signals, indicating the possibility of resolving them. Furthermore, this indicates that XAS, and other core-level spectroscopic techniques, offer the appealing prospect of directly probing the effects of selective chemical substitution and its ability to affect chemical control over excited-state molecular dynamics.

Electronic relaxation and dissociation dynamics in formaldehyde: pump wavelength dependence.

The effect of the incident UV pump wavelength on the subsequent excited state dynamics, electronic relaxation, and ultimate dissociation of formaldehyde is studied using first principles simulation and Coulomb explosion imaging (CEI) experiments. Transitions in a vibronic progression in the à ← X̃ absorption band are systematically prepared using a tunable UV source which generates pulses centered at 304, 314, 329, and 337 nm. We find, both via ab initio simulation and experimental results, that the rate of excited state decay and subsequent dissociation displays a prominent dependence on which vibronic transition in the absorption band is prepared by the pump. Our simulations predict that nonadiabatic transition rates and dissociation yields will increase by a factor of >100 as the pump wavelength is decreased from 337 to 304 nm. The experimental results and theoretical simulations are in broad agreement and both indicate that the dissociation rate plateaus rapidly after ≈2 ps following an ultrafast sub-ps rise.

A perturbative approximation to DFT/MRCI: DFT/MRCI(2).

We introduce a perturbative approximation to the combined density functional theory and multireference configuration interaction (DFT/MRCI) approach. The method, termed DFT/MRCI(2), results from the application of quasi-degenerate perturbation theory (QDPT) and the Epstein-Nesbet partitioning to the DFT/MRCI Hamiltonian matrix. The application of QDPT obviates the need to diagonalize the large DFT/MRCI Hamiltonian; electronic energies are instead obtained as the eigenvalues of a small effective Hamiltonian, affording an orders of magnitude savings in the computational cost. Most importantly, the DFT/MRCI(2) approximation is found to be of excellent accuracy, furnishing excitation energies with a root mean squared deviation from the canonical DFT/MRCI values of less than 0.03 eV for an extensive test set of organic molecules.

A quantum molecular movie: polyad predissociation dynamics in the VUV excited 3pσ2Σu state of NO2.

The optical formation of coherent superposition states, a wavepacket, can allow the study of zeroth-order states, the evolution of which exhibit structural and electronic changes as a function of time: this leads to the notion of a molecular movie. Intramolecular vibrational energy redistribution, due to anharmonic coupling between modes, is the molecular movie considered here. There is no guarantee, however, that the formed superposition will behave semi-classically (e.g. Gaussian wavepacket dynamics) or even as an intuitively useful zeroth-order state. Here we present time-resolved photoelectron spectroscopy (TRPES) studies of an electronically excited triatomic molecule wherein the vibrational dynamics must be treated quantum mechanically and the simple picture of population flow between coupled normal modes fails. Specifically, we report on vibronic wavepacket dynamics in the zeroth-order 3pσ2Σu Rydberg state of NO2. This wavepacket exemplifies two general features of excited state dynamics in polyatomic molecules: anharmonic multimodal vibrational coupling (forming polyads); nonadiabatic coupling between nuclear and electronic coordinates, leading to predissociation. The latter suggests that the polyad vibrational states in the zeroth-order 3p Rydberg manifold are quasi-bound and best understood to be scattering resonances. We observed a rapid dephasing of an initially prepared 'bright' valence state into the relatively long-lived 3p Rydberg state whose multimodal vibrational dynamics and decay we monitor as a function of time. Our quantum simulations, based on an effective spectroscopic Hamiltonian, describe the essential features of the multimodal Fermi resonance-driven vibrational dynamics in the 3p state. We also present evidence of polyad-specificity in the state-dependent predissociation rates, leading to free atomic and molecular fragments. We emphasize that a quantum molecular movie is required to visualize wavepacket dynamics in the 3pσ2Σu Rydberg state of NO2.

Propagative block diagonalization diabatization of DFT/MRCI electronic states.

We present a framework for the calculation of diabatic states using the combined density functional theory and multireference configuration interaction (DFT/MRCI) method. Due to restrictions present in the current formulation of the DFT/MRCI method (a lack of analytical derivative couplings and the inability to use non-canonical Kohn-Sham orbitals), most common diabatization strategies are not applicable. We demonstrate, however, that diabatic wavefunctions and potentials can be reliably calculated at the DFT/MRCI level of theory using a propagative variant of the block diagonalization diabatization method (P-BDD). The proposed procedure is validated via the calculation of diabatic potentials for LiH and the simulation of the vibronic spectrum of pyrazine. In both cases, the combination of the DFT/MRCI and P-BDD methods is found to correctly recover the non-adiabatic coupling effects of the problem.

Substituent effects on nonadiabatic excited state dynamics: Inertial, steric, and electronic effects in methylated butadienes.

The photochemical dynamics of double-bond-containing hydrocarbons is exemplified by the smallest alkenes, ethylene and butadiene. Chemical substituents can alter both decay timescales and photoproducts through a combination of inertial effects due to substituent mass, steric effects due to substituent size, and electronic (or potential) effects due to perturbative changes to the electronic potential energy surface. Here, we demonstrate the interplay of different substituent effects on 1,3-butadiene and its methylated derivatives using a combination of ab initio simulation of nonadiabatic dynamics and time-resolved photoelectron spectroscopy. The purely inertial effects of methyl substitution are simulated through the use of mass 15 "heavy-hydrogen" atoms. As expected from both inertial and electronic influences, the excited-state dynamics is dominated by pyramidalization at the unsubstituted carbon sites. Although the electronic effects of methyl group substitution are weak, they alter both decay timescales and branching ratios by influencing the initial path taken by the excited wavepacket following photoexcitation.

Directing excited state dynamics via chemical substitution: A systematic study of π-donors and π-acceptors at a carbon-carbon double bond.

Functional group substituents are a ubiquitous tool in ground-state organic chemistry often employed to fine-tune chemical properties and obtain desired chemical reaction outcomes. Their effect on photoexcited electronic states, however, remains poorly understood. To help build an intuition for these effects, we have studied ethylene, substituted with electron acceptor (cyano) and/or electron donor (methoxy) substituents, both theoretically and experimentally: using ab initio quantum molecular dynamics and time-resolved photoelectron spectroscopy. Our results show the consistent trend that photo-induced ethylenic dynamics is primarily localized to the carbon with the greater electron density. For doubly substituted ethylenes, the trend is additive when both substituents are located on opposite carbons, whereas the methoxy group (in concert with steric effects) dominates when both substituents are located on a single carbon atom. These results point to the development of rules for structure-dynamics correlations; in this case, a novel mechanistic ultrafast photochemistry for conjugated carbon chains employing long-established chemical concepts.

Dynamics at Conical Intersections.

The nonadiabatic coupling of electronic and vibrational degrees of freedom is the defining feature of electronically excited states of polyatomic molecules. Once considered a theoretical curiosity, conical intersections (CIs) are now generally accepted as being the dominant source of coupled charge and vibrational energy flow in molecular excited states. Passage through CIs leads to the conversion of electronic to vibrational energy, which drives the ensuing photochemistry, isomerization being a canonical example. It has often been remarked that the CI may be thought of as a transition state in the excited state. As such, we expect that both the direction and the velocity of approach to the CI will matter. We explore this suggestion by looking for dynamical aspects of passage through CIs and for analogies with well-known concepts from ground-state reaction dynamics. Great progress has been made in the development of both experimental techniques and ab initio dynamics simulations, to a degree that direct comparisons may now be made. Here we compare time-resolved photoelectron spectroscopy results with on-the-fly ab initio multiple spawning calculations of the experimental observables, thereby validating each. We adopt a phenomenological approach and specifically concentrate on the excited-state dynamics of the C=C bond in unsaturated hydrocarbons. In particular, we make use of selective chemical substitution (such as replacing an H atom by a methyl group) so as to alter the inertia of certain vibrations relative to others, thus systematically varying (mass-weighted) directions and velocities of approach to a CI. Chemical substituents, however, may affect both the nuclear and electronic components of the total wave function. The former, which we call an inertial effect, influences the direction and velocity of approach. The latter, which we call a potential effect, modifies the electronic structure and therefore the energetic location and topography of the potential energy surfaces involved. Using a series of examples, we discuss both types of effects. We argue that there is a need for dynamical pictures and simple models of nonadiabatic dynamics at CIs and hope that the phenomenology presented here will help inspire such developments.

Site-Selective Isomerization of Cyano-Substituted Butadienes: Chemical Control of Nonadiabatic Dynamics.

The photochemistry of organic chromophores generally involves the co-evolution of the electronic and nuclear degrees of freedom. To obtain a specific and predetermined photochemical reaction outcome, chemical substitution can be used to selectively alter the underlying electronic potential energy surfaces to favor a particular reaction pathway. We show using ab initio simulation that the substitution of s- trans-1,3-butadiene with a cyano group can effectively "direct" a molecular wavepacket to particular regions of the seam of conical intersection and either favor or inhibit the photoinitiated cis-trans isomerization. The substituent is able to effect this control due to the formation of transient charge-separated electronic structures that arise during the nonadiabatic dynamical process. The atomic site at which this charge develops can be selectively stabilized (or destabilized) depending on the location of the cyano substituent and gives rise to a single dominant decay pathway. This work aims to demonstrate how the application of known electron density effects to ultrafast dynamics may be used to obtain desired photochemical reactions and properties.

Efficient calculation of X-ray absorption spectra using Chebyshev-Slepian filter diagonalisation.

The efficient, yet accurate, simulation of X-ray absorption spectra represents a significant challenge for ab initio electronic structure methods. Conventional approaches involve the explicit calculation of all core-excited states spanning the energy range of interest, even though only a small number of these states will contribute appreciably to the spectrum. We here report a different approach, based on a time-independent Chebyshev filter diagonalization scheme, which allows for the X-ray absorption spectrum to be computed without the explicit calculation of the core-excited eigenstates. Furthermore, in a subsequent postprocessing calculation, selected peaks may be analyzed via the calculation of natural transition orbitals, if desired. The scheme presented here is based on a refinement of the time-independent Chebyshev filter diagonalization approach. Previous formulations of this method have been characterized by a requirement for significant "user input" via the (sometimes unintuitive) tuning of various numerical parameters. To circumvent this, we introduce a new class of filters based on discrete prolate spheroidal sequences. We demonstrate that the resulting method, which we term Chebyshev-Slepian filter diagonalization, makes filter diagonalization essentially a black-box procedure. The Chebyshev-Slepian filter diagonalization method is implemented at the second-order algebraic diagrammatic construction level of theory and validated through the calculation of the X-ray absorption spectra of trifluoroacetonitrile and 1,4-benzoquinone.

Electron transfer in photoexcited pyrrole dimers.

Following on from previous experimental and theoretical work [Neville et al., Nat. Commun. 7, 11357 (2016)], we report the results of a combined electronic structure theory and quantum dynamics study of the excited state dynamics of the pyrrole dimer following excitation to its first two excited states. Employing an exciton-based analysis of the Ã(π3s/σ*) and B̃(π3s/3p/σ*) states, we identify an excited-state electron transfer pathway involving the coupling of the Ã(π3s/σ*) and B̃(π3s/3p/σ*) states and driven by N-H dissociation in the B̃(π3s/3p/σ*) state. This electron transfer mechanism is found to be mediated by vibronic coupling of the B̃ state, which has a mixed π3s/3p Rydberg character at the Franck-Condon point, to a high-lying charge transfer state of the πσ* character by the N-H stretch coordinate. Motivated by these results, quantum dynamics simulations of the excited-state dynamics of the pyrrole dimer are performed using the multiconfigurational time-dependent Hartree method and a newly developed model Hamiltonian. It is predicted that the newly identified electron transfer pathway will be open following excitation to both the Ã(π3s/σ*) and B̃(π3s/3p/σ*) states and may be the dominant relaxation pathway in the latter case.

The simulation of X-ray absorption spectra from ground and excited electronic states using core-valence separated DFT/MRCI.

We present an extension of the combined density functional theory (DFT) and multireference configuration interaction (MRCI) method (DFT/MRCI) [S. Grimme and M. Waletzke, J. Chem. Phys. 111, 5645 (1999)] for the calculation of core-excited states based on the core-valence separation (CVS) approximation. The resulting method, CVS-DFT/MRCI, is validated via the simulation of the K-edge X-ray absorption spectra of 40 organic chromophores, amino acids, and nucleobases, ranging in size from CO2 to tryptophan. Overall, the CVS-DFT/MRCI method is found to yield accurate X-ray absorption spectra (XAS), with consistent errors in peak positions of ∼2.5-3.5 eV. Additionally, we show that the CVS-DFT/MRCI method may be employed to simulate XAS from valence excited states and compare the simulated spectra to those computed using the established wave function-based approaches [ADC(2) and ADC(2)x]. In general, each of the methods yields excited state XAS spectra in qualitative and often quantitative agreement. In the instances where the methods differ, the CVS-DFT/MRCI simulations predict intensity for transitions for which the underlying electronic states are characterized by doubly excited configurations relative to the ground state configuration. Here, we aim to demonstrate that the CVS-DFT/MRCI approach occupies a specific niche among numerous other electronic structure methods in this area, offering the ability to treat initial states of arbitrary electronic character while maintaining a low computational cost and comparatively black box usage.

Ultrafast X-Ray Spectroscopy of Conical Intersections.

Ongoing developments in ultrafast x-ray sources offer powerful new means of probing the complex nonadiabatically coupled structural and electronic dynamics of photoexcited molecules. These non-Born-Oppenheimer effects are governed by general electronic degeneracies termed conical intersections, which play a key role, analogous to that of a transition state, in the electronic-nuclear dynamics of excited molecules. Using high-level ab initio quantum dynamics simulations, we studied time-resolved x-ray absorption (TRXAS) and photoelectron spectroscopy (TRXPS) of the prototypical unsaturated organic chromophore, ethylene, following excitation to its S_{2}(ππ^{*}) state. The TRXAS, in particular, is highly sensitive to all aspects of the ensuing dynamics. These x-ray spectroscopies provide a clear signature of the wave packet dynamics near conical intersections, related to charge localization effects driven by the nuclear dynamics. Given the ubiquity of charge localization in excited state dynamics, we believe that ultrafast x-ray spectroscopies offer a unique and powerful route to the direct observation of dynamics around conical intersections.

A general approach for the calculation and characterization of x-ray absorption spectra.

We present a general approach for the calculation and assignment of X-ray absorption spectra based on electronic wavepacket propagations performed using explicitly time-dependent electronic structure calculations. Such calculations have the appeal of yielding the entire absorption spectrum for the cost of a single set of electronic wavepacket propagations, obviating the need to explicitly calculate large numbers of core-excited states. The spectrum can either be calculated from the Fourier transform of the time-dependent dipole moment or from the Fourier transform of the wavepacket autocorrelation function. We propose that calculating the absorption spectrum using the latter approach will generally be the preferred option. This method has two important advantages. First, the autocorrelation functions can be obtained for twice the propagation time, resulting in a halving of the computational effort required to calculate the spectrum relative to the time-dependent dipole moment approach. Second, using the tools of filter diagonalisation, the autocorrelation functions may be used to determine the time-independent final core-excited states underlying the peaks of interest in the spectrum. The proposed scheme is validated by calculating and characterizing the X-ray absorption spectra of benzene and trifluoroacetonitrile at the time-dependent second-order algebraic diagrammatic construction level of theory.

Vacuum ultraviolet excited state dynamics of the smallest ring, cyclopropane. I. A reinterpretation of the electronic spectrum and the effect of intensity borrowing.

Cyclopropane, the smallest organic ring compound, exhibits complex spectroscopy and excited state dynamics. In Paper I, we reinterpret the vacuum ultraviolet (VUV) electronic absorption spectrum of cyclopropane via ab initio computation. The first two bands in the VUV spectrum are simulated using wavepacket propagations employing the multiconfigurational time-dependent Hartee method and a newly parameterized linear vibronic coupling Hamiltonian. The parameters of the model Hamiltonian are obtained directly from high level multireference configuration interaction calculations. An analysis of the results, with an emphasis on previously neglected vibronic coupling effects, reveals that these vibronic coupling terms must be included in order to account for strong intensity borrowing effects. This treatment dramatically changes the assignment of much of the VUV spectrum, with intensity borrowing by the optically dark A 2 ' (σ3p x /3p y ) and A 1 ' (σ3p x /3p y ) states from the E'(σ3p x /3p y ) state being found to give rise to almost all the spectral intensities below 8 eV. This is in stark contrast to previous studies, which attributed the first two bands to transitions to the E'(σ3p x /3p y ) state. This highlights the limitations of assigning spectral features based solely on calculated electronic excitation energies and oscillator strengths. Furthermore, we address the significant but infrequently discussed difficulties involved in determining the electronic character of a wavepacket produced in the pump step of ultrafast pump-probe experiments for systems exhibiting strong vibronic coupling.

Excited state non-adiabatic dynamics of the smallest polyene, trans 1,3-butadiene. II. Ab initio multiple spawning simulations.

The excited state non-adiabatic dynamics of the smallest polyene, trans 1,3-butadiene (BD), has long been the subject of controversy due to its strong coupling, ultrafast time scales and the difficulties that theory faces in describing the relevant electronic states in a balanced fashion. Here we apply Ab Initio Multiple Spawning (AIMS) using state-averaged complete active space multistate second order perturbation theory [SA-3-CAS(4/4)-MSPT2] which describes both static and dynamic electron correlation effects, providing a balanced description of both the initially prepared bright 11Bu (ππ*) state and non-adiabatically coupled dark 21Ag state of BD. Importantly, AIMS allows for on-the-fly calculations of experimental observables. We validate our approach by directly simulating the time resolved photoelectron-photoion coincidence spectroscopy results presented in Paper I [A. E. Boguslavskiy et al., J. Chem. Phys. 148, 164302 (2018)], demonstrating excellent agreement with experiment. Our simulations reveal that the initial excitation to the 11Bu state rapidly evolves via wavepacket dynamics that follow both bright- and dark-state pathways as well as mixtures of these. In order to test the sensitivity of the AIMS results to the relative ordering of states, we considered two hypothetical scenarios biased toward either the bright 1Bu or the dark 21Ag state. In contrast with AIMS/SA-3-CAS(4/4)-MSPT2 simulations, neither of these scenarios yields favorable agreement with experiment. Thus, we conclude that the excited state non-adiabatic dynamics in BD involves both of these ultrafast pathways.