# al-split-t1-5Q-echo-160-phase-cycling-p1.in # oQ -> -5Q -> 1Q -> -1Q # approach: 160-phase cycling the first pulse, the third pulse and the receiver # spin-5/2 echo amplitude optimization # versus the first-pulse duration # in three-pulse split-t1 phase-modulated 5Q-MAS sequence, # the second-pulse duration p2 = 4 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.5 sw 1.0e6/tsw np 34 crystal_file rep10 gamma_angles 10 proton_frequency 800e6 start_operator I1z detect_operator I1c verbose 1101 variable rf 90000 variable rf3 90000 variable p2 4 variable p3 2 variable N1 20 variable N3 8 variable deltap1 5 variable deltap3 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 [expr $par(ph31) - 90] reset prop 2 } } 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 # 1.34229531e-06 -1.11022302e-15 # 0.000382444094 0 # 0.0078618998 7.21644966e-15 # 0.064464131 7.21644966e-15 # 0.288092813 6.88338275e-15 # 0.773671356 1.08246745e-14 # 1.37628568 8.21565038e-15 # 1.79294308 1.18238752e-14 # 1.8958684 1.63757896e-14 # 1.79814074 1.4974133e-14 # 1.65550539 8.9928065e-15 # 1.54217929 1.94289029e-16 # 1.47824636 -5.48172618e-15 # 1.4771626 -5.66213743e-15 # 1.53269598 -8.10462808e-15 # 1.59733862 -6.605827e-15 # 1.64046333 -1.80411242e-15 # 1.67429122 -1.37390099e-15 # 1.68721399 2.77555756e-17 # 1.64358884 0 # 1.58359163 -4.23272528e-15 # 1.62691618 -5.74540415e-15 # 1.83649461 -4.71844785e-15 # 2.13566974 2.77555756e-15 # 2.37606122 6.80011603e-15 # 2.48355882 1.15272375e-14 # 2.49982262 1.25732758e-14 # 2.46411229 1.24622535e-14 # 2.37173255 8.37524494e-15 # 2.23997292 9.71445147e-15 # 2.10332842 8.43769499e-15 # 1.93797142 1.03805853e-14 # 1.68091988 1.07691633e-14 # END