# al-split-t1-5Q-antiecho-160-phase-cycling-p2.in # oQ -> 5Q -> 1Q -> -1Q # approach: 160-phase cycling the first pulse, the third pulse and the receiver # spin-5/2 antiecho amplitude optimization # versus the second-pulse duration # in three-pulse split-t1 phase-modulated 5Q-MAS sequence, # the first-pulse duration p1 = 15 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 p1 15 variable p3 2 variable N1 20 variable N3 8 variable deltap1 -5 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 # 0.0510608567 -4.29517533e-15 # 0.568198454 -4.19109192e-15 # 1.48049325 -7.10542736e-15 # 1.38235231 -6.66133815e-15 # -0.491862054 -4.85722573e-15 # -2.34574461 3.88578059e-15 # -2.66048054 3.44516082e-15 # -1.85343817 -1.41553436e-15 # -1.23579646 -1.60982339e-15 # -1.80990131 4.5519144e-15 # -2.57058204 8.93729535e-15 # -2.65773322 1.66533454e-15 # -2.21188773 -3.85802501e-15 # -0.576324816 -4.74620343e-15 # 1.99675571 -3.07393e-15 # 4.49232299 -8.46545056e-15 # 5.36255371 -8.82627305e-15 # 2.86848469 -8.8817842e-16 # -0.922828507 1.63064007e-15 # -2.62436675 -1.66533454e-16 # -2.06102322 -4.88498131e-15 # -1.27580563 -7.88258347e-15 # -0.253584248 -5.13478149e-15 # 0.985279132 1.02435421e-15 # 1.14610862 2.58126853e-15 # 0.450953683 -1.11022302e-15 # -0.54588908 -7.63278329e-16 # -2.23443906 9.85322934e-16 # -3.32110451 -4.82947016e-15 # -3.09230707 -5.71764858e-15 # -2.9075078 -2.08166817e-15 # -2.98173125 4.82947016e-15 # -2.21555186 9.98853777e-15 # END