Global Investigation on the Potential Energy Surface of CH3CN: Application of the Scaled Hypersphere Search Method

M←NCCH3, M-η2-(NC)-CH3, and CN-M-CH3 Prepared by Reactions of Ce, Sm, Eu, and Lu Atoms with Acetonitrile: Matrix Infrared Spectra and Theoretical Calculations

The reactions of laser-ablated Ce, Sm, Eu, and Lu atoms with acetonitrile were studied by matrix infrared spectra in a neon matrix, and M←NCCH₃, M-η²-(NC)-CH₃, and CN-M-CH₃ were identified with isotopic substitution and quantum chemical calculations. The major product is the insertion complex (CN-M-CH₃), while the end-on and side-on complexes (M←NCCH₃ and M-η²-(NC)-CH₃) are also trapped in the matrix. The CCN antisymmetric stretching mode for Ce-η²-(NC)-CH₃ was observed at 1536.9 cm^-1, which is much lower than the same modes observed for other lanthanides. NBO analysis reveals that Ce exhibits a remarkable 4f-orbital contribution in bonding to N and to C, reconfirming an active 4f-orbital contribution of cerium in bonding in the side-on complex, while the 4f contributions of Sm and Eu to the M-N and M-C bonds are much lower and the 4f orbital of Lu is not involved in bonding.

A Photoionization Reflectron Time-of-Flight Mass Spectrometric Study on the Detection of Ethynamine (HCCNH2 ) and 2H-Azirine (c-H2 CCHN)

Ices of acetylene (C₂ H₂ ) and ammonia (NH₃ ) were irradiated with energetic electrons to simulate interstellar ices processed by galactic cosmic rays in order to investigate the formation of C₂ H₃ N isomers. Supported by quantum chemical calculations, experiments detected product molecules as they sublime from the ices using photoionization reflectron time-of-flight mass spectrometry (PI-ReTOF-MS). Isotopically-labeled ices confirmed the C₂ H₃ N assignments while photon energies of 8.81 eV, 9.80 eV, and 10.49 eV were utilized to discriminate isomers based on their known ionization energies. Results indicate the formation of ethynamine (HCCNH₂ ) and 2H-azirine (c-H₂ CCHN) in the irradiated C₂ H₂ :NH₃ ices, and the energetics of their formation mechanisms are discussed. These findings suggest that these two isomers can form in interstellar ices and, upon sublimation during the hot core phase, could be detected using radio astronomy.

High performance global exploration of isomers and isomerization channels on quantum chemical potential energy surface of H5 C2 NO2

High performance global exploration of isomers and isomerization channels on the quantum chemical potential energy surface (PES) is performed for H₅ C₂ NO₂ by using the scaled hypersphere search-anharmonic downward distortion following (SHS-ADDF) method. A multi-node operation, NeoGRRM, has achieved high performance exploration calculations for the large system by submitting SHS-ADDF sub-jobs into many cores in parallel and unifying the results of sub-jobs into the total lists of the main-job. Global exploration of equilibrium (EQ) and transition-state structures at the level of B3LYP/6-31G(d) gave 3210 EQs and 23278 TSs. Nine compounds were found in the low energy regions of 0-100 kJ/mol; the lowest energy compound is N-methylcarbamic acid, the second is methyl carbamate, and the third is glycine (the most fundamental amino acid). Interconversion pathways between the conformers of each of the low energy compounds were surveyed. Isomerization channels around glycine were explored in detail. The lowest energy barriers around some of the EQs turned to be negative after zero-point energy corrections. This indicates that those structures cannot exist as independent structures because they spontaneously collapse into more stable structures. The global PES search showed various interesting dissociating channels which indicate synthon reaction pathways in the reverse directions.

Infrared Spectra of CH3CN→M, M-η2-(NC)-CH3, CH3-MNC Prepared by Reactions of Laser-Ablated Fe, Ru, and Pt Atoms with Acetonitrile in Excess Argon

Reactions of laser-ablated Fe, Ru, and Pt atoms with acetonitrile have been carried out in excess argon, and the products identified in the matrix spectra. CH₃CN→Fe and Fe-η²-(NC)-CH₃ observed in the original deposition spectra converted to CH₃-FeNC on uv irradiation. CH₃CN→Ru, the only product detected in the Ru system, dissociated on uv irradiation, but was partly reproduced on subsequent visible irradiation and annealing. Similar behavior was found for CH₃CN→Pt. The major products (CH₃-FeNC, CH₃CN→Ru, and CH₃CN→Pt) are the most stable constituents in the previously proposed reaction path for Group 4, 5, 6, and 7 metal atoms and acetonitrile, parallel with the previous results. The Group 8 metal π-coordination products are weakly bound complexes due to limited back-donation to the π*-orbitals of CH₃CN. Calculations show that the Fe insertion product has a much less bent structure than the Ru analogue, in line with its higher s-character from the first row transition-metal to the C-Fe bond, and the group 8 metal methylidenes are not agostically distorted. The Pt to N bond in CH₃CN→Pt is the strongest of all the metals we have investigated owing in large part to its higher electron affinity, which prevents nitrogen lone pair density from entering the pi* orbitals of the C-N group.

Study of Potential Energy Surfaces towards Global Reaction Route Mapping

The potential energy surface (PES) is just a theoretical construct based on the Born-Oppenheimer approximation, but it underlies various phenomena, including molecular vibrations, collisional ionizations, and chemical reactions. This account describes how a new idea for global reaction route mapping (GRRM), which had seemed to be impossible for chemical systems with more than three atoms, was born and has been developed during the course of the study of the PES. GRRM has pioneered new fields of chemistry. Furthermore, techniques for GRRM are still developing, and GRRM is further extending its application to various areas of chemistry and chemical physics.

Potential Energy Surface-Based Automatic Deduction of Conformational Transition Networks and Its Application on Quantum Mechanical Landscapes of d-Glucose Conformers

This paper describes our approach that is built upon the potential energy surface (PES)-based conformational analysis. This approach automatically deduces a conformational transition network, called a conformational reaction route map (r-map), by using the Scaled Hypersphere Search of the Anharmonic Downward Distortion Following method (SHS-ADDF). The PES-based conformational search has been achieved by using large ADDF, which makes it possible to trace only low transition state (TS) barriers while restraining bond lengths and structures with high free energy. It automatically performs sampling the minima and TS structures by simply taking into account the mathematical feature of PES without requiring any a priori specification of variable internal coordinates. An obtained r-map is composed of equilibrium (EQ) conformers connected by reaction routes via TS conformers, where all of the reaction routes are already confirmed during the process of the deduction using the intrinsic reaction coordinate (IRC) method. The postcalculation analysis of the deduced r-map is interactively carried out using the RMapViewer software we have developed. This paper presents computational details of the PES-based conformational analysis and its application to d-glucose. The calculations have been performed for an isolated glucose molecule in the gas phase at the RHF/6-31G level. The obtained conformational r-map for α-d-glucose is composed of 201 EQ and 435 TS conformers and that for β-d-glucose is composed of 202 EQ and 371 TS conformers. For the postcalculation analysis of the conformational r-maps by using the RMapViewer software program we have found multiple minimum energy paths (MEPs) between global minima of ¹C₄ and ⁴C₁ chair conformations. The analysis using RMapViewer allows us to confirm the thermodynamic and kinetic predominance of ⁴C₁ conformations; that is, the potential energy of the global minimum of ⁴C₁ is lower than that of ¹C₄ (thermodynamic predominance) and that the highest energy of those of all the TS structures along a route from ⁴C₁ to ¹C₄ is lower than that of ¹C₄ to ⁴C₁ (kinetic predominance).

An automated transition state search using classical trajectories initialized at multiple minima

TS finding using iterative TSSCDS with trajectories initialized at different minima.

Infrared spectra of CH2=M(H)NC, CH3-MNC, and eta2-M(NC)-CH3 produced by reactions of laser-ablated group 5 metal atoms with acetonitrile

Methylidene isocyanides, methyl isocyanides, and eta(2)-nitrile-pi-complexes are observed in the matrix IR spectra from reactions of Group 5 metals with acetonitrile isotopomers. The primary isocyanide products with no trace of cyanide complexes are consistent with the reaction path proposed in the analogous Zr study. The major products (CH(2)=Ta(H)NC, CH(3)-NbNC, eta(2)-Nb(NC)-CH(3), and eta(2)-V(NC)-CH(3)) after codeposition and reaction of metal with CH(3)CN clearly show the increasing preference for the higher oxidation-state complex on going down the group column, and the subsequent photochemistry provides further information for molecular rearrangements. The Group 5 metal methylidene isocyanides exhibit more agostic distortion than the Zr counterparts and are comparable to the previously studied Group 5 metal methylidene hydrides and halides. The computed structures and observed frequencies indicate that the effects of metal conjugation (C=Ta-N=C:) are minor.

Insight into Global Reaction Mechanism of [C2, H4, O] System from ab Initio Calculations by the Scaled Hypersphere Search Method

A detailed computational study is performed on the singlet potential energy surface (PES) for possible isomerization and dissociation reactions of CH(3)CHO at the DFT (B3LYP/6-311++G(d,p)) and CCSD(T)/cc-pVTZ//B3LYP/6-311++G(d,p) levels. The pathways around the equilibrium structures can be discovered by the scaled hypersphere search (SHS) method, which enables us to make a global analysis of the PES for a given chemical composition. Fourteen isomers inclusive of 11 single-molecules and three "non-stabilized" oxygen-based ylides, 5 energetically favored complexes, and 79 interconversion transition states have been found on the singlet PES. Four lowest lying isomers with thermodynamic stability are also kinetically stable with respect to metastable intermediates. It was revealed that vinyl alcohols, which could be generated by the tautomerization of acetaldehyde, could undergo dissociation to form acetylene and water. In addition, recombination channels between some fragments, such as H(2)CO + (1)CH(2) and (1)CHOH + (1)CH(2), are energetically accessible via collision complex or oxygen-based ylides. Most of available unimolecular decompositions are found to be responsible for favorable hydrogen abstraction processes.

State-Resolved Dynamics of the CN(B2Σ+) and CH(A2Δ) Excited Products Resulting from the VUV Photodissociation of CH3CN

Fourier transform visible spectroscopy, in conjunction with VUV photons produced by a synchrotron, is employed to investigate the photodissociation of CH3CN. Emission is observed from both the CN(B2Sigma+-X2Sigma+) and CH(A2Delta-X2Pi) transitions; only the former is observed in spectra recorded at 10.2 and 11.5 eV, whereas both are detected in the 16 eV spectrum. The rotational and vibrational temperatures of both the CN(B2Sigma+) and CH(A2Delta) radical products are derived using a combination of spectral simulations and Boltzmann plots. The CN(B2Sigma+) fragment displays a bimodal rotational distribution in all cases. Trot(CN(B2Sigma+)) ranges from 375 to 600 K at lower K' and from 1840 to 7700 K at higher K' depending on the photon energy used. Surprisal analyses indicate clear bimodal rotational distributions, suggesting CN(B2Sigma+) is formed via either linear or bent transition states, respectively, depending on the extent of rotational excitation in this fragment. CH(A2Delta) has a single rotational distribution when produced at 16 eV, which results in Trot(CH(A2Delta))=4895+/-140 K in v'=0 and 2590+/-110 K in v'=1. From thermodynamic calculations, it is evident that CH(A2Delta) is produced along with CN(X2Sigma+)+H2. These products can be formed by a two step mechanism (via excited CH3* and ground state CN(X2Sigma+)) or a process similar to the "roaming" atom mechanism; the data obtained here are insufficient to definitively conclude whether either pathway occurs. A comparison of the CH(A2Delta) and CN(B2Sigma+) rotational distributions produced by 16 eV photons allows the ratio between the two excited fragments at this energy to be determined. An expression that considers the rovibrational populations of both band systems results in a CH(A2Delta):CN(B2Sigma+) ratio of (1.2+/-0.1):1 at 16 eV, thereby indicating that production of CH(A2Delta) is significant at 16 eV.

Conversion pathways between a fullerene and a ring among C20 clusters by a sphere contracting walk method: Remarkable difference in local potential energy landscapes around the fullerene and the ring

Conversion pathways from a fullerene to a ring and vice versa among C(20) clusters have been explored. A new technique has been developed for discovering multistep reaction pathways, which can be searched by finding constrained energy minima on spheres whose surfaces are contracting from the starting point to the terminal. Local landscapes that appeared on the pathways were found to be considerably different around the fullerene and the ring, from which one can make a new interpretation for no generation of the C(20) fullerene in laser vaporization of graphite.