Unexpected course of reaction between (E)-2-aryl-1-cyano-1-nitroethenes and diazafluorene: why is there no 1,3-dipolar cycloaddition?

Abstract Reactions between (E)-2-aryl-1-cyano-1-nitroethenes and diazafluorene lead to acyclic 2,3-diazabuta-1,3-diene derivatives, instead of the expected pyrazoline systems. DFT calculations suggest that this is a consequence of formation of zwitterionic structure in the first stage of the reaction. It must be noted that this is a specific property of the (E)-2-aryl-1-cyano-1-nitroethenes group, in contrast to most other conjugated nitroalkenes. Graphical abstract Electronic supplementary material The online version of this article (doi:10.1007/s00706-016-1893-5) contains supplementary material, which is available to authorized users.

(E)-2-Phenyl-1-cyano-1-nitroethene and their aryl-substituted analogs (ACNE) were prepared for the first time in the first half of the twentieth century [7]. At present, several compounds from this group are known [8][9][10][11]. However, their chemical properties are not well known. Some compounds have been described very recently. In particular, some examples of participation of ACNE in Diels-Alder reactions as dienophiles [12][13][14][15] as well as heterodienes [8,16,17] are explored. Additionally, some examples of catalyzed 1,3-dipolar cycloadditions between ACNE and trimethylsilyl azide are also analyzed [18,19]. Unfortunately, this method does not yield stable adducts, because the primary reaction products decompose partially under reaction conditions. Actually, any examples of thermal (non-catalyzed) 1,3-dipolar cycloadditions involving ACNE are unknown. This work is an attempt to fill this gap and is a continuation of our systematic study about participation of CNA in cycloaddition reactions [15,[20][21][22][23][24]. In particular, we have decided to shed light on the reactions between the homogenous series of ACNE (1a-1d) and diazafluorene (2), as a model allenyl-type 1,3dipole with [(C -)-(N ? ):N functional group. Theoretically, these reactions should give attractive, from a practical point of view, nitrofunctionalized pyrazoline Electronic supplementary material The online version of this article (doi:10.1007/s00706-016-1893-5) contains supplementary material, which is available to authorized users.
& Radomir Jasiński radomir@chemia.pk.edu.pl systems (Scheme 1). In addition to experimental studies, we also performed comprehensive quantum chemical studies to understand better the nature and molecular mechanism of these reactions.

Results and discussion
Firstly, we prepared the reaction components. For this purpose, we applied methodologies described earlier (see the experimental section). Next, the analysis of the interaction between the addends in the expected cycloaddition course was performed. Global and local electronic properties of (E)-2-aryl-1-cyano-1-nitroethenes had been analyzed in detail previously [17,25]. It was discovered that all these compounds were characterized by a high global electrophilicity (in the case of compounds 1a-1d and were in the range of 3.14-3.68eV). In comparison (Table 1), diazafluorene has evidently weaker electrophilic nature (x = 1.84 eV). Additionally, 2 is characterized by a relatively high global nucleophilcity (more than 3.5 eV). In consequence, the interaction between CNA 1a-1d and 2 may be considered a polar one [26]. Subsequently, the local electronic properties of addends have been analyzed. As established earlier, in all CNA 1a-1d, most of the electrophilic center is located at the b-position of the nitrovinyl moiety [17,25]. On the other hand, most of the nucleophilic center in the [CNN moiety of 2 is located on the terminal nitrogen atom. In polar cycloadditions, a reaction course is controlled by the nature of the local nucleophileelectrophile interaction. Therefore, we have assumed that the reaction channel A should be preferred.
To verify the quantum chemical simulations, we performed experimental tests of the reactions of interest. In the first step, we explored reactions involving ACNE 1a.
It was found that this reaction proceeds in MeNO 2 solution under mild conditions and yields a dark brown solid. HPLC analysis of the post-reaction mixture shows the existence of one reaction product, which was isolated by crystallization from ethanol. Its constitution was established by means of elemental analysis as well as spectral techniques. It was found unexpectedly that the results of the elemental analysis were fundamentally different from those of the expected adduct. Next, in the IR spectrum, any bands from NO 2 as well as CN groups were not observed. On the other hand, the MS spectrum gives a molecular ion, which suggests that the molecular weight is less than that in the expected adduct. This suggest the absence of the C(CN)NO 2 moiety in the molecule. Thus, we established that it is not a five-membered heterocyclic product, but a 2,3-diazabuta-1,3-diene derivative 5a. This was fully confirmed by the single crystal X-ray structure determination ( Fig. 1 as well as Supplementary Material). For a full characteristic of this compound, 1 H and 13 C NMR spectra were also recorded (see ''Experimental'' section).
This phenomenon can be explained when assuming that the primary reaction product 3a-3d spontaneously decomposed with: C(CN)NO 2 carbene elimination according to the retro-[4 ? 1]-cycloaddition scheme (path C on Scheme 2). Theoretically, the highly substituted fivemember heterocycles may have decomposed with the ring Scheme 1 opening via carbene elimination. It should be highlighted, however, that any cases of this type of decomposition of nitropyrazoline systems had not been previously described. Some of these heterocycles decomposed under mild conditions, but via completely different mechanisms [28][29][30][31][32]. Alternatively, it may be assumed that in the first reaction stage, a zwitterionic intermediate I is formed. The zwitterionic structure of I is probably stabilized by a pushpull electronic effect, which is determined by the presence of two EWG groups (NO 2 and CN) in the terminal position. A similar effect has been recently explored [11] in detail in the case of a molecule, which has similar structural moieties. In the next step, it is converted via C(CN)NO 2 carbene elimination into a 2,3-diazabuta-1,3-diene derivative (Scheme 2). It should be noted at this point that the possibility of the existence of zwitterionic structures on the paths of reactions involving CNA has been recently described in the case of interactions between diarylnitrones and 1-EWG-substituted 1-nitroethenes [33,34] as well as between thiocarbonylylides and nitroethene [35,36].
To confirm this hypothesis, we have performed DFT simulation of theoretically possible paths for a model reaction involving CNA 1c. It was found that, in contrast to the reaction 1c ? 2 ? I, both cycloadditions leading finally to pyrazoline systems (3c and 4c) should be treated as forbidden from a kinetic point of view. In particular, for reactions 1c ? 2 ? 3c(4c), Gibbs free energy of activation is higher than 141 kJ/mol, whereas in reaction 1c ? 2 ? I about 128 kJ/mol. In consequence, our

Scheme 2
Unexpected course of reaction between (E)-2-aryl-1-cyano-1-nitroethenes and diazafluorene:… 911 calculations may support mechanism via paths D ? E, which is illustrated in Scheme 2. The key physical parameters of all considered transition states are presented in Fig. 2. It should be noted that the zwitterionic nature of I, as well as the charge distribution on its molecule, was confirmed by a population Mulliken analysis (see the GEDT value). Finally, it should be noted that the described reaction course is a specific property of the ACNE group. In comparison, a similar reaction between (E)-3,3,3-trichloro-1nitroprop-1-ene (6) (which has similar electrophilicity to ACNE [22]) and diazafluorene proceeds in MeNO 2 solution under mild conditions and gives pyrazoline 7 with quantitative yield (Scheme 3). Interestingly, this adduct is easily decomposed into 8 under relatively mild conditions. This decomposition proceeded via the HCl elimination stage.

Experimental
The melting points were determined on a Boetius apparatus and are uncorrected. Elemental analyses were performed on a Perkin-Elmer PE-2400 CHN apparatus. IR spectra were recorded on a Bio-Rad spectrophotometer in CCl 4 solution. The 1 H NMR (500 MHz) and 13 C NMR (125 MHz) spectra were recorded on a Bruker AMX 500 spectrometer. Liquid chromatography (HPLC) was done using a Knauer apparatus equipped with a UV/Vis detector. For monitoring of the reaction progress, LiChrospher 18-RP 5 lm column (4 9 240 mm) and 75% methanol as the eluent at a flow rate of 1.0 cm 3 /min were used.
The experimental details and final atomic parameters for 5a have been deposited with the Cambridge Crystallographic Data Centre as supplementary material (CCDC ID 1448932). Copies of the data can be obtained free of charge on request via http://www.ccdc.cam.ac.uk/conts/retrieving. html (or from the Cambridge Crystallographic Data Centre, 12, Union Road, Cambridge CB2 1EZ, UK; fax: ?44 1223 336033).

Reactions between CNA and diazafluorene: general procedure
A mixture of appropriate nitroalkene (1a-1d or 6, 0.012 mol) and diazafluorene (0.010 mol) in 5 cm 3 of the appropriate solvent was stirred at room temperature for 12 h. The solvent was evaporated in vacuo to dryness and the semisolid residue was recrystallized firstly from ethanol and then from cyclohexane.

Quantum chemical calculations
All calculations reported in this paper were performed on ''Zeus'' supercomputer in the ''Cyfronet'' computational center in Cracow. Global and local electronic properties were estimated on the basis of structures obtained-according to Domingo suggestions [47,48]-on the basis of B3LYP/6-31G(d) calculations. For this purpose, we have used structures created by the standard procedure [49]. In particular, the electronic chemical potentials (l) and chemical hardness (g) were evaluated in terms of oneelectron energies of FMO (E HOMO and E LUMO ) using the equations: Next, the values of l and g were then used for the calculation of global electrophilicity (x) [47,48] according to the formula: x ¼ l 2 =2g; and the global nucleophilicity (N) [50] can be expressed in terms of the equation: N ¼ E HOMOðdiazafluoreneÞ À E HOMOðtetracyanoetheneÞ : The local nucleophilicity (N k ) [51] condensed to atom k was calculated using global nucleophilicity N and Parr function P k - [52] according to the formula: For the simulation of the reaction paths, hybrid functional B3LYP with the 6-31 ??G(d), basis set included in the GAUSSIAN 09 package [53] was used. It was found previously that calculations using B3LYP functional illustrate well the structure of transition states in polar 1,3-dipolar cycloadditions involving conjugated nitroalkenes [22,23,30]. The critical points on the reaction paths were localized in an analogous manner as in the case of the previously analyzed cycloadditions of diaryldiazomethanes with nitroacetylene [54]. In particular, for structure optimization of the reactants and the reaction products, the Berny algorithm was applied. First-order saddle points were localized using the QST2 procedure. The transition states were verified by diagonalization of the Hessian matrix and by an analysis of the intrinsic reaction coordinates (IRC).
All calculations were carried out for the simulated presence of nitromethane as the reaction medium. For this purpose, PCM [55] was used. For optimized structures, the thermochemical data for the temperature T = 298 K and pressure p = 1 atm were computed using vibrational analysis data. Global electron density transfer (GEDT) [56] was calculated according to the formula: where q A is the net Mulliken charge and the sum is taken over all the atoms of the dipolarophile. Indexes of r-bonds development (l) were calculated according to the formula [30]: where r A-B TS is the distance between the reaction centers A and B at the TS and r A-B P is the same distance at the corresponding product.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.