# al-zfilter-amplitude-100-cycling-p3.in # approach: 100-phase cycling the first and the second pulses and the receiver # spin-5/2 echo and antiecho amplitude optimization # versus the first-pulse duration # in three-pulse z-filtered amplitude-modulated 5Q-MAS sequence, # the first-pulse duration p1 = 4 micro seconds # the third-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 rep10 gamma_angles 10 proton_frequency 800e6 start_operator I1z detect_operator I1c verbose 1101 variable rf 90000 variable rf3 93000 variable p1 4 variable p2 4 variable NA 10 variable NB 10 variable deltapA -5 variable deltapB 5 } proc pulseq {} { global par maxdt $par(tsw) acq $par(phREC) pulse $par(p1) $par(rf) $par(phA) pulse $par(p2) $par(rf) $par(phB) for {set i 1} {$i < $par(np)} {incr i} { pulse $par(tsw) $par(rf3) $par(phC) acq [expr $par(phREC) - 90] } } proc main {} { global par set par(phC) 0 for {set jB 0} {$jB < $par(NB)} {incr jB} { set par(phB) [expr $jB*360./$par(NB)] 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(deltapB)*$par(phB)] 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 # -1.09948769 0.000436225738 # -1.96013135 0.00139506572 # -2.42810097 0.00205625649 # -2.48781819 0.00155288552 # -2.25489229 -0.00086506078 # -1.90627592 -0.00591030116 # -1.5882431 -0.014089065 # -1.36123207 -0.0250674266 # -1.20815254 -0.0371071752 # -1.07955653 -0.0473163043 # -0.931503781 -0.0527760618 # -0.741483418 -0.0518421789 # -0.51497676 -0.0448978559 # -0.285198742 -0.0343340147 # -0.0939007841 -0.0238552521 # 0.0398455182 -0.017273945 # END