# na-spam-echo-selective-12-cycling-p2.in # approach: 12-phase cycling the first pulse and the receiver # spin-3/2 echo amplitude optimization # versus the second-pulse duration # in three-pulse spam-echo -3Q-MAS sequence, # the first-pulse duration p1 = 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 rep10 gamma_angles 10 proton_frequency 800e6 start_operator I1z detect_operator I1c verbose 1101 variable rf 90000 variable rf3 9300 variable p1 4 variable p3 2 variable NA 12 variable deltapA 3 } proc pulseq {} { global par maxdt $par(tsw) pulse $par(p1) $par(rf) $par(phA) store 10 reset acq $par(phREC) for {set i 1} {$i < $par(np)} {incr i} { reset prop 10 pulse [expr ($i)*$par(tsw)] $par(rf) x pulse $par(p3) $par(rf3) -x acq [expr $par(phREC) - 90] } } 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 23Na: " puts [resfreq 23Na $par(proton_frequency)] } # SIMP # NP=17 # SW=4000000 # TYPE=FID # DATA # 0 0 # 0.135236263 -5.55111512e-17 # 0.488908227 -1.66533454e-16 # 1.02313553 -1.66533454e-16 # 1.61231779 -4.4408921e-16 # 2.13335276 -4.4408921e-16 # 2.49863914 -1.66533454e-16 # 2.67167155 2.22044605e-16 # 2.66395758 3.33066907e-16 # 2.51966833 6.10622664e-16 # 2.29696575 4.99600361e-16 # 2.05252661 1.06165077e-15 # 1.83145293 9.43689571e-16 # 1.66180846 1.2490009e-15 # 1.55268571 6.66133815e-16 # 1.49580481 8.8817842e-16 # 1.47108365 8.32667268e-16 # END