Benchmarking Ground-State Geometries and Vertical Excitation Energies of a Selection of P-Type Semiconducting Molecules with Different Polarity

Theoretical Study on the Optoelectronic Properties of Merocyanine-Dyes

Merocyanines, as prototypes of highly polar π-conjugated molecules, have been intensively investigated for their self-assembly and optoelectronic properties, both experimentally and theoretically. However, an accurate description of their structural and electronic properties remains challenging for quantum-chemical methods. We assessed several theoretical approaches, TD-DFT, GW-BSE, STEOM-DLPNO-CCSD, and CASSCF/NEVPT2-FIC for their reliability in reproducing optoelectronic properties of a series of donor/acceptor (D/A) merocyanines, focusing on the first excitation energy. Additionally, we tested an all-electron perturbative method based on time-dependent coupled-perturbed density functional theory, denoted as TDCP-DFT. Particular focus was set on direct and indirect solvent effects, which affect excited-state energies by electrostatic interaction and molecular geometry. The molecular configuration space was sampled at the semiempirical tight-binding level. Our results corroborate previous investigations, showing that the S0 - S1 excitation energy strongly depends on the merocyanine molecular structure and the dielectric constant of the solvent. We found significant effects of the polar solution environment on the geometry of the merocyanines, which strongly affect the calculated excitation energies. Taking these effects into account, the best agreement between calculated and measured excitation energies was obtained with TDCP-DFT and GW-BSE. We also calculated excitation energies of molecular crystals at the TDCP-DFT level and compared the results to the corresponding monomers.

Quantum mechanical assessment on the optical properties of capsanthin conformers

As optical properties, the ultraviolet-visible (UV-Vis) absorption spectra of capsanthin-based red natural dye are a decisive parameter for their usage in various applications. Thus, accurately predicting the maximum UV-Vis wavelength (λ max ) values is critical in designing dye-conjugated material. Extensive metadynamics simulations were carried out to generate capsanthin conformers at various levels of the extended tight-binding method. Benchmarking the time-dependent density-functional theory (TD-DFT) methods help understand the results of a particular functional and allows a comparison between results obtained with different functional. The long-range correction (LC) scheme in LC-TD-DFT-D4/ωB97X/def2-SVP has been found to reproduce the experimentalλ max , and exhibited the effect of conformational changes to the calculated wavelengths. On the other hand, an inexpensive yet efficient LC-TD-DFTB method reproduced the experimentalλ max insensitive to conformational changes.

Computational Approaches for Organic Semiconductors: From Chemical and Physical Understanding to Predicting New Materials

While a complete understanding of organic semiconductor (OSC) design principles remains elusive, computational methods─ranging from techniques based in classical and quantum mechanics to more recent data-enabled models─can complement experimental observations and provide deep physicochemical insights into OSC structure-processing-property relationships, offering new capabilities for in silico OSC discovery and design. In this Review, we trace the evolution of these computational methods and their application to OSCs, beginning with early quantum-chemical methods to investigate resonance in benzene and building to recent machine-learning (ML) techniques and their application to ever more sophisticated OSC scientific and engineering challenges. Along the way, we highlight the limitations of the methods and how sophisticated physical and mathematical frameworks have been created to overcome those limitations. We illustrate applications of these methods to a range of specific challenges in OSCs derived from π-conjugated polymers and molecules, including predicting charge-carrier transport, modeling chain conformations and bulk morphology, estimating thermomechanical properties, and describing phonons and thermal transport, to name a few. Through these examples, we demonstrate how advances in computational methods accelerate the deployment of OSCsin wide-ranging technologies, such as organic photovoltaics (OPVs), organic light-emitting diodes (OLEDs), organic thermoelectrics, organic batteries, and organic (bio)sensors. We conclude by providing an outlook for the future development of computational techniques to discover and assess the properties of high-performing OSCs with greater accuracy.

One- and two-electron reduction of triarylborane-based helical donor-acceptor compounds

One-electron chemical reduction of 10-(dimesitylboryl)-N,N-di-p-tolylbenzo[c]phenanthrene-4-amine (3-B(Mes)2-[4]helix-9-N(p-Tol)2) 1 and 13-(dimesitylboryl)-N,N-di-p-tolyldibenzo[c,g]phenanthrene-8-amine (3-B(Mes)2-[5]helix-12-N(p-Tol)2) 2 gives rise to monoanions with extensive delocalization over the annulated helicene rings and the boron p z orbital. Two-electron chemical reduction of 1 and 2 produces open-shell biradicaloid dianions with temperature-dependent population of the triplet states due to small singlet-triplet gaps. These results have been confirmed by single-crystal X-ray diffraction, EPR and UV/vis-NIR spectroscopy, and DFT calculations.

A combined Monte Carlo/DFT approach to simulate UV-vis spectra of molecules and aggregates: Merocyanine dyes as a case study

The combination of a Monte Carlo (MC) sampling of the configurational space with time dependent-density functional theory (TD-DFT) to estimate vertical excitations energies has been applied to compute the absorption spectra of a family of merocyanine dyes in both their monomeric and dimeric forms. These results have been compared to those obtained using a static DFT/TD-DFT approach as well as to the available experimental spectra. Though suffering of the limitations related to the use of DFT and TD-DFT for this type of systems, our data clearly show that the classical MC sampling provides a suitable alternative to classical molecular dynamics to explore the structural flexibility of these donor-acceptor (D-π-A) chromophores enabling a realistic description of the potential energy surface of both their monomers and aggregates (here dimers) and thus of their spectra. Overall, the combination of MC sampling with quantum mechanics (TD-DFT) calculations, carried out in implicit dioxane solvent on random snapshots, provides a workable compromise to solve the combined challenge of accuracy and time-consuming problem not only for merocyanines momers, but also for their dimers, up to now less investigated. Indeed, the simulated absorption spectra fairly agree with the experimental ones, suggesting the general reliability of the method.

Twisting versus Delocalization in CAAC- and NHC-Stabilized Boron-Based Biradicals: The Roles of Sterics and Electronics

Twisted boron-based biradicals featuring unsaturated C2 R2 (R=Et, Me) bridges and stabilization by cyclic (alkyl)(amino)carbenes (CAACs) were recently prepared. These species show remarkable geometrical and electronic differences with respect to their unbridged counterparts. Herein, a thorough computational investigation on the origin of their distinct electrostructural properties is performed. It is shown that steric effects are mostly responsible for the preference for twisted over planar structures. The ground-state multiplicity of the twisted structure is modulated by the σ framework of the bridge, and different R groups lead to distinct multiplicities. In line with the experimental data, a planar structure driven by delocalization effects is observed as global minimum for R=H. The synthetic elusiveness of C2 R2 -bridged systems featuring N-heterocyclic carbenes (NHCs) was also investigated. These results could contribute to the engineering of novel main group biradicals.

The ortho-benzyne cation is not planar

A recent review on the photoionisation of the C₆H₄ isomer ortho-benzyne suggests that bands reported in earlier photoelectron spectra might be due to side products or contaminations, while computations raise doubts, whether the cation has a planar geometry. We therefore reinvestigate the photoionisation of ortho-benzyne, generated by pyrolysis from benzocyclobutenedione, by photoion mass-selected threshold photoelectron (ms-TPE) spectroscopy using synchrotron radiation. The experiments are accompanied by a theoretical study that investigates the structure of the ortho-benzyne cation systematically as a function of the computational method, up to CASPT2(11,14) ab initio computations. Our study leads to a re-evaluation of the ionisation energy of ortho-benzyne. It reveals that the ortho-benzyne cation has indeed a twisted C₂ geometry rather than a C_2v structure. A vertical ionisation energy IE_vert of 9.77 eV and an adiabatic ionisation energy of IE_ad = 9.56 eV are computed for ortho-benzyne. A Franck-Condon simulation of the photoelectron spectrum based on the CASPT2 results and including three electronic states of the cation is in agreement with the experiment and yields IE_ad = 9.51 eV (+50 meV/-100 meV). Since this value is in contrast with previous work, the ionisation energy has to be revised based on our study. Computational methods based on density functional theory give a reasonable description of the cationic ground state, but fail for the corresponding excited electronic states that are indispensible for a proper assignment of the photoelectron spectrum.

The electronic character of PTCDA thin films in comparison to other perylene-based organic semi-conductors: ab initio-, TD-DFT and semi-empirical computations of the opto-electronic properties of large aggregates

Perylene-based compounds are promising materials for opto-electronic thin film devices but despite intense investigations, important details of their electronic structure are still under debate. For perylene-3,4,9,10-tetracarboxylic dianhydrid (PTCDA), the theoretical models predict a different relative energetic order of Frenkel and Charge Transfer (CT) states. Additionally, while one model indicates strong differences between PTCDA on one hand and other perylene-based compounds on the other, recent ab initio computations indicate electronic properties of all perylene-based compounds to resemble each other. Finally, the models disagree about the amount of mixing between CT and Frenkel states. Definitive answers to these questions are difficult because the approaches use various approximations. Up to date, the ab initio based methods employ rather small model systems and neglect environmental effects. In the present work, we improve our former approach by analyzing the effects of the various simplifications. In more detail, we increase the size of the model systems, include environmental effects and investigate the influence of exciton-phonon couplings on the absorption spectrum. The computations for larger aggregates were performed with the ZINDO/S approach, because benchmark computations show that it provides accurate vertical excitation energies for Frenkel, as well as CT states.

Electronic Structure and Excited-State Dynamics of an Arduengo-Type Carbene and its Imidazolone Oxidation Product

We describe an investigation of the excited-state dynamics of isolated 1,3-di-tert-butyl-imidazoline-2-ylidene (tBu₂ Im, C₁₁ H₂₀ N₂ , m/z=180), an Arduengo-type carbene, by time- and frequency-resolved photoionization using a picosecond laser system. The energies of several singlet and triplet excited states were calculated by time-dependent density functional theory (TD-DFT). The S₁ state of the carbene deactivates on a 100 ps time scale possibly by intersystem crossing. In the experiments we observed an additional signal at m/z=196, that was assigned to the oxidation product 1,3-di-tert-butyl-imidazolone, tBu₂ ImO. It shows a well-resolved resonance-enhanced multiphoton ionization (REMPI) spectrum with an origin located at 36951 cm^-1 . Several low-lying vibrational bands could be assigned, with a lifetime that depends strongly on the excitation energy. At the origin the lifetime is longer than 3 ns, but drops to 49 ps at higher excess energies. To confirm formation of the imidazolone we also performed experiments on benzimidazolone (BzImO) for comparison. Apart from a redshift for BzImO the spectra of the two compounds are very similar. The TD-DFT values display a very good agreement with the experimental data.

The dimer-approach to characterize opto-electronic properties of and exciton trapping and diffusion in organic semiconductor aggregates and crystals

We present the recently developed dimer approach which seems to include all main effects determining the photo-physics of organic semiconductor aggregates.

A theoretical description of charge reorganization energies in molecular organic P-type semiconductors

Charge transport properties of materials composed of small organic molecules are important for numerous optoelectronic applications. A material's ability to transport charges is considerably influenced by the charge reorganization energies of the composing molecules. Hence, predictions about charge-transport properties of organic materials deserve reliable statements about these charge reorganization energies. However, using density functional theory which is mostly used for the predictions, the computed reorganization energies depend strongly on the chosen functional. To gain insight, a benchmark of various density functionals for the accurate calculation of charge reorganization energies is presented. A correlation between the charge reorganization energies and the ionization potentials is found which suggests applying IP-tuning to obtain reliable values for charge reorganization energies. According to benchmark investigations with IP-EOM-CCSD single-point calculations, the tuned functionals provide indeed more reliable charge reorganization energies. Among the standard functionals, ωB97X-D and SOGGA11X yield accurate charge reorganization energies in comparison with IP-EOM-CCSD values. © 2016 Wiley Periodicals, Inc.