On the Question of Zwitterionic Intermediates in the [3+2] Cycloaddition Reactions between Aryl Azides and Ethyl Propiolate

The molecular mechanism of the [3+2] cycloaddition reactions between aryl azides and ethyl propiolate was evaluated in the framework of the Molecular Electron Density Theory. It was found that independently of the nature of the substituent within the azide molecule, the cycloaddition process is realized via a polar but single-step mechanism. All attempts of localization as postulated earlier by Abu-Orabi and coworkers’ zwitterionic intermediates were not successful. At the same time, the formation of zwitterions with an “extended” conformation is possible on parallel reaction paths. The ELF analysis shows that the studied cycloaddition reaction leading to the 1,4-triazole proceeds by a two-stage one-step mechanism. It also revealed that both zwitterions are created by the donation of the nitrogen atom’s nonbonding electron densities to carbon atoms of ethyl propiolate.

The most universal methodology of the preparation of five-membered heterocycles is [3+2] cycloaddition reactions (32CA reactions) with the participation of different types of three-atom components (TACSs) [18][19][20].The preparation of the 1,2,3-triazole molecular segment is also possible in this way due to 32CA between azides and alkynes.
Full regioselective synthesis of 1,2,3-triazoles via the 32CA reaction and under milder conditions is possible via the catalytic protocol, including the presence of cooper or rhutenium catalysts.For example, benzyl azide reacts with simple alkynes at 60 • C. The reaction is realized with full regioselectivity (Scheme 3) [33].
Molecules 2023, 28, x FOR PEER REVIEW 2 of 24 Scheme 1.The regioselectivity of the 32CA reactions between benzyl azide and phenoxyethyne.
Relatively easier are analogous 32CAs with the participation of cyclooctyne analogs [30,31].This is, however, an incidental example, determined by the specific, strained nature of the alkyne.Similar, atypical reactivity is also observed in the case of azidefullerene C20 cycloadditions (Scheme 2) [32].
Scheme 2. The 32CA reactions between azides and strained alkynes.
Full regioselective synthesis of 1,2,3-triazoles via the 32CA reaction and under milder conditions is possible via the catalytic protocol, including the presence of cooper or rhutenium catalysts.For example, benzyl azide reacts with simple alkynes at 60 °C.The reaction is realized with full regioselectivity (Scheme 3) [33].
Scheme 3. The rhutenium-catalyzed 32CA reactions between benzyl azide and simple alkynes.Some years ago, Abu-Orabi and coworkers [34] described an interesting case of the 32CA reactions between aryl azides (1a-e) and ethyl propiolate (2) (Scheme 4).Independently of the absence of a catalyst, all analyzed reactions were realized under relatively mild conditions and (most importantly) with full regioselectivity.The authors explained this phenomenon by assuming a stepwise mechanism via a zwitterionic intermediate.The presence of the zwitterionic intermediates is, however, not confirmed in any way.Scheme 4. The regioselectivity and reaction conditions of the 32CAs between aryl azides 1a-e and ethyl propiolate 2. Scheme 3. The rhutenium-catalyzed 32CA reactions between benzyl azide and simple alkynes.Some years ago, Abu-Orabi and coworkers [34] described an interesting case of the 32CA reactions between aryl azides (1a-e) and ethyl propiolate (2) (Scheme 4).Independently of the absence of a catalyst, all analyzed reactions were realized under relatively mild conditions and (most importantly) with full regioselectivity.The authors explained this phenomenon by assuming a stepwise mechanism via a zwitterionic intermediate.The presence of the zwitterionic intermediates is, however, not confirmed in any way.
Molecules 2023, 28, x FOR PEER REVIEW 2 of 24 Scheme 1.The regioselectivity of the 32CA reactions between benzyl azide and phenoxyethyne.
Relatively easier are analogous 32CAs with the participation of cyclooctyne analogs [30,31].This is, however, an incidental example, determined by the specific, strained nature of the alkyne.Similar, atypical reactivity is also observed in the case of azidefullerene C20 cycloadditions (Scheme 2) [32].
Scheme 2. The 32CA reactions between azides and strained alkynes.
Full regioselective synthesis of 1,2,3-triazoles via the 32CA reaction and under milder conditions is possible via the catalytic protocol, including the presence of cooper or rhutenium catalysts.For example, benzyl azide reacts with simple alkynes at 60 °C.The reaction is realized with full regioselectivity (Scheme 3) [33].
Scheme 3. The rhutenium-catalyzed 32CA reactions between benzyl azide and simple alkynes.Some years ago, Abu-Orabi and coworkers [34] described an interesting case of the 32CA reactions between aryl azides (1a-e) and ethyl propiolate (2) (Scheme 4).Independently of the absence of a catalyst, all analyzed reactions were realized under relatively mild conditions and (most importantly) with full regioselectivity.The authors explained this phenomenon by assuming a stepwise mechanism via a zwitterionic intermediate.The presence of the zwitterionic intermediates is, however, not confirmed in any way.It should be underlined that, at this moment, the actual state of knowledge shed new light on the earlier interpretation of the 32CA reaction mechanism as one-step "concerted" independently of the nature of starting molecules.In particular, at this moment, several different mechanisms can be assumed within the exploration of the 32CA mechanisms; in particular, (i) polar mechanisms (the one-step mechanism, two-stage one-step mechanism, and stepwise zwitterionic mechanism) [35,36]) or (ii) non-polar mechanism (the one-step mechanism, two-stage one-step mechanism asynchronous mechanism, and stepwise biradical mechanism) [37,38].Which of them is implemented in practice is determined by the interplay of two independent factors-the nature of electronic interactions [39] and the power of steric interactions [40].Next, Domingo generally undermined the term "pericyclic" regarding all 32CA reactions based on the recent progress of the exploration of reaction profiles via BET and ELF techniques [41,42].Therefore, the mechanistic aspects of the title processes evidently require deeper exploration.
Molecules 2023, 28, x FOR PEER REVIEW 3 of 24 It should be underlined that, at this moment, the actual state of knowledge shed new light on the earlier interpretation of the 32CA reaction mechanism as one-step "concerted" independently of the nature of starting molecules.In particular, at this moment, several different mechanisms can be assumed within the exploration of the 32CA mechanisms; in particular, (i) polar mechanisms (the one-step mechanism, two-stage one-step mechanism, and stepwise zwitterionic mechanism) [35,36]) or (ii) non-polar mechanism (the one-step mechanism, two-stage one-step mechanism asynchronous mechanism, and stepwise biradical mechanism) [37,38].Which of them is implemented in practice is determined by the interplay of two independent factors-the nature of electronic interactions [39] and the power of steric interactions [40].Next, Domingo generally undermined the term "pericyclic" regarding all 32CA reactions based on the recent progress of the exploration of reaction profiles via BET and ELF techniques [41,42].Therefore, the mechanistic aspects of the title processes evidently require deeper exploration.
In the initial step, we decided to analyze the nature of the intermolecular interactions between respective pairs of reagents (Table 1, Figures 1 and 2).It was found that ethyl propiolate is characterized by global electrophilicity equal to 1.48 eV.Therefore, according to the unique Domingo electrophilicity scale [43,44], this component should be treated as a moderate electrophile.The analysis of the distribution of local electrophilicities subsequently shows that the most activated reaction center is located on the terminal carbon atom of acetylenic moiety (ω 1 = 0.67 eV).In comparison, the local electrophilicity at the second carbon atom of the acetylenic moiety is equal to 0.04 eV.Therefore, in light of the Molecular Electron Density Theory [45], the course of the addition reaction with the participation of ethyl propiolate should be determined by the attack of nucleophilic agents on the terminal position of the ethyne moiety.On the other hand, the considered organic azides are characterized by global nucleophilicities in the range of 2.12-3.45eV.Relatively stronger nucleophilicity properties are observed in the case of EDG-substituted azides, whereas weaker nucleophilicities are assigned to EWG-substituted azides (Table 1).It is interesting that local nucleophilicities on both important reaction centers within the NNN moiety are not significantly different and exist in the range of 0.5-0.7 eV.This suggests that in light of electronic effects, both regioisomeric cycloaddition paths should be kinetically allowed.In the initial step, we decided to analyze the nature of the intermolecular interactions between respective pairs of reagents (Table 1, Figures 1 and 2).It was found that ethyl propiolate is characterized by global electrophilicity equal to 1.48 eV.Therefore, according to the unique Domingo electrophilicity scale [43,44], this component should be treated as a moderate electrophile.The analysis of the distribution of local electrophilicities subsequently shows that the most activated reaction center is located on the terminal carbon atom of acetylenic moiety (ω1 = 0.67 eV).In comparison, the local electrophilicity at the second carbon atom of the acetylenic moiety is equal to 0.04 eV.Therefore, in light of the Molecular Electron Density Theory [45], the course of the addition reaction with the participation of ethyl propiolate should be determined by the attack of nucleophilic agents on the terminal position of the ethyne moiety.On the other hand, the considered organic azides are characterized by global nucleophilicities in the range of 2.12-3.45eV.Relatively stronger nucleophilicity properties are observed in the case of EDG-substituted azides, whereas weaker nucleophilicities are assigned to EWG-substituted azides (Table 1).It is interesting that local nucleophilicities on both important reaction centers within the NNN moiety are not significantly different and exist in the range of 0.5-0.7 eV.This suggests that in light of electronic effects, both regioisomeric cycloaddition paths should be kinetically allowed.In the initial step, we decided to analyze the nature of the intermolecular interactions between respective pairs of reagents (Table 1, Figures 1 and 2).It was found that ethyl propiolate is characterized by global electrophilicity equal to 1.48 eV.Therefore, according to the unique Domingo electrophilicity scale [43,44], this component should be treated as a moderate electrophile.The analysis of the distribution of local electrophilicities subsequently shows that the most activated reaction center is located on the terminal carbon atom of acetylenic moiety (ω1 = 0.67 eV).In comparison, the local electrophilicity at the second carbon atom of the acetylenic moiety is equal to 0.04 eV.Therefore, in light of the Molecular Electron Density Theory [45], the course of the addition reaction with the participation of ethyl propiolate should be determined by the attack of nucleophilic agents on the terminal position of the ethyne moiety.On the other hand, the considered organic azides are characterized by global nucleophilicities in the range of 2.12-3.45eV.Relatively stronger nucleophilicity properties are observed in the case of EDG-substituted azides, whereas weaker nucleophilicities are assigned to EWG-substituted azides (Table 1).It is interesting that local nucleophilicities on both important reaction centers within the NNN moiety are not significantly different and exist in the range of 0.5-0.7 eV.This suggests that in light of electronic effects, both regioisomeric cycloaddition paths should be kinetically allowed.Within the next research step, we explored theoretically possible channels of reactions between aryl azides and ethyl propiolate, starting from the model addition 1c+2.The results of the wb97xd/6-311+G(d)(PCM) calculations show that from a qualitative point of view, enthalpy profiles of reactions leading to 3-phenyl-5-carboethoxy-1,2,3pyrazoles (A) or 3-phenyl-4-carboethoxy-1,2,3-pyrazoles (B) are very similar.In both cases, between the valley of the individual reagents and the valley of the respective cycloadduct, two critical points were detected.These are connected by the existence of the pre-reaction molecular complex (MC) and transition state (TS), respectively.
Interactions between addend molecules initially lead to the formation of the prereaction molecular complex MC (MCA and MCB for paths A and B, respectively).This transformation of the reaction system is accompanied by reductions in enthalpy of 3.2 and 5.0 kcal/mol for MCA and MCB, respectively (Table 2, Figure 3).At the same time, however, a substantial reduction in the entropy of the reaction system is observed.Consequently, Gibbs free energies of the formation of MCs exhibit positive values.This excludes the possibility of the existence of MCs as thermodynamically stable intermediates.In the framework of MCs, substructures of addends adopt specific orientations, which stimulate stabilization via the coulombic interactions (Figure 4).In general, the relative orientation of the reaction centers at this stage determines the further observed regioselectivity.Therefore, localized intermediates can be interpreted as orientation complexes [46].Within MCs, any new bonds are not formed (Table 3).Key interatomic distances of N3-C4 and C5-N1 exist beyond the area of the range typical for new C-N bonds within transition states [47][48][49].These types of MCs were experimentally observed earlier regarding 32CAs of ozone with ethene and ethyne [50,51].It should be underlined that both MCs do not exhibit the nature of the charge-transfer complexes [52,53].This was confirmed by analysis of global electron density transfer (GEDT [54]) values (GEDT = 0.00e).Within the next research step, we explored theoretically possible channels of reactions between aryl azides and ethyl propiolate, starting from the model addition 1c+2.The results of the wb97xd/6-311+G(d)(PCM) calculations show that from a qualitative point of view, enthalpy profiles of reactions leading to 3-phenyl-5-carboethoxy-1,2,3-pyrazoles (A) or 3-phenyl-4-carboethoxy-1,2,3-pyrazoles (B) are very similar.In both cases, between the valley of the individual reagents and the valley of the respective cycloadduct, two critical points were detected.These are connected by the existence of the pre-reaction molecular complex (MC) and transition state (TS), respectively.
Interactions between addend molecules initially lead to the formation of the prereaction molecular complex MC (MCA and MCB for paths A and B, respectively).This transformation of the reaction system is accompanied by reductions in enthalpy of 3.2 and 5.0 kcal/mol for MCA and MCB, respectively (Table 2, Figure 3).At the same time, however, a substantial reduction in the entropy of the reaction system is observed.Consequently, Gibbs free energies of the formation of MCs exhibit positive values.This excludes the possibility of the existence of MCs as thermodynamically stable intermediates.In the framework of MCs, substructures of addends adopt specific orientations, which stimulate stabilization via the coulombic interactions (Figure 4).In general, the relative orientation of the reaction centers at this stage determines the further observed regioselectivity.Therefore, localized intermediates can be interpreted as orientation complexes [46].Within MCs, any new bonds are not formed (Table 3).Key interatomic distances of N3-C4 and C5-N1 exist beyond the area of the range typical for new C-N bonds within transition states [47][48][49].These types of MCs were experimentally observed earlier regarding 32CAs of ozone with ethene and ethyne [50,51].It should be underlined that both MCs do not exhibit the nature of the charge-transfer complexes [52,53].This was confirmed by analysis of global electron density transfer (GEDT [54]) values (GEDT = 0.00e).The further transformation of MCs on both considered reaction paths leads to TS (TSA and TSB on reaction paths A and B, respectively).Kinetic factors clearly favored the formation of the TSA structure, which is connected to the 3-phenyl-5-carboethoxy-1,2,3-triazole 3c.This adduct was detected experimentally in the post-reaction mixture.It is interesting that the competitive reaction path B should be treated as disfavored, but not forbidden from a kinetic point of view, because the activation barrier in this case is only 1.1 kcal/mol lower than in the case of path A. However, one should wonder whether the authors of the work [34] carefully analyzed the post-reaction mass.The published procedure only mentions the evaporation of the solvent and recrystallization of the postreaction mass.No detailed composition analyses were performed using, e.g., the HPLC technique [55][56][57].We believe that the second regioisomer is likely formed there as a minority, with a yield of several percent.This is very likely because the isolated product is produced with a yield of 80-90% and the vast majority of 32CA reactions occur with a total yield close to 95-99%.This hypothesis correlates well with the above-described analysis of the local reactivity of title addends.
Both optimized TSs are characterized by imaginary frequencies (−484.72 and −524.70 cm −1 , respectively) and are structurally similar.In particular, within both TSs, key interatomic distances (N3-C4 and C5-N1) are substantially reduced in comparison to the respective MCs (Table 4).However, the kinetically favored TSA transition state is slightly less unsymmetrical.It should be underlined that both considered TSs exhibit evidently polar natures.This was confirmed by an examination of GEDT values (0.31 and 0.37e for TSA and TSB, respectively).The IRC analysis confirms without any doubt that optimized TSs are connected directly to the valleys of the respective MCs and respective products (Figure 3).All attempts at the localization of hypothetical zwitterionic intermediates were not successful.Therefore, postulated reaction channels C and D cannot be realized in practice.At the same time, however, we detected the possibility of the formation of other zwitterionic structures as a consequence of the interactions between azide 1c and alkyne 2 (Scheme 6).These are not postulated intermediates I1C and I2D, but other, different adducts I3E and I4F (Figures 5 and 6) characterized by an "extended" conformation [58,59].This type of conformation excludes the possibility of the cyclization of zwitterions to cycloadducts.Their conversion to triazoles is realized via dissociation to individual reagents, and, in the next step, via one-step cycloaddition according to the A or B path.Similar dissociations are observed during the rotation around the NNCC bond.The formation of zwitterions I3E and I4F are realized via transition states TSE and TSF, respectively.Within the considered TSs, only one new single bond is formed-N3-C4 in the case of TSE and N1-C5 in the case of TSF.It should be underlined, however, that the formation of intermediates I3E and I4F is not favored from a thermodynamic point of view, because the Gibbs free energies of the formation of I3E and I4F are substantially positive.Next, from the kinetic point of view, both these transformations should be considered forbidden due to the very high activation barriers in comparison to reaction paths A and B.
In a similar way, we examined similar reactions with the participation of azides substituted by different types of EDG or EWG groups.In all cases, we detected analogous molecular mechanisms as in the case of the addition 1c+2.Therefore, the proposed mechanism can be assumed to be general for a defined group of 32CA reactions.
In the last step, we decided to characterize the mechanism of 32CA of phenylazide 1c with ethyl propiolate 2 as a representative model with the bonding evolution theory (BET) [60] study along the more favorable reaction path.In Table 4 and Scheme 7, detailed BET data of the critical points of the reaction are presented.The choice requirement is the change in the electronic structure, such as the creation or disappearance of a basin compared to the previous point.The BET study of 716 points along the reaction path revealed twelve phases following intrinsic reaction coordinates.The most important observations can be summed up as follows: (1) In phases I-II, disynaptic basins V(N1,N2) and V ′ (N1,N2) go through topological changes leading to the creation of a new monosynaptic basin V(N2).First, in P1, disynaptic basin V ′ (N1,N2) disappears, transferring its electron density into V(N1,N2) and integrating 4.20e.Next, in P2, monosynaptic basin V(N2) representing nonbonding electron density is created with a population of 0.62e originating from V(N1,N2).The energy increases by 9.4 kcal/mol.(2) In phases III-VI, events leading to the creation of the first new single bond can be observed.During phase III, the integration of disynaptic basins V(N1,N2) and V(N1,N2) decreases by 0.81e and 0.42e, respectively, increasing the population of monosynaptic basin V(N2).In P3, a pseudoradical center is created on C5 represented by the V(C5) monosynaptic basin, integrating 0.14e.An increase of 9.9 kcal/mol in energy can be observed, and the GEDT is −0.08e.Next, a short-phase V starts with the creation of a pseudoradical center on C4 integrating 0.05e; at the same time, an increase in the population in basins V(C5) and V(N2) can be seen.In P5, monosynaptic basin V(C4) disappears, with its electron density being transferred to the monosynaptic basin V(N3) and integrating 3.48e.(3) Phase VII, d(N1-C5) = 2.104 Å and d(N3-C4) = 1.878Å, starts with the creation of the first new C-N bond by donation of the nonbonding electron density of N3 to C4.The disynaptic basin V(N3,C4) is created with an initial population of 1.45e, while the V(N3) monosynaptic basin's integration decreases to 2.06e, as shown by structures P6 ′ and P6 in Figure 6.(4) In P7, a new monosynaptic basin V ′ (N1) integrating 0.38e is created with its population originating from V(N1) now integrating 3.45e.Meanwhile, a transfer of population from disynaptic basins V(C4,C5) and V ′ (C4,C5) of the C4-C5 double bond to V(N3,C4) and V(C5) can be observed.The GEDT is 0.14e, and the energy decreases by 15.3 kcal/mol.( 5) Phase IX starts with the creation of a monosynaptic basin V ′ (N3), integrating 0.50e, and at the same time the population of the disynaptic basin V(N3,C4) decreases to 1.33e.Next, at P9, basin V ′ (N3) disappears transferring its population back to V(N3,C4), integrating now 2.02e.Additionally, an increase in population of the monosynaptic basin V ′ (N1) of 0.30e originating from V(N1) can be seen.The energy decreases by 20.2 kcal/mol and the GEDT increases to 0.24e.(6) The second new C-N bond is created at P10, d(N1-C5) = 1.720Å and d(N3-C4) = 1.561Å, by shearing of the pseudoradical center C5 and nonbonding electron density of N1.The monosynaptic basin V(C5) and V ′ (N1), both integrating 0.75e, merge to create the new disynaptic basin V(N1,C5) with a population of 1.50e, as shown by structures P10 ′ and P10 in Figure 6. ( 7) Phase XII starts with the disappearance of V ′ (C4,C5) with the transfer of its population to V(C4,C5), now integrating 3.67e.A decrease in energy of 27.3 kcal/mol and a GEDT of 0.29e can be observed.(8) The 32CA of phenylazide 1c with alkyne 2 takes place via a two-stage one-step mechanism of the attack of nucleophilic N3 of the azide on the most electrophilic carbon C4 of the alkyne.
In the same way, the process of N3-C4 zwitterion creation has been evaluated (Table 5, Scheme 8).Analysis of hte 264 points along the IRC revealed four phases, of which the most important aspects are: (1) In P1, a monosynaptic basin V(C5) integrating 0.24e is created with its population originating from two disynaptic basins V(C4,C5) and V ′ (C4,C5), as shown by P1 ′ and P1 in Figure 7.The energy increases by 17.0 kcal/mol and the GEDT is −0.03e.(2) At P2, the most significant topological change can be observed, and a new bond is created by sharing the N3 nonbonding electron density with C4.The monosynaptic basin V(N3) is depopulated and a new disynaptic basin V(N3,C4) is created, integrating 1.68e, as shown by structures P2 ′ and P2 in Figure 7.The GEDT increases to −0.31e, and another significant rise in energy of 17.7 kcal/mol can be seen.( 3) The last phase IV starts with the creation of monosynaptic basin V ′ (N3) integrating 0.69e.The energy decreases by 0.9 kcal/mol, while the GEDT increases to −0.51e.The energy of the final structure of the zwitterion I3E is 0.5 kcal/mol lower and the GEDT reaches −0.53e.In a similar way, we examined similar reactions with the participation of azides substituted by different types of EDG or EWG groups.In all cases, we detected analogous molecular mechanisms as in the case of the addition 1c+2.Therefore, the proposed mechanism can be assumed to be general for a defined group of 32CA reactions.
In the last step, we decided to characterize the mechanism of 32CA of phenylazide 1c with ethyl propiolate 2 as a representative model with the bonding evolution theory (BET) [60] study along the more favorable reaction path.In Table 4 and Scheme 7, detailed BET data of the critical points of the reaction are presented.The choice requirement is the change in the electronic structure, such as the creation or disappearance of a basin compared to the previous point.In the same way, the process of N3-C4 zwitterion creation has been evaluated (Table 5, Scheme 8).Analysis of hte 264 points along the IRC revealed four phases, of which the most important aspects are:      The analysis of the mechanism of formation of less favorable zwitterion N1-C5 (Table 6, Scheme 9) revealed ten phases during which the most important topological changes are:

Figure 1 .
Figure 1.The distribution of local electrophilicities along the azide 1c molecule.

Figure 1 .
Figure 1.The distribution of local electrophilicities along the azide 1c molecule.

Figure 1 .
Figure 1.The distribution of local electrophilicities along the azide 1c molecule.

Figure 2 .
Figure 2. The distribution of local nucleophilicities along the ethyl propiolate 2 molecule.

Figure 2 .
Figure 2. The distribution of local nucleophilicities along the ethyl propiolate 2 molecule.

Figure 3 .
Figure 3. Gibbs free energy profiles for all explored paths of the reaction between phenyl azide 1c and ethyl propiolate 2 in the ethanolic solution according to the wb97xd/6-311+G(d)(PCM) calculations.

Figure 3 .
Figure 3. Gibbs free energy profiles for all explored paths of the reaction between phenyl azide 1c and ethyl propiolate 2 in the ethanolic solution according to the wb97xd/6-311+G(d)(PCM) calculations.

Figure 4 .
Figure 4. Views of critical structures of the 32CA reaction between phenyl azide 1c and ethyl propiolate 2 in the ethanolic solution according to the wb97xd/6-311+G(d)(PCM) calculations.

Figure 5 .
Figure 5. Views of critical structures of the creation of zwitterions via reaction between phenyl azide 1c and ethyl propiolate 2 in the ethanolic solution according to the wb97xd/6-311+G(d)(PCM) calculations.

Figure 5 .
Figure 5. Views of critical structures of the creation of zwitterions via reaction between phenyl azide 1c and ethyl propiolate 2 in the ethanolic solution according to the wb97xd/6-311+G(d)(PCM) calculations.

Figure 6 .
Figure 6.ELF valence basins attractors of structures of the most important topological changes, creation of N3-C4 and N1-C5 single bonds, during the reaction of phenyl azide 1c and ethyl propiolate 2 along the reaction path leading to 3c.Structures with apostrophes represent the first structure before the critical point.The formation of zwitterions I3E and I4F are realized via transition states TSE and TSF, respectively.Within the considered TSs, only one new single bond is formed-N3-C4 in the case of TSE and N1-C5 in the case of TSF.It should be underlined, however, that the formation of intermediates I3E and I4F is not favored from a thermodynamic point of view, because the Gibbs free energies of the formation of I3E and I4F are substantially positive.Next, from the kinetic point of view, both these transformations should be considered forbidden due to the very high activation barriers in comparison to reaction paths A and B.In a similar way, we examined similar reactions with the participation of azides substituted by different types of EDG or EWG groups.In all cases, we detected analogous molecular mechanisms as in the case of the addition 1c+2.Therefore, the proposed mechanism can be assumed to be general for a defined group of 32CA reactions.In the last step, we decided to characterize the mechanism of 32CA of phenylazide 1c with ethyl propiolate 2 as a representative model with the bonding evolution theory (BET)[60] study along the more favorable reaction path.In Table4and Scheme 7, detailed BET data of the critical points of the reaction are presented.The choice requirement is the change in the electronic structure, such as the creation or disappearance of a basin compared to the previous point.

Figure 6 . 24 Scheme 7 .
Figure 6.ELF valence basins attractors of structures of the most important topological changes, creation of N3-C4 and N1-C5 single bonds, during the reaction of phenyl azide 1c and ethyl propiolate 2 along the reaction path leading to 3c.Structures with apostrophes represent the first structure before the critical point.Molecules 2023, 28, x FOR PEER REVIEW 15 of 24

( 1 ) 7 . 24 Scheme 8 .
Scheme 7. Simplified representation of the molecular mechanism of 32CA of phenyl azide 1c and ethyl propiolate 2 by Lewis-like structures based on the topological analysis of the ELF along the reaction path leading to 3c.

Figure 7 .
Figure 7. ELF valence basins attractors of structures of the most important topological changes, creation of negative charge on C5 and creation of N3-C4 single bond, during the reaction of phenyl azide 1c and ethyl propiolate 2 along the reaction path leading to I3E.Structures with apostrophes represent the first structure before the critical point.

Scheme 8 . 24 Scheme 8 .
Scheme 8. Simplified representation of the molecular mechanism of the reaction of phenyl azide 1c and ethyl propiolate 2 by Lewis-like structures based on the topological analysis of ELF along the reaction path leading to I3E.

Figure 7 .
Figure 7. ELF valence basins attractors of structures of the most important topological changes, creation of negative charge on C5 and creation of N3-C4 single bond, during the reaction of phenyl azide 1c and ethyl propiolate 2 along the reaction path leading to I3E.Structures with apostrophes represent the first structure before the critical point.

Figure 7 .
Figure 7. ELF valence basins attractors of structures of the most important topological changes, creation of negative charge on C5 and creation of N3-C4 single bond, during the reaction of phenyl azide 1c and ethyl propiolate 2 along the reaction path leading to I3E.Structures with apostrophes represent the first structure before the critical point.

( 1 )Scheme 9 .
Scheme 9. Simplified representation of the molecular mechanism of the reaction of phenyl azide 1c and ethyl propiolate 2 by Lewis-like structures based on the topological analysis of the ELF along the reaction path leading to I4F.

Figure 8 .
Figure 8. ELF valence basins attractors of structures of the most important topological changes, creation of negative charge on C4, and creation of N1-C5 single bond during the reaction of phenyl azide 1c and ethyl propiolate 2 along the reaction path leading to I4F.Structures with apostrophes represent the first structure before the critical point.

Scheme 9 .Scheme 9 .
Scheme 9. Simplified representation of the molecular mechanism of the reaction of phenyl azide 1c and ethyl propiolate 2 by Lewis-like structures based on the topological analysis of the ELF along the reaction path leading to I4F.

Figure 8 .
Figure 8. ELF valence basins attractors of structures of the most important topological changes, creation of negative charge on C4, and creation of N1-C5 single bond during the reaction of phenyl azide 1c and ethyl propiolate 2 along the reaction path leading to I4F.Structures with apostrophes represent the first structure before the critical point.

Figure 8 .
Figure 8. ELF valence basins attractors of structures of the most important topological changes, creation of negative charge on C4, and creation of N1-C5 single bond during the reaction of phenyl azide 1c and ethyl propiolate 2 along the reaction path leading to I4F.Structures with apostrophes represent the first structure before the critical point.

Table 1 .
Global and local nucleophilicities at azide molecules 1a-c (in eV).

Table 1 .
Global and local nucleophilicities at azide molecules 1a-c (in eV).

Table 1 .
Global and local nucleophilicities at azide molecules 1a-c (in eV).

Table 3 .
Key parameters of critical structures of the 32CA reaction between aryl azides 1a-e and ethyl propiolate 2 in the ethanolic solution according to the wb97xd/6-311+G(d)(PCM) calculations.

Table 4 .
ELF valence basin populations of the IRC points, MC 1c-3c, defining the twelve different phases characterizing the reaction of the phenyl azide 1c and ethyl propiolate 2. The stationary points 1c, 2, MCA, TSA, and 3c are also included.Distances are given in angstroms, Å, electron populations in an average number of electrons, [e], relative energies in kcal•mol −1 , and GEDT values in an average number of electrons, [e].

Table 5 .
ELF valence basin populations of the IRC points, MCE-I3E, defining the four different phases characterizing the reaction of the phenyl azide 1c and ethyl propiolate 2. The stationary points 1c, 2, MCE, TSE, and I3E are also included.Distances are given in angstroms, Å, electron populations in the average number of electrons, [e], relative energies in kcal•mol −1 , and GEDT values in an average number of electrons, [e].

Table 6 .
ELF valence basin populations of the IRC points, MCF-I4F, defining the ten different phases characterizing the reaction of the phenyl azide 1c and ethyl propiolate 2. The stationary points 1c, 2, MCF, TSF, and I4F are also included.Distances are given in angstroms, Å, electron populations in the average number of electrons, [e], relative energies in kcal•mol-1, and GEDT values in an average number of electrons, [e].