Molecular dynamics and metadynamics simulations of [2 + 2] photocycloaddition

Triplet state mechanism of [2 + 2] photocycloaddition forming a cyclobutane ring from two ethylenes is investigated in the context of photocatalysis. High‐level ab initio calculations are combined with ab initio adiabatic molecular dynamics and ab initio metadynamics for rare events modeling. In a photocatalytic scheme, a reactant reaches the triplet state either via intersystem crossing (ISC) or triplet sensitization. The model system adopts a biradical structure, which represents energy intersection with the ground state. The system either completes cyclization or undergoes fragmentation into two olefinic units. The potential and free energy surfaces of the cyclobutane/ethylenes system are mapped with multireference approaches describing possible reaction pathways. To obtain a full picture of a double bond photoreactivity, ab initio adiabatic dynamical calculations were used to estimate reaction yields and to model the effects of excess energy. The potential use of density functional theory based approaches for [2 + 2] photocycloaddition was investigated for future simulations and design of realistic photocatalytic systems. Dynamical aspects of [2 + 2] photocycloaddition via a triplet state manifold are investigated by combining ab initio multireference methods and ab initio molecular dynamics and metadynamics. The reaction pathways are studied for a model system of two ethylenes forming a cyclobutane ring to provide a basis for further studies on design of photocatalytic systems.     

Polystyrene models. II. Ab initio study of isobutylbenzene

The structure and stability of the various conformations of isobutylbenzene are studied using ab initio molecular orbital theory. The calculations show that coupling between the structural units is important. The results indicate that complete geometry optimization of the stable and transition structures of isobutylbenzene produce significant changes in geometrical parameters and charge distributions of this molecule when compared with the corresponding results obtained using the rigid-rotor approximation. These changes are particularly noticeable in one of the gauche conformations and in transition structures of isobutylbenzene generated by the phenyl group rotation. For polystyrene, these results present evidence that there is a strong coupling between the chain-backbone folding and the rotation of the phenyl group. Multidimensional potential energy surfaces are displayed using a topological representation. © 1992 John Wiley & Sons, Inc.     

Theoretical Study on the Long-lived Complexes for the Na＋I2 Collision System

For the Na I2 collision system, theoretical study is performed on the QCISD(T) level by using ab initio method. The ab initio potential energy surfaces are got and on them the long-lived complexes are found and optimized. These results verify the crossed molecule beam experimental phenomenon and the detailed geometry structures are given for the first time. The role of the complexes in the reaction path is also described in detail.     

Transition between icosahedral and cuboctahedral nanoclusters of lead

We have used ab initio methods to study the possible transition between icosahedral (ico) and cuboctahedral (fcc) structures in lead nanoclusters of sizes up to 309 atoms. Spontaneous fcc-to-ico transition in Pb(13) was observed in the ab initio molecular dynamics (MD) simulations at various temperatures. The transition path can be described predominantly by an angular variable s, which can, generally be applied to the similar transitions in clusters of larger sizes and was observed to follow the Mackay model. We have calculated the two-dimensional energy surface that describes the transition in Pb(13) and found a barrierless fcc-to-ico transition path, which is consistent with the observed spontaneous transition in the ab initio MD simulations. The atomic displacements in the transition were identified as one of the vibrational eigenmodes of these two Pb(13) clusters. For clusters of larger sizes (Pb(n), where n = 55, 147, and 309), the possible transitions following similar paths were determined not to be barrierless and the sizes of the barriers were determined by the ab initio elastic band method.     

Nickel and copper doped palladium clusters from a first-principles perspective

A complete study on the evolution of structures and the variation of the energy properties of MPd_n−1 (M = Ni and Cu; n = 2-13) clusters is presented. The study was performed employing auxiliary density functional theory. The obtained results were compared with the results of Pd_n clusters studied with the same methodology. For each cluster size, several structures were studied to determine the lowest energy structures. The initial structures for the geometry optimization were taken along ab initio Born-Oppenheimer molecular dynamics trajectories. Different potentials energy surfaces were studied. All cluster structures were fully optimized without any symmetry restriction. Stable structures, frequencies, spin multiplicities, averaged bond lengths, spin density plots, different energy properties, dipole and magnetic moments as well as charge transfers are reported. This investigation indicates that the palladium clusters doped with a Ni atom are the most stable and potentially the most chemical active ones.     

Ab initio and direct quasiclassical-trajectory study of the F+CH4-->HF+CH3 reaction

We present an electronic structure and dynamics study of the F+CH4-->HF+CH3 reaction. CCSD(T)/aug-cc-pVDZ geometry optimizations, harmonic-frequency, and energy calculations indicate that the potential-energy surface is remarkably isotropic near the transition state. In addition, while the saddle-point F-H-C angle is 180 degrees using MP2 methods, CCSD(T) geometry optimizations predict a bent transition state, with a 153 degrees F-H-C angle. We use these high-quality ab initio data to reparametrize the parameter-model 3 (PM3) semiempirical Hamiltonian so that calculations with the improved Hamiltonian and employing restricted open-shell wave functions agree with the higher accuracy data. Using this specific-reaction-parameter PM3 semiempirical Hamiltonian (SRP-PM3), we investigate the reaction dynamics by propagating quasiclassical trajectories. The results of our calculations using the SRP-PM3 Hamiltonian are compared with experiments and with the estimates of two recently reported potential-energy surfaces. The trajectory calculations using the SRP-PM3 Hamiltonian reproduce quantitatively the measured HF vibrational distributions. The calculations also agree with the experimental HF rotational distributions and capture the essential features of the excitation function. The results of the SRP semiempirical Hamiltonian developed here clearly improve over those using the two prior potential-energy surfaces and suggest that reparametrization of semiempirical Hamiltonians is a promising strategy to develop accurate potential-energy surfaces for reaction dynamics studies of polyatomic systems.     

Potential energy surfaces of an adenine-thymine base pair and its methylated analogue in the presence of one and two water molecules: molecular mechanics and correlated ab initio study

Potential energy surfaces of monohydrated and dihydrated adenine-thymine and 9-methyladenine-1-methylthymine base pairs were examined by the molecular dynamics/quenching technique using the Cornell et al. force field (J. Am. Chem. Soc. 1995, 117, 5179). Long runs of molecular dynamics/quenching calculations allowed us to evaluate the free energy surface. The most stable and populated structures found were fully reoptimized at the correlated ab initio level employing the resolution of identity M?ller-Plesset method. A systematic study of the base pairs' microhydration using both the empirical and the high-level correlated ab initio approaches is presented for the first time. We show that the occurrence of water molecules and their gradually increasing number as well as the methylation of the bases favor stacked structures over the planar hydrogen-bonded ones. These results based on the correlated ab initio calculations are in the excellent agreement with data obtained from our previous empirical potential molecular dynamics study (Kabelác et al. Chem.-Eur. J. 2001, 7, 2067).     

Predissociation via conformational change: photodissociation of N,N-dimethylnitrosamine in the S(1) state

We investigated the photodissociation mechanism of N,N-dimethylnitrosamine (CH(3))(2)NNO (DMN) by ab intio quantum chemical methods. Inspired by an earlier study we calculated two-dimensional potential energy surfaces of the S(1) state of DMN in its planar and pyramidal conformations. While the planar molecular geometry appears to possess no direct dissociation channel, the pyramidal configuration is dissociative yielding the products NO + (CH(3))(2)N. Using wave packet dynamics on the planar S(1) potential energy surface the experimental absorption spectrum was well reproduced which gives indirect but strong support for the nondissociative nature of this surface. The transition from the planar to the pyramidal conformation of DMN was then investigated by an ab initio molecular dynamics method which revealed the time evolution of the geometrical parameters of the molecule up to the dissociation of the N-N bond. This occurs about 90 fs after photon excitation. The calculated minimum energy path along the N-N coordinate and the structural changes of the molecule along this coordinate provided a detailed picture of this indirect dissociation or, more specific, predissociation process via conformational change.     

Reaction path potential for complex systems derived from combined ab initio quantum mechanical and molecular mechanical calculations

Combined ab initio quantum mechanical and molecular mechanical calculations have been widely used for modeling chemical reactions in complex systems such as enzymes, with most applications being based on the determination of a minimum energy path connecting the reactant through the transition state to the product in the enzyme environment. However, statistical mechanics sampling and reaction dynamics calculations with a combined ab initio quantum mechanical (QM) and molecular mechanical (MM) potential are still not feasible because of the computational costs associated mainly with the ab initio quantum mechanical calculations for the QM subsystem. To address this issue, a reaction path potential energy surface is developed here for statistical mechanics and dynamics simulation of chemical reactions in enzymes and other complex systems. The reaction path potential follows the ideas from the reaction path Hamiltonian of Miller, Handy and Adams for gas phase chemical reactions but is designed specifically for large systems that are described with combined ab initio quantum mechanical and molecular mechanical methods. The reaction path potential is an analytical energy expression of the combined quantum mechanical and molecular mechanical potential energy along the minimum energy path. An expansion around the minimum energy path is made in both the nuclear and the electronic degrees of freedom for the QM subsystem internal energy, while the energy of the subsystem described with MM remains unchanged from that in the combined quantum mechanical and molecular mechanical expression and the electrostatic interaction between the QM and MM subsystems is described as the interaction of the MM charges with the QM charges. The QM charges are polarizable in response to the changes in both the MM and the QM degrees of freedom through a new response kernel developed in the present work. The input data for constructing the reaction path potential are energies, vibrational frequencies, and electron density response properties of the QM subsystem along the minimum energy path, all of which can be obtained from the combined quantum mechanical and molecular mechanical calculations. Once constructed, it costs much less for its evaluation. Thus, the reaction path potential provides a potential energy surface for rigorous statistical mechanics and reaction dynamics calculations of complex systems. As an example, the method is applied to the statistical mechanical calculations for the potential of mean force of the chemical reaction in triosephosphate isomerase.     

羟基负离子与一氟甲烷双分子亲核取代反应的量子化学研究

本文用LCAO-MO-SCF ab initio方法, 对OH^-+CH3F→CH3OH+F^-反应进行了过渡态理论及前线轨道理论的量子化学研究, 以4-31G为基组, 计算了反应进程的势能曲线,得到了过渡态的几何构型, 并用MP2方法进行了电子库仑相关效应的校正, 反应活化能的计算值与实验数据较为一致。用前线轨道理论对反应中分子重新组合过程进行了轨道分析, 较全面地解释了该反应的机理。     

Modeling water exchange on an aluminum polyoxocation

For the first time, water exchange on a polymeric complex has been modeled using a combination of gas-phase ab initio calculations and molecular dynamics (MD) simulations. The GaO4Al12(OH)24(H2O)12(7+)aq ion (GaAl12) was chosen because high-quality experimental data exist, including an activation enthalpy (+63 +/- 7 kJ/mol) and an activation volume (+3 +/- 1 cm3/mol). We took a two-step approach. First, the local solvent structure and the initial states for reaction were inferred from the molecular dynamics simulations. Second, we used this information to evaluate initial-state structures in the ab initio calculations. The energy differences between the initial and transition states from the ab initio calculations varied from +59 kJ/mol to +53 kJ/mol depending upon details, closely approximating the activation enthalpy.     

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

Technical details of a new global mapping technique for finding equilibrium (EQ) and transition structures (TS) on potential energy surfaces (PES), the scaled hypersphere search (SHS) method (Ohno, K.; Maeda, S. Chem. Phys. Lett. 2004, 384, 277), are presented. On the basis of a simple principle that reaction pathways are found as anharmonic downward distortions of PES around an EQ point, the reaction pathways can be obtained as energy minima on the scaled hypersphere surface, which would have a constant energy when the potentials are harmonic. Connections of SHS paths between each EQ are very similar to corresponding intrinsic reaction coordinate (IRC) connections. The energy maximum along the SHS path reaches a region in close proximity to the TS of the reaction pathway, and the subsequent geometry optimization from the SHS maximum structure easily converges to the TS. The SHS method, using the one-after-another algorithm connecting EQ and TS, considerably reduces the multidimensional space to be searched to certain limited regions around the pathways connecting each EQ with the neighboring TS. Applications of the SHS method have been made to ab initio surfaces of formaldehyde and propyne molecules to obtain systematically five EQ and nine TS for formaldehyde and seven EQ and 32 TS for propyne.     

Multiple time scale Hartree–Fock molecular dynamics

This is the first application of a rigorous, established multiple time-step method to ab initio molecular dynamics. The resulting algorithm is conceptually simple and easy to implement, but very effective. It translates the large mass differences present in ab initio molecular dynamics into substantial savings in computer time while retaining high accuracy. This transforms ab initio molecular dynamics from a desirable but prohibitively expensive possibility into a viable method, at least for short-time phenomena in small systems or for otherwise inaccessibly complicated potential energy surfaces. © 1993 John Wiley & Sons, Inc.     

A test of the new “integrated MO + MM” (IMOMM) method for the conformational energy of ethane and n-butane

Test calculations of the newly developed “Integrated Molecular Orbital + Molecular Mechanics” (IMOMM) method were performed for the optimized equilibrium and transition structures and energies of ethane and n-butane. In this method, the total energy of a large molecular system is expressed as a sum of the MO energy of the small “model” system and a modified MM energy of the “real” system, and full geometry optimization is carried out using the gradient of this total energy. Various schemes of partition of the system into the MO part and the MM part, including some not intended in the original design of the method, were examined and compared with the pure ab initio MO and the pure MM results. In most reasonable partition schemes, the IMOMM method can reproduce the pure ab initio and the pure MM geometries and energies quite well. © 1996 John Wiley & Sons, Inc.     

Calculating geochemical reaction pathways--exploration of the inner-sphere water exchange mechanism in Al(H2O)6(3+)(aq) + nH2O with ab Initio calculations and molecular dynamics

We have simulated exchange of inner-sphere and bulk water molecules for different sizes of Al3+(aq) clusters, Al(H2O)63+ + nH2O for n = 0, 1, 6, or 12, with ab initio and molecular dynamics simulations, in order to understand how robust the ab initio method is for identifying hydrolytic reaction pathways of particular importance to geochemistry. In contrast to many interfacial reactions, this particular elementary reaction is particularly simple and well-constrained by experiment. Nevertheless, we find that a rich array of parallel reaction pathways depend sensitively on the details of the solvation sphere and structure and that larger clusters are not necessarily better. Inner-sphere water exchange in Al3+(aq) may occur through two Langford-Gray dissociative pathways, one in which the incoming and outgoing waters are cis, the other in which they are trans to one another. A large majority of exchanges in the molecular dynamics simulations occurred via the trans mechanism, in contrast to the predictions of the ab initio method. In Al(H2O)63+ + H2O, the cis mechanism has a transition state of 84.3 kJ/mol, which is in good agreement with previous experimental and ab initio results, while the trans mechanism has only a saddle point with two negative frequencies, not a transition state, at 89.7 kJ/mol. In addition to the exchange mechanisms, dissociation pathways could be identified that were considerably lower in energy than experiment and varied considerably between 60 and 100 kJ/mol, depending on the particular geometry and cluster size, with no clear relation between the two. Ab initio calculations using large clusters with full second coordination spheres (n = 12) were unable to find dissociation or exchange transition states because the network of hydrogen bonds in the second coordination sphere was too rigid to accommodate the outgoing inner-sphere water. Our results indicate that caution should surround ab initio simulation of complicated dynamic processes such as hydrolysis, ion exchange, and interfacial reactions that involve several steps. Dynamic methods of simulation need to accompany static methods such as ab initio calculation, and it is best to consider simulated pathways as hypotheses to be tested experimentally rather than definitive properties of the reaction.     

Smooth solvation method for d-orbital semiempirical calculations of biological reactions. 1. Implementation

The present paper describes the extension of a recently developed smooth conductor-like screening model for solvation to a d-orbital semiempirical framework (MNDO/d-SCOSMO) with analytic gradients that can be used for geometry optimizations, transition state searches, and molecular dynamics simulations. The methodology is tested on the potential energy surfaces for separating ions and the dissociative phosphoryl transfer mechanism of methyl phosphate. The convergence behavior of the smooth COSMO method with respect to discretization level is examined and the numerical stability of the energy and gradient are compared to that from conventional COSMO calculations. The present method is further tested in applications to energy minimum and transition state geometry optimizations of neutral and charged metaphosphates, phosphates, and phosphoranes that are models for stationary points in transphosphorylation reaction pathways of enzymes and ribozymes. The results indicate that the smooth COSMO method greatly enhances the stability of quantum mechanical geometry optimization and transition state search calculations that would routinely fail with conventional solvation methods. The present MNDO/d-SCOSMO method has considerable computational advantages over hybrid quantum mechanical/molecular mechanical methods with explicit solvation, and represents a potentially useful tool in the arsenal of multi-scale quantum models used to study biochemical reactions.     

Derivation of class II force fields. VIII. Derivation of a general quantum mechanical force field for organic compounds

A class II valence force field covering a broad range of organic molecules has been derived employing ab initio quantum mechanical "observables." The procedure includes selecting representative molecules and molecular structures, and systematically sampling their energy surfaces as described by energies and energy first and second derivatives with respect to molecular deformations. In this article the procedure for fitting the force field parameters to these energies and energy derivatives is briefly reviewed. The application of the methodology to the derivation of a class II quantum mechanical force field (QMFF) for 32 organic functional groups is then described. A training set of 400 molecules spanning the 32 functional groups was used to parameterize the force field. The molecular families comprising the functional groups and, within each family, the torsional angles used to sample different conformers, are described. The number of stationary points (equilibria and transition states) for these molecules is given for each functional group. This set contains 1324 stationary structures, with 718 minimum energy structures and 606 transition states. The quality of the fit to the quantum data is gauged based on the deviations between the ab initio and force field energies and energy derivatives. The accuracy with which the QMFF reproduces the ab initio molecular bond lengths, bond angles, torsional angles, vibrational frequencies, and conformational energies is then given for each functional group. Consistently good accuracy is found for these computed properties for the various types of molecules. This demonstrates that the methodology is broadly applicable for the derivation of force field parameters across widely differing types of molecular structures. Copyright 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1782-1800, 2001     

Computational strategies for evaluating barrier heights for gas-phase reactions of lithium enolates

Gas-phase activation energies were calculated for three lithium enolate reactions by using several different ab initio and density functional theory (DFT) methods to determine which levels of theory generate acceptable results. The reactions included an aldol-type addition of an enolate to an aldehyde, a proton transfer from an alcohol to a lithium enolate, and an S(N)2 reaction of an enolate with chloromethane. For each reaction, the calculations were performed for both the monomeric and dimeric forms of the lithium enolate. It was found that transition state geometry optimization with B3LYP followed by single point MP2 calculations generally provided acceptable results compared to higher level ab initio methods.     

Photodissociation of acetic acid in the gas phase: an ab initio study

Photodissociation of acetic acid in the gas phase was investigated using ab initio molecular orbital methods. The stationary structures on the ground-state potential energy surfaces were mainly optimized at the MP2 level of theory, while those on the excited-state surfaces were determined by complete active space SCF calculations with a correlation-consistent basis set of cc-pVDZ. The reaction pathways leading to different photoproducts are characterized on the basis of the computed potential energy surfaces and surface crossing points. The calculations reproduce the experimental results well and provide additional insight into the mechanism of the ultraviolet photodissociation of acetic acid and related compounds.