# al-split-t1-5Q-antiecho-selective-80-phase-cycling-p3.in # oQ -> 5Q -> (5Q, 1Q, -3Q) -> -1Q # approach: 80-phase cycling the first pulse, the third pulse and the receiver # spin-5/2 antiecho amplitude optimization # versus the third-pulse duration # in three-pulse split-t1 phase-modulated 5Q-MAS sequence, # the first-pulse duration p1 = 15 micro seconds # the second-pulse duration p2 = 8.5 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 rep10 gamma_angles 10 proton_frequency 800e6 start_operator I1z detect_operator I1c verbose 1101 variable rf 90000 variable rf3 9000 variable p1 15 variable p2 8.5 variable N1 20 variable N3 4 variable deltap1 -5 variable deltap3 2 } proc pulseq {} { global par maxdt $par(tsw) acq $par(ph31) pulse $par(p1) $par(rf) $par(ph1) pulse $par(p2) $par(rf) x for {set i 1} {$i < $par(np)} {incr i} { pulse $par(tsw) $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 # -0.00614629303 -4.71844785e-16 # -0.0243737015 1.66533454e-16 # -0.0542301932 -1.94289029e-16 # -0.0952216659 2.35922393e-16 # -0.146949295 -2.08166817e-16 # -0.20911255 2.77555756e-16 # -0.281365653 8.32667268e-17 # -0.36316571 2.49800181e-16 # -0.453769068 -8.8817842e-16 # -0.55235374 -4.30211422e-16 # -0.658119522 1.38777878e-17 # -0.770320559 -9.71445147e-17 # -0.888300885 -4.16333634e-17 # -1.01151624 -7.21644966e-16 # -1.13944063 -5.689893e-16 # -1.2714173 -3.46944695e-16 # -1.40665567 -3.1918912e-16 # -1.54436361 7.63278329e-16 # -1.68377368 3.74700271e-16 # -1.82401679 -5.55111512e-17 # -1.96406026 4.4408921e-16 # -2.10278836 -6.9388939e-17 # -2.23910058 3.95516953e-16 # -2.37198657 -4.85722573e-17 # -2.50060316 -2.56739074e-16 # -2.62428227 -6.24500451e-17 # -2.74245049 -3.46944695e-16 # -2.85459994 -6.97358837e-16 # -2.96037071 -8.11850587e-16 # -3.05960406 -5.67254577e-16 # -3.15227936 -5.67254577e-16 # -3.23845233 -1.28196065e-15 # -3.31828167 -1.80064297e-15 # END