Light-Induced Excited-State Spin Trapping in Tetrazole-Based Spin Crossover Systems

Ab initio calculation of magnetic anisotropy and thermal spin transition in the variable temperature crystal conformations of [Co(terpy)2]2

The structure-property correlation of [Co(terpy)2]2+, which shows a spin crossover at 270 K, has been computationally investigated based on its variable temperature crystal structures. Among the employed DFT functionals, only the re-parametrized hybrid B3LYP* functional could describe the correct spin transition temperature. Explicit consideration of metal-ligand sigma bonding with dynamic electron correlation is found to be necessary for an accurate determination of the SCO temperature with multi-reference calculations. The metal-ligand axial bond distances are found to be the most significant internal coordinates in deciding SCO. A small structural change along the axial distance causes a change in the t2g orbital splitting pattern and a reorientation of the magnetization axes at the SCO temperature. The complex shows an unusual triaxial magnetic anisotropy, with an easy axis of magnetization developing at higher temperatures. The strong coupling of low-frequency wagging motion of the two terpyridine ligands with the spin states of the complex provides an effective pathway for the relaxation of magnetization, resulting in a small magnetic anisotropy barrier.

Optical Phenomena in Molecule-Based Magnetic Materials

Since the last century, we have witnessed the development of molecular magnetism which deals with magnetic materials based on molecular species, i.e., organic radicals and metal complexes. Among them, the broadest attention was devoted to molecule-based ferro-/ferrimagnets, spin transition materials, including those exploring electron transfer, molecular nanomagnets, such as single-molecule magnets (SMMs), molecular qubits, and stimuli-responsive magnetic materials. Their physical properties open the application horizons in sensors, data storage, spintronics, and quantum computation. It was found that various optical phenomena, such as thermochromism, photoswitching of magnetic and optical characteristics, luminescence, nonlinear optical and chiroptical effects, as well as optical responsivity to external stimuli, can be implemented into molecule-based magnetic materials. Moreover, the fruitful interactions of these optical effects with magnetism in molecule-based materials can provide new physical cross-effects and multifunctionality, enriching the applications in optical, electronic, and magnetic devices. This Review aims to show the scope of optical phenomena generated in molecule-based magnetic materials, including the recent advances in such areas as high-temperature photomagnetism, optical thermometry utilizing SMMs, optical addressability of molecular qubits, magneto-chiral dichroism, and opto-magneto-electric multifunctionality. These findings are discussed in the context of the types of optical phenomena accessible for various classes of molecule-based magnetic materials.

Interpenetration Phenomena via Anion Template Effects in Fe(II) and Co(II) Coordination Networks with a Bis-(1,2,4-triazole) Ligand

Seven new coordination networks, [Fe(tbbt)3](BF4)2 (1), [Co(tbbt)3](BF4)2 (2), [Fe(tbbt)3](ClO4)2 (3), [Co(tbbt)3](ClO4)2 (4), [Fe(NCS)2(tbbt)2] (5), [Co(NCS)2(tbbt)2] (6), and [Fe(H2O)2(tbbt)2]Br2·2H2O (7), were synthesized with the linker 1,1'-(trans-2-butene-1,4-diyl)bis-1,2,4-triazole (tbbt) and structurally investigated. The structure of complexes 1-4 is composed of three interpenetrating, symmetry-related 3D networks. Each individual 3D network forms a primitive, nearly cubic lattice (pcu) with BF4- or ClO4- anions present in the interstitial spaces. The structure of compounds 5 and 6 is composed of two-dimensional sql layers, which are parallel to each other in the AB stacking type. These layers are interpenetrated by one-dimensional chains, both having the same formula unit, [M(NCS)2(tbbt)2] (M = Fe, Co). The structure of compound 7 consists of parallel, two-dimensional sql layers in the ABCD stacking type. The interpenetration in 1-6 is not controlled by π-π-interactions between the triazole rings or C=C bonds, as could have been expected, but by (triazole)C-H⋯F4B, C-H⋯O4Cl, and C-H⋯SCN anion hydrogen bonds, which suggests a template effect of the respective non-coordinated or coordinated anion for the interpenetration. In 7, the (triazole)C-H⋯Br anion interactions are supplemented by O-H⋯O and O-H⋯Br hydrogen bonds involving the aqua ligand and crystal water molecules. It is evident that the coordinated and non-coordinated anions play an essential role in the formation of the networks and guide the interpenetration. All iron(II) coordination networks are colorless, off-white to yellow-orange, and have the metal ions in the high-spin state down to 77 K. Compound 5 stays in the high spin state even at temperatures down to 10 K.

Toward High-Temperature Light-Induced Spin-State Trapping in Spin-Crossover Materials: The Interplay of Collective and Molecular Effects

Spin-crossover (SCO) materials display many fascinating behaviors including collective phase transitions and spin-state switching controlled by external stimuli, e.g., light and electrical currents. As single-molecule switches, they have been fêted for numerous practical applications, but these remain largely unrealized-partly because of the difficulty of switching these materials at high temperatures. We introduce a semiempirical microscopic model of SCO materials combining crystal field theory with elastic intermolecular interactions. For realistic parameters, this model reproduces the key experimental results including thermally induced phase transitions, light-induced spin-state trapping (LIESST), and reverse-LIESST. Notably, we reproduce and explain the experimentally observed relationship between the critical temperature of the thermal transition, T_1/2, and the highest temperature for which the trapped state is stable, T_LIESST, and explain why increasing the stiffness of the coordination sphere increases T_LIESST. We propose strategies to design SCO materials with higher T_LIESST: optimizing the spin-orbit coupling via heavier atoms (particularly in the inner coordination sphere) and minimizing the enthalpy difference between the high-spin (HS) and low-spin (LS) states. However, the most dramatic increases arise from increasing the cooperativity of the spin-state transition by increasing the rigidity of the crystal. Increased crystal rigidity can also stabilize the HS state to low temperatures on thermal cycling yet leave the LS state stable at high temperatures following, for example, reverse-LIESST. We show that such highly cooperative systems offer a realistic route to robust room-temperature switching, demonstrate this in silico, and discuss material design rationale to realize this.

Assessing the Performances of CASPT2 and NEVPT2 for Vertical Excitation Energies

Methods able to simultaneously account for both static and dynamic electron correlations have often been employed, not only to model photochemical events but also to provide reference values for vertical transition energies, hence allowing benchmarking of lower-order models. In this category, both the complete-active-space second-order perturbation theory (CASPT2) and the N-electron valence state second-order perturbation theory (NEVPT2) are certainly popular, the latter presenting the advantage of not requiring the application of the empirical ionization-potential-electron-affinity (IPEA) and level shifts. However, the actual accuracy of these multiconfigurational approaches is not settled yet. In this context, to assess the performances of these approaches, the present work relies on highly accurate (±0.03 eV) aug-cc-pVTZ vertical transition energies for 284 excited states of diverse character (174 singlet, 110 triplet, 206 valence, 78 Rydberg, 78 n → π*, 119 π → π*, and 9 double excitations) determined in 35 small- to medium-sized organic molecules containing from three to six non-hydrogen atoms. The CASPT2 calculations are performed with and without IPEA shift and compared to the partially contracted (PC) and strongly contracted (SC) variants of NEVPT2. We find that both CASPT2 with IPEA shift and PC-NEVPT2 provide fairly reliable vertical transition energy estimates, with slight overestimations and mean absolute errors of 0.11 and 0.13 eV, respectively. These values are found to be rather uniform for the various subgroups of transitions. The present work completes our previous benchmarks focused on single-reference wave function methods ( J. Chem. Theory Comput. 2018, 14, 4360; J. Chem. Theory Comput. 2020, 16, 1711), hence allowing for a fair comparison between various families of electronic structure methods. In particular, we show that ADC(2), CCSD, and CASPT2 deliver similar accuracies for excited states with a dominant single-excitation character.

CASPT2 molecular geometries of Fe(II) spin-crossover complexes

Using fully internally contracted (FIC)-CASPT2 analytical gradients, geometry optimizations of spin-crossover complexes are reported. This approach is tested on a series of Fe(II) complexes with different sizes, ranging from 13 to 61 atoms. A combination of active space and basis set choices are employed to investigate their role in determining reliable molecular geometries. The reported strategy demonstrates that a wave function-based level of theory can be used to optimize the geometries of metal complexes in reasonable times and enables one to treat the molecular geometry and electronic structure of the complexes using the same level of theory. For a series of smaller Fe(II) SCO complexes, strong field ligands in the LS state result in geometries with the largest differences between DFT and CASPT2; however, good agreement overall is observed between DFT and CASPT2. For the larger complexes, moderate sized basis sets yield geometries that compare well with DFT and available experimental data. We recommend using the (10e,12o) active space since convergence to a minimum structure was more efficient than with truncated active spaces despite having similar Fe-ligand bond distances.

Quantum dynamics simulations of the thermal and light-induced high-spin to low-spin relaxation in Fe(bpy)3 and Fe(mtz)6

TBC

Spin Crossover of Thiophosgene via Multidimensional Heavy-Atom Quantum Tunneling

The spin-crossover reaction of thiophosgene has drawn broad attention from both experimenters and theoreticians as a prime example of radiationless intramolecular decay via intersystem crossing. Despite multiple attempts over 20 years, theoretical predictions have typically been orders of magnitude in error relative to the experimentally measured triplet lifetime. We address the T₁ → S₀ transition by the first application of semiclassical golden-rule instanton theory in conjunction with on-the-fly electronic-structure calculations based on multireference perturbation theory. Our first-principles approach provides excellent agreement with the experimental rates. This was only possible because instanton theory goes beyond previous methods by locating the optimal tunneling pathway in full dimensionality and thus captures "corner cutting" effects. Since the reaction is situated in the Marcus inverted regime, the tunneling mechanism can be interpreted in terms of two classical trajectories, one traveling forward and one backward in imaginary time, which are connected by particle-antiparticle creation and annihilation events. The calculated mechanism indicates that the spin crossover is sped up by many orders of magnitude due to multidimensional quantum tunneling of the carbon atom even at room temperature.

Spin-crossover complexes in nanoscale devices: main ingredients of the molecule-substrate interactions

Spin-crossover complexes embedded in nanodevices experience effects that are absent in the bulk that can modulate, quench and even suppress the spin-transition. In this work we explore, by means of state-of-the-art quantum chemistry calculations, different aspects of the integration of SCO molecules on active nanodevices, such as the geometry and energetics of the interaction with the substrate, extension of the charge transfer between the substrate and SCO molecule, impact of the applied external electric field on the spin-transition, and sensitivity of the transport properties on the local conditions of the substrate. We focus on the recently reported encapsulation of Fe(II) spin-crossover complexes in single-walled carbon nanotubes, with new measurements that support the theoretical findings. Even so our results could be useful to many other systems where SCO phenomena take place at the nanoscale, the spin-state switching is probed by an external electric field or current, or the substrate is responsible for the quenching of the SCO mechanism.

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)₂(NCS)₂] 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.

The Role of Vibrational Anharmonicity in the Computational Study of Thermal Spin Crossover

Spin crossover in transition metal complexes can be studied in great detail with computational chemistry. Over the years, the understanding has grown that the relative stability of high-spin (HS) versus low-spin (LS) states is a subtle balance of many factors that all need to be taken into account for a reliable description. Among the different contributions, the zero-point energy (ZPE) and the entropy play key roles. These quantities are usually calculated assuming a harmonic oscillator model for the molecular vibrations. We investigated the impact of including anharmonic corrections on the ZPE and the entropy and indirectly on the critical temperature of spin crossover. As test systems, we used a set of ten Fe(II) complexes and one Fe(III) complex, covering different coordination modes (mono-, bi-, and tri-dentate ligands), decreasing coordination number upon spin crossover, coordination by second- and third-row atoms, and changes in the oxidation state. The results show that the anharmonicity has a measurable effect, but it is in general rather small, and tendencies are not easily recognized. As a conclusion, we put forward that for high precision results, one should be aware of the anharmonic effects, but as long as computational chemistry is still struggling with other larger factors like the influence of the environment and the accurate determination of the electronic energy difference between HS and LS, the anharmonicity of the vibrational modes is a minor concern.

Spin trapping and flipping in FeCO through relativistic electron dynamics

Transition metal compounds are very versatile, and their characteristics can differ profoundly depending on their electronic structure. Compounds in which a spin transition from a low-spin to a high-spin state can be achieved through means of an optical excitation are particularly intriguing, as a controlled spin-flip opens promising avenues in areas such as sensing, information technology, molecular switches and energy technology. The fundamental mechanisms in spin crossover and spin transitions remain unanswered, due to the complexity of electronic structure and interplay of relativistic effects. Presented here is a new approach that allows the first direct study of spin flip dynamics through a mapping of spin-mixed to spin-pure states. The method is applied to FeCO and addresses the spin-flip dynamics during a spin transition. Wave packets that combine different spin states are generated through optical excitation and relevant mechanisms in optically triggered spin transitions are discussed.

Electron transport through a spin crossover junction. Perspectives from a wavefunction-based approach

The present paper reports the application of a computational framework, based on the quantum master equation, the Fermi's golden Rule, and conventional wavefunction-based methods, to describe electron transport through a spin crossover molecular junction (Fe(bapbpy) (NCS)₂, 1, bapbpy = N-(6-(6-(Pyridin-2-ylamino)pyridin-2-yl)pyridin-2-yl)-pyridin-2-amine). This scheme is an alternative to the standard approaches based on the relative position and nature of the frontier orbitals, as it evaluates the junction's Green's function by means of accurate state energies and wavefunctions. In the present work, those elements are calculated for the relevant states of the high- and low-spin species of 1, and they are used to evaluate the output conductance within a given range of bias- and gate-voltages. The contribution of the ground and low-lying excited states to the current is analyzed, and inspected in terms of their 2S + 1 Ms-states. In doing so, it is shown the relevance of treating not only the ground state in its maximum-Ms projection, as usually done in most computational-chemistry packages, but the whole spectrum of low-energy states of the molecule. Such improved representation of the junction has a notable impact on the total conductivity and, more importantly, it restores the equivalence between alpha and beta transport, which means that no spin polarization is observed in the absence of Zeeman splitting. Finally, this work inspects the strong- and weak-points of the suggested theoretical framework to understand electron transport through molecular switchable materials, identifies a pathway for future improvement, and offers a new insight into concepts that play a key role in spintronics.

Diffusion Monte Carlo Perspective on the Spin-State Energetics of [Fe(NCH)6](2.)

The energy difference between the high spin and the low spin state of the model compound [Fe(NCH)6](2+) is investigated by means of Diffusion Monte Carlo (DMC), where special attention is dedicated to analyzing the effect of the fix node approximation on the accuracy of the results. For this purpose, we compare several Slater-Jastrow and multireference Slater-Jastrow trial wave functions. We found that a Slater-Jastrow trial wave function constructed with the generalized Kohn-Sham orbitals from hybrid DFT represents the optimal choice. This is understood by observing that hybrid functionals account for the subtle balance between exchange and correlation effects and the respective orbitals accurately describe the ligand-metal hybridization as well as the charge reorganization accompanying the spin transition. Finally the DMC results are compared with those obtained by Hartree-Fock, DFT, CASSCF, and CASPT2. While there is no clear reference value for the high spin-low spin energy difference, DMC and high level CCSD(T) calculations agree within around 0.3 eV.

TD-DFT study of the light-induced spin crossover of Fe(iii) complexes

Two light-induced spin-crossover Fe(iii) compounds have been studied with time-dependent density functional theory (TD-DFT) to investigate the deactivation mechanism and the role of the ligand-field states as intermediates in this process.

Insight into the light-induced spin crossover of [Fe(bpy)3]2+ in aqueous solution from molecular dynamics simulation of d–d excited states

Lifetimes of triplet d–d states were evaluated through molecular dynamics simulations to gain insight into relaxation dynamics of aqueous [Fe(bpy)₃]²⁺.

Spin-crossover in phenylazopyridine-functionalized Ni–porphyrin: trans–cis isomerization triggered by π–π interactions

To elucidate the light-induced spin-crossover mechanism of the PAPy-functionalized Ni(ii)–porphyrin, a DFT/CASSCF/CASPT2 study has been performed to determine the most stable cis and trans conformers and to characterize the excitation that triggers the process.

Towards an accurate and computationally-efficient modelling of Fe(ii)-based spin crossover materials

A theoretical approach is proposed to accurately calculate the LS–HS energy gap of SCO complexes in the solid state.

Spectroscopic, structural and magnetic investigations of iron(II) complexes based on 1-isopropyl- and 1-isobutyl-substituted tetrazole ligands

Two partly new [Fe(intz)6](BF4)2 complexes with n = 3 and 4 (i3tz = 1-isopropyl-1H-tetrazole, i4tz = 1-isobutyl-1H-tetrazole) were synthesized and characterised by X-ray powder diffraction, magnetic susceptibility measurements, vibrational, electronic and 57Fe-Mößbauer spectroscopy as well as DSC. The [Fe(i3tz)6](BF4)2 complex was re-investigated and shows an incomplete spin transition at Tc ∼ 109 K, while the [Fe(i4tz)6](BF4)2 features a complete but rather gradual spin transition with T½ = 223 K. In the lack of suitable crystals of [Fe(intz)6](BF4)2 with n = 3 and 4, we synthesized crystals of [Ni(intz)6](BF4)2 with n = 3 and 4, determined their X-ray crystal structures, and proved them to be homeotypic with the respective Fe-complexes by X-ray powder diffraction. DSC measurements showed an endothermic peak for the [Fe(i3tz)6](BF4)2 around T = 260 K not corresponding to a spin transition and suggesting a structural phase transition at this temperature. A well-developed peak at Tp = 225 K matches the spin-transition temperature T½ = 223 K for [Fe(i4tz)6](BF4)2.

Theoretical Investigation of the Electronic Structure of Fe(II) Complexes at Spin-State Transitions

The electronic structure relevant to low spin (LS)↔high spin (HS) transitions in Fe(II) coordination compounds with a FeN₆ core are studied. The selected [Fe(tz)₆]²⁺ (1) (tz = 1H-tetrazole), [Fe(bipy)₃]²⁺ (2) (bipy = 2,2′-bipyridine), and [Fe(terpy)₂]²⁺ (3) (terpy = 2,2′:6′,2″-terpyridine) complexes have been actively studied experimentally, and with their respective mono-, bi-, and tridentate ligands, they constitute a comprehensive set for theoretical case studies. The methods in this work include density functional theory (DFT), time-dependent DFT (TD-DFT), and multiconfigurational second order perturbation theory (CASPT2). We determine the structural parameters as well as the energy splitting of the LS–HS states (ΔE_HL) applying the above methods and comparing their performance. We also determine the potential energy curves representing the ground and low-energy excited singlet, triplet, and quintet d⁶ states along the mode(s) that connect the LS and HS states. The results indicate that while DFT is well suited for the prediction of structural parameters, an accurate multiconfigurational approach is essential for the quantitative determination of ΔE_HL. In addition, a good qualitative agreement is found between the TD-DFT and CASPT2 potential energy curves. Although the TD-DFT results might differ in some respect (in our case, we found a discrepancy at the triplet states), our results suggest that this approach, with due care, is very promising as an alternative for the very expensive CASPT2 method. Finally, the two-dimensional (2D) potential energy surfaces above the plane spanned by the two relevant configuration coordinates in [Fe(terpy)₂]²⁺ were computed at both the DFT and CASPT2 levels. These 2D surfaces indicate that the singlet–triplet and triplet–quintet states are separated along different coordinates, i.e., different vibration modes. Our results confirm that in contrast to the case of complexes with mono- and bidentate ligands, the singlet–quintet transitions in [Fe(terpy)₂]²⁺ cannot be described using a single configuration coordinate.

Molecular spin crossover phenomenon: recent achievements and prospects

Recently we assisted a strong renewed interest in the fascinating field of molecular spin crossover complexes by (1) the emergence of nanosized spin crossover materials through direct synthesis of coordination nanoparticles and nanopatterned thin films as well as by (2) the use of novel sophisticated high spatial and temporal resolution experimental techniques and theoretical approaches for the study of spatiotemporal phenomena in cooperative spin crossover systems. Besides generating new fundamental knowledge on size-reduction effects and the dynamics of the spin crossover phenomenon, this research aims also at the development of practical applications such as sensor, display, information storage and nanophotonic devices. In this critical review, we discuss recent work in the field of molecule-based spin crossover materials with a special focus on these emerging issues, including chemical synthesis, physical properties and theoretical aspects as well (223 references).