Error Control and Automatic Detection of Reference Active Spaces in Many-Body Expanded Full Configuration Interaction.

We present a wide-reaching revamp of the generalized many-body expanded full configuration interaction (MBE-FCI) method. First, we outline how to automatize the selection of reference active spaces, whereby the inherent bias introduced through a manual identification is reduced, also within the context of traditional complete active space methods. Second, we allow for the use of compact orbital clusters as expansion objects, which works to circumvent the unfavorable scaling with the number of orbitals included in the space complementary to the reference orbitals. Finally, we present a new algorithm for ensuring that many-body expansions can be efficiently terminated while conservatively accounting for resulting errors. These developments are all tested on a variety of molecular systems and different orbital representations to illustrate the abilities of our algorithm to produce correlation energies within predetermined error bounds, significantly broadening the overall applicability of the MBE-FCI method.

MBE-CASSCF Approach for the Accurate Treatment of Large Active Spaces.

We present a novel implementation of the complete active space self-consistent field (CASSCF) method that makes use of the many-body expanded full configuration interaction (MBE-FCI) method to incrementally approximate electronic structures within large active spaces. On the basis of a hybrid first-order algorithm employing both Super-CI and quasi-Newton strategies for the optimization of molecular orbitals, we demonstrate both computational efficacy and high accuracy of the resulting MBE-CASSCF method. We assess the performance of our implementation on a set of established numerical tests before applying MBE-CASSCF in the investigation of the triplet-quintet spin gap of iron(II) porphyrin with active spaces as large as 50 electrons in 50 orbitals.

Symmetrization of Localized Molecular Orbitals.

We present a novel algorithm for (i) detecting approximate symmetries inherently present among spatially localized molecular orbitals and (ii) enforcing these in numerically exact manners by means of unitary optimization techniques. The vast potential of our algorithm to compress a full set of molecular orbitals into only a minimal set of symmetry-unique orbitals is demonstrated, starting from localized bases of either Pipek-Mezey or Foster-Boys orbitals. Comparisons of results based on either of these two localization procedures indicate that Foster-Boys molecular orbitals can be spanned by a smaller number of symmetry-unique orbitals on average, making these outstanding candidates for the exploitation of general, (non-)Abelian point-group symmetries in a range of local correlation methods. As an illustration of said compressibility, our algorithm is capable of identifying a mere 14 symmetry-unique orbitals for the buckminsterfullerene in the highly symmetric Ih molecular point group, corresponding to only 1.7% of its total 840 molecular orbitals in a standard double-ζ basis set. The present work thus marks an important advancement in the exploitation of point-group symmetry within local correlation methods, for which the appropriate adaption of symmetry uniqueness among orbitals has the potential to yield unprecedented speed-ups.

Photophysical properties of the triangular [Au(HN[double bond, length as m-dash]COH)]3 complex and its dimer.

Rate constants for radiative and non-radiative transitions of the [Au(HN[double bond, length as m-dash]COH)]3 complex and its dimer were calculated within the Herzberg-Teller approximation based on quantum mechanical principles. A high triplet quantum yield was estimated for the monomer. Internal conversion (IC) was found to be the major competing process to the intersystem crossing (ISC) from the lowest excited singlet state (S1) to the lowest triplet state (T1). ISC and IC from the spin-mixed T[combining tilde]1 state also dominate the triplet relaxation process resulting in a negligible phosphorescence quantum yield for the monomer. The IC and ISC rate constants of the dimer are considerably smaller due to much lower Franck-Condon factors. For the dimer a triplet quantum yield of 0.71 was estimated using the extended multi-configuration quasi-degenerate second-order perturbation theory (XMCQDPT2) method to calculate the transition energies. Fluorescence is the major competing process to the ISC relaxation of the S1 state of the dimer. The ISC and IC processes are insignificant for the relaxation of the T1 state, resulting in unity phosphorescence quantum yield. The high triplet and phosphorescence quantum yields of the [Au(HN[double bond, length as m-dash]COH)]3 dimer make it and its higher oligomers potential candidates as dopants for phosphorescent organic light emitting diodes and as down-converters in solid-state lighting systems.

Correction: Calculation of vibrationally resolved absorption and fluorescence spectra of the rylenes.

Correction for 'Calculation of vibrationally resolved absorption and fluorescence spectra of the rylenes' by Jonas Greiner et al., Phys. Chem. Chem. Phys., 2020, 22, 2379-2385.

Calculation of vibrationally resolved absorption and fluorescence spectra of the rylenes.

A generating function method was used to simulate the vibrationally resolved absorption and emission spectra of perylene, terrylene and quaterrylene. This method operates on the basis of adiabatic excitation energies and electronic ground and excited state vibrational frequencies. These parameters were calculated using density functional theory with the PBE0 functional for perylene and terrylene and with the BH-LYP functional for quaterrylene. The vertical excitation energies of the lower excited states were calculated using functionals with differing amounts of Hartree-Fock exchange. The optimal functional for each molecule was chosen by comparing these energies to literature excitation energies. Using this technique the calculated absorption spectra and the calculated emission spectrum of perylene were found to be in excellent agreement with the literature experimental spectra after introducing a shift and a scaling factor. The most prominent bands of the absorption spectra were assigned to their respective vibronic transitions.