# Thermally-Induced Spin Crossover and LIESST Effect in the Neutral [FeII(Mebik)2(NCX)2] Complexes: Variable-Temperature Structural, Magnetic, and Optical Studies ($X = S, Se$; Mebik = bis(1-methylimidazol-2-yl)ketone)

Siddhartha De, Lise-Marie Chamoreau, Hasnaa El Said, Yanling Li, Alexandrine Flambard, Marie-Laure Boillot, Subrata Tewary, Gopalan Rajaraman, Rodrigue Lescouëzec

► To cite this version:

Siddhartha De, Lise-Marie Chamoreau, Hasnaa El Said, Yanling Li, Alexandrine Flambard, et al. Thermally-Induced Spin Crossover and LIESST Effect in the Neutral [FeII(Mebik)2(NCX)2] Complexes: Variable-Temperature Structural, Magnetic, and Optical Studies ($X = S, Se$; Mebik = bis(1-methylimidazol-2-yl)ketone). Frontiers in Chemistry, 2018, 6, pp.326. 10.3389/fchem.2018.00326. hal-01867908

HAL Id: hal-01867908

https://hal.sorbonne-universite.fr/hal-01867908

Submitted on 4 Sep 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Thermally-Induced Spin Crossover and LIESST Effect in the Neutral $[Fe^{II}(Me_{2}bik)_{2}(NCX)_{2}]$ Complexes: Variable-Temperature Structural, Magnetic, and Optical Studies ($X = S, Se; Me_{2}bik = bis(1-methylimidazol-2-yl)ketone$)

OPEN ACCESS

Edited by:

Matilde Fondo,
Universidade de Santiago de
Compostela, Spain

Reviewed by:

Jose Manuel Seco,
University of the Basque Country
(UPV/EHU), Spain
Joan Cano,
Universitat de València, Spain
Emilio Pardo,
Universitat de València, Spain

*Correspondence:

Marie-Laure Boillot
marie-laure.boillot@u-psud.fr
Gopalan Rajaraman
rajaraman@chem.iitb.ac.in
Rodrigue Lescouëzec
rodrigue.lescouenezec
@sorbonne-universite.fr

Specialty section:

This article was submitted to
Inorganic Chemistry,
a section of the journal
Frontiers in Chemistry

Received: 29 May 2018

Accepted: 11 July 2018

Published: 21 August 2018

Citation:

De S, Chamoreau L-M, El Said H, Li Y, Flambard A, Boillot M-L, Tewary S, Rajaraman G and Lescouëzec R (2018) Thermally-Induced Spin Crossover and LIESST Effect in the Neutral $[Fe^{II}(Me_{2}bik)_{2}(NCX)_{2}]$ Complexes: Variable-Temperature Structural, Magnetic, and Optical Studies ($X = S, Se; Me_{2}bik = bis(1-methylimidazol-2-yl)ketone$). *Front. Chem.* 6:326.
doi: 10.3389/fchem.2018.00326

Siddhartha De<fnref n="1" />, Lise-Marie Chamoreau<fnref n="1" />, Hasnaa El Said<fnref n="1" />, Yanling Li<fnref n="1" />, Alexandrine Flambard<fnref n="1" />, Marie-Laure Boillot²*, Subrata Tewary³, Gopalan Rajaraman³* and Rodrigue Lescouëzec<fnref n="1" />*

Two new iron(II) neutral complexes of bis(1-methylimidazol-2-yl)ketone ($Me_{2}bik$) with molecular formula $[Fe^{II}(Me_{2}bik)_{2}(NCS)_{2}]$ (1) and $[Fe^{II}(Me_{2}bik)_{2}(NCSe)_{2}]$ (2) have been synthesized and characterized by magnetic measurements, single-crystal X-ray diffraction, and solid state UV-vis spectroscopy. The temperature dependent magnetic susceptibility measurements of crystalline samples of both compound show the occurrence of a gradual spin transition centered at $T_{1/2} = 260$ K and 326 K, respectively. The crystal structures of both compounds were determined at different temperatures, below and above the transition, in order to detect the structural changes associated with the spin transition. The main structural modifications, when passing from the low-spin to the high-spin form, consist of an important lengthening of the Fe-N(Mebik) and Fe-N (C-S/Se) distances (by ca. 0.20 and 0.18 Å, respectively) and a noticeable variation of the N-Fe-N angles, leading to a more distorted [Fe-N6] octahedron. The spin-transition phenomenon also affects the optical properties, with significant decrease of the intensity of the Metal-to-Ligand charge transfer band upon increasing the temperature. Finally, both complexes exhibit a light-induced excited spin-state trapping under laser light irradiation at low temperature. DFT calculations were also carried out on these complexes in order to rationalize the theoretically predicted magnetic and optical behavior with those of the experimental one. The results clearly highlights the dramatic alteration of the magneto-structural behavior of the tris-chelate $[Fe^{II}(Me_{2}bik)_{3}]^{+}$ spin-crossover complex upon substituting one Mebik with NCS and NCSe ligands.

Keywords: spin crossover, iron(II) complexes, photomagnetism, DFT, B3LYP/B3LYP*

<footnote n="1">¹ Sorbonne Université, CNRS, Institut Parisien de Chimie Moléculaire, Paris, France, ² Institut de Chimie Moléculaire et des Matériaux d'Orsay, CNRS, Bat 420, Univ. Paris-Sud, Université Paris-Saclay, Orsay, France, ³ Department of Chemistry, Indian Institute of Technology Bombay, Mumbai, India</footnote>
## INTRODUCTION

The rational design of switchable molecular materials featuring stimuli-responsiveness is attracting strong research efforts because of the potential use as molecular sensors, switches, actuators, or memories in future molecular devices (Sauvage and Amendola, 2001; Feringa, 2011; Shepherd et al., 2013; Ferrando-Soria et al., 2017). The spin-crossover (SCO) complexes represent an emblematic class of switchable molecular systems. In these complexes, the application of thermal variation, pressure, light-irradiation, electric, or magnetic fields can induce an electronic change between a low-spin (LS) and a high-spin (HS) state, which leads to important changes in both optical and magnetic properties (Halcrow, 2013). The Fe(II) SCO complexes are of particular interest as the electronic reorganization involves a transition between a diamagnetic low-spin state ($S = 0, t_{2g}^6$) and a paramagnetic high-spin state ($S = 2, t_{2g}^4 e_{g}^2$). The spin-state change is also accompanied by a structural reorganization, in particular with a significant lengthening of the Fe-N(ligand) distances (ca. 0.2 Å). The important structural reorganization have been shown to favor the Light-Induced Spin-State Trapping (LIESST) and photomagnetic effects are often observed in Fe(II) SCO complexes (Létard, 2006).

In recent years, we have been interested by the use of the bis(1-R-imidazol-2-yl)ketone ligands (with R = methyl, ethyl, vinyl), named “Rbik,” for designing switchable photomagnetic molecules (Scheme 1). From a structural point of view, the Rbik ligands are β di-imine, which can be compared to the well-known α di-imine ligands: L-L = 2,2'-bipyridine (bipy) and 1,10-phenanthroline (phen). The presence of the ketone group between the two imidazolyl donor groups confers some aromaticity to these ligands, which are thus π-acceptor ligands, like the bipy and phen. However Rbik ligands form six-membered chelate rings leading to weaker ligand field than for the five membered chelate rings observed in similar bipy and phen complexes. Thus, whereas $[Fe^{II}(L-L)_{3}]^{2+}$ are a low-spin complexes, the $[Fe^{II}(Rbik)_{3}]^{2+}$ are SCO complexes (De et al., 2018). The high-spin six-membered rings $[Fe^{II}(Rbik)_{3}]^{2+}$ complexes are also more labile at room temperature in comparison to the inert $[Fe^{II}(L-L)_{3}]^{2+}$, which makes the preparation of bis-chelate complexes more easy. Taking advantage of such feature, we synthesized and studied $[Fe(Rbik)_{2}(NC-)]_{2}$ SCO systems, where one of the Rbik ligand is replaced by two cyanido metallo-ligand, abbreviated (NC-). For example, we investigated two tetranuclear rhombus $\{Fe_{2}M_{2}\}$ (M = Fe, Mo) switchable complexes containing the $\{Fe^{II}(Rbik)_{2}(NC-))_{2}\}$ subunit (NC- represents cyanido metallo ligand: $[Fe^{III}(Tp)(CN)_{3}]^{-}$ or $[Mo^{V}(CN)_{8}]^{3-}$) (Mondal et al., 2012, 2014). Both complexes show gradual thermally-induced spin transition on the Fe(II) centers and photomagnetic effect (LIESST) at low temperature. In the $\{[Fe^{III}(Tp)(CN)_{3}]_{2}[Fe^{II}(Me_{2}bik)_{2}]_{2}\}$ $[Fe^{III}(Tp)(CN)_{3}]_{2}18H_{2}O \cdot 4CH_{3}OH$ compound, the transition is centered at ca. $T_{1/2} = 330$ K whereas in the $\{[Mo(CN)_{8}]_{2}[Fe(Me_{2}bik)_{2}]_{2}(H_{2}MeIm)_{2}5H_{2}O \cdot CH_{3}CN$ compound, it occurs at a slightly higher temperature, ca. $T_{1/2} = 350$ K. In both cases the transition is significantly shifted toward lower temperature (ca. 250 K for the two compounds)

and becomes more gradual upon solvent removal. A striking difference between the two $\{Fe_{2}M_{2}\}$ rhombus complexes and the $[Fe^{II}(Me_{2}bik)_{3}](BF_{4})_{2}$ monometallic complex arises from the photomagnetic behavior, which was investigated using different laser diode in the visible and near IR range (404, 532, 605, 750, 808, 900 nm). Whereas the $\{Fe_{2}Fe_{2}\}$ compound show a significant LIESST effect in the 700–900 nm range (only in the dehydrated phase), the $\{Fe_{2}Mo_{2}\}$ shows the strongest increase of magnetization at 405 nm, whereas LIESST effect is much less significant at 800 nm. In contrast, the $[Fe(Me_{2}bik)_{3}](BF_{4})_{2}$ monometallic complex shows a maximal effect at 635 nm (De et al., 2018).

In order to further extend these seminal investigation, the present contribution focuses on the switchable properties of the mononuclear SCO complexes $[Fe^{II}(Me_{2}bik)_{2}(NCX)_{2}]$ ($X = S, Se$), which can be considered as model compounds for our previously reported systems as the group NC- is here a simple inorganic ligand. We thus present here a structural, spectroscopic, (photo)magnetic study of these two new spin crossover compounds: $[Fe^{II}(Me_{2}bik)_{2}(NCS)_{2}]$ (1) and $[Fe^{II}(Me_{2}bik)_{2}(NCSe)_{2}]$ (2).

<figure><img src="image_2.png" /><figcaption>SCHEME 1 | bis(1-methylimidazol-2-yl)ketone, $Me_{2}bik$.</figcaption></figure>

## RESULTS AND DISCUSSION

### Synthesis

Both compounds were prepared similarly, by simply reacting in-situ prepared $[Fe(Me_{2}bik)_{2}(S)_{2}]$ complex ($S = solvent$) with a stoichiometric amount of NCX salts in air and at room temperature (see details in the experimental part).

### Description of the Crystal Structure

The single crystal X-ray diffraction structure of 1 and 2 were determined at 200(2), 300(2), and 400(2) K. The compounds are isostructural and crystallize in the monoclinic C2/c space group (Z = 4). Selected crystallographic data are reported in Table 1 and Tables S1, S2.

The structures of 1 and 2 consist of mononuclear $[Fe(Me_{2}bik)_{2}(NCX)_{2}]$ complexes (Figure 1), which are linked by Van der Waals interactions and pseudo-hydrogen bonds. No solvent molecules are observed in the crystal lattice, which is coherent with the TGA data (see Supplementary information). In both compounds, the asymmetric unit is made of half a complex, which is related to the other half through a $C_{2}$ symmetry axis. A view of the crystal packing along the c axis is shown in Figure 2. In the complexes, the iron(II) atom is
<table><caption>**TABLE 1** | Cell parameters and volume cells of 1 and 2 obtained at 100, 200 and 400 K.</caption><thead><tr><th></th><th colspan="3">[Fe(Me<sup>bik</sup>)<sub>2</sub>(NCS)<sub>2</sub>]</th><th colspan="3">[Fe(Me<sup>bik</sup>)<sub>2</sub>(NCSe)<sub>2</sub>]</th></tr></thead><tbody><tr><td>T / [K]</td><td>200 (2)</td><td>300 (2)</td><td>400 (2)</td><td>200 (2)</td><td>300 (2)</td><td>400 (2)</td></tr><tr><td>a [Å]</td><td>7.7745 (2)</td><td>7.6476 (2)</td><td>7.6806 (3)</td><td>7.9156 (10)</td><td>7.9181 (10)</td><td>7.8489 (4)</td></tr><tr><td>b [Å]</td><td>13.7943 (3)</td><td>14.4341 (4)</td><td>14.5138 (6)</td><td>13.8415 (2)</td><td>13.9977 (2)</td><td>14.5090 (6)</td></tr><tr><td>c [Å]</td><td>22.3611 (5)</td><td>22.3945 (6)</td><td>22.5939 (10)</td><td>22.5668 (4)</td><td>22.6609 (4)</td><td>22.9089 (9)</td></tr><tr><td>β [°]</td><td>100.564 (1)</td><td>99.280 (2)</td><td>99.281 (2)</td><td>99.7010 (10)</td><td>99.333 (10)</td><td>98.819 (3)</td></tr><tr><td>Volume [Å<sup>3</sup>]</td><td>2357.45 (10)</td><td>2439.69 (11)</td><td>2485.68 (18)</td><td>2437.15 (6)</td><td>2478.38 (6)</td><td>2578.0 (2)</td></tr></tbody></table>

coordinated by two thio- or seleno-cyanate nitrogen atoms (for 1 or 2, respectively), which are in cis position and four Me\textit{bik} nitrogen atoms, leading to a distorted octahedral [N\textit{6}] coordination sphere. Selected interatomic distances and angles for 1 and 2 are listed in Tables S3, S4 and selected data of the coordination spheres for 1 and 2 are given in Table 2.

The comparison of the structures at different temperature points to the occurrence of a thermally-induced spin crossover in both compounds, which is also revealed by magnetometry measurements (vide infra). We examine first the structural data at 200 K.

The average Fe-N(NCS) and Fe-N(Me\textit{bik}) bond lengths measured at 200 K are equal respectively in 1 (1.967 Å) and very similar in 2 (respectively 1.959 and 1.962 Å). These distances are very homogeneous with an average deviation from the orthogonality of τ 0.006° (1) and 0.008° (2) (Table 2). These bond lengths are typical for low-spin octahedral Fe(II) complexes. For example, they are close to those previously reported in the low-spin [Fe(phen)<sub>2</sub>(NCS)<sub>2</sub>] (Real et al., 1992) and [Fe(bipy)<sub>2</sub>(NCS)<sub>2</sub>] (Konno and Mikami-Kido, 1991) complexes (Table 3). However in these last compounds, the bond deviation, τ, are higher because the Fe-NCX bonds are notably shorter. The distortion from the octahedral geometry in 1 and 2 is more apparent in the deviation from orthogonality of the twelve N-Fe-N angles involving the nitrogen donor atoms (see Table 2). The measured values, ca. 23°, are slightly above that measured in the low-spin [Fe(Me\textit{bik})<sub>3</sub>](BF<sub>4</sub>)<sub>2</sub> complex (13.5°). The distortions of the Fe(II) coordination sphere were also analyzed by continuous shape measurements using the SHAPE program (Casanova et al., 2005; Llunell et al., 2005). The reported output (shape factor) allows assessing the matching between an idealized polyhedron and actual coordination sphere: the lower the shape factor is, the better the matching between the actual coordination sphere and the idealized polyhedron. The analyses lead to numerical values of 0.85 for both 1 and 2 at 200 K (0 would correspond to a perfect octahedron), a value which is a bit larger than that observed in the low-spin [Fe(Me\textit{bik})<sub>3</sub>](BF<sub>4</sub>)<sub>2</sub> complex (0.42) (De et al., 2018).

Concerning the geometry of the ligands, the β-di-imine Me\textit{bik} bidentate ligands are not planar in contrast with the bipy and phen α-di-imine ligands. However we showed in a previous study that the presence of the CO group between the two imidazolyl groups confers some aromaticity to the Me\textit{bik} ligand and thus, a π-acceptor character. (De et al., 2018) Here the dihedral angle between the imidazolyl rings are 24.9 (1) and 24.1° (2) at 200 K.

These values are quite larger than those observed in the low-spin [Fe(Me\textit{bik})<sub>3</sub>](BF<sub>4</sub>)<sub>2</sub> (between ca. 4 and 11°) (De et al., 2018). Finally the NCX ligands are linear (Tables S3, S4) however the Fe-N-C angle deviates significantly from linearity (ca. 165° in 1 and 2) as previously observed in related [Fe(phen)<sub>2</sub>(NCX)<sub>2</sub>] complexes (Konno and Mikami-Kido, 1991; Real et al., 1992; MacLean et al., 2003).

The X-ray structure analyses carried out at higher temperature (300 and 400 K) reveal that complexes 1 and 2 undergo structural changes that are to be correlated with a spin-crossover process. The most obvious change is the significant increase of the metal-ligand bond lengths occurring upon heating. The average Fe-N coordination bonds increase by ca. 0.190 and 0.177 Å in 1 and 2, respectively, between 200 and 400 K. Such an increase is within the expected range for a spin crossover occurring in a {Fe<sup>II</sup>N<sub>6</sub>} system (Table 3) (Konno and Mikami-Kido, 1991; Real et al., 1992; MacLean et al., 2003). It is correlated to the increase of electronic population in the antibonding e<sup>*</sup> orbitals and the decrease of electronic population in the t<sub>2g</sub> orbital, which decreases the π-backbonding on going from the LS to the HS state. It is worth noticing that the distance variations are larger for the Fe-N(Me\textit{bik}) bonds than for the Fe-NCX ones. As suggested in previous studies, this may be correlated to the better π-acceptor character of the Me\textit{bik} ligand in comparison to the NCX (Konno and Mikami-Kido, 1991). However, it is not clear why this difference is less pronounced in the Me\textit{bik} complexes than in the α-di-imine complexes (Table 3). The significant increase of the Fe-N coordination bonds upon heating is also accompanied by an increase of the distortion in the coordination spheres, as revealed by the significant increase of the ∑ value and the S factor (Table 3). This behavior is characteristic of the conversion of low-spin Fe<sup>II</sup> state to a high-spin Fe<sup>II</sup> state.

The intermolecular interactions are also affected by the spin-state change, in particular, some intermetallic distances tend to decrease upon heating. Actually the analysis of the unit-cell parameters shows for both compounds a moderate thermal volumic expansion (>5%) that is related to an anisotropic and opposite variation of the axes. Beside the expansion along the b (5%) and c axes (1%), we observe a slight contraction along the a axis (1–2%). This contraction is correlated to the decrease of the smallest Fe···Fe distance (Table S5) between 200 and 400 K, from 7.774 to 7.681 Å in 1 and, from 7.916 to 7.849 Å in 2 (Table 1). This feature has already been observed in related [Fe(L-L)<sub>2</sub>(NCS)<sub>2</sub>] complexes (L-L are α-diimine ligands) and it is due to a scissor-like distortion of the complex
<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | View of the low-spin crystal structure of 1 (**Left**) and 2 (**Right**) at 200 K with the atom numbering of the coordination sphere and the NCX ligands. All hydrogen atoms are omitted for clarity (Fe: golden, C: gray, N: blue, O: red; S: yellow, Se: orange).</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**FIGURE 2** | View of the pseudo H-bond interaction (dotted green lines) in 1 at 200 K (view along *c* axis).</figcaption></figure>

during the spin-state change: the angle between the plans containing the L-L ligands decreases as the temperature increases (Shepherd et al., 2012).

In contrast, the pseudo hydrogen bonds between the S or Se atoms and a hydrogen atom of a $Me_{2}bik$ ligand belonging to a neighboring complex (Figure 2, Supplementary Figure S1),
<table><caption>**TABLE 2** | Selected structural data of the coordination sphere in 1 and 2 (Distance are given in Å and angles in degrees).</caption><thead><tr><th></th><th colspan="3">[Fe(Mebik)₂(NCS)₂]</th><th colspan="3">[Fe(Mebik)₂(NCSe)₂]</th></tr></thead><tbody><tr><td>T (K)</td><td>200</td><td>300</td><td>400</td><td>200</td><td>300</td><td>400</td></tr><tr><td>&lt;Fe-N(L)&gt;<sup>a</sup></td><td>1.967</td><td>2.141</td><td>2.169</td><td>1.960</td><td>1.989</td><td>2.140</td></tr><tr><td>Fe-NCX</td><td>1.967</td><td>2.108</td><td>2.133</td><td>1.961</td><td>1.974</td><td>2.132</td></tr><tr><td>&lt;N-Fe-N&gt;<sup>a</sup> bite angle</td><td>88.98</td><td>84.03</td><td>83.08</td><td>89.11</td><td>88.11</td><td>83.74</td></tr><tr><td>[Fe-N-C(X)]<sub>av</sub> angle</td><td>164.90</td><td>159.40</td><td>158.60</td><td>165.42</td><td>165.10</td><td>159.20</td></tr><tr><td>Σ|90 - α|<sup>b</sup></td><td>23.41</td><td>33.80</td><td>37.26</td><td>23.36</td><td>26.64</td><td>34.74</td></tr><tr><td>τ<sup>c</sup></td><td>0.0058</td><td>0.0317</td><td>0.0375</td><td>0.0075</td><td>0.0033</td><td>0.0083</td></tr><tr><td>Shape factor S (OC-6)</td><td>0.085</td><td>0.188</td><td>0.240</td><td>0.085</td><td>0.089</td><td>0.205</td></tr></tbody></table>

<table><caption>**TABLE 3** | Metal-Ligand distances in related Fe(II) SCO complexes.</caption><thead><tr><th rowspan="2">Compounds</th><th colspan="2">High temperature</th><th colspan="2">Low temperature</th><th colspan="2">ΔR upon HS-to-LS conversion</th><th rowspan="2">References</th></tr><tr><th>[Fe-NCX]</th><th>[Fe-N(L)]</th><th>[Fe-NCX]</th><th>[Fe-N(L)]</th><th>ΔR [Fe-NCX]</th><th>ΔR [Fe-N(L)]</th></tr></thead><tbody><tr><td>[Fe(phen)₂(NCS)₂]</td><td>2.057</td><td>2.206</td><td>1.958</td><td>2.009</td><td>0.099</td><td>0.197</td><td>Real et al., 1992</td></tr><tr><td>[Fe(bipy)₂(NCS)₂]</td><td>2.053</td><td>2.173</td><td>1.945</td><td>1.966</td><td>0.108</td><td>0.207</td><td>Konno and Mikami-Kido, 1991</td></tr><tr><td><sup>a</sup>[Fe(btz)₂(NCS)₂]</td><td>2.064</td><td>2.170</td><td>1.948</td><td>1.973</td><td>0.116</td><td>0.197</td><td>Real et al., 1992</td></tr><tr><td>[Fe(Mebik)₂(NCS)₂]</td><td>2.133</td><td>2.169</td><td>1.967</td><td>1.967</td><td>0.166</td><td>0.202</td><td>This paper</td></tr><tr><td>[Fe(Mebik)₂(NCSe)₂]</td><td>2.132</td><td>2.140</td><td>1.961</td><td>1.960</td><td>0.171</td><td>0.180</td><td>This paper</td></tr><tr><td>[Fe(phen)₂(NCSe)₂]</td><td>2.080</td><td>2.187</td><td>1.933</td><td>1.982</td><td>0.147</td><td>0.205</td><td>MacLean et al., 2003</td></tr></tbody></table>

<table><caption>**TABLE 4** | Pseudo-hydrogen bonds (Å) in 1 and 2.</caption><thead><tr><th></th><th>r(vdw)<sup>a</sup></th><th>200(2)</th><th>300(2)</th><th>400(2)</th></tr></thead><tbody><tr><td>S(1)-H (7)</td><td>3.00</td><td>2.85</td><td>2.94</td><td>2.97</td></tr><tr><td>Se(1)-H (7)</td><td>3.10</td><td>2.96</td><td>2.99</td><td>3.02</td></tr></tbody></table>

increase when the temperature increases (Table 4). This might be correlated to the decrease of the Fe(1)-N(1)-C(1) angle (ca. 6° in both 1 and 2) between 200 and 400 K. This angle variation might allow in part the absorption of the strain caused by the increase of Fe-N bond lengths accompanied by the spin-transition (Pap et al., 2011a,b,c; Scheja et al., 2015). The bending also indicates that the Fe-N=C=X mesomeric form (X = S or Se) coexists with the Fe-N=C-X⁻, which is expected to be linear. This is also reflected in the C-S and C-Se bond distances, which are intermediate between a double and a triple bond (ca. 1.62 and 1.78 Å respectively).

Overall the structural changes in the coordination sphere of 1 and 2 point to a full (or almost complete) spin crossover between 200 and 400 K. However the main structural changes occur between 200 and 300 K in 1 whereas in 2 they occur between 300 and 400 K. These structural data indicate that the spin transition in 1 should be centered at a lower temperature (ca. < 300 K) than that in 2 (ca. > 300 K), which is consistent with the

weaker ligand field induced by the NCS ligand as compared to the NCSe one (Nakano et al., 2004; Ross et al., 2011; Klingele et al., 2012).

## Optical Properties

The solid-state absorption spectra were collected on crystallites of the compounds dispersed in KBr pellets within a temperature range of 125–375 K. Similarly to the [Fe(Mebik)₃]²⁺ spin-crossover complex, 1 and 2 exhibit an asymmetric broad absorption band centered near 615 nm, with several contributions included those centered at ~ 570 and 655 nm (1) and 550 and 650 (2) (Figure 3 and Supplementary Figure S2). These bands responsible for the deep blue color of the complexes are ascribed to metal-to-ligand charge transfer (MLCT) absorption as confirmed by TD-DFT calculations (see below). Upon heating, the spectral evolution is similar in both compounds: (i) a pseudo isosbestic point appears; (ii) a decrease of intensity and a high-energy shift of the MLCT absorption are observed. These features match to the trends generally observed in this class of spin-crossover complexes. (Mondal et al., 2014; De et al., 2018) In fact, the larger metal-ligand bond lengths in the high-spin FeII complexes result in a weaker overlap between metal-centered and ligand-centered orbitals as compared to the low-spin state. Therefore, the intensity of the MLCT band strongly decreases upon converting the LS into HS species.

<footnote><sup>a</sup> <i>Fe-N</i> > and <i>N-Fe-N</i> > are the averaged values over the different coordination bonds and bite angles, respectively. <sup>b</sup>Sum of the deviation from the orthogonality of the twelve N-Fe-N angles. <sup>c</sup>τ = Σ (<i>Fe-N</i> - <i>Fe-N</i> >).</footnote>

<footnote><sup>a</sup>btz = 2,2'-bi-4,5-dihydrothiazine; The [Fe-N] distances are averaged values.</footnote>

<footnote><sup>a</sup>Sum of van der Waals radii.</footnote>
<figure><img src="image_5.png" /><figcaption>FIGURE 3 | Temperature dependence of the UV-Vis spectra of 1 between 125 and 375 K.</figcaption></figure>

The intense bands observed below 450 nm are ascribed to intra-ligand charge transfer transitions. On the low energy side of the spectrum, no other bands were clearly evidenced in the solid-state but a very weak unresolved absorption is observed in solution (see Supplementary Figures S2, S3). This weak band could be due to a ligand-field transition.

## Magnetic Properties

The magnetic properties of both complexes were investigated in the 4.0–400 K range by measuring the temperature dependence of the χT product (χ being the molar magnetic susceptibility) upon heating and cooling (at 2 K/min). The χT curves of 1 and 2 exhibit sigmoidal shapes, which are typical of the occurrence of a spin-crossover. In 1, the transition occurs between 170 and 380 K with an estimated transition temperature T₁/₂ of ca. 250 K. The significant variation of the χT value (from 0.08 cm³ mol⁻¹ K at 150 K to 3.25 cm³ mol⁻¹ K at 400 K) points to a complete thermally-induced S = 0 (t₂₉g) ⇌ S = 2 (t₂₉g e₂g) spin crossover. The magnetic behavior of 2 follows a similar trend but the spin transition starts at higher temperature (ca. 240 K) and seems almost complete at 400 K. The χT value increases from ca. 0.10 cm³ mol⁻¹ K at 200 K up to ca. 3.30 cm³ mol⁻¹ K at 400 K. The estimated transition temperature is ca. T₁/₂ ≈ 330 K (Figure 4).

Overall these magnetic data agree well with the structural investigation at different temperature and the optical spectroscopic studies at variable temperature. The thiocyanate derivative (1) exhibits a LS ⇌ HS spin transition, occurring below room temperature whereas the observed transition in selenocyanate derivative is shifted above room temperature. This is coherent with a stronger ligand field in the [Fe(Mebik)₂(NCSe)₂] complex than in the [Fe(Mebik)₂(NCS)₂] one. This agrees well with previous studies on related α-diimine-based complexes such as [Fe(phen)₂(NCS)₂] and [Fe(phen)₂(NCSe)₂] (Real et al., 1992; MacLean et al., 2003).

<figure><img src="image_6.png" /><figcaption>FIGURE 4 | Temperature Plot of χT vs. T for 1 and 2 under an applied magnetic external field of 1 T.</figcaption></figure>

In these compounds, the spin transition was quite abrupt, whereas in the present case the transition is much more gradual. This difference can be explained by the weaker intermolecular interactions involved in 1 and 2 as compared to those observed in [Fe(phen)₂(NCX)₂], which show π-π interaction beside the pseudo hydrogen bonds (Real et al., 2003; Gütlich et al., 2005) (MacLean et al., 2003).

The cooperativity between the spin-crossover complexes in the crystal lattice of 1 and 2 can be quantified by analyzing the χT vs. T curves through a regular solution model (expressed by the following equation): (Slichter and Drickamer, 1972).

$$ \ln \left[ \frac{(1 - n_{HS})}{(n_{HS} - f_{HS})} \right] = \frac{[\Delta H + \Gamma(f_{HS} + 1 - 2n_{HS})]}{RT} - \frac{\Delta S}{R} \quad (1) $$

where ΔH and ΔS are the enthalpy and entropy variation induced by the spin transition, Γ parameter is the cooperativity factor associated with the spin crossover (fHS is the residual HS molar fraction at low temperature and nHS represents the HS molar fraction).

The analyses lead to ΔH = 22.2 (1) 20.1 kJ mol⁻¹ (2), ΔS = 85.3 (1) 61.7 (2) J K⁻¹ mol⁻¹, and Γ = 1.5 (1) 3.3 (2) kJ mol⁻¹, respectively. The enthalpy and entropy variations are somehow lower than those obtained in the [Fe(Mebik)₃](BF₄)₂ complex but they remain in the expected range for Fe(II) SCO complexes. ((Nakamoto et al., 2001; Martínez et al., 2009); (Scott et al., 2013; Kumar et al., 2015)).

In contrast, the Γ values are moderate but higher in 1 and 2 than in the [Fe(Rbik)₃](BF₄)₂ complexes where only very weak intermolecular interaction were found. The larger value measured in 2 agrees well with the crystallographic data mentioned above, which revealed a stronger pseudo H-bond in compound 2 in comparison to compound 1. Finally, it is worth noticing that both spin transition are perfectly reversible: the χT vs. T curve obtained upon cooling superimposed perfectly with
that measured upon heating. The absence of solvent molecule in the crystal lattices of **1** and **2** likely accounts for this reversibility.

## The Photomagnetic Behavior

The photo-sensitivity of both compounds was investigated at 20 K by using laser diodes (with a power of ca. 7 mW/cm²) in the range 400–1,300 nm (405, 532, 635, 808, 900, 1,313 nm). Both compounds show Light-Induced Excited Spin-State Trapping (LIESST) effect in almost the all energy range (Figure 5). However significant differences are observed depending on the irradiation wavelengths. In both cases, the strongest effect is obtained at 900 nm. LIESST effects involving near IR irradiations are not so common, however some examples have been reported (Glijer et al., 2008; Buron-Le Cointe et al., 2012; Marino et al., 2014). Irradiation at 635 nm in the center of the MLCT band also leads to a LIESST effect but with a lower conversion that could be due to a smaller light penetration. Surprisingly, the effect of the 808 nm laser diode is remarkably low and it leads to a very poor conversion as compared to the irradiation at lower (900 nm) and higher (635 nm) energy.

The relative efficiency of the different wavelengths is similar in both compounds, however the conversion rates are overall higher in **2** than in **1**. In particular the irradiation at 900 nm leads to an almost quantitative spin conversion in **2**: the $χ_{M}T$ value reaches 2.83 cm³ mol¹ K after 30 min in comparison to 1.46 cm³ mol⁻¹ K in **1**. When comparing these results to those obtained on our previously reported $[Fe^{II}(Me_{2}bik)_{2}(N-)]_{2}$ complexes, we observe that the nature of the N- donors groups allows to tune the energy of the most efficient wavelength promoting the LIESST effect. In the $[Fe(Me_{2}bik)_{3}](BF_{4})_{2}$ complex, the optimal LIESST effect is observed at 635 nm (De et al., 2018). When the N- donors are the metallo-ligands, $[Mo(CN)_{8}]^{3-}$ or $[Fe(Tp)(CN)_{3}]^{-}$ (in the cyanide-bridged rhombus molecules of formula $\{[Mo(CN)_{8}]_{2}[Fe(Me_{2}bik)_{2}]_{2}\}^{(HMeIm)}_{2} \cdot 5H_{2}O \cdot CH_{3}CN$, where $^{MeImH}^{+}$ = N-methyl-imidazolium cation, and $\{[Fe^{III}(Tp)(CN)_{3}]_{2}[Fe^{II}(Me_{2}bik)_{2}]_{2}\} \cdot [Fe^{III}(Tp)(CN)_{3}]_{2} \cdot 18H_{2}O \cdot 4CH_{3}OH$, the most efficient effect are observed at 405 and 808 nm, respectively (Mondal et al., 2013, 2014) These differences should be related to the different electronic structures of the compounds and the nature of the accessible excited-states involved in the LIESST effect. The occurrence of reverse LIESST effect in some of these compounds was also proposed to account for a partial decrease of the photomagnetic response (Moussa et al., 2007). The rationalizing of the wavelength dependence on the LIESST effect would be interesting but it requires a set of advanced physical measurements that is beyond the scope of the present report.

It is also worth noting that the relaxation temperature, $T_{LIESST}$, are significantly higher (ca. 60 K) in all the previously reported complexes based on the $[Fe^{II}(Me_{2}bik)(N-)]_{2}$ units in comparison to those observed in **1** and **2**. Here, both complexes exhibit a rapid relaxation if the light is switched off, revealing that the $T_{LIESST}$ are below 20 K. This unusual behavior might be associated to the existence of a different relaxation process in **1** and **2**. Recently, Heinze and McCusker showed that the existence of low-lying $^{3}T_{1}$ state might open efficient relaxation pathway from the MLCT state to the $^{1}A_{1}$ ground state in Fe$^{II}$ complexes containing

pyridine-based ligand connected by keto groups (Jamula et al., 2014; Mengel et al., 2015).

## Theoretical Calculations

As observed in our earlier studies on the $[Fe^{(R)bik}_{3}](BF_{4})_{2}$ complexes featuring spin crossover properties (De et al., 2018), the choice of a suitable exchange-correlation functional within DFT method which could accurately predict both the ground state structure and SCO characteristics is challenging. For example, although the spin-state energies ($ΔE_{HS-LS}$) is in accordance with the SCO properties, hybrid B3LYP functional fails to correctly predict the ground state. In contrast, a modified version where percentage Hartree-Fock (%HF) exchange is scaled to 15% (B3LYP*) and OPBE predict the ground state correctly but the predicted spin-state splitting are too high to suggest SCO characteristics. Cooperative effects in SCO molecules play a very important role and there are methods available within the DFT to estimate such effects comprehensively (Rackwitz et al., 2013; Scott et al., 2013; Paulsen, 2016). In this present work, we have probed the energetic of spin-state rearrangements and the excited state properties within the intramolecular framework. The PCM (Tomasi et al., 2005) solvent model was employed to counter the effects of solvent molecules around the metal complex. While this qualitatively address the SCO effect in solution, the effect observed in solid-state are still missing.

As %HF exchange is found to play a role in dictating the spin-state splitting, we decided to vary the %HF exchange in B3LYP to assess and understand how this influence the SCO characteristics of $[Fe(bik)_{2}(NCX)_{2}]$ complexes. Additionally for comparative purpose, calculations are also performed using OPBE functional. The optimized energies of complexes **1** and **2** (Table 5) reveal that the B3LYP functional predicts a high-spin ground state in both the cases, whereas all other variation attempted (BLYP, B3LYP*, OPE, etc.) suggest a low spin as the ground state. The difference in energy between the high-spin and the low-spin state defined as $ΔE_{HS-LS}$ is computed varying the % HF exchange using BLYP functional. A linear dependency of $ΔE_{HS-LS}$ on the %HF exact exchange is observed for both complexes **1-2** (Figure 6) as reported earlier for other SCO complexes. (Grimme, 2006; Bowman and Jakubikova, 2012; Kepp, 2013) The $ΔE_{HS-LS}$ energies in case of B3LYP and BLYP (with 25%HF exchange) are estimated to be -12.5 and -37.3 kJ/mol in **1**, and -9.2 and -34.2 kJ/mol in **2**. Calculations reveal that as the % HF exchange increases, the $ΔE_{HS-LS}$ gap decreases with HS becoming ground state from 20% HF onwards. The $ΔE_{HS-LS}$ energies summarized in these cases as, 15.1 (15%), 45.8 (10%), and 70.4 (0%) kJ/mol in **1**, whereas the same estimates are 19.0 (15%), 49.5 (10%), and 73.1 (0%) kJ/mol in **2**. On the other hand, OPBE functional estimated the $ΔE_{HS-LS}$ energies are, 14.0 and 19.2 kJ/mol for complexes **1** and **2**, respectively. It is important to note that the $ΔE_{HS-LS}$ energies are slightly higher in **2** as compare to complex **1**. This is in agreement with the greater covalency associated with the selenium atom as compare to sulfur. Lastly the computed energies using B3LYP (%HF = 20), B3LYP* (%HF = 15) and OPBE functional estimates the $ΔE_{HS-LS}$ gap below 20 kJ/mol for both **1** and **2**. These values are coherent
<figure><img src="image_7.png" /><figcaption>FIGURE 5 | $\chi_{MT}$ vs. irradiation time of 1 (Left) and 2 (Right) measured at 20 K for different wavelengths (ca. 5 mW/cm$^{2}$).</figcaption></figure>

<table><caption>TABLE 5 | DFT computed $\Delta E_{HS-LS}$ (in kJ/mol) obtained with different %HF exchange values (0–25%) using BLYP and including OPBE.</caption><thead><tr><th></th><th>HF-0%</th><th>HF-10%</th><th>HF-15%</th><th>HF-20%</th><th>HF-25%</th><th>OPBE</th></tr></thead><tbody><tr><td>1</td><td>70.4</td><td>45.8</td><td>15.1</td><td>-12.5</td><td>-37.3</td><td>14.0</td></tr><tr><td>2</td><td>73.1</td><td>49.5</td><td>19.0</td><td>-9.2</td><td>-34.2</td><td>19.2</td></tr></tbody></table>

with the occurrence of a spin-crossover phenomena (Ye and Neese, 2010b) in both complexes (1 and 2) and in line with the experimental observations where transition are measured at $T_{1/2}$ of 260 and 326 K for 1 and 2, respectively.

To understand how temperature influences the geometrical parameters, we have performed geometry optimization starting from the X-ray structures of 1 and 2 obtained at different temperatures. Due to the difference in the $T_{1/2}$ observed, it is clear from the estimated structural parameters that the high-spin optimized structure at 300 K of 1 is identical with the one computed at 400 K, while the 300 K structure in case of 2 is identical with the low-spin optimized structure at 200 K (Table 2). Here the optimized spin-state energies ($\Delta E_{HS-LS}$) of 1 and 2 are considered for discussion.

The optimized geometries of the low-spin and high-spin geometries of 1 and 2 are shown in Figure 7 and some selected bond parameters are summarized in Table 6. The choice of DFT functional plays an important role on the optimized bond parameters. Here using B3LYP functionals, metal-ligand bond distances (in the structure optimized in gas phase at 0 K) are overall in very good agreement with the experimental values (Table 2). They appear slightly overestimated as compare to the experimental X-ray structures, in particular for the Fe-N(bik) distances. The significant increase of the Fe-N bond length upon heating is well reproduced by the theoretical calculation. In fact the bond distance variation associated with the Fe-N(bik) and

Fe-NCX bonds are estimated to be 0.22 and 0.15 Å in 1, and 0.21 and 0.16 Å in 2, respectively. These values are close to the experimental ones (Table 2), in particular, the larger variation in the Fe-N(bik) bond lengths as compare to the Fe-NCX bonds is also observed in the optimized structure. This further supports the fact that the $Me_{2}bik$ ligands offer suitable $\pi$-acceptor pathway as compare to the NCX ligand.

The distortion of the octahedral coordination sphere observed in the experimental X-ray structures is also reflected in the calculated ones. For example, the computed axial ($N-Fe-N)_{av}$ angles decrease by ca. 5° between the LS and HS states, as observed in the experiment (Table 2). It is worth noting that optimized geometries are more symmetric compare to the X-ray structures for 1 and 2. For example the average axial and bite ($N-Fe-N)_{av}$ bond angles in the computed structures are close to the 180° and 90°, respectively. This reveals that solid-state effects contribute to the distortion in the coordination sphere in the actual structure. This is particularly clear for the bending of the NCX ligands, which are involved in pseudo hydrogen interaction in the crystal lattice: whereas the computed Fe-N-C angle is close to linearity, the experimental ones vary between ca. 165° (LS state) and ca. 159° (HS state).

## Thermodynamic Parameters and Molecular Orbital Analysis

Thermodynamic parameters estimated from the frequency calculations on top of the optimized coordinates are summarized in Table 7. The contribution from vibrational modes is factored in the energy as zero-point corrected energies are used for this purpose. The high-spin state is always favored by the vibronic entropy; i.e., strong field ligands that energetically favors the low-spin state suffers larger loss to the vibronic entropy. The estimates of the calculated entropy ($\Delta S$) and enthalpy ($\Delta H$) are 88.3 and 84.9 J/mol-K, and -18.0 and -14.2 kJ/mol (B3LYP)
<figure><img src="image_8.png" /><figcaption>**FIGURE 6** | DFT computed spin-state energies ($\Delta E_{HS-LS}$) vs. HF exchange percentage employing BLYP functional for complexes 1 and 2.</figcaption></figure>

for 1-2, respectively. These estimates of the calculated entropy are on a higher note but in fair agreement with the experimental data derived from the magnetic measurements (see above). The calculated $\Delta H$ values are slightly underestimated as compared to the experimental data obtained on the solid-state measurements. They are also relatively smaller as compared to the estimates from $[Fe(Rbik)_{3}](BF_{4})_{2}$ complexes (De et al., 2018), and these further manifests that the $^{Rbik}$ligand exerts stronger ligand fields as compared to the NCX ($X = S/Se$) ligand due to relatively stronger $\pi$-acceptor abilities. The OPBE calculated values of $\Delta S$ and $\Delta H$ parameters follows the footsteps of the B3LYP calculated results whereas the B3LYP* computed results are unsatisfactory. The differences in the entropy ($T\Delta S$) and the enthalpy ($\Delta H$) contributions sum up to lead to the difference in the values of free energy component ($\Delta G_{el}$), 43.9 kJ/mol and 39.3 kJ/mol for 1 and 2 respectively. Calculated $\Delta E_{vib}$ values on the other hand are on the range of typically iron(II) octahedral complexes as observed earlier (De et al., 2018).

The nature of the spin ground state is determined by the orbital splitting and the spin pairing energy. Therefore, to gain further insights into the nature of the bonding, the d-based orbitals of complexes 1 and 2 at their respective high-spin state (400 K) were plotted using the eigenvalue orbital energies (Figure 8). For 1 and 2 HS structures, the B3LYP calculated electronic configuration yields the following, $(d_{xz})^{2}(d_{xy})^{1}(d_{yz})^{1}(d_{x2-y2})^{1}(d_{z2})^{1}$. The splitting of d orbitals is estimated to be 2.685 and 2.541 eV for complexes 1 and 2, respectively. These values are much higher as compared to $[Fe(Rbik)_{3}](BF_{4})_{2}$ (~1.5 eV) complexes (De et al., 2018). Beside, the $t_{2g}$ set of bonding orbitals are non-degenerate in 1 and 2 in contrast with the previous observations. This is due to the different $\sigma$-bonding interaction pattern of the $^{Rbik}$and the NCX ligands. The splitting patterns for complexes 1 and 2 are consistent with each other. There is a slight difference in the energy of crystal splitting due to the greater covalency of the

metal-ligand bonds in the selenium derivative as compared to the sulfur one. The extended delocalization of $\pi$-clouds from the $d_{xz/yz}$ orbitals of the iron(II) center to the imidazole group is significantly high in the NCSe analog (2) as compared to the NCS analog (1) which can also be identified from the shorter Fe-N($^{Me}_{2}bik$) bond distances and the larger computed $\Delta E_{HS-LS}$ energy gap. Finally, it is worth noticing that there is no C-H···O interactions between the methyl hydrogen and the ketone oxygen atoms, in contrast with the situation in the $[Fe(Rbik)_{3}](BF_{4})_{2}$ complexes (De et al., 2018). This is correlated to the weaker $\pi$-bonding interaction in $[Fe(Rbik)_{2}(NCX)_{2}]$ complexes as compared to the $[Fe(Rbik)_{3}](BF_{4})_{2}$ analogs. Hence the hydrogen atoms of the C-H bond of the methyl group in complexes 1 and 2 are less polarized to promote such interaction.

The spin density calculations also reflect somehow the bonding interaction in the complexes. The iron(II) center of 1 and 2 high spin structures shows almost identical spin density values of 3.830 and 3.826, respectively (Figure 9). The smaller spin-density value of the iron(II) in 2 reflects the relatively larger delocalization (back bonding) as compared to 1. The net spin densities at the N-donor center for N-bik and NCX ligand are 0.020 and 0.013 for both 1 and 2, respectively. The larger spin densities at the N-bik reflects the greater $\pi$-acceptor ability of the organic ligand as compared to the NCX ligands. Although the plotted spin density isosurface of complex 1-2 and $[Fe(Rbik)_{3}](BF_{4})_{2}$ complexes (De et al., 2018) are nearly identical and the quoted spin density values are small, the difference observed between the two set of complexes clearly indicates significant structural and electronic alteration upon ligand modifications owing to variations in the SCO properties.

## Optical Properties

In order to rationalize the optical properties of complexes 1 and 2, TD-DFT calculations using B3LYP functional have been
<figure><img src="image_9.png" /><figcaption>**FIGURE 7** | B3LYP optimized structures along with some selected bond parameters of complexes 1-2 at 200 and 400 K. Colour code: Blue = N; Red = O; Yellow = S/Se; Gray = Carbon. Hydrogen atoms are omitted for clarity.</figcaption></figure>

performed on the X-ray structures. The computed absorption spectra are shown in **Figure 10**, **Supplementary Figure S4**. Calculated spectra of complexes **1** and **2** in the low-spin state show intense bands near ca. 611 nm and 596 nm, respectively, which are consistent with the experimentally observed intense band that were presumably assigned to MLCT bands (see above). Actually the theoretical calculations reveal that these bands consist of a series of transitions LLCT and MLCT as it has been noted earlier in other Fe(II) low-spin complexes. (Mengel et al., 2015) The different contributions of the absorption spectra of **1** and **2** are listed in **Table 8** and the corresponding orbital diagrams are given in **Supplementary Figure S6** in ESI.

## CONCLUSION

The current report presents a straightforward synthesis of two new air-stable Fe(II) spin-crossover complexes based on β-diimine ligands, namely [Fe(Mebik)₂(NCS)₂], **1**, and

<table><caption>**TABLE 6** | Selected bond parameters based on B3LYP optimized geometries for 1 and 2.</caption><thead><tr><th rowspan="2">B3LYP<br>T (K)</th><th colspan="2">[FeMe(bik)₂(NCS)₂]</th><th colspan="2">[FeMe(bik)₂(NCSe)₂]</th></tr><tr><th>LS</th><th>HS</th><th>LS</th><th>HS</th></tr></thead><tbody><tr><td>d[Fe-N(L)]<sub>av</sub> (Å)</td><td>2.022</td><td>2.237</td><td>2.021</td><td>2.235</td></tr><tr><td>d[(Fe-NCX)]<sub>av</sub>(Å)</td><td>1.985</td><td>2.139</td><td>1.980</td><td>2.144</td></tr><tr><td>axial (N-Fe-N)<sub>av</sub> angle (°)</td><td>179.4</td><td>172.9</td><td>179.9</td><td>175.3</td></tr><tr><td>bite (N-Fe-N)<sub>av</sub> angle (°)</td><td>87.9</td><td>86.6</td><td>87.9</td><td>86.9</td></tr><tr><td>[Fe-N-C(X)]<sub>av</sub> angle (°)</td><td>178.2</td><td>174.2</td><td>178.7</td><td>175.8</td></tr></tbody></table>

[Fe(Mebik)₂(NCSe)₂], **2**, which are reminiscent of the well-known [FeII(phen)₂(NCX)₂] complexes. The magnetic studies reveal the occurrence of weakly cooperative spin transition. This is coherent with the structural analyses, which show only weak intermolecular interactions (pseudo-hydrogen bonds) between the SCO complexes. The transitions are perfectly reversible upon heating and cooling, which is in agreement with
<table><caption>TABLE 7 | Comparison of the calculated thermodynamic parameters of 1 and 2.</caption><thead><tr><th rowspan="2">Δ = HS - LS<sup>[b]</sup></th><th colspan="3">[Fe(bik)<sub>2</sub>(NCS)<sub>2</sub>]</th><th colspan="3">[Fe(bik)<sub>2</sub>(NCSe)<sub>2</sub>]</th></tr><tr><th>B3LYP</th><th>B3LYP*</th><th>OPBE</th><th>B3LYP</th><th>B3LYP*</th><th>OPBE</th></tr></thead><tbody><tr><td>ΔH<sub>el</sub> (kJ/mol)</td><td>-18.0</td><td>62.3</td><td>12.9</td><td>-14.2</td><td>70.7</td><td>12.9</td></tr><tr><td>ΔH<sub>el</sub> <sup>[a]</sup> (kJ/mol)</td><td>-20.1</td><td>7.1</td><td>7.9</td><td>-10.0</td><td>35.1</td><td>12.9</td></tr><tr><td>ΔG<sub>el</sub> (kJ/mol)</td><td>-43.9</td><td>53.9</td><td>-17.1</td><td>-39.3</td><td>39.7</td><td>-13.1</td></tr><tr><td>ΔG<sub>el</sub> <sup>[a]</sup> (kJ/mol)</td><td>-42.8</td><td>-15.1</td><td>-22.6</td><td>-33.4</td><td>11.7</td><td>-13.1</td></tr><tr><td>ΔE<sub>vib</sub> (kJ/mol)</td><td>-5.4</td><td>-7.9</td><td>-6.2</td><td>-5.0</td><td>-3.7</td><td>-6.2</td></tr><tr><td>ΔE<sub>vib</sub> <sup>[a]</sup> (kJ/mol)</td><td>-8.4</td><td>-8.8</td><td>-5.8</td><td>-3.3</td><td>13.8</td><td>-6.2</td></tr><tr><td>ΔS (J/mol·K)</td><td>88.3</td><td>88.1</td><td>101.7</td><td>84.9</td><td>103.7</td><td>101.6</td></tr><tr><td>ΔS <sup>[a]</sup> (J/mol·K)</td><td>74.5</td><td>74.9</td><td>102.5</td><td>78.6</td><td>79.5</td><td>101.6</td></tr></tbody></table>

<figure><img src="image_10.png" /><figcaption>FIGURE 8 | Energy splitting of d-based orbital's for (A) complex 1 (B) complex 2 in the high spin structures.</figcaption></figure>

the absence of volatile solvent molecule in the crystal lattices of both compounds. Actually, both compounds do not significantly lose their crystallinity upon heating. This has allowed the structural determination of the LS and HS states in both cases, revealing the lengthening of the metal-ligand distances and the higher distortion of the coordination sphere in the HS state.

The spin transition are centered at 240 and 330 K, in 1 and 2 respectively, in agreement with the higher ligand field induced by the NCSe ligand compared to the NCS one. The theoretical calculations performed on the electronic structures of both low-spin and high-spin structures of 1 and 2 support the occurrence of a spin-transition. They also enlighten the bonding frame in the two complexes: the M-NCSe bond is

more covalent than the M-NCS one. In addition, the back-donation from the Fe to the $Me_{2}bik$ ligand is also more efficient in the NCSe complex and it leads to a stronger ligand field. Both effects favor a higher spin-transition temperature in the NCSe derivative, as previously observed in $[Fe^{II}(phen)_{2}(NCX)_{2}]$ derivatives. Several functionals and variation of % HF exchange were explored to accurately reproduce the HS-LS gap for these two complexes. These attempts reveal that lower % HF exchange such as the one present in B3LYP* functional lead to superior results for estimating the gap while it worsens the estimation of the thermodynamic parameters. In other words, a universal functional, which could accurately compute spin-state splitting, thermodynamic parameters and geometries is unfounded.

<footnote>[a] B3LYP-D2 dispersion corrected values; [b] difference between the HS and LS values computed are listed here.</footnote>
<figure><img src="image_11.png" /><figcaption>FIGURE 9 | B3LYP computed spin density plot for 1 (left) and 2 (right) in the high spin geometry.</figcaption></figure>

Interestingly, a LIESST effect is observed at low temperature in both complexes. However, in contrast with other related $[Fe^{II}(Rbik)_{2}(NC-)]^{2+}$ SCO complexes, the relaxation temperature are very low here, and the relaxation of the magnetization occurs when the light is switched off. The experiments also reveal a remarkable dependence of the photomagnetic effect on the irradiation wavelengths. In the present case, a LIESST effect is observed when irradiating the sample in the MLCT band (near 600 nm) but with a weak efficiency. The most efficient photo-switching of the magnetization occurs at 900 nm, in the near IR.

Finally, it is worth noting that the comparison of these results with previous photomagnetic measurements on related $[Fe(Rbik)_{2}(NC-)]^{2+}$ complexes show that the most efficient wavelength can be varied by changing the nature of the -N donor ligands or metallo-ligands. The rationalizing of this dependency is difficult. At that stage, advanced physical measurements would be useful, for example, to identify the nature of the excited states involved in the LIESST effect. This might be a necessary would be of very high value to gain control on the photomagnetic effect by rational chemical design.

# MATERIALS AND METHODS

## Synthesis

The $Me_{2}bik$ ligand was synthesized according to a literature procedure (Lucas et al., 2000). Commercially available chemicals including KSCN, KScCN, and $Fe(SO_{4})·7H_{2}O$ were purchased from commercial sources and used without prior purification.

## Preparation of $[Fe(Me_{2}bik)_{2}(NCS)_{2}]$ (1) and $[Fe(Me_{2}bik)_{2}(NCSe)_{2}]$ (2)

A solution of $Me_{2}bik$ (2.5 mmol) in water (5 mL) was added to a solution of $Fe(SO_{4})·7H_{2}O$ (1.2 mmol) in the same solvent (30 mL). The resulting deep blue solution was stirred for few

minutes. Then a freshly prepared solution of KSCN/KScCN (16 mmol) in water (5 mL) was quickly added to the blue reaction mixture. A green precipitate immediately appeared. The turbid mixture was allowed to stir for 30 min. The precipitate was isolated by water suction, washed vigorously with distilled water and dried in vacuum. The crude product was dissolved in a minimum amount of DMSO at 60°C and the resulting bluish gray solution (ca. 10 mL) was stirred for 30 min and filtered. Slow evaporation of the filtrate under ambient condition afforded in 2 weeks dark blue plate like crystals of 1 and 2 that were suitable for single-crystal X-ray diffraction. Yields: 81% (1), 86% (2). IR (ATR): (1) $v_{max}/cm^{-1}$: 3147, 3127, 2959, 2100 (NC_{NCS}), 2066 (NC_{NCS}), 1639 (CO_{bik}), 1520, 1484, 1417, 1287, 1169, 897 $cm^{-1}$. (2)!: $v_{max}/cm^{-1}$: 3144, 3112, 2953, 2097 (NC_{NCSe}), 2069 (NC_{NCSe}), 1628 (CO_{bik}), 1518, 1482, 1416, 1288, 1170, 892 $cm^{-1}$. Anal. Calcd for $C_{20}H_{20}FeN_{10}O_{2}Se_{2}$: C, 37.17; H, 3.12; N, 21.68 Found: C, 37.06; H, 3.27 N, 21.31. Anal. Calcd for $C_{20}H_{20}FeN_{10}O_{2}S_{2}$: C, 43.48; H, 3.65; N, 25.36 Found: C, 43.10; H, 3.70; N, 25.01.

## X-Ray Data Collection and Structure Refinement

A single crystal of 1 and 2 compound were selected, mounted onto a cryoloop, and transferred in a cold nitrogen gas stream. Intensity data were collected with a BRUKER Kappa-APEXII diffractometer with graphite-monochromated Mo-Kα radiation (λ = 0.71073 Å). Data collection were performed with APEX2 suite (BRUKER). Unit-cell parameters refinement, integration, and data reduction were carried out with SAINT program (BRUKER). SADABS (BRUKER) was used for scaling and multi-scan absorption corrections. In the WinGX suite of programs (Farrugia, 1999), the structure were solved with SHELXT-14 (Palatinus and Chapuis, 2007) program and refined by full-matrix least-squares methods using SHELXL-14 (Sheldrick, 2015). The Bilbao Crystallographic Server (de la Flor et al., 2016)
<figure><img src="image_12.png" /><figcaption>**FIGURE 10** | TDDFT computed absorption spectra of complexes low-spin state (Top) and high-spin state (Bottom) of 1 (red) and 2 (black) depicting important electronic transitions.</figcaption></figure>

was used to compare the crystal structure of 1 measured at 200 K to others because it was recorded in a different setting.

CCDC 1844949-1844955 contain the supplementary crystallographic data for this paper. These data can be obtained free of charge from The Cambridge Crystallographic Data Centre via www.ccdc.cam.ac.uk/data_request/cif.

## Analytical and Spectroscopic Measurements

Elemental analyses for C, H, and N were performed on a Perkin-Elmer 240C analyser at the ISCN (Gif sur Yvette, France).

FTIR spectroscopic data were carried out on a Vertex 70 Bruker instrument working in the ATR mode and collected in the 400–4000 cm⁻¹ range at room temperature (with a 4 cm⁻¹ resolution).

Solid-state UV-vis spectra were measured in the range of 350–1200 nm on a CARY 5000 double-beam spectrophotometer equipped with the Eurolabo variable-temperature cell (21525, quartz windows) and Specac temperature controller. The

<table><caption>**TABLE 8** | Computed absorption spectra for complexes 1 and 2 along the various contributions to these transitions and their assignments (see Supplementary Figure S6 in ESI for orbital diagram).</caption><thead><tr><th></th><th>λ (nm)</th><th>E (cm<sup>-1</sup>)</th><th>Osc. Strength</th><th>Nature of the transition <sup>a</sup></th></tr></thead><tbody><tr><td rowspan="3">High-spin</td><td>298</td><td>33616</td><td>0.374</td><td>π → π<sup>*</sup> (N<sup>bik</sup> → N<sup>bik</sup>)</td></tr><tr><td>322</td><td>31089</td><td>0.031</td><td>LMCT (NCS → Fe)</td></tr><tr><td>402</td><td>24897</td><td>0.000</td><td>MLCT (Fe → N<sup>bik</sup>)</td></tr><tr><td rowspan="2">[FeII(Me bik)₂(NCS)₂]</td><td>488</td><td>20495</td><td>6.83E-05</td><td>LL'CT (NCS → N<sup>bik</sup>)</td></tr><tr><td>904</td><td>11062</td><td>0.029</td><td>d-d (Fe d<sub>xz</sub> → Fe d<sub>yz</sub>)</td></tr><tr><td rowspan="5">FeII(Me bik)₂(NCSe)₂</td><td>320</td><td>31256</td><td>0.029</td><td>MLCT (Fe → N<sup>bik</sup>)</td></tr><tr><td>340</td><td>29407</td><td>0.054</td><td>LMCT (NCS → Fe)</td></tr><tr><td>435</td><td>22981</td><td>0.005</td><td>LL'CT (NCS → N<sup>bik</sup>)</td></tr><tr><td>588</td><td>16997</td><td>0.006</td><td>LL'CT (NCS → N<sup>bik</sup>)</td></tr><tr><td>1023</td><td>9777</td><td>0.084</td><td>d-d (Fe → Fe)</td></tr><tr><td rowspan="5">Low-spin</td><td>329</td><td>30383</td><td>0.081</td><td>π → π<sup>*</sup> (N<sup>bik</sup> → N<sup>bik</sup>)</td></tr><tr><td>354</td><td>28273</td><td>0.026</td><td>MLCT (Fe → N<sup>bik</sup>)</td></tr><tr><td>448</td><td>22301</td><td>0.011</td><td>LL'CT (NCS → N<sup>bik</sup>)</td></tr><tr><td>506</td><td>19765</td><td>0.005</td><td>LL'CT (NCS → N<sup>bik</sup>)</td></tr><tr><td>611</td><td>16377</td><td>0.045</td><td>MLCT (Fe → N<sup>bik</sup>)</td></tr><tr><td rowspan="5">FeII(Me bik)₂(NCSe)₂</td><td>351</td><td>28457</td><td>0.023</td><td>LL'CT (NCS → N<sup>bik</sup>)</td></tr><tr><td>429</td><td>23303</td><td>0.005</td><td>MLCT (Fe → N<sup>bik</sup>)</td></tr><tr><td>478</td><td>20931</td><td>0.015</td><td>LL'CT (NCS → N<sup>bik</sup>)</td></tr><tr><td>596</td><td>16769</td><td>0.074</td><td>MLCT (Fe → N<sup>bik</sup>)</td></tr></tbody></table>

measurements were performed on KBr pellets. ~2 mg of fresh crystallites of compounds were dispersed without any grinding in *ca.* 99 mg of KBr, this latter being previously ground. This sample preparation aimed at minimizing the formation of crystalline defects that could alter the SCO characteristics.

## Magnetic Measurements

Magnetic susceptibility data were collected using a Quantum Design SQUID magnetometer (MPMS-5S Model) calibrated against a standard palladium sample. The magnetic susceptibility values were corrected from the diamagnetism of the molecular constituents and of the sample holder. The measurement was carried out on fresh crystals of **1** and **2** in the temperature range 4–400 K. Magnetic field of 1T was used for the study.

Photomagnetic measurements were carried out by using a sample holder equipped with an optical fiber. In a typical experiment 0.5 mg of finely ground crystals were deposited on an adhesive tape. The sample was separated from the end of the fiber by 5.5 cm. All the irradiations were carried out at 20 K to minimize the temperature variation induced by light. Both the samples were irradiated with laser diodes of wavelength 405, 532, 635, 808, 900, and 1313 nm (ca. 5–10 mW cm⁻²). The experimental data were corrected in a similar way to that of bulk magnetic susceptibility measurement.

## Theoretical Details

DFT calculations were performed on the X-ray structural coordinates using the Gaussian 09 (Frisch et al., 2009) suite programs for complexes **1-2** in 200, 300, and 400 K structures.

<footnote><sup>a</sup>MLCT, LMCT and LL'CT denoting transitions with a main character of metal-to-ligand, ligand-to-metal and ligand-to-ligand charge transfer.</footnote>
In order to investigate the dependence of the spin state energetic on the amount of exact exchange in the B3LYP functional, the %HF exchange is varied systematically from 0 to 25% also by adjusting the DFT exchange accordingly (Becke, 1988) and geometry optimization followed by the frequency calculation is performed in the DMSO solvent environment using PCM solvent model (Tomasi et al., 2005). Additional frequency calculations are also performed on top of the optimized coordinates by incorporating the dispersion correction (D2) recommended by Grimme (2006) Additional frequency calculations were performed to extract the dispersion corrected thermodynamic parameters. All calculations employ Ahlrichs polarized triple-ζ valence (TZVP) basis set for iron, nitrogen, sulfur, selenium and single-ζ valence (SVP) basis set for rest of the atoms (Schäfer et al., 1992, 1994). The geometry optimization and frequency calculations were also performed using OPBE(Perdew et al., 1997; Handy and Cohen, 2001) functional as this is shown to be superior in estimating thermodynamic quantities (Kepp, 2013). Beside, the time dependent DFT (TDDFT) calculations were performed on the low spin and high spin structures of complexes **1** and **2** using ORCA 3.0.3 suite (Neese, 2012) by employing B3LYP (Becke, 1988) as functional along with def2-TZVP (Schäfer et al., 1992, 1994; Weigend and Ahlrichs, 2005) basis set for iron, sulfur and selenium, and TZVP (Schäfer et al., 1992, 1994; Weigend and Ahlrichs, 2005) basis set for the rest of the atoms also by incorporating the solvent effects using COSMO (Tomasi et al., 2005; Ye and Neese, 2010b) method. TDDFT calculations were performed

in order to understand the absorption properties of these complexes.

## AUTHOR CONTRIBUTIONS

SD and HS prepared the samples and made the measurements. AF provided some help in the NMR measurements. YL has made some of the magnetic measurements. L-MC carried out the XRD analyses. RL supervised the work and wrote the draft. ST carried out the theoretical calculation and participate to the writing process. GR supervised the theoretical study and contributed to writing the article. M-LB supervised part of the work and the UV-vis measurement, and actively contributed to the improvement of the article.

## ACKNOWLEDGMENTS

The authors thank the Centre National de la Recherche Scientifique (CNRS, France) for its financial support and the Ministère de l'Enseignement Supérieur et de la Recherche (France). GR would like to thank SERB (EMR/2014/000247) for financial support.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2018.00326/full#supplementary-material

## REFERENCES

* Becke, A. D. (1988). Density-functional exchange-energy approximation with correct asymptotic behavior. *Phys. Rev.* A 38, 3098-3100. doi: 10.1103/PhysRevA.38.3098

* Bowman, D. N., and Jakubikova, E. (2012). Low-spin versus high-spin ground state in pseudo-octahedral iron complexes. *Inorg. Chem.* 51, 6011-6019. doi: 10.1021/ic202344w

* Buron-Le Cointe, M., Hébert, J., Baldé, C., Moisan, N., Toupet, L., and Guionneau, P. (2012). Intermolecular control of thermoswitching and photoswitching phenomena in two spin-crossover polymorphs. *Phys. Rev. B* 85, 064114-064122. doi: org/10.1103/PhysRevB.85.064114

* Casanova, D., Llunell, M., Alemany, P., and Alvarez, S. (2005). The Rich stereochemistry of eight-vertex polyhedra: a continuous shape measures study. *Chem. Eur. J.* 11, 1479-1494. doi: 10.1002/chem.200400799

* de la Flor, G., Orobengoa, D., Tasci, E., Perez-Mato, J. M., and Aroyo, M. I. (2016). Comparison of structures applying the tools available at the Bilbao crystallographic server. *J. Appl. Crystallogr.* 49, 653-664. doi: 10.1107/S1600576716002569

* De, S., Tewary, S., Garnier, D., Li, Y., Gontard, G., Lisnard, L., et al. (2018). Solution and solid-state study of the spin-crossover $[Fe^{II}(R-bik)_{3}](BF_{4})_{2}$ complexes ($R$ = Me, Et, Vinyl): solution and solid-state study of the spin-crossover $[Fe^{II}(R-bik)_{3}](BF_{4})_{2}$ complexes ($R$ = Me, Et, Vinyl). *Eur. J. Inorg. Chem.* 2018, 414-428. doi: 10.1002/ejic.201701013

* Farrugia, L. J. (1999). *WinGX* suite for small-molecule single-crystal crystallography. *J. Appl. Crystallogr.* 32, 837-838. doi: 10.1107/S0021889899006020

* Feringa, B. L. (2011). *Molecular Switches*. Weinheim: Wiley-VCH-Verl.

* Ferrando-Soria, J., Vallejo, J., Castellano, M., Martínez-Lillo, J., Pardo, E., Cano, J., et al. (2017). Molecular magnetism, quo vadis? A historical perspective from a coordination chemist viewpoint*. *Coord. Chem. Rev.* 339, 17-103. doi: 10.1016/j.ccr.2017.03.004

* Frisch, M.-J., Trucks, G.-W., Schlegel, H. B., Scuseria, G.-E., Robb, M.-A., Cheeseman, J.-R., et al. (2009). *Gaussian 09, Revision A.02*. Wallington, CT: Gaussian Inc.

* Glijer, D., Hébert, J., Trzop, E., Collet, E., Toupet, L., Cailleau, H., et al. (2008). Photoinduced phenomena and structural analysis associated with the spin-state switching in the $[Fe^{II}(DPEA)(NCS)_{2}]$ complex. *Phys. Rev. B* 78, 134112-134120. doi: 10.1103/PhysRevB.78.134112

* Grimme, S. (2006). Semiempirical GGA-type density functional constructed with a long-range dispersion correction. *J. Comput. Chem.* 27, 1787-1799. doi: 10.1002/jcc.20495

* Gütlich, P., Ksenofontov, V., and Gaspar, A. B. (2005). Pressure effect studies on spin crossover systems. *Coord. Chem. Rev.* 249, 1811-1829. doi: 10.1016/j.ccr.2005.01.022

* Halcrow, M. A. (2013). *Spin-Crossover Materials: Properties and Applications*. Chichester: Wiley.

* Handy, N. C., and Cohen, A. J. (2001). Left-right correlation energy. *Mol. Phys.* 99, 403-412. doi: 10.1080/00268970010018431

* Jamula, L. L., Brown, A. M., Guo, D., and McCusker, J. K. (2014). Synthesis and characterization of a high-symmetry ferrous polypyridyl complex: approaching the $^{5}$T$_{2}$/$^{3}$T$_{1}$ crossing point for $Fe^{II}$. *Inorg. Chem.* 53, 15-17. doi: 10.1021/ic402407k

* Kepp, K. P. (2013). Consistent descriptions of metal-ligand bonds and spin-crossover in inorganic chemistry. *Coord. Chem. Rev.* 257, 196-209. doi: 10.1016/j.ccr.2012.04.020

* Klingele, J., Kaase, D., Klingele, M. H., and Lach, J. (2012). Synthesis and mononuclear complexes of the bis-bidentate ligand 2,5-di(2-pyridyl)-1,3,4-thiadiazole (dptd): spin crossover in $[Fe(II)(dptd)_{2}(NCS)_{2}]$ and $[Fe(II)(dptd)_{2}(NCBH_{3})_{2}] \cdot H_{2}O$. *Dalton Trans.* 41, 1397-1406. doi: 10.1039/c1dt11396h

* Konno, M., and Mikami-Kido, M. (1991). Temperature- or pressure-induced structure changes of a spin crossover $Fe(II)$ complex; $[Fe(bpy)_{2}(NCS)_{2}]$. *Bull. Chem. Soc. Jpn.* 64, 339-345. doi: 10.1246/bcsj.64.339
Kumar, K. S., Salitros, I., Heinrich, B., Fuhr, O., and Ruben, M. A. (2015). A charge neutral iron (II) complex with an above room temperature spin crossover (SCO) and hysteresis loop. *Mater. Chem.* 3, 11635–11644. doi: 10.1039/C5TC02079D

Létard, J.-F. (2006). Photomagnetism of iron(ii) spin crossover complexes—the T (LIESST) approach. *J. Mater. Chem.* 16, 2550–2559. doi: 10.1039/B603473J

Llunell, M., Casanova, D., Cirera, J., Bofill, J. M., Alemany, P., Alvarez, S. et al. (2005). *SHAPE*, v2.1. University of Barcelona and The Hebrew University of Jerusalem: Barcelona.

Lucas, P., El Mehdi, N., Ho, H. A., Bélanger, D., and Breau, L. (2000). Expedient synthesis of symmetric aryl ketones and of ambient-temperature molten salts of imidazole. *Synthesis* 2000, 1253–1258. doi: 10.1055/s-2000-6416

MacLean, E. J., McGrath, C. M., O'Connor, C. J., Sangregorio, C., Seddon, J. M. W., Sinn, E., et al. (2003). Structural study of the thermal and photochemical spin states in the spin crossover complex $[Fe(phen)_{2}(NCSe)_{2}]$. *Chem. Eur. J.* 9, 5314–5322. doi: 10.1002/chem.200304865

Marino, A., Chakraborty, P., Servol, M., Lorenc, M., Collet, E., and Hauser, A. (2014). The role of ligand-field states in the ultrafast photophysical cycle of the prototypical iron(II) spin-crossover complex $[Fe(pzt)_{6}](BF_{4})_{2}$. *Angew. Chem. Int. Ed.* 53, 3863–3867. doi: 10.1002/anie.201310884

Martínez, V., Gaspar, A. B., Muñoz, M. C., Bukin, G. V., Levchenko, G., and Real, J. A. (2009). Synthesis and characterisation of a new series of bistable iron (II) spin-crossover 2D metal-organic frameworks. *Chem. Eur. J.* 15, 10960–10971. doi: 10.1002/chem.200901391

Mengel, A. K., Förster, C., Breivogel, A., Mack, K., Ochsmann, J. R., Laquai, F., et al. (2015). A heteroleptic push-pull substituted iron(II) Bis(tridentate) complex with low-energy charge-transfer states. *Chem. Eur. J.* 21, 704–714. doi: 10.1002/chem.201404955

Mondal, A., Li, Y., Chamoreau, L.-M., Seuleiman, M., Rechignat, L., Bousseksou, A., et al. (2014). Photo- and thermo-induced spin crossover in a cyanide-bridged $[Mo^{2}Fe^{II}]$ rhombus molecule. *Chem. Commun.* 50, 2893–2895. doi: 10.1039/C3CC49164A

Mondal, A., Li, Y., Herson, P., Seuleiman, M., Boillot, M.-L., Rivière, E., et al. (2012). Photomagnetic effect in a cyanide-bridged mixed-valence $[Fe^{II}Fe^{III2}]$ molecular square. *Chem. Commun.* 48:5653. doi: 10.1039/c2cc17835d

Mondal, A., Li, Y., Seuleiman, M., Julve, M., Toupet, L., Buron-Le Cointe, M., et al. (2013). On/Off photoswitching in a cyanide-bridged $[Fe_{2}Co_{2}]$ magnetic molecular square. *J. Am. Chem. Soc.* 135, 1653–1656. doi: 10.1021/ja3087467

Moussa, N. O., Trzop, E., Mouris, S., Zein, S., Molnár, G., Gaspar, A. B., et al. (2007). Wavelength selective light-induced magnetic effects in the binuclear spin crossover compound $[[Fe(bt)(NCS)_{2}]_{2}(bpym)]$. *Phys. Rev. B* 75, 1–8. doi: 10.1103/PhysRevB.75.054101

Nakamoto, T., Tan, Z. C., and Sorai, M. (2001). Heat capacity of the spin crossover complex $[Fe(2-pic)_{3}]Cl_{2}·MeOH$: a spin crossover phenomenon with weak cooperativity in the solid state. *Inorg. Chem.* 40, 3805–3809. doi: 10.1021/ic010073z

Nakano, K., Suemura, N., Kawata, S., Fuyuhiro, A., Yagi, T., Nasu, S., et al. (2004). Magnetic behavior and Mossbauer spectra of spin-crossover pyrazolate bridged dinuclear diiron(II) complexes: x-ray structures of high-spin and low-spin $[(Fe(CNBH_{3})(py))_{2}(\mu-bpyz)_{2}]$. *Dalton Trans.* 0, 982–988. doi: 10.1039/b314211f

Neese, F. (2012). The ORCA program system: the ORCA program system. *Wiley Interdiscip. Rev. Comput. Mol. Sci.* 2, 73–78. doi: 10.1002/wcms.81

Palatinus, L., and Chapuis, G. (2007). SUPERFLIP – a computer program for the solution of crystal structures by charge flipping in arbitrary dimensions. *J. Appl. Crystallogr.* 40, 786–790. doi: 10.1107/S0021889807029238

Pap, J. S., Bányai, V., Szilvási, D. S., Kaizer, J., Speier, G., and Giorgi, M. (2011a). Influence of meridional N3-ligands on supramolecular assembling and redox behavior of carboxylatocopper(II) complexes. *Inorg. Chem. Commun.* 14, 1767–1772. doi: org/10.1016/j.inoche.2011.08.005

Pap, J. S., Kripli, B., Bányai, V., Giorgi, M., Korecz, L., Gajda, T., et al. (2011b). Tetra-, penta- and hexacoordinate Copper(II) complexes with N3 donor isoindoline-based ligands: characterization and SOD-like activity. *Inorg. Chim. Acta* 376, 158–169. doi: 10.1016/j.ica.2011.06.001

Pap, J. S., Kripli, B., Giorgi, M., Kaizer, J., and Speier, G. (2011c). Redox properties of cobalt(II) complexes with isoindoline-based ligands. *Trans. Met. Chem.* 36, 481–487. doi: 10.1007/s11243-011-9493-z

Paulsen, H. (2016). Periodic density functional calculations in order to assess the cooperativity of the spin transition in

$fe(phen)2(NCS)2$. *Magnetochemistry* 2:14. doi: 10.3390/magnetochemistry2010014

Perdew, J. P., Burke, K., and Ernzerhof, M. (1997). Generalized gradient approximation made simple. *Phys. Rev. Lett.* 78, 1396–1396. doi: 10.1103/PhysRevLett.78.1396

Rackwitz, S., Klopper, W., Schünemann, V., and Wolny, J. A. (2013). Quantification of intramolecular cooperativity in polynuclear spin crossover Fe(ii) complexes by density functional theory calculations. *Phys. Chem. Chem. Phys.* 15:15450. doi: 10.1039/c3cp52240g

Real, J. A., Gallois, B., Granier, T., Suez-Panama, F., and Zarembowitch, J. (1992). Comparative investigation of the spin-crossover compounds $Fe(btz)_{2}(NCS)_{2}$ and $Fe(phen)_{2}(NCS)_{2}$ (where btz = 2,2'-bi-4,5-dihydrothiazine and phen = 1,10-phenanthroline). Magnetic properties and thermal dilatation behavior and crystal structure of $Fe(btz)_{2}(NCS)_{2}$ at 293 and 130 K. *Inorg. Chem.* 31, 4972–4979. doi: 10.1021/ic00050a013

Real, J. A., Gaspar, A. B., Niel, V., and Muñoz, M. C. (2003). Communication between Iron (II) building blocks in cooperative spin transition phenomena. *Coord. Chem. Rev.* 236, 121–141. doi: 10.1016/S0010-8545(02)00220-5

Ross, T. M., Moubaraki, B., Wallwork, K. S., Batten, S. R., and Murray, K. S. (2011). A temperature-dependent order-disorder and crystallographic phase transition in a 0D Fe(II) spin crossover compound and its non-spin crossover Co(II) isomorph. *Dalton Trans.* 40, 10147–10155. doi: 10.1039/c1dt10961h

Sauvage, J.-P., and Amendola, V. (2001). *Molecular Machines and Motors*. Berlin: Springer.

Schäfer, A., Horn, H., and Ahlrichs, R. (1992). Fully optimized contracted Gaussian basis sets for atoms Li to Kr. *J. Chem. Phys.* 97, 2571–2577. doi: 10.1063/1.463096

Schäfer, A., Huber, C., and Ahlrichs, R. (1994). Fully optimized contracted Gaussian basis sets of triple zeta valence quality for atoms Li to Kr. *J. Chem. Phys.* 100, 5829–5835. doi: 10.1063/1.467146

Scheja, A., Baabe, D., Menzel, D., Pietzonka, C., Schweyen, P., and Bröring, M. (2015). Spin crossover and valence tautomerism in neutral homoleptic iron complexes of Bis(pyridylimino)isoindolines. *Chem. Eur. J.* 21, 14196–14204. doi: 10.1002/chem.201501551

Scott, H. S., Ross, T. M., Chilton, N. F., Gass, I. A., Moubaraki, B., Chastanet, G., et al. (2013). Crown-linked dipyridylamino-triazine ligands and their spin-crossover Iron (II) derivatives: magnetism, photomagnetism and cooperativity. *Dalton Trans.* 42, 16494–16509. doi: 10.1039/c3dt51839f

Sheldrick, G. M. (2015). Crystal structure refinement with SHELXL. *Acta Crystallogr. Sect. C Struct. Chem.* 71, 3–8. doi: 10.1107/S2053229614024218

Shepherd, H. J., Gural'skiy, I. A., Quintero, C. M., Tricard, S., Salmon, L., Molnár, G., et al. (2013). Molecular actuators driven by cooperative spin-state switching. *Nat. Commun.* 4:2607. doi: 10.1038/ncomms3607

Shepherd, H. J., Palamarciu, T., Rosa, P., Guionneau, P., Molnár, G., Létard, J.-F., et al. (2012). Antagonism between extreme negative linear compression and spin crossover in $[Fe(dpp)_{2}(NCS)_{2}]_{2}·py$. *Angew. Chem. Int. Ed.* 51, 3910–3914. doi: 10.1002/anie.201108919

Slichter, C. P., and Drickamer, H. G. (1972). Pressure-induced electronic changes in compounds of iron. *J. Chem. Phys.* 56, 2142–2160. doi: 10.1063/1.1677511

Tomasi, J., Mennucci, B., and Cammi, R. (2005). Quantum mechanical continuum solvation models. *Chem. Rev.* 105, 2999–3094. doi: 10.1021/cr9904009

Weigend, F., and Ahlrichs, R. (2005). Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: design and assessment of accuracy. *Phys. Chem. Chem. Phys.* 7, 3297–3305. doi: 10.1039/b508541a

Ye, S., and Neese, F. (2010b). Accurate modeling of spin-state energetics in spin-crossover systems with modern density functional theory. *Inorg. Chem.* 49, 772–774. doi: 10.1021/ic902365a

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2018 De, Chamoreau, El Said, Li, Flambard, Boillot, Tewary, Rajaraman and Lescouëzec. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.