Accurate state energetics in spin-crossover systems using pure density functional theory.

The energy difference between different spin states of systems with transition metals is an outstanding challenge for electronic structure calculation methods. The small energy difference between high- and low-spin states in spin-crossover systems makes most post-Hartree-Fock or density functional theory-based methods provide inaccurate values. A test case of twenty systems showing spin transitions has been used to evaluate the accuracy of a new family of training meta-GGA (Generalized Gradient Approximation) functionals. One of the functionals of this new family provides comparable or even better values than the best functional reported so far for this type of system, the TPSSh hybrid meta-GGA functional, but without having to use the exact exchange term. It also improves the results obtained with the r2SCAN meta-GGA functional, which was the best alternative to the TPSSh hybrid functional. This makes it possible to calculate the spin energetics of any kind of compound, especially large systems or periodic structures where the exact exchange requires large computational resources.

Benchmarking Periodic Density Functional Theory Calculations for Spin-State Energies in Spin-Crossover Systems.

Spin energetics is one of the biggest challenges associated with energy calculations for electronic structure methods. The energy differences of the spin states in spin-crossover compounds are very small, making them one of the most difficult systems to calculate. Few methods provide accurate results for calculating these energy differences. In addition, studies have usually focused on calculating energetics of single molecules, while spin-crossover properties are usually experimentally studied in the solid phase. In this paper, we have used periodic boundary conditions employing methods based on density functional theory to calculate the high- and low-spin energy differences for a test case of 20 extended systems. Compounds with different metals and ligands have been selected, and the results indicate that a semiquantitative description of the energy differences can be obtained with the combination of geometry optimization using the PBE functional including many-body dispersion approach and the use of meta-GGA functionals, such as r2SCAN but especially KTBM24, for the energy calculation. Other hybrid functionals, such as TPSSh, give generally good results, but the calculation of the exact exchange with periodic boundary conditions involves a huge increase in computer time and computational resources. It makes the proposed nonhybrid functional approach (KTBM24//PBE+MB) a great advantage for the study of periodic systems.

Thermal Jahn-Teller Distortion Changes and Slow Relaxation of Magnetization in Mn(III) Schiff Base Complexes.

The impact that the anion and alkyl group has on the electronic structures and magnetic properties of four mononuclear Mn(III) complexes is explored in [Mn(salEen-Br)2]Y (salEen-Br = 2-{[2-(ethylamino)ethylimino]methyl}-4-Br-phenol; Y = ClO4- 1 and BF4-·1/3CH2Cl2 2) and [Mn(salBzen-Br)2]Y (salBzen-Br = 2-{[2-(benzylamino)ethylimino]methyl}-4-Br-phenol; Y = ClO4- 3 and BF4- 4). X-ray structures of [Mn(salEen-Br)2]ClO4·0.45C6H14 1-hexane, [Mn(salEen-Br)2]BF4·0.33CH2Cl2·0.15C6H14 2-dcm-hexane, and 3-4 reveal that they crystallize in ambient conditions in the monoclinic P21/c space group. Lowering the temperature, 2-dcm-hexane uniquely exhibits a structural phase transition toward a monoclinic P21/n crystal structure determined at 100 K with the unit cell trebling in size. Remarkably, at room temperature, the axially elongated Jahn-Teller axis in 2-dcm-hexane is poorly defined but becomes clearer at low temperature after the phase transition. Magnetic susceptibility measurements of 1-4 reveal that only 3 and 4 show slow relaxation of magnetization with Δeff/kB = 27.9 and 20.7 K, implying that the benzyl group is important for observing single-molecule magnet (SMM) properties. Theoretical calculations demonstrate that the alkyl group subtly influences the orbital levels and therefore very likely the observed SMM properties.

New series of mononuclear ß-diketonate cerium(III) field induced single-molecule magnets.

Five new ß-diketonate Ce3+ mononuclear complexes, [Ce(Btfa)3(H2O)2] (1), [Ce(Btfa)3(phen)] (2), [Ce(Btfa)3(bipy)] (3), [Ce(Btfa)3(terpy)] (4) and [Ce(Btfa)3(bathophen)(DMF)] (5), where Btfa- = 4,4,4-trifluoro-1-phenyl-1,3-butanedionate, phen = 1,10-phenanthroline, bipy = 2,2'-bipyridyl, terpy = 2,2':6',2''-terpyridine and bathophen = 4,7-diphenyl-1,10-phenanthroline, have been synthesized and structurally characterized through X-ray diffraction of single crystals. The central Ce3+ atom displays a coordination number of 8 for 1, 2 and 3 and of 9 for 4 and 5. Under a 0 T external magnetic field, none of the given compounds exhibits single molecule magnet (SMM) behaviour. However, a small magnetic field, between 0.02 and 0.1 T, is enough for all the compounds to exhibit slow relaxation of the magnetization. A comprehensive magnetic analysis, with experimental magnetic data and ab initio calculations, was undertaken for all the complexes, and the study highlights the significance of the different spin relaxation mechanisms that must be considered for a Ce3+ lanthanide ion.

Increasing the Magnetic Blocking Temperature of Single-Molecule Magnets.

The synthesis of single-molecule magnets (SMMs), magnetic complexes capable of retaining magnetization blocking for a long time at elevated temperatures, has been a major concern for magnetochemists over the last three decades. In this review, we describe basic SMMs and the different approaches that allow high magnetization-blocking temperatures to be reached. We focus on the basic factors affecting magnetization blocking, magnetic axiality and the height of the blocking barrier, which can be used to group different families of complexes in terms of their SMM efficiency. Finally, we discuss several practical routes for the design of mono- and polynuclear complexes that could be applied in memory devices.

Improving spin crossover characteristics in heteroleptic [FeIII(qsal-5-I)(qsal-5-OMe)]A complexes.

A family of heteroleptic spin crossover (SCO) [FeIII(qsal-5-I)(qsal-5-OMe)]A·sol (qsal-5-X = 5-X-2-[(8-quinolylimino)methyl]phenolate; A = NO3-1 sol = 2MeOH, NCS-2 sol = 0.75MeOH·1.3H2O, BF4-3 sol = MeOH, OTf-4, sol = MeOH) complexes have been synthesized. Most of the complexes exhibit gradual SCO, with the exception of NCS, which is principally high spin. In contrast, the OTf complex shows an abrupt hysteretic SCO (35 K) after solvent loss. The magnetic properties of this complex are significantly improved in comparison to the related homoleptics, [Fe(qsal-I)2]OTf 5 (hysteresis, 8 K) and [Fe(qsal-5-OMe)2]OTf·CH2Cl26 (gradual SCO). Structural studies reveal that slight changes in the crystal packing cause stronger interactions improving the cooperativity. These findings are supported by DFT calculations using the r2SCAN functional in which the calculated structures show that SCO from the LS to the HS state causes pronounced scissoring of the 1D π-π chains and substantial changes in their relative orientation following loss of MeOH.

Fe-Gd Ferromagnetic Cyclic Coordination Cluster [FeIII4GdIII4(teaH)8(N3)8(H2O)] with Magnetic Anisotropy─Theory and Experiment.

The synthesis, structural, and magnetic characterization of [FeIII4LnIII4(teaH)8(N3)8(H2O)] (Ln = Gd and Y) and the previously reported isostructural Dy analogue are discussed. The commonly held belief that both FeIII and GdIII can be regarded as isotropic ions is shown to be an oversimplification. This conclusion is derived from the magnetic data for the YIII analogue in terms of the zero-field splitting seen for FeIII and from the fact that the magnetic data for the new GdIII analogue can only be fit employing an additional anisotropy term for the GdIII ions. Furthermore, the Fe4Gd4 ring shows slow relaxation of magnetization. Our analysis of the experimental magnetic data employs both density functional theory as well as the finite-temperature Lanczos method which finally enables us to provide an almost perfect fit of magnetocaloric properties.

The spin magnetic order of Con+ (n ≤ 5) clusters.

The magnetism of transition metal clusters has been for decades a complicated puzzle, with experimental results disagreeing with calculations performed within the density functional theory formalism. In this work, we provide a key to this puzzle by investigating the lowest-energy spin states of cobalt cluster, Con+ (n ≤ 5), using CASSCF/NEVPT2 calculations with very large active spaces. The geometries as well as the spin configurations adopted by the clusters in their ground-state are known from experiments, making Con+ clusters an ideal model system for theoretical investigation. Here, using the experimentally known geometries determined by far-infrared spectroscopy as inputs, we calculated the lowest-energy spin configurations of the clusters, revealing that the CASSCF/NEVPT2 formalism correctly predicts the preferred electronic configuration of the clusters known experimentally. This is in contrast to the widely used density functional theory, with results that depend on the selected exchange-correlation functional. The reasons for the failure of density functional theory, in opposition to CASSCF/NEVPT2, are discussed, providing a solid framework for investigating other transition metal and transition metal oxide clusters.

Redox-Programmable Spin-Crossover Behaviors in a Cationic Framework.

Metal-organic frameworks (MOFs) provide versatile platforms to construct multi-responsive materials. Herein, by introducing the neutral tetradentate ligand and the linear dicyanoaurate(I) anion, we reported a rare cationic MOF [FeII(TPB){AuI(CN)2}]I·4H2O·4DMF (TPB = 1,2,4,5-tetra(pyridin-4-yl)benzene) with hysteretic spin-crossover (SCO) behavior near room temperature. This hybrid framework with an open metal site (AuI) exhibits redox-programmable capability toward dihalogen molecules. By means of post-synthetic modification, all the linear [AuI(CN)2]- linkers can be oxidized to square planar [AuIII(CN)2X2]- units, which results in the hysteretic SCO behaviors switching from one-step to two-step for Br2 and three-step for I2. More importantly, the stepwise SCO behaviors can go back to one-step via the reduction by l-ascorbic acid (AA). Periodic DFT calculations using various SCAN-type exchange-correlation functionals have been employed to rationalize the experimental data. Hence, these results demonstrate for the first time that switchable one-/two-/three-stepped SCO dynamics can be manipulated by chemical redox reactions, which opens a new perspective for multi-responsive molecular switches.

Insights into the Spin Dynamics of Mononuclear Cerium(III) Single-Molecule Magnets.

Four novel CeIII mononuclear complexes of formulas [Ce(ntfa)3(MeOH)2] (1), [Ce(ntfa)3(5,5'-Me2bipy)] (2), [Ce(ntfa)3(terpy)] (3), and [Ce(ntfa)3(bipy)2] (4), where ntfa = 4,4,4-trifluoro-1-(naphthalen-2-yl)butane-1,3-dionato, 5,5'-Me2bipy = 5,5'-dimethyl-2,2'-dipyridyl, terpy = 2,2':6',2″-terpyridine, and bipy = 2,2'-bipyridine, have been synthesized and structurally characterized with CeIII displaying coordination numbers of 8, 8, 9, and 10, respectively. Magnetic measurements indicate that all the complexes show a field-induced single-ion magnet behavior under a small applied dc field. The magnetic analysis shows the relevance of the different spin relaxation mechanisms in the magnetic relaxation of the CeIII compounds, with special emphasis on the local-mode process. Multiconfigurational calculations were also performed to get more information on the axiality of the compounds.

Dinuclear Fluoride Single-Bridged Lanthanoid Complexes as Molecule Magnets: Unprecedented Coupling Constant in a Fluoride-Bridged Gadolinium Compound.

A new synthetic method allows isolating fluoride-bridged complexes Bu4N{[M(3NO2,5Br-H3L1,1,4)]2(µ-F)} (M = Dy, 1; M = Ho, 2; M = Gd, 3) and Bu4N{[Dy(3Br,5Cl-H3L1,2,4)]2(µ-F)}·2H2O, 4·2H2O. The crystal structures of 1·5CH3C6H5,·2·2H2O·0.75THF, 3, and 4·2H2O·2THF show that all of them are dinuclear compounds with linear single fluoride bridges and octacoordinated metal centers. Magnetic susceptibility measurements in the temperature range of 2-300 K reveal that the GdIII ions in 3 are weakly antiferromagnetically coupled, and this constitutes the first crystallographically and magnetically analyzed gadolinium complex with a fluoride bridge. Variable-temperature magnetization demonstrates a poor magnetocaloric effect for 3. Alternating current magnetic measurements for 1, 2, and 4·2H2O bring to light that 4·2H2O is an SMM, 1 shows an SMM-like behavior under a magnetic field of 600 Oe, while 2 does not show relaxation of the magnetization even under an applied magnetic field. In spite of this, 2 is the first fluoride-bridged holmium complex magnetically analyzed. DFT and ab initio calculations support the experimental magnetic results and show that apparently small structural differences between 1 and 4·2H2O introduce important changes in the dipolar interactions, from antiferromagnetic in 1 to ferromagnetic in 4·2H2O.

Spin-Phonon Coupling and Slow-Magnetic Relaxation in Pristine Ferrocenium.

We report the spin dynamic properties of non-substituted ferrocenium complexes. Ferrocenium shows a field-induced single-molecule magnet behaviour in DMF solution while cobaltocene lacks slow spin relaxation neither in powder nor in solution. Multireference quantum mechanical calculations give a non-Aufbau orbital occupation for ferrocenium with small first excitation energy that agrees with the relatively large measured magnetic anisotropy for a transition metal S=1/2 system. The analysis of the spin relaxation shows an important participation of quantum tunnelling, Raman, direct and local-mode mechanisms which depend on temperature and the external field conditions. The calculation of spin-phonon coupling constants for the vibrational modes shows that the first vibrational mode, despite having a low spin-phonon constant, is the most efficient process for the spin relaxation at low temperatures. In such conditions, vibrational modes with higher spin-phonon coupling constants are not populated. Additionally, the vibrational energy of this first mode is in excellent agreement with the experimental fitted value obtained from the local-mode mechanism.

Magnetic anisotropy in YbIII complex candidates for molecular qubits: a theoretical analysis.

The magnetic properties of mononuclear YbIII complexes have been explored by using multiconfigurational CASPT2/RASSI calculations. Such complexes, in particular the case of [Yb(trensal)] complex, have been proposed as molecular qubits due to their spin dynamics properties. We have verified the accuracy of the theoretical approach to study such systems by comparing with experimental magnetic data. In order to have a wide overview of the magnetic properties of mononuclear YbIII complexes, we have considered simple charged and neutral models, [Yb(H2O)n]3+ and [Yb(OH)3(H2O)n-3], for many coordination modes. Thus, the results for more than 100 models allow extraction of some conclusions about the best ligand distributions in the coordination sphere to tailor the magnetic properties. Some low coordination, between 3 and 5, complexes that have no experimental magnetic data have been studied computationally to check if they can present high magnetic anisotropy.

Metal-Organic Nanocapsules with Functionalized s-Heptazine Ligands.

A metalloorganic capsule was synthesized where the ligand is a derivative of heptazine with three carboxylic groups that are coordinated to CuII cations, forming paddle-wheel motifs. Each nanocapsule is neutral, with 12 CuII centers and 8 ligands adopting a rhombicuboctahedron shape. It has almost 3 nm diameter, and the main intermolecular interactions in the solid are π··· π stacking between the C6N7 heptazine moieties. The nanocapsules can form monolayers deposited on graphite as observed by atomic force microscopy, which confirms their stability in solution.

Dysprosium-based complexes with a flat pentadentate donor: a magnetic and ab initio study.

The reactivity of the well-known pentadentate N3O2 Schiff base H2L (2,6-bis(2-hydroxyphenyliminomethyl)pyridine) towards a lanthanoid metal, in this case DyIII, has been investigated for the first time. This reactivity markedly depends on the pH of the medium and, accordingly, two different complexes, [Dy(HL)(NO3)2]·H2O (1·H2O) and [Dy(L)(NO3)(EtOH)(H2O)]·2H2O (2·2H2O), could be isolated from dysprosium(iii) nitrate and H2L. In addition, reaction of H2L with dysprosium(iii) chloride in methanol yields [Dy(HL')2][Dy(L)(Cl2)] (3), where H2L' ((6-(2-hydroxyphenyliminomethyl)-2-methoxyhydroxymethyl)pyridine) is an N2O2 hemiacetal donor derived from the partial hydrolysis of the H2L ligand, and subsequent addition of the methanol solvent to the carbonyl group. This latter reaction has been firstly observed for a lanthanoid metal. Single crystal X-ray diffraction studies of 1·1.15Py·0.3CH3C6H5, 2·2H2O and 3 show that the Schiff base is acting as a nearly flat pentadentate donor in all the cases, this behaviour being independent of the deprotonation degree of the phenolic oxygen atoms, both mono- or bisdeprotonated. Complexes 1·1.15Py·0.3CH3C6H5 and 2·2H2O show DyN3O6 cores, with distorted geometries closer to spherical tricapped trigonal prism or spherical capped square antiprism for 1·1.15Py·0.3CH3C6H5 and 2·2H2O, respectively. In the case of 3, the [Dy(HL')2]+ cation shows a dysprosium ion in an N4O4 triangular dodecahedron environment, while the [Dy(L)(Cl2)]- anion displays a DyN3O2Cl2 core with distorted pentagonal bipyramidal geometry. Moreover, attempts to dilute 1·H2O with yttrium yielded single crystals of (Et3NH)[Dy0.09Y0.91(L)(NO3)2] (4), where the Schiff base shows a similar pentadentate coordination mode. Dynamic magnetic studies of 1·H2O, 2·2H2O and 3 show that 2·2H2O and 3 present field-induced slow relaxation of the magnetisation, with Ueff barriers of 46.1(9) and 31.0(7) K for 2·2H2O and 3, respectively, while 1·H2O does not exhibit frequency-dependent peaks of the out of phase susceptibility, even in the presence of an external dc magnetic field. By contrast, the dilute sample 4 behaves as a SIM at zero dc field, with an energy barrier of ca. 49 K. Ab initio calculations using CASSCF methods including spin-orbit effects qualitatively support the obtained magnetic results, indicating that axiality is not the only factor that should be taken into account in order to increase effective energy barriers.

Single molecule magnets of cobalt and zinc homo- and heterometallic coordination polymers prepared by a one-step synthetic procedure.

The synthesis of 1D cobalt and zinc monometallic and heterometallic coordination polymers (CPs) was carried out applying one-pot synthetic methods by using either supercritical carbon dioxide or ethanol as the solvent. A collection of four 1D CPs were thus obtained by the combination of a metal (or a mixture of metals) with the linker 1,4-bis(4-pyridylmethyl)benzene. The used metallic complexes were zinc and cobalt hexafluoroacetylacetonate, which can easily incorporate pyridine ligands in the coordination sphere of the metal centre. Independently of the used solvent, the precipitated phases involving Zn(ii), i.e., homometallic CP of Zn(ii) and bimetallic CP of Zn(ii)/Co(ii), were isostructural. Contrarily, homometallic CPs of Co(ii) were precipitated as an isostructural phase of Zn(ii) or with a different structure, depending on the used solvent. All the structures were resolved by XRD using synchrotron radiation. In addition, the magnetic properties of the new CPs involving Co(ii) were studied. Remarkably, at low temperatures with the application of an external field, they acted as field-induced single molecule magnets.

Modular supramolecular dimerization of optically tunable extended aryl viologens.

Cucurbit[8]uril (CB[8]) mediated assembly of extended aryl viologens (EVs) into optically tunable dimers is reported for the first time. We show that the modular design and synthesis of a new class of π-conjugated viologen derivatives with rigid aromatic or heteroaromatic bridging units as well as electron donating molecular recognition motifs enable their self-assembly into 2 : 2 complexes with CB[8]. The quantitative dimerization process involving these two molecular components in an aqueous solution enables excimer-like interactions between closely packed charged guests giving rise to distinct spectroscopic behavior. The nature of these dimers (CB[8]2·(EV[X]R)2) in the ground and excited states was characterized by NMR, isothermal titration calorimetry, and steady-state spectroscopic measurements.

Dinuclear CoIIYIIIvs. tetranuclear CoY complexes: the effect of increasing molecular size on magnetic anisotropy and relaxation dynamics.

A new CoII2YIII2 complex with the formula [{Co(µ-L)Y(NO3)}2(µ-CO3)2]·2CH3OH·2H2O (where H2L = N,N',N''-trimethyl-N,N''-bis(2-hydroxy-3-methoxy-5-methylbenzyl)diethylenetriamine) has been prepared and its structure solved by single-crystal X-ray diffraction. The tetranuclear structure is formed by the connection of two [Co(µ-L)Y(NO3)] dinuclear units through two carbonate bridging ligands, which exhibit a µ3-κ2-O,O':κ-O:κ-O'' tetradentate coordination mode. The CoII ion exhibits a slightly distorted octahedral CoN3O3 coordination environment. From direct-current magnetic data a large and positive axial anisotropy parameter was extracted (D = +80.6 cm-1) and its sign unambiguously confirmed by HFEPR spectra and ab initio calculations. The extracted D value is rather larger than those previously reported for the analogous CoIIYIII dinuclear complexes, which agrees with the fact that the CoII ion in the CoII2YIII2 complex exhibits a lower distortion from the octahedral geometry in this family of CoIInYIIIn complexes. Dynamic ac magnetic measurements show that the reported compound presents field-induced slow relaxation for magnetization reversal, through a combination of direct and Raman processes. Magnetic measurements on the diluted magnetic counterpart (Zn/Co = 10/1) show the persistence of these processes, pointing out their single-ion origin. The Raman relaxation process for the CoII2YIII2 complex is faster that those observed for its CoIIYIII dinuclear counterparts. This fact and the existence of the persistent direct process at low temperature could be attributed to the former molecule being larger and more flexible than the latter ones.

Hexagonal Bipyramidal Dy(III) Complexes as a Structural Archetype for Single-Molecule Magnets.

Single-molecule magnets (SMMs), are regarded as excellent nanomaterials for high-density information storage and quantum computing. The local symmetry of the crystal field for the metal ion plays an important role in pursuing a high-performance SMM. Herein, two highly stable distorted hexagonal bipyramidal (quasi- D6 h) Dy complexes exhibiting slow relaxation of the magnetization are reported. A hexagonal bipyramidal Dy model complex with 18-crown-6 was also designed to study the relationship between magnetic anisotropy and symmetry. The combined experimental and theoretical results indicate that quantum tunneling is highly dependent on the local symmetries of the crystal field. The magnetic anisotropy becomes much stronger when the symmetry is closer to a standard D6 h geometry. These results support the conclusion that the hexagonal bipyramidal geometry is a viable one for the design of new classes of SMMs.

Anion-π Interactions in Computer-Aided Drug Design: Modeling the Inhibition of Malate Synthase by Phenyl-Diketo Acids.

Human infection by Mycobacterium tuberculosis (Mtb) continues to be a global epidemic. Computer-aided drug design (CADD) methods are used to accelerate traditional drug discovery efforts. One noncovalent interaction that is being increasingly identified in biological systems but is neglected in CADD is the anion-π interaction. The study reported herein supports the conclusion that anion-π interactions play a central role in directing the binding of phenyl-diketo acid (PDKA) inhibitors to malate synthase (GlcB), an enzyme required for Mycobacterium tuberculosis virulence. Using density functional theory methods (M06-2X/6-31+G(d)), a GlcB active site template was developed for a predictive model through a comparative analysis of PDKA-bound GlcB crystal structures. The active site model includes the PDKA molecule and the protein determinants of the electrostatic, hydrogen-bonding, and anion-π interactions involved in binding. The predictive model accurately determines the Asp 633-PDKA structural position upon binding and precisely predicts the relative binding enthalpies of a series of 2-ortho halide-PDKAs to GlcB. A screening model was also developed to efficiently assess the propensity of each PDKA analog to participate in an anion-π interaction; this method is in good agreement with both the predictive model and the experimental binding enthalpies for the 2-ortho halide-PDKAs. With the screening and predictive models in hand, we have developed an efficient method for computationally screening and evaluating the binding enthalpy of variously substituted PDKA molecules. This study serves to illustrate the contribution of this overlooked interaction to binding affinity and demonstrates the importance of integrating anion-π interactions into structure-based CADD.