Deciphering the Reactive Pathways of Competitive Reactions inside Carbon Nanotubes

Nanoscale control of chemical reactivity, manipulation of reaction pathways, and ultimately driving the outcome of chemical reactions are quickly becoming reality. A variety of tools are concurring to establish such capability. The confinement of guest molecules inside nanoreactors, such as the hollow nanostructures of carbon nanotubes (CNTs), is a straightforward and highly fascinating approach. It mechanically hinders some molecular movements but also decreases the free energy of translation of the system with respect to that of a macroscopic solution. Here, we examined, at the quantum mechanics/molecular mechanics (QM/MM) level, the effect of confinement inside CNTs on nucleophilic substitution (SN2) and elimination (syn-E2 and anti-E2) using as a model system the reaction between ethyl chloride and chloride. Our results show that the three reaction mechanisms are kinetically and thermodynamically affected by the CNT host. The size of the nanoreactor, i.e., the CNT diameter, represents the key factor to control the energy profiles of the reactions. A careful analysis of the interactions between the CNTs and the reactive system allowed us to identify the driving force of the catalytic process. The electrostatic term controls the reaction kinetics in the SN2 and syn/anti-E2 reactions. The van der Waals interactions play an important role in the stabilization of the product of the elimination process.


Introduction
Host-guest chemistry is emerging as a novel way to manipulate reactivity by confining reactive systems (guests) inside containers (hosts). Reaction mechanisms are affected by confinement inside nanoscale reactors such as calixarenes, cucurbiturils, cyclodextrins, metalorganic frameworks (MOF), and zeolites [1]. Confinement acts on the potential energy surface (PES) topology, which is responsible for the kinetics (catalyzing or inhibiting reaction channels) but also changes the thermodynamics of the reaction because of the different interactions that the supramolecular host establishes with the reactants/products of the reaction [2][3][4][5]. The confinement inside a chemical host can have two major effects: the first occurs when the shape/volume of the host cavity (partly) restrains the guest system, and the second occurs when new chemical interactions between the host and guest are established.
Carbon nanotubes (CNTs) are nanosystems formed by trivalent carbon atoms covalently linked to form hollow cylindrical structures. The inside of the cylindrical shape allows them to host molecules [6,7]. CNTs are characterized by great chemical inertia. A crucial difference between CNTs and other nanocontainers is that their low reactivity and high stability make possible reactions under harsh conditions that would destroy other hosts. Experimentally, CNTs have already been used as containers for preparative chemical reactions [8][9][10][11][12][13][14][15][16][17][18][19].
For example, CNT nanoreactors were used to control the regioselectivity of aromatic halogenation reactions [10] and azide-alkyne cycloadditions [13] or to modulate the catalytic activity of metal or metal oxide nanoparticles encapsulated inside the cavity of the CNTs [8,11].
We previously studied the reaction between methyl chloride and a chloride anion (the prototypical S N 2 reaction) inside the confined environment of CNTs of different sizes [39].
Here, considering ethyl chloride as a reactant, we evaluated the effect of CNT confinement in competitive reactions.
Ethyl chloride can react with the chloride anion via either an S N 2 or an E2 mechanism. Scheme 1 shows the S N 2 mechanism (top) and the two pathways of the E2 mechanism, which occur via different stereochemical paths that depend on the syn (syn-E2) and anti (anti-E2) periplanar alignment of the leaving chloride anions and β-proton abstraction. Inside CNTs of different diameters, the requisite of periplanarity can be hard to satisfy and affect the possible reaction mechanisms. In this work, we used computational approaches to understand how the competitive S N 2 and E2 mechanisms are affected by the CNT environment. An energy decomposition analysis was also performed to decipher the effects of the host-CNT environment on the reactivity of the guest systems.
For example, CNT nanoreactors were used to control the regioselectivity of a halogenation reactions [10] and azide-alkyne cycloadditions [13] or to modulate t lytic activity of metal or metal oxide nanoparticles encapsulated inside the cavit CNTs [8,11].
We previously studied the reaction between methyl chloride and a chlorid (the prototypical SN2 reaction) inside the confined environment of CNTs of differe [39]. Here, considering ethyl chloride as a reactant, we evaluated the effect of C finement in competitive reactions.
Ethyl chloride can react with the chloride anion via either an SN2 or an E2 mec Scheme 1 shows the SN2 mechanism (top) and the two pathways of the E2 mec which occur via different stereochemical paths that depend on the syn (syn-E2) (anti-E2) periplanar alignment of the leaving chloride anions and β-proton abs Inside CNTs of different diameters, the requisite of periplanarity can be hard to and affect the possible reaction mechanisms. In this work, we used computatio proaches to understand how the competitive SN2 and E2 mechanisms are affecte CNT environment. An energy decomposition analysis was also performed to deci effects of the host-CNT environment on the reactivity of the guest systems. Scheme 1. SN2, and syn-and anti-E2 reaction mechanisms.

Materials and Methods
The gas-phase reaction profiles were calculated at the DFT level with the functional [44]. The functional was selected on the basis of previous studies [45, basis set used was 6-311++G(2df,2p) [47].
The coordinates of nanotubes 24 Å long were obtained with Nanotube Bui plemented in VMD [48]. CNT lengths greater than this value do not significant the reactivity of guest systems [39].
The reaction profiles inside the CNTs were calculated using the ONIOM available in Gaussian16 [49,50]. Mechanical embedding and electrostatic embeddi considered in the QM/MM calculations [49,[51][52][53]. The high-level layer (i.e., the guest system of ethyl chloride and chloride anion) was described at the M 311++G(2df,2p) level of theory. The low-level layer (i.e., the host CNTs) was d Scheme 1. S N 2, and synand anti-E2 reaction mechanisms.

Materials and Methods
The gas-phase reaction profiles were calculated at the DFT level with the M06-2X functional [44]. The functional was selected on the basis of previous studies [45,46]. The basis set used was 6-311++G(2df,2p) [47].
The coordinates of nanotubes 24 Å long were obtained with Nanotube Builder implemented in VMD [48]. CNT lengths greater than this value do not significantly affect the reactivity of guest systems [39].
The reaction profiles inside the CNTs were calculated using the ONIOM method available in Gaussian16 [49,50]. Mechanical embedding and electrostatic embedding were considered in the QM/MM calculations [49,[51][52][53]. The high-level layer (i.e., the reactive guest system of ethyl chloride and chloride anion) was described at the M06-2X/6-311++G(2df,2p) level of theory. The low-level layer (i.e., the host CNTs) was described using the Universal Force Field (UFF) method [54]. Partial atomic charges calculated using the QEq scheme [55] were used to calculate electrostatic interactions. This methodology was effective to describe the energetics and the geometries of reactions of molecules confined inside CNTs [39], although some approximations are required to reduce the complexity of the system, such as the elimination of counter ions and solvent molecules inside the CNTs.
The structure's minima and saddle points were fully optimized, and frequency calculations determined the nature of these critical points.
The effect of confinement on the reaction path (∆E con f ) was calculated as the difference between the QM/MM energy (reacting molecules in the CNTs) and the QM model system energy (reacting molecules in the gas phase).
Considering the scheme implemented in the ONIOM method, ∆E con f can be decomposed as where ∆E vdW and ∆E elect are the energy differences between host and guest of the van der Waals and electrostatic interactions (comprising both Coulomb interactions and the polarization effects due to the electrostatic embedding scheme). ∆E strain(host) and ∆E strain(guest) measure the energy strain of the host and the guest molecules induced upon the formation of the complex. These terms were calculated similarly to the Activation Strain Model (ASM) proposed by van Zeist and Bickelhaupt [56]. Within this formulation, the effect of confinement on the kinetics and thermodynamics of the S N 2, syn-E2, and anti-E2 mechanisms is decomposed into a sum of contributions with a well-defined chemical meaning, providing a better comprehension of the role of the CNT cavity in the reactive system.

Potential Energy Surfaces for SN2, anti-E2, and syn-E2 Reactions in the Gas Phase
Preliminary DFT calculations were carried out to calculate, at the M06-2X/6-311++G(2df,2p) level of theory, the stationary points of the potential energy surfaces of the three competitive mechanisms in the gas phase ( Figure 1). fined inside CNTs [39], although some approximations are required to reduce th plexity of the system, such as the elimination of counter ions and solvent molecule the CNTs.
The structure's minima and saddle points were fully optimized, and freque culations determined the nature of these critical points.
The effect of confinement on the reaction path (∆ ) was calculated as the ence between the QM/MM energy (reacting molecules in the CNTs) and the QM system energy (reacting molecules in the gas phase).
Considering the scheme implemented in the ONIOM method, ∆ can be posed as where ∆ and ∆ are the energy differences between host and guest of der Waals and electrostatic interactions (comprising both Coulomb interactions polarization effects due to the electrostatic embedding scheme). ∆ (ℎ ∆ ( ) measure the energy strain of the host and the guest molecules i upon the formation of the complex. These terms were calculated similarly to the tion Strain Model (ASM) proposed by van Zeist and Bickelhaupt [56].
Within this formulation, the effect of confinement on the kinetics and thermo ics of the SN2, syn-E2, and anti-E2 mechanisms is decomposed into a sum of contri with a well-defined chemical meaning, providing a better comprehension of the the CNT cavity in the reactive system.

Potential Energy Surfaces for SN2, anti-E2, and syn-E2 Reactions in the Gas Phase
Preliminary DFT calculations were carried out to calculate, at the M0 311++G(2df,2p) level of theory, the stationary points of the potential energy surface three competitive mechanisms in the gas phase ( Figure 1). The approaching nucleophile Cl − forms a preliminary complex (Rx) with ethyl chloride. From this complex, the system can follow three main paths: S N 2, syn-E2, and anti-E2.
The S N 2 mechanism, characterized by the lowest transition state energy (17.6 kcal mol −1 ), is the energetically favorite path and leads to a product (Pd SN2 ) isoenergetic with Rx, since the nucleophile and the leaving group are the same species (Cl − ). Alternatively, the system can undergo an anti-E2 mechanism, with an activation barrier of 30.8 kcal mol −1 , or a syn-E2 mechanism, with an activation barrier of 40.6 kcal mol −1 . The two elimination pathways are endothermic, since Pd E2 is 6.9 kcal mol −1 higher in energy than Pd SN2 , and lead to the formation of an ethylene molecule.
These results were benchmarked (Table 1) against "gold standard" calculations, reported in the literature at the Coupled Cluster CCSD and CCSD(T) levels, with infinite basis set extrapolation, with contributions of inner-shell correlation, scalar relativistic effects, and first-order spin-orbit coupling (CCSD(T)/CBS) [45]. Table 1. Critical point energies (kcal mol −1 ) of the S N 2, and synand anti-E2 mechanisms. Energies were calculated with the M06-2X functional with the 6-311++G(2df,2p) basis set and reference CCSD(T)/CBS values obtained by Bento [45].
The SN2 mechanism, characterized by the lowest transition state energy (17.6 kcal mol −1 ), is the energetically favorite path and leads to a product (PdSN2) isoenergetic with Rx, since the nucleophile and the leaving group are the same species (Cl -). Alternatively, the system can undergo an anti-E2 mechanism, with an activation barrier of 30.8 kcal mol −1 , or a syn-E2 mechanism, with an activation barrier of 40.6 kcal mol −1 . The two elimination pathways are endothermic, since PdE2 is 6.9 kcal mol −1 higher in energy than PdSN2, and lead to the formation of an ethylene molecule.
These results were benchmarked (Table 1) against "gold standard" calculations, reported in the literature at the Coupled Cluster CCSD and CCSD(T) levels, with infinite basis set extrapolation, with contributions of inner-shell correlation, scalar relativistic effects, and first-order spin-orbit coupling (CCSD(T)/CBS) [45]. Table 1. Critical point energies (kcal mol −1 ) of the SN2, and syn-and anti-E2 mechanisms. Energies were calculated with the M06-2X functional with the 6-311++G(2df,2p) basis set and reference CCSD(T)/CBS values obtained by Bento [45].

Potential Energy Surfaces for SN2, anti-E2, and syn-E2 Reactions Confined in CNTs
We investigated the potential energy surfaces for SN2, and syn-and anti-E2 mechanisms inside CNTs of different diameters (Figure 2), namely, CNT(6,6), CNT(7,7), CNT (8,8), CNT (9,9), CNT(10,10), and CNT(12,12) ( Figure S1). The critical points corresponding to the reactant complex (Rx) and the SN2, and synand anti-E2 transition states (TSs) and products (Pds) were fully optimized inside the different CNTs. The activation (∆ ‡ ) and reaction (∆ ) energies of the reactions inside the different-diameter CNTs are reported in Tables 2 and 3. The interactions between the CNT "nanoreactor" and the guest molecules were significantly large. Confinement could lead to either a decrease (catalysis) or increase (inhibition) in the activation energies (Table 2). The critical points corresponding to the reactant complex (Rx) and the S N 2, and synand anti-E2 transition states (TSs) and products (Pds) were fully optimized inside the different CNTs. The activation (∆E ‡ ) and reaction (∆E) energies of the reactions inside the different-diameter CNTs are reported in Tables 2 and 3. The interactions between the CNT "nanoreactor" and the guest molecules were significantly large. Confinement could lead to either a decrease (catalysis) or increase (inhibition) in the activation energies (Table 2). The interaction with the CNTs could also affect the thermodynamics of the reactions (Table 3) and, in particular, the E2 profile that was characterized by reactants and products with different geometries, which implies that locally, an energy difference is established among the isoenergetic reactants/products of the gas phase.
To better understand the effect of confinement in the substitution and elimination mechanisms, we discuss the two mechanisms separately in the next sections.

Substitution Mechanism (S N 2)
The computed activation energies for the S N 2 reaction inside the different tubes are reported in Table 2. A comparison with activation energy ∆E ‡ in the gas phase (17.6 kcal mol −1 ) is shown in Figure 3, indicating the tendency of the reaction to be catalyzed or inhibited by the confinement within CNTs of different sizes. (Table 3) and, in particular, the E2 profile that was char with different geometries, which implies that locally, among the isoenergetic reactants/products of the gas p To better understand the effect of confinement i mechanisms, we discuss the two mechanisms separate The computed activation energies for the SN2 rea reported in Table 2. A comparison with activation ene mol −1 ) is shown in Figure 3, indicating the tendency o hibited by the confinement within CNTs of different s  Figure 3 shows that activation barrier ∆ ‡ incre (6,6), ∆∆ ‡ = 6.0 kcal mol −1 ) and for the larger tubes (C of 4.2 and 11.5 kcal mol −1 , respectively). On the contra by CNT (7,7), CNT (8,8), and CNT (9,9), showing decre 2.3, and 3.2 kcal mol −1 . Illustrative pictures are provide  Figure 3 shows that activation barrier ∆E ‡ increased for the narrowest tube (CNT (6,6), ∆∆E ‡ = 6.0 kcal mol −1 ) and for the larger tubes (CNT (10,10) and CNT (12,12), ∆∆E ‡ of 4.2 and 11.5 kcal mol −1 , respectively). On the contrary, the S N 2 reaction was catalyzed by CNT(7,7), CNT (8,8), and CNT (9,9), showing decreases in the activation energy of 0.9, 2.3, and 3.2 kcal mol −1 . Illustrative pictures are provided in Figure 4. An energy decomposition analysis of these terms (∆∆E ‡ in Table 4) was carried out to shed light on the contributions that affect the S N 2 reaction confined inside CNTs.  . Geometries of Rx and TS of the SN2 mechanism inside (6,6), (9,9), and (12,12) CNTs.  The distortion of the host, ΔEstrain(host), did not contribute to the modification of the kinetics of the SN2 reaction inside the CNTs, as it was close to 0 kcal mol −1 for all the tubes. This is not surprising because of the rigidity of the CNT walls. The energy term related to the distortion of the guest molecules, i.e., ΔEstrain(guest), was also small. This effect can be associated with the rigidity of the guest system inside the CNT cavity. In general, the distortion of the reactive complex, imposed by the presence of the tube, is larger in Rx than in the transition state, because the geometry of the TS is less flexible (than the reactant complex); therefore, the imposed strain is lower. The higher destabilization of Rx compared with that of the TS gives a net stabilization energy.
In the largest tube (CNT(12,12)), the available volume was large enough to allow the reacting molecules to maintain the original arrangement observed in the gas-phase system; thus, the ΔEstrain(guest) term was zero.
When we analyze the chemical interactions occurring in the host/guest system, we can decompose it in two terms: the vdW term (ΔEvdW) and the electrostatic term (ΔEelect). Of course, the two terms act synergistically.
From the decomposition analysis, it appears that the electrostatic interaction governs the CNT effect on the SN2 kinetics. For CNT(6,6), CNT(10,10), and CNT (12,12), ΔEelect was detrimental (5.5, 4.2, and 11.6 kcal mol −1 ), while for CNT(7,7), CNT (8,8), and CNT (9,9), the electrostatic interaction between the tube and the guest contributed to the catalyzing effect by −0.7, −1.9, and −4.0 kcal mol −1 . This term takes into consideration the polarization induced on the tube by the reacting molecules.  The distortion of the host, ∆E strain(host) , did not contribute to the modification of the kinetics of the S N 2 reaction inside the CNTs, as it was close to 0 kcal mol −1 for all the tubes. This is not surprising because of the rigidity of the CNT walls. The energy term related to the distortion of the guest molecules, i.e., ∆E strain(guest) , was also small. This effect can be associated with the rigidity of the guest system inside the CNT cavity. In general, the distortion of the reactive complex, imposed by the presence of the tube, is larger in Rx than in the transition state, because the geometry of the TS is less flexible (than the reactant complex); therefore, the imposed strain is lower. The higher destabilization of Rx compared with that of the TS gives a net stabilization energy.
In the largest tube (CNT (12,12)), the available volume was large enough to allow the reacting molecules to maintain the original arrangement observed in the gas-phase system; thus, the ∆E strain(guest) term was zero.
When we analyze the chemical interactions occurring in the host/guest system, we can decompose it in two terms: the vdW term (∆E vdW ) and the electrostatic term (∆E elect ). Of course, the two terms act synergistically.

Elimination Mechanism (E2)
CNTs of different diameters have catalyzing/inhibiting effects also on the elimination mechanism. Figure 5 compares, with reference to the gas phase, the variation in the activation barrier of the two E2 mechanisms (anti-E2 and syn-E2) inside CNTs of different diameters.

Elimination Mechanism (E2)
CNTs of different diameters have catalyzing/inhibiting effects also on the elimination mechanism. Figure 5 compares, with reference to the gas phase, the variation in the activation barrier of the two E2 mechanisms (anti-E2 and syn-E2) inside CNTs of different diameters. The effect of confinement in the tubes is quite different for the two paths (anti-E2 and syn-E2 mechanisms). In fact, the anti-E2 mechanism is less sensitive to confinement. The main reason is the different orientation of the chlorine anion in the deprotonation step ( Figure 6).
In the anti-E2 mechanism, the antiperiplanar geometry of the C-H and C-Cl bonds allows the TS to reorient inside the tube and to optimize its interaction with the tube walls (in particular with the chlorine atoms), compared with the syn disposition, which is blocked. Figure 6. Geometries of Rx and TS in the anti-E2 mechanism inside (6,6), (9,9), and (12,12) CNTs.
Two situations are particularly interesting: (i) inside CNT (9,9), the geometry of the TS was a perfect fit for the tube, and the interaction was stabilizing (catalysis); (ii) inside CNT (12,12), the tube was too large to interact at the same time with the two chlorines of the TS, and the reaction was slowed (inhibition). The effect of confinement in the tubes is quite different for the two paths (anti-E2 and syn-E2 mechanisms). In fact, the anti-E2 mechanism is less sensitive to confinement. The main reason is the different orientation of the chlorine anion in the deprotonation step ( Figure 6).
CNTs of different diameters have catalyzing/inhibiting effects also on the elimination mechanism. Figure 5 compares, with reference to the gas phase, the variation in the activation barrier of the two E2 mechanisms (anti-E2 and syn-E2) inside CNTs of different diameters. The effect of confinement in the tubes is quite different for the two paths (anti-E2 and syn-E2 mechanisms). In fact, the anti-E2 mechanism is less sensitive to confinement. The main reason is the different orientation of the chlorine anion in the deprotonation step ( Figure 6).
In the anti-E2 mechanism, the antiperiplanar geometry of the C-H and C-Cl bonds allows the TS to reorient inside the tube and to optimize its interaction with the tube walls (in particular with the chlorine atoms), compared with the syn disposition, which is blocked. Figure 6. Geometries of Rx and TS in the anti-E2 mechanism inside (6,6), (9,9), and (12,12) CNTs.
Two situations are particularly interesting: (i) inside CNT (9,9), the geometry of the TS was a perfect fit for the tube, and the interaction was stabilizing (catalysis); (ii) inside CNT (12,12), the tube was too large to interact at the same time with the two chlorines of the TS, and the reaction was slowed (inhibition). Figure 6. Geometries of Rx and TS in the anti-E2 mechanism inside (6,6), (9,9), and (12,12) CNTs.
In the anti-E2 mechanism, the antiperiplanar geometry of the C-H and C-Cl bonds allows the TS to reorient inside the tube and to optimize its interaction with the tube walls (in particular with the chlorine atoms), compared with the syn disposition, which is blocked.
Two situations are particularly interesting: (i) inside CNT (9,9), the geometry of the TS was a perfect fit for the tube, and the interaction was stabilizing (catalysis); (ii) inside CNT (12,12), the tube was too large to interact at the same time with the two chlorines of the TS, and the reaction was slowed (inhibition).
On the other hand, the syn-E2 mechanism (Figure 7) has a well-defined conformation inside the CNT that allows both chlorine atoms to interact with the inner surface of the tube. On the other hand, the syn-E2 mechanism (Figure 7) has a well-defined conformation inside the CNT that allows both chlorine atoms to interact with the inner surface of the tube. Figure 7. Geometries of Rx and TS in the syn-E2 mechanism inside (6,6), (8,8), and (12,12) CNTs.
Inside CNT (6,6), in the TS, a distortion of the approach direction of the chlorine anion involved in deprotonation (C-H-Cl angles of 151° vs. 165° in the gas phase) was observed, with consequent destabilization. Inside CNT(7,7), CNT (8,8), and CNT (9,9), the chlorine anion was free enough to maintain its ideal angle of approach. The CNT walls fit well the TS structure, in particular in the case of CNT (8,8), in which the activation barrier was lowered by 8.3 kcal mol −1 .
Again, for the larger CNT(10,10) and CNT (12,12), the stabilization of the TS was lower than that of Rx, and the reaction was globally inhibited.
The analysis of the energy contributions due to confinement inside the tubes (Tables 5 and 6) showed that the electrostatic interactions between the CNTs and the guest system were once again responsible for the catalysis/inhibition of the anti-E2 and syn-E2 mechanisms inside the tubes. Particularly interesting is the value of ΔEstrain(guest) in the case of the syn-E2 mechanism inside CNT (6,6). A considerable value of 4.7 kcal mol −1 appeared and quantitatively described the distortion due to the lack of space for the TS inside the tube, which we previously described qualitatively.   Inside CNT (6,6), in the TS, a distortion of the approach direction of the chlorine anion involved in deprotonation (C-H-Cl angles of 151 • vs. 165 • in the gas phase) was observed, with consequent destabilization. Inside CNT(7,7), CNT (8,8), and CNT (9,9), the chlorine anion was free enough to maintain its ideal angle of approach. The CNT walls fit well the TS structure, in particular in the case of CNT (8,8), in which the activation barrier was lowered by 8.3 kcal mol −1 .
Again, for the larger CNT(10,10) and CNT (12,12), the stabilization of the TS was lower than that of Rx, and the reaction was globally inhibited.
The analysis of the energy contributions due to confinement inside the tubes (Tables 5 and 6) showed that the electrostatic interactions between the CNTs and the guest system were once again responsible for the catalysis/inhibition of the anti-E2 and syn-E2 mechanisms inside the tubes. Particularly interesting is the value of ∆E strain(guest) in the case of the syn-E2 mechanism inside CNT (6,6). A considerable value of 4.7 kcal mol −1 appeared and quantitatively described the distortion due to the lack of space for the TS inside the tube, which we previously described qualitatively.  Since the elimination mechanisms generate a product different from the reactant (in the S N 2 mechanism, product and reactant are the same), the effect of the confinement on the reaction thermodynamics was also investigated.
As reported in Figure 8 and Table 7, CNTs provide the stabilization of the reaction product depending on the CNT size. Since the elimination mechanisms generate a product differ the SN2 mechanism, product and reactant are the same), the effe the reaction thermodynamics was also investigated.
The electrostatic term is always important, but in the prese played by van der Waals interactions, since the conjugated syst interacts with the newly formed alkene via π-π interactions (t QM/MM model, are included in the van der Waals term). The us which likely increases the EvdW term, can potentially alter the pro voring the elimination channels.
The electrostatic term is always important, but in the present case, a crucial role is played by van der Waals interactions, since the conjugated system of the wall strongly interacts with the newly formed alkene via π-π interactions (these interactions, in the QM/MM model, are included in the van der Waals term). The use of multi-walled CNTs, which likely increases the E vdW term, can potentially alter the products of the reaction, favoring the elimination channels.

Conclusions
A hybrid QM/MM method was used to describe the confinement effects on the competing S N 2, syn-E2, and anti-E2 mechanisms of the reaction of ethyl chloride with chloride anion inside CNTs of different diameters. Confinement can have kinetic effects (either catalysis or inhibition), depending on the CNT diameter. The size of the nanoreactor represents the key factor to control the reaction profiles of the reactions. In general, when the geometry of the TS fits well with the size of the tube, the interaction is stabilizing (catalysis), similar to what happens with enzymes in their catalytic sites. When the tube is too narrow or too broad, the reactions are inhibited. In general, shape complementarity is a crucial factor for the interaction between carbon nanomaterials and molecules [57][58][59][60][61][62].
The selection of an appropriate substrate and of a CNT of the proper size may drive competitive reactions toward a target structure. For example, it may favor the elimination products against the substitution products by stabilizing the generated alkene or impeding the Walden inversion, typical of the S N 2 channel.
The energy decomposition analysis of confinement effects in terms of strain of the guest and host systems, vdW, and electrostatic interactions demonstrated that the electrostatic term is the dominant contributor for both catalysis and inhibition. Nonetheless, in the case of elimination products, van der Waals contributions become important to modify the reaction thermodynamics.
To summarize, in this work, we show that the confinement of molecules inside CNTs can significantly alter the pathways of chemical reactions. The most straightforward use of these nanoreactors is to catalyze reactions, but they can also be used to inhibit the reactivity of molecules trapped inside the nanotube cavity. In this way, highly reactive molecules, such as polyynes or π-conjugated [63][64][65][66][67][68] polymers, can be confined and protected inside CNTs of suitable dimensions.

Data Availability Statement:
The data presented in this study are available upon request from the corresponding authors.