DOI:
10.1039/C9SC00962K
(Edge Article)
Chem. Sci., 2019,
10, 6482-6493
Doubly zwitterionic, di-reduced, highly electron-rich, air-stable naphthalenediimides: redox-switchable islands of aromatic–antiaromatic states†
Received
25th February 2019
, Accepted 21st May 2019
First published on 21st May 2019
Abstract
The di-reduced state of the naphthalene moiety and its congeners have long captivated chemists as it is elusive to stabilize these intrinsically reactive electron-rich π-systems and for their emergent multifaceted properties. Herein we report the synthesis and isolation of two-electron (2e−) reduced, highly electron-rich naphthalenediimides (NDIs). A doubly zwitterionic structure is observed for the first time in a naphthalene moiety and validated by single crystal X-ray crystallography and spectroscopic methods. The synthesis avoids hazardous reducing agents and offers an easy, high-yielding route to bench-stable di-reduced NDIs. Notably, we realized high negative first oxidation potentials of up to −0.730 V vs. Fc/Fc+ in these systems, which establish these systems to be one of the strongest ambient stable electron donors. The study also provides the first insights into the NMR spectra of the di-reduced systems revealing a large decrease in diatropicity of the naphthalene ring compared to its 2e− oxidized form. The NICS, NICS-XY global ring current, gauge-including magnetically induced current (GIMIC) and AICD ring current density calculations revealed switching of the antiaromatic and aromatic states at the naphthalene and the imide rings, respectively, in the di-reduced system compared to the 2e− oxidized form. Notably, the substituents at the phosphonium groups significantly tune the antiaromatic–aromatic states and donor ability, and bestow an array of colors to the di-reduced systems by virtue of intramolecular through-space communication with the NDI scaffold. Computational studies showed intramolecular noncovalent interactions to provide additional stability to these unprecedented doubly zwitterionic systems.
Introduction
Insight into the stability of multi-electron accumulated states is crucial to develop super-reductants, conductors, and electrode materials in batteries, solar cells, etc.1–4 Also, realization of super-electron rich systems is of value as their excited states can be applied towards metal-free organic transformations.5 However, π-conjugated systems exhibiting high negative oxidation potentials are intrinsically unstable under ambient conditions.6 There is also great interest to discern the effects of electron accumulation on the ring currents of these π-systems.7 As seen recently, intriguing electronic scenarios evolve if aromatic and antiaromatic states co-exist in π-systems.8 Understanding and control over these states may lead to realization of a rare class of bench-stable electron-rich π-systems. Presently, rational design concepts are lacking that can facilitate multi-electron injection, restrict their high reactivity and address their isolation challenges.
Isolation of a two-electron reduced naphthalene dianion and its π-extended systems is a challenge being pursued over last several decades. In 1965, a naphthalene dianion was realized in the solution state9a and later isolated in the solid state,9b while a π-extended acenaphthalene dianion10 was realized in solution (Scheme 1a). All these dianions are extremely reactive and are handled under inert conditions. Subsequently, arylenediimide π-scaffolds11 like naphthalenediimide (NDI) and perylenediimide (PDI) attracted huge interest due to their ability to accept two-electrons (2e−).12–14 While the 1e− reduced NDI radical ion has been isolated,15,16 isolation of the highly electron-rich di-reduced NDI remains elusive (Scheme 1a). In an elegant synthetic development, the 2e− reduced PDI dianion has recently been isolated.14 PDI with its π-extended structure can aid significantly in charge stabilization vis-à-vis the smaller NDI congener. Although through-space charge delocalization has been elegantly formulated in NDI cyclophanes, ambient stabilization of 2e− reduced individual NDIs has not been possible.17 Further, there have been considerable efforts to apply the photo-excited states of these electron-rich systems.18,19 Ground-state ambient stabilization and their stoichiometric usage is anticipated to open-up new vistas in controlling electron transfer reactions. We envisaged to explore the effect of 2e− uptake and release as well as the tunability of the aromatic naphthalene π-backbone and the nonaromatic imide rings. Thus, redox-assisted inter-conversion between the aromatic and nonaromatic states can be envisioned.
![image file: c9sc00962k-s1.tif](/image/article/2019/SC/c9sc00962k/c9sc00962k-s1.gif) |
| Scheme 1 (a) Naphthalene-based reported di-reduced species. (b) Pictorial and molecular representation of π-conjugated and betaine type di-zwitterions and (c) molecular structures of di-reduced naphthalenediimides. TMEDA = tetramethylethylenediamine. | |
Herein we report the first synthesis and isolation of highly electron-rich di-reduced NDIs. The doubly zwitterionic structure employed in the naphthalene scaffold affords an appealing mode of charge distribution and ambient stabilization (Scheme 1b). A facile eco-friendly synthetic protocol was rationalized from the electrochemical data and leads to strong electron donors (1–5) with the first oxidation potential ranging from −0.724 to −0.263 V vs. Fc/Fc+ (Scheme 1c). The single crystal X-ray structures for the di-reduced systems were compared with the corresponding 2e− oxidized dications.
Theoretical calculations showed that additional electrons are delocalized over the NDI unit and the carbonyl O atoms embrace bulk of the negative charge. On the other hand, the phosphonium groups contain majority of the positive charge. Such an arrangement of the opposite charges at two different sites within the molecule endows it with di-zwitterionic character. From the resonating structures, a hybrid of a π-conjugated zwitterion and a π-betaine zwitterion can be obtained as shown in Scheme 1b.
Results and discussion
Synthesis
Herein we report the first synthetic protocol towards the isolation of ambient stable NDI-based di-reduced compounds (1–5). We observed that second reduction potential of the dicationic NDI (2e− oxidized form) plays a key role during the synthesis of di-reduced compounds. EWGs at the axial- and the core-positions of the NDI led to a significant contrast between the first and second reduction potentials, which significantly favoured the one-step synthesis of 3–5. The reduction potentials and the corresponding LUMO levels of the parent dicationic compounds are listed in Scheme 2.
![image file: c9sc00962k-s2.tif](/image/article/2019/SC/c9sc00962k/c9sc00962k-s2.gif) |
| Scheme 2 Reduction potentials against Fc/Fc+ in DCM. The LUMO energies were calculated as ELUMO = [−5.1(Fc) − ER1] eV against vacuum on the basis of CV experiments. | |
We also realized a direct one-pot synthesis of 3 and 4 along with the corresponding radical ions when NDIs comprised of electron withdrawing groups (EWGs) at the axial positions (Scheme 3b). UV-vis absorption spectroscopy revealed the formation of di-reduced NDIs along with the radical ions (ESI Fig. S1†). The exclusive formation of 5 was observed when phosphenes were substituted with p-CF3Ph EWGs (Scheme 3c). Finally, we synthesized the corresponding dicationic compounds via one-electron oxidation of the radical ion or two-electron oxidation of the di-reduced compounds (Scheme 3d).
![image file: c9sc00962k-s3.tif](/image/article/2019/SC/c9sc00962k/c9sc00962k-s3.gif) |
| Scheme 3 Synthesis of di-reduced NDIs. Reagents and conditions: (i) PR3, Et3N, 100 °C, 30 min. and (ii) Na2S in MeOH, MeCN, RT, 30 min. (a) Stepwise synthesis. (b) In situ synthesis of radical and di-reduced species. (c) In situ synthesis of di-reduced 5 and (d) synthesis of dications from radical ions as well as di-reduced compounds. Reagents and conditions: (i) and (ii) NOBF4, RT, 15 min. | |
Importantly, 32+ and 52+ revealed exceptionally low LUMO levels of −5.15 and −5.03 eV, respectively. Interestingly, the LUMO of 32+ is 0.25 eV lower than that of the ultra-electron deficient NDI dication,16 and establishes itself as the lowest LUMO arylenediimide to be reported to date. Thus, 32+ is one of the strongest electron acceptors to be isolated.14,16a,20,21
Earlier we demonstrated the plausible mechanism for the in situ synthesis of diphosphonium substituted NDI radical ions.16b We anticipated that due to the contrast between the first and second reduction potentials in the dicationic systems embracing EWGs in the axial or core-positions, the radical ions would be able to accept an electron from the triethylamine to form in situ direduced compounds. To confirm this possibility, we performed a reaction of 12+ and 42+ with triethylamine at 100 °C. We observed gradual conversion of 42+ to 4, while 12+ only formed the corresponding radical ion (ESI Fig. S2†). This confirmed that triethylamine can further reduce the radical ion and form the direduced species if the dication is substituted with EWGs. On the basis of our earlier findings and this result, we proposed the probable mechanism for the formation of the direduced compounds (ESI Scheme S1†).
NMR spectroscopy
In the presence of an external magnetic field, aromatic systems sustain a diatropic ring current while antiaromatic systems sustain a paratropic ring current.22 The 1H NMR spectra of 1 shows a significant upfield chemical shift of all the protons compared to those of 12+ (Fig. 1). In the case of 1, the proton signal Ha corresponding to the proton of the naphthalene unit appears at 7.33 ppm (J = 17 Hz), which is shifted upfield by 1.61 ppm compared to the same protons of 12+ that appear at 8.94 ppm (J = 11 Hz). This signifies a significant contribution of the paratropic ring current in the naphthalene ring and decrease in diatropicity. This upfield 1H NMR chemical shift is observable even for the protons not directly attached to NDI scaffold. For example, the protons Hf and Hg of the ethyl groups linked to the phosphonium groups in 1 appear at 2.34 and 1.13 ppm, respectively, while in 12+ they appear at 2.75 and 1.30 ppm, respectively. The axial group protons, He (5.30 ppm) and Hb–d (7.28–7.15 ppm), of 1 also show measurable upfield chemical shifts compared to the same protons in 12+ appearing at 5.41 and 7.48–7.27 ppm, respectively.
![image file: c9sc00962k-f1.tif](/image/article/2019/SC/c9sc00962k/c9sc00962k-f1.gif) |
| Fig. 1
1H NMR spectra of (a) 1 and (b) 12+ in DCM at room temperature. The large upfield chemical shifts in 1 compared to 12+ are shown by red lines. (*) denotes solvent residual peaks of toluene (*, from crystallization), DCM (*) and water (*). | |
Such an upfield chemical shift in the 1H NMR spectra was also observed for 2–5. The phenylene/naphthalene protons in 2 (ESI Fig. S25†), 4 (ESI Fig. S36†) and 5 (ESI Fig. S42†) appear from 7.79, 7.78 and 7.92 ppm, which is significantly shifted upfield compared to 22+ (ESI Fig. S28†), 42+ (ESI Fig. S39†) and 52+ (ESI Fig. S44†), in which these protons begin from 8.42, 8.35 and 8.45 ppm, respectively. In 2, the CH protons of the phenyl diisopropyl axial groups are 0.31 ppm shifted upfield compared to the CH protons of 22+. Likewise, in 5 the methyl protons of the axial mesityl groups displayed 0.11–0.07 ppm upfield shifts compared to 52+. In 3 (ESI Fig. S31†) with a triethyl phosphonium group at the NDI-core, the phenylene/naphthalene peaks start from 7.50 ppm, which is 1.39 ppm upfield shifted compared to 32+ (ESI Fig. S34†), which appears at 8.89 ppm. The CH2 and CH3 protons of the phosphonium groups of 3 appear upfield shifted to 2.35 and 1.15 ppm respectively, compared to 2.75 and 1.31 ppm, respectively, in 32+.
The 31P NMR spectra of these di-reduced compounds also exhibited significant upfield chemical shifts of up to 11 ppm. Compounds 1, 2, 3, 4 and 5 showed 31P peaks at 32.90, 25.45, 33.25, 25.87 and 24.89, respectively. In the corresponding 12+, 22+, 32+, 42+ and 52+ this peak appears at 44.07, 31.09, 43.53, 31.12 and 30.45 ppm, respectively. The 13C NMR of the di-reduced compounds showed measurable upfield chemical shifts in the 13C signals compared to the dications. These significant upfield 1H, 31P and 13C NMR chemical shifts suggest efficient delocalization of the two additional electrons in the di-reduced systems.
FT-IR spectroscopy
The stretching frequencies of the imide carbonyl bonds of 1–5 are shifted by ∼50–120 cm−1 to lower frequencies compared to the dicationic compounds (ESI Fig. S11†). The imide carbonyl asymmetric and symmetric stretching frequencies for 1 were found to be 1648 cm−1 and 1617 cm−1, respectively, while for 12+ they were 1715 and 1667 cm−1, respectively. Likewise, for 2 the respective frequencies were 1630 and 1558 cm−1, respectively, compared to 1716 and 1664 cm−1 for 22+. Compound 32+ showed the corresponding frequencies of 1714 and 1664 cm−1 while for 3 it was observed to be 1619 and 1556 cm−1, respectively. The asymmetric and symmetric frequencies for 42+ and 52+ were found to be 1716 and 1658 cm−1 and 1720 and 1660 cm−1, respectively. They were found to be 1619 and 1558 cm−1 and 1634 and 1555 cm−1 for 4 and 5, respectively.
X-ray crystallography
To have the first insight into the structural characteristics of the di-reduced compounds, we sought to crystallize 1 and 2. X-ray suitable single crystals of 1 and 2 were grown from DCM
:
toluene (2
:
1) solution.
Further, single crystals of 12+ and 22+ were grown from DCM: toluene (2
:
1) solution at room temperature and the structural parameters compared with 1 and 2.23
Crystallographic studies of 1 and 12+.
1 crystallizes in the monoclinic space group “C2” with one molecule per asymmetric unit (Fig. 2a and b). We observed significant increment in the C
O bond length of 1 compared to 12+ and the longest C
O bond was found to be 1.252 Å, suggesting its partial double bond character. In addition, bond lengthening in the annular C5–C5′ and the transverse bonds of the naphthalene moiety along with shortening in the P1–C2, C2–C3, C1–C6 and C4–C7 bonds was observed indicating delocalization of the two additional electrons over the NDI scaffold (ESI Table S1†). The shortest intramolecular P⋯O non-bonded distance was found to be 2.806 Å, which is slightly shorter than that of 12+ (vide infra).24 The bond lengthening of the C
O groups and shortening of the P⋯O distance suggest significant polarization of the positive and negative charges over the P atom of the phosphonium group and the O atom of the imide groups, respectively. A π-conjugated betaine-type di-zwitterionic structure is therefore pertinent from the crystallographic insight.
![image file: c9sc00962k-f2.tif](/image/article/2019/SC/c9sc00962k/c9sc00962k-f2.gif) |
| Fig. 2 ORTEP representations of: [1] (a) top and (b) side view. Selected bond lengths (Å): C1–C2 1.456(4), C2–C3 1.371(4), C3–C4 1.423(4), C1–C5′ 1.405(4), C4–C5 1.409(4), C5–C5′ 1.437(4), C6–O1 1.252(4), C7–O2 1.244(4), N1–C6 1.394(4), N1′–C7 1.406(4), C2–P1 1.795(4), P1⋯O1, P1′⋯O1′ 2.806 and 2.818. ORTEP representations of [2]: (c) top and (d) side views. Selected bond lengths (Å): C1–C2 1.444(4), C2–C3 1.376(3), C3–C4 1.421(4), C1–C5′ 1.403(4), C4–C5 1.412(4), C5–C5′ 1.437(3), C6–O1 1.244(4), C7–O2 1.237(5), N1–C6 1.390(4), N1′–C7 1.415(5), C2–P1 1.779(3), P1⋯O1 2.658. ORTEP representations of (12+)·2BF4−: (e) top and (f) side views. Selected bond lengths (Å): C1–C2 1.396(2), C2–C3 1.413(2), C3–C4 1.373(2), C1–C5′ 1.407(2), C4–C5 1.410(2), C5–C5′ 1.410(2), C6–O1 1.212(2), C7–O2 1.219(2), N1–C6 1.392(2), N1′–C7 1.386(1), C2–P1 1.827(1), P1⋯O1 2.811, B⋯Ct 3.524, F1⋯C1–5 3.113–3.451. (g) Top view and (h) side view of crystal structure of (22+)·2BF4−. Selected bond lengths (Å): C1–C2 1.394(4), C2–C3 1.404(4), C3–C4 1.376(4), C1–C5′ 1.412(4), C4–C5 1.394(4), C5–C5′ 1.416(4), C6–O1 1.213(4), C7–O2 1.208(4), N1–C6 1.388(4), N1′–C7 1.404(4), C2–P1 1.824(4), P1⋯O1 2.738, B⋯Ct 3.565, F1⋯C1–5 2.968–4.010. H atoms have been removed for clarity. Thermal ellipsoids are shown at 50% probability. | |
12+ crystallizes in the triclinic space group “P
” with a half molecule per asymmetric unit (Fig. 2e and f). The imide C
O adjacent to the phosphonium group shows a bond distance of 1.212 Å, while the C
O bond situated remote to the phosphonium group had a distance of 1.219 Å. The intramolecular P⋯O non-bonded distance was found to be 2.811 Å. The BF4− counter anions show strong anion–π interaction with a small shift towards the imide region of the NDI scaffold. The distance between the F atom of anion and the C atoms of the naphthalene moiety is in the range of 3.113–3.451 Å and the distance from the centroid of naphthalene (Ct) to the B atom of the anion was found to be 3.524 Å.
Electrochemistry
Cyclic voltammetry (CV) and differential pulse voltammetry (DPV) experiments were performed for the di-reduced and dicationic forms in DCM (ESI Fig. S3 and S4†). The CV and DPV studies of 1–5 exhibited two reversible oxidation peaks corresponding to the sequential two-step ET process (Fig. 3). Notably, the first oxidation potentials of 1 and 3 are shifted to more negative values than that of the second reduction potentials of the corresponding dicationic compounds (Table 1).
![image file: c9sc00962k-f3.tif](/image/article/2019/SC/c9sc00962k/c9sc00962k-f3.gif) |
| Fig. 3 CV and DPV of di-reduced molecules (a) 1, (b) 2, (c) 3, (d) 4, (e) and 5. Conditions: 5 × 10−4 M in DCM; reference electrode, Ag/AgCl; working and auxiliary electrodes, Pt with 0.1 M Bu4NPF6 and (Fc/Fc+); 298 K; scan rate, 200 mV s−1. | |
Table 1 Oxidation potentials of 1–5 determined by CV studies against Fc/Fc+ in DCM. HOMO level energies and UV-vis absorption spectroscopya
Mol. |
Potential (V) vs. Fc/Fc+ |
Exp. λabsmax [nm] (ε, L mol−1 cm−1) |
E
OX1 (ΔEp) |
E
OX2 (ΔEp) |
HOMO |
ΔEp is the potential difference between the anodic and cathodic peaks. HOMO energies were calculated as EHOMO = [−5.1 (Fc) − EOX1] eV against vacuum on the basis of CV experiments.
|
1
|
−0.73 (0.08) |
−0.18 (0.08) |
−4.36 |
572 (21 200) |
2
|
−0.67 (0.11) |
−0.23 (0.11) |
−4.42 |
593 (21 600) |
3
|
−0.47 (0.11) |
0.06 (0.11) |
−4.62 |
575 (18 100) |
4
|
−0.55 (0.09) |
−0.12 (0.09) |
−4.54 |
584 (31 800) |
5
|
−0.49 (0.18) |
−0.12 (0.18) |
−4.61 |
591 (33 900) |
The most negatively shifted oxidation potentials at −0.724 and −0.186 V vs. Fc/Fc+ was observed for 1 corresponding to the 1st and 2nd oxidation potential, respectively. The HOMO was determined to be −4.376 eV from the CV experiment. Compound 2 showed oxidation potentials at −0.573 and −0.140 V corresponding to a HOMO of −4.527 eV. On the other hand, 3, 4 and 5 exhibited the 1st and 2nd oxidation potentials at −0.474 and 0.060 V, −0.554 and −0.122 V, and −0.263 and 0.095 V, respectively. The corresponding HOMO levels for 3–5 were found to be −4.626, −4.546 and −4.837 eV, respectively. Therefore, 1 is the strongest electron donor followed by 2, 4, 3 and 5. Notably, 1 is one of the strongest ambient stable electron donors to be isolated to date.25 Notably, in this class of electron donors the HOMO level can be modulated by ∼0.5 eV by altering the substituents at the axial and the core positions of the NDI scaffold.
CV and DPV studies of 12+–52+ demonstrate two quasi-reversible reduction peaks corresponding to sequential two-step electron transfer processes. The reduction potentials of compounds 12+–52+ clearly showed the effect of the EWGs at the axial- and the core-substituents (Scheme 2). On comparison of the reduction potential data, it is clear that EWG substituents at the NDI-axial groups play a major role in decreasing both the first and second reduction potentials, while EWGs at the NDI-core predominantly affect the second reduction potential and significantly reduces the ER1 − ER2 potential gap. Moreover, the plot of the square root of the scan rate with respect to the current peak of both the redox waves suggested its thermodynamic stability and reversibility (ESI Fig. S5†).
UV-vis spectroscopy
The UV-vis absorption spectra of the di-reduced compounds in DCM exhibit characteristic spectroscopic features of the NDI dianion. All di-reduced compounds exhibited intense absorption bands in the 412–592 nm regions with the absorption coefficient (ε) reaching 3.3 × 104 L mol−1 cm−1 (Fig. 4a and Table 1). Compound 1 absorbs at 412, 437, 523 and 571 nm, while 2 shows red-shifted absorption bands at 418, 440, 541 and 592 nm. The λmax of 2 undergoes a 21 nm red-shift compared to 1 suggesting the role of the axial and core substituents in substantially influencing the frontier molecular orbital energies. A similar trend is observed on going from 3–5. Fig. 4b shows the characteristics absorption bands of the other two oxidized states of 1, i.e. the radical ion 1˙+ and the dication 12+.
![image file: c9sc00962k-f4.tif](/image/article/2019/SC/c9sc00962k/c9sc00962k-f4.gif) |
| Fig. 4 (a) Normalized UV-vis-NIR spectra of 1–5 in DCM. (b) Spectra showing the three-state spectra of 1, 1˙+ and 12+. | |
Moreover, we could succeed in tuning the colour of the di-reduced compounds by changing the axial and core functional groups of the NDI (Fig. 4a). Compound 1 which has triethylphosphonium group at the NDI core and benzyl groups at the axial positions gives brownish-red colour, while on replacement of axial benzyl with p-CF3 benzyl in 3 the colour changes to orange-red.
On the other hand, 2 having phenyl phosphonium groups at the core with diisopropylphenyl as the axial substituents revealed a purple-blue colour, which can be further modulated to light purple colour in 4 and light blue in 5. Next, we followed the gradual formation of 4 from the radical ion 4˙+ in the presence of a mild reducing agent like tetrabutyl ammonium cyanide in DMF by UV-vis-NIR absorption spectroscopy (ESI Fig. S6†). The complete conversion to the di-reduced compound was confirmed by the complete disappearance of the characteristic absorption bands of the radical ion 4˙+ at 469, 630 and 724 nm and the appearance of the characteristic bands for 4 at 436, 499, 536 and 584 nm.
Stability in the solution state
The solution stability of a representative compound 4 was determined in different solvents through UV-vis-NIR absorption spectroscopy. The λmax absorption band was followed to calculate the half-life time (t1/2). In toluene
:
DCM (97
:
3 v/v), the t1/2 was found to be 6.1 months, which is the highest stability amongst any reported di-reduced NDI molecules (ESI Fig. S7†). In DCM the t1/2 was found to be 51.4 days. On the other hand, in polar aprotic solvents like THF, MeCN, DMF and DMSO the t1/2 ranges between ∼3 and 4 days. These results confirmed that in nonpolar solvents the di-reduced compound has significantly higher stability. It is to be noted that all the experiments were done under ambient conditions. The trace impurities including water present in polar solvents possibly reduce the stability of the direduced compounds. We also examined stability of 4 in the solid-state by exposing it to open air for several months and found that there are no significant changes in the absorption spectrum albeit formation of a negligible amount of the radical ion.
NICS computations
Next, we performed nucleus independent chemical shift (NICS) calculations. The negative NICS values indicate a diatropic ring current, while positive NICS values indicate a paratropic ring current.26 NICS (0) and NICS (1) values result from the calculated magnetic shielding tensor for a “dummy” atom located in the center of the ring and 1 Å above/below the plane of the ring system under examination, respectively. The NICS values of 1 and 2 are compared with the corresponding NICS values of 12+ and 22+, respectively (Table 2). NICS (0) and NICS (1) values of −8.3 and −10.3, and −7.4 and −9.6 ppm for the naphthalene ring of 12+ and 22+, respectively, indicate that the aromaticity is intact.
Table 2 NICS (0) and NICS (1) parameters for dication and di-reduced compounds
![image file: c9sc00962k-u1.tif](/image/article/2019/SC/c9sc00962k/c9sc00962k-u1.gif)
|
|
NICS (0), ppm |
NICS (1), ppm |
1
|
2
|
12+
|
22+
|
1
|
2
|
12+
|
22+
|
A |
4.0 |
4.9 |
−8.3 |
−7.4 |
0.8 |
1.6 |
−10.3 |
−9.3 |
B |
−4.3 |
−3.5 |
4.7 |
5.2 |
−5.9 |
−5.4 |
0.3 |
0.7 |
However, the positive NICS (0) values of 4.0 ppm and 4.9 ppm for the naphthalene ring in the case of the di-reduced systems 1 and 2, respectively, indicate mild antiaromatic nature. On the other hand, the NICS (1) values for the same was found to be 0.8 and 1.6 ppm respectively, which indicates the nonaromatic nature of the naphthalene moiety.
In contrast, the imide regions of 12+ and 22+ comprising of the six-membered rings shows positive NICS (0) values of 4.7 and 5.2 ppm, respectively, indicating antiaromaticity. The NICS (1) values for the same were found to be 0.3 and 0.7 ppm, respectively, indicating nonaromatic nature. However in 1 and 2, the NICS (0) turns negative with values of −4.3 and −3.5 ppm, indicating its weakly aromatic nature, while NICS (1) turns negative with with values of −5.9 and −5.4, suggesting its aromatic nature. These findings indicate that in the di-reduced compounds, the naphthalene moiety nurtures antiaromatic/nonaromatic states and the imide rings provide the aromatic/nonaromatic state. Thus, a 2e− redox reaction can switch the antiaromatic–aromatic sites in these electron-rich and electron-deficient systems.
NICS-XY scan
To comprehend the global ring current, we performed the NICS-XY scans for the di-reduced systems and compared with the corresponding dications.27 We performed NICS-XY scan for 1 and 12+ in different directions keeping NICS probes (dummy atoms) at 1.7 Å above the molecular plane under examination at 0.1 Å intervals and scanned along the X-axis (Fig. 5). The NICS-XY scan of 1 along the H′s of the naphthalene ring shows a NICSπ,zz value of −1.0 ppm at the centre, which indicates the nonaromatic nature of the naphthalene ring. It also shows a value of −4.9 ppm at the carbon centre bearing the H atom of naphthalene indicating a local diatropic current at these carbon centres (Fig. 5c). Further, the NICS-XY scan along the imide N atoms shows NICSπ,zz values of −1.0 ppm indicating the nonaromatic nature at the centre of the naphthalene ring. This imide N to imide N NICS-XY scan also shows a value of −10.8 ppm at the centre of the imide rings indicating significant diatropic current and the aromatic nature of the imide rings (Fig. 6a). The O to O NICS-XY scan shows a value of −6.6 ppm at the periphery of the imide ring also suggesting the aromatic nature of the imide rings (Fig. 5b). The NICS-XY scan of the dication molecule 12+ along one of the Hs of the naphthalene ring to another H of the same shows NICSπ,zz values of −20.8 ppm indicating a diatropic current within the boundaries, suggesting aromatic nature. The imide N to imide N and imide O to imide O NICS-XY scans also show NICSπ,zz values of −20.8 ppm indicating diatropic current at the centre and a value of −2.2 ppm at the imide ring which indicates the weakly antiaromatic nature of the imide ring. Although, NICS is a decent tool for interpretation of the aromaticity of a variety of cyclic compounds yet it has been criticized due to its unreliability to describe the aromaticity of the polycyclic structures.28 This is because NICS calculates the induced magnetic field at single point in space of a ring also affected by the σ-electrons while aromaticity is a property of global and local ring currents induced by a single ring and the rings associated with that ring.29
![image file: c9sc00962k-f5.tif](/image/article/2019/SC/c9sc00962k/c9sc00962k-f5.gif) |
| Fig. 5 NICS-XY scans of 12+ and 1. (a) Along imide N to imide N (solid; red = 1; blue = 12+), (b) imide O to imide O (hollow; red = 1; blue = 12+) and (c) along the diagonal naphthalene hydrogen; H to H (red = 1; blue = 12+). | |
![image file: c9sc00962k-f6.tif](/image/article/2019/SC/c9sc00962k/c9sc00962k-f6.gif) |
| Fig. 6 AICD isosurface plots of (a) 1, and (b) 12+. The induced ring current vectors are plotted on the AICD isosurface to denote the diatropic and paratropic ring currents. The AICD plots are plotted at an isosurface value of 0.05. | |
AICD calculations
To obtain further insight into the ring current density as well as delocalization of the π-electron clouds, we performed AICD calculations.30 According to the AICD calculations, a high critical isosurface value (CIV) reflects strong delocalization of the π-electrons, and a low CIV represents weak conjugation at a critical point.30b The CIVs have been assigned at the points of less conjugation for understanding the delocalization of electrons in both the dicationic and di-reduced systems. The AICD calculations of 12+ and 22+ show clockwise ring current on the naphthalene moiety indicating its aromatic nature, while in the imide region several of the current vectors were found to be disoriented which along with the anticlockwise current flow indicated antiaromatic nature of the imide rings (Fig. 6b and ESI Fig. S8b†).
The isosurface of the naphthalene annular bond of 12+ gets broken at a CIV of 0.073 and for 22+ it broke at 0.069 whereas, the side bonds rupture at 0.079–0.081 for 12+ and 0.076–0.099 for 22+. These values confirm their diatropic nature. In contrast, in the imide part the C–N bond ruptures at 0.054 to 0.063 indicating less delocalization. The CIV value of the carbonyl group adjacent to the phosphonium group was found to be 0.092 for 12+ and 0.089 for 22+ which is significantly lower than the CIV (0.120 and 0.113) value of the carbonyl group remote to the phosphonium group, respectively. This clearly suggests the bond lengthening of carbonyl groups adjacent to the phosphonium group due to the strong P⋯O non-bonding interactions as validated from crystallographic studies (Fig. 6b).
On the other hand, in 1 and 2 the ring currents in the naphthalene unit are oriented in the anti-clockwise direction, which suggests antiaromaticity. In the imide rings, several of the current vectors were found to be disoriented along with current vectors oriented in a clockwise manner suggesting weakly aromatic nature (Fig. 6a and ESI Fig. S8a†). The isosurface of the naphthalene annular bond of 1 and 2 ruptures at a CIV value of 0.071 and 0.059 while the peripheral bonds were found to break at 0.062–0.072 and 0.059–0.072, respectively. These low CIV values of the naphthalene unit of di-reduced compounds showed slightly less electron delocalization compared to the dicationic molecules. There is a small increment in the CIV in the imide region in one of the C–N bonds and a small decrease in the CIV of other imide bonds in the di-reduced molecules. Furthermore, the CIV value of carbonyl group adjacent to phosphonium group was found to be 0.059 and 0.056 for 1 and 2, respectively, while that of the carbonyl group remote to the phosphonium group was found to be 0.084 and 0.082. These CIV values are significantly lower than the CIV values found in the corresponding dicationic diphosphonium carbonyl groups. This corroborates lengthening in the bond lengths. Therefore, there is a swapping in the electron delocalization pattern between the naphthalene and the imide rings in both the dicationic and di-reduced systems.
Current density maps and integral analysis
The detailed information about the magnetically induced ring current densities,31 current strength and current pathways in both the direduced and dicationic compounds were investigated using the gauge-including magnetically induced current (GIMIC)32 method. The induced ring currents were visualized using ParaView 5.6.0.33 On careful examination of the ring currents it has been found that two weakly interacting ring currents; one corresponding to the diatropic ring current at the periphery of the molecule and one to the paratropic ring current in the center of the naphthalene ring of the NDI moiety circulate in both the direduced and dicationic compounds (Fig. 7 and ESI Fig. S9†).
![image file: c9sc00962k-f7.tif](/image/article/2019/SC/c9sc00962k/c9sc00962k-f7.gif) |
| Fig. 7 The magnetically induced current density of (a) 1 and (b) 12+ calculated in a plane placed at 1.0 Bohr (1 Bohr = 0.529 Å) above the molecular plane. A magnet was placed along the X axis. Diatropic currents are assumed to circle clockwise and the paratropic ones circle anticlockwise. | |
To estimate the accurate ring current strength and current pathways we performed integration analysis of the key bonds by placing a plane at each of the bonds. The integration analysis of the current density passing through the plane of the C1–C2, C2–C3, C3–C4 and C4–C5 bonds in 12+ shows a strong contribution of diatropic ring current of 16.37 to 17.97 nA T−1. It also shows a paratropic current contribution of −6.40 to −12.00 nA T−1 due to local current around carbon atoms and paratropic currents circling the naphthalene rings.
The integration of current strength passing a plane at the C5–C5′ bond yields a diatropic ring current of 13.04 nA T−1 and also a paratropic current contribution of 13.46 nA T−1. The net diatropic ring current of 4.74–11.6 nA T−1 in 12+ validates the aromatic nature of the naphthalene ring (Table 3). The integration analysis of the C1–C6 bond corresponding to the imide region shows diatropic and paratrpoic ring currents of 6.76 nA T−1 and −19.57 nA T−1, respectively. A net current of −12.81 nA T−1 confirms antiaromatic nature, while the C6–N1, C7′–N1 and C4′–C7′ bonds associated with the imide ring show a diatropic current of 7.24–7.87 nA T−1 and paratropic contribution of −6.84 to −7.37 nA T−1. This corresponds to a net current of ∼0.5 nA T−1 and suggests non-aromatic nature. The net current of the C
O bonds of the dication also shows antiaromatic nature (Table 3).
Table 3 The diatropic, paratropic and net current strengths (current strength susceptibility in nA T−1) calculated at the B3LYP/6311++G(d,p) level. The numbering of the molecules is given in the figure
![image file: c9sc00962k-u2.tif](/image/article/2019/SC/c9sc00962k/c9sc00962k-u2.gif)
|
Bond |
Diatropic |
Paratropic |
Net current |
1
|
12+
|
1
|
12+
|
1
|
12+
|
C1–C2 |
8.76 |
17.97 |
−9.31 |
−6.40 |
−0.54 |
11.56 |
C2–C3 |
8.43 |
16.37 |
−17.57 |
−13.34 |
−9.13 |
3.03 |
C3–C4 |
8.17 |
17.13 |
−10.63 |
−7.71 |
−2.45 |
9.42 |
C4–C5 |
8.88 |
17.02 |
−16.70 |
−11.98 |
−7.81 |
5.03 |
C1–C5′ |
8.49 |
16.74 |
−16.44 |
−12.00 |
−7.94 |
4.74 |
C5–C5′ |
7.08 |
13.04 |
−16.03 |
−13.46 |
−8.95 |
−0.42 |
C1–C6 |
9.36 |
6.76 |
−17.92 |
−19.57 |
−8.56 |
−12.81 |
C4′–C7′ |
9.94 |
7.37 |
−4.45 |
−6.84 |
5.49 |
0.52 |
C6–N1 |
11.31 |
7.87 |
−5.73 |
−7.37 |
5.57 |
0.49 |
C7′–N1′ |
10.67 |
7.24 |
−5.24 |
−6.89 |
5.43 |
0.34 |
C6–O1 |
11.00 |
12.81 |
−18.36 |
−18.63 |
−7.35 |
−5.82 |
C7–O2 |
11.00 |
12.03 |
−16.05 |
−16.83 |
−5.04 |
−4.80 |
C2–P1 |
7.54 |
7.97 |
−6.51 |
−6.53 |
1.02 |
1.44 |
On the other hand, the diatropic ring current in C1–C2, C2–C3, C3–C4, C4–C5 and C5–C5′ bonds of 1 was found to be considerably moderated with values of 7.08–8.88 nA T−1, while the paratropic current contribution was found to be strong with a value of −9.31 to −17.57 nA T−1. The net current was found to be −0.54 to −8.95 nA T−1, which corroborates that this naphthalene ring possesses non-aromatic regions to antiaromatic regions (Table 3). The diatropic and paratropic ring contributions at the plane of the C1–C6 bond were found to be 9.36 nA T−1 and −17.92 nA T−1. This corresponds to a net current of −8.56 nA T−1, which suggests antiaromaticity at this position. However the C6–N1, C7′–N1 and C4′–C7′ bonds, which are also the part of the imide ring, yield diatropic 9.94–11.31 nA T−1 and paratropic −6.84 to −7.37 nA T−1 current contribution.
This leads to a net diatropic current with a value of 5.43 to 5.57 nA T−1 and suggests weak aromaticity in the imide ring in 1. We also calculated the current contributions for 2 and 22+ which are given in ESI Table S3.† The integration analysis also validated the switching of the aromatic and antiaromatic regions in the dicationic and the direduced compounds.
In essence, NICS, NICS-XY scan, AICD plots and GIMIC calculations suggest that there is a significant degree of switching in the aromaticity/antiaromaticity/nonaromaticity character in the naphthalene and the imide rings as a result of the two additional electrons in the di-reduced systems.
FMO, ESP and NPA calculations
The ground state geometry optimization of 1 and 2 and 12+ and 22+ was performed by using the DFT/B3LYP method with the 6-311++G(d,p) basis set using the IEFPCM model in DCM. The axial groups in all the cases were replaced by the methyl group to save computational cost. The calculated HOMO and LUMO energies of 1 were found to be −4.145 and −1.436 eV, respectively, while for 2 they were −4.225 and −1.628 eV, respectively. For 12+, the HOMO and LUMO levels were −8.372 and −4.903 eV, respectively, and for 22+ they were −8.063 and −4.894 eV, respectively. The calculated HOMO energy levels of di-reduced and LUMO energy levels of dications are in excellent correlation with electrochemically obtained HOMO–LUMO values (ESI Fig. S10†). The HOMO surface was found to be delocalized over the NDI moiety including the carbonyl oxygen in both 1 and 2. The LUMO of 1 was predominantly delocalized over the NDI while in 2 it delocalizes over the two phosphonium groups with minor delocalization over the NDI unit. In the case of 12+ and 22+, the LUMO was found to be delocalized over the NDI π-surface, while the HOMO was found to be delocalized over the phenyl phosphonium groups in 22+ and over the NDI π-surface in 12+. The theoretically calculated HOMO and LUMO energy levels of 4, 42+, 5 and 52+ were also found to be in good agreement with the experimental data (ESI Fig. S8†).
The 3D electrostatic surface potential (ESP) distribution of 22+ and 2 further supports the above findings. In 22+, a highly electron deficient region over the NDI surface is evident, while 2 shows a highly electron rich region delocalized over the NDI surface and also the triphenyl rings of the phosphonium groups (Fig. 8).
![image file: c9sc00962k-f8.tif](/image/article/2019/SC/c9sc00962k/c9sc00962k-f8.gif) |
| Fig. 8 ESP maps of 2 and 22+. The contours are color-coded from red (electron-rich) to blue (electron-deficient). | |
Natural population analysis (NPA)34 also revealed existence of higher negative charge on the oxygen atoms adjacent to phosphonium groups in the dicationic as well as in the di-reduced compounds compared to the oxygen atoms remote to phosphonium groups. Likewise, the phosphorus atoms retain most of the positive charge (ESI Table S4†).
NBO analysis
Natural bond orbital (NBO) analysis35 shows the orbital interactions between the non-bonding orbitals of the imide oxygen atoms and the antibonding orbitals of the P–C bonds
(ESI Fig. S11†). A donor–acceptor type interaction with a stabilization energy of E(2) = 2.24 kcal mol−1 for 22+ and E(2) = 4.44 kcal mol−1 for 2 per P⋯O interaction was estimated using second-order perturbation theory. This stronger P⋯O interaction in 2 further supports greater distribution of the negative charge through the imide carbonyl groups due to the additional two electrons in the di-reduced compounds.
AIM analysis
Atoms in molecule (AIM) analysis36 also validate strong P⋯O interaction. An analysis of the topological properties of the electron density and determination of the bond critical point (BCP), which is accompanied by a bond path between the P atom of the phosphonium group and O atom of the imide carbonyl (P⋯O) was performed for 1, 2 and its dicationic precursors.
The 2D critical bond paths are shown in Fig. 9 and the corresponding 3D contour plots are shown in ESI Fig. S12†. The computed electron density parameters at the BCP, ∇2ρ(rb) P⋯O, E(rb)P⋯O, G(rb)P⋯O, and V(rb)P⋯O for the P⋯O interactions are summarized in Table 3. The electron density, ρ(rb), values at the BCP of the P⋯O are in the range of H-bonding interactions (ρH-bond ≈ 0.002–0.04 a.u.).31 The ρ(rb) for 1 (0.021 a.u.) and for 2 (0.025 a.u.) are higher than that for 12+ (0.018 a.u.) and 22+ (0.018 a.u.) validating a stronger P⋯O interaction in the di-reduced compounds. The Laplacian density, ∇2ρ(rb), is also in the range of strong H-bond interactions.
![image file: c9sc00962k-f9.tif](/image/article/2019/SC/c9sc00962k/c9sc00962k-f9.gif) |
| Fig. 9 2D contour plots of (a) 1 and (b) 12+ showing the bond path between P and O. | |
The positive 0.000 value of E(rb) for molecule 22+ indicates its dominant electrostatic nature and the negative E(rb) values for 2 indicate its partial covalent nature [E(rb) value at a specific BCP indicates whether the interaction is electrostatic dominant E(rb) > 0 or covalent dominant E(rb) < 0]. The positive value of G(rb) and negative value of V(rb) support its stable, bound stationary states (Table 4). Inter-locked non-covalent interactions can thus play a significant role in electron delocalization.37
Table 4 AIM parameters for the dications and di-reduced compoundsa
|
Bond |
ρ(rb)b |
∇2ρ(rb)c |
E(rb)d |
G(rb)e |
V(rb)f |
Theoretical topological properties at the BCP; values from the 6-311++G(d,p) basis set.
The electron density at the BCP.
Laplacian of electron density.
The electron energy density.
Kinetic energy electron density.
Potential energy electron density. All are in a.u.
|
1
|
P⋯O |
0.021 |
0.051 |
−0.001 |
0.013 |
−0.014 |
12+
|
P⋯O |
0.018 |
0.050 |
−0.000 |
0.012 |
−0.013 |
2
|
P⋯O |
0.025 |
0.055 |
−0.001 |
0.015 |
−0.016 |
22+
|
P⋯O |
0.018 |
0.051 |
0.000 |
0.012 |
−0.012 |
Conclusions
We demonstrated isolation of two-electron (2e−) reduced, highly electron-rich, bench-stable NDIs. A doubly zwitterionic structure is observed in a naphthalene moiety and validated by single crystal X-ray crystallography. This new genre of compounds was achieved in high-yields via a solvent-free synthetic protocol. The di-reduced systems endow exceptional inherent stability with a half-life time of more than six months in toluene under ambient conditions. Notably, high negative first oxidation potentials up to −0.730 V vs. Fc/Fc+ and HOMO levels extending to −4.360 eV were realized.
The study offers the first insights into the NMR spectra of the di-reduced systems revealing large upfield chemical shifts for the NDI-core atoms and the adjoining groups compared to their 2e− oxidized form, suggesting a large decrease in the diatropicity of the naphthalene ring. The NICS, NICS-XY, GIMIC and AICD calculations revealed redox switching of the antiaromatic and aromatic states at the naphthalene and the imide rings, respectively, in the di-reduced system compared to the corresponding rings of the 2e− oxidized form. The substituents at the axial- and core-positions tune the antiaromatic–aromatic states and the electron donor ability. Further an array of diverse colors was achieved for the di-reduced systems. The through-space non-covalent interactions are key elements assisting this tunability. Isolation of this new class of ambient stable systems should have fascinating implications in controlling electron transfer reactions and development of new switchable materials.
Conflicts of interest
There are no conflicts to declare.
Acknowledgements
PM acknowledges financial support under Swarna Jayanti Fellowship (DST/SJF-02/CSA-02/2013-14), DST-FIST and DST-PURSE. PM thanks AIRF, JNU for the instrumentation facilities. S. K., J. S. and K. M. acknowledge UGC for their research fellowship. We thank S. Subrahmanyam, K. Reddy and Prof. V. G. Anand, IISER Pune, for their help in AICD calculations. We thank Prof. Dr R. Herges and F. Roehricht, Otto-Diels-Institut für Organische Chemie Christian-Albrechts-Universität zu Kiel Otto-Hahn Germany, for their help with AICD calculations. We thank Prof. A. Stanger, Technion–Israel Institute of Technology, for his help in AROMA calculations. We also thank H. Fliegl, R. Bast, S. G. Patra for their help in GIMIC calculations and K. Moreland for help with the ParaView visualization.
Notes and references
-
(a) A. B. Buades, V. S. Arderiu, D. Olid-Britos, C. Viñas, R. Sillanpää, M. Haukka, X. Fontrodona, M. Paradinas, C. Ocal and F. Teixidor, J. Am. Chem. Soc., 2018, 140, 2957 CrossRef CAS PubMed;
(b) A. V. Zabula, A. S. Filatov, S. N. Spisak, A. Yu. Rogachev and M. A. Petrukhina, Science, 2011, 333, 1008 CrossRef CAS PubMed.
-
(a) S. Goswami, J. N. Nelson, T. Islamoglu, Y.-L. Wu, O. K. Farha and M. R. Wasielewski, Chem. Mater., 2018, 30, 2488 CrossRef CAS;
(b) J. A. Murphy, J. Org. Chem., 2014, 79, 3731 CrossRef CAS PubMed.
-
(a) J.-i. Yamada, H. Akutsu, H. Nishikawa and K. Kikuchi, Chem. Rev., 2004, 104, 5057 CrossRef CAS PubMed;
(b) C. Rovira, Chem. Rev., 2004, 104, 5289 CrossRef CAS PubMed.
-
(a) J. E. Coughlin, Z. B. Henson, G. C. Welch and G. C. Bazan, Acc. Chem. Res., 2014, 47, 257 CrossRef CAS PubMed;
(b) D. P. Dubal, O. Ayyad, V. Ruiza and P. Gómez-Romero, Chem. Soc. Rev., 2015, 44, 1777 RSC;
(c) H. Huang, L. Yang, A. Facchetti and T. J. Marks, Chem. Rev., 2017, 117, 10291 CrossRef CAS PubMed;
(d) M. Stępień, E. Gońka, M. Żyła and N. Sprutta, Chem. Rev., 2017, 117, 3479 CrossRef PubMed.
-
(a) I. Ghosh, T. Ghosh, J. I. Bardagi and B. König, Science, 2014, 346, 725 CrossRef CAS PubMed;
(b) I. Ghosh, L. Marzo, A. Das, R. Shaikh and B. König, Acc. Chem. Res., 2016, 49, 1566 CrossRef CAS PubMed.
-
(a) N. L. Holy, Chem. Rev., 1974, 74, 243 CrossRef CAS;
(b) A. Rajka, A. Safronov, S. Rajka and J. Wongsriratanakul, J. Am. Chem. Soc., 2000, 122, 3351 CrossRef;
(c) Y. Matsuo and E. Nakamura, J. Am. Chem. Soc., 2005, 127, 8457 CrossRef CAS PubMed;
(d) T. Kupfer, H. Braunschweig and K. Radacki, Angew. Chem., Int. Ed., 2015, 54, 15084 CrossRef CAS PubMed;
(e) H. Oshima, A. Fukazawa, T. Sasamori and S. Yamguchi, Angew. Chem., Int. Ed., 2015, 54, 7636 CrossRef CAS PubMed.
-
(a) T. M. Krygowski, H. Szatylowicz, O. A. Stasyuk, J. Dominikowska and M. Palusiak, Chem. Rev., 2014, 114, 6383 CrossRef CAS PubMed;
(b) R. Gershoni-Poranne and A. Stanger, Chem. Soc. Rev., 2015, 44, 6597 RSC;
(c) Z. Zeng, X. Shi, C. Chi, J. T. L. Navarrete, J. Casado and J. Wu, Chem. Soc. Rev., 2015, 44, 6578 RSC.
-
(a) T. Sarma, G. Kim, S. Sen, W.-Y. Cha, Z. Duan, M. D. Moore, V. M. Lynch, Z. Zhang, D. Kim and J. L. Sessler, J. Am. Chem. Soc., 2018, 140, 12111 CrossRef CAS PubMed;
(b) T. Y. Gopalakrishna, J. S. Reddy and V. G. Anand, Angew. Chem., Int. Ed., 2014, 53, 10984 CrossRef CAS PubMed;
(c) C. K. Frederickson, L. N. Zakharov and M. M. Haley, J. Am. Chem. Soc., 2016, 138, 16827 CrossRef CAS PubMed.
-
(a) J. Smid, J. Am. Chem. Soc., 1965, 87, 655 CrossRef CAS;
(b) J. J. Brooks, W. Rhine and G. D. Stucky, J. Am. Chem. Soc., 1972, 94, 7346 CrossRef CAS.
- L. D. Kershner, J. M. Gaidis and H. H. Freedman, J. Am. Chem. Soc., 1972, 94, 985 CrossRef CAS.
-
(a) M. Al Kobaisi, S. V. Bhosale, K. Latham, A. M. Raynor and S. V. Bhosale, Chem. Rev., 2016, 116, 11685 CrossRef PubMed;
(b) C. Huang, S. Barlow and S. R. Marder, J. Org. Chem., 2011, 76, 2386 CrossRef CAS PubMed;
(c) S.-L. Suraru and F. Würthner, Angew. Chem., Int. Ed., 2014, 53, 7428 CrossRef CAS PubMed.
-
(a) G. Andric, J. F. Boas, A. M. Bond, G. D. Fallon, K. P. Ghiggino, C. F. Hogan, J. A. Hutchison, M. A.-P. Lee, S. J. Langford, J. R. Pilbrow, G. J. Troup and C. P. Woodward, Aust. J. Chem., 2004, 57, 1011 CrossRef CAS;
(b) S. Guha and S. Saha, J. Am. Chem. Soc., 2010, 132, 17674 CrossRef CAS PubMed;
(c) V. V. Roznyatovskiy, D. M. Gardner, S. W. Eaton and M. R. Wasielewski, Org. Lett., 2014, 16, 696 CrossRef CAS PubMed;
(d) M. R. Ajayakumar, G. Hundal and P. Mukhopadhyay, Chem. Commun., 2013, 49, 7684 RSC.
-
(a) E. Shirman, A. Ustinov, N. Ben-Shitrit, H. Weissman, M. A. Iron, R. Cohen and B. Rybtchinski, J. Phys. Chem. B, 2008, 112, 8855 CrossRef CAS PubMed;
(b) M. A. Iron, R. Cohen and B. Rybtchinski, J. Phys. Chem. A, 2011, 115, 2047 CrossRef CAS PubMed.
- S. Seifert, D. Schmidt and F. Würthner, Chem. Sci., 2015, 6, 1663 RSC.
- G. Belangér-Chabot, A. Ali and F. P. Gabbaï, Angew. Chem., Int. Ed., 2017, 56, 9958 CrossRef PubMed.
-
(a) S. Kumar, M. R. Ajayakumar, G. Hundal and P. Mukhopadhyay, J. Am. Chem. Soc., 2014, 136, 12004 CrossRef CAS PubMed;
(b) S. Kumar and P. Mukhopadhyay, Green Chem., 2018, 20, 4620 RSC.
-
(a) S. T. Schneebeli, M. Frasconi, Z. Liu, Y. Wu, D. M. Gardner, N. L. Strutt, C. Cheng, R. Carmieli, M. R. Wasielewski and J. F. Stoddart, Angew. Chem., Int. Ed., 2013, 52, 13100 CrossRef CAS PubMed;
(b) Y. Wu, M. Frasconi, D. M. Gardner, P. R. McGonigal, S. T. Schneebeli and M. R. Wasielewski, Angew. Chem., Int. Ed., 2014, 53, 9476 CrossRef CAS PubMed.
- N. T. La Porte, J. F. Martinez, S. Chaudhuri, S.-H. Victor, S. Batista and M. R. Wasielewski, Coord. Chem. Rev., 2018, 361, 98 CrossRef CAS.
-
(a) N. T. La Porte, J. F. Martinez, S. Hedström, B. Rudshteyn, B. T. Phelan, C. M. Mauck, R. M. Young, V. S. Batista and M. R. Wasielewski, Chem. Sci., 2017, 8, 3821 RSC;
(b) D. Gosztola, M. P. Niemczyk, W. Svec, A. S. Lukas and M. R. Wasielewski, J. Phys. Chem. A, 2000, 104, 6545 CrossRef CAS.
-
(a) S. Kumar, J. Shukla, Y. Kumar and P. Mukhopadhyay, Org. Chem. Front., 2018, 5, 2254 RSC;
(b) M. Kivala and F. Diederich, Acc. Chem. Res., 2009, 42, 235 CrossRef CAS PubMed;
(c) A. Nowak-Król, K. Shoyama, M. Stolte and F. Würthner, Chem. Commun., 2018, 54, 13763 RSC;
(d) H. Zhylitskaya and M. Stępień, Org. Chem. Front., 2018, 5, 2395 RSC;
(e) F. N. Miros and S. Matile, ChemistryOpen, 2016, 5, 219 CrossRef CAS PubMed.
-
(a) J. Míšek, A. V. Jentzsch, S.-i. Sakurai, D. Emery, J. Mareda and S. Matile, Angew. Chem., Int. Ed., 2010, 49, 7680 CrossRef PubMed;
(b) J. Chang, Q. Ye, K.-W. Huang, J. Zhang, Z.-K. Chen, J. Wu and C. Chi, Org. Lett., 2012, 14, 2964 CrossRef CAS PubMed;
(c) W. Yue, W. Jiang, M. Böckmann, N. L. Doltsinis and Z. Wang, Chem.–Eur. J., 2014, 20, 5209 CrossRef CAS PubMed;
(d) Y. Kumar, S. Kumar, S. K. Keshri, J. Shukla, S. S. Singh, T. S. Thakur, M. Denti, A. Facchetti and P. Mukhopadhyay, Org. Lett., 2016, 18, 472 CrossRef CAS PubMed.
-
(a) R. Breslow, Acc. Chem. Res., 1973, 6, 393 CrossRef CAS;
(b) K. B. Wiberg, Chem. Rev., 2001, 101, 1317 CrossRef CAS PubMed.
- CCDC 1822819 (1), CCDC 1822817 (12+), CCDC 1822818 (2) and CCDC 1822820 (22+) contain the supplementary crystallographic data for this paper.†.
- S. S. Batsanov, Inorg. Mater., 2001, 37, 871 CrossRef CAS.
- J. Broggi, T. Terme and P. Vanelle, Angew. Chem., Int. Ed. Engl., 2014, 53, 384 CrossRef CAS PubMed.
-
(a) P. v. R. Schleyer, C. Maerker, A. Dransfeld, H. Jiao and N. J. R. v. E. Hommes, J. Am. Chem. Soc., 1996, 118, 6317 CrossRef CAS PubMed;
(b) Z. Chen, C. S. Wannere, C. Corminboeuf, R. Puchta and P. V. Schleyer, Chem. Rev., 2005, 105, 3842 CrossRef CAS PubMed.
- R. Gershoni-Poranne and A. Stanger, Chem.–Eur. J., 2014, 20, 5673 CrossRef CAS PubMed.
-
(a) R. Gershoni-Poranne and A. Stanger, Chem. Soc. Rev., 2015, 44, 6597 RSC;
(b) P. Lazzeretti, Prog. Nucl. Magn. Reson. Spectrosc., 2000, 36, 1 CrossRef CAS.
-
(a) P. Lazzeretti, Phys. Chem. Chem. Phys., 2004, 6, 217 RSC;
(b) J.-i. Aihara, Chem. Phys. Lett., 2002, 365, 34 CrossRef CAS;
(c) A. Stanger, Chem. Commun., 2009, 1939 RSC;
(d) N. S. Mills and K. B. Llagostera, J. Org. Chem., 2007, 72, 9163 CrossRef CAS PubMed;
(e) T. K. Dickens, R. B. Mallion and S. Radenković, J. Phys. Chem. A, 2019, 123, 1445 CrossRef CAS PubMed;
(f) P. B. Karadakov, P. Hearnshaw and K. E. Horner, J. Org. Chem., 2016, 81, 11346 CrossRef CAS PubMed.
-
(a) R. Herges and D. Geuenich, J. Phys. Chem. A, 2001, 105, 3214 CrossRef CAS;
(b) D. Geuenich, K. Hess, F. Köhler and R. Herges, Chem. Rev., 2005, 105, 3758 CrossRef CAS PubMed;
(c) H. Fliegl, J. Jusélius and D. Sundholm, J. Phys. Chem. A, 2016, 120, 5658 CrossRef CAS PubMed.
-
(a) D. Sundholm, H. Fliegl and R. J. F. Berger, Wiley Interdiscip. Rev.: Comput. Mol. Sci., 2016, 6, 639 CAS;
(b) H. Fliegl, D. Sundholm, S. Taubert, J. Jusélius and W. Klopper, J. Phys. Chem. A, 2009, 113, 8668 CrossRef CAS PubMed;
(c) S. Pelloni and P. Lazzeretti, J. Phys. Chem. A, 2013, 117, 9083 CrossRef CAS PubMed.
-
(a) J. Jusélius, D. Sundholm and J. Gauss, J. Chem. Phys., 2004, 121, 3952 CrossRef PubMed;
(b) H. Fliegl, S. Taubert, O. Lehtonen and D. Sundholm, Phys. Chem. Chem. Phys., 2011, 13, 20500 RSC;
(c) S. Taubert, D. Sundholm and J. Jusélius, J. Chem. Phys., 2011, 134, 054123 CrossRef PubMed;
(d) D. Sundholm, R. J. F. Berger and H. Fliegl, Phys. Chem. Chem. Phys., 2016, 18, 15934 RSC.
-
(a)
J. Ahrens, B. Geveci and C. Law, ParaView: An End-User Tool for Large Data Visualization, Visualization Handbook, Elsevier, 2005, ISBN-13: 978-0123875822 Search PubMed;
(b)
U. Ayachit, The ParaView Guide: A Parallel Visualization Application, Kitware, 2015, ISBN 978-1930934306 Search PubMed.
- A. E. Reed and R. B. Weinstock, J. Chem. Phys., 1985, 83, 735 CrossRef CAS.
-
E. D. Glendening, J. K. Badenhoop, A. E. Reed, J. E. Carpenter, J. A. Bohmann, C. M. Morales, C. R. Landis and F. Weinhold, NBO 6.0, Theoretical Chemistry Institute, University of Wisconsin, Madison, WI, 2013 Search PubMed.
-
R. F. W. Bader, Atoms in Molecules: A Quantum Theory, Claren-don Press, Oxford, U.K., 1990 Search PubMed.
- H. Huang, L. Yang, A. Facchetti and T. J. Marks, Chem. Rev., 2017, 117, 10291 CrossRef CAS PubMed.
Footnote |
† Electronic supplementary information (ESI) available. CCDC 1822817–1822820. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/c9sc00962k |
|
This journal is © The Royal Society of Chemistry 2019 |