# na-shifted-echo-echo-phase-filtering-p2.in # approach: filtering the desired coherence (-3Q) after the first pulse # and (1Q) after the second pulse # spin-3/2 echo amplitude optimization # versus the second-pulse duration # in three-pulse shifted-echo phase-modulated -3Q-MAS sequence, # the first-pulse duration p1 = 4 micro seconds # the third-pulse duration p3 = 2 micro seconds spinsys { channels 23Na nuclei 23Na quadrupole 1 1 1e6 1 0 0 0 } par { spin_rate 5000 variable tsw 0.25 sw 1.0e6/tsw np 17 crystal_file rep100 gamma_angles 10 proton_frequency 800e6 start_operator 96*I1z detect_operator I1c verbose 1101 variable rf 90000 variable rf3 93000 variable p1 4 variable p3 2 } proc pulseq {} { global par maxdt $par(tsw) matrix set 1 coherence {-3} matrix set 2 coherence {1} pulse $par(p1) $par(rf) x store 10 reset acq -y for {set i 1} {$i < $par(np)} {incr i} { reset prop 10 filter 1 pulse [expr ($i)*$par(tsw)] $par(rf) x filter 2 pulse $par(p3) $par(rf3) x acq -y } } proc main {} { global par fsave [fsimpson] $par(name).fid puts "Larmor frequency (Hz) of 23Na: " puts [resfreq 23Na $par(proton_frequency)] } # SIMP # NP=17 # SW=4000000 # TYPE=FID # DATA # 0 0 # 0.00266047223 -7.30771317e-19 # 0.040473719 -6.86363939e-18 # 0.188538494 -3.66833489e-17 # 0.531256559 -8.72235634e-17 # 1.12216689 -1.53729775e-16 # 1.95733183 -3.48473229e-16 # 2.97155439 -6.04913785e-16 # 4.05516759 -8.66851094e-16 # 5.08119738 -8.96907503e-16 # 5.93169717 -1.16722628e-15 # 6.51653318 -1.20075084e-15 # 6.78354426 -1.27144605e-15 # 6.72208875 -1.20228414e-15 # 6.36176417 -1.18456718e-15 # 5.76649732 -1.08083753e-15 # 5.02372271 -1.01290893e-15 # END