# al-spam-echo-3Q-selective-12-cycling-p3.in # approach: 12-phase cycling the first pulse and the receiver # spin-5/2 echo amplitude optimization # versus the third-pulse duration # in three-pulse spam-echo 3Q-MAS sequence, # the first-pulse duration p1 = 4 micro seconds # the second-pulse duration p2 = 4 micro seconds spinsys { channels 27Al nuclei 27Al quadrupole 1 1 1e6 1 0 0 0 } par { spin_rate 5000 variable tsw 0.25 sw 1.0e6/tsw np 17 crystal_file rep10 gamma_angles 10 proton_frequency 800e6 start_operator I1z detect_operator I1c verbose 1101 variable rf 90000 variable rf3 9300 variable p1 4 variable p2 4 variable NA 12 variable deltapA -3 } proc pulseq {} { global par maxdt $par(tsw) acq $par(phREC) pulse $par(p1) $par(rf) $par(phA) pulse $par(p2) $par(rf) x for {set i 1} {$i < $par(np)} {incr i} { pulse $par(tsw) $par(rf3) x acq [expr $par(phREC) - 90] } } proc main {} { global par for {set jA 0} {$jA < $par(NA)} {incr jA} { set par(phA) [expr $jA*360./$par(NA)] set par(phREC) [expr $par(deltapA)*$par(phA)] 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.01960994 -4.21884749e-15 # -1.16326606 -4.4408921e-15 # -1.29689289 -3.77475828e-15 # -1.41875902 -3.77475828e-15 # -1.5275976 -3.33066907e-15 # -1.62269124 -3.55271368e-15 # -1.70392475 -3.77475828e-15 # -1.7717808 -2.66453526e-15 # -1.82727102 -2.66453526e-15 # -1.87181497 -2.88657986e-15 # -1.90709388 -9.99200722e-16 # -1.93490994 -2.10942375e-15 # -1.95707444 -1.33226763e-15 # -1.97533371 -9.99200722e-16 # -1.99132573 -8.8817842e-16 # -2.00655016 -9.99200722e-16 # END