# al-zfilter-5Q-amplitude-40p2p3-filtering-p3.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 third-pulse duration # in three-pulse z-filtered amplitude-modulated 5Q-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 40*I1z detect_operator I1c verbose 1101 variable rf 90000 variable rf3 9300 variable p1 4 variable p2 4 variable nb1 0 variable nb2 0 } proc pulseq {} { global par maxdt $par(tsw) matrix set 1 coherence {-1} matrix set 2 coherence {1} matrix set 3 coherence {-5 5} matrix set 10 coherence {4} matrix set 20 coherence {-4} matrix set 30 coherence {0} acq -y pulse $par(p1) $par(rf) x filter $par(nb1) pulse $par(p2) $par(rf) x filter $par(nb2) for {set i 1} {$i < $par(np)} {incr i} { pulse $par(tsw) $par(rf3) x acq -y } } 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 # -0.0611014363 2.64934828e-05 # -0.122301417 0.000106297699 # -0.183711003 0.000239698501 # -0.245337073 0.000427241591 # -0.307025807 0.000670161845 # -0.368513999 0.000970098579 # -0.429531686 0.00132806552 # -0.489873623 0.00174336513 # -0.549403583 0.00221341692 # -0.608011984 0.00273477192 # -0.665578167 0.00330450592 # -0.72197349 0.00392098789 # -0.777091582 0.00458376334 # -0.830866818 0.00529299093 # -0.883262364 0.00604903596 # -0.934240549 0.00685246243 # END