Semiempirical UNO–CAS and UNO–CI: Method and Applications in Nanoelectronics

Tetracene Dimers: A Platform for Intramolecular Down- and Up-conversion

Photon energy conversion can be accomplished in many different ways, including the two opposing manners, down-conversion (i.e., singlet fission, SF) and up-conversion (i.e., triplet-triplet annihilation up-conversion, TTA-UC). Both processes have the potential to help overcome the detailed balance limit of single-junction solar cells. Tetracene, in which the energies of the lowest singlet excited state and twice the triplet excited state are comparable, exhibits both down- and up-conversion. Here, we have designed meta-diethynylphenylene- and 1,3-diethynyladamantyl-linked tetracene dimers, which feature different electronic coupling, to characterize the interplay between intramolecular SF (intra-SF) and intramolecular TTA-UC (intra-TTA-UC) via steady-state and time-resolved absorption and fluorescence spectroscopy. Furthermore, we have used Pd-phthalocyanine as a sensitizer to enable intra-TTA-UC in the two dimers via indirect photoexcitation in the near-infrared part of the solar spectrum. The work is rounded off by temperature-dependent measurements, which outline key aspects of how thermal effects impact intra-SF and intra-TTA-UC in different dimers.

Hybrid Density Functional Tight Binding (DFTB)─Molecular Mechanics Approach for a Low-Cost Expansion of DFTB Applicability

The density functional-based tight binding (DFTB) method has seen a rise in adoption for materials modeling, as it offers significant improvement in scalability with accuracy comparable to the density functional theory (DFT) when good parameterizations exist. The cost reduction in DFTB compared to DFT is achieved by the pre-parameterization of the elements of the Hamiltonian matrix as well as the repulsion potential between all pairs of atoms. Parameterization for new systems with accuracies competitive with DFT in specific applications requires specialized manpower and computational resources. This prevents the application of the DFTB method to systems for which it was not parameterized. In this work, we explore an approach to address the problem of missing parameters of DFTB by modeling the interactions with missing Slater-Koster parameters with an interatomic interaction potential. When the distance between two atoms modeled at the force-field level is sufficiently large, the approach results in accurate structural and electronic properties. The resulting calculation is therefore a hybrid between DFTB and molecular mechanics, a pure DFTB for atoms with a complete set of interatomic parameterizations, and a mix between DFTB and molecular mechanics for atoms with a missing interatomic parameterization. The approach is expected to be particularly useful for hybrid materials and interfaces. The method is tested on the examples of 2D materials, mixed oxides, and a large-scale calculation of an oxide-oxide interface.

Profiling charge transport: A new computational approach

To maintain life, charge transfer processes must be efficient to allow electrons to migrate across distances as large as 30-50 Å within a timescale from picoseconds to milliseconds, and the free-energy cost should not exceed one electron volt. By employing local ionization and local affinity energies, we calculated the pathway for electron and electron-hole transport, respectively. The pathway is then used to calculate both the driving force and the activation energy. The electronic coupling is calculated using configuration interaction procedure. When the charge acceptor is not known, as in oxidative stress, the charge transport terminals are found using Monte-Carlo simulation. These parameters were used to calculate the rate described by Marcus theory. Our approach has been elaborately explained using the famous androstane example and then applied to two proteins: electron transport in azurin protein and hole-hopping migration route from the heme center of cytochrome c peroxidase to its surface. This model gives an effective method to calculate the charge transport pathway and the free-energy profile within 0.1 eV from the experimental measurements and electronic coupling within 3 meV.

Carbon Nanodots from an In Silico Perspective

Carbon nanodots (CNDs) are the latest and most shining rising stars among photoluminescent (PL) nanomaterials. These carbon-based surface-passivated nanostructures compete with other related PL materials, including traditional semiconductor quantum dots and organic dyes, with a long list of benefits and emerging applications. Advantages of CNDs include tunable inherent optical properties and high photostability, rich possibilities for surface functionalization and doping, dispersibility, low toxicity, and viable synthesis (top-down and bottom-up) from organic materials. CNDs can be applied to biomedicine including imaging and sensing, drug-delivery, photodynamic therapy, photocatalysis but also to energy harvesting in solar cells and as LEDs. More applications are reported continuously, making this already a research field of its own. Understanding of the properties of CNDs requires one to go to the levels of electrons, atoms, molecules, and nanostructures at different scales using modern molecular modeling and to correlate it tightly with experiments. This review highlights different in silico techniques and studies, from quantum chemistry to the mesoscale, with particular reference to carbon nanodots, carbonaceous nanoparticles whose structural and photophysical properties are not fully elucidated. The role of experimental investigation is also presented. Hereby, we hope to encourage the reader to investigate CNDs and to apply virtual chemistry to obtain further insights needed to customize these amazing systems for novel prospective applications.

Artificial intelligence-enhanced quantum chemical method with broad applicability

High-level quantum mechanical (QM) calculations are indispensable for accurate explanation of natural phenomena on the atomistic level. Their staggering computational cost, however, poses great limitations, which luckily can be lifted to a great extent by exploiting advances in artificial intelligence (AI). Here we introduce the general-purpose, highly transferable artificial intelligence-quantum mechanical method 1 (AIQM1). It approaches the accuracy of the gold-standard coupled cluster QM method with high computational speed of the approximate low-level semiempirical QM methods for the neutral, closed-shell species in the ground state. AIQM1 can provide accurate ground-state energies for diverse organic compounds as well as geometries for even challenging systems such as large conjugated compounds (fullerene C60) close to experiment. This opens an opportunity to investigate chemical compounds with previously unattainable speed and accuracy as we demonstrate by determining geometries of polyyne molecules-the task difficult for both experiment and theory. Noteworthy, our method's accuracy is also good for ions and excited-state properties, although the neural network part of AIQM1 was never fitted to these properties.

Comprehensive Computational Investigation of the Barton-Kellogg Reaction for Both Alkyl and Aryl Systems

The course of the Barton-Kellogg (BK) reaction for alkyl- and aryl-substituted substrates has been investigated at the DLPNO-CCSD(T)/def2-TZVPP//ωB97X-D/def2-TZVPP level of theory, with results compared to available experimental kinetic data. Through comparison with the unsubstituted parent system, the preference for the formation of 1,3,4-dihydrothiadiazole over the isomeric 1,2,3-dihydrothiadiazole was observed to result from reduced steric repulsion in the relevant transition-state structure. Nitrogen extrusion [retro-(3 + 2)-cycloaddition] from the intermediate dihydrothiadiazole was found to be the rate-determining step. The barrier for this process was, however, significantly lower for aromatic substrates, which is consistent with the difficulty in isolating aryl-substituted dihydrothiadiazoles. The electronic structure of the transient thiocarbonyl ylide was also investigated, highlighting the contradictory results from wave-function theory- and density functional theory-based methods. Correlation of unrestricted natural orbital eigenvalues with previous experimental models suggested that the dipole intermediates possess low diradical character and are therefore considered to be closed-shell species. Exergonic conrotatory electrocyclization of the dipole led to sterically congested thiirane products, even for very bulky systems (di-t-butyl). These results complement the recent work of Mlostoń et al. Finally, DLPNO-CCSD(T)//ωB97X-D was found to be a reliable method for estimating the feasibility of the BK reaction, which should assist experimentalists in the selection of viable substrates.

Propagation of Holes and Electrons in Metal-Organic Frameworks

Charge transport in two zinc metal-organic frameworks (MOFs) has been investigated using periodic semiempirical molecular orbital calculations with the AM1* Hamiltonian. Restricted Hartree-Fock calculations underestimate the band gap using Koopmans theorem (ca. 2 eV compared to the experimental value of 2.8 eV). However, it almost doubles when the constraint on the wave function to remain spin-restricted is removed and the energies of the UHF Natural Orbitals are used. Charge-transport simulations using propagation of the electron- or hole-density in imaginary time allow charge-transport paths and mechanisms to be determined. The calculated relative mobilities in the directions of the three crystal axes agree with experimental expectations, but the absolute values are not reliable using the current technique. Hole-mobility along the crystal c-axis (along the metal stacks) is found to be 13 times higher in the zinc MOF with anthracene linker (Zn-ANMOF-74) than in the other directions, whereas the factor is far smaller (1.7) for electron mobility. Directional preferences are far less distinct in the equivalent structure with phenyl linkers (Zn-MOF-74). The imaginary-time simulation technique does not give quantitative mobilities. The simulations reveal a change in mechanism between the different directions: Coherent polaron migration is observed along the stacks but tunneling hops between them.

Semiempirical Quantum-Chemical Methods with Orthogonalization and Dispersion Corrections

We present two new semiempirical quantum-chemical methods with orthogonalization and dispersion corrections: ODM2 and ODM3 (ODM x). They employ the same electronic structure model as the OM2 and OM3 (OM x) methods, respectively. In addition, they include Grimme's dispersion correction D3 with Becke-Johnson damping and three-body corrections E ABC for Axilrod-Teller-Muto dispersion interactions as integral parts. Heats of formation are determined by adding explicitly computed zero-point vibrational energy and thermal corrections, in contrast to standard MNDO-type and OM x methods. We report ODM x parameters for hydrogen, carbon, nitrogen, oxygen, and fluorine that are optimized with regard to a wide range of carefully chosen state-of-the-art reference data. Extensive benchmarks show that the ODM x methods generally perform better than the available MNDO-type and OM x methods for ground-state and excited-state properties, while they describe noncovalent interactions with similar accuracy as OM x methods with a posteriori dispersion corrections.

Dehydrative π-extension to nanographenes with zig-zag edges

Zig-zag nanographenes are promising candidates for the applications in organic electronics due to the electronic properties induced by their periphery. However, the synthetic access to these compounds remains virtually unexplored. There is a lack in efficient and mild strategies origins in the reduced stability, increased reactivity, and low solubility of these compounds. Herein we report a facile access to pristine zig-zag nanographenes, utilizing an acid-promoted intramolecular reductive cyclization of arylaldehydes, and demonstrate a three-step route to nanographenes constituted of angularly fused tetracenes or pentacenes. The mild conditions are scalable to gram quantities and give insoluble nanostructures in close to quantitative yields. The strategy allows the synthesis of elusive low bandgap nanographenes, with values as low as 1.62 eV. Compared to their linear homologues, the structures have an increased stability in the solid-state, even though computational analyses show distinct diradical character. The structures were confirmed by X-ray diffraction or scanning tunneling microscopy.

Semiempirical Quantum-Chemical Orthogonalization-Corrected Methods: Benchmarks for Ground-State Properties

The semiempirical orthogonalization-corrected OMx methods (OM1, OM2, and OM3) go beyond the standard MNDO model by including additional interactions in the electronic structure calculation. When augmented with empirical dispersion corrections, the resulting OMx-Dn approaches offer a fast and robust treatment of noncovalent interactions. Here we evaluate the performance of the OMx and OMx-Dn methods for a variety of ground-state properties using a large and diverse collection of benchmark sets from the literature, with a total of 13035 original and derived reference data. Extensive comparisons are made with the results from established semiempirical methods (MNDO, AM1, PM3, PM6, and PM7) that also use the NDDO (neglect of diatomic differential overlap) integral approximation. Statistical evaluations show that the OMx and OMx-Dn methods outperform the other methods for most of the benchmark sets.

Semiempirical Quantum-Chemical Orthogonalization-Corrected Methods: Theory, Implementation, and Parameters

Semiempirical orthogonalization-corrected methods (OM1, OM2, and OM3) go beyond the standard MNDO model by explicitly including additional interactions into the Fock matrix in an approximate manner (Pauli repulsion, penetration effects, and core–valence interactions), which yields systematic improvements both for ground-state and excited-state properties. In this Article, we describe the underlying theoretical formalism of the OMx methods and their implementation in full detail, and we report all relevant OMx parameters for hydrogen, carbon, nitrogen, oxygen, and fluorine. For a standard set of mostly organic molecules commonly used in semiempirical method development, the OMx results are found to be superior to those from standard MNDO-type methods. Parametrized Grimme-type dispersion corrections can be added to OM2 and OM3 energies to provide a realistic treatment of noncovalent interaction energies, as demonstrated for the complexes in the S22 and S66×8 test sets.

The Electronic Structure of Amorphous Carbon Nanodots

We have studied hydrogen-passivated amorphous carbon nanostructures with semiempirical molecular orbital theory in order to provide an understanding of the factors that affect their electronic properties. Amorphous structures were first constructed using periodic calculations in a melt/quench protocol. Pure periodic amorphous carbon structures and their counterparts doped with nitrogen and/or oxygen feature large electronic band gaps. Surprisingly, descriptors such as the elemental composition and the number of sp(3)-atoms only influence the electronic structure weakly. Instead, the exact topology of the sp(2)-network in terms of effective conjugation defines the band gap. Amorphous carbon nanodots of different structures and sizes were cut out of the periodic structures. Our calculations predict the occurrence of localized electronic surface states, which give rise to interesting effects such as amphoteric reactivity and predicted optical band gaps in the near-UV/visible range. Optical and electronic gaps display a dependence on particle size similar to that of inorganic colloidal quantum dots.

"Tschitschibabin type biradicals": benzenoid or quinoid?

Our findings provide a better understanding of the discrepancies related to the discussion in the literature on the ground state of Tschitschibabin's hydrocarbon. A series of phenylene bridged bisnitroxide molecules were designed and studied for the comparison. The simple theoretical and experimental methodologies have been developed and utilized to understand the singlet biradicaloids which exist in semi-quinoid form and exhibits characteristics of biradicaloid and quinoid form simultaneously.

Doping the Buckminsterfullerene by Substitution: Density Functional Theory Studies of C59X (X = B, N, Al, Si, P, Ga, Ge, and As)

The heterofullerenes C59X (X = B, N, Al, Si, P, Ga, Ge, and As) were investigated by quantum chemistry calculations based on density functional theory. These hybrid cages can be seen as doping the buckminsterfullerene by heteroatom substitution. The geometrical structures, relative stabilities, electronic properties, vibrational frequencies, dielectric constants, and aromaticities of the doped cages were studied systemically and compared with those of the pristine C60cage. It is found that the doped cages with different heteroatoms exhibit various electronic, vibrational, and aromatic properties. These results imply the possibility to modulate the physical properties of these fullerene-based materials by tuning substitution elements.

Doped polycyclic aromatic hydrocarbons as building blocks for nanoelectronics: a theoretical study

Density functional theory (DFT) and semiempirical UHF natural orbital configuration interaction (UNO-CI) calculations are used to investigate the effect of heteroatom substitution at the central position of a model polycyclic aromatic hydrocarbon. The effects of the substitution on structure, strain, electronic and spectral properties, and aromaticity of the compounds are discussed.

A π-stacked porphyrin-fullerene electron donor-acceptor conjugate that features a surprising frozen geometry

A "frozen" electron donor-acceptor array that bears porphyrin and fullerene units covalently linked through the ortho position of a phenyl ring and the nitrogen of a pyrrolidine ring, respectively, is reported. Electrochemical and photophysical features suggest that the chosen linkage supports both through-space and through-bond interactions. In particular, it has been found that the porphyrin singlet excited state decays within a few picoseconds by means of a photoinduced electron transfer to give the rapid formation of a long-lived charge-separated state. Density functional theory (DFT) calculations show HOMO and LUMO to be localized on the electron-donating porphyrin and the electron-accepting fullerene moiety, respectively, at this level of theory. More specifically, semiempirical molecular orbital (MO) configuration interaction (CI) and unrestricted natural orbital (UNO)-CI methods shed light on the nature of the charge-transfer states and emphasize the importance of the close proximity of donor and acceptor for effective electron transfer.