K.
Keller
a,
I.
Ritsch
a,
H.
Hintz
b,
M.
Hülsmann
b,
M.
Qi
b,
F. D.
Breitgoff
a,
D.
Klose
a,
Y.
Polyhach
*a,
M.
Yulikov
*a,
A.
Godt
*b and
G.
Jeschke
a
aLaboratory of Physical Chemistry, Department of Chemistry and Applied Biosciences, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland. E-mail: yevhen.polyhach@phys.chem.ethz.ch; maxim.yulikov@phys.chem.ethz.ch
bFaculty of Chemistry and Center for Molecular Materials (CM2), Bielefeld University, Universitätsstraße 25, 33615 Bielefeld, Germany. E-mail: godt@uni-bielefeld.de
First published on 14th September 2020
Determining distributed exchange couplings is important for understanding the properties of synthetic magnetic molecules. Such distributions can be determined from pulsed dipolar spectroscopy (PDS) data, but this is challenging due to the similar influence of both exchange and dipolar couplings on such data. In this work we introduce two models that aim to identify these two contributions to the spin–spin couplings from frequency-domain PDS data of shape-persistent molecules having either two Cu(II) ions, or a Cu(II) ion and a nitroxide radical as the paramagnetic moieties. The first model assumes correlated Lorentzian or Gaussian exchange and dipole–dipole coupling distributions whose parameters are the model's unknowns. The second model relies on prior knowledge of the distance distribution and by performing Tikhonov regularization along the exchange coupling dimension yields the latter distribution model-free. Both models were able to differentiate between the absence and the presence of exchange interaction, to determine the coupling regime (ferro- or antiferromagnetic) and to estimate the distribution shape. In contrast, calculations within the exchange resilient model of the neural network analysis implemented in DeerAnalysis2018 were not able for our data to identify exchange couplings and return correct distance distributions. However, the generic model was able to identify and separate the strongly curved intermolecular background in the relaxation-induced dipolar modulation enhancement (RIDME) experiments. Our analysis revealed that in such systems exchange coupling may be present up to at least 3.3 nm in π-conjugated systems involving Cu(II)–PyMTA, while it is negligible for distances r ≥ 4.5 nm between Cu(II) ions and r ≥ 3.8 nm between a Cu(II) ion and an unpaired electron of a nitroxide radical. Disruption of the π-conjugation between the ligand of the Cu(II) complex and the nitroxide leads to negligible exchange coupling at distances r ≥ 2.6 nm in the corresponding [Cu(II)–TAHA]–nitroxide ruler. Overall, for cases with known distance distributions, the presented analysis techniques allow to determine distributions of exchange couplings from PDS data.
Paramagnetic Cu(II) ions are interesting for determining distances between specific sites, either naturally occurring metal binding sites of metalloproteins17–20 or selected sites of interest, where Cu(II) spin tags21–23 can be attached by site-directed spin labelling.24–27 Significant progress in performing and analysing Cu(II)-based PDS experiments has been achieved over the past few years.28–33 Exchange interactions are often observed for Cu(II) ions that are bound to conjugated π-systems,34–39 or in strongly coupled systems, such as molecular magnets.7,40,41 The exchange coupling may give additional information on the system of interest, yet, the dipolar and exchange couplings can become difficult to disentangle.3,42–44
Approaches to determine exchange coupling strengths based on the analysis of the EPR lineshape have been reported.45–48 In cases, where the isotropic tumbling limit is achieved in liquid solution, the separation of isotropic exchange and anisotropic dipole–dipole coupling is straightforward because anisotropic contributions are averaged to zero as for example demonstrated for Cu(II) compounds in ref. 34 and 35. Note however that fast dynamics of the molecules in solution averages distribution of the exchange couplings for different possible conformations of the magnetic molecule as well as for dynamical shape distortions within each conformer, so that only the mean value can be accessed in this regime.
In the other limiting case, i.e. in the rigid limit, in which molecules are randomly oriented and immobilized in a glassy solution, both exchange and dipole–dipole couplings influence the shape of the EPR spectrum. If the half-field transition – which is independent from isotropic exchange coupling – can be detected, it may be a convenient measure to separate the exchange and dipolar contributions.42,49 Generally, separating and accurately determining both couplings in EPR lineshape analysis is a complex task. All parameters must be explicitly taken into account which influence the resonance frequencies of the coupled spins, such as principal values and relative orientations of all relevant magnetic tensors (g-tensor, hyperfine coupling tensor etc.). The problem is complicated further if any of the parameters in question are distributed. In such cases, independent knowledge of the molecular structure is helpful as it can, for instance, provide information on the dipole–dipole couplings, and leave only the distribution of exchange couplings as an unknown.39 Finally, if the line width exceeds the spin–spin coupling strength, the latter remains unresolved in EPR spectra in both the fast tumbling and the rigid limit regimes, thus making the EPR lineshape analysis approach inapplicable for determining spin–spin interactions.
For many spin systems the spin–spin couplings can be separated well from other interactions, such as g and hyperfine anisotropies by PDS. Hence, knowledge of the principal values of the corresponding magnetic tensors as well as their mutual orientations is no longer strictly required to be able to analyze the spin–spin couplings. Another strength of PDS experiments is that they can provide detailed information on the distribution of spin–spin couplings. For instance, in biological applications this is used to obtain site-to-site distance distributions from precisely measured dipole–dipole coupling distributions.1,4–6 In synthetic magnetic molecules it is also important to understand how broadly the spin–spin couplings are distributed. Such magnetic molecules are often designed to have a fixed geometry, in order to obtain defined magnetic properties. This implies a rather narrow distribution of spin–spin distances, and, hence, dipolar couplings. However, since the exchange interaction strongly depends on the degree of orbital overlap, it can vary substantially even upon small distortions of the molecule structure. Thus, while for such magnetic molecules the assumption of narrow spin–spin distance distribution is often valid, the distribution of exchange couplings could still be quite broad.36,37 Furthermore, stiffness of magnetic molecules facilitates the so called orientation-selection effect,50 where, depending on the excitation scheme of a particular experiment (e.g. microwave pulse positions, frequencies, excitation bandwidths and shape of the pulses etc.) different parts of the overall spin–spin coupling distribution can be attenuated or amplified reflecting the molecular geometry. Thus, orientation selection data are richer in information. However, processing of such data is much more complicated as it requires to consider orientations of hyperfine and g tensors.
In some practical cases, e.g. molecular magnets or certain types of biradicals used for DNP, the dipolar and exchange coupling values may be too large to be detected with PDS. However, measurements of weaker spin–spin couplings on analogous magnetic molecules appear useful for understanding the principles of the magnetic molecule design. Spin–spin coupling can be controlled by varying the distance between the two paramagnetic moieties in the molecule.51
PDS experiments are usually performed in frozen glassy solvents and, thus, they can reveal distributions of exchange couplings for different conformations of the magnetic molecules, representing a snapshot of the dynamic processes in solution corresponding to the glass transition temperature.52 Such a determination in PDS experiments is, however, still quite challenging due to the mentioned combination of distributions of dipole–dipole and exchange couplings along with orientation selection effects. Quite substantial work was done to demonstrate how orientation selection effects in the PDS experiments can be reduced by averaging over a set of pump and detection positions in Cu(II) and nitroxide EPR spectra.29,33,53–56 Such (partially) orientation-averaged PDS data are good test examples to study whether and how accurately the distributions of exchange and dipolar couplings can be disentangled and simultaneously determined. In the most general case, where both dipolar and exchange couplings are broadly distributed and correlated to each other, information content of the data may be too small for unambiguous analysis. However, many practical examples allow for simplifications, which then make the overall analysis more reliable.
Any PDS approach relies on the detection of time traces that oscillate with frequencies determined by the spin–spin coupling distribution. It is limited by the transverse relaxation time and thus exhibits higher resolution than the inhomogeneously broadened EPR spectrum. Fourier transform of such time traces yields an exchange-distorted Pake-like dipolar frequency pattern. If its singularities are well resolved, their dependence on exchange and dipolar interaction57,58 can be used to extract the individual couplings.59 Extraction of spin–spin couplings has also been performed by simulation of the experimental form factor and taking the spin–spin couplings as well as excited orientations into account.36,37,60,61 In some of the cases singularities of the Pake-like patterns were taken as initial guess for the spin–spin couplings in the simulation of the full time traces. Similar simulations, with an analytical expression of the spin–spin couplings in time domain, have been presented for time traces with negligible orientation-selection and presuming that the spin–spin distance is known.38
We have recently encountered a curious case of exchange interaction in shape-persistent molecules featuring Cu(II)–Cu(II) or Cu(II)–nitroxide radical spin pairs, wherein the strength of spin–spin exchange interaction depends on the types of paramagnetic moieties involved. In case of the Cu(II) ion complex, it also depends on the type of ligand. In this contribution, we thus address the problem of determining exchange coupling and distance distributions with four different shape-persistent molecules involving Cu(II) ions (see Fig. 1). We apply the chelators PyMTA or TAHA and consider Cu(II)–Cu(II) as well as Cu(II)–nitroxide rulers with both types of chelators. We attempt to extract electron spin–spin couplings from orientation-averaged frequency-domain PDS data using three different approaches. The paper is structured as follows: first, we describe the studied spin systems, their synthesis and their characterization by CW EPR and DFT calculations. Second, we present UV/VIS absorption spectra for the ruler compounds and the corresponding ruler precursors before Cu(II)-loading. Next, we describe the time traces obtained by Relaxation-Induced Dipolar Modulation Enhancement (RIDME)62,63 and Double Electron Electron Resonance (DEER)64,65 experiments on these compounds. Fourth, on the example of the [Cu–PyMTA]–[Cu–PyMTA] ruler, we test the applicability of our analysis routines and of analysis by neural networks as available in DeerAnaylsis2018.66 Finally, we apply the best performing data analysis to all other studied systems.
Fig. 1 Structural formulae of the Cu(II)–Cu(II) and Cu(II)–nitroxide rulers and the corresponding complexes Cu(II)–PyMTA and Cu(II)–TAHA used in this study. The effective spacer is marked violet. The displayed coordination of the Cu(II) ions shows the maximum potential coordination bonds. See Section 4.1 for the question about the numbers of coordinating atoms. The Cu–TAHA moiety is partially protonated with about 1.9 protons per moiety.56 Due to the structural similarity of PyMTA to diethylenetriaminepentaacetic acid (DTPA) the Cu–PyMTA moiety is expected to be completely deprotonated at pH 7 as reported for the complex Cu(II)–DTPA.75 |
(1) |
The first two terms describe the orientation-dependent electron Zeeman interaction for the A and B spins. The next two terms describe the hyperfine interaction between the electron spins and non-zero nuclear spins. The last term describes electron–electron spin interactions Ĥee. In the weak coupling approximation, which applies to an overwhelming majority of spin pairs because of the strong anisotropy of the g and Cu(II) hyperfine tensors and relatively small spin–spin coupling strengths, this term can be simplified to:
Ĥee = JŜA,zŜB,z + ωddŜA,zŜB,z, | (2) |
The dipolar coupling ωdd is dependent on the angle θ between the spin–spin vector and the external magnetic field
(3) |
It is also dependent on the effective g-value of each spin at a given orientation (ϕ, θ)
(4) |
(5) |
The dependency of the dipolar coupling constant on the angle θ between the spin–spin axis and the external magnetic field can lead to aforementioned orientation selectivity effects in experimental dipolar evolution data if not all orientations are excited uniformly.50 Such data can be described by introducing a geometry factor ξ(θ), which accounts for excited fractions of orientations in the signal evolution and depends on sample geometry and the applied microwave pulses setup:58
(6) |
In cases where all orientations can be excited uniformly, or the time-domain signal can be completely averaged over all orientations, the signal in frequency domain corresponds to the Pake-like pattern for a given spin–spin distance rAB and exchange coupling constant J.57,58 In fact, orientation selectivity effects are often observed as distortions of the Pake-like pattern and they may complicate data analysis significantly as knowledge of the molecular geometry is required in such cases. The strong anisotropy of Cu(II) spectra and insufficient excitation bandwidth of microwave pulses often lead to orientation selection effects.36,54,67–70 Such effects can be included into data analysis and give additional information on the geometry of the system.36,69,70 However, to simplify the problem of disentangling exchange and dipole–dipole couplings, we aimed to acquire data with negligible orientation selection by averaging over several detection fields.33,53–55
Extraction of the spin–spin interactions between two paramagnetic centres with broad EPR spectra posed the need to develop PDS experiments that can cover a large bandwidth. In recent years two main techniques developed in this respect: the use of ultrawideband (UWB) inversion pulses in combination with the DEER experiment71 or inversion by stochastic longitudinal relaxation events instead of an inversion pulse in the RIDME sequence.62,63,72
The (UWB-)DEER and RIDME pulse sequences are shown in Fig. 2. In both experiments electron spin coherence is created by the first π/2 pulse of the initial preparation period (π/2–τ–π sequence). After this initial preparation period, a modulation to the local field at the position of the observed spin A is induced by inversion of a second spin B. This causes a shift of the resonance frequency of the A spin. The change in local field is proportional to the spin–spin coupling, which in turn depends on the inverse cube of the A–B distance and on orientation for the dipole–dipole coupling, as well as on the exchange coupling constant. Thus, depending on the inversion efficiency pinv of the coupled B spin, a fraction of the electron spin echo signal will oscillate with the sum of the dipolar frequency ωdd and the exchange coupling constant J. This fraction λ is often called modulation depth. The evolution of the A spin with changed frequency results in a phase gain of ±(ωdd + J)t, which is then detected using a final refocusing pulse. An important difference between the two pulse sequences is the nature of the inversion of the B spins.
Fig. 2 (a) Four-pulse DEER65 and (b) five-pulse RIDME63 sequences. In (a) the dipolar evolution signal is detected as the intensity variation of the refocused primary echo with respect to the shift of the time position of the pump pulse (blue). In (b) the so-called refocused virtual echo is detected (see ref. 63) and its intensity varies with the change of the position of the mixing block, formed by the last two π/2 pulses. |
In the DEER experiment the B spins are inverted by a pump pulse applied at a second frequency (or with a frequency-modulated pulse, whose excitation band does not overlap with the detection band) and the inversion efficiency is given by the excitation profile of the pump pulse. In the absence of orientation selection, the modulation depth of the DEER experiment is equal to the inversion efficiency of the pump pulse.
In the RIDME experiment the magnetization is stored in the direction of the external field in the form of a polarization grating for a certain time Tmix. During this mixing time Tmix longitudinal relaxation events lead to stochastic inversion of B spins. After the mixing time, the magnetization is transferred back into the transverse plane, where it is eventually detected. The modulation depth is thus determined by the time allowed for such B spin flips. To first approximation, the modulation builds up according to the probability of an odd number of spin flips during the mixing block:73,74
(7) |
Q-band data were acquired either on a commercial Bruker Elexsys E580 Q-band spectrometer or on a home-built spectrometer based on a Keysight Arbitrary Waveform Generator (AWG)76 with a high-power Q-band extension.77 A helium flow cryostat (ER 4118 CF, Oxford Instruments) was used to adjust the measurement temperature to 20 K. The commercial spectrometer, equipped with a home-built cavity operating at about 34.5 GHz,78,79 was used for the metal–metal RIDME experiments. The home-built spectrometer was combined with either the same resonator for the orthogonal-RIDME experiments or with a broadband resonator with a pent loop-gap design80 for the DEER experiments. For the pent loop-gap resonator the coupler position was set to provide the lowest possible Q-value, which was previously estimated to be about 50, corresponding to a bandwidth of about 700 MHz.80
Cu(II)–Cu(II) RIDME data were acquired at 20 K using the refocused five-pulse sequence shown in Fig. 2(a) with π/2-pulses set to 12 ns and π-pulses to 24 ns. The sequence was . The initial time t was set to d3,0 − d1, with d1 = 400 ns and d3,0 = 280 ns (negative initial t = −120 ns before the position of the primary echo), in order to measure both sides of the region around zero dipolar evolution time. The interpulse delay d2 was adjusted for each sample. If not stated otherwise, the mixing time was set to Tmix = 40 μs for the ruler 13, 15 and 3. To remove echo crossing and phase offsets, an eight-step phase cycle was used.63 ESEEM averaging81 was performed. To reduce orientation-selection effects, RIDME data were acquired at different observer fields29,33,53–55 over the full Cu(II) spectrum using a spacing of 150 G. Resulting time traces were normalized to unit intensity at the RIDME zero-time and summed weighted by the spectral intensity of each observer position. With this approach most observer orientations are expected to be excited uniformly. Orientation selection from the inverted spin is assumed to be negligible in RIDME experiments because relaxation events can occur in the full spectral range.
Cu(II)–Cu(II) DEER data were recorded at 20 K with the dead-time free 4-pulse DEER sequence .65 The delay τ1 was set to 400 ns and τ2 was adjusted for the different compounds. All observer pulses were 12 ns long and set 100 MHz off the center of the resonator mode f0. The pump pulse was a sech/tanh (HS) pulse of order 6.82,83 The length was 100 ns, the bandwidth 800 MHz and the frequency flank facing the observer pulse was offset from the observer frequency by 100 MHz.33 HS pulses were compensated for the resonator profile and non-linearity of the excitation chain.76 Two pulse setup versions were employed, one in which the observer pulses were set to higher frequency than the pump pulse, and one for which the observer pulses were lower in frequency. Similar to the RIDME experiments, measurements were performed with the observer pulses at the maximum of the field swept spectrum and at several offsets from the maximum in steps of 150 G, respectively. The DEER traces from the individual field positions were added corresponding to the spectral intensity of each observer position. Finally, the summed traces from both setups were normalized and added together. With such an averaging procedure, RIDME and DEER traces agree sufficiently well to indicate (see Supporting information part 1 (ESI1†)) that orientation selection effects are strongly suppressed.
Cu(II)–nitroxide RIDME data were acquired at 20 K using the resonance offset-compensated five-pulse RIDME sequence (2:2:2:2:1 length ratio of the frequency-swept pulses) following the previously reported procedure.56 To cope with orientation-selection effects the pulses had a bandwidth of 400 MHz, which was centered within the nitroxide spectrum. The π/2-pulses and the long π-pulse were 150 ns long, the short π-pulse was 75 ns long. Quarter-sine flanks were used with rise-times of 30 ns for the long pulses, and 15 ns for the short pulse. The interpulse delays were set to d1 = 6100 ns, d3,0 = 5980 ns, and d2 = 3500 ns for ruler 21, d1 = 4100 ns, d3,0 = 3980 ns, and d2 = 3500 ns for ruler 4, and d1 = 4100 ns, d3,0 = 3980 ns, and d2 = 5500 ns for ruler 23. These choices of d1 suppress the effect of a low-frequency artefact, which was reported for nitroxide detected RIDME.56 The longer d1 for ruler 21 was required to enable reproducible fitting of the RIDME background, which is less stable for data without pronounced oscillations. The mixing time was set to Tmix = 400 μs. To remove echo crossing and phase offsets, an eight-step phase cycle was used.632H ESEEM averaging was performed with the stationary echo protocol.56,81
The parameters of g and A tensors were calculated using different functionals and basis sets to test for variations throughout the DFT calculations. g-Tensor calculations were performed using (ma-)def2-TZVP and the hybrid density functional B3LYP in combination with the RIJCOSX approximation or using PBE094 without and with the RIJCOSX approximation. In addition, for the Cu(II)-ion a core properties basis set with three polarization functions CP(PPP)95 for the copper transition metal center was tested. The A tensor was calculated using the PBE0 functional in combination with the EPR II or III96 basis sets for light atoms as well as the B3LYP functional in combination with the EPR II basis set for light atoms. On the Cu(II)-ion either the (ma)-def2-TZVPP or the CP(PPP) basis set was applied. The results of the DFT computations are shown in ESI1.†
In the first step all experimental traces were analyzed with the DeerAnalysis97 open-source package (version 2018) whereby the standard functionality for processing dipolar evolution data was utilized (Tikhonov regularization in the time domain using the standard K(r,t) kernel which assumes only dipolar coupling as given in eqn (8). Where appropriate, g-values for computing the dipolar coupling constants ωdd (eqn (3)) were adjusted.
Second, the data were analyzed using the generic and exchange resilient neural networks as implemented in DeerAnalysis2018.66
Two further custom computational approaches were developed – the coupled exchange and distance distributions model and the regularization in the exchange-coupling domain procedure. These are described below.
The coupled exchange and distance distributions model was parameterized with 4 adjustable parameters: mean and standard deviation of distance and exchange-coupling distributions (see ESI1,† Fig. S8). Computations were performed in frequency domain using a point-dipole approximation applied to a weighted average position of the (distributed) electron spin at both ends of the molecule. Such a point-dipole approximation is permitted if the two spin density distributions have axial symmetry with coinciding symmetry axes. Spin coupling spectra obtained from experimental time-domain form-factors by Fourier transformation were used as simulation targets. Model spectra were obtained in the following way. A distance axis representing the spin–spin distance r was defined and used for generating model distributions of distances P(r) and exchange couplings J(r) for chosen combination of rmean, σr and Jmean, σJ values. Since P(r) and J(r) are coupled, each point of the distance distribution corresponds to a particular value from the isotropic exchange coupling distribution J(r) (ESI,† Fig. S8). For each point in P(r), an ideal Pake pattern was computed based on eqn (3); each single frequency from the pattern was combined with the corresponding (isotropic) exchange coupling value from J(r) thus yielding a kernel function for a particular value of rAB. Averaging over all points of the distance distribution P(r) resulted in a simulated spin–spin coupling spectrum for a given rmean, σr and Jmean, σJ combination.
Direct fitting to experimental frequency spectra in the 4-parameter space (i.e. rmean, σr and Jmean, σJ) was prone to become stuck in local minima of the error surface. Therefore, an iterative grid-search approach was adopted. In order to identify a global minimum, initial computation was always performed for a sufficiently broad range of model parameters using a rough grid for each parameter. Subsequently, parameter variation ranges were reduced while grids were made finer in order to find the best matching model for a given experimental data set. Two representative data sets in the presence and absence of exchange couplings ([Cu–PyMTA]–[Cu–PyMTA] 1n, n = 3, 5) were subject to analysis within the coupled exchange and distance distribution model. Details of the grid-search procedure for each case are given in ESI1.†
(8) |
(9) |
Prior to applying this computational approach to experimental data, it was subjected to an extensive validation using simulated data sets. When mathematical problems are solved by kernel inversion without regularization, the rows of the kernel must be pairwise orthogonal to ensure a stable solution. The pure dipolar evolution kernel K(t,r) computed with (8) does not fulfill this condition. The solution must thus be regularized. On the contrary, a hypothetical purely exchange coupling kernel K(t,J) could be constructed such that the pairwise orthogonality requirement is fulfilled thus allowing for computing correct exchange coupling distributions, potentially without regularization. However, determining the character of P(J) (i.e. ferromagnetic or antiferromagnetic) would be impossible in this way as the entries corresponding to ferro- and antiferromagnetic couplings of the same strength would be exactly the same according to eqn (9) for cosine being an even function. The two situations are distinguishable only in the presence of the dipolar contribution of comparable magnitude, whose sign is known. Thus the presence of dipolar coupling contributions in the K(t,J)p(r) kernel on one hand makes the rows of the kernel less (pairwise) orthogonal resulting in less stable solutions requiring regularization. On the other hand distinguishing between the ferro- and antiferromagnetic regimes becomes possible. Results of the calibrating simulations are summarized in ESI1.†
The Cu(II)–nitroxide rulers were synthesized via a reaction sequence, in which the nitroxide is attached to the spacer first and then the ligand moiety (Scheme 1).56 In this way, synthesis effort needed for the variation of the metal ion chelating moiety is kept to a minimum. Important details about this reaction sequence can be found in a recent publication which describes the synthesis of the short [Cu–PyMTA]–nitroxide and [Cu–TAHA]–nitroxide rulers 21 and 4.56 The long [Cu–PyMTA]–nitroxide ruler 23 was assembled largely in the same way, but the linking of the nitroxide moiety and the spacer was different. For the assembly of the short Cu–nitroxide rulers the nitroxide was attached via amide bond formation between a spacer with a terminal amino group and 2,2,5,5-tetramethyl-3-pyrrolin-1-oxyl-3-carboxylic acid (6a), the latter in form of an N-acylpyridinium chloride 6b (method A).56 In case of the assembly of the long Cu–nitroxide ruler 23, the nitroxide was first transformed into the N-(4-iodophenyl)amide 7, which was attached to the oligoPE 83via a Sonogashira–Hagihara coupling reaction (method B; Scheme 1). Method B was applied when the oligoPE backbone was found to be harmed under the conditions of method A. Later it was realized that this was an effect of long reaction times and could be avoided by limiting the reaction time to 1–2 h. Method B has its pros and cons. On the one hand it reduces the number of consecutive steps to be performed with the spacer and avoids to handle very small amounts of a highly water sensitive acylation reagent in cases where the amount of oligoPE is small. On the other hand, the Sonogashira–Hagihara coupling reaction is occasionally accompanied by reduction of a few percent of the nitroxide to hydroxylamine.56
Scheme 1 Synthesis of [Cu(II)–PyMTA]–nitroxide ruler 23. (a) (1) DMAP, CH2Cl2, (2) thionyl chloride, 0 °C to rt, 65 min; 6b was not isolated; (b) rt, 65 min, 98%; (c) Pd2(dba)3, CuI, PPh3, piperidine, THF, rt, 45 h, 77%; (d) K2CO3, MeOH, CH2Cl2, rt, 14.5 h, 91%; (e) (1) Pd2(dba)3, CuI, PPh3, iPr2NH, THF, rt, 45 h, (2) 2-methylbut-3-yn-2-ol, rt, 22 h, 44%; the addition of 2-methylbut-3-yn-2-ol serves to trap and thus polarly tag residual 4-iodo-PyMTA ethyl ester (11) for its easy chromatographic removal. (f) nBu4NF, THF, rt, 10 min; (g) [Cu(phen)(PPh3)2]NO3·0.5CH2Cl2, toluene, 67% in reference to compound 12; (h) (1) NaOH, EtOH, H2O, 18 h, rt, (2) pH adjusted to ∼7 through addition of HCl in H2O, rt, 80%. The protonation degree at pH 7 given with the formula of PyMTA–nitroxide 16 is based on the data for diethylenetriaminepentaacetic acid (DTPA) which is structurally similar to PyMTA.75 (i) CuCl2, NaOD, D2O, rt. For further details see ESI2.† Abbrevations: DMAP = 4-(dimethylamino)pyridine, phen = phenanthroline, THF = tetrahydrofuran, TIPS = triisopropylsilyl, TMS = trimethylsilyl, rt = room temperature. |
The [Cu–PyMTA]–[Cu–PyMTA] rulers 1n were prepared in the same way as the corresponding [Gd–PyMTA]–[Gd–PyMTA] rulers,99 only using CuCl2 instead of GdCl3 in the final step. When aiming at linear rulers with two identical spin labels, as is the case for the Cu(II)–Cu(II) rulers 1n and 3, an oxidative alkyne dimerization (Glaser coupling) of a building block of the type ligand–(EP)nE–H which has a terminal ethyne group giving a ruler of the type ligand–(EP)nB(PE)n–ligand with B standing for a butadiyne unit (B = EE), is most convenient in respect to the number of synthetic steps because the spacer length is doubled in one step. The butadiyne unit makes the spacer only slightly more bendable.52 This strategy has been exploited for the synthesis of rulers with two nitroxide units.102 However, the PyMTA ester moiety turned out to be incompatible with the simultaneous presence of Cu(II) and oxygen,103i.e. with the typical conditions for such a dimerization. Under these conditions, the aminomethyl substituent at the pyridine ring is transformed into a formyl group.103 Although being a slow side reaction, it matters because also the dimerization of long spacers proceeds rather slowly and it is difficult to remove the side product. Therefore, the short synthetic pathway was no option for the synthesis of long [Gd–PyMTA]–[Gd–PyMTA] rulers. In contrast, the TAHAester moiety proved inert under these conditions and, therefore, the oxidative dimerization was employed in the synthesis of the [Cu–TAHA]–[Cu–TAHA] ruler 3 (Scheme 2): the TAHA unit was attached to half of the final spacer via Sonogashira–Hagihara coupling of spacer building block 8199 with 4-iodo-TAHA ethyl ester (17) and the alkyne protecting trimethylsilyl group of coupling product 18 was removed. Glaser coupling of alkyne 19 completed the scaffold assembly. Despite the inertness of TAHA ester 17 against the simultaneous presence of Cu(II) salts and oxygen, we trapped all copper and palladium ions with metal scavengers during work-up of the reaction mixtures because column chromatography can be insufficient to remove the metal ions completely and even traces of residual metal ions can cause uncontrolled reactions with terminal alkyne groups. To avoid ester hydrolysis (there are six ester groups per molecule) under the basic conditions of desilylation, we applied water free potassium carbonate, thereby enabling strict exclusion of water.
Scheme 2 Synthesis of [Cu–TAHA]–[Cu–TAHA] ruler 3. (a) PdCl2(PPh3)2, CuI, iPr2NH, THF, rt, 21.5 h, 82%; (b) K2CO3, EtOH, THF, rt, 19 h, 95%; (c) PdCl2(PPh3)2, CuI, iPr2NH, THF, in air, rt, 23 h, 84%; (d) Bu4NF, THF, rt, 75 min; (e) [Cu(phen)(PPh3)2]NO3·0.5CH2Cl2, toluene, 30 °C, 69.5 h, 78% in reference to compound 20; (f) (1) NaOH, H2O, rt, (2) HCl in H2O, rt, 54%; the average protonation degree of TAHA–TAHA 23 is 1.7 protons per TAHA calculated with the pKa values reported for TAHA.101 (g) CuCl2, NaOD, D2O, rt. For further details see ESI2.† Abbreviations: phen = phenanthroline, THF = tetrahydrofuran, TIPS = triisopropylsilyl, TMS = trimethylsilyl, rt = room temperature. |
After scaffold assembly, as in the synthesis of the other rulers, the triisopropylsilyl groups protecting the ethyne units in the side chains were removed and the deprotected ethyne moieties were used to attach the PEG side chains by the Cu(I) catalyzed alkyne–azide cycloaddition of alkyne 21 with PEG–azide 14 giving PEGylated compound 22. The late attachment of the polar PEG side chains allowed us to use standard extraction and standard column chromatography on silica gel for purification of all products up to this step. The obtained PEGylated compound 22 was saponified to ruler precursor TAHA–TAHA 23 (denoted ‘3-pre’ in the UV/Vis experiments) containing the free ligands ready for complexation. By addition of an aqueous solution of copper(II) chloride the ruler [Cu–TAHA]–[Cu–TAHA] 3 was obtained.
When preparing the rulers, about 5% less than the calculated required amount of copper salt was added to the ruler precursor in order to avoid the presence of unchelated copper ions and therefore copper ions which are spectroscopically different from the chelated ones.
Analysis of the Cu(II) EPR spectra in combination with DFT calculations was performed in order to estimate the spin density distribution at the Cu(II) complexes, and estimate the fraction of spin density transferred to the spacer's conjugated π-system. Fig. 3 shows the X-band CW EPR spectra of Cu–PyMTA (a) and Cu–TAHA (b) and their associated rulers at 140 K. Lineshape fits are presented in ESI1.†
The EPR spectra of Cu–PyMTA and [Cu–PyMTA]–[Cu–PyMTA] rulers 13 and 15 overlay well. In particular, the positions of three out of the four low-field peaks due to the Cu(II) hyperfine splitting in the parallel orientation of the g-frame are well resolved, and these positions are the same for Cu–PyMTA and [Cu–PyMTA]–[Cu–PyMTA] rulers 13 and 15. The spectra of the [Cu–PyMTA]–nitroxide rulers 21 and 23 were obtained with a microwave power that heavily saturated the nitroxide radicals. Still the nitroxide line dominates the EPR lineshape around g ≈ 2, which makes it difficult to interpret the contribution of the Cu(II) ion to this region visualized by a grey background in Fig. 3(b). However, in the parallel orientation the same three peaks due to Cu(II) hyperfine splitting are clearly visible at the same positions as for the other Cu–PyMTA compounds. Somewhat stronger broadening is observed for these three lines of the [Cu–PyMTA]–nitroxide rulers. Importantly for our analysis, the strong hyperfine anisotropy in the Cu(II) EPR spectrum is due to the dipolar character of the hyperfine coupling, which appears due to the non-spherical distribution of the electron spin density around the Cu(II)-nucleus. A strong difference in the electron spin density distribution between Cu–PyMTA and [Cu–PyMTA]–[Cu–PyMTA] ruler would result in a clear difference of this hyperfine splitting, which is not observed for the [Cu–PyMTA]–[Cu–PyMTA] rulers 13 and 15. Thus, also in case of these rulers the major part of the electron spin density must be limited to the atoms of the Cu–PyMTA moiety. On the other hand, the lineshape of [Cu–PyMTA]–[Cu–PyMTA] ruler 11, which has the shortest Cu–Cu distance, strongly deviates from the one of the [Cu–PyMTA]–[Cu–PyMTA] rulers 13 and 15 featuring longer spacers. This indicates that the two Cu(II) centres in this compound are significantly exchange-coupled. The broadening due to pure dipole–dipole interaction at the anticipated Cu(II)–Cu(II) mean distance of about 2.1 nm would not be sufficient to change the CW EPR lineshape of ruler 11 that strong. Spin-density transfer into the ruler π-system is an explanation of the EPR lineshape change for the ruler 11. Also, a large electron–electron coupling is consistent with a strongly reduced modulation depth in PDS experiments (see ESI1†), as large couplings cannot be excited by the applied microwave pulses.104 Calculation of a broad unknown distance distribution along with simultaneous determination of a distribution of exchange couplings would be highly ambiguous, since any distribution of dipolar couplings can be substituted by an appropriate distribution of exchange couplings. The DFT-based estimates of the spin density distributions in such samples would not provide sufficient accuracy of the dipole–dipole coupling estimate to stabilize the fitting procedure. Even more importantly, a set of thermally excitable geometries would need to be considered in order to obtain an estimate of the distribution of dipole–dipole couplings. We thus refrained from detailed analysis of [Cu–PyMTA]–[Cu–PyMTA] ruler 11 by PDS experiments.
The Cu–PyMTA spectrum was analysed previously33 and is best reproduced using a rhombic g tensor (gx = 2.048, gy = 2.078, gz = 2.253) as well as a rhombic A tensor (Ax ∼ 25 MHz, Ay ∼ 75 MHz, Az ∼ 510 MHz). Note however that the broad lines render extraction of the perpendicular tensor components ambiguous (see ESI† for details).
Concerning the TAHA-containing compounds, the EPR spectrum for the [Cu(II)–TAHA]–nitroxide ruler 4 does almost perfectly overlay with the corresponding EPR spectrum of Cu–TAHA if the contribution of the nitroxide is neglected. The lineshape simulation for Cu–TAHA gives axial g and A tensors (gx,y = 2.066, gz = 2.284, Ax,y = 23 MHz, Az = 498 MHz). In the case of [Cu–TAHA]–[Cu–TAHA] ruler 3 some broadening and a slight shift of the parallel Cu-hyperfine components occur relative to the other two TAHA-based compounds. Thus, for simulating the spectrum of [Cu–TAHA]–[Cu–TAHA] ruler 3 it is not sufficient to only adjust strains, but a small shift in gz and Az is needed (see ESI1†). Nevertheless, since this shift is small we can assume that the spin density transfer into the effective spacer must be weak in the series of [Cu–TAHA]-based compounds.
The clearly visible components of the Cu(II) hyperfine quartet in the parallel region with peak positions determined by gz and Az serve as a fingerprint for the local Cu(II) coordination according to the analysis by Peisach and Blumberg.105 Analyzing gz and Az of Cu–PyMTA and Cu–TAHA accordingly shows (ESI1,† Fig. S24) that both complexes fall into an ambiguous region with three or two directly coordinating nitrogens, with a tendency more towards three coordinating nitrogen atoms for Cu–PyMTA and more towards two coordinating nitrogen atoms for Cu–TAHA.
The degree of spin delocalization was also studied by DFT calculations of Cu–PyMTA and Cu–TAHA. Cu–PyMTA is assumed to be fully deprotonated at pH 7 because a complete deprotonation is reported for the structurally similar complex Cu(II)–DTPA with DTPA being diethylenetriaminepentaacetic acid.75 For Cu–TAHA on average 1.7 protons are anticipated.56 Nevertheless, in solution an equilibrium of different protonation states is expected and thus geometries for the complexes with different numbers of protons were considered. The optimized geometries were sensitive to the protonation degree, which also led to differences in the tensor symmetries and calculated values of g, AN as well as ACu (see ESI1†). The g and AN parameters are relatively stable for different functionals and basis sets, but ACu varies strongly and in some cases is qualitatively wrong for the TZVPP basis, indicating that the degrees of freedom in this basis, particularly in the core region, are insufficiently accurate to model core properties such as metal hyperfine couplings. For some of the protonation states, the calculated g and ACu tensor symmetries match the experimental data qualitatively. However, we did not find sufficient quantitative agreement to differentiate between different structures. For Cu–PyMTA and [Cu–TAHA]–[Cu–TAHA] ruler 3, adding the largest DFT calculated 14N hyperfine values (∼50 MHz) to the lineshape simulation is in agreement with the experimental data, i.e. this 14N hyperfine coupling can be accommodated as an unresolved coupling in the lineshape fits. Similarly, for Cu–TAHA an experimental upper limit of about 40 MHz is found for the z-component of the nitrogen hyperfine coupling or an Aiso of 35 MHz (in DFT the largest 14N hyperfine coupling is ∼44 MHz, see ESI1†). This suggests that no spin density on nitrogen atoms beyond 13% for Cu–PyMTA and [Cu–TAHA]–[Cu–TAHA] ruler 3 and beyond 10% for Cu–TAHA is expected. This is in agreement with the Mulliken population analysis that shows throughout all studied cases a spin density of 60–80% on the Cu(II)-ion and no significant spin density on atoms beyond the first coordination sphere. Spin density distributions are shown in ESI1,† Fig. S25 and S26. This indicates that the spin density remains approximately within a sphere with a radius of about 2–3 Å around the Cu(II) ion and strongly peaks at the ion. Compared to a spin–spin distance of 2.4 nm or larger such an extension in space is rather small. Furthermore, this distribution is within the most narrow distance distribution widths computed from the PDS data of these rulers (see the following sections). Thus, taken together with the results of the CW EPR lineshape analysis, a point-dipole approximation with localization of the spin density at the copper ions appears to be a good approximation and will be assumed in the following analysis of PDS data for this series of compounds.
Note that Cu(II) coordination can be fairly variable in terms of changes in bond lengths and bond angles. This can lead to a distribution in g-tensor value and, more importantly, of orientations.106 Here we strongly reduce the effect of orientations by the RIDME-based effective broad-band measurements and by the detection position averaging. Additionally, the well resolved spectral features in the Cu(II) EPR spectra also indicate rather narrow g-values distributions, which in turn suggest rather low complex plasticity. Furthermore, we did not observe any measurable contribution of the free Cu(II) EPR signal in the CW EPR spectra of studied model compounds. This suggests rather low Kd for the Cu(II) chelates studied: Kd here cannot be higher than 1–2 μM, and more realistically it lies yet lower, somewhere in the nM range.
The UV/Vis spectra of the rulers and the ruler precursors are shown in Fig. 4. As predicted, the long wavelength cutoff in the extinction spectra (which corresponds to the lowest energy excitation) correlates with the number of spacer units n for each ruler type: the longer the ruler, the longer the cutoff wavelength. This is in agreement with UV/Vis studies of conjugated polymers.107–109 A list of the lowest energy peak wavelengths can be found in ESI1.† Note that the near UV spectral range (λ = 300 nm and below) is shown, but not analyzed in the following, and the spectra in this range do not vary significantly between different rulers.
Fig. 4 UV/Vis spectra of the ruler precursors (‘-pre’) and corresponding rulers, normalized to the maximum of the lowest energy peak of each individual trace (indicated by cross). Ruler precursors and rulers (a) with TAHA ligands; (b) with one PyMTA ligand and one nitroxide moiety; (c) with two PyMTA ligands; all maxima and shifts are tabulated in ESI1.† |
In the case of the TAHA-based compounds (ruler precursors 3-pre and 4-pre and the corresponding rulers 3 and 4) loading of the ruler precursors with Cu(II) does not affect the long-wavelength cutoff area of the UV/Vis spectra significantly. In the case of PyMTA-based compounds (ruler precursors 1n-pre and 2n-pre with n = 1, 3, 5, and the corresponding rulers 1n and 2n) we observed weak red-shifts of the long wavelength cutoff upon loading the ruler precursors with Cu(II). This red-shift is larger for the ruler precursor/ruler pairs with shorter spacer lengths. The largest absolute shifts are observed for the pairs 11-pre/11 (Fig. 4(c), red-shift of 6.8 nm) and 21-pre/21 (Fig. 4(c), red-shift of 4.3 nm). For the ruler precursor/ruler pairs with intermediate spacer lengths the red-shift is smaller, but still clearly visible (ruler precursors 13-pre and 23-pre and the corresponding rulers 13 and 23, Fig. 4(b and c)). The red-shift of the long wavelength cutoff for the pair ruler precursor 13-pre/ruler 13 with two PyMTA moieties, one at each end of the spacer, is slightly larger than the corresponding red-shift for the case of the pair ruler precursor 23-pre/ruler 23 with only one PyMTA moiety and one nitroxide moiety at the other end of the spacer. Note that loading of the ruler precursor corresponds to the addition of two Cu(II) metal ions in the first case, and of only one metal ion in the latter case. For the pair with the longest spacer, ruler precursor 15-pre/ruler 15, the red-shift is nearly invisible.
The maximum strength of the spin–spin couplings for the studied rulers is below 20 MHz (see the following sections). This indicates that the electron spin density of the Cu(II) ion is only weakly coupled to the π-system of the effective spacer in the electronic ground state. However, first, in the UV/Vis spectra we observe excitation of the entire conjugated π-system, not only the unpaired electron, which is relevant for EPR. Second, the UV/Vis spectra depend not only on the electronic configuration of the ground state, but also on the electronic configuration of the excited state, where stronger delocalization of the unpaired electron of the Cu(II) ion cannot be excluded a priori. Thus, the overall picture given by the UV/Vis data is in line with the assumption of the weak exchange between paramagnetic centres in these compounds.
Fig. 5 PDS data for the studied ruler compounds (see ESI1† for the primary data). (a and c) Background-corrected RIDME form factors F(t) and (b and d) apparent distance distributions P(r) computed neglecting exchange coupling. (a and b) Molecular rulers with one paramagnetic center being Cu(II) and one being nitroxide: [Cu–PyMTA]–nitroxide rulers 2n (n = 1, 3) and [Cu–TAHA]–nitroxide ruler 4. (c and d) Both paramagnetic centers are Cu(II): [Cu–PyMTA]–[Cu–PyMTA] rulers 1n (n = 3, 5), [Cu–TAHA]–[Cu–TAHA] ruler 3. |
Given the identical architecture of the PyMTA-based compounds with different spacer lengths, it is excluded that the shorter spacers exhibit higher flexibility, and the corresponding compounds possess truly broader distributions of spin–spin distances to which such a faster dampening could be related. Indeed, the measured distance distribution widths should be comparable to those found for a series of stiff rod-like nitroxide biradicals exhibiting pure dipole–dipole character of spin–spin couplings.52 Neither is a higher flexibility of these spacers supported by earlier distance measurements on the same [M–PyMTA]–[M–PyMTA] compounds using M = Gd(III)71,110 or Mn(II)81 as metal centers. These measurements resulted in metal-to-metal distance distributions of similar widths centered around 3.4 nm and 4.7 nm for the short and long spacer, respectively. In contrast, it was shown earlier37 that such a dampening can result from distributed J-couplings. In conjugated π-systems – as present in the [Cu–PyMTA]–nitroxide and the [Cu–PyMTA]–[Cu–PyMTA] rulers – non-negligible contributions from exchange couplings have been observed for spin–spin distances up to 3.6 nm.45 The exponential decay of the exchange coupling with distance or a number of structurally identical linker units would be able to explain its absence for the compounds with longer spacers.44,51
Very similar qualitative conclusions can be made for the present series of compounds simply based on the apparent distance distributions presented in Fig. 5, and on the nearly identical CW EPR spectra within the series of PyMTA-based and TAHA-based compounds. First, the Cu(II) CW EPR spectra of the PyMTA-based rulers 21, 23, 13 and 15 are essentially the same, which strongly suggests that all Cu–PyMTA containing rulers show the same distribution of conformations of the Cu–PyMTA moiety, and the same holds true for all Cu–TAHA containing rulers. Since the PyMTA-based rulers exhibit narrow apparent distance distributions when the spacer is long and broad ones when it is short, it is reasonable to assume that the apparent short spin–spin distances in the latter case result from distortion of the dipolar spectra by exchange interaction, whereas the true distribution of spin–spin distances is narrow. Second, the apparent distance distributions show either distances that are close to the expected spin–spin distance or shorter. However, the spin–spin distances which are predicted based on the compound structures and which are experimentally observed in the case of the longer spacers, correspond to the nearly closest placement of the Cu(II) ion with respect to the partner Cu(II) or nitroxide spin. Therefore, yet significantly shorter spin–spin distances are not feasible with the given structures of the ruler compounds. Since apparent short distances in the distance distributions in Fig. 5 correspond to stronger spin–spin couplings than expected from dipolar interaction, we conclude that exchange interaction is the reason for the broad apparent distance distributions.
If we interpret the widths of the apparent distance distributions in Fig. 5 as a measure of the strengths of exchange couplings, some further useful comparisons can be made. In the case of [Cu–TAHA]–nitroxide ruler 4 and [Cu–TAHA]–[Cu–TAHA] ruler 3, the three single bonds between the N atoms being in contact with the Cu(II) ion and the effective spacer are expected to decouple the conjugated π-system of the spacer from the unpaired electron of the Cu(II) ion. This is supported by the UV/Vis data described above that do not show a red-shift in the absorption spectra for the rulers 3/4 involving at least one TAHA ligand as compared to corresponding ruler precursors. Thus, the range of spacer lengths for which we expect exchange coupling is presumed to be much shorter. However, while UV/Vis spectra can reveal conjugation of the Cu(II) ion with the π-system of the effective spacer, as found for the rulers with Cu(II)–PyMTA moieties, they would not be able to detect weak exchange coupling through space or through saturated bonds, which likely takes place for the rulers with Cu–TAHA moieties. The assumption of weak exchange coupling when Cu–TAHA is present is supported by the observation of a narrow single peak in the apparent distance distribution of [Cu–TAHA]–nitroxide ruler 4. In contrast, [Cu–TAHA]–[Cu–TAHA] ruler 3 shows several additional peaks contributing to a broad background next to a relatively narrow peak at 3.3 nm, the latter being in agreement with the distance anticipated based on the molecular structure. The comparison of the PDS data for the two Cu(II)–nitroxide rulers 21 and 4 allows us to conclude that coupling of the Cu(II)–TAHA spin with the spacer π-system is weaker than for Cu(II)–PyMTA. The same conclusion can be made by comparing the PDS data of [Cu–PyMTA]–[Cu–PyMTA] ruler 13 and [Cu–TAHA]–[Cu–TAHA] ruler 3. The broad distance peak obtained for [Cu–TAHA]–[Cu–TAHA] ruler 3 overlaps with the more narrow peak at the anticipated spin–spin distance, while for [Cu–PyMTA]–[Cu–PyMTA] ruler 13 the anticipated distance peak is not well visible, rather a broad exchange-distorted distribution is observed. This indicates that the mean exchange coupling and the width of the exchange coupling distribution is smaller for 3 as compared to 13 in agreement with the interruption of the conjugation between the π-system of the effective spacer and the Cu(II) ion in the Cu–TAHA moiety.
Furthermore, based on stronger distortions of the apparent distance distribution, we can state that the exchange interaction between two Cu–TAHA complexes (ruler 3) is stronger than between one Cu–TAHA complex and one nitroxide moiety (ruler 4). Note that the expected spin–spin distance in [Cu–TAHA]–[Cu–TAHA] ruler 3 is even longer than in the [Cu–TAHA]–nitroxide ruler 4, which makes the last conclusion even stronger. The comparison of [Cu–PyMTA]–[Cu–PyMTA] with [Cu–TAHA]–[Cu–TAHA] is also in line with stronger coupling to the spacer conjugated π-system for Cu–PyMTA moieties as compared to the Cu–TAHA moieties. To summarize these observations, the Cu–TAHA complex exhibits a weaker exchange coupling with the nitroxide moiety than the Cu–PyMTA complex and the exchange coupling is weaker in Cu(II)–nitroxide than in Cu(II)–Cu(II) pairs. Thus, the propensity to induce exchange interaction with the partner moiety in the studied series of three paramagnetic units has the order nitroxide < Cu–TAHA < Cu–PyMTA.
Possible distortions of the PDS data can be excluded as follows. Note that apart from the contribution of spectral diffusion to the RIDME experiment, which results in a faster, strongly curved background shape, orientation-averaged RIDME data are semi-quantitatively concurrent with orientation-averaged DEER data of all measurements in this work (see ESI1,† Fig. S1 and S2). This strongly indicates that the additional peaks in some of the cases result neither from incomplete ESEEM averaging nor from orientation selection effects. These conclusions are further supported by RIDME background reference measurements for Cu–TAHA as well as Cu–PyMTA shown in ESI1,† Fig. S7. Aggregation of the rulers appears unlikely considering the solvent polarity. Moreover, the nitroxide–nitroxide DEER measurements with the Cu(II)–nitroxide rulers 21 and 4 did not reveal any dipolar modulation above the noise level as presented in ESI1.† Note further that for the Cu(II)–Cu(II) and Cu(II)–nitroxide rulers 15, 23 and 4 we see pure dipolar couplings corresponding to the anticipated narrow spin–spin distance distributions and hence no indications for aggregation.
We cannot fully exclude some weak residual contributions from orientation selection. However, such a contribution cannot explain the observed broad apparent distance distributions (see additional analysis in ESI1†). Finally, we considered whether the partially charged end groups in the TAHA moiety might influence spacer conformation and thus cause additional peaks in the distance distribution. With a Bjerrum length of 7.1 Å in water and a dielectric constant of about 60 of the water/glycerol mixtures111 that we use, one expects that the Coulomb interaction between two unit charges drops below thermal energy at about 9 Å. Since the distances between copper ions are much larger in all ruler compounds, we disregard such electrostatic effects.
In the following we report on the performance of our data analysis routines to extract exchange couplings, assuming that all mentioned distortions are weak and do not change the overall picture.
The apparent distance distribution of the shorter compound 13 is broad and bimodal. The longer-distance peak, centered at about 3 nm, roughly agrees with the anticipated distance for this ruler, but it is too broad, while the second peak is also broad and shifted towards shorter distances. This is taken as indication of the presence of distributed exchange couplings which may, in combination with dipolar couplings, result in the broadly distributed total spin–spin coupling frequencies. Such frequencies may or may not coincide with dipolar couplings expected for the ruler 13 in the absence of exchange coupling. One may further speculate that the presence of the distance band at around 3 nm in this case may correspond to a fraction of rulers with weaker exchange couplings distributed around zero. Based on this preliminary inspection of corresponding experimental data, [Cu–PyMTA]–[Cu–PyMTA] compounds 15 and 13 could serve as test cases for situations of the absence of exchange coupling (for 15) and the presence of distributed exchange couplings (for 13).
Fig. 6(a) shows spin coupling spectra corresponding to the best solutions of the grid search procedure for both distribution types used. Except for the high-frequency regions of the spectra marked by arrows, we observe an excellent match between both grid search results and experimental data. Furthermore, spectra obtained from the grid search analysis nearly perfectly coincide with the spectrum corresponding to the solution obtained using regularization in the distance domain and a standard dipolar K(t,r) kernel.
The most likely reason accounting for most of the aforementioned deviations in the high-frequency regions between the modelled and experimental spectra of 15 is residual orientation selection, which might be present in the measured orientation-averaged data. Indeed, for a single, relatively narrowly distributed distance one expects to observe a typical Pake-pattern like frequency spectrum with the perpendicular- and the parallel-to-field singularities well seen and correctly weighted. In the case of compound 15, the position and the shape of the perpendicular-to-field features (or “horns”) of the spectrum match well. The parallel-to-field features (or “wings”) are underrepresented in the experimental data compared to all simulations as well as to the Tikhonov regularization fit. This indicates that these parallel orientations of the interspin axis with respect to the magnetic field were under-sampled during data acquisition.
The spin–spin coupling situation in the shorter compound 13 is more complex. The spectrum (Fig. 6(b)) is comprised of well-pronounced moderately broadened horns at about ±1.75 MHz and a very broad contribution extending beyond ±12 MHz. Such a frequency content together with the very stiff nature of the spacer indicates a distribution of exchange couplings.
For Lorentzian distributions, the initial search using a large rough parameter grid revealed a global minimum. Subsequent searches with finer parameter grids around that point yield the best-fit solution with rmean = 2.53 nm, σr = 0.09 nm, Jmean = −5.38 MHz and σJ = 0.04 MHz, (Fig. 6(b)). Contrary to this, the initial search using a large rough parameter grid and Gaussian distributions revealed two solutions with similar quality (rmsd figures). Subsequent fine grid searches in the vicinity of these solutions yielded the following two results featuring somewhat higher and somewhat lower exchange coupling values: rmean = 2.79 nm, σr = 0.36 nm, Jmean = −8.19 MHz and σJ = 0.1 MHz and rmean = 2.5 nm, σr = 0.16 nm, Jmean = −5.5 MHz, σJ = 0.08 MHz. This latter solution with shorter mean distance and smaller exchange coupling constant is very close to the result obtained using Lorentzian distributions. Details of the grid search analysis of compound 13 are given in ESI1.† Both solutions corresponding to the Gausssian distributions are very similar in quality to each other. However, in terms of describing the experimental data both solutions using the Gaussian distributions are worse than the solution obtained with the Lorentzian distributions that seemingly describes the experimental data in a nearly perfect way. Note that all solutions either with the Lorentzian or Gaussian distribution shapes underestimate the mean interspin distance expected from the structure of 13. This fact, together with the discussed deviations between the model solutions in terms of the interspin distance and the exchange coupling values may be attributed to the use of simplified distribution functions and the model assumption of a strict correlation between the dipole–dipole and exchange coupling distributions – which may have led to error compensation during the computations. The better fit of the experimental data with Lorentzian shaped distributions as well as the finding of a single best solution in that case indicates most probably that the underlying distributions are better described by slowly decaying tails implied by the Lorentzian shape. Nevertheless, all grid search calculations agree qualitatively in the description of the coupling situation in [Cu–PyMTA]–[Cu–PyMTA] compound 13: a moderate ferromagnetic coupling is predicted. This prediction was subject to further analysis using the approach of regularization in the exchange-coupling domain as is discussed in the next section. Note finally that in the case of a narrow distance distribution (compound 15) the absence of significant exchange coupling was recognized by this fitting approach, despite the potential possibility to describe part of the couplings with distributed exchange.
We thus conclude that both investigated models are able to recognize the presence or absence of distributed exchange couplings and to estimate the strength of distributed exchange couplings as long as other distortions, such as orientation selection or multi-spin contributions, can be ruled out or suppressed. However, for a precise, quantitative extraction of exchange coupling parameters, precise knowledge of the spin–spin distance distribution would be required.
When the spin–spin distance distribution can be estimated by some independent method, the use of regularization in the exchange domain appears advantageous, since it is not restricting the shape of the exchange coupling distribution by an ad hoc model. Due to the uncorrelated dipolar and exchange coupling distributions, this approach is also more stable with respect to small distortions of the spin–spin coupling spectrum, which may result from orientation averaging imperfections or any other experimental artefacts.
Fig. 8 shows that the generic neural network is able to capture correctly the strongly curved RIDME background in the primary data of both rulers 13 and 15. Furthermore, the associated distance distributions match those obtained using the standard K(t,r) kernel and Tikhonov regularization well: in case of ruler 15 that exhibits no exchange coupling this match is perfect, while in the case of ruler 13 in which (distributed) exchange coupling is present apparent distance distributions obtained by both methods agree reasonably well. On the contrary, the exchange resilient network failed when applied to these RIDME data. In case of the shorter compound 13, the background decay was apparently well described by the neural network analysis. However, the resulting distance distribution – a single relatively broad peak at around 2.5 nm – strongly contradicts the interspin distance expected for this compound based on its chemical structure. In case of the longer compound 15, the background decay calculation failed altogether resulting in an unrealistic distance distribution.
Similar results were obtained when analyzing the other compounds as shown in ESI1.† We thus conclude that the exchange resilient network is not capable of realistically approximating the distributed exchange couplings in the presented PDS data, probably because the networks were trained with single-valued exchange couplings. This problem may not be easy to fix by extending the training set, since the problem of simultaneous uncorrelated distribution of dipole–dipole and exchange couplings does not have a unique solution unless further restraints are introduced. Analysis of all the other compounds was therefore performed by regularization in the exchange-coupling domain with the K(t,J)p(r) kernel.
A spin pair consisting of a Cu(II) ion coupled with a nitroxide radical represented here by the [Cu–PyMTA]–nitroxide rulers 2n (n = 1, 3, Fig. 9(a and b)) behaves very similar to the spin pair consisting of two coupled Cu(II) ions represented by the [Cu–PyMTA]–[Cu–PyMTA] rulers 1n (n = 3, 5) as discussed above. Primary PDS data of the long [Cu–PyMTA]–nitroxide ruler 23 exhibit long lasting oscillations extending beyond the experimentally chosen detection window of 5 μs thus indicating that the underlying interspin distance distribution must be very narrow (Fig. 9(e)). Direct processing of these data using a standard dipolar K(t,r) kernel resulted in a narrow peak centered at 3.77 nm matching the prediction of 3.80 nm based on the molecular structure. The calculation of an exchange coupling distribution in this case reveals essentially its absence, as expected (Fig. 9(c)). Primary PDS data for the shorter [Cu–PyMTA]–nitroxide compound 21 show a very strongly damped oscillation, (Fig. 9(a)). Direct processing of this trace with the K(t,d) kernel expectedly results in a broad apparent distance distribution (Fig. 9(d)), being incompatible with the molecular structure. The calculation of the exchange coupling distribution using a narrow distance peak at 2.37 nm (estimated on the basis of the molecular structure) yields a broad band in the exchange coupling domain with a mean coupling value of about −2 MHz, (Fig. 9(c)).
For the [Cu–TAHA]-based rulers, the conjugated π-system of the spacer is not strongly coupled with the Cu(II) ion. It appears, however, that some exchange coupling mechanism exists for the Cu–TAHA complex studied here. The absence of exchange interaction is nicely observed in [Cu–TAHA]–nitroxide ruler 4 (Fig. 9i–l). For this compound, structure calculation suggests a mean interspin distance of about 2.7 nm which is in good agreement with the experimentally obtained narrow distance peak centered at 2.64 nm and the extraction of exchange coupling results in the narrow coupling band centered at zero. In contrast to this, for [Cu–TAHA]–[Cu–TAHA] compound 3 our analysis indicates a fraction of strained molecular structures with negligible exchange coupling and another fraction of such structures with distributed exchange coupling of ferromagnetic type. In comparison to [Cu–PyMTA]–[Cu–PyMTA] compound 13, the fraction with significant exchange couplings is smaller and the mean exchange coupling for this fraction is slightly weaker: ∼−2.5 MHz for [Cu–TAHA]–[Cu–TAHA] compound 3) instead of ∼−4 MHz for [Cu–PyMTA]–[Cu–PyMTA] compound 13. The presence of such a non-zero mean exchange coupling in [Cu–TAHA]–[Cu–TAHA] compound 3 indicates that there is a geometric arrangement that allows for coupling. As a speculation, this might be related to the spatial arrangement of those carboxylic groups of TAHA, which are not directly coordinated to the Cu(II) ion and facilitate through-space exchange coupling mechanism.
Finally, we should mention as an outlook that a comparison of the presented here PDS-based analysis and a more common CW EPR-based analysis of the exchange couplings in solution,113 based on the Kivelson's rotational averaging description,114 or in a frozen glassy state would be of significant interest. Similar studies in the frozen glassy state have been successfully performed to analyse systems with only through-space dipolar couplings present.115 However, when performing such a comparison one would have to face the problem that applicability range for simultaneous use of both techniques would be rather narrow: the lower limit for the detection of spin–spin couplings in CW EPR is usually around 10 MHz for nitroxide biradicals, which are characterized by narrower EPR lines, as compared to Cu(II). For the Cu(II) case this lower limit might get even larger, while the upper limit for detecting the spin–spin couplings with PDS is hardware dependent, and it usually lies somewhere between 50 MHz and 100 MHz.
In addition to this, we would like to mention that in contrast to the EPR data in the frozen state, the solution EPR reports on the rotationally averaged exchange coupling. In other words, this experiment would provide a single exchange coupling value, which might make a wrong impression that also in the frozen state all biradicals are characterized by the same coupling. The low-temperature PDS data we reported here clearly show that this is an oversimplification of the actual situation. Unravelling how does the exchange averaging happen in solution based on the ambient temperature solution CW EPR and low temperature PDS data would be an interesting study as well, however this would imply yet another level of complexity for the experimental work and for the analysis.
Analysis performed with DEERNet neural networks, as available in DeerAnalysis2018, was able to recover the strongly curved background in RIDME experiments using the generic model. However, within the exchange resilient model the neural network analysis was not able to identify the exchange couplings and provide the correct distance distribution for our data, likely because the present form of approximating exchange by a single J value in DEERNet appears to be unable to treat such cases with broadly distributed exchange couplings. Our work shows that approaches to determine exchange distributions in stiff synthetic molecules most likely require that the distance distribution and thus the distribution of dipolar couplings is known. Provided that the molecular geometry is known and the spin density distributions at both spin centres can be estimated, it should be possible to compute the dipole–dipole contribution independently and thus to determine the distribution of exchange couplings, e.g. by the here presented one-dimensional Tikhonov regularization.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d0cp03105d |
This journal is © the Owner Societies 2020 |