# na-zfilter-amplitude-p2half-12-cycling-p2.in # approach: 12-phase cycling the first pulse, the second pulse and the receiver # spin-3/2 echo and antiecho amplitude optimization # versus the second-pulse duration # in three-pulse z-filtered amplitude-modulated 3Q-MAS sequence, # the first-pulse duration p1 = 4 micro seconds # the third-pulse duration p3 = 2 micro seconds spinsys { channels 23Na nuclei 23Na quadrupole 1 1 1e6 1 0 0 0 } par { spin_rate 5000 variable tsw 0.25 sw 1.0e6/tsw np 17 crystal_file rep10 gamma_angles 10 proton_frequency 800e6 start_operator 2*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}\ { 60 0 90}\ {120 0 270}\ {180 0 90}\ {240 0 270}\ {300 0 90}\ { 0 180 90}\ { 60 180 270}\ {120 180 90}\ {180 180 270}\ {240 180 90}\ {300 180 270}} { 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 23Na: " puts [resfreq 23Na $par(proton_frequency)] } # SIMP # NP=17 # SW=4000000 # TYPE=FID # DATA # 0 0 # -0.740892714 9.99200722e-16 # -1.26310063 6.66133815e-16 # -1.57285154 7.77156117e-16 # -1.70312398 7.77156117e-16 # -1.69172165 1.22124533e-15 # -1.57380998 1.08246745e-15 # -1.38648084 9.679757e-16 # -1.1741252 9.43689571e-16 # -0.985056587 5.55111512e-16 # -0.859144239 3.33066907e-16 # -0.814950495 2.49800181e-16 # -0.846110081 -1.52655666e-16 # -0.929975788 -2.11636264e-16 # -1.04264292 -3.60822483e-16 # -1.17036647 -5.82867088e-16 # -1.31087705 -6.66133815e-16 # END