# zfilter_p1S.in # spin-x/2 central-transition echo and antiecho # amplitude optimization (x = 3, 5, 7 or 9) # with the first pulse p1 in z-filtered MQMAS sequence, # OQ -> ±pQ -> OQ -> -1Q # all the coherences belonging to the same MQ coherence transfer pathway # are considered # author: R. HAJJAR spinsys { channels 27Al nuclei 27Al quadrupole 1 2 5e6 1 0 0 0 } par { spin_rate 5000 variable tsw 0.25 sw 1.0e6/tsw np 17 crystal_file rep100 gamma_angles 10 proton_frequency 800e6 start_operator I1z verbose 1101 variable rf 90000 variable rf3 9300 variable p2 4.0 variable p3 9.0 } proc pulseq {} { global par maxdt $par(tsw) #============= SPIN 3/2 ============= #----- ±3Q coherences ----- # matrix set 1 coherence {-3 3} #----- 0Q coherence ----- # matrix set 2 coherence {0} #----- -1Q coherence of the central-line ----- # matrix set detect elements {{2 3}} #============= SPIN 5/2 ============= #----- ±3Q coherences ----- matrix set 1 coherence {-3 3} #----- ±5Q coherences ----- # matrix set 1 coherence {-5 5} #----- 0Q coherence ----- matrix set 2 coherence {0} #----- -1Q coherence of the central-line ----- matrix set detect elements {{3 4}} #============= SPIN 7/2 ============= #----- ±3Q coherences ----- # matrix set 1 coherence {-3 3} #----- ±5Q coherences ----- # matrix set 1 coherence {-5 5} #----- ±7Q coherences ----- # matrix set 1 coherence {-7 7} #----- 0Q coherence ----- # matrix set 2 coherence {0} #----- -1Q coherence of the central-line ----- # matrix set detect elements {{4 5}} #============= SPIN 9/2 ============= #----- ±3Q coherences ----- # matrix set 1 coherence {-3 3} #----- ±5Q coherences ----- # matrix set 1 coherence {-5 5} #----- ±7Q coherences ----- # matrix set 1 coherence {-7 7} #----- ±9Q coherences ----- # matrix set 1 coherence {-9 9} #----- 0Q coherence ----- # matrix set 2 coherence {0} #----- -1Q coherence of the central-line ----- # matrix set detect elements {{5 6}} #========================================================== acq for {set i 1} {$i < $par(np)} {incr i} { pulse $par(tsw) $par(rf) x store 2 filter 1 pulse $par(p2) $par(rf) x filter 2 pulse $par(p3) $par(rf3) x acq -y reset prop 2 } } proc main {} { global par fsave [fsimpson] $par(name).fid puts "Larmor frequency (Hz) of 27Al: " puts [resfreq 27Al $par(proton_frequency)] }