A multi-descriptor analysis of substituent effects on the structure and aromaticity of benzene derivatives: π-Conjugation versus charge effects.

This work provides a detailed multi-component analysis of aromaticity in monosubstituted (X = CH3, C H 2 - , C H 2 + , NH2, NH-, NH+, OH, O-, and O+) and para-homodisubstituted (X = CH3, CH2, NH2, NH, OH, and O) benzene derivatives. We investigate the effects of substituents using single-reference (B3LYP/DFT) and multireference (CASSCF/MRCI) methods, focusing on structural (HOMA), vibrational (AI(vib)), topological (ELFπ), electronic (MCI), magnetic (NICS), and stability (S0-T1 splitting) properties. The findings reveal that appropriate π-electron-donating and π-electron-accepting substituents with suitable size and symmetry can interact with the π-system of the ring, significantly influencing π-electron delocalization. While the charge factor has a minimal impact on π-electron delocalization, the presence of a pz orbital capable of interacting with the π-electron delocalization is the primary factor leading to a deviation from the typical aromaticity characteristics observed in benzene.

Stability and Reactivity of the Phenalene and Olympicene Isomers.

The phenalene (triangulene) and olympicene molecules belong to the polycyclic aromatic hydrocarbon class, which have attracted substantial technological interest due to their unique electronic properties. Electronic structure calculations serve as a valuable tool in investigating the stability and reactivity of these molecular systems. In the present work, the multireference calculations, namely, the complete active space second-order perturbation theory and multireference averaged quadratic coupled cluster (MR-AQCC), were employed to study the reactivity and stability of phenalene and olympicene isomers, as well as their modified structures where the sp3-carbon at the borders were removed. The harmonic oscillator model of aromaticity (HOMA) and the nucleus-independent chemical shift as geometric and magnetic indexes calculated with density functional theory were utilized to assess the aromaticity of the studied molecules. These indexes were compared with properties such as the excitation energy and natural orbitals occupation. The reactivity analyzed using the HOMA index combined with MR-AQCC revealed the radical character of certain structures as well as the weakening of their aromaticity. Moreover, the results suggest that the removal of sp3-carbon atoms and the addition of hydrogen atoms did not alter the π network and the excitation energies of the phenalene molecules.

Quantification of the Ionic Character of Multiconfigurational Wave Functions: The Qat Diagnostic.

The complete active space self-consistent field (CASSCF) method is a cornerstone in modern excited-state quantum chemistry providing the starting point for most common multireference computations. However, CASSCF, when used with a minimal active space, can produce significant errors (>2 eV) even for the excitation energies of simple hydrocarbons if the states of interest possess ionic character. After illustrating this problem in some detail, we present a diagnostic for ionic character, denoted as Q at, that is readily computed from the transition density. A set of 11 molecules is considered to study errors in vertical excitation energies. State-averaged CASSCF obtains a mean absolute error (MAE) of 0.87 eV for the 34 singlet states considered. We highlight a strong correlation between the obtained errors and the Q at diagnostic, illustrating its power to predict problematic cases. Conversely, using multireference configuration interaction with single and double excitations and Pople's size extensivity correction (MR-CISD+P), excellent results are obtained with an MAE of 0.11 eV. Furthermore, correlations with the Q at diagnostic disappear. In summary, we hope that the presented diagnostic will facilitate reliable and user-friendly multireference computations on conjugated organic molecules.

A comprehensive analysis of charge transfer effects on donor-pyrene (bridge)-acceptor systems using different substituents.

The alternant polycyclic aromatic hydrocarbon pyrene has photophysical properties that can be tuned with different donor and acceptor substituents. Recently, a D (donor)-Pyrene (bridge)-A (acceptor) system, DPA, with the electron donor N,N-dimethylaniline (DMA), and the electron acceptor trifluoromethylphenyl (TFM), was investigated by means of time-resolved spectroscopic measurements (J. Phys. Chem. Lett. 2021, 12, 2226-2231). DPA shows great promise for potential applications in organic electronic devices. In this work, we used the ab initio second-order algebraic diagrammatic construction method ADC(2) to investigate the excited-state properties of a series of analogous DPA systems, including the originally synthesized DPAs. The additionally investigated substituents were amino, fluorine, and methoxy as donors and nitrile and nitro groups as acceptors. The focus of this work was on characterizing the lowest excited singlet states regarding charge transfer (CT) and local excitation (LE) characters. For the DMA-pyrene-TFM system, the ADC(2) calculations show two initial electronic states relevant for interpreting the photodynamics. The bright S1 state is locally excited within the pyrene moiety, and an S2 state is localized ~0.5 eV above S1 and characterized as a donor to pyrene CT state. HOMO and LUMO energies were employed to assess the efficiency of the DPA compounds for organic photovoltaics (OPVs). HOMO-LUMO and optical gaps were used to estimate power conversion and light-harvesting efficiencies for practical applications in organic solar cells. Considering the systems using smaller D/A substituents, compounds with the strong acceptor NO2 substituent group show enhanced CT and promising properties for use in OPVs. Some of the other compounds with small substituents are also found to be competitive in this regard.

Structural stability and the low-lying singlet and triplet states of BN-n-acenes, n = 1-7.

The chemical stability and the low-lying singlet and triplet excited states of BN-n-acenes (n = 1-7) were studied using single reference and multireference methodologies. From the calculations, descriptors such as the singlet-triplet splitting, the natural orbital (NO) occupations and aromaticity indexes are used to provide structural and energetic analysis. The boron and nitrogen atoms form an isoelectronic pair of two carbon atoms, which was used for the complete substitution of these units in the acene series. The structural analysis confirms the effects originated from the insertion of a uniform pattern of electronegativity difference within the molecular systems. The covalent bonds tend to be strongly polarized which does not happen in the case of a carbon-only framework. This effect leads to a charge transfer between neighbor atoms resulting in a more strengthened structure, keeping the aromaticity roughly constant along the chain. The singlet-triplet splitting also agrees with this stability trend, maintaining a consistent gap value for all molecules. The BN-n-acenes molecules possess a ground state with monoconfigurational character indicating their electronic stability. The low-lying singlet excited states have charge transfer character, which proceeds from nitrogen to boron.

Non-Kasha fluorescence of pyrene emerges from a dynamic equilibrium between excited states.

Pyrene fluorescence after a high-energy electronic excitation exhibits a prominent band shoulder not present after excitation at low energies. The standard assignment of this shoulder as a non-Kasha emission from the second-excited state (S2) has been recently questioned. To elucidate this issue, we simulated the fluorescence of pyrene using two different theoretical approaches based on vertical convolution and nonadiabatic dynamics with nuclear ensembles. To conduct the necessary nonadiabatic dynamics simulations with high-lying electronic states and deal with fluorescence timescales of about 100 ns of this large molecule, we developed new computational protocols. The results from both approaches confirm that the band shoulder is, in fact, due to S2 emission. We show that the non-Kasha behavior is a dynamic-equilibrium effect not caused by a metastable S2 minimum. However, it requires considerable vibrational energy, which can only be achieved in collisionless regimes after transitions into highly excited states. This strict condition explains why the S2 emission was not observed in some experiments.

Excitonic and charge transfer interactions in tetracene stacked and T-shaped dimers.

Extended quantum chemical calculations were performed for the tetracene dimer to provide benchmark results, analyze the excimer survival process, and explore the possibility of using long-range-corrected (LC) time-dependent second-order density functional tight-biding (DFTB2) for this system. Ground- and first-excited-state optimized geometries, vertical excitations at relevant minima, and intermonomer displacement potential energy curves (PECs) were calculated for these purposes. Ground-state geometries were optimized with the scaled-opposite-spin (SOS) second-order Møller-Plesset perturbation (MP2) theory and LC-DFT (density functional theory) and LC-DFTB2 levels. Excited-state geometries were optimized with SOS-ADC(2) (algebraic diagrammatic construction to second-order) and the time-dependent approaches for the latter two methods. Vertical excitations and PECs were compared to multireference configuration interaction DFT (DFT/MRCI). All methods predict the lowest-energy S0 conformer to have monomers parallel and rotated relative to each other and the lowest S1 conformer to be of a displaced-stacked type. LC-DFTB2, however, presents some relevant differences regarding other conformers for S0. Despite some state-order inversions, overall good agreement between methods was observed in the spectral shape, state character, and PECs. Nevertheless, DFT/MRCI predicts that the S1 state should acquire a doubly excited-state character relevant to the excimer survival process and, therefore, cannot be completely described by the single reference methods used in this work. PECs also revealed an interesting relation between dissociation energies and the intermonomer charge-transfer interactions for some states.

Dynamics of benzene excimer formation from the parallel-displaced dimer.

Excimers play a key role in a variety of excited-state processes, such as exciton trapping, fluorescence quenching, and singlet-fission. The dynamics of benzene excimer formation in the first 2 ps after S1 excitation from the parallel-displaced geometry of the benzene dimer is reported here. It was simulated via nonadiabatic surface-hopping dynamics using the second-order algebraic diagrammatic construction (ADC(2)). After excitation, the benzene rings take ∼0.5-1.0 ps to approach each other in a parallel-stacked structure of the S1 minimum and stay in the excimer region for ∼0.1-0.4 ps before leaving due to excess vibrational energy. The S1-S2 gap widens considerably while the rings visit the excimer region in the potential energy surface. Our work provides detailed insight into correlations between nuclear and electronic structure in the excimer and shows that decreased ring distance goes along with enhanced charge transfer and that fast exciton transfer happens between the rings, leading to the equal probability of finding the exciton in each ring after around 1.0 ps.

Electronic structure theory gives insights into the higher efficiency of the PTB electron-donor polymers for organic photovoltaics in comparison with prototypical P3HT.

The electron donor poly-thienothiophene-benzodithiophene (PTB) polymer series displays remarkable properties that lead to more efficient bulk heterojunction (BHJ) organic solar cells. In this work, the ground and four excited states (bright S 1 and dark S 2-S 4) of three different members of the PTBn (n = 1, 6, 7) series were studied and compared with the prototypical poly(3-hexylthiophene) (P3HT) donor polymer. Time-dependent density functional theory was employed to investigate oligomers of similar sizes (∼50 Å). Charge alternation electron accumulation and depletion regions of the four transitions are concentrated on the inner units, thereby favoring interaction with the electron acceptor in a BHJ. The bright S 1 transition energies of PTBn are about 0.2 eV lower as compared to P3HT, thereby allowing a better match of their levels with the typical C60-type acceptor moiety in a BHJ. Side chains play a minor role in the electronic spectrum (less than ∼0.1 eV). The most efficient PTB7 transfers more electronic charge from its electron-rich benzodithiophene subunit to its electron-deficient thieno[3,4-b] thiophene subunit as compared to PTB1 and PTB6. We show that the dipolar effect, a partial concentration of negative and positive charges on the different parts of the donor polymer that favors charge separation, is more pronounced in PTBn polymers and typically an order of magnitude larger as compared to P3HT. These effects are conspicuous for the most efficient polymer of the series, PTB7, with its fluorine substituent shown to play a crucial role.

Interplay between Aromaticity and Radicaloid Character in Nitrogen-Doped Oligoacenes Revealed by High-Level Multireference Methods.

Aromaticity is a multivariable concept in organic chemistry that plays a central role for understanding the structure, stability, and reactivity of polycyclic aromatic hydrocarbons (PAHs). Several types of PAHs are characterized as singlet biradicaloid species and their chemical stability is intimately linked to the degree of aromatic character. In this study, theoretically designed routes to tune the biradical character (and thereby its chemical stability) of nitrogen-substituted octacenes have been investigated on the basis of the high-level multireference averaged quadratic coupled-cluster MR-AQCC method necessary for the appropriate description of polyradicaloid systems. The influence of nitrogen centers on the aromaticity of octacene is probed through structural (HOMA) and electron localization (ELF) indices by comparing the N- against NH-doping cases. These analyses reveal that the aromaticity and biradical character of octacene is only slightly affected by replacing one pair of CH groups with N atoms, i.e., by N-doping. However, a significant aromatic stabilization can be obtained when NH-doping is applied at the inner octacene rings; this is also accompanied by an overall decrease of the open-shell character, as evidenced by the gradual quenching of the unpaired electrons and increase in the singlet-triplet splittings when the NH doping groups are moved toward the center of the octacene molecule. Our findings aid in the rational design of new PAH compounds with balanced biradicaloid character and chemical stability which is important, e.g., for practical applications in organic solar cells based on the singlet-fission mechanism.

How to efficiently tune the biradicaloid nature of acenes by chemical doping with boron and nitrogen.

Acenes are fascinating polyaromatic compounds that combine impressive semiconductor properties with an open-shell character by varying their molecular sizes. However, the increasing chemical instabilities related to their biradicaloid structures pose a great challenge for synthetic chemistry. Modifying the π-bond topology through chemical doping allows modulation of the electronic properties of graphene-related materials. In spite of the practical importance of these techniques, remarkably little is known about the basic question - the extent of the radical character created or quenched thereby. In this work, we report a high-level computational study on two acene oligomers doubly-doped with boron and nitrogen, respectively. These calculations demonstrate precisely which the chemical route is in order to either quench or enhance the radical character. Moving the dopants from the terminal rings to the central ones leads to a remarkable variation in the biradicaloid character (and thereby also in the chemical stability). This effect is related to a π-charge transfer involving the dopants and the radical carbon centers at the zigzag edges. This study also provides specific guidelines for a rational design of large polyaromatic compounds with enhanced chemical stability.

Singlet La and Lb Bands for N-Acenes (N = 2-7): A CASSCF/CASPT2 Study.

In this work CASPT2 calculations of polyacenes (from naphthalene to heptacene) were performed to find a methodology suitable for calculations of the absorption spectra, in particular of the La (B2u state) and Lb (B3u state) bands, of more extended systems. The effect of the extension of the active space and of freezing σ orbitals was investigated. The MCSCF excitation energy of the B2u state is not sensitive to the size of the active space used. However, the CASPT2 results depend strongly on the amount of σ orbitals frozen reflecting the ionic character of the B2u state. On the other hand, the excitation energies of the B3u state are much more sensitive to the size of the active space used in the calculations reflecting its multiconfigurational character. We found a good agreement with experimental data for both bands by including 14 electrons in 14 π orbitals in the active space followed by the CASPT2(14,14) perturbation scheme in which both σ and π orbitals are included.

Thermochemical and Kinetics of Hydrazine Dehydrogenation by an Oxygen Atom in Hydrazine-Rich Systems: A Dimer Model.

The kinetics of the reaction of N2H4 with oxygen depends sensitively on the initial conditions used. In oxygen-rich systems, the rate constant shows a conventional positive temperature dependence, while in hydrazine-rich setups the dependence is negative in certain temperature ranges. In this study, a theoretical model is presented that adequately reproduces the experimental results trend and values for hydrazine-rich environment, consisting of the hydrogen abstraction from the hydrazine (N2H4) dimer by an oxygen atom. The thermochemical properties of the reaction were computed using two quantum chemical approaches, the coupled cluster theory with single, double, and noniterative triple excitations (CCSD(T)) and the M06-2X DFT approach with the aug-cc-pVTZ and the maug-cc-pVTZ basis sets, respectively. The kinetic data were calculated with the improved canonical variational theory (ICVT) using a dual-level methodology to build the reaction path. The tunneling effects were considered by means of the small curvature tunneling (SCT) approximation. Potential wells on both sides of the reaction ((N2H4)2 + O → N2H4···N2H3 + OH) were determined. A reaction path with a negative activation energy was found leading, in the temperature range of 250-423 K, to a negative dependence of the rate constant on the temperature, which is in good agreement with the experimental measurements. Therefore, the consideration of the hydrazine dimer model provides significantly improved agreement with the experimental data and should be included in the mechanism of the global N2H4 combustion process, as it can be particularly important in hydrazine-rich systems.

The electronic states of a double carbon vacancy defect in pyrene: a model study for graphene.

The electronic states occurring in a double vacancy defect for graphene nanoribbons have been calculated in detail based on a pyrene model. Extended ab initio calculations using the MR configuration interaction (MRCI) method have been performed to describe in a balanced way the manifold of electronic states derived from the dangling bonds created by initial removal of two neighboring carbon atoms from the graphene network. In total, this study took into account the characterization of 16 electronic states (eight singlets and eight triplets) considering unrelaxed and relaxed defect structures. The ground state was found to be of (1)Ag character with around 50% closed shell character. The geometry optimization process leads to the formation of two five-membered rings in a pentagon-octagon-pentagon (5-8-5) structure. The closed shell character increases thereby to ∼70%; the analysis of unpaired density shows only small contributions confirming the chemical stability of that entity. For the unrelaxed structure the first five excited states ((3)B3g, (3)B2u, (3)B1u, (3)Au and (1)Au) are separated from the ground state by less than 2.5 eV. For comparison, unrestricted density functional theory (DFT) calculations using several types of functionals have been performed within different symmetry subspaces defined by the open shell orbitals. Comparison with the MRCI results gave good agreement in terms of finding the (1)Ag state as a ground state and in assigning the lowest excited states. Linear interpolation curves between the unrelaxed and relaxed defect structures also showed good agreement between the two classes of methods opening up the possibilities of using extended nanoflakes for multistate investigations at the DFT level.

Absorption and fluorescence spectra of poly(p-phenylenevinylene) (PPV) oligomers: an ab initio simulation.

The absorption and fluorescence spectra of poly(p-phenylenevinylene) (PPV) oligomers with up to seven repeat units were theoretically investigated using the algebraic diagrammatic construction method to second order, ADC(2), combined with the resolution-of-the-identity (RI) approach. The ground and first excited state geometries of the oligomers were fully optimized. Vertical excitation energies and oscillator strengths of the first four transitions were computed. The vibrational broadening of the absorption and fluorescence spectra was studied using a semiclassical nuclear ensemble method. After correcting for basis set and solvent effects, we achieved a balanced description of the absorption and fluorescence spectra by means of the ADC(2) approach. This fact is documented by the computed Stokes shift along the PPV series, which is in good agreement with the experimental values. The experimentally observed band width of the UV absorption and fluorescence spectra is well reproduced by the present simulations showing that the nuclear ensemble generated should be well suitable for consecutive surface hopping dynamics simulations.

Hydrogen abstraction from the hydrazine molecule by an oxygen atom.

Thermochemical and kinetics properties of the hydrogen abstraction from the hydrazine molecule (N2H4) by an oxygen atom were computed using high-level ab initio methods and the M06-2X DFT functional with aug-cc-pVXZ (X = T, Q) and maug-cc-pVTZ basis sets, respectively. The properties along the reaction path were obtained using the dual-level methodology to build the minimum energy path with the potential energy surface obtained with the M06-2X method and thermochemical properties corrected with the CCSD(T)/CBS//M06-2X/maug-cc-pVTZ results. The thermal rate constants were calculated in the framework of variational transition-state theory. Wells on both sides of the reaction (reactants and products) were found and considered in the chemical kinetics calculations. Additionally, the product yields were investigated by means of a study of the triplet and singlet surfaces of the N2H4 + O → N2H2 + H2O reaction.

A multireference configuration interaction study of the photodynamics of nitroethylene.

Extended multireference configuration interaction with singles and doubles (MR-CISD) calculations of nitroethylene (H2C═CHNO2) were carried out to investigate the photodynamical deactivation paths to the ground state. The ground (S0) and the first five valence excited electronic states (S1-S5) were investigated. In the first step, vertical excitations and potential energy curves for CH2 and NO2 torsions and CH2 out-of-plane bending starting from the ground state geometry were computed. Afterward, five conical intersections, one between each pair of adjacent states, were located. The vertical calculations mostly confirm the previous assignment of experimental spectrum and theoretical results using lower-level calculations. The conical intersections have as main features the torsion of the CH2 moiety, different distortions of the NO2 group and CC, CN, and NO bond stretchings. In these conical intersections, the NO2 group plays an important role, also seen in excited state investigations of other nitro molecules. Based on the conical intersections found, a photochemical nonradiative deactivation process after a π-π* excitation to the bright S5 state is proposed. In particular, the possibility of NO2 release in the ground state, an important property in nitro explosives, was found to be possible.