Clean Absorption Mode NMR Data Acquisition Based on Time-Proportional Phase Incrementation
Total Page:16
File Type:pdf, Size:1020Kb
Supporting Information
Clean absorption Mode NMR data acquisition based on time-proportional phase incrementation
Yibing Wu‡ • Arindam Ghosh‡ • Thomas Szyperski
Yibing Wu‡ • Arindam Ghosh‡ • Thomas Szyperski ( )
Department of Chemistry, The State University of New York at Buffalo, Buffalo, NY 14260,
E-mail: [email protected].
‡: These authors made equal contributions to this work and should both be considered as first authors.
S1 I. c+- sampling (Forward TPPI)
II. c- sampling (Backward TPPI)
IIII. Dual TPPI sampling III.1 cADD- sampling (Sum of forward and backward TPPI interferograms) III.2 cSUB- sampling (Difference of forward and backward TPPI interferograms) III.3 Measurement of phase shift III.4 Simulations to verify theory of dual TPPI data acquisition
IV. Generalization to multiple chemical shift evolution periods IV.1 Multidimensional NMR
(a) c+ c+ -sampling for K = 1
(b) c c -sampling for K = 1
ADD ADD (c) c c -sampling for K = 1
IV.2 RD projection NMR
(a) c+ c+ -sampling for K = 1
(b) c c -sampling for K = 1
(c) c+ c+ + c c-sampling for K = 1
(d) cADD cADD -sampling for K = 1
V. Comparative contour plots and cross sections
V.1 RD 2D CC(CON)HN recorded with forward TPPI in RD dimension
V.2 RD 2D CC(CON)HN recorded with ‘double’ dual TPPI in RD dimension
S2 In this Supporting Information, we describe in detail the theory of TPPI sampling schemes alluded to in the paper, using the following definition to characterize the schemes: the time evolution of chemical shift can be sampled with N time domain data point (indexed by q = 1,2,3…,N) employing Time
Proportional Phase Incrementation (TPPI) as c (t) : cos( (q 1)t (q 1) ) or, equivalently, as , with representing a chemical shift, ‘±’ indicating c (t) : cos((q 1)t / 2 (q 1) ) forward (‘+’) and backward (‘’) sampling, denoting the TPPI phase increment and representing a phase shift which gives rise to a dispersive frequency domain peak component. In the following, equations are derived for the real time domain signal c±(t).
The real part of complex Fourier transformation (FT) of c±(t) (or equivalently a real FT) yields a frequency domain peak located at /t)(see below) and its quadrature image peak (‘quad peak’) located at /t), which is discarded(see below), both consisting of an absorptive peak component and a dispersive peak component. Here, represents the ‘TPPI frequency shift’. The four peak components can be written as elements of a vector F. With
A+ denoting the absorptive peak component of the peak at + , D+ denoting the dispersive peak component of the peak at + A denoting the absorptive peak component of the peak at , and D denoting the dispersive peak component of the peak at , the frequency domain signal is then proportional to
A D A D A D Re(F [c±(t)]) F= (S1). C A D
A D A D FC denotes the complex FT and = represents a ‘coefficient vector’ which provides the relative intensities of absorptive (A+) and dispersive (D+) peak components located at frequency + , as well as the relative intensities of absorptive (A) and dispersive (D) peak components located at the quad peak frequency – .
S3 I. c+ -sampling (Forward TPPI)
The forward TPPI sampled time domain signal is proportional to
S (t) c (t) cos(q 1)t (q 1) cos (q 1)t cos t t 1 ei t ei ei t ei 2
(S2),
where c+(t) is the discretely forward TPPI sampled interferogram, t = (q1)t =0,t,2t,3t,…,Nt and = /t denotes the “TPPI frequency shift”. With 1/t = 2•SW, where SW is the spectral width in Hz and = /2, the TPPI shift for the chemical shift is given by = •SW rad s-1. Eq. S2 yields, after discarding the modulations giving rise to the discarded quad peak:
1 1 i S (t) coseit sin e 2 ei t (S3). 2 2
After FT, one obtains for c+-sampling:
1 A cosΦ 2 (S4), 1 D sin Φ 2 yielding peak components A+ and D+ scaled with cos and sin, respectively.
Phase correction of the frequency domain spectrum requires both the real and the imaginary part. Complex FT of the time domain signal given in Eq. S2 is given by
1 i t i i t i 1 i 1 i FC e e e e e (A) i(D) e (A) i(D) 2 2 2 1 1 cos(A) sin (D) isin (A) i cos(D) cos(A) sin (D) isin (A) i cos(D) 2 2 1 i cos(A) sin (D) cos(A) sin (D) sin (A) cos(D) sin (A) cos(D) 2 2 R iI
(S5),
S4 where R and I are the real and the complex parts of the spectrum. The well known phase correction requires combining the real and imaginary parts according to cosR+sinI, which yields: cos R sin I 1 cos2 (A) cos sin (D) cos2 (A) cos sin (D) 2 1 sin 2 (A) sin cos (D) sin 2 (A) sin cos (D) (S6). 2 1 1 1 (A) cos 2(A) sin 2(D) 2 2 2
Hence, the dispersive peak component of the desired frequency domain peak at is eliminated while the dispersive peak component of the discarded quad peak is doubled. A first order phase correction requires that the ‘pivot point’ [2] is located at .
II. c -sampling (Backward TPPI)
The backward TPPI sampled time domain signal is proportional to
S (t) c (t) cos(q 1)t (q 1) cos (q 1)t cos t t 1 ei t ei ei t ei 2
(S7). which yields, after discarding the modulations giving rise to the discarded quad peak,
1 1 i S (t) coseit sin e 2 eit (S8). 2 2
After FT, one obtains for c-sampling:
1 A cosΦ 2 (S9), 1 D sin Φ 2
that is, the sign of the dispersive peak component D+ is inverted when compared with c+-sampling.
S5 III. Dual TPPI sampling
III.1 cADD-sampling (Sum of forward and backward TPPI interferograms)
Addition of St) in Eq. S3 and St) in Eq. S8 yields for such dual TPPI sampling
1 1 i 1 1 i SADD (t) S (t) S (t) cosei t sin e 2 ei t cosei t sin e 2 ei t 2 2 2 2 cosei t ei t cosei t
(S10).
After FT, one obtains for SADD (t) -sampling:
A ADD cosΦ D (S11), ADD 0 that is, the dispersive peak components D+ cancel and the absorptive peak component A+ is scaled by cos.
III.2 cSUB-sampling (Difference of forward and backward TPPI interferograms)
Subtraction of St) in Eq. S8 from St) in Eq. S3 yields for such dual TPPI sampling
1 1 i 1 1 i SSUB (t) S (t) S (t) cos ei t sin e 2 ei t cosei t sin e 2 ei t 2 2 2 2 i i sin e 2 ei t ei t sin e 2 ei t
(S12).
S6 After FT, one obtains for SSUB (t) -sampling:
A SUB 0 D (S13), SUB sin
that is, the absorptive peak components A+ cancel and the dispersive signal at D+ is scaled by sin. A zero order phase correction of /2, or a Hilbert transformation [1], converts D+ into A+, so that
A SUB sin D (S14). SUB 0
III.3 Measurement of phase shift
A A Inspection of Eqs. S11 and S14 reveals that measurement of the ratio of ADD and SUB , which can readily be measured by frequency domain signal integration, enables one to measure the phase shift according to
A SUB tan A ADD (S15).
For experiments designed to encode an NMR parameter in this enables one to measure .
III.4 Simulations to verify theory of dual TPPI data acquisition
Numerical simulations were performed using the program package MATLAB to verify equations provided in sections I-III. Constant amplitude time domain cosine modulations with oscillation frequency of 2,000 Hz and phase shift = 300 were generated. The TPPI frequency shift ( ) was chosen to be 6,000 Hz. Prior to performing complex FT, the time domain data was subjected to apodization using cosine bell window function.
Figure S1 shows the real part of the frequency domain spectrum obtained by complex FT of the simulated time domain signal. With a TPPI frequency shift of 6,000 Hz, the peak corresponding to frequency Hz appears at 8,000 Hz, while the quad peak appears at 8,000 Hz.
S7
Figure S1. Simulated spectra obtained by employing TPPI sampling with = 6,000 Hz. The desired and the quad peaks (on the left, indicated with ‘Quad peak’) appear at 8,000 and 8,000 Hz, respectively (see text).
In Figs, S1a and S1b, the spectra corresponding to forward TPPI sampling (Eq. S2) are shown without and with phase correction, respectively. As predicted by Eq. S6, the phase correction eliminates the phase of the desired peak while that of the quad peak is doubled.
Fig. S1c shows the spectrum obtained when using backward TPPI sampling (Eq. S7), where the phases of dispersive components of both the desired and the quad peak are reversed when compared to the forward TPPI sampled spectrum (Fig. S1a).
Fig. S1d shows the spectrum obtained after adding the forward and backward TPPI sampling spectra shown in Figs. S1a and S1c, where the dispersive peak components are eliminated (Eq. S10).
Fig. S1e shows the spectrum obtained after forming the difference subtraction of the forward and backward TPPI sampling spectra shown in Figs. S1a and S1c, where the absorptive peak components are eliminated (Eq. S12).
Fig. S1f shows the spectrum obtained after a 90o zero-order phase correction of the spectrum shown in Fig. S1e. The ratio of peak volumes of the desired (or quad) peak in Fig. S1e and those in this spectrum is given by tan (Eq. S15) and enables measurement of Ф.
S8 IV. Generalization to multiple chemical shift evolution periods
IV.1 Multidimensional NMR
TPPI sampling of indirect dimension can readily be generalized to K+1 indirect dimensions of a multi- dimensional NMR experiment, where K+1 chemical shifts 0, 1,…,K, that are associated with phase shifts 0, 1,…K, are sampled using either forward, backward or dual TPPI manner. For c±-sampling, the corresponding interferogram is obtained by K-fold tensor product formation:
K c (tK ) c (tK 1) c (t0 ) c (t j ) (S16). j0
In the following, the generalization of TPPI based sampling of multidimensional NMR spectra is derived for two indirect dimensions (K = 1), considering that the majority of practical applications will be for K = 1. Generalization to higher dimensionality (K > 1) is straightforward.
(a) c+ c+ –sampling for K = 1:
Starting with Eq. S3, the time domain signal for two indirect dimensions is proportional to
i i 1 i t i t 1 i t i t S (t ) S (t ) c (t ) c (t ) cos e 1 1 1 sin e 2 e 1 1 1 cos e 0 0 0 sin e 2 e 0 0 0 1 0 1 0 1 1 0 0 2 2 i 1 i t i t i t i t cos cos e 1 1 1 e 0 0 0 cos sin e 1 1 1 e 2 e 0 0 0 1 0 1 0 4 i i i i t i t i t i t sin cos e 2 e 1 1 1 e 0 0 0 sin sin e 2 e 1 1 1 e 2 e 0 0 0 1 0 1 0
(S17).
After FT, one obtains
Re(Fc (S (t1) S (t0 ))) 1 cos cos (A )(A ) cos sin (A )(D ) sin cos (D )(A ) sin sin (D )(D ) 4 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 (S18).
S9 The intensity of the desired peak component (A1+)(A0+), which is absorptive in both dimensions and is located at (+1+ 1 ,+ 0 ), is modulated with cos(cos(. In general, for K+1 dimensions the
cos j relative intensity of the desired peak component is modulated with j , where j represents the phase shift in the jth dimension. The other terms in Eq. S18 represent peak components which are dispersive in either one or in both dimensions.
(b) c c–sampling for K = 1:
Starting with Eq. S8, the time domain signal for two indirect dimensions is proportional to
i i 1 i t i t 1 i t i t S (t ) S (t ) c (t ) c (t ) cos e 1 1 1 sin e 2 e 1 1 1 cos e 0 0 0 sin e 2 e 0 0 0 1 0 1 0 1 1 0 0 2 2 i 1 i t i t i t i t cos cos e 1 1 1e 0 0 0 cos sin e 1 1 1e 2 e 0 0 0 1 0 1 0 4 i i i i t i t i t i t sin cos e 2 e 1 1 1e 0 0 0 sin sin e 2 e 1 1 1e 2 e 0 0 0 1 0 1 0
(S19).
After FT, one obtains
Re(Fc (S (t1) S (t0 ))) 1 cos cos (A )(A ) cos sin (A )(D ) sin cos (D )(A ) sin sin (D )(D ) 4 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 (S20).
(c) cADD cADD –sampling for K = 1:
Dual TPPI for both indirect dimensions requires all possible combinations of forward and backward TPPI sampling, that is, starting with Eq. S10, the time domain signal is proportional to
S10 ADD ADD ADD S (t1,t0 ) S (t1) S (t0 ) (c (t1) c (t1)) (c (t0 ) c (t0 ))
c (t1) c (t0 ) c (t1) c (t0 ) c (t1) c (t0 ) c (t1) c (t0 ) (S21).
c c c c c c c c
The first and the second terms in Eq. S21 are represented, respectively, by Eqs. S17 and S19. It is then straightforward to show that the time domain signal corresponding to the third term in Eq. S21 is proportional to
i i 1 i t i t 1 i t i t S (t ) S (t ) c (t ) c (t ) cos e 1 1 1 sin e 2 e 1 1 1 cos e 0 0 0 sin e 2 e 0 0 0 1 0 1 0 1 1 0 0 2 2 i 1 i t i t i t i t cos cos e 1 1 1e 0 0 0 cos sin e 1 1 1e 2 e 0 0 0 1 0 1 0 4 i i i i t i t i t i t sin cos e 2 e 1 1 1e 0 0 0 sin sin e 2 e 1 1 1 e 2 e 0 0 0 1 0 1 0
(S22).
After FT, one obtains
Re(Fc (S (t1) S (t0 ))) 1 cos cos (A )(A ) cos sin (A )(D ) sin cos (D )(A ) sin sin (D )(D ) 4 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 (S23).
The time domain signal represented by the fourth term in Eq. S21 is proportional to
S11 i i 1 i t i t 1 i t i t S (t ) S (t ) c (t ) c (t ) cos e 1 1 1 sin e 2 e 1 1 1 cos e 0 0 0 sin e 2 e 0 0 0 1 0 1 0 1 1 0 0 2 2 i 1 i t i t i t i t cos cos e 1 1 1 e 0 0 0 cos sin e 1 1 1 e 2 e 0 0 0 1 0 1 0 4 i i i i t i t i t i t sin cos e 2 e 1 1 1 e 0 0 0 sin sin e 2 e 1 1 1 e 2 e 0 0 0 1 0 1 0
(S24).
After FT, one obtains
Re(Fc (S (t1) S (t0 ))) 1 cos cos (A )(A ) cos sin (A )(D ) sin cos (D )(A ) sin sin (D )(D ) 4 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 (S25).
Addition of the spectra represented by Eqs. S18, S20, S23 and S25 give rise to a spectrum exhibiting only the absorptive peak component scaled by cos1cos0.
In general, for K+1 dimensions the relative intensity of the absorptive peak component obtained using
cos multiple dual TPPI is proportional to j . j
IV.2 RD projection NMR
TPPI can be applied to an arbitrary sub-set of K+1 chemical shift evolution periods which are jointly sampled in RD NMR, where K+1 chemical shifts 0, 1,…K , which are associated with phase shifts 0,
1,…K, are measured as linear combinations (0+ 0 ±(11+ 1 ±…±KK+ K . The scaling factors j enable one to achieve different maximal evolution times for the different jointly sampled shifts, while proper adjustment of the TPPI frequency shifts enables one to shift peaks originating from j different linear combinations of chemical shifts to different spectral regions, thereby avoiding additional spectral overlap arising from the projection.
S12 Choosing c-- or c+-sampling for each of the shift evolution periods, the resulting interferogram can be obtained by K-fold tensor product formation according to:
K c ( K ,t) c ( K1,t) c ( 0 ,t) c ( j ,t) (S26), j0
where t t0 t1 /1 t2 / 2 t j / j tK / K . In the following, the generalization of TPPI based sampling of RD NMR spectra is derived for two jointly sampled chemical shifts, where 0 is detected in quadrature and 1 is encode in an in-phase splitting (K = 1), considering that the majority of practical applications will be for K = 1. Generalization to more than two jointly sampled chemical shifts (K > 1) is straightforward.
(a) c+ c+ –sampling for K = 1:
Starting with Eq. S3, the time domain signal is proportional to
1 i t i i t i 1 i t i i t i S ( ,t) S ( ,t) c ( ,t) c ( ,t) e 1 1 e 1 e 1 1 e 1 e 0 0 e 0 e 0 0 e 0 1 0 1 0 2 2
1 i t i( ) i t i( ) e 1 0 1 0 e 1 0 e 1 0 1 0 e 1 0 4 i t i( ) i t i( ) e 1 0 1 0 e 1 0 e 1 0 1 0 e 1 0
(S27).
FT reveals that the first two terms give rise to the chemical shift doublet with components located at
0+ 1 + 0 and 1 + 0 , with phase shifts 1+0 and 1+0, respectively. The locations of these peak components provide the desired linear combinations 0 and 0 of the jointly sampled chemical shifts 0 and The third and the fourth terms give rise to discarded quadrature image peaks. The phase associated with a given linear combination of chemical shifts is then given by an analogous linear combination of the associated phase shifts.
(b) c c –sampling for K = 1:
Starting with Eq. S8, the time domain signal is proportional to
S13 1 i t i i t i 1 i t i i t i S ( ,t) S ( ,t) c ( ,t) c ( ,t) e 1 1 e 1 e 1 1 e 1 e 0 0 e 0 e 0 0 e 0 1 0 1 0 2 2
1 i t i( ) i t i( ) e 1 0 1 0 e 1 0 e 1 0 1 0 e 1 0 4 i t i( ) i t i( ) e 1 0 1 0 e 1 0 e 1 0 1 0 e 1 0
(S28).
FT reveals that the same shift doublet components are obtained as with forward sampling (Eq. S27), while the phases are of opposite signs.
(c) c c+ c c –sampling for K = 1:
Addition of Eqs. S27 and S28 yields
S (1,t) S (0 ,t) S (1,t) S (0 ,t) c (1,t) c ( 0 ,t) c (1,t) c (0 ,t)
1 i t i( ) i t i( ) i t i( ) i t i( ) e 1 0 1 0 e 1 0 e 1 0 1 0 e 1 0 e 1 0 1 0 e 1 0 e 1 0 1 0 e 1 0 4 i t i( ) i t i( ) i t i( ) i t i( ) e 1 0 1 0 e 1 0 e 1 0 1 0 e 1 0 e 1 0 1 0 e 1 0 e 1 0 1 0 e 1 0
1 i t i t cos( )e 1 0 1 0 cos( )e 1 0 1 0 2 1 0 1 0
i10 1 0 t i10 1 0 t cos(1 0 )e cos(1 0 )e
(S29).
Inspection of Eq. S29 reveals that for obtaining a clean absorption mode RD NMR spectrum using dual TPPI, it is sufficient to add two spectra, i.e., one recorded with forward sampling of both chemical shift evolution periods and one recorded with backward sampling for both periods. However, the intensities of the peaks forming the peak doublet at 0+ 1 + 0 and 1 + 0 , are then scaled differently by cos(1+0) and cos(1+0), respectively.
Eq. S29 can be extended for an arbitrary number of projected chemical shifts N ,...,1 , where j j j with associated phase shifts N,N-1,…,1 yield
S (t, N ,...,0) S (t,N ,..., 0) 1 cos ... ei N N 1 ...1 0 t ei N N 1 ...1 0 t 2N N N 1 1 0
S14
(S30).
Inspection of Eq. S30 reveals that, in contrast to clean absorption mode GFT NMR spectra (devoid of quadrature peaks) acquired using the hypercomplex method [7], clean absorption mode RD experiments can be obtained by recording only two RD NMR spectra, independent of the number of jointly sampled chemical shift evolution periods.
ADD ADD (d) c c –sampling for K = 1:
The imbalance of intensities of the peak pairs forming the chemical shift doublet (Eq. S29) can be eliminated by using dual TPPI for both evolution periods. This requires recoding of all possible combinations of forward and backward sampling for the two periods. Starting with Eq. S10, the resulting time domain signal is proportional to
ADD ADD ADD S (1,0 ,t) S (1,t) S (0 ,t) (c (1,t) c (1,t)) (c (0 ,t) c (0 ,t))
c (1,t) c (0 ,t) c (1,t) c ( 0 ,t) c (1,t) c (0 ,t) c (1,t) c (0 ,t) (S31).
c c c c c c c c
The first and the second terms in Eq. S31 are represented, respectively, by Eqs. S27 and S28.
The time domain signal corresponding to the third term in Eq. S31 is proportional to
S15 1 i t i i t i 1 i t i i t i S ( ,t) S ( ,t) c ( ,t) c ( ,t) e 1 1 e 1 e 1 1 e 1 e 0 0 e 0 e 0 0 e 0 1 0 1 0 2 2
1 i t i( ) i t i( ) e 1 0 1 0 e 1 0 e 1 0 1 0 e 1 0 4 i t i( ) i t i( ) e 1 0 1 0 e 1 0 e 1 0 1 0 e 1 0
(S32), and the time domain signal corresponding to the fourth term in Eq. S31 is proportional to
1 i t i i t i 1 i t i i t i S ( ,t) S ( ,t) c ( ,t) c ( ,t) e 1 1 e 1 e 1 1 e 1 e 0 0 e 0 e 0 0 e 0 1 0 1 0 2 2
1 i t i( ) i t i( ) e 1 0 1 0 e 1 0 e 1 0 1 0 e 1 0 4 i t i( ) i t i( ) e 1 0 1 0 e 1 0 e 1 0 1 0 e 1 0
(S33).
Addition of two RD NMR spectra represented by Eqs. S32 and S33 yields a clean absorption mode spectrum. However, the intensity imbalance of the two components of the shift doublet is ‘swapped’ when compared to Eq. S29:
S (1,t) S (0 ,t) S (1,t) S (0 ,t) c (1,t) c (0 ,t) c (1,t) c (0 ,t)
1 i t i( ) i t i( ) i t i( ) i t i( ) e 1 0 1 0 e 1 0 e 1 0 1 0 e 1 0 e 1 0 1 0 e 1 0 e 1 0 1 0 e 1 0 4 i t i( ) i t i( ) i t i( ) i t i( ) e 1 0 1 0 e 1 0 e 1 0 1 0 e 1 0 e 1 0 1 0 e 1 0 e 1 0 1 0 e 1 0
1 i t i t cos( )e 1 0 1 0 cos( )e 1 0 1 0 2 1 0 1 0
i101 0 t i10 1 0 t cos(1 0 )e cos(1 0 )e
(S34).
Addition of the spectra represented by Eqs. S29 and S32 yields:
ADD i10 1 0 t i10 1 0 t i10 1 0 t i10 1 0 t S (1, 0 ,t) cos 1 cos 0 e e e e (S35). S16 FT reveals that using ‘double’ dual TPPI sampling yields a clean absorption mode RD NMR spectrum in which the components of a given chemical shift doublet have the same intensity and are scaled by cos1cos0. In general, for K+1 jointly sampled chemical shifts, the relative intensity is proportional to
cos j . j
If the quadrature detection of 0 is accomplished by using the hypercomplex method, the same results are obtained.
V. Comparative contour plots and cross sections
Contour plots of the RD 2D CC(CON)HN spectra alluded to in the paper are shown in this section.
V.1 RD 2D CC(CON)HN recorded with forward TPPI in RD dimension
S17 Figure S2. a Contour plot of spectrum obtained after FT, with the discarded quadrature peaks being indicated. The desired peaks are within the dashed red box. b Contour plot of desired peaks only, with peaks located at (13C)(13C). Residual dispersive components which cannot be removed by a zero- or first-order phase correction are apparent. The horizontal dashed horizontal line at 56 ppm indicates the carrier position for evolution of (13C). The cross section chosen for Fig. 3e is indicated by a vertical grey dashed line. The digital
1 13 resolutions in 2( H) and 1( C) are, respectively, 7.8 and 24 Hz/pt. V.2 RD 2D CC(CON)HN recorded with ‘double’ dual TPPI in RD dimension
S18 Figure S3. Same as in Fig. S2, except that ‘double’ dual TPPI sampling was used (see text). The elimination of dispersive components is apparent. The cross section chosen for Fig. 3f is indicated by a vertical grey dashed line.
S19