Search

We have investigated the scope and limitations of a method for predicting the regioisomer distribution in electrophilic aromatic substitution reactions that are under kinetic control. This method is based on calculation of the relative stabilities of the sigma-complex intermediates using density functional theory. Predictions from this method can be used quantitatively for halogenations; it agreed to an accuracy of about 1 kcal/mol with experimental observations in 10 of the 11 investigated halogenation reactions. For nitrations, the method gave useful predictions for heterocyclic substrates. The method failed for nitration of monosubstituted benzenes, and we expect that more elaborate model systems, including explicit solvent molecules, will be necessary to obtain quantitatively useful predictions for such cases. For Lewis acid promoted Friedel-Crafts acylations, the method can be expected to give qualitatively correct predictions, that is, to point out the dominating isomer. For substrates where the regioisomeric outcome is highly dependent on the reaction conditions, the method can only be of qualitative use if the concentration of the free Lewis acid is high during the reaction. We have also compared the predictive capacity of the method to that of a modern reactivity index, the average local ionization energy, I(r). The latter method is found to predict the regisolectivity in halogenations and nitrations qualitatively correctly if the positions for the I(r) minima (I(S,min)) are not too sterically hindered but fails for qualitative predictions of F-C reactions. The downscaled I(S,min) values also perform well for the quantitative prediction of regioisomer distributions of halogenations. The accuracy is slightly lower than that for the new method.

Ab initio embedded-cluster calculations have been performed for the CeO2(110) surface using temperature induced structures from molecular dynamics (MD) snapshots. As a first step towards understanding how temperature induced distortions of the surface structure influence the surface oxygen reactivity, the energy cost of removing an O atom from the surface was calculated for 41 snapshots from the MD simulation at 300 K. The quantum mechanical embedded-cluster calculations show that already at 300 K the dynamics causes significant fluctuations (root mean square of 0.37 eV) in the O vacancy formation energy (Evac) while the distribution of the two excess electrons associated with the vacancy is virtually unaffected by the surface dynamics and remains localized on the two Ce ions close to the vacancy. It is also found that the quantum mechanical Evac fluctuations can be reproduced by oxygen vacancy calculations using only the relaxed shell-model force field (FF) itself and the MD geometries. Using the FF as the interaction model, the effect of raising the temperature to 750 K and the effect of doping with Ca were investigated for the oxygen vacancy formation. [ABSTRACT FROM AUTHOR]

A combined molecular dynamics (MD)+quantum mechanics (QM) method for studying processes on ionic surfaces is presented. Through the combination of classical MD and ab initio embedded-cluster calculations, this method allows the modeling of surface processes involving both the structural and dynamic features of the substrate, even for large-scale systems. The embedding approach used to link the information from the MD simulation to the cluster calculation is presented, and rigorous tests have been carried out to ensure the feasibility of the method. The electrostatic potential and electron density resulting from our embedded-cluster model have been compared with periodic slab results, and confirm the satisfying quality of our embedding scheme as well as the importance of applying embedding in our combined MD+QM approach. We show that a highly accurate representation of the Madelung potential becomes a prerequisite when the embedded-cluster approach is applied to temperature-distorted surface snapshots from the MD simulation. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]