# al-zfilter-amplitude-selective-40-cycling-p2.in # approach: 40-phase cycling the first and the third pulses and the receiver # spin-5/2 echo and antiecho amplitude optimization # versus the second-pulse duration # in three-pulse z-filtered amplitude-modulated 5Q-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 rep10 gamma_angles 10 proton_frequency 800e6 start_operator I1z detect_operator I1c verbose 1101 variable rf 90000 variable rf3 9300 variable p1 4 variable p3 2 variable NA 10 variable NC 4 variable deltapA -5 variable deltapC 1 } proc pulseq {} { global par maxdt $par(tsw) pulse $par(p1) $par(rf) $par(phA) store 10 reset acq $par(phREC) for {set i 1} {$i < $par(np)} {incr i} { reset prop 10 pulse [expr ($i)*$par(tsw)] $par(rf) $par(phB) pulse $par(p3) $par(rf3) $par(phC) acq [expr $par(phREC) - 90] } } proc main {} { global par set par(phB) 0 for {set jC 0} {$jC < $par(NC)} {incr jC} { set par(phC) [expr $jC*360./$par(NC)] for {set jA 0} {$jA < $par(NA)} {incr jA} { set par(phA) [expr $jA*360./$par(NA)] set par(phREC) [expr $par(deltapA)*$par(phA) + $par(deltapC)*$par(phC)] 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.00146850995 -3.21811744e-05 # -0.00449274115 -4.73841965e-05 # -0.0387916398 4.1912177e-06 # -0.129617481 0.000198123197 # -0.287204301 0.00058513071 # -0.48715967 0.00114387263 # -0.678990635 0.0017719667 # -0.812741983 0.00232321776 # -0.863192817 0.00267050368 # -0.836466046 0.00275536828 # -0.759305509 0.00259770728 # -0.662539156 0.00227351175 # -0.569492959 0.00188307865 # -0.491283248 0.0015152778 # -0.425713113 0.00120755583 # -0.359148753 0.0009253913 # END