Alkylated Benzodithienoquinolizinium Salts as Possible Non-Fullerene Organic N-Type Semiconductors: An Experimental and Theoretical Study

Three photobicyclized benzodithienoquinolizinium tetrafluoroborates (BPDTQBF4) were prepared and evaluated by UV–Vis and fluorescence spectral, electrochemical analysis, and by theoretical calculations as possible organic n-type semiconductors. Evaluation and comparison of their LUMO levels, HOMO-LUMO energy gaps as monomeric and π-stacked dimers with those of other materials, suggest their potential as organic n-type semiconductors. Calculations of their relative charge carrier mobilities confirmed this potential for one derivative with a long (C-14) alkyl chain appended to the polycyclic planar π-system.


Introduction
The development of new and more efficient organic semiconductors is an active area of research with a plethora of applications. They are normally classified as p-type (electrondonating or hole-conducting), n-type (electron-accepting or electron-conducting), and ambipolar (hole-and electron-conducting) semiconductors. The preparation and the study of electron-donating or p-type semiconductors have received considerable attention [1]. Compared to them, the development of their n-type counterparts has been left behind due in part to their greater atmospheric sensitivity and the rather limited number of molecular scaffolds employed for their synthesis. Research has been mainly concentrated on fullerene derivatives because of their high electron affinity and electron mobility [2][3][4][5][6]. Due to their ease of polymer intercalation and efficient electron transport, fullerene derivatives became an important class of electron acceptors for highly efficient organic solar cells, particularly when coupled with high-performance, electron-donating polymers [7]. Table 1 shows a comparison of some properties of fullerene derivatives and small molecules, when they are used as n-type semiconductors in organic solar cells. Table 1. LUMO energy levels (obtained by electrochemical measurements), electron mobility values, donor/acceptor ratio and power conversion efficiency (PCE%) of fullerene derivatives and some organic small molecules used as n-type semiconductors.

Acceptor.
LUMO Energy (eV) µ e (cm 2 [8] −3.75 1.75 × 10 −4 1:0.7 1.14 PyC61BM [8] −3.72 2.13 × 10 −4 1:0.7 1.95 oo-PDI [9] −3.90 -1:1 8.12 PDI-3 [10] −3.63 -1:2 0.96 PDI-1 [10] −3.67 -1:2 0.13 NDITz [11] −3.99 0.15 --NDINI [12] −4.20 1.75 --C8-NDTI [13] −4.00 0.05 --Despite these advantages, however, recent studies have suggested that degradation of organic solar cells may be due to their high photo-and oxygen-sensitivity, and their tendency to form macroscopic aggregates [14][15][16]. Fullerene and its derivatives are weak absorbers of visible light [17,18]. This feature results in poor light-harvesting properties, thereby substantially reducing the potential of the photocurrent generation of organic solar cells in the UV-Visible range of the solar spectrum. These limitations and the poor synthetic flexibility shown by fullerene have led to efforts to develop other organic n-type semiconductors with a non-fullerene scaffold. Research in this field has concentrated on various oligomers [19,20] and polymeric materials [21,22], with few examples of small molecules. A variety of oligomeric N-hydrogenated/N-methylated pyridinium salts have been used as cathode interfacial layers to increase the electron-mobility without modifying the work function of the electrode [23,24]. Oligothiophenes [25][26][27] have exhibited decreased LUMO energies, together with increased electron mobilities. Many of these properties can be anticipated theoretically, allowing for the identification of potentially interesting materials for the development of n-type organic semiconductors. An example of such computational studies is the theoretical calculation of HOMO-LUMO energy gaps of thiophenebased discotic systems, their charge-transfer rates, and relative charge-mobilities [28][29][30][31]. Theoretical structural and electronic properties of these kinds of compounds have been studied [32][33][34][35]. A series of ten 2-naphthyl and 2-anthrylbithiophene derivatives were synthesized and studied as possible n-type semiconductors [36]. Spectroscopic measurements and X-ray analyses were carried, together with theoretical calculations with DFT methods. Some experimental properties were in good agreement with the theoretical predictions. Theoretical and experimental studies with monomeric and dimeric DAE molecules have been carried out [37]. Theoretical calculations showed that the π extension of the dimer produced a decrease in the LUMO energy by comparison with the monomer. Low values for the internal reorganization energies (λ) were obtained for the dimer when compared with the monomer, leading to a higher semiconductor performance of the former.
Polycyclic aromatic hydrocarbons (PAHs) have been proposed as semiconductors due to their ability to produce a discotic stacking, allowing an efficient charge-transport through their vertical intermolecular axis [38][39][40]. Polycyclic pyridinium salts are positively charged PAHs [41,42] and they constitute a class of electron-deficient planar πsystem whose versatile synthesis can include other electron-donor or acceptor groups. Recently, pyridinium-based molecules have been used in dye-sensitized solar cells [43,44]. A pyridinium fragment, appended to a polyelectrolyte system, enabled efficient electron transport/collection, giving rise to a PCE value of 16.14% in an organic solar cell [45].
The above features of pyridinium-and thiophene-based extended π-systems and their possible formation of a discotic stacking have attracted our attention as a new class of electron-deficient systems with potential use in solar cells. To the best of our knowledge, photocyclized pyridinium salts have never been proposed as potential n-type organic semiconductors. Their synthesis by photocyclization of readily obtained N-arylpyridinium salts has been described for more than thirty years [46,47]. More recently, some of their derivatives such as 2-phenyl-benzo [8,9]quinolizino [4,5,6,7-fed]-phenanthridinylium tetrafluoroborate (PQPBF 4 ) have been shown to form columnar π-stacked aggregates in the solid phase [48]. Although such aggregation in acetonitrile has been disputed [49], spectral evidence for its formation in water has been obtained in our laboratories [50]. Thus, planar electron-deficient π-systems like PQP + or its dithienoquinolizium analog BPDTQ + should be good candidates for promising n-type semiconductors. The incorporation of side-chains of variable size into PQP + derivatives has led to materials with different discotic stackings in the solid state [51,52], suggesting their use in supramolecular assemblies on electronic surfaces [53,54].
In line with other theoretical studies [55] that not only shed light on their properties, but also help identify new interesting compounds by comparing their properties with those of known materials, we investigate in the present communication the properties of polycyclic pyridinium salts as candidates for the development of n-type organic semiconductors. Besides their electron-deficient character and their potential ability to form columnar aggregates, other structural features were included in their synthesis, which were above-mentioned as leading to better n-type organic semiconductors. They include the presence of dithieno substituents to the pyridinium core, and of alkyl sidechains of variable size, all of which should be expected to improve the charge-mobility in these systems. The effect of these structural modifications on a discotic arrangement of these molecules was assessed by calculations of their electronic behavior, injection, and charge-transport properties, HOMO and LUMO energies, and supported by experimental determination of their thermal stability and oxidation/reduction potentials through electrochemical measurements.

Synthesis of Polycyclic Heteroaromatic Salts
Melting points were measured with a capillary Microthermal apparatus and were not corrected. 1 H NMR spectra were obtained with a Bruker Avance 400 MHz and the samples were prepared in acetone-d 6 . The UV-Vis spectra were recorded with a Cary 50 spectrometer. Fluorescence spectra were recorded with a Perkin Elmer LS55 spectrofluorimeter.

Electrochemical Measurements
Cyclic voltammetry experiments were carried out in a PalmSens 4 potentiostat/galvanostat/ impedance analyzer, using a one-compartment cell with a conventional three-electrode arrangement. As the working electrode, a glassy carbon electrode (CH Instruments, TX, USA) was used; as a counter electrode, we used a Pt wire (CH Instruments, TX, USA) and as a reference, we used a Ag/AgCl (3 M KCl) electrode. This last electrode was sealed and separated using a glass tube connected to the solution through a platinum bridge, working as a Lugging capillary. This arrangement avoids any moisture contamination to the working solution at the timescale of the CV experiments [57]. Tetrabutylammonium perchlorate (TBPA) 0.1 M was employed as the supporting electrolyte in 1 mM solutions of the tetrafluoroborate salts in DMF, with scan rates of 100 mV/s.

Thermogravimetric Measurements
The thermogravimetric analysis was performed using a NETZSCH TG 209 F1 Libra thermoanalyzer at a constant heating rate of 3 K/min in a range of 30-500 • C. A constant purging gas flow of 20 mL/min nitrogen was applied with a protective gas flow of 8 mL/min nitrogen. Data analysis was carried out with Proteus Software (version 6.1).

Computational Details
The optimized structures of the monomeric BPDTQ cations and of their tetrafluoroborate salts, and the various conformations of their dimers (parallel and face to face) were calculated with Gaussian 09 [58] employing the hybrid DFT-B3LYP and M062X functionals, with the 6-31G (d) basis set. The electronic spectra of monomeric and dimeric molecules in methanol were obtained by TD-DFT methods with the 6-31G (d) basis set. Solvent effects were mimicked with the polarized continuum model (PCM). Intermolecular non-covalent interactions within the various conformations of the dicationic dimers were depicted with the NCIPLOT software (version 3.0) [59].
Intermolecular orbital overlap for holes and electrons were calculated using the Koopman approximation.

Synthesis of the Polycyclic Heteroaromatic Salts 4a-c
The synthetic route to the heteroaromatic salts 4a-c is shown in Scheme 1. als, with the 6-31G (d) basis set. The electronic spectra of monomeric and dimeric molecules in methanol were obtained by TD-DFT methods with the 6-31G (d) basis set. Solvent effects were mimicked with the polarized continuum model (PCM). Intermolecular noncovalent interactions within the various conformations of the dicationic dimers were depicted with the NCIPLOT software (version 3.0) [59].
Intermolecular orbital overlap for holes and electrons were calculated using the Koopman approximation.

Synthesis of the Polycyclic Heteroaromatic Salts 4a-c
The synthetic route to the heteroaromatic salts 4a-c is shown in Scheme 1.

Spectroscopic Characterization of Tetrasubtituted and Photocyclized Pyridinium Salts
The tetrasubtituted pyridinium salts 3a-c showed two absorption maxima 327 and 312 nm, respectively, in chloroform ( Figure 1). When the absorption spectra of these salts were measured in methanol as the solvent, two absorption maxima were also registered at 378 and 385 nm (7S). The spectra were not affected by the size of the alkyl substituent of the N-phenyl group, but there existed a hypsochromic shift of the maxima passing from the chloroform to methanol medium (−15 and −7 nm), probably because of the major stabilization of the ground state of the molecules with the polar protic solvent. The basic condensation of 2-thienylethanone with benzaldehyde to form propenone (1), followed by acid-catalyzed cyclization of this product with benzaldehyde to give the pyrylium tetrafluoroborate (2), are classical procedures reported previously [56]. Conversion of (2) into the corresponding N-aryl-2,6-bis(2-thienyl)-4-phenyl pyrydinium tetrafluoroborates 3a-c by reaction with 4-alkylanilines also followed a standard procedure, described in the literature [46]. The photocatalytic bicyclization of N-aryl-2,4,6-trisubstituted pyridinium salts has been previously described [47] and the procedure was adapted to the preparation of tetrafluoroborate salts 4a-c from the corresponding pyridinium salts 3a-c.

Spectroscopic Characterization of Tetrasubtituted and Photocyclized Pyridinium Salts
The tetrasubtituted pyridinium salts 3a-c showed two absorption maxima 327 and 312 nm, respectively, in chloroform ( Figure 1). When the absorption spectra of these salts were measured in methanol as the solvent, two absorption maxima were also registered at 378 and 385 nm (7S). The spectra were not affected by the size of the alkyl substituent of the N-phenyl group, but there existed a hypsochromic shift of the maxima passing from the chloroform to methanol medium (−15 and −7 nm), probably because of the major stabilization of the ground state of the molecules with the polar protic solvent.
The optical band gap (Eg) for compounds 4a-c compounds could be estimated by the equation where λ onset is the low-energy absorption edge obtained from the UV-Visible spectra in chloroform ( Figure 2a) employing a geometric method [60]. From λ onset values of 476, 475, and 477 nm for compounds 4a, 4b, and 4c, respectively, E g values of 2.61, 2.60, and 2.60 eV were obtained. Upon cyclization, the extended heterocyclic π−system exhibited several new bands around 360, 430, and 470 nm, in addition to the major absorption around 330 nm, which consisted of at least two resolved bands ( Figure 2). The optical band gap (Eg) for compounds 4a-c compounds could be estimated by the equation where λ(onset) is the low-energy absorption edge obtained from the UV-Visible spectra in chloroform ( Figure 2a) employing a geometric method [60]. From λ(onset) values of 476, 475, and 477 nm for compounds 4a, 4b, and 4c, respectively, Eg values of 2.61, 2.60, and 2.60 eV were obtained. Unlike salts 3, the rigid polycyclic scaffold of 4a-c prevented non-radiative decay processes of their excited species, leading to an emission band around 480 nm ( Figure 3).  The optical band gap (Eg) for compounds 4a-c compounds could be estimated by the equation where λ(onset) is the low-energy absorption edge obtained from the UV-Visible spectra in chloroform ( Figure 2a) employing a geometric method [60]. From λ(onset) values of 476, 475, and 477 nm for compounds 4a, 4b, and 4c, respectively, Eg values of 2.61, 2.60, and 2.60 eV were obtained. Unlike salts 3, the rigid polycyclic scaffold of 4a-c prevented non-radiative decay processes of their excited species, leading to an emission band around 480 nm ( Figure 3). The comparison of the UV-Vis and fluorescence spectra of 4a-c in chloroform and in methanol (Figures 2 and 3) revealed some differences, but are not conclusive as to the possible formation of aggregates in the more hydrophilic methanolic medium. Discotic aggregation in methanol has been suggested for polycyclic systems analogous to 4 on the basis of dynamic light-scattering measurements [51]. Following this view of a postulated face-to-face aggregation in the solid phase of analogous highly planar hydrophobic struc- The comparison of the UV-Vis and fluorescence spectra of 4a-c in chloroform and in methanol (Figures 2 and 3) revealed some differences, but are not conclusive as to the possible formation of aggregates in the more hydrophilic methanolic medium. Discotic aggregation in methanol has been suggested for polycyclic systems analogous to 4 on the basis of dynamic light-scattering measurements [51]. Following this view of a postulated face-to-face aggregation in the solid phase of analogous highly planar hydrophobic structures [48], we cannot discard the formation of salt 4 aggregates in methanol. A more detailed investigation of this possibility, and of its consequences on the absorption spectra of these compounds and on their electron-deficiency was carried out with theoretical tools in Sections 3.4 and 3.5.

Electrochemical Studies
Cyclic voltammetric studies of pyridinium salts 3 and of their cyclized derivatives 4 were carried out to gain information on the potentials of these compounds as electrondeficient materials, and on the effect of the bicyclization of 3 on these potentials.
The nature of the alkyl substituent R in series 3a-c and 4a-c, did not alter the obtained voltammograms. Figure 4 presents the cyclic voltammogram of pyridinium tetrafluoroborate 3a. The same experimental conditions were applied in obtaining the voltametric profile of the bicyclized derivative 4a shown in Figure 5. The comparison of the UV-Vis and fluorescence spectra of 4a-c in chloroform and in methanol (Figures 2 and 3) revealed some differences, but are not conclusive as to the possible formation of aggregates in the more hydrophilic methanolic medium. Discotic aggregation in methanol has been suggested for polycyclic systems analogous to 4 on the basis of dynamic light-scattering measurements [51]. Following this view of a postulated face-to-face aggregation in the solid phase of analogous highly planar hydrophobic structures [48], we cannot discard the formation of salt 4 aggregates in methanol. A more detailed investigation of this possibility, and of its consequences on the absorption spectra of these compounds and on their electron-deficiency was carried out with theoretical tools in Sections 3.4 and 3.5.

Electrochemical Studies
Cyclic voltammetric studies of pyridinium salts 3 and of their cyclized derivatives 4 were carried out to gain information on the potentials of these compounds as electrondeficient materials, and on the effect of the bicyclization of 3 on these potentials.
The nature of the alkyl substituent R in series 3a-c and 4a-c, did not alter the obtained voltammograms. Figure 4 presents the cyclic voltammogram of pyridinium tetrafluoroborate 3a. The same experimental conditions were applied in obtaining the voltametric profile of the bicyclized derivative 4a shown in Figure 5.  The potential of the first reduction process can be related to the electron affinity (EA) of the compounds. For this, it is necessary to relate the electrochemical potential to the vacuum level. Using Equation (1), it is possible to obtain the reduction potential (E red ) relative to the vacuum level (E vac ), the onset reduction potential (E red ), and the Ag/AgCl reference electrode (E Ag/AgCl ) [61,62].
In order to produce a good ohmic contact between an n-type semiconductor and a cathode, and a favorable electron injection into the former, the difference between its LUMO and the cathode work function must be ideally smaller than 0.3 eV [67]. The estimated LUMO energies of the tetrafluoroborate salts 4a-c (≈−3.6 eV), though slightly less negative than those of the other n-type organic semiconductors of Table 3, fulfils this requirement for commonly used electrodes such as Mg, with a work function of −3.66 eV.

Thermogravimetric Measurements
The stability of BPDTQ tetrafluoroborates 4a-c was next evaluated by thermogravimetric analysis. Figure 6 reproduces the TGA of compound 4c, showing the loss in mass of the salt with the increased temperature. As can be seen, the BPDTQ salts were thermally very stable, only starting to decompose and lose mass at temperatures around 350 • C, thus ensuring their use in solar cells.

Theoretical Studies
In the previous sections, the potential of the photocyclized BPDTQ salts as promising organic n-type semiconductors was established by experimental estimates of their HOMO and LUMO energies, which fell within the range of frontier-orbital energies of other organic n-type semiconductors and led to an adequate ohmic contact between them and a metal cathode.
In order to gain a deeper insight into their behavior, we carried out theoretical calculations with three main purposes: (1) to have a clearer picture of the frontier orbitals and of the HOMO → LUMO transition of the isolated cations BPDTQ + ; (2) to investigate the effect of π-stacking of these extended planar structures on their properties as n-type electron-accepting species; and (3) to have an insight into the electron-injection process from the BPTDQ + cation into a metallic surface and to evaluate charge mobility in these species.

Monomeric BPDTQ + Cations
Frontier-orbital calculations of the optimized structures of BPDTQ + cations of 4a-c were carried out with the hybrid DFT-B3LYP and M062X functionals, employing the 6-31G (d) basis set. Transition energies were obtained with TD-DFT calculations. Solvent effects were simulated with the PCM approach. Figure 7 reproduces the HOMO and LUMO of the 4a cation. More detailed results from the frontier molecular orbital analyses

Theoretical Studies
In the previous sections, the potential of the photocyclized BPDTQ salts as promising organic n-type semiconductors was established by experimental estimates of their HOMO and LUMO energies, which fell within the range of frontier-orbital energies of other organic n-type semiconductors and led to an adequate ohmic contact between them and a metal cathode.
In order to gain a deeper insight into their behavior, we carried out theoretical calculations with three main purposes: (1) to have a clearer picture of the frontier orbitals and of the HOMO → LUMO transition of the isolated cations BPDTQ + ; (2) to investigate the effect of π-stacking of these extended planar structures on their properties as n-type electronaccepting species; and (3) to have an insight into the electron-injection process from the BPTDQ + cation into a metallic surface and to evaluate charge mobility in these species.

Monomeric BPDTQ + Cations
Frontier-orbital calculations of the optimized structures of BPDTQ + cations of 4a-c were carried out with the hybrid DFT-B3LYP and M062X functionals, employing the 6-31G (d) basis set. Transition energies were obtained with TD-DFT calculations. Solvent effects were simulated with the PCM approach. Figure 7 reproduces the HOMO and LUMO of the 4a cation. More detailed results from the frontier molecular orbital analyses are given in the Supplementary Materials. The HOMO-LUMO transition involves a charge-transfer process from the benzoand diethieno-rings to the quinolizinium and the 2-phenyl substituent. Table 4 lists some calculated values for the 4a-c cations, which can be compared to the experimentally-derived values of Tables 4 and 5. In general, the nature of the appended alkyl chain did not affect the calculated values, so only the results for the 7tetradecyl BPDTQ + cation are given in the Table.   Table 4. Calculated values of frontier-orbital energies and corresponding absorption of the longestwavelength band of the 7-tetradecyl BPDTQ + cation in methanol, together with its ionization potential IP and electron affinity EA.   The HOMO-LUMO transition involves a charge-transfer process from the benzo-and diethieno-rings to the quinolizinium and the 2-phenyl substituent. Table 4 lists some calculated values for the 4a-c cations, which can be compared to the experimentally-derived values of Tables 4 and 5. In general, the nature of the appended alkyl chain did not affect the calculated values, so only the results for the 7-tetradecyl BPDTQ + cation are given in Table 4. Table 4. Calculated values of frontier-orbital energies and corresponding absorption of the longestwavelength band of the 7-tetradecyl BPDTQ + cation in methanol, together with its ionization potential IP and electron affinity EA. A comparison of these values with those of Tables 4 and 5 revealed a significant departure of theory from the experimentally-derived values. Calculations with the B3LYP/6-31g(d) method led to better results than with the DFT/M062X method, but still only tolerably reproduced the experimental values. LUMO levels in particular, were significantly lower than the experimental value of ≈−3.6 eV. As pointed out at the end of Section 3.3, caution is required when employing a HOMO-LUMO terminology based on different methods of calculation [68]. Calculated values of HOMO and LUMO energy levels may therefore differ from those estimated by optical or electrochemical measurements. By adopting a different method of calculation, with the assumption that the LUMO energy E LUMO = −EA [69], the calculated LUMO energy of cation 4c draws closer to the electrochemically-derived experimental value of −3.6 eV.

π-Stacking in Dimeric BPDTQ + Cations
Although the spectra of Figures 2 and 3 did not offer sufficient evidence for the formation of aggregates in more polar methanol, reports of π-stacking in dimeric PQP + [48,51] and in BPDTQ + cations [50] led us to investigate the possible structures and electronic properties of the dimers of the 4a-c cations.
In the absence of data from crystal packing, we considered two forms of parallel stacking, one in which the two delocalized ring systems were aligned in the same direction (a face-to-face (FF) alignment) and another where they were aligned in opposite directions (face-to-tail (FT) alignment). The two conformations are reproduced in Figure 8 for the hexyl derivative 4a. 31g(d) method led to better results than with the DFT/M062X method, but still only tolerably reproduced the experimental values. LUMO levels in particular, were significantly lower than the experimental value of  -3.6 eV. As pointed out at the end of Section 3.3, caution is required when employing a HOMO-LUMO terminology based on different methods of calculation [68]. Calculated values of HOMO and LUMO energy levels may therefore differ from those estimated by optical or electrochemical measurements. By adopting a different method of calculation, with the assumption that the LUMO energy ELUMO = −EA [69], the calculated LUMO energy of cation 4c draws closer to the electrochemically-derived experimental value of −3.6 eV.

π-Stacking in Dimeric BPDTQ + Cations
Although the spectra of Figures 2 and 3 did not offer sufficient evidence for the formation of aggregates in more polar methanol, reports of π-stacking in dimeric PQP + [48,51] and in BPDTQ + cations [50] led us to investigate the possible structures and electronic properties of the dimers of the 4a-c cations.
In the absence of data from crystal packing, we considered two forms of parallel stacking, one in which the two delocalized ring systems were aligned in the same direction (a face-to-face (FF) alignment) and another where they were aligned in opposite directions (face-to-tail (FT) alignment). The two conformations are reproduced in Figure 8 for the hexyl derivative 4a.

Face-to-face (FF) alignment
Face-to-tail (FT) alignment  Table 5 lists some calculated properties of these dimers: the distance between their π-stacked planes, their dipole moments, the interaction energies between the two monomers, and their longest-wavelength λmax values.
A non-covalent interaction (NCI) analysis was performed for all dimers of Table 5, and the results are summarized in Figure 9. More detailed results from the NCI analyses are given in the Supplementary Materials.  Table 5 lists some calculated properties of these dimers: the distance between their πstacked planes, their dipole moments, the interaction energies between the two monomers, and their longest-wavelength λ max values.
A non-covalent interaction (NCI) analysis was performed for all dimers of Table 5, and the results are summarized in Figure 9. More detailed results from the NCI analyses are given in the Supplementary Materials.
As can be seen, interactions are governed by π − π charge repulsions between the two cationic monomers, shown in red, which are compensated by a large, green iso-surface of attractive van-der-Waals interactions between the two polycyclic systems [70].
Interestingly, face-to-face dimers present an additional non-covalent interaction between the parallel alkyl chains, giving rise to an additional attraction between the two monomers. Such quasi-covalent interactions between parallel chains have been described previously [71]. This "fastener effect" between alkyl chains in solid samples should be responsible for an increased conductivity and electron mobility in their resulting films. More detailed NCI analyses are shown the in Supplementary Materials (Figure S10).
The dipole moments listed in Table 5 point to significant differences between the FF and the FT conformations, giving rise to correspondingly distinct interaction energies. Relatively non-polar FT dimers should be favored in non-polar media, whereas an FF alignment should predominate in polar solvents, giving rise to J-aggregates. In the latter media, the non-covalent attraction between the alkyl chains should provide an additional stabilization to these conformers. Table 5 lists some calculated properties of these dimers: the distance between their π-stacked planes, their dipole moments, the interaction energies between the two monomers, and their longest-wavelength λmax values.
A non-covalent interaction (NCI) analysis was performed for all dimers of Table 5, and the results are summarized in Figure 9. More detailed results from the NCI analyses are given in the Supplementary Materials. Figure 9. Face-to-face (FF) and face-to-tail (FT) π-stacked dimers of cations from 4a-c, with noncovalent interactions obtained from NCI analyses.
A noteworthy result from Table 5 is the bathochromic shift of the longest-wavelength λ max value of the dimers when compared with the corresponding values in Table 4 for the monomeric cations. Though small (<10 nm), using the B3LYP functional, this red shift points to a HOMO-LUMO gap for the discotic π-stacked aggregates of compounds 4a-c smaller than the ones listed in Table 4 for the monomeric cations. Following previous suggestions of the formation of columnar aggregates in the solid phase [48] for analogous systems, leading to materials with different discotic stackings [51,52], our calculations point to a similar possibility for polycyclic systems 4a-c, suggesting their use as n-type semiconductors in supramolecular assemblies.

Charge Mobility in Dimeric BPDTQ + Cations
Charge mobility is an important factor to be considered in organic semiconductors. Transport properties are commonly discussed with the aid of Marcus theory [55], though other approaches are also found in the literature [72][73][74][75][76][77][78].
Having assumed the two parallel alignments for our cationic dimers, we estimated the charge transport in these discotic systems employing the Marcus theory [55,79,80] and Einstein's equation for a one-dimensional diffusion process [81,82].
According to the latter, the diffusion µ of a charge or a hole between two parallel surfaces separated by a distance l is given by Equation (4) where e is the electron charge; k B is the Boltzmann's constant; and k ET is the charge-transfer rate constant between the two contiguous molecules. The k ET value can be obtained from Equation (5), derived from Marcus theory where t is the electronic coupling between the two contiguous molecules and λ is the reorganization energy of the system. The first term t depends on the distance between the parallel monomers and can be calculated using different methodologies [83]. In the present work, we employed Koopman's methodology using the HOMO and LUMO energies of the interacting monomers [84][85][86], and the hopping process where a dimeric system of two Scheme 2. Hole/electron exchange between two molecules M in a hopping process.
The reorganization energy λ of Equation (5) comprises an inner and an outer contribution. The first arises from rapid changes in the molecular geometries of the donor (D) and acceptor (A) species, and the latter from nuclear polarization/relaxation of the surrounding medium upon an electron transfer. In our case, this outer contribution can be neglected because our system is a static crystal pack [87]. Thus, the reorganization energy λ (inner) can be obtained from Equations (6), (7), and (10) for holes, and Equations (8) λ = λ  λ D (10) Table 6 lists the obtained values for parameters t and λ, μ, and kET, for the mobility of electrons and holes in all cationic dimers from the BPDTQ + salts. Table 6. Calculated values of the electronic coupling t, the reorganization energy λ, the chargetransfer rate constant kET, and the relative mobility μ for all cationic dimers from 4a-c. An inspection of Table 6 revealed differences among the three dimers of 4a-c and between each pair of FF/FT conformations. These variations point to possible differences in the behavior of these π-stacked dimers as semiconductors.

Dimer
With the exception of 4c dimers, the relative hole (μh) and electron (μe) mobilities of all dimeric species are of the same order of magnitude, with μh < μe for the FF conformations, and μh >> μe for the FT conformations. This amounts to predicting an ambipolar, or even a p-character for semiconductors obtained from these dimeric molecules. Dimers 4c present a distinct behavior. For both conformers FF and FT, the relative hole mobility μh is much smaller than the electron mobility μe, anticipating for these dimeric species an n-type semiconductor behavior. Such low hole mobilities are a result of internal hole reorganization energies λh almost four times greater than those of electrons λe. The fact that Scheme 2. Hole/electron exchange between two molecules M in a hopping process.
The reorganization energy λ of Equation (5) comprises an inner and an outer contribution. The first arises from rapid changes in the molecular geometries of the donor (D) and acceptor (A) species, and the latter from nuclear polarization/relaxation of the surrounding medium upon an electron transfer. In our case, this outer contribution can be neglected because our system is a static crystal pack [87]. Thus, the reorganization energy λ (inner) can be obtained from Equations (6), (7), and (10) for holes, and Equations (8)- (10) for electrons, where E(M) is the total energy of species M and the superscripts A and D refer to the electron-accepting and -donor molecules. Table 6 lists the obtained values for parameters t and λ, µ, and k ET , for the mobility of electrons and holes in all cationic dimers from the BPDTQ + salts. Table 6. Calculated values of the electronic coupling t, the reorganization energy λ, the charge-transfer rate constant k ET , and the relative mobility µ for all cationic dimers from 4a-c. An inspection of Table 6 revealed differences among the three dimers of 4a-c and between each pair of FF/FT conformations. These variations point to possible differences in the behavior of these π-stacked dimers as semiconductors.
With the exception of 4c dimers, the relative hole (µ h ) and electron (µ e ) mobilities of all dimeric species are of the same order of magnitude, with µ h < µ e for the FF conformations, and µ h >> µ e for the FT conformations. This amounts to predicting an ambipolar, or even a p-character for semiconductors obtained from these dimeric molecules. Dimers 4c present a distinct behavior. For both conformers FF and FT, the relative hole mobility µ h is much smaller than the electron mobility µ e , anticipating for these dimeric species an n-type semiconductor behavior. Such low hole mobilities are a result of internal hole reorganization energies λ h almost four times greater than those of electrons λ e . The fact that dimers 4c, with the longest C 14 chain, show a strong n-type character, at variance with analogs 4a and 4b suggests that the size of the 7-substituent alkyl chain plays an important role in the charge-transport and mobility of these dimers. This echoes previous reports of changes in the hopping mobility of discotic alkylated hexabenzocoronenes with the increased size of the alkyl substituents [88]. Thus, trapping and recombination in these discotic materials are retarded by the increased size of the alkyl chains, a mechanism that may also operate for materials derived from 4c.

Conclusions
Photobicyclization of N-aryl-2,6-dithienyl-4-phenyl pyridinium salts led to the formation of benzodithienoquinolizinium tetrafluoroborates (BPDTQBF 4 ) that were investigated by spectral, electrochemical, and theoretical tools as possible n-type semiconductors for solar cells.
Their UV-Vis spectra and their voltamograms in DMF allowed for the estimation of their LUMO levels and HOMO-LUMO energy gaps, which had values compatible with those of other n-type semiconductors. Their UV-Vis spectra in chloroform and in methanol as well as their fluorescence in the same solvents did not show any conclusive evidence of aggregation, as might be surmised from reports of analogous systems in the literature. Nevertheless, theoretical calculations of the monomeric cations 4a-c and of their π-stacked dimers in two different conformations (face-to-face, FF, and face-to-tail, FT) indicated that aggregation could lead to smaller HOMO-LUMO gaps, with a beneficial impact on the behavior of these materials as n-type semiconductors. In order to substantiate this expectation and obtain a deeper view on the many aspects that govern this complex behavior, preliminary calculations on their relative charge mobilities were performed. Estimates of relative charge mobilities for the two conformations of the cationic dimers from 4a-c, based on Einstein's diffusion equation and Marcus' theory, led to inconclusive results, with the characterization of dimers from cation 4c, with a long (C 14 ) alkyl chain appended to the polycyclic core as a probable, potentially interesting n-type semiconductor.
In conclusion, in the present communication, the preparation of a new family of alkylated benzodithienoquinolizinium salts was described. Their electronic properties, estimated from spectroscopic and electrochemical measurements in solution, coupled with theoretical calculations in the gas-phase of their monomers and dimeric species, suggest their potential use as organic n-type semiconductors. In order to confirm their potentiality, measurements in the solid phase, and with thin films are necessary, aside from calculations based on their crystal structures. We will address these necessary issues in a future work, with the most promising compound of the series.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/ma14216239/s1, Figure S1: NMR spectrum of compound 3a in acetone-d 6 . Figure S2: NMR spectrum of compound 3b in acetone-d 6 . Figure S3: NMR spectrum of compound 3c in acetone-d 6 . Figure S4: NMR spectrum of compound 4a in acetone-d 6 . Figure S5: NMR spectrum of compound 4b in acetone-d 6 . Figure S6: NMR spectrum of compound 4c in acetone-d 6 . Figure S7: Absorption spectra of pyridinium salts 3a, 3b, and 3c in methanol (ca. 9 × 10 −4 M). Figure S8: Cyclic voltammograms obtained for tetrasubtituted pyridinium salts 3b and 3c in DMF (c = 1mM), with TBPA 0.1 M as the supporting electrode and a scan rate of 100 mV/s. Figure S9: Cyclic voltammograms obtained for BPDTQ salts 4b and 4c in DMF (c = 1 mM) with TBPA 0.1 M as the supporting electrode and a scan rate of 100 mV/s. Figure S10: NCI 3D representation of isosurfaces and analysis of the s((λ 2 )ρ) graph for the face-to-face (FF) and face-to-tail (FT) dimers. The color-scale code for the 3D representations is shown. Figure S11: Surfaces of HOMO (a) and LUMO (b) of the FT dimer of 4a + . Figure S12: Surfaces of HOMO (a) and LUMO (b) of the FF dimer of 4a + . Figure S13: Surfaces of HOMO (a) and LUMO (b) of 4b + . Figure S14: Surfaces of HOMO (a) and LUMO (b) of the FT dimer of 4b + . Figure S15: Surfaces of HOMO (a) and LUMO (b) of the FF dimer of 4b + . Figure S16