A DFT computational study of the molecular mechanism of [3 + 2] cycloaddition reactions between nitroethene and benzonitrile N-oxides

DFT calculations were performed to shed light on the molecular mechanism of [3 + 2] cycloadditions of simple conjugated nitroalkenes to benzonitrile N-oxides. In particular, it was found that these processes proceed by a one-step mechanism through asynchronous transition states. According to the latest terminology, they should be considered polar but not stepwise processes.

With this in mind, we have found that the molecular mechanism of [3 + 2] cycloadditions of nitrile N-oxides to conjugated nitroalkenes requires comprehensive theoretical study using DFT methods. It should be emphasized that the aforementioned reactions have not yet been investigated in this manner, even though such work is important from theoretical and practical points of view. Due to the strong electrophilicity (ω) of the conjugated nitroalkenes (ω > 1.5 eV [20][21][22]), it is in fact very likely that these reactions proceed through a zwitterionic intermediate (paths A and B in Scheme 2). At the same time, it is possible that zwitterionic structures with Bextendedĉ onformations are created along the competing paths of the cycloaddition reaction (paths C and D in Scheme 2). Such a path for addition has recently been analyzed for the reactions of nitroacetylene with a series of allenyl-type TACs [23]. All of the above options deserve detailed consideration.
Recently, we performed theoretical studies of a number of different cycloaddition processes involving conjugated nitroalkenes [24][25][26][27][28][29][30], and subsequently carried out experiments examining the reaction selectivity [27,28,[30][31][32], activation parameters [33], and kinetic effects of the solvent and substituents [29,34] for these processes. This comprehensive approach provided good insight into the key properties of the transition states involved. We concluded that DFT calculations should be capable of determining the molecular mechanism of these [3 + 2] cycloadditions. Therefore, in the work reported in the present paper, in order to allow general conclusions to be drawn, benzonitrile N-oxide as well as analogs of it containing substituents with different electronic properties were tested as model TACs. In particular, we decided to (a) probe the nature of the interaction between the addents in an elementary cycloaddition step and (b) run simulations of theoretically possible reaction paths of nitroethene with various nitrile N-oxides. As an extension of those studies, we also present theoretical studies of similar reactions involving the extremely electrophilic [13] 1,1dinitroethene molecule here.

Computational methods
Global reactivity indices (electronic chemical potential μ, chemical hardness η, global electrophilicity ω, global nucleophilicity N) were estimated according to the equations recommended by Parr [33] and Domingo [35][36][37]. In particular, the electronic chemical potentials and chemical hardnesses of the reactants studied here were evaluated in terms of the one-electron energies of the frontier molecular orbitals using the following equations [35,36]: The values of μ and η were then used to calculate ω according to the formula N can be expressed [37] as The local electrophilicity (ω k ) [38] concentrated on atom k was calculated by projecting the index ω onto any reaction center k in the molecule using the Parr function P + k : The local nucleophilicity (N k ) [38] concentrated on atom k was calculated using the global nucleophilicity N and the Parr function P − k according to the formula The global and local electronic properties of the reactants considered in this work are collected in Tables 1 and 2.
To localize the transition states (TSs), the hybrid B3LYP functional and the 6-31G(d) basis set included in the GAUSSIAN 09 software package were applied. We subsequently also performed analogous computation at the more advanced B3LYP/6-31 + G(d) and B3LYP/6-311G(d) levels of theory.
The structures corresponding to critical points on the PES along the reaction paths were localized in an analogous manner to that used in a previous analysis of [3 + 2] cycloadditions of (Z)-C,N-diphenylnitrone with 1,1dinitroethene [13]. First-order saddle points were localized using the QST2 and Berny procedures. The transition states were verified by diagonalizing the Hessian matrix and analyzing the intrinsic reaction coordinates (IRC). The polarizable continuum model (PCM) [39], in which the cavity is created as a series of overlapping spheres, was used to calculate solvent effects. Calculations of all critical structures were performed for a temperature T = 298 K and pressure p = 1 atm. The global electron density transfer (GEDT) [40] was calculated according to the formula where q A is the net Mulliken charge, and the sum is performed over all the atoms of nitroalkene.
Indices of σ-bond development (l) were calculated according to the formula [14] where r TS A-B is the distance between the reaction centers A and B at the TS and r P A-B is the corresponding distance at the product.
The kinetic parameters as well as important properties of the critical structures are displayed in Tables 3, 4, 5, and 6.

Analysis of interactions between addents
First, we decided to shed some light on the nature of the interactions between the addents. To do this, we used the electronic properties of the reactants, which were estimated using equations defined based on conceptual density functional theory [41]. A similar approach was successfully used recently to explain the paths followed by a number of different biomolecular processes (see for example [42][43][44][45]). In this theory, nitroethene is classified as a strong electrophile (ω > 1.5 eV) [36]. Its electrophilicity can be enhanced by introducing a second electron-withdrawing group (EWG) at position 1 of the nitrovinyl fragment. Accordingly, the presence of a Cl substituent increases the value of ω for nitroalkene to 2.88 eV, and that of the nitro group to 3.56 eV.
The electronic properties of N-oxides of aromatic nitriles vary widely. The global electrophilicity of benzonitrile N-oxide 2d is 1.46 eV, which makes it a moderate electrophile. However, gradually increasing the electron-donating power of a substituent at the 4-position on the N-oxide phenyl ring causes a gradual change in its e l e c t r o n i c p r o p e r t i e s . I n p a r t i c u la r, fo r t h e 4dimethylamino-substituted N-oxide 2a, ω is below 1 eV. This means that 2a will exhibit only marginally electrophilic properties, but will show strong nucleophilic power, as indicated by the value of N (>3.8 eV). Replacing the dimethylamino group with a strongly electronwithdrawing nitro group results in a dramatic change in the properties of the N-oxide. In particular, the 4-nitrosubstituted N-oxide 2h is characterized by strong electrophilic properties (ω > 3 eV)-stronger than nitroethene (!).
We then analyzed the local reactivity for different pairs of reactants. We found that the oxygen atom on the CNO fragment is a strongly nucleophilic reaction center in all of the Noxides. In turn, the strongest electrophilic reaction center is Table 1 Global and local electronic properties of nitroethene (1) and its selected 1-substituted derivatives (5-7)

R
Global properties Local properties always the β atom of the nitrovinyl fragment in the nitroalkenes. If we assume that these centers govern the reaction path, then the products of the cycloadditions should be 4nitroisoxazolines. However, this conclusion conflicts with the experimental data, because only the corresponding 5nitroisoxazolines are formed in this process. It should be noted at this point that only extremely electrophilic addents were considered for further detailed study of the mechanistic aspects of the corresponding cycloadditions. This allowed us to get a good idea of the mechanism associated with the [3 + 2] cycloaddition processes without having to calculate the full reaction paths for all such processes.

Kinetic aspects
DFT calculations show that the [3 + 2] cycloaddition of nitroethene (1) to N-oxide 2d in weakly polar toluene proceeds by a one-step mechanism (Fig. 1, Table 3). In the first step of the process, the enthalpy of the reacting system decreases due to the formation of a pre-reaction complex (MC). Table 2 Global and local electronic properties of benzonitrile N-oxide and its 4-R-substituted analogs (2a-h)

Substituent
Global properties Local properties Depending on the reaction path, this complex can be converted to regioisomeric isoxazoline 3d or 4d. Each of these conversions entails overcoming an activation barrier. The kinetically favored conversion is that associated with path B, leading to the formation of a 5-nitro-substituted adduct, which-as shown by experimental studies-does indeed form during this reaction [6,10]. So, contrary to expectations arising from regioselectivity, based on the analysis of the stationary states of the addends (see the preceding paragraph), a DFT-based exploration of the reaction pathways can accurately determine the regiochemistry of cycloaddition. However, all attempts to find a reaction path that leads to the adduct through an acyclic intermediate were unsuccessful. We also failed to obtain any stable structures of hypothetical zwitterions in an Bextendedĉ onformation. Similar studies were also performed for reactions involving representative N -substituted analogs of benzonitrile N-oxide. In the calculations of the reaction paths for the parent reaction system involving nitroethene   and unsubstituted benzonitrile N-oxide, several basis sets were applied (Table 3). It was found that calculations using the simple 6-31G(d) basis set gave practically identical (from a mechanistic point of view) results to calculations using higher levels of theory. In particular, in each instance, DFT calculations indicated a one-step reaction mechanism in which the favored path leads to a product with the nitro group at the 5-position on the heterocyclic ring. Thus, only the B3LYP/6-31G(d) level of theory was applied in subsequent investigations. We found that, in quantitative terms, the energy profiles for the cycloadditions 1 + 2a and 1 + 2h are the same as that for 2d + 1 cycloaddition, although they would be somewhat different quantitatively. In particular, introducing the (electron-donating) dimethylamino substituent at the 4-position of the phenyl ring of the N-oxide reduces the activation barrier, whereas the introduction of the (electron-withdrawing) nitro group increases the aforementioned barrier. However, reaction path B is always kinetically favored.
We then decided to investigate the effect of solvent polarity on the course of the reaction. We found that replacing toluene with a more polar medium (water) prompts an increase in the activation barrier along both regioisomeric paths.
As an extension of our studies, we performed a similar analysis of the energy profile of a similar cycloaddition involving 1,1-dinitroethene 7 and N-oxide 2d (Scheme 3), which has not yet been studied experimentally. Our previous work [13] indicated that nitroalkene 7 will react with C,Ndiphenylnitrone via a stepwise zwitterionic mechanism. As in the case of 1 + 2 cycloaddition, we considered both of the theoretically possible regioisomeric reaction routes (A and B; see Scheme 3).
The results indicated that, in toluene, the energy profile of 7 + 2d cycloaddition is qualitatively the same as that of 2d + 1 cycloaddition. However, due to the strongly electrophilic character of nitroalkene 7, these processes take place much more quickly than in the reaction involving nitroethene 1 ( Table 4). The regioselectivity of this cycloaddition is different from that of 1 + 2d cycloaddition. In particular, from a kinetic point of view, the path leading to the 4-nitro-substituted adduct is favored. All attempts to find paths leading to cycloadducts through an acyclic intermediate step were unsuccesful, as were efforts to find Fig. 1 Energetic profiles for cycloaddition between nitroethene (1) and benzonitrile N-oxide (2d) in toluene according to B3LYP/6-31G(d) (PCM) calculations paths leading to hypothetical zwitterions with an Bextended^conformation. Very similar results were seen for calculations of 7 + 2d cycloaddition in the presence of a more polar medium (water).
The studies described above indicate that there is an important difference between the paths of the reactions of model nitriles and imine N-oxides with nitroethene 1 and dinitroethene 7. The former TAC reacts with both of those nitroalkenes (despite their rather different global electrophilicities) by a one-step mechanism, while the latter reacts with nitroethene 1 according to a single-step mechanism but with 1,1-dinitroethene 7 by a two-step zwitterionic mechanism [13].
Key structures B3LYP/6-31G (d) calculations showed that two new σbonds always form in both TSs of the 1 + 2d reaction (Fig. 2). These bonds are formed between the atoms O1 and C5 and between C3 and C4, although the bond with the β carbon of the nitroalkene substructure forms faster. Thus, the structures of the TSs correlate well with the data obtained upon analyzing the local electrophilicity in nitroalkenes. The energetically favored TS along path B is even more asynchronous (see the Δl values in Table 5). Both of the TSs considered are polar, as demonstrated by their dipole moments and GEDT values ( Table 5).
The synchronicity of TS A can be controlled to some extent by altering the substituent on the phenyl ring of TAC. In particular, introducing an electron-donating group increases the asynchronicity of the TS, while adding an electron-withdrawing group reduces its asynchronicity. The type of substituent present does not appear to influence the asynchronicity of TS B .
Next, we analyzed the effect of the polarity of the reaction medium on the characteristics of the TSs. We found that replacing the weakly polar solvent toluene with a strongly polar medium (nitromethane or water) increases the asynchronicity of the TSs. At the same time, they increase in polarity. However, these changes are not significant enough to induce a stepwise zwitterionic mechanism.
A similar analysis was performed for the TSs in the reaction between dinitroethene 7 and TAC 2d. In accordance with our expectations, these structures were found to be far more asynchronous and strongly polar (Fig. 3, Table 6) than the corresponding TSs in the 1 + 2d reaction. Our attempts to optimize stable zwitterionic structures that could hypothetically form in the 7 + 2d reaction failed.

Conclusions
Regardless of the basis set applied, our B3LYP results clearly indicate that [3 + 2] cycloadditions of simple nitroalkenes to arylonitrile N-oxides proceed via a one-step mechanism. According to Domingo's terminology [43], this mechanism should be interpreted as polar. DFT calculations also showed that the favored reaction path leads to an adduct with a nitro group on C5, which agrees well with experimental observations. The transition-state synchronicity can be controlled to a certain degree by changing the polarity of the reaction medium and the nature of the substituent on the N-oxide phenyl Fig. 2 Views of the TSs that occur during the cycloaddition of nitroethene (1) to benzonitrile N-oxide (2d) in toluene, as derived via B3LYP/6-31G(d) (PCM) calculations Scheme 3 Fig. 3 Views of the TSs that occur during the cycloaddition between 1,1dinitroethene (7) and benzonitrile N-oxide (2d) in toluene, as derived via B3LYP/6-31G(d) (PCM) calculations ring. However, this is not enough to induce a switch to a twostep reaction mechanism with a zwitterionic intermediate. Nor was such a mechanism found to be possible in an analogous cycloaddition with a much stronger electrophilic component (1,1-dinitroethene).