Global Mapping of Equilibrium and Transition Structures on Potential Energy Surfaces by the Scaled Hypersphere Search Method: Applications to ab Initio Surfaces of Formaldehyde and Propyne Molecules

Quantum Chemistry Study of H+(H2O)8: A Global Search for Its Isomers by the Scaled Hypersphere Search Method, and Its Thermal Behavior

The structures of the protonated water cluster H+(H2O)8 have been globally explored by the scaled hypersphere search method. On the Hartree-Fock potential energy surface 174 isomers were found, among which 168 were computed to be minima at the B3LYP/6-31+G** level, and their energies were further refined at the level of MP2/6-311++G(3df,2p). The global minimum on the potential energy surface computed at the B3LYP/6-31+G** level shows a cagelike structure with the "Eigen" motif, while the lowest-free-energy isomer has a five-membered-ring structure at 170 K and a chain form at 273 K. The present results are well in line with previous experimental findings. In addition, the ADMP (atom-centered density matrix propagation) simulation indicates that the extra proton in the lowest-free-energy isomer (170 K), which has a five-membered ring and the "Zundel" feature, is often in an asymmetrical hydrogen bond.

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.

Structures of Water Octamers (H2O)8: Exploration on Ab Initio Potential Energy Surfaces by the Scaled Hypersphere Search Method

The potential energy surface (PES) of water octamers has been explored by the scaled hypersphere search method. Among 164 minima on the PES (based on MP2/6-311++G(3df,2p)//B3LYP/6-311+G(d,p) calculations), the cubic structure with D2d symmetry has been confirmed to be the global minimum. In a thermodynamic simulation using these 164 structures, the cubic structure with S4 symmetry has the highest population at low temperature, though double rings can become dominant as temperature going up, in good accord with a recent Monte Carlo simulation using an empirical potential. A transition temperature from cubic to noncubic has significantly been underestimated when potential energy data of B3LYP/6-311+G(d,p) calculations are employed in the simulation. This serious discrepancy between the MP2 and the B3LYP results suggests an importance of dispersion interactions for discussions on thermodynamics of water octamers.

Effect of confinement by porous carbons on the unimolecular decomposition of formaldehyde

As part of an effort to understand the effect of confinement by porous carbons on chemical reactions, we have carried out density functional theory calculations on the unimolecular decomposition of formaldehyde within graphitic carbons. Our results show that the interactions with the carbon walls result in a lowering of the reaction barrier. For larger pores, there is also a shift of the equilibrium towards the formation of carbon monoxide and hydrogen at low temperatures. This trend is reversed for small pore sizes.

Global Reaction Route Mapping on Potential Energy Surfaces of Formaldehyde, Formic Acid, and Their Metal-Substituted Analogues

Global reaction route mapping of equilibrium structures, transition structures, and their connections on potential energy surface (PES) has been done for MCHO (M = H, Li, Na, Al, Cu) and HCO2M (M = H, Li). A one-after-another technique based on the scaled hypersphere search method has been successfully applied to exploring unknown chemical structures, transition structures, and reaction pathways for organometallic systems. Upon metal substitution, considerable changes of stable structures, reaction pathways, and relative heights of transition structures have been discovered, though some features are similar among the analogues. Al and Cu atoms were found to behave as very strong scissors to cut the CO double bond in MCHO. Energy profiles of the CO insertion into Li-H and Li-CH3 bonds were found to be very similar, especially around the structures where the Li atom is not directly connected with the methyl group, which indicates little effects of alkyl substitution on the reaction route topology.

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.

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

An extensive quantum chemical study of the potential energy surface (PES) for all possible isomerization and dissociation reactions of CH3CN is reported at the DFT (B3LYP/6-311++G(d,p)) and CCSD(T)/ cc-pVTZ//B3LYP/6-311++G(d,p) levels of theory. 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 potential energy surface for a given chemical composition in combination with a downhill-walk algorithm. Seventeen equilibrium structures and 59 interconversion transition states have been found on the singlet PES. The four lowest lying isomers with thermodynamic stability are also kinetically stable with the lowest conversion barriers of 49.69-101.53 kcal/mol at the CCSD(T)/cc-pVTZ//B3LYP/6-311++G(d,p) level, whereas three-membered-ring isomers c-CH2NCH, c-CH2CNH, and c-CHNHCH can be considered as metastable intermediates which can further convert into the low-lying chain-like isomers and higher lying acyclic isomers with the lowest conversion energies of 21.70-59.99 kcal/mol. Thirteen available dissociation channels depending on the different initial isomers have been identified. A prediction can be made for the possible mechanism explaining the migration of a hydrogen atom in competition with the CC bond dissociation. Several new energetically accessible pathways are found to be responsible for the migration of the hydrogen atom. The present results demonstrate that the SHS method is an efficient and powerful technique for global mapping of reaction pathways on PESs.