# al-zfilter-phase-72-p1p2-cycling-p2.in # approach: 72-phase cycling the first pulse, the second pulse and the receiver # spin-5/2 echo amplitude optimization # versus the second-pulse duration # in three-pulse z-filtered phase-modulated +3Q-MAS sequence, # the first-pulse duration p1 = 4 micro seconds # the third-pulse duration p3 = 2 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 I1z detect_operator I1c verbose 1101 variable rf 90000 variable rf3 93000 variable p1 4 variable p3 2 } proc pulseq {} { global par maxdt $par(tsw) pulse $par(p1) $par(rf) $par(ph1) store 10 reset acq $par(ph31) for {set i 1} {$i < $par(np)} {incr i} { reset prop 10 pulse [expr ($i)*$par(tsw)] $par(rf) $par(ph2) pulse $par(p3) $par(rf3) x acq $par(ph31) } } proc main {} { global par foreach p {{ 0 0 270}\ { 30 0 180}\ { 60 0 90}\ { 90 0 0}\ {120 0 270}\ {150 0 180}\ {180 0 90}\ {210 0 0}\ {240 0 270}\ {270 0 180}\ {300 0 90}\ {330 0 0}\ { 0 60 90}\ { 30 60 0}\ { 60 60 270}\ { 90 60 180}\ {120 60 90}\ {150 60 0}\ {180 60 270}\ {210 60 180}\ {240 60 90}\ {270 60 0}\ {300 60 270}\ {330 60 180}\ { 0 120 270}\ { 30 120 180}\ { 60 120 90}\ { 90 120 0}\ {120 120 270}\ {150 120 180}\ {180 120 90}\ {210 120 0}\ {240 120 270}\ {270 120 180}\ {300 120 90}\ {330 120 0}\ { 0 180 90}\ { 30 180 0}\ { 60 180 270}\ { 90 180 180}\ {120 180 90}\ {150 180 0}\ {180 180 270}\ {210 180 180}\ {240 180 90}\ {270 180 0}\ {300 180 270}\ {330 180 180}\ { 0 240 270}\ { 30 240 180}\ { 60 240 90}\ { 90 240 0}\ {120 240 270}\ {150 240 180}\ {180 240 90}\ {210 240 0}\ {240 240 270}\ {270 240 180}\ {300 240 90}\ {330 240 0}\ { 0 300 90}\ { 30 300 0}\ { 60 300 270}\ { 90 300 180}\ {120 300 90}\ {150 300 0}\ {180 300 270}\ {210 300 180}\ {240 300 90}\ {270 300 0}\ {300 300 270}\ {330 300 180}} { set par(ph1) [lindex $p 0] set par(ph2) [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 # 0.0669092292 0.000161430688 # 0.730481963 0.000855611005 # 2.56927108 0.00201021433 # 5.44325004 0.00423557914 # 8.79504752 0.00691861288 # 12.2537555 0.00363133773 # 15.4486849 -0.0149298883 # 17.3723154 -0.048372061 # 16.4554277 -0.0800882972 # 11.7213591 -0.0889152325 # 3.92387627 -0.0675532705 # -4.65568572 -0.0271360207 # -11.6934966 0.0156177302 # -16.1180385 0.0524782629 # -18.0381578 0.0836126217 # -17.8665962 0.107853442 # END