Open Access Article
Iqbal
Javed
*abc,
Ayub
Khurshid
*de,
Muhammad Nadeem
Arshad
f and
Yue
Wang
*a
aState Key Laboratory of Supramolecular Structure and Materials, College of Chemistry, Jilin University, Changchun 130012, P. R. China. E-mail: javedkhattak79@gmail.com; yuewang@jlu.edu.cn; Tel: +46-704763801
bDepartment of Chemistry, Umeå University, SE-901 87 Umeå, Sweden
cDepartment of Chemistry & Biochemistry, University of Agriculture, Faisalabad, Pakistan
dDepartment of Chemistry, College of Science, King Faisal University, Al-Ahsa 31982, Kingdom of Saudi Arabia. E-mail: kayub@kfu.edu.sa
eDepartment of Chemistry, COMSATS Institute of Information Technology, Abbottabad, KPK, Pakistan 22060
fCentre of Excellence for Advanced Materials Research (CEAMR), King Abdulaziz University, P.O. Box 80203, Jeddah 21589, Saudi Arabia
First published on 27th November 2013
Quinacridone diimines 1–10 were synthesized by the condensation of anilines with alkyl substituted quinacridones (QA). Photophysical and electrochemical properties of the compounds were investigated. Unconventional behavior of absorption spectra suggested a decrease in π_conjugation within the QA skeleton as well as lack of extended π_conjugation between the QA skeleton and the N-phenyl rings. A computational study of compounds 1–10, a variable temperature 1H NMR study of compounds 2, 7 and 10 (for instance), and single crystal X-ray analysis of 2, 3, 6, 7, 8 and 10 indicated that the anomalous behavior is due to the buckled, non-planar structure of the quinacridones. Moreover the molecules existed in two interconvertible geometric isomeric forms at different temperatures. Molecular orbital calculations were performed at B3LYP/6-31+G(d), B3PW91/6-31G(d) and PBEPBE/6-31G(d) levels of theory at B3PW91/6-31G(d) optimized structures for both isomers of all compounds (1–10); the results obtained are in close agreement with the experimentally determined values.
The literature reveals that a great deal of effort has been devoted to the synthesis and characterization of various QA derivatives; however little attention has been paid to functionalization of the carbonyl group.10 To fully exploit the potential applications of such functional compounds in organic material based devices, it is necessary to design and synthesize some new types of QA derivatives with different molecular structures. A careful understanding of the relationship between molecular structures, conformations and optical and electronic properties in solution and condensed states will lead to the development of new strategies for the preparation of high-performance organic optical and electronic materials.
A series of QA diimines were synthesized by the condensation of alkyl substituted quinacridones with aromatic amines. The introduction of the aniline base backbone generally not only enhances the conducting and optoelectronic properties of the derivatives but also broadens the absorption band due to extended π-conjugation but this requires coplanarity, at least to within 15°.11 Physical properties are closely related to structure but after measuring the photophysical properties of the synthesized compounds the coplanarity was questionable as they show anomalous behavior in their physical properties.
For in-depth structural exploration of these compounds, single crystals have been grown at different temperatures for some compounds and single crystal X-ray analysis is carried out. Density Functional Theory (DFT) calculations are also performed in an attempt to attain a certain level of predictability in these systems. To get an unclouded picture we embarked on the study of the variable temperature 1H NMR spectrum of compounds 2, 7 and 10 for instance. We also evaluated the relationship between the structure, conformation and optoelectronic properties of the QA diimines.
On the basis of single crystal X-ray analysis, molecular modeling calculations and the variable temperature 1H NMR spectrum study, herein we report that these QA diimines exist in different geometrical isomeric forms at different temperatures with the buckled QA skeleton in a non-planar fashion.
The CnQA (n = 4,8), N,N-di(n-butyl)-1,3,8,10-tetramethylquinacridone (TMQA) and N,N-di(n-butyl)-2,9-difluoroquinacridone (DFDBQA) were synthesized according to the similar method reported in the literature.5
:
1). The desired product was obtained as a red solid in 50% yield.
The same procedure was adopted for compounds 2–10 and reaction conditions are tabulated in Table S1 (ESI†).
O and N–H, therefore the solubility of QA is increased. The soluble alkyl substituted QAs were treated with anilines to form a new series of QA derivatives. These compounds are soluble in common organic solvents, such as toluene, chloroform and tetrahydrofuran (THF) etc. The chemical structures were verified by 1H NMR, mass spectrometry and elemental analysis. Structures of some of the compounds were further confirmed by single crystal X-ray analysis.
It is worth noting that λmax of the majority of these compounds is blue shifted with respect to their parent non-iminated compounds (alkyl substituted QA), which is quite opposite to the expected red shifting as a result of extended π-conjugation. A generally accepted concept, for the relationship between structure and absorption, requires maximum planarity for producing the red shift because of increase of π-conjugation. From the observed absorption spectra it can be deduced that the N-phenyl ring might not be in the plane of the quinacridone nucleus, so it is not involved in extended π-conjugation. This conjecture was confirmed by single crystal X-ray analysis. Single crystal structure analysis of six compounds reveals that the QA nucleus is not planar but slightly buckles to a “butterfly” conformation. In addition, the N-phenyl rings are twisted out of the plane of the central ring system and prefer to be in conjugation with the nitrogen lone pair and not in extended π-conjugation. The QA skeleton is buckled to a large extent in compounds 5 and 6, which results in large blue shifts in the UV-Vis spectrum and this is confirmed by the crystal structure of 6 where the QA skeleton is buckled by 31.68° (Fig. 8). Compound 7 is the only molecule which exhibits a slight red shift in λmax by 17 nm as each N-phenyl ring is substituted by five fluorine atoms, which enhances intermolecular interactions.
Absorption spectra were also recorded in acetone and toluene in order to observe the effect of the polarity of solvents on absorption wavelengths (Fig. 2). Generally, QAs form more stable aggregates in less polar solvents and, as a result, the absorption maxima are red shifted. On the other hand, the absorption maxima are blue shifted in more polar solvents, probably due to the loss of aggregation. The solvent effect, i.e., the maximum absorption wavelength shifts to a longer wavelength with decreasing polarity, was not observed in regularity. This means that different substituents obviously affect the absorption of substituted QA.
Absorption spectra of these derivatives as thin films on quartz from CHCl3 solution were also recorded (Fig. 1). The λmax is slightly red shifted and this spectral displacement on going from solution to the solid state is in accordance with the molecular exciton theory of Kasha.13 It is also important to note that the spectral shape in the solid state is quite similar to that in solution and the numbers of absorption bands are the same indicating that the molecular nature is not changed.
![]() | ||
| Fig. 3 Quinacridone diimines in (a) sterically hindered form, (b) the “(A)” isomer, (c) the “(B)” isomer and (d) the monosubstituted product. | ||
QA carbonyl groups are less reactive towards nucleophilic attack than aliphatic carbonyl groups as these are the part of a conjugation system. The second carbonyl group left after mono-substitution was found to be more reactive than the first one. This is inferred from the fact that the mono-substituted product (Fig. 3d) was very difficult to obtain and the reaction generally yielded a mixture of starting material and the disubstituted product with very small or negligible amounts of mono-substituted product. A proposed explanation is that when the first imine is formed after mono-substitution, the QA system buckles and moves the second carbonyl group out of conjugation, making it more reactive than the first one. After disubstitution the QA nucleus is buckled and the N-phenyl ring is out of the plane of the QA nucleus. This buckling and non-planarity of the N-phenyl ring are evident in X-ray structures of molecules.
A variable temperature 1H NMR study of 2, 7 and 10 gives an evidence of different geometrical isomers (ESI†).
A variable temperature 1H NMR study of compound 2 revealed interesting molecular dynamics, one at lower temperature (200–230 K) and the other at higher temperatures (298–335 K). The 1H NMR spectrum of 2 (in chloroform-d3) at ambient temperature shows two broad singlets at δ 3.50 and 4.30, with a relative ratio of 1.0
:
0.58 and are attributed to ZZ(2B) (major) and EE(2A) (minor) isomers, respectively (vide infra). The broad signals become sharp when the temperature is reduced however at further lower temperatures, broadening of the signals is also observed (Fig. 5 and ESI†). The signal at δ 3.43 is broadened at 233 K and it splits at 223 K into two signals at δ 3.08 and 3.72 in a relative ratio of 1
:
1. The coalescence temperature is expected to be somewhere between 233 and 223 K. These new signals cannot be attributed to the EZ isomer because at a further lower temperature (213 K) the signal at δ 4.33 also splits in a similar fashion into two signals at δ 4.24 and 4.48 in a relative ratio of 1
:
1. The splitting of both peaks at lower temperatures indicates some process other than E–Z interconversion, mainly because of the fact that both isomers show a similar splitting pattern but different peaks appear in the NMR spectrum. The energy of the dynamic process is calculated to be 10.2 kcal mol−1 which is much lower than expected for the E–Z isomerization around a double bond.15 Since the signals due to N–CH2 are affected at lower temperature, we believe that the variable NMR behavior may be associated with some conformational process involved around N–CH2. At room temperature, the N-alkyl chains are free rotating however at lower temperature, free rotation of N-alkyl chains may be restricted due to interaction with the peri-hydrogens. Due to the restricted rotation, both hydrogens on N–CH2 (Ha and Hb) become diastereotopic due to the absence of the symmetry plane in the molecule, and behave differently in NMR measurements. Due to the restricted rotation, one alkyl chain is restricted above the plane of quinacridone whereas the other is restricted below the plane of the quinacridone. Some signals in the aromatic region (δ 7.5 and 8.4) are also broadened which may be assigned to the peri-hydrogens. We realized that the NMR signals at δ 3.46 and 4.34 are quite sharp at a temperature range of 263–283; however broadening is also observed at 303 (vide supra), which clearly points to some dynamic process just above room temperature. Since chloroform-d3 boils at 333 K, we have carried out high temperature NMR measurements in benzene. The 1H NMR spectrum of 2 (in benzene-C6D6) at ambient temperature shows two sharp singlets at δ 3.26 and 3.89, with a relative ratio of 1.0
:
0.58, which are attributed to ZZ(2B) (major) and EE(2A) (minor) isomers, respectively. However the dynamic processes in benzene proceed at different rates compared to chloroform. For example, the signals at δ 3.26 and 3.89 start broadening at 310 K and coalesced at 330 K. This indicates that the dynamic process is quite fast in benzene compared to the rate in chloroform. Since the EZ isomer is quite stable over the temperature range of ca. 320–335 K (the temperature range under study), it is a quite stable species on the potential energy surface. These findings are consistent with the results from the computational analysis where predicted relative energy of the EZ isomer is not very different from the EE isomer. The energy difference between ZZ(2B) and ZE is relatively large whereas the difference in relative energies between EZ and EE(2A) isomers is small (vida infra).
The variable temperature NMR behavior of compounds 7 and 10 is also analyzed (ESI†). The room temperature NMR spectrum (in benzene) of 7 shows a single peak at δ 4.38 which is a composite signal from the coalesced NCH2 groups of the isomer mixture due to rapid interconversion on the NMR time-scale at room temperature. These observations point towards the fact is that the ZZ–EE isomerization is much faster in 7 compared to 2. At ambient temperature the interconversion rate is fast on the NMR time scale. Occurrence of the ZZ(7B) isomer of compound 7 cannot be observed in the variable temperature NMR spectrum at any temperature under study. Not only the signal at δ 4.38 is broad at 298 K, but also the signals at δ 7.8–8.8. At 260 K, two broad singlets at δ 3.70 and 4.40 are also observed which coalesce at 275 K. The signals at δ 3.70 and 4.40 may be attributed to the EZ isomer because two signals are expected for the unsymmetrical isomer and they both coalesce at the same temperature (275 K) and are converted into a broad singlet at δ 4.38. The coalescence temperature is believed to be around 275 K. At further lower temperature (245 K), the peaks at δ 3.70 and 4.40 start splitting into two peaks with equal intensity. This process is very similar to the restricted conformation process described above for 2. However the splitting of the peak at δ 4.38 is not observed at temperatures as low as 230 K, which indicates that the restriction in conformational flexibility is achieved earlier for the Z component of the EZ isomer or for the ZZ(7B) isomer. This is very much consistent with the structural analysis. For the ZZ(7B) isomer (or the Z component of the EZ isomer), the phenyl ring of the imine part is bent towards the alkyl chain which pushes the latter towards the peri-hydrogens and therefore causes more steric hindrance for rotation. Therefore based on the structural analysis, the order of restricted conformational flexibility should be ZZ(B) > EZ > EE(A) and this is also observed experimentally. The high temperature NMR spectrum of compound 7 (in benzene) from 298–338 K is also shown in the ESI.† No significant change in the NMR of 7 could be observed in this temperature range, which indicates that the dynamic (E–Z) isomerization processes complete at 298 K (vide supra). Quite contrary to variable temperature of 2, no significant difference in the rates of dynamic processes is observed with the change of solvent from chloroform to benzene.
The variable temperature NMR spectrum of compound 10 is more complicated probably because in this case a number of conformational and geometrical isomerization processes are taking place simultaneously (ESI†). Compound 10 is structurally different than compounds 2 and 7 in two ways; it contains a long more flexible alkyl chain and a cyano group on the phenyl ring. The electron withdrawing effect of cyano probably speeds up the imino nitrogen inversion process and lowers the barrier to EZ isomerisation. The proton NMR spectral analysis of compound 10 reveals some flips of the chemical shifts. At room temperature, the signal at δ 3.54 is attributed to the EE(10A) isomer. EE–EZ isomerization is observed at 260 K, however the EZ isomer has two peaks, one at 3.04 and the other at 3.50, which is quite opposite to the behavior of 2 and 7 in variable temperature NMR. The protons due to the Z-component of the EZ isomer appear downfield compared to the E fragment, whereas in 2 and 7 a reverse phenomenon is observed. At lower temperature, each of the peaks split into two, which is attributed to restricted conformation (vide supra).
The rate constant, k, for the interconversion of the isomers at the coalescence temperature (Tc) 330 K for the broad NCH2 signals in benzene-d6 was calculated from the Gutowsky–Holm equation, and the imine interconversion barrier (ΔG‡) was calculated to be 15.42 kcal mol−1 by using the Eyring equation.
To further confirm the “(A)” and “(B)” isomers of QA diimines we grew crystals of these compounds at ambient and low temperatures. Importantly, we got single crystals of 3, 6, 8 and 10 at low temperature with “(B)” conformation and of 2 and 7 with “(A)” conformation at ambient temperature, which shows that QA diimines exist in two isomeric forms and more likely interconvert through intermediates.
In order to evaluate the stability of the EZ isomer relative to EE(A) and ZZ(B) isomers, we have calculated relative energies of all these isomers for compounds 7 and 10. For both these molecules, the EE isomer is the most stable isomer however the EZ isomer is unstable relative to the EE isomer by 2.55 and 1.54 kcal mol−1 for 7 and 10, respectively. The ZZ isomers for 7 and 10 lie relative to the EE isomer at 5.9 and 5.33 kcal mol−1, respectively (see Tables S6 and S7, ESI†).
Molecular orbitals of both isomers ((A) and (B)) are analyzed and MO energies and isodensity plots of HOMO − 2 to LUMO + 2 are shown in Fig. 6. Molecular orbitals of both isomers have strong resemblance among themselves regarding isodensity however their energies are quite different (see Table S3, ESI†). For 1a, the HOMO (−5.0 eV at B3PW91) is located primarily on the QA skeleton. The next highest occupied MOs (HOMO − 1 and HOMO − 2) are 0.43 and 0.53 eV below the HOMO and are centered mainly on the phenyl ring of the imine group. The lowest unoccupied orbital (LUMO) lies 2.99 eV above the HOMO. The LUMO and the next higher orbitals (LUMO + 1) are centered mainly on the QA core, however LUMO + 2 is exclusively based on the QA core. The internal orientation of the phenyl ring of imines has little effect on the isodensity of all of these orbitals.
A comparison of eigenvalues of HOMOs and LUMOs of QAs 1–10 (A) and (B) is shown in Table S3 (ESI†) at three different levels of theory. B3PW91/6-31G(d) and B3LYP/6-31+G(d) methods are quite good at predicting the energy of HOMOs (see Table S4, ESI†) in comparison with the experimental methods, which is consistent with the fact that hybrid functionals are better at predicting the HOMO eigenvalues than those without hybrid functionals (PBEPBE);16,17 although the former (B3PW91) slightly underestimates the energies, the latter (B3LYP) overestimates them. The simulated LUMO eigenvalues are much higher (for both hybrid and non-hybrid functionals) than the experimental values and this discrepancy is because of the poor description of virtual orbitals by time independent DFT methods. Hybrid functionals although significantly increase the accuracy of the HOMO eigenvalues, however they are less accurate in computing the HOMO–LUMO gap when compared with non-hybrid functionals PBEPBE (see Table S3, ESI†). The difference between theoretical and experimental HOMO–LUMO gaps decreases with increasing size of hydrocarbons.18 QAs are big enough to show the accurate HOMO–LUMO gap through calculations, however we think that the shape of the molecule and planarity are also key features. The QA skeleton is not even planar and it shows buckling effect.
The bond distances in this aromatic ring are shorter than the rg value for benzene (1.399 ± 0.001 Å) and this is also attributed to the fluorine substituent at the para position of the aromatic ring. It is also worth noting that the QA skeleton is not planar, but is buckled by 16.51(2)° and the angle between the QA body and the N-phenyl ring is 126.32(7)°.
As far as molecular packing is concerned molecules are connected to each other by N-phenyl rings through hydrogen bonding in such a way that the QA skeleton seems to be embedded perpendicularly between N-phenyl rings of different molecules.
![]() | ||
| Fig. 9 Crystal structure of compound 3 showing C–H⋯F hydrogen bonding interaction and unit cell packing. | ||
A contortion is observed in the N-phenyl ring, which is quite larger than the previous compounds due to the electron withdrawing character of the five fluorine atoms. The angles are slightly larger than expected for pure sp2 hybridization. The bond distances in this aromatic ring are shorter than the rg value of benzene (1.399 ± 0.001 Å) and this may be attributed to fluorine atoms substituted at the aromatic ring. In this molecule, the QA skeleton is also not planar, and is buckled by 20.52(2)° and the angle between the QA body and the N-phenyl ring is 125.71(8)°.
As far as molecular interactions are concerned each molecule is connected to four other molecules through two types of hydrogen bonds. Atoms C(14) and C(20) in the molecule at (x, y, z) act as hydrogen-bond donors, via atoms H(14) and H(20b) to atoms F(1) and F(3) at (1 + x, y, z) and (x, 1 + y, 1 + z) respectively to form the dimers (Fig. 12). The dH⋯A distance in C14–H14⋯F1 and C2–H20b⋯F3 is 2.59 Å and 2.53 Å respectively.
As far as molecular packing is concerned molecules are connected to each other by N-phenyl rings through hydrogen bonding in such a way that the QA skeleton seems to be embedded perpendicularly between N-phenyl rings of different molecules.
As this single crystal was grown at low temperature, it exists in the ‘ZZ(10B)’ form of conformation (Fig. 13). There is no π⋯π stacking or intermolecular hydrogen bonding involved between molecules. The only interactions which might be involved are van der Waals forces. The molecule is buckled by 22.74(1)° and the N-phenyl ring is twisted out of the plane of the central core by 51.5(9)°. This large buckling and twisting effect is due to the presence of the electron withdrawing ‘CN’ group at the N-phenyl ring. As far as molecular packing is concerned the central core of the molecules is located at the edges of the unit cell while the N-phenyl ring is embedded at the center of the unit cell.
Footnote |
| † Electronic supplementary information (ESI) available: 1H NMR spectra, cyclic voltammetric and crystallographic data. CCDC 923157 (2), 923158 (3), 923159 (6), 923160 (7), 923161 (8) and 923162 (10). For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/c3nj00477e |
| This journal is © The Royal Society of Chemistry and the Centre National de la Recherche Scientifique 2014 |