Phase space geometry of dynamics passing through saddle coupled with spatial rotation

Nonlinear reaction dynamics through a rank-one saddle is investigated for many-particle system with spatial rotation. Based on the recently developed theories of the phase space geometry in the saddle region, we present a theoretical framework to incorporate the spatial rotation which is dynamically coupled with the internal vibrational motions through centrifugal and Coriolis interactions. As an illustrative simple example, we apply it to isomerization reaction of HCN with some nonzero total angular momenta. It is found that no-return transition state (TS) and a set of impenetrable reaction boundaries to separate the "past" and "future" of trajectories can be identified analytically under rovibrational couplings. The three components of the angular momentum are found to have distinct effects on the migration of the "anchor" of the TS and the reaction boundaries through rovibrational couplings and anharmonicities in vibrational degrees of freedom. This method provides new insights in understanding the origin of a wide class of reactions with nonzero angular momentum.     

One-dimensional quantum description of the bending vibrations of HCN/CNH for high values of the total angular momentum

For high values of the quantum number of the total angular momentum J (up to J = 20), quantum mechanical eigenstates (eigenvalues and eigenfunctions) are calculated by the method of Gatti et al. (J. Mol. Spectrosc. 181 (1997) 403) for the bending deformations of HCN and CNH. In particular, we have examined the l-type resonances in highly excited rovibrational states within the framework of a one-dimensional model, i.e. along the reaction pathway for the isomerization reaction HCN/CNH. The potential energy surface used is that of Bowman et al. (J. Chem. Phys. 99 (1993) 308).     

Insights into the Thermal and Photochemical Reaction Mechanisms of Azidoacetonitrile. Spectroscopic and MS-CASPT2 Calculations

This work studies the photochemical and thermal decompositions of azidoacetonitrile (N₃CH₂CN) from both the experimental and theoretical points of view. The data of the photochemical experiments are taken from the literature, while the thermal decomposition have been carried out by us. In addition, we have performed ab initio calculations of the multiconfigurational type [complete active space self-consistent field (CASSCF) and the multistate multireference perturbation theory (MS-CASPT2)]. It is found that the first step of both type of decompositions is N₂ elimination and formation of closed shell singlet nitrene. Afterwards, the nitrene tends to rapidly rearrange into formimidoyl cyanide (HNCHCN). As both reactions progress, the imine isomerizes into formimidoyl isocyanide (HNCHNC). The photoisomerization of the imine takes places thorough a conical intersection, while the same reaction on the ground electronic state occurs via a conventional transition state. The last step of the global reaction is decomposition of the imines into HCN and CNH. In photochemical conditions, the conjunction of the imines and its dissociation products (HCN and CNH) yields adenine     

The role of long-range forces in the determination of translational kinetic energy release. Loss of C4H4+ from the benzene and pyridine cations

Kinetic energy release distributions (KERDs) for the benzene ion fragmenting into C 4H 4 (+) and C 2H 2 have been recorded by double-focusing mass spectrometry in the metastable energy window and by a retarding field experiment up to an energy of 5 eV above the fragmentation threshold. They are compared with those resulting from the HCN loss reaction from the pyridine ion. Both reactions display a similar variation of the kinetic energy release as a function of the internal energy: the average release is smaller than statistically expected, with a further restriction of the phase-space sampling for the C 5H 5N (+) dissociation. Ab initio calculations of the potential-energy profile have been carried out. They reveal a complicated reaction mechanism, the last step of which consists in the dissociation of a weakly bound ion-quadrupole or ion-dipole complex. The KERDs have been analyzed by the maximum entropy method. The fraction of phase space effectively sampled by the pair of fragments has been determined and is similar for both dissociations. Both reactions are constrained by the square root of the released translational energy, epsilon (1/2). This indicates that in the latter stage of the dissociation process, the reaction coordinate is adiabatically decoupled from the bath of the bound degrees of freedom. For the C 6H 6 (+) fragmentation, the analysis of the experimental results strongly suggests that, just as for a spherically symmetric interaction potential, the translational motion is confined to a plane. For the dissociation of the pyridine ion, the main dynamical constraint is also a restriction to a two-dimensional subspace. This dimensionality reduction of the translational phase space is due to the fact that the Hamiltonian of both weakly bound complexes contains a cyclic coordinate.     

Controlled subnanosecond isomerization of HCN to CNH in solution

We report a study of control of the HCN-->CNH isomerization in a liquid Ar solution. We show, using molecular dynamics simulations, nearly complete conversion from HCN to CNH can be achieved in solution on the subnanosecond time scale without requiring laser pulse shaping or molecular alignment. The mechanism of the isomerization reaction involves multiphoton rovibrational excitation on the ground electronic state potential energy surface coupled with rapid rovibrational relaxation in solution. The results demonstrate the important role of rotation-vibration coupling in multiphoton excitation of small molecules and constitute the first realistic computational demonstration of fast, robust, and high-yield laser field manipulation of solution-phase molecular processes.     

Collision-induced dissociation of HS-(HCN): unsymmetrical hydrogen bonding in a proton-bound dimer anion

The energy-resolved competitive collision-induced dissociation of the proton-bound complex [HS.H.CN](-) is studied in a guided ion beam tandem mass spectrometer. H(2)S and HCN have nearly identical gas-phase acidities, and therefore, the HS(-) + HCN and the CN(-) + H(2)S product channels exhibit nearly the same threshold energies, as expected. However, the HS(-) + HCN channel has a cross section up to a factor of 50 larger than CN(-) + H(2)S at higher energies. The cross sections are modeled using RRKM theory and phase space theory. The complex dissociates to HS(-)+ HCN via a loose transition state, and it dissociates to CN(-) + H(2)S via a tight transition state. Theoretical calculations show that the proton-transfer potential energy surface has a single minimum and that the hydrogen bonding in the complex is strongly unsymmetrical, with an ion-molecule complex of the form HS(-)..HCN rather than CN(-)..H(2)S or an intermediate structure. The requirement for proton transfer before dissociation and curvature along the reaction path impedes the CN(-) + H(2)S product channel.     

Nonempirical statistical theory for atomic evaporation from nonrigid clusters: applications to the absolute rate constant and kinetic energy release

A high energy atomic cluster undergoing frequent structural isomerization behaves like a liquid droplet, from which atoms or molecules can be emitted. Even after evaporation, the daughter cluster may still keep changing its structure. We study the dynamics of such an evaporation process of atomic evaporation. To do so, we develop a statistical rate theory for dissociation of highly nonrigid molecules and propose a simple method to calculate the absolute value of classical phase-space volume for a potential function that has many locally stable basins. The statistical prediction of the final distribution of the released kinetic energy is also developed. A direct application of the Rice-Ramsperger-Kassed-Marcus (RRKM) theory to this kind of multichannel chemical reaction is prohibitively difficult, unless further modeling and/or assumptions are made. We carry out a completely nonempirical statistical calculation for these dynamical quantities, in that nothing empirical is introduced like remodeling (or reparametrization) of artificial potential energy functions or recalibration of the phase-space volume referring to other "empirical" values such as those estimated with the molecular dynamics method. The so-called dividing surface is determined variationally, at which the flux is calculated in a consistent manner with the estimate of the phase-space volume in the initial state. Also, for the correct treatment of a highly nonrigid cluster, the phase-space volume and flux are estimated without the separation of vibrational and rotational motions. Both the microcanonical reaction rate and the final kinetic energy distribution thus obtained have quite accurately reproduced the corresponding quantities given by molecular dynamics calculations. This establishes the validity of the statistical arguments, which in turn brings about the deeper physical insight about the evaporation dynamics.     

Communication: transition state theory for dissipative systems without a dividing surface

Transition state theory is a central cornerstone in reaction dynamics. Its key step is the identification of a dividing surface that is crossed only once by all reactive trajectories. This assumption is often badly violated, especially when the reactive system is coupled to an environment. The calculations made in this way then overestimate the reaction rate and the results depend critically on the choice of the dividing surface. In this Communication, we study the phase space of a stochastically driven system close to an energetic barrier in order to identify the geometric structure unambiguously determining the reactive trajectories, which is then incorporated in a simple rate formula for reactions in condensed phase that is both independent of the dividing surface and exact.     

Is N-protonated hydrogen isocyanide, H2NC+, an observable interstellar species?

Ab initio molecular orbital theory is used to examine the singlet and triplet potential energy surfaces for the CH2N+ system. The results confirm those of earlier studies which suggested that the singlet H2NC+ isomer could be formed via the corresponding triplet isomer. Also, it is shown that the reaction HCN+ + H2 might lead to this metastable isomer without invoking the triplet species. The best test of the hypothesis that this molecule can be formed by gas phase, ion molecule reactions and may be an important precursor in the interstellar synthesis of HCN and HNC is to search for it in space. To this end, theoretical predictions are made of its rotational frequencies and its vibrational frequencies and intensities to serve as a guide to laboratory spectroscopists and radioastronomers.     

Ni+ reactions with aminoacrylonitrile, a species of potential astrochemical relevance

The reaction of aminoacrylonitrile, a species of astrochemical interest, with Ni(+)((2)D(5/2)) was investigated by means of mass spectrometry techniques and density functional theory calculations. The dominant fragmentations in the MIKE spectrum correspond to the loss of [C2,N,H3], HCN, and NH3, the loss of H2 being very minor. The structure and bonding of the different aminoacrylonitrile-Ni(+) complexes were investigated at the B3LYP/6-311G(d,p) level of theory. The same approach was employed in our survey of the corresponding potential energy surface. This survey indicates that the [C2,N,H3] neutral product can be formed either as ketenimine (CH2CNH) or acetonitrile. The formation of the latter is significantly more exothermic but involves slightly higher activation barriers; so very likely, both isomers are produced along the reaction process. The lost of HNC is not competitive with the loss of HCN, because when the former is formed the products lie higher in energy and the corresponding mechanisms involve energy barriers above the entrance channel. The loss of NH3 is associated with the formation of a complex between cyanoacetylene, HCCCN, which is very abundant in the interstellar media, and Ni(+).     

Formation of surface CN from the coupling of C and N atoms on Pt(111)

Reflection absorption infrared spectroscopy (RAIRS) and temperature-programmed desorption have been used to study the coupling of C and N atoms on Pt(111) to form surface CN. This reaction underlies the important synthesis of HCN from methane and ammonia over platinum catalysts. Since CH4 and NH3 do not thermally dissociate on Pt(111) under ultrahigh vacuum conditions, we used CH3I and electron bombardment of NH3 to generate reactive surface species. Surface CN is formed at a temperature of 500 K from surface Nads and Cads atoms. The presence of surface CN is detected by HCN desorption and through the reaction of hydrogen with CNads to form a surface >CNH2 (aminocarbyne) species, which has a characteristic RAIR spectrum.     

Ni(+) reactions with aminoacetonitrile, a potential prebiological precursor of glycine

The gas-phase reactions between Ni(+) ((2)D(5/2)) and aminoacetonitrile, a molecule of prebiological interest as possible precursor of glycine, have been investigated by means of mass spectrometry techniques. The mass-analyzed ion kinetic energy (MIKE) spectrum reveals that the adduct ions [NC--CH(2)--NH(2), Ni(+)] spontaneously decompose by loosing HCN, H(2), and H(2)CNH, the loss of hydrogen cyanide being clearly dominant. The structures and bonding characteristics of the aminoacetonitrile-Ni(+) complexes as well as the different stationary points of the corresponding potential energy surface (PES) have been theoretically studied by density functional theory (DFT) calculations carried out at B3LYP/6-311G(d,p) level. A cyclic intermediate, in which Ni(+) is bisligated to the cyano and the amino group, plays an important role in the unimolecular reactivity of these ions, because it is the precursor for the observed losses of HCN and H(2)CNH. In all mechanisms associated with the loss of H(2), the metal acts as hydrogen carrier favoring the formation of the H(2) molecule. The estimated bond dissociation energy of aminoacetonitrile-Ni(+) complexes (291 kJ mol(-1)) is larger than those measured for other nitrogen bases such as pyridine or pyrimidine and only slightly smaller than that of adenine.     

HCN Adsorption on the Fe(100), Fe(111) and Fe(110) Surfaces: a Density Functional Theory Study

Twenty kinds of adsorptions of HCN on the Fe(100), Fe(111) and Fe(110) surfaces at the 1/4 monolayer coverage are found using the density functional theory. For Fe(100), the adsorption energy of the most stable configuration where the HCN locates at the fourfold site with the C-N bonded to four Fe atoms is 1.928 eV. The most favored adsorption structure for HCN on Fe(111) is f-η3(N)-h-η3(C), in which the C-N bond is almost parallel to the surface, and the adsorption energy is 1.347 eV. On Fe(110), the adsorption energy in the most stable configuration in which HCN locates at the two long-bridge sites is 1.777 eV. The adsorption energy of the parallel orientation for HCN is larger than that of the perpendicular configuration. The binding mechanism of HCN on the Fe(100), Fe(111) and Fe(110) surfaces is also analyzed by Mulliken charge population and the density of states in HCN. The result indicates that the configurations in which the adsorbed HCN becomes the non-linear are beneficial to the formation of the addition reaction for hydrogen. The nature that the introduction of Fe into the catalyst could increase the catalytic activity of the bimetallic catalyst in the addition reaction of hydrogen for nitriles is revealed.     

C-N coupling on transition metal surfaces: a density functional theory study

We have investigated the formation of C-N bonds from individual atoms and single hydrogenated moieties on a series of transition metals. These reactions play a role in HCN formation at high oxygen coverage, also known as Andrussow oxidation, and they are fundamental to understand the ability of other materials to form part of alloys where Pt is the major component. Dehydrogenations take place quite easily under these high oxygen conditions and thus, the C+N, HC+N, and N+CH recombinations to form HCN or its isomer CNH might represent the rate-limiting steps for the reaction. For all the metals in the present study we have found that the activation energy for the reactions between H(x)C and NH(y) (x,y = 0,1) involved in C-N formation follow a linear relationship with the adsorption energy of the N atom. This is due to the common nature of all these transition states, where N-containing fragments get activated from three-fold hollow sites to bridge positions. The slopes of the linear dependence, though, depend on the valence of the N fragment, i.e., smaller slopes are found for NH moieties with respect to N ones.     

Theoretical investigation of the dissociation dynamics of vibrationally excited vinyl bromide on an ab initio potential-energy surface obtained using modified novelty sampling and feed-forward neural networks

The reaction dynamics of vibrationally excited vinyl bromide have been investigated using classical trajectory methods on a neural network potential surface that is fitted to an ab initio database of 12 122 configuration energies obtained from electronic structure calculations conducted at the MP4(SDQ) level of theory using a 6-31G(d,p) basis set for the carbon and hydrogen atoms and Huzinaga's (43334334) basis set augmented with split outer s and p orbitals (4332143214) and a polarization f orbital with an exponent of 0.5 for the bromine atom. The sampling of the 12-dimensional configuration hyperspace of vinyl bromide prior to execution of the electronic structure calculations is accomplished by combining novelty-sampling methods, chemical intuition, and trajectory sampling on empirical and neural network surfaces. The final potential is obtained using a two-layer feed-forward neural network comprising 38 and 1 neurons, respectively, with hyperbolic tangent sigmoid and linear transfer functions in the hidden and output layers, respectively. The fitting is accomplished using the Levenberg-Marquardt algorithm with early stopping and Bayesian regularization methods to avoid overfitting. The interpolated potentials have a standard deviation from the ab initio results of 0.0578 eV, which is within the range generally regarded as "chemical accuracy" for the purposes of electronic structure calculations. It is shown that the potential surface may be easily and conveniently transferred from one research group to another. The files required for transfer of the vinyl bromide surface can be obtained from the Electronic Physics Auxiliary Publication Service. Total dissociation rate coefficients for vinyl bromide are obtained at five different excitation energies between 4.50 and 6.44 eV. Branching ratios into each of the six open reaction channels are computed at 24 vibrational energies in the range between 4.00 and 6.44 eV. The distribution of vibrational energies in HBr formed via three-center dissociation from vinyl bromide is determined and compared with previous theoretical and experimental results. It is concluded that the combination of ab initio electronic structure calculations, novelty sampling with chemical intuition and trajectories on empirical analytic surfaces, and feed-forward neural networks provides a viable framework in which to execute purely ab initio molecular-dynamics studies on complex systems with multiple open reaction channels.     

Synthesis and modification of carbon nanohorns structure for hyperthermic application

Single-wall carbon carcass nanoparticles (nanohorns) (CNHs) are synthesized by two methods: electric arc and electron evaporation of graphite in the inert atmosphere. Distinctions in the structures of obtained materials are revealed using electron microscopy and Raman spectroscopy. The chemical structure of the CNH surface is investigated by X-ray photoelectron and NEXAFS spectroscopy during oxidation. It is found that oxidation causes the destruction of CNH agglomerates and weakly affects the structure of graphene nets. However, these changes are sufficient for an increase in the infrared radiation absorption by the dispersion of nanohorns in water. It is shown that the efficiency of 808 nm laser heating of a CNH dispersion depends on the synthesis method and chemical modification of nanoparticles, which enables their potential use for local hyperthermia of cells of living organisms in cancer therapy.     

The topology of energy hypersurfaces II. Reaction topology in euclidean spaces

By introducing equivalence classes of critical points of potential energy hypersurfaces a unique topological space, Reaction Topology (^3N E, T _C ) is defined over an Euclidean nuclear configuration space ^3N E for a system of N nuclei and k electrons. Relations between the topological concepts of molecular structure and reaction mechanism are analyzed. Topological equivalences between Euclidean and Riemannian representations of nuclear configuration spaces are exploited for the analysis of quantum mechanical reaction networks of all possible chemical reactions over the given potential energy hypersurface.