The Highly Conducting Spin-Crossover Compound Combining Fe ( III ) Cation Complex with TCNQ in a Fractional Reduction State . Synthesis , Structure , Electric and Magnetic Properties

Three systems [Fe(III)(sal2-trien)](TCNQ)n·X (n = 1, 2, X = MeOH, CH3CN, H2O) showing spin-crossover transition, conductivity and ferromagnetic coupling were synthesized and studied by X-ray diffraction, Montgomery method for resistivity, SQUID magnetometry and X-band EPR. Spin-spin interactions between local magnetic moments of Fe(III) ions and electron spins of organic TCNQ network were discovered and discussed within the framework of intermolecular superexchange coupling.


Introduction
The synthesis and investigation of multifunctional molecule materials combining conductivity and magnetism in the same crystal lattice attract considerable attention, because the interplay of these properties may lead to novel behavior [1][2][3][4][5][6].Until recently, research in this direction has been focused on the family of the quasi-two-dimensional (super)conductors based on the radical cation salts of bis(ethylenedithio)tetrathiafulvalene (BEDT-TTF) and its derivatives with paramagnetic metal complex anions of different nature [2][3][4][5][6].In such materials, conductivity is associated with mobile electrons in organic layers, whereas magnetism usually originates from localized spins of transition metal ions in insulating counterion layers.In particular, salts of BEDT-TTF and its selenium-substituted derivative bis(ethylenedithio)tetraselenafulvalene (BETS) have been shown to combine (super)conducting and paramagnetic [7][8][9][10] and even antiferromagnetic [11,12] and ferromagnetic properties [13].Moreover, interaction between localized spins in insulating magnetic layers and itinerant spins in conducting organic layers was found to lead to new fascinating phenomena such as field-induced superconductivity observed on λ-(BETS) 2 FeCl 4 [14] and κ-(BETS) 2 FeBr 4 [15].

Synthesis
The compound [Fe(III)(sal 2 -trien)](TCNQ) 2 •CH 3 OH (1•MeOH) were obtained by mixing CH 3 CN solution of TCNQ with the solutions of [Fe(III)(sal 2 -trien)]I•1.5H 2 O and LiTCNQ in a mixture of acetonitrile/methanol (see Materials and Methods).The resulting solution was left to stand in a refrigerator overnight.The black plate crystals of 1•MeOH were formed (Figure S1). Figure S2 displays the thermogram of 1•MeOH.With increasing temperature, a weight loss of 4.43% was observed in the temperature range 70-150 • C with endothermic peak at 126.7 • C, which is assigned to the loss of crystallization methanol molecule (calc.3.77%).As this takes place, the ions with m/e 29 (CHO) and 31 (CH 3 O) relating to the fragments of methanol molecule are observed in the mass spectrum.On heating above 200 • C, the complex begins to decompose (two DSC-peaks at 210.1 and 222.6 • C).
The compound [Fe(III)(sal 2 -trien)](TCNQ)•H 2 O (3) in the polycrystalline form was prepared by metathesis reaction between [Fe(III)(sal 2 -trien)]I•1.5H 2 O and Li(TCNQ) in methanol.The recrystallization of 3 from acetonitrile yielded the crystals of [Fe(III)(sal 2 -trien)](TCNQ)•CH 3 CN (2) (Materials and Methods).The thermogram of 2 demonstrates a mass loss of 6.38% in the temperature range 70-100 • C with endothermic peak at 93 • C, which corresponds to the loss of lattice acetonitrile molecule (calc.6.27%).(Figure S3).In the mass spectrum recorded in the gas phase, the peaks are observed at m/z 41, 15 and 26 relating to CH 3   The structure of 1•MeOH has been investigated at several different temperatures between 100 and 385 K.The asymmetric unit of the triclinic P1 unit cell of initial 1•MeOH crystal includes two independent TCNQ •δ− radical anions (denoted as I and II), one [Fe(sal 2 -trien)] + cation and one methanol molecule in general positions (Figure 1).The crystal structure consists of layers of TCNQ radical anions parallel to the ab-plane which alternate with mixed Fe(sal 2 -trien)-solvent layers along the c-axis (Figure 2).The structure of 1•MeOH has been investigated at several different temperatures between 100 and 385 K.The asymmetric unit of the triclinic P1 unit cell of initial 1•MeOH crystal includes two independent TCNQ •δ− radical anions (denoted as I and II), one [Fe(sal2-trien)] + cation and one methanol molecule in general positions (Figure 1).The crystal structure consists of layers of TCNQ radical anions parallel to the ab-plane which alternate with mixed Fe(sal2-trien)-solvent layers along the c-axis (Figure 2).In the 100 K structure, charge state of both independent TCNQ •δ− radical anions I and II is slightly different and close to −0.4 and −0.6, respectively, according to calculations on the base of intramolecular bond lengths using Kistenmacher's empirical formula [41] (Table S1).π-stacking of TCNQ anions leads to a formation of infinite stacks along the [1 1 0] direction.The stacks are composed of tetrads …-(II-I-I-II)-(II-I-I-II)-… (Figure 2) with ring-over-bond overlap inside the tetrad due to longitudinal shift of adjacent anions and ring-over-ring one between the tetrads   In the 100 K structure, charge state of both independent TCNQ •δ− radical anions I and II is slightly different and close to −0.4 and −0.6, respectively, according to calculations on the base of intramolecular bond lengths using Kistenmacher's empirical formula [41] (Table S1).π-stacking of TCNQ anions leads to a formation of infinite stacks along the [1 1 0] direction.The stacks are composed of tetrads …-(II-I-I-II)-(II-I-I-II)-… (Figure 2) with ring-over-bond overlap inside the tetrad due to longitudinal shift of adjacent anions and ring-over-ring one between the tetrads In the 100 K structure, charge state of both independent TCNQ •δ− radical anions I and II is slightly different and close to −0.4 and −0.6, respectively, according to calculations on the base of intramolecular bond lengths using Kistenmacher's empirical formula [41] (Table S1).π-stacking of TCNQ anions leads to a formation of infinite stacks along the  S2) that is also characteristic for the LS state.The coordination FeN 4 O 2 octahedrons are much more regular in LS state than in HS one [24][25][26]34].There is no π-stacking between phenolate ligands of the nearest cations.Most of the cation-anion hydrogen contacts are formed with TCNQ II that is in agreement with higher charge of the latter in comparison with TCNQ I. (longitudinal shift is absent, Figure S4).Interplane separations in the stack at 100 K are 3.274(4) Å (I…I), 3.19(4) Å (I…II) and 3.297( 6) Å (II…II), the angle between mean planes of I and II is 2.37(5)°.
There are a lot of intrastack C…C contacts shorter than 3.4 Å, while adjacent stacks interact only by weak hydrogen contacts of С-Н…N type with H…N distances of 2.58-2.73Å.
The [Fe(III)(sal2-trien)] + cations are located between the TCNQ layers and connected by hydrogen contacts with surrounding TCNQ and methanol molecules.The Fe(III) ion has a slightly distorted octahedral coordination geometry formed by two oxygen and four nitrogen atoms of ligand.Key bond distances and angles in the Fe(III)(sal2-trien) + cation at 100 K are listed in Table S2.Mean values of the Fe-O, Fe-Nimine and Fe-Namine bond distances in the FeN4O2 chromophore (1.885(8), 1.936(8) and 2.007(2) Å, respectively) and shape of the cation with the dihedral angle α between two phenoxy groups of 53.87(6)° (Figure 3a) correspond to the low-spin state of the Fe(III) ion.The diagonals of the octahedron at Fe centre are close to linear (N-Fe-O and N-Fe-N angles are 174°-177°, Table S2) that is also characteristic for the LS state.The coordination FeN4O2 octahedrons are much more regular in LS state than in HS one [24][25][26]34].There is no π-stacking between phenolate ligands of the nearest cations.Most of the cation-anion hydrogen contacts are formed with TCNQ II that is in agreement with higher charge of the latter in comparison with TCNQ I.There are two N-Hcat…NII contacts with H…N distances of 2.  In order to follow the spin-crossover transition in the crystal of 1•MeOH, additional X-ray diffraction experiments were performed at several temperatures between 100 and 385 K.The X-ray study showed single crystal to single crystal phase transition with solvent loss at annealing the crystal at 350 K in nitrogen stream.According to crystal structure refinement, the methanol site is fully occupied at 295 K [refining the occupancy gives value of 0.988(4)].Process of solvent removal has started at the end of the 24-h X-ray experiment at 340 K [refined MeOH occupancy value is 0.956 (6)].It has been empirically found that annealing of the crystal at 350 K during one day transforms the crystal completely into a new phase with chemical formula [Fe(III)(sal2-trien)](TCNQ)2 (1), i.e., totally free of solvent.All further diffraction experiments were done for the new phase 1, overall sequence of temperature changes being following: warming to 350 and 385 K, further cooling to 295 and 220 K and again warming to 260 and 325 K with full data collection at all these temperatures.The treatment to carry out X-ray study at 400 K was unsuccessful due to melting of the crystal.Thus, temperatures of solvent loss and crystal decay in the X-ray experiments appeared to be about 50-80 K lower than the temperatures of corresponding peaks in DSC spectrum that is associated with the distinctly different conditions of heating in X-ray and thermogravimetric experiments.Symmetry of the crystal (triclinic P 1 ) persists upon the transition, but solvent loss causes significant distortion of the lattice (Materials and Methods).Relatively large temperature-dependent structural changes have been found in both cationic magnetic and anionic conducting layers.They are accumulated in Table 1 and Figure 4.All the structural parameters investigated show In order to follow the spin-crossover transition in the crystal of 1•MeOH, additional X-ray diffraction experiments were performed at several temperatures between 100 and 385 K.The X-ray study showed single crystal to single crystal phase transition with solvent loss at annealing the crystal at 350 K in nitrogen stream.According to crystal structure refinement, the methanol site is fully occupied at 295 K [refining the occupancy gives value of 0.988(4)].Process of solvent removal has started at the end of the 24-h X-ray experiment at 340 K [refined MeOH occupancy value is 0.956 (6)].It has been empirically found that annealing of the crystal at 350 K during one day transforms the crystal completely into a new phase with chemical formula [Fe(III)(sal 2 -trien)](TCNQ) 2 (1), i.e., totally free of solvent.All further diffraction experiments were done for the new phase 1, overall sequence of temperature changes being following: warming to 350 and 385 K, further cooling to 295 and 220 K and again warming to 260 and 325 K with full data collection at all these temperatures.The treatment to carry out X-ray study at 400 K was unsuccessful due to melting of the crystal.Thus, temperatures of solvent loss and crystal decay in the X-ray experiments appeared to be about 50-80 K lower than the temperatures of corresponding peaks in DSC spectrum that is associated with the distinctly different conditions of heating in X-ray and thermogravimetric experiments.Symmetry of the crystal (triclinic P1) persists upon the transition, but solvent loss causes significant distortion of the lattice (Materials and Methods).Relatively large temperature-dependent structural changes have been found in both cationic magnetic and anionic conducting layers.They are accumulated in Table 1 and Figure 4.All the structural parameters investigated show discontinuous changes between 340 and 350 K at the transition from solvate to solvent-free phase.In the magnetic subsystem, the Fe-N bond lengths and dihedral angle α between two phenoxy groups in the [Fe(III)(sal 2 -trien)] + cation, both of which are sensitive to spin state of the complex, show linear temperature dependences in the 100-340 K range and notably increase at 350 K (Figure 4a).An abrupt growth of these parameters reflects transition of Fe(III) from LS (S = 1/2) to HS state (S = 5/2).The mean Fe-N imine and Fe-N amine bond length values increase from 1.94 and 2.01 Å, typical for the LS Fe(III), to 2.00 and 2.08 Å, respectively, at 385 K. Taking into account corresponding average Fe-N distances in high-spin state, 2.12 and 2.20 Å [26,34], one can presume that about one third of iron centers are HS at this temperature.Stronger distortion of the FeN 4 O 2 octahedron at 385 K (diagonal angles are 164 • -177 • , Table S2) is an additional evidence of partial transition to HS state.It should be emphasized that in the temperature range 100-340 K, in which the solvate phase 1•MeOH exists, the Fe-N bonds correspond to LS state and do not show any notable growth.In the solvent-free phase 1 a reversible character of SCO transition is observed at further temperature cycling.The Fe-N bond lengths return to initial, typical for LS Fe(III) values on cooling to 220 K and go again to HS values on a new cycle of warming (empty marks in Figure 4a) but in the phase 1 these changes are smooth, without jump.The angle α in 1 at low temperatures is larger than it was in 1•MeOH (Figure 4a) reflecting appearance of an additional space in the crystal lattice of the solvent-free phase.
Table 1.The average Fe-N imine and Fe-N amine bond length values and dihedral angle α between two phenoxy groups in [Fe(III)(sal 2 -trien)] + cation in 1•MeOH and 1, the charges (δ) of different TCNQ •δ− species I and II estimated from Kistenmacher's empirical formula (Table S1) [41] and then scaled to make the total charge on TCNQ pair exactly equal to −1 and intrastack I-I, I-II and II-II interplanar separations at different temperatures.In contrast to the magnetic subsystem, changes in the conducting TCNQ layer occur in the entire temperature range investigated.Two independent TCNQ •δ− radical anions show tendency to a leveling of their charge state on warming (Figure 4b).Charges on TCNQ •δ− I and II radicals were calculated at different temperatures using the Kistenmacher formula (Table S1) and then scaled to make the total charge on the TCNQ pair exactly equal to −1 according to stoichiometry of the compound.The resulted scaled values are used in Table 1 and Figure 4b.Difference in charge on TCNQ •δ− I and II observed at 100 K fully disappears at 385 K and both of the anions become half-charged.Evolution of the charge is not uniform in 100-340 K range showing minimum for I and maximum for II at about room temperature.There is a jump at 350 K equalizing the charges.Further cooling leads to a linear divergence of charges to the starting values.Temperature evolution of the interplane separations along the TCNQ stack is not similar for three types of interactions: I-II and II-II distances grow on warming whereas I-I has a minimum at 220 K (Figure 4c).Again there are clear jumps on the heating curves in Figure 4c
In order to investigate temperature behavior of the structure, additional X-ray diffraction experiments were carried out on the crystal 2 at 220, 295 and 325 K.There was no indication of spin transition in the 100-325 K range because the mean Fe-N imine and Fe-N amine bond length values as well as dihedral angle α between two phenoxy groups in the [Fe(III)(sal 2 -trien)] + cation remain almost unchanged (Table 2).At 350 K the diffraction quality sharply decreases, apparently, due to solvent loss.Disappearance of the long range order in the crystal makes impossible further structural investigation of the spin-crossover transition.

Conductivity and Magnetic Properties
The normalized dc resistance, R(T)/R(300 K), was measured by a standard four-probe method for the single crystal [Fe(III)(sal 2 -trien)](TCNQ) 2 •CH 3 OH (1•MeOH) along the TCNQ stacks.It demonstrated a semiconducting type behavior (Figure 6).The placement of the electrodes on the crystal is shown in Figure S1.The value of the conductivity at room temperature is 1.5 Ω −1 •cm −1 .Below 110 K the data points were well described by the exponential law R theor (T) = exp(∆E/kT) with the energy gap ∆E = 0.05 eV (inset, Figure 6).
The conductivity at room temperature for crystals 1•MeOH was 1-2 order of magnitude higher, and the activation energy was an order less, than those reported for conducting SCO complexes of Fe(II) and Co(II) with TCNQ in a fractional reduction state [32,33,40].This difference was associated with the structure of the TCNQ stacks.In 1•MeOH, the stacks are rather regular by charge distribution and intermolecular separation, while in the other structures the stacks are subdivided into the pronounced TCNQ triads.It should be noted, that the complex [29] possesses a record conductivity among known conducting SCO complexes.In contrast to 1•MeOH, complex 2 was practically an insulator (σ 300K = 10 −8 Ω −1 •cm −1 ), that is characteristic of many TCNQ salts with a fully reduced TCNQ.The low conductivity is associated primarily with a strong dimerization of TCNQ radicals and absence of a long-range stacking.

Conductivity and Magnetic Properties
The normalized dc resistance, R(T)/R(300 K), was measured by a standard four-probe method for the single crystal [Fe(III)(sal2-trien)](TCNQ)2•CH3OH (1•MeOH) along the TCNQ stacks.It demonstrated a semiconducting type behavior (Figure 6).The placement of the electrodes on the crystal is shown in Figure S1.The value of the conductivity at room temperature is 1.5 Ω −1 •cm −1 .Below 110 K the data points were well described by the exponential law Rtheor(T) = exp(ΔE/kT) with the energy gap ΔE = 0.05 eV (inset, Figure 6).The bulk static magnetic susceptibility χ of 1•MeOH was measured on a polycrystalline sample in the temperature range of 2-400 K (Figure 7a).The entire range was scanned twice at rates 5 K/min in the interval 20-400 K and 0.5 K/min at the interval 2-20 K.The exposure time at the highest temperature 400 K was 10 min.Two cycles of measurements at external field values 0.1 T and 4.0 T provided coinciding data points above 20 K.The data of the latter cycle were omitted in Figure 7.The value of spin response below 100 K indicated that ~98% of Fe(III) ions exist in LS state.Spin response of conducting TCNQ sublattice appeared to be negligible within less than 10% (χ ≤ 2 × 10 −2 cm 3 •K•mol −1 ) [42,43].This is in agreement with earlier measured spin response of TCNQ layers in conducting Fe(II) SCO system [40].A temperature independent contribution (TIP) to the spin susceptibility was also found insignificant for the description of the low temperature evolution of χ(T) in Figure 7a.A total magnetization, M(B), measured at 2 K (Figure 7b), was described by a Brillouin curve for S = 1/2 (solid line).Thus, a paramagnetic response of 1•MeOH below 100 K is related to the local moments S = 1/2 of LS Fe(III).
Sharp growth of χT above ~150 K is consistent with the structural changes in TCNQ sublattices (Figure 4b,c) and deviation of the resistivity from the exponential law (Figure 6).It arises due to SCO transition and can be modeled by replacing S = 1/2 to S = 5/2 magnetic moments via Boltzmann activation mechanism [40].A respective curve is shown in Figure 7 (inset).The SCO parameters were T* = 410 K and ∆T 1/2 = 108 K for the mid point and the width of the transition at the levels ±1/2 of the midpoint.In contrast to the diffraction data, the behavior of the magnetic susceptibility was found completely reversible.This can be understood in terms of different experimental conditions (measurement time window: minutes vs. hours, atmosphere: helium vs. nitrogen, etc.).A fully realized HS state was not achieved, because it rests above the stability threshold for this structure as it follows from the thermogravimetric analysis (Figure S2).We concluded that solvate methanol did not leave the structure while heating-cooling cycles in the magnetic measurements.A calculated value, χT = 4.4 cm 3 •K•mol −1 , seemed a reliable estimate for total the spin response in HS state even though a complete was not achieved experimentally (Figure 7a, inset).A continuous transition curve allows determining mutual concentrations of S = 1/2 and S = 5/2 moments at any temperature within the SCO range.
provided coinciding data points above 20 K.The data of the latter cycle were omitted in Figure 7.The value of spin response below 100 K indicated that ~98% of Fe(III) ions exist in LS state.Spin response of conducting TCNQ sublattice appeared to be negligible within less than 10% (χ ≤ 2  10 −2 cm 3 •K•mol −1 ) [42,43].This is in agreement with earlier measured spin response of TCNQ layers in conducting Fe(II) SCO system [40].A temperature independent contribution (TIP) to the spin susceptibility was also found insignificant for the description of the low temperature evolution of χ(T) in Figure 7a.A total magnetization, M(B), measured at 2 K (Figure 7b), was described by a Brillouin curve for S = 1/2 (solid line).Thus, a paramagnetic response of 1•MeOH below 100 K is related to the local moments S = 1/2 of LS Fe(III).Sharp growth of χT above ~150 K is consistent with the structural changes in TCNQ sublattices (Figure 4b,c) and deviation of the resistivity from the exponential law (Figure 6).It arises due to SCO transition and can be modeled by replacing S = 1/2 to S = 5/2 magnetic moments via Boltzmann The EPR spectrum of 1•MeOH at 100 K was an intensive signal with anisotropic g-factor: g 1 = 2.219(9), g 2 = 2.171(0), g 3 = 1.961 (5).The evolution of the lineshape while heating up to 370 K is shown in Figure 8a.Temperature dependence of the g-components is presented in Figure 8b.The total spectrum was simulated at the temperatures of the measurement by using a standard lineshape model with anisotropic g-factor (WINEPR).The g-parameters were taken from the respective trend lines in Figure 8.An individual single line corresponding to the central g component has been extracted and further analyzed in terms of a relaxation rate.Qualitatively the linewidth behavior, ∆B(T), can be traced in Figure 8.Here we are showing the peak positions B p+ and B p− for low-and high-field halves of the first derivative signal (∆B = B p− − B p+ ).Therefore, we were able to avoid a g-strain effect and only analyze relaxation.
Magnetochemistry 2017, 3, 9 9 of 26 activation mechanism [40].A respective curve is shown in Figure 7 (inset).The SCO parameters were T* = 410 K and ΔT1/2 = 108 K for the mid point and the width of the transition at the levels ±1/2 of the midpoint.In contrast to the diffraction data, the behavior of the magnetic susceptibility was found completely reversible.This can be understood in terms of different experimental conditions (measurement time window: minutes vs. hours, atmosphere: helium vs. nitrogen, etc.).A fully realized HS state was not achieved, because it rests above the stability threshold for this structure as it follows from the thermogravimetric analysis (Figure S2).We concluded that solvate methanol did not leave the structure while heating-cooling cycles in the magnetic measurements.A calculated value, χT = 4.4 cm 3 •K•mol −1 , seemed a reliable estimate for total the spin response in HS state even though a complete was not achieved experimentally (Figure 7a, inset).A continuous transition curve allows determining mutual concentrations of S = 1/2 and S = 5/2 moments at any temperature within the SCO range.
The EPR spectrum of 1•MeOH at 100 K was an intensive signal with anisotropic g-factor: g1 = 2.219( 9), g2 = 2.171(0), g3 = 1.961 (5).The evolution of the lineshape while heating up to 370 K is shown in Figure 8a.Temperature dependence of the g-components is presented in Figure 8b.The total spectrum was simulated at the temperatures of the measurement by using a standard lineshape model with anisotropic g-factor (WINEPR).The g-parameters were taken from the respective trend lines in Figure 8.An individual single line corresponding to the central g component has been extracted and further analyzed in terms of a relaxation rate.Qualitatively the linewidth behavior, ΔB(T), can be traced in Figure 8.Here we are showing the peak positions Bp+ and Bp− for low-and high-field halves of the first derivative signal (ΔB = Bp− − Bp+).Therefore, we were able to avoid a gstrain effect and only analyze relaxation.The observed signal is a typical X-band spectrum of LS Fe(III) with S = 1/2 in six-coordinated position.When a geometry of the ligand changes (for example a dihedral angle α in Figure 4a), the symmetry of hexa-coordination would also change, from rhombic distortions to axial symmetry.Axial symmetry arises in HS state due to Jahn-Teller effect.In various 4N-coordinated chelate sites the g1 value reaches 3.0 following the ascending linear trend when coordination geometry gradually The observed signal is a typical X-band spectrum of LS Fe(III) with S = 1/2 in six-coordinated position.When a geometry of the ligand changes (for example a dihedral angle α in Figure 4a), the symmetry of hexa-coordination would also change, from rhombic distortions to axial symmetry.Axial symmetry arises in HS state due to Jahn-Teller effect.In various 4N-coordinated chelate sites the g 1 value reaches 3.0 following the ascending linear trend when coordination geometry gradually changes to axial symmetry [44].In other words, g 1 is a sensitive probe of ligand geometry as well as a solvent presence.On contrary to expected, g 1 and g 2 components in our experiment demonstrated weakly descending trend with the temperature.This is in an indicative agreement with unchanged F-N am and Fe-N im distances in Figure 4a.This also implies that the coordination geometry in LS complexes [Fe(III)(sal 2 -trien)] remains sustainable up to 370 K (63% LS, 37% HS).Variation of the dihedral angle between two phenoxy arms within ~5• does not cause a systematic effect on electronic structure of the metal cation.The HS configuration arises due to spontaneous thermal activation, without a precursor.The X-band spectrum for HS Fe(III), S = 5/2, usually shows a strong line at g = 6 (900 G < B < 1400 G) and a complimentary weak absorption at g ≈ 2.0 [45][46][47].This is because the S = 5/2 multiplet ( 6 A 1 ) forms three Kramers' doublets that are separated by energies significantly larger than the microwave quantum, ~0.3 cm −1 for X-band.For higher rhombicities (E/D) an additional line arises at g = 4.3 due to quantum-mechanically mixed states.We detected none of the three lines in the range 100-370 K. Since the contribution of HS Fe(III) magnetic moments in 1•MeOH was determined by SQUID measurements, this is an indication of fast spin relaxation.As we show further, this signal was observed in the insulating compounds 2 and 3. Hence, fast relaxation of S = 5/2 moments is caused by spin-spin interactions with mobile spins in TCNQ sublattice.
The concentration of LS [Fe(III)(sal 2 -trien)] complexes was independently verified by EPR.According to Schumacher-Slichter method a double integrated intensity of the total EPR spectrum, I EPR , is proportional to the respective spin concentration [48].This method works best for spin S = 1/2 and detection of relative changes.The absolute values determined at high temperatures and/or for broad spectra are less reliable.The product I EPR T, corresponding to the intensity I EPR (T)/I EPR (100 K) in Figure 7, has decreased 1.48 times by 370 K giving the estimate 67%.This is in good agreement with 63% portion of S = 1/2 moments obtained from SQUID data at the same temperature.
Local magnetic moments of LS Fe(III) can be considered as a useful probe for studying spin-spin interactions in the process of spin crossover.Indeed, the local environment remains stable as it follows from the behavior of g-tensor, whereas the linewidth experiences an abrupt broadening.The value of partial linewidth, ∆B, at 100 K was 18.8 G, that corresponds to a spin-lattice relaxation in an individual isolated complex LS [Fe(III)(sal 2 -trien)] at finite temperature.Figure 9 demonstrates a sharp growth of ∆B at increasing relative concentration of the moments S = 5/2, n/N(%).The value of n/N was extracted from the magnetic susceptibility data in Figure 7a (inset).This growth was fitted well by a simple expression ∆B = ∆B 0 + k(n/N), where ∆B 0 = −46 G is a base parameter relating to a threshold concentration (4.8%) and k = 21 G/% is a broadening factor.Weak and gradual broadening at the concentrations below 5% was reasonably affiliated with a trivial spin-phonon mechanism (spin-lattice relaxation).Above 5% the HS Fe(III) magnetic moments become a dominant source of relaxation.This is evidently a spin-spin relaxation.It looks credible that the relaxation rate in a reservoir A (S = 1/2) is proportional to the spin concentration in a reservoir B (S = 5/2).The question arises if this relaxation is of dipole-dipole nature, or a result of a weak superexchange coupling?The dipole-dipole could be a favorable mechanism once the dramatic broadening had been confirmed for other [Fe(III)(sal 2 -trien)] complexes.So far, we do not have the proofs.Exchange coupling might not seem a realistic scenario, unless the counterion system of TCNQ in a fractional reduction state, possessing highly mobile delocalized electrons.The abrupt broadening of the conduction electron EPR signal was observed in the SCO system [Fe(II){HC(pz) 3 } 2 ](TCNQ) 3 [40].Note, that LS Fe(II) cation does not possess a magnetic moment.The broadening in the SCO range was related to the spin concentration of HS Fe(II) ions with S = 2 [40].Owing to the fractional reduction state TCNQ triads, tetrads or even dimers might serve as efficient mediator of spin-spin interactions, when a direct overlapping of d-orbitals does not exist.
Absence of a direct exchange or superexchange does not eliminate a dynamic local spin density on ligands.In turn, ligand shell closely interacts with CN groups in TCNQ via a network of short contacts.This is a "bottleneck" junction, which facilitates cross relaxation and may affect the broadening factor k [49].One more argument that is favorable is the absence of EPR signal from TCNQ sublattice in 1•MeOH.A typical signal in TCNQ semiconductors has width in the range 1-10 G that makes it very detectable even at low concentrations.It was shown for the system [Fe(II){HC(pz) 3 } 2 ](TCNQ) 3 that at spin S = 2 concentrations above 5% in the SCO subsystem this signal rapidly spreads [40].If so, there is no question why we did not observe it in 1•MeOH where 100% of S = 1/2 moments were presented at all temperatures below the SCO.As soon as delocalized spins in anion sublattice of 1•MeOH were not detected by EPR, due to fast cross relaxation, and by SQUID, due to its negligible contribution to the total susceptibility, the only source of information about their reaction to SCO was an electron scattering.Below 110 K the resistance, R(T)/R(300 K), perfectly obeys thermal activation mechanism (Figure 6).One would expect that this mechanism also works at the higher temperatures, unless the structure of TCNQ stacks change substantially.Indeed, the gap of Δ = 0.05 eV corresponds to 600 K (kT), which is higher than T* = 410 K. Hence, theoretically the resistance was expected to follow the same trend in absence of SCO, Rtheor(T).However, the electron scattering can be affected by spin crossover due to local rearrangements inside the TCNQ tetrads I-II and II-II in the vicinity of HS complexes [Fe(III)(sal2trien)] (see Figure 4).For quantitative analysis we plotted a relative deviation of the actual resistance from its theoretical estimate, ΔR, against the concentration of the HS complexes, ΔR = (R(n/N) − Rtheor)/Rtheor, that was shown in Figure 9.The data were described by a logarithmic dependence ΔR = A•ln(B•n/N), where A = 0.3 and B = 3.0.Note, that the argument 3n/N = n/(N/3) indicates that changes of resistivity occur in a reduced scale, ~(N/3).It is also worth to note, that the same type dependence was found for the in-plane components of the resistivity, ρ||a and ρ||b, in the conducting SCO complex [Fe{HC(pz)3}2](TCNQ)3 [40].Moreover, it was shown in [40] that the transverse transport, ρ(T), did not react to the changes in the cation layers.Therefore, we suggested that in 1•MeOH all the observed changes of the resistance took place for the in-plane transport.Comparing data on Figure 9 one can conclude, that the major effect of SCO on transport properties occurs at lower concentrations n/N < 5%, whereas spin-spin depolarization is effective at n/N >> 5%.Similar effect was observed in [40].This is comprehensible and has a qualitative explanation.At the very low concentrations of HS [Fe(III)(sal2-trien)] complexes, the local displacements inside and/or between the TCNQ tetrads behave as randomly distributed defects.As the temperature approaches the midpoint of the As soon as delocalized spins in anion sublattice of 1•MeOH were not detected by EPR, due to fast cross relaxation, and by SQUID, due to its negligible contribution to the total susceptibility, the only source of information about their reaction to SCO was an electron scattering.Below 110 K the resistance, R(T)/R(300 K), perfectly obeys thermal activation mechanism (Figure 6).One would expect that this mechanism also works at the higher temperatures, unless the structure of TCNQ stacks change substantially.Indeed, the gap of ∆ = 0.05 eV corresponds to 600 K (kT), which is higher than T* = 410 K. Hence, theoretically the resistance was expected to follow the same trend in absence of SCO, R theor (T).However, the electron scattering can be affected by spin crossover due to local rearrangements inside the TCNQ tetrads I-II and II-II in the vicinity of HS complexes [Fe(III)(sal 2 -trien)] (see Figure 4).For quantitative analysis we plotted a relative deviation of the actual resistance from its theoretical estimate, ∆R, against the concentration of the HS complexes, ∆R = (R(n/N) − R theor )/R theor , that was shown in Figure 9.The data were described by a logarithmic dependence ∆R = A•ln(B•n/N), where A = 0.3 and B = 3.0.Note, that the argument 3n/N = n/(N/3) indicates that changes of resistivity occur in a reduced scale, ~(N/3).It is also worth to note, that the same type dependence was found for the in-plane components of the resistivity, ρ ||a and ρ ||b , in the conducting SCO complex [Fe{HC(pz) 3 } 2 ](TCNQ) 3 [40].Moreover, it was shown in [40] that the transverse transport, ρ ⊥ (T), did not react to the changes in the cation layers.Therefore, we suggested that in 1•MeOH all the observed changes of the resistance took place for the in-plane transport.Comparing data on Figure 9 one can conclude, that the major effect of SCO on transport properties occurs at lower concentrations n/N < 5%, whereas spin-spin depolarization is effective at n/N >> 5%.Similar effect was observed in [40].This is comprehensible and has a qualitative explanation.At the very low concentrations of HS [Fe(III)(sal 2 -trien)] complexes, the local displacements inside and/or between the TCNQ tetrads behave as randomly distributed defects.As the temperature approaches the midpoint of the transition, the anion sublattice turns into a modulated structure or to a certain extent into an amorphous one, which depends on heating/cooling rates and cooperative effects.Hence, a scattering rate gets less sensitive.On contrary, the efficiency of a cross relaxation 1/2 → 5/2 depends on a capacity of the spin reservoir of S = 5/2, i.e., concentration of the HS complexes within the interaction range.Thus, a momentum scattering and spin depolarization of conduction electrons in 1•MeOH are decoupled to a large extent.In other words, this allows spin manipulation by an external source.It was found in [40], that for delocalized spin moments in the conducting layer depolarization is even more efficient, ∆B ~(n/N) 2 .Once we obtain magnetically ordered local moments, we would get polarized conduction electrons.
The structures [Fe(III)(sal 2 -trien)](TCNQ)•CH 3 CN ( 2) and [Fe(III)(sal 2 -trien)](TCNQ)•H 2 O (3) could be viewed as the reference systems for distinguishing the effects of conduction electrons and solvent molecules in the magnetic and resonance properties of 1•MeOH.Solvent effects are responsible for the difference in their SCO transitions.The static magnetic susceptibility of 2 and 3 was measured in the range of 2-400 K four times, in two sequential cycles.Heating and cooling rates were taken the same as for 1•MeOH.The measurements begun at the room temperature in the external field B = 4.0 T. The virgin curves for the product χT are shown in Figures 10 and 11 (( ) data points) for 2 and 3, respectively.Upon cooling down to 2.0 K we measured the field dependences of the total magnetization by passing a complete field loop from +4.0 T → -5.0 T →+5.0 T →+0.1 T. The exposure at liquid helium temperature lasted several hours.Then the magnetic susceptibility was measured second time at the field value 0.1 T (∆).The measurements at heating were finished at 360 K, where the sample was exposed for about 10 min.A second cycle started at 360 K by measurements at cooling down to 2.0 K (∇) and further heating up to 360 K ( ).Repeating field measurements were not performed.The χT data for the heating curve of the cycle 2 in 3 are not shown in Figure 11, because we had to interrupt experiment due to technical reasons.We believe this was not critical for the discussion and conclusions.(□) data points) for 2 and 3, respectively.Upon cooling down to 2.0 K we measured the field dependences of the total magnetization by passing a complete field loop from +4.0 T → -5.0 T →+5.0 T →+0.1 T. The exposure at liquid helium temperature lasted several hours.Then the magnetic susceptibility was measured second time at the field value 0.1 T (Δ).The measurements at heating were finished at 360 K, where the sample was exposed for about 10 min.A second cycle started at 360 K by measurements at cooling down to 2.0 K () and further heating up to 360 K (○).Repeating field measurements were not performed.The χT data for the heating curve of the cycle 2 in 3 are not shown in Figure 11, because we had to interrupt experiment due to technical reasons.We believe this was not critical for the discussion and conclusions.It is worth to explain, how repeated measurements at the high and low fields, such as 4.0 T and 0.1 T, allow accounting magnetic impurities.If a sample contains impurities (for example iron oxide nanoparticles from the solvent) with a contribution χIMP, it is practically difficult to subtract it from the total magnetic response in experiment, χexp = χIMP + χsample(T).In an ideal case, when χexp coincides with the expected theoretical value, χsample theor , the experimental data are entirely related to the spin system of the studied compound.This took place for 1•MeOH, so we did not present the high field data.In other cases, proofs would be required.The absolute values of χexpT, 1.25 and 1.75 cm 3 •K•mol −1 , coincide for both 2 and 3.That means the absence of extrinsic magnetic contributions.We also asserted the absence of single radicals TCNQ •− as paramagnetic defects.Corresponding EPR signals [50] had negligible intensities, both in 2 and in 3 (Figure S8).These arguments as well as the measurement results allowed us to conclude, that the spin susceptibilities of 2 and 3 below 100 K were associated with the cation sublattice, namely with the magnetic moments of [Fe(III)(sal2-trien)] + complexes.The abrupt growth of χT above 100 K confirms SCO transition for both systems, which is in agreement with that for 1•MeOH.Meanwhile, in both systems, 2 and 3, the transition showed a thermal hysteresis.Intricate hysteresis of spin susceptibility is typical for mononuclear iron (III) SCO complexes [51].The absolute values of χT contained a significant quantity of residual HS complexes at temperatures below 100 K.The measurements in 2 yielded 80% of LS and 20% of HS complexes at 100 K in the cycle 1, and 69% of LS and 31% of HS complexes in the cycle 2 (Figure 10).For 3 the χT data were fitted with 76% of LS and 24% of HS complexes in the first cycle, and 68% to 32% in the second cycle, respectively (Figure 11).This is close to the balance ratio 2:1, which is evidently determined by a complete loss of solvent (2/3 complexes in LS state and 1/3-in HS state).Field dependences shown in Figures 10 and 11 also confirm, that considerable ferric entity remains in the HS state at T = 2 K. Magnetization curves were successfully fitted with 82% and 18% for 2 (solid line in Figure 10b, and 86% and 14% for 3 (solid line in Figure 11b).Note, that these numbers were obtained from original samples (before heating).Different volumes of residual HS fractions and invariable gap between the temperature curves χT in Figure 10   It is worth to explain, how repeated measurements at the high and low fields, such as 4.0 T and 0.1 T, allow accounting magnetic impurities.If a sample contains impurities (for example iron oxide nanoparticles from the solvent) with a contribution χ IMP , it is practically difficult to subtract it from the total magnetic response in experiment, χ exp = χ IMP + χ sample (T).In an ideal case, when χ exp coincides with the expected theoretical value, χ sample theor , the experimental data are entirely related to the spin system of the studied compound.This took place for 1•MeOH, so we did not present the high field data.
In other cases, proofs would be required.The absolute values of χ exp T, 1.25 and 1.75 cm 3 •K•mol −1 , measured respectively for 2 and 3 at 100 K, were found considerably higher than their theoretical estimates for the LS state.The question is, whether this was a total spin response of SCO system, or it included the extrinsic component χ IMP .Magnetic response of the solid particles usually does not depend on the external field of a measurement, whereas the actual signal χ sample is proportional to its strength.By performing measurement at a higher field, we reduced contribution of the extrinsic part, ~χIMP /B.Thus, a divergence of χ exp T data obtained at the various fields let us distinguish and account for external impurities, should they present in the sample.The negative side of the high field measurements is a descending trend, that arises due to insufficient population of upper spin state at higher Zeeman splitting (see the data points ( ) below 15 K in Figures 10 and 11).Therefore, a temperature behavior below 20 K in Figures 10 and 12 was not discussed for the data points obtained at 4.0 T. In the range near 100 K the values χ exp T measured at a broad span 0.1-4.0G pretty much coincide for both 2 and 3.That means the absence of extrinsic magnetic contributions.We also asserted the absence of single radicals TCNQ  S8).These arguments as well as the measurement results allowed us to conclude, that the spin susceptibilities of 2 and 3 below 100 K were associated with the cation sublattice, namely with the magnetic moments of [Fe(III)(sal 2 -trien)] + complexes.The abrupt growth of χT above 100 K confirms SCO transition for both systems, which is in agreement with that for 1•MeOH.Meanwhile, in both systems, 2 and 3, the transition showed a thermal hysteresis.Intricate hysteresis of spin susceptibility is typical for mononuclear iron (III) SCO complexes [51].The absolute values of χT contained a significant quantity of residual HS complexes at temperatures below 100 K.The measurements in 2 yielded 80% of LS and 20% of HS complexes at 100 K in the cycle 1, and 69% of LS and 31% of HS complexes in the cycle 2 (Figure 10).For 3 the χT data were fitted with 76% of LS and 24% of HS complexes in the first cycle, and 68% to 32% in the second cycle, respectively (Figure 11).This is close to the balance ratio 2:1, which is evidently determined by a complete loss of solvent (2/3 complexes in LS state and 1/3-in HS state).Field dependences shown in Figures 10 and 11 also confirm, that considerable ferric entity remains in the HS state at T = 2 K. Magnetization curves were successfully fitted with 82% and 18% for 2 (solid line in Figure 10b, and 86% and 14% for 3 (solid line in Figure 11b).Note, that these numbers were obtained from original samples (before heating).Different volumes of residual HS fractions and invariable gap between the temperature curves χT in Figure 10  The concentration of LS Fe(III)(sal2-trien) complexes in 2 and 3 was also verified by EPR.The measurements were carried out on the same samples that had passed through SQUID experiments.The absolute values, IEPR, at 100 K were found close to 2/3 of the theoretical estimate to within ~20% accuracy.The relative intensities IEPR(T)/IEPR(100 K) plotted in Figures 10 and 11 (right ordinate axes) were obtained in the heating regime.In the range 95-325 K respective spin concentration of LS Fe(III) decreased with the temperature to 60% of its original value for 2, and to 70% within 95-305 K range for 3. Temperature evolution of the lineshapes is shown in the Figures S7 and S8.The extracted gparameters are presented in Figure 12 ((a) for 2 and (b) for 3).The character of anisotropy, as well as the weak temperature dependence, were discussed earlier for 1•MeOH.A common peculiar feature, axial anisotropy, was described by two parameters g1 = g2 = g and g3 = g||.At 95 K they were 2.172(3) and 1.957(8) for 2, and 2.142(6) and 2.031(6) for 3.As the complexes underwent a thermal SCO, the EPR signal broadened, while the g-parameters remained unchanged.A partial peak-to-peak linewidth, ΔB, corresponding to g component was extracted and analyzed as a function of temperature like it had been done for 1•MeOH.The guiding lines for Bp+ and Bp− are shown in Figure 12.It was somewhat unclear, why at presence of 1/3 (33%) complexes in HS state in 2 and 3 the EPR signal of S = 1/2 moments has the linewidth similar to that at 5% of HS complexes in 1•MeOH.Consensus becomes apparent if we assume spatial inhomogeneity.Spin-spin relaxation becomes ineffective if the fraction with S = 1/2 moments of LS Fe(III) and the HS fraction are separated in the cation layers.Besides that, insulating TCNQ dimers in 2 and 3 are unable to facilitate expanded interactions by mediating cross-relaxation via conduction electrons.At the same time, the LS fraction undergoes SCO transition, while the HS fraction remains unchanged.Indistinct n/N estimates together with less effective line broadening in 2, did not result in finding an appropriate correlation between EPR linewidth and concentration of Fe(III) ions in HS state.However, a credible linkage between line broadening and χT growth was found in the system 3.The transition was extrapolated by Boltzmann model.The exponential growth was superimposed by an ascending foothill segment (Figure S9).This feature arises because the thermal hysteresis is a kinetic effect, and spin states appear in result of cooperative interactions [52].Spin concentration of S = 5/2, x, was obtained by solving the equation at every temperature point of the measurement, χexp(T) = (1 − x)•χ1/2 theor + x•χ5/2 theor , where χ1/2 theor and χ5/2 theor -theoretical values of the magnetic susceptibilities at 100% concentrations of LS or HS complexes, respectively.This procedure might seem arguable, because each data set was obtained in the different conditions, due to thermal irreversibility (spin fractionation).However, we believe The concentration of LS Fe(III)(sal 2 -trien) complexes in 2 and 3 was also verified by EPR.The measurements were carried out on the same samples that had passed through SQUID experiments.The absolute values, I EPR , at 100 K were found close to 2/3 of the theoretical estimate to within ~20% accuracy.The relative intensities I EPR (T)/I EPR (100 K) plotted in Figures 10 and 11 (right ordinate axes) were obtained in the heating regime.In the range 95-325 K respective spin concentration of LS Fe(III) decreased with the temperature to 60% of its original value for 2, and to 70% within 95-305 K range for 3. Temperature evolution of the lineshapes is shown in the Figures S7 and S8.The extracted g-parameters are presented in Figure 12 ((a) for 2 and (b) for 3).The character of anisotropy, as well as the weak temperature dependence, were discussed earlier for 1•MeOH.A common peculiar feature, axial anisotropy, was described by two parameters g 1 = g 2 = g ⊥ and g 3 = g || .At 95 K they were 2.172(3) and 1.957(8) for 2, and 2.142(6) and 2.031(6) for 3.As the complexes underwent a thermal SCO, the EPR signal broadened, while the g-parameters remained unchanged.A partial peak-to-peak linewidth, ∆B, corresponding to g ⊥ component was extracted and analyzed as a function of temperature like it had been done for 1•MeOH.The guiding lines for B p+ and B p− are shown in Figure 12.
The ∆B values at 100 K for 2 and 3 equaled 43.0 G and 45.3 G, which was broader than 18.8 G for 1•MeOH.Meanwhile, in 1•MeOH it quickly reached the value 46 G at HS Fe(III) concentration of 5%.It was somewhat unclear, why at presence of 1/3 (33%) complexes in HS state in 2 and 3 the EPR signal of S = 1/2 moments has the linewidth similar to that at 5% of HS complexes in 1•MeOH.Consensus becomes apparent if we assume spatial inhomogeneity.Spin-spin relaxation becomes ineffective if the fraction with S = 1/2 moments of LS Fe(III) and the HS fraction are separated in the cation layers.
Besides that, insulating TCNQ dimers in 2 and 3 are unable to facilitate expanded interactions by mediating cross-relaxation via conduction electrons.At the same time, the LS fraction undergoes SCO transition, while the HS fraction remains unchanged.Indistinct n/N estimates together with less effective line broadening in 2, did not result in finding an appropriate correlation between EPR linewidth and concentration of Fe(III) ions in HS state.However, a credible linkage between line broadening and χT growth was found in the system 3.The transition was extrapolated by Boltzmann model.The exponential growth was superimposed by an ascending foothill segment (Figure S9).This feature arises because the thermal hysteresis is a kinetic effect, and spin states appear in result of cooperative interactions [52].Spin concentration of S = 5/2, x, was obtained by solving the equation at every temperature point of the measurement, theor , where χ 1/2 theor and χ 5/2 theor -theoretical values of the magnetic susceptibilities at 100% concentrations of LS or HS complexes, respectively.This procedure might seem arguable, because each data set was obtained in the different conditions, due to thermal irreversibility (spin fractionation).However, we believe that the linkage between EPR linewidth and spin concentration of S = 5/2 moments remains valid, if the latter was extracted from the χT data at heating in the first measurement cycle (data points (∆) in Figure 11).Indeed, it follows from Figure 13 that the value ∆B LS increases proportionally to n/N above the threshold of 23%.Similar dependence took place in 1•MeOH right at the beginning of the reversible SCO transition.Hence, a broadening part of the dependence in 3, ∆B ~k•(n/N) (solid line in Figure 13a), where k = 30.1 G/%, occurs due to spin-spin relaxation.This contribution above the threshold is induced by thermally activated HS complexes with S = 5/2, homogeneously appearing inside the LS fraction.For reference, in Figure S10 we present the temperature dependences ∆B LS (T) in 2 and 3 with respective fitting parameters in Table S4.The threshold is consistent with spatial in homogeneity of a spin ensemble, when spin probes of S = 1/2 do not exist inside the intact fraction of S = 5/2 moments.The residual (intact) HS fraction remains spatially separated.EPR spectra from Fe(III) ions in the HS state in 2 and 3 are shown in Figure S8.Very weak temperature dependences of g-factor with g = 5.6(0) at 300 K were observed for 2, and with g = 4.06-for 3, respectively.Line broadening in the SCO range, ∆B HS (T), is shown in Figure 13b.The linewidth was found to follow a linear dependence ∆B HS = ∆B HS 0 + k HS •T, where ∆B HS 0 = 5.66•10 2 G, and k HS = 1.7 G/K for 2, and ∆B HS 0 = 5.09•10 2 G, and k HS = 0.91 G/K for 3. A proportional broadening ∆B HS (T) ~T was observed in the entire range of SCO in 3, whereas in 2 it begun at higher temperatures, which correlates with the data in Figure 12a and in further discussed Figure 14a.
An assumption of a phase fractionation, promoted by a solvent loss, implies a partition of the total spin susceptibility to the fractional contributions and further analysis of their individual temperature behavior.The magnetic susceptibility of the LS phase in 2 was reconstructed by using the expression, χ LS T = (χ exp (T) -y•χ 5/2 theor )•T/(1 -y), where y is the concentration of the moments S = 5/2 in the HS phase.Similarly, the expression χ HS T = (χ exp (T) − (1 − y)•χ 1/2 theor )•T/y described the contribution of the HS phase.Here we took g = 2.00 and 2.17 for calculated values of χ 1/2 theor and χ 5/2 theor , respectively.A comparative plot "χ LS T vs. T" in Figure 14 demonstrates, that the thermal cycling affects a foothill domain of the SCO transition curve, leaving its steep slope unchanged.EPR line broadening correlates with the steep part, but remains insensitive to a foothill region (see also Figure S9).This is consistent with the threshold in 3 (Figure 13), that was attributed to the cooperative interactions [52].The broadening is caused by spin-spin interactions between S = 1/2 spin probes with magnetic moments S = 5/2 in the neighboring HS [Fe(III)(sal 2 -trien)] + complexes, which do not belong to the residual HS fraction (HS phase).After several cycles of measurements, the amount of [Fe(III)(sal 2 -trien)] + complexes, taking part in the SCO transition, approached 2/3 of total.
in Figure 13a), where k = 30.1 G/%, occurs due to spin-spin relaxation.This contribution above the threshold is induced by thermally activated HS complexes with S = 5/2, homogeneously appearing inside the LS fraction.For reference, in Figure S10 we present the temperature dependences ΔBLS(T) in 2 and 3 with respective fitting parameters in Table S4.The threshold is consistent with spatial in homogeneity of a spin ensemble, when spin probes of S = 1/2 do not exist inside the intact fraction of S = 5/2 moments.The residual (intact) HS fraction remains spatially separated.EPR spectra from Fe(III) ions in the HS state in 2 and 3 are shown in Figure S8.Very weak temperature dependences of g-factor with g = 5.6(0) at 300 K were observed for 2, and with g = 4.06-for 3, respectively.Line broadening in the SCO range, ΔBHS(T), is shown in Figure 13b.The linewidth was found to follow a linear dependence ΔBHS = ΔB HS 0 + kHS•T, where ΔB HS 0 = 5.66•10 2 G, and kHS = 1.7 G/K for 2, and ΔB HS 0 = 5.09•10 2 G, and kHS = 0.91 G/K for 3. A proportional broadening ΔBHS(T) ~ T was observed in the entire range of SCO in 3, whereas in 2 it begun at higher temperatures, which correlates with the data in Figure 12a and in further discussed Figure 14a.
An assumption of a phase fractionation, promoted by a solvent loss, implies a partition of the total spin susceptibility to the fractional contributions and further analysis of their individual temperature behavior.The magnetic susceptibility of the LS phase in 2 was reconstructed by using the expression, χLST = (χexp(T) -y•χ5/2 theor )•T/(1 -y), where y is the concentration of the moments S = 5/2 in the HS phase.Similarly, the expression χHST = (χexp(T) − (1 − y)•χ1/2 theor )•T/y described the contribution of the HS phase.Here we took g = 2.00 and 2.17 for calculated values of χ1/2 theor and χ5/2 theor , respectively.A comparative plot "χLST vs. T" in Figure 14 demonstrates, that the thermal cycling affects a foothill domain of the SCO transition curve, leaving its steep slope unchanged.EPR line broadening correlates with the steep part, but remains insensitive to a foothill region (see also Figure S9).This is consistent with the threshold in 3 (Figure 13), that was attributed to the cooperative interactions [52].The broadening is caused by spin-spin interactions between S = 1/2 spin probes with magnetic moments S = 5/2 in the neighboring HS [Fe(III)(sal2-trien)] + complexes, which do not belong to the residual HS fraction (HS phase).After several cycles of measurements, the amount of [Fe(III)(sal2-trien)] + complexes, taking part in the SCO transition, approached 2/3 of total.The total spin response χT of 2 demonstrated unusual peaks at T = 4.11 K shown in Figure 10.System 3 also showed a significant growth of χT, though with no maximum, at temperatures T < 10 K (Figure 11).In this study, we attributed the observed peaks to the HS phase, χHST, as shown in Figure 14b.We incline to discuss their ferromagnetic nature rather than metastable high spin trapping.Our opinion is based on the reasoning below.Enhancement of a spin response at the low temperatures is often a sign of a thermal-induced metastable spin-state trapping (TIESST) [53].Indeed, some mononuclear Fe(III) complexes, especially those with two-step spin-crossover, exhibit a metastable spin-state HS* emerging after thermal quench from the true HS state of a complete SCO transition.A key feature is the rate of cooling, which determines the amount of quenched fraction.A drop of χHS*T to the original LS level is, in fact, a relaxation of a metastable phase, while heating above the temperature, that corresponds to a respective energy barrier.A width of the transition depends on a heating rate, and usually consists tens of degrees in the transition range from 50 to 100 K [48].A distinct feature of TIESST phase is the butte-like shape of the χT curve, where the flat top may reach ~50% of the total HS response, depending on a freezing rate.Here are the details, that contradict a mechanism of quenched HS* spinstate for 2: (1) gradual refrigeration down to 2.0 K took 5 h (including measurement time), which is ~100 times slower than the flash freezing in SQUID chamber; (2) the difference (χHST-χLST) at maximum exceeded 100% of the χLST gain at 353 K; (3) χHST curves had no plateau and further relaxation drop, associated with decaying HS*spin-state; (4) gradual 50% decrease of χHST between 4 and 10 K was unlikely caused by released dynamics (thermal relaxation) of sal2-trien ligand; (5) maximum value of χT in Figure 14b did grow at cycle 2, whereas LS phase has diminished.However, since the complete SCO transition has not been reached, we cannot reject this scenario completely.Importantly, the TIESST mechanism does not imply magnetic coupling in the HS* system.The total spin response χT of 2 demonstrated unusual peaks at T = 4.11 K shown in Figure 10.System 3 also showed a significant growth of χT, though with no maximum, at temperatures T < 10 K (Figure 11).In this study, we attributed the observed peaks to the HS phase, χ HS T, as shown in Figure 14b.We incline to discuss their ferromagnetic nature rather than metastable high spin trapping.Our opinion is based on the reasoning below.
Enhancement of a spin response at the low temperatures is often a sign of a thermal-induced metastable spin-state trapping (TIESST) [53].Indeed, some mononuclear Fe(III) complexes, especially those with two-step spin-crossover, exhibit a metastable spin-state HS* emerging after thermal quench from the true HS state of a complete SCO transition.A key feature is the rate of cooling, which determines the amount of quenched fraction.A drop of χ HS* T to the original LS level is, in fact, a relaxation of a metastable phase, while heating above the temperature, that corresponds to a respective energy barrier.A width of the transition depends on a heating rate, and usually consists tens of degrees in the transition range from 50 to 100 K [48].A distinct feature of TIESST phase is the butte-like shape of the χT curve, where the flat top may reach ~50% of the total HS response, depending on a freezing rate.Here are the details, that contradict a mechanism of quenched HS* spin-state for 2: (1) gradual refrigeration down to 2.0 K took 5 h (including measurement time), which is ~100 times slower than the flash freezing in SQUID chamber; (2) the difference (χ HS T-χ LS T) at maximum exceeded 100% of the χ LS T gain at 353 K; (3) χ HS T curves had no plateau and further relaxation drop, associated with decaying HS*spin-state; (4) gradual 50% decrease of χ HS T between 4 and 10 K was unlikely caused by released dynamics (thermal relaxation) of sal 2 -trien ligand; (5) maximum value of χT in Figure 14b did grow at cycle 2, whereas LS phase has diminished.However, since the complete SCO transition has not been reached, we cannot reject this scenario completely.Importantly, the TIESST mechanism does not imply magnetic coupling in the HS* system.
The alternative scenario suggests the magnetic exchange interactions in the HS phase.A maximum in χT curve at helium temperatures is a characteristic feature of ferromagnetic or ferrimagnetic coupling in many molecular magneto-active systems, including HS Fe(III) complexes [52,54].Ferromagnetic interactions were also found in numerous TCNQ based compounds with metallo-complex counerions [55].Ferro-and ferrimagnetic coupling was reported for networks of transition metals, bridged by TCNQ −1 radicals M(TCNQ) 2 (M = Mn, Fe, Co, Ni), among which was the compound Mn(TCNQ) 2 with T c = 44 K [56,57].
In the crystal field approach, the structural characteristics of 1•MeOH, 2 and 3 deny exchange coupling between the magnetic moments of Fe(III) in neighboring complexes.Moreover, the shortest distances d Fe-Fe = 7.070 Å (350 K, 1•MeOH) and 7.729 Å (325 K, 2) are not sufficient for the effective spin-spin relaxation via dipole-dipole mechanism.Therefore, magnetic and resonance behavior of the Fe(III) magnetic moments, found in the current study, implies the key role of TCNQ molecules as a mediator.Schematic structural arrangements of the TCNQ molecules and [Fe(  (2).Thus, an interplay between π bonding and intramolecular distortion could be the driving force for the exchange coupling.Note that that the acetonitrile molecule in 2 acts as an "anchor" due to short contact with N amine (N3, Figure 15B).Therefore, the solvent removal may release the distortions.However, a key role belongs to the electrons in the TCNQ sublattice.Hopping electrons of semiconducting TCNQ layers in 1•MeOH could serve as an efficient spin reservoir for cross relaxation between S = 1/2 and S = 5/2 local moments (Figure 15A).High conductivity along TCNQ stacks is capable to provide long-range spin-spin interactions.Localized spins of adjacent radicals TCNQ •− in weakly conducting compounds 2 and 3 suggest intrinsic interactions in the form of singlet-triplet splitting or exchange coupling.Taking into consideration π-stacking between phenoxy groups of neighboring ligands in Figure 15B, one can consider an alternating spin chain with the units or [-5/2-1/2-1/2-5/2-] appearing upon passing the SCO transition.There were several drawbacks.We could not simulate the temperature dependences in Figure 14b by a Heisenberg model for the individual four spin unit with two exchange constant J 1 , J 2 and one ZFS parameter D. We also did not find a coercivity in magnetization curves in Figures 11 and 12. Though the latter can be understood, as the measurements were performed before heating-cooling cycles.Thus, the discussion of physical properties in this study was extended with somewhat hypothetical consideration of fundamental mechanisms.We believe, that holistically, this will draw attention to the important aspects of spin-spin interactions and electronic transport in [Fe(III)(sal 2 -trien)](TCNQ) 2 (1•MeOH) system as well as an interplay between ferro/ferrimagnetic coupling and SCO in [Fe(III)(sal 2 -trien)](TCNQ) (2, 3).While the first compound has a record conductivity up-to-date and reversible SCO transition, the compound [Fe(III)(sal 2 -trien)](TCNQ)•CH 3 CN is the first system, where ferromagnetic coupling was triggered by the SCO transition.In other words, in the system [Fe(III)(sal 2 -trien)](TCNQ) n switchable magnetic moments of isolated metallo-complexes coexist and interact with spin system of organic network in full or fractional reduction state.In turn, the organic network allows electronic/spin transport.The fact that the spin interactions are switchable, makes this conducting system a prospective candidate for molecular spintronics.

X-ray Crystallography
Single crystal X-ray diffraction experiments were carried out on a Bruker SMART APEX2 CCD diffractometer (Bruker AXS Advanced X-ray Solutions GmbH, Karlsruhe, Germany) (for 1•MeOH at 100, 220 K and 2 at 100 K) and an Oxford Diffraction Gemini-R CCD diffractometer (Oxford Diffraction, Oxford, Oxfordshire, United Kingdom) [for all other temperatures, λ(MoKα) = 0.71073 Å, graphite monochromator, ω-scan mode].Multi temperature experiment procedure was used without control of cooling/warming rate; the average rate was roughly estimated as 2-4 K/min.Data collection for 1•MeOH and 1 proceeded during 24 h at each temperature.Crystal of 1•MeOH had begun to lose solvent at the end of the 340 K experiment; in this reason some latest X-ray frames were excluded from the data list for correct refinement of the solvate structure.Then the crystal was warmed to 350 K and kept at this temperature during one day before data collection to obtain the new phase 1 totally free of solvent.The structures were solved by the direct method and refined by the full-matrix least-squares technique against F 2 in an anisotropic approximation for all non-hydrogen atoms.Hydrogen atoms were localized from the Fourier synthesis of the electron density and refined in the isotropic approximation.MeOH sites in 1•MeOH were finally refined as fully occupied at 295 and 340 K though occupancy refinement gave values of 0.988(4) and 0.956 (6), respectively.All calculations were performed using SHELXTL PLUS 5.0 and SHELX-2016 program packages [59].Selected crystallographic data and refinement parameters are given in Table 3.The full data of studies are available at the Cambridge Crystallographic Data Centre.Distortion of the lattice due to solvent loss transforms acute unit cell of 1•MeOH into obtuse one for 1, the latter is given in non-standard setting both in Table 3 and cif -files for direct comparison of the structure data.
[110] direction.The stacks are composed of tetrads . . .-(II-I-I-II)-(II-I-I-II)-. . .(Figure 2) with ring-over-bond overlap inside the tetrad due to longitudinal shift of adjacent anions and ring-over-ring one between the tetrads (longitudinal shift is absent, Figure S4).Interplane separations in the stack at 100 K are 3.274(4) Å (I. . .I), 3.19(4) Å (I. . .II) and 3.297(6) Å (II. . .II), the angle between mean planes of I and II is 2.37(5) • .There are a lot of intrastack C. . .C contacts shorter than 3.4 Å, while adjacent stacks interact only by weak hydrogen contacts of C-H. . .N type with H. . .N distances of 2.58-2.73Å.The [Fe(III)(sal 2 -trien)] + cations are located between the TCNQ layers and connected by hydrogen contacts with surrounding TCNQ and methanol molecules.The Fe(III) ion has a slightly distorted octahedral coordination geometry formed by two oxygen and four nitrogen atoms of ligand.Key bond distances and angles in the Fe(III)(sal 2 -trien) + cation at 100 K are listed in Table S2.Mean values of the Fe-O, Fe-N imine and Fe-N amine bond distances in the FeN 4 O 2 chromophore (1.885(8), 1.936(8) and 2.007(2) Å, respectively) and shape of the cation with the dihedral angle α between two phenoxy groups of 53.87(6) • (Figure 3a) correspond to the low-spin state of the Fe(III) ion.The diagonals of the octahedron at Fe centre are close to linear (N-Fe-O and N-Fe-N angles are 174 • -177 • , Table 35 and 2.52 Å and three C-Hcat…NII contacts with H…N = 2.44, 2.53 and 2.65 Å and only two C-H…NI contacts to TCNQ I with longer H…N distances of 2.69 and 2.73 Å.The [Fe(III)(sal2-trien)] + cations are also hydrogen bonded to the methanol solvent molecules (N-Hcat…Osolv, C-Hcat…Osolv and Ocat…H-Osolv contacts with H…O distances of 2.04, 2.46 and 1.80 Å, respectively).

of 26 Figure 4 .
Figure 4. Temperature dependences of mean Fe-Namine and Fe-Nimine bond length values and dihedral angle α between two phenoxy arms in the cation (a); TCNQ •δ− I and II charges (b) and interplane separations I-I, I-II and II-II in TCNQ stack (c).Starting point was at 100 K, warming from 100 to 385 K was followed by cooling from 385 to 220 K (see arrows) and new warming from 220 to 325 K (points corresponding to second warming are shown by empty marks).The data points obtained while warming in the range 100 K to 340 K belong to the solvate phase 1•MeOH.The parameters obtained while warming above 350 K and further cooling correspond to the solvent-free phase 1.

Figure 4 .
Figure 4. Temperature dependences of mean Fe-N amine and Fe-N imine bond length values and dihedral angle α between two phenoxy arms in the cation (a); TCNQ •δ− I and II charges (b) and interplane separations I-I, I-II and II-II in TCNQ stack (c).Starting point was at 100 K, warming from 100 to 385 K was followed by cooling from 385 to 220 K (see arrows) and new warming from 220 to 325 K (points corresponding to second warming are shown by empty marks).The data points obtained while warming in the range 100 K to 340 K belong to the solvate phase 1•MeOH.The parameters obtained while warming above 350 K and further cooling correspond to the solvent-free phase 1.

Figure 4 .
Figure 4. Temperature dependences of mean Fe-Namine and Fe-Nimine bond length values and dihedral angle α between two phenoxy arms in the cation (a); TCNQ •δ− I and II charges (b) and interplane separations I-I, I-II and II-II in TCNQ stack (c).Starting point was at 100 K, warming from 100 to 385 K was followed by cooling from 385 to 220 K (see arrows) and new warming from 220 to 325 K (points corresponding to second warming are shown by empty marks).The data points obtained while warming in the range 100 K to 340 K belong to the solvate phase 1•MeOH.The parameters obtained while warming above 350 K and further cooling correspond to the solvent-free phase 1.

Figure 6 .
Figure 6.Temperature dependence of the normalized resistance, R(T)/R(300 K) for 1•MeOH, in the logarithmic scale.Inset-the plot of R(T)/R(300 K) vs. scaled reciprocal temperature, 1000/T.Solid line is the best fit curve in the range 78 K < T < 110 K, Rtheor(T) = exp(ΔE/kT), with the energy gap ΔE = 0.05 eV.Zoomed window-a deviation of the experimental data from the theoretical curve starting above 105 K.

Figure 6 .
Figure 6.Temperature dependence of the normalized resistance, R(T)/R(300 K) for 1•MeOH, in the logarithmic scale.Inset-the plot of R(T)/R(300 K) vs. scaled reciprocal temperature, 1000/T.Solid line is the best fit curve in the range 78 K < T < 110 K, R theor (T) = exp(∆E/kT), with the energy gap ∆E = 0.05 eV.Zoomed window-a deviation of the experimental data from the theoretical curve starting above 105 K.

Figure 7 .
Figure 7. Temperature dependences (logarithmic scale) of the product χT measured for 1•MeOH by SQUID magnetometer (B = 1 kG, heating (∆) and cooling (∇) regimes) and the relative spin concentration I EPR (T)/I EPR (100 K) measured by EPR (B = 3 kG, heating ( ) regime) (a).Inset-detailed evolution of the χT in the range of spin crossover transition between LS Fe(III) S = 1/2 and HS Fe(III) S = 5/2.Solid line is a simulation by a Boltzmann distribution.Field dependence of the magnetization, M(B), measured at T = 2.0 K (b).Solid line denotes the Brillouin function for S = 1/2.

Figure 8 .
Figure 8. Temperature evolution of the EPR spectrum of 1•MeOH in the range 100 K to 370 K (a).Temperature dependence of the g-values (g1 (◊), g2 (□), g3 (○)) measured while heating in the range of SCO in 1•MeOH: (b).Solid straight lines are linear fits.Solid bold lines connecting squares (■) designate positions of the peaks Bp+ and Bp− for a simulated central line of the spectrum corresponding to g2.The peak-to-peak linewidth of the central line is ΔB = Bp+ − Bp−.

Figure 8 .
Figure 8. Temperature evolution of the EPR spectrum of 1•MeOH in the range 100 K to 370 K (a).Temperature dependence of the g-values (g 1 (♦), g 2 ( ), g 3 ( )) measured while heating in the range of SCO in 1•MeOH: (b).Solid straight lines are linear fits.Solid bold lines connecting squares ( ) designate positions of the peaks B p+ and B p− for a simulated central line of the spectrum corresponding to g 2 .The peak-to-peak linewidth of the central line is ∆B = B p+ − B p− .
large extent.In other words, this allows spin manipulation by an external source.It was found in[40], that for delocalized spin moments in the conducting layer depolarization is even more efficient, ΔB ~ (n/N)2 .Once we obtain magnetically ordered local moments, we would get polarized conduction electrons.The structures [Fe(III)(sal2-trien)](TCNQ)•CH3CN (2) and [Fe(III)(sal2-trien)](TCNQ)•H2O (3) could be viewed as the reference systems for distinguishing the effects of conduction electrons and solvent molecules in the magnetic and resonance properties of 1•MeOH.Solvent effects are responsible for the difference in their SCO transitions.The static magnetic susceptibility of 2 and 3 was measured in the range of 2-400 K four times, in two sequential cycles.Heating and cooling rates were taken the same as for 1•MeOH.The measurements begun at the room temperature in the external field B = 4.0 T. The virgin curves for the product χT are shown in Figures10 and 11 (

Figure 12 .
Figure 12.Temperature dependences of the g-parameters for 2 (a) and 3 (b) extracted from LS Fe(III) EPR signal while measurements in heating regime: g (◊), g|| (□).Solid lines connecting squares (■) fit the positions of the peaks Bp+ and Bp− for simulation the partial spectral line corresponding to g, ΔB= Bp+ − Bp− .

Figure 12 .
Figure 12.Temperature dependences of the g-parameters for 2 (a) and 3 (b) extracted from LS Fe(III) EPR signal while measurements in heating regime: g ⊥ (♦), g || ( ).Solid lines connecting squares ( ) fit the positions of the peaks B p+ and B p− for simulation the partial spectral line corresponding to g ⊥ , ∆B = B p+ − B p− .

Figure 14 .
Figure 14.Temperature dependences of the product χLST for the LS (SCO) fraction in 2 (measurements in heating regime, cycle 1 (Δ), cycle 2 (•)) (a).Temperature dependence of the EPR linewidth for g component measured at heating in cycle 1 (□).Temperature dependences of the product χHST for the HS (residual) fraction (b).

Figure 14 .
Figure 14.Temperature dependences of the product χ LS T for the LS (SCO) fraction in 2 (measurements in heating regime, cycle 1 (∆), cycle 2 ( )) (a).Temperature dependence of the EPR linewidth for g ⊥ component measured at heating in cycle 1 ( ).Temperature dependences of the product χ HS T for the HS (residual) fraction (b).

Figure 15 .
Figure 15.Schematic structural arrangements of the TCNQ molecules and [Fe(sal2-trien)] + cation complexes in 1•MeOH (A) and 2 (B).Dotted lines represent N…N, C…C short contacts of less than the sum of the van der Waals radii; (C) Mutual arrangement of coordination Fe(III) octahedra with neighbor TCNQ units II in 1•MeOH.Orbitals are presented schematically.In the crystal field approach, the structural characteristics of 1•MeOH, 2 and 3 deny exchange coupling between the magnetic moments of Fe(III) in neighboring complexes.Moreover, the shortest distances dFe-Fe = 7.070 Å (350 K, 1•MeOH) and 7.729 Å (325 K, 2) are not sufficient for the effective spin-spin relaxation via dipole-dipole mechanism.Therefore, magnetic and resonance behavior of the magnetic moments, found in the current study, implies the key role of TCNQ molecules as a mediator.Schematic structural arrangements of the TCNQ molecules and [Fe(sal2-trien)] + complexes in 1•MeOH (A) and 2 (B) are shown in Figure 15.Dotted lines represent short contacts N…N, C…C, having distances less than the sum of their van der Waals radii including effect of a high spin Fe dilation.In 1•MeOH at ambient conditions the distances Namine…NCN were 3.128 Å (N4,

Figure 15 .
Figure 15.Schematic structural arrangements of the TCNQ molecules and [Fe(sal 2 -trien)] + cation complexes in 1•MeOH (A) and 2 (B).Dotted lines represent N. . .N, C. . .C short contacts of less than the sum of the van der Waals radii; (C) Mutual arrangement of coordination Fe(III) octahedra with neighbor TCNQ units II in 1•MeOH.Orbitals are presented schematically.

Table 2 .
The average Fe-N imine and Fe-N amine bond length values and dihedral angle α between two phenoxy groups in [Fe(III)(sal 2 -trien)] + cation in 2 and interplanar separations d in the TCNQ dimer at different temperatures.