# shiftedantiecho_p2.in # spin-x/2 central-transition antiecho # amplitude optimization (x = 3, 5, 7 or 9) # with the second pulse p2 in phase-modulated shifted-antiecho 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 15000 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 100000 variable rf3 10000 variable p1 4 variable p3 15 } proc pulseq {} { global par maxdt $par(tsw) #============= SPIN 3/2 ============= #----- 3Q coherence ----- # matrix set 1 coherence {3} #----- 1Q coherences ----- # matrix set 2 coherence {1} #----- -1Q coherence of the central-line ----- # matrix set detect elements {{2 3}} #============= SPIN 5/2 ============= #----- -3Q coherences ----- matrix set 1 coherence {-3} #----- 5Q coherence ----- # matrix set 1 coherence {5} #----- 1Q coherences ----- matrix set 2 coherence {1} #----- -1Q coherence of the central-line ----- matrix set detect elements {{3 4}} #============= SPIN 7/2 ============= #----- -3Q coherences ----- # matrix set 1 coherence {-3} #----- -5Q coherences ----- # matrix set 1 coherence {-5} #----- 7Q coherence ----- # matrix set 1 coherence {7} #----- 1Q coherences ----- # matrix set 2 coherence {1} #----- -1Q coherence of the central-line ----- # matrix set detect elements {{4 5}} #============= SPIN 9/2 ============= #----- -3Q coherences ----- # matrix set 1 coherence {-3} #----- -5Q coherences ----- # matrix set 1 coherence {-5} #----- -7Q coherences ----- # matrix set 1 coherence {-7} #----- 9Q coherence ----- # matrix set 1 coherence {9} #----- 1Q coherences ----- # matrix set 2 coherence {1} #----- -1Q coherence of the central-line ----- # matrix set detect elements {{5 6}} #========================================================== pulse $par(p1) $par(rf) x store 10 reset acq -y for {set i 1} {$i < $par(np)} {incr i} { reset prop 10 filter 1 pulse [expr ($i)*$par(tsw)] $par(rf) x filter 2 pulse $par(p3) $par(rf3) x acq -y } } proc main {} { global par fsave [fsimpson] $par(name).fid puts "Larmor frequency (Hz) of 27Al: " puts [resfreq 27Al $par(proton_frequency)] }