Ashlyn A.
Kamin
a,
Ian P.
Moseley
b,
Jeewhan
Oh
c,
E. J.
Brannan
a,
Paige M.
Gannon
a,
Werner
Kaminsky
a,
Joseph M.
Zadrozny
b and
Dianne J.
Xiao
*a
aDepartment of Chemistry, University of Washington, Seattle, Washington 98195, USA. E-mail: djxiao@uw.edu
bDepartment of Chemistry, Colorado State University, Fort Collins, Colorado 80523, USA
cDepartment of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, USA
First published on 27th March 2023
Redox-active tetraoxolene ligands such as 1,4-dihydroxybenzoquinone provide access to a diversity of metal–organic architectures, many of which display interesting magnetic behavior and high electrical conductivity. Here, we take a closer look at how structure dictates physical properties in a series of 1D iron–tetraoxolene chains. Using a diphenyl-derivatized tetraoxolene ligand (H2Ph2dhbq), we show that the steric profile of the coordinating solvent controls whether linear or helical chains are exclusively formed. Despite similar ligand environments, only the helical chain displays temperature-dependent valence tautomerism, switching from (FeII)(Ph2dhbq2−) to (FeIII)(Ph2dhbq3˙−) at temperatures below 203 K. The stabilization of ligand radicals leads to exceptionally strong magnetic exchange coupling (J = −230 ± 4 cm−1). Meanwhile, the linear chains are more amenable to oxidative doping, leading to Robin–Day class II/III mixed-valency and an increase in electrical conductivity by nearly three orders of magnitude. While previous studies have focused on the effects of changing metal and ligand identity, this work highlights how altering the metal–ligand connectivity can be a similarly powerful tool for tuning materials properties.
Extended metal–dioxolene materials are typically accessed by fusing two or more dioxolene units to form a single conjugated polytopic ligand. Well-known examples include 1,4-dihydroxybenzoquinone (H2dhbq),16,17 which contains two fused dioxolene units (also known as tetraoxolene) (Fig. 1), and 2,3,6,7,10,11-hexahydroxytriphenylene (HHTP),10 which contains three fused dioxolene units. We have been interested in studying extended 1D metal–tetraoxolene chains with the goal of identifying structure–property relationships that can be translated to higher dimensional 2- and 3D networks. Simple 1D chains are ideally suited for probing how changes to the metal, ligand, and metal–ligand connectivity influence electronic and magnetic coupling across the metal–tetraoxolene linkages. However, a barrier to the systematic synthesis and structural modification of such materials is the complex coordination chemistry of tetraoxolene ligands. For example, while H2dhbq and its substituted derivatives (H2R2dhbq, where R = F, Cl, Br, NO2, etc.) have been shown to bind metals to form 1-, 2-, and 3D networks (Fig. 1),16 how to select for each structure remains poorly understood. In the case of iron, 3D cubic frameworks,11 2D honeycomb sheets,8,18–25 1D linear chains,49,50 and 1D zig-zag chains6,26 have all been reported.
Fig. 1 (a) Binding modes for H2R2dhbq ligands and the resulting metal–organic structures. (b) Redox series of R2dhbqn− ligands. The ligand contains a radical in the R2dhbq3˙− state. |
Here we show how new linear and helical iron–tetraoxolene chain architectures can be achieved using a more sterically encumbered ligand, H2Ph2dhbq (2,5-dihydroxy-3,6-dipheny-1,4-benzoquinone). Both linear and helical structures can be predictably accessed simply by adjusting the steric bulk of the coordinating solvent. Despite identical bridging ligands, switching the chain structure dramatically alters how charges are distributed across the metal–tetraoxolene linkage. Unlike the linear chain, which remains (FeII)(Ph2dhbq2−) at all observed temperatures, the helical chain undergoes valence tautomerization, adopting a (FeIII)(Ph2dhbq3˙−) configuration at temperatures below 203 K. The stabilization of ligand radicals leads to exceptionally strong antiferromagnetic metal–ligand exchange coupling (J = −230 ± 4 cm−1). Further, we show that the linear chains are more easily oxidatively doped, leading to orders-of-magnitude improvement in electrical conductivity. Together, this work illustrates how changing the metal–ligand connectivity can be a powerful, yet often overlooked, tool to alter materials properties.
The greater steric bulk of the diphenyl substituted ligand provides an opportunity to probe the nuances of 1D chain formation without interference from competing 2D and 3D phases. With H2Ph2dhbq in hand, we were able to carry out a comprehensive survey of different synthesis conditions that led to two important structural findings (Fig. 2). First, we found that the steric bulk of the solvent employed during synthesis dictates the overall chain architecture. Second, we observed that the initial oxidation states of the metal and ligand starting materials greatly impact crystallite size and compositional purity. As we will describe in subsequent sections, these differences in chain architecture lead to significant differences in the physical properties of these chains, including the electrical conductivity and magnetism.
The first surprising structural outcome we observed was that distinct helical and linear chain architectures could be achieved by subtly tuning the steric bulk of the coordinating solvent (Fig. 2). Specifically, heating H2Ph2dhbq and FeCl2 in dimethylformamide (DMF) at 120 °C under an inert atmosphere provides the helical cis-Fe(Ph2dhbq)(DMF)2 chain (Fig. 2a and b). In this structure, the quasi-octahedral iron sites are coordinated by two cis DMF molecules and bridged by bis-bidentate Ph2dhbq ligands to form a helix. To our knowledge, these represent the first chiral 1D chains synthesized from H2dhbq or its substituted derivatives. On the other hand, changing the solvent from DMF to a slightly bulkier dimethylacetamide (DMA) provides trans-Fe(Ph2dhbq)(DMA)2 (Fig. 2c). While the iron sites remain quasi-octahedral in this structure, the bulkier DMA molecules now adopt a trans configuration, leading to linear rather than helical chains.
Given that both FeII and H2Ph2dhbq are redox-active, we next probed how the initial metal and ligand oxidation states impact chain formation. Previous reports have shown that changing the initial ligand oxidation state can lead to improved crystallinity or even entirely new phases.33,34 We hypothesized that similar chain structures could be achieved by combining the fully reduced tetrahydroxy ligand, H4Ph2dhbq, with an FeIII source, as FeIII can bind the ligand as well as serve as an in situ oxidant. Indeed, heating H4Ph2dhbq and FeCl3 in DMF under an inert atmosphere led to the formation of single crystals that are structurally and compositionally identical to those obtained from H2Ph2dhbq/FeII, but with significantly larger crystallite sizes (Fig. 2a). Given the improved crystallinity, all subsequent characterization of cis-Fe(Ph2dhbq)(DMF)2 will focus on samples prepared using H4Ph2dhbq/FeIII.
In contrast to the helical chain, which remains compositionally pure despite changes to the initial metal and ligand oxidation states, the linear chain is more readily over-oxidized. Heating H4Ph2dhbq and FeBr3 at 120 °C in DMA under an inert atmosphere leads to oxidative doping of the 1D chain and incorporation of additional charge-balancing Br− anions. Single crystal X-ray crystallography confirmed the formation of the desired linear chain geometry, and detailed elemental analysis revealed the formula to be trans-Fe(Ph2dhbq)(DMA)2Br0.55 (Fig. 2d). Although the bromine atoms are not observed directly in the crystal structure, it is likely that they reside in the infinite channels along the crystallographic c-axis and their contribution to the diffraction pattern was removed along with that of disordered solvent using SQUEEZE during refinement of the crystal structure.35
Close analysis of the bond lengths in cis-Fe(Ph2dhbq)(DMF)2 at 100 K and 250 K suggests the chain undergoes valence tautomerism, switching between (FeIII)(Ph2dhbq3˙−) at low temperatures and (FeII)(Ph2dhbq2−) at high temperatures (Fig. 3a). At 100 K, the average Fe–O bond length is 1.991(13) Å, consistent with high-spin FeIII. In addition, the average C–O bond length at this temperature is 1.305(13) Å, which is longer than the ∼1.24–1.28 Å typically observed in first-row transition metal structures containing delocalized R2dhbq2− ligands.16 In contrast, at 250 K, the average Fe–O bond lengths is 2.09(3) Å which is clearly lengthened relative to the Fe–O bond length observed at 100 K. While we were unable to observe a statistically significant change in the average C–O bond length between these temperatures due to lower data quality at 250 K, the average C–O bond length of 1.28(3) Å is in the range expected of R2dhbq2− ligands. Overall, the bond lengths observed at 100 K are more characteristic of FeIII metal centers bound to R2dhbq3˙− radical ligands, while the bond lengths at 250 K resemble previous structures containing FeII metal centers and R2dhbq2− ligands.6,16 While iron-based valence tautomerism is relatively rare,5 similar behavior was recently observed by Harris and coworkers in a structurally related Fe-based chain synthesized from chloranilic acid, [Me4N]2[Fe(Cl2dhbq)Cl2].6 The authors of this earlier work hypothesized that the FeIII state is stabilized by the presence of anionic, π-donating Cl− ligands. In contrast, our system contains only simple DMF ligands, suggesting that the relative arrangement of the R2dhbqn− ligands may be the critical factor, rather than the properties of any additional coordinating ligands.
Unlike cis-Fe(Ph2dhbq)(DMF)2, no structural changes indicative of valence tautomerism are observed in trans-Fe(Ph2dhbq)(DMA)2 (Fig. 3b). The Fe–O and C–O bond lengths are consistent with high-spin FeII metal centers and Ph2dhbq2− ligands at all temperatures. Single crystal X-ray diffraction data for trans-Fe(Ph2dhbq)(DMA)2 were collected at 250 K using a synchrotron X-ray source. The Fe–O bond lengths in trans-Fe(Ph2dhbq)(DMA)2 average 2.096(13) Å and C–O bond lengths average 1.264(3) Å.
Finally, the average bond lengths in oxidatively doped trans-Fe(Ph2dhbq)(DMA)2Br0.55 are in between those expected for FeII and FeIII (Fig. 3b).21 The average Fe–O bond length is 2.034(1) Å and the average C–O bond length is 1.277(3) Å at 250 K. Neither of these bond lengths change appreciably (by more than ±0.8%) when the crystal is cooled to 100 K (Table S3†). The C–O bond lengths match the bond lengths observed for cis-Fe(Ph2dhbq)(DMF)2 at 250 K, suggesting that despite the additional positive charge on the iron centers, the ligand is still in the Ph2dhbq2− oxidation state. No significant changes are observed in the νC–O region of the IR spectrum, consistent with metal-centered oxidation (Fig. S6 and S7†).
For cis-Fe(Ph2dhbq)(DMF)2, the Mössbauer spectrum at 240 K shows a single quadrupole doublet with an isomer shift (δ) of 1.12 mm s−1 and quadrupole splitting (|ΔEQ|) of 2.48 mm s−1, which are characteristic of high-spin, octahedral FeII (Fig. 4a). As the sample is cooled to 90 K, the spectrum displays parameters that are more consistent with high-spin, octahedral FeIII (δ = 0.64 mm s−1 and |ΔEQ| = 1.30 mm s−1).6,21 In between these two temperatures, it is possible to resolve three peaks which can be fit by two quadrupole doublets representing a mixture of FeII and FeIII species in roughly equal proportions (Table S5†). Together, these observations are fully consistent with valence tautomerization between (FeII)(Ph2dhbq2−) at high temperatures and (FeIII)(Ph2dhbq3˙−) at low temperatures, as suggested by our structural analysis.
For trans-Fe(Ph2dhbq)(DMA)2, the Mössbauer spectrum does not demonstrate any temperature-dependent change in the iron oxidation state. Instead, the data reveals a single quadrupole doublet with δ = 1.17 mm s−1 and |ΔEQ| = 1.43 mm s−1 at 240 K, which splits into two slightly offset quadrupole doublets at 90 K (δ = 1.22, 1.24 and |ΔEQ| = 1.64, 2.17 mm s−1) (Fig. 4b). All three doublets can unambiguously be assigned to high-spin, octahedral FeII. Again, these assignments agree with the bond lengths observed by single crystal X-ray diffraction.
Unlike the sharp quadrupole doublets observed in the previous two chains, the Mössbauer spectra of oxidatively doped trans-Fe(Ph2dhbq)(DMA)2Br0.55 shows an asymmetric quadrupole doublet at 240 K that broadens as the temperature is lowered (Fig. 4c). The spectrum at 240 K can be fit approximately by a single quadrupole doublet with δ = 0.79 mm s−1 and |ΔEQ| = 1.59 mm s−1. The presence of an averaged quadrupole doublet with an isomer shift that is in between what is typically observed for high-spin FeII and FeIII suggests that the charge is delocalized at this temperature, with electron transfer between the iron centers occurring at a rate that is faster than the Mössbauer time scale (∼10−7 s).21,36,37 This is in agreement with the crystallographically measured average Fe–O bond lengths, which also reside in between typical FeII and FeIII ranges (Fig. 3b). Decreasing the temperature to 90 K (the lowest temperature accessible on our instrument) slows down the rate of charge transfer, leading to increased line broadening.36 Overall, the Mössbauer data obtained for trans-Fe(Ph2dhbq)(DMA)2Br0.55 closely matches the behavior observed in previously reported Robin–Day class II/III mixed-valent iron complexes,36,37 where the extent of electronic coupling lies in an intermediate regime in between fully localized and delocalized.38 The presence of iron-based mixed-valency is further supported by UV-vis-NIR diffuse reflectance spectroscopy. The spectrum for trans-Fe(Ph2dhbq)(DMA)2Br0.55 displays a broad absorption band in the near-IR region with νmax ∼ 4200 cm−1 (Fig. S8†). We tentatively attribute this feature, which is not present in the undoped samples, to intervalence charge transfer (IVCT) between FeII and FeIII centers. The observed band has an asymmetric shape with a sharper cutoff at lower energies, consistent with a class II/III assignment.11,38 Similar IVCT bands in the near-IR region have been observed in other mixed-valence iron complexes and extended materials.36,39
Fig. 5 (a) Dc magnetic susceptibility (χMT) measurements for cis-Fe(Ph2dhbq)(DMF)2 (purple squares) and trans-Fe(Ph2dhbq)(DMA)2 (pink circles) under an applied field of 1000 G. The solid grey lines are best fits with Jcis = −230 cm−1 and Jtrans = +0.30 cm−1. Inset: zoom in of the valence tautomerization temperature (T1/2) of the cis chain. (b) Variable-temperature zero-field dc magnetic susceptibility (χM′) of cis-Fe(Ph2dhbq)(DMF)2 and trans-Fe(Ph2dhbq)(DMA)2. The grey lines are linear fits that yield correlation lengths of 12 for the cis chain and 4 for the trans chain. Data were collected at zero dc field with an applied ac field frequency of 1 Hz. See main text and ESI† for fitting details. |
The presence of the ligand radical, Ph2dhbq3˙−, at lower temperatures should give rise to strong metal–ligand magnetic coupling. Indeed, below the valence tautomerization transition of cis-Fe(Ph2dhbq)(DMF)2, the value of χMT exceeds the value expected for uncoupled high-spin FeIII (S = 5/2) with Ph2dhbq3˙− radicals (S = 1/2). This discrepancy likely indicates that magnetic coupling between the metal ions and ligand radicals is enabled as soon as valence tautomerization takes place. The best simulations of this data suggest a strong antiferromagnetic metal–radical exchange coupling constant of J = −230 ± 4 cm−1. Due to the differing numbers of unpaired electrons on the metal versus the ligand, the χMT increases with decreasing temperature, reaching a maximum at 13 K (χMTmax = 150 cm3 K mol−1). For strongly magnetically coupled 1D systems, this increase in susceptibility is indicative of an increasing intrachain spin correlation length (2ξ) between domain walls.40 Below 13 K, the magnetic behavior becomes governed by defect sites that limit the chain growth.41 As this chain growth becomes limited, χMT decreases precipitously, reaching 34 cm3 K mol−1 at 1.8 K.
To the best of our knowledge, the magnitude of the exchange coupling constant observed in cis-Fe(Ph2dhbq)(DMF)2 (J = −230 ± 4 cm−1) is the largest reported for any 1D chain, and is almost three times that observed for [Me4N]2[Fe(Cl2dhbq)Cl2] (J = −81 cm−1).6,42 The relative increase in the magnitude of the exchange coupling constant observed in chains constructed from Ph2dhbq versus Cl2dhbq may be attributable to improved orbital overlap between the metal ion and ligand radicals.43 While the average Fe–O bond lengths can be assigned to FeIII in both cases, the Fe–O bond length in cis-Fe(Ph2dhbq)(DMF)2 is 1.991(13) Å, which is shorter than the value of 2.027(2) Å observed in the Cl2dhbq chain,6 which should promote stronger magnetic coupling.
In contrast to the cis chain, dc magnetic susceptibility measurements of trans-Fe(Ph2dhbq)(DMA)2 exhibit no signatures of valence tautomerism. Instead, the susceptibility data reflect the expected behavior for a chain of ferromagnetically-coupled FeII centers (Fig. 5). Best fits of the data suggest a modest ferromagnetic exchange coupling, with J = +0.30 ± 0.05 cm−1. The trans-Fe(Ph2dhbq)(DMA)2Br0.55 chain exhibits similar magnetic behavior to trans-Fe(Ph2dhbq)(DMA)2, neither of which display the valence tautomerization behavior of the cis chain (Fig. S10 and S11†). Based on the Mössbauer and UV-vis-NIR analysis, the dc susceptibility of trans-Fe(Ph2dhbq)(DMA)2Br0.55 was fit assuming contributions from both S = 2 and S = 5/2 ions (see ESI† for further discussion). We note that the simulation of magnetic data in one-dimensional mixed valence systems is far from trivial.44 Nevertheless, best fits using this method yielded an exchange coupling value (J = +0.64 ± 0.01 cm−1) similar in magnitude to trans-Fe(Ph2dhbq)(DMA)2, yet still small relative to cis-Fe(Ph2dhbq)(DMF)2.
Magnetic hysteresis was not observed for any of the chains (Fig. S18–S20†), and charging curves reveal negligible interchain interactions owing to the absence of a low-field inflection in magnetization with increasing field (Fig. S21–S23†).
The in-phase ac susceptibility (χM′), measured as a function of temperature at a fixed ac frequency (the ‘zero-field dc susceptibility’), reveals the effective lengths of the chains in these solids (Fig. 5b and S10†) when plotted as ln(χM′T) versus 1/T.40,45 Here, the linear region was fit (see ESI† for details) to yield the domain wall energy (Δζ). These values are considerably larger for the cis chain versus the trans chains, with 37 cm−1 for cis-Fe(Ph2dhbq)(DMF)2 and 6 cm−1 for trans-Fe(Ph2dhbq)(DMA)2. The extension of the chains with increased exchange coupling follows expectations for 1D molecular magnetic materials.40,45 Likewise, trans-Fe(Ph2dhbq)(DMA)2Br0.55 exhibits a larger domain wall energy relative to the undoped chain (12 cm−1), which is expected given the increased exchange coupling.
The correlation length (n) was determined as described in the ESI† and gives the average number of contiguous metal ions coupled via the exchange interaction (J). As expected, for the cis chain we obtain n = 12 while trans-Fe(Ph2dhbq)(DMA)2 yields n = 4 (Table S12†). Interestingly, the doped trans-Fe(Ph2dhbq)(DMA)2Br0.55 yields the same correlation length as the undoped chain (n = 4), despite possessing larger exchange coupling as well as a larger domain wall energy. Overall, the magnetic characterization demonstrates clear differences between each of the three chains.
Two-probe pressed pellet conductivities for all chains were collected at ambient temperature in air immediately following removal of the samples from the glovebox. The electrical conductivities of cis-Fe(Ph2dhbq)(DMF)2 and trans-Fe(Ph2dhbq)(DMA)2 were found to be similar, with measured values of 2.31 ± 0.13 × 10−8 and 6.0 ± 1.5 × 10−9 S cm−1, respectively. In contrast, the electrical conductivity of mixed-valence trans-Fe(Ph2dhbq)(DMA)2Br0.55 was found to be 2.5 ± 1.3 × 10−6 S cm−1, nearly three orders of magnitude higher. We note that pressed pellet conductivity measurements are strongly affected by grain boundaries, particularly in highly anisotropic systems such as the 1D chains reported here. Nevertheless, the dramatic improvement suggests that the conductivity of these chains may be further improved through controlled oxidative doping. A more comprehensive study on the conductivity of these materials, including single crystal measurements and controlled oxidation and reduction experiments, will be the subject of future work.
Footnote |
† Electronic supplementary information (ESI) available. CCDC 2215980–2215984. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d2sc06392a |
This journal is © The Royal Society of Chemistry 2023 |