# na-shifted-echo-amplitude-selective-24-cycling-p1.in # approach: 24-phase cycling the first pulse, the third pulse and the receiver # spin-3/2 echo and antiecho amplitude optimization # versus the first-pulse duration # in three-pulse shifted-echo amplitude-modulated 3Q-MAS sequence, # 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 } proc pulseq {} { global par maxdt $par(tsw) acq $par(ph31) for {set i 1} {$i < $par(np)} {incr i} { pulse $par(tsw) $par(rf) $par(ph1) store 2 pulse $par(p2) $par(rf) x pulse $par(p3) $par(rf3) $par(ph3) acq $par(ph31) reset prop 2 } } proc main {} { global par foreach p {{ 0 0 270}\ { 60 0 90}\ {120 0 270}\ {180 0 90}\ {240 0 270}\ {300 0 90}\ { 0 90 90}\ { 60 90 270}\ {120 90 90}\ {180 90 270}\ {240 90 90}\ {300 90 270}\ { 0 180 270}\ { 60 180 90}\ {120 180 270}\ {180 180 90}\ {240 180 270}\ {300 180 90}\ { 0 270 90}\ { 60 270 270}\ {120 270 90}\ {180 270 270}\ {240 270 90}\ {300 270 270}} { set par(ph1) [lindex $p 0] set par(ph3) [lindex $p 1] set par(ph31) [lindex $p 2] 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 # 6.8899869e-07 3.05311332e-16 # 2.12223522e-05 0 # 0.000151178414 6.66133815e-16 # 0.000581766162 6.66133815e-16 # 0.00157611006 0 # 0.00338016096 -4.4408921e-16 # 0.00610749622 -2.22044605e-16 # 0.00965282028 -1.33226763e-15 # 0.0136854856 -4.4408921e-16 # 0.0177409185 2.22044605e-16 # 0.0213780437 -9.99200722e-16 # 0.0243322783 -6.66133815e-16 # 0.0265915289 0 # 0.0283614981 -1.22124533e-15 # 0.029944974 -1.66533454e-16 # 0.0316025874 1.52655666e-16 # END