On the Mechanism of the Synthesis of Nitrofunctionalised ∆ 2 -Pyrazolines via [3+2] Cycloaddition Reactions between α -EWG-Activated Nitroethenes and Nitrylimine TAC Systems On the Mechanism of the Synthesis of Nitrofunctionalised Δ 2 -Pyrazolines via [3+2] Cycloaddition Reactions between α -EWG-Activated Nitroethenes and Nitrylimine TAC Systems

: We investigated the reactivity of different substituted nitrylimine-type three atom components (TACs) in [3+2] cycloaddition (32CAs) reactions with electrophilically activated nitroethenes within molecular electron density theory (MEDT). In parallel research, the molecular mechanism of the considered transformation was examined through analysis of all possible reaction channels and full optimization of all critical structures. In particular, the existence of zwitterionic intermediates on reaction paths was veriﬁed. On the basis of the bonding evolution theory (BET), the mechanism of the 32CA reaction between C,N-diphenylnitrylimine and (E)-2-phenyl-1-cyano-1-nitroethene should be treated as a one-step two-stage mechanism. Abstract: We investigated the reactivity of different substituted nitrylimine-type three atom components (TACs) in [3+2] cycloaddition (32CAs) reactions with electrophilically activated nitroethenes within molecular electron density theory (MEDT). In parallel research, the molecular mechanism of the considered transformation was examined through analysis of all possible reaction channels and full optimization of all critical structures. In particular, the existence of zwitterionic intermediates on reaction paths was verified. On the basis of the bonding evolution theory (BET), the mechanism of the 32CA reaction between C,N-diphenylnitrylimine and (E)-2-phenyl-1-cyano-1-nitroethene should be treated as a one-step two-stage mechanism.


Introduction
This work is a continuation of our comprehensive, experimental, and theoretical study regarding the synthesis and reactivity of (E)-2-aryl-1-cyano-1-nitroethenes (ACN). Some examples of this group of conjugated nitroalkenes have been known since the second half of the 20th century [1,2]. In recent years, similar-type compounds have been synthetized [3,4]. Unfortunately, knowledge of their chemical properties is still limited. Recently, we detected a series of interesting chemical properties of ACNs regarding their participation in cycloaddition processes. For example, despite the high reactivity of the shielded reaction sites, ACNs react rapidly with cyclopentadiene even at r.t., yielding mixtures of respective endo-nitro and exo-nitro cycloadducts [5] (Scheme 1). In similar reactions involving the mixture of methylcyclopentadienes, of which many are possible, stereoisomeric products are formed [6]. In contrast, less sterically crowded 2-aryl-1-nitroethenes react with cyclopentadiene at temperatures up to 80 • C [7].

Introduction
This work is a continuation of our comprehensive, experimental, and theoretical study regarding the synthesis and reactivity of (E)-2-aryl-1-cyano-1-nitroethenes (ACN). Some examples of this group of conjugated nitroalkenes have been known since the second half of the 20th century [1,2]. In recent years, similar-type compounds have been synthetized [3,4]. Unfortunately, knowledge of their chemical properties is still limited. Recently, we detected a series of interesting chemical properties of ACNs regarding their participation in cycloaddition processes. For example, despite the high reactivity of the shielded reaction sites, ACNs react rapidly with cyclopentadiene even at r.t., yielding mixtures of respective endo-nitro and exo-nitro cycloadducts [5] (Scheme 1). In similar reactions involving the mixture of methylcyclopentadienes, of which many are possible, stereoisomeric products are formed [6]. In contrast, less sterically crowded 2-aryl-1nitroethenes react with cyclopentadiene at temperatures up to 80 °C [7].
[3+2] cycloaddition (32CA) reactions of ACNs and N-methylazomethine ylide proceed under mild conditions and produce nitropyrrolidines [8]. In contrast, in analogous 32CAs involving diazafluorene, zwitterionic acyclic adducts were formed in the first reaction step. These zwitterions converted spontaneously under the reaction conditions into azine molecular systems [9]. At the same time, other nitroalkenes characterized by similar electrophilicity as ACNs reacted with diazafluorene to produce ∆ 1 -nitropyrazolines (Scheme 2).

Scheme 2. The zwitterionic intermediate in reaction of diazafluorene and ACNs.
Non-catalyzed 32CA reactions between ACNs and nitrile N-oxides proceed unexpectedly with the participation of a nitrile bond, instead of resulting in the expected >C=C< moiety of nitroalkene [10] (Scheme 3). There are evidently rare cases of 32CAs in the CN bond, which is generally recognized as inactive in the 32CA processes. There is currently no published work regarding the participation of ACNs in 32CA processes involving nitrylimine TAC systems. Our work initiates comprehensive research in this area. We aim to shed light on the selectivity and molecular mechanism of model processes involving different substituted diarylnitrylimines (DNI) and ACNs. This analysis will be help to guide further experimental study in the presented area. Non-catalyzed 32CA reactions between ACNs and nitrile N-oxides proceed unexpectedly with the participation of a nitrile bond, instead of resulting in the expected >C=C< moiety of nitroalkene [10] (Scheme 3). There are evidently rare cases of 32CAs in the CN bond, which is generally recognized as inactive in the 32CA processes.

Scheme 2. The zwitterionic intermediate in reaction of diazafluorene and ACNs.
Non-catalyzed 32CA reactions between ACNs and nitrile N-oxides proceed unexpectedly with the participation of a nitrile bond, instead of resulting in the expected >C=C< moiety of nitroalkene [10] (Scheme 3). There are evidently rare cases of 32CAs in the CN bond, which is generally recognized as inactive in the 32CA processes. There is currently no published work regarding the participation of ACNs in 32CA processes involving nitrylimine TAC systems. Our work initiates comprehensive research in this area. We aim to shed light on the selectivity and molecular mechanism of model processes involving different substituted diarylnitrylimines (DNI) and ACNs. This analysis will be help to guide further experimental study in the presented area. There is currently no published work regarding the participation of ACNs in 32CA processes involving nitrylimine TAC systems. Our work initiates comprehensive research in this area. We aim to shed light on the selectivity and molecular mechanism of model processes involving different substituted diarylnitrylimines (DNI) and ACNs. This analysis will be help to guide further experimental study in the presented area.

Computational Details
The global reactivity descriptors of the addends, namely electronic potential µ, chemical hardness η, global electrophilicity ω and global nucleophilicity N, were approximated in pursuance of the equations defined on the basis of conceptual density functional theory (CDFT) according to the equations recommended by Parr [11] and Domingo [12,13]. In the calculation we used the correlation-exchange functional B3LYP together with the basic level set of 6-31G(d) in the gas phase [12][13][14][15].
The electronic chemical potentials (µ) and chemical hardness (η) were evaluated in terms of one-electron energies of FMO (E HOMO and E LUMO ) using the following equations [11][12][13][14][15][16]: where E HOMO and E LUMO may be approached in terms of the one-electron energies of the frontier MOs respectively HOMO and LUMO. Next, values of µ and η were then used to calculate a global electrophilicity (ω) according to the formula [12,13]: The global nucleophilicity (N) can be presented as follow [13]: where E HOMO (TCE) is the HOMO energy for tetracyanoethylene (TCE); is the reference, because it presents the lowest HOMO (E HOMO (TCE) = −9.368 eV). The local electrophilicity (ω k ) and the local nucleophilicity (N k ) concentrated on atom k was calculated based on global properties and the Parr function (P k + or P k − ), according to the formulas [17]: For localization of the transition states (TSs) the wb97xd/6-311+g(d) level of theory was applied [18]. All transition states were verified by diagonalization of the Hessian matrix and by analysis of the intrinsic reaction coordinates (IRC). For the simulation of solvents effect, the polarizable continuum model (PCM) [19] was used. Calculations of all critical structures were performed for the temperature T = 298 K and pressure p = 1 atm. Global electron density transfer (GEDT) [20] was calculated according to the formula: where q A is the net Mulliken charge and the sum takes over all the atoms of nitroalkene. Indexes of σ-bonds development (I) were calculated according to formula [21]: 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 same distance at the corresponding product.
Organics 2022, 3, FOR PEER REVIEW 4 pyrazolines 4 can be formed in the consequence of interactions between addends. Firstly, we decide to lead an exploration of the nature of these interactions.

Global Reactivity
The Conceptual Density Functional Theory (CDFT) should be considered as a strong implement that helps comprehension the reactivity of components in polar processes such as cycloaddition. All of indices was determined via calculations based on B3LYP/6-31G(d) theory level in the gas phase. It is very useful for estimate electrophilicity and nucleophilicity of the reagents [35][36][37]. Due to this, such indices as global reactivity indices: electronic chemical potential μ, chemical hardness η, global electrophilicity ω, global nucleophilicity N were determined and given in Table 1.  The Conceptual Density Functional Theory (CDFT) should be considered as a strong implement that helps comprehension the reactivity of components in polar processes such as cycloaddition. All of indices was determined via calculations based on B3LYP/6-31G(d) theory level in the gas phase. It is very useful for estimate electrophilicity and nucleophilicity of the reagents [35][36][37]. Due to this, such indices as global reactivity indices: electronic chemical potential µ, chemical hardness η, global electrophilicity ω, global nucleophilicity N were determined and given in Table 1. Analysis of the electronic chemical potential µ, can define the direction of the electron density flux between reagents in determined path of reaction. In case of the model reaction between nitrylimine 1a and nitroethene 2a, the electronic chemical potential µ [35,36] of 1a, −3.37 eV, is significantly higher than this of 2a, −5.28 eV, (Table 1). That means, in a course of model reaction 1a with 2a, the electron density flux will take place from nitrylimine 1a to nitroethene 2a. Similar results can be observed for the other reagents, as long as all the nitrylimine components 1b-e can be characterized by substantially higher value of electronic chemical potential than nitroethanes 2b-c (Table 1).
Calculated index of electrophilicity [35,36] ω of C,N-diphenylnitrylimine (1a) is 1.50 eV. In turn, the calculated index of nucleophilicity [13] N for this compound is 3.86 eV ( Table 1). These values give the conclusion that model nitrylimine 1a acts similar to moderate electrophile and strong nucleophile in a polar reaction.
In turn, calculated phenyl index of nucleophilicity N for this compound is 2.87 eV ( Table 1). These values provide the conclusion that model nitroalkene 2a acts similar to strong electrophile and also moderate nucleophile in a polar reaction.
Very recently, organic reactions have been classified as forward electron density flux (FEDF) and reverse electron density flux (REDF) reactions, depending on the direction of the flux of the electron density [38,40]. Non-polar reactions are classified as null electron density flux (NEDF) reactions [38,41]. Thus, the reactions involving nucleophilic nitrylimines 1a, 1b and 1d and electrophilic nitroalkenes 2a-c are classified as FEDF in agreement with the CDFT analysis. On the other hand, the non-polar 32CA reactions of 1c or 1e with all of the nitrylimine components 1a-e should be classified as NEDF.

Local Reactivity
The regioselectivity of polar processes of non-symmetric reagents can be specified through interaction between the most electrophilic center of the electrophile and the most nucleophilic ones of the nucleophile. For this purpose, electrophilic P k + and nucleophilic P k − Parr functions, derived from the changes of spin electron density reached via the GEDT process from the nucleophile to the electrophile, can be used as a powerful tool in the study of the local reactivity [42][43][44]. According to the nucleophilic P k − Parr functions of nitrylimine components 1a,c,e and the electrophilic P k + Parr functions of nitroethanes 2a-c in order to characterize the most nucleophilic and electrophilic centers of the species involved in this polar 32CA reaction were analyzed ( Figure 1). 1). It follows that, the 32CA reaction of nitrylimines 1c and 1e with nitroethenes 2a-c will not have a polar character [13,16].
Very recently, organic reactions have been classified as forward electron density flux (FEDF) and reverse electron density flux (REDF) reactions, depending on the direction of the flux of the electron density [38,40]. Non-polar reactions are classified as null electron density flux (NEDF) reactions [38,41]. Thus, the reactions involving nucleophilic nitrylimines 1a, 1b and 1d and electrophilic nitroalkenes 2a-c are classified as FEDF in agreement with the CDFT analysis. On the other hand, the non-polar 32CA reactions of 1c or 1e with all of the nitrylimine components 1a-e should be classified as NEDF.

Local Reactivity
The regioselectivity of polar processes of non-symmetric reagents can be specified through interaction between the most electrophilic center of the electrophile and the most nucleophilic ones of the nucleophile. For this purpose, electrophilic Pk + and nucleophilic Pk − Parr functions, derived from the changes of spin electron density reached via the GEDT process from the nucleophile to the electrophile, can be used as a powerful tool in the study of the local reactivity [42][43][44]. According to the nucleophilic Pk -Parr functions of nitrylimine components 1a,c,e and the electrophilic Pk + Parr functions of nitroethanes 2ac in order to characterize the most nucleophilic and electrophilic centers of the species involved in this polar 32CA reaction were analyzed ( Figure 1). Analysis of the electrophilic Pk + Parr functions of (E)-2-phenyl-1-cyano-1-nitroethenes (2a) indicates that the most electrophilic center is situated on carbon atom α, Pk + = 1.41 eV (Figure 1). The presence of amino or nitro group at para position of the phenyl ring of nitroethenes 2b and 2c does not cause significant changes in local reactivity. It means that for nitroethenes 2b and 2c the most electrophilic center is also located in atom α of carbon, Analysis of the electrophilic P k + Parr functions of (E)-2-phenyl-1-cyano-1-nitroethenes (2a) indicates that the most electrophilic center is situated on carbon atom α, P k + = 1.41 eV (Figure 1). The presence of amino or nitro group at para position of the phenyl ring of nitroethenes 2b and 2c does not cause significant changes in local reactivity. It means that for nitroethenes 2b and 2c the most electrophilic center is also located in atom α of carbon, P k + = 1.44 eV (2b) and 1.27 eV (2c), respectively ( Figure 1). It implies that center will react with the most nucleophilic α center of nitrylimines 1a,c,e.
On the other hand, the analysis of the nucleophilic P k − Parr functions of the C,Ndiphenylnitrylimine (1a) indicates that the carbon atom of -N=N=C-fragment constitutes the most nucleophilic center at molecule, presenting the maximum values at P k − = 0.32 eV (Figure 1). Same as in the previous example the presence of amino group at para position of the benzene ring for nitroethenes 2c and 2e does not significantly change the values of local reactivities. Therefore, nitrylimines 2c and 2e show a similar reactivity when compared with C,N-diphenylnitrylimine (1a) and presenting the maximum values at P k − = 0.33 (2c) and 0.36 (2e) eV, respectively ( Figure 1).
Summarizing, according to CDFT theory [12,14] the polar reactions of nitrylimines 1a,c,e with nitroalkenes 2a-c will be realized through the interaction of α atom of carbon for 1a,c,e with the carbon atom of -N=N=C-fragment for 2a-c. Therefore, the more preferred reaction channel for analyzed [3+2] cycloaddition reactions are forming 5-nitrosubstituted ∆ 2 -isoxazolines 4a,c,e according to path B. The regioselectivity of 32CA between nitrylimines 1a,c,e and nitroalkenes 2a-c cannot be determined using similar approach due to the non-polar nature of these processes [45].

Reaction Profiles
The mechanistic study regarding the title reactions we initiated from the exploration of the model process involving parent nitrylimine 1a and parent CNA 2a. It was found, that in toluene solution, the nature of energy profiles of both considered reaction channels are qualitatively similar (Figure 2). In particular, between valleys of individual reagents and products, two critical points were localized. These are connected with the existence of pre-reaction molecular complexes (MC) as well as transition states (TS). The first transformation of the reaction system is a formation of MC. This stage is characterized by reduction of the enthalpy of the reaction system about 11-12 kcal/mol (Table 2). However, the entropy factors do not determine the possibility of the existence of MCs as stable intermediates (∆G > 0 kcal/mol). Within MCs, addends are oriented relatively, for the achieving maximally good coulombic interactions between substructures (Table 3, Figure 3). These are not a charge-transfer complexes, which was confirmed by GEDT analysis [20]. Within MCs, reaction centers are oriented for the fashion, which determine further, positive interactions between substructures. Key interatomic distances (Tables 4 and 5) are characterized by values which are beyond of the typical range for formed sigma-bonds in transition states. The further conversion of MCs, lead on both considered paths to area of the existence of TSs. This requires an overcome of the energetical barrier about 8kcal/mol. It should be mentioned, that both considered paths should be permitted, as allowed from the kinetic point of view. Next, small difference between energies of the activation suggest rather low reaction regioselectivity. Within TSs, great amount of the charge transfer between substructures are observed (0.67e and 0.39e for TSA and TSB, respectively). This confirms expected previously, polar nature of considered cycloadditions. Next, in the framework of TS, the key interatomic distances are reduced substantially, in the comparison for analogous distances within MCs. The formation of new sigma-bonds proceeds evidently according to the asynchronous manner (Tables 4 and 5). This is typical for polar 32Cas processes involving unsymmetrical substituted components. The asynchronicity of TSs are however not sufficient for the extort the stepwise mechanism of the cycloaddition. The IRC analysis connects without any doubts both TSs, with respective MCs and respective products.
In the second part of our research, we decided to perform a similar study regarding to the 1a+2a process for the simulated presence of polar solvent-nitromethane. It was found, that both, the nature of energy profiles, as well as the quantitative description of all critical points, are close to results obtained for the simulated presence of toluene. Parallel, we detected two alternative channels of the addition reactions between considered reagents, which are not available in the toluene solution (paths C and D) (Scheme 5). In both cases, the first stage of the reaction is-similarly as in the case of 32CAs process-the formation of pre-reaction MC complexes. Their nature is close, to observed for cycloaddition reactions, but their further transformations proceed via different mechanistic scheme. In particular, conversion of mentioned MCs, leads directly to transition states TSC and TCD, respectively. Within these TSs, the distance between reaction centers C5-N1 are reduced to values, which can be considered as the stage of the formation of new sigma-bond. At the same time, the distances C3-C4 are characterized by values which are beyond of the typical range for formed sigma-bonds in transition states. Both localized TSs exhibit evidently polar nature, which is confirmed by analysis of GEDT values (Tables 4 and 5). The further conversion of TSC and TSD leads directly for acyclic adducts Z1 and Z2, respectively. This was confirmed by IRC analysis. Optimized adducts exhibit zwitterionic nature, which is evidently right due to GEDT values (GEDT > 1e). The direct transformation of these zwitterions into cycloadduct is however impossible, as the consequence of their conformation. The theoretically possible transformations from Z1 or Z2 to pyrazoline systems must proceed via dissociation in to individual reagents, and, in next step, via 32CA processes described above. Only, transformation between Z1 and Z2 zwitterions are possible, according to the simple rotation around C5-N1 bond.
Lastly, we examined the influence of the substituents' nature into reaction course. It was found that, the molecular reaction mechanism in all cases is close to observed for 32CAs involving diarylnitrylimine 1a and ACN 2a. In particular, in any case the stepwise zwitterionic mechanism was not detected. Only, the qualitative description of energy profile as well as key structures were changing in some range.    In the second part of our research, we decided to perform a similar study regarding to the 1a+2a process for the simulated presence of polar solvent-nitromethane. It was found, that both, the nature of energy profiles, as well as the quantitative description of all critical points, are close to results obtained for the simulated presence of toluene. Parallel, we detected two alternative channels of the addition reactions between considered reagents, which are not available in the toluene solution (paths C and D) (Scheme 5). In both cases, the first stage of the reaction is-similarly as in the case of 32CAs process-the   of their conformation. The theoretically possible transformations from Z1 or Z2 to pyra zoline systems must proceed via dissociation in to individual reagents, and, in next step via 32CA processes described above. Only, transformation between Z1 and Z2 zwitterion are possible, according to the simple rotation around C5-N1 bond. Lastly, we examined the influence of the substituents' nature into reaction course. I was found that, the molecular reaction mechanism in all cases is close to observed fo 32CAs involving diarylnitrylimine 1a and ACN 2a. In particular, in any case the stepwis zwitterionic mechanism was not detected. Only, the qualitative description of energy pro file as well as key structures were changing in some range.

BET Analysis of the 32CA between C,N-diphenylnitrylimine 1a and (E)-2-phenyl-1-cyano-1-nitroethene 2a
In order to carefully investigate the bonding changes occurring in the 32CA of the 1 and 2a, we decided to conduct a BET analysis. Scheme 6 represents the molecular mecha nism by Lewis-like structures resulted from the ELF topological analysis. The most signif icant ELF basin populations together with attractor positions for a reaction leading to 4a are gathered in Table 6 and Figure 4.  In order to carefully investigate the bonding changes occurring in the 32CA of the 1a and 2a, we decided to conduct a BET analysis. Scheme 6 represents the molecular mechanism by Lewis-like structures resulted from the ELF topological analysis. The most significant ELF basin populations together with attractor positions for a reaction leading to 4aa are gathered in Table 6 and Figure 4.

BET Study of the Creation of Acyclic Adduct Z1
We also decided to analyze the bonding changes along the formation acyclic addu Z1 by doing the BET analysis. The most significant ELF basin populations and attrac positions for structures participating in the studied reaction are collected in Table 7 a Figure 5. The proposed molecular mechanism of the 32CA reaction between 1a and 2a. The bonding changes along this 32CA between C,N-diphenylnitrylimine 1a and 2a are characterized by nine different topological phases ( Table 4). The ELF picture of the first structure MCB is very similar to those of the divided reagents. Phase II begins at P1 1 , where it can be observed that the V(C4,C5) and V'(C4,C5) disynaptic basins current in previous MCB point, have combined into one new V(C4,C5) disynaptic basin with initial population 3.45 e. Along Phase III, which begins at P2 1 , the new V(N2) monosynaptic basin is created with population of 0.94 e. Creation of this monosynaptic basin is related to depopulation of V(N1,N2), V(N2,C3) and V'(N2,C3) disynaptic basins and the growth of the population of V(C3) monosynaptic basin. At this point, we found the TS structure of the 32CA between 1a and 2a (TSB, d(C3-C4) = 2.176 Å and d(C5-N1) = 2.646 Å). Phase IV initiates at the structure P3 1 , where can be noticed the creation of a V(C4) monosynaptic basin with 0.16 e, which is connected with depopulation of V(C4,C5) disynaptic basin. Creation of this monosynaptic basin is associated with formation a pseudoradical [46] center at C4 carbon atom (Figure 4).
Along Phase V, which begins at P4 1 , we noticed that pseudoradical centers located on C3 and C4 atoms combined into C3-C4 bonding area with an starting value of 1.41 e. The electron density changes mean that creation of the C3-C4 bond starts with a distance d(C3,C4) = 2.014 Å, by divvying the non-bonding electron densities of the two C3 and C4 centers. Phase VI, starts at P5 1 , and is related to formation a new pseudoradical center at C5 carbon atom with value of 0.41 e (Figure 4 and Table 6). Formation of a new V(C5) monosynaptic basin causes depopulation of V(C4,C5) disynaptic basin. In Phase VII, the V(N2,C3) and V'(N2,C3) disynaptic basins, have merged and a new V(N2,C3) disynaptic basin was created with 3.17 e. This shift is related to disruption of the N2-C3 bond and formation a partial double bond (Scheme 6). At P7 1 begins Phase VIII, which is associated with creation a V'(N2) monosynaptic basin with starting value of 0.48 e. This shift is caused by division of V(N1) monosynaptic basin, present in P6 1 , for two new V(N1) and V'(N1) monosynaptic basins and depopulation of V(N1,N2) disynaptic basin. Subsequently we observed in the last Phase IX, which is located between P8 1 and 4aa (Table 6). In this phase, the creation of the second C5-N1 bond follows through the connection of the non-bonding electron densities of the two V(C5) and V(N1) monosynaptic basins. The C5-N1 bond was formed with starting distance of d(C5,N1) = 1.750 Å.
Organics 2022, 3, FOR PEER REVIEW 14 On the grounds of the BET study, we gather that: (i) the 32CA reaction between C,Ndiphenylnitrylimine 1a and (E)-2-phenyl-1-cyano-1-nitroethene 2a can be depicted by nine topologically various phases; (ii) the activation energy of this reaction, 4.0 kcal·mol −1 , is related to the formation of C4 pseudoradical center and lone pair at N2 nitrogen atom; (iii) creation of first C3-C4 bond follows in Phase V by way of merging two C3 and C pseudoradical centers; (iv) creation of the C5-N1 bond occurs in the last Phase IX by merging the C5 pseudoradical center and N1 non-bonding lone pair; (v) it is worth nothing that when the C5-N1 bond begins to form, the first bond is fully formed. According to that, we can conclude that the 32CA of the 1a and 2a proceeds according to one-step two-stage mechanism.  On the grounds of the BET study, we gather that: (i) the 32CA reaction between C,Ndiphenylnitrylimine 1a and (E)-2-phenyl-1-cyano-1-nitroethene 2a can be depicted by nine topologically various phases; (ii) the activation energy of this reaction, 4.0 kcal·mol −1 , is related to the formation of C4 pseudoradical center and lone pair at N2 nitrogen atom; (iii) creation of first C3-C4 bond follows in Phase V by way of merging two C3 and C pseudoradical centers; (iv) creation of the C5-N1 bond occurs in the last Phase IX by merging the C5 pseudoradical center and N1 non-bonding lone pair; (v) it is worth nothing that when the C5-N1 bond begins to form, the first bond is fully formed. According to that, we can conclude that the 32CA of the 1a and 2a proceeds according to one-step two-stage mechanism.

BET Study of the Creation of Acyclic Adduct Z1
We also decided to analyze the bonding changes along the formation acyclic adduct Z1 by doing the BET analysis. The most significant ELF basin populations and attractor positions for structures participating in the studied reaction are collected in Table 7 and Figure 5. Scheme 6. The proposed molecular mechanism of the 32CA reaction between 1a and 2a.

BET Study of the Creation of Acyclic Adduct Z1
We also decided to analyze the bonding changes along the formation acyclic adduct Z1 by doing the BET analysis. The most significant ELF basin populations and attractor positions for structures participating in the studied reaction are collected in Table 7 and Figure 5. The bonding changes along the formation of acyclic adduct Z1 are described by eight various phases ( Table 7). The first point of the IRC is the molecular complex MCC, which presents ELF picture of separated reagents. We observed that the V(N1) monosynaptic and V(N2,C3), V'(N2,C3) and V(C4,C5) disynaptic basins show a slight larger ELF valence basin populations than those shown in the substrates. In turn, in the event of V(N1,N2) and V'(C4,C5) disynaptic and V(C3) monosynaptic basins, we observed the lower values of the ELF basin populations. Phase II begins at structure P1 2 . In this area, the two V(C4,C5) and V'(C4,C5) disynaptic basins merged into one new V(C4,C5) disynaptic basin with starting population of 3.49 e. This topological shift is connected with break of C4-C5 bond in (E)-2phenyl-1-cyano-1-nitroethene 2a. Along Phase III, the C3 monosynaptic basin disappear and we observed the increase in the population of V(N2,C3) disynaptic basin to about 3.81 e. Along Phase IV, the V(N2,C3) and the V'(N2,C3) disynaptic basins merged into one original V(N2,C3) disynaptic basin integrating 5.87 e. At P3 2 , we notice the TS of the formation of adduct Z1 (TSC, d(C4-N1) = 2.062 Å). Phase V begins at P4 2 and is related to create a V(C5) monosynaptic basin integrating 0.56 e. This is connected with a creation of pseudoradical center at C5 carbon atom and depopulation of V(C4,C5) disynaptic basin. In Phase VI, we observed the formation a new V(C4,N1) disynaptic basin with population of 1.31 e, created by depopulation of V(N1) monosynaptic and V(C4,C5) disynaptic basins. This topological change is associated with formation of the C4-N1 single bond ( Figure 5). The next Phase VII begins at P6 2 and is associated with creation of a second V'(C5) monosynaptic basin with population of 0.37 e. Forming of pseudoradical center at C5 carbon atom is associated with decrease the population of V(C4,C5) disynaptic basin. The final Phase VIII, placed between P7 2 and Z1 is related to formation a V(N2) monosynaptic basin through the depopulation of V(N2,C3) disynaptic basin. starting population of 3.49 e. This topological shift is connected with break of C4-C in (E)-2-phenyl-1-cyano-1-nitroethene 2a. Along Phase III, the C3 monosynaptic ba appear and we observed the increase in the population of V(N2,C3) disynaptic about 3.81 e. Along Phase IV, the V(N2,C3) and the V'(N2,C3) disynaptic basins into one original V(N2,C3) disynaptic basin integrating 5.87 e. At P32, we notice th the formation of adduct Z1 (TSC, d(C4-N1) = 2.062 Å ). Phase V begins at P42 and is to create a V(C5) monosynaptic basin integrating 0.56 e. This is connected with a of pseudoradical center at C5 carbon atom and depopulation of V(C4,C5) disynapti In Phase VI, we observed the formation a new V(C4,N1) disynaptic basin with pop of 1.31 e, created by depopulation of V(N1) monosynaptic and V(C4,C5) disynaptic This topological change is associated with formation of the C4-N1 single bond (Fi The next Phase VII begins at P62 and is associated with creation of a second V'(C osynaptic basin with population of 0.37 e. Forming of pseudoradical center at C5 atom is associated with decrease the population of V(C4,C5) disynaptic basin. T Phase VIII, placed between P72 and Z1 is related to formation a V(N2) monosynapt through the depopulation of V(N2,C3) disynaptic basin. On the basis of BET study, we may say that: (i) the molecular mechanism of mation of acyclic adduct Z1 can be featured by eight various phases; (ii) the ac On the basis of BET study, we may say that: (i) the molecular mechanism of the formation of acyclic adduct Z1 can be featured by eight various phases; (ii) the activation energy of this reaction, 8.1 kcal·mol −1 , is mostly related to breaking of the N2-C3 double bond and formation of C5 pseudoradical center; (iii) formation of C4-N1 single bond takes place in Phase VI by the depopulation of V(N1) monosynaptic and V(C4,C5) disynaptic basins.

Conclusions
The DFT computational study shows, that despite of high electrophilic nature of nitroalkenes, the [3+2] cycloaddition reactions between diarylnitrylimines and 2-aryl-1cyano-1-nitroethenes proceed via single-step mechanism. In a polar nitromethane solution, the cycloaddition process can however compete, with the formation of zwitterionic structures characterized by "extended" conformations. These intermediates cannot be cyclized in to heterocyclic systems via simple one-step reactions, because the key reaction sites are localized on too long relative distance. Their conversion into pyrazoline systems is possible only via the stage of dissociation into individual reagents, and subsequently, one-step cycloaddition.
The topological analysis of the bonding changes associated with 32CA reaction between C,N-diphenylnitrylimine 1a and (E)-2-phenyl-1-cyano-1-nitroethene 2a, can be described by nine topologically different phases. Formation of the first C-C single bond follows by merging of two pseudoradical centers. In turn, the second C-N bond is formed in the last phase of the reaction when the first bond is fully formed. According to that, we can conclude that the 32CA of the 1a and 2a progress according to one-step two-stage mechanism.