# na-cogwheel_24_0_4_19_7.in # approach: selection of the desired coherences (-3Q, 3Q) after the first pulse # and (0Q) after the second pulse with cogwheel phase cycling # N = 24, wA = 0, wB = 4, wC = 19, and wRec = 7 # Cog24(0, 4, 19, 7) # spin-3/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 spinsys { channels 23Na nuclei 23Na quadrupole 1 1 1e6 1 0 0 0 } par { spin_rate 5000 variable tsw 0.25 sw 1.0e6/tsw np 17 crystal_file rep100 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 24 variable wA 0 variable wB 4 variable wC 19 variable wRec 7 } 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 23Na: " puts [resfreq 23Na $par(proton_frequency)] } # SIMP # NP=17 # SW=4000000 # TYPE=FID # DATA # 0 0 # 3.91902455e-06 3.07393e-15 # 0.000117094148 2.69229083e-15 # 0.000795529255 2.22044605e-15 # 0.00287802458 6.66133815e-16 # 0.00725731211 -1.11022302e-15 # 0.0144295742 -5.55111512e-16 # 0.0242623665 -2.77555756e-15 # 0.0361384856 -2.55351296e-15 # 0.049323588 -4.21884749e-15 # 0.0632269656 -5.55111512e-15 # 0.0773291382 -4.10782519e-15 # 0.090853889 -4.21884749e-15 # 0.102514048 -3.99680289e-15 # 0.110643958 -3.16413562e-15 # 0.113745863 -3.3029135e-15 # 0.111140994 -2.87964097e-15 # END