# Bi-cog140(0,10,1,71).in # spin-9/2 anti echo and echo amplitude optimization # versus the first-pulse duration # in three-pulse amplitude-modulated z-filter 7QMAS, # the second-pulse duration p2 = 4 micro seconds # the third-pulse duration p3 = 2 micro seconds 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 N 140 variable wA 0 variable wB 10 variable wC 1 variable wRec 71 } proc pulseq {} { global par maxdt $par(tsw) acq -x 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 for {set j 0} {$j < $par(N)} {incr j} { set par(phA) [expr $j*$par(wA)*360./$par(N)] set par(phB) [expr $j*$par(wB)*360./$par(N)] set par(phC) [expr $j*$par(wC)*360./$par(N)] set par(phREC) [expr $j*$par(wRec)*360./$par(N)] 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 #2.34355457e-09 8.20143642e-09 #1.6490509e-09 9.07336339e-09 #-1.65426246e-07 1.02934568e-08 #-1.02707119e-05 1.65200711e-08 #-0.00018705453 1.07328959e-07 #-0.00166694176 8.19750758e-07 #-0.00934973644 3.9964337e-06 #-0.0383235178 1.33189303e-05 #-0.124828184 3.36938159e-05 #-0.336393408 7.11500505e-05 #-0.766014137 0.000136484498 #-1.50177413 0.000253959624 #-2.58981003 0.00047187132 #-4.01914543 0.000864282518 #-5.73449142 0.00151077698 #-7.65903325 0.00245690995 #END