Matija Zlatar

Theoretical study of the magnetic anisotropy and magnetic tunnelling in mononuclear Ni(II) complexes with potential molecular magnet behavior

Magnetic molecules that present a slow decay of their magnetization (molecular magnets) are very interesting both from a fundamental and applied points of view. While many approaches focus strongly on finding systems with strong magnetic anisotropy giving rise to large spin-reversal barriers, less is known on the behavior of magnetic tunnelling, which is also a fundamental component of molecular magnet behavior. In this work, we propose a model to describe both the spin-reversal barrier and magnetic tunnelling in Ni(II) trigonal bipyramidal complexes, which could be easily extended to other transition-metal systems. Based on this model, we show the criteria that lead to the optimal complexes to find molecular magnet behavior. We test our proposal with multi-reference configuration-interaction (MRCI) and ligand-field-density-functional-theory (LF-DFT) first-principles calculations applied over several families of mononuclear Ni(II) complexes. As a salient result, we find that the complex [NiCl3(Hdabco)2]+ (dabco is 1,4-diazabicyclo[2.2.2]-octane) displays both a very large magnetic anisotropy energy, 524 cm−1, and a small tunnelling splitting, 0.2 cm−1, when compared to other systems containing the same metal. We expect molecular magnet behaviour to be observed when small magnetic fields are employed to disrupt tunnelling. These values are reached due to the choice of ligands that favor a complete destruction of the Jahn–Teller distortions through spin–orbit coupling and an unquenched orbital momentum.

Density functional theory study of the Jahn-Teller effect in cobaltocene*

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-of-plane 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.

Origin of the Zero-Field Splitting in Mononuclear Octahedral MnIV Complexes: A Combined Experimental and Theoretical Investigation

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.

Treatment of the multimode Jahn-Teller problem in small aromatic radicals.

The family of the Jahn-Teller (JT) active hydrocarbon rings, C(n)H(n) (n = 5-7), was analyzed by the means of multideterminantal density functional theory (DFT) approach. The multimode problem was addressed using the intrinsic distortion path (IDP) method, in which the JT distortion is expressed as a linear combination of all totally symmetric normal modes in the low symmetry minimum energy conformation. Partitioning of the stabilization energy into the various physically meaningful terms arising from Kohn-Sham DFT has been performed to get further chemical insight into the coupling of the nuclear movements and the electron distribution.

Assessment of TD-DFT and LF-DFT for study of d - d transitions in first row transition metal hexaaqua complexes.

Herein, we present the systematic, comparative computational study of the d - d transitions in a series of first row transition metal hexaaqua complexes, [M(H2O)6](n+) (M(2+/3+) = V (2+/3+), Cr(2+/3+), Mn(2+/3+), Fe(2+/3+), Co(2+/3+), Ni(2+)) by the means of Time-dependent Density Functional Theory (TD-DFT) and Ligand Field Density Functional Theory (LF-DFT). Influence of various exchange-correlation (XC) approximations have been studied, and results have been compared to the experimental transition energies, as well as, to the previous high-level ab initio calculations. TD-DFT gives satisfactory results in the cases of d(2), d(4), and low-spin d(6) complexes, but fails in the cases when transitions depend only on the ligand field splitting, and for states with strong character of double excitation. LF-DFT, as a non-empirical approach to the ligand field theory, takes into account in a balanced way both dynamic and non-dynamic correlation effects and hence accurately describes the multiplets of transition metal complexes, even in difficult cases such as sextet-quartet splitting in d(5) complexes. Use of the XC functionals designed for the accurate description of the spin-state splitting, e.g., OPBE, OPBE0, or SSB-D, is found to be crucial for proper prediction of the spin-forbidden excitations by LF-DFT. It is shown that LF-DFT is a valuable alternative to both TD-DFT and ab initio methods.

Magnetic anisotropy in ‘scorpionate’ first-row transition-metal complexes: a theoretical investigation

In this work we have analyzed in detail the magnetic anisotropy in a series of hydrotris(pyrazolyl)borate (Tp(-)) metal complexes, namely [VTpCl](+), [CrTpCl](+), [MnTpCl](+), [FeTpCl], [CoTpCl], and [NiTpCl], and their substituted methyl and tert-butyl analogues with the goal of observing the effect of the ligand field on the magnetic properties. In the [VTpCl](+), [CrTpCl](+), [CoTpCl], and [NiTpCl] complexes, the magnetic anisotropy arises as a consequence of out-of-state spin-orbit coupling, and covalent changes induced by the substitution of hydrogen atoms on the pyrazolyl rings does not lead to drastic changes in the magnetic anisotropy. On the other hand, much larger magnetic anisotropies were predicted in complexes displaying a degenerate ground state, namely [MnTpCl](+) and [FeTpCl], due to in-state spin-orbit coupling. The anisotropy in these systems was shown to be very sensitive to perturbations, for example, chemical substitution and distortions due to the Jahn-Teller effect. We found that by substituting the hydrogen atoms in [MnTpCl](+) and [FeTpCl] by methyl and tert-butyl groups, certain covalent contributions to the magnetic anisotropy energy (MAE) could be controlled, thereby achieving higher values. Moreover, we showed that the selection of ion has important consequences for the symmetry of the ground spin-orbit term, opening the possibility of achieving zero magnetic tunneling even in non-Kramers ions. We have also shown that substitution may also contribute to a quenching of the Jahn-Teller effect, which could significantly reduce the magnetic anisotropy of the complexes studied.

A Practical Computational Approach to Study Molecular Instability Using the Pseudo-Jahn-Teller Effect.

Vibronic coupling theory shows that the cause for spontaneous instability in systems presenting a nondegenerate ground state is the so-called pseudo-Jahn-Teller effect, and thus its study can be extremely helpful to understand the structure of many molecules. While this theory, based on the mixing of the ground and excited states with a distortion, has been long studied, there are two obscure points that we try to clarify in the present work. First, the operators involved in both the vibronic and nonvibronic parts of the force constant take only into account electron-nuclear and nuclear-nuclear interactions, apparently leaving electron-electron repulsions and the electrons kinetic energy out of the chemical picture. Second, a fully quantitative computational appraisal of this effect has been up to now problematic. Here, we present a reformulation of the pseudo-Jahn-Teller theory that explicitly shows the contributions of all operators in the molecular Hamiltonian and allows connecting the results obtained with this model to other chemical theories relating electron distribution and geometry. Moreover, we develop a practical approach based on Hartree-Fock and density functional theory that allows quantification of the pseudo-Jahn-Teller effect. We demonstrate the usefulness of our method studying the pyramidal distortion in ammonia and its absence in borane, revealing the strong importance of the kinetic energy of the electrons in the lowest a2″ orbital to trigger this instability. The present tool opens a window for exploring in detail the actual microscopic origin of structural instabilities in molecules and solids.

Role of spin state and ligand charge in coordination patterns in complexes of 2,6-diacetylpyridinebis(semioxamazide) with 3d-block metal ions: A density functional theory study

We report here a systematic computational study on the effect of the spin state and ligand charge on coordination preferences for a number of 3d-block metal complexes with the 2,6-diacetylpyridinebis(semioxamazide) ligand and its mono- and dianionic analogues. Our calculations show excellent agreement for the geometries compared with the available X-ray structures and clarify some intriguing experimental observations. The absence of a nickel complex in seven-coordination is confirmed here, which is easily explained by inspection of the molecular orbitals that involve the central metal ion. Moreover, we find here that changes in the spin state lead to completely different coordination modes, in contrast to the usual situation that different spin states mainly result in changes in the metal-ligand bond lengths. Both effects result from different occupations of a combination of π- and σ-antibonding and nonbonding orbitals.

Benchmarking density functional tight binding models for barrier heights and reaction energetics of organic molecules

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.

Resolving the origin of the multimode Jahn–Teller effect in metallophthalocyanines

A detailed Density Functional Theory (DFT) analysis was performed in order to study the multimode Jahn-Teller (JT) problem in the electronic ground state of manganese phthalocyanine (MnPc). A comparison with the magnesium phthalocyanine ion (MgPc-) and the phthalocyanine trianion (Pc3-), also prone to the JT effect, is presented. Our results clarify the origin and provide the microscopic insight into the symmetry breaking process. The JT distortion is highly influenced by the coordination of phthalocyanine to the MnII ion, and occurs over the whole system, while the MgPc- complex ion possesses mainly ligand-based instability.