SPC52dlsw: 2D large DQ F1 spectral width SPC5 pulse program for TopSpin2.1

*** Outline ***

• Code for Avance III spectrometers with topSpin2.1 operating system
• Example: ³¹P in VPI-5 zeolite with AV700
• References

Code for Avance III spectrometers with topSpin2.1 operating system

;spc52dlsw (TopSpin 2.0)

;2D SQ-DQ correlation experiment with SPC5 sequence
;for the original C7 sequence see: Lee et al. Chem Phys Lett 242, 304-309, 1995
;see Hohwy, Rienstra, Jaroniec and Griffin, JCP 110, 7983, 1999
;using Mei Hong trick: J. Magn. Reson. 136, 86-91 (1999)
;by JOS 03/03/03, modified by HF 14.5.07

;Avance III version
;parameters:
;d1 : recycle delay

;p1 : f1 presaturation pulses and detection pulse at pl1

;pl2  : =120dB, not used
;pl11 : for SPC5 recoupling sequence B1=5*cnst31 in Hz

;cnst31 : spinning speed
;l0  : number of composite SPC5 mq excitation and reconversion cycles multiple of 5
;l3  : number of 2/5 rotor revolution increments for DQ evolution, usually 1 (why else would I use this?)
;l20 : # of pulses in saturation pulse train
;in0 : inf1=l3*(2s/cnst31)/5, t1 increment, multiples of 2/5th of a rotation period
;ns : 16*n
;FnMode : undefined
;mc2 : STATES-TPPI
;nd0 : 1
;WDW : F1 QSINE 3,  F2 QSINE 2 or EM
;zgoptns :-Dpresat or blank

;$COMMENT=SQ-DQ experiment with SPC5 sequence, optimised for large sweep width
;$CLASS=Solids
;$DIM=2D
;$TYPE=direct excitation
;$SUBTYPE=homonuclear correlation
;$OWNER=Bruker

define pulse pul360
  "pul360=(1s/cnst31)/5"    ;360° pulse
define pulse pul90
  "pul90=(0.25s/cnst31)/5"  ; 90° pulse
define pulse pul270
  "pul270=(0.75s/cnst31)/5" ;270° pulse

  "d31=(1s/cnst31)"

define loopcounter count    ;for STATES-TPPI procedure
  "count=td1/2"             ;and STATES cos/sin procedure

define delay del360
  "del360=(1s/cnst31)/5"
  "in0=l3*(2s/cnst31)/5"
  "inf1=l3*(2s/cnst31)/5"
  "l2=0"

define delay dwellf1
  "dwellf1=in0"

;cnst11 : to adjust t=0 for acquisition, if digmod = baseopt
"acqt0=1u*cnst11"

#include <rot_prot.incl>
                            ;protect for too slow rotation

  ze                        ;acquire into a cleared memory
1 d31
  dwellf1

#ifdef presat               ;set with -Dpresat
pres, d20                   ;delay between saturation pulses
  (p1 pl1 ph1):f1           ;saturation loop if required
  lo to pres times l20
#endif /* presat */

2 d1                        ;recycle delay
  1m rpp11                  ;reset the phase ph11 pointer to the first element
  1m rpp12                  ;reset the phase ph12 pointer to the first element
  1m rpp13                  ;reset the phase ph13 pointer to the first element
  1m rpp14                  ;reset the phase ph14 pointer to the first element

3 (pul90 pl11 ph11 ipp13 ipp14):f1    ;SPC5 DQ excitation
                            ;switch to SPC5 RF condition
                            ;increment reconversion phase ph13 and ph14 pointers
  pul360 ph12 ipp12         ;increment phase ph12 pointer
  pul270 ph11 ipp11         ;increment phase ph11 pointer
  lo to 3 times l0          ;l0 = multiple of 5 

                            ;double-quantum evolution period
4 del360 ipp13 ipp14        ;increment reconversion pulse phase ph13 and ph14 
  del360                    ;pointers to the next phase in the lists
  lo to 4 times l2          ;del360 = (1s/cnst31)/5

5 pul90  ph13               ;SPC5 DQ reconversion
  pul360 ph14 ipp14         ;increment phase ph14 pointer
  pul270 ph13 ipp13         ;increment phase ph13 pointer
  lo to 5 times l0          ;l0 = multiple of 5 

  (p1 pl1 ph5):f1           ;detection pulse, flip into the xy plane
  2u
  gosc ph31                 ;gosc does not loop to 1
                            ;start ADC with ph31 signal routing
                            ;DQ filtering (four phase cycling):
  40u ip13*16384            ;increments all phases of ph13 by 90°
  40u ip14*16384            ;increments all phases of ph14 by 90°
  lo to 1 times ns          ;next scan

  100m wr #0 if #0 zd       ;delay for disk I/O, store signal,
                            ;increase FID number
                            ;delete memory data
                            ;do not perform dummy scans
                            ;with next acquisition

  1m ip11*8192              ;increments all phases of ph11 by 45°,
                            ;90° phase for DQ coherence
  1m ip12*8192              ;increments all phases of ph12 by 45°,
                            ;90° phase for DQ coherence
  lo to 1 times 2           ;t1 quadrature detection

6 1m iu2                    ;increment counter l2 in 2*Tr/5
  lo to 6 times l3

  ;1m rp11                  ;reset all phases of ph11, ph12, ph13, and ph14 
  ;1m rp12                  ;to their original values, i.e. to the values they 
  ;1m rp13                  ;had before the first ip11, ip12, ip13, and ip14
  ;1m rp14                  ;in case of STATES remove semicolon at beginning of the 4 lines

  lo to 1 times count       ;count = td1/2

HaltAcqu, 1m
exit


ph1=   0                    ;for saturation pulse

ph11= (65536)     0 13107 26214 39322 52429 32768 45875 58982  6554 19661
ph12= (65536) 32768 45875 58982  6554 19661     0 13107 26214 39322 52429
ph13= (65536) 16384 29491 42598 55706  3277 49152 62259  9830 22937 36044
ph14= (65536) 49152 62259  9830 22937 36044 16384 29491 42598 55706  3277

ph5=   0 0 0 0 2 2 2 2 1 1 1 1 3 3 3 3
ph31 = 0 2 0 2 2 0 2 0 1 3 1 3 3 1 3 1    ;ph31 = ph5 + 2*ph13