Density Functional Tight Binding (DFTB) models are two to three orders of magnitude faster than ab initio and Density Functional Theory (DFT) methods and therefore are particularly attractive in applications to large molecules and condensed phase systems. To establish the applicability of DFTB models to general chemical reactions, we conduct benchmark calculations for barrier heights and reaction energetics of organic molecules using existing databases and several new ones compiled in this study. Structures for the transition states and stable species have been fully optimized at the DFTB level, making it possible to characterize the reliability of DFTB models in a more thorough fashion compared to conducting single point energy calculations as done in previous benchmark studies. The encouraging results for the diverse sets of reactions studied here suggest that DFTB models, especially the most recent third-order version (DFTB3/3OB augmented with dispersion correction), in most cases provide satisfactory description of organic chemical reactions with accuracy almost comparable to popular DFT methods with large basis sets, although larger errors are also seen for certain cases. Therefore, DFTB models can be effective for mechanistic analysis (e.g., transition state search) of large (bio)molecules, especially when coupled with single point energy calculations at higher levels of theory.

Please note that technical editing may introduce minor changes to the text and/or graphics contained in the manuscript submitted by the author(s) which may alter content, and that the standard Terms & Conditions and the ethical guidelines that apply to the journal are still applicable. In no event shall the RSC be held responsible for any errors or omissions in these Accepted Manuscript manuscripts or any consequences arising from the use of any information contained in them.

The aim of this work was to determine and understand the origin of the electronic properties of Mn(IV) complexes, especially the zero-field splitting (ZFS), through a combined experimental and theoretical investigation on five well-characterized mononuclear octahedral Mn(IV) compounds, with various coordination spheres (N6, N3O3, N2O4 in both trans (trans-N2O4) and cis configurations (cis-N2O4) and O4S2). High-frequency and -field EPR (HFEPR) spectroscopy has been applied to determine the ZFS parameters of two of these compounds, MnL(trans-N2O4) and MnL(O4S2). While at X-band EPR, the axial-component of the ZFS tensor, D, was estimated to be +0.47 cm(-1) for MnL(O4S2), and a D-value of +2.289(5) cm(-1) was determined by HFEPR, which is the largest D-magnitude ever measured for a Mn(IV) complex. A moderate D value of -0.997(6) cm(-1) has been found for MnL(trans-N2O4). Quantum chemical calculations based on two theoretical frameworks (the Density Functional Theory based on a coupled perturbed approach (CP-DFT) and the hybrid Ligand-Field DFT (LF-DFT)) have been performed to define appropriate methodologies to calculate the ZFS tensor for Mn(IV) centers, to predict the orientation of the magnetic axes with respect to the molecular ones, and to define and quantify the physical origin of the different contributions to the ZFS. Except in the case of MnL(trans-N2O4), the experimental and calculated D values are in good agreement, and the sign of D is well predicted, LF-DFT being more satisfactory than CP-DFT. The calculations performed on MnL(cis-N2O4) are consistent with the orientation of the principal anisotropic axis determined by single-crystal EPR, validating the calculated ZFS tensor orientation. The different contributions to D were analyzed demonstrating that the d-d transitions mainly govern D in Mn(IV) ion. However, a deep analysis evidences that many factors enter into the game, explaining why no obvious magnetostructural correlations can be drawn in this series of Mn(IV) complexes.

In a combined experimental and theoretical study we characterize dissociative electron attachment (DEA) to, and electronically excited states of, Fe(CO). Both are relevant for electron-induced degradation of Fe(CO). The strongest DEA channel is cleavage of one metal-ligand bond that leads to production of Fe(CO). High-resolution spectra of Fe(CO) reveal fine structures at the onset of vibrational excitation channels. Effective range R-matrix theory successfully reproduces these structures as well as the dramatic rise of the cross section at very low energies and reveals that virtual state scattering dominates low-energy DEA in Fe(CO) and that intramolecular vibrational redistribution (IVR) plays an essential role. The virtual state hypothesis receives further experimental support from the rapid rise of the elastic cross section at very low energies and intense threshold peaks in vibrational excitation cross sections. The IVR hypothesis is confirmed by our measurements of kinetic energy distributions of the fragment ions, which are narrow (∼0.06 eV) and peak at low energies (∼0.025 eV), indicating substantial vibrational excitation in the Fe(CO) fragment. Rapid IVR is also revealed by the yield of thermal electrons, observed in two-dimensional (2D) electron energy loss spectroscopy. We further measured mass-resolved DEA spectra at higher energies, up to 12 eV, and compared the bands observed there to resonances revealed by the spectra of vibrational excitation cross sections. Dipole-allowed and dipole/spin forbidden electronic transitions in Fe(CO)-relevant for neutral dissociation by electron impact-are probed using electron energy loss spectroscopy and time-dependent density functional theory calculations. Very good agreement between theory and experiment is obtained, permitting assignment of the observed bands.

Electron induced chemistry of metal-containing precursor molecules is central in focused electron beam induced deposition (FEBID). While some elementary processes leading to precursor decomposition were quantitatively characterized, data for neutral dissociation is missing. We provide this data for the model precursor Pt­(PF3)4 by using the available cross sections for electronic excitation and characterizing fragmentation of the excited states theoretically by TDDFT. The potential energy curves for a number of states visible in the experimental electron energy loss spectra are dissociative, either directly or via conical intersections, indicating that the quantum yield for dissociation is close to 100%. Taking into account typical electron energy distribution at the FEBID spot reveals that the importance of neutral dissociation exceeds that of dissociative electron attachment, which has been so far considered to be the dominant decomposition process. We thus established neutral dissociation as an important, albeit often neglected, channel for FEBID using Pt­(PF3)4. The calculations revealed a number of other phenomena that can play a role in electron induced chemistry of this compound, e.g., a considerable increase of bond dissociation energy with sequential removal of multiple ligands.

A detailed discussion of the potential energy surface of bis(cyclopentadienyl)cobalt(II), cobaltocene, is given. Vibronic coupling coefficients are calculated using density functional theory (DFT). Results are in good agreement with experimental findings. On the basis of our calculation, there is no second-order Jahn-Teller (JT) effect as predicted by group theory. The JT distortion can be expressed as a linear combination of all totally symmetric normal modes of the low-symmetry, minimum-energy conformation. The out-ofplane ring deformation is the most important mode. The JT distortion is analyzed by seeking the path of minimal energy of the adiabatic potential energy surface.

A general model for the analysis of the Adiabatic Potential Energy Surfaces (APES) of the molecules that are subject to the multimode Jahn-Teller effect is presented. The method utilizes the information obtained by DFT calculations on a distorted stationary point on the APES. The essence of the model is to express the distortion along a model minimal energy path called Intrinsic Distortion Path (IDP), projecting the geometry of the system on the normal modes of the either high-symmetry (HS) or low symmetry (LS) nuclear configuration. This allows us to determine the significance of all of the involved normal modes along a relevant particular path of distortion, and the direct calculation of the vibronic coupling constants. The IDP analysis is illustrated by the discussion of the multimode H ⊗ (g + 2h) JT effect in fullerene cations (C(60)(+)) giving a deep insight into the origin and the mechanism of vibronic coupling in fullerene based molecules.

No abstract