Distance measurement between trityl radicals by pulse dressed electron paramagnetic resonance with phase modulation
- 1Department of Chemistry and Applied Biosciences, Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
- 2Faculty of Chemistry and Center for Molecular Materials (CM), Bielefeld University, Universitätsstrasse 25, 33615 Bielefeld, Germany
Correspondence: Nino Wili (firstname.lastname@example.org)
Distance measurement in the nanometre range is among the most important applications of pulse electron paramagnetic resonance today, especially in biological applications. The longest distance that can be measured by all presently used pulse sequences is determined by the phase memory time Tm of the observed spins. Here we show that one can measure the dipolar coupling during strong microwave irradiation by using an appropriate frequency- or phase-modulation scheme, i.e. by applying pulse sequences in the nutating frame. This decouples the electron spins from the surrounding nuclear spins and thus leads to significantly longer relaxation times of the microwave-dressed spins (i.e. the rotating frame relaxation times T1ρ and T2ρ) compared to Tm. The electron–electron dipolar coupling is not decoupled as long as both spins are excited, which can be implemented for trityl radicals at Q-band frequencies (35 GHz, 1.2 T). We show results for two bis-trityl rulers with inter-electron distances of about 4.1 and 5.3 nm and discuss technical challenges and possible next steps.
Pulsed dipolar electron paramagnetic resonance (EPR) spectroscopy emerged as a powerful tool to measure distance distributions between electron spins in the nanometre range (Jeschke, 2012). This information is particularly useful when studying molecules and molecule assemblies that are intrinsically disordered or partially disordered or otherwise hard to crystallize and difficult to study with NMR or cryo-EM alone, e.g. certain membrane proteins (Bordignon and Bleicken, 2018) or protein-RNA complexes (Duss et al., 2014). The distance information is encoded in the magnetic dipole–dipole coupling between the electron spins, which depends on the inverse cubed distance, r−3. A plethora of different techniques have been introduced, most notably double electron–electron resonance (DEER) (Milov et al., 1984; Pannier et al., 2000), double quantum coherence (DQC) (Borbat and Freed, 1999), the single frequency for refocusing (SIFTER) (Jeschke et al., 2000), and relaxation-induced dipolar modulation enhancement (RIDME) (Kulik et al., 2001; Milikisyants et al., 2009). The limiting factor for all these pulse sequences is the electron phase memory time Tm, which determines the maximal dipolar evolution time and thus the longest distance that can be measured. In many cases, the phase memory time can be prolonged by deuterating the solvent or even the whole protein (Georgieva et al., 2010; Ward et al., 2010; Schmidt et al., 2016). However, such an approach is costly and is rarely feasible; e.g. it is very difficult for membrane proteins in a lipid bilayer and impossible for in-cell work.
In recent years, several groups tried to use dynamical decoupling sequences based on multiple refocusing pulses (also known as Carr–Purcell sequences) in order to prolong the coherence times (Borbat et al., 2013; Spindler et al., 2015). Although shaped pulses significantly improved the fidelity of EPR experiments, pulse frequency band overlap and non-uniform inversion are still a problem in these sequences and can lead to artifacts, which may be corrected if traces with sufficient signal-to-noise ratio and only moderately decaying background can be acquired (Breitgoff et al., 2017). Nevertheless, the improvements in Tm so far are of the order of a factor of 2, which only marginally (though sometimes decisively) improves the longest attainable distance.
Here we propose a sequence where the longest dipolar evolution time is, in principle, limited by the rotating frame relaxation time T2ρ, which is often much longer than Tm (for a discussion of T2ρ vs. the more familiar T1ρ, vide infra). The complete dipolar evolution takes place during strong microwave irradiation. This decouples the electron spins from the surrounding nuclei (Jeschke and Schweiger, 1997) while the electron–electron coupling is still active. The spin manipulation during the strong microwave irradiation is achieved by short intervals of sinusoidal phase modulation. The frequency of this modulation needs to match the Rabi or nutation frequency of the spin-locking irradiation.
The latter approach was discovered more than once in the history of magnetic resonance. It traces back to investigations of Redfield on “rotary saturation” (Redfield, 1955). Hoult introduced the related idea of longitudinal field modulation for nutation frequency selective pulses to magnetic resonance imaging (MRI) (Hoult, 1979). Grzesiek and Bax picked up Hoult's idea but used a phase-modulation scheme instead and applied it to homonuclear mixing in solution-state NMR (Grzesiek and Bax, 1995). They termed the technique “Audio-frequency NMR in a nutating frame”, because their phase-modulation (PM) frequency is in the audible range, and the pulse sequences effectively take place in a frame that nutates with the Rabi frequency of the spin lock. Independently, Jeschke used longitudinal field modulation during a spin lock for pulse EPR (Jeschke, 1999) and used the term “dressed EPR”, because the spins are dressed by the microwave field during the spin lock. This term is borrowed from quantum optics (Cohen-Tannoudji et al., 1992). The idea of dressed EPR originated in artifacts in hyperfine-decoupled electron–nuclear double resonance (ENDOR) spectra, which appear if the radio-frequency coil is not aligned perfectly perpendicular to the static field (Jeschke and Schweiger, 1997). Much later, it was also realized that field modulation should also prolong Rabi oscillations in the presence of inhomogeneous microwave fields (Saiko et al., 2018). Recently, the quantum information processing community picked up the idea of dressing electron spins in order to prolong coherence times (Laucht et al., 2016, 2017; Cohen et al., 2017). During the writing of this paper, Chen and Tycko came up with the idea of phase modulation during a spin lock independently again and used it for slice selection during off-resonance spin locks in solid-state, DNP-enhanced MRI (Chen and Tycko, 2020).
Here we combine the ideas of applying pulse sequences to dressed spins (Grzesiek and Bax, 1995; Jeschke, 1999) with the one of prolonging coherence times as a means of improving distance distribution resolution or prolonging distance range in pulsed dipolar EPR spectroscopy. To test the method, we used two bis-trityl rulers in which two trityl radicals are connected by a rather stiff linker. Linker length and residual flexibility are known (Godt et al., 2006; Jeschke et al., 2010). The chosen trityl radical is structurally closely related to the Finland trityl radical and has similar EPR spectroscopic properties (Hintz et al., 2019). The narrow EPR spectrum of the used trityl radical makes it particularly amenable to single-frequency techniques for measurements of the dipole–dipole coupling (Reginsson et al., 2012) in a regime that is analogous to the one of homonuclear NMR experiments. Note that the sequence presented in this work relies on the narrow spectrum of the trityl radicals. We do not expect it to work with the much more commonly used nitroxide radicals.
The article is organized as follows: first, we review mathematically, in the language of the magnetic resonance community, what happens to all the interactions in the spin Hamiltonian if we apply a strong microwave field. In order to do this, we will introduce a nutating frame description. Then we explain how an appropriate phase-modulation scheme leads to “pulses” in the nutating frame. In the results section we show the synthesis of the bis-trityl rulers and present the application of a dressed-spin echo experiment to such rulers to measure the dipolar coupling between two trityl radicals.
We use the following convention for operators: no prime refers to the laboratory frame and one prime to the electron-spin rotating frame, i.e. the interaction frame with the Zeeman Hamiltonian of the electrons. Two primes refer to the nutating frame, which is obtained with an additional interaction frame transformation with the pulse Hamiltonian. We will usually only denote the Hamiltonian with primes and not all operators. If we mention axes in the text, we will explicitly use the primes, but we will omit them in mathematical formulas.
2.1 Averaging of interactions by strong continuous microwave irradiation
In order to understand the observations in this work, we need to study the influence of strong microwave irradiation on the different interactions present in the spin system. The spin Hamiltonian of a system with two coupled electrons () in a bath of nuclei is given in the rotating frame by
The first term is the microwave Hamiltonian, which is given in the electron-spin rotating frame by
The Rabi or nutation frequency is denoted by ω1, which depends on the microwave amplitude B1 and the gyromagnetic ratio of the electron, γe. We assume a constant microwave phase and neglect the influence of the microwaves on the nuclear spins. In the following, we will apply an interaction frame transformation (IAT) with to all other terms and use first-order average Hamiltonian theory to gain physical insight. The new frame is referred to as the nutating frame. The nutating frame Hamiltonian is based on spin operators for dressed electron spins and bare nuclear spins. For mathematical details please consult the Supplement.
If we choose the nutating frame frequency ωPM equal to the Rabi frequency, ωPM=ω1, the irradiation term is completely absorbed into the frame. In a real experiment with an ensemble of spins, ω1 will be distributed due to microwave inhomogeneities; thus, we will always have a remaining contribution of
The dressed-spin offset Ωd will be distributed over the sample, but as a molecule is by orders of magnitude smaller than the microwave wavelength. Ωd will be the same for all electron spins within one molecule.
As usual, the influence of a small g anisotropy and of an inhomogeneous static magnetic field B0 is captured in offset terms in the rotating frame
We neglect any tilt of the electron spin quantization axis due to strong g anisotropy, which is a good approximation for trityl and other organic radicals. The first-order average Hamiltonian after an IAT with vanishes, i.e.
In pulse EPR, the spectral width is often much larger than the Rabi frequency. In this case, the first-order approximation will be poor. It is, however, not a poor approximation for trityl radicals with our setup. For simplicity, we will mostly neglect the effect of resonance offsets ΩS, 1 and ΩS, 2.
The most important term in the context of distance measurements is the electron–electron coupling Hamiltonian, which contains dipolar and exchange (J) contributions
where μ0 is the vacuum permeability, μB is the Bohr magneton, g1 and g2 are the g factors of the two electron spins, and θ is the angle between the external magnetic field and the interspin vector with length r12. The exchange contribution is often but not always negligible in pulse EPR-based distance measurements. The prefactor of the dipolar coupling contains the distance information and is given by
This amounts to 52.04 MHz for r=1 nm. After transformation to the nutating frame, we obtain
The electron–electron dipolar coupling is not averaged to zero, but only scaled by a factor of . It is also tilted such that the unique axis of the coupling Hamiltonian points along the spin-lock axis (; in the NMR literature, often a tilted frame is used). In other words, the two dressed spins are still dipole–dipole coupled with half the original coupling strength and with an inverted sign of the interaction. This result is well-known in solid-state NMR (Rhim et al., 1970), where it is used to generate “magic echoes”. The isotropic J coupling is unaffected if both spins are irradiated. Note however that the difference of the resonance frequencies of the two dressed spins is much smaller than the one of the bare spins, as remarked upon already by Grzesiek and Bax (1995). The difference in relative magnitude of the exchange coupling and resonance frequency difference can lead to a different manifestation of the exchange coupling in the spectra. If both the dressed-spin offsets as well as the spin states of the two spins are the same, exchange coupling has no influence on the evolution. This is analogous to the situation of magnetically equivalent nuclei in liquid-state NMR. This different averaging of dipolar and exchange contributions might be exploited experimentally to distinguish the two contributions.
The term contains all electron–nucleus (hyperfine) couplings. If the Rabi frequency of the irradiation is much larger than all hyperfine couplings and nuclear Zeeman frequencies, this term also averages to zero in the nutating frame, i.e.
Equations (9) and (10) might appear to be irrelevant to distance measurements between electrons, but they are not. The terms and do not commute if nuclear–nuclear flip–flop terms are present, even if the hyperfine coupling is purely secular (no electron spin echo envelope modulation effect). For example, for the flip–flop terms in , . A simple spin echo sequence on the electron spins thus does not completely refocus the hyperfine coupling – the result is a dephasing of the electron spins, sometimes loosely referred to as “relaxation”. In principle, this dephasing stems from coherent evolution, but since the nuclear spin bath is usually very large, it is computationally very expensive to simulate a real system. Accordingly, most theoretical studies treat the internuclear couplings phenomenologically using effective flip rates (Klauder and Anderson, 1962). The situation during microwave irradiation of an electron–nuclear spin system has many parallels with heteronuclear decoupling in solid-state NMR (Ernst, 2003), where one distinguishes between the “real” transverse relaxation time due to incoherent dynamics, T2, and the effective relaxation time that is measured with a spin echo, T2′, and has large coherent contributions. Of course, in EPR, the coupling strengths and Rabi frequencies are several orders of magnitude higher than in NMR.
In EPR measurements of organic radicals at sufficiently low temperatures, usually at 50 K and below, the hyperfine and nuclear–nuclear couplings dominate the dephasing (Brown, 1979). In this case, averaging the hyperfine coupling to zero should drastically increase the dephasing time, because commutes with all remaining terms containing electron spin operators. At the same time, according to Eq. (8), the effective dipolar coupling is scaled by a factor . If the gain in dephasing time is larger than a factor of 2, it should – in principle – be possible to measure longer dipolar dephasing traces and thus longer distances. As we shall see later, the effective scaling factor may be even more favourable (), as the flip–flop terms in the electron–electron dipolar Hamiltonian may be truncated for bare spins but can be significant for dressed spins.
The immediate next question is then how one can measure the dipolar coupling during a spin-lock pulse. We propose using a phase-modulation scheme that we discuss in the next section.
2.2 Pulse dressed spin resonance with phase-modulated pulses
For simplicity and illustration, we first look at an isolated electron (spin 1∕2) in a static magnetic field B0 along the laboratory-frame z axis. If we irradiate this system with a linearly polarized electromagnetic field with frequency ωmw and amplitude 2B1, the Hamiltonian in angular frequency units is given by
with . We include an arbitrary phase ϕmw(t), which we will use later to generate dressed-spin PM pulse sequences. As usual, we now go into a rotating frame with frequency ωmw. If we neglect the time-dependent terms (rotating wave approximation, RWA), we obtain
with the offset , which is also used in Eq. (4). The main effect of the time-dependent terms is a Bloch–Siegert shift, i.e. just a small correction of ΩS. We can choose the PM as
with a modulation amplitude aPM, a modulation frequency ωPM, and a modulation phase ϕPM. The phase ϕ0 is what one would conventionally call the phase of the microwave pulse applied to the bare spins, i.e. for . Likewise, ϕPM is the phase of the PM pulse that is applied to the dressed spins. We use ϕ0=0 for the following discussion. For small modulation amplitudes, aPM≪1, we can use the approximations cos (ϕmw)≈1 and sin (ϕmw)≈ϕmw and obtain a truncated rotating-frame Hamiltonian
For a hard pulse, i.e. ω1≫ΩS, we can now apply a second interaction frame transformation with , use the RWA again, and obtain the dressed rotating frame Hamiltonian
with the dressed-spin offset , already introduced in Eq. (3), and a dressed-spin nutation (Rabi) frequency of ω1aPM∕2. Again, the RWA implies that we neglect a Bloch–Siegert shift, now for the dressed spins, which would introduce a correction to Ωd. The whole situation is analogous to the rotating frame Hamiltonian in Eq. (12), but with an exchange of axes.
Some words of caution: first, in EPR unlike in NMR, the hard pulse limit will often not be fulfilled. In a first step, one can use an interaction frame transformation with the whole effective nutation field, . For the sake of intuitive clarity, we will not do this for the qualitative discussion. Second, one can easily choose a large aPM, such that the RWA leading from Eqs. (14) to (15) is seriously invalid. This was recognized already in Grzesiek and Bax (1995) and studied separately in Laucht et al. (2016). In our study, imperfection of the RWA is visible in nutation curves, but the final results do not seem to be affected. The problem might be alleviated by using an appropriate frequency or amplitude modulation in order to generate a circularly polarized field in the rotating frame.
There are two alternatives to the phase-modulation schemes. One could equivalently formulate the dressed-spin resonance as a frequency modulation. Phase and frequency modulations are physically equivalent, but we prefer the phase modulation because the description of frequency modulation involves a time-dependent offset/detuning and thus a “wobbling” frame, which makes it harder to keep track of relative phases of coherences. Instead of any microwave/radio-frequency modulation, one could also use a modulation of the magnetic field along the laboratory frame z direction (Jeschke, 1999). Depending on the setup, the relative phase of the modulation can be locked to the phase of the driving field or not. If an arbitrary waveform generator setup is available, phase modulation may be preferable, as it does not require modulation coils and a radiofrequency amplifier and makes synchronization of bare-spin and dressed-spin pulses much easier. However, the amplitude of the phase pulses depends on the Rabi frequency itself in the case of phase modulation. By using an external oscillating field, this dependence would vanish.
2.3 Pulse sequence
The pulse sequence used to measure the dipolar coupling in this work is the dressed-spin primary echo sequence shown in Fig. 1. It can be readily understood with results from the previous sections. For dipolar measurements, one chooses τ1=τ2 and constant TSL. The first π∕2 pulse generates electron coherence. Since we deal with trityl radicals, the excitation can be nearly uniform on our setup. The spins are then locked with a spin-lock pulse that is 90∘ phase shifted with respect to the first pulse. Let us assume that the spin lock and the coherences are along x′. For free dressed-spin evolution, i.e. in the absence of phase modulation, we can assume the following Hamiltonian during the spin lock in the nutating frame:
where we recall that . Note that ω1 is inhomogeneous over the sample but is the same within each pair of spins.
The Hamiltonian in Eq. (16) is analogous to the one in the rotating frame, but to a very good approximation the offsets are the same for both dressed electron spins. Additionally, all hyperfine couplings vanish. The phase-modulation pulses act on the dressed spins in the nutating frame. We can thus generate a dressed-spin echo by a phase–pulse sequence . A third π∕2 pulse is needed that rotates any refocused dressed-spin coherence back to the axis. The magnetization resulting from this backrotation is locked again until it is detected by the remaining echo sequence. A very similar sequence was already demonstrated with z-modulation pulses (Jeschke, 1999), albeit not for dipole–dipole coupled electron spins.
The dressed-spin echo is needed to refocus microwave field inhomogeneities (i.e. a distribution of ω1 and thus also Ωd). The dipolar part of the Hamiltonian is unaffected by the PM-π pulse, because this pulse inverts both spin operators at the same time. With effects of the other terms being refocused, it is sufficient to only keep the dipolar part during the periods τ1 and τ2:
At the start of the period τ1, the system is in the state (or along y′′, depending on the phase ϕ2). For τ1=τ2, this evolves according to
The z′′ terms are then flipped to , are transferred to bare-spin coherence at the end of the microwave pulse, and are then detected by the echo. The other terms do not contribute to the detected signal. The factor of 3∕4 has two contributions. A factor of is due to the spin lock and the partial averaging of the dipolar coupling; see Eq. (8). A factor of 3∕2 is due to the strong coupling regime in the dressed frame, because the dressed electron spins are equivalent. This scaling by a factor of 3∕2 for trityl biradicals has been observed before at short distances with established single-frequency techniques (Meyer et al., 2018), where it results from the dipole–dipole coupling being much larger than the mean difference of the resonance frequencies of two trityl radicals. In conclusion, we expect that the dressed spin echo intensity oscillates with 3∕4 of the dipolar coupling, which for a fixed or narrowly distributed distance will manifest in a Pake pattern because the measurements are conducted in frozen solution.
The timing τSL of the read-out echo does, in principle, affect the resulting dipolar spectrum, because it acts as a filter with the signal intensity scaling with cos ωddτSL. However, for short interpulse delays and long distances, such filtering should be negligible. If necessary, a SIFTER-type read-out sequence could be used, which refocuses both the offsets and the dipolar couplings.
It is noteworthy that, in principle, a normal two-pulse echo on the bare spins with non-selective pulses would be sufficient to measure the dipolar coupling. In practice, this approach is usually much inferior to the DQC and SIFTER sequences, because the phase memory time is of the same order of magnitude as the dipolar oscillations, echo decay is not monoexponential and contains other contributions, and dead time is significant. The combination of these complications makes it very difficult to separate the dipolar oscillation. Under the spin lock, the relaxation is sufficiently slowed down, such that the dipolar evolution is clearly distinguishable, and the dead time in a PM-pulse sequence is nearly zero. If the dead time becomes too large for the relevant dipolar oscillations, one could, in principle, apply the known dead-time free single-frequency pulse sequences DQC and SIFTER also as a phase–pulse sequence in the nutating frame.
2.4 Expected limitations
The derivation of the modulation formula in Eq. (18) depends on the condition that ω1 is much bigger than all other frequencies present in the system. Especially for the bare-spin resonance offsets, this approximation is not fulfilled very well. In principle, one could account for the different offsets analytically, but this is rather tedious and does not provide much additional insight. We will present numerical simulations in the results section to illustrate the deviations.
All measurements were performed on a home-built Q-band spectrometer equipped with a Keysight M8190A arbitrary waveform generator operating at 8 GS s−1 and an ADC with a sampling frequency of 2 GHz (SP Devices ADQ412) (Doll, 2016). The highly flexible software made it straightforward to implement the pulse sequences with PM pulses, in contrast to commercial analogues. Microwave pulses were amplified with a travelling wave tube (TWT) amplifier with 150 W nominal output power (Applied Systems Engineering). A home-built Q-band loop-gap resonator for 1.6 mm tubes was used (Tschaggelar et al., 2017). Note that the long spin-lock pulses with full power can be dangerous for the receiver, since much of the power is reflected by the overcoupled resonator. We did not take any special measures beyond the receiver protection switch (Doll and Jeschke, 2017). However, we are planning to install an additional limiter or a slow switch that could take more power. Since the spin-lock pulses are rather long, a slow switch could be used for most of the time, while the fast switch could be used for the transient times of the pulses to still provide the small dead time.
As model compounds, we used bis-trityl rulers with electron–electron distances of about 4.1 and 5.3 nm. The synthesis is discussed in Sect. 4.1. The bis-trityl rulers were dissolved in ortho-terphenyl (OTP) or its perdeuterated analogue dOTP providing solutions of different concentrations. More details are given in each figure and the Supplement.
Measurements were conducted at 50 K using a liquid helium flow cryostat. We did not systematically test the optimal temperature for each measurement. However, it is likely that higher temperatures would allow for shorter shot repetition times without dramatically changing the dephasing times.
Frequency-domain spectra were measured with chirp echoes and subsequent Fourier transform instead of field sweeps (Doll and Jeschke, 2014). Chirp pulses covered a range of 300 MHz symmetrically around the centre of the spectrum. The powder spectrum was simulated with the EasySpin library (Stoll and Schweiger, 2006).
The two-pulse dephasing time Tm was measured with a sequence echo with tπ=2tπ. Different pulse lengths were used to check whether instantaneous diffusion contributes to coherence loss. Similarly to previous findings by Meyer et al. (2018), it was found that flip angles of π∕2 or 3π∕2 for the second pulse gave higher echo intensities than an angle of π. More details are given in the Supplement.
The rotating frame relaxation time T1ρ was measured with the sequence in Fig. 1 in the absence of any phase-modulation pulses and variable TSL and with τSL=200 ns. Interestingly, T1ρ is significantly different when measured with a simple spin-locked echo with the sequence lock echo. More details are given in the results section and the Supplement.
The rotating frame relaxation time T2ρ for the mono-trityl was measured with the sequence in Fig. 1 including the phase-modulation pulses and fixed TSL and with τSL=200 ns. In the case of the bis-trityls, it is impossible to measure T2ρ independently of the dipolar coupling. Where applicable, we mention the decay rate of the “intramolecular background” for comparison.
All decay rates were obtained by fitting a stretched exponential of the functional form
to the relaxation curves, where t=2τ and T=Tm for the two-pulse echo decay, t=2τ1 and T=T2ρ for dressed echo decays, and t=TSL and T=T1ρ for the longitudinal rotating frame relaxation time.
The Rabi frequency ω1 was measured with a nutation experiment echo. As a control, we performed a dressed-spin resonance experiment with the sequence in Fig. (4) but only one PM pulse with low amplitude and variable frequency. This also yields the ω1 spectrum (see the Supplement). A similar experiment with z modulation was demonstrated in Jeschke (1999).
When the Rabi spectrum is known, one can set the value of the PM frequency ωPM. One then needs to choose a value for the modulation amplitude aPM and set up the PM pulse lengths. This can be achieved with a PM nutation experiment. Again, one uses the basic sequence in Fig. 1, with one pulse only with now fixed ωPM. One then observes the echo intensity as a function of the PM pulse length. That way the optimal PM pulse length can be determined. When choosing aPM=0.3, we observed only slight Bloch–Siegert shift-related oscillations in the PM pulse nutation traces while achieving a PM π-pulse length of 40–42 ns.
The dressed echo cannot be detected directly, because τ1 and τ2 are both indirect variables. Only the actual echo at the end of the microwave pulse sequence is digitized continuously. In order to optimize indirect detection, we checked that the last PM pulse in Fig. 1 is applied at the correct position. We observed that the position seems to be nearly perfectly predictable by setting , where tπ∕2 refers to the length of the PM-π∕2 pulse. We observed crossing dressed-spin echoes when changing interpulse delays in the PM pulse sequence, similar to what is known in microwave multi-pulse sequences in pulse EPR. Interestingly, the position of some unwanted echoes depends on the choice of τ0. Nevertheless, all these unwanted echoes can be suppressed by phase cycling the initial phases ϕPM of the PM pulses, ϕ2−4.
A step-by-step guide to setting up the sequence is provided in the Supplement.
The synthesis of bis-trityl rulers 1 and 2 is presented in Fig. 2. They were assembled from the rod-like building blocks 6 equipped with amino groups at both ends and trityl acid chloride 8. The latter was prepared from the corresponding trityl acid 7 (also named mono-trityl) using a procedure that has been described for the corresponding conversion of the structurally related Finland trityl radical (Shevelev et al., 2014). To achieve a complete conversion of the building blocks 6, trityl acid chloride 8 was used in excess. Leftover trityl acid chloride 8 is hydrolysed upon workup and the resulting trityl acid 7 is easily removed by filtration through silica gel. The building blocks 6 were obtained through a sequence of alkynyl–aryl coupling reactions (Sahoo et al., 2010; Qi et al., 2016; Ritsch et al., 2019) and a final oxidative alkyne dimerization. Oxidative alkyne dimerization is a very efficient way to obtain rod-like spacers with the same functional groups at both ends. Although this gives a butadiyne moiety, the spacer is still rather stiff and therefore the spin–spin distance sufficiently well-defined (Godt et al., 2006; Jeschke et al., 2010).
4.2 Relaxation of mono-trityl 7
As a reference, we measured the spectrum and the relaxation properties of the mono-trityl 7 in dOTP; see Fig. 3. As visible in panel b, T1ρ is orders of magnitude larger than Tm. Unfortunately, our TWT prevents us from using spin-lock pulses of more than 40 µs, meaning that uncertainty in T1ρ is rather large. Nevertheless, fitting a single stretched exponential to each curve yields values of Tm=2.9 µs and T1ρ≈930 µs. As mentioned above, the distance measurements based on dressed-spin echoes are limited by the transverse rotating frame relaxation time T2ρ rather than the longitudinal one T1ρ. The blue curve in panel b shows the dressed echo decay, indicating that , with a fitted value of T2ρ=13.1 µs.
While conceptually simple, the large difference between T1ρ and T2ρ was rather surprising to us. We are not aware of any example in the literature where T2ρ is discussed in-depth in the context of EPR, although there are several discussions in NMR and MRI (Michaeli et al., 2004). It remains unclear what the limiting contribution to T2ρ is. In analogy to solid-state NMR, residual coupling terms of the hyperfine interactions certainly contribute. An additional contribution would be the remaining intermolecular dipolar couplings, but then we would expect a strong dependence on the concentration, which we did not observe. Another factor that will definitely contribute is the noise of the driving field (Cohen et al., 2017). The noise (phase and amplitude) of the TWT during spin lock will not be refocused by the dressed echo. It is hard to quantify this contribution, since we do not have high-power amplifiers with different noise figures. In the future, we might investigate the influence of artificially added driving noise on T2ρ.
The large difference between T1ρ and T2ρ is unfortunate, because our proposed sequence will be limited by the latter. Nevertheless, one might come up with a sequence that will be limited by the former, longer relaxation time, and thus we measured T1ρ also for the bis-trityl rulers.
4.3 Bis-trityl 1, r≈4.1 nm
The results for bis-trityl 1 are shown in Fig. 4. The chirp echo FT-EPR spectrum is shown in panel a. The spectrum consists of a slightly asymmetric line with an FWHM of 16 MHz. The theoretical excitation profile of a 4 and an 8 ns microwave pulse are overlaid, showing that the whole spectrum can be excited almost uniformly with rectangular pulses.
The relaxation measurements for Tm and T1ρ are displayed in panel b, and they show the same trends as in the case of the mono-trityl. Note that the Tm measurement displayed was done with 100/200 ns pulses. Otherwise, the dipolar oscillations are already strongly visible in the two-pulse echo decay. It is immediately clear that the rotating frame relaxation time T1ρ is much longer than the phase-memory time, T1ρ≫Tm. The phase memory time is about 3.3 µs, while after 40 µs of spin lock, the echo intensity is still more than 90 % of its maximal value. A naive fit with a stretched exponential yields T1ρ≈560 µs.
The modulation of the dressed-spin echo is displayed in panel c. Clear oscillations are visible in the primary data. Since we do not currently have a model for the background, we fitted a stretched exponential to the data. This background is very similar to the T2ρ measurement of the mono-trityl (14.3 µs vs. 13.1 µs decay constant), which also means that it decays much faster than T1ρ. Note that not only intermolecular dipolar couplings from remote spins contribute to the background. Transverse relaxation of dressed spins with time constant T2ρ also contributes because we do not perform a constant time experiment. Regarding modulation depth we would have expected it to be unity, which is clearly not seen in our experiments. We suspect that imperfections in the dressed-spin π pulse lead to an unmodulated background, which cannot be removed by phase cycling. The phenomenon is similar to reduced instantaneous diffusion for a Hahn echo if the flip angle of the π pulse is reduced.
After background correction by division and a Fourier transform, we obtain the spectrum in Fig. 4d. The spectrum is a nice Pake pattern with the characteristic singularities at 1 and 2 times the dressed-spin dipolar frequency. The singularities appear at the expected positions. The splitting parameter d can be calculated from the expected distance of 4.1 nm, but it is scaled by a factor of 3∕4 as discussed above.
4.4 Bis-trityl 2, r≈ 5.3 nm
The analogous data of bis-trityl 2 are displayed separately in Fig. 5. The chirp echo FT EPR spectrum looks essentially the same as for bis-trityl 1, with the same slight asymmetry and an FWHM of 16 MHz.
The two-pulse microwave echo decay is slightly faster for bis-trityl 2 (2.9 µs vs. 3.3 µs). Again, it is difficult to really quantify a decoherence time that is not influenced by residual echo envelope contributions from intramolecular electron–electron coupling. Even with 100/200 ns pulses the excitation profile of the π pulse is still larger than the dipolar coupling, and some dipolar contribution to the echo envelope function is expected. The signal decay of dressed-spin polarization under the spin lock (T1ρ=730 µs) is again much slower and comparable to the case of bis-trityl 1. Unfortunately, the dipolar oscillations in panel c are not as clear as in the case of shorter distances. Also, the background is already rather fast compared to the dipolar frequencies (14.5 µs decay constant). In the dipolar spectrum, panel d, it becomes clear that this case is more complicated, because additional singularities appear at around . These features must result from the breakdown of some approximation that we have made in our theoretical description. Most likely they are due to the finite strength of the spin lock compared to the inhomogeneous spectral width. For two spins with different bare-spin resonance offsets, both the direction and magnitude of the effective field in the rotating frame differ. Accordingly, the two dressed spins have different resonance frequencies and quantization axes. Unless the dipole–dipole coupling is much larger than the frequency difference, it is significantly perturbed. In order to give a more quantitative explanation, we will show simplified numerical simulations in the following.
4.5 Numerical simulations
In order to understand the deviation of our experimental results from the theoretical expectation based on first-order average Hamiltonian theory (especially in the case of bis-trityl 2), we performed simplified numerical simulations. In principle, one could simulate the complete sequence, including the time-dependent phase during the phase pulses. We chose a simplified route: we start with both spins along z′ and then calculate the expectation value of during the spin lock using the Hamiltonian
In order to refocus the nutation of the spins around the effective field, we invert the phase of the irradiation in the middle of the spin lock, such that
This emulates the effect of the dressed refocusing (phase) pulse. With this choice, the evolution consists of two periods with time-independent Hamiltonians, which is straightforward to calculate on a computer.
In our implementation, which is available online, the parameters Ω1, Ω2, r, and θ are drawn in Monte Carlo fashion from their respective distributions (Gaussian for the first three, P(θ)=sin (θ) with for the latter). Statistical independence of the parameters is assumed. It is not unlikely that this assumption is at least partially wrong, since the respective orientation of the trityl moieties is restricted by the rigid linker. Although we have implemented simulations with a distance distribution, we do not consider such cases here, but rather assume fixed values of r. Additionally, all the simulations shown here assume on-resonance irradiation in the sense that the mean values of Ω1 and Ω2 are 0.
Some illustrative simulations are shown in Fig. 6. For each parameter set, we display the numerical simulation in the time and frequency domain as solid lines and show the analytical dipolar powder pattern (scaled by 3∕4) as dashed lines on top. In panel a, we show simulations assuming infinitely narrow EPR lines. In this case, the numerical and analytical results are the same. Panel b shows a simulation where we assume a FWHM of 16 MHz for both offset distributions (denoted by ΓΩ). For the case of r=5.3 nm, the simulation qualitatively reproduces the experimental results for bis-trityl (2), especially regarding the singularities in the dipolar spectra. For r=4.1 nm, the experimental results actually look better than the simulation if one regards the additional singularities at as an artifact. In this case, simulations with ΓΩ=8 MHz are actually closer to the experimental results (see panel c). This might suggest that the difference in offsets of bis-trityl (1) is smaller than the EPR spectrum might suggest. Either hyperfine and dipolar couplings significantly contribute to the linewidth of the EPR spectrum, or the offsets are not completely uncorrelated in reality.
In order to guide future developments, we also simulated traces assuming ΓΩ=8, but with significantly larger microwave strengths of MHz; see panel d. Compared to panels b and c, these simulations already show much better defined dipolar spectra. In conclusion, the simulations confirm that at least some of the artifacts in the experimental results are due to the finite size of the electron spin nutation frequency. The contribution of the artifacts becomes larger for larger offset differences and smaller dipolar couplings.
We showed that it is possible to measure the dipolar coupling between trityl radicals during a spin lock by using short intervals of phase modulations, i.e. by a dressed-spin echo generated with PM pulses. The relaxation during the spin lock is much slower compared to a simple two-pulse echo decay. The phenomena can be conceptually understood by describing the spin lock in a nutating frame and using average Hamiltonian theory. For an electron–electron distance of ≈4.1 nm, the experimental spectra agree very well with the theoretical expectations that assume a microwave Rabi frequency much larger than all other interactions in the system. For a distance of ≈5.3 nm, additional singularities appear in the dipolar spectrum. While the spin dynamics underlying these additional contributions can be understood by numerical simulations, they might seriously complicate data analysis in terms of distance distributions and have to be addressed in the future if the sequence should be used in application work. Additionally, we showed a profound difference between the longitudinal and transverse rotating frame relaxation times, T1ρ and T2ρ. In our case, the latter is much smaller than the former and unfortunately limits the distance measurements by the sequence introduced here. Preliminary results with the OX063 trityl and its partially deuterated analogue OX71 in different solvent compositions (not shown) revealed that even bare-spin relaxation at low temperatures and low concentrations is complicated to understand, let alone dressed-spin relaxation with characteristic times T2ρ and T1ρ. We are planning to investigate this in more detail and to compare the different relaxation times also at different temperatures. Note that in dOTP, the Tm values of the slow relaxing component of nitroxides (the relaxation of nitroxides in dOTP can be described by a sum of two stretched exponentials) can still be bigger than the T2ρ times measured here for trityl radicals (Soetbeer et al., 2018).
Since there are still significant artifacts present in the dipolar spectra when measuring longer distances, we refrained from a systematic analysis of signal-to-noise ratio and a comparison with existing pulse sequences.
Nevertheless, we are confident that the presented obstacles can be overcome. First, it might very well be possible to come up with a dressed pulse sequence that measures the dipolar coupling with an observation time limited by T1ρ instead of T2ρ. This appears feasible because, unlike the sum of dressed-spin polarizations of the two spins, their difference is affected by dipolar coupling. This fact is used in cross-polarization in solid-state NMR and oscillatory behaviour of magnetization transfer in the rotating frame has been studied in the context of heteronuclear correlation spectroscopy (Müller and Ernst, 1979). Second, the ratio of Rabi frequency to offsets could be reduced by going to a lower field. While in principle we could have done the experiments at X-band frequencies, our TWT in this range can only generate pulses of up to 15 µs. The Rabi frequencies generated by our setup are already rather high (≈100 MHz compared to ≈50 MHz in most commercial setups), but several groups around the world are working on micro resonators (Anders and Lips, 2019; Sidabras et al., 2019; Narkowicz et al., 2008; Blank et al., 2017), which generally give higher conversion factors and could be used to generate higher Rabi frequencies. If these difficulties can be overcome, pulse dressed electron paramagnetic resonance could significantly expand the measurable distance range, at least for trityl radicals.
The supplement related to this article is available online at: https://doi.org/10.5194/mr-1-75-2020-supplement.
NW designed the research with input from GJ. HH synthesized the bis-trityl rulers under the supervision of AG. NW carried out all measurements, data analysis, and simulations with help from AV. The EPR part of the manuscript was written by NW and edited by GJ and AG. The synthetic part was written by HH and AG.
The authors declare that they have no conflict of interest.
We thank Lukas Schreder, who carried out a research project in a very independent manner and sparked NW's interest in dressed electron spin resonance. We thank Jan Henrik Ardankjær-Larsen for providing OX063 and OX071 for preliminary relaxation studies. Matthias Ernst is acknowledged for helpful discussions and critically reading the theoretical part of the article.
This research has been supported by the ETH Zürich (grant no. ETH-48 16-1) and the Deutsche Forschungsgemeinschaft (SPP 1601 (grant no. GO 555/6-2)).
This paper was edited by Robert Bittl and reviewed by Thomas Prisner and Jack H. Freed.
Borbat, P. P., Georgieva, E. R., and Freed, J. H.: Improved sensitivity for long-distance measurements in biomolecules: Five-pulse double electron-electron resonance, J. Phys. Chem. Lett., 4, 170–175, https://doi.org/10.1021/jz301788n, 2013. a
Bordignon, E. and Bleicken, S.: New limits of sensitivity of site-directed spin labeling electron paramagnetic resonance for membrane proteins, BBA-Biomembranes, 1860, 841–853, https://doi.org/10.1016/j.bbamem.2017.12.009, 2018. a
Breitgoff, F. D., Soetbeer, J., Doll, A., Jeschke, G., and Polyhach, Y. O.: Artefact suppression in 5-pulse double electron electron resonance for distance distribution measurements, Phys. Chem. Chem. Phys., 19, 15766–15779, https://doi.org/10.1039/C7CP01488K, 2017. a
Brown, I. M.: Electron spin-echo studies of relaxation processes in molecular solids, in: Time Domain Electron Spin Resonance, edited by: Kevan, L. and Schwartz, R., John Wiley and Sons, Inc, New York, p. 200, 1979. a
Chen, H.-Y. and Tycko, R.: Slice Selection in Low-Temperature, DNP-Enhanced Magnetic Resonance Imaging by Lee-Goldburg Spin-Locking and Phase Modulation, J. Magn. Reson., 313, 106715, https://doi.org/10.1016/j.jmr.2020.106715, 2020. a
Cohen-Tannoudji, C., Dupont-Roc, J., and Grynberg, G.: Atom-Photon Interactions, Wiley-VCH, Berlin, 1992. a
Duss, O., Yulikov, M., Jeschke, G., and Allain, F. H.-T.: EPR-aided approach for solution structure determination of large RNAs or protein-RNA complexes, Nat. Commun., 5, 3669, https://doi.org/10.1038/ncomms4669, 2014. a
Georgieva, E. R., Ramlall, T. F., Borbat, P. P., Freed, J. H., and Eliezer, D.: The Lipid-binding Domain of Wild Type and Mutant a-Synuclein, J. Biol. Chem., 285, 28261–28274, https://doi.org/10.1074/jbc.m110.157214, 2010. a
Grzesiek, S. and Bax, A.: Audio-Frequency NMR in a Nutating Frame. Application to the Assignment of Phenylalanine Residues in Isotopically Enriched Proteins, J. Am. Chem. Soc., 117, 6527–6531, https://doi.org/10.1021/ja00129a016, 1995. a, b, c, d, e
Hintz, H., Vanas, A., Klose, D., Jeschke, G., and Godt, A.: Trityl Radicals with a Combination of the Orthogonal Functional Groups Ethyne and Carboxyl: Synthesis without a Statistical Step and EPR Characterization, J. Org. Chem., 84, 3304–3320, https://doi.org/10.1021/acs.joc.8b03234, 2019. a, b
Jeschke, G.: Coherent superposition of dressed spin states and pulse dressed electron spin resonance, Chem. Phys. Lett., 301, 524–530, https://doi.org/10.1016/S0009-2614(99)00041-X, 1999. a, b, c, d, e, f
Jeschke, G.: DEER Distance Measurements on Proteins, Annu. Rev. Phys. Chem., 63, 419–46, https://doi.org/10.1146/annurev-physchem-032511-143716, 2012. a
Jeschke, G., Pannier, M., Godt, A., and Spiess, H.: Dipolar spectroscopy and spin alignment in electron paramagnetic resonance, Chem. Phys. Lett., 331, 243–252, https://doi.org/10.1016/S0009-2614(00)01171-4, 2000. a
Jeschke, G., Sajid, M., Schulte, M., Ramezanian, N., Volkov, A., Zimmermann, H., and Godt, A.: Flexibility of shape-persistent molecular building blocks composed of p-phenylene and ethynylene units, J. Am. Chem. Soc., 132, 10107–10117, https://doi.org/10.1021/ja102983b, 2010. a, b
Kulik, L. V., Dzuba, S. A., Grigoryev, I. A., and Tsvetkov, Y. D.: Electron dipole-dipole interaction in ESEEM of nitroxide biradicals, Chem. Phys. Lett., 343, 315–324, https://doi.org/10.1016/S0009-2614(01)00721-7, 2001. a
Laucht, A., Simmons, S., Kalra, R., Tosi, G., Dehollain, J. P., Muhonen, J. T., Freer, S., Hudson, F. E., Itoh, K. M., Jamieson, D. N., McCallum, J. C., Dzurak, A. S., and Morello, A.: Breaking the rotating wave approximation for a strongly driven dressed single-electron spin, Phys. Rev. B, 94, 1–5, https://doi.org/10.1103/PhysRevB.94.161302, 2016. a, b
Laucht, A., Kalra, R., Simmons, S., Dehollain, J. P., Muhonen, J. T., Mohiyaddin, F. A., Freer, S., Hudson, F. E., Itoh, K. M., Jamieson, D. N., McCallum, J. C., Dzurak, A. S., and Morello, A.: A dressed spin qubit in silicon, Nat. Nanotechnol., 12, 61–66, https://doi.org/10.1038/nnano.2016.178, 2017. a
Meyer, A., Jassoy, J. J., Spicher, S., Berndhäuser, A., and Schiemann, O.: Performance of PELDOR, RIDME, SIFTER, and DQC in measuring distances in trityl based bi- and triradicals: Exchange coupling, pseudosecular coupling and multi-spin effects, Phys. Chem. Chem. Phys., 20, 13858–13869, https://doi.org/10.1039/c8cp01276h, 2018. a, b
Michaeli, S., Sorce, D. J., Idiyatullin, D., Ugurbil, K., and Garwood, M.: Transverse relaxation in the rotating frame induced by chemical exchange, J. Magn. Reson., 169, 293–299, https://doi.org/10.1016/j.jmr.2004.05.010, 2004. a
Milikisyants, S., Scarpelli, F., Finiguerra, M. G., Ubbink, M., and Huber, M.: A pulsed EPR method to determine distances between paramagnetic centers with strong spectral anisotropy and radicals: The dead-time free RIDME sequence, J. Magn. Reson., 201, 48–56, https://doi.org/10.1016/j.jmr.2009.08.008, 2009. a
Milov, A. D., Ponomarev, A. B., and Tsvetkov, Y. D.: Electron-electron double resonance in electron spin echo: Model biradical systems and the sensitized photolysis of decalin, Chem. Phys. Lett., 110, 67–72, https://doi.org/10.1016/0009-2614(84)80148-7, 1984. a
Narkowicz, R., Suter, D., and Niemeyer, I.: Scaling of sensitivity and efficiency in planar microresonators for electron spin resonance, Rev. Sci. Instrum., 79, 084702, https://doi.org/10.1063/1.2964926, 2008. a
Pannier, M., Veit, S., Godt, A., Jeschke, G., and Spiess, H. W.: Dead-Time Free Measurement of Dipole-Dipole Interactions between Electron Spins, J. Magn. Reson., 142, 331–340, https://doi.org/10.1006/jmre.1999.1944, 2000. a
Qi, M., Hülsmann, M., and Godt, A.: Spacers for Geometrically Well-Defined Water-Soluble Molecular Rulers and Their Application, J. Org. Chem., 81, 2549–2571, https://doi.org/10.1021/acs.joc.6b00125, 2016. a
Reginsson, G. W., Kunjir, N. C., Sigurdsson, S. T., and Schiemann, O.: Trityl radicals: Spin labels for nanometer-distance measurements, Chem.-Eur. J., 18, 13580–13584, https://doi.org/10.1002/chem.201203014, 2012. a
Ritsch, I., Hintz, H., Jeschke, G., Godt, A., and Yulikov, M.: Improving the Accuracy of Cu(II)-Nitroxide RIDME in the Presence of Orientation Correlation Evaluated with Water-soluble Cu(II)-Nitroxide Rulers, Phys. Chem. Chem. Phys., 21, 9810–9830, https://doi.org/10.1039/C8CP06573J, 2019. a, b
Sahoo, D., Thiele, S., Schulte, M., Ramezanian, N., and Godt, A.: Polar tagging in the synthesis of monodisperse oligo(p-phenyleneethynylene)s and an update on the synthesis of oligoPPEs, Beilstein J. Org. Chem., 6, 20–24, https://doi.org/10.3762/bjoc.6.57, 2010. a
Saiko, A. P., Fedaruk, R., and Markevich, S. A.: Suppression of electron spin decoherence in Rabi oscillations induced by an inhomogeneous microwave field, J. Magn. Reson., 290, 60–67, https://doi.org/10.1016/j.jmr.2018.02.003, 2018. a
Schmidt, T., Wälti, M. A., Baber, J. L., Hustedt, E. J., and Clore, G. M.: Long Distance Measurements up to 160 Å in the GroEL Tetradecamer Using Q-Band DEER EPR Spectroscopy, Angew. Chem. Int. Edit., 55, 15905–15909, https://doi.org/10.1002/anie.201609617, 2016. a
Shevelev, G. Y., Krumkacheva, O. A., Lomzov, A. A., Kuzhelev, A. A., Rogozhnikova, O. Y., Trukhin, D. V., Troitskaya, T. I., Tormyshev, V. M., Fedin, M. V., Pyshnyi, D. V., and Bagryanskaya, E. G.: Physiological-temperature distance measurement in nucleic acid using triarylmethyl-based spin labels and pulsed dipolar EPR spectroscopy, J. Am. Chem. Soc., 136, 9874–9877, https://doi.org/10.1021/ja505122n, 2014. a
Sidabras, J. W., Duan, J., Winkler, M., Happe, T., Hussein, R., Zouni, A., Suter, D., Schnegg, A., Lubitz, W., and Reijerse, E. J.: Extending electron paramagnetic resonance to nanoliter volume protein single crystals using a self-resonant microhelix, Science Advances, 5, eaay1394, https://doi.org/10.1126/sciadv.aay1394, 2019. a
Soetbeer, J., Hülsmann, M., Godt, A., Polyhach, Y., and Jeschke, G.: Dynamical decoupling of nitroxides in o-terphenyl: a study of temperature, deuteration and concentration effects, Phys. Chem. Chem. Phys., 20, 1615–1628, https://doi.org/10.1039/C7CP07074H, 2018. a
Spindler, P. E., Waclawska, I., Endeward, B., Plackmeyer, J., Ziegler, C., and Prisner, T. F.: Carr-Purcell Pulsed Electron Double Resonance with Shaped Inversion Pulses, J. Phys. Chem. Lett., 6, 4331–4335, https://doi.org/10.1021/acs.jpclett.5b01933, 2015. a
Tschaggelar, R., Breitgoff, F. D., Oberhänsli, O., Qi, M., Godt, A., and Jeschke, G.: High-Bandwidth Q-Band EPR Resonators, Appl. Magn. Reson., 48, 1273–1300, https://doi.org/10.1007/s00723-017-0956-z, 2017. a
Ward, R., Bowman, A., Sozudogru, E., El-Mkami, H., Owen-Hughes, T., and Norman, D. G.: EPR distance measurements in deuterated proteins, J. Magn. Reson., 207, 164–167, https://doi.org/10.1016/j.jmr.2010.08.002, 2010. a
Wili, N.: Distance measurements between trityl radicals by pulse dressed electron paramagnetic resonance with phase modulation: Raw Data, Processing Scripts, Simulations (Version submitted version), Data set, Zenodo, https://doi.org/10.5281/zenodo.3703053, 2020. a