# al-spam-echo_-5Q-selective-20-cycling-p1.in # approach: 20-phase cycling the first pulse and the receiver # spin-5/2 antiecho amplitude optimization # versus the first-pulse duration # in three-pulse spam-echo 5Q-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 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 I1z detect_operator I1c verbose 1101 variable rf 90000 variable rf3 9300 variable p2 4 variable p3 2 variable NA 20 variable deltapA 5 } proc pulseq {} { global par maxdt $par(tsw) acq $par(phREC) for {set i 1} {$i < $par(np)} {incr i} { pulse $par(tsw) $par(rf) $par(phA) store 2 pulse $par(p2) $par(rf) x pulse $par(p3) $par(rf3) -x acq [expr $par(phREC) - 90] reset prop 2 } } proc main {} { global par 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)] 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.83786022e-010 2.22044605e-015 # 1.40611439e-007 8.8817842e-016 # 5.11944722e-006 1.66533454e-016 # 6.32019497e-005 5.03069808e-016 # 0.000427261325 -3.88578059e-016 # 0.00195823428 1.55431223e-015 # 0.00681986019 0 # 0.0193225289 6.66133815e-016 # 0.0465516878 1.77635684e-015 # 0.0983120906 1.55431223e-015 # 0.186047289 -6.66133815e-016 # 0.32074566 -8.8817842e-016 # 0.510274412 -1.11022302e-015 # 0.756883234 -1.77635684e-015 # 1.05563753 -2.22044605e-015 # 1.39426804 -1.77635684e-015 # END