# shiftedantiecho_p1.in # spin-x/2 central-transition antiecho # amplitude optimization (x = 3, 5, 7 or 9) # with the first pulse p1 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 p2 1.75 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}} #========================================================== 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)] }