Extraction of giant Dzyaloshinskii-Moriya interaction from ab initio calculations: First-order spin-orbit coupling model and methodological study

Physical origin of the anisotropic exchange tensor close to the first-order spin-orbit coupling regime and impact of the electric field on its magnitude

This article follows earlier studies on the physical origin of magnetic anisotropy and the means of controlling it in polynuclear transition metal complexes. The difficulties encountered when focusing a magnetic field on a molecular object have led to consider the electric field as a more appropriate control tool. It is therefore fundamental to understand what governs the sensitivity of magnetic properties to the application of an electric field. We have already studied the impact of the electric field on the isotropic exchange coupling and on the Dzyaloshinskii-Moriya interaction (DMI). Here, we focus on the symmetric exchange anisotropy tensor. In order to obtain significant values of anisotropic interactions, we have carried out this study on a model complex that exhibits first-order spin-orbit coupling. We will show that (i) large values of the axial parameter of symmetric exchange can be reached when close to the first-order spin-orbit coupling regime, (ii) both correlated energies and wave functions must be used to achieve accurate values of the symmetric tensor components when the DMI is non-zero, and (iii) finally, an interferential effect between the DMI and the axial parameter of symmetric exchange occurs for a certain orientation of the electric field, i.e., the latter decreases in magnitude as the former increases. While DMI is often invoked as being involved in magneto-electric coupling, isotropic exchange and the symmetrical anisotropic tensor also contribute. Finally, we provide a recipe for generating significant anisotropic interactions and a significant change in magnetic properties under an electric field.

The resolution of the weak-exchange limit made rigorous, simple and general in binuclear complexes

The correct interpretation of magnetic properties in the weak-exchange regime has remained a challenging task for several decades. In this regime, the effective exchange interaction between local spins is quite weak, of the same order of magnitude or smaller than the various anisotropic terms, which generates a complex set of levels characterized by spin mixing. Although the model multispin Hamiltonian in the absence of local orbital momentum, , is considered good enough to map the experimental energies at zero field and in the strong-exchange limit, theoretical works pointed out limitations of this simple model. This work revives the use of ĤMS from a new theoretical perspective, detailing point-by-point a strategy to correctly map the computational energies and wave functions onto ĤMS, thus validating it regardless of the exchange limit. We will distinguish two cases, based on experimentally characterized dicobalt(II) complexes from the literature. If centrosymmetry imposes alignment of the various rank-2 tensors constitutive of ĤMS in the first case, the absence of any symmetry element prevents such alignment in the second case. In such a context, the strategy provided herein becomes a powerful tool to rationalize the experimental magnetic data, since it is capable of fully and rigorously extracting the multispin model without any assumption on the orientation of its constitutive tensors. Furthermore, the strategy allows to question the use of the spin Hamiltonian approach by explicitly controlling the projection norms on the model space, which is showcased in the second complex where local orbital momentum could have occurred (distorted octahedra). Finally, previous theoretical data related to a known dinickel(II) complex is reinterpreted, clarifying initial wanderings regarding the weak exchange limit.

Impact of the electric field on isotropic and anisotropic spin Hamiltonian parameters

One may obviously think that the best way to control magnetic properties relies on using a magnetic field. However, it is not convenient to focus a magnetic field on a small object, whereas it is much easier to do so with an electric field. Magnetoelectric coupling allows one to control the magnetization with the electric field and the polarization with the magnetic field and could therefore provide a solution to this problem. This paper aims at quantifying the impact of the electric field on both the isotropic magnetic exchange and the Dzyaloshinskii–Moriya interaction in the case of a binuclear system of S = 1/2 spins. This study follows previous studies that showed that very high Dzyaloshinskii–Moriya interaction, i.e., the antisymmetric exchange, can be generated when close to first order spin orbit coupling. We will, therefore, explore this regime in a model Cu(II) complex that exhibits a quasi-degeneracy of the [Formula: see text] and d xy orbitals. This situation is indeed the one that allows us to obtain the largest spin orbit couplings in transition metal complexes. We will show that both the magnetic exchange and the Dzyaloshinskii–Moriya interaction are very sensitive to the electric field and that it would therefore be possible to modulate and control magnetic properties by the electric field. Finally, rationalizations of the obtained results will be proposed.

Antisymmetric Exchange in a Real Copper Triangular Complex

The antisymmetric exchange, also known as the Dzyaloshinskii-Moriya interaction (DMI), is an effective interaction that may be at play in isolated complexes (with transition metals or lanthanides, for instance), nanoparticles, and highly correlated materials with adequate symmetry properties. While many theoretical works have been devoted to the analysis of single-ion zero-field splitting and to a lesser extent to symmetric exchange, only a few ab initio studies deal with the DMI. Actually, it originates from a subtle interplay between weak electronic interactions and spin-orbit couplings. This article aims to highlight the origin of this interaction from theoretical grounds in a real tri-copper(II) complex, capitalizing on previous methodological studies on bi-copper(II) model complexes. By tackling this three-magnetic-center system, we will first show that the multispin model Hamiltonian is appropriate for trinuclear (and likely for higher nuclearity) complexes, then that the correct application of the permutation relationship is necessary to explain the outcomes of the ab initio calculations, and finally, that the model parameters extracted from a binuclear model transfer well to the trinuclear complex. For a more theory-oriented purpose, we will show that the use of a simplified structural model allows one to perform more demanding electronic structure calculations. On this simpler system, we will first check that the previous transferability is still valid, prior to performing more advanced calculations on the derived two-magnetic-center model system. To this end, we will explain in detail the physics of the DMI in the copper triangle of interest, before advocating further theory/experiment efforts.

On the role of dynamic electron correlation in non-orthogonal configuration interaction with fragments

Two different approaches have been implemented to include the effect of dynamic electron correlation in the Non-Orthogonal Configuration Interaction for Fragments (NOCI-F) method. The first is based on shifting the diagonal matrix elements of the NOCI matrix, while the second incorporates the dynamic correlation explicitly in the fragment wave functions used to construct the many-electron basis functions of the NOCI. The two approaches are illustrated for the calculation of the electronic coupling relevant in singlet fission and the coupling of spin moments in organic radicals. Comparison of the calculated diabatic couplings, the NOCI energies and wave functions shows that dynamic electron correlation is not only efficiently but also effectively incorporated by the shifting approach and can largely affect the coupling between electronic states. Also, it brings the NOCI coupling of the spin moments in close agreement with benchmark calculations.

Zero-Field Splitting Calculations by Multiconfiguration Pair-Density Functional Theory

Zero-field splitting (ZFS) is a fundamental molecular property that is especially relevant for single-molecule magnets (SMMs), electron paramagnetic resonance spectra, and quantum computing. Developing a method that can accurately predict ZFS parameters can be very powerful for designing new SMMs. One of the challenges is to include external correlation in an inherently multiconfigurational open-shell species for the accurate prediction of magnetic properties. Previously available methods depend on expensive multireference perturbation theory calculations to include external correlation. In this paper, we present spin-orbit-inclusive multiconfiguration and multistate pair-density functional theory (MC-PDFT) calculations of ZFSs; these calculations have a cost comparable to complete-active-space self-consistent field (CASSCF) theory, but they include correlation external to the active space. We found that combining a multistate formulation of MC-PDFT, namely, compressed-state multistate pair-density functional theory, with orbitals optimized by weighted-state-averaged CASSCF, yields reasonably accurate ZFS results.