key: cord-0058690-cb499qpu authors: Rosi, Marzio; Pacifici, Leonardo; Skouteris, Dimitrios; Caracciolo, Adriana; Casavecchia, Piergiorgio; Falcinelli, Stefano; Balucani, Nadia title: A Computational Study on the Insertion of N((2)D) into a C—H or C—C Bond: The Reactions of N((2)D) with Benzene and Toluene and Their Implications on the Chemistry of Titan date: 2020-08-20 journal: Computational Science and Its Applications - ICCSA 2020 DOI: 10.1007/978-3-030-58808-3_54 sha: 2a1a7bee2ab301917fc3d65a72c61df442ea1433 doc_id: 58690 cord_uid: cb499qpu The reactions between nitrogen atoms in their first electronically excited state (2)D with benzene and toluene have been characterized by electronic structure calculations of the stationary points along the minimum energy path. We focused our attention, in particular, to the channels leading to the imidogen radical for the first reaction implying the insertion of nitrogen into a C—H bond and to the NCH(3) radical for the second reaction implying the insertion of nitrogen into a C—C bond. The minima along these reaction paths have been characterized using different ab initio methods in order to find a reasonable compromise between chemical accuracy and computational costs. Our results suggest that, while for geometry optimizations even relatively low level calculations are adequate, for energies higher level of calculations are necessary in order to obtain accurate quantitative results, in particular when strong correlation effects are present. information provided by the NASA/ESA/ASI Cassini-Huygens mission [5] and the observations performed with the ALMA interferometer [6] . Among the species identified by Cassini Ion Neutral Spectrometer (INMS), benzene shows a relatively important mole fraction, being 1.3 Â 10 −6 at 950 km [7] . Solid benzene has been identified also on the surface of Titan [8] . Toluene is efficiently produced in Titan's atmosphere since C 6 H 5 , the main product of the photodissociation of benzene, reacts with the radicals most abundant in Titan's atmosphere, which are H and CH 3 [9] . However, molecular dinitrogen, in its closed shell electronic configuration, cannot react with hydrocarbons in the atmosphere of Titan because of the presence of relatively high activation energy barriers. In the upper atmosphere of Titan, however, molecular dinitrogen can be converted into atomic nitrogen or ions by energetic processes [10] or by the interaction with Extreme Ultra-Violet (EUV) radiation. The dissociation of molecular dinitrogen induced by dissociative photoionization, galactic cosmic ray absorption, N 2 + dissociative recombination, or dissociation induced by EUV photons produces atomic nitrogen in its electronic ground state 4 S and in the first excited 2 D state in similar amounts [10] . While atomic nitrogen in its 4 S ground state exhibits very low reactivity with closed shell molecules, atomic nitrogen in its first electronically excited 2 D state shows a significant reactivity with several molecules identified in the atmosphere of Titan. Atomic nitrogen in its excited 2 D state is metastable but it shows a radiative lifetime long enough to react in binary collisions with other constituents of the upper atmosphere of Titan (6.1 Â 10 4 s and 1.4 Â 10 5 s for the 2 D 3/2 and 2 D 5/2 state, respectively) [11] [12] [13] [14] [15] [16] [17] [18] [19] . We have already investigated the reactions of atomic nitrogen in its excited 2 D state with aliphatic hydrocarbons, like CH 4 , C 2 H 2 , C 2 H 4 , C 2 H 6 , in laboratory experiments by the crossed molecular beam technique and by ab initio and kinetics calculations [20] [21] [22] [23] [24] [25] . More recently, with the same approach, we started the investigation of the reaction between N( 2 D) and aromatic hydrocarbons like benzene [26] . In this study we will report preliminary calculations on the interaction of N( 2 D) with benzene and toluene paying attention to the computational requirements. Indeed, the ab initio study of these systems is computationally very challenging and a compromise between chemical accuracy and computational resources is necessary. The reactive channels of the N( 2 D) + C 6 H 6 and C 6 H 5 CH 3 systems leading to the insertion of atomic nitrogen into a C-H or C-C bond were investigated by locating the lowest stationary points at the B3LYP level of theory [27, 28] , in conjunction with the 6 − 311 + G** basis set [29, 30], on the doublet ground state potential energy surface. The basis set superposition error (BSSE) has been estimated using the counterpoise method [31, 32] . At the same level of theory we have computed the harmonic vibrational frequencies in order to check the nature of the stationary points, i.e. minimum if all the frequencies are real, saddle point if there is only one imaginary frequency. The assignment of the saddle points was performed using intrinsic reaction coordinate (IRC) calculations [33, 34] The geometry of all the species was optimized without any symmetry constraints considering for all the investigated species the electronic ground state. In order to check the accuracy of the computed geometries, we have optimized all the minima using a more extended basis set. We have optimized all the minima at the B3LYP level [27, 28] with the correlation consistent aug-cc-pVTZ basis set [35] . For all the stationary points, the energy was computed also at the higher level of calculation CCSD(T) [36] [37] [38] using the same correlation consistent aug-cc-pVTZ basis set [35] and the B3LYP/aug-cc-pVTZ optimized geometries, following a well established computational scheme [39] [40] [41] [42] [43] [44] [45] [46] [47] [48] [49] . The accuracy of the employed approach, in particular as far as basis set completeness is concerned, has been recently investigated [50] . Both the B3LYP and the CCSD(T) energies were corrected to 0 K by adding the zero point energy correction computed using the scaled harmonic vibrational frequencies evaluated at B3LYP level. The energy of N( 2 D) was estimated by adding the experimental [51] separation N( 4 S) -N( 2 D) of 55 kcal mol −1 to the energy of N( 4 S) at all levels of calculation. All calculations were done using Gaussian 09 [52] while the analysis of the vibrational frequencies was performed using Molekel [53, 54]. Figure 1 reports the optimized geometries of the stationary points localized along the reactive channel leading from N( 2 D) + benzene to the insertion product of N( 2 D) into a C-H bond and finally to the formation of the imidogen radical NH. We have found four minima along this reaction channel connected by three transition states: the first minimum (1) is only an adduct where N is slightly interacting with the benzene ring; the second one (2) shows a C atom of the ring bonded to an H and the N atoms; the third minimum (3) shows the nitrogen atom bonded to two adjacent carbon atoms and then the fourth minimum (4) shows the insertion of the nitrogen atom into the C-H bond with the formation of an NH group. Figure 1 shows the main geometrical parameters optimized at B3LYP level using both the 6 − 311 + G** and the aug-cc-pVTZ basis set. The last results are reported in parentheses. Some of these optimized structures were previously reported [55] . We can notice that there is a good agreement between the geometries optimized with the 6 − 311 + G** and the aug-cc-pVTZ basis set, the differences being within 0.01Å. There is only one relevant difference represented by the N-C distance in species 1: this distance indeed is 3.044 Å with the 6 − 311 + G** basis set and 3.186 Å with the aug-cc-pVTZ basis set. However, this is an adduct and the energy dependence from the distance in this case is weak, so it is very difficult to compute the correct bond length. In Fig. 2 we have reported the minimum energy path for the channel leading to species 4 and then to the formation of the imidogen radical. The approach of the nitrogen to the benzene ring leads to the adduct 1 which is more stable than the reactants by 30.1 kcal/mol. This species, through the transition state TS12 which is still well under the reactants gives rise to species 2 where we have a true C-N bond. Species 2, through a very small barrier of only 0.1 kcal/mol evolves towards species 3 which shows the nitrogen bridging two adjacent carbon atoms. It is necessary to overcome a barrier of 48.4 kcal/mol for species 3 to evolve to the more stable species 4 which represents the lowest minimum in this path. Species 4 shows the insertion of the nitrogen into a C-H bond. The dissociation of the C-N bond leads to the formation of the imidogen radical in a process endothermic by 101.2 kcal/mol. The imidogen radical can be formed also by direct abstraction of the hydrogen by the N( 2 D). If the nitrogen approaches the benzene ring from the hydrogen side, the adduct 5 is formed which is 30.0 kcal/mol more stable than the reactants. This adduct, through a transition state which is 22.7 kcal/mol under the reactants gives the products. The direct abstraction of hydrogen implies only one step for the formation of the imidogen radical, while the reaction path through species 4 implies several steps which lie however at lower energies. In order to state which one is the preferred path, kinetic calculations are necessary and they will be performed in the future. Optimized geometries at B3LYP/6 − 311 + G** and B3LYP/aug-cc-pVTZ (in parentheses) level of the stationary points considered along the minimum energy paths leading from N( 2 D) + C 6 H 6 to C 6 H 5 NH. Bond lengths in Å. Figure 3 reports the optimized geometries of the stationary points localized along the reactive channel leading from N( 2 D) + toluene to the insertion product of N( 2 D) into a C-C bond and finally to the formation of the NCH 3 radical. We have found three minima along this reaction channel connected by two transition states: the first minimum 1 is only an adduct where N is slightly interacting with the toluene ring, being the C-N distance as high as 3.027 Å; the second one 2 shows the nitrogen atom bonded to two adjacent carbon atoms and then the third minimum 3 shows the insertion of the nitrogen atom into the C-C bond with the formation of an NCH 3 group. Figure 3 shows the main geometrical parameters optimized at B3LYP level using both the 6 − 311 + G** and the aug-cc-pVTZ basis set. The last results are reported in parentheses. We were not able to optimize the initial adduct at the B3LYP/aug-cc-pVTZ level. From Fig. 3 we can notice that there is a good agreement between the geometries optimized with the 6 − 311 + G** and the aug-cc-pVTZ basis set, the differences being within 0.01Å. We can also notice that the presence of the methyl substituent in toluene, with respect to benzene, changes the electronic distribution around the carbon atom, preventing the formation of the equivalent of species 2 reported in Fig. 1. In Fig. 4 we have reported the minimum energy path for the channel leading to species 3 and then to the formation of the NCH 3 radical. The approach of the nitrogen to the ipso carbon of the toluene ring leads to the adduct 1 which is more stable than the reactants by 30.2 kcal/mol. This species, through the transition state TS12 which is still well under the reactants gives rise to species 2 which shows the nitrogen bridging the ipso and the ortho carbon atoms. It is necessary to overcome a barrier of 47.5 kcal/mol for species 2 to evolve to the more stable species 3 which represents the lowest minimum in this path. Species 3 shows the insertion of the nitrogen into a C-C bond. The dissociation of the C-N bond leads to the formation of the NCH 3 radical in a process endothermic by 88.5 kcal/mol. Comparing Fig. 2 and Fig. 4 we can notice that, except for the presence of minimum 2 in Fig. 2 which, however, can evolve to minimum 3 through an almost negligible barrier height, the insertions of N( 2 D) into a C-H or a C-C bond are very similar processes. Moreover, this seems to be a realistic path to the formation of a imidogen or a NCH 3 radical. Fig. 3 . Optimized geometries at B3LYP/6 − 311 + G** and B3LYP/aug-cc-pVTZ (in parentheses) level of the stationary points considered along the minimum energy path leading from N( 2 D) + toluene to C 6 H 5 NCH 3 . Bond lengths in Å. Fig. 2 and Fig. 4 all the species, including the transition states, are in energy well under the reactants, some doubts could emerge due to the low level of accuracy of the performed calculations. For this reason, we performed more accurate calculations using the CCSD(T) method in conjunction with a more extended basis set as the aug-cc-pVTZ. We optimized all the stationary points (minima and saddle points) at the B3LYP/aug-cc-pVTZ level and at these optimized geometries we performed CCSD(T)/aug-cc-pVTZ calculations. The energies of main steps along the minimum energy path for the reaction N( 2 D) + benzene ! C 6 H 5 + NH and N( 2 D) + toluene ! C 6 H 5 + NCH 3 are reported in Table 1 . For the main processes we reported both the enthalpy differences and the barrier heights computed at 0 K. In Table 1 , the first entry is the energy computed at B3LYP/6 − 311 + G**, the second entry the relative energy computed at B3LYP/aug-cc-pVTZ level and the third entry the energy computed at CCSD(T)/aug-cc-pVTZ//B3LYP/aug-cc-pVTZ level. We were not able to optimize the geometry of species 1 for the toluene reaction at the B3LYP/aug-cc-pVTZ level. For this reason, in Table 1 we have reported also the reaction from the reactants to species 2. Comparing the values of Table 1 we can notice that the B3LYP/6 − 311 + G** and the B3LYP/aug-cc-pVTZ energies differ always by a small amount of energy (around 1 kcal/mol). This difference is essentially due to the basis set superposition error computed with the counterpoise method which is 1.1 kcal/mol [31, 32]. Therefore, comparing the geometries and the energies, we can conclude that it is worthless to use a more extended basis set, with respect to the 6 − 311 + G**, to optimize the geometries at B3LYP level. Comparing the B3LYP and CCSD(T) energies we have a complete different situation, being the difference sometimes higher than ten kcal/mol. In particular, the difference between B3LYP and CCSD(T) energies is relevant for reactions involving species 1: this is an expected result since species 1 is an adduct which is very difficult to describe accurately at B3LYP level; moreover, in species 1 there is a N( 2 D) atom almost unperturbed and it is well known that the B3LYP method cannot describe correctly the nitrogen excited 2 D state. For the reactions considered in this context, the situation however does not change dramatically going from B3LYP to CCSD(T) results because the energies of all the species investigated are still well under the reactants, but for other reaction this point could be very relevant. The study at ab initio level of the reactions N( 2 D) + C 6 H 6 ! C 6 H 5 + NH and N( 2 D) + C 6 H 5 CH 3 ! C 6 H 5 + NCH 3 , performed using different methods in order to find a reasonable compromise between chemical accuracy and computational costs Table 1 . Enthalpy changes and barrier heights (kcal/mol, 0 K) computed at the B3LYP/6 − 311 + G**, B3LYP/aug-cc-pVTZ and CCSD(T)/aug-cc-pVTZ//B3LYP/aug-cc-pVTZ level of theory for the minimum energy path for the reactions N( 2 D) + C 6 H 6 ! C 6 H 5 + NH and N( 2 D) + C 6 H 5 CH 3 ! C 6 H 5 + NCH 3 . First entry B3LYP/6 − 311 + G**, second entry B3LYP/aug-cc-pVTZ and third entry CCSD(T)/aug-cc-pVTZ//B3LYP/aug-cc-pVTZ Titan's atmosphere and climate The nitrogen chemistry of Titan's upper atmosphere revealed Titan: Surface, Atmosphere and Magnetosphere Elementary reactions of N atoms with hydrocarbons: first steps towards the formation of prebiotic N-containing molecules in planetary atmospheres Titan from Cassini-Huygens Mapping vinyl cyanide and other nitriles in Titan's atmosphere using ALMA Formation and distribution of benzene on Titan Detection and mapping of hydrocarbon deposits on Titan The photochemical production of aromatics in the atmosphere of Titan Energy deposition and primary chemical products in Titan's upper atmosphere Critical review of N, N + , N 2 + , N ++ and N 2 ++ main production processes and reactions of relevance to Titan's atmosphere Nitrogen fixation by photochemistry in the atmosphere of Titan and implications for prebiotic chemistry Elementary reactions and their role in gas-phase prebiotic chemistry Formation of nitrogenated organic aerosols in the Titan upper atmosphere An improved algorithm for reaction path following Reaction path following in mass-weighted internal coordinates Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen Many-body perturbation theory and coupled cluster theory for electron correlation in molecules Quadratic configuration interaction. A general technique for determining electron correlation energies Full configuration-interaction and state of the art correlation calculations on water in a valence double-zeta basis with polarization functions Experimental detection of theoretically predicted N 2 CO SSOH and HSSO radicals: an experimental and theoretical study of [S 2 OH] 0/+/-species The intermolecular potential in NO-N 2 and (NO-N 2 ) + systems: implications for the neutralization of ionic molecular aggregates Observation of organosulfur products (thiovinoxy, thioketene and thioformyl) in crossed-beam experiments and low temperature rate coefficients for the reaction S( 1 D) + C 2 H 4 Crossed-beam and theoretical studies of the S( 1 D) + C 2 H 2 reaction The proton affinity and gas-phase basicity of sulfur dioxide Low temperature kinetics, crossed beam dynamics and theoretical studies of the reaction S( 1 D) + CH 4 and low temperature kinetics of S( 1 D) + C 2 H 2 Theoretical study of reactions relevant for atmospheric models of Titan: rocarbons Dimerization of methanimine and its charged species in the atmosphere of Titan and interstellar/cometary ice analogs Stereoselectivity in autoionization reactions of hydrogenated molecules by metastable gas atoms: the role of electronic couplings Vanadium hydroxide cluster ions in the gas phase: bond-forming reactions of doubly-charged negative ions by SO 2 -promoted V-O activation An experimental and theoretical investigation of 1-butanol pyrolysis MOLEKEL 4.3, Swiss Center for Scientific Computing MOLEKEL: an interactive molecular graphics tool A computational study of the reaction N( 2 D) + C 6 H 6 leading to pyridine and phenylnitrene Acknowledgments. This work has been supported by MIUR "PRIN 2015" funds, project "STARS in the CAOS (Simulation Tools for Astrochemical Reactivity and Spectroscopy in the