;hncacbtr.fa
;3D TROSY-HN(CA)CB
;Yang and Kay, J. Am. Chem. Soc. 1999, 121, 2571-2575.
;Bruker Avance/Xwin-nmr version. This program requires XWIN-NMR 2.5+
;Written up by F. Abildgaard, NMRFAM (abild@nmrfam.wisc.edu)
;
;  $Id: hncacbtr.fa,v 1.3 2000/09/25 18:27:52 abild Exp abild $
;
;  Disclaimer: This pulse program is provided "as is" for your
;  information. Support for the use of this pulse program is only
;  provided to users of the National Magnetic Resonance Facility
;  at Madison (NMRFAM). Users of this pulse program employ it at
;  their own risk. Neither NMRFAM nor University of Wisconsin-Madison
;  are liable for any physical or other damage incurred during the
;  use of this pulse program.
;
;f1: 1H, f2: 15N, f3: 13C, f4: 2H (channel assignments may be changed below)
;o1p: 4.7 ppm
;o2p: 118 ppm,
;o3p: 43 ppm
;o4p: 4.5 ppm
;
;d1: relaxation delay = d1+0.1s
;p1 90 H at pl1
;p2 90 N at pl2
;p5 180 dgr. CO rectangular shaped (spnam5) at power sp5, offset 176-43,
;   semi-selective: field strength=dNu/sqrt(3),	dNu=(176ppm-43ppm)*bf3.
;   p5 43.2 us at 600 MHz.
;p6  90 dgr. Cab at pl6 for 90 dgr. pulse at highest field possible. 
;p7 180 dgr. Ca rectangular shaped (spnam7) at power sp7, offset 58-43 ppm,
;   semi-selective: field strength=dNu/sqrt(3),	dNu=(176ppm-58ppm)*bf3.
;   p7 48.6 us at 600 MHz.
;p8 180 dgr. shaped Cab refocusing REBURP (spnam8) pulse (non-selective for
;   aliphatic). Centered at 43 ppm (offset 0). p8 ~350us. Adjust phase
;   program ph22 for max. intensity.
;p9 180 dgr. shaped CO/Cg G3 (spnam9) pulse at power sp9, offset 150-43 ppm.
;   p9 300 us at 600 MHz.
;p11 90 H1 shaped (spnam1, EBURP-1) H2O pulse at power sp1, offset 0.
;   p11 ~ 7 ms. Check for a possible phase difference between hard H1 and
;   soft H1 pulses and set the phase program ph21 accordingly.
;p25 90 dgr. 2H pulse at pl15   
;pl0 120dB
;
;H2 Waltz-16 (cpdprg5) decoupling, 90 dgr. pulse (PCPD) at pl15
;N15 evolution:
;  in10=in30, SW(N)=1/(2*in10), typ. 30-40 ppm
;  l4 complex points; max. is (d10/in10)+1
;  Process as echo-antiecho.
;  Chemical shift axis is reversed.
;C13 evolution:
;  SW(C)=1/(2*in0), typ. 55 ppm
;  l6 complex points.
;  set cnst0 to 0 (preferably) or 1, 2, 3 .. to make d0 the smallest possible
;  positive delay. cnst0=0 gives (90,-180) phase distortion in F1.
;  cnst0=1 gives (270,-540) phase distortion in F1 (use LP to correct).
;  Chemical shift axis is reversed.
;ns=8, 16, ..., ds=8, 16,...
;
;Recommendations for gradients:
;gpz1:   7%
;gpy2:  30%
;gpy3: -44%
;gpz4:  30%
;gpx6:  35%
;gpz7:  15%
;gpx8:  54% adjust for magic angle
;gpz8:  30%
;gpx9:  44%
;gpy10: 30%
;gpx11: 54% adjust for magic angle
;gpz11: 30%
;gpnam1: sine.50
;gpnam2: sine.100
;gpnam3: sine.100
;gpnam4: sine.10
;gpnam6: sine.10
;gpnam7: sine.50
;gpnam8: sine.100
;gpnam9: sine.50
;gpnam10: sine.50
;gpnam11: sine.20
;
;Define one or more of the following options to tailor this pulse program
; to your specific needs.
;
;#define ONE_D		; uncomment for 1D experiment
#define N15_EVOL	; comment out for 2D w/o N15 evolution
#define C13_EVOL	; comment out for 2D w/o C13 evolution
;#define H2_DEC		; uncomment if 2H labeled sample
#define OPTIM_P19	; uncomment if you want to optimize p19 (GRAD11)
#define EXPTCORR	; uncomment if you want "expt" to report
;			; the correct expt time (works with XWIN-NMR 2.x)
;
;Define channel assignments:
#define H f1
#define N f2
#define C f3
#define D f5
;
;You shouldn't have to worry about anything beyond this point :-)
;
;sanity checks
;
#ifdef ONE_D
#undef N15_EVOL
#undef C13_EVOL
#endif
;
aqseq 321
;
define delay DELTA
define delay TAUA
define delay TAUA1
define delay TAUA9
define delay TAUA10
define delay TAUB
define delay TAUC
define delay TAUC4
define delay TAUC6
define delay TN
define delay TN2
define delay CEN_HN1
define delay CEN_CN1
define pulse H1_90
define pulse H1_180
define pulse H1_S90
define pulse N15_90
define pulse N15_180
define pulse CO_180
define pulse CAB_90
define pulse CA_180
define pulse CA_REBP
define pulse CO_G3
define pulse GRAD1
define pulse GRAD2
define pulse GRAD3
define pulse GRAD4
define pulse GRAD6
define pulse GRAD7
define pulse GRAD8
define pulse GRAD9
define pulse GRAD10
define pulse GRAD11

"d11=100m"			;disk i/o
"d12=10u"			;power switching etc.
"d13=5u"			;just a short delay
"d14=60u"			;ip,id etc
"d16=300u"			;gradient recovery
"d17=50u"			;short gradient recovery
"H1_90=p1"
"H1_180=H1_90*2"
"H1_S90=p11"
"N15_90=p2"
"N15_180=N15_90*2"
"CO_180=p5"
"CAB_90=p6"
"CA_180=p7"
"CA_REBP=p8"
"CO_G3=p9"
;
"GRAD1=400u"
"GRAD2=1.0m"
"GRAD3=750u"
"GRAD4=80u"
"GRAD6=80u"
"GRAD7=500u"
"GRAD8=3.000m"
"GRAD9=400u"
"GRAD10=400u"
#ifndef OPTIM_P19
"p19=305u"			; configurable: set to the optimum value on your instrument
#endif
"GRAD11=p19"

"TAUA=2.2m"			; configurable: ~1/(4*JNH)
"TAUA1=TAUA-GRAD1-d16-d13"
"TAUA9=TAUA-GRAD9-d16-d13"
"TAUA10=TAUA-GRAD10-d16-d13"
"TAUB=12.0m"			; configurable: usually 12ms
"TAUC=6.8m"			; configurable: usually 6.8-7.1ms for Cb only,
				; or 3.5 ms for Ca and Cb
#ifdef H2_DEC
  "TAUC4=TAUC-GRAD4-d17-p25-d12-d13*2"
  "TAUC6=TAUC-GRAD6-d17-p25-d12-d13*2"
#else
  "TAUC4=TAUC-GRAD4-d17-d13*2"
  "TAUC6=TAUC-GRAD6-d17-d13*2"
#endif
"TN=12.0m"			; configurable: usually 12ms
"d10=TN"
"d30=d13"
"TN2=TN-d13-GRAD8-d16-CO_180-d30"
"DELTA=GRAD11+d16+d13*2"
"CEN_HN1=N15_90-H1_90"
"CEN_CN1=(N15_180-CA_180)/2"

#ifdef C13_EVOL
  "d0=((cnst0*2+1)*in0-CAB_90*1.273-CO_G3-d12*2-d13*2)/2"
#endif

#ifdef EXPTCORR
  "d31=2*(GRAD1+TAUA1+TAUB+GRAD4+TAUC4+GRAD6+TAUC6+TAUA9+GRAD9+TAUA10+GRAD10)+GRAD2+GRAD3+GRAD7+GRAD8+TN2+DELTA+GRAD11"
#endif

#ifdef H2_DEC
#define H2_DEC_ON d12 pl15:D \n p25:D ph1 \n d13 cpds5:D
#define H2_DEC_OFF d13 do:D \n p25:D ph3 d12
#else
#define H2_DEC_ON d13
#define H2_DEC_OFF d13
#endif

#include <Avance.incl>
#include <Grad.incl>

1 ze
  d11 LOCKDEC_ON
2 d13
  d14 H2_LOCK
  d11 LOCKH_OFF
  3m
  d14
3 d14
  d14
  d14
  d14
  d14
  d14
4 d14
  d14
  d14
  d14
5 d14
  d14
  d14
6 d13
#ifdef EXPTCORR
#include <faexptcorr.incl>
#endif  
  d1 pl1:H pl2:N
  d13 LOCKH_ON
  d13 UNBLKGRAMP
  d13 H2_PULSE
; INEPT transfer from H to N
  (d13 d12 pl0 H1_S90:sp1 ph21 d13 d12 pl1):H	;Selective 90 H1 pulse phase -y
  (H1_90 ph0):H
  d13
  GRAD1:gp1		; 0.4ms, 5G/cm, z, sine.50
  d16
  TAUA1
  (CEN_HN1 H1_180 ph0):H (N15_180 ph0):N
  d13
  TAUA1
  GRAD1:gp1		; 0.4ms, 5G/cm, z, sine.50
  d16
  (H1_90 ph19):H	; phase y on Bruker DMX
  d13
  GRAD2:gp2		; 1.0ms, 10G/cm, z, sine.100
  d16 
; INEPT transfer from N to Ca
  (N15_90 ph10):N
  TAUB pl0:C
  (CEN_CN1 CA_180:sp7 ph0):C (N15_180 ph0):N
  TAUB
  (N15_90 ph1):N
  d13
  GRAD3:gp3		; 700us, 10G/cm, z, sine.100
  d16 pl6:C
  H2_DEC_ON
; Ca-Cb transfer
  (CAB_90 ph11):C
  TAUC4
  H2_DEC_OFF
  d13
  GRAD4:gp4		; 80us, 20G/cm, z, sine.10
  d17 pl0:C
  (CA_REBP:sp8 ph22):C  ; adjust ph22 for max. signal
  d13
  GRAD4:gp4		; 80us, 20G/cm, z, sine.10
  d17
  H2_DEC_ON
  TAUC4 pl6:C
  (CAB_90 ph12):C
; Begin C13 evolution
#ifdef C13_EVOL
  (d13 d12 pl0):C
  d0
  (CO_G3:sp9 ph0):C
  d0
  (d13 d12 pl6):C
#else
   d13
#endif
; End C13 evolution  
  (CAB_90 ph1):C
; Ca-Cb back-transfer
  TAUC6
  H2_DEC_OFF
  d13
  GRAD6:gp6		; 80us, 25G/cm, x, sine.10
  d17 pl0:C
  (CA_REBP:sp8 ph22):C  ; adjust ph22 for max. signal
  d13
  GRAD6:gp6		; 80us, 25G/cm, x, sine.10
  d17
  H2_DEC_ON
  TAUC6 pl6:C
  (CAB_90 ph0):C
  H2_DEC_OFF
  GRAD7:gp7		; 0.5ms, 10G/cm, z, sine.50
  d16 pl0:C
; Begin constant time evolution on N
  (N15_90 ph13):N
  d10
  (N15_180 ph14 TN2):N (CEN_CN1 CA_180:sp7 ph0):C
  d13
  GRAD8:gp8*EA*-1	; 3.000 ms, -/+30G/cm, ma, sine.100
  d16 pl0:C
  (CO_180:sp5 ph0):C
  d30
; End constant time
; Sensitivity enhanced coherence transfer from N to H
  (N15_90 ph15):N
  (H1_90 ph0):H
  d13
  GRAD9:gp9		; 0.4ms, 3.1G/cm, z, sine.50
  d16
  TAUA9
  (CEN_HN1 H1_180 ph0):H (N15_180 ph0):N
  d13
  TAUA9
  GRAD9:gp9		; 0.4ms, 3.1G/cm, z, sine.50
  d16
  (H1_90 ph1):H
  d13
  (N15_90 ph1):N
  d13
  GRAD10:gp10		; 0.4ms, 5.35G/cm, z, sine.50
  TAUA10
  d16
  (CEN_HN1 H1_180 ph0):H (N15_180 ph0):N
  d13
  TAUA10
  GRAD10:gp10		; 0.4ms, 5.35G/cm, z, sine.50
  d16
  (N15_90 ph0):N
  d13
  (H1_90 ph0):H
  DELTA
  (H1_180 ph0):H
  d13
  GRAD11:gp11		; 305 us, 30G/cm, ma, sine.20
  d13
  d16 BLKGRAMP

  go=2 ph31

#ifdef ONE_D
  d11 wr #0 H2_LOCK
#else
  d11 wr #0 if #0 zd H2_LOCK
#endif
  d13 LOCKH_OFF
#ifdef N15_EVOL
  3m igrad EA
  d14 ip15
  d14 ip15
  lo to 3 times 2
  d14 dd10
  d14 id30
  d14 ip13
  d14 ip13
  d14 ip31
  d14 ip31
  lo to 4 times l4
  d14 rd10
  d14 rd30
#else
  d14*10
  3m
#endif
#ifdef C13_EVOL
  d14 ip11
  d14 ip12
  lo to 5 times 2
  d14 id0
  d14 ip31
  d14 ip31
  lo to 6 times l6
#endif
  d14 LOCKDEC_OFF

exit

ph0=0
ph1=1
ph2=2
ph3=3
ph10=0 2
ph11=0 0 2 2
ph12={1}*4 {3}*4
ph13=1
ph14={0}*4 {2}*4
ph15=0
ph19=1			; phase y on Bruker DMX
ph21=(360) 270		; phase ph19+180
ph22=(360) 0		; phase x plus any phase diff.
ph31=0 2 2 0