Dinitrogen Activation Mediated by the (P2PPh)Fe Complex: Electronic Structure, Dimerization Mechanism, and Magnetic Coupling.

Herein, we report the estimation of the extent of dinitrogen activation by different charged and structural forms of (P2PPh)Fe biomimetic catalysts, which, in the presence of light, exhibit significant yield in the N2-to-NH3 conversion. Complete active space self-consistent field (CASSCF) calculations have been used to determine the electronic structure of different reduced forms of the mononuclear complexes: the neutral (P2PPh)Fe(N2)2 adduct and the anionic [(P2PPh)Fe(N2)]- and [(P2PPh)Fe(N2)]2- complexes. These calculations also revealed that the extent of reduction of a dinitrogen molecule reaches up to one electron (N21-) due to the back-bonding from the Fe center, in agreement with the changes observed in the vibration frequency of the N-N bond in these complexes. In addition, the energy profile of the dimerization of the mononuclear (P2PPh)Fe(N2)2 complex to the dinuclear mono-N2-bridged [(P2PPh)Fe]2(µ-N2) complex has been determined by means of density functional theory (DFT) calculations. A three-step mechanism has been proposed for the dimerization, favored by both kinetics and thermodynamics criteria. Finally, the magnetic coupling constant in the diiron (µ-N2) complex is estimated from CASSCF/NEVPT2 calculations. Such a dinuclear complex presents a strong antiferromagnetic coupling resulting from the interaction between two S = 1 d6 Fe2+ ions, bridged by a highly activated dinitrogen molecule (N22-) with two electrons on the π* orbitals.

How complex-surface interactions modulate the spin transition of Fe(II) SCO complexes supported on metallic surfaces?

The deposition of a prototypical spin-crossover [Fe(phen)2(NCS)2] complex on Au(111), Cu(111) and Ag(111) surfaces has been investigated by means of periodic DFT+U calculations, with the aim of understanding how different metallic surfaces affect the spin state switching. Our results show that adsorption is metal- and spin-dependent, with different preferred adsorption sites for the different surfaces and spin states. For the three considered surfaces adsorption energies are larger in the LS state than in the HS one, which increases the transition enthalpy by 58.7 kJ mol-1 for Cu(111), 14.6 kJ mol-1 for Au(111) and 9.6 kJ mol-1 for Ag(111) with respect to the free molecule. There is a clear correlation between this effect and the extent of the complex-surface interaction, which can be established from adsorption energies, surface-complex distances and charge density difference plots as: Cu(111) > Au(111) > Ag(111). Therefore, a stronger interaction with the surface produces a larger energy difference between two spin states, making the spin transition less probable to occur. Finally, our calculations show that it would be possible to probe the spin-state of the deposited molecules from the STM images, in line with the recent experimental results.

Deposition of the Spin Crossover FeII -Pyrazolylborate Complex on Au(111) Surface at the Molecular Level.

The interaction at the molecular level of the spin-crossover (SCO) FeII ((3,5-(CH3 )2 Pz)3 BH)2 complex with the Au(111) surface is analyzed by means of rPBE periodic calculations. Our results show that the adsorption on the metallic surface enhances the transition energy, increasing the relative stability of the low spin (LS) state. The interaction indeed is spin-dependent, stronger for the low spin than the high spin (HS) state. The different strength of the Fe ligand field at low and high temperature manifests on the nature, spatial extension and relative energy of the states close to the Fermi level, with a larger metal-ligand hybridization in the LS state. This feature is of relevance for the differential adsorption of the LS and HS molecules, the spin-dependent conductance, and for the differences found in the corresponding STM images, correctly reproduced from the density of states provided by the rPBE calculations. It is expected that this spin dependence will be a general feature of the SCO molecule-substrate interaction, since it is rooted in the different ligand field of Fe site at low and high temperatures, a common hallmark of the FeII SCO complexes. Finally, the states involved in the LIESST phenomenon has been identified through NEVPT2 calculations on a model reaction path. A tentative pathway for the photoinduced LS→HS transition is proposed, that does not involve the intermediate triplet states, and nicely reproduces both the blue laser wavelength required for the activation, and the wavelength of the reverse HS → LS transition.

Theoretical inspection of the spin-crossover [Fe(tzpy)2(NCS)2] complex on Au(100) surface.

We explore the deposition of the spin-crossover [Fe(tzpy)2(NCS)2] complex on the Au(100) surface by means of density functional theory (DFT) based calculations. Two different routes have been employed: low-cost finite cluster-based calculations, where both the Fe complex and the surface are maintained fixed while the molecule approaches the surface; and periodic DFT plane-wave calculations, where the surface is represented by a four-layer slab and both the molecule and surface are relaxed. Our results show that the bridge adsorption site is preferred over the on-top and fourfold hollow ones for both spin states, although they are energetically close. The LS molecule is stabilized by the surface, and the HS-LS energy difference is enhanced by about 15%-25% once deposited. The different Fe ligand field for LS and HS molecules manifests on the composition and energy of the low-lying bands. Our simulated STM images indicate that it is possible to distinguish the spin state of the deposited molecules by tuning the bias voltage of the STM tip. Finally, it should be noted that the use of a reduced size cluster to simulate the Au(100) surface proves to be a low-cost and reliable strategy, providing results in good agreement with those resulting from state-of-the-art periodic calculations for this system.

Light-Induced Control of the Spin Distribution on Cu⁻Dithiolene Complexes: A Correlated Ab Initio Study.

Metal dithiolene complexes-M(dmit)2-are key building blocks for magnetic, conducting, and optical molecular materials, with singular electronic structures resulting from the mixing of the metal and dmit ligand orbitals. Their use in the design of magnetic and conducting materials is linked to the control of the unpaired electrons and their localized/delocalized nature. It has been recently found that UV⁻Vis light can control the spin distribution of some [Cu(dmit)2]-2 salts in a direct and reversible way. In this work, we study the optical response of these salts and the origin of the differences observed in the EPR spectra under UV⁻Vis irradiation by means of wave function-based quantum chemistry methods. The low-lying states of the complex have been characterized and the electronic transitions with a non-negligible oscillator strength have been identified. The population of the corresponding excited states promoted by the UV⁻Vis absorption produces significant changes in the spin distribution, and could explain the changes observed in the system upon illumination. The interaction between neighbor [Cu(dmit)2]-2 complexes is weakly ferromagnetic, consistent with the relative orientation of the magnetic orbitals and the crystal packing, but in disagreement with previous assignments. Our results put in evidence the complex electronic structure of the [Cu(dmit)2]-2 radical and the relevance of a multideterminantal approach for an adequate analysis of their properties.

Light-Induced Spin Transitions in Copper-Nitroxide-Based Switchable Molecular Magnets: Insights from Periodic DFT+U Calculations.

The electronic structure and magnetic interactions of three members of the breathing crystal Cu(hfac)2 LR family (hfac=hexafluoroacetylacetonato, LR =pyrazole-substituted nitronyl nitroxides with R=Me, Et, Pr, iPr, Bu ), mainly Cu(hfac)2 LPr (1), Cu(hfac)2 LBu ⋅0.5 C8 H18 (2) and Cu(hfac)2 LBu ⋅0.5 C8 H10 (3), have been analyzed by means of periodic plane-wave based DFT+U calculations. These CuII -nitroxide based molecular magnets display thermally and optically induced switchable behavior and light-induced excited spin state trapping phenomena. The calculations confirm the presence of temperature-dependent exchange interaction within the spin triads formed by the nitroxide-copper(II)-nitroxide units, in line with the changes observed in the effective magnetic moment. Moreover, they quantify the interchain interaction mediated by the terminal nitroxide group of two spin triads in neighboring polymer chains. This interaction competes with the exchange interaction within the spin triads at high temperature, and introduces 1D exchange channels that do not coincide with the polymeric chains. The density of states reveal that the low-lying conduction states potentially involved in the UV/Vis transitions are located on the nitroxide radicals, the hfac groups and the Cu atoms. Then, the density of states is almost independent of the solvent and the R group. This suggests the possibility of light-induced spin switching for other members of this family. The 500 nm band of the low-temperature phase can be ascribed to a ligand-to-metal charge transfer transition between the nitroxide and Cu bands.

Electronic Structure and Magnetic Interactions in the Radical Salt [BEDT-TTF]2[CuCl4].

The magnetic behavior and electric properties of the hybrid radical salt [BEDT-TTF]2[CuCl4] have been revisited through extended experimental analyses and DDCI and periodic DFT plane waves calculations. Single crystal X-ray diffraction data have been collected at different temperatures, discovering a phase transition occurring in the 250-300 K range. The calculations indicate the presence of intradimer, interdimer, and organic-inorganic π-d interactions in the crystal, a magnetic pattern much more complex than the Bleaney-Bowers model initially assigned to this material. Although this simple model was good enough to reproduce the magnetic susceptibility data, our calculations demonstrate that the actual magnetic structure is significantly more intricate, with alternating antiferromagnetic 1D chains of the organic BEDT-TTF+ radical, connected through weak antiferromagnetic interactions with the CuCl42- ions. Combination of experiment and theory allowed us to unambiguously determine and quantify the leading magnetic interactions in the system. The density-of-states curves confirm the semiconductor nature of the system and the dominant organic contribution of the valence and conduction band edges. This general and combined approach appears to be fundamental in order to properly understand the magnetic structure of these complex materials, where experimental data can actually be fitted from a variety of models and parameters.

Evaluation of the Giant Ferromagnetic π-d Interaction in Iron-Phthalocyanine Molecule.

The interaction between itinerant π and localized d electrons in metal-phthalocyanines, namely, Jπd interaction, is considered as responsible for the giant negative magnetoresistance observed in several phthalocyanine-based conductors, among many other important physical properties. Despite the fundamental and technological importance of this on-site intramolecular interaction, its giant ferromagnetic nature has been only recently demonstrated by the experiments conducted by Murakawa et al. in the neutral radical [Fe(Pc)(CN)2]·2CHCl3 ( Phys. Rev. B 2015 , 92 , 054429 ). In this article, we present the theoretical evaluation of this interaction combining wave function-based electronic calculations on isolated Fe(Pc)(CN)2 molecules and density functional theory-based periodic calculations on the crystal. Our calculations confirm the ferromagnetic nature of the π-d interaction, with a coupling constant as large as Jπd/kB = 570 K, in excellent agreement with the experiments, and the presence of intermolecular antiferromagnetic interactions driven by the π-π overlap of neighboring phthalocyaninato molecules. The analysis of the wave function of the ground state of the Fe(Pc)(CN)2 molecule provides the clues of the origin of this giant ferromagnetic π-d interaction.

Analysis of the Magnetic Exchange Interactions in Yttrium(III) Complexes Containing Nitronyl Nitroxide Radicals.

We report a combined theoretical and experimental investigation of the exchange interactions governing the magnetic behavior of a series of nitronyl nitroxide (NIT)-based Y(III) complexes, i.e., Y(hfac)3(NIT-R)2 with R = PhOPh (1), PhOEt (2), and PhOMe (3a, 3b). Even though some of these complexes or their Dy(III) parents were previously described in the literature [ Zhao et al. Transition Met. Chem. 2006 , 31 , 593 ; Bernot et al. J. Am. Chem. Soc. 2009 , 131 , 5573 ], their synthesis procedure as well as their structural and magnetic properties were completely reconsidered. Depending on the nature of R and the crystallization conditions, Y(hfac)3(NIT-R)2 units can be organized as supramolecular dimers or linear or orthogonal chains. Such structural diversity within the series induces extremely different magnetic behaviors. The observed behaviors are rationalized by state-of-the-art wave function-based quantum-chemical approaches (CASSCF/DDCI) that demonstrate the existence of not only effective intramolecular interactions between the NIT-R radical ligands of an isolated Y(hfac)3(NIT-R)2 molecule but also intermolecular interactions between NIT-R moieties belonging to different Y(hfac)3(NIT-R)2 units. These results are supported by the use of spin Hamiltonian models going beyond the basic Bleaney-Bowers formalism to properly fit the experimental magnetic data. Finally, the microscopic mechanisms behind the evidenced intramolecular exchange interactions are elucidated through the inspection of the calculated wave functions. In particular, whereas the role of Y orbitals was already proposed, we herein demonstrate the contribution of the hfac- ancillary ligands in mediating the magnetic interactions between the NIT radicals.

Highly efficient perturbative + variational strategy based on orthogonal valence bond theory for the evaluation of magnetic coupling constants. Application to the trinuclear Cu(ii) site of multicopper oxidases.

A new strategy based on orthogonal valence-bond analysis of the wave function combined with intermediate Hamiltonian theory has been applied to the evaluation of the magnetic coupling constants in two AF systems. This approach provides both a quantitative estimate of the J value and a detailed analysis of the main physical mechanisms controlling the coupling, using a combined perturbative + variational scheme. The procedure requires a selection of the dominant excitations to be treated variationally. Two methods have been employed: a brute-force selection, using a logic similar to that of the CIPSI approach, or entanglement measures, which identify the most interacting orbitals in the system. Once a reduced set of excitations (about 300 determinants) is established, the interaction matrix is dressed at the second-order of perturbation by the remaining excitations of the CI space. The diagonalization of the dressed matrix provides J values in good agreement with experimental ones, at a very low-cost. This approach demonstrates the key role of d → d* excitations in the quantitative description of the magnetic coupling, as well as the importance of using an extended active space, including the bridging ligand orbitals, for the binuclear model of the intermediates of multicopper oxidases. The method is a promising tool for dealing with complex systems containing several active centers, as an alternative to both pure variational and DFT approaches.

Metal-Metal Interactions in Trinuclear Copper(II) Complexes [Cu3(RCOO)4(H2TEA)2] and Binuclear [Cu2(RCOO)2(H2TEA)2]. Syntheses and Combined Structural, Magnetic, High-Field Electron Paramagnetic Resonance, and Theoretical Studies.

The trinuclear [Cu3(RCOO)4(H2TEA)2] copper(II) complexes, where RCOO(-) = 2-furoate (1), 2-methoxybenzoate (2), and 3-methoxybenzoate (3, 4), as well as dimeric species [Cu2(H2TEA)2(RCOO)2]·2H2O, have been prepared by adding triethanolamine (H3TEA) at ambient conditions to hydrated Cu(RCOO)2 salts. The newly synthesized complexes have been characterized by elemental analyses, spectroscopic techniques (IR and UV-visible), magnetic susceptibility, single crystal X-ray structure determination and theoretical calculations, using a Difference Dedicated Configuration Interaction approach for the evaluation of magnetic coupling constants. In 1 and 2, the central copper atom lies on an inversion center, while in the polymorphs 3 and 4, the three metal centers are crystallographically independent. The zero-field splitting parameters of the trimeric compounds, D and E, were derived from high-field, high-frequency electron paramagnetic resonance spectra at temperatures ranging from 3 to 290 K and were used for the interpretation of the magnetic data. It was found that the dominant interaction between the terminal and central Cu sites J12 is ferromagnetic in nature in all complexes, even though differences have been found between the symmetrical or quasi-symmetrical complexes 1-3 and non-symmetrical complex 4, while the interaction between the terminal centers, J23, is negligible.

Mechanism of Magnetostructural Transitions in Copper-Nitroxide-Based Switchable Molecular Magnets: Insights from ab Initio Quantum Chemistry Calculations.

The gradual magnetostructural transition in breathing crystals based on copper(II) and pyrazolyl-substituted nitronyl nitroxides has been analyzed by means of DDCI quantum chemistry calculations. The magnetic coupling constants (J) within the spin triads of Cu(hfac)2L(Bu)·0.5C8H18 have been evaluated for the X-ray structures reported at different temperatures. The coupling is strongly antiferromagnetic at low temperature and becomes ferromagnetic when the temperature increases. The intercluster magnetic coupling (J') is antiferromagnetic and shows a marked dependence on temperature. The magnetostructural transition can be reproduced using the calculated J values for each structure in the simulation of the magnetic susceptibility. However, the µ(T) curve can be improved nicely by considering the coexistence of two phases in the transition region, whose ratio varies with temperature corresponding to both the weakly and strongly coupled spin states. These results complement a recent VT-FTIR study on the parent Cu(hfac)2L(Pr) compound with a gradual magnetostructural transition.

Electronic structure aspects of the complete O2 transfer reaction between Ni(II) and Mn(II) complexes with cyclam ligands.

This work explores the electronic structure aspects involving the complete intermolecular O2 transfer between Ni(ii) and Mn(ii) complexes, both containing N-tetramethylated cyclams (TMC). The energy of the low-lying states of reactants, intermediates and products is established at the CASSCF level and also the DDCI level when possible. The orthogonal valence bond analysis of the wave functions obtained from CASSCF and DDCI calculations indicates the dominant superoxide nature of all the adducts participating in the reaction, and consequently that the whole reaction can be described as the transfer of the superoxide O2(-) between Ni(ii) and Mn(ii) complexes, without any additional change in the electronic structure of the fragments.

Improving the calculation of magnetic coupling constants in MRPT methods.

The magnetic coupling in transition metal compounds with more than one unpaired electron per magnetic center has been studied with multiconfigurational perturbation theory. The usual shortcomings of these methodologies (severe underestimation of the magnetic coupling) have been overcome by describing the Slater determinants with a set of molecular orbitals that maximally resemble the natural orbitals of a high-level multiconfigurational reference configuration interaction calculation. These orbitals have significant delocalization tails onto the bridging ligands and largely increase the coupling strengths in the perturbative calculation.

Magnetic coupling constants of self-assembled Cu(II) [3×3] grids: alternative spin model from theoretical calculations.

This paper reports a theoretical analysis of the electronic structure and magnetic properties of a ferromagnetic Cu(II) [3×3] grid. A two-step strategy, combining calculations on the whole grid and on binuclear fragments, has been employed to evaluate all the magnetic interactions in the grid. The calculations confirm an S = 7/2 ground state, which is in accordance with the magnetisation versus field curve and the thermal dependence of the magnetic moment data. Only the first-neighbour coupling terms present non-negligible amplitudes, all of them in agreement with the structure and arrangement of the Cu 3d magnetic orbitals. The results indicate that the dominant interaction in the system is the antiferromagnetic coupling between the ring and the central Cu sites (J3 = J4 ≈ -31 cm(-1)). In the ring two different interactions can be distinguished, J1 = 4.6 cm(-1) and J2 = -0.1 cm(-1), in contrast to the single J model employed in the magnetic data fit. The calculated J values have been used to determine the energy level distribution of the Heisenberg magnetic states. The effective magnetic moment versus temperature plot resulting from this ab initio energy profile is in good agreement with the experimental curve and the fitting obtained with the simplified spin model, despite the differences between these two spin models. This study underlines the role that the theoretical evaluations of the coupling constants can play on the rationalisation of the magnetic properties of these complex polynuclear systems.

On the reaction mechanism of the complete intermolecular O2 transfer between mononuclear nickel and manganese complexes with macrocyclic ligands.

The recently described intermolecular O2 transfer between the side-on Ni-O2 complex [(12-TMC)Ni-O2](+) and the manganese complex [(14-TMC)Mn](2+), where 12-TMC and 14-TMC are 12- and 14-membered macrocyclic ligands, 12-TMC=1,4,7,10-tetramethyl-1,4,7,10-tetraazacyclododecane and 14-TMC=1,4,8,11-tetramethyl-1,4,8,11-tetraazacyclotetradecane, is studied by means of DFT methods. B3LYP calculations including long-range corrections and solvent effects are performed to elucidate the mechanism. The potential energy surfaces (PESs) compatible with different electronic states of the reactants have been analyzed. The calculations confirm a two-step reaction, with a first rate-determining bimolecular step and predict the exothermic character of the global process. The relative stability of the products and the reverse barrier are in line with the fact that no reverse reaction is experimentally observed. An intermediate with a µ-η(1):η(1)-O2 coordination and two transition states are identified on the triplet PES, slightly below the corresponding stationary points of the quintet PES, suggesting an intersystem crossing before the first transition state. The calculated activation parameters and the relative energies of the two transition sates and the products are in very good agreement with the experimental data. The calculations suggest that a superoxide anion is transferred during the reaction.

Voltage-induced modulation of the magnetic exchange in binuclear Fe(III) complex deposited on Au(111) surface.

We present a complete computational study devoted to the deposition of a magnetic binuclear complex on a metallic surface, aimed to obtain insight into the interaction of magnetically coupled complexes with their supporting substrates, as well as their response to external electrical stimuli applied through a surface-molecule-STM molecular junction-like architecture. Our results not only show that the deposition is favorable in two of the four studied orientations, but also, that the magnetic coupling is only slightly perturbed once the complex is adsorbed. We observe that the effects of the applied bias voltage on the magnetic coupling strongly depend on the molecule orientation with respect to the surface and the voltage polarity. Further analysis shows that this behavior is attributable to the stabilization/destabilization of the d-type singly occupied orbitals of the iron centers, reinforced by the strong local electric fields and induced charge densities only present in certain orientations of the deposited molecule and applied voltage polarity.

Correction: Computational demonstration of isomer- and spin-state-dependent charge transport in molecular junctions composed of charge-neutral iron(II) spin-crossover complexes.

Correction for 'Computational demonstration of isomer- and spin-state-dependent charge transport in molecular junctions composed of charge-neutral iron(II) spin-crossover complexes' by Nicolás Montenegro-Pohlhammer, et al., Dalton Trans., 2023, 52, 1229-1240, https://doi.org/10.1039/D2DT02598A.

On the controversial fitting of susceptibility curves of ferromagnetic Cu(II) cubanes: insights from theoretical calculations.

This paper reports a theoretical analysis of the electronic structure and magnetic properties of a tetranuclear Cu(II) complex, [Cu(4) (HL)(4)], which has a 4+2 cubane-like structure (H(3) L=N,N'-(2-hydroxypropane-1,3-diyl)bis(acetylacetoneimine)). These theoretical calculations indicate a quintet (S=2) ground state; the energy-level distribution of the magnetic states confirm Heisenberg behaviour and correspond to an S(4) spin-spin interaction model. The dominant interaction is the ferromagnetic coupling between the pseudo-dimeric units (J(1) =22.2 cm(-1)), whilst a weak and ferromagnetic interaction is found within the pseudo-dimeric units (J(2) =1.4 cm(-1)). The amplitude and sign of these interactions are consistent with the structure and arrangement of the magnetic Cu 3d orbitals; they accurately simulate the thermal dependence of magnetic susceptibility, but do not agree with the reported J values (J(1) =38.4 cm(-1), J(2) =-18.0 cm(-1)) that result from the experimental fitting. This result is not an isolated case; many other polynuclear systems, in particular 4+2 Cu(II) cubanes, have been reported in which the fitted magnetic terms are not consistent with the geometrical features of the system. In this context, theoretical evaluation can be considered as a valuable tool in the interpretation of the macroscopic behaviour, thus providing clues for a rational and directed design of new materials with specific properties.

Multidimensional network structures and versatile magnetic properties of intermolecular compounds of a radical-anion ligand, [1,2,5]thiadiazolo[3,4-f][1,10]phenanthroline 1,1-dioxide.

The crystal structures and magnetic properties of seven kinds of [1,2,5]thiadiazolo[3,4-f][1,10]phenanthroline 1,1-dioxide (tdapO2) radical-anion salts, namely, K·tdapO2, K·tdapO2·0.5MeCN, K·(tdapO2)2, Rb·(tdapO2)2, Cs7·(tdapO2)6·ClO4, (NH4)2·tdapO2·I, and Hppda·tdapO2·MeCN, were investigated. Single-crystal X-ray analyses of these radical-anion salts revealed formation of π-stacking columns and the presence of intercolumnar coordination bonding or hydrogen bonding. The intermolecular magnetic coupling constants in these salts range from strong antiferromagnetic (J/kB = -310 K) to ferromagnetic (J/kB = 24 K). Ab initio calculations performed on the nearest-neighbor radical pairs in the π-stacking columns suggested that the magnetic interactions are strongly governed by the overlap between the two anionic radical species and well explain the observed ferromagnetic and antiferromagnetic interactions. In addition, calculations of a hypothetical oxygen-less tdap analogue suggested that the presence of oxygen in tdapO2 significantly reduces the hopping integral and enhances the probability of ferromagnetic interaction.