# al-zfilter-amplitude-24-cycling-p3.in # approach: 24-phase cycling the first pulse, the third pulse and the receiver # spin-5/2 echo and antiecho amplitude optimization # versus the third-pulse duration # in three-pulse z-filtered amplitude-modulated 3Q-MAS sequence, # the first-pulse duration p1 = 4 micro seconds # the second-pulse duration p2 = 4 micro seconds spinsys { channels 27Al nuclei 27Al quadrupole 1 2 3e6 0 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 1.5*I1z detect_operator I1c verbose 1101 variable rf 90000 variable rf3 93000 variable p1 4 variable p2 4 } proc pulseq {} { global par maxdt $par(tsw) acq $par(ph31) pulse $par(p1) $par(rf) $par(ph1) pulse $par(p2) $par(rf) x for {set i 1} {$i < $par(np)} {incr i} { pulse $par(tsw) $par(rf3) $par(ph3) acq $par(ph31) } } proc main {} { global par foreach p {{ 0 0 270}\ { 60 0 90}\ {120 0 270}\ {180 0 90}\ {240 0 270}\ {300 0 90}\ { 0 90 0}\ { 60 90 180}\ {120 90 0}\ {180 90 180}\ {240 90 0}\ {300 90 180}\ { 0 180 90}\ { 60 180 270}\ {120 180 90}\ {180 180 270}\ {240 180 90}\ {300 180 270}\ { 0 270 180}\ { 60 270 0}\ {120 270 180}\ {180 270 0}\ {240 270 180}\ {300 270 0}} { set par(ph1) [lindex $p 0] set par(ph3) [lindex $p 1] set par(ph31) [lindex $p 2] set g [fsimpson] if [info exists f] { fadd $f $g funload $g } else { set f $g } } fsave $f $par(name).fid funload $f puts "Larmor frequency (Hz) of 27Al: " puts [resfreq 27Al $par(proton_frequency)] } # SIMP # NP=17 # SW=4000000 # TYPE=FID # DATA # 0 0 # -1.38038835 0.00071214647 # -4.37482241 0.00379881024 # -9.32779806 0.010683085 # -15.0957419 0.0196573447 # -19.6477665 0.0253225347 # -21.0729772 0.0245022591 # -18.5176089 0.0210463118 # -12.664157 0.0222334276 # -5.53727452 0.0320310278 # 0.427219144 0.0489107043 # 3.59469252 0.0677671236 # 3.83551344 0.0816099222 # 2.33608939 0.0828011788 # 0.747159283 0.0671382676 # 0.27190484 0.0387126652 # 1.16957631 0.00849555533 # END