# al-split-t1-echo-selective-48-phase-cycling-p1.in # oQ -> 3Q -> (5Q, 1Q, -3Q) -> -1Q # approach: 48-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 3Q-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 9000 variable p2 4 variable p3 2 variable N1 12 variable N3 4 variable deltap1 -3 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 # -0.000685209338 1.22124533e-15 # -0.0133621968 -4.4408921e-16 # -0.073438793 1.22124533e-15 # -0.177974954 5.55111512e-17 # -0.276409656 -2.49800181e-16 # -0.350511249 -1.11022302e-15 # -0.418988113 -1.11022302e-15 # -0.489066376 -2.22044605e-16 # -0.540587541 -7.77156117e-16 # -0.541199235 -6.66133815e-16 # -0.487321649 -6.9388939e-17 # -0.401648033 -3.05311332e-16 # -0.314194558 7.49400542e-16 # -0.238447304 8.32667268e-16 # -0.165826822 3.05311332e-16 # -0.0959536221 1.94289029e-16 # -0.0478289791 8.32667268e-16 # -0.0370582527 9.99200722e-16 # -0.057632351 2.22044605e-16 # -0.0914969318 -3.33066907e-16 # -0.10386823 -2.22044605e-16 # -0.0909332183 -7.35522754e-16 # -0.0881315657 -5.55111512e-17 # -0.0976315833 -3.60822483e-16 # -0.116685574 3.60822483e-16 # -0.141716467 9.71445147e-17 # -0.173295533 -4.51028104e-16 # -0.218664133 -4.64905892e-16 # -0.249605394 -1.94289029e-16 # -0.253860842 -1.11022302e-16 # -0.218680498 0 # -0.131366224 -3.88578059e-16 # -0.0476156066 1.66533454e-16 # END