Setting up non-uniformly sampled spectra/NUS guide for Bruker according to Arrowsmith group in Toronto: Difference between revisions
Line 104: | Line 104: | ||
It can be called by typing '''sparse''' from the command prompt. | It can be called by typing '''sparse''' from the command prompt. | ||
[[Image:TopspinNUSSnap0.png|thumb|right| | [[Image:TopspinNUSSnap0.png|thumb|right|700px|Figure 1. Changing the pulse program and calling sparse.py script]] | ||
After the initial greeting, you will be guided through a few steps as shown below. | After the initial greeting, you will be guided through a few steps as shown below. |
Revision as of 20:52, 11 December 2009
__NORICHEDITOR__
Setting up acquisition of non-uniformly sampled 3D spectra on Bruker spectrometers with TopSpin
__NORICHEDITOR__
A standardized way to set non-uniform sampling (NUS) for TopSpin is in the making, however, while it is not yet included in the current distribution, this page will focus on the way Non-uniform sampling of 3D spectra is implemented by the Arrowsmith group in Toronto (University Health Network and University of Toronto). It is tailored for the subsequent processing with MDDGUI software, described elsewhere on this site.
To set up 3D NUS experiments according to this method one needs the following:
- Modified pulse program(s)
- A Jython script for generating the sampling schedule
__NORICHEDITOR__
Pulse program modifications
__NORICHEDITOR__
First of all, one has to set up a regular experiment with all pulses and delays calibrated as usual. Once this step is finished, one needs to replace the pulse program by a modified one.
This section briefly describes the essential pulse sequence modifications. The 15N-NOESY-HSQC pulse program is used as an example. In most cases one just need to copy these changes into relevant pulse program.
- Additional definitions should come with other definitions in the file:
;********Added for NUS******* define delay short "short=50u" "l2=td2/2" "l3=td1/2" ;^^^^^^^^^^^^^^^^^^^^^^^^^^^^
- Main if statement determines whether an Echo/AntiEcho or Real/Imaginary pair is recorded or skipped. The information is taken from the vdlist file which is generated by a separate script before running the experiments. If the current value of vd is 1u, the pair is skipped (goto 100 statement), if it is 3u, it is recorded. This addition replaces the text commented out and comes before the d1 delay:
;1 ze ; d11 pl16:f3 ;2 d11 do:f3 ;3 d12 ;********Added for NUS******* 1 ze d11 pl16:f3 2 d11 do:f3 short*2 3 short*5 4 1u if "vd < 2u" goto 100 short*6 5 short 6 d11 ;^^^^^^^^^^^^^^^^^^^^^^^^^^^^ d1
- The final modification comes at the end of the pulse program and uses explicit looping (old style Bruker pulse programming). Note that F1PH and F2EA statements are commented out and the phase and delay increments or decrements are stated explicitly. These increments or decrements need to be executed even if the particular fids are skipped
go=2 ph31 cpd3:f3 ;********Added for NUS******* d11 do:f3 wr #0 if #0 zd ;---------------------------- 100 short igrad EA ; even if the pair is skipped, phases have to be increment appropriately. short ip5*2 ; especially important if Complex or States-TPPI mode used for this dimension lo to 3 times 2 ; end of Echo-AntiEcho loop short id10 ; increment delay and phases for 15N short ip3*2 short ip6*2 short ip31*2 short ivd ; check the next value of vd lo to 4 times l2 ; go back to beginng with the new vd value short rd10 ; reset 15N delays and phases short rp3 short rp6 short rp31 short ip8 ; increment 1H-indirect phases short ip9 lo to 5 times 2 short id0 ; increment 1H-indirect delay. lo to 6 times l3 1m ;^^^^^^^^^^^^^^^^^^^^^^^^^^^^ ; d11 do:f3 mc #0 to 2 ; F1PH(rd10 & rp3 & rp6 & rp31 & ip8 & ip9, id0) ; F2EA(igrad EA & ip5*2, id10 & ip3*2 & ip6*2 & ip31*2) exit
A good number of pulse programs are already modified and are available upon request from the Arrowsmith group (agutmana@uhnres.utoronto.ca):
|
|
__NORICHEDITOR__
Generating the sampling schedule and running the experiment
__NORICHEDITOR__
Once the experiment is otherwise set up and the pulse program modifications are in place, the last step before running the experiment is to generate the sampling schedule. Depending on whether your computer runs Linux or Windows, the generator differs slightly (mainly in the usage of forward and backslashes). Both scripts are available upon request from the Arrowsmith group (together with the above pulse programs). The corresponding script should be renamed sparse.py and put in /opt/topspin/stan/nmr/..... on your computer.
It can be called by typing sparse from the command prompt.
After the initial greeting, you will be guided through a few steps as shown below.
- All vdlist files should be copied to a standard directory. If the suggested directory is wrong, please modify.
- Main relevant parameters are presented in this window. The level of "sparsing" allows to fine tune the acquisition time. It should generally be 30% or more. T2sp parameters dictate the exponential biasing of the sampling and only need to be approximate. In the Arrowsmith group we use 20ms for 13C, 35 ms for 1H and 50ms for 15N. Usually the values are guessed correctly, but you can modify them if desired.
- Concluding window shows the final information on the number of FIDs to be recorded. Note that the experimental time shown by expt command will be wrong and should be multiplied by the "level of sparsing" parameter to give a correct estimate.
- To run the experiment type zg or multizg as the case may be. To process the spectrum you will need the vdlist file in addition to the usual parameter files.