Search

The fate of thymine upon excitation by ultraviolet radiation has been the subject of intense debate over the past three decades. Today, it is widely believed that its ultrafast excited state decay stems from a radiationless transition from the bright ${\pi}{\pi}^*$ state to a dark $n{\pi}^*$ state. However, conflicting theoretical predictions have made the experimental data difficult to interpret. Here we simulate the ultrafast dynamics in thymine at the highest level of theory to date, performing wavepacket dynamics with a new coupled cluster method. Our simulation confirms an ultrafast ${\pi}{\pi}^*$ to $n{\pi}^*$ transition (${\tau} = 41 \pm 14$ fs). Furthermore, the predicted oxygen-edge X-ray absorption spectra agree quantitatively with the experimental results. Our simulation also predicts an as-yet uncharacterized photochemical pathway: a ${\pi}{\sigma}^*$ channel that leads to hydrogen dissociation at one of the two N-H bonds in thymine. Similar behavior has been identified in other heteroaromatic compounds, including adenine, and several authors have speculated that a similar pathway may exist in thymine. However, this was never confirmed theoretically or experimentally. This prediction calls for renewed efforts to experimentally identify or exclude the presence of this channel., Comment: 42 pages, 23 figures

Simulations of photochemical reaction dynamics have been a challenge to the theoretical chemistry community for some time. In an effort to determine the predictive character of current approaches, we predict the results of an upcoming ultrafast diffraction experiment on the photodynamics of cyclobutanone after excitation to the lowest lying Rydberg state (S$_2$). A picosecond of nonadiabatic dynamics is described with ab initio multiple spawning. We use both time dependent density functional theory and equation-of-motion coupled cluster for the underlying electronic structure theory. We find that the lifetime of the S$_2$ state is more than a picosecond (with both TDDFT and EOM-CCSD). The predicted UED spectrum exhibits numerous structural features, but weak time dependence over the course of the simulations.

Directly imaging structural dynamics involving hydrogen atoms by ultrafast diffraction methods is complicated by their low scattering cross-sections. Here we demonstrate that megaelectronvolt ultrafast electron diffraction is sufficiently sensitive to follow hydrogen dynamics in isolated molecules. In a study of the photodissociation of gas phase ammonia, we simultaneously observe signatures of the nuclear and corresponding electronic structure changes resulting from the dissociation dynamics in the time-dependent diffraction. Both assignments are confirmed by ab initio simulations of the photochemical dynamics and the resulting diffraction observable. While the temporal resolution of the experiment is insufficient to resolve the dissociation in time, our results represent an important step towards the observation of proton dynamics in real space and time.

Electrocyclic reactions are characterized by the concerted formation and cleavage of both {\sigma} and {\pi} bonds through a cyclic structure. This structure is known as a pericyclic transition state for thermal reactions and a pericyclic minimum in the excited state for photochemical reactions. However, the structure of the pericyclic geometry has yet to be observed experimentally. We use a combination of ultrafast electron diffraction and excited state wavepacket simulations to image structural dynamics through the pericyclic minimum of a photochemical electrocyclic ring-opening reaction in the molecule {\alpha}-terpinene. The structural motion into the pericyclic minimum is dominated by rehybridization of two carbon atoms, which is required for the transformation from two to three conjugated {\pi} bonds. The {\sigma} bond dissociation largely happens after internal conversion from the pericyclic minimum to the electronic ground state. These findings may be transferrable to electrocyclic reactions in general., Comment: Combined manuscript and supplementary information

We develop a quartic-scaling implementation of coupled-cluster singles and doubles based on low-rank tensor hypercontraction (THC) factorizations of both the electron repulsion integrals (ERIs) and the doubles amplitudes. This extends our rank-reduced coupled-cluster method to incorporate higher-order tensor factorizations. The THC factorization of the doubles amplitudes accounts for most of the gain in computational efficiency as it is sufficient, in conjunction with a Cholesky decomposition of the ERIs, to reduce the computational complexity of most contributions to the CCSD amplitude equations. Further THC factorization of the ERIs reduces the complexity of certain terms arising from nested commutators between the doubles excitation operator and the two-electron operator. We implement this new algorithm using graphical processing units (GPUs) and demonstrate that it enables CCSD calculations for molecules with 250 atoms and 2500 basis functions using a single computer node. Further, we show that the new method computes correlation energies with comparable accuracy to the underlying RR-CCSD method.

We present the first investigation of excited state dynamics by resonant Auger-Meitner spectroscopy (also known as resonant Auger spectroscopy) using the nucleobase thymine as an example. Thymine is photoexcited in the UV and probed with X-ray photon energies at and below the oxygen K-edge. After initial photoexcitation to a {\pi}{\pi}* excited state, thymine is known to undergo internal conversion to an n{\pi}* excited state with a strong resonance at the oxygen K-edge, red-shifted from the ground state {\pi}* resonances of thymine (see our previous study Wolf et al., Nat. Commun., 2017, 8, 29). We resolve and compare the Auger-Meitner electron spectra associated both with the excited state and ground state resonances, and distinguish participator and spectator decay contributions. Furthermore, we observe simultaneously with the decay of the n{\pi}* state signatures the appearance of additional resonant Auger-Meitner contributions at photon energies between the n{\pi}* state and the ground state resonances. We assign these contributions to population transfer from the n{\pi}* state to a {\pi}{\pi}* triplet state via intersystem crossing on the picosecond timescale based on simulations of the X-ray absorption spectra in the vibrationally hot triplet state. Moreover, we identify signatures from the initially excited {\pi}{\pi}* singlet state which we have not observed in our previous study.

The de Broglie wave nature of matter is a paradigmatic example of fundamental quantum physics and enables precise measurements of forces, fundamental constants and even material properties. However, even though matter-wave interferometry is nowadays routinely realized in many laboratories, this feat has remained an outstanding challenge for the vast class of native polypeptides, the building blocks of life, which are ubiquitous in biology but fragile and difficult to handle. Here, we demonstrate the quantum wave nature of gramicidin, a natural antibiotic composed of 15 amino acids. Femtosecond laser desorption of a thin biomolecular film with intensities up to 1~TW/cm$^2$ transfers these molecules into a cold noble gas jet. Even though the peptide's de Broglie wavelength is as tiny as 350~fm, the molecular coherence is delocalized over more than 20 times the molecular size in our all-optical time-domain Talbot-Lau interferometer. We compare the observed interference fringes for two different interference orders with a model that includes both a rigorous treatment of the peptide's quantum wave nature as well as a quantum chemical assessment of its optical properties to distinguish our result from classical predictions. The successful realization of quantum optics with this polypeptide as a prototypical biomolecule paves the way for quantum-assisted molecule metrology and in particular the optical spectroscopy of a large class of biologically relevant molecules.

Singlet fission is a spin-allowed exciton multiplication process in organic semiconductors that converts one spin-singlet exciton to two triplet excitons. It offers the potential to enhance solar energy conversion by circumventing the Shockley-Queisser limit on efficiency. Recently, the mechanism of the primary singlet fission process in pentacene and its derivatives have been extensively investigated, however, the nature of the primary ultrafast process in singlet fission is still a matter of debate. Here, we study the singlet fission process in a pentacene film by employing a combination of transient-grating (TG) and two-dimensional (2D) electronic spectroscopy complemented by quantum chemical and nonadiabatic dynamics calculations. The high sensitivity of heterodyne detected TG spectroscopy enabled us to capture the vibrational coherence and to show that it mediates the transition from the singlet excited electronic state to the triplet-pair state. This coherent process is further examined by 2D electronic spectroscopy. Detailed analysis of the experimental data reveals that significant vibronic couplings of a few key modes in the low- and high-frequency region connect the excited singlet and triplet-pair states. Based on quantum chemical calculations, we identify these key intermolecular rocking modes along the longitudinal molecular axis between the pentacene molecules. They play the essential role of an electronic bridge between the singlet and triplet-pair states. Along with high-frequency local vibrations acting as tuning modes, these rocking motions drive the ultrafast dynamics at the multidimensional conical intersection in the singlet fission process.

The maturation of analytical derivative theory over the past few decades has enabled classical electronic structure theory to provide accurate and efficient predictions of a wide variety of observable properties. However, classical implementations of analytical derivative theory take advantage of explicit computational access to the approximate electronic wavefunctions in question, which is not possible for the emerging case of hybrid quantum/classical methods. Here, we develop an efficient Lagrangian-based approach for analytical first derivatives of hybrid quantum/classical methods using only observable quantities from the quantum portion of the algorithm. Specifically, we construct the key first-derivative property of the nuclear energy gradient for the recently-developed multistate, contracted variant of the variational quantum eigensolver (MC-VQE) within the context of the ab initio exciton model (AIEM). We show that a clean separation between the quantum and classical parts of the problem is enabled by the definition of an appropriate set of relaxed density matrices, and show how the wavefunction response equations in the quantum part of the algorithm (coupled-perturbed MC-VQE or CP-MC-VQE equations) are decoupled from the wavefunction response equations and and gradient perturbations in the classical part of the algorithm. We explore the magnitudes of the Hellmann-Feynman and response contributions to the gradients in quantum circuit simulations of MC-VQE+AIEM and demonstrate a quantum circuit simulator implementation of adiabatic excited state dynamics with MC-VQE+AIEM., Comment: 26 pages, 3 figures, 4 tables

We have investigated the structural dynamics in photoexcited 1,2-diiodotetrafluoroethane molecules (C2F4I2) in the gas phase experimentally using ultrafast electron diffraction and theoretically using FOMO-CASCI excited state dynamics simulations. The molecules are excited by an ultra-violet femtosecond laser pulse to a state characterized by a transition from the iodine 5p orbital to a mixed 5p|| hole and CF2 antibonding orbital, which results in the cleavage of one of the carbon-iodine bonds. We have observed, with sub-Angstrom resolution, the motion of the nuclear wavepacket of the dissociating iodine atom followed by coherent vibrations in the electronic ground state of the C2F4I radical. The radical reaches a stable classical (non-bridged) structure in less than 200 fs., Comment: 13 pages, 11 figures

We develop an extension of the variational quantum eigensolver (VQE) algorithm - multistate, contracted VQE (MC-VQE) - that allows for the efficient computation of the transition energies between the ground state and several low-lying excited states of a molecule, as well as the oscillator strengths associated with these transitions. We numerically simulate MC-VQE by computing the absorption spectrum of an ab initio exciton model of an 18-chromophore light-harvesting complex from purple photosynthetic bacteria., Comment: Updated with addl references and SI case study of H-aggregate BChl-a stack

We studied the reaction dynamics of a proposed prebiotic reaction theoretically. The chemical process involves acetone cyanohydrin or formalcyanohydrin reacting with hydrosulfide in an aqueous environment. Rate constants and populations of reactant and product bimolecular geometric orientations for the reactions were obtained by using density functional theory for the energies, transition state theory for the rates and matrix exponentiation as well as the hybrid tau-leaping algorithm for the population dynamics. The role of including the solvent explicitly versus implicitly was also investigated. We found that adding water or hydrogen sulfide molecules explicitly lowers the activation energy barrier and leads to a more efficient reaction pathway. In particular, hydrogen sulfide was a better proton donor. Further, when adding the solvent explicitly a synchronous reaction mechanism was observed while a concerted reaction mechanism was seen when using an implicit solvent model. Finally, we studied the role of including more than one reactant and product bimolecular orientation geometry in the dynamics. Including all pathways, reactant to reactant, product to product, reactant to product and product to reactant reduced the reaction half-time by two orders of magnitude and was therefore determined to be fundamental for an accurate description of the dynamics. Overall, we found that most reactions which involve formalcyanohydrin occur more rapidly or at the same speed as reactions which involve acetone cyanohydrin at room temperature.

We derive the crossing conditions at conical intersections between electronic states in coupled cluster theory, and show that if the coupled cluster Jacobian matrix is nondefective, two (three) independent conditions are correctly placed on the nuclear degrees of freedom for an inherently real (complex) Hamiltonian. Calculations using coupled cluster theory on an $2 {^{1}}A' / 3 {^{1}}A'$ conical intersection in hypofluorous acid illustrate the nonphysical artifacts associated with defects at accidental same-symmetry intersections. In particular, the observed intersection seam is folded about a space of the correct dimensionality, indicating that minor modifications to the theory are required for it to provide a correct description of conical intersections in general. We find that an accidental symmetry allowed $1 {^{1}}A" / 2 {^{1}}A"$ intersection in hydrogen sulfide is properly described, showing no artifacts as well as linearity of the energy gap to first order in the branching plane., Comment: 9 pages and 4 figures

The conductor-like polarization model (C-PCM) with switching/Gaussian smooth discretization is a widely used implicit solvation model in chemical simulations. However, its application in quantum mechanical calculations of large-scale biomolecular systems can be limited by computational expense of both the gas phase electronic structure and the solvation interaction. We have previously used graphical processing units (GPUs) to accelerate the first of these steps. Here, we extend the use of GPUs to accelerate electronic structure calculations including C-PCM solvation. Implementation on the GPU leads to significant acceleration of the generation of the required integrals for C-PCM. We further propose two strategies to improve the solution of the required linear equations: a dynamic convergence threshold and a randomized block-Jacobi preconditioner. These strategies are not specific to GPUs and are expected to be beneficial for both CPU and GPU implementations. We benchmark the performance of the new implementation using over 20 small proteins in solvent environment. Using a single GPU, our method evaluates the C-PCM related integrals and their derivatives more than 10X faster than a conventional CPU based implementation. Our improvements to the linear solver provide a further 3X acceleration. The overall calculations including C-PCM solvation require typically 20-40% more effort than their gas phase counterparts for moderate basis set and molecule surface discretization level. The relative cost of the C-PCM solvation correction decreases as the basis sets and/or cavity radii increase. Therefore description of solvation with this model should be routine. We also discuss applications to the study of the conformational landscape of an amyloid fibril., Comment: 36 pages, 12 figures

Hybrid quantum mechanical-molecular mechanical (QM/MM) simulations are widely used in studies of enzymatic catalysis. Until recently, it has been cost prohibitive to determine the asymptotic limit of key energetic and structural properties with respect to increasingly large QM regions. Leveraging recent advances in electronic structure efficiency and accuracy, we investigate catalytic properties in catechol O-methyltransferase, a representative example of a methyltransferase critical to human health. Using QM regions ranging in size from reactants-only (64 atoms) to nearly one-third of the entire protein (940 atoms), we show that properties such as the activation energy approach within chemical accuracy of the large-QM asymptotic limits rather slowly, requiring approximately 500-600 atoms if the QM residues are chosen simply by distance from the substrate. This slow approach to asymptotic limit is due to charge transfer from protein residues to the reacting substrates. Our large QM/MM calculations enable identification of charge separation for fragments in the transition state as a key component of enzymatic methyl transfer rate enhancement. We introduce charge shift analysis that reveals the minimum number of protein residues (ca. 11-16 residues or 200-300 atoms for COMT) needed for quantitative agreement with large-QM simulations. The identified residues are not those that would be typically selected using criteria such as chemical intuition or proximity. These results provide a recipe for a more careful determination of QM region sizes in future QM/MM studies of enzymes., Comment: 27 pages, 10 figures

Multiple time-scale algorithms exploit the natural separation of time-scales in chemical systems to greatly accelerate the efficiency of molecular dynamics simulations. Although the utility of these methods in systems where the interactions are described by empirical potentials is now well established, their application to ab initio molecular dynamics calculations has been limited by difficulties associated with splitting the ab initio potential into fast and slowly varying components. Here we show that such a timescale separation is possible using two different schemes: one based on fragment decomposition and the other on range separation of the Coulomb operator in the electronic Hamiltonian. We demonstrate for both water clusters and a solvated hydroxide ion that multiple time-scale molecular dynamics allows for outer time steps of 2.5 fs, which are as large as those obtained when such schemes are applied to empirical potentials, while still allowing for bonds to be broken and reformed throughout the dynamics. This permits computational speedups of up to 4.4x, compared to standard Born-Oppenheimer ab initio molecular dynamics with a 0.5 fs time step, while maintaining the same energy conservation and accuracy., Comment: 28 pages, 7 Figures

We recently showed that the Dephasing Representation (DR) provides an efficient tool for computing ultrafast electronic spectra and that further acceleration is possible with cellularization [M. \v{S}ulc and J. Van\'i\v{c}ek, Mol. Phys. 110, 945 (2012)]. Here we focus on increasing the accuracy of this approximation by first implementing an exact Gaussian basis method, which benefits from the accuracy of quantum dynamics and efficiency of classical dynamics. Starting from this exact method, the DR is derived together with ten other methods for computing time-resolved spectra with intermediate accuracy and efficiency. These methods include the Gaussian DR, an exact generalization of the DR, in which trajectories are replaced by communicating frozen Gaussian basis functions evolving classically with an average Hamiltonian. The newly obtained methods are tested numerically on time correlation functions and time-resolved stimulated emission spectra in the harmonic potential, pyrazine S0/S1 model, and quartic oscillator. Numerical results confirm that both the Gaussian basis method and the Gaussian DR increase the accuracy of the DR. Surprisingly, in chaotic systems the Gaussian DR can outperform the presumably more accurate Gaussian basis method, in which the two bases are evolved separately., Comment: 15 pages, 7 figures

Configuration-space matrix elements of N-body potentials arise naturally and ubiquitously in the Ritz-Galerkin solution of many-body quantum problems. For the common specialization of local, finite-range potentials, we develop the eXact Tensor HyperContraction (X-THC) method, which provides a quantized renormalization of the coordinate-space form of the N-body potential, allowing for a highly separable tensor factorization of the configuration-space matrix elements. This representation allows for substantial computational savings in chemical, atomic, and nuclear physics simulations, particularly with respect to difficult "exchange-like" contractions., Comment: Third version of the manuscript after referee's comments. In press in PRL. Main text: 4 pages, 2 figures, 1 table; Supplemental material (also included): 14 pages, 2 figures, 2 tables