Benchmark theoretical study of the ionization energies, electron affinities and singlet–triplet energy gaps of azulene, phenanthrene, pyrene, chrysene and perylene

Stabilization of phenanthrene anions in helium nanodroplets

It has been debated for years if the polycyclic aromatic hydrocarbon phenanthrene exists in its anionic form, or, in other words, if its electron affinity (EA) is positive or negative. In this contribution we confirm that the bare phenanthrene anion Ph⁻ created in a binary collision with an electron at room temperature has a lifetime shorter than microseconds. However, the embedding of neutral phenanthrene molecules in negatively charged helium nanodroplets enables the formation of phenanthrene anions by charge transfer processes and the stabilization of the latter in the ultracold environment. Gentle shrinking of the helium matrix of phenanthrene-doped HNDs by collisions with helium gas makes the bare Ph⁻ visible by high-resolution mass spectrometry. From these and previous measurements we conclude, that the EA of phenanthrene is positive and smaller than 24.55 meV.

Phenanthrene anions are stabilized in the ultracold environment of helium nanodroplets. Gentle shrinking of the helium matrix by collisions with helium gas makes the bare phenanthrene anion visible by high-resolution mass spectrometry.

Phenanthrene: establishing lower and upper bounds to the binding energy of a very weakly bound anion

Quite a few molecules do not form stable anions that survive the time needed for their detection; their electron affinities (EA) are either very small or negative. How does one measure the EA if the anion cannot be observed? Or, at least, can one establish lower and upper bounds to their EA? We propose two approaches that provide lower and upper bounds. We choose the phenanthrene (Ph) molecule whose EA is controversial. Through competition between helium evaporation and electron detachment in HenPh- clusters, formed in helium nanodroplets, we estimate the lower bound of the vertical detachment energy (VDE) of Ph- as about -3 meV. In the second step, Ph is complexed with calcium whose electron affinity is just 24.55 meV. When CaPh- ions are collided with a thermal gas of argon, one observes Ca- product ions but no Ph-, suggesting that the EA of Ph is below that of Ca.

A Mountaineering Strategy to Excited States: Highly Accurate Energies and Benchmarks for Bicyclic Systems

Pursuing our efforts to define highly accurate estimates of the relative energies of excited states in organic molecules, we investigate, with coupled-cluster methods including iterative triples (CC3 and CCSDT), the vertical excitation energies of 10 bicyclic molecules (azulene, benzoxadiazole, benzothiadiazole, diketopyrrolopyrrole, furofuran, phthalazine, pyrrolopyrrole, quinoxaline, tetrathiafulvalene, and thienothiophene). In total, we provide aug-cc-pVTZ reference vertical excitation energies for 91 excited states of these relatively large systems. We use these reference values to benchmark various wave function methods, i.e., CIS(D), EOM-MP2, CC2, CCSD, STEOM-CCSD, CCSD(T)(a)*, CCSDR(3), CCSDT-3, ADC(2), ADC(2.5), and ADC(3), as well as some spin-scaled variants of both CC2 and ADC(2). These results are compared to those obtained previously on smaller molecules. It turns out that while the accuracy of some methods is almost unaffected by system size, e.g., CIS(D) and CC3, the performance of others can significantly deteriorate as the systems grow, e.g., EOM-MP2 and CCSD, whereas others, e.g., ADC(2) and CC2, become more accurate for larger derivatives.

In the quest for a stable triplet state in small polyaromatic hydrocarbons: an in silico tool for rational design and prediction

Combining the roles of spin frustration and geometry of odd and even numbered rings in polyaromatic hydrocarbons (PAHs), we design small molecules that show exceedingly small singlet-triplet gaps and stable triplet ground states. Furthermore, a computationally efficient protocol with a model spin Hamiltonian is shown to be capable of qualitative agreement with respect to high level multireference calculations and therefore, can be used for fast molecular discovery and screening.

Tuning the HOMO-LUMO Energy Gap of Small Diamondoids Using Inverse Molecular Design

Functionalized diamondoids show great potential as building blocks for various new optoelectronic applications. However, until now, only simple mono and double substitutions were investigated. In this work, we considered up to 10 and 6 sites for functionalization of the two smallest diamondoids, adamantane and diamantane, respectively, in search for diamondoid derivatives with a minimal and maximal HOMO-LUMO energy gap. To this end, the energy gap was optimized systematically using an inverse molecular design methodology based on the best-first search algorithm combined with a Monte Carlo component to escape local optima. Adamantane derivatives were found with HOMO-LUMO gaps ranging from 2.42 to 10.63 eV, with 9.45 eV being the energy gap of pure adamantane. For diamantane, similar values were obtained. The structures with the lowest HOMO-LUMO gaps showed apparent push-pull character. The push character is mainly formed by sulfur or nitrogen dopants and thiol groups, whereas the pull character is predominantly determined by the presence of electron-withdrawing nitro or carbonyl groups assisted by amino and hydroxyl groups via the formation of intramolecular hydrogen bonds. In contrast, maximal HOMO-LUMO gaps were obtained by introducing numerous electronegative groups.

Accurate Ionization Potentials and Electron Affinities of Acceptor Molecules I. Reference Data at the CCSD(T) Complete Basis Set Limit

In designing organic materials for electronics applications, particularly for organic photovoltaics (OPV), the ionization potential (IP) of the donor and the electron affinity (EA) of the acceptor play key roles. This makes OPV design an appealing application for computational chemistry since IPs and EAs are readily calculable from most electronic structure methods. Unfortunately reliable, high-accuracy wave function methods, such as coupled cluster theory with single, double, and perturbative triples [CCSD(T)] in the complete basis set (CBS) limit are too expensive for routine applications to this problem for any but the smallest of systems. One solution is to calibrate approximate, less computationally expensive methods against a database of high-accuracy IP/EA values; however, to our knowledge, no such database exists for systems related to OPV design. The present work is the first of a multipart study whose overarching goal is to determine which computational methods can be used to reliably compute IPs and EAs of electron acceptors. This part introduces a database of 24 known organic electron acceptors and provides high-accuracy vertical IP and EA values expected to be within ±0.03 eV of the true non-relativistic, vertical CCSD(T)/CBS limit. Convergence of IP and EA values toward the CBS limit is studied systematically for the Hartree-Fock, MP2 correlation, and beyond-MP2 coupled cluster contributions to the focal point estimates.

Rational computing of energy levels for organic electronics: the case of 2-benzylidene-1,3-indandiones

To compute ionization potential and electron affinity (HOMO & LUMO levels by tradition), diffuse functions are redundant, while range-separated hybrid DFT is not always better than B3LYP. Solvent modelling is essential and can be cheap with CPCM.

Towards an understanding of the singlet–triplet splittings in conjugated hydrocarbons: azulene investigated by anion photoelectron spectroscopy and theoretical calculations

Electron density overlaps are correlated with singlet triplet splittings: azulene as a test case.