scholarly journals Magnetic Exchange Interactions in Binuclear and Tetranuclear Iron (III) Complexes Described by Spin-Flip DFT and Heisenberg Effective Hamiltonians

Author(s):  
Saikiran Kotaru ◽  
Maristella Alessio ◽  
Anna I. Krylov

Low-energy spectra of single-molecule magnets (SMMs) are often described by the Heisenberg Hamiltonian. Within this formalism, exchange interactions between magnetic centers determine the ground-state multiplicity and energy separation between the ground and excited states. In this contribution, we extract exchange coupling constants (J) for a set of iron (III) binuclear and tetranuclear complexes from all-electron calculations using non-collinear spin-flip time-dependent density functional theory (NC-SF-TDDFT). For the series of binuclear complexes with J-values ranging from -6 to -132 cm−1 , our benchmark calculations using the short-range hybrid LRC-ωPBEh functional and 6-31G(d,p) basis set agree well (mean absolute error of 4.7 cm−1) with the experimentally derived values. For the tetranuclear SMMs, the computed J constants are within 6 cm−1 from the values extracted from the experiment. We explore the range of applicability of the Heisenberg model by analyzing the radical character in the binuclear iron (III) complexes using natural orbitals (NO) and their occupations. On the basis of the number of effectively unpaired electrons and the NO occupancies, we attribute larger errors observed in strongly anti-ferromagnetic species to an increased ionic character. The results illustrate the efficiency of the spin-flip protocol for computing the exchange couplings and the utility of the NO analysis in assessing the validity of effective spin Hamiltonians.

2021 ◽  
Author(s):  
Saikiran Kotaru ◽  
Maristella Alessio ◽  
Anna I. Krylov

Low-energy spectra of single-molecule magnets (SMMs) are often described by the Heisenberg Hamiltonian. Within this formalism, exchange interactions between magnetic centers determine the ground-state multiplicity and energy separation between the ground and excited states. In this contribution, we extract exchange coupling constants (J) for a set of iron (III) binuclear and tetranuclear complexes from all-electron calculations using non-collinear spin-flip time-dependent density functional theory (NC-SF-TDDFT). For the series of binuclear complexes with J-values ranging from -6 to -132 cm−1 , our benchmark calculations using the short-range hybrid LRC-ωPBEh functional and 6-31G(d,p) basis set agree well (mean absolute error of 4.7 cm−1) with the experimentally derived values. For the tetranuclear SMMs, the computed J constants are within 6 cm−1 from the values extracted from the experiment. We explore the range of applicability of the Heisenberg model by analyzing the radical character in the binuclear iron (III) complexes using natural orbitals (NO) and their occupations. On the basis of the number of effectively unpaired electrons and the NO occupancies, we attribute larger errors observed in strongly anti-ferromagnetic species to an increased ionic character. The results illustrate the efficiency of the spin-flip protocol for computing the exchange couplings and the utility of the NO analysis in assessing the validity of effective spin Hamiltonians.


2013 ◽  
Vol 91 (9) ◽  
pp. 866-871 ◽  
Author(s):  
Silvia Gómez-Coca ◽  
Eliseo Ruiz

The exchange coupling constants of a Mn14 complex constituted by two weakly coupled Mn7 moieties were calculated using two different density functional theory (DFT) approaches: the Perdew–Burke–Ernzerhof (PBE) functional with a numerical basis set and the hybrid Becke, three-parameter Lee–Yang–Parr (B3LYP) functional employed with a Gaussian basis set. The sign and relative strength of the exchange coupling constants calculated with both methods were consistent; as expected, the values calculated with the PBE functional were slightly overestimated, as corroborated by comparison with the experimental magnetic susceptibility curve. Both methods gave a ground spin configuration of S = 3/2 for the Mn7 moiety, which was weakly antiferromagnetically coupled with the other Mn7 fragment, leading to an S = 0 ground spin configuration for the entire Mn14 complex.


2021 ◽  
Author(s):  
Yannick J. Franzke ◽  
Jason M. Yu

We present a highly efficient implementation of the electron-nucleus hyperfine coupling matrix within one-electron exact two-component (X2C) theory. The complete derivative of the X2C Hamiltonian is formed, i.e. the derivatives of the unitary decoupling transformation are considered. This requires solution of the response and Sylvester equations, consequently increasing the computational costs. Therefore, we apply the diagonal local approximation to the unitary decoupling transformation (DLU). The finite nucleus model is employed for both the scalar potential and the vector potential. Two-electron picture-change effects are modeled with the (modified) screened-nuclear spin--orbit approach. Our implementation is fully integral direct and OpenMP-parallelized. An extensive benchmark study regarding the Hamiltonian, the basis set, and the density functional approximation is carried out for a set of 12--17 transition-metal compounds. The error introduced by DLU is negligible and the DLU-X2C Hamiltonian accurately reproduces its four-component ``fully'' relativistic parent results. Functionals with a large amount of Hartree--Fock exchange such as CAM-QTP-02 and omega-B97X-D are generally favorable. The pure density functional r2SCAN performs remarkably and even outperforms the common hybrid functionals TPSSh and CAM-B3LYP. Fully uncontracted basis sets or contracted quadruple-zeta bases are required for accurate results. The capability of our implementation is demonstrated for [Pt(C6Cl5)4]- with more than 4700 primitive basis functions and four rare-earth single molecule magnets: [La(OAr*)3]-, [Lu(NR2)3]-, [Lu(OAr*)3]-, and [TbPc2]-. Here, the spin--orbit DLU-X2C Hamiltonian results in an excellent agreement with the experimental findings of all Pt, La, Lu, and Tb molecules.


2021 ◽  
Author(s):  
Yannick J. Franzke ◽  
Jason M. Yu

We present a highly efficient implementation of the electron-nucleus hyperfine coupling matrix within one-electron exact two-component (X2C) theory. The complete derivative of the X2C Hamiltonian is formed, i.e. the derivatives of the unitary decoupling transformation are considered. This requires solution of the response and Sylvester equations, consequently increasing the computational costs. Therefore, we apply the diagonal local approximation to the unitary decoupling transformation (DLU). The finite nucleus model is employed for both the scalar potential and the vector potential. Two-electron picture-change effects are modeled with the (modified) screened-nuclear spin--orbit approach. Our implementation is fully integral direct and OpenMP-parallelized. An extensive benchmark study regarding the Hamiltonian, the basis set, and the density functional approximation is carried out for a set of 12--17 transition-metal compounds. The error introduced by DLU is negligible and the DLU-X2C Hamiltonian accurately reproduces its four-component ``fully'' relativistic parent results. Functionals with a large amount of Hartree--Fock exchange such as CAM-QTP-02 and omega-B97X-D are generally favorable. The pure density functional r2SCAN performs remarkably and even outperforms the common hybrid functionals TPSSh and CAM-B3LYP. Fully uncontracted basis sets or contracted quadruple-zeta bases are required for accurate results. The capability of our implementation is demonstrated for [Pt(C6Cl5)4]- with more than 4700 primitive basis functions and four rare-earth single molecule magnets: [La(OAr*)3]-, [Lu(NR2)3]-, [Lu(OAr*)3]-, and [TbPc2]-. Here, the spin--orbit DLU-X2C Hamiltonian results in an excellent agreement with the experimental findings of all Pt, La, Lu, and Tb molecules.


2021 ◽  
Vol 7 (2) ◽  
pp. 22
Author(s):  
Itziar Oyarzabal ◽  
Estitxu Echenique-Errandonea ◽  
Eider San Sebastián ◽  
Antonio Rodríguez-Diéguez ◽  
José Manuel Seco ◽  
...  

New dinuclear MII-LnIII complexes of general formulas [Cu(µ-L)(µ-OAc)Ln(NO3)2]·CH3CN·H2O (LnIII = Gd (1), Tb (2), Dy (3) and Er (4)), [Ni(CH3CN)(µ-L)(µ-OAc)Ln(NO3)2]·CH3CN (LnIII = Nd (5), Gd (6), Tb (7), Dy (8), Er (9) and Y (10)) and [Co(CH3CN)(µ-L)(µ-OAc)Ln(NO3)2]·CH3CN (LnIII = Gd (11), Tb (12), Dy (13), Er (14) and Y (15)) were prepared from the compartmental ligand N,N′-dimethyl-N,N′-bis(2-hydroxy-3-formyl-5-bromo-benzyl)ethylenediamine (H2L). In all these complexes, the transition metal ions occupy the internal N2O2 coordination site of the ligand, whereas the LnIII ions lie in the O4 external site. Both metallic ions are connected by an acetate bridge, giving rise to triple mixed diphenoxido/acetate bridged MIILnIII compounds. Direct current (dc) magnetic measurements allow the study of the magnetic exchange interactions between the 3d and 4f metal ions, which is supported by density functional theory (DFT) theoretical calculations for the GdIII-based counterparts. Due to the weak ferromagnetic exchange coupling constants obtained both experimentally and theoretically, the magneto-thermal properties of the less anisotropic systems (compounds 1 and 6) are also studied. Alternating current (ac)magnetic measurements reveal the occurrence of slight frequency dependency of the out-of-phase signal for complexes 8, 9 and 13, while complex 15 displays well-defined maximums below ~6 K.


2011 ◽  
Vol 10 (03) ◽  
pp. 325-348 ◽  
Author(s):  
ANNA IGNACZAK

Detailed comparative analysis of properties of the tert-butyl radical and cation is performed using 14 density functional (DFT) methods combined with double-zeta and triple-zeta quality Gaussian basis sets with polarization and diffuse functions. Stability of different conformers is discussed. Structural parameters, dipole moment, adiabatic ionization potential (IP), inversion barrier and isotropic hyperfine coupling constants are examined and compared to values obtained at the standard MP2 level and to experimental data available. All methods indicate that that the CC bond in the radical is longer than in the cation by about 0.033 Å. The IP values are found to be very sensitive to the method used and range from 612 to 709 kJ/mol, but majority oscillate around 646÷656 kJ/mol. Calculated inversion barrier for the radical is higher than the experimental estimate of 2.68 kJ/mol; with the 6-311++G** basis set and most DFT methods it is predicted in the range 3.86÷4.82 kJ/mol. All DFT methods predict for the out-of-plane CC3 bending mode of the radical the frequency around 260 cm-1, while in the cation the corresponding frequency is higher by about 180 cm-1.


Molecules ◽  
2020 ◽  
Vol 25 (17) ◽  
pp. 3931 ◽  
Author(s):  
Kacper Rzepiela ◽  
Aneta Buczek ◽  
Teobald Kupka ◽  
Małgorzata A. Broda

We report on the density functional theory (DFT) modelling of structural, energetic and NMR parameters of uracil and its derivatives (5-halogenouracil (5XU), X = F, Cl, Br and I) in vacuum and in water using the polarizable continuum model (PCM) and the solvent model density (SMD) approach. On the basis of the obtained results, we conclude that the intramolecular electrostatic interactions are the main factors governing the stability of the six tautomeric forms of uracil and 5XU. Two indices of aromaticity, the harmonic oscillator model of aromaticity (HOMA), satisfying the geometric criterion, and the nuclear independent chemical shift (NICS), were applied to evaluate the aromaticity of uracil and its derivatives in the gas phase and water. The values of these parameters showed that the most stable tautomer is the least aromatic. A good performance of newly designed xOPBE density functional in combination with both large aug-cc-pVQZ and small STO(1M)−3G basis sets for predicting chemical shifts of uracil and 5-fluorouracil in vacuum and water was observed. As a practical alternative for calculating the chemical shifts of challenging heterocyclic compounds, we also propose B3LYP calculations with small STO(1M)−3G basis set. The indirect spin–spin coupling constants predicted by B3LYP/aug-cc-pVQZ(mixed) method reproduce the experimental data for uracil and 5-fluorouracil well.


Sign in / Sign up

Export Citation Format

Share Document