Understanding the Electronic Structure of Magnetic Trinuclear Complexes Based on the Tris-Dioxolene Triphenylene Non-Innocent Bridging Ligand, a Theoretical Study.

A complete theoretical analysis using first the simple Hückel model followed by more sophisticated multi-reference calculations on a trinuclear Ni(II) complex (Tp#Ni3 HHTP), bearing the non-innocent bridging ligand HHTP3- , is carried out. The three semiquinone moieties of HHTP3- couple antiferromagnetically and lead to a single unpaired electron localized on one of the moieties. The calculated exchange coupling integrals together with the zero-field parameters allow, when varied within a certain range, reproducing the experimental data. These results are generalized for two similar other trinuclear complexes containing Ni(II) and Cu(II). The electronic structure of HHTP3- turns out to be independent of both the chemical nature and the geometry of the metal ions. We also establish a direct correlation between the geometrical and the electronic structures of the non-innocent ligand that is consistent with the results of calculations. It allows experimentalists to get insight into the magnetic behavior of this type of complexes by an analysis of their X-ray structure.

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.

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 dx2-y2 and dxy 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.

Difficulty of the evaluation of the barrier height of an open-shell transition state between closed shell minima: The case of small C4n rings.

C4n cyclacenes exhibit strong bond-alternation in their equilibrium geometry. In the two equivalent geometries, the system keeps an essentially closed-shell character. The two energy minima are separated by a transition state suppressing the bond-alternation, where the wave function is strongly diradical. This paper discusses the physical factors involved in this energy difference and possible evaluations of the barrier height. The barrier given as the energy difference between the restricted density functional theory (DFT)/B3LYP for the equilibrium and the broken symmetry DFT/B3LYP of the transition state is either negative or small, in contradiction with the most reliable Wave Function Theory calculations. The minimal (two electrons in two molecular orbitals) Complete Active Space self-consistent field (CASSCF) overestimates the barrier, and the subsequent second-order perturbation cancels it. Due to the collective character of the spin-polarization effect, it is necessary to perform a full π CASSCF + second-order perturbation to reach a reasonable value of the barrier, but this type of treatment cannot be applied to large molecules. DFT procedures treating on an equal foot the closed-shell and open-shell geometries have been explored, such as Mixed-Reference Spin-Flip Time-dependent-DFT and a new spin-decontamination proposal, namely, DFT-dressed configuration interaction, but the results still depend on the density functional. M06-2X without or with spin-decontamination gives the best agreement with the accurate wave function results.

Trinuclear Cyanido-Bridged [Cr2 Fe] Complexes: To Be or not to Be a Single-Molecule Magnet, a Matter of Straightness.

Trinuclear systems of formula [{Cr(LN3O2Ph )(CN)2 }2 M(H2 LN3O2R )] (M=MnII and FeII , LN3O2R stands for pentadentate ligands) were prepared in order to assess the influence of the bending of the apical M-N≡C linkages on the magnetic anisotropy of the FeII derivatives and in turn on their Single-Molecule Magnet (SMM) behaviors. The cyanido-bridged [Cr2 M] derivatives were obtained by assembling trans-dicyanido CrIII complex [Cr(LN3O2Ph )(CN)2 ]- and divalent pentagonal bipyramid complexes [MII (H2 LN3O2R )]2+ with various R substituents (R=NH2 , cyclohexyl, S,S-mandelic) imparting different steric demand to the central moiety of the complexes. A comparative examination of the structural and magnetic properties showed an obvious effect of the deviation from straightness of the M-N≡C alignment on the slow relaxation of the magnetization exhibited by the [Cr2 Fe] complexes. Theoretical calculations have highlighted important effects of the bending of the apical C-N-Fe linkages on both the magnetic anisotropy of the FeII center and the exchange interactions with the CrIII units.

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

The Dzyaloshinskii-Moriya interaction is expected to be at the origin of interesting magnetic properties, such as multiferroicity, skyrmionic states, and exotic spin orders. Despite this, its theoretical determination is far from being established, neither from the point of view of ab initio methodologies nor from that of the extraction technique to be used afterward. Recently, a very efficient way to increase its amplitude has been demonstrated near the first-order spin-orbit coupling regime. Within the first-order regime, the anisotropic spin Hamiltonian involving the Dzyaloshinskii-Moriya operator becomes inappropriate. Nevertheless, in order to approach this regime and identify the spin Hamiltonian limitations, it is necessary to characterize the underlying physics. To this end, we have developed a simple electronic and spin-orbit model describing the first-order regime and used ab initio calculations to conduct a thorough methodological study.

How to create giant Dzyaloshinskii-Moriya interactions? Analytical derivation and ab initio calculations on model dicopper(II) complexes.

This paper is a theoretical "proof of concept" on how the on-site first-order spin-orbit coupling (SOC) can generate giant Dzyaloshinskii-Moriya interactions in binuclear transition metal complexes. This effective interaction plays a key role in strongly correlated materials, skyrmions, multiferroics, and molecular magnets of promising use in quantum information science and computing. Despite this, its determination from both theory and experiment is still in its infancy and existing systems usually exhibit very tiny magnitudes. We derive analytical formulas that perfectly reproduce both the nature and the magnitude of the Dzyaloshinskii-Moriya interaction calculated using state-of-the-art ab initio calculations performed on model bicopper(II) complexes. We also study which geometrical structures/ligand-field forces would enable one to control the magnitude and the orientation of the Dzyaloshinskii-Moriya vector in order to guide future synthesis of molecules or materials. This article provides an understanding of its microscopic origin and proposes recipes to increase its magnitude. We show that (i) the on-site mixings of 3d orbitals rule the orientation and magnitude of this interaction, (ii) increased values can be obtained by choosing more covalent complexes, and (iii) huge values (∼1000 cm-1) and controlled orientations could be reached by approaching structures exhibiting on-site first-order SOC, i.e., displaying an "unquenched orbital momentum."

Cyano-Bridged Fe(II)-Cr(III) Single-Chain Magnet Based on Pentagonal Bipyramid Units: On the Added Value of Aligned Axial Anisotropy.

A cyano-bridged Fe(II)-Cr(III) single-chain magnet designed to ensure a parallel orientation of the axial anisotropy of the building units is reported. This ferromagnetic chain compound consists of a pentagonal bipyramid Fe(II) complex with Ising-type anisotropy and a dicyanide Cr(III) complex interlinked through their apical positions. It is characterized by an energy gap for the magnetization reversal of Δeff/ kB = 113 K and exhibits magnetic hysteresis with a coercive field of 1400 Oe at 2 K which positions this compound among the very few examples of SCMs with spin reversal barriers above 100 K. The quite remarkable performances of this single-strand SCM are attributed to the alignment of the local anisotropy axes, which is supported by ab initio modeling. A discrete Cr2Fe complex based on the same building units and behaving as a SMM in zero field is also reported.

Pentagonal Bipyramid FeII Complexes: Robust Ising-Spin Units towards Heteropolynuclear Nanomagnets.

Pentagonal bipyramid FeII complexes have been investigated to evaluate their potential as Ising-spin building units for the preparation of heteropolynuclear complexes that are likely to behave as single-molecule magnets (SMMs). The considered monometallic complexes were prepared from the association of a divalent metal ion with pentadentate ligands that have a 2,6-diacetylpyridine bis(hydrazone) core (H2 LN3O2R ). Their magnetic anisotropy was established by magnetometry to reveal their zero-field splitting (ZFS) parameter D, which ranged between -4 and -13 cm-1 and was found to be modulated by the apical ligands (ROH versus Cl). The alteration of the D value by N-bound axial CN ligands, upon association with cyanometallates, was also assessed for heptacoordinated FeII as well as for related NiII and CoII derivatives. In all cases, N-coordinated cyanide ligands led to large magnetic anisotropy (i.e., -8 to -18 cm-1 for Fe and Ni, +33 cm-1 for Co). Ab initio calculations were performed on three FeII complexes, which enabled one to rationalize the role of the ligand on the nature and magnitude of the magnetic anisotropy. Starting from the pre-existing heptacoordinated complexes, a series of pentanuclear compounds were obtained by reactions with paramagnetic [W(CN)8 ]3- . Magnetic studies revealed the occurrence of ferromagnetic interactions between the spin carriers in all the heterometallic systems. Field-induced slow magnetic relaxation was observed for mononuclear FeII complexes (Ueff /kB up to 53 K (37 cm-1 ), τ0 =5×10-9  s), and SMM behavior was evidenced for a heteronuclear [Fe3 W2 ] derivative (Ueff /kB =35 K and τ0 =4.6 10-10  s), which confirmed that the parent complexes were robust Ising-type building units. High-field EPR spectroscopic investigation of the ZFS parameters for a Ni derivative is also reported.

From Positive to Negative Zero-Field Splitting in a Series of Strongly Magnetically Anisotropic Mononuclear Metal Complexes.

A series of mononuclear [M(hfa)2(pic)2] (Hhfa = 1,1,1,5,5,5-hexafluoro-2,4-pentanedione; pic = 4-methylpyridine; M = FeII, CoII, NiII, ZnII) compounds were obtained and characterized. The structures of the complexes have been resolved by single-crystal X-ray diffraction, indicating that, apart from the zinc derivative, the complexes are in a trans configuration. Moreover, a dramatic lenghthening of the Fe-N distances was observed, whereas the nickel(II) complex is almost perfectly octahedral. The magnetic anisotropy of these complexes was thoroughly studied by direct-current (dc) magnetic measurements, high-field electron paramagnetic resonance, and infrared (IR) magnetospectroscopy: the iron(II) derivative exhibits an out-of-plane anisotropy (DFe = -7.28 cm-1) with a high rhombicity, whereas the cobalt(II) and nickel(II) complexes show in-plane anisotropy (DCo ∼ 92-95 cm-1; DNi = 4.920 cm-1). Ab initio calculations were performed to rationalize the evolution of the structure and identify the excited states governing the magnetic anisotropy along the series. For the iron(II) complex, an out-of-phase alternating-current (ac) magnetic susceptibility signal was observed using a 0.1 T dc field. For the cobalt(II) derivative, the ac magnetic susceptibility shows the presence of two field-dependent relaxation phenomena: at low field (500 Oe), the relaxation process is beyond single-ion behavior, whereas at high field (2000 Oe), the relaxation of magnetization implies several mechanisms including an Orbach process with Ueff = 25 K and quantum tunneling of magnetization. The observation by µ-SQUID magnetization measurements of hysteresis loops of up to 1 K confirmed the single-ion-magnet behavior of the cobalt(II) derivative.

Ising-type Magnetic Anisotropy and Slow Relaxation of the Magnetization in Four-Coordinate Amido-Pyridine FeII Complexes.

A family of four-coordinate FeII complexes formed with N,N'-chelating amido-pyridine ligands was synthesized, and their magnetic properties were investigated. These distorted tetrahedral complexes exhibit significant magnetic anisotropy with zero-field splitting parameter D ranging between -17 and -12 cm-1. Ab initio calculations enabled identification of the structural factors that control the nature of the magnetic anisotropy and the rationalization of the variation of D in these complexes. It is shown that a reduced N-Fe-N angle involving the chelating nitrogen atoms of the ligands is at the origin of the negative D value and that the torsion between the two N-Fe-N planes imposed by steric hindrances further increases the |D| value. Field-induced slow relaxation of magnetization was observed for the three compounds, and a single-molecule magnet behavior with an energy barrier for magnetization flipping (Ueff) of 27 cm-1 could be evidenced for one of them.

An extended DFTB-CI model for charge-transfer excited states in cationic molecular clusters: model studies versus ab initio calculations in small PAH clusters.

We present an extension of the constrained density functional tight binding scheme combined with configuration interaction (DFTB-CI) to efficiently compute excited states of molecular cluster cations and their oscillator strengths from the ground state. The present extension consists of generalizing the initial model, relying on configurations with holes in the monomer HOMOs only, to configurations involving sub-HOMO holes, allowing for the description of higher excited states. The extended scheme is benchmarked on selected energy pathways with respect to available ab initio and new CASPT2 reference calculations on the benzene, naphthalene and pyrene dimer cations. The ability of the model to describe the potential energy surfaces and the transition dipole moments is discussed. The vertical electronic absorption spectra of the three dimer cations are calculated and compared with the theoretical litterature and available experimental data. Finally, the electronic absorption spectra of low energy isomers of the trimer and tetramer pyrene cluster cations are also predicted.

Electrically switchable magnetic molecules: inducing a magnetic coupling by means of an external electric field in a mixed-valence polyoxovanadate cluster.

Herein we evaluate the influence of an electric field on the coupling of two delocalized electrons in the mixed-valence polyoxometalate (POM) [GeV14 O40 ](8-) (in short V14 ) by using both a t-J model Hamiltonian and DFT calculations. In absence of an electric field the compound is paramagnetic, because the two electrons are localized on different parts of the POM. When an electric field is applied, an abrupt change of the magnetic coupling between the two delocalized electrons can be induced. Indeed, the field forces the two electrons to localize on nearest-neighbors metal centers, leading to a very strong antiferromagnetic coupling. Both theoretical approaches have led to similar results, emphasizing that the sharp spin transition induced by the electric field in the V14 system is a robust phenomenon, intramolecular in nature, and barely influenced by small changes on the external structure.

In search of a rational dressing of intermediate effective Hamiltonians.

The intermediate effective Hamiltonians are designed to provide M exact energies and the components of the corresponding eigenvectors in the N-dimensional model space, with N > M. The effective Hamiltonian is not entirely defined by these N × M conditions, and several dressings of the Hamiltonian matrix in the model space are possible. Some of them lead to unreliable N - M roots associated with the intermediate model space. This defect appears dramatically when one refers to the weak separability property, namely, the fact that in a noninteracting A···B problem where the model space only involves excitations on A, the consideration of the excitations on B should not affect the spectrum of A. We suggest variants that should maintain the physical meaning of the intermediate roots. Numerical comparisons illustrate the relevance of this proposal.

Singly occupied MOs in mono- and diradical conjugated hydrocarbons: comparison between variational single-reference, π-fully correlated and Hückel descriptions.

This work compares three descriptions of the unpaired electrons distribution in conjugated monoradical and diradical hydrocarbons involving one or two methylene groups attached to an aromatic skeleton. The first one is the simple Hückel topological Hamiltonian, the singly occupied molecular orbitals (SOMO) of which may be analytically obtained. The second one is the restricted open-shell self-consistent field (ROHF-SCF) method. The so-obtained distribution of the unpaired electrons on the skeleton appears deeply different from that predicted by the Hückel Hamiltonian, being more strongly localized on the external methylene groups. More elaborate methods treat all π electrons in the π valence molecular orbitals (MOs) through a full valence π complete active space self-consistent field (CASSCF) treatment. The distributions of the unpaired electrons (given by the natural MOs of occupation number close to 1) are surprisingly similar to those predicted by the Hückel model. The spin density distributions, including spin polarization effects, can be improved by further configuration interactions involving one hole-one particle excitations and compared with the experimental hyperfine coupling constant ratios. This comparison confirms the lack of delocalization of the magnetic orbitals defined from the self-consistent single-reference treatment. We show that, provided correct SOMO are used, a single excitation CI performed on top of a single reference gives accurate spin densities. Finally, a rationalization of the role of the dynamic correlation in correcting the excessive localization of the unpaired electron(s) at the ROHF level on the exocyclic methylene group(s) is given, attributing it to the dynamic charge polarization of the charge transfer configurations between methylene and the aromatic frame.

Rationalization of the behavior of M2(CH3CS2)4I (M = Ni, Pt) chains at room temperature from periodic density functional theory and ab initio cluster calculations.

The electrical conductivities and plausible charge-ordering states in the room temperature (r.t.) phase for MMX chains [Ni(2)(dta)(4)I](∞) and [Pt(2)(dta)(4)I](∞) (dta = CH(3)CS(2)(-)) have been analyzed with periodic density functional theory (DFT) and correlated ab initio calculations combined with the effective Hamiltonian theory. Periodic DFT calculations show a more delocalized nature of the ground state in [Pt(2)(dta)(4)I](∞) compared to [Ni(2)(dta)(4)I](∞), which features a rather large energy gap between the occupied and empty bands, and charge polarized dimer units. A larger electrical conductivity for the Pt chain can be expected, especially because the Fermi level lies within a band with contributions from Pt and I orbitals. Electronic structure parameters extracted from ab initio cluster calculations show that the large difference between the observed conductivities at 300 K for Ni and Pt compounds, of 3 orders of magnitude, cannot be explained from the parameters extracted from an embedded M(2)(dta)(4)I(2) dimer fragment alone. When tetramer fragments are considered, we observe that the interdimer transfer integral (t) between neighboring M(2) units connected by an iodine atom at correlated level is comparable in both chains. On the other hand, the energy to transfer an electron from a dimer to the neighboring one (Coulomb repulsion U) is three times larger in the Ni compound with respect to the Pt chain, in line with the poor conductivity of the former. The electronic structure of the M(4)(dta)(8)I(3) fragment points to an alternate charge-polarization state for Ni and an average valence state for Pt when the r.t. X-ray structure is considered.

Ab initio study of the influence of structural parameters on the potential energy surfaces of spin-crossover Fe(II) model compounds.

In spin-crossover (SCO) compounds exhibiting a light induced excited spin state trapping (LIESST) effect, the thermodynamic T(1∕2) and kinetic T(LIESST) temperature values depend on the features of the potential energy surfaces (PES) of the two lowest singlet and quintet states but also on vibrational contributions, collective effects, such as electrostatics, for instance, spin-orbit couplings to a lesser extent, etc. In this work, the question of the link between the shape of the PES of SCO compounds exhibiting a LIESST effect and their first coordination sphere structure is addressed from wave function theory based ab initio calculations. Fe(II) complexes based on model ligands suited to reproduce the main characteristics of the PES of such compounds are distorted to emphasize selectively the role played by the metal-ligand distances and the ligand-metal-ligand angles. The studied angular deformations are those usually observed in many Fe(L)(2)(NCS)(2) complexes. It is shown that the larger the deformation between the low spin and high spin equilibrium geometries, the higher the energy barrier from the high spin state and the weaker the energy difference between the bottom of the wells. These results corroborate observations made by experimentalists on a large number of complexes. While the PES features only constitutes one of the contributions to these temperatures, it is worth noticing that, relating T(1∕2) to the energy difference between the bottoms of the singlet and quintet wells and the T(LIESST) to the energy barrier from the quintet bottom well, the same slope of the empirical law T(LIESST) = -0.3T(1∕2)+T(0) is observed.

Designing magnetic organic lattices from high-spin polycyclic units.

This work addresses the conception of purely organic magnetic materials by properly bridging high-spin polycyclic hydrocarbons A and B, through covalent ligands L. The strategy varies two degrees of freedom that govern the magnetic character of the A-L--B sequence, namely, the bridge response to spin polarization and the relative signs of spin density on carbon atoms to which the bridge is attached. Topological prescriptions based on Ovchinnikov's rule are proposed to predict ground-state spin multiplicities of various A-L-B sets. The relevance of these guiding principles is essentially confirmed through DFT calculations on dimers connected by conjugated bridges. The transferability of interunit magnetic couplings to larger assemblies is further checked, the building blocks tending to maintain their high-spin character whatever the environment. Such local designs open the way to periodic lattices of ferromagnetic, antiferromagnetic, ferrimagnetic, or paramagnetic materials.

Chemical tuning of spin clock transitions in molecular monomers based on nuclear spin-free Ni(ii).

We report the existence of a sizeable quantum tunnelling splitting between the two lowest electronic spin levels of mononuclear Ni complexes. The level anti-crossing, or magnetic "clock transition", associated with this gap has been directly monitored by heat capacity experiments. The comparison of these results with those obtained for a Co derivative, for which tunnelling is forbidden by symmetry, shows that the clock transition leads to an effective suppression of intermolecular spin-spin interactions. In addition, we show that the quantum tunnelling splitting admits a chemical tuning via the modification of the ligand shell that determines the crystal field and the magnetic anisotropy. These properties are crucial to realize model spin qubits that combine the necessary resilience against decoherence, a proper interfacing with other qubits and with the control circuitry and the ability to initialize them by cooling.

Theoretical design of high-spin polycyclic hydrocarbons.

High-spin organic structures can be obtained from fused polycyclic hydrocarbons, by converting selected peripheral HC(sp(2)) sites into H(2)C(sp(3)) ones, guided by Ovchinnikov's rule. Theoretical investigation is performed on a few examples of such systems, involving three to twelve fused rings, and maintaining threefold symmetry. Unrestricted DFT (UDFT) calculations, including geometry optimizations, confirm the high-spin multiplicity of the ground state. Spin-density distributions and low-energy spectra are further studied through geometry-dependent Heisenberg-Hamiltonian diagonalizations and explicit correlated ab initio treatments, which all agree on the high-spin character of the suggested structures, and locate the low-lying states at significantly higher energies. In particular, the lowest-lying state of lower multiplicity is always found to be higher than kT at room temperature (at least ten times higher). Simplification of the ferromagnetic organization based on sets of semilocalized nonbonding orbitals is proposed. Molecular architectures are thus conceived in which the ferromagnetically-coupled unpaired electrons tally up to one third of the involved conjugated carbons. Connecting such building blocks should provide bidimensional materials endowed with robust magnetic properties.