Ab Initio and Quantum Chemical Topology studies on the isomerization of HONO to HNO2. Effect of the basis set in QCT

The article focus on the isomerization of nitrous acid HONO to hydrogen nitryl HNO₂. Density functional (B3LYP) and MP2 methods, and a wide variety of basis sets, have been chosen to investigate the mechanism of this reaction. The results clearly show that there are two possible paths: 1) Uncatalysed isomerisation, trans‐HONO → HNO₂, involving 1,2‐hydrogen shift and characterized by a large energetic barrier 49.7 ÷ 58.9 kcal/mol, 2) Catalysed double hydrogen transfer process, trans‐HONO + cis‐HONO → HNO₂ + cis‐HONO, which displays a significantly lower energetic barrier in a range of 11.6 ÷ 18.9 kcal/mol. Topological analysis of the Electron Localization Function (ELF) shows that the hydrogen transfer for both studied reactions takes place through the formation of a ‘dressed’ proton along the reaction path. 1 Use of a wide variety of basis sets demonstrates a clear basis set dependence on the ELF topology of HNO₂. Less saturated basis sets yield two lone pair basins, V₁(N), V₂(N), whereas more saturated ones (for example aug‐cc‐pVTZ and aug‐cc‐pVQZ) do not indicate a lone pair on the nitrogen atom. Topological analysis of the Electron Localizability Indication (ELI‐D) at the CASSCF (12,10) confirms these findings, showing the existence of the lone pair basins but with decreasing populations as the basis set becomes more saturated (0.35e for the cc‐pVDZ basis set to 0.06e for the aug‐cc‐pVTZ). This confirms that the choice of basis set not only can influence the value of the electron population at the particular atom, but can also lead to different ELF topology. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

A Comparative Study of Hydrogen Bonding Using Density Functional Theory

The hydrogen‐bond energies of water dimer and water‐formaldehyde complexes have been studied using density functional theory (DFT). Basis sets up to aug‐cc‐pVXZ (X=D, T, Q) were used. It was found that counterpoise corrected binding energies using the aug‐cc‐pVDZ basis set are very close to those predicted with the aug‐cc‐pVQZ set. Comparative studies using various DFT functionals on these two systems show that results from B3LYP, mPW1PW91 and PW91PW91 functionals are in better agreements with those predicted using high‐level ab initio methods. These functionals were applied to the study of hydrogen bonding between guanine (G) and cytosine (C), and between adenine (A) and thy mine (T) base pairs. With the aug‐cc‐pVDZ basis set, the predicted binding energies of base pairs are in good agreement with the most elaborate ab initio results.     

Computational NMR coupling constants: Shifting and scaling factors for evaluating 1JCH

Optimized shifting and/or scaling factors for calculating one‐bond carbon–hydrogen spin–spin coupling constants have been determined for 35 combinations of representative functionals (PBE, B3LYP, B3P86, B97‐2 and M06‐L) and basis sets (TZVP, HIII‐su3, EPR‐III, aug‐cc‐pVTZ‐J, ccJ‐pVDZ, ccJ‐pVTZ, ccJ‐pVQZ, pcJ‐2 and pcJ‐3) using 68 organic molecular systems with 88 ¹J_CH couplings including different types of hybridized carbon atoms. Density functional theory assessment for the determination of ¹J_CH coupling constants is examined, comparing the computed and experimental values. The use of shifting constants for obtaining the calculated coupling improves substantially the results, and most models become qualitatively similar. Thus, for the whole set of couplings and for all approaches excluding those using the M06 functional, the root‐mean‐square deviations lie between 4.7 and 16.4 Hz and are reduced to 4–6.5 Hz when shifting constants are considered. Alternatively, when a specific rovibrational contribution of 5 Hz is subtracted from the experimental values, good results are obtained with PBE, B3P86 and B97‐2 functionals in combination with HIII‐su3, aug‐cc‐pVTZ‐J and pcJ‐2 basis sets. Copyright © 2013 John Wiley & Sons, Ltd.     

ChemInform Abstract: Smallest Deltahedra Silicon Dicarbide: C2Si32‐.

A total of 20 singlet and 36 triplet C₂Si₃^2‐ isomers are obtained by quantum chemical calculations at the B3LYP/aug‐cc‐pVTZ level.     

Improving on equilibrium isotropic nuclear shielding constants

Isotropic nuclear shielding constants at the equilibrium molecular structure σ_eq and zero‐point vibrational corrections (ZPVCs) to σ_eq are evaluated using the B3LYP/aug‐cc‐pVTZ level of theory, as well as the KT2/aug‐cc‐pVTZ level of theory. Various scaling factors and systematic corrections are obtained by linear regression to experimental shielding constants. Comparisons of the scaled and systematically corrected equilibrium and vibrationally averaged shielding constants reveal that, at the 99% confidence level, the ZPVCs via second‐order perturbation theory do not improve the agreement of B3LYP/aug‐cc‐pVTZ and KT2/aug‐cc‐pVTZ calculated shielding constants with experiment. This holds true when the same analysis is applied to CCSD(T)/aug‐cc‐pCV[TQ]Z calculated σ_eq of Teale et al. [Journal of Chemical Physics 2013, 138, 024111]. In addition, at the 99% confidence level, B3LYP/aug‐cc‐pVTZ and KT2/aug‐cc‐pVTZ scaled and systematically corrected shielding constants are found to be statistically no different from CCSD(T)/aug‐cc‐pCV[TQ]Z calculated σ_eq. The use of scaling factors and systematic corrections could thus provide a cheaper but yet reasonably accurate alternative for the study of nuclear shielding constants of larger systems.     

Compatibility of correlation‐consistent basis sets with a hybrid Hartree–Fock/density functional method

The present study examines the feasibility of combining the correlation‐consistent basis sets developed by Dunning and coworkers with the hybrid Hartree–Fock/density functional method B3LYP. Furthermore, extrapolation to the complete basis set (CBS) limit minimizes errors due to the presence of an incomplete basis set and can act as a rigorous test of the limitations of the B3LYP method. Equilibrium geometries, energies, and harmonic vibrational frequencies were determined for a series of well‐studied, yet computationally challenging, small inorganics and their respective ions. The results were then extrapolated to the CBS limit, where applicable, and compared to experiment. It was found that a union between the hybrid Hartree–Fock/density functional B3LYP method and Dunning's augmented correlation‐consistent basis sets gave results that were comparable to molecular orbital methods that explicitly account for electron correlation. Furthermore, the minimum basis set necessary to attain reasonable results for the systems studied was aug‐cc‐pVTZ. Upgrading to the aug‐cc‐pVQZ level and subsequent extrapolation to the CBS limit further improved the overall agreement with the experiment. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 207–216, 1999     

Comparative study on formamide–water complex

The hydrogen bonding of 1:1 complexes formed between formamide and water molecule have been investigated systematically using Hartree–Fock (HF), hybrid density functional theory (B3LYP), and post‐Hartree–Fock (MP2 and CCSD(T)) methods with range of basis sets 6‐31G(d), cc‐pVXZ (X = D, T, Q) and aug‐cc‐pVYZ (Y = D, T). Three stable structures are considered on the potential energy surface of formamide and water system. The optimized geometric parameters and interaction energies for various isomers at different levels are estimated. The IR frequencies, intensities, and frequency shifts are reported. This study shows that B3LYP/aug‐cc‐pVDZ method gives better performance for formamide‐water complexes. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010.     

Highly accurate benchmark calculations of the interaction energies in the complexes C6H6···C6X6 (X = F,Cl, Br,and I)

Different from the case of the benzene dimer, the differences between the interaction energies are always less than 0.50 kcal/mol for face‐to‐face eclipsed, face‐to‐face staggered, and parallel‐displaced configurations of all investigated complexes C₆H₆···C₆X₆ (X = F, Cl, Br, and I). Hence, it is a great challenge for quantum chemists to accurately calculate the interaction energies for the three configurations of the complexes C₆H₆···C₆X₆. This work demonstrates that results obtained with the PBE0 density functional combined with the D3 dispersion correction (PBE0‐D3) and the basis set def2‐TZVPP are in excellent agreement with the estimates of the coupled‐cluster singles, doubles, and perturbative triples [CCSD(T)] complete basis set (CBS) limit. The other finding in this study is that, in comparison with the gold‐standard CCSD(T)/CBS benchmark, the spin‐component scaled (SCS) zeroth‐order symmetry‐adapted perturbation theory (SAPT0), when paired with the basis set aug‐cc‐pVDZ, performs also very well, and its performance is even better than that of the PBE0‐D3/def2‐TZVPP method or the conventional SAPT/aug‐cc‐pVQZ method. The findings of this study are very significant because both PBE0‐D3/def2‐TZVPP and SCS‐SAPT0/aug‐cc‐pVDZ can deal with the systems with more than 200 atoms.     

An analytical approach for computing Franck‐Condon integrals of harmonic oscillators with arbitrary dimensions

We have developed an analytical approach for computing Franck‐Condon integrals (FCIs) of harmonic oscillators (HOs) with arbitrary dimensions in which the mode‐mixing Duschinsky effect is taken into account. A general formula of FCIs of HOs was obtained and was applied to study the photoelectron spectroscopy of vinyl alcohol and ovalene (C₃₂H₁₄). The equilibrium geometries, harmonic vibrational frequencies and normal modes of vinyl alcohol, ovalene, and their cations were computed at the B3LYP/aug‐cc‐pVTZ or the B3LYP/6‐31G(d) level, from which Franck‐Condon factors were calculated and photoelectron spectra were simulated. The adiabatic ionization energies of vinyl alcohol were also computed by extrapolating the CCSD(T) energies to the complete basis set limit with aug‐cc‐pVXZ (X = D, T, Q, 5). The simulated photoelectron spectra of both vinyl alcohol and ovalene are in agreement with the experiments. The computed adiabatic ionization energies of syn‐ and anti‐vinyl alcohol are in consistent with the experiment within 0.008 eV and 0.014 eV, respectively. We show, for the first time, that the analytical approach of computing FCIs is also efficient and promising for the studies of vibronic spectra of macrosystems. © 2012 Wiley Periodicals, Inc.     

Computational thermochemistry of glycolaldehyde

We use a variant of the focal point analysis to refine estimates of the relative energies of the four low‐energy torsional conformers of glycolaldehyde. The most stable form is the cis‐cis structure which enjoys a degree of H‐bonding from hydroxyl H to carbonyl O; here dihedral angles τ₁ (O?C? C? O) and τ₂ (C? C? O? H) both are zero. We optimized structures in both CCSD(T)/aug‐cc‐pVDZ and aug‐cc‐pVTZ; the structures agree within 0.01 Å for bond lengths and 1.0 degrees for valence angles, but the larger basis brings the rotational constants closer to experimental values. According to our extrapolation of CCSD(T) energies evaluated in basis sets ranging to aug‐cc‐pVQZ the trans‐trans form (180°, 180°) has a relative energy of 12.6 kJ/mol. The trans‐gauche conformer (160°, ±75°) is situated at 13.9 kJ/mol and the cis‐trans form (0°, 180°) at 18.9 kJ/mol. Values are corrected for zero point vibrational energy by MP2/aug‐cc‐pVTZ frequencies. Modeling the vibrational spectra is best accomplished by MP2/aug‐cc‐pVTZ with anharmonic corrections. We compute the Watsonian parameters that define the theoretical vibrational‐rotational spectra for the four stable conformers, to assist the search for these species in the interstellar medium. Six transition states are located by G4 and CBS‐QB3 methods as well as extrapolation using energies for structures optimized in CCSD(T)/aug‐cc‐pVDZ structures. We use two isodesmic reactions with two well‐established thermochemical computational schemes G4 and CBS‐QB3 to estimate energy enthalpy and Gibbs energy of formation as well as the entropy of the gas phase system. Our extrapolated electronic energies of species appearing in the isodesmic reactions produce independent values of thermodynamic quantities consistent with G4 and CBS‐QB3. © 2013 Wiley Periodicals, Inc.     

Benchmarking density functionals in conjunction with Grimme's dispersion correction for noble gas dimers (Ne2, Ar2, Kr2, Xe2, Rn2)

Eleven exchange‐correlational functionals of different types corrected for dispersion by Grimme's D3 correction in conjunction with the aug‐cc‐pVTZ basis set were tested on the following noble gas (Ng) dimers: Ne₂, Ar₂, Kr₂, Xe₂, and Rn₂. For comparison, the D2 and D3BJ corrections were probed with the B3LYP functional. From post‐HF wavefunction methods, CCSD(T) theory was also included. The investigated properties involved potential energy curves, equilibrium bond distances, and interaction energies. The B3LYP‐D3, B3LYP‐D3BJ, and PBE0‐D3 functionals performed overall best for bond distances, while B3LYP‐D3 and B97‐D3 performed best for interaction energies. The importance of fortunate error cancellations was seen in the often reduced agreement with reference data upon correction for BSSE. As several functionals performed well selectively for some noble gases (and poorly for others), we also analysed the performance on the Ng₂ dimers individually and recommended DFT‐D3 functionals for the calculation of large clusters of each Ng.     

Structures,energetics, and isomerism of [Be,C,O,S]: Stability of triply bonded sulfur

Multiply bonded sulfur has continued to attract attention both experimentally and theoretically. Triply sulfur‐bonded compounds are still rare, due to either the lack of suitable generation precursors or the conversion instability toward doubly sulfur‐bonded structures. A detailed computational study was performed on the structures and stability of various [Be,C,O,S] isomers at the coupled cluster singles doubles (triple excitations) (CCSD(T))/aug‐cc‐pVTZ//B3LYP/6‐311+G(d)+ZPVE level to predict intrinsically stable isomers with triply bonded sulfur. The molecular orbital, bond distance, and harmonic vibrational frequency analysis were carried out at aug‐cc‐pVTZ‐B3LYP, M06‐2X, and CCSD(T) levels to investigate the bonding nature of linear structures. It was shown that two low‐lying isomers are linear SBeCO 01 (0.0 kcal/mol) and SBeOC 02 (15.7 kcal/mol), both of which possess the SBe triple bonding. The Lewis acid–base association of SBe + CO can barrierlessly form 01 and 02, with the former more abundant, while the insertion reaction of SCO + Be might generate more 02 than 01 via photochemical processes. By contrast, formation of the SC‐bearing isomer SCBeO 04 (39.4 kcal/mol) seems unlikely due to its higher energy and less kinetic competition than that of 01 and 02, via either simple association or insertion reactions. The new stable isomers SBeCO 01 and SBeOC 02 add to the number of SBe triply bonded species. Their unique structures and varied branching ratios under association and insertion processes deserve future experimental study. © 2013 Wiley Periodicals, Inc.     

Conventional strain energies of azetidine and phosphetane: Can density functional theory yield reliable results?

The conventional strain energies for azetidine and phosphetane are determined within the isodesmic, homodesmotic, and hyperhomodesmotic models. Optimum equilibrium geometries, harmonic vibrational frequencies, and corresponding electronic energies and zero‐point vibrational energies are computed for all pertinent molecular systems using self‐consistent field theory, second‐order perturbation theory, and density functional theory and using the correlation consistent basis sets cc‐pVDZ, cc‐pVTZ, and cc‐pVQZ. Single point fourth‐order perturbation theory, CCSD, and CCSD(T) calculations using the cc‐pVTZ and the cc‐pVQZ basis sets are computed using the MP2/cc‐pVTZ and MP2/cc‐pVQZ optimized geometries, respectively, to ascertain the contribution of higher order correlation effects and to determine if the quadruple‐zeta valence basis set is needed when higher order correlation is included. In the density functional theory study, eight different functionals are used including B3LYP, wB97XD, and M06‐2X to determine if any functional can yield results similar to those obtained at the CCSD(T) level. © 2012 Wiley Periodicals, Inc.     

Structure and NMR spectra of cyclophanes with unsaturated bridges (cyclophenes)

The calculated structures of several known and hypothetical cyclophanes with ethylene bridges (cyclophenes) are reported together with experimental and calculated values of their NMR parameters. Of the exchange‐correlation functionals and basis sets used in this work, only the ωB97X‐D/6‐311++G(2d,2p) and ωB97X‐D/cc‐pVQZ yielded values of the C_sp3–C_sp3 bond length close to the experimental data, although significant differences still remain. As far as the NMR parameters are concerned, except for close‐lying signals, chemical shifts and coupling constants calculated at the ωB97X‐D/cc‐pVQZ level reproduce in most cases the experimental trends. Contrary to the calculations of geometries, an agreement between the values of the NMR parameters obtained at ωB97X‐D/cc‐pVQZ level and the experimental ones is the poorest compared with that of the ωB97X‐D/6‐311++G(2d,2p) one. Taking into account that the results of the different calculations show the same qualitative trends in most cases, we believe that they correctly describe the structure and properties of the hypothetical molecules studied here. Copyright © 2012 John Wiley & Sons, Ltd.     

A Theoretical Study of the Kinetics of the Hydrogen Atom Abstraction Reactions from Cyclopropane by H,O (3P), and Cl (2P3/2) Atoms and OH Radicals

The rate constants of the H‐abstraction reactions from cyclopropane by H, O (³P), Cl (²P_3/2), and OH radicals have been calculated over the temperature range of 250?2500 K using two different levels of theory. Calculations of optimized geometrical parameters and vibrational frequencies are performed using the MP2 method combined with the cc‐pVTZ basis set and the 6–311++G(d,p) basis set. Single‐point energy calculations have been carried out with the highly correlated ab initio coupled cluster method in the space of single, double, and triple (perturbatively) electron excitations CCSD(T) using either the cc‐pVTZ, aug‐cc‐pVTZ, and aug‐cc‐pVQZ basis sets or the 6–311++G(3df,3pd) basis set. The CCSD(T) calculated potential energies have been extrapolated to the complete basis limit (CBS) limit. The Full Configuration Interaction (FCI) energies have been also estimated using the continued‐fraction approximation as proposed by Goodson (J. Chem. Phys., 2002, 116, 6948–6956). Canonical transition‐state theory combined with an Eckart tunneling correction has been used to predict the rate constants as a function of temperature using two kinetic models (direct abstraction or complex mechanism) at two levels of theory (CCSD(T)‐cf/CBS//MP2/cc‐pVTZ and CCSD(T)‐cf/6–311++G(3df,3pd)//MP2/6–311++G(d,p)). The calculated kinetic parameters are in reasonable agreement with their literature counterparts for all reactions. In the light of these trends, the use of the Pople‐style basis sets for studying the reactivity of other systems such as larger cycloalkanes or halogenated cycloalkanes is recommended because the 6–311++G(3df,3pd) basis set is less time consuming than the aug‐cc‐pVQZ basis set. Based on our calculations performed at the CCSD(T)‐cf/CBS//MP2/cc‐pVTZ level of theory, the standard enthalpy of formation at 298 K for the cyclopropyl radical has been reassessed and its value is (290.5 ± 1.6) kJ mol^?1.     

Highly Accurate CCSD(T) and DFT–SAPT Stabilization Energies of H‐Bonded and Stacked Structures of the Uracil Dimer

The CCSD(T) interaction energies for the H‐bonded and stacked structures of the uracil dimer are determined at the aug‐cc‐pVDZ and aug‐cc‐pVTZ levels. On the basis of these calculations we can construct the CCSD(T) interaction energies at the complete basis set (CBS) limit. The most accurate energies, based either on direct extrapolation of the CCSD(T) correlation energies obtained with the aug‐cc‐pVDZ and aug‐cc‐pVTZ basis sets or on the sum of extrapolated MP2 interaction energies (from aug‐cc‐pVTZ and aug‐cc‐pVQZ basis sets) and extrapolated ΔCCSD(T) correction terms [difference between CCSD(T) and MP2 interaction energies] differ only slightly, which demonstrates the reliability and robustness of both techniques. The latter values, which represent new standards for the H‐bonding and stacking structures of the uracil dimer, differ from the previously published data for the S22 set by a small amount. This suggests that interaction energies of the S22 set are generated with chemical accuracy. The most accurate CCSD(T)/CBS interaction energies are compared with interaction energies obtained from various computational procedures, namely the SCS–MP2 (SCS: spin‐component‐scaled), SCS(MI)–MP2 (MI: molecular interaction), MP3, dispersion‐augmented DFT (DFT–D), M06–2X, and DFT–SAPT (SAPT: symmetry‐adapted perturbation theory) methods. Among these techniques, the best results are obtained with the SCS(MI)–MP2 method. Remarkably good binding energies are also obtained with the DFT–SAPT method. Both DFT techniques tested yield similarly good interaction energies. The large magnitude of the stacking energy for the uracil dimer, compared to that of the benzene dimer, is explained by attractive electrostatic interactions present in the stacked uracil dimer. These interactions force both subsystems to approach each other and the dispersion energy benefits from a shorter intersystem separation.     

Global exploration of isomers and isomerization channels on the quantum chemical potential energy surface of H3CNO3

Global exploration of isomers and isomerization channels on the quantum chemical potential energy surface (PES) is performed for H₃CNO₃ using the Scaled Hypersphere Search‐Anharmonic Downward Distortion Following (SHS‐ADDF) method. The molecular formula of H₃CNO₃ includes functional groups of CH₃, OH, NH₂, COOH, NO, NO₂, and NO₃, which are very important in connection with amino acids and NOx. Geometrical structures and interconversion pathways are disclosed after 18719781 force calculations and 534726 Hessian calculations at the level of B3LYP/6‐31G(d). The explored results are confirmed to be valid, especially for the important lower energy regions, by re‐optimization at the higher level of B3LYP/6‐311++G(d,p). A global reaction route‐mapping using SHS‐ADDF demonstrates the entire view and undeveloped landscapes on PES of H₃CNO₃. Typical compounds of H₃CNO₃, aminoxy formic acid, hydroxycarbamic acid, aminoperformic acid, hydroxymethyl nitrite, nitromethanol, methyl nitrate, methyl peroxynitrite, and dioxaziridine, are well separated from others by very high energy‐barriers. The stable‐most conformer of H₃CNO₃ is difficult to be determined, because of seven structures existing with nearly the same energies within 5.7 kJ/mol at the level of CCSD(T)/aug‐cc‐pVTZ. © 2017 Wiley Periodicals, Inc.     

Thermochemistry for enthalpies and reaction paths of nitrous acid isomers

Recent studies show that nitrous acid, HONO, a significant precursor of the hydroxyl radical in the atmosphere, is formed during the photolysis of nitrogen dioxide in soils. The term nitrous acid is largely used interchangeably in the atmospheric literature, and the analytical methods employed do not often distinguish between the HONO structure (nitrous acid) and HNO₂ (nitryl hydride or isonitrous acid). The objective of this study is to determine the thermochemistry of the HNO₂ isomer, which has not been determined experimentally, and to evaluate its thermal and atmospheric stability relative to HONO. The thermochemistry of these isomers is also needed for reference and internal consistency in the calculation of larger nitrite and nitryl systems. We review, evaluate, and compare the thermochemical properties of several small nitric oxide and hydrogen nitrogen oxide molecules. The enthalpies of HONO and HNO₂ are calculated using computational chemistry with the following methods of analysis for the atomization, isomerization, and work reactions using closed‐ and open‐shell reference molecules. Three high‐level composite methods G3, CBS‐QB3, and CBS‐APNO are used for the computation of enthalpy. The enthalpy of formation, ΔH^o_f(298 K), for HONO is determined as ?18.90 ± 0.05 kcal mol^?1 (?79.08 ± 0.2 kJ mol^?1) and as ?10.90 ± 0.05 kcal mol^?1 (?45.61 ± 0.2 kJ mol^?1) for nitryl hydride (HNO₂), which is significantly higher than values used in recent NO_x combustion mechanisms. H‐NO₂ is the weakest bond in isonitrous acid; but HNO₂ will isomerize to HONO with a similar barrier to the HO? NO bond energy; thus, it also serves as a source of OH in atmospheric chemistry. Kinetics of the isomerization is determined; a potential energy diagram of H/N/O₂ system is presented, and an analysis of the triplet surface is initiated. © 2007 Wiley Periodicals, Inc. Int J Chem Kinet 39: 378–398, 2007     

Isomerization study of C5H5NO molecules

The complex potential energy surface (PES) for the isomerization of C₅H₅NO species, including 18 isomers and 23 interconversion transition states, is probed theoretically at the B3LYP/6‐311++G(d,p) and MP2//B3LYP/6‐311++G(d,p) levels of theory. The geometries and relative energies for various stationary points were determined. The zero‐point vibrational energy (ZPVE) corrections have been made to calculate the reliable energy. We predicted a six‐membered ring structure as a global minima isomer I, which is 118.49 and 131.48 kcal · mol^?1 more stable than the least stable, four‐ and three‐membered ring isomer VIII at B3LYP and MP2//B3LYP levels of theory, respectively. The isomers and interconversion transition states have verified by frequency calculation. The intrinsic reaction coordinates (IRC) calculations have been performed to confirm that each transition state is linked by the desired reactants and products. The isomer stability has been studied using relative energies, chemical hardness, and chemical potential. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

On the mechanism of intramolecular nitrogen‐atom hopping in the carbon chain of C6N radical: A Plausible 3c−4e crossover Long‐Bond

Linear isomers of C₆N radical differ in the position of the nitrogen atom in the carbon chain of C₆N. Reaction routes, involving intramolecular nitrogen atom insertion at varying position in the carbon chain of C₆N, are analyzed for the isomerisation between linear isomers of C₆N. Through an automated and systematic search performed with global reaction route mapping of the potential energy surface, thermal isomerisation pathways for C₆N radical are proposed based on the computations carried out at CASSCF/aug‐cc‐pVTZ, and CCSD(T)/6‐311++G(d,p)//B3LYP/6‐311++G(d,p) levels of the theory. Notably, a high lying linear isomer, centrosymmetric with respect to the nitrogen atom, is observed to be stabilized by a unique crossover three center‐four electron π long bond between the carbon atoms that are spatially separated by a nitrogen atom in a natural bond orbital. This long bond is concluded to be responsible for the predicted thermal isomerisation to be more feasible than the dissociation during the isomerisation pathway of a linear isomer of C₆N. © 2014 Wiley Periodicals, Inc.