First-principles calculations of zero-field splitting parameters

Charge and Spin Transfer Dynamics in a Weakly Coupled Porphyrin Dimer

The dynamics of electron and spin transfer in the radical cation and photogenerated triplet states of a tetramethylbiphenyl-linked zinc-porphyrin dimer were investigated, so as to test the relevant parameters for the design of a single-molecule spin valve and the creation of a novel platform for the photogeneration of high-multiplicity spin states. We used a combination of multiple techniques, including variable-temperature continuous wave EPR, pulsed proton electron-nuclear double resonance (ENDOR), transient EPR, and optical spectroscopy. The conclusions are further supported by density functional theory (DFT) calculations and comparison to reference compounds. The low-temperature cw-EPR and room-temperature near-IR spectra of the dimer monocation demonstrate that the radical cation is spatially localized on one side of the dimer at any point in time, not coherently delocalized over both porphyrin units. The EPR spectra at 298 K reveal rapid hopping of the radical spin density between both sites of the dimer via reversible intramolecular electron transfer. The hyperfine interactions are modulated by electron transfer and can be quantified using ENDOR spectroscopy. This allowed simulation of the variable-temperature cw-EPR spectra with a two-site exchange model and provided information on the temperature-dependence of the electron transfer rate. The electron transfer rates range from about 10.0 MHz at 200 K to about 53.9 MHz at 298 K. The activation enthalpies Δ‡H of the electron transfer were determined as Δ‡H = 9.55 kJ mol-1 and Δ‡H = 5.67 kJ mol-1 in a 1:1:1 solvent mixture of CD2Cl2/toluene-d8/THF-d8 and in 2-methyltetrahydrofuran, respectively, consistent with a Robin-Day class II mixed valence compound. These results indicate that the interporphyrin electronic coupling in a tetramethylbiphenyl-linked porphyrin dimer is suitable for the backbone of a single-molecule spin valve. Investigation of the spin density distribution of the photogenerated triplet state of the Zn-porphyrin dimer reveals localization of the triplet spin density on a nanosecond time scale on one-half of the dimer at 20 K in 2-methyltetrahydrofuran and at 250 K in a polyvinylcarbazole film. This establishes the porphyrin dimer as a molecular platform for the formation of a localized, photogenerated triplet state on one porphyrin unit that is coupled to a second redox-active, ground-state porphyrin unit, which can be explored for the formation of high-multiplicity spin states.

Control of the geometry and anisotropy driven by the combination of steric and anion coordination effects in CoII complexes with N6-tripodal ligands: the impact of the size of the ligand on the magnetization relaxation time

Four mononuclear CoII complexes of formula [Co(L)(SCN)2(CH3OH)0.5(H2O)0.5]·1.5H2O·0.75CH3OH (1), [Co(L1)Cl2]·H2O·2CH3CN (2), [Co(L1)(SCN)2]·1.5H2O·CH3OH (3) and [Co(L1)]ClO4·2CH3OH (4) were prepared from the N6-tripodal Schiff base ligands (S)P[N(Me)NC(H)2-Q]3 (L) and (S)P[N(Me)NC(H)1-ISOQ]3 (L1), where Q and ISOQ represent quinolyl and isoquinolyl moieties, respectively. In 1, the L ligand does not coordinate to the CoII ion in a tripodal manner but using a new N,N,S tridentate mode, which is due to the fact that the N6-tripodal coordination promotes a strong steric hindrance between the quinolyl moieties. However, L1 can coordinate to the CoII ions either in a tripodal manner using CoII salts with poorly coordinating anions to give 4 or in a bisbidentate fashion using CoII salt-containing medium to strongly coordinating anions to afford 2 and 3. In the case of L1, there is no steric hindrance between ISOQ moieties after coordination to the CoII ion. The CoII ion exhibits a distorted octahedral geometry for compounds 1-3, with the anions in cis positions for the former and in trans positions for the two latter compounds. Compound 4 shows an intermediate geometry between an octahedral and trigonal prism but closer to the latter one. DC magnetic properties, HFEPR and FIRMS measurements and ab initio calculations demonstrate that distorted octahedral complexes 1-3 exhibit easy-plane magnetic anisotropy (D > 0), whereas compound 4 shows large easy-axis magnetic anisotropy (D < 0). Comparative analysis of the magneto-structural data underlines the important role that is played not only by the coordination geometry but also the electronic effects in determining the anisotropy of the CoII ions. Compounds 2-3 show a field-induced slow relaxation of magnetization. Despite its large easy-axis magnetic anisotropy, compound 4 does not show significant slow relaxation (SMR) above 2 K under zero applied magnetic fields, but its magnetic dilution with ZnII triggers SMR at zero field. Finally, it is worth remarking that compounds 2-4 show smaller relaxation times than the analogous complexes with the tripodal ligand bearing in its arms pyridine instead of isoquinoline moieties, which is most likely due to the increase of the molecular size in the former one.

Structural, optical and magnetic properties of a new metal-organic CoII-based complex

A mononuclear cobalt(ii) complex [C5H8N3]2[CoCl4(C5H7N3)2] (I) was synthesized and structurally characterized. Single crystal X-ray diffraction analysis indicates that monometallic Co(ii) ions acted as coordination nodes in a distorted octahedral geometry, giving rise to a supramolecular architecture. The latter is made up of a ½ unit form composed of an anionic element [Co0.5Cl2(C5H7N3)]- and one 2-amino-4-methylpyrimidinium cation [C5H8N3]+. The crystalline arrangement of this compound adopts the sandwich form where inorganic parts are sandwiched between the organic sheets following the [100] direction. More information regarding the structure hierarchy has been supplied based on Hirshfeld surface analysis; the X⋯H (X = N, Cl) interactions play a crucial role in stabilizing the self-assembly process of I, complemented by the intervention of π⋯π electrostatic interaction created between organic entities. Thermal analyses were carried out to study the thermal behavior process. Static magnetic measurements and ab initio calculations of compound I revealed the easy-axis anisotropy character of the central Co(ii) ion. Two-channel field-induced slow-magnetic relaxation was observed; the high-frequency channel is characterized by underbarrier relaxation with U eff = 16.5 cm-1, and the low-frequency channel involves a direct relaxation process affected by the phonon-bottleneck effect.

Auxiliary Rather Than Dominant. The Role of Direct Dy-S Coordination in Single-Molecule Magnet Unveiled via ab initio Study

The role of Dy-S coordination in a single-molecule magnet (SMM) is investigated via an ab initio study in a group of mononuclear structures. The SMM performance of this group is well interpreted via a concise criterion consisting of long quantum tunneling of magnetization (QTM) time τQTM and high effective barrier for magnetic reversal Ueff. The best SMMs in the selected group, i.e., 1Dy (CCDC refcode: PUKFAF) and 2Dy (CCDC refcode: NIKSEJ), are just those holding the longest τQTM and the highest Ueff simultaneously. Further analysis based on the crystal field model and ab initio magneto-structural exploration indicates that the influence of Dy-S coordination on the SMM performance of 1Dy is weaker than that of axial Dy-O coordination. Thus, Dy-S coordination is more likely to play an auxiliary role rather than a dominant one. However, if placed at the suitable equatorial position, Dy-S coordination could provide important support for good SMM performance. Consequently, starting from 1Dy, we built two new structures where Dy-S coordination only exists at the equatorial position and two axial positions are occupied by strong Dy-O/Dy-F coordination. Compared to 1Dy and 2Dy, these new ones are predicted to have significantly longer τQTM and higher Ueff, as well as a nearly doubled blocking temperature TB. Thus, they are probable candidates of SMM having clearly improved performance.

Comprehensive Exploration of Optics and Magnetism in a Hydrothermally Synthesized Nickel(II) Complex

The hydrothermal reaction of 2-MeIm (2-MeIm: 2-methylimidazole) with nickel sulfate hexahydrate in methanol afforded a mononuclear complex formulated as [Ni(SO4)(2-MeIm)2(H2O)3]·CH3OH (1). The title compound was described by X-ray single-crystal diffraction, thermal assessment, IR, and UV-vis spectroscopy. The crystal structure of 1 is composed of segregated [Ni(SO4)(2-MeIm)2(H2O)3] neutral entities and a solvent methanol molecule. Two (2-MeIm) ligands, a sulfate group, and a water molecule reside in the equatorial positions of the vertices in this 6-fold coordination. Two aqua ligands lay in the apical positions, resulting in a subtly distorted octahedral framework, as was supported by spectroscopic analysis. The complex's self-assembly is firmly governed by robust O-H···O/N-H···O interactions. Further details on these bonds have been furnished via Hirshfeld surface scrutiny and 2D fingerprint plots. As proven by TGA/DSC analysis, raising the temperature of 1 above 60 °C instigates progressive decomposition stages, which culminates in the production of metal oxide as the ultimate product at 700 °C. The optical analysis suggests the dielectric nature of the material with large direct and indirect gap energies of 5.25 and 4.96 eV, respectively. The results of magnetic studies suggest that 1 undergoes a transition to a magnetically ordered state below 6 K.

Identification, Characterization, and Electronic Structures of Interconvertible Cobalt-Oxygen TAML Intermediates

The reaction of Li[(TAML)CoIII]·3H2O (TAML = tetraamido macrocyclic tetraanionic ligand) with iodosylbenzene at 253 K in acetone in the presence of redox-innocent metal ions (Sc(OTf)3 and Y(OTf)3) or triflic acid affords a blue species 1, which is converted reversibly to a green species 2 upon cooling to 193 K. The electronic structures of 1 and 2 have been determined by combining advanced spectroscopic techniques (X-band electron paramagnetic resonance (EPR), electron nuclear double resonance (ENDOR), X-ray absorption spectroscopy/extended X-ray absorption fine structure (XAS/EXAFS), and magnetic circular dichroism (MCD)) with ab initio theoretical studies. Complex 1 is best represented as an S = 1/2 [(Sol)(TAML•+)CoIII---OH(LA)]- species (LA = Lewis/Brønsted acid and Sol = solvent), where an S = 1 Co(III) center is antiferromagnetically coupled to S = 1/2 TAML•+, which represents a one-electron oxidized TAML ligand. In contrast, complex 2, also with an S = 1/2 ground state, is found to be multiconfigurational with contributions of both the resonance forms [(H-TAML)CoIV═O(LA)]- and [(H-TAML•+)CoIII═O(LA)]-; H-TAML and H-TAML•+ represent the protonated forms of TAML and TAML•+ ligands, respectively. Thus, the interconversion of 1 and 2 is associated with a LA-associated tautomerization event, whereby H+ shifts from the terminal -OH group to TAML•+ with the concomitant formation of a terminal cobalt-oxo species possessing both singlet (SCo = 0) Co(III) and doublet (SCo = 1/2) Co(IV) characters. The reactivities of 1 and 2 at different temperatures have been investigated in oxygen atom transfer (OAT) and hydrogen atom transfer (HAT) reactions to compare the activation enthalpies and entropies of 1 and 2.

Influence of the phonon-bottleneck effect and low-energy vibrational modes on the slow spin-phonon relaxation in Kramers-ions-based Cu(II) and Co(II) complexes with 4-amino-3,5-bis-(pyridin-2-yl)-1,2,4-triazole and dicyanamide

Two new complexes, bis-[4-amino-3,5-bis-(pyridin-2-yl)-1,2,4-triazole-κ2N2,N6]bis-(dicyanamide-κN8)copper(II), [Cu(abpt)2(dca)2] (1) and bis-[4-amino-3,5-bis-(pyridin-2-yl)-1,2,4-triazole-κ2N2,N6]bis-(dicyanamide-κN8)cobalt(II), [Co(abpt)2(dca)2] (2), have been prepared and magneto-structurally characterised. Single crystal studies of both complexes have shown that their crystal structures are molecular, in which the central atoms are six-coordinated in the form of a distorted octahedron by two bidentate abpt and two monodentate dca ligands. Even if both complexes have the same composition and crystallize in the same P1̄ space group, they are not isostructural. Both structures contain strong intermolecular N-H⋯N hydrogen bonds and π-π stacking interactions. IR spectra are consistent with the solved structures of both complexes and confirmed the terminal character of the dca ligands and the bidentate coordination of the abpt ligands. The analysis of the magnetic properties showed that both complexes exhibit field-induced slow spin-phonon relaxation. In both complexes, the slow spin-phonon relaxation is influenced by a severe phonon-bottleneck effect that affects the direct process, a dominant relaxation mechanism at low temperatures in both complexes. The phonon-bottleneck effect in 1 was suppressed by simply reducing the crystallite size, and further analysis of the field dependence of the relaxation time yielded the characteristic energy of vibrational modes of 11 cm-1 participating in the Raman process at low magnetic fields. The analysis of magnetic properties and ab initio calculations confirmed that 2 represents a system with a moderate uniaxial anisotropy yielding an average energy barrier of 82 cm-1 (from all four nonequivalent Co(II) sites in the structure of 2).

Unravelling the Highest Magnetic Anisotropy Among all the n-Shells in 20 Metallocene

Single-molecule magnets (SMMs) with a large magnetization reversal barrier are predominated by the lanthanide systems due to their strong spin-orbit coupling (SOC). However, the transition metals have also emerged as potential contenders and the largest magnetic anisotropy has been identified for a cobalt system among any d-series-based SMMs (Bunting et al. Science 2018, 362, eaat7319). In this work, we have explored the magnetic anisotropy in highly axial ligand field systems of metallocene, having different d-subshell (3d4, 4d4, and 5d4). The wave function-based multireference methods including static and dynamic electron correlations have been employed to investigate the zero-field splitting (ZFS) parameters. Here, we report exceptionally large magnetic anisotropy for a 5d complex of [WCp2]0 with the highest energy barrier that is nearly twice as high as the previous record value for the Co complex. We have also observed that the axial ZFS parameter (D) is increasing down the group in the order of 3d < 4d < 5d, pertaining to a large SOC.

Can Iron-Porphyrins Behave as Single-Molecule Magnets?

The study of magnetic properties, especially the magnetic anisotropy of iron-porphyrin complexes employing multiconfigurational methods, is quite challenging due to many strongly correlated electrons in nearly degenerate orbitals. However, a prerequisite for observing the magnetic anisotropy and slow magnetization relaxation, the zero-field splitting parameter, D, was experimentally observed decades ago for halide-based axially ligated penta-coordinate Fe(III)-porphyrins. In these complexes, the signs of D were reported mostly as positive; in a few cases, inconclusive signs of the D parameter were also mentioned. However, no ab initio calculations have been reported to shed light on this. Deciphering the electronic structure of these penta-coordinated complexes employing the complete active space self-consistent field method and N-electron valence second-order perturbation theory, we confirm the positive D values. However, a negative D value is highly desired to observe the single-molecule magnet properties without an external magnetic field, which we observed in the Fe(II)-porphyrin complexes with axial imidazole ligands instead of halide ligands. The detailed analysis of the multireference wave functions unravels the role of axial ligands in determining the sign and magnitude of the D parameters.

Zero-Field Splitting Tensor of the Triplet Excited States of Aromatic Molecules: A Valence Full-π Complete Active Space Self-Consistent Field Study

A method to predict the D tensor in the molecular frame with multiconfigurational wave functions in large active space was proposed, and the spin properties of the lowest triplets of aromatic molecules were examined with full-π active space; such calculations were challenging because the size of active space grows exponentially with the number of π electrons. In this method, the exponential growth of complexity is resolved by the density matrix renormalization group (DMRG) algorithm. From the D tensor, we can directly determine the direction of the magnetic axes and the ZFS parameters, D- and E-values, of the phenomenological spin Hamiltonian with their signs, which are not usually obtained in ESR experiments. The method using the DMRG-CASSCF wave function can give correct results even when the sign of D- and E-values is sensitive to the accuracy of the prediction of the D tensor and existing methods fail to predict the correct magnetic axes.

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.

Synthesis and photoinduced behavior of DPP-anchored nitronyl nitroxides: a multifaceted approach

Understanding and controlling spin dynamics in organic dyes is of significant scientific and technological interest. The investigation of 2,5-dihydropyrrolo[4,3-c]pyrrolo-1,4-dione derivatives (DPPs), one of the most widely used dyes in many fields, has so far been limited to closed-shell molecules. We present a comprehensive joint experimental and computational study of DPP derivatives covalently linked to two nitronyl nitroxide radicals (DPPTh-NN2). Synthesis, single crystal X-ray diffraction study, photophysical properties, magnetic properties established using steady-state and pulse EPR, fast spin dynamics, and computational modelling using density functional theory and ab initio methods of electronic structure and spectroscopic properties of DPPTh-NN2 are presented. The single-crystal X-ray diffraction analysis of DPPTh-NN2 and computational modeling of its electronic structure suggest that effective conjugation along the backbone leads to noticeable spin-polarization transfer. Calculations using ab initio methods predict a weak exchange interaction of radical centers through a singlet ground state of DPPTh with a small singlet-triplet splitting (ΔEST) of about 25 cm-1 (∼0.07 kcal mol-1). In turn, a strong ferromagnetic exchange interaction between the triplet state of DPPTh chromophore and nitronyl nitroxides (with J ∼ 250 cm-1) was predicted.

Probing the dynamic behaviour and magnetic identification of seven coordinated Mn(II) complexes: a combined AIMD and multi-reference approach

We present an in-depth solution phase dynamics of rare seven coordinated pentagonal bipyramidal Mn(II) complexes, together with their binding affinity anticipated using ab initio molecular dynamics (AIMD) simulations and density functional theory (DFT). Moreover, the simulations at different temperatures (25 °C and 90 °C) interpret the rigidity and stability of the ligands with Mn(II) ions. An intuitive approach for modulating the easy plane magnetic anisotropy of the mononuclear Mn(II) complex has been revealed by this work. In this regard, we have performed an extensive theoretical study based on the ab initio CASSCF/NEVPT2 method, exhibiting the presence of an easy plane magnetic anisotropy with a positive value of axial zero-field splitting (ZFS) parameter D. The complex's magnetic properties and electronic relaxation reveal that the rhombic ZFS term (E) can be modulated as the symmetry around the Mn(II) ion varies. The magnitude of the D-value increased with a more symmetrical equatorial ligand as found in the order of [Mn(pydpa)(H2O)] > [Mn(cbda)(H2O)]- > [Mn(dpaaa)(H2O)]- > [Mn(dpasam)(H2O)]-. Furthermore, we found that substituting the equatorial oxygen atom with heavier S and Se-donor atoms switches the sign of magnetic anisotropy for the Mn(II) complexes. The magnitude of the D-value increased when the energy levels of the ground state (GS) and the first excited state (ES) decreased. The observed magneto-structural correlation reveals that shortening the distance of the axial water molecule (Mn-O(w)) increases the D-value by an order of magnitude for the symmetrical [Mn(pydpa)(H2O)] complex. Overall, the combined analysis of solution phase dynamics of Mn(II) complexes and their magnetic characterization opens up new avenues in coordination chemistry, molecular magnetism, spin-crossover materials, and catalysis.

Zero-Field SMM Behavior Triggered by Magnetic Exchange Interactions and a Collinear Arrangement of Local Anisotropy Axes in a Linear Co3II Complex

A new linear trinuclear Co(II)3 complex with a formula of [{Co(μ-L)}2Co] has been prepared by self-assembly of Co(II) ions and the N3O3-tripodal Schiff base ligand H3L, which is obtained from the condensation of 1,1,1-tris(aminomethyl)ethane and salicylaldehyde. Single X-ray diffraction shows that this compound is centrosymmetric with triple-phenolate bridging groups connecting neighboring Co(II) ions, leading to a paddle-wheel-like structure with a pseudo-C3 axis lying in the Co-Co-Co direction. The Co(II) ions at both ends of the Co(II)3 molecule exhibit distorted trigonal prismatic CoN3O3 geometry, whereas the Co(II) at the middle presents an elongated trigonal antiprismatic CoO6 geometry. The combined analysis of the magnetic data and theoretical calculations reveal strong easy-axis magnetic anisotropy for both types of Co(II) ions (|D| values higher than 115 cm-1) with the local anisotropic axes lying on the pseudo-C3 axis of the molecule. The magnetic exchange interaction between the middle and ends Co(II) ions, extracted by using either a Hamiltonian accounting for the isotropic magnetic coupling and ZFS or the Lines' model, was found to be medium to strong and antiferromagnetic in nature, whereas the interaction between the external Co(II) ions is weak antiferromagnetic. Interestingly, the compound exhibits slow relaxation of magnetization and open hysteresis at zero field and therefore SMM behavior. The significant magnetic exchange coupling found for [{Co(μ-L)}2Co] is mainly responsible for the quenching of QTM, which combined with the easy-axis local anisotropy of the CoII ions and the collinearity of their local anisotropy axes with the pseudo-C3 axis favors the observation of SMM behavior at zero field.

Zero-field splitting parameters within exact two-component theory and modern density functional theory using seminumerical integration

An efficient implementation of zero-field splitting parameters based on the work of Schmitt et al. [J. Chem. Phys. 134, 194113 (2011)] is presented. Seminumerical integration techniques are used for the two-electron spin-dipole contribution and the response equations of the spin-orbit perturbation. The original formulation is further generalized. First, it is extended to meta-generalized gradient approximations and local hybrid functionals. For these functional classes, the response of the paramagnetic current density is considered in the coupled-perturbed Kohn-Sham equations for the spin-orbit perturbation term. Second, the spin-orbit perturbation is formulated within relativistic exact two-component theory and the screened nuclear spin-orbit (SNSO) approximation. The accuracy of the implementation is demonstrated for transition-metal and diatomic main-group compounds. The efficiency is assessed for Mn and Mo complexes. Here, it is found that coarse integration grids for the seminumerical schemes lead to drastic speedups while introducing clearly negligible errors. In addition, the SNSO approximation substantially reduces the computational demands and leads to very similar results as the spin-orbit mean field Ansatz.

State-Interaction Approach for Evaluating g-Tensors within EOM-CC and RAS-CI Frameworks: Theory and Benchmarks

Among various techniques designed for studying open-shell species, electron paramagnetic resonance (EPR) spectroscopy plays an important role. The key quantity measured by EPR is the g-tensor, describing the coupling between an external magnetic field and molecular electronic spin. One theoretical framework for quantum chemistry calculations of g-tensors is based on response theory, which involves substantial developments that are specific to the underlying electronic structure models. A simplified and easier-to-implement approach is based on the state-interaction scheme, in which perturbation is included by considering a small number of states. We describe and benchmark the state-interaction approach using equation-of-motion coupled-cluster and restricted-active-space configuration interaction wave functions. The analysis confirms that this approach can deliver accurate results and highlights caveats of applying it, such as a choice of the reference state, convergence with respect to the number of states used in calculations, etc. The analysis also contributes toward a better understanding of challenges in calculations of higher-order properties using approximate wave functions.

SOiCISCF: Combining SOiCI and iCISCF for Variational Treatment of Spin-Orbit Coupling

It has recently been shown that the SOiCI approach [Zhang, N.; J. Phys.: Condens. Matter 2022, 34, 224007], in conjunction with the spin-separated exact two-component relativistic Hamiltonian, can provide very accurate fine structures of systems containing heavy elements by treating electron correlation and spin-orbit coupling (SOC) on an equal footing. Nonetheless, orbital relaxations/polarizations induced by SOC are not yet fully accounted for due to the use of scalar relativistic orbitals. This issue can be resolved by further optimizing the still real-valued orbitals self-consistently in the presence of SOC, as done in the spin-orbit coupled CASSCF approach [Ganyushin, D.; et al. J. Chem. Phys. 2013, 138, 104113] but with the iCISCF algorithm [Guo, Y.; J. Chem. Theory Comput. 2021, 17, 7545-7561] for large active spaces. The resulting SOiCISCF employs both double group and time reversal symmetries for computational efficiency and the assignment of target states. The fine structures of p-block elements are taken as showcases to reveal the efficacy of SOiCISCF.

Unravelling the effect of paramagnetic Ni2+ on the 13C NMR shift tensor for carbonate in Mg2-xNixAl layered double hydroxides by quantum-chemical computations

Structural disorder and low crystallinity render it challenging to characterise the atomic-level structure of layered double hydroxides (LDH). We report a novel multi-step, first-principles computational workflow for the analysis of paramagnetic solid-state NMR of complex inorganic systems such as LDH, which are commonly used as catalysts and energy storage materials. A series of 13CO32--labelled Mg2-xNixAl-LDH, x ranging from 0 (Mg2Al-LDH) to 2 (Ni2Al-LDH), features three distinct eigenvalues δ11, δ22 and δ33 of the experimental 13C chemical shift tensor. The δii correlate directly with the concentration of the paramagnetic Ni2+ and span a range of |δ11 - δ33| ≈ 90 ppm at x = 0, increasing to 950 ppm at x = 2. In contrast, the isotropic shift, δiso(13C), only varies by -14 ppm in the series. Detailed insight is obtained by computing (1) the orbital shielding by periodic density-functional theory involving interlayer water, (2) the long-range pseudocontact contribution of the randomly distributed Ni2+ ions in the cation layers (characterised by an ab initio susceptibility tensor) by a lattice sum, and (3) the close-range hyperfine terms using a full first-principles shielding machinery. A pseudohydrogen-terminated two-layer cluster model is used to compute (3), particularly the contact terms. Due to negative spin density contribution at the 13C site arising from the close-by Ni2+ sites, this step is necessary to reach a semiquantitative agreement with experiment. These findings influence future NMR investigations of the formally closed-shell interlayer species within LDH, such as the anions or water. Furthermore, the workflow is applicable to a variety of complex materials.

Dynamic Magnetic Properties of Germole-ligated Lanthanide Sandwich Complexes

The first germole-ligated single-molecule magnets are reported, with contrasting properties found for the near-linear sandwich complexes [(η8 -COT)Ln(η5 -CpGe ]- , where Ln=Dy (1Dy ) or Er (1Er ), COT is cyclo-octatetraenyl and CpGe is [GeC4 -2,5-(SiMe3 )2 -3,4-Me2 ]2- . Whereas 1Er has an energy barrier of 120(1) cm-1 in zero applied field and open hysteresis loops up to 10 K, the relaxation in 1Dy is characterized by quantum tunneling within the ground state.

Synthesis and isolation of a triplet bismuthinidene with a quenched magnetic response

Large Spin-Orbit Coupling (SOC) is an intrinsic property of the heavy-elements that directly affects the electronic structures of the compounds. Herein we report the synthesis and characterization of a mono-coordinate bismuthinidene featuring a rigid and bulky ligand. All magnetic measurements (SQUID, NMR) point to a diamagnetic compound. However, multiconfigurational quantum chemical calculations predict the ground state of the compound to be dominated (76%) by a spin-triplet. The apparent diamagnetism is explained by an extremely large SOC induced positive zero-field-splitting of more than 4500 cm^-1 that leaves the M_S = 0 magnetic sublevel thermally isolated in the electronic ground state.

On the Single-Molecule Magnetic Behavior of Linear Iron(I) Arylsilylamides

The rational design of 3d-metal-based single-molecule magnets (SMM) requires a fundamental understanding of their intrinsic electronic and structural properties and how they translate into experimentally observable features. Here, we determined the magnetic properties of the linear iron(I) silylamides K{crypt}[FeL₂] and [KFeL₂] (L = -N(Dipp)SiMe₃; crypt = 4,7,13,16,21,24-Hexaoxa-1,10-diazabicyclo[8.8.8]hexacosan). For the former, slow-relaxation of the magnetization with a spin reversal barrier of U_eff = 152 cm^-1 as well as a closed-waist magnetic hysteresis and magnetic blocking below 2.5 K are observed. For the more linear [KFeL₂], in which the potassium cation is encapsulated by the aryl substituents of the amide ligands, the relaxation barrier and the blocking temperature increase to U_eff = 184 cm^-1 and T_B = 4.5 K, respectively. The increase is rationalized by a more pronounced axial anisotropy in [KFeL₂] determined by dc-SQUID magnetometry. The effective relaxation barrier of [KFeL₂] is in agreement with the energy spacing between the ground and first-excited magnetic states, as obtained by field-dependent IR-spectroscopy (178 cm^-1), magnetic measurements (208 cm^-1), as well as theoretical analysis (212 cm^-1). In comparison with the literature, the results show that magnetic coercivity in linear iron(I) silylamides is driven by the degree of linearity in conjunction with steric encumbrance, whereas the ligand symmetry is a marginal factor.

Computational NMR of the iron pyrazolylborate complexes [Tp2Fe]+ and Tp2Fe including solvation and spin-crossover effects

Transition metal complexes have important roles in many biological processes as well as applications in fields such as pharmacy, chemistry and materials science. Paramagnetic nuclear magnetic resonance (pNMR) is a valuable tool in understanding such molecules, and theoretical computations are often advantageous or even necessary in the assignment of experimental pNMR signals. We have employed density functional theory (DFT) and the domain-based local pair natural orbital coupled-cluster method with single and double excitations (DLPNO-CCSD), as well as a number of model improvements, to determine the critical hyperfine part of the chemical shifts of the iron pyrazolylborate complexes [Tp₂Fe]⁺ and Tp₂Fe using a modern version of the Kurland-McGarvey theory, which is based on parameterising the hyperfine, electronic Zeeman and zero-field splitting interactions via the parameters of the electron paramagnetic resonance Hamiltonian. In the doublet [Tp₂Fe]⁺ system, the calculations suggest a re-assignment of the ¹³C signal shifts. Consideration of solvent via the conductor-like polarisable continuum model (C-PCM) versus explicit solvent molecules reveals C-PCM alone to be insufficient in capturing the most important solvation effects. Tp₂Fe exhibits a spin-crossover effect between a high-spin quintet (S = 2) and a low-spin singlet (S = 0) state, and its recorded temperature dependence can only be reproduced theoretically by accounting for the thermal Boltzmann distribution of the open-shell excited state and the closed-shell ground-state occupations. In these two cases, DLPNO-CCSD is found, in calculating the hyperfine couplings, to be a viable alternative to DFT, the demonstrated shortcomings of which have been a significant issue in the development of computational pNMR.

Simulating Spin-Orbit Coupling with Quasidegenerate N-Electron Valence Perturbation Theory

We present the first implementation of spin-orbit coupling effects in fully internally contracted second-order quasidegenerate N-electron valence perturbation theory (SO-QDNEVPT2). The SO-QDNEVPT2 approach enables the computations of ground- and excited-state energies and oscillator strengths combining the description of static electron correlation with an efficient treatment of dynamic correlation and spin-orbit coupling. In addition to SO-QDNEVPT2 with the full description of one- and two-body spin-orbit interactions at the level of two-component Breit-Pauli Hamiltonian, our implementation also features a simplified approach that takes advantage of spin-orbit mean-field approximation (SOMF-QDNEVPT2). The accuracy of these methods is tested for the group 14 and 16 hydrides, 3d and 4d transition metal ions, and two actinide dioxides (neptunyl and plutonyl dications). The zero-field splittings of group 14 and 16 molecules computed using SO-QDNEVPT2 and SOMF-QDNEVPT2 are in good agreement with the available experimental data. For the 3d transition metal ions, the SO-QDNEVPT2 method is significantly more accurate than SOMF-QDNEVPT2, while no substantial difference in the performance of two methods is observed for the 4d ions. Finally, we demonstrate that for the actinide dioxides the results of SO-QDNEVPT2 and SOMF-QDNEVPT2 are in good agreement with the data from previous theoretical studies of these systems. Overall, our results demonstrate that SO-QDNEVPT2 and SOMF-QDNEVPT2 are promising multireference methods for treating spin-orbit coupling with a relatively low computational cost.

Change in the Electronic Structure of the Cobalt(II) Ion in a One-Dimensional Polymer with Flexible Linkers Induced by a Structural Phase Transition

A new 1D-coordination polymer [Co(Piv)₂(NH₂(CH₂)₆NH₂)]_n (1, Piv is Me₃CCO₂^- anion) was obtained, the mononuclear fragments {Co(O₂CR)₂} within which are linked by μ-bridged molecules of hexamethylenediamine (NH₂(CH₂)₆NH₂). For this compound, two different monoclinic C2/c (α-1) and P2/n (β-1) phases were found at room temperature by single-crystal X-ray diffraction analysis, with a similar structure of chains and their packages in unit cells. The low-temperature phase (γ-1) of crystal 1 at 150 K corresponds to the triclinic space group P-1. As the temperature decreases, the structural phase transition (SPT) in the α-1 and β-1 crystals is accompanied by an increase in the crystal packing density caused by the rearrangements of both H-bonds and the nearest ligand environment of the cobalt atom ("octahedral CoN₂O₄ around the metal center at room temperature" → "pseudo-tetrahedral CoN₂O₂ at 150 K"). The SPT was confirmed by DSC in the temperature range 210-150 K; when heated above 220 K, anomalies in the behavior of the heat flow are observed, which may be associated with the reversibility of SPT; endo effects are observed up to 300 K. The SPT starts below 200 K. At 100 K, a mixture of phases was found in sample 1: 27% α-1 phase, 61% γ-1 phase. In addition, at 100 K, 12% of the new δ-1 phase was detected, which was identified from the diffraction pattern at 260 K upon subsequent heating: the a,b,c-parameters and unit cell volume are close to the structure parameters of γ-1, and the values of the α,β,γ-angles are significantly different. Further heating leads to a phase transition from δ-1 to α-1, which both coexist at room temperature. According to the DC magnetometry data, during cooling and heating, the χ_MT(T) curves for 1 form a hysteresis loop with ~110 K, in which the difference in the χ_MT values reaches 9%. Ab initio calculations of the electronic structure of cobalt(II) in α-1 and γ-1 have been performed. Based on the EPR data at 10 K and the ab initio calculations, the behavior of the χ_MT(T) curve for 1 was simulated in the temperature range of 2-150 K. It was found that 1 exhibits slow magnetic relaxation in a field of 1000 Oe.

Computational design of magnetic molecules and their environment using quantum chemistry, machine learning and multiscale simulations

Having served as a playground for fundamental studies on the physics of d and f electrons for almost a century, magnetic molecules are now becoming increasingly important for technological applications, such as magnetic resonance, data storage, spintronics and quantum information. All of these applications require the preservation and control of spins in time, an ability hampered by the interaction with the environment, namely with other spins, conduction electrons, molecular vibrations and electromagnetic fields. Thus, the design of a novel magnetic molecule with tailored properties is a formidable task, which does not only concern its electronic structures but also calls for a deep understanding of the interaction among all the degrees of freedom at play. This Review describes how state-of-the-art ab initio computational methods, combined with data-driven approaches to materials modelling, can be integrated into a fully multiscale strategy capable of defining design rules for magnetic molecules.

Field-Induced Slow Magnetic Relaxation in CoII Cyclopropane-1,1-dicarboxylates

New Co^II substituted malonate field-induced molecular magnets {[Rb₆Co₃(cpdc)₆(H₂O)₁₂]∙6H₂O}_n (1) and [Cs₂Co(cpdc)₂(H₂O)₆]_n (2) (where cpdc²⁻ stands for cyclopropane-1,1-dicarboxylic acid dianions) were synthesized. Both compounds contain mononuclear bischelate fragments {Co^II(cpdc)₂(H₂O)₂}²⁻ where the quasi-octahedral cobalt environment (CoO₆) is complemented by water molecules in apical positions. The alkali metal atoms play the role of connectors between the bischelate fragments to form 3D and 2D polymeric structures for 1 and 2, respectively. Analysis of dc magnetic data using the parametric Griffith Hamiltonian for high-spin Co^II supported by ab initio calculations revealed that both compounds have an easy axis of magnetic anisotropy. Compounds 1 and 2 exhibit slow magnetic relaxation under an external magnetic field (H_DC = 1000 and 1500 Oe, respectively).

Magnetic exchange and valence delocalization in a mixed valence [Fe2+Fe3+Te2]+ complex: insights from theory and interpretations of magnetic and spectroscopic data

A mixed valence binuclear Fe^2.5+-Fe^2.5+ (Robin-Day Class III) transition metal complex, [Fe^2.5+μTe₂Fe^2.5+]^1-, composed of two FeN₂Te₂ pseudo-tetrahedral units with μ-bridging Te^2- ligands was reported to exist in an unprecedented S = 3/2 ground state (Nature Chemistry, https://doi.org/10.1038/s41557-021-00853-5). For this and the homologous complexes containing Se^2- and S^2-, the Anderson-Hasegawa double exchange spin-Hamiltonian was broadly used to interpret the corresponding structural, spectroscopic and magnetic data. First principles multireference ab initio calculations are used here to simulate magnetic and spectroscopic EPR data; analysis of the results affords a rationale for the stabilization of the S = 3/2 ground state of the Fe₂ pair. Complete Active Space Self-Consistent Field (CASSCF) calculations and dynamical correlation accounted for by means of N-Electron Valence Perturbation Theory to Second Order (NEVPT2) reproduce well the g-factors determined from simulations of X-band EPR spectra. A crucial technical tool to achieve these results is: (i) use of a localized orbital formulation of the many-particle problem at the scalar-relativistic CASSCF step; (ii) choice of state averaging over states of a given spin (at the CASCI/NEVPT2 step); and (iii) accounting for spin-orbit coupling within the non-relativistic Born-Oppenheimer (BO) many-particle basis using Quasi-Degenerate Perturbation Theory (QDPT). The inclusion of the S = 5/2 spin manifold reproduced the observed increase in the magnetic susceptibility (χT) in the high temperature range (T > 100 K), which is explained by thermal population of the S = 5/2 excited state at energy 160 cm^-1 above the S = 3/2 ground state. Theoretical values of χT from experimentally reported data points in the temperature range from 3 to 30 K were further computed and analyzed using a model which takes spin-phonon coupling into account. The model considerations and the computational protocols of this study are generally applicable to any Class I/II mixed valence dimer. The work can potentially stimulate further experimental and theoretical work on bi- and oligonuclear transition metal complexes of importance to bioinorganic chemistry and life sciences.

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.

Slow magnetic relaxation in dinuclear Co(III)-Co(II) complexes containing a five-coordinated Co(II) centre with easy-axis anisotropy

Two air-stable Co(III)-Co(II) mixed-valence complexes of molecular formulas [Co^IICo^III(L)(DMAP)₃(CH₃COO)]·H₂O·CH₃OH (1) and [Co^IICo^III(L)(4-Pyrrol)₃ (CH₃COO)]·0.5CH₂Cl₂ (2) (H₄L = 1,3-bis-(5-methyl pyrazole-3-carboxamide) propane; DMAP = 4-dimethylaminopyridine; and 4-Pyrrol = 4-pyrrolidinopyridine) were synthesized and characterized by single-crystal X-ray crystallography, high-field electron paramagnetic resonance (HFEPR) spectroscopy, and magnetic measurements. Both complexes possess one five-coordinated paramagnetic Co(II) ion and one six-coordinated Co(III) ion with octahedral geometry. Direct-current magnetic susceptibility and magnetization measurements show the easy-axis magnetic anisotropy that is also confirmed by low-temperature HFEPR measurements and theoretical calculations. Frequency- and temperature-dependent alternating-current magnetic susceptibility measurements reveal their field-assisted slow magnetic relaxation, which is a characteristic behavior of single-molecule magnets (SMMs), caused by the individual Co(II) ion. The effective energy barrier of complex 1 (49.2 cm^-1) is significantly higher than those of the other dinuclear Co(III)-Co(II) SMMs. This work hence presents the first instance of the dinuclear Co(III)-Co(II) single-molecule magnets with a five-coordinated environment around the Co(II) ion.

A trust-region augmented Hessian implementation for state-specific and state-averaged CASSCF wave functions

In this work, we present a one-step second-order converger for state-specific (SS) and state-averaged (SA) complete active space self-consistent field (CASSCF) wave functions. Robust convergence is achieved through step restrictions using a trust-region augmented Hessian (TRAH) algorithm. To avoid numerical instabilities, an exponential parameterization of variational configuration parameters is employed, which works with a nonredundant orthogonal complement basis. This is a common approach for SS-CASSCF and is extended to SA-CASSCF wave functions in this work. Our implementation is integral direct and based on intermediates that are formulated in either the sparse atomic-orbital or small active molecular-orbital basis. Thus, it benefits from a combination with efficient integral decomposition techniques, such as the resolution-of-the-identity or the chain-of-spheres for exchange approximations. This facilitates calculations on large molecules, such as a Ni(II) complex with 231 atoms and 5154 basis functions. The runtime performance of TRAH-CASSCF is competitive with the other state-of-the-art implementations of approximate and full second-order algorithms. In comparison with a sophisticated first-order converger, TRAH-CASSCF calculations usually take more iterations to reach convergence and, thus, have longer runtimes. However, TRAH-CASSCF calculations still converge reliably to a true minimum even if the first-order algorithm fails.

Harnessing Colossal Magnetic Anisotropy in Sandwiched 3d2-Metallocenes

Single-molecule magnets are gaining attention in recent years with the growing focus on achieving higher barriers of magnetization reversal. Metallocenes, owing to their unique sandwiched structure, assure themselves as plausible molecular systems for the development of novel single-molecule magnets (SMMs). Here in this work, we have explicitly investigated metallocenes of first-row transition elements, along with their one-electron-oxidized (cationic) and -reduced (anionic) analogues, for their magnetic anisotropies by adopting multireference ab initio calculations. Herein, we report a high magnetic anisotropy for 3d² systems among all 3d-metallocenes.

Giant Rotational Magnetocaloric Effect in Ni(en)(H2O)4·2H2O: Experiment and Theory

An experimental study of the rotational magnetocaloric effect in Ni(en)(H2O)4SO4∙2H2O (en = ethylenediamine) single crystal is presented. The study was carried out at temperatures above 2 K and was associated with adiabatic crystal rotation between the easy plane and hard axis in magnetic fields up to 7 T. The magnetocaloric properties of the studied system were investigated by isothermal magnetization measurement. The experimental observations were completed with ab initio calculations of the anisotropy parameters. A large rotational magnetic entropy change ≈12 Jkg−1K−1 and ≈16.9 Jkg−1K−1 was achieved in 5 T and 7 T, respectively. The present study suggests a possible application of this material in low-temperature refrigeration since the adiabatic rotation of the single crystal in 7 T led to a cooldown of the sample from the initial temperature of 4.2 K down to 0.34 K. Finally, theoretical calculations show that S = 1 Ni(II)-based systems with easy-plane anisotropy can have better rotational magnetocaloric properties than costly materials containing rare-earth elements in their chemical structures.

SOiCI and iCISO: combining iterative configuration interaction with spin-orbit coupling in two ways

The near-exact iCIPT2 approach for strongly correlated systems of electrons, which stems from the combination of iterative configuration interaction (iCI, an exact solver of full CI) with configuration selection for static correlation and second-order perturbation theory (PT2) for dynamic correlation, is extended to the relativistic domain. In the spirit of spin separation, relativistic effects are treated in two steps: scalar relativity is treated by the infinite-order, spin-free part of the exact two-component (X2C) relativistic Hamiltonian, whereas spin-orbit coupling (SOC) is treated by the first-order, Douglas-Kroll-Hess-like SOC operator derived from the same X2C Hamiltonian. Two possible combinations of iCIPT2 with SOC are considered, i.e., SOiCI and iCISO. The former treats SOC and electron correlation on an equal footing, whereas the latter treats SOC in the spirit of state interaction, by constructing and diagonalizing an effective spin-orbit Hamiltonian matrix in a small number of correlated scalar states. Both double group and time reversal symmetries are incorporated to simplify the computation. Pilot applications reveal that SOiCI is very accurate for the spin-orbit splitting (SOS) of heavy atoms, whereas the computationally very cheap iCISO can safely be applied to the SOS of light atoms and even of systems containing heavy atoms when SOC is largely quenched by ligand fields.

Dielectric-Optical Switches: Photoluminescent, EPR, and Magnetic Studies on Organic-Inorganic Hybrid (azetidinium)2MnBr4

A new organic-inorganic hybrid, AZEMnBr, has been synthesized and characterized. The thermal differential scanning calorimetry, differential thermal analysis, and thermogravimetric analyses indicate one structural phase transition (PT) at 346 and 349 K, on cooling and heating, respectively. AZEMnBr crystallizes at 365 K in the orthorhombic, Pnma, structure, which transforms to monoclinic P2₁/n at 200 K. Due to the X-ray diffraction studies, the anionic MnBr₄^2- moiety is discrete. The azetidinium cations show dynamical disorder in the high-temperature phase. In the proposed structural PT, the mechanism is classified as an order-disorder type. The structural changes affect the dielectric response. In this paper, the multiple switches between low- and high- dielectric states are presented. In addition, it was also observed that the crystal possesses a mutation of fluorescent properties between phase ON and OFF in the PT's point vicinity. We also demonstrate that EPR spectroscopy effectively detects PTs in structurally diverse Mn(II) complexes. AZEMnBr compounds show DC magnetic data consistent with the S = 5/2 spin system with small zero-field splitting, which was confirmed by EPR measurements and slow magnetic relaxation under the moderate DC magnetic field typical for a single-ion magnet behavior. Given the above, this organic-inorganic hybrid can be considered a rare example of multifunctional materials that exhibit dielectric, optical, and magnetic activity.

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.

Hexacoordinate high-spin Fe(III) complexes composed of pentadentate amino-type ligand and pseudohalido coligands

A new series of Fe(III) mononuclear complexes of [Fe(Lam)(X)] type {where X = Cl (1), NCSe (2), NCS (3), N3 (4), NCO (5)} have been prepared and characterized in detail....

Modulation of the magnetic dynamics of pentagonal-bipyramidal Co(ii) complexes by fine-tuning the coordination microenvironment

Dynamic magnetic behaviours of a series of Co(ii) SIMs with pentagonal-bipyramidal geometry have been modulated by an alteration of the ligand field effect.

Magnetically bistable cobalt-dioxolene complexes with a tetradentate N-donor base

We describe a family of cobalt-dioxolene complexes exhibiting a high diversity of magnetic properties: from field-induced single-ion magnetism to thermally induced valence-tautomerism.

First-Principles Study of Optical Absorption Energies, Ligand Field and Spin-Hamiltonian Parameters of Cr3+ Ions in Emeralds

Herein, we study the electronic structure, energies, and vibronic structure of optical d-d transitions of Cr³⁺ ions doped in beryl (Be₃Si₆Al₂O₁₈:Cr³⁺, emerald). A computational protocol is developed that combines periodic density functional theory (for modeling of the bulk crystalline lattice of emerald) and the multireference configuration interaction complete active space self-consistent field method supplemented with n-electron valence second-order perturbation theory (for the calculation of the energy levels, wave functions, and spin-Hamiltonian and ligand-field parameters of the trigonal Cr³⁺ centers in the [CrO₆]^9- clusters embedded in an extended point charge field). Ligand-field parameters were extracted from mapping the effective ligand-field Hamiltonian onto the full many-particle Hamiltonian from one side and from a direct fit to energies of computed d-d transitions on the other side. These have been analyzed using ab initio ligand-field theory. The quality of the theoretical predictions is critically assessed through a detailed comparison with the available experimental data.

Rigid versions of PDTA4- incorporating a 1,3-diaminocyclobutyl spacer for Mn2+ complexation: stability, water exchange dynamics and relaxivity

Rigid derivatives of the acyclic ligand PDTA4- (H₄PDTA = propylenediamine-N,N,N',N'-tetraacetic acid) were prepared by functionalization of a 1,3-diaminocyclobutyl spacer. The new ligands contain either four acetate groups attached to the central scaffold (H4L1) or incorporate pyridyl (H2L2) or propylamide (H2L3) units replacing two of the carboxylate groups. The ligand protonation constants and the stability constants of their Mn²⁺ complexes were determined using potentiometric and spectrophotometric titrations. The stability of the [Mn(L1)]^2- complex was found to be significantly higher than that of the flexible [Mn(PDTA)]^2- derivative (log K_MnL = 10.78 and 10.01, respectively). A detailed study of the ¹H Nuclear Magnetic Relaxation Dispersion (NMRD) profiles and ¹⁷O NMR measurements evidence that the [Mn(L1)]^2- and [Mn(L2)] complexes display a hydration equilibrium in solution involving a seven-coordinate species with an inner-sphere water molecule and a six-coordinate species that lacks a coordinated water molecule. As a result the ¹H relaxivities of these complexes are somewhat lower than that of [Mn(EDTA)]^2- and related systems. The introduction of propylamide groups in [Mn(L3)] shifts the hydration equilibrium to the seven-coordinate species, which results in a ¹H relaxivity (r_1p = 3.7 mM^-1 s^-1 at 22 MHz and 25 °C) exceeding that of [Mn(EDTA)]^2- (r_1p = 3.3 mM^-1 s^-1 at 22 MHz and 25 °C). The parameters that control the relaxivities in this family of complexes were determined by simultaneous fitting of the experimental ¹H NMRD and ¹⁷O NMR data (transverse relaxation rates and chemical shifts), with the aid of computational studies performed at the DFT and CASSCF/NEVPT2 levels. These studies provide detailed insight of the parameters that control the efficiency of these relaxation agents at the molecular level.

Modern physical methods for the molecular design of single-molecule magnets

Many paramagnetic metal complexes have emerged as unique magnetic materials (single-molecule magnets), which behave as conventional magnets at the single-molecule level, thereby making it possible to use them in modern devices for data storage and processing. The rational design of these complexes, however, requires a deep understanding of the physical laws behind a single-molecule magnet behaviour, the mechanisms of magnetic relaxation that determines the magnetic properties and the relationship of these properties with the structure of single-molecule magnets. This review focuses on the physical methods providing such understanding, including different versions and various combinations of magnetometry, electron paramagnetic and nuclear magnetic resonance spectroscopy, optical spectroscopy and X-ray diffraction. Many of these methods are traditionally used to determine the composition and structure of new chemical compounds. However, they are rarely applied to study molecular magnetism.

The bibliography includes 224 references.

Understanding the Effect of the Electron Spin Relaxation on the Relaxivities of Mn(II) Complexes with Triazacyclononane Derivatives

Investigating the relaxation of water ¹H nuclei induced by paramagnetic Mn(II) complexes is important to understand the mechanisms that control the efficiency of contrast agents used in diagnostic magnetic resonance imaging (MRI). Herein, a series of potentially hexadentate triazacyclononane (TACN) derivatives containing different pendant arms were designed to explore the relaxation of the electron spin in the corresponding Mn(II) complexes by using a combination of ¹H NMR relaxometry and theoretical calculations. These ligands include 1,4,7-triazacyclononane-1,4,7-triacetic acid (H₃NOTA) and three derivatives in which an acetate group is replaced by sulfonamide (H₃NO2ASAm), amide (H₂NO2AM), or pyridyl (H₂NO2APy) pendants. The analogue of H₃NOTA containing three propionate pendant arms (H₃NOTPrA) was also investigated. The X-ray structure of the derivative containing two acetate groups and a sulfonamide pendant arm [Mn(NO2ASAm)]^- evidenced six-coordination of the ligand to the metal ion, with the coordination polyhedron being close to a trigonal prism. The relaxivities of all complexes at 20 MHz and 25 °C (1.1-1.3 mM^-1 s^-1) are typical of systems that lack water molecules coordinated to the metal ion. The nuclear magnetic relaxation profiles evidence significant differences in the relaxivities of the complexes at low fields (<1 MHz), which are associated with different spin relaxation rates. The zero field splitting (ZFS) parameters calculated by using DFT and CASSCF methods show that electronic relaxation is relatively insensitive to the nature of the donor atoms. However, the twist angle of the two tripodal faces that delineate the coordination polyhedron, defined by the N atoms of the TACN unit (lower face) and the donor atoms of the pendant arms (upper face), has an important effect in the ZFS parameters. A twist angle close to the ideal value for an octahedral coordination (60°), such as that in [Mn(NOTPrA)]^-, leads to a small ZFS energy, whereas this value increases as the coordination polyhedron approaches to a trigonal prism.

Field-induced single-ion magnet based on a quasi-octahedral Co(II) complex with mixed sulfur-oxygen coordination environment

Synthesis and characterization of structure and magnetic properties of the quasi-octahedral complex (pipH₂)[Co(TDA)₂] 2H₂O (I), (pipH₂²⁺ = piperazine dication, TDA^2- = thiodiacetic anion) are described. X-ray diffraction studies reveal the first coordination sphere of the Co(II) ion, consisting of two chelating tridentate TDA ligands with a mixed sulfur-oxygen strongly elongated octahedral coordination environment. SQUID magnetometry, frequency-domain Fourier-transform (FD-FT) THz-EPR spectroscopy, and high-level ab initio SA-CASSCF/NEVPT2 quantum chemical calculations reveal a strong "easy-plane" type magnetic anisotropy (D ≈ +54 cm^-1) of complex I. The complex shows field-induced slow relaxation of magnetization at an applied DC field of 1000 Oe.

Spin-free formulation of the multireference driven similarity renormalization group: A benchmark study of first-row diatomic molecules and spin-crossover energetics

We report a spin-free formulation of the multireference (MR) driven similarity renormalization group (DSRG) based on the ensemble normal ordering of Mukherjee and Kutzelnigg [J. Chem. Phys. 107, 432 (1997)]. This ensemble averages over all microstates of a given total spin quantum number, and therefore, it is invariant with respect to SU(2) transformations. As such, all equations may be reformulated in terms of spin-free quantities and they closely resemble those of spin-adapted closed-shell coupled cluster (CC) theory. The current implementation is used to assess the accuracy of various truncated MR-DSRG methods (perturbation theory up to third order and iterative methods with single and double excitations) in computing the constants of 33 first-row diatomic molecules. The accuracy trends for these first-row diatomics are consistent with our previous benchmark on a small subset of closed-shell diatomic molecules. We then present the first MR-DSRG application on transition-metal complexes by computing the spin splittings of the [Fe(H₂O)₆]²⁺ and [Fe(NH₃)₆]²⁺ molecules. A focal point analysis (FPA) shows that third-order perturbative corrections are essential to achieve reasonably converged energetics. The FPA based on the linearized MR-DSRG theory with one- and two-body operators and up to a quintuple-ζ basis set predicts the spin splittings of [Fe(H₂O)₆]²⁺ and [Fe(NH₃)₆]²⁺ to be -35.7 and -17.1 kcal mol^-1, respectively, showing good agreement with the results of local CC theory with singles, doubles, and perturbative triples.