# al-zfilter-3Q-amplitude-24p2p3-filtering-p1.in # approach: filtering the desired coherences after the first pulse # and after the second pulse, three times # spin-5/2 echo and antiecho amplitude optimization # versus the first-pulse duration # in three-pulse z-filtered amplitude-modulated 3Q-MAS sequence, # the second-pulse duration p2 = 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 24*I1z detect_operator I1c verbose 1101 variable rf 90000 variable rf3 9300 variable p2 4 variable p3 2 variable nb1 0 variable nb2 0 } proc pulseq {} { global par maxdt $par(tsw) matrix set 1 coherence {-1 5} matrix set 2 coherence { 1 -5} matrix set 3 coherence {-3 3} matrix set 10 coherence {-4} matrix set 20 coherence {4} matrix set 30 coherence {0} acq -y for {set i 1} {$i < $par(np)} {incr i} { pulse $par(tsw) $par(rf) x store 2 filter $par(nb1) pulse $par(p2) $par(rf) x filter $par(nb2) pulse $par(p3) $par(rf3) x acq -y reset prop 2 } } proc main {} { global par for {set j 1} {$j < 4} {incr j} { set par(nb1) $j set par(nb2) [expr 10*$j] 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 # -2.42725708e-05 -6.42900018e-07 # 0.000273712464 -3.78464053e-07 # 0.00017614564 8.42643559e-06 # -0.00579187364 4.66168477e-05 # -0.0303434426 0.000132112772 # -0.0843402926 0.000258275216 # -0.163948528 0.000406007366 # -0.254569825 0.000582606448 # -0.341062215 0.000822116822 # -0.413162873 0.00113864158 # -0.472308564 0.00150509699 # -0.526914262 0.00187469212 # -0.57716916 0.00220325684 # -0.613424903 0.00245655438 # -0.626540505 0.00260664668 # -0.614468691 0.00263111477 # END