A theoretical study of the time-resolved x-ray absorption spectrum of the photoionized BT-1T cation.

The time-resolved x-ray absorption spectrum of the BT-1T cation (BT-1T+) is theoretically simulated in order to investigate the charge transfer reaction of the system. We employ both trajectory surface hopping and quantum dynamics to simulate the structural evolution over time and the changes in the state populations. To compute the static x-ray absorption spectra (XAS) of the ground and excited states, we apply both the time-dependent density functional theory and the coupled cluster singles and doubles method. The results obtained are in good agreement between the methods. It is, furthermore, found that the small structural changes that occur during the reaction have little effect on the static XAS. Hence, the tr-XAS can be computed based on the state populations determined from a nuclear dynamics simulation and one set of static XAS calculations, utilizing the ground state optimized geometry. This approach can save considerable computational resources, as the static spectra need not to be calculated for all geometries. As BT-1T is a relatively rigid molecule, the outlined approach should only be considered when investigating non-radiative decay processes in the vicinity of the Franck-Condon point.

New Implementation of an Equation-of-Motion Coupled-Cluster Damped-Response Framework with Illustrative Applications to Resonant Inelastic X-ray Scattering.

We present an implementation of a damped response framework for calculating resonant inelastic X-ray scattering (RIXS) at the equation-of-motion coupled-cluster singles and doubles (CCSD) and second-order approximate coupled-cluster singles and doubles (CC2) levels of theory in the open-source program eT. This framework lays the foundation for future extension to higher excitation methods (notably, the coupled-cluster singles and doubles with perturbative triples, CC3) and to multilevel approaches. Our implementation adopts a fully relaxed ground state and different variants of the core-valence separation projection technique to address convergence issues. Illustrative results are compared with those obtained within the frozen-core core-valence separated approach, available in Q-Chem, as well as with experiment. The performance of the CC2 method is evaluated in comparison with that of CCSD. It is found that, while the CC2 method is noticeably inferior to CCSD for X-ray absorption spectra, the quality of the CC2 RIXS spectra is often comparable to that of the CCSD level of theory, when the same valence excited states are probed. Finally, we present preliminary RIXS results for a solvated molecule in aqueous solution.

Core spectroscopy of oxazole.

We have measured, analyzed, and simulated the ground state valence photoelectron spectrum, x-ray absorption (XA) spectrum, x-ray photoelectron (XP) spectrum as well as normal and resonant Auger-Meitner electron (AE) spectrum of oxazole at the carbon, oxygen, and nitrogen K-edge in order to understand its electronic structure. Experimental data are compared to theoretical calculations performed at the coupled cluster, restricted active space perturbation theory to second-order and time-dependent density functional levels of theory. We demonstrate (1) that both N and O K-edge XA spectra are sensitive to the amount of dynamical electron correlation included in the theoretical description and (2) that for a complete description of XP spectra, additional orbital correlation and orbital relaxation effects need to be considered. The normal AE spectra are dominated by a singlet excitation channel and well described by theory. The resonant AE spectra, however, are more complicated. While the participator decay channels, dominating at higher kinetic energies, are well described by coupled cluster theory, spectator channels can only be described satisfactorily using a method that combines restricted active space perturbation theory to second order for the bound part and a one-center approximation for the continuum.

Efficient implementation of molecular CCSD gradients with Cholesky-decomposed electron repulsion integrals.

We present an efficient implementation of ground and excited state coupled cluster singles and doubles (CCSD) gradients based on Cholesky-decomposed electron repulsion integrals. Cholesky decomposition and density fitting are both inner projection methods, and, thus, similar implementation schemes can be applied for both methods. One well-known advantage of inner projection methods, which we exploit in our implementation, is that one can avoid storing large V3O and V4 arrays by instead considering three-index intermediates. Furthermore, our implementation does not require the formation and storage of Cholesky vector derivatives. The new implementation is shown to perform well, with less than 10% of the time spent calculating the gradients in geometry optimizations. Furthermore, the computational time per optimization cycle is significantly lower compared to other implementations based on an inner projection method. We showcase the capabilities of the implementation by optimizing the geometry of the retinal molecule (C20H28O) at the CCSD/aug-cc-pVDZ level of theory.

The Second-Order-Polarization-Propagator-Approximation (SOPPA) in a four-component spinor basis.

A theoretical framework for understanding molecular structures is crucial for the development of new technologies such as catalysts or solar cells. Apart from electronic excitation energies, however, only spectroscopic properties of molecules consisting of lighter elements can be computationally described at a high level of theory today since heavy elements require a relativistic framework, and thus far, most methods have only been derived in a non-relativistic framework. Important new technologies such as those mentioned above require molecules that contain heavier elements, and hence, there is a great need for the development of relativistic computational methods at a higher level of accuracy. Here, the Second-Order-Polarization-Propagator-Approximation (SOPPA), which has proven to be very successful in the non-relativistic case, is adapted to a relativistic framework. The equations for SOPPA are presented in their most general form, i.e., in a non-canonical spin-orbital basis, which can be reduced to the canonical case, and the expressions needed for a relativistic four-component SOPPA are obtained. The equations are one-index transformed, giving more compact expressions that correspond to those already available for the four-component RPA. The equations are ready for implementation in a four-component quantum chemistry program, which will allow both linear response properties and excitation energies to be calculated relativistically at the SOPPA level.

RPA(D) and HRPA(D): Two new models for calculations of NMR indirect nuclear spin-spin coupling constants.

In this article, the RPA(D) and HRPA(D) models for the calculation of linear response functions are presented. The performance of the new RPA(D) and HRPA(D) models is compared to the performance of the established RPA, HRPA, and SOPPA models in calculations of indirect nuclear spin-spin coupling constants using the CCSD model as a reference. The doubles correction offers a significant improvement on both the RPA and HRPA models; however, the improvement is more dramatic in the case of the RPA model. For all coupling types investigated in this study, the results obtained using the HRPA(D) model are comparable in accuracy to those given by the SOPPA model, while requiring between 30% and 90% of the calculation time needed for SOPPA. The RPA(D) model, while of slightly lower accuracy compared to the CCSD model than HRPA(D), offered calculation times of only approximately 25% of those required for SOPPA for all the investigated molecules. © 2018 Wiley Periodicals, Inc.