Synthesis and Combined Experimental and Theoretical Characterization of Dihydro-tetraaza-acenes.

We present a combined experimental and theoretical study of electronic and optical properties of dihydro-tetraaza-acenes (DHTAn). Using solvent-free condensation, we are able to synthesize not only DHTA5 but also the longer DHTA6 and DHTA7 molecules. We then investigate their gas-phase electronic structures by means of ab initio density functional calculations employing an optimally tuned range-separated hybrid functional. By comparing with the parent linear oligoacenes (nA) and based on computed ionization potentials and electron affinities, we predict DHTAn molecules to be more stable than acenes of the same length, where we expect DHTAn molecules to be persistent at least up to n = 7 rings. We further exploit the analogy with nA by analyzing the entire intramolecular π-band structure of the DHTAn molecules. This clearly reveals that the additional two electrons donated by the dihydropyrazine group are delocalized over the entire molecule and contribute to its π-electron system. As a consequence, the symmetry of the frontier orbitals of DHTAn differs from that of the parent nA molecule. This also affects the UV-vis absorption spectra which have been measured for DHTA5, 6, and 7 dissolved in dimethyl sulfoxide and analyzed by means of excited state calculations within a time-dependent density functional theory framework.

Introduction

In the field of organic electronics, pentacene (5A) is presumably one of the most investigated organic molecules covered by a vast literature ranging from device physics[1,2] over to fundamental studies on electronic[3−7] and optical properties.[8−11] However, pentacene suffers from long-term stability issues because of its photo-oxidation and low thermal stability.[12] To overcome these limitations and to provide alternative molecules with improved film morphology and stability, nitrogen-containing derivatives of pentacene, N-heteroacenes, have been shown to be promising.[13−19] The synthesis of larger heteroacenes and the investigation of their properties have already a long-standing history.[16] As early as in 1901, Hinsberg[20] has presented a comprehensive investigation of such multiunit nitrogeneous ring systems, thereby discovering interesting properties regarding their chemical stability or fluorescence properties. In the 1960s, Kummer and Zimmermann have systematically studied the electronic properties of linear diaza- and tetra-azaacenes compared with that of the parent oligoacenes.[21] More recently, Nuckolls and co-workers have demonstrated the fabrication of a thin-film transistor based on a hydrogenated diaza-pentacene,[14] whereas Winkler and Houk have discussed the usage of nitrogen-rich oligoacenes as n-channel transistors based on density functional calculations.[15]

We decided to focus on 5-14-dihydro-5,7,12,14-tetraazapentacene derivatives (DHTAn) depicted in Scheme because of two fundamental features inherent to their chemical structure which offer a number of opportunities to manipulate and control the electronic properties, the stability, and the supramolecular arrangement in the solid.[22] First, these N-heteroacenes can be viewed as constituted of dihydropyrazine units which add two more electrons to the π system and thereby stabilize the molecule.[23] Second, the presence of H-donor (dihydropyrazine, N–H) and H-acceptor sites (pyrazine, N) leads to well-organized intermolecular H-bonding networks, forming a stable arrangement in the solid state.[24] The first documented molecule of that sort, 5-14-dihydro-5,7,12,14-tetraazapentacene (DHTA5), has in fact been synthesized as early as 1890[25] and also been investigated in the seminal work by Hinsberg who called this molecule homofluorindin.[20] Its structural and electrochemical properties have already been reported,[26,27] and its electronic band structure has been revealed by ultraviolet photoemission spectroscopy.[28] Driven by new efficient synthesis routes,[22] its fundamental properties and its potential for organic devices have been assessed.[17,18,23,29−31] However, extended DHTAn for n > 5 have not been reported so far where the presence of extra phenyl rings should affect the stability of their oxidized and/or reduced forms. Also, a systematic study on the electronic structure of the whole DHTAn family of molecules using a state-of-the-art and predictive ab initio electronic structure method is still lacking.

Scheme 1

Chemical Structures of the Investigated Molecules

The left column shows representatives of the oligoacene series denoted as nA, while the right column depicts dihydro-tetraaza-acenes (DHTAn) with n = 5, 6, 7 rings.

In this work, extended DHTAn have been synthesized (n = 6, 7) and fully characterized. Using a combined experimental and theoretical approach, we investigate their electronic and optical properties. Theoretically, we perform density functional theory (DFT) calculations for the gas-phase oligoacene (nA) and dihydro-tetraaza-acene (DHTAn) molecules. In addition to a standard exchange–correlation functional, a generalized gradient approximation (GGA) as parameterized by Perdew, Burke, and Ernzerhof (PBE),[32] we employ a range-separated hybrid (RSH) density functional, with an asymptotically exact and short-range fractional Fock exchange, where the parameters are tuned nonempirically, per system, on the basis of exact physical constraints, also known as optimally tuned RSH (OT-RSH).[33−36] This method leads to highly reliable values for the ionization potentials (IPs) and the electron affinities (EAs) for DHTAn which we compare to the parent nA series. Moreover, we also present a comparison of the π bands of DHTAn with those of nA with an emphasis on the changes induced by the pyrazine and dihydropyrazine groups. We envision that the insights gained by our study will be not only important for further exploiting DHTAn molecules for various applications but also useful for designing other novel N-heteroacenes.

Experimental Section

Synthesis and Characterization of DHTA6

A well-mixed solid of phenazine-2,3-diol (m = 800 mg, 3.72 mmol, 1 equiv), 2,3-diaminonaphthalene (m = 626 mg, 3.96 mmol, 1.05 equiv), and benzoic acid (m = 4.9 g, 40.0 mmol, 10 equiv) was heated to 400 °C for 10 min. Then, EtOH was added to the hot mixture. This dark-colored suspension in EtOH (v = 50 mL) was sonicated (at 60 °C, 1 h), filtrated, and washed with dimethyl sulfoxide (DMSO) and acetone. The desired product was afforded after 3 times sublimations under vacuum (8 mbar, 450 °C, 15 min) as a dark violet powder (m = 712 mg, 2.31 mmol, 58% yield). 1H NMR (400 MHz, DMSO-d6): δ (ppm) = 10.17 (br s, 2H), 7.78–7.76 (m, 2H), 7.53–7.48 (m, 4H), 7.14–7.12 (m, 2H), 6.86 (s, 2H), 6.63 (s, 2H). Solid 13C NMR (100 MHz): δ (ppm) = 145.3, 144.1, 138.8, 130.7, 128.6, 123.0, 108.2, 100.8, 100.0. HRMS (MALDI-TOF) m/z: [M + H]+ for C22H15N4+ calcd, 335.1291; found, 335.1290; error < 1 ppm.

Synthesis and Characterization of DHTA7

A well-mixed solid of 2,5-dihydroxy-1,4-benzoquinone (m = 445 mg, 3.18 mmol, 1 equiv), 2,3-diaminonaphthalene (m = 1020 mg, 6.45 mmol, 2.03 equiv), and benzoic acid (m = 7.2 g, 59.0 mmol, 18 equiv) was heated to 360 °C for 15 min. Then, EtOH was added to the hot mixture. This dark-colored suspension in EtOH (v = 50 mL) was sonicated (at 60 °C, 1 h) and filtrated thrice. The desired product was afforded after 3 times sublimations under vacuum (8 mbar, 450 °C, 15 min) as a dark violet powder (m = 745 mg, 1.94 mmol, 61% yield). 1H NMR (400 MHz, DMSO-d6): δ (ppm) = 10.55 (br s, 2H), 8.37 (s, 2H), 8.03–8.00 (m, 2H), 7.57–7.55 (m, 2H), 7.43–7.41 (m, 2H), 7.19–7.17 (m, 2H), 7.00 (s, 2H), 6.66 (s, 2H). Solid 13C NMR (100 MHz): δ (ppm) = 146.0, 138.3, 132.2, 126.1, 123.1, 108.2, 99.1. HRMS (MALDI-TOF) m/z: [M + H]+ for C26H17N4+ calcd, 385.1447; found, 385.1447; error < 2 ppm.

Computational Details

Within the generalized Kohn–Sham formalism,[35,37] RSH functionals consist of a partition of the Coulomb interaction as[38−40]Here, r is the interelectron coordinate, erf is the error function, and α, β, and γ are parameters, which in principle may be freely chosen or determined empirically.[38] The first term in eq is treated within Hartree–Fock (HF) theory, while the second one is treated within a standard semilocal (sl) approximation.[41] The range-separation parameter γ controls which of the two terms dominate at a given distance. It has been repeatedly shown that γ can be strongly system- and size-dependent;[33,35,42−44] thus, one uniform value for γ is not sufficient. Therefore, we prefer to determine it separately for each system by enforcing a nonempirical condition: we make use of the IP theorem[45−48] and choose γ such that the difference between the energy of the highest occupied molecular orbital (HOMO) level and the IP is minimized. Specifically, this is accomplished by minimizing the following target function[42]Here, the εγ,α are the HOMOs of the (N + i) electron molecular systems and i is an integer number. The IPγ,α(N + i) are the corresponding IPs, which are calculated from the total energy difference between the (N + i) electron and the N + i – 1 electron system. In eq , when i = 0, the HOMO of the neutral system is adjusted with the IP, whereas when i = 1, the IP of the anionic system, that is, the EA of the neutral system, is considered, such that its difference from the lowest unoccupied molecular orbital (LUMO) eigenvalue is minimized.[35] Note that γ has to be tuned for each choice of α separately.

As can be seen from eq , the parameters α and β control the behavior of the Fock term for r → 0 and r → ∞, leading to and , respectively. The asymptotic behavior of the xc functional for r → ∞ has been shown to be crucial for obtaining accurate HOMO–LUMO gaps and thereby introduces another constraint on the parameters. For a gas-phase molecule, enforcing α + β = 1 leads to the correct asymptotic behavior.[49] As a consequence, the semilocal contribution in the long range vanishes and α controls the amount of nonlocal Fock exchange in the short range. This can be seen in the expression for the exchange–correlation energy of the RSH functional, obtained with this constraint[34]where sl,x and sl,c denote semilocal exchange and correlation, respectively, and HF,x denotes nonlocal Fock exchange.

We have also simulated the optical absorption spectrum of molecules in solution. To this end, we have taken into account the dielectric constant of the solvent, DMSO, for which we take an experimental literature value of ε = 2.18. Enforcing the correct asymptotic behavior of the exchange potential in the solvent, , requires to choose . Thus, the exchange–correlation energy is generalized to[36]Note that eq reduces to eq in the gas-phase case where ε = 1. Using the functional form of eq , we have computed the spin-restricted vertical excitation energies within the framework of time-dependent DFT (TD-DFT) using the Casida formalism.[50]

All electronic structure calculations have been performed using the ab initio computational chemistry software package NWChem.[51] The Gaussian-based DFT module has been used with a cc-PVTZ basis set. For the computations, we have utilized the standard PBE exchange–correlation functional[32] as well as the long-range-corrected RSH functional (LC-wPBE). Note that the preceding geometry relaxations have been performed using the PBE functional and all molecular geometries are found to be planar.

Results

Synthesis

DHTA5 was prepared according to the procedure reported in the literature[24] using solvent-free condensation between the commercially available 2,5-dihydroxy-p-benzoquinone 1 and o-diaminobenzene 2 (2 equiv) (Scheme ). This reaction occurs at high temperature in the presence of excess PhCO2H as the catalyst to afford the target pentacycle that could be easily isolated by filtration in 80% yield. DHTA6 could be isolated as a violet solid in 58% yield by reacting similarly the commercially available 2,3-diaminonaphthalene 4 and 2,3-dihydroxyphenazine 3, which could be obtained in a previous work from 1 and 2 (1 equiv).[22] The same procedure using precursors 1 and 4 (2 equiv) leads to the formation of DHTA7 in 61% yield (see the Experimental Section).

Scheme 2

Synthesis of DHTA5,[24] DHTA6, and DHTA7

Influence of Short-Range Exchange

For our gas-phase DFT calculations, we have tuned the range-separation parameter γ by minimizing the target function defined in eq by taking into account, both, the neutral molecule (i = 0) and the anion (i = 1). To judge the influence of α, the fraction of short-range HF exchange, on the orbital energies, we have varied α from 0 to 0.5 and optimally tuned γ for each value of α. The results for DHTA5 are plotted in Figure . First, we notice that compared to the GGA-PBE result of 1.77 eV, the HOMO–LUMO gap is considerably widened to 5.60 eV in a symmetric manner. The OT-RSH value of the gap, in fact also the absolute HOMO and the LUMO position, is quite insensitive to the value of α. This is because the IP and the EA are mainly governed by the asymptotic behavior of the respective orbital that depends on the long-range behavior of the exchange–correlation potential which is independent of α. As expected, the optimal range-separation parameter γopt reduces with increasing α. Such a behavior can be rationalized because the range 1/γopt at which full Fock exchange sets in can be extended to larger distances if the amount of the Fock exchange at short range, governed by α, is increased.[34,36]

Figure 1

(a) Orbital energies of DHTA5 computed within PBE-GGA and the OT-RSH functional for α-values between 0 and 0.5; the corresponding optimal γ-values are also shown (right axis). (b) Kohn–Sham orbitals resulting from the OT-RSH calculations for α = 0.25; the labels are explained in the text.

We have color-coded the orbital energies according to the symmetry of the respective orbitals that are depicted in panel (b) of Figure . All π orbitals with no node in the x,z-plane, denoted as π, and the respective orbital energies are plotted in red, while π orbitals with a node in the x,z-plane, denoted as π, are plotted in orange. Here, m counts the nodes along the long molecular axis (x-axis). Orbitals of σ symmetry are shown in green, specifically two states denoted as σa and σb which are centered around the pyrazine groups of the molecule. In contrast to the π states, these orbitals exhibit a pronounced dependence on the amount of short-range HF exchange. This is the typical behavior for comparably localized states resulting from the growing compensation of the self-interaction error with increasing α.[36] Note that for the remaining part of the paper, we have set α = 0.25 because this value has been shown to yield the best agreement with quasiparticle calculations for other pentacene derivatives[36] as well as for the organic pigment PTCDA.[34]

It is worth mentioning that the two types of π orbitals introduced in Figure , namely, π and π, are sometimes also termed linking and apex π bands, respectively, referring to the position of the carbon atoms with the strongest contribution.[52] For instance, the HOMO of DHTA5, denoted as π6,0, has six nodes along the long molecular axis and it belongs to the linking band. Interestingly, the HOMO of DHTA5 has the same nodal structure as the LUMO of pentacene.[53] This is because DHTA5 has two π electrons more compared to pentacene owing to the two N–H groups which each contribute one more π electron. Thus, overall DHTA5 exhibits 24 π electrons which lead to 12 occupied states, 5 of which are in the apex band (π0,1, π1,1, ..., π4,1) and 7 in the linking band (π0,0, π1,0, ..., π6,0).

Frontier Orbitals

Naively, one may conclude that 24 π electrons for DHTA5, as compared to 22 π electrons for 5A, may destabilize DHTA5 relative to 5A following the so-called 4n Hückel rule.[23] This is, however, not the case as can be seen from a comparison of the frontier orbitals, that is, the HOMO and LUMO energies, between DHTAn and the parent nA molecules which is the focus of this subsection.

In Figure , we have depicted the IP and the EA, computed from the OT-RSH functional with α = 0.25, for oligoacenes from naphthalene (n = 2) up to heptacene (n = 7) as black triangles and for the DHTAn for n = 5, 6, 7, 9, 11 as red circles. For the oligoacenes up to pentacene, we have also included experimental data (blue plus signs) for the IP[54] and the EA.[55] Both, experiment and calculation, nicely show the reduction of the IP and the increase in EA, respectively, in the series from naphthalene (n = 2), anthracene (n = 3), tetracene (n = 4), and pentacene (n = 5). The OT-RSH results slightly underestimate IP with a maximum deviation of 0.3 eV compared to the experiment, while the EA values all agree within ±0.15 eV. The further decrease in IP (increase in EA) computed for hexacene (n = 6) and heptacene (n = 7) is in fact the reason for the well-known chemical instability of linear acenes longer than pentacene.

Figure 2

(a) IPs and EAs for nA and DHTAn vs number of rings n on a reciprocal axis. Experimental IP and EA values for nA (blue plus signs) are from refs (54) and (55). The OT-RSH-computed values for nA and DHTAn are shown as black triangles and red circles, respectively. (b) Optimal range-separation parameter γopt corresponding to the computational results of panel (a).

We now turn to the OT-RSH results for the DHTAn molecules represented as red circles in Figure a. When comparing nA with DHTAn of the same length, we clearly observe a stabilization of, both, the IP and the EA. For instance for n = 5, the IP increases by 0.3 eV and the EA reduces by almost 0.5 eV because of the introduction of the pyrazine and dihydropyrazine groups. We predict that this trend is continued for the longer molecules such that one can expect chemically stable dihydro-tetraaza-acenes up to 7 rings. While our computed EA for DHTAn has the expected trend as a function of n, we observe a peculiar behavior of the IP for which the OT-RSH results predict an almost constant value of 6.6 eV between n = 5 and n = 7 before it considerably reduces for n = 9 and n = 11. We believe that this is related to the underlying structure of the π bands of linear oligoacenes which will be the focus of the next subsection.

Before we close this section, we briefly point out that, as depicted in panel (b) of Figure , the optimal range-separation parameter γopt also follows a clear trend as a function of the molecular length n, where almost the same values for nA and DHTAn have been computed. The longer the molecule, the smaller γopt has to be chosen; in other words, the larger the distance dividing the short from the long range, 1/γopt has to be chosen. Interestingly, we find a linear trend as a function of the inverse length n with a saturation value of γopt = 0.07 bohr for a hypothetical infinitely long acene chain.

π Bands

For pentacene, it has been shown that the orbital energies of the 11 occupied π orbitals can be interpreted as intramolecular band structure grouped into the linking band containing six and the apex band containing five orbitals, respectively. This intramolecular band structure could also be measured by means of angle-resolved photoemission spectroscopy (ARPES).[52] Here, we show that the π orbitals of dihydro-tetraaza-acenes strongly resemble those of the parent oligoacenes and can be equally represented as the intramolecular band structure.

Clearly, to establish an intramolecular energy dispersion relation, orbital energy E versus momentum k, a momentum value k needs to be assigned to a given orbital. Now, for a molecule of finite length, the momentum k is not a good quantum number. For molecules with a quasiperiodicity such as nA and DHTAn, we can nevertheless define a quasimomentum value by computing the Fourier transform of the orbital and assigning the k-value to the maximum intensity in Fourier space. Such theoretical momentum maps are depicted in Figure a for the LUMO + 1, LUMO, and HOMO of 5A (top to bottom) and in Figure b for the LUMO, HOMO, and HOMO – 1 of DHTA5, respectively, where the k-values corresponding to the maximum intensity in momentum space are indicated by the vertical dashed lines. It should be noted that the theoretical momentum maps for 5A have already been experimentally confirmed by analyzing the photoelectron angular distribution within ARPES experiments.[52,53,56−58] Using the technique of photoemission tomography,[59] such measurements are being planned also for DHTAn in future.

Figure 3

(a) Theoretical momentum maps of the LUMO + 1 (π7,0), LUMO (π6,0), and HOMO (π4,1) of 5A with the k-value at maximum intensity indicated as vertical dashed lines. (b) Same as (a) for the LUMO (π7,0), HOMO (π6,0), and HOMO – 1 (π4,1) of DHTA5. (c) Orbital energies of nA and DHTAn molecules plotted vs the respective k-value defined in panel (a). Filled and open symbols correspond to the linking and apex π bands denoted as π and π, respectively. The solid and dashed black lines correspond to the π-band structure of a hypothetical infinite polyacene chain, where empty bands (dashed lines) have been rigidly shifted to match the empty orbitals of the acenes.

As already mentioned earlier, the LUMO and HOMO belong to the linking band with seven and six nodes perpendicular to the long axis, respectively, and no nodal plane parallel to the long molecular axis. Roughly speaking, given the length of DHTA5, L ≈ 12.1 Å, the k-values are simply determined by , where m is the number of nodes. This leads to 1.55 and 1.82 Å–1 for the HOMO (π6,0) and LUMO (π7,0), respectively, while the more precise values determined from the momentum maps, from which also the momentum spread Δk is apparent, are 1.51 and 1.70 Å–1. Similarly, the HOMO – 1 (π4,1) would be assigned a k-value of 1.04 Å–1 from simple node counting arguments, while the momentum map yields 1.18 Å–1.

When plotting the computed orbital energies against the momentum values k determined from the momentum maps, the quasiband structure, or equivalently, the intramolecular band structure can be determined.[52,60] Several interesting observations can be made when comparing the intramolecular band structures of nAs with those from the DHTAn molecules as depicted in Figure c. First, all π orbitals of the oligoacenes are nicely threaded on the energy bands of a hypothetical, infinitely long polyacene chain shown as solid and dashed lines. Note that polyacene actually turns out to be gapless with valence and conduction bands crossing slightly before the Brillouin zone boundary (indicated by a vertical line), but for the sake of comparison, we have shifted the conduction bands of polyacene (dashed lines) such that they match the empty states of the oligoacenes in energy. Second, we observe that all orbital energies of DHTAn are, on average, shifted by ≈1 eV to higher binding energies and are slightly more scattered. We also notice a systematic increase in the k-values from nA to DHTAn as indicated by the brown lines connecting black and red circles for n = 5. We attribute these effects to the presence of the pyrazine and dihydropyrazine groups that are responsible for the dipole moment of the molecule for which we obtain 3.49, 3.84, and 3.21 debye for n = 5, 6, and 7, respectively. As a result of the dipole moment, the electrostatic potential exhibits a step along the long molecular axis, leading to an asymmetric orbital structure and a concomitant deviation from the simple particle-in-a-boxlike band dispersion. The potential step is also responsible for a stronger spatial confinement of the orbital and thereby explains the increased k-values. A third observation has already been mentioned earlier but can be well-confirmed by plotting the band structure, namely, that the symmetry of the frontier orbitals changes when going from nA to DHTAn. While the HOMOs of nA belong to the apex band, π, the HOMOs of DHTAn, for n = 5, 6, 7, are rather part of the linking band π. This results from the presence of two N–H groups that add two more electrons to the π system which fill up one more π orbital. Similarly, the LUMOs of DHTAn have in fact the analogous orbital structure as the LUMO + 1’s of the oligoacenes. Last but not the least, the analysis of the π-band structure of DHTAn allows us also to explain the aforementioned peculiar behavior of the IP (the HOMO orbital energy), which remains at an almost constant value of 6.6 eV for n = 5, 6, 7 but reduces significantly for the 9- and 11-ring molecules. The reason is that in the case of DHTA9 and DHTA11, the topmost orbital of the apex band, π, is predicted to become the HOMO rather than the highest orbital of the linking band as is the case for n = 5, 6, 7.

UV–Vis Absorption Spectra

We have also measured UV–vis absorption spectra of DHTA5,[22] DHTA6, and DHTA7 which are depicted as solid lines in Figure . All three molecules show two main absorption bands, one in the UV near 300 nm and another in the visible range between 450 and 650 nm. Note that the results for DHTA6 and DHTA7 are thus similar to those already observed earlier for the parent system DHTA5.[22] With increasing molecular length, we observe an overall red shift of the absorption features, a phenomenon that is already well-known for oligoacenes and for other aza-acenes.[21] This can be clearly seen for the peak in the UV, but also the onset of the first excitation shifts to lower energies, which is in line with a reduction of the HOMO–LUMO gap.

Figure 4

Solid lines show the experimental UV–vis absorption of compounds DHTAn dissolved in DMSO, whereas the dashed lines show the corresponding calculated TD-DFT spectra.

To analyze the experimental observations, we have simulated the vertical excitations energies of DHTAn within the framework of TD-DFT.[50] To account for the solvent, we have considered its dielectric constant, ε = 2.18, in the expression for the exchange correlation energy given by eq . Compared to the gas-phase calculations with ε = 1, the dielectric screening in the environment of the solvent leads to a reduction of the single-particle HOMO–LUMO gap from the gas-phase values of 5.60, 5.19, and 5.11 eV for DHTA5, 6, and 7, respectively (see Figure ), to 3.80, 3.54, and 3.49 eV in the solvent. Focusing on the lowest energy transition, the excitations energies resulting from the TD-DFT calculation turn out to be slightly smaller than those of the HOMO–LUMO gaps, namely, 2.81, 2.61, and 2.62 eV for DHTA5, 6, and 7, respectively. Within our TD-DFT approach and using an OT-RSH functional, this difference between the HOMO–LUMO gap and the optical excitation energy may indeed be interpreted as exciton binding energies with an accuracy comparable to GW-BSE calculations.[61] Thus, we predict exciton binding energies to range between 0.9 and 1.0 eV, which are typical values to be expected for such molecules in a solvent. Note that a reduction of the exciton binding energy is expected for the molecules in the solid state.[8,9,62,63]

Figure 4 shows the lowest 10 excitations of DHTA5 with their corresponding oscillator strengths and a Gaussian broadening of 0.1 eV, which are represented by a black dashed line. Apart from an overall shift of approximately 0.5 eV to higher energies, the simulated spectrum compares quite favorably to the experimental result. With the help of TD-DFT, we can thereby identify the two major absorption bands separated by about 2 eV to arise primarily from the HOMO–LUMO transition and the HOMO – 1/LUMO + 1 transition, respectively. Owing to the symmetry and nodal structure of the involved orbitals (HOMO–LUMO = π6,0/π7,0 and HOMO – 1/LUMO + 1 = π4,1/π5,1), both excitations exhibit a transition dipole parallel to the long molecular axis with calculated oscillator strengths of 0.62 and 2.03 whose ratio fits excellently to the experimental spectrum. It is important to note that the optical absorption of DHTA5 therefore distinctly differs from pentacene for which the lowest-lying transition is polarized perpendicular to the molecule.

Finally, also for DHTA6 and DHTA7, the calculations predict the lowest energy transition to be polarized parallel to the long molecular axis. Moreover, also the second major absorption bands of DHTA6 and DHTA7 around 4.34 and 4.26 eV, respectively, mainly arise from the HOMO–LUMO + 2 transition polarized parallel to the long axis. Compared to the experiment, a similar shift of about 0.5 eV for both absorption bands is observed. One reason for this discrepancy could be the fact that the effect of the solvent is taken into account only by its dielectric constant, while for a more accurate description of the spectra a more realistic description of the molecule–solvent interactions may be necessary. Also note that the fine structure apparent in the lowest energy absorption band in the experiment is likely due to a vibrational progression which is not taken into account in the theoretical spectra.

Conclusions and Outlook

Building on the already reported synthesis root for DHTA5, in this work, we have succeeded to synthesize and isolate the longer DHTA6 and DHTA7 molecules which we have characterized, both, experimentally and by means of density functional calculations with a nonempirically tuned RSH functional. We have computed the gas-phase electronic structure of linear oligoacene (nA) derivatives in which the benzene rings left and right to the central ring are replaced by rings containing a pyrazine and a dihydropyrazine group, respectively. Despite the fact that such dihydro-tetraaza-acenes (DHTAn, compare Scheme ) contain 4n π electrons and are thus nominally expected to be unstable, our computed IPs and EAs reveal the opposite, namely, an increase in IP and a decrease in EA when going from nA to DHTAn for a given number of rings n. This conclusion supports previous works based on the resonance energies, extra cyclic resonance energy, and nucleus independent chemical shift data of 4n π electron dihydro-aza-pentacenes which are stabilized by aromaticity and conjugations of the dihydropyrazine moieties.[23] We further predict that all π orbitals of DHTAn strongly resemble the one of the parent nA molecules; however, the symmetry of the frontier orbitals changes owing to the addition of two electrons to the π system.

We have also computed photoemission momentum maps for DHTAn which, in future studies, should be tested against experimental observations and could confirm the orbital energies and spatial structure of the frontier orbitals. As a first step toward this direction, we have measured the UV–vis absorption spectra of DHTA5, 6, and 7 in solution and compared them to TD-DFT results. This allows us to identify the origin and polarization of the main optical transitions which are distinctly different from the parent oligoacene series. Moreover, preliminary investigations on the intermolecular packing of DHTAn have emphasized the pronounced role of the molecule’s dipole moment and the importance of hydrogen bonding for these DHTAn molecules. In view of applications in organic electronics, the intramolecular electronic structure together with the expected intermolecular bonding forces makes DHTAn molecules an interesting class of molecules that certainly deserve further investigations regarding thin-film morphologies and bulk electronic structure.