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.

Few-femtosecond electronic and structural rearrangements of CH4+ driven by the Jahn-Teller effect

The Jahn-Teller effect (JTE) is central to the understanding of the physical and chemical properties of a broad variety of molecules and materials. Whereas the manifestations of the JTE in stationary properties of matter are relatively well studied, the study of JTE-induced dynamics is still in its infancy, largely owing to its ultrafast and non-adiabatic nature. For example, the time scales reported for the distortion of CH 4 + from the initial T d geometry to a nominal C 2 v relaxed structure range from 1.85 fs over 10 ± 2 fs to 20 ± 7 fs. Here, by combining element-specific attosecond transient-absorption spectroscopy and quantum-dynamics simulations, we show that the initial electronic relaxation occurs within 5 fs and that the subsequent nuclear dynamics are dominated by the Q2 scissoring and Q1 symmetric stretching modes, which dephase in 41 ± 10 fs and 13 ± 3 fs, respectively. Significant structural relaxation is found to take place only along the e-symmetry Q2 mode. These results demonstrate that CH 4 + created by ionization of CH 4 is best thought of as a highly fluxional species that possesses a long-time-averaged vibrational distribution centered around a D 2 d structure. The methods demonstrated in our work provide guidelines for the understanding of Jahn-Teller driven non-adiabatic dynamics in other more complex systems.

Machine Learning Seams of Conical Intersection: A Characteristic Polynomial Approach

The machine learning of potential energy surfaces (PESs) has undergone rapid progress in recent years. The vast majority of this work, however, has been focused on the learning of ground state PESs. To reliably extend machine learning protocols to excited state PESs, the occurrence of seams of conical intersections between adiabatic electronic states must be correctly accounted for. This introduces a serious problem, for at such points, the adiabatic potentials are not differentiable to any order, complicating the application of standard machine learning methods. We show that this issue may be overcome by instead learning the coordinate-dependent coefficients of the characteristic polynomial of a simple decomposition of the potential matrix. We demonstrate that, through this approach, quantitatively accurate machine learning models of seams of conical intersection may be constructed.

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.

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.

Removing the Deadwood from DFT/MRCI Wave Functions: The p-DFT/MRCI Method

The combined density functional theory and multireference configuration interaction (DFT/MRCI) method is a powerful tool for the calculation of excited electronic states of large molecules. There exists, however, a large amount of superfluous configurations in a typical DFT/MRCI wave function. We show that this deadwood may be effectively removed using a simple configuration pruning algorithm based on second-order Epstein-Nesbet perturbation theory. The resulting method, which we denote p-DFT/MRCI, is shown to result in orders of magnitude saving in computational timings, while retaining the accuracy of the original DFT/MRCI method.

Vacuum Ultraviolet Excited State Dynamics of the Smallest Ketone: Acetone

We combined tunable vacuum-ultraviolet time-resolved photoelectron spectroscopy (VUV-TRPES) with high-level quantum dynamics simulations to disentangle multistate Rydberg-valence dynamics in acetone. A femtosecond 8.09 eV pump pulse was tuned to the sharp origin of the A₁(n3d_yz) band. The ensuing dynamics were tracked with a femtosecond 6.18 eV probe pulse, permitting TRPES of multiple excited Rydberg and valence states. Quantum dynamics simulations reveal coherent multistate Rydberg-valence dynamics, precluding simple kinetic modeling of the TRPES spectrum. Unambiguous assignment of all involved Rydberg states was enabled via the simulation of their photoelectron spectra. The A₁(ππ*) state, although strongly participating, is likely undetectable with probe photon energies ≤8 eV and a key intermediate, the A₂(nπ*) state, is detected here for the first time. Our dynamics modeling rationalizes the temporal behavior of all photoelectron transients, allowing us to propose a mechanism for VUV-excited dynamics in acetone which confers a key role to the A₂(nπ*) state.

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 CO₂ 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.

Efficient Solution of the Electronic Eigenvalue Problem Using Wavepacket Propagation

We report how imaginary time wavepacket propagation may be used to efficiently calculate the lowest-lying eigenstates of the electronic Hamiltonian. This approach, known as the relaxation method in the quantum dynamics community, represents a fundamentally different approach to the solution of the electronic eigenvalue problem in comparison to traditional iterative subspace diagonalization schemes such as the Davidson and Lanczos methods. In order to render the relaxation method computationally competitive with existing iterative subspace methods, an extended short iterative Lanczos wavepacket propagation scheme is proposed and implemented. In the examples presented here, we show that by using an efficient wavepacket propagation algorithm the relaxation method is, at worst, as computationally expensive as the commonly used block Davidson-Liu algorithm, and in certain cases, significantly less so.

Beyond structure: ultrafast X-ray absorption spectroscopy as a probe of non-adiabatic wavepacket dynamics

The excited state non-adiabatic dynamics of polyatomic molecules, leading to the coupling of structural and electronic dynamics, is a fundamentally important yet challenging problem for both experiment and theory. Ongoing developments in ultrafast extreme vacuum ultraviolet (XUV) and soft X-ray sources present new probes of coupled electronic-structural dynamics because of their novel and desirable characteristics. As one example, inner-shell spectroscopy offers localized, atom-specific probes of evolving electronic structure and bonding (via chemical shifts). In this work, we present the first on-the-fly ultrafast X-ray time-resolved absorption spectrum simulations of excited state wavepacket dynamics: photo-excited ethylene. This was achieved by coupling the ab initio multiple spawning (AIMS) method, employing on-the-fly dynamics simulations, with high-level algebraic diagrammatic construction (ADC) X-ray absorption cross-section calculations. Using the excited state dynamics of ethylene as a test case, we assessed the ability of X-ray absorption spectroscopy to project out the electronic character of complex wavepacket dynamics, and evaluated the sensitivity of the calculated spectra to large amplitude nuclear motion. In particular, we demonstrate the pronounced sensitivity of the pre-edge region of the X-ray absorption spectrum to the electronic and structural evolution of the excited-state wavepacket. We conclude that ultrafast time-resolved X-ray absorption spectroscopy may become a powerful tool in the interrogation of excited state non-adiabatic molecular dynamics.

Conformer-resolved quantum dynamics study of the photodissociation of 3-pyrroline

A model Hamiltonian based on the vibronic coupling model is developed to describe the excited state dynamics of 3-pyrroline. With the use of the method of improved relaxation in conjunction with the MCTDH wavepacket propagation algorithm, vibrational eigenstates corresponding to both the axial and equatorial conformers of 3-pyrroline are calculated and subsequently used in a conformer-resolved study of the photodissociation of 3-pyrroline following excitation to its S1(3s/πσ*) and S2(3px) states. In analogy with ammonia, the excited state dynamics of both conformers of 3-pyrroline are found to be dominated by the (quasi-) planarization of the molecule in its electronically excited states and predominantly diabatic behavior of dissociation mediated by a conical intersection between the S1 and S0 states.