Molecular electronic spin qubits from a spin-frustrated trinuclear copper complex.

The trinuclear copper(ii) complex [Cu3(saltag)(py)6]ClO4 (H5saltag = tris(2-hydroxybenzylidene)triaminoguanidine) was synthesized and characterized by experimental as well as theoretical methods. This complex exhibits a strong antiferromagnetic coupling (J = -298 cm-1) between the copper(ii) ions, mediated by the N-N diazine bridges of the tritopic ligand, leading to a spin-frustrated system. This compound shows a T2 coherence time of 340 ns in frozen pyridine solution, which extends to 591 ns by changing the solvent to pyridine-d5. Hence, the presented compound is a promising candidate as a building block for molecular spintronics.

Quantum technology is one of the challenging fields of the current century due to potential major advances in information technology. 1 Qubits are the basic building blocks of quantum computers and in a straightforward approach can be realized by nuclear or electronic systems. 2 Although nuclear spins generally show longer coherence times than electronic spins, the latter have been shown to be suitable in terms of molecular magnets. 3 In fact, electronic spins related to magnetic molecules offer important advantages, as their structural and electronic properties can be tuned by chemical means and are easier to manipulate by employing electric fields. 4 Along this line, homotrinuclear complexes of paramagnetic metal ions have been attracting interest for several reasons. For magnetically anisotropic metal ions with large spins this includes cases which have been investigated with regard to their spin interaction in the triangle, 5 tuning single-molecule magnets 6 and single-molecule toroic behavior. 7 In contrast to ferromagnetically coupled molecular triangles, 8 an antiferromagnetic coupling causes situations which Kahn described as ranging from ''competing spin interactions'' to ''degenerate frustration''. 9 For molecules with an equilateral triangle topology and strict C 3 symmetric arrangement, antiferromagnetic interactions lead to specific properties. In the case of ions with integer spin like Ni(II) this results in a diamagnetic ground state (S = 0), 10 whereas systems based on non-integer spin such as Cu(II) ions give rise to the magnetic phenomenon of two degenerate spin ground states, the so-called spin frustration. 9,11 This situation in molecular systems is predicted to give rise to spin-electric coupling that allows manipulating the molecular spin states. 12 For practical implementation a large antiferromagnetic isotropic exchange interaction to provide an energetically well separated chiral spin ground state and small antisymmetric exchange interactions are favorable molecular properties. 13 However, this combination is hardly present in most complexes reported so far, 14 although such systems are highly demanded for quantum computing and spintronic applications. 15 In this work, we present the synthesis and characterization of the trinuclear Cu(II) complex [Cu 3 (saltag)(py) 6 ]ClO 4 (Cu 3 saltag), which is based on the central tritopic triaminoguanidine derived ligand (H 5 saltag = tris(2-hydroxybenzylidene)triaminoguanidine). The compound was obtained by a stepwise synthetic approach utilizing three equivalents of Cu(II) perchlorate, one equivalent of the ligand hydrochloride (H 5 saltagÁHCl), six equivalents of NEt 3 as a base, and finally an excess of pyridine to saturate the coordination sphere of the copper centers (see the ESI †). Cu 3 saltag crystallizes as Cu 3 saltagÁ1/2([Hpy]ClO 4 )ÁH 2 O in the trigonal space group P31c (Table S1, ESI †). The unit cell contains two crystallographically independent trinuclear monocationic complexes [Cu 3 (saltag)(py) 6 ] + ( Fig. 1 and Fig. S1, ESI †), one pyridinium cation, three perchlorate anions, and two co-crystallized water molecules. Each complex cation is centered on one of the two distinct crystallographic 3-fold axes, so that within the independent cations all Cu centers are equivalent. The general arrangement of the individual components along the two crystallographic 3-fold axes is illustrated in Fig. S2 (ESI †). The phase purity of the bulk material was confirmed by X-ray powder diffraction ( The distorted trigonal bipyramidal coordination environment is confirmed by continuous shape measures, which give the best fit for a trigonal bipyramid (Cu1: S(D 3h ) = 1.694; Cu2: S(D 3h ) = 1.705). 16 The three perchlorate anions present in the crystal structure, which counterbalance the charge of the monocationic complex molecules and the co-crystallized pyridinium cation, are located with the chlorine atom on special positions on either of the crystallographic 3-fold axes with one of the Cl-O bonds aligned along these axes. For all three perchlorate anions this leads to ClO 4 Á Á ÁCn1 ionic interactions with the central carbon atom of the triaminoguanidine ligand ( Fig. S5 and Table S3, ESI †). In the case of the complex molecule with Cu1 two perchlorate anions are embedded in the two propeller shaped cavities formed by the pyridine co-ligands (Fig. S6, ESI †). On the other hand, for the second complex molecule (Cu2) only one perchlorate anion is present (Fig. S7, ESI †), with the remaining void on the opposite site filled with co-crystallized water molecules (Fig. S2, ESI †). DFT studies for these arrangements reveal remarkably high electrostatic interaction energies with the cationic trinuclear complexes ([Cu 3 (saltag)(py) 6 ] + Á Á ÁClO 4 À binding energies: Cu1, 69 and 63; Cu2, 60 kcal mol À1 ) (Table S4, ESI †).
These attractive energies contain a small amount of pyÁ Á ÁClO 4 À interaction (Cu1: 17 and 14; Cu2: 13 kcal mol À1 ; see Table S5 and Fig. S8, ESI †). Interestingly, structural optimization studies with DFT showed that the perchlorate anions are most probably the reason for stabilizing the copper(II) ions in a trigonal bipyramidal coordination sphere in the crystalline material (Fig. S9, ESI †). This is confirmed by the optimized structure obtained without perchlorate anions being present, for which the trinuclear complex exhibits a somewhat more square pyramidal coordination environment (Fig. S9, ESI †). Additionally, also p-p interactions between the pyridine co-ligands of the crystallographically independent complex molecules are observed (Fig. S5, ESI †). Hence, each trinuclear cationic complex is linked to three neighboring complexes of the other crystallographic type, forming an infinite twodimensional double layer within the crystallographic ab plane ( Fig. S10 and S11, ESI †).
Temperature-dependent magnetic susceptibility measurements for Cu 3 saltag (Fig. 2) revealed a room temperature w M T value of 0.66 cm 3 K mol À1 , which is significantly lower than the spin-only value expected for three independent Cu(II) ions with S = 1/2 (1.13 cm 3 K mol À1 ). Upon lowering the temperature, the w M T value further decreases and reaches a plateau at roughly 80 K. This behavior is indicative of a very strong antiferromagnetic intramolecular magnetic exchange between the three copper(II) centers. At temperatures below 10 K, the w M T value again decreases, which indicates the presence of very weak intermolecular interactions. The data were fitted to the spin Hamiltonian given in eqn (1) using the software PHI 17 (see the ESI †). The spin coupling scheme in combination with the resulting magnetic states is depicted in Fig. S12 (ESI †). The best fit was obtained with the parameters given in Table S6 (ESI †), revealing a strong antiferromagnetic exchange ( J = À297.8(5) cm À1 ) and weak antiferromagnetic intermolecular interactions (zJ = À0.048(3) cm À1 ). Moreover, including antisymmetric exchange in the Hamiltonian does not improve the quality of the obtained fit (Table S7 and Fig. S13, ESI †). In fact, intermolecular interactions better describe the low temperature  behavior. The upper limit for a possible contribution from antisymmetric exchange can be estimated to be |G z | o 1 cm À1 (|G z /J| o 0.003). Although the effective value for G z is expected to be significantly lower, even the given upper limit is already well below the usually observed range. 14,18 Broken-symmetry DFT (BS-DFT) calculations were performed to get further insights into the magnetic exchange (Fig. S14, ESI †). The BS-DFT results obtained for pairwise interactions in zinc-replaced model structures verify the strong antiferromagnetic exchange in both complex cations of Cu 3 saltag ( J Cu1Cu1 = À281 cm À1 ; J Cu2Cu2 = À274 cm À1 ; J av = À278 cm À1 ; see Table S8, ESI †) and are in good agreement with the experimental value of À297.8(5) cm À1 . Alternatively, the exchange coupling can also be directly derived from the energy difference DE = E Q À E D = 3/2J, assuming that both the quartet and doublet states are correctly represented by the single determinate solutions for the HS and BS states (DE = E HS À E BS ), respectively, which, however, slightly underestimates the experimental values (Cu1: À255; Cu2: À249 cm À1 ; cf. Table S9, ESI †). 19 The strong antiferromagnetic exchange within the complex cation is mediated by the N-N diazine bridges, which are formed upon deprotonation of the ligand. Spin-density plots of the high-spin and broken-symmetry states as well as the magnetic orbitals obtained by orbital transformation 20 are depicted in Fig. S15-S18 (ESI †). The donor atoms of the [N 2 O] binding pockets of the saltag 5À ligand show a significant spin polarization, whereas the two pyridine nitrogen atoms of the co-ligands are significantly less spin polarized.
Continuous wave X-band ESR measurements at 6 K of a powdered crystalline sample of Cu 3 saltag revealed axial type anisotropy for the exclusively populated magnetic ground state (see Fig. S19 and Table S10, ESI †). A suitable fit for the data is obtained using a g factor with the two principal components g > = 2.04 and g 8 = 2.19 together with a strain parameter for g 8 of 0.20. The latter is indicative of a distribution of individual g 8 values due to minor differences in the coordination environment of the copper(II) ions along their magnetic z axis. This observation appears reasonable based on the presence of two different complex cations in the structure and the rigid ligand plane provided by the saltag 5À ligand, the latter being in contrast to the rather mobile monodentate pyridine co-ligands. Similarly, for X-band ESR spectra of frozen solutions of Cu 3 saltag in pyridine (pyridine-d 5 ) at 3 K (see Fig. S19, ESI †), an axial signal was observed, which can be fitted with a rather similar g > of 2.05 (2.04) and but requires a slightly larger g 8 of 2.29 (2.30) as compared to the solid state spectra. However, also for the spectra of the frozen solutions an identical strain parameter g 8 of 0.20 (0.25) was necessary to reproduce the experimental results in the fit. Nevertheless, upon dissolving Cu 3 saltag in pyridine (pyridine-d 5 ) we expect the intermolecular p-p as well as the ionic perchlorate interactions to be strongly reduced. Therefore, it is tempting to attribute these changes in interactions to the observed difference in the g 8 values between the bulk material and pyridine solutions. The latter is indicative of a slight variation of the copper(II) coordination environments between both cases.
Ab initio calculations were carried out for both independent copper(II) centers (Cu1 and Cu2) to investigate their magnetic anisotropy ( Fig. S20 and Table S11, ESI †). For both centers a 2 D multiplet with a 2 A 0 ground state is obtained as expected for copper(II) in a trigonal bipyramidal coordination sphere (for relative CASSCF and CASPT2 energies see Table S12 (ESI †); for spin-orbit coupled energies see Table S13, ESI †). Since the 2 A 0 ground state does not show a first-order spin-orbit contribution, any magnetic anisotropy is solely based on the second-order spin-orbit contributions of the excited 2 E 0 and 2 E 00 states. The magnetic anisotropy of both copper(II) centers shows a strong rhombic distortion (g x a g y a g z , see Table S14 and Fig. S21 and S22, ESI †). The hard axis of magnetization g z (Cu1/Cu2: 2.019/2.019) at both metal centers coincides with the Cun-On1 bond vector (angle Cu1/Cu2: 1.31/1.91), whereas the intermediate axis g y (Cu1/Cu2: 2.255/2.250) is nearly parallel to the Cun-Nn2 bond vector (angle Cu1/Cu2: 9.01/10.31) and the easy axis of magnetization g x is parallel to the molecular C 3 axis (Cu1/Cu2: 2.413/2.413). Utilizing the single-ion anisotropies to simulate the magnetic susceptibility of the trinuclear cationic complexes (see the ESI †) shows good agreement with the experiment (see Fig. S23, ESI †), but slightly too large w M T values. Such a shift in w M T values is based on an overestimation in the corresponding g factors and is a known phenomenon. 22 The calculated g factors for the molecular framework are given in Table S15 (ESI †). In fact, the situation becomes somewhat different when going from the single-ion anisotropy to the anisotropy of the trinuclear framework. For the latter, an axial system (Cu1: g > = 2.158 and g 8 = 2.425; Cu2: g > = 2.158 and g 8 = 2.423) is confirmed, where the g 8 axis coincides with the C 3 axis of the molecular framework.
For perspective applications, spin relaxation times are a crucial feature of the system. They can usually be enhanced via spatial separation of the spin centers. 23 Therefore, we performed spin-relaxation time measurements via echo-detected X-band ESR experiments in frozen pyridine and pyridine-d 5 solutions of Cu 3 saltag (see the ESI †). The data were fitted according to eqn (2) to determine the T 2 coherence times depending on the used solvent. In pyridine a T 2 time of 340 ns was obtained (see Fig. S24, ESI †), whereas in pyridine-d 5 this coherence time is significantly longer (591 ns; see Fig. 3). Fitting the echo decay gives a stretch parameter of x = 0.73, suggesting that the coherence time is limited by spectral diffusion. 24 This could be due to ligand rotation and libration processes related to the flexible pyridine ligands. 25 Such processes could possibly be suppressed by replacing them with sterically demanding ligands, further increasing the coherence time.
Y(t) = Y(0) exp(À(T 2 /2t) x ) (2) In conclusion, the spin-frustrated trinuclear copper(II) complex Cu 3 saltag with large antiferromagnetic exchange was synthesized. The magnetic properties were characterized by experimental and theoretical methods. Contributions from antisymmetric exchange are estimated to be well below |G z /J| o 0.003, making the system favorable for future practical implementation. Moreover, the compound shows a T 2 coherence time of 591 ns in a frozen pyridine-d 5 solution. These findings make Cu 3 saltag a promising candidate for further investigations of its applicability as a qubit and of spin-electric couplings.
We acknowledge DESY (Hamburg, Germany), a member of the Helmholtz Association HGF, for the provision of experimental facilities. Parts of this research were carried out at beamline P11, PETRA III. The authors would like to acknowledge networking support by the MolSpin COST Action CA15128. JM would like to thank EPSRC and Magdalen College, Oxford, for a doctoral studentship.

Conflicts of interest
There are no conflicts to declare. Fig. 3 Decay of the Hahn-echo intensity for Cu 3 saltag in a frozen pyridine-d 5 solution measured at 3 K. The solid red line shows the best fit according to eqn (2). The fit assumes that the oscillation in the signal is due to the ESEEM effect dominated by a single harmonic at the deuterium Zeeman frequency. 7