Synthesis, Solution, and Solid State Properties of Homological Dialkylated Naphthalene Diimides—A Systematic Review of Molecules for Next-Generation Organic Electronics

This systematic study aimed at finding a correlation between molecular structure, solubility, self-assembly, and electronic properties of a homological series of N-alkylated naphthalene diimides (NDIs). NDIs are known for their n-type carrier mobility and, therefore, have potential in the field of organic electronics, photovoltaics, and sensors. For the purpose of this study, nine symmetrical N,N′-dialkylated naphthalene diimides (NDIC3-NDIC11) were synthesized in the reaction of 1,4,5,8-naphthalenetetracarboxylic dianhydride with alkylamines ranging from propyl- to undecyl-. The NDIs were characterized by spectroscopic (NMR, UV-Vis, FTIR), microscopic, and thermal methods (TGA and DSC), and X-ray diffraction (XRD). Our experimental study, extensively referring to findings reported in the literature, indicated that the NDIs revealed specific trends in spectroscopic and thermal properties as well as solubility and crystal morphology. The solubility in good solvents (chloroform, toluene, dichlorobenzene) was found to be the highest for the NDIs substituted with the medium-length alkyl chains (NDIC5–NDIC8). Systematic FTIR and XRD studies unraveled a distinct parity effect related to the packing of NDI molecules with odd or even numbers of methylene groups in the alkyl substituents. The NDIs with an even number of methylene groups in the alkyl substituents revealed low-symmetry (P1−) triclinic packing, whereas those with an odd number of carbon atoms were generally monoclinic with P21/c symmetry. The odd–even parity effect also manifested itself in the overlapping of the NDIs’ aromatic cores and, hence, the π-π stacking distance (dπ-π). The odd-numbered NDIs generally revealed slightly smaller dπ-π values then the even-numbered ones. Testing the NDIs using standardized field-effect transistors and unified procedures revealed that the n-type mobility in NDIC6, NDIC7, and NDIC8 was 10- to 30-fold higher than for the NDIs with shorter or longer alkyl substituents. Our experimental results indicate that N,N′-alkylated NDIs reveal an optimum range of alkyl chain length in terms of solution processability and charge transport properties.


Introduction
The N,N -dialkylated naphthalene diimides (NDIs) represent the group of aromatic diimides also referred to as rylene dyes (Figure 1). The NDIs and their larger-core counterparts, perylenediimides (PDIs, Figure 1), are the most extensively studied rylene dyes. NDIs are molecules with a remarkable planar structure built up of fused benzene rings ( Figure 1). They exhibit high chemical, thermal, and oxidative stability. The high stability in conjunction with their redox activity [1], spectroscopic [2,3], and electronic properties [4,5] make NDIs interesting for materials scientists and supramolecular chemists. NDIs reveal Historically, NDIs were first mentioned as early as in the 1930s in the patents and articles on their use as various dyes [10][11][12][13][14]. The use of NDIs as dyes continued until the mid-1970s. At the turn of the 1970s and 1980s, the NDIs attracted interest from biochemists due to their ability to achieve intercalative binding of DNA [15]. Later in the 1980s, NDIs gained the interest of materials scientists due to their potential in organic field-effect transistors (OFETs), electroluminescent diodes (OLEDs), and other, also flexible, electronic devices, energy storage materials, and sensorics. NDI-based OFETs typically reveal electron mobility oscillating around 0.2 cm 2 × V −1 × s −1 . The OFET based on the N,N′-dicyclohexyl-2,6-naphthalene NDI revealed electron mobility reaching as high as 6 cm 2 × V −1 × s −1 [7,[16][17][18]. NDIs have been successfully used to enhance the electronic performance of polymerbased OFETs [19][20][21]. Lithium-ion batteries with NDI-based electrodes revealed discharge potentials ranging from 2.3 to 2.9 V [22]. NDIs were successfully used as acceptors in organic light-emitting diodes (OLEDs) [23] or in flexible photovoltaic cells (OPV) with power conversion efficiency of 7.65% [24]. NDIs were also used as molecular sensors, playing the role of acceptors in host-guest-type complexes or as an aromatic electron acceptor in synthetic ion channels [25,26].
Interest in the synthesis and applications of NDIs is still on the rise despite decades of preceding research [27][28][29][30][31]. There are novel applications reported in the fields of organic electronics, modern power sources, or power storage units [29,30,32]. NDIs are still studied for their photophysical properties [3]. Nonetheless, in the literature, consistent systematic investigations of homological series of NDIs are still relatively rare [33]. Such a survey for symmetrical NDIs is necessary in order to support the engineering of NDIbased materials. The research papers typically focus on selected properties of arbitrary selected NDIs. In our previous paper, we have already faced the problem of selecting NDI derivatives for specific purposes in the fabrication of OFETs [21]. Motivated by the issues we have experienced in our previous research, herein, we extensively survey-both in Historically, NDIs were first mentioned as early as in the 1930s in the patents and articles on their use as various dyes [10][11][12][13][14]. The use of NDIs as dyes continued until the mid-1970s. At the turn of the 1970s and 1980s, the NDIs attracted interest from biochemists due to their ability to achieve intercalative binding of DNA [15]. Later in the 1980s, NDIs gained the interest of materials scientists due to their potential in organic field-effect transistors (OFETs), electroluminescent diodes (OLEDs), and other, also flexible, electronic devices, energy storage materials, and sensorics. NDI-based OFETs typically reveal electron mobility oscillating around 0.2 cm 2 × V −1 × s −1 . The OFET based on the N,N -dicyclohexyl-2,6-naphthalene NDI revealed electron mobility reaching as high as 6 cm 2 × V −1 × s −1 [7,[16][17][18]. NDIs have been successfully used to enhance the electronic performance of polymer-based OFETs [19][20][21]. Lithium-ion batteries with NDI-based electrodes revealed discharge potentials ranging from 2.3 to 2.9 V [22]. NDIs were successfully used as acceptors in organic light-emitting diodes (OLEDs) [23] or in flexible photovoltaic cells (OPV) with power conversion efficiency of 7.65% [24]. NDIs were also used as molecular sensors, playing the role of acceptors in host-guest-type complexes or as an aromatic electron acceptor in synthetic ion channels [25,26].
Interest in the synthesis and applications of NDIs is still on the rise despite decades of preceding research [27][28][29][30][31]. There are novel applications reported in the fields of organic electronics, modern power sources, or power storage units [29,30,32]. NDIs are still studied for their photophysical properties [3]. Nonetheless, in the literature, consistent systematic investigations of homological series of NDIs are still relatively rare [33]. Such a survey for symmetrical NDIs is necessary in order to support the engineering of NDI-based materials. The research papers typically focus on selected properties of arbitrary selected NDIs. In our previous paper, we have already faced the problem of selecting NDI derivatives for specific purposes in the fabrication of OFETs [21]. Motivated by the issues we have experienced in our previous research, herein, we extensively survey-both in terms of literature review and experimentally-the physical properties of homological series of symmetrically N,Ndialkylated NDIs. The NDIs in this study were synthesized, purified, and analyzed in a unified, consistent way which, therefore, enabled their direct comparison. In addition to the experimental data, for the sake of comparison, here, we also provide and discuss experimental results previously reported in the literature. This extensive study unveils the characteristic trends in the physical properties of N,N -dialkylated NDIs. We demonstrate here that NDIs substituted with moderately long alkyl chains reveal the highest solubility in good solvents and the highest charge carrier mobility. Calorimetric analysis of the full array of n-propyl-through n-undecyl-substituted NDIs indicated that a stepwise extending of the alkyl substituents leads to a gradual fading of aromatic-related and increase in alkylrelated features with a relatively well-distinguished range of molecules with intermediate properties. Based on such observations, we demonstrate and discuss the patterns of properties characteristic of a homological series of N,N -dialkylated NDIs.

Remarks on Synthesis: Towards Increasing the Synthesis Scale
The general synthetic route of NDIs is shown in Figure 2.
terms of literature review and experimentally-the physical properties of homological series of symmetrically N,N′-dialkylated NDIs. The NDIs in this study were synthesized, purified, and analyzed in a unified, consistent way which, therefore, enabled their direct comparison. In addition to the experimental data, for the sake of comparison, here, we also provide and discuss experimental results previously reported in the literature. This extensive study unveils the characteristic trends in the physical properties of N,N′-dialkylated NDIs. We demonstrate here that NDIs substituted with moderately long alkyl chains reveal the highest solubility in good solvents and the highest charge carrier mobility. Calorimetric analysis of the full array of n-propyl-through n-undecyl-substituted NDIs indicated that a stepwise extending of the alkyl substituents leads to a gradual fading of aromatic-related and increase in alkyl-related features with a relatively well-distinguished range of molecules with intermediate properties. Based on such observations, we demonstrate and discuss the patterns of properties characteristic of a homological series of N,N′-dialkylated NDIs.

Remarks on Synthesis: Towards Increasing the Synthesis Scale
The general synthetic route of NDIs is shown in Figure 2. The synthetic protocols of NDIs reported in the literature involve heating of 1,4,5,8naphthalenetetracarboxylic dianhydride with an amine in a solvent at elevated temperature. Depending on the temperature (ranging typically from 100 to 220 °C) and reactivity of the amine, the reaction lasted from a few to 48 h [34][35][36][37]. After the reaction, the crude product was usually precipitated with methanol or ethanol, and separated by filtration. The separated crude product was purified by column chromatography, sublimation, or recrystallization.
Most of the synthetic protocols are typically based on the 1-4 mmol scale [34,35,[38][39][40]. Due to extensive analytical needs in this study, larger quantities were required. That motivated a search for an alternative synthetic route permitting at least 4 g (7-11 mmol) scale. Experimental verification of the reported synthetic protocols indicated that the most scale-limiting factors were the viscosity of the reaction medium and the efficiency of final NDI isolation. Typically, DMF [34,35], quinoline [41,42], or pyridine [22,36] were used as reaction media. The imidazole-toluene mixture was found a less reported reaction medium [37,43]. The use of molten imidazole rather than quinoline as a reaction medium was first reported by Langhals in the synthesis of PDIs [43]. Both the imidazole and quinoline were found to increase the yield, especially in the case of less reactive, sterically hindered amines. Both imidazole and quinoline can increase the solubility of the anhydride ( Figure  2). The use of imidazole permits decreasing the reaction temperature from approx. 180-230 °C (typical range for quinoline) down to 100-130 °C [43], which points towards a catalytic influence of imidazole on the reaction. Despite an unambiguously positive impact on the reaction, the role of imidazole was not finally clarified. Toluene lowers viscosity and, therefore, facilitates stirring and prevents splashing of the reaction mixture [8,[44][45][46]. The synthetic protocols of NDIs reported in the literature involve heating of 1,4,5,8naphthalenetetracarboxylic dianhydride with an amine in a solvent at elevated temperature. Depending on the temperature (ranging typically from 100 to 220 • C) and reactivity of the amine, the reaction lasted from a few to 48 h [34][35][36][37]. After the reaction, the crude product was usually precipitated with methanol or ethanol, and separated by filtration. The separated crude product was purified by column chromatography, sublimation, or recrystallization.
Most of the synthetic protocols are typically based on the 1-4 mmol scale [34,35,[38][39][40]. Due to extensive analytical needs in this study, larger quantities were required. That motivated a search for an alternative synthetic route permitting at least 4 g (7-11 mmol) scale. Experimental verification of the reported synthetic protocols indicated that the most scale-limiting factors were the viscosity of the reaction medium and the efficiency of final NDI isolation. Typically, DMF [34,35], quinoline [41,42], or pyridine [22,36] were used as reaction media. The imidazole-toluene mixture was found a less reported reaction medium [37,43]. The use of molten imidazole rather than quinoline as a reaction medium was first reported by Langhals in the synthesis of PDIs [43]. Both the imidazole and quinoline were found to increase the yield, especially in the case of less reactive, sterically hindered amines. Both imidazole and quinoline can increase the solubility of the anhydride (Figure 2). The use of imidazole permits decreasing the reaction temperature from approx. 180-230 • C (typical range for quinoline) down to 100-130 • C [43], which points towards a catalytic influence of imidazole on the reaction. Despite an unambiguously positive impact on the reaction, the role of imidazole was not finally clarified. Toluene lowers viscosity and, therefore, facilitates stirring and prevents splashing of the reaction mixture [8,[44][45][46]. Moreover, the imidazole-toluene mixture facilitates NDI isolation: toluene-unlike the aforementioned DMF, quinoline, or pyridine-can be easily evaporated. Imidazole can be effectively removed during washing with methanol. The yields for 5 mmol-scale synthesis in the imidazole-toluene mixture were reported to reach 80%. These advantages underpinned the decision to choose the imidazole-toluene mixture in our study. Because of the reaction scale, conventional multiple washing on a filter funnel was replaced with continuous washing with methanol in the Soxhlet apparatus. Additional purification of NDIC3 and NDIC4 (short-column chromatography/sublimation) was necessary due to contamination of NDIs with monoimides resulting probably from the low boiling point of n-propylamine and n-butylamine (see Table 1); evaporation of the amines probably caused a variation in the reagents' ratio over time. Final NDIs were obtained with 35-88% yields ( Table 1). The highest yields, exceeding 80%, were reached for NDIC5, NDIC6, and NDIC9. The yields were found, to some extent, correlated with the amines' boiling points. In principle, the amines with boiling points in the range from 100 to 220 • C permitted the higher yields (Table 1).

UV-Vis Spectroscopy: Solubility of the NDIs
The UV-Vis spectra of NDIs in 0.1 mM chloroform solutions, their description, and molar absorptivity (ε) data are presented in Figure S1 and Table S3 in the Supplementary Materials. Generally, all NDIs showed the same characteristic UV-Vis absorption features, regardless of the length of the alkyl chains. In the UV-Vis spectra of the NDIs, three characteristic sharp absorption bands at 342, 360, and 380 nm were identified. Such bands are characteristic of symmetrical NDIs and are related to the π-π* z-polarized transition [47,48]. The position of each band is practically constant irrespective of the alkyl chains' length. Such absorption characteristics provide the perfect opportunity to quantify the solubility of NDIs in different solvents based on an analysis of the UV-Vis spectra.
It is generally known that NDIs are soluble in the lipophilic solvents like chloroform, dichloromethane, or toluene, or in polar aprotic solvents like DMSO or DMF [49][50][51][52]. The solubility of NDIs in cyclohexane, n-hexane, or n-heptane is significantly lower. In our study, we tested and verified the solubility of NDIs in solvents commonly used in the fabrication of organic electronic films [5,51,53]. The solubility of NDIs was, therefore, tested in chloroform, 1,2-dichlorobenzene, and toluene-good solvents for NDIs, as well as n-heptane-a poor solvent [49,53,54]. Our solubility tests indicated that the solubility of NDIC5 through NDIC8 and of NDIC10 in chloroform exceeded 100 mg/mL. The solubility of the NDIC3, NDIC4, NDIC9, and NDIC11 in chloroform was 46.9 ± 1.5, 65.5 ± 1.6, 47 ± 1.0, and 15.15 ± 0.6 mg/mL, respectively. Considering applications of NDIs in, for instance, printed organic electronics, that solubility in the case of all NDIs exceeds usual needs since typically used concentrations oscillate within 2-20 mg/mL [5,51,53]. Conversely, the solubility of NDIs in n-heptane ( Figure 3 and Table S4) was two orders of magnitude lower than in chloroform. Generally, solubility in n-heptane was found to increase with increasing alkyl chain length. The solubility of NDIC11 was found approx. 6-fold higher than that of NDIC3. Plotting the solubility against alkyl chain length ( Figure 3a) seems to reveal two distinct trends: a sharp power-like increase from NDIC3 through NDIC6, and a quasi-linear trend starting from a local minimum observed for NDIC7 up to NDIC11. fold higher than that of NDIC3. Plotting the solubility against alkyl chain length ( Figure  3a) seems to reveal two distinct trends: a sharp power-like increase from NDIC3 through NDIC6, and a quasi-linear trend starting from a local minimum observed for NDIC7 up to NDIC11. The solubility of NDIs in toluene and 1,2-dichlorobenzene revealed different trends. In both these cases, maximum solubility was observed for NDIC6 (Figure 3b,c and Tables S5 and S6). NDIC5 through NDIC8 reveal the best solubility in toluene and 1,2-dichlorobenzene, hence distinguishing them in the whole range of NDIs. Extending the alkyl chains beyond 8 carbon atoms significantly reduced the solubility. Despite the NDIC3 has the shortest alkyl substituents among the tested NDIs it revealed relatively good solubility in 1,2-dichlorobenzene. This is probably due to its unique crystal packing that is discussed later in the article.

Characteristic Features of NMR Spectra of the Symmetrically N-Substituted NDIs
Analysis of the 13 C NMR spectra is intuitive and the spectral features can be unambiguously assigned to NDIs' molecular structures. Briefly: the 13 C NMR spectra reveal characteristic signal of the imide carbonyl group at 162 ppm, three signals of the aromatic carbons at 126-131 ppm, and a set of peaks characteristic for the alkyl chains on the NDIs' molecules. The collection of 13 C NMR data is provided in Table S1 in the Supplementary Materials. All NMR spectra are provided as a supplementary file to this article.
The 1 H NMR spectra (Table S2 in (Table S3). Most of the 1 H NMR spectral features are in agreement with the literature [21,34,38,39,[55][56][57][58][59][60][61]. However, the multiplicity of the signal at 4.2 ppm originating from the protons of methylene group adjacent to the imide nitrogen atoms ( Figure 4) often seems misinterpreted in the literature [38,55,56,59,61]. Linear alkyl chains and the symmetry of the NDIs' molecules intuitively suggest that coupling of the protons of the methylene group adjacent to the nitrogen atom with the two protons of the subsequent methylene group should result in a simple triplet. Hence, in the literature, the signals appearing within 4.19-4.24 ppm are typically considered triplets [38,55,59,61]. Detailed analysis of NDIC3 and NDIC5 spectra revealed, however, that the middle peak in the discussed range is flanked by two additional peaks ( Figure 4a). In other NDIs, the flanking peaks partially overlap with the middle peak, which appears as middle peak broadening ( Figure 4b). Thus, the signals appearing from 4.19 to 4.24 ppm should be considered a second-order spin system; that is, and, therefore, they are multiplets. The multiplet character of the signal at 4.19-4.24 ppm is similar to the AA'XX' spin system characteristic of, e.g., aryl n-propyl sulfones [62]. On the other hand, in the case of a multiplet with broadened middle line description, assigning them as triplets is conditionally acceptable. In the case when all five lines in a multiplet are visible, referring to them as a "multiplet" or "AA'XX' spin system" is definitely more appropriate. Moreover, magnetic The solubility of NDIs in toluene and 1,2-dichlorobenzene revealed different trends. In both these cases, maximum solubility was observed for NDIC6 ( Figure 3b,c and Tables S5 and S6). NDIC5 through NDIC8 reveal the best solubility in toluene and 1,2dichlorobenzene, hence distinguishing them in the whole range of NDIs. Extending the alkyl chains beyond 8 carbon atoms significantly reduced the solubility. Despite the NDIC3 has the shortest alkyl substituents among the tested NDIs it revealed relatively good solubility in 1,2-dichlorobenzene. This is probably due to its unique crystal packing that is discussed later in the article.

Characteristic Features of NMR Spectra of the Symmetrically N-Substituted NDIs
Analysis of the 13 C NMR spectra is intuitive and the spectral features can be unambiguously assigned to NDIs' molecular structures. Briefly: the 13 C NMR spectra reveal characteristic signal of the imide carbonyl group at 162 ppm, three signals of the aromatic carbons at 126-131 ppm, and a set of peaks characteristic for the alkyl chains on the NDIs' molecules. The collection of 13 C NMR data is provided in Table S1 in the Supplementary Materials. All NMR spectra are provided as a supplementary file to this article.
The 1 H NMR spectra (Table S2 in the Supplementary Materials) include a characteristic singlet from the four aromatic naphthyl protons at 8.74-8.78 ppm and a set of signals corresponding to n-alkyl chains. The signals of the alkyl chains appear in the range of 0.87-1.81 ppm and at 4.19-4.21 ppm (Table S3). Most of the 1 H NMR spectral features are in agreement with the literature [21,34,38,39,[55][56][57][58][59][60][61]. However, the multiplicity of the signal at 4.2 ppm originating from the protons of methylene group adjacent to the imide nitrogen atoms ( Figure 4) often seems misinterpreted in the literature [38,55,56,59,61]. Linear alkyl chains and the symmetry of the NDIs' molecules intuitively suggest that coupling of the protons of the methylene group adjacent to the nitrogen atom with the two protons of the subsequent methylene group should result in a simple triplet. Hence, in the literature, the signals appearing within 4.19-4.24 ppm are typically considered triplets [38,55,59,61]. Detailed analysis of NDIC3 and NDIC5 spectra revealed, however, that the middle peak in the discussed range is flanked by two additional peaks ( Figure 4a). In other NDIs, the flanking peaks partially overlap with the middle peak, which appears as middle peak broadening ( Figure 4b). Thus, the signals appearing from 4.19 to 4.24 ppm should be considered a second-order spin system; that is, and, therefore, they are multiplets. The multiplet character of the signal at 4.19-4.24 ppm is similar to the AA'XX' spin system characteristic of, e.g., aryl n-propyl sulfones [62]. On the other hand, in the case of a multiplet with broadened middle line description, assigning them as triplets is conditionally acceptable. In the case when all five lines in a multiplet are visible, referring to them as a "multiplet" or "AA'XX' spin system" is definitely more appropriate. Moreover, magnetic non-equivalence of the two methylene protons neighboring the nitrogen atom causes increased spectral complexity of the adjacent, following the methylene group (the second from the N-atom, see Figure 4c). non-equivalence of the two methylene protons neighboring the nitrogen atom causes increased spectral complexity of the adjacent, following the methylene group (the second from the N-atom, see Figure 4c).

Effects of Alkyl Chain Length on FTIR Spectra of Solid NDIs
Full spectra are provided in Figure S2, and the positions of the absorption bands are listed in Table S7 in the Supplementary Materials. The NDIs' spectra are dominated by absorption bands characteristic of aromatic, alkyl, and imide moieties [39,63]. Selected characteristic bands are shown in Figure 5. Expectedly, the intensity of the bands corresponding to the motion of alkyl chains monotonically increases with the increasing number of carbon atoms in the alkyl chains. However, the variation in the alkyl-related bands is not only related to their intensity. Inspection of the spectra in the range of 2800 to 3000 cm −1 reveals an interesting transition: the position of the peak visible at υ~2940 cm −1 for short-alkyl chain-NDIs drops by δυ ≈ 15 cm −1 down to ~2925 cm −1 in approximately the middle range of the studied NDIs (Figure 5a). That transition can be likely attributed to a so-called gamma-gauche effect reported, for instance, for tetrathiafulvalene-based molecules [64]. Briefly, the gamma-gauche effect is related to the formation of specific packing of longer alkyl chains in regular crystal arrays. In order, however, to ultimately determine the nature of the transition observed here, a separate solid-state NMR study would be necessary [65,66]. The length of alkyl chains also exerts an influence on spectral features corresponding to the aromatic moieties of the NDI molecules. The absorption peaks related to aromatic C=C stretching (1560-1600 cm −1 Figure 5b) and aromatic out-of-plane bending (740-780 cm −1 Figure 5d) are subtly shifted towards the lower υ-numbers for NDIC3; these bands are also broader and less intense than observed in the NDIs with longer alkyl chains ( Figure 5). Another interesting and, to the best of our knowledge, sofar unreported feature is related to the C-N stretching bands (1320-1350 cm −1 , Figure 5c). The imide C-N stretching peak positions reveal a clear odd-even parity effect. Such effects related to the orientation of aromatic cores were recently reported for a homological series of double N,N′-bis(n-alkyl)-naphthalenediimides [33]. Here, in the case of symmetrical NDIs, the C-N stretching peaks corresponding to NDIs with odd numbered alkyl substituents are located at higher υ whereas those for the NDIs with an even number of carbon atoms in alkyl chains tend to be lower (Figure 5c). The odd-even δυ difference averages approximately 11 cm −1 ( Figure 5). There are two exceptions from the odd-even rule: (i) NDIC6, where the band maximum is located at υ higher than excepted, and (ii) NDIC3 that shows a noticeably different spectral profile within 1320 < υ < 1350 cm −1 corresponding to the imide C-N stretching band ( Figure 5). The variation in spectral features related to the aromatic cores of the NDIs can be explained on the basis of crystal packing discussed in the next section.

Effects of Alkyl Chain Length on FTIR Spectra of Solid NDIs
Full spectra are provided in Figure S2, and the positions of the absorption bands are listed in Table S7 in the Supplementary Materials. The NDIs' spectra are dominated by absorption bands characteristic of aromatic, alkyl, and imide moieties [39,63]. Selected characteristic bands are shown in Figure 5. Expectedly, the intensity of the bands corresponding to the motion of alkyl chains monotonically increases with the increasing number of carbon atoms in the alkyl chains. However, the variation in the alkyl-related bands is not only related to their intensity. Inspection of the spectra in the range of 2800 to 3000 cm −1 reveals an interesting transition: the position of the peak visible at υ~2940 cm −1 for short-alkyl chain-NDIs drops by δ υ ≈ 15 cm −1 down to~2925 cm −1 in approximately the middle range of the studied NDIs (Figure 5a). That transition can be likely attributed to a so-called gamma-gauche effect reported, for instance, for tetrathiafulvalene-based molecules [64]. Briefly, the gamma-gauche effect is related to the formation of specific packing of longer alkyl chains in regular crystal arrays. In order, however, to ultimately determine the nature of the transition observed here, a separate solid-state NMR study would be necessary [65,66]. The length of alkyl chains also exerts an influence on spectral features corresponding to the aromatic moieties of the NDI molecules. The absorption peaks related to aromatic C=C stretching (1560-1600 cm −1 Figure 5b) and aromatic out-of-plane bending (740-780 cm −1 Figure 5d) are subtly shifted towards the lower υ-numbers for NDIC3; these bands are also broader and less intense than observed in the NDIs with longer alkyl chains ( Figure 5). Another interesting and, to the best of our knowledge, so-far unreported feature is related to the C-N stretching bands (1320-1350 cm −1 , Figure 5c). The imide C-N stretching peak positions reveal a clear odd-even parity effect. Such effects related to the orientation of aromatic cores were recently reported for a homological series of double N,N -bis(n-alkyl)-naphthalenediimides [33]. Here, in the case of symmetrical NDIs, the C-N stretching peaks corresponding to NDIs with odd numbered alkyl substituents are located at higher υ whereas those for the NDIs with an even number of carbon atoms in alkyl chains tend to be lower (Figure 5c). The odd-even δ υ difference averages approximately 11 cm −1 ( Figure 5). There are two exceptions from the odd-even rule: (i) NDIC6, where the band maximum is located at υ higher than excepted, and (ii) NDIC3 that shows a noticeably different spectral profile within 1320 < υ < 1350 cm −1 corresponding to the imide C-N stretching band ( Figure 5). The variation in spectral features related to the aromatic cores of the NDIs can be explained on the basis of crystal packing discussed in the next section.

Crystal Structure of NDIs
Research performed over the past 30 years provided extensive crystallographic data on various NDIs [7,19,[55][56][57]67]. A survey of the related literature indicates that the research is mostly related to short-chained NDIs. Data for longer chain NDIs are less represented, so here we fill this gap. In order to complete the existing crystallographic dataset in the range from NDIC3 to NDIC11, we have determined the crystal structure of NDIC7, NDIC9, and NDIC11 (see chapter SI-4 in the Supplementary Materials). Our data, together with the crystal data of NDIC3 [55], NDIC4 [56], NDIC5 [57], NDIC6 [7], NDIC8 [19], and NDIC10 [67] reported previously in literature are collected in Table 2. The results acquired by the single-crystal X-ray diffraction technique were additionally compared with the crystallographic data determined from powder patterns ( Table 2, numbers in parentheses). Selected crystal data and structure refinement details of the above NDIs are extensively reported in the Supplementary Materials. Similarly, as observed in the FTIR spectra ( Figure 5), there is an interesting odd-even parity effect related to the alkyl chain length of the NDIs. The NDIs with an even number of methylene groups in the alkyl chains form typically triclinic unit cells with translational symmetry only (P1 ̅ , see Table 2). The crystals

Crystal Structure of NDIs
Research performed over the past 30 years provided extensive crystallographic data on various NDIs [7,19,[55][56][57]67]. A survey of the related literature indicates that the research is mostly related to short-chained NDIs. Data for longer chain NDIs are less represented, so here we fill this gap. In order to complete the existing crystallographic dataset in the range from NDIC3 to NDIC11, we have determined the crystal structure of NDIC7, NDIC9, and NDIC11 (see chapter SI-4 in the Supplementary Materials). Our data, together with the crystal data of NDIC3 [55], NDIC4 [56], NDIC5 [57], NDIC6 [7], NDIC8 [19], and NDIC10 [67] reported previously in literature are collected in Table 2. The results acquired by the single-crystal X-ray diffraction technique were additionally compared with the crystallographic data determined from powder patterns ( Table 2, numbers in parentheses). Selected crystal data and structure refinement details of the above NDIs are extensively reported in the Supplementary Materials. Similarly, as observed in the FTIR spectra ( Figure 5), there is an interesting odd-even parity effect related to the alkyl chain length of the NDIs. The NDIs with an even number of methylene groups in the alkyl chains form typically triclinic unit cells with translational symmetry only (P − 1, see Table 2). The crystals of NDIs substituted with odd-numbered alkyl chains reveal higher symmetrythey form typically P2 1 /c monoclinic cells ( Table 2). The odd-even parity effect also manifests itself in the π-π stacking distance (d π-π , see Table S8). The odd-numbered NDIs generally reveal slightly smaller d π-π values (3.2 Å on average) then the even-numbered ones (3.3 Å on average). In terms of crystal structure the NDIC3 is the most different NDI in the whole NDIC3-NDIC11 range. The NDIC3 has an orthorhombic (Pbca) unit cell and it was the only NDI in the study showing non-parallel stacking of the aromatic cores. That non-parallel arrangement explains the difference observed in the C-N stretching band in the FTIR spectrum of the NDIC3 ( Figure 5). The parameters of the unit cells and their formulae were used to Miller-index the powder patterns (see chapter SI-4.2 in the Supplementary Materials). Room-temperature crystal structures of NDI powders (and films) were found identical with the crystal structures of the single crystals. Table 2. Crystallographic data of NDIs determined from single crystal and powder X-ray diffraction patterns. Unit cell parameters determined from the powder patterns are given in brackets. The explanation for why d π-π for NDIC3 was not calculated is given in the text.

Crystal Morphologies of NDIs
One of the valid issues in modern technologies including fields of, among others, organic electronics, photovoltaics, or sensor systems is the fabrication of thin films. On a technological scale, the films are often fabricated by sublimation [21,68]. Nevertheless, however, a challenge is to fabricate films via scalable solution-based technologies like inkjet or screen printing or spray-or meniscus-guided coating [19]. It is, therefore, important to know the solubility-as provided here in Section 2.2-as well as to know the crystal morphologies resulting from isotropic crystallization-i.e., crystallization unaffected by any kind of external force. Here in Figure 6 we provide the polarizing optical morphologies of NDI films crystallized from chloroform at 25 • C on glass substrates. Generally, the NDIs in the study form two kind of morphologies: spherulitic/acicular for NDIs with shorter aliphatic chains and lamellar/feather-like morphologies for NDIs with longer aliphatic chains. NDIC3 and NDIC4 showed radiating domains with spherulitic morphology, whereas, in the case of NDIC4 on the images of the thicker part of the layer, needle and acicular structures could be also distinguished. Crystallization of NDIs C5, C6, and C7 resulted in fan-type sheets with fine lamellar textures. The NDIs C9, C10, and C11 show lamellar feather-like structures. An inspection of the micrographs in Figure 6 generally indicates that NDIs tend to dewet the glass leading to issues with film continuity. It seems that NDIC5, NDIC6, and NDIC7 formed the most continuous films after crystallization from chloroform.
Molecules 2023, 28, x FOR PEER REVIEW 9 of 21 C7 resulted in fan-type sheets with fine lamellar textures. The NDIs C9, C10, and C11 show lamellar feather-like structures. An inspection of the micrographs in Figure 6 generally indicates that NDIs tend to dewet the glass leading to issues with film continuity. It seems that NDIC5, NDIC6, and NDIC7 formed the most continuous films after crystallization from chloroform.

Thermal Properties of Solid NDIs
NDIs are relatively thermally stable molecules. TGA measurements indicated that non-oxidative thermal decomposition of NDIs begins above 300 °C-i.e., beyond the melting point. Onsets of thermal decomposition (T−3, T−5, and T−10, see Figure 7 and Table S9 in the Supplementary Materials) of the NDIs were found increasing along with the increasing number of carbon atoms in the alkyl chains. Interestingly, the onset of thermal decompositions of NDIC7 was negatively deviated from a generally linear trend. That negative deviation was, to some extent, qualitatively similar to a deviation of NDIC7's solubility in heptane (Figure 3a). The TGA profiles of NDIs (see Figure S5 in the Supplementary Materials) suggest that the thermal decomposition of NDIs is a single-step process.

Thermal Properties of Solid NDIs
NDIs are relatively thermally stable molecules. TGA measurements indicated that non-oxidative thermal decomposition of NDIs begins above 300 • C-i.e., beyond the melting point. Onsets of thermal decomposition (T −3 , T −5 , and T −10 , see Figure 7 and Table  S9 in the Supplementary Materials) of the NDIs were found increasing along with the increasing number of carbon atoms in the alkyl chains. Interestingly, the onset of thermal decompositions of NDIC7 was negatively deviated from a generally linear trend. That negative deviation was, to some extent, qualitatively similar to a deviation of NDIC7's solubility in heptane (Figure 3a). The TGA profiles of NDIs (see Figure S5 in the Supplementary Materials) suggest that the thermal decomposition of NDIs is a single-step process.
Further investigations of the NDIs' thermal properties were performed with differential scanning calorimetry (DSC). Some DSC studies of NDIs were already reported in literature [56,67,69]. In an attempt to find the general pattern of the phase transitions, the whole range of the NDIs in this work was analyzed in a unified way. The recorded raw experimental DSC curves are shown in Figure S6. Phase transition points and transition enthalpies are collected in Tables S10 and S11 where they are compared with the available literature data. In order to provide a clear overview of melting enthalpies and phase transition points for all NDIs, the data are graphically presented in Figure 8. The transition enthalpies were calculated from integrals over corresponding endo-or exotherms. Then, the calculated enthalpies were charted against temperature (Figure 8). The transition enthalpies in Figure 8 are marked in the shades of red or blue to distinguish between endothermic (red) or exothermic (blue) phase transitions. As indicated in the color scale bar, darker shades indicate higher enthalpy (Figure 8). Heating in the range from −20 • C to 250 • C causes the NDIs to reversibly transform two or three times before the melting (isotropization). Most notably, the isotropization temperature of NDIs decreases monotonically with increasing alkyl chain lengths (Figure 8). Unlike double N,N -bis(n-alkyl)-naphthalenediimides, and despite a regular odd-even pattern of crystal unit cells (see Table S4), NDIs do not show any kind of odd-even parity effect on melting or crystallization. However, and quite interestingly, the transformation points preceding the isotropization form two distinct trends. In the range of NDIC3-NDIC7, a trend roughly following the line from 100 • C to 150 • C is observed. In the range of NDIC7-NDIC11, another set of transitions following the line connecting −10 • C (for NDIC7) and 100 • C (for NDIC11) can be discerned. NDIC7 and NDIC8 revealed the most complex phase behavior. Based on the melting enthalpies and reversibility, the phase transitions preceding isotropization were identified as mostly crystal-crystal polymorphic transformations (Tables S10 and S11 in the Supplementary Materials). In some cases, such as NDIC4, the phase transitions preceding isotropization can be attributed to the crystal to liquid crystal transition [67].  Further investigations of the NDIs' thermal properties were performed with differential scanning calorimetry (DSC). Some DSC studies of NDIs were already reported in literature [56,67,69]. In an attempt to find the general pattern of the phase transitions, the whole range of the NDIs in this work was analyzed in a unified way. The recorded raw experimental DSC curves are shown in Figure S6. Phase transition points and transition enthalpies are collected in Tables S10 and S11 where they are compared with the available literature data. In order to provide a clear overview of melting enthalpies and phase transition points for all NDIs, the data are graphically presented in Figure 8. The transition enthalpies were calculated from integrals over corresponding endo-or exotherms. Then, the calculated enthalpies were charted against temperature (Figure 8). The transition enthalpies in Figure 8 are marked in the shades of red or blue to distinguish between endothermic (red) or exothermic (blue) phase transitions. As indicated in the color scale bar, darker shades indicate higher enthalpy (Figure 8). Heating in the range from −20 °C to 250 °C causes the NDIs to reversibly transform two or three times before the melting (isotropization). Most notably, the isotropization temperature of NDIs decreases monotonically with increasing alkyl chain lengths (Figure 8). Unlike double N,N′-bis(n-alkyl)-naphthalenediimides, and despite a regular odd-even pattern of crystal unit cells (see Table S4), NDIs do not show any kind of odd-even parity effect on melting or crystallization. However, and quite interestingly, the transformation points preceding the isotropization form two distinct trends. In the range of NDIC3-NDIC7, a trend roughly following the line from 100 °C to 150 °C is observed. In the range of NDIC7-NDIC11, another set of transitions following the line connecting −10 °C (for NDIC7) and 100 °C (for NDIC11) can be discerned. NDIC7 and NDIC8 revealed the most complex phase behavior. Based on the melting enthalpies and reversibility, the phase transitions preceding isotropization were identified as mostly crystal-crystal polymorphic transformations (Tables S10 and S11 in the Supplementary Materials). In some cases, such as NDIC4, the phase transitions preceding isotropization can be attributed to the crystal to liquid crystal transition [67].

Field-Effect Electron Mobility
Solution-processed, linear alkyl-substituted NDIs exhibit relatively high field-effect electron mobilities (µ e ), reportedly reaching almost 0.2 cm 2 V −1 s −1 [17,42,50,70,71]. The µ e for NDIC8 (studied in this work) can reach 0.16 cm 2 V −1 s −1, whereas longer-chain NDIs tend to reveal lower mobilities: µ e for NDIC12 and NDIC18 (not included in the NDI range in this work) were reported to be 0.01 and 0.005 cm 2 V −1 s −1 , respectively [42,72]. Despite a proven potential in the field of electronic and energy materials as n-type semiconductors [54], a survey of the literature reveals that a complete comparative study of their electronic properties is still lacking. Analysis of the scattered published data indicates, for instance, that the effect of the alkyl chain length cannot be unequivocally identified based on previous studies mainly because of different methodologies used in different works [51,70]. Hence, in our work, we performed unified tests of electron field mobility using organic field-effect transistors (OFETs) as diagnostic devices. For the sake of comparison, we have used completely unified and fully reproducible experimental procedures: the NDI films were deposited on commercially available OFET substrates with prepatterned gold electrodes. In order to avoid uncontrolled phase transitions and dewetting of NDI films (see Section 2.7) the OFETs were neither thermally nor solvent-vapor annealed. The solvent (chloroform) was removed by placing the OFET in a mild vacuum. Obviously, because of unequilibrated crystal morphology and a mismatch between gold work function (−5.1 eV) and HOMO/LUMO levels of NDIs (approx. −7.2 eV/−3.6 eV, respectively), our methodology did not allow reaching the highest possible µ e or the lowest possible threshold voltage (V th ) numbers [20,70,71,73]. On the other hand, however, our approach enabled a direct comparison between all NDIs in the studied range. The data collected in Table 3 suggest that the V th ranges from 40 to 55 V and reveal no clear trend. The recorded V th can be considered typical for molecular, NDI-like materials [70,71]. Unlike the V th , µ e numbers reveal the optimum range of the alkyl chain length. Comparison of the data in Table 3 indicates that the range of studied NDIs' µ e number reaches a maximum for NDIC6, NDIC7, and NDIC8. For these three NDIs, the µ e is in the range of 10 −4 cm 2 V −1 s −1 whereas, for the NDIs substituted with chains shorter or longer than in the abovementioned range, the µ e tends to be lower by one order of magnitude (Table 3).

Synthesis
In the study, we adopted the following synthetic protocol: a mixture of 1,4,5,8naphthalenetetracarboxylic dianhydride (4.0 g, 15 mmol), imidazole (8.2 g, 120 mmol), and appropriate amine (45 mmol in the case of n-pentylamine to n-undecylamine and 120 mmol of n-butylamine and n-propylamine) in toluene (75 mL) were stirred and heated under reflux at 150-155 • C (oil bath temperature) for about 12 h. The reaction mixture was cooled to room temperature and evaporated to dryness. The obtained residue was suspended in 96% ethanol (~150 mL) and gradually added to vigorously stirred, cold 2N hydrochloric acid (~300 mL). The resulting suspension was stirred for about 1 h and filtered directly through a Soxhlet thimble. The thimble, together with the crude product, was placed in a Soxhlet extractor body (200 mL, Quickfit ® (DWK Life Sciences Ltd., Stoke-on-Trent, UK)) equipped with a 500 mL round bottom flask and reflux condensers open to the air. In order to remove any colored impurities and starting materials, the crude product was extracted with methanol (300 mL) for 4-8 h. The heating was controlled to maintain the extraction rate at the level of approx. one drop of methanol per second. The extraction was continued 1 h beyond the point at which the washings became colorless (or almost colorless). Afterwards, the extraction solvent was changed to chloroform (300 mL) and the extraction was continued until all material soluble in chloroform was transferred to the flask. The chloroform solution of the NDI was then evaporated to dryness, and the residue was treated with methanol (50 mL), filtered, and dried in air. In the case of NDIC3 and NDIC4, additional purification via short-column chromatography and sublimation was required. Hence, NDIC3 and NDIC4 were additionally dissolved in a small amount of CH 2 Cl 2 , passed through a pad of silica gel (20 g, 70-230 mesh, Sigma-Aldrich), and finally washed with approx. 150 mL of CH 2 Cl 2 . Evaporation and sublimation (200 • C, 1.2·10 −2 mmHg) resulted in pure NDIC3 and NDIC4.

Infrared Spectroscopy
The infrared spectra (FTIR) were recorded for isotropic samples in the diamond ATR cell. In the study, we used the Bruker Vertex 70v Fourier transform IR spectrometer (Billerica, MA, USA). The spectra were collected as 64 scans at room temperature over the range of 4000 to 400 cm −1 with 2 cm −1 resolution. The measurements, the instrument control, and initial data processing were performed using OPUS software (v. 7.0 Bruker Optics, Ettlingen, Germany). The spectra were processed and visualized using OriginPro 2019.

UV-Vis Spectroscopy and Solubility Studies
The UV-Vis absorption spectra of the NDIs were measured in the chloroform (0.10 mM). The spectra were collected using a Perkin Elmer Lambda 20 spectrophotometer with a spectral resolution of 1 nm (Waltham, MA, USA). The NDIs were dissolved in chloroform and measured in quartz cuvettes with 10 mm or 1 mm path length at 20 • C. Data processing was performed using UV WinLab 2.70 software. For the purpose of data presentation, the UV-Vis spectra were exported to OriginPro 2019 as ASCII files. In order to test the solubility, 100.0 mg of each of the NDIs was dissolved in 1.0 mL of chloroform. For NDIs that dissolved completely, their solubility was considered as >100 mg/mL [49]. In case of incomplete dissolution, the solubility was determined quantitatively. The quantitative solubility determination was performed for selected NDIs in chloroform and all NDIs in n-heptane, toluene, and 1,2-dichlorobenzene. In order to determine the solubility, the NDIs were volumetrically dissolved in n-heptane, toluene, 1,2-dichlorobenzene, or chloroform to obtain stock solutions with concentrations of 0.040 mg/mL, 0.400 mg/mL, or 0.400 mg/mL, or 3.000 mg/mL, respectively. Then, the stock solutions were diluted to obtain five standard solutions. Concentrations of standard solutions were: 0.0200 mg/mL, 0.0133 mg/mL, 0.0100 mg/mL, 0.0080 mg/mL, and 0.0067 mg/mL in n-heptane; 0.2000 mg/mL, 0.1330 mg/mL, 0.1000 mg/mL, 0.080 mg/mL, and 0.0667 mg/mL in toluene and 1,2-dichlorobenzene; or 1.500 mg/mL, 1.000 mg/mL, 0.750 mg/mL, 0.600 mg/mL, and 0.500 mg/mL in chloroform. For each standard solution, a UV-Vis spectrum in the λ-range from 230 nm to 400 nm was measured. The absorbances at 376 nm (for n-heptane solutions), 281 nm (toluene), 363 nm (1,2-dichlorobenzene), or 265 nm (chloroform) were calibrated against the concentration. Saturated solutions were prepared by mixing NDIs (in amounts exceeding solubility) with 1 mL of each of the above solvents. The mixtures were stirred for 24 h at 25 • C. Afterwards, the saturated solutions were filtered through a 0.2 µm PTFE syringe filter and their absorbance at the corresponding wavelength was measured. The concentrations of the saturated solutions were determined from the absorbance using the calibration plots obtained for the corresponding standard solution series. The validity of the Lambert-Beer law was confirmed by measuring the absorbance of almost saturated solutions and checking whether they follow the calibration line. No deviations from the Lambert-Beer law were observed in the whole analyzed range.

NMR Spectroscopy
The 1 H and 13 C NMR spectra were measured on the 600 MHz Bruker Avance spectrometer in CDCl 3 (~15 mg/mL in the case of NDIC11 and~30-40 mg/mL in the case of other NDIs) as the solvent. Chemical shifts (δ, ppm) were determined relative to CDCl 3 signals (7.26 ppm for 1 H NMR and 77.16 ppm for 13 C NMR spectra).

Thermogravimetric Analysis and Differential Scanning Calorimetry
Thermogravimetric analysis (TGA) was performed using the Mettler Toledo (Columbus, OH, USA) TGA/DSC1 system, calibrated with indium and zinc as standards. The tests were carried out in a temperature range of 25-600 • C with a heating rate of 10 • C/min under an argon atmosphere at a flowrate of 60 mL/min. The analyzed samples (3-11 mg) were placed in the ceramic crucibles (70 µL). Differential scanning calorimetry (DSC) measurements were performed using the Mettler Toledo DSC1 coupled with the Huber TC 100 intracooler. The instrument was calibrated using indium standard (T m = 156.6 • C, ∆H m = 28.45 J/g). NDI samples (2-12 mg) were sealed in 40 µL aluminum pans and heated under a constant nitrogen purge (60 mL/min) from −60 • C to the temperature above the isotropization points of NDIs. After the isotropization, the samples were thermally equilibrated for 2 min and then cooled down to −60 • C. The samples were heated and cooled at a rate of 10 • C/min. Data processing was performed using the generic Mettler Toledo STAR e software and visualized using OriginPro 2019.

Polarized Optical Microscopy
All micrographs were taken with crossed polarizers using the Olympus BX60 (Tokyo, Japan) equipped with a ×10 magnification lens and a Panasonic DMC-G2 camera (Kadoma, Japan).

X-ray Diffraction Measurements
The single-crystal diffraction (SCD) data were collected on a New Xcalibur EosS2 diffractometer (Rigaku, Tokyo, Japan). The patterns were recorded using monochromatic MoK α radiation (λ = 0.71073 Å) for NDIC11 and NDIC7 and AgK α radiation (λ = 0.56087 Å) for NDIC9. The data were corrected for Lorentz polarization as well as for absorption effects [74]. The structures were solved by direct methods using the OLEX program with a plugin using the SHELXT program [75], which used an intrinsic phasing solution method, and refined with the SHELXL using least squares minimization [76]. All non-hydrogen atoms were refined anisotropically, while all hydrogen atoms were placed in idealized positions and refined as the 'riding model'. The isotropic NDIs samples were packed in 1 mm Mark capillaries (Capillary Tube Supplies Ltd., London, UK). The XRD measurements were performed using the Ganesha 300 XL+ (SAXSLAB ApS, Kopenhagen, Denmark). The instrument was equipped with a microfocus copper source, generic parabolic multilayer mirror optics, and two x/y JJ X-ray scatterless slits providing parallel monochromatic (λ = 1.5482 Å), point-collimated beams. The diffraction patterns were recorded on the Pilatus 300 k detector (Dectris AG, Baden-Daettwil, Switzerland) with 1800 s exposure. The camera length and, therefore, the 2θ-range were calibrated using silver behenate. Conversion to one-dimensional diffraction profiles (powder patterns) was performed by integration of 2D patterns over 180 • azimuthal angle (χ) using the Datasqueeze 3.0.9 program. Intensity normalization and background correction were performed numerically using OriginPro 2019.

Determination of Field-Effect Electron Mobility
Field-effect electron mobility was characterized using thin-film organic field-effect transistors (OFETs) as diagnostic devices. OFETs with bottom gate, bottom contacts configuration were used. The NDI active films were spin-coated from 5 mg/mL chloroform solutions at 2500 rpm for 30 s. In order to avoid phase transitions or sublimation, the OFETs were not thermally annealed or vacuum dried. All operations and measurements related to OFETs were performed under nitrogen atmosphere in the glove box system. The OFET measurements were performed with the Keithley 2614B dual-source meter (Tektronix Inc., Beaverton, OR, USA). Output characteristics were measured in the range from −10 to 80 V with 5 V gate-source step and 2 V drain-source step. The electron mobility was extracted from transfer characteristics in the saturation regime according to Equation (1).
where I DS is the drain current, W is the channel width, L is the channel length (see Section 3.1), C i is the capacitance of gate dielectric, and V g is the gate voltage.

Conclusions
Herein, we present a study on a homological series of NDIs N-substituted with linear aliphatic chains with 3 to 11 carbon atoms. The NDIs were analyzed spectroscopically with UV-VIS, NMR, and FTIR. The research was supplemented with studies on crystal structure and morphology performed by, respectively, X-ray diffraction (both powder and singlecrystal) and optical microscopy. Our systematic approach permitted unraveling interesting trends in the whole homological series as well as some peculiar features of specific NDIs. We started our studies with a determination of solubility in solvents commonly used in the fabrication of organic electronics films. It turned out that NDIs in the range of NDIC5-NDIC8 revealed the best of all solubility in chloroform, dichlorobenzene, or toluene. Extending the length of alkyl substituents beyond 8 carbon atoms caused a 2-or 3-fold decrease in the solubility in each of the solvents. Generally, the NDIs with shorter alkyl substituents also reveal solubility significantly lower than that observed in the mediumlength alkyl substituents. Observations on solubility already at this stage make the range of NDIC5-NDIC8 distinct from the viewpoint of practical applications. The effect of alkyl chain length on the solubility of NDIs in n-heptane is different. Extending the alkyl substituent chain length resulted in a gradual increase in the solubility in n-heptane from 0.005 mg/mL (NDIC3) to 0.600 mg/mL (NDIC11). However, the solubility of short-and long-alkyl-chained NDIs in n-heptane follow slightly different trends. The solubility of NDIs with substituents up to 6 carbon atoms increases faster with the number of carbon atoms. Starting from the the n-heptyl substituted NDIC7, we observed a different, linear trend in solubility with increasing length of the alkyl substituents. Our solution NMR studies indicated that the multiplicity of the methylene protons adjacent to the nitrogen atom is a non-first-order AA'XX' system and not a simple triplet, as is typically reported in literature. This observation suggests more complex than intuitively perceptible interactions of alkyl chains with imide moieties of the NDIs. The FTIR analysis of the solid NDIs revealed interesting spectral features related to the aromatic cores of the NDIs. The imide C-N stretching peak positions revealed a clear odd-even parity effect related to the length and number of methylene groups in the alkyl substituents. For the NDIs with odd numbered alkyl substituents, the C-N stretching peaks were found at wavenumbers distinctly higher than the peaks of the even-numbered NDIs. Further FTIR analysis revealed that NDIs can be divided into two distinct subgroups: the short-and the long-chained. The NDIs in the middle of the analyzed range (i.e., NDIC5-NDIC7) reveal some additional characteristic spectral features. Crystal packing of NDIs revealed a clear odd-even parity effect, thus explaining the spectroscopy findings. The NDIs with an even number of methylene groups in the alkyl substituents revealed low-symmetry (P − 1) triclinic packing, whereas those with an odd number of carbon atoms were generally monoclinic with P2 1 /c symmetry. The odd-even parity effect manifests itself in an overlapping of the NDIs' aromatic cores and the π-π stacking distance (d π-π ). The odd-numbered NDIs generally reveal slightly smaller d π-π values then the even-numbered ones. Solid NDIs revealed relatively high thermal stability that was found increasing from 320 • C to 420 • C with increasing substituent chain length. Interestingly, here again, some deviation from a linear trend was observed in the middle range of the alkyl chain length. The NDIs revealed remarkable patterns of thermally induced phase transitions. All the NDIs show a series of solid-state transitions prior to the isotropization. We observed two distinct patterns (kinds) of these pre-melting transitions for short-and long-chained NDIs. The NDIC7, exactly in the middle of the studied range, revealed transitions of both kinds. The NDIs from the middle of the analyzed range also stand out in terms of electronic properties. The field-effect electron mobilities determined for the NDIC6-NDIC8 range in normalized OFETs were 10-to 30-fold higher than for NDIs with shorter or longer alkyl substituents. All these observations may suggest that balanced aromatic and alkyl interactions are the key to achieve the best solution processability and electronic performance.

Supplementary Materials:
The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/molecules28072940/s1, UV-VIS spectra and supplementary UV-VIS data, FTIR spectra and supplementary FTIR data, solubility data, crystallographic analysis of NDIC7, NDIC9 and NDIC11, exemplary powder diffraction data, TGA and DSC data, NMR data and exemplary OFET characteristics are included in the supplementary file. All NMR spectra are provided in a separate file. New crystal structures of NDIC7, NDIC9, and NDIC11 were deposited in the CCDC database under deposition numbers, 2242898, 2242899, and 2242897, respectively (CCDC 2242897-2242899). Table S1: 13 C NMR data of the studied NDIs; Table S2: 1 H NMR data of the studied NDIs; Table S3: Wavelength (λ) and molar absorptivity (ε) of the main bands of NDIs; Tables S4-S6: Solubilities of NDIs in n-heptane, toluene and 1,2-dichlorobenzene; Table S7: Positions and intensities of the main FTIR absorption bands of the NDIs; Table S8: π-π stacking planes; Table S9: Thermal decomposition temperatures of NDIs: temperatures of 3, 5 and 10 % weight loss; Table S10: DSC data for the first heating and cooling cycles for the studied NDIs, together with literature data; Table S11: DSC data for the second heating cycle for the studied NDIs, together with literature data; Figure S1: UV-VIS spectra of the NDIs in the study; Figure S2: Series of FTIR spectra of all investigated NDIs; Figure S3: Unit cells of NDIC6 and NDIC7; Figure S4: Exemplary one-dimensional I(q) profiles of NDIC6: simulated in Mercury software and experimental; Figure S5: TGA curves of NDIs; Figure S6: DSC plots for the 1st heating and 1st cooling cycles, and 2nd heating cycle of NDIs; Figure S7: Exemplary output current-voltage characteristics for NDC7 [77].

Author Contributions:
The NDIs were synthesized, purified, and spectroscopically (UV-Vis, NMR) characterized by W.G. Solubility experiments were performed by W.G. and M.S. DSC, TGA, and powder X-ray diffraction measurements and analysis were performed by D.C. with support from A.K. The FTIR measurements were carried out by M.S. and D.C. The morphology analysis, OFET characterization, and OFET parameters' extraction were performed by K.J. The project was designed and supervised by W.G., A.K. and D.C. The manuscript was written by A.K. and W.G. with support from all the coauthors. Funding for the project was acquired by A.K. All authors have read and agreed to the published version of the manuscript.