Search

Time-resolved near-edge X-ray absorption fine structure (TR-NEXAFS) spectroscopy is a powerful technique for studying photochemical reaction dynamics with femtosecond time resolution. In order to avoid ambiguity in TR-NEXAFS spectra from nonadiabatic dynamics simulations, core- and valence-excited states must be evaluated on equal footing and those valence states must also define the potential energy surfaces used in the nonadiabatic dynamics simulation. In this work, we demonstrate that hole-hole Tamm-Dancoff-approximated density functional theory (hh-TDA) is capable of directly simulating TR-NEXAFS spectroscopies. We apply hh-TDA to the excited-state dynamics of acrolein. We identify two pre-edge features in the oxygen K-edge TR-NEXAFS spectrum associated with the S2 (ππ*) and S1 (nπ*) excited states. We show that these features can be used to follow the internal conversion dynamics between the lowest three electronic states of acrolein. Due to the low, O(N2) apparent computational complexity of hh-TDA and our GPU-accelerated implementation, this method is promising for the simulation of pre-edge features in TR-NEXAFS spectra of large molecules and molecules in the condensed phase.

We study the excited state absorption (ESA) properties of the four DNA bases (thymine, cytosine, adenine, and guanine) by different single reference quantum mechanical methods, i.e. equation of motion coupled cluster singles and doubles (EOM-CCSD), singles, doubles and perturbative triples (EOM-CC3), and time-dependent density functional theory (TD-DFT), with the long-range corrected CAM-B3LYP functional. Preliminary results at the Tamm-Dancoff (TDA) CAM-B3LYP level using the maximum overlap method (MOM) are reported for Thymine. In the gas phase, the three methods predict similar One Photon Absorption (OPA) spectra, which are also consistent with the experimental results and with the most accurate computational studies available in the literature. The ESA spectra are then computed for the pp states (one for pyrimidine, two for purines) associated with the lowest energy absorption band, and for the close-lying np state. The EOM-CC3, EOM-CCSD and CAM-B3LYP methods provide similar ESA spectral patterns, which are also in qualitative agreement with literature RASPT2 results. Once validated in the gas phase, TD-CAM-B3LYP has been used to compute the ESA in chloroform, including solvent effect by the polarizable continuum model (PCM). The predicted OPA and ESA spectra in chloroform are very similar to those in the gas phase, most of the bands shifting by less than 0.1 eV, with a small increase of the intensities and a moderate destabilization of the np state. Finally, ESA spectra have been computed from the minima of the lowest energy pp state, and are consistent with the available experimental transient absorption spectra of the nucleosides in solution, providing a final validation of our computational approach.

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange–correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear–electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an “open teamware” model and an increasingly modular design. This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

In this article, we use two extensively studied systems, a retinal model system and azobenzene, to explore the use of coupled cluster models for describing ground and singlet excited state potential energy surfaces of photoswitchable systems. While not being suitable for describing nuclear dynamics of photoisomerization, coupled cluster models have useful attributes, such as the inclusion of dynamical correlation, their black box nature, and the systematic improvement offered by truncation level. Results for the studied systems show that when triple excitations (here through the CC3 model) are included, ground and excited state potential energy surfaces for isomerization paths may reliably be generated, also for states of doubly excited character. For ground state equilibrium cis- and trans-azobenzene, the molecular geometry and basis set is seen to significantly impact the vertical excitation energies for the two lowest excited states. Efficient implementations of coupled cluster models can therefore constitute valuable tools for investigating photoswitchable systems and can be used for preliminary black box studies to gather information before more complicated excited state dynamics approaches are pursued.

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 third-order algebraic-diagrammatic construction method for studies of electron detachment processes within the electron propagator framework [IP-ADC(3)] was extended to treat the properties of molecular states with a detached electron using the intermediate state representation (ISR) formalism. The second-order ISR(2) equations for the one-particle (transition) density matrix have been derived and implemented as an extension of the IP-(U)ADC(3) method available in the Q-CHEM program. As a first systematic test of the present IP-(U)ADC(3)/ISR(2) method, the dipole moments of various electronic states of closed- and open-shell molecules have been computed and compared to full configuration interaction (FCI) results. The present study employing FCI benchmarks also provides the first rigorous estimates for the accuracy of electron detachment energies obtained using the IP-ADC(3) method.

The third-order non-Dyson algebraic-diagrammatic construction approach to the electron propagator [IP-ADC(3)] is extended using the intermediate state representation (ISR) formalism, allowing the wave functions and properties of molecular states with detached electron to be studied. The second-order ISR equations [ISR(2)] for the one-particle (transition) density matrix have been derived and implemented in the Q-CHEM program. The approach is completely general and enables evaluation of arbitrary one-particle operators and interpretation of electron detachment processes in terms of density-based quantities. The IP-ADC(3)/ISR(2) equations were implemented for Ŝz-adapted intermediate states, allowing open-shell molecules to be studied using unrestricted Hartree-Fock references. As a first test for computations of ground state properties, dipole moments of various closed- and open-shell molecules have been computed by means of electron detachment from the corresponding anions. The results are in good agreement with experimental data. The potential of IP-ADC(3)/ISR(2) for the interpretation of photoelectron spectra is demonstrated for the galvinoxyl free radical.

We present a new and efficient implementation of the closed shell coupled cluster singles and doubles with perturbative triples method (CC3) in the electronic structure program eT. Asymptotically, a ground state calculation has an iterative cost of 4nV4nO3 floating point operations (FLOP), where nV and nO are the number of virtual and occupied orbitals, respectively. The Jacobian and transpose Jacobian transformations, required to iteratively solve for excitation energies and transition moments, both require 8nV4nO3 FLOP. We have also implemented equation of motion (EOM) transition moments for CC3. The EOM transition densities require recalculation of triples amplitudes, as nV3nO3 tensors are not stored in memory. This results in a noniterative computational cost of 10nV4nO3 FLOP for the ground state density and 26nV4nO3 FLOP per state for the transition densities. The code is compared to the CC3 implementations in CFOUR, DALTON, and PSI4. We demonstrate the capabilities of our implementation by calculating valence and core excited states of l-proline.

We present a computational study of the one-photon and excited-state absorption (ESA) from the two lowest energy excited states of uracil in the gas phase: an nπ* dark state (1n) and the lowest energy bright ππ* state (1π). The predictions of six different linear response electronic structure methods, namely, TD-CAM-B3LYP, EOM-CCSD, EOM-CC3, ADC(2), ADC(2)-x, and ADC(3) are critically compared. In general, the spectral shapes predicted by TD-CAM-B3LYP, EOM-CCSD, EOM-CC3, and ADC(3) are fairly similar, though the quality of TD-CAM-B3LYP slightly deteriorates in the high-energy region. By computing the spectra at some key structures on different potential energy surfaces (PES), that is, the Franck-Condon point, the 1n minimum, and structures representative of different regions of the 1π PES, we obtain important insights into the shift of the ESA spectra, following the motion of the wavepacket on the excited-state PES. Though 1π has larger ESA than 1n, some spectral regions are dominated by these latter signals. Aside from its methodological interest, we thus obtain interesting indications to interpret transient absorption spectra to disentangle the photoactivated dynamics of nucleobases.