# Bi-5Q-referencephase100.in # spin-9/2 anti echo and echo amplitude optimization # versus the first-pulse duration # in three-pulse amplitude-modulated z-filter 5QMAS, # the second-pulse duration p2 = 4 micro seconds # the third-pulse duration p3 = 2 micro seconds # nested phase cycling with reference receiver phase spinsys { channels 209Bi nuclei 209Bi 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 10 variable NB 10 variable deltapA 5 variable deltapB -5 } 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(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 209Bi: " puts [resfreq 209Bi $par(proton_frequency)] } #SIMP #NP=17 #SW=4000000 #TYPE=FID #DATA #0 0 #4.61302917e-08 -3.3244744e-11 #1.87996659e-05 -1.32728332e-08 #0.00052904163 -3.77624033e-07 #0.00464067582 -3.47431777e-06 #0.0199315351 -1.66029355e-05 #0.0496937985 -5.15526566e-05 #0.0702562123 -0.000119548847 #0.00919378003 -0.00022651231 #-0.262294714 -0.000362491501 #-0.912864601 -0.000474415822 #-2.09324212 -0.000444723107 #-3.842429 -0.000122955456 #-6.00874027 0.000582298588 #-8.24982805 0.00161335347 #-10.1222539 0.00274894426 #-11.2006649 0.00367532504 #END