# Co-3Q-digitizedphase48.in # spin-7/2 anti echo and echo amplitude optimization # versus the first-pulse duration # in three-pulse amplitude-modulated z-filter 3QMAS, # the second-pulse duration p2 = 4 micro seconds # the third-pulse duration p3 = 2 micro seconds # nested phase cycling with digitized phase for receiver spinsys { channels 59Co nuclei 59Co 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 p2 4 variable p3 2 variable NA 6 variable NC 8 variable deltapA 3 variable deltapC 1 } proc pulseq {} { global par maxdt $par(tsw) acq -y for {set i 1} {$i < $par(np)} {incr i} { pulse $par(tsw) $par(rf) $par(phA) store 2 pulse $par(p2) $par(rf) $par(phB) pulse $par(p3) $par(rf3) $par(phC) acq [expr $par(phREC) - 90] reset prop 2 } } 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 59Co: " puts [resfreq 59Co $par(proton_frequency)] } #SIMP #NP=17 #SW=4000000 #TYPE=FID #DATA #0 0 #0.000560253803 -4.5893597e-07 #0.017895857 -1.23404915e-05 #0.127162105 -7.3849188e-05 #0.477281695 -0.000230260428 #1.24310804 -0.000492184077 #2.53059563 -0.000817261337 #4.28700234 -0.00113038157 #6.29189604 -0.00136871622 #8.24876923 -0.00151651041 #9.91385426 -0.00160285769 #11.1837145 -0.00167590739 #12.096944 -0.00177372045 #12.7642481 -0.00189287628 #13.2908759 -0.00196884029 #13.7335299 -0.00189322548 #14.0867933 -0.0015666893 #END