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.

A theoretical framework for understanding molecular structures is crucial for the development of new technologies such as catalysts or solar cells. Apart from electronic excitations energies however, only spectroscopic properties of molecules consisting of lighter elements can be computationally described at high level of theory today, since heavy elements require a relativistic framework and most methods have only been derived in a non-relativistic one so far. Important new technologies like the above mentioned require molecules that contain heavier elements and hence there is a great need for the development of relativistic computational methods at higher level of accuracy.

No abstract

We have measured, analyzed and simulated the ground state valence photoelectron (PES), X-ray absorption (XAS), X-ray photoelectron (XPS) as well as normal and resonant Auger-Meitner electron (AES) spectra of oxazole at the carbon, oxygen and nitrogen K-edge in order to understand its electronic structure. Experimental data is 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 that both N and O K-edge XAS are sensitive to the amount of dynamical electron correlation included in the theoretical description, and that for a complete description of XPS, additional orbital correlation and orbital relaxation effects need to be considered. The normal AES are dominated by a singlet excitation channel and well described by theory. The resonant AES, 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.

We present an efficient implementation of ground and excited state CCSD gradients based on Cholesky-decomposed electron repulsion integrals. Cholesky decomposition, like density-fitting, is an inner projection method, 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. The computational time spent per optimization cycle are furthermore found to be significantly lower compared to other implementations based on an inner projection method. We illustrate the capabilities of the implementation by optimizing the geometry of the retinal molecule (C20H28O) at the CCSD/aug-cc-pVDZ level of theory.