Recent advances in electronic structure theory and their influence on the accuracy of ab initio potential energy surfaces

Recent advances in electronic structure theory and the availability of high speed vector processors have substantially increased the accuracy of ab initio potential energy surfaces. The recently developed atomic natural orbital approach for basis set contraction has reduced both the basis set incompleteness and superposition errors in molecular calculations. Furthermore, full CI calculations can often be used to calibrate a CASSCF/MRCI approach that quantitatively accounts for the valence correlation energy. These computational advances also provide a vehicle for systematically improving the calculations and for estimating the residual error in the calculations. Calculations on selected diatomic and triatomic systems will be used to illustrate the accuracy that currently can be achieved for molecular systems. In particular, the F+H2 yields HF+H potential energy hypersurface is used to illustrate the impact of these computational advances on the calculation of potential energy surfaces

Density functional theory

Density functional theory (DFT) finds increasing use in applications related to biological systems. Advancements in methodology and implementations have reached a point where predicted properties of reasonable to high quality can be obtained. Thus, DFT studies can complement experimental investigations, or even venture with some confidence into experimentally unexplored territory. In the present contribution, we provide an overview of the properties that can be calculated with DFT, such as geometries, energies, reaction mechanisms, and spectroscopic properties. A wide range of spectroscopic parameters is nowadays accessible with DFT, including quantities related to infrared and optical spectra, X-ray absorption and Mössbauer, as well as all of the magnetic properties connected with electron paramagnetic resonance spectroscopy except relaxation times. We highlight each of these fields of application with selected examples from the recent literature and comment on the capabilities and limitations of current methods

A computational study of the heterogeneous synthesis of hydrazine on Co3Mo3N

Periodic and molecular density functional theory calculations have been applied to elucidate the associative mechanism for hydrazine and ammonia synthesis in the gas phase and hydrazine formation on Co3Mo3N. We find that there are two activation barriers for the associative gas phase mechanism with barriers of 730 and 658 kJ/mol, corresponding to a hydrogenation step from N2 to NNH2 and H2NNH2 to H3NNH3, respectively. The second step of the mechanism is barrierless and an important intermediate, NNH2, can also readily form on Co3Mo3N surfaces via the Eley–Rideal chemisorption of H2 on a pre-adsorbed N2 at nitrogen vacancies. Based on this intermediate a new heterogeneous mechanism for hydrazine synthesis is studied. The highest relative barrier for this heterogeneous catalysed process is 213 kJ/mol for Co3Mo3N containing nitrogen vacancies, clearly pointing towards a low-energy process for the synthesis of hydrazine via a heterogeneous catalysis route

Is Mn-Bound Substrate Water Protonated in the S2 State of Photosystem II?

In spite of great progress in resolving the geometric structure of the water-splitting Mn4OxCa cluster in photosystem II, the binding sites and modes of the two substrate water molecules are still insufficiently characterized. While time-resolved membrane-inlet mass spectrometry measurements indicate that both substrate water molecules are bound to the oxygen-evolving complex (OEC) in the S2 and S3 states (Hendry and Wydrzynski in Biochemistry 41:13328–13334, 2002), it is not known (1) if they are both Mn-bound, (2) if they are terminal or bridging ligands, and (3) in what protonation state they are bound in the different oxidation states Si (i = 0, 1, 2, 3, 4) of the OEC. By employing 17O hyperfine sublevel correlation (HYSCORE) spectroscopy we recently demonstrated that in the S2 state there is only one (type of) Mn-bound oxygen that is water exchangeable. We therefore tentatively identified this oxygen as one substrate ‘water’ molecule, and on the basis of the finding that it has a hyperfine interaction of about 10 MHz with the electron spin of the Mn4OxCa cluster, we suggest that it is bound as a Mn–O–Mn bridge within a bis-μ2 oxo-bridged unit (Su et al. in J Am Chem Soc 130:786–787, 2008). Employing pulse electron paramagnetic resonance, 1H/2H Mims electron-nuclear double resonance and 2H-HYSCORE spectroscopies together with 1H/2H-exchange here, we test this hypothesis by probing the protonation state of this exchangeable oxygen. We conclude that this oxygen is fully deprotonated. This result is discussed in the light of earlier reports in the literature

Coupling and uncoupling mechanisms in the methoxythreonine mutant of cytochrome P450cam: a quantum mechanical/molecular mechanical study

The Thr252 residue plays a vital role in the catalytic cycle of cytochrome P450cam during the formation of the active species (Compound I) from its precursor (Compound 0). We investigate the effect of replacing Thr252 by methoxythreonine (MeO-Thr) on this protonation reaction (coupling) and on the competing formation of the ferric resting state and H2O2 (uncoupling) by combined quantum mechanical/molecular mechanical (QM/MM) methods. For each reaction, two possible mechanisms are studied, and for each of these the residues Asp251 and Glu366 are considered as proton sources. The computed QM/MM barriers indicate that uncoupling is unfavorable in the case of the Thr252MeO-Thr mutant, whereas there are two energetically feasible proton transfer pathways for coupling. The corresponding rate-limiting barriers for the formation of Compound I are higher in the mutant than in the wild-type enzyme. These findings are consistent with the experimental observations that the Thr252MeO-Thr mutant forms the alcohol product exclusively (via Compound I), but at lower reaction rates compared with the wild-type enzyme

Orbital interactions and chemical reactivity of metal particles and metal surfaces

A review is presented with 101 refs. on chem. bonding to metal surfaces and small metal particles demonstrating the power of symmetry concepts to predict changes in chem. bonding. Ab-initio calcns. of chemisorption to small particles, as well as semiempirical extended Hueckel calcns. applied to the study of the reactivity of metal slabs are reviewed. On small metal particles, classical notions of electron promotion and hybridization are found to apply. The surroundings of a metal atom (ligands in complexes, other metal atoms at surfaces), affect bonding and reactivity through the prehybridization they induce. A factor specific for large particles and surfaces is the required localization of electrons on the atoms involved in the metal surface bond. At the surface, the bond energy is found to relate to the grou8p orbital local d. of states at the Fermi level. The use of this concept is extensively discussed and illustrated for chemisorption of CO and dissocn. of NO on metal surfaces. A discussion is given of the current decompn. schemes of bond energies and related concepts (exchange (Pauli)-repulsion, polarization, charge transfer). The role of non-orthogonality of fragment orbitals and of kinetic and potential energy for Pauli repulsion and (orbital) polarization is analyzed. Numerous examples are discussed to demonstrate the impact of those concepts on chem. bonding theor

Improving the study of proton transfers between amino acid sidechains in solution: choosing appropriate DFT functionals and avoiding hidden pitfalls

We have studied the influence of implicit solvent models, inclusion of explicit water molecules, inclusion of vibrational effects, and density functionals on the quality of the predicted pK a of small amino acid side chain models. We found that the inclusion of vibrational effects and explicit water molecules is crucial to improve the correlation between the computed and the experimental values. In these micro-solvated systems, the best agreement between DFT-computed electronic energies and benchmark values is afforded by BHHLYP and B97-2. However, approaching experimental results requires the addition of more than three explicit water molecules, which generates new problems related to the presence of multiple minima in the potential energy surface. It thus appears that a satisfactory ab initio prediction of amino acid side chain pK a will require methods that sample the configurational space in the presence of large solvation shells, while at the same time computing vibrational contributions to the enthalpy and entropy of the system under study in all points of that surface. Pending development of efficient algorithms for those computations, we strongly suggest that whenever counterintuitive protonation states are found in a computational study (e.g., the presence of a neutral aspartate/neutral histidine dyad instead of a deprotonated aspartate/protonated histidine pair), the reaction profile should be computed under each of the different protonation micro-states by constraining the relevant N–H or O–H bonds, in order to avoid artifacts inherent to the complex nature of the factors contributing to the pK a