J.
Vallejo
a,
E.
Pardo
*a,
M.
Viciano-Chumillas
a,
I.
Castro
a,
P.
Amorós
b,
M.
Déniz
c,
C.
Ruiz-Pérez
c,
C.
Yuste-Vivas
a,
J.
Krzystek
d,
M.
Julve
a,
F.
Lloret
a and
J.
Cano
*ae
aInstitut de Ciència Molecular (ICMOL), Universitat de València, 46980 Paterna, València, Spain. E-mail: joan.cano@uv.es; Emilio.Pardo@uv.es
bInstitut de Ciència del Materials (ICMUV), Universitat de València, 46980 Paterna, València, Spain
cLaboratorio de Rayos X y Materiales Moleculares, Departamento de Física, Facultad de Ciencias (Sección Física), Universidad de La Laguna, 38201 Tenerife, Spain
dNational High Magnetic Field Laboratory (NHMFL), Florida State University, Tallahassee, Florida 32310, USA
eFundació General de la Universitat de València (FGUV), Spain
First published on 13th February 2017
A vast impact on molecular nanoscience can be achieved using simple transition metal complexes as dynamic chemical systems to perform specific and selective tasks under the control of an external stimulus that switches “ON” and “OFF” their electronic properties. While the interest in single-ion magnets (SIMs) lies in their potential applications in information storage and quantum computing, the switching of their slow magnetic relaxation associated with host–guest processes is insufficiently explored. Herein, we report a unique example of a mononuclear cobalt(II) complex in which geometrical constraints are the cause of easy and reversible water coordination and its release. As a result, a reversible and selective colour and SIM behaviour switch occurs between a “slow-relaxing” deep red anhydrous material (compound 1) and its “fast-relaxing” orange hydrated form (compound 2). The combination of this optical and magnetic switching in this new class of vapochromic and thermochromic SIMs offers fascinating possibilities for designing multifunctional molecular materials.
Molecules with a single slow-relaxing metal centre, referred to as single-ion magnets (SIMs), are the simplest model systems for fundamental research on magnetic anisotropy and magnetic relaxation dynamics, from both experimental and theoretical points of view.32,33 Being a non-cooperative property of an isolated molecule, the SIM behaviour depends on subtle and small structural variations in the metal coordination environment that directly determine the nature and magnitude of the local single-ion magnetic anisotropy, e.g. peripheral ligand modifications.16 SIMs are also particularly appealing because the slow magnetic relaxation effects can be controlled with structural modulation and an external stimulus (e.g., a magnetic field). These features together with their small size and easy handling and addressing to surfaces make these simple switchable molecules candidates for new magnetic devices in molecular electronics and spintronics.
The “entatic state” is a common term in biology, but is only rarely used in other disciplines such as bioinorganic or coordination chemistry.34–36 It is typically found in the activity of some enzymes, specifically those controlling redox processes. An “entatic state” is related to a constrained geometry halfway between the ideal geometries at the beginning and end of a reaction. The energetic requirements for a chemical reaction, i.e. the energy barriers, are minimized in these cases and the reaction is favoured. In other cases, the uptake and release of the substrate molecules are controlled by the entatic states. An illustrative example is provided by iron and manganese superoxide dismutases (SODs), where a distortion of the trigonal bipyramidal geometry (TBPY-5) of the active site provides easy entry for the charged superoxide (O2−) substrate. After evolution to hydrogen peroxide or molecular oxygen, the active site adopts a strongly distorted octahedral (OC-6) geometry that largely limits the space available to accommodate the generated neutral substrates (H2O2 and O2), which are minimized in these cases easily released because the Coulomb attraction between these substrates and the metal ion vanishes. Recently, the entatic state concept has also been indirectly used for “local sterically constrained environments” to improve new hybrid solar thermal fuels through the preservation of the robust cyclability and stability of the material.37
With the design and synthesis of multifunctional molecular materials for magnetic (or optical) switching in mind,38–41 we have developed a new strategy based on mononuclear cobalt(II) complexes with β,β′-dimethyl-substituted aromatic α,α′-diimines, such as 6,6′-dimethyl-2,2′-bipyridine (dmbpy) and 2,9-dimethyl-1,10-phenanthroline (dmphen). This class of sterically hindered, chelating ligands causes important structural distortions of the metal environment that can ultimately lead to a large single-ion magnetic anisotropy.17 Along this line, we report herein the synthesis, X-ray structures, thermal and sorption behaviours, spectroscopic and magnetic properties, and theoretical calculations of two new mononuclear cobalt(II) complexes with formulas [CoII(dmbpy)2](ClO4)2 (1) and [CoII(dmbpy)2(H2O)](ClO4)2 (2). This pair of anhydrous (1) and hydrated (2) forms exhibits a unique vapour- and thermochromic-induced SIM switching in the solid state, which is accompanied by reversible changes in their colour and magnetic properties. This switching behaviour is related to the coordination and release of a water molecule through vapour adsorption and heating, reminiscent of classical cobalt(II)-substituted molecular sieves (e.g. aluminophosphates or silicates like the well-known silica gel).42 The reversibility of this process is connected to a large distortion in both complexes, which is associated with entatic states induced by the steric constrains caused by the bulky chelating ligands.
Each cobalt atom in 1 is four-coordinate in a somewhat distorted tetrahedral (T-4) CoN4 environment, whereas it is five-coordinate in 2 adopting a CoN4Ow coordination sphere that is intermediate between a square pyramid (SPY-5) and a trigonal bipyramid (TBPY-5) (Fig. 1). The so-called Berry pseudo-rotation parameter, defined as τ = |α − β|/60, where α and β are the largest LCoL angles, is 0.44 for 2 (τ is equal to 1 and 0 for ideal TBP and SQP geometries).43 Nevertheless, this parameter is not the most adequate way to describe the highly distorted environment in 2, which is closer to a distorted asymmetric T path with a N(2)Co(1)N(4) angle (θ) equal to 165.8° and two different γ angles, 115.2 and 139.4° for OwCo(1)N(4) and OwCo(1)N(1), respectively. The tetrahedron in 1 is defined by four imine-nitrogen atoms from two almost perpendicularly placed chelating dmbpy ligands (Fig. 1), with the angle (α) between the Co(1)N(1)N(2) and Co(1)N(3)N(4) planes equal to 88.7(2)° and close to those found in a few very similar complexes (86.5° and 81.5°).44,45 As in these other cases, the bite angles from the two chelating dmbpy ligands [N(1)–Co(1)–N(2) = 83.5(2)° and N(3)–Co(1)–N(4) = 83.02(2)°] and the remaining free interbond angles [N–Co–N = 115.7(2)–131.4(2)°] in 1 are smaller and larger, respectively, than those for an ideal tetrahedron (109.5°). In contrast, it is very difficult to assign the axial and equatorial positions in 2 due to the highly distorted metal environment intermediate between a SQP and TBP (Fig. 1). Overall, a “plier” distortion in 1 leads to pseudo D2d symmetry,46 whereas the symmetry of 2 is halfway between C4v and D3h. This rare geometry for 2 is a consequence of the high geometrical restrictions that avoid a relaxing geometry when the water molecule is removed. As a result, continuous shape measurements47 (CShM) relative to the tetrahedron (T-4) and square plane (SP-4) on CoN4 coordination spheres are similar for 1 (5.2 and 25.7°) and for the dehydrated form of 2 (2deh, 9.4 and 18.3°). Moreover, quasi-coincident deviations from the pathway connecting the two ideal polyhedra were found for 1 and 2deh (23.9 and 22.4°, respectively).
This reversibility is linked to the highly distorted geometries of the metal environments in 1 and 2, which clearly favour the easy and reversible entry and release of a water molecule into the cobalt coordination sphere. As with the iron and manganese SODs, there are two entatic states that minimize the energy barrier for the entry and release of the water molecule, converting this into a reversible process. The structural modifications of the CoII complex accompanying the dehydration/rehydration process are quite small (Fig. 1 and S1, ESI†) – as was already established from the continuous shape measurements – which accounts for the modest value of the energy barriers for each step in the dehydration/rehydration process and explains its reversibility. Accordingly, in agreement with the experimental observation and from the free energies calculated in the gas phase, the hydrated form is the most stable with an energy barrier of just 0.18 kcal mol−1 for the hydration process (Fig. 3). Even if the energy barrier for the dehydration is larger (2.09 kcal mol−1), it is small enough to remove the coordinated water molecule with mild heating.
The successive heating and exposure to air of both a polycrystalline sample and a solution of 2 were followed up with various spectroscopic techniques corroborating the reversible uptake and release of a water molecule from the metal coordination sphere. The XRPD patterns of the dehydrated and rehydrated solid samples, in perfect agreement with the theoretical patterns of 1 and 2 (Fig. S2, ESI†), showed significant shifts of the peaks. The time-dependent XRPD patterns showed that, after a short induction period of approximately 10 min, half and complete conversions were reached in 20 and 60 min, respectively (Fig. 2), retaining the crystallinity of the sample in both phases after five dehydration/rehydration cycles.
The UV-Vis spectra in the solid state and in solution of the dehydrated and rehydrated samples matched with those of 1 and 2, supporting the assumption that the changes upon the reversible dehydration/rehydration process occur in the metal coordination sphere (Fig. 4 and S3a, ESI†). Both compounds showed the same electronic absorption spectra for the high energy region (200–300 nm), with two very intense and sharp intraligand bands. Each of the two compounds also exhibits a weaker and asymmetric band in the visible region although the intensity and energy of this band differs moderately, reflecting the changes in the metal environment. The deconvolution of the visible-region band in the solution spectra showed that it is composed of a set of electronic transitions (Table S3, ESI†). The fact that some of the electronic transitions in both compounds have the same energy suggests the presence of a mixture of the hydrated and dehydrated species. Nevertheless, the most intense electronic transitions have lower energies and are largely more intense in 1 (548, 572 and 588 nm) than in 2 (527 nm). These features agree with d–d transitions in the T-4 and SPY-5 ligand fields for 1 and 2, respectively. In this respect, it is well known that more intense and less energetic electronic transitions are provided by a less symmetric and weaker T-4 ligand field. Despite a slight shift in energy from the experimental data, TDDFT calculations in acetonitrile solution on both complexes supported the experimental observation, reproducing the relative positions and intensities of the electronic transitions in the discussed energy region. Besides, the TDDFT calculations predicted only two d–d electronic transitions for the hydrated and dehydrated species in the 500–600 nm region (Table S3, ESI†), as can be easily seen from the involved natural transition orbitals (Fig. 5). This is a confirmation that the experimental spectra in solution correspond to a mixture of the hydrated and dehydrated species with one of them being predominant, rather than to a single compound.
Fig. 5 Perspective views of the natural transition orbitals (NTOs) involved in the theoretical electronic excitations shown in Table S3† for 1 (left) and 2 (right). NTOs have been used because they provide an easier visualization of the electronic transitions compared to molecular orbitals. The isodensity surfaces correspond to a cut-off value of 0.03 e bohr−3. In this picture, one electron is promoted from the orbital at the left side to the other one at the right side. |
In our case, complex 1 belongs to the HFEPR-silent category, in qualitative agreement with the large negative D (−57.0 cm−1) and very low E/D (0.004) values estimated from NEVPT2 calculations. Unlike in 1, EPR spectra were observed for 2 in the 52–320 GHz frequency range which correspond to those expected for an S = 3/2 spin state with a positive D and sizeable rhombic ratio E/D. The pattern does not depend on frequency, thus the condition hν ≪ |D| (ν is the Larmor frequency) is satisfied for the whole frequency range. Fig. 6 shows one such spectrum at 56 GHz and low temperature (4.5 K) together with simulations that assume the magnitude of |D| obtained from NEVPT2 calculations. It is clear that the sign of D is positive. The spectrum at a much higher frequency (319 GHz) and its simulations are shown in Fig. S4 in the ESI.† An analysis of the multifrequency set of resonances (Fig. S5, ESI†) yields a moderate rhombicity (E/D = 0.125) of the zfs tensor together with an axial Landé g-tensor: g⊥ = 2.37 and g∥ = 2.16. These results agree with those obtained from the NEVPT2 calculations (E/D = 0.198), with the rhombicity E/D factor being slightly overestimated by the theory, as usually occurs in most cobalt(II) cases.
Because of the limited access to the HFEPR technique, successive dehydration/rehydration processes on a sample of 2 were followed using X-band EPR spectroscopy at 4.0 K. Whereas 2 shows a rhombic signal that could be analysed as a doublet effective spin (Seff = 1/2), with g1 = 2.10, g2 = 4.05, and g3 = 5.54 in solution (Fig. 4) and in the solid state (Fig. S3b, ESI†), no signal was recorded for 1 yet again. To be exact, despite all the care taken, and in accordance with the easy hydration of 1 observed from the XPRD studies and UV-Vis electronic spectroscopy, a zoom of the EPR spectra of 1 in solution confirmed partial hydration to afford 2 but in a very low amount (less than 2 per cent). These gi components for an effective spin agree well with those deduced from the g-tensor and E/D values found in the HFEPR study for S = 3/2 using D = +28 cm−1 obtained from NEVPT2 calculations (gx = 5.58, gy = 3.88, and gz = 2.07).49 Like for other spectroscopic techniques, the EPR spectra successively matched up with those of 1 and 2 after several dehydration/rehydration cycles, confirming the reversibility of the uptake/release of a water molecule.
The direct current (dc) magnetic susceptibility behaviour of 1 and 2 in the form of χMT versus T and M versus H/T plots (χM and M being the molar magnetic susceptibility and magnetisation, and T and H the absolute temperature and the applied magnetic field) revealed the occurrence of large single-ion magnetic anisotropies (Fig. 7 and S6, ESI†). The values of χMT at room temperature were 2.72 (1) and 2.73 cm3 mol−1 K (2). They are higher than those expected for an isotropic S = 3/2 spin moment [χMT = 1.875 cm3 mol−1 K with g = 2.0] but in agreement with those expected for one high-spin d7 CoII ion with unquenched orbital momentum. Upon cooling, χMT continuously decreases for both compounds to reach 2.06 (1) and 1.68 cm3 mol−1 K (2) at 5 K, values which suggest a different magnitude of the zfs in them. Isothermal magnetization curves in the form M versus H/T do not superimpose when a zfs or weak magnetic couplings are present, which occurred in 1 and 2 (Fig. S6, ESI†). However, if the zfs is very large, the two split Kramers doublets are well separated in energy, and a small or no temperature influence on the M vs. H/T magnetization isotherms is to be expected. In such a case, it is not possible to accurately extract the values of the zfs parameters from the magnetization curves, but the experiment allowed us to conclude that the zfs is significant in 1 and 2 and probably greater in the former compound.
The magnetic behaviour of OC-6 cobalt(II) complexes is usually induced by a first order SOC. Consequently, the use of the SOC Hamiltonian introduced through T-P isomorphism [eqn (S1) and (S2), ESI†] is the most appropriate option in the analysis of their magnetic properties. The α, λ, Δ and TIP parameters in this Hamiltonian are the orbital reduction factor, spin–orbit constant, energy level splitting of the orbital moment through the lowering of the molecular symmetry and temperature-independent paramagnetism, respectively.50 The least-squares fit of the magnetic susceptibility data of 2 from 5 to 300 K is: α = 1.32, λ = −118 cm−1, Δ = +805 cm−1, TIP = 246 × 10−6 cm3 mol−1 and F = 1.2 × 10−5 (F is the agreement factor defined as ∑[Pexp − Pcalcd]2/∑[Pexp]2, with P being the physical property under study). Among the six Kramers doublets coming from the coupling between the spin (S = 3/2) and orbit (L = 1) momenta in a distorted octahedron, the two lowest doublets are usually the unique populated states below 50 K; therefore, a phenomenological approach based on the zfs of an S = 3/2 state that considers only these doublets [as summarized in the spin Hamiltonian of eqn (S2) (ESI†)] is also justified at low temperatures. Above 50 K, a TIP term, however, should be introduced into this picture to take into account the depopulation of the remaining higher energy Kramers doublets. This approach including both terms is always efficient under 100 K but only occasionally up to room temperature. The experimental magnetic susceptibility of 2, in the temperature range 5–100 K, was analyzed with the VPMAG program.51 The least-squares fit gave: g⊥ = 2.41, g∥ = 2.07, D = +35.8 cm−1 and TIP = 4.37 × 10−6 cm3 mol−1 with F = 3.0 × 10−6. Although usually the thermal dependence of the magnetic susceptibility cannot unambiguously determine the sign of the D parameter, this was possible in 1 – and also in 2 – because of the large anisotropy in the g-factors. Both methods, T-P isomorphism and zfs, were appropriate and provided the same energy gap between the two lowest Kramers doublets (71.6 cm−1, see the Magnetic measurements section in the ESI†). It is remarkable that the axial and rhombic zfs parameters and the components of the Landé g-tensor obtained from magnetic susceptibility and HFEPR spectroscopy are similar to the theoretical ones [D = +27.6 cm−1, E/D = 0.198, g⊥ = 2.42 (gx = 2.50 and gy = 2.34) and g∥ = 2.12].
The electronic structure in T-4 Co(II) complexes is different than it is in those of OC-6 because the ground state is a non-degenerate orbital term and the SOC Hamiltonian in eqn (S1) (ESI†) is inadequate. Hence, the zfs comes from a 2nd order SOC and, as a result, it is generally smaller than in the OC-6 complexes. However, geometrical distortions can transform the ground state into a magnetic degenerated or quasi-degenerated orbital term increasing the SOC contributions to the zfs. In this sense, the magnetic susceptibility data of 1 from 5 to 300 K were perfectly simulated (F = 2.9 × 10−6) through the spin Hamiltonian described by eqn (S2)† and with the following best-fit parameters: g⊥ = 2.11, g∥ = 2.67, D = −50.6 cm−1, and TIP = 529 × 10−6 cm3 mol−1 (Fig. 7). Once more, they are similar to those derived through NEVPT2 calculations [D = −57.0 cm−1, E/D = 0.004, g⊥ = 2.10 (gx = 2.10 and gy = 2.10) and g∥ = 2.77].
A first order SOC cannot occur in an ideal T-4 environment because the ground state corresponds to a singlet orbital term. In such cases, the axial zfs is smaller than in the previous geometries. Some geometrical changes, however, such as a plier distortion, can convert the GS into a doublet orbital term, significantly increasing the axial zfs by means of a 1st order SOC (Fig. S7, ESI†). This is the case for 1 where the first excited quartet state (Q1) is close to the GS (1400 cm−1) and it is almost solely responsible for the very large negative axial zfs (Table S4, ESI†).
Despite 1 being distorted, the orientation of its zfs tensor is similar to that shown for the molecular axis derived from the d orbitals in an ideal T-4 symmetry (Fig. S8, ESI†). This is not the case for 2 which is highly distorted, and there is no coincidence even with the orientation observed in 1 (Fig. S8, ESI†). This fact should be the cause of the very different axial and rhombic zfs of 1 and 2.
Fig. 8 Temperature dependence of the alternating current (a) in- and (b) out-of-phase molar magnetic susceptibilities. (c) Arrhenius plot in the high temperature region of 1 under an applied static field of 1000 G with a ±5.0 G oscillating field in the frequency range of 100–10000 Hz (from black to green). The solid line in (c) is the best-fit curve through a model with two Orbach relaxation processes (Table S5†). |
Deconvolution of the thermal dependence of χ′′M allowed us to accurately calculate the maxima of χ′′M at a given frequency (τ = 1/2πν, τ being the relaxation time). The plot of ln(τ) vs. 1/T follows an Arrhenius law characteristic of a thermally activated mechanism with two different relaxation processes [τ = {(1/τ01) exp(Ea1/kBT) + (1/τ02) exp(Ea2/kBT)}−1] (solid lines in Fig. 8c and S11, ESI†). Each one of these two relaxation processes is predominant in different temperature regions; i.e. one of them relaxes faster in a particular 1/T region although both coexist at any temperature. What is responsible for the curvature in the Arrhenius plot is that there is a temperature region where both processes are competing without there being a clear winner. The data processing with both processes rivalling is adequate; however, to individually analyse each area using only one relaxation process is not acceptable because the presence of another relaxation misrepresents the predominant relaxation process. The calculated values of the first pre-exponential factor (τ0) and activation energy (Ea) (Table S5, ESI†) obtained this way are consistent with those found for other previously reported SIMs of OC-6 cobalt(II) complexes.52 Lower values of Ea were, however, found for the second relaxation process such as is observed for other CoII systems showing a double relaxation process.49 Like in that paper, analysis considering quantum tunnelling or Raman mechanisms for the second and faster relaxation processes was ambiguous, leading to disparate and even inconsistent values for the parameter that determines the direct, quantum tunnelling, and Raman mechanisms (see the Magnetic measurements section in the ESI†). A model with two Orbach processes, which could be related to parallel and transversal relaxation times, seems to account for the magnetic behaviour observed in 1.
Nevertheless, against the double relaxation process, small values of the α parameter were extracted from almost perfect semicircle Cole–Cole plots (Table S5 and Fig. S12, ESI†), discarding the spin-glass behaviour (α ≥ 0.4). This fact allowed us to conclude that two observed relaxation processes occur in a unique cobalt(II) complex. Similarly to the majority of cobalt(II) complexes with first order SOC, the energy barrier calculated from the experimental D value [Ea = −D(S2 − 1/4) = −2D = 101.2 cm−1] is much larger than the experimental one derived from the Arrhenius behaviour.52
In summary, structural changes in 1 and 2, associated with a reversible dehydration/rehydration process, are responsible for the reversal of the sign of D and the disappearance of the SIM behaviour. It was originally believed that only negative D values permitted slow magnetic relaxation effects and such effects would not be possible for positive D values.53 The dynamic magnetic behaviours exhibited by 1 and 2 match this logical assumption. However, since about a dozen examples of Co(II) SIMs with large and positive values of D can already be found in the literature, explaining the magnetic switch between 1 and 2 taking into account only the D sign is naïve.17,18,23,49,54,55 On the other hand, the magnetic relaxation through a high energy barrier induced by an axial zfs could be assisted by extra energy investment from phonons emitted by the solid network from a vibrational relaxation.54 Hence, the synthesis and investigation of novel cobalt(II) complexes are essential to understanding this singular magnetic behaviour, although an improvement of the current theoretical models or the development of new ones that answer the open questions will be a landmark in the field.
Footnote |
† Electronic supplementary information (ESI) available: Preparation methods and physical characterization data. Crystallographic refinement and computational details. Additional figures (Fig. S1–S12) and tables (Tables S1–S5). CCDC 952077 and 938463. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/c6sc05188j |
This journal is © The Royal Society of Chemistry 2017 |