# al-split-t1-echo-96-phase-cycling-p2.in # oQ -> 3Q -> 1Q -> -1Q # approach: 96-phase cycling the first pulse, the third pulse and the receiver # spin-5/2 echo amplitude optimization # versus the second-pulse duration # in three-pulse split-t1 phase-modulated 3Q-MAS sequence, # the first-pulse duration p1 = 9 micro seconds # the third-pulse duration p3 = 2 micro seconds spinsys { channels 27Al nuclei 27Al quadrupole 1 1 3e6 1 0 0 0 } par { spin_rate 5000 variable tsw 0.25 sw 1.0e6/tsw np 34 crystal_file rep100 gamma_angles 10 proton_frequency 800e6 start_operator I1z detect_operator I1c verbose 1101 variable rf 90000 variable rf3 90000 variable p1 9 variable p3 2 variable N1 12 variable N3 8 variable deltap1 -3 variable deltap3 2 } proc pulseq {} { global par maxdt $par(tsw) pulse $par(p1) $par(rf) $par(ph1) store 10 reset acq $par(ph31) 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) $par(ph3) acq [expr $par(ph31) - 90] } } proc main {} { global par for {set j1 0} {$j1 < $par(N1)} {incr j1} { set par(ph1) [expr $j1*360./$par(N1)] for {set j3 0} {$j3 < $par(N3)} {incr j3} { set par(ph3) [expr $j3*360./$par(N3)] set par(ph31) [expr $par(deltap1)*$par(ph1) + $par(deltap3)*$par(ph3) ] 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 27Al: " puts [resfreq 27Al $par(proton_frequency)] } # SIMP # NP=17 # SW=4000000 # TYPE=FID # DATA # 0 0 # -3.07985342 -3.71794609e-15 # -9.85579627 -3.57353036e-15 # -15.9036315 5.85642645e-15 # -18.3668067 4.60742555e-15 # -17.1645919 6.9388939e-15 # -14.3785118 5.44009282e-15 # -12.2391223 3.88578059e-15 # -11.5118452 4.4408921e-16 # -11.3399754 5.55111512e-17 # -10.1810732 -4.38538095e-15 # -7.48105973 -2.33146835e-15 # -4.38104139 -3.94129174e-15 # -2.77190901 -2.77555756e-15 # -4.06327694 -2.24126273e-15 # -8.09039716 -1.47104551e-15 # -12.8560399 -1.19348975e-15 # -15.9486759 -3.77475828e-15 # -16.1749676 -3.71924713e-15 # -13.973069 -1.49880108e-15 # -10.8096564 -3.88578059e-16 # -8.08663483 2.49800181e-15 # -6.33751458 3.44169138e-15 # -5.30053911 3.88578059e-15 # -4.58284062 2.72004641e-15 # -4.2680421 1.69309011e-15 # -4.92113574 4.4408921e-16 # -6.91609465 1.11022302e-16 # -9.71827688 -1.33226763e-15 # -11.9888904 -3.5249581e-15 # -12.5786063 -3.05311332e-15 # -11.3827632 -4.26741975e-15 # -9.29662225 -2.61943245e-16 # -7.35497174 1.11022302e-16 # END